【BZOJ3782】上学路线 组合数+容斥+CRT
【BZOJ3782】上学路线
Description
Input
Output
Sample Input
3 0
1 1
2 2
Sample Output
HINT
题解:从(0,0)走到(n,m)的总方案数=C(n+m,n)。
依旧考虑容斥,先将点排序,用f[i]表示从(0,0)走到(x[i],y[i]),途中不经过其它障碍的方案数,那么如果j在i的左下方,则f[i]-=f[j]*(从j走到i的方案数)。
然而1019663265不是质数?分解质因数的1019663265=3*5*6793*10007,分别求解,再用中国剩余定理合并即可。
EXGCD还能写错~
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
ll mod,ans;
ll n,m,P;
int T;
ll f[210],jc[1000010],jcc[1000010],ine[1000010];
struct node
{
ll x,y;
}p[210];
bool cmp(node a,node b)
{
return (a.x==b.x)?(a.y<b.y):(a.x<b.x);
}
ll C(ll a,ll b)
{
if(a<b) return 0;
if(!b) return 1;
if(a<mod&&b<mod) return jc[a]*jcc[a-b]%mod*jcc[b]%mod;
return C(a%mod,b%mod)*C(a/mod,b/mod)%mod;
}
ll calc()
{
memset(f,0,sizeof(f));
int i,j;
jc[0]=jcc[0]=1;
ine[1]=1;
for(i=2;i<mod;i++) ine[i]=(mod-(mod/i)*ine[mod%i]%mod)%mod;
for(i=1;i<mod;i++) jc[i]=jc[i-1]*i%mod,jcc[i]=jcc[i-1]*ine[i]%mod;
for(i=1;i<=T;i++)
{
f[i]=C(p[i].x+p[i].y,p[i].x);
for(j=1;j<i;j++)
if(p[i].x>=p[j].x&&p[i].y>=p[j].y) f[i]=(f[i]-f[j]*C(p[i].x-p[j].x+p[i].y-p[j].y,p[i].x-p[j].x)%mod+mod)%mod;
}
return f[T];
}
void exgcd(ll a,ll b,ll &x,ll &y)
{
if(!b)
{
x=1,y=0;
return ;
}
exgcd(b,a%b,x,y);
ll tmp=x;
x=y,y=tmp-a/b*x;
}
ll work(ll a,ll c,ll b)
{
ll x,y;
exgcd(a,b,x,y),x=(x+b)%b;
return x*a%P*c%P;
}
int main()
{
scanf("%lld%lld%d%lld",&n,&m,&T,&P);
int i;
for(i=1;i<=T;i++) scanf("%lld%lld",&p[i].x,&p[i].y);
p[++T].x=n,p[T].y=m;
sort(p+1,p+T+1,cmp);
if(P==1000003)
{
mod=P,printf("%lld",calc());
return 0;
}
ll a1,a2,a3,a4;
mod=3,a1=calc();
mod=5,a2=calc();
mod=6793,a3=calc();
mod=10007,a4=calc();
ans=(ans+work(P/3,a1,3))%P;
ans=(ans+work(P/5,a2,5))%P;
ans=(ans+work(P/6793,a3,6793))%P;
ans=(ans+work(P/10007,a4,10007))%P;
printf("%lld",(ans+P)%P);
return 0;
}//3 4 3 1000003 3 0 1 1 2 2
【BZOJ3782】上学路线 组合数+容斥+CRT的更多相关文章
- BZOJ3782 上学路线 【dp + Lucas + CRT】
题目链接 BZOJ3782 题解 我们把终点也加入障碍点中,将点排序,令\(f[i]\)表示从\((0,0)\)出发,不经过其它障碍,直接到达\((x_i,y_i)\)的方案数 首先我们有个大致的方案 ...
- bzoj3782上学路线(Lucas+CRT+容斥DP+组合计数)
传送门:https://www.lydsy.com/JudgeOnline/problem.php?id=3782 有部分分的传送门:https://www.luogu.org/problemnew/ ...
- bzoj3782上学路线
题意:从n*m网格图的左下角走到右上角(n,m<=10^10),有t个坐标不能经过(t<=200),只能向上向右走,问有多少种不同的走法,对p取模, p只有两种取值,1000003(质数) ...
- Codeforces 100548F - Color (组合数+容斥)
题目链接:http://codeforces.com/gym/100548/attachments 有n个物品 m种颜色,要求你只用k种颜色,且相邻物品的颜色不能相同,问你有多少种方案. 从m种颜色选 ...
- BZOJ5306 [HAOI2018]染色 【组合数 + 容斥 + NTT】
题目 为了报答小 C 的苹果, 小 G 打算送给热爱美术的小 C 一块画布, 这块画布可 以抽象为一个长度为 \(N\) 的序列, 每个位置都可以被染成 \(M\) 种颜色中的某一种. 然而小 C 只 ...
- 【BZOJ4710】[Jsoi2011]分特产 组合数+容斥
[BZOJ4710][Jsoi2011]分特产 Description JYY 带队参加了若干场ACM/ICPC 比赛,带回了许多土特产,要分给实验室的同学们. JYY 想知道,把这些特产分给N 个同 ...
- cf997C. Sky Full of Stars(组合数 容斥)
题意 题目链接 \(n \times n\)的网格,用三种颜色染色,问最后有一行/一列全都为同一种颜色的方案数 Sol Orz fjzzq 最后答案是这个 \[3^{n^2} - (3^n - 3)^ ...
- BZOJ3782 上学路线
设障碍个数为,\(obs\)则一般的容斥复杂度为\(O(2^{obs})\).但因为这个题是网格图,我们可以用DP解.设\(f[i]\)表示不经过任何障碍到达第\(i\)个障碍的方案数,转移时枚举可以 ...
- HDU - 5201 :The Monkey King (组合数 & 容斥)
As everyone known, The Monkey King is Son Goku. He and his offspring live in Mountain of Flowers and ...
随机推荐
- hdu 4353 统计点在三角形内的个数
Finding Mine Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Tota ...
- javaweb学习总结(十五)——JSP基础语法(转)
任何语言都有自己的语法,JAVA中有,JSP虽然是在JAVA上的一种应用,但是依然有其自己扩充的语法,而且在JSP中,所有的JAVA语句都可以使用. 一.JSP模版元素 JSP页面中的HTML内容称之 ...
- GridView从行寻找到该控件,以及从该控件获知该行
原文发布时间为:2008-08-01 -- 来源于本人的百度文章 [由搬家工具导入] GridView从行寻找到该控件: protected void Button3_Click(object sen ...
- WEB学习-兼容问题
css选择器 儿子选择器 (IE7开始兼容,IE6不兼容.) div>p{ color:red; } div的儿子p.和div的后代p的截然不同. 能够选择: <div> <p ...
- 关于Red5整合springMVC提示scope not found 的错误
https://www.cnblogs.com/qgc88:
- AC日记——猴子 cogs 2043
2043. 猴子 ★★ 输入文件:monkeya.in 输出文件:monkeya.out 简单对比时间限制:1 s 内存限制:256 MB [题目描述] 有n只猴子,第一只尾巴挂在树上 ...
- VUE之命令行报错:Expected indentation of 4 spaces but found 6
使用vue时候,经常被一大片警告惊呆了,这是webpack默认的语法检查插件ESLint在做警告, [ESLint是一个语法规则和代码风格的检查工具,可以用来保证写出语法正确.风格统一的代码] 但是我 ...
- Java 界面编程【03】事件监听
当你把界面都设计好了,总需要添加相应的执行动作给组件,在JAVA中有相应的时间处理机制,叫做“监听器”,给组件添加相应执行动作的过程叫做“注册”,其中的“监听器”是一个接口,里面包含了相应的执行函数, ...
- HDU 5733 tetrahedron(计算几何)
题目链接 tetrahedron 题目大意 输入一个四面体求其内心,若不存在内心则输出"O O O O" 解题思路 其实这道题思路很简单,只要类推一下三角形内心公式就可以了. 至于 ...
- 洛谷——P2658 汽车拉力比赛
P2658 汽车拉力比赛 题目描述 博艾市将要举行一场汽车拉力比赛. 赛场凹凸不平,所以被描述为M*N的网格来表示海拔高度(1≤ M,N ≤500),每个单元格的海拔范围在0到10^9之间. 其中一些 ...