結果
問題 | No.1693 Invasion |
ユーザー |
|
提出日時 | 2021-10-01 21:31:08 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 17 ms / 2,000 ms |
コード長 | 1,466 bytes |
コンパイル時間 | 598 ms |
コンパイル使用メモリ | 71,296 KB |
実行使用メモリ | 5,632 KB |
最終ジャッジ日時 | 2024-07-19 10:35:44 |
合計ジャッジ時間 | 1,474 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 21 |
ソースコード
#include <iostream> #include <vector> using namespace std; typedef long long ll; const ll MOD = 998244353; const int INF = 10000000; long long modpow(long long x, long long n){ long long res = 1; long long mul = x; while(n){ if(n % 2) res = res * mul % MOD; mul = mul * mul % MOD; n /= 2; } return res; } struct Combination { int maxn; std::vector<long long> fact; std::vector<long long> ifact; Combination(int maxn) : maxn(maxn){ fact.resize(maxn + 1); ifact.resize(maxn + 1); fact[0] = 1; for(int i = 1; i <= maxn; i++) fact[i] = fact[i - 1] * i % MOD; ifact[maxn] = modpow(fact[maxn], MOD - 2); for(int i = maxn - 1; i >= 0; i--) ifact[i] = ifact[i + 1] * (i + 1) % MOD; } long long com(int n, int k){ if(k > n || k < 0) return 0; if(k == n || k == 0) return 1; return (fact[n] * ifact[k] % MOD) * ifact[n - k] % MOD; } }; int main() { int n, m; cin >> n >> m; int a[102]; for(int i = 0; i < n; i++) cin >> a[i]; int dp[100005]; dp[0] = 0; for(int i = 1; i <= m; i++) dp[i] = INF; for(int i = 0; i < n; i++){ int d[100005]; for(int j = 0; j <= m; j++) d[j] = dp[j]; for(int j = a[i]; j <= m; j++) d[j] = min(d[j], d[j - a[i]] + 1); swap(dp, d); } ll ans = 0; Combination com(m); for(int i = 0; i <= m; i++){ ans = (ans + com.com(m - dp[i], i - dp[i])) % MOD; } cout << ans << endl; }