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 23 results. Next

A099265 Partial sums of A056272.

Original entry on oeis.org

1, 3, 8, 23, 75, 277, 1132, 4977, 22979, 109451, 531456, 2610931, 12917683, 64181625, 319695980, 1594859885, 7963472187, 39784944799, 198827606704, 993846943839, 4968361974491, 24839192686973, 124188113975628, 620917025694793, 3104514504312595, 15522360665856147, 77611167795714752
Offset: 1

Views

Author

Nelma Moreira, Oct 10 2004

Keywords

Comments

Density of regular language L{0}* over {0, 1, 2, 3, 4, 5} (i.e., the number of strings of length n), where L is described by regular expression with c = 5: Sum_{i=1..c} Product_{j=1..i} (j(1+...+j)*), where "Sum" stands for union and "Product" for concatenation. I.e., L = L((11* + 11*2(1 + 2)* + ... + 11*2(1 + 2)*3(1 + 2 + 3)*4(1 + 2 + 3 + 4)*5(1 + 2 + 3 + 4 + 5)*)0*).

Crossrefs

Programs

  • Maple
    with (combinat):seq(sum(sum(stirling2(k, j),j=1..5), k=1..n), n=1..23); # Zerinvary Lajos, Dec 04 2007
  • PARI
    a(n) = sum(m=1, n, sum(i=1, 5, stirling(m, i, 2))) \\ Petros Hadjicostas, Mar 10 2021

Formula

a(5,n) = (1/96)*5^n + (1/8)*3^n + (1/3)*2^n + (3/8)*n - 15/32.
a(n) = Sum_{m=1..n} Sum_{i=1..5} S(m,i), where S(m,i) = A008277(m,i) (i.e., partial sum of the sum of Stirling numbers of second kind S(n,i) for i = 1..5).
For c = 5, a(c,n) = g(1,c)*n + Sum_{k=2..c} g(k,c)*k*(k^n - 1)/(k - 1), where g(1,1) = 1, g(1,c) = g(1,c-1) + (-1)^(c-1)/(c-1)! for c > 1, and g(k,c) = g(k-1, c-1)/k for c > 1 and 2 <= k <= c.
G.f.: x*(-1 + 19*x^3 - 24*x^2 + 9*x)/((3*x-1)*(2*x-1)*(5*x-1)*(x-1)^2). [Maksym Voznyy (voznyy(AT)mail.ru), Jul 28 2009]

Extensions

Name and Formula section edited by Petros Hadjicostas, Mar 10 2021
More terms from Michel Marcus, Jan 05 2025

A007051 a(n) = (3^n + 1)/2.

Original entry on oeis.org

1, 2, 5, 14, 41, 122, 365, 1094, 3281, 9842, 29525, 88574, 265721, 797162, 2391485, 7174454, 21523361, 64570082, 193710245, 581130734, 1743392201, 5230176602, 15690529805, 47071589414, 141214768241, 423644304722, 1270932914165, 3812798742494, 11438396227481
Offset: 0

Views

Author

Keywords

Comments

