P3358 最长k可重区间集问题

• 问题描述：给定n个开区间，从中选取若干个，使得每个数的被覆盖次数不超过k，且选取的区间长度和最大；
• 建图方法：将区间上的端点视为图中的点，对于每个区间，左端点向右端点连权值为区间长度，容量为1（保证了每个区间只选一次）的边，源点向所有的左端点连权值为0，容量为1的边，所有的右端点向汇点连权值为0，容量为1的边，对于每个不重叠的区间，左边区间的右端点向右边区间的左端点连权值为0，容量为1的边，这样就形成了若干组串联关系；现在考虑并联，流并联的两条链等价于覆盖次数+2，现在问题就变成了选取k条权值最大的串联链，因此再建立一个super源向源点连容量为k的边；跑最大费用最大流即可。
• 总结：依旧延续了流表达限制条件的方法，将区间拆成左右端点限制了区间的选取次数，通过类比于电路图中的串并联关系（串联电路覆盖为1，i条串联最大覆盖数为i）将覆盖不超过k这一约束条件转换成不超过k条“串联电路”。

#include <iostream>
#include <map>
#include <cstdio>
#include <algorithm>
#include <queue>
#include <cstring>
using namespace std;
const int maxn=1105,maxm=500*500*2,inf=0x3f3f3f3f;
int cnt=1,last[maxn],cur[maxn],pos=1,inq[maxn],st=1100,n,k,l[maxn],r[maxn],len,ml,ans,d[maxn];
bool vis[maxn];
struct edge{
    int v,w,next,f;
}e[maxm];
inline void add(int u,int v,int w,int f)
{
    e[++cnt].v=v;
    e[cnt].next=last[u];
    e[cnt].w=w;
    e[cnt].f=f;
    last[u]=cnt;
}
bool spfa()
{
    queue<int>q;
    memset(inq,0,sizeof(inq));
    for(int i=1;i<=st;i++) d[i]=inf;
    d[0]=0;
    q.push(0);
    while(!q.empty())
    {
        int u=q.front();q.pop();
        inq[u]=0;
        for(int i=last[u];i;i=e[i].next)
        {
            int v=e[i].v,w=e[i].w,f=e[i].f;
            if(d[v]>d[u]+w&&f)
            {
                d[v]=d[u]+w;
                if(!inq[v])
                {
                    inq[v]=1;
                    q.push(v);
                }
            }
        }
    }
    return d[st]!=inf;
}
int dfs(int u,int dis)
{
    vis[u]=1;
    if(u==st||!dis) return dis;
    for(int i=cur[u];i;i=e[i].next)
    {
        cur[u]=e[i].next;
        int v=e[i].v,w=e[i].w,f=e[i].f;
        if(d[v]==d[u]+w&&f&&!vis[v])
        {
            int di=dfs(v,min(dis,f));
            if(di>0)
            {
                e[i].f-=di;
                e[i^1].f+=di;
                return di;
            }
        }
    }
    return 0;
}
map<int,int>p;
void dinic()
{
    while(spfa())
    {
        for(int i=0;i<=st;i++)
        cur[i]=last[i];
        memset(vis,0,sizeof(vis));
        while(int mi=dfs(0,inf))
        {
            ans+=d[st]*mi;
            ml+=mi;;
            memset(vis,0,sizeof(vis));
        }
    }
}

int main()
{

   cin>>n>>k;
   add(0,1,0,k);
   add(1,0,0,0);
   for(int i=1;i<=n;i++)
   {
       cin>>l[i]>>r[i];
       len=r[i]-l[i];
       if(!p[l[i]]) p[l[i]]=++pos;
       if(!p[r[i]]) p[r[i]]=++pos;
       add(p[l[i]],p[r[i]],-len,1);
       add(p[r[i]],p[l[i]],len,0);
       add(1,p[l[i]],0,1);
       add(p[l[i]],1,0,0);
       add(p[r[i]],st,0,1);
       add(st,p[r[i]],0,0);
   }
   for(int i=1;i<=n;i++)
   for(int j=1;j<=n;j++)
   if(r[i]<=l[j])
   {
       add(p[r[i]],p[l[j]],0,1);
       add(p[l[j]],p[r[i]],0,0);
   }
   dinic();
   cout<<-ans;
   return 0;
}