cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

A031149 Numbers whose square with its last digit deleted is also a square.

Original entry on oeis.org

0, 1, 2, 3, 4, 7, 13, 16, 19, 38, 57, 136, 253, 487, 604, 721, 1442, 2163, 5164, 9607, 18493, 22936, 27379, 54758, 82137, 196096, 364813, 702247, 870964, 1039681, 2079362, 3119043, 7446484, 13853287, 26666893, 33073696, 39480499, 78960998, 118441497, 282770296
Offset: 1

Views

Author

Keywords

Comments

Square root of A023110(n).
For the first 4 terms, the square has only one digit, but in analogy to the sequences in other bases (A204502, A204512, A204514, A204516, A204518, A204520, A004275, A055793, A055792), it is understood that deleting this digit yields 0.
From Robert Israel, Feb 16 2016: (Start)
Solutions x to x^2 = 10*y^2 + j, j in {0,1,4,6,9}, in increasing order of x.
j=0 occurs only for x=0.
Let M be the 2 X 2 matrix [19, 60; 6, 19].
Solutions of x^2 = 10*y^2 + 1 are (x,y)^T = M^k (1,0)^T for k >= 0.
Solutions of x^2 = 10*y^2 + 4 are (x,y)^T = M^k (2,0)^T for k >= 0.
Solutions of x^2 = 10*y^2 + 6 are (x,y)^T = M^k (4,1)^T and M^k (16,5)^T for k >= 0.
Solutions of x^2 = 10*y^2 + 9 are (x,y)^T = M^k (3,0)^T, M^k (7,2)^T, M^k (13,4)^T for k >= 0.
Since (1,0)^T <= (2,0)^T <= (3,0)^T <= (4,1)^T <= (7,2)^T <= (13,4)^T <= (16,5)^T <= (19,6)^T = M (1,0)^T (element-wise) and M has positive entries, we see that the terms always occur in this order, for successive k.
The eigenvalues of M are 19 + 6*sqrt(10) and 19 - 6*sqrt(10).
From this follow my formulas below and the G.f. (End)

Examples

			364813^2 = 133088524969, 115364^2 = 13308852496.
		

References

  • R. K. Guy, Neg and Reg, preprint, Jan 2012. [From N. J. A. Sloane, Jan 12 2012]

Crossrefs

Programs

  • Maple
    for i from 1 to 150000 do if (floor(sqrt(10 * i^2 + 9)) > floor(sqrt(10 * i^2))) then print(floor(sqrt(10 * i^2 + 9))) end if end do;
  • Mathematica
    fQ[n_] := IntegerQ@ Sqrt@ Quotient[n^2, 10]; Select[ Range[ 0, 40000000], fQ] (* Harvey P. Dale, Jun 15 2011 *) (* modified by Robert G. Wilson v, Jan 16 2012 *)
  • PARI
    s=[]; for(n=0, 1e7, if(issquare(n^2\10), s=concat(s,n))); s \\ Colin Barker, Jan 17 2014; typo fixed by Zak Seidov, Jan 31 2014

Formula

Appears to satisfy: a(n)=38a(n-7)-a(n-14) which would require a(-k) to look like 3, 2, 1, 4, 7, 13, 16, 57, 38, 19, 136, ... for k>0. - Henry Bottomley, May 08 2001
Empirical g.f.: x^2*(1 + 2*x + 3*x^2 + 4*x^3 + 7*x^4 + 13*x^5 + 16*x^6 - 19*x^7 - 38*x^8 - 57*x^9 - 16*x^10 - 13*x^11 - 7*x^12 - 4*x^13) / ((1 - 38*x^7 + x^14)). - Colin Barker, Jan 17 2014
a(n) = 38*a(n-7) - a(n-14) for n>15 (conjectured). - Colin Barker, Dec 31 2017
With e1 = 19 + 6*sqrt(10) and e2 = 19 - 6*sqrt(10),
a(2+7k) = (e1^k + e2^k)/2,
a(3+7k) = e1^k + e2^k,
a(4+7k) = (3/2) (e1^k + e2^k),
a(5+7k) = (2+sqrt(10)/2) e1^k + (2-sqrt(10)/2) e2^k,
a(6+7k) = (7/2+sqrt(10)) e1^k + (7/2-sqrt(10)) e2^k,
a(7+7k) = (13/2+2 sqrt(10)) e1^k + (13/2-2 sqrt(10)) e2^k,
a(8+7k) = (8+5 sqrt(10)/2) e1^k + (8-5 sqrt(10)/2) e2^k. - Robert Israel, Feb 16 2016

