CF #576 Div.1 E. Rectangle Painting 2 //最大流(最小点覆盖)

最新推荐文章于 2020-08-07 14:30:25 发布

RSHS

最新推荐文章于 2020-08-07 14:30:25 发布

阅读量237

点赞数

分类专栏：网络流

本文链接：https://blog.csdn.net/qq_41730604/article/details/100749074

版权

网络流专栏收录该内容

16 篇文章 0 订阅

订阅专栏

题意：
在这里插入图片描述

思路：
每个点只需要把宽或者长对应的染黑，那么就是最小点覆盖问题。
离散化，将所有边都延长割成网状，割成一个个小矩形，左闭右开统计大矩阵由小矩阵组成。

最小点覆盖 $\rightarrow$ 二分图最大匹配 $\rightarrow$ 最大流

建图:长 $\rightarrow$ 宽的边权为 $i n f$ ，其他边权为对应的长度

源点 $\rightarrow$ 离散化后的一段段x $\rightarrow$ 离散化后的一段段y $\rightarrow$ 汇点

跑最大流求最小覆盖

~~不愧*2500~~

/*   Author : Rshs
 *   Data : 2019-09-11-20.16
 */
#include<bits/stdc++.h>
using namespace std;
#define FI first
#define SE second
#define int long long
#define MP make_pair
#define PII pair<int,int>
#define SZ(a) (int)a.size()
const double pai = acos(-1);
const int mod = 1e9+7;
const int MX = 1e5+5;
//****************************************************
const int max_n = 500;
struct no{int to,cap,rev;}; //arc
vector<no>g[max_n]; //图
int level[max_n]; //到起点的距离
int iter[max_n]; //当前弧，在其之前的边已经没用了
void addarc(int s,int e,int c){
	g[s].push_back((no){e,c,(int)g[e].size()});
	g[e].push_back((no){s,0,(int)g[s].size()-1});
}
//更新层次，即level
void bfs(int s){
	memset(level,-1,sizeof(level));
	level[s]=0;
	queue<int>q;
	q.push(s);
	while(!q.empty()){
		int now=q.front();q.pop();
		for(int i=0;i<(int)g[now].size();i++){
			no &arc=g[now][i];
			if(level[arc.to]!=-1||arc.cap<=0) continue;
			level[arc.to]=level[now]+1;
			q.push(arc.to);
		}
	}
}
//寻找增广路
int dfs(int v,int t,int f){
	if(v==t) return f;
	for(iter[v];iter[v]<(int)g[v].size();iter[v]++){
		no &arc=g[v][iter[v]];
		if(arc.cap<=0||level[arc.to]!=level[v]+1) continue;
		int d=dfs(arc.to,t,min(f,arc.cap));
		if(d>0) {
			arc.cap=arc.cap-d;
			g[arc.to][arc.rev].cap+=d;
			return d;
		}
	}
	return 0;
}
int Dinic(int s,int t){
	int re=0;
	while(1){
		bfs(s);
		memset(iter,0,sizeof(iter));
		if(level[t]==-1) return re;
		int f;
		while((f=dfs(s,t,LLONG_MAX))>0)
			re=re+f;
	}
	return re;
}
//****************************************************
int a[MX],b[MX],x[MX],y[MX];
vector<int>lx,ly;
int vis[200][200];
int main(){
    int n,m;cin>>n>>m;
    for(int i=1;i<=m;i++){
        cin>>a[i]>>b[i]>>x[i]>>y[i];
        x[i]++,y[i]++;
        lx.push_back(a[i]);lx.push_back(x[i]);
        ly.push_back(b[i]),ly.push_back(y[i]);
    }
    sort(lx.begin(),lx.end());sort(ly.begin(),ly.end());
    lx.erase(unique(lx.begin(),lx.end()),lx.end());
    ly.erase(unique(ly.begin(),ly.end()),ly.end());
    int sx=SZ(lx),sy=SZ(ly);int SS=0,TT=sx+sy+1;
    for(int i=1;i<sx;i++) addarc(SS,i,lx[i]-lx[i-1]); //左闭又开，S连边
    for(int i=1;i<sy;i++) addarc(i+sx,TT,ly[i]-ly[i-1]);
    for(int i=1;i<=m;i++){
        int l1,l2,r1,r2;
        l1=lower_bound(lx.begin(),lx.end(),a[i])-lx.begin()+1;
        l2=lower_bound(lx.begin(),lx.end(),x[i])-lx.begin()+1;
        r1=lower_bound(ly.begin(),ly.end(),b[i])-ly.begin()+1;
        r2=lower_bound(ly.begin(),ly.end(),y[i])-ly.begin()+1;
        for(int j=l1;j<l2;j++){
            for(int k=r1;k<r2;k++){
                vis[j][k]=1;
            }
        }
    }
    for(int j=1;j<sx;j++){
        for(int k=1;k<sy;k++){
            if(vis[j][k]) addarc(j,k+sx,LLONG_MAX/10000);
        }
    }
    cout<<Dinic(SS,TT);
    return 0;
}