#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define endl '\n'
const int N=5e5+5,mod=998244353;
ll n,m;
string s,st;
ll dp[1050][1050];
ll fac[N],invfac[N];
ll quick_pow(ll a,ll b){
	if(a==0)return 0;
	ll ans=1;
	while(b){
		if(b&1)ans=ans*a%mod;
		a=a*a%mod;
		b=b>>1;
	}
	return ans%mod;
}
ll inv(ll x){
	return quick_pow(x,mod-2);
}
void init(){
	fac[0]=invfac[0]=1;
	for(int i=1;i<N;i++){
		fac[i]=fac[i-1]*i%mod;
	}
	invfac[N-1]=inv(fac[N-1]);
	for(int i=N-2;i>=1;i--){
		invfac[i]=invfac[i+1]*(i+1)%mod;
	}
}
ll C(ll a,ll b){
	if(b<a)return 0;
	return fac[b]*invfac[b-a]%mod*invfac[a]%mod;
}
void solve(){
	ll k;
	cin>>n>>m>>k;
	cin>>s;
	cin>>st;
	ll cnt=0;
	for(int i=0;i<n;i++){
		if(s[i]!=st[i])cnt++;
	}
	for(int i=0;i<=m;i++){
		for(int j=0;j<=n;j++){
			dp[i][j]=0;
		}
	}
	dp[0][0]=1;
	for(int i=1;i<=m;i++){
		for(int j=0;j<=n;j++){
			for(int u=0;u<=j;u++){
				if(k-u>n-j)continue;
				if(j+k-2*u<0||j+k-2*u>n)continue;
				dp[i][j+k-2*u]+=dp[i-1][j]*C(u,j)%mod*C(k-u,n-j)%mod;
				dp[i][j+k-2*u]%=mod;
			}
		}
	}
	cout<<dp[m][cnt]*inv(C(cnt,n))%mod<<endl;
	return ;
}
int main(){
	ios::sync_with_stdio(0);
	cin.tie(0);cout.tie(0);
	int _=1;
	init();
	cin>>_;
	while(_--)
	solve();
	return 0;
}