network_static.ipynb (80341B)
1 { 2 "cells": [ 3 { 4 "cell_type": "code", 5 "execution_count": 1, 6 "id": "4c9bab78", 7 "metadata": {}, 8 "outputs": [], 9 "source": [ 10 "import itertools\n", 11 "import random\n", 12 "import collections\n", 13 "import os\n", 14 "from datetime import datetime\n", 15 "\n", 16 "from scipy.optimize import curve_fit\n", 17 "import numpy as np\n", 18 "import pandas as pd\n", 19 "\n", 20 "from community import community_louvain\n", 21 "import networkx as nx\n", 22 "\n", 23 "import matplotlib.pyplot as plt\n", 24 "import matplotlib.cm as cm\n", 25 "import matplotlib \n", 26 "\n", 27 "matplotlib.rcParams.update({'font.size': 15})\n", 28 "%matplotlib inline\n", 29 "\n", 30 "import plotly\n", 31 "import plotly.graph_objs as pgo\n", 32 "\n", 33 "from nxplot_3d import get_trace_nodes_3D_legend\n", 34 "from nxplot_3d import get_trace_edges_3D\n", 35 "from nxplot_3d import plot_3D\n", 36 "\n", 37 "from export_import import export_data\n", 38 "from export_import import read_data\n", 39 "\n", 40 "from csv_to_network import get_graph\n", 41 "\n", 42 "path = \"../data/static/\"" 43 ] 44 }, 45 { 46 "cell_type": "markdown", 47 "id": "0dd9a67b", 48 "metadata": {}, 49 "source": [ 50 "# Fetch Data & make Network" 51 ] 52 }, 53 { 54 "cell_type": "code", 55 "execution_count": 19, 56 "id": "9ec366d6", 57 "metadata": {}, 58 "outputs": [], 59 "source": [ 60 "req = pd.read_csv(path + 'my_req.csv', sep='|', keep_default_na=False)\n", 61 "DG = get_graph(req, cut=0, direction=\"\")\n", 62 "list_node_deg = sorted(list(dict(DG.degree()).items()), key=lambda x: x[1], reverse=True)" 63 ] 64 }, 65 { 66 "cell_type": "markdown", 67 "id": "95be6588", 68 "metadata": {}, 69 "source": [ 70 "# Degree Distribution" 71 ] 72 }, 73 { 74 "cell_type": "code", 75 "execution_count": 57, 76 "id": "cbbcddda", 77 "metadata": {}, 78 "outputs": [ 79 { 80 "name": "stderr", 81 "output_type": "stream", 82 "text": [ 83 "/tmp/ipykernel_30438/714574286.py:48: RuntimeWarning: divide by zero encountered in log\n", 84 " ax[0].plot(k, np.exp(linear(np.log(k), popt[0]+1, popt[1])),\\\n", 85 "/tmp/ipykernel_30438/714574286.py:48: RuntimeWarning: divide by zero encountered in log\n", 86 " ax[0].plot(k, np.exp(linear(np.log(k), popt[0]+1, popt[1])),\\\n", 87 "/tmp/ipykernel_30438/714574286.py:48: RuntimeWarning: divide by zero encountered in log\n", 88 " ax[0].plot(k, np.exp(linear(np.log(k), popt[0]+1, popt[1])),\\\n" 89 ] 90 }, 91 { 92 "data": { 93 "image/png": "iVBORw0KGgoAAAANSUhEUgAAAsgAAAFgCAYAAACmDI9oAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABVeUlEQVR4nO3deXxU1f3/8deZTDbCDklkIMoOCYtRIiguqKggigoqgktBUbR16Wartmpt1a9W609rtSpu1B20VkEjVK2AooBssiPIYiCSRHbInpzfH0NiEpKQZZI7M/f9fDzy0Llz5+ZzHfjMxzOfc46x1iIiIiIiIn4epwMQEREREQkmKpBFRERERCpQgSwiIiIiUoEKZBERERGRClQgi4iIiIhUoAJZRERERKQCFcgiIiIiIhWoQBYRERERqcDrdAANZYyJA/4JFAJzrbWvOxySiIgrKR+LSLgJqhFkY8xLxphsY8zqKsdHGmM2GGM2GWPuPHx4LPCOtfYG4KJmD1ZEJIwpH4uImwVVgQxMA0ZWPGCMiQCeBs4HUoAJxpgUoAuQcfi0kmaMUUTEDaahfCwiLhVULRbW2vnGmK5VDg8GNllrNwMYY94CLga240/KK6il0DfGTAGmAMTFxQ3q27dv4AMXEamHpUuX/mitjXc6jtooH4uIG9SUj4OqQK5BZ34amQB/Ih4CPAk8ZYy5AJhV04uttVOBqQBpaWl2yZIlTRiqiMjRGWO2OR1DAykfi0hYqSkfh0KBXC1r7SHgWqfjEBFxO+VjEQk3wdaDXJ0dQFKFx10OH6szY8xoY8zUffv2BTQwERGXUT4WEVcIhQL5a6CXMaabMSYKGA/MrM8FrLWzrLVT2rRp0yQBioi4hPKxiLhCULVYGGPeBM4EOhpjtgN/sta+aIy5BZgDRAAvWWvXOBimSNgrKipi+/bt5OfnOx1KSIuJiaFLly5ERkY6HUq9KR+LNI7yaHCpbz4OqgLZWjuhhuPpQHpDr2uMGQ2M7tmzZ0MvIeIq27dvp1WrVnTt2hVjjNPhhCRrLbt27WL79u1069bN6XDqTflYpHGUR4NHQ/JxKLRYNJq+0hOpn/z8fDp06KCk3gjGGDp06KDRoyqUj8UtlEeDR0PysSsKZBGpPyX1xtN/QxF3Uw4IHvV9L1xRIGvWtIhIcFA+FpFQ4IoCWV/pibjPmWeeiTaiCD7KxyKhwe051BUFsoiIiIhIXQXVKhYiEroys3LIyMwmyZeAL/GIbe3rbevWrVx44YWsXr0agL/97W8cPHiQuXPnMmTIED777DP27t3Liy++yOmnn05eXh7XXnst33zzDX379iUvL6/RMYiINKdA5lHl0MZxRYGsZYVEmlZmVg7TZqRTXFKKN8LDpHGjAlIk16S4uJjFixeTnp7On//8Zz755BOeeeYZWrRowbp161i5ciUnnnhik/1+aTjlY5HqNWceVQ49Ole0WKjnTaRpZWRmU1xSSlyLGIpLSsnIzG7S3zd27FgABg0axNatWwGYP38+V199NQADBw5k4MCBTRqDNIzysUj1mjOPKocenSsKZBFpWkm+BLwRHg7l5uON8JDkS2j0Nb1eL6WlpeWPK65fGR0dDUBERATFxcWN/l0iIk4LdB5VDm0cFcgi0mi+xHgmjRvFiGGDA/a1YGJiItnZ2ezatYuCggI++OCDWs8/44wzeOONNwBYvXo1K1eubHQMIiLNJdB5VDm0cVzRgywiTc+XGB/QfrnIyEjuvfdeBg8eTOfOnenbt2+t5//85z/n2muvJTk5meTkZAYNGhSwWEREmkMg86hyaOMYa63TMTS5CpNCbti4caPT4YgEvXXr1pGcnOx0GGGhuv+Wxpil1to0h0JyVEPzcaBXSRFpasqjwac++dgVLRaaFCIiEhwako/LZvfPmbeYaTPSyczKacIIRURcUiA32A8/wOOPOx2FiIirZWRmE7l/P2krlzTLKikiIiqQa9O6NbzwAuzd63QkIiKuleRLoDgujuXJx9M69yBJiR2cDklEwpwK5NrExcGqVdC2rdORiIi4li8xnokTRnPOOacyZeWX+P50D7hg/oyIOEcF8tF4PDBhAnz7rdORiIi4li8xniEn9KPFC8/7By7uuENFsog0GVcUyMaY0caYqfv27WvYBZKT4aGHAhuUiIgLNToft2wJ6enw2WewaVNggxMROcwVBXKjV7G47TaYNQu2bAlsYCIiLhOQVYXat4dFi6BXL1i+vF4vzczKYdHyNVoJQ0Rq5YoCudHatoV774Xvv3c6EhERAX/72+7dMHo0HN7962i0XJyI1JUK5Lq67TY44wzIzXU6EhFphIyMDM466yxSUlLo168ff//732s8d/bs2fTp04eePXvy8MMPlx//+9//Tv/+/enXrx9PPPFEM0Qt1WrfHmbPhl//Go6yjS74l4srLiklrkWMlosTaaTrrruOhIQE+vfv32zXas6crAK5Ph57zD8xRERCltfr5bHHHmPt2rUsXLiQp59+mrVr1x5xXklJCTfffDMfffQRa9eu5c0332Tt2rWsXr2a559/nsWLF/PNN9/wwQcfsEm9sM7p3x9mzoQZM456apIvAW+Eh0O5+XgjPCT5EpohQJHwNGnSJGbPnt1s12runKwCuT6uuQZef92/gYiINKmNGzfStWvX8kRXVFREamoqGRkZjbpup06dOPHEEwFo1aoVycnJ7Nix44jzFi9eTM+ePenevTtRUVGMHz+e999/n3Xr1jFkyBBatGiB1+tl2LBhvPvuu42KSRppyBB45RXYuhWWLavxNF9iPJPGjWLEsMFMGjdKW1ZL2GuqPApwxhln0L59+0Zfp67Xau6crAK5PhIT/UXy3/7mdCQiYa9Xr15MmTKFOXPmAPDUU09x0UUXkZSUVONrTj/9dFJTU4/4+eSTT6o9f+vWrSxfvpwhQ4Yc8dyOHTsq/a4uXbqwY8cO+vfvz+eff86uXbvIzc0lPT09IB82EgCrV8MFF8C6dTWeUrZcnIpjcYOG5NFAqm9Ork1z52Rvo68QAowxo4HRPXv2bPzFfv97WLKk8dcRCSX33Qd//vNPj8v+DqSl/XTsT3/yn+fz/fQty4knwtKlMGUKPP/8T+fu2OE/7yj69+/PJ598wu7du3nxxRdZtGhRred//vnndbsf4ODBg1x66aU88cQTtG7dus6vS05O5o477uC8884jLi6O1NRUIiIi6vx6twtoPq7qwgvhr3+FESPg88/huOMC/ztEGipE8ug555zDzp07jzj+4IMPcvHFFx/191VUn5zcUE2Vk10xghyQZYXKdO4Mo0b5/7CKuMV99/k3ZSj7GTTI/1Px2H33+c/NzPzpWNnfk6lTK59bh6QO0Lt3bzZs2MB9993H7bffzubNm/nb4W9wbr75ZnKrTJqt62hFUVERl156KVdddRVjx46t9nd37ty50ijE9u3b6dy5MwCTJ09m6dKlzJ8/n3bt2tG7d+863Y8EOB9X52c/888V+e67Rl9KS8JJQIVIHv3kk09YvXr1ET/1LY4hsCPIzZ2TXTGCHHA7dsB558HGjf5Z1CLSJHr06MGyZcvYt28fTzzxBK+99hrHH3884B8BbtGiRaXz6zJaYa1l8uTJJCcn85vf/KbG80466SQ2btzIli1b6Ny5M2+99RZvHF5OLDs7m4SEBL7//nveffddFi5c2Ii7lIC7+WZ/AfHsszB+vH+pznoqWxKuuKQUb4RHPcsSsuqbRwMpkCPIzZ2TXTGCHHBdu8LFF0Mty0OJSONFRkbSunVrHn74YTweD6tWrWLgwIHs378fY0yDrrlgwQJeffVV/ve//5WPZqSnp5c/P2rUKDIzM/F6vTz11FOMGDGC5ORkxo0bR79+/QC49NJLSUlJYfTo0Tz99NO0bUABJs1g/Xp/T/KhQ/V+qZaEk3DRFHkUYMKECZxyyils2LCBLl268OKLLwb8WmX5GGj2nGysi/ayT0tLs0sC1T+8aROcfLJ/1nTLloG5pkiQWLduHcnJyU6HAcCxxx7Ltm3bMMbw4osvsn79erxeL23atOHOO+90Oryjqu6/pTFmqbU2rYaXuEJA83FNSkvhuutg507/UnBRUXV+qUaQpbGUR4NPffKxWiwaqmdP+OILiItzOhKRsLV161aOO+648lGOyZMnOxyRhBSPB154AR5/3F8s10PZknAZmdkk+RJUHEvIUh5tGLVYNEafPv42i4MHnY5EJCx17dq1WWZBSxjzeuF3v/OPIt99t783uY60JJyEA+XRhlGB3BjGwJdf+ieCiIhI8OrQAT7+2L/ChYtaC0WkYVQgN9bdd/u3oM7LczoSERGpSatW8NFHkJ4Ojz7a6MtpCTiR8KYe5MYaONA/We+TT2D0aKejEQkYa22jZjiL/7+hBJH27eG//230GsmawCd1pTwaPOqbj10xgmyMGW2Mmbpv376m+QUzZqg4lrASExPDrl27VOA1grWWXbt2ERMT43QoQaXJ8/HR+Hxw+unwyCNweA3V+tIScFIXyqPBoyH52BUjyNbaWcCstLS0G5rkF0RGwrRp/p7kiROb5FeINKcuXbqwfft2cnL09XFjxMTE0KVLF6fDCCpNno/ratQoOOccaNPGv1ZyPST5EvBGeDiUm483wkOSL6GJgpRQpjwaXOqbj11RIDeL3r3hqqvgyiv9BbNICIuMjKRbt25OhyHSdPr3h/ff93/7N2cOnHBCnV+qJeCkLpRHQ5srWiyaxdCh0L07vP6605GIiEhdDBni70nu16/eK1toCTiR8KYCOZDuvRe+/97pKEREpK5SU6GkBAYPhnXrGn05rW4hEh7UYhFIw4b5f4qK1GYhIhIqYmPhlltgxAj4/HM47rgGXUarW4iED40gB9ratZCWVu9tTUVExEETJ8Ltt8MllzQ4f2t1C5HwoQI50JKTISYG3n3X6UhERKQ+brsNZs3yr0jUgM2ftLqFSPhQgRxoxsA998ADD2g7UxGRUNOlC7z3Hpx7Lhw6VK+Xlq1uMWLYYLVXiIQ49SA3hQsugE2boKDAP5osIiKh4+KL/UvAXXopzJwJUVF1fqkvMb5SYZyZlaPl4ERCkEaQm4Ix8KtfwebNGkUWEQk1Hg+88AK0aAFPPdXgy5RN2pszbzHTZqRrZQuREKICualYC1dc4V+AXkREQovXC2++6e9L3rmzQYMdmrQnErpUIDcVY+APf4D779cosohIKIqO9hfKN94Id9xR71yuSXsioUsFclMaNw5+/BE++8zpSEREpKFefhnS0+Hhh+v1Mk3aEwldIVsgG2O6G2NeNMa843QsNYqI8C/3NmSI05GIiDSZkMjHjdG+vX9L6hkz4Icf6vXSqltSa6c9kdDgSIFsjHnJGJNtjFld5fhIY8wGY8wmY8ydtV3DWrvZWju5KeMMSCLr1w9WrIDlywMWl4hIoIRKPnaczwdLl8Ixx/hzegNo0p5I6HBqBHkaMLLiAWNMBPA0cD6QAkwwxqQYYwYYYz6o8tPkjVwBTWRr18JddwUuOBGRwJlGkOfjoOHxwI4dMHIkfPhhvV+uSXsiocORAtlaOx/YXeXwYGDT4ZGIQuAt4GJr7Spr7YVVfpo8qwQ0kU2c6C+SFy8OXIAiIgEQCvk4qHTp4l8j+dprYd68er1Uk/ZEQkcw9SB3BjIqPN5++Fi1jDEdjDHPAicYY2ocnjXGTDHGLDHGLMnJqfsocEATWVSUfwb08883/BoiIs0nqPJx0BkyxL8E3Ntv1+tl1U3aU0+ySHAK2Z30rLW7gJvqcN5UYCpAWlpandfoKUtkAdsB6YYb/F/PiYiEmabOx0Fp+HD/z7p1/mU9+/at08sq7rRX1spXXFKKN8KjlS5EgkgwVWw7gKQKj7scPtZoxpjRxpip+/btq9frqs4+bqjMrBwWrdlI9sLFcGetc11ERIJB0OXjoLV8OZx3HmzbVu+XqidZJHgFU4H8NdDLGNPNGBMFjAdmBuLC1tpZ1topbdq0CcTl6qXiZL9/LVxNyYsvwpo1zR6HiEg9hGU+bhJXXgm33w7nngtZWfV6qXqSRYKXU8u8vQl8BfQxxmw3xky21hYDtwBzgHXADGttyFeSFUcI8ryR7JhwNfzf/zkdlogI4K583GRuuw2uvx6+/75eL6vakwyoH1kkSDjSg2ytnVDD8XQgPdC/zxgzGhjds2fPQF/6qKqOEET+8jb4xU1QWOifvCci4iA35eMm9fvfQ2kpPPkkTJ4McXF1ellZT7L6kUWCSzC1WDQZJ7/SqzpC0KlHN5gzB7whOz9SRKTBwq7Foqply+Cyy/yDIPWgfmSR4OKKAtlpR0z2Ky2FQYNgyxZnAxMRkcDxeOCFFyAmBq6+GkpK6vxS9SOLBBcVyE7weGDUKPjrX52OREREAsnr9a+RPHgw2LqvZKd+ZJHg4ooCOSiXFfr1r/2LzG/f7nQkIiLNJijzcaDFxPhXtvj2W3jggToXymXfNgLlqx9Nm5GuIlnEAa4okIOy561jR3j0UTh0qNJh7aokIuEsKPNxUznmGHjrLXj44Xq9TP3IIs7TTDEnXXcd7N8Pe/ZAu3aaxSwiEk7at4f//hdOPx06dIApU+r0MvUjizhPBbLT7r/f//Xb3/5WadTgUG4+GZnZKpBFREKZzwcffwy7d9f9JYf7kTMys4mOjiwfQdbngUjzcUWLRVD3vP3yl/DSS5CTo1EDEQl7QZ2Pm0r37pCWBn/4A3zwQZ1e4kuMJ8mXQPqnX6kXWcQBriiQg7rnrUsXGDcOHn/8iFnMGi0QkXAT1Pm4qV18sb+1bt68Op2uXmQR56jFIhjccw8cHk0p21VJRETCzJAh/iXgLr8cFiyAXr1qPb3it4pg2bvvAJlZOfqMEGkGrhhBDnqdO/tXtZgzx+lIRESkKQ0fDrNn+9sujrL8W9m3ikNOSAYLi1asU6uFSDNxRYEcEj1v+/f7d17av9/pSEREmkxI5OOmduKJcOCAvy9527ZaT/UlxtO2dSswRq0WIs3IFQVySPS89ewJI0bA00/XeIrWSBaRUBcS+bg5tG0LEyfCuedCVlatp9bUaiEiTccVBXLI+OMf4ZlnoLj4iKfK1kjWbGYRkTBx221w1VVwxRW1tluo1UKk+alADibJyfDNN+A9cu6kZjOLiIShe++F11/3F8j5+TWeplYLkealAjnYtG4Nt9wCeXmVDmuNZBGRMGSMf6L2yy/DmDFQWFjjqWq1EGk+KpCDTUQEbN8OL7xQ6bDWSBYRCWMTJ0JMjH+ydklJtaeo1UKk+biiQA65WdN33w2PPAIFBZUO+xLjGXJCPxXHIhKyQi4fNxev179G8u7d8OqrNZ6mVguR5uGKAjnkZk2npcFpp8HKlU5HIiISUCGXj5tTTAzMmgU/+xlkZtY4cU+tFiJNzxUFckh680046aSjLiQvIiJhJDbW35c8fjw8/HC1p6jVQqTpqUAOZvfcA6+84nQUIiLSnIyBt97yz0V55plqT6naalFQWMS8hStUJIsEiArkYDZ8ODz4YI0TNkREJEz5fPDxx/Dcc1BDv3ZZq8X+A4c4mJvHuk3bNJIsEiAqkIPZsGGQmAjTp9d6mnbYExEJQ927w7Jl0LKlf438KspaLfr2PI64FrG0aRWnSXsiAaICOZgZAw884P9nDbTDnohIGPN4YNMm/5bU8+Yd8bQvMZ5hJ6cSExXJ/gOHKC4uIjo60oFARcKLKwrkkF5WaNgw/2SNnOoLX+2wJyKhJKTzsVP69PFP3L78cli69IinfYnxjBp+CtaWgvGQ/ulXGiwRaSRXFMghv6zQ3Ln+0YNqVrTQDnsiEkpCPh87ZfhwmDoV/v3vap8uKCjC642kTas4TdgTCQCv0wFIHZx5pr/NYtYsuOiiSk+V9aBlZGaT5EvQJiIiIuHqkkv8P0uXQseOcNxx5U9VN2Fvy/eZ2nlVpIFcMYIc8ozx7653//3VjiJrhz0RERdZtMj/rWJWVvmhqhP2WsRGcygvn1Xrv3MwUJHQpRHkUDFmjL9QtrbWSXsiIhLmfvEL/7yUESP8LXht2wI/TdjbuCWDnF17AFi0bC0D+vbQAIpIPWkEOVR4PHDxxTB7tnbXExFxu3vv9bdbZGZWOuxLjGdIagox0dF0bN8OjNHkbZEGUIEcan71K/+IwVFobWQRkTBmDNx3H/TuDY8/DoWF5U8NSO5BXGwMRUXFgGXvvgP6LBCpJxXIoSQiAv7wB38vci20NrKIiIvMnw9XX12+62pZP/KQE5LBwqIV6/RZIFJPKpBDzVVXwdat1a6FWUZrI4uIuITX618jedcuuPHG8hY8X2I8bVu3AmOIivJqwp5IPalADjWRkbBgAZx4Yo2naG1kEREXiYmB99+H/v0rzVHx535Lzq495BcUsGjZWo0ii9SRKwrksNu5qVMnePVVWLas2qfLvl4bMWyw1sAUkaASdvk4WLRs6Z+jsmSJvyeZyhP22rZpRUFRkUaRRerIFQVyWO7ctH+/f4JGDbQ2sogEo7DMx8GkSxd46il45hnAP2EvOiqSvfsOaBRZpB5cUSCHpcmT/SMFK1Y4HYmIiAQLnw8+/hgefBCmT9coskgDqUAOVbGxcPvt8NFHdTpdy76JiLhE9+7+NfMHDAAqjyLn5Rcwf9EKVqz91uEgRYKbdtILZb/+dZ121Stb9q24pBRvhEd9ySIi4a5/f/8/b7kF3+WXMyQ1hbkLl1NUVExefiHvps8loUM7fRaI1EAjyKHMGP/6l7/4Ra2nadk3ERGXGjMGLr+cEwsPEumNoNRaPMZQVFyiVguRWqhADnWpqfD227BxY42naNk3ERGXGj4cpk4lfuI1jEtLJjYmmlJbSkFhkVotRGqhFotQ17o13HIL/N//wcsvV3tK2bJvGZnZJPkS9JWaiIibXHIJdOlC/0GDOMN6mbtohVotRI5CI8jh4Lbb/JP2KiwQX5WWfRMRcbG0NPjhB8761c9pm3dIrRYiR6ECORy0awf//CccPFin07WihYiIC/l8RI0dw23TX6SdLaHUllJYVKS1kUWqoQI5XOTnQ+/ekJFR62llK1rMmbeYaTPSlRRFRNzk3nuJHXEeUz6dSWxMjNZGFqmBCuRwERMDV10FjzxS62la0UJExMWMgccfx/7jSWIiPBzYtUdrI4tUQwVyOLn9dnj9ddi5s8ZTtKKFiIjLeTwcc/wAxm38hqv/8zoR1pKXX8j09z9RkSxymArkcHLMMfDii+CteXGSshUtRgwbrA1DRERcrNVdd9Ay7xCXzpyOLS0lv7CId9PnqvVOBBXI4WfMGMjKgl27ajxFK1qIiIjvuCQOvvEGnbMz6bdhDREej1a1EDkspAtkY8wlxpjnjTHTjTHnOR1P0PjHP+Cxx+p0qla0EJHGUi4OXQNPOpHd773PltRBtNy3RxuIiBzmWIFsjHnJGJNtjFld5fhIY8wGY8wmY8ydtV3DWvuetfYG4CbgiqaMN6TccQc89xzs3l3raVrRQkSUi+X4QcdzxkkDufH15zj16y/UjyyCsyPI04CRFQ8YYyKAp4HzgRRggjEmxRgzwBjzQZWfirPL7j78OgHo1g0uusi/NnIttKKFiKBcLMCAfr14a9IvGD7/v6SuXEJ+YRFvf/CZBk7EtRzbatpaO98Y07XK4cHAJmvtZgBjzFvAxdbah4ALq17DGGOAh4GPrLXLqvs9xpgpwBSAY489NnA3EOwefti/u14ttKKFiDRXLj58njvzcQjwJcZz1jWX81JuLpe//Sqr+g4kz8Cq9d9pvoq4kmMFcg06AxV3utgODKnl/FuBc4A2xpie1tpnq55grZ0KTAVIS0ureS/mcJOYCGvW+H/Gjav2lLIVLTIys8uL40XL15DkS1BCFHG3gOdicHE+DhGpKb3ZdO5wnmjXEU+ppVPWDxw4mOt0WCKOCLYCuV6stU8CTzodR9AyBm69FS64AOLiqj3FlxiPLzG+vB+5uKQUb4RHS8CJSJ0pF4ePoScNYOWG72j77XpueG0qr3g9rOjWhdSU3k6HJtKsgm0Vix1AUoXHXQ4faxRjzGhjzNR9+/Y19lKhJSUFzjgDnq12MKcS9SOLSAVNkovBxfk4RPgS4zl10AByuvXknYvHc81LT/O/51/ThD1xnWArkL8GehljuhljooDxwMzGXtRaO8taO6VNmzaNDjDk3H03vP022Nq/zVQ/sohU0CS5GFyej0PEgOQeRHojWNV3ALPOuZBeq1dowp64jpPLvL0JfAX0McZsN8ZMttYWA7cAc4B1wAxr7RqnYgwLxx8PCxb42y1qoR32RNxJuViq8iXGM3bUmUR6I1iSOpj/nTacxPVrmf3W+yqSxTWcXMViQg3H04H0QP4uY8xoYHTPnj0DednQkZcHEybAO+9AdHSNp5X1I4uIezRnLgbl41CRmtKbTZu3s2DpKgB6bt3E8envMO3AISbdcKU+KyTsBVuLRZNw/Vd6LVtCcTH86191Ol2764lIU3F9Pg4hQ08aQItY/6DKx2ecy3ddezJ+6t9Z+MXXDkcm0vRcUSAL/l7khx6CoqJaT9PueiIiAv5vFS+/8GyiIr1gDO+NuIQtx3Zn4+Jl+myQsOeKAlmzpoFTT4WRI2FH7RPRtZqFiDQl5ePQkprSmwmXnIvXG4H1ePjg3NH8GBlD9t1/gpISp8MTaTKuKJD1ld5hzzwDXbrUmtS0moWINCXl49CTmtKbIcenVDoW9+mnZF5y6VFXSBIJVa4okKWCa66B6dNrfFqrWYiISFUV+5GLIyN5cfx1FC1fzrarfqYiWcKSCmS3mTQJHnwQSktrPMWXGM+QE/oBaLKeiIiU9yN7IyIAKIiOYepVN7L0QIE2EZGw5IoCWT1vFZx3nn/b6XffrfU0TdYTkaagfBy6UlN6c/7ZJ5c/zm0Rx+dpp/L1E8+w92+PORiZSOC5okBWz1sFxsDjj/t7kWuhyXoi0hSUj0Pb2UPTSE2uvIZ1VotWRD70ELz5pkNRiQSeKwpkqeLUU2HAANi0qcZTNFlPRESqc+6wwcRER5U/3tW+I09dcT2FN98CH33kYGQigaMC2a1mz4Yrr6xxcoUm64mISHV8ifFccdFwIjw/lRA7EzvxzLhrWeutebdWkVCiAtmtxoyBQ4fgv/+t8ZSyyXq+xHjtriciIuVSU3ozavgplY5t7ZTEy4vWsPuSsbBsmUORiQSGKwpkTQqphscDf/wjPPDAUU/VhD0RCRTl4/BRXT9ycUkps+I6UDRiJKxf71BkIo3nigJZk0JqcMUV8PTTRz1NE/ZEJFCUj8NL1X5kgBW9+/PusBEUn3MO7N3rTGAijeSKAllqEBEBPXrA1Km1nlY2YW//gUMUFxcRHR3ZTAGKiEgwK+tH9nojKh1fOCCN96b8Etq2rXXdfZFgpQLZ7aKi4KGH4MsvazzFlxjPqOGnYG0pGA/pn36lNgsREQH8/chXjTkPj8dUOr7ARvHJC6/6V07SSLKEGBXIbhcZCXfdBfffX+tpBQVFeL2RtGkVpzYLERGpJDWlNxcMH3rE8Q+37yIjqStceCHk5jZ/YCIN1KAC2RjTNsBxNClNCjmKiRPh229h27YaT9G6yCLBSflYgsXZQ9MYdnJq5YPG8Pd+Q9jdMQFuu82RuEQawtga1sGt9UXGzLTWXlTlWBdr7faARdYE0tLS7JIlS5wOIzjl50NMTK2nZGblkJGZTXR0JAUFRST5ErQ+skgDGGOWWmvTAnQt5WMJKu/Nmce8hSsqHfPaUiaePYT+QwaB1+ufAyMSBGrKxw1tsZhvjPm/ChcfDHzc0OAkCMTEwB/+AN98U+MpvsR4knwJpH/6lZZ8EwkeyscSVC4ZMYykTpW/ZSw2Hl6et4SsW38JN91U4yZVIsGiQQWytfZvQE9jzFhjzDjgRWBsQCOT5tehAzz4YK2naMk3keCifCzB6OzTBh1xrLTU8mRiDw4tWgx33KEiWYJavQpkY8w8Y8z/M8ZcBTwG3AfcBJxhrV3XBPFJc7rpJpg3D9bV/FaqF1kkOCgfSzBLTel9ZD8ykOuN4qFR4zn4zr9hwYLmD0ykjrz1PP9m4AQgDbgBSAIKgAeMMSustc8HOD5pTnFx/jaLlSshObnaU3yJ8UwaN4qMzGz1IIs4S/lYgtolI4YBHNGPfCi2BX+ZcBMjPTGcnZkJPp8D0YnUrl4FsrV2NbAaeLXsmDGmB3Ai/kQtoe6Xv/T/s6SkxkkUvsR4fInxZGblsGj5GhXKIg5QPpZQcMmIYbRuFccHnyyo1FFRFBXFR7Pnceq0vxP917/ChAnOBSlSjaO2WBhjIo0xvzLGzDbGLDbGvGeMucUY0wrAWvudtfZta+0fmj7chtGyQvX05pswZUqtp2Rm5TBtRrom64k0MWNMvDFmujFmtzEmzxjzlTFmdHXnKh9LMDp7aBoXnnPqEceLI7w8edGVFN16G3z4oQORidSsLj3I/wD+HxABfA2UAPcD24wxlzZhbAFjrZ1lrZ3Spk0bp0MJDSNHwnvvwdatNZ6iyXoizebvwCXAbOAVIA54zxgTkkNuysfuVO0ayUBmfCeeufRnFN53n/+bS5EgUZcC+XLgXmvtudbam621l+LvdXsMeN0Yc2GTRijNr107uPFGePjhGk/RZD2RZnMO8Btr7ZXW2huttQOBqcBDDsclUi+XjBhWbZG8xXcsT0y6FQoKYO3a5g9MpBp1KZAt8EWlA9YetNY+CDwOPNAUgYnDfvObGifqwU+T9YackMygAX2aMTAR1+kIVN1R42EgyRhznAPxiDRYTUXyDz/u4bU77qPkrLNg/frmD0ykiroUyLOBMTU8NwdQdRSOOnb0bwv63Xe1nrZ05QYWrVinPmSRplVa5fF2wADtHIhFpFEuGTGs2oGVpR07M/3Ucyk862zYts2ByER+UpcCeTNwrTHmQWNM1WR8OrAm8GFJUNi1C046CXburPZp9SGLNJvHD695PNkYczJQ1sBrnAxKpKGuHjuS5J5HfgHy9fEn8WHqyXw/7dVqXiXSfOpSIN8MtATuAr43xnxojHneGPM/4HpgUhPGJ07q2BGuvhoee6zap6v2IUdHR7Jo+RqNJIsE1t+BPOBK4HlgAVD2l+xRY8xfjDGXGWN6G2NUMEvImHLVJXRL6nTE8fknD+PvpjU/vv0O7N3b/IGJAMbWYatHY0wScHyVnx74C+x9wCrgG+Aba+0LTRZtI6WlpdklS6q28kmttm+HgQPh22/9BXMVmVk5ZGRmcygvj7lfLgPjISYqkknjRmltZJEaGGOWWmvTGvC6BGDA4Z+Bh/+ZAsTiny+SZ61tGchYm4rysZR58qUZbMn44YjjYz6eyYmH9tDyi8/9G1mJNIGa8nGdtpq21mZYaz+w1j5orR1nre0DtAZOAe7AXyCfAPwtkEFLEOjSBWbNghqWZPIlxpPkS2Dul8s4mJtPXl4++YVFarcQaQLW2mxr7afW2iestddZa0/C/w1fX+AKlIMlBN123bhq2y3eG34ha62XXeecC4WFDkQmblanArk61tpca+0ia+3Uw8u/nWatbRvA2CRYnHqqfxH3PXuqfTojMxtjPER4PJSUloIt1bJvIs3E+n1rrX3HWnuf0/GINER17RbW42H6ReNZRRQ7v93oUGTiVg0ukMVlZs6EJ5+s9qkkXwLRUZG0iI2hZYsYxo46U+0VIiJSL7ddN47OiZVb+UojInh/5BieeWcO+x54EOrQFioSCK4okLW1aQDcdRc89RTs33/EU2VrIo86+xRuumYMCR3aabKeiFRL+Vhqc/tNV1U7ca+gsIi9z07l4C23qkiWZuGKAllbmwZAr15w3nnw5pvVPu1LjGfICf0AmDYjnTnzFmttZBE5gvKxHM1t1407okguiI5h6lU3kvvuf2rd5VUkUFxRIEuAPPssTJlS6ylaG1lERBrrtuvGkdC+baVjuS3iePrKG3lz4UremzPPmcDENVQgS921agUffwzTptV4StW1kTVZT0REGuKuWyceUSTvb92GxSeezK6XX+GDyTc7E5i4ggpkqZ+EBPjDHyAvr9qny/qRh5yQXO1WoiIiInV1160Tadv6yDWQf2wfzxnTX2XmL+90ICpxAxXIUj+pqZCWBi++WOtpS1duYNGKdepDFhGRRrnhyoupukfkzsROvDhhMme98DRv/O5PzgQmYU0FstTfPffA4sU1Pl3WhxwV5eVQXj6r1n/XjMGJiEg48SXGc/uNV9K2deVNIr/v0pWXr7iO9dbLky/NcCg6CVcqkKX+TjoJXnmlxqV2/H3Hlpxde8gvKGDRsrUaRRYRkQbzJcbzp19PZtjJqZWObzmuOwfjWnLGY//Hqw8+7kxwEpZUIEvD/PADDB4MRUVHPOVLjGdIagox0dF0bN8OjNFqFiIi0miXjBh2xPwW6/Gwtnc/LnzkL7z8yNMORSbhRgWyNEynTv5VLd54o9qnByT3IC42hqKiYq1mISIiAXP12JFHFMlfpw5m7tCzuOCRPzNz5n8dikzCidfpACSE3X03/PzncPXVEBFR6amy1SwyMrNJ8iVo62kREQmYq8eOpGVcLPMWrig/Nv/kYWw+tjvbl69ja/Zubrt+vHMBSsjTCLI03Flnwfjx1W4/DT/trqfiWEREAq26dovtviSSdnzP+Q/ew6OPTXUoMgkHKpCl4YyBP//Z34dcWup0NCIi4jJXjx1J58SOlY5t79SFPW3aceGzj/OXh59xKDIJdSqQpfEuvxz+859qn8rMymHR8jVaxUJERJrE7TddVWnHPevxMP2i8RR5Izlt1jv8QUWyNEDIFsjGmGRjzLPGmHeMMT93Oh5X++1v4YEHjlj2LTMrh2kz0kn/31c8++p/WLH2W4cCFJGmpHwsTrvr1omV2i1KIyJ49dKf8elpwyk+cJA7H3zKwegkFDlSIBtjXjLGZBtjVlc5PtIYs8EYs8kYU+v+kdbaddbam4BxwKlNGa8cxejR/uL4gw8qHc7IzKagsIjcvHwO5ubzbvpcjSSLBBnlYwkXV48dyehzf/rjVxwZSW6LOC6Z8x7npf9HRbLUi1MjyNOAkRUPGGMigKeB84EUYIIxJsUYM8AY80GVn4TDr7kI+BBIb97wpRJj/FtPp6ZWOpzkS8DaUkpKS4nweMB4tB6ySPCZhvKxhImzh6Yx8fLzK21N/eHwC0jeuJbTPpvNb+9/UgM1UieOFMjW2vnA7iqHBwObrLWbrbWFwFvAxdbaVdbaC6v8ZB++zkxr7fnAVTX9LmPMFGPMEmPMkpwc/aVoMoMGQWEhfPNN+SFfYjxjR51JyxYxREdFUVJSQsaOLCUnkSCifCzhJjWlN7ffeGX549wWcTx7zc8ZsnwRCT9k8uizb+hzSI4qmHqQOwMZFR5vP3ysWsaYM40xTxpjnqOWEQtr7VRrbZq1Ni0+XsuNNalFi+CWWyr1Iqem9GbsqDOxWPLy8/lq+WpeeGOWkpNIcFM+lpDmS4zndzddiTk8lLy/dRseven37EzsRJv9e3n8+bccjlCCXTAVyPVirZ1rrb3NWnujtVZ7SwaDK66ArCyYN6/S4YIC/3bUHo8Hg6GgqEitFiJhRPlYgpEvMZ7bb5xQ/rgoKoro/Hx+9fzj9F63SqtbSK2CqUDeASRVeNzl8LFGM8aMNsZM3bdvXyAuJzWJiIC77oLHHqt0OMmXQHRUJNZaSq0FLNHRkc7EKCJ1oXwsYaFsJLlMQUwML19xLePfexPfhrX86W/POxidBLNgKpC/BnoZY7oZY6KA8cDMQFzYWjvLWjulTZs2gbic1Obqq+Ff/6p0yJcYz/UTRnPyif2JjYkiIsJL+qdfqc1CJHgpH0vYqFokf9+lK69e9jPOnf9f9h88xN+efd3B6CRYObXM25vAV0AfY8x2Y8xka20xcAswB1gHzLDWrnEiPmmEyEj/ZL2//KXSYV9iPMf6EomOiqZNqziKS0rVZiESBJSPxQ3KimSPx9+TvLF7b567+iZi8vMoWLuO9+bMO8oVxG28TvxSa+2EGo6n0wRLBBljRgOje/bsGehLS3XatYPnn4fzz4eTTio/nORLwBvh4VBuPt4ID0m+BAeDFBFQPhb38CXG89spE3j02TcA/457vbZuYsxH7/KPiAgALhkxzMkQJYgEU4tFk9FXes0sOhp+/3v/7noV+BLjmTRuFENOSK6045GIuIfysTjJlxjPsJNTyx+vSh7I3KFncdMrz7D0k8958qUZzgUnQcUVBbI44PrrYft22Lv3iKeWrtzAohXrmDYjXX3IIiLSrC4ZMYzOiR3LH88/eRhLBw6ix7bv2JLxAw/941+1vFrcwhUFsmZNOyA2FpYsgbZtKx3OyMymuKSUuBYx6kMWcSHlYwkGt990Fe1btyp//N8zR/JNv1SSv13L3p1Z3P3Icw5GJ8HAFQWyvtJziDEwZgysW1d+SH3IIu6mfCzB4p5fX1epSAY4fu0KJs14mfwDB/nDX59lxdpvHYpOnOaKAlkcdNJJ8OCD5Q/L+pBHDBvMqOGnkJGZrTYLERFxxD2/vo7knseVP54x+gqKvJFc/e5r5Ofm8a+3P1KR7FIqkKVp3XILzJkDGzeWH/IlxpPkSyD906+YM2+xepFFRMQxU666hImXnw9AaUQEr176M/a1ak10YQEAr7wzW59RLuSKAlk9bw5q3dq/mkV25V7jVeu+41BePpGRXvUii7iI8rEEo9SU3uWrKxVHRvLe+WPxlJZy6uLPsaWlPPrsG/zvyyUORynNyRUFsnreHHbjjXDyyXDgAACZWTksWrGW/IICfty9B6xVL7KISygfS7C6euxIhp2cSqQ3ovzYqV8vYPgXnwAw6+MFTH39PYeik+bmigJZgsBjj/nXRobDo8WG+A7tiIqMpHOneGdjExERwb8E3CN/vIW42BhyW8Tx7DU/5+RlCxn69RcArNu0jfsff8nhKKU5qECW5nHttTB9OmzfXr6SRV5eAYVFxWzbkaU+ZBERCRoP/P5G4mJj2N+6Dc9e83PyYmLLn9u9/wB/ePgZB6OT5qACWZpHfLy/SH7kkfKVLPr2PI64FrG0iI3mwKFcPvpsoYpkEREJCg/8/ka6JXViV/uOLB8wiJNWLCbl2zUA5BUU8rsHntIKF2HMFQWyJoUEid/9Di64APhpu09vhIfsH/eQX1DImm8388Ibs1Qki4Qx5WMJJbddN658a+qsjgmMf+9Num/9DoDikhItAxfGXFEga1JIkDjmGBg2DBYsAPxF8pDUFCK9ERhjAENeQYFWtBAJY8rHEmouGTGMiZefT0ZSV1697GdMmvEyvp07yp//4OMFDkYnTcUVBbIEkZwcGD0afvwRgAHJPYiOjsJai7WWgsIioqMjHQ5SRETkJ6kpvbn9xivZ2L03r1z2M/a0aVf+3K69+3nt3dkORidNQQWyNK+kJLj8cnjiCcA/inzqoAHEREcRFxtLdFQkWTm7nY1RRESkCl9iPL+76Uq29OxDYWQUk994nnZ7/Z9XS1dt4O5Hp6rdIoyoQJbmd+ed8OyzcPAg4B9Fjo2JJq8gn8KiIhYtW6s+ZBERCTq+xHhGDR9KidfLxu69uemVZ2h50L/G/6HDW1NrQ5HwoAJZml+3brBkCbRsCfzUixwTHU3H9u3AGPUhi4hIUDp7aBqDBvRh/snDWDbgRKa8/hyekpLy52d9vEAjyWHAFQWyZk0HoeOOg4cegv37Af8oclxsDEVFxXgjPNpZTyRMKR9LOLh67EgmXn4+C0aNYcboKyiNiMCUlpY/r9UtQp8rCmTNmg5CxsDq1fDPfwKUr4085IRkBg3o43BwItJUlI8lXKSm9Ob+399I5JDB9P5uA9e/8TwRxcXlz7/y9kdqFwxhriiQJUj98Y/w+ONw6FD5oaUrN7BoxTrtrCciIiHhtuvGse/ENIq9Xq76z2vlI8kW+Oe0d50NThpMBbI4JyUFzjoL5s4FICMzm+KSUuJaxFBcUqo+ZBERCQl3/uo6Zt9wK3G5hzh7wf/Kjx/Kz+eRZ17VgE8I8jodgLjcG2+Ax///aUm+BLwRHg7l5qsPWUREQsrvf3MDj5WUkLNnH1GFBRRGRoEx/JC9m8eee5PThxzPJSOGOR2m1JFGkMVZHg+88AJMnao+ZBERCWm//d1NRLRtw4T33mT4F5+UHy+1lnkLV2hDkRCiAlmcd/zxcP/9UFAA+PuQv1y6mmdf/Y9mAYuISEi5eeKlzLrwMk5etpChX1fehnrpqg28N2eeQ5FJfbiiQNayQkHupJOgXz945RUyMrMpKCwiNy+fg7n5TH//ExXJImFE+VjCnS8xnsm/mcLcv/yV4V98Ur7bXpl5C1eoSA4BriiQtaxQCLj3Xti1iyRfAtaWUlJairWW/MIi3k2fqwkOImFC+VjcwJcYz2W3TOaLV95iT9v2tNm3t9Lz8xauYOrr7zkSm9SNKwpkCQFDh8Kdd+Jr3ZIzh55IxOGJe96ICDAerWghIiIh56KLzuO8Psfy2+f+Ro+tmyo9t27TNo0kBzEVyBI8Vq+m6PjjWfT1KqKjIzHGEBUVRUxUpFa0EBGRkHT++DF8c8+fmThjGl0yMyo998Xib/QNaZBSgSzBo18/8lvE0XvpQtq1aU2ruBb079ONSeNG4UuMdzo6ERGRBjntzt/w40MPc26FlS0ASkotT7wwXUVyEFKBLMHDGArvuJMzPvmQQwdziY6KJDG+HfMWrtBEPRERCWndbruZH6c+T8uDB2i7d0/58aLiEp577X19zgUZbRQiQaXD+HHs3bGdkaeeyMGSEj789CustSxZuR6A1JTeDkcoIiLSMGefehKef71Cyr+n84/rbuNgy1YA7D94iH+9/REftV/IkEEpnD00zeFIRSPIElyMoe3tv2Wwt5Ss7N1Ya4mI8GBtKUtXbnA6OhERkUY5c+o/yBh2Fje9+gwxebmVnsvevYdZHy/QChdBQAWyBB9rYfJkhmRuwRgoLi7BWti2faf6tEREJOQN+vdb5J52Gl1/yKj2+XWbtmnXPYepQJbgYwz88Y90eWEq3bt0wuuNoH3b1ni9Xi33JiIioc8Yer0zndGP/JnLivbRqV3rI05ZumqDimQHuaJA1s5NoSfztDPYn7mTlgu/orTUUlRcgjfCo+XeREKc8rHIT3zxHTh1zXJ+/8VHpHTvcsTzS1dtYOrr77Fo+Rp9g9rMXFEga+em0JOx80femnwrewam0rJFLAkd25IY357sXXuO/mIRCVrKxyIVeDzw5puwaxc3zJ/NoP5HTkRft2kb/5k9l3+89A7/+3KJA0G6kysKZAk9Sb4E9h/TifZrV9Pp+y1s3pbJmm+38Mo7s5UgREQkfMTEwHvvQVwcV486k0ED+hxxSkFhMfmFhcz6eIF232smKpAlKPkS45k0bhSnxrdm9P8+BMBai7WWj+ct1ldNIiISPlq1gieegN27uXr39mqL5DLaG6B5qECWoOVLjKfHH++g454fOS5jS/lxT4Qm64mISBgyBh56iKuztjDx8vMZnJrCcZ0TjzhtztyFDgTnLiqQJbhFRhJ9zz1c+mMGMVGRtIiNJS42WpP1REQk/Ph88PHH8OCDpH6zlAkXn8u40cOJ8FQu13bm7FG7YRPTTnoS/CZPpssNN3Brzi4yMrOJjo4sH0H2JcY7HJyIiEgAde8Os2fD5s2A/3Nu1PBTmPXxgkqnffDJl7Rv21o7zDYRFcgS/LxeWLsW33PPkX3jzbybPpfSUv/k37GjzlRyEBGR8NK/v//nqadg4EDOPuMMvt30PRu2/LSxiLWW196dA6DPwSagFgsJDd26UfTGm3z27DQOHMrjUF4eBw7l8faszzRhT0REwlNyMlx2GSxbxkUjTicionLZVlJSyuv/+a8m7TUBFcgSEjL3H+STwadz1qcfVTqeV5DPqnXfORSViIhIExo+HJ57Di64AN/+vYw6+5QjTikuLmH6zE81WBRgKpAlJGRkZrNw8GlElBQTUVxcftwYA8bBwERERJrSmDEwbRr4fJw9NI3R556Kx1P5gy+/oJBV6zVYFEgqkCUkJPkSiGjThpevvMG/DA4QEx1Fm1YtGdC3h8PRiYiINKERI8BauPBCzu6RxDWXjsRTZWWLzd9nOhRceNIkPQkJvsR4rp8wmlVrNzL0umvI+OczHDimM0m+BK1kISIi4a9lS0hLgxEjSJ07l/8ldiTjh5/2BNi0dTuZWTn6TAwQjSBLyPAlxjPirKG0mngNKe9MZ8gJ/SolgsysHBYtX6M+LBERCU9/+hOccQaMGcOxnSrvB1BaalmwZJVDgYWfkC6QjTFxxpglxpgLnY5FmtEvf+nft37rVsBfGM+Zt5AX3pjFnHmLmTYjXUWySDNTPhZpBsb4t6R+/HGGDh54RBvAd9t2OBFVWHKkQDbGvGSMyTbGrK5yfKQxZoMxZpMx5s46XOoOYEbTRClBq107+Mc/oKSEzKwcps1IZ/6ib9h38CCRkV6KS0q1FbVIHSkfi4QYjwdSU/HNn8sv5vwbU1pa/lRWzm7tsBcgTo0gTwNGVjxgjIkAngbOB1KACcaYFGPMAGPMB1V+Eowx5wJrAVVCbnTllewsKOJ/7/+XfQcOYozBWsvBQ7l4Izzailqk7qahfCwSekaPxhcB49Lf8U/gOyz9f19pXeQAcGSSnrV2vjGma5XDg4FN1trNAMaYt4CLrbUPAUd8ZWeMOROIw5+884wx6dba0qrnSXjKzMph882/5Ni8fJaeP5bComKMgRNSejH0pAGapCBSR8rHIiEqJobo9A9JPvkUTln6FV+lDQX8m4dMn/kpCR3a6bOwEYKpB7kzkFHh8fbDx6plrf2jtfZXwBvA8zUlY2PMlMN9cUtyctSXGg4ys3L46LOFfDzoVAZ9s4RWB/YDEB0VRVLnRCUEkcZTPhYJBa1a0eaLzym4YjxRhQXlh/MLCjVhr5GCqUBuEGvtNGvtB7U8P9Vam2atTYuPV+EU6sp6jjdt3c7+Vm1YNnAQZ371GQAejyE6OtLhCEXcS/lYxAHt2jH83NO4YfpLDP36i/LDK9Z8y4q132p1pwYKpnWQdwBJFR53OXys0Ywxo4HRPXv2DMTlxEEZmdkUl5TSpnVLCn7cw2dnjSS6qBBvRASFRcW8mz5PXyuJNJ7ysUgI8SXGs/appzjvisvJi4ll5fEnYfDw9gefUVpaisfj4fILzyI1pbfToYaMYBpB/hroZYzpZoyJAsYDMwNxYWvtLGvtlDZt2gTicuKgJF8C3ggPhYXFtG3dkiEjz6Zbr+70WfMNxcUlHDiUy5df62slkUZSPhYJMSkjz+GHV17lok8/pKMtxmLJzcsnv6CQ3Lx8ps/8VCPJ9eDUMm9vAl8BfYwx240xk621xcAtwBxgHTDDWrvGifgkePkS45k0bhQjhg3m+gmjGZDcg30/7mbC+28Qk5cLwL6DhxyOUiR0KB+LhI++F19I7vLlnDX6PIYmtMVUeE59yfXjSIFsrZ1gre1krY201nax1r54+Hi6tba3tbaHtfbBQP0+Y8xoY8zUffv2BeqS4iBfYjxDTugHwLQZ6WyLjmVN736csWg+xhjatIzT/yWL1JHysUh48XU9liGtYxh5+y10+3FnpeeWrdqgz8c6CqYWiyajr/TCU1k/crs2rVlw3oUMXbmE1tFRrN20TbvpiQQp5WORZtCjBxHPT+Wmt16k056fPgvzCwpZtf47BwMLHa4okCU8lfUjH8rNJy/pOL5+bQae6GjiWsRoNz0REXG3MWOI/NujjFm5uNLhH7J3ORRQaAmmVSxE6qWsHzkjM9u/c15xMR3/8ic+uuRKvC1itZueiIi428SJbD62N+1mzqbIG8nBlq1Y9+1WVqz9VitaHIUrRpDV8xa+yvqRfYnx4PVyTJSXUZtWMmhAH6dDE5FqKB+LNK8BKT1J27iGm159hpi8XAqLi3n7g8/UhngUriiQ1fMW/so2EPnPoNPoOv0NFi1crj5kkSCkfCzSvHyJ8fj++iBbu/fihjeeJ7LQv+ybVrSonSsKZAl/ZRP2co7tyvqefenww3Zy8/P56LOFzJm7kMysHDKzcrSjkIiIuE5qvz78cOfd/JDYCV9WJqAVLY5GPcgSFsom7OXm5/PWxRP8Bw/lsXrDZtZu3MKCpavweiLAGLwRHiaNG6Xd9kRExDWGDjmef1x6Nfn5BRy/ZgWrUgayav13+iysgStGkNXzFv7KJux1S/IRFenl4s/SOXnZQv+TFgoKCikoKiIqysuhvHwtcyPiEOVjEWf4EuMZ1L8PEbaUoUsWcNkHb3PggDbWqokrCmT1vLnHjp05FJeUsKpHX875/BM8JSWUWovX6yUiwkPOrj3kFxSwaNlafbUk4gDlYxHnDD1pAC1at+Kl8ZPx7dxB4mOP6rOwBq4okMUd/OseG9q1acW2rj3Z07YdJ69fSaQ3gtR+vTh10ABioqPp2L4dGKN1kkVExFV8ifEM7NuDwpgYXrjmJkoLCli9egOA5ulUoR5kCRv+dY8tu/cexFrLh2dfQOuD+ykuKWHdt1u5aMRpxMXGUFRUjDfCo3WSRUTEdVrFtcAYQ16Llsw8fwwX79zJnqeeZpptSXFJqebpHOaKEWT1vLmDLzGelJ5d8XgMxsC247qzqu8AOhUXgjEUFBQxadwoRgwbrL/8Ig5RPhZx1oDkHrRp2ZLo6CjatGxJn57HEnf//fRa9AVgKSgs0jesuKRAVs+bO2Rm5bB201ZKSy3WgrWWnt9v5mdTn4DSEjIys1i17juSfAkqjkUconws4ixfYjzXXzmai849jeuvHE3ikMFsnvoCI9+bTqclizhwKJdDeXlOh+k4tVhI2CjrQU7o2I59+w/Ss2sXBl02kpivPqXXyqV8dSgXgEXL13L9laNVJIuIiCv5EuMrfQbuO7YbcyfdjD14AGstc79cRt8ex7n6c9IVI8jiDmVrIRcWFtMqrgXnn3Uyqf36kHXjzxn28YcYwGAoKNLXRyIiImWSfAlkdu3Gxu69OfeLT4jftsX1y6FqBFnCRtlayBmZ2ZXaKFpOuILPNmzBWgsGoiMjNUFPRETkMF9iPGNHncnbH3xGVod4Jk77J6/GtWRA3x6uHUVWgSxhperXRgC+YxI4657f0e2pZ1neuz+J8R0cik5ERCQ4pab0Jit7N/OBzyIMVz37GCvSBuC77CKnQ3OEK1osNGta8Hjo9fLzeOd+xlfLV/PCG7O01qOIA5SPRYLXgOQeREdFMr/v8bw9+nK+2LLTtZ+VriiQNWva3TKzcvho3mI+Pv1czp07BywcystzfX+ViBOUj0WCly8xniGpKcRER7PzpJMxB/bjHTMG9u51OrRm54oCWdwrMyuHF96cxZpvt7A0+XhaHdyPb8f3FBUXa7tpERGRKspGkffuO8BubyRbvDEUjBgBublOh9asVCBLWMrMymHOvIV89NlC8vIL8HgMREbyjym/YeexXYnv0F7bTYuIiFRRNoocFRlJXFwLPho9jv2JneCaa5wOrVlpkp6EnbJR430HDmItgAXAGAMdO3LGoi/I6nws2d17ajULERGRKhIT2lNYVExBYRHGGHY88H/EmxIO78IFnvAfXw3/OxTXycjMLv9LbYwhIiKChI7t6BTfkRP798YX354Rcz8ipXdX5i1cwYq13x5xjcysHBYtX6MWDBERcZ2CgiLiWsTStnVL4lrEUlBiYcAAmDqVQxMnsWjZ6rD/fNQIsoSdJF8C0VGR5OUXYK2ltBSycvYAsCMrh4j2Xbh78yY2/3smGZ2PZcnK9YB/iRvwF8fTZqRTXFKKN8LDpHGjXLsOpIiIuE+SL4GYqEiKS0qJifpp74AfzjmP0oceoTjz90y75Iqw/nx0RYFsjBkNjO7Zs6fToUgz8CXGc/2E0Xz02UI2bd1OaWkphUXFGOP/ZqjE6+XTU8+hx9ZNZB7bldJSy/yF3wD+/2veu+8AxSWlxLWI4VBuPhmZ2WGbAESam/KxSPCraeOt7/fnMv+GXzHp6b+yrVMXVg3sE7afj8b6mzRdIS0tzS5ZssTpMKSZlPUi791/kCP+nFsLxpQ/jI2OpqDI/5WSN8Ljb1s2RiPI0iSMMUuttWlOx+Ek5WOR0FP2uVq8cyd5MbG0j4lh4rWXhfRnZE35WD3IErbKZuLGxkTTplUcUZGRJHVKICY6ilat4kjZsZVJH/2bFrGxxMREYa31F8cYhpyYwohhg1Uci4iIHFb2uVrSvgPtO3bgZ/98hNyXpjkdVpNQgSxhbUByD6IjIykqLiEuNoazTxtEq7gWeCO87O3dl77rV9F+RwYHDuUClrz8AopLisFCdHQkGZnZYT8RQUREpK4GJPcgLjaGouISZk64lqRHH2LXa284HVbAuaIHWVzO/PTPhA7tyvuqDuXl8WnaqQybO5vXx14NQGlpKQcP5fH54pXkFRQQ1yKWmKhIjSSLiIjwU3/yqvXfschaXpv4C8b+6lf8MHQonbp3dTq8gNEIsoQ1/0YghvZtWwOmfMLdkBP6kZWzhy8Gn8aBuJYcXjAZr9f//4wWW95yUVxSqg1FREREDvMlxtO2dSswhr3J/Xjm9/fz/d6DsH2706EFjApkCWtJvgS8ER4O5ebjjfCUL1WzYu23/JD9I3nRscwccQkt8vxbaNpSizEGa0sBOJSbB1htKCIiIlJBxc9XYmPoticbTjoJ1q8Pi70E1GIhYa26pWpWrP2WV/89m9LSw6PGRUXc+fRDPH79rynw+Th98EC+/mY91lqKiksoLi5x+C5ERESCS9XP14TEeHj4YYrPOYd3b/g1u9u0D+mVoDSCLGGvrKWi7C/ouo3bKi37VhwZyaIThjD8y/8Bhrz8QgAiIjx4jKGkVC0WIiIiVVX9fGXiRLZPuJpBcz8mrkVMSLcoqkAW10nudRw/zdzzm3vKWaSuWk6LXT/yQ9aPAFhrsViiIyOJjo4M+a+LREREmlrU7b/lszHjidm6hbiCvPIWxVBru3BFi4V2bpKKEjq0o2WLWPIKCvAYQ0qvrlhg+riJHDQefvwhG2Pg+OSeHBPfgcSE9qR/+pW2nhYJAOVjkfDmS4xn0hUXYP94Nx0/fZ/oWyeRmZXDtBnpIfU56ooRZGvtLGvtlDZt2jgdigSBjMxsvF4vx8R3IK5FC/r27Epyz65s6NOfjnt20fLgAayFqKgoRpx5MgUFReVbT4fy10UiwUD5WCT8+RLj6Tz1GaL79oVLL2X7th0h9znqigJZpKLqVrZI8iUQHR3J4OWLOOvLzzDGHG7FqHklDBEREamBxwMvvggJCXTN3R9yn6OuaLEQqai6lS0ys3IY0LcH6y69gp/95U6yr5tMVvZuMjvkHHF+9q49zFu4gsT4dsTFxpb/Ra94PREREdfzeuFf/yKhtJQrly3jq+4pJPfuGhKfkyqQxZV8ifHlf0Ezs3J44c1Z7DtwEGthefJA2rz+Gv89aySLlq/l+itHl5/vXyJujn8Cn7W0iI0lOsoLFjAmZHqrREREmssP2zKIfPIfHNujD+mXXEFCh3ZB/zmpFgtxvYzMbAoKizDGv7LFB8MvZO7QszEYCoqKKvVKlS0R5/GUrYJhKSgsoqCoKKR6q0RERJrL93sP8sYNv6LP+lUMmTMrJD4nVSCL6yX5EoiOiixfGzkvLo743dmcuGIx0ZGRlXqlknsdhzGmfJMRMERHRRIdGRlSvVUiIiLNJcmXQFHr1rwy5de0yM8jqVNwjx6DWixE8CXGc/2E0axa/x1YSExoz4FlKxhyy43s+uuf6VTha6DUlN6AfyRZPcgiIiJHV2kuz03X4MveCV/Mh8suczq0GqlAFlfJzMopL2QBVq3/jgMHcwFoFdeCxIT2ZOXs5kC7juxMPZH8Rx8j65e/ZOOW7WTn7KZlXCzHxHdg2MmpAEd8TZS9aw8ZmdlER0dSUFCkgllERATKPwszMrPx7ttDws03Q2wsXHCBw5FVTwWyuEbFhcrBUlxSwsFDeZW2nTYGyh5+1+8kLv/gbf7VqZf/iXKbWLB0FV5PxOHjFiyUlJZyMDeP2Oho8goKiGsRS0xUpCbtiYiI61XdLOSGaf8ifuLP4J134IwznA7vCOpBFtfIyMwuX6i8oLCIgoLCKhtO+4vjsmM7Ezrxj2tvrVIc+xUUFJZPzCubpBcR4SnfntpaizfCo0l7IiIiVP4MLi4pZfMxSfDWWxAT43Ro1VKBLK5RccOP6KhIoqOjsFXOMYZKx6IKC/n5v/6Jt6io0nnR0VHlE/PKJumVlJRijMFgMMaU/1+yJu2JiIjbVbvp1tlnw+DB8Mc/wvr1TodYiVosxDWqbvgBtfQgVzz25Wwm7Pqer4acUd6DPCC5B0Cla1XsPVYPsoiIyE+q26SrXO/ecN558PnncNxxzgVZgQpkcZWyDT/KJusN6Nuj/PGq9d+xact2DubmUVJSyqCBfUjo0I6MyVMYeNftZFxwEUvXb+b7zGw2Z2QSFRlJt2M7lf9lH3JCv/LrhsIi6CIiIs2p4iZdZTKzcsgYmEbfKTfS5txz4auvoEMHhyL8iQpkcZ2qEwVGDT+Fmf/9gr37D1aasLfm2820bNECr9fLWd368s2nn3OgbXsgj9179wOwesNmWsTGEhcbzajhp5D+6Vfl19XkPBERkZpV/Dz+tE0XbrznXjq0bet0WEAI9yAbY840xnxujHnWGHOm0/FI6Kg6UWDdxm3+nfSqnGctFBQWEtcihvfOvYgDca3wlJRUc0Vbfp2K19XkPHEL5WMRaYiqn8eb+p8AGRkwZgzk5joamyMFsjHmJWNMtjFmdZXjI40xG4wxm4wxdx7lMhY4CMQA25sqVgk/VScKJPc6zr+TXpXzjIHoqCj/RLzoSK76z+ucuGppNVc05dc5YgKCSJBTPhYRp1Q7ce/YY6F1a7j0UigsdCw2U/Er5Wb7pcacgT+ZvmKt7X/4WATwLXAu/gT7NTABiAAeqnKJ64AfrbWlxphE4P9Za6862u9NS0uzS5YsCdyNSMiquGFIxR7kAwdzj+xBPnze5mmv0fexv/LsXQ/SoWO78h7kst30KvY2a3Ke1MYYs9Ram+Z0HKB8LCLOqvZzs7gYLr8cOnaE559v0t9fUz52pAfZWjvfGNO1yuHBwCZr7WYAY8xbwMXW2oeAC2u53B4guqYnjTFTgCmHHx40xmyo8HQbYN9R/r0j8GOtN1Q3Fa/ZmHNreq6641WPBeP9Hu28utxXdcfqcu8Nu9ffXF+XGGui9/boz9fnvaz6OND321TvbXBM00b5uJHnhsKf4drirO95wZeP6xZjQ85VPq7+ePP/WX7hhaPFWJ3G52NrrSM/QFdgdYXHlwEvVHh8DfBULa8fCzwHTAfObGAMU+vw70sCdL9TA3FuTc9Vd7zqsWC836OdV5f7Otq91nS/em+D772t73vZ1Pfb1O9tsPwoHwf2fQ2mP8P1uV/l4/q/t/V8P0Pqva3ve9nU99vc+ThkV7Gw1r4LvNvIy8yqw78HSn2uWdu5NT1X3fGqx4Lxfo92Xl3uq7pjdb33QNB727Dz6nq/4fzehgXl46Med/LPcH2uqXx89OeUj+v+ONTe20oc6UEGOPyV3gf2p563U4D7rLUjDj++C8D6v9JzjDFmiQ2SXsHm4Kb7ddO9gu5XaqZ8HJzcdL9uulfQ/YaCYFrm7WuglzGmmzEmChgPzHQ4JoCpTgfQzNx0v266V9D9St0pHwcHN92vm+4VdL9Bz6lVLN4EzsTftJ0F/Mla+6IxZhTwBP6Z0i9Zax9s9uBERFxE+VhE5EiOtViIiIiIiASjYGqxEBERERFxnApkEREREZEKVCCLiIiIiFSgArmejDFxxph/GWOeN8YcdTvVUGaM6W6MedEY847TsTQHY8wlh9/X6caY85yOp6kZY5KNMc8aY94xxvzc6Xia2uG/u0uMMbXtBCchRPk4fCkfh7dQyMcqkAFjzEvGmGxjzOoqx0caYzYYYzYZY+48fHgs8I619gbgomYPtpHqc6/W2s3W2snORBoY9bzf9w6/rzcBVzgRb2PV837XWWtvAsYBpzoRb2PU8+8twB3AjOaNUupL+Vj5GJSPnYi3McIxH6tA9psGjKx4wBgTATwNnA+kABOMMSlAFyDj8GklzRhjoEyj7vcaDqZR//u9+/DzoWga9bhfY8xFwIdAevOGGRDTqOO9GmPOBdYC2c0dpNTbNJSPlY9/onwcGqYRZvlYBTJgrZ0P7K5yeDCw6fD/tRcCbwEXA9vxJ2UIwf9+9bzXkFef+zV+fwU+stYua+5YA6G+76+1dqa19nwg5L6erue9ngmcDFwJ3GCMCbm/u26hfKx8jPKx8nEQ8DodQBDrzE8jE+BPxEOAJ4GnjDEX0DT7qjuh2ns1xnQAHgROMMbc5fQ2swFU03t7K3AO0MYY09Na+6wTwTWBmt7fM/F/RR1NaI5YVKfae7XW3gJgjJkE/GitLXUgNmk45WPlY+Xj0BPS+VgFcj1Zaw8B1zodR3Ow1u7C3//lCtbaJ/F/4LqCtXYuMNfhMJqVtXaa0zFI4Cgfhy/l4/AX7Pk4KIe1g8QOIKnC4y6Hj4UjN90r6H7D+X7ddK9u4qb31U33CrrfcL7fkL5XFcg1+xroZYzpZoyJAsYDMx2Oqam46V5B9xvO9+ume3UTN72vbrpX0P2G8/2G9L2qQAaMMW8CXwF9jDHbjTGTrbXFwC3AHGAdMMNau8bJOAPBTfcKut9wvl833aubuOl9ddO9gu43nO83HO/VWGudjkFEREREJGhoBFlEREREpAIVyCIiIiIiFahAFhERERGpQAWyiIiIiEgFKpBFRERERCpQgSwiIiIiUoEKZBERERGRClQgi9SBMWa8MabAGBPpdCwiIm6mfCzNQQWySN0cD6y11hY5HYiIiMspH0uTU4EsUjfHA8udDkJERJSPpempQBapm1RgRdkDY8zZxphdxpgnjDERjkUlIuI+qSgfSxNTgSxyFMaYeKAThxOyMeZW4APgTmvtr6y1JQ6GJyLiGsrH0ly8TgcgEgKOByyw1hjzPDAWGGWtnetoVCIi7qN8LM1CBbLI0aUCu4F3gXhgiLV2k6MRiYi4UyrKx9IMVCCLHN3xgAFOB85RMhYRcYzysTQLY611OgaRoGaMWQn8G+gJDAMGWWtznI1KRMR9lI+luWiSnkgtjDHRQF9gJTAF+BF42xijb19ERJqR8rE0JxXIIrVLASKBVdbaPGAM0A943NGoRETcR/lYmo0KZJHaHQ8cAr4DsNZuA64Afm6MmeRgXCIibqN8LM1GPcgiIiIiIhVoBFlEREREpAIVyCIiIiIiFahAFhERERGpQAWyiIiIiEgFKpBFRERERCpQgSwiIiIiUoEKZBERERGRClQgi4iIiIhU8P8BZGZSOW6OMxsAAAAASUVORK5CYII=\n", 94 "text/plain": [ 95 "<Figure size 720x360 with 2 Axes>" 96 ] 97 }, 98 "metadata": { 99 "needs_background": "light" 100 }, 101 "output_type": "display_data" 102 }, 103 { 104 "data": { 105 "image/png": "iVBORw0KGgoAAAANSUhEUgAAAsgAAAFgCAYAAACmDI9oAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAB/YElEQVR4nO3deXxU1dnA8d+ZmSSTPexLgEAEWQyKgriVuFQrqIjgUrXWiopdtMXWVm3fulVF61KhxbauULXWWgviCtpSDaCooAgIsogJi2wJMNkzM3fO+8fNJDPJZJ+ZO5N5vnzmQ3LnLueScPLk3Oc8R2mtEUIIIYQQQphsVjdACCGEEEKIWCIBshBCCCGEEAEkQBZCCCGEECKABMhCCCGEEEIEkABZCCGEEEKIABIgCyGEEEIIEUACZCGEEEIIIQJIgCyEEEIIIUQAh9UN6CylVDrwZ8ANvKe1/rvFTRJCiIQk/bEQoruJqRFkpdSzSqkDSqmNTbZPVkptUUptV0rdXr95BvCK1noWcGHUGyuEEN2Y9MdCiEQWUwEysBCYHLhBKWUHHgemAGOAK5RSY4BBwK763YwotlEIIRLBQqQ/FkIkqJhKsdBaFymlhjbZPBHYrrXeAaCUegmYBuzG7JTX0Uqgr5S6AbgBID09ffyoUaPC3/AQXC7YuxeidDkhRBxZu3Ztqda6j9XtaE136o9bs3s3aA2DB1vdEiGEFVrqj2MqQG5BLo0jE2B2xCcBfwTmK6XOB15v6WCt9ZPAkwATJkzQa9asiWBTG9XVQf/+8MYb5t9CCOGnlCqxug2dFJf9cWt274Zjj4V33oGePa1ujRAi2lrqj+MhQA5Ja10FzLS6HS1JSYE77jBHkiVAFkJ0Z7HeH7dm0CC4+mpYtw7OOsvq1gghYkU8BMh7gMCHX4Pqt7WbUmoqMHX48OHhbFebfvEL8PmiekkhhIikuO2PWzN3rvm31qCUpU0RQsSIWJukF8onwAil1DClVDJwOfBaR06gtX5da31DdnZ2RBrYEq8X8vPNUWQhhOgG4rY/bsu998JTT1ndCiFErIipEWSl1D+AM4DeSqndwF1a62eUUjcBywA78KzW+gsLm9luDgeMHQtvvglXXml1a4ToGI/Hw+7du6mtrbW6KXHL6XQyaNAgkpKSrG5Kh3W3/rgtZ5wB114L110HdrvVrRHdgfShsaWj/XFMBcha6yta2P4W8FZnz2vlI73p02HRIgmQRfzZvXs3mZmZDB06FCXPnTtMa01ZWRm7d+9m2LBhVjenw7pjf9yab30LeveGV1+Fiy+2ujWiO5A+NHZ0pj+OhxSLLrPykd6FF8KAAVG/rBBdVltbS69evaRj7ySlFL169ZLRoyZiNcVCKbj7bjM1TohwkD40dnSmP46pEeTuqHdv+NOfZPKHiE/SsXeN/PvFl3PPNfvq8nLIyrK6NaI7kD4gdnT0a5EQI8hKqalKqSddFs2WW7IEfvxjSy4thBAxxer+uC2vvAJXhEwuEUIkkoQIkK1+pHf88WanK4/uhOiYU0891eomiDCzuj9uy9Sp8OmnsGGD1S0RouukD+08SbGIgiFDYNgwKCqSQvSieypcUEiJq/liRHnZeRTNLOr0eT/44IOuNEuIDnM64Wc/g4cfhuees7o1IlFIHxp7JECOkh/+0MxrE6I7KnGV0Cet2VL2ITv8jsjIyKCyspL33nuPu+++m969e7Nx40bGjx/PCy+8IPl9IiJ+9CNYutTqVohEIn1o7EmIFItYyHm7/nq46CLLLi9E3Pvss8+YO3cumzZtYseOHaxatcrqJolOiIX+uC09esCll8Lnn1vdEiHCR/rQjkmIADlWct6uugo++cTSJggRtyZOnMigQYOw2WyMGzeO4uJiq5skOiFW+uO27NsHZ54Jhw9b3RIhwkP60I5JiAA5VgwZYi4aIoTouJSUlIaP7XY7Xpn1KiJo0CBzwt5f/2p1S4QID+lDO0YC5CiaMcMMkLW2uiVCCCHa8stfwh//CLLWixCJRybpRdH48TByJLhckJNjdWuECJ+87LwWZ2ALEa/GjoXXX4eAgTchIkL60NijdAIMZyqlpgJThw8fPmvbtm1WN0eIuLB582ZGjx5tdTPiXqh/R6XUWq31BIuaZKl4648NA/72N/jBD8But7o1Ip5IHxp7OtIfJ0SKRSxNCikuNmdHCyFEIoql/rg9bDZ44gl47TWrWyKEiKaECJBjyaBB8N57UNK10oZCCCGiQCm49Vb4/e9l/ogQiUQC5ChzOMyZ0YsXW90SIYSIH5UVXp5/opjKiujPvL/oIigrg/Xro35pIYRFJEC2wGWXmTU2hRBCtM+q5aVs+NTFquWlUb+23Q5r1sBxx0X90kIIi0iAbIHJk+HBB61uhRBCxIfKCi+rV5SRlZ3E6hVllowiZ2fDY4/Bpk1Rv7QQwgIJESDH4tKmTzwBL71kdSuEECK6OtMfr1peis/QJCXb8BnaklFkgOpqePhhSy4thIiyhAiQY3HWdHY2PP+81a0QQojo6mh/7B89TnGaNdZSnHbLRpF//GNYsgR27476pYUQUZYQAXIsOu88WLECysutbokQQsQu/+ix3a4AsNuVZaPIPXvCddfB++9H/dJCiCiTANkiWVlmNQvJZxMiuq699lr69u1LQUFBi/vMmzePgoICjjnmGObOndvsfcMwOP7447ngggsi2FIBsLO4Gp/PHEn2v3w+c7sVHnoIvvc9Sy4tRExoTx/aXkuXLmXkyJEMHz6cB1uYnNVafxzJvliWmrbQ3/9udQuESDzXXHMNN910E1dffXXI9zdu3MhTTz3Fxx9/THJyMpMnT+aCCy5g+PDhDfvMmzeP0aNHUy6PgCJu1uz8sJ6vssLL4hd3M/3KQWRkdvxHoFKwcCEcOQI33xzWpgkRF9rqQ9vLMAxuvPFG3n33XQYNGsSJJ57IhRdeyJgxYxr2aas/jmRfLCPIFjIMuPJKqK21uiVCxKZt27YxdOhQtm/fDoDH42HcuHHs2rWr0+csLCykZ8+eLb6/efNmTjrpJNLS0nA4HJx++uksWrSo4f3du3fz5ptvcv3113e6DcI64SgXN368OZJcVxfGhgkRAVb0oe318ccfM3z4cPLz80lOTubyyy9nyZIlQfu01h9Hui+WANlCdrs52eM//7G6JULEphEjRnDDDTewbNkyAObPn8+FF17I4MGDm+07adIkxo0b1+z1nw7+BysoKGDFihWUlZVRXV3NW2+9FfTD5Oabb+ahhx7CZpPuM96Eq1zc2LFmTWSZaC1iXUf60HBpb1+8Z8+eoHYMGjSIPXv2BO3TWn8c6b44IXr4WCzz5jd9uqyqJ+LD3dyNCviztv5P4La7uRuAgQxs2Dae8QDcwA1B+37DN+26bkFBAVu2bOHQoUM888wz3Hbbbdxxxx3N9luxYgXr1q1r9jr77LM7dJ+jR4/mtttu4zvf+Q6TJ09m3Lhx2O1mBYU33niDvn37Mn78+A6dUzSysj8OZ7m4X/8aYvBHiohhsd6Hnn322RQUFDR7NR3VbUu4+mJouT+ORl+cEDnIWuvXgdcnTJgwy+q2NDV9uuQii/hwd/2fpjS62bZQHfeT9X866uijj+bxxx/n7rvv5pe//CUVFRV4PJ5m+02aNImKiopm2x955JEOd8zXXXcd1113HQC/+c1vGDRoEACrVq3itdde46233qK2tpby8nKuuuoqXnjhhQ7fV6Kyqj9uqVzcaWf17lQucmGh+XK7ITk53K0V3VGs96EdfdrWkvb2xbm5uUFP53bv3k1ubm6z40L1x9HoixMiQI5lQ4fCJ59Y3QohYtdRRx3Fp59+isvlYu7cubzzzjuMGzeu2X4rVqwI2zUPHDhA37592blzJ4sWLWL16tUAPPDAAzzwwAMAvPfeezzyyCMSHMeJhnJxzsZycXW15ijyudP6d+6cq+D2282SnULEqvb2oeHS3r74xBNPZNu2bXz99dfk5uby0ksv8eKLLzbbL1R/nJOTE/G+OCFSLGLdxx/D/PlWt0KI2JSUlERWVhYPPvggNpuNdevWdalzv+KKKzjllFPYsmULgwYN4plnngHgvPPO45tvzFGbiy++mDFjxjB16lQef/xxcnJywnAnwkqRKBd38smwdy+sXBnGhgoRZtHqQzvK4XAwf/58zj33XEaPHs1ll13GMcccA8RGf6y0bj60311NmDBBr1mzxupmNLNxI5x/PhQXmyWEhIgFmzdvZvTo0VY3A4AhQ4ZQUlKCUorrrruOp556Km4myYX6d1RKrdVaT7CoSTEhVvvjjvrrX+Gtt+C116xuiYg10ofGno70x4n3rxODjjkGUlLg00+tbokQsae4uJi8vDxU/W+PzzzzTEJ27CL6Kiu8PP9EcavVLn7wAzjjDEigsSYRZ6QP7Rz5F4oBSpmT9eQxnRDNDR06NKz5xUK0V3tqJqemwi9+ATt2RLFhQnSA9KGdIwFyjJgzB2bPtroVQgghoGM1k10uOPFEaFLCVQgRxyRAjhF2Ozz2GNQvdiOEEMJCHamZnJ0NV18N8+ZFsYFCiIiSADmGbN8O//qX1a0QQojE1lLN5NZGkX/+c3jmGVk8RIjuIiEC5FheSS/QjBmyqp4QonuLh/64oWayvbFmclujyHl58Mor4HRGq5VCiEhKiABZa/261vqG7Oxsq5vSqsJC2LkTSju/+qkQQsS0eOiPO1sz+cwzYdkyqKuLUkOFEBEjK+nFkKQksxayjEAIIYR1Zs3Ob3Ofygovi1/czfQrBwUtVf3nP5uDHNdeG8kWCiEiLSFGkONJRYU5WU8I0XELFy5sWH1JiEhqqQTcrbfCww+Dz2dRw4ToAulDG0mAHGMyMuDuuyXNQsSn9iysEEnSuXc/hQsKyZub1+xVuKDQsja1VgLuzDMhPR3efdey5ok4Jn1o7JAAuQ0+ojsMkJoK55wDr78e1csKERbtWViho/7whz9QUFBAQUEBc+fOpbi4mIKCgob3H3nkEe6++25eeeUV1qxZw/e+9z3GjRtHTU1N2NogrFPiKqFPWh/6pPWhZ1bPho9LXCWWtam1EnBKwRtvmP24EB0lfWjskAC5FcUUM57xePBE9brTp8Obb0b1kkJ0WUcWVmivtWvXsmDBAj766CNWr17NU089xeHDh0Pue8kllzBhwgT+/ve/s27dOlJTU7t8fRE7vClevrjsCyoGVFjajvaUgOvfH/79b/jkE6taKeKR9KGxRQLkVgxlKL3oxQIWRPW6l1wCL7wQ1UsK0WUdWVihvVauXMn06dNJT08nIyODGTNmyJKpCcpR52DY/4ax4+wdHB4a+gd8NLS3BNzBg3D//Va0UMQr6UNjiwTIbZjDHH7H76gheo8aUlJg9Wr48MOoXVKILunMwgqddeTIEXwBM6Bqa2vDfg0RmzL3ZjLirREon7KsDW2VgPPnkF5ysZcPP4Qvv7SsqSKOSB8aeyRAbsNEJnIHd1BLdL+BtmyBP/0pqpcUotM6s7BCe0yaNIlXX32V6upqqqqqWLx4MVOmTOHAgQOUlZVRV1fHG2+80bB/ZmYmFRXWPoIXkZVWlkbOzhwqCyqZwxw0OqrXnzU7n3vnFTR7+UvD+XNIP1tdyo03wpIlUW2eiFPSh8YeqYPcDj/kh2xjGzZsZBOd4vbTpsHtt4PbDcnJUbmkEJ0WOKrWdHtXnHDCCVxzzTVMnDgRgOuvv54TTzyRO++8k4kTJ5Kbm8uoUaMa9r/mmmv40Y9+RGpqKh9++GHC59B1B3nZeSEn5I0oHcFLvMRBDvIoj2KLgfGepjmkN/+2N9k58mNWtE360NijtI7ub99WmjBhgl6zZk2njr2O6xjIQO7l3jC3qmWnnQZ33QXf+U7ULilEg82bNzN69GirmxH3Qv07KqXWaq0nWNSkmNCV/tjvMIeZylRmMINf8Iswtazzli3ZxwfvlZKW7qC6ysupZ/QGZ3/WrYPbbrO6dSLapA+NPR3pj63/lTtO3MVd/Jk/s5/9Ubvm4sVw9tlRu5wQQsSVHvTgHd7hh/yQPeyhiirL2tJSDumgXC8PPQTl5ZY1TQjRCRIgt9MQhvB9vs/DPBy1a/buDU8+CYYRtUsKIURcSSONdNL5K3/lHM7hEIcsaUdLOaS7vyrl3HPNvlwIET/iNkBWSuUrpZ5RSr0SqWs0XdHmLu7iN/wmUpdrxmaDv/xFqlkIIWJbNPrjttzDPZzKqUxiErvZHfXrt1TdYse2KkYP+4bt22TtaSHiiSWzB5RSzwIXAAe01gUB2ycD8wA78LTW+sGWzqG13gFcF8kO2T8buW9/J+dO608PerCTnTzHc9zMzZG6bJDp081Ui299KyqXE0IkmHjpj9tiw8YjPEI/+rGZzQxiUFSv769i0dSyJfv4+q39TD/Phs/XH1vcDksJkVis+q+6EJgcuEEpZQceB6YAY4ArlFJjlFJjlVJvNHn1jXQDW1rRJpts5jCHTWyKdBMAmDEDFi2CBJpLKYSIroXEeH/cEb/iV5zDOTzGY6yha5MAuyrw58gbSyo4cYLGJwPJQsQFSwJkrXURNEsUmwhs11rv0Fq7gZeAaVrrDVrrC5q8DkS6jS2taJNNNrdyK7/lt5FuAgBjx8K774Kyri6+EKIbi4f+uDPyyWcKU/gP/7GsDYE/R3pk1eFyeXnrLcuaI4TogFh62JML7Ar4fHf9tpCUUr2UUn8FjldK/bqV/W5QSq1RSq05ePBguxpy5hNTeHHJx3xV8SXr96/nq4oveXHJx5z5hDnIciM3spvdUZkMopS5st7//hfxSwkhhF/M9MedNY1p/Jt/cyVXUkRRRK8VStOqFs5UO+NG7+OBB+RxoBDxIJYC5A7RWpdprX+ktT5Ka/1AK/s9qbWeoLWe0KdPn3adW301CodKxm634bA5zL9VMuors3ZeKql8xEf0oEd4bqYNu3bBz34WlUsJIUSHRbI/7opCCnmf95nIRCqI7upgoapajD7KxcB+lSx4vDgiSwgLIcInlgLkPcDggM8H1W/rMqXUVKXUky6Xq137Z1cfhdJ2kjyZDS+l7WRXH9V4ThRXc3VUHt+dfDKUlsK2bRG/lBAJoba2lokTJ3LcccdxzDHHcNdddzXbZ9euXZx55pmMGTOGY445hnnz5rXrvW4iZvrjrhrNaJw4uZiL+R2/i9rS1KGqWgAcN2IXaz6q6vISwkJY7dprr6Vv374UFBS0vXMr2tOftrRPRPtirbUlL2AosDHgcwewAxgGJAOfA8eE85rjx4/X7THksSF6/BPjm72GPDYkaL+X9Ev6RH2i9mlfq+erKPfo5/76ta4o97Tr+qH88Ida//73nT5ciA7btGmT1U2IGJ/PpysqKrTWWrvdbj1x4kT94YcfBu3zzTff6LVr12qttS4vL9cjRozQX3zxRZvvNRXq3xFYoy3qe0O9Yrk/Dpe9eq8+Th+nb9I3aUMbUb22X0W5R99x80adke7WN83c2qWfCSL2dec+VGut33//fb127Vp9zDHHdOk87elPW9qnI32x1h3rjy0ZQVZK/QP4EBiplNqtlLpOa+0FbgKWAZuBl7XWX1jRvva6lEvx4GExi1vdz18urisjBr/9Lcyc2enDhYhL27ZtY+jQoWzfvh0Aj8fDuHHj2LVrVxtHtk4pRUZGRsM5PR4PqslM2AEDBnDCCScAkJmZyejRo9mzZ0+b78Wb7tIft6U//Xmf99nKVjay0ZI2rFpeikJz0gkuPvikl4wii4iLVB8KUFhYSM+ePbt8nvb0py3tE8m+2JI6yFrrK1rY/hYQ9jm+SqmpwNThw4eH9bw2bMxlbqvLmzYtF3faWb3JyOz4P/ugQbBqlflxFFL3hIgJI0aM4IYbbmDZsmUMHz6c+fPnc+GFFzJ48OAWj5k0aRIVFc3zTR955BHODli73TAMxo8fz/bt27nxxhs56aSTWjxncXExn332Wch9WnsvHnSX/rg9sslmKUtRKBaykEu4hAwyonLtwEl7p044zCN/yefdpdsZNzGHd17bx/QrB3XqZ4MQrelMHxpO7e2P/drTn7a0T7j74ljKQY4YrfXrWusbsrOz27V/XnYeB6sPNnvlZec12/d0TudczuUrvgp5rpbKxXXGX/4C//pXpw8XokvuvtusquJ/rV1rvgK33X23ue/AgY3bxo83t91wQ/C+33zTvusWFBSwZcsWDh06xDPPPMNtt93W6v4rVqxg3bp1zV5NO2O73c66devYvXs3H3/8MRs3hh5VrKys5OKLL2bu3LlkZWW1+z0RWkf743BTKHz4WMlKvs23KSU6o7iBk/bS0wwun/YNdgxeXriry08YRXyIlz707LPPpqCgoNlryZIlHb7n9vbH0L7+tKV9ItEXy6+rIRTN7FhJoLWs5VIuZStbSSGlYXvTMj8pTnuXRpFnzDCD5J/8pMOHCtFld9/d2HkHCrWITaiO+8knzVdHHX300Tz++OPcfffd/PKXv0QpxY033khKSgppaWncd999Qft3dMQiJyeHM888k6VLlzabbOLxeLj44ov53ve+x4wZM9r9nog9hQsKKXGVNHyu0bhOczFkxBB25OygP/0jev3ASXsAuf1c7N7rpKKyioEDuvaEUcSHeOlD//Of8BUfaG9/3J7+tKV9ItUXy//EMJjIRMYylr/yV2Yzu2F7w4iBs7HMT12tOYp87rSOd8bnnmvmIR86BGFI+xEiLhx11FF8+umnuFwu5s6dy6OPPso111zDiSeeyCWXXNJs/xUrVrR5zoMHD5KUlEROTg41NTW8++67zUZVtNZcd911jB49ml/84hftfk/EphJXCX3SgvPT+n7Wl5IDJfQ8vyc+fNgi+FA11FLUZ0yqIVkdIS/PhafK2+mfDUK0pqN9aDi1pz9uT3/a0j6R7IsTIsUiGmWF7ud+HuCBoHzkUGV+fD5ze2ekp8OyZebfQiSKpKQksrKyePDBB7HZbHzxxReMHz8et9tNWlpap865d+9ezjzzTI499lhOPPFEzjnnHC644IKG98877zxWrFjB888/z/Llyxk3bhzjxo3jrfpl0FatWtXie6J10S7z1pa0bWloNCdxEh/xUdSuW1nhZeTQb/h4XW+8RuMTRqmPLMItEn0owBVXXMEpp5zCli1bGDRoEM8880ynztNaf3reeefxzTfftLhPJPtipUON7XdTEyZM0GvWrInY+T/jM8YxDkXk1oWurTVzlk47LWKXEAKAzZs3M3r0aKubAcCQIUMoKSlBKcWbb77Ja6+9BsBPfvITjjvuOItb17pQ/45KqbVa6wkWNSkmRLo/bipvbl6zEWSAg9UHKbm5hDd4g5nM5O/8ne/wnYi3Z9mSfXzwXikvvjqME451ccLYciorvKRn2Lnp9hGSatENSB8aezrSH8v/wDA6nuP5C3/hUi6lN70jco3qapgyBfbulZFkkRiKi4vJy8trKMN2/vnnc/7551vcKtHdXMAFLGYx3+f7fMZn5JAT0ev5nzBOOXM3aaleKis0NdUGrsNuSbUQYSV9aOdIgBxm61nPV3zFIzwSkfP37AknnWSmWsi8IJEIhg4d2q48NiG66lt8i81sJoUUiiiikMKIXSswJ3nFCqir87Jq6RbsNpmwJ8JL+tDOkRzkMLuDO3iWZ9nN7k4dX1nh5fknilvNQ5s+HRa3vjaJEELEJKtykNtbvtOJExcurud67uKuqCxNvWsX/OJmX9hKggohui4hfj3VWr8OvD5hwoRZkb7WQAZyAzfwKq9yEzd1+Hj/qnt9+ztbfMR22WVQv3CMEELElWj2x4E6Ur4zhxxWspIpTOEAB5jPfOzYI9a2KZO9/OTHcOBQOkPT60hx2vnw/VL27Krmsh8MkZFkISyQECPI0TaHOZ0KjpuuutfSKHLv3jBsmDnqIEQkJdIk3kiQf7/41Ze+/I//0Z/+ER9F/nhFKScfX8q6L3IAsyRoVaXB+rWygEi8kz4gdnT0ayEBcgTYsLGSldzADR06riOr7v3tbzBnTldbKkTLnE4nZWVl0sF3ktaasrIynE6n1U0RnZRFFndxF/vZz/f4HhU0X/AgHHYWV3NCwWHOOmUPlRVeyl0eamsN7DYlpd/imPShsaMz/XFCPLdRSk0Fpg4fPjxq1zye47mUS/mMzzie49vcv6Or7k2fDoWF8PjjYJNfc0QEDBo0iN27d3Pw4EGrmxK3nE4ngwYNsroZMcWK/rir+tOfDDI4i7N4i7foQ/NycV3hn7D36aewdCmMP8YsAZeW7qBaFhCJW9KHxpaO9sdSBzmC5jOft+r/tMVfEzMtvTEYrq7ycuoZvVvsGAsKzGUnTz01bE0WQkSB1EGOfn/cVRrNndzJUpbyER9FZNW9Awdg5EjN9ZdvoUeOxm5XuN0+jhxyc+u9o+g3QJ5GCBFuUgfZAjdwA1/xFW7cJJPc6r6Bq+413d6Sxx+HIUPC0lQhhBD1ChcUUuIqabY9LzuPxTMXo1DsZjeDCO/Tgb594Vsn17D6055ccM4hAGqqDepqfby8cBc//fWIsF5PCNEyCZAjKJlkHuMxvuEbBjCg1RX2Amtittfpp8PWraA1qMgt3ieEEAmlxFUSctW9ElcJgxjE53zOuZzLYhZzCqeE9donTyjlby/kUFnhxefTVJR7sNlgx/ZKKiu8UtFCiCiR7NUI02gu5ELe5M0un6tpjWSt4TvfgY0bu3xqIYQQ7XQcx7GQhUxjGm/zdljP/X/3DmHr11ncO6+AwrP70LNXMrlD0sjKSpKKFkJEkQTIEaZQ3MEd/B//hw9fl87lr5Hs7ySVMifrLVoUjpYKIYRor8lM5jVeC8vgR1OlpXDiiZpV7zVO3HYk2Xjn9X3s31sb9usJIZpLiADZqpWb/C7kQtJI4yVe6vQ5WqqRPGOGrKonhIgfVvfH4XQyJzOf+XzJl/yZP4ftvL17g6fOzYbNWdjtZv5cYC6yECLyEiJA1lq/rrW+ITs725LrKxTP8ixncVanz9FSjeRTT4UbbzTTLYQQItZZ3R9HQhpp/JE/8ht+E7ZFRU4/9RAfrOlDRblZF7lpLrIQIrISIkCOBaMZTSWVvMM7HT62pRrJlRVe7Ha49lrYvz/cLRZCiMSUl53HweqDzV552Xkh9x/CEFaykv/wH27hlrC0Ye5fBnDFVU5++TvJRRbCCjIdNoqOcISZzGQb20gjrd3H+UeP7U7zUZvdrqir1Q3F4//zH7j3Xli5MlItF0KIxFE0s6jDx/SmN8tZzja24cOHBw8ppHS6DUrB3Lmwb2/wAIk/F3ncxBypiyxEBMkIchRNYAKncArzmd+h4wJrJPtfPl9jjeQzzoAvvoB9+yLQaCGEEA0KFxSSNzev2atwQSEZZHA8x/N3/s5kJlNOeZeupTWcfLJm1x6n5CILEWUyghxl93Ef05nOLdyCHXu7jmmrRnJKCkyZAkuWwA9/GI5WCiGECKW1Gsl+V3Ilq1nNGZzB27xNP/p16lpKwUknuPhgTW8G9N0ldZGFiCIZQY6yUYxiHevaHRy3189+BqNHh/WUQgghOsGOnfnMZxrTmMvcLp3rmed7c/BINj/4qeQiCxFNCREgx1pZoWSSuZiL2Uf4ciJOPhlOPBFqpUSmECKGxVp/HCkKxV3cxRzmsL7+T2dkZMDTT4PXEzoXWeoiCxEZCREgW1VWqKVctdMXnM4QhnAf94XlOv4V9q64wse//x2WUwohRER0xzJvrVEovuIrzuEcVtK5mdTTpsEnqw5RWWmTXGQhokSSlyKotVy1f/NvRjGKW7iFYQzr0nX8K+wdNaSCxYuz+d73unQ6IYQQHbS/cj95c5uXgcvLzqNoZhEZZDCDGSxkIedxXofP/8wCJ4anF2ecckBykYWIAvkfZZE+9OFu7qaEki4FyIEr7JVX7uXdd7OoqVGkpoaxsUIIIQAz4A2ckBeotcl753AOb/AGCoVGo1Aduu7Tf8vilFOyeHt5X1b9dx8fvFdKWrqDygov8x/cxk23j5AgWYgwkv9NFvopP8WHDxcusunc40Z/jWSn005qipcrLqmgsjJLAmQhhIiAlmokhxo9bmoiEwG4gzvIJptf8st2X3f4cLOk5+PzDWpLG3ORPR4fu0vc/O/tA0y9bGC7zyeEaF1C5CDHsud4jh/wg04dG2qFvQE9doGWZUiFECJW3cANPMuz3MZtHVqa+uGHIW9AmblwlF1hGJqaagObDVb896AsQS1EGEmAbLHLuZy1rGU1qzt8bMMKe/bGFfYqK20cfbTC4wl3S4UQQrTXlrItbDiwgT3le5otKDKYwaxgBe/zPi/wQrvPOXQobP2yjuLdaVRWeDlU6sZnaEDhdvuY/+A2CZKFCBMJkCMoLzuPg9UHm73yshsfxTlxchd38Rt+06GRBAi9wp4zxSAn20NRx1dKFUIIESZuw43D5sBus9MnrU/Dy5+T3Ite/Jf/ciVXspnN1FDTrvOeWDiYdVvzufXeUWRmOxgwKJXcIak4U+3sLqnhf28fiORtCZEwJAc5glrKVWvqGq6hL307fP6WVthL6w+LF8O3v93hUwohhOiEppP3DJ8BQJItqcVj0kkH4HEeZz3reY3XyCGn1etMngy33gp//EO5+QTR2TzV4swpfWXCnhBdJCPIMcCBg/M5n3/xrw6PIody6aWQN8Tg+SeK5XGbEEJEQdHMIkpuLml45WblMrbvWEb1HtXmsX/kj4xjHKdzOnvZ2+q+SpkB8r/+ndLwBDEw1aK62mD52/vDdFdCJK6E+BVTKTUVmDp8+HCrm9IiheJBHsSGjUu4pEvnGj4cjh15kHffcNG3v5Nzp/UPUyuFEKJr4qE/joQtZVtwG24MnxFU8cJfJ3ke83iYhznAAQYwoNVzXX45TJ+eTkZGAZUVXh69Zws9eydTWeGl/IiHlctLOWtKPxlFFqILEmIEOR5WbrJhYw5z+C2/xUvXRn0rK7w8Nlfz8fpcVq8ok1FkIUTMiIf+OBLayklWKG7lVo7jOG7lVj7jsxbPlZQEpaXwwAONk7UBqiq9OByKGhlFFqLLEiJAjhfnci796MfbvN2l86xaXsrgAVVs3p6F4dWsWl4aphYKIYRoj6aTtA2fgdfnbTUn2W8iEzmXc3mf91vcp3dvePRRWLPGjc9HQ5qFzwdoKHpXyr4J0RXy/CWGKBRv8AaZZHb6HP7ayEOH2PB4Fa7KVFavKOO0s3rL4zYhhIiSppO08+bmBa2011bKRU96cimX8i7vchzHNTt/Rgb88IfwzaEh/P7BxjQLu11x+JCbynKvLB4iRBfICHKMySSTJSzhaZ7u1PH+x20Oh+KMU8owDDs+Q0aRhRAilrSVcnEWZ/EBHzCWsRzkYMhz/PSnsHo1FL1bKouHCBFmEiDHoHzy+S2/pYKKDh8bWBt57MhScjIr8fnM7UIIIWKPf1GRpguLXLvgWhSKyUzmQR5sVuWof39Yswb27KqWihZChJk8c49BYxnLOZzDH/gDd3FXh44NrI2sNZxwArz6KvTq6eX5J4qZfuUgSbUQQogoa61Osn802c+filHiKkGheJ3XOZdzOchBHuZhbAFjW243/PP1fJYsAe2TihZChIuMIMeoe7iHXezq0jmUguOPNxcNWbW8lA2fuiTVQgghLNCVOskDGUgRRWxnO3vYE/ReSgpkZ8Ozz0pFCyHCSQLkGJVPPk/zdLuXH23JjBnwyiua1SvKyMpOkrJvQggRwzyGJ2S6xbQF01jCEgYykEd5lGoa0+ZuvdWsaPH1V9VS0UKIMJEAOYbVUstIRnZpJPnssyHFUYvXo0lKtsmEPSGEiAGBZeD8JeC8PjOIddgczSbw+dMzNJr1rOcczuEwhwE46SQ47zw475J8brtvFJnZDgYMSiV3SCrpmQ5qa3z87+0Dlt2rEPFIAuQY5sTJ9/ge93BPp8/h9XgpnLCD1DQ7AClOu4wiCyGExQJTLvzpFmP7jiXJbtZJrvXW4jbczUaTz1pwFgtYwMmcTCGFHOEIAH/+MwwYACv/KxUthAgHCZBj3K3cyhKWsIUtnTp+1fJSSnY5efm1XADsdoXXo5n/4DbpLIUQIgaEGk32aR82ZWs2mlziKsGGjUd4hId4iGyyG1ZfvfBCWPaOlooWQoSBBMgxrgc9+Bt/6/TiITuLq+ndo47N29IpK/VRWeGlptrgm101kmohhBAxINRocrI9GafDSa23lhpPTcNosn8k+fQFpzOFKbhwUUABn/Ipl18OG7cN4N55BUGpFhlZDtCwcnmpDIwI0U4SIMeB8ziPWmrZytYOHztrdj4P/XUMZ59j54QzxjR0mn37OyXVQgghYpzWGqUUYE7gM7TBvsp9fLDrA/Lm5nHhggt5kAeZzGT6X/E/Nm+GTz+VihZCdJUURowRhQsKg2pk+vmXHV3KUl7jNZaytFPn/8EPoLy8sdN0Ou14qrysWl7KudP6d7X5QgghwsBfL9lfJ1mj8a8PopQCTUPNZH/KxUVcRA968NPkn/KHP63F4UhqWDTKn2ahlMJnwLqPj3DhZblW3Z4QcSOuA2Sl1EXA+UAW8IzW+h1rW9R5Ja6ShuLwTbcDXM/1PMIjvM/7nM7pHT7/ZZeZOWmP3lNGijN4wt5pZ/WWIvJCiE7rTn2x1YpmFgGQNzePPml92HBgAw6bgxqPWfJTo6nx1KDRbDiwAcNnmKXgsvP4dOanMB2WlL/L9O+dgzOl+cIhlZVeKiu80ucL0QbLUiyUUs8qpQ4opTY22T5ZKbVFKbVdKXV7a+fQWr+qtZ4F/Aj4biTba7Vkkvkdv2MOczp9jplX17Jlezp2u/m4zm5XUvZNiAQnfXFs8k/c80/a02i0NoeSW0q5OGruUZz20mlc//AWLrjnE1YuPyhpFkJ0kpU5yAuByYEblFJ24HFgCjAGuEIpNUYpNVYp9UaTV9+AQ39bf1y3dgVX8Hf+3unjbdSyaVs2lRXehpfXq1m5XMr/CJHAFiJ9cczxT9wLnLSXmpSKQjXsEypQXrt1Lc6sZ/nohaN4eeNaDJ8OWjjEn2YhhGidZQGy1roIONRk80Rgu9Z6h9baDbwETNNab9BaX9DkdUCZfg+8rbX+NNR1lFI3KKXWKKXWHDx4MLI3FWF27CSTzM/4GT58HT7+94/1Zs+BHtz1aAH3zjNfhWf3oa7WJ6PIQiSoaPXF0L3642hpbSTZTymFQjWUg8vNtpF+zFL2efs0q2ahFA1pFkKIlsVaFYtcCFo2bnf9tpb8FDgbuEQp9aNQO2itn9RaT9BaT+jTp3mOb7zJJJM1rOEf/KPDxw4dCmeeCbt3m59XVnhlCWohRChh74uh+/XH0eAfST518Kn0z+iPXdkbFhPxB8o+7WvITfaXg6uZ+Fv22x9n0fIv+MbYiw8fVZVe7A5FVaXB0iV7rbwtIWJeXGfpa63/CPzR6naEg3/mcqjtgRSKOczhWq7lUi4lmeQOXefllxs/looWQohw6E59caxqbfKeri9z4f/bY3jQ2V+zzlXHM89dSH6/dHZVHibZSDP30LBm1WEuuWqwRXcjROyLtQB5DxD4P3ZQ/bYuUUpNBaYOHz68q6eKGH/n1x5ncAancRpf8iXHcmyHrnPkCFx1FfzjRXP0WCpaCCFCiEhfDPHRH8eykGXgmvBv8x0ezIr/jWfNTyZyiesZvD3KyTgwAG3TVFZ4+Xp7JcOGZ0S1/ULEi1hLsfgEGKGUGqaUSgYuB17r6km11q9rrW/Izs7ucgNjxfM8z1jGNiwx2l7Z2bBtGyx86gg+QwdVtJAlqIUQ9SLSF0P37I+jKVTKBRA0ea/B4I8gczf91vwOX7lGH3Tg9rrRGrSGhY8XR7fxQsQRK8u8/QP4EBiplNqtlLpOa+0FbgKWAZuBl7XWX1jVxlh3D/fwCI906BilYPp0eHupHZ+PoIoWsgS1EIlH+uL4FBgoB+YlN3PaQxz+bBp2kkh256ABw/ABmtIDbr7eXhnNZgsRN1TT2bDdUcAjvVnbtm2zujlhs4UtfItvsZWt9KBHu4/75BP405/guecat5mLiGzBblMYPs0td42UVAshIkQptVZrPcHqdlihu/bHVgrMS/YYHnPVPF1f6cinoLI/J3gmM9o1DbtOwalzUEphNxz07pvCPY8VWHsDQliopf441lIsIqK7PtIbyUgu4iIe4qEOHXfiicHBMTRO2EtKtsniIUKIiOmu/bGVAkvBNSsDZ9PgqKPy4+9jJ4lUowfKp7AZdnT9KPL+vbXWNV6IGJUQAXJ3did3dniiHsCiRfDkk+bH/nJvTSfsSS6yEELEvjbTLZIr2LppBB+7P6LSsY8a+2EMvHiVG42Ph5/+0JqGCxHDEiJAVkpNVUo96XK5rG5K2A1mMN/lu3zCJx06LjsbnnrK/Ng/eixLUAshIq0798dWK5pZ1GzyHgAOD5w8jyOffLdxFBmFXScBioptqQz/oCBkRQwhElVCBMjx/kivcEEheXPzmr0KFxQCcJjDTGEKO9jR/nMWwo4dsGsX7CyuDpqwV+7yUO7ysmNbVaRuSQiRoOK9P451TUeTG4x/kn2HUvg89Q0qHfs4nLIDn/JiKDc2w07fDWcz9NOhlrVbiFgjs7DiQImrhD5pzVed8i8s0ote/JSfchd38TzPt+ucSUlw2WWwbh3Mmp0f9N6yJfv471v7yR+R3uW2CyGEiL6imUUULijkg10fYGgDnOVw00j67Z+DDTuZngEobcdcpNrG0avO48PaeaQsS+HEgSey8gcrrb4FISyVECPIieAX/IJ3eIctbGn3MX/5C0ydGrxNlp8WQojuwZ9y0TCSbDN454PBvNznZ1TbD3EkudgcRaYOp5FF1oHBuMe7WTVpFac+f6q1jRfCYgkRICdCzlsmmaxmNUdzdIeO+8lPoDQg1ViqWQghIikR+uNYEhQkKyCpip4rH8Wm7aQYWaAVWoHCwen7fwOrgSPw4VkfkvJoSkMqnxCJJiEC5ETJeRvGMF7kRT7io3Yfc+AAvP66+bFUsxBCRFqi9MexJChIPvVRytZ/B+1NCZqsp1AMqBlHijcbXgd2gHu0mw92fSBBskhICREgJ5JaarmN29o9G3n6dFi82Py4aTULANdhD8vf3h+JpgohhIgSf5Cc0v9r3Mc8x6e+D5tP1sPOCWXXmAcsBz4BY6jBipoVEiSLhCMBchzwF4Fv+srLzmu27w/4AfvYxzu8065zn38+bNwIXm/zahaHSt2463ysXyuPQoUQIt4VzSyi9re1TPrJP8jLzEDpxsl6Nu0wJ+uVnxd8UDpwNawwVuC8zymBskgYUsUiDhTNLGr3vg4c3Md9vMM7nMu5be6fkwPbt4PNFlzNwr/0dE6PJAxDU1nhlaWnhRCiGyiaWcSA55bg6v0qF6Wfh89xgGzPYNCQ6utB75oxlKZuMnfeCNQCV0Ddv+tY8bU5mtyRn0tCxKOEGEFOtEkhF3Mxj/Jou9Mstm6FRx4J3iaT9YQQkZBo/XGs+uvN0xj02V9ItWWQ4jMn6ylsKGycs/fe4J23Ay8Bh8xPV+xcQfaDkkMuureECJATbVKIQrGTnZzMyXhpe4JdTg7cfz/U1Zmfh5qs9+H7pTw7f4dM2BNCdEmi9cexaupUKK/Q1B4eTprRq74ash3Q5HiH0LtmTPABuwEXMBWYCOV15djusUnKhei2EiJATkSDGUwaaSxkYZv79u8PxxwDy5ebn4daerqq0mD9WpeMJAshRDdgs8ETTyYz9vgc+vZJpWePNHzKiw8DUM1Hkf2KgBOBM0GjWbFzhQTKoluSADmOtbYEtULxAA9wD/dQS22b55o+Hd57z/w41NLTtbUGdpuSsm9CCNFNTJsGORlHqKhUVJR7UNhaH0UGcxR5AXBU/YvGQFmCZNGdyKyrONbWEtQnczI/5+dUUIETZ6vn+tnPwFH/3RBq6ekP3islLd1BdZWXVctLOXda//DchBBCCMt4k/P5fFsN3zr+KzyeJKoqvLgND0rb+c7eObyYf1Hzg6oxg2QDyAX2mR/7g2SZwCe6g06NICulcsLcjohK5Ekhv+AXuHFTQUWr+yUlwQsvwGefBW+XxUOEiG3SH4uu+P734aNPkik75KCi3INPa2zaHC3J9uaS5R4c+kCj/u8TgSuBZPNTmcAnuovOplg813SDUmpQF9sSMYk+KeQO7uBRHm1zvx074Pnng7cF5iMbhubIITdej1S1ECKGSH8sOi07G35yo509hwbTs1cyGRlmcKyUOW3vnP2/a/0ESzCrW/wASDM3yQQ+0R10NkAuUkrN8X+ilJoIvBueJolwu4M7+BN/4iAHW91vxgxYtAh0QHW4wHzkQ6VuqqsMaqoNdhZXR7jVQoh2kv5YdMnPfw7DhrjwejUV5fVPB+t/DvSvPZbe9lZ+39LAm8AGwB64WfKSRXzrVICstX4EGK6UmqGUugx4BpgR1paJsBnGMK7kSh7kwVb3KygwUy2+/LJx26zZ+dw7r4Db7htFZraD/gOdZGY7uOLaIRFutRCiPaQ/Fl01YADc+9BAeg7sh7KBI0mBWcQIreGPA1ei79JkpWS1fJLVQAVwEdC3cbMEySJedShAVkq9r5T6g1Lqe8CjwN3Aj4BCrfXmCLRPtKIjS1DfxV3cwi2tnk8pWLcORo9u/p4sHCJEbJH+WIRTSQncfX9PvB4znS5wnamPV5grhLhudzFpyKTWT7QN+D4QkLoseckiHnV0BPlG4DNgAvB7zP8CmcB9SqlZYW6baEPRzCJKbi4JeuVl51HiKmlW+m3GghkoFE/yZJvnfeih4M9lop4QMUn6YxE2Y8ZA/tA6ivf1RSmwBaRLuN2ar7dXAubPnVaD5C+AxcB3gYCYuLyuHHWPktFkETc6FCBrrTdqrZ/XWv9ca32G1roHcDnwHjAsEg0UHeMv/db0VeIqIZVU/o//40u+bPH41FSYOxe2bWvc1nThEADXYQ/L394fwTsRQrRG+mMRbo//NZ31W/rRo2cymVlJOByqIVBe8KevG/bzB8kKFfpEO4AnMGsmN8nKkJQLES/aDJCVUklKqZuVUkuVUh8rpV5VSt2klMoE0Fp/pbX+l9b6N5FvbudIWSFTDjn8kl9yB3e0uI/NBhddBIsXN25runDIoVI37jof69cm9r+nENGglOqjlPqnUuqQUqpGKfWhUmpqqH2lPxZdcdpp8MMf7MUwdEPJN199ObeyUg/79zYuOlU0swjfXb6W85IrMCOMq4GTg9+S1fdEPGjPCPKfgD9gzk/9BLP64b1AiVLq4gi2LWykrFCjn/JTNrOZUlrOIZ4+HV5/vfFz/0S9ppP1DENLmoUQkTcPc+rTUsySbunAq0qpK6xsVGdJfxy7lIJf3TmYI3V59OjZWPLNP1D8/JPFzY5pNS/Zh/kdewLw7eC3pMqFiHXtCZAvBe7UWp+jtb5Ra30xZq7bo8DflVIXRLSFIqzSSGM96+lN7xb3OfNMeOON0O/JZD0hou5s4Bda6yu11j/UWh8LPAk8YHG7RDeUkgL/ftXJjpK0ZiXfirdVhxwUKZpZ1PJIcjnmqnvZQErztyVIFrGqPQGyBlYGbdC6Umt9P/AYcF8kGiYix4aNH/JD/sf/Qr7vcJiLhvyvydtNJ+s5kmy88/q+oMduQoiw6w2sabLtQWCwUqp5yRohusBmg/vmJPNl8cCQJd9e/9eekMf5R5JD5iXXAIswnz+fSVC9ZDCDZMfvHGG8CyG6rj0B8lJgegvvLQNGhq85oqvaW/rtdE7n1/waHVjLJ0BxMdx/f/C2ppP1aqoN6mp9vLxwVyRuRQjRyNfk892YYUsPC9oiurmrroLtO2xUVia1WPItFH9esl3ZQ+/gA3oCV9FsNNnQBuqeFib9CWGB9vzKtgO4WSlVCTyitT4c8N4kzKIuIkYUzSxq136Xczm/5/e8xmtMY1qz9889F665Bg4dgp49zW2Bk/V8PnMSh80GO7ZXUlnhJSNTRgCEiJDHlFIfY/a3XwBb67dLRCHCzumE+3+7h292QrkLlJ2GyXput2b/3lr6DXC2eLz3Ti/ZD2ZTXlce/IYP+DcwBXNp6mcwR5UDqHsUdmXHe6fMbxHWUlqHHkFs2EGpMhpHKaqAIuAb4Kj61/la642RbGS4TJgwQa9Z0/RJZeL6D//hCEe4hEtCvj99urn89Pe/3/y9ZUv28cF7paSlO6iu8nLqGb05d1r/CLdYiO5BKbVWaz2hnfs+BhQAYzHXKPN32gpYDnwArK9/bdNtdeoxQvrj2HfJ9CoG9dxNakod5S6v+R2nYdjRafzyrlFtHh8ySPYbDOzCHKYLEQtnpWThul0qnYjIa6k/bnPIT2vdSyk1GDgu4FWIGRzbgBVKqQ3A58DnWuunw9pyETFnczYazR72kEtus/f/8pfG0eNAgbnIhqGprvLy4fulnHZWbxlFFiLMtNY/93+slOqLGSiPBY6t//sWIBUzcK4BMixopuiGyl0ePvy6FycW7DQ3NJms11Z/7w9wQwbKuzBrJM8EXgKalNX3Lyyi74qL3/dEN9SuaEZrvQvz27mhtoFSKg2zc/YHzcdjLjApAXIMKVxQSImrpNn2vOw8imYWsZa1XMZlfMmXJJMctE///vDiizBtGqSnN25vyEV2KlxHPNTW+ACDVctLZRRZiAjSWh8A/lv/AkAppYARmAFzgUVNE93Qs8/lMHqUj2NH7iEt1YfXYwarWsPSJXu55KrBbZzB5LrdheN3DgzdJJ+iHHgXM3J4GdjZ/FgJkoVVOrrUdAOtdbXW+iOt9ZP15d++pbXOCWPbRBi0trIewAQmMJKRPMVTIY9fsACWLQve5s9FLnd5qCj3oBTU1hrs2FYV6dsRQjShTVu11q9ore+2uj2i+xg0CI4ZeYQDZWkY3uAg9aNWJuuF4r3TG7oU3CbMvOQzWz5WlqgWVuh0gCy6jznM4X7up4rmAe706cGr6kHjwiGFZ/ehZ69kcoek0bNXMvkj0psdL4QQIn5dd/URCkY3L+VZU2V0eKEo1+2u0KPBXwN/A5KBFlKbpV6yiLaECJBladPWHc/xvMzLpJLa7L2LLoI33wS3O3h705rIKU47q1eUycp6QohWSX8cX2bNzmdX6RA2fdWr3TWR26Lv0qFLwaUB5wKnhj5OgmQRTQkRIMvSpm37Ft/iRV7kEMGPzQYOhFWrICkpeP+mNZEBjhzyMO/+rRIkCyFaJP1x/ElWh/n0i354PO2vidwW753e5ktUHwGexZzVdHbo41bsXEH2g/K9IyIvIQJkYdpStoUNBzaw4cAG9pTvIW9uHnlz8xp+Iy+iiN/z+2bH9ekD//1v8LbAmsiVFV4Olbpx1/nYu7tWlp8WQohuZOJEH+lpPor35GALGPj110TurJBLVFdgLk29t+Xj/BUuhIgkCZC7ucCV9Wo8NQ3bnQ5nswl7d3InT/EU3/BN0DkOHzZrIfsC1vLy5yHfO6+A2+4bRUamHbsDbHbFB++XyiiyEEJ0EzfcnM+Pb6jFmZpEZlYStoDI4R9PN6+S1BEh85JrMZfDyQO+S4v1ttQ9SkaTRcRIgNzNFc0souTmEkpuLiE3K5exfccytu9YRvVuPhNiEIO4lmt5kReDto8YAb17w+rVoa+xankpVZUGqj45rarSK6PIQgjRjfTrdYTRw49wqMwIGizZsa0qLAMiISfv7QLqCLk0tV95XTmO30n9fRF+EiCLIHOYwy3c0mz79OmwaFHz/SsrzEVC6uoMlA1sCtx1mmWv7e3SozchhBCxY9bsfGxpg3hnVT52R2N6g2HA8rf3t3Jk+zULkn3Aq5jpFue1fJyhDQmSRdhJgCyCJJPMx3zMbGYHbf/xj+Gmm5rvHzh6rFAopfD5NHW1mpcX7opSq4UQQkRan5zD7C9N5WBZ8HDumg8Oh+0aIStcLMNcpiwF6BH6OEMbkpcswkoCZNHMGMbwT/7JetY3bBswACoqYPfu4H13Flfj8fjw+TRerw+P14euf/y2Y3ul5CILIUQ3MXJMGhPHHWLd5n5B2w+XucPa13vvDHEuD2ZO8kyglQVbZSRZhIsEyAkkcMJe4CsvOy9ov0wy+TW/5v/4v6DtCxbA000WEp81O5+5C47nT8+fwJ+eP4Gplwykd99khgxLIysrSXKRhRCim9hZXM2Yow4wsG9F0Hafr/M1kVsSciR5K/AWZk7y0NDHyUiyCBcJkBNI4IS9wFfRzKJm+/6IH9GPftRR17Btxozmq+oFClw8xDA01VVmfrKMIgshRPybNTufAQPtjD7qEIdczqD3ulITuSUhayV/CbwC5LR+rATJoqvkWUSCKlxQ2FDeLVBedh5FM4tIIYWneZrDHCaZZBSKU06Bffvgq6/gqKOan7Nh8RCnwnXEQ22NDzBYtbyUc6e18kxMCCFEXOg30Mn+A4rF747gqgu/ICXZABprIvcb4GzjDB3jH8AJCniL6/8+BjMv+dPQx6p7FJOGTAo5CCREW2QEOUGVuEoa6iAHvpoGzZOZzFKWAmC3w6uvQt++oc/pXzyk3OWhotyDUlBba7BjW1WE70YIIUQ0zJqdzzlTshg6yMWmr3oHvff8k8URu27IlIu9wLfqXy1YsXOF5CWLTpEAWbTqVm7lN/wGH+bMuxNPhC++CL2vf/GQwrP70LNXMv1zU0lKspE7JLwjCkIIIayzfq2L40fv5/Mv++A1Gkd2i7dVRzSlznunNzhIPoS5NHUBMLHl4yQvWXSGBMii1SWoZzADBw7+xb8AczLGlClmqkUogXnIlRVe3HU+Vi6XPGQhhOgu+g10kjvAw5RJX2O3NdYu1jp8NZFb4r3TG7w8dSXm0tQbACetRjUSJIuOiNsAWSk1Win1V6XUK0qpH1vdnnjmNtw4bA4cNgd2mz0o3UKheJqnmVj/63lyshkgL1kS+lz+PGQwV9RzOBQ11UbEO00hhHWkP04ss2bnc9TIDPr1qWJbcU90wPoen64OX03klrhudwUHyXVADXAacDmQ1PKxEiSL9rIkQFZKPauUOqCU2thk+2Sl1Bal1Hal1O2tnUNrvVlr/SPgMsz/FiJCjuM4kkjifd4HWl5VDxrzkA+VuvEZ2lySVJuP5IQQsUf6Y9EZ+SPS6dEzmQ1b+/D17uyG7YcPeaLyxNB1u6t5TvL/gCrg+5ijyS2QIFm0h1UjyAuByYEblFJ24HFgCjAGuEIpNUYpNVYp9UaTV9/6Yy4E3sSsjCg6ILAmsuEz8Pq8eH1ekmyhf/X+hm+4iquooYYpU+DnPw993lmz87ntvlFkZjsYMCiV3CGpDBiUimFoSbMQIjYtRPpj0UE7i6uprTEYN2Y/n23q1zCK7PXosNdEbon3Tm/w8tQ+YAmwkxbrJPtJkCzaYkmArLUuwkyvDzQR2K613qG1dgMvAdO01hu01hc0eR2oP89rWuspwPdaupZS6gal1Bql1JqDBw9G6pbiTmBN5NysXMb2HUuyPRmPz9MsH7lwQSETmciJnMif+TMZGXDKKVDSvEocEFDuzW52QHa7wmdoWTREiBgk/bHoDH+axVFDjlDrdlB2JLXhvVXLy6I6IBIUJAP8B7Ne8ligV8vHSZAsWhNLOci5wK6Az3fXbwtJKXWGUuqPSqknaGXEQmv9pNZ6gtZ6Qp8+fcLX2m4oMBc5MB/ZX/rtPu7jz/wZA4NFi+CWW0Kfx59mUVnhpdzlYd83tXi9mp3F1VG8GyFEF0h/LNqUPyKdnr2SuXTyFnr3qGnYrjW88a9votqWZkEygB34ATCw5eMkSBYtidvigFrr94D3LG5Gt5CXnUeJqwTDZwRtb5puMYYxfM7n2LEzdSrcfDPU1EBqatBuzJqd3/DxsiX7+O9b+yk8u48sFiJENyX9cWLaWVxNdZWX5CQf67f0YVD/Cnpm1wKwcnkpF1w6kIzM6IUZ+i5N4YJCVuxcYW5YB1QDV2I+A9kd+jhZUESEEksjyHuAwQGfD6rf1mVKqalKqSddLpkoFoo/3cKfauF/jeo9qtm+6aRzDdfg632A8ePhnXdaPq+/5FtWdhKrV0T3kZsQokukPxZtmjU7nx69krHZwO2xs25zv4b3tDYHSKKtaGZRcIWLrcA/gIO0GvE0BNVC1IulAPkTYIRSaphSKhmzWMtr4Tix1vp1rfUN2dnZbe8sWqVQZJHFHOZw770wqnkM3cCfi5yUbJMcZCHii/THol36DXRisysKRhzk693ZVFY3PnkseveAJQMjrttdTBoyqXHDHsxScN8FJrR8nKRbiEBWlXn7B/AhMFIptVspdZ3W2gvcBCwDNgMva61bWLNNRFLgwiFNJ+sB/B//x/M8z6DTSkhPB4+n+TkCFwwBSHHaZRRZiBgk/bHoilmz80nPcJDqNBiVX8ae/ZkN73m9WDYw0mwkGeBt4BSgsOXjJEgWflZVsbhCaz1Aa52ktR6ktX6mfvtbWuujtdZHaa3vD9f15JFe+/hLv9V4aoK2Ox3OoMl6/ejHr/k1X/IlF18MRSHStgIrWRiG5sghN546H/Mf3CZBshAxRPpj0VVz5o/lmHHZnHbCHkYOCy6IsnljuUWtClEr+Qjm0tQjgVamxEiQLCC2UiwiRh7ptU9HcpF/yS/5Dt/hvIvcLF7c/FyBlSwOlbqprjKoqPDyza4aSbUQIoFJf9w95Q5xYrPD1uKefP5lY4WSA/tqLWyVWSs5SBXwNLAPGE6LUZAEySIhAmQRGf/gH6yacQuLF2OumBdg1ux87p1X0LBoSO++yRheTc9eyZJqIYQQ3cz6tS7Q0DO7hs8298MwzADTdcjL/r3WBsn6Lh2cbuGvCHcCcAUtLk2t7lFkPyi/yCWqhAiQ5ZFeZExnOptGLub7dxTjdofex59qUeHyojWUl3tkwp4QCUz64+6p30AnySk2eveooXdODVuKeza897c/f21hy0yu213NayX/C6gArgZSQxwElNeVN8y/EYklIQJkeaTXeYET9vyT9fwT9lJJ5Q7uYM0PZ1FR0fxY/0Q9u11RW2sOMdfVamw2JaPIQiQo6Y+7p8CSb8eP2U9tXWP945IdNTHT3wflJGvM2izrAaOFAzBLwEnKReJJiABZdIx/sl7TCXv+yXqBE/au5VrO+vznnHaaRjf55dw/euw6bJa5UPX9i+uIW0aRhRCim/GPIuf2q+SEMfsxAlLvli7Za13DAjTLSQazqKEXc0GR3i0fK0FyYpEAWTTjn6zXdMJeqMl6SSRx23GTOVJXyxebgiPkncXVeL26YfTYH0DX1Wrcbp8sPS2EEN1I4CjyvrJ0XvvviIb31qw6bGHLgum7dPN0Cx+wEXNp6hYXVTeDZEm5SAwJESBLzltkKQV6+iLuX7Q5aPus2fkUnt2H3n2TGTIsjSHD0sjukYRSMHBQatCS1EKIxCD9cffWb6CTrJwk+veqoqommb0H0wGoqzViJs3Cr1mQvB4z5eK7gLPl41bsXIHjd9FbQltYIyECZMl565qmC4f485H9v0XbsHHnzMGsOvbPGE0SuQLLvZW7PFSUe1A22LqpwvKZzUKI6JP+uHvzD4ykptkYN3o/n20yl5/2+axbNKQ1zYLkbcBfgFpaTbcwtCEVLrq5hAiQRde4DTcOmyPoZbfZG/KQAW46bhKDz9/Av2reCDrWX+7t3nkFFJ7dh569ksnMSkJreHnhrmjfihBCiAjbsa2K2hofo4aV0b93FVqD16stXTSkNc2C5BogBbgKOLHl48rrymUkuRuTAFm0Ki87D8Nn4PV5qfHUUO2pptpTjcfwBC1BrVCMvnkZ3/z1wpDn8Ve0cCTZqKr0YrfDts2VMooshBDdTO4QJ8oGycmaE47ZT2V1MmD9oiGt0Xfp4AoXdcAC4CTgjJaPM7Qhk/e6KQmQRauKZhY1TNRLsieRlpRGWlIaqUmp2G32oIoW06c4WbD4MM/xXLPz+Cta1FQboMFms6G1llFkIYToZvyLhvh8UOe28fLbo6iqToqJRUNa473TGxwkuzCXpk4H7C0cVE+C5O4nIQJkmRQSXrXeWmo8NbgNd1B95Af2nE3x+mxu2f8QlVQGHeOvaFFR7kFrMAyNUood2ytjbuKGECJypD/u/vzl3gBSkn0cPfQQ67eYy08/O3+HlU1rU7MguRp4EzPl4ju0GihLukX3khABskwKCS+tNbp+rU6P4cHQBvsq97F633v0mzqPk6vPYi5zg47xT9zo2SuZ3CGp5A5JpX+uE4Vi+dv7LbgLIYQVpD/u/mbNzic7p3H95uNGHWDzjt64PTZ2F9fG9CgymEFys7zkOiAbs1ZycujjDG1IkNyNJESALLrGv3CIPxdZ09hxKKVQqIaJe57x8/h97s/YwIZm5wmsaFFZ4eVQqRt3nc98HCeEEKLbsNlVw+JQWRluzjn1a2w282fH0/O+srBl7RcUJBvAK8AhzMl7LZCc5O5DftURbSqaWQRA3tw8+qT1YcOBDXgMD0o17wS0N5kzBw9n69Z/4s324gj4Fguse7x/by0P3fElffulYBiaygovGZny7SiEEN1Bv4FO9u+tRddX/hw8oIKd32SR26+Cb3bV8fX2SoYNz7C2ke2g79I4fufA0Ob8Gd4E+ta/mYI5shyCukc1H4UWcUVGkEW7BY4ka7SZaqF1UKCsHG4mTIAlb3kYzWj2sCfkuV5euIu6Wh+1dYYsOy2EEN3MrNn59BvgxBYQZaz7si9bi3sA8JeHt1nUso7z3uklKyWrccMBzCD5x0Cflo+TkeT4JgGyaDf/EtS5Wbkk25NJTUolNSkVoGHS3p7yPazOvI0f/34ZtRtquZd7m51n/95atm2uxO6AqgoDR5KN1SvKZLKeEEJ0I/0GOrHZG4PEE8bsZ93mfmgNVZWaTz86ZGHrOsZ1u6t5kPxf4GpgUMvHydLU8SshAmSZNR1eTWsj+7SvIS/Z0AYVQ/+Op8cX6FWaV3iF7WwPOv7lhbvQWmNTNjRQU+2VUWQhEoT0x4lj1ux80jMaU+dy+1XgsGt27jUDzQXziy1qWee4bm/yPbsBeBWzVnIrVuxcEaEWiUhKiABZZk2HV9HMIk4dfCr9M/pjU+a3kKr/Y1M2krIOo759J7ZqB0/yJPaAujiVFV6+3l4FCgyfmaJRUe7F69XsLK626paEEFEi/XFimTN/LJlZZpCsFEw5/SsG9zdX1PP5YNP6+PpFqdmCIl8B/wYygTEtHyfLUsefhAiQRfi1lG7hdDgBMNZfxqGXH2QGM0glla/5GjAXDMnMcjBoSBq5g1MZkJtKUpKNkyb1DJrEJ4QQonsITLPISPNQvCebA4fM9Lyn5sV2XeRQvHd6mTRkUvDGZMw6yS2MJpfXlTesPCvigwTIIiJsQ1ZTs/EcvF74F//iJm4CWi71tu7jIzz/RLHkIQshRDczZ/5YHI7GILmqJom1GwcA4K7VMV8XOZSimUXBVSrKMFfdGw+cHvqYna6drNi5QmolxwkJkEVEqJyd2Ht8w8qVcAM38AVfsJKVzJqdz73zCrh3XgG33TeKzGwHffqnUHbQzedrDkseshBCdEPpAWU8R+WXse9gOoddKQA8/2SxRa3quqCR5HJgAVDS+jGGNiTlIg5IgCy6LNmejNfnDXoZPoNh33mdujpIIYV7uIe7uTvouFXLS/EZmtpqA69XY3iRahZCCNENzZk/FqfTDDmSHJqCow/y1S6z5NvOr+J3/knRzKLgILkGKAaOAS6jxaWpy+vKI9420TUyzi+6JC87r8XtRXf+BgCt4Sp1Fd/m2w3vV1Z4Wb2iDEeSjUOlbgDcbo27zmDV8lLOndY/8o0XQggRNaOOzaJ4exWuwx7GH7OvYaU9f9pdvC4W5V9MK6ju8ZfAKMxV914i5IIi6h5FVkpW8+oYIiYkxAiylBWKHP9kvaYvf4cxdSqsXQt27PSgB7/hN/jwNYweV1d68flo6CirqwwZRRaiG5P+OHHNmp3PSZN6AmCzwe79mXy2uS9aw9Iley1uXdcFVbcwgEXAQVrMSQZzJFlykmNTQgTIUlbIOsccA4sWmR+nkcZ/+S8v8zI7i6vxes0Sb2COMgN43Bqvxye5yEJ0U9IfJ7Z1Hx9p6O+zM+r4bFN/3B4bH6+Mn0VDWuK90xscJGvgLcwFRbKBnNDHGf71uEVMkV9bRFgVLiikxFXC/sr91Bl1UHYivPo3HkgZg13ZseXbmHnhTMpnl7N8SRlLl+zFMMBn6IZOs6bKkJrIQgjRDdnsCpvNTKvIynAzqH85m7b35oRjDsR1moWf905z0Cf7wezGPGMDGAKcDfwdcxW+JrIfzJZUixiTECPIInpKXCUcqjlkBscAA9dA/3VQl45Go4oVar9iAxvYWVxNRmYSmVkOlA1sdvMX7qHD06UmshBCdEP9BjpJcTaOsh4/ej9HKlLw+ehWTw6bLU29AXgH+D5msNyEpFrEHgmQRdgULihkT/keKt2VjRttGi65ElKq0PVDxH3e7MMJnMD1s4dx77wCCs/uQ89eyWRmJaFoTLcQQgjRvcyanc/Qo9IaPu/Ts4YzJu5Ca9i8sXtVdmg2IvwFsBhICr2/oQ3UPUpKwMUICZBF2JS4SrDb7ChU8BulI+DF14I2PciDPMIjQdUsqiq92O2wbXNlXBaOF0II0Tavxxf0uasihUXvHM2Bb7pfv6/v0sEjyTswl6c+BTgu9DHldeWoe5SMKFtMAmQReT12wK5T0EcGN2yaxjQe4iH+s3wXPkNTU22ABpvNhtaalxfusrDBQgghImVXcU3Q51kZdRg+G59/kdYtKxi5bnc1X5p6K3AGZqDcAkMbEiRbSAJkEXl2A0a+Dl9e1LBpNKO5gAt4r3hDfTULD1qDYWiUUuzYXtktO0ohhEh0o47NwpnaGH4oZeYif7apH4te3G1hyyKnaGZRcIWLMsxV944Hjm35OAmSrSMBsgirZHsymhBJxAUvgTsDhWpYXORu7uaY2fsbcpBzh6Q2vLKykrrVhA0hhBCmWbPzyc4JTsQ9ashh+vWu5sP3DlvUqsjz3ukNTrcoB57FzE3OhKbZiX4SJFtDAmQRViN7jSQjOYNkezLJ9uSGfGT7iP+Scuaj9MvoR4mrhMIFheSRx/VczxfF+xpWUip3edi3p5byIx52bKuy+G6EEEJEgr/cW8PnNjj1+D24PTY2re++5c5ct7uCR5JrMcvAnQ58lxaL7xrawHmfk7y5eRQuKIx8Q0ViBMiyclN05GXncbD6ID1Te9I/oz/9M/pjUzbSk9IZ138cfT5/APvGq+iT1ocSVwkAhznM/bO/zXXzMhoqWuj6Mej8EenW3pAQIuykPxZglntLSg4OQXw++Odbo3nonu6ZZuHnvdPLpCGTgie0vwW4MZemTgl9XJ1Rx07XTlbsXCFBchQkRIAsKzdFR6hlp3OzchnVexQASdmlHPr43KBjetKTG7mRu7mbygovH75fis+n8fngg/dLJQ9ZiG5G+mMBZppFcpMA2WaD0UeV8cnn/fh6e2ULR3YPRTOL8N3lQ9+lzRFlH+bS1HuAPm0fv2LnCikHF2EJESALa20p28KGAxvY3f9JXFvHsX7X1+wp39PwG/At3MJbvMWryzdRVWmYv1VrOFTqZt79WyVIFkKIbshmb550O/bog3y1M4c/3F9iQYusEbRE9bvAbszqFj1aP65hpT4RERIgi4hzG24cNgdG8iEY/W88B4ZiaIMPdn1A3tw8pi6YytKK99j2vo26OgOUWc3CZ8A3u2r539sh1uUUQggR1+bMH0uKMzhITnV6OeuUEqorvAk1OBIUJAPUAdcA/Vo/TkaRI0cCZBE1Wmts06/HlrsOhcJuszfkIx9c3pMjlTVoND5fcBWMd9/YJwuHCCFEN9SzdwqqyUDysEEuKquTue/2LQkXJDf4FFiKuTR135aPkVHkyJEAWURUXnYehs/A6/Oawa9X4fv7q2hvEm7Dzdq9a9np2sm8t1+gttaD1/DhM4LPYRjIwiFCCNEN9RvobBYgA6z9oj8rV2ezdMne6DfKQpOGTCLFnmKOJm8G/g6U0mJ1CwDnfU6ZtBcBEiCLiCqaWURuVi7J9mRzg92A2hzY8e2g/d4e8CsWjDqDtecswFBeNMFLkW7dVCGjyEII0c3Mmp1Pekbz6O/40fv5/Mu+/OfNgwnV9xfNLKL2t7UNlS7YizmB7xrMRUVCqDPqpLJFBEiALKLCbbgbS9qMWgxfTm+2j702nYL3LkaFWGfE55NRZCGE6I5CTdbr26ua7Mxatpf05IUnE2fCXqCimUWNS1QvAiYBp7W8/8qdK6PRrIQhAbKIOH+aRcMKe6MXw+FhzfYbc2Q6Dk8K5rdl8w5Tlp8WQojuZ878sTgczfv8M0/aSf6gI+zYWpVQo8iBGoLkQ5hLU48CWpiXp9EyihxGEiCLiPOnWWQkZ5gbehTDD85ptl/vulHU2l1UJe/H7axE+b87FSgFaWkOWX5aCCG6ofRMR7Nc5KwMN2VHUtlXmpawo8hg/gwFoAJ4BnABBYSM4GQUOXwkQBZRM7LXyMblp78+HT68Oej9dwbexgv5U3npqOl8M/wTtK1+tp4GraGy3MvO4uroN1wIIUREzZk/luNOzGkWJJdXJrN6XS5fb6tK6CeIQavu2YCxhFyaWkaRw0cCZBFVyfZkM9UirQxWz4Ym+cZKKRSKrX3+gaF92Ov/89vtZpB84WUDo99oIYQQETdrdj5ZOUlBQfJReYcpr0xm38G0hH6C+K0h3yLFXr8GtQ/4J1CNWQauydLUMoocHhIgi6jIy87jYPVBeqb2xK7sJA3YCjYv7BuHqv8DmAuK+AxO/fwX2LUdXR9Aa23WUZaJekII0X35R5IdSebPBLsNjht1gM07erFtc4XFrbOOv7pFw6Q9H7AE+AxwB++rm448iU6RAFlERdHMIkpuLqHk5hJOHXwqAzL7Yy9YDN9MQNf/AbPaRZI3k5TDQ/D5Gmsi+3zma8c2magnhBDd2YWXDUT7aBhJLhhRSuGEXZRXeKxtWAwomllEVkpW44Z1QDJwLdCzcbOssNd1cR0gK6XSlVJrlFIXWN0W0X7+YNnz/q8YctY7JNuTSUtKa3gdV3EZtY7D1CW50A4vJHsbJuw5HLaEfswmRKyS/liEy8sLd6G1bniCaLdrDrlSeeOtbL7eXmlt42KA63ZX8IY6zJHka4AB5qbyunLJRe4iSwJkpdSzSqkDSqmNTbZPVkptUUptV0rd3o5T3Qa8HJlWikhTClxLZ6PLjgra3qvmaJS2k2SkY/M58HkVun7dkNpagx3bqixorRDdk/THItbs3VOLr0mWQFqqh0839WPuAzutaVSMCRpFBjNAfhO4BLCbm2TxkK6xagR5ITA5cINSyg48DkwBxgBXKKXGKKXGKqXeaPLqq5Q6B9gEHIh240X4GJV98H1xUdC294bcyT9HzGDp+CuY//wJcMlaVO9K0jPs5PRIJn9EujWNFaJ7Woj0xyKGzJk/lgG5TmwBEUqa08vwIYdZ9VGPhK2JHKjZKDLAFuCvgEHDSPLHez6OYqu6F0sCZK11EWbZ60ATge1a6x1aazfwEjBNa71Ba31Bk9cB4AzgZOBKYJZSKq7TRRJV6til+DZPC9pW663FbbjZU76H4Q8XUL6kHxVHaqiorGPPkX28uORjhj9SIL8ZCxEG0h+LWNRvoJPklOBvo3GjD3DgUBrPPVFsTaNijAqxoBYeIBO4HBhvLkMtPys7J5Y6sVwgsETB7vptIWmt/09rfTPwIvCU1v6H8MGUUjfU58WtOXjwYDjbK8LAOXw1umIAusacUFDrrcVX/6U0tEHPPadgMxw4jSx8Ni8pRiYOnIw+fBElrsQtHC9EhEl/LCw1a3Y+PXolB40iZ2fWMfXMryjeVi2jyJil30KqwHwudCowCVbtWkXe3Dzy5uZJsNwBsRQgd4rWeqHW+o1W3n9Saz1Baz2hT58+0WyaaIehPXNJ/sUoahx7qfZUNwTHfj1rjibFl4nSdsxBKRsObyrplXnWNFgI0SLpj0U4hRpFrq2zs/jdo3n2T19z9y82JnSg3KyiRaDDwLOAA3zax07XTna6dkpecgfEUoC8Bxgc8Pmg+m1dppSaqpR60uUKkbMjLFU0s4ivr68m++OHGD9gfMNKezZlw6ZsFPWbQ7X9EIeSv+JQj6840vdrPMkVFB/9gtVNF6I7k/5YWG7W7Hyyc5KCtjlTDGw2zX/fS+Xgfjf/eDqxnyS6bnc1LiDSVBXwP6AHcB4NEZ/kJbdPLAXInwAjlFLDlFLJmBk0r4XjxFrr17XWN2RnS13AWJSTA+X//RGf7ijBbbjRaHzah0/7GH1kOjZtRysDVWEjtbI3HsPAt+1o9pTvkUdGQkSG9MciZh0/ej+fbe6H1rDty6qEL/02MXdi6ztUYOYlXwEkSV5ye1lV5u0fwIfASKXUbqXUdVprL3ATsAzYDLystf7CivaJ6EpLA+fRq9Bbzmv2Xu+6Udiwk2rkkOkZQEpVJunuPvSuHoOhDcqqyyQXWYgukP5YxDKbXWGzgd3eOCFt8IByBvappM5j1jNb8Kdii1oXG4pmFjFpyCTsyk6yPZlke3LwDl7MAowVwKXmJhlFbpvDiotqra9oYftbwFvhvp5Saiowdfjw4eE+tQiTtLFLqdl4Nox7Lmj7OwNvA8Bp5DB95zMk4yFV53C4507O2ncXn/R7HNKsaLEQ3YP0xyKW9RvopK7WR021QV2tgVJgGFB44i48XjNoLit1s2m9izHHJu5TiaKZRRQuKGwYMNrpalIvWmM+A8oxP61LqYtm8+JSLKVYRIw80ot9aSe8Dhdd0+L7Y45Mx66TcOpM0DB27yUMrTydUUcuilobhRBdJ/2x6IhZs/O5d14BIwsyyemZjKpff1prWPTOSA6UmSMkz/xxh5XNjAn+VWpLbi5pOS/5CDAUmAUTFk2IXuPiUEIEyCL2De2ZC9snQ8mkZu85jRxGll+AzZeM0jYUNhxeJz67wQjXeSR5My1osRBCiGjxB8oOhxkgKwWj8sv4bHM/AGprdMLnIgdqNS+5GHgP1n5nLce/fnyUWhR/JEAWMaFoZhHKNRT7pzc2K37eMHrsy0Jhw3xWBCmeLGyGg5xdJ5E3N4/hDxdwxS//TGWF14I7EEIIEWmjjs0iNc2OzQZjhpexZ38GrgpztFRGkRv585Jb9DnwOmzosyFqbYo3CREgS1mh+JB+7LsYW89Be4NT43vXjSLFl4kNB6DqX2DDho0kjnadz+Ej1fTccwpp+0dx4b2/khm6QsQo6Y9FV8yanc/IgkyycpJITvJx+sRdKGUOmhwu8yZ0XeSmimYWtb7DVjCWGZzw6gn8i39Fp1FxJCECZMl5iw8Vc7ZyyrieLD/LjV3ZG+oh/yf31+xJ/QyND//osUmR5EjFrhyMPXwlR7vOo9pRxtBD5/BNWdOVc4UQsUD6Y9FVs2bnc9KknigbjMg7gt2mqa0zK1o8/2SxtY2LMS3mIgfYuH8jP+fn/Jk/R6FF8SMhAmQRP5YuhTPPbL79gHMDblUFTdIvbJ4kVJKdvPJJKG3HsLlR2saQA+dEp8FCCCGibv1aF2jw+WDNxv6s39IXgOJt1ZJmF2Bi7kTsyt7qPr79Plawgsd4jEd4JEoti30SIIuYkp4ODz8M2hf8rdm3dizJOiNgizmSrLTitV/9BLBRp6oA8NprGVx6pnSSQgjRTQUuQ33cqANs3NYbj8eG1rD87f0Wty52FM0s4tTBp7a6j6ENfrDgB6xkJWdzNkb9n0SXEAGy5LzFD7sdnn8e7HtOC1pR74BzAx5bNaDN7Rh4qUVpG6fP+zV2XxIZ3r4obQflQ2Fn1fJSq29HCNGE9MciHGbNzqdHr2SUgpysOgb2rWTzjl5A/eiyaFA0s6jNVIuP93xMP/oxjnH8iT9xJVdSR2LXSk6IAFly3uLL9Ong3H4FGcnmiLG/zJvdlwz1NS5s2LFjft7nyGiSySTFl0WmZwBJnkyUtrOzuNrK2xBChCD9sQiXXn1SUDaw2eGU4/eQP/gIABlZlqyBFtPaWo66zqgjb24eeXPz+MOf/sDr21+n/57+VFARpRbGnoQIkEV8mT4dPJvOo2dqT+zKzpgj07HV5xf78ODDi0YDCq+tmn1p66jJPMThHl/jSXLx3pif81XhHGbNzrf6VoQQQkRI/oh0evZKJiPDQXaGG6/Xxt6D6Rw64La6aTGnPaPIfdL60CetD31T+jL6f6PxHPLwF/4SpRbGHgmQRcw57jj4+vNBlNxcQm5WLn3dY7CrJLz2amqSDlHjOITGh1dVA4pDjmJsVXaMLDc+rRi/+Te8dclyq29DCCFEBO0srsbr1VSUm/NNXBUpFH0ymNKDbin3FoJ/FFk1+eO34cCGhtfG/RvxLPFwC7fwKZ+yk50tnbbbSogAWXLe4otSUFEBy5aZny8f9H+8cvSlvHL0pbyQP5WN2a9Q6diH21ZNki+NY1wX46YaasBmT6Zn7QjO+b+bpRayEDFI+mMRLoHl3mx2GDKwHK0Vu/ZlSrm3EIpmFmFXdlKTUoNefh7DE/RyG27OXHAmn/AJ3+JbbGKTha2PvoQIkCXnLf7s3w+33tp8e4o3i5HlF+C11eH0ZaPRJOFEaRs2lx1nbRYaH6MPT5dayELEIOmPRTg1lHszzMGVcaP3s3FrHyn31gKHzYHX5w16+Smlgl4AJa4SfsgPmcMczuRMVrPaqqZHXUIEyCL+nHIK7NsHfetOQaGo89ZR561jdH0+coqRBdr/eEjRwzOMTPdAs4qFTePQKVILWQghurl+A52kOM2lp+12xYihh/j2KcVoDfPu3ypBchMTcyfSP6N/0MvPXzXK/wLYU76HwgWFXMVVLGQhaaRZ1fSok6meIibZ7XDRRTDC/hKf3Na4fcoPH8ZOEqlGDwB0fa1GGzYMvJSnfIXuZ2ArczTUQs7IlG9zIYTojmbNzuepeTvYuaOammoDrQ3q3A4+25yDUgdYtbyUc6f1b/tECSLU8tPO+5zUGaFLutltdkpcJQBMYQoAt3ALJ3AC3+N7kWtoDJARZBGz7rkHfvSjxs8LFxSydOCtDTnIZSnbOJJsThzQaJw6y1w/pAp8WV4MQ3PGb3+I8z5nQ/kayUsWQojuZdbsfG67bxSZ2Q4GDEplSF4Sn23qT0VVCste2ysT9trQL6MfyfbkoAl7fm7D3TCK7Hct13I7t/NH/hjNZkadBMgiZvXvDytXwoED5uf+32J7143Chp1UI4dszyAUtoZUi0zPAFJdOaRW5WBTdvrVjUGjG8rX+M8hhBCi+1i1vBSfobHbFR63m4IRpXy2qQ91tZqXF+6yunkxLS87D8Nn1JdPDaZQQaPIAMdwDCtZyeM8zlu8Fc2mRlVCPHtWSk0Fpg4fPtzqpogOuvqed/EO+h+Zp/2dPeV7AHhnoJlzke0ewuXF/8SLG6XgiGMndpJYNORa6uwuVLLCqZ3gs/IOhBCBpD8WkbCzuBqfD8pdHirKPYwddYAl/xmBzwfbNlewf28t/QY4rW5mTCqaWUThgkI+2PUBhg5eYlqjg0aR/SkaeeTxER+RRRZb2MJwhmPHbkXzIyYhRpBl1nT80iMXYWyaSp+0PthtdhQKmzK/bSftvxVQoHygwenLwqbtjKv4LgDqKoVvqETHQsQS6Y9FJMyanc+98wooPLsPPXsl07+f4orzN2Ozgc+HjCK3oWhmEblZuQ2pFjZla3iFGkUGyCEHGzZ+y2+5jMuopXulsiTECLKIX87R73Honw/jrcoEzDI0WmucRg79a48FNDZt/taaavSk2l5KdvVRkAO+1T7qTq+DbWYBdDAfFxUuKAyZapGXnRdyAoMQQoj40LB4iMuD16t4uyifcyd9zY5tlTJpuw152XnsKd+DRqN1cLqF23Czv3J/yONe4AWu4irO4zxe5VWyyIpGcyMuIUaQRfyypVQz8lfXY3dWA+B0OElNSmVs+SX4Uqo45PyKspTtlKVsp9Kxjy9y/t2QgsEmQAGjzdqP/vqPJa6ShpzkwJfkJwshRHybNTufwrP74Ey1k5Sksds1W77uhcNhY9XyUqubF9P8o8hA0AiyfxQ5sGZyoBRSeImXOJET2ce+aDY5ouRXKRHzUnrvoWLreJL7bsJtuAHoVXM0Ph+kGj1I8WWhNdQ5yuldNyr44DcAD1R7qhs27XTtZJdrF0n2JJJsSYzq3eQYIYQQcWvHtipqawx8Pjh+zD7eXTWMghEH2bGtyuqmxby87Dx2unY2G0FWSmFog7y5eSGPKZpZxO/5PRrNr/gVP+EnDGNYtJodETKCLGKe9ibz1V8fYkTmsYztO5axfceyafR8FiyYzJUXn4Cq/+12U86/+U/urxtylAHYC6meVMil4bdgfykbh82Bx+ex6K6EEEJEQv6IdJyp5uIh/XtXM3hABR6fk/wR6VY3Lea1tBy102FOcGzr6atCkUcek5jEBjZYdRthISPIIqblZedR4tpO0uDP2fXZGNLGvtOwvbLCy4fvl6J8NhQwruxqdqZ+RIHrYj7oM486uwuAmswauAR8f/JBwARd/6jy2r1rG69X/9ux5CMLIUR8ChxBBjhjolkv/6utMoLcHv50xFC2lG1peJLrZ/iMoAoXN3ETvenN2ZzNalbH7UhyQgTIUlYofvn/w813wMcfn8ZzNze+t2zJPqoqDbOQhbaRpNM5e+/vcPpycCXt4tNeC8wddwP7gBOhrWXk+6T1AZB8ZCEiRPpjEWn5I9LZsbWyIUi22eCt94czaKQsO90e/TL6NfwsDLRu3zrchhuHrXno2PRn5uVcTgEFDGUoRzhCDjmRam7EJESKhZQVin9XXAG3BSw57R89rqszUIr6xUIgxzuEWnWEkeUXkGIEfL2XAydDiIWCGoRaRUgIEV7SH4tIazqC7PPBsEGHeO4FqYPcHnnZeRysPtjsFSowbk0BBXjwMJGJPMdzEWpt5CTECLKIf716QXk57NgB+fnmqklVlQYKRZJKwRuQO5Hm60WtOsIxR2Y0jiIfAJ6AEAsFNdBoPt37acNqQs77nPTL6AdIyoUQQsQL/wiyu86HsoHXoxk7ysXnrw3ms8/g+OOtbmFsa+lnXd7cPPZVmlUqar21DRP5NJo95XvIm5vX7GdlMsksYQmTmUwppfyCX0T+BsIkIUaQRffwz3/CI4+YH+8srsbj8WH4NIahoX6paVAk6TQ8qrb5KHINMAVIDX1+/wiyfzKfLFEthBDxZ8e2KmprDbQGr1ejNXjcBtPPO2h10+Kaf0lqr8+LTzcuwmVTNuw2e4s/K0czmpWs5A3e4BCHotnkLpEAWcSNGTPg1VfNx2WzZuczd8HxTL1kIL37JpPTIxkDD15Viw9vw6p6xxyZEXwSBUyyovVCCCGiIX9EurmaXq4Tu12RlKRQNrjgAoOjj4ZD8ROjxRR/neSxfceSbE9uVuGiNYMZzHKWk0UWj/M4XmI/H1xSLETcOPpo6NEDPv4YTj7Z3LazuBqfDw65qrBhA21DYSPFyKLWcbh5XeQi4MeYk/XKg9/yp1b4Hxu5DTdfln4pdZKFECKONPxcKHXjMzRKKdCwfq2LNRtzObDfx2kn7GT6lYNkZb0OMitLlWD4jKDtSbYkAPZX7m+1VrIbN6/xGu/yLv/gH6S29Eg3Bsh3hogrzzwDRx3V+Pms2flAY27U8WXXMKb0u6zv+XfW9ny2+QkqMYPj/kA5DakUAGlJadR4aszOFDNQljrJQggRX2bNzqeywsuj92yhZ+9k7HaFYWg8Hh9Xf9/L8SfYcNRW0rd/KedO6291c+OKP784b25eyEoXXp835HZ/6kUaabzO61zN1UxmMstYhpPYnDwpKRYirpx0EpSUgA4x2S7Fm83wI5OptpdytOv84PzjQCuArch3vxBCdFOrlpfiMzR2uzngYbcrfIZm02f7OXqYi81f9WX1ijIqK2L/UX93k0wyL/IiP+NnpJCCD1/bB1lARpBF3LnkEnjjDSgoaNyW5MnkvJ0PkeruhSt5D0m+VKbuns/rg26izu7CaeRwysGfNS4gMgYYDfrfjZG2f+GQwCU23Ya7YSERx+/M/y4Om4OJuROlqoUQQsQof5pF0wB4/VoX35pQiUESrsMelr+9nwsvy7WolfHLn2rRVKhScFvKtlDjqWmWepGXncfJM09mClNYzGKO4qhmx1pJAmQRV5SCiy6CxYuDA+Rh+y6gZ91RaDRpRg+0VvT2juS4Q1eypu8TjDkynaGVpzcuILINmIyZarGv8Tw2ZUNrTWpSKjWeGjQahSI1qTFPyuvzSlULIYSIYf70u0D+tIt+/XwcOFDLzj1OVi4v5awp/SQXuYNaKwXXlNtwN1S5CFTiKiGXXH7MjymkkDd5k3GMi0RzOyUhHjIrpaYqpZ50uVxWN0WEwfTpsGhR4+eVFV6GHpyMv8xbipGD05eFxscY1wyy6wYzsvwCqu2ljaXfPJipFt+25h6ESFTSHwur+NMuAL7Z6+DdVcOoqjRY/vZ+i1uW2H7Mj5nLXK7juphKt0iIAFlWbupevvUtuP32xjzkVctLcfjSADNEtuHAhgPQOLST0/b/Epu2Y9jcwaXf1gJbgs+ttUaj8fq8DZP3hBDhI/2xsEpgdYt+vSpJT/XwVUkP1q+VX9bCJdQqfIbPaKhy0ZJLuZQP+RAvXlawIkqtbZ08UxBxx2430yw2feHl4/d28fX2SpLsDmw2RY27FjtJgMKVtButNLk1EzicVIJdJ5Pu7cto10V8kbOIOlywBhgM7DLP7Q+K3Ya74XoaTa23tl21HoUQQsSmptUtzjnjCP9b1Z9TTtpBZYVX0izCIFTqhfM+Jx6fhw0HNgRt9y/O5ZdMMlvZymVcxhzmMJOZEW1rWxJiBFl0P//5D1z+XR+frz2C64gXVf/Hpuz4Uy0yPbmk+LIAhdOXRYanP3aSSTN6BS8gcj5wdOvX06HKZgghhIgrgdUtRo2o5Hsz9uAzNKuWl1rdtG7L6/PisDmavby+5hVEjuZo3ud97uEeHuIhC1rbSAJkEZdOmuhl21d2auqS8HrM5aa9Xh9K2xv2sZNMsjcDhSLFyCZZp+NPuxhQc0LjyZYDZ7V+PY2mzltHnbcOhSIvu/lEBCGEELEtsLpFdaUXj9vDh2t7sbO42uqmdVv+YLjpK1TFCzCD5FWssrw+sjxPEHFp7QelDM9zUurqRcHRZZx6Rm/OndafPz2wja2bKnA4bHi9PkaP6csjqWdz3Be3MbB6Aj7lwaaTUPUDwgqF82snNe4aGA1sNrfblPm7o0/7GD9gPAerD1Jys1SuEEKIeNa0ukVdHeTnw0OnW9SgBNAvo1/IxUMOVh9s8ZhccvkZP+MjPuJpnubP/JkkWs9jDjcZQRZxp7LCy+oVZRSecpiB/epIcdpZvaKMr7dXsm1zJbb6QWSbHbZtrqTnkWMZWD0eXb/2u8ZLbs0EstyDgfo8qJeBLy26ISGEEJZISYHZs+GRR6xuiQhlLGPZwx5mMINqojvKrxIpt3LChAl6zZo1VjdDdNGyJfv44L1S0tIdVNfYSHJoPG5zSehDpW4cDhsajdvjRWuNDy8KOz7lBRR27UADe1I/5o3BP208cR6QAXzR8TbZVWNqR3pyOq7bZVa0aJlSaq3WeoLV7bCS9MciVrhccNxxsGEDZGZa3ZruJ29uHodqDgVNfgcwfAanDj61zUW3PHgaJuy9wAthb19L/bGkWIi4E5hD9s/Xh3DM0S4KRrood5lBstfrM0vAaRsKsGNH48OukwBb/cxZTU93k1V73MC5mMtQe0Jf25964S8HB+YIdIojpWGfKndVGO9WCCFEJGVnw7ZtkBTdJ/gJIy87jz3le7Db7EHbnQ5nuxbdSiKJ53iOAxygllqOcIT+9I9UcxtIgCziTmAO2dBjYenSLO6dN7hhm7+Mz1eHviLd3Z8a5wFs2NhSMJ+8ddfjw4sNB4uGXBt84r3ATuAkYGVUbkUIIUQMcDjg6qvhscegVy+rW9O9FM0sIm9uXofzkAPZsNGf/ixiEb/klyxjGSMYEe6mNrmmEHHswgvhnXegpqZxm7+Mj8OXChochhOlbeRtuzz0giGB/gcMiVrzhRBCxAClzCD5L3+xuiWiNTOYwW/4DadzOp/yaUSvJQGyiGt9+piTK+rqzM/9E/gcSTaSvZloZeDwZGAog2zXKDyqFgC3rbpx2elAZcCL0b0HIYQQ1vvVr2D+/OABFxF7rud6HudxlrM8oteRFAsR92bNgsOHzY9XLS/F6/FRUe4FrbDhQCsfzpregI1szyAM5QGlcdsqOebIDD7ttaDhXE4jh1MO/owPfjWPusUuqAy+lk83Xydeo6n2mLNrFQqNRt3TuEKQfwJfuq8X51XexVN33SArNgkhRIwZPdp8Kvnll3D88Va3pvvbUraFGk8NeXOD1xXIy85rc+LedKYDsIxlVFHFDEI8Ee4i+Skt4t7evXDssebfO4urqaky8Lg1tvpvb6XtDQtaOnDi0CmgzWC2d92ooHONOTKdoZWn43pzF58WLoC3OtaWwIl7Sim01g0T+EYeuJCMQ2NYtbyUc6dFfoKBEEKIjnnySfNvrc20CxEeedl5zSbk1XhqcDqczXKT2zNxz68PfZjJTMooYxazwtJWP0mxEHFvwADIy4MVK+CKa4eQmZNEn37JKBR2u7kEtQ8PoOtfoDHwJdfwfr85pCWlkZaURqqRw8jyC6i2lzLyswtIycuGnPa3I9mezPgB40m2Jzd7L8WbzQjXedQmHWb1ijIqK5ovsSmEEMJ6N94Ir7xidSu6l6KZRZTcXBL0ys3KZVTvUW0f3IoTOIH3eZ8HeIAneCJMrTXFbYCslDpDKbVCKfVXpdQZVrdHWGv6dFi0qHGCXoWrflGQ+owIhQNQ9S8AGw5PetBEvdFHpjdO4jPsHPPUDAiuStNpRx+eik3b8dk8+AzNquWl4TmxEDFA+mPRnZxzDjz0kDmKLGLfCEawkpVMZnJYz2tJgKyUelYpdUAptbHJ9slKqS1Kqe1KqdvbOI3GzBB1Arsj1VYRHy6/HMaMMVi9ogy7XVFba0bGOmBJ6UbmqLLNSGGU60JSvFmkeLMZWX4BbpuZS+y2VTNy2wWkHMmGLhaOT/FmM/zI5IZz+1f+k1FkEQukPxYi2IUXQkUFvP++1S0R7TWQgeSR1/aOHWDVCPJCCA71lVJ24HFgCjAGuEIpNUYpNVYp9UaTV19ghdZ6CnAbcE+U2y9izFFHwbABB6mrA9fhFlb5aEJhI8XI4ujDFzaM8GplAKCVYZaCS5kB53WtbUcfnooKOLfdrmQUWcSShUh/LEQDmw0eeADc7rb3Fd2XZUtNK6WGAm9orQvqPz8FuFtrfW79578G0Fo/0MZ5koEXtdaXtPD+DcAN9Z+OBLYEvJ0NuNr4uDcQjkgm8Jxd2bel90Jtb7otFu+3rf3ac18A2dmZg3s57M50m83RPAm4BRofbptZqsKhU5u971U1lCft6cjpPNhIojGXgyxPrnnu+vRn5VMeAK9RW+Wq2LUj8B6Qr22o7R35PNz3G6n/t3la6+ZV8y0i/XGn942H7+HW2tnR/drdH9Px+43Xr22obbH4vRx7/XFvRmKn+c9rAzelQX1DW20Mpev9sdbakhcwFNgY8PklwNMBn38fmN/K8TOAJ4B/Amd0sg1PtuPjNWG63yfDsW9L74Xa3nRbLN5vW/u1577auteW7le+trH3te3o1zLS9xvpr22svKQ/Du/XNZa+hztyv9Ifd/xr28GvZ1x9bTv6tYz0/Ua7P47bMm9a60XAoi6e5vV2fBwuHTlna/u29F6o7U23xeL9trVfe+4r1Lb23ns4yNe2c/u1936789e2W5D+uM3tVn4Pd+Sc0h+3/Z70x+3/PN6+tkHiPsUi0pRSa7TWE6xsQzQl0v0m0r2C3K9omfTHsSmR7jeR7hXkfuNBLJV5+wQYoZQaVp/HdjnwmsVtAnjS6gZEWSLdbyLdK8j9ivaT/jg2JNL9JtK9gtxvzLNkBFkp9Q/gDMyk7f3AXVrrZ5RS5wFzMavPPqu1vj/qjRNCiAQi/bEQQjRnWYqFEEIIIYQQsSiWUiyEEEIIIYSwnATIQgghhBBCBJAAWQghhBBCiAASIHeQUipdKfU3pdRTSqnvWd2eSFJK5SulnlFKvWJ1W6JBKXVR/df1n0qp71jdnkhTSo1WSv1VKfWKUurHVrcn0ur/765RSl1gdVtEeEh/3H1Jf9y9xUN/LAEyoJR6Vil1QCm1scn2yUqpLUqp7Uqp2+s3zwBe0VrPAi6MemO7qCP3qrXeobW+zpqWhkcH7/fV+q/rj4DvWtHerurg/W7WWv8IuAw4zYr2dkUH/98C3Aa8HN1Wio6S/lj6Y5D+2Ir2dkV37I8lQDYtBCYHblBK2YHHgSnAGOAKpdQYYBCwq343I4ptDJeFtP9eu4OFdPx+f1v/fjxaSAfuVyl1IfAm8FZ0mxkWC2nnvSqlzgE2AQei3UjRYQuR/lj640bSH8eHhXSz/lgCZEBrXQQcarJ5IrC9/rd2N/ASMA3YjdkpQxz++3XwXuNeR+5XmX4PvK21/jTabQ2Hjn59tdavaa2nAHH3eLqD93oGcDJwJTBLKRV3/3cThfTH0h8j/bH0xzHAYXUDYlgujSMTYHbEJwF/BOYrpc4nMuuqWyHkvSqlegH3A8crpX5t9TKzYdTS1/anwNlAtlJquNb6r1Y0LgJa+vqegfmIOoX4HLEIJeS9aq1vAlBKXQOUaq19FrRNdJ70x9IfS38cf+K6P5YAuYO01lXATKvbEQ1a6zLM/K+EoLX+I+YP3ISgtX4PeM/iZkSV1nqh1W0Q4SP9cfcl/XH3F+v9cUwOa8eIPcDggM8H1W/rjhLpXkHutzvfbyLdayJJpK9rIt0ryP125/uN63uVALllnwAjlFLDlFLJwOXAaxa3KVIS6V5B7rc7328i3WsiSaSvayLdK8j9duf7jet7lQAZUEr9A/gQGKmU2q2Uuk5r7QVuApYBm4GXtdZfWNnOcEikewW53+58v4l0r4kkkb6uiXSvIPfbne+3O96r0lpb3QYhhBBCCCFihowgCyGEEEIIEUACZCGEEEIIIQJIgCyEEEIIIUQACZCFEEIIIYQIIAGyEEIIIYQQASRAFkIIIYQQIoAEyEIIIYQQQgSQAFmIdlBKXa6UqlNKJVndFiGESGTSH4tokABZiPY5DtiktfZY3RAhhEhw0h+LiJMAWYj2OQ74zOpGCCGEkP5YRJ4EyEK0zzhgnf8TpdRZSqkypdRcpZTdslYJIUTiGYf0xyLCJEAWog1KqT7AAOo7ZKXUT4E3gNu11jdrrQ0LmyeEEAlD+mMRLQ6rGyBEHDgO0MAmpdRTwAzgPK31e5a2SgghEo/0xyIqJEAWom3jgEPAIqAPcJLWerulLRJCiMQ0DumPRRRIgCxE244DFDAJOFs6YyGEsIz0xyIqlNba6jYIEdOUUuuBfwPDgdOB8Vrrg9a2SgghEo/0xyJaZJKeEK1QSqUAo4D1wA1AKfAvpZQ8fRFCiCiS/lhEkwTIQrRuDJAEbNBa1wDTgWOAxyxtlRBCJB7pj0XUSIAsROuOA6qArwC01iXAd4EfK6WusbBdQgiRaKQ/FlEjOchCCCGEEEIEkBFkIYQQQgghAkiALIQQQgghRAAJkIUQQgghhAggAbIQQgghhBABJEAWQgghhBAigATIQgghhBBCBJAAWQghhBBCiAASIAshhBBCCBHg/wEV2fLTr6b5mwAAAABJRU5ErkJggg==\n", 106 "text/plain": [ 107 "<Figure size 720x360 with 2 Axes>" 108 ] 109 }, 110 "metadata": { 111 "needs_background": "light" 112 }, 113 "output_type": "display_data" 114 } 115 ], 116 "source": [ 117 "def linear(x, a, ε):\n", 118 " return -a*x + ε\n", 119 "\n", 120 "def curve_fit_log(xdata, ydata):\n", 121 " xdata_log = np.log(xdata)\n", 122 " ydata_log = np.log(ydata)\n", 123 " popt_log, pcov_log = curve_fit(linear, xdata_log, ydata_log)\n", 124 " ydata_fit = np.exp(linear(xdata_log, *popt_log))\n", 125 " return (popt_log, pcov_log, ydata_fit)\n", 126 "\n", 127 "def get_dist(G, directed=False):\n", 128 " degree_sequence = sorted(list(dict(G.degree()).values()), reverse=True)\n", 129 " degreeCount = collections.Counter(degree_sequence)\n", 130 " k, p_k = zip(*degreeCount.items())\n", 131 " k = np.array(k)\n", 132 " p_k = np.array(p_k)/sum(p_k)\n", 133 " CDD = np.cumsum(p_k)\n", 134 " if directed==True:\n", 135 " degree_sequence_in = sorted(list(dict(DG.in_degree()).values()), reverse=True)\n", 136 " degreeCount_in = collections.Counter(degree_sequence_in)\n", 137 " k_in, p_k_in = zip(*degreeCount_in.items())\n", 138 " k_in = np.array(k_in)\n", 139 " p_k_in = np.array(p_k_in)/sum(p_k_in)\n", 140 " CDD_in = np.cumsum(p_k_in)\n", 141 " \n", 142 " degree_sequence_out = sorted(list(dict(DG.out_degree()).values()), reverse=True)\n", 143 " degreeCount_out = collections.Counter(degree_sequence_out)\n", 144 " k_out, p_k_out = zip(*degreeCount_out.items())\n", 145 " k_out = np.array(k_out)\n", 146 " p_k_out = np.array(p_k_out)/sum(p_k_out)\n", 147 " CDD_out = np.cumsum(p_k_out)\n", 148 " return (k, p_k, CDD), (k_in, p_k_in, CDD_in), (k_out, p_k_out, CDD_out)\n", 149 " return (k, p_k, CDD)\n", 150 " \n", 151 "(k, p_k, CDD), (k_in, p_k_in, CDD_in), (k_out, p_k_out, CDD_out) = get_dist(DG, directed=True)\n", 152 "\n", 153 "popt, pcov, CDD_fit = curve_fit_log(k[:-1], CDD[:-1])\n", 154 "popt_in, pcov_in, CDD_fit_in = curve_fit_log(k_in[:-1], CDD_in[:-1])\n", 155 "popt_out, pcov_out, CDD_fit_out = curve_fit_log(k_out[:-1], CDD_out[:-1])\n", 156 "\n", 157 "\n", 158 "def plot_pk_CDD(k, p_k, CDD, CDD_fit, popt,\\\n", 159 " direction=\"u\", node_color='slategrey',\\\n", 160 " alpha=0.8, fit_color='r', marker='.',\\\n", 161 " p_n=[0, 1]):\n", 162 "\n", 163 " ax[0].loglog(k, p_k, marker, alpha=alpha, c=node_color , label=direction)\n", 164 " ax[0].plot(k, np.exp(linear(np.log(k), popt[0]+1, popt[1])),\\\n", 165 " c=fit_color, lw=1, linestyle='--', label=fr\"$γ_{direction[0]} = {round(popt[0], 2)+1}$\")\n", 166 " \n", 167 " ax[0].set_ylabel(fr'$p_k$', fontsize=15)\n", 168 " ax[0].set_xlabel(fr'$k$', fontsize=15)\n", 169 " ax[0].legend(loc='best')\n", 170 " ax[0].set_ylim((1e-5, 1))\n", 171 " \n", 172 "\n", 173 " ax[1].loglog(k, CDD, marker, alpha=alpha, c=node_color, label=direction)\n", 174 " ax[1].loglog(k[:-1], CDD_fit, linestyle='--',c=fit_color, lw=1,\\\n", 175 " label=fr\"$γ_{direction[0]}-1 = {round(popt[0], 2)}$\")\n", 176 " ax[1].legend(loc='best')\n", 177 " ax[1].set_ylabel(r'$P_k$', fontsize=15)\n", 178 " ax[1].set_xlabel(r'$k$', fontsize=15)\n", 179 " ax[1].set_ylim((1e-5, 1))\n", 180 "\n", 181 "fig, ax= plt.subplots(1, 2, figsize=[10, 5])\n", 182 "plot_pk_CDD(k, p_k, CDD, CDD_fit, popt, marker='.', direction=\"und\")\n", 183 "fig.tight_layout()\n", 184 "fig.savefig('../pres/pics/dist_u.png', dpi=300)\n", 185 "\n", 186 "fig, ax= plt.subplots(1, 2, figsize=[10, 5])\n", 187 "plot_pk_CDD(k_in, p_k_in, CDD_in, CDD_fit_in, popt_in,\\\n", 188 " direction='in', node_color='green', fit_color='lime', marker='s')\n", 189 "plot_pk_CDD(k_out, p_k_out, CDD_out, CDD_fit_out, popt_out,\\\n", 190 " direction='out', node_color='slateblue', fit_color='blue', marker='^')\n", 191 "fig.tight_layout()\n", 192 "fig.savefig('../pres/pics/dist_d.png', dpi=300)" 193 ] 194 } 195 ], 196 "metadata": { 197 "kernelspec": { 198 "display_name": "Python 3", 199 "language": "python", 200 "name": "python3" 201 }, 202 "language_info": { 203 "codemirror_mode": { 204 "name": "ipython", 205 "version": 3 206 }, 207 "file_extension": ".py", 208 "mimetype": "text/x-python", 209 "name": "python", 210 "nbconvert_exporter": "python", 211 "pygments_lexer": "ipython3", 212 "version": "3.10.5" 213 } 214 }, 215 "nbformat": 4, 216 "nbformat_minor": 5 217 }