(鼠尾草)
@缓存函数
def欧拉数(n,k,m):#欧拉数
如果n==0:如果k==0,则返回1,否则为0
返回((m*(n-k)+m-1)*欧拉数(n-1,k-1,m)+
(m*k+1)*欧拉数(n-1,k,m)
@缓存函数
def BS(n,m):#广义标度伯努利数
return(加(加(欧拉数(n,j,m)*二项式(j,n-k)
对于j in(0..n))/((-m)^k*(k+1))对于k in(0..n))
C=mul(滤波器(λp:(n+1)%p==0)或(n%(p-1)==0,素数(n+2))
返回C*BS(n,2)