%I#15 2020年5月22日16:21:34
%S 3224331257776168071000001610513712935378247593751419857,
%电话:24760994084101515363264343431188137620511492430000028629151,
%电话:3913539345435425218756939577923516890224199115856201
%素因子提升到五次幂的N个数。
%F总和{k>=1}1/a(k)=zeta(5)/zeta(10)-1=A157291-1.-_Amiram Eldar,2020年5月22日
%e 7776=32*243=2^5*3^5,所以素因子2和3被提升到五次方。
%t选择[范围@41^5,并集[Last/@FactorInteger@#]=={5}&](*_Robert G.Wilson v_*)
%t休息[Select[Range[100],SquareFreeQ]^5](*_Vaclav Kotesovec_,2020年5月22日*)
%o(PARI)allpwrfact(n,p)=\所有素数因子提升到幂p{local(x,j,ln,y,flag);对于(x=4,n,y=Vec(factor(x));ln=长度(y[1]);flag=0;对于(j=1,ln),如果(y[2][j]==p,flag++););如果(flag==ln,print1(x“,”);)}
%Y A000584的适当子集。
%Y A329332第5列的非单位术语,按升序排列。
%Y参考A157291。
%K容易,不是
%O 4,1号机组
%A Clino Hilliard,2006年1月25日
%E来自_Robert G.Wilson的更多条款v_,2006年1月26日
|