题意
小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;
}