$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ł的更多相关文章

  1. BZOJ3738 [Ontak2013]Kapitał 【扩展Lucas】

    题目链接 BZOJ3738 题解 复习 同上 但是为了消去因子\(10\),处理\(2^k\)的时候,乘回\(2^{k_1}\)时,应同时计算\(5^{k_2}\) 如果\(k_1 \ge k_2\) ...

  2. bzoj AC倒序

    Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...

  3. BZOJ3734 : [Ontak2013]Miny

    将所有炸弹按坐标排序 x<-y连边表示x爆炸了y也会爆炸 如果是DAG则直接拓扑排序+DP求出每个点出发能走到的最左端和最右端的点 有环则SCC缩点后再拓扑 用线段树优化建图的过程 边数$O(n ...

随机推荐

  1. BZOJ 1455

    STL的基本用法 (居然能空间卡过去= =!!!) #include <cstdio> #include <ext/pb_ds/priority_queue.hpp> #inc ...

  2. intellij idea 如何更改编辑器文本字体和大小

    换上了intellij idea之后,第一件事就是想要改变下文字字体,因为在我这个27寸的2k分辨率的屏幕上,文字显然太小了. intellij idea字体设值分成两部分,一部分是UI部分字体字号设 ...

  3. 【转】基于LDA的Topic Model变形

    转载自wentingtu 基于LDA的Topic Model变形最近几年来,随着LDA的产生和发展,涌现出了一批搞Topic Model的牛人.我主要关注了下面这位大牛和他的学生:David M. B ...

  4. static_cast dynamic_cast const_cast reinterpret_cast总结对比

    [本文链接] http://www.cnblogs.com/hellogiser/p/static_cast-dynamic_cast-const_cast-reinterpret_cast.html ...

  5. 29.调整数组顺序使奇数位于偶数前面[ReOrderArray]

    [题目] 输入一个整数数组,调整数组中数字的顺序,使得所有奇数位于数组的前半部分,所有偶数位于数组的后半部分.要求时间复杂度为O(n). [分析] 如果不考虑时间复杂度,最简单的思路应该是从头扫描这个 ...

  6. [Android Pro] Android 之使用LocalBroadcastManager解决BroadcastReceiver安全问题

    参考博客: http://blog.csdn.net/t12x3456/article/details/9256609 http://blog.csdn.net/lihenair/article/de ...

  7. JUC回顾之-ScheduledThreadPoolExecutor底层实现原理和应用

    项目中经常使用定时器,比如每隔一段时间清理下线过期的F码,或者应用timer定期查询MQ在数据库的配置,根据不同version实现配置的实时更新等等.但是timer是存在一些缺陷的,因为Timer在执 ...

  8. WMI

    https://wiki.jenkins-ci.org/display/JENKINS/Windows+slaves+fail+to+start+via+DCOM#Windowsslavesfailt ...

  9. wp8 入门到精通 ---时间

    DateTime.Now.ToShortTimeString()DateTime dt = DateTime.Now;dt.ToString();//2005-11-5 13:21:25dt.ToFi ...

  10. Laravel 4 系列入门教程(一)

    默认条件 本文默认你已经有配置完善的PHP+MySQL运行环境,懂得PHP网站运行的基础知识.跟随本教程走完一遍,你将会得到一个基础的包含登录的简单blog系统,并将学会如何使用一些强大的Larave ...