「SNOI2019」通信 分治建图

根据题意 每个点可以直接与S,T相连 也可以和前面的哨站相连 暴力建边的话 有n²条边

要用分治优化建边:

类似于归并排序 先对每一层分为左半边与右半边 对每一半都拿出来先排序去重后 直接排成一条链建边

if (l == r) {
                return ;
        }
        int mid = (l + r) >> ;
        solve(l, mid), solve(mid + , r);
        int cnt = ;
        for (int i = l; i <= r; i++) {
                aa[++cnt] = a[i];
        }
        sort(aa + , aa +  + cnt);
        cnt = unique(aa + , aa + cnt + ) - aa - ;
        int now = MCMF::MAXP;
        for (int i = ; i < cnt; i++) {
                MCMF::addedge(now + i, now + i + , inf, aa[i + ] - aa[i]);
                MCMF::addedge(now + i + , now + i, inf, aa[i + ] - aa[i]);
        }

然后对于区间内的每个数 前半边的出后半边的入

        for (int i = l; i <= r; i++) {
                int aim = lower_bound(aa + , aa + cnt + , a[i]) - aa;
                if (i <= mid) {
                        MCMF::addedge(now + aim, i + n, , );
                } else {
                        MCMF::addedge(i, now + aim, , );
                }
        }

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef long long JQK;
const int inf = INT_MAX / ;int n, a[], aa[];
void solve(int l, int r) {
        if (l == r) {
                return ;
        }
        int mid = (l + r) >> ;
        solve(l, mid), solve(mid + , r);
        int cnt = ;
        for (int i = l; i <= r; i++) {
                aa[++cnt] = a[i];
        }
        sort(aa + , aa +  + cnt);
        cnt = unique(aa + , aa + cnt + ) - aa - ;
        int now = MCMF::MAXP;
        for (int i = ; i < cnt; i++) {
                MCMF::addedge(now + i, now + i + , inf, aa[i + ] - aa[i]);
                MCMF::addedge(now + i + , now + i, inf, aa[i + ] - aa[i]);
        }
        for (int i = l; i <= r; i++) {
                int aim = lower_bound(aa + , aa + cnt + , a[i]) - aa;
                if (i <= mid) {
                        MCMF::addedge(now + aim, i + n, , );
                } else {
                        MCMF::addedge(i, now + aim, , );
                }
        }
        MCMF::MAXP += cnt;
}
int main() {
        int x;
        int s, t;
        ll W;
        RR(n), RR(x);
        W = x;
        MCMF::MAXP =  * n + , s =  * n + , t =  * n + ;
        for (int i = ; i <= n; i++) {
                RR(a[i]);
                MCMF::addedge(s, i, , ), MCMF::addedge(i, t, , W), MCMF::addedge(i + n, t, , );
        }
        solve(, n);
        MCMF::init(s, t);
        cout << MCMF::MCMF() << "\n";
        return ;
}