%I#19 2020年10月13日05:45:08

%S 6472915625466561176491000000715614826809752953611390625，

%电话：241375694704588185766121113379904148035889308915776594823321，

%电话：72900000088750368112914679691544804416183656252565726409301093638435187437614750104241

%素因子提升到六次方的N个数。

%H Nathaniel Johnston，n的表，a（n）表示n=1..5000</a>

%F a（n）=A005117（n+1）^6.-_Nathaniel Johnston，2011年6月21日

%F和{n>=1}1/a（n）=zeta（6）/zeta（12）-1=A269404-1.-_Amiram Eldar，2020年10月13日

%p for n from 2 to 100 do if（numtheory[issqrfree]（n））then printf（“%d，”，n^6）：fi:od:#_Nathaniel Johnston_，2011年6月21日

%t选择[范围@37^6，并集[Last/@FactorInteger@#]=={6}&]（*_Robert G.Wilson v_*）

%t选择[Range[2,37]，SquareFreeQ]^6（*_Amiram Eldar_，2020年10月13日*）

%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“，”）；）}

%A001014的Y子集。A030516的超集。

%Y A329332第6列的非单位术语按升序排列。

%Y参考A005117、A269404。

%K容易，不是

%O 1,1号机组

%A Clino Hilliard，2006年1月25日

%E更多条款摘自Robert G.Wilson v_，2006年1月26日