• POJ1830 开关问题 高斯消元


    这题是给定N个灯的初始和最终状态,再给定一些关系,这些关系说明按某个开关可能影响其他的灯的开关情况,可以将这种关系视为一种取反的关系。

    对于这题我们假设一组数据:

    3
    0 1 0
    1 1 0
    1 2
    2 3
    1 3
    3 2
    0 0

    对于以上的数据,我们用矩阵来表示整个操作的过程和状态 

       0       1

    S =   1           E =    1        

       0         0

    我们可以很显然的知道,某个等的亮灭情况将体现出各个开关的复合结果,因此我们可以得到这样的方程组,模二加等价于异或操作:

    E(a) = xa*A11 ^ xb*A12 ^ xc*A13 ^ S(a);

    E(b) = xa*A21 ^ xb*A22 ^ xc*A23 ^ S(b);

    E(c) = xa*A31 ^ xb*A32 ^ xc*A33 ^ S(c);

    其中S(a)表示初始状态,E(a)表示a号灯最终的情况,xa表示a号开关是否按下,A13表示3号开关是否影响1号灯的亮灭,异或操作体现了操作的具体影响。

    将S项移到方程的左边,我们抽象出系数矩阵,然后构造出增广矩阵,然后就是进行高斯消元,高斯消元的意思可以看作是符合一系列两个逻辑关系的情况下,对变元要求的变换。只不过这里用的全部都是异或操作,最后我们进行秩的判定,如果矩阵的秩和邻接矩阵的秩不相等的话那么就无解,如果相同的话,那么说明达到理想灯的亮灭情况下至少需要确定的状态的开关数为矩阵的秩。那么解就是剩下开关的任意选择了,即2^X(每个开关对应关或者是开)。

    代码如下:

    #include <cstdlib>
    #include <cstring>
    #include <cstdio>
    #include <algorithm>
    using namespace std;
    
    int N, s[50], e[50];
    
    void swap(int &a, int &b)
    {
        int t = a;
        a = b;
        b = t;
    }
    
    struct Matrix
    {
        int a[50][50];    
        void rswap(int x, int y)
        {
            for (int j = 1; j <= N+1; ++j) {
                swap(a[x][j], a[y][j]);  
            }
        }
        void relax(int x, int y) 
        {
            for (int j = 1; j <= N+1; ++j) {
                a[y][j] ^= a[x][j];
            }
        }
        void init(){
            memset(a, 0, sizeof (a));
            for (int i = 1; i <= N; ++i) {
                a[i][i] = 1;
            }
        }
    }M;
    
    void Gauss()
    {
        int ptr, t, i, j;
        for (i = 1, j = 1; i <= N, j <= N; ++i, ++j) {
            for (ptr = i; ptr <= N; ++ptr) {
                if (M.a[ptr][j]) {
                    // 当第i行的j列不为零
                    break;
                }
            }
            if (ptr == N+1) {
                --i;
                continue;
            }
            if (ptr != i) {
                M.rswap(ptr, i);
            }
            for (int k = i + 1; k <= N; ++k) {
                if (!M.a[k][j]) {
                    continue;
                }
                M.relax(i, k);
            }
        }
        for (int k = i; k <= N; ++k) {
            if (M.a[k][N+1] != 0) {
                puts("Oh,it's impossible~!!");
                return;
            }
        }
        printf("%d\n", 1<<(N-i+1));
    }
    
    int main()
    {
        int T, a, b;
        scanf("%d", &T);
        while (T--) {
            scanf("%d", &N);
            M.init();
            for (int i = 1; i <= N; ++i) {
                scanf("%d", &s[i]);
            }
            for (int i = 1; i <= N; ++i) {
                scanf("%d", &e[i]);
                M.a[i][N+1] = s[i]^e[i];
            }
            while (scanf("%d %d", &a, &b), a | b) {
                M.a[b][a] = 1;
            }
            Gauss();
        }
        return 0;
    }
  • 相关阅读:
    DotnetCore 使用Jwks验证JwtToken签名
    HashCode
    C# RedisRateLimiter
    Centos7 使用Dockerfile 制作自己的Dotnetcore程序镜像
    ES6 HttpApplication Middleware
    请转发!简单2分钟制作无接触式小区进出微信登记表!全免费!数据安全!所有数据均存在创建人登录的QQ腾讯文档里!
    理解虚基类、虚函数与纯虚函数的概念
    不无道理
    乔布斯:不要为明天忧虑!
    【心态不好,人生易老】
  • 原文地址:https://www.cnblogs.com/Lyush/p/2609271.html
Copyright © 2020-2023  润新知