BZOJ3738 : [Ontak2013]Kapitał
$C_{N+M}^N=\frac{(N+M)!}{N!M!}$
考虑求出$ans\bmod 10^9$的值
$10^9=2^9\times5^9$
以$2^9$为例,先预处理出$1$..$2^9$中不是2的倍数的数的前缀积s[],显然$n!\bmod 2^9$有着长度为$2^9$的循环节
将答案表示成$a\times2^b$的形式,$a$与$2^9$互质,可以直接逆元,b直接相减即可
cal(n).a=s[n%512]*pow(s[512],n/512)*cal(n/2).a
cal(n).b=n/2+cal(n/2).b
如此递归计算即可
答案中末尾0的个数为min(2的个数,5的个数)
以$2^9$为例,除以10相当于乘上5的逆元,同时2的个数减1
分别算出答案后再用中国剩余定理合并即可
#include<cstdio>
typedef long long ll;
ll n,m,k,x,y,P,B,s[2000000],res[2],del,ans,i,T=1;
ll exgcd(ll a,ll b){
if(!b)return x=1,y=0,a;
ll d=exgcd(b,a%b),t=x;
return x=y,y=t-a/b*y,d;
}
ll rev(ll a,ll P){exgcd(a,P);while(x<0)x+=P;return x%P;}
ll pow(ll a,ll b,ll P){ll t=1;for(;b;b>>=1LL,a=a*a%P)if(b&1LL)t=t*a%P;return t;}
struct Num{
ll a,b;
Num(){a=1,b=0;}
Num(ll _a,ll _b){a=_a,b=_b;}
Num operator*(Num x){return Num(a*x.a%P,b+x.b);}
Num operator/(Num x){return Num(a*rev(x.a,P)%P,b-x.b);}
}now[2];
Num cal(ll n){return n?Num(s[n%P]*pow(s[P],n/P,P)%P,n/B)*cal(n/B):Num(1,0);}
void pre(){for(i=s[0]=1;i<P;i++)if(i%B)s[i]=s[i-1]*i%P;else s[i]=s[i-1];s[P]=s[P-1];}
int main(){
scanf("%lld%lld%lld",&n,&m,&k);
while(k--)T*=10;
B=2,P=512,pre();
now[0]=cal(n+m)/cal(n)/cal(m);
del=now[0].b;
B=5,P=1953125,pre();
now[1]=cal(n+m)/cal(n)/cal(m);
if(del>now[1].b)del=now[1].b;
while(del--)P=512,now[0]=now[0]/Num(5,1),P=1953125,now[1]=now[1]/Num(2,1);
B=2,P=512,res[0]=now[0].a*pow(B,now[0].b,P)%P;
B=5,P=1953125,res[1]=now[1].a*pow(B,now[1].b,P)%P;
ans=(1953125LL*rev(1953125,512)%T*res[0]%T+512LL*rev(512,1953125)%T*res[1]%T)%T;
while(ans*10<T)putchar('0'),T/=10;
return printf("%lld",ans),0;
}
BZOJ3738 : [Ontak2013]Kapitał的更多相关文章
- BZOJ3738 [Ontak2013]Kapitał 【扩展Lucas】
题目链接 BZOJ3738 题解 复习 同上 但是为了消去因子\(10\),处理\(2^k\)的时候,乘回\(2^{k_1}\)时,应同时计算\(5^{k_2}\) 如果\(k_1 \ge k_2\) ...
- bzoj AC倒序
Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...
- BZOJ3734 : [Ontak2013]Miny
将所有炸弹按坐标排序 x<-y连边表示x爆炸了y也会爆炸 如果是DAG则直接拓扑排序+DP求出每个点出发能走到的最左端和最右端的点 有环则SCC缩点后再拓扑 用线段树优化建图的过程 边数$O(n ...
随机推荐
- unity3d的四元数 Quaternion
原地址:http://www.cnblogs.com/88999660/archive/2013/04/02/2995074.html 今天准备学习和研究下unity3d的四元数 Quaternion ...
- C++模板常用使用方法介绍
转载:http://developer.51cto.com/art/201002/182202.htm C++编程语言中的模板应用在一定程度上大大提高了程序开发的效率.我们在这篇文章中为大家详细讲解一 ...
- Coursera台大机器学习课程笔记3 – 机器学习的分类和机器学习的可能性
第三讲比较简单,参考:http://www.cnblogs.com/HappyAngel/p/3466527.html 第四讲很抽象,尤其是第四个视频,目的仍然是为了证明机器学习是可能的,不过这个博主 ...
- Python 正则表达式_re模块_使用compile加速
使用compile加速 compile( rule [,flag] ) 将正则规则编译成一个Pattern对象,以供接下来使用. 第一个参数是规则式,第二个参数是规则选项. 返回一个Pattern对象 ...
- 用Python操纵MySQL
本例用Python操纵MySQL,从指定文件读取数据,并对数据进行处理,处理之后批量插入MySQL. 贴上代码: # -*- coding: gbk -*- import re import MySQ ...
- Redis和Memcache的区别分析
1. Redis中,并不是所有的数据都一直存储在内存中的,这是和Memcached相比一个最大的区别. 2. Redis不仅仅支持简单的k/v类型的数据,同时还提供list,set,hash等数据结构 ...
- linux安装setup工具
如果你的Linux系统是最小化安装的,可能会没有setup命令工具,环境是centos 5.8 安装setup命令工具的步骤. 安装setuptool #yum install setuptool 系 ...
- 在Python脚本中判断Python的版本
引自:http://segmentfault.com/q/1010000000127878 如果是给人读,用 sys.version,如果是给机器比较,用 sys.version_info,如果是判断 ...
- 【转】windows下安装和调用curl的方法
本文转自:http://1316478764.iteye.com/blog/2100778 curl是利用URL语法在命令行方式下工作的开源文件传输工具.它支持很多协议:FTP, FTPS, HTTP ...
- Android 和iOS 中关于View 的一点知识
View的概念和方法十分重要,这里将对Android 和iOS中出现的,关于视图的一些知识点进行总结,预计文章会比较长,要许多时间慢慢补充. 先转载一部分资料,感谢原作者! 原链接为:http://b ...