JZOJ 5612 T3

Description

有个人要从($0$,$0$)走到($n$,$m$)，假设他现在在($x$,$y$)，那么他可以走到($x$+$1$,$y$+$1$),($x$+$1$,$y$-$1$)或($x$+$2$,$y$)，求方案数模100003。

Data Constraint

$n,m$<=$10^{18}$

Solution

首先变换坐标系，逆时针旋转45度，则原题目变成，每次可以从($x$,$y$)走到($x$,$y+1$)、($x+1$,$y$)或($x+1$,$y+1$)，求从($0$,$0$)走到(${n-m}\over 2$,${n+m}\over 2$)的方案数。

枚举先向右走然后再向上走的路径总共有多少个，定义只有这样的两段路径才能被($x$,$y$)走到($x+1$,$y+1$)的一步替代，每一个这样的路径都能被替代或者不被替代两种选择，所以有

$$Ans=\sum_{i=1}^{{n-m}\over 2}C_{{n+m}\over2}^iC_{{n-m}\over 2}^i2^i$$
由于模数比较小，直接用 $Lucas$定理拆位统计答案即可。

Code

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>

#define fo(i,j,l) for(int i=j;i<=l;++i)
#define fd(i,j,l) for(int i=j;i>=l;--i)

using namespace std;
typedef long long ll;
const ll N=1e6,mo=1e5+3,P=mo+10;

ll jc[N],ny[N],n,m;
ll f[5][P][3],m2[N];
ll fn[5],fm[5];

inline ll ksm(ll o,ll t)
{
    ll y=1;
    for(;t;t>>=1,o=o*o%mo)
    if(t&1)y=y*o%mo;
    return y;
}

ll C(ll a,ll b)
{return a<b?0:(a<mo?jc[a]*ny[b]%mo*ny[a-b]%mo:C(a/mo,b/mo)*C(a%mo,b%mo)%mo);}

int main()
{
    jc[0]=ny[0]=1;
    m2[0]=1;
    fo(i,1,mo)m2[i]=m2[i-1]*2%mo;
    fo(i,1,mo-1)jc[i]=jc[i-1]*i%mo;
    ny[mo-1]=ksm(jc[mo-1],mo-2);
    fd(i,mo-2,1)ny[i]=ny[i+1]*(i+1)%mo;
    cin>>n>>m; ll op=n;
    n=(n+m)>>1; m=(op-m)>>1;
    for(;n;n=n/mo)fn[++fn[0]]=n%mo;
    for(;m;m=m/mo)fm[++fm[0]]=m%mo;
    ll ans=1;
    fo(l,1,fn[0]){
        ll lj=0;
        fo(i,0,min(fm[l],fn[l]))lj=(lj+C(fm[l],i)*C(fn[l],i)%mo*m2[i])%mo;
        ans=ans*lj%mo;
    }
    cout<<ans;
}