hdu5001（概率dp）

题目连接：http://acm.hdu.edu.cn/showproblem.php?pid=5001

题意：一个人随即从一个点出发，到达邻接点的概率相同，求出走d步都不会到达1~n点的每一点i的概率。

分析：每次删掉i点去计算d步后到达其他点的概率，加起来的都是没有到达i点的概率。

一开始算出x步到达每点的概率，然后d步后（1-dp[x][i]）^d，错了，因为算每点到达的概率时已经算了到达i点时的概率了。

网络赛时这道水题做不出，好殇。。。。。

#include <cstdio>

#include <cstring>

#include <string>

#include <cmath>

#include <iostream>

#include <algorithm>

#include <queue>

#include <cstdlib>

#include <stack>

#include <vector>

#include <set>

#include <map>

#define LL long long

#define mod 100000000

#define inf 0x3f3f3f3f

#define eps 1e-9

#define N 100010

#define FILL(a,b) (memset(a,b,sizeof(a)))

#define lson l,m,rt<<1

#define rson m+1,r,rt<<1|1

using namespace std;

double dp[][];

vector<int>g[];

int main()

{

    int t,n,m,d;

    scanf("%d",&t);

    while(t--)

    {

        scanf("%d%d%d",&n,&m,&d);

        for(int i=;i<=n;i++)g[i].clear();

        for(int i=;i<=m;i++)

        {

            int u,v;

            scanf("%d%d",&u,&v);

            g[u].push_back(v);

            g[v].push_back(u);

        }

        for(int i=;i<=n;i++)

        {

            FILL(dp,);

            for(int j=;j<=n;j++)

                if(j!=i)dp[][j]=1.0/n;

            for(int x=;x<=d;x++)

            {

                for(int j=;j<=n;j++)

                {

                    if(i==j)continue;

                    for(int k=,sz=g[j].size();k<sz;k++)

                    {

                        int v=g[j][k];

                        if(v!=i)dp[x][v]+=dp[x-][j]/sz;

                    }

                }

            }

            double res=;

            for(int j=;j<=n;j++)

                if(i!=j)res+=dp[d][j];

            printf("%.10lf\n",res);

        }

    }

}

巴特西

hdu5001（概率dp）

最新文章

热门文章