OFFSET
1,1
COMMENTS
Obviously b(n) = 100^n - 10^n + 1 = (91, 9901, 999001, 99990001, ...) is a subsequence. Are { b(2), b(4), b(6), b(8) } the only terms of this sequence that are prime? - M. F. Hasler, Mar 30 2008. Answer: The smallest prime in this sequence that is not of the form b(n) is A054216(155) = 811451682377384625400019885321 [Max Alekseyev, Oct 08 2008]. See A145381 for further prime terms.
Other subsequences are c(n) = ( 10^(6n) - 2*10^(5n) - 10^(3n) - 2*10^n + 1 )/3 (n>=2), d(n) = (33/101)*(100^(404n+71)+1)+10^(404n+71) (n>=0) and e(n) = (33/101)*(100^(404n-71)+1)+10^(404n-71) (n>=1). Primes among these include c(10), c(14) and d(0). - M. F. Hasler, Oct 09 2008
A positive integer m is in this sequence if and only if m^2 == -1 (mod 10^k + 1) where k is the number of decimal digits in m. Note that k cannot be odd, since in this case 11 divides 10^k + 1 while -1 is not a square modulo 11. - Max Alekseyev, Oct 09 2008
REFERENCES
Luca, Florian, and Pantelimon Stănică. "Perfect Squares as Concatenation of Consecutive Integers." The American Mathematical Monthly 126.8 (2019): 728-734.
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..1000
FORMULA
a(n) = sqrt(A054215(n)). - Max Alekseyev, May 14 2007
EXAMPLE
'8242' + '8242-1' gives 82428241 which is 9079^2.
Leading zeros are not allowed, which is why c(1)=266327 is not in this sequence although c(1)^2 = 070930 070929.
PROG
(PARI) isA054216(n)={ 1==[1, -1]*divrem(n^2, 10^(#Str(n^2)\2)) & #Str(n^2)%2==0 }
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Patrick De Geest, Feb 15 2000
EXTENSIONS
More terms from Max Alekseyev, May 14 2007
Several corrections and additions from M. F. Hasler, Oct 09 2008
STATUS
approved