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

A016169 a(n) = 7^n - 6^n.

Original entry on oeis.org

0, 1, 13, 127, 1105, 9031, 70993, 543607, 4085185, 30275911, 222009073, 1614529687, 11664504865, 83828316391, 599858908753, 4277376525367, 30411820662145, 215703854542471, 1526853641242033, 10789535445362647
Offset: 0

Views

Author

Keywords

Comments

a(n) is also the number of n-digit numbers whose smallest decimal digit is 3. - Stefano Spezia, Nov 15 2023

Crossrefs

Programs

  • Magma
    [n le 2 select n-1 else 13*Self(n-1) -42*Self(n-2): n in [1..31]]; // G. C. Greubel, Nov 10 2024
    
  • Maple
    a:=n->sum(6^(n-j)*binomial(n,j),j=1..n): seq(a(n), n=0..30); # Zerinvary Lajos, Apr 18 2009
  • Mathematica
    Table[7^n-6^n,{n,0,30}] (* or *) LinearRecurrence[{13,-42},{0,1},30] (* Harvey P. Dale, Apr 25 2020 *)
  • SageMath
    A016169=BinaryRecurrenceSequence(13,-42,0,1)
    [A016169(n) for n in range(41)] # G. C. Greubel, Nov 10 2024

Formula

G.f.: x/((1-6*x)*(1-7*x)).
E.g.f.: exp(7*x) - exp(6*x). - Mohammad K. Azarian, Jan 14 2009
a(0)=0, a(n) = 7*a(n-1) + 6^(n-1). - Vincenzo Librandi, Feb 09 2011
a(0)=0, a(1)=1, a(n) = 13*a(n-1) - 42*a(n-2). - Vincenzo Librandi, Feb 09 2011

A081200 6th binomial transform of (0,1,0,1,0,1,...), A000035.

Original entry on oeis.org

0, 1, 12, 109, 888, 6841, 51012, 372709, 2687088, 19200241, 136354812, 964249309, 6798573288, 47834153641, 336059778612, 2358521965909, 16540171339488, 115933787267041, 812299450322412, 5689910849522509, 39848449432985688, 279034513462540441, 1953718431395986212
Offset: 0

Views

Author

Paul Barry, Mar 11 2003

Keywords

Comments

Binomial transform of A081199.
Conjecture (verified up to a(9)): Number of collinear 5-tuples of points in a 5 X 5 X 5 X ... n-dimensional cubic grid. - Ron Hardin, May 24 2010
a(n) is also the total number of words of length n, over an alphabet of seven letters, of which one of them appears an odd number of times. See the Lekraj Beedassy, Jul 22 2003, comment on A006516 (4-letter case), and the Balakrishnan reference there. For the 2-, 3-, 5-, 6- and 8-letter case analogs see A131577, A003462, A005059, A081199, A081201 respectively. - Wolfdieter Lang, Jul 17 2017

Examples

			The a(2) = 12 words of length 2 over {A, B, C, D, E, F, G} with say, A, appearing an odd number of times (that is once) are: AB, AC, AD, AE, AF, AG; BA, CA, DA, EA, FA, GA. - _Wolfdieter Lang_, Jul 17 2017
		

Crossrefs

Cf. A000035, A003462, A005059, A006516, A081199, A081201 (binomial transform, and 8-letter analog), A121213, A131577.
Apart from offset same as A016161.

Programs

  • Magma
    [7^n/2-5^n/2: n in [0..25]]; // Vincenzo Librandi, Aug 07 2013
  • Mathematica
    CoefficientList[Series[x / ((1 - 5 x) (1 - 7 x)), {x, 0, 20}], x] (* Vincenzo Librandi, Aug 07 2013 *)
    LinearRecurrence[{12,-35},{0,1},30] (* Harvey P. Dale, Feb 07 2014 *)
  • Sage
    [lucas_number1(n,12,35) for n in range(0, 21)] # Zerinvary Lajos, Apr 27 2009
    

Formula

a(n) = 12*a(n-1) - 35*a(n-2), a(0) = 0, a(1) = 1.
G.f.: x/((1-5*x)*(1-7*x)).
a(n) = 7^n/2 - 5^n/2.
a(n) = Sum_{k=0..n-1} 7^k * 5^(n-k-1), with a(0)=0. - Reinhard Zumkeller, Aug 01 2010
a(n) = A121213(n)/2. - Reinhard Zumkeller, Aug 01 2010
E.g.f.: exp(5*x)*(exp(2*x) - 1)/2. - Stefano Spezia, Jun 19 2021

