我们枚举第一个经过的坏点,然后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. Tensorflow_入门学习_1

    1.0 TensorFlow graphs Tensorflow是基于graph based computation: 如: a=(b+c)∗(c+2) 可分解为 d=b+c e=c+2 a=d∗e ...

  2. shell脚本自动部署及监控

    一.shell脚本部署nginx反向代理和三个web服务 1 对反向代理服务器进行配置 #!/bin/bash #修改用户交互页面 用户输入参数执行相应的参数 #安装epel扩展包和nginx fun ...

  3. 【转】Python 访问 HDFS

    1.前言 hdfs , Hadoop Distributed File System.Hadoop的分布式文件系统,安全行和扩展性没得说. 访问HDFS的方式有以下几种: 命令行方式:FS Shell ...

  4. java mongodb 增删改查 工具类

    package com.jttx.demo;   import com.mongodb.*; import com.mongodb.util.JSON;   import java.net.Unkno ...

  5. CF-1110 (2019/02/08)

    CF-1110 A. Parity 快速幂的思想,考虑最后一位即可 #include <bits/stdc++.h> using namespace std; typedef long l ...

  6. 【树状数组 离散化】bzoj1573: [Usaco2009 Open]牛绣花cowemb

    解方程题! Description Bessie学会了刺绣这种精细的工作.牛们在一片半径为d(1 <= d <= 50000)的圆形布上绣花. 它们一共绣了N (2 <= N < ...

  7. [JOYOI] 1035 棋盘覆盖

    题目限制 时间限制 内存限制 评测方式 题目来源 1000ms 131072KiB 标准比较器 Local 题目描述 给出一张nn(n<=100)的国际象棋棋盘,其中被删除了一些点,问可以使用多 ...

  8. 【DB_MySQL】 limit——取查询结果的子集

    语法:select * from student limit beginIndex,length; 这里结果集的下标同数组一样从0开始,beginIndex表示起始位置,length表示从beginI ...

  9. win32应用程序和win32控制台应用程序的区别

          win32应用程序是有窗体的(当然也可以没有),有Windows消息循环机制的.而win32控制台应用程序只是在控制台下运行的程序,类似以前dos的程序. 后续答案: Win32 Appl ...

  10. laravel服务容器(IOC控制反转,DI依赖注入),服务提供者,门脸模式

    laravel的核心思想: 服务容器: 容器:就是装东西的,laravel就是一个个的对象 放入:叫绑定 拿出:解析 使用容器的目的:这里面讲到的是IOC控制反转,主要是靠第三方来处理具体依赖关系的解 ...