题目地址:点击打开链接
【分析】
我们设总的方案数为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;
}