Extensions

4 initial terms added by M. F. Hasler, Jan 15 2012
a(40) from Robert G. Wilson v, Jan 15 2012

A031150 Appending a digit to n^2 gives another perfect square.

Original entry on oeis.org

1, 2, 4, 5, 6, 12, 18, 43, 80, 154, 191, 228, 456, 684, 1633, 3038, 5848, 7253, 8658, 17316, 25974, 62011, 115364, 222070, 275423, 328776, 657552, 986328, 2354785, 4380794, 8432812, 10458821, 12484830, 24969660, 37454490
Offset: 1

Views

Author

Keywords

Comments

Square root of 'Squares from A023110 with last digit removed'.
One could include an initial '0', and even list it with multiplicity 3 or 4, since 00, 01, 04 and 09 are all perfect squares: In analogy to corresponding sequences for other bases, this sequence could be defined as sqrt(floor[A023110/10]), see A204512 [base 8], A204517 (base 7), A204519 (base 6), A204521 [base 5], A001353 [base 3], A001542 [base 2]. (For bases 4 and 9, the corresponding sequence contains all integers.) - M. F. Hasler, Jan 16 2012

Examples

			5^2 = 25 and 16^2 = 256, so 5 is in the sequence.
115364^2 = 13308852496, 364813^2 = 133088524969.
		

References

  • R. K. Guy, Neg and Reg, preprint, Jan 2012.

Crossrefs

See A202303 for the resulting squares.

Programs

  • Maple
    for i from 1 to 150000 do if (floor(sqrt(10 * i^2 + 9)) > floor(sqrt(10 * i^2))) then print(i) end if end do;
  • Mathematica
    CoefficientList[Series[(x^10 + 2 x^9 + 4 x^8 + 5 x^7 + 18 x^6 + 12 x^5 + 6 x^4 + 5 x^3 + 4 x^2 + 2 x + 1)/(x^14 - 38 x^7 + 1), {x, 0, 50}], x] (* Vincenzo Librandi, Oct 19 2013 *)
    LinearRecurrence[{0,0,0,0,0,0,38,0,0,0,0,0,0,-1},{1,2,4,5,6,12,18,43,80,154,191,228,456,684},40] (* Harvey P. Dale, Jun 09 2017 *)

Formula

G.f.: x*(x^10+2*x^9+4*x^8+5*x^7+18*x^6+12*x^5+6*x^4+5*x^3+4*x^2+2*x+1) / (x^14-38*x^7+1). - Colin Barker, Jan 30 2013

A030686 Smallest nontrivial extension of n^2 which is a square.

Original entry on oeis.org

16, 49, 900, 169, 256, 361, 4900, 6400, 8100, 10000, 12100, 1444, 16900, 19600, 22500, 25600, 28900, 3249, 36100, 40000, 44100, 48400, 52900, 57600, 62500, 67600, 72900, 78400, 84100, 90000, 96100, 102400, 108900, 115600, 122500
Offset: 1

Views

Author

Keywords

Comments

Nontrivial extension means appending at least one digit even if the number is already a square.

Crossrefs

See also A023110 = A031149^2 and A202303 = A031150^2 for a related concept, and cross-references there (and in links) for the analog in bases other than 10. - M. F. Hasler, Sep 28 2014

Formula

a(n) = A030687(n)^2. - M. F. Hasler, Sep 28 2014
a(n) = A030666(n^2). - Alonso del Arte, Apr 01 2020
Showing 1-3 of 3 results.