ZOJ 3624 Count Path Pair (组合计数，乘法逆元)

乘法逆元

题意： 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;
}