a[n_] := (2*n-1)!! * Hypergeometric1F1[-n, 1-2*n, -2]; a[1] = 0;
twounifll[{}]:={{}}; twounifll[set:{i_, ___}]:=Join@@Function[s, Prepend[#, s]&/@twounifll[Complement[set, s]]]/@Table[{i, j}, {j, If[i==1, Select[set, 2<#<Last[set]&], Select[set, #>i+1&]]}];
Table[Length[twounifll[Range[n]]], {n, 0, 14, 2}] (*
Gus Wiseman, Feb 27 2019 *)