【JZOJ A组】登山

Description

恶梦是一个登山爱好者，今天他来到了黄山。
俗话说的好，不走回头路。所以在黄山，你只能往前走，或者往上走。并且很显然的是，当你走到山脊的时候，你不能够往上走，你只能往前走一步再往上走。
抽象一点而言就是，你可以把黄山视为一个N * N格点图，恶梦从(0,0)开始出发，要走到(N,N)。当他走到位置(x,y)的时候，它可以往(x + 1,y),或(x,y+1)走。
并且当他走到(x,x)的时候，由于他已经处在了山脊上，所以他不能够往(x,x+1)方向上走。
当恶梦兴致勃勃准备开始爬山的时候，他的同伴告诉他，黄山由于年久失修，有一些位置出现了大坑，不能走。恶梦觉得更刺激了，但他想先知道他能有多少种方式走到黄山顶。
由于这个数字很大，所以你只需要将答案对10^9 + 7取模输出即可。

Input

第一行包括两个整数N,C,分别表示你可以把黄山视作一个N * N的格点图，并且黄山上面有C个位置出现了大坑。
接下来的C行，每行包括两个整数X,Y,表示X,Y这个位置不能走。保证X>=Y,也就是说(X,Y)必然在山上。
保证这C个点互不相同。

Output

输出只有一个整数Ans,表示恶梦爬上山顶的路径数对10^9+7取模的值。

Sample Input

输入1：
5 2
5 0
1 1
输入2：
7 4
6 5
5 3
2 1
7 1

Sample Output

输出1：
27
输出2：
34

Data Constraint

对于30%的数据,保证N<=5000
对于另外20%的数据，保证C=0
对于另外20%的数据，保证C=1
对于100%的数据，保证N<=100000,C<=1000
保证对于(0,0),(N,N)不存在障碍点。

思路

很像卡特兰数的思想

用总的减去不合法的

对于每一个坑，可以做一个关于直线y=x+1的对称点，那么从起点到那个点的距离也是O（1）可求的

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#define f first
#define S second
#define ll long long
using namespace std;
typedef pair<int,int> P;
const int maxm=1005,maxn=200015,mod=int(1e9)+7;
P b[maxm];
int fc[maxn],Rv[maxn],f[maxn],n,C;
int power(int x,int c)
{
	ll b=1;
	while(c)
	{
		if(c&1) b=b*x*1ll%mod;
		c>>=1; x=(long long)x*x*1ll%mod;
	}
	return b;
}
int getc(int n,int m)
{
	if (m>n) return 0;
	return fc[n]*1ll*Rv[m]%mod*Rv[n-m]%mod;
}
int normal(int s,int t,int x,int y)
{
	if (x<s||y<t) return 0;
	return getc(x+y-s-t,x-s);
}
int walk(int s,int t,int x,int y)
{
	return (normal(s,t,x,y)-normal(s,t,y-1,x+1)+mod)%mod;
}
int main()
{
	scanf("%d%d",&n,&C);
	fc[0]=1;
	for(int i=1; i<=n*2+5; i++) fc[i]=fc[i-1]*1ll*i%mod;
	Rv[n*2+5]=power(fc[n*2+5],mod-2);
	for(int i=n*2+5; i; i--) Rv[i-1]=Rv[i]*1ll*i%mod;
	for(int i=1; i<=C; i++) scanf("%d%d",&b[i].f,&b[i].S);
	b[++C]=P(n,n);
	sort(b+1,b+C+1);
	for(int i=1; i<=C; i++)
	{
		f[i]=walk(0,0,b[i].f,b[i].S);
		for(int j=1; j<i; j++)
		if (b[j].S<=b[i].S)
			f[i]=(f[i]-f[j]*1ll*walk(b[j].f,b[j].S,b[i].f,b[i].S)%mod+mod)%mod;
	}
	printf("%d\n",f[C]);
}