A016170 Expansion of 1/((1-6*x)*(1-8*x)).

Original entry on oeis.org

1, 14, 148, 1400, 12496, 107744, 908608, 7548800, 62070016, 506637824, 4113568768, 33271347200, 268347559936, 2159841173504, 17357093552128, 139326933401600, 1117436577120256, 8956419276406784, 71752914167922688
Offset: 0

Views

Author

Keywords

Crossrefs

Programs

  • Magma
    [n le 2 select 14^(n-1) else 14*Self(n-1) -48*Self(n-2): n in [1..31]]; // G. C. Greubel, Nov 10 2024
    
  • Maple
    A016170:=n->4*8^n-3*6^n: seq(A016170(n), n=0..30); # Wesley Ivan Hurt, May 03 2017
  • Mathematica
    CoefficientList[Series[1/((1-6x)(1-8x)),{x,0,30}],x] (* or *) LinearRecurrence[{14,-48},{1,14},30] (* Harvey P. Dale, Dec 08 2011 *)
  • PARI
    Vec(1/((1-6*x)*(1-8*x))+O(x^99)) \\ Charles R Greathouse IV, Sep 24 2012
    
  • SageMath
    A016170=BinaryRecurrenceSequence(14,-48,1,14)
    [A016170(n) for n in range(31)] # G. C. Greubel, Nov 10 2024

Formula

a(n) = Sum_{k=1..n} 2^(n-1)*3^(n-k)*binomial(n,k). - Zerinvary Lajos, Sep 24 2006
From R. J. Mathar, Sep 18 2008: (Start)
a(n) = 4*8^n - 3*6^n = A081201(n+1).
Binomial transform of A081033. (End)
a(n) = 8*a(n-1) + 6^n. - Vincenzo Librandi, Feb 09 2011
a(0)=1, a(1)=14, a(n) = 14*a(n-1) - 48*a(n-2). - Harvey P. Dale, Dec 08 2011
E.g.f.: 4*exp(8*x) - 3*exp(6*x). - G. C. Greubel, Nov 10 2024

A081202 8th binomial transform of (0,1,0,1,0,1,....), A000035.

Original entry on oeis.org

0, 1, 16, 193, 2080, 21121, 206896, 1979713, 18640960, 173533441, 1602154576, 14701866433, 134294124640, 1222488408961, 11099284691056, 100571785292353, 909893629141120, 8222275592839681, 74233110849544336, 669726411243809473, 6038936596379658400, 54430221633714537601
Offset: 0

Views

Author

Paul Barry, Mar 11 2003

Keywords

Comments

Binomial transform of A081201.
From Wolfdieter Lang, Jul 17 2017: (Start)
For a combinatorial interpretation of a(n) with special 9-letter words of length n see the comment in A081200 on the 7-letter analog.
The binomial transform of {a(n)}_{n >=0} is A081203, the 10-letter analog.
(End)

Crossrefs

Apart from offset same as A016178.

Programs

Formula

a(n) = 16*a(n-1) - 63*a(n-2), a(0)=0, a(1)=1.
G.f.: x/((1-7*x)*(1-9*x)).
a(n) = (9^n - 7^n)/2.
E.g.f.: exp(7*x)*(exp(2*x) - 1)/2. - Stefano Spezia, Jul 23 2024

A081203 9th binomial transform of (0,1,0,1,0,1,.....), A000035.

Original entry on oeis.org

0, 1, 18, 244, 2952, 33616, 368928, 3951424, 41611392, 432891136, 4463129088, 45705032704, 465640261632, 4725122093056, 47800976744448, 482407813955584, 4859262511644672, 48874100093157376, 490992800745259008
Offset: 0

Views

Author

Paul Barry, Mar 11 2003

Keywords

Comments

Binomial transform of A081202.
From Wolfdieter Lang, Jul 17 2017: (Start)
For a combinatorial interpretation of a(n) with special 10-letter words of length n see the comment in A081200 on the 7-letter analog.
The binomial transform of {a(n)}_{n >= 0} is {0, A016190}, the 11-letter analog.
(End)

Crossrefs

Apart from the first term, identical to A016186.

