A101738 Indices of primes in sequence defined by A(0) = 41, A(n) = 10*A(n-1) + 81 for n > 0.
0, 1, 3, 6, 23, 171, 172, 217, 438, 2562, 2781, 2879, 3714, 4391, 10819, 14591, 25053, 27768, 40167, 180991, 193906
Offset: 1
Examples
49991 is prime, hence 3 is a term.
References
- Klaus Brockhaus and Walter Oberschelp, Zahlenfolgen mit homogenem Ziffernkern, MNU 59/8 (2006), pp. 462-467.
Links
Programs
-
Mathematica
Select[Range[0, 1000], PrimeQ[(450*10^# - 81)/9] &] (* Robert Price, May 30 2015 *)
-
PARI
a=41;for(n=0,1500,if(isprime(a),print1(n,","));a=10*a+81)
-
PARI
for(n=0,1500,if(isprime((450*10^n-81)/9),print1(n,",")))
Formula
a(n) = A103001(n) - 1.
Extensions
More terms from Herman Jamke (hermanjamke(AT)fastmail.fm), Jan 01 2008
a(17)-a(19) from Kamada data by Ray Chandler, May 01 2015
a(20)-a(21) from Robert Price, May 30 2015
Comments