// xtqqwq
#include<bits/stdc++.h>
#define pb push_back
#define fi first
#define se second
#define mp make_pair
using namespace std;
typedef pair<int, int> pii;
typedef long long ll;
template<typename T>bool chkmax(T&x,T y) { return (y>x) ? x=y,1 : 0;}
template<typename T>bool chkmin(T&x,T y) { return (y<x) ? x=y,1 : 0;}
template<typename T>void inline read(T&x) {
    int f=1;x=0;char s=getchar();
    while (s<'0'||s>'9') { if (s=='-') f=-1;s=getchar();}
    while (s<='9'&&s>='0') x=x*10+(s^48),s=getchar();
    x*=f;
}
const int N=1e5+5,P=998244353;
int n,m,op;
int inline power(int a,int b) {
	int res=1;
	while (b) {
		if (b&1) res=(ll)res*a % P;
		a=(ll)a*a % P;
		b>>=1;
	}
	return res;
}
namespace Task1{
	int f[N];
	int find(int x) { return x==f[x] ? x : f[x]=find(f[x]);}
	void inline merge(int x,int y) {
		f[find(x)]=find(y);
	}
	set<int>s[N];
	void inline work() {
		for (int i=1;i<=n;i++) f[i]=i;
		for (int i=1;i<n;i++) {
			int a,b;scanf("%d%d",&a,&b);
			s[a].insert(b),s[b].insert(a);
		}
		for (int i=1;i<n;i++) {
			int a,b;scanf("%d%d",&a,&b);
			if (s[a].count(b)) merge(a,b);
		}
		int cnt=0;
		for (int i=1;i<=n;i++)
			if (find(i)==i) cnt++;
		printf("%d\n",power(m,cnt));
	}
}
void inline add(int&x,int y) {
	x+=y;
	if (x>=P) x-=P;
}
namespace Task2{
	vector<int>g[N];
	int fa[N],f[N][2],h[2],I,ans;
	void upd(int u,int v) {
		//cut
		for (int i=0;i<2;i++)
			add(h[i],1ll*f[u][i]*f[v][1] % P);
		//link
		for (int i=0;i<2;i++)
			for (int j=0;j+i<2;j++)
				add(h[i+j],1ll*f[u][i]*f[v][j] % P*I % P);
		for (int i=0;i<2;i++)
			f[u][i]=h[i],h[i]=0;
	}
	void dfs(int u) {
		f[u][0]=1,f[u][1]=1;
		for (int v: g[u]) {
			if (v==fa[u]) continue;
			fa[v]=u;
			dfs(v);
			upd(u,v);
		}
	}
	void inline work() {
		I=(1ll-m+P)*power(m,P-2) % P*power(n,P-2) % P;
		for (int i=1;i<n;i++) {
			int u,v;read(u),read(v);
			g[u].pb(v),g[v].pb(u);
		}
		int s=1ll*power(m,n)*power(n,n-2) % P;
		dfs(1);
		int ans=f[1][1];
		ans=1ll*ans*s % P;
		printf("%d\n",ans);
	}
}
namespace Task3{
	typedef vector<int>Poly;
	#define pb push_back
	const int S=8e5+5,P=998244353,G=3;
	int A[S],rev[S],mod,inv[S],fact[S],infact[S];
	int lim=1,len=0,W[20][S];
	int inline power(int a,int b,int Mod=P) {
		int res=1;
		while (b) {
			if (b&1) res=(ll)res*a % Mod;
			a=(ll)a*a % Mod;
			b>>=1;
		}
		return res;
	}
	int Gi=power(G,P-2,P),inv2=power(2,P-2,P);
	void inline NTT(int c[],int lim,int o) {
		for (int i=0;i<lim;i++)
			if (i<rev[i]) swap(c[i],c[rev[i]]);
		for (int k=1,t=0;k<lim;k<<=1,t++) {
			for (int i=0;i<lim;i+=(k<<1)) {
				for (int j=0;j<k;j++) {
					int u=c[i+j],v=(ll)c[i+k+j]*W[t][j] % P;
					c[i+j]=u+v>=P ? u+v-P : u+v;
					c[i+j+k]=u-v<0 ? u-v+P : u-v;
				}
			}
		}
		if (o==-1) {
			reverse(c+1,c+lim);
			int inv=power(lim,P-2,P);
			for (int i=0;i<lim;i++)
				c[i]=(ll)c[i]*inv % P;
		}
	}
	void inline setN(int n) {
		lim=1,len=0;
		while (lim<n) lim<<=1,len++;
		for (int i=0;i<lim;i++)
			rev[i]=(rev[i>>1]>>1)|((i&1)<<(len-1));
	}
	Poly inline NTT(Poly a,int o) {
		int n=a.size();
		for (int i=0;i<n;i++) A[i]=a[i];
		NTT(A,lim,o);
		a.clear();
		for (int i=0;i<lim;i++) a.push_back(A[i]),A[i]=0;
		return a;
	}
	Poly operator+(const Poly a,const Poly b)  {
		Poly c(max(a.size(),b.size()));
		for (int i=0;i<c.size();i++) {
			if (i<a.size()) {
				c[i]+=a[i];if (c[i]>=P) c[i]-=P;
			}
			if (i<b.size()) {
				c[i]+=b[i];if (c[i]>=P) c[i]-=P;
			}
		}
		return c;
	}
	Poly operator-(const Poly a,const Poly b)  {
		Poly c(max(a.size(),b.size()));
		for (int i=0;i<c.size();i++) {
			if (i<a.size()) {
				c[i]+=a[i];if (c[i]>=P) c[i]-=P;
			}
			if (i<b.size()) {
				c[i]-=b[i];if (c[i]<0) c[i]+=P;
			}
		}
		return c;
	}
	Poly inline mul (Poly a,Poly b,int newn=-1) {
		if (newn==-1) newn=a.size()+b.size()-1;
		setN(a.size()+b.size()-1);
		Poly c=NTT(a,1),d=NTT(b,1);
		for (int i=0;i<lim;i++) c[i]=(ll)c[i]*d[i] % P;
		d=NTT(c,-1);d.resize(newn);
		return d;
	}
	Poly inline reverse(Poly a) {
		int n=a.size()-1;
		for (int i=0;i<n-i;i++) swap(a[i],a[n-i]);
		return a;
	}
	Poly inline dx(Poly a) {
		int n=a.size();
		Poly b;b.resize(n-1);
		for (int i=0;i<b.size();i++) b[i]=a[i+1]*(i+1ll) % P;
		return b;
	}
	Poly inline F(Poly a) {
		int n=a.size();
		Poly b;b.resize(n+1);
		for (int i=1;i<b.size();i++) b[i]=(ll)a[i-1]*inv[i] % P;
		return b;
	}
	 
