login
A183174
Numbers k such that (10^(2k+1) - 6*10^k - 1)/3 is prime.
45
1, 3, 7, 61, 90, 92, 269, 298, 321, 371, 776, 1567, 2384, 2566, 3088, 5866, 8051, 9498, 12635, 24512, 32521, 43982
OFFSET
1,2
COMMENTS
a(23) > 10^5. - Robert Price, Jan 29 2016
REFERENCES
C. Caldwell and H. Dubner, "Journal of Recreational Mathematics", Volume 28, No. 1, 1996-97, pp. 1-9.
FORMULA
a(n) = (A077775(n) - 1)/2.
MATHEMATICA
Do[If[PrimeQ[(10^(2n + 1) - 6*10^n - 1)/3], Print[n]], {n, 3000}]
PROG
(PARI) for(n=1, 1e3, if(ispseudoprime((10^(2*n+1)-6*10^n-1)/3), print1(n", "))) \\ Charles R Greathouse IV, Jul 15 2011
KEYWORD
nonn,base
AUTHOR
Ray Chandler, Dec 28 2010
EXTENSIONS
a(21)-a(22) from Robert Price, Jan 29 2016
STATUS
approved