%I#8 2018年5月18日09:32:06
%S 0,0,1,1,2,1,4,3,3,2,6,3,8,4,5,6,11,5,13,6,8,8,17,7,15,9,13,8,11,7,23,
%电话:12,14,12,17,10,27,14,18,13,31,10,33,16,19,18,37,14,33,16,25,19,42,15,
%U 31,20,29,23,48,14,50,25,30,27,38,17,55,27,36,21,59
%N非完美幂(A007916)的数量小于N且相对素数为N。
%a(21)=8个小于21且相对素数为21的正整数是{2,5,10,11,13,17,19,20}。
%t表格[Length[Select[Range[2,n],And[GCD@@FactorInteger[#][[All,2]]==1,GCD[n,#]==1]&]],{n,50}]
%o(PARI)a(n)=总和(k=2,n-1,!i幂(k)&&(gcd(n,k)==1));\\_米歇尔·马库斯,2018年5月15日
%Y参考A000005、A0000010、A000961、A001597、A005117、A0071916、A008683、A073311、A139555、A304326、A304362、A304574、A304575、A304576。
%K非n
%O 1.5
%A _Gus Wiseman_,2018年5月14日
|