A216162 -id:A216162 - cp's OEIS frontend

A006454 Solution to a Diophantine equation: each term is a triangular number and each term + 1 is a square.

0, 3, 15, 120, 528, 4095, 17955, 139128, 609960, 4726275, 20720703, 160554240, 703893960, 5454117903, 23911673955, 185279454480, 812293020528, 6294047334435, 27594051024015, 213812329916328, 937385441796000, 7263325169820735, 31843510970040003, 246739243443988680

Offset: 0

N. J. A. Sloane, Jeffrey Shallit

Alternative definition: a(n) is triangular and a(n)/2 is the harmonic average of consecutive triangular numbers. See comments and formula section of A005563, of which this sequence is a subsequence. - Raphie Frank, Sep 28 2012
As with the Sophie Germain triangular numbers (A124174), 35 = (a(n) - a(n-6))/(a(n-2) - a(n-4)). - Raphie Frank, Sep 28 2012
Sophie Germain triangular numbers of the second kind as defined in A217278. - Raphie Frank, Feb 02 2013
Triangular numbers m such that m+1 is a square. - Bruno Berselli, Jul 15 2014
From Vladimir Pletser, Apr 30 2017: (Start)
Numbers a(n) which are the triangular number T(b(n)), where b(n) is the sequence A006451(n) of numbers n such that T(n)+1 is a square.
a(n) also gives the x solutions of the 3rd-degree Diophantine Bachet-Mordell equation y^2 = x^3 + K, with y = T(b(n))*sqrt(T(b(n))+1) = A285955(n) and K = T(b(n))^2 = A285985(n), the square of the triangular number of b(n) = A006451(n).
Also: This sequence is a subsequence of A000217(n), namely A000217(A006451(n)). (End)

			From _Raphie Frank_, Sep 28 2012: (Start)
35*(528 - 15) + 0 = 17955 = a(6),
35*(4095 - 120) + 3 = 139128 = a(7),
35*(17955 - 528) + 15 = 609960 = a(8),
35*(139128 - 4095) + 120 = 4726275 = a(9). (End)
From _Raphie Frank_, Feb 02 2013: (Start)
a(7) = 139128 and a(9) = 4726275.
a(9) = (2*(sqrt(8*a(7) + 1) - 1)/2 + 3*sqrt(a(7) + 1) + 1)^2 - 1 = (2*(sqrt(8*139128 + 1) - 1)/2 + 3*sqrt(139128 + 1) + 1)^2 - 1 = 4726275.
a(9) = 1/2*((3*(sqrt(8*a(7) + 1) - 1)/2 + 4*sqrt(a(7) + 1) + 1)^2 + (3*(sqrt(8*a(7) + 1) - 1)/2 + 4*sqrt(a(7) + 1) + 1)) = 1/2*((3*(sqrt(8*139128 + 1) - 1)/2 + 4*sqrt(139128 + 1) + 1)^2 + (3*(sqrt(8*139128 + 1) - 1)/2 + 4*sqrt(139128 + 1) + 1)) = 4726275. (End)
From _Vladimir Pletser_, Apr 30 2017: (Start)
For n=2, b(n)=5, a(n)=15
For n=5, b(n)=90, a(n)= 4095
For n = 3, A006451(n) = 15. Therefore, A000217(A006451(n)) = A000217(15) = 120. (End)

Edward J. Barbeau, Pell's Equation, New York: Springer-Verlag, 2003, p. 17, Exercise 1.2.
Allan J. Gottlieb, How four dogs meet in a field, and other problems, Technology Review, Jul/August 1973, pp. 73-74.
Vladimir Pletser, On some solutions of the Bachet-Mordell equation for large parameter values, to be submitted, April 2017.
Jeffrey Shallit, personal communication.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Vladimir Pletser, Table of n, a(n) for n = 0..1000 (first 60 terms from Vincenzo Librandi)
M.A. Bennett and A. Ghadermarzi, Data on Mordell's curve.
Michael A. Bennett and Amir Ghadermarzi, Mordell's equation : a classical approach, arXiv:1311.7077 [math.NT], 2013.
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.
Jeffrey Shallit, Letter to N. J. A. Sloane, Oct. 1975.
K. B. Subramaniam, Almost Square Triangular Numbers, The Fibonacci Quarterly, Vol. 37, No. 3 (1999), pp. 194-197.
Eric Weisstein's World of Mathematics, Mordell Curve.
Index entries for linear recurrences with constant coefficients, signature (1,34,-34,-1,1).

