题目链接

题目思路

这个\(dp\)比较巧妙

设\(dp[i][j]\)表示在第\(i\)次操作后，还有\(j\)个数不一样的答案

那么最后输出\(dp[k][0]\)

转移方程即每次枚举选了s个与最终态不同的灯，m-s个与最终态相同的灯操作

那么\(dp\)方程即为

$ dp[i+1][j-s+m-s]=dp[i+1][j-s+m-s]+dp[i][j]c(j,s)c(n-j,m-s);$

代码

#include<bits/stdc++.h>
#define fi first
#define se second
#define debug cout<<"I AM HERE"<<endl;
using namespace std;
typedef long long ll;
const int maxn=1e2+5,inf=0x3f3f3f3f,mod=998244353;
const double eps=1e-6;
int n,k,m;
char s[maxn],t[maxn];
ll dp[maxn][maxn];
ll fac[maxn],finv[maxn];
ll qpow(ll a,ll b){
    ll ans=1,base=a;
    while(b){
        if(b&1) ans=ans*base%mod;
        base=base*base%mod;
        b=b>>1;
    }
    return ans;
}
void init(int n){
    fac[0]=finv[0]=1;
    for(int i=1;i<=n;i++){
        fac[i]=fac[i-1]*i%mod;
    }
    finv[n]=qpow(fac[n],mod-2);
    for(int i=n-1;i>=1;i--){
        finv[i]=finv[i+1]*(i+1)%mod;
    }
}
ll c(ll a,ll b){
    if(a<b) return 0;
    ll ans=fac[a]*finv[b]%mod*finv[a-b]%mod;
    return ans;
}
signed main(){
    init(100);
    int _;scanf("%d",&_);
    while(_--){
        scanf("%d%d%d",&n,&k,&m);
        scanf("%s %s",s+1,t+1);
        memset(dp,0,sizeof(dp));
        int cnt=0;
        for(int i=1;i<=n;i++){
            if(s[i]!=t[i]){
                cnt++;
            }
        }
        dp[0][cnt]=1;
        for(int i=0;i<=k-1;i++){
            for(int j=0;j<=n;j++){
                for(int s=0;s<=min(j,m);s++){
                    if(j-s+m-s<0) continue;
                    if(j-s+m-s>n) continue;
                    dp[i+1][j-s+m-s]=(dp[i+1][j-s+m-s]+dp[i][j]*c(j,s)%mod*c(n-j,m-s)%mod)%mod;
                }
            }
        }
        printf("%lld\n",dp[k][0]);
    }
    return 0;
}