cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

A001189 Number of degree-n permutations of order exactly 2.

Original entry on oeis.org

0, 1, 3, 9, 25, 75, 231, 763, 2619, 9495, 35695, 140151, 568503, 2390479, 10349535, 46206735, 211799311, 997313823, 4809701439, 23758664095, 119952692895, 618884638911, 3257843882623, 17492190577599, 95680443760575, 532985208200575, 3020676745975551
Offset: 1

Views

Author

Keywords

Comments

Number of set partitions of [n] into blocks of size 2 and 1 with at least one block of size 2. - Olivier Gérard, Oct 29 2007
For n>=2, number of standard Young tableaux with height <= n - 1. That is, all tableaux (A000085) but the one with just one column. - Joerg Arndt, Oct 24 2012

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Column k=1 of A143911, column k=2 of A080510, A182222. - Alois P. Heinz, Oct 24 2012
Column k=2 of A057731. - Alois P. Heinz, Feb 14 2013

Programs

  • Magma
    m:=30; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Exp(x + x^2/2) -Exp(x) )); [0] cat [Factorial(n+1)*b[n]: n in [1..m-2]]; // G. C. Greubel, May 14 2019
    
  • Maple
    a:= proc(n) option remember; `if`(n<3, [0$2, 1][n+1],
          a(n-1) +(n-1) *(1+a(n-2)))
        end:
    seq(a(n), n=1..30);  # Alois P. Heinz, Oct 24 2012
    # alternative:
    A001189 := proc(n)
        local a,prs,p,k ;
        a := 0 ;
        for prs from 1 to n/2 do
            p := product(binomial(n-2*k,2),k=0..prs-1) ;
            a := a+p/prs!;
        end do:
        a;
    end proc:
    seq(A001189(n),n=1..13) ; # R. J. Mathar, Jan 04 2017
  • Mathematica
    RecurrenceTable[{a[1]==0,a[2]==1,a[n]==a[n-1]+(1+a[n-2])(n-1)},a[n],{n,25}] (* Harvey P. Dale, Jul 27 2011 *)
  • PARI
    {a(n) = sum(j=1,floor(n/2), n!/(j!*(n-2*j)!*2^j))}; \\ G. C. Greubel, May 14 2019
    
  • Sage
    m = 30; T = taylor(exp(x +x^2/2) - exp(x), x, 0, m); a=[factorial(n)*T.coefficient(x, n) for n in (0..m)]; a[1:] # G. C. Greubel, May 14 2019

Formula

E.g.f.: exp(x + x^2/2) - exp(x).
a(n) = A000085(n) - 1.
a(n) = b(n, 2), where b(n, d)=Sum_{k=1..n} (n-1)!/(n-k)! * Sum_{l:lcm{k, l}=d} b(n-k, l), b(0, 1)=1 is the number of degree-n permutations of order exactly d.
From Henry Bottomley, May 03 2001: (Start)
a(n) = a(n-1) + (1 + a(n-2))*(n-1).
a(n) = Sum_{j=1..floor(n/2)} n!/(j!*(n-2*j)!*(2^j)). (End)

A008307 Table T(n,k) giving number of permutations of [1..n] with order dividing k, read by antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 4, 1, 1, 1, 10, 3, 2, 1, 1, 26, 9, 4, 1, 1, 1, 76, 21, 16, 1, 2, 1, 1, 232, 81, 56, 1, 6, 1, 1, 1, 764, 351, 256, 25, 18, 1, 2, 1, 1, 2620, 1233, 1072, 145, 66, 1, 4, 1, 1, 1, 9496, 5769, 6224, 505, 396, 1, 16, 3, 2, 1, 1, 35696, 31041, 33616, 1345, 2052, 1, 56, 9, 4, 1, 1
Offset: 1

Views

Author

Keywords

Comments

Solutions to x^k = 1 in Symm_n (the symmetric group of degree n).

