%I #23 Mar 01 2021 10:15:59
%S 1,0,1,3,10,40,153,921,5144,30717,230748,1766056,14052445,116580521,
%T 897876519,7657321097,75743979608,788733735080,7569825650083,
%U 75242386295617,831978453306391,9444103049405370,120064355466770831,1579842230380587833
%N Number of partitions of {1, 2, ..., 2n} into pairs whose differences are primes.
%e a(3) = 3: {{1,6}, {2,4}, {3,5}}, {{1,4}, {2,5}, {3,6}}, {{1,3}, {2,5}, {4,6}}. - _Alois P. Heinz_, Nov 15 2016
%p b:= proc(s) option remember; `if`(s={}, 1, (j-> add(`if`(i<j
%p and isprime(j-i), b(s minus {i, j}), 0), i=s))(max(s)))
%p end:
%p a:= n-> b({$1..2*n}):
%p seq(a(n), n=0..12); # _Alois P. Heinz_, Nov 15 2016
%t b[s_] := b[s] = If[s == {}, 1, Function[j, Sum[If[i < j && PrimeQ[j - i], b[s ~Complement~ {i, j}], 0], {i, s}]][Max[s]]];
%t a[n_] := b[Range[2n]];
%t Table[Print[n, " ", a[n]]; a[n], {n, 0, 18}] (* _Jean-François Alcover_, Mar 01 2021, after _Alois P. Heinz_ *)
%Y Cf. A000341.
%K nonn,hard
%O 0,4
%A _David W. Wilson_
%E a(0), a(14)-a(18) from _Alois P. Heinz_, Nov 15 2016
%E a(19)-a(23) from _Bert Dobbelaere_, Feb 20 2020
|