	Poly polyInv(Poly a) {
		int n=a.size();
		if (n==1) { Poly b;b.push_back(power(a[0],P-2,P));return b;}
		Poly b=a;b.resize((n+1)>>1);
		b=polyInv(b);
		setN(2*n);
		a=NTT(a,1),b=NTT(b,1);
		for (int i=0;i<lim;i++)
			b[i]=(ll)b[i]*(2ll-(ll)a[i]*b[i] % P+P) % P;
		b=NTT(b,-1);
		b.resize(n);
		return b;
	}
	//注意必须保证 n>=m
	void inline div (Poly a,Poly b,Poly&Q,Poly&R) {
		int n=a.size()-1,m=b.size()-1;
		Poly ar=reverse(a),br=reverse(b);
		ar.resize(n-m+1),br.resize(n-m+1);
		Q=reverse(mul(ar,polyInv(br),n-m+1));
		R=a-mul(b,Q);R.resize(m);
	}
	Poly t(1,1);
	Poly inline ln(Poly a) {
		Poly b=F(mul(dx(a),polyInv(a)));
		b.resize(a.size());
		return b;
	}
	Poly exp(Poly a) {
		int n=a.size();
		if (n==1) return t;
		Poly b=a;b.resize((n+1)>>1);
		b=exp(b);b.resize(n);
		Poly c=a-ln(b);
		(c[0]+=1) %=P;
		b=mul(b,c,a.size());
		return b;
	}
	void cdq(int l,int r) {
		if (r-l<=0) return;
		int mid=(l+r)>>1,len=r-l;
		cdq(l,mid);
		//Do sth
		cdq(mid+1,r);
	}
	void inline preInv(int n) {
		inv[1]=1;
		for (int i=2;i<=n;i++)
			inv[i]=((ll)P-P/i)*inv[P % i] % P;
	}
	void inline factPrework(int n) {
		fact[0]=infact[0]=1;
		for (int i=1;i<=n;i++) fact[i]=(ll)fact[i-1]*i % P;
		infact[n]=power(fact[n],P-2);
		for (int i=n-1;i;i--) infact[i]=infact[i+1]*(i+1ll) % P;
	}
	//用到的最大的 n
	void inline init(int n) {
		preInv(n);
		factPrework(n);
		setN(2*n);
		for (int k=1,t=0;k<lim;k<<=1,t++) {
			int wn=power(G,(P-1)/(k<<1));
			W[t][0]=1;
			for (int j=1;j<k;j++) W[t][j]=(ll)W[t][j-1]*wn % P;
		}
	}
	int I;
	void inline work() {
		init(1e5+1);
		I=(1ll-m+P)*power(m,P-2) % P*power(n,P-2) % P*power(n,P-2) % P;
		Poly A(n+1,0);
		int v=1;
		for (int i=1;i<=n;i++) {
			A[i]=1ll*infact[i-1]*v % P*(i==1 ? 1 : power(i,i-2)) % P*i % P;
			v=1ll*v*I % P;
		}
		Poly B=exp(A);
		int s=1ll*power(m,n)*power(n,n-2) % P*power(n,n-2) % P;
		int ans=1ll*B[n]*s % P*fact[n] % P;
		printf("%d\n",ans);
	}
}
int main() {
	scanf("%d%d%d",&n,&m,&op);
	if (op==0) Task1::work();
	else if (op==1) Task2::work();
	else if (op==2) Task3::work();
	return 0;
}