#include #include using namespace std; using namespace atcoder; typedef long long int ll; typedef long double ld; typedef vector vi; typedef vector vl; typedef vector vvl; typedef vector vvvl; typedef vector vvvvl; typedef vector vb; typedef vector vvb; typedef vector vvvb; typedef vector vvvvb; typedef pair pl; typedef pair ppl; typedef pair pppl; typedef pair pppppl; #define rep(i,a,b) for(int i=(a);i<(b);i++) #define rrep(i,a,b) for(int i=(b)-1;i>=(a);i--) #define all(a) begin(a),end(a) #define sz(a) (int)(a).size() #define F first #define S second #define bs(A,x) binary_search(all(A),x) #define lb(A,x) (ll)(lower_bound(all(A),x)-A.begin()) #define ub(A,x) (ll)(upper_bound(all(A),x)-A.begin()) #define cou(A,x) (ll)(upper_bound(all(A),x)-lower_bound(all(A),x)) templateusing min_priority_queue=priority_queue,greater>; templatebool chmax(T&a,T b){if(abool chmin(T&a,T b){if(b vm; typedef vector vvm; typedef vector vvvm; typedef vector vvvvm; ostream&operator<<(ostream&os,mint a){os<>(istream&is,mint&a){int x;is>>x;a=mint(x);return is;} //*/ templateostream&operator<<(ostream&os,pairp){os<istream&operator>>(istream&is,pair&p){is>>p.F>>p.S;return is;} templateostream&operator<<(ostream&os,vectorv){rep(i,0,sz(v))os<istream&operator>>(istream&is,vector&v){for(T&in:v)is>>in;return is;} int main(){ cin.tie(0)->sync_with_stdio(0); cin.exceptions(cin.failbit); ll N,M,E;cin>>N>>M>>E; vm ex(2e6,1),re(2e6); rep(i,1,2e6)ex[i]=i*ex[i-1]; rep(i,0,2e6)re[i]=1/ex[i]; rep(k,1,N+1){ if(M<=k){cout<<1<1){ if(k==N){ ll m=M-N*(l-1); if(1<=m&&m=m+N*(l-1))p++; } else{ mint q=(mint(N-k)/(N-k+1)).pow(M-k*(l-1))/pow_mod(N-k+1,l-1,mod)*re[l-2],r=0; rep(m,1,k)if(k*(l-1)+m<=M)r+=(mint(N-k+1)/(N-k)).pow(m)*ex[M-m-k*(l-1)+l-2]*re[M-m-k*(l-1)]; p+=q*r; } } ans+=p*pow_mod(l,E,mod); } cout<