HZOJ visit

对于前30%的数据，可以考虑dp，f[i][j][k]表示时间为i，在i，j位置的方案数，枚举转移即可。要注意的是可以走到矩阵外。

对于另外30%数据，考虑推一下式子，设向右走y步，左z，上s，下x。那么y-z=n,s-x=m。所以我们枚举s就可以求得sxzy，步数确定之后就比较简单了，显然答案为

$∑C_T^s*C_{T-s}^x*C_{T-s-x}^z*C_{T-s-x-z}^y$,化减得ans=∑$\frac{T!}{s!x!z!y!}$,由于mod是质数，直接逆元干就行了。

对于100%数据，难点就在于p不是质数，开始我想着分解质因数，然后T了。

考虑一下如何求组合数%合数？（合数为若干质数乘积）

将合数质因数分解，用卢卡斯求出组合数mod每个质数，然后发现其实就是线性同余方程组，CRT解即可。

 #include<iostream>
 #include<cstring>
 #include<cstdio>
 #include<queue>
 #define int LL
 #define ma(x) memset(x,0,sizeof(x))
 #define LL long long
 using namespace std;
 int T,mod,n,m;
 int s,x,z,y;
 LL f[][][];
 LL jc[];
 int cnt[];
 int prime[],num;
 bool isprime[];
 LL exgcd(int a,int b,int &x,int &y)
 {
     if(!b){x=,y=;return a;}
     int gcd=exgcd(b,a%b,x,y),t=x;
     x=y,y=t-a/b*y;
     return gcd;
 }
 void ss()
 {
     const int N=;
     for(int i=;i<=N;i++)isprime[i]=;
     for(int i=;i<=N;i++)
     {
         if(isprime[i])prime[++num]=i;
         for(int j=;j<=num&&i*prime[j]<=N;j++)
         {
             isprime[i*prime[j]]=;
             if(i%prime[j]==)break;
         }
     }
 }
 void add(int x,int nu)
 {
     for(int i=;prime[i]*prime[i]<=x;i++)
     while(x%prime[i]==){cnt[prime[i]]+=nu;x/=prime[i];}
     cnt[x]+=nu;
 }
 bool pdprime(int x)
 {
     for(int i=;i*i<=x;i++)
         if(x%i==)return ;
     return ;
 }
 int fj[],cntt;
 LL poww(LL a,int b,int mod)
 {
     LL ans=;
     while(b)
     {
         if(b&)ans=ans*a%mod;
         a=a*a%mod;
         b=b>>;
     }
     return ans;
 }
 LL C(int n,int m,int mod)
 {
     if(m>n)return ;
     return jc[n]*poww(jc[m],mod-,mod)%mod*poww(jc[n-m],mod-,mod)%mod;
 }
 LL Lucas(int n,int m,int mod)
 {
     if(m>n)return ;
     if(!m)return ;
     return Lucas(n/mod,m/mod,mod)*C(n%mod,m%mod,mod)%mod;
 }
 int CRT(int W[],int B[],int k)
 {
     int x,y,a=,m,n=;
     for(int i=;i<=k;i++)
         n*=W[i];
     for(int i=;i<=k;i++)
     {
         m=n/W[i];
         exgcd(W[i],m,x,y);
         a=(a+y*m*B[i])%n;
     }
     return a>?a:(a+n);
 }
 int w[],b[];
 signed main()
 {
 //    freopen("out.doc","w",stdout);

     ss();
     cin>>T>>mod>>n>>m;int tem=mod;
     for(int i=;prime[i]*prime[i]<=tem;i++)
     while(tem%prime[i]==)
     {
         fj[++cntt]=prime[i];
         tem/=prime[i];
     }
     if(tem>)fj[++cntt]=tem;
 //    for(int i=1;i<=cntt;i++)cout<<fj[i]<<" ";puts("");
     if(T<=)
     {
         f[][][]=;
         for(int i=;i<=T;i++)
             for(int j=;j<=n+;j++)
                 for(int k=;k<=m+;k++)
                     f[i][j][k]=((f[i-][j][k-]+f[i-][j][k+])%mod+(f[i-][j-][k]+f[i-][j+][k])%mod)%mod;
         cout<<f[T][n+][m+]%mod<<endl;
         return ;
     }
     if(pdprime(mod))
     {
         if(n<)n=-n;
         if(m<)m=-m;
         jc[]=;for(int i=;i<=;i++)jc[i]=jc[i-]*i%mod;
         LL ans=;
         for(int s=m;s<=T;s++)
         {
             x=s-m;
             if((T-s-x+n)%!=)continue;
             y=(T-s-x+n)/;
             z=y-n;
             if(y<||z<)break;
             ans=(ans+jc[T]*poww(jc[s],mod-,mod)%mod*poww(jc[x],mod-,mod)%mod*poww(jc[z],mod-,mod)%mod*poww(jc[y],mod-,mod)%mod)%mod;
         }
         cout<<ans<<endl;
     }
     else
     {
         if(n<)n=-n;
         if(m<)m=-m;
         jc[]=;for(int i=;i<=;i++)jc[i]=jc[i-]*i%mod;
         LL ans=;
         for(int s=m;s<=T;s++)
         {
             x=s-m;
             if((T-s-x+n)%!=)continue;
             y=(T-s-x+n)/;
             z=y-n;
             if(y<||z<)break;
             LL tem=;
             ma(w),ma(b);
             for(int i=;i<=cntt;i++)
                 w[i]=fj[i],b[i]=Lucas(T,s,fj[i]);
             tem=tem*CRT(w,b,cntt)%mod;
             ma(w),ma(b);
             for(int i=;i<=cntt;i++)
                 w[i]=fj[i],b[i]=Lucas(T-s,x,fj[i]);
             tem=tem*CRT(w,b,cntt)%mod;
             ma(w),ma(b);
             for(int i=;i<=cntt;i++)
                 w[i]=fj[i],b[i]=Lucas(T-s-x,z,fj[i]);
             tem=tem*CRT(w,b,cntt)%mod;
             ans=(ans+tem)%mod;
         }
         cout<<ans%mod<<endl;
     }
 }