【SCOI2010】【数论】【扩展欧几里徳】字符串

这是一道很有思维难度的题。

首先将模型就行变换，假设一个人在一个矩形中，站在(0,0)这个点上，只能向右走或向上走，走到(n,m)这个点的方案数，就为所有生成字符串的方案数(可以想象成向右走就取了一个1，向上走就取了一个0)。

得到总方案数以后，只要我们求出所有不合法的方案数即可。不难发现，所有不合法的方案一定穿过了y = x这条直线，所以我们做(0,0)关于y = x + 1的对称点(-1,1),这样从(-1,1)到(n,m)的方案数，即为所有不合法的方案数，减去这一部分可以了，而从一个点到另一个点的方案数可以使用组合数来求。

推导过程为:

考虑到模数为一个质数，所以我们可以使用乘法逆元的方法来加速(暴力也可以求)，这样这道题就完美地解决了。

代码:

#include<cstdio>
#include<cstring>
using namespace std;
const int mo = 20100403;
long long ans,n,m;
long long x,y;
void init()
{
	freopen("bzoj1856.in","r",stdin);
	freopen("bzoj1856.out","w",stdout);
}

void readdata()
{
	scanf("%lld%lld",&n,&m);
}

void extend_gcd(long long a,long long b,long long &x,long long &y)
{
	if(!b)
	{
		x = 1;
		y = 0;
	}
	else
	{
		extend_gcd(b,a % b,y,x);
		y -= x * (a / b);
	}
}

void solve()
{
	long long tmp1 = n + 1 - m;
	for(long long i = 1;i <= n + m;i++)
	{
		tmp1 = (long long)(tmp1 * i) % mo;;
	}
	long long tmp2 = n + 1;
	for(int i = 1;i <= n;i++)
	{
		tmp2 = (long long)(tmp2 * i) % mo;
	}
	for(int i = 1;i <= m;i++)
	{
		tmp2 = (long long)(tmp2 * i) % mo;
	}
	extend_gcd(tmp2,mo,x,y);
	ans = (tmp1 * (x + mo) % mo) % mo;
	printf("%lld\n",ans);
}

int main()
{
	init();
	readdata();
	solve();
	return 0;
}