結果

問題 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
}
}
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0