题意: 4个点A(0,0),B(p,0),C(m,q),D(m,n),保证m>p,n>q,求从A走到D和从B走到C两条路径不相交的走法的种数,只能up和right。
解法: 直接求不相交路径不好求,那么就求所有的走法,再减去所有的相交路径。而相交的路径种数,其实等价于从A走到C和从B走到D的走法种数,因为从A到C和从B到D路径必相交,并且A和B只要在交点交换一下路径,变成A去D,B去C,就是A~D和B~C相交情况下的路径了。
answer = C(n+m,n) * C(m-p+q,q) - C(m+q,q) * C(m+n-p,n)
= ( (n+m)! * (m-p+q)! - (m+q)! * (m+n-q)! ) / ( n! * m! * q! * (m-p)! ) % MOD
除法部分求一下逆元即可
关于逆元:
X / a = ? %MOD
--> x * (a^-1) = ? % MOD
求逆元:
a * (a^-1) = 1 % MOD
--> a*(a^-1) + k * MOD = 1
exgcd求解即可。
/* Created Time: Tuesday, November 05, 2013 PM11:33:46 CST */
#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
const int M = 201000;
const int MOD = 1e8+7;
typedef long long lld;
int n,m,p,q,jie[M];
int exgcd(int a,int b,int &x,int &y) {
if (!b) {
x = 1,y = 0;
return a;
} else {
int c = exgcd(b,a%b,y,x);
y -= a/b*x;
return c;
}
}
void init() {
jie[0] = 1;
for (int i = 1; i < M; i ++)
jie[i] = (lld)jie[i-1]*i%MOD;
}
void work() {
int ans = 0;
ans = (lld)jie[n+m]*jie[m-p+q]%MOD;
ans = ((lld)ans-(lld)jie[m-p+n]*jie[m+q]%MOD)%MOD;
ans = (ans+MOD)%MOD;
int inv,d;
int a = (lld)jie[n]*jie[m]%MOD*jie[q]%MOD*jie[m-p]%MOD;
exgcd(a,MOD,inv,d);
inv = (inv+MOD)%MOD;
ans = (lld)ans*inv%MOD;
printf("%d\n",ans);
}
int main() {
init();
while (~scanf("%d%d%d%d",&m,&n,&p,&q))
work();
return 0;
}