結果
問題 | No.1693 Invasion |
ユーザー |
|
提出日時 | 2022-10-18 00:31:44 |
言語 | Rust (1.83.0 + proconio) |
結果 |
AC
|
実行時間 | 118 ms / 2,000 ms |
コード長 | 1,744 bytes |
コンパイル時間 | 14,072 ms |
コンパイル使用メモリ | 376,940 KB |
実行使用メモリ | 83,456 KB |
最終ジャッジ日時 | 2024-06-28 13:04:20 |
合計ジャッジ時間 | 14,853 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 21 |
ソースコード
const INF: usize = 1usize << 60;const MOD: usize = 998244353;fn fact_inv(size: usize) -> (Vec<usize>, Vec<usize>) {let mut facts = vec![0usize; size+1];let mut fact_invs = vec![0usize; size+1];let mut invs = vec![0usize; size+1];facts[0] = 1;facts[1] = 1;fact_invs[0] = 1;fact_invs[1] = 1;invs[1] = 1;for i in 2..=size {facts[i] = facts[i-1] * i % MOD;invs[i] = MOD - invs[MOD % i] * (MOD / i) % MOD;fact_invs[i] = fact_invs[i - 1] * invs[i] % MOD;}(facts, fact_invs)}fn comb(n: usize, r: usize, facts: &Vec<usize>, fact_invs: &Vec<usize>) -> usize {facts[n] * fact_invs[r] % MOD * fact_invs[n-r] % MOD}fn main() {let mut nm = String::new();std::io::stdin().read_line(&mut nm).ok();let nm: Vec<usize> = nm.trim().split_whitespace().map(|s| s.parse().unwrap()).collect();let n = nm[0];let m = nm[1];let mut a = String::new();std::io::stdin().read_line(&mut a).ok();let a: Vec<usize> = a.trim().split_whitespace().map(|s| s.parse().unwrap()).collect();let mut dp = vec![vec![INF; m+1]; n+1];dp[0][0] = 0;for i in 0..n {for j in 0..=m {dp[i+1][j] = dp[i+1][j].min(dp[i][j]);if j+a[i] <= m {dp[i+1][j+a[i]] = dp[i+1][j+a[i]].min(dp[i][j] + 1);}if j >= a[i] {dp[i+1][j] = dp[i+1][j].min(dp[i+1][j-a[i]] + 1);}}}let (facts, fact_invs) = fact_inv(m+10);let mut result = 1usize;for i in 1..=m {if dp[n][i] == INF { continue; }result += comb(m-dp[n][i], i-dp[n][i], &facts, &fact_invs);result %= MOD;}println!("{}", result);}