login

Revision History for A005246

(Bold, blue-underlined text is an addition; faded, red-underlined text is a deletion.)

Showing entries 1-10 | older changes
a(n) = (1 + a(n-1)*a(n-2))/a(n-3), a(0) = a(1) = a(2) = 1.
(history; published version)
#118 by Alois P. Heinz at Thu Jul 27 12:37:47 EDT 2023
STATUS

editing

approved

#117 by Alois P. Heinz at Thu Jul 27 12:37:17 EDT 2023
COMMENTS

From Eric Simon Jacob, Jul 20 2023: (Start)

G.f.: (1 + x - 3*x^2 - 2*x^3) / (x^4 - 4*x^2 + 1).

a(n) = 0*a(n-1) + 4*a(n-2) + 0*a(n-3) - 1*a(n-4) for n > 3, a(0) = 1, a(1) = 1, a(2) = 1, a(3) = 2.

A = (989*sqrt(3)+1713-(2340*sqrt(3)+4053)*sqrt(2-sqrt(3)))/(9360*sqrt(3)+16212)

B = ((45*sqrt(3)+78)*sqrt(2-sqrt(3))+19*sqrt(3)+33)/(180*sqrt(3)+312)

C = ((2*sqrt(3)+3)*sqrt(2+sqrt(3))+3*sqrt(3)+5)/(8*sqrt(3)+12)

D = (3*sqrt(3)+5-(2*sqrt(3)+3)*sqrt(2+sqrt(3)))/(8*sqrt(3)+12)

A ~ -0.0237470898486668

B ~ 0.235071955253854

C ~ 0.877300480441941

D ~ -0.0886253458471277

a(n) = ((-1)^n*A + B)*sqrt(2+sqrt(3))^n + (C + (-1)^n*D)*sqrt(2-sqrt(3))^n.

(End)

#116 by Alois P. Heinz at Thu Jul 27 09:27:09 EDT 2023
STATUS

proposed

editing

Discussion
Thu Jul 27
09:32
Alois P. Heinz: "G.f.: (1 + x - 3*x^2 - 2*x^3)/(1 - 4*x^2 + x^4)." is the first formula in the formula section ... there is no need to duplicate this in the form "(1 + x - 3*x^2 - 2*x^3) / (x^4 - 4*x^2 + 1)" ... first read then write ...
12:36
Alois P. Heinz: no need for duplicates ... no need for overly complicated formulas ...
#115 by Jon E. Schoenfield at Thu Jul 27 05:49:02 EDT 2023
STATUS

editing

proposed

Discussion
Thu Jul 27
09:27
Alois P. Heinz: 0*a(n-1) = 0 ...
#114 by Jon E. Schoenfield at Thu Jul 27 05:47:48 EDT 2023
COMMENTS

For n > 2, a(n) is the smallest integer > a(n-1) such that sqrt(3)*a(n) is closer to and greater than an integer than sqrt(3)*a(n-1). I.e., a(n) is the smallest integer > a(n-1) such that frac(sqrt(3)*a(n)) < frac(sqrt(3)*a(n-1)). - Benoit Cloitre, Jan 20 2003

This sequence is a particular case of the following situation: a(0)=1, a(1)=a, a(2)=b with the recurrence relation a(n+3)=(a(n+2)*a(n+1)+q)/a(n) where q is given in Z to have Q=(a*b^2+q*b+a+q)/(a*b) itself in Z. The g.f . is f: f(z)=(1+a*z+(b-Q)*z^2+(a*b+q-a*Q)*z^3)/(1-Q*z^2+z^4); so we have the linear recurrence: a(n+4)=Q*a(n+2)-a(n). The general form of a(n) is given by: a(2*m) = Sum_{p=sum0..floor(m/2)} (-1)^p*binomial(m-p,p)*Q^(m-2*p), + (b-Q)*Sum_{p=0..floor((m-1)/2))+(b-Q)*sum(} (-1)^p*binomial(m-1-p,p)*Q^(m-1-2*p), and a(2*m+1) = a*Sum_{p=0..floor((m-1)/2)) and a(2*m+1)=a*sum(} (-1)^p*binomial(m-p,p)*Q^(m-2*p), + (a*b+q-a*Q)*Sum_{p=0..floor((m-1)/2))+(a*b+q-a*Q)*sum(} (-1)^p*binomial(m-1-p,p)*Q^(m-1-2*p),p=0..floor((m-1)/2)). - Richard Choulet, Feb 24 2010

