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-7 of 7 results.

A303694 Array read by antidiagonals: T(n,k) is the number of noncrossing partitions up to rotation composed of n blocks of size k.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 2, 3, 1, 1, 1, 1, 3, 7, 6, 1, 1, 1, 1, 3, 11, 19, 14, 1, 1, 1, 1, 4, 17, 52, 86, 34, 1, 1, 1, 1, 4, 25, 102, 307, 372, 95, 1, 1, 1, 1, 5, 33, 187, 811, 1936, 1825, 280, 1, 1, 1, 1, 5, 43, 300, 1772, 6626, 13207, 9143, 854, 1
Offset: 0

Views

Author

Andrew Howroyd, Apr 28 2018

Keywords

Comments

Also, the number of unlabeled planar k-gonal cacti having n polygons.
The number of noncrossing partitions counted distinctly is given by A070914(n,k-1).

Examples

			Array begins:
==================================================================
n\k| 1   2    3     4      5       6       7        8        9
---+--------------------------------------------------------------
0  | 1   1    1     1      1       1       1        1        1 ...
1  | 1   1    1     1      1       1       1        1        1 ...
2  | 1   1    1     1      1       1       1        1        1 ...
3  | 1   2    2     3      3       4       4        5        5 ...
4  | 1   3    7    11     17      25      33       43       55 ...
5  | 1   6   19    52    102     187     300      463      663 ...
6  | 1  14   86   307    811    1772    3412     5993     9821 ...
7  | 1  34  372  1936   6626   17880   40770    82887   154079 ...
8  | 1  95 1825 13207  58385  191967  518043  1213879  2558305 ...
9  | 1 280 9143 93496 532251 2141232 6830545 18471584 44121134 ...
...

Links

Crossrefs

Columns 2..7 are A002995(n+1), A054423, A054362, A054365, A054368, A054371.
Cf. A070914, A209805, A303864, A303929.

Programs

  • Mathematica
    T[0, _] = 1;
