OFFSET
0,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..75
FORMULA
a(n) = Sum_{k=0..n} abs( A118686(n,k) ).
MATHEMATICA
g[n_]:= If[PrimeQ[n], n, 1]; p[n_]:= p[n]= If[n==0, 1, g[n]*p[n-1]];
A119489= Flatten[Join[{{1}}, Table[Apply[Plus, Abs[Reverse[ CoefficientList[Product[x-p[n], {n, 0, m}], x]]]], {m, 0, 30}]]]
PROG
(Magma)
R<x>:=PowerSeriesRing(Rationals(), 50);
g:= func< n | IsPrime(n) select n else 1 >;
p:=[1] cat [n le 1 select 1 else g(n)*Self(n-1): n in [1..50]];
T:= func< n, k | k eq 0 select 1 else Coefficient(R!( (&*[x-p[j+1]: j in [0..n-1]]) ), n-k) >;
[(&+[Abs(T(n, k)): k in [0..n]]): n in [0..30]]; // G. C. Greubel, Mar 31 2024
(SageMath)
def g(n): return n if is_prime(n) else 1
def p(n): return 1 if n==0 else g(n)*p(n-1)
def T(n, k): return 1 if k==0 else ( product(x-p(j) for j in range(n)) ).series(x, n+2).list()[n-k] # T = A118686
[sum(abs(T(n, k)) for k in range(n+1)) for n in range(31)] # G. C. Greubel, Mar 31 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Roger L. Bagula, May 25 2006
EXTENSIONS
Terms a(12) onward added by G. C. Greubel, Mar 31 2024
STATUS
approved