OFFSET
1,3
COMMENTS
a(33) > 2*10^5. - Robert Price, Mar 16 2014
All terms are verified primes (i.e., not merely probable primes). - Robert Price, Mar 16 2014
REFERENCES
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
LINKS
H. C. Williams and C. R. Zarnke, Some prime numbers of the forms 2*3^n+1 and 2*3^n-1, Math. Comp., 26 (1972), 995-998.
MATHEMATICA
Do[ If[ PrimeQ[ 10*3^n + 1], Print[n]], {n, 0, 6810}]
PROG
(Magma) [n: n in [0..3567] | IsPrime(10*3^n + 1) ]; // Vincenzo Librandi, Sep 26 2012
(PARI) is(n)=ispseudoprime(10*3^n+1) \\ Charles R Greathouse IV, Feb 20 2017
CROSSREFS
KEYWORD
hard,nonn
AUTHOR
EXTENSIONS
More terms from Robert G. Wilson v, Sep 07 2000
a(1)=0 added and typo in Mathematica program fixed by Vincenzo Librandi, Sep 26 2012
a(22)-a(32) from Robert Price, Mar 16 2014
STATUS
approved