刷题总结——随机图（ssoi）

题目：

随机图
(random.cpp/c/pas)
【问题描述】
BG 为了造数据，随机生成了一张�个点的无向图。他把顶点标号为1~�。
根据BG 的随机算法，对于一个点对�, �(1 ≤ � < � ≤ �)，它有�‰(千分之�)的概率成为
这张图中的一条边� − �。不同的边出现的概率是相互独立的。
为了保证数据的强度，BG要求生成的图中至少要有一个大于等于 4的连通块。于是他
想知道，在不做任何改进的情况下，根据他的算法造出强数据的概率是多少。

【输入】
输入文件名为random.in。
共一行两个整数�, �，中间用一个空格隔开，表示图的点数和边在图上的概率。

【输出】
输出文件名为random.out。
输出共一行一个实数，表示所求的概率。结果保留4 位小数。 【输入输出样例1】

random.in random.out
3 123

0.0000

【输入输出样例2】
random.in random.out
4 500

0.5938

题解：�

概率DP方程

�[�+1][�−1][�+1]，乘上系�·(1−�)�

代码：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<cctype>
#include<cstring>
#include<string>
#include<algorithm>
using namespace std;
const int N=;
double dp[N][N][N],p,ans=;
double c[N][];
int n;
inline void pre()
{
  for(int i=;i<=n;i++)  c[i][]=i;
  for(int i=;i<=n;i++)
    for(int j=;j<=n;j++)
      c[i][j]=c[i-][j]+c[i-][j-];
}
int main()
{
  scanf("%d",&n);
  cin>>p;
  p=p/;
  pre();
  dp[][][]=;
  for(int i=;i<=n;i++)
    for(int j=;j<=n/;j++)
      for(int k=;k<=n/;k++)
      {
        int l=(i-*j-*k);
        dp[i+][j][k]+=dp[i][j][k]*pow(-p,i);
        dp[i+][j+][k]+=dp[i][j][k]*pow(-p,i-)*l*p;
        dp[i+][j][k+]+=dp[i][j][k]*c[l][]*pow(-p,i-)*p*p;
        dp[i+][j-][k+]+=dp[i][j][k]*pow(-p,i-)**j*p;
        dp[i+][j-][k+]+=dp[i][j][k]*pow(-p,i-)*j*p*p;
      }
  for(int j=;j<=n/;j++)
    for(int k=;k<=n/;k++)
      ans-=dp[n][j][k];
  printf("%0.4lf",ans);
  return ;
}

� ≤ 1,000。