求大组合数

众所周知 ,所以我们只需要求出m!和n!(m-n)!。

但是由于n和m都偏大，直接乘一定会爆掉（喜闻乐见(●'◡'●)）。所以我们首先想到的是边乘边模，但是尝试了一组数据后我们神奇地发现答案不对。

但是我们知道另一种定理，n!%p=n!边乘边模的逆元。

所以这时候我们可以利用逆元来解决这个问题。

而逆元可以运用扩展欧几里得的方法求得。

下面我们给出欧几里得算法的证明:

首先我们设k为gcd(a,b)，则a=km，b=kn。

则a%b=a-c*b=km-c*kn=(m-cn)k

gcd(b,a%b)=gcd(kn,(m-cn)k)

由于k为a，b的最大公约数，所以n与m-cn互质，所以gcd(b,a%b)=gcd(a,b)=k。

知道了欧几里得算法，我们可以对其进行延伸得到了扩展欧几里得算法：

对于任意两个互质的a,b，总有gcd(a,b)=ax+by,扩展欧几里得算法可以用来求解

x，y。求法：

    根据欧几里得算法可知gcd(a,b)=gcd(b,a%b)。

    则bx‘+(a%b)y'=gcd(a,b)

    将a%b=a-(a/b)*b带入得

    ay'+b(x'-(a/b)*y')=gcd(a,b)

    对于a，b而言，他们对应的x,y则分别为y',(x'-(a/b)*y')。

而当b=0时，a*1+b*0=gcd(a,b)。

知道了这些后，此题就变得非常简单了，但是仍有需要注意的地方。

一：要防止逆元为负数的情况发生。

二：要注意因子中含有p的倍数，这样变成边模后结果就会变成0.

对于第一种情况，我们很好解决。对于第二种情况，我们可以先记录因子可以整除p几次后记录，将能整除的除去，剩下的就乘进去取模。

如果分母中的因子可以整除p的次数和与分子中的相等的话，就照原来的方法输出，否则输出0。

------------------------------------------------------------------------------------------------------------------------------------------------------------------------

代码如下

#include<iostream>
#include<cstdio>
#define LL long long
using namespace std;
LL extgcd(LL a,LL b,LL& x,LL& y)
{
    if(b!=0)
    {
        LL d=extgcd(b,a%b,y,x);
        y-=(a/b)*x;
        return d;
    }
    else
    {
       x=1,y=0;
       return a;
    }
}
int main()
{
    LL m,n,p;
    cin>>m>>n>>p;
    LL a=1,b=1,c=1;
    LL ka=0,kb=0,kc=0;
    for(int i=1;i<=m;i++)
    {
        int temp=i;
         while(temp%p==0)
          {
            temp/=p;
            ka++;
          }
          a=(a*temp)%p;
    }
     for(int i=1;i<=n;i++)
     {
        int temp=i;
         while(temp%p==0)
          {
            temp/=p;
            kb++;
          }
          b=(b*temp)%p;
     }
     for(int i=1;i<=m-n;i++)
     {
        int temp=i;
         while(temp%p==0)
          {
            temp/=p;
              kc++;
          }
          c=(c*temp)%p;
     }
    LL x1,x2,y;
    extgcd(b,p,x1,y);
    x1=(x1%p+p)%p; 
    extgcd(c,p,x2,y);
    x2=(x2%p+p)%p;
    a=((a*x1)%p)*x2%p;
    if(ka==kb+kc) cout<<a;
    else cout<<0;
}