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;
}