|
%I #23 Feb 01 2021 19:00:05
%S 16,96,376,1160,3121,7532,16754,34796,68339,127952,229956,398688,
%T 669781,1094076,1742710,2713604,4139111,6195712,9115304,13199072,
%U 18833449,26509260,36843322,50603884,68740107,92414192,123039628,162323200,212312453,275448380
%N Four-fold convolution of primes with themselves.
%H Alois P. Heinz, <a href="/A014344/b014344.txt">Table of n, a(n) for n = 0..1000</a>
%F G.f.: ((1/x)*Sum_{k>=1} prime(k)*x^k)^4. - _Ilya Gutkovskiy_, Mar 10 2018
%p b:= proc(n, k) option remember; `if`(k=1, ithprime(n+1),
%p add(b(j, floor(k/2))*b(n-j, ceil(k/2)), j=0..n))
%p end:
%p a:= n-> b(n, 4):
%p seq(a(n), n=0..35); # _Alois P. Heinz_, Mar 10 2018
%t b[n_, k_] := b[n, k] = If[k==1, Prime[n+1], Sum[b[j, Floor[k/2]] b[n-j, Ceiling[k/2]], {j, 0, n}]];
%t a[n_] := b[n, 4];
%t a /@ Range[0, 35] (* _Jean-François Alcover_, Nov 16 2020, after _Alois P. Heinz_ *)
%o (PARI) my(N = 50, x = 'x + O('x^N)); Vec(((1/x)*sum(k=1, N, prime(k)*x^k))^4) \\ _Michel Marcus_, Mar 10 2018
%Y Cf. A000040, A014342, A014343.
%Y Column k=4 of A340991.
%K nonn
%O 0,1
%A _N. J. A. Sloane_
|
