【BZOJ 3907】网格(Catalan数)
题目链接
这个题推导公式跟\(Catalan\)数是一样的,可得解为\(C_{n+m}^n-C_{n+m}^{n+1}\)
然后套组合数公式\(C_n^m=\frac{n!}{m!(n-m)!}\)
用阶乘分解的方法对分子和分母分解质因数然后指数相减,最后把剩下的高精度乘起来就行了,这样就避免了高精除法。可以用快速幂,但我太lan了,就直接暴力乘起来了。
说一下怎么阶乘分解,直接对每个数分解质因数的时间复杂度是\(O(n\sqrt{n})\),这显然是不可忍受的。
于是,考虑先用线筛求出\(1~n\)之间所有质数,然后枚举所有质数\(p\),\(1~n\)中所有\(p\)的倍数包含一个质因子\(p\),所有\(p^2\)的倍数包含一个质因子\(p\),...
所以,\(n!\)中\(p\)的指数为\(\sum_{i=1}^{\left\lfloor \log_p n\right\rfloor}\left\lfloor n/p^i \right\rfloor\)
枚举\(p\)求就行了,时间复杂度\(O(n\ log\ n)\)。
Code:
#include <cstdio>
#include <cstring>
const int MAXN = 10010;
int prime[MAXN], v[MAXN], c[MAXN], d[MAXN];
int n, cnt, m;
const int MOD = 10000;
struct Int{
int s[MAXN];
Int(int x){ memset(s, 0, sizeof s); s[++s[0]] = x % 10, x /= 10; }
void mul(int x){
s[1] *= x;
for(int i = 2; i <= s[0]; ++i){
s[i] = s[i] * x + s[i - 1] / MOD;
s[i - 1] %= MOD;
}
while(s[s[0]] >= MOD){
s[++s[0]] = s[s[0] - 1] / MOD;
s[s[0] - 1] %= MOD;
}
}
void cut(Int x){
for(int i = 1; i <= x.s[0]; ++i)
s[i] -= x.s[i];
for(int i = 1; i <= x.s[0]; ++i){
if(s[i] < 0){
s[i] += MOD;
--s[i + 1];
}
}
if(!s[s[0]]) --s[0];
}
void print(){
printf("%d", s[s[0]]);
for(int i = s[0] - 1; i; --i)
printf("%04d", s[i]);
}
}a(1), b(1);
int main(){
scanf("%d%d", &n, &m);
for(int i = 2; i <= n + m; ++i){
if(!v[i]){
v[i] = i;
prime[++cnt] = i;
}
for(int j = 1; j <= cnt; ++j){
if(prime[j] > v[i] || prime[j] * i > n + m) break;
v[prime[j] * i] = prime[j];
}
}
for(int i = 1; i <= cnt; ++i)
for(int j = prime[i]; j <= n + m; j *= prime[i])
c[i] += (n + m) / j;
for(int i = 1; i <= cnt && prime[i] <= n; ++i)
for(int j = prime[i]; j <= n; j *= prime[i])
c[i] -= n / j;
for(int i = 1; i <= cnt && prime[i] <= m; ++i)
for(int j = prime[i]; j <= m; j *= prime[i])
c[i] -= m / j;
for(int i = 1; i <= cnt; ++i)
for(int j = 1; j <= c[i]; ++j)
a.mul(prime[i]);
for(int i = 1; i <= cnt; ++i)
for(int j = prime[i]; j <= n + m; j *= prime[i])
d[i] += (n + m) / j;
for(int i = 1; i <= cnt && prime[i] <= n + 1; ++i)
for(int j = prime[i]; j <= n + 1; j *= prime[i])
d[i] -= (n + 1) / j;
for(int i = 1; i <= cnt && prime[i] <= m - 1; ++i)
for(int j = prime[i]; j <= m - 1; j *= prime[i])
d[i] -= (m - 1) / j;
for(int i = 1; i <= cnt; ++i)
for(int j = 1; j <= d[i]; ++j)
b.mul(prime[i]);
a.cut(b);
a.print();
return 0;
}
【BZOJ 3907】网格(Catalan数)的更多相关文章
- BZOJ 3907: 网格 [Catalan数 高精度]
3907: 网格 Time Limit: 1 Sec Memory Limit: 256 MBSubmit: 402 Solved: 180[Submit][Status][Discuss] De ...
- BZOJ 3907: 网格
Description 求不跨过直线 \(y=x\) ,到达 \((n,m)\) 的方案数. Sol 组合数学+高精度. 这个推导过程跟 \(Catalan\) 数是一样的. 答案就是 \(C^{n+ ...
- bzoj 3907: 网格 组合数学
3907: 网格 Time Limit: 1 Sec Memory Limit: 256 MBSubmit: 13 Solved: 7[Submit][Status][Discuss] Descr ...
- BZOJ 3907: 网格( 组合数 + 高精度 )
(0,0)->(n,m)方案数为C(n,n+m), 然后减去不合法的方案. 作(n,m)关于y=x+1的对称点(m-1,n+1), 则(0,0)->(m-1,n+1)的任意一条路径都对应( ...
- bzoj 3907 网格 bzoj2822 [AHOI2012]树屋阶梯——卡特兰数(阶乘高精度模板)
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3907 https://www.lydsy.com/JudgeOnline/problem.p ...
- BZOJ 3907: 网格【组合数学】
Description 某城市的街道呈网格状,左下角坐标为A(0, 0),右上角坐标为B(n, m),其中n >= m.现在从A(0, 0)点出发,只能沿着街道向正右方或者正上方行走,且不能经过 ...
- [Catalan数三连]网格&有趣的数列&树屋阶梯
如何让孩子爱上打表 Catalan数 Catalan数是组合数学中一个常出现在各种计数问题中的数列. 以比利时的数学家欧仁·查理·卡塔兰 (1814–1894)的名字来命名. 先丢个公式(设第n项为$ ...
- BZOJ 1485: [HNOI2009]有趣的数列 [Catalan数 质因子分解]
1485: [HNOI2009]有趣的数列 Description 我们称一个长度为2n的数列是有趣的,当且仅当该数列满足以下三个条件: (1)它是从1到2n共2n个整数的一个排列{ai}: (2)所 ...
- BZOJ 1856: [Scoi2010]字符串 [Catalan数]
1856: [Scoi2010]字符串 Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 1418 Solved: 790[Submit][Status][ ...
随机推荐
- Java Set集合(HashSet、TreeSet)
什么是HashSet?操作过程是怎么样的? 1.HashSet底层实际上是一个HashMap,HashMap底层采用了哈希表数据结构 2.哈希表又叫做散列表,哈希表底层是一个数组,这个数组中每一个元素 ...
- [POJ3585]Accumulation Degree
题面 \(\text{Solution:}\) 有些题目不仅让我们做树型 \(\text{dp}\) ,而且还让我们换每个根分别做一次, 然后这样就愉快的 \(\text{TLE}\) 了,所以我们要 ...
- 使用vue和web3创建你的第一个以太坊APP
欢迎回到这个很牛的教程系列的第2部分,在教程中我们亲手构建我们的第一个分布式应用程序. 在第二部分中,我们将介绍VueJS和Vuex的核心概念,并引入web3js以与metamask进行交互. 如果你 ...
- C++内置类型如何存放于计算机内存中
摘要:内置类型的机器实现.字/字节/比特.内存 一.概念 计算机以比特序列存储数据,每个比特非0即1,如:00011011011100010110010000111011... 二.计算机以块来处理内 ...
- Tomcat 中如何给 web 项目配置虚拟目录的方法
为什么要给 web 项目配置虚拟目录? 初学 JavaWeb 时,会发现只要我们把 web 项目放到 Tomcat 的 webapps 目录下,再通过 http://localhost:8080/项目 ...
- css如何选择相同class下的第一个class元素和最后一个元素?
如图,如果像选择类名为 class="exerciseInfo" 中的第一个和最后一个div,做法如下: 选择第一个类名: .exerciseInfo: nth-of-type ...
- php数据缓存到文件类设计
// 自定义缓存类 class Cache_Filesystem { // 缓存写保存 function set ($key, $data, $ttl) { //打开文件为读/写模式 $h = fop ...
- kudu介绍及安装配置
kudu介绍及安装配置 介绍 Kudu 是一个针对 Apache Hadoop 平台而开发的列式存储管理器.Kudu 共享 Hadoop 生态系统应用的常见技术特性: 它在 commodity har ...
- COGS1752. [BOI2007]摩基亚Mokia CDQ
CDQ的板子题 #include<cstdio> #include<cstring> #include<iostream> #include<algorith ...
- COGS 930. [河南省队2012] 找第k小的数 主席树
主席树裸板子 #include<cstdio> #include<iostream> #include<algorithm> #define MAXN 100005 ...