%I#34 2020年2月25日08:05:56

%S 1,1,1,1,2,1,3,1,3,2,6,1,7,2,4,3,10,11,2,6,4,14,12,5,10,5,19,1,20，

%电话：6,11,7,15,3,25,8,14,6,28,2,29,8,12,10,32,3,28,7,19,11,37,4,26,10,22，

%U 14,42,2,43,14,20,15,32,5,48,15,27,8,51,6,52,17,41,6,57,12,33,20

%N在N的约化残渣体系中非素数（复合数和1）的个数。

%C在n=20、21…时与A039776不同。。。

%H Michael De Vlieger，n的表，n的a（n）=1..10000</a>

%H Abhijit A J，A.Satyanarayana Reddy，<A href=“https://arxiv.org/abs/1907.09908“>模数n的单位集合中的非犯罪数量，arXiv:1907.09908[math.GM]，2019。

%H Abhijit A.J.和A.Satyanarayana Reddy，<A href=“http://www.indianmathsociety.org.in/mathstudent-part-1-2019.pdf#page=153“>模n单位集合中的非素数，《数学学生》，第88卷，第1-2期（2019年），第147-152页。

%F a（n）=A036997（n）+1.-_Peter Luschny_，2010年10月22日

%F a（n）=A000010（n）-（A000720（n）-A001221（n））。

%在n=10时，我们看到10以下的数与10互素是1，3，7，9。去掉3和7这两个素数，我们只剩下两个数字，1和9。因此a（10）=2。

%e当n=100时，φ（100）=40，φ（100-（pi（100）-A001221（100））=17，因此a（100）=17。

%p A048864:=n->nops（选择（k->gcd（k，n）=1，删除（isprime，[$1..n]））；#_Peter Luschny2010年10月22日

%t阵列[EulerPhi@#-（PrimePi@#-PrimeNu@#）&，82]（*_Michael De Vlieger_，2016年7月3日*）

%t表[长度[Select[Range[n]，GCD[n，#]==1&&Not[PrimeQ[#]]&]]，{n，80}]（*_Alonso del Arte_，2017年10月2日*）

%o（PARI）a（n）=eulerphi（n）-（素数（n）-omega（n））；\\_Indranil Ghosh，2017年4月27日

%o（Python）

%o从sympy输入到pient，primepi，primefactors

%o定义a（n）：返回值（n）-（primepi（n）-len（primefactors（n）））#_Indranil Ghosh，2017年4月27日

%Y参见A039776、A000010、A000720、A001221、A037228、A072022、A072023、A074915。

%K nonn公司

%O 1,5型

%Abos Elemer（_L）_

%E将第二个公式转换为方程式，在示例中添加逗号-R.J.Mathar_，2010年10月23日