BZOJ P1059 [ZJOI2007]矩阵游戏——solution
Posted Der Barde, Nietzsche
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了BZOJ P1059 [ZJOI2007]矩阵游戏——solution相关的知识,希望对你有一定的参考价值。
1059: [ZJOI2007]矩阵游戏
Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 4604 Solved: 2211
[Submit][Status][Discuss]
Description
小Q是一个非常聪明的孩子,除了国际象棋,他还很喜欢玩一个电脑益智游戏——矩阵游戏。矩阵游戏在一个N
*N黑白方阵进行(如同国际象棋一般,只是颜色是随意的)。每次可以对该矩阵进行两种操作:行交换操作:选择
矩阵的任意两行,交换这两行(即交换对应格子的颜色)列交换操作:选择矩阵的任意行列,交换这两列(即交换
对应格子的颜色)游戏的目标,即通过若干次操作,使得方阵的主对角线(左上角到右下角的连线)上的格子均为黑
色。对于某些关卡,小Q百思不得其解,以致他开始怀疑这些关卡是不是根本就是无解的!!于是小Q决定写一个程
序来判断这些关卡是否有解。
Input
第一行包含一个整数T,表示数据的组数。接下来包含T组数据,每组数据第一行为一个整数N,表示方阵的大
小;接下来N行为一个N*N的01矩阵(0表示白色,1表示黑色)。
Output
输出文件应包含T行。对于每一组数据,如果该关卡有解,输出一行Yes;否则输出一行No。
Sample Input
2
2
0 0
0 1
3
0 0 1
0 1 0
1 0 0
2
0 0
0 1
3
0 0 1
0 1 0
1 0 0
Sample Output
No
Yes
【数据规模】
对于100%的数据,N ≤ 200
Yes
【数据规模】
对于100%的数据,N ≤ 200
通过对行列的移动使黑点移动到正对角线上,
对于一种可行的步骤,至少先移动行与先移动列是一样的;
于是:
如果存在方案,则存在先完成所有列移动再移动行的方案;
讨论这一类方案:
在列操作之后,每行中有黑点的列已经确定了,则该行可以通过移动而贡献到对角线上哪些点,也已经确定了;
但每行只能贡献对角线上的一个点;(注意)
于是有了模糊的二分图匹配的思路;
然后讨论在行操作之前的列操作:
她其实没什么用;
因为她大概可以改变对角线上某一点永远不能被行移动贡献的情况;
但这只是一种拆东墙补西墙的行为,没有意义;
于是直接忽略列操作;
二分图匹配;
在行号和用列好表示的对角线上点间建边
对行和每行的黑点的列编号建边即可;
代码:
1 #include<cstdio> 2 #include<cstring> 3 #define INF 2147483000 4 using namespace std; 5 int n,m,S,T; 6 struct ss{ 7 int next,to,w,cp; 8 }x[5000010]; 9 int first[4010],num; 10 int dep[4010]; 11 int que[1000000]; 12 void bui_(int ,int ,int ); 13 void build(int ,int ,int ); 14 int bfs(); 15 int dfs(int ,int ); 16 int main() 17 { 18 int TT,i,j,k,ans,add; 19 scanf("%d",&TT); 20 while(TT--){ 21 memset(first,0,sizeof(first));num=0;ans=0; 22 scanf("%d",&n); 23 S=0,T=(n<<1)+1; 24 for(i=1;i<=n;i++){ 25 bui_(S,i,1);bui_(n+i,T,1); 26 for(j=1;j<=n;j++){ 27 scanf("%d",&k); 28 if(k) 29 bui_(i,n+j,1); 30 } 31 } 32 while(bfs()) 33 while(add=dfs(S,INF)) 34 ans+=add; 35 if(ans==n)printf("Yes\n"); 36 else printf("No\n"); 37 } 38 return 0; 39 } 40 void bui_(int f,int t,int d){ 41 build(f,t,d);x[num].cp=num+1; 42 build(t,f,0);x[num].cp=num-1; 43 } 44 void build(int f,int t,int d){ 45 x[++num].next=first[f]; 46 x[num].to=t; 47 x[num].w=d; 48 first[f]=num; 49 } 50 int bfs(){ 51 memset(dep,0,sizeof(dep)); 52 int h=0,t=1,j; 53 que[t]=S; 54 while(h<t){ 55 ++h; 56 for(j=first[que[h]];j;j=x[j].next) 57 if(x[j].w&&x[j].to!=S&&!dep[x[j].to]){ 58 dep[x[j].to]=dep[que[h]]+1; 59 que[++t]=x[j].to; 60 } 61 } 62 if(dep[T])return 1; 63 return 0; 64 } 65 int dfs(int now,int min){ 66 int j,re=0; 67 if(now==T) 68 return min; 69 for(j=first[now];j;j=x[j].next) 70 { 71 if(x[j].w&&dep[x[j].to]==dep[now]+1){ 72 re=dfs(x[j].to,min<x[j].w?min:x[j].w); 73 if(re){ 74 x[j].w-=re; 75 x[x[j].cp].w+=re; 76 return re; 77 } 78 dep[x[j].to]=0; 79 } 80 } 81 return re; 82 }
以上是关于BZOJ P1059 [ZJOI2007]矩阵游戏——solution的主要内容,如果未能解决你的问题,请参考以下文章
bzoj 1059: [ZJOI2007]矩阵游戏 二分图匹配