(PARI)isok(n)=总和(k=0,n,!是_A005153号(二项式(n,k))==n-1;
(Python)
从itertools导入计数,islice
来自数学导入梳
来自sympy导入因子
定义A334084型_gen(startvalue=1):#术语生成器>=startvalue
对于计数中的n(max(startvalue,1)):
对于范围(1,n)中的k:
c=梳(n,k)
l=(~c&c-1).bit_length()
如果l>0:
P=(1<<l+1)-1
对于因子(c>>l).items()中的p,e:
如果p>1+p:
打破
P*=(P**(e+1)-1)//(P-1)
其他:
打破
其他:
产量n