OFFSET
1,2
COMMENTS
Numbers k such that A331973(k) = 1.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
f[p_, e_] := Module[{b = IntegerDigits[e, 2]}, m = Length[b]; Product[If[b[[j]] > 0, 1 + p^(2^(m - j)), 1], {j, 1, m}]]; is[1] = 0; is[n_] := Times @@ f @@@ FactorInteger[n] - n;
seq[max_] := Module[{v = Table[0, {max}], i}, Do[i = is[k] + 1; If[i <= max, v[[i]]++], {k, 1, max^2}]; -1 + Position[v, 1] // Flatten];
seq[500]
PROG
(PARI) s(n) = {my(f = factor(n), b); prod(i=1, #f~, b = binary(f[i, 2]); prod(k=1, #b, if(b[k], f[i, 1]^(2^(#b-k)) + 1, 1))) - n; }
lista(nmax) = {my(v = vector(nmax+1)); for(k=1, nmax^2, i = s(k) + 1; if(i <= nmax+1, v[i] += 1)); for(i = 1, nmax+1, if(v[i] == 1, print1(i-1, ", "))); }
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Mar 11 2023
STATUS
approved