OFFSET
1,4
PROG
(PARI) a(n)=sum(k=1, n, if(numerator(bernfrac(2*n))%k, 0, 1))
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Jun 21 2003
STATUS
approved