Examples

			Array begins:
  1,   1,    1,    1,    1,     1,    1,     1, ...
  1,   2,    1,    2,    1,     2,    1,     2, ...
  1,   4,    3,    4,    1,     6,    1,     4, ...
  1,  10,    9,   16,    1,    18,    1,    16, ...
  1,  26,   21,   56,   25,    66,    1,    56, ...
  1,  76,   81,  256,  145,   396,    1,   256, ...
  1, 232,  351, 1072,  505,  2052,  721,  1072, ...
  1, 764, 1233, 6224, 1345, 12636, 5761, 11264, ...
		

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 257.
  • J. D. Dixon, B. Mortimer, Permutation Groups, Springer (1996), Exercise 1.2.13.

Crossrefs

Programs

  • Maple
    A:= proc(n,k) option remember; `if`(n<0, 0, `if`(n=0, 1,
           add(mul(n-i, i=1..j-1)*A(n-j,k), j=numtheory[divisors](k))))
        end:
    seq(seq(A(1+d-k, k), k=1..d), d=1..12); # Alois P. Heinz, Feb 14 2013
    # alternative
    A008307 := proc(n,m)
        local x,d ;
        add(x^d/d, d=numtheory[divisors](m)) ;
        exp(%) ;
        coeftayl(%,x=0,n) ;
        %*n! ;
    end proc:
    seq(seq(A008307(1+d-k,k),k=1..d),d=1..12) ; # R. J. Mathar, Apr 30 2017
  • Mathematica
    t[n_ /; n >= 0, k_ /; k >= 0] := t[n, k] = Sum[(n!/(n - d + 1)!)*t[n - d, k], {d, Divisors[k]}]; t[, ] = 1; Flatten[ Table[ t[n - k, k], {n, 0, 12}, {k, 1, n}]] (* Jean-François Alcover, Dec 12 2011, after given formula *)

Formula

T(n+1,k) = Sum_{d|k} (n)_(d-1)*T(n-d+1,k), where (n)_i = n!/(n - i)! = n*(n - 1)*(n - 2)*...*(n - i + 1) is the falling factorial.
E.g.f. for n-th row: Sum_{n>=0} T(n,k)*t^n/n! = exp(Sum_{d|k} t^d/d).

Extensions

More terms from Vladeta Jovovic, Apr 13 2001

A005388 Number of degree-n permutations of order a power of 2.

Original entry on oeis.org

1, 1, 2, 4, 16, 56, 256, 1072, 11264, 78976, 672256, 4653056, 49810432, 433429504, 4448608256, 39221579776, 1914926104576, 29475151020032, 501759779405824, 6238907914387456, 120652091860975616, 1751735807564578816, 29062253310781161472, 398033706586943258624
Offset: 0

Views

Author

Keywords

Comments

Differs from A053503 first at n=32. - Alois P. Heinz, Feb 14 2013

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Example 5.2.10.

Crossrefs

Programs

  • Magma
    R:=PowerSeriesRing(Rationals(), 40);
    f:= func< x | Exp( (&+[x^(2^j)/2^j: j in [0..14]]) ) >;
    Coefficients(R!(Laplace( f(x) ))); // G. C. Greubel, Nov 17 2022
    
  • Maple
    a:= proc(n) option remember; `if`(n<0, 0, `if`(n=0, 1,
           add(mul(n-i, i=1..2^j-1)*a(n-2^j), j=0..ilog2(n))))
        end:
    seq(a(n), n=0..25);  # Alois P. Heinz, Feb 14 2013
  • Mathematica
    max = 23; CoefficientList[ Series[ Exp[ Sum[x^2^m/2^m, {m, 0, max}]], {x, 0, max}], x]*Range[0, max]! (* Jean-François Alcover, Sep 10 2013 *)
  • SageMath
    def f(x): return exp(sum(x^(2^j)/2^j for j in range(15)))
    def A005388_list(prec):
        P. = PowerSeriesRing(QQ, prec)
        return P( f(x) ).egf_to_ogf().list()
    A005388_list(40) # G. C. Greubel, Nov 17 2022

Formula

E.g.f.: exp(Sum_{m>=0} x^(2^m)/2^m).
E.g.f.: 1/Product_{k>=1} (1 - x^(2*k-1))^(mu(2*k-1)/(2*k-1)), where mu() is the Moebius function. - Seiichi Manyama, Jul 06 2024
Showing 1-3 of 3 results.