OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..138 (terms below 10^10)
MATHEMATICA
usigma[1] = 1; usigma[n_] := Times @@ (1 + Power @@@ FactorInteger[n]); seq = {}; u1 = usigma[1]; u2 = usigma[2]; Do[u3 = usigma[n]; If[u3 == 2*u2 - u1, AppendTo[seq, n - 2]]; u1 = u2; u2 = u3, {n, 3, 10^4}]; seq (* Amiram Eldar, Jan 04 2020 *)
PROG
(PARI) usigma(n) = sumdiv(n, d, if(gcd(d, n/d)==1, d));
lista(nn) = {u = vector(nn, n, usigma(n)); v = vector(#u-1, n, u[n+1] - u[n]); w = vector(#v-1, n, v[n+1] - v[n]); x = select(y->y==0, w, 1); } \\ Michel Marcus, Aug 12 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Jason Earls, Nov 01 2001
EXTENSIONS
More terms from Michel Marcus, Aug 12 2018
STATUS
approved