SPOJ-DETER3 Find The Determinant III

求矩阵行列式的模板题

因为需要取模,并且有除法,可以在除法完成之后再取模

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
using namespace std;
typedef long long ll;
const int N=200+10;
ll mod;
ll a[N][N],b[N][N];
ll det(int n)
{
	ll ret=1;
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
			b[i][j]=a[i][j];
	for(int i=0;i<n;i++)
	{
		for(int j=i+1;j<n;j++)
		{
			while(b[j][i])
			{
				ll t=b[i][i]/b[j][i];
				for(int k=i;k<n;k++)
				{
					b[i][k]=(b[i][k]-t*b[j][k]%mod+mod)%mod;
					swap(b[i][k],b[j][k]);
				}
				ret=(-ret+mod)%mod;
			}
		}
		ret=ret*b[i][i]%mod;
	}
	return ret;
}
int main()
{
    int n;
    while(~scanf("%d%lld",&n,&mod))
    {
    	for(int i=0;i<n;i++)
    		for(int j=0;j<n;j++)
    			scanf("%lld",&a[i][j]);
    	printf("%lld\n",det(n));
    }
    return 0;
}

猜你喜欢

转载自blog.csdn.net/wl16wzl/article/details/83573837
今日推荐