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.

Previous Showing 11-12 of 12 results.

A246219 Number of endofunctions on [n] where the largest cycle length equals 9.

Original entry on oeis.org

40320, 4032000, 268329600, 15328051200, 823379356800, 43413497487360, 2298646350000000, 123818753182924800, 6837492347913427200, 388855261570122547200, 22836250136299660220160, 1386932177757615390720000, 87175183432121364413433600, 5671938409008942797114572800
Offset: 9

Views

Author

Alois P. Heinz, Aug 19 2014

Keywords

Crossrefs

Column k=9 of A241981.

Programs

  • Maple
    with(combinat):
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
          add((i-1)!^j*multinomial(n, n-i*j, i$j)/j!*
          b(n-i*j, i-1), j=0..n/i)))
        end:
    A:= (n, k)-> add(binomial(n-1, j-1)*n^(n-j)*b(j, min(j, k)), j=0..n):
    a:= n-> A(n, 9) -A(n, 8):
    seq(a(n), n=9..25);

Formula

a(n) ~ (9*exp(7129/2520) - 8*exp(761/280)) * n^(n-1). - Vaclav Kotesovec, Aug 21 2014

A246220 Number of endofunctions on [n] where the largest cycle length equals 10.

Original entry on oeis.org

362880, 43908480, 3448811520, 228012744960, 13954338478080, 827512686000000, 48753634065776640, 2895879112057451520, 174984885490926551040, 10817178515493080290560, 686533182382689959116800, 44833266187415969387604480, 3016487768851293040555130880
Offset: 10

Views

Author

Alois P. Heinz, Aug 19 2014

Keywords

Comments

In general, number of endofunctions on [n] where the largest cycle length equals k is asymptotic to (k*exp(H(k)) - (k-1)*exp(H(k-1))) * n^(n-1), where H(k) is the harmonic number A001008/A002805, k>=1. - Vaclav Kotesovec, Aug 21 2014

Crossrefs

Column k=10 of A241981.

Programs

  • Maple
    with(combinat):
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
          add((i-1)!^j*multinomial(n, n-i*j, i$j)/j!*
          b(n-i*j, i-1), j=0..n/i)))
        end:
    A:= (n, k)-> add(binomial(n-1, j-1)*n^(n-j)*b(j, min(j, k)), j=0..n):
    a:= n-> A(n, 10) -A(n, 9):
    seq(a(n), n=10..25);

Formula

a(n) ~ (10*exp(7381/2520) - 9*exp(7129/2520)) * n^(n-1). - Vaclav Kotesovec, Aug 21 2014
Previous Showing 11-12 of 12 results.