https://www.luogu.com.cn/problem/solution/P3358
洛谷大佬们题解写的很棒,我就不献丑了。
说下我的理解吧,可以把这个限制过程当初物理中的电路并联与串联。流量就是电流。
我们让源点产生k电流。这样保证流经每个点的电流不超过k。
每个区间我们给他变成左开右闭区间。
仔细体会下这个过程,很有意思的!
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int MAXN = 5001;
const int MAXM = 50001;
struct MCMF
{
int n, s, t, cnt = 1;
long long maxflow=0, mincost=0;
int dis[MAXN], head[MAXN], incf[MAXN], pre[MAXN];//dis表示最短路,incf表示当前增广路上最小流量,pre表示前驱
bool vis[MAXN];
struct EDGE {
int nxt, to, cost, flow;
}ee[MAXM << 1];
inline void AD(int from, int to, int flow, int cost) {
ee[++cnt].nxt = head[from];
ee[cnt].to = to;
ee[cnt].cost = cost;
ee[cnt].flow = flow;
head[from] = cnt;
}
inline void add(int x, int y, int flow, int cost){
AD(x,y,flow,cost);AD(y,x,0,-cost);
//双向边
}
inline bool spfa()
{
queue <int> q;
for(int i=1;i<=n;i++)dis[i]=0x3f3f3f3f,vis[i]=0;
q.push(s);
dis[s] = 0;vis[s] = 1;
incf[s] = 1 << 30;
while(!q.empty()) {
int u = q.front();
vis[u] = 0;
q.pop();
for(register int i = head[u]; i; i = ee[i].nxt) {
if(!ee[i].flow) continue;//没有剩余流量
int v = ee[i].to;
if(dis[v] > dis[u] + ee[i].cost) {
dis[v] = dis[u] + ee[i].cost;
incf[v] = min(incf[u], ee[i].flow);//更新incf
pre[v] = i;
if(!vis[v]) vis[v] = 1, q.push(v);
}
}
}
if(dis[t] == 0x3f3f3f3f) return 0;
return 1;
}
void gao()
{
while(spfa()) {//如果有增广路
int x = t;
maxflow += incf[t];
mincost += (long long)dis[t] * incf[t];
int i;
while(x != s) {//遍历这条增广路,正向边减流反向边加流
i = pre[x];
ee[i].flow -= incf[t];
ee[i^1].flow += incf[t];
x = ee[i^1].to;
}
}
}
//先初始化
void init(int nn,int S,int T)
{
n=nn,s=S,t=T;
cnt=1;
for(int i=0;i<=n;i++)head[i]=0;
mincost=maxflow=0;
}
}mc;
const int M =1e5+7;
int li[M],sz;
struct node{
int l,r,len;
}p[M];
int main()
{
ios::sync_with_stdio(false);
cin.tie(0);
int n,k;
cin>>n>>k;
for(int i=1;i<=n;i++){
int l,r;
cin>>l>>r;
p[i]=node{l,r-1,r-l};//左开右闭区间
li[++sz]=l,li[++sz]=r-1;
}
sort(li+1,li+1+sz);
sz=unique(li+1,li+1+sz)-(li+1);
int s=sz+1,t=sz+2;
mc.init(sz+2,s,t);
int inf=0x3f3f3f3f;
mc.add(s,1,k,0);//首末流量设为k,限制了流向所有点的流量不超过k即选某个点不超过k次
mc.add(sz,t,k,0);
for(int i=1;i<sz;i++)
mc.add(i,i+1,inf,0);
for(int i=1;i<=n;i++){
p[i].l=lower_bound(li+1,li+1+sz,p[i].l)-li;
p[i].r=lower_bound(li+1,li+1+sz,p[i].r)-li;
mc.add(p[i].l,p[i].r,1,-p[i].len);//建负边,跑最大费用流
}
mc.gao();
cout<<-mc.mincost<<endl;
return 0;
}