题目链接🔗:361. 观光奶牛 - AcWing题库
分析:
本题属于SPFA判断负环的综合应用 但是要从题面转换到SPFA判断环还是有很大难度的。
首先明确方法:由于答案上下界是确定的(0~1000)所以可以采用二分答案。设每次二分的答案为mid 则有 :
设一个环 S 的边权为 t1,t2,t3…点权为 f1,f2,f3…
若 mid<=ans,即存在一个环S使得 mid<=∑fi∑ti,变换一下:∑(mid∗ti−fi)<=0否则,则 mid>ans
每次 check 的时候,一条 u 指向 v,边权为 w 的边权变为:fu-mid*tv。我们只需检查这个图是否存在正环即可。
上代码~还有不懂请看注释 (注:正解确实就是这样 但可能由于是Py语言 会TLE三个点)
# 读入点权和边权
L,P = map(int,input().split())
lst = [0]
for i in range(L) : lst.append(int(input()))
Map = [[0 for i in range(L+1)]for j in range(L+1)]
for i in range(P) :
a,b,c = map(int,input().split())
Map[a][b] = c
def check(mid) : # spfa
flag = [1 for i in range(L+1)] # 队列里是否有这个元素
cnt = [0 for i in range(L+1)] # 计数
queue = [i for i in range(1,L+1)]
dis = [0 for i in range(L+1)]
while queue :
node = queue.pop(0)
flag[node] = 0
for i in range(1,L+1) :
if Map[node][i] :
new_dis = dis[node] + lst[node] - Map[node][i]*mid
if dis[i] < new_dis : # 注意这里是求最长路,因为是要判断是否存在正环
dis[i] = new_dis
cnt[i] = cnt[node] + 1
if cnt[i] >= L : return True # 存在正环
if not flag[i] : queue.append(i); flag[i] = 1 # 如果这个点不在队列里,加入队列
return False
left,right,eps = 0,1010,10**(-4) # 二分答案
while right - left > eps :
mid = (left+right)/2
if check(mid) : left = mid
else : right = mid
print('%.2f'%right)