結果
問題 | No.2527 H and W |
ユーザー | anonymously |
提出日時 | 2023-11-03 23:13:55 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 650 ms / 2,000 ms |
コード長 | 912 bytes |
コンパイル時間 | 1,594 ms |
コンパイル使用メモリ | 166,708 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-09-25 21:26:20 |
合計ジャッジ時間 | 5,240 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 23 |
ソースコード
#include <bits/stdc++.h> using namespace std; const long MODULO=998244353L; long modinv(long a) { long b=MODULO, u=1L, v=0L; while(b){ long t=a/b; a-=(t*b); swap(a, b); u-=(t*v); swap(u, v); } u%=MODULO; if(u<0)u+=MODULO; return u; } long modplus(long a,long b){ long res=(a+b)%MODULO; if(res<0)res+=MODULO; return res; } long modtimes(long a,long b){ return (a*b)%MODULO; } long comb(long n,long r){ if(r*2>n){ return comb(n,n-r); }else{ long nCr=1L; for(long i=1;i<=r;i++){ nCr=modtimes(modtimes(nCr,n-r+i),modinv(i)); } return nCr; } } int main(){ long h,w,k; cin >> h >> w >> k; long ans=0L; for(long i=max(1L,(k+w-1)/w);i<=min(h,k);i++){ if(k%i==0){ long j=k/i; //i行、j列を残して選択し、赤くする ans=modplus(ans,modtimes(comb(h,i),comb(w,j))); } } cout << ans << endl; return 0; }