hdoj 5852 Intersection is not allowed! 【组合数学】

本文链接：https://blog.csdn.net/chenzhenyu123456/article/details/52235185

本文介绍了一种解决不相交路径计数问题的方法，利用行列式的计算来找到从多个起点到对应终点的所有可能路径数量。通过编程实现，具体展示了如何使用C++代码来解决这个问题。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

hdoj 5852 Intersection is not allowed!
题意： $n*n$ 的方格上， $k$ 个起点、 $k$ 个终点，第 $i$ 个起点对应第 $i$ 个终点，每次只可以向下或者向右移动一步，问你不相交的路径有多少条？

看过朝鲜爷题解后，ORZ。每个起点对于所有终点，求出方案数。列出一个行列式，求解行列式的值。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <string>
#include <cstring>
#include <map>
#include <vector>
#include <queue>
#include <stack>
#include <set>
#include <cstdlib>
#define ll o<<1
#define rr o<<1|1
#define CLR(a, b) memset(a, (b), sizeof(a))
using namespace std;
typedef long long LL;
typedef pair<LL, LL> pii;
const int MAXN = 2e5 + 10;
const int INF = 1e7 + 10;
const int MOD = 1e9 + 7;
void add(LL &x, LL y) { x += y; x %= MOD; }
LL a[110][110];
int X1[110], Y1[110], X2[110], Y2[110];
LL fac[MAXN], idv[MAXN];
LL pow_mod(LL a, int n) {
    LL ans = 1;
    while(n) {
        if(n & 1) {
            ans = ans * a % MOD;
        }
        a = a * a % MOD;
        n >>= 1;
    }
    return ans;
}
LL C(int n, int m) {
    return fac[n] * idv[n - m] % MOD * idv[m] % MOD;
}
LL Solve(int n) {
    LL ans = 1; int T = 0;
    for(int i = 1; i <= n; i++) {
        for(int j = i + 1; j <= n; j++) {
            //a[j][i] a[i][i];
            int x = i, y = j;
            for(; a[y][i] != 0 ;) {
                LL ta = a[x][i] / a[y][i];
                for(int k = i; k <= n; k++) {
                    a[x][k] = ((a[x][k] - ta * a[y][k]) % MOD + MOD) % MOD;
                }
                swap(x, y);
            }
            if(x != i) {
                for(int k = i; k <= n; k++) {
                    swap(a[i][k], a[x][k]);
                }
                T ^= 1;
            }
        }
        ans = ans * a[i][i] % MOD;
        if(ans == 0) return 0;
    }
    if(T & 1) {
        ans = (-ans + MOD) % MOD;
    }
    return ans;
}
int main()
{
    fac[0] = 1;
    for(int i = 1; i <= MAXN - 1; i++) {
        fac[i] = fac[i - 1] * i % MOD;
    }
    for(int i = 0; i <= MAXN - 1; i++) {
        idv[i] = pow_mod(fac[i], MOD - 2);
    }
    int t; scanf("%d", &t);
    while(t--) {
        int n, k; scanf("%d%d", &n, &k);
        for(int i = 1; i <= k; i++) {
            scanf("%d", &Y1[i]);
            X1[i] = 1;
        }
        for(int i = 1; i <= k; i++) {
            scanf("%d", &Y2[i]);
            X2[i] = n;
        }
        for(int i = 1; i <= k; i++) {
            for(int j = 1; j <= k; j++) {
                if(Y1[i] <= Y2[j]) {
                    a[i][j] = C(X2[j] - X1[i] + Y2[j] - Y1[i], X2[j] - X1[i]);
                }
                else {
                    a[i][j] = 0;
                }
            }
        }
        printf("%lld\n", Solve(k));
    }
    return 0;
}