OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..47 (terms below 10^11)
EXAMPLE
11 is in the sequence since uphi(11) * usigma(11) = 10 * 12 = uphi(12) * usigma(12) = 6 * 20 = 120.
MATHEMATICA
usigma[n_] := If[n == 1, 1, Times @@ (1 + Power @@@ FactorInteger[n])];
uphi[n_] := (Times @@ (Table[#[[1]]^#[[2]] - 1, {1}] & /@ FactorInteger[n]))[[1]]; u[n_] := uphi[n]*usigma[n]; aQ[n_] := u[n] == u[n + 1]; Select[Range[10^6], aQ]
PROG
(PARI) A191414(n) = {my(f = factor(n)); prod(i = 1, #f~, f[i, 1]^(2*f[i, 2])-1); }
lista(kmax) = {my(a1 = 1, a2); for(k = 2, kmax, a2 = A191414(k); if(a1 == a2, print1(k-1, ", ")); a1 = a2); } \\ Amiram Eldar, Nov 09 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Dec 29 2017
STATUS
approved