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::(); let sum_r = r.iter().sum::(); let dp_v = calc_dp(&v, n); let dp_r = calc_dp(&r, m); let cum_dp_v = dp_v.cum(|cum, dpi| cum + dpi, ModInt::new(0)); let mut ans = ModInt::::new(0); for i in 1..=sum_r { ans += dp_r[i] * (cum_dp_v[(b * i + 1).min(sum_v + 1)] - cum_dp_v[(a * i).min(sum_v)]); } println!("{}", ans); } fn calc_dp(a: &Vec, n: usize) -> Vec> { let sum_a = a.iter().sum::(); let mut dp = vec![vec![ModInt::::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() } use cumulative::*; mod cumulative { pub trait Cumulative where T: Copy, F: Fn(T, T) -> T, { fn cum(&self, op: F, e: T) -> Vec; fn cum_rev(&self, op: F, e: T) -> Vec; } impl Cumulative for [T] where T: Copy, F: Fn(T, T) -> T, { //cum[i]:=op[0,i) fn cum(&self, op: F, e: T) -> Vec { let n = self.len(); let mut res = vec![e; n + 1]; for i in 0..n { res[i + 1] = op(res[i], self[i]); } res } //cum_rev[i]:=op[i,n) fn cum_rev(&self, op: F, e: T) -> Vec { let n = self.len(); let mut res = vec![e; n + 1]; for i in (0..n).rev() { res[i] = op(res[i + 1], self[i]); } res } } } const MOD: usize = 1000000007; use modint::*; mod modint { use std::fmt; use std::ops; #[derive(Copy, Clone, PartialEq, Eq)] pub struct ModInt { pub val: usize, } impl ModInt { 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 From for ModInt { fn from(value: usize) -> Self { Self { val: value % MOD } } } impl fmt::Display for ModInt { fn fmt(&self, f: &mut fmt::Formatter<'_>) -> fmt::Result { write!(f, "{}", self.val) } } impl fmt::Debug for ModInt { fn fmt(&self, f: &mut fmt::Formatter<'_>) -> fmt::Result { write!(f, "{}", self.val) } } impl ops::Neg for ModInt { type Output = Self; fn neg(self) -> Self::Output { Self { val: (MOD - self.val) % MOD, } } } impl ops::Add for ModInt { type Output = Self; fn add(self, rhs: Self) -> Self::Output { Self { val: (self.val + rhs.val) % MOD, } } } impl ops::AddAssign for ModInt { fn add_assign(&mut self, rhs: Self) { *self = *self + rhs; } } impl ops::Mul for ModInt { type Output = Self; fn mul(self, rhs: Self) -> Self::Output { Self { val: self.val * rhs.val % MOD, } } } impl ops::MulAssign for ModInt { fn mul_assign(&mut self, rhs: Self) { *self = *self * rhs; } } impl ops::Sub for ModInt { 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 ops::SubAssign for ModInt { fn sub_assign(&mut self, rhs: Self) { if self.val < rhs.val { self.val += MOD; } *self = *self - rhs; } } impl ops::Div for ModInt { type Output = Self; fn div(self, rhs: Self) -> Self { assert!(rhs.val != 0); self * rhs.pow(MOD - 2) } } impl ops::DivAssign for ModInt { fn div_assign(&mut self, rhs: Self) { *self = *self / rhs } } }