bzoj 3738: [Ontak2013]Kapitał

题意：

求$C_{n+m}^n$去掉所有末尾的0后对$10^k$取模的结果。

题解：

显然是扩展lucas的题了。
去掉0可以看作除$10^x$后在取余。
那么就将2,5的因子都提出来，最后再补回去。
其他就是普通的扩展lucas了。
code：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#define LL long long
using namespace std;
LL n,m,k,MOD=1LL,k2,k5,P,K;
LL pre[2][2000010];
LL pow(LL a,LL b,LL mod)
{
    LL ans=1;
    while(b)
    {
        if(b&1) ans=ans*a%mod;
        a=a*a%mod;b>>=1;
    }
    return ans;
}
LL inv(LL a,LL pi,LL pk)
{
    if(!a) return 1;
    return pow(a,pk-(pk/pi)-1,pk);
}
LL mul(LL n,LL pi,LL pk)
{
    if(!n) return 1;
    LL ans;
    if(pi==2) ans=pow(pre[0][pk],n/pk,pk)*pre[0][n%pk];
    if(pi==5) ans=pow(pre[1][pk],n/pk,pk)*pre[1][n%pk];
    return ans*mul(n/pi,pi,pk)%pk;
}
LL num(LL n,LL m,LL p)
{
    LL k=0;
    for(LL i=n;i;i/=p) k+=i/p;
    for(LL i=m;i;i/=p) k-=i/p;
    for(LL i=n-m;i;i/=p) k-=i/p;
    return k;
}
LL lucas(LL n,LL m,LL mod,LL pi,LL pk)
{
    if(m>n) return 0;
    LL a=mul(n,pi,pk),b=mul(m,pi,pk),c=mul(n-m,pi,pk);
    LL ans;
    if(pi==2) ans=a*inv(b,pi,pk)%pk*inv(c,pi,pk)%pk*inv(pow(5,k5,pk),pi,pk)%pk*pow(P,K,pk)%pk;
    if(pi==5) ans=a*inv(b,pi,pk)%pk*inv(c,pi,pk)%pk*inv(pow(2,k2,pk),pi,pk)%pk*pow(P,K,pk)%pk;
    return ans*(mod/pk)%mod*inv(mod/pk,pi,pk)%mod;
}
void print(LL ans)
{
    LL x=ans,p=0;
    while(x) p++,x/=10;
    for(int i=1;i<=k-p;i++) printf("0");
    printf("%lld",ans);
}
int main()
{
    scanf("%lld %lld %lld",&n,&m,&k);n+=m;
    LL ans=0;
    for(int i=1;i<=k;i++) MOD*=10LL;
    LL pk2=pow(2,k,MOD),pk5=pow(5,k,MOD);
    pre[0][0]=pre[1][0]=1;
    for(int i=1;i<=2000000;i++)
    {
        if(i%2) pre[0][i]=pre[0][i-1]*i%pk2;else pre[0][i]=pre[0][i-1];
        if(i%5) pre[1][i]=pre[1][i-1]*i%pk5;else pre[1][i]=pre[1][i-1];
    }
    k2=num(n,m,2);k5=num(n,m,5);
    if(k2>k5) P=2,K=k2-k5;
    else P=5,K=k5-k2;
    ans=(lucas(n,m,MOD,2,pow(2,k,MOD))+lucas(n,m,MOD,5,pow(5,k,MOD)))%MOD;
    print(ans);
}