%I#40 2017年9月15日23:29:25

%S 1,1,1,1,6,1,18,1,3,2,1,1,1,2,2,1,2,2,2,2,211,6,2,1,3,1,2,6,2,4,2,1,1，

%温度4,2,3,1,1,2,2,3,2,2,1,1,1,1,1,1,2,12,2,2,2,1,4,2,2,2,1,5,1,1，

%U 4,2,2,2,3,1,7,2,1,2,2,6,7,1,2,2,8

%N a（N）=gcd（s1，s2），其中s1是奇数之和，s2是N的Collatz（3x+1）轨迹中偶数之和。

%C前10^6项a（n）的统计：

%C类+------+-----------------+------------+

%C||字数||

%C||这样||

%C|n|gcd（s1，s2）=n|百分比|

%C类+------+-----------------+------------+

%C|1|401614|40.16%|

%C|2|305471|30.54%|

%C|3 | 44381 | 4.44%|

%C|4|76228|7.62%|

%C|5|15966|1.60%|

%C|6|34514|3.45%|

%C|7|8969|0.90%|

%C|8|19156|1.92%|

%C|9|4941|0.49%|

%C | 10 | 12212 | 1.22%|

%C|11|3316|0.33%|

%C|12|8234|0.82%|

%C|>12|64998|6.50%|

%C类+------+-----------------+------------+

%当n趋于无穷大时，第三列的值似乎会无限振荡。

%C记录：1、6、18、21、23、93、187、560、1730、5098、10552、11060、11657、31072、32468、306770、793906、1956888、3107101、12210181等。；它们出现在第1、5、7、18、133、147、186、270、839、5090、5244、5488、23255、62132、113624、153341、793842、6849034、9321240、12210146等处。-Robert G.Wilson v_，2016年10月3日

%H Michel Lagneau，n的表，a（n）表示n=1..10000</a>

%H Robert G.Wilson v，<a href=“/A277068/A277068.txt”>第一次出现a（n）</a>

%e a（5）=6，因为5的Collatz轨迹是5->16->8->4->2->1=>s1=5+1=6，s2=16+8+4+2=30，gcd（6,30）=6。

%p nn：=10^7：

%p代表从1到100的n do：

%pm:=n:s1:=0:s2:=0:

%i的p从1到nn，而（m<>1）do：

%p如果irem（m，2）=0

%那么p

%p s2:=s2+m:m:=m/2：

%p其他

%ps1:=s1+m:m:=3*m+1：

%p fi:

%操作说明：

%p x：=gcd（s1+1，s2）：打印f（`%d，`，x）：

%日期：

%t排序规则[n_]：=嵌套WhileList[If[OddQ[#]，3#+1，#/2]&，n，#>1&]；f[n_]：=块[{c=Collatz@n}，GCD[Plus@@Select[c，OddQ]，Plus@@Set[c，EvenQ]]；数组[f，86]（*_Robert G.Wilson v_，2016年10月3日*）

%o（PARI）a（n）={my（se=0）；my（so=0））；while（n！=1，if（n%2，so+=n；n=3*n+1，se+=n，n=n/2）；）；gcd（se，so+1）；}\\米歇尔·马库斯，2016年10月3日

%Y参见A213909、A213916、A271973。

%K nonn公司

%O 1,5型

%2016年9月28日，拉格瑙市