Number of ordered trees with n edges and height at most 4.
Number of palindromic structures using a maximum of three different symbols. - Marks R. Nester
Number of compositions of all even natural numbers into n parts <= 2 (0 is counted as a part), see example. - Adi Dani, May 14 2011
Consider the mapping f(a/b) = (a + 2*b)/(2*a + b). Taking a = 1, b = 2 to start with, and carrying out this mapping repeatedly on each new (reduced) rational number gives the sequence 1/2, 4/5, 13/14, 40/41, ... converging to 1. The sequence contains the denominators = (3^n+1)/2. The same mapping for N, i.e., f(a/b) = (a + N*b)/(a+b) gives fractions converging to N^(1/2). - Amarnath Murthy, Mar 22 2003
Second binomial transform of the expansion of cosh(x). - Paul Barry, Apr 05 2003
The sequence (1, 1, 2, 5, ...) = 3^n/6 + 1/2 + 0^n/3 has binomial transform A007581. - Paul Barry, Jul 20 2003
Number of (s(0), s(1), ..., s(2n+2)) such that 0 < s(i) < 6 and |s(i) - s(i-1)| = 1 for i = 1, 2, ..., 2n+2, s(0) = 1, s(2n+2) = 1. - Herbert Kociemba, Jun 10 2004
Density of regular language L over {1,2,3}^* (i.e., number of strings of length n in L) described by regular expression 11*+11*2(1+2)*+11*2(1+2)*3(1+2+3)*. - Nelma Moreira, Oct 10 2004
Sums of rows of the triangle in A119258. - Reinhard Zumkeller, May 11 2006
Number of n-words from the alphabet A = {a,b,c} which contain an even number of a's. - Fung Cheok Yin (cheokyin_restart(AT)yahoo.com.hk), Aug 30 2006
Let P(A) be the power set of an n-element set A. Then a(n) = the number of pairs of elements {x,y} of P(A) for which either 0) x and y are disjoint and for which x is not a subset of y and y is not a subset of x, or 1) x = y. - Ross La Haye, Jan 10 2008
a(n+1) gives the number of primitive periodic multiplex juggling sequences of length n with base state <2>. - Steve Butler, Jan 21 2008
a(n) is also the number of idempotent order-preserving and order-decreasing partial transformations (of an n-chain). - Abdullahi Umar, Oct 02 2008
Equals row sums of triangle A147292. - Gary W. Adamson, Nov 05 2008
Equals leftmost column of A071919^3. - Gary W. Adamson, Apr 13 2009
A010888(a(n))=5 for n >= 2, that is, the digital root of the terms >= 5 equals 5. - Parthasarathy Nambi, Jun 03 2009
Let A be the Hessenberg matrix of order n, defined by: A[1,j]=1, A[i,i]:=5, (i>1), A[i,i-1]=-1, and A[i,j]=0 otherwise. Then, for n>=1, a(n-1)=(-1)^n*charpoly(A,2). - Milan Janjic, Jan 27 2010
Let A be the Hessenberg matrix of order n, defined by: A[1,j]=1, A[i,i]:=6, (i>1), A[i,i-1]=-1, and A[i,j]=0 otherwise. Then, for n>=1, a(n)=(-1)^(n-1)*charpoly(A,3). - Milan Janjic, Feb 21 2010
It appears that if s(n) is a rational sequence of the form s(1)=2, s(n)= (2*s(n-1)+1)/(s(n-1)+2), n>1 then s(n)=a(n)/(a(n-1)-1).
Form an array with m(1,n)=1 and m(i,j) = Sum_{k=1..i-1} m(k,j) + Sum_{k=1..j-1} m(i,k), which is the sum of the terms to the left of m(i,j) plus the sum above m(i,j). The sum of the terms in antidiagonal(n-1) = a(n). - J. M. Bergot, Jul 16 2013
From Peter Bala, Oct 29 2013: (Start)
An Engel expansion of 3 to the base b := 3/2 as defined in A181565, with the associated series expansion 3 = b + b^2/2 + b^3/(2*5) + b^4/(2*5*14) + .... Cf. A034472.
More generally, for a positive integer n >= 3, the sequence [1, n - 1, n^2 - n - 1, ..., ( (n - 2)*n^k + 1 )/(n - 1), ...] is an Engel expansion of n/(n - 2) to the base n/(n - 1). Cases include A007583 (n = 4), A083065 (n = 5) and A083066 (n = 6). (End)
Diagonal elements (and one more than antidiagonal elements) of the matrix A^n where A=(2,1;1,2). - David Neil McGrath, Aug 17 2014
From M. Sinan Kul, Sep 07 2016: (Start)
a(n) is equal to the number of integer solutions to the following equation when x is equal to the product of n distinct primes: 1/x = 1/y + 1/z where 0 < x < y <= z.
If z = k*y where k is a fraction >= 1 then the solutions can be given as: y = ((k+1)/k)*x and z = (k+1)*x.
Here k can be equal to any divisor of x or to the ratio of two divisors.
For example for x = 2*3*5 = 30 (product of three distinct primes), k would have the following 14 values: 1, 6/5, 3/2, 5/3, 2, 5/2, 3, 10/3, 5, 6, 15/2, 10, 15, 30.
As an example for k = 10/3, we would have y=39, z=130 and 1/39 + 1/130 = 1/30.
Here finding the number of fractions would be equivalent to distributing n balls (distinct primes) to two bins (numerator and denominator) with no empty bins which can be found using Stirling numbers of the second kind. So another definition for a(n) is: a(n) = 2^n + Sum_{i=2..n} Stirling2(i,2)*binomial(n,i).
(End)
a(n+1) is the smallest i for which the Catalan number C(i) (see A000108) is divisible by 3^n for n > 0. This follows from the rule given by Franklin T. Adams-Watters for determining the multiplicity with which a prime divides C(n). We need to find the smallest number in base 3 to achieve a given count. Applied to prime 3, 1 is the smallest digit that counts but requires to be followed by 2 which cannot be at the end to count. Therefore the number in base 3 of the form 1{n-1 times}20 = (3^(n+1) + 1)/2 + 1 = a(n+1)+1 is the smallest number to achieve count n which implies the claim. - Peter Schorn, Mar 06 2020
Let A be a Toeplitz matrix of order n, defined by: A[i,j]=1, if ij; A[i,i]=2. Then, for n>=1, a(n) = det A. - Dmitry Efimov, Oct 28 2021
a(n) is the least number k such that A065363(k) = -(n-1), for n > 0. - Amiram Eldar, Sep 03 2022

Examples

			From _Adi Dani_, May 14 2011: (Start)
a(3)=14 because all compositions of even natural numbers into 3 parts <=2 are
for 0: (0,0,0)
for 2: (0,1,1), (1,0,1), (1,1,0), (0,0,2), (0,2,0), (2,0,0)
for 4: (0,2,2), (2,0.2), (2,2,0), (1,1,2), (1,2,1), (2,1,1)
for 6: (2,2,2).
(End)
		

References

  • J. M. Borwein, D. H. Bailey and R. Girgensohn, Experimentation in Mathematics, A K Peters, Ltd., Natick, MA, 2004. x+357 pp. See p. 47.
  • Adi Dani, Quasicompositions of natural numbers, Proceedings of III congress of mathematicians of Macedonia, Struga Macedonia 29 IX -2 X 2005 pages 225-238.
  • R. K. Guy, Unsolved Problems in Number Theory, Springer, 1st edition, 1981. See section E11.
  • M. R. Nester (1999). Mathematical investigations of some plant interaction designs. PhD Thesis. University of Queensland, Brisbane, Australia. [See A056391 for pdf file of Chap. 2]
  • P. Ribenboim, The Book of Prime Number Records. Springer-Verlag, NY, 2nd ed., 1989, p. 60.
  • P. Ribenboim, The Little Book of Big Primes, Springer-Verlag, NY, 1991, p. 53.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

