%I #15 Jan 12 2021 10:10:45
%S 1,0,2,2,1,6,1,8,6,6,14,6,18,14,18,24,23,30,35,38,46,54,55,74,72,90,
%T 100,106,128,136,152,178,185,216,238,252,302,308,359,390,420,478,512,
%U 564,628,668,745,810,871,974,1035,1140,1238,1336,1459,1586,1700,1868,1993,2168,2354,2512,2751,2930,3177,3418,3677,3960
%N Expansion of Product_{k>=1} (1 + x^prime(k))^2.
%C Number of partitions of n into distinct prime parts, with 2 types of each part.
%C Self-convolution of A000586. - _Ilya Gutkovskiy_, Jan 19 2018
%H Vaclav Kotesovec, <a href="/A280245/b280245.txt">Table of n, a(n) for n = 0..10000</a> (terms 0..1000 from Seiichi Manyama)
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/PrimePartition.html">Prime Partition</a>
%H <a href="/index/Par#partN">Index entries for related partition-counting sequences</a>
%F G.f.: Product_{k>=1} (1 + x^prime(k))^2.
%F log(a(n)) ~ 2*Pi*sqrt(n/(3*log(n/2))). - _Vaclav Kotesovec_, Jan 12 2021
%e a(5) = 6 because we have [5], [5'], [3, 2], [3', 2], [3, 2'], [3', 2'].
%t nmax = 67; CoefficientList[Series[Product[(1 + x^Prime[k])^2, {k, 1, nmax}], {x, 0, nmax}], x]
%Y Cf. A000009, A000586, A000607, A022567, A070215.
%K nonn
%O 0,3
%A _Ilya Gutkovskiy_, Dec 29 2016