|
%I #41 Nov 03 2023 10:58:29
%S 31,191,257,367,3061,110503,1145393
%N Numbers k such that (36^k + 1)/37 is prime.
%C All such numbers k are prime.
%C Note that a(6) = 110503 corresponds to (36^110503 + 1)/37, which is only a probable prime with 171975 digits.
%C The primes corresponding to the terms of this sequence have 1 as their last digit and an even number as their next-to-last digit. - _Iain Fox_, Dec 08 2017
%H J. Brillhart et al., <a href="http://dx.doi.org/10.1090/conm/022">Factorizations of b^n +- 1</a>, Contemporary Mathematics, Vol. 22, Amer. Math. Soc., Providence, RI, 3rd edition, 2002.
%H H. Dubner and T. Granlund, <a href="http://www.cs.uwaterloo.ca/journals/JIS/VOL3/DUBNER/dubner.html">Primes of the Form (b^n+1)/(b+1)</a>, J. Integer Sequences, 3 (2000), #P00.2.7.
%H H. Lifchitz, <a href="http://www.primenumbers.net/Henri/us/MersFermus.htm">Mersenne and Fermat primes field</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/Repunit.html">Repunit</a>.
%t Do[ p=Prime[n]; If[ PrimeQ[ (36^p + 1)/37 ], Print[p] ], {n, 1, 9592} ]
%o (PARI) is(n)=isprime((36^n+1)/37) \\ _Charles R Greathouse IV_, Feb 17 2017
%Y Cf. A000978 = numbers n such that (2^n + 1)/3 is prime. Cf. A007658, A057171, A057172, A057173, A057175, A001562, A057177, A057178, A057179, A057180, A057181, A057182, A057183, A057184, A057185, A057186, A057187, A057188, A057189, A057190, A057191, A071380, A071381, A071382, A084741, A084742, A065507, A126659, A126856, A185240.
%K hard,more,nonn
%O 1,1
%A _Robert Price_, Sep 15 2013
%E a(6) = 110503 (posted by Lelio R. Paula on primenumbers.net) from _Paul Bourdelais_, Dec 08 2017
%E a(7) from _Paul Bourdelais_, Nov 03 2023
|
