码迷,mamicode.com
首页 > 其他好文 > 详细

BZOJ P1059 [ZJOI2007]矩阵游戏——solution

时间:2017-05-25 20:48:06      阅读:167      评论:0      收藏:0      [点我收藏+]

标签:page   names   esc   struct   com   pac   define   stat   using   

1059: [ZJOI2007]矩阵游戏

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 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

Sample Output

No
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

标签:page   names   esc   struct   com   pac   define   stat   using   

原文地址:http://www.cnblogs.com/nietzsche-oier/p/6905545.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!