%I #11 Mar 01 2022 07:25:03
%S 1,2,4,10,24,50,140,368,1152,3682,11784,39902,134612,463066,1635092,
%T 5818384,20684072,73693068,266943648,967762792,3533666568,13036452946,
%U 48102671884,178315730764,661567489568,2450447537226,9123572154720,34201574126260
%N a(n) is the constant term in expansion of Product_{k=1..n} (x^prime(k) + 1/x^prime(k))^2.
%p p:= proc(n) option remember; `if`(n=0, 1,
%p p(n-1)*(x^ithprime(n)+1/x^ithprime(n))^2)
%p end:
%p a:= n-> coeff(p(n), x, 0):
%p seq(a(n), n=0..30); # _Alois P. Heinz_, Jan 26 2022
%t p[n_] := p[n] = If[n == 0, 1, p[n - 1]*(x^Prime[n] + 1/x^Prime[n])^2];
%t a[n_] := Coefficient[p[n], x, 0];
%t Table[a[n], {n, 0, 30}] (* _Jean-François Alcover_, Mar 01 2022, after _Alois P. Heinz_ *)
%o (PARI) a(n) = polcoef (prod(k=1, n, (x^prime(k) + 1/x^prime(k))^2), 0); \\ _Michel Marcus_, Jan 21 2022
%Y Cf. A000980, A022894, A022896, A047653.
%K nonn
%O 0,2
%A _Ilya Gutkovskiy_, Jan 20 2022