#include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using ll = long long; const int inf=1<<30; const ll INF=1LL<<62; typedef pair P; typedef pair PP; const ll MOD=998244353; template //関数ポインタ struct lazy_segtree { public: lazy_segtree() : lazy_segtree(0) {} explicit lazy_segtree(int n) : lazy_segtree(std::vector(n, e())) {} explicit lazy_segtree(const std::vector& v) : _n(int(v.size())) { log = ceil_pow2(_n); size = 1 << log; d = std::vector(2 * size, e()); lz = std::vector(size, id()); for (int i = 0; i < _n; i++) d[size + i] = v[i]; for (int i = size - 1; i >= 1; i--) { update(i); } } void set(int p, S x) {//a[p]=xの操作を行う aのp番目にxを代入する assert(0 <= p && p < _n); p += size; for (int i = log; i >= 1; i--) push(p >> i); d[p] = x; for (int i = 1; i <= log; i++) update(p >> i); } S get(int p) {//aのp番目の要素を返す a[p] assert(0 <= p && p < _n); p += size; for (int i = log; i >= 1; i--) push(p >> i); return d[p]; } // op( a[l],a[l+1],a[l+2],...,a[r-1] )を返す l==rの時はe()を返す S prod(int l, int r) { assert(0 <= l && l <= r && r <= _n); if (l == r) return e();//範囲が0の時 l += size; r += size; for (int i = log; i >= 1; i--) { if (((l >> i) << i) != l) push(l >> i); if (((r >> i) << i) != r) push((r - 1) >> i); } S sml = e(), smr = e();//関数eには長さ0に対応する情報を書く while (l < r) { if (l & 1) sml = op(sml, d[l++]);// if (r & 1) smr = op(d[--r], smr); l >>= 1; r >>= 1; } return op(sml, smr); } S all_prod() { return d[1]; }//op(a[0],...,a[n-1])を計算する void apply(int p, F f) {//a[p] = f(a[p]) assert(0 <= p && p < _n); p += size; for (int i = log; i >= 1; i--) push(p >> i); d[p] = mapping(f, d[p]); for (int i = 1; i <= log; i++) update(p >> i); } void apply(int l, int r, F f) { //i=l,l+1,l+2,...,r-1についてa[i] = f(a[i]) assert(0 <= l && l <= r && r <= _n); if (l == r) return; l += size; r += size; for (int i = log; i >= 1; i--) { if (((l >> i) << i) != l) push(l >> i); if (((r >> i) << i) != r) push((r - 1) >> i); } { int l2 = l, r2 = r; while (l < r) { if (l & 1) all_apply(l++, f); if (r & 1) all_apply(--r, f); l >>= 1; r >>= 1; } l = l2; r = r2; } for (int i = 1; i <= log; i++) { if (((l >> i) << i) != l) update(l >> i); if (((r >> i) << i) != r) update((r - 1) >> i); } } template int max_right(int l) { return max_right(l, [](S x) { return g(x); }); } template int max_right(int l, G g) { assert(0 <= l && l <= _n); assert(g(e())); if (l == _n) return _n; l += size; for (int i = log; i >= 1; i--) push(l >> i); S sm = e(); do { while (l % 2 == 0) l >>= 1; if (!g(op(sm, d[l]))) { while (l < size) { push(l); l = (2 * l); if (g(op(sm, d[l]))) { sm = op(sm, d[l]); l++; } } return l - size; } sm = op(sm, d[l]); l++; } while ((l & -l) != l); return _n; } template int min_left(int r) { return min_left(r, [](S x) { return g(x); }); } template int min_left(int r, G g) { assert(0 <= r && r <= _n); assert(g(e())); if (r == 0) return 0; r += size; for (int i = log; i >= 1; i--) push((r - 1) >> i); S sm = e(); do { r--; while (r > 1 && (r % 2)) r >>= 1; if (!g(op(d[r], sm))) { while (r < size) { push(r); r = (2 * r + 1); if (g(op(d[r], sm))) { sm = op(d[r], sm); r--; } } return r + 1 - size; } sm = op(d[r], sm); } while ((r & -r) != r); return 0; } private: int _n, size, log; std::vector d; std::vector lz; // 2^x がnを超えるときのxを返す int ceil_pow2(int n) { int x = 0; while ((1U << x) < (unsigned int)(n)) x++; return x; } //opに従ってアップデートする opにはs[l,m)とs[m,r)からs[l,r)を計算する手法を書く void update(int k) { d[k] = op(d[2 * k], d[2 * k + 1]); } void all_apply(int k, F f) { d[k] = mapping(f, d[k]);//更新クエリにより, s[l,r)がどのように変化するかを表すのがmapping if (k < size) lz[k] = composition(f, lz[k]);//クエリを処理する順番 先にlz[k]を行った後で次にfを行う } void push(int k) { all_apply(2 * k, lz[k]); all_apply(2 * k + 1, lz[k]); lz[k] = id();//何もしない更新クエリ id } }; using S=ll; using F=ll; S op(S s_l, S s_r) { //s[l,r)のsum についてのものと s[l,r)のlenについてのもの return min(s_l, s_r); } S e() {//長さ区間が0に対する情報 return INF; } S mapping(F f, S s) { //更新クエリによってs[l,r]がどのように変化するか // sumは変わるけど lenは変わらない return f+s; } F composition(F f, F g) { //クエリを適用する順番はg,f 引数はf,gの順で書く(仕様) return f+g; } F id() { //何もしない 更新クエリ(b,c)=(1,0)の時が何もしないを表す return 0; } ll solve(vector R,vector B,ll K){ ll N=R.size(),M=B.size(); vector> balls; for(int i=0;i dp(N+M+1); int zero=N; dp.set(zero,0); for(auto e:balls){ ll num=e.first; int color=e.second; if(color==1){ //赤の場合 dp.apply(0,zero,num);//[0,zero), dp.apply(zero,N+M+1,-num);//[zero,N+M+1) zero--; }else{ //青が来た場合 ll sv=dp.prod(zero,zero+1);//[zero,zero+1) dp.apply(0,zero+1,-num);//[0,zero+1) dp.apply(zero+1,N+M+1,num);//[zero+1,N+M+1) dp.set(zero+1,min(sv,dp.prod(zero+1,zero+2))); zero++; } } return dp.prod(zero,zero+1)/K; } int main(){ ll N,M,K; cin>>N>>M>>K; vector B(N),R(M); for(int i=0;i>B[i]; } for(int i=0;i>R[i]; } if(N=M map,vector>> mp; for(int i=0;i b=e.second.first; vector r=e.second.second; if(b.size()0){ ans+=solve(r,b,K); } } cout<