%I#31 2022年6月28日15:28:47

%S 0,0,0,1,0,1,1,0,2,0,1,1,1,1,0,02,2,0,2,1,0,1,1,1,1,1,0,1，

%温度0,1,2,0,1,1,0,0,1,1,1,2,1,1,1,2,0,1,0,1,1,1,1,1,0,1,11,0,1,3,0,11,1，

%U 0,1,1,1,2,0,0,1,1,1,0,4,0,1,0,11,0,1,11,0,2,1,1,1,1,1,1,1,0,2,3,0,10,1,1,1,0,1,3,0,1,1,1,0,1,0,10,0,2,1,1,1

%N的素因式分解中所有4k+3素数的指数的位或。

%C A001481（2个平方和的数字）给出了该序列中偶数项的位置，而其补码A022544（不是2个平方总和的数字）则给出了奇项的位置。

%如果我们用普通的方法将4k+3素数的指数相加，而不是按位排序（A003986），我们将得到序列A065339。关于这两个序列不同的位置，请参见A260730。

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

%F如果n<3，a（n）=0；此后，对于任何偶数n:a（n）=a（n/2），对于任何n，其最小素因子（A020639）的形式为4k+1:a（n=a（A032742（n）），否则[当A020639（n）的形式是4k+3]a（n。

%F其他身份。对于所有n>=0：

%F A229062（n）=1-A000035（a（n））。[减少模2并进行补码，序列给出A001481的特征函数。]

%F a（n）=a（A097706（n））。[结果仅取决于形式4k+3的素因子。]

%F a（n）=A267116（A097706（n））。

%F a（n）=A267113（A267099（n））。

%e对于n=21=3^1*7^1，我们计算A003986（1,1）=1，因此a（21）=1。

%e对于n=63=3^2*7^1，我们计算A003986（2,1）=A003986=1,2）=3，因此a（63）=3。

%t表[BitOr@@（映射[Last，FactorInteger@n/.{p_，_}/；MemberQ[{0，1，2}，Mod[p，4]]->Nothing]），{n，0，120}]（*_Michael De Vlieger_，2016年2月7日*）

%o（方案）（定义（A260728 n）（cond（（<n 3）0）（（偶数？n）（A260428（/n 2）））（（=1（模（A020639 n）4））（A267028（A032742 n））））；；A003986bi实现按位或（请参见A003986）。

%o（Python）

%o从functools导入reduce

%o从操作员导入或_

%o来自症状输入因子int

%o def A260728（n）：返回reduce（or _，（e for p，e in factorint（n）.items（）if p&3==3），0）#_Chai Wah Wu_，2022年6月28日

%Y参见A000035、A001481、A022544、A003986、A020639、A028234、A032742、A067029、A097706、A229062、A260730。

%Y另请参阅A267113、A267116和A267099。

%Y在n=21时第一次与A065339不同，其中a（21）=1，而A065339-（21）=2。

%K nonn公司

%0、10

%A _Antti Karttunen，2015年8月12日