标签：Box return int CF578F fa Mirror inline find

description

CF578F

solution

• 考虑转化题目的要求
  1.对于任意一条边，都存在一条从界垂直射入的光线，经过反射穿过这条边。
  当图中有环时，环内的边一定不满足条件，而在不存在环时感性理解一下就能满足条件
  2.从任意一条界垂直射入的光线经过反射，从相邻的一条界射出；
  对于这个条件，因为光路可逆，所以我们需要将图中的边两两分组，并用镜子将他们围住
• 多画几张图后，你就会发现：将图中点黑白染色后，答案一定包括一个黑色点或者白色点的生成树，并且再确定这棵树后，剩下的镜子选法只有一种
• 于是上矩阵树定理即可

code

#include<bits/stdc++.h>
using namespace std;
const int N=310;
int n,m,p,fa[N*N];
char s[N][N];
inline int id(int x,int y){return (x-1)*(m+1)+y;}
inline int find(int x){return fa[x]==x?x:fa[x]=fa[fa[x]]=fa[fa[fa[x]]]=find(fa[fa[fa[x]]]);}
inline void merge(int x,int y){x=find(x);y=find(y);x!=y?fa[x]=y:0;} 
int a[2][N<<2][N<<2],tot[2],bel[N*N];
inline int add(int p,int x,int y){
	++a[p][bel[x]][bel[x]];++a[p][bel[y]][bel[y]];
	--a[p][bel[x]][bel[y]];--a[p][bel[y]][bel[x]]; 
}
inline int inv(int a,int b=p-2){
	int ret=1;
	while(b){
		if(b&1) ret=1ll*ret*a%p;
		a=1ll*a*a%p;b>>=1;
	}
	return ret;
}
inline int matrix_tree(int pt,int x){
	int ans=1;
	for(int j=1;j<=x;++j){
		for(int i=j+1;i<=x;++i){
			if(a[pt][i][j]){
				int d=1ll*a[pt][j][j]*inv(a[pt][i][j])%p;
				for(int k=j;k<=x;++k) a[pt][j][k]=(a[pt][j][k]-1ll*d*a[pt][i][k]%p+p)%p;
				swap(a[pt][i],a[pt][j]);ans=-ans;
			}
		}
		ans=1ll*ans*a[pt][j][j]%p;
		if(ans<0) ans+=p;
	}
	return ans;
}
int main(){
	scanf("%d%d%d",&n,&m,&p);
	for(int i=1;i<=(n+1)*(m+1);++i) fa[i]=i;
	for(int i=1;i<=n;++i){
		scanf("%s",s[i]+1);
		for(int j=1;j<=m;++j){
			if(s[i][j]=='/') merge(id(i+1,j),id(i,j+1));
			else if(s[i][j]!='。') merge(id(i,j),id(i+1,j+1));
		}
	}
	for(int i=1;i<=n+1;++i)
		for(int j=1;j<=m+1;++j)
			if(fa[id(i,j)]==id(i,j)){
				int wh=(i+j)&1;
				bel[id(i,j)]=++tot[wh];
			}
	for(int i=1;i<=n;++i){
		for(int j=1;j<=m;++j){
			if(s[i][j]=='*'){
				add((i+1+j)&1,find(id(i+1,j)),find(id(i,j+1)));
				add((i+j)&1,find(id(i,j)),find(id(i+1,j+1)));
			}
		}
	}
	printf("%d\n",(matrix_tree(0,tot[0]-1)+matrix_tree(1,tot[1]-1))%p);
	return 0;
}

标签：Box,return,int,CF578F,fa,Mirror,inline,find
来源： https://www.cnblogs.com/tqxboomzero/p/13860284.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。