login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A204516
Numbers such that floor(a(n)^2 / 7) is a square.
20
0, 1, 2, 3, 8, 16, 45, 127, 254, 717, 2024, 4048, 11427, 32257, 64514, 182115, 514088, 1028176, 2902413, 8193151, 16386302, 46256493, 130576328, 261152656, 737201475, 2081028097, 4162056194, 11748967107, 33165873224, 66331746448
OFFSET
1,3
COMMENTS
Or: Numbers whose square, with its last base-7 digit dropped, is again a square (where for the first 3 terms, dropping the digit is meant to yield zero).
FORMULA
G.f. = (x + 2*x^2 + 3*x^3 - 8*x^4 - 16*x^5 - 3*x^6 )/(1 - 16*x^3 + x^6).
floor(a(n)^2 / 7) = A204517(n)^2.
MATHEMATICA
LinearRecurrence[{0, 0, 16, 0, 0, -1}, {0, 1, 2, 3, 8, 16, 45}, 30] (* or *) CoefficientList[Series[ (x+2x^2+3x^3-8x^4-16x^5-3x^6)/(1-16x^3+x^6), {x, 0, 30}], x] (* Harvey P. Dale, Apr 22 2023 *)
PROG
(PARI) b=7; for(n=0, 2e9, issquare(n^2\b) & print1(n", "))
CROSSREFS
Cf. A031149 (base 10), A204502 (base 9), A204514 (base 8), A204518 (base 6), A004275=sqrt(A055808) (base 4), A001075=sqrt(A055793) (base 3), A001541=sqrt(A055792) (base 2).
Sequence in context: A051573 A292853 A363204 * A331679 A277346 A005648
KEYWORD
nonn
AUTHOR
M. F. Hasler, Jan 15 2012
STATUS
approved