关闭

高斯消元

107人阅读 评论(0) 收藏 举报
#include<stdio.h>
#include<algorithm>
#include<iostream>
#include<string.h>
#include<math.h>
using namespace std;
int a[72][72],x[72],r,c,mod;
int gcd(int a,int b)
{
    return b?gcd(b,a%b):a;
}
int lcm(int a,int b)
{
    return a/gcd(a,b)*b;
}
int pow_mod(int a,int b)
{
    int res=1;
    while(b)
    {
        if(b&1) res=res*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return res;
}
void Gauss()
{
    int k,col,max_r,i,j,LCM,ta,tb,temp;
    for(k=col=0;k<r&&col<c;col++)  //枚举行,列从0开始
    {
        max_r=k;
        for(i=k+1;i<r;i++)
            if(abs(a[i][col])>abs(a[max_r][col])) max_r=i;   //找到这一列里面绝对值最大的
        if(!a[max_r][col]) continue;   //最大的都是0,那么这一列不用考虑了
        if(max_r!=k)
            for(i=col;i<=c;i++) swap(a[max_r][i],a[k][i]);  //否则交换
        for(i=k+1;i<r;i++)
            if(a[i][col]!=0)
            {
                LCM=lcm(abs(a[i][col]),abs(a[k][col]));
                ta=LCM/abs(a[i][col]);
                tb=LCM/abs(a[k][col]);
                if(a[i][col]*a[k][col]<0) tb=-tb;
                for(j=col;j<=c;j++) a[i][j]=((a[i][j]*ta%mod-a[k][j]*tb%mod)+mod)%mod;
            }
        k++;
    }
    //for(i=k;i<r;i++) if(a[i][col]) return -1;
    for(i=c-1;i>=0;i--)
    {
        temp=a[i][col];
        for(j=i+1;j<c;j++)
            if(a[i][j]!=0) temp=((temp-a[i][j]*x[j]%mod)+mod)%mod;
        x[i]=temp*pow_mod(a[i][i],mod-2)%mod;
    }
}
int main()
{
    int i,j,t,v;
    scanf("%d",&t);
    char s[72];
    while(t--)
    {
        memset(a,0,sizeof(a));
        scanf("%d%s",&mod,s);
        r=c=strlen(s);
        for(i=0;i<r;i++) a[i][c]=(s[i]=='*'?0:s[i]-96);
        for(i=0;i<r;i++)
        {
            v=a[i][0]=1;
            for(j=1;j<c;j++)
            {
                v=v*(i+1)%mod;
                a[i][j]=v;
            }
        }
        Gauss();
        for(i=0;i<r;i++)
            printf("%d%c",x[i],i+1==r?'\n':' ');
    }
}









// 模2时的高斯消元
#include<stdio.h>
#include<algorithm>
#include<iostream>
#include<string.h>
using namespace std;
int a[35][35],b[35],r,c,mod=2;
int Gauss()
{
    int k,col,max_r,i,j;
    for(k=col=0;k<r&&col<c;col++)  //枚举行,列从0开始
    {
        max_r=k;
        for(i=k+1;i<r;i++)
            if(abs(a[i][col])>abs(a[max_r][col]) )max_r=i;   //找到这一列里面绝对值最大的
        if(!a[max_r][col]) continue;   //最大的都是0,那么这一列不用考虑了
        if(max_r!=k)
            for(i=col;i<=c;i++) swap(a[max_r][i],a[k][i]);  //否则交换
        for(i=k+1;i<r;i++)
            if(a[i][col]!=0)
                for(j=col;j<=c;j++) a[i][j]^=a[k][j];
        k++;
    }
    for(i=k;i<r;i++)
        if(a[i][col]) return 0;
    if(k<c) return 1<<(c-k);
    return 1;
}
int main()
{
    int i,t,v,u;
    scanf("%d",&t);
    while(t--)
    {
        memset(a,0,sizeof(a));
        scanf("%d",&r);  c=r;
        for(i=0;i<r;i++) scanf("%d",&b[i]);
        for(i=0;i<r;i++)
        {
            scanf("%d",&v);
            a[i][c]=b[i]^v;
            a[i][i]=1;
        }
        while(scanf("%d%d",&u,&v)&&u) a[v-1][u-1]=1;
        v=Gauss();
        if(v) printf("%d\n",v);
        else printf("Oh,it's impossible~!!\n");
    }
}

0
0

查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:6341次
    • 积分:552
    • 等级:
    • 排名:千里之外
    • 原创:50篇
    • 转载:4篇
    • 译文:0篇
    • 评论:0条