nn=120;rad[x_]:=rad[x]=倍@FactorInteger[x][[All,1]];
f[x_]:=倍@@Prime@位置[Reverse@IntegerDigit[x,2],1][[All,1];
c[_]:=错误;c[1]=正确;q[_]:=1;a[0]=1;r[_]:=1;
执行[If[PrimeQ[#],
而[c[集[k,#^q[#]]],q[#]++],
而[Or[c[r[#]]!可除[#,rad[r[#]]],r[#]++];k=r[#]]&[f[i]];集合[{a[i],c[k]},{k,True}],{i,nn}];
数组[a,nn+1,0]
|