枚举LCP以及下一位变小成什么,统计出剩下的有几个可以在原位置。
然后枚举剩下的至少有几个在原位置,容斥计算答案。
时间复杂度$O(n^3)$。
#includetypedef long long ll;const int N=70,B=10000,MAXL=30;int n,i,j,a[N],b[N],v[N],fun;inline int max(int a,int b){return a>b?a:b;}struct Num{ int a[MAXL],len,fu; Num(){len=1,fu=a[1]=0;} Num operator+(Num b){ Num c; c.len=max(len,b.len)+2; int i; for(i=1;i<=c.len;i++)c.a[i]=0; if(fu==b.fu){ for(i=1;i<=len;i++)c.a[i]=a[i]; for(i=1;i<=b.len;i++)c.a[i]+=b.a[i]; for(i=1;i<=c.len;i++)if(c.a[i]>=B)c.a[i+1]++,c.a[i]-=B; while(c.len>1&&!c.a[c.len])c.len--; c.fu=fu; }else{ bool flag=0; if(len==b.len){ for(i=len;i;i--)if(a[i]!=b.a[i]){ if(a[i]>b.a[i])flag=1; break; } }else{ if(len>b.len)flag=1; } if(flag){ for(i=1;i<=len;i++)c.a[i]=a[i]; for(i=1;i<=b.len;i++)c.a[i]-=b.a[i]; for(i=1;i<=c.len;i++)if(c.a[i]<0)c.a[i+1]--,c.a[i]+=B; while(c.len>1&&!c.a[c.len])c.len--; c.fu=fu; }else{ for(i=1;i<=b.len;i++)c.a[i]=b.a[i]; for(i=1;i<=len;i++)c.a[i]-=a[i]; for(i=1;i<=c.len;i++)if(c.a[i]<0)c.a[i+1]--,c.a[i]+=B; while(c.len>1&&!c.a[c.len])c.len--; c.fu=b.fu; } } return c; } Num operator-(Num b){ b.fu^=1; return *this+b; } Num operator*(Num b){ Num c; c.len=len+b.len+2; c.fu=fu^b.fu; int i,j; for(i=1;i<=c.len;i++)c.a[i]=0; for(i=1;i<=len;i++)for(j=1;j<=b.len;j++){ c.a[i+j-1]+=a[i]*b.a[j]; if(c.a[i+j-1]>=B){ c.a[i+j]+=c.a[i+j-1]/B;c.a[i+j-1]%=B; if(c.a[i+j]>=B)c.a[i+j+1]+=c.a[i+j]/B,c.a[i+j]%=B; } } while(c.len>1&&!c.a[c.len])c.len--; return c; } void write(){ printf("%d",a[len]); for(int i=len-1;i;i--)printf("%04d",a[i]); } void set(int x){ fu=0; len=1; a[1]=x; }}f[N],C[N][N],ans,tmp;void solve(int x,int y){ int i,j,now; for(i=0;i fun)return; now=fun-now; int ret=0; for(i=x+1;i