【网络流24题】最长k可重区间集问题

（网络流24题大多需要spj，所以需要一个有spj的oj，本系列代码均在www.oj.swust.edu.cn测试通过）
这道题可以说是经典中的经典了，我们可以这样想，由于要求包含一个点的区间最多不超过k个，那么我们就可以只给这个序列k的流量，一旦选择了一条线段就在这条线段的区间范围内减少一点流量，这样就能保证没有任何一点上面的区间会超过k个，那么具体怎么实现呢，由于我们不清楚区间的大小，所以我们要先离散化区间的左右端点，然后每一个位置向它下一个位置连一条INF，费用为0（这条线段不选），然后每一条线段的左端点向其右端点连一条1，费用为区间长度，跑一边最大费用最大流就行了。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<iostream>
#include<iomanip>
#include<ctime>
#include<cmath>
#include<algorithm>
using namespace std;
#define INF 100000000
struct bian
{
    int l,r,f,v;
}a[1000000];
int fir[1000000];
int nex[1000000];
int fro[10000];
int tot=1;
int s=0,t=9999;
int dis[10000];
bool pd[10000];
void add_edge(int l,int r,int f,int v)
{
    a[++tot].l=l;
    a[tot].r=r;
    a[tot].f=f;
    a[tot].v=v;
    nex[tot]=fir[l];
    fir[l]=tot;
}
bool spfa()
{
    memset(dis,0x8f,sizeof(dis));
    static int dui[1000000];
    int top=1,my_final=2;
    dui[top]=s;
    pd[s]=true;
    dis[s]=0;
    while(top<my_final)
    {
        int u=dui[top++];
        for(int o=fir[u];o;o=nex[o])
        {
            if(a[o].f && dis[a[o].r]<dis[u]+a[o].v)
            {
                dis[a[o].r]=dis[u]+a[o].v;
                fro[a[o].r]=o;
                if(!pd[a[o].r])
                {
                    pd[a[o].r]=true;
                    dui[my_final++]=a[o].r;
                }
            }
        }
        pd[u]=false;
    }
    if(dis[t]==0x8f8f8f8f) return false;
    return true;
}
int cost=0;
void add_flow()
{
    int mid=t;
    int temp=2147483647;
    while(mid!=s)
    {
        temp=min(temp,a[fro[mid]].f);
        mid=a[fro[mid]^1].r;
    }
    mid=t;
    cost+=temp*dis[t];
    while(mid!=s)
    {
        a[fro[mid]].f-=temp;
        a[fro[mid]^1].f+=temp;
        mid=a[fro[mid]^1].r;
    }
}
struct point
{
    int x,tim,num;
}xian[1000000];
bool cmp1(point a,point b)
{
    return a.x<b.x;
}
bool cmp2(point a,point b)
{
    return a.tim<b.tim;
}
int top=0;
int main()
{
    int n,k;
    scanf("%d%d",&n,&k);
    for(int i=1;i<=n;i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        xian[++top].x=x;
        xian[top].tim=top;
        xian[++top].x=y;
        xian[top].tim=top;
    }
    sort(xian+1,xian+1+top,cmp1);
    int wz=0;
    for(int i=1;i<=top;i++)
    {
        if(i==1 || xian[i].x!=xian[i-1].x) wz++;
        xian[i].num=wz;
    }
    for(int i=1;i<wz;i++)
    {
        add_edge(i,i+1,INF,0);
        add_edge(i+1,i,0,0);
    }
    add_edge(s,1,k,0);
    add_edge(1,s,0,0);
    add_edge(wz,t,k,0);
    add_edge(t,wz,0,0);
    sort(xian+1,xian+1+top,cmp2);
    for(int i=1;i<=top;i+=2)
    {
        int val=xian[i+1].x-xian[i].x;
        add_edge(xian[i].num,xian[i+1].num,1,val);
        add_edge(xian[i+1].num,xian[i].num,0,-val);
    }
    while(spfa()) add_flow();
    cout<<cost;
    return 0;
}