結果
問題 | No.1043 直列大学 |
ユーザー |
|
提出日時 | 2024-09-10 11:33:10 |
言語 | Rust (1.83.0 + proconio) |
結果 |
AC
|
実行時間 | 198 ms / 2,000 ms |
コード長 | 4,214 bytes |
コンパイル時間 | 14,162 ms |
コンパイル使用メモリ | 383,576 KB |
実行使用メモリ | 77,696 KB |
最終ジャッジ日時 | 2024-09-10 11:33:27 |
合計ジャッジ時間 | 16,291 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 28 |
ソースコード
use proconio::input;fn main() {input! {n:usize,m:usize,v:[usize;n],r:[usize;m],a:usize,b:usize,}let sum_v = v.iter().sum::<usize>();let sum_r = r.iter().sum::<usize>();let dp_v = calc_dp(&v, n);let dp_r = calc_dp(&r, m);let cum_dp_v = dp_v.iter().scan(ModInt::new(0), |cum, dpi_v| {*cum += *dpi_v;Some(*cum)}).collect::<Vec<_>>();let mut ans = ModInt::<MOD>::new(0);for i in 1..=sum_r {ans += dp_r[i] * (cum_dp_v[(b * i).min(sum_v)] - cum_dp_v[(a * i - 1).min(sum_v)]);}println!("{}", ans);}fn calc_dp(a: &Vec<usize>, n: usize) -> Vec<ModInt<MOD>> {let sum_a = a.iter().sum::<usize>();let mut dp = vec![vec![ModInt::<MOD>::new(0); sum_a + 1]; n + 1];dp[0][0] = ModInt::new(1);for i in 0..n {for j in 0..=sum_a {dp[i + 1][j] = dp[i + 1][j] + dp[i][j];if j + a[i] <= sum_a {dp[i + 1][j + a[i]] = dp[i + 1][j + a[i]] + dp[i][j];}}}dp[n].clone()}const MOD: usize = 1000000007;use modint::*;mod modint {use std::fmt;use std::ops;#[derive(Copy, Clone, PartialEq, Eq)]pub struct ModInt<const MOD: usize> {pub val: usize,}impl<const MOD: usize> ModInt<MOD> {pub fn new(val: usize) -> Self {Self { val: val % MOD }}pub fn pow(mut self, mut e: usize) -> Self {let mut res = Self::new(1);while 0 < e {if e & 1 != 0 {res *= self;}self *= self;e >>= 1;}res}}impl<const MOD: usize> From<usize> for ModInt<MOD> {fn from(value: usize) -> Self {Self { val: value % MOD }}}impl<const MOD: usize> fmt::Display for ModInt<MOD> {fn fmt(&self, f: &mut fmt::Formatter<'_>) -> fmt::Result {write!(f, "{}", self.val)}}impl<const MOD: usize> fmt::Debug for ModInt<MOD> {fn fmt(&self, f: &mut fmt::Formatter<'_>) -> fmt::Result {write!(f, "{}", self.val)}}impl<const MOD: usize> ops::Neg for ModInt<MOD> {type Output = Self;fn neg(self) -> Self::Output {Self {val: (MOD - self.val) % MOD,}}}impl<const MOD: usize> ops::Add for ModInt<MOD> {type Output = Self;fn add(self, rhs: Self) -> Self::Output {Self {val: (self.val + rhs.val) % MOD,}}}impl<const MOD: usize> ops::AddAssign for ModInt<MOD> {fn add_assign(&mut self, rhs: Self) {*self = *self + rhs;}}impl<const MOD: usize> ops::Mul for ModInt<MOD> {type Output = Self;fn mul(self, rhs: Self) -> Self::Output {Self {val: self.val * rhs.val % MOD,}}}impl<const MOD: usize> ops::MulAssign for ModInt<MOD> {fn mul_assign(&mut self, rhs: Self) {*self = *self * rhs;}}impl<const MOD: usize> ops::Sub for ModInt<MOD> {type Output = Self;fn sub(mut self, rhs: Self) -> Self::Output {if self.val < rhs.val {self.val += MOD;}Self {val: (self.val - rhs.val) % MOD,}}}impl<const MOD: usize> ops::SubAssign for ModInt<MOD> {fn sub_assign(&mut self, rhs: Self) {if self.val < rhs.val {self.val += MOD;}*self = *self - rhs;}}impl<const MOD: usize> ops::Div for ModInt<MOD> {type Output = Self;fn div(self, rhs: Self) -> Self {assert!(rhs.val != 0);self * rhs.pow(MOD - 2)}}impl<const MOD: usize> ops::DivAssign for ModInt<MOD> {fn div_assign(&mut self, rhs: Self) {*self = *self / rhs}}}