卡特兰数(高精度乘法压位)
889. 满足条件的01序列
- 题目
- 提交记录
- 讨论
- 题解
- 视频讲解
给定 n
个 0 和 n 个 1,它们将按照某种顺序排成长度为 2n 的序列,求它们能排列成的所有序列中,能够满足任意前缀序列中 0 的个数都不少于 1
的个数的序列有多少个。
输出的答案对 109+7
取模。
输入格式
共一行,包含整数 n
。
输出格式
共一行,包含一个整数,表示答案。
数据范围
1≤n≤105
输入样例:
3
输出样例:
5
* 假定在一个二维矩阵的地图里,向右走为0,向上走为1,则要求能排列成的所有序列中,
* 能够满足任意前缀序列中 0 的个数都不少于 1 的个数的序列,所有点都满足为或者在
* (0,0) 到 (n,n) 的直线上或者在下方。只要有点越过直线 y=x+1 ,再到达点(n,n),
* 则这种排列就是不符合要求的。所以不符合要求的走法一定满足从(0,0) 到 (n-1,n+1) 的
* 路线关于 y=x+1 对称 (因为(0,0) 到 (n,n+1) 的路线一定是不合法的,并且(n,n)与
* (n-1,n+1) 关于 y=x+1 对称。)。那么我们只要 C(2*n,n) - C(2*n,n-1) 就是答案。
*
* C(2*n,n) - C(2*n,n-1) =C(2*n,n)/(n+1);()
*
* 卡特兰数:
* C(2*n,n) - C(2*n,n-1) =(2*n)! / (n! * n! * (n+1))
/*** 假定在一个二维矩阵的地图里,向右走为0,向上走为1,则要求能排列成的所有序列中,* 能够满足任意前缀序列中 0 的个数都不少于 1 的个数的序列,所有点都满足为或者在* (0,0) 到 (n,n) 的直线上或者在下方。只要有点越过直线 y=x+1 ,再到达点(n,n),* 则这种排列就是不符合要求的。所以不符合要求的走法一定满足从(0,0) 到 (n-1,n+1) 的* 路线关于 y=x+1 对称 (因为(0,0) 到 (n,n+1) 的路线一定是不合法的,并且(n,n)与* (n-1,n+1) 关于 y=x+1 对称。)。那么我们只要 C(2*n,n) - C(2*n,n-1) 就是答案。* * C(2*n,n) - C(2*n,n-1) =C(2*n,n)/(n+1);()* * 卡特兰数:* C(2*n,n) - C(2*n,n-1) =(2*n)! / (n! * n! * (n+1))
*/#include <iostream>using namespace std;typedef long long LL;const int mod=1e9+7,maxn=1e5+10;
int fac[2*maxn],infac[maxn];int qmi(int a,int b,int p)
{int res=1;while(b){if(b&1)res=(LL)res*a%mod;a=(LL)a*a%mod;b>>=1;}return res;
}void init(int n)
{fac[0] = infac[0] = 1; //初始化for(int i=1;i<=2*n;++i)fac[i]=(LL)fac[i-1]*i%mod;for(int i=1;i<=n;++i)infac[i]=(LL)infac[i-1]*qmi(i,mod-2,mod)%mod; //必须要用qmi求逆元
}int main()
{int n;cin >> n;init(n);cout << (LL)fac[2*n]*infac[n]%mod*infac[n]%mod *qmi(n+1,mod-2,mod)%mod << endl;return 0;
}
130. 火车进出栈问题
- 题目
- 提交记录
- 讨论
- 题解
- 视频讲解
一列火车 n
节车厢,依次编号为 1,2,3,…,n
。
每节车厢有两种运动方式,进栈与出栈,问 n
节车厢出栈的可能排列方式有多少种。
输入格式
输入一个整数 n
,代表火车的车厢数。
输出格式
输出一个整数 s
表示 n
节车厢出栈的可能排列方式数量。
数据范围
1≤n≤60000
输入样例:
3
输出样例:
5
* 卡特兰数:
* C(2*n,n) - C(2*n,n-1) =(2*n)! / (n! * n! * (n+1));
*
* 这个题与满足条件的01序列解法一样,只不过这个题是输出全部方案数,不取模,
* 输出高精度数据;由于数据较大,所以我们用质因数分解的方法来求解。
* 最后由于高精度数据相乘是主要的算法时间消耗,因此我们优化高精度数据相乘,
* 用LL来存储结果,每次相乘我们以 M=1e9 作为进制数,(这样两个M相乘也不会爆LL);
* 这样能减少许多乘法步骤;这种方法称为高精度压位。
/*** 假定在一个二维矩阵的地图里,向右走为0,向上走为1,则要求能排列成的所有序列中,* 能够满足任意前缀序列中 0 的个数都不少于 1 的个数的序列,所有点都满足为或者在* (0,0) 到 (n,n) 的直线上或者在下方。只要有点越过直线 y=x+1 ,再到达点(n,n),* 则这种排列就是不符合要求的。所以不符合要求的走法一定满足从(0,0) 到 (n-1,n+1) 的* 路线关于 y=x+1 对称 (因为(0,0) 到 (n,n+1) 的路线一定是不合法的,并且(n,n)与* (n-1,n+1) 关于 y=x+1 对称。)。那么我们只要 C(2*n,n) - C(2*n,n-1) 就是答案。* * C(2*n,n) - C(2*n,n-1) =C(2*n,n)/(n+1);()* * 卡特兰数:* C(2*n,n) - C(2*n,n-1) =(2*n)! / (n! * n! * (n+1));* * 这个题与满足条件的01序列解法一样,只不过这个题是输出全部方案数,不取模,* 输出高精度数据;由于数据较大,所以我们用质因数分解的方法来求解。* 最后由于高精度数据相乘是主要的算法时间消耗,因此我们优化高精度数据相乘,* 用LL来存储结果,每次相乘我们以 M=1e9 作为进制数,(这样两个M相乘也不会爆LL);* 这样能减少许多乘法步骤;这种方法称为高精度压位。
*/#include <iostream>
#include <algorithm>using namespace std;typedef long long LL;const LL M=1e9;
const int N = 12e4+10;
int p[N],sum[N],num=0;
bool hs[N];void get_primer(int n) //获得n以内的素数
{for(int i=2;i<=n;++i){if(hs[i]==0)p[num++]=i;for(int j=0;p[j]<=n/i;++j){hs[p[j]*i]=1;if(i%p[j]==0)break;}}
}int cal(int n,int p) //计算n! 内有多少个因子p
{int cnt=0;while (n){cnt+=n/p;n/=p;}return cnt;
}void get_pow(int a,int b) //将C(a,b) 进行质因子分解
{for(int i=0;i<num;++i){int cnt=cal(a,p[i]) -cal(b,p[i])*2;sum[i]=cnt;}
}void mul(vector<LL> &A,int b) //高精度与整数相乘
{LL d=0;for(int i=0;i<A.size();++i){d+=A[i]*b;A[i]=d%M;d/=M;}while(d){A.push_back(d%M);d/=M;}
}void print(vector<LL> &res) //输出高精度数据
{cout << res.back(); // 最后一位不需要输出九位for(int i=res.size()-2;i>=0;--i)printf("%09lld",res[i]); //因为是按着1e9作为进制数,所以中间的位要输出九位,cout << endl; //不足九位的补零
}int main()
{int n;cin >> n;get_primer(2*n);get_pow(2*n,n);int r=n+1;for(int i=0;i<num && r>1;++i){while(r%p[i]==0){sum[i]-=1;r/=p[i];}}vector<LL> res;res.push_back(1);for(int i=0;i<num;++i)for(int j=0;j<sum[i];++j)mul(res,p[i]);print(res);return 0;}
卡特兰数(高精度乘法压位)相关推荐
- 2017 西安网络赛A Tree(树上静态查询,带权并查集,矩阵乘法压位,好题)
题目链接 题意: 给出 \(n(n \leq 3000)\) 个结点的一棵树,树上每个结点有一个 \(64 \times 64\) 的 \(0,1\)矩阵,每个结点上的矩阵是根据输入的 \(seed\ ...
- 高精度(压位+判负数+加减乘+读写)
本算法目前属于还处于测试状态,欢迎Hack! struct gj{bool fu; //是否是负数int tt,mod; //高精的长度int s[40005]; //压位用的数组inline gj( ...
- 刷题总结——regular words(hdu1502 dp+高精度加法+压位)
题目: Problem Description Consider words of length 3n over alphabet {A, B, C} . Denote the number of o ...
- 【算法专题】卡特兰数
卡特兰数 1. 概述 卡特兰数:首先这个一个数,很多问题的结果都是卡特兰数,比如2016年全国三卷数学选择题压轴题让求解的就是卡特兰数,问题如下: 首先是结论:卡特兰数为: C 2 n n n + 1 ...
- 【卡特兰数】【递推】【高精+压位】JZOJ·WZK打雪仗
题目大意: 有2n个人站在一个圈里,有n人要去其他n人那(两人连一线),线不能与其他线碰在一起,让你求方案数. 思路: 卡特兰数,因为位数大,高精压位走起 CodeCodeCode: #include ...
- (压位)高精度乘法,高精度加法,高精度减法,高精度除法
大值 乘 大值: 题目链接:https://www.luogu.com.cn/problem/P1303 # include <iostream> # include <string ...
- 火车进出栈问题(高精度+压位+质因数分解消除除法)
题目描述 一列火车n节车厢,依次编号为1,2,3,-,n. 每节车厢有两种运动方式,进栈与出栈,问n节车厢出栈的可能排列方式有多少种. 输入格式 输入一个整数n,代表火车的车厢数. 输出格式 输出一个 ...
- C语言丨小 学 数 学(二):高精度乘法
计算机,顾名思义,是用于高速计算的电子计算机器.它可以进行数值计算,又可以进行逻辑计算,还具有存储记忆功能.是能够按照程序运行,自动.高速处理海量数据的现代化智能电子设备.[1]但是,计算机的功能毕竟 ...
- 阿里云 超级码力在线编程大赛初赛 第2场 题目4. 小栖的金字塔(超级卡特兰数+除法求模/乘法逆元)
文章目录 1. 题目 2. 解题 1. 题目 来源:https://tianchi.aliyun.com/oj/15165469968503404/76745683739284070 2. 解题 按道 ...
最新文章
- powerdesigner箭头如何画_用Scratch编程画几何图形:如何画多边形
- mysql录入foreigen错误_Python MySQLdb 使用utf-8 编码插入中文数据问题
- oracle11g 读写,oracle11g pysical standby开启临时读写
- Java源码研究之object to json string debug
- 学弟学妹们,C#为啥学不好?是因为你没真正理解面向对象的思想!
- SGU 327 Yet Another Palindrome(状态压缩DP)
- 初中女生数学不好能学计算机,初中数学成绩好的女生,多半有这3个“小特征”,学渣可装不来!...
- C 标准库中输出到字符串、到文件的相关函数
- python绘制emoji_使用Emoji表情拼成汉字
- Linux版本怎么设置指纹,Linux上使用Thinkpad指纹识别
- svn的安装linux
- 微分方程的数值解法—四阶龙格库塔(Runge-Kutta)的推导以及C代码实现
- 机器学习之集成学习(七)随机森林scikit-learn库
- HTML Rendering Error
- Tableau的安装与下载
- Altera FPGA 配置方式概述
- vue中实现拖拽调整顺序功能
- idea的卸载和安装
- python卡方拟合优度检验_如何使用Python中的科学库执行卡方拟合优度检验?
- 联想笔记本电脑——开启BIOS/CPU开启虚拟化(VT)
热门文章
- 基于华为云自定义模板的图片操作演示(框选、拖拽)
- Richard Feynman 《What Do You Care What Other People Think》
- 多看阅读怎么横屏阅读
- Python学习笔记(一):Pandas中的数据结构
- 银行动态口令令牌工作原理简析
- 特征筛选———贷款用户是否会逾期
- 12V 电路采用电阻分压,产生 ± 6V的电源
- CT主控计算机功能不包括,​项目四 计算机X线体层成像 任务二CT的基本硬件与CT的类型...
- DAZ3D STUDIO 中复杂背景道具的简化
- 大数据任务调度工具 Apache DolphinScheduler