hdu 4592（高斯消元好题）

最新推荐文章于 2021-11-04 18:06:34 发布
最新推荐文章于 2021-11-04 18:06:34 发布
阅读量1.2k 收藏
点赞数
分类专栏： North--数论 North--高斯消元 北门的智慧——数论 文章标签： 高斯消元
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_27599517/article/details/52505600
版权
该博客介绍了如何使用求逆打表法解决数据范围大的高斯消元问题。由于普通高斯消元在大数据下效率低下，博主建议采用求逆矩阵的方法。题目链接指向HDU 4592，说明该题在特定情况下（4,5时）存在多解。代码实现中注意避免处理0,1矩阵以防止超时，最终得到的是十进制解决方案。" 108773708,10048130,Python多协程与队列爬取豆瓣图书Top250,"['Python', '编程语言', '网络爬虫']
摘要由CSDN通过智能技术生成

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4592；

题目分析：  本题其实就是一个普通的高斯消元题，但是本题数据范围巨大，所以不能用普通的高斯消元来写，只能采用求逆打表法，进行求解。只有满秩矩阵才能有逆元，通过消元可以看出，只有4,5时时多解的，所以其他就打表好了，但不能打0,1矩阵，还是会超时，最后打出的是十进制数就好办了。

代码如下：


#include <set>
#include <map>
#include <stack>
#include <queue>
#include <math.h>
#include <vector>
#include <string>
#include <utility>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <functional>

using namespace std;
int A2[4][2]=
{
    {3,2},{3,1},
    {2,3},{1,3}
};

int A3[9][3]=
{
    {5,1,6},{0,2,7},{5,4,3},
    {1,3,1},{2,7,2},{4,6,4},
    {6,1,5},{7,2,0},{3,4,5}
};

int A6[36][6]=
{
    {43,8,55,10,44,40},{3,20,53,17,14,4},{44,42,7,32,55,34},{13,21,56,1,59,17},{48,10,43,34,28,8},{53,4,59,20,13,5},
    {8,28,2,27,34,44},{20,54,5,59,21,14},{42,43,10,48,2,55},{21,53,20,3,16,59},{10,27,40,55,42,28},{4,14,16,54,17,13},
    {55,2,48,42,27,10},{53,5,56,1,59,17},{7,10,28,40,48,32},{56,20,14,5,3,1},{43,40,7,32,55,34},{59,16,3,21,54,20},
    {10,27,42,48,2,55},{17,59,1,56,5,53},{32,48,40,28,10,7},{1,3,5,14,20,56},{34,55,32,7,40,43},{20,54,21,3,16,59},
    {44,34,27,2,28,8},{14,21,59,5,54,20},{55,2,48,10,43,42},{59,16,3,20,53,21},{28,42,55,40,27,10},{13,17,54,16,14,4},
    {40,44,10,55,8,43},{4,14,17,53,20,3},{34,55,32,7,42,44},{17,59,1,56,21,13},{8,28,34,43,10,48},{5,13,20,59,4,53}
};

int A7[49][7]=
{
    {109,65,14,84,88,16,96},{108,34,27,2,28,40,112},{3,20,53,81,110,84,56},{107,0,107,8,119,42,28},{96,20,86,69,59,21,14},
    {27,34,108,32,28,10,7},{91,65,56,21,13,4,3},{65,99,21,86,68,56,16},{34,119,32,7,42,108,40},{20,54,69,91,5,86,84},
    {0,0,8,28,34,107,42},{20,54,81,109,80,53,21},{34,119,2,112,42,27,10},{65,99,84,53,17,14,4},{14,21,59,5,54,68,88},
    {27,32,107,8,119,42,28},{53,69,88,17,99,5,110},{107,8,119,34,0,34,119},{86,81,13,68,99,80,59},{108,2,107,8,119,42,28},
    {56,84,110,80,54,17,13},{84,86,5,91,5,86,84},{2,7,8,27,8,7,2},{81,91,17,96,17,91,81},{8,28,34,107,34,28,8},
    {69,109,68,3,68,109,69},{32,112,8,108,8,112,32},{21,53,80,109,80,53,21},{88,68,54,5,59,21,14},{28,42,119,8,107,32,27},
    {110,5,99,17,88,69,53},{119,34,0,34,119,8,107},{59,80,99,68,13,81,86},{28,42,119,8,107,2,108},{13,17,54,80,110,84,56},
    {16,56,68,86,21,99,65},{40,108,42,7,32,119,34},{84,86,5,91,69,54,20},{42,107,34,28,8,0,0},{21,53,80,109,81,54,20},
    {10,27,42,112,2,119,34},{4,14,17,53,84,99,65},{96,16,88,84,14,65,109},{112,40,28,2,27,34,108},{56,84,110,81,53,20,3},
    {28,42,119,8,107,0,107},{14,21,59,69,86,20,96},{7,10,28,32,108,34,27},{3,4,13,21,56,65,91}
};

