BZOJ4031: [HEOI2015]小Z的房间【矩阵树定理】

4031: [HEOI2015]小Z的房间

裸的矩阵树定理。

#include<cstdio>
#include<algorithm>
using namespace std;
const int MOD=1e9;
int n,m,cnt,Ans,mp[105][105],a[105][105];char ch[15];
void Add(int x,int y){a[x][x]++,a[y][y]++,a[x][y]--,a[y][x]--;}
int Gauss(){
	Ans=1;
	for(int i=1;i<cnt;i++){
		for(int j=i+1;j<cnt;j++)
		while(a[j][i]){
			int t=a[i][i]/a[j][i];
			for(int k=i;k<cnt;k++) a[i][k]=(a[i][k]-1ll*t*a[j][k]%MOD+MOD)%MOD;
			swap(a[i],a[j]);Ans=-Ans;
		}
		Ans=1ll*Ans*a[i][i]%MOD;
	}
	return (Ans+MOD)%MOD;
}
int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++){
		scanf("%s",ch+1);
		for(int j=1;j<=m;j++) if(ch[j]=='.') mp[i][j]=++cnt;
	}
	for(int i=1;i<=n;i++)
	for(int j=1;j<=m;j++)
	if(mp[i][j]){
		if(mp[i][j-1]&&j>1) Add(mp[i][j],mp[i][j-1]);
		if(mp[i-1][j]&&i>1) Add(mp[i][j],mp[i-1][j]);
	}
	printf("%d\n",Gauss());
	return 0;
}