SPOJ (不知道题号是多少) DETER3 - Find The Determinant III

%%%%%%%%%Jin Bin

欧几里得消元

相当好用

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
ll a[205][205],p;
int n;
ll Gauss(){
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++)
	a[i][j]=(a[i][j]%p+p)%p;
	ll ans=1;
	for(int i=1;i<=n;i++){
		for(int j=i+1;j<=n;j++)
		while(a[j][i]){
			ll t=a[i][i]/a[j][i];
			for(int k=i;k<=n;k++)
			(a[i][k]-=a[j][k]*t)%=p;
			for(int k=i;k<=n;k++)
			swap(a[i][k],a[j][k]);
			ans=-ans;
		}
		if(!a[i][i])return 0;
		ans=ans*a[i][i]%p;
	}
	return (ans+p)%p;
}
int main(){
	//freopen("a.in","r",stdin);
	//freopen("a.out","w",stdout);
	while(~scanf("%d%lld",&n,&p)){
		for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
		scanf("%lld",&a[i][j]);
		printf("%lld\n",Gauss());
	}
	return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值