[ZJOI2015] 地震后的幻想乡

关注

原创最新推荐文章于 2022-09-12 20:34:57 发布 1k 阅读 AI 写同款· GEO 优化 ›

【BZOJ】3925 [Zjoi2015]地震后的幻想乡状压+期望DP||定积分题目传送门这题是真的神仙题……整整花了我两个礼拜来理解这题首先这题据我了解有三种做法：纯OI做法、积分+数学推导、直接积分请做好一定的心理准备，接下来的东西可能有点难理解~~（好像不是一点点的难吧……）~~ 1.纯OI做法首先我们根据期望的线性性，可以得到ans=∑x=0mw(x)×p(x)ans=\sum_{x=0}^mw(x)\times p(x)ans=∑x=0mw(x)×p(x... 阅读详情

P3255【ZJOI2015 Day1】地震后的幻想乡

时间限制 : 20000 MS 空间限制 : 262144 KB

问题描述

        傲娇少女幽香是一个很萌很萌的妹子，而且她非常非常地有爱心，很喜欢为幻想乡的人们做一些自己力所能及的事情来帮助他们。

        这不，幻想乡突然发生了地震，所有的道路都崩塌了。现在的首要任务是尽快让幻想乡的交通体系重新建立起来。幻想乡一共有n个地方，那么最快的方法当然是修复n-1条道路将这n个地方都连接起来。

        幻想乡这n个地方本来是连通的，一共有m条边。现在这m条边由于地震的关系，全部都毁坏掉了。每条边都有一个修复它需要花费的时间，第i条边所需要的时间为ei。地震发生以后，由于幽香是一位人生经验丰富，见得多了的长者，她根据以前的经验，知道每次地震以后，每个ei会是一个0到1之间均匀分布的随机实数。并且所有ei都是完全独立的。

        现在幽香要出发去帮忙修复道路了，她可以使用一个神奇的大魔法，能够选择需要的那n-1条边，同时开始修复，那么修复完成的时间就是这n-1条边的ei的最大值。当然幽香会先使用一个更加神奇的大魔法来观察出每条边ei的值，然后再选择完成时间最小的方案。

        幽香在走之前，她想知道修复完成的时间的期望是多少呢？

输入格式

        第一行两个数n,m，表示地方的数量和边的数量。其中点从1到n标号。

        接下来m行，每行两个数a,b，表示点a和点b之间原来有一条边。

        这个图不会有重边和自环。

输出格式        一行输出答案，四舍五入保留6位小数。



样例输入样例输入1：

5 4

1 2

1 5

4 3

5 3



样例输入2：

7 10

2 5

1 4

5 6

3 7

5 3

7 1

6 2

4 6

4 3

7 6

样例输出样例输出1：

0.800000



样例输出2：

0.617027

提示提示：

        （以下内容与题意无关，对于解题也不是必要的。）

        对于n个[0,1]之间的随机变量x1,x2,...,xn，第k小的那个的期望值是。



样例解释：

        对于第一个样例，由于只有4条边，幽香显然只能选择这4条，那么答案就是4条边的ei中最大的数的期望，由提示中的内容，可知答案为0.8。



数据范围：

        对于所有数据：n<=10, m<=n(n-1)/2, n,m>=1。

        对于15%的数据：n<=3。

        另有15%的数据：n<=10, m=n。

        另有10%的数据：n<=10, m=n(n-1)/2。

        另有20%的数据：n<=5。

        另有20%的数据：n<=8。

前面的分析这位大神分析的非常到位：







那么最后答案是什么呢？1m+1∑i=0mf[all][i](icnt[all])




代码：

#include<cstdio>
#include<iostream>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<queue>
#include<iomanip>
#define ll long long
using namespace std;
const int inf=1000000000;
template <typename T>
inline void _read(T& x){
    char t=getchar();bool sign=true;
    while(t<'0'||t>'9'){if(t=='-')sign=false;t=getchar();}
    for(x=0;t>='0'&&t<='9';t=getchar())x=x*10+t-'0';
    if(!sign)x=-x;
}
int n,m;
int cnt[1<<11],Set[15];
ll c[55][55],f[1<<11][55],g[1<<11][55];
int bitcnt(int x){
	int temp=0;
	while(x){
		temp++;
		x-=x&(-x);
	}
	return temp;
}
int main(){
	int i,j,k,s,s1;
	cin>>n>>m;
	for(i=1;i<=m;i++){
		int x,y;
		_read(x);_read(y);
		x--;y--;
		Set[x]|=(1<<y);
		Set[y]|=(1<<x);
	}
	for(i=0;i<=m;i++){
		for(j=0;j<=i;j++){
			if(j==0)c[i][j]=1;
			else c[i][j]=c[i-1][j]+c[i-1][j-1];
		}
	}
	int all=(1<<n)-1;
	for(s=1;s<=all;s++){
		int size=bitcnt(s);
		if(size==1){
			g[s][0]=1;
			continue;
		}
		for(i=0;i<n;i++){
			if(s&(1<<i)){
				cnt[s]+=bitcnt(Set[i]&s);
			}
		}
		cnt[s]>>=1;
		int p=s&(-s);
		for(s1=(s-1)&s;s1;s1=(s1-1)&s){
			if(s1&p){
				for(i=0;i<=cnt[s1];i++){
					for(j=0;j<=cnt[s-s1];j++){
						f[s][i+j]+=g[s1][i]*c[cnt[s-s1]][j];
					}
				}
			}
		}
		for(i=0;i<=cnt[s];i++){
			g[s][i]=c[cnt[s]][i]-f[s][i];
		}
	}
	double ans=0.0;
	for(i=0;i<=m;i++){
		ans+=(double)f[all][i]/(double)c[cnt[all]][i];
	}
	ans/=double(m+1);
	cout<<fixed<<setprecision(6)<<ans;
}