OFFSET
1,1
PROG
(PARI) {s=0; h=1; forprime(p=3, 2400000, q=(p+1)/2; s=s+sum(k=h, q, 1/k^3); if(numerator(s)%p==0, print1(p, ", ")); h=q+1)} - Klaus Brockhaus, Nov 10 2006
CROSSREFS
KEYWORD
more,nonn
AUTHOR
Alexander Adamchuk, Nov 07 2006
EXTENSIONS
a(5) from Klaus Brockhaus, Nov 10 2006
STATUS
approved