Search: a096231 -id:a096231
|
|
A000931
|
|
Padovan sequence (or Padovan numbers): a(n) = a(n-2) + a(n-3) with a(0) = 1, a(1) = a(2) = 0.
(Formerly M0284 N0102)
|
|
+10
246
|
|
|
1, 0, 0, 1, 0, 1, 1, 1, 2, 2, 3, 4, 5, 7, 9, 12, 16, 21, 28, 37, 49, 65, 86, 114, 151, 200, 265, 351, 465, 616, 816, 1081, 1432, 1897, 2513, 3329, 4410, 5842, 7739, 10252, 13581, 17991, 23833, 31572, 41824, 55405, 73396, 97229, 128801, 170625
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,9
|
|
COMMENTS
|
Number of compositions of n into parts congruent to 2 mod 3 (offset -1). - Vladeta Jovovic, Feb 09 2005
a(n) is the number of compositions of n into parts that are odd and >= 3. Example: a(10)=3 counts 3+7, 5+5, 7+3. - David Callan, Jul 14 2006
Referred to as N0102 in R. K. Guy's "Anyone for Twopins?" - Rainer Rosenthal, Dec 05 2006
Zagier conjectures that a(n+3) is the maximum number of multiple zeta values of weight n > 1 which are linearly independent over the rationals. - Jonathan Sondow and Sergey Zlobin (sirg_zlobin(AT)mail.ru), Dec 20 2006
Starting with offset 6: (1, 1, 2, 2, 3, 4, 5, ...) = INVERT transform of A106510: (1, 1, -1, 0, 1, -1, 0, 1, -1, ...). - Gary W. Adamson, Oct 10 2008
Starting with offset 7, the sequence 1, 2, 2, 3, 4, 5, 7, 9, 12, 16, 21, 28, ... is called the Fibonacci quilt sequence by Catral et al., in Fib. Q. 2017. - N. J. A. Sloane, Dec 24 2021
Triangle A145462: right border = A000931 starting with offset 6. Row sums = Padovan sequence starting with offset 7. - Gary W. Adamson, Oct 10 2008
Starting with offset 3 = row sums of triangle A146973 and INVERT transform of [1, -1, 2, -2, 3, -3, ...]. - Gary W. Adamson, Nov 03 2008
a(n+5) corresponds to the diagonal sums of "triangle": 1; 1; 1,1; 1,1; 1,2,1; 1,2,1; 1,3,3,1; 1,3,3,1; 1,4,6,4,1; ..., rows of Pascal's triangle (A007318) repeated. - Philippe Deléham, Dec 12 2008
With offset 3: (1, 0, 1, 1, 1, 2, 2, ...) convolved with the tribonacci numbers prefaced with a "1": (1, 1, 1, 2, 4, 7, 13, ...) = the tribonacci numbers, A000073. (Cf. triangle A153462.) - Gary W. Adamson, Dec 27 2008
a(n) is also the number of strings of length (n-8) from an alphabet {A, B} with no more than one A or 2 B's consecutively. (E.g., n = 4: {ABAB,ABBA,BABA,BABB,BBAB} and a(4+8) = 5.) - Toby Gottfried, Mar 02 2010
p(n):=A000931(n+3), n >= 1, is the number of partitions of the numbers {1,2,3,...,n} into lists of length two or three containing neighboring numbers. The 'or' is inclusive. For n=0 one takes p(0)=1. For details see the W. Lang link. There the explicit formula for p(n) (analog of the Binet-de Moivre formula for Fibonacci numbers) is also given. Padovan sequences with different inputs are also considered there. - Wolfdieter Lang, Jun 15 2010
Equals the INVERTi transform of Fibonacci numbers prefaced with three 1's, i.e., (1 + x + x^2 + x^3 + x^4 + 2x^5 + 3x^6 + 5x^7 + 8x^8 + 13x^9 + ...). - Gary W. Adamson, Apr 01 2011
When run backwards gives (-1)^n*A050935(n).
a(n) is the top left entry of the n-th power of the 3 X 3 matrix [0, 0, 1; 1, 0, 1; 0, 1, 0] or of the 3 X 3 matrix [0, 1, 0; 0, 0, 1; 1, 1, 0]. - R. J. Mathar, Feb 03 2014
Figure 4 of Brauchart et al., 2014, shows a way to "visualize the Padovan sequence as cuboid spirals, where the dimensions of each cuboid made up by the previous ones are given by three consecutive numbers in the sequence". - N. J. A. Sloane, Mar 26 2014
a(n) is the number of closed walks from a vertex of a unidirectional triangle containing an opposing directed edge (arc) between the second and third vertices. Equivalently the (1,1) entry of A^n where the adjacency matrix of digraph is A=(0,1,0;0,0,1;1,1,0). - David Neil McGrath, Dec 19 2014
Number of compositions of n-3 (n >= 4) into 2's and 3's. Example: a(12)=5 because we have 333, 3222, 2322, 2232, and 2223. - Emeric Deutsch, Dec 28 2014
The Hoffman (2015) paper "offers significant evidence that the number of quantities needed to generate the weight-n multiple harmonic sums mod p is" a(n). - N. J. A. Sloane, Jun 24 2016
a(n) gives the number of compositions of n-5 into odd parts where the order of the 1's does not matter. For example, a(11)=4 counts the following compositions of 6: (5,1)=(1,5), (3,3), (3,1,1,1)=(1,3,1,1)=(1,1,3,1)=(1,1,1,3), (1,1,1,1,1,1). - Gregory L. Simay, Aug 04 2016
For n > 6, a(n) is the number of maximal matchings in the (n-5)-path graph, maximal independent vertex sets and minimal vertex covers in the (n-6)-path graph, and minimal edge covers in the (n-5)-pan graph and (n-3)-path graphs. - Eric W. Weisstein, Mar 30, Aug 03, and Aug 07 2017
a(2n + 5) + 2n - 4, n > 2, is the number of maximal subsemigroups of the monoid of order-preserving mappings on a set with n elements.
a(n + 6) + n - 3, n > 3, is the number of maximal subsemigroups of the monoid of order-preserving or reversing mappings on a set with n elements.
(End)
Has the property that the largest of any four consecutive terms equals the sum of the two smallest. - N. J. A. Sloane, Aug 29 2017 [David Nacin points out that there are many sequences with this property, such as 1,1,1,2,1,1,1,2,1,1,1,2,... or 2,3,4,5,2,3,4,5,2,3,4,5,... or 2,2,1,3,3, 4,1,4, 5,5,1,6,6, 7,1,7, 8,8,1,9,9, 10,1,10, ... (spaces added for clarity), and a conjecture I made here in 2017 was simply wrong. I have deleted it. - N. J. A. Sloane, Oct 23 2018]
a(n) is also the number of maximal cliques in the (n+6)-path complement graph. - Eric W. Weisstein, Apr 12 2018
a(n+8) is the number of solus bitstrings of length n with no runs of 3 zeros. - Steven Finch, Mar 25 2020
Named after the architect Richard Padovan (b. 1935). - Amiram Eldar, Jun 08 2021
Shannon et al. (2006) credit a French architecture student Gérard Cordonnier with the discovery of these numbers.
For n >= 3, a(n) is the number of sequences of 0s and 1s of length (n-2) that begin with a 0, end with a 0, contain no two consecutive 0s, and contain no three consecutive 1s. - Yifan Xie, Oct 20 2022
|
|
REFERENCES
|
A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, p. 47, ex. 4.
Minerva Catral, Pari L. Ford, Pamela E. Harris, Steven J. Miller, Dawn Nelson, Zhao Pan, and Huanzhong Xu, Legal Decompositions Arising from Non-positive Linear Recurrences, Fib. Quart., 55:3 (2017), 252-275. [Note that there is an earlier version of this paper, with only five authors, on the arXiv in 2016. Note to editors: do not merge these two citations. - N. J. A. Sloane, Dec 24 2021]
Richard K. Guy, "Anyone for Twopins?" in D. A. Klarner, editor, The Mathematical Gardner. Prindle, Weber and Schmidt, Boston, 1981, pp. 10-11.
Silvia Heubach and Toufik Mansour, Combinatorics of Compositions and Words, CRC Press, 2010.
A. G. Shannon, P. G. Anderson and A. F. Horadam, Properties of Cordonnier, Perrin and Van der Laan numbers, International Journal of Mathematical Education in Science and Technology, Volume 37:7 (2006), 825-831. See P_n.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
Ian Stewart, L'univers des nombres, "La sculpture et les nombres", pp. 19-20, Belin-Pour La Science, Paris, 2000.
Steven J. Tedford, Combinatorial identities for the Padovan numbers, Fib. Q., Vol. 57, No. 4 (2019), pp. 291-298.
Hans van der Laan, Het plastische getal. XV lessen over de grondslagen van de architectonische ordonnantie. Leiden, E.J. Brill, 1967.
Don Zagier, Values of zeta functions and their applications, in First European Congress of Mathematics (Paris, 1992), Vol. II, A. Joseph et al. (eds.), Birkhäuser, Basel, 1994, pp. 497-512.
|
|
LINKS
|
David Applegate, Marc LeBrun and N. J. A. Sloane, Dismal Arithmetic, J. Int. Seq. 14 (2011) # 11.9.8.
Richard K. Guy, Anyone for Twopins?, in D. A. Klarner, editor, The Mathematical Gardner. Prindle, Weber and Schmidt, Boston, 1981, pp. 2-15. [Annotated scanned copy, with permission]
Dov Jarden, Recurring Sequences, Riveon Lematematika, Jerusalem, 1966. [Annotated scanned copy] See p. 90.
Steven J. Miller and Alexandra Newlon, The Fibonacci Quilt Game, arXiv preprint arXiv:1909.01938 [math.NT], 2019. Also Fib. Q., Vol. 58, No. 2 (2020), pp. 157-168. (See Fig. 2, The "Fibonacci Quilt" sequence.)
Richard Padovan, Dom Hans van der Laan and the Plastic Number, Chapter 74, pp. 407-419, Volume II of K. Williams and M. J. Ostwald (eds.), Architecture and Mathematics from Antiquity to the Future, DOI 10.1007/978-3-319-00143-2_27, Springer International Publishing Switzerland, 2015.
Eric Weisstein's World of Mathematics, Pan Graph.
Eric Weisstein's World of Mathematics, Path Graph.
|
|
FORMULA
|
G.f.: (1-x^2)/(1-x^2-x^3).
a(n) is asymptotic to r^n / (2*r+3) where r = 1.3247179572447... = A060006, the real root of x^3 = x + 1. - Philippe Deléham, Jan 13 2004
a(n)^2 + a(n+2)^2 + a(n+6)^2 = a(n+1)^2 + a(n+3)^2 + a(n+4)^2 + a(n+5)^2 (Barniville, Question 16884, Ed. Times 1911).
a(n+5) = a(0) + a(1) + ... + a(n).
a(n) = central and lower right terms in the (n-3)-th power of the 3 X 3 matrix M = [0 1 0 / 0 0 1 / 1 1 0]. E.g., a(13) = 7. M^10 = [3 5 4 / 4 7 5 / 5 9 7]. - Gary W. Adamson, Feb 01 2004
G.f.: 1/(1 - x^3 - x^5 - x^7 - x^9 - ...). - Jon Perry, Jul 04 2004
a(n+4) = Sum_{k=0..floor((n-1)/2)} binomial(floor((n+k-2)/3), k). - Paul Barry, Jul 06 2004
a(n+3) is diagonal sum of A026729 (as a number triangle), with formula a(n+3) = Sum_{k=0..floor(n/2)} Sum_{i=0..n-k} (-1)^(n-k+i)*binomial(n-k, i)*binomial(i+k, i-k). - Paul Barry, Sep 23 2004
a(n+3) = Sum_{k=0..floor(n/2)} binomial((n-k)/2, k)(1+(-1)^(n-k))/2. - Paul Barry, Sep 09 2005
The sequence 1/(1-x^2-x^3) (a(n+3)) is given by the diagonal sums of the Riordan array (1/(1-x^3), x/(1-x^3)). The row sums are A000930. - Paul Barry, Feb 25 2005
a(n+5) corresponds to the diagonal sums of A030528. The binomial transform of a(n+5) is A052921. a(n+5) = Sum_{k=0..floor(n/2)} Sum_{k=0..n} (-1)^(n-k+i)*binomial(n-k, i)binomial(i+k+1, 2k+1). - Paul Barry, Jun 21 2004
r^(n-1) = (1/r)*a(n) + r*(n+1) + a(n+2), where r = 1.32471... is the real root of x^3 - x - 1 = 0. Example: r^8 = (1/r)*a(9) + r*a(10) + a(11) = (1/r)*2 + r*3 + 4 = 9.483909... - Gary W. Adamson, Oct 22 2006
a(n) = (r^n)/(2r+3) + (s^n)/(2s+3) + (t^n)/(2t+3) where r, s, t are the three roots of x^3-x-1. - Keith Schneider (schneidk(AT)email.unc.edu), Sep 07 2007
a(n) = -k*a(n-1) + a(n-2) + (k+1)a(n-2) + k*a(n-4), n > 3, for any value of k. - Gary Detlefs, Sep 13 2010
a(0) + a(2) + a(4) + a(6) + ... + a(2*n) = a(2*n+3).
a(0) + a(3) + a(6) + a(9) + ... + a(3*n) = a(3*n+2)+1.
a(0) + a(5) + a(10) + a(15) + ... + a(5*n) = a(5*n+1)+1.
a(0) + a(7) + a(14) + a(21) + ... + a(7*n) = (a(7*n) + a(7*n+1) + 1)/2. (End)
a(n+3) = Sum_{k=0..floor((n+1)/2)} binomial((n+k)/3,k), where binomial((n+k)/3,k)=0 for noninteger (n+k)/3. - Nikita Gogin, Dec 07 2012
a(n) = the k-th difference of a(n+5k) - a(n+5k-1), k>=1. For example, a(10)=3 => a(15)-a(14) => 2nd difference of a(20)-a(19) => 3rd difference of a(25)-a(24)... - Bob Selcoe, Mar 18 2014
Construct the power matrix T(n,j) = [A^*j]*[S^*(j-1)] where A=(0,0,1,0,1,0,1,...) and S=(0,1,0,0,...) or A063524. [* is convolution operation] Define S^*0=I with I=(1,0,0,...). Then a(n) = Sum_{j=1...n} T(n,j). - David Neil McGrath, Dec 19 2014
If x=a(n), y=a(n+1), z=a(n+2), then x^3 + 2*y*x^2 - z^2*x - 3*y*z*x + y^2*x + y^3 - y^2*z + z^3 = 1. - Alexander Samokrutov, Jul 20 2015
For the sequence shifted by 6 terms, a(n) = Sum_{k=ceiling(n/3)..ceiling(n/2)} binomial(k+1,3*k-n) [Doslic-Zubac]. - N. J. A. Sloane, Apr 23 2017
a(2n) = 2*a(n-1)*a(n) + a(n)^2 + a(n+1)^2, for n > 8.
a(2n-1) = 2*a(n)*a(n+1) + a(n-1)^2, for n > 8.
a(2n+1) = 2*a(n+1)*a(n+2) + a(n)^2, for n > 7. (End)
0*a(0) + 1*a(1) + 2*a(2) + ... + n*a(n) = n*a(n+5) - a(n+9) + 2. - Greg Dresden and Zi Ye, Jul 02 2021
2*a(n) = a(n+2) + a(n-5) for n >= 5.
3*a(n) = a(n+4) - a(n-9) for n >= 9.
4*a(n) = a(n+5) - a(n-9) for n >= 9. (End)
|
|
EXAMPLE
|
G.f. = 1 + x^3 + x^5 + x^6 + x^7 + 2*x^8 + 2*x^9 + 3*x^10 + 4*x^11 + ...
|
|
MAPLE
|
A000931 := proc(n) option remember; if n = 0 then 1 elif n <= 2 then 0 else procname(n-2)+procname(n-3); fi; end;
A000931:=-(1+z)/(-1+z^2+z^3); # Simon Plouffe in his 1992 dissertation; gives sequence without five leading terms
a[0]:=1; a[1]:=0; a[2]:=0; for n from 3 to 50 do a[n]:=a[n-2]+a[n-3]; end do; # Francesco Daddi, Aug 04 2011
|
|
MATHEMATICA
|
CoefficientList[Series[(1-x^2)/(1-x^2-x^3), {x, 0, 50}], x]
a[0]=1; a[1]=a[2]=0; a[n_]:= a[n]= a[n-2] + a[n-3]; Table[a[n], {n, 0, 50}] (* Robert G. Wilson v, May 04 2006 *)
LinearRecurrence[{0, 1, 1}, {1, 0, 0}, 50] (* Harvey P. Dale, Jan 10 2012 *)
Table[RootSum[-1 -# +#^3 &, 5#^n -6#^(n+1) +4#^(n+2) &]/23, {n, 0, 50}] (* Eric W. Weisstein, Nov 09 2017 *)
|
|
PROG
|
(Haskell)
a000931 n = a000931_list !! n
a000931_list = 1 : 0 : 0 : zipWith (+) a000931_list (tail a000931_list)
(PARI) {a(n) = if( n<0, polcoeff(1/(1+x-x^3) + x * O(x^-n), -n), polcoeff( (1 - x^2)/(1-x^2-x^3) + x * O(x^n), n))}; /* Michael Somos, Sep 18 2012 */
(Magma) I:=[1, 0, 0]; [n le 3 select I[n] else Self(n-2) + Self(n-3): n in [1..60]]; // Vincenzo Librandi, Jul 21 2015
(Sage)
P.<x> = PowerSeriesRing(ZZ, prec)
return P( (1-x^2)/(1-x^2-x^3) ).list()
(GAP) a:=[1, 0, 0];; for n in [4..50] do a[n]:=a[n-2]+a[n-3]; od; a; # G. C. Greubel, Dec 30 2019
(Python)
def aupton(nn):
alst = [1, 0, 0]
for n in range(3, nn+1): alst.append(alst[n-2]+alst[n-3])
return alst
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy,nice
|
|
AUTHOR
|
|
|
EXTENSIONS
|
Deleted certain dangerous or potentially dangerous links. - N. J. A. Sloane, Jan 30 2021
|
|
STATUS
|
approved
|
|
|
|
|
A007283
|
|
a(n) = 3*2^n.
(Formerly M2561)
|
|
+10
239
|
|
|
3, 6, 12, 24, 48, 96, 192, 384, 768, 1536, 3072, 6144, 12288, 24576, 49152, 98304, 196608, 393216, 786432, 1572864, 3145728, 6291456, 12582912, 25165824, 50331648, 100663296, 201326592, 402653184, 805306368, 1610612736, 3221225472, 6442450944, 12884901888
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,1
|
|
COMMENTS
|
Same as Pisot sequences E(3, 6), L(3, 6), P(3, 6), T(3, 6). See A008776 for definitions of Pisot sequences.
Also least number m such that 2^n is the smallest proper divisor of m which is also a suffix of m in binary representation, see A080940. - Reinhard Zumkeller, Feb 25 2003
Length of the period of the sequence Fibonacci(k) (mod 2^(n+1)). - Benoit Cloitre, Mar 12 2003
Total number of Latin n-dimensional hypercubes (Latin polyhedra) of order 3. - Kenji Ohkuma (k-ookuma(AT)ipa.go.jp), Jan 10 2007
Number of different ternary hypercubes of dimension n. - Edwin Soedarmadji (edwin(AT)systems.caltech.edu), Dec 10 2005
For n >= 1, a(n) is equal to the number of functions f:{1, 2, ..., n + 1} -> {1, 2, 3} such that for fixed, different x_1, x_2,...,x_n in {1, 2, ..., n + 1} and fixed y_1, y_2,...,y_n in {1, 2, 3} we have f(x_i) <> y_i, (i = 1,2,...,n). - Milan Janjic, May 10 2007
a(n) written in base 2: 11, 110, 11000, 110000, ..., i.e.: 2 times 1, n times 0 (see A003953). - Jaroslav Krizek, Aug 17 2009
Numbers containing the number 3 in their Collatz trajectories. - Reinhard Zumkeller, Feb 20 2012
a(n-1) gives the number of ternary numbers with n digits with no two adjacent digits in common; e.g., for n=3 we have 010, 012, 020, 021, 101, 102, 120, 121, 201, 202, 210 and 212. - Jon Perry, Oct 10 2012
If n > 1, then a(n) is a solution for the equation sigma(x) + phi(x) = 3x-4. This equation also has solutions 84, 3348, 1450092, ... which are not of the form 3*2^n. - Farideh Firoozbakht, Nov 30 2013
a(n) is the upper bound for the "X-ray number" of any convex body in E^(n + 2), conjectured by Bezdek and Zamfirescu, and proved in the plane E^2 (see the paper by Bezdek and Zamfirescu). - L. Edson Jeffery, Jan 11 2014
If T is a topology on a set V of size n and T is not the discrete topology, then T has at most 3 * 2^(n-2) many open sets. See Brown and Stephen references. - Ross La Haye, Jan 19 2014
Comment from Charles Fefferman, courtesy of Doron Zeilberger, Dec 02 2014: (Start)
Fix a dimension n. For a real-valued function f defined on a finite set E in R^n, let Norm(f, E) denote the inf of the C^2 norms of all functions F on R^n that agree with f on E. Then there exist constants k and C depending only on the dimension n such that Norm(f, E) <= C*max{ Norm(f, S) }, where the max is taken over all k-point subsets S in E. Moreover, the best possible k is 3 * 2^(n-1).
The analogous result, with the same k, holds when the C^2 norm is replaced, e.g., by the C^1, alpha norm (0 < alpha <= 1). However, the optimal analogous k, e.g., for the C^3 norm is unknown.
For the above results, see Y. Brudnyi and P. Shvartsman (1994). (End)
Also, coordination sequence for (infinity, infinity, infinity) tiling of hyperbolic plane. - N. J. A. Sloane, Dec 29 2015
The average of consecutive powers of 2 beginning with 2^1. - Melvin Peralta and Miriam Ong Ante, May 14 2016
For n > 1, a(n) is the smallest Zumkeller number with n divisors that are also Zumkeller numbers (A083207). - Ivan N. Ianakiev, Dec 09 2016
Also, for n >= 2, the number of length-n strings over the alphabet {0,1,2,3} having only the single letters as nonempty palindromic subwords. (Corollary 21 in Fleischer and Shallit) - Jeffrey Shallit, Dec 02 2019
Also, a(n) is the minimum link-length of any covering trail, circuit, path, and cycle for the set of the 2^(n+2) vertices of an (n+2)-dimensional hypercube. - Marco Ripà, Aug 22 2022
The finite subsequence a(3), a(4), a(5), a(6) = 24, 48, 96, 192 is one of only two geometric sequences that can be formed with all interior angles (all integer, in degrees) of a simple polygon. The other sequence is a subsequence of A000244 (see comment there). - Felix Huber, Feb 15 2024
A level 1 Sierpiński triangle is a triangle. Level n+1 is formed from three copies of level n by identifying pairs of corner vertices of each pair of triangles. For n>2, a(n-3) is the radius of the level n Sierpiński triangle graph. - Allan Bickle, Aug 03 2024
|
|
REFERENCES
|
Jason I. Brown, Discrete Structures and Their Interactions, CRC Press, 2013, p. 71.
T. Ito, Method, equipment, program and storage media for producing tables, Publication number JP2004-272104A, Japan Patent Office (written in Japanese, a(2)=12, a(3)=24, a(4)=48, a(5)=96, a(6)=192, a(7)=384 (a(7)=284 was corrected)).
Kenji Ohkuma, Atsuhiro Yamagishi and Toru Ito, Cryptography Research Group Technical report, IT Security Center, Information-Technology Promotion Agency, JAPAN.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
|
|
LINKS
|
|
|
FORMULA
|
G.f.: 3/(1-2*x).
a(n) = 2*a(n - 1), n > 0; a(0) = 3.
a(n) = Sum_{k = 0..n} (-1)^(k reduced (mod 3))*binomial(n, k). - Benoit Cloitre, Aug 20 2002
a(n) = abs(b(n) - b(n+3)) with b(n) = (-1)^n*A084247(n). (End)
|
|
MAPLE
|
|
|
MATHEMATICA
|
|
|
PROG
|
(PARI) a(n)=3*2^n
(Haskell)
a007283 = (* 3) . (2 ^)
a007283_list = iterate (* 2) 3
(Scala) (List.fill(40)(2: BigInt)).scanLeft(1: BigInt)(_ * _).map(3 * _) // Alonso del Arte, Nov 28 2019
(Python)
|
|
CROSSREFS
|
Coordination sequences for triangular tilings of hyperbolic space: A001630, A007283, A054886, A078042, A096231, A163876, A179070, A265057, A265058, A265059, A265060, A265061, A265062, A265063, A265064, A265065, A265066, A265067, A265068, A265069, A265070, A265071, A265072, A265073, A265074, A265075, A265076, A265077.
Subsequence of the following sequences: A029744, A029747, A029748, A029750, A362804 (after 3), A364494, A364496, A364289, A364291, A364292, A364295, A364497, A364964, A365422.
Row sums of (5, 1)-Pascal triangle A093562 and of (1, 5) Pascal triangle A096940.
|
|
KEYWORD
|
easy,nonn,changed
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|
|
A001630
|
|
Tetranacci numbers: a(n) = a(n-1) + a(n-2) + a(n-3) + a(n-4), with a(0)=a(1)=0, a(2)=1, a(3)=2.
(Formerly M0795 N0301)
|
|
+10
49
|
|
|
0, 0, 1, 2, 3, 6, 12, 23, 44, 85, 164, 316, 609, 1174, 2263, 4362, 8408, 16207, 31240, 60217, 116072, 223736, 431265, 831290, 1602363, 3088654, 5953572, 11475879, 22120468, 42638573, 82188492, 158423412, 305370945, 588621422, 1134604271, 2187020050
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,4
|
|
COMMENTS
|
Also (with a different offset), coordination sequence for (4,infinity,infinity) tiling of hyperbolic plane. - N. J. A. Sloane, Dec 29 2015
Apparently for n>=2 the number of 1-D walks of length n-2 using steps +1, +3 and -1, avoiding consecutive -1 steps. - David Scambler, Jul 15 2013
For n > 1, a(n) is the number of ways to tile a skew double-strip of n-1 cells with one extra initial cell, using squares and all possible "dominoes". Here is the skew double-strip corresponding to n=12, with 11 cells:
___ ___ ___ ___ ___
| | | | | |
___ _|___|___|___|___|___|
| | | | | | |
|___|___|___|___|___|___|,
and here are the three possible "domino" tiles:
___ ___
| | | |
_| _| |_ |_ _______
| | | | | |
|___|, |___|, |_______|.
As an example, here is one of the a(12) = 609 ways to tile the skew double-strip of 11 cells:
___ _______ _______
| | | | | |
_____|_ |___|_ _|_ _| _|
| | | | | |
|_______|___|___|___|___|. (End)
|
|
REFERENCES
|
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
|
|
LINKS
|
|
|
FORMULA
|
G.f.: -x^2*(1+x)/(-1+x+x^2+x^3+x^4). [Simon Plouffe in his 1992 dissertation]
a(n) = 2*a(n-1) - a(n-5) with n>4, a(0)=a(1)=0, a(2)=1, a(3)=2, a(4)=3. - Vincenzo Librandi, Dec 21 2010
|
|
EXAMPLE
|
G.f. = x^2 + 2*x^3 + 3*x^4 + 6*x^5 + 12*x^6 + 23*x^7 + 44*x^8 + 85*x^9 + ...
|
|
MAPLE
|
a:= proc(n) option operator; local M; M := Matrix(4, (i, j)-> if (i=j-1) or j=1 then 1 else 0 fi)^n; M[1, 4]+M[1, 3] end; seq (a(n), n=0..34); # Alois P. Heinz, Aug 01 2008
|
|
MATHEMATICA
|
RecurrenceTable[{a[0] == a[1] == 0, a[2] == 1, a[3] == 2, a[n] == a[n - 1] + a[n - 2] + a[n - 3] + a[n - 4]}, a, {n, 35}] (* or *) a = {0, 0, 1, 2}; Do[AppendTo[a, a[[-1]] + a[[-2]] + a[[-3]] + a[[-4]]], {35}]; a (* Bruno Berselli, Jan 29 2013 *)
CoefficientList[Series[- x^2 * (1 + x)/(- 1 + x + x^2 + x^3 + x^4), {x, 0, 35}], x] (* Vincenzo Librandi, Jan 29 2013 *)
LinearRecurrence[{1, 1, 1, 1}, {0, 0, 1, 2}, 40] (* Harvey P. Dale, Aug 25 2013 *)
|
|
PROG
|
(Magma) I:=[0, 0, 1, 2]; [n le 4 select I[n] else Self(n-1)+ Self(n-2) + Self(n-3) + Self(n-4): n in [1..40]]; // Vincenzo Librandi, Jan 29 2013
(PARI) concat([0, 0], Vec(-x^2*(1+x)/(-1+x+x^2+x^3+x^4) + O(x^50))) \\ Michel Marcus, Dec 30 2015
|
|
CROSSREFS
|
Coordination sequences for triangular tilings of hyperbolic space: A001630, A007283, A054886, A078042, A096231, A163876, A179070, A265057, A265058, A265059, A265060, A265061, A265062, A265063, A265064, A265065, A265066, A265067, A265068, A265069, A265070, A265071, A265072, A265073, A265074, A265075, A265076, A265077.
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|
|
A179070
|
|
a(1)=a(2)=a(3)=1, a(4)=3; thereafter a(n) = a(n-1) + a(n-3).
|
|
+10
43
|
|
|
1, 1, 1, 3, 4, 5, 8, 12, 17, 25, 37, 54, 79, 116, 170, 249, 365, 535, 784, 1149, 1684, 2468, 3617, 5301, 7769, 11386, 16687, 24456, 35842, 52529, 76985, 112827, 165356, 242341, 355168, 520524, 762865, 1118033, 1638557, 2401422, 3519455, 5158012, 7559434
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,4
|
|
COMMENTS
|
Also (essentially), coordination sequence for (2,4,infinity) tiling of hyperbolic plane. - N. J. A. Sloane, Dec 29 2015
Column sums of shifted (1,2) Pascal array:
1 1 1 1 1 1 1 1 1
......2 3 4 5 6 7
............2 5 9
.................
----------------- +
1 1 1 3 4 5 8 ...
a(n+1) is the number of multus bitstrings of length n with no runs of 2 0's. - Steven Finch, Mar 25 2020
For n >= 5, a(n) gives the number of ways to tile the following board of length n-3 with squares and trominos:
._ _
|_|_|
|_|_|_ _ _ _ _
|_|_|_|_|_|_|_| ... . (End)
|
|
LINKS
|
|
|
FORMULA
|
G.f.: x - x^2*(1+2*x^2) / ( -1+x+x^3 ). - R. J. Mathar, Oct 30 2011
|
|
MATHEMATICA
|
|
|
PROG
|
(Haskell)
a179070 n = a179070_list !! (n-1)
a179070_list = 1 : zs where zs = 1 : 1 : 3 : zipWith (+) zs (drop 2 zs)
|
|
CROSSREFS
|
Coordination sequences for triangular tilings of hyperbolic space: A001630, A007283, A054886, A078042, A096231, A163876, A179070, A265057, A265058, A265059, A265060, A265061, A265062, A265063, A265064, A265065, A265066, A265067, A265068, A265069, A265070, A265071, A265072, A265073, A265074, A265075, A265076, A265077.
|
|
KEYWORD
|
easy,nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|
|
A054886
|
|
Layer counting sequence for hyperbolic tessellation by cuspidal triangles of angles (Pi/3,Pi/3,0) (this is the classical modular tessellation).
|
|
+10
37
|
|
|
1, 3, 6, 10, 16, 26, 42, 68, 110, 178, 288, 466, 754, 1220, 1974, 3194, 5168, 8362, 13530, 21892, 35422, 57314, 92736, 150050, 242786, 392836, 635622, 1028458, 1664080, 2692538, 4356618, 7049156, 11405774, 18454930, 29860704, 48315634, 78176338, 126491972
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
COMMENTS
|
The layer sequence is the sequence of the cardinalities of the layers accumulating around a ( finite-sided ) polygon of the tessellation under successive side-reflections; see the illustration accompanying A054888.
Equivalently, coordination sequence for (3,3,infinity) tiling of hyperbolic plane. - N. J. A. Sloane, Dec 29 2015
Equivalently, spherical growth series for modular group.
Also, number of sequences of length n with terms 1, 2, and 3, with no adjacent terms equal, and no three consecutive terms (1, 2, 3) or (3, 2, 1). - Pontus von Brömssen, Jan 03 2022
|
|
REFERENCES
|
P. de la Harpe, Topics in Geometric Group Theory, Univ. Chicago Press, 2000, p. 156.
|
|
LINKS
|
|
|
FORMULA
|
G.f.: (1+2*x+2*x^2+x^3)/(1-x-x^2) = (x^2+x+1)*(1+x)/(1-x-x^2).
a(n) = 2*F(n+2) for n >= 2, with F(n) the n-th Fibonacci number (cf. A000045).
E.g.f.: 2*exp(x/2)*(5*cosh(sqrt(5)*x/2) + 3*sqrt(5)*sinh(sqrt(5)*x/2))/5 - 1 - x. - Stefano Spezia, Apr 18 2022
|
|
MATHEMATICA
|
Join[{1, 3}, 2Fibonacci[Range[4, 40]]] (* Harvey P. Dale, Jan 06 2012 *)
|
|
PROG
|
(PARI) my(x='x+O('x^50)); Vec((1+2*x+2*x^2+x^3)/(1-x-x^2)) \\ G. C. Greubel, Aug 06 2017
|
|
CROSSREFS
|
Coordination sequences for triangular tilings of hyperbolic space: A001630, A007283, A054886, A078042, A096231, A163876, A179070, A265057, A265058, A265059, A265060, A265061, A265062, A265063, A265064, A265065, A265066, A265067, A265068, A265069, A265070, A265071, A265072, A265073, A265074, A265075, A265076, A265077.
|
|
KEYWORD
|
nonn,easy,nice
|
|
AUTHOR
|
Paolo Dominici (pl.dm(AT)libero.it), May 23 2000
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|
|
A078042
|
|
Expansion of (1-x)/(1+x-x^2+x^3).
|
|
+10
30
|
|
|
1, -2, 3, -6, 11, -20, 37, -68, 125, -230, 423, -778, 1431, -2632, 4841, -8904, 16377, -30122, 55403, -101902, 187427, -344732, 634061, -1166220, 2145013, -3945294, 7256527, -13346834, 24548655, -45152016, 83047505, -152748176, 280947697, -516743378, 950439251, -1748130326, 3215312955
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
COMMENTS
|
Absolute values give coordination sequence for (3,infinity,infinity) tiling of hyperbolic plane. - N. J. A. Sloane, Dec 29 2015
a(n) is the upper left entry of the n-th power of the 3 X 3 matrix M = [-2, -2, 1; 1, 1, 0; 1, 0, 0]; a(n) = M^n [1, 1]. - Philippe Deléham, Apr 19 2023
|
|
LINKS
|
|
|
FORMULA
|
a(n) = -a(n-1) + a(n-2) - a(n-3) for n > 2; a(0)=1, a(1)=-2, a(2)=3. - Harvey P. Dale, Jun 01 2012
|
|
MATHEMATICA
|
CoefficientList[Series[(1-x)/(1+x-x^2+x^3), {x, 0, 40}], x] (* or *) LinearRecurrence[{-1, 1, -1}, {1, -2, 3}, 40] (* Harvey P. Dale, Jun 01 2012 *)
|
|
PROG
|
(Magma) [n le 3 select -n*(-1)^n else -Self(n-1)+Self(n-2)-Self(n-3): n in [1..50]]; // Vincenzo Librandi, Dec 30 2015
|
|
CROSSREFS
|
Coordination sequences for triangular tilings of hyperbolic space: A001630, A007283, A054886, A078042, A096231, A163876, A179070, A265057, A265058, A265059, A265060, A265061, A265062, A265063, A265064, A265065, A265066, A265067, A265068, A265069, A265070, A265071, A265072, A265073, A265074, A265075, A265076, A265077.
|
|
KEYWORD
|
sign,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|
|
A182097
|
|
Expansion of 1/(1-x^2-x^3).
|
|
+10
29
|
|
|
1, 0, 1, 1, 1, 2, 2, 3, 4, 5, 7, 9, 12, 16, 21, 28, 37, 49, 65, 86, 114, 151, 200, 265, 351, 465, 616, 816, 1081, 1432, 1897, 2513, 3329, 4410, 5842, 7739, 10252, 13581, 17991, 23833, 31572, 41824, 55405, 73396, 97229, 128801, 170625, 226030, 299426, 396655, 525456, 696081, 922111, 1221537, 1618192, 2143648, 2839729, 3761840, 4983377, 6601569, 8745217
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,6
|
|
COMMENTS
|
Number of compositions (ordered partitions) into parts 2 and 3. - Joerg Arndt, Aug 21 2013
a(n) is the top left entry of the n-th power of any of the 3X3 matrices [0, 1, 1; 0, 0, 1; 1, 0, 0], [0, 1, 0; 1, 0, 1; 1, 0, 0], [0, 1, 1; 1, 0, 0; 0, 1, 0] or [0, 0, 1; 1, 0, 0; 1, 1, 0]. - R. J. Mathar, Feb 03 2014
Conjectured values of d(n), the dimension of a Z-module in MZV(conv). See the Waldschmidt link. - Michael Somos, Mar 14 2014
Shannon et al. (2006) call these the Van der Laan numbers. - N. J. A. Sloane, Jan 11 2022
|
|
REFERENCES
|
A. G. Shannon, P. G. Anderson and A. F. Horadam, Properties of Cordonnier, Perrin and Van der Laan numbers, International Journal of Mathematical Education in Science and Technology, Volume 37:7 (2006), 825-831. See R_n.
Michel Waldschmidt, "Multiple Zeta values and Euler-Zagier numbers", in Number theory and discrete mathematics, International conference in honour of Srinivasa Ramanujan, Center for Advanced Study in Mathematics, Panjab University, Chandigarh, (Oct 02, 2000).
|
|
LINKS
|
Michel Waldschmidt, Multiple Zeta values and Euler-Zagier numbers, Slides, Number theory and discrete mathematics, International conference in honour of Srinivasa Ramanujan, Center for Advanced Study in Mathematics, Panjab University, Chandigarh, (Oct 02, 2000).
|
|
FORMULA
|
G.f.: 1 / (1 - x^2 - x^3).
a(n) = a(n-2) + a(n-3) for all n in Z.
a(n+4) = A164001(n) if n > 1. - (End)
|
|
EXAMPLE
|
G.f. = 1 + x^2 + x^3 + x^4 + 2*x^5 + 2*x^6 + 3*x^7 + 4*x^8 + 5*x^9 + ...
|
|
MATHEMATICA
|
a[ n_] := If[n < 0, SeriesCoefficient[ (1 + x) / (1 + x - x^3), {x, 0, -n}], SeriesCoefficient[ 1 / (1 - x^2 - x^3), {x, 0, n}]]; (* Michael Somos, Dec 13 2013 *)
CoefficientList[Series[1/(1-x^2-x^3), {x, 0, 60}], x] (* or *) LinearRecurrence[ {0, 1, 1}, {1, 0, 1}, 70] (* Harvey P. Dale, Dec 04 2014 *)
|
|
PROG
|
(PARI) {a(n) = if( n<0, polcoeff( (1 + x) / (1 + x - x^3) + x * O(x^-n), -n), polcoeff( 1 / (1 - x^2 - x^3) + x * O(x^n), n))}; /* Michael Somos, Dec 13 2013 */
(PARI) Vec(1/(1-x^2-x^3) + O(x^99)) \\ Altug Alkan, Sep 02 2016
(Magma) m:=50; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!(1/(1-x^2-x^3))); // G. C. Greubel, Aug 11 2018
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|
|
A163876
|
|
Number of reduced words of length n in Coxeter group on 3 generators S_i with relations (S_i)^2 = (S_i S_j)^6 = I.
|
|
+10
27
|
|
|
1, 3, 6, 12, 24, 48, 93, 180, 351, 684, 1332, 2592, 5046, 9825, 19128, 37239, 72498, 141144, 274788, 534972, 1041513, 2027676, 3947595, 7685400, 14962368, 29129580, 56711106, 110408373, 214949232, 418475259, 814711182, 1586125572, 3087958512
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
COMMENTS
|
Also, coordination sequence for (6,6,6) tiling of hyperbolic plane. - N. J. A. Sloane, Dec 29 2015
The initial terms coincide with those of A003945, although the two sequences are eventually different.
Computed with MAGMA using commands similar to those used to compute A154638.
|
|
LINKS
|
|
|
FORMULA
|
G.f.: (x^6 + 2*x^5 + 2*x^4 + 2*x^3 + 2*x^2 + 2*x + 1)/(x^6 - x^5 - x^4 - x^3 - x^2 - x + 1).
G.f.: (1+x)*(1-x^6)/(1-2*x+2*x^6-x^7). - G. C. Greubel, Apr 25 2019
|
|
MATHEMATICA
|
CoefficientList[Series[(1+x)*(1-x^6)/(1-2*x+2*x^6-x^7), {x, 0, 40}], x] (* G. C. Greubel, Aug 06 2017, modified Apr 25 2019 *)
|
|
PROG
|
(PARI) x='x+O('x^40); Vec((x^6+2*x^5+2*x^4+2*x^3+2*x^2+2*x+1)/(x^6-x^5- x^4-x^3-x^2-x+1)) \\ G. C. Greubel, Aug 06 2017
(Magma) R<x>:=PowerSeriesRing(Integers(), 40); Coefficients(R!( (1+x)*(1-x^6)/(1-2*x+2*x^6-x^7) )); // G. C. Greubel, Apr 25 2019
(Sage) ((1+x)*(1-x^6)/(1-2*x+2*x^6-x^7)).series(x, 40).coefficients(x, sparse=False) # G. C. Greubel, Apr 25 2019
|
|
CROSSREFS
|
Coordination sequences for triangular tilings of hyperbolic space: A001630, A007283, A054886, A078042, A096231, A163876, A179070, A265057, A265058, A265059, A265060, A265061, A265062, A265063, A265064, A265065, A265066, A265067, A265068, A265069, A265070, A265071, A265072, A265073, A265074, A265075, A265076, A265077.
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|
|
A265057
|
|
Coordination sequence for (2,3,7) tiling of hyperbolic plane.
|
|
+10
27
|
|
|
1, 3, 5, 7, 9, 12, 16, 20, 24, 28, 33, 40, 48, 57, 67, 78, 92, 109, 129, 152, 178, 209, 246, 290, 342, 402, 472, 555, 653, 769, 905, 1064, 1251, 1471, 1731, 2037, 2396, 2818, 3314, 3898, 4586, 5395, 6346, 7464, 8779, 10327, 12148, 14290, 16809, 19771, 23256, 27356, 32179, 37852, 44524, 52372
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
LINKS
|
Index entries for linear recurrences with constant coefficients, signature (-1, 0, 1, 1, 1, 1, 1, 0, -1, -1).
|
|
FORMULA
|
G.f.: (x^6+x^5+x^4+x^3+x^2+x+1)*(x^2+x+1)*(x+1)^2/(x^10+x^9-x^7-x^6-x^5-x^4-x^3+x+1).
|
|
MATHEMATICA
|
CoefficientList[Series[(x^6 + x^5 + x^4 + x^3 + x^2 + x + 1) (x^2 + x + 1) (x + 1)^2/(x^10 + x^9 - x^7 - x^6 - x^5 - x^4 - x^3 + x + 1), {x, 0, 60}], x] (* Vincenzo Librandi, Dec 30 2015 *)
|
|
PROG
|
(PARI) x='x+O('x^50); Vec((x^6+x^5+x^4+x^3+x^2+x+1)*(x^2+x+1)*(x+1)^2/(x^10+x^9-x^7-x^6-x^5-x^4-x^3+x+1)) \\ G. C. Greubel, Aug 06 2017
|
|
CROSSREFS
|
Coordination sequences for triangular tilings of hyperbolic space: A001630, A007283, A054886, A078042, A096231, A163876, A179070, A265057, A265058, A265059, A265060, A265061, A265062, A265063, A265064, A265065, A265066, A265067, A265068, A265069, A265070, A265071, A265072, A265073, A265074, A265075, A265076, A265077.
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|
|
A265058
|
|
Coordination sequence for (2,3,8) tiling of hyperbolic plane.
|
|
+10
27
|
|
|
1, 3, 5, 7, 9, 12, 16, 21, 27, 33, 40, 49, 61, 76, 94, 116, 142, 174, 214, 264, 326, 401, 493, 606, 745, 917, 1129, 1390, 1710, 2103, 2587, 3183, 3917, 4820, 5931, 7297, 8977, 11045, 13590, 16722, 20575, 25315, 31147, 38322, 47151, 58015, 71382, 87828, 108062, 132958, 163590, 201280, 247654
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
LINKS
|
Index entries for linear recurrences with constant coefficients, signature (0, 0, 1, 0, 1, 0, 1, 0, 0, -1).
|
|
FORMULA
|
G.f.: (x+1)^2*(x^2+x+1)*(x^6+x^4+x^2+1)/(x^10-x^7-x^5-x^3+1).
|
|
MATHEMATICA
|
CoefficientList[Series[(x + 1)^2 (x^2 + x + 1) (x^6 + x^4 + x^2 + 1)/(x^10 - x^7 - x^5 - x^3 + 1), {x, 0, 60}], x] (* Vincenzo Librandi, Dec 30 2015 *)
|
|
PROG
|
(PARI) x='x+O('x^50); Vec((x+1)^2*(x^2+x+1)*(x^6+x^4+x^2+1)/(x^10-x^7-x^5-x^3+1)) \\ G. C. Greubel, Aug 06 2017
|
|
CROSSREFS
|
Coordination sequences for triangular tilings of hyperbolic space: A001630, A007283, A054886, A078042, A096231, A163876, A179070, A265057, A265058, A265059, A265060, A265061, A265062, A265063, A265064, A265065, A265066, A265067, A265068, A265069, A265070, A265071, A265072, A265073, A265074, A265075, A265076, A265077.
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
Search completed in 0.037 seconds
|