70分逆元求助
查看原帖
70分逆元求助
282080
RefreshinglyNaive楼主2020/7/20 12:59
#include <bits/stdc++.h>
#define maxn 2000005
#define mod 20100403
using namespace std;
long long a[maxn],inv[maxn];
void exgcd(long long a,long long b,long long &x,long long &y)
{
	if(b==0) x=1,y=0;
	else exgcd(b,a%b,y,x), y=y-(a/b)*x;
}
long long C(long long x,long long y)
{
	long long xx,yy,x1,x2;
	exgcd(a[y],mod,xx,yy); x1=(xx%mod+mod)%mod;
	exgcd(a[x-y],mod,xx,yy); x2=(xx%mod+mod)%mod;
	return a[x]*x1%mod*x2%mod;
}
int main()
{
	int n,m;
	while(cin>>n>>m)
	{
		a[0]=1;
		for(int i=1;i<=n+m;i++)
			a[i]=a[i-1]*i%mod;
		cout<<C(n+m,m)-C(n+m,m-1)<<endl;
	}
	return 0;
} 

求助为什么70分

2020/7/20 12:59
加载中...