OFFSET
1,1
COMMENTS
A prime p is in the sequence if p is of the form 660*n + 551.
PROG
(PARI) forstep(p=11, 49391, 12, if(isprime(p)&&!Mod(numerator(bernfrac(p-1)), (p-1)/2)==0, print1(p, ", ")));
CROSSREFS
KEYWORD
nonn
AUTHOR
Arkadiusz Wesolowski, Nov 17 2013
STATUS
approved