codeforces560EGerald and Giant Chess（组合数学）

题目链接：

http://codeforces.com/contest/560/problem/E

题目大意：

从（1,1）走到(h,w)，要求不能经过题目中所给的坏点，问最后有几种走法。

思路：

如果从(1,1)走到（x，y），可以知道有C(x-1+y-1,x-1)种。

然后考虑将走到坏点的情况减掉。

定义一个f[i]为从出发点到坏点i且不经过其他坏点的路径走法。对于某个坏点i来说，如果他的左上角有一个坏点j，那么我必然要将他减掉，那么到i的走法就是C（dx+dy,dx）-f[j]*从j到i的走法。

所以我只要先将坏点按照x,y从小到大排序，然后枚举坏点i，将其前面的坏点全部减去就好了。而且根据我们对f[i]的定义，全部减去的过程中是不会有重复的。

代码：

#include<stdio.h>
#include<string.h>
#include<math.h>
#include<algorithm>
#define ll __int64
#define mod 1000000007
using namespace std;
ll f[2005];
ll fac[200005],inv[200005];
struct node{
	int x,y;
}p[2005];
ll quick(ll a,ll b){
	ll ans=1;
	while(b){
		if(b%2)ans=ans*a%mod;
		a=a*a%mod;
		b=b/2;
	}
	return ans%mod;
}
void extgcd(ll a,ll b,ll& d,ll& x,ll& y){  
    if(!b){  
        d=a;x=1;y=0;  
    }else{  
        extgcd(b,a%b,d,y,x);  
        y-=x*(a/b);  
    }  
}  
ll mod_inv(ll a,ll n){//a在模n意义下的逆元  
    ll d,x,y;  
    extgcd(a,n,d,x,y);  
    return d==1?(x+n)%n:-1;  
}  
  
bool cmp(node a,node b){
	if(a.x!=b.x)return a.x<b.x;
	return a.y<b.y;
}
void init(){
	fac[0]=fac[1]=inv[0]=inv[1]=1;
	for(int i=2;i<=200000;i++){
		fac[i]=fac[i-1]*i%mod;
	//	inv[i]=quick(fac[i],mod-2)%mod;
		inv[i]=mod_inv(fac[i],mod)%mod;
		}
}
ll CC(ll a,ll b){
	ll ans;
	ans=fac[a]%mod;
	ans=(ans*inv[b])%mod;
	ans=(ans*inv[a-b])%mod;
	return ans%mod;
}
int main()
{
	ll h,w,n;
	int i,j,k;
	init();
	scanf("%I64d%I64d%I64d",&h,&w,&n);
	for(i=1;i<=n;i++){
		scanf("%d%d",&p[i].x,&p[i].y);
	}
	sort(p+1,p+1+n,cmp);
	memset(f,0,sizeof(f));
	p[n+1].x=h;
	p[n+1].y=w;
	ll as=CC(h+w-2,h-1);;
	for(i=1;i<=n+1;i++){
		int dx,dy;
		dx=p[i].x-1;
		dy=p[i].y-1;
		ll tmp=0;
		tmp=CC(dy+dx,dx);
		for(j=1;j<i;j++){
			if(p[i].x>=p[j].x&&p[i].y>=p[j].y){
			int xx=p[i].x-p[j].x;
				int yy=p[i].y-p[j].y;
				tmp-=f[j]*CC(xx+yy,xx);
				tmp=(tmp+mod)%mod;
			}
		}
		f[i]=tmp;
	}
	printf("%I64d\n",(f[n+1]+mod)%mod);
}
/*
100000 100000 4
50001 50001
50000 50000
50000 50001
50001 50000
*/