图

Description

这个世界是一个有向图，图中有n个点m条边且无重边无自环，每秒第i条边出现的概率是p[i]/100，一开始Samjia在1点，每一秒假设Samjia在点x上，那么Samjia要从存在的边中选一条来走，不可以不走，如果不存在可以走的边，那么Samjia就会掉出这个世界，假设Samjia绝顶聪明，问最后Samjia可以成功到达n的概率是多少。
n<=50

Solution

首先，假设你像Samjia一样绝顶聪明，你可以预先知道每个点走到n的概率。
那么对于当前的状态，你肯定是优先走到概率最大的点，如果不能走到就走到次大的。
但是我们并不像Samjia一样绝顶聪明，于是我们可能需要猜测每个点的相对大小。
暴力的做法我们可以枚举相对大小然后列方程高斯消元，然后判断是否合法。
但是这样做显然不够优秀，我们可以考虑调整。
先随机一个相对大小，消元解出一个新的相对大小，如果一样就证明找到了解，否则就用这个新的相对大小迭代进行。
复杂度我也不知道是多少反正能跑过_ (:з」∠) _

Code

#include <cmath>
#include <vector>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;

typedef double db;

const int N=55;

int n,m,x,y,p,id[N],Id[N];
db A[N][N],ans,eps=1e-6;
bool vis[N];
vector<int> v[N],to[N];

#define pb(a) push_back(a)

struct node{
    int id,to;
    db p;
}a[N];

bool cmp(node x,node y) {return x.id<y.id;}
bool cmp1(node x,node y) {return x.p>y.p;}

void gauss() {
    fo(i,1,n) {
        fo(j,i,n)
            if (fabs(A[j][i])>fabs(A[i][i])) 
                fo(k,i,n+1) swap(A[j][k],A[i][k]);
        if (fabs(A[i][i])<eps) continue;
        fo(j,i+1,n) {
            db d=A[j][i]/A[i][i];
            fo(k,i,n+1) A[j][k]=A[i][k]*d-A[j][k];
        }
    }
    fd(i,n,1) {
        if (fabs(A[i][i])<eps) continue;
        fd(j,i-1,1) {
            db d=A[j][i]/A[i][i];
            fo(k,j,n+1) A[j][k]=A[i][k]*d-A[j][k];
        }
    }
}

void get_id() {
    fo(i,1,n-1) a[i].p=A[i][n+1]/A[i][i],a[i].id=i;
    sort(a+1,a+n,cmp1);
    fo(i,1,n-1) Id[a[i].id]=i;
}

bool same() {
    fo(i,1,n-1) if (Id[i]!=id[i]) return 0;
    return 1;
}

void solve() {
    fo(i,1,n) fo(j,1,n+1) A[i][j]=0;
    fo(i,1,n) {
        if (i==n) {
            A[i][i]=1.0;
            A[i][n+1]=1.0;
            continue;
        }
        for(int j=0;j<to[i].size();j++) {
            a[j].id=id[to[i][j]];
            a[j].to=to[i][j];
            a[j].p=v[i][j]/100.0;
        }
        sort(a,a+to[i].size(),cmp);
        db p=1;
        for(int j=0;j<to[i].size();j++) {
            int x=a[j].to;
            A[i][x]=p*a[j].p;
            p=p*(1-a[j].p);
        }
        A[i][i]=-1.0;
    }
    gauss();
    get_id();
    if (same()) {
        printf("%.7lf\n",A[1][n+1]/A[1][1]);
        return;
    }
    fo(i,1,n-1) id[i]=Id[i];
    solve();
}

int main() {
    freopen("graph.in","r",stdin);
    freopen("graph.out","w",stdout);
    scanf("%d%d",&n,&m);
    fo(i,1,m) {
        scanf("%d%d%d",&x,&y,&p);
        to[x].pb(y);v[x].pb(p);
    }
    solve();
    return 0;
}