【bzoj1300】大数计算器

题意：

求C(n,m) 如果C(n,m)的位数<=12 那么直接输出 否则按XXX...XXXXXXXXX的形式输出

题解：

这题之前打过 但是昨天又想出一种新的做法 先说下新的做法吧- -

________________（坑爹的博客园出现了一些bug 不能显示回车 so 我是华丽的分割线）________________

首先要知道位数 显然可以用res存 如果是乘的res就加log(10,x) 除的就减掉

如果res<13那么暴力算 直接输出就完了

如果res>=13 就要知道前3位 和答案模10^10的值

前3位很简单 把res的整数部分边成2 10^res的整数部分就是前3位

后面要算C(n,m)%10^10的值 我们发现10^10=(2^10)*(5^10)

就把要乘或除的数a 转换成a=a'*(2^x)*(5^y) (a',10)=1

把2和5的次方分离出来算 那么就能保证运算的数都和10^10互质 就能用乘法逆元了

_________________________________________________________________________________

但是这有个问题- -

10^10*10^10可能会爆long long

对这问题有两种办法解决

1.高精 不解释

2.分别算模2^10 和模5^10 最后线性模方程算答案

鉴于这两种方法都听难打 而且bzoj还把这题变成高富帅题 于是我放弃治疗了- -

_________________________________________________________________________________

正解：

我们可以把每个数都分解质因数 因为n和m都比较小 所以可以把每个质数的指数都存下来 最后统一算就完了

_________________________________________________________________________________

代码：

 #include <cstdio>
 #include <cmath>
 using namespace std;
 typedef long long ll;
 const ll mo=1000000000000ll;
 ll n,m,ans=,num,sum[],bo[],primer[],ns;
 double save;
 ll power(ll a,ll b){
             ll res=;
             for (;b;b>>=){
                 if (b&) res=res*a%mo;
                 a=a*a%mo;
             }
             return res;
 }
 void makepr(ll t){
      for (ll i=;i<=t;i++){
          if (!bo[i]) primer[++ns]=i;
          for (ll j=;j<=ns && i*primer[j]<=t;j++){
              bo[i*primer[j]]=;
              if (i%primer[j]==) break;
          }
      }
 }
 int main(){
     freopen("bz1300.in","r",stdin);
     freopen("bz1300.out","w",stdout);
     ll i,j;
     scanf("%I64d%I64d\n",&n,&m);
     makepr(n);
     for (i=;i<=ns;i++){
         for (j=primer[i];j<=n;j*=primer[i]) sum[i]+=n/j;
         for (j=primer[i];j<=m;j*=primer[i]) sum[i]-=m/j;
         for (j=primer[i];j<=n-m;j*=primer[i]) sum[i]-=(n-m)/j;
     }
     for (i=;i<=ns;i++)
     if (sum[i]){
                 ans=(ans*power(primer[i],sum[i]))%mo;
                 save+=log10(primer[i])*sum[i];
     }
     if (save<) printf("%I64d",ans);
     else{
          long double ans2=pow(,save-static_cast<int>(std::floor(save))+);
          printf("%I64d...%I64d",static_cast<ll>(floor(ans2)),ans%(mo/));
     }
     fclose(stdin);
     fclose(stdout);
 }