OFFSET
1,1
COMMENTS
A prime p is in the sequence if p is of the form 660*n + 109.
EXAMPLE
109 is in the sequence because B(110) = (5 * 157 * 76493 * C)/1518 (where C is some large, unfactored composite number), the numerator of which is not divisible by 110/2 = 5 * 11.
97 is not in the sequence because B(98) = (7^2 * 2857 * 3221 * C)/6, the numerator of which is divisible by 98/2 = 49 = 7^2.
MATHEMATICA
Select[12Range[864] + 1, PrimeQ[#] && Not[Divisible[Numerator[Bernoulli[# + 1]], (# + 1)/2]] &] (* Alonso del Arte, Nov 17 2013 *)
PROG
(PARI) forstep(p=1, 46309, 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