抄题解者的求救！！！

#include <bits/stdc++.h>

#define i64 long long
#define rep(i,l,r) for(int i=(l);i<=(r);i++)
#define fdn(i,r,l) for(int i=(r);i>=(l);i--)
#define pii pair<int,int>
using namespace std;

typedef long long ll;
typedef double db;
typedef __int128 i128;

const int N=1001;
const int M=201;
const ll mod=1e9+7;

int n,m,K;
string a,b;
ll dp[2][M][M][2];

int main()
{
#ifndef ONLINE_JUDGE
    freopen("in.in","r",stdin);
    freopen("out.out","w",stdout);
#endif
    // ios::sync_with_stdio(false);cin.tie(nullptr);cout.tie(nullptr);
    cin>>n>>m>>K;
    cin>>a>>b; a=' '+a,b=' '+b;
    dp[1][0][0][0]=dp[0][0][0][0]=1;
    memset(dp,0,sizeof(dp));
    bool u=1;
    for(int i=1;i<=n;i++,u^=1) rep(j,1,m) rep(k,1,K)
    {
        if(a[i]==b[j])
        {
            dp[u][j][k][0]=(dp[u^1][j][k][0]+dp[u^1][j][k][1])%mod;
            dp[u][j][k][1]=(dp[u^1][j-1][k][1]+\
                        (dp[u^1][j-1][k-1][1]+dp[u^1][j-1][k-1][0])%mod)%mod;
        }
        else
        {
            dp[u][j][k][0]=(dp[u^1][j][k][0]+dp[u^1][j][k][1])%mod;
            dp[u][j][k][1]=0;
        }
    }
    cout<<(dp[n&1][m][K][0]+dp[n&1][m][K][1])%mod<<endl;
}