[go: up one dir, main page]

login
Search: a015531 -id:a015531
     Sort: relevance | references | number | modified | created      Format: long | short | data
Triangle read by rows, giving the numbers T(n,m) = binomial(n+1, m+1); or, Pascal's triangle A007318 with its left-hand edge removed.
+10
48
1, 2, 1, 3, 3, 1, 4, 6, 4, 1, 5, 10, 10, 5, 1, 6, 15, 20, 15, 6, 1, 7, 21, 35, 35, 21, 7, 1, 8, 28, 56, 70, 56, 28, 8, 1, 9, 36, 84, 126, 126, 84, 36, 9, 1, 10, 45, 120, 210, 252, 210, 120, 45, 10, 1, 11, 55, 165, 330, 462, 462, 330, 165, 55, 11, 1, 12, 66, 220, 495, 792, 924, 792
OFFSET
0,2
COMMENTS
T(n,m) is the number of m-faces of a regular n-simplex.
An n-simplex is the n-dimensional analog of a triangle. Specifically, a simplex is the convex hull of a set of (n + 1) affinely independent points in some Euclidean space of dimension n or higher, i.e., a set of points such that no m-plane contains more than (m + 1) of them. Such points are said to be in general position.
Reversing the rows gives A074909, which as a linear sequence is essentially the same as this.
From Tom Copeland, Dec 07 2007: (Start)
T(n,k) * (k+1)! = A068424. The comment on permuted words in A068424 shows that T is related to combinations of letters defined by connectivity of regular polytope simplexes.
If T is the diagonally-shifted Pascal matrix, binomial(n+m, k+m), for m=1, then T is a fundamental type of matrix that is discussed in A133314 and the following hold.
The infinitesimal matrix generator is given by A132681, so T = LM(1) of A132681 with inverse LM(-1).
With a(k) = (-x)^k / k!, T * a = [ Laguerre(n,x,1) ], a vector array with index n for the Laguerre polynomials of order 1. Other formulas for the action of T are given in A132681.
T(n,k) = (1/n!) (D_x)^n (D_t)^k Gf(x,t) evaluated at x=t=0 with Gf(x,t) = exp[ t * x/(1-x) ] / (1-x)^2.
[O.g.f. for T ] = 1 / { [ 1 - t * x/(1-x) ] * (1-x)^2 }. [ O.g.f. for row sums ] = 1 / { (1-x) * (1-2x) }, giving A000225 (without a leading zero) for the row sums. Alternating sign row sums are all 1. [Sign correction noted by Vincent J. Matsko, Jul 19 2015]
O.g.f. for row polynomials = [ (1+q)**(n+1) - 1 ] / [ (1+q) -1 ] = A(1,n+1,q) on page 15 of reference on Grassmann cells in A008292. (End)
Given matrices A and B with A(n,k) = T(n,k)*a(n-k) and B(n,k) = T(n,k)*b(n-k), then A*B = C where C(n,k) = T(n,k)*[a(.)+b(.)]^(n-k), umbrally. The e.g.f. for the row polynomials of A is {(a+t) exp[(a+t)x] - a exp(a x)}/t, umbrally. - Tom Copeland, Aug 21 2008
A007318*A097806 as infinite lower triangular matrices. - Philippe Deléham, Feb 08 2009
Riordan array (1/(1-x)^2, x/(1-x)). - Philippe Deléham, Feb 22 2012
The elements of the matrix inverse are T^(-1)(n,k)=(-1)^(n+k)*T(n,k). - R. J. Mathar, Mar 12 2013
Relation to K-theory: T acting on the column vector (-0,d,-d^2,d^3,...) generates the Euler classes for a hypersurface of degree d in CP^n. Cf. Dugger p. 168 and also A104712, A111492, and A238363. - Tom Copeland, Apr 11 2014
Number of walks of length p>0 between any two distinct vertices of the complete graph K_(n+2) is W(n+2,p)=(-1)^(p-1)*Sum_{k=0..p-1} T(p-1,k)*(-n-2)^k = ((n+1)^p - (-1)^p)/(n+2) = (-1)^(p-1)*Sum_{k=0..p-1} (-n-1)^k. This is equal to (-1)^(p-1)*Phi(p,-n-1), where Phi is the cyclotomic polynomial when p is an odd prime. For K_3, see A001045; for K_4, A015518; for K_5, A015521; for K_6, A015531; for K_7, A015540. - Tom Copeland, Apr 14 2014
Consider the transformation 1 + x + x^2 + x^3 + ... + x^n = A_0*(x-1)^0 + A_1*(x-1)^1 + A_2*(x-1)^2 + ... + A_n*(x-1)^n. This sequence gives A_0, ..., A_n as the entries in the n-th row of this triangle, starting at n = 0. - Derek Orr, Oct 14 2014
See A074909 for associations among this array, the Bernoulli polynomials and their umbral compositional inverses, and the face polynomials of permutahedra and their duals (cf. A019538). - Tom Copeland, Nov 14 2014
From Wolfdieter Lang, Dec 10 2015: (Start)
A(r, n) = T(n+r-2, r-1) = risefac(n,r)/r! = binomial(n+r-1, r), for n >= 1 and r >= 1, gives the array with the number of independent components of a symmetric tensors of rank r (number of indices) and dimension n (indices run from 1 to n). Here risefac(n, k) is the rising factorial.
As(r, n) = T(n+1, r+1) = fallfac(n, r)/r! = binomial(n, r), r >= 1 and n >= 1 (with the triangle entries T(n, k) = 0 for n < k) gives the array with the number of independent components of an antisymmetric tensor of rank r and dimension n. Here fallfac is the falling factorial. (End)
The h-vectors associated to these f-vectors are given by A000012 regarded as a lower triangular matrix. Read as bivariate polynomials, the h-polynomials are the complete homogeneous symmetric polynomials in two variables, found in the compositional inverse of an e.g.f. for A008292, the h-vectors of the permutahedra. - Tom Copeland, Jan 10 2017
For a correlation between the states of a quantum system and the combinatorics of the n-simplex, see Boya and Dixit. - Tom Copeland, Jul 24 2017
LINKS
Paul Barry, On the f-Matrices of Pascal-like Triangles Defined by Riordan Arrays, arXiv:1805.02274 [math.CO], 2018.
L. Boya and K. Dixit, Geometry of density states, arXiv:808.1930 [quant-phy], 2017.
V. Buchstaber, Lectures on Toric Topology, Trends in Mathematics - New Series, Information Center for Mathematical Sciences, Vol. 10, No. 1, 2008. p. 7.
D. Dugger, A Geometric Introduction to K-Theory [From Tom Copeland, Apr 11 2014]
Atli Fannar Franklín, Pattern avoidance enumerated by inversions, arXiv:2410.07467 [math.CO], 2024. See pp. 2, 12.
Atli Fannar Franklín, Anders Claesson, Christian Bean, Henning Úlfarsson, and Jay Pantone, Restricted Permutations Enumerated by Inversions, arXiv:2406.16403 [cs.DM], 2024. See p. 4.
B. Grünbaum and G. C. Shephard, Convex polytopes, Bull. London Math. Soc. (1969) 1 (3): 257-300.
G. Hetyei, Meixner polynomials of the second kind and quantum algebras representing su(1,1), arXiv preprint arXiv:0909.4352 [math.QA], 2009, p. 4 (Added by Tom Copeland, Oct 01 2015).
Wolfdieter Lang, On Generating functions of Diagonals Sequences of Sheffer and Riordan Number Triangles, arXiv:1708.01421 [math.NT], August 2017.
Wikipedia, Simplex
FORMULA
T(n, k) = Sum_{j=k..n} binomial(j,k) = binomial(n+1, k+1), n >= k >= 0, else 0. (Partial sum of column k of A007318 (Pascal), or summation on the upper binomial index (Graham et al. (GKP), eq. (5.10). For the GKP reference see A007318.) - Wolfdieter Lang, Aug 22 2012
E.g.f.: 1/x*((1 + x)*exp(t*(1 + x)) - exp(t)) = 1 + (2 + x)*t + (3 + 3*x + x^2)*t^2/2! + .... The infinitesimal generator for this triangle has the sequence [2,3,4,...] on the main subdiagonal and 0's elsewhere. - Peter Bala, Jul 16 2013
T(n,k) = 2*T(n-1,k) + T(n-1,k-1) - T(n-2,k) - T(n-2,k-1), T(0,0)=1, T(1,0)=2, T(1,1)=1, T(n,k)=0 if k<0 or if k>n. - Philippe Deléham, Dec 27 2013
T(n,k) = A193862(n,k)/2^k. - Philippe Deléham, Jan 29 2014
G.f.: 1/((1-x)*(1-x-x*y)). - Philippe Deléham, Mar 13 2014
From Tom Copeland, Mar 26 2014: (Start)
[From Copeland's 2007 and 2008 comments]
A) O.g.f.: 1 / { [ 1 - t * x/(1-x) ] * (1-x)^2 } (same as Deleham's).
B) The infinitesimal generator for T is given in A132681 with m=1 (same as Bala's), which makes connections to the ubiquitous associated Laguerre polynomials of integer orders, for this case the Laguerre polynomials of order one L(n,-t,1).
C) O.g.f. of row e.g.f.s: Sum_{n>=0} L(n,-t,1) x^n = exp[t*x/(1-x)]/(1-x)^2 = 1 + (2+t)x + (3+3*t+t^2/2!)x^2 + (4+6*t+4*t^2/2!+t^3/3!)x^3+ ... .
D) E.g.f. of row o.g.f.s: ((1+t)*exp((1+t)*x)-exp(x))/t (same as Bala's).
E) E.g.f. for T(n,k)*a(n-k): {(a+t) exp[(a+t)x] - a exp(a x)}/t, umbrally. For example, for a(k)=2^k, the e.g.f. for the row o.g.f.s is {(2+t) exp[(2+t)x] - 2 exp(2x)}/t.
(End)
From Tom Copeland, Apr 28 2014: (Start)
With different indexing
A) O.g.f. by row: [(1+t)^n-1]/t.
B) O.g.f. of row o.g.f.s: {1/[1-(1+t)*x] - 1/(1-x)}/t.
C) E.g.f. of row o.g.f.s: {exp[(1+t)*x]-exp(x)}/t.
These generating functions are related to row e.g.f.s of A111492. (End)
From Tom Copeland, Sep 17 2014: (Start)
A) U(x,s,t)= x^2/[(1-t*x)(1-(s+t)x)] = Sum_{n >= 0} F(n,s,t)x^(n+2) is a generating function for bivariate row polynomials of T, e.g., F(2,s,t)= s^2 + 3s*t + 3t^2 (Buchstaber, 2008).
B) dU/dt=x^2 dU/dx with U(x,s,0)= x^2/(1-s*x) (Buchstaber, 2008).
C) U(x,s,t) = exp(t*x^2*d/dx)U(x,s,0) = U(x/(1-t*x),s,0).
D) U(x,s,t) = Sum[n >= 0, (t*x)^n L(n,-:xD:,-1)] U(x,s,0), where (:xD:)^k=x^k*(d/dx)^k and L(n,x,-1) are the Laguerre polynomials of order -1, related to normalized Lah numbers. (End)
E.g.f. satisfies the differential equation d/dt(e.g.f.(x,t)) = (x+1)*e.g.f.(x,t) + exp(t). - Vincent J. Matsko, Jul 18 2015
The e.g.f. of the Norlund generalized Bernoulli (Appell) polynomials of order m, NB(n,x;m), is given by exponentiation of the e.g.f. of the Bernoulli numbers, i.e., multiple binomial self-convolutions of the Bernoulli numbers, through the e.g.f. exp[NB(.,x;m)t] = (t/(e^t - 1))^(m+1) * e^(xt). Norlund gave the relation to the factorials (x-1)!/(x-1-n)! = (x-1) ... (x-n) = NB(n,x;n), so T(n,m) = NB(m+1,n+2;m+1)/(m+1)!. - Tom Copeland, Oct 01 2015
From Wolfdieter Lang, Nov 08 2018: (Start)
Recurrences from the A- and Z- sequences for the Riordan triangle (see the W. Lang link under A006232 with references), which are A(n) = A019590(n+1), [1, 1, repeat (0)] and Z(n) = (-1)^(n+1)*A054977(n), [2, repeat(-1, 1)]:
T(0, 0) = 1, T(n, k) = 0 for n < k, and T(n, 0) = Sum_{j=0..n-1} Z(j)*T(n-1, j), for n >= 1, and T(n, k) = T(n-1, k-1) + T(n-1, k), for n >= m >= 1.
Boas-Buck recurrence for columns (see the Aug 10 2017 remark in A036521 also for references):
T(n, k) = ((2 + k)/(n - k))*Sum_{j=k..n-1} T(j, k), for n >= 1, k = 0, 1, ..., n-1, and input T(n, n) = 1, for n >= 0, (the BB-sequences are alpha(n) = 2 and beta(n) = 1). (End)
T(n, k) = [x^k] Sum_{j=0..n} (x+1)^j. - Peter Luschny, Jul 09 2019
EXAMPLE
The triangle T(n, k) begins:
n\k 0 1 2 3 4 5 6 7 8 9 10 11 ...
0: 1
1: 2 1
2: 3 3 1
3: 4 6 4 1
4: 5 10 10 5 1
5: 6 15 20 15 6 1
6: 7 21 35 35 21 7 1
7: 8 28 56 70 56 28 8 1
8: 9 36 84 126 126 84 36 9 1
9: 10 45 120 210 252 210 120 45 10 1
10: 11 55 165 330 462 462 330 165 55 11 1
11: 12 66 220 495 792 924 792 495 220 66 12 1
... reformatted by Wolfdieter Lang, Mar 23 2015
Production matrix begins
2 1
-1 1 1
1 0 1 1
-1 0 0 1 1
1 0 0 0 1 1
-1 0 0 0 0 1 1
1 0 0 0 0 0 1 1
-1 0 0 0 0 0 0 1 1
1 0 0 0 0 0 0 0 1 1
- Philippe Deléham, Jan 29 2014
From Wolfdieter Lang, Nov 08 2018: (Start)
Recurrence [Philippe Deléham]: T(7, 3) = 2*35 + 35 - 15 - 20 = 70.
Recurrence from Riordan A- and Z-sequences: [1,1,repeat(0)] and [2, repeat(-1, +1)]: From Z: T(5, 0) = 2*5 - 10 + 10 - 5 + 1 = 6. From A: T(7, 3) = 35 + 35 = 70.
Boas-Buck column k=3 recurrence: T(7, 3) = (5/4)*(1 + 5 + 15 + 35) = 70. (End)
MAPLE
for i from 0 to 12 do seq(binomial(i, j)*1^(i-j), j = 1 .. i) od;
MATHEMATICA
Flatten[Table[CoefficientList[D[1/x ((x + 1) Exp[(x + 1) z] - Exp[z]), {z, k}] /. z -> 0, x], {k, 0, 11}]]
CoefficientList[CoefficientList[Series[1/((1 - x)*(1 - x - x*y)), {x, 0, 10}, {y, 0, 10}], x], y] // Flatten (* G. C. Greubel, Nov 22 2017 *)
PROG
(PARI) for(n=0, 20, for(k=0, n, print1(1/k!*sum(i=0, n, (prod(j=0, k-1, i-j))), ", "))) \\ Derek Orr, Oct 14 2014
(Sage)
Trow = lambda n: sum((x+1)^j for j in (0..n)).list()
for n in (0..10): print(Trow(n)) # Peter Luschny, Jul 09 2019
KEYWORD
easy,nonn,tabl
AUTHOR
Zerinvary Lajos, Dec 02 2007
EXTENSIONS
Edited by Tom Copeland and N. J. A. Sloane, Dec 11 2007
STATUS
approved
a(p^k) = p^k - p^(k-1) + p^(k-2) - ... +- 1, and then extend by multiplicativity.
+10
42
1, 1, 2, 3, 4, 2, 6, 5, 7, 4, 10, 6, 12, 6, 8, 11, 16, 7, 18, 12, 12, 10, 22, 10, 21, 12, 20, 18, 28, 8, 30, 21, 20, 16, 24, 21, 36, 18, 24, 20, 40, 12, 42, 30, 28, 22, 46, 22, 43, 21, 32, 36, 52, 20, 40, 30, 36, 28, 58, 24, 60, 30, 42, 43, 48, 20, 66, 48, 44, 24, 70, 35
OFFSET
1,3
COMMENTS
For more information see the Comments in A061020.
a(n) is the number of integers j such that 1 <= j <= n and gcd(n,j) is a perfect square. For example, a(12) = 6 because |{1,4,5,7,8,11}|=6 and the respective GCDs with 12 are 1,4,1,1,4,1, which are squares. - Geoffrey Critzer, Feb 16 2015
If m is squarefree (A005117), then a(m) = A000010(m) where A000010 is the Euler totient function. - Michel Marcus, Nov 08 2017
Also it appears that the primorials (A002110) is the sequence of indices of minimum records for a(n)/n, and these records are A038110(n)/A060753(n). - Michel Marcus, Nov 09 2017
Also called rho(n). When rho(n) | n, then n is called k-imperfect, with k = n/rho(n), cf. A127724. - M. F. Hasler, Feb 13 2020
REFERENCES
P. J. McCarthy, Introduction to Arithmetical Functions, Springer Verlag, 1986, page 25.
LINKS
László Tóth, A survey of the alternating sum-of-divisors function, arXiv:1111.4842 [math.NT], 2011-2014.
FORMULA
a(n) = abs(A061020(n)).
a(n) = n*Sum_{d|n} lambda(d)/d, where lambda(n) is A008836(n). - Enrique Pérez Herrero, Sep 23 2012
Dirichlet g.f.: zeta(s - 1)*zeta(2*s)/zeta(s). - Geoffrey Critzer, Feb 25 2015
From Michel Marcus, Nov 05 2017: (Start)
a(2^n) = A001045(n+1);
a(3^n) = A015518(n+1);
a(5^n) = A015531(n+1);
a(7^n) = A015552(n+1);
a(11^n) = A015592(n+1). (End)
a(p^k) = p^k - a(p^(k - 1)) for k > 0 and prime p. - David A. Corneth, Nov 09 2017
a(n) = Sum_{d|n, d is a perfect square} phi(n/d), where phi(k) is the Euler totient function. - Daniel Suteu, Jun 27 2018
a(p^k) = A071324(p^k), for k >= 0 and prime p. - Michel Marcus, Aug 11 2018
Sum_{k=1..n} a(k) ~ Pi^2 * n^2 / 30. - Vaclav Kotesovec, Feb 07 2019
G.f.: Sum_{k>=1} lambda(k)*x^k/(1 - x^k)^2. - Ilya Gutkovskiy, May 23 2019
a(n) = Sum_{i=1..n} A010052(gcd(n,i)). - Ridouane Oudra, Nov 24 2019
a(p^k) = round(p^(k+1)/(p+1)). - M. F. Hasler, Feb 13 2020
MAPLE
a:= n-> mul(add(i[1]^(i[2]-j)*(-1)^j, j=0..i[2]), i=ifactors(n)[2]):
seq(a(n), n=1..100); # Alois P. Heinz, Nov 03 2017
MATHEMATICA
Table[Length[Select[Range[n], IntegerQ[GCD[n, #]^(1/2)] &]], {n, 72}] (* Geoffrey Critzer, Feb 16 2015 *)
a[n_] := n*DivisorSum[n, LiouvilleLambda[#]/#&]; Array[a, 72] (* Jean-François Alcover, Dec 04 2017, after Enrique Pérez Herrero *)
f[p_, e_] := Sum[(-1)^(e-k)*p^k, {k, 0, e}]; a[1] = 1; a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 100] (* Amiram Eldar, Jan 01 2020 *)
PROG
(Haskell)
a206369 n = product $
zipWith h (a027748_row n) (map toInteger $ a124010_row n) where
h p e = sum $ take (fromInteger e + 1) $
iterate ((* p) . negate) (1 - 2 * (e `mod` 2))
-- Reinhard Zumkeller, Feb 08 2012
(PARI) a(n) = sum(k=1, n, issquare(gcd(n, k)));
(PARI) ak(p, e)=my(s=1); for(i=1, e, s=s*p + (-1)^i); s
a(n)=my(f=factor(n)); prod(i=1, #f~, ak(f[i, 1], f[i, 2])) \\ Charles R Greathouse IV, Dec 27 2016
(PARI) a(n) = sumdiv(n, d, eulerphi(n/d) * issquare(d)); \\ Daniel Suteu, Jun 27 2018
(PARI) apply( {A206369(n)=vecprod([f[1]^(f[2]+1)\/(f[1]+1)|f<-factor(n)~])}, [1..99]) \\ M. F. Hasler, Feb 13 2020
(Python)
from math import prod
from sympy import factorint
def A206369(n): return prod((lambda x:x[0]+int((x[1]<<1)>=p+1))(divmod(p**(e+1), p+1)) for p, e in factorint(n).items()) # Chai Wah Wu, Mar 05 2024
CROSSREFS
Cf. A027748 row, A124010, A206475 (first differences).
Cf. A078429.
Cf. A127724 (k-imperfect), A127725 (2-imperfect), A127726 (3-imperfect).
KEYWORD
nonn,mult
AUTHOR
N. J. A. Sloane, Feb 06 2012
STATUS
approved
Scaled Chebyshev U-polynomials evaluated at i*sqrt(5)/2. Generalized Fibonacci sequence.
+10
29
1, 5, 30, 175, 1025, 6000, 35125, 205625, 1203750, 7046875, 41253125, 241500000, 1413765625, 8276328125, 48450468750, 283633984375, 1660422265625, 9720281250000, 56903517578125, 333118994140625, 1950112558593750, 11416157763671875, 66831351611328125, 391237546875000000
OFFSET
0,2
COMMENTS
a(n) gives the length of the word obtained after n steps with the substitution rule 0->11111, 1->111110, starting from 0. The number of 1's and 0's of this word is 5*a(n-1) and 5*a(n-2), resp.
a(n) / a(n-1) converges to (5 + (3 * sqrt(5))) / 2 as n approaches infinity. (5 + (3 * sqrt(5))) / 2 can also be written as phi^2 + (2 * phi), phi^3 + phi, phi + sqrt(5) + 2, (3 * phi) + 1, (3 * phi^2) - 2, phi^4 - 1 and (5 + (3 * (L(n) / F(n)))) / 2, where L(n) is the n-th Lucas number and F(n) is the n-th Fibonacci number as n approaches infinity. - Ross La Haye, Aug 18 2003, on another version
Pisano period lengths: 1, 3, 3, 6, 1, 3, 24, 12, 9, 3, 10, 6, 56, 24, 3, 24,288, 9, 18, 6, ... - R. J. Mathar, Aug 10 2012
LINKS
Martin Burtscher, Igor Szczyrba, Rafał Szczyrba, Analytic Representations of the n-anacci Constants and Generalizations Thereof, Journal of Integer Sequences, Vol. 18 (2015), Article 15.4.5.
A. F. Horadam, Special properties of the sequence W_n(a,b; p,q), Fib. Quart., 5.5 (1967), 424-434. Case n->n+1, a=0,b=1; p=5, q=5.
Tanya Khovanova, Recursive Sequences
W. Lang, On polynomials related to powers of the generating function of Catalan's numbers, Fib. Quart. 38 (2000) 408-419. Eqs.(39) and (45),rhs, m=5.
Eric Weisstein's World of Mathematics, Horadam Sequence
FORMULA
a(n) = 5*(a(n-1) + a(n-2)), a(-1)=0, a(0)=1.
a(n) = S(n, i*sqrt(5))*(-i*sqrt(5))^n with S(n, x) := U(n, x/2), Chebyshev's polynomials of the 2nd kind, A049310.
G.f.: 1/(1 - 5*x - 5*x^2).
a(n) = (1/3)*Sum_{k=0..n} binomial(n, k)*Fibonacci(k)*3^k. - Benoit Cloitre, Oct 25 2003
a(n) = ((5 + 3*sqrt(5))/2)^n(1/2 + sqrt(5)/6) + (1/2 - sqrt(5)/6)((5 - 3*sqrt(5))/2)^n. - Paul Barry, Sep 22 2004
(a(n)) appears to be given by the floretion - 0.75'i - 0.5'j + 'k - 0.75i' + 0.5j' + 0.5k' + 1.75'ii' - 1.25'jj' + 1.75'kk' - 'ij' - 0.5'ji' - 0.75'jk' - 0.75'kj' - 1.25e ("jes"). - Creighton Dement, Nov 28 2004
a(n) = Sum_{k=0..n} 4^k*A063967(n,k). - Philippe Deléham, Nov 03 2006
G.f.: G(0)/(2-5*x), where G(k)= 1 + 1/(1 - x*(9*k-5)/(x*(9*k+4) - 2/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 17 2013
From Ehren Metcalfe, Nov 18 2017: (Start)
With F(n) = A000045(n), L(n) = A000032(n), beta = (1-sqrt(5))/2:
a(2*n-1) = 5^n*F(4*n)/3 = (5^(n-1/2)*L(4*n) - 2*5^(n-1/2)*beta^(4*n))/3.
a(2*n) = 5^n*L(4*n+2)/3 = (5^(n+1/2)*F(4*n+2) + 2*5^n*beta^(4*n+2))/3.
a(n) = round 5^((n+1)/2)*F(2*(n+1))/3.
a(n) = round 5^(n/2)*L(2*(n+1))/3. (End)
MAPLE
a[0]:=0:a[1]:=1:for n from 2 to 50 do a[n]:=5*a[n-1]+5*a[n-2]od: seq(a[n], n=1..33); # Zerinvary Lajos, Dec 14 2008
MATHEMATICA
LinearRecurrence[{5, 5}, {1, 5}, 30] (* G. C. Greubel, Jan 16 2018 *)
PROG
(Sage) [lucas_number1(n, 5, -5) for n in range(1, 22)] # Zerinvary Lajos, Apr 24 2009
(PARI) x='x+O('x^30); Vec(1/(1 - 5*x - 5*x^2)) \\ G. C. Greubel, Jan 16 2018
(Magma) I:=[1, 5]; [n le 2 select I[n] else 5*Self(n-1) + 5*Self(n-2): n in [0..30]]; // G. C. Greubel, Jan 16 2018
KEYWORD
nonn,easy
AUTHOR
Wolfdieter Lang, Aug 11 2000
STATUS
approved
a(n) = 7*a(n-1) + 8*a(n-2), a(0) = 0, a(1) = 1.
+10
26
0, 1, 7, 57, 455, 3641, 29127, 233017, 1864135, 14913081, 119304647, 954437177, 7635497415, 61083979321, 488671834567, 3909374676537, 31274997412295, 250199979298361, 2001599834386887, 16012798675095097, 128102389400760775, 1024819115206086201, 8198552921648689607
OFFSET
0,3
COMMENTS
A linear 2nd order recurrence. A Jacobsthal number sequence.
Binomial transform of A053573 (preceded by zero). - Paul Barry, Apr 09 2003
Second binomial transform of A080424. Binomial transform of A053573, with leading zero. Binomial transform is 0,1,9,81,729,....(9^n - 0^n)/9. Second binomial transform is 0,1,11,111,1111,... (A002275: repunits). - Paul Barry, Mar 14 2004
Number of walks of length n between any two distinct nodes of the complete graph K_9. Example: a(2)=7 because the walks of length 2 between the nodes A and B of the complete graph ABCDEFGHI are: ACB, ADB, AEB, AFB, AGB, AHB and AIB. - Emeric Deutsch, Apr 01 2004
Unsigned version of A014990. - Philippe Deléham, Feb 13 2007
The ratio a(n+1)/a(n) converges to 8 as n approaches infinity. - Felix P. Muga II, Mar 09 2014
LINKS
Jean-Paul Allouche, Jeffrey Shallit, Zhixiong Wen, Wen Wu, and Jiemeng Zhang, Sum-free sets generated by the period-k-folding sequences and some Sturmian sequences, arXiv:1911.01687 [math.CO], 2019.
M. Dukes and C. D. White, Web Matrices: Structural Properties and Generating Combinatorial Identities, arXiv:1603.01589 [math.CO], 2016.
Dale Gerdemann, Fractal generated from (7,8) recursion, YouTube Video, Dec 5, 2014.
FORMULA
From Paul Barry, Apr 09 2003: (Start)
a(n) = (8^n - (-1)^n)/9.
a(n) = J(3*n)/3 = A001045(3*n)/3. (End)
From Emeric Deutsch, Apr 01 2004: (Start)
a(n) = 8^(n-1) - a(n-1).
G.f.: x/(1-7*x-8*x^2). (End)
a(n) = Sum_{k = 0..n} A106566(n,k)*A099322(k). - Philippe Deléham, Oct 30 2008
a(n) = round(8^n/9). - Mircea Merca, Dec 28 2010
From Peter Bala, May 31 2024: (Start)
G.f: A(x) = x/(1 - x^2) o x/(1 - x^2), where o denotes the black diamond product of power series as defined by Dukes and White. Cf. A054878.
The black diamond product A(x) o A(x) is the g.f. for the number of walks of length n between any two distinct nodes of the complete graph K_81.
Row 8 of A062160. (End)
E.g.f.: exp(-x)*(exp(9*x) - 1)/9. - Elmo R. Oliveira, Aug 17 2024
EXAMPLE
G.f. = x + 7*x^2 + 57*x^3 + 455*x^4 + 3641*x^5 + 29127*x^6 + 233017*x^7 + ...
MAPLE
seq(round(8^n/9), n=0..25); # Mircea Merca, Dec 28 2010
MATHEMATICA
k=0; lst={k}; Do[k=8^n-k; AppendTo[lst, k], {n, 0, 5!}]; lst (* Vladimir Joseph Stephan Orlovsky, Dec 11 2008 *)
LinearRecurrence[{7, 8}, {0, 1}, 30] (* Harvey P. Dale, Mar 04 2016 *)
PROG
(Sage) [lucas_number1(n, 7, -8) for n in range(0, 21)] # Zerinvary Lajos, Apr 24 2009
(Magma) [Round(8^n/9): n in [0..30]]; // Vincenzo Librandi, Jun 24 2011
(PARI) x='x+O('x^30); concat([0], Vec(x/(1-7*x-8*x^2))) \\ G. C. Greubel, Dec 30 2017
KEYWORD
nonn,easy
STATUS
approved
Expansion of x^3/(1 - 2*x + x^3 - 2*x^4) = x^3/( (1-2*x)*(1+x)*(1-x+x^2) ).
+10
26
0, 0, 0, 1, 2, 4, 7, 14, 28, 57, 114, 228, 455, 910, 1820, 3641, 7282, 14564, 29127, 58254, 116508, 233017, 466034, 932068, 1864135, 3728270, 7456540, 14913081, 29826162, 59652324, 119304647, 238609294, 477218588, 954437177, 1908874354, 3817748708
OFFSET
0,5
COMMENTS
A transform of the Jacobsthal numbers. A059633 is the equivalent transform of the Fibonacci numbers.
Paul Curtz, Aug 05 2007, observes that the inverse binomial transform of 0,0,0,1,2,4,7,14,28,57,114,228,455,910,1820,... gives the same sequence up to signs. That is, the extended sequence is an eigensequence for the inverse binomial transform (an autosequence).
The round() function enables the closed (non-recurrence) formula to take a very simple form: see Formula section. This can be generalized without loss of simplicity to a(n) = round(b^n/c), where b and c are very small, incommensurate integers (c may also be an integer fraction). Particular choices of small integers for b and c produce a number of well-known sequences which are usually defined by a recurrence - see Cross Reference. - Ross Drewe, Sep 03 2009
LINKS
M. Bernstein and N. J. A. Sloane, Some canonical sequences of integers, arXiv:math/0205301 [math.CO], 2002. [Link to arXiv version]
M. Bernstein and N. J. A. Sloane, Some canonical sequences of integers, Linear Alg. Applications, 226-228 (1995), 57-72; erratum 320 (2000), 210. [Link to Lin. Alg. Applic. version together with omitted figures]
N. J. A. Sloane, Transforms
FORMULA
a(n) = 2a(n-1) - a(n-3) + 2a(n-4).
a(n) = Sum_{k=0..floor(n/2)} binomial(n-k,k)*A001045(k).
a(n) = Sum_{k=0..n} binomial((n+k)/2,k)*A001045((n-k)/2)*(1+(-1)^(n-k))/2.
a(3n) = A015565(n), a(3n+1) = 2*A015565(n), a(3n+2) = 4*A015565(n). - Paul Curtz, Nov 30 2007
From Paul Curtz, Dec 16 2007: (Start)
a(n+1) - 2a(n) = A131531(n).
a(n) + a(n+3) = 2^n. (End)
a(n) = round(2^n/9). - Ross Drewe, Sep 03 2009
9*a(n) = 2^n + (-1)^n - 3*A010892(n). - R. J. Mathar, Mar 24 2018
MAPLE
A010892 := proc(n) op((n mod 6)+1, [1, 1, 0, -1, -1, 0]) ; end proc:
A113405 := proc(n) (2^n-(-1)^n)/9 -A010892(n-1)/3; end proc: # R. J. Mathar, Dec 17 2010
MATHEMATICA
CoefficientList[Series[x^3/(1-2x+x^3-2x^4), {x, 0, 40}], x] (* or *) LinearRecurrence[{2, 0, -1, 2}, {0, 0, 0, 1}, 40] (* Harvey P. Dale, Apr 30 2011 *)
PROG
(PARI) a(n)=2^n\/9 \\ Charles R Greathouse IV, Jun 05, 2011
(Magma) [Round(2^n/9): n in [0..40]]; // Vincenzo Librandi, Aug 11 2011
CROSSREFS
From Ross Drewe, Sep 03 2009: (Start)
Other sequences a(n) = round(b^n / c), where b and c are very small integers:
A001045 b = 2; c = 3
A007910 b = 2; c = 5
A016029 b = 2; c = 5/3
A077947 b = 2; c = 7
abs(A078043) b = 2; c = 7/3
A007051 b = 3; c = 2
A015518 b = 3; c = 4
A034478 b = 5; c = 2
A003463 b = 5; c = 4
A015531 b = 5; c = 6
(End)
KEYWORD
easy,nonn
AUTHOR
Paul Barry, Oct 28 2005
EXTENSIONS
Edited by N. J. A. Sloane, Dec 13 2007
STATUS
approved
Scaled Chebyshev U-polynomials evaluated at i*sqrt(6)/2. Generalized Fibonacci sequence.
+10
15
1, 6, 42, 288, 1980, 13608, 93528, 642816, 4418064, 30365280, 208700064, 1434392064, 9858552768, 67757668992, 465697330560, 3200729997312, 21998563967232, 151195763787264, 1039165966526976, 7142170381885440
OFFSET
0,2
COMMENTS
a(n) gives the length of the word obtained after n steps with the substitution rule 0->1^6, 1->(1^6)0, starting from 0. The number of 1's and 0's of this word is 6*a(n-1) and 6*a(n-2), resp.
LINKS
Martin Burtscher, Igor Szczyrba, and Rafał Szczyrba, Analytic Representations of the n-anacci Constants and Generalizations Thereof, Journal of Integer Sequences, Vol. 18 (2015), Article 15.4.5.
A. F. Horadam, Special properties of the sequence W_n(a,b; p,q), Fib. Quart., 5.5 (1967), 424-434. Case n->n+1, a=0,b=1; p=6, q=6.
Tanya Khovanova, Recursive Sequences
Wolfdieter Lang, On polynomials related to powers of the generating function of Catalan's numbers, Fib. Quart. 38 (2000) 408-419. Eqs.(39) and (45),rhs, m=6.
FORMULA
a(n) = 6*a(n-1) + 6*a(n-2); a(0)=1, a(1)=6.
a(n) = S(n, i*sqrt(6))*(-i*sqrt(6))^n with S(n, x) := U(n, x/2), Chebyshev's polynomials of the 2nd kind, A049310.
G.f.: 1/(1-6*x-6*x^2).
a(n) = Sum_{k=0..n} 5^k*A063967(n,k). - Philippe Deléham, Nov 03 2006
MATHEMATICA
Join[{a=0, b=1}, Table[c=6*b+6*a; a=b; b=c, {n, 100}]] (* Vladimir Joseph Stephan Orlovsky, Jan 16 2011 *)
LinearRecurrence[{6, 6}, {1, 6}, 40] (* Harvey P. Dale, Nov 05 2011 *)
PROG
(Sage) [lucas_number1(n, 6, -6) for n in range(1, 21)] # Zerinvary Lajos, Apr 24 2009
(Magma) I:=[1, 6]; [n le 2 select I[n] else 6*Self(n-1)+6*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Nov 14 2011
(PARI) x='x+O('x^30); Vec(1/(1-6*x-6*x^2)) \\ G. C. Greubel, Jan 24 2018
KEYWORD
nonn,easy
AUTHOR
Wolfdieter Lang, Aug 11 2000
STATUS
approved
Expansion of x/(1 - 5*x - 4*x^2).
+10
14
0, 1, 5, 29, 165, 941, 5365, 30589, 174405, 994381, 5669525, 32325149, 184303845, 1050819821, 5991314485, 34159851709, 194764516485, 1110461989261, 6331368012245, 36098688018269, 205818912140325, 1173489312774701, 6690722212434805
OFFSET
0,3
COMMENTS
First differences give A122690(n) = {1, 4, 24, 136, 776, 4424, 25224, ...}. Partial sums of a(n) are {0, 1, 6, 35, 200, ...} = (A123270(n) - 1)/8. - Alexander Adamchuk, Nov 03 2006
For n >= 2, a(n) equals the permanent of the (n-1) X (n-1) tridiagonal matrix with 5's along the main diagonal, and 2's along the superdiagonal and the subdiagonal. - John M. Campbell, Jul 19 2011
Pisano period lengths: 1, 1, 8, 1, 4, 8, 48, 1, 24, 4, 40, 8, 42, 48, 8, 2, 72, 24, 360, 4, ... - R. J. Mathar, Aug 10 2012
LINKS
Lucyna Trojnar-Spelina and Iwona Włoch, On Generalized Pell and Pell-Lucas Numbers, Iranian Journal of Science and Technology, Transactions A: Science (2019), 1-7.
FORMULA
a(n) = 5*a(n-1) + 4*a(n-2).
a(n) = Sum_{k=0..floor((n-1)/2)} C(n-k-1, k)*4^k*5^(n-2*k-1). - Paul Barry, Apr 23 2005
a(n) = Sum_{k=0..(n-1)} A122690(k). - Alexander Adamchuk, Nov 03 2006
a(n) = 2^(n-1)*Fibonacci(n, 5/2) = (2/i)^(n-1)*ChebyshevU(n-1, 5*i/4). - G. C. Greubel, Dec 26 2019
MAPLE
seq( simplify((2/I)^(n-1)*ChebyshevU(n-1, 5*I/4)), n=0..20); # G. C. Greubel, Dec 26 2019
MATHEMATICA
LinearRecurrence[{5, 4}, {0, 1}, 30] (* Vincenzo Librandi, Nov 12 2012 *)
Table[2^(n-1)*Fibonacci[n, 5/2], {n, 0, 30}] (* G. C. Greubel, Dec 26 2019 *)
PROG
(Sage) [lucas_number1(n, 5, -4) for n in range(0, 22)] # Zerinvary Lajos, Apr 24 2009
(Magma) [n le 2 select n-1 else 5*Self(n-1)+4*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Nov 12 2012
(PARI) x='x+O('x^30); concat([0], Vec(x/(1-5*x-4*x^2))) \\ G. C. Greubel, Jan 01 2018
(GAP) a:=[0, 1];; for n in [3..30] do a[n]:=5*a[n-1]+4*a[n-2]; od; a; # G. C. Greubel, Dec 26 2019
KEYWORD
nonn,easy
STATUS
approved
a(n) = 6*a(n-1) + 7*a(n-2), a(0) = 0, a(1) = 1.
+10
13
0, 1, 6, 43, 300, 2101, 14706, 102943, 720600, 5044201, 35309406, 247165843, 1730160900, 12111126301, 84777884106, 593445188743, 4154116321200, 29078814248401, 203551699738806, 1424861898171643, 9974033287201500, 69818233010410501, 488727631072873506
OFFSET
0,3
COMMENTS
Number of walks of length n between any two distinct nodes of the complete graph K_8. Example: a(2)=6 because the walks of length 2 between the nodes A and B of the complete graph ABCDEFGH are: ACB, ADB, AEB, AFB, AGB and AHB. - Emeric Deutsch, Apr 01 2004
The ratio a(n+1)/a(n) converges to 7 as n approaches infinity. - Felix P. Muga II, Mar 09 2014
LINKS
Jean-Paul Allouche, Jeffrey Shallit, Zhixiong Wen, Wen Wu, Jiemeng Zhang, Sum-free sets generated by the period-k-folding sequences and some Sturmian sequences, arXiv:1911.01687 [math.CO], 2019.
FORMULA
a(n) = 6*a(n-1) + 7*a(n-2).
From Emeric Deutsch, Apr 01 2004: (Start)
G.f.: x/(1-6*x-7*x^2).
a(n) = 7^(n-1) - a(n-1). (End)
a(n) = (7^n - (-1)^n)/8. - Rolf Pleisch, Jul 06 2009
a(n) = round(7^n/8). - Mircea Merca, Dec 28 2010
E.g.f. exp(3*x)*sinh(4*x)/4. - Elmo R. Oliveira, Aug 17 2024
EXAMPLE
G.f. = x + 6*x^2 + 43*x^3 + 300*x^4 + 2101*x^5 + 14706*x^6 + 102943*x^7 + ...
MAPLE
seq(round(7^n/8), n=0..25); # Mircea Merca, Dec 28 2010
MATHEMATICA
k=0; lst={k}; Do[k=7^n-k; AppendTo[lst, k], {n, 0, 5!}]; lst (* Vladimir Joseph Stephan Orlovsky, Dec 11 2008 *)
Table[(7^n - (-1)^n)/8, {n, 0, 30}] (* G. C. Greubel, Dec 2017 *)
PROG
(PARI) {a(n) = if ( n<0, 0, (7^n - (-1)^n) / 8)};
(Sage) [lucas_number1(n, 6, -7) for n in range(0, 21)] # Zerinvary Lajos, Apr 24 2009
(Magma) [Round(7^n/8): n in [0..30]]; // Vincenzo Librandi, Jun 24 2011
KEYWORD
nonn,easy
STATUS
approved
a(n) = 5*a(n-1) + 6*a(n-2), a(0) = 0, a(1) = 1.
+10
12
0, 1, 5, 31, 185, 1111, 6665, 39991, 239945, 1439671, 8638025, 51828151, 310968905, 1865813431, 11194880585, 67169283511, 403015701065, 2418094206391, 14508565238345, 87051391430071, 522308348580425
OFFSET
0,3
COMMENTS
Number of walks of length n between any two distinct vertices of the complete graph K_7. Example: a(2)=5 because the walks of length 2 between the vertices A and B of the complete graph ABCDEFG are ACB, ADB, AEB, AFB and AGB. - Emeric Deutsch, Apr 01 2004
Pisano period lengths: 1, 1, 2, 2, 2, 2, 14, 2, 2, 2, 10, 2, 12, 14, 2, 2, 16, 2, 18, 2, ... - R. J. Mathar, Aug 10 2012
Sum_{i=0..m} (-1)^(m+i)*6^i, for m >= 0, gives all terms after 0. - Bruno Berselli, Aug 28 2013
The ratio a(n+1)/a(n) converges to 6 as n approaches infinity. Also A053524, A080424, A051958. - Felix P. Muga II, Mar 09 2014
LINKS
Jean-Paul Allouche, Jeffrey Shallit, Zhixiong Wen, Wen Wu, Jiemeng Zhang, Sum-free sets generated by the period-k-folding sequences and some Sturmian sequences, arXiv:1911.01687 [math.CO], 2019.
Ji Young Choi, A Generalization of Collatz Functions and Jacobsthal Numbers, J. Int. Seq., Vol. 21 (2018), Article 18.5.4.
F. P. Muga II, Extending the Golden Ratio and the Binet-de Moivre Formula, March 2014; Preprint on ResearchGate.
FORMULA
a(n) = 5*a(n-1) + 6*a(n-2).
From Paul Barry, Apr 20 2003: (Start)
a(n) = (6^n - (-1)^n)/7.
G.f.: x/((1-6*x)*(1+x)).
E.g.f.: (exp(6*x) - exp(-x))/7. (End)
a(n) = 6^(n-1) - a(n-1). - Emeric Deutsch, Apr 01 2004
a(n+1) = Sum_{k=0..n} binomial(n-k, k)*5^(n-2*k)*6^k. - Paul Barry, Jul 29 2004
a(n) = round(6^n/7). - Mircea Merca, Dec 28 2010
a(n) = (-1)^(n-1)*Sum_{k=0..n-1} A135278(n-1,k)*(-7)^k) = (6^n - (-1)^n)/7 = (-1)^(n-1)*Sum_{k=0..n-1} (-6)^k. Equals (-1)^(n-1)*Phi(n,-6), where Phi is the cyclotomic polynomial when n is an odd prime. (For n > 0.) - Tom Copeland, Apr 14 2014
EXAMPLE
G.f. = x + 5*x^2 + 31*x^3 + 185*x^4 + 1111*x^5 + 6665*x^6 + 39991*x^7 + ...
MAPLE
seq(round(6^n/7), n=0..25); # Mircea Merca, Dec 28 2010
MATHEMATICA
k=0; lst={k}; Do[k = 6^n-k; AppendTo[lst, k], {n, 0, 5!}]; lst (* Vladimir Joseph Stephan Orlovsky, Dec 11 2008 *)
CoefficientList[Series[x / ((1 - 6 x) (1 + x)), {x, 0, 50}], x] (* Vincenzo Librandi, Mar 26 2014 *)
LinearRecurrence[{5, 6}, {0, 1}, 30] (* Harvey P. Dale, May 12 2015 *)
PROG
(Sage) [lucas_number1(n, 5, -6) for n in range(21)] # Zerinvary Lajos, Apr 24 2009
(Magma) [Floor(6^n/7-(-1)^n/7): n in [0..30]]; // Vincenzo Librandi, Jun 24 2011
(PARI) x='x+O('x^30); concat([0], Vec(x/((1-6*x)*(1+x)))) \\ G. C. Greubel, Jan 24 2018
(PARI) a(n) = round(6^n/7); \\ Altug Alkan, Sep 05 2018
CROSSREFS
Partial sums are in A033116. Cf. A014987.
KEYWORD
nonn,easy
STATUS
approved
Number of closed walks of length n on the complete graph on 6 nodes from a given node.
+10
12
1, 0, 5, 20, 105, 520, 2605, 13020, 65105, 325520, 1627605, 8138020, 40690105, 203450520, 1017252605, 5086263020, 25431315105, 127156575520, 635782877605, 3178914388020, 15894571940105, 79472859700520
OFFSET
0,3
LINKS
Ji Young Choi, A Generalization of Collatz Functions and Jacobsthal Numbers, J. Int. Seq., Vol. 21 (2018), Article 18.5.4.
Christopher R. Kitching, Henri Kauhanen, Jordan Abbott, Deepthi Gopal, Ricardo Bermúdez-Otero, and Tobias Galla, Estimating transmission noise on networks from stationary local order, arXiv:2405.12023 [cond-mat.stat-mech], 2024. See p. 48.
FORMULA
G.f.: (1 - 4*x)/(1 - 4*x - 5*x^2).
a(n) = (5^n + 5*(-1)^n)/6.
a(n) = 5^(n-1) - a(n-1), a(0) = 1. - Jon E. Schoenfield, Feb 08 2015
MATHEMATICA
k=0; lst={k}; Do[k=5^n-k; AppendTo[lst, k], {n, 1, 5!}]; lst (* Vladimir Joseph Stephan Orlovsky, Dec 11 2008 *)
CoefficientList[Series[(1 - 4*x)/(1 - 4*x - 5*x^2), {x, 0, 50}], x] (* or *) Table[(5^n + 5*(-1)^n)/6, {n, 0, 30}] (* G. C. Greubel, Dec 30 2017 *)
PROG
(PARI) for(n=0, 30, print1((5^n + 5*(-1)^n)/6, ", ")) \\ G. C. Greubel, Dec 30 2017
(Magma) [(5^n + 5*(-1)^n)/6: n in [0..30]]; // G. C. Greubel, Dec 30 2017
CROSSREFS
Cf. A109502.
Cf. sequences with the same recurrence form: A001045, A078008, A097073, A115341, A015518, A054878, A015521, A109499, A015531. - Vladimir Joseph Stephan Orlovsky, Dec 11 2008
KEYWORD
nonn,easy
AUTHOR
Mitch Harris, Jun 30 2005
EXTENSIONS
Corrected by Franklin T. Adams-Watters, Sep 18 2006
Edited by Jon E. Schoenfield, Feb 08 2015
STATUS
approved

Search completed in 0.025 seconds