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.

A030662 Number of combinations of n things from 1 to n at a time, with repeats allowed.

Original entry on oeis.org

1, 5, 19, 69, 251, 923, 3431, 12869, 48619, 184755, 705431, 2704155, 10400599, 40116599, 155117519, 601080389, 2333606219, 9075135299, 35345263799, 137846528819, 538257874439, 2104098963719, 8233430727599, 32247603683099, 126410606437751, 495918532948103
Offset: 1

Views

Author

Donald Mintz (djmintz(AT)home.com)

Keywords

Comments

Add terms of an increasingly bigger diamond-shaped part of Pascal's triangle:
.......................... 1
............ 1 .......... 1 1
.. 1 ...... 1 1 ........ 1 2 1
. 1 1 =5 . 1 2 1 =19 .. 1 3 3 1 =69
.. 2 ...... 3 3 ........ 4 6 4
............ 6 ......... 10 10
.......................... 20
- Ralf Stephan, May 17 2004
The prime p divides a((p-1)/2) for p in A002144 (Pythagorean primes). - Alexander Adamchuk, Jul 04 2006
Also, number of square submatrices of a square matrix. - Jono Henshaw (jjono(AT)hotmail.com), Apr 22 2008
Partial sums of A051924. - J. M. Bergot, Jun 22 2013
Number of partitions with Ferrers diagrams that fit in an n X n box (excluding the empty partition of 0). - Michael Somos, Jun 02 2014
Also number of non-descending sequences with length and last number are less or equal to n, and also the number of integer partitions (of any positive integer) with length and largest part are less or equal to n. - Zlatko Damijanic, Dec 06 2024

Examples

			G.f. = x + 5*x^2 + 19*x^3 + 69*x^4 + 251*x^5 + 923*x^6 + 3431*x^7 + ...
		

Crossrefs

Column k=2 of A047909.
Central column of triangle A014473.
Right-hand column 2 of triangle A102541.

Programs

  • Magma
    [(n+1)*Catalan(n)-1: n in [1..40]]; // G. C. Greubel, Apr 07 2024
  • Maple
    seq(sum((binomial(n,m))^2,m=1..n),n=1..23); # Zerinvary Lajos, Jun 19 2008
    f:=n->add( add( binomial(i+j,i), i=0..n),j=0..n); [seq(f(n),n=0..12)]; # N. J. A. Sloane, Jan 31 2009
  • Mathematica
    Table[Sum[Sum[(2n-i-j)!/(n-i)!/(n-j)!,{i,1,n}],{j,1,n}],{n,1,20}] (* Alexander Adamchuk, Jul 04 2006 *)
    a[n_] := 2*(2*n-1)!/(n*(n-1)!^2)-1; Table[a[n], {n, 1, 26}] (* Jean-François Alcover, Oct 11 2012, from first formula *)
  • PARI
    a(n)=binomial(2*n,n)-1 \\ Charles R Greathouse IV, Jun 26 2013
    
  • Python
    from math import comb
    def a(n): return comb(2*n, n) - 1
    print([a(n) for n in range(1, 27)]) # Michael S. Branicky, Jul 11 2023
    
  • Sage
    def a(n) : return binomial(2*n,n) - 1
    [a(n) for n in (1..26)] # Peter Luschny, Apr 21 2012
    

Formula

