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.

A087193 H(p)/p where p runs through the primes and H(k) is the k-th central hexanomial coefficient (A018901).

Original entry on oeis.org

3, 9, 156, 3431, 2280921, 64092034, 55747783602, 1700642242677, 1657887524047959, 54732141299289779730, 1783584256683646551447, 63884853139612229737722392, 71016623651822742997810429944, 2380864745882038026563515929162, 2701273375177028344436110369387929
Offset: 1

Views

Author

Benoit Cloitre, Oct 19 2003

Keywords

Crossrefs

Cf. A000040 (primes), A018901 (central hexanomials).

Programs

  • Mathematica
    f[n_] := Max[CoefficientList[Expand[Sum[x^k, {k, 0, 5}]^n], x]]; Table[f[p]/p, {p, Prime[Range[15]]}] (* Amiram Eldar, Apr 25 2025 *)
  • PARI
    \\ here b(n) is A018901.
    b(n) = {if(n==0, 1, sum(k=0, 5*n\12, (-1)^k*binomial(n,k)*binomial(n + 5*n\2 - 6*k - 1, n - 1)))}
    a(n) = {my(p=prime(n)); b(p)/p} \\ Andrew Howroyd, Jan 08 2020

Formula

a(n) = A018901(prime(n))/prime(n). - Andrew Howroyd, Jan 08 2020

Extensions

Terms a(13) and beyond from Andrew Howroyd, Jan 08 2020

A063260 Sextinomial (also called hexanomial) coefficient array.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 2, 3, 4, 5, 6, 5, 4, 3, 2, 1, 1, 3, 6, 10, 15, 21, 25, 27, 27, 25, 21, 15, 10, 6, 3, 1, 1, 4, 10, 20, 35, 56, 80, 104, 125, 140, 146, 140, 125, 104, 80, 56, 35, 20, 10, 4, 1, 1, 5, 15, 35, 70, 126, 205, 305, 420, 540, 651, 735, 780
Offset: 0

Views

Author

Wolfdieter Lang, Jul 24 2001

Keywords

Comments

The sequence of step width of this staircase array is [1,5,5,...], hence the degree sequence for the row polynomials is [0,5,10,15,...]=A008587.
The column sequences (without leading zeros) are for k=0..5 those of the lower triangular array A007318 (Pascal) and for k=6..9: A062989, A063262-4. Row sums give A000400 (powers of 6). Central coefficients give A063419; see also A018901.
This can be used to calculate the number of occurrences of a given roll of n six-sided dice, where k is the index: k=0 being the lowest possible roll (i.e., n) and n*6 being the highest roll.

Examples

			The irregular table T(n, k) begins:
n\k 0 1 2  3  4  5  6  7  8  9 10 11 12 13 14 15
1:  1
2:  1 1 1  1  1  1
3:  1 2 3  4  5  6  5  4  3  2  1
4:  1 3 6 10 15 21 25 27 27 25 21 15 10  6  3  1
...reformatted - _Wolfdieter Lang_, Oct 31 2015
		

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, pp. 77,78.

Crossrefs

The q-nomial arrays for q=2..5 are: A007318 (Pascal), A027907, A008287, A035343 and for q=7: A063265, A171890, A213652, A213651.
Columns for k=0..9 (with some shifts) are: A000012, A000027, A000217, A000292, A000332, A000389, A062989, A063262, A063263, A063264.

Programs

  • Maple
    #Define the r-nomial coefficients for r = 1, 2, 3, ...
    rnomial := (r,n,k) -> add((-1)^i*binomial(n,i)*binomial(n+k-1-r*i,n-1), i = 0..floor(k/r)):
    #Display the 6-nomials as a table
    r := 6:  rows := 10:
    for n from 0 to rows do
    seq(rnomial(r,n,k), k = 0..(r-1)*n)
    end do;
    # Peter Bala, Sep 07 2013
  • Mathematica
    Flatten[Table[CoefficientList[(1 + x + x^2 + x^3 + x^4 + x^5)^n, x], {n, 0, 25}]] (* T. D. Noe, Apr 04 2011 *)
  • PARI
    concat(vector(5,k,Vec(sum(j=0,5,x^j)^k)))  \\ M. F. Hasler, Jun 17 2012

