Vijos 3747 随机图

【问题描述】 　　

随机生成一张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;
}