给定实直线 L L L 上 n n n 个开区间组成的集合 I I I,和一个正整数 k k k,试设计一个算法,从开区间集合 I I I 中选取出开区间集合 S⊆I S \subseteq I S⊆I,使得在实直线 L L L 的任何一点 x x x,S S S 中包含点 x x x 的开区间个数不超过 k k k。且 ∑z∈S∣z∣ \sum\limits_{z \in S} | z | z∈S∑∣z∣ 达到最大。这样的集合 S S S 称为开区间集合 I I I 的最长 k k k 可重区间集。∑z∈S∣z∣ \sum\limits_{z \in S} | z | z∈S∑∣z∣ 称为最长 k k k 可重区间集的长度。
对于给定的开区间集合 I I I 和正整数 k k k,计算开区间集合 I I I 的最长 k k k 可重区间集的长度。
输入格式
文件的第 1 1 1 行有 2 2 2 个正整数 n n n 和 k k k,分别表示开区间的个数和开区间的可重迭数。
接下来的 n n n 行,每行有 2 2 2 个整数 li l_i li 和 ri r_i ri,表示开区间的左右端点坐标,注意可能有 li>ri l_i > r_i li>ri,此时请将其交换。
输出格式
输出最长 k k k 可重区间集的长度。
样例
样例输入
4 2
1 7
6 8
7 10
9 13
样例输出
15
数据范围与提示
1≤n≤500,1≤k≤3 1 \leq n \leq 500, 1 \leq k \leq 3 1≤n≤500,1≤k≤3
构建一个图 保证满流是满足条件的
首先 源向起点连边权为k的权值 然后 最末端向汇点连边权为Inf 费用为0的边
每个区间从左端点向右端点连边 权值为1 表示这个区间只能使用它覆盖一次 费用为这个区间的长度 然后跑最大费用最大流即可求出答案
#include<queue>
#include<vector>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define inf 0x3f3f3f3f
#define N 1100
using namespace std;
vector<int> aq;
inline char gc(){
static char now[1<<16],*S,*T;
if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
return *S++;
}
inline int read(){
int x=0;char ch=gc();
while(ch<'0'||ch>'9') ch=gc();
while(ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=gc();}
return x;
}
int h[N],pre[N],flag[N],f[N],path[N],l[N],r[N],n,k,num=1,T;
struct node{
int y,z,next,c;
}data[N*N];
inline void insert1(int x,int y,int z,int c){
data[++num].y=y;data[num].z=z;data[num].next=h[x];h[x]=num;data[num].c=c;
data[++num].y=x;data[num].z=0;data[num].next=h[y];h[y]=num;data[num].c=-c;
}
inline bool spfa(){
memset(flag,0,sizeof(flag));memset(f,128,sizeof(f));queue<int>q;q.push(0);flag[0]=1;f[0]=0;memset(pre,-1,sizeof(pre));
while(!q.empty()){
int x=q.front();q.pop();flag[x]=0;
for (int i=h[x];i;i=data[i].next){
int y=data[i].y,z=data[i].z,c=data[i].c;
if (f[x]+c>f[y]&&z){
f[y]=f[x]+c;pre[y]=x;path[y]=i;
if (!flag[y]) flag[y]=1,q.push(y);
}
}
}if (pre[T]==-1) return 0;else return 1;
}
int main(){
freopen("loj.in","r",stdin);
n=read();k=read();
for (int i=1;i<=n;++i){
l[i]=read();r[i]=read();if(l[i]>r[i]) swap(l[i],r[i]);aq.push_back(l[i]);aq.push_back(r[i]);
}sort(aq.begin(),aq.end());insert1(0,1,k,0);int nn=unique(aq.begin(),aq.end())-aq.begin();T=nn+1;
for (int i=1;i<=n;++i){
int cost=r[i]-l[i];
int ll=lower_bound(aq.begin(),aq.begin()+nn,l[i])-aq.begin()+1,rr=lower_bound(aq.begin(),aq.begin()+nn,r[i])-aq.begin()+1;
insert1(ll,rr,1,cost);
}insert1(nn,T,inf,0);
for (int i=1;i<nn;++i) insert1(i,i+1,inf,0);int ans=0;
while(spfa()){
int minn=inf,now=T;
while(now) minn=min(minn,data[path[now]].z),now=pre[now];now=T;
while(now){ans+=data[path[now]].c*minn;data[path[now]].z-=minn;data[path[now]^1].z+=minn;now=pre[now];}
}printf("%d",ans);
return 0;
}