【例题】【数学期望】NKOJ3853 constandine

NKOJ3853 constandine
时间限制 : - MS 空间限制 : 165536 KB
评测说明 : 2s

问题描述
古月族是一个神奇的种族，有一天古月腾和古月豪在一起看《地狱神探(constandine)》。
剧中，小康是一个驱魔师。有一天，他找到了一本驱魔指南，上面记录了许多咒语，每个咒语由许多字符组成，例如：“!?####T_Tabc”。
仪式一开始，需要呼告，咒语总是从某个字符串st开始，最后总是以某个字符串en结束。咒语之间存在着一种古老的联系，如念完咒语i后，接着只能念k条规定的咒语(保证从任何一个咒语开始念都能念到en，并且从st开始，能念到任何咒语，后面的咒语不能念到前面的，即不能倒着念)，否则会破坏自然的平衡，反噬自身，变成恶魔。每两个咒语连着念会产生一定的魔法值Wij。
念到一句咒语时（还没念到到en），会面临k种选择来完成整个仪式，并且每种选择的概率为1/k。
小康为了提前准备仪式，以防走火入魔，他想让你帮帮他，算出整个咒语的魔法值的期望值是多少？
四舍五入，保留两位小数。

输入格式
第一行为两个正整数n，m，表示有n条咒语(保证字符串长度不超过5)，有m种联系。

第二行到第m+1行，每行两个字符串，一个正整数，a，b，c，代表咒语a与咒语b有联系，且连续念产生的魔法值为c（输入出现的第一个不同的字符串为st，k==0的字符串为en）。

输出格式
输出整个咒语的魔法值的期望值，四舍五入保留两位小数。

样例输入
4 4
ah aaaaa 1
ah hjh 2
aaaaa hjh 3
hjh ht 4

样例输出
7.00

提示
【数据范围】

对于20%的数据，n<=100

对于另20%的数据，n<=1000

对于另20%的数据，n<=10000

对于100%的数据，n<=100000，m<=2*n

节点到e的期望值为各下家到e的期望乘上该节点到该下家的概率，该题中概率为1/k;

#include<cstdio>
#include<vector>
#include<map>
#include<iostream>
using namespace std;
#define END fclose(stdin),fclose(stdout);return 0
const int needn=100003;

int n,m;

//..........................................................
inline void in_(int & d)
{
    char t=getchar();
    while(t<'0'||t>'9') t=getchar();
    for(d=0;!(t<'0'||t>'9');t=getchar()) d=(d<<1)+(d<<3)+t-'0';
}
//..........................................................
map<string,int> ma;
int s=1,tot_,e;
int chu[needn];

int addma(const string & a)
{
    if(ma.find(a)==ma.end()) return ma[a]=++tot_;
    else return ma[a];
}
//..........................................................
int tot,fi[needn],la[needn<<1],en[needn<<1],len[needn<<1];

void add(int a,int b,int c)
{
    //cout<<a<<" "<<b<<endl;
    tot++;
    la[tot]=fi[a];
    fi[a]=tot;
    en[tot]=b;
    len[tot]=c;
}
//..........................................................
long long cnt,all;
double ans[needn];

double dfs(int x)
{
    if(ans[x]!=-1.0) return ans[x];
    if(x==e) return 0;
    double anss=0.0;
    int son=0;
    for(int t=fi[x],y;t;t=la[t])
    {
        y=en[t];
        son++;
        anss+=dfs(y)+(double)len[t];
    }
    return ans[x]=anss/son;
}
//..........................................................

int main()
{
    scanf("%d%d",&n,&m);
    string a,b;int c,x,y;
    for(int i=1;i<=m;i++) 
    {
        cin>>a>>b;in_(c);
        x=addma(a);y=addma(b);
        add(x,y,c);
        chu[x]++;
    }
    for(int i=1;i<=n;i++) 
    {
        ans[i]=-1.0;
        if(chu[i]==0)
        {
            e=i;
            break;
        } 
    }
    printf("%.2lf",dfs(s));
}