(网络流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;
}