cp's OEIS Frontend

A131631 Supersubfactorials: partial product of positive subfactorials (A000166).

%I A131631 #8 May 07 2020 06:13:01
%S A131631 1,2,18,792,209880,389117520,5771780174160,770509566129663360,
%T A131631 1028600220910021528728960,15104551945968674840127424147200,
%U A131631 2661646219535110627933754465838408595200
%N A131631 Supersubfactorials: partial product of positive subfactorials (A000166).
%C A131631 This is to subfactorials (A000166, rencontres numbers, or derangements) as superfactorials (A000178) are to factorials (A000142).
%F A131631 a(n) = Product_{k=2..n} A000166(k).
%F A131631 a(n) ~ c * n^(n^2/2 + n + 5/12) * (2*Pi)^((n+1)/2) / (A * exp(3*n^2/4 + 2*n - 13/12)), where A = A074962 = 1.2824271291... is the Glaisher-Kinkelin constant and c = 1.2517384488693662195086340541087053383189277225386098721341690164735... . - _Vaclav Kotesovec_, Jul 11 2015
%e A131631 a(2) = 1.
%e A131631 a(3) = 1 * 2 = 2.
%e A131631 a(4) = 1 * 2 * 9 = 18 = 2 * 3^2.
%e A131631 a(5) = 1 * 2 * 9 * 44 = 792 = 2^3 * 3^2 * 11.
%e A131631 a(6) = 1 * 2 * 9 * 44 * 265 = 209880 = 2^3 * 3^2 * 5 * 11 * 53.
%e A131631 a(7) = 1 * 2 * 9 * 44 * 265 * 1854 = 389117520 = 2^4 * 3^4 * 5 * 11 * 53 * 103.
%t A131631 Table[Product[k!*Sum[(-1)^j/j!, {j,0,k}], {k,2,n}], {n,2,15}] (* _Vaclav Kotesovec_, Jul 11 2015 *)
%Y A131631 Cf. A000142, A000166, A000178.
%K A131631 easy,nonn
%O A131631 2,2
%A A131631 _Jonathan Vos Post_, Sep 01 2007