洛谷传送门
BZOJ传送门
题目描述
你突然有了一个大房子,房子里面有一些房间。事实上,你的房子可以看做是一个包含 n∗m n ∗ m 个格子的格状矩形,每个格子是一个房间或者是一个柱子。在一开始的时候,相邻的格子之间都有墙隔着。
你想要打通一些相邻房间的墙,使得所有房间能够互相到达。在此过程中,你不能把房子给打穿,或者打通柱子(以及柱子旁边的墙)。同时,你不希望在房子中有小偷的时候会很难抓,所以你希望任意两个房间之间都只有一条通路。现在,你希望统计一共有多少种可行的方案。
输入输出格式
输入格式:
第一行两个数分别表示 n n 和。
接下来 n n 行,每行个字符,每个字符都会是’.’或者’ ∗ ∗ ’,其中’.’代表房间,’*’代表柱子。
输出格式:
一行一个整数,表示合法的方案数
输入输出样例
输入样例#1:
2 2
..
..
输出样例#1:
4
输入样例#2:
2 2
*.
.*
输出样例#2:
0
说明
对于前20%的数据, n,m≤3 n , m ≤ 3
对于前50%的数据, n,m≤5 n , m ≤ 5
对于前100%的数据, n,m≤9 n , m ≤ 9
有40%的数据保证, min(n,m)≤3 m i n ( n , m ) ≤ 3
有30%的数据保证,不存在柱子
解题分析
矩阵树定理裸题…
做这道题之前, 我们要了解关于基尔霍夫矩阵和行列式的一些定理:
基尔霍夫矩阵的行列式就是对应图的生成树数量。
这意味着我们可以把生成树计数问题转换为矩阵上的问题。
交换矩阵中两行, 行列式为相反数。
这意味着我们可以任意交换矩阵中的两行,而仍可以轻松计算出行列式的值。
矩阵中一行所有元素减去另一行对应元素的相同倍数, 行列式不变。
这意味着我们可以像高斯消元一样消去某一列中的元素。
一个上三角矩阵的行列式等于其对角线上的元素之积
有了这一条和之前的性质, 我们就可以将矩阵消成上三角矩阵, 并计算出其行列式来了。
这道题因为需要模一个合数, 不能直接暴力消元, 我们就可以利用以上2、 3条性质, 用辗转相除的方法消元。
代码如下:
#include <cstdio>
#include <cstring>
#include <cctype>
#include <algorithm>
#include <cmath>
#include <cstdlib>
#define R register
#define IN inline
#define gc getchar()
#define W while
#define MX 205
#define ll long long
#define MOD 1000000000ll
bool mp[MX][MX];
int id[MX][MX];
ll mat[MX][MX];
int n, m, cnt;
IN ll Gauss()
{
ll ans = 1;
ll mul;
for (R int i = 1; i < cnt; ++i)
{
for (R int j = i + 1; j < cnt; ++j)
{
W (mat[j][i])
{
mul = mat[i][i] / mat[j][i];
for (R int k = i; k < cnt; ++k)
{
mat[i][k] = ((mat[i][k] - mul * mat[j][k]) % MOD + MOD) % MOD;
std::swap(mat[i][k], mat[j][k]);
}
ans = -ans;
}
}
if(!mat[i][i]) return 0;
ans = (ans * mat[i][i] + MOD) % MOD;
}
return (ans % MOD + MOD) % MOD;
}
int main(void)
{
char buf[500];
scanf("%d%d", &n, &m);
for (R int i = 1; i <= n; ++i)
{
scanf("%s", buf + 1);
for (R int j = 1; j <= m; ++j)
{
mp[i][j] = (buf[j] == '.');
}
}
for (R int i = 1; i <= n; ++i)
{
for (R int j = 1; j <= m; ++j)
{
if(mp[i][j])
{
id[i][j] = ++cnt;
if(mp[i - 1][j])
{
++mat[id[i - 1][j]][id[i - 1][j]];
++mat[cnt][cnt];
mat[id[i - 1][j]][cnt] = -1;
mat[cnt][id[i - 1][j]] = -1;
}
if(mp[i][j - 1])
{
++mat[id[i][j - 1]][id[i][j - 1]];
++mat[cnt][cnt];
mat[id[i][j - 1]][cnt] = -1;
mat[cnt][id[i][j - 1]] = -1;
}
}
}
}
printf("%lld", Gauss());
}