P4111 [HEOI2015]小 Z 的房间 图的生成树个数 高斯消元计算行列式

题意

求一个图的生成树个数

分析

拉普拉斯生成树？？？

生成树矩阵=度数矩阵-邻接矩阵
生成树的个数等于：生成树矩阵降一维后求行列式的值

采用高斯消元的方式计算当前行列式的值

注意：本题中每一个矩形都算图中的一个点，计算行列式的矩阵需要开$n^2$

代码

#include <iostream>
#include <string>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <queue>
#include <iomanip>
#include <map>
#include <cstdio>
#include <stack>
#include <set>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int ,int > pii;
#define endl '\n'
ll gcd(ll a, ll b){
    return b == 0 ? a : gcd(b, a % b);
}
void input(){
    freopen("in.txt", "r", stdin);
    freopen("out.txt", "w", stdout);
}
inline int read(){
	int x=0,f=1;char c=getchar();
	while(c<'0'||c>'9') {if(c=='-') f=-1;c=getchar();}
	while (c>='0'&&c<='9') x=(x<<3)+(x<<1)+(c^48),c=getchar();
	return x*f;
}
const int N = 110;
const int mod = 1e9;
int A[N][N], id[N][N], n, m;
char mp[N][N];
void add(int id1, int id2){
    A[id1][id1]++, A[id2][id2]++, A[id1][id2]--, A[id2][id1]--;
}
// 求一个图的生成树个数，（度数矩阵-邻接矩阵）随机降一维的行列式值
// 求行列式的值
int gauss(int n){
    ll ans = 1;
    for(int i = 1; i < n; i++)
		for(int j = i + 1; j <= n; j++){
			while(A[j][i]){
				int l = A[i][i] / A[j][i];
				for(int k = 1; k <= n; k++) 
					A[i][k] = (A[i][k] - 1ll*A[j][k] * l % mod + mod) % mod;
				for(int k = 1; k <= n; k++) swap(A[i][k],A[j][k]);
				ans *= -1;
			}
		}
	for(int i = 1; i <= n; i++) ans=(ans * A[i][i] % mod + mod) %mod;
    return ans;
}
int main(){
    ios::sync_with_stdio(0); cin.tie(0); cout.tie(0);
    cin>>n>>m;
    int cnt = 0;
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= m; j++){
            cin>>mp[i][j];
            if(mp[i][j] == '.') id[i][j] = ++cnt;
        }
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= m; j++){
            if(mp[i][j] == '.' && mp[i+1][j] == '.') add(id[i][j], id[i+1][j]);
            if(mp[i][j] == '.' && mp[i][j+1] == '.') add(id[i][j], id[i][j+1]);
        }
    cnt--;
    
    cout<<gauss(cnt)<<endl;
    return 0;
}