Cf. sqrt(a(n) + 1) = A006452(n + 1) = A216162(2n + 2) and (sqrt(8a(n) + 1) - 1)/2 = A006451.
Cf. A217278, A124174, A216134. - Raphie Frank, Feb 02 2013
Subsequence of A182334.
Cf. A001109, A245031.
Cf. A285955, A285985, A000217, A006451, A081119, A054504. - Vladimir Pletser, Apr 30 2017

I:=[0,3,15,120,528,4095]; [n le 6 select I[n] else 35*(Self(n-2) - Self(n-4)) + Self(n-6): n in [1..30]]; // Vincenzo Librandi, Dec 21 2015

A006454:=-3*z*(1+4*z+z**2)/(z-1)/(z**2-6*z+1)/(z**2+6*z+1); # conjectured (correctly) by Simon Plouffe in his 1992 dissertation
restart: bm2:=-1: bm1:=0: bp1:=2: bp2:=5: print ('0,0','1,3','2,15'); for n from 3 to 1000 do b:= 8*sqrt((bp1^2+bp1)/2+1)+bm2; a:=b*(b+1)/2; print(n,a); bm2:=bm1; bm1:=bp1; bp1:=bp2; bp2:=b; end do: # Vladimir Pletser, Apr 30 2017

Clear[a]; a[0] = a[1] = 1; a[2] = 2; a[3] = 4; a[n_] := 6a[n - 2] - a[n - 4]; Array[a, 40]^2 - 1 (* Vladimir Joseph Stephan Orlovsky, Mar 03 2011 *)
LinearRecurrence[{1,34,-34,-1,1},{0,3,15,120,528},30] (* Harvey P. Dale, Feb 18 2023 *)

concat(0, Vec(3*x*(1 + 4*x + x^2) / ((1 - x)*(1 - 6*x + x^2)*(1 + 6*x + x^2)) + O(x^30))) \\ Colin Barker, Apr 30 2017

a(n) = A006451(n)*(A006451(n)+1)/2.
a(n) = A006452(n)^2 - 1. - Joerg Arndt, Mar 04 2011
a(n) = 35*(a(n-2) - a(n-4)) + a(n-6). - Raphie Frank, Sep 28 2012
From Raphie Frank, Feb 01 2013: (Start)
a(0) = 0, a(1) = 3, and a(n+2) = (2x + 3y + 1)^2 - 1  = 1/2*((3x + 4y + 1)^2 + (3x + 4y + 1)) where x = (sqrt(8*a(n) + 1) - 1)/2 = A006451(n) =  1/2*(A216134(n + 1) + A216134(n - 1)) and y = sqrt(a(n) + 1) = A006452(n + 1) = 1/2*(A216134(n + 1) - A216134(n - 1)).
Note that A216134(n + 1) = x + y, and A216134(n + 3) = (2x + 3y + 1) + (3x + 4y + 1) = (5x + 7y + 2), where A216134 gives the indices of the Sophie Germain triangular numbers. (End)
a(n) = (1/64)*(((4 + sqrt(2))*(1 -(-1)^(n+1)*sqrt(2))^(2* floor((n+1)/2)) + (4 - sqrt(2))*(1 + (-1)^(n+1)*sqrt(2))^(2*floor((n+1)/2))))^2 - 1. - Raphie Frank, Dec 20 2015
From Vladimir Pletser, Apr 30 2017: (Start)
Since b(n) = 8*sqrt(T(b(n-2))+1)+ b(n-4) = 8*sqrt((b(n-2)*(b(n-2)+1)/2)+1)+ b(n-4), with b(-1)=-1, b(0)=0, b(1)=2, b(2)=5 (see A006451) and a(n) = T(b(n)) (this sequence), we have:
a(n) = ((8*sqrt((b(n-2)*(b(n-2)+1)/2)+1)+ b(n-4))*(8*sqrt((b(n-2)*(b(n-2)+1)/2)+1)+ b(n-4)+1)/2). (End)
From Colin Barker, Apr 30 2017: (Start)
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) + 34*a(n-2) - 34*a(n-3) - a(n-4) + a(n-5) for n > 4.
(End)
a(n)  = (A001109(n/2+1) - 2*A001109(n/2))^2 - 1 if n is even, and (A001109((n+3)/2) - 4*A001109((n+1)/2))^2 - 1 if n is odd (Subramaniam, 1999). - Amiram Eldar, Jan 13 2022

