高斯消元专题

最新推荐文章于 2018-08-19 14:16:40 发布
置顶 最新推荐文章于 2018-08-19 14:16:40 发布
阅读量3.6k 收藏 1
点赞数 1
分类专栏： 数论 文章标签： ini sun
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qiqijianglu/article/details/7945197
版权

由于老师的讲课，我又重拾高斯消元，做了一天的高斯，感觉老师讲得很好，模板也不错，我懒得自己写，就把同僚的模板拿过来用了，反正基本原理都懂了。

http://poj.org/problem?id=1222

高斯消元模板题，只有唯一解。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;
const int maxn=30;
int a[maxn][maxn+1],x[maxn];
int equ,var,free_num;
void Debug()
{
   for(int i=0;i<equ;i++)
   {
      for(int j=0;j<var+1;j++)
        cout<<a[i][j]<<" ";
      cout<<endl;
   }
}
int gcd(int a,int b)
{
     if(a<0) return gcd(-a,b);
     if(b<0) return gcd(a,-b);
     return b==0?a:gcd(b,a%b);
}
int Gauss()
{
    int k,col=0;
    for(k=0;k<equ&&col<var;k++,col++)
    {
        int mx=k;
        for(int i=k+1;i<equ;i++)
          if(a[i][col]>a[mx][col]) mx=i;
        if(mx!=k)
           for(int i=k;i<var+1;i++)  swap(a[k][i],a[mx][i]);
        if(!a[k][col]) { k--;continue; }
        for(int i=k+1;i<equ;i++)
           if(a[i][col]!=0)
           {
               int lcm=a[k][col]/gcd(a[k][col],a[i][col])*a[i][col];
               int ta=lcm/a[i][col],tb=lcm/a[k][col];
               if(a[i][col]*a[k][col] < 0)  tb = -tb;
               for(int j=col;j<var+1;j++)
                 a[i][j]=((a[i][j]*ta)%2 - (a[k][j]*tb)%2+2)%2;
           }
     }
    // Debug();
     //cout<<col<<endl<<endl;
     for(int i=k;i<equ;i++)
        if(a[i][var]) return -1;
     for(int i=0,j;i<equ;i++)
         if(!a[i][i])
         {
              for(j=i+1;j<var;j++)
                 if(a[i][j]) break;
              if(j>=var) break;
              for(int r=0;r<equ;r++) swap(a[r][i],a[r][j]);
         }

         //唯一解
      for(int i=k-1;i>=0;i--)
      {
          int tmp=a[i][var]%2;
          for(int j=i+1;j<var+1;j++)
             if(a[i][j]) tmp=(tmp-a[i][j]*x[j]%2+2)%2;
          x[i]=(tmp/a[i][i])%2;
      }
}
void init()
{
     memset(a,0,sizeof(a));
     memset(x,0,sizeof(x));
     for(int i=0;i<5;i++)
     for(int j=0;j<6;j++)
      {
           if(i!=0) a[i*6+j][(i-1)*6+j]=1;
           if(i!=4) a[i*6+j][(i+1)*6+j]=1;
           if(j!=0) a[i*6+j][i*6+j-1]=1;
           if(j!=5) a[i*6+j][i*6+j+1]=1;
           a[i*6+j][i*6+j]=1;
      }
}
int main(int argc, char *argv[])
{
    int T,ca=1;
    scanf("%d",&T);
    equ=30,var=30;
    while(T--)
    {
        init();
        for(int i=0;i<30;i++) scanf("%d",&a[i][30]);
        Gauss();
        printf("PUZZLE #%d\n",ca++);
        for(int i=0;i<30;i++)
        {
            if(i%6!=0) putchar(' ');
            printf("%d",x[i]);
            if(i%6==5) puts("");
        }
    }
    return EXIT_SUCCESS;
}

 

fzu1704  http://acm.fzu.edu.cn/problem.php?pid=1704

