BJ模拟（1） D1T1 Delight for a Cat

  Delight for a Cat

题目背景：

thoj22

分析：

首先，如果你没有学习过NOI2008的志愿者招募，请先去A掉那道题这里

首先，假设这只猫在所有天数都是睡觉，我们稍后再决策将哪些天数改为吃东西。

我们用 X_i = 0 表示这只猫在第i天睡觉， X_i = 1 表示吃东西。

则有下列不等式：

m_e ≤ X₁ + X₂+ X₃ +...+ X_k ≤ k − m_s

m_e ≤ X₂ + X₃+ X₄ +...+ X_k+1 ≤ k − m_s

...

m_e ≤ X_n-k+1 + X_n-k+2+...+ X_n ≤ k − m_s

我们令 0 ≤ Y₁ ≤ k – m_s − m_e ，则上述不等式可以转化为以下

的等式：

X₁ + X₂+ X₃ +...+ X_k + Y₁ = k − m_s

X₂ + X₃+ X₄ +...+ X_k+1 + Y₂ = k − m_s

...

X_n−k+1 + X_n−k+2 +...+ X_n +Y_n−k+1 = k − m_s

将这些等式两两相减，得：

X₁ − X_k+1 + Y₁− Y₂=0

X₂ − X_k+2 + Y₂− Y₃=0

...

X_n−k − X_n + Y_n−k − Y_n−k+1 = 0

(k − m_s) – X₁ − X₂ −...− X_k − Y₁ = 0

X_n−k+1 + X_n−k+2 +...+ X_n +Y_n−k+1 − (k − m_s) = 0

我们把等式左边加起来，发现恰好为0，满足流量平衡。

把上述每个等式看成一个点，若在第j个等式有-Y_i，第k个等式有Y_i，则向连一条从j到k的容量为k–m_s–m_e的边，若在第j个等式有-X_i，在第k个等式有X_i则连一条j到k容量为1的边，然后对于正的常数c出现在第j个等式，从原点到j连一条容量为c的边，若对于负的常数-c出现在第k个等式，从k向汇点连一条容量为c的边，因为这一题还需要有费用，所以在相应的边上加上费用限制即可，最后跑一遍费用流。

Source：

#include 
   
   
    
    
#include 
    
    
     
     
#include 
     
     
      
      
#include 
      
      
       
       
#include 
       
       
         #include 
        
          #include 
         
           using namespace std; inline char read() { static const int IN_LEN = 1024 * 1024; static char buf[IN_LEN], *s, *t; if (s == t) { t = (s = buf) + fread(buf, 1, IN_LEN, stdin); if (s == t) return -1; } return *s++; } template 
          
            inline bool R(T &x) { static char c; static bool iosig; for (iosig = false, c = read(); !isdigit(c); c = read()) { if (c == -1) return false; if (c == '-') iosig = true; } for (x = 0; isdigit(c); c = read()) x = (x << 3) + (x << 1) + (c ^ '0'); if (iosig) x = -x; return true; } const int OUT_LEN = 1024 * 1024; char obuf[OUT_LEN], *oh = obuf; inline void writechar(char c) { if (oh == obuf + OUT_LEN) fwrite(obuf, 1, oh - obuf, stdout), oh = obuf; *oh++ = c; } template 
           
             inline void W(T x) { static int buf[30], cnt; if (!x) writechar(48); else { if (x < 0) writechar('-'), x = -x; for (cnt = 0; x; x /= 10) buf[++cnt] = x % 10 + 48; while (cnt) writechar(buf[cnt--]); } } inline void flush() { fwrite(obuf, 1, oh - obuf, stdout); } const long long INF = 1e18; const int MAXN = 1000000 + 10; int s, t; long long ans1, sum, ans; long long dis[MAXN], f[MAXN]; bool vis[MAXN]; struct node { int to, rev; long long w, price; node(int to, long long w, long long price, int rev) : to(to), w(w), price(price), rev(rev) {} node() {} }; vector 
            
              edge[MAXN]; inline void create(int x, int y, long long w, long long c) { edge[x].push_back(node(y, w, c, edge[y].size())); if (x == y) edge[x].back().rev++; edge[y].push_back(node(x, 0, -c, edge[x].size() - 1)); } inline long long dfs(int cur, long long low, int t) { if (cur == t) return ans1 += low * sum, low; vis[cur] = true; int res = 0, delta; for (int p = 0; p < edge[cur].size(); ++p) { node &e = edge[cur][p]; if (e.w > 0 && !e.price && !vis[e.to]) { delta = dfs(e.to, min(low - res, e.w), t); if (delta) { e.w -= delta, edge[e.to][e.rev].w += delta, res += delta; if (res == low) break; } } } return res; } inline bool spfa_slf(int s, int t) { for (int i = s; i <= t; ++i) dis[i] = -INF; deque 
             
               q; dis[t] = 0, q.push_back(t); while (!q.empty()) { int j = q.front(); q.pop_front(); for (int p = 0; p < edge[j].size(); ++p) { node &e = edge[j][p]; if (edge[e.to][e.rev].w && dis[j] - e.price > dis[e.to]) { dis[e.to] = dis[j] - e.price; if (!q.empty() && dis[e.to] > dis[q.front()]) q.push_front(e.to); else q.push_back(e.to); } } } for (int i = s; i <= t; ++i) for (int p = 0; p < edge[i].size(); ++p) edge[i][p].price += dis[edge[i][p].to] - dis[i]; sum += dis[s]; return dis[s] > -INF; } inline void cost_flow(int s, int t) { while (spfa_slf(s, t)) do memset(vis, false, sizeof(vis)); while(dfs(s, INF, t)); } int n, k, ms, me; int ss[MAXN], ee[MAXN]; int main() { // freopen("in.in", "r", stdin); R(n), R(k), R(ms), R(me), me = k - me; for (int i = 1; i <= n; ++i) R(ss[i]); for (int i = 1; i <= n; ++i) R(ee[i]), ans += ee[i]; s = 0, t = n - k + 2 + 1; create(s, 1, me, 0); create(n - k + 2, t, me, 0); for (int i = 1; i <= n; ++i) { int u = max(1, i - k + 1), v = min(i + 1, n - k + 2); create(u, v, 1, ss[i] - ee[i]); } for (int i = 1; i <= n - k + 1; ++i) create(i, i + 1, me - ms, 0); cost_flow(s, t); W(ans1 + ans); flush(); return 0; }