GCPC2017E Perpetuum Mobile题解

GCPC2017E Perpetuum Mobile题解

一张图，图上有n$(2\le n\le 800)$个转换器，m$(0\le m\le 4000)$个转换关系。如果转换器a,b之间有一个权值为c$(0<c\le 5.0，精度最高为4位小数)$的转换关系，表明向a注入1的能量，会从b产出c的能量(翻了c倍)。那么显然，如果该图上有一个权值之积大于1的环，那么这张图违背了能量守恒，是不合法的
题目保证给定数据中环权值之积$p\le 0.9$或$p\ge 1.1$

思路：
图中最多有800个点，4000条边，想要暴力对每一个环算权值的积不太现实。更不用说权值的精度可以到小数点后4位，直接乘对精度的要求也很大（不过既然权值之积有范围，可能精度的问题也不是很大）。

分析：我们需要判断对于一个环，其权值的积是否大于1，即：
$$\begin{array}{r}\prod _{i=0}^n{c}_i>1\end{array}$$
两侧取对数
$$\begin{array}{r}\sum _{i=0}^{n}ln{c}_i>0\end{array}$$

也就是说，可以通过取对数将权值的积转换为权值对数的和。也即我们可以对权值取对数在建图，判断图上是否有正权环。

AC代码：SPFA判正环，注意判负环用的是最短路，判正环就要用最长路了。然后就是边的权值dis数组这些平常开longlong的变量别忘了改成double

#include <bits/stdc++.h>
using namespace std;

const int maxn = 810;
const int maxm = 4010;
int n, m;

struct edge {
	int next, to;
	double w;
}e[maxm << 1];
int tot;
int head[maxn];

void add(int u, int v, double w)
{
	e[++tot].w = w;
	e[tot].to = v;
	e[tot].next = head[u];
	head[u] = tot;
}

int vis[maxn];
double dis[maxn];
queue<int> q;

int cnt[maxn];

bool spfa(int s)
{
	for (int i = 0; i <= maxn; i++) {
		dis[i] = INT64_MIN;
	}
	memset(vis, 0, sizeof(vis));
	while (!q.empty())
		q.pop();

	q.push(s);
	dis[s] = 0;
	vis[s] = 1;
	cnt[s] = 1;
	while (!q.empty()) {
		int u = q.front();
		q.pop();
		vis[u] = 0;
		for (int i = head[u]; i; i = e[i].next) {
			int v = e[i].to;
			if (dis[v] < dis[u] + e[i].w) {
				dis[v] = dis[u] + e[i].w;
				cnt[v] = cnt[u] + 1;
				if (cnt[v] > n) {
					return true;
				}
				if (!vis[v]) {
					q.push(v);
					vis[v] = 1;
				}
			}
		}
	}
	return false;
}

int main()
{
	scanf("%d%d", &n, &m);
	for (int i = 0; i < m; i++)
	{
		int u, v; double w;
		scanf("%d%d%lf", &u, &v, &w);
		w = log(w);
		add(u, v, w);
	}

	int hasPosRing = spfa(1);
	if (hasPosRing)
		printf("inadmissible");
	else
		printf("admissible");

	return 0;
}