Nothing Special   »   [go: up one dir, main page]

login
Search: a154105 -id:a154105
     Sort: relevance | references | number | modified | created      Format: long | short | data
Hex (or centered hexagonal) numbers: 3*n*(n+1)+1 (crystal ball sequence for hexagonal lattice).
(Formerly M4362)
+10
283
1, 7, 19, 37, 61, 91, 127, 169, 217, 271, 331, 397, 469, 547, 631, 721, 817, 919, 1027, 1141, 1261, 1387, 1519, 1657, 1801, 1951, 2107, 2269, 2437, 2611, 2791, 2977, 3169, 3367, 3571, 3781, 3997, 4219, 4447, 4681, 4921, 5167, 5419, 5677, 5941, 6211, 6487, 6769
OFFSET
0,2
COMMENTS
The hexagonal lattice is the familiar 2-dimensional lattice in which each point has 6 neighbors. This is sometimes called the triangular lattice.
Crystal ball sequence for A_2 lattice. - Michael Somos, Jun 03 2012
Sixth spoke of hexagonal spiral (cf. A056105-A056109).
Number of ordered integer triples (a,b,c), -n <= a,b,c <= n, such that a+b+c=0. - Benoit Cloitre, Jun 14 2003
Also the number of partitions of 6n into at most 3 parts, A001399(6n). - R. K. Guy, Oct 20 2003
Also, a(n) is the number of partitions of 6(n+1) into exactly 3 distinct parts. - William J. Keith, Jul 01 2004
Number of dots in a centered hexagonal figure with n+1 dots on each side.
Values of second Bessel polynomial y_2(n) (see A001498).
First differences of cubes (A000578). - Cecilia Rossiter (cecilia(AT)noticingnumbers.net), Dec 15 2004
Final digits of Hex numbers (hex(n) mod 10) are periodic with palindromic period of length 5 {1, 7, 9, 7, 1}. Last two digits of Hex numbers (hex(n) mod 100) are periodic with palindromic period of length 100. - Alexander Adamchuk, Aug 11 2006
All divisors of a(n) are congruent to 1, modulo 6. Proof: If p is an odd prime different from 3 then 3n^2 + 3n + 1 = 0 (mod p) implies 9(2n + 1)^2 = -3 (mod p), whence p = 1 (mod 6). - Nick Hobson, Nov 13 2006
For n>=1, a(n) is the side of Outer Napoleon Triangle whose reference triangle is a right triangle with legs (3a(n))^(1/2) and 3n(a(n))^(1/2). - Tom Schicker (tschicke(AT)email.smith.edu), Apr 25 2007
Number of triples (a,b,c) where 0<=(a,b)<=n and c=n (at least once the term n). E.g., for n = 1: (0,0,1), (0,1,0), (1,0,0), (0,1,1), (1,0,1), (1,1,0), (1,1,1), so a(1)=7. - Philippe Lallouet (philip.lallouet(AT)wanadoo.fr), Aug 20 2007
Equals the triangular numbers convolved with [1, 4, 1, 0, 0, 0, ...]. - Gary W. Adamson and Alexander R. Povolotsky, May 29 2009
From Terry Stickels, Dec 07 2009: (Start)
Also the maximum number of viewable cubes from any one static point while viewing a cube stack of identical cubes of varying magnitude.
For example, viewing a 2 X 2 X 2 stack will yield 7 maximum viewable cubes.
If the stack is 3 X 3 X 3, the maximum number of viewable cubes from any one static position is 19, and so on.
The number of cubes in the stack must always be the same number for width, length, height (at true regular cubic stack) and the maximum number of visible cubes can always be found by taking any cubic number and subtracting the number of the cube that is one less.
Examples: 125 - 64 = 61, 64 - 27 = 37, 27 - 8 = 19. (End)
The sequence of digital roots of the a(n) is period 3: repeat [1,7,1]. - Ant King, Jun 17 2012
The average of the first n (n>0) centered hexagonal numbers is the n-th square. - Philippe Deléham, Feb 04 2013
A002024 is the following array A read along antidiagonals:
1, 2, 3, 4, 5, 6, ...
2, 3, 4, 5, 6, 7, ...
3, 4, 5, 6, 7, 8, ...
4, 5, 6, 7, 8, 9, ...
5, 6, 7, 8, 9, 10, ...
6, 7, 8, 9, 10, 11, ...
and a(n) is the hook sum Sum_{k=0..n} A(n,k) + Sum_{r=0..n-1} A(r,n). - R. J. Mathar, Jun 30 2013
a(n) is the sum of the terms in the n+1 X n+1 matrices minus those in n X n matrices in an array formed by considering A158405 an array (the beginning terms in each row are 1,3,5,7,9,11,...). - J. M. Bergot, Jul 05 2013
The formula also equals the product of the three distinct combinations of two consecutive numbers: n^2, (n+1)^2, and n*(n+1). - J. M. Bergot, Mar 28 2014
The sides of any triangle ABC are divided into 2n + 1 equal segments by 2n points: A_1, A_2, ..., A_2n in side a, and also on the sides b and c cyclically. If A'B'C' is the triangle delimited by AA_n, BB_n and CC_n cevians, we have (ABC)/(A'B'C') = a(n) (see Java applet link). - Ignacio Larrosa Cañestro, Jan 02 2015
a(n) is the maximal number of parts into which (n+1) triangles can intersect one another. - Ivan N. Ianakiev, Feb 18 2015
((2^m-1)n)^t mod a(n) = ((2^m-1)(n+1))^t mod a(n) = ((2^m-1)(2n+1))^t mod a(n), where m any positive integer, and t = 0(mod 6). - Alzhekeyev Ascar M, Oct 07 2016
((2^m-1)n)^t mod a(n) = ((2^m-1)(n+1))^t mod a(n) = a(n) - (((2^m-1)(2n+1))^t mod a(n)), where m any positive integer, and t = 3(mod 6). - Alzhekeyev Ascar M, Oct 07 2016
(3n+1)^(a(n)-1) mod a(n) = (3n+2)^(a(n)-1) mod a(n) = 1. If a(n) not prime, then always strong pseudoprime. - Alzhekeyev Ascar M, Oct 07 2016
Every positive integer is the sum of 8 hex numbers (zero included), at most 3 of which are greater than 1. - Mauro Fiorentini, Jan 01 2018
Area enclosed by the segment of Archimedean spiral between n*Pi/2 and (n+1)*Pi/2 in Pi^3/48 units. - Carmine Suriano, Apr 10 2018
This sequence contains all numbers k such that 12*k - 3 is a square. - Klaus Purath, Oct 19 2021
The continued fraction expansion of sqrt(3*a(n)) is [3n+1; {1, 1, 2n, 1, 1, 6n+2}]. For n = 0, this collapses to [1; {1, 2}]. - Magus K. Chu, Sep 12 2022
REFERENCES
M. Gardner, Time Travel and Other Mathematical Bewilderments. Freeman, NY, 1988, p. 18.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
LINKS
G. L. Alexanderson and John E. Wetzel, Dissections of a tetrahedron, J. Combinatorial Theory Ser. B 11 (1971), 58--66. MR0303412 (46 #2549). See p. 58.
B. T. Bennett and R. B. Potts, Arrays and brooks, J. Austral. Math. Soc., 7 (1967), 23-31 (see p. 30).
B. T. Bennett and R. B. Potts, Arrays and brooks, J. Austral. Math. Soc., 7 (1967), 23-31. [Annotated scanned copy]
Aran Bingham, Commutative n-ary Arithmetic, University of New Orleans Theses and Dissertations, Paper 1959, 2015.
J. H. Conway and N. J. A. Sloane, Low-Dimensional Lattices VII: Coordination Sequences, Proc. Royal Soc. London, A453 (1997), 2369-2389 (pdf).
M. Gardner & N. J. A. Sloane, Correspondence, 1973-74
R. K. Guy, The strong law of small numbers. Amer. Math. Monthly 95 (1988), no. 8, 697-712.
R. K. Guy, The strong law of small numbers. Amer. Math. Monthly 95 (1988), no. 8, 697-712. [Annotated scanned copy]
Md. Towhidul Islam, Extending triangle
G. S. Kazandzidis, On a Conjecture of Moessner and a General Problem, Bull. Soc. Math. Grèce, Nouvelle Série - vol. 2, fasc. 1-2, pp. 23-30 (1961).
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992
B. K. Teo and N. J. A. Sloane, Magic numbers in polygonal and polyhedral clusters, Inorgan. Chem. 24 (1985), 4545-4558.
Eric Weisstein's World of Mathematics, Hex Number
Eric Weisstein's World of Mathematics, Nexus Number
Eric Weisstein's World of Mathematics, Outer Napoleon Triangle.
FORMULA
a(n) = 3*n*(n+1) + 1, n >= 0 (see the name).
a(n) = (n+1)^3 - n^3 = a(-1-n).
G.f.: (1 + 4*x + x^2) / (1 - x)^3. - Simon Plouffe in his 1992 dissertation
a(n) = 6*A000217(n) + 1.
a(n) = a(n-1) + 6*n = 2a(n-1) - a(n-2) + 6 = 3*a(n-1) - 3*a(n-2) + a(n-3) = A056105(n) + 5n = A056106(n) + 4*n = A056107(n) + 3*n = A056108(n) + 2*n = A056108(n) + n.
n-th partial arithmetic mean is n^2. - Amarnath Murthy, May 27 2003
a(n) = 1 + Sum_{j=0..n} (6*j). E.g., a(2)=19 because 1+ 6*0 + 6*1 + 6*2 = 19. - Xavier Acloque, Oct 06 2003
The sum of the first n hexagonal numbers is n^3. That is, Sum_{n>=1} (3*n*(n-1) + 1) = n^3. - Edward Weed (eweed(AT)gdrs.com), Oct 23 2003
a(n) = right term in M^n * [1 1 1], where M = the 3 X 3 matrix [1 0 0 / 2 1 0 / 3 3 1]. M^n * [1 1 1] = [1 2n+1 a(n)]. E.g., a(4) = 61, right term in M^4 * [1 1 1], since M^4 * [1 1 1] = [1 9 61] = [1 2n+1 a(4)]. - Gary W. Adamson, Dec 22 2004
Row sums of triangle A130298. - Gary W. Adamson, Jun 07 2007
a(n) = 3*n^2 + 3*n + 1. Proof: 1) If n occurs once, it may be in 3 positions; for the two other ones, n terms are independently possible, then we have 3*n^2 different triples. 2) If the term n occurs twice, the third one may be placed in 3 positions and have n possible values, then we have 3*n more different triples. 3) The term n may occurs 3 times in one way only that gives the formula. - Philippe Lallouet (philip.lallouet(AT)wanadoo.fr), Aug 20 2007
Binomial transform of [1, 6, 6, 0, 0, 0, ...]; Narayana transform (A001263) of [1, 6, 0, 0, 0, ...]. - Gary W. Adamson, Dec 29 2007
a(n) = (n-1)*A000166(n) + (n-2)*A000166(n-1) = (n-1)floor(n!*e^(-1)+1) + (n-2)*floor((n-1)!*e^(-1)+1) (with offset 0). - Gary Detlefs, Dec 06 2009
a(n) = A028896(n) + 1. - Omar E. Pol, Oct 03 2011
a(n) = integral( (sin((n+1/2)x)/sin(x/2))^3, x=0..Pi)/Pi. - Yalcin Aktar, Dec 03 2011
Sum_{n>=0} 1/a(n) = Pi/sqrt(3)*tanh(Pi/(2*sqrt(3))) = 1.305284153013581... - Ant King, Jun 17 2012
a(n) = A000290(n) + A000217(2n+1). - Ivan N. Ianakiev, Sep 24 2013
a(n) = A002378(n+1) + A056220(n) = A005408(n) + 2*A005449(n) = 6*A000217(n) + 1. - Ivan N. Ianakiev, Sep 26 2013
a(n) = 6*A000124(n) - 5. - Ivan N. Ianakiev, Oct 13 2013
a(n) = A239426(n+1) / A239449(n+1) = A215630(2*n+1,n+1). - Reinhard Zumkeller, Mar 19 2014
a(n) = A243201(n) / A002061(n + 1). - Mathew Englander, Jun 03 2014
a(n) = A101321(6,n). - R. J. Mathar, Jul 28 2016
E.g.f.: (1 + 6*x + 3*x^2)*exp(x). - Ilya Gutkovskiy, Jul 28 2016
a(n) = (A001844(n) + A016754(n))/2. - Bruce J. Nicholson, Aug 06 2017
a(n) = A045943(2n+1). - Miquel Cerda, Jan 22 2018
a(n) = 3*Integral_{x=n..n+1} x^2 dx. - Carmine Suriano, Apr 10 2018
a(n) = A287326(A000124(n), 1). - Kolosov Petro, Oct 22 2018
From Amiram Eldar, Jun 20 2020: (Start)
Sum_{n>=0} a(n)/n! = 10*e.
Sum_{n>=0} (-1)^(n+1)*a(n)/n! = 2/e. (End)
G.f.: polylog(-3, x)*(1-x)/x. See the Simon Plouffe formula above, and the g.f. of the rows of A008292 by Vladeta Jovovic, Sep 02 2002. - Wolfdieter Lang, May 08 2021
a(n) = T(n-1)^2 - 2*T(n)^2 + T(n+1)^2, n >= 1, T = triangular number A000217. - Klaus Purath, Oct 11 2021
a(n) = 1 + 2*Sum_{j=n..2n} j. - Klaus Purath, Oct 19 2021
a(n) = A069099(n+1) - A000217(n). - Klaus Purath, Nov 03 2021
From Leo Tavares, Dec 03 2021: (Start)
a(n) = A005448(n) + A140091(n);
a(n) = A001844(n) + A002378(n);
a(n) = A005891(n) + A000217(n);
a(n) = A000290(n) + A000384(n+1);
a(n) = A060544(n-1) + 3*A000217(n);
a(n) = A060544(n-1) + A045943(n).
a(2*n+1) = A154105(n).
(End)
EXAMPLE
G.f. = 1 + 7*x + 19*x^2 + 37*x^3 + 61*x^4 + 91*x^5 + 127*x^6 + 169*x^7 + 217*x^8 + ...
From Omar E. Pol, Aug 21 2011: (Start)
Illustration of initial terms:
.
. o o o o
. o o o o o o o o
. o o o o o o o o o o o o
. o o o o o o o o o o o o o o o o
. o o o o o o o o o o o o
. o o o o o o o o
. o o o o
.
. 1 7 19 37
.
(End)
From Klaus Purath, Dec 03 2021: (Start)
(1) a(19) is not a prime number, because besides a(19) = a(9) + P(29), a(19) = a(15) + P(20) = a(2) + P(33) is also true.
(2) a(25) is prime, because except for a(25) = a(12) + P(38) there is no other equation of this pattern. (End)
MAPLE
A003215:=n->3*n*(n+1)+1; seq(A003215(n), n=0..100); # Wesley Ivan Hurt, Mar 28 2014
MATHEMATICA
FoldList[#1 + #2 &, 1, 6 Range@ 50] (* Robert G. Wilson v, Feb 02 2011 *)
LinearRecurrence[{3, -3, 1}, {1, 7, 19}, 47] (* Robert G. Wilson v, Jul 06 2013 *)
PROG
(PARI) {a(n) = 3*n*(n+1) + 1};
(Haskell)
a003215 n = 3 * n * (n + 1) + 1 -- Reinhard Zumkeller, Oct 22 2011
(Maxima) makelist(3*n*(n+1)+1, n, 0, 30); /* Martin Ettl, Nov 12 2012 */
(Magma) [3*n*(n+1)+1: n in [0..50]]; // G. C. Greubel, Nov 04 2017
(Python) [3*n*(n+1)+1 for n in range(47)] # Michael S. Branicky, Jan 07 2021
CROSSREFS
Column k=3 of A080853, and column k=2 of A047969.
See also A220083 for a list of numbers of the form n*P(s,n)-(n-1)*P(s,n-1), where P(s,n) is the n-th polygonal number with s sides.
Cf. A287326(A000124(n), 1).
Cf. A008292.
Cf. A154105.
KEYWORD
nonn,easy,nice,changed
EXTENSIONS
Partially edited by Joerg Arndt, Mar 11 2010
STATUS
approved
a(n) = 2*n*(4*n - 1).
+10
20
0, 6, 28, 66, 120, 190, 276, 378, 496, 630, 780, 946, 1128, 1326, 1540, 1770, 2016, 2278, 2556, 2850, 3160, 3486, 3828, 4186, 4560, 4950, 5356, 5778, 6216, 6670, 7140, 7626, 8128, 8646, 9180, 9730, 10296, 10878, 11476, 12090, 12720, 13366, 14028, 14706
OFFSET
0,2
COMMENTS
Even hexagonal numbers.
Number of edges in the join of two complete graphs of order 3n and n, K_3n * K_n - Roberto E. Martinez II, Jan 07 2002
Bisection of A000384. Also, this sequence arises from reading the line from 0, in the direction 0, 6, ..., in the square spiral whose vertices are the triangular numbers A000217. Perfect numbers are members of this sequence because a(A134708(n)) = A000396(n). Also, positive members are a bisection of A139596. - Omar E. Pol, May 07 2008
FORMULA
a(n) = C(4*n,2), n>=0. - Zerinvary Lajos, Jan 02 2007
O.g.f.: 2*x*(3+5*x)/(1-x)^3. - R. J. Mathar, May 06 2008
a(n) = 8*n^2 - 2*n. - Omar E. Pol, May 07 2008
a(n) = a(n-1) + 16*n - 10 (with a(0)=0). - Vincenzo Librandi, Nov 20 2010
E.g.f.: (8*x^2 + 6*x)*exp(x). - G. C. Greubel, Jul 18 2017
From Vaclav Kotesovec, Aug 18 2018: (Start)
Sum_{n>=1} 1/a(n) = 3*log(2)/2 - Pi/4.
Sum_{n>=1} (-1)^n / a(n) = log(2)/2 + log(1+sqrt(2))/sqrt(2) - Pi / 2^(3/2). (End)
a(n) = A154105(n-1) - A016754(n-1). - Leo Tavares, May 02 2023
MAPLE
[seq(binomial(4*n, 2), n=0..43)]; # Zerinvary Lajos, Jan 02 2007
MATHEMATICA
s=0; lst={s}; Do[s+=n++ +6; AppendTo[lst, s], {n, 0, 7!, 16}]; lst (* Vladimir Joseph Stephan Orlovsky, Nov 16 2008 *)
Table[2*n*(4*n - 1), {n, 0, 50}] (* G. C. Greubel, Jul 18 2017 *)
PolygonalNumber[6, Range[0, 90, 2]] (* or *) LinearRecurrence[{3, -3, 1}, {0, 6, 28}, 50] (* Harvey P. Dale, Jan 21 2023 *)
PROG
(Magma) [2*n*(4*n-1): n in [0..50]]; // Vincenzo Librandi, Apr 25 2011
(PARI) a(n)=2*n*(4*n-1) \\ Charles R Greathouse IV, Jun 17 2017
CROSSREFS
KEYWORD
nonn,easy
EXTENSIONS
More terms from Erich Friedman
STATUS
approved
a(n) = n^3 + sum((-1)^j*a(j)); for j=1 to n-1; a(1)=1.
+10
4
1, 7, 33, 37, 135, 91, 309, 169, 555, 271, 873, 397, 1263, 547, 1725, 721, 2259, 919, 2865, 1141, 3543, 1387, 4293, 1657, 5115, 1951, 6009, 2269, 6975, 2611, 8013, 2977, 9123, 3367, 10305, 3781, 11559, 4219, 12885, 4681, 14283, 5167, 15753, 5677, 17295
OFFSET
1,2
COMMENTS
1 followed by interleaving of A154105 and 3*A154106. - Klaus Brockhaus, Jan 04 2009
FORMULA
G.f.: x*(1 + 7*x + 30*x^2 + 16*x^3 + 39*x^4 + x^5 + 2*x^6)/((1+x)^3*(1-x)^3). - Klaus Brockhaus, Jan 04 2009
From Walter Carlini, Jan 12 2009: (Start)
a(n) = 3n^2 - 3n + 1 if n is 1 or an even number;
a(n) = 9n^2 - 21n + 15 if n is any odd number other than 1. (End)
EXAMPLE
a(1)=1, a(2)=2^3-a(1)=8-1=7, a(3)=3^3+a(2)-a(1)=27+7-1=33, a(4)=64-33+7-1=37, a(5)=125+37-33+7-1=135, a(6)=216-135+37-33+7-1=91, etc.
PROG
(Magma) S:=[ 1 ]; for n in [2..45] do Append(~S, n^3 + &+[ (-1)^j*S[j]: j in [1..n-1] ]); end for; S; // Klaus Brockhaus, Jan 04 2009
CROSSREFS
The third of a family of sequences that includes A153284 and A153285.
Cf. A154105 (12*n^2 + 18*n + 7), A154106 (12*n^2 + 22*n + 11). - Klaus Brockhaus, Jan 04 2009
KEYWORD
easy,nonn
AUTHOR
Walter Carlini, Dec 23 2008, Jan 01 2009
EXTENSIONS
Extended beyond a(30) by Klaus Brockhaus, Jan 04 2009
G.f. corrected by Klaus Brockhaus, Oct 15 2009
STATUS
approved
a(n) = 3*n^4 + 6*n^3 - 3*n + 1.
+10
2
1, 7, 91, 397, 1141, 2611, 5167, 9241, 15337, 24031, 35971, 51877, 72541, 98827, 131671, 172081, 221137, 279991, 349867, 432061, 527941, 638947, 766591, 912457, 1078201, 1265551, 1476307, 1712341, 1975597, 2268091, 2591911, 2949217, 3342241, 3773287, 4244731
OFFSET
0,2
COMMENTS
If gcd(n,7) = gcd(n+1,7) = gcd(2*n+1,7) = 1 then a(n) == 0 (mod 7) (E. Picutti, see References).
REFERENCES
Ettore Picutti, Sul numero e la sua storia, Feltrinelli Economica, 1977, p. 208.
FORMULA
G.f.: (1 + 2*x + 66*x^2 + 2*x^3 + x^4)/(1-x)^5.
a(n) = a(-n-1) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5).
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4) + 6*12.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + 6*A008594(n-1).
a(n) = 2*a(n-1) - a(n-2) + 6*A003154(n).
a(n) = a(n-1) + 6*A007588(n).
a(n) = 1 + 6*A062392(n).
a(n) = 7*A000540(n)/A000330(n) = A154105(A000096(n-1)) for n > 0.
Sum_{i=0..n} a(i) = (3*n^5 + 15*n^4 + 20*n^3 - 3*n + 5)/5.
a(n) = 7*(3*n^2 + 3*n - 1)*(Sum_{k=1..n} k^6)/(5*Sum_{k=1..n} k^4), n > 0. - Gary Detlefs, Oct 18 2011
MATHEMATICA
Table[3 n^4 + 6 n^3 - 3 n + 1, {n, 0, 40}] (* Vincenzo Librandi, Mar 26 2013 *)
LinearRecurrence[{5, -10, 10, -5, 1}, {1, 7, 91, 397, 1141}, 40] (* Harvey P. Dale, Jul 12 2022 *)
PROG
(Magma) [3*n^4+6*n^3-3*n+1: n in [0..31]];
CROSSREFS
Subsequence of A003215.
KEYWORD
nonn,easy
AUTHOR
Bruno Berselli, Oct 25 2010 - Oct 29 2010
EXTENSIONS
Formula, program and crossref added by Bruno Berselli, Aug 22 2011
STATUS
approved
a(n) = 12*n^2 + 18*n.
+10
1
0, 30, 84, 162, 264, 390, 540, 714, 912, 1134, 1380, 1650, 1944, 2262, 2604, 2970, 3360, 3774, 4212, 4674, 5160, 5670, 6204, 6762, 7344, 7950, 8580, 9234, 9912, 10614, 11340, 12090, 12864, 13662, 14484, 15330, 16200, 17094, 18012, 18954, 19920
OFFSET
0,2
COMMENTS
For n>=3, a(n) is the second Zagreb index of the double-wheel graph DW[n]. The second Zagreb index of a simple connected graph g is the sum of the degree products d(i) d(j) over all edges ij of g.
The double-wheel graph DW[n] consists of two cycles C[n], whose vertices are connected to an additional vertex.
The M-polynomial of the double-wheel graph DW[n] is M(DW[n],x,y) = 2*n*x^3*y^3 + 2*n*x^3*y^{2*n}.
LINKS
E. Deutsch and Sandi Klavzar, M-polynomial and degree-based topological indices, Iranian J. Math. Chemistry, 6, No. 2, 2015, 93-102.
FORMULA
G.f.: 6*x*(5-x)/(1-x)^3.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3).
a(n) = 6*A014106(n).
a(n) = A152746(n+1) - 6 = A154105(n) - 7. - Omar E. Pol, May 08 2018
EXAMPLE
a(3) = 162. Indeed, the double-wheel graph DW[3] has 6 edges with end-point degrees 3,3 and 6 edges with end-point degrees 3,6. Then the second Zagreb index is 6*9 + 6*18 = 162.
MAPLE
seq(12*n^2+18*n, n = 0 .. 50);
MATHEMATICA
Table[12 n^2 + 18 n, {n, 0, 45}] (* Vincenzo Librandi, Nov 09 2016 *)
PROG
(Magma) [12*n^2+18*n: n in [0..40]]; // Vincenzo Librandi, Nov 09 2016
(PARI) a(n)=12*n^2+18*n \\ Charles R Greathouse IV, Nov 09 2016
CROSSREFS
First bisection of A277978.
After 0, subsequence of A255265.
KEYWORD
nonn,easy
AUTHOR
Emeric Deutsch, Nov 08 2016
STATUS
approved

Search completed in 0.009 seconds