OFFSET
0,3
COMMENTS
a(n) > 0 is a square such that a(n) - 1 is a product of powers. - Michel Lagneau, Feb 16 2012
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 0..654
M. F. Hasler, Truncated squares, OEIS wiki, Jan 16 2012
Giovanni Lucca, Integer Sequences and Circle Chains Inside a Circular Segment, Forum Geometricorum, Vol. 18 (2018), 47-55.
Giovanni Lucca, Circle chains inside the arbelos and integer sequences, Int'l J. Geom. (2023) Vol. 12, No. 1, 71-82.
Index entries for linear recurrences with constant coefficients, signature (35,-35,1).
FORMULA
From Colin Barker, Sep 15 2014: (Start)
a(n) = 35*a(n-1) - 35*a(n-2) + a(n-3) for n > 3.
G.f.: -x*(9*x^2 - 26*x + 1) / ((x-1)*(x^2 - 34*x + 1)). (End)
a(n) = c*k^n + 1/2 + o(1) with k = 17+sqrt(288) = 33.97... and c = 17/4 - sqrt(18). - Charles R Greathouse IV, May 07 2015
a(n) = (4 + 2*(17 + 12*sqrt(2))^(1-n) + (34 - 24*sqrt(2))*(17 + 12*sqrt(2))^n)/8 for n > 0. - Colin Barker, Mar 02 2016
EXAMPLE
a(2) = 9 because 9 = 3^2 = 1001_2 and 100_2 = 4 = 2^2.
MATHEMATICA
LinearRecurrence[{35, -35, 1}, {0, 1, 9, 289}, 25] (* Paolo Xausa, Jul 22 2024 *)
PROG
(PARI) concat(0, Vec(-x*(9*x^2-26*x+1)/((x-1)*(x^2-34*x+1)) + O(x^100))) \\ Colin Barker, Sep 15 2014
(PARI) is(n)=issquare(n) && issquare(n\2) \\ Charles R Greathouse IV, May 07 2015
CROSSREFS
KEYWORD
nonn,base,easy
AUTHOR
Henry Bottomley, Jul 14 2000
STATUS
approved