(完结撒花)网络流二十三题 ————(二十三）、P3357 最长k可重线段集问题 费用流并联串联问题

先说下为啥完结撒花。。

刷完感觉网络流无敌了！！！

https://blog.csdn.net/bjfu170203101/article/details/109062331

这道题的升级版，懒得再说一遍了。

不过这题有几个要注意的点：

1：坐标距离会爆longlong；

2：1 1 2 2 这种  (x1+1==x2)的线段对结果没有影响，即：这些线段可以随便取，先把这些线段取完扔掉。

3：1 1 1 2 这种与x轴垂直的线段是与x=1相交的，要注意处理这种情况：

比较好的处理方法是：坐标扩域。

比如x轴有1-5，则扩为:

1122334455.

正常x1+1<x2的线段 连边： x1*2 -> x2*2(左开右闭)  你会发现刚好线段影响的数是x1+1 -> x2范围内的 符合左开右闭区间的性质

若x1==x2，则直接连：x1*2-1    ->   x2*2  这样刚好只影响x1这个点。

 

#include<bits/stdc++.h> 
using namespace std;
typedef long long ll;
const int N=1000005;
const int M=1000005;
const ll inf=0x7fffffff;
struct MCMF{
	struct EDGE {
		ll nxt, to, cost, flow;
	}ee[M << 1];
	ll n,s,t;
	ll cur[N],head[N],vis[N];
	ll dis[N];
	ll mincost=0,maxflow=0;
	inline void AD(int from, int to, ll flow, ll cost) {
		ee[++cnt].nxt = head[from];
		ee[cnt].to = to;
		ee[cnt].cost = cost;
		ee[cnt].flow = flow;
		head[from] = cnt;
	}
	int cnt=1;
	inline void add(int u,int v,ll flow,ll cost){
	    AD(u,v,flow,cost);
	    AD(v,u,0,-cost);
//	   cout<<u<<"  ->  "<<v<<endl;
	}
	inline int spfa()
	{
	    for(int i=1;i<=n;++i)dis[i]=inf,vis[i]=0;
	    queue<ll>q;q.push(s);
	    dis[s]=0;vis[s]=1;
	    while(!q.empty()){
	        int u=q.front();
	        vis[u]=0;q.pop();
	        for(int i=head[u];i;i=ee[i].nxt)
	        {
	            int v=ee[i].to;
	            if(!ee[i].flow)continue;
	            if(ee[i].flow&&dis[v]>dis[u]+ee[i].cost){
	                dis[v]=dis[u]+ee[i].cost;
	                if(!vis[v]){
	                    vis[v]=1;
	                    q.push(v);
	                }
	            }
	        }
	    }
	    return dis[t]!=inf;
	}
	inline ll dfs(int u,ll flow)
	{
	    if(u==t)return flow;
	    ll rest=flow;
	    vis[u]=1;
	    for(int i=cur[u];i&&rest;i=ee[i].nxt){
	        cur[u]=i;
	        int v=ee[i].to;
	        ll w=ee[i].cost,flow=ee[i].flow;
	        if(ee[i].flow&&(dis[v]==dis[u]+w)&&!vis[v]){
	            ll k=dfs(v,min(flow,rest));
	            if(k){
	                mincost+=k*w;
	                ee[i].flow-=k;
	                ee[i^1].flow+=k;
	                rest-=k;
	            }
	        }
	    }
	    vis[u]=0;
	    return flow-rest;
	}
	void gao(){
		while(spfa()){
	        for(int i=1;i<=n;++i)cur[i]=head[i];
	        maxflow+=dfs(s,inf);
	    }
	}
	void init(int nn,int ss,int tt){
		n=nn,s=ss,t=tt;
		maxflow=mincost=0;
		for(int i=0;i<=n;i++)head[i]=0;
	}
}mc;
int s,t,nm,n,m,k,x,y;
void path(int ip,int p,int t){//打印路径 
	while(p!=t){
	//	cout<<p<<" -> ";
		for(int i=mc.head[p];i;i=mc.ee[i].nxt){
			if(i&1)continue;//cout<<p<<"  -   "<<mc.ee[i].to<<"     "<<mc.ee[i].flow<<endl;
			if(mc.ee[i^1].flow){
				mc.ee[i^1].flow--;
				int to=mc.ee[i].to;
			//	cout<<p<<" -> "<<to<<"  =  ";
				if(to>n*m){
					if(to-n*m-p==1)cout<<ip<<" "<<1<<endl;
					else cout<<ip<<" "<<0<<endl;
				}
				p=to;
				break;
			}
		}
	}
}

ll li[7777];
struct node{
	ll l,r,len,op;
}p[7777];
ll dis(ll x1,ll y1,ll x2,ll y2){
	return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
int main()
{
	ios::sync_with_stdio(false);
  	cin.tie(0);
  	int n,k,sz=0;
  	cin>>n>>k;	
	ll ans=0;
  	for(int i=1;i<=n;i++){
  		ll x1,x2,y1,y2;
  		cin>>x1>>y1>>x2>>y2;
  		ll ds=dis(x1,y1,x2,y2);
  		if(x1>x2)swap(x1,x2),swap(y1,y2);
  		if(x1+1==x2){
  			ans+=ds;
  			continue;
		}
		if(x1==x2){
			li[++sz]=x1*2-1,li[++sz]=x2*2;
			p[i]=node{x1*2-1,x2*2,ds,1};
			continue;
		}
  		li[++sz]=x1*2,li[++sz]=x2*2;
		p[i]=node{x1*2,x2*2,ds,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);
	for(int i=1;i<sz;i++)
		mc.add(i,i+1,inf,0);
//	cout<<"----- "<<endl;
	for(int i=1;i<=n;i++){
		if(p[i].op==0)continue;
		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);
	}
//	cout<<"------------------"<<endl;
	mc.add(s,1,k,0);
	mc.add(sz,t,inf,0);
	mc.gao();
	cout<<ans-mc.mincost<<endl;
	return 0;
}