A176398 -id:A176398 - cp's OEIS frontend

A261327 a(n) = (n^2 + 4) / 4^((n + 1) mod 2).

1, 5, 2, 13, 5, 29, 10, 53, 17, 85, 26, 125, 37, 173, 50, 229, 65, 293, 82, 365, 101, 445, 122, 533, 145, 629, 170, 733, 197, 845, 226, 965, 257, 1093, 290, 1229, 325, 1373, 362, 1525, 401, 1685, 442, 1853, 485, 2029, 530, 2213, 577, 2405, 626, 2605, 677

Offset: 0

Paul Curtz, Aug 15 2015

Using (n+sqrt(4+n^2))/2, after the integer 1 for n=0, the reduced metallic means are b(1) = (1+sqrt(5))/2, b(2) = 1+sqrt(2), b(3) = (3+sqrt(13))/2, b(4) = 2+sqrt(5), b(5) = (5+sqrt(29))/2, b(6) = 3+sqrt(10), b(7) = (7+sqrt(53))/2, b(8) = 4+sqrt(17), b(9) = (9+sqrt(85))/2, b(10) = 5+sqrt(26), b(11) = (11+sqrt(125))/2 = (11+5*sqrt(5))/2, ... . The last value yields the radicals in a(n) or A013946.
b(2) = 2.41, b(3) = 3.30, b(4) = 4.24, b(5) = 5.19 are "good" approximations of fractal dimensions corresponding to dimensions 3, 4, 5, 6: 2.48, 3.38, 4.33 and 5.45 based on models. See "Arbres DLA dans les espaces de dimension supérieure: la théorie des peaux entropiques" in Queiros-Condé et al. link. DLA: beginning of the title of the Witten et al. link.
Consider the symmetric array of the half extended Rydberg-Ritz spectrum of the hydrogen atom:
   0,    1/0,     1/0,     1/0,     1/0,      1/0,      1/0,     1/0, ...
-1/0,      0,     3/4,     8/9,   15/16,    24/25,    35/36,   48/49, ...
-1/0,   -3/4,       0,    5/36,    3/16,   21/100,      2/9,  45/196, ...
-1/0,   -8/9,   -5/36,       0,   7/144,   16/225,     1/12,  40/441, ...
-1/0, -15/16,   -3/16,  -7/144,       0,    9/400,    5/144,  33/784, ...
-1/0, -24/25, -21/100, -16/225,  -9/400,        0,   11/900, 24/1225, ...
-1/0, -35/36,    -2/9,   -1/12,  -5/144,  -11/900,        0, 13/1764, ...
-1/0, -48/49, -45/196, -40/441, -33/784, -24/1225, -13/1764,       0, ... .
The numerators are almost A165795(n).
Successive rows: A000007(n)/A057427(n), A005563(n-1)/A000290(n), A061037(n)/A061038(n), A061039(n)/A061040(n), A061041(n)/A061042(n), A061043(n)/A061044(n), A061045(n)/A061046(n), A061047(n)/A061048(n), A061049(n)/A061050(n).
A144433(n) or A195161(n+1) are the numerators of the second upper diagonal (denominators: A171522(n)).
c(n+1) = a(n) + a(n+1) = 6, 7, 15, 18, 34, 39, 63, 70, 102, 111, ... .
c(n+3) - c(n+1) = 9, 11, 19, 21, 29, 31, ... = A090771(n+2).
The final digit of a(n) is neither 4 nor 8. - Paul Curtz, Jan 30 2019

Colin Barker, Table of n, a(n) for n = 0..1000
Diogo Queiros-Condé, Jean Chaline, and Jacques Dubois, Le monde des fractales La Nature trans-échelles, 478p., ellipses, Paris, 2015, page 220.
T. A. Witten, Jr. and L. M. Sander, Diffusion-Limited Aggregation, a Kinetic Critical Phenomenom, Phys. Rev. Lett., Vol. 47 (Nov 09 1981), pp. 1400-1403.
Index entries for linear recurrences with constant coefficients, signature (0,3,0,-3,0,1).

Cf. A000007, A000290, A001622, A002522, A005563, A010685, A010698, A013946, A014176, A057427, A061035-A061050, A078370, A087475, A090771, A098316, A098317, A098318, A144433, A168077, A171621, A176398, A176439, A176458, A176522, A176537, A195161.

Cf. A069894.

