OFFSET
1,1
EXAMPLE
a(3) = 67 is a term because it is prime and 67^2 = 4489 is the concatenation of 44 and 2*44+1=89.
MAPLE
dcat:=proc(a, b) a*10^(1+ilog10(b))+b end proc:
f:= proc(t) local s;
if not issqr(t) then return NULL fi;
s:=sqrt(t);
if isprime(s) then return s fi
end proc:
map(f, [seq(dcat(x, 2*x+1), x=1..5*10^7)]);
CROSSREFS
KEYWORD
nonn,base
AUTHOR
J. M. Bergot and Robert Israel, Jul 21 2022
EXTENSIONS
More terms from Jinyuan Wang, Jul 21 2022
STATUS
approved