大=100 000s=0S=设置([])看不见=1seen(v)=if(v<big,bit test(s,v),set search(s,v))参见(v)=如果(v<大,s=位或(s,2^v),s=集合并集(s,集合([v]));while(可见(不可见),不可见++)b(n)={local(x='x);subst(lift(Mod(1,2)*subst)(Pol(binary(n),x),x,1+x)),x(2)});其他(p)={参见(p);我的(bp=b(p));for(v=看不见,oo,如果(!seen(v)&b(p*v)==bp*b(v),回报(v);););}v=1;对于(n=1,oo,打印(n“”v);v=其他(v));退出