OFFSET
1,3
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
PROG
(PARI) a(n) = my(v = binomial(4*n, 2*n)/binomial(2*n, n)); bigomega(numerator(v)) - bigomega(denominator(v)); \\ Michel Marcus, Sep 30 2013
(PARI) vp(n, p)=my(s); while(n\=p, s+=n); s
a(n)=my(s); forprime(p=2, 4*n, s+=vp(4*n, p)-3*vp(2*n, p)+2*vp(n, p)); s \\ Charles R Greathouse IV, Sep 30 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved