[CF578F]Mirror Box

Description

有一个n*m的网格,每个格子是一个斜45°摆放的镜子
现在有一些格子上的镜子不见了,你需要求出有多少种摆放镜子的方法,满足:
1:从网格某一个边缘的中点射进来的光线会从这条边的一个临边射出
2:网格图中每一条网格线都存在1中的一条光线,满足这条光线经过这个网格线
在这里插入图片描述n,m<=100,保证空格子的数量<=200

Solution

考虑把所有格点黑白染色,我们可以发现,一种合法的方案对应黑/白格子上的一棵生成树
在这里插入图片描述
简单证明的话,第二个条件相当于图中无环,第一个条件相当于我们可以把边缘分组,然后每一组用一条封闭的线围住
空格子的数量<=200的话,我们把所有边缩在一起,剩余的点数<=401
直接矩阵树算生成数个数即可

Code

#include <cstdio>
#include <cstring>
#include <algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;

typedef long long ll;

const int N=105,M=505;

int n,m,Mo,fa[N*N],id[N*N],tot[2];
char a[N][N];

int Id(int x,int y) {return x*(m+1)+y;}
int get(int x) {return fa[x]==x?x:fa[x]=get(fa[x]);}

void merge(int x,int y) {
	x=get(x);y=get(y);
	fa[y]=x;
}

int pwr(int x,int y) {
	int z=1;
	for(;y;y>>=1,x=(ll)x*x%Mo)
		if (y&1) z=(ll)z*x%Mo;
	return z;
}

struct Graph{
	int a[M][M];
	void add(int x,int y) {
		a[x][x]++;a[y][y]++;
		a[x][y]--;a[y][x]--;
	}
	int DET(int n) {
		int ret=1;
		fo(i,1,n) {
			if (!a[i][i])
				fo(j,i+1,n)
					if (a[j][i]) {
						fo(k,i,n) swap(a[i][k],a[j][k]);
						ret=-ret;break;
					}
			if (!a[i][i]) break;
			int inv=pwr(a[i][i],Mo-2);
			fo(j,i+1,n) {
				int w=(ll)a[j][i]*inv%Mo;
				fo(k,i,n) (a[j][k]-=(ll)w*a[i][k]%Mo)%=Mo;
			}
		}
		fo(i,1,n) ret=(ll)ret*a[i][i]%Mo;
		return ret;
	}
}G[2];

int main() {
	scanf("%d%d%d",&n,&m,&Mo);
	fo(i,1,n) scanf("%s",a[i]+1);
	fo(i,0,(n+1)*(m+1)) fa[i]=i;
	fo(i,1,n) fo(j,1,m) {
		if (a[i][j]=='/') merge(Id(i,j-1),Id(i-1,j));
		if (a[i][j]=='\\') merge(Id(i-1,j-1),Id(i,j));
	}
	fo(i,0,n) fo(j,0,m) {
		int x=Id(i,j);
		if (fa[x]==x) id[x]=++tot[(i+j)&1];
	}
	fo(i,1,n)
		fo(j,1,m)
			if (a[i][j]=='*') {
				G[(i+j)&1].add(id[get(Id(i-1,j-1))],id[get(Id(i,j))]);
				G[((i+j)&1)^1].add(id[get(Id(i,j-1))],id[get(Id(i-1,j))]);
			}
	int ans=(G[0].DET(tot[0]-1)+G[1].DET(tot[1]-1))%Mo;
	printf("%d\n",(ans+Mo)%Mo);
	return 0;
}
©️2020 CSDN 皮肤主题: 编程工作室 设计师: CSDN官方博客 返回首页
实付0元
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、C币套餐、付费专栏及课程。

余额充值