Formula

G.f. for row n: (Sum_{j=0..5} x^j)^n.
G.f. for column k: (x^(ceiling(k/5)))*N6(k, x)/(1-x)^(k+1) with the row polynomials from the staircase array A063261(k, m) and with N6(6,x) = 5 - 10*x + 10*x^2 - 5*x^3 + x^4.
T(n, k) = 0 if n=-1 or k<0 or k >= 5*n + 1; T(0, 0)=1; T(n, k) = Sum_{j=0..5} T(n-1, k-j) else.
T(n, k) = Sum_{i = 0..floor(k/6)} (-1)^i*binomial(n,i)*binomial(n+k-1-6*i,n-1) for n >= 0 and 0 <= k <= 5*n. - Peter Bala, Sep 07 2013
T(n, k) = Sum_{i = max(0,ceiling((k-2*n)/3)).. min(n,k/3)} binomial(n,i)*trinomial(n,k-3*i) for n >= 0 and 0 <= k <= 5*n. - Matthew Monaghan, Sep 30 2015

Extensions

More terms and corrected recurrence from Nicholas M. Makin (NickDMax(AT)yahoo.com), Sep 13 2002

A077042 Square array read by falling antidiagonals of central polynomial coefficients: largest coefficient in expansion of (1 + x + x^2 + ... + x^(n-1))^k = ((1-x^n)/(1-x))^k, i.e., the coefficient of x^floor(k*(n-1)/2) and of x^ceiling(k*(n-1)/2); also number of compositions of floor(k*(n+1)/2) into exactly k positive integers each no more than n.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 2, 1, 1, 0, 1, 3, 3, 1, 1, 0, 1, 6, 7, 4, 1, 1, 0, 1, 10, 19, 12, 5, 1, 1, 0, 1, 20, 51, 44, 19, 6, 1, 1, 0, 1, 35, 141, 155, 85, 27, 7, 1, 1, 0, 1, 70, 393, 580, 381, 146, 37, 8, 1, 1, 0, 1, 126, 1107, 2128, 1751, 780, 231, 48, 9, 1, 1, 0, 1, 252, 3139
Offset: 0

Views

Author

Henry Bottomley, Oct 22 2002

Keywords

Comments

From Michel Marcus, Dec 01 2012: (Start)
A pair of numbers written in base n are said to be comparable if all digits of the first number are at least as big as the corresponding digit of the second number, or vice versa. Otherwise, this pair will be defined as uncomparable. A set of pairwise uncomparable integers will be called anti-hierarchic.
T(n,k) is the size of the maximal anti-hierarchic set of integers written with k digits in base n.
For example, for base n=2 and k=4 digits:
- 0 (0000) and 15 (1111) are comparable, while 6 (0110) and 9 (1001) are uncomparable,
- the maximal antihierarchic set is {3 (0011), 5 (0101), 6 (0110), 9 (1001), 10 (1010), 12 (1100)} with 6 elements that are all pairwise uncomparable. (End)

Examples

			Rows of square array start:
  1,    0,    0,    0,    0,    0,    0, ...
  1,    1,    1,    1,    1,    1,    1, ...
  1,    1,    2,    3,    6,   10,   20, ...
  1,    1,    3,    7,   19,   51,  141, ...
  1,    1,    4,   12,   44,  155,  580, ...
  1,    1,    5,   19,   85,  381, 1751, ...
  ...
Read by antidiagonals:
  1;
  0, 1;
  0, 1, 1;
  0, 1, 1, 1;
  0, 1, 2, 1, 1;
  0, 1, 3, 3, 1, 1;
  0, 1, 6, 7, 4, 1, 1;
  ...
		

Crossrefs

Programs

Formula

By the central limit theorem, T(n,k) is roughly n^(k-1)*sqrt(6/(Pi*k)).
T(n,k) = Sum{j=0,h/n} (-1)^j*binomial(k,j)*binomial(k-1+h-n*j,k-1) with h=floor(k*(n-1)/2), k>0. - Michel Marcus, Dec 01 2012

