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-10 of 15 results. Next

A016160 Duplicate of A005062.

Original entry on oeis.org

0, 1, 11, 91, 671, 4651, 31031, 201811, 1288991, 8124571, 50700551
Offset: 0

Views

Author

Keywords

A005061 a(n) = 4^n - 3^n.

Original entry on oeis.org

0, 1, 7, 37, 175, 781, 3367, 14197, 58975, 242461, 989527, 4017157, 16245775, 65514541, 263652487, 1059392917, 4251920575, 17050729021, 68332056247, 273715645477, 1096024843375, 4387586157901, 17560804984807, 70274600998837, 281192547174175, 1125052618233181
Offset: 0

Views

Author

Keywords

Comments

Number of 2 X n binary arrays with a path of adjacent 1's from top row to bottom row, see A359576. - R. H. Hardin, Mar 21 2002
Number of binary vectors (x_1, x_2, ..., x_{2n}) such that in at least one of the disjoint pairs (x_1, x_2), (x_3, x_4), ..., (x_{2n-1}, x_{2n}) both x_{2i-1} and x_{2i} are both 1. Equivalently, number of solutions (x_1, ..., x_n) to the equation x_1*x_2 + x_3*x_4 + x_5*x_6 + ... +x_{2n-1}*x_{2n} = 1 in base-2 lunar arithmetic. - N. J. A. Sloane, Apr 23 2011
a(n)/4^n is the probability that two randomly selected (with replacement) subsets of [n] will have at least one element in common if the probability of selection is equal for all subsets. - Geoffrey Critzer, May 09 2009
This sequence is also the second column of the Sheffer triangle A143495 (3-restricted Stirling2 numbers). (See the e.g.f. given below.) - Wolfdieter Lang, Oct 08 2011
Also, the number of numbers with at most n digits whose largest digit equals 3. See A255463 for the first differences (i.e., ...with exactly n digits...). - M. F. Hasler, May 03 2015
If 2^k | n then a(2^k) | a(n). - Bernard Schott, Oct 08 2020
a(n) is the number of ordered n-tuples with elements from {0,1,2,3} in which any of these elements, say 0, appears at least once. For example, a(2)=7 since 01,10,02,20,03,30,00 are the ordered 2-tuples that contain 0. - Enrique Navarrete, Apr 05 2021
a(n) is the number of n-digit numbers whose smallest decimal digit is 6. - Stefano Spezia, Nov 15 2023

Examples

			G.f. = x + 7*x^2 + 37*x^3 + 175*x^4 + 781*x^5 + 3367*x^6 + 14197*x^7 + ...
		

Crossrefs

Cf. A001047, A002250, A005060, A005062, A143495, A255463 (first differences), A359576.
Array column A047969(n-1, 3), or triangle's subdiagonal A047969(n+2, n-1), for n >= 1.

Programs

Formula

a(n) = 4*a(n-1) + 3^(n-1) for n>=1. - Xavier Acloque, Oct 20 2003
Binomial transform of A001047. - Ross La Haye, Sep 17 2005
From Mohammad K. Azarian, Jan 14 2009: (Start)
G.f.: 1/(1-4*x)-1/(1-3*x).
E.g.f.: exp(4*x)-exp(3*x). (End)
a(n) = 2^n * Sum_{i=0...n} binomial(n,i)*(2^i-1)/2^i. - Geoffrey Critzer, May 09 2009
a(n) = 7*a(n-1) - 12*a(n-2) for n>=2. - Bruno Berselli, Jan 25 2011
From Joe Slater, Jan 15 2017: (Start)
a(n) = 3*a(n-1) + 4^(n-1) for n>=0.
a(n+1) = Sum_{k=0..n} 4^(n-k) * 3^k. (End)
a(n) = -a(-n) * 12^n for all n in Z. - Michael Somos, Jan 22 2017

A047969 Square array of nexus numbers a(n,k) = (n+1)^(k+1) - n^(k+1) (n >= 0, k >= 0) read by upwards antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 5, 7, 1, 1, 7, 19, 15, 1, 1, 9, 37, 65, 31, 1, 1, 11, 61, 175, 211, 63, 1, 1, 13, 91, 369, 781, 665, 127, 1, 1, 15, 127, 671, 2101, 3367, 2059, 255, 1, 1, 17, 169, 1105, 4651, 11529, 14197, 6305, 511, 1, 1, 19, 217, 1695, 9031
Offset: 0

Views

Author

Keywords

Comments

