C: =f->子({x=1,y=1},f);
#查找CA中由规则定义的0到M代的ON单元数
#当n-1为奇数时,如果nbd中的ON单元数为
#其中nbd由多项式或Laurent级数f(x,y)定义。
奇数CA:=proc(f,M)全局C;局部n,a,i,f2,p;
f2:=简化(展开(f))mod 2;
a: =[];p: =1;
对于从0到M的n,做a:=[op(a),C(p)];p: =扩展(p*f2)模块2;日期:
l打印([seq(a[i],i=1..nops(a)]);
结束;
f: =1/(x*y)+1/x+1/y+y+x/y+x+x*y;
OddCA(f,130);
|