【TC SRM 594】FoxAndGo3 （二分图最大独立集）

传送门

---

题解：

注意到一个事实，白棋与和它相邻的一个空格子在终止情况下不可能同时为空。

我们连边表示两个位置不可能同时为空，显然这是一个二分图。

为了最大化空的位置，我们需要在这里面找一个最大独立集。

直接建图然后Dinic求一个二分图最大匹配就行了。

---

代码：

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

using std::string;

class FoxAndGo3{
	private:
		static cs int N=5e3+7,M=55;
		int S,T,tot,n,m,id[M][M];
		struct edge{int to,cap,rev;};
		typedef std::vector<edge>::iterator iter;
		std::vector<edge> G[N];iter cur[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});
		}
		
		int lev[N];
		inline bool BFS(){
			memset(lev+1,0,sizeof(int)*tot);
			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;if(e.to==T)return true;
					cur[e.to]=G[e.to].begin();q.push(e.to);
				}
			}
			return false;
		}
		
		int dfs(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=dfs(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 ans;
				}
			}
			return ans;
		}
		
		inline int Flow(){
			int flow=0;
			while(BFS())flow+=dfs(S,0x3f3f3f3f);
			return flow;
		}
	public:
		FoxAndGo3(){}
		int maxEmptyCells(std::vector<string> s){
			n=s.size();m=s[0].size();
			S=1,T=tot=2;int ans=0;
			for(int re i=0;i<n;++i)
			for(int re j=0;j<m;++j)
			if(s[i][j]!='x')id[i][j]=++tot,++ans;
			for(int re i=0;i<n;++i)
			for(int re j=0;j<m;++j)if(id[i][j]){
				if(s[i][j]=='o'){
					adde(S,id[i][j],1);
					for(int re k=0;k<4;++k){
						static cs int dx[]={0,1,0,-1};
						static cs int dy[]={1,0,-1,0};
						int x=i+dx[k],y=j+dy[k];
						if(!~x||x==n||!~y||y==m||s[x][y]!='.')continue;
						adde(id[i][j],id[x][y],1);
					}
				}
				else adde(id[i][j],T,1);
			}
			return ans-Flow();
		}
};

#ifdef zxyoi

FoxAndGo3 Solver;

signed main(){
	std::cout<<Solver.maxEmptyCells(
{"o.o.o",
 ".ox..",
 "oxxxo",
 "..x..",
 "o.o.o"}
	)<<"\n";
	
	return 0;
}

#endif