高斯消元练习

最新推荐文章于 2024-05-02 17:34:07 发布
最新推荐文章于 2024-05-02 17:34:07 发布
阅读量531 收藏
点赞数
分类专栏： algorithm_数论 文章标签： 数学
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/theArcticOcean/article/details/50686001
版权
POJ 1222 EXTENDED LIGHTS OUT
开关灯问题：每一个开关对四周5点领域都有影响，状态变化：如果灯是亮着的，熄灭；如果灯是熄灭的，亮起。把每一个开关的开闭看做一个未知数x，影响领域看做系数a，每一盏灯现在的状态是y，问题就是 . 等式解释：如果灯是关着的，那么 必须等于0，维持状态；如果灯是开着的， 必须等于1，改变状态。
对于“ 5点领域”换个角度看，只有周围的5盏灯会对其产生影响，所以他们的当前系数是1。有：
假设有规模为3*3的9盏灯：
0 1 0
1 0 0
1 1 1
增广矩阵对应的是：
1 1 0 1 0 0 0 0 0 0
1 1 1 0 1 0 0 0 0 1
0 1 1 0 0 1 0 0 0 0
1 0 0 1 1 0 1 0 0 1
0 1 0 1 1 1 0 1 0 0
0 0 1 0 1 1 0 0 1 0
0 0 0 1 0 0 1 1 0 1
0 0 0 0 1 0 1 1 1 1
0 0 0 0 0 1 0 1 1 1 

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
using namespace std;
const int N=905;
int gcd(int a,int b){
    return b==0?a:gcd(b,a%b);
}
// a[][]：增广矩阵 n: 方程数 m：变量数
void Gauss(int a[N][N],int n,int m,int &r,int &c){
    for(r=0,c=0;r<n&&c<m;r++,c++){
        int maxi=r;
        for(int i=r+1;i<n;i++){
            if(abs(a[i][c])>abs(a[maxi][c])) maxi=i;
        }
        if(maxi!=r){
            for(int j=r;j<m+1;j++){ 
                swap(a[r][j],a[maxi][j]);
            }
        }
        if(a[r][c]==0){
            r--;
            continue;
        }
        for(int i=r+1;i<n;i++){
            if(a[i][c]!=0){
                int x=abs(a[i][c]);
                int y=abs(a[r][c]);
                int lcm=x/gcd(x,y)*y;
                int tx=lcm/x;
                int ty=lcm/y;
                if(a[i][c]*a[r][c]<0) ty=-ty;
                for(int j=c;j<m+1;j++){
                    a[i][j]=(a[i][j]*tx-a[r][j]*ty+2)%2;
                }
            }
        }
    }
}
int reback(int a[][N],int n,int m,int r,int c,int x[]){
    for(int i=r-1;i>=0;i--){
        int t=a[i][c];
        for(int j=i+1;j<c;j++)  t-=a[i][j]*x[j];
        x[i]=t/a[i][i]%2;
        x[i]=(x[i]+2)%2;
    }
    return 0;
}

int a[N][N];
int x[N];

int main(){
    int t;
    int ca=1;
    scanf("%d",&t);
    while(t--){
        memset(a,0,sizeof(a));
        for(int i=0;i<5;i++){
            for(int j=0;j<6;j++){
                a[6*i+j][6*i+j]=1;
                if(i>0) a[6*i+j][6*(i-1)+j]=1;
                if(i<4) a[6*i+j][6*(i+1)+j]=1;
                if(j>0) a[6*i+j][6*i+j-1]=1;
                if(j<5) a[6*i+j][6*i+j+1]=1;
            }
        }
        for(int i=0;i<5;i++){
            for(int j=0;j<6;j++) scanf("%d",&a[6*i+j][30]);
        }
        int r,c;
        Gauss(a,30,30,r,c);
        reback(a,30,30,r,c,x);
        printf("PUZZLE #%d\n",ca++);
        for(int i=0;i<30;i++){
            if((i+1)%6==0) printf("%d\n",x[i]);
            else printf("%d ",x[i]);
        }
    }
    return 0;
}

