ลองกำหนดลำดับของรากที่สองจำนวนเต็ม ครั้งแรก (1) = 1 จากนั้น (n) เป็นจำนวนเต็มบวกที่เล็กที่สุดไม่ได้เห็นมาก่อนเช่นว่า

sqrt(a(n) + sqrt(a(n-1) + sqrt(... + sqrt(a(1)))))

เป็นจำนวนเต็ม ตัวอย่างบางส่วน:

a (2) คือ 3 เพราะมันเป็นจำนวนเต็มที่น้อยที่สุดที่sqrt(a(2) + sqrt(a(1))) = sqrt(a(2) + 1)เป็นจำนวนเต็มและ 3 ไม่เคยเกิดขึ้นในลำดับก่อนหน้านี้

a (3) คือ 2 เพราะมันเป็นจำนวนเต็มที่น้อยที่สุดที่sqrt(a(3) + sqrt(a(2) + sqrt(a(1)))) = sqrt(a(3) + 2)เป็นจำนวนเต็มและ 2 ไม่เคยเกิดขึ้นมาก่อนในลำดับ

a (4) คือ 7 เพราะsqrt(a(4) + 2)เป็นจำนวนเต็ม เราไม่สามารถมี (4) = 2 เพราะ 2 เกิดขึ้นแล้วในลำดับของเรา

เขียนโปรแกรมหรือฟังก์ชั่นที่กำหนดพารามิเตอร์nส่งคืนลำดับของตัวเลข a (1) ไปยัง a (n)

ลำดับเริ่มต้น 1,3,2,7,6,13,5, ....

^{แหล่งที่มาของลำดับนี้มาจากคำถาม Math.SEนี้}

เนื้อเรื่องขององค์ประกอบ 1,000 ชิ้นแรกในลำดับ:

Python 2 , 80 ไบต์

s=[]
exec'x=q=1\nwhile(x in s)+q%1:x+=1;q=(v+x)**.5\nv=q;s+=x,;'*input()
print s

ลองออนไลน์!

Haskell , 103 87 ไบต์

ไม่มีประสิทธิภาพอย่างน่ากลัว แต่ไม่พึ่งพาเลขคณิตจุดลอยตัว นี่a(x) = sqrt(f(x)+a(x-1))คือลำดับผู้ช่วยที่ทำให้การคำนวณง่ายขึ้น

a 0=0
a x=[k|k<-[1..],m<-[k^2-a(x-1)],m>0,notElem m$f<$>[1..x-1]]!!0
f x=(a x)^2-a(x-1)

ลองออนไลน์!

Python 2 , 87 ไบต์

t,=s=1,
for n in~-input()*s:
 while(n in s)+(t+n)**.5%1:n+=1
 s+=n,;t=(t+n)**.5
print s

ลองออนไลน์!

-3 ต้องขอบคุณนาย Xcoder
-5 ขอบคุณที่OVS

MATL , 30 27 ไบต์

lXHiq:"`@ymH@+X^1\+}8MXHx@h

ลองออนไลน์! หรือดูจอแสดงผลกราฟิก (ใช้เวลาสักครู่; หมดเวลาสำหรับอินพุตเกินประมาณ60)

คำอธิบาย

l          % Push 1. This is the array that holds the sequence, initialized to
           % a single term. Will be extended with subsequent terms
XH         % Copy into clipboard H, which holds the latest result of the 
           % "accumulated" square root
