%I#33 2021年9月22日13:15:03

%S 0,0,0,1,0,0 0,0，0,1,0，0，0,0，

%温度1,0,1,1,0,0,1,10,0,0，0,00,01,0,02,0,1,0,1,0，

%U 0,1,0,0,1,1,1,1,1,0,0,0，1,0,1,0，0,01,0,0

%N素数=1模3除以N。

%C a（n）=2的第一个例子是n=91；a（n）=3的第一个实例是n＝1729。1729是著名的拉马努扬出租车号码——见A001235_Harvey P.Dale_，2013年6月25日

%H Antti Karttunen，n的表，n=1..10000的a（n）</a>

%H S.R.Finch和Pascal Sebah，<a href=“http://arXiv.org/abs/math.NT/0604465“>平方和立方模n，arXiv:math/0604465[math.NT]，2006-2016。

%F如果p=1（mod 3），则a（p^e）=1的加法，否则为0。

%F From _Antti Karttunen，2017年7月10日：（开始）

%F a（1）=0；对于n>1，（（A020639（n）mod 3）mod 2）+a（A028234（n））。

%F a（n）=A001221（n）-A005090（n）-A079978（n）。

%F（结束）

%p A005088:=程序（n）

%p局部a，pe；

%p a：=0；

%ifactors（n）[2]do中pe的p

%p如果modp（op（1，pe），3）=1，则

%pa:=a+1；

%p end if；

%p端do：

%p a；

%p end程序：#_R.J.Mathar_，2020年5月19日

%t连接[{0}，表[Count[Transpose[FactorInteger[n]][[1]，_？（型号[#-1,3]==0&）]，{n，2100}]]（*哈维·P·戴尔，2021年9月22日*）

%t数组[DivisorSum[#，1&，And[PrimeQ@#，Mod[#，3]==1]&]&，91]（*_Michael De Vlieger_，2017年7月11日*）

%o（PARI）a（n）=我的（f=系数（n）[，1]）；总和（i=1，#f，f[i]%3==1）\\查尔斯·格里特豪斯IV_，2017年1月16日

%o（方案）（定义（A005088 n）（如果（=1 n）0（+（模数（A020639 n）3）2）（A0050088（A028234 n）））_Antti Karttunen，2017年7月10日

%Y参见A001221、A005070、A005090、A020639、A028234、A079978。

%Y参考A121940（第一个数字有n个此类因素）。

%K nonn公司

%O 1,91型

%A _N.J.A.斯隆_