%I#67 2023年3月31日14:14:35

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

%温度0,0,1,0,0,0，1,0,1,0，0,01,0,1,1,0,10,0:0,0,1，0,0，

%U 0,0,1,0,10,0,00,0.0,0,1,0

%如果N是素数幂p^k（k>=0），则为N1，否则为0。

%C单位或素数幂的特征函数p^k（k>=1）。素数幂的特征函数p^k（k>=0）_Daniel Forgues_，2009年3月3日

%C部分金额见A065515_Reinhard Zumkeller_，2009年11月22日

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

%特征函数的索引项</a>

%F Dirichlet生成函数：1+ppzeta（s）。这里ppzeta（s）=和{p素数}和{k>=1}1/（p^k）^s_富兰克林·亚当斯·沃特斯，2005年9月11日

%F a（n）=0^（A119288（n）-1）_Reinhard Zumkeller_，2006年5月13日

%F a（A000961（n））=1；a（A024619（n））=0.-_Reinhard Zumkeller_2011年11月17日

%F a（n）=如果A001221（n）<=1，则为1，否则为0.-_Reinhard Zumkeller_，2015年11月28日

%F如果n>=2，a（n）=A069513（n）.-_Jeppe Stig Nielsen，2016年2月2日

%F猜想：对于所有n>1.-，a（n）=（n-A048671（n））/A000010（n）_韦林·亚涅夫（Velin Yanev），2021年3月10日【猜测是真的】

%p A010055:=程序（n）

%p如果n=1，则

%p 1；

%p其他

%p如果nops（ifactors（n）[2]）=1，则

%第1页；

%p其他

%p 0；

%p end if；

%p结束if；

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

%t A010055[n_]：=布尔[PrimeNu[n]<=1]；A010055/@Range[20]（*_Enrique Pérez Herrero_，2011年5月30日*）

%t{1}~加入~表格[Boole@PrimePowerQ@n，{n，2，105}]（*_Michael De Vlieger_，2016年2月2日*）

%o（PARI）用于（n=1120，打印1（ω（n）<=1，“，”））

%o（哈斯克尔）

%o a010055 n=如果a001221 n<=1，则1为0

%o--_Reinhard Zumkeller_，2015年11月28日，2013年3月19日，2011年11月17日

%o（Python）

%o来自sympy导入因子

%o定义A010055（n）：返回int（len（primefactors（n））<=1）#_Chai Wah Wu_，2023年3月31日

%Y参考A069513（如果n是素数幂p^k（k>=1），则为1，否则为0。）

%Y参考A268340。

%Y参考A100995。

%Y参考A001221、A000961、A024619、A065515、，

%K nonn公司

%O 1,1号机组

%A _N.J.A.斯隆_

%E更多条款来自_Charles R Greathouse IV_，2008年3月12日

%E由_Daniel Forgues_编辑，2009年3月2日

%关于Galois油田的E Comment由_Franklin T.Adams-Waters_移至A069513，2009年11月2日