If each row started with an initial 0 (i.e., a(n,k) = (n+1)^k - n^k) then each row would be the binomial transform of the preceding row. - Henry Bottomley, May 31 2001
a(n-1, k-1) is the number of ordered k-tuples of positive integers such that the largest of these integers is n. - Alford Arnold, Sep 07 2005
From Alford Arnold, Jul 21 2006: (Start)
The sequences in A047969 can also be calculated using the Eulerian Array (A008292) and Pascal's Triangle (A007318) as illustrated below: (cf. A101095).
1 1 1 1 1 1
1 1 1 1 1 1
-----------------------------------------
1 2 3 4 5 6
1 2 3 4 5
1 3 5 7 9 11
-----------------------------------------
1 3 6 10 15 21
4 12 24 40 60
1 3 6 10
1 7 19 37 61 91
-----------------------------------------
1 4 10 20 35 56
11 44 110 220 385
11 44 110 220
1 4 10
1 15 65 175 369 671
----------------------------------------- (End)
From Peter Bala, Oct 26 2008: (Start)
The above remarks of Alford Arnold may be summarized by saying that (the transpose of) this array is the Hilbert transform of the triangle of Eulerian numbers A008292 (see A145905 for the definition of the Hilbert transform). In this context, A008292 is best viewed as the array of h-vectors of permutohedra of type A. See A108553 for the Hilbert transform of the array of h-vectors of type D permutohedra. Compare this array with A009998.
The polynomials n^k - (n-1)^k, k = 1,2,3,..., which give the nonzero entries in the columns of this array, satisfy a Riemann hypothesis: their zeros lie on the vertical line Re s = 1/2 in the complex plane. See A019538 for the connection between the polynomials n^k - (n-1)^k and the Stirling polynomials of the simplicial complexes dual to the type A permutohedra.
(End)
Empirical: (n+1)^(k+1) - n^(k+1) is the number of first differences of length k+1 arrays of numbers in 0..n, k > 0. - R. H. Hardin, Jun 30 2013
a(n-1, k-1) is the number of bargraphs of width k and height n. Examples: a(1,2) = 7 because we have [1,1,2], [1,2,1], [2,1,1], [1,2,2], [2,1,2], [2,2,1], and [2,2,2]; a(2,1) = 5 because we have [1,3], [2,3], [3,1], [3,2], and [3,3] (bargraphs are given as compositions). This comment is equivalent to A. Arnold's Sep 2005 comment. - Emeric Deutsch, Jan 30 2017

