|
%I M2337 #33 Sep 08 2022 08:44:33
%S 0,1,3,4,7,9,12,18,22,102,112,157,162,289,619,763,1389,1783,1882,3294,
%T 3567,13297,14932,18954,19612,23598,33882,66874,70546,86568,187626,
%U 190738
%N Numbers k such that 10*3^k + 1 is prime.
%C a(33) > 2*10^5. - _Robert Price_, Mar 16 2014
%C All terms are verified primes (i.e., not merely probable primes). - _Robert Price_, Mar 16 2014
%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
%H H. C. Williams and C. R. Zarnke, <a href="https://doi.org/10.1090/S0025-5718-1972-0314747-X">Some prime numbers of the forms 2*3^n+1 and 2*3^n-1</a>, Math. Comp., 26 (1972), 995-998.
%t Do[ If[ PrimeQ[ 10*3^n + 1], Print[n]], {n, 0, 6810}]
%o (Magma) [n: n in [0..3567] | IsPrime(10*3^n + 1) ]; // _Vincenzo Librandi_, Sep 26 2012
%o (PARI) is(n)=ispseudoprime(10*3^n+1) \\ _Charles R Greathouse IV_, Feb 20 2017
%K hard,nonn
%O 1,3
%A _N. J. A. Sloane_
%E More terms from _Robert G. Wilson v_, Sep 07 2000
%E a(1)=0 added and typo in Mathematica program fixed by _Vincenzo Librandi_, Sep 26 2012
%E a(22)-a(32) from _Robert Price_, Mar 16 2014
|
