%I #25 Jul 05 2020 08:23:09
%S 1,4,5,11,26,29,38,65,166,490,641,2300,9440,44741,65296,161930,586240
%N Numbers n such that 8*9^n +1 is prime.
%C The even members of A005538, halved. - _Jeppe Stig Nielsen_, Oct 13 2019
%t Do[If[PrimeQ[8*9^n + 1], Print[n]], {n, 0, 1000}]
%t Select[Range[3000], PrimeQ[8 9^# + 1] &] (* _Vincenzo Librandi_, Sep 28 2012 *)
%o (PARI) is(n)=ispseudoprime(8*9^n+1) \\ _Charles R Greathouse IV_, Mar 01 2017
%Y Cf. A005538.
%K nonn,more
%O 1,2
%A _Robert G. Wilson v_, Aug 22 2000
%E a(12) from _Vincenzo Librandi_, Sep 28 2012
%E a(13)-a(16) from _Robert Price_, Feb 22 2016
%E a(17) from _Paul S. Vanderveen_, Jul 04 2020
|