码迷,mamicode.com
首页 > 其他好文 > 详细

poj_1222_高斯消元

时间:2015-07-20 23:51:48      阅读:251      评论:0      收藏:0      [点我收藏+]

标签:高斯消元   数学   

第一次学习使用高斯消元,将灯板化为线性方程组,进行求解。

/*#########################################################################
# File Name: poj_1222.cpp
# Author: CaoLei
# Created Time: 2015/7/20 15:48:04
#########################################################################*/

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <set>
#include <queue>
#include <map>
using namespace std;
#define MAX(x,y) (((x)>(y))?(x):(y))
#define MIN(x,y) (((x)<(y))?(x):(y))
#define N 500010
#define pi acos(-1.0)
#define inf 100000000
typedef long long ll;
typedef unsigned long long ull;
typedef vector<double> vec;
typedef vector<vec> mat;
const double eps=1e-8;
int a[35][35];
int ans[35];
void gauss(){
    for(int i=0;i<30;i++){
        int tmp=0,k=0;
        for(int j=0;j<30;j++){
            if(a[j][i]>tmp) k=j;
        }
        for(int j=0;j<=30;j++){
            swap(a[i][j],a[k][j]);
        }
        for(int j=0;j<30;j++){
            if(i!=j&&a[j][i])
              for(int k=i+1;k<=30;k++){
                  a[j][k]=a[j][k]^a[i][k];
              }
        }
    }
}

int main(){
    freopen("in.txt","r",stdin);
    int t;
    scanf("%d",&t);
    for(int cnt=1;cnt<=t;cnt++){
        printf("PUZZLE #%d\n",cnt);
        memset(a,0,sizeof(a));
        for(int i=0;i<30;i++){
            scanf("%d",&a[i][30]);
            ans[i]=0;
        }
        for(int i=0;i<30;i++){
            a[i][i]=1;
            if(i%6!=0) a[i-1][i]=1;
            if(i%6!=5) a[i+1][i]=1;
            if(i>5) a[i-6][i]=1;
            if(i<24) a[i+6][i]=1;
        }
        gauss();
        for(int i=0;i<30;i++){
            ans[i]=a[i][30];
        }
        for(int i=0;i<30;i++){
            if((i+1)%6==0)
                printf("%d\n",ans[i]);
            else printf("%d ",ans[i]);
        }
    }
    return 0;
}



版权声明:本文为博主原创文章,未经博主允许不得转载。

poj_1222_高斯消元

标签:高斯消元   数学   

原文地址:http://blog.csdn.net/bryant03/article/details/46971339

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!