int A8[64][8]=
{
    {220,138,7,128,199,170,108,40},{214,81,13,64,237,1,238,68},{27,0,27,32,107,168,199,130},{237,81,54,16,14,5,3,1},
    {183,138,108,8,112,160,192,128},{216,0,216,4,214,21,227,65},{107,138,176,2,183,128,119,34},{59,81,224,1,227,85,54,20},
    {138,219,10,192,42,171,130,108},{81,219,17,224,65,3,69,238},{0,0,32,112,136,172,42,199},{81,219,65,56,21,13,4,3},
    {138,219,130,28,168,176,32,192},{0,0,4,14,17,53,84,227},{138,219,136,7,130,192,162,119},{81,219,80,3,84,213,65,54},
    {7,10,28,160,172,2,171,170},{13,17,54,16,14,5,3,1},{27,32,107,136,183,10,172,168},{54,65,213,68,59,20,13,5},
    {108,130,171,34,220,40,176,160},{216,4,214,17,237,80,53,21},{176,136,108,8,112,160,192,128},{224,80,56,5,53,64,213,85},
    {128,192,160,112,136,172,42,199},{64,224,16,216,84,14,65,237},{32,112,136,172,10,183,136,107},{16,56,68,214,69,59,21,14},
    {8,28,34,107,162,220,168,112},{4,14,17,53,80,237,17,214},{2,7,8,27,42,112,130,183},{1,3,5,14,17,53,84,227},
    {199,42,172,136,112,160,192,128},{237,65,14,84,216,16,224,64},{107,136,183,10,172,136,112,32},{14,21,59,69,214,68,56,16},
    {112,168,220,162,107,34,28,8},{214,17,237,80,53,17,14,4},{183,130,112,42,27,8,7,2},{227,84,53,17,14,5,3,1},
    {170,171,2,172,160,28,10,7},{1,3,5,14,16,54,17,13},{168,172,10,183,136,107,32,27},{5,13,20,59,68,213,65,54},
    {160,176,40,220,34,171,130,108},{21,53,80,237,17,214,4,216},{128,192,160,112,8,108,136,176},{85,213,64,53,5,56,80,224},
    {108,130,171,42,192,10,219,138},{238,69,3,65,224,17,219,81},{199,42,172,136,112,32,0,0},{3,4,13,21,56,65,219,81},
    {192,32,176,168,28,130,219,138},{227,84,53,17,14,4,0,0},{119,162,192,130,7,136,219,138},{54,65,213,84,3,80,219,81},
    {40,108,170,199,128,7,138,220},{68,238,1,237,64,13,81,214},{130,199,168,107,32,27,0,27},{1,3,5,14,16,54,81,237},
    {128,192,160,112,8,108,138,183},{65,227,21,214,4,216,0,216},{34,119,128,183,2,176,138,107},{20,54,85,227,1,224,81,59}
};
const int maxn=405;
int a[maxn][maxn*2+1],x[maxn];
int equ,var;
int free_num,ans=100000000;
void Debug(){
    int i,j;
    for(i=0;i<equ;i++){
        for(j=0;j<var+1;j++)
        cout<<a[i][j]<<" ";
        cout<<endl;
    }
    cout<<endl;
}

