标签：head 匹配 int Day2 Alice ecnt 雅礼 len 2017
description
给一个\(n*m\)的棋盘，'.'为可通行，'#'为障碍。Alice选择一个起始点，Bob先手从该点往四个方向走一步，Alice再走，不能走走过的点，谁不能动谁就输了。
问Alice选择哪些出发点能赢。
solution
棋盘黑白染色
分类后建二分图
1.Alice选择关键匹配点(所有的最大匹配都包含该点)，Bob赢。
2.否则，Alice选择(含不包含该点的匹配），Alice赢。
证明：如果为2的话（假如Alice选的起点为S），存在一个最大匹配是不包含S的，所以对于该匹配Bob所到的一定是匹配点，接着Alice一直走非匹配点，Bob只能走到匹配点……
最后一定结束在匹配点，否则会有新的增广路。
做法肯定能想到暴力删每个点跑一次，判断是否为关键匹配点。
不过复杂度太爆炸，因此先跑一次最大流得到匹配边。其中没有被选的点已经确定为非关键匹配点了，它所连出的边（非匹配边）所到的另一个集合的点再回到它所在的集合中的匹配点，那这个匹配点一定是非关键匹配点（因为可以由它替换）。这个利用dfS配残量len=1就可以判断，记得从\(t\)搜的时候要判反边的len。
code
戳我
#include<bits/stdc++.h>
using namespace std;
const int N=1e3+5;
const int M=1e6+5;
const int inf=1e9;
char mp[N][N];
int n,m,dir[5][2]={{0,-1},{0,1},{1,0},{-1,0}};
int s,t,dis[M],gap[M],id[N][N],nxt[M],to[M],head[M],len[M],ecnt=1;
bool vis[M];
void add_edge(int u,int v,int z) {
nxt[++ecnt]=head[u];to[ecnt]=v;len[ecnt]=z;head[u]=ecnt;
nxt[++ecnt]=head[v];to[ecnt]=u;len[ecnt]=0;head[v]=ecnt;
}
queue<int> Q;
void BFS() {
for(int i=s;i<=t;i++) dis[i]=-1;
Q.push(t);gap[dis[t]=0]++;
while(!Q.empty()) {
int u=Q.front(); Q.pop();
for(int i=head[u];i;i=nxt[i]) {
int v=to[i];
if(dis[v]!=-1)continue;
gap[dis[v]=dis[u]+1]++;Q.push(v);
}
}
}
int dfs(int u,int flow) {
if(u==t) return flow;
int used=0;
for(int i=head[u];i;i=nxt[i]) {
int v=to[i];
if(len[i]&&dis[u]==dis[v]+1) {
int tmp=dfs(v,min(len[i],flow-used));
if(tmp) {len[i]-=tmp;len[i^1]+=tmp;used+=tmp;}
if(used==flow)return used;
}
}
--gap[dis[u]];
if(!gap[dis[u]]) dis[s]=t;
dis[u]++,gap[dis[u]]++;
return used;
}
int mxflow=0;
void ISAP() {
for(BFS();dis[s]<t;mxflow+=dfs(s,inf));
}
void Build() {
int tc=0;
s=0;t=n*m+1;
for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)id[i][j]=++tc;
for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) {
if(((i+j)&1)&&mp[i][j]=='.') {
for(int d=0;d<4;d++) {
int x=i+dir[d][0],y=j+dir[d][1];
if(id[x][y]&&mp[x][y]=='.')add_edge(id[i][j],id[x][y],1);
}
add_edge(s,id[i][j],1);
}
else add_edge(id[i][j],t,1);
}
}
void dfs1(int u) {
vis[u]=1;
for(int i=head[u];i;i=nxt[i]) if(!vis[to[i]]&&len[i]) dfs1(to[i]);
}
void dfs2(int u) {
vis[u]=1;
for(int i=head[u];i;i=nxt[i]) if(!vis[to[i]]&&len[i^1]) dfs2(to[i]);
}
bool ans[N][N];
void solve() {
ISAP();dfs1(s);
int cnt=0;
for(int i=1;i<=n;i++)for(int j=1;j<=m;j++) {
if(((i+j)&1)&&vis[id[i][j]]&&mp[i][j]=='.') {ans[i][j]=1;cnt++;}
vis[id[i][j]]=0;
}
dfs2(t);
for(int i=1;i<=n;i++)for(int j=1;j<=m;j++) {
if((i+j)&1)continue;
if(mp[i][j]=='.'&&vis[id[i][j]]) {ans[i][j]=1;cnt++;}
}
printf("%d\n",cnt);
for(int i=1;i<=n;i++)for(int j=1;j<=m;j++) {
if(ans[i][j]) printf("%d %d\n",i,j);
}
}
int main() {
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++) scanf("%s",mp[i]+1);
Build();
solve();
return 0;
}
标签：head,匹配,int,Day2,Alice,ecnt,雅礼,len,2017 来源： https://www.cnblogs.com/bestime/p/16107352.html
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。