結果

問題 No.1043 直列大学
ユーザー Yukino DX.
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff # 

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
        }
    }
}
0