OFFSET
1,3
MATHEMATICA
mx = 120; s = {0, 1} ~Join~ Select[Range[2, mx], GCD @@ (Last /@ FactorInteger[#]) > 1 &]; Union[ Reap[Do[v = s[[i]] + s[[j]]; If[v <= mx, Sow@v], {i, Length@s}, {j, i}]][[2, 1]]] (* Giovanni Resta, Sep 19 2019 *)
PROG
(PARI) upto(n)={my(p=(1 + x + sum(k=2, sqrtint(n), sum(e=2, logint(n, k), x^(k^e))) + O(x*x^n))^2); select(i->polcoef(p, i), [0..n])} \\ Andrew Howroyd, Sep 22 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Elan Roth, Sep 18 2019
STATUS
approved