%I#10 2012年3月30日18:56:08

%S 1,2,3,4,5,6,7,8,9,10,11,12,13,15,16,17,18,19,20,22,23,24,25,27,28,31，

%电话：32,33,34,36,37,39,40,43,44,47,48,51,52,55,57,58,60,63,64,67,68,71,72，

%U 73,76,79,80,82,85,87,88,92,95,96,97,99100102103108111112119123单位

%N a（N）=最小正整数>a（N-1），而不是a（i）^2+a（j）^2+a（k）^2对于1<=i<j<k<=N。

%H Klaus Brockhaus，n表，n=1..10000时的a（n）</a>

%o（PARI）{m=70；z=300；a=[]；a2=Set（[]）；a3=Set（[]）；s=1；for（n=1，m，p=s；while（setsearch（a3，p，0），p++）；if（p>z，break）；new=vector（#a2）；for（k=1，#a2，q=eval（a2[k]）+p^2；if（#a）；对于（k=1，#a，q=a[k]^2+p^2；如果（q<=z，new[k]=q））；a2=setunion（a2，setminus（Set（new），Set（0））；a=concat（a，p）；s=p+1；print1（p，“，”）}/*Klaus Brockhaus，2008年7月3日*/

%Y参考A004440。【摘自R.J.Mathar_，2008年10月23日】

%K nonn公司

%O 1,2号机组

%百灵鸟金伯利_

%E由Klaus Brockhaus扩展，2008年7月3日