hdu5001 鞍山赛区网络赛 E题

来源：互联网发布：周立功单片机北京代理编辑：程序博客网时间：2024/04/28 15:45

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

每轮模拟后把特定点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 0x3f3f3f3fint 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;}

0 0