求解的个数，+高精

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <memory.h>
using namespace std;
const int maxn=110;
const int inf=0x3fffffff;
int a[maxn][maxn+1];
int equ,var,free_num;
int n,m;
inline int gcd(int a,int b)
{
     return b==0?a:gcd(b,a%b);
}
int LCM(int a,int b)
{
    return (a/gcd(a,b))*b;
}
void multi(int x)
{
    int str[40];
    memset(str,0,sizeof(str));
    str[0]=1;
    int c=1;
    for(int i=1;i<=x;i++)
    {
        int g=0;
        for(int j=0;j<c;j++)
        {
            str[j]=str[j]*2+g;
            g=str[j]/10;
            str[j]=str[j]%10;
        }
        if(g!=0) str[c++]=g;
    }
    for(int i=c-1;i>=0;i--)
    cout<<str[i];
    puts("");
}
void Gauss()
{
    int k,col=0;
    for(k=0;k<equ&&col<var;k++,col++)
    {
        int mx=k;
        for(int i=k+1;i<equ;i++)
          if(a[i][col]>a[mx][col]) mx=i;
        if(mx!=k)
           for(int i=k;i<var+1;i++)  swap(a[k][i],a[mx][i]);
        if(a[k][col]==0) { k--;continue; }
        for(int i=k+1;i<equ;i++)
           if(a[i][col]!=0)
           {
               int lcm=LCM(a[i][col],a[k][col]);
               int ta=lcm/a[i][col],tb=lcm/a[k][col];
               if(a[i][col]*a[k][col] < 0)  tb = -tb;
               for(int j=col;j<var+1;j++)
                 a[i][j]=((a[i][j]*ta)%2 - (a[k][j]*tb)%2+2)%2;
           }
     }

     for(int i=k;i<equ;i++)
        if(a[i][col]!=0) {puts("0");return ;}
    multi(var-k);
}
int main()
{
    int t,count=1;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&equ,&var);
        memset(a,0,sizeof(a));
        for(int i=0;i<equ;i++)
        scanf("%d",&a[i][var]);
        for(int i=0;i<var;i++)
        {
            int u,v;
            scanf("%d",&u);
            while(u--)
            {
                scanf("%d",&v);
                v--;
                a[v][i]=1;
            }
        }
        Gauss();
    }
    return EXIT_SUCCESS;
}


hdu3364求多少组解http://acm.hdu.edu.cn/showproblem.php?pid=3364

 

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;
const int maxn=51;
const int inf=0x3fffffff;
int a[maxn][maxn+1],b[maxn][maxn+1];
int equ,var,free_num;
int n,m;
int gcd(int a,int b)
{
     if(a<0) return gcd(-a,b);
     if(b<0) return gcd(a,-b);
     return b==0?a:gcd(b,a%b);
}
void Gauss()
{
    int k,col=0;
    for(k=0;k<equ&&col<var;k++,col++)
    {
        int mx=k;
        for(int i=k+1;i<equ;i++)
          if(a[i][col]>a[mx][col]) mx=i;
        if(mx!=k)
           for(int i=k;i<var+1;i++)  swap(a[k][i],a[mx][i]);
        if(!a[k][col]) { k--;continue; }
        for(int i=k+1;i<equ;i++)
           if(a[i][col]!=0)
           {
               int lcm=a[k][col]/gcd(a[k][col],a[i][col])*a[i][col];
               int ta=lcm/a[i][col],tb=lcm/a[k][col];
               if(a[i][col]*a[k][col] < 0)  tb = -tb;
               for(int j=col;j<var+1;j++)
                 a[i][j]=((a[i][j]*ta)%2 - (a[k][j]*tb)%2+2)%2;
           }
     }

     for(int i=k;i<equ;i++)
        if(a[i][col]) {puts("0");return ;}
    printf("%I64d\n",1LL<<(var-k));

}
int main()
{
    int t,count=1;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&equ,&var);
        memset(b,0,sizeof(b));
        for(int i=0;i<var;i++)
        {
            int q;
            scanf("%d",&q);
            while(q--)
            {
                int num;
                scanf("%d",&num);
                num--;
                b[num][i]=1;
            }
        }
        printf("Case %d:\n",count++);
        int q;
        scanf("%d",&q);
        while(q--)
        {
         memcpy(a,b,sizeof(b));
         for(int i=0;i<equ;i++)
         scanf("%d",&a[i][var]);
         Gauss();
        }
    }
    //system("pause");
    return EXIT_SUCCESS;
}


 


http://poj.org/problem?id=1681

动最少的砖，枚举自由元，找1的个数最少的那组解。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <stdio.h>
using namespace std;
const int inf=0x3fffffff;

const int maxn=15*15;
int a[maxn][maxn+1],x[maxn],equ,var,n;
 char s[20][20];

