%I#20 2016年10月4日13:34:29
%S 0,0,0,1,0,1,1,0,2,0,1,1,1,2,0,1,2,2,0,2,2,2,1,1,1,2,2,
%温度1,3,0,1,2,0,1,0,3,3,0,0,1,1,2,3,02,2,3,4,0,12,2,2,2,1,0,3,
%U 3,3,0,3,3,1,0,2,3,4,2,0,2,2,3,0,1,3,12,2,2,0,1,4,4,0,12,4
%N a(N)是算法的运行次数。设置b_0=n,如果是素数或1或0,则停止;否则,设置c0=n(!=n)的最大除数;设置b1=c0-b0/c0。使用b_1运行。
%C对于大n:Sum_{k=1..n}a(k)~n*log(n)/2-n/2(推测)。
%H Yuriy Sibirmovsky,n的表,n=1..1000的a(n)</a>
%H Yuriy Sibirmovsky,<a href=“/A276771/A276771.png”>n=1..20000时sum_{k=1..n}a(k)的绘图</a>
%e对于n=14:b_0=14,不是素数或1或0。c0=7。b_1=7-2=5。5是质数。
%e简而言之:14->{7,2}->5。运行次数a(14)=1。
%t牛米=100;
%t a=表[0,{n,1,Nm}];
%t做[b0=n;
%t j=0;
%t当[PrimeQ[b0]==假&&b0=1&&b0=0,c=反向[除数[b0]];
%t b1=c[2]]-b0/c[2]];
%t b0=b1;j++];
%tα[[n]]=j,{n,1,Nm};
%时间
%o(PARI)停止(n)=(n<=1)||isprime(n);
%o a(n)={b=n;nb=0;while(!stop(b),d=除数(b);c=d[#d-1];b=c-b/c;nb++;);nb;}\\米歇尔·马库斯,2016年9月19日
%Y参考A000040,A276770。
%K nonn公司
%O 1,12个
%A_Yuriy Sibirmovsky,2016年9月17日
|