POJ 1830 开关问题
大意：给出所有灯的开始和结束状态，开关之间的影响关系，求出有多少种方案能使得状态改变成功
和上一题有点类似。用数学表达：
如果存在自由变元，答案则是2^ans
如果方案唯一，则是2^0。总之就是1<<ans 

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int N=30;
int gcd(int a,int b){
    return b==0?a:gcd(b,a%b);
}
void Gauss(int a[N][N],int n,int m,int &r,int &c){
    for(r=0,c=0;r<n&&c<m;r++,c++){
        int maxi=r;
        for(int i=r+1;i<n;i++){
            if(a[i][c]>a[maxi][c]) maxi=i;
        }
        if(maxi!=r){
            for(int j=r;j<m+1;j++){
                swap(a[r][j],a[maxi][j]);
            }
        }
        if(a[r][c]==0){
            r--;
            continue;
        }
        for(int i=r+1;i<n;i++){
            if(a[i][c]&&a[r][c]){
                for(int j=c;j<m+1;j++){
                    a[i][j]=(a[i][j]-a[r][j]+2)%2;
                }
            }
        }
    }
}
// -2 没有整数解 -1：无解  0: 唯一解 m-r：（自由变元个数）无穷多个解
int reback(int a[][N],int n,int m,int r,int c,int x[],bool tag[]){
    for(int i=r;i<n;i++) if(a[i][c]!=0) return -1;
    if(r<m){
        memset(tag,0,sizeof(tag));
        for(int i=r-1;i>=0;i--){
            int dex=0;
            int cnt=0;
            for(int j=0;j<m;j++){
                if(a[i][j]!=0&&tag[j]==0){
                    cnt++;
                    dex=j;
                }
            }
            if(cnt>1) continue;
            int t=a[i][m];
            for(int j=0;j<m;j++){
                if(a[i][j]!=0&&j!=dex){
                    t-=a[i][j]*x[j];
                }
            }
            x[dex]=t/a[i][dex]%2;
            x[dex]=(x[dex]+2)%2;
            tag[dex]=1;
        }
        return m-r;
    }
    for(int i=r-1;i>=0;i--){
        int t=a[i][c];
        for(int j=i+1;j<c;j++)  t-=a[i][j]*x[j];
        if(t%a[i][i]!=0) return -2;
        x[i]=t/a[i][i]%2;
        x[i]=(x[i]+2)%2;
    }
    return 0;
}
int e1[N],e2[N];
int a[N][N];
int x[N];
bool tag[N];
int main()
{
    //freopen("cin.txt","r",stdin);
    int k;
    cin>>k;
    while(k--){
        memset(a,0,sizeof(a));
        int N;
        scanf("%d",&N);
        for(int i=0;i<N;i++) scanf("%d",&e1[i]);
        for(int i=0;i<N;i++) scanf("%d",&e2[i]);
        for(int i=0;i<N;i++) {
            a[i][N]=(e2[i]-e1[i]+2)%2;
            a[i][i]=1;
        }
        int p1,p2;
        while(scanf("%d%d",&p1,&p2)){
            if(p1+p2==0) break;
            a[p2-1][p1-1]=1;
        }
        int r,c;
        Gauss(a,N,N,r,c);
        int g=reback(a,N,N,r,c,x,tag);
        if(g<0) printf("Oh,it's impossible~!!\n");
        else printf("%d\n",1<<g);
    }
    return 0;
}

POJ 2947 Widget Factory
大意：周期性时间问题，建立模线性方程组。每一个装饰物的加工应该是需要一定人手的，且有一定的工作时间。设每一个装饰物的加工时间为 , 有
例如例一： 

#include <iostream>
#include <cstdio>
#include <map>
#include <algorithm>
#include <cstring>
using namespace std;
const int N=305;

int gcd(int a,int b){
    return b==0?a:gcd(b,a%b);
}

