【BZOJ4031】【HEOI2015】小Z的房间（Matrix-Tree定理，线性同余迭代消元）

Description

你突然有了一个大房子，房子里面有一些房间。事实上，你的房子可以看做是一个包含n*m个格子的格状矩形，每个格子是一个房间或者是一个柱子。在一开始的时候，相邻的格子之间都有墙隔着。

你想要打通一些相邻房间的墙，使得所有房间能够互相到达。在此过程中，你不能把房子给打穿，或者打通柱子（以及柱子旁边的墙）。同时，你不希望在房子中有小偷的时候会很难抓，所以你希望任意两个房间之间都只有一条通路。现在，你希望统计一共有多少种可行的方案。

Solution

Matrix-Tree定理模板题。
定理本身没什么好说的，记住就行反正也不会证明
这题的模数比较恶心:$10^9$
开始打了个乘法逆元处理除法，发现$10^9$根本TM没有逆元。。。
于是就不会做了，orz了一发zhou888，发现他写了个类似辗转相除法的东西消元，感觉非常巧妙，具体可以看看代码。

Code

/************************************************
 * Au: Hany01
 * Date: Apr 4th, 2018
 * Prob: [BZOJ4031][HEOI2015] 小z的房间
 * Email: hany01@foxmail.com
************************************************/

#include<bits/stdc++.h>

using namespace std;

typedef long long LL;
typedef pair<int, int> PII;
#define File(a) freopen(a".in", "r", stdin), freopen(a".out", "w", stdout)
#define rep(i, j) for (register int i = 0, i##_end_ = (j); i < i##_end_; ++ i)
#define For(i, j, k) for (register int i = (j), i##_end_ = (k); i <= i##_end_; ++ i)
#define Fordown(i, j, k) for (register int i = (j), i##_end_ = (k); i >= i##_end_; -- i)
#define Set(a, b) memset(a, b, sizeof(a))
#define Cpy(a, b) memcpy(a, b, sizeof(a))
#define x first
#define y second
#define pb(a) push_back(a)
#define mp(a, b) make_pair(a, b)
#define ALL(a) (a).begin(), (a).end()
#define SZ(a) ((int)(a).size())
#define INF (0x3f3f3f3f)
#define INF1 (2139062143)
#define Mod (1000000000)
#define debug(...) fprintf(stderr, __VA_ARGS__)
#define y1 wozenmezhemecaia

template <typename T> inline bool chkmax(T &a, T b) { return a < b ? a = b, 1 : 0; }
template <typename T> inline bool chkmin(T &a, T b) { return b < a ? a = b, 1 : 0; }

inline int read()
{
    register int _, __; register char c_;
    for (_ = 0, __ = 1, c_ = getchar(); c_ < '0' || c_ > '9'; c_ = getchar()) if (c_ == '-') __ = -1;
    for ( ; c_ >= '0' && c_ <= '9'; c_ = getchar()) _ = (_ << 1) + (_ << 3) + (c_ ^ 48);
    return _ * __;
}

#define idx(i, j) ((i - 1) * m + j)

inline LL Pow(LL a, LL b) {
    LL Ans = 1;
    for ( ; b; b >>= 1, (a *= a) %= Mod) if (b & 1) (Ans *= a) %= Mod;
    return Ans;
}

struct Line { LL a[120]; };
struct Matrix { Line a[120]; }A;
int id[10][10], n, m, tot;
char s[10][10];

Line operator + (Line A, Line B) {
    Line C;
    For(i, 1, tot) C.a[i] = (A.a[i] + B.a[i]) % Mod;
    return C;
}

Line operator - (Line A, Line B) {
    Line C;
    For(i, 1, tot) C.a[i] = (A.a[i] - B.a[i]) % Mod;
    return C;
}

Line operator * (Line A, LL d) {
    Line C;
    For(i, 1, tot) C.a[i] = A.a[i] * d % Mod;
    return C;
}

inline void insert(int u, int v) {
    A.a[u].a[v] = Mod - 1, A.a[v].a[u] = Mod - 1, ++ A.a[u].a[u], ++ A.a[v].a[v];
}

inline LL calc(Matrix A, int n)
{
    LL Ans = 1;
    For(i, 1, n)
    {
        if (!A.a[i].a[i]) {
            register int mark = 0;
            For(j, i + 1, n) if (A.a[j].a[i]) {
                swap(A.a[i], A.a[j]), mark = 1;
                break;
            }
            //assert(mark);
        }
        For(j, i + 1, n) if (A.a[j].a[i]) {
            int x = A.a[i].a[i], y = A.a[j].a[i];
            while (y)
            {
                int t = x / y;
                x %= y, swap(x, y);
                A.a[i] = A.a[i] - A.a[j] * t;
                swap(A.a[i], A.a[j]), Ans = -Ans;
            }
            if (!A.a[i].a[i]) swap(A.a[i], A.a[j]), Ans = -Ans;
        }
        //assert(A.a[i].a[i]);
        (Ans *= A.a[i].a[i]) %= Mod;
    }
    return (Ans % Mod + Mod) % Mod;
}

int main()
{
#ifdef hany01
    File("bzoj4031");
#endif

    //Input
    n = read(), m = read();
    For(i, 1, n) scanf("%s", s[i] + 1);
    For(i, 1, n) For(j, 1, m) if (s[i][j] == '.') id[i][j] = ++ tot;

    //Build the graph
    For(i, 1, n) For(j, 1, m)
        if (s[i][j] == '.') {
            if (j != m && s[i][j + 1] == '.') insert(id[i][j], id[i][j + 1]);
            if (i != n && s[i + 1][j] == '.') insert(id[i][j], id[i + 1][j]);
        }

    //Calc
    printf("%lld\n", calc(A, tot - 1));

    return 0;
}
//天长地久有时尽，此恨绵绵无绝期。
//    -- 白居易《长恨歌》