editing
approved
editing
approved
Alois P. Heinz, <a href="/A284943/b284943.txt">Table of n, a(n) for n = 1..3313</a>
b:= proc(n) option remember; nops(ifactors(n)[2])=1 end:
a:= proc(n) option remember; `if`(n=0, 0, add(a(n-j)+
`if`(b(j), ceil(2^(n-j-1)), 0), j=1..n))
end:
seq(a(n), n=1..33); # Alois P. Heinz, Aug 07 2019
approved
editing
proposed
approved
editing
proposed
(PARI) x='x+O('x^34); concat([0], Vec(sum(k=2, 34, (1\omega(k))*x^k*(1 - x)^2/(1 - 2*x)^2))) \\ Indranil Ghosh, Apr 06 2017
proposed
editing
editing
proposed
nmax = 33; Rest[CoefficientList[Series[Sum[Floor[1/PrimeNu[k]] x^k (1 - x)^2/(1 - 2 x)^2, {k, 2, 2 nmax}], {x, 0, nmax}], x]]
nmax = 33; Rest[CoefficientList[Series[Sum[Floor[1/PrimeNu[k]] x^k (1 - x)^2/(1 - 2 x)^2, {k, 2, 2 nmax}], {x, 0, nmax}], x]]
proposed
editing
editing
proposed