(PARI) sopf(n, s, fac, i)=fac=factor(n); for(i=1, matsize(fac)[1], s=s+fac[i, 1]); return(s);
for(n=1, 23, print(sopf(n!+1)))
(PARI) sopf(n)= { local(f, s=0); f=factor(n); for(i=1, matsize(f)[1], s+=f[i, 1]); return(s) }
{ for (n=1, 60, write("b064011.txt", n, " ", sopf(n! + 1)) ) } \\
Harry J. Smith, Sep 06 2009