#include <bits/stdc++.h>
using namespace std;
#define int long long
const int N=55,M=(1<<10)+5,mods=998244353;
int pows(int a,int b){
	if(b==0)return 1;
	int res=pows(a,b>>1);
	res=res*res%mods;
	if(b&1)res=res*a%mods;
	return res;
}
struct node{
	int x,y,w;
}g[N];
int n,m,qs,sz,f[N][N][M],pw[N],cnt[M],jc[M],ny[M],bh[N],idx,dp[N][M][M],C[N][N];
vector<int>zj[M],dy[N];
signed main(){
	vector<int>jl;
	auto fd=[&](int x){
		return lower_bound(jl.begin(),jl.end(),x)-jl.begin();
	};
	cin>>n>>m>>qs;
	for(int i=1;i<=qs;i++){
		cin>>g[i].x>>g[i].y>>g[i].w;
		if(!bh[g[i].y])bh[g[i].y]=++idx;
		jl.push_back(g[i].x);
	}
	for(int i=0;i<(1<<idx);i++){
		for(int j=i;j;j=j-1&i)zj[i].push_back(j);
		zj[i].push_back(0);
		cnt[i]=__builtin_popcount(i);
	}
	sort(jl.begin(),jl.end());
	jl.resize(unique(jl.begin(),jl.end())-jl.begin());
	for(int i=1;i<=qs;i++)dy[fd(g[i].x)].push_back(i);
	sz=jl.size();
	for(int i=0;i<=n;i++)pw[i]=pows(i,m-sz);
	for(int i=0;i<=n;i++)C[i][0]=1;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			C[i][j]=(C[i-1][j]+C[i-1][j-1])%mods;
		}
	}
	f[0][n-idx][(1<<idx)-1]=1;
	for(int i=0;i<n;i++){
		for(int j=0;j<=n;j++){
			for(int k=0;k<(1<<idx);k++){
				for(auto s:zj[k^(1<<idx)-1]){
    				int tmp=pw[j-cnt[k]-cnt[s]];
    				for(int u=0;u<sz;u++){
    					int sm=j-cnt[k]-cnt[s],op=0;
    					for(auto d:dy[u]){
    						if((s>>bh[g[d].y]-1&1^1)&&(k>>bh[g[d].y]-1&1^1)&&g[d].w>i+1)sm--;
    						if(g[d].w==i+1)op=bh[g[d].y];
    					}
    					if(op){
    						if(k>>op-1&1)sm=0;
    						else if(s>>op-1&1)sm=0;
    						else sm=1;
    					}
    					tmp=tmp*sm%mods;
    				}
    				if(cnt[s]&1)dp[j][k][s]=-tmp;
    				else dp[j][k][s]=tmp;
				}
				for(int t=0;t<idx;t++){
					for(auto s:zj[k^(1<<idx)-1]){
						if(s>>t&1)dp[j][k][s]+=dp[j][k][s^1<<t];
					}
				}
				for(auto s:zj[k^(1<<idx)-1])dp[j][k][s]%=mods;
			}
		}
		for(int j=0;j<=n-idx;j++){
			for(int k=0;k<(1<<idx);k++){
				if(!f[i][j][k])continue;
			    for(int l=0;l<=j;l++){
			    	for(auto t:zj[k]){
			    		int ans=0,sl=0;
			    		for(int r=0;r<=l&&n-i-j+l-r>=0;r++){
			    			int sm=dp[n-i-j+l-r][k^t][t];
			    			if(r&1)ans-=sm*C[l][r]%mods;
			    			else ans+=sm*C[l][r]%mods;
			    		}
						ans%=mods;
			    		f[i+1][j-l][k^t]+=f[i][j][k]*ans%mods*pows(i+1,l+cnt[t])%mods*C[j][l]%mods;
					}
			    }
			}
		}
		for(int j=0;j<=n-idx;j++){
			for(int k=0;k<(1<<idx);k++){
				f[i+1][j][k]%=mods;
			//	cout<<i+1<<" "<<j<<" "<<k<<" "<<f[i+1][j][k]<<endl;
			}
		}
		cerr<<i<<endl;
	}
	cout<<(f[n][0][0]+mods)%mods;
}