组合数取模

最新推荐文章于 2021-02-28 11:55:45 发布

甦萌

最新推荐文章于 2021-02-28 11:55:45 发布

阅读量218

点赞数

分类专栏：组合数文章标签：组合数取模

本文链接：https://blog.csdn.net/lanyanzhiji123asd/article/details/87357274

版权

组合数专栏收录该内容

1 篇文章 0 订阅

订阅专栏

刷了几道组合数取模的题，模板都是一样的 Lucas定理

关于Lucas定理 https://blog.csdn.net/cai_haiq/article/details/75954298 这个博客讲得很好

下面是模板

#include<bits/stdc++.h>
using namespace std;
#define ll long long int
const int INF = 0x3f3f3f3f;
const ll p = 1e9+7;
ll nn,mm,kk;

ll qpow(ll n, ll m, ll p)
{
	ll res=1;
	while(m)
	{
		if(m%2!=0)
		{
			res*=n;
			res%=p;
		}
		n*=n;
		n%=p;
		m/=2;
	}
	return res;
}

ll C(ll n,ll m,ll p)
{
	if(m>n) return 0;
	ll a=1,b=1;
	while(m)
	{
		a*=n;
		a%=p;
		n--; 
		b*=m;
		m--;
		b%=p;
	}
	return a*qpow(b,p-2,p)%p;
}

ll Lucas(ll n,ll m,ll p)
{
	if(m==0)
		return 1;
	return Lucas(n/p,m/p,p)*C(n%p,m%p,p)%p;
	
	
}

int main()
{
	std::ios::sync_with_stdio(false);
	int t;
	while(cin>>nn>>mm)
	{
		
		if(nn<mm)swap(nn,mm);

			ll b=Lucas(nn+mm-4,nn-2,p)%p;
			cout<<b<<endl;

		
	} 

	return 0;
}