luogu P4111 [HEOI2015]小Z的房间

背景：

今天马刺被掘金$19pts$逆转。

题目传送门：

https://www.luogu.org/problemnew/show/P4111

题意：

给定一幅地图（有空地和墙），现在你要走遍所有的空地，求方案数。

思路：

基尔霍夫矩阵&矩阵树定理
搞出$A,D$矩阵，求出$K$矩阵后高斯消元即可。
这道题有点奇怪，可以不打逆元，因为模数不是质数，因此可能没有逆元。

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define mod (int)(1e9)
#define LL long long
using namespace std;
	int D[110][110],A[110][110],K[110][110],id[20][20];
	int n,m,tot=0;
void calc(int x,int y)
{
	D[x][x]++,D[y][y]++;A[x][y]++,A[y][x]++;
}
int Matrix_tree(int tot)
{
	int ans=1;
	for(int i=1;i<=tot;i++)
	{
		for(int j=i+1;j<=tot;j++)
			while(K[j][i])
			{
				int t=K[i][i]/K[j][i];
				for(int k=i;k<=tot;k++) K[i][k]=(K[i][k]-(LL)K[j][k]*t%mod+mod)%mod;
				for(int k=i;k<=tot;k++) swap(K[i][k],K[j][k]);
				ans=-ans;
			}
		ans=((LL)ans*K[i][i]%mod+mod)%mod;
	}
	return ans;
}
int main()
{
	char s[20];
	scanf("%d %d",&n,&m);
	for(int i=1;i<=n;i++)
	{
		scanf("%s",s+1);
		for(int j=1;j<=m;j++)
			if(s[j]=='.') id[i][j]=++tot;
	}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			if(id[i][j])
			{
				if(id[i-1][j]) calc(id[i][j],id[i-1][j]);
				if(id[i][j-1]) calc(id[i][j],id[i][j-1]);
			}
	for(int i=1;i<=tot;i++)
		for(int j=1;j<=tot;j++)
			K[i][j]=D[i][j]-A[i][j];
	printf("%d",Matrix_tree(tot-1));
}