void Gauss(int a[N][N],int n,int m,int &r,int &c){
    for(r=0,c=0;r<n&&c<m;r++,c++){
        int maxi = r;
        for(int i=r+1; i<n; i++)
            if(a[i][c] > a[maxi][c])
                maxi = i;
        if(maxi != r)
        {
            for(int i=r; i<m+1; i++)
                swap(a[r][i],a[maxi][i]);
        }
        if(a[r][c] == 0)
        {
            r--;
            continue;
        }
        for(int i=r+1; i<n; i++)
        {
            if(a[i][c] != 0)
            {
                int LCM = a[i][c]/gcd(a[i][c],a[r][c])*a[r][c];
                int tx = LCM / a[i][c];
                int ty = LCM / a[r][c];
                for(int j=c; j<m+1; j++)
                    a[i][j] = ((a[i][j] % 7 * tx % 7 - a[r][j] % 7 * ty % 7) % 7 + 7) % 7;
            }
        }
    }
}
//-2 没有整数解 这一情况除去，模线性方程中没有整除的概念
//  -1：无解  0: 唯一解 m-r：（自由变元个数）无穷多个解
int reback(int a[][N],int n,int m,int r,int c,int x[],bool tag[]){
    for(int i=r;i<n;i++) if(a[i][c]!=0) return -1;
    if(r<m) return m-r;
    for(int i=r-1;i>=0;i--){
        int t=a[i][c];
        for(int j=i+1;j<c;j++)  {
            t-=a[i][j]*x[j]%7;
            t=(t%7+7)%7;
        }
        //if(t%a[i][i]!=0) return -2;
        //int ni,y;
        //exgcd(t,a[i][i],ni,y);
        //ni=(ni%7+7)%7;
        //x[i]=t*ni%7;
        while(t%a[i][i]!=0) t+=7;
        x[i]=t/a[i][i]%7;
        if(x[i]<3) x[i]+=7;
    }
    return 0;
}

int a[N][N];
int x[N];
bool tag[N];
struct StrCmp{
    bool operator()(const char *s1,const char *s2)const{
         return strcmp(s1,s2)<0;
    }
};
map<const char *,int,StrCmp> mp;
int main()
{
    //freopen("cin.txt","r",stdin);
    mp["SUN"]=0;
    mp["MON"]=1;
    mp["TUE"]=2;
    mp["WED"]=3;
    mp["THU"]=4;
    mp["FRI"]=5;
    mp["SAT"]=6;
    int n,m;
    while(cin>>n>>m){
        if(n+m==0) break;
        memset(a,0,sizeof(a));
        int k,dex;
        char s1[10],s2[10];
        for(int i=0;i<m;i++){
            scanf("%d%s%s",&k,s1,s2);
            a[i][n]=((mp[s2]-mp[s1]+1)%7+7)%7; // fired day counts
            for(int j=0;j<k;j++){
                scanf("%d",&dex);
                a[i][dex-1]=(a[i][dex-1]+1)%7;
            }
        }
        int r,c;
        Gauss(a,m,n,r,c);
        int g=reback(a,m,n,r,c,x,tag);
        if(g>0) printf("Multiple solutions.\n");
        else if(g<0) printf("Inconsistent data.\n");
        else {
            for(int i=0;i<n-1;i++) printf("%d ",x[i]);
            printf("%d\n",x[n-1]);
        }
    }
    return 0;
}







