编辑
经核准的
P： =proc（n）局部a，x；a： =1；x： =0；而a^2<trunc（n/2）
如果frac（sqrt（n-a^2））=0，则x:=x+1；fi；a： =a+1；od；x；结束时间：
seq（P（i），i=1..100）#保罗·拉瓦2018年3月12日
检验过的
提出
return-int（not（any（（e-1 if p==2 else e）&1 for p，e in f）or n&1））+（（（m:=prod（1 if p==2其他（e+1 if p&3==1其他（e+1）&1）for p，e in f））+#柴华武2022年9月8日
（Python）
来自数学导入产品
来自sympy导入因子
定义A025441号（n） ：
f=因子（n）.items（）