Formula

a(n) = 3*a(n-1) - 1.
Binomial transform of Chebyshev coefficients A011782. - Paul Barry, Mar 16 2003
From Paul Barry, Mar 16 2003: (Start)
a(n) = 4*a(n-1) - 3*a(n-2) for n > 1, a(0)=1, a(1)=2.
G.f.: (1 - 2*x)/((1 - x)*(1 - 3*x)). (End)
E.g.f.: exp(2*x)*cosh(x). - Paul Barry, Apr 05 2003
a(n) = Sum_{k=0..floor(n/2)} binomial(n, 2*k)*2^(n-2*k). - Paul Barry, May 08 2003
This sequence is also the partial sums of the first 3 Stirling numbers of second kind: a(n) = S(n+1, 1) + S(n+1, 2) + S(n+1, 3) for n >= 0; alternatively it is the number of partitions of [n+1] into 3 or fewer parts. - Mike Zabrocki, Jun 21 2004
For c=3, a(n) = (c^n)/c! + Sum_{k=1..c-2}((k^n)/k!*(Sum_{j=2..c-k}(((-1)^j)/j!))) or = Sum_{k=1..c} g(k, c)*k^n where g(1, 1) = 1, g(1, c) = g(1, c-1) + ((-1)^(c-1))/(c-1)! for c > 1, and g(k, c) = g(k-1, c-1)/k for c > 1 and 2 <= k <= c. - Nelma Moreira, Oct 10 2004
The i-th term of the sequence is the entry (1, 1) in the i-th power of the 2 X 2 matrix M = ((2, 1), (1, 2)). - Simone Severini, Oct 15 2005
If p[i]=fibonacci(2i-3) and if A is the Hessenberg matrix of order n defined by: A[i,j]=p[j-i+1], (i<=j), A[i,j]=-1, (i=j+1), and A[i,j]=0 otherwise. Then, for n>=1, a(n-1)= det A. - Milan Janjic, May 08 2010
INVERT transform of A001519: [1, 1, 2, 5, 13, 34, ...]. - Gary W. Adamson, Jun 13 2011
a(n) = M^n*[1,1,1,0,0,0,...], leftmost column term; where M = an infinite bidiagonal matrix with all 1's in the superdiagonal and (1,2,3,...) in the main diagonal and the rest zeros. - Gary W. Adamson, Jun 23 2011
a(n) = M^n*[1,1,1,0,0,0,...], top entry term; where M is an infinite bidiagonal matrix with all 1's in the superdiagonal, (1,2,3,...) as the main diagonal, and the rest zeros. - Gary W. Adamson, Jun 24 2011
a(n) = A201730(n,0). - Philippe Deléham, Dec 05 2011
a(n) = A006342(n) + A006342(n-1). - Yuchun Ji, Sep 19 2018
From Dmitry Efimov, Oct 29 2021: (Start)
a(2*m+1) = Product_{k=-m..m} (2+i*tan(Pi*k/(2*m+1))),
a(2*m) = Product_{k=-m..m-1} (2+i*tan(Pi*(2*k+1)/(4*m))),
where i is the imaginary unit. (End)

A007581 a(n) = (2^n+1)*(2^n+2)/6.

Original entry on oeis.org

1, 2, 5, 15, 51, 187, 715, 2795, 11051, 43947, 175275, 700075, 2798251, 11188907, 44747435, 178973355, 715860651, 2863377067, 11453377195, 45813246635, 183252462251, 733008800427, 2932033104555, 11728128223915, 46912504507051, 187650001250987
Offset: 0

Views

Author

Keywords

Comments

Number of palindromic structures using a maximum of four different symbols. - Marks R. Nester
Dimension of the universal embedding of the symplectic dual polar space DSp(2n,2) (conjectured by A. Brouwer, proved by P. Li). - J. Taylor (jt_cpp(AT)yahoo.com), Apr 02 2004.
Apart from initial term, same as A124303. - Valery A. Liskovets, Nov 16 2006
Hankel transform is := [1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, ...]. - Philippe Deléham, Dec 04 2008
a(n) is also the number of distinct solutions (avoiding permutations) to the equation: XOR(A,B,C)=0 where A,B,C are n-bit binary numbers. - Ramasamy Chandramouli, Jan 11 2009
The rank of the fundamental group of the Z_p^n - cobordism category in dimension 1+1 for the case p=2 (see paper below). The expression for any prime p is (p^(2n-1)+p^(n+1)-p^(n-1)+p^2-p-1)/(p^2-1). - Carlos Segovia Gonzalez, Dec 05 2012
The number of isomorphic classes of regular four coverings of a graph with respect to the identity automorphism (S. Hong and J. H. Kwak). - Carlos Segovia Gonzalez, Aug 01 2013
The density of a language with four letters (N. Moreira and R. Reis). - Carlos Segovia Gonzalez, Aug 01 2013

References

  • P. Li, On the Brouwer Conjecture for Dual Polar Spaces of Symplectic Type over GF(2). Preprint.
  • M. R. Nester (1999). Mathematical investigations of some plant interaction designs. PhD Thesis. University of Queensland, Brisbane, Australia. [See A056391 for pdf file of Chap. 2]
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

