A233906 - OEIS

%I #9 Dec 18 2013 13:39:20

%S 2,1499,1783,9719,9311,67883,134947,203317,189433,560171,438533,

%T 943849,640973,578827,2172383,28687,1505657,7595033,2822971,1242379,

%U 22899523,9232219,5730031,12336083,3487607,35451433,12174803,10234079,84459019,68736683,44671169,85507057

%N Primes of the form (3^k mod k^3) + 1, in order of increasing k.

%H K. D. Bajpai, <a href="/A233906/b233906.txt">Table of n, a(n) for n = 1..1000</a>

%e 1499 is in the sequence because (3^13 mod 13^3) + 1 = 1499 which is prime.

%e 9719 is in the sequence because (3^29 mod 29^3) + 1 = 9719 which is prime.

%p KD := proc() local a; a:=3^n mod n^3 + 1; if isprime(a) then RETURN (a); fi; end: seq(KD(), n=1..1000);

%Y Cf. A000040 (prime numbers).

%Y Cf. A007519 (primes congruent to 1 mod 8).

%K nonn

%O 1,1

%A _K. D. Bajpai_, Dec 17 2013