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

(http://www.elijahqi.win/2017/12/12/loj-6014%E3%80%8C%E7%BD%91%E7%BB%9C%E6%B5%81-24-%E9%A2%98%E3%80%8D%E6%9C%80%E9%95%BF-k-%E5%8F%AF%E9%87%8D%E5%8C%BA%E9%97%B4%E9%9B%86/%20%E2%80%8E)
题目描述

给定实直线 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 l​i​​ 和 ri r_i r​i​​，表示开区间的左右端点坐标，注意可能有 li>ri l_i > r_i l​i​​>r​i​​，此时请将其交换。
输出格式

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

4 2
1 7
6 8
7 10
9 13

样例输出

15

数据范围与提示

1≤n≤500,1≤k≤3 1 \leq n \leq 500, 1 \leq k \leq 3 1≤n≤500,1≤k≤3

构建一个图 保证满流是满足条件的

首先 源向起点连边权为k的权值 然后 最末端向汇点连边权为Inf 费用为0的边

每个区间从左端点向右端点连边 权值为1 表示这个区间只能使用它覆盖一次 费用为这个区间的长度 然后跑最大费用最大流即可求出答案

#include<queue>
#include<vector>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define inf 0x3f3f3f3f
#define N 1100
using namespace std;
vector<int> aq;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0;char ch=gc();
    while(ch<'0'||ch>'9') ch=gc();
    while(ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=gc();}
    return x;
}
int h[N],pre[N],flag[N],f[N],path[N],l[N],r[N],n,k,num=1,T;
struct node{
    int y,z,next,c;
}data[N*N];
inline void insert1(int x,int y,int z,int c){
    data[++num].y=y;data[num].z=z;data[num].next=h[x];h[x]=num;data[num].c=c;
    data[++num].y=x;data[num].z=0;data[num].next=h[y];h[y]=num;data[num].c=-c;
}
inline bool spfa(){
    memset(flag,0,sizeof(flag));memset(f,128,sizeof(f));queue<int>q;q.push(0);flag[0]=1;f[0]=0;memset(pre,-1,sizeof(pre));
    while(!q.empty()){
        int x=q.front();q.pop();flag[x]=0;
        for (int i=h[x];i;i=data[i].next){
            int y=data[i].y,z=data[i].z,c=data[i].c;
            if (f[x]+c>f[y]&&z){
                f[y]=f[x]+c;pre[y]=x;path[y]=i;
                if (!flag[y]) flag[y]=1,q.push(y);
            }
        }
    }if (pre[T]==-1) return 0;else return 1;
}
int main(){
    freopen("loj.in","r",stdin);
    n=read();k=read();
    for (int i=1;i<=n;++i){
        l[i]=read();r[i]=read();if(l[i]>r[i]) swap(l[i],r[i]);aq.push_back(l[i]);aq.push_back(r[i]);
    }sort(aq.begin(),aq.end());insert1(0,1,k,0);int nn=unique(aq.begin(),aq.end())-aq.begin();T=nn+1;
    for (int i=1;i<=n;++i){
        int cost=r[i]-l[i];
        int ll=lower_bound(aq.begin(),aq.begin()+nn,l[i])-aq.begin()+1,rr=lower_bound(aq.begin(),aq.begin()+nn,r[i])-aq.begin()+1;
        insert1(ll,rr,1,cost);
    }insert1(nn,T,inf,0);
    for (int i=1;i<nn;++i) insert1(i,i+1,inf,0);int ans=0;
    while(spfa()){
        int minn=inf,now=T;
        while(now) minn=min(minn,data[path[now]].z),now=pre[now];now=T;
        while(now){ans+=data[path[now]].c*minn;data[path[now]].z-=minn;data[path[now]^1].z+=minn;now=pre[now];}
    }printf("%d",ans);
    return 0;
}