udebug通过但UVaWA,求助
查看原帖
udebug通过但UVaWA,求助
660645
pengpengXiuEr楼主2022/6/22 23:55
#include<cstdio>
#include<cstring>
#include<vector>
#include<map>//强制性 
#include<algorithm>
using namespace std;
int T,n,m,f[32],k[32],p,q;//蓝题 
vector<int> v[32];
map<vector<int>,long double> ma;
int head(int x){
	return f[x]==x?x:f[x]=head(f[x]);
}
void sor(int x){
	int b=v[x].size();
	for(int i=0;i!=b;i++){
		f[i]=v[x][i];
	}
	sort(f+1,f+b);
	v[x].clear();
	for(int i=0;i!=b;i++){
		v[x].push_back(f[i]);
	}
	return;
}
long double dp(long double u,int w){
	if(ma[v[w]]){
		long double c=ma[v[w]]*u;//注意引用 
		v[w].clear();
	    return c;
	}
	if(v[w].size()==1){
	    v[w].clear();
		return 0; 
	} 
	long double ans=0;
	for(int i=v[w].size()-1;i;i--){
		for(int j=0;j!=v[w].size();j++){
			if(i==j){
				v[w+1][0]+=v[w][j];
				continue;
			}
			v[w+1].push_back(v[w][j]);
		}
		ans+=dp(u*v[w][i]/(n-v[w][0]),w+1)+u*v[w][i]/(n-v[w][0]);
	}
	ans+=u*(v[w][0]-1)/(n-v[w][0]);
	ma[v[w]]=ans/u;
	v[w].clear();
	return ans;
}
signed main(){
	scanf("%d",&T);
	for(int r=1;r<=T;r++){
		memset(k,0,sizeof(k));
		v[0].clear();
		ma.clear();
		scanf("%d%d",&n,&m);
		for(int i=1;i<=n;i++){
			f[i]=i;
		}
		while(m--){
			scanf("%d%d",&p,&q);
			if(p>q)p^=q^=p^=q;
			f[q]=p;
		}
		for(int i=1;i<=n;i++){
			k[head(i)]++;
		}
		for(int i=1;i<=n;i++){
			if(k[i])v[0].push_back(k[i]);
		}
		sor(0);
		printf("Case %d: %.16Lf\n",r,dp(1.0,0));//不能共存的节日
	}
	return 0;
} 
2022/6/22 23:55
加载中...