def run():
V = 5
names = ['USD', 'EUR', 'GBP', 'CHF', 'CAD']
rates = [[1, 0.741, 0.657, 1.061, 1.005],
[1.349, 1, 0.888, 1.433, 1.366],
[1.521, 1.126, 1, 1.614, 1.538],
[0.942, 0.698, 0.619, 1, 0.953],
[0.995, 0.732, 0.650, 1.049, 1]]
g = EdgeWeightedDigraph(V)
for v in range(V):
for w in range(V):
e = DirectedEdge(v, w, -math.log(rates[v][w]))
g.add_edge(e)
# print(g)
spt = BellmanFordSP(g, 0)
# print(spt)
print(spt.negative_cycle().queue)
print(spt.has_negative_cycle())
if spt.has_negative_cycle():
print('hit')
stake = 1000.0
for e in spt.negative_cycle().queue:
print(e)
print(f'{stake} {names[e.tail()]}')
stake *= math.exp(-e.weight())
print(f'= {stake}\n {names[e.head()]}')
else:
print('No arbitrage opportunity')
def run():
with open("../resources/jobsPC.txt", ) as f:
values = "".join(f.readlines()).splitlines()
print(values)
n = int(values[0])
source, sink = 2 * n, 2 * n + 1
durations = list()
for x in values[1:]:
durations.append(float(x.split(' ')[0]))
precedent = list()
print(values[1:])
for prec in values[1:]:
if prec.split(' ')[-1] == '0':
precedent.append([])
else:
precedent.append("".join(prec.split(' ')[-1]).split(' '))
print(n, durations, precedent)
g = EdgeWeightedDigraph(2 * n + 2)
for i in range(n):
g.add_edge(DirectedEdge(source, i, 0.0))
g.add_edge(DirectedEdge(i + n, sink, 0.0))
g.add_edge(DirectedEdge(i, i + n, durations[i]))
m = len(precedent[i])
for j in range(m):
p = int(precedent[i][j])
g.add_edge(DirectedEdge(n + i, p, 0.0))
print(g)
lp = AcyclicLP(g, source)
print('job start finish')
print('----------------')
for i in range(n):
print(f'{i}, {lp.dist_to(i)}, {lp.dist_to(i+n)}')
print(f'finish time {lp.dist_to(sink)}')
def run_nc():
with open("../resources/tinyEWDnc.txt", ) as f:
values = "".join(f.readlines()).splitlines()
V, E = int(values[0]), int(values[1])
g = EdgeWeightedDigraph(V)
for line in values[2:]:
vertices = line.split(' ')
v, w, weight = int(vertices[0]), int(vertices[1]), float(vertices[2])
e = DirectedEdge(v, w, weight)
g.add_edge(e)
s = 0
sp = BellmanFordSP(g, s)
if sp.has_negative_cycle():
print('negative cycle')
q = sp.negative_cycle()
while not q.empty():
print(q.get())
else:
for t in range(g.get_V()):
if sp.has_path_to(t):
print(f'{s} to {t} ({sp.dist_to(t)}) ', end="")
q = sp.path_to(t)
while not q.empty():
print(q.get())
print()
else:
print(f'{s} to {t} no path\n')
def main():
g = EdgeWeightedDigraph(8)
with open("../resources/tinyEWD.txt", ) as f:
for line in f.readlines():
vertices = " ".join(line.splitlines()).split(' ')
if len(vertices) < 3:
continue
else:
v, w, weight = int(vertices[0]), int(vertices[1]), float(
vertices[2])
e = DirectedEdge(v, w, weight)
g.add_edge(e)
print(g)
finder = EdgeWeightedDirectedCycle(g)
print(finder)
if finder.has_cycle():
print('cycle')
while not finder.cycle().empty():
e = finder.cycle().get()
print(f'{e} ')
print()
else:
print('No directed cycle')
def __init__(self, g):
self._marked = [False for _ in range(g.get_V())]
self._on_stack = [False for _ in range(g.get_V())]
self._edge_to = [DirectedEdge() for _ in range(g.get_V())]
self._cycle = None
for v in range(g.get_V()):
if not self._marked[v]:
self.__dfs(g, v)
assert self.__check()
def main():
with open("../resources/tinyEWD.txt", ) as f:
values = "".join(f.readlines()).splitlines()
V, E = int(values[0]), int(values[1])
g = EdgeWeightedDigraph(V)
for line in values[2:]:
vertices = line.split(' ')
v, w, weight = int(vertices[0]), int(vertices[1]), float(
vertices[2])
e = DirectedEdge(v, w, weight)
g.add_edge(e)
dijkstra = DijkstraAllPairsSP(g)
print(dijkstra)
def main():
g = EdgeWeightedDigraph(8)
with open("../resources/tinyEWD.txt", ) as f:
for line in f.readlines():
vertices = " ".join(line.splitlines()).split(' ')
if len(vertices) < 3:
continue
else:
v, w, weight = int(vertices[0]), int(vertices[1]), float(
vertices[2])
e = DirectedEdge(v, w, weight)
g.add_edge(e)
print(g)
print(str(g))
def main():
with open("../resources/tinyEWDAG.txt", ) as f:
values = "".join(f.readlines()).splitlines()
V, E = int(values[0]), int(values[1])
g = EdgeWeightedDigraph(V)
for line in values[2:]:
vertices = line.split(' ')
v, w, weight = int(vertices[0]), int(vertices[1]), float(vertices[2])
e = DirectedEdge(v, w, weight)
g.add_edge(e)
s = 5
asp = AcyclicSP(g, s)
for t in range(g.get_V()):
if asp.has_path_to(t):
print(f'{s} to {t} ({asp.dist_to(t)}) ', end="")
q = asp.path_to(t)
while not q.empty():
print(q.get())
print()
else:
print(f'{s} to {t} no path\n')