From 8549749ccfdb985af8a223ba63d5cd1c06ee8e4a Mon Sep 17 00:00:00 2001 From: kilian Date: Wed, 20 May 2026 12:41:01 +0200 Subject: [PATCH] Dateien nach "ILP" hochladen --- ILP/caffeinesynthesis.py | 51 ++++++++++++++++++++-------------------- 1 file changed, 25 insertions(+), 26 deletions(-) diff --git a/ILP/caffeinesynthesis.py b/ILP/caffeinesynthesis.py index 87bfad5..c925a89 100644 --- a/ILP/caffeinesynthesis.py +++ b/ILP/caffeinesynthesis.py @@ -2,31 +2,31 @@ import gurobipy as gp from gurobipy import GRB, Model, quicksum HYPEREDGES = { - 1: ([], ["Xanthine"]) - 2: (["Xanthine"], ["p_{0,0}"]) - 3: (["Xanthine"], ["p_{0,1}"]) - 4: (["Xanthine"], ["p_{0,2}"]) - 5: (["p_{0,0}"], ["p_{0,3}"]) - 6: (["p_{0,0}"], ["p_{0,4}"]) - 7: (["p_{0,1}"], ["p_{0,3}"]) - 8: (["p_{0,1}"], ["p_{0,5}"]) - 9: (["p_{0,2}"], ["p_{0,4}"]) - 10: (["p_{0,2}"], ["p_{0,5}"]) - 11: (["p_{0,3}"], ["Caffeine"]) - 12: (["p_{0,4}"], ["Caffeine"]) - 13: (["p_{0,5}"], ["Caffeine"]) - 14: (["Caffeine"], []) + 1: ([], ['Xanthine']), + 2: (['Xanthine'], ['p_{0,0}']), + 3: (['Xanthine'], ['p_{0,1}']), + 4: (['Xanthine'], ['p_{0,2}']), + 5: (['p_{0,0}'], ['p_{0,3}']), + 6: (['p_{0,0}'], ['p_{0,4}']), + 7: (['p_{0,1}'], ['p_{0,3}']), + 8: (['p_{0,1}'], ['p_{0,5}']), + 9: (['p_{0,2}'], ['p_{0,4}']), + 10: (['p_{0,2}'], ['p_{0,5}']), + 11: (['p_{0,3}'], ['Caffeine']), + 12: (['p_{0,4}'], ['Caffeine']), + 13: (['p_{0,5}'], ['Caffeine']), + 14: (['Caffeine'], []) } #Similyrity of Nodes NMR to Measured NMR VERTICES = { - 1: ([0, 0], ["Xanthine"]) - 2: ([0, 0], ["p_{0,0}"]) - 3: ([0, 0], ["p_{0,1}"]) - 4: ([1, 0], ["p_{0,2}"]) - 5: ([0, 0], ["p_{0,3}"]) - 6: ([0, 0], ["p_{0,4}"]) - 7: ([0, 1], ["p_{0,5}"]) - 8: ([0, 0], ["Caffeine"]) + 1: ([0, 0], ['Xanthine']), + 2: ([0, 0], ['p_{0,0}']), + 3: ([0, 0], ['p_{0,1}']), + 4: ([1, 0], ['p_{0,2}']), + 5: ([0, 0], ['p_{0,3}']), + 6: ([0, 0], ['p_{0,4}']), + 7: ([0, 1], ['p_{0,5}']), + 8: ([0, 0], ['Caffeine']), } FIXED_FLOWS = { 1: 1, @@ -50,9 +50,8 @@ def build_model(name, hyperedges, vertices, excluded_support=None): model.addConstr(x[e_id] == value, name=f"fixed_flow_{e_id}") for e_id in hyperedges: - - model.addGenConstrIndicator(b[e_id], 0, x[e_id] == 0, name=f"unused_implies_zero_{e_id}") - model.addConstr(x[e_id] >= b[e_id], name=f"used_implies_positive_flow_{e_id}") + model.addGenConstrIndicator(b[e_id], 0, x[e_id] == 0, name=f"unused_implies_zero_{e_id}") + model.addConstr(x[e_id] >= b[e_id], name=f"used_implies_positive_flow_{e_id}") if excluded_support: model.addConstr(quicksum(b[e_id] for e_id in excluded_support) <= len(excluded_support) - 1, name="different_hyperedges",) @@ -77,7 +76,7 @@ def print_solution(title, flow_solution, binary_solution, hyperedges): print(f"Number of used hyperedges: {len(binary_solution)}") def main(): - model, x, b = build_model("HypergraphFlow", HYPEREDGES, NODES) + model, x, b = build_model("HypergraphFlow", HYPEREDGES, VERTICES) model.optimize() if model.status != GRB.Status.OPTIMAL: print("No optimal solution found for the first model.")