nsfRun[n_]:=模块[{i=n},While[!SquareFreeQ[i],i++];i-n]
sfRun[n_]:=模块[{i=n},而[SquareFreeQ[i],i++];i-n]
sfBlockSearch[i_]:=模块[{searching=True,j=i,r,s},而[searching,r=nsfRun[j];s=sfRun[j+r];如果[s<2,j+=r+s,搜索=False]];j+r+s]
nsfsfPairQ[i_]:=nsfRun[i]==1&&sfRun[i+1]==1
nsfsfEndQ[i_]:=nsfRun[i]==1&&sfRun[i+1]>1
nsfsfRun[i_]:=模块[{searching=True,count,j=i,s,e},j=sfBlockSearch[j];当[searching,count=0;s=j;当[nsfsfPairQ[j],count++;j+=2];e=j;如果[count==0||!nsfsfEndQ[j],j=sfBlockSearch[j]、searching=False]];{s,e,count}]
a271145[{low_,high_},b_]:=模块[{i=低,k,k3,list=表[{},b]},而[i<=高,k=nsfsfRun[i];k3=最后[k]/2;如果[list[[k3]]=={},list[[k3]]=k[[1]]-2];i=k【2】】;列表]
a271145[{0,10000},3](*计算a(1),a(2),a
|