|
（SageMath）
P.<P>=QQ[]
定义多流（nn，pol=P（0），kk=0）：
如果kk>=nn：
产量pol-1
其他：
对于sx范围内的ii（二项式（nn，kk）+1）：
对于多流中的xx（nn，pol+2*ii*p^kk*（1-p）^（nn-kk），kk+1）：
产量xx
定义计算（nn）：
解决方案=设置（）
对于多流中的pol（nn）：
rootlist=[xx[0]表示pol.roots中的xx（环=QQbar）]
对于rootlist中的root：
如果root.real（）==root且0<=root<=1：
解决方案+=集合（[root]）
返回长度（解决方案）
| 