a(n) = A000984(n) - 1.
a(n) = 2*A001700(n-1) - 1.
a(n) = 2*(2*n-1)!/(n!*(n-1)!)-1.
a(n) = Sum_{k=1..n} binomial(n, k)^2. - Benoit Cloitre, Aug 20 2002
a(n) = Sum_{j=0..n} Sum_{i=j..n+j} binomial(i, j). - Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de), Jul 23 2003
a(n) = Sum_{i=0..n-1} Sum_{j=0..n-1} binomial(i+j, i). - N. J. A. Sloane, Jan 31 2009
Also for n>1: a(n)=(2*n)!/(n!)^2-1. - Hugo Pfoertner, Feb 10 2004
a(n) = Sum_{j=1..n} Sum_{i=1..n} (2n-i-j)!/((n-i)!*(n-j)!). - Alexander Adamchuk, Jul 04 2006
a(n) = A115112(n) + 1. - Jono Henshaw (jjono(AT)hotmail.com), Apr 22 2008
G.f.: Q(0)*(1-4*x)/x - 1/x/(1-x), where Q(k)= 1 + 4*(2*k+1)*x/( 1 - 1/(1 + 2*(k+1)/Q(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 11 2013
D-finite with recurrence: n*a(n) +2*(-3*n+2)*a(n-1) +(9*n-14)*a(n-2) +2*(-2*n+5)*a(n-3)=0. - R. J. Mathar, Jun 25 2013
0 = a(n)*(+16*a(n+1) - 70*a(n+2) + 68*a(n+3) - 14*a(n+4)) + a(n+1)*(-2*a(n+1) + 61*a(n+2) - 96*a(n+3) + 23*a(n+4)) + a(n+2)*(-6*a(n+2) + 31*a(n+3) - 10*a(n+4)) + a(n+3)*(-2*a(n+3) + a(n+4)) for all n in Z. - Michael Somos, Jun 02 2014
From Ilya Gutkovskiy, Jan 25 2017: (Start)
O.g.f.: (1 - x - sqrt(1 - 4*x))/((1 - x)*sqrt(1 - 4*x)).
E.g.f.: exp(x)*(exp(x)*BesselI(0,2*x) - 1). (End)
a(n) = 3*n*Sum_{k=1..n} (-1)^(k+1)/(2*n+k)*binomial(2*n+k,n-k). - Vladimir Kruchinin, Jul 29 2025
a(n) = n * binomial(2*n, n) * Sum_{k = 1..n} 1/(k*binomial(n+k, k)). - Peter Bala, Aug 05 2025

A014473 Pascal's triangle - 1: Triangle read by rows: T(n, k) = A007318(n, k) - 1.

Original entry on oeis.org

0, 0, 0, 0, 1, 0, 0, 2, 2, 0, 0, 3, 5, 3, 0, 0, 4, 9, 9, 4, 0, 0, 5, 14, 19, 14, 5, 0, 0, 6, 20, 34, 34, 20, 6, 0, 0, 7, 27, 55, 69, 55, 27, 7, 0, 0, 8, 35, 83, 125, 125, 83, 35, 8, 0, 0, 9, 44, 119, 209, 251, 209, 119, 44, 9, 0, 0, 10, 54, 164, 329, 461, 461, 329, 164, 54, 10, 0
Offset: 0

Views

Author

Keywords

Comments

Indexed as a square array A(n,k): If X is an (n+k)-set and Y a fixed k-subset of X then A(n,k) is equal to the number of n-subsets of X intersecting Y. - Peter Luschny, Apr 20 2012

Examples

			Triangle begins:
   0;
   0, 0;
   0, 1,  0;
   0, 2,  2,  0;
   0, 3,  5,  3,  0;
   0, 4,  9,  9,  4,  0;
   0, 5, 14, 19, 14,  5, 0;
   0, 6, 20, 34, 34, 20, 6, 0;
   ...
Seen as a square array read by antidiagonals:
  [0] 0, 0,  0,  0,   0,   0,   0,    0,    0,    0,    0,     0, ... A000004
  [1] 0, 1,  2,  3,   4,   5,   6,    7,    8,    9,   10,    11, ... A001477
  [2] 0, 2,  5,  9,  14,  20,  27,   35,   44,   54,   65,    77, ... A000096
  [3] 0, 3,  9, 19,  34,  55,  83,  119,  164,  219,  285,   363, ... A062748
  [4] 0, 4, 14, 34,  69, 125, 209,  329,  494,  714, 1000,  1364, ... A063258
  [5] 0, 5, 20, 55, 125, 251, 461,  791, 1286, 2001, 3002,  4367, ... A062988
  [6] 0, 6, 27, 83, 209, 461, 923, 1715, 3002, 5004, 8007, 12375, ... A124089
		

Crossrefs

Triangle without zeros: A014430.
Related: A323211 (A007318(n, k) + 1).
A000295 (row sums), A059841 (alternating row sums), A030662(n-1) (central terms).
Columns include A000096, A062748, A062988, A063258.
Diagonals of A(n, n+d): A030662 (d=0), A010763 (d=1), A322938 (d=2).

Programs

  • Haskell
    a014473 n k = a014473_tabl !! n !! k
    a014473_row n = a014473_tabl !! n
    a014473_tabl = map (map (subtract 1)) a007318_tabl
    -- Reinhard Zumkeller, Apr 10 2012
    
  • Magma
    [Binomial(n,k)-1: k in [0..n], n in [0..12]]; // G. C. Greubel, Apr 08 2024
    
  • Maple
    with(combstruct): for n from 0 to 11 do seq(-1+count(Combination(n), size=m), m = 0 .. n) od; # Zerinvary Lajos, Apr 09 2008
    # The rows of the square array:
    Arow := proc(n, len) local gf, ser;
    gf := (x - 1)^(-n - 1) + (-1)^(n + 1)/(x*(x - 1));
    ser := series(gf, x, len+2): seq((-1)^(n+1)*coeff(ser, x, j), j=0..len) end:
    for n from 0 to 9 do lprint([n], Arow(n, 12)) od; # Peter Luschny, Feb 13 2019
  • Mathematica
    Table[Binomial[n,k] -1, {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Apr 08 2024 *)
  • SageMath
    flatten([[binomial(n,k)-1 for k in range(n+1)] for n in range(13)]) # G. C. Greubel, Apr 08 2024

Formula

G.f.: x^2*y/((1 - x)*(1 - x*y)*(1 - x*(1 + y))). - Ralf Stephan, Jan 24 2005
T(n,k) = A109128(n,k) - A007318(n,k), 0 <= k <= n. - Reinhard Zumkeller, Apr 10 2012
T(n, k) = T(n-1, k-1) + T(n-1, k) + 1, 0 < k < n with T(n, 0) = T(n, n) = 0. - Reinhard Zumkeller, Jul 18 2015
If seen as a square array read by antidiagonals the generating function of row n is: G(n) = (x - 1)^(-n - 1) + (-1)^(n + 1)/(x*(x - 1)). - Peter Luschny, Feb 13 2019
From G. C. Greubel, Apr 08 2024: (Start)
T(n, n-k) = T(n, k).
Sum_{k=0..floor(n/2)} T(n-k, k) = A129696(n-2).
Sum_{k=0..floor(n/2)} (-1)^k*T(n-k, k) = b(n-1), where b(n) is the repeating pattern {0, 0, -1, -2, -1, 1, 1, -1, -2, -1, 0, 0}_{n=0..11}, with b(n) = b(n-12). (End)

Extensions

More terms from Erich Friedman

A124051 Quasi-mirror of A062196 formatted as a triangular array.

Original entry on oeis.org

3, 6, 8, 10, 30, 15, 15, 80, 90, 24, 21, 175, 350, 210, 35, 28, 336, 1050, 1120, 420, 48, 36, 588, 2646, 4410, 2940, 756, 63, 45, 960, 5880, 14112, 14700, 6720, 1260, 80, 55, 1485, 11880, 38808, 58212, 41580, 13860, 1980, 99, 66, 2200, 22275, 95040, 194040, 199584, 103950, 26400, 2970, 120
Offset: 0

Views

Author

Zerinvary Lajos, Nov 03 2006

Keywords

Examples

			Triangle begins as:
   3;
   6,    8;
  10,   30,    15;
  15,   80,    90,    24;
  21,  175,   350,   210,     35;
  28,  336,  1050,  1120,    420,     48;
  36,  588,  2646,  4410,   2940,    756,     63;
  45,  960,  5880, 14112,  14700,   6720,   1260,    80;
  55, 1485, 11880, 38808,  58212,  41580,  13860,  1980,   99;
  66, 2200, 22275, 95040, 194040, 199584, 103950, 26400, 2970, 120;
		

Crossrefs

Columns k: A000217(n+2) (k=0), A002417(n+1) (k=1), A001297(n) (k=2), A105946(n-2) (k=3), A105947(n-3) (k=4), A105948(n-4) (k=5), A107319(n-5) (k=6).
Diagonals: A005563(n+1) (k=n), A033487(n) (k=n-1), A027790(n) (k=n-2), A107395(n-3) (k=n-3), A107396(n-4) (k=n-4), A107397(n-5) (k=n-5), A107398(n-6) (k=n-6), A107399(n-7) (k=n-7).
Sums: A322938(n+1) (row).

Programs

  • Magma
    A124051:= func< n,k | Binomial(n+1,n-k+1)*Binomial(n+3,n-k+1) >;
    [A124051(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Feb 07 2025
    
  • Maple
    for n from 0 to 10 do seq(binomial(n,i-1)*binomial(n+2,n+1-i), i=1..n ) od;
  • Mathematica
    A124051[n_, k_]:= Binomial[n+1,n-k+1]*Binomial[n+3,n-k+1];
    Table[A124051[n,k], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Feb 07 2025 *)
  • SageMath
    def A124051(n,k): return binomial(n+1,n-k+1)*binomial(n+3,n-k+1)
    print(flatten([[A124051(n,k) for k in range(n+1)] for n in range(13)])) # G. C. Greubel, Feb 07 2025

Formula

From G. C. Greubel, Feb 07 2025: (Start)
T(n, k) = binomial(n+1, n-k+1)*binomial(n+3, n-k+1).
T(2*n, n) = (1/2)*A000894(n) + (5/2)*[n=0].
Sum_{k=0..n} (-1)^k*T(n, k) = (1/2)*( (1+(-1)^n)*(-1)^(n/2)*A286033((n+4)/2) + (1-(-1)^n)*((-1)^((n+1)/2)*A000108((n+1)/2) - 1) ). (End)
Showing 1-3 of 3 results.