%I#21 2022年2月9日09:04:56
%S 2,3,2,5,7,7,2,3,3,11,7,13,11,4,2,7,7,19,3,5,4,23,7,5,17,11,29,13,
%电话31,2,7,5,6,7,37,23,8,3,41,4,43,4,4,3,47,7,3,10,17,53,7,8,11,7,
%U 59、13、61、6、5、2、9、19、67、5、13、17、71、7、73、41、4、23、9、6、79、3、4、83、4、11、47
%N对于有理数p/q,设f(p/q)=p+q的不同素数因子(A008472)之和除以p+q不同素数因数(A001221)的个数;a(n)是通过迭代f获得的,从n/1开始,直到达到整数,或者如果没有达到整数,则a(n)=0。
%C A247462给出了达到a(n)所需的迭代次数_Reinhard Zumkeller,2014年9月17日
%H Reinhard Zumkeller,n的表,n=1..10000的a(n)</a>
%H P.Schogt,<a href=“http://arxiv.org/abs/1211.6583“>野生数字问题:数学还是虚构?</a>,arXiv预印本arXiv:1211.6583[math.HO],2012.-发件人:N.J.A.Sloane,2013年1月3日
%e f(5/1)=5/2,f(5/2)=7,因此a(5)=7。
%t nxt[n_]:=模块[{s=分子[n]+分母[n]},总计[Transpose[FactorInteger[s]][[1]]/PrimeNu[s]];表[NestWhile[nxt,nxt[n],!整数Q[#]&],{n,90}](*哈维·P·戴尔,2013年3月15日*)
%o(哈斯克尔)
%o导入数据。比率(%),分子,分母)
%o a058977=分子。直到(==1)。分母)f。f、。来自Integral
%o其中f x=a008472 z%a001221 z
%o其中z=分子x+分母x
%o——Reinhard Zumkeller,2014年8月29日
%o(PARI)f2(p,q)=我的(f=系数(p+q)[,1]~);vecsum(f)/#f;
%o f1(r)=f2(分子(r),分母(r));
%o循环(列表)={my(v=Vecrev(列表));对于(i=2,#v,如果(v[i]==v[1],返回(1)););}
%o a(n)={my(ok=0,m=f2(n,1),list=list());while(分母(m)!=1,m=f1(m);listput(列表,m);if(循环(列表),return(0););return
%Y参考A058971、A058972。
%Y参考A008472、A001221。
%Y参见A247462和A247468。
%不,简单,好
%O 1,1号机组
%A _N.J.A.Sloane,2001年1月14日
%E更多来自马修·康罗伊的条款,2001年4月18日