%I#23 2021年5月27日17:16:19

%S 441693108911971449161718811953220522772541270927932961，

%电话：3069324933813465371739333969422142574389447345574653，

%电话：476149775229530153135445573358415929598562376321641764896633

%N个数有多个因子分解成S-素数。定义见A054520和A057948。

%与3模4同余的k>=4素因子（具有多重性）的C数，其中k-1不相等_Charlie Neder_，2018年11月3日

%H Eric M.Schmidt，n的表，n=1..10000的a（n）</a>

%e 2205位于S={1,5,9，…4i+1，…}，2205=5*9*49=5*21^2；5、9、21和49是S-素数（A057948）。

%o（鼠尾草）

%o定义A057949_list（绑定）：

%o数字项=（bound-1）//4+1

%o M=[1]*数字项

%o对于范围（1，数字项）中的k：

%o如果M[k]==1：

%o k功率=k

%o当kpower<numterms:

%o阶跃=4*k功率+1

%o对于范围内的j（kpower，numterms，step）：

%o M[j]*=4*k+1

%o kpower=4*kpower*k+kpower+k

%o#现在M[k]包含p^e项的乘积，其中p是S-素数

%o#和e是最大的，使得p^e除以4*k+1

%o如果M[k]>4*k+1]，则返回范围内k的[4*k+1（数字项）

%o#_Eric M.Schmidt，2016年12月11日

%o（PARI）ok（n）={if（n%4==1，my（f=因子（n））；my（s=[f[i，2]|i<-[1..#f~]，f[i、1]%4==3]）；vecsum（s）>=4&&vecmax（s）<vecsum

%Y参考A054520、A057948、A057950。

%Y参考A343826（仅1路）、A343827（正好2路）和A343828（正好3路）。

%K nonn公司

%O 1，1

%A _Jud McCranie，2000年10月14日

%E偏移由_Eric M.Schmidt_校正，2016年12月11日