Search: id:A000254 Results 1-1 of 1 results found. %I A000254 M2902 N1165 %S A000254 0,1,3,11,50,274,1764,13068,109584,1026576,10628640,120543840, %T A000254 1486442880,19802759040,283465647360,4339163001600,70734282393600, %U A000254 1223405590579200,22376988058521600,431565146817638400 %N A000254 Unsigned Stirling numbers of first kind, s(n+1,2): a(n+1)=(n+1)*a(n)+n!. %C A000254 Number of permutations of n elements with exactly two cycles. %C A000254 a(n)=number of cycles in all permutations of [n]. Example: a(3)=11 because the permutations (1)(2)(3), (1)(23), (12)(3), (13)(2), (132), (123) have 11 cycles alltogether. - Emeric Deutsch (deutsch(AT)duke.poly.edu), Aug 12 2004 %C A000254 The sum of the top levels of the last column over all deco polyominoes of height n. A deco polyomino is a directed column-convex polyomino in which the height, measured along the diagonal, is attained only in the last column. Example: a(2)=3 because the deco polyominoes of height 2 are the vertical and horizontal dominoes, the levels of their last columns being 2 and 1, respectively. - Emeric Deutsch (deutsch(AT)duke.poly.edu), Aug 12 2006 %C A000254 a(n) is divisible by n for all composite n >= 6. a(2n) is divisible by (2n+1). - Leroy Quet May 20 2007 %C A000254 For n >= 2 the determinant of the n-1 X n-1 matrix M(i,j) = i + 2 for i = j and 1 otherwise (i,j = 1..n-1). E.g. for n = 3 the determinant of [(3, 1), (1, 4)]. See 53rd Putnam Examination, 1992, Problem B5. - Franz Vrabec (franz.vrabec(AT)planetuniqa.at), Jan 13 2008, Mar 26 2008 %C A000254 Contribution from Eric Desbiaux (moongerms(AT)wanadoo.fr), Jan 07 2009: (Start) %C A000254 The numerator of the fraction when we sum (without simplification) the terms in the harmonic sequence. %C A000254 (1+1/2=2/2+1/2= 3/2; 3/2+1/3=9/6+2/6= 11/6; 11/6+1/4=44/24+6/24= 50/24; ...) %C A000254 The denominator of this fraction is n! A000142. %C A000254 (End) %C A000254 Contribution from Johannes W. Meijer (meijgia(AT)hotmail.com), Oct 20 2009: (Start) %C A000254 The asymptotic expansion of the higher order exponential integral E(x, m=2,n=1) ~ exp(-x)/x^2*(1 - 3/x + 11/x^2 - 50/x^3 + 274/x^4 - 1764/ x^5 + 13068/x^6 - ...) leads to the sequence given above. See A163931 and A028421 for more information. %C A000254 (End) %C A000254 Contribution from Tom Woodward (twoodward(AT)macalester.edu), Nov 12 2009: (Start) %C A000254 a(n)=number of permutations of [n+1] containing exactly 2 cycles. Example: a(2)=3 %C A000254 because the permutations (1)(23), (12)(3), (13)(2) are the only permutations of [3] %C A000254 with exactly 2 cycles. (End) %D A000254 M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 833. %D A000254 E. Barcucci, A. Del Lungo and R. Pinzani, "Deco" polyominoes, permutations and random generation, Theoretical Computer Science, 159, 1996, 29-42. %D A000254 A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, identities 186-190. %D A000254 N. Bleistein and R. A. Handelsman, Asymptotic Expansions of Integrals, Dover Publications, 1986, see page 2. MR0863284 (89d:41049) %D A000254 L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 217. %D A000254 F. N. David, M. G. Kendall and D. E. Barton, Symmetric Function and Allied Tables, Cambridge, 1966, p. 226. %D A000254 N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence). %D A000254 N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence). %H A000254 T. D. Noe, Table of n, a(n) for n=0..100 %H A000254 Leroy Quet, Home Page (listed in lieu of email address) %H A000254 M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy]. %H A000254 INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 31 %H A000254 J. Scholes, 53rd Putnam 1992, Problem B5. %F A000254 Let P(n,X)=(X+1)(X+2)(X+3)...(X+n); then a(n) is the coefficient of X; or a(n)=P'(n,0) - Benoit Cloitre (benoit7848c(AT)orange.fr), May 09 2002 %F A000254 Sum_{n>0} a(n)x^n/n!^2 = exp(x)(Sum_{n>0}(-1)^(n+1)x^n/(n*n!)). - Michael Somos Mar 24 2004. Corrected by Warren Smith, Feb 12 2006. %F A000254 a(n) is coefficient of x^(n+2) in (-log(1-x))^2, multiplied by (n+2)!/ 2. %F A000254 Also a(n) = n!*Sum 1/i, i=1..n = n!*H(n), H(n) = harmonic number = A001008/ A002805. %F A000254 a(n) ~ 2^(1/2)*pi^(1/2)*log(n)*n^(1/2)*e^-n*n^n - Joe Keane (jgk(AT)jgk.org), Jun 06 2002 %F A000254 E.g.f.: log(1-x)/(x-1). (= (log(1-x))^2/2 if offset 1). - Michael Somos Feb 05 2004 %F A000254 a(n)=a(n-1)(2n-1)-a(n-2)(n-1)^2, if n>1. - Michael Somos Mar 24 2004 %F A000254 a(n)=A081358(n)+A092691(n). - Emeric Deutsch (deutsch(AT)duke.poly.edu), Aug 12 2004 %F A000254 a(n) = n!*Sqrt[Sum[Sum[1/(i*j), {i, 1, n}], {j, 1, n}]] - Alexander Adamchuk (alex(AT)kolmogorov.com), Oct 26 2004 %F A000254 a(n) = n!*Sum_{k=1..n} (-1)^(k+1)*binomial(n, k)/k. - Vladeta Jovovic (vladeta(AT)eunet.rs), Jan 29 2005 %F A000254 p^2 divides a(p-1) for prime p>3. a(n) = Sum[ 1/i, {i,1,n}] / Product[ 1/i, {i,1,n}]. - Alexander Adamchuk (alex(AT)kolmogorov.com), Jul 11 2006 %F A000254 For n>=1, a(n)=n!*sum(1/(n-k),k=0..n-1); [From Milan R. Janjic (agnus(AT)blic.net), Dec 14 2008] %e A000254 (1-x)^-1 * (-log(1-x)) = x + 3/2*x^2 + 11/6*x^3 + 25/12*x^4 + ... %p A000254 A000254 := proc(n) option remember; if n<=1 then 1 else n*A000254(n-1)+(n-1)!; fi; end; %p A000254 a:=n->sum(n!/k, k=1..n): seq(a(n), n=0..19); - Zerinvary Lajos (zerinvarylajos(AT)yahoo.com), Jan 22 2008 %t A000254 Table[ (PolyGamma[ m ]+EulerGamma) (m-1)!, {m, 1, 24} ] - from W. Meeussen (wouter.meeussen(AT)pandora.be) %t A000254 Table[ n!*HarmonicNumber[n], {n, 0, 19}] (from Robert G. Wilson v (rgwv(AT)rgwv.com), May 21 2005) %t A000254 Table[Sum[1/i,{i,1,n}]/Product[1/i,{i,1,n}],{n,1,30}] - Alexander Adamchuk (alex(AT)kolmogorov.com), Jul 11 2006 %o A000254 (MuPAD) A000254 := proc(n) begin n*A000254(n-1)+fact(n-1) end_proc: A000254(1) := 1: %o A000254 (PARI) a(n)=if(n<0,0,(n+1)!/2*sum(k=1,n,1/k/(n+1-k))) %o A000254 sage: [stirling_number1(i,2) for i in xrange(1,22)] - Zerinvary Lajos (zerinvarylajos(AT)yahoo.com), Jun 27 2008 %Y A000254 Cf. A000399, A000774, A004041, A024167, A046674, A049034, A008275 (Stirling1 triangle). %Y A000254 Cf. A081358, A092691, A151881. %Y A000254 With signs: A081048. %Y A000254 Column 1 in triangle A008969. %Y A000254 Cf. A121633. %Y A000254 Sequence in context: A162477 A115081 A103466 this_sequence A081048 A065048 A024335 %Y A000254 Adjacent sequences: A000251 A000252 A000253 this_sequence A000255 A000256 A000257 %K A000254 nonn,easy,nice,new %O A000254 0,3 %A A000254 N. J. A. Sloane (njas(AT)research.att.com). Search completed in 0.003 seconds