fa[n_]：=fa[n]=系数整数[n]；tes1[n]：=并集@表[IntegerQ[n/（fa[n][[i，1]]^fa[n][[i、2]]-1）]，{i，长度[fa[n4]}]=={True}；G[n]：=并集@表[Mod[n/fa[n][[i，1]]^（fa[n][[i，2]]），fa[n][[i，1]]]==1，{i，长度[fa[n]]}]=={True}；选择[Range[40000]，tes1[#]和&G[#]&]
（*或*）
ok[n_]：=n==1||块[{f=因子整数@n，p}，p=功率@@@f；和@@IntegerQ/@（n/（p-1））&&Union@Mod[n/p，First/@f]=={1}]；选择[Range@80000，ok](*乔瓦尼·雷斯塔2017年3月10日*）
