组合数C(M,N)

求一个数的组合数C(M,N)
C(M,N)=12… * m/n * (n-1) * (n-2) *… *(n-m+1)
意思就是C(3,10)=1 * 2 * 3/8 *9 *10
当M和N很大的时候就会爆精度，这个时候需要用到去除质因子来得到正确的答案。
先将分母里面的质因子数量记住，然后拿分子里面质因子数量来抵消，最后把剩下的乘起来就行了。

#include<stdio.h>
#include<string.h>
#include<math.h>
#include<iostream>
#include<algorithm>
#include<vector>
#include<map>
#include<queue>
using namespace std;
typedef long long ll;
int a[200005],b[200005],c[200005];
ll kpow(int a,int n)//快速幂 
{
 ll ans=1;
 while(n)
 {
  if(n%2==1)ans*=a;
  a*=a;
  n/=2;
 }
 return ans;
}
int main()
{
    int m,n;
    scanf("%d%d",&m,&n);
    int k=0;
    for(int i=2;i<=n;i++)//素筛求m里的素数，存到a数组中
    {
        if(b[i]==0)
        {
            a[k++]=i;
            for(int j=i+i;j<=n;j+=i)
                b[j]=1;
        }
    }
    for(int i=n;i>=n-m+1;i--)//分母拆分质因子
    {
        int co=i;
        for(int j=0;j<k;j++)
        {
            while(co%a[j]==0)
            {
                co/=a[j];
                c[j]++;//存入数量 
            }
            if(co==1)
                break;
        }
    }
    for(int i=2;i<=m;i++)//分字拆分质因子
    {
        int co=i;
        for(int j=0;j<k;j++)
        {
            while(co%a[j]==0)
            {
                co/=a[j];
                c[j]--;//相互抵消 
            }
            if(co==1)
                break;
        }
    }
    ll ans=1;
    for(int i=0;i<k;i++)
    {
        if(c[i]!=0)
         ans=ans*kpow(a[i],c[i]);
    }
    printf("%lld\n",ans);
}