OFFSET
1,5
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..16384
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..65537
MATHEMATICA
Array[Block[{w = {{#}}}, While[w[[-1]] != {1}, w = Join[w, {Union@ Flatten[# - #/FactorInteger[#][[All, 1]] & /@ w[[-1]] ]}]]; Length@ Union@ Flatten@ w - DivisorSigma[0, #]] &, 93] (* Michael De Vlieger, May 09 2020, after Robert G. Wilson v at A332809 *)
PROG
(PARI)
up_to = 20000;
A332809list(up_to) = { my(v=vector(up_to)); v[1] = Set([1]); for(n=2, up_to, my(f=factor(n)[, 1]~, s=Set([n])); for(i=1, #f, s = setunion(s, v[n-(n/f[i])])); v[n] = s); apply(length, v); }
v332809 = A332809list(up_to);
A332809(n) = v332809[n];
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, May 09 2020
STATUS
approved