A063419 Central sextinomial coefficients.

Original entry on oeis.org

1, 6, 146, 4332, 135954, 4395456, 144840476, 4836766584, 163112472594, 5542414273884, 189456975899496, 6507792553644256, 224442843729333276, 7766945604528200460, 269557528994032024080, 9378595792117360310832
Offset: 0

Views

Author

Wolfdieter Lang, Jul 24 2001

Keywords

Comments

Largest coefficient of (Sum_{j=0..5} x^j)^(2*n). a(n)= A018901(2*n).
The exponent of 2 in prime factorization of a(n) for n=1,2,...,7 is in A023520. - Roman Witula, Oct 07 2012
Central coefficients in triangle A063260. - Zagros Lalo, Sep 25 2018

References

  • Shara Lalo and Zagros Lalo, Polynomial Expansion Theorems and Number Triangles, Zana Publishing, 2018, ISBN: 978-1-9995914-0-3, pp. 605, 606.
  • R. Witula and D. Slota, Central trinomial coefficients and convolution type identities, Congressus Numerantium 201 (2010), 109-126.

Crossrefs

Central q-nomial coefficients (appearing once) for q=2..5: A000984, A002426, A005721, A005191. For q=7: A025012.
Cf. A063260.

Programs

  • GAP
    Concatenation([1],List([1..15],n->Sum([0..Int(5*n/6)],k->(-1)^k*Binomial(2*n,k)*Binomial(7*n-6*k-1,2*n-1)))); # Muniru A Asiru, Sep 26 2018
  • Mathematica
    Table[Sum[(-1)^k*Binomial[2*n, k]*Binomial[7*n - 6*k - 1, 2*n - 1], {k, 0, Floor[5*n/6]}], {n, 0, 50}] (* G. C. Greubel, Mar 02 2017 *)
  • PARI
    concat([1], for(n=1,25, print1(sum(k=0,floor(5*n/6),(-1)^(k)*binomial(2*n,k)* binomial(7*n-6*k-1, 2*n-1)), ", "))) \\ G. C. Greubel, Mar 02 2017
    

Formula

a(n) = A063260(2*n, 5*n)= [x^(5*n)](Sum_{j=0..5} x^j)^(2*n).
a(n) = Sum_{k=0..floor(5*n/6)} ((-1)^(k)*binomial(2*n,k)*binomial(7*n-6*k-1, 2*n-1)). - Warut Roonguthai, May 22 2006
2*Pi*a(n) = Integral[-Pi;Pi] (sin(6*x)/sin(x))^(2*n) dx. The proof of this fact is in the Witula/Slota paper. - Roman Witula, Oct 07 2012
The Almkvist-Zeilberger algorithm in EKHAD establishes the following recurrence: -31104*(2*n+5)*(2*n+3)*(2*n+1)*(7*n+19)*(5*n+11)*(7*n+20)*(7*n+13)*(n+2)*(n+1)*a(n)+ 864*(7*n+20)*(2*n+5)*(2*n+3)*(n+2)*(25480*n^5+ 223496*n^4+755066*n^3+1223233*n^2+946889*n+279936)*a(n+1)- 6*(5*n+6)*(2*n+5)*(7*n+6)*(499359*n^6+ 6777015*n^5+38079431*n^4+113390385*n^3+18872398*n^2+ 166469280*n+60800544)*a(n+2)+ 5*(5*n+14)*(5*n+13)*(5*n+12)*(7*n+12)*(5*n+11)*(5*n+6)*(7*n+13)*(7*n+6)*(n+3)*a(n+3) = 0. - Doron Zeilberger, Apr 02 2013
a(n) ~ sqrt(3) * 36^n / sqrt(35*Pi*n). - Vaclav Kotesovec, Dec 09 2021

A061676 Triangle T(n,k) of number of ways of throwing k standard dice to produce a total of n.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 4, 6, 4, 1, 1, 5, 10, 10, 5, 1, 0, 6, 15, 20, 15, 6, 1, 0, 5, 21, 35, 35, 21, 7, 1, 0, 4, 25, 56, 70, 56, 28, 8, 1, 0, 3, 27, 80, 126, 126, 84, 36, 9, 1, 0, 2, 27, 104, 205, 252, 210, 120, 45, 10, 1, 0, 1, 25, 125, 305, 456, 462, 330, 165, 55, 11, 1
Offset: 1

Views

Author

Henry Bottomley, Apr 01 2002

Keywords

Examples

			Rows start:
1;
1,1;
1,2,1;
1,3,3,1;
1,4,6,4,1;
1,5,10,10,5,1;
0,6,15,20,15,6,1;
0,5,21,35,35,21,7,1;
etc.
T(8,2)=5 since 8 =2+6 =3+5 =4+4 =5+3 =6+2.
		

Crossrefs

First 21 terms as A007318 (see formula). Cf. A001592, A069713.
Cf. A030528 (2-sided dice), A078803 (3-sided), A213887 (4-sided), A213888 (5-sided).

Programs

  • Maple
    pts := 6; # A213889 and A061676
    g := 1/(1-t*z*add(z^i,i=0..pts-1)) ;
    for n from 1 to 13 do
        for k from 1 to n do
            coeftayl(g,z=0,n) ;
            coeftayl(%,t=0,k) ;
            printf("%d ",%) ;
        end do:
        printf("\n") ;
    end do: # R. J. Mathar, May 28 2025

Formula

T(n, k)=T(n-1, k-1)+T(n-2, k-1)+T(n-3, k-1)+T(n-4, k-1)+T(n-5, k-1)+T(n-6, k-1) starting with T(0, 0)=1. T(n, k)=T(7k-n, k); if n>6k or n6k-6, T(n, k)=C(7k-n-1, k-1); T([7k/2], k)=A018901(k).

A273975 Three-dimensional array written by antidiagonals in k,n: T(k,n,h) with k >= 1, n >= 0, 0 <= h <= n*(k-1) is the coefficient of x^h in the polynomial (1 + x + ... + x^(k-1))^n = ((x^k-1)/(x-1))^n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 2, 3, 2, 1, 1, 3, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 3, 4, 3, 2, 1, 1, 3, 6, 7, 6, 3, 1, 1, 4, 6, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 3, 4, 5, 4, 3, 2, 1, 1, 3, 6, 10, 12, 12, 10, 6, 3, 1, 1, 4, 10
Offset: 1

Views

Author

Andrey Zabolotskiy, Nov 10 2016

Keywords

Comments

Equivalently, T(k,n,h) is the number of ordered sets of n nonnegative integers < k with the sum equal to h.
From Juan Pablo Herrera P., Nov 21 2016: (Start)
T(k,n,h) is the number of possible ways of randomly selecting h cards from k-1 sets, each with n different playing cards. It is also the number of lattice paths from (0,0) to (n,h) using steps (1,0), (1,1), (1,2), ..., (1,k-1).
Shallow diagonal sums of each triangle with fixed k give the k-bonacci numbers. (End)
T(k,n,h) is the number of n-dimensional grid points of a k X k X ... X k grid, which are lying in the (n-1)-dimensional hyperplane which is at an L1 distance of h from one of the grid's corners, and normal to the corresponding main diagonal of the grid. - Eitan Y. Levine, Apr 23 2023

Examples

			For first few k and for first few n, the rows with h = 0..n*(k-1) are given:
k=1:  1;  1;  1;  1;  1; ...
k=2:  1;  1, 1;  1, 2, 1;  1, 3, 3, 1;  1, 4, 6, 4, 1; ...
k=3:  1;  1, 1, 1;  1, 2, 3, 2, 1;  1, 3, 6, 7, 6, 3, 1; ...
k=4:  1;  1, 1, 1, 1;  1, 2, 3, 4, 3, 2, 1; ...
For example, (1 + x + x^2)^3 = 1 + 3*x + 6*x^2 + 7*x^3 + 6*x^4 + 3*x^5 + x^6, hence T(3,3,2) = T(3,3,4) = 6.
From _Eitan Y. Levine_, Apr 23 2023: (Start)
Example for the repeated cumulative sum formula, for (k,n)=(3,3) (each line is the cumulative sum of the previous line, and the first line is the padded, alternating 3rd row from Pascal's triangle):
  1  0  0 -3  0  0  3  0  0 -1
  1  1  1 -2 -2 -2  1  1  1
  1  2  3  1 -1 -3 -2 -1
  1  3  6  7  6  3  1
which is T(3,3,h). (End)
		

Crossrefs

k-nomial arrays for fixed k=1..10: A000012, A007318, A027907, A008287, A035343, A063260, A063265, A171890, A213652, A213651.
Arrays for fixed n=0..6: A000012, A000012, A004737, A109439, A277949, A277950, A277951.
Central n-nomial coefficients for n=1..9, i.e., sequences with h=floor(n*(k-1)/2) and fixed n: A000012, A000984 (A001405), A002426, A005721 (A005190), A005191, A063419 (A018901), A025012, (A025013), A025014, A174061 (A025015), A201549, (A225779), A201550. Arrays: A201552, A077042, see also cfs. therein.
Triangle n=k-1: A181567. Triangle n=k: A163181.

Programs

  • Mathematica
    a = Table[CoefficientList[Sum[x^(h-1),{h,k}]^n,x],{k,10},{n,0,9}];
    Flatten@Table[a[[s-n,n+1]],{s,10},{n,0,s-1}]
    (* alternate program *)
    row[k_, n_] := Nest[Accumulate,Upsample[Table[((-1)^j)*Binomial[n,j],{j,0,n}],k],n][[;;n*(k-1)+1]] (* Eitan Y. Levine, Apr 23 2023 *)

Formula

T(k,n,h) = Sum_{i = 0..floor(h/k)} (-1)^i*binomial(n,i)*binomial(n+h-1-k*i,n-1). [Corrected by Eitan Y. Levine, Apr 23 2023]
From Eitan Y. Levine, Apr 23 2023: (Start)
(T(k,n,h))_{h=0..n*(k-1)} = f(f(...f(g(P))...)), where:
(x_i)_{i=0..m} denotes a tuple (in particular, the LHS contains the values for 0 <= h <= n*(k-1)),
f repeats n times,
f((x_i){i=0..m}) = (Sum{j=0..i} x_j)_{i=0..m} is the cumulative sum function,
g((x_i){i=0..m}) = (x(i/k) if k|i, otherwise 0)_{i=0..m*k} is adding k-1 zeros between adjacent elements,
and P=((-1)^i*binomial(n,i))_{i=0..n} is the n-th row of Pascal's triangle, with alternating signs. (End)
From Eitan Y. Levine, Jul 27 2023: (Start)
Recurrence relations, the first follows from the sequence's defining polynomial as mentioned in the Smarandache link:
T(k,n+1,h) = Sum_{i = 0..s-1} T(k,n,h-i)
T(k+1,n,h) = Sum_{i = 0..n} binomial(n,i)*T(k,n-i,h-i*k) (End)

A270918 Largest coefficient of (1+x+...+x^n)^(2*n).

Original entry on oeis.org

1, 2, 19, 580, 38165, 4395456, 786588243, 202384723528, 70886845397481, 32458256583753952, 18832730699014127291, 13507852690353224821652, 11738630472138500287398379, 12155701820213424461220851360, 14790850878997102285050287114419
Offset: 0

Views

Author

Vaclav Kotesovec, Mar 26 2016

Keywords

Crossrefs

Programs

  • Mathematica
    Table[Max[CoefficientList[Expand[Sum[x^k, {k, 0, n}]^(2n)], x]], {n, 0, 20}]
  • PARI
    a(n) = vecmax(Vec((sum(k=0,n,x^k))^(2*n))); \\ Michel Marcus, Apr 01 2016

Formula

a(n) ~ exp(2) * sqrt(3/Pi) * n^(2*n - 3/2).

Extensions

Typo in formula corrected by Vaclav Kotesovec, Dec 10 2021
Showing 1-7 of 7 results.