OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..1204 (terms 1..75 from Harry J. Smith)
MATHEMATICA
g[1]={1, 1}; g[n_] := {Times @@ ((f = FactorInteger[n])[[;; , 2]] + 1), Times @@ (1 + Power @@@ f)}; s={}; g1={0, 0}; Do[g2=g[n]; If[g1==g2, AppendTo[s, n-1]]; g1=g2, {n, 1, 50000}]; s (* Amiram Eldar, Jun 19 2019 *)
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); } for(n=1, 10^6, if(usigma(n)==usigma(n+1) && numdiv(n)==numdiv(n+1), print1(n, ", ")))
(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; d=0; for (m=1, 10^9, us=usigma(m + 1); ud=numdiv(m + 1); if(s==us && d==ud, write("b064348.txt", n++, " ", m); if (n==100, break)); s=us; d=ud ) } \\ Harry J. Smith, Sep 12 2009
CROSSREFS
KEYWORD
nonn
AUTHOR
Jason Earls, Oct 15 2001
STATUS
approved