A row of the array in A278984.

Programs

Formula

a(n) = (3*2^(n-1) + 2^(2*n-1) + 1)/3.
a(n) = Sum_{k=1..4} Stirling2(n, k). - Winston Yang (winston(AT)cs.wisc.edu), Aug 23 2000
Binomial transform of 3^n/6 + 1/2 + 0^n/3, i.e., of A007051 with an extra leading 1. a(n) = binomial(2^n+2, 2^n-1)/2^n. - Paul Barry, Jul 19 2003
a(n) = C(2+2^n, 3)/2^n = a(n-1) + 2^(n-1) + 4^(n-3/2) = A092055(n)/A000079(n). - Henry Bottomley, Feb 19 2004
Second binomial transform of A001045(n-1) + 0^n/2. G.f.: (1-5*x+5*x^2)/((1-x)*(1-2*x)*(1-4*x)). - Paul Barry, Apr 28 2004
a(n) is the top entry of the vector M^n*[1,1,1,1,0,0,0,...], where M is an infinite bidiagonal matrix with M(r,r)=r, r >= 1, as the main diagonal, M(r,r+1)=1, and the rest zeros. ([1,1,1,...] is a column vector and transposing gives the same in terms of a leftmost column term.) - Gary W. Adamson, Jun 24 2011
a(0)=1, a(1)=2, a(2)=5, a(n) = 7*a(n-1) - 14*a(n-2) + 8*a(n-3). - Harvey P. Dale, Jul 24 2011
E.g.f.: (exp(2*x) + 1/3*exp(4*x) + 2/3*exp(x))/2 = G(0)/2; G(k)=1 + (2^k)/(3 - 6/(2 + 4^k - 3*x*(8^k)/(3*x*(2^k) + (k+1)/G(k+1)))); (continued fraction). - Sergei N. Gladkovskii, Dec 08 2011

A278984 Array read by antidiagonals downwards: T(b,n) = number of words of length n over an alphabet of size b that are in standard order.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 4, 2, 1, 1, 8, 5, 2, 1, 1, 16, 14, 5, 2, 1, 1, 32, 41, 15, 5, 2, 1, 1, 64, 122, 51, 15, 5, 2, 1, 1, 128, 365, 187, 52, 15, 5, 2, 1, 1, 256, 1094, 715, 202, 52, 15, 5, 2, 1, 1, 512, 3281, 2795, 855, 203, 52, 15, 5, 2, 1, 1, 1024, 9842, 11051, 3845, 876, 203, 52, 15, 5, 2, 1
Offset: 1

Views

Author

Joerg Arndt and N. J. A. Sloane, Dec 05 2016

Keywords

Comments

We study words made of letters from an alphabet of size b, where b >= 1. We assume the letters are labeled {1,2,3,...,b}. There are b^n possible words of length n.
We say that a word is in "standard order" if it has the property that whenever a letter i appears, the letter i-1 has already appeared in the word. This implies that all words begin with the letter 1.
Let X be the random variable that assigns to each permutation of {1,2,...,b} (with uniform distribution) its number of fixed points (as in A008290). Then T(b,n) is the n-th moment about 0 of X, i.e., the expected value of X^n. - Geoffrey Critzer, Jun 23 2020

Examples

			The array begins:
1,.1,..1,...1,...1,...1,...1,....1..; b=1, A000012
1,.2,..4,...8,..16,..32,..64,..128..; b=2, A000079
1,.2,..5,..14,..41,.122,.365,.1094..; b=3, A007051 (A278985)
1,.2,..5,..15,..51,.187,.715,.2795..; b=4, A007581
1,.2,..5,..15,..52,.202,.855,.3845..; b=5, A056272
1,.2,..5,..15,..52,.203,.876,.4111..; b=6, A056273
...
The rows tend to A000110.
		

Crossrefs

Rows 1 through 16 of the array are: A000012, A000079, A007051 (or A124302), A007581 (or A124303), A056272, A056273, A099262, A099263, A164863, A164864, A203641-A203646.
The limit of the rows is A000110, the Bell numbers.
See A278985 for the words arising in row b=3.
Cf. A203647, A137855 (essentially same table).

Programs

  • Maple
    with(combinat);
    f1:=proc(L,b) local t1;i;
    t1:=add(stirling2(L,i),i=1..b);
    end:
    Q1:=b->[seq(f1(L,b), L=1..20)]; # the rows of the array are Q1(1), Q1(2), Q1(3), ...
  • Mathematica
    T[b_, n_] := Sum[StirlingS2[n, j], {j, 1, b}]; Table[T[b-n+1, n], {b, 1, 12}, {n, b, 1, -1}] // Flatten (* Jean-François Alcover, Feb 18 2017 *)

Formula

The number of words of length n over an alphabet of size b that are in standard order is Sum_{j = 1..b} Stirling2(n,j).

A124302 Number of set partitions with at most 3 blocks; number of Dyck paths of height at most 4; dimension of space of symmetric polynomials in 3 noncommuting variables.

Original entry on oeis.org

1, 1, 2, 5, 14, 41, 122, 365, 1094, 3281, 9842, 29525, 88574, 265721, 797162, 2391485, 7174454, 21523361, 64570082, 193710245, 581130734, 1743392201, 5230176602, 15690529805, 47071589414, 141214768241, 423644304722, 1270932914165, 3812798742494, 11438396227481
Offset: 0

