OFFSET
1,2
COMMENTS
FORMULA
a(n) <= 2^(n-1).
a(p) = 2^(p-1) for a prime p.
MATHEMATICA
q[n_] := UnsameQ @@ Differences[Divisors[n]]; seq[len_] := Module[{s = Table[0, {len}], c = 0, n = 1, i}, While[c < len, i = DivisorSigma[0, n]; If[i <= len && s[[i]] == 0 && q[n], c++; s[[i]] = n]; n++]; s]; seq[20]
PROG
(PARI) isA060683(k) = my(d=divisors(k)); #Set(vector(#d-1, k, d[k+1]-d[k])) == #d-1; \\ Michel Marcus at A060683
list(len) = {my(s = vector(len), c = 0, n = 1, i); while(c < len, i = numdiv(n); if(i <= len && s[i] == 0 && isA060683(n), c++; s[i] = n); n++); s; }
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Sep 16 2025
STATUS
approved