int Gauss(){
    int k,col=0;
    for(k=0;k<equ&&col<var;++k,++col) {
        int max_r=k;
        for(int i=k+1;i<equ;++i)
        if(a[i][col]>a[max_r][col])
        max_r=i;
        if(max_r!=k){
            for(int i=k;i<var+1;++i)
            swap(a[k][i],a[max_r][i]);
        }
        if(a[k][col]==0){
            k--;
            continue;
        }
        for(int i=k+1;i<equ;++i) {
            if(a[i][col] != 0) {
                for(int j=col;j<var+1;++j)
                a[i][j]^=a[k][j];
            }
        }
    }
    for(int i=k;i<equ;++i)
    if(a[i][col]!=0)return -1;
//    Debug();
    for(int i=0;i<equ;++i)
    if(!a[i][i]) {
        int j;
        for(j=i+1;j<var;++j)
        if(a[i][j]) break;
        if(j==var) break;
        for(int k=0;k<equ;++k)
        swap(a[k][i],a[k][j]);
    }
    int sum=0x3f3f3f3f;
    if (var-k>0){
        for(int num=0;num<(1<<(var-k));num++){
            int minn=0;
            memset(x,0,sizeof(x));
            for(int ans=0;ans<var-k;ans++){
                if((1<<ans)&num){
                    x[var-ans-1]=1;
                    minn++;
                }
            }
            for(int i=k-1;i>=0;--i)
            {
                int tmp=a[i][var];
                for(int j=i+1;j<var;++j)
                if(a[i][j]!=0)
                tmp^=x[j];
                x[i]=tmp;
                minn+=x[i];
            }
            sum=min(minn,sum);
        }
    }
    if (var-k==0){
        sum=0;
        for(int i=k-1;i>=0;--i)
        {
            int tmp=a[i][var];
            for(int j=i+1;j<var;++j)
            if(a[i][j]!=0)
            tmp^=a[i][j]*x[j];
            x[i]=tmp;
        }
        for(int i=0;i<var;i++)
        sum+=x[i];
    }
    return sum;
}
void init(int n){
    equ=n*n;
    var=equ;
    memset(a,0,sizeof(a));
    memset(x,0,sizeof(x[0])*30);
    for(int i=0; i<n; i++)
    for(int j=0; j<n; j++){
        if(i!=0) a[i*n+j][(i-1)*n+j]=1;
        if(i!=n-1) a[i*n+j][(i+1)*n+j]=1;
        if(j!=0) a[i*n+j][i*n+j-1]=1;
        if(j!=n-1) a[i*n+j][i*n+j+1]=1;
        a[i*n+j][i*n+j]=1;
    }
}
int main(){
    int cas,n,X;
    scanf("%d",&cas);
    while(cas--){
        scanf("%d",&n);
        if(n==1){
            scanf("%d",&X);
            printf("%d\n",X);
            continue;
        }
        if(n==4||n==5){
            init(n);
            for(int i=0;i<n;i++){
                scanf("%d",&X);
                int r=(i+1)*n-1,l=i*n;
                while(r>=l){
                    a[r][var]=X&1;
                    X>>=1;
                    r--;
                }
            }
            int ans=Gauss();
            printf("%d\n",ans);
            continue;
        }
        int ans=0;
        int t[80];
        for(int i=0; i<n; i++)
        scanf("%d",&t[i]);
        equ=n*n;
        memset(x,0,sizeof(x[0])*(equ+4));
        if(n==2)
        for(int j=0; j<n; j++)
        if(t[j])
        for(int i=0; i<equ; i++)
        x[i]=x[i]^(A2[i][j]&t[j]);
        else if(n==3)
        for(int j=0;j<n;j++)
        if(t[j])
        for(int i=0; i<equ; i++)
        x[i]=x[i]^(A3[i][j]&t[j]);
        else if(n==6)
        for(int j=0; j<n; j++)
        if(t[j])
        for(int i=0; i<equ; i++)
        x[i]=x[i]^(A6[i][j]&t[j]);
        else if(n==7){
            for(int j=0; j<n; j++)
            if(t[j]){
                for(int i=0; i<equ; i++)
                    x[i]=x[i]^(A7[i][j]&t[j]);
            }
        }
        else if(n==8){
            for(int j=0;j<n;j++)
            if(t[j]){
                for(int i=0; i<equ; i++){
                    x[i]=x[i]^(A8[i][j]&t[j]);
                }
            }
        }
        for(int i=0; i<equ; i++){
            int t=0;
            while(x[i]){
                t^=(x[i]&1);
                x[i]>>=1;
            }
            ans+=t;
        }
        printf("%d\n",ans);
    }
    return 0;
}


