題目連結:http://poj.org/problem?id=1222
題意:開關是四連通的,每按一個就會翻轉自己以及附近的四個格(假如有)。問需要翻轉幾個,使他們都變成關。
把每一個燈看作一個未知量,可以有30個未知量,給每一個未知量從0~29編號,再列30個方程,每一個方程中描述每一個燈被翻轉後其相鄰等的狀態變化,倒著考慮,把全關狀態翻成當前狀態就行了。
轉移矩陣是這樣的,很好理解。
#include <bits/stdc++.h>
using namespace std; const int maxn = ;
int equ, var;
int a[maxn][maxn];
int x[maxn];
int free_x[maxn];
int free_num; int gauss() {
int max_r, col, k;
free_num = ;
for(k = , col = ; k < equ && col < var; k++, col++) {
max_r = k;
for(int i = k + ; i < equ; i++) {
if(abs(a[i][col]) > abs(a[max_r][col]))
max_r = i;
}
if(a[max_r][col] == ) {
k--;
free_x[free_num++] = col;
continue;
}
if(max_r != k) {
for(int j = col; j < var + ; j++)
swap(a[k][j], a[max_r][j]);
}
for(int i = k + ; i < equ; i++) {
if(a[i][col] != ) {
for(int j = col; j < var + ; j++) {
a[i][j] ^= a[k][j];
}
}
}
}
for(int i = k; i < equ; i++) {
if(a[i][col] != )
return -;
}
if(k < var)
return var - k;
for(int i = var - ; i >= ; i--) {
x[i] = a[i][var];
for(int j = i + ; j < var; j++) {
x[i] ^= (a[i][j] & x[j]);
}
}
return ;
} int main() {
// freopen("in", "r", stdin);
equ = , var = ;
int T, _ = ;
scanf("%d", &T);
while(T--) {
memset(a, , sizeof(a));
memset(x, , sizeof(x));
memset(free_x, , sizeof(free_x));
for(int i = ; i < ; i++) {
scanf("%d", &a[i][var]);
}
for(int i = ; i < ; i++) {
a[i][i] = ;
if(i % != ) a[i-][i] = ;
if(i % != ) a[i+][i] = ;
if(i > ) a[i-][i] = ;
if(i < ) a[i+][i] = ;
}
int v = gauss();
int cnt = ;
printf("PUZZLE #%d\n", _++);
for(int i = ; i < ; i++) {
for(int j = ; j < ; j++) {
printf("%d ", x[cnt++]);
}
printf("\n");
}
}
return ;
}