这类问题一般就是反转一个位置,那么和它相邻的其他位置也跟着反转,求全部反转成一样的最少反转次数。
在构造增广矩阵时,要将每个开关看成一个变元,并给每个开关列一个方程。即若有n个开关,就要列n个方程,每个方程含有n个变元。若开关i需要反转,则将第i个方程的解设为1(a[i][var]=1),否则为0(a[i][var]=0)。若反转开关j会使开关i同时反转,则在系数矩阵中a[i][j] = 1.
构造完增广矩阵后,用高斯消元求解。
若方程组无解,则返回无解。
若有唯一解,带回求出解,x[i](i=0~var-1)中1的个数就是所需的反转次数
若有多解,枚举解(解的个数为1<<(自由变元个数)),找出其中反转次数最小的。
poj--1753 Flip Game
给定4乘4的01矩阵,反转成全0或者全1。求最小的反转次数。
http://blog.csdn.net/fad_further/article/details/23107597
poj -- 3185 The Water Bowls
有一个长度为20的01序列,要将其翻成全0,每次翻动也会将其两侧的字符反转。求最小反转次数。
http://blog.csdn.net/fad_further/article/details/24549669
poj--1681 Painter's Problem
给一个n*n的矩阵将其全部反转成‘y’,求最少的反转次数
#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <queue>
#include <stack>
#include <cstdio>
#include <vector>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define lson l, m, rt << 1
#define rson m + 1, r, rt << 1 | 1
#define pb push_back
#define MP make_pair
typedef long long ll;
typedef pair<int, int> PII;
const int inf = 0x3f3f3f3f;
const int mod = 1000000007;
const int maxn = 400 + 10;
int dxy[4][2] = { {1, 0}, {0, 1}, {-1, 0}, {0, -1} };
int n, t;
char str[20];
int equ, var, a[maxn][maxn], x[maxn], free_num;
bool free_x[maxn];
void init(){
memset(a, 0, sizeof(a));
memset(x, 0, sizeof(x));
memset(free_x, true, sizeof(free_x));
equ = var = n * n;
}
void Debug(){
puts("debug");
for(int i=0; i<equ; i++)
for(int j=0; j<=var; j++)
printf("%d%c", a[i][j], j == var ? '\n' : ' ');
puts("end");
}
void Guass(){
int col, row, mx;
// Debug();
for(col=0, row=0; col<var&&row<equ; col++, row++){
mx = row;
for(int i=row+1; i<equ; i++)
if(abs(a[i][col]) > abs(a[mx][col])) mx = i;
if(a[mx][col] == 0){
row--;
continue;
}
if(mx != row)
for(int i=col; i<=var; i++)
swap(a[mx][i], a[row][i]);
for(int i=row+1; i<equ; i++){
if(a[i][col] == 0) continue;
for(int j=col; j<=var; j++)
a[i][j] ^= a[row][j];
}
}
// Debug();
for(int i=row; i<equ; i++)
if(a[i][var] != 0){
puts("inf");
return;
}
if(row == equ){
for(int i=var-1; i>=0; i--){
x[i] = a[i][var];
for(int j=i+1; j<var; j++)
if(a[i][j] && x[j])
x[i] ^= 1;
}
int ans = 0;
for(int i=0; i<var; i++)
ans += x[i];
printf("%d\n", ans);
return;
}
for(int i=row-1; i>=0; i--){
int idx;
free_num = 0;
for(int j=0; j<var; j++)
if(a[i][j] && free_x[j]) {
free_num++;
idx = j;
}
if(free_num > 1) continue;
free_x[idx] = false;
}
int num = var - row, xx[maxn], cnt = 0, ans = inf;
for(int i=var-1; i>=0 && cnt<num; i--)
if(free_x[i]) xx[cnt++] = i;
for(int i=0; i<(1 << num); i++){
for(int j=0; j<num; j++)
if(i & (1 << j))x[xx[j]] = 1;
else x[xx[j]] = 0;
for(int j=row-1; j>=0; j--){
int tmp;
for(tmp=0; tmp<var; tmp++)
if(a[j][tmp]) break;
x[tmp] = a[j][var];
for(int k=tmp+1; k<=var; k++)
if(a[j][k] && x[k])
x[tmp] ^= 1;
}
int sum = 0;
for(int j=0; j<var; j++) sum += x[j];
ans = min(sum, ans);
}
printf("%d\n", ans);
}
int main(){
scanf("%d", &t);
while(t--){
scanf("%d", &n);
init();
for(int i=0; i<n; i++){
scanf("%s", str);
for(int j=0; str[j]; j++){
int tmp = i * n + j;
if(str[j] == 'w') a[tmp][var] = 1;
for(int k=0; k<4; k++){
int xx = i + dxy[k][0];
int yy = j + dxy[k][1];
if(xx < 0 || xx >= n || yy < 0 || yy >= n) continue;
a[tmp][xx * n + yy] = 1;
}
}
}
for(int i=0; i<equ; i++) a[i][i] = 1;
Guass();
}
return 0;
}
poj--1830 开关问题
给n个开关,并给出他们之间的相互联系求有多少种到达指定状态的方法。
方法数就是1 << 自由变元的个数
#include <queue>
#include <stack>
#include <cstdio>
#include <vector>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
const int maxn = 40;
const int inf = 0x3f3f3f3f;
int equ, var;
int a[maxn][maxn];
int x[maxn];
bool free_x[maxn];
int free_num;
int T, n;
void Debug(){
puts("debug:");
for(int i=0; i<equ; i++)
for(int j=0; j<=var; j++)
printf("%d%c", a[i][j], j == var ? '\n' : ' ');
puts("end");
}
void init(){
memset(a, 0, sizeof(a));
memset(x, 0, sizeof(x));
memset(free_x, true, sizeof(free_x));
}
void Guass(){
int col = 0, row = 0, mx;
// Debug();
for(; col<var&&row<equ; row++, col++){
mx = row;
for(int i=row+1; i<equ; i++)
if(abs(a[i][col]) > abs(a[mx][col])) mx = i;
if(a[mx][col] == 0){
row--;
continue;
}
if(mx != row){
for(int i=col; i<=var; i++)
swap(a[mx][i], a[row][i]);
}
for(int i=row+1; i<equ; i++){
if(a[i][col] == 0) continue;
for(int j=col; j<=var; j++)
a[i][j] ^= a[row][j];
}
}
for(int i=row; i<equ; i++)
if(a[i][var] != 0) {
puts("Oh,it's impossible~!!");
return;
}
printf("%d\n", (1 << (var - row)));
}
int main(){
scanf("%d", &T);
while(T--){
init();
scanf("%d", &n);
int x, y, s[maxn], t[maxn];
for(int i=0; i<n; i++){
scanf("%d", &s[i]);
}
for(int i=0; i<n; i++){
scanf("%d", &t[i]);
}
while(scanf("%d%d", &x, &y) == 2){
if(!x && !y) break;
a[y - 1][x - 1] = 1;
}
equ = var = n;
for(int i=0; i<n; i++){
a[i][var] = (s[i] ^ t[i]);
a[i][i] = 1;
}
Guass();
}
return 0;
}