Count Path Pair

题目地址：点击打开链接

【分析】

我们设总的方案数为A，然后f，g相交的方案数为B。那么答案就是A-B了。

A很好计算，A=C(m+n,n)*C(m-p+q,q)

从图中我们可以看出，f,g在E点相交，等价于f从E点穿过。所以相交的方案数B就等价于“A到C，B到D的方案数”。

那么，B=C(m+q,q)*C(m-p+n,n)

至此，问题完美解决

【代码】

/****************************
    ID:Ciocio
    LANG:C++
    DATE:2014-1-26
    TASK:Count Path Pair
****************************/
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <iostream>
#include <algorithm>
#include <queue>

using namespace std;

#define MODER 100000007
#define MAXN 200000

typedef long long LL;
LL m,n,p,q,Plus,Sub;
LL f[MAXN+10];

LL _inv(LL a){
	int b=MODER-2;
	LL rt=1;
	a%=MODER;
	while(b){
		if(b&1) rt=(rt*a)%MODER;
		b>>=1;
		a=(a*a)%MODER;
	}
	return rt;
}

LL C(LL a,LL b){
	LL rt=f[a];
	rt=(rt*(_inv(f[b]*f[a-b]%MODER)))%MODER;
	return rt;
}

void _init(){
	Plus=C(m+n,n)*C(m-p+q,q)%MODER;
	Sub=C(m+q,q)*C(m-p+n,n)%MODER;
}

void _solve(){
	cout<<(Plus-Sub+MODER)%MODER<<endl;
}

int main(){
	f[1]=1;
	for(int i=2;i<=MAXN;i++)
	    f[i]=(f[i-1]*i)%MODER;
	while(cin>>m>>n>>p>>q){
	    _init();
	    _solve();
	}
	return 0;
}