高斯消元法求解开关问题

这类问题一般就是反转一个位置，那么和它相邻的其他位置也跟着反转，求全部反转成一样的最少反转次数。

在构造增广矩阵时，要将每个开关看成一个变元，并给每个开关列一个方程。即若有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;
}