题目大意

从 \((0,0)\) 开始,每次只可走上下左右一个单位长度,可走重复路,求第 \(T\) 步正好走到 \((n,m)\) 的方案数。

答案要求对 \(MOD\) 取模,\(MOD\) 保证是几个不同质数的乘积。

数据范围

\(1\leq T\leq 100000;-T\leq n,m\leq T;1\leq MOD\leq 10^9+7\)

样例输入

4 10

2 2

样例输出

6

思路

枚举向左走了 \(l\) 步,则向右走了 \(r\) 步,向上走了 \(u\) 步,向下走了 \(d\) 步,依据高考数学中的平均分组问题,答案就是:

\[\frac{A_T^T}{A_l^l\ A_r^r\ A_u^u\ A_d^d}=C_T^l\ C_{T-l}^r\ C_{T-l-r}^u
\]

但是这个式子比较繁琐不好计算,但是聚聚 \(\texttt{skyh}\) 给出了一个简单的式子:

\[C_T^{\frac{T-n-m}{2}}\ C_T^{\frac{T-\vert n-m\vert}{2}}
\]

问了数奥的同学该式子的正确性,请大家学习

显然使用 \(Lucas\) 定理求解。

噫,好,我会了!

然后一顿怒操作发现样例输出 \(0\)。这是因为在模 \(10\) 的意义下样例是没有逆元的。所以我们需要将模数分解质因子,在每一个质因子下的模意义下求出答案,最后用 \(CRT\) 合并答案即可。

代码

数论全家桶

#include <bits/stdc++.h>
#define int long long
using namespace std;
const int maxn=1e5+10;
int T,Mod,N,M;
int res[maxn],fac[maxn],inv[maxn]; inline int read(){
int x=0;bool fopt=1;char ch=getchar();
for(;!isdigit(ch);ch=getchar())if(ch=='-')fopt=0;
for(;isdigit(ch);ch=getchar())x=(x<<3)+(x<<1)+ch-48;
return fopt?x:-x;
} void exgcd(int a,int b,int &x,int &y){
if(!b)return x=1,y=0,void();
exgcd(b,a%b,x,y);
int z=x;x=y;y=z-(a/b)*y;
} inline int qpow(int x,int b,int p){
int ans=1,base=x;
while(b){
if(b&1)ans=ans*base%p;
base=base*base%p;
b>>=1;
}
return ans;
} bool Miller_Rabin(int n){
if(n==2)
return true;
for(int i=1;i<=30;i++){
int a=rand()%(n-2)+2;
if(qpow(a,n,n)!=a)
return false;
}
return true;
} int v[maxn];
inline void Get(int x){
int u=sqrt(x);
for(int i=2;i<=u;i++){
if(x%i!=0)continue;
if(Miller_Rabin(i)){
v[++v[0]]=i;
x/=i;
}
}
if(x>1)v[++v[0]]=x;
} inline int C(int n,int m,int p){
if(n<m)return 0;
if(!m||n==m)return 1;
return fac[n]*inv[m]%p*inv[n-m]%p;
} int Lucas(int n,int m,int p){
if(n<m)return 0;
if(!m)return 1;
return C(n%p,m%p,p)*Lucas(n/p,m/p,p)%p;
} inline int CRT(int n){
int M=1,ans=0;
for(int i=1;i<=n;i++)
M*=v[i];
for(int i=1;i<=n;i++){
int m=M/v[i];
int x=0,y=0;
exgcd(m,v[i],x,y);
ans+=res[i]*m*(x<0?x+v[i]:x);
}
return ans%M;
} signed main(){
srand(time(0));
T=read();Mod=read();N=read();M=read();
Get(Mod);
for(int i=1;i<=v[0];i++){
int p=v[i];
fac[0]=fac[1]=inv[0]=inv[1]=1;
for(int j=2;j<=T;j++)
inv[j]=(p-p/j)*inv[p%j]%p;
for(int j=2;j<=T;j++){
inv[j]=inv[j-1]*inv[j]%p;
fac[j]=fac[j-1]*j%p;
}
res[i]=Lucas(T,(T-N-M)/2,p)*Lucas(T,(T-abs(N-M))/2,p)%p;
}
printf("%lld\n",CRT(v[0]));
return 0;
}

【组合计数】visit的更多相关文章

  1. bzoj 2281 [Sdoi2011]黑白棋(博弈+组合计数)

    黑白棋(game) [问题描述] 小A和小B又想到了一个新的游戏. 这个游戏是在一个1*n的棋盘上进行的,棋盘上有k个棋子,一半是黑色,一半是白色. 最左边是白色棋子,最右边是黑色棋子,相邻的棋子颜色 ...

  2. BZOJ 4555: [Tjoi2016&Heoi2016]求和 [分治FFT 组合计数 | 多项式求逆]

    4555: [Tjoi2016&Heoi2016]求和 题意:求\[ \sum_{i=0}^n \sum_{j=0}^i S(i,j)\cdot 2^j\cdot j! \\ S是第二类斯特林 ...

  3. BZOJ 4555: [Tjoi2016&Heoi2016]求和 [FFT 组合计数 容斥原理]

    4555: [Tjoi2016&Heoi2016]求和 题意:求\[ \sum_{i=0}^n \sum_{j=0}^i S(i,j)\cdot 2^j\cdot j! \\ S是第二类斯特林 ...

  4. 【BZOJ5491】[HNOI2019]多边形(模拟,组合计数)

    [HNOI2019]多边形(模拟,组合计数) 题面 洛谷 题解 突然特别想骂人,本来我考场现切了的,结果WA了几个点,刚刚拿代码一看有个地方忘记取模了. 首先发现终止态一定是所有点都向\(n\)连边( ...

  5. [总结]数论和组合计数类数学相关(定理&证明&板子)

    0 写在前面 0.0 前言 由于我太菜了,导致一些东西一学就忘,特开此文来记录下最让我头痛的数学相关问题. 一些引用的文字都注释了原文链接,若侵犯了您的权益,敬请告知:若文章中出现错误,也烦请告知. ...

  6. 【BZOJ5323】[JXOI2018]游戏(组合计数,线性筛)

    [BZOJ5323][JXOI2018]游戏(组合计数,线性筛) 题面 BZOJ 洛谷 题解 显然要考虑的位置只有那些在\([l,r]\)中不存在任意一个约数的数. 假设这样的数有\(x\)个,那么剩 ...

  7. 【BZOJ5305】[HAOI2018]苹果树(组合计数)

    [BZOJ5305][HAOI2018]苹果树(组合计数) 题面 BZOJ 洛谷 题解 考虑对于每条边计算贡献.每条边的贡献是\(size*(n-size)\). 对于某个点\(u\),如果它有一棵大 ...

  8. 【BZOJ3142】[HNOI2013]数列(组合计数)

    [BZOJ3142][HNOI2013]数列(组合计数) 题面 BZOJ 洛谷 题解 唯一考虑的就是把一段值给分配给\(k-1\)天,假设这\(k-1\)天分配好了,第\(i\)天是\(a_i\),假 ...

  9. 【BZOJ4005】[JLOI2015] 骗我呢(容斥,组合计数)

    [BZOJ4005][JLOI2015] 骗我呢(容斥,组合计数) 题面 BZOJ 洛谷 题解 lalaxu #include<iostream> using namespace std; ...

  10. 【BZOJ4830】[HNOI2017]抛硬币(组合计数,拓展卢卡斯定理)

    [BZOJ4830][HNOI2017]抛硬币(组合计数,拓展卢卡斯定理) 题面 BZOJ 洛谷 题解 暴力是啥? 枚举\(A\)的次数和\(B\)的次数,然后直接组合数算就好了:\(\display ...

随机推荐

  1. [程序员代码面试指南]递归和动态规划-换钱的最少货币数(DP,完全背包)

    题目描述 给定arr,arr中所有的值都为正数且不重复.每个值代表一种面值的货币,每种面值的货币可以使用任意张,再给定一个整数aim,求组成aim的最少货币数. 解题思路 dp[i][j]表示只用第0 ...

  2. 跨平台框架与React Native基础

    跨平台框架 什么是跨平台框架? 这里的多个平台一般是指 iOS 和 Android . 为什么需要跨平台框架? 目前,移动开发技术主要分为原生开发和跨平台开发两种.其中,原生应用是指在某个特定的移动平 ...

  3. Vue iview可收缩多级菜单的实现

    递归组件实战 views/layout.vue <template> <div class="layout-wrapper"> <Layout cla ...

  4. bash运行脚本的几种方式

    转载自https://www.jianshu.com/p/ba6efda13e23 转载地址:http://www.jquerycn.cn/a_8354 bash shell 脚本执行的方法有多种,本 ...

  5. hystrix源码之AbstractCommand

    AbstractCommand HystrixCommand和HystrixObservableCommand的父类.每个command对应一个HystrixCommandKey.HystrixThr ...

  6. RXJAVA之变换操作

    RXJAVA提供了以下变换操作,对Observable的消息进行变换操作: 1.window 定期将来自Observable的数据分拆成一些Observable窗口,然后发射这些窗口,而不是每次发射一 ...

  7. FreeSWITCH 处理Refer盲转时,UUI传递不对(没有将SIP 消息头Refer-To中的User-to-User传递给B-Leg)

    运行环境:     CentOS 7.6     FreeSWICH 1.6.18   一.问题场景:     FreeSWITCH收到REFER命令后,重新发起的INVITE消息中的 "U ...

  8. 结构体排序中sort的自定义函数cmp()

    水题链接 #include<iostream> #include<cstdio> #include<algorithm> using namespace std; ...

  9. Spring学习(一)初识Spring

    什么是Spring 定义:Spring 是一个轻量级的 DI / IoC 和 AOP 容器的开源框架,目的为了简化java开发. DI:注入 IOC:控制反转 AOP:面向切面编程 原理:利用了jav ...

  10. [HarekazeCTF2019]Avatar Uploader 1 &&

    [HarekazeCTF2019]Avatar Uploader 1 这是一个文件上传的题目,但是这导体是通过满足条件来获取flag的. 他有两个函数,一个是getimagesize,还有一个是FIL ...