void init()
{
    equ=var=n*n;
    for(int i=0;i<n;i++)
       scanf("%s",s[i]);
    memset(a,0,sizeof(a));
    for(int i=0;i<var;i++) a[i][var]= (s[i/n][i%n]!='y'),a[i][i]=1;
    for(int i=0;i<n;i++)
    for(int j=0;j<n;j++)
    {
       /* if(i!=0) a[(i-1)*n+j][i*n+j]=1;
        if(i!=n-1) a[(i+1)*n+j][i*n+j]=1;
        if(j!=0) a[i*n+j-1][i*n+j]=1;
        if(j!=n-1) a[i*n+j+1][i*n+j]=1;*/
             if(i-1>=0) a[i*n+j][(i-1)*n+j]=1;
             if(i+1<n) a[i*n+j][(i+1)*n+j]=1;
             if(j-1>=0) a[i*n+j][i*n+j-1]=1;
             if(j+1<n) a[i*n+j][i*n+j+1]=1;
    }
}
void debug()
{
    for(int i=0;i<equ;i++)
    {
        for(int j=0;j<=var;j++)
          cout<<a[i][j]<<" ";
        cout<<endl;
    }
    cout<<endl;
}
int gcd(int a,int b)
{
    if(a<0) return gcd(-a,b);
    if(b<0) return gcd(a,-b);
    return b==0?a:gcd(b,a%b);
}

int gauss()
{
    int k,col;
    for(k=col=0;k<equ&&col<var;k++,col++)
    {
        int mx=k;
        for(int i=k+1;i<equ;i++)
           if(a[i][col]>a[mx][col]) mx=i;
        if(k!=mx){
            for(int i=k;i<var+1;i++)
               swap(a[k][i],a[mx][i]);
        }
        if(!a[k][col]){
            k--;continue;
        }
        for(int i=k+1;i<equ;i++)
        if(a[i][col])
        {
            int lcm=a[k][col]/gcd(a[k][col],a[i][col])*a[i][col];
            int ta=lcm/a[i][col],tb=lcm/a[k][col];
            for(int j=col;j<var+1;j++)
               a[i][j]=((a[i][j]*ta-a[k][j]*tb)%2+2)%2;
        }
    }
   // debug();
    for(int i=k;i<equ;i++)
      if(a[i][col]) return inf;

    for(int i=0,j;i<equ;i++)
       if(!a[i][i])
       {
           for(j=i+1;j<var;j++)
             if(a[i][j]) break;
           if(j>=var) break;
           for(int r=0;r<equ;r++) swap(a[r][i],a[r][j]);
       }

 /*   int ret=inf,lim= (1<<(var-k));

   for(int j=0;j<lim;j++)
   {
     for(int i=0;i<(var-k);i++)
       if(j&(1<<i)) x[equ-1-i]=1;
       else x[equ-1-i]=0;
     for(int i=k-1;i>=0;i--)
     {
        int tmp=a[i][var];
        for(int j=i+1;j<var;j++)
          tmp=((tmp-a[i][j]*x[j])%2+2)%2;
        x[i]=tmp/a[i][i]%2;
     }
     int cnt=0;
     for(int i=0;i<var;i++)
        if(x[i]) cnt++;
     ret=min(cnt,ret);
     if(ret==0) break;
   }*/     // 二进制枚举
     int Min=inf;
     for(int j=0;j<(1<<(equ-k));j++)
     {
        int tmp=j,p=equ-1;
        while(tmp) x[p--]=tmp%2,tmp>>=1;
        for(int i=k-1;i>=0;i--)
        {
          int tmp=a[i][var]%2;
          for(int j=i+1;j<var;j++)
             if(a[i][j]) tmp=(tmp-a[i][j]*x[j]%2+2)%2;
          x[i]=(tmp/a[i][i])%2;
        }
        tmp=0;
        for(int i=0;i<var;i++) tmp+=x[i];
        Min=min(Min,tmp);
        if(Min==0)
        break;
     }   //全部枚举
   return Min;
}
int main()
{
    int T;

    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        init();
       // debug();
        int re=gauss();
        //cout<<re<<endl;
        if(re==inf) puts("inf");
        else printf("%d\n",re);
    }

    return 0;
}


http://poj.org/problem?id=1753

同上，不贴代码了

http://poj.org/problem?id=3185

没什么特别的

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=3858

枚举四种图案，没什么特别的

http://acm.hdu.edu.cn/showproblem.php?pid=4200

