poj 1681

最新推荐文章于 2020-12-31 13:40:43 发布

转载最新推荐文章于 2020-12-31 13:40:43 发布 · 92 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/zhaozhe/archive/2011/05/25/2057258.html

本文详细介绍了一种解决线性方程组的有效方法——高斯消元法，并提供了完整的C++实现代码。通过该方法可以有效地求解由多个线性方程组成的方程组，适用于各种科学计算和工程应用中。

高斯消元

代码：

#include<iostream>
#include<fstream>

using namespace std;

int move[5][2]={0,1,0,-1,1,0,-1,0,0,0};
int n,m;
int a[300][300];

void solve(){
    int i,j,k,s;
    j=0;
    for(i=0;i<n&&j<n;i++,j++)
    {
        for(k=i;k<n;k++)
            if(a[k][j])
                break;
        if(k==n){
            i--;
            continue;
        }
        if(k!=i)
        {
            for(s=j;s<=n;s++)
                swap(a[i][s],a[k][s]);
        }
        for(k=i+1;k<n;k++)
            if(a[k][j])
                for(s=j;s<=n;s++)
                    a[k][s]=a[k][s]^a[i][s];
    }
    k=i;
    for(i=k;i<n;i++)
        if(a[i][n])
        {
            cout<<"inf"<<endl;
            return;
        }
    s=0;
    for(i=k-1;i>=0;i--)
    {
        for(j=i+1;j<n;j++)
            a[i][n]^=(a[i][j]&&a[j][n]);
        if(a[i][n]) s++;
    }
    cout<<s<<endl;
}

void read(){
//    ifstream cin("in.txt");
    int i,j,k,s,t;
    int cas;
    cin>>cas;
    char c;
    while(cas--)
    {
        cin>>m;
        n=m*m;
        memset(a,0,sizeof(a));
        for(i=0;i<m;i++)
            for(j=0;j<m;j++)
            {
                cin>>c;
                if(c=='y')
                    a[i*m+j][n]=0;
                else
                    a[i*m+j][n]=1;
            }
        for(i=0;i<m;i++)
            for(j=0;j<m;j++)
            {
                for(k=0;k<5;k++)
                {
                    s=i+move[k][0];
                    t=j+move[k][1];
                    if(s>=0&&s<m&&t>=0&&t<m)
                        a[i*m+j][s*m+t]=1;
                }
            }
        solve();
    }

}

int main(){
    read();
    return 0;
}

转载于:https://www.cnblogs.com/zhaozhe/archive/2011/05/25/2057258.html