bzoj 3782: 上学路线 dp+中国剩余定理+lucas定理

题意

小C所在的城市的道路构成了一个方形网格，它的西南角为(0,0)，东北角为(N,M)。小C家住在西南角，学校在东北角。现在有T个路口进行施工，小C不能通过这些路口。小C喜欢走最短的路径到达目的地，因此他每天上学时都只会向东或北行走；而小C又喜欢走不同的路径，因此他问你按照他走最短路径的规则，他可以选择的不同的上学路线有多少条。由于答案可能很大，所以小C只需要让你求出路径数mod P的值。
1<=N,M<=10^10
0<=T<=200
p=1000003或p=1019663265

分析

可以把所有障碍加进来排序后dp，设f[i]表示不经过任何障碍到达点i的方案，然后用全部方案减去不合法的方案即可。不合法的方案可以枚举第一次经过哪个障碍。
组合数可以用lucas定理来算。
第二个模数不是质数，可以拆开分别算后用CRT合并。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;

typedef long long LL;

const int N=205;

int n,f[N],ny[1000005],jc[1000005];
struct data{LL x,y;}a[N];

bool cmp(data a,data b)
{
    return a.x<b.x||a.x==b.x&&a.y<b.y;
}

void updata(int &x,int y,int mo)
{
    x+=y;x-=x>=mo?mo:0;
}

int get_c(LL n,LL m,int mo)
{
    if (n<m) return 0;
    if (!n) return 1;
    if (n<mo&&m<mo) return (LL)jc[n]*ny[m]%mo*ny[n-m]%mo;
    return (LL)get_c(n/mo,m/mo,mo)*get_c(n%mo,m%mo,mo)%mo;
}

int solve(int mo)
{
    jc[0]=jc[1]=ny[0]=ny[1]=1;
    for (int i=2;i<mo;i++) jc[i]=(LL)jc[i-1]*i%mo,ny[i]=(LL)(mo-mo/i)*ny[mo%i]%mo;
    for (int i=2;i<mo;i++) ny[i]=(LL)ny[i-1]*ny[i]%mo;
    f[1]=1;
    for (int i=2;i<=n;i++)
    {
        f[i]=get_c(a[i].x+a[i].y,a[i].y,mo);
        for (int j=2;j<i;j++)
            if (a[j].x<=a[i].x&&a[j].y<=a[i].y) updata(f[i],mo-(LL)f[j]*get_c(a[i].x+a[i].y-a[j].x-a[j].y,a[i].x-a[j].x,mo)%mo,mo);
    }
    return f[n];
}

void exgcd(LL &x,LL &y,LL a,LL b)
{
    if (!b) {x=1;y=0;return;}
    exgcd(x,y,b,a%b);
    LL tmp=x;x=y;y=tmp-(a/b)*y;
}

int merge(LL x1,LL mo1,LL x2,LL mo2)
{
    LL x,y,mo=mo1*mo2;
    exgcd(x,y,mo1,mo2);
    return ((x*(x2-x1)%mo*mo1+x1)%mo+mo)%mo;
}

int main()
{
    LL x,y,w;
    scanf("%lld%lld%d%lld",&x,&y,&n,&w);
    for (int i=1;i<=n;i++) scanf("%lld%lld",&a[i].x,&a[i].y);
    a[++n]=(data){x,y};a[++n]=(data){0,0};
    sort(a+1,a+n+1,cmp);
    if (w==1000003) printf("%d\n",solve(w));
    else
    {
        int ans=solve(3);
        ans=merge(ans,3,solve(5),5);
        ans=merge(ans,15,solve(6793),6793);
        ans=merge(ans,15*6793,solve(10007),10007);
        printf("%d",ans);
    }
    return 0;
}