[Numerator(1+n^2/4): n in [0..60]]; // Vincenzo Librandi, Aug 15 2015

A261327:=n->numer((4 + n^2)/4); seq(A261327(n), n=0..60); # Wesley Ivan Hurt, Aug 15 2015

LinearRecurrence[{0, 3, 0, -3, 0, 1}, {1, 5, 2, 13, 5, 29}, 60] (* Vincenzo Librandi, Aug 15 2015 *)
a[n_] := (n^2 + 4) / 4^Mod[n + 1, 2]; Table[a[n], {n, 0, 52}] (* Peter Luschny, Mar 18 2022 *)

vector(60, n, n--; numerator(1+n^2/4)) \\ Michel Marcus, Aug 15 2015

Vec((1+5*x-x^2-2*x^3+2*x^4+5*x^5)/(1-x^2)^3 + O(x^60)) \\ Colin Barker, Aug 15 2015

a(n)=if(n%2,n^2+4,(n/2)^2+1) \\ Charles R Greathouse IV, Oct 16 2015

[(n*n+4)//4**((n+1)%2) for n in range(60)] # Gennady Eremin, Mar 18 2022

[numerator(1+n^2/4) for n in (0..60)] # G. C. Greubel, Feb 09 2019

a(n) = numerator(1 + n^2/4). (Previous name.) See A010685 (denominators).
a(2*k) = 1 + k^2.
a(2*k+1) = 5 + 4*k*(k+1).
a(2*k+1) = 4*a(2*k) + 4*k + 1.
a(4*k+2) = A069894(k). - Paul Curtz, Jan 30 2019
a(-n) = a(n).
a(n+2) = a(n) + A144433(n) (or A195161(n+1)).
a(n) = A168077(n) + period 2: repeat 1, 4.
a(n) = A171621(n) + period 2: repeat 2, 8.
From Colin Barker, Aug 15 2015: (Start)
a(n) = (5 - 3*(-1)^n)*(4 + n^2)/8.
a(n) = n^2/4 + 1 for n even;
a(n) = n^2 + 4 for n odd.
a(n) = 3*a(n-2) - 3*a(n-4) + a(n-6) for n>5.
G.f.: (1 + 5*x - x^2 - 2*x^3 + 2*x^4 + 5*x^5)/ (1 - x^2)^3. (End)
E.g.f.: (5/8)*(x^2 + x + 4)*exp(x) - (3/8)*(x^2 - x + 4)*exp(-x). - Robert Israel, Aug 18 2015
Sum_{n>=0} 1/a(n) = (4*coth(Pi)+tanh(Pi))*Pi/8 + 1/2. - Amiram Eldar, Mar 22 2022

New name by Peter Luschny, Mar 18 2022

A221874 Numbers m such that 10*m^2 + 6 is a square.

Original entry on oeis.org

1, 5, 43, 191, 1633, 7253, 62011, 275423, 2354785, 10458821, 89419819, 397159775, 3395598337, 15081612629, 128943316987, 572704120127, 4896450447169, 21747674952197, 185936173675435, 825838944063359, 7060678149219361, 31360132199455445

Offset: 1

Bruno Berselli, Jan 28 2013

The Diophantine equation 10*x^2 + k = y^2, for |k| < 10, has integer solutions with the following k values:
k =  1, the nonnegative x values are in A084070;
k = -1,               "                 A097315;
k =  4,               "                 2*A084070;
k = -4,               "                 2*A097315;
k =  6,               "                 this sequence;
k = -6,               "                 A221875;
k =  9,               "                 A075836;
k = -9,               "                 A052454.
a(n+1)/a(n) tends alternately to (sqrt(2)+sqrt(5))^2/3 and (2*sqrt(2)+sqrt(5))^2/3; a(n+2)/a(n) tends to A176398^2.

Vincenzo Librandi, Table of n, a(n) for n = 1..1000
Index entries for linear recurrences with constant coefficients, signature (0,38,0,-1).

Cf. A097315, A084070, A075836, A052454, A129556, A221875.

Subsequence of A031150.

m:=22; R:=PowerSeriesRing(Integers(), m); Coefficients(R!((1+x)*(1+4*x+x^2)/((1-6*x-x^2)*(1+6*x-x^2))));

A221874:=proc(q)
local n;
for n from 1 to q do if type(sqrt(10*n^2+6),integer) then print(n);
fi; od; end:
A221874(100000000000000000); # Paolo P. Lava, Feb 11 2013

LinearRecurrence[{0, 38, 0, -1}, {1, 5, 43, 191}, 22]

makelist(expand(((-5*(-1)^n+2*sqrt(10))*(3+sqrt(10))^(2*floor(n/2))-(5*(-1)^n+2*sqrt(10))*(3-sqrt(10))^(2*floor(n/2)))/10), n, 1, 22);

G.f.: x*(1+x)*(1+4*x+x^2)/((1-6*x-x^2)*(1+6*x-x^2)).
a(n) = ((-5*(-1)^n+2*t)*(3+t)^(2*floor(n/2)) - (5*(-1)^n+2*t)*(3-t)^(2*floor(n/2)))/10, where t=sqrt(10).
a(n) = 2*A129556(n) + 1.
a(n)*a(n-3) - a(n-1)*a(n-2) = -36 + 12(-1)^n.

A179237 a(0) = 1, a(1) = 2; a(n+1) = 6*a(n) + a(n-1) for n>1.

Original entry on oeis.org

1, 2, 13, 80, 493, 3038, 18721, 115364, 710905, 4380794, 26995669, 166354808, 1025124517, 6317101910, 38927735977, 239883517772, 1478228842609, 9109256573426, 56133768283165, 345911866272416, 2131604965917661, 13135541661778382, 80944854936587953

Offset: 0

Gary W. Adamson, Jul 04 2010

a(n)/a(n-1) converges to 1/(sqrt(10) - 3) = 6.16227766017... = A176398.

			a(5) = 3038 = 6*a(5) + a(4) = 6*493 + 80.
a(5) = term (1,1) in M^5 where M^5 = [3038, 4215, 4215, 5848].

Colin Barker, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (6,1).

I:=[1,2]; [n le 2 select I[n] else 6*Self(n-1)+Self(n-2): n in [1..40]]; // Vincenzo Librandi, Oct 13 2015

CoefficientList[Series[(-1 + 4 x)/(-1 + 6 x + x^2), {x, 0, 33}], x] (* Vincenzo Librandi, Oct 13 2015 *)

Vec((-1+4*x)/(-1+6*x+x^2) + O(x^40)) \\ Colin Barker, Oct 13 2015

Let M = the 2x2 matrix [2,3; 3,4]. a(n) = term (1,1) in M^n.
G.f.: ( -1+4*x ) / ( -1+6*x+x^2 ). a(n) = A005668(n) + A015451(n). - R. J. Mathar, Jul 06 2012
a(n) = ((3-sqrt(10))^n*(1+sqrt(10))+(-1+sqrt(10))*(3+sqrt(10))^n)/(2*sqrt(10)). - Colin Barker, Oct 13 2015
a(n) = Sum_{k=0..n-2} A168561(n-2,k)*6^k + 2 * Sum_{k=0..n-1} A168561(n-1,k)*6^k, n>0. - R. J. Mathar, Feb 14 2024
a(n) = A005668(n+1) - 4*A005668(n). - R. J. Mathar, Feb 14 2024

Corrected by R. J. Mathar, Jul 06 2012

A244593 Decimal expansion of z_c = phi^5 (where phi is the golden ratio), a lattice statistics constant which is the exact value of the critical activity of the hard hexagon model.

Original entry on oeis.org

1, 1, 0, 9, 0, 1, 6, 9, 9, 4, 3, 7, 4, 9, 4, 7, 4, 2, 4, 1, 0, 2, 2, 9, 3, 4, 1, 7, 1, 8, 2, 8, 1, 9, 0, 5, 8, 8, 6, 0, 1, 5, 4, 5, 8, 9, 9, 0, 2, 8, 8, 1, 4, 3, 1, 0, 6, 7, 7, 2, 4, 3, 1, 1, 3, 5, 2, 6, 3, 0, 2, 3, 1, 4, 0, 9, 4, 5, 1, 2, 2, 4, 8, 5, 3, 6, 0, 3, 6, 0, 2, 0, 9, 4, 6, 9, 5, 5, 6, 8, 7, 4, 2

Offset: 2

Jean-François Alcover, Jul 01 2014

Essentially the same digit sequence as A239798, A019863 and A019827. - R. J. Mathar, Jul 03 2014

The minimal polynomial of this constant is x^2 - 11*x - 1. - Joerg Arndt, Jan 01 2017

			11.09016994374947424102293417182819058860154589902881431067724311352630...

Steven R. Finch, Mathematical Constants, Cambridge University Press, 2003, Section 5.12.1 Phase transitions in Lattice Gas Models, p. 347.
Alfred S. Posamentier, Math Charmers, Tantalizing Tidbits for the Mind, Prometheus Books, NY, 2003, pages 138-139.
David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See p. 83.

Eric Weisstein's MathWorld, Hard Hexagon Entropy Constant.
Wikipedia, Hard Hexagon Model.
D. W. Wood and R. W. Turnbull, z^2-11z-1 as an algebraic invariant for the hard-hexagon model, 1988 J. Phys. A: Math. Gen. 21 L989.

Cf. A019863, A019827, A085850, A239798.

Cf. A001622, A014176, A098316, A098317, A098318, A176398, A176439, A176458 A176522, A176537 (metallic means).

RealDigits[GoldenRatio^5, 10, 103] // First

(5*sqrt(5)+11)/2 \\ Charles R Greathouse IV, Aug 10 2016

Equals ((1 + sqrt(5))/2)^5 = (11 + 5*sqrt(5))/2.
Equals phi^5 = 11 + 1/phi^5 = 3 + 5*phi, an integer in the quadratic number field Q(sqrt(5)). - Wolfdieter Lang, Nov 11 2023
Equals lim_{n->infinity} S(n, 5*(-1 + 2*phi))/ S(n-1, 5*(-1 + 2*phi)), with the S-Chebyshev polynomials (see A049310). - Wolfdieter Lang, Nov 15 2023

A222390 Nonnegative integers m such that 10m(m+1)+1 is a square.

Original entry on oeis.org

0, 3, 15, 132, 588, 5031, 22347, 191064, 848616, 7255419, 32225079, 275514876, 1223704404, 10462309887, 46468542291, 397292260848, 1764580902672, 15086643602355, 67007605759263, 572895164628660, 2544524437949340, 21754929612286743, 96624921036315675

Offset: 1

Bruno Berselli, Feb 18 2013

a(n+1)/a(n) tends alternately to (7+2*sqrt(10))/3 and (13+4*sqrt(10))/3; a(n+2)/a(n) tends to A176398^2.

Subsequence of A014601.

Bruno Berselli, Table of n, a(n) for n = 1..500
Index entries for linear recurrences with constant coefficients, signature (1,38,-38,-1,1).

Cf. nonnegative integers m such that k*m*(m+1)+1 is a square: A001652 (k=2), A001921 (k=3), A001477 (k=4), A053606 (k=5), A105038 (k=6), A105040 (k=7), A053141 (k=8), this sequence (k=10), A105838 (k=11), A061278 (k=12), A104240 (k=13); A105063 (k=17), A222393 (k=18), A101180 (k=19), A077259 (k=20) [incomplete list].

Cf. A221875.

m:=22; R:=PowerSeriesRing(Integers(), m); [0] cat Coefficients(R!(3*(1+4*x+x^2)/((1-x)*(1-6*x-x^2)*(1+6*x-x^2))));

I:=[0,3,15,132,588]; [n le 5 select I[n] else Self(n-1) +38*Self(n-2)-38*Self(n-3)-Self(n-4)+Self(n-5): n in [1..25]]; // Vincenzo Librandi, Aug 18 2013

LinearRecurrence[{1, 38, -38, -1, 1}, {0, 3, 15, 132, 588}, 23]
CoefficientList[Series[3 x (1 + 4 x + x^2)/((1 - x) (1 - 6 x - x^2) (1 + 6 x - x^2)), {x, 0, 25}], x] (* Vincenzo Librandi, Aug 18 2013 *)

makelist(expand(-1/2+((5+(-1)^n*sqrt(10))*(3-sqrt(10))^(2*floor(n/2))+(5-(-1)^n*sqrt(10))*(3+sqrt(10))^(2*floor(n/2)))/20), n, 1, 23);

x='x+O('x^30); concat([0], Vec(3*x*(1+4*x+x^2)/((1-x)*(1-6*x-x^2)*(1+6*x-x^2)))) \\ G. C. Greubel, Jul 15 2018

G.f.: 3*x*(1+4*x+x^2)/((1-x)*(1-6*x-x^2)*(1+6*x-x^2)).
a(n) = a(-n+1) = a(n-1)+38*a(n-2)-38*a(n-3)-a(n-4)+a(n-5).
a(n) = -1/2+((5+t*(-1)^n)*(3-t)^(2*floor(n/2))+(5-t*(-1)^n)*(3+t)^(2*floor(n/2)))/20, where t=sqrt(10).
2*a(n)+1 = A221875(n).

A261391 a(n) = n^5 + 5n^3 + 5n.

Original entry on oeis.org

0, 11, 82, 393, 1364, 3775, 8886, 18557, 35368, 62739, 105050, 167761, 257532, 382343, 551614, 776325, 1069136, 1444507, 1918818, 2510489, 3240100, 4130511, 5206982, 6497293, 8031864, 9843875, 11969386, 14447457, 17320268, 20633239, 24435150, 28778261, 33718432, 39315243, 45632114

Offset: 0

Raphael Ranna, Aug 17 2015

Also numbers of the form (n-th metallic mean)^5 - 1/(n-th metallic mean)^5, see link to Wikipedia.

Colin Barker, Table of n, a(n) for n = 0..1000
Wikipedia, Metallic mean.
Index entries for linear recurrences with constant coefficients, signature (6,-15,20,-15,6,-1).

Cf. A001622, A014176, A098316, A098317, A098318, A176398, A176439, A176458, A176522.

Array[#^5 + 5 #^3 + 5 # &, 34] (* Michael De Vlieger, Aug 18 2015 *)
Table[n^5 + 5*n^3 + 5*n, {n,0, 50}] (* G. C. Greubel, Aug 21 2015 *)
LinearRecurrence[{6,-15,20,-15,6,-1},{0,11,82,393,1364,3775},40] (* Harvey P. Dale, May 07 2018 *)

concat(0, Vec(x*(11*x^4+16*x^3+66*x^2+16*x+11)/(x-1)^6 + O(x^100))) \\ Colin Barker, Aug 18 2015

a(n) = ( (n+sqrt(n^2+4))/2 )^5 - 1/( (n+sqrt(n^2+4))/2 )^5.
a(n) = -a(-n) = 6*a(n-1)-15*a(n-2)+20*a(n-3)-15*a(n-4)+6*a(n-5)-a(n-6). - Colin Barker, Aug 18 2015
G.f.: x*(11*x^4+16*x^3+66*x^2+16*x+11) / (x-1)^6. - Colin Barker, Aug 18 2015
E.g.f.: (x^5 + 15*x^4 + 70*x^3 + 120*x^2 + 71*x + 11)*e^x. - G. C. Greubel, Aug 21 2015

Offset changed from 1 to 0, initial 0 added and b-file adapted from Bruno Berselli, Aug 25 2015

A218991 Power floor sequence of 3+sqrt(10).

Original entry on oeis.org

6, 36, 221, 1361, 8386, 51676, 318441, 1962321, 12092366, 74516516, 459191461, 2829665281, 17437183146, 107452764156, 662153768081, 4080375372641, 25144406003926, 154946811396196, 954825274381101, 5883898457682801

Offset: 0

Clark Kimberling, Nov 12 2012

See A214992 for a discussion of power floor sequence and the power floor function, p1(x) = lim_{n->oo} a(n,x)/x^n. The present sequence is a(n,r), where r = 3+sqrt(10), and the limit p1(r) = 5.815421188487681054332319082...

See A218992 for the power floor function, p4. For comparison with p1, we have lim_{r->oo} p4(r)/p1(r) = (3+sqrt(10))/5 = 1.23245553....

			a(0) = floor(r) = 6, where r = 3+sqrt(10);
a(1) = floor(6*r) = 36;
a(2) = floor(36*r) = 221.

Clark Kimberling, Table of n, a(n) for n = 0..250
Index entries for linear recurrences with constant coefficients, signature (7,-5,-1).

Cf. A214992, A005668, A015451, A218992.

Cf. A176398 (3+sqrt(10)).

[IsZero(n) select Floor(r) else Floor(r*Self(n)) where r is 3+Sqrt(10): n in [0..20]]; // Bruno Berselli, Nov 22 2012

x = 3 + Sqrt[10]; z = 30; (* z = # terms in sequences *)
f[x_] := Floor[x]; c[x_] := Ceiling[x];
p1[0] = f[x]; p2[0] = f[x]; p3[0] = c[x]; p4[0] = c[x];
p1[n_] := f[x*p1[n - 1]]
p2[n_] := If[Mod[n, 2] == 1, c[x*p2[n - 1]], f[x*p2[n - 1]]]
p3[n_] := If[Mod[n, 2] == 1, f[x*p3[n - 1]], c[x*p3[n - 1]]]
p4[n_] := c[x*p4[n - 1]]
t1 = Table[p1[n], {n, 0, z}]  (* A218991 *)
t2 = Table[p2[n], {n, 0, z}]  (* A005668 *)
t3 = Table[p3[n], {n, 0, z}]  (* A015451 *)
t4 = Table[p4[n], {n, 0, z}]  (* A218992 *)

a(n) = floor(r*a(n-1)), where r=3+sqrt(10), a(0) = floor(r).
a(n) = 7*a(n-1) - 5*a(n-2) - a(n-3).
G.f.: (6 - 6*x - x^2)/(1 - 7*x + 5*x^2 + x^3).
a(n) = ((5+sqrt(10))*(3-sqrt(10))^(n+2) + (5-sqrt(10))*(3+sqrt(10))^(n+2)+2)/12. - Bruno Berselli, Nov 22 2012

A218992 Power ceiling sequence of 3+sqrt(10).

Original entry on oeis.org

7, 44, 272, 1677, 10335, 63688, 392464, 2418473, 14903303, 91838292, 565933056, 3487436629, 21490552831, 132430753616, 816075074528, 5028881200785, 30989362279239, 190965054876220, 1176779691536560, 7251643204095581

Offset: 0

Clark Kimberling, Nov 12 2012

See A214992 for a discussion of power ceiling sequence and the power ceiling function, p4(x) = limit of a(n,x)/x^n. The present sequence is a(n,r), where r = 3+sqrt(10), and the limit p4(r) = 7.16724801485749657...

See A218991 for the power floor function, p1(x); for comparison of p1 and p4, we have limit(p4(r)/p1(r) = (3+sqrt(10))/5 = 1.23245553...

			a(0) = ceiling(r) = 7, where r = 3+sqrt(10);
a(1) = ceiling(7*r) = 44;
a(2) = ceiling(44*r) = 272.

Clark Kimberling, Table of n, a(n) for n = 0..250
Index entries for linear recurrences with constant coefficients, signature (7,-5,-1).

Cf. A214992, A005668, A015451, A218991.

Cf. A176398 (3+sqrt(10)).

[IsZero(n) select Ceiling(r) else Ceiling(r*Self(n)) where r is 3+Sqrt(10): n in [0..20]]; // Bruno Berselli, Nov 22 2012

(See A218991.)
LinearRecurrence[{7,-5,-1},{7,44,272},20] (* Harvey P. Dale, Sep 22 2016 *)

a(n) = ceiling(r*a(n-1)), where r=3+sqrt(10), a(0) = ceiling(r).
a(n) = 7*a(n-1) - 5*a(n-2) - a(n-3).
G.f.: (7 - 5*x - x^2)/(1 - 7*x + 5*x^2 + x^3).
a(n) = ((5+sqrt(10))*(3-sqrt(10))^(n+3)+(5-sqrt(10))*(3+sqrt(10))^(n+3)-10)/60. [Bruno Berselli, Nov 22 2012]

A261540 a(n) = n^7 + 7n^5 + 14n^3 + 7*n.

Original entry on oeis.org

0, 29, 478, 4287, 24476, 101785, 337434, 946043, 2333752, 5206581, 10714070, 20633239, 37597908, 65378417, 109216786, 176222355, 275832944, 420346573, 625528782, 911300591, 1302512140, 1829807049, 2530582538, 3450050347, 4642403496, 6172093925, 8115226054

Offset: 0

Raphael Ranna, Aug 24 2015

Also numbers of the form (n-th metallic mean)^7 - 1/(n-th metallic mean)^7, see link to Wikipedia.

Raphael Ranna, Table of n, a(n) for n = 0..100
Wikipedia, Metallic mean
Index entries for linear recurrences with constant coefficients, signature (8,-28,56,-70,56,-28,8,-1).

Cf. A001622, A014176, A098316, A098317, A098318, A176398, A176439, A176458, A176522, A261391.

[n^7 + 7*n^5 + 14*n^3 + 7*n: n in [0..30]]; // Vincenzo Librandi, Aug 24 2015

Table[n^7 + 7 n^5 + 14 n^3 + 7 n, {n, 0, 30}] (* Bruno Berselli, Aug 24 2015 *)
LinearRecurrence[{8, -28, 56, -70, 56, -28, 8, -1}, {0, 29, 478, 4287, 24476, 101785, 337434, 946043}, 30] (* Vincenzo Librandi, Aug 24 2015 *)

a(n)=n^7+7*n^5+14*n^3+7*n \\ Charles R Greathouse IV, Aug 24 2015

[n^7+7*n^5+14*n^3+7*n for n in (0..30)] # Bruno Berselli, Aug 24 2015

a(n) = -a(-n) = ( (n+sqrt(n^2+4))/2 )^7 - 1/( (n+sqrt(n^2+4))/2 )^7.

G.f.: x*(29 + 246*x + 1275*x^2 + 1940*x^3 + 1275*x^4 + 246*x^5 + 29*x^6)/(1 - x)^8. - Bruno Berselli, Aug 24 2015

Offset changed from 1 to 0 and initial 0 added by Bruno Berselli, Aug 25 2015

A261574 a(n) = n(n^2 + 3)(n^6 + 6n^4 + 9n^2 + 3).

Original entry on oeis.org

0, 76, 2786, 46764, 439204, 2744420, 12813606, 48229636, 153992264, 432083484, 1092730090, 2537720636, 5489037036, 11179326964, 21624372014, 40001698260, 71163830416, 122319408236, 203920464114, 330799604044, 523606640180, 810600392196, 1229857906486

Offset: 0

Raphael Ranna, Aug 24 2015

Also numbers of the form (n-th metallic mean)^9 - 1/(n-th metallic mean)^9, see link to Wikipedia.

Colin Barker, Table of n, a(n) for n = 0..1000 (first 100 terms from Raphael Ranna)
Wikipedia, Metallic mean
Index entries for linear recurrences with constant coefficients, signature (10,-45,120,-210,252,-210,120,-45,10,-1).

Cf. A001622, A014176, A098316, A098317, A098318, A176398, A176439, A176458, A176522, A261391, A261540.

[n*(n^2+3)*(n^6+6*n^4+9*n^2+3): n in [0..25]]; // Bruno Berselli, Aug 25 2015

Table[n (n^2 + 3) (n^6 + 6 n^4 + 9 n^2 + 3), {n, 0, 25}] (* Bruno Berselli, Aug 25 2015 *)

concat(0, Vec(2*x*(38*x^8 +1013*x^7 +11162*x^6 +43907*x^5 +69200*x^4 +43907*x^3 +11162*x^2 +1013*x +38) / (x -1)^10 + O(x^50))) \\ Colin Barker, Aug 25 2015

a(n) = -a(-n) = ( (n+sqrt(n^2+4))/2 )^9-1/( (n+sqrt(n^2+4))/2 )^9.

G.f.: 2*x*(38*x^8 +1013*x^7 +11162*x^6 +43907*x^5 +69200*x^4 +43907*x^3 +11162*x^2 +1013*x +38) / (x -1)^10. - Colin Barker, Aug 25 2015

Formula in Name by Colin Barker, Aug 25 2015

Offset changed from 1 to 0 and initial 0 added by Bruno Berselli, Aug 25 2015

cp's OEIS Frontend

A261327 a(n) = (n^2 + 4) / 4^((n + 1) mod 2).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

PARI

Python

Sage

Formula

Extensions

A221874 Numbers m such that 10*m^2 + 6 is a square.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Maxima

Formula

A179237 a(0) = 1, a(1) = 2; a(n+1) = 6*a(n) + a(n-1) for n>1.

Views

Author

Keywords

Comments

Examples

Links

Programs

Magma

Mathematica

PARI

Formula

Extensions

A244593 Decimal expansion of z_c = phi^5 (where phi is the golden ratio), a lattice statistics constant which is the exact value of the critical activity of the hard hexagon model.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A222390 Nonnegative integers m such that 10*m*(m+1)+1 is a square.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Magma

Mathematica

Maxima

PARI

Formula

A261391 a(n) = n^5 + 5*n^3 + 5*n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

A222390 Nonnegative integers m such that 10m(m+1)+1 is a square.

A261391 a(n) = n^5 + 5n^3 + 5n.

A261540 a(n) = n^7 + 7n^5 + 14n^3 + 7*n.

A261574 a(n) = n(n^2 + 3)(n^6 + 6n^4 + 9n^2 + 3).