|
%I #22 Aug 10 2019 09:00:29
%S 1,1,5,31,217,1661,13721,121703,1157857,11826121,129877645,1535504015,
%T 19546846441,267633414517,3932330905361,61806788736551,
%U 1035452546213441,18421374554192017,346790652640704725,6885640002624595007,143771244649798115257
%N Number of permutations of [2n] with exactly n cycles such that the elements of each cycle form an integer interval.
%C All terms are odd.
%H Alois P. Heinz, <a href="/A287899/b287899.txt">Table of n, a(n) for n = 0..448</a>
%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Permutation">Permutation</a>
%F a(n) = A084938(2n,n).
%F a(n) = [x^n] (1/(1 - x/(1 - x/(1 - 2*x/(1 - 2*x/(1 - 3*x/(1 - 3*x/(1 - ...))))))))^n, a continued fraction. - _Ilya Gutkovskiy_, Sep 29 2017
%F a(n) ~ exp(1) * n * n!. - _Vaclav Kotesovec_, Sep 29 2017
%e a(2) = 5: (1)(2,3,4), (1)(2,4,3), (1,2)(3,4), (1,2,3)(4), (1,3,2)(4).
%p b:= proc(n, i) option remember; `if`(n=0 or i=1, n!,
%p add(b(n-j, i-1)*j!, j=0..n))
%p end:
%p a:= n-> b(n$2):
%p seq(a(n), n=0..25);
%t Table[SeriesCoefficient[1/(1 + ContinuedFractionK[-Floor[(k + 1)/2]*x, 1, {k, 1, n}])^n, {x, 0, n}], {n, 0, 20}] (* _Vaclav Kotesovec_, Sep 29 2017 *)
%t Table[SeriesCoefficient[Sum[k!*x^k, {k, 0, n}]^n, {x,0,n}], {n,0,25}] (* _Vaclav Kotesovec_, Aug 10 2019 *)
%Y Cf. A084938, A088218 (analog for set partitions).
%K nonn
%O 0,3
%A _Alois P. Heinz_, Jun 02 2017
|
