【问题描述】
随机生成一张n个点的无向图。把他顶点标号为1~n。
对于一个点对i,j(1 <= i < j <= n),他有千分之p的概率成为这张图中的一条边(i,j)。不同的边出现的概率是相互独立的。
那么生一个至少含有一个大于等于4的连通块的图的概率是多少。
【输入格式】
一行两个整数n,p,表示图的点数和边在图上的概率。
【输出格式】
输出一行一个实数,表示所求的概率。结果保留4位小数。
【输入样例】
6 100
【输出样例】
0.1184
【样例解释】
时间限制:1秒 内存限制:128M
【数据范围】
4<=n<=100
嗯,一道概率+动态规划的好题。
题目要求生成一个至少含有一个大于等于4的连通块的图的概率是多少。如果我们从正面入手,直接求解,会发现很难往下做(“至少”并没有表明到底有多少);我们可以从反面入手,所谓正难则反易,这么看答案就很容易得出了——它等于1减去“所有连通块都小于等于3”的概率。
状态函数定义:dp(i,j,k) 表示“前i个点中,含有2个点的连通分量数为j,含有3个点的连通分量数为k,含有1个点的连通分量数为i-2*j-3*k”的概率。下一步是状态转移,用刷表法,免不了分类讨论:
情况1:第i+1个点不与前面i个点中的任意一个点相连,则状态转移方程为dp(i+1,j,k)+=dp(i,j,k)*(1-p)^i
情况2:第i+1个点与前面i个点中的某一个孤立点连一条边,则状态转移方程为dp(i+1,j+1,k)+=dp(i,j,k)*C(i-2*j-3*k,1)*p*(1-p)^(i-1)
情况3:第i+1个点与前面i个点中的某两个孤立点各连一条边,则状态转移方程为dp(i+1,j,k+1)+=dp(i,j,k)*C(i-2*j-3*k,2)*p^2*(1-p)^(i-2)
情况4:第i+1个点与前面i个点中的某个含有2个点的连通分量连一条边,则状态转移方程为dp(i+1,j-1,k+1)+=dp(i,j,k)*2*C(j,1)*p*(1-p)^(i-1)
情况5: 第i+1个点与前面i个点中的某个含有2个点的连通分量中的两个点分别连一条边,则状态转移方程为dp(i+1,j-1,k+1)+=dp(i,j,k)*C(j,1)*p^2*(1-p)^(i-2)
边界分析:dp(0,0,0)=1
答案分析:Ans=1-∑dp(n,x,y) (0<=x<=n/2 && 0<=y<=n/3 && 2*x+3*y<=n)
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#define dp(x,y,z) dp[x][y][z]
using namespace std;
int n,p;
double P,dp[105][55][35],qkpow1[105],qkpow2[105];
void ready(int nn)
{
qkpow1[0]=qkpow2[0]=1;
for(int i=1;i<=nn;i++)
{
qkpow1[i]=qkpow1[i-1]*P;
qkpow2[i]=qkpow2[i-1]*(1-P);
}
return;
}
void init()
{
scanf("%d%d",&n,&p);
P=(double)p/1000;
ready(n);
return;
}
void DP()
{
memset(dp,0,sizeof(dp));
dp[0][0][0]=1;
for(int i=0;i<n;i++)
for(int j=0;2*j<=i;j++)
for(int k=0;2*j+3*k<=i;k++)
{
int l=i-2*j-3*k;
dp(i+1,j,k)+=dp(i,j,k)*qkpow2[i];
dp(i+1,j+1,k)+=dp(i,j,k)*l*P*qkpow2[i-1];
dp(i+1,j,k+1)+=dp(i,j,k)*l*(l-1)/2*qkpow1[2]*qkpow2[i-2];
if(j>=1)
dp(i+1,j-1,k+1)+=dp(i,j,k)*(2*j*P*qkpow2[i-1]+j*qkpow1[2]*qkpow2[i-2]);
}
double Ans=1;
for(int x=0;2*x<=n;x++)
for(int y=0;2*x+3*y<=n;y++)
Ans-=dp(n,x,y);
printf("%0.4lf\n",Ans);
return;
}
int main()
{
// freopen("in.txt","r",stdin);
// freopen("out.txt","w",stdout);
init();
DP();
return 0;
}