Better description from Harvey P. Dale, Jan 28 2001
More terms from Larry Reeves (larryr(AT)acm.org), Feb 07 2001
Minor edits by N. J. A. Sloane, Oct 24 2009

A266504 a(n) = 2*a(n - 2) + a(n - 4) with a(0) = a(1) = 2, a(2) = 1, a(3) = 3.

Original entry on oeis.org

2, 2, 1, 3, 4, 8, 9, 19, 22, 46, 53, 111, 128, 268, 309, 647, 746, 1562, 1801, 3771, 4348, 9104, 10497, 21979, 25342, 53062, 61181, 128103, 147704, 309268, 356589, 746639, 860882, 1802546, 2078353, 4351731, 5017588, 10506008, 12113529, 25363747, 29244646, 61233502

Offset: 0

Raphie Frank, Dec 30 2015

This sequence gives all x in N | 2*x^2 - 7(-1)^x = y^2. The companion sequence to this sequence, giving y values, is A266505.
A266505(n)/a(n) converges to sqrt(2).
Alternatively, 1/4*(3*A002203(floor[n/2]) - A002203(n-(-1)^n)), where A002203 gives the Companion Pell numbers, or, in Lucas sequence notation, V_n(2, -1).
Alternatively, bisection of A266506.
Alternatively, A048654(n -1) and A078343(n + 1) interlaced.
Alternatively, A100525(n-1), A266507(n), A038761(n) and A253811(n) interlaced.
Let b(n) = (a(n) - a(n)(mod 2))/2, that is b(n) = {1, 1, 0, 1, 2, 4, 4, 9, 11, 23, 26, 55, 64, ...}. Then:
A006452(n) = {b(4n+0) U b(4n+1)} gives n in N such that n^2 - 1 is triangular;
A216134(n) = {b(4n+2) U b(4n+3)} gives n in N such that n^2 + n + 1 is triangular (indices of Sophie Germain triangular numbers);
A216162(n) = {b(4n+0) U b(4n+2) U b(4n+1) U b(4n+3)}, sequences A006452 and A216134 interlaced.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (0,2,0,1).

Cf. A000129, A001333, A002203, A002965, A006451, A006452, A002965, A038761, A038762, A048654, A048655, A054490, A078343, A098586, A098790, A100525, A101386, A135532, A216134, A216162, A253811, A255236, A266504, A266505, A266507.

I:=[2,2,1,3]; [n le 4 select I[n] else 2*Self(n-2)+Self(n-4): n in [1..70]]; // Vincenzo Librandi, Dec 31 2015

LinearRecurrence[{0, 2, 0, 1}, {2, 2, 1, 3}, 70] (* Vincenzo Librandi, Dec 31 2015 *)
Table[SeriesCoefficient[(1 - x) (2 + 4 x + x^2)/(1 - 2 x^2 - x^4), {x, 0, n}], {n, 0, 41}] (* Michael De Vlieger, Dec 31 2015 *)

Vec((1-x)*(2+4*x+x^2)/(1-2*x^2-x^4) + O(x^50)) \\ Colin Barker, Dec 31 2015