LINKS

T. D. Noe, <a href="/A005246/b005246.txt">Table of n, a(n) for n = 0..500</a>

FORMULA

Limit_{n->infinityoo} a(2n+1)/a(2n) = (3+sqrt(3))/3 = 1.5773502...; lim_{n->infinityoo} a(2n)/a(2n-1) = (3+sqrt(3))/2 = 2.3660254.... - Benoit Cloitre, Aug 07 2002

a(n) = a(2 - n) for all n in Z. - Michael Somos, Nov 15 2006

For n > 2: a(n) = a(n-1) + Sum_{k=1..floor((n-1)/2)} a(2*k). - Reinhard Zumkeller, Dec 16 2007

Discussion
Thu Jul 27
05:49
Jon E. Schoenfield: Okay, I’m done editing for now. :-)
#113 by Jon E. Schoenfield at Thu Jul 27 05:37:25 EDT 2023
COMMENTS

In the closed form formula (sqrt(2+sqrt(3))^n)=((sqrt(6)+sqrt(2))/2)^n;(-sqrt(2+sqrt(3))^n)=((-sqrt(6)-sqrt(2))/2)^n;(sqrt(2-sqrt(3))^n)=((sqrt(6)-sqrt(2))/2)^n ;(-sqrt(2-sqrt(3))^n)=((sqrt(2)-sqrt(6))/2)^n. - Tim Monahan, Jul 07 2011

From Tim Monahan, Jul 07 2011: (Start)

In the closed-form formula,

sqrt(2+sqrt(3))^n = ((sqrt(6)+sqrt(2))/2)^n;

-sqrt(2+sqrt(3))^n = ((-sqrt(6)-sqrt(2))/2)^n;

sqrt(2-sqrt(3))^n = ((sqrt(6)-sqrt(2))/2)^n;

-sqrt(2-sqrt(3))^n = ((sqrt(2)-sqrt(6))/2)^n.

(End)

STATUS

proposed

editing

#112 by Jon E. Schoenfield at Thu Jul 27 05:31:24 EDT 2023
STATUS

editing

proposed

#111 by Jon E. Schoenfield at Thu Jul 27 05:24:04 EDT 2023
COMMENTS

Integer solutions to the equation floor(sqrt(3)*x^2) = x*floor(sqrt(3)*x). - Benoit Cloitre, Mar 18 2004

For n > 2, a(n) is the smallest integer > a(n-1) such that sqrt(3)*a(n) is closer to and greater than an integer than sqrt(3)*a(n-1). iI.e. , a(n) is the smallest integer > a(n-1) such that frac(sqrt(3)*a(n))<frac(sqrt(3)*a(n-1)). - Benoit Cloitre, Jan 20 2003

a(n) for n > 1 are the integer square roots of (floor(m^2/3)+1), where the values of m are given by A143643. Also see A082630. - Richard R. Forberg, Nov 14 2013

a(n) = 0*a(n-1) + 4*a(n-2) + 0*a(n-3) - 1*a(n-4) for n > 3, a(0) = 1, a(1) = 1, a(2) = 1, a(3) = 2.

STATUS

proposed

editing

Discussion
Thu Jul 27
05:26
Jon E. Schoenfield: Thanks … but the same g.f. Is already given in the Formula section … and I don’t understand the rationale for including terms with a coefficient of 0.
05:31
Jon E. Schoenfield: It seems to me that the remaining content in your contribution (reorganized to put the definitions of A, B, C, and D at the end) belongs in the Formula section, not the Comments section (provided that it doesn’t duplicate existing content).
#110 by Eric Simon Jacob at Thu Jul 27 05:18:28 EDT 2023
STATUS

editing

proposed

#109 by Eric Simon Jacob at Thu Jul 20 13:16:11 EDT 2023
COMMENTS

a(n) = 0.*a(n-1) + 4.*a(n-2) + 0.*a(n-3) -1.*a(n-4) for n>3, a(0) = 1, a(1) = 1, a(2) = 1, a(3) = 2.

Discussion
Thu Jul 27
05:18
Eric Simon Jacob: i did the correction Jon E. Schoenfield.