Gauss poj 1830

原创 2016年05月30日 18:20:02

POJ1830 开关问题

巧用位运算。
对于有效true的影响:0->1 , 1->0 ;
无效false的影响:0->0 , 1->1 ;
true ^ 0->1 , true^1->0 ;
false^0->0 , false^1->1 ;
已知结果矩阵,和初始矩阵,和关系矩阵。初始矩阵^(关系矩阵*变元) == 结果矩阵 变换->
变元*(关系矩阵)=初始矩阵^结果矩阵
变元的值域【0,1】。
方程解的个数 1<<变元个数

#include <iostream>
#include <cmath>
#include <cstdio>
#include <cstdlib>
using namespace std;
const int maxn = 50;
int a[maxn][maxn]; // zeng guang
int x[maxn] ; // ans_x
int free_x[maxn] ;
inline int gcd(int a , int b) {
    return b == 0 ? a : gcd(b,a%b) ;
    /*int t ;
    while(b){
        t = b ;
        b = a%b ;
        a = t ;
    }
    return a;*/
}
inline int lcm(int a ,int b){
    return a/gcd(a,b)*b;
}
inline int Abs(int a ){
    return a > 0 ? a : -a;
}
void debug(int equ , int var){
    for(int i = 0 ; i < equ ; ++i){
        for(int j = 0 ;j <= var ; ++j){
            cout << a[i][j] << " " ;
        }
        cout << endl;
    }
}
int Gauss(int equ , int var) {
    int max_r ;
    for(int i = 0 ; i <= var ; ++i){
        free_x[i] = false ;
        x[i] = 0 ;
    }
    //cout << equ << " " << var << endl;
    int row , col ;
    for(row = 0 ,col = 0 ; row < equ && col < var ; ++row , ++col) {
        max_r = row ;
        for(int i = row+1 ; i < equ ; ++i){         // find the max element in this col
            if(Abs(a[max_r][col]) < Abs(a[i][col])) max_r = i;
        }
        if(max_r != row){
            for(int j = col ; j <= var ; ++j){  // change
                swap(a[max_r][j],a[row][j]) ;
            }
        }
        if(!a[row][col]) {                  // 0
            row--;
            //cout << "123" <<endl;
            continue ;
        }
        //debug(3,3) ;
        for(int i = row+1 ; i < equ ; ++i){
            if(a[i][col]){// sub a[i][col] -> 0
                    /*int tmp = lcm(abs(a[i][col]),abs(a[row][col])) ;
                    int up = tmp/(abs(a[row][col])) ;
                    int down = tmp/(abs(a[i][col])) ;
                    //cout << tmp << " " << up << " " << down << endl;
                    if(a[i][col]*a[row][col]<0) up = -up ;
                for(int j = col ; j <= var; ++j){
                    a[i][j] = a[i][j]*down - a[row][j]*up ;
                }*/
                for(int j = 0 ; j <= var ; ++j){
                    a[i][j] ^= a[row][j] ;
                }
            }
        }
    //system("pause") ;
    }
    // col == var  , [row,equ) == 皆是0000...an
    for(int i = row ; i < equ ; ++i){
        if(a[i][col] != 0) {
            return -1 ;         // no solve
        }
    }
    if(row < var){ //[row , equ -1 ] -> (000.000)  have free_x ,at least var - row . cal the x[]
            //cout << row << " xxx " << endl;
        return var - row ;
    }
    return 0;
}
int s[maxn] , e[maxn] ;
void init(int var) {
    for(int i = 0 ;i < var ; ++i){
            cin >> s[i] ;
        }
        for(int i = 0 ;i < var ;++i){
            cin >> e[i] ;
            a[i][var] = s[i]^e[i] ;
        }
        for(int i = 0 ;i < var ; ++i){
            for(int j = 0 ;j < var ; ++j){
                a[i][j] = i==j?1:0;
            }
        }
}
int main(){
    int t ;
    cin >> t ;
    while(t--){
        int var ;
        cin >> var ;
        init(var) ;
        int i , j ;
        while(cin >> i >> j && (i|j)){          // I influence J
            a[j-1][i-1] = true ;
        }
        int flag = Gauss(var,var) ;
        if(flag == -1){
            cout << "Oh,it's impossible~!!" <<endl;
        }
        else cout << (1<<flag) <<endl;
    }
    return 0;
}
版权声明:本文为博主原创文章,未经博主允许不得转载。

poj 1830 高斯消元模版题

开关问题 Time Limit: 1000MS Memory Limit: 30000K Total Submissions: 6067 Accepted: 2309 Desc...

poj 1830 开关问题(guass)

这题用到了高斯消元法,有很多要注意的地方: a[i][j]表示第 j 个开关对第 i 个开关的影响,1表示有影响,0表示没有影响; 对矩阵消元的时候我们用对角线上的元素去把这列其他不为0的元素消去, ...

Guass消元 poj 1830

代码是拷别人的。。。。 举个例子 当n=3时 灯1  (a11*x1)^(a12*x2)^(a13*x3)=sta[0]^end[0] 灯2  (a21*x1)^(a22*x2)^(a23*x...

poj1830 开关问题 高斯消元

Language: Default 开关问题 Time Limit: 1000MS   Memory Limit: 30000K Total Submissio...

POJ 1830

#include #include #include #include using namespace std; const int maxn = 33; int a[maxn][maxn]; in...
  • zz_1215
  • zz_1215
  • 2012年07月23日 00:21
  • 545

poj 1830 开关问题

开关问题 Time Limit: 1000MS   Memory Limit: 30000K Total Submissions: 9328   Accepted:...
  • yjt9299
  • yjt9299
  • 2017年11月27日 21:02
  • 13

poj 1830 开关问题(高斯消元法)

Description 有N个相同的开关,每个开关都与某些开关有着联系,每当你打开或者关闭某个开关的时候,其他的与此开关相关联的开关也会相应地发生变化,即这些相联系的开关的状态如果原来为开就变为...

poj 1830开关问题

先送上参考学习网址: http://www.cppblog.com/menjitianya/archive/2014/06/08/207226.html

poj 1830 高斯消元

题意:有N个相同的开关,每个开关都与某些开关有着联系,每当你打开或者关闭某个开关的时候,其他的与此开关相关联的开关也会相应地发生变化,即这些相联系的开关的状态如果原来为开就变为关,如果为关就变为开。你...
  • bllsll
  • bllsll
  • 2017年11月02日 09:34
  • 34

POJ 1830 开关问题 高斯消元 异或方程

题意:给出N个开关,当改变一个开关z
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:Gauss poj 1830
举报原因:
原因补充:

(最多只允许输入30个字)