a(n) = 1/sqrt(8)*(+sqrt(2)*(1+sqrt(2))^(floor(n/2)-(-1)^n)*(-1)^n - 3*(1-sqrt(2))^(floor(n/2)-(-1)^n) + sqrt(2)*(1-sqrt(2))^(floor(n/2)-(-1)^n)*(-1)^n + 3*(1+sqrt(2))^(floor(n/2)-(-1)^n)).
a(n) = 1/4*((3*((1+sqrt(2))^floor(n/2)+(1-sqrt(2))^floor(n/2))) - (-1)^n*((1+sqrt(2))^(floor(n/2)-(-1)^n)+(1-sqrt(2))^(floor(n/2)-(-1)^n))).
a(2n) = (+sqrt(2)*(1+sqrt(2))^(n-1) - 3 *(1-sqrt(2))^(n-1) + sqrt(2)*(1-sqrt(2))^(n-1) + 3*(1 + sqrt(2))^(n-1))/sqrt(8) = A048654(n -1).
a(2n) = 1/4*((3*((1+sqrt(2))^n+(1-sqrt(2))^n)) - ((1+sqrt(2))^(n-1)+(1-sqrt(2))^(n-1))) = A048654(n -1).
a(2n + 1) = (-sqrt(2)*(1+sqrt(2))^(n+1) - 3 *(1-sqrt(2))^(n+1) - sqrt(2)*(1-sqrt(2))^(n+1) + 3*(1+sqrt(2))^(n+1))/sqrt(8) = A078343(n + 1).
a(2n + 1) =1/4*((3*((1+sqrt(2))^n+(1-sqrt(2))^n)) + ((1+sqrt(2))^(n+1)+(1-sqrt(2))^(n+1))) =  A078343(n + 1).
a(4n + 0) = 6*a(4n - 4) - a(4n - 8) = A100525(n-1).
a(4n + 1) = 6*a(4n - 3) - a(4n - 7) = A266507(n).
a(4n + 2) = 6*a(4n - 2) - a(4n - 6) = A038761(n).
a(4n + 3) = 6*a(4n - 1) - a(4n - 5) = A253811(n).
sqrt(2*a(n)^2 - 7(-1)^a(n))*sgn(2*n - 1) = A266505(n).
(a(2n + 1) + a(2n))/2 = A002203(n), where A002203 gives the companion Pell numbers.
(a(2n + 1) - a(2n))/2 = A000129(n), where A000129 gives the Pell numbers.
(a(2n+2) + a(2n+1))*2 = A002203(n+2)
(a(2n+2) - a(2n+1))*2 = A002203(n-1).
G.f.: (1-x)*(2+4*x+x^2) / (1-2*x^2-x^4). - Colin Barker, Dec 31 2015

A073325 a(n) = least k > 0 such that prime(k) == n (mod k).

Original entry on oeis.org

1, 2, 3, 4, 75, 9, 79, 18, 17, 10, 19, 20, 91, 22, 23, 41, 83, 24, 16049, 43, 2711, 94, 25, 26, 95, 198, 449, 452, 99, 50, 451, 48, 453, 1072, 447, 54, 16043, 55, 2719, 56, 459, 57, 101, 472, 100371, 62, 105, 102, 103, 104, 467, 110, 107, 65, 109, 63, 115, 118, 117

Offset: 1

Labos Elemer, Jul 30 2002

nonn

First appearance of n-1 in A004648. Are all positive integers present in A004648 and hence in this sequence? - Zak Seidov, Sep 02 2012

			a(4) = 75 as prime(75) = 379 == 4 (mod 75).
a(44) = 100371 since prime(100371) = 1304867 == 44 (mod 100371) and prime(k) <> 44 (mod k) for k < 100371.

Zak Seidov, Table of n, a(n) for n = 1..301

Cf. A000040, A002808, A004648, A073324, A073326.

nn = 60; f[x_] := Mod[Prime[x], x]; t = Table[0, {nn}]; k = 0; While[Times @@ t == 0, k++; n = f[k]; If[n <= nn && t[[n]] == 0, t[[n]] = k]]; Join[{1}, t]
lk[n_]:=Module[{k=1},While[Mod[Prime[k],k]!=n,k++];k]; Array[lk,60,0] (* Harvey P. Dale, Nov 29 2013 *)

