|
%I #17 Jul 02 2018 22:38:07
%S 0,0,0,0,0,0,720,5760,25920,86400,237600,570240,1235520,892045440,
%T 13348249200,106757164800,604924594560,2722120577280,10344007402560,
%U 34479959558400,24928970490633600,546446134633639680,6281586217487489040,50248618811434961280
%N Number of degree-n permutations of order exactly 7.
%H Alois P. Heinz, <a href="/A153760/b153760.txt">Table of n, a(n) for n = 1..200</a>
%F E.g.f.: -exp(x)+exp(x+1/7*x^7). - _Alois P. Heinz_, Feb 16 2013
%p a:= proc(n) option remember;
%p `if`(n<7, 0, a(n-1)+(1+a(n-7))*(n-1)!/(n-7)!)
%p end:
%p seq(a(n), n=1..30); # _Alois P. Heinz_, Feb 16 2013
%t Rest[CoefficientList[Series[-Exp[x] + Exp[x + 1/7*x^7], {x, 0, 25}], x]*Range[0, 25]!] (* _G. C. Greubel_, Aug 27 2016 *)
%o (PARI) seq(n)={Vec(serlaplace(-exp(x + O(x*x^n)) + exp(x + x^7/7 + O(x*x^n))),-n)} \\ _Andrew Howroyd_, Jul 02 2018
%Y Column k=7 of A057731. - _Alois P. Heinz_, Feb 16 2013
%K nonn
%O 1,7
%A _Franz Vrabec_, Jan 01 2009
%E More terms from _Alois P. Heinz_, Feb 16 2013
|
