%I#17 2022年11月17日16:54:21

%S 1,2,5,29,4857，10734549562233501，612433，15,8，17353，12124753，

%电话：13,14,20,16,18,22,24,251201,26,412357,28556233,37,302269,39,32，

%U 35,523929,4015438641,82,45,65,345381,78,50,36,38,42,44,46,48,51,3,9,21,27,33,54,55,91

%N a（1）=1，a（2）=2；对于n>2，a（n）=未出现的最小正数，该正数与a（n-2）^2+a（n-1）^2具有公因数。

%C大多数术语都集中在a（n）=n线或其上方，导致前5000个术语中有51个固定点。然而，有些项要大得多，因为前两项的平方和是素数。

%C猜想：序列是正整数的置换。

%H Michael De Vlieger，n的表，a（n）表示n=1..10000</a>

%H Scott R.Shannon，<a href=“/A358444/A358444.png”>前10000个术语的图像，其中a（n）在n的10%以内。绿线为a（n）=n。

%H Michael De Vlieger，a（n）的对数散点图，n=1..24857。a（24858）是质数因子2345424289569907866042152579118178340801^2+24922^2的倍数。

%H Michael De Vlieger，n表，n=1..65536，a（n）</a>。

%e a（5）=4等于a（3）^2+a（4）^2=25+841=866，4是与866共用一个因子的最小未使用数。

%e a（9）=539562233501作为a（7）^2+a（8）^2=100+539562233401=539566233501，这是一个质数。

%t nn=120；c[_]=错误；q[_]=1；做[Set[{a[i]，c[i]、q[i]}、{i、True、2}]、{i，2}]；i=a[1]^2；j=a[2]^2；Do[k=i+j；s=FactorInteger[k][[All，1]]；Do[（m=q[#]；While[c[#m]，m++]；q[#]=m；如果[#m<k，k=#m]）&[s[[n]]]，{n，长度[s]}]；集合[{a[n]，c[k]，i，j}，{k，真，j，k^2}]，{n，3，nn}]；数组[a，nn]（*_Michael De Vlieger_，2022年11月17日*）

%Y参见A337136、A347594、A098550、A336957。

%K nonn公司

%O 1,2号机组

%A _斯科特·R·香农，2022年11月16日