OFFSET
1,1
EXAMPLE
441 is in the sequence because 441^2 = 194481 and the prefix 44 is in the middle of the decimal expansion of 194481.
MATHEMATICA
lst={}; Do[a=IntegerDigits[n^2]; b=Length[a]; c=IntegerLength[(n-Mod[n, 10])/10]; If[EvenQ[b-c]&&FromDigits[Take[a, {(b-c)/2+1, (b+c)/2}]]==(n-Mod[n, 10])/10, AppendTo[lst, n]], {n, 23, 5*10^6}]; lst
PROG
(PARI) ok(n)={my(t=n^2, d=logint(n, 10), b=(1+logint(t, 10)-d)/2 ); d>0 && frac(b)==0 && t\10^b%10^d==n\10}
{ for(n=1, 10^6, if(ok(n), print1(n, ", "))) } \\ Andrew Howroyd, Mar 02 2020
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Michel Lagneau, May 28 2014
EXTENSIONS
Terms a(23) and beyond from Andrew Howroyd, Mar 02 2020
STATUS
approved