「网络流 24 题」最长 k 可重区间集

#6014. 「网络流 24 题」最长 k 可重区间集

内存限制：256 MiB 时间限制：1000 ms 标准输入输出

题目类型：传统 评测方式：文本比较

上传者： 匿名

题目描述

给定实直线 L L L 上 n n n 个开区间组成的集合 I I I，和一个正整数 k k k，试设计一个算法，从开区间集合 I I I 中选取出开区间集合 S⊆I S \subseteq I S⊆I，使得在实直线 L L L 的任何一点 x x x，S S S 中包含点 x x x 的开区间个数不超过 k k k。且 ∑z∈S∣z∣ \sum\limits_{z \in S} | z | z∈S∑​∣z∣ 达到最大。这样的集合 S S S 称为开区间集合 I I I 的最长 k k k 可重区间集。∑z∈S∣z∣ \sum\limits_{z \in S} | z | z∈S∑​∣z∣ 称为最长 k k k 可重区间集的长度。

对于给定的开区间集合 I I I 和正整数 k k k，计算开区间集合 I I I 的最长 k k k 可重区间集的长度。

输入格式

文件的第 1 1 1 行有 2 2 2 个正整数 n n n 和 k k k，分别表示开区间的个数和开区间的可重迭数。
接下来的 n n n 行，每行有 2 2 2 个整数 li l_i li​ 和 ri r_i ri​，表示开区间的左右端点坐标，注意可能有 li>ri l_i > r_i li​>ri​，此时请将其交换。

输出格式

输出最长 k k k 可重区间集的长度。

样例

样例输入

4 2
1 7
6 8
7 10
9 13

样例输出

15

 

题意：区间覆盖费用流。注意垂直的时候处理。可以先把所有的坐标*2，然后左端点-1，不影响最终答案！

 

#include<bits/stdc++.h>
using namespace std;
#define rep(i,a,n) for (int i=a;i<n;i++)
#define per(i,a,n) for (int i=n-1;i>=a;i--)
#define oo cout<<"!!!"<<endl;
typedef long long ll;
typedef unsigned long long ull;
#define ms(s) memset(s, 0, sizeof(s))
const int inf = 0x3f3f3f3f;
//head

const int maxn = 1e6+11;

int ver[maxn],edge[maxn],cost[maxn],nxt[maxn],head[maxn],cnt;
int _edge[maxn];
int d[maxn],incf[maxn],pre[maxn],v[maxn];
int n,m,s,t,ans,maxflow;

void add(int x,int y,int z,int c)
{
    ver[++cnt] = y, edge[cnt] = z;cost[cnt] = c; nxt[cnt]  = head[x];head[x] = cnt;
    ver[++cnt] = x; edge[cnt] = 0;cost[cnt] = -c;nxt[cnt] = head[y];head[y] = cnt;
}

bool spfa()
{
    queue<int>q;
    memset(d,0xcf,sizeof d);
    ms(v);
    q.push(s);d[s] = 0;v[s] = 1;
    incf[s] = inf;

    while(q.size())
    {
        int x = q.front();q.pop();v[x] = 0;
        for(int i = head[x];i;i=nxt[i])
        {
            if(!edge[i])continue;
            int y = ver[i];
            if(d[y] < d[x] + cost[i])
            {
                d[y] = d[x] + cost[i];
                incf[y] = min(incf[x],edge[i]);
                pre[y] = i;
                if(!v[y])q.push(y),v[y] = 1;
            }
        }
    }
    if(d[t] == 0xcfcfcfcf)return false;
    return true;
}


void update()
{
    int x = t;
    while(x != s)
    {
        int i = pre[x];
        edge[i] -= incf[t];
        edge[i^1] += incf[t];
        x = ver[i^1];
    }
    maxflow += incf[t];
    ans += d[t] * incf[t];
}

int w[maxn],a[maxn],l[maxn],r[maxn];

int main()
{
    int n,k;
    cin >> n >> k;
    int tot = 0;
    cnt = 1;
    rep(i,1,n+1)
    {
        int x1,x2,y1,y2;
        scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
        w[i] = sqrt((ll)(x1-x2)*(x1-x2) + (ll)(y1-y2)*(y1-y2));
        l[i] = x1 << 1, r[i] = x2 << 1;
        if(x1 == x2)l[i]--;
        if(l[i] > r[i])swap(l[i],r[i]);
        a[++tot] = l[i],a[++tot] = r[i];
    }
    sort(a+1,a+1+tot),tot = unique(a+1,a+1+tot)-a-1;
    s = tot+1,t = s+1;
    add(s,1,k,0),add(tot,t,k,0);
    rep(i,1,tot)add(i,i+1,inf,0);
    rep(i,1,n+1)
    {
        int L = lower_bound(a+1,a+1+tot,l[i]) - a;
        int R = lower_bound(a+1,a+1+tot,r[i]) - a;
        add(L,R,1,w[i]); 
    } 
    while(spfa())update();   
    cout << ans << endl;

}