poj2345


题目描述:

一个工人搞几个暖气,最后都打开

题解:

%2等于1的高斯消元

重点:

用^代替加法.

代码:
#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <cmath>
#include <ctype.h>
#include <limits.h>
#include <cstdlib>
#include <algorithm>
#include <vector>
#include <queue>
#include <map>
#include <stack>
#include <set>
#include <bitset>
#define CLR(a) memset(a, 0, sizeof(a))
#define REP(i, a, b) for(int i = a;i < b;i++)
#define REP_D(i, a, b) for(int i = a;i <= b;i++)

typedef long long ll;

using namespace std;

const int maxn = 250 + 100;
int a[maxn][maxn];
int x[maxn];
int n;
int equ, var;

int gauss()
{
    for(int k = 0, col=0;k<equ&&col<var;k++,col++)
    {
        int max_k = k;
        for(int i=k+1;i<equ;i++)
        {
            if(a[max_k][col]==1)
            {
                break;
            }
            if(a[i][col]==1)
            {
                max_k = i;
                break;
            }
        }
        if(a[max_k][col]==0)
        {
            return 0;
        }
        if(max_k!=k)
        {
            swap(x[k], x[max_k]);
            for(int i = col;i<var;i++)
            {
                swap(a[k][i], a[max_k][i]);
            }
        }
        for(int i=0;i<equ;i++)
        {
            if(i!=k)
            {
                if(a[i][col]==1)
                {
                    a[i][col] = 0;
                    for(int j=col+1;j<var;j++)
                    {
                        a[i][j] ^= a[k][j];
                    }
                    x[i] ^= x[k];
                }
            }
        }
    }
    return 1;
}
void solve()
{
    CLR(a);
    REP(i, 0, n)
    {
        int d;
        while(scanf("%d", &d)&&d!=-1)
        {
            d--;
            a[d][i] = 1;
        }
    }
    REP(i, 0, n)
    {
        x[i] = 1;
    }
    equ = n;
    var = n;
    gauss();
    vector<int> ans;
    REP(i, 0, n)
    {
        if(x[i])
        {
            ans.push_back(i+1);
        }
    }
    REP(i, 0, ans.size())
    {
        if(i==ans.size()-1)
        {
            printf("%d\n", ans[i]);
        }
        else
        {
            printf("%d ", ans[i]);
        }
    }
}

int main()
{
    //freopen("10Jin.txt", "r", stdin);
    //freopen("10Jout.txt", "w", stdout);
    while(scanf("%d", &n) != EOF)
    {
        solve();
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值