Examples

			Array a begins:
  [n\k][0  1   2    3    4   5  6  ...
  [0]   1  1   1    1    1   1  1  ...
  [1]   1  3   7   15   31  63  ...
  [2]   1  5  19   65  211  ...
  [3]   1  7  37  175  ...
  ...
Triangle T begins:
  n\m   0   1    2     3     4      5      6      7      8     9  10 ...
  0:    1
  1:    1   1
  2:    1   3    1
  3:    1   5    7     1
  4:    1   7   19    15     1
  5:    1   9   37    65    31      1
  6:    1  11   61   175   211     63      1
  7:    1  13   91   369   781    665    127      1
  8:    1  15  127   671  2101   3367   2059    255      1
  9:    1  17  169  1105  4651  11529  14197   6305    511     1
  10:   1  19  217  1695  9031  31031  61741  58975  19171  1023   1
  ...  - _Wolfdieter Lang_, May 07 2021
		

References

  • J. H. Conway and R. K. Guy, The Book of Numbers, Copernicus Press, NY, 1996, p. 54.

Crossrefs

Cf. A047970.
Cf. A009998, A108553 (Hilbert transform of array of h-vectors of type D permutohedra), A145904, A145905.
Row n sequences of array a: A000012, A000225(k+1), A001047(k+1), A005061(k+1), A005060(k+1), A005062(k+1), A016169(k+1), A016177(k+1), A016185(k+1), A016189(k+1), A016195(k+1), A016197(k+1).
Column k sequences of array a: (nexus numbers): A000012, A005408, A003215, A005917(n+1), A022521, A022522, A022523, A022524, A022525, A022526, A022527, A022528.
Cf. A343237 (row reversed triangle).

Programs

  • Mathematica
    Flatten[Table[n = d - e; k = e; (n + 1)^(k + 1) - n^(k + 1), {d, 0, 100}, {e, 0, d}]] (* T. D. Noe, Feb 22 2012 *)
  • Maxima
    T(n,m):=if m=0 then 1 else sum(k!*(-1)^(m+k)*stirling2(m,k)*binomial(n+k-1,n),k,0,m); /* Vladimir Kruchinin, Jan 28 2018 */

Formula

From Vladimir Kruchinin: (Start)
O.g.f. of e.g.f of rows of array: ((1-x)*exp(y))/(1-x*exp(y))^2.
T(n,m) = Sum_{k=0..m} k!*(-1)^(m+k)*Stirling2(m,k)*C(n+k-1,n), T(n,0)=1.(End)
From Wolfdieter Lang, May 07 2021: (Start)
T(n,m) = a(n-m,m) = (n-m+1)^(m+1) - (n-m)^(m+1), n >= 0, m = 0, 1,..., n.
O.g.f. column k of the array: polylog(-(k+1), x)*(1-x)/x. See the Peter Bala comment above, and the Eulerian triangle A008292 formula by Vladeta Jovovic, Sep 02 2002.
E.g.f. of e.g.f. of row of the array: exp(y)*(1 + x*(exp(y) - 1))*exp(x*exp(y)).
O.g.f. of triangle's exponential row polynomials R(n, y) = Sum_{m=0} T(n, m)*(y^m)/m!: G(x, y) = exp(x*y)*(1 - x)/(1 - x*exp(x*y))^2. (End)

A016129 Expansion of 1/((1-2*x)*(1-6*x)).

Original entry on oeis.org

1, 8, 52, 320, 1936, 11648, 69952, 419840, 2519296, 15116288, 90698752, 544194560, 3265171456, 19591036928, 117546237952, 705277460480, 4231664828416, 25389989101568, 152339934871552, 914039609753600, 5484237659570176, 32905425959518208, 197432555761303552
Offset: 0

Views

Author

Keywords

Crossrefs

Row sums of A100851.
Sequences with gf 1/((1-n*x)*(1-6*x)): A000400 (n=0), A003464 (n=1), this sequence (n=2), A016137 (n=3), A016149 (n=4), A005062 (n=5), A053469 (n=6), A016169 (n=7), A016170 (n=8), A016172 (n=9), A016173 (n=10), A016174 (n=11), A016175 (n=12).

Programs

Formula

a(n) = A071951(n+2, 2) = 9*(2*3)^(n-1) - (2*1)^(n-1) = (2^(n-1))*(3^(n+1)-1), n>=0. - Wolfdieter Lang, Nov 07 2003
From Lambert Klasen (lambert.klasen(AT)gmx.net), Feb 05 2005: (Start)
G.f.: 1/((1-2*x)*(1-6*x)).
E.g.f.: (-exp(2*x) + 3*exp(6*x))/2.
a(n) = (6^(n+1) - 2^(n+1))/4. (End)
a(n)^2 = A144843(n+1). - Philippe Deléham, Nov 26 2008
a(n) = 8*a(n-1) - 12*a(n-2). - Philippe Deléham, Jan 01 2009
a(n) = det(|ps(i+2,j+1)|, 1 <= i,j <= n), where ps(n,k) are Legendre-Stirling numbers of the first kind (A129467). - Mircea Merca, Apr 06 2013

A016189 a(n) = 10^n - 9^n.

Original entry on oeis.org

0, 1, 19, 271, 3439, 40951, 468559, 5217031, 56953279, 612579511, 6513215599, 68618940391, 717570463519, 7458134171671, 77123207545039, 794108867905351, 8146979811148159, 83322818300333431, 849905364703000879, 8649148282327007911, 87842334540943071199, 890581010868487640791
Offset: 0

Views

Author

Keywords

Comments

Almost all numbers contain any given sequence of digits (in any base) [Theorem 143 of Hardy and Wright]. a(7) = 5217031, more than 52% of the numbers < 10^7 contain any given nonzero decimal digit. - Frank Ellermann, May 30 2001
a(n) gives the number of integers from 0 to 10^n-1 which contain (at least) any one given decimal digit except 0. - Michael Taktikos, Aug 24 2004
These are the numerators of a(n)=(integral{x=0 to 0.2} (1-0.5*x)^n dx). E.g., a(3)=3439/20000. The denominators are b(n)=5*(n+1)*10^n. E.g., b(3)=20000. - Al Hakanson (hawkuu(AT)excite.com), Feb 22 2004
Binomial transforms of sequences defined by a(n)=(C+1)^n-C^n are the sequences (C+2)^n-(C+1)^n. The binomial transform of this here is in A016195, for example. - R. J. Mathar, Nov 27 2008
First differences are given in A088924. - M. F. Hasler, May 04 2015

References

  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 5th ed., Oxford Univ. Press, 1979, th. 143

Crossrefs

Base 2: A000225, 3: A001047, 4: A005061, 5: A005060, 6: A005062, base 7: A016169, 8: A016177, 9: A016185 11: A016195 12: A016197.
Equals A155671 - 1.

Programs

Formula

G.f.: x/((1-9x)(1-10x)).
a(0) = 0, a(1) = 1, then a(n+1) = 9*a(n) + 10^n.
a(n) = 19*a(n-1) - 90*a(n-2), n > 1; a(0)=0, a(1)=1. - Philippe Deléham, Jan 01 2009
E.g.f.: e^(10*x) - e^(9*x). - Mohammad K. Azarian, Jan 14 2009

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

A193685 5-Stirling numbers of the second kind.

Original entry on oeis.org

1, 5, 1, 25, 11, 1, 125, 91, 18, 1, 625, 671, 217, 26, 1, 3125, 4651, 2190, 425, 35, 1, 15625, 31031, 19981, 5590, 740, 45, 1, 78125, 201811, 170898, 64701, 12250, 1190, 56, 1, 390625, 1288991, 1398097, 688506, 174951, 24150, 1806, 68, 1, 1953125, 8124571, 11075670, 6906145, 2263065, 416451, 44016, 2622, 81, 1, 9765625, 50700551, 85654261, 66324830, 27273730, 6427575, 900627, 75480, 3675, 95, 1
Offset: 0

Views

Author

Wolfdieter Lang, Oct 06 2011

Keywords

Comments

This is the lower triangular Sheffer matrix (exp(5*x),exp(x)-1). For Sheffer matrices see the W. Lang link under A006232 with references, and the rules for the conversion to the umbral notation of S. Roman's book.
The general case is Sheffer (exp(r*x),exp(x)-1), r=0,1,..., corresponding to r-Stirling2 numbers with row and column offsets 0. See the Broder link for r-Stirling2 numbers with offset [r,r].
a(n,m), n >= m >= 0, gives the number of partitions of the set {1.2....,n+5} into m+5 nonempty distinct subsets such that 1,2,3,4 and 5 belong to distinct subsets.
a(n,m) appears in the following normal ordering of Bose operators a and a* satisfying the Lie algebra [a,a*]=1: (a*a)^n (a*)^5 = Sum_{m=0..n} a(n,m)*(a*)^(5+m)*a^m, n >= 0. See the Mikhailov papers, where a(n,m) = S(n+5,m+5,5).
With a->D=d/dx and a*->x we also have
(xD)^n x^5 = Sum_{m=0..n} a(n,m)*x^(5+m)*D^m, n >= 0.

Examples

			n\m  0    1    2   3  4  5 ...
0    1
1    5    1
2   25   11    1
3  125   91   18   1
4  625  671  217  26  1
5 3125 4651 2190 425 35  1
...
5-restricted S2: a(1,0)=5 from 1,6|2|3|4|5, 2,6|1|3|4|5,
3,6|1|2|4|5, 4,6|1|2|3|5 and 5,6|1|2|3|4.
Recurrence: a(4,2) = (5+2)*a(3,2)+ a(3,1) = 7*18 + 91 = 217.
Normal ordering (n=1): (xD)^1 x^5 = Sum_{m=0..1} a(1,m)*x^(5+m)*D^m = 5*x^5 + 1*x^6*D.
a(2,1) = Sum_{j=0..1} S1(5,5-j)*S2(7-j,6) = 1*21 - 10*1 = 11.
		

Crossrefs

Cf. A196834 (row sums), A196835 (alternating row sums).
Columns: A000351 (m=0), A005062 (m=1), A019757 (m=2), A028165 (m=3), ...

Programs

  • Mathematica
    a[n_, m_] := Sum[ StirlingS1[5, 5-j]*StirlingS2[n+5-j, m+5], {j, 0, Min[5, n-m]}]; Flatten[ Table[ a[n, m], {n, 0, 10}, {m, 0, n}] ] (* Jean-François Alcover, Dec 02 2011, after Wolfdieter Lang *)

Formula

E.g.f. of row polynomials s(n,x):=Sum_{m=0..n} a(n,m)*x^m: exp(5*z + x(exp(z)-1)).
E.g.f. of column no. m (with leading zeros):
exp(5*x)*((exp(x)-1)^m)/m!, m >= 0 (Sheffer).
O.g.f. of column no. m (without leading zeros):
1/Product_{j=0..m} (1-(5+j)*x), m >= 0. (Compute the first derivative of the column e.g.f. and compare its Laplace transform with the partial fraction decomposition of the o.g.f. x^(m-1)/Product_{j=0..m} (1-(5+j)*x). This works for every r-restricted Stirling2 triangle.)
Recurrence: a(n,m) = (5+m)*a(n-1,m) + a(n-1,m-1), a(0,0)=1, a(n,m)=0 if n < m, a(n,-1)=0.
a(n,m) = Sum_{j=0..min(5,n-m)} S1(5,5-j)*S2(n+5-j,m+5), n >= m >= 0, with S1 and S2 the Stirling1 and Stirling2 numbers A008275 and A048993, respectively (see the Mikailov papers).
Dobinski-type formula for the row polynomials: R(n,x) = exp(-x)*Sum_{k>=0} k*(4+k)^(n-1)*x^(k-1)/k!. - Peter Bala, Jun 23 2014

A257286 a(n) = 5*6^n - 4*5^n.

Original entry on oeis.org

1, 10, 80, 580, 3980, 26380, 170780, 1087180, 6835580, 42575980, 263268380, 1618672780, 9907349180, 60420657580, 367406757980, 2228854610380, 13495197974780, 81581539411180, 492540994279580, 2970504754739980, 17899322473752380
Offset: 0

Views

Author

M. F. Hasler, May 03 2015

Keywords

Comments

First differences of 6^n - 5^n = A005062.
a(n-1) is the number of numbers with n digits having the largest digit equal to 5. Or, equivalently, number of n-letter words over a 6-letter alphabet {a,b,c,d,e,f}, which must not start with the first letter of the alphabet, and in which the last letter of the alphabet must appear.

Crossrefs

Cf. A005062.
Coincides with A125373 only for the first terms.

Programs

  • Magma
    [5*6^n-4*5^n: n in [0..20]]; // Vincenzo Librandi, May 04 2015
  • Mathematica
    Table[5 6^n - 4 5^n, {n, 0, 30}] (* Vincenzo Librandi, May 04 2015 *)
  • PARI
    a(n)=5*6^n-4*5^n
    

Formula

a(n) = 11 a(n-1) - 30 a(n-2).
G.f.: (1-x)/((1-5*x)*(1-6*x)). - Vincenzo Librandi, May 04 2015
E.g.f.: exp(5*x)*(5*exp(x) - 4). - Stefano Spezia, Nov 15 2023

A028165 Expansion of 1/((1-5x)*(1-6x)*(1-7x)*(1-8x)).

Original entry on oeis.org

1, 26, 425, 5590, 64701, 688506, 6906145, 66324830, 616252901, 5580303586, 49508360265, 432061044870, 3720287489101, 31681154472266, 267320885100785, 2238337148081710, 18621251375573301, 154069635600426546
Offset: 0

Views

Author

Keywords

Comments

This is the column m=2 sequence (without leading zeros) of the Sheffer triangle (exp(5*x), exp(x)-1) of the 5-restricted Stirling2 numbers A193685. For a proof see the column o.g.f. formula there. - Wolfdieter Lang, Oct 07 2011

Crossrefs

Programs

Formula

If we define f(m,j,x) = Sum_{k=j..m} binomial(m,k)*Stirling2(k,j)*x^(m-k) then a(n-3) = f(n,3,5), (n >= 3). - Milan Janjic, Apr 26 2009
a(n) = 26*a(n-1) - 251*a(n-2) + 1066*a(n-3) - 1680*a(n-4), n >= 4. - Vincenzo Librandi, Mar 19 2011
a(n) = 15*a(n-1) - 56*a(n-2) + 6^(n+1) - 5^(n+1), a(0)=1, a(1)=26. - Vincenzo Librandi, Mar 19 2011
E.g.f.: (d^3/dx^3)(exp(5*x)*((exp(x)-1)^3)/3!). See the Sheffer triangle comment above. - Wolfdieter Lang, Oct 07 2011
a(n) = -125*5^n/6 + 108*6^n - 343*7^n/2 + 256*8^n/3. - R. J. Mathar, Jun 23 2013

A087749 Primes p such that 6^p - 5^p is composite.

Original entry on oeis.org

3, 7, 17, 19, 29, 31, 37, 41, 43, 47, 53, 59, 67, 71, 73, 79, 89, 97, 101, 103, 107, 109, 113, 127, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 241, 251, 257, 263, 269, 271, 277, 281
Offset: 1

Views

Author

Cino Hilliard, Oct 26 2003

Keywords

Crossrefs

Cf. A005062.

Programs

  • Mathematica
    Select[Prime[Range[60]],!PrimeQ[6^#-5^#]&] (* Harvey P. Dale, Apr 26 2011 *)
  • PARI
    apmb(a,b,n) = { forprime(x=2,n, y=a^x-b^x; if(!ispseudoprime(y), print1(x","); ) ) }

Extensions

Corrected by Andrew S. Plewe, Apr 24 2007 (all the terms were wrong)
Checked by Ray Chandler, Apr 25 2007
Showing 1-10 of 15 results. Next