这题在比赛的时候没做出来实在是不应该。本来也是用对数处理的,但是后面写的七零八落的。
这题也可以直接求,记录每个点最少消耗的流量,剩余流量为初始流量减去这个值,逐步向下迭代。
当然这里可以去求他的剩余值,根据公式 最后的流量 L = I * (1-p1) * (1-p2) * ... * (1-pn) 我们对两边同时取对数的话,那么我们就将乘法化成了加法,并且直接求一个最长路就可以了。最后再拿总流量减去最大的剩余量即可。
代码如下:
#include <cstdio> #include <cstring> #include <cstdlib> #include <queue> #include <algorithm> #include <cmath> #define eps 1e-6 using namespace std; int head[50005], N, idx, S, T, P, q[10000005], visit[50005]; double dis[50005]; struct Edge { int v, next; double fee; }e[2500005]; void insert(int x, int y, int fee) { ++idx; e[idx].v = y, e[idx].fee = log(1 - fee/100.); e[idx].next = head[x], head[x] = idx; } void spfa() { int front = 0, tail = 1, pos; q[tail] = S; fill(dis+1, dis+N+1, 0); dis[S] = log(1.*P); while (front != tail) { pos = q[++front]; visit[pos] = 0; for (int u = head[pos]; u != -1; u = e[u].next) { if (dis[pos] + e[u].fee - dis[e[u].v] > eps) { dis[e[u].v] = dis[pos] + e[u].fee; if (!visit[e[u].v]) { // spfa中防止多次入队 q[++tail] = e[u].v; visit[e[u].v] = 1; } } } } if (dis[T] == 0) { puts("IMPOSSIBLE!"); } else { printf("%.2lf\n", P - exp(dis[T])); } } int main() { int M, y, fee; while (scanf("%d", &N) == 1) { memset(head, -1, sizeof (head)); memset(visit, 0, sizeof (visit)); idx = -1; for (int i = 1; i <= N; ++i) { scanf("%d", &M); for (int j = 1; j <= M; ++j) { scanf("%d %d", &y, &fee); insert(i, y, fee); } } scanf("%d %d %d", &S, &T, &P); spfa(); } return 0; }