高斯消元模板

int gauss(int r, int c)
{
bool flag=false;
int coe=1;
int i=0,t=0;
for(int j=0; j<c; ++j)
{
int index=i;
for(int k=i; k<r; ++k)
if(a[k][j]>0)
{
index=k;
break;
}
if(a[index][j])
{
if(index != i)
{
for(int k=j; k<c; ++k)
swap(a[i][k],a[index][k]);
flag = !flag;
}
for(int k=i+1; k<r; ++k)
if(a[k][j])
{
coe=(coe*a[i][j])%MOD;
++ t;
for(int l=c-1; l>=j; --l)
{
a[k][l]=(a[k][l]*a[i][j]-a[i][l]*a[k][j])%MOD;
if(a[k][l]<0)
a[k][l]+=MOD;
}
}
++i;
}
}
for(i=1; i<MOD; ++i)
if((coe*i)%MOD==1)
break;
int result=i;
for(i=0; i<r; ++i)
result=(result*a[i][i])%MOD;
if(flag)
result =MOD-result;
return result;
}