1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57
   | #include<bits/stdc++.h> using namespace std;   void read(int &x) {     x=0;int f=1;char ch=getchar();     for(;!isdigit(ch);ch=getchar()) if(ch=='-') f=-f;     for(;isdigit(ch);ch=getchar()) x=x*10+ch-'0';x*=f; }   void print(int x) {     if(x<0) putchar('-'),x=-x;     if(!x) return ;print(x/10),putchar(x%10+48); } void write(int x) {if(!x) putchar('0');else print(x);putchar('\n');}   #define lf double #define ll long long    #define pii pair<int,int > #define vec vector<int >   #define pb push_back #define mp make_pair #define fr first #define sc second   #define data asd09123jdf02i3h   #define FOR(i,l,r) for(int i=l,i##_r=r;i<=i##_r;i++)   const int maxn = 1e6+10; const int inf = 1e9; const lf eps = 1e-8; const int mod = 998244353;   int n,m,k,s[5002][5002];   int qpow(int a,int x) {     int res=1;     for(;x;x>>=1,a=1ll*a*a%mod) if(x&1) res=1ll*res*a%mod;     return res; }   int main() {     read(n),read(m),read(k);     s[0][0]=1;     for(int i=1;i<=k;i++)          for(int j=1;j<=i;j++)             s[i][j]=(1ll*s[i-1][j]*j%mod+s[i-1][j-1])%mod;     int ans=0,a=1,b=1,p=qpow(m,mod-2);     for(int i=1;i<=k;i++) {         a=1ll*a*n%mod,n--;         b=1ll*b*p%mod;         ans=(ans+1ll*s[k][i]*a%mod*b%mod)%mod;     }write(ans);     return 0; }
   |