Views

Author

Mike Zabrocki, Oct 25 2006

Keywords

Comments

Row sums of triangle in A056241. - Philippe Deléham, Oct 30 2006
Row sums of triangle in A147746. - Philippe Deléham, Dec 04 2008
Hankel transform is := [1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, ...]. - Philippe Deléham, Dec 04 2008
Number of nonisomorphic graded posets with 0 and 1 and uniform Hasse graph of rank n with no 3-element antichain. (Uniform used in the sense of Retakh, Serconek and Wilson. Graded used in Stanley's sense that every maximal chain has the same length n.) - David Nacin, Feb 26 2012
Number of Dyck paths of length 2n and height at most 4. - Ira M. Gessel, Aug 06 2012

Examples

			There are 15 set partitions of {1,2,3,4}, only {{1},{2},{3},{4}} has more than 3 blocks, so a(4) = 14.
G.f. = 1 + x + 2*x^2 + 5*x^3 + 14*x^4 + 41*x^5 + 122*x^6 + 365*x^7 + ...
		

References

  • R. Stanley, Enumerative combinatorics, Vol. 1, Cambridge University Press, Cambridge, 1997, pp. 96-100.

Crossrefs

Essentially the same as A007051.

Programs

  • Magma
    I:=[1, 1, 2]; [n le 3 select I[n] else  4*Self(n-1) - 3*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Dec 25 2012
    
  • Maple
    a:= proc(n); if n<3 then [1,1,2][n+1]; else 4*a(n-1)-3*a(n-2); fi; end:
    # Mike Zabrocki, Oct 25 2006
    with(GraphTheory): G:=PathGraph(5): A:= AdjacencyMatrix(G): nmax:=27; for n from 0 to 2*nmax do B(n):=A^n; b(n):=B(n)[1,1]; od: for n from 0 to nmax do a(n):=b(2*n) od: seq(a(n),n=0..nmax);
    # Johannes W. Meijer, May 29 2010
  • Mathematica
    a=Exp[x]-1; Range[0, 20]! CoefficientList[Series[1+a+a^2/2+a^3/6, {x,0,20}],x]
    Join[{1}, LinearRecurrence[{4, -3}, {1, 2}, 20]] (* David Nacin, Feb 26 2012 *)
    CoefficientList[Series[1 / (1 - x / (1 - x / (1 - x / (1 - x)))), {x, 0, 30}], x] (* Vincenzo Librandi, Dec 25 2012 *)
    Table[Sum[StirlingS2[n,k],{k,0,3}],{n,0,30}] (* Robert A. Russell, Mar 29 2018 *)
  • PARI
    {a(n) = if( n<1, n==0, (3^(n-1) + 1) / 2)}; /* Michael Somos, Apr 03 2014 */
  • Python
    def a(n, adict={0:1, 1:1, 2:2}):
        if n in adict:
            return adict[n]
        adict[n]=4*a(n-1) - 3*a(n-2)
        return adict[n] # David Nacin, Mar 04 2012
    

Formula

O.g.f.: (q^2 - 3*q + 1)/(3*q^2 - 4*q + 1) = Sum_{k=0..3} (q^k/Product_{i=1..k} (1-i*q)).
a(n) = 4*a(n-1) - 3*a(n-2); a(0) = 1, a(1) = 1, a(2) = 2, a(n) = Sum_{k=1..3} A008277(n,k).
Inverse binomial transform of A007581. - Philippe Deléham, Oct 30 2006
a(n) = Sum_{k=0..n} A056241(n,k), n >= 1. - Philippe Deléham, Oct 30 2006
a(0) = 1, a(n) = (3^(n-1) + 1)/2 for n >= 1, see A007051. - Philippe Deléham, Oct 30 2006
E.g.f.: (2 + 3*exp(x) + exp(3x))/6.
G.f.: 1 / (1 - x / (1 - x / (1 - x / (1 - x)))). - Michael Somos, May 03 2012
G.f.: 1 + x + 3*x^2*U(0)/2 where U(k) = 1 + 2/(3*3^k + 3*3^k/(1 - 18*x*3^k/ (9*x*3^k - 1/U(k+1)))); (continued fraction, 4-step). - Sergei N. Gladkovskii, Nov 01 2012
G.f.: 1+x*G(0) where G(k) = 1 + 2*x/( 1-2*x - x*(1-2*x)/(x + (1-2*x)*2/G(k+1) )); (recursively defined continued fraction). - Sergei N. Gladkovskii, Dec 10 2012
a(n) = Sum_{k=0..3} Stirling2(n,k). - Robert A. Russell, Mar 29 2018
G.f.: Sum_{j=0..k} A248925(k,j)*x^j / Product_{j=1..k} 1-j*x with k=3. - Robert A. Russell, Apr 25 2018

A056273 Word structures of length n using a 6-ary alphabet.

Original entry on oeis.org

1, 1, 2, 5, 15, 52, 203, 876, 4111, 20648, 109299, 601492, 3403127, 19628064, 114700315, 676207628, 4010090463, 23874362200, 142508723651, 852124263684, 5101098232519, 30560194493456, 183176170057707, 1098318779272060, 6586964947803695, 39510014478620232, 237013033135668883
Offset: 0

Views

Author

Keywords

Comments

Set partitions of the n-set into at most 6 parts; also restricted growth strings (RGS) with six letters s(1),s(2),...,s(6) where the first occurrence of s(j) precedes the first occurrence of s(k) for all j < k. - Joerg Arndt, Jul 06 2011
Permuting the alphabet will not change a word structure. Thus aabc and bbca have the same structure.
Density of regular language L over {1,2,3,4,5,6}^* (i.e., number of strings of length n in L) described by regular expression with c=6: Sum_{i=1..c} Product_{j=1..i} (j(1+...+j)*) where Sum stands for union and Product for concatenation. - Nelma Moreira, Oct 10 2004
Word structures of length n using an N-ary alphabet are generated by taking M^n* the vector [(N 1's),0,0,0,...], leftmost column term = a(n+1). In the case of A056273, the vector = [1,1,1,1,1,1,0,0,0,...]. As the vector approaches all 1's, the leftmost column terms approach A000110, the Bell sequence. - Gary W. Adamson, Jun 23 2011
From Gary W. Adamson, Jul 06 2011: (Start)
Construct an infinite array of sequences representing word structures of length n using an N-ary alphabet as follows:
.
1, 1, 1, 1, 1, 1, 1, 1, ...; N=1, A000012
1, 2, 4, 8, 16, 32, 64, 128, ...; N=2, A000079
1, 2, 5, 14, 41, 122, 365, 1094, ...; N=3, A007051
1, 2, 5, 15, 51, 187, 715, 2795, ...; N=4, A007581
1, 2, 5, 15, 52, 202, 855, 3845, ...; N=5, A056272
1, 2, 5, 15, 52, 203, 876, 4111, ...; N=6, A056273
...
The sequences tend to A000110. Finite differences of columns reinterpreted as rows generate A008277 as a triangle: (1; 1,1; 1,3,1; 1,7,6,1; ...). (End)

Examples

			For a(4) = 15, the 7 achiral patterns are AAAA, AABB, ABAB, ABBA, ABBC, ABCA, and ABCD; the 8 chiral patterns are the 4 pairs AAAB-ABBB, AABA-ABAA, AABC-ABCC, and ABAC-ABCB.
		

References

  • M. R. Nester (1999). Mathematical investigations of some plant interaction designs. PhD Thesis. University of Queensland, Brisbane, Australia. [See A056391 for pdf file of Chap. 2]

Crossrefs

A row of the array in A278984 and A320955.
Cf. A056325 (unoriented), A320936 (chiral), A305752 (chiral).

Programs

  • GAP
    List([0..25],n->Sum([0..6],k->Stirling2(n,k))); # Muniru A Asiru, Oct 30 2018
    
  • Magma
    [(&+[StirlingSecond(n, i): i in [0..6]]): n in [0..30]]; // Vincenzo Librandi, Nov 07 2018
  • Maple
    egf := (265+264*exp(x)+135*exp(x*2)+40*exp(x*3)+15*exp(x*4)+exp(6*x))/6!:
    ser := series(egf,x,30): seq(n!*coeff(ser,x,n),n=0..22); # Peter Luschny, Nov 06 2018
  • Mathematica
    Table[Sum[StirlingS2[n,k],{k,0,6}],{n,0,30}] (* or *) LinearRecurrence[ {16,-95,260,-324,144},{1,1,2,5,15,52},30] (* Harvey P. Dale, Jun 05 2015 *)
  • PARI
    Vec((1 - 15*x + 81*x^2 - 192*x^3 + 189*x^4 - 53*x^5)/((1-x)*(1-2*x)*(1-3*x)*(1-4*x)*(1-6*x)) + O(x^30)) \\ Michel Marcus, Nov 07 2018
    

Formula

a(n) = Sum_{k=0..6} Stirling2(n, k).
For n > 0, a(n) = (1/6!)*(6^n + 15*4^n + 40*3^n + 135*2^n + 264). - Vladeta Jovovic, Aug 17 2003
From Nelma Moreira, Oct 10 2004: (Start)
For n > 0 and c = 6:
a(n) = (c^n)/c! + Sum_{k=0..c-2} ((k^n)/k!*(Sum_{j=2..c-k}(((-1)^j)/j!))).
a(n) = Sum_{k=1..c} (g(k, c)*k^n) where g(1, 1) = 1; g(1, c) = g(1, c-1) + ((-1)^(c-1))/(c-1)! if c>1. For 2 <= k <= c: g(k, c) = g(k-1, c-1)/k if c>1. (End)
G.f.: (1 - 15*x + 81*x^2 - 192*x^3 + 189*x^4 - 53*x^5)/((1-x)*(1-2x)*(1-3x)*(1-4x)*(1-6x)). - Maksym Voznyy (voznyy(AT)mail.ru), Jul 26 2009 [corrected by R. J. Mathar, Sep 16 2009] [Adapted to offset 0 by Robert A. Russell, Nov 06 2018]
G.f.: Sum_{j=0..k} A248925(k,j)*x^j / Product_{j=1..k} 1-j*x with k=6. - Robert A. Russell, Apr 25 2018
E.g.f.: (265 + 264*exp(x) + 135*exp(x*2) + 40*exp(x*3) + 15*exp(x*4) + exp(6*x))/6!. - Peter Luschny, Nov 06 2018

Extensions

a(0)=1 prepended by Robert A. Russell, Nov 06 2018

A198982 T(n,k) = number of n X k 0..5 arrays with values 0..5 introduced in row major order and no element equal to any horizontal or vertical neighbor.

Original entry on oeis.org

1, 1, 1, 2, 4, 2, 5, 34, 34, 5, 15, 481, 1835, 481, 15, 52, 8731, 146286, 146286, 8731, 52, 202, 174454, 12662226, 53082012, 12662226, 174454, 202, 855, 3603244, 1112962873, 19622872903, 19622872903, 1112962873, 3603244, 855, 3845, 75251971
Offset: 1

Views

Author

R. H. Hardin, Nov 01 2011

Keywords

Comments

Number of colorings of the grid graph P_n X P_k using a maximum of 6 colors up to permutation of the colors. - Andrew Howroyd, Jun 26 2017

Examples

			Table starts
.....1...........1.................2........................5
.....1...........4................34......................481
.....2..........34..............1835...................146286
.....5.........481............146286.................53082012
....15........8731..........12662226..............19622872903
....52......174454........1112962873............7267830860056
...202.....3603244.......98102456246.........2692353648978984
...855....75251971.....8651794282083.......997397244990907738
..3845..1577395861...763087851014929....369492074075459555844
.18002.33105096904.67305520316532514.136880688981914387733120
...
Some solutions with all values from 0 to 5 for n=6, k=4:
..0..1..0..1....0..1..0..1....0..1..0..1....0..1..0..1....0..1..0..1
..1..2..3..0....1..2..3..0....1..2..3..0....1..2..3..0....1..2..3..0
..0..3..0..1....0..3..0..1....0..3..0..1....0..3..0..1....0..3..0..1
..2..0..3..0....1..0..3..0....1..0..3..0....1..0..3..0....1..0..3..0
..1..2..0..2....4..1..0..1....3..1..0..4....3..4..2..1....3..4..2..5
..4..5..2..1....5..4..2..4....5..2..3..1....1..5..1..5....5..1..5..2
		

Crossrefs

Columns 1-7 are A056272(n-1), A198976, A198977, A198978, A198979, A198980, A198981.
Main diagonal is A198975.
Cf. A207997 (3 colorings), A198715 (4 colorings), A198906 (5 colorings), A222281 (labeled 6 colorings), A198723 (7 colorings), A198914 (8 colorings), A207868 (unlimited).

A056324 Number of reversible string structures with n beads using a maximum of five different colors.

Original entry on oeis.org

1, 1, 2, 4, 11, 32, 116, 455, 1993, 9134, 43580, 211659, 1041441, 5156642, 25640456, 127773475, 637624313, 3184387574, 15910947980, 79521737939, 397510726681, 1987259550002, 9935420646296, 49674470817195, 248364482308833, 1241798790172214
Offset: 0

Views

Author

Keywords

Comments

A string and its reverse are considered to be equivalent. Permuting the colors will not change the structure. Thus aabc, cbaa and bbac are all considered to be identical.
Number of set partitions of an unoriented row of n elements with five or fewer nonempty subsets. - Robert A. Russell, Oct 28 2018
There are nonrecursive formulas, generating functions, and computer programs for A056272 and A305751, which can be used in conjunction with the formula. - Robert A. Russell, Oct 28 2018
From Allan Bickle, Jun 02 2022: (Start)
a(n) is the number of (unlabeled) 5-paths with n+7 vertices. (A 5-path with order n at least 7 can be constructed from a 5-clique by iteratively adding a new 5-leaf (vertex of degree 5) adjacent to an existing 5-clique containing an existing 5-leaf.)
Recurrences appear in the papers by Bickle, Eckhoff, and Markenzon et al. (End)

Examples

			For a(4)=11, the 7 achiral patterns are AAAA, AABB, ABAB, ABBA, ABCA, ABBC, and ABCD.  The 4 chiral pairs are AAAB-ABBB, AABA-ABAA, AABC-ABCC, and ABAC-ABCB.
		

References

  • M. R. Nester (1999). Mathematical investigations of some plant interaction designs. PhD Thesis. University of Queensland, Brisbane, Australia. [See A056391 for pdf file of Chap. 2]

Crossrefs

Cf. A032122.
Column 5 of A320750.
Cf. A056272 (oriented), A320935 (chiral), A305751 (achiral).
The numbers of unlabeled k-paths for k = 2..7 are given in A005418, A001998, A056323, A056324, A056325, and A345207, respectively.
The sequences above converge to A103293(n+1).

Programs

  • Mathematica
    Ach[n_, k_] := Ach[n, k] = If[n<2, Boole[n==k && n>=0], k Ach[n-2,k] + Ach[n-2,k-1] + Ach[n-2,k-2]] (* A304972 *)
    k=5; Table[Sum[StirlingS2[n,j]+Ach[n,j],{j,0,k}]/2,{n,0,40}]  (* Robert A. Russell, Oct 28 2018 *)
    LinearRecurrence[{11, -34, -16, 247, -317, -200, 610, -300}, {1, 1, 2, 4, 11, 32, 116, 455, 1993}, 40] (* Robert A. Russell, Oct 28 2018 *)

Formula

Use de Bruijn's generalization of Polya's enumeration theorem as discussed in reference.
G.f.: (1-10x+25x^2+32x^3-196x^4+149x^5+225x^6-321x^7+85x^8)/((1-x)*(1-2x)*(1-3x)*(1-5x)*(1-2x^2)*(1-5x^2)). - Colin Barker, Nov 24 2012 [Adapted to offset 0 by Robert A. Russell, Nov 07 2018]
From Robert A. Russell, Oct 28 2018: (Start)
a(n) = (A056272(n) + A305751(n)) / 2.
a(n) = A056272(n) - A320935(n) = A320935(n) + A305751(n).
a(n) = Sum_{j=0..k} (S2(n,j) + Ach(n,j)) / 2, where k=5 is the maximum number of colors, S2 is the Stirling subset number A008277, and Ach(n,k) = [n>=0 & n<2 & n==k] + [n>1]*(k*Ach(n-2,k) + Ach(n-2,k-1) + Ach(n-2,k-2)).
a(n) = A000007(n) + A057427(n) + A056326(n) + A056327(n) + A056328(n) + A056329(n). (End)
For n>8, a(n) = 11*a(n-1) - 34*a(n-2) - 16*a(n-3) + 247*a(n-4) - 317*a(n-5) - 200*a(n-6) + 610*a(n-7) - 300*a(n-8). - Muniru A Asiru, Oct 30 2018
From Allan Bickle, Jun 04 2022: (Start)
a(n) = 5^n/240 + 3^n/24 + 2^n/12 + 13*5^(n/2)/120 + 2^(n/2)/6 + 5/16 for n>0 even;
a(n) = 5^n/240 + 3^n/24 + 2^n/12 + 5^((n+1)/2)/24 + 2^((n+1)/2)/12 + 5/16 for n>0 odd. (End)

Extensions

Terms added by Robert A. Russell, Oct 30 2018
a(0)=1 prepended by Robert A. Russell, Nov 07 2018

A164864 Number of ways of placing n labeled balls into 10 indistinguishable boxes; word structures of length n using a 10-ary alphabet.

Original entry on oeis.org

1, 1, 2, 5, 15, 52, 203, 877, 4140, 21147, 115975, 678569, 4213530, 27641927, 190829797, 1381367941, 10448276360, 82285618467, 672294831619, 5676711562593, 49344452550230, 439841775811967, 4005444732928641, 37136385907400125, 349459367068932740
Offset: 0

Views

Author

Alois P. Heinz, Aug 28 2009

Keywords

Crossrefs

Programs

  • Maple
    # First program:
    a:= n-> ceil(2119/11520*2^n +103/1680*3^n +53/3456*4^n +11/3600*5^n +6^n/1920 +7^n/15120 +8^n/80640 +10^n/3628800): seq(a(n), n=0..25);
    # second program:
    a:= n-> add(Stirling2(n, k), k=0..10): seq(a(n), n=0..25);
  • Mathematica
    Table[Sum[StirlingS2[n,k],{k,0,10}],{n,0,30}] (* Harvey P. Dale, Nov 22 2023 *)

Formula

a(n) = Sum_{k=0..10} Stirling2 (n,k).
a(n) = ceiling(2119/11520*2^n +103/1680*3^n +53/3456*4^n +11/3600*5^n +6^n/1920 +7^n/15120 +8^n/80640 +10^n/3628800).
G.f.: (148329*x^9 -613453*x^8 +855652*x^7 -596229*x^6 +240065*x^5 -59410*x^4 +9177*x^3 -862*x^2 +45*x-1) / ((10*x-1) *(8*x-1) *(7*x-1) *(6*x-1) *(5*x-1) *(4*x-1) *(3*x-1) *(2*x-1) *(x-1)).
a(n) <= A000110(n) with equality only for n <= 10.

A198528 T(n,k) is the number of n X k 0..4 arrays with values 0..4 introduced in row major order and each element equal to no more than two horizontal and vertical neighbors.

Original entry on oeis.org

1, 2, 2, 5, 15, 5, 15, 193, 193, 15, 52, 3660, 16714, 3660, 52, 202, 81844, 1877316, 1877316, 81844, 202, 855, 1948672, 222590953, 1084539825, 222590953, 1948672, 855, 3845, 47436498, 26670041125, 634586561196, 634586561196, 26670041125
Offset: 1

Views

Author

R. H. Hardin, Oct 26 2011

Keywords

Examples

			Table starts:
.....1............2...................5......................15
.....2...........15.................193....................3660
.....5..........193...............16714.................1877316
....15.........3660.............1877316..............1084539825
....52........81844...........222590953............634586561196
...202......1948672.........26670041125.........371815743708461
...855.....47436498.......3201911378187......217885196778717007
..3845...1163606279.....384557171168810...127683385189755792564
.18002..28617909415...46189600128813487.74824145653256981522691
.86472.704465305625.5547962760669962158
Some solutions for n=6 and k=4:
..0..0..0..0....0..0..0..0....0..0..0..0....0..0..0..0....0..0..0..0
..0..1..1..0....0..1..1..0....0..1..1..0....0..1..1..0....0..1..1..0
..0..0..0..0....0..0..0..0....0..0..0..0....0..0..0..0....0..0..0..0
..2..2..1..2....2..2..1..2....2..2..1..2....2..2..1..2....2..2..1..2
..2..2..3..2....0..3..1..2....3..0..0..3....2..0..1..0....0..3..4..4
..1..4..0..4....3..2..4..4....3..1..4..3....1..3..0..0....3..3..0..4
		

Crossrefs

Main diagonal is A198521.
Showing 1-10 of 23 results. Next