login
A204514
Numbers such that floor(a(n)^2 / 8) is again a square.
20
0, 1, 2, 3, 6, 17, 34, 99, 198, 577, 1154, 3363, 6726, 19601, 39202, 114243, 228486, 665857, 1331714, 3880899, 7761798, 22619537, 45239074, 131836323, 263672646, 768398401, 1536796802, 4478554083, 8957108166, 26102926097, 52205852194, 152139002499, 304278004998, 886731088897
OFFSET
1,3
COMMENTS
Or: Numbers whose square, with its last base-8 digit dropped, is again a square. (Except maybe for the 3 initial terms whose square has only 1 digit in base 8.)
See A204504 for the squares resulting from truncation of a(n)^2, and A204512 for their square roots. - M. F. Hasler, Sep 28 2014
FORMULA
G.f. = (x^2 + 2*x^3 - 3*x^4 - 6*x^5)/(1 - 6*x^2 + x^4).
a(n) = sqrt(A055872(n)). - M. F. Hasler, Sep 28 2014
a(2n) = A001541(n-1). a(2n+1) = A003499(n-1). - R. J. Mathar, Feb 05 2020
MAPLE
A204514 := proc(n) coeftayl((x^2+2*x^3-3*x^4-6*x^5)/(1-6*x^2+x^4), x=0, n); end proc: seq(A204514(n), n=1..30); # Wesley Ivan Hurt, Sep 28 2014
MATHEMATICA
CoefficientList[Series[(x^2 + 2*x^3 - 3*x^4 - 6*x^5)/(x (1 - 6*x^2 + x^4)), {x, 0, 30}], x] (* Wesley Ivan Hurt, Sep 28 2014 *)
LinearRecurrence[{0, 6, 0, -1}, {0, 1, 2, 3, 6}, 40] (* Harvey P. Dale, Nov 23 2022 *)
PROG
(PARI) b=8; for(n=0, 1e7, issquare(n^2\b) & print1(n", "))
(PARI) A204514(n)=polcoeff((x + 2*x^2 - 3*x^3 - 6*x^4)/(1 - 6*x^2 + x^4+O(x^(n+!n))), n-1, x)
CROSSREFS
Cf. A031149=sqrt(A023110) (base 10), A204502=sqrt(A204503) (base 9), A204516=sqrt(A055859) (base 7), A204518=sqrt(A055851) (base 6), A204520=sqrt(A055812) (base 5), A004275=sqrt(A055808) (base 4), A001075=sqrt(A055793) (base 3), A001541=sqrt(A055792) (base 2).
Sequence in context: A181185 A241945 A018284 * A338594 A078344 A024498
KEYWORD
nonn,easy
AUTHOR
M. F. Hasler, Jan 15 2012
STATUS
approved