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");
}
}