The Evaluation of Determinant(求行列式mod一个数的值）

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<string>
#define LL long long
using namespace std;
LL n,m,A[][],p[],pos,d[],r[],len,B[][];
bool vd[]={};
void prime()
{
    pos=;
    for(int i=;i<;i++)
    {
        if(!vd[i])
        {
            if(i>) p[pos++]=i;
            for(int j=(i<<);j<;j+=i)
                vd[i]=;
        }
    }
}
void deal(LL k)
{
    len=;
    for(int i=;i<pos&&k!=;i++)
    {
        if(k%p[i]==)
        {
            while(k%p[i]==)
            {
                d[len++]=p[i];
                k/=p[i];
            }
        }
    }
}
LL exp(LL a,LL b,LL mod)
{
    LL ans=;
    while(b)
    {
        if(b&)ans=ans*a%mod;
        a=a*a%mod;b>>=;
    }
    return ans;
}
void ex_gcd(LL a,LL b,LL &dd,LL &x,LL &y)
{
    if(b==)
        x=,y=,dd=a;
    else
    {
        ex_gcd(b,a%b,dd,y,x);
        y-=x*(a/b);
    }
}
LL gauss(LL mod)
{
    bool flag=;
    LL ans=;
    for(int i=;i<n;i++)
        for(int j=;j<n;j++)
            B[i][j]=A[i][j];
    for(int k=;k<n-;k++)
    {
        LL max_b=B[k][k];int bin=k;
        for(int i=k+;i<n;i++)
            if(B[i][k]>max_b)
                max_b=B[i][k],bin=i;
        if(bin!=k)
        {
            for(int i=k;i<n;i++)
                swap(B[bin][i],B[k][i]);
            flag^=;
        }
        if(B[k][k]<)B[k][k]+=mod;
        LL Ni,y,dd;
        ex_gcd(B[k][k],mod,dd,Ni,y);
        Ni%=mod;
        if(Ni<)Ni+=mod;
        for(int j=k+;j<n;j++)
        {
            B[k][j]=B[k][j]*Ni%mod;
            if(B[k][j]<)B[k][j]+=mod;
            for(int i=k+;i<n;i++)
            {
                B[i][j]=(B[i][j]-(B[k][j]*B[i][k])%mod)%mod;
                if(B[i][j]<)B[i][j]+=mod;
            }
        }
        ans*=B[k][k];
        ans%=mod;
        if(ans<)ans+=mod;
    }
    ans*=B[n-][n-];
    ans%=mod;
    if(flag)ans=-ans;
    if(ans<)ans+=mod;
    return ans;
}

LL china_remain()
{
    LL a,b,c,c1,c2,x,y,dd,N;
    a=d[],c1=r[];
    if(c1==)c1=d[];
    for(int i=;i<len;i++)
    {
        b=d[i],c2=r[i];
        ex_gcd(a,b,dd,x,y);
        c=c2-c1;
        LL b1=b/dd;
        x=((c/dd*x)%b1+b1)%b1;
        c1=a*x+c1;
        a=a*b1;
    }
    return c1%m;
}
int main()
{
    prime();
    while(cin>>n>>m)
    {
        deal(m);
        for(int i=;i<n;i++)
            for(int j=;j<n;j++)
                cin>>A[i][j];
        if(m==)
        {
            cout<<<<endl;
            continue;
        }
        for(int i=;i<len;i++)
        {
            r[i]=gauss(d[i]);
        }
        cout<<china_remain()<<endl;
    }
    return ;
}