牛客 跑路ing(java解题）

题目描述

vigoss18 辞职成功终于逃出了公司，但是没过太久，公司就发现vigoss18 的所作所为，于是派人来把他抓 回去。 vigoss18
必须一直跑路，躲避公司的围捕。可以抽象的看成一个有向图，图中可能存在重边和自环。 刚开始他站在位置1，每单位时间vigoss18
必须从目前站的位置，等概率选择一条边然后移动到对应的节
点上去或者不动（如果当前节点有t条边，则有1/(t+1)的概率选择一条边移动或者原地不动），可以认为每次需 要花费1 单位时间。
他就这样一直跑一直跑，过了很长很长的时间… 公司把你派出来寻找vigoss18，如果能抓到他，你将能升官发财赢取白富美走向人生巅峰。
但是你精力有限，不是太走的开身，所以写了一个程序，来计算vigoss18 在每个位置的概率，可以认为过 了很长时间以后，vigoss18
在每个位置的概率是收敛的。所以你需要告诉上司，他最可能在哪个位置(概率 最大的那个位置)。
你的上司并不想知道过程，他只想知道结果，所以你只需要告诉他这个概率最大是多少即可。

输入描述:

多组输入，保证绝大部分为小数据。 每组输入第一行n m(1<=n<=100,1<=m<=10000)，表示n个点m条有向边。
接下来m行，每行u v(1<=u,v<=n)，表示有一条有向边从u连向v

输出描述:

算出vigoss18在所有位置的概率，并输出其中的最大值即可。 你的答案与标准答案的误差应保持在1e-6以内。

输入

3 3 1 2 2 3 3 1

输出

0.333333333

思路
第一步 定义dp状态

发现dp[i]根本不够，无法做出此题 一定需要再来个数组存储当前时间的概率 所以定义dp[ i ][ j ]定义为当前 i
单位时间在 j节点的概率

第二步 定义状态方程

显然此题需要存储边，又因为是有重边和自环，我们直接定义一个二维数组存储这些边，用于根据当前节点算出概率1/(t+1)

ArrayList<Integer>[] t = new ArrayList[n+1];
			for(int i=1;i<=n;i++){
				t[i] = new ArrayList<>();
			}
			for(int i=0;i<m;i++){
				int u = sc.nextInt();
				int v = sc.nextInt();
				t[u].add(v);
			}

初始化dp[ 0 ] [ 1 ] 概率为1（因为从 1 开始） 因为概率会趋于收敛 所以 i 的最大值 200多可以，不行就加
然后就是循环遍历更新dp状态 dp方程为 当前节点 i 单位时间概率 + 当前节点 i - 1单位时间概率 dp[ i ][ j ] +=
dp[ i - 1 ][ j ] * ( 1.0 / ( g[ j ].size() + 1)); 当前节点连接的另一个节点 i
单位时间概率 + 当前节点 i - 1单位时间概率 dp[ i ][ to ] += dp[ i - 1 ] [ j ] * ( 1.0
/ ( num + 1));

double[][] dp = new double[201][n+1];
			dp[0][1] = 1;
			for(int i=1;i<=200;i++){
				for(int j=1;j<=n;j++){
					int num = t[j].size();
					dp[i][j] += dp[i-1][j]*(1.0/(num+1));
					for(int to : t[j]){
						dp[i][to] += dp[i-1][j]*(1.0/(num+1));
					}
				}
			}

代码：

import java.util.ArrayList;
import java.util.Scanner;

public class Main {
	public static void main(String[] args){
		Scanner sc = new Scanner(System.in);
		while(sc.hasNext()){
			int n = sc.nextInt();
			int m = sc.nextInt();
			ArrayList<Integer>[] t = new ArrayList[n+1];
			for(int i = 1; i <= n; i ++){
				t[i] = new ArrayList<>();
			}
			for(int i = 0; i < m; i ++){
				int u = sc.nextInt();
				int v = sc.nextInt();
				t[u].add(v);
			}
			double[][] dp = new double[201][n+1];
			dp[0][1] = 1;
			for(int i = 1; i <= 200; i ++){
				for(int j = 1; j <= n; j ++){
					int num = t[j].size();
					dp[i][j] += dp[i-1][j]*(1.0/(num+1));
					for(int to : t[j]){
						dp[i][to] += dp[i-1][j]*(1.0/(num+1));
					}
				}
			}
			double res = 0.0;
			for(int i = 1; i <= n; i ++){
				res = Math.max(res,dp[200][i]);
			}
			System.out.printf("%.9f",res);
			System.out.println();
		}
	}
}