Programs

  • Magma
    [10^n/2 - 8^n/2: n in [0..25]]; // Vincenzo Librandi, Aug 07 2013
  • Mathematica
    CoefficientList[Series[x / ((1 - 8 x) (1 - 10 x)), {x, 0, 20}], x] (* Vincenzo Librandi, Aug 07 2013 *)
    LinearRecurrence[{18,-80},{0,1},20] (* Harvey P. Dale, Aug 05 2018 *)

Formula

a(n) = 18*a(n-1)-80*a(n-2), a(0)=0, a(1)=1.
G.f.: x/((1-8*x)*(1-10*x)).
a(n) = 10^n/2 - 8^n/2.

A162590 Polynomials with e.g.f. exp(x*t)/csch(t), triangle of coefficients read by rows.

Original entry on oeis.org

0, 1, 0, 0, 2, 0, 1, 0, 3, 0, 0, 4, 0, 4, 0, 1, 0, 10, 0, 5, 0, 0, 6, 0, 20, 0, 6, 0, 1, 0, 21, 0, 35, 0, 7, 0, 0, 8, 0, 56, 0, 56, 0, 8, 0, 1, 0, 36, 0, 126, 0, 84, 0, 9, 0, 0, 10, 0, 120, 0, 252, 0, 120, 0, 10, 0, 1, 0, 55, 0, 330, 0, 462, 0, 165, 0, 11, 0, 0, 12, 0, 220, 0, 792, 0, 792, 0
Offset: 0

Views

Author

Peter Luschny, Jul 07 2009

Keywords

Comments

Comment from Peter Bala (Dec 06 2011): "Let P denote Pascal's triangle A070318 and put M = 1/2*(P-P^-1). M is A162590 (see also A131047). Then the first column of (I-t*M)^-1 (apart from the initial 1) lists the row polynomials for" A196776(n,k), which gives the number of ordered partitions of an n set into k odd-sized blocks. - Peter Luschny, Dec 06 2011
The n-th row of the triangle is formed by multiplying by 2^(n-1) the elements of the first row of the limit as k approaches infinity of the stochastic matrix P^(2k-1) where P is the stochastic matrix associated with the Ehrenfest model with n balls. The elements of a stochastic matrix P give the probability of arriving in a state j given the previous state i. In particular the sum of every row of the matrix must be 1, and so the sum of the terms in the n-th row of this triangle is 2^(n-1). Furthermore, by the properties of Markov chains, we can interpret P^(2k) as the (2k)-step transition matrix of the Ehrenfest model and its limit exists and it is again a stochastic matrix. The rows of the triangle divided by 2^(n-1) are the even rows (second, fourth, ...) and the odd rows (first, third, ...) of the limit matrix P^(2k). - Luca Onnis, Oct 29 2023

Examples

			Triangle begins:
  0
  1,  0
  0,  2,  0
  1,  0,  3,  0
  0,  4,  0,  4,  0
  1,  0, 10,  0,  5,  0
  0,  6,  0, 20,  0,  6,  0
  1,  0, 21,  0, 35,  0,  7,  0
  ...
  p[0](x) = 0;
  p[1](x) = 1
  p[2](x) = 2*x
  p[3](x) = 3*x^2 +  1
  p[4](x) = 4*x^3 +  4*x
  p[5](x) = 5*x^4 + 10*x^2 +  1
  p[6](x) = 6*x^5 + 20*x^3 +  6*x
  p[7](x) = 7*x^6 + 35*x^4 + 21*x^2 + 1
  p[8](x) = 8*x^7 + 56*x^5 + 56*x^3 + 8*x
.
Cf. the triangle of odd-numbered terms in rows of Pascal's triangle (A034867).
p[n] (k), n=0,1,...
k=0:  0, 1,  0,   1,    0,     1, ... A000035, (A059841)
k=1:  0, 1,  2,   4,    8,    16, ... A131577, (A000079)
k=2:  0, 1,  4,  13,   40,   121, ... A003462
k=3:  0, 1,  6,  28,  120,   496, ... A006516
k=4:  0, 1,  8,  49,  272,  1441, ... A005059
k=5:  0, 1, 10,  76,  520,  3376, ... A081199, (A016149)
k=6:  0, 1, 12, 109,  888,  6841, ... A081200, (A016161)
k=7:  0, 1, 14, 148, 1400, 12496, ... A081201, (A016170)
k=8:  0, 1, 16, 193, 2080, 21121, ... A081202, (A016178)
k=9:  0, 1, 18, 244, 2952, 33616, ... A081203, (A016186)
k=10: 0, 1, 20, 301, 4040, 51001, ... ......., (A016190)
.
p[n] (k), k=0,1,...
p[0]: 0,  0,   0,    0,    0,     0, ... A000004
p[1]: 1,  1,   1,    1,    1,     1, ... A000012
p[2]: 0,  2,   4,    6,    8,    10, ... A005843
p[3]: 1,  4,  13,   28,   49,    76, ... A056107
p[4]: 0,  8,  40,  120,  272,   520, ... A105374
p[5]: 1, 16, 121,  496, 1441,  3376, ...
p[6]: 0, 32, 364, 2016, 7448, 21280, ...
		

Crossrefs

Cf. A119467.

Programs

  • Maple
    # Polynomials: p_n(x)
    p := proc(n,x) local k;
    pow := (n,k) -> `if`(n=0 and k=0,1,n^k);
    add((k mod 2)*binomial(n,k)*pow(x,n-k),k=0..n) end;
    # Coefficients: a(n)
    seq(print(seq(coeff(i!*coeff(series(exp(x*t)/csch(t), t,16),t,i),x,n), n=0..i)), i=0..8);
  • Mathematica
    p[n_, x_] := Sum[Binomial[n, 2*k-1]*x^(n-2*k+1), {k, 0, n+2}]; row[n_] := CoefficientList[p[n, x], x] // Append[#, 0]&; Table[row[n], {n, 0, 12}] // Flatten (* Jean-François Alcover, Jun 28 2013 *)
    n = 15; "n-th row"
    mat = Table[Table[0, {j, 1, n + 1}], {i, 1, n + 1}];
    mat[[1, 2]] = 1;
    mat[[n + 1, n]] = 1;
    For[i = 2, i <= n, i++, mat[[i, i - 1]] = (i - 1)/n ];
    For[i = 2, i <= n, i++, mat[[i, i + 1]] = (n - i + 1)/n];
    mat // MatrixForm;
    P2 = Dot[mat, mat];
    R1 = Simplify[
      Eigenvectors[Transpose[P2]][[1]]/
       Total[Eigenvectors[Transpose[P2]][[1]]]]
    R2 = Table[Dot[R1, Transpose[mat][[k]]], {k, 1, n + 1}]
    even = R1*2^(n - 1) (* Luca Onnis, Oct 29 2023 *)

Formula

p_n(x) = Sum_{k=0..n} (k mod 2)*binomial(n,k)*x^(n-k).
E.g.f.: exp(x*t)/csch(t) = 0*(t^0/0!) + 1*(t^1/1!) + (2*x)*(t^2/2!) + (3*x^2+1)*(t^3/3!) + ...
The 'co'-polynomials with generating function exp(x*t)*sech(t) are the Swiss-Knife polynomials (A153641).

A102728 Array read by antidiagonals: T(n, k) = ((n+1)^k-(n-1)^k)/2.

Original entry on oeis.org

0, 0, 1, 0, 1, 0, 0, 1, 2, 1, 0, 1, 4, 4, 0, 0, 1, 6, 13, 8, 1, 0, 1, 8, 28, 40, 16, 0, 0, 1, 10, 49, 120, 121, 32, 1, 0, 1, 12, 76, 272, 496, 364, 64, 0, 0, 1, 14, 109, 520, 1441, 2016, 1093, 128, 1, 0, 1, 16, 148, 888, 3376, 7448, 8128, 3280, 256, 0, 0, 1, 18, 193, 1400, 6841
Offset: 0

Views

Author

Lambert Klasen (lambert.klasen(AT)gmx.net) and Gary W. Adamson, Feb 07 2005

Keywords

Comments

