OFFSET
1,1
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..10000 (first 800 terms from Harry J. Smith)
PROG
(PARI) usigma(n, s=1, fac, i)=fac=factor(n); for(i=1, matsize(fac)[1], s=s*(1+fac[i, 1]^fac[i, 2])); return(s); j=[]; for(n=1, 50000, if(usigma(n)==usigma(n+1), j=concat(j, n))); j
(PARI) usigma(n)= { local(f, s=1); f=factor(n); for(i=1, matsize(f)[1], s*=1 + f[i, 1]^f[i, 2]); return(s) }
{ n=0; s=0; for (m=1, 10^9, us=usigma(m+1); if(s==us, write("b064125.txt", n++, " ", m); if (n==800, break)); s=us ) } \\ Harry J. Smith, Sep 08 2009
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Jason Earls, Sep 10 2001
STATUS
approved