|
|
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
|
|
|
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
|
|
CROSSREFS
|
|
|
KEYWORD
|
easy,nonn,base
|
|
AUTHOR
|
C. Ronaldo (aga_new_ac(AT)hotmail.com), Jan 15 2005
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|