OFFSET
1,2
LINKS
Harry J. Smith, Table of n, a(n) for n = 1..600
PROG
(PARI) sopfr(n, s, fac, i) = fac=factor(n); for(i=1, matsize(fac)[1], s=s+fac[i, 1]*fac[i, 2]); return(s); aa=-1; j=[]; for(n=1, 150000, x=sopfr(n); if(issquare(x), b=x; if(b>aa, aa=b; j=concat(j, n)))); j
(PARI) sopfr(n)= { local(f, s=0); f=factor(n); for(i=1, matsize(f)[1], s+=f[i, 1]*f[i, 2]); return(s) } { n=0; a=-1; for (m=1, 10^9, if (issquare(b=sopfr(m)), if(b>a, a=b; write("b064463.txt", n++, " ", m); if (n==600, break))) ) } \\ Harry J. Smith, Sep 15 2009
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Jason Earls, Oct 02 2001
STATUS
approved