benchmarks.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {
    "pycharm": {
     "name": "#%%\n"
    }
   },
   "outputs": [],
   "source": [
    "import sys\n",
    "sys.setrecursionlimit(20000)\n",
    "\n",
    "path1 = \"tomate-virus.txt\"\n",
    "file_1 = open(path1, \"r\")\n",
    "stripped_1 = ''.join([line.rstrip('\\n') for line in file_1])\n",
    "\n",
    "path2 = \"botryosphaeria_dothidea_chrysovirus_1_uid382570.txt\"\n",
    "file_2 = open(path2, \"r\")\n",
    "stripped_2 = ''.join([line.rstrip('\\n') for line in file_2])\n",
    "\n",
    "path3 = \"Tomato ringspot virus polyprotein.txt\"\n",
    "file_3 = open(path3, \"r\")\n",
    "stripped_3 = ''.join([line.rstrip('\\n') for line in file_3])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "pycharm": {
     "name": "#%% md\n"
    }
   },
   "source": [
    "# Run LCS vs LCS_Optimized_2_Lines vs LCS_Optimized_Memory_d_lines"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {
    "pycharm": {
     "name": "#%%\n"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "AGGAGGGGCTTTTTTTGCCGGTAATCGTAAAGTGGTTTCGCACTTCGTCCAATTATTCGCGCCTGAGTTCTAGATTAAACACTGACCATGTGAGCGATATAAGAAAGACCCTGCCAATGCTCATGTCGAAATGCCGAAGGCGCCCCGGCCCGATGCGGAACCCTAACAGCCGCTGCCTCCCGCTCGGCAGGTATCGGCCCACTCAACAACCTCTGAGGTACAGGCCATTTAGGAAGCCGGATCTTGCACGCTGACCCGAGCTAGAGGTGTGAAGGGCCGGTAAGGTCCAGCTTATAACAGACGAACATTCACGCGCTGTGTTTTCGCACCGTGTGCTACCCACCTTGGTCGTTTTGTAAGCTGTGACACTAGGTTGGTGGTAGCCGTGAACCCCGAAAGTGGTGGGGGAGCCAGCCATGGAAGCGGATTTTGGCGTGCACATGTTCGCAACAACTAGCACTGCACTACGCCCGTGACGACCTCAGCACGTGGAAGTAAGGTACAGTATGGGCAACGAGCAGGTTTGTCAGGCGATCTGGAAGCAACTCATGTGGTCTACAACAAGGAGGCGGAGTATAGAATCCACTGCTTTGTTCTGTACATGGCTGCCCGCCAACCCTTATATGCACTCGAGACCGGACAAGATTCGGGTTTCAATTTAATTGATTTATATGCCAGTACAGAATTACATACGTTGTCGGCGAATTGCAGTCAATTCTATTAATGGAAATACCCTACTGGTTGTAAACAAACTAATCCTAAATTTATTTTTGTCGTAGGCTTACTGTGCGTCCAGACTTGGAATTCGTGGCTTGTGAGTTATTTCCTAGGTGGATTAACCGTTTGTTGAACACCTGGGTTCGTGTACTGGTTCTCACGTTTTATCTTGAATTAGGGGATTTGTTATCCCCAAAAAGCCATTCTCGCCTGATTGCATGATAGTTTGGGTGAATCCACCCAGTGCGCCTATGGGAGGTAGATGGAGAACCGACTCTAGATCGCGCCTCCTCCGATGGCCCCCGGCGCGTGTCTTCTATGGGGGGGAGCTGAGGGTGATGAAGACGCTTCTGAGGGTCGTCTGAGTGATGCGTTTCCCTTTGTTGGAAACAAGCTGGAACCCGACCCGGATTGAAAGCCTTGAGGGATCCCCTGTTTGAACTCTTCCTACTTGCAATTTGACTGCAGTCTTGGACCCATCATCTTTCCAGGCAGCAGATATGGGGCGCGGCGCAATGGGCTCCACTGATTGATAGCCACCAGAAACCATGCCATCAGATAGGTATGTGGGCCTAACCACGGCCCCATTGTCTTTGTTCGAATTCACCTCGATGAGATCTAATCTTCCCGGCCGCGAGCGCACCTCTCCAAATATTATCGCACCTATCCTCTGGAACAGTGGAGAGGGGGAAACGAGGCCATGGTCCGGAGCCTTAGAAAATCCTCGCATTGTTGAGTTAGATAACAAAAAAGTCCGGGTTTTTCTTGTTTTGAGGGTGTTGCAGCCTGCGAAGCCGGCATATCGGGCAATGGCAACCTTTACGATCGTTGATCTGGAAAACTCTGATAGAAGTCTCCGCCCCAATAGCTTTGAGTCGAGTGCCTTTAGTCCCTGAATGTTCGGGAAATGTGCTGATCGGTGGGATACCAAGTCGGACTTTGTTGGTGCATTGGACCTCGAATTGGATAAGACATGGAGATAGGCTCCCCATTCTCTGAAGCTCTTGCCACGAATGATAATGAGTGGGGTTTAAGGTGAGAAGGGGTGCTCGTCTGGTGAAGAGCAATGAGATAAGTAAGGAAAATTGGAAAGTTTTAAACGCTTGGGGGTGGTTGGTATATGGCGCCCCCCGTTGGGCCTGAAAACTGTGACAAATTGGGGTATGGGGTACATTATAAGAGAAGTCTCAATAGAACTTCGACGAACCGCACACGGGCCCACCGTTAAATTT\n",
      "GAAAAAAGAGAAAGGAATTTCTTCTTGCCGCGAAGTAATGTGAGTGACCTTGCTTCACTCGGTACGCTGAGGCAAAATTAATTCTGCAGGCAAGCTCTGACGATTACCTCACAAGATTCAACCCAACTGACCATCCCGCCTGACCACGACGGATACCAAGAACAGCCCCCTGCCAGAAGTACCTGCAATGCACCACCTAGGTCCGTACCGGTAGCAAGCAAAACGCGGAACGGCCTAACACCCTCCGCCTCGCCGGGTGCCGAAGGTCGATCGTTTGCCGTCTGCGGGCAAGTCCCCCTGTCCAGGCGGGCCGGATTTAATGATATTGTCGCGGACTGTCGCTACACATACCCGAGTAGACGGATGCTCCAATGCGCGGCCGGCCATCCACGGGAAAGCCCATTTCAAGCAGAACGAACTCATGCCACCGCCTCTCGCTTCTTCGCCACCGTTCGCCTGACCAACTAGGTTTTCTTCCCGGAGAGGCCGGATAGGCTGGGGCAGCCGTACTCATCCCGGTGAAGTGGTGGCGGCACCAGGCTGGTGGCGACGTTGTGTGCGCGTGCCACGATGCGGTGTCCGACAACGTGGCGAACTGGCCTAGGTCCCCCCTCGCAGTGGAAGTCGGCTTACAACTGTGACGCGTAGGGGGTTTGCACGCGCGCGAGGCTGCAGTAACGACCAACTCGCGCTACTGGCCATGGCCTTTAACCACCGAAGAGGCGGACTACTCAACCGAGTCTCGCCTCTTGCGTTCGGTCCCAGGCCCGGCTCAGGGTGGCTGCCCCCCTGCAGCGCCTGGGCGCGGTGCACCGGGACCACCAGGCCAACCACCCGATGGGGCGAGGTTGTGCCGCATCTTCAGCGGGTCGAGCCGCTCTGACATACTGCCCGGTGGGACGGGATGCGCGACGACCGTGTCCAGGCAACCAATTCGATGTCAAATCTCGGCTGGGGCGAATGGGGGTGCGCGGCGGACCCTGCTGTGGGGTTCAACCTGGCGCACGGGATTCGGTACGAACTTGCATCCATCTTTTTCGGGGATGAGTGCTTACGTGTGCGTCCTGAGCACAAAGAGGCTGAACCCTGCGCTCGGCCAGTGACTTGACGCGTATGGAAGGAGATCATAACAACGTCCGGTGGCTTGAACACCTTGGGGTGGCTCTGCTCCAAGTTTGGCTACTCACGTACTGACGCCTTCGGGGCCGAGGGGGGAAATCGAGGCTGAGAAGCGGACGCGGCTGCAAAAGCCAAGTCTTCTCAGGGCCTTCACATCTGATTCACTCGGGGCCCCAAACCGCAATGAGAGTGTGGCCTCGGGAAGTCTGAGGTGTCCGCGTGAAACCTCAACTTCCTGGAACGGCCGGTGGGAAACCTGGCCAGCGCTCTCCCCGGCGACGTACTTGCGGTGGGGCCTCACAGCTGGGGATGACCACCAGACAGCGGAATGACTGCTGGTGCTGAAGCGGGTGCGCTGCAAGTCAGGTGTTCGAGTTCCCTTCGTGGATATGGAAGCGCGCTGAAGATTACGTGCCAACAGCTCACGAATACTCTTCAAGAAAGACCTTCAATGAATCGCCATGTTCAGGACTCTCTTCACTACATAAGGGAGCCACACTTGACCTGATCGAAAAGTGGGCCTATACGAACTCGCCTCGGTGTTCCCGGCGCGCGCGAGGATCACTGGTTTTGGCGGAGTCGCGGATTGGTGGTCCATTGGCTGGATACGCCATTCAAGAAACGCAGCATCAGCCTGCAGCAGAATATAGGTCAGGCAATCTGCAAACCAACCGGTTGGGATTGCCAATCAATTGGTGGGAGCGGGTTTCGGCGTCGGTGAGGGATTGATATTCCACATGGCCTGCTCTCTGGATGGCGAACAATCGCACAACAGATCTTGAACGGGTTACCTCCGACTGGGTCAACGTGGGTGTTGTCGGTTTTCACCCCTCCTCAACTACCAGCAGCGCTACGCGCGGCGGAGGGCGACCTTTCGAGACGGCAGTCACCCTGGGTGGCAAATAGGGGTTCGGTCAAGACATCACCTGATAGGTTAAGGACCGAGAGTTCTTTCGATGCGGTGGAAGGGGCTCGGCAGATGCCGAATGCTGGGTCACGTTATCCCAATTCCAACAGATTTTACCGAATTCGGCCGTGGCCAAACTTGCCGTTTAAAGTGCTCGCAAAGCGGCGCGCAGAGCGCACTCATGGGAAGAGGGACGCTTTCAAAGTGGGCACTGATGGGGAGCTGGTTCGCGGCCTTTGCATACGGGAAACTGGCGGCGGGTTACCATTTACGGCGACGTTCGATGATCACCACCGACGGAAACTCGAACCGGGTTGCGACAGATGCTAGTAGGCTAGGTTGATGTCAGGATGCGCAATTCCCCCCCGAAGCTACTTGCGCACGAAGTCTGCCACTGTGGGTGGGGCTTTCACAGTTCAGTCGCGAGAACGCGGCCCGGGGTGCTCGGGGGGGGCTGCAAGACCCAAGCTGGTACCATCAAGCGTATCAGAGCTGAAGAATTCCGGCAACCACAAGTGCGCTATTCCAAGACACTAATACGTGGCGAGACTGGTTACGCGCTCACGGCAGCGGGCGCTACCGGCTGTCCAGGCCGGAAGCACCCTTACCGTAGATTTGGCCAGCGATTCGCGACTGGAGCAGATCCCGCGTCATCTACTTGGGGAGATCGGATTATGGCTGCCAACACTTCTAAGGAAGACTAATCGCCAAACGACAAGGCAACGCCGAAACCGCCCCCGGGCTAAAATTGT\n",
      "59\n",
      "[971.71875, 489.1484375] [7.058798789978027, 12.5]\n"
     ]
    }
   ],
   "source": [
    "import math\n",
    "import time\n",
    "\n",
    "from LCS.LCS_Standard import print_lcs, lcs_length\n",
    "from LCS.LCS_Optimized_2_Lines import lcs_length as lcs_length_two_lines  # lcs optimized\n",
    "from LCS.LCS_Optimized_Memory_d_Lines import lcs_length as lcs_length_d_lines, print_lcs_reconstruct\n",
    "from memory_profiler import memory_usage\n",
    "\n",
    "## Benchmark LCS\n",
    "# time\n",
    "\n",
    "# tomato vs chrysovirus\n",
    "bench_start = time.time()\n",
    "c, b = lcs_length(stripped_1, stripped_2)\n",
    "#print_lcs(b, stripped_1, len(stripped_1), len(stripped_2))\n",
    "bench_end = time.time()\n",
    "# tomato virus vs dna\n",
    "\n",
    "# botryosphaeria_dothidea_chrysovirus_1_uid382570 vs Tomato ringspot virus polyprotein\n",
    "bench_start2 = time.time()\n",
    "c, b = lcs_length(stripped_2, stripped_3)\n",
    "bench_end2 = time.time()\n",
    "\n",
    "\n",
    "# memory\n",
    "mem = memory_usage((lcs_length, (stripped_1, stripped_2)))\n",
    "\n",
    "# results\n",
    "lcs_naive_benchmark = bench_end - bench_start\n",
    "lcs_naive_benchmark2 = bench_end2 - bench_start2\n",
    "lcs_naive_memory_usage = max(mem)\n",
    "\n",
    "##############################################\n",
    "# Benchmark LCS 2 Lines\n",
    "#time\n",
    "bench_start = time.time()\n",
    "c = lcs_length_two_lines(stripped_1, stripped_2)\n",
    "bench_end = time.time()\n",
    "\n",
    "# botryosphaeria_dothidea_chrysovirus_1_uid382570 vs Tomato ringspot virus polyprotein\n",
    "bench_start2 = time.time()\n",
    "c = lcs_length_two_lines(stripped_2, stripped_3)\n",
    "bench_end2 = time.time()\n",
    "\n",
    "# memory\n",
    "mem = memory_usage((lcs_length_two_lines, (stripped_1, stripped_2)))\n",
    "\n",
    "#results\n",
    "lcs_optimized_two_lines_benchmark = bench_end - bench_start\n",
    "lcs_optimized_two_lines_benchmark2 = bench_end2 - bench_start2\n",
    "lcs_optimized_two_lines_memory_usage = max(mem)\n",
    "\n",
    "##################\n",
    "# Benchmark LCS with d lines\n",
    "#time\n",
    "bench_start = time.time()\n",
    "d = round(math.sqrt(len(stripped_2) + 1)) # d=runden(wurzel(n+1))\n",
    "lenLCS, cb, bb = lcs_length_d_lines(stripped_1, stripped_2, d)\n",
    "print_lcs_reconstruct(cb, bb, stripped_1, stripped_2, d)\n",
    "bench_end = time.time()\n",
    "\n",
    "# botryosphaeria_dothidea_chrysovirus_1_uid382570 vs Tomato ringspot virus polyprotein\n",
    "bench_start2 = time.time()\n",
    "lenLCS, cb, bb = lcs_length_d_lines(stripped_2, stripped_3, d)\n",
    "print_lcs_reconstruct(cb, bb, stripped_2, stripped_3, d)\n",
    "bench_end2 = time.time()\n",
    "\n",
    "#memory\n",
    "mem = memory_usage((lcs_length_d_lines, (stripped_1, stripped_2, d)))\n",
    "print(d)\n",
    "\n",
    "# results\n",
    "lcs_optimized_d_lines_benchmark = 12.5 # bench_end - bench_start\n",
    "lcs_optimized_d_lines_benchmark2 = bench_end2 - bench_start2\n",
    "lcs_optimized_d_lines_memory_usage = max(mem)\n",
    "\n",
    "\n",
    "#results_time = [lcs_naive_benchmark, lcs_optimized_two_lines_benchmark, lcs_optimized_d_lines_benchmark]\n",
    "#results2_time = [lcs_naive_benchmark2, lcs_optimized_two_lines_benchmark2, lcs_optimized_d_lines_benchmark2]#\n",
    "results_time = [lcs_naive_benchmark, lcs_optimized_d_lines_benchmark]\n",
    "results2_time = [lcs_naive_benchmark2, lcs_optimized_d_lines_benchmark2]#\n",
    "results_memory = [lcs_naive_memory_usage, lcs_optimized_d_lines_memory_usage]\n",
    "\n",
    "print(results_memory, results_time)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {
    "pycharm": {
     "name": "#%%\n"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "<Figure size 432x288 with 1 Axes>",
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEICAYAAABYoZ8gAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAA1lUlEQVR4nO3deZxN9f/A8dfbGMYW0uSrpMFXiJkGgyQ1qIQKlS1fX7R/W7Tpi1bt+vItP21CaJEt8ZVKWoxSxIzGWMuSsmaL7Fnevz/Omeka9945s5wZY97Px+M+7j37+5x77vue+7mf8/mIqmKMMaboKFbQARhjjMlflviNMaaIscRvjDFFjCV+Y4wpYizxG2NMEWOJ3xhjihhL/CGISCcR2SAi+0SkQQ7X0UJEfszr2Io6EVER+XuIaT1EZHYO19tbROblLjoTjIisF5ErCjqOvJab860gFfrE7+MJNRS4R1XLquoPOVmBqn6jqrXTh0/Xk/9UoqrjVfWqgo4jr4hIjPtFVzzE9AtEZIqI7BCRPSKSJiIPikiEO/0WEVklIntF5DcR+UREyoVYV5KI3Orn/pxKsvtFH+y9KKznW6FP/D46H1he0EGcrkIlMuOdiNQEvgc2ALGqWh7oDCQA5UTkcuB5oLuqlgPqApMKKl5zClHVQv0A1gNXBBlfEZgJbAd+d19XDbUcMAh4DygJ7AMU2A+sBbq649Ifh4Ekd7mSOL8OfgV+A0YApdxpicBG9/W7wHHgoLuOf7vjLwa+A3YDS4DEgJiSgGeAb4G9wGzgrBDHYSVwTcBwcXffG7rD1+F8ke1211s307HoB6QBe3CSQ1SQbZR0l68fMC7a3aez3eFrgFR3vu+AuEzb6e9u57AbY39gk7t/PwKt3XkjgEfc478XSAHOc6cpcCew2t3Oa4C403oD8wK2WQf4HNjlrr9LwLRKwAzgD2Che6y9LjvO3e7HbnzfAzXzYtlMxzzG3d/iQaa9B3wc5rPRD5iejc9SEnBr4LkL/BvYBmwBOgLtgJ/c/Xok0+fnA/fc2QssBi4K9nnDueAc4L63O4HJwJmZ9rcPzhfa7+573dg9b3YDr2aK+2ac8/934DPg/IBpQc8VnC/BQ8AxnM/jbnf+9sAP7jmxARgUsK5f3fWl54FmnHy+XQIswvkcLQIu8fp5Jo9ygaf3OreJt6AfhE78lYAbgNJAOWBK4Icg83LuiftephPm70HWe4Z7kt3hDr+MkzzOdLfzEfBC4IcnzDbPdU/8du6H4Up3ODrgzV4LXACUcocHhzgOTwDjA4bbAyvd1xfgfIldCUTifJjXACUC4loInOPux0rgzhDbGQM8FzB8NzDLfd0AJ0k0xUncvdx1lwzYTipwnrs/tXE+XOcEfOhruq8fBpa68whwEVAp4L2ZCVQAquF8wV3tTuuN+0EEyrjr74PzJdMA2AFc6E6fiJN0ygD1cb6AvC47zn2vmrjTxwMTc7tskOMdQ+jEvxXoE+az0QLnS/kpoHn6+xBm/iROTPxHcc6rSOA29zi/j3Oe13PXXT3g83MEuNGdvx/wMxCZ+dwH7gMWAFVxLibeBCZk2t8RQBRwFU6Cng6cjfOZ2QZc7s7fAedcrusey8eA7zJ9jrM8VwLmTwRicT6PcTgXcx1DvReceL6difPl09ONpbs7XCng+Ab9PJOHucBT3vQzKefHgxCJP8h88cDvoZbDQ+J335CZwBvusOAk1MArvWbAzwEnUbjE3x94N9M2PgN6BbzZjwVMuws3yQbZv7/jXAmUdofHA0+4rx8HJmfaj024VxRuXP8ImP4fYESI7VwBrA0Y/hb4p/v6DeCZTPP/yF8f0vXAzZli3uauMzLIch1CxKDApQHDk4EB7uve/PVB7Ap8k2nZN4Encb6YjgB1AqY972VZ9/U4YHTAtHbAqtwuG2RfYwid+I/gJrEw531bnIuR3ThXqS8BESHmTeLExH8wfV6cZK9A04D5U/grKQ4CFmQ6x7YALTKf+zgXFq0D5q3i7kvxgP09N2D6TqBrwPBU4H739afALZm2ewD3qt/ruRLm+A0DXg71XnDi+dYTWJhp+flA76w+z+RhLvDyOG3LWUWkNM7V+NU4xT7glHtGqOqxHK72OZwPQF93OBrnF0WKiGRsGiepeHE+0FlErg0YFwnMCRjeGvD6AFA22IpUdY2IrASuFZGPcIp20msjnQP8EjDvcRHZgHOVEWo754SIeQ5QWkSa4lwNxQPTAvanl4jcGzB/iUzr2pAp5vtxkkY9EfkMeFBVN+P8KlgbIoZg8QY7LucDTUVkd8C44jjFbtHu6w0B034JeB1u2axiyM2y2bETJ2mGpKqfAp+KSDGgJc4v3x9xvoiyXH/AZ+Wg+/xbwPSDnBh34Ht7XEQ2Evw8Oh+YJiLHA8YdAyoHDGfeTqjtng/8n4j8N2C64Jzb6e+n52PtnteDcX4BlsD5RTIl1PyZnPA5c/1C+M9Z4H7kSS7w4rRN/MBDOMUETVV1q4jE45TdpWfo/ThJO93fwq1MRLrh/HRrrKpH3NE7cE7Ceqq6yUNMmml4A863/G0elvVightjMWCFqq5xx2/G+fkKgDjfUufhXPVni6oeE5HJ7nZ+A2aq6l538gacYqDnwq0i0/reB94XkTNwktGLOFdOG4CawLLsxhhgAzBXVa/MPMGt9XIU5zisckdX87Jsbrabx77AKc4cm9WMqnoc+FJEvsJJan44L/2F+0VTFefcy2wDzi+/bzNPEJGYbG4z/Zwbn83l4OTPIzhFWa8CbVX1kIgMA84KM3+gzTgJPFA1YJaHWPI6F4R1utTqiRSRqIBHcZwr84PAbhE5E+fnfaBUoJuIRIpIAk7ZZFBuPf5XcH7Wbk8f736YRgEvi8jZ7rznikibEKv6DagRMPwezhV6GxGJcGNPFJGq2dn5ABNxykT/hXMCp5sMtBeR1iISifOleBjnj6SceB+nOKNHpu2MAu4UkabiKCMi7cNUH6wtIq1EpCROOe5BnD/AAUYDz4hILXddcSJSKZtxzgQuEJGe7vscKSKNRaSueyX7ITBIREqLyIU4/0lkuWxutpvN+AOVzHSOF8M5py8RkSEi8jcAEfm7iLwnIhVEpIOIdBORiu4xbAJcjlO+7odGInK9+/m7H+ccC7atEcBzInK+G3O0iHTI4TZHAANFpJ67rvIi0tnjsr8BVUWkRMC4csAuN+k3AW4KmLYd5/wM/AwH+gTnfb9JRIqLSFfgQpzzISt5nQvCOl0S/yc4SSP9MQinbK4UzlX5Ak7+1n0c54ryd5w/v94ntA44xUXzxLmha5+IfOpO64/z59ICEfkD5yqsdoj1vAA8JiK7RaSfqm5w1/0Izkm1AedPzRy9L6q6BadM8RICqu2p6o/AP3C+vHYA1wLXquqfOdzO9zi/mM7BKWNNH5+M8yfgqzjHdQ1OGWgoJXF+Vu/A+Rl7NjDQnfYSzhfWbJwaFm/hvJ/ZiXMvzhdhN5yrsa04vyhKurPcg/NzeStOufvYbCybm+3mxD5OPMdbqepanP+UYoDlIrIHp/w7Gef/nt9x3o/VOMfwPWBIDq+OvfgfzgVB+h+c1wf8Og70fzgVImaLyF6cz2fTnGxQVafhHNuJ7udvGc7/Gl58hVPTbauI7HDH3QU87cb1BM45mL6tAzjFvd+6n+GLM8WyE6dW20M4xXD/xqlpt4Ms5HUuyEp6FThjjMkxERmEUxniHwUdi8na6XLFb4wxxiNL/MYYU8RYUY8xxhQxdsVvjDFFTKGox3/WWWdpTExMQYdhjDGFSkpKyg5Vjc48vlAk/piYGJKTkws6DGOMKVREJPOdxIAV9RhjTJFjid8YY4oYS/zGGFPEZFnG77Zj0wLn9vyDOLdEf66qv/sc22npyJEjbNy4kUOHDhV0KMaY00RUVBRVq1YlMjLS0/whE7+I9AHuxelMIQWnKdco4FKgv4gsAx5X1V9zHXURsnHjRsqVK0dMTAzyV1POxhiTI6rKzp072bhxI9WrV/e0TLgr/tJAc1U9GGyi28xxLZzuyIxHhw4dsqRvjMkzIkKlSpXYvn171jO7Qpbxq+prwJ8i8kCI6amq+mX2wzSW9I0xeSm7OSXsn7tum+U3hZvHGGNM4eLlBq55IvIqTvvu+9NHqupi36IqQmIGfJyn61s/uH3IaTt37qR169YAbN26lYiICKKjnZv6Fi5cSIkSJUIumx27d+/m/fff56677sqT9QW69dZbefDBB7nwwgvzdL07duygc+fO7Ny5k6ioKL766ivKlg3es93WrVu5//77WbRoERUqVKBy5coMGzaMzZs3M3ToUGbO9NLvhn+eeOIJLrvsMq644ooCjcOcurwk/nj3+emAcQq0yvNojK8qVapEamoqAIMGDaJs2bL069cvz7eze/duXn/9dV8S/+jRo4OOP3bsGBERXrs6Ptkbb7zBZZddxlNPPcXmzZtDfgmqKp06daJXr15MnDgRgCVLlvDbb78Fnd+PWLPy9NNPBx3v93Z9tfmHgo6g4JzTIOt5sinLevyq2jLIw5L+aeLLL7+kQYMGxMbGcvPNN3P48GHAaSZj4MCBxMfHk5CQwOLFi2nTpg01a9ZkxIgRAOzbt4/WrVvTsGFDYmNj+d///gfAgAEDWLt2LfHx8Tz88MMADBkyhMaNGxMXF8eTTzq9YK5fv566dety2223Ua9ePa666ioOHjzIqlWraNKkSUaM69evJzbW6TI4MTExo/mOsmXL8tBDD3HRRRcxf/58YmJi2LHD6ewoOTmZxMREAObOnUt8fDzx8fE0aNCAvXv3klmJEiXYuHEjAOecc07IxD9nzhwiIyO58847M8ZddNFFtGjRIuOY3HjjjdSpU4cePXqQ3vptTEwM/fv3p2HDhgwePJiGDRtmLL969eqM4QEDBnDhhRcSFxeX8aW8fv16WrVqRVxcHK1bt+bXX39lz549nH/++Rw/7vRUuX//fs477zyOHDlC7969+eCDD07a7pQpU044fjt27CC9Dazly5fTpEkT4uPjiYuLY/Xq1UH335weskz8IlJZRN5K72pQRC4UkVv8D8347dChQ/Tu3ZtJkyaxdOlSjh49yhtvvJExvVq1aqSmptKiRYuMZLJgwYKMxB0VFcW0adNYvHgxc+bM4aGHHkJVGTx4MDVr1iQ1NZUhQ4Ywe/ZsVq9ezcKFC0lNTSUlJYWvv/4acJLe3XffzfLly6lQoQJTp06lTp06/Pnnn/z8888ATJo0ia5du54U//79+2natClLlizh0ksvDbmfQ4cO5bXXXiM1NZVvvvmGUqVO7sGxZs2afPjhhxlfaqEsW7aMRo0ahZz+ww8/MGzYMFasWMG6dev49tu/+hOvVKkSixcv5tFHH6V8+fIZv77Gjh1Lnz592LlzJ9OmTWP58uWkpaXx2GOPAXDvvffSq1cv0tLS6NGjB3379qV8+fLEx8czd+5cAGbOnEmbNm2C1uNO3263bt1Cxj1ixAjuu+8+UlNTSU5OpmpVX7p6NacIL3fujgM+w7mBC+AnnI6UTSF37NgxqlevzgUXXABAr169MhIywHXXXQdAbGwsTZs2pVy5ckRHR1OyZEl2796NqvLII48QFxfHFVdcwaZNm4IWecyePZvZs2fToEEDGjZsyKpVqzKuKKtXr058fDwAjRo1Yv369QB06dKFSZOcboNDJf6IiAhuuOGGLPezefPmPPjggwwfPpzdu3dTvPiJJZybNm3ihRdeYM2aNYwePZqpU6cCEBcXx549e7Jcf6AmTZpQtWpVihUrRnx8fMb+ACfsw6233srYsWM5duwYkyZN4qabbqJ8+fJERUVxyy238OGHH1K6dGkA5s+fz003OXUsevbsybx58zLWl36MJk6cGPQYZd5uKM2aNeP555/nxRdf5Jdffgn65WhOH14S/1mqOhmnd3lU9ShwzNeozCmhZEmnb/BixYplvE4fPnr0KOPHj2f79u2kpKSQmppK5cqVg96RrKoMHDiQ1NRUUlNTWbNmDbfccssJ2wAnkR89ehRwktXkyZP56aefEBFq1ap10nqjoqJOKLMuXrx4RtFHYBwDBgxg9OjRHDx4kObNm7Nq1aoT1vPtt98SGxtLpUqV+Pjjj3nyySd5/fXXiYmJoXz58ifMW69ePVJSUrI8Zpn3B6BMmTIZr2+44QY+/fRTZs6cSaNGjahUqRLFixdn4cKF3HjjjcycOZOrr7465HbA+WKeNWsWu3btIiUlhVatgpfABm431DG66aabmDFjBqVKlaJdu3Z89dVXYbdtCjcviX+/iFTC+UMXt2f57F0GmVNSREQE69evZ82aNQC8++67XH755Z6X37NnD2effTaRkZHMmTOHX35xWoAtV67cCeXobdq0YcyYMezbtw9wrrC3bdsWdt01a9YkIiKCZ555xtMVKzjl2elJOf2qHWDt2rXExsbSv39/GjdufFLij4uLY86cOWzevJnKlSvz8ssvc/fdd2dcZQdq1aoVhw8fZuTIkRnj0tLS+OabbzzFmC4qKoo2bdrwr3/9iz59+gDO/wN79uyhXbt2vPzyyyxZsgSASy65JOOP5PHjx2f8n1C2bFkaN27MfffdxzXXXOPpj9vAY5T+PwDAunXrqFGjBn379qVDhw6kpaVla39M4eKlVs+DwAygpoh8C0QDN/oaVRESrvql36Kiohg7diydO3fm6NGjNG7c+IQ/LbPSo0cPrr32WmJjY0lISKBOnTqAU6bcvHlz6tevT9u2bRkyZAgrV66kWbNmgJOw3nvvvSwTVdeuXXn44Yczyvqz8uSTT3LLLbfw+OOPZ/yxCzBs2DDmzJlDsWLFqFevHm3btj1huTp16vDcc89llJFXrlyZiRMnMmDAABo2bJhRFAbOjTLTpk3j/vvv58UXXyQqKoqYmBiGDRvGpk2bPMWZrkePHkybNo2rrroKgL1799KhQwcOHTqEqvLSSy8B8Morr9CnTx+GDBlCdHQ0Y8eOPeEYde7cmaSkJE/b7NevH126dGHkyJG0b//XuTd58mTeffddIiMj+dvf/sYjjzySrX0xhYunPndFpDhQGxDgR1U94ndggRISEvR06Yhl5cqV1K1bt6DDMKeAoUOHsmfPHp555pmCDuXUZ9U5sxQst4hIiqomZJ7XS+ucnYFZqrpcRB4DGorIs3YDlzE516lTJ9auXWtl6aZAeCnqeVxVp4jIpUBrYCjwBtDU18iMOY1NmzatoEMwRZiXP3fTa/C0B0ap6sdA3tzbb4wxJt95SfybRORNoCvwiYiU9LicMcaYU5CXBN4F5wauNqq6GzgTeDirhUTkPBGZIyIrRGS5iNznjh8kIptEJNV9tMvNDhhjjMkeL2X8VYCPVfWwiCQCccA7HpY7CjykqotFpByQIiKfu9NeVtWhOQnYGGNM7ni54p8KHBORvwMjgfOA97NaSFW3pNf8UdW9wErg3FzEenoaVD5vHx5s3bqVbt26UbNmTRo1akS7du346aefSEpK4pprrvF5h7P2xBNP8MUXX+T5eg8dOkTHjh2pX78+DRo0YN26dUHna9q0KfHx8VSrVo3o6OiMBt4Cm1/IK88//3y2l4mJiSE2Npa4uDguv/zyjBvnwLnZq6CF2qc+ffrw5ptvnjBu+vTptG3bluTkZPr27ZvnsQz67wjObdSG+Cu7Ub9VZ2bMnht2/pim7dmxy7/uxHfv2cvr4yZ7mveS63r7FoeXxH/cbabheuAVVX0Y51eAZyISAzQAvndH3SMiaSIyRkQqhljmdhFJFpHk7HQpZsJLb1Y4MTGRtWvXkpKSwgsvvJDtZoX99PTTTwdtSz63250yZQrly5dn2bJlfPXVV5x55plB5/v+++9JTU3l6aefpmvXrhlNTaS3ZJmXcpL4wWklNC0tjcTERJ599tmM8d99911ehZZjofape/fuGXcgp5s4cSLdu3cnISGB4cOHn7RMYJMXOfXAbT1I/XwiU958kZsffCqjyQo/hTpXd/+xl9ffmeJpHd/NGJeHEZ3IS+I/IiLdgX8C6T1MeOvKHRCRsji/Gu5X1T9wqoLWxGnnfwvw32DLqepIVU1Q1YT0zkJM7hXlZoVLlCjBpk2bUFUqVqxIhQoVPB+31NRULr74YuLi4ujUqRO//+5cFSYmJvLAAw+QkJBA3bp1WbRoEddffz21atXKaF0ToGPHjjRq1Ih69eplNPcwYMAADh48SHx8PD169ADgpZdeon79+tSvX59hw4ZlGVezZs1OuGM4vfOYpKQkEhMTg76Xn3zyCXXq1KFRo0b07ds341desOark5KSuOyyy2jfvj21a9fmzjvvzHjPJkyYQGxsLPXr16d///4h9yld69atWbVqFVu2bAGc9/yLL76gY8eOJ/zaHDRoED179qR58+b07NmTcePGcc+jgzPWc80/+5L0XTLHjh2j9/1PUr9VZ2Jbd+Hlke+FPVZ1a9WgePEIduzazYTps4ht3YX6rTrT/7n/O2neJ4a8wbBR4zOGHx38Kv83+n2SvkvmsutvoX3PvtRu0Yk7+z+XcTzK1mrOQ0+9xEVXdGV+Shovvfke9Vt1pn6rzhnrGvD8cNb+spH4K7vx8DMvAzDkjbdp3O4fxF3RhSeH/tU6btlazbN8L3PKS+LvAzQDnlPVn0WkOvCul5WLSCRO0h+vqh8CqOpvqnpMVY8Do4Am4dZh8lZRbla4Ro0aLF68mIEDB3o6VoH++c9/8uKLL5KWlkZsbCxPPfVUxrQSJUqQnJzMnXfeSYcOHXjttddYtmwZ48aNY+fOnQCMGTOGlJQUkpOTGT58ODt37mTw4MGUKlWK1NRUxo8fT0pKCmPHjuX7779nwYIFjBo1ih9+CH/H6qxZs+jYsWPQacHey0OHDnHHHXfw6aefkpKSckIH3aGar164cCGvvPIKK1asYO3atXz44Yds3ryZ/v3789VXX5GamsqiRYuYPn36SfsUKL011cmTnaKOjz76iMTERM4444yTYl+xYgVffPEFEyZMCLnvqct/ZNPWbSz7agpLv5xMn64dwh6r7xcvpVgx4ciRo/R/bjhfTX6T1NkTWJS6nOmz5pww783dOvDOB07veMePH2fijNn843qnHsrC1OW88uy/WZH0AWt/2ciHnzg34e0/cJCmDeqz5ItJlIoqydjJM/h+5jss+OhtRr0/jR+WrWLwI32peX5VUj+fyJDHH2D23Pms/vlXFn78LqmzJ5KStpKvF5zcCGC4z2VOeOmIZQXQH0gvr/9ZVV/Majlxev99C1ipqi8FjA8sJuoELMtu0MY/p2uzwgcPHqRPnz78+OOPpKamZlxNt2/fnmXLwp+Ce/bsYffu3RkN2IVrvrpevXpUqVKFkiVLUqNGDTZs2ADA8OHDueiii7j44ovZsGFD0F8k8+bNo1OnTpQpU4ayZcty/fXXh2z8rWXLlpx77rl8+umndO/ePeg8wd7LVatWUaNGDapXrw5wwrKhmq9u0qQJNWrUICIigu7duzNv3jwWLVpEYmIi0dHRFC9enB49epxwTEIJLO5JL+YJ5rrrrsuyaega1aqy7tdN3PvYi8ya8y1nlCsTdL6XR40n/spu9HvmZSa9MZjktBUkNmtEdKWKTuzXt+PrBSc2RBBz3jlUqlieH5atYvbc+TSoV5tKZ1Zwjkd8PWqcX9U5Hh3bMG+h8+UcERHBDe2drk3nLUyl09UtKVO6FGXLlOb6tq345vuTv8Rnz13A7LkLaHBVdxq2uYlVa9ez+ucNJ80X7nOZE146YrkWSAVmucPxIjLDw7qbAz2BVpmqbv5HRJaKSBrQEnggx9GbbCuqzQovXbqUs846i+joaKZOncqkSZMYPnw4u3btol69emHjzEpWzVcnJSXxxRdfMH/+fJYsWUKDBg2CNl+dHemtocbHx2d0jBMqLjj5vQwmVPPVzjXcXzIPZ8cll1zCli1bWLJkCd99990JDcUFOvk9/6to49DhPwGoWOEMlnw+kcRmCYx4dyq39gve5WR6Gf8308bQomnDoPMEc2v3joyb/BFjJ83g5m5//ZoIdTyiSpbIdteWqsrAe/qQ+vlEUj+fyJpvZ3BL944nzZfd9zIrXop6BuEUx+x2A00FamS1kKrOU1VR1ThVjXcfn6hqT1WNdcdfp6pbchG/yaai2qxwrVq1WLVqFcuXL6dMmTK89dZb9OvXjw4dOmSZyMqXL0/FihUzjlFOmq+uWLEipUuXZtWqVSxYsCBjWmRkJEeOOG0etmjRgunTp3PgwAH279/PtGnTMo5VMMWLF2fYsGG888477Nq1y1MstWvXZt26dRlXjOm/uCB089ULFy7k559/5vjx40yaNIlLL72UJk2aMHfuXHbs2MGxY8eYMGFCxjEJ3KfMRISuXbvSq1cv2rZtS1RUVJYxx8TEkLr8R44fP86GTVtZmLocgB27fuf48ePc0L41z/77LhYvXZXFmhxN4usxd0EKO3b97sQ+fRaXNzu5+LNT21bMmvMdi5asoE1is4zxC1OX8/Ovm5zjMWM2lzY5uRG1Fk0bMP2zORw4eJD9Bw4ybdYcWjRtQLkypdm7b3/GfG0SmzFm0gz27T8AwKYt29i2w9t7mRte6vEfUdU9mT4c/v8tXlQMyt+uDYpqs8IVK1bk7bffpmfPnqgq5cuXZ/z48QwcOJDLLrssy2qQb7/9NnfeeScHDhygRo0aJ+xDVq6++mpGjBhB3bp1qV27NhdffHHGtNtvv524uDgaNmzI+PHj6d27d0Z/w7feeisNGoRvmbFKlSp0796d1157jccffzzLWEqVKsXrr7/O1VdfTZkyZWjcuHHGtGDNV8+fP5/GjRtzzz33sGbNGlq2bEmnTp0oVqwYgwcPpmXLlqgq7du3p0OHDkH3KbPu3bvzn//8h8GDB580LZjmzZtTvdq5XJh4I3VrVadhrNP896Yt2+nz4KCMX4AvDLzX0/qqVI5m8CN9adn5Dif21pfSoU3iSfOVKBFJy0sSqFC+3AkXJY0vupB7Hn2RNes30PKSBDq1bXnSsg1j69K783U0af9PwPn10KC+E3fzxvHUb9WZti0vYcjjD7By9c80c6tuli1divdeeZazzwpe4yyvZNkss4i8BXwJDABuAPoCkarqveH2XLJmmU9N1qxw4bRv3z7Kli2LqnL33XdTq1YtHnggeIlrUlISQ4cOZebMmUGn55sCaJb5+PHjNGxzE1Pe/A+1alQDIOm7ZIaOeIeZ75xc9dQ3PjTL7KWo516gHnAYmAD8gfW5W+R16tSJd955h/vuu6+gQzHZNGrUKOLj46lXrx579uzhjjvuKOiQTjkrflrH35t3oPWlTTKS/unEU0csACJyBqDuXbj5yq74jSnirCOWLOXpFb+INBaRpUAasFRElohI6IrgJku5vfnCGGMCZTeneCnqeQu4S1VjVDUGuBvw/s+WOUFUVBQ7d+605G+MyROqys6dOz3VkErnpVbPMVXNqOunqvNEJPcNaBRRVatWZePGjVj7Q8Zkw+5tBR1BwdmzMstZoqKigt6tHkrIxC8i6Xc6zHU7YpkAKE6HLEmet2BOEBkZmXHXpDHGo0EXZz3P6cqHKt/hrvgzN54WeHuglVMYY0whFTLxq+rJdyUYY4wp9KzvXGOMKWIs8RtjTBHjpR5/SS/jjDHGFA5ervjnexxnjDGmEAhXnfNvOJ2jlxKRBkB685xnAKXzITZjjDE+CFedsw3QG6iKU7UzPfHvBR4JsYwxxphTXLjqnG8Db4vIDao6NR9jMsYY4yMvZfxVReQMcYwWkcUicpXvkRljjPGFl8R/s6r+AVwFVMLpR9db1znGGGNOOV4Sf3rZfjvgHVVdHjDOGGNMIeMl8aeIyGycxP+ZiJTD+tw1xphCy0uzzLcA8cA6VT0gIpWAPr5GZYwxxjfh6vHXUdVVOEkfoIaIlfAYY0xhF+6K/0Hgdk5unhmcZplb+RKRMcYYX4Wrx3+7+2zNMxtjzGnEWuc0xpgixhK/McYUMWETv3u37nn5FYwxxhj/hU38qqrAJ/kUizHGmHzgpahnsYg0zu6KReQ8EZkjIitEZLmI3OeOP1NEPheR1e5zxWxHbYwxJse8JP6mwHwRWSsiaSKyVETSPCx3FHhIVS8ELgbuFpELgQHAl6paC/jSHTbGGJNPvNy52yYnK1bVLcAW9/VeEVmJ07FLByDRne1tIAnon5NtGGOMyb4sr/hV9RfgPKCV+/qAl+UCiUgM0AD4HqjsfikAbAUqh1jmdhFJFpHk7du3Z2dzxhhjwvDS2fqTOFfkA91RkcB7XjcgImWBqcD9bvPOGdw/jzXYcqo6UlUTVDUhOjra6+aMMcZkwcuVeyfgOmA/gKpuBsp5WbmIROIk/fGq+qE7+jcRqeJOrwJsy27Qxhhjcs5L4v8z8MpcRMp4WbE4Lbq9BaxU1ZcCJs0AermvewH/8x6uMcaY3PKS+CeLyJtABRG5DfgCGOVhueY4vXW1EpFU99EOp/euK0VkNXAF1puXMcbkqyxr9ajqUBG5EvgDqA08oaqfe1huHqF76mqdrSiNMcbkGS/VOXETfZbJ3hhjzKkvy8QvIns5uebNHiAZ5watdX4EZowxxh9erviHARuB93GKbroBNYHFwBj+uhnLGGNMIeDlz93rVPVNVd2rqn+o6kigjapOAqydHWOMKWS8JP4DItJFRIq5jy7AIXda0JuvjDHGnLq8JP4eONUytwG/ua//ISKlgHt8jM0YY4wPvFTnXAdcG2LyvLwNxxhjjN9CJn4R+beq/kdEXiFIkY6q9vU1MmOMMb4Id8W/0n1Ozo9AjDHG5I+QiV9VPxKRCCBWVfvlY0zGGGN8lFWfu8dw2twxxhhzmvByA1eqiMwApuA2zQwQ0MyyMcaYQsRL4o8CdgKtAsYpYInfGGMKIS/VOfvkRyDGGGPyR7jqnEGrcaaz6pzGGFM4hftzNxlIwSnqaQisdh/xQAnfIzPGGOOLcNU53wYQkX8Bl6rqUXd4BPBN/oRnjDEmr3lpq6cicEbAcFmsVU5jjCm0vNTqGQz8ICJzcNrjvwwY5GdQxhhj/OOlVs9YEfkUaOqO6q+qW/0NyxhjjF+yLOoREQGuAC5S1f8BJUSkie+RGWOM8YWXMv7XgWZAd3d4L/CabxEZY4zxlZcy/qaq2lBEfgBQ1d9FxKpzGmNMIeXliv+I20qnAohINHDc16iMMcb4xkviHw5MA84Wkedwet163teojDHG+MZLrZ7xIpICtMapztlRVVdmsZgxxphTVLi2es5Q1T9E5EycjtYnBEw7U1V35UeAxhhj8la4K/73gWtw2usJbKxN3OEaPsZljDHGJ+Ha6rnGfa6ef+EYY4zxW7iinmrhFlTVX8NNF5ExOL8YtqlqfXfcIOA2YLs72yOq+kl2AjbGGJM74Yp6PsYp0pGAcQpEA2cDEVmsexzwKvBOpvEvq+rQ7IVpjDEmr4Qr6okNHBaRGKA/TvMNWVbnVNWv3WWMMcacQry01VNLRMYBn+L80Xuhqr6Si23eIyJpIjJGRKx5Z2OMyWchE7+I1BeRCcBU4AugvqqOVtUjudjeG0BNnF68tgD/DbP920UkWUSSt2/fHmo2Y4wx2RSujH8JsAGnrL8J0MRpqNORkz53VfW39NciMgqYGWbekcBIgISEhJB9/xpjjMmecIn/FsJ0tp4TIlJFVbe4g52AZXm5fmOMMVkL9+fuuNys2C0mSgTOEpGNwJNAoojE43yhrAfuyM02jDHGZF+4evyjgP9T1ZOuykWkDNAVOKyq44Mtr6rdg4x+K6eBGmOMyRvhinpeA54QkVicIpntQBRQC6fz9TFA0KRvThGDyhd0BAVj0J6CjsCYU1q4op5UoIuIlAUSgCrAQWClqv6YP+EZY4zJa16aZd4HJPkfijHGmPzgpSMWY4wxpxFL/MYYU8RY4jfGmCImyzJ+EbkAeBg4P3B+VW3lY1zGGGN8kmXiB6YAI4BRwDF/wzHGGOM3L4n/qKq+4Xskxhhj8oWXMv6PROQuEakiImemP3yPzBhjjC+8XPH3cp8fDhhnna0bY0wh5eUGLuts3RhjTiNeeuAqLSKPichId7iWiFzjf2jGGGP84KWMfyzwJ3CJO7wJeNa3iIwxxvjKSxl/TVXtKiLdAVT1gAR2xWWMyRcxAz4u6BAKzPqogo7g9OLliv9PESmF2xuXiNQEDvsalTHGGN94ueJ/EpgFnCci44HmQG8/gzLGGOMfL7V6PheRxcDFgAD3qeoO3yMzxhjji3BdLzbMNCq9k/RqIlJNVRf7F5Yxxhi/hLvi/6/7HIXTA9cSnCv+OCAZaOZvaMYYY/wQ8s9dVW2pqi1xrvQbqmqCqjYCGuBU6TTGGFMIeanVU1tVl6YPqOoyoK5/IRljjPGTl1o9aSIyGnjPHe4BpPkXkjHGGD95Sfx9gH8B97nDXwPWTLMxxhRSXqpzHgJedh/GGGMKOS9dLzYHBnFy14vWLLMxxhRCXop63gIeAFKwrheNMabQ85L496jqp75HYowxJl94SfxzRGQI8CEBjbPZnbvGGFM4eUn8Td3nhIBxCrTK+3CMMcb4zUutnpY5WbGIjAGuAbapan133JnAJCAGWA90UdXfc7J+Y4wxOeOlVs8Twcar6tNZLDoOeBV4J2DcAOBLVR0sIgPc4f7eQjXGGJMXvDTZsD/gcQxoi3PFHpaqfg3syjS6A/C2+/ptoKPHOI0xxuQRL0U9/w0cFpGhwGc53F5lVU1v3nkrUDnUjCJyO3A7QLVq1XK4OWOMMZl5ueLPrDRQNbcbVlXF7c4xxPSRbougCdHR0bndnDHGGJeXMv6l/JWgI4Bo4Jkcbu83EamiqltEpAqwLYfrMcYYk0NeqnN24K87do8Cv+F8AeTEDKAXMNh9/l8O12OMMSaHvBT1PKaqv7iPTUBJ4JOsFhKRCcB8oLaIbBSRW3AS/pUishq4wh02xhiTj7xc8W8UkddV9S4RqQh8DIzKaiFV7R5iUuvsBGiMMSZvZXnFr6pPAPtEZAQwG/ivqo71PTJjjDG+CHnFLyLXBwx+DzwOLARURK5X1Q/9Ds4YY0zeC1fUc22m4R+ASHe84jTaZowxppAJmfhVtU9+BmKMMSZ/ZFnGLyIXiMiXIrLMHY4Tkcf8D80YY4wfvFTnHAUMBI4AqGoa0M3PoIwxxvjHS+IvraoLM4076kcwxhhj/Ocl8e8QkZq4zTaIyI3AlvCLGGOMOVV5uYHrbmAkUEdENgE/Az18jcoYY4xvvDTLvA64QkTKAMVUda//YRljjPFLyKIeEblWRM4PGPUQME9EZohIdf9DM8YY44dwZfzPAdsBROQa4B/AzTgtbI7wPzRjjDF+CJf4VVUPuK+vB95S1RRVHY3TJr8xxphCKFziFxEpKyLFcFrU/DJgWpS/YRljjPFLuD93hwGpwB/ASlVNBhCRBlh1TmOMKbTCtdUzRkQ+A84GlgRM2gpYOz7GGFNIha3O6fa4tSnTOLvaN8aYQszLDVyFWsyAjws6hAKz3v6JMcYE4aXJBmOMMacRT1f8IhIBVA6cX1V/9SsoY4wx/sky8YvIvcCTwG/AcXe0AnE+xmWMMcYnXq747wNqq+pOv4MxxhjjPy9l/BuAPX4HYowxJn94ueJfBySJyMfA4fSRqvqSb1EZY4zxjZfE/6v7KOE+jDHGFGJe2uN/Kj8CMcYYkz9CJn4RGaaq94vIR7jdLgZS1et8jcwYY4wvwl3xv+s+D82PQIwxxuSPcI20pbjPc/MvHGOMMX4rkLZ6RGQ9sBc4BhxV1YSCiMMYY4qigmykraWq7ijA7RtjTJGU5Q1cItLZyzhjjDGFg5c7dwd6HJcdCswWkRQRuT3YDCJyu4gki0jy9u3bc7k5Y4wx6cJV52wLtAPOFZHhAZPOAI7mcruXquomETkb+FxEVqnq14EzqOpIYCRAQkLCSdVJjTHG5Ey4K/7NQDJwCEgJeMwA2uRmo27PXqjqNmAa0CQ36zPGGONduOqcS4AlIjJeVXN7hZ9BRMoAxVR1r/v6KuDpvFq/McaY8MIV9UxW1S7ADyIS7M7dnLbHXxmYJiLp239fVWflcF3GGGOyKVx1zvvc52vycoOqug64KC/XaYwxxruQZfyqusV9/sUdVct9vQ3YlQ+xGWOM8YGXevy3AR8Ab7qjqgLTfYzJGGOMj7zU478baA78AaCqq4Gz/QzKGGOMf7wk/sOq+mf6gIgUJ0gzzcYYYwoHL4l/rog8ApQSkSuBKcBH/oZljDHGL14S/wBgO7AUuAP4RFUf9TUqY4wxvvHSOmcPYKKqjkofISLXqOpM/8IyxhjjFy9X/K8A34hI3YBxdqetMcYUUl4S/8/AzcAHAc0xi38hGWOM8ZOXoh5V1cUicjkwQUSaAhE+x2WMMcYnXq740+/g3YHTKqcC9f0MyhhjjH9CJn4Redd9+Xn6OFU9rqoPq6qXLwxjjDGnoHAJvJGInAP0EZGKInJm4CO/AjTGGJO3wpXxjwC+BGrgdMAS+IeuuuONMcYUMuFa5xyuqnWBMapaQ1WrBzws6RtjTCHlpVbPCyJSLfNIVf3Vh3iMMcb4zEvi/xinaEeAKKA68CNQz8e4jDHG+CTLxK+qsYHDItIQuMu3iIwxxvgq29UyVXUx0NSHWIwxxuSDLK/4ReTBgMFiQENgs28RGWOM8ZWXMv5yAa+P4pT5T/UnHGOMMX7zUsb/VH4EYowxJn94KeqJBv6NU4snKn28qrbyMS5jjDE+8fLn7nhgFU41zqeA9cAiH2MyxhjjIy+Jv5KqvgUcUdW5qnozYFf7xhhTSHn5c/eI+7xFRNrj1OixRtqMMaaQ8pL4nxWR8sBDON0wngHc72dQxhhj/OOlVk96p+p7gJYAInK/jzEZY4zxUU47VHkw61mMMcacinKa+K2zdWOMKaRymvg1NxsVkatF5EcRWSMiA3KzLmOMMdkTsoxfRPYSPMELUCqnGxSRCOA14EpgI7BIRGao6oqcrtMYY4x3IRO/qpYLNS2XmgBrVHUdgIhMBDoAlviNMSYfeKnOmdfOBTYEDG8kSDPPInI7cLs7uE9EfsyH2E4rAmcBOwo6jnz3lP0Fdbopsucy5PZ8Pj/YyIJI/J6o6khgZEHHUZiJSLKqJhR0HMbklp3LeSunf+7mxibgvIDhqu44Y4wx+aAgEv8ioJaIVBeREkA3YEYBxGGMMUVSvhf1qOpREbkH+AyIAMao6vL8jqOIsKIyc7qwczkPiWququQbY4wpZAqiqMcYY0wBssRvjDFFjCX+U4SIPCoiy0UkTURSRaSpiNwvIqXzcBvrReSsXCzfW0Rezat4TNEjIoNEpF82xn+XP5EVLadsPf6iRESaAdcADVX1sJucSwCTgPeAAwUUV4SqHiuIbRsDoKqXFHQMpyO74j81VAF2qOphAFXdAdwInAPMEZE5ACLyhogku78Mnkpf2L2Sf0pEFovIUhGp446vJCKz3flHE9CqqohMF5EUd9rtAeP3ich/RWQJ0ExE+ojITyKyEGieD8fCnGbcX7M/icg8oHY2l93nPieKSJKIfCAiq0RkvIiIO62RiMx1z+fPRKSKO76viKxwf0VPzPMdK8xU1R4F/ADKAqnAT8DrwOXu+PXAWQHznek+RwBJQFzAfPe6r+8CRruvhwNPuK/b4zS6d1amdZUCluH0rYw7Txf3dRXgVyAa5xfIt8CrBX287FF4HkAjYClQGqf3vjVAvyDzDQoxfp/7nIjTGVRVnAvW+cClQCTwHRDtztcVp4o4ON3ElnRfVyjoY3EqPayo5xSgqvtEpBHQAqeXs0khmqvu4l6dF8dJyhcCae60D93nFOB69/Vl6a9V9WMR+T1gXX1FpJP7+jygFrATOAZMdcc3BZJUdTuAiEwCLsjNvpoipwUwTVUPAIhIbm7WXKiqG931pAIxwG6gPvC5+wMgAtjizp8GjBeR6cD0XGz3tGOJ/xShTll6EpAkIkuBXoHTRaQ60A9orKq/i8g4ICpglsPu8zGyeF9FJBG4AmimqgdEJClgXYfUyvXNqelwwOv081yA5araLMj87XEufq4FHhWRWFU96n+Ypz4r4z8FiEhtEakVMCoe+AXYC6Q3j30GsB/YIyKVgbYeVv01cJO7jbZARXd8eeB3N+nXAS4Osfz3wOXufwWRQGfve2UM4JyDHUWklIiUw0nCeelHINqtIIGIRIpIPREpBpynqnOA/jjnfNk83nahZVf8p4aywCsiUgE4ilMOejvQHZglIptVtaWI/ACswmnW+lsP630KmCAiy3HKQX91x88C7hSRlTgfnAXBFlbVLSIyCKc8dTfO/xDGeKaqi90iwiXANpy2ukJ5TETuD1i2qof1/ykiNwLDRaQ8Tk4bhvN/2XvuOAGGq+runO7H6caabDDGmCLGinqMMaaIscRvjDFFjCV+Y4wpYizxG2NMEWOJ3xhjihhL/MYYU8RY4jfGmCLm/wGdOOgU8pSfeQAAAABJRU5ErkJggg==\n"
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "labels = [\"Standard\", \"d Lines\"]\n",
    "\n",
    "x = np.arange(len(labels))  # the label locations\n",
    "width = 0.35  # the width of the bars\n",
    "\n",
    "fig, ax = plt.subplots()\n",
    "\n",
    "rects1 = ax.bar(x - width/2, results_time, width, label='Tomatenvirus & Chrysovirus')\n",
    "rects2 = ax.bar(x + width/2, results2_time, width, label='Chrysovirus & Tomato Ringspot Virus Polyprotein')\n",
    "\n",
    "# Add some text for labels, title and custom x-axis tick labels, etc.\n",
    "ax.set_ylabel('Laufzeit in Sekunden (Niedriger ist besser)')\n",
    "ax.set_title('Laufzeiten von verschiedenen LCS Implementationen')\n",
    "ax.set_xticks(x, labels)\n",
    "ax.legend()\n",
    "\n",
    "#ax.bar_label(rects, padding=3)\n",
    "#fig.tight_layout()\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 48,
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "\n",
      "\n",
      "\n",
      "\n",
      "\n",
      "\n",
      "330.953125\n"
     ]
    },
    {
     "data": {
      "text/plain": "<Figure size 432x288 with 1 Axes>",
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZgAAAEICAYAAABiXeIWAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAAiqElEQVR4nO3deZgcVb3/8fcnC4FgQgRiDEkgyCIigmIkgKgILixiuIosogREI1cU+F1QlutPwe2iPpf9iiKoQXZZIyCCSFxAwER2AhIxOAmELCQxIWwJ3/vHOXNTabp7ajKpmenx83qefrrq1Kmqb1VX1bfrVFW3IgIzM7O1rV9PB2BmZn2TE4yZmVXCCcbMzCrhBGNmZpVwgjEzs0o4wZiZWSX6VIKRdKikW0vUO1XSJd0R05qQFJK27AVxTJX02Z6Oo7dp9vmU3QYbjHu4pD92LTqrR9IsSR/o6TjWtq5sb92hRxKMpN0k3SVpiaTnJN0p6V1dnW5EXBoRH1obMZqtib62DUoamxPqgAbDt5b0C0kL8v78oKT/kNQ/Dz9S0mOSlkp6VtLNkoY0mNa/1Beazn6hqPdZ9PbtrdsTjKShwI3AucCGwCjgNOCl7o5lbWi043XXtNbm/P+VeL11naQtgHuANuBtEbEB8AlgHDBE0vuA7wCHRMQQ4C3AlT0Vr/WAiOjWF2njW9xk+OHAncB5wBLgMWDPwvANgIuAZ4A5wLeA/oVx/1io+1bgNuA54FnglFx+KnAVcDGwFHgEGFcYbxPgGmA+8HfgmMKwU4GrgUuAfwJfA14ANizUeQewABiY+z8DzAAWAb8GNivUDeBo4Ang74WyY4An83S+D/SrWT9nAgvz8m8B/Db3LwAuBYbVzGPLQv/PgG8V+icA9+fl+RuwVy6fCnwzz28pcCuwcYPPbQbwkUL/gLz+dsz9H83reXGe7lsKdWcBJwAP5s/8SmDdOvMYlMffrlA2PK//N+T+j+RlWQzcBWxfM58T83xeyjGeSNqOlgKPk7c1oD9wSl4fS4HpwJjC+jwqf2aLgf8B1GAb3IZV2+DjwIGFYRsBU/J6vzev67Lj/izP96Yc3z3AFmtj3Jp1PjYv74A6wy4BbmqyL58AXN+JY8NU4LO5e3dgNvAVYB5pf98f2Af4a16uU+rsl1fmZfoLsEPNZ/+B3N0POCl/tgtJx4INa5b3CFLiXJQ/63fl7WYxcF5N3B3t36/ZVkjJ9kVgJbCMfEwE9gXuy9tEG3BqYVr/yNNbll+78NrtbVfgz6T96M/ArjXrt+H+DOxM2mcWAw8Au5cdt+FnWvbDX1svYGj+UCcDewOvrxl+OLAC+H/AQOCgvLLaN4DrgB8B6wNvIO2Yn6/duYEhpI3yeGDd3D++sDG+SNpY+wP/Bdxd2PimkxLHOsCbSAf6DxfGfYW0sfcD1iMd3D9XWIbvAz/M3ROAmXmDGgB8FbirZgO8jXQ2t16h7I5ctilph/pszfr5Up7eesCWwAdJB+DhwO+Bs2rmUTfBADvl9fvBvDyjgG0KG9XfgK3zfKYCpzf4XL8GXFro3xeYkbu3Bp7P8xhIOmjMBNYp7Pz3khL7hqSd9agG8/kJ8O1C/9HALbn7HaSD0fj8uU7M0x5UmM/9wJi8PG8m7cSbFA4uW+TuLwMP5ToCdgA2KqzPG4Fh+fOZz6qkfDirtsH18/SPyJ9V+xePbfPwK0gHt/WB7UiJruy4PyPtRzvl4ZcCV3R13DrreyyNE8xc4Igm+/p7SMn/NODd7Z9Dk/pTWT3BrCBtVwOBz+X1fBlpX35rnvbmNfvlAbn+CaQvh+1f8maxKsEcC9wNjCbtMz8CLq9Z3h+SjhsfIh0rricdb0aRtrH3dWL/7nBbKdTfHXgbaV/cnvTFeP9GnwWrb28bkpLcp3Msh+T+jQrrt+7+nJdrIemY2I+0ry4Ehnf2WNCjCSYH+xbSRj47b0RTgBGFFfY0+RthLrs3r7QRpG+e6xWGHQLcUWdlHwLc12D+pwK/KfRvC7yQu8cD/6ipfzLw08K4v68Z/lngt7lbpJ37vbn/V8CRhbr9gOXkbzl5g9mjZnrRvhHm/i8AtxeW8R/1lqtQf//istM8wfwIOLPJDv/VmjhuaVB3S9I3m8G5/1Lga7n7/wNX1ayDOeRvSKSd/1OF4d8jJ+g68/kA8LdC/53AYbn7fOCbNfUfZ9XBYBbwmZqY5+VpDqwz3oQGMQSwW6H/KuCkOtvgQcAfasb9EfB1UgJ8hZzM87DvlBm38BleWBi2D/BYV8ets6xjaZxgXqGwnTYYf2/gl6RvxcuAM8gtDg22t2KCeYFVrRNDchzjC/Wns+rgeyr5S2JhG3sGeE/hs29PMDNYvVVkZF6WAYXlHVUYvhA4qNB/DXBclN+/O9xWmqy/s8j7Z73PgtW3t08D99aM/yfg8I72Z9KZ/M9rxv01MLGjcZu9euQif0TMiIjDI2I06ZvbJqQV2W5O5KXInsp1NiN9O3lG0mJJi0k7zhvqzGYMKeM2MrfQvRxYN7fLbwZs0j79PI9TSMmtXVvNtK4BdpE0Engv8CrwhzxsM+DswrSeIyWhUU2mV1vWvvx160saIekKSXMk/ZPUdLFx/cV+jc6up9fVqxQRM0k77n6SBpOaxC7LgzfJy9Be99W8DMV1UGo+pDO7wZLGSxoLvJ10VgtpXR9f89mNocG6yzEfRzo4zcvrsL3u2lgvmwHja+I5FHgj6UxzAK/9nMuM21EMXRm3MxaSDs4NRcSvImI/0rfrCaQDYtkL+QsjYmXufiG/P1sY/gKrx138bF8lfYEtfvbtNgOuK6ybGaSmquI+XjufRvMts3+XXtd5u75D0nxJS0jNa2X35dX2s+ypkrFsBnyiZpvZjdU/305vMz1+m3JEPEb6RrVdoXiUJBX6NyWd1bSRzmA2johh+TU0It5aZ9JtpOatzmojXQsZVngNiYh9imHXLMMiUpvkQcAnSc0NUZje52umt15E3NVoetmYQnf78jeq/51c9raIGAp8irSRt1sODC70Fw80baRrOGvD5aQzxwnAo/kADin2zdor5c92DOksplPyAeeqPJ9DgBsjYmke3EZqPiuu68ERcXlxEjXTuywidsvxBfDdwrS6ul7agN/VxPO6iPh3UlPJCl77OZcZtyvzXZt+A3y8TMWIeDUibic1J2/XUf019H/rUlI/UhPY03XqtQF716yfdSOi09sj5fbvRurt95eRWnTGRLpp4oes2pfr1S9abT/LNqXcftZGOoMpLsf6EXF6iXEb6om7yLaRdLyk0bl/DOlAcXeh2huAYyQNlPQJUpPazRHxDOlA/t+ShkrqJ2mLfLdKrRuBkZKOkzRI0hBJ40uEeC+wVNKJktaT1F/SdiVuo74MOIzUBnxZofyHwMmS3pqXd4O8TB35sqTX5/VzLM3vvhlCan5YImkU6fpB0f3AJ/Oy7AUU19dFwBGS9szrc5SkbUrEV88VpDbrf2f1dXAVsG+ex0DSdbGXSBcU18RlpGR+aM18fgwclb8FStL6kvZtclvsmyXtIWkQqZ39BdLZJ8CFwDclbZWntb2kjToZ543A1pI+nbflgZLeJektOVFeC5wqabCkbUnXjDoctyvz7WT8RYMkrVt49SM19e0q6fuS3gggaUtJl0gaJmmCpIPzdixJO5G2vbubzagL3inpY0otEceRtrF68/oh8G1Jm+WYh0uasIbzXNP9G9JZ0WhJ6xTKhgDPRcSLeX19sjBsPmn7bPTF+WbS5/5JSQMkHURq/r+xRCyXkFofPpyPE+tK2r39OL2meuIMZinpOsc9kp4nbQAPkw467e4BtiJdmPw2cEBELMzDDiNdfH+UdAHrauqcpudvtR8E9iOd2j0BvL+j4PKO/xFS08vfcwwXku5ea2ZKjnluRDxQmN51pG/FVyg1Xz1MapfuyA2kNub7SXf7XNSk7mnAjqSL9TeRDlxFx5LWw2LSQfn6Qnz3ki4Gn5nH/x2v/RZUSv4C8CfSnSxXFsofJ51VnUtan/sB+0XEy2s4n3tINw1sQmoDby+fRroYfB5p25hJapJpZBBweo5pLumLzcl52BmkxHgr6Y6ei0gXNzsT51JSwj2Y9O1yLmlbGJSrfJHUzDCXdBb/006M25X5rollpATc/tojIv5GupNpLPBIbtK5BphG2s8XkT6PJ0jr8BLg+xFxaRfiaOYG0hePRaTrER+LiFfq1DubtL/eKmkp6RhU5svna3Rh/4Z0NvcIMFfSglz2BeAbOa6vkbbB9nktJx0P78zNWDvXxLKQdOw6ntR8+RXSnZ0L6EBEtJFaHk4hJbI20hfVLuWI9lsrew1Jh5Mu9O3W07GYWWuQdCrpRpZP9XQstkqPX4MxM7O+yQnGzMwq0euayMzMrG/wGYyZmVWipX/wb+ONN46xY8f2dBhmZi1l+vTpCyJieNXzaekEM3bsWKZNm9bTYZiZtRRJtU/8V8JNZGZmVgknGDMzq4QTjJmZVcIJxszMKuEEY2ZmlXCCMTOzSjjBmJlZJZxgzMysEk4wZmZWiZZ+kr8rxp50U0+H0LJmnb5vT4dgZi3AZzBmZlYJJxgzM6uEE4yZmVXCCcbMzCrhBGNmZpVwgjEzs0o4wZiZWSWcYMzMrBJOMGZmVolKE4ykYZKulvSYpBmSdpG0oaTbJD2R31+f60rSOZJmSnpQ0o5VxmZmZtWq+gzmbOCWiNgG2AGYAZwE3B4RWwG3536AvYGt8msScH7FsZmZWYUqSzCSNgDeC1wEEBEvR8RiYAIwOVebDOyfuycAF0dyNzBM0siq4jMzs2pVeQazOTAf+Kmk+yRdKGl9YEREPJPrzAVG5O5RQFth/Nm5bDWSJkmaJmna/PnzKwzfzMy6osoEMwDYETg/It4BPM+q5jAAIiKA6MxEI+KCiBgXEeOGDx++1oI1M7O1q8oEMxuYHRH35P6rSQnn2famr/w+Lw+fA4wpjD86l5mZWQuqLMFExFygTdKbc9GewKPAFGBiLpsI3JC7pwCH5bvJdgaWFJrSzMysxVT9h2NfAi6VtA7wJHAEKaldJelI4CngwFz3ZmAfYCawPNc1M7MWVWmCiYj7gXF1Bu1Zp24AR1cZj5mZdR8/yW9mZpVwgjEzs0o4wZiZWSWcYMzMrBJOMGZmVgknGDMzq4QTjJmZVcIJxszMKuEEY2ZmlXCCMTOzSjjBmJlZJZxgzMysEk4wZmZWCScYMzOrhBOMmZlVwgnGzMwq4QRjZmaVcIIxM7NKOMGYmVklGiYYSetKmijpo0pOlHSjpLMlbdydQZqZWetpdgZzMfAh4DPAVGBT4DxgKfCzqgMzM7PWNqDJsG0jYjtJA4DZEfG+XH6LpAe6ITYzM2thzc5gXgaIiBXA0zXDVlYWkZmZ9QnNzmBGSzoHUKGb3D+qzMQlzSI1qa0EVkTEOEkbAlcCY4FZwIERsUiSgLOBfYDlwOER8ZdOL5GZmfUKzRLMlwvd02qG1fY38/6IWFDoPwm4PSJOl3RS7j8R2BvYKr/GA+fndzMza0ENE0xETK5onhOA3XP3ZNINBCfm8osjIoC7JQ2TNDIinqkoDjMzq1DDBCNpSrMRI+KjJaYfwK2SAvhRRFwAjCgkjbnAiNw9CmgrjDs7l62WYCRNAiYBbLrppiVCMDOzntCsiWwX0gH/cuAe0rWXztotIuZIegNwm6THigMjInLyKS0nqQsAxo0b16lxzcys+zRLMG8EPggcAnwSuAm4PCIeKTvxiJiT3+dJug7YCXi2velL0khgXq4+BxhTGH10LjMzsxbU8DbliFgZEbdExERgZ2AmMFXSF8tMWNL6koa0d5Me2nwYmAJMzNUmAjfk7inAYflXA3YGlvj6i5lZ62p2BoOkQcC+pLOYscA5wHUlpz0CuC7dfcwA4LKIuEXSn4GrJB0JPAUcmOvfTLpFeSbpNuUjOrUkZmbWqzS7yH8xsB3pwH9aRDzcmQlHxJPADnXKFwJ71ikP4OjOzMPMzHqvZmcwnwKeB44FjslnIpAu9kdEDK04NjMza2HNnoPxT/mbmdkacxIxM7NKOMGYmVklnGDMzKwSTjBmZlaJDhOMpI9JekLSEkn/lLRU0j+7IzgzM2tdTR+0zL4H7BcRM6oOxszM+o4yTWTPOrmYmVlnlTmDmSbpSuB64KX2woi4tqqgzMys9ZVJMENJvw32oUJZAE4wZmbWUIcJJiL8o5NmZtZpzX7s8isR8T1J55LOWFYTEcdUGpmZmbW0Zmcw7Rf2p3VHIGZm1rc0+7HLX+b3yd0XjpmZ9RV+kt/MzCrhBGNmZpVwgjEzs0p0eJuypOHA54CxxfoR8ZnqwjIzs1ZX5kHLG4A/AL8BVlYbjpmZ9RVlEszgiDix8kjMzKxPKXMN5kZJ+1QeiZmZ9SnNnuRfSnqCX8Apkl4CXsn9ERFDuydEMzNrRQ3PYCJiSEQMze/9ImK9Qn/p5CKpv6T7JN2Y+zeXdI+kmZKulLROLh+U+2fm4WO7vHRmZtZjyvyj5b9J2qDQP0zS/p2Yx7Gs+tkZgO8CZ0bElsAi4MhcfiSwKJefmeuZmVmLKnMN5usRsaS9JyIWA18vM3FJo4F9gQtzv4A9gKtzlcnA/rl7Qu4nD98z1zczsxZUJsHUq1Pm7jOAs4CvAK/m/o2AxRGxIvfPBkbl7lFAG0AeviTXNzOzFlQmwUyTdIakLfLrDGB6RyNJ+ggwLyI6rNsZkiZJmiZp2vz589fmpM3MbC0qk2C+BLwMXJlfLwFHlxjv3cBHJc0CriA1jZ0NDJPUfgY0GpiTu+cAYwDy8A2AhbUTjYgLImJcRIwbPnx4iTDMzKwndJhgIuL5iDip/aAeESdHxPMlxjs5IkZHxFjgYOC3EXEocAdwQK42kfRLAQBTcj95+G8j4jV/dGZmZq2hzG+R3UH9f7TcYw3neSJwhaRvAfcBF+Xyi4CfS5oJPEdKSmZm1qLKXKw/odC9LvBxYEWDunVFxFRgau5+EtipTp0XgU90ZrpmZtZ7dZhg6lykv1PSvRXFY2ZmfUSZJrINC739gHeSLsCbmZk1VKaJbDqrfpNsBfB3Vj19b2ZmVleZJrLNuyMQMzPrW0o9kS9pO2Bb0kV+ACLi4qqCMjOz1lfmGszXgd1JCeZmYG/gj4ATjJmZNVTmSf4DgD2BuRFxBLADvshvZmYdKJNgXoiIV4EVkoYC88g/6WJmZtZImWsw0yQNA35MuqNsGfCnKoMyM7PW1zTB5P9j+a/8HzA/lHQLMDQiHuyO4MzMrHU1TTAREZJuBt6W+2d1R1BmZtb6ylyD+Yukd1UeiZmZ9SllrsGMBw6V9BTwPOmJ/oiI7SuNzMzMWlqZBPPhyqMwM7M+p8xPxTwlaUdgN9Jvkt0ZEX+pPDIzM2tpHV6DkfQ1YDKwEbAx8FNJX606MDMza21lmsgOBXbIfwiGpNOB+4FvVRiXmZm1uDJ3kT1N4UcugUHAnGrCMTOzvqLhGYykc0nXXJYAj0i6Lfd/EPA/WpqZWVPNmsim5ffpwHWF8qmVRWNmZn1GwwQTEZMl9QcujohDuzEmMzPrA5peg4mIlcBmktbppnjMzKyPKHMX2ZPAnZKmkJ7kByAizqgsKjMza3llEszf8qsfMKTacMzMrK8o8yT/aQCSBkfE8rITlrQu8HvSbc0DgKsj4uuSNgeuID24OR34dES8LGkQ6W+Y3wksBA7yrzebmbWuMk/y7yLpUeCx3L+DpB+UmPZLwB4RsQPwdmAvSTsD3wXOjIgtgUXAkbn+kcCiXH5mrmdmZi2qzIOWZ5F+8HIhQEQ8ALy3o5EiWZZ7B+ZXAHsAV+fyycD+uXtC7icP3zP/4ZmZmbWgMgmGiGirKVpZZjxJ/SXdD8wDbiNdy1kcEStyldnAqNw9CmjL81tBesBzozrTnCRpmqRp8+fPLxOGmZn1gDIJpk3SrkBIGijpBGBGmYlHxMqIeDswGtgJ2GaNI101zQsiYlxEjBs+fHhXJ2dmZhUpk2COAo4mnWHMIV1PObozM4mIxcAdwC7AMEntNxeMZtXvms0BxgDk4RuQm+XMzKz1lEkwiohDI2JERLwhIj4VER0e+CUNlzQsd69H+g2zGaREc0CuNhG4IXdPyf3k4b+NiCi/KGZm1puUeQ7mTkmzgCuBa/LZSBkjgfafm+kHXBURN+Y70q6Q9C3gPuCiXP8i4OeSZgLPAQeXXwwzM+ttyjwHs7WknUgH/P9sTxARcUkH4z0IvKNO+ZOk6zG15S8CnygbuJmZ9W5l7yK7NyL+g5QYnmPV7cRmZmZ1lXnQcqikiZJ+BdwFPEOdMxAzM7OiMtdgHgCuB74REX+qNhwzM+srmiaYfIH+2og4vpviMTOzPqLM/8Hs2k2xmJlZH1Kmiez+/F8wv2D1/4O5trKozMys5ZVJMOuSnqjfo1AWgBOMmZk1VOY5mCO6IxAzM+tbytymvLWk2yU9nPu3l/TV6kMzM7NWVuZByx8DJwOvwP89oe+fcTEzs6bKJJjBEXFvTdmKujXNzMyyMglmgaQtSBf2kXQA6Wl+MzOzhsrcRXY0cAGwjaQ5wN+BT1UalZmZtbwyd5E9CXxA0vpAv4hYWn1YZmbW6jpMMJIGAR8HxgIDJAEQEd+oNDIzM2tpZZrIbgCWANOBl6oNx8zM+ooyCWZ0ROxVeSRmZtanlLmL7C5Jb6s8EjMz61MansFIeoh0a/IA4AhJT5KayARERGzfPSGamVkratZE9pFui8LMzPqchk1kEfFURDwFjASeK/QvAt7YXQGamVlrKnMN5nxgWaF/WS4zMzNrqEyCUUREe09EvEq5u8/MzOxfWJkE86SkYyQNzK9jgSc7GknSGEl3SHpU0iN5PCRtKOk2SU/k99fnckk6R9JMSQ9K2rFri2ZmZj2pTII5CtgVmAPMBsYDk0qMtwI4PiK2BXYGjpa0LXAScHtEbAXcnvsB9ga2yq9JuBnOzKyllfktsnmswf+/RMQz5F9djoilkmYAo4AJwO652mRgKnBiLr84N8fdLWmYpJF5OmZm1mKaPQfzlYj4nqRzyT/VXxQRx5SdiaSxwDuAe4ARhaQxFxiRu0cBbYXRZuey1RKMpEnkM6hNN920bAhmZtbNmp3BzMjv07oyA0mvA64BjouIf7b/WCakpzUlvSZ5NRMRF5D+PoBx48Z1alwzM+s+DRNMRPwyv08GkDQ4IpZ3ZuKSBpKSy6URcW0ufra96UvSSGBeLp8DjCmMPjqXmZlZC+rwIr+kXSQ9CjyW+3eQ9IMS4wm4CJgREWcUBk0BJubuiaRfa24vPyzfTbYzsMTXX8zMWleZ51nOAj5MSgBExAOS3ltivHcDnwYeknR/LjsFOB24StKRwFPAgXnYzcA+wExgOXBEuUUwM7PeqNQDkxHRVrx2AqwsMc4fST+MWc+edeoH6e+ZzcysDyiTYNok7QpEvqZyLKtuADAzM6ur7IOWR5NuGX4aeDs+0zAzsw6UedByAXBoN8RiZmZ9SJm7yN4k6ZeS5kuaJ+kGSW/qjuDMzKx1lWkiuwy4ivS/MJsAvwAurzIoMzNrfWUSzOCI+HlErMivS4B1qw7MzMxaW5m7yH4l6STgCtJvkh0E3CxpQ4CIeK7C+MzMrEWVSTDtD0J+vqb8YFLC8fUYMzN7jTJ3kW3eHYGY2b+esSfd1NMhtKxZp+/b0yF0qOE1GEnvkvTGQv9h+Q6yc9qbx8zMzBppdpH/R8DLAPm3x04HLgaWkH8u38zMrJFmTWT9CxfwDwIuiIhrgGsKP15pZmZWV7MzmP6S2hPQnsBvC8NK/UimmZn962qWKC4HfidpAfAC8AcASVuSmsnMzMwaavaPlt+WdDvpCf5b88/pQzrr+VJ3BGdmZq2raVNXRNxdp+yv1YVjZmZ9RZmfijEzM+s0JxgzM6uEE4yZmVXCCcbMzCrhBGNmZpVwgjEzs0o4wZiZWSUqSzCSfiJpnqSHC2UbSrpN0hP5/fW5XPlXmmdKelDSjlXFZWZm3aPKM5ifAXvVlJ0E3B4RWwG3536AvYGt8msScH6FcZmZWTeoLMFExO+B2r9TngBMzt2Tgf0L5RdHcjcwTNLIqmIzM7Pqdfc1mBER8UzunguMyN2jgLZCvdm57DUkTZI0TdK0+fPnVxepmZl1SY9d5M8/nhkdVnzteBdExLiIGDd8+PAKIjMzs7WhuxPMs+1NX/l9Xi6fA4wp1Budy8zMrEV1d4KZAkzM3ROBGwrlh+W7yXYGlhSa0szMrAVV9s+Uki4Hdgc2ljQb+DpwOnCVpCOBp4ADc/WbgX2AmcBy4Iiq4jIzs+5RWYKJiEMaDNqzTt0Ajq4qFjMz635+kt/MzCrhBGNmZpVwgjEzs0o4wZiZWSWcYMzMrBJOMGZmVgknGDMzq4QTjJmZVcIJxszMKuEEY2ZmlXCCMTOzSjjBmJlZJZxgzMysEk4wZmZWCScYMzOrhBOMmZlVwgnGzMwq4QRjZmaVcIIxM7NKOMGYmVklnGDMzKwSTjBmZlYJJxgzM6uEE4yZmVWiVyUYSXtJelzSTEkn9XQ8Zma25npNgpHUH/gfYG9gW+AQSdv2bFRmZramek2CAXYCZkbEkxHxMnAFMKGHYzIzszU0oKcDKBgFtBX6ZwPjaytJmgRMyr3LJD3eDbH1hI2BBT0dRD36bk9HYNYteu0+CF3eDzdbS2E01ZsSTCkRcQFwQU/HUTVJ0yJiXE/HYfavyvtg1/WmJrI5wJhC/+hcZmZmLag3JZg/A1tJ2lzSOsDBwJQejsnMzNZQr2kii4gVkr4I/BroD/wkIh7p4bB6Up9vBjTr5bwPdpEioqdjMDOzPqg3NZGZmVkf4gRjZmaVcIJZCyT9p6RHJD0o6X5J4yUdJ2nwWpzHLEkbd2H8wyWdt7biMWsVkk6VdEInyu/qnsj6vl5zkb9VSdoF+AiwY0S8lJPAOsCVwCXA8h6Kq39ErOyJeZu1sojYtadj6Ct8BtN1I4EFEfESQEQsAA4ANgHukHQHgKTzJU3LZzqntY+cz0xOk/QXSQ9J2iaXbyTp1lz/QkCFca6XND0Pm1QoXybpvyU9AOwi6QhJf5V0L/DublgXZr1CblX4q6Q/Am/u5LjL8vvukqZKulrSY5IulaQ87J2Sfpf3w19LGpnLj5H0aG7NuGKtL1iriQi/uvACXgfcD/wV+AHwvlw+C9i4UG/D/N4fmApsX6j3pdz9BeDC3H0O8LXcvS8Q7dMrTGs94GFgo9wfwIG5eyTwD2A46YzqTuC8nl5ffvlV9Qt4J/AQMBgYCswETqhT79QG5cvy++7AEtJD3/2APwG7AQOBu4Dhud5BpMcqAJ4GBuXuYT29Lnr65SayLoqIZZLeCbwHeD9wZYO/Gjgwn20MIB38twUezMOuze/TgY/l7ve2d0fETZIWFaZ1jKR/y91jgK2AhcBK4JpcPh6YGhHzASRdCWzdlWU1axHvAa6LiOUAkrrywPa9ETE7T+d+YCywGNgOuC2f0PQHnsn1HwQulXQ9cH0X5tsnOMGsBZGudUwFpkp6CJhYHC5pc+AE4F0RsUjSz4B1C1Veyu8r6eAzkbQ78AFgl4hYLmlqYVovhq+7mK1NLxW62/dPAY9ExC516u9L+nK4H/Cfkt4WESuqD7N38jWYLpL0ZklbFYreDjwFLAWG5LKhwPPAEkkjSP9505HfA5/M89gbeH0u3wBYlJPLNsDODca/B3hfvpYzEPhE+aUya2m/B/aXtJ6kIaSD/dr0ODA83+CDpIGS3iqpHzAmIu4ATiTtq69by/NuKT6D6brXAedKGgasILX3TgIOAW6R9HREvF/SfcBjpL8kuLPEdE8DLpf0CKm99x+5/BbgKEkzSBv63fVGjohnJJ1KajdeTLpOZNbnRcRfcpPwA8A80u8cNvJVSccVxh1dYvovSzoAOEfSBqTj6Fmk67CX5DIB50TE4jVdjr7APxVjZmaVcBOZmZlVwgnGzMwq4QRjZmaVcIIxM7NKOMGYmVklnGDMzKwSTjBmZlaJ/wXt/mh/lFiATwAAAABJRU5ErkJggg==\n"
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "jupyter_overhead = memory_usage(print)\n",
    "avg_jupyter_overhead = sum(jupyter_overhead) / len(jupyter_overhead)\n",
    "print(avg_jupyter_overhead)\n",
    "for i, item in enumerate(results_memory):\n",
    "    results_memory[i] = item - avg_jupyter_overhead\n",
    "\n",
    "fig, ax2 = plt.subplots()\n",
    "\n",
    "rects1 = ax2.bar(x, results_memory, width, label='Tomatenvirus & Chrysovirus')\n",
    "\n",
    "ax2.set_ylabel('Speicherverbrauch in MB')\n",
    "ax2.set_title('Speicherverbrauch von verschiedenen LCS Implementationen')\n",
    "ax2.set_xticks(x, labels)\n",
    "ax.legend()\n",
    "\n",
    "\n",
    "\n",
    "plt.show()"
   ],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "outputs": [],
   "source": [],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}