bzoj1300: [LLH邀请赛]大数计算器

传送门
最后的十二位我们大力质因数分解之后快速幂爆出来。
前面已看只有三位于是考虑大力用log来爆过去。
我们计算出答案的log值，然后截取ans-int(ans)+2的log值。
然后我们大力pow一波然后取整爆过去你就AC辣。

#include<bits/stdc++.h>
#define mo 1000000000000ll
#define ll long long
#define N 1000005
using namespace std;
int n,m,tot;
int num[N],pri[N],fl[N],mn[N];
ll ans=1;
double dig;
ll mul(ll x,ll y){
    ll s=0;
    for (;y;y/=2,x=(x+x)%mo)
        if (y&1) s=(s+x)%mo;
    return s; 
}
ll Pow(ll x,ll y){
    ll s=1;
    for (;y;y/=2,x=mul(x,x))
        if (y&1) s=mul(s,x);
    return s;
}
void add(int x,int f){
    for (;x^1;x/=pri[mn[x]]) num[mn[x]]+=f;
}
void pre(int n){
    for (int i=2;i<=n;i++){
        if (!fl[i]) pri[++tot]=i,mn[i]=tot;
        for (int j=1;j<=tot&&i*pri[j]<=n;j++){
            fl[i*pri[j]]=1;
            mn[i*pri[j]]=j;
            if (i%pri[j]==0) break;
        }
    }
}
int main(){
    scanf("%d%d",&n,&m);
    pre(n);
    for (int i=m+1;i<=n;i++) add(i,1);
    for (int i=1;i<=n-m;i++) add(i,-1);
    for (int i=1;i<=tot;i++)
        if (num[i]){
            ans=mul(ans,Pow(pri[i],num[i]));
            dig+=log10(pri[i])*num[i];
        }
    if (dig<12) printf("%lld",ans);
    else{
        long double ans2=pow(10,dig-int(dig)+2);
        printf("%lld...%lld",(ll)ans2,ans%1000000000);
    }
}