hdu5001 鞍山赛区网络赛 E题

最新推荐文章于 2015-09-04 23:34:27 发布

cqsh3vj2

最新推荐文章于 2015-09-04 23:34:27 发布

阅读量425

点赞数

分类专栏：概率dp

本文链接：https://blog.csdn.net/cqsh3vj2/article/details/39253799

版权

概率dp 专栏收录该内容

13 篇文章 0 订阅

订阅专栏

从一个状态到下一个状态是存在一个转移矩阵的，但是转移矩阵应该不如直接模拟的复杂度低

每轮模拟后把特定点k清0，最后把k点以外的各个点的概率加和

原本考虑到矩阵的收敛性，添加了一个收敛性判断（精度要求比较高才能过），1.8s

后来去掉收敛性判断直接暴力反而更快，1.2s

可能是因为d本身比较小，精度要求有些高，没有明显的区别

#include<iostream>
#include<algorithm>
#include<iostream>
#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<vector>
#include<queue>
#include<cmath>
using namespace std;
#define ll __int64
#define inf 0x3f3f3f3f
int dir[4][2]={0,1,0,-1,1,0,-1,0};

int n,m,d;
vector<int>bians[60];
double map[60][2];

int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        int i,j,k,p,num;
        int ta,tb;
        int te,le;
        int flag;
        double ans;
        scanf("%d%d%d",&n,&m,&d);
        for(i=1;i<=n;i++)
            bians[i].clear();
        for(i=1;i<=m;i++)
        {
            scanf("%d%d",&ta,&tb);
            bians[ta].push_back(tb);
            bians[tb].push_back(ta);
        }
        for(k=1;k<=n;k++)
        {
            for(i=1;i<=n;i++)
                map[i][0]=1.0/(double)n;
            map[k][0]=0;
            for(p=1;p<=d;p++)
            {
                te=p%2;
                le=(p+1)%2;
                for(i=1;i<=n;i++)
                    map[i][te]=0;
                for(i=1;i<=n;i++)
                {
                    num=bians[i].size();
                    for(j=0;j<num;j++)
                        map[bians[i][j]][te]+=map[i][le]/(double)num;
                }
                map[k][te]=0;
            }
            ans=0;
            for(i=1;i<=n;i++)
                ans+=map[i][te];
            printf("%.10lf\n",ans);
        }
    }
    return 0;
}