%I#16 2023年12月10日17:37:14
%S 0,1,1,2,1,1,1,1,2,1,1,2,1,1,2,1,1,1,1,4,1,2,2,1,1,2,1,2,2,1,1,1,4,1,1,
%在1,1,1,1,1,2,1,1,2,1,1,2,1,4,2,1,2,1,2,1,2,1,2,1,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,1,1,2,2,1,1,1,1,2,2,2,2,1,1,1,2,
%U 1,1,1,1,1,1,2,1,1,1,4,1,1,2,1,1,2,1,1,1,1,1,2,1,2,1,1,1,1,4,1,2,2,1,1
%当从N开始并用A334870迭代时,第一个费米-狄拉克因子(形式为p^(2^k),A050376)的素数达到N,其中a(1)=0。
%H Antti Karttunen,n表,n=1..65537的a(n)</a>
%e对于n=27,当用A334870迭代时,我们得到路径27->18->9,其中9是遇到的第一个素数幂,其指数形式为2^k,为9=3^2,因此a(27)=2。查看二叉树A334860或A334866了解这些路径的运行方式。
%e对于n=900,当使用A334870迭代时,我们得到路径900->30->15->10->5,5^1最终是一个素数幂,其指数为2的幂,因此a(900)=1。请注意,900是该序列中第一个未在A333634中列出的位置1。
%o(PARI)
%o A209229(n)=(n位和(n,n-1));
%o A302777(n)=A209229(i功率(n));
%o A334870(n)=如果(issquare(n),sqrtint(n)),my(c=芯(n)、m=n);对于素数(p=2,如果(!(c%p),m/=p;断裂,m*=p));(m) );
%o A335428(n)={而(n>1&&!A302777(n),n=A334870(n));isprimpower(n);};
%o(PARI)
%o速度更快,A209229和A302777如上图所示:
%o A335428(n)=如果(1==n,0,while(!A302777(n),if(issquarefree,n),return(1));如果(发行方(n),n=平方(n)、n/=核心(n));isprimpower(n));
%Y参见A050376、A209229、A302777、A334860、A33486、A334870和A334872。
%K nonn公司
%O 1,4型
%2020年6月26日,安蒂·卡图宁