51nod 连环进阶题 1118-1119-1120

51nod 1118 1119 1120 机器人走方格问题

Question1(for 1118 1119):

M * N的方格，一个机器人从左上走到右下，只能向右或向下走。有多少种不同的走法？由于方法数量可能很大，只需要输出Mod 10^9 + 7的结果。

Question2(for 1120):

N * N的方格，从左上到右下画一条线。一个机器人从左上走到右下，只能向右或向下走。并要求只能在这条线的上面或下面走，不能穿越这条线，有多少种不同的走法？由于方法数量可能很大，只需要输出Mod 10007的结果。

 

针对第一个问题：

很好想到甚至小学的时候或许就接触推导过，就是最经典的dp O(∩_∩)O 应该对大家没压力叭~

左上角为f[1][1],右下角为f[m][n],f[i][1]=1,f[1][i]=1;f[i][j]=f[i-1][j]+f[i][j-1]，这样就应该AC了1118

但蒟蒻表示没看见1118就先做了1119，后来就直接帖啦~qwq

那不急这就讲1119，进阶版辣，数据范围变到了n,m<=1000000

这就完全不能用这么传统的dp辣，显然要用到数论方面的知识喽~

其实也不难，如果我们知道向右走是怎样走的，那我们自然就知道了向下怎么走，对吗？

所以只要知道向右走的走法就是最后的答案。辣么当然答案就是C(min(n-1,m-1),m+n-2)喽~

这应该很显然吧。

 

 

 

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>
#define ll long long
using namespace std;
const int mod=1e9+7;
ll exgcd(ll a,ll b,ll &x,ll &y){
	if(!a && !b) return -1;
	if(!b){x=1;y=0;return a;}
	ll d=exgcd(b,a%b,y,x);
	y-=a/b*x;
	return d;
}
ll reverse(ll n,ll mo){
	ll x,y;
	ll d=exgcd(n,mo,x,y);
	if(d==1) return(x%mo+mo)%mo;
	else return -1;
}
ll C(ll m,ll n){
	ll s1=1,s2=1;
	for(int i=n;i>=n-m+1;i--) s1=s1*i%mod;
	for(int i=1;i<=m;i++) s2=s2*i%mod;
	return s1*reverse(s2,mod)%mod;
}
int main()
{
	ll n,m,ans;
	scanf("%lld%lld",&n,&m);
	ans=C(min(n-1,m-1),m+n-2);
	printf("%lld\n",ans);
	return 0;
}

 

 

 

再看第二个问题:

这应该略有一些难度，会用到lucas定理，最近刚做到过那就用到啦。

而且表示崩溃的是，这题不仅要用lucas定理，还是卡特兰数。

那又要用到卡特兰数，鉴于数据范围这么大，还得预处理，蒟蒻崩溃啊啊啊！

还是坚持下来了(虽然本蒟蒻什么都不知道)，因为只能一直在上面或一直在下面走，所以答案就是卡特兰数*2

卡特兰数又是多少呢，就是C(n+n,n)-C(n+n,n-1)注意有模所以在减之前先加上模数。

至于lucas定理，比较高深，我会再做讲解qwq

 

 

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>
#define mod 10007
using namespace std;
int product[mod],inv[mod];
inline int C(int n,int m){
	if(n<m) return 0;
	return product[n]*inv[m]%mod*inv[n-m]%mod;
}
int lucas(int n,int m){
	return m==0 ? 1 : C(n%mod,m%mod)*lucas(n/mod,m/mod)%mod;
}
int main()
{
	int n;
	scanf("%d",&n);
	n--;
	product[0]=1;
	for(int i=1;i<mod;i++) product[i]=product[i-1]*i%mod;
	inv[1]=1;
	for(int i=2;i<mod;i++) inv[i]=(mod-mod/i)*inv[mod%i]%mod;
	inv[0]=1;
	for(int i=1;i<mod;i++) inv[i]=inv[i-1]*inv[i]%mod;
	printf("%d\n",(lucas(n+n,n)+mod-lucas(n+n,n-1))%mod*2%mod);
	return 0;
}