网络流24题——15.最长k可重区间集问题

题目链接

https://www.luogu.org/problemnew/show/P3358

最长k可重区间集问题

先离散化，要去重，设去完重之后点数为size，每个点向后面一个点连流量为k，费用为0的边，当然，S向1连，size向T连，然后每条线段的左端点向右端点连流量为1，费用为区间大小的边，然后跑最大费用最大流就行了。
我们考虑每个点有两个决策，一个是走到下一个点，表示不取任何一条线段，贡献为0。另一个是走到这条线段的右端点，表示选择这条线段，显然边都是从左向右的，所以一条增广路中选择的线段都不相交，然后最多只能选k条增广路，所以S到1流量控制为k，然后中间的点保证每次都能流就行，所以流量设为k。

// luogu-judger-enable-o2
#include<cstdio>
#include<algorithm>
#include<cctype>
#include<cstring>
#include<iostream>
#define LL long long
#define INF (2139062143)
#define N (1005)
using namespace std;
int n,k,tot,S,T,all,h,t,ans;
int head[N],cur[N],co[N],nxt[N<<1],f[N<<1],fare[N<<1],a[N<<1],from[N],que[N<<2],dis[N],data[N<<1],fnd[N<<1]; 
bool vis[N];
template <typename T> void read(T&t) {
    t=0;
    bool fl=true;
    char p=getchar();
    while (!isdigit(p)) {
        if (p=='-') fl=false;
        p=getchar();
    }
    do {
        (t*=10)+=p-48;p=getchar();
    }while (isdigit(p));
    if (!fl) t=-t;
}
struct node{
    int data,id;
}q[N];
inline bool cmp(node a,node b){
    return a.data<b.data;
}
inline void add(int x,int y,int z,int t){
    //printf("%d %d %d %d\n",x,y,z,t);
    a[++tot]=y,f[tot]=z,fare[tot]=t,nxt[tot]=head[x],head[x]=tot;
}
inline int spfa(){
    memset(dis,-127,sizeof(dis));
    h=t=0;
    que[++t]=S;
    vis[S]=1;
    dis[S]=0;
    while (h<t){
        int u=que[++h];
        for (int p=head[u];~p;p=nxt[p]){
            //printf("%d %d %d %d %d\n",u,a[p],dis[u],dis[a[p]],fare[p]);
            if (f[p]&&dis[u]+fare[p]>dis[a[p]]){
                dis[a[p]]=dis[u]+fare[p];
                from[a[p]]=p;
                if (!vis[a[p]]){
                    vis[a[p]]=1;
                    que[++t]=a[p];
                }
            }
        }
        vis[u]=0;
    } 
    //for (int i=1;i<=all+2;i++) printf("%d ",dis[i]);
    //puts("");
    return dis[T];
}
int main(){
    read(n),read(k);
    for (int i=1;i<=2*n;i++){
        read(data[i]);
        q[i].data=data[i];
        q[i].id=i;
    }
    sort(q+1,q+2*n+1,cmp);
    fnd[q[1].id]=++all;
    for (int i=2;i<=2*n;i++){
        if (q[i].data==q[i-1].data) fnd[q[i].id]=all;
        else fnd[q[i].id]=++all;
    }

    tot=1,S=all+1,T=all+2;
    memset(head,-1,sizeof(head));
    for (int i=1;i<=n;i++){
        add(fnd[i*2-1],fnd[i*2],1,data[i*2]-data[i*2-1]);
        add(fnd[i*2],fnd[i*2-1],0,-data[i*2]+data[i*2-1]);
    }
    for (int i=1;i<all;i++){
        add(i,i+1,k,0);
        add(i+1,i,0,0);
    }
    add(S,1,k,0);
    add(1,S,0,0);
    add(all,T,k,0);
    add(T,all,0,0);
    while (spfa()!=dis[0]){
        if (dis[T]==0) break;
        ans+=dis[T];
        int now=T;
        while (now!=S){
            f[from[now]]--;
            f[from[now]^1]++;
            now=a[from[now]^1];
        }   
    }
    printf("%d",ans);
    return 0;
}