洛谷 - P1492 猩猩散步

点击打开链接

 

POINT：

 

答案为C（n+m,m）或者C(n+m,n)。

然后因为完全模拟为超时。JAVA大数乘法超时。

那么先把分母分子全部质因数分解。

把答案化为质因数的个数。

然后作高精度低乘

 

 

#include <iostream>
#include <stdio.h>
#include <algorithm>
#include <string.h>
#include <vector>
#include <queue>
using namespace std;
#define LL long long
const LL N = 50002*2;
LL p[N];
LL k=0;
LL isp[N];
void getprime()//质数筛
{
    for(LL i=2;i<N;i++) isp[i]=1;
    for(LL i=2;i<N;i++){
        if(isp[i]){
            p[++k]=i;
            for(LL j=2*i;j<N;j+=i){
                isp[j]=0;
            }
        }
    }
}
LL cnt[N];
 
void add(LL x)//分子的分解质因数
{
    for(LL i=1;i<=k&&p[i]<=x;i++){
        if(x%p[i]==0){
            cnt[i]++;
            x=x/p[i];
            while(x%p[i]==0)
            {
                cnt[i]++;
                x=x/p[i];
            }
        }
    }
}
 
void devide(LL x)//分母的分解质因数
{
    if(x==1) return;
    for(LL i=1;i<=k&&p[i]<=x;i++){
        if(x%p[i]==0){
            cnt[i]--;
            x=x/p[i];
            while(x%p[i]==0)
            {
                cnt[i]--;
                x=x/p[i];
            }
        }
    }
}
 
LL ans[22];
 
LL mod = 100000;
 
void mul(LL it)//高精度低乘
{
    LL jw=0;
    for(LL i=1;i<=20;++i)
    {
        ans[i]=ans[i]*it+jw;
        jw=ans[i]/mod;
        ans[i]%=mod;
    }
}
 
int main()
{
    getprime();
    LL m,n;
    scanf("%lld %lld",&m,&n);
    for(LL i=m+n;i>=m+1;i--){//分子
        add(i);
    }
    for(LL i=n;i>=2;i--){//分母
        devide(i);
    }
    ans[1]=1;
    for(LL i=k;i>=1;i--){//结果
        while(cnt[i]--){
            mul(p[i]);
        }
    }
    for(LL i=20;i>=1;i--){
        printf("%05lld",ans[i]);
        if(i&1) printf("\n");
    }
 
 
}