求矩阵行列式的模板题
因为需要取模,并且有除法,可以在除法完成之后再取模
#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;
}