hdoj5698

果然以前不想搞的东西,今天他妈全来了,我要爆炸,除了说操。。。。真是欲哭无泪啊。。。。。

//这道题目卡在逆元了。。。。

//利用逆元计算1/(n!(m-n)!)
//对于正整数a,m如果有ax≡1(modm),那么把这个同余方程中x的最小正整数解叫做a模m的逆元。

problem1:
//这道题为什么要有乘法逆元呢?

模运算与基本四则运算有些相似,但是除法例外。其规则如下:
若a≡b (% p),则对于任意的c,都有(a + c) ≡ (b + c) (%p);

若a≡b (% p),则对于任意的c,都有(a * c) ≡ (b * c) (%p);

若a≡b (% p),c≡d (% p),则 (a + c) ≡ (b + d) (%p),(a - c) ≡ (b - d) (%p),

(a * c) ≡ (b * d) (%p),(a / c) ≡ (b / d) (%p);

problem2:
怎么求逆元:
//逆元一般用【扩展欧几里得】算法来求得,但是这里的m是素数,那么还可以根据【费马小定理】得到逆元为 a^m-2 mod m。

#include<cstdio>
#include<algorithm>
typedef long long LL;
using namespace std;
const int maxn = 1e6+7;
const int mod = 1e9+7;
long long fac[maxn];


long long qpow(long long a,long long b)
{
    long long ans=1;
    a%=mod;
    while(b)
    {
        if(b%2)
            ans=ans*a%mod;
        a=a*a%mod;
        b/=2;
    }
    return ans;
}

long long C(long long n,long long m)
{
    if(m>n||m<0)return 0;
    long long s1=fac[n];
    long long s2=fac[n-m]*fac[m]%mod;
  //  printf("%lld %lld
",s1,s2);
    return s1*qpow(s2,mod-2)%mod;
}

int n,m;
int main()
{
    fac[0]=1;
    for(int i=1;i<maxn;i++)
        fac[i]=fac[i-1]*i%mod;
    LL a,b;
    while(scanf("%d%d",&n,&m)!=EOF)
        printf("%lld
",C(n+m-4,n-2));
}
原文地址:https://www.cnblogs.com/keyboarder-zsq/p/5934524.html