poj 2065 高斯消元

最新推荐文章于 2023-05-23 16:47:36 发布

NeighThorn

最新推荐文章于 2023-05-23 16:47:36 发布

阅读量385

点赞数 1

分类专栏： poj题解模板

本文链接：https://blog.csdn.net/neighthorn/article/details/51387642

版权

poj题解同时被 2 个专栏收录

94 篇文章 0 订阅

订阅专栏

模板

54 篇文章 0 订阅

订阅专栏

萌萌哒的高斯消元模板题
恩，这道题题目很鬼畜，￣へ￣，请各位拿好速效救心丸>o<
首先给出cas**case数目
每个case一行
首先输入p**此题再mod p下进行
然后输入一个字符串,字符串n个字母就有n个未知数
然后就有了n个方程
（a1*1^0+a2*1^1+……+an*1^n-1)%p=str[1]%p
.
.
.
.
.
.
（a1*n^0+a2*n^1+……+an*n^n-1)%p=str[n]%p
(⊙v⊙)嗯
然后你就把gauss模板一摆，就万事大吉啦~(≧▽≦)/~

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
const int maxn=70+10;
int cas,p,cnt,a[maxn],len,coefficient[maxn][maxn],ans[maxn];//每个方程的系数 
char str[maxn];
void gauss(){
    for(int i=1,j=1;i<=cnt&&j<=cnt;i++,j++){//j枚举第几个主元，i枚举第几个方程 
        int k;
        for(k=i;k<=cnt;k++)
            if(coefficient[k][j])
                break; //寻找用第k方程来消去其他方程的第j主元
        for(int s=1;s<=cnt+1;s++)
            swap(coefficient[i][s],coefficient[k][s]);
        for(int k=1;k<=cnt;k++){
            if(k!=i&&coefficient[k][j]){
                int b1=coefficient[k][j],b2=coefficient[i][j];
                for(int s=1;s<=cnt+1;s++)
                    coefficient[k][s]=((coefficient[k][s]*b2-coefficient[i][s]*b1)%p+p)%p;
            }
        }
    } 
    for(int i=cnt;i>0;i--){//第i个方程 倒着枚举可以保证下面枚举ans是都是已经求出的 
        int temp=coefficient[i][cnt+1];
        for(int j=i+1;j<=cnt;j++)
            temp=((temp-coefficient[i][j]*ans[j])%p+p)%p;
        while(temp%coefficient[i][i]!=0)
            temp+=p;
        ans[i]=(temp/coefficient[i][i])%p;
    }
}
int main(){
    scanf("%d",&cas);
    while(cas--){
        cnt=0,scanf("%d",&p),cin>>str,len=strlen(str);
        for(int i=0;i<len;i++){
            if(str[i]=='*')
                a[++cnt]=0;
            else
                a[++cnt]=(str[i]-'a'+1)%p;
            coefficient[cnt][len+1]=a[cnt];
        }
        for(int i=1;i<=cnt;i++)
            for(int j=1;j<=cnt;j++)
                if(j==1)
                    coefficient[i][1]=1;
                else
                    coefficient[i][j]=coefficient[i][j-1]*i%p;
        gauss();
        for(int i=1;i<cnt;i++)
            cout<<ans[i]<<" ";
        cout<<ans[cnt]<<endl;
    }
    return 0;
}