login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A102567 Numbers k such that the concatenation of k with itself is a biperiod square. 49
13223140496, 20661157025, 29752066116, 40495867769, 52892561984, 66942148761, 82644628100, 183673469387755102041, 326530612244897959184, 510204081632653061225, 734693877551020408164 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
Also, numbers N associated with A106497.
Also, numbers k such that k concatenated with k-1 gives the product of two numbers which differ by 2. E.g., 13223140496//13223140495 = 36363636363 * 36363636365, where // denotes concatenation. - Giovanni Resta and Franklin T. Adams-Watters, Nov 13 2006
REFERENCES
Andrew Bridy, Robert J. Lemke Oliver, Arlo Shallit, and Jeffrey Shallit, The Generalized Nagell-Ljunggren Problem: Powers with Repetitive Representations, Experimental Math, 28 (2019), 428-439.
R. Ondrejka, Problem 1130: Biperiod Squares, Journal of Recreational Mathematics, Vol. 14:4 (1981-82), 299. Solution by F. H. Kierstead, Jr., JRM, Vol. 15:4 (1982-83), 311-312.
LINKS
Dr Barker, Can Numbers Like These Be Square?, YouTube video, 2023.
Andrew Bridy, Robert J. Lemke Oliver, Arlo Shallit, and Jeffrey Shallit, The Generalized Nagell-Ljunggren Problem: Powers with Repetitive Representations, preprint arXiv:1707.03894 [math.NT], July 14 2017.
EXAMPLE
13223140496 concatenated with 13223140496 is 1322314049613223140496 = 36363636364^2.
40495867769 is in the sequence because writing it twice gives the square number 4049586776940495867769 = 63636363637^2.
MAPLE
with(numtheory): Digits:=50:for d from 1 to 35 do tendp1:=10^d+1: tendp1fact:=ifactors(tendp1)[2]: n:=mul(piecewise(tendp1fact[i][2] mod 2=1, tendp1fact[i][1], 1), i=1..nops(tendp1fact)):for i from ceil(sqrt((10^(d-1))/n)) to floor(sqrt((10^d-1)/n)) do printf("%d, ", n*i^2) od: od:
MATHEMATICA
A102567L[n_] := Catenate@Table[Module[{fac = FactorInteger[10^k + 1], min}, If[Max@fac[[All, -1]] == 1, {}, min = Times @@ Cases[fac, {a_, _?OddQ} :> a]; Table[min s^2, {s, Ceiling@Sqrt[10^(k - 1)/min], Floor@Sqrt[(10^k - 1)/min]}]]], {k, n}]; A102567L[30] (* JungHwan Min, Dec 11 2016 *)
A102567Q = IntegerQ@Sqrt@FromDigits[Join[#, #] &@IntegerDigits[#]] & (* JungHwan Min, Dec 11 2016 *)
PROG
(Python)
from itertools import count, islice
from sympy import sqrt_mod
def A102567_gen(): # generator of terms
for j in count(0):
b = 10**j
a = b*10+1
for k in sorted(sqrt_mod(0, a, all_roots=True)):
if a*b <= k**2 < a*(a-1):
yield k**2//a
A102567_list = list(islice(A102567_gen(), 10)) # Chai Wah Wu, Feb 19 2024
CROSSREFS
Sequence in context: A204096 A113639 A059122 * A016932 A016992 A017160
KEYWORD
easy,nonn,base
AUTHOR
C. Ronaldo (aga_new_ac(AT)hotmail.com), Jan 15 2005
EXTENSIONS
Entry revised by N. J. A. Sloane, Nov 14 2006 and also Nov 27 2006
Definition edited and reference added by William Rex Marshall, Nov 12 2010
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 18 22:56 EDT 2024. Contains 370952 sequences. (Running on oeis4.)