[hdu6391]Lord Li's problem

首先发现结果与需要改变的具体位置无关，只和需要改变的位置的个数有关，因此设f[i][j]表示选取了i个数字异或结果有j个1，只要分析接下来选择的数和这j个1有几个重合即可：

1. 三个数字全部重合，即$f[i][j+3]+=f[i-1][j]\cdot c(n-j,3)$

2. 有两个数字重合，即$[i][j+1]+=f[i-1][j]\cdot c(n-j,2)\cdot j$

3. 有一个数字重合，即$f[i][j-1]+=f[i-1][j]\cdot (n-j)\cdot c(j,2)$

4. 有0个数字重合，即$f[i][j-3]+=f[i-1][j]\cdot c(j,3)$

然而这并不正确，因为不能重复选择，因此去掉($c(n,3)-(i-2))\cdot f[i-2][j]$，同时这里我们考虑了操作的顺序和操作的位置，再除以$n!\cdot c(n,m)$才是最终答案。

#include<bits/stdc++.h>
using namespace std;
#define mod 19260817
int n,m,t,sum,inv[101],fac[101],f[101][101];
char s1[101],s2[101];
int c(int n,int m){
    return 1LL*fac[n]*inv[n-m]%mod*inv[m]%mod;
}
int cc(int n,int m){
    return 1LL*inv[n]*fac[n-m]%mod*fac[m]%mod;
}
int main(){
    inv[0]=inv[1]=fac[0]=fac[1]=1;
    for(int i=2;i<=40;i++)fac[i]=1LL*fac[i-1]*i%mod;
    for(int i=2;i<=40;i++)inv[i]=1LL*(mod-mod/i)*inv[mod%i]%mod;
    for(int i=2;i<=40;i++)inv[i]=1LL*inv[i-1]*inv[i]%mod;
    while (scanf("%d%d",&n,&m)!=EOF){
        if ((!n)&&(!m))return 0;
        scanf("%s%s",s1,s2);
        sum=0;
        for(int i=0;i<n;i++)
            if (s1[i]!=s2[i])sum++;
        memset(f,0,sizeof(f));
        f[0][0]=1;
        for(int i=0;i<=m;i++)
            for(int j=0;j<=n;j++){
                if (i>1)f[i][j]=(f[i][j]-f[i-2][j]*(i-1LL)%mod*(c(n,3)-i+2)%mod+mod)%mod;
                for(int k=-1;k<3;k++)
                    if ((n>j+k)&&(j+k>1))
                        f[i+1][j+2*k-1]=(f[i+1][j+2*k-1]+1LL*c(n-j,k+1)*c(j,2-k)*f[i][j])%mod;
            }
        printf("Case #%d: %d\n",++t,1LL*f[m][sum]*cc(n,sum)%mod*inv[m]%mod);
    }
}