%I#22 2015年3月13日05:00:41
%S 1,578471694303933883479607867766958139162783255664111328,
%电话:2361072467213112574934426225149868841245028132331973768,
%电话:2052521249005626564663947536410504249801128210084
%N编号N,使A033493(N)/N为整数。
%C如果A033493(n)/n=M是偶数,则2*n是序列的成员,而A033492(2*n)/(2*n)=M/2+1。
%Ca(31)>10^7_Derek Orr_,2015年3月12日
%C倒数之和似乎很快收敛到1.0208…-Derek Orr_,2015年3月12日
%H Eric Weisstein的数学世界,<a href=“http://mathworld.wolfram.com/CollatzProblem.html“>Collatz问题</a>
%H维基百科,<a href=“http://en.wikipedia.org/wiki/Collatz_conjustic网站“>科拉茨猜想</a>
%H<a href=“/index/3#3x1”>与3x+1(或Collatz)问题相关的序列的索引条目</a>
%ta033493[n_]:=块[{f},f[1]=1;f[x_Integer?OddQ]:=3 x+1;f[x_Integer?EvenQ]:=x/2-1+Plus@@FixedPointList[f,n]];选择[Range[10^5],IntegerQ[a033493[#]/#]&](*_Michael De Vlieger_,2015年2月9日,在_Alonso del Arte_ at a033493*之后)
%o(PARI)Tsum(n)=s=n;而(n!=1,如果(n==Mod(0,2),n=n/2;s+=n);如果(n==Mod(1,2)&&n=1,n=3*n+1;s+=n));秒
%o表示(n=1,10^6,如果(类型(Tsum(n)/n)==“t_INT”,打印1(n,“,”))
%Y参考A033493。
%K nonn,更多,难
%O 1、2
%A _Derek Orr_,2015年2月7日
|