最长 k 可重区间集问题

本题其实就是求最大权不相交路径，可以用费用流来建模。

方法如下：

首先建立源点ST和汇点ED还要建立另一个间接汇点SS。由ST向SS连流量为k，权值为0的边。要把一条线段拆成两个端点来处理。再按照左端点排序。由SS向每个做端点连流量为1，权值为0的边，如果选显然有流量流过。这时再把所有点的右端点向ED连流量为1，权值为0的边。对于线段i，如果其右端点小于其右边点的左端点，说明两者显然不会重合，这时就由其的右端点向满足条件的点的左端点连权值为0，流量为1的边。最后每个线段的左端点向右端点连权值为长度，流量为1的边，跑一边最大费用最大流就是答案.

代码如下：

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <vector>
#include <queue>

#define M 101011
#define INF 0x7fffffff/3
#define ST 2001
#define ED 2003
#define SS 2002

using namespace std;

int n,m;
int num=0,ans=0;
int cost[M],p[M],flag[M],pre[M];
struct node
{
    int a;
    int b;
    int w;
    int f;
    int next;
}str1[M];
int q[4*M];

void add(int x,int y,int z,int fl)
{
    str1[num].a=x;
    str1[num].b=y;
    str1[num].w=z;
    str1[num].f=fl;
    str1[num].next=p[x];
    p[x]=num++;

    swap(x,y);

    str1[num].a=x;
    str1[num].b=y;
    str1[num].w=-z;
    str1[num].f=0;
    str1[num].next=p[x];
    p[x]=num++;
}

bool spfa(int src)
{
    memset(flag,0,sizeof(flag));
    for(int i=0;i<=ED;i++)
        cost[i]=INF;
    cost[src]=0;
    flag[src]=1;
    pre[src]=-1;
    int l=0,r=0;
    q[r++]=src;
    while(l<r)
    {
        int k=q[l++];
        flag[k]=0;
        for(int i=p[k];i!=-1;i=str1[i].next)
        {
            int v=str1[i].b;
            if(cost[v]>cost[k]+str1[i].w&&str1[i].f>0)
            {
                cost[v]=cost[k]+str1[i].w;
                pre[v]=i;
                if(!flag[v])
                {
                    flag[v]=1;
                    q[r++]=v;
                }
            }
        }
    }
    if(cost[ED]==INF)
        return false;
    return true;
}

void start()
{
    int minx=INF;
    for(int i=pre[ED];i!=-1;i=pre[str1[i^1].b])
        minx=min(minx,str1[i].f);
    for(int i=pre[ED];i!=-1;i=pre[str1[i^1].b])
    {
        str1[i].f-=minx;
        str1[i^1].f+=minx;
    }
    ans+=cost[ED]*minx;
}

void spfaflow()
{
    ans=0;
    while(spfa(ST))
        start();
}

struct node1
{
    int fx;
    int fy;
    int fw;
}str2[M];

bool cmp(const node1 &a,const node1 &b)
{
    if(a.fx==b.fx)
        return a.fy<b.fy;
    return a.fx<b.fx;
}

void init()
{
    scanf("%d%d",&n,&m);
    memset(p,-1,sizeof(p));
    for(int i=1;i<=n;i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        str2[i].fx=x;
        str2[i].fy=y;
        str2[i].fw=y-x;
    }
    sort(str2+1,str2+n+1,cmp);
    //for(int i=1;i<=n;i++)
    //    cout<<str2[i].fx<<" "<<str2[i].fy<<" "<<str2[i].fw<<endl;
    add(ST,SS,0,m);
    for(int i=1;i<=n;i++)
    {
        add(SS,i,0,1);
        add(i+n,ED,0,1);
        add(i,i+n,-str2[i].fw,1);
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=i+1;j<=n;j++)
        {
            if(str2[i].fy<=str2[j].fx)
                add(i+n,j,0,1);
        }
    }
}

int main()
{
    //freopen("T21.in","r",stdin);
    //freopen("T21.out","w",stdout);

    init();
    spfaflow();
    printf("%d\n",-ans);

    return 0;
}