leader_win
关注
专栏目录
hdu 3359 高斯消元法
_rabbit的专栏
03-07 1485
Kind of a Blur Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 1226    Accepted Submission(s): 356 Problem Description Image blurring o
HDU4418——高斯消元+DP
Dr.Cell--CellTech
03-11 1404
HDU4418 题目大意：一个人站在编号为0..n-1的数轴上，以p[i]的概率走i步(i∈[1,m]，且∑p[i]=1)给定n,m,起点s,终点e,方向d,求出这人走到终点的步数期望。 题目难点：列方程求解；不要用pascal 题目思路： 首先按一般思路先列出f数组表示期望（废话……） 显然f[e]=0,f[i]=∑((f[x+j]+j)*p[j]) 然后把第二条式子变形，则得到：
HDU 4592 Boring Game(高斯消元好题)
代码改变世界
09-11 1338
传送门 Boring GameTime Limit: 8000/4000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others)Total Submission(s): 827    Accepted Submission(s): 80Problem Description Zero likes to play a boring
高斯消元基础题
ACdreamer
05-31 4073
题目：EXTENDED LIGHTS OUT   题意：5*6的一个由灯组成的方阵 操作一个灯 周围的上下左右四个灯会发生相应变化 即由灭变亮 由亮变灭 问如何操作使灯全灭。 本题方法：猛戳这里 #include #include #include #include using namespace std; const int N=105; int eq
hdu 3364 Lanterns 高斯消元，模板题
冷月之殇
09-17 2590
#include #include #include #include #include #include #include #include using namespace std; const int maxn=55; struct matrix{ int f[maxn][maxn]; }e,g; int find(matrix a,int m,int n)//高斯消元
【例题】【高斯消元(异或)】
ANG
12-19 1989
1、 NKOJ1985 中央暖气 时间限制 : 10000 MS 空间限制 : 65536 KB问题描述 冬天来了，但URAL大学的暖气系统还没启动。这个暖气系统包括许多阀门，只有所有阀门都打开了才能供应暖气。大学里有一些技术员，他们每人负责一个或多个阀门，有可能存在一个阀门由多个人负责的情况。当一个技术员接到启动暖气系统的指示后他就会把所有自己负责的阀门都转动一次。这意味着如果一个阀门
HDU3359(高斯消元浮点数版）
creatorx的博客
06-12 361
Kind of a BlurTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 2568 Accepted Submission(s): 694Problem Description Image blurring occurs whe
期望与概率DP专题/高斯消元
Qianyri的博客
08-19 1398
高斯消元 总结：概率和期望问题的求解都要依靠递推式，概率问题按照递推式进行dp，期望问题根据递推式化简的难易程度通过化简再dp或高斯消元求解，该类问题大部分仍属于各个类型的dp求解问题 2018.9.3-概率与期望dp学习 Discovering Gold 一排1到n的格子，每个格子上有黄金 ai ,你最开始在 1 号,每一次投骰子决定到哪一个格子，超出1~n范围则重新投掷，你到了哪个格子...
高斯消元配合概率dp-图上随机游走模型
qq_35577488的博客
01-03 1019
题目大意: 给你一张无向连通图,从第111点随机游走到nnn点。花费为经过的边的编号的和.问你如何安排边的编号使得期望花费最小化. n≤500,m≤n2n \leq 500,m \leq n^2n≤500,m≤n2 题目思路: 很容易想到，求出每条边期望经过次数，然后贪心安放即可. 但是mmm太大，无法求解.考虑先求每个点期望经过次数. 考虑到nnn点就停止了，所以不考虑nnn的贡献。 显然有转移: f1=1+∑(1,j)∈Efjdujf_1=1+\sum_{(1,j)\in E}^{}\frac{f_j}
hdu 3359 /poj 3999 Kind of a Blur(高斯消元模板题)
TUT我好菜啊
08-14 864
Kind of a Blur Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 2242    Accepted Submission(s): 609 Problem Description Image blurri
高斯消元好题选讲
wljoi的博客
10-24 283
例1：[SDOI2006]线性方程组 传送门 例2：[JSOI2008]球形空间产生器 传送门 例3：Barracuda 传送门 例4 传送门 Solution 例1： 看似是一个高斯消元，但普通的程序交上去会被卡成808080分WAWAWA。 我们看下面的一组hack数据： 2 0 2 3 0 0 0 化简2y+3=0,0=02y+3=0,0=02y+3=0,0=0 易得该方程组有无数解，但是...
HDU5088（高斯消元）
cq_phqg
11-04 929
Revenge of Nim II Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 443    Accepted Submission(s): 157 Problem Description Nim is a m
hdu 3915 高斯消元
weixin_30555515的博客
11-01 92
http://acm.hdu.edu.cn/showproblem.php?pid=3915 这道题目是和博弈论挂钩的高斯消元。本题涉及的博弈是nim博弈，结论是：当先手处于奇异局势时（几堆石子数相互异或为0），其必败。 思路在这里，最后由于自由变元能取1、0两种状态，所以，最终答案是2^k，k表示自由变元的个数。 #include <cstdio> #inclu...
【高斯消元】[NOIP2004]虫食算 —— 正解
weixin_30614109的博客
01-30 229
题目描述：所谓虫食算，就是原先的算式中有一部分被虫子啃掉了，需要我们根据剩下的数字来判定被啃掉的字母。来看一个简单的例子： +43a9865a045008468a663344445506978 其中a号代表被虫子啃掉的数字。 根据算式，我们很容易判断：第一行的两个数字分别是5和3，第二行的数字是5。 现在，我们对问题做两个限制： 首先，我们只考虑加法的...
HDU 5755 高斯消元
weixin_30906701的博客
07-30 98
http://acm.hdu.edu.cn/showproblem.php?pid=5755 题意：一个N*M的棋盘，每个点是0,1,2中的一个，在取模3的含义下，翻转一个这个+2，周围的4个加1.。。。 思路：感慨一下知识面太狭隘了。。学习了一下高斯消元，这个题目其实就是高斯消元的一个经典题目。。每个结点的翻转次数都是一个变量。。猛地一看时间复杂度是n^3m^3但是高斯消元里面变...
高斯消元专题训练
HeartWine的博客
11-04 488
文章目录POJ 1222 EXTENDED LIGHTS OUT POJ 1222 EXTENDED LIGHTS OUT 题意： 给定一个555行666列的010101矩阵，每一个小格子代表一个开关，000表示开关关闭，111表示开关闭合，当按下一个开关后，其也会影响该开关上下左右的开关的状态。现在给定一个这样的01矩阵，输出一种方案，使得所有开关都变成闭合状态。 思路：555行666列，相当于一共是303030个未知数，可以根据题意列出303030个方程，然后高斯消元，因为都是010101，所以是解异或
hdu 3976 高斯消元 模板
转角的守望
11-22 482
Electric resistance Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 655    Accepted Submission(s): 346 Problem Description Now give you
hdu 4592 Boring Game
浮躁的人，终究一事无成！
09-03 1743
Boring Game Time Limit: 8000/4000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others) Total Submission(s): 544    Accepted Submission(s): 91 Problem Description Zero likes to pla
HDU 3364 (高斯消元)
morejarphone~
12-08 428
求出每种状态下的自由元个数n,答案就是2^n. #include #include #include #include #include using namespace std; #define maxn 333 #define free Free int a[maxn][maxn]; int aa[maxn][maxn]; int n, m; int free[maxn*maxn]
HDU动态规划：46题实战分析与背包问题详解
动态规划是计算机科学中解决优化问题的一种重要算法策略，特别是在处理具有重叠子问题和最优子结构特性的算法中。在HDU竞赛中，动态规划...熟练掌握动态规划思想和技巧，能够帮助参赛者在这些挑战性问题中取得好成绩。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值