stop=110000; for(n=0,59,k=1; while(k

from sympy import prime, nextprime
def A073325(n):
    p, m = prime(n), n
    while p%m != n-1:
        p = nextprime(p)
        m += 1
    return m # Chai Wah Wu, Mar 18 2023

a(n) = Min{x; Mod[A000040(x), x]=n} = Min{x; A004648[x]=n}.

Definition revised by N. J. A. Sloane, Aug 12 2009

A216162 merged into this sequence by T. D. Noe, Sep 07 2012

A217278 Sequences A124174 and A006454 interlaced.

Original entry on oeis.org

0, 0, 1, 3, 10, 15, 45, 120, 351, 528, 1540, 4095, 11935, 17955, 52326, 139128, 405450, 609960, 1777555, 4726275, 13773376, 20720703, 60384555, 160554240, 467889345, 703893960, 2051297326, 5454117903, 15894464365, 23911673955, 69683724540, 185279454480

Offset: 0

Raphie Frank, Sep 29 2012

a(2n) and 2*a(2n) + 1 are triangular.

a(2n + 1) is triangular and a(2n + 1)/2 is the harmonic mean of consecutive triangular numbers (therefore, a(2n + 1) + 1 is square).

			a(18) = 35*(52326 - 1540) + 45 = 1777555,
a(19) = 35*(139128 - 4095) + 120 = 4726275.

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

Cf. (sqrt(8a(2n) + 1) - 1)/2 = A216134(n) = A216162(2n + 1).
Cf. sqrt(a(2n+1) + 1) = A006452(n + 1) = A216162(2n + 2).
Cf. (sqrt(8a(2n+1) + 1) - 1)/2 = A006451(n).

LinearRecurrence[{0,1,0,34,0,-34,0,-1,0,1},{0,0,1,3,10,15,45,120,351,528},40] (* Harvey P. Dale, Aug 04 2019 *)

concat([0,0], Vec(-x^2*(3*x^5+x^4+12*x^3+9*x^2+3*x+1)/((x-1)*(x+1)*(x^2-2*x-1)*(x^2+2*x-1)*(x^4+6*x^2+1)) + O(x^100))) \\ Colin Barker, Jun 23 2015

a(n) = 35*(a(n-4) - a(n-8)) + a(n-12).
lim n --> infinity a(2n)/a(2n - 1) = (3 + sqrt(8))/2.
From Raphie Frank, Dec 21 2015: (Start)
a(2*n + 1) = 1/64*(((4 + sqrt(2)) * (1 - (-1)^(n+1) * sqrt(2))^(2*floor((n+1)/2)) + (4 - sqrt(2)) * (1+(-1)^(n+1) * sqrt(2))^(2*floor((n+1)/2))))^2 - 1.
a(2*n + 2) = 1/2*(3*(a(2*n + 1)) + sqrt((a(2*n + 1)) + 1) * sqrt(8*(a(2*n + 1)) + 1) + 1).
(End)

A266505 a(n) = 2*a(n - 2) + a(n - 4) with a(0) = -1, a(1) = 1, a(2) = 3, a(3) = 5.

Original entry on oeis.org

-1, 1, 3, 5, 5, 11, 13, 27, 31, 65, 75, 157, 181, 379, 437, 915, 1055, 2209, 2547, 5333, 6149, 12875, 14845, 31083, 35839, 75041, 86523, 181165, 208885, 437371, 504293, 1055907, 1217471, 2549185, 2939235, 6154277, 7095941, 14857739, 17131117, 35869755, 41358175, 86597249, 99847467

Offset: 0

Raphie Frank, Dec 30 2015

a(n)/A266504(n) converges to sqrt(2).
Alternatively, bisection of A266506.
Alternatively, A135532(n) and A048655(n) interlaced.
Alternatively, A255236(n-1), A054490(n), A038762(n) and A101386(n) interlaced.
Let b(n) = (a(n) - (a(n) mod 2))/2, that is b(n) = {-1, 0, 1, 2, 2, 5, 6, 13, 15, 32, 37, 78, 90, ...}. Then:
A006451(n) = {b(4n+0) U b(4n+1)} gives n in N such that triangular(n) + 1 is square;
A216134(n) = {b(4n+2) U b(4n+3)} gives n in N such that triangular(n) follows form n^2 + n + 1 (twice a triangular number + 1).

G. C. Greubel, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (0,2,0,1).

I:=[-1,1,3,5]; [n le 4 select I[n] else 2*Self(n-2)+Self(n-4): n in [1..70]]; // Vincenzo Librandi, Dec 31 2015

a:=proc(n) option remember; if n=0 then -1 elif n=1 then 1 elif n=2 then 3 elif n=3 then 5 else 2*a(n-2)+a(n-4); fi; end:  seq(a(n), n=0..50); # Wesley Ivan Hurt, Jan 01 2016

LinearRecurrence[{0, 2, 0, 1}, {-1, 1, 3, 5}, 70] (* Vincenzo Librandi, Dec 31 2015 *)
Table[SeriesCoefficient[(-1 + 3 x) (1 + x)^2/(1 - 2 x^2 - x^4), {x, 0, n}], {n, 0, 42}] (* Michael De Vlieger, Dec 31 2015 *)

my(x='x+O('x^40)); Vec((-1+3*x)*(1+x)^2/(1-2*x^2-x^4)) \\ G. C. Greubel, Jul 26 2018

G.f.: (-1 + 3*x)*(1 + x)^2/(1 - 2*x^2 - x^4).
a(n) = (-(1+sqrt(2))^floor(n/2)*(-1)^n - sqrt(8)*(1-sqrt(2))^floor(n/2) - (1-sqrt(2))^floor(n/2)*(-1)^n + sqrt(8)*(1+sqrt(2))^floor(n/2))/2.
a(n) = 3*(((1+sqrt(2))^floor(n/2)-(1-sqrt(2))^floor(n/2))/sqrt(8)) - (-1)^n*(((1+sqrt(2))^(floor(n/2)-(-1)^n)-(1-sqrt(2))^(floor(n/2)-(-1)^n))/sqrt(8)).
a(n) = (3*A000129(floor(n/2)) - A000129(n-(-1)^n)), where A000129 gives the Pell numbers.
a(n) = sqrt(2*A266504(n)^2 - 7*(-1)^A266504(n))*sgn(2*n-1), where A266504 gives all x in N such that 2*x^2 - 7*(-1)^x = y^2. This sequence gives associated y values.
a(2n) = (-(1 + sqrt(2))^n - sqrt(8)*(1 - sqrt(2))^n - (1 - sqrt(2))^n + sqrt(8)*(1 + sqrt(2))^n)/2 = a(2n) = A135532(n).
a(2n) = 3*(((1+sqrt(2))^n-(1-sqrt(2))^n)/sqrt(8)) - (((1+sqrt(2))^(n-1)-(1-sqrt(2))^(n-1))/sqrt(8)) = A135532(n).
a(2n+1) = (+(1 + sqrt(2))^n - sqrt(8)*(1 - sqrt(2))^n + (1 - sqrt(2))^n + sqrt(8)*(1 + sqrt(2))^n)/2 = a(2n + 1) = A048655(n).
a(2n+1) = 3*(((1+sqrt(2))^n-(1-sqrt(2))^n)/sqrt(8)) + (((1+sqrt(2))^(n+1)-(1-sqrt(2))^(n+1))/sqrt(8)) = A048655(n).
a(4n + 0) = 6*a(4n - 4) - a(4n - 8) = A255236(n-1).
a(4n + 1) = 6*a(4n - 3) - a(4n - 7) = A054490(n).
a(4n + 2) = 6*a(4n - 2) - a(4n - 6) = A038762(n).
a(4n + 3) = 6*a(4n - 1) - a(4n - 5) = A101386(n).
(sqrt(2*(a(2n + 1) )^2 + 14*(-1)^floor(n/2)))/2 = A266504(n).
(a(2n + 1) + a(2n))/8 = A000129(n), where A000129 gives the Pell numbers.
a(2n + 1) - a(2n) = A002203(n), where A002203 gives the companion Pell numbers.
(a(2n + 2) + a(2n + 1))/2 = A000129(n+2).
(a(2n + 2) - a(2n + 1))/2 = A000129(n-1).

A211202 Positive numbers n such that Lambda_n = A002336(n) is divisible by n.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 12, 15, 16, 18, 20, 21, 22, 23, 24

Offset: 1

Raphie Frank, Feb 18 2013

nonn

Observations:
For all n in this sequence to n = 24, then y = Lambda_n/n follows form: y = (x^2 + x^k) - (floor[z^2/4]) or y = (x^2 + x^k) + (floor[z^2/4]); k = 1 or 2 and z = 0, 1, 3, 6 or 7. y (= A222786) gives the average number of spheres/dimension of the laminated lattice Kissing numbers in A222785.
e.g. Where T_x is the x-th triangular number = (1/2*(x^2 + x)), 2*T_x is the x-th pronic number = (x^2 + x) = floor[(2*x + 1)^2/4], and S_x is the x-th square = (x^2) = floor[(2*x)^2/4]:
For k = 1, z = 0 or 1, then n = {1, 4, 6, 8, 15, 20, 24}, x = {1, 2, 3, 5, 12, 29, 90}, and y = 2*T_x = {2, 6, 12, 30, 156, 870, 8190}.
For k = 2, z = 0 or 1, then n = {1, 5, 7, 23}, x = {1, 2, 3, 45}, and y = 2*T_x + 2*T_(-x) = 2*S_x = {2, 8, 18, 4050}.
For k = 1, z = 3, then n = {3, 7, 12, 16}, x = {2, 4, 7, 16}, and y = 2*T_x - 2*T_1 = {4, 18, 54, 270}.
For k = 1, z = 6, then n = {2, 18}, x = {3, 20}, and y = 2*T_x - S_3 = {3, 411}.
For k = 1, z = 7, then n = {5, 7, 8, 21}, x = {4, 5, 6, 36}, and y = 2*T_x - 2*T_3 = {8, 18, 30, 1320}.
For k = 1, z = 7, then n = {6, 7, 12, 22}, x = {0, 2, 6, 47}, and y = 2*T_x + 2*T_3 = {12, 18, 54, 2268}.
For the special case where k = 1 and z = 0 or 1, then all associated x values follow form (A216162(n) - A216162(n - 2)) [type 1] or (A216162(n) - A216162(n - 1)) [type II] for some n in N. Type II x values = {1, 2, 5, 90} (= A215797(n+1)) are associated with the positive Ramanujan-Nagell triangular numbers = {1, 3, 15, 4095} (= A076046(n+1)) by the formula 1/2*(x^2 + x) = T_x.

			Lambda_6/6 = 72/6 = 12, so 6 is in this sequence.
Lambda_12/12 = 648/12 = 54, so 12 is in this sequence.
Lambda_18/18 = 7398/18 = 411, so 18 is in this sequence.
Lambda_24/24 = 196560/24 = 8190, so 24 is in this sequence.
But...
Lambda_19/19 = 10668/19 = 561.47368..., so 19 is not in this sequence.

Cf. A222785, A222786, A002336, A216162

cp's OEIS Frontend

A006454 Solution to a Diophantine equation: each term is a triangular number and each term + 1 is a square.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

Extensions

A266504 a(n) = 2*a(n - 2) + a(n - 4) with a(0) = a(1) = 2, a(2) = 1, a(3) = 3.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A073325 a(n) = least k > 0 such that prime(k) == n (mod k).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

Extensions

A217278 Sequences A124174 and A006454 interlaced.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A266505 a(n) = 2*a(n - 2) + a(n - 4) with a(0) = -1, a(1) = 1, a(2) = 3, a(3) = 5.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

A211202 Positive numbers n such that Lambda_n = A002336(n) is divisible by n.

Views

Author

Keywords

Comments

Examples

Crossrefs