果然以前不想搞的东西,今天他妈全来了,我要爆炸,除了说操。。。。真是欲哭无泪啊。。。。。
//这道题目卡在逆元了。。。。
//利用逆元计算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));
}