Consider a 2 X 2 matrix M = [N, 1] / [1, N]. The n-th row of the array contains the values of the non-diagonal elements of M^k, k=0,1,.... (Corresponding diagonal entry = non-diagonal entry + (N-1)^k.) Table:
N: row sequence g.f. cross references.
0: (1^n-(-1)^n)/2 x/((1+1x)(1-1x)) A000035
1: (2^n-0^n)/2 x/(1-2x) A000079
2: (3^n-1^n)/2 x/((1-1x)(1-3x)) A003462
3: (4^n-2^n)/2 x/((1-2x)(1-4x)) A006516
4: (7^n-3^n)/2 x/((1-3x)(1-5x)) A005059
5: (6^n-4^n)/2 x/((1-4x)(1-6x)) A016149
6: (7^n-5^n)/2 x/((1-5x)(1-7x)) A016161 A081200
7: (8^n-6^n)/2 x/((1-6x)(1-8x)) A016170 A081201
8: (9^n-7^n)/2 x/((1-7x)(1-9x)) A016178 A081202
9: (10^n-8^n)/2 x/((1-8x)(1-10x)) A016186 A081203
10: (11^n-9^n)/2 x/((1-9x)(1-11x)) A016190
11: (12^n-10^n)/2 x/((1-10x)(1-12x)) A016196
...
Characteristic polynomial x^2-2nx+(n^2-1) has roots n+-1, so if r(n) denotes a row sequence, r(n+1)/r(n) converges to n+1.
Columns follow polynomials with certain binomial coefficients:
column: polynomial
0: 0
1: 1
2: 2n
3: 3n^2+ 1 (see A056107)
4: 4n^3+ 4n (= 8*A006003(n))
5: 5n^4+ 10n^2+ 1
6: 6n^5+ 20n^3+ 6n
7: 7n^6+ 35n^4+ 21n^2+ 1
8; 8n^7+ 56n^5+ 56n^3+ 8n
9: 9n^8+ 84n^6+126n^4+ 36n^2+ 1
10: 10n^9+ 120n^7+252n^5+120n^3+ 10n
11: 11n^10+165n^8+462n^6+330n^4+ 55n^2+ 1

Examples

			Array begins:
0,1,0,1,0,1...
0,1,2,4,8,16...
0,1,4,13,40,121...
0,1,6,28,120,496...
0,1,8,49,272,1441...
...
		

Programs

  • PARI
    MM(n,N)=local(M);M=matrix(n,n);for(i=1,n, for(j=1,n,if(i==j,M[i,j]=N,M[i,j]=1)));M for(k=0,12, for(i=0,k,print1((MM(2,k-i)^i)[1,2],","))) T(n, k) = ((n+1)^k-(n-1)^k)/2 for(k=0,10, for(i=0,10,print1(T(k,i),","));print()) for(k=0,10, for(i=0,10,print1(((k+1)^i-(k-1)^i)/2,","));print()) for(k=0,10, for(i=0,10,print1(polcoeff(x/((1-(k-1)*x)*(1-(k+1)*x)),i),","));print())

A105373 Square array by antidiagonals of number of straight lines with n points in a k-dimensional hypercube with n points on each edge.

Original entry on oeis.org

1, 1, 1, 1, 6, 1, 1, 28, 8, 1, 1, 120, 49, 10, 1, 1, 496, 272, 76, 12, 1, 1, 2016, 1441, 520, 109, 14, 1, 1, 8128, 7448, 3376, 888, 148, 16, 1, 1, 32640, 37969, 21280, 6841, 1400, 193, 18, 1, 1, 130816, 192032, 131776, 51012, 12496, 2080, 244, 20, 1, 1, 523776
Offset: 1

Views

Author

Henry Bottomley, Apr 02 2005

Keywords

Examples

			Rows start:
  1,  1,   1,   1,    1,     1, ...;
  1,  6,  28, 120,  496,  2016, ...;
  1,  8,  49, 272, 1441,  7448, ...;
  1, 10,  76, 520, 3376, 21280, ...;
  1, 12, 109, 888, 6841, 51012, ...;
  etc.
T(5,3)=109 because in a 5 X 5 X 5 cube there are 25 columns, 25 linear rows in one direction, 25 linear rows in another direction, 5 short diagonals in each of 6 directions and 4 long diagonals; and 3*25 + 6*5 + 4 = 109.
		

Crossrefs

See A102728. Rows essentially include A000012, A006516, A005059, A016149 or A081199, A016161 or A081200, A016170 or A081201, A016178 or A081202 etc. Columns essentially include A000012, A005843, A056107, A105373.

Formula

T(1, k)=1. For n>1: T(n, k) = ((n+2)^k-n^k)/2 = (n+2)*T(n, k-1)+n^(k-1) = A102728(k, n+1).
Showing 1-8 of 8 results.