我们枚举第一个经过的坏点,然后DP即可。

状态转移方程不是难点,难点在于组合数的处理。

将狼踩尽的博客中有很详细的证明过程,但是我只记住了结论

$n=a_1 * p^k+a_2*p^k-1...$

$m=b_1 * p^k+b_2*p^k-1...$

$C(_{m}^{n})=C(_{b_1}^{a_1})*...$

大概的意思就是转化成$p$进制下的每一位做组合数,那么我们就可以预处理阶乘以及它的逆元进行计算。

所以说Lucas只能跑过$10^5$当质数很大的时候就放弃。

如果不是质数,那么可以分解质因数,每一个因数做一次Lucas,然后用CRT合并。

以前留下的大坑终于补完了,EXGCD和CRT终于明白了(真是弱(。・・)ノ)

突然发现namespace写起来挺好用的,以后挂链就用它了

注意 1.n<m时候需要判定,因为%意义下会有小的情况产生。

2.随时取模

#include <map>
#include <cmath>
#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define F(i,j,k) for (int i=j;i<=k;++i)
#define D(i,j,k) for (int i=j;i>=k;--i)
#define ll long long
#define mp make_pair namespace Subtask1{
const int p=1000003;
ll fac[p],inv[p];
void Shaker()
{
int i;
for (fac[0]=1,i=1;i<p;++i)
fac[i]=fac[i-1]*i%p;
for (inv[1]=1,i=2;i<p;++i)
inv[i]=(p-p/i)*inv[p%i]%p;
for (inv[0]=1,i=1;i<p;++i)
(inv[i]*=inv[i-1])%=p;
}
ll C(ll n,ll m)
{
if (n<m) return 0;
if (n<p&&m<p) return fac[n]*inv[m]*inv[n-m]%p;
return C(n%p,m%p)*C(n/p,m/p)%p;
}
} namespace Subtask2{
const int p=1019663265;
int pri[]={3,5,6793,10007};
ll fac[4][10007],inv[4][10007],k1=339887755,k2=407865306,k3=673070820,k4=618502650;
void Shaker()
{
int i,j;
for (j=0;j<4;++j)
{
int np=pri[j];
for (fac[j][0]=1,i=1;i<np;++i)
fac[j][i]=fac[j][i-1]*i%np;
for (inv[j][1]=1,i=2;i<np;++i)
inv[j][i]=(np-np/i)*inv[j][np%i]%np;
for (inv[j][0]=1,i=1;i<np;++i)
(inv[j][i]*=inv[j][i-1])%=np;
}
}
ll C(ll n,ll m,int j)
{
int np=pri[j];
if (n<m) return 0;
if (n<np&&m<np) return fac[j][n]*inv[j][m]*inv[j][n-m]%np;
return C(n%np,m%np,j)*C(n/np,m/np,j)%np;
}
ll C(ll n,ll m)
{
ll r1=C(n,m,0),r2=C(n,m,1),r3=C(n,m,2),r4=C(n,m,3);
return (r1*k1+r2*k2+r3*k3+r4*k4)%p;
}
} struct Point{
ll x,y;
void read(){scanf("%lld%lld",&x,&y);}
}points[220]; bool cmp(Point a,Point b)
{return a.x==b.x?a.y<b.y:a.x<b.x;} ll ksm(ll a,ll b,ll p)
{
ll ret=1;
for (;b;b>>=1,(a*=a)%=p)
if (b&1) (ret*=a)%=p;
return ret;
} ll f[220],t,n,m; void solve1()
{
using namespace Subtask1;
Shaker();
F(i,1,t)
{
f[i]=C(points[i].x+points[i].y,points[i].x);
F(j,1,i-1) if (points[j].y<=points[i].y)
{
(f[i]-=f[j]*C(points[i].x-points[j].x+points[i].y-points[j].y,points[i].x-points[j].x)%p);
f[i]+=p; f[t]%=p;
}
}
} void solve2()
{
using namespace Subtask2;
Shaker();
F(i,1,t)
{
f[i]=C(points[i].x+points[i].y,points[i].x);
F(j,1,i-1) if (points[j].y<=points[i].y)
{
(f[i]-=f[j]*C(points[i].x-points[j].x+points[i].y-points[j].y,points[i].x-points[j].x)%p);
f[i]%=p;f[i]+=p; f[t]%=p;
}
}
} void Finout()
{
freopen("in.txt","r",stdin);
freopen("wa.txt","w",stdout);
} int main()
{
ll p;
scanf("%lld%lld%lld%lld",&n,&m,&t,&p);
F(i,1,t) points[i].read();++t;points[t].x=n;points[t].y=m;
sort(points+1,points+t+1,cmp);
if (p==1000003) solve1(); else solve2();
cout<<f[t]<<endl;
}

  

BZOJ 3782 上学路线 ——动态规划 Lucas定理 中国剩余定理的更多相关文章

  1. BZOJ 3782: 上学路线 [Lucas定理 DP]

    3782: 上学路线 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 192  Solved: 75[Submit][Status][Discuss] ...

  2. 【bzoj3782】上学路线 dp+容斥原理+Lucas定理+中国剩余定理

    题目描述 小C所在的城市的道路构成了一个方形网格,它的西南角为(0,0),东北角为(N,M).小C家住在西南角,学校在东北角.现在有T个路口进行施工,小C不能通过这些路口.小C喜欢走最短的路径到达目的 ...

  3. bzoj 3782 上学路线 卢卡斯定理 容斥 中国剩余定理 dp

    LINK:上学路线 从(0,0)走到(n,m)每次只能向上或者向右走 有K个点不能走求方案数,对P取模. \(1\leq N,M\leq 10^10 0\leq T\leq 200\) p=10000 ...

  4. 卢卡斯定理&&中国剩余定理

    卢卡斯定理(模数较小,且是质数) 式子C(m,n)=C(m/p,n/p)*C(m%p,n%p)%p 至于证明(我也不会QAQ,只要记住公式也该就好了). 同时卢卡斯定理一般用于组合数取模上 1.首先当 ...

  5. BZOJ 3782 上学路线

    首先这个题需要dp.dp[i]=C(x[i]+y[i],x[i])-Σdp[j]*C(x[i]-x[j]+y[i]-y[j],x[i]-x[j])(x[i]>=x[j],y[i]>=y[j ...

  6. 【bzoj1951】[Sdoi2010]古代猪文 费马小定理+Lucas定理+中国剩余定理

    题目描述 求  $g^{\sum\limits_{k|n}C_{n}^{\frac nk}}\mod 999911659$ 输入 有且仅有一行:两个数N.G,用一个空格分开. 输出 有且仅有一行:一个 ...

  7. Ceizenpok’s formula Gym - 100633J 扩展Lucas定理 + 中国剩余定理

    http://codeforces.com/gym/100633/problem/J 其实这个解法不难学的,不需要太多的数学.但是证明的话,我可能给不了严格的证明.可以看看这篇文章 http://ww ...

  8. hdu 5446 Unknown Treasure Lucas定理+中国剩余定理

    Unknown Treasure Time Limit: 1500/1000 MS (Java/Others)    Memory Limit: 131072/131072 K (Java/Other ...

  9. 【题解】P2480 [SDOI2010]古代猪文 - 卢卡斯定理 - 中国剩余定理

    P2480 [SDOI2010]古代猪文 声明:本博客所有题解都参照了网络资料或其他博客,仅为博主想加深理解而写,如有疑问欢迎与博主讨论✧。٩(ˊᗜˋ)و✧*。 题目描述 猪王国的文明源远流长,博大精 ...

随机推荐

  1. 【排序】插入排序:最稳定:时间复杂度O(n^2)

    想象着自己在玩扑克的时候抓牌,每抓到一张牌,按照从小到大的顺序排序. 如果第二张的点数小于第一张,就交换这两张牌,默认每次抓牌之前,前面的已经排好序了. 再来一张牌,与第二张比较,如果小于第二张,交换 ...

  2. C#调用CMD程序

    最近写了两个小程序都要调用Windows自带的命令行程序,一个是调用Openfiles.exe查询正在编辑的共享文档,一个是调用DiskPart.exe查询硬盘状态.两种命令行程序调用有点不同,记录一 ...

  3. mongodb主从配置信息查看与确认

    在local库中不仅有主从日志 oplog集合,还有一个集合用于记录主从配置信息 system.replset: > use local > show collections > d ...

  4. SAP CRM中间件下载equipment时遇到的一个错误

    在CRM开发系统上进行equipment下载,发现不工作.调试发现错误信息在下图定96行的WHEN default分支抛出的: MESSAGE ID 'AZ' ... 通过阅读源代码发现,ERP端支持 ...

  5. Xcode5 如何添加一个Github/Repository 并且Checkout

    1. 添加一个Account  也就是添加一个 Repository. In Xcode, choose Xcode > Preferences, and click Accounts. Pre ...

  6. 卷积网络中的通道(Channel)和特征图

    转载自:https://www.jianshu.com/p/bf8749e15566 今天介绍卷积网络中一个很重要的概念,通道(Channel),也有叫特征图(feature map)的. 首先,之前 ...

  7. bootstrap 两端对齐的导航

    您可以在屏幕宽度大于768px时,通过在分别使用.nav .nav-tabs或.nav .nav-pills的同时使用class.nav-justified,让标签式或胶囊式导航菜单与父元素等宽,在更 ...

  8. dedeCMS数据库字段详细介绍

    dede_addonarticle 附加文章表 aid int(11) 文章编号 typeid int(11) 分类栏目编号 body mediumtext 文章内容 dede_addonflash ...

  9. shell脚本,检查给出的字符串是否为回文

    [root@localhost wyb]# .sh #!/bin/bash #检查给出的字符串是否为回文 read -p "Please input a String:" numb ...

  10. non-JRMP server at remote endpoint

    #在相应的domain的domain.xml文件添加下面红色设置,并重启domain <admin-service system-jmx-connector-name="system& ...