[hdu5304]Eastest Magical Day Seep Group‘s Summer

Eastest Magical Day Seep Group's Summer

题解

看到题面时应该很容易想到矩阵树定理，但是由于是一个条边的连通图，并不能直接用矩阵树来处理。

很容易发现，任意一个个点条边的连通图，一定是一棵基环树，于是我们可以考虑缩点。

考虑状压，表示一根起点为的链状态为的图的构成方案数，可以通过背包来进行更新。

之后枚举每个，再计算出每个状态为的环的构成方案数，将所有在中的点缩点后再将该图的基尔霍夫矩阵建出来，计算出它的行列式值，加上它的贡献即可。由于每个环有两个方向，还需要除个2。

时间复杂度，按道理说是可能T的，但还是卡过去了。

源码

#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<vector>
#include<queue>
using namespace std;
#define MAXN 20
typedef long long LL;
const LL INF=0x7f7f7f7f;
const LL mo=998244353;
const LL inv2=499122177;
typedef pair<int,int> pii;
template<typename _T>
_T Fabs(_T x){return x<0?-x:x;}
template<typename _T>
void read(_T &x){
	_T f=1;x=0;char s=getchar();
	while(s>'9'||s<'0'){if(s=='-')f=-1;s=getchar();}
	while('0'<=s&&s<='9'){x=(x<<3)+(x<<1)+(s^48);s=getchar();}
	x*=f;
}
bool Zero(LL x){return x==0;}
LL qkpow(LL a,int s){
	LL t=1;
	while(s){
		if(s&1)t=1ll*a*t%mo;
		a=1ll*a*a%mo;s>>=1;
	}
	return t;
}
LL inv(LL x){return qkpow(x,mo-2);}
LL add(LL x,LL y){return x+y>=mo?x+y-mo:x+y;}
struct matrix{
	int c[MAXN][MAXN];int n;
	matrix(){memset(c,0,sizeof(c));n=0;}
	void Guass(){
		for(int i=1,mx=1;i<=n;i++,mx=i){ 
			for(int j=i;j<=n;j++)if(!Zero(c[j][i])){mx=j;break;}
			if(mx!=i)for(int j=i;j<=n;j++)swap(c[mx][j],c[i][j]);
			if(Zero(c[i][i]))continue;
			for(int j=i+1;j<=n;j++){
				LL tmp=1ll*c[j][i]*inv(c[i][i])%mo;if(Zero(tmp))continue;
				for(int k=i;k<=n;k++)c[j][k]=add(c[j][k],mo-1ll*c[i][k]*tmp%mo);
			}
		}
	}
	LL query(){
		LL res=1;
		for(int i=1;i<=n;i++)res=1ll*res*c[i][i]%mo;
		return res;
	}
	void clear(){memset(c,0,sizeof(c));n=0;}
	void print(){
		puts("\n");
		for(int i=1;i<=n;i++,puts(""))
			for(int j=1;j<=n;j++)
				printf("%d ",c[i][j]);
		puts("\n");
	}
};
int mp[MAXN][MAXN],id[MAXN],cnt;
LL sum[20][(1<<17)+5];
struct Kirchhoff{
	matrix ansMat;
	LL solve(int s,int nn){
		ansMat.clear();cnt=1;
		for(int i=1;i<=nn;i++)if(s&(1<<i-1))id[i]=cnt; 
		for(int i=1;i<=nn;i++)if(!(s&(1<<i-1)))id[i]=++cnt;
		for(int i=1;i<=nn;i++)
			for(int j=i+1;j<=nn;j++){
				if(!mp[i][j])continue;int u=i,v=j;
				ansMat.c[id[u]][id[u]]++;ansMat.c[id[v]][id[v]]++;
				ansMat.c[id[u]][id[v]]--; ansMat.c[id[v]][id[u]]--;
			}
		if(nn-cnt<2)return 0;
		ansMat.n=cnt-1;ansMat.Guass();
		return (ansMat.query()+mo)%mo;
	}
}Tree;
signed main(){
	int n,m;
	while(scanf("%d %d",&n,&m)!=EOF){
		memset(mp,0,sizeof(mp));
		for(int i=1;i<=n;i++)
			for(int j=0;j<(1<<n);j++)
				sum[i][j]=0;
		for(int i=1;i<=m;i++){
			int u,v;read(u);read(v);
			mp[u][v]=mp[v][u]=1;
		}
		for(int i=1;i<=n;i++){
			sum[i][1<<i-1]=1;
			for(int s=(1<<i-1);s<(1<<n);s++)
				if((((1<<i-1)-1)&s)==0){
					for(int j=i;j<=n;j++)if(sum[j][s])
						for(int k=i+1;k<=n;k++)
							if(((1<<k-1)&s)==0&&mp[j][k])
								sum[k][s|(1<<k-1)]=add(sum[k][s|(1<<k-1)],sum[j][s]);
					//if(sum[i][s])printf("%d %d:%d\n",i,s,sum[i][s]);
				}
		}
		LL ans=0;
		for(int s=1;s<(1<<n);s++){
			LL tmp=0;int k=0;
			for(int i=1;i<=n;i++)if((1<<i-1)&s){k=i;break;}
			for(int i=k+1;i<=n;i++)
				if(((1<<i-1)&s)&&mp[i][k])
					tmp=add(tmp,sum[i][s]);
			if(!tmp)continue;//printf("%d:%d\n",s,tmp);
			ans=add(ans,1ll*Tree.solve(s,n)*tmp%mo*inv2%mo);
		}
		printf("%lld\n",ans);
	}
	return 0;
}

谢谢！！！