T[n_, k_] := (DivisorSum[n, EulerPhi[n/#] Binomial[k #, #]&] + DivisorSum[ GCD[n-1, k], EulerPhi[#] Binomial[n k/#, (n-1)/#]&])/(k n) - Binomial[k n, n]/(n (k-1) + 1);
Table[T[n-k, k], {n, 0, 12}, {k, n, 1, -1}] // Flatten (* Jean-François Alcover, May 22 2018 *)
  • PARI
    T(n,k)={if(n==0, 1, (sumdiv(n,d,eulerphi(n/d)*binomial(k*d,d)) + sumdiv(gcd(n-1,k), d, eulerphi(d)*binomial(n*k/d, (n-1)/d)))/(k*n) - binomial(k*n,n)/(n*(k-1)+1))}

Formula

T(n,k) = ((Sum_{d|n} phi(n/d)*binomial(k*d,d)) + (Sum_{d|gcd(n-1,k)} phi(d) * binomial(n*k/d, (n-1)/d)))/(k*n) - binomial(k*n,n)/(n*(k-1)+1) for n > 0.
T(n,k) ~ A070914(n,k-1)/(n*k) for fixed k > 1.

A006082 Number of labeled projective plane trees (or "flat" trees) with n nodes.

Original entry on oeis.org

1, 1, 1, 2, 3, 6, 12, 27, 65, 175, 490, 1473, 4588, 14782, 48678, 163414, 555885, 1913334, 6646728, 23278989, 82100014, 291361744, 1039758962, 3729276257, 13437206032, 48620868106, 176611864312, 643834562075, 2354902813742, 8640039835974, 31791594259244
Offset: 1

Views

Author

N. J. A. Sloane

Keywords

Comments

Also, the number of noncrossing partitions up to rotation and reflection composed of n-1 blocks of size 2. - Andrew Howroyd, May 03 2018

References

  • R. W. Robinson, personal communication.
  • R. W. Robinson, Efficiency of power series operations for graph counting, AGRC Grant, Math. Dept., Univ. Newcastle, Australia, 1982.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Links

Crossrefs

Column k=2 of A302828 and A303929.
Cf. A006079, A006080, A006081.
Cf. A002995 (noncrossing partitions into pairs up to rotations only), A126120, A001405, A185100.

Programs

  • Mathematica
    u[n_, k_, r_] := (r*Binomial[k*n + r, n]/(k*n + r));
e[n_, k_] := Sum[ u[j, k, 1 + (n - 2*j)*k/2], {j, 0, n/2}]
c[n_, k_] := If[n == 0, 1, (DivisorSum[n, EulerPhi[n/#]*Binomial[k*#, #]&] + DivisorSum[GCD[n-1, k], EulerPhi[#]*Binomial[n*k/#, (n-1)/#]&])/(k*n) - Binomial[k*n, n]/(n*(k - 1) + 1)];
T[n_, k_] := (1/2)*(c[n, k] + If[n == 0, 1, If[OddQ[k], If[OddQ[n], 2*u[ Quotient[n, 2], k, (k + 1)/2], u[n/2, k, 1] + u[n/2 - 1, k, k]], e[n, k] + If[OddQ[n], u[Quotient[n, 2], k, k/2]]]/2]) /. Null -> 0;
a[n_] := T[n, 2];
Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Jun 14 2018, after Andrew Howroyd and A303929 *)
  • PARI
    \\ from David Broadhurst, Apr 06 2022, added by N. J. A. Sloane, Apr 06 2022
{A006082(n)=my(c(n)=binomial(2*n,n));
if(n<2,1,n--;(c(n)+if(n%2,2*n*(n+2),(n+1)^2)*c(n\2)
+(n+1)*sumdiv(n,d,if(d>2,eulerphi(d)*c(n/d))))/(4*n*(n+1)));}

Formula

a(n) = A006080(n) - A006081(n) + A126120(n-2). [Stockmeyer] [Corrected by Andrey Zabolotskiy, Apr 06 2021]
a(n) = (2 * A002995(n) + A126120(n-2) + A001405(n-1)) / 4 for n > 1. - Andrey Zabolotskiy, May 24 2018
There is a compact formula from David Broadhurst - see the Pari code - N. J. A. Sloane, Apr 06 2022.
a(n) ~ 2^(2*n-4) / (sqrt(Pi) * n^(5/2)). - Vaclav Kotesovec, Jun 01 2022

Extensions

a(25) and a(26) from Robert W. Robinson, Oct 17 2006
a(27) and beyond from Andrew Howroyd, May 03 2018

A369929 Array read by antidiagonals: T(n,k) is the number of achiral noncrossing partitions composed of n blocks of size k.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 2, 3, 1, 1, 1, 1, 3, 3, 6, 1, 1, 1, 1, 3, 5, 7, 10, 1, 1, 1, 1, 4, 5, 16, 12, 20, 1, 1, 1, 1, 4, 7, 18, 31, 30, 35, 1, 1, 1, 1, 5, 7, 31, 35, 102, 55, 70, 1, 1, 1, 1, 5, 9, 34, 64, 136, 213, 143, 126, 1
Offset: 0

Views

Author

Andrew Howroyd, Feb 07 2024

Keywords

Comments

T(n,2*k-1) is the number of achiral noncrossing k-gonal cacti with n polygons.

Examples

			Array begins:
===============================================
n\k| 1  2   3   4    5    6    7    8     9 ...
---+-------------------------------------------
0  | 1  1   1   1    1    1    1    1     1 ...
1  | 1  1   1   1    1    1    1    1     1 ...
2  | 1  1   1   1    1    1    1    1     1 ...
3  | 1  2   2   3    3    4    4    5     5 ...
4  | 1  3   3   5    5    7    7    9     9 ...
5  | 1  6   7  16   18   31   34   51    55 ...
6  | 1 10  12  31   35   64   70  109   117 ...
7  | 1 20  30 102  136  296  368  651   775 ...
8  | 1 35  55 213  285  663  819 1513  1785 ...
9  | 1 70 143 712 1155 3142 4495 9304 12350 ...
...

Links

Crossrefs

Columns are: A000012, A001405(n-1), A047749 (k=3), A369930 (k=4), A143546 (k=5), A143547 (k=7), A143554 (k=9), A192893 (k=11).
Cf. A070914, A303694, A303929, A361236, A361239, A370060, A370062.

Programs

  • PARI
    \\ u(n,k,r) are Fuss-Catalan numbers.
u(n,k,r) = {r*binomial(k*n + r, n)/(k*n + r)}
e(n,k) = {sum(j=0, n\2, u(j, k, 1+(n-2*j)*k/2))}
T(n, k)={if(n==0, 1, if(k%2, if(n%2, 2*u(n\2, k, (k+1)/2), u(n/2, k, 1) + u(n/2-1, k, k)), e(n, k) + if(n%2, u(n\2, k, k/2)))/2)}

Formula

T(n,k) = 2*A303929(n,k) - A303694(n,k).
T(n,2*k-1) = 2*A361239(n,k) - A361236(n,k).

A302828 Array read by antidiagonals: T(n,k) = number of noncrossing path sets on k*n nodes up to rotation and reflection with each path having exactly k nodes.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 4, 2, 1, 1, 3, 21, 22, 3, 1, 1, 6, 111, 494, 201, 6, 1, 1, 10, 604, 9400, 18086, 2244, 12, 1, 1, 20, 3196, 157040, 1141055, 794696, 29096, 27, 1, 1, 36, 16528, 2342480, 55967596, 161927208, 38695548, 404064, 65, 1
Offset: 0

Views

Author

Andrew Howroyd, May 01 2018

Keywords

Examples

			Array begins:
=======================================================
n\k| 1  2     3        4           5              6
---+---------------------------------------------------
0  | 1  1     1        1           1              1 ...
1  | 1  1     1        2           3              6 ...
2  | 1  1     4       21         111            604 ...
3  | 1  2    22      494        9400         157040 ...
4  | 1  3   201    18086     1141055       55967596 ...
5  | 1  6  2244   794696   161927208    23276467936 ...
6  | 1 12 29096 38695548 25334545270 10673231900808 ...
...

Links

Crossrefs

Columns 2..4 are A006082(n+1), A303330, A303867.
Row n=1 is A005418(k-2).
Cf. A303839, A303864, A303868, A303929.

Programs

  • Mathematica
    nmax = 10; seq[n_, k_] := Module[{p, q, h, c}, p = 1 + InverseSeries[ x/(k*2^(k - 3)*(1 + x)^k) + O[x]^n, x]; h = p /. x -> x^2 + O[x]^n; q = x*D[p, x]/p; c = Integrate[((p - 1)/k + Sum[EulerPhi[d]*(q /. x -> x^d + O[x]^n), {d, 2, n}])/x, x] + If[OddQ[k], 0, 2^(k/2 - 2)*x*h^(k/2)]; If[k == 1, 2/(1 - x) + O[x]^n, 3/2 + c + If[OddQ[k], h + x^2*2^(k - 3)*h^k + x*2^((k - 1)/2)*h^((k + 1)/2), If[k == 2, x*h, 0] + h/(1 - 2^(k/2 - 1)*x*h^(k/2))]/2]/2];
Clear[col]; col[k_] := col[k] = CoefficientList[seq[nmax, k], x];
T[n_, k_] := col[k][[n + 1]];
Table[T[n - k, k], {n, 0, nmax}, {k, n, 1, -1}] // Flatten (* Jean-François Alcover, Jul 04 2018, after Andrew Howroyd *)
  • PARI
    seq(n,k)={ \\ gives gf of k'th column
my(p=1 + serreverse( x/(k*2^(k-3)*(1 + x)^k) + O(x*x^n) ));
my(h=subst(p,x,x^2+O(x*x^n)), q=x*deriv(p)/p);
my(c=intformal( ((p-1)/k + sum(d=2,n,eulerphi(d)*subst(q,x,x^d+O(x*x^n))))/x) + if(k%2, 0, 2^(k/2-2)*x*h^(k/2)));
if(k==1, 2/(1-x) + O(x*x^n), 3/2 + c + if(k%2, h + x^2*2^(k-3)*h^k + x*2^((k-1)/2)*h^((k+1)/2), if(k==2,x*h,0) + h/(1-2^(k/2-1)*x*h^(k/2)) )/2)/2;
}
Mat(vector(6, k, Col(seq(7, k))))

A082938 Number of solid 2-trees with 2n+1 edges.

Original entry on oeis.org

1, 1, 1, 2, 5, 13, 49, 201, 940, 4643, 24037, 127859, 696365, 3858759, 21704863, 123619126, 711787259, 4137614454, 24256010068, 143271593982, 852001881614, 5097719884665, 30670572676389, 185466705697057
Offset: 0

Views

Author

N. J. A. Sloane, May 26 2003

Keywords

Comments

Also, the number of noncrossing partitions up to rotation and reflection composed of n blocks of size 3. - Andrew Howroyd, May 03 2018

Links

Crossrefs

Column k=3 of A303929.
Cf. A047749, A054423.

Programs

  • Mathematica
    u[n_, k_, r_] := (r*Binomial[k*n + r, n]/(k*n + r));
e[n_, k_] := Sum[ u[j, k, 1 + (n - 2*j)*k/2], {j, 0, n/2}]
c[n_, k_] := If[n == 0, 1, (DivisorSum[n, EulerPhi[n/#]*Binomial[k*#, #] &] + DivisorSum[GCD[n-1, k], EulerPhi[#]*Binomial[n*k/#, (n-1)/#] &])/(k*n) - Binomial[k*n, n]/(n*(k - 1) + 1)];
T[n_, k_] := (1/2)*(c[n, k] + If[n == 0, 1, If[OddQ[k], If[OddQ[n], 2*u[ Quotient[n, 2], k, (k + 1)/2], u[n/2, k, 1] + u[n/2 - 1, k, k]], e[n, k] + If[OddQ[n], u[Quotient[n, 2], k, k/2]]]/2]) /. Null -> 0;
a[n_] := T[n, 3];
Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Jun 14 2018, after Andrew Howroyd and A303929 *)

Formula

a(n) = (A047749(n)+A054423(n))/2. - Vladeta Jovovic, Sep 11 2004
a(n) ~ 3^(3*n - 1/2) / (sqrt(Pi) * n^(5/2) * 2^(2*n + 3)). - Vaclav Kotesovec, Jun 01 2022

Extensions

More terms from Vladeta Jovovic, Sep 11 2004

A303870 Number of noncrossing partitions up to rotation and reflection composed of n blocks of size 4.

Original entry on oeis.org

1, 1, 1, 3, 8, 34, 169, 1019, 6710, 47104, 342772, 2566209, 19621256, 152669854, 1205358482, 9636786366, 77890590994, 635628049370, 5231328157060, 43382605871299, 362225044991368, 3043083681629249, 25708398651274529, 218296978274674435, 1862280135781609982
Offset: 0

Views

Author

Andrew Howroyd, May 01 2018

Keywords

Links

Crossrefs

Column k=4 of A303929.
Cf. A054362.

Programs

  • Mathematica
    u[n_, k_, r_] := (r*Binomial[k*n + r, n]/(k*n + r));
e[n_, k_] := Sum[ u[j, k, 1 + (n - 2*j)*k/2], {j, 0, n/2}]
c[n_, k_] := If[n == 0, 1, (DivisorSum[n, EulerPhi[n/#]*Binomial[k*#, #] &] + DivisorSum[GCD[n - 1, k], EulerPhi[#]*Binomial[n*k/#, (n - 1)/#] &])/(k*n) - Binomial[k*n, n]/(n*(k - 1) + 1)];
T[n_, k_] := (1/2)*(c[n, k] + If[n == 0, 1, If[OddQ[k], If[OddQ[n], 2*u[Quotient[n, 2], k, (k + 1)/2], u[n/2, k, 1] + u[n/2 - 1, k, k]], e[n, k] + If[OddQ[n], u[Quotient[n, 2], k, k/2]]]/2]) /. Null -> 0;
a[n_] := T[n, 4];
Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Jun 14 2018, after _Andrew Howroyd and A303929 *)

Formula

a(n) ~ 2^(8*n - 5/2) / (sqrt(Pi) * n^(5/2) * 3^(3*n + 3/2)). - Vaclav Kotesovec, Jun 01 2022

A303871 Number of noncrossing partitions up to rotation and reflection composed of n blocks of size 5.

Original entry on oeis.org

1, 1, 1, 3, 11, 60, 423, 3381, 29335, 266703, 2507232, 24177705, 238003111, 2383370158, 24217426745, 249182213284, 2592138293117, 27225668134063, 288405507217589, 3078471666603235, 33085393411436772, 357782389095170193, 3890765813426578535, 42527471172438573757
Offset: 0

Views

Author

Andrew Howroyd, May 01 2018

Keywords

Links

Crossrefs

Column k=5 of A303929.
Cf. A054365.

Programs

  • Mathematica
    u[n_, k_, r_] := (r*Binomial[k*n + r, n]/(k*n + r));
e[n_, k_] := Sum[ u[j, k, 1 + (n - 2*j)*k/2], {j, 0, n/2}]
c[n_, k_] := If[n == 0, 1, (DivisorSum[n, EulerPhi[n/#]*Binomial[k*#, #] &] + DivisorSum[GCD[n - 1, k], EulerPhi[#]*Binomial[n*k/#, (n - 1)/#] &])/(k*n) - Binomial[k*n, n]/(n*(k - 1) + 1)];
T[n_, k_] := (1/2)*(c[n, k] + If[n == 0, 1, If[OddQ[k], If[OddQ[n], 2*u[Quotient[n, 2], k, (k + 1)/2], u[n/2, k, 1] + u[n/2 - 1, k, k]], e[n, k] + If[OddQ[n], u[Quotient[n, 2], k, k/2]]]/2]) /. Null -> 0;
a[n_] := T[n, 5];
Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Jun 14 2018, after Andrew Howroyd and A303929 *)

Formula

a(n) ~ 5^(5*n - 1/2) / (sqrt(Pi) * n^(5/2) * 2^(8*n + 9/2)). - Vaclav Kotesovec, Jun 01 2022
Showing 1-7 of 7 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.