大致题意: 有\(n\)个点和\(m\)条边,每条边的权值是一个\(0\sim1\)的随机实数,要你用\(n-1\)条边将图联通,问这\(n-1\)条边中边权最大值的期望最小值。
提示
这题应该是一道比较难的\(DP\)题吧。
首先,我们需要注意到提示中的一句话:
\(Hint\)
对于\(n\)个\([0,1]\)之间的随机变量\(x_1,x_2,...,x_n\),第\(k\)小的那个的期望值是\(\frac k{n+1}\)。
其实,这就很明显在提示我们,只要求出这\(n-1\)条边中权值最大边的期望最小排名,我们就能轻松求出它的期望最小值。
于是我们就能形成一个比较基础的思路了。
基础思路
首先,应该不难想到,可以用\(P_i\)来表示权值最大边排名为\(i\)的概率,而答案就是\(\sum_{i=1}^mi·P_i\)(这应该还是比较显然的)。
但问题在于这样很难进行状态转移。
我们可以将原式进行进一步的转化,则不难发现,答案其实就相当于\(\sum_{i=1}^m\sum_{j=i}^mP(j)\)。
由于\(\sum_{i=1}^mP(i)=1\),所以上面式子中的\(\sum_{j=i}^mP(j)\)其实又等价于\(1-\sum_{j=1}^{i-1}P(j)\)。
考虑\(\sum_{j=1}^{i-1}P(j)\)表示的是加入\(i-1\)条边能使原图联通的概率,所以\(1-\sum_{j=1}^{i-1}P(j)\)其实就相当于加入\(i-1\)条边不能使原图联通的概率。
不难发现,这比求原式要容易多了。
转化后的问题
整理一下上面的思路,就可以得到一个这样的问题:设\(P'(i)\)为加入\(i\)条边后原图不连通的概率,求\(\sum_{i=0}^{m-1}P'(i)\)。
注意到\(n\le10\)这极小的范围,我们首先想到的自然是乱搞状压\(DP\)。
状压\(DP\)
考虑用\(f_{i,t,0/1}\)来表示连接\(i\)条边后使点集\(j\)不连通/联通的方案数,注意,这里记录的是方案数而不是概率,这样可以防止出现精度误差(毕竟要保留\(6\)位小数)。
我们可以用\(Size_i\)和\(line\_tot_i\)分别表示点集\(i\)中点与边的数量。
首先,我们要知道一个比较基础的知识,即在点集\(i\)中选择\(t\)条边的方案数应为\(C_{line\_tot_i}^t\),即\(f_{i,t,0}+f_{i,t,1}=C_{line\_tot_i}^t\)。
既然这样,只要求出\(f_{i,t,0}\)和\(f_{i,t,1}\)中的一个,我们就可以轻松求出另一个的值了。
比较显然\(f_{i,t,0}\)相对而言要更容易求一些。
状态转移
我们来考虑如何进行状态转移。
考虑在当前点集中任意选取一点\(k\)作为定点(\(k\)取什么点无所谓,为了方便推荐选择\(lowbit(i)\),可以\(O(1)\)求出),然后枚举含有点\(k\)的点集\(j\),进行转移。
不难想到,如果要让点集\(j\)变成\(i\),自然需要加入\(i\)^\(j\)这个点集。
我们可以在保证点集\(j\)联通的情况下,枚举\(x,y(0\le x\le line\_tot_j,0\le y\le line\_tot_{i\text{^}j})\),表示两个点集中的分别选择的边数。
由于选\(x\)条边使点集\(j\)联通的方案数为\(f_{j,x,1}\),而\(y\)条边的选择方案数共有\(C_{line\_tot_{i\text{^}j}}^y\),所以可以得出转移方程:
\]
这样就可以不重不漏地统计到所有情况了。
最后的答案
根据一开始得出的结论,最后的答案就是:
\]
还有,千万注意精度啊。
一些细节:关于\(Size\)与\(line\_tot\)数组的求法
关于\(Size\)和\(line\_tot\)两个数组,其实有一种比较巧妙的求法。(会的大佬可以直接跳过这一部分)
- \(Size_i=Size_{i>>1}+(i\&1)\)
- 首先,由于\((i>>1)<i\),所以\(Size_{i>>1}\)一定在求\(Size_i\)之前已经求解过。
- 其次,由于\((i>>1)\)与\(i\)在二进制下只相差最后一位,所以只要比较最后一位(即\(i\&1\))就可以求出结果。
- \(line\_tot_i=\frac{\sum Size_{Next_j\&i}}2\)
- 注意,上面式子中的\(j\)表示点集\(i\)中的一个点,而\(Next_j\)表示与\(j\)有边相连的点的点集。
- 首先,由于\(Next_j\&i\le i\),所以\(Size_{Next_j\&i}\)肯定已经求解过。
- 其次,将与\(j\)相连的点与点集\(i\)中的点\(\&\)一下,得到的结果就是点集\(i\)中与\(j\)相连的点的点集,它的\(Size\)就是点集\(i\)中与\(j\)有边相连的点的数量,即以点\(j\)为一个端点的边数。
- 但这样会把每条边统计两次,因此最后除以\(2\)即可。
代码
#include<bits/stdc++.h>
#define N 10
#define M 50
#define LL long long
using namespace std;
int n,m,Next[N+5],line_tot[1<<N+1],Size[1<<N+1];LL f[1<<N+1][M+5][2],C[M+5][M+5];
class Class_FIO
{
private:
#define Fsize 100000
#define tc() (A==B&&(B=(A=Fin)+fread(Fin,1,Fsize,stdin),A==B)?EOF:*A++)
char ch,*A,*B,Fin[Fsize];
public:
Class_FIO() {A=B=Fin;}
inline void read(int &x) {x=0;while(!isdigit(ch=tc()));while(x=(x<<3)+(x<<1)+(ch&15),isdigit(ch=tc()));}
}F;
inline void Init() {for(register int i=C[0][0]=1,j;i<=M;++i) for(C[i][0]=C[i][i]=j=1;j<i;++j) C[i][j]=C[i-1][j-1]+C[i-1][j];}//预处理组合数
int main()
{
register int i,j,k,x,y,Full;register double ans=0.0;
for(Init(),F.read(n),F.read(m),i=1;i<=m;++i) F.read(x),F.read(y),Next[x]|=(1<<y-1),Next[y]|=(1<<x-1);//用Next[i]存储与i有边相连的点的点集
for(i=1,Full=(1<<n)-1;i<=Full;++i)//枚举每一个点集i
{
if((Size[i]=Size[i>>1]+(i&1))==1) {f[i][0][1]=1;continue;}//如果只有一个点,跳过
for(j=1;j<=n;++j) if(i&(1<<j-1)) line_tot[i]+=Size[Next[j]&i];line_tot[i]>>=1;//统计边数
for(k=i&-i,j=i&(i-1);j;j=i&(j-1)) if(j&k) for(x=0;x<=line_tot[j];++x) for(y=0;y<=line_tot[i^j];++y) f[i][x+y][0]+=f[j][x][1]*C[line_tot[i^j]][y];//状态转移
for(j=0;j<=line_tot[i];++j) f[i][j][1]=C[line_tot[i]][j]-f[i][j][0];//根据f[i][j][0]得到f[i][j][1]的值
}
for(i=0;i<=m;++i) ans+=1.0*f[Full][i][0]/C[line_tot[Full]][i];ans/=m+1;//统计答案
return printf("%.6lf",ans),0;
}