diff --git a/Exercise sheet 4/exercise_sheet_04.ipynb b/Exercise sheet 4/exercise_sheet_04.ipynb index 5e4508d..5037f6f 100644 --- a/Exercise sheet 4/exercise_sheet_04.ipynb +++ b/Exercise sheet 4/exercise_sheet_04.ipynb @@ -137,7 +137,7 @@ "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAb4AAAEuCAYAAADx63eqAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABIgElEQVR4nO3dd1gU59oG8HsbLIqAVBHsDRsaG0SioigqRY1IIEcTezcxOcYkxmPKMeZoYosxdo2JJSBYIoINpNhQQcWKiooBlS5SZJct8/3hxyaEIsjuzszu87surkt2h9lnFfee9523CBiGYUAIIYQYCSHbBRBCCCH6RMFHCCHEqFDwEUIIMSoUfIQQQowKBR8hhBCjQsFHCCHEqFDwEUIIMSoUfIQQQowKBR8hhBCjQsFHCCHEqFDwEUIIMSoUfIQQQowKBR8hhBCjQsFHCCHEqFDwEUIIMSoUfIQQQowKBR8hhBCjQsFHCCHEqFDwEUIIMSoUfIQQQowKBR8hhBCjQsFHCCHEqIjZLoAQwp68EjnCkzORmlWEIpkSFlIxXJpZILC3M2zMTdkujxCdEDAMw7BdBCFEv1IyCvFzXBri7+YCAORKteY5qVgIBoBnJzvMGdQePVpYsVMkITpCwccDdFVOtGl3YjqWRaVCplShtv/9AgEgFYuw2McFE9xb660+QnSNgo/D6KqcaNvL0LuNMoX61Qf/PzOJEIt9OlP4EYNBwcdRdFVOtC0loxDBWxNRplBVejwvYiVk6SlQK2QQNW4KC/cANOkxvNIxZhIRQme4w9XZSo8VE6IbFHwcRFflRBdm7ErCydvZVS6kynMfQdK0OQRiCRT5Gcjauwj2gV/DtFl7zTECATC8iwM2Teij56oJ0T4a1ckxKRmFWBaVWmPold6KR+HZ36EqyoWocVPY+H4EaYtuKFOosSwqFa7OVnRVTqrIK5Ej/m5utb0HJnat/vadAAIIoHz2tFLwMQwQeycX+SVyuq9MeI+Cj2N+jkuDTKmq9rmyh1fwLG4n7EZ/BpPmHaEqKaj0vEypwoa4NLoqJ1WEJ2fW+nz+8Q0ovR4DRimHiUM7mLWr+jskABB+ORMzB7bTUZWE6AcFH4fUdlUOAM/P7IGlx7swdXIBAIib2FZ6nq7K2cMwDNRqNVQqFZRKpebrn99X95g+jkmSdIXcxLnG+m2Gz4H1sJmQP06F7M/rEIgkVY6RKdVIfVqsy79GQvSCgo9DarsqZ9QqyJ+mway9Gx5vmg5GVY5GHdxhNXgKhJK/Qk6fV+VqtVqnH9xcCY26HiMQCCAWizVfIpGo0vfVPaaLY0xMTGBmZlbp+ZvPrAFZ7f+eAqEI0hZdUXozFsVXomDRZ1SVY4pkCh39NhGiPxR8HJKaVVRpysLfqUoLAbUSL+6chcOEFRAIRcjd/y2enwtF00Hva46TKdXYsf8YzmxN1PmHPcMwr/xA1teHvVQqrdN5dFWfSCSCUMjdFQAfhF7B/atP6nawWg3ls6fVPmUhrdoSJIRvKPg4pEimrPE5wf+36pr09ofY3Prln/uOqRJ8AGBp5wjvPt46+XD/+/dc/qAnlbk0s4CpOKvKhZWqtBCyRykwa98PArEJZOlXUXo7Hrb+C6ucQyIEOjVroq+SCdEZCj4OsZDW/M8hkppD9I97ejXp2qEt3gvqqaWqiCEY19sZa6LvVn1CIEDxlaPIP74BYNQQW9qjqdd0NOroXuVQhUKBn/49Huqp72PChAmwsLDQQ+WEaB8FH4fUdFVewbz7UBQnH4FZ296ASIzipD/QqH3fSsdIxUK4ONJVOanM1twUgzraVZnHJ2pkiWbjl7/y5wUCwNvVGUEjvsPGjRuxePFiBAUFYfbs2ejRo4cOKydE+6ivikPG9a551B0AWHoEw8SxAx5vmYknW2fBxKEdLPsHVTqGATCuV+3nIcZprmd7SMWi1/pZqViEuZ7tMWTIEISFheHmzZtwcnKCn58f+vfvj127dkEme8XoGUI4glZu4ZiaVteoC1pdg7yKtlcFUiqViIyMxMaNG5GcnIyJEydi1qxZaN++fdUTEcIR1OLjmIZelc/xpA8cUrMJ7q2x2KczzCQiCAS1HysQvFyjs7al8MRiMUaPHo1jx44hMTERQqEQ/fv3h7e3Nw4ePAilsuYBW4SwhVp8HERrdRJdu5ZZiA1xaYi9kwsBXk6DqVCx88fgTnaY49m+3kvgyWQyhIeHY+PGjXj06BGmT5+OadOmwcnJSavvgZDXRcHHUbQ7A9GH/BI5wi9nIvVpMc4nX4WlmQRjvdwxrpd29nq8du0aNm3ahJCQEAwePBizZ8/GkCFDaCoMYRUFH4fp8qqckH9atWoVMjIysHbtWq2fu7i4GHv27MHGjRtRVlaGWbNmYeLEibCxsdH6axHyKhR8PPD3q/ITcafRvVM7DO7VSWtX5YQAQGhoKMLDwxEWFqaz12AYBufPn8fGjRsRERGB0aNHY/bs2XBzc4PgVTcdCdES6m/gARtzU8wc2A5rgnqi3eOT+FdbJWYObEehR7TKyckJjx8/1ulrCAQCzfSHtLQ0dOvWDePHj0evXr2wZcsWlJSU6PT1CQEo+HhHKpXSfCmiE05OTsjMrH37Im2ytbXFwoULce/ePSxfvhxRUVFo2bIl5s2bh5s3b+qtDmJ8KPh4hoKP6Erz5s2RlZUFtbruo4m1QSgUYvjw4Th06BBSUlJgY2MDb29vDBw4EL///jvkcrle6yGGj4KPZyj4iK6YmprCysoKOTk5rNXQokULfPPNN0hPT8f8+fOxfft2tGzZEosWLcLDhw9Zq4sYFgo+nqHgI7qk7+7OmkgkEgQEBCA6OhoJCQmQy+Xo27cvfHx8EBERAZVKxXaJhMco+HiGgo/okrOzs84HuNRXp06dsHr1amRkZOCdd97Bt99+i7Zt22LZsmXIyspiuzzCQxR8PEPBR3RJHyM7X5eZmRkmTZqECxcu4ODBg0hPT0fnzp0RFBSEuLg40MwsUlcUfDxDwUd0iStdna/Sq1cvbN26Fenp6RgwYADmzp2Lrl27Yt26dSgsLGS7PMJxFHw8Y2pqSsFHdIaLXZ21sbS0xLx583Djxg1s2rQJ58+fR5s2bTBt2jQkJyezXR7hKAo+nqEWH9ElLnd11kYgEGimP6SmpqJdu3YICAhAv379sGPHDrx48YLtEgmHUPDxjFQqpXlNRGf40tVZGwcHByxatAj379/HV199hQMHDqBly5b46KOPkJqaynZ5hAMo+HiGWnxEl/jW1VkbkUgEX19fHDlyBElJSWjcuDE8PT01u8grFAq2SyQsoeDjGQo+oksWFhZgGAZFRUVsl6JVrVu3xrJly/Dnn39i5syZ+Pnnn9GqVSssWbIEGRkZbJdH9IyCj2co+IguCQQC3t7nqwsTExPN9Ifo6Gg8f/4cPXv21Owir+/l2gg7KPh4hoKP6JqzszPv7/PVRZcuXbBu3Tr8+eef8Pf3x+LFi9GhQwd8//33yM3NZbs8okMUfDxDwUd0zZBbfNVp3Lgxpk2bhqSkJPz++++4ffs2OnTogPHjx+PMmTM0Md4AUfDxDAUf0TVjC74KAoEA/fr1wy+//IIHDx6gT58+mDp1Knr06IENGzYY3H1PY0bBxzMUfETXDGFKQ0NZW1vj448/RmpqKtauXYvY2Fi0bt0as2bNQkpKCtvlkQai4OMZCj6ia4Y0paGhBAKBZvrDzZs34eTkBD8/P80u8vR/kZ8o+HiGgo/omrF2db6Ko6MjlixZgocPH+Kzzz7Dnj170LJlSyxcuBBpaWlsl0fqgYKPZyj4iK5RV2ftxGKxZvrD+fPnIRAI0L9/f3h7e+PgwYNQKpVsl0heQcDQkCVeKS0thb29PUpLS9kuhRgolUoFMzMzlJSUwMTEhO1yeEEmk2H//v3YsGEDHj16hOnTp2P69Olo3rw526WRalCLj2cqdmeg6xWiKyKRCA4ODnj69CnbpfCGVCrF+PHjcfbsWURFRSE7OxvdunXT7CJPE+O5hYKPZ8RiMYRCIXWnEJ2i7s7X5+rqqmn5DRs2DAsWLICLiwtWr16NgoICtssjoODjJbrPR3SNRnY2XJMmTTBr1ixcvXoVO3fuxJUrV9CuXTvNLvLUa8MeCj4eouAjukYjO7WnYvDLrl27cO/ePXTt2hXjx49H7969sXXrVrpfzwIKPh6i4CO6Rl2dumFra4uFCxfi7t27WL58OaKiotCyZUvMmzcPN2/eZLs8o0HBx0MUfETXqKtTt4RCoWb6Q0pKCmxsbODt7a3ZRZ42m9YtCj4eouAjukZdnfrj7OyMb775Bunp6Zg/fz62b9+Oli1bYtGiRUhPT2e7PIMkZrsAUn8UfETXqKtT/yQSCQICAhAQEIC7d+9i06ZN6NOnD9zc3DB79myMHDkSIpHotc6dVyJHeHImUrOKUCRTwkIqhkszCwT2doaNuamW3wn30QR2HhowYAD+97//4a233mK7FGKgysrK0LRpU5SVlUEgELBdjtEqKytDaGgoNm7ciOzsbMyYMQNTp06Fg4NDnX4+JaMQP8elIf7uy/0F5cq/5hNKxUIwADw72WHOoPbo0cJKB++Am6irk4cqJrEToitmZmZo3Lgx8vLy2C7FqJmZmWmmPxw4cADp6elwcXFBcHAw4uLiap0SsTsxHcFbE3HydjbkSnWl0AMA2f8/duJWNoK3JmJ3YrqO3w13UPDxEHV1En2g+3zc0qtXL2zZsgXp6el46623MHfuXHTt2hXr1q1DYWFhpWN3J6ZjWdRtlClUeFWfHsMAZQoVlkXdNprwo+DjIQo+og/Ozs50n4+DLC0tMW/ePNy4cQObNm3C+fPn0aZNG0ybNg3JyclIySjEsqhUlCkqt/BUZcXI2f8t/lwVgMwNk1F6M67S82UKNZZFpeJaZqH+3gxLKPh4iIKP6AO1+LhNIBBopj+kpqaiXbt2CAgIQODX21CmqLqkYcGJjRCIJHD+YDds/T9B/okNKM99VOkYmVKFDXGGv8USBR8PUfARfaDg4w8HBwcsWrQIF1JuQe3QGUDlAUnqchle3DkHq4ETIDQxg7RFVzRq74bSm7GVjmMYIPZOLvJLDHseIQUfD1HwEX2grk7+OXj1KUSiqh/ryoLHEAiFkFg7aR6T2LeB4h8tPuBlZIZfNux/dwo+HqLgI/pALT7+Sc0qqjJ6EwDUijIITBtVekxo2gjq8rIqx8qUaqQ+LdZZjVxAwcdDFHxEHyj4+KdIVv12ZUKJGRh55ZBj5C8gNDGr4TwKrdfGJRR8PETBR/SBujr5x0Ja/WJcYmsnMGoVFAV/XciU5zyExK5VDeeR6KQ+rqDg4yEKPqIPTZs2RXl5OW2bwyMd7c0hFlSduCc0kaJRpzdReHoP1OUyyDJv4UXaBTTuOrjKsVKxEC6OTfRRLmso+HiIgo/og0AgQPPmzam7kweKi4vx448/4vtZb0OpUlV7jLX3HDDKcmT+NB55h3+AjfccmFTT4mMAjOvlrOOK2UWLVPMQBR/Rl4ruzo4dO7JdCqnG48eP8dNPP2Hbtm0YPHgwQn7dhh33xDh5O7vKii0isyawD/hPrecTCIDBnewMfuFqavHxEAUf0Rca4MJN165dw8SJE9G9e3e8ePECFy9eRFhYGNzd3THXsz2k4tfbxUEqFmGOZ3stV8s9FHw8RMFH9IWCjzsYhsHx48fh7e2NESNGwMXFBWlpaVi3bh3atm2rOa5HCyss9nGBmaR+H+9mEiEW+7jA1dlKy5VzD3V18hAFH9EXZ2dn3L17l+0yjJpcLsfvv/+O1atXg2EYLFiwAO+++y5MTWvujpzg3hoAsCwqFTJl7QtVCwQvW3qLfVw0P2foKPh4iIKP6IuTkxNiY2NffSDRumfPnmHz5s346aef0KVLF/zwww/w9vau8/6IE9xbw9XZChvi0hB7JxcCvJycXkEsUEMkEmNwJzvM8WxvFC29ChR8PETBR/SFujr17+HDh1i7di127doFPz8/REVFoUePHq91LldnK2ya0Af5JXKEX85E6tNiFMkUePrnA5TnPETo/z42+IEs1aHg4yEKPqIvzs7OFHx6cuHCBaxatQoxMTGYNm0arl27Bmdn7UwrsDE3xcyB7TTfJycLMXHiUtiYf66V8/MNDW7hIalUCrncsFdPJ9zQrFkz5ObmQqmsfiks0jBqtRp//PEHBgwYgHfeeQf9+/dHeno6VqxYobXQq46rqyvS09NRVFSks9fgMmrx8RC1+Ii+iMVi2NraIisrS6cfxMbmxYsX+O2337B69WpYWlrik08+QUBAAMRi/XwkSyQS9OzZE5cuXYKXl5deXpNLqMXHQxR8RJ/oPp/25OTk4KuvvkLr1q0RFRWFbdu24eLFiwgKCtJb6FVwd3fHhQsX9PqaXEHBx0MUfESfaLHqhktNTcWMGTPQqVMnZGVlISEhAYcPH8bAgQPrPEpT29zc3JCYmMjKa7ONgo+HKPiIPlGL7/UwDIP4+Hj4+/tj4MCBaN68Oe7cuYPNmzfDxcWF7fI0LT6mtkl+Boru8fGQqakpBR/RGwq++lEqlQgPD8eqVavw/PlzLFiwAPv27YOZWfV737HF2dkZIpEI6enpaNOmDdvl6BUFHw+JxWKo1WoolUq93xcgxsfZ2Rk3btxguwzOKy4uxrZt2/Djjz+iZcuWWLJkCfz8/CAUcrNjTSAQaFp9xhZ83PwXIbUSCAQ0pYHoDbX4apeZmYlPP/0UrVu3RmJiIvbt24eEhASMGjWKs6FXwVjv83H7X4XUiO7zEX2h4Kve1atX8d5778HV1RXl5eVISkpCaGgo+vXrx3ZpdWasIzsp+HiKgo/oi5OTEzIzM41yEMQ/MQyDY8eOYejQofD19UW3bt1w//59rF27lpfdhb1798a1a9eMrveIbhDxFAUf0Rdzc3OYmpri2bNnsLa2ZrscVsjlcuzduxerV6+GQCDAJ598guDgYJiYmLBdWoOYm5ujQ4cOSElJ4VVLtaGoxcdTFHxEn4y1u7OgoADfffcd2rRpg5CQEKxevRopKSl4//33eR96Fdzc3Iyuu5OCj6co+Ig+Gdsk9gcPHuCDDz5A+/btcffuXRw7dgzHjx/HsGHDWJtwrivu7u5GN8CFgo+nKPiIPhlLiy8xMRGBgYHo168fzM3NcePGDezcuROurq5sl6Yzxtjio3t8PEXBR/TJkINPpVLh8OHDWLVqFR4/foyPP/4Yv/zyC8zNzdkuTS9cXFyQl5eH3Nxc2NnZsV2OXlDw8RQFH9EnZ2dnJCcns12GVr148QK//vorVq9ejaZNm2LhwoV4++23jW5RCKFQiL59++LixYvw9fVluxy9oK5OnqLgI/pUMaXBEGRnZ+PLL79E69atcfz4cezYsQMXLlxAYGCg0YVeBWO7z0fBx1MUfESfDKGr8/bt25g+fTpcXFyQm5uLM2fO4NChQxgwYIDBDVipL2O7z2eclzcGgIKP6JOzszMvg69ih4SVK1fi0qVLmDt3Lu7evWs097Lqys3NDRcvXoRareb8MmvaQMHHU7RWJ9EnW1tblJSUoKysjHO7DFRHoVAgPDwcK1euRGlpKRYsWICwsDBe1M4GOzs72NjY4M6dO+jcuTPb5egcBR9PUYuP6JNAIICjoyOePHmCdu3asV1OjYqKirBt2zasXbsWbdu2xTfffAMfHx+jaMU0VMV9PmMIPvpt4CkKPqJvXO7uzMjIwMKFC9GmTRtcunQJBw4cQFxcHKe3BeIaY7rPR78RPEXBR/SNiyM7r1y5ggkTJqBHjx5QqVS4fPkyfv/9d/Tp04ft0njHmEZ2UvDxFAUf0TeujOxkGAZHjx6Fl5cX/P394erqigcPHmD16tVo1aoV2+XxVo8ePXDv3j2UlpayXYrO0T0+nqLgI/rm7OyMR48esfb6crkce/bswapVqyCRSLBgwQIEBQUZzGLRbDM1NUX37t2RlJSEQYMGsV2OTlGLj6co+Ii+sdXVmZ+fj2XLlqF169YICwvDunXrcOXKFbz33nsUelpmLBvTUvDxFAUf0Td9d3Xev38f8+bNQ4cOHXD//n2cPHlS08Vp7BPOdcXNzc0o7vNR8PEUBR/RN32N6jx//jwCAgLg7u4OS0tL3Lx5Ezt27EC3bt10/trGrmKAC8MwbJeiU3SPj6dMTU0p+IheOTo6IisrCyqVCiKRSKvnVqlU+OOPP7By5UpkZWXh448/xq+//mo0OyRwRevWraFUKpGZmYkWLVqwXY7OUPDxFLX4iL6ZmJigadOmyMnJgaOjo1bOWVpaip07d2LNmjWwtbXFJ598grffflvrwUrqRiAQaO7zUfARzqHgI/qWVyKHtUcQPj14E5LGmbCQiuHSzAKBvZ1hY25ar3NlZWVh/fr12Lx5MwYMGIBff/0V/fv3p3t3HFBxn2/cuHFsl6IzFHw8RcFH9CUloxA/x6Uh/m4uyjt44XSGHEAOAEAqzsKa6Lvw7GSHOYPao0cLq1rPdevWLaxevRoHDhzAu+++i3PnzqFDhw66fxOkztzd3fH111+zXYZOUfDxFAUf0YfdielYFpUKmVIFhgEgrPyRIVOqAQAnbmUj4W4eFvu4YIJ760rHMAyD2NhYrFq1CsnJyZg3bx7u3r0LW1tbPb0LUh99+/bFlStXoFAoIJFI2C5HJyj4eIqCj+jay9C7jTKF+pXHMgxQplBhWdRtAMAE99ZQKBTYt28fVq1aBZlMhgULFmD//v2QSqW6Lp00gIWFBVq3bo3r16+jV69ebJejExR8PEXBR3QpJaMQy6JSqw09RcFjPNk+D41dPGDr/0ml58oUanwbeRvX4iIRsmEF2rdvj6VLl2LkyJG0WDSPVNznM9Tgo99EnqLgI7r0c1waZEpVtc8VnNgEU8ea78uVlSsR/VSIQ4cO4dSpU/D19aXQ4xlDX8GFfht5ioKP6EpeiRzxd3NR3Rzm0lvxEEobQ9qqR40/LxAKUWrRGq06dtVhlUSXDH2LIgo+nqLgI7oSnlz9epxq+QsUnt6DpkOmvvIcAgDhl7m1hRGpu65du+LJkyd49uwZ26XoBAUfT0kkEiiVSqjVrx54QEh9pGYVQa6s+ntVmLAL5j28Ibawe+U5ZEo1Up8W66I8ogcikQi9e/fGxYsX2S5FJ2hwC08JBAJIpVLI5XKYmZmxXQ7hEbVajYKCAuTk5CA7Oxs5OTmV/pxo0hNoUnlfu/LsB5A9SoHj5B/r/DpFMoWWKyf6VLFu5/Dhw9kuReso+HisoruTgo/IZLIqAVZdqOXk5CAvLw8WFhawt7eHvb09HBwcNH/u3bs38osdkZz/j/P/eR3K59nI3DAZAMCUywBGjad582sMQwupYc4BMxZubm7YsmUL22XoBAUfj9F9PsPFMAyePXtWa4D9/c8ymazaIHNycsIbb7xR6XFbW9ta97Fj4u/jRvTdSt2d5j2Ho3HngZrviy4egPJ5NqyHz632HGKBGq0s6eOFz9zc3DB16lQwDGNwS8nRbyaPUfDxi1wuR25ubp2CLDc3F40bN64SZg4ODujRo0eVgLO0tNTah9O43s5YE3230mNCiRSQ/DXxXCCRQiA2gaiRZbXnUKnU+HriCJxy64WgoCCMGjUKTZo00Up9RD8cHR1hbm6OtLQ0g1tWjoKPxyj42MUwDJ4/f16nIMvOzsaLFy9gZ2dXJbQcHBzQvXv3So/b2dnB1LR+Cz9ri625KQZ1tMPJ29nVTmkAAKsB42v8eYEAGN7dCd/fu4U//vgDe/fuxZw5c+Dt7Y2goCD4+vpS9zxPVNznM7TgEzCGvuOgAcorkSM8OROrd4TAtbcbnB1sXnuVfFKZQqHQhNWrgiw3NxdSqbTaIKvuz02bNuVNl1FKRiGCtp6HrA7Llf2TmUSE0BnucHW20jyWn5+PgwcPIiQkBElJSfDz80NwcDC8vb1r7XYl7Fq9ejUePHiA9evXs12KVlHw8cjfV8kHUOkejFQsBAPUeZV8Y8EwDIqKiuoUZDk5OSguLoatrW2dgsze3t5g151Uq9Xwmvk1Hln3glpY90EqZhIhFvt0rrJQ9d9lZWUhPDwcoaGhuHXrFsaMGYPg4GAMHjwYYjF1QnHJuXPn8OGHHyIpKYntUrSKgo8nqqySXwOBAJCKRdWukm8oFAoF8vLy6jzwQyKR1CnIHBwc0LRpU1peC8Dnn3+O06dPY+rynfjh5H2d/d5lZGRg3759CAkJwZ9//omAgAAEBwfjrbfeon8HDigrK4OtrS3y8vIMqnuago8H6rNKfoW6XHlzBcMwKCkpqXOr7Pnz57CxsalTkNnZ2aFRo0Zsv0Ve2bhxI9auXYuzZ8/C1tYW1zILsSEuDbF3ciHAX1sRAX/1NAzuZIc5nu0rdW/W1/379xEaGoqQkBAUFBQgMDAQwcHB6NevH2+6iA1R3759sXbtWnh4eLBditZQ8HFcSkYhgrcmokzx14LBf66qvDMyoyxHkzd8YO09q9Lj1d1r0RelUon8/Pw6t8oEAkGlkYu1hZq1tTVEIpHe35MxiIiIwIwZM3DmzBm0a9eu0nP5JXKEX85E6tNiFMkUsJBK4OLYBON6af/e8q1btzQhWF5ejuDgYAQHB8PV1ZVCUM/mzZuHtm3b4t///jfbpWgNBR/HzdiVVOvoOnW5DJk/TYB94NeQtuxW6TmBABjexQGbJvTRSi2lpaV1DrJnz56hadOmdWqV2dvbo3Hjxlqpkby+S5cuwcfHB5GRkejXrx/b5QB42Rtw9epVTQhKpVJNCLq4uLBdnlHYvXs3Dh8+jH379rFditZQ8HFYXokcHitOVbtuYoWS6zF4fmYvms/aVu2VsKlYiHOfDan2ilylUqGgoKDOw/HVanWdW2U2NjY0UIFHHjx4gLfeegubNm3CqFGj2C6nWgzD4MKFCwgJCcG+fftgb2+PoKAgBAUFoW3btmyXZ7Du3buHoUOH4tGjR2yXojUUfBy2Kf4+1vxjBY1/ytr7BaQtutY4r0osUKO3SRaaP79VJcgKCgpgaWlZ5+H45ubm1M1kgPLz8+Hh4YEPPvgAc+dWvxIL16hUKpw5cwYhISHYv38/2rRpg+DgYAQGBsLZ2Znt8gwKwzCwtbXFjRs34OjoyHY5WkHBx2EfhV7BoatPanxe+TwHjzdNQ/OZWyCxalbjca2ZbIyyL6wSZLa2tpBIaD1FYyaTyTB06FD0798f33//PdvlvBalUolTp04hJCQEhw4dQrdu3RAcHIxx48bB3t6e7fIMgo+PD2bMmIExY8awXYpWUPBx2JRfL+FUak6NzxeeDYEs/SqajV9e63m8XOyxfWJfbZdHeE6tViM4OBhCoRB79+41iOkDcrkcJ06cQEhICCIjI9G3b18EBwfj7bffhrW1Ndvl8dY333yDsrIyLF9e+2cNX/D/N92AWUhrv0dWeuMUzLsNqcN5qFVHqvr000+RlZWFnTt3GkToAYCpqSn8/f2xZ88ePHnyBDNmzMDRo0fRpk0b+Pn5Yffu3Sgupn0C68vd3d2gdmQ3jN92A+XSzAKm4ur/iWSZt6EqyUcjl7dqPYdULISLIy0OTCr76aefEBkZiUOHDhns6jONGjVCYGAgwsPDkZGRgeDgYISEhMDZ2Rnjxo1DWFgYXrx4wXaZvNCvXz8kJSVBpVK9+mAeoODjsHG9a75JX3ojBo069ofQtPbJ2QyAcb3oZj/5y6FDh/C///0PUVFRRtP9Z2FhgQkTJuDIkSN4+PAhRo4ciS1btqB58+YYP348IiIiIJfL2S6Ts5o2bQonJyfcvHmT7VK0goKPwypWya9uIKXNiHmw9V9Q688LBC9X1KCFq0mFxMRETJ8+HYcPH0abNm3YLocV1tbWmDp1Kk6ePIk7d+7Aw8MDP/zwAxwdHTFlyhScOHECSqWS7TI5x83NDYmJiWyXoRUUfBw317M9pOLXW6VEKhZhjmd7LVdE+CotLQ1vv/02du7ciT59tLOoAd85ODhgzpw5SEhIwLVr19C9e3csWbIEzZs3x+zZsxEfH28w3XsNZUj3+Sj4OK5HCyssHNYOUJXX6+cYhRwTXZuwslwZ4Z68vDyMHDkSX3/9NXx9fdkuh5OcnZ3x8ccf48KFC0hMTETLli0xf/58tGzZEh999BESExNhzIPgDanFR9MZOI5hGEyZMgWPJC3wpNmbkCvVdVolf6RjGX7/7xxER0eja9eu+iuYcE5ZWRm8vLzg6emJ7777ju1yeOf27duaJdPkcrlmtZiePXsa1YIOCoUCTZs2xZMnT2BhYcF2OQ1CwcdxW7ZswY8//ogLFy7gQaGyXqvk7927FwsXLsSpU6fQqVMn1t4DYY9KpUJgYCDMzMywa9cug5m2wAaGYXDt2jWEhIQgNDQUEokEwcHBCAoKQpcuXdguTy8GDBiAr7/+Gl5eXmyX0iAUfBxWsWjwmTNnKgVXfVbJ37lzJ5YsWYK4uLgqq+0Tw/fRRx8hJSUFx44dg6kpDXLSFoZhcOnSJU0I2tjYaELQkP+fLVy4EFZWVli8eDHbpTQIBR9H5eXloU+fPli1ahUCAgIadK4tW7bgu+++Q3x8PFq1aqWlCgnXrV27Flu3bsXZs2dhZWXFdjkGS61W4+zZswgJCUF4eDhatWqFoKAgvPPOO2jRogXb5WlVeHg4fvvtNxw+fJjtUhqEgo+DVCoVfHx84Orqih9++EEr51y/fj3WrFmD+Ph4WsTXCOzfvx/z58/HuXPn0LJlS7bLMRpKpRKxsbEIDQ3FwYMH0aVLFwQFBSEwMBAODg5sl9dgmZmZ6NWrF7Kzs3l9f5OCj4O+/PJLJCQkIDo6Wqtb+6xatQqbN29GfHy8wayyTqo6d+4cxowZg+PHj+ONN95guxyjVV5ejhMnTiA0NBQRERHo06cPgoODMXbsWF4vHODk5IQzZ87weh4o3enmmMjISOzYsQMhISFa389uwYIFmDRpEry8vJCTU/Pi14S/7t69i7Fjx+K3336j0GOZiYkJ/Pz8sGvXLjx9+hSzZ8/G8ePH0aZNG/j6+mLXrl0oKipiu8x6c3d35/20Bgo+Dnnw4AGmTJmC0NBQNGtW8zZDDfHFF18gMDAQQ4cORX5+vk5eg7AjJycHI0eOxLJlyzBixAi2yyF/Y2ZmhoCAAISFhSEzMxPjx49HWFgYWrRogbFjx2Lfvn28WTfUzc2N9xPZqauTI8rKytC/f39MmjQJ8+fP1+lrMQyDRYsW4eTJk4iJiaGBDwagtLQUQ4YMwfDhw/Hf//6X7XJIHT179gyHDh1CSEgILly4gJEjRyI4OBgjRozg7CjchIQEfPrpp7xu9VHwcQDDMJg6dSpevHiB33//XS83jRmGwb///W+cO3cOJ0+e5P2EVGOmUqkwduxYWFlZYefOnbwedGDMcnJysH//foSGhuLatWsYPXo0goKC4OXlxakNo0tLS2Fvb4+CggLOhvOrUPBxwNatW7F27VpcuHAB5ubmentdhmEwb948zTwvfb420Q6GYfDBBx8gNTUVUVFRMDExYbskogWPHz9GWFgYQkNDkZaWhoCAAAQHB2PAgAEQiV5v7V5teuONN7Bp0ya4ubmxXcproeBjWVJSEkaOHInTp0/DxcVF76+vVqsxc+ZMpKWlITIyEo0a1b7NEeGWlStX4rfffsPp06dhaWnJdjlEBx4+fIh9+/YhJCQE2dnZCAwMRHBwMNzd3Vlr3c+aNQtdunTBhx9+yMrrNxQNbmFRfn4+xo0bh40bN7ISegAgFAqxefNmtGjRAmPGjIFMJmOlDlJ/oaGh+PHHHxEVFUWhZ8DatGmDzz77DFeuXEFsbCxsbGwwZcoUtGnTBp9++ikuX76s98Wz+T6yk1p8LFGpVPD19UW3bt2wcuVKtsuBSqXChAkTUFRUhAMHDvC2795YnD59GgEBAYiOjoarqyvb5RA9YxgG169f1yyeLRKJEBQUhODgYL0sSn/79m34+fnh/v37On8tXaDgY8lXX32FuLg4xMTEaH2+3utSKBQIDg6GSqVCWFgYp26ok7/cvn0bnp6e2L17N4YNG8Z2OYRlDMMgKSkJoaGhCA0NhZWVlWbd0PbtdbMfp1qthrW1Ne7duwc7OzudvIYuUfCxICoqCjNmzEBSUpLO5uu9rvLycowbNw5SqRR79+7lTCiTl7KysvDmm2/iq6++wqRJk9guh3CMWq3GuXPnNOuGOjs7Izg4GO+8847Wl64b4jMa3UdNh6CpE4pkSlhIxXBpZoHA3lUXy+caCj49e/jwIdzd3bF//3689dZbbJdTLblcjtGjR8PW1ha//vorJ0aREaCkpASenp4YNWoUvvzyS7bLIRynVCoRHx+PkJAQHDhwAC4uLggODkZgYGCDLrhTMgrxc1waom89hQAMVPjr86FiezTPTnaYM6g9erSwavgb0QEKPj0qKyuDh4cHJk6cqPNJ6g1VVlYGPz8/tGrVCtu2baN93FimVCoxZswYODg4YNu2bTRXj9RLeXk5oqOjERISgoiICPTq1QtBQUEICAiAjY1Nnc+zOzEdy6JSIVOq6rQh9mIfF0xwb93wN6BlFHx6NHXqVJSWluptknpDlZaWYuTIkejatSs2bNjAi5oNEcMwmD17NtLT0xEREUH3XkmDlJWV4ejRowgNDcWxY8fg4eGBoKAgjBkzptbRwS9D7zbKFOoaj/knM4kQi306cy78KPj0ZNu2bVizZo3eJ6k3VHFxMby9vdGvXz+sXbuWwo8Fy5cvR2hoKBISEtCkSRO2yyEGpKSkBBEREQgNDUVsbCyGDBmC4OBg+Pn5oXHjxprjUjIKEbw1EWUKleaxouQIlF6PQXluOhp3HgRbv4+rfQ0ziQihM9zh6myl67dTZxR8epCcnIyRI0ciISGBtfl6DVFYWIihQ4diyJAhWLFiBYWfHu3duxeLFi3C+fPn0bx5c7bLIQassLBQs27o+fPnK60b+mHYDZy8nV2pe/PFnXOAQICyh5fBKMprDD6BABjexQGbJvTR0zt5NQo+HcvPz0efPn3www8/YNy4cWyX89oKCgowZMgQ+Pv7Y+nSpWyXYxRiY2MRFBSE2NhYvczNIqRCbm4uDhw4gJCQEKSkPoDlxJ/ACKof5PYsYRdURXk1Bh8AmIqFOPfZEM6M9qQRCzqkUqkwfvx4BAQE8Dr0AMDa2honT57EwYMH8e2337JdjsG7efMmgoODERoaSqFH9M7Ozg4zZ85EbGwsPt9yCMIG9vIIAIRfztROcVpAk7R06L///S/KysqwfPlytkvRCjs7O0RHR8PT0xOmpqZYuHAh2yUZpCdPnsDX1xerV6/G4MGD2S6HGLknLwBVA9tIMqUaqU+LtVRRw1Hw6UhUVBS2b9+OpKQkg5oE3qxZM8TExGDQoEEwMTHh/LQMvikuLoavry9mzJiB8ePHs10OISiSKbV0HoVWzqMNhvOJzCEPHz7E5MmTsX//fs6tzKINTk5OOHXqlCb8Zs+ezXZJBkGhUCAwMBD9+vXDokWL2C6HEACAhVQ7MWEh5c40HAo+LSsrK0NAQAAWLVrE2ZVZtKFly5aIiYmBp6cnTExMMHXqVLZL4rWKuXoikQg///wzjZwlnOHSzAKm4izIlZXn7zFqFVDxxajBKMsBoQgCYdVBMFKxEC6O3JmKQ8GnZfPmzUPHjh2Noguwbdu2iImJweDBg2FiYoL33nuP7ZJ4a9myZbhy5Qri4+MNqmuc8N+43s5YE323yuPPz4bg+dnfNd+X3oyFpce7sBpQtYueATCul7Muy6wX+h+mRdu2bcP58+dx8eJFo7li79ChA06ePAkvLy+YmJggKCiI7ZJ457fffsP27dtx/vx5Xi1uQIyDrbkpBnW0qzKPz2rA+GpD7p8EAmBwJzvOTGUAKPi0Jjk5GYsWLcLp06eN7sOrc+fOOH78OIYNGwYTExO8/fbbbJfEGzExMVi4cCHi4uIM8n4wMQxzPdvj9L28Siu31JVULMIcT91sj/S6aB6fFnBhJ3W2de/eHUePHsWsWbNw5MgRtsvhhevXr+Pdd99FWFgYOnfuzHY5hNSoRwsrLPZxgVRcv8h4uVanC6eWKwMo+BqsYufysWPH8n6SekO98cYbiIiIwJQpU3DixAm2y+G0zMxM+Pr6Yt26dRg4cCDb5RDySuPdWsH+8WmIGBVedSdHIHi5RicXF6gGKPgabOnSpSgtLTWYSeoN1a9fPxw6dAgTJkxAbGws2+VwUlFREXx9fTFv3jwEBwezXQ4hdRIaGoq8xEMImeGO4V0cYCoWVmkBSsVCmIqFGN7FAaEz3DkZegCt1dkgR48exbRp05CUlARHR0e2y+GU+Ph4BAYG4sCBAwY9raO+FAoFfHx80LFjR6xfv95oBkERfsvKykKPHj1w5MgR9O3bFwCQXyJH+OVMpD4tRpFMAQupBC6OTTCuF+3AbrAqdlIPDw/HgAED2C6Hk6Kjo/Gvf/0LERERcHNzY7sc1jEMg8mTJ6OgoAAHDhygaQuEFxiGwdtvv42uXbti2bJlbJejFdTV+RpkMhnGjRuHzz//nEKvFkOHDsXOnTsxatQoJCcns10O67755hvcunULv//+O4Ue4Y09e/bg/v37+PLLL9kuRWuoxfcapk2bhuLiYoSEhFBXVR388ccfmDlzJo4fP44ePXqwXQ4rduzYgWXLluHcuXNwcHBguxxC6uTJkyfo2bMnjh49it69e7NdjtbQZWc9bd++HefOnTOqSeoNNXr0aJSXl2PEiBGIjo42um12jh8/ji+++ALx8fEUeoQ3GIbBzJkzMWvWLIMKPYCCr16Sk5Px+eefG+Uk9YYKDAyEQqGAt7c3Tp06hU6dOrFdkl5cvXoV7733Hg4ePGg075kYht9++w1//vkn9u/fz3YpWkfBV0cFBQVGP0m9of71r3+hvLwcQ4cORVxcHNq1a8d2STr1559/ws/PDxs2bICHhwfb5RBSZ48fP8bChQtx4sQJmJiYsF2O1lHw1YFaraZJ6loyadIklJeXw8vLC/Hx8WjVqhXbJelEYWEhfHx88O9//5t+ZwivMAyD6dOnY+7cuejZsyfb5egEBV8dLF26FCUlJTRJXUtmzJiB8vJyDBkyBPHx8XB25s6q7dpQXl6OsWPHwsvLCx9//DHb5RBSL7/88guePn2KL774gu1SdIZGdb7CsWPHMHXqVJqkrgOrVq3C5s2bER8fbzB/twzD4P3330dpaSnCwsIgElXdm4wQrsrIyECvXr0QExMDV1dXtsvRGWrx1SI9PR0TJ05EeHi4wXwwc8mCBQsgl8vh5eWFuLg42Nvbs11Sgy1ZsgRpaWmIiYmh0CO8wjAMpk2bhvnz5xt06AEUfDWiSer68cUXX0Aul2Po0KGIjY2FjY0N2yW9ti1btiA0NBTnzp1Do0aN2C6HkHrZtm0b8vPz8dlnn7Fdis5RV2cNpk+fjqKiIpqkrgcMw2DRokU4efIkYmJiYGVlxXZJ9RYVFYWpU6fi9OnTaN+eW3uPEfIqjx49Qp8+fRAbG4tu3bqxXY7OUYuvGjt27MCZM2dokrqeCAQC/O9//4NcLsfw4cNx8uRJWFhYsF1WnSUnJ2PSpEk4fPgwhR7hHYZhMHXqVCxYsMAoQg+gFl8Vly9fxvDhw5GQkECbg+oZwzCYN28eUlJScOzYMV4sEpCeng4PDw+sX7+edp4nvLRp0yb88ssvOHv2rNGsIUvB9zcFBQXo3bs3vv/+ewQGBrJdjlFSq9WYOXMm0tLSEBkZyel7Zc+ePYOHhwdmzZqFDz/8kO1yCKm3hw8fom/fvkhISECXLl3YLkdvKPj+n1qthp+fH1xcXLB69Wq2yzFqarUakyZNQlZWFg4fPgypVMp2SVXI5XJ4e3ujT58+WLVqFdvlEFJvarUaXl5e8PHxwcKFC9kuR68o+P7fN998g+joaJw6dQoSiYTtcoyeUqnEhAkTUFxcjAMHDsDUlDsbW6rVaowfPx5KpRKhoaEQCml3L8I/69evx549e3DmzBmjm3pDwQeapM5VCoUCQUFBUKvVCAsL48wFyeeff44zZ84gOjqak61RQl7l/v37cHNzw9mzZ41y8XSjv1StmKQeEhJCoccxEokEISEhUKlUmhYW2zZu3IiDBw/ijz/+oNAjvKRWqzF58mR88cUXRhl6gJEHX8Uk9c8++4wmqXOUiYkJwsLCUFRUhEmTJkGlUrFWy+HDh7F06VIcPXqU1xPtiXH76aefoFarMX/+fLZLYY1Rd3VOnz4dhYWF2LdvH83X47iysjL4+fmhVatW2LZtm97vq126dAk+Pj6IjIxEv3799PrahGjLvXv38Oabb+L8+fPo0KED2+WwxmhbfBWT1Hfs2EGhxwNmZmY4fPgw0tLSMHfuXOjzeu3BgwcYPXo0tm/fTqFHeEulUmHy5MlYsmSJUYceYKTBd/nyZXz22Wc4cOAAmjRpwnY5pI4aN26MyMhIXL16FR999JFewi8/Px8jR47Ef/7zH4waNUrnr0eIrvz4448QCoX44IMP2C6FdUbX1VlQUIA+ffpgxYoVNEmdpwoLCzF06FAMGTIEK1as0FmLvaysDMOGDYOHhwdWrFihk9cgRB/u3LkDDw8PXLhwAe3atWO7HNYZVfCp1Wr4+/ujU6dONEmd5woKCjBkyBD4+/tj6dKlWj+/Wq1GUFAQxGIx9uzZQ3P1CG+pVCq89dZbmDBhAubOnct2OZxgHAuz/b9vv/0WRUVFdPVuAKytrXHy5EkMHjwYpqam+M9//qPV8y9cuBA5OTk4ceIEhR7htdWrV0MqlWL27Nlsl8IZRhN8x44dw+bNm5GUlMSZidCkYezs7BAdHQ1PT0+YmppqbdmldevWISoqCmfPnuXUijGE1NetW7ewYsUKXLp0iS7g/sYogi89PR2TJk1CWFgYTVI3MM2aNUNMTAwGDRoEExOTBs9NOnjwIFasWIGzZ8/C2tpaS1USon9KpRKTJk3Ct99+izZt2rBdDqcYfPBVTFL/9NNPaZK6gXJycsKpU6cwaNAgmJqaYtasWa91nsTERMyYMQPHjh1D69attVskIXq2cuVKWFpaYubMmWyXwjkGP7hlxowZePbsGU1SNwIPHjyAp6cnvv76a0yZMqXSc3klcoQnZyI1qwhFMiUspGK4NLNAYG9n2JibIi0tDQMGDMD27dvh4+PD0jsgRDtu3LiBwYMHIykpCa1atWK7HM4x6BbfL7/8gtOnT9NO6kaibdu2iImJweDBgyGRSPDee+8hJaMQP8elIf5uLgBArlRrjpeKs7Am+i7ebG2J05v/g2+++YZCj/CeQqHApEmT8N1331Ho1cBgg+/KlSv49NNPkZCQQJPUjUiHDh1w4sQJDB06FMlFjXA8uzFkShWq69eQ/X8Ixt3Lh8jrIzRyddVztYRo34oVK2BjY4Np06axXQpnGWRX57Nnz9C7d28sX74c77zzDtvlEBas2H8OG85nQSD5a1Rm1p7PIX9yBwLhy73HRE1s4DRjs+Z5M4kQi306Y4J7a32XS4hWXLt2DV5eXrh8+TJatGjBdjmcZXAtPrVajQkTJmD06NEUekYqJaMQO1OKKoVeBWvvWWjSY3i1P1emUGNZVCpcna3g6myl4yoJ0S6FQoGJEydixYoVFHqvYHATO5YtW4aioiJ8//33bJdCWPJzXBpkytfbvkimVGFDXJqWKyJE97777js4Ojpi8uTJbJfCeQbV4jt+/Dg2bdpEk9SNWF6JHPF3c6u9pwcAhXG/ojDuV0isnWA18D1IW1W+r8cwQOydXOSXyGFjTpPXCT9cvXoVP//8M65cuUID+eqAN8H3quHojx49wsSJE7Fv3z6apG7EwpMza3yu6eDJkNi0gEAkQentBOTsXwrHyesgaVr590UAIPxyJmYOpMV8CfeVl5dj4sSJ+OGHH+Dk5MR2ObzA+eCry3D0Ae1tkPTbMixcuBADBw5kq1SiYyqVCsXFxSgqKqryVfF4RJ415Grban/etHknzZ/Nu3uh9FY8yu4nQdLHv9JxMqUaqU+LdfpeCNGWb7/9Fi1btsT777/Pdim8weng252YjmVRqa8cjh59OxvCPlNg37+Hniskr8IwDORyeY1B9aqvvx9XVlYGc3NzWFhY1PhVbmZX9+IEAgDV94kWyRTa+QsgRIeSk5OxadMmXL16lbo464Gzwfcy9G6jTKF+9cECIdQCIb47ehsCAWg4uhao1WqUlJTUOaBqCzMAsLS0rBJSTZo0qfS9g4NDjYHWpEkTNG7c+JUL7X4UegWZV59UfT+yEsif3IG0ZXdAKELp7QTIM27A2mt6tecxNzG4cV/EwMjlckycOBFr1qxB8+bN2S6HVzgZfCkZhVgWlVol9JSF2cg/sQHlj1MBsQSNO3mg6dAZmnlZNBz9ZX9/XcLqVceUlpaiUaNGNYZUxZednR3atWtX63H63OHApZkFTMVZlbrEAYBRq1CYsBuKgkxAIITExhl2Y/8DiY1z1ZOoFAjZtAqP95fDz88PPj4+sLe319M7IKRu/vvf/6J9+/b417/+xXYpvMPJCewzdiXh5O3sKt2b2fu+gqiRFWxGzIVaVors0P/AvMdwWPQZpTlGIACGd3HApgl99Fz162MYBi9evHjtkPr7l1KprLUr8FWtroovc3NziEQitv9q6i2vRA6PFaeqBF99mIqFODzVFefjTiIiIgLR0dHo2rUr/Pz84O/vj27dulG3EmHVpUuX4Ofnh5SUFDRr1oztcniHcy2+2oajK59nw6K3HwRiE4jMTWDWpjcUeX9WOkafw9GVSmW1oVTf+1fFxcUwMTF5ZUhZWlqiRYsWtR4jlUqN+kPZ1twUgzraVXvhVBcCATC4kx06tXZCp0mTMGnSJMjlcsTHxyMiIgL+/i8HwlSEYMVegIToi0wmw8SJE/Hjjz9S6L0mzrX4NsXfx5rou9VesRdfiYI88zasR8yFWlaCnNAvYTVgAhp16l/pOKlYiI+Hdax2ODrDMJDJZA0aZFHxJZPJamwx1fR4dcc1adKE5h1qUUpGIYK3JqJMUf9J7GYSEUJnuNfYVc4wDG7evIkjR44gIiICN27cgJeXF/z8/ODr6wsHB4cGVk9I7T7//HOkpaUhLCzMqC9yG4JzwfdR6BUcqmZwAgAo8jKQF7ES5TkPAUaNxt28YOP7UbX/+LbFD+Dw8Fi1ISYSiV4rpP55bOPGjekXj6PqNTjq/73OWp25ubk4evQoIiIicPLkSXTq1An+/v7w9/eHq6sr/X4QrUpMTMSYMWNw7do1uu/cAJwLvim/XsKp1JwqjzOMGo83TkWTniNg0W8s1Ioy5Ef+CImNE5oOnlLleJcmCsxzFVcbXNQ1ZRxeNR2mgkAASMUiLPZxadCI4PLyciQkJCAiIgIRERFQKpWaLtHBgwdDKpW+9rkJKSsrwxtvvIGlS5ciMDCQ7XJ4jXPBV1OLT/XiOTLXjUeLj0IhlDYGALy4ex6FCbvQfNqGKse/3dMJa4J66rpcwnHXMguxIS4NsXdyIcBfcz+Bl13iDF7e05vj2V6rI4EZhsHt27c1XaIpKSkYMmSIpkuUVhci9fXJJ58gIyMDoaGhbJfCe5wb3FLTcHRRI0uILR1QfCUKFm5jwZSXoeR6DCT2baqcQyoWwsWR9uAjgKuzFTZN6IP8EjnCL2ci9WkximQKWEglcHFsgnG9nHUyCEogEKBLly7o0qULPv30U+Tl5eHYsWOIiIjAwoUL0aFDB01rsGfPntQlSmp17tw57NmzB9evX2e7FIPAuRZfbcPRy7MfoCB6CxQ5DwGhCNKW3WHtPRuixlaVjjMVC3HusyG0yDDhJIVCgdOnT2u6RGUymSYEhwwZAjMzM7ZLJBzy4sUL9OzZE8uXL8fYsWPZLscgcC74gJrn8dUFH+fxEePFMAzu3LmDiIgIHDlyBFeuXIGnpyf8/f3h6+tLK3IQfPzxx8jOzsbevXvZLsVgcDL4dDkcnRAuKygo0HSJHj9+HG3bttW0Bnv16kVdokbm9OnTCAoKwvXr12FjY8N2OQaDk8EH6G84OiFcpVAocPbsWU2XaGlpKXx9feHv7w8vLy80atSI7RKJDpWWlqJHjx5YtWoVRo8ezXY5BoWzwQfofzg6IVx29+5dTZdocnIyBg4cqOkSdXauZs1Rwmsffvghnj17hl27drFdisHhdPAB7A1HJ4TLnj17huPHjyMiIgLHjh1Dq1atNF2ivXv3fuUuFoTb4uLiMH78eFy/fh3W1tZsl2NwOB98FfQ9HJ0QvlAqlTh37pymNVhYWKjpEh06dCgaN27MdomkHkpKSuDq6op169bBz8+P7XIMEm+CjxBSN2lpaZoQvHjxIgYMGAB/f3/4+fmhRYsWbJdHXmHu3Ll48eIFfvnlF7ZLMVgUfIQYsMLCQhw/fhxHjhzB0aNH4ezsrAnBvn37Upcox8TExGDSpEm4fv06rKys2C7HYFHwEWIklEolEhMTNaNE8/PzNV2iw4YNg7m5OdslGrWioiK4urpi48aNGDlyJNvlGDQKPkKM1P379zVriV64cAEeHh6a1mCrVq3YLs/ozJw5EyqVCtu2bWO7FINHwUcIwfPnz3HixAkcOXIEUVFRcHR01IRgv379IBKJ2C7RoJ04cQLTp0/HtWvXYGlpyXY5Bo+CjxBSiUqlQmJioqY1mJOTAx8fH/j7+8Pb2xtNmtAC8Nr0/PlzdO/eHdu2bYO3tzfb5RgFCj5CSK0ePnyoCcHz58/jzTff1Gy227p1a7bL471p06ZBJBJh8+bNbJdiNCj4CCF1VlxcrOkSjYyMhL29vaZL1N3dnbpE6+no0aOYPXs2rl+/Ti1pPaLgI4S8FpVKhYsXL2pag0+fPoWPjw/8/PwwfPhwWFhYsF0ipxUWFqJ79+7YuXMnvLy82C7HqFDwEUK04tGjR5oQPHv2LNzc3DRdom3btmW7PM6ZPHkyzMzMsGHDBrZLMToUfIQQrSspKcHJkyc1XaLW1taaLtE333wTYrGY7RJZFRkZiQ8++ADXrl2j+ZMsoOAjhOiUWq3GpUuXNK3BzMxMjBw5En5+fhgxYoTRDd9/9uwZunfvjt27d8PT05PtcowSBR8hRK8yMjI0IXj69Gn069dPs7NE+/bt2S5P595//31YWVlh3bp1bJditCj4CCGsKS0tRXR0NCIiIhAZGQlLS0tNl6iHh4fBdYn+8ccfWLBgAVJSUmjXDBZR8BFCOEGtViM5OVnTGnz06BFGjBih6RJt2rQp2yU2SH5+Prp3747Q0FAMGDCA7XKMGgUfIYSTMjMzERkZiYiICCQkJKB3796aLtGOHTuyXV69jR8/Hvb29lizZg3bpRg9Cj5CCOe9ePECMTExmn0Gzc3NNVMlPDw8IJFI2C6xVgcOHMDnn3+Oq1evolGjRmyXY/Qo+AghvKJWq3HlyhVNl+iDBw8wfPhw+Pv7Y8SIEbC2tma7xEpyc3Ph6uqK8PBweHh4sF0OAQUfIYTnnjx5oukSjYuLwxtvvKHpEu3UqRMEAgGr9QUFBaFFixZYuXIlq3WQv1DwEUIMRllZGU6dOqXpEjUzM9OE4IABA3TSJZpXIkd4ciZSs4pQJFPCQiqGSzMLBPZ2xqmjh7FkyRJcuXIFZmZmWn9t8noo+AghBolhGFy9elUTgvfu3YO3tzf8/f0xcuRI2NjYNOj8KRmF+DkuDfF3cwEAcqVa85xULISaYVB2PwnL3x+C8SPfatBrEe2i4COEGIWnT58iKioKEREROHXqFHr06KFpDXbu3LleXaK7E9OxLCoVMqUKtX6CMgzMTMRY7OOCCe6tG/weiHZQ8BFCjI5MJkNsbKymNSiRSDQhOHDgQJiYmNT4sy9D7zbKFOoaj/knM4kQi306U/hxBAUfIcSoMQyDa9euaUIwNTUVw4YN03SJ2tnZaY5NyShE8NZElClUf/28UoH8ExsgS78KtawEYitHNB30Psza9an0OmYSEUJnuMPV2Upfb43UgIKPEEL+Jjs7G5GRkThy5AhiYmLQrVs3zTJq667IcPJ2dqXuTXW5DEUX9sO8+1CILO1Qdj8JeYd/QPMp6yG2ctAcJxAAw7s4YNOEPtW8KtEnCj5CCKmBTCZDfHw8IiIiEHEiFoIxywDRq0eGPtk+D5Ye76KxS+V5e6ZiIc59NgQ25qa6KpnUgZDtAgghhKukUimGDx+O9evXY9HWP2Bi8urQU5U+g6LgMUzsWlZ5TgAg/HKmDiol9UHBRwghdZCaVYxyVe3HMCol8g6vhHl3L0hsWlR5XqZUI/VpsY4qJHVFwUcIIXVQJFPW+jzDqJF3ZBUgEsN62KxazqPQdmmknij4CCGkDiykNe8NyDAM8qPWQVVaCLu3v4BAVPOxFlJuL6htDCj4CCGkDlyaWcBUXP1HZsHxn6HIz4D9uC8hlNQ8cEUqFsLFsYmuSiR1RKM6CSGkDvJK5PBYcarS0mQAoHyeg8cbpwAiCQRCkeZx6xFzYd51cKVjaVQnN9TcHieEEKJha26KQR3tqszjE1vao9XnR1758wIBMLiTHYUeB1BXJyGE1NFcz/aQikWvPrAaUrEIczzba7ki8joo+AghpI56tLDCYh8XmEnq99H5cq1OF1qujCOoq5MQQuqhYqHpuuzOIBC8bOnR7gzcQoNbCCHkNVzLLMSGuDTE3smFAC8np1eQioVg8PKe3hzP9tTS4xgKPkIIaYD8EjnCL2ci9WkximQKWEglcHFsgnG9nGkgC0dR8BFCCDEqNLiFEEKIUaHgI4QQYlQo+AghhBgVCj5CCCFGhYKPEEKIUaHgI4QQYlQo+AghhBgVCj5CCCFGhYKPEEKIUaHgI4QQYlQo+AghhBgVCj5CCCFGhYKPEEKIUaHgI4QQYlQo+AghhBgVCj5CCCFGhYKPEEKIUaHgI4QQYlQo+AghhBgVCj5CCCFGhYKPEEKIUfk/rSfBGNhI43EAAAAASUVORK5CYII=\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAb4AAAEuCAYAAADx63eqAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABMsElEQVR4nO3dd1iT5/4G8Dsh7CGoIIhMEZAtQxTE3draarXWOqq2ta1Yte5Wre2p1tGt1j1aR9W2Wmedp+6NiiAgQ2QPAUVEhgRI8p4/LBxjEnZ4Mr6f6+I65ychuTk/4M7zvs/gcRzHgRBCCNESfNYBCCGEkNZExUcIIUSrUPERQgjRKlR8hBBCtAoVHyGEEK1CxUcIIUSrUPERQgjRKlR8hBBCtAoVHyGEEK1CxUcIIUSrUPERQgjRKlR8hBBCtAoVHyGEEK1CxUcIIUSrUPERQgjRKlR8hBBCtAoVHyGEEK1CxUcIIUSrUPERQgjRKlR8hBBCtAoVHyGEEK1CxUcIIUSrCFgHIEQVFZZVYt+tHCTll6BEKIKZgQDu1mYYGdAJ7Uz0WccjhDQDj+M4jnUIQlRFTHYx1p1PwYXkhwCASpGk9nMGAj44AH3dLDGljwt87czZhCSENAsVHyH/2hWRgWXHkyAUiVHXbwWPBxgIdLBwsDvG9XBstXyEkJZBlzoJQU3pJaKiWlLvYzkOqKgWY9nxRACg8iNEzdCIj2i9mOxijN4SgYpqsdS/l9w6gvK4M6h6mAHjrn3Q/vVZMl9rqKuDPZN6wKeTeSulJYQ0F83qJFpv3fkUCEVimX8XmLRDm5BRMPF5SeHXCkVirD+fosx4hJAWRsVHtFphWSUuJD+Ue0/PyC0ERq49wTc0U/j1HAecu/sQj8oqlZiSENKSqPiIVtt3K6fZz8EDsC+q+c9DCGkdVHxEqyXll0gtWWgKoUiCpLzSFkpECFE2mtVJtBLHcSgoKEBaTn6LPF+JsLpFnocQonxUfESjcRyH/Px8xMfHIyEhQeo/eTweOgybB1h5Nvt1Iq9dwpKsk+jRowe6d++ONm3atEB6Qogy0HIGohE4jkNeXp5MuSUkJIDP58PT0xOenp7w8PCo/U8rKytsupiGlaeT5V7u5CRiQCJG8eXfIS59hHavfgLwdcDj60g9Tl+Hh5c7VsMo8yoiIiJw69YtODg4oEePHrUfHh4e0NHRkXkNQkjro+IjaoXjONy/f19uwenq6koV2/MFp0hhWSVCvzsrt/iKL+3Gkyt/SP1bm9AxMA97R+rf9AV8XJ3Xv3YPz+rqaty5cwcRERG1H3l5eQgKCqotwuDg4DpzEUKUh4qPqCSO45Cbmyu34PT19eUWnKWlZZNea9LOSJxKLKhzmzJFeDxgkEcHbBwXWOfjHj16hBs3btQW4Y0bN2BhYSE1KvTz84Oenl6TvgdCSMNR8RGmOI5DTk6O3IIzNDSUW3Dt27dv0QyKdm5piKbu3CKRSJCcnCw1Krx37x58fX2lytDOzg48Hq/RuQghilHxkVbBcRyys7PlFpyRkZHce3Dt2rVrtXyN2auzhqEuHwsHd22xvTrLysoQGRlZW4TXrl2Djo6OVBEGBATA2Ni4RV6PEG1FxUdaFMdxyMrKkim4xMREGBsby5RbaxdcXRp8OgMAA13ln87AcRwyMzOlRoVxcXFwc3OTulfo6upKo0JCGoGKjzSJRCJBdna2zOgtISEBpqamcguubdu2rGPXKzanGOvPp+Dc3Yfg4dni9BoGAj6qRWIYP0nHrs/HM9mYWigU4vbt21JlWFJSguDg4Noy7N69OywsLFo9GyHqgoqP1EkikSArK0um4BITE2FmZia34DThj+6jskrsi8pBUl4pSoTVMDPQhbuNKV7r2g7dPLrgypUrcHV1ZR0TAJCfn4/r16/XFmFkZCQ6deokdYnU09MTAgEt2yUEoOIj/5JIJMjMzJRbcObm5nILztzcnHVsJhYuXIji4mKsW7eOdRS5RCIR4uPjpUaFOTk5CAwMlLpEam1tzToqIUxQ8WkZiUSCjIwMmYJLSkqChYWF3IKjXUik3b9/H56enkhNTVWLy7cA8PjxY6nlFNevX0ebNm1kllPo6+uzjkqI0lHxaSiJRIL09HS5BdeuXTu5BWdmpvj4HSJtwoQJ8PT0xLx581hHaRKJRIJ79+5JjQqTk5Ph7e0tVYYODg40cYZoHCo+NScWi5Geni4zi/Lu3bto3769TMF17dqVCq4FREdHY8iQIUhPT4euri7rOC2ivLwct27dklpOwXGcVBEGBgbCxMSEdVRCmoWKT02IxWKkpaXJLTgrKyu5BWdqaso6tkbr168fPvroI4wdO5Z1FKWoWXv5/KgwJiYGXbp0kSpDV1dX8Pl0whlRH1R8KkYsFiM1NVWm4JKTk9GhQ4facnu+4OgdOBt///03vv76a9y8eVNrLgdWVlYiJiZGqgwfP34ss5xCXe59Eu1ExceISCSSW3D37t2DtbW1zFZd7u7uVHAqRiKRwM3NDVu3bkVYWBjrOMwUFBTILKewsbGRGhV6e3vTcgqiMtSm+ArLKrHvVg6S8ktQIhTBzEAAd2szjAzoVLsrvioSiURISUmRW3AdO3aUW3C0JZX6WLduHU6fPo2DBw+yjqIyxGKxzHKKrKwsBAQESJWhjY0N66hES6l88cVkF2Pd+RRcSH4IAFLHxxgI+OAA9HWzxJQ+LvC1M2cTEs+OoklNTUV8fLxUwaWkpMDW1lZuwRkZGTHLS1pGeXk5HBwccP36dXTu3Jl1HJVVXFyMmzdvSpWhiYmJVBF269YNBgYGrKMSLaDSxdfgvRN5gIFA+XsnAs8KLiUlRWaZQEpKCjp16iRTcG5ublRwGm7+/Pl4+vQpVq9ezTqK2uA4DikpKVJFmJSUBC8vL6kydHR01Jr7p6T1qGzxsd4tv6qqSm7Bpaamws7OTm7BGRoaNvt1ifrJycmBj48P0tLStHY3m5bw9OlTmeUUYrFYZjkFzVYmzaWSxSfvfDROVI1H/6yHMOM2JMIyCMxtYNFnAgw7Sx8A2tjz0aqqqnDv3j2ZgktLS4O9vb1Mwbm6ulLBERnvvPMOunXrhrlz57KOojFqzmp8cTmFs7OzVBm6u7vTcgrSKCpZfPJOxJZUCVFyfT9MvAdCp40lKlIjUfj3D+g4cS0E5h1qH6foROyqqiokJyfLFFx6ejocHBzkFhzdbyANFRkZiTfffBNpaWk0e1GJqqqqEBsbK1WGhYWF6N69u9Q+pKpy1BVRTSpXfIVllQj97qzUJBZF7v86DW1Cx8DYPVTq33X5wHyPMmTe/f9Ek4yMDDg6OsotONqfkLSE3r17Y+rUqRg1ahTrKFrl4cOHUsspbt68iQ4dOsgsp2Cxw466zkbXdCpXfBsvpGLl6eR6i09c/hg56yei48TV0G1nJ/1JURU6PY5B/46S2sXeXbp0oYIjSnXw4EF8++23iIiIoAkZDInFYiQmJkqNCjMyMuDv7y81KrS1tVVaBnWZja6tVK74Zu6JxqHb9+t8DCcW4cHeryCwsEG7V6bJfcxwP1usHOWnhISEyCcWi+Hq6oqdO3ciJCSEdRzynCdPnsgspzA0NJQaFfr7+7fI/XtVnI1OpKnczYgSoajOz3OcBIVHfwJ0BGj70uQ6nqe6paMRUicdHR3MmDEDK1asoOJTMW3atMHAgQMxcOBAAM8mzqSmptYe0fTnn38iISEBHh4eUmXo7OzcqNF7Y2ajcxxQUS3GsuOJAEDl14rUasTHcRweHf8ZoicFsBq5CHxdxZcuacRHWCgtLYWjoyMiIyPh5OTEOg5phIqKCkRFRUmNCoVCoVQRBgUFKTzdRN5sdAAQV5Ti0fGfIcyIBt/QDBZ93oWxZ1+pxzR2NjppHpUb8blbm0FfkC/3Hl/Rf9eh+lE2OoxeWmfpGQj4cLehtT6k9ZmammLixIlYvXo1Vq5cyToOaQRDQ0OEhoYiNPT/k+VycnJqJ84sWrQI0dHRcHJykirDrl27gs/nY935FAhFYpnnLfpnA3g6uuj0yS5UFaThwb7F0LVygp6lQ+1jhCIx1p9PkZmNTpRD5UZ8imZ1ip48QO6GiYCOLnh8ndp/b/vKVJh49pN6rL6Aj6vz+tOsKcJEVlYW/Pz8kJGRQWcfapjq6mqZ5RQPHjyAf0gfZPp9BAmk1xNKqoTIXjUaHT9cB922zybTFB75CTqm7WDR9z2px9LfrdajciO+9ib66ONqKbOOT9DGCg7zj9b79Twe0M/Nkn54CDP29vZ4+eWX8euvv2LWrFms45AWpKuri4CAAAQEBGDq1KkAgMLCQizdfx1ZGbKjPVFRLnh8fm3pAYCulRMqs+JkHssDsC8qB+G9ac9XZVPJ7Q6m9nWBgUCn/gfKYSDQwZS+Li2ciJDGmT17Nn7++WeIRHVP1iLqr3379oB5R4gh+zdLUl0Bnr70Xr18fSNIqipkHisUSZCUV6q0nOT/VLL4fO3MsXCwOwx1Gxfv2V6d7nSDmDDXvXt32Nra4tChQ6yjkFagaDY6X9cQXKV0yXGVT8HXk79sgmajtw6VLD7g2dTehYO7wlBXB/XNJubxns2KaqkNqglpCbNnz8aKFStYxyCtwMxA/l0jQVtbcBIxqotya/+t6kE6dJ+b2CL9PK2/u4w2UtniA56V355JPTDIowP0BXwYCKTj8iQi6PA4DPLogD2TelDpEZUybNgw5OXlISIignUUomTPZqPL/jnl6xnAyK0nii/thqRKCGFOAp6mXIfxCxPyAJqN3ppUblanIo/KKrEvKgdJeaUoEVbDzEAX4qIs3PprHa6c/S/reITItXLlSkRERGDPnj2soxAlqmuP4Yas4wNoVmdrUpvik6eqqgq2tra4ceMGLRYmKqmkpAROTk6IioqCg4P8y1tE/YlEIgxYtAcZ1WbgNeGIJEWnyhDlUOlLnfXR09PDqFGjsGvXLtZRCJHLzMwM7777LtasWcM6ClGSpKQkhIaGQhx3HAZ6NBtdHah18QHA+PHjsXPnTqjxwJVouOnTp2Pbtm0oLaWp6ppEIpHg559/Rq9evTBhwgRcPLgTX77mQbPR1YDaF1/37t3B4/Fw/fp11lEIkcvR0RH9+/fHtm3bWEchLSQjIwMDBgzAnj17cO3aNUydOhV8Pr9Rs9HBSSCAhGajM6D2xcfj8WpHfYSoqtmzZ2PVqlUQi2V39yDqg+M4/PrrrwgKCsIrr7yCS5cuoUuXLlKPqW82uoGAD30BH6EOpig5uBgvO0svcCfKp9aTW2pkZGQgKCgIubm50NPTYx2HELl69OiBefPmYfjw4ayjkCbIy8vDRx99hNzcXPz222/w9vau92vkzUZ3tzHFW/7PTmCfPXs2ysrKsHnz5lb4DkgNjSg+AOjTpw9mzZqFYcOGsY5CiFx79+7FmjVrcOnSJdZRSCP9+eefmDFjBsLDw/HFF1+02Bvs4uJiuLu74+TJk/Dz82uR5yT105ji++WXX3DixAns37+fdRRC5BKJROjcuTP27duHoKAg1nFIAxQWFmLq1KmIjY3Fjh070L179xZ/jY0bN+LPP//EuXPnGnXoLWk6tb/HV2PkyJE4ffo0ioqKWEchRC6BQIDp06fTOX1q4ujRo/Dx8UHHjh0RFRWllNIDgA8//BBFRUU4cOCAUp6fyNKYER8AjBo1Cv369cPkyZNZRyFEridPnsDJyQkxMTGws7NjHYfIUVJSglmzZuHs2bPYtm0b+vbtq/TXPHv2LD788EMkJCTAwMBA6a+n7TRmxAcAEyZMoNmdRKW1adMG48ePx9q1a1lHIXKcPXsWPj4+0NHRQWxsbKuUHgD0798fvr6+dDWglWjUiK+6uhqdOnXClStX4OJCuyAQ1ZSWlobu3bsjIyMDJiYmrOMQAE+fPsX8+fNx4MABbN68GYMHD271DKmpqQgODkZcXBxsbGxa/fW1iUaN+HR1dTF69GjawoyoNGdnZ/Tu3Rvbt29nHYUAiIiIQLdu3VBYWIjY2FgmpQcAnTt3xgcffIDPP/+cyetrE40a8QHArVu38PbbbyMlJYVmSBGVdfnyZbz//vu4e/cu+E3Y1Jg0X2VlJRYvXoxff/0Va9euxciRI1lHQklJCdzc3HDkyBEEBtKG1cqicb9x/v7+0NfXx9WrV1lHIUSh0NBQmJub4+jRo6yjaKWYmBh0794dd+7cQUxMjEqUHvBsU/MlS5Zg5syZtP+wEmlc8dEWZkQd8Hg8OqGdAZFIhOXLl2PgwIGYNWsWDh8+DGtra9axpLz//vsoLy+nMxyVSOMudQJAVlYWunXrhtzcXJoaTFRWdXU1nJ2dcfjwYfj7+7OOo/Hu3r2Ld999F8bGxti2bRvs7e1ZR1Lo4sWLGD9+PBITE2FkRHt5tjSNG/EBgL29PXx9fXHs2DHWUQhRSFdXF5988glNYVcyiUSC1atXIzQ0FOPGjcOpU6dUuvQAoHfv3ujevTt+/PFH1lE0kkaO+ABg+/btOHjwIA4fPsw6CiEKPX78GM7Ozrhz5w5sbW1Zx9E4mZmZeP/991FRUYEdO3bA1dWVdaQGy8jIQEBAAGJiYtCpUyfWcTSKRo74AGDEiBG4cOECCgsLWUchRCELCwuMGzcO69atYx1Fo9QcHxQYGIhBgwbh8uXLalV6wLNzHCdPnoz58+ezjqJxNHbEBwBjx45FaGgopk6dyjoKIQqlpKSgZ8+eyMjIgLGxMes4aq/m+KCcnBzs3LmzQccHqaqysjK4ublh//796NGjB+s4GkNjR3zAsy3MfvvtN9YxCKmTi4sLQkJC6Ge1BezZswd+fn7o1q0bbty4odalBwAmJiZYvnw5Zs6cCYlEwjqOxtDoEZ9IJIKdnR3Onz8PNzc31nEIUejChQuYNGkSEhMTaUF7Ezx69AhTp07F7du38dtvvyntJAUWJBIJevTogenTp2PcuHGs42gEjf4NEwgEGDt2LK3pIyqvd+/eMDY2xokTJ1hHUTvHjh2Dj48PbGxsEB0drVGlBwB8Ph+rVq3CggULUF5ezjqORtDoER8A3L59G8OGDUNaWhq9kyYqbdeuXdi2bRvOnDnDOopaYHF8EEtjx46Fi4sLvv76a9ZR1J7GN4Gvry9MTU1x+fJl1lEIqdPbb7+NpKQkxMTEsI6i8s6dOwcfHx/w+fxWPT6Ipe+++w7r1q1DZmYm6yhqT+OLj8fj0SQXohb09PQwbdo0WtBeh6dPn2LGjBkYN24c1q9fjy1btsDU1JR1rFZhZ2eHTz75BPPmzWMdRe1p/KVOAMjNzYW3tzdyc3NhaGjIOg4hChUVFaFz585ISEigM9leEBERgXfffRcBAQFYu3Yt2rZtyzpSq3v69Cnc3d3x+++/o1evXqzjqC2NH/EBgK2tLQIDA/H333+zjkJIndq2bYsxY8Zg/fr1rKOojKqqKixcuBBvvPEGli5dit9//10rSw8AjIyM8O2339LyhmbSiuIDQCc2ELUxY8YMbNq0CRUVFayjMBcTE4OgoCDExsaq1PFBLI0ZMwZ6enrYsWMH6yhqSysudQLPdkDo1KkTkpOTYWVlxToOIXUaMmQIhgwZgkmTJrGOwoRIJML333+PlStX4vvvv8d7771HB0s/58aNGxg2bBju3r2rNfc4W5LWjPhMTEwwdOhQ/PHHH6yjEFKvWbNmYdWqVVp5GGlycjLCwsJw5swZREZG4v3336fSe0H37t0xcOBALF++nHUUtaQ1xQfQ5U6iPvr16wc9PT3897//ZR2l1dQcHxQSEoJ33nkHp06dgoODA+tYKuubb77B5s2bkZaWxjqK2tGaS50AIBaLYW9vj1OnTsHDw4N1HELqtGPHDuzevRv//PMP6yhKp87HB7G0dOlSREdHY//+/ayjqBWtGvHp6OjgnXfeoVEfUQujR49GXFwc7ty5wzqK0nAch61btyIwMBAvv/wyLl26RKXXCHPmzMGtW7dw/vx51lHUilaN+AAgLi4OgwcPRmZmJm1hRlTe0qVLkZ6ejl9//ZV1lBaXl5eHSZMmITs7G7/99ht8fHxYR1JLe/fuxfLly3Hr1i3o6OiwjqMWtO4vv7e3N9q3b0/vkIhamDx5Mg4cOICCggLWUVrU3r174efnB19fX9y4cYNKrxlGjhwJMzMzjXxzpCxaN+IDgBUrViAuLg7btm1jHYWQeoWHh8PGxgaLFi1iHaXZHj16hGnTpiE6Oho7duxAcHAw60gaISoqCoMHD8bdu3fRpk0b1nFUnlYWX15eHjw8PJCbmwsjIyPWcQipU2JiIvr27YvMzEwYGBiwjtNkx44dw6RJkzBy5EgsX76cfvda2AcffAALCwv8+OOPrKOoPK0sPgB49dVXMX78eIwdO5Z1FELqNXjwYIwYMQIffPAB6yiNVlJSgtmzZ+PMmTPYunUr+vXrxzqSRsrPz4eXlxeuXbuGLl26sI6j0rTuHl8NWtNH1MmsWbOwcuVKtVvQXnN8EPBs+zEqPeWxtrbGp59+irlz57KOovK0dsT39OlT2Nra0i74RC1wHAcfHx/89NNPePnll1nHqVdFRQUWLFiAv/76C5s3b8Zrr73GOpJWqKyshIeHBzZt2oSBAweyjqOytHbEZ2RkhGHDhtEWZkQt8Hi82lGfqrt+/Tq6deuGgoICxMbGUum1In19ffz444+YOXMmRCIR6zgqS2uLD3h2uZMOqCXqYuzYsYiOjkZCQgLrKHLVHB80dOhQfP311/jjjz/Qrl071rG0zrBhw2BlZYXNmzezjqKytPZSJ/Bsb0BHR0ccO3YM3t7erOMQUq/FixcjNzdX5f6oxcbGYsKECbCzs8OWLVtgbW3NOpJWi42NxUsvvYSkpCRYWFiwjqNytLr4AGDBggUQi8X4/vvvWUchpF4PHjyAm5sbkpOTYWlpyToORCIRfvjhB6xYsQLfffcdnaSgQiZPngwDAwOsWrWKdRSVo/XFl5CQgJdeeglZWVm03Q9RCx9++CEcHBzw5ZdfMs2RnJyMd999F4aGhti2bRudpKBiHj58CA8PD1y6dAnu7u6s46gUrb7HBwAeHh6wsbHB2bNnWUchpEFmzpyJ9evXo7KyksnrSyQSrFmzBiEhIRg7dixOnz5NpaeCLC0tsWDBAsyePZt1FJWj9cUH0Jo+ol68vLzg4+PDZEZyZmYmXnrpJezevRtXr17FJ598Qpu9q7Bp06YhJSUFJ06cYB1FpdBPLIAxY8bg77//RllZGesohDRIay9o5zgO27ZtQ2BgIAYOHIjLly/T8UFqQE9PDz/99BNmz56N6upq1nFUBhUfACsrK/Tq1QsHDx5kHYWQBhk0aBBEIlGrXKLPz8/H0KFDsWrVKpw5cwYLFiyAQCBQ+uuSlvH666/Dzs4O69evZx1FZVDx/YvW9BF1wuPxMHPmTKUvaN+7dy98fX3h6+uLmzdv0vFBaojH42HlypVYunQpCgsLWcdRCVo/q7NGRUUFbG1tERcXB1tbW9ZxCKlXRUUFHB0dcfHiRbi5ubXocxcVFWHq1KmIiorCb7/9RscHaYBp06aB4zisW7eOdRTmaMT3L0NDQ7z55pv4/fffWUchpEEMDQ0RHh7e4uu0jh8/Dm9vb1hZWSE6OppKT0MsXrwYf/31F+7cucM6CnM04nvOxYsXMXXqVMTGxtIiXKIW8vPz0bVrV6SkpDR7e7CSkhLMmTMHp06dwrZt2+gkBQ20evVqHDlyBP/8849W/42jEd9zevXqhdLSUsTExLCOQkiDWFtbY9iwYdi0aVOznuf8+fPw9fUFx3GIjY2l0tNQH3/8MXJzc3HkyBHWUZiiEd8LvvjiCzx9+hQrVqxgHYWQBomNjcWrr76KyLi7OBxXgKT8EpQIRTAzEMDd2gwjAzqhnYm+3K99/vigTZs24fXXX2/l9KS1nTx5Ep988gnu3LkDfX35PxeajorvBXfv3kXfvn2RnZ1NU7aJWojJLsbYZdtR2dYFOjo6qBRJaj9nIOCDA9DXzRJT+rjA18689nM3btzAhAkT4Ofnh3Xr1tFJClrktddeQ//+/TFnzhzWUZig4pMjODgYixcvxiuvvMI6CiF12hWRgWXHkyCsFqOuX2QeDzAQ6GDhYHe87d8RX3/9NbZs2YLVq1dj1KhRrZaXqIa7d++iV69eiI+Ph5WVFes4rY6KT45169bhypUrNMOTqLRnpZeIimpJ/Q/+l74OD/yYQ+jM3ceWLVtgY2OjxIRElc2aNQtPnz5t9v1hdUTFJ0dhYSFcXFyQlZUFMzMz1nEIkRGTXYzRWyJQUS2W+/nqolzc/3UajN1D0X7IXKnP6fI47Ps4FL52dE6bNnv8+DHc3d3xzz//wNfXl3WcVkWzOuVo3749+vTpg/3797OOQohc686nQCiSX3oAUPTPRujbdJH7ORF42HAhVVnRiJqwsLDAokWLMHPmzFbb81VVUPEpMGHCBDqxgaikwrJKXEh+CEV/q8oTLoBvYAwDB/nv4jkOOHf3IR6VsTnWiKiOjz76CIWFhVq3TzEVnwKvv/46YmJikJWVxToKIVL23cpR+DlJ5VMUX9oNi/4f1PkcPAD7ohQ/D9EOAoEAq1atwty5cyEUClnHaTVUfAro6+vjrbfewu7du1lHIURKUn6J1JKF5xVf3AkT35chMLOs8zmEIgmS8kqVEY+omQEDBsDHx6fFt75TZVR8dai53Klt17+JaisRiuT+e1VBGoSZMTALeqNBz5Nf9AQSScNnhBLN9eOPP+LHH39EXl4e6yitgmZ11oHjOLi4uGDPnj0IDAxkHYdoMY7jkJycjEuXLmFzXCUeGDnKPKbk5mEUX/wNPD3DZ19TJQQ4CXTb2cHm/Z9lHi9JvYai46vg7u4ODw8PqQ8nJyfo6Ogo+9siKuSzzz5DYWEhtm7dyjqK0lHx1WPRokV4/Pgxfv5Z9g8HIcoiEokQExODS5cu4dKlS7h8+TIMDAwQFhYGgferuFpqgSqx9K+upFoIrrKi9v8uuXEAoicFaDtoKnSM2kg91kDAx6yXXDHatz0SExORkJCAhISE2v9eUFCALl26yBSii4sLdHV1W+V/A9K6SkpK4ObmhqNHjyIgIIB1HKWi4qtHSkoKQkNDkZOTQ7/wRGkqKipw48aN2qKLiIhAp06dEBYWVvthb28P4NmsztDvziq8z1ej+NJuiIrzZNbxAYC+gI+r8/or3MOzvLwcSUlJtYVY85GdnY3OnTtLlWHXrl3h6uoKAwOD5v8PQZjasmULduzYgUuXLmn06Q1UfA0QEhKCzz//nDbwJS2muLgYV69excWLF3Hp0iXExMTA09OztuRCQ0PRvn17hV8/aWckTiUWKFzSUBceDxjk0QEbxzX+8n1FRQWSk5OlRokJCQlIS0uDvb29zAjR3d0dRkZGjQ9JmBCLxQgICMCCBQs0eis7Kr4G2LhxI86dO4c9e/awjkLUVF5eXu1o7tKlS0hNTUVQUBDCwsLQu3dv9OjRA8bGxg1+vvp2bqmLoa4O9kzqAZ9O5o3+WkWqqqqQkpIiM0K8d+8ebGxsZEaIXbt2pV2RVNSFCxcwYcIEJCUlwdDQkHUcpaDia4CioiI4OTkhMzMT5ubmrOMQFcdxHFJSUqSKrqioCL169aod0fn7+0NPT69Zr9OUvToNdflYOLgrxvVwbNZrN5RIJEJ6erpMISYlJaFt27YyI8SuXbuibdu2rZKNKPbWW2/B19cXX375JesoSkHF10AjRozAq6++ig8//JB1FKJixGIx4uLiai9bXr58GQKBQOr+nIeHB/j8ll89tCsiA0uOxkNYLQGvjud//nSG1iq9ukgkEmRmZspMqklISICxsXFtCT5fipaWlhp930mVpKenIzAwELGxsbC1tWUdp8VR8TXQoUOHsHLlSly4cIF1FMJYZWUlbt68WTuau3r1KqytrWsvW4aFhcHBwaHV/kgPHjcZlZ37IkfSBjw8W5xeo+Y8vn5ulpjS16VFL28qA8dxyM3NlRkhJiQkgM/ny4wQPTw8YGNjQ4WoBJ9//jmys7M1cutGKr4GqqqqQseOHREZGQlHR0fWcUgrKikpwdWrV2uLLioqCu7u7rWjuV69ejE70ywxMRF9+vRBeno6hJwA+6JykJRXihJhNcwMdOFuY4q3/BWfwK4uOI7DgwcP5BaiUCiUW4h2dnZKGWVri9LSUri7u2P//v3o0aMH6zgtioqvEaZOnQobGxt88cUXrKMQJSooKJC6P5ecnIzAwMDaouvZsydMTU1ZxwTwbHchNzc3LFy4kHUUZgoLC+WuRSwuLpa5XNq1a1danN8I27dvx8aNG3H16lWNehNBxdcIERERmDBhAu7evUuXVjQEx3FIT0+XKroHDx4gJCSk9rJlQEAA9PVVb8SUlpaG7t27IzU1FW3atKn/C7TMkydPZJZdJCQk4MGDB3B1dZUZIXbu3JnW6r5AIpEgODgYM2fOxDvvvMM6Touh4msEjuPg5uaGnTt3Ijg4mHUc0gQSiQR37tyRKjqO46Qmonh5eanFiCA8PByWlpZYunQp6yhqpaysrHZx/vPFmJOTU7s4//mRoqurq0q+8WktV65cwejRo5GUlNSoJTeqjIqvkZYsWYKCggKsXbuWdRTSAFVVVYiMjJSaiNK+fXuponN2dla7EXxOTg58fHyQnJxc50J30nA1i/NfHCFmZGTIXZzv5uamNYvzx4wZA1dXVyxevJh1lBZBxddI6enp6N69O3Jzc5u9Dou0vLKyMly7dq226G7evIkuXbrUXrbs1asXrK2tWcdsthkzZkAgEOCnn35iHUXjNWZxfs1oUVXuAbeUrKwsdOvWDdHR0bVb56kzKr4m6N27N+bMmYM33mjY8S9EeR4+fIjLly/XFl1iYiK6detWO5oLCQnRuPtfBQUF6Nq1K+7cuYOOHTuyjqO1RCIR0tLSZC6ZJiUloV27dnIL0cLCgnXsJvvqq6+QnJyMP/74g3WUZqPia4ItW7bg5MmT2L9/P+soWiczM7N2ofilS5dw//59hISE1BZdUFCQxm+WPH/+fJSWlmLdunWsoxA5Xlyc//yHiYmJ3KUXlpZ1HxysCsrLy+Hu7o4///wToaGhKCyrxL5bOUjKL0GJUAQzAwHcrc0wMkD1l89Q8TVBcXExHBwckJGRodbv4FSdRCJBYmKi1ESUyspKqYXiPj4+ajERpaUUFRWhS5cuiIqKgoODA+s4pBE4jkNOTo7MTNP4+HgIBAKpkaGqLs7fvXs3vv91D4Lf+xIX7hUCgNQpITUbJvR1s8SUPi7wtTNnE7QeVHxN9Pbbb2PAgAEIDw9nHUVjVFdXIyoqSuoMOnNzc6mJKF26dFGpPwStbfHixcjMzNSKw0K1BcdxKCgokCrDmnKsrKxUuDifxe/BzmsZ+PJgNKCjC0Dx66vaFnkvouJroiNHjuDbb7/FlStXWEdRW+Xl5YiIiKgtuhs3bsDZ2Vmq6Oge1v+VlpbC2dkZV65cgaurK+s4pBW8uDi/5qOkpATu7u4yhejo6Ki0KyDqsCl6Q1HxNVF1dTVsbW1x7do1dO7cmXUctVBUVCQ1ESUuLg6+vr61ly1DQkLo0nEdvv/+e0RHR2vE5ALSPMXFxUhMTJQpxQcPHsDNzU1mUk1zF+crOgar8MiPEGbEQFIthI6xBcx6jICp7yCpxyjjGKzmouJrhunTp6Nt27ZYtGgR6ygqKTs7W+r+XFZWFnr06FE7mgsODtbY875aWkVFBZycnHDq1Cl4e3uzjkNU1POL85//yM3NrV2c//xHly5dGrQ4X9HBx1UPM6Fr0RE8gS6qH2Uj//cFsBq5CPrWLrWPac7Bx8pCxdcMN2/exJgxY3Dv3j2tvu8EPLtPkZSUJFV05eXlUpct/fz8IBAIWEdVS2vWrMGZM2dw6NAh1lGIGqprcb6Dg4PMxJrnF+cXllUi9LuzUpNY5Kl+lIOC3xfAYuAkGHcNk/qcvoCPq/P6q8xsTyq+ZuA4Dh4eHvj1118REhLCOk6rEolEuH37ttREFCMjo9rLlmFhYXBzc9P6NwQtobKyEi4uLjhw4ACCgoJYxyEapKqqCvfu3ZNZi3jv3j107NgRHh4e4Hm8jAQdJ4g4+b/Lj/67HuVxZ8CJKqHXoTM6vPMt+HrSV3IMBHzMeskV4b1V47YQFV8zLV++HFlZWdi4cSPrKEpVUVGB69ev1xZdREQE7O3tpUZ0dnZ2rGNqpC1btmD//v04efIk6yhESzy/OH/trVKkiNrW+XhOIkZlbhKEWXFo0+Mt8HRkr+wM97PFylF+SkrcOFR8zVSzlc/9+/c1aiPbx48f48qVK7VFFxMTA29v79qSCw0NRbt27VjH1HgikQhubm7Yvn07wsLC6v8CQlrYxB03cTbpQYMe++jkWui2t4dZ4FCZzw1wt8Kv76rGFQu64dJM9vb28PHxwdGjRzFixAjWcZrs/v37Uvfn0tLSEBwcjLCwMCxduhTBwcEaszO7Ovnzzz/RqVMnKj3CjJlBI2pCIoHocZ6C51GdI5+o+FrA+PHjsXPnTrUpPo7jcO/ePamiKy4uRq9evRAWFobx48fD39+fziZjTCKRYPny5Vi1ahXrKESLuVubQV+QLzO5RVxeDGFmDAxduoMn0IMw4zbKEy+g/ZBPZZ7DQMCHu43qbNxNlzpbQElJCezt7ZGSkqKSR8SIxWLExMRITUTR09OTuj/XtWtXjTphWRPs378f3333Ha5fv06ThAgzimZ1ip8+wcOD36DqQTrASSBoYwXTgCEw9XtF5jloVqeGGjNmDLr17A0z35eZb9oqFApx8+ZNqTPobG1tpYqO9nlUbRzHISAgAIsWLcLQobL3SwhpTYrW8TUErePTUDHZxfjPniuIeVANfX39Vt+09cmTJ7h69Wpt0UVHR6Nr1661JderVy+12P2d/N/x48cxf/583L59m0bihDlFO7c0BO3cooGe7V+XBKFIXOe7oZbctDU/P1/q/ty9e/cQFBRUW3Q9e/aEiYlJs16DsMNxHEJCQjBz5kyMGjWKdRxCAGjWXp00uaUZGvODwHFARbUYy44nAkCDfxA4jkNaWlptyV28eBGPHj1CaGgowsLCsH79egQEBNBp8Brk3LlzKCoqwltvvcU6CiG1av5mLTueBGG1GHWNmOh0Bg2laOifv3s+Ku/fBY//bId0HdN2sJ20SeoxdQ39xWIx7ty5IzWi4/F4UvfnvLy86PKXBhswYADGjx+P9957j3UUQmTE5hRj1paTSK80gp6uLoRybu30c7PElL4uKnV583lUfE2k6GZv/u75MPbqJ7ND+fOev9lbWVmJyMhIqYkoVlZWUkXn5OREs/q0xLVr12r3f6XlJERV9erVC9M//RyPLdyQlFeKEmE1zAx04W5jirf8Vf8EdrrU2QSFZZW4kPywSTOcgGeXPf+Jz0OvAa/g9vXLcHNzQ1hYGCZOnIitW7eiQ4cOLRuYqI1ly5Zh3rx5VHpEZeXl5SE+Ph5vvDJAbXerouJrgn23cur8fPH5HSg+vwO6bW1h3ns8DBx8ZB7DA4fuo6fj+MG9MDMzU1ZUokaio6MRHR2Nffv2sY5CiEKHDx/Ga6+9pralB1DxNUlSfonCIzos+r0P3XZ24OnoojzxIh7sXwKb91dD18JG6nFi6IAz60ilR2otW7YMc+fOhYGBAesohCh04MABTJ48mXWMZqEZEk1QIhQp/Jx+Rzfw9Y3AE+jCxHsA9G27oiI1UsHzVCsrIlEzCQkJuHTpEiZNmsQ6CiEKFRUV4fr16xg0SPEcBnVAxdcEjdq0lccDFEz8VaVNWwlb33zzDWbMmEEbgROVduTIEQwcOFDtf06p+Jrg2aatsv/TSYRlqEi7BU5UBU4iRln8OVRm34Ghk7/MY1Vt01bCTmpqKk6cOIGpU6eyjkJInQ4cOIA333yTdYxmo+UMTVDXpq0P9i5CdVEOwONDt10nmIeNg6FTN5nnULVNWwk7kyZNQocOHbBkyRLWUQhRqKysDLa2tsjMzIS5uTnrOM1Ck1uaoL2JPvq4Wsqs49MxagOb91bW/wQSCXza61LpEWRnZ2Pfvn1ITk5mHYWQOh0/fhwhISFqX3oAXepssql9XWAg0GnS1+rq8HD1l6/w8ccfo7i4uGWDEbXy448/YuLEiSp5nBUhzztw4IDanDlaHyq+JvK1M8fCwe4w1G3c/4SGunx8NdQL8ZdOgOM4eHp6Yu/evaArztqnoKAAO3fuxJw5c1hHIaROQqEQJ0+e1Jgjsqj4mmFcD0csHNwVhro6AFf3RtU83rM9Omt2Kjc3N8fGjRuxd+9eLF68GK+//joyMzNbKTlRBStWrMCYMWNgY2NT/4MJYejUqVPw8/ODlZUV6ygtgoqvmcb1cMQvYzxRlX4Lejo8GLww29NAwIe+gI9BHh2wZ1IPmZ3KQ0NDER0djdDQUAQEBOCnn36CSKR4nSDRDEVFRfjll1/w2WefsY5CSL00ZTZnDZrV2QJ+/vlnREREYO2W7dgXldPkTVtTUlIwefJkFBUVYfPmzQgMVJ0Ti0nLWrRoEbKzs/Hrr7+yjkJInaqrq2FjY4Po6GjY2dmxjtMiqPiaieM4eHh4YNOmTejdu3eLPN+uXbvw6aefYvTo0ViyZAlMTWm9nyYpKSmBs7Mzrl27hi5durCOQ0idzpw5g88//xzXr19nHaXF0KXOZjp//jx0dHQQFhbWIs/H4/Ewfvx43LlzB0+ePIGnpycOHz7cIs9NVMOGDRvw8ssvU+kRtbB//36NuswJ0Iiv2d5++2307t0b06ZNU8rznzt3DuHh4fDy8sKaNWtga2urlNchrePp06dwdnbG6dOn4eXlxToOIXWSSCSwtbXFxYsXNeqNGo34miEvLw+nTp3C+PHjlfYa/fr1Q2xsLLy9veHr64s1a9ZALBbX/4VEJf3yyy/o2bMnlR5RCxEREbC0tNSo0gNoxNcsS5cuRVZWFjZv3twqr5eYmIjw8HBUVlZi8+bN8PX1bZXXJS2jsrISnTt3xqFDh2jiElELc+fOhYmJCRYtWsQ6SouiEV8TicVibN68GR9//HGrvWbXrl1x/vx5fPTRR3jppZfw2Wefoby8vNVenzTPjh074OXlRaVH1ALHcRp5fw+g4muyY8eOoWPHjujWTXYDamXi8/n48MMPERcXh9zcXHh7e+PkyZOtmoE0nkgkwrfffosvvviCdRRCGuT27dvQ0dGBt7c36ygtjoqviTZs2NCqo70XdejQAbt378aGDRswdepUjB49Gvn5+czykLr98ccfsLe3R69evVhHIaRBavbm5PF4rKO0OCq+JkhLS0NkZCTefvtt1lEwaNAgxMXFwdHRET4+Pti8eTMkkrq3TyOtSyKRYPny5Vi4cCHrKIQ0mKbt1vI8Kr4m2LRpE959910YGhqyjgIAMDIywrfffovTp09j69at6N27NxISEljHIv86cOAAzMzMMHDgQNZRCGmQpKQkPHnyBEFBQayjKAUVXyNVVlZi27ZtCA8PZx1Fho+PD65cuYKxY8eiT58++PLLLyEUClnH0mocx2Hp0qVYuHChRl4yIprpwIEDGD58OPh8zawIzfyulGjfvn3w8/NT2XUtOjo6mDJlCmJiYpCYmAhvb2+cPXuWdSytdfz4cXAch9dff511FEIaTJPO3pOH1vE1Uq9evTBnzhwMHz6cdZQGOXLkCKZNm4a+ffvip59+ogNPWxHHcejZsydmz56tEveDCWmIjIwMdO/eHffv34dAIGAdRyloxNcIsbGxyMjIwJAhQ1hHabAhQ4YgPj4e7dq1g6enJ3bs2EGH3raSs2fPori4WKPfORPNc/DgQQwdOlRjSw+g4muUDRs24KOPPlK7HwgTExOsWLECx48fx+rVqzFw4EDcu3ePdSyNt2zZMixYsAA6OjqsoxDSYJo8m7MGXepsoNLSUtjb2+POnTtqvVG0SCTCmjVrsGzZMsycOROfffYZ9PT0WMfSOFevXsU777yD5ORk6Orqso5DSIPk5+eja9euyM/Ph75+/WeIqisa8TXQrl270L9/f7UuPQAQCASYNWsWoqKicP36dfj5+eHy5cusY2mcZcuWYd68eVR6RK0cOnQIgwcP1ujSA6j4GoTjOOY7tbQ0e3t7/P3331iyZAlGjx6NSZMm4fHjx6xjaYSoqCjcvn0b7733HusohDSKNlzmBKj4GuTq1asQCoXo378/6ygtisfjYcSIEYiPj4euri48PT3x559/0uSXZlq2bBnmzp0LAwMD1lEIabCioiJcv34dr7zyCusoSkf3+Bpg3Lhx8Pf3x+zZs1lHUapr165h0qRJ6NSpE9avXw8nJyfWkdROfHw8+vfvj7S0NBgbG7OOQ0iD/fbbbzh06BAOHDjAOorS0YivHg8fPsTRo0e14rJVz549ERUVhT59+iAoKAjff/89qqurWcdSK9988w1mzpxJpUfUjqYeQSQPjfjq8f333yMhIQHbt29nHaVVpaam4uOPP0ZBQQE2b96M4OBg1pFUXmpqKoKDg5Gamoo2bdqwjkNIg5WVlaFjx47IysqCubk56zhKRyO+OkgkEmzatAlTpkxhHaXVde7cGf/9738xb948DBs2DNOmTUNJSQnrWCrt22+/xZQpU6j0iNo5ceIEQkJCtKL0ACq+Ov3zzz8wNzfX2B3K68Pj8TB27FjEx8dDKBTCw8MDBw4coMkvcmRnZ2P//v2YMWMG6yiENJqm7835IrrUWYc33ngDQ4YMwYcffsg6ikq4ePEiwsPD4erqirVr18LOzo51JJUxffp06Ovr44cffmAdhZBGEQqFsLa2RnJyMqysrFjHaRU04lMgKysLly9fxpgxY1hHURm9e/fG7du3ERAQgG7duuHnn3+GWCxmHYu5/Px87Nq1C3PmzGEdhZBGO336NHx9fbWm9AAqPoW2bNmCd955h2bnvUBfXx//+c9/cOXKFRw6dAjBwcGIjo5mHYupFStWYOzYsbC2tmYdhZBG05ZF68+jS51yVFVVwcHBAWfOnIGHhwfrOCqL4zhs374d8+fPx7hx47B48WKYmJiwjtWqHj16hC5duuD27duwt7dnHYeQRhGJRLC2tkZ0dLRW3bqgEZ8chw4dgpubG5VePXg8Ht5//33cuXMHDx8+hJeXF44dO8Y6VqtavXo1hg8fTqVH1NKFCxfg7OysVaUH0IhPrn79+mHy5MkYNWoU6yhq5fTp05g8eTL8/f3x888/w8bGhnUkpSopKYGzszOuXbuGLl26sI5DSKNNnToVdnZ2mD9/PusorYpGfC9ITExEYmKi2pywrkoGDhyIuLg4dOnSBT4+PtiwYQMkEgnrWEqzfv16DBo0iEqPqCWJRIKDBw9q3f09gEZ8MmbMmAETExMsW7aMdRS1dufOHYSHh4PjOGzevBleXl6sI7Wop0+fwsnJCWfOnNG4741oh6tXryI8PBxxcXGso7Q6GvE9p7y8HLt27cKkSZNYR1F7Xl5euHTpEt59913069cPn3/+OSoqKljHajFbtmxBaGgolR5RW9o4m7MGFd9z/vzzT4SEhMDBwYF1FI3A5/MRHh6O2NhYpKamwtvbG6dPn2Ydq9kqKyvxww8/YOHChayjENIkHMdR8ZFnNO2wWVVhY2ODPXv2YPXq1fjwww8xfvx4PHjwgHWsJtu+fTu8vb0REBDAOgohTRITEwMejwcfHx/WUZig4vvXzZs3UVRUhEGDBrGOorEGDx6M+Ph4dOjQAd7e3ti6dava7ftZXV2N7777Dl988QXrKIQ02f79+zFixAjweDzWUZig4vvXhg0bEB4eDh0dHdZRNJqxsTF+/PFHnDx5Ehs2bEC/fv1w9+5d1rEa7I8//oC9vT1CQ0NZRyGkybT5MidAxQcAePz4MQ4ePIiJEyeyjqI1unXrhoiICLz55psIDQ3F4sWLUVlZyTpWncRiMb755hsa7RG1lpSUhOLiYnTv3p11FGao+ADs2LEDgwcPhqWlJesoWkVHRwfTp09HdHQ0oqOj4evriwsXLrCOpdCBAwfQpk0bDBgwgHUUQprs4MGDGD58OPh87f3zr73f+b84jsPGjRtpUgtDdnZ2OHToEL755huMGzcOH3zwAYqKiljHksJxHJYtW4aFCxdq7X0Rohlq7u9pM60vvnPnzkFXV5fu2aiA4cOHIz4+HkZGRvD09MTu3btVZvJLzR6kr7/+OuMkhDRdZmYmMjMzERYWxjoKU1pffDVLGOhdvGowMzPDmjVrcPjwYfzwww8YNGgQUlNTmWbiOA5Lly7F559/Tj8nRK0dPHgQQ4cOhUAgYB2FKa0uvvv37+P06dMYN24c6yjkBd27d8fNmzfx0ksvITg4GN988w2qq6uZZDl79iyKi4u1/vIQUX/aPpuzhlYX3y+//IJRo0bBzMyMdRQih66uLj799FPcvHkTFy9ehL+/P65du9bqOWpGe7TUhaiz/Px8xMXFYeDAgayjMKe1xScSibBlyxaa1KIGnJyccPz4cXzxxRcYMWIEpkyZguLi4lZ57StXriAjIwNjxoxpldcjRFkOHz6MV199Ffr6+qyjMKe1xXf06FHY2dnB19eXdRTSADweD6NGjUJ8fDzEYjE8PT3x119/KX3yy7JlyzB//nzo6uoq9XUIUTa6zPl/Wnss0aBBgzBu3DiMHz+edRTSBJcvX0Z4eDicnJywbt06pWwsfuvWLbzxxhtITU2ld8lErT1+/BiOjo7Izc2FiYkJ6zjMaeWILyUlBVFRURg5ciTrKKSJevXqhejoaPTs2RMBAQFYsWIFRCJRi77GsmXLMHfuXCo9ovaOHDmC/v37U+n9SyuLb9OmTXj//fdhYGDAOgppBj09PSxcuBDXrl3DsWPHEBwcjFu3brXIc8fHx+PKlSv46KOPWuT5CGGJLnNK07pLnUKhEHZ2doiIiEDnzp1ZxyEthOM47Ny5E5999hlGjx6NJUuWwNTUtN6vKyyrxL5bOUjKL0GJUAQzAwHcrc1w9pdl6ObRBQsWLGiF9IQoT1lZGWxtbZGRkQELCwvWcVSC1q1i/OuvvxAQEEClp2F4PB4mTJiAwYMHY+7cufD09MTatWsxdOhQuY+PyS7GuvMpuJD8EABQKZLUfk5PJw+V1kNgYm2DmOxi+NqZt8a3QIhSnDhxAj179qTSe47WjfhCQkIwb948vPHGG6yjECU6d+4cwsPD4e3tjdWrV8PW1rb2c7siMrDseBKEIjHq+unn8QADgQ4WDnbHuB6Oyg9NiBKMGTMG/fr1w6RJk1hHURladY8vJiYG2dnZeO2111hHIUrWr18/xMbGwtPTE35+fli3bh3EYvG/pZeIiuq6Sw8AOA6oqBZj2fFE7IrIaJXchLQkoVCIEydO0Bv9F2jViG/y5MmwtbXFl19+yToKaUUJCQkIDw/HU4P2KAv+CJXi///IZ/30ltRjOVEVTLsNRtuXJ0v9u6GuDvZM6gGfTuatEZmQFnHs2DF89913uHjxIusoKkVriq+kpASOjo6Ij4+HjY0N6ziklUkkEry67ACSyvTBU3AOmaRKiJw142A1chEM7L2kPsfjAYM8OmDjuMDWiEtIi5g4cSJ8fX0xY8YM1lFUitZc6ty1axcGDBhApaelip5WI6PKWGHpAcDTu1egY9QG+naeMp/jOODc3Yd4VKbap8QTUkMkEuHvv//G8OHDWUdROVpRfBzH1R4/RLTTvls59T6mLO4MjL36Kzx6iAdgX1T9z0OIKrh48SKcnJxgb2/POorK0Yriu3z5Mqqrq9GvXz/WUQgjSfklUksWXiR68gCV2Xdg7D1A4WOEIgmS8kqVEY+QFkeL1hXTinV8GzZswOTJk+kQUS1WIqx7O7OyO2eh38kDuubW9TwPmzMBCWkMiUSCgwcP4ty5c6yjqCSNL74HDx7g+PHjWLduHesopJVxHIf09HRERkYiJaEU4CsutfI7Z9Gmx1sKP1/jxOH98P/52ebYzs7OcHJyqv1wdHSkbfCISrh+/TosLCzg6urKOopK0vji27p1K958803atUDDcRyH3NxcREZG4ubNm4iMjERkZCQMDQ0RGBgIa9/ByKsGquXMYRbmJEJc9ghG7r3qfA0DAR/TPxqLsPYjkJaWhvT0dMTExODQoUNIT09HdnY22rZtK1OINR+dOnWiw2xJq6DLnHXT6OUMYrEYLi4u2Lt3L4KCgljHIS2ooKCgttxqyk4ikSAoKAiBgYG1HzWzeAvLKhH63Vm59/kenVwLrroS7YfMqfM19QV8XJ3XH+1M5J/WIBaLcf/+/dpSfPHj4cOHsLOzk1uKzs7OaN++PV2OJ83GcRxcXFywf/9++Pn5sY6jkjS6+I4fP46vvvoKN2/eZB2FNENRURFu3bolNZIrLS1FYGCgVNHZ2dnVWRyTdkbiVGJBvTu2yNMS6/iEQiEyMzNlCrGmKKuqquQWYs1/pyNlSEPcvn0bI0aMQEpKCr2RUkCji2/IkCEYPnw4Jk6cyDoKaaDS0lJERUVJldyDBw/g7+9fW3BBQUFwdnZu9C91THYxRm+JQEW1uNG5WmPnlidPnsgtxPT0dGRkZMDY2FhuIdZMWdfT01NaNqI+/vOf/6CiogI//PAD6ygqS2OLLzMzEwEBAcjKyoKRkRHrOESOp0+f4vbt21KXK7OysuDr6ytVcq6uri12b+z/e3UqXtrwIkNdPhYO7sp0o2qO41BQUCBTiDUf9+/fR4cOHRTeX7SxsQG/jsX7RHN4eXlhy5Yt6NmzJ+soKktji2/hwoUoLy/HqlWrWEchACorKxEXFyc1+eTevXvw8PCQumTp4eEBXV1dpWbRxNMZqqurkZOTo/D+4pMnT+Dg4KDw/iJN/tIMd+/eRf/+/ZGdnU1vdOqgkcVXVVUFe3t7nD9/Hu7u7qzjaB2RSISEhASpy5Xx8fFwcXGRuifn4+MDfX35E0WULTanGOvPp+Dc3Yfg4dni9BoGAj44AP3cLDGlr4tGbExdXl6OjIwMhfcX+Xy+wvuLjo6OMDQ0ZP0tkAb45ptvkJubi7Vr17KOotI0svj27NmDTZs24ezZs6yjaDyJRIK7d+9KXa6MjY2FnZ2d1OVKPz8/lbzk/KisEvuicpCUV4oSYTXMDHThbmOKt/w7KZy9qWk4jkNRUZHC+4tZWVmwsLBQeH+xU6dOEAg0fmWUWggKCsJ3332H/v37s46i0jSy+Pr27YupU6di5MiRrKNoFI7jkJaWJnW5MioqCpaWllKXK/39/WFmZsY6LmkhEokE9+/fl1uK6enpePDgAWxtbeWWopOTE6ysrGh2YSvIysqCv78/8vPz6Y1IPTSu+BISEjBw4EBkZmYq/V6RJuM4Djk5OTILwo2NjWXWyrVt25Z1XMJQZWUlsrKyFN5frKiokHtvseaD3iS1jFWrViE2NhZbt25lHUXlaVzxTZ8+HW3atMGSJUtYR1ErBQUFUgVXs/axpuSCgoIQEBAAa+u697Ik5EUlJSVyCzEtLQ0ZGRkwNDRUeH/RwcGBlmk0UO/evTFv3jy89tprrKOoPI0qvvLyctjb2+P27duws7NjHUdlPXr0CLdu3ZIazZWXl8ssCO/UqRNdoiJKxXEcHjx4oPD+Ym5uLqysrBTeX+zYsSPNXsSzN65ubm7Iz8+n/WIbQKOKb8uWLTh69CgOHz7MOorKKCkpkVkQXlhYKLMg3MnJiUqOqByRSIScnByF6xcfP34Me3t7hesX27Ztq3E/14Vlldh3KwdJ+SUoEYpgZiDAk8xElMb8g327t7OOpxY0pvg4jkNAQACWL1+OV155hXUcJp4+fYro6Gipy5U5OTm1C8JrRnOurq70LplohKdPn8pdplFTlADqvL+oijONFYnJLsa68ym4kPwQAKT2neVJqqGjI8AAD2tM6eMCXztzRinVg8YU3/Xr1zF27Fjcu3dPK/6oV1ZWIjY2VupyZUpKCjw9PaUuV3p4eNAML6KVOI7D48ePFZZiZmYmzM3NFd5ftLOzU5nfHU3cdIEljSm+9957D56envj0009ZR2lx1dXVMgvCExIS4OrqKnW50svLi9mCcELUjUQiQV5ensL7iwUFBejYsaPC+4sdOnRolcuo6rrNnirTiOIrKipC586dce/ePbRv3551nGYRi8VyF4Q7ODhIXa709fVVq8s0hKibqqoqZGVlKby/WF5eDkdHR4X3F9u0adPsDPVtrF6ecAHFV/6AuOQhdIwt0O61mTCw8wLQOhurqyuNKL4VK1YgOjoaO3fuZB2lUTiOQ2pqqtTlyujoaFhZWUldrvT394epqSnruISQ55SWlsq9jFrzoaenp7AUHRwcGjT7sq6jtCrSo/HoxGpYvjEPeh1dIS4rAgAITJ+9+W+Jo7Q0ldoXn0Qigbu7O7Zv346QkBDWcRTiOA7Z2dlSJXfr1i2YmppKXa709/enBeGEqDmO4/Dw4UOF9xdzcnJgaWmp8P5ix44d8bhCpPDwZADI3zkXxj4vw9T3ZYU56js8WVupffGdPn0ac+bMwe3bt1Vq2nJ+fr7MgnA+n4+goKDa0VxAQAA6dOjAOiohpJWJRCLk5uYqvL9YVFSEjgMmgPMaDI4vuwMVJxEj68cRMA97B2Ux/4ATV8GoSw+Y95sIvu7/S85AwMesl1wR3rtza357Kk9tik/e2hV3azMcX/cfDO4fhsmTJzPL9ujRI6mCi4yMREVFhdQ9uaCgIHTs2FGlypkQopoqKiowZdcNnEsrk/t5Uekj5K57F3rWLrB86z/g8XXwcP9S6Nt7w6LPBKnHDvezxcpRfq2QWn2oxlzdOtS1dkVfJw9Ch7dgqW+DntnFrbJ25cmTJzILwh89eoSAgAAEBgZi7NixWLlyJRwdHankCCFNYmhoCJ6eEQD5xcf7d1RnGjAEApNnt0ZMg4bhydU9MsVXIqxWalZ1pNLFV9/alUoxB55AD2fvPcLV9IgWX7tSXl4usyA8NzcXfn5+CAoKwhtvvIElS5agS5cuWrF2kBDSeswMFP951jEwgY5pw2awmxnQZv0vUtnia8zaFY4DKqrFWHY8EQCaVH6VlZWIiYmRulyZlpYGLy8vBAYGYsCAAZg/fz7c3d1VZlErIURzuVubQV+Qr3Byi4n3QJTeOgpD5wBAR4DSyMMwcgmSeoyBgA93G5oR/iKVvMdX19oVUXEBHv2zHlW5SYBAF8ZuobAYOAk8vg6Ahq1dqa6uRnx8vFTJJSYmws3NTWZBOO0MTwhhobCsss5ZnZxYhKLTm1GecAE8gS6M3cNg0e998AT//5tFszrlU8niq2vtSsHer6BjZI52r0yFRFiOgj1fwMR3EMwChwKQXbtSsyD8+XtysbGxcHR0lFkQbmho2JrfJiGE1Kmuv4X1oXV8iqncNbvCskpcSH6o8P/RoicFMAt4HTyBHnRM9GDoFIDqwqzaz3MccDohH9PmzEfszauIjo6GtbV1bcGNHDkS/v7+MDExaaXviBBCmmZqXxdculeocOeWuhgIdDClr4sSUqk/lSu+fbdy6vy8WeBQlCdchL69NyTCMlSkRcI8bJzUY8RiMQpNXfDVVy/B398fFhYWyoxMCCFK4WtnjoWD3Zu4V6c7bVemgMoVX1J+icJr2gBgYOeNstv/RfaKtwFOAmOvATB07Sn1GI4vgE3XQAwY4KfktIQQolw1k/XodIaWo3Jz8EuEIoWf4zgJCvb+B0ZuIbCfsx+dZvwOibAMxee3yXkeWrtCCNEM43o4Ys+kHhjk0QH6Aj4MBNJ/ug0EfOgL+Bjk0QF7JvWg0quHyo346lq7IqkohbjkIUz9XwdPoAsdgS5MfAai+OJOWPSb+MLz0NoVQojm8Olkjo3jAvGorBL7onKQlFeKEmE1zAx04W5jirf8O9HszQZSueKra+2KjlEbCNp0QGn0cZgFvwmuqgJlcWega+Uk9Thau0II0VTtTPRp781mUrnlDPWtXakqSEPR6c2ofpAO8HVgYO+Nti9/DB1j89rH0NoVQgghiqjciK+9iT76uFoqXLui18EZ1u98q/DreTygn5sllR4hhBC5VG5yC/Bs7YqBQKdJX0trVwghhNRFJYuvZu2KoW7j4tHaFUIIIfVRuUudNWjtCiGEEGVQucktL4rNKcb68yk4d/cheACEz016MRDwweHZPb0pfV1opEcIIaReKl98NWjtCiGEkJagNsVHCCGEtASVnNxCCCGEKAsVHyGEEK1CxUcIIUSrUPERQgjRKlR8hBBCtAoVHyGEEK1CxUcIIUSrUPERQgjRKlR8hBBCtAoVHyGEEK1CxUcIIUSrUPERQgjRKlR8hBBCtAoVHyGEEK1CxUcIIUSrUPERQgjRKlR8hBBCtAoVHyGEEK1CxUcIIUSrUPERQgjRKlR8hBBCtMr/ADx050/3OYAKAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] @@ -152,42 +152,6 @@ "nx.draw(example_graph,with_labels=True)" ] }, - { - "cell_type": "code", - "execution_count": 15, - "id": "5f082e4e", - "metadata": {}, - "outputs": [ - { - "data": { - "text/plain": [ - "array([[0. , 0.33, 0. , 0.33, 0.33, 0. , 0. , 0. , 0. ],\n", - " [0.33, 0. , 0.33, 0. , 0.33, 0. , 0. , 0. , 0. ],\n", - " [0. , 0.5 , 0. , 0. , 0. , 0.5 , 0. , 0. , 0. ],\n", - " [0.33, 0. , 0. , 0. , 0.33, 0. , 0.33, 0. , 0. ],\n", - " [0.2 , 0.2 , 0. , 0.2 , 0. , 0.2 , 0. , 0.2 , 0. ],\n", - " [0. , 0. , 0.25, 0. , 0.25, 0. , 0. , 0.25, 0.25],\n", - " [0. , 0. , 0. , 0.5 , 0. , 0. , 0. , 0.5 , 0. ],\n", - " [0. , 0. , 0. , 0. , 0.25, 0.25, 0.25, 0. , 0.25],\n", - " [0. , 0. , 0. , 0. , 0. , 0.5 , 0. , 0.5 , 0. ]])" - ] - }, - "execution_count": 15, - "metadata": {}, - "output_type": "execute_result" - } - ], - "source": [ - "#h = [[0 for _ in range (9)] for _ in range(9)]\n", - "n = example_graph.number_of_nodes()\n", - "h = np.zeros((n, n))\n", - "for x in range(n): \n", - " for k in example_graph.neighbors(x):\n", - " h[x,k] = round(1/example_graph.degree(x),2)\n", - "\n", - "h" - ] - }, { "cell_type": "markdown", "id": "952c97cd", @@ -213,7 +177,7 @@ }, { "cell_type": "code", - "execution_count": 33, + "execution_count": 4, "id": "449da919", "metadata": { "deletable": false, @@ -239,7 +203,7 @@ }, { "cell_type": "code", - "execution_count": 34, + "execution_count": 5, "id": "604ce3f4", "metadata": { "deletable": false, @@ -287,7 +251,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 6, "id": "f52fcbd5", "metadata": { "deletable": false, @@ -302,14 +266,34 @@ "task": false } }, - "outputs": [], + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEICAYAAABI7RO5AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAgOUlEQVR4nO3de7RVdd3v8fdHLoIEqbA1BAz0bEFQI9uJhCnl0URL1PRJvFGnREoyMiu0HlPzjOHpkPlYBKEhpHk7XpI6HK9Jmg8aGyLdSCgq6hYUxEe84AXwe/5YE1os1l57Ttxz7wX78xpjjTXnb/5+v/mdc8D67nn7TUUEZmZmae3U1gGYmdn2xYnDzMwyceIwM7NMnDjMzCwTJw4zM8ukY1sH0Bp69eoV/fv3b+swzMy2KwsWLHg1ImpKy9tF4ujfvz/19fVtHYaZ2XZF0vPlyn2qyszMMnHiMDOzTJw4zMwsk3ZxjcPMtl/r16+nsbGRd999t61D2WF16dKFvn370qlTp1T1nTjMrKo1NjbSvXt3+vfvj6S2DmeHExGsWbOGxsZGBgwYkKqNT1WZWVV799136dmzp5NGTiTRs2fPTEd0ThxmVvWcNPKVdf86cZiZWSa+xmGWs5EjR25Tu7lz57ZoHLZtXn/9dW688Ua+9a1vMXfuXCZPnsyf/vSnVlv/zJkzOfroo9lrr70A+MY3vsH555/P4MGDM/XTkrH7iMPMrILXX3+dX//617muY8OGDU0umzlzJitWrNg8f+2112ZOGi0t1yMOSccA/wF0AK6NiCtKlg8CrgMOBn4UEZOT8oHALUVV9wEujoirJF0CnA2sTpZdFBFz8twOsw/DRw7bt0mTJvHMM88wdOhQOnXqRLdu3Tj55JNpaGjgU5/6FDfccAOSWLBgAeeffz5vvfUWvXr1YubMmfTu3ZtFixYxfvx41q1bx7777suMGTPYbbfdGDlyJJ/5zGd45JFHOP744xk5cuRW7R955BHq6+s5/fTT6dq1K/PmzWPUqFFMnjyZuro67r77bi666CI2btxIr169eOCBB/jb3/7GxIkTeeedd+jatSvXXXcdAwcObNF9klvikNQBmAIcBTQC8yXNjogni6q9BpwHnFDcNiKWAkOL+nkJuLOoyi82JRmzanHH0pWttq6TBvZutXW1d1dccQUNDQ0sWrSIuXPnMnr0aBYvXsxee+3FiBEjeOSRRxg2bBjf/va3ueuuu6ipqeGWW27hRz/6ETNmzOCss87il7/8JUcccQQXX3wxl156KVdddRVQOJr5y1/+wvr16zniiCPKtv/Vr361OVEUW716NWeffTYPPfQQAwYM4LXXXgNg0KBBPPTQQ3Ts2JH777+fiy66iNtvv71F90meRxyHAMsi4lkASTcDo4HNiSMiVgGrJB1XoZ8jgWciouxgW2bWfrXF9aNDDjmEvn37AjB06FCWL1/OrrvuSkNDA0cddRQAGzdupHfv3qxdu5bXX3+dI444AoCxY8dyyimnbO7rK1/5CgBLly4t276SRx99lMMPP3zzsxe77747AGvXrmXs2LE8/fTTSGL9+vXbvK1NyTNx9AFeLJpvBIZtQz+nAjeVlE2QdBZQD3wvIv6rtJGkccA4gL333nsbVmtmtrWdd95583SHDh3YsGEDEcGQIUOYN2/eFnXXrl1bsa9u3boBNNm+kogoexvtv//7v/O5z32OO++8k+XLl29zcq0kz8RR7sbgyNSB1Bk4HriwqHgq8NOkr58CPwf+x1YripgOTAeoq6vLtF4z2z60xvWj7t278+abb1asM3DgQFavXs28efMYPnw469ev56mnnmLIkCHstttuPPzww3z2s5/l+uuv33z0kbZ9U+sfPnw45557Ls8999zmU1W77747a9eupU+fPkDhwnoe8kwcjUC/ovm+wIom6jZlFLAwIl7ZVFA8LekaoPXuizOzdqdnz56MGDGCAw44gK5du7LnnntuVadz587cdtttnHfeeaxdu5YNGzYwceJEhgwZwqxZszZfHN9nn3247rrrMrX/6le/yvjx4zdfHN+kpqaG6dOnc9JJJ/HBBx+wxx57cN999/GDH/yAsWPHcuWVV/L5z38+l32iiHz+GJfUEXiKwjWKl4D5wGkRsbhM3UuAt0oveCfXRe6JiOuKynpHxMpk+rvAsIg4tVIsdXV14Rc5Wd58cTwfS5YsYf/992/rMHZ45fazpAURUVdaN7cjjojYIGkCcA+F23FnRMRiSeOT5dMkfYzCdYoewAeSJgKDI+INSbtQuCPrnJKufyZpKIVTVcvLLDczsxzl+hxH8nzFnJKyaUXTL1M4hVWu7TqgZ5nyM1s4TDMzy8BPjpuZWSZOHGZmlokTh5mZZeLEYWZmmXhYdTPbrrT0bc/VdmvzyJEjy45N1ZS2GOrdRxxmZpaJE4eZWTPefvttjjvuOD7xiU9wwAEHcMstt3DZZZfx6U9/mgMOOIBx48ax6WHqkSNH8t3vfpfDDz+c/fffn/nz53PSSSdRW1vLj3/8YwCWL1/OoEGDGDt2LAcddBAnn3wy69at22q99957L8OHD+fggw/mlFNO4a233gLg7rvvZtCgQRx22GHccccdrbcjEk4cZmbNuPvuu9lrr734xz/+QUNDA8cccwwTJkxg/vz5NDQ08M4772xxqqhz58489NBDjB8/ntGjRzNlyhQaGhqYOXMma9asAQoj4o4bN47HH3+cHj16bPWyqFdffZXLL7+c+++/n4ULF1JXV8eVV17Ju+++y9lnn80f//hHHn74YV5++eVW3RfgxGFm1qwDDzyQ+++/nx/+8Ic8/PDDfPSjH+XBBx9k2LBhHHjggfz5z39m8eJ/jaZ0/PHHb243ZMgQevfuzc4778w+++zDiy8WBg3v168fI0aMAOCMM87gr3/96xbrfPTRR3nyyScZMWIEQ4cOZdasWTz//PP885//ZMCAAdTW1iKJM844o5X2wr/44riZWTP2228/FixYwJw5c7jwwgs5+uijmTJlCvX19fTr149LLrmEd999d3P9TUOv77TTTlsMw77TTjttfk1s6ZDopfMRwVFHHcVNN235VolFixaVHU69NfmIw8ysGStWrGCXXXbhjDPO4IILLmDhwoUA9OrVi7feeovbbrstc58vvPDC5tFub7rpJg477LAtlh966KE88sgjLFu2DIB169bx1FNPMWjQIJ577jmeeeaZzW1bm484zGy70ha3zz7xxBN8//vfZ6eddqJTp05MnTqVP/zhDxx44IH079+fT3/605n73H///Zk1axbnnHMOtbW1fPOb39xieU1NDTNnzmTMmDG89957AFx++eXst99+TJ8+neOOO45evXpx2GGH0dDQ0CLbmVZuw6pXEw+rbq3Bw6rnY0ccVn358uV88YtfbPUf/EqyDKvuU1VmZpaJE4eZWSvr379/VR1tZOXEYWZVrz2cUm9LWfevE4eZVbUuXbqwZs0aJ4+cRARr1qyhS5cuqdv4riozq2p9+/alsbGR1atXt3UoO6wuXbrQt2/Zl7GW5cRhZlWtU6dODBgwoK3DsCI+VWVmZpnkmjgkHSNpqaRlkiaVWT5I0jxJ70m6oGTZcklPSFokqb6ofHdJ90l6OvneLc9tMDOzLeWWOCR1AKYAo4DBwBhJg0uqvQacB0xuopvPRcTQkgdQJgEPREQt8EAyb2ZmrSTPI45DgGUR8WxEvA/cDIwurhARqyJiPrA+Q7+jgVnJ9CzghBaI1czMUsozcfQBXiyab0zK0grgXkkLJI0rKt8zIlYCJN97lGssaZykekn1vhvDzKzl5Jk4yo37m+VG7BERcTCFU13nSjo8y8ojYnpE1EVEXU1NTZamZmZWQZ6JoxHoVzTfF1iRtnFErEi+VwF3Ujj1BfCKpN4AyfeqFonWzMxSyTNxzAdqJQ2Q1Bk4FZidpqGkbpK6b5oGjgY2DewyGxibTI8F7mrRqM3MrKLcHgCMiA2SJgD3AB2AGRGxWNL4ZPk0SR8D6oEewAeSJlK4A6sXcGfylquOwI0RcXfS9RXArZK+DrwAnJLXNpiZ2dZyfXI8IuYAc0rKphVNv0zhFFapN4BPNNHnGuDIFgzTzMwy8JPjZmaWiceqMtvBVHoT4cVnfnmb+rzs+tubXLY9vI3Q+6Rl+YjDzMwy8RGHWTtS6a/k9sr7JDsfcZiZWSZOHGZmlokTh5mZZeLEYWZmmThxmJlZJk4cZmaWiROHmZll4sRhZmaZOHGYmVkmThxmZpaJE4eZmWXixGFmZpk4cZiZWSZOHGZmlokTh5mZZeLEYWZmmeSaOCQdI2mppGWSJpVZPkjSPEnvSbqgqLyfpAclLZG0WNJ3ipZdIuklSYuSz7F5boOZmW0ptzcASuoATAGOAhqB+ZJmR8STRdVeA84DTihpvgH4XkQslNQdWCDpvqK2v4iIyXnFbmZmTcvziOMQYFlEPBsR7wM3A6OLK0TEqoiYD6wvKV8ZEQuT6TeBJUCfHGM1M7OUmk0ckuolnStpt4x99wFeLJpvZBt+/CX1Bz4JPFZUPEHS45JmNBWXpHFJ7PWrV6/OulozM2tCmiOOU4G9KJxqulnSFyQpRbtydSJLcJI+AtwOTIyIN5LiqcC+wFBgJfDzcm0jYnpE1EVEXU1NTZbVmplZBc0mjohYFhE/AvYDbgRmAC9IulTS7hWaNgL9iub7AivSBiapE4Wk8fuIuKMonlciYmNEfABcQ+GUmJmZtZJU1zgkHUThL/v/TeHH/GTgDeDPFZrNB2olDZDUmcKRy+yU6xPwW2BJRFxZsqx30eyJQEOaPs3MrGU0e1eVpAXA6xR+yCdFxHvJosckjWiqXURskDQBuAfoAMyIiMWSxifLp0n6GFAP9AA+kDQRGAwcBJwJPCFpUdLlRRExB/iZpKEUTnstB87JssFmZvbhpLkd95SIeLa4QNKAiHguIk6q1DD5oZ9TUjataPplCqewSv2V8tdIiIgzU8RsZmY5SXOq6raUZWZm1g40ecQhaRAwBPiopOIjix5Al7wDMzOz6lTpVNVA4IvArsCXisrfBM7OMSYzM6tiTSaOiLgLuEvS8IiY14oxmZlZFat0quoHEfEz4DRJY0qXR8R5uUZmZmZVqdKpqiXJd31rBGJmZtuHSqeq/ph8z9pUJmkn4CNFw3+YmVk7k2aQwxsl9ZDUDXgSWCrp+/mHZmZm1SjNcxyDkyOMEyg8zLc3hae6zcysHUqTODolAw6eANwVEevJOMqtmZntONIkjt9QGBOqG/CQpI9TGODQzMzaoWbHqoqIq4Gri4qel/S5/EIyM7NqlmZ03J2BLwP9S+pfllNMZmZWxdKMjnsXsBZYALzXTF0zM9vBpUkcfSPimNwj2Q6NHDlym9rNnTu3ReOA6oqlWnifmOUjTeL4T0kHRsQTuUdThe5YurLJZa+ue7/F+zxpYO8ml7VmLJXiMLP2LU3iOAz4qqTnKJyqEhARcVCukW0HLrv+9rYOYbNqiqU1VUqm5/3mphbv0wnVLF3iGJV7FGZmtt1o9jmOiHge6Ad8Pplel6admZntmNKMVfUT4IfAhUlRJ+CGPIMyM7PqlebI4UTgeOBtgIhYAXRP07mkYyQtlbRM0qQyywdJmifpPUkXpGkraXdJ90l6OvneLU0sZmbWMtIkjvcjIkjGp0pGyW2WpA7AFArXSAYDYyQNLqn2GnAeMDlD20nAAxFRCzyQzJuZWStJkzhulfQbYFdJZwP3A9ekaHcIsCwino2I94GbgdHFFSJiVUTMB9ZnaDsa2PSOkFkUBl80M7NWkmasqsmSjqIwsOFA4OKIuC9F332AF4vmG4FhKeOq1HbPiFiZxLZS0h7lOpA0DhgHsPfee6dcrZmZNSfN7bgkiSJNsiimcl21QttC5YjpwHSAuro6DwNvZtZCmkwckt6kwo91RPRopu9GCrfxbtIXWJEyrkptX5HUOzna6A2sStmnmZm1gErvHO8OIOky4GXgegpHAqeT7q6q+UCtpAHAS8CpwGkp46rUdjYwFrgi+b4rZZ9mZtYC0pyq+kJEFF+bmCrpMeBnlRpFxAZJE4B7gA7AjIhYLGl8snyapI8B9UAP4ANJE0leVVuubdL1FRQu2H8deAE4Je3GmpnZh5cmcWyUdDqFO5sCGANsTNN5RMyh8J7y4rJpRdMvUzgNlaptUr4GODLN+s3MrOWluR33NODfgFeSzymkP+VkZmY7mDS34y6n5PkLMzNrvzxYoZmZZeLEYWZmmThxmJlZJs1e45B0fpnitcCCiFjU4hGZmVlVS3PEUQeMpzB+VB8K4z+NBK6R9IP8QjMzs2qU5jmOnsDBEfEWbH6x023A4cACmnkQ0MzMdixpjjj2Bt4vml8PfDwi3gHeyyUqMzOrWmmOOG4EHpW0aUyoLwE3JS90ejK3yMzMrCqleQDwp5L+HzCCwiCH4yOiPll8ep7BmZlZ9Un1Pg7g7xSGNe8IIGnviHght6jMzKxqpbkd99vATyiMU7WRwlFHAAflG5qZmVWjNEcc3wEGJqPSmplZO5fmrqoXKTzwZ2ZmluqI41lgrqT/S9HttxFxZW5RmZlZ1UqTOF5IPp2Tj5mZtWNpbse9tDUCMTOz7UOTiUPSVRExUdIfKdxFtYWIOD7XyMzMrCpVOuK4Pvme3BqBmJnZ9qHJu6oiYkEyOTQi/lL8AYam6VzSMZKWSlomaVKZ5ZJ0dbL8cUkHJ+UDJS0q+rwhaWKy7BJJLxUtOzbrRpuZ2bZLczvu2DJlX22ukaQOwBRgFDAYGCNpcEm1UUBt8hkHTAWIiKURMTQihgKfAtYBdxa1+8Wm5RExJ8U2mJlZC6l0jWMMcBowQNLsokXdgTQPAx4CLIuIZ5P+bgZGs+XAiKOB30VEUBhIcVdJvSNiZVGdI4FnIuL5VFtkZma5qnSN4z+BlUAv4OdF5W8Cj6fouw+Fhwc3aQSGpajTJ1nvJqcCN5W0myDpLKAe+F5E/FfpyiWNo3AUw957750iXDMzS6PSNY7nI2JuRAwvucaxMCI2pOhb5brNUkdSZ+B44P8ULZ8K7EvhOstKtkxqxfFPj4i6iKirqalJEa6ZmaVR6VTVXyPiMElvsuUPvoCIiB7N9N0I9Cua70thhN0sdUYBCyPilU0FxdOSrgH+1EwcZmZV446lK5uv1IJOGti7xftsMnFExGHJd/dt7Hs+UCtpAPAShVNOp5XUmU3htNPNFE5jrS25vjGGktNUJddATgQatjE+y8HIkSO3qd3cuXNbNA6rfv63srWLz/zyNrW77PrbWziSytIMq74v0BgR70kaSWE49d9FxOuV2kXEBkkTgHuADsCMiFgsaXyyfBowBzgWWEbhzqmvFa13F+Ao4JySrn8maSiFo6DlZZabWZWo9Nf1q+veb3LZtvSZx1/WVl6asapuB+ok/TfgtxSOEm6k8INfUXKr7JySsmlF0wGc20TbdUDPMuVnpojZclTpx+C835Tex/Dh+/QPwo6ptf9K3h5sL/skzXMcHyQXw08EroqI7wL+n2xm1k6lSRzrk2c6xvKvC9Gd8gvJzMyqWZrE8TVgOPA/I+K55GL3DfmGZWZm1SrNsOpPAucVzT8HXJFnUGZmVr0qPcdxa0T8m6QnKD+s+kG5RmZmZlWp0hHHd5LvL7ZGIGZmtn2o9ADgpvsjTwJujYiXWickMzOrZmkujvcA7pX0sKRzJe2Zd1BmZla9mk0cEXFpRAyh8KDeXsBfJN2fe2RmZlaV0hxxbLIKeJnCuzj2yCccMzOrds0mDknflDQXeIDCuznO9h1VZmbtV5qxqj4OTIyIRTnHYmZm24E0DwBOao1AzMxs+5DlGoeZmZkTh5mZZePEYWZmmThxmJlZJk4cZmaWiROHmZll4sRhZmaZ5Jo4JB0jaamkZZK2eh5EBVcnyx+XdHDRsuWSnpC0SFJ9Ufnuku6T9HTyvVue22BmZlvKLXFI6gBMAUYBg4ExkgaXVBsF1CafccDUkuWfi4ihEVFXVDYJeCAiaikMg+IHFM3MWlGeRxyHAMsi4tmIeB+4GRhdUmc08LsoeBTYVVLvZvodDcxKpmcBJ7RgzGZm1ow8E0cf4MWi+cakLG2doPAekAWSxhXV2XPTS6aS77Ij9UoaJ6leUv3q1as/xGaYmVmxPBOHypSVvru8Up0REXEwhdNZ50o6PMvKI2J6RNRFRF1NTU2WpmZmVkGeiaMR6Fc03xdYkbZORGz6XgXcSeHUF8Arm05nJd+rWjxyMzNrUp6JYz5QK2mApM7AqcDskjqzgbOSu6sOBdZGxEpJ3SR1B5DUDTgaaChqMzaZHgvcleM2mJlZiTTv49gmEbFB0gTgHqADMCMiFksanyyfBswBjgWWAeuAryXN9wTulLQpxhsj4u5k2RXArZK+DrwAnJLXNpiZ2dZySxwAETGHQnIoLptWNB0U3mVe2u5Z4BNN9LkGOLJlIzUzs7T85LiZmWXixGFmZpk4cZiZWSZOHGZmlokTh5mZZeLEYWZmmThxmJlZJk4cZmaWiROHmZll4sRhZmaZOHGYmVkmThxmZpaJE4eZmWXixGFmZpk4cZiZWSZOHGZmlokTh5mZZeLEYWZmmThxmJlZJk4cZmaWSa6JQ9IxkpZKWiZpUpnlknR1svxxSQcn5f0kPShpiaTFkr5T1OYSSS9JWpR8js1zG8zMbEsd8+pYUgdgCnAU0AjMlzQ7Ip4sqjYKqE0+w4CpyfcG4HsRsVBSd2CBpPuK2v4iIibnFbuZmTUtzyOOQ4BlEfFsRLwP3AyMLqkzGvhdFDwK7Cqpd0SsjIiFABHxJrAE6JNjrGZmllKeiaMP8GLRfCNb//g3W0dSf+CTwGNFxROSU1szJO1WbuWSxkmql1S/evXqbdwEMzMrlWfiUJmyyFJH0keA24GJEfFGUjwV2BcYCqwEfl5u5RExPSLqIqKupqYmY+hmZtaUPBNHI9CvaL4vsCJtHUmdKCSN30fEHZsqRMQrEbExIj4ArqFwSszMzFpJnoljPlAraYCkzsCpwOySOrOBs5K7qw4F1kbESkkCfgssiYgrixtI6l00eyLQkN8mmJlZqdzuqoqIDZImAPcAHYAZEbFY0vhk+TRgDnAssAxYB3wtaT4COBN4QtKipOyiiJgD/EzSUAqntJYD5+S1DWZmtrXcEgdA8kM/p6RsWtF0AOeWafdXyl//ICLObOEwzcwsAz85bmZmmThxmJlZJk4cZmaWiROHmZll4sRhZmaZOHGYmVkmThxmZpaJE4eZmWXixGFmZpk4cZiZWSZOHGZmlokTh5mZZeLEYWZmmThxmJlZJk4cZmaWiROHmZll4sRhZmaZOHGYmVkmThxmZpaJE4eZmWWSa+KQdIykpZKWSZpUZrkkXZ0sf1zSwc21lbS7pPskPZ1875bnNpiZ2ZZySxySOgBTgFHAYGCMpMEl1UYBtclnHDA1RdtJwAMRUQs8kMybmVkryfOI4xBgWUQ8GxHvAzcDo0vqjAZ+FwWPArtK6t1M29HArGR6FnBCjttgZmYlOubYdx/gxaL5RmBYijp9mmm7Z0SsBIiIlZL2KLdySeMoHMUAvCVp6bZsxA6iF/BqWwdRZbxPtuZ9srX2vk8+Xq4wz8ShMmWRsk6athVFxHRgepY2OypJ9RFR19ZxVBPvk615n2zN+6S8PE9VNQL9iub7AitS1qnU9pXkdBbJ96oWjNnMzJqRZ+KYD9RKGiCpM3AqMLukzmzgrOTuqkOBtclpqEptZwNjk+mxwF05boOZmZXI7VRVRGyQNAG4B+gAzIiIxZLGJ8unAXOAY4FlwDrga5XaJl1fAdwq6evAC8ApeW3DDsSn7LbmfbI175OteZ+UoYhMlw7MzKyd85PjZmaWiROHmZll4sSxg5LUT9KDkpZIWizpO20dU7WQ1EHS3yX9qa1jqQaSdpV0m6R/Jv9ehrd1TNVA0neT/zsNkm6S1KWtY6oWThw7rg3A9yJif+BQ4NwyQ760V98BlrR1EFXkP4C7I2IQ8Am8b5DUBzgPqIuIAyjcpHNq20ZVPZw4dlARsTIiFibTb1L4MejTtlG1PUl9geOAa9s6lmogqQdwOPBbgIh4PyJeb9OgqkdHoKukjsAubP0cWrvlxNEOSOoPfBJ4rI1DqQZXAT8APmjjOKrFPsBq4Lrk9N21krq1dVBtLSJeAiZTuOV/JYVnzO5t26iqhxPHDk7SR4DbgYkR8UZbx9OWJH0RWBURC9o6lirSETgYmBoRnwTexiNOk7yuYTQwANgL6CbpjLaNqno4cezAJHWikDR+HxF3tHU8VWAEcLyk5RRGXP68pBvaNqQ21wg0RsSmo9HbKCSS9u6/A89FxOqIWA/cAXymjWOqGk4cOyhJonDeeklEXNnW8VSDiLgwIvpGRH8KFzr/HBHt+q/IiHgZeFHSwKToSODJNgypWrwAHCppl+T/0pH4poHN8hwd19rWCOBM4AlJi5KyiyJiTtuFZFXq28Dvk3HhniUZ+qc9i4jHJN0GLKRwh+Lf8fAjm3nIETMzy8SnqszMLBMnDjMzy8SJw8zMMnHiMDOzTJw4zMwsEycOMzPLxInDzMwyceIwayOSLpF0QVvHYZaVE4eZmWXixGG2jST1T96Yd03yprh7JXVNlp2fvDmuQdLEojY/krRU0v3AwKLyMyT9TdIiSb+R1KGJdT4o6ahk+nJJV+e7lWZb81hVZh9OLTAmIs6WdCvwZUlLKIz3NAwQ8Jikv1D4Q+1UCu9G6UhhHKQFkvYHvgKMiIj1kn4NnA78rsz6fgJcJmmPpJ/j8908s605cZh9OM9FxKJkegHQH+gJ3BkRbwNIugP4LIXEcWdErEvKZyftjgQ+BcwvDMRKV2BVuZVFxEPJaK3nAyMjYmMO22RWkROH2YfzXtH0Rgo/+qpQv9yoogJmRcSFza1M0oFAb+DV5JXAZq3O1zjMWt5DwAnJuxy6AScCDyflJ0rqKqk78KWk/gPAycnpJyTtLunjpZ1K6g38nsKb6d6W9IVW2BazrfiIw6yFRcRCSTOBvyVF10bE3wEk3QIsAp6nkEyIiCcl/Ri4V9JOwHrg3KQOSbtdKLyF7nsRsUTST4H/BdzTKhtlVsTv4zAzs0x8qsrMzDJx4jAzs0ycOMzMLBMnDjMzy8SJw8zMMnHiMDOzTJw4zMwsk/8PA1ayPmeD2FoAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "def chain_Q_histogram(graph,start,k):\n", " '''Produce a histogram (a Numpy array of length equal to the number of \n", " nodes of graph) of the states visited (excluding initial state) by the \n", " Q Markov chain in the first k steps when started at start.'''\n", - " # YOUR CODE HERE\n", - " raise NotImplementedError()\n", + " n = graph.number_of_nodes()\n", + " number_of_visits = np.zeros(n)\n", + " x = start\n", + " \n", + " for _ in range(k):\n", + " x = sample_proposal(graph, x)\n", + " number_of_visits[x] += 1\n", + " \n", + " return number_of_visits\n", "\n", "def transition_matrix_Q(graph):\n", " '''Construct transition matrix Q from graph as two-dimensional Numpy array.'''\n", @@ -317,17 +301,26 @@ " Q = np.zeros((n, n))\n", " for x in range(n): \n", " for k in example_graph.neighbors(x):\n", - " Q[x,k] = round(1/example_graph.degree(x),2)\n", + " Q[x,k] = 1/example_graph.degree(x)\n", " return Q\n", "\n", "# Compare histogram and stationary distribution in a plot\n", - "# YOUR CODE HERE\n", - "raise NotImplementedError()" + "x_start = 1\n", + "k = 100000\n", + "\n", + "plt.figure()\n", + "x_list = [i + 1 for i in range(example_graph.number_of_nodes())]\n", + "plt.bar(x_list, chain_Q_histogram(example_graph, x_start, k)/k, color=\"lightblue\", label=\"sampled\")\n", + "plt.scatter(x_list, stationary_distributions(transition_matrix_Q(example_graph)), s=200, marker=\"_\", color=\"black\", label=\"theoretical\")\n", + "plt.ylabel(\"visiting density\")\n", + "plt.xlabel(\"node $x$\")\n", + "plt.legend()\n", + "plt.show()" ] }, { "cell_type": "code", - "execution_count": null, + "execution_count": 7, "id": "812cd2ef", "metadata": { "deletable": false,