#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 ;

}