模意义下的高斯消元

求解线性方程组是一个基础的问题,但如果是在模意义下,事情就没有那么简单了。现在给出一个方程组和模数p(是一个质数),如何解出这个方程组呢(保证有唯一一组整数解)?首先,因为我们使用了整形数组来存储数据,计算中就不允许出现小数点,我们可以把两行各乘上一个数后(要保证主元被消去)再进行相减,而且相减的过程要在模意义下进行。注意,设相减的结果为x,如果只是令x=x\ mod\ p,是存在问题的,因为x可能是一个负数。那么令x=(x+p)\ mod\ p是否可行呢?如果x是一个绝对值很大的负数,加上p之后还是负数,该怎么处理呢?我们可以先把x缩小到(-p,p)的范围。所以,我们要令x=(x\ mod\ p+p)\ mod\ p。接下来考虑回代的问题。回代的过程显然要用到除法操作,但模意义下是不能使用除法的,于是,我们可以使用乘法逆元。因为保证了模数p是一个质数,所以我们可以根据费马小定理得知:x^{p-2}即为x在模p意义下的逆元。代码如下。

#include<cstdio>
#define maxn 110
#define r register
using namespace std;
typedef long long ll;
int n,p,maxi;
ll tmp,ans[maxn],a[maxn][maxn];
int read()
{
	r char ch=getchar();r int in=0;
	while(ch>'9'||ch<'0') ch=getchar();
	while(ch>='0'&&ch<='9') in=(in<<3)+(in<<1)+ch-'0',ch=getchar();
	return in;
}
ll ksm(r ll x,r int y)
{
	if(!y) return 1;
	r ll ret=ksm(x,y>>1);
	if(y&1) return ret*ret%p*x%p;
	return ret*ret%p;
}
int main()
{
	n=read(),p=read();
	for(r int i=1;i<=n;i++)
		for(r int j=1;j<=n+1;j++)
			a[i][j]=read();
	for(r int i=1;i<=n;i++)
	{
		if(!a[i][i])//主元不能为0
		{
			maxi=0;
			for(r int j=i+1;j<=n&&!maxi;j++)
				if(a[j][i]) maxi=j;
			if(!maxi) continue;//如果一整列都为0,不需要消元
			for(r int j=i;j<=n+1;j++)
				tmp=a[maxi][j],a[maxi][j]=a[i][j],a[i][j]=tmp;
		}
		for(r int j=i+1;j<=n;j++)
		{
			tmp=a[j][i];
			if(!tmp) continue;//已经为0,不需要消元
			for(r int k=i;k<=n+1;k++)
				a[j][k]=((a[j][k]*a[i][i]-a[i][k]*tmp)%p+p)%p;
		}
	}
	for(r int i=n;i;i--)
	{
		for(r int j=i+1;j<=n;j++)
			a[i][n+1]=((a[i][n+1]-ans[j]*a[i][j])%p+p)%p;
		ans[i]=a[i][n+1]*ksm(a[i][i],p-2)%p;
	}
	for(r int i=1;i<=n;i++) printf("%lld ",ans[i]);
	return 0;
}

 

  • 3
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值