BZOJ 3503 CQOI 2014 和谐矩阵 高斯消元

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

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

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

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

CODE:

#include <cstdio>#include <cstring>#include <iostream>#include <algorithm>#define MAX 50using 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;}

当一个人真正觉悟的一刻,他放弃追寻外在世界的财富,而开始追寻他内心世界的真正财富

BZOJ 3503 CQOI 2014 和谐矩阵 高斯消元

相关文章:

你感兴趣的文章:

标签云: