OFFSET
0,1
PROG
(PARI) bouay2(n) = bernfrac(n-1)+1/n for (i=1, 500, if(denominator(bouay2(prime(i)))==30, print(i)))
CROSSREFS
KEYWORD
nonn
AUTHOR
mohammed bouayoun (bouyao(AT)wanadoo.fr), Feb 10 2004
STATUS
approved