poj3557 Map Generator(概率dp）及其扩展

对于原题：

假如说n个点要联通，我们从反面开始考虑，考虑不联通的情况。比如说1所在的连同块大小为j（j>=1&&j<=n-1)则其他点不管他们是否连通，这j个点与剩下点n－j个点必然不连边，则dp[n]=sigma(dp[j]*C[n-1][j-1]*dp[j]*(1-p)^(j*(n-j))).

同理dp[i]也可以按dp［n］来表示，dp［i］表示i个点连通的概率。

下面考虑其扩展情况：

首先确定顶点的个数N和一个概率p，然后遍历所有的点对(i,j)，如果生成的0到1之间的随机实数小于p，那么就会有一条边连接这两个顶点。那么，请计算整个无向图有K个连通块的概率是多少？

同理k＝1时按上述dp来求解，k不等于1时，dp[i][j]表示i个点有j个连同块的概率。j不等于1时,考虑1所在的连同块有f个点，那么剩下i－f个点则分成j－1个的连同块，而这f个点与i－f个点不连边。

dp[i][j]=sigma(dp[f][1]*C[i-1][f-1]*(1-p)^(f*(i-f))).

j等于1时3557处理即可。

3557代码：

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
double dp[50];
double c[22][22];
void init(void){
    int i,j;
    c[1][1]=1;
    for(i=1;i<=20;i++){
        c[i][0]=1;
    }
    for(i=2;i<=20;i++){
        for(j=1;j<i;j++){
            c[i][j]=c[i-1][j-1]+c[i-1][j];
        }
        c[i][i]=1;
    }
}
int main()
{
    int i,j,n;
    double p;
    init();
    while(scanf("%d%lf",&n,&p)!=EOF){
        memset(dp,0,sizeof(dp));
        for(i=1;i<=n;i++){
            for(j=1;j<=i-1;j++){
                dp[i]+=c[i-1][i-j-1]*pow(1-p,((i-j)*j)*1.0)*dp[i-j];
            }
            dp[i]=1-dp[i];
        }
        printf("%.5lf\n",dp[n]);
    }
}

扩展题代码（如果是BIT的学生，则在网络教室的琪露诺算法教室里有，题目叫：

图生成器

下面是代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<iostream>
using namespace std;
double dp[22][22];
int zuhe[22][22];
void init(void){
    int i,j;
    for(i=0;i<=20;i++){
        for(j=0;j<=i;j++){
            if(j==0) zuhe[i][j]=1;
            else zuhe[i][j]=zuhe[i-1][j]+zuhe[i-1][j-1];
        }
    }
}
int main()
{
    int i,j,n,k,f;
    double p;
    init();
    while(cin>>n>>p>>k){
        memset(dp,0,sizeof(dp));
        dp[1][1]=1;
        for(i=2;i<=n;i++){
            for(j=1;j<=i&&j<=k;j++){
                if(j==1){
                    for(f=1;f<=i-1;f++){
                        dp[i][j]+=zuhe[i-1][f-1]*dp[f][1]*pow(1-p,f*(i-f));
                    }
                    dp[i][j]=1-dp[i][j];
                }
                else{
                    for(f=1;f<=i-1;f++){
                        dp[i][j]+=dp[f][1]*zuhe[i-1][f-1]*dp[i-f][j-1]*pow(1-p,f*(i-f));
                    }
                }
            }
        }
        printf("%.3lf\n",dp[n][k]);
    }
}