|
|
PROG
|
(Perl) use bigint; use ntheory ":all"; say "$_ ", ramanujan_tau(10 ** $_) for 0..19; # Dana Jacobsen, Sep 05 2015
(PARI) taup(p, e)={
if(e==1,
(65*sigma(p, 11)+691*sigma(p, 5)-691*252*sum(k=1, p-1, sigma(k, 5)*sigma(p-k, 5)))/756
,
my(t=taup(p, 1));
sum(j=0, e\2,
(-1)^j*binomial(e-j, e-2*j)*p^(11*j)*t^(e-2*j)
)
)
};
a(n)=taup(5, n)*taup(2, n) \\ Charles R Greathouse IV, Sep 06 2015
|