theArcticOcean
关注
专栏目录
高斯消元练习题集
Miserable_ccf的博客
11-01 1876
1.POJ-1222 EXTENDED LIGHTS OUT 题意 开关问题，给一个矩阵，每次翻转，上下左右都会一起翻转，问你翻转哪些位置，可以把灯全部关上。 思路 异或运算下的高斯消元 典型的开关问题，和 POJ 1830 开关问题 是一样的属于 XOR 类型的消元。 /*************************** *author:ccf *source:poj-1222- EXTE...
高斯消元训练
而濡木染
09-06 447
EXTENDED LIGHTS OUT Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 6443 Accepted: 4229 Description In an extended version of the game Lights Out, is a puzzle with 5
高斯消元基础题
ACdreamer
05-31 4104
题目：EXTENDED LIGHTS OUT   题意：5*6的一个由灯组成的方阵 操作一个灯 周围的上下左右四个灯会发生相应变化 即由灭变亮 由亮变灭 问如何操作使灯全灭。 本题方法：猛戳这里 #include #include #include #include using namespace std; const int N=105; int eq
高斯消元习题
cheng__yu_的博客
07-23 2502
高斯消元P3389 【模板】高斯消元法 （浮点数模板题） P3389 【模板】高斯消元法 （浮点数模板题） 链接：https://www.luogu.com.cn/problem/P3389 #include <bits/stdc++.h> #define ll long long using namespace std; const int maxn=110; int n; double a[maxn][maxn]; double x[maxn]; int guess(double a[]
矩阵高斯消元解方程组
打渔还是晒网 —— stone
05-10 1396
数论 - 高斯消元算法
林深不见鹿 的博客
07-27 2185
1.高斯消元 (1)定义 高斯消元法是求解线性方阵组的一种算法，它也可用来求矩阵的秩，以及求可逆方阵的逆矩阵。它通过逐步消除未知数来将原始线性系统转化为另一个更简单的等价的系统。它的实质是通过初等行变化将线性方程组的增广矩阵转化为行阶梯矩阵. (2)初等行变换 初等行变换 定义：所谓数域P上矩阵的初等行变换是指下列3种变换： (1）以P中一个非零的数乘矩阵的某一行 (2）把矩阵的某一行的c倍加到另一行，这里c是P中的任意一个数 (3）互换矩阵中两行的位置 一般来说，一个矩阵经过初等行变换后就变成了另一个矩
【数学】高斯消元
最新发布
OI博客
05-02 856
讲解基本高斯消元
高斯消元
Nismilesucc的学习札记
08-20 2514
1.高斯消元 在O(n^3)下求解一个包含n个方程组n个未知数的多元线性方程组 n元线性方程组解的三种情况：无穷多组解、唯一解、无解 例题1 解线性方程组 输入 3 1.00 2.00 -1.00 -6.00 2.00 1.00 -3.00 -9.00 -1.00 -1.00 2.00 7.00 输出 1.00 -2.00 3.00 代码 例题2 解异或线性方程组 输入 3 1 1 0 1 0 1 1 0 1 0 0 1 输出 1 0 0 代码 ...
高斯消元入门⑥-配合求解概率/期望dp-HDU4870
qq_35577488的博客
01-02 517
前言: 对于期望dp/概率dp，我们有状态以及状态转移方程.我们将所有的状态看成未知数，状态转移方程看成线性方程。那么就得到了一个线性方程组. 适用情形: 状态转移有后效性时 优化点: 大多数情况下它会是一个稀疏矩阵(因为一个转移方程所牵涉到的状态比较少).所以在 消元的时候，遇到0元素就continue. 这个trick可以将复杂度从三次方优化到近似二次方.(前提是要是一个稀疏矩阵.) 题目一:HDU4780 题目大意: 你有两个数A,BA,BA,B。 每次你会选择min(A,B)min(A,B)min(
NOIP ACM 高斯消元练习+数据题解
04-01
本资料集包含了针对这些竞赛的高斯消元练习及对应的数据题解，旨在帮助参赛者深入理解和熟练应用高斯消元法。 一、高斯消元法概述 高斯消元法由数学家卡尔·弗里德里希·高斯提出，其核心思想是通过一系列行初等...
gaosi.rar_方程组_高斯消元_高斯消元法
09-20
高斯消元法是线性代数中一种重要的求解线性方程组的方法，由19世纪德国...在"gaosi.txt"文件中，可能包含具体示例或练习题目，通过应用上述理论步骤，可以解出方程组的解。这个文件是学习和实践高斯消元法的宝贵资源。
高斯消元入门②-异或线性方程组-POJ1222
qq_35577488的博客
12-31 552
题目大意: 给你一个5∗65 * 65∗6大小的二维01矩阵代表灯的开关状态。改变任意一个灯的状态，会导致其上下左右的灯的状态的改变。问一种灯的开关方案使得最终所有灯都关闭. 题目思路: 首先，一个灯只会被按0次或1次。两次等于没按。不同灯按下的效果是叠加起来的. 解释:想象按下某个位置,等效于原矩阵异或上一个新矩阵。 例如按下(2,2)(2,2)(2,2)产生的效果. A2,2=[010111010]A_{2,2} = \left[\begin{array}{lll} 0 & 1 & 0
algorithm 题集三 (16.05.24)
theArcticOcean
05-24 5492
练习的故事仍然继续下去，今天写下的题目难度适中，记录自己近段时间的训练成果。涉及数学，动态规划，字符串问题。poj 2159 Ancient Cipher(规律)http://poj.org/problem?id=2159 大意：密码加密方法，两种加密方式，第一种，替换加密：Substitution cipher changes all occurrences of each letter to
费马数
theArcticOcean
12-03 3731
形如的数字被称之为费马数，费马数有可能是素数，有可能不是，费马数的素因子都具备这样的形式： 对于所有的非负整数n，有这样的定理：推导：应用数学归纳法，对于n=1，有  现在做一个假设：  有：        由费马数的特殊的表达式可以猜想任意两个不同的费马数一定是互质的。事实上这也是费马数的一个性质，由前面的定理可以知道，，从中随便拿一个数字来和研究，设， 那么有，所以， 因为费马数是奇
快速乘取模
theArcticOcean
12-30 3639
问题来源于：http://acm.nefu.edu.cn/JudgeOnline/problemshow.php?problem_id=609 简单的讲：求解a*b%c 其中0 分析：这样的式子和a^b%c很像，所以可以用类似于二进制取模的方法来做。即，将b写成二进制来看，然后拆开相加：
毕达哥拉斯三元组
theArcticOcean
02-17 2938
毕达哥拉斯三元组：x^2+y^2=z^2 如果有(x,y,z)=1，那么这个这个毕达哥拉斯三元组称为原本的。 原本的毕达哥拉斯三元组满足： x=m^2-n^2 y=2mn z=m^2+n^2 （m,n互素，且m,n奇偶性不同） fzu 1669 Right-angled Triangle http://acm.fzu.edu.cn/problem.php?pid=1669
c(n,m) mod p 1 Lucas 定理
theArcticOcean
09-20 2843
普通的组合数C(n,m)在数据较小的情况下可以先用杨辉三角存储组合值，取模的话再%p即可。但是如果n,m很大，组合的结果自然很多，pascal自然不能完成任务，这样的取模问题可以使用数论里的Lucas定理来解决。 数论Lucas定理是用来求 c(n,m) mod p的值,p是素数（从n取m组合，模上p）。 描述为: Lucas(n,m,p)=Cm(n%p,m%p)* Lucas(n/p,m/
素因子分解 （快速筛法&&试除法）
theArcticOcean
06-15 2776
素因子分解的算法有很多，费马因子分解：比试除法更加高效，是计算机中广泛使用的很多更有效的因子分解算法的基础。二次筛法和数域筛法用于数百位的十进制的大数字。而数字越大数域筛法更好。现在暂时仅仅写了最基础的试除法，更好的算法还等着我去学习~~ #include #include using namespace std; const int N=5e6; //有些特殊的数字素因子会特别大,这里筛选出
数值计算期末复习：牛顿迭代法、高斯消元与插值问题详解
本资源是一份数值计算期末复习的练习题与答案集合，主要涵盖了计算题、填空题以及判断题。以下是对这些题目的详细解析： 一、计算题 1. 题目要求证明某个区间内一给定方程有唯一实根，并利用牛顿迭代法求解，使得...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值