【TC SRM 575】TheTilesDivOne（最大流）

传送门

现在题面看不了，题意可以在TheTileDivTwo看，数据范围是把横纵坐标扩展到47

Div2传送门

---

题解：

注意到不管怎么旋转，两个白色格子一定不同行。

对于偶数行的白色格子，源点向它连容量为1的边，奇数行的白色格子，它向汇点连容量为1的边。

然后偶数行的白色格子向所有相邻的黑色格子连边，所有黑色格子向相邻的奇数行白色格子连边。

所有点拆开成入点和出点，限制流量为1。

跑最大流即为答案

---

代码：

#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const

using std::string;

class TheTilesDivOne{
	private:
		static cs int N=5e3+7,M=55;
		int S,T,n,m,tot;
		int in[M][M],out[M][M];
		
		struct edge{int to,cap,rev;};
		std::vector<edge> G[N]; 
		typedef std::vector<edge>::iterator iter;
		iter cur[N];int lev[N];
		inline void adde(int u,int v,int val){
			G[u].push_back((edge){v,val,G[v].size()});
			G[v].push_back((edge){u,0,G[u].size()-1});
		}
		
		inline bool BFS(){
			memset(lev,0,sizeof(int)*(tot+1));
			std::queue<int> q;q.push(S);lev[S]=1;
			cur[S]=G[S].begin();
			while(!q.empty()){
				int u=q.front();q.pop();
				for(auto e:G[u])if(e.cap&&!lev[e.to]){
					lev[e.to]=lev[u]+1;cur[e.to]=G[e.to].begin();
					if(e.to==T)return true;
					q.push(e.to);
				}
			}
			return false;
		}
		
		int Dinic(int u,int flow){
			if(u==T)return flow;
			int ans=0;
			for(iter &e=cur[u];e!=G[u].end();++e)
			if(e->cap&&lev[e->to]==lev[u]+1){
				int delta=Dinic(e->to,std::min(flow-ans,e->cap));
				if(delta){
					e->cap-=delta;
					G[e->to][e->rev].cap+=delta;
					if((ans+=delta)==flow)return flow;
				}
			}
			return ans;
		}
		
		inline int Flow(){
			int ans=0;
			while(BFS())ans+=Dinic(S,0x3f3f3f3f);
			return ans;
		}
		
		inline int col(int x,int y){
			if(~(x^y)&1)return 2;
			return (x&1)<<1|1;
		}
	public:
		TheTilesDivOne(){}
		int find(std::vector<string> s){
			n=s.size();m=s[0].size();
			S=tot=1,T=++tot;
			for(int re i=0;i<n;++i)
			for(int re j=0;j<m;++j)
			in[i][j]=++tot,out[i][j]=++tot;
			for(int re i=0;i<n;++i)
			for(int re j=0;j<m;++j)if(s[i][j]=='.'){
				adde(in[i][j],out[i][j],1);
				switch(col(i,j)){
					case 3:adde(out[i][j],T,1);break;
					case 1:adde(S,in[i][j],1);
					case 2:{
						static cs int dx[]={0,1,0,-1};
						static cs int dy[]={1,0,-1,0};
						for(int re d=0;d<4;++d){
							int x=i+dx[d],y=j+dy[d];
							if(0<=x&&x<n&&0<=y&&y<m&&
								s[x][y]=='.'&&col(x,y)==col(i,j)+1)
								adde(out[i][j],in[x][y],1);
						}
						break;
					}
				}
			}
			return Flow();
		}
};

#ifdef zxyoi 

TheTilesDivOne Solver;

signed main(){
	std::cout<<Solver.find(
		{"X.....XXX.XX..XXXXXXXXX...X.XX.XX....X",
 ".XXXX..X..XXXXXXXX....XX.X.X.X.....XXX",
 "....XX....X.XX..X.X...XX.X..XXXXXXX..X",
 "XX.XXXXX.X.X..X..XX.XXX..XX...XXX.X..."}
	)<<"\n";
	return 0;
}

#endif