iq:"       % Input n. Do the following n-1 times
  `        %   Do...while
    @      %     Push interaton index k, starting at 1. This is the candidate
           %     to being the next term of the sequence
    y      %     Push copy of array of terms found so far
    m      %     Ismbmer? True if k is in the array
    H      %     Push accumulated root
    @+     %     Add k
    X^     %     Square root
    1\     %     Modulo 1. This gives 0 if k gives an integer square root
    +      %     Add. Gives nonzero if k is in the array or doesn't give an
           %     integer square root; that is, if k is invalid.
           %   The body of the do...while loop ends here. If the top of the
           %   stack is nonzero a new iteration will be run. If it is zero that
           %   means that the current k is a new term of the sequence
  }        %   Finally: this is executed after the last iteration, right before
           %   the loop is exited
    8M     %     Push latest result of the square root
    XH     %     Copy in clipboard K
    x      %     Delete
    @      %     Push current k
    h      %     Append to the array
           % End do...while (implicit)
           % Display (implicit)

Mathematica, 104 ไบต์

(s=f={i=1};Do[t=1;While[!IntegerQ[d=Sqrt[t+s[[i]]]]||!f~FreeQ~t,t++];f~(A=AppendTo)~t;s~A~d;i++,#-1];f)&  

ลองออนไลน์!

ลำดับของรากที่สองนั้นน่าสนใจมาก ...
และแสดงรูปแบบที่คล้ายกัน

1,2,2,3,3,4,3,5,3,6,4,4,5,4,6,5,5,6,6,7,4,7,5,7,6, 8,4,8,5,8,6,9,5,9,6,10,5,10,6,11,5,11,6,12,6,13,6,14,7,7, 8,7,9,7,10,7,11,7,12,7,13,7,14,8,8,9,8,10 ...

นี่คือความแตกต่างของลำดับหลัก

Python 2 , 117 115 112 102 99 87 ไบต์

t,=r=1,;exec"x=1\nwhile(t+x)**.5%1or x in r:x+=1\nr+=x,;t=(t+x)**.5;"*~-input();print r

ลองออนไลน์!

ใช้t=(t+x)**.5ตรรกะจากคำตอบของ Erik

JavaScript (ES7), 89 82 77 76 ไบต์

i=>(g=k=>(s=(++n+k)**.5)%1||u[n]?g(k):i--?[u[n]=n,...g(s,n=0)]:[])(n=0,u=[])

การสาธิต

let f =

i=>(g=k=>(s=(++n+k)**.5)%1||u[n]?g(k):i--?[u[n]=n,...g(s,n=0)]:[])(n=0,u=[])

console.log(JSON.stringify(f(10)))

จัดรูปแบบและแสดงความคิดเห็น

i => (                             // given i = number of terms to compute
  u = [],                          // u = array of encountered values
  g = p =>                         // g = recursive function taking p = previous square root
    (s = (++n + p) ** .5) % 1      // increment n; if n + p is not a perfect square,
    || u[n] ?                      // or n was already used:
      g(p)                         //   do a recursive call with p unchanged
    :                              // else:
      i-- ?                        //   if there are other terms to compute:
        [u[n] = n, ...g(s, n = 0)] //     append n, set u[n] and call g() with p = s, n = 0
      :                            //   else:
        []                         //     stop recursion
  )(n = 0)                         // initial call to g() with n = p = 0

R , 138 105 99 99 ไบต์

function(n){for(i in 1:n){j=1
while(Reduce(function(x,y)(y+x)^.5,g<-c(T,j))%%1|j%in%T)j=j+1
T=g}
T}

ลองออนไลน์!

-33 ไบต์โดยใช้เคล็ดลับที่ฉลาดของ Tfeldsqrt()%%1ในขณะที่ลูป

-6 ไบต์โดยใช้ T แทน F

คำตอบเดิม 138 ไบต์:

function(n,l={}){g=function(L)Reduce(function(x,y)(y+x)^.5,L,0)
for(i in 1:n){T=1
while(g(c(l,T))!=g(c(l,T))%/%1|T%in%l)T=T+1
l=c(l,T)}
l}

ลองออนไลน์!

Husk , 21 ไบต์

!¡oḟȯΛ±sFo√+Som:`-N;1

ลองออนไลน์!

อย่างไร?

!¡oḟȯΛ±sFo√+Som:`-N;1    Function that generates a list of prefixes of the sequence and indexes into it
                   ;1    The literal list [1]
 ¡                       Iterate the following function, collecting values in a list
  oḟȯΛ±sFo√+Som:`-N        This function takes a prefix of the sequence, l, and returns the next prefix.
                `-N      Get all the natural numbers that are not in l.
            Som:         Append l in front each of these numbers, generates all possible prefixes.
    ȯΛ±sFo√+               This predicate tests if sqrt(a(n) + sqrt(a(n-1) + sqrt(... + sqrt(a(1))))) is an integer.
        F                Fold from the left
         o√+             the composition of square root and plus
       s                 Convert to string
    ȯΛ±                  Are all the characters digits, (no '.')
  oḟ                     Find the first list in the list of possible prefixes that satisfies the above predicate
!                        Index into the list