连动题

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;
const int maxn=101;
const int inf=0x3fffffff;
int a[maxn][maxn+1],x[maxn];
int equ,var,free_num;
int n,m;
int gcd(int a,int b)
{
     if(a<0) return gcd(-a,b);
     if(b<0) return gcd(a,-b);
     return b==0?a:gcd(b,a%b);
}
int Gauss()
{
    int k,col=0;
    for(k=0;k<equ&&col<var;k++,col++)
    {
        int mx=k;
        for(int i=k+1;i<equ;i++)
          if(a[i][col]>a[mx][col]) mx=i;
        if(mx!=k)
           for(int i=k;i<var+1;i++)  swap(a[k][i],a[mx][i]);
        if(!a[k][col]) { k--;continue; }
        for(int i=k+1;i<equ;i++)
           if(a[i][col]!=0)
           {
               int lcm=a[k][col]/gcd(a[k][col],a[i][col])*a[i][col];
               int ta=lcm/a[i][col],tb=lcm/a[k][col];
               if(a[i][col]*a[k][col] < 0)  tb = -tb;
               for(int j=col;j<var+1;j++)
                 a[i][j]=((a[i][j]*ta)%2 - (a[k][j]*tb)%2+2)%2;
           }
     }

     for(int i=k;i<equ;i++)
        if(a[i][var]) return inf;
     for(int i=0,j;i<equ;i++)
         if(!a[i][i])
         {
              for(j=i+1;j<var;j++)
                 if(a[i][j]) break;
              if(var==j) break;
              for(int r=0;r<equ;r++) swap(a[r][i],a[r][j]);
         }

     int ret=inf;
     long long lim= (1<<(var-k));

   for(long long  j=0;j<lim;j++)
   {
     for(int i=0;i<(var-k);i++)
       if(j&(1<<i)) x[equ-1-i]=1;
       else x[equ-1-i]=0;
     for(int i=k-1;i>=0;i--)
     {
        int tmp=a[i][var];
        for(int j=i+1;j<var;j++)
          tmp=((tmp-a[i][j]*x[j])%2+2)%2;
        x[i]=tmp/a[i][i]%2;
     }
     int cnt=0;
     for(int i=0;i<var;i++)
        if(x[i]) cnt++;
     ret=min(cnt,ret);
     if(ret==0) break;
   }

     return ret;
}
void  init()
{
     memset(a,0,sizeof(a));
     memset(x,0,sizeof(x));
     for(int i=0;i<n;i++)
     {
         scanf("%d",&a[i][var]);
         a[i][i]=1;
     }
     for(int i=0;i<n;i++)
     {
         int p,q;
        if(i-m>=0) p=i-m;else p=0;
        if(i+m<n) q=i+m;else q=n-1;
        for(int j=p;j<=q;j++)
           a[j][i]=1;
     }
     return ;
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&n,&m);
        //cout<<"*"<<endl;
        equ=n;var=n;
        init();
        int ans=Gauss();
        if(ans==inf) puts("impossible");
        else printf("%d\n",ans);
    }
    return 0;
}


http://poj.org/problem?id=2947

同余题

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;
const int maxn=301;
const int inf=0x3fffffff;
int a[maxn][maxn+1],x[maxn];
int equ,var,free_num;
char str[8][4]={"MON","TUE","WED","THU","FRI","SAT","SUN"};
__int64 exGcd(__int64 a,__int64 b, __int64 &xx, __int64 &yy)
{
       if(b == 0)
       {
       xx = 1;
       yy = 0;
       return a;
       }
       __int64 r = exGcd(b, a % b, xx, yy);
       __int64 t = xx;
       xx = yy;
       yy = t - a / b * yy;
    return r;
}
int gcd(int a,int b)
{
     if(a<0) return gcd(-a,b);
     if(b<0) return gcd(a,-b);
     return b==0?a:gcd(b,a%b);
}
void Gauss()
{

    int k,col=0;
    __int64 xx,yy;
    for(k=0;k<equ&&col<var;k++,col++)
    {
        int mx=k;
        for(int i=k+1;i<equ;i++)
          if(abs(a[i][col])>abs(a[mx][col])) mx=i;
        if(mx!=k)
           for(int i=k;i<var+1;i++)  swap(a[k][i],a[mx][i]);
        if(!a[k][col]) { k--;continue; }
        for(int i=k+1;i<equ;i++)
           if(a[i][col]!=0)
           {
               int lcm=a[k][col]/gcd(a[k][col],a[i][col])*a[i][col];
               int ta=lcm/a[i][col],tb=lcm/a[k][col];
               if(a[i][col]*a[k][col] < 0)  tb = -tb;
               for(int j=col;j<var+1;j++)
                 a[i][j]=((a[i][j]*ta)%7 - (a[k][j]*tb)%7+7)%7;
           }
     }
     for(int i=k;i<equ;i++)
        if(a[i][col]!=0) {puts("Inconsistent data.");return ;}
        if(var>k){puts("Multiple solutions.");return ;}
    for(int i=0,j;i<equ;i++)
         if(a[i][i]==0)
         {
              for(j=i+1;j<var;j++)
                 if(a[i][j]) break;
              //if(var==j) break;
              if(j>=var) break;
              for(int r=0;r<equ;r++) swap(a[r][i],a[r][j]);
         }

      for(int i=k-1;i>=0;i--)
      {
          int tmp=a[i][var]%7;
          for(int j=i+1;j<var;j++)
          if(a[i][j]) tmp=(tmp-a[i][j]*x[j]%7+7)%7;
       __int64 d=exGcd(a[i][i],7,xx,yy);
       if(tmp%d!=0){puts("Inconsistent data.");return ;}
       else{
              x[i]=tmp/d*xx%7+7;
              x[i]=x[i]%(7/d);
            }
            while(x[i]<3) x[i]=x[i]+7;
      }

      for(int i=0;i<var;i++)
      {
          if(i==0) cout<<x[i];
          else cout<<" "<<x[i];
      }
      puts("");

}


void  init()
{
     memset(a,0,sizeof(a));
     memset(x,0,sizeof(x));
     int d;
     char s1[4],s2[4];
     for(int i=0;i<equ;i++)
     {
         scanf("%d%s%s",&d,s1,s2);
         int q1,q2;
         for(int j=0;j<8;j++)
         if(!strcmp(s1,str[j])){q1=j;break;}
         for(int j=0;j<8;j++)
         if(!strcmp(s2,str[j])){q2=j;break;}
         while(d--)
         {
             int num;
             scanf("%d",&num); num--;a[i][num]++;
             a[i][num]%=7;
         }
         a[i][var]=(q2-q1+1+7)%7;
     }
     return ;
}
int main()
{
    while(scanf("%d%d",&var,&equ)!=EOF)
    {
        if(var==0&&equ==0) break;
        init();
        Gauss();
    }
    return 0;
}


http://poj.org/problem?id=2065

题目大意：给你一个素数P(P<=30000)和一串长为n的字符串str[]。字母'*'代表0，字母a-z分别代表1-26，这n个字符所代表的数字分别代表f(1)、f(2)....f(n)。
定义: f (k) = ∑0<=i<=n-1aiki (mod p) (1<=k<=n,0<=ai<P)
求a0、a1.....an-1。题目保证肯定有唯一解

解题思路：高斯消元。根据上面的公式显然可以列出有n个未知数的n个方程式：
a0*1^0 + a1*1^1+a2*1^2+........+an-1*1^(n-1) = f(1)
a0*2^0 + a1*2^1+a2*2^2+........+an-1*2^(n-1) = f(2)
..............
a0*n^0 + a1*n^1+a2*n^2+........+an-1*n^(n-1) = f(n)

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <math.h>
using namespace std;
const int maxn=71;
const int inf=0x3fffffff;
int a[maxn][maxn+1],x[maxn];
int equ,var,p;
__int64 exGcd(__int64 a,__int64 b, __int64 &xx, __int64 &yy)
{
       if(b == 0)
       {
       xx = 1;
       yy = 0;
       return a;
       }
       __int64 r = exGcd(b, a % b, xx, yy);
       __int64 t = xx;
       xx = yy;
       yy = t - a / b * yy;
    return r;
}
inline int gcd(int a,int b)
{
     return b==0?a:gcd(b,a%b);
}
int LCM(int a,int b)
{
    return (a/gcd(a,b))*b;
}
void Gauss()
{

    int k,col=0;
    __int64 xx,yy;
    for(k=0;k<equ&&col<var;k++,col++)
    {
        int mx=k;
        for(int i=k+1;i<equ;i++)
          if(abs(a[i][col])>abs(a[mx][col])) mx=i;
        if(mx!=k)
           for(int i=k;i<var+1;i++)  swap(a[k][i],a[mx][i]);
        if(a[k][col]==0) { k--;continue; }
        for(int i=k+1;i<equ;i++)
           if(a[i][col]!=0)
           {
               int lcm=LCM(a[i][col],a[k][col]);
               int ta=lcm/a[i][col],tb=lcm/a[k][col];
               if(a[i][col]*a[k][col] < 0)  tb = -tb;
               for(int j=col;j<var+1;j++)
                 a[i][j]=((a[i][j]*ta)%p - (a[k][j]*tb)%p+p)%p;
           }
     }

      for(int i=0,j;i<equ;i++)
         if(a[i][i]==0)
         {
              for(j=i+1;j<var;j++)
                 if(a[i][j]!=0) break;
              if(j>=var) break;
              for(int r=0;r<equ;r++) swap(a[r][i],a[r][j]);
         }

      for(int i=k-1;i>=0;i--)
      {
          int tmp=a[i][var]%p;
          for(int j=i+1;j<var;j++)
          if(a[i][j]) tmp=(tmp-a[i][j]*x[j]%p+p)%p;
          __int64 d=exGcd(a[i][i],p,xx,yy);

              x[i]=tmp/d*xx%p+p;
              x[i]=x[i]%(p/d);
      }

      for(int i=0;i<var;i++)
      {
          if(i==0) cout<<x[i];
          else cout<<" "<<x[i];
      }
      puts("");

}

char s[71];
void  init()
{
     memset(a,0,sizeof(a));
     memset(x,0,sizeof(x));
     scanf("%s",s);
     equ=var=strlen(s);
     for(int i=0;i<equ;i++)
     {
         if(s[i]=='*') a[i][var]=0;
         else a[i][var]=s[i]-'a'+1;
         a[i][var]%=p;
     }
     for(int i=0;i<equ;i++)
     {
         int tt=1;
         for(int j=0;j<var;j++)
       {
         a[i][j]=tt;
         tt*=(i+1);tt%=p;

       }
     }
     return ;
}

int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&p);
        init();
        Gauss();
    }
    return 0;
}


 

qiqijianglu
关注
  • 1
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
高斯消元法专题
10-13 832
首先贡献一道求方程组解的高斯消元法： 1.判断方程组解的问题： 无解：系数矩阵的秩 有唯一解：系数矩阵的秩==增广矩阵的秩，且等于n（方程组数） 有无穷解：系数矩阵的秩==增广矩阵的秩，且小于n http://www.bnuoj.com/bnuoj/problem_show.php?pid=11952 代码：http://blog.csdn.net/qq415200973/arti
（ 数论专题 ）【 高斯消元 】
石旭先森的博客
08-07 296
（ 数论专题 ）【 高斯消元 】 作用：解一元多次方程。 例如：x+y=2， 2x+3y=5 . 解得x=1，y=1 这里介绍的是高斯-约旦消元法。 相对于传统的高斯消元，约旦消元法的精度更好、代码更简单，没有回带的过程。 约旦消元法大致思路如下： 1.选择一个尚未被选过的未知数作为主元，选择一个包含这个主元的方程。 2.将这个方程主元的系数化为1。 3.通过加减消元，消掉其它方程中的这个未知数。 4.重复以上步骤，直到把所有式子变成形如: a1+b0+c*0……=d 我们用矩阵表示每.
高斯消元
NYIST_TC_LYQ的博客
08-09 463
高斯消元  上周打完多校队内开了个会商量了一下每个人的学习计划和接下来的比赛状态。于是这个专题开始了，还有几道题待做，先把题解写一下吧。   我是看大白书学的，挺好理解的，不过前提是要学过线性代数，这样理解起来比较容易，如果没学过可能会有点困难。鄙人愚见：高斯消元类似小学学过的解二元方程组，那么方程组至少要有两个以上才会有解，这很好理解的吧。但求解多元方程组就不那么容易了，而线性代数里面的求解
矩阵乘+高斯消元专题小结【高斯消元】
Senbonzakura
05-20 1090
http://acm.hust.edu.cn/vjudge/contest/view.action?cid=20851#overview I-P 其它题太水了就不写了，这里说一下O和P。 【我太懒不想写】 OP这两个题是高斯消元求期望的题。 对于P题，求从@点到$点的步数期望，先从@点BFS一遍找到所有可以到达的点，对于可达点每个点是一个变量。 因为终点有任意多而起点只有一个，这里可以
高斯消元（非常详细） 高斯消元（非常详细）
最新发布
06-19
高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元...
gaosixiaoyuanfa.rar_高斯消元
09-23
这个名为“gaosixiaoyuanfa.rar_高斯消元”的压缩包文件包含了一份详细的高斯消元法课程资料，通常以PPT的形式呈现，可以帮助学习者深入理解并掌握这一算法。 高斯消元法的基本思想是通过一系列行初等变换，将系数...
高斯消元附解方程神器
04-09
先输入一个数n表示方程元数。 然后再输入一个 n*(n+1)的矩阵 具体见我的博客
gauss_seq.c.zip_高斯消元
09-23
此外，高斯-约旦消元（Gauss-Jordan elimination）是高斯消元的一种扩展，它不仅将矩阵转化为上三角形，还能进一步将其化为单位矩阵，从而直接得到解。 总结来说，`gauss_seq.c`是一个用C语言实现的高斯消元法串行...
高斯消元模版深度认识完整版
06-19
高斯消元模版
高斯消元模板
ACM之梦
05-05 1997
mod 2 方程写法： int Gauss(int n, int m){ int r = 0, c = 0; while(r < n && c < m){ int id = r; for(int i = r + 1; i < n; i++) if(maze[i][c] > maze[id][c])
poj1830(高斯消元解mod2方程组)
weixin_30418341的博客
09-24 155
题目链接：http://poj.org/problem?id=1830 题意：中文题诶～ 思路：高斯消元解 mod2 方程组 有 n 个变元，根据给出的条件列 n 个方程组，初始状态和终止状态不同的位置对应的方程右边常数项为１，状态相同的位置对于的方程组右边的常数项为０．然后用高斯消元解一下即可．若有唯一解输出１即可，要是存在 k 个变元，则答案为 1 << k, 因...
期望与概率DP专题/高斯消元
Qianyri的博客
08-19 1360
高斯消元 总结：概率和期望问题的求解都要依靠递推式，概率问题按照递推式进行dp，期望问题根据递推式化简的难易程度通过化简再dp或高斯消元求解，该类问题大部分仍属于各个类型的dp求解问题 2018.9.3-概率与期望dp学习 Discovering Gold 一排1到n的格子，每个格子上有黄金 ai ,你最开始在 1 号,每一次投骰子决定到哪一个格子，超出1~n范围则重新投掷，你到了哪个格子...
高斯消元解方程组 Gauss-Jordan elimination
willsongrui的专栏
05-02 1250
#include #include #include #include #include using namespace std; const int MAXN=50; int a[MAXN][MAXN];//增广矩阵 int x[MAXN];//解集 bool free_x[MAXN];//标记是否是不确定的变元 /* void Debug(void) { int i, j
poj 2947 高斯消元求解方程组取模
渐近自由
08-04 594
题意： (mat[1][1]*x[1] + mat[1][2]*x[2] + … + mat[1][n]*x[n])%7 =mat[1][n+1] (mat[2][1]*x[1] + mat[2][2]*x[2] + … + mat[2][n]*x[n])%7 =mat[2][n+1] … … (mat[m][1]*x[1] + mat[m][2]*x[2] +
扩展欧几里得&&中国剩余定理 解线性同余方程组 专题
Shirley
10-04 4183
poj1061  扩展欧几里得解方程 #include #include #include using namespace std; __int64 exGcd(__int64 a,__int64 b, __int64 &xx, __int64 &yy); int main() { __int64 x,y,m,n,L,xx,yy; scanf("%I64d%I64d%
hdu4405 Aeroplane chess---金华网络赛 期望
Shirley
09-22 3752
Aeroplane chess Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 24    Accepted Submission(s): 20 Problem Description Hzz loves aeroplan
poj3270 Cow Sorting-------置换群
Shirley
08-31 3118
Cow Sorting Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 4746   Accepted: 1730 Description Farmer John's N (1 ≤ N ≤ 10,000) cows are lined up to be milked
matlab 高斯消元
11-12
高斯消元是一种线性代数中常用的求解线性方程组的方法，Matlab中也提供了相应的函数来实现高斯消元。下面是一个简单的例子： 假设有如下线性方程组： ``` 2x + y - z = 8 -3x - y + 2z = -11 -2x + y + 2z = -3 ``...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值