\(\mathcal{Description}\)

Link.

有 \(n\) 个位置，从左至右编号 \(1\sim n\)。在第 \(i\) 个位置放一座塔的代价为 \(c_i\)，一个位置可以放任意数量的塔。给定 \(m\) 个要求，第 \(i\) 个表示 \([l_i,r_i]\) 内至少有 \(d_i\) 座塔。求最小的代价和。

\(n\le10^3\)，其余参数 \(\le10^4\)。

\(\mathcal{Solution}\)

经历了逝量的 whk 学习，我学会了背套路。（

原问题可以写成线规：

\[\begin{aligned}
\min~~~~&z=\sum_{i=1}^nc_i(x_i-x_{i-1})\\
\text{s.t.}~~~~&x_i\ge x_{i-1}\\
&x_{r_i}-x_{l_i-1}\ge d_i
\end{aligned},
\]

而它显然等价于：

\[\begin{aligned}
\min~~~~z'=\sum_{i=1}^n(c_i-c_{i+1})x_i+\sum_{i=1}^nI\max\{x_{i-1}-x_i,0\}+\sum_{i=1}^mI\max\{d_i-x_{r_i}+x_{l_i-1},0\}
\end{aligned}.
\]

其中 \(I\) 是一个足够大的常数。

我们的论文套路是，对于以下线规：

\[\min~~~~z=\sum_ub_ux_u+\sum_{\lang u,v\rang}c_{uv}\max\{0,x_v-x_u-w_{uv}\},
\]

其对应最小费用流模型：

\(b_u\) 表示 \(u\) 点的流出量-流入量。

为了流量守恒，对于 \(b_u>0\) 的 \(u\)，连 \(\lang S,u,b_u,0\rang\)；

对于 \(b_u<0\) 的 \(u\)，连 \(\lang u,T,-b_u,0\rang\)；
\(c_{uv}\) 表示 \(\lang u,v\rang\) 的容量，\(w_{uv}\) 表示 \(\lang u,v\rang\) 的费用，所以连边 \(\lang u,v,c_{uv},w_{uv}\rang\)。

该图的最小费用最大流的费用的相反数就是答案。

点的系数出减入，max 外是容量，max 内是费用，减数连向被减数，最小费用相反数！

但 whk 是落后的，我们来证明。令 \(f_{uv}\) 表示 \(\lang u,v\rang\) 的流量，考虑在上述条件下的最小费用流形式：

\[\begin{aligned}
\min~~~~&z=\sum_{\lang u,v\rang}c_{uv}f_{uv}\\
\text{s.t.}~~~~&-f_{uv}\ge-c_{uv}\\
&\sum_v f_{vu}-\sum_vf_{uv}=-b_u
\end{aligned}.
\]

令 \(p_{uv}\) 表示第一类限制的对偶，\(q_u\) 表示第二类限制的对偶，得到：

\[\begin{aligned}
\max~~~~&z=\sum_u -b_uq_u+\sum_{\lang u,v\rang}-c_{uv}p_{uv}\\
\text{s.t.}~~~~&q_v-q_u-p_{uv}\le w_{uv}
\end{aligned}.
\]

那么这个 \(p_{uv}\) 非常的自由，直接给它定成最优，所以

\[\min~~~~z=\sum_ub_uq_u+\sum_{\lang u,v\rang}c_{uv}\max\{0,q_v-q_u-w_{uv}\}
\]

的相反数就是答案。 \(\square\)

\(\mathcal{Code}\)

总之得写势能 Dijkstra。

/*+Rainybunny+*/

#include <bits/stdc++.h>

#define rep(i, l, r) for (int i = l, rep##i = r; i <= rep##i; ++i)

#define per(i, r, l) for (int i = r, per##i = l; i >= per##i; --i)

typedef std::pair<int, int> PII;

#define fi first

#define se second

const int MAXN = 1e3, MAXM = 1e4, IINF = 0x3f3f3f3f;

int n, m, c[MAXN + 5], l[MAXM + 5], r[MAXM + 5], d[MAXM + 5];

namespace FlowGraph {

const int MAXND = MAXN + 3, MAXEG = MAXN * 2 + MAXM;

int ecnt = 1, S, T, head[MAXND + 5], curh[MAXND + 5];

int hig[MAXND + 5], dis[MAXND + 5];

bool instk[MAXND + 5];

struct Edge { int to, flw, cst, nxt; } graph[MAXEG * 2 + 5];

inline void link(const int s, const int t, const int f, const int c) {

    graph[++ecnt] = { t, f, c, head[s] }, head[s] = ecnt;

    graph[++ecnt] = { s, 0, -c, head[t] }, head[t] = ecnt;

}

inline bool spfa() {

    static bool inq[MAXND + 5]; static std::queue<int> que;

    rep (i, 0, T) dis[i] = IINF;

    dis[S] = 0, que.push(S);

    while (!que.empty()) {

        int u = que.front(); que.pop(), inq[u] = false;

        for (int i = head[u], v; i; i = graph[i].nxt) {

            if (graph[i].flw && dis[u] + graph[i].cst < dis[v = graph[i].to]) {

                dis[v] = dis[u] + graph[i].cst;

                if (!inq[v]) inq[v] = true, que.push(v);

            }

        }

    }

    return dis[T] != IINF;

}

inline bool dijkstra() {

    static std::priority_queue<PII, std::vector<PII>, std::greater<PII> > heap;

    rep (i, 0, T) hig[i] += dis[i], dis[i] = IINF;

    heap.push({ dis[S] = 0, S });

    while (!heap.empty()) {

        PII p(heap.top()); heap.pop();

        if (dis[p.se] != p.fi) continue;

        for (int i = head[p.se], v; i; i = graph[i].nxt) {

            int d = p.fi + graph[i].cst + hig[p.se] - hig[v = graph[i].to];

            if (graph[i].flw && dis[v] > d) heap.push({ dis[v] = d, v });

        }

    }

    return dis[T] != IINF;

}

inline PII augment(const int u, int iflw) {

    if (u == T) return { iflw, 0 };

    PII ret(0, 0); instk[u] = true;

    for (int &i = curh[u], v; i; i = graph[i].nxt) {

        if (graph[i].flw && !instk[v = graph[i].to]

          && dis[v] == dis[u] + hig[u] - hig[v] + graph[i].cst) {

            PII tmp(augment(v, std::min(iflw, graph[i].flw)));

            ret.fi += tmp.fi, ret.se += graph[i].cst * tmp.fi + tmp.se;

            iflw -= tmp.fi, graph[i].flw -= tmp.fi, graph[i ^ 1].flw += tmp.fi;

            if (!iflw) break;

        }

    }

    if (ret.fi) instk[u] = false;

    return ret;

}

inline PII dinic() {

    PII ret(0, 0);

    for (spfa(); dijkstra(); ) {

        rep (i, 0, T) curh[i] = head[i], instk[i] = false;

        PII tmp(augment(S, IINF));

        ret.fi += tmp.fi, ret.se += tmp.se;

    }

    return ret;

}

} using namespace FlowGraph;

int main() {

    scanf("%d %d", &n, &m);

    rep (i, 1, n) scanf("%d", &c[i]);

    rep (i, 1, m) scanf("%d %d %d", &l[i], &r[i], &d[i]);

    S = n + 1, T = n + 2;

    rep (i, 0, n) {

        if (c[i] > c[i + 1]) link(S, i, c[i] - c[i + 1], 0);

        else link(i, T, c[i + 1] - c[i], 0);

    }

    rep (i, 0, n - 1) link(i + 1, i, IINF, 0);

    rep (i, 1, m) link(r[i], l[i] - 1, IINF, -d[i]);

    printf("%d\n", -dinic().se);

    return 0;

}

巴特西

Solution -「ZJOI 2013」「洛谷 P3337」防守战线

\(\mathcal{Description}\)

\(\mathcal{Solution}\)

\(\mathcal{Code}\)

最新文章

热门文章