设为首页 加入收藏

TOP

BZOJ 3503 CQOI 2014 和谐矩阵 高斯消元
2015-07-24 05:04:03 来源: 作者: 【 】 浏览:6
Tags:BZOJ 3503 CQOI 2014 和谐 矩阵 高斯

题目大意:给出m和n,求出一种方案使得每一个点和周围的四个点的1的个数为偶数。


思路:根据题意可以列出m*n个异或方程,然后组成异或方程组。解这个异或方程组然后输出任意一个解就可以了。

PS:值得注意的是,全是0肯定是一个解,显然你不能输出这个解。所以你需要让一个或一些自由元的值为1,至于怎么做,随便yy就行了。

PS2:这个题的样例吞掉了空格,然而又是SPJ,所以就是wa。。然后我wa了一下午。。


CODE:


#include 
  
   
#include 
   
     #include 
    
      #include 
     
       #define MAX 50 using namespace std; const int dx[] = {0,1,-1,0,0}; const int dy[] = {0,0,0,1,-1}; int num[MAX][MAX]; int arr[MAX * MAX][MAX * MAX],solves; bool ans[MAX]; int represent[MAX]; void Gauss() { int cnt = m * n,now = 1; for(int i = 1; i <= cnt; ++i) { for(int j = now + 1; j <= cnt; ++j) if(arr[j][i] > arr[now][i]) for(int k = 1; k <= cnt; ++k) swap(arr[now][k],arr[j][k]); if(!arr[now][i]) { ans[i] = 1; for(int j = 1; j <= cnt; ++j) if(arr[j][i]) arr[i][j] = 0,ans[i][cnt + 1] ^= 1; continue; } for(int j = now + 1; j <= cnt; ++j) if(arr[j][i]) for(int k = i; k <= cnt; ++k) arr[j][k] ^= arr[now][k]; represent[now] = i; ++now; } for(int i = now; i; --i) { int solve = represent[i]; for(int j = 1; j <= now - 1; ++j) if(arr[j][solve]) arr[j][solve] = 0,arr[j][cnt + 1] ^= arr[i][cnt + 1]; ans[solve] = arr[i][cnt + 1]; } } int main() { cin >> m >> n; for(int i = 1; i <= m; ++i) for(int j = 1; j <= n; ++j) num[i][j] = ++cnt; for(int i = 1; i <= m; ++i) for(int j = 1; j <= n; ++j) for(int k = 0; k <= 4; ++k) { int fx = i + dx[k],fy = j + dy[k]; if(!fx || !fy || fx > m || fy > n) continue; arr[num[i][j]][num[fx][fy]] = 1; } Gauss(); for(int i = 1; i <= m; ++i) { for(int j = 1; j <= n; ++j) printf("%d",ans[num[i][j]]); puts(""); } return 0; }
     
    
   
  


】【打印繁体】【投稿】【收藏】 【推荐】【举报】【评论】 【关闭】 【返回顶部
分享到: 
上一篇Codeforces Round #290 (Div. 2)(.. 下一篇C++拾遗--构造函数(二)

评论

帐  号: 密码: (新用户注册)
验 证 码:
表  情:
内  容: