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

A034266 Partial sums of A027818.

Original entry on oeis.org

0, 1, 15, 99, 435, 1485, 4257, 10725, 24453, 51480, 101530, 189618, 338130, 579462, 959310, 1540710, 2408934, 3677355, 5494401, 8051725, 11593725, 16428555, 22940775, 31605795, 43006275, 57850650, 76993956, 101461140, 132473044, 171475260, 220170060, 280551612
Offset: 0

Views

Author

Keywords

References

  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, N. Y., 1964, pp. 194-196.

Crossrefs

a(n)=f(n, 6) where f is given in A034261.
Cf. A093564 ((7, 1) Pascal, column m=8).
Cf. similar sequences listed in A254142.

Programs

  • GAP
    List([0..35], n-> (7*n+1)*Binomial(n+6,7)/8); # G. C. Greubel, Aug 29 2019
  • Magma
    [0] cat [(7*n+8)*Binomial(n+7, 7)/8: n in [0..30]]; // Vincenzo Librandi, Mar 20 2015
    
  • Maple
    f:=n->(7*n+8)*binomial(n+7, 7)/8; [seq(f(n),n=-1..40)]; # N. J. A. Sloane, Mar 25 2015
  • Mathematica
    CoefficientList[Series[x(1+6x)/(1-x)^9, {x, 0, 30}], x] (* Vincenzo Librandi, Mar 20 2015 *)
    Table[(7*n+1)*Binomial[n+6,7]/8, {n,0,35}] (* G. C. Greubel, Aug 29 2019 *)
  • PARI
    lista(nn) = for (n=0, nn, print1((7*n+1)*binomial(n+6,7)/8, ", ")); \\ Michel Marcus, Mar 20 2015
    
  • Sage
    [(7*n+1)*binomial(n+6,7)/8 for n in (0..35)] # G. C. Greubel, Aug 29 2019
    

Formula

a(n) = (7*n+1)*binomial(n+6, 7)/8.
G.f.: x*(1+6*x)/(1-x)^9.
E.g.f.: x*(8! +262080*x +383040*x^2 +210000*x^3 +52080*x^4 +6216*x^5 + 344*x^6 +7*x^7)*exp(x)/8!

Extensions

Better description from Barry E. Williams, Jan 25 2000
Corrected and extended by N. J. A. Sloane, Apr 21 2000
More terms from Michel Marcus, Mar 20 2015

A062190 Coefficient triangle of certain polynomials N(5; m,x).

Original entry on oeis.org

1, 1, 6, 1, 14, 21, 1, 24, 84, 56, 1, 36, 216, 336, 126, 1, 50, 450, 1200, 1050, 252, 1, 66, 825, 3300, 4950, 2772, 462, 1, 84, 1386, 7700, 17325, 16632, 6468, 792, 1, 104, 2184, 16016, 50050, 72072, 48048, 13728, 1287, 1, 126, 3276, 30576, 126126, 252252, 252252, 123552, 27027, 2002
Offset: 0

Views

Author

Wolfdieter Lang, Jun 19 2001

Keywords

Comments

The e.g.f. of the m-th (unsigned) column sequence without leading zeros of the generalized (a=5) Laguerre triangle L(5; n+m,m)= A062138(n+m,m), n >= 0, is N(5; m,x)/(1-x)^(2*(m+3)), with the row polynomials N(5; m,x) := Sum_{k=0..m} a(m,k)*x^k.

Examples

			Triangle begins as:
  1;
  1,   6;
  1,  14,   21;
  1,  24,   84,    56;
  1,  36,  216,   336,    126;
  1,  50,  450,  1200,   1050,    252;
  1,  66,  825,  3300,   4950,   2772,     462;
  1,  84, 1386,  7700,  17325,  16632,    6468,    792;
  1, 104, 2184, 16016,  50050,  72072,   48048,  13728,   1287;
  1, 126, 3276, 30576, 126126, 252252,  252252, 123552,  27027,  2002;
  1, 150, 4725, 54600, 286650, 756756, 1051050, 772200, 289575, 50050, 3003;
		

Crossrefs

Family of polynomials (see A062145): A008459 (c=1), A132813 (c=2), A062196 (c=3), A062145 (c=4), A062264 (c=5), this sequence (c=6).
Columns k: A028557 (k=1), A104676 (k=2), A104677 (k=3), A104678 (k=4), A104679 (k=5), A104680 (k=6).
Diagonals: A000389 (k=n), A027818 (k=n-1), A104670 (k=n-2), A104671 (k=n-3), A104672 (k=n-4), A104673 (k=n-5), A104674 (k=n-6).
Cf. A003516 (row sums), A113894 (main diagonal).

Programs

  • Magma
    A062190:= func< n,k | Binomial(n,k)*Binomial(n+5,k) >;
    [A062190(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Feb 28 2025
    
  • Maple
    A062190 := proc(m,k)
        add( (binomial(m, j)*(2*m+5-j)!/((m+5)!*(m-j)!))*(x^(m-j))*(1-x)^j,j=0..m) ;
        coeftayl(%,x=0,k) ;
    end proc: # R. J. Mathar, Nov 29 2015
  • Mathematica
    NN[5, m_, x_] := x^m*(2*m+5)!*Hypergeometric2F1[-m, -m, -2*m-5, (x-1)/x]/((m+5)!*m!); Table[CoefficientList[NN[5, m, x], x], {m, 0, 8}] // Flatten (* Jean-François Alcover, Sep 18 2013 *)
    A062190[n_,k_]:= Binomial[n,k]*Binomial[n+5,k];
    Table[A062190[n,k], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Feb 28 2025 *)
  • SageMath
    def A062190(n,k): return binomial(n,k)*binomial(n+5,k)
    print(flatten([[A062190(n,k) for k in range(n+1)] for n in range(13)])) # G. C. Greubel, Feb 28 2025

Formula

T(m, k) = [x^k]N(5; m, x), with N(5; m, x) = ((1-x)^(2*(m+3)))*(d^m/dx^m)(x^m/(m!*(1-x)^(m+6))).
N(5; m, x) = Sum_{j=0..m} ((binomial(m, j)*(2*m+5-j)!/((m+5)!*(m-j)!))*(x^(m-j))*(1-x)^j).
N(5; m, x)= x^m*(2*m+5)! * 2F1(-m, -m; -2*m-5; (x-1)/x)/((m+5)!*m!). - Jean-François Alcover, Sep 18 2013
T(n, k) = binomial(n, k)*binomial(n+5, k). - G. C. Greubel, Feb 28 2025

A245334 A factorial-like triangle read by rows: T(0,0) = 1; T(n+1,0) = T(n,0)+1; T(n+1,k+1) = T(n,0)*T(n,k), k=0..n.

Original entry on oeis.org

1, 2, 1, 3, 4, 2, 4, 9, 12, 6, 5, 16, 36, 48, 24, 6, 25, 80, 180, 240, 120, 7, 36, 150, 480, 1080, 1440, 720, 8, 49, 252, 1050, 3360, 7560, 10080, 5040, 9, 64, 392, 2016, 8400, 26880, 60480, 80640, 40320, 10, 81, 576, 3528, 18144, 75600, 241920, 544320
Offset: 0

Views

Author

Reinhard Zumkeller, Aug 30 2014

Keywords

Comments

row(0) = {1}; row(n+1) = row(n) multiplied by n and prepended with (n+1);
A111063(n+1) = sum of n-th row;
T(2*n,n) = A002690(n), central terms;
T(n,0) = n + 1;
T(n,1) = A000290(n), n > 0;
T(n,2) = A011379(n-1), n > 1;
T(n,3) = A047927(n), n > 2;
T(n,4) = A192849(n-1), n > 3;
T(n,5) = A000142(5) * A027810(n-5), n > 4;
T(n,6) = A000142(6) * A027818(n-6), n > 5;
T(n,7) = A000142(7) * A056001(n-7), n > 6;
T(n,8) = A000142(8) * A056003(n-8), n > 7;
T(n,9) = A000142(9) * A056114(n-9), n > 8;
T(n,n-10) = 11 * A051431(n-10), n > 9;
T(n,n-9) = 10 * A049398(n-9), n > 8;
T(n,n-8) = 9 * A049389(n-8), n > 7;
T(n,n-7) = 8 * A049388(n-7), n > 6;
T(n,n-6) = 7 * A001730(n), n > 5;
T(n,n-5) = 6 * A001725(n), n > 5;
T(n,n-4) = 5 * A001720(n), n > 4;
T(n,n-3) = 4 * A001715(n), n > 2;
T(n,n-2) = A070960(n), n > 1;
T(n,n-1) = A052849(n), n > 0;
T(n,n) = A000142(n);
T(n,k) = A137948(n,k) * A007318(n,k), 0 <= k <= n.

Examples

			.  0:   1;
.  1:   2,  1;
.  2:   3,  4,   2;
.  3:   4,  9,  12,    6;
.  4:   5, 16,  36,   48,    24;
.  5:   6, 25,  80,  180,   240,   120;
.  6:   7, 36, 150,  480,  1080,  1440,    720;
.  7:   8, 49, 252, 1050,  3360,  7560,  10080,   5040;
.  8:   9, 64, 392, 2016,  8400, 26880,  60480,  80640,  40320;
.  9:  10, 81, 576, 3528, 18144, 75600, 241920, 544320, 725760, 362880.
		

Crossrefs

Programs

  • Haskell
    a245334 n k = a245334_tabl !! n !! k
    a245334_row n = a245334_tabl !! n
    a245334_tabl = iterate (\row@(h:_) -> (h + 1) : map (* h) row) [1]
  • Mathematica
    Table[(n!)/((n - k)!)*(n + 1 - k), {n, 0, 9}, {k, 0, n}] // Flatten (* Michael De Vlieger, Sep 10 2017 *)

Formula

T(n,k) = n!*(n+1-k)/(n-k)!. - Werner Schulte, Sep 09 2017

A093564 (7,1) Pascal triangle.

Original entry on oeis.org

1, 7, 1, 7, 8, 1, 7, 15, 9, 1, 7, 22, 24, 10, 1, 7, 29, 46, 34, 11, 1, 7, 36, 75, 80, 45, 12, 1, 7, 43, 111, 155, 125, 57, 13, 1, 7, 50, 154, 266, 280, 182, 70, 14, 1, 7, 57, 204, 420, 546, 462, 252, 84, 15, 1, 7, 64, 261, 624, 966, 1008, 714, 336, 99, 16, 1, 7, 71, 325, 885
Offset: 0

Views

Author

Wolfdieter Lang, Apr 22 2004

Keywords

Comments

The array F(7;n,m) gives in the columns m>=1 the figurate numbers based on A016993, including the 9-gonal numbers A001106, (see the W. Lang link).
This is the seventh member, d=7, in the family of triangles of figurate numbers, called (d,1) Pascal triangles: A007318 (Pascal), A029653, A093560-3, for d=1..6.
This is an example of a Riordan triangle (see A093560 for a comment and A053121 for a comment and the 1991 Shapiro et al. reference on the Riordan group). Therefore the o.g.f. for the row polynomials p(n,x):=Sum_{m=0..n} a(n,m)*x^m is G(z,x)=(1+6*z)/(1-(1+x)*z).
The SW-NE diagonals give A022097(n-1) = Sum_{k=0..ceiling((n-1)/2)} a(n-1-k,k), n >= 1, with n=0 value 6. Observation by Paul Barry, Apr 29 2004. Proof via recursion relations and comparison of inputs.

Examples

			Triangle begins
  [1];
  [7,  1];
  [7,  8,  1];
  [7, 15,  9,  1];
  ...
		

References

  • Kurt Hawlitschek, Johann Faulhaber 1580-1635, Veroeffentlichung der Stadtbibliothek Ulm, Band 18, Ulm, Germany, 1995, Ch. 2.1.4. Figurierte Zahlen.
  • Ivo Schneider: Johannes Faulhaber 1580-1635, Birkhäuser, Basel, Boston, Berlin, 1993, ch. 5, pp. 109-122.

Crossrefs

Row sums: A000079(n+2), n>=1, 1 for n=0, alternating row sums are 1 for n=0, 6 for n=2 and 0 otherwise.
The column sequences give for m=1..9: A016993, A001106 (9-gonal), A007584, A051740, A051877, A050403, A027818, A034266, A055994.
Cf. A093565 (d=8).

Programs

  • Haskell
    a093564 n k = a093564_tabl !! n !! k
    a093564_row n = a093564_tabl !! n
    a093564_tabl = [1] : iterate
                   (\row -> zipWith (+) ([0] ++ row) (row ++ [0])) [7, 1]
    -- Reinhard Zumkeller, Sep 01 2014
  • Maple
    N:= 20: # to get the first N rows
    T:=Matrix(N,N):
    T[1,1]:= 1:
    for m from 2 to N do
    T[m,1]:= 7:
    T[m,2..m]:= T[m-1,1..m-1] + T[m-1,2..m];
    od:
    for m from 1 to N do
    convert(T[m,1..m],list)
    od; # Robert Israel, Dec 28 2014

Formula

a(n, m)=F(7;n-m, m) for 0<= m <= n, otherwise 0, with F(7;0, 0)=1, F(7;n, 0)=7 if n>=1 and F(7;n, m):=(7*n+m)*binomial(n+m-1, m-1)/m if m>=1.
Recursion: a(n, m)=0 if m>n, a(0, 0)= 1; a(n, 0)=7 if n>=1; a(n, m)= a(n-1, m) + a(n-1, m-1).
G.f. column m (without leading zeros): (1+6*x)/(1-x)^(m+1), m>=0.
T(n, k) = C(n, k) + 6*C(n-1, k). - Philippe Deléham, Aug 28 2005
exp(x) * e.g.f. for row n = e.g.f. for diagonal n. For example, for n = 3 we have exp(x)*(7 + 15*x + 9*x^2/2! + x^3/3!) = 7 + 22*x + 46*x^2/2! + 80*x^3/3! + 125*x^4/4! + .... The same property holds more generally for Riordan arrays of the form ( f(x), x/(1 - x) ). - Peter Bala, Dec 22 2014

A128629 A triangular array generated by moving Pascal sequences to prime positions and embedding new sequences at the nonprime locations. (cf. A007318 and A000040).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 4, 6, 4, 1, 1, 4, 9, 10, 5, 1, 1, 6, 10, 16, 15, 6, 1, 1, 5, 18, 20, 25, 21, 7, 1, 1, 8, 15, 40, 35, 36, 28, 8, 1, 1, 9, 27, 35, 75, 56, 49, 36, 9, 1
Offset: 1

Views

Author

Alford Arnold, Mar 29 2007

Keywords

Comments

The array can be constructed by beginning with A007318 (Pascal's triangle) placing each diagonal on a prime row. The other rows are filled in by mapping the prime factorization of the row number to the known sequences on the prime rows and multiplying term by term.

Examples

			Row six begins 1 6 18 40 75 126 ... because rows two and three are
1 2 3 4 5 6 ...
1 3 6 10 15 21 ...
The array begins
1 1 1 1 1 1 1 1 1 A000012
1 2 3 4 5 6 7 8 9 A000027
1 3 6 10 15 21 28 36 45 A000217
1 4 9 16 25 36 49 64 81 A000290
1 4 10 20 35 56 84 120 165 A000292
1 6 18 40 75 126 196 288 405 A002411
1 5 15 35 70 126 210 330 495 A000332
1 8 27 64 125 216 343 512 729 A000578
1 9 36 100 225 441 784 1296 2025 A000537
1 8 30 80 175 336 588 960 1485 A002417
1 6 21 56 126 252 462 792 1287 A000389
1 12 54 160 375 756 1372 2304 3645 A019582
1 7 28 84 210 462 924 1716 3003 A000579
1 10 45 140 350 756 1470 2640 4455 A027800
1 12 60 200 525 1176 2352 4320 7425 A004302
1 16 81 256 625 1296 2401 4096 6561 A000583
1 8 36 120 330 792 1716 3432 6435 A000580
1 18 108 400 1125 2646 5488 10368 18225 A019584
1 9 45 165 495 1287 3003 6435 12870 A000581
1 16 90 320 875 2016 4116 7680 13365 A119771
1 15 90 350 1050 2646 5880 11880 22275 A001297
1 12 63 224 630 1512 3234 6336 11583 A027810
1 10 55 220 715 2002 5005 11440 24310 A000582
1 24 162 640 1875 4536 9604 18432 32805 A019583
1 16 100 400 1225 3136 7056 14400 27225 A001249
1 14 84 336 1050 2772 6468 13728 27027 A027818
1 27 216 1000 3375 9261 21952 46656 91125 A059827
1 20 135 560 1750 4536 10290 21120 40095 A085284
		

Crossrefs

Cf. A064553 (second diagonal), A080688 (second diagonal resorted).

Programs

  • Maple
    A128629 := proc(n,m) if n = 1 then 1; elif isprime(n) then p := numtheory[pi](n) ; binomial(p+m-1,p) ; else a := 1 ; for p in ifactors(n)[2] do a := a* procname(op(1,p),m)^ op(2,p) ; od: fi; end: # R. J. Mathar, Sep 09 2009

Extensions

A-number added to each row of the examples by R. J. Mathar, Sep 09 2009

A108681 a(n) = (n+1)*(n+2)^2*(n+3)*(n+4)*(n+5)*(2*n+3)/720.

Original entry on oeis.org

1, 15, 98, 420, 1386, 3822, 9240, 20196, 40755, 77077, 138138, 236600, 389844, 621180, 961248, 1449624, 2136645, 3085467, 4374370, 6099324, 8376830, 11347050, 15177240, 20065500, 26244855, 33987681, 43610490, 55479088, 70014120, 87697016, 109076352, 134774640
Offset: 0

Views

Author

Emeric Deutsch, Jun 18 2005

Keywords

Comments

Kekulé numbers for certain benzenoids.

References

  • S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (p. 232, # 4).

Programs

  • Maple
    G:=factor(sum(a(n)*z^n,n=0..infinity)); series(G,z=0,37);
  • Mathematica
    Table[(n+1)(n+2)^2(n+3)(n+4)(n+5)(2n+3)/720,{n,0,30}] (* or *) LinearRecurrence[{8,-28,56,-70,56,-28,8,-1},{1,15,98,420,1386,3822,9240,20196},30] (* Harvey P. Dale, Sep 23 2017 *)

Formula

G.f.: (1+x)*(1+6*x)/(1-x)^8.
From Amiram Eldar, Jun 02 2022: (Start)
Sum_{n>=0} 1/a(n) = 20*Pi^2 - 3072*log(2)/7 + 4531/42.
Sum_{n>=0} (-1)^n/a(n) = 768*Pi/7 - 10*Pi^2 - 256*log(2)/7 - 9227/42. (End)
a(n) = A027818(n)+A027818(n-1). - R. J. Mathar, Jul 22 2022
Showing 1-6 of 6 results.