bzoj3907 网格 & bzoj2822 [AHOI2012]树屋阶梯——卡特兰数+高精度
题目:bzoj3907:https://www.lydsy.com/JudgeOnline/problem.php?id=3907
bzoj2822:https://www.lydsy.com/JudgeOnline/problem.php?id=2822
bzoj3907:
从网格图的角度很好想啊,就是像定义一样,一下子就出来了 ans = C(n+m,n) - C(n+m,m-1);
那么从01串的角度呢?这个是0和1的个数不同的问题呢...
看到这篇博客:https://blog.csdn.net/u014097230/article/details/44244793
里面有很多很好的小例子!其中的 hdu 1133,不正可以转化成这道题吗!
于是就解决了,得到的答案和上面那种方法得到的一样:ans = C(n+m,n) - C(n+m,m-1);
然而难道是卡特兰数的题都太容易写了,所以非得来个高精度?!一时间不想做了...
直到看到了这篇博客:https://www.cnblogs.com/Tunix/p/4354348.html
怎么会有这么优美的高精度求组合数写法!!!顿时有了信心,马上抄了一遍成功写好了!
我会牢牢记住这篇美丽的高精度的(话说好像还是第一次写高精度求组合数啊)。
bzoj2822:就是求卡特兰数第 n 项,因为正好符合了那个递推式;
也要高精度,就用同一个板子A了(令 m = n 就好啦)。
代码如下:
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
int const maxn=,mod=1e8;
int n,m,pri[maxn],mn[maxn],cnt,p[maxn];
ll a[],b[];
void init()
{
mn[]=;
for(int i=;i<=n+m;i++)
{
if(!mn[i])pri[++cnt]=i,mn[i]=i;
for(int j=;j<=cnt&&i*pri[j]<=n+m;j++)
{
mn[i*pri[j]]=pri[j];
if(i%pri[j]==)break;
}
}
}
ll pw(int a,int b)
{
ll ret=;
for(;b;b>>=,a*=a)
if(b&)ret*=a;
return ret;
}
void mul(ll a[],ll k)
{
ll x=;
for(int i=;i<=a[];i++)
{
a[i]=a[i]*k+x;
x=a[i]/mod;
a[i]%=mod;
}
while(x)a[++a[]]+=x%mod,x/=mod;
}
void getc(ll a[],int n,int m)
{
memset(p,,sizeof p);
for(int i=;i<=n;i++)p[i]++;
for(int i=;i<=m;i++)p[i]--;
for(int i=;i<=n-m;i++)p[i]--;
for(int i=n;i>=;i--)
{
if(mn[i]==i) mul(a,pw(i,p[i]));
else p[mn[i]]+=p[i],p[i/mn[i]]+=p[i];//转化到质因子上!太妙了!
}
}
void dec(ll a[],ll b[])
{
for(int i=;i<=a[];i++)
{
if(a[i]<b[i])a[i]+=mod,a[i+]--;
a[i]-=b[i];
}
while(a[]&&!a[a[]])a[]--;
}
void print(ll a[])
{
printf("%lld",a[a[]]);
for(int i=a[]-;i;i--)
printf("%08lld",a[i]);
printf("\n");
}
int main()
{
scanf("%d%d",&n,&m);
init();
a[]=a[]=b[]=b[]=;
getc(a,n+m,n);
getc(b,n+m,m-);
dec(a,b);
print(a);
return ;
}
bzoj3907 网格 & bzoj2822 [AHOI2012]树屋阶梯——卡特兰数+高精度的更多相关文章
- BZOJ2822[AHOI2012]树屋阶梯——卡特兰数+高精度
题目描述 暑假期间,小龙报名了一个模拟野外生存作战训练班来锻炼体魄,训练的第一个晚上,教官就给他们出了个难题.由于地上露营湿气重,必须选择在高处的树屋露营.小龙分配的树屋建立在一颗高度为N+1尺(N为 ...
- BZOJ2822:[AHOI2012]树屋阶梯(卡特兰数,高精度)
Description 暑假期间,小龙报名了一个模拟野外生存作战训练班来锻炼体魄,训练的第一个晚上,教官就给他们出了个难题.由于地上露营湿气重,必须选择在高处的树屋露营.小龙分配的树屋建立在一颗高度为 ...
- bzoj 3907 网格 bzoj2822 [AHOI2012]树屋阶梯——卡特兰数(阶乘高精度模板)
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3907 https://www.lydsy.com/JudgeOnline/problem.p ...
- bzoj2822[AHOI2012]树屋阶梯(卡特兰数)
2822: [AHOI2012]树屋阶梯 Time Limit: 1 Sec Memory Limit: 128 MBSubmit: 879 Solved: 513[Submit][Status] ...
- [bzoj2822][AHOI2012]树屋阶梯 (卡特兰数+分解质因数+高精度)
Description 暑假期间,小龙报名了一个模拟野外生存作战训练班来锻炼体魄,训练的第一个晚上,教官就给他们出了个难题.由于地上露营湿气重,必须选择在高处的树屋露营.小龙分配的树屋建立在一颗高度为 ...
- BZOJ 2822: [AHOI2012]树屋阶梯 [Catalan数 高精度]
2822: [AHOI2012]树屋阶梯 Time Limit: 1 Sec Memory Limit: 128 MBSubmit: 779 Solved: 453[Submit][Status] ...
- bzoj 2822 [AHOI2012]树屋阶梯 卡特兰数
因为规定n层的阶梯只能用n块木板 那么就需要考虑,多出来的一块木板往哪里放 考虑往直角处放置新的木板 不管怎样,只有多的木板一直扩展到斜边表面,才会是合法的新状态,发现,这样之后,整个n层阶梯就被分成 ...
- P2532 [AHOI2012]树屋阶梯 卡特兰数
这个题是一个卡特兰数的裸题,为什么呢?因为可以通过划分来导出递推式从而判断是卡特兰数,然后直接上公式就行了.卡特兰数的公式见链接. https://www.luogu.org/problemnew/s ...
- 【BZOJ 2822】[AHOI2012]树屋阶梯 卡特兰数+高精
这道题随便弄几个数就发现是卡特兰数然而为什么是呢? 我们发现我们在增加一列时,如果这一个东西(那一列)他就一格,那么就是上一次的方案数,并没有任何改变,他占满了也是,然后他要是占两格呢,就是把原来的切 ...
随机推荐
- java虚拟机(八)--java性能监控与故障处理工具
问题定位: 除了个人经验,知识,工具也是很重要的,通过数据进行问题分析,包括:运行日志.异常堆栈.GC日志.线程快照(threaddump/javacore文件 ).堆转储快照(heapdump/hp ...
- 让元素div消失在视野中
让元素div消失在视野中1.position:absolute/relative/fixed + 方位 top/bottom/left/right: -9999px2.display:none3.vi ...
- Gym - 101550A(Artwork 倒序+并查集)
题目: 思路: 1.对输入数据离线,先把所有的黑线都画出来,统计一下剩余的白色连通块的个数,dfs过程将一个连通块放到一个集合中. 2.倒着往前消去黑线,如果当前的块A是白块就看他的四周有没有白块:有 ...
- Linux无人值守安装系统
yum install dhcp -yvim /etc/dhcp/dhcpd.conf--------------------------------------------- allow booti ...
- 59.关系型与document类型数据模型对比
现假设有如下两个类: class Department(object): def __init__(self, dept_id, name, desc, employees=[]): self.dep ...
- LINUX-光盘
cdrecord -v gracetime=2 dev=/dev/cdrom -eject blank=fast -force 清空一个可复写的光盘内容 mkisofs /dev/cdrom > ...
- python之cookbook-day01
第一章:数据结构和算法 1.1 解压序列赋值给多个变量 >>> p = (4, 5) >>> x, y = p >>> x 4 >>& ...
- reading/writing files in Python
file types: plaintext files, such as .txt .py Binary files, such as .docx, .pdf, iamges, spreadsheet ...
- java诗词横版--转为竖版
import java.util.Arrays; /* 诗句横版转成竖版输出 */public class PrintPoem { public static void main(String[] a ...
- phpcms 短信替换
后台表单向导文件路径: [/www/wwwroot/phpcms/phpcms/modules/formguide/templates/formguide_info_list.tpl.php] pub ...