OFFSET
1,1
COMMENTS
There are 5 primes in the first 6100 terms of this sequence, see A283561.
LINKS
Robert Israel, Table of n, a(n) for n = 1..365
Eric Weisstein's World of Mathematics, Consecutive Number Sequences
MAPLE
A[1]:= 2:
x:= 2:
for n from 2 to 20 do
if issqr(x+1) then x:= x+2 else x:= x+1 fi;
A[n]:= A[n-1]*10^(1+ilog10(x))+x;
od:
seq(A[n], n=1..20); # Robert Israel, Mar 13 2017
MATHEMATICA
cns[n_]:=FromDigits[Flatten[IntegerDigits[Table[k+Floor[1/2+Sqrt[k]], {k, 1, n}]]]]
Module[{nn=20, sqf}, sqf=Select[Range[2, nn], !IntegerQ[Sqrt[#]]&]; Table[FromDigits[Flatten[ IntegerDigits/@ Take[sqf, n]]], {n, Length[sqf]}]] (* Harvey P. Dale, Sep 07 2024 *)
PROG
(PARI) a(n)=my(s=""); for(k=1, n, s=Str(s, (sqrtint(4*k)+1)\2 + k)); eval(s) \\ Charles R Greathouse IV, Mar 10 2017
CROSSREFS
KEYWORD
base,nonn,changed
AUTHOR
XU Pingya, Mar 10 2017
STATUS
approved