{ "cells": [ { "cell_type": "markdown", "id": "ea6f6697", "metadata": {}, "source": [ "# Core Distribution Samplers" ] }, { "cell_type": "code", "execution_count": 1, "id": "c4dd2a04", "metadata": { "execution": { "iopub.execute_input": "2026-06-13T10:26:25.028405Z", "iopub.status.busy": "2026-06-13T10:26:25.028031Z", "iopub.status.idle": "2026-06-13T10:26:26.524725Z", "shell.execute_reply": "2026-06-13T10:26:26.524266Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "lal.MSUN_SI != Msun\n" ] } ], "source": [ "import warnings\n", "warnings.filterwarnings('ignore', 'Wswiglal-redir-stdio')\n", "\n", "from gwGenealogy.utils import (\n", " sample_uniform_1d,\n", " sample_loguniform_1d,\n", " sample_gaussian_1d,\n", " sample_lognormal_1d,\n", " sample_powerlaw_1d,\n", " sample_maxwellian_1d,\n", " sample_beta_1d,\n", ")\n", "from gwGenealogy.utils import set_rcparams\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "set_rcparams()" ] }, { "cell_type": "markdown", "id": "ffdb930a", "metadata": {}, "source": [ "## Draw samples from all seven distributions" ] }, { "cell_type": "code", "execution_count": 2, "id": "52285cd2", "metadata": { "execution": { "iopub.execute_input": "2026-06-13T10:26:26.526130Z", "iopub.status.busy": "2026-06-13T10:26:26.525995Z", "iopub.status.idle": "2026-06-13T10:26:26.529076Z", "shell.execute_reply": "2026-06-13T10:26:26.528703Z" } }, "outputs": [], "source": [ "uniform = sample_uniform_1d(5000, low=0, high=10, seed=0)\n", "logunif = sample_loguniform_1d(5000, low=1, high=1000, seed=1)\n", "gauss = sample_gaussian_1d(5000, mean=30, std=5, seed=2)\n", "lognorm = sample_lognormal_1d(5000, mean=1.0, sigma=0.5, seed=3)\n", "powerlaw = sample_powerlaw_1d(5000, beta=-2.35, xmin=5, xmax=100, seed=4)\n", "maxwell = sample_maxwellian_1d(5000, sigma=265, seed=5)\n", "beta_samp = sample_beta_1d(5000, a=1.4, b=3.6, seed=6)" ] }, { "cell_type": "markdown", "id": "a0c40fca", "metadata": {}, "source": [ "## Histograms" ] }, { "cell_type": "code", "execution_count": 3, "id": "a3f11c6b", "metadata": { "execution": { "iopub.execute_input": "2026-06-13T10:26:26.530121Z", "iopub.status.busy": "2026-06-13T10:26:26.530059Z", "iopub.status.idle": "2026-06-13T10:26:26.873415Z", "shell.execute_reply": "2026-06-13T10:26:26.873075Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAABhoAAAL5CAYAAABGsP6IAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjksIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvJkbTWQAAAAlwSFlzAAAPYQAAD2EBqD+naQAAvZxJREFUeJzs/VuMHWeeGHhGJhMtGihSFLOkLfS0t1dM8mUMs6TkBbItA4Mi6XqqMVDFyz6MFwvMkNleQ94BbCSLetiXBcTKhN4Ew51UGxgssEBTZDfWkIHFFJPtedB4tENmVpkPheGQSWnQhV5ZKl4q2W6LXJK5+Ic6UidPxrnHucSJ3w84PMk4cb6IEyfO90V8/+8ysbGxsZEAAAAAAAB0YbKbNwEAAAAAAAg0AAAAAAAAPdGjAQAAAAAA6JpAA5TIvXv3knGwuro67F0AYIzKlX5QVgGMT3khTwcYrLKXG/22Oqb1YgINUKJMaGlpKRkX58+fH/YuAFTauJUr/aCsAhif8kKeDjAY41Ju9Nv5MawXE2igrx49epScOnUqOXToUDIxMZE+4v/1P6arV69urvPKK68kJ06c6Dr6GducmZlJFhcXG2Z4sQ/xei/bGaT4THHMFhYWtr0WnyNei+e5ubn0WPZrH+K4tUq/nf2ZnZ1Npqenk0uXLvVlX4H+5uORd8bve3l5eeCHWrnS/3Kl9jiXoUzp17rKKiiXuMaP33OUV5G/xCP+n13rR5nVrzytn1rd2wyrvMiWZ/lpXBs028d+3bPEdiPdrGVofN/x/9hGLXk6lM+o3ofEPmT7k1efFPtWVP3WKInPEJ8rPlM8D+s+o991YL2mf6Lq5dIGDMDDhw834nTbt29f0/XOnTu3ceXKlZ62tba2lm7r5MmTDV+L59inPXv2bCwtLW2MuuPHj6f7nHe8FhYWtq1b5GeK4xjbiUccu2bfT6f7Mzs7m/u5gPLm44OiXOlPudJJnj9KZUq/1g3KKhh98/PzafmUl6fEa/H7juv+ovO0QWh2bzOs8iLK4MhL65fFMc7bz37es0QeHcen9hFpN1vf/QeUy6jdh4TIZ2KfIk9pJMqfMpY77Xz2Zp+7n/cZ/a4DKyL92YqXSwINDO5ka5EJh/jx9vMHFhnGKBVO7bh27VpuprSyspIe03aXF3WT06gw6GZ/Iq1h3TQB/cnHB0m5Umy50kmeP0plSr/WzSirYLRFuRSPqIhqpF/B03HXqLyIe7Y4nvF6rbiuzxp1DeqeJbYZ93jxHJV69ftUT54O5TRq9yGR38Qj9qtRJXTkN+NUgVxbpjb6Lvp5n9Hv8qSo9E9WvFwydBKVsra2luzbty8pkxjXLm/ctlgeXa3qZcsG3TW8m/05efJk2q0wus4BlNE4lSujpJMypV/rZpRVMLpiWIMYmuDKlSvJnj17Gq4Xv/1mr9NZeXH48OHc47l37970ufbavt/3LLHN+fn59ByIYTqOHz/edH15OlCUyHPiPqB2mL6q6+d9Rr/Lk6LS31vxckmgAUZcZDh5GVMsb1S5FRf+165dG8De9b4/8dnGbkw6oHAu3vtfroySTsqUfq1bS1kFoycb8zh+n+0EfE+fPj2Q/apCeRGVLg8fPtz2WpbH1lbUjNo9S5CnA0Xdh0Rlcqgff7+q+nmf0e/yZJjl1fExqhebGvYOQDPRQuns2bNpxh4/vA8//DD98WU/8sgEaieYifUig79582b62srKymY6Fy9eTJ8fPHiQTjxTG4HO1E5gdv/+/XTitXPnzm3JeOKGJrYTy8+cOZMui32J7UZazfY3IpsRJY1IZZaJ3LhxY9vnyETEtFFGl20jT2wnjsEgdbs/MVFOFM4R8QXGU6u8tb7S6MiRI2lPgVgv8s/IQy5fvpzmudHioxfKlcblyijppEzp17q1lFUweuKaOvt9tiOu1eM+oF7kC5FWTMgYZVSovy6vLTvq7zHyltemG2VZ1koxyrbYj6wivtU6je5tut33du6n2rkPyRPrx/G9fv36UO5Z4nNGWtHbIq9Fai15OlRHv+9DIr+JuozYTuSteWl3s1+DrntqpzxppdNyo1P9Lk+KTn+1quXSsMduojp6GUs7xniLR/2kLDHhWP2ybP28bTUbSy6W148RF+Op5Y2VFvM81E4SE/+vXa/Z/kaa7X6ObNy/RpMh1U/CVvtZip6Lotk4er3sT7/mlABGY2zUdvPWbBLJ2jEsIx/N8pXIK+rH31auFFeujNIcDZ2UKf1at56yCkZPNhFnr/lUfZ4QeWSz+4W81/LuPRqtF/lJJ+s0Sr/Tfe/0fqqd8iLy8mwc6ki7vpwexD1L/RjY8dxqYk15OpTPqN2H1OePkZfFPtau12iOhlGse+qkPGlUr9bP+4x+lydFpn+y4uWSoZMohYjKRvS2PoIckcGILtfrNIoaEeGIUNanHxHc2G79WGyR/kcffbQZcY7WRVmXuVb7G5Hrdj9HRFQjojzOsu9qXMajA7rLW6OlTeQDta1IovdZtjxagRQ5vrZyhU4oq2D0ZL0TeikboiyKciaeMxcuXEhbIeaNxdxsSIVa8f68a9vaoTXaWafZdjvd907LvXbuQ2K/ovVltKSN1piHDh1Ktz9Isd3aMbDjOVr+NuvpIk+H8Tfo+5CsPigbPaOI/RpU3VOn5UkjVai/aseJipdLAg2UQmTEkYHWX2QX1S0rMtVGBUJkCnmT2cS2s8KmvtBptr+NludlKLEsr0BrVcgNOnPqZX+y9+Z1ZQfKrZO8NW9S5ez//ZifQbkyujopU/q1bqN9UlbB6MibeLhTUc7Uz/GQ/d7zyp5sm62WR3rx/qhUqK24iQqfbPiEdtYpet87uZ9qdB+SJ9aLgEPsTwQbsu0P4p4lb5iS2I/Yh0YVZPJ0GH+Dvg/JhlDKKu2L2K9B1T11Wp400km50al+lydFpn+u4uWSQAOl0ejCvldZq5tmNw6d3Gi0er2TzxHb7eZzR+bUrwy+G+3szzhEboFvWsR0k7fWVkzU5wvtVLh0Q7lSTp2UcUWvq6yC0VE7z0Gj32u0HI0KnVdeeSWZmJhIx8CurcSJSpVsnoJYv7b1ZjZGdTciL4kWpzE2cwQSYttRztUGFNpZp5lu9r3f9yFZ5VleRdkg71myyrJWE3fK02G8DPs+JJsDNMqevPxlVOueiioLu62/KkK/68B6TX9fhcolgQYGplGr/VrRDaxRxtSvTKPbiGGr/Wn0eiefI9ZttH/NXovjHF3iBqnX/SnD5KRAa9kFdKd56+nTpze7/maiW2+zCbGUK8WWK6OkkzKlX+vmUVbB6MiGGGp00x6/95jYMirzYwiIrAK8fnLLqEyJiv6LFy9u9igo4r4j0nn48GG6f1GWRR6T13uh1TrNdLrvRZUXUSmXNwREdh9XO3xSP+9ZIrAR+9JIq/JOng7jZZD3Id0MoTTKdU9FlIX9vs/odx1YEemfUi4lUz19C9CB6CrUrAvZsFriZ5lFo9ZQsU/DughtVokWhWGzWe+bjf/WD93uT/b5RqkHBtCduDjP8stO89bIAyIfiYvzrNLoyJEjTS/wlSudayc4Mwo6KVP6tW4tZRWMnmx+gBiDOsqaZtfr2XVmfYOmuDeJ4EM0durler++YiKrrIpyKntEgCMqIKKci/+3s04zRe17N+VFVOblbTM7DrUtgPt5z9JqP+I6Io88HcbPoO9DGon8L4bOyeq+IuA96nVPRZUn/b7P6HcdWBHpryqX9GhgcLIutHEz0CiCGhOkDFoUKhGtzZsELduvdrv/Fi1uhhp1VYubkLxJ5GpvWgap2/2JQlaQAcZDXEjXjh/aSd4aeUVcvEUaUdESj1YX98qVYsuVUdJJmdKvdWspq2A0RVkRlTqtJt9s1OshypGoWKivWKnNI+rvXRoNh1G/vHayztpeGLVBiXbWaaSbfS+qvIi8MiqkGh3n2nu6ft6zxHVG3nebDfmRN052kKfD+Bn0fUirfYm8ub532qjWPRVVnvT7PqPfdWBFpH9SuSTQwOBEphUX09EVqz7DjQz1xo0b6Y8yT1xs513UN1oeOln+4Ycfpsvre1zEhX5kJnkXqc0itZ3ub7NuyY26Tsd+xfHKurZlolCM41xfeR9j08a4tN3K9rvRvna6P5mIGA96mCegc61ap0R+GXl57QVqJ3lrvC8rH+ICLyoBWk0+plwptlzpJM/vtVwpskzp17q1lFUwuqLCO64lIz/Km2Qx8ryo/MhrpZk3HnakEXlFXt4alRB560fa9flZlH3160aleG0lfDvrZOr3p9N97/T+pFl5EZVoZ8+e3fK+2JeoiIqyvfaerp/3LDEkVjaEViauIWJb8nQYH6N2HxKv1Q4Rl6e2J8Oo1z11Wp402qci7zOiHKgvC/pdB9Zp+nn7eEG5lExsbGxstH3UoQCRqUTENDKVrPtyZEh5GWpkdrFuZFbxvvjRRwYcBUHt8sgQYnm0aKpfP37okSHUrx/r1o/RmkWPp6enk7W1tW37lV24ZjcxkX5tl7p297fV58guziO9yLia/Uzjgj6ixtk+RyQ+L2CTTXDUaXe42L+sIM1a/8R+xneXV3i2uz+Z+LyNvn9g+CJvisqE2gvq+E1HHpBdwEYlbHbBmJdftcpbM43GpY7txUV5o0pg5Upx5UoneX435Uo/y5R+rRuUVTD6suv0bBilLF+Ka/X4fUeeWD8sRVbGZXMjhCwvyO4tsnuJ2vwjGkhlw/LEOlEB8dFHH6VpZ+vH/sRzbWVMvJ6ln1VoNVsn794i0o9ttrvv8Zn7cR+Slb3Z/sf6te8f1D1L7X5kFWXZMWpEng7lMGr3IbGdCDhn24x8JvKTRnUZ8VpsPy8vG6W6p17Lk6ylf9H3GbE8r9dav+vAOkk/bx8fVbxcEmiAEReZTQREihgKKSucm2VwgxYFUWTeQLXFxVXkB7XdlLObh+xCPG+YBjo37uVKPyirgCoqsrwYpbJFng7Uch9S/nKj7PcqM2NULzY57B0AmovIZ6Nud52KCrthTWydJ6Lzo57hA/0XrWIif6ofCzVr5XL9+vXc8TLpzjiXK/2grAKqqsjyYlTKFnk6UMt9yHiUG2W+V7k6ZvViAg0w4rKubK3GKm9HpDNKEy9nky0BNFM7sRu9G+dypR+UVUBVFVlejErZIk8HOuE+pBzlRpnvVZbGrF5MoAFKIDKe+onOOjVqLYFj8qMY627Uo8tA/2Vjc8Z4mHliLNRxuvgaBeNYrvSDsgqouiLKi1EpW+TpQD33IeUvN8p8r3JpDOvFzNEAJerSF5PMdFvZFhV49cOSDEsMgXL58mUVh8C2bqORz8UYlZkYqzICDaMyzuc4GadypR+UVQDFlBejULbI04Fm3IeUt9wo673K6pjWiwk0QIlERjQOY7eNy+cAKDv5sWMDUIXyouz7D1A28t1qHh+BBgAAAAAAoGtT3b+Vei9evEj+4i/+Itm1a1cyMTHhAAFd29jYSB4/fpz87u/+bjI5aTqdcaf8AIqg7KgWZQdQBGVHtSg7gH6WHQINBYogw9/8m3+zyCSBivvzP//z5Pd+7/eGvRv0mfIDKJKyoxqUHUCRlB3VoOwA+ll2CDQUKHoyZAd59+7dbb3n+fPn6USXMfHljh07itwdRpTvvHq6+c7X19fTwGWWrzDeOi0/5CPV4zuvpk6/d2VHNcuOf/tv/23y/e9/373EX5Nf5nNcHJNGlB3VknffIX+gGecHnZQdAg0FyoZLisy6k0DDd77znXR9gYZq8J1XTy/fuWHYqqHT8kM+Uj2+82rq9ntXdlRD9j27l9hKfpnPcXFM2s1TqN59h/yBZpwfdFJ2jESgYXFxMbl//34yPT2dtto6ceJEcvLkyY7TefToUXL27NnkzJkzDd8f61y8eDF9vnfvXvLgwYPkwoULuevHfsQM4JFePMf6S0tL6XvjGQAAAAAAqm7ogYa5ubm0O/jCwsKWCv4IAJw7d66tNE6dOpXs3bs3/fvq1atpYCBPBAjOnz+fbmvPnj3pstXV1eTQoUNpoOHKlStb1o99iCBIPDLHjx9Prl271tVnBQAAAACAcTPUQENU8l+6dCmdqbpWBAKi8r/dQEMWIIgeB5FeI9GToTbIEKKnQiyLAMTy8nIaSMjs27cvDVrcuHEj/TsCILWvAwAAAABA1Q010BDDD0VFf71sWfRO6GYIpUYivXjE8Ey1suBBBCxqAwnRS2J+fr6w7QMAAAAAwLiZHObGowdB9BTIE70Oih6iKLYVwyHlbSvkvQYAAAAAAIxoj4YY6qjRUETRm+DmzZuFbq9R4CKGcApHjhxp+Hrsy+HDh3N7YNQ7evRosmPHji3Lbt261XD29hcvXqTPVIPvvHpafecHDx7MfQ8AAAAAlMHQAg0xMXMz0cug1TpFDuEU26ufEyJ6OMTcDTE3Q7wWPTBi7ogYYqlRT4zw9OnTZHJya2eRO3fu5K4blY+xnbt37257D+PJd149rb7zJ0+e5L4HAAAAAMpgqD0aRkEED+IRwYPaSaJDFmDIRO+LmBw6ltfP81BrZWUl2b17d1vbj1bLUfm4f//+bb0gGE++8+pp9Z3fvn1727L19fW0ZxcAAAAAjLqhBRrqK/XrDao3w6lTp9IeDXmTTtf3cMiCDdHLodlE1VGR2EnQIFo4d/oeys13Xj2dfufyAwAAAADKYmTH6olhRloFI4oIMly4cCE3oNBINmRS0RNVAwAAAABAGQ116KQIJERAoVGPhph8uV+iV0JM/jw/P98wCBGTVccwSHka7TeD9aMPPkm+erx9fPvw6q6Xko/fedtXAoydf/gv/l3yxaP/mExN/TpJkoktr8n7AIBm1xC/+cunua+5hgCKpL4GqmeogYbTp08nN2/ebPh6zIXQD5cuXUqmp6e3BRlieda7YXV1NXfC5yzAEEEKhi+CDF+sfz3s3QAYqN/85ZPkN3/1PGYAceQB+mRxcTG5f/9+et8Q87PFvUmjoVOLSCNbN2t0dejQoY56XkO71xBfrOc31AIokvoaqJ6hBhqi10BU7seFdO0wSTE5czYfQtEi7dheXk+G2nkh4gZgYWFh2zoxN0Oo8kX/KEalJyeS5LVdO9O/v3z8dfJiY+C7ADBw8j6A/pibm0tmZma23A9EkCAaHbV7H9BJGrFu9LiubegU90mxPOaTg6K5hgAGRX4D1THUQEMEEqJC/+LFi1suwOPvK1eubJujIS7UQ7QGypMFChoNaxRDIcXFejahc+17stcyMXdD/YV99HKIfc3btyoZxah0BBk+ffdY+vdb710fuf0D6IfXdr0k7wMoWFzzRyX/xsbWlitxj9JuL4NO0oiGUHFvUd+bOtaJdesbZUF3DcQ2kmfPniUP/tM3vSHdPwGDIr+B6hhqoCFEpX10E46K/6xLcVTw53UpjgvsvXv3blse741AQVzQZ/+PyZpj3dpAQbQgivXioj9PbbAjthX/j32Jv+N94fr168ns7Gwhn73sRKUBABg3cf+Qd72fLYsezq2GUOokjbiHye416kXwIV5z/8G4NBADAMbX0AMNodGEzPUaTcycN8RRnkY9IRqJAIOuyo2JSgMAMG6ih0Gjiv24P4gGTa0CDZ2kEcGEaChVO19cJoIQggwU00DspbRHw9RUVAFMpEPeAgCMXaABAABgFEQPgkZzxUWP6Zs3bxaaRgQcItgQPamjt3c8wtmzZzf/hl4biP2P5/+L5M6dO8mBAweSHTt2OKAAQOEEGgAAAGrmb2skeiO0WqebNKLn9qlTp9KeEK+88kraiyGGbG1nbobTp08nO3fuTCYmJjaX3bp1q7Lf5/Pnz5MXL16kz2RzhGy0OC5b16uKKp4rBw8e3LasSp8fgP4SaAAAABiiCChEoCGeY/6GGDKp3R4NUUn49OnTLcui5XpVRcXxgwcPkrt37yaTk5NJlcVQSdlzHI9Gx6V2vSqdO1U8V548ySYH33ocAKAIAg1j5kcffJJO+pUnxuH8+J23B75PAABQBq16ELTqrdBtGidOnNgcOimGXYq/I+AwMzOT9nZolma854033jAcTk3gJSqO9+/fX/ljMjX16zgi6bwMcTwaHZfa9WJopaqo4rly+/btbcvW19fTId0AoFcCDWMmggxfrH897N0AAICxE62fYz6FItOIoEL0ZqidHDomi47JoeO1mCh6aWmpYXrREjsqSatSUdoOxySTDac1kZ4fjY/L1vWqxLmSVO47B6B/BBrG1OTEN5N+hS8ff528yIbdBAAAGoreAxEMaNQb4fDhw4WmEQGFjY3tF+vnzp1L1718+bJvCwCAkVeNgQgrKIIMn757LH1kAQcAAKD15MoxfFEjMczRINIIWS8HAAAYdXo0AKViHhIAoJ9iGKPoZRC9CWrnRlheXk6fjx8/Xmgas7Oz6fK8dGP5mTNnev5MAADQb3o0AKWchyTv0WgidACAdkWFf/QkuHjx4pblCwsL6cTL9RMzx4TN8eg2jfh/zMVQ3wNidXU1nathfn7elwcAwMjTowEoJfOQAAD9EpX/i4uL6UTM09PTydraWhoMyBvKKIIGe/fu7TqNbPLnCELUiuBFpAEAAGUg0ACUeh6S8NZ719MeDQAARWm3J8HKykrPaUSwYWlpqe19AwCAUWPoJAAAAAAAoGt6NAClmvj5y8d6LgAAAADAKBFoAEZ+4mcAAAAAYHQJNACF9TbIvLrrpeTjd97uy8TPtdsAAAAAAIZPoAEY+d4GtRM/AwBAVQyycQ8AQC8EGirKBSv96G0Q8ye82HBsGYzFxcXk/v37yfT0dLK2tpacOHEiOXnyZF/SePToUXLx4sX0+d69e8mDBw+SCxcuNNxeEfsGAGAoUQCgLAQaKsoFK/3obfDWe9fNqcBAzM3NJTMzM8nCwsLmsqjMjwDAuXPnCk0jggvnz59P19uzZ0+6bHV1NTl06FAaPLhy5Urh+wYAUEvjHmCcRCPFqD+opZcWlJ9AQ8W5YAXKJir5L126lGxsbO0+ExX7UfnfTmV+J2lET4baIEOYnZ1Nl0UAYnl5OTl+/Hhh+wYAUE/jHmCcxEgIgxqKGRicyQFuixG+YI1H/WS7AKNoaWkpreivly27evVqoWnE3xEkqJcFF2p7NBSxbwAAAOMoei18b/fOLY9oAAuMBz0aACiV6EGQV5kfotfBtWvXWs6H0Eka+/btS27evJm7XoghkYrcNwAAYPz1Oq9bNo9cvD/Sid7Vkcb8/HwyqvImsDcEM4wPgQYo0cTZo7QvMCwxGXPWm6De3r17c4MCvaQRwYE8cSEfjhw5Uui+AQAA463Xed2yeeSiR3Xtstdffz25cePGtnnkysC8DVB+Ag1QoomzR2lfYBji4rmZ6DXQap0i0ghxUR/rZjcCRaR79OjRZMeOHVuW3bp1a9t6GzXPz58/37b022WMi/hOX7x44butmGbf+8GDB3PXB6pDpRyUUxHzun300UdpGqdOndps6BT3G/F3DNcaDaCiZ3aZmLcByk+gAUo4cfYo7QtUUQyRFI9oKVQ7SXSvnj59mkxObp0+6c6dO9vWe/bs2eZz9nreMsZHVDZHC7e7d+9uO0eo5vf+5MmT3PWB6lApB+XUzrxurYZQOnz4cO59SPSiDu00nBoVMTpDPfUcUE4CDdDhxNmjMIbgKO0LDFKrSv12LqiLSCNaDsUNQu0NQBHprqysJLt372653tTUr6PtcjI1NZUcOHCg4TLGR7RUj8rm/fv3b+v1QjW/99u3b29bf319fbOCARhf/aqU00MCBqOIed3i/Q8fPsxNO3oyNEp/FJm3AcaHQAO0mAMhLriBcoiWv712EW6VRgQZLly40FaX5k73LSoS26lEnqh5/nb9b5eqiB5P0aK93XOEan7vzg2ohn5VyukhAYPRr3ndoidE3HNcv3695bq1Q7bGEE7Rs/pXv/pVUqxehnY1LOyoMIQrBzsYslWgAWqYAwFGX7TyiQvoRr0Gohtxv9KICddi8uf5+fm+7RsAwCAZtgQGp6j54mqDFhFgiAmg432fffZZW0O71g/ZmvWiLHKIzl6GdjUs7OgwhCtPOhiyVaABWsyB0OwCvJ/0rij++NV+l3ktwSiH06dPN23lc+LEib6kEZOtTU9PbwsyxPKsd0MR+wYAUMTwR+32zDZsCZRX9JiO+5MIMsR9SUwmHfPItRo6qXbI1n4N0dnL0K6GhR0dhnDldgdDtgo0QIs5EIZF7wrHj8ZDF8VFdFxM17bWifFIQ6NuyL2kEctj3byeDLUtjorYNwCAbhn+CEZbEfO6NUo37lXW1tbSYEM8Nxu2tX5Yxv4M0dnL0K6GhR0lhnClXqPfdHF9ooC+9a743u6dWx6D7l0xLscv/qb8orI+Jke7ePHiluULCwtp6536i/eZmZn00W0a0R15bm4uvViPoZPiEf+PR/RQqL2A73TfAACKEPcH9fcM7h2gfGIY1l7uGaLhU4h7FoBB06OBSnQbLvOQOaPQu6LMao9fEZPkMRqi0n5xcTG9gI7hjCIIEBX/UclfLy7U87r0tZtGBBMi2BA9FfJEEKHbfQMAKEKre5zOJ2IF+qWIed2i10Lc41y7dm3L8uy+Z3V1taC9BWifQANjS7dhGG+NJmTOG3+0lzQiUNCvfQMAAKqliHndIpCQNzRSFsBoNUcDQD8INAzZP/wX/y754tF//OuJbiYatr43sW378oYVit4NEXgAAAAAGJYi5nWLdaIndb2sh8OZM2cK3WeAdgg0DNlv/vJJ8pu/im6szbuymhi4t27DhswBAAAAhq12XrfaYVibzTlX39N6aWkpOXv2bPLhhx9urh/DvcYQrufOnTNsKzAUAg0jIiapjbHkW7W+b3e9ouT1pIjtwqjPy+E8BQAAYBznnIthkyLIEO+P16N3RAQaIl1zwwHDItAwIl7b9VJbE9YOemJbPSkoE/NyAAAAUIU55yLAED0bAEaFQANtqe1J0WwuBBiGZufioM/T/PlUNpJnz54l39vzVfJv/unfH+j+AABAEcwbCAA0I9BAW2p7UkAZ5uUYlma9gKam6gMQAABQDs2uc2N5NoRpbYOfUbpOBwD6S6ABoA8GPZ8KAAAM+jq3NvDQ72F9AYDRJtAAI2BYk27r/tw/W+dTWU6+WNebAQCA8brObXQfo5ENAFTPSAQaFhcXk/v37yfT09PJ2tpacuLEieTkyZMdp/Po0aPk7NmzyZkzZ5q+v5PtFbVvMIqTbpvsGwAA6Fbe0EgxhJLeDQBQPUMPNMzNzSUzMzPJwsLC5rKozH/w4EFy7ty5ttI4depUsnfv3vTvq1evpoGGIrZXxL7BoCbdjpZD3Y6LapgfAAAAAKCUgYbV1dXk0qVLycbG1sHLo2L/0KFDbVfmX7lyJX2+d+9eml4R2ytq32BQk25H9+RuWw5tHeZHCyQAAAAAoCSBhqWlpWR2dnbb8mxZ9E4ocpiiTrY36H3rRW1L9kGM689oyevxUMS4qL30kKjK3BDDmlsDAAAAAEbJUAMNy8vLuZX5Yc+ePcm1a9cKrczvZHuD3rde9NKSnfLr17io/TyvxmVuiHH5HAAAAABQ2kBDDHV0/Pjx3NdizoWbN28ObXu97NvRo0eTHTt2bFl269at3HU3ap6fP3++bWmzZd/9zu/ULN8qXus0vfaXNf8k2Xr/8F/8u+Q3f5nfav2733kp+df/5O8mvet+/waz3nax/osXL+re1316ve1fe+fVl4+fpIGHb3o5LG8u276Nzrb7zdwQL23ZRnfptas/6dV+jrzfYP7v/FsHDx7ctqyY/QMAAACAMQ40PHr0qOnr0Wug1Tr92l6v+/b06dNkcnJyy7I7d+7krvvs2bPN52yddpe9/w9ebbqfnabX7rJ2P8cXj/5j8pu/yq8sbZVeu3rZv0GslyeCDDGh+N27dzfPk17S62X/2j2v/quPPk+/y296OWwNHvVyvuz9GzuS/+7Hv7dlG92k165+pVf7OWq1+zmePHmSe54AAAAAQBkMtUfDuFpZWUl2797d1rpTU7+OtsvJ1NRUcuDAgY6WtauXbbS73WbvzWu13s3naHe7g1zvwX96nvyf/zT+br+3RrRUjyDD/v37N3u+9PL9FvE5Wq33vT1fJVNTjXunFHm+FHXet7vdQaTXar3bt29vW7a+vp72ngIAgGExHyAAMPKBhugV0EyRvRk63V6v+xaVx/VDJzUyUfP87Xu+Xdp8WbvaTa+X7TZ+72u7diafvnusbu6Abj5Hu9sd3Hp5Lfzb+WzRk2HredLL99v/z/tv/unfH8B2iz7ve9m/YtPL/503V8z+AQBA98wHCACUvkdDDC2zb9++kdzeoPeN0fNq3Xj8WWufb+YYAAAAGK/7nXZeAwCqa6iBhug5EJX2jXoNHD58eGjbG/S+US4fv/P2tmXf9tag7H70wSfJV5uTXW+/scr7/gEAYFy43gUAShVoOH36dHLz5s2Gr584cWJo2xv0vgGjI4IMgkYAAAAAUIJAw6lTp5JLly6lPQRq50VYXl5On48fPz607Q163xiNVusx/BFkvpnIfOfmuWFoLAAAABj8hPQZowzA6Joc5sajsv7kyZPJxYsXtyxfWFhIrly5sm1S5pmZmfTRapLmRkMedbK9TveN8rZar32oSKZWNpF5PLKAAwAAADC4CelrH42GOQaGb+iTQUel/eLiYnL+/Plkeno6WVtbS+bm5tJK/npRub93795ty+O99+7dS1ZXVzf/f+3atXTdpaWlrrfXybqULxqe9V6obbWeMcEZnZxDAADQKz2uARrXyRhlAEbf0AMNYX5+vq31VlZWcpdHL4N+bK/TdSlPNDyv1Tp0ew4BAECvzBMG0HhC+mjw514cRttIBBrIp+V0cZr1UCi694IxBEdfN9/RIM8hAACqS49rAKCMBBpGmJbTxWlWgVw039vo6+Y7GuQ5BABAdelxDQCUkUDDCNJyupyMITj6fEcAAAAAUDyBhhFUtZbTeZOeNdNqiJthMYbg6A8Z5jsCAAAAgOIJNDB0Jj2jlqGnAAAAAKBcBBoYmUmo8yY9q9+vqIRmPBkyDAAAAADKSaCBkWlR3mrSswh+DGO/GIxRHA4LAKiuxcXF5P79+8n09HSytraWnDhxIjl58mRf07h06VK6XmZmZiY5d+5cT58DAAAGQaCBobQoH1aviaqpPc6jPscFAMComJubSyv5FxYWNpdFkODBgwdtV/x3ksajR4+SU6dOpY9s/dXV1fT/hw8fTmZnZwv7bIzuPHXuiwCAMhNoYNMgK5+Nw+84AwCMoqjgj54FGxtbx+yMAMChQ4faCjR0msbZs2eTffv2bVkeAYl47Nmzp+fPxOgxTx0AMG4EGhgo4/AP7zib46I5vWwAgLC0tJTbgyBbdvXq1ZZDKHWSRgQl4v+1QyaF48ePJw8fPvSljLm8eeqK7k0OADAIAg0MlCF7hneczXHRnF42AEBYXl5uOFRR9C64du1ay0BDJ2lcvHgxXRY9GqieVvPUAQCUhUADUGl62QAAte7du5f2Jsizd+/e5ObNm4WmEetGkCHmaYjhlkJMIG0iaAAAykSgAag0vWwAgExU9jcTPQ9ardNpGjF0UgQlIsgwPz+/uTwmgl5ZWUmHYWrm9OnTyc6dO5OJiYnNZbdu3Uqq6vnz58mLFy/S59GWzd+xMZB9HexxGexnG/9zpTgHDx7ctqxKnx+A/hJoAAAAGIIs4BBDLdUHFGLi6OjVEAGHRr0jskrCp0+fbll2586dpKqi4jgm0b57924yOTmZjKpnz55tPg/i+xrkcRn0Zxv3c6VIT548yT0OAFAEgQYYsH/4L/5d8sWj/5hMTf06SZKJzYmIAQAYruht0EtvhU7TyNaNoZPq52jI/h8BiGaBhitXriRvvPFGsmPHjpb7VgUReImK4/3794/0MfnmXuB5MjU1lRw4cGCsjsugP9u4nytFun379rZl6+vr6ZBuANArgQYYsN/85ZPkN38V3VN1UQUAKJNo/dzrpM15aTQKTsTymMOhmWiJHZWkVakobUc5jkk21NXEwPZzcMdl8J9tvM+V/qryZwegWAINFRKt5t967/rm3wzX5ESSvLZrZ9sTEwMA0H9RuR/BgEa9EQ4fPlxoGtlE0I3WBQCAMhBoqJAXG0nyxboAw6h4bddLyafvHhv2bgAAUDe58s2bNxsekxMnThSaRgyL9NFHHzVct53ABgAADFs1ZjyquGgl/73dO3MfWtADAMC3YvLl1dXVbb0JYsLm0Gy+hG7SOH/+fLpe/RBJ8f4wNzfn6wEAYOTp0VABH7/z9rB3AQAASiGCACdPnkwuXryYLCwsbC6Pv2Pi5fr5FGZmZtLntbW1rtKIoZPm5+fTgMK1a9c2l589ezZdPjs727fPCgAARRFoAAAAqBHBgMXFxbS3wfT0dBpEiEBABA/qRdBg7969PaURAYhLly6lPSEirZjf4cKFC7nrAgDAKBJogIFPxP2kUsf8Rx98knxV95lNRg4AjLroTdCOlZWVntMI586dSx8AAFBGAg0wAFWeiDuCDFX97PRXtBK9f//+ZivRmFiz05afnaYRY2jHUBZnzpxpuF6kEcNcxDrxHGNuLy0tpe+NZwAAAIBxI9AAfZQ/2fZG8uzZs+S738l7bXxNTiTJa7t2bllmMnK6FUNPxJjYteNeRwV/DDXRbmvQTtLIhrIIV69eTYMIjcT7I4ARj9qxumvH3QYAABgXRjIAgkADDHgi7ufPnyd37txJDhw4UKljH0GGT989NuzdYAysrq6m41hvbGxsWR4Bg0OHDrUVaOg0jRhnO0TvhHhfMzGpZwQibty4kf4dwYsINAAAAIwjIxkAQaABgFKJ4YdiSKJ62bLocdBqCKUi0mgkej50MiY3AADAODCSAVSbQAMApbK8vJwbJAh79uxJhyhqFSQoIg0AAAC+ZSQDqLbJYe8AdOrLx18nb713PX3E30C1xPBF2XwJ9WL5zZs3B5JGu8MzxTMAAADAONOjgdJ5sZEkX6wLMEAVPXr0qOnr0Ruh1TpFpNFMTAZ9/vz5dG6GmOshek/EvA8xz0PM2dDM0aNHkx07dmxZduvWrW3rbdQ8x7wv9Uu/Xca4iO/0xYsXvtuKafa9Hzx4MHd9AACAYRBooDRe3fVSV68BDFIWYMjERNAxOXQsX1tba/rep0+fJpOTWzsbxuTx9Z49e7b5nL2et4zxEZXNEcS6e/futnOEan7vT548yV0fYNR6otfft338zttD2ycAoH8EGigNF6TlvKkwvBVFit4GzbTTE6GINJqpDTLUBhuil0OrSaZXVlaS3bt3t9zG1NSvo+1yMjU1lRw4cKDhMsZHtFSPyub9+/dv6/VCNb/327dvb1t/fX294bBwAIOmJzoAVItAA1A4NxUMS7T8bTU80SDSqJel12qS6ahIbKcSeaLm+dv1v12qIno8RYv2ds8Rqvm9OzeAUZDX2zwaH8U9AgAwvgQagMIY3opBiB4JEQxo1Bvh8OHDA0kjz6lTp9KJpqNnQp5G2wQAGOee6NHb2Tx7ADDeBBqAwhjeikE4ffp0cvPmzYavx1wIg0gjz+rqam5viCzAcOTIka7SBQAAABhlZhMEoFSi10BU6NfPpbC8vLw5H8Ig0sgTwyLF8Ej1Ym6GRvM3AAAAAJSdQAMApRJBgKjQv3jx4pblCwsLyZUrV7ZN9jwzM5M+ekkjkwUmGg2BdOHChWRubm7LsghoxHaapQsAjK8fffBJOnRQ7SPmLAAAqOTQST/84Q+T//6//+/7uzdA38VNTdzcZH9DP/Wr7IhK+8XFxeT8+fPJ9PR0sra2llbw5020HJX7e/fu7SmNWCfmXoigQfb/6LkQ6S4tLW3ZVgQrIp34O94Trl+/nszOzhZ8FABwj0IZfPX4ifkJoKSUMwB9CDREhcq/+lf/Kvmv/+v/uoPkgVHzYiNxo8PA9LPsmJ+fb2u9RhMzd5JGBA/aFQGG2uADAP3jHoUymZxIktd27dyy7NVdLw1tf4DWlDMAfZoM+u7du8lPf/rT5Lvf/W46zvTu3bs7eTswRM1uYtzg0E/KDgCUM/BNkOHTd485FFAy7mcACg40xBATP/nJT9K/f/vb3yaXL19Oh4Q4cuRI8uMf/7jdZIAh+fidtx17Bk7ZAYByBoCycj8D0IdAQxZkCC+//HJy9uzZ9O/PPvssef/995P79+8nZ86cSd54442kUzFGdrw/GyP7xIkTuWNk95pGNiFnDIGxb9++lmlGGjGmdnyueI7ASgyHEZOBGhYDYLhlBwAoZwDoJ+XMaM87WTtKg8aVULKhk+p9/vnnyaVLlzYr369evZpW7s/MzCSnT59ua2ilmDAz1q8d/zoq+B88eJAOz9SOdtOIQEHsYzwaic+SvSfeHwGMeGSOHz+ejtEHwPDKDgBQzgAwDO5nhsu8kzAGgYYLFy6kvQHCH/3RH6UVRNFDYGNjI62Yj8r+N998c8vQStG7oNmwSvH+qGyKNGpFwODQoUNtBRo6SSMCDREoyOvNEEGFeL12/VgvWtreuHEj/TuCF/F+AIZXdgCAcgaAQXA/Mzry5paM3g0ReABKFmjIKoeWl5fTCqKocP/oo4+2dCOrHx4jKo1iaIx//s//ecM0Y0iietmyrJVrq/1qN40YoqNRb4So7Prwww+3LNu7d28yPz/fdPsADLbsAADlDACD4H5mdOQNjRRDKH2x/vVQ9gfoIdAQw1vE3Ac/+9nP0laoUSHUSqt1ouIpL0gQ9uzZkwYFWgUaOkkjeiTkiWBEvBbrA1CcfpQdAKCcAWAQ3M8A9CHQEJX5N2/ezB2bLlr+14+p/dOf/jT5kz/5k6aBgmwoozyRZt72ekkjb714fwyNVDu/Q71ojRvpHD58uGFQo9bRo0eTHTt2bFl269at3HU3ap6fP3+elNe3n+Tbz5G3jDgWL168cEwqpNXv/ODBg9uWjctvph9lBwAoZwAo+/1MzAcaI2/E0LHROCsawHZyHxRBkBimNp6jbiuGBI+hntxLASMfaPjud7+7bVkMbxGZ4crKSvLZZ5+lGfAPfvCD9LWY2yBEJpcnMsJmondBq3WKSOP8+fPJlStXcl+LTDpej8w+WuJG74mY9yHWz5vnIfP06dNkcnJyy7I7d+7krvvs2bPN50brlEHe5xiXz1a0CDLEuXX37t1t5wnjqdVv4cmTJ7nnyTgouuwAAOUMAGW/n4nhu2dmZrY0eo26p6graGe+0qjrivqqeH82Okc0ko06qwg0NKrnAhiJQEN9lDYb3uLYsWOb/48xtbPMNSb3zCb4HFUxZNKRI0cavp4FGGp7REShEcujUGkkCpu845VnaurX0XY5mZqaSg4cOJCUVd7nGJfPVrRoqR5Bhv3792/r+cJ4avVbuH379rZl6+vraQuZshvHsgOA0aGcAaBs5UwEBC5dupTOYVcrggYRKGgn0BA9GWqDDCECHrEsAhDRULbR6B8AQw80TExMtFwn5kNod/LOVvMhtOqJUEQakfk2mhw65GXukVHH+5pNVB2Vx+1WIE/UPJe70vnbT/Lt58hbRoieDJ2cJ5RbN7/zcTk3ii47AEA5wyj70QefJF893tpb9cvHJiqtPRYxeWutV3e9lDvJK4zr/UxMMJ03LHe2rFl9UybWiUd9I9gsuBA9GgQagJEJNMS4cjHGW3QJCzEm3Q9/+MOGCcXr7URd2xXdxZoNT9RrGhFBjs/X6Tay9duZqBqgaoZddgAw3pQzjLoIMnyxLrDQyIuNxPEhqXo5E70NGs3/GQ1q26lvirqpvLkjsga5UR8GMDKBhp/97Gfpc0RII9OMKG59t67ajCzWP3v2bEcbj/c1yvyiN0JMvtyvNCKC3CzIcOrUqbRwiWGQ8si0AYZTdgBQXcoZymJyIkle27VzW8v9qsr77NG7IQIPULVyJuqaGvU2iOFz8wII9RqNzhGNakOzYcLD0aNHN3vRx+eLuUZ/9atfJd3LjtFGOlz04Axru9URxzXmkHR8q+vgwYPbljU6H1oOnRRR1Ii0RlT3o48+Sop0+vTpphlozIXQrzQigtxs6KXInPMCEVmAoVWmDVBl/Sw7AEA5w6iLIMOn7347hnvV5Q2NFEMo6f1B1cqZVkN8Rz1VO0OJNxKNaiONVr0sIrAQQzrXzyVZu6wTz54923y+c+dOV2mUabtVEkGGqAvt5fyg3J482TokZHZedD1HQ1S4xyTIRYteAzEBTmSitZX+EQQI7Ywn120azSLIWaESk+jUi6h2MNQHwHDKDgBQzgDQb2W7n4l6sHjE/Ayt5jSN0Tuyia6zIMP+/fu7nitwaurXkVIyNTWVHDhwoKs0yrTdKini/KDcbt++vW3Z+vp62gOr68mgf/KTn7Rc55e//GXyxhtvtJtkWtEfFfoXL17cUqkff+dljDMzM+lz7WQ3naYR2okOX7hwIZmbm0ujwbW9HGI77WTaAPSn7ACAjHIGgDKVM63qknrpzRANcaMOq535RKPCuLbSOFqq1y/rTDZp9sSAK6OHtd1q6f38YNw0OhfaDjS0IyrhL1++3NF7otJ+cXExOX/+fDI9PZ0GEaKCPy9jjAw5L1rSSRpZOvFoNPlOtk4EKyKd+Dt6QITr1683fV/VxTib0QU2+xugH2UHALRLOQNAWcqZGKKm2XyizYIM0WDW6BvAMOUGGt5///103oM//uM/3lzWak6CiLpmlfGdmp+fb2u9RhMzd5JG5uHDhy3XiQBDbY8GWovJvIyzCdU06LIDgGpRzgAwDuVM1DVl83/mpXf48OGO0otGt7GfndaLjWuj19pJ6PPmhwEGHGj4l//yXyaff/75lsw1eglEZtcoshqz1DfKKBl/kYF38xowPpQdAChnACirQd3PnD59Og1oNHLixIm204o5S2Nkj/ogQyyvUu8GjV5hhAMNMRdBfUYZGevPf/7znnsJMJ5EiQFlBwD9pJwBYBzKmRjmKAIB0Xuhds6GmMg5m4u0HbF+pJHXk6GXuR7KJK9ha/RuiMADMCKBhpdffjl91GpnCKHayZgBqBZlBwDKGQDKalD3MxFIiDlFY26H2vfG3zEHaf2E0TMzM5u9KzIxXFPMKRppxdBJtcGF7LWqNnqNIZQM6Q3D0fZk0K+//noh6wBQHcoOAJQzAJRVv+5nIqCwuLiYBgli6KMIIkRwIAIQ9SLwsHfv3m3DK0VAIXpG5NEQGBjpQMMvf/nLzb/feOON9PnP/uzP0swxupadOXMm+fGPf9yfvQSglJQdAChnACirft7PtDt588rKyrZltb0bAEbFZLsr/uEf/mEaaY2Iabh+/XoaQY2MNSaY+Z//5/85+dM//dN+7isAJaPsAEA5A0BZuZ8B6EOPhhgTLjLYTAQd9u3bl1y+fDn9/7Fjx5L333+/g00DMO6UHQAoZwAoK/czAH3o0VBvdXU1d+w4AFB2ADAM7lEAUM4AjHig4ZVXXtn8O4ZNmpiYSIdOqrWxsVHs3gFQasoOAJQzAJSV+xmAPgQa7t69u2X2+pj1/gc/+MHmss8//3xLBgwAyg4A+kk5A4ByBqBkgYa5ubnk8OHDyYEDB5Ll5eXkww8/3Ozd8Ad/8AfpuHWPHj3q574CUDLKDgCUMwCUlfsZgD5MBv36668nN2/eTH7xi1+kk0C//PLLm6+dOnUqfcRwSgCg7ABgENyjAKCcAShZoCHz5ptvbvn/sWPHNv/+0z/902L2CoCxouwAQDkDQFm5nwEocOikdiwtLRWZHAAVoOwAQDkD1fXl46+Tt967vuXxow8+GfZuQdvczwB00aPh/fffTy5fvpz7WszPcO/evU6SA6AClB0AKGeARl5sJMkX6187QIws9zMABQcafvrTnyaXLl1KJ4SOORrq3b9/P3nw4EG7yQFQAcoOAJQzQJ5Xd72U27shAg8wKtzPAPQh0BC9FVoFEk6fPt3BpgEYd8oOAJQzQJ6P33l727IYNknvBkaJ+xmAPszRcOLEiZbrLCwsdLBpAMadsgMA5QwAZeV+BqAPPRpiDoZWPvvss+T111/vYPMAjDNlBwBlLWcWFxfT4WGnp6eTtbW1tLLp5MmTA0vj0KFDycrKSsf7DUBx3M8A9KFHw7lz59IJcD7//POG6ywtLXWwaQDGnbIDgDKWM3Nzc5s9tufn59M04hFz1g0ijXjv6upqx/sNQLHczwD0oUdDZK4RyT1//nyyZ8+eZO/evelzJl6LsesAQNkBwCD04x4lKvgjGLCxsXVG2ggYRC+D2GY/04j33rx5s6N9BqA/1IUB9CHQcO3ateTw4cPJT37yk/QCvl50CW41WTQA1aLsAKBs5Uz0Opidnd22PFt29erVlsMf9ZLG5cuXkzNnzujRADAC3M8A9CHQsG/fvuTnP/9503VOnz7dwaYBGHfKDgDKVs4sLy/nBglC9JaISqdWgYZu04g5HS5cuNDREE0A9I/7GYA+BBo+/PDDlutEV2AAUHYAMAj9uEeJoZaOHz+e+1r0mmhnWKNu0oghk6JCq3bop3ZEIGXnzp3JxMTE5rJbt24lVfX8+fPkxYsX6fNwZMNlbQxxH0bxuIze8SvPMSnOwYMHty2r0ufvhrowgD4EGt58882W67z++usdbBqAcafsAKBM5UzM6dBMBAFardNtGjFkUjcNt6KS8OnTp1uW3blzJ6mqqDiO4bLu3r2bTE5ODnz7z54923wepe9h2MdlFI9fWY5JkZ48eZJ7HGjM/QxAHwIN4Ze//GU60Vq0wImL4P/mv/lv0uV/8Ad/kLak+cEPftBJcgBUgLIDAOVMc9mQSd24cuVK8sYbbyQ7duxwov114CUqjvfv3z+UYzI19evYi2Rqaio5cODAyHwnwz4uo3j8ynJMinT79u1ty9bX13PnuOFb7mcACg40/OIXv0iOHTuWdgH+2c9+tuW1P/zDP0z+5E/+JM184yIXAJQdAPRb0fcorYYtatVboZs0uh0yKRMtsaOStCoVpaN/TLIhrCZG7jspx7ky2ONXjmPSX1X+7O1QFwbQh0BDXLivrKxsdj3+oz/6oy2v/+QnP0nef/99gQYAlB0ADMSg71FimJUIChSZRrdDJjF8P/rgk+Srx1uHovny8ddD2x+geOrCAPoQaIiLd3MwANCJfpYdMczE/fv3k+np6WRtbS05ceJEcvLkyb6mEa1Qz549m5w5c6bpekXsGwDDKWeiZ0EEAxqVA4cPHy4sjatXr6Y9Gubm5rask00WnS2PQES3PR7onwgyfLEusADjTF0YQB8CDd/97ne3/H9jY2PbOlGpAgD9Ljui4mVmZmZLC9CozI9KnXPnzhWexqlTpzbHro1KoQg09HPfABheORNzz2UV/XkiTy8qjQhC5wWioyyJAMTS0lLb+83wTE4kyWu7dm5Z9uqul4a2P0Bx1IUB9CHQEJMk/ft//++T73//++n/JyaysRO/EV2SAaDfZUdUvFy6dGlbZVJU7B86dKityvxO04iJNsO9e/fS9/Vz3wAYbjkTweXIy6PnQW0vguXl5fQ55oMYRBqURwQZPn332LB3A+gDdWEAfZqjIcYSjdY3R44cSYeCiNadUekSLW3iAvrGjRsdbBqAcdePsiPeNzs7u215tix6HLQapqiINAaZLgCDK2ciCBB59cWLF7f0Tou/I/BcP4RR9GILse1u06jXaNglAAZLXRhAHwINcTEc3X+jG+/8/Hy6LOvKG/+PzBcA+l12RGvQvMr8bHvXrl1rWZlfRBqDTBeAwd6jRDAg5ts5f/785nw7sY28PDz2IRter9s0MtELIsqKCEyH6A0X8zkYQglgONSFldeXj79O3nrv+rZh7T5+5+2h7ROMu7YDDSFaC8WF729/+9v0gj4uqN98883+7R0ApVd02RGtVBsNORFpNxsTu8g0+pHu0aNHkx07dmxZduvWrW3rbdQ8P3/+fNvSb5cxLuI7ffHihe+2Ypp97wcPHsxdv4r6dY+SBS5aWVlZ6TmNTAyxZ5i90fSjDz5JJ3+ur8QCxp+6sHJ6sZEkX6zLp2FkAw2Zl19+OTl2zBiUAAy27Iixrlu1OGq1ThFp9Cvdp0+fJpOTk1uW3blzZ9t6z54923zOXs9bxviIyuYYSiXGCa4/R6jm9/7kyZPc9avMPQr9FEEGFVZQbcqZcoheC/UiMByBB2CEAg2//OUv05udaCUUEd3du3f3b88AGAvKjvZFq9h2ytapqV9H2+VkamoqOXDgQMNljI9oqR6Vzfv379/W64Vqfu+3b9/etv76+nruED7jTjnDIE1OfDP5c6tKLWB8KGfKJW9opBhCSbAYRiDQEDcsZ8+e3RwntFZMuhbjhf7+7/9+v/YPgBLqV9nRavLMdnoiFJFGv9KNisR2KpEnap6/Xf/bpSqix1O0aG/3HKGa33uVzg33KAxLBBk+fVfvfhh3yhmAzjXte//+++8nr7zySjqR2euvv56OdRrDXsRz/P/nP/952rPh3Xff7WLTAIyjYZYd0euuVYX/INIYZLoAVeMeBQDlDECJejR8+OGHyXvvvZf87Gc/Sycki7Ho6sWEa5cvX05++tOfJtPT08k/+2f/rN/7C8AIG0TZEZX1UWnfqNfA4cOHB5LGINMF4BvuUQDoJ+UMQMGBhl/84hfJwsJCOlZ0tD5tJCqQoiLp9OnTaevU48ePJ9///vd72B0AympQZUe87+bNmw1fj6GZBpHGINMFwD0KVElM3BpjqtfPhZE39joURV0YQB8CDdESNcbPblZRVN+C86OPPkpbsUYr1U4tLi4m9+/fT1u2rq2tpRUxJ0+e7EsasXx2djY5c+ZM+nzv3r30s0ZL03jux74BVMGgyo5Tp04lly5dSvPt2qGIlpeX0+cIXAwijUGmC8Dg71GA4XmxkZi4lYFTzgD0YY6GqFCP8bQ7EZUnDx8+7HgH5ubm0udoBTs/P5/ePMQjKmr6kUYMaRHBg0OHDiUTExPJzMxMsrq6mhtkKGLfAKpiUGVHvCcCvhcvXtyyPPLqmBeifh6EyOfj0Usa9RM6Nxoeqdt0AWhtkPcowHBEr4Xv7d655TE54dtgMJQzAH3o0RBDWXSj0wqUqOCPSvuNjY1tFTIRCIihNYpOIz5b9Ga4ceNG+nf0UMhrYVrEvgFUyaDKjhCV9hE0Pn/+/GaPswgO5/U4i/T37t3bUxqxTvSAi7Ih+/+1a9fSdOsD1Z2kC8BoljPAcOQNjRRDKH2x/vVQ9odqUc4A9CHQEC39u5FXkdNMVM7E8EX1smVXr15tWTHTaRqxj9E7YRD7BlAlgyo7Mu3k5SHmjOg1jQgyd6LddAEY3XIGgGpRzgD0Yeik+lb8/cqUY8zqRhHjaHkUrUUHkcYg0wUYV4MqOwCoJuUMAMoZgIr0aGg0ZnUjMQxFo4kxo+XRzZs3+5ZGDH8Rrx0+fDi350Iv+3b06NFkx44dW5bdunUrd92Nmufnz583TJPxEd/zixcvfN8V0up3fvDgwW3LypgfDKrsAKCalDMAKGcAShZoiPGlP/vss45bDWVjV7cjm1Czkeg10GqdbtKICq0YNzvmZoh5FqLnQsy5EJ8568HQ6749ffo0mZzc2lnkzp07ues+e/Zs87nROoyXCDLEeXj37t1t5wnjqdXv/MmTJ7nnSdkMouwAoLqUMwAoZwBKFmgI7fQmKOPwF1mAIRO9FmJy6FgeE3YWIcYD3717d1vrTk39OtouJ1NTU8mBAwcK2T6jLVqqR5Bh//7923q+MJ5a/c5v3769bdn6+nopx5Qe17IDgNGgnAFAOQNQokBDDCXUzUV8DEPUrugV0EyrXgXdplEbZKgNNkQvh2yC5173LSqP261Anqh5VulcHdGToZPzhHLr5ndexnNjEGUHANWlnAFAOQMwunLHbWk0CXIr3b4vTwwt06rCv6g0sv1ud4LnIvYNYNyMQtkBwPhSzgCgnAEoWY+Gjz76qKvEOn1fVNY3mgQ0eg2008q1kzROnTqVTvIcQxvlqU2niH0DqJJBlR0AVJNyBgDlzPD96INPkq8eb51n8MvHXw9tf4DRMdSZaE+fPp1W/DcS8yYUmUZMOJo35nkWUDhy5Eih+wYAAAAA4yKCDF+sf73l8WJj2HsFJFUPNEQPg6j8r5/zYHl5eXPuhCLTiPkX8oZHirkZ6udvKGLfAAAAAGDcTE4kyfd279zyeHXXS8PeLaCqgYaorI/K/4sXL25ZvrCwkFy5cmXbPAgzMzPpo9s0Lly4kMzNzW1ZL4IJ8d76dTvdNwAAAACogtd27Uw+fffYlsfH77w97N0CRm2OhkGKSvvFxcXk/PnzyfT0dLK2tpYGA6KSv15U7ucNfdRuGvH+CBTEa/F3NjTS9evXk9nZ2Z72DQAAAAAAqmjogYYwPz/f1nqNJnHuJI0IMCwtLRW+bwAAAJSfiU4BAEoaaAAAAIBRmugUAID2CTQAAABAzkSnMQZ5LROdAgDkE2gAAACABhOdAgDQ2mQb6wAAAAAAAOQSaAAAAAAAALom0AAAAAAAAHTNHA0AAAAAwFj78vHXyVvvXd+y7NVdLyUfv/P20PYJxolAAwAAAAAw1l5sJMkX618PezdgbAk0AAAAACNNS2SgW9FrIS9PicADUByBBgAAAGCkaYkMdCtvaKQYQknvBiiWQAMAAAAwkrREBoByEGgAAAAARpKWyABQDpPD3gEAAAAAAKC8BBoAAAAAAICuGToJAAAAAAZocXExuX//fjI9PZ2sra0lJ06cSE6ePNlxOo8ePUrOnj2bnDlzpqv3AxRFoAEAAAAABmRubi6ZmZlJFhYWNpdFoOHBgwfJuXPn2krj1KlTyd69e9O/r169mgYaAIZJoAEAAAAABmB1dTW5dOlSsrGxsWV5BB0OHTrUdqDhypUr6fO9e/fS9ACGzRwNAAAAADAAS0tLyezs7Lbl2bLonQBQRgINAAAAADAAy8vLyb59+3Jf27NnT3Lt2jXfA1BKAg0AAAAAMAAx1FE2t0K9WH7z5k3fA1BK5mgAAAAAgD579OhR09ejR0OrdYpw9OjRZMeOHenfMVfE06dPk1/96ldtvjubW2Ijef78eVJe4/I5+iuOzYsXLxyjCjt48OC2ZY1+MwINAAAAAFAREViYnJzcUml49+7dLcsaefbs2ebznTt3krIal8/RbxFkePDgQdvnB+PnyZMnuedFHoEGAAAAAOiz6LHQzCB6M4SVlZVk9+7dW4IM+/fv3+zl0MzU1K/jXcnU1FRy4MCBpKzG5XP0W6fnB+Pn9u3b25atr6/nDgEn0AAAAAAAQxYtxxtNFF2kqDCurTSOlur1yxqb2Hwud8XzuHyO/uvs/KAKGp0LAg0AAABU0o8++CT56vHWIQG+fPz10PYHqEavhggoNOrRcPjw4YHvE0ARBBoAAACopAgyfLEusAAMzunTp5ObN282fP3EiRO+DqCUzOIBAABApU1OJMn3du/c8nh110vD3i1gDJ06dSpZXV3dNh/D8vJy+nz8+PEh7RlAb/RoAAAAoNJe27Uz+fTdY8PeDaACIpBw8uTJ5OLFi8nCwsLm8vj7ypUr2yaMnpmZSZ/X1tZy08sCFo2GYwIYFIEGAAAAABiQCCgsLi4m58+fT6anp9MgwtzcXBqAqBeBh717925bHu+9d+9e2jsi+/+1a9fSdZeWlgbyOQBqCTQAAADUiQqg+/fvb1YAxZjZeRVARaQRrVGjZWs8R6VRtEq9cOFCx9sDoDzm5+fbWm9lZSV3eW1vCIBRINAAAABQI1qVxlAVtZU4ESSIAMC5c+cKTSOCC9EKNdbLhsuI1qmHDh1KAw3R6pVi/OiDT9LJn2t9+dhE0AAARRBoAAAA+GtRyX/p0qVkY2NjyzGJQEBU/rcTaOgkjWyM7toxuWdnZ9NlEYCIyUFNDFqMCDJ8sS6wME4iUPTWe9e3LItJvD9+5+2h7RMAVNXksHcAAABgVMS41lHRXy9bdvXq1ULTiL8j+FAvCy7o0VC8yYkk+d7unVseUTlN+bzYSNLgUe2jvtcKADAYejQAAAD8tehBkBckCNHrICbabDV3Qidp7Nu3L7l582bueiGGWqJYr+3amXz67jGHtcTyAkPRuyECDwDAcAg0AFBKg5yks5N1Y3lULp05cyZ9jkk9o2VrjMEdzwCMtsi3Gw1VtHfv3tygQC9pRNCh0fBL4ciRI023dfr06WTnzp3JxMTE5rJbt24lVfX8+fPkxYsX6fN2WS30RoPXq3pcyuf/9X/5O9uW/b2Ff5t8sf6k7e933I5JOw4ePLhtWZU+P+QxBBsUR6ABgNIZ5CSdna4byyIoEY9MVDY1qkgCYHREULiZ6GXQap0i0ggRnI51W5VrUUn49OnTLcvu3LmTVFVUHEdZfPfu3WRycutIwc+ePdt8rtoxanZcxkWn328Vjkm9J0+e5B4HqLJsCDagdwINAJTKoCfp7HR7MQRG9Ga4ceNG+ncEJEziCUAnYuileMT8DLWTROeJdd54441kx44dDvJfB16i4nj//v3bjsnU1K9jjWRqaio5cOBApY5Xs+MyLjr9fqtwTOrdvn1727L19fW0pxVUjSHYoHgCDQCUSjsTbLYaQqmTNDrdXtyozc/Pd/y5ABi+VpX67fREKCKNU6dOpeVPO0MCRkvsqCStSkVpOxofk2x4qYlKHq/xP1c6/37H/5i0VuXPTrV9/M7b25a99d51vRugB9XoHwjA2IgWntFToNkEm0WmUcT2ABgPMcxKq0BCr2lEkOHChQttDwUIAACjQI8GAEpl0JN0dru9GHIpXjt8+HBuj4g8R48e3daqLG9Cz2+nsqydwK+6E1xWQRUnrKT5925Cz/6JIEAEAxr1Roh8vV9pnD9/Pp38Wc84AADKZiQCDTFh5v3795Pp6elkbW0tHc+6nW7C3aQRF/YXL15Mn6PyKG4AosVQ3rqRRlQOxVjb8RzrRxfmeG88AzDek3R2s70oV6KiKMqQaI0aPSJiLocYQ7tRz4hMTOZZPxlh3mSGeZMdVnmCyyqo4oSVNP/eTejZP6dPn24atI78vR9pxHxAcS9TH2SI5Xo3AAAw6oYeaJibm0tmZmbSSTVrL7zjpqrdC+p204jKoKj8ifWy7srR4jQqgCLQEJVAteL9EcCIRyZatRomA4BGsgBDbbkRAetYHoHwZlZWVpLdu3d3NdlhlSe4rIIqTlhJ8+/dhJ79E0MXReV+3DvUDnEUgePQqJdbL2nE8lg3rydDO3M6AABApQMNUckfF+AbG9lwD9+IQEBU/rcTaOgkjejJUBtkCNFTIZZFACIu8Gsv+qPlaVQO3bhxI/07KonaubEAYDwm6exme3llV5QdUc60mqi63ckIv53qsHYCv2+Wfvn4SfL3Fv6HLeu/uuul3MnOKBcTVlZTJ9+7IFQxIs+OvDq7d8jE39Ewqb5siAZPoTaY3Eka0Ws6Gk5lZUVt+ZK9BrTvy8dfpxO61nItBABjHmiI4Yfyxq3OlrWqkOk0jfg7HvUtSrPgQVz01wYSYuxt46MClEf0RGs1PFGRabS7brZO9IjrdGjATr3YSJIv1r/u6zYAxl3cF0Sv5qj4z4ZmjQr/vDw8ggZx39BtGtGYKQIK0XgqT22gAmjNtRDQKwFLKGGgIXoQNJogMy7Y26mQ6SSNqOjJGys1a1HUaMI2AKo7SWcn68ZQGVFZFEMg5elnORMt9fIukONmG4DOtdvgqFGe324arYbVA9rjWggoioAllDDQEJUxjYYiilZBzSZQ6yaNRnMrxPBL4ciRIw1fj3SiMqlRUKPW0aNHt3Vdv3XrVu66GzXPMQ4v4y++55jc0fddHa1+5wcPHty2zPkxOpN0drJulBd5PRyyAEOjcqYIeUMjxbABejcAUEX/8F/8u+Q3f/l0WwCe8eVaCOiVgCWUNNDQahztaEHaap0i0siGX4p168fVjoqh6OqcTewZvSdi3ofoBt1sqIynT5+m4+nWunPnTu66z54923xutA7jJYIMcW7F5I715wnjqdXv/MmTJ7nnCaMxSWcn60YPurwhLmLYvtDO3EMAQO9+85dPki/Wt19jAUAjApZQ4h4NoyAqiuKRN7FbFmDIRGVSTA4dy5t1cY7u07t3725r+1NTv462y8nU1FRy4MCBHj4JZREt1SPIsH//fpM2VkSr3/nt27e3LVtfX88d75nBT9LZyboXLlxIx9+OAHZtL4d4b96+AQD9NTmRJK/t2tmyxSoAACUNNLSqbGmnJ0IRaURL1agQypsLIq/laVQ4RS+HZhNVx7BJ9UMnNTJR89zueyi/6MnQyXlCuXXzO3dujM4knZ2sG9uKAES8Fn/H8H7h+vXrbQ29BwAUK4IMn757zGEFAKhqj4YYWqbZ8ERFpBFBhmh92slQFll67UxUDUD5J+nsdN0IMNT2aAAAAAAYd0MdID4qY7IJMvN6I8Tky/1KI1qlxqScjSqOIggR8zE00mibAAAAAABQJUMNNJw+fXpzWIk8MRdCP9KIST1j6Iv6IEMsrx1TO2+ojSzAEEEKAAAAAACouqEGGqLXQFTo18+lEJMzZ/MhFJ1GLI9183oy1KYRwyLF8Ej1Ym6G0MlwSwAAAAAAMK6GGmiIIEBU6F+8eHHL8phIMyberJ/seWZmJn10m0b0fIgJOmMSzxg6KR7x/3hEz4fa+Rxi7oZYXisCGrGdvH0DAAAAAIAqGvpk0FFpv7i4mFb6x3BGEQSICv68iZajcj9vOKN204hgQgQbaodIqg9O1G4r/h/pxN/Z8EzXr19PZmdnC/jkAAAAAMCo+/Lx18lb713fsuzVXS8lH7/z9tD2CUbN0AMNodGEzPVWVlZ6SiMCEJ2IAMPS0lJH7wEAAAAAxseLjST5Yv3rYe8GjLSRCDQAAAAAAIyS6LWQ17shAg/AVgINAAAAAAB18oZGiiGU9G6AEZsMGgAAAAAAKDeBBgAAAAAAoGsCDQAAAAAAQNcEGgAAAAAAgK4JNAAAAAAAAF2b6v6tAEAZfPn46+St965vWfbqrpeSj995e2j7BAAwzGuh737nd5L3/8GrvgQAKIhAAwCMuRcbSfLF+tfD3g0AgBG6FtrwbQBAgQQaAGBMRa+FvBZ9cbMNADDuXAsB/aTnOGwl0AAAYypvaKQYNkDvBgCgClwLAf2k5zhsJdAAAAAAANAGvaUgn0ADAAAAAEAb9JaCfAINAAAAAAA9Mm8DVSbQAAAAQGn96INPkq8eP/nr/20kz549Sx78p+dD3isAqsi8DVSZQAMAAAClFUGGL9a/HvZuUDJfPn6S/FcffZ5MTf06SZKJzXHX84ZEAWjFvA0g0AAAAMAYmJxIktd2vZT2aJiaijZ1E7kVP5C1Ov7NX0XPF71fgN6ZtwEEGgAAABgDr+3amfyP5/+L5M6dO8mBAweSHTt2DHuXGEHfBp82NoNS0bshAg8AQPcMnQQAAABUqtXx8+fPN4NSf2/hfzD8FrScA+fbyY4B8gg0AAAAAABbmAMH6IRAAwAAAKVoSZtH61qKEufSW+9db7meSaOp5hw4O7csMwcOUE+gAQAqKO9G2k0zAKNES1oGLeZp+GLdsDBQL4IMn757zIEBmhJoAIAKciMNQJlb0ubRupZutXvuREMNk0YDQD6BBgCo+I20m2YAxqUlbUzwC91OEN1K9AbV4wEA8gk0AEDFb6TdNAMAAAC9mOzp3QAAAAAAQKXp0QAAAAAA0AcxVG30Iq8f0rbdYdugLAQaAAAAAAD6ICaRN78LVSDQAAA0bGmTR+sbAACA1vdNefdcEXiAcSTQAACktLQBAAAoRt7QSNGwS+8GxpVAAwBUXF5Lmzxa3wAAAAB5BBoAoOLanYRM6xsAAIDemSCacSTQAAAAAAAwxGFr4//1c+aZH48yEWgAADqi9Q0AAEAxw9bWBhwEHygzgQYAoCNa3wBQtB998Eny1eMn2wLbADDuw9bmlYHNgg8wqgQaAICBtb5plG6780QAMJ6igkVFCgBV1G7wIQLw0egLRpVAAwAw1NY3xiIFIDM5kSSv7drZMtANw2QYSWAY917RiEtQnlEm0AAAdK3d4EMeY5ECUC+CDJ++e8yBoXTDSAJA1Qk0AACFancYpF56Q7QbzMgYngkA6FVe75q8oUyaXae4JgF6pVcVo0qgAQAoxViktRfUWhECAKM6lIk5R4Bh9qoS7KTSgYbFxcXk/v37yfT0dLK2tpacOHEiOXnyZN/S6Ne63fgP/+E/JMnf2FNYeoy+gwcPJk+ePElu37497F1hQPzO+0PZsadyN/CNLqi/t3vrWN6tghTttigcpV4Tyo5q8r0P1yiXM7VOnz6d3LlzJxlVg67s8LtxXIZ1rtRea8Tf9XOO6PnAsMuZRuSb49erqlmws9P58ZwflCrQMDc3l8zMzCQLCwubyyJjfPDgQXLu3LnC0+jXugAMjrKjWppNAtqqkqpZkCLvIrv+9U60Sq/dfW7XuLRUGpfPwXgZ5XKmbNqt7MgqZqGs8q41auccya5J9NBkGOUM498oq1Wws90haqEXExsbG3Vxr8FZXV1NDh06lNTvQqPlvabRr3Uz6+vrycsvv5z89re/TXbv3t3GEUiS1/+v/89k42/sSb63+6Xk03ePt/Ueyu1v/a2/tdlyZseOHcPeHQagm995N/lJVYxb2dHN963saF+reSDa1azXRNHp5fkP/+GLJE6f733vfxeXb21tt51t1Fbk91Lh30vvj7whJ2o/wyAnhR21oEen1wzKjvEvZ/K+79dffz3t0TDo68p2f/dZS8tGlR1F/+5dazsugz5X2i07mpV3teV29psZVBmo7Bj/cqbV993st5Cdt4O+JqN97V5LNxuittl6tfchr+7aqREQTcuOofZoWFpaSmZnZ7ctz5ZdvXq1ZZevTtLo17oADI6yg37MA9FMUZXsPbUi+uuhsr5Yf9K0cqLTbdS2Jm62fj97fzQbXiJvyKt2tRtEqd+vdva5l220u88MzyiXM6Ok0zHoa1t2t6qYhTJpN99ut4dmq4AE5Tfocobx1ipv6bQ3xLb8Z8t9yDcNnqCRoQYalpeXczPGsGfPnuTatWstM8ZO0ujXugAMjrKDXhVdkdtuer1UROf1aCiqQr3ZvBfdBEfa7f3RzvASrSa6KyKI0kg7x6DXbbQKZvzV8XeTjf+03nO6jE8504t2W1130lOhPkjYTKvKDhh3nZ73RQXbGT2DLmcYb0X81pvNhZfeh+zcnSQTkyORL41aD2RGKNBw79695Pjx/GFE9u7dm9y8ebPQNPq1bibrmhbd1Oq7m/1P/9P/lJvWiyf/KdmY/J3k2dfPkocPH+auw3iK79vQSdXQ6nf+d/7O39m27Pnz5+nzEEe3G1njVnZ0U34oO8rp//Ff/a2u3/v222+nXdr/zaef5pYdWd7S7jb+jx/+f5L7f5l/gT79nZeSPz4723K9Vu/tdLsvT/3/Nj9H/P3sd54l3fjy8dPNv//iq7/a9vpru36nZRqtjkER22grvcnfSZLJqdxrBmVHNcuZvLLjz//8z5P//D//z5OJiYmmZcf/96v7W861Wn/xVZIc/r/96/TvRus08t1dv5P8v/9x4999rUHe77jWdlzKeq48+/ovkxdPniYv4rf59V91mcbvuO8YYYMuZ1rdd3z22Wfp8rzfQnY+qrMaT82uubPr4bgPefz3/9u0Z0Mv+VLttUYvml2nFLWNKvpHb/3vk//T3/k/5L7WyX3H0AINjx49avp6RGBbrdNJGv1at9bjx4/T57t37+Zm9M38eazzf2+6CmPm1VdfHfYuMGDd/M4jX4lx7xjfsqOX8kPZUT2DKjv2ni/fe9sVv5t+70u722jXq0udTeio7BjPcqZR2fHs2bPkf/1f/9eO7j2KPG9HtSxyre24VPlccd8xugZdznRy39HstzCqeT39tXk9/KtiJhcv+hp5WNsYRytJkvy3Xbyv/r5jqD0axs3v/u7vpi2Kdu3ataVFEUCnIiocGXbkK4w/5QdQBGVHtSg7gCIoO6pF2QH0s+wYWqAhIqy9RGg7TaNf69aanJxMfu/3fq/pewHapSdDNcqOoPwAiqLsGN9ypp6yAyiKsmN8y5l6yg6gn2XHZDKiHjx40DLzLDKNfq0LwOAoOwCoajkDQPkNupwBKNJQAw2R8UUG2CgCe/jw4ULT6Ne6AAyOsgOAqpYzAJTfoMsZgEoEGk6fPp3cu3ev4esnTpwoNI1+rQvA4Cg7AKhqOQNA+Q26nAEYlImNmL1hSJaXl9PM7+HDh1u6dTVa3msa/Vq3W4uLi8n9+/eT6enpZG1tLU335MmTPaXJaIrvdnZ2Njlz5kz6HBcES0tLaUuDeKb84rs8e/Zs+h03+x373feuymWH86d65UEn37nzY3zLA+dBtcuZ//K//C+Tf/SP/pHri7/+fV28eDF9jvwzWvNeuHBBvlhzbCJPibxldXU1PX/m5+dzz9GqlhmHDh1KVlZWcl+r6jFh/MuZZpz3qMugUBtDdvLkyY35+fkty44fP75x5cqVbevu27cvffSSRr/W7dS5c+c2FhYWtqW9tLTUc9qMntnZ2QjobXnE9035RT4Rv+d4xPfaLH/wuy/2uFet7HD+VK886OQ7d36Mb3ngPKhmORP//8/+s/8sXd/1xTcePnyY/h7iObOyspIenzheVc4Xs2NTv2zPnj2VPza1svw5T5XOF6pZzuRx3lebugz6Yag9GjqNoEbrg7179ybXrl3rOo1+rtuuaF0Sn6X+0DdaTvmdOnUqOXLkSHLjxo1k37596Xl0/PjxYe8WBYpWdTMzM8mVK1dy8wi/++JVqexw/lSvPOjkO3d+jG954DwYrlEoZ6Lnk+uLb5w/fz7tvVDfQjeOW7wWxz/LT6uWL166dCmZm5vbcgyyMufq1avp+RRlThWPTe3ni95m8dzuZx/3Y0L1yplaznsy6jIoVF/CF7SMGkeLxjytWixRTvUtjBg/a2trTX+/fvf0wvlTvfKgk+/c+TG+5YHzANcXrVv0Zr0aavPXquWLcQyi98K1a9dyW/DH61U9Nplo9R09FvKqQKp6TKg25z0Z1xoUaaiTQVdVjJmXtSipFy108qLUQLn53eP8oV95hvylXPr13ToPqqlK33t8zpiToV7Ww6H2tSodlxA9X2I89voectlxiNfrl1Xl2GQtvqM3TCNVPCbgvKddzhU6IdAwpG5J0e0tTyy/efPmwPeJwYjuidG1OZ6pFr97nD90Uh50kmfIX8qlX9+t86CaqvS9R2VvVKbXy/LRGJauiselkRgyKYIvMYxbraodmzg/IojQbFLcqh0TCM572uVcoRMCDUOYzb2ZuABqtQ7lExf5MXZsPJ87dy59jnEWI8Nm/Pnd4/yhk/KgkzxD/lIu/fpunQfV5Hv/xtLSUvp7iDy16sclypJovR9zM8Rx+eyzz7b0Zqjisbl8+XLT8eqreEzAeU+7nCt0aqrjdwAdiwmZspufEN2az5w5ky6PCZsAqAblAUCxwznEI1rtN2uxXhXRcn9+fj6tGIpecxHIjmNTG2yoklZDJgEAxdKjYcBaXQBrLTGeaoMMtcGGaHUU3ZoZb373OH/opDzoJM+Qv5RLv75b50E1+d6TzZb7tS3WHZdvjkEEHKJ8qe01V6Vj086QSVU7JpBx3tMu5wqdEmgYMTGEgtY41ZBNOGZyMfzu6YXzp3rlQSffufOjXPr13ToPqmncv/cIMkRr9bwAbpWPS/0xCjFkX9WOTashk9o1TscE2uW8x7lCtwydNARxoRIZd6MWE4cPHx74PtHfC/xoRbSyspL7eqNzgfHid4/zh07Kg07yDPlLufTru3UeVFNVv/eoOI/Jn6PVfp6qHZfotRCTFtcHrLMJjrMJs6tybKKHYHzmubm5LcuzSZ2z5QsLC+nxqMIxgXrOe9rlXKETAg1DcPr06c2LnEbjN5OMXbfdetnFbNwkMf787nH+0El50EmeIX8pl359t86Daqri9x5zD0xPT28LMsTyrHdD1Y5Lq/Kldo6GKhyb6MmQ15shAgxxrGK4rVpVOCZQz3lPu5wrdMLQSUNq0RgXOPXjPcZEZiHG0mR8xEVu3nAY2VjcnXb3ppz87nH+0El50EmeIX8pl359t86Daqra9x6fKz5rXk+G2mNQteMSnyevt1xW5pw5c6ayx6YdjglV5LzHuUJfbDAUJ0+e3Jifn9+y7Pjx4xtXrlzxjYyZhw8fbpw7d27LspWVlY09e/b4vsdIfKeRpS4tLTVcx++eXjh/qlcedPKdOz/GtzxwHlSb64tvra2tbezbty/NR+P3E4/4Ox55v58q5YtxbOLzRjlTuyzyovpyp2rHpv5zN6oCqeoxodqc9wTXGhRpIv7pTwiDVhYXF5P79++nXX/X1tbSLplFTFjF6IkWQzGWbIxtF+Nzh5i8rrYbM+UU32t8p9EyLJ7jO46WYDEmbn237OB3Ty+cP9UrDzr5zp0f41seOA+qx/XFdjMzM5v5Zp5o0V+fl1YpX6wtX+LvOFYxVJAy45thtaJ3R9aDMM6TmHehPm+u0vkCGed9dbnWoB8EGgAAAAAAgK6ZowEAAAAAAOiaQAMAAAAAANA1gQYAAAAAAKBrAg0AAAAAAEDXBBoAAAAAAICuCTQAAAAAAABdE2gAAAAAAAC6JtAAAAAAAAB0TaABAAAAAADomkADAAAAAADQNYEGAAAAAACgawINAAAAAABA1wQaAAAAAACArgk0AAAAAAAAXRNoAAAAAAAAuibQAH20uLiYHDp0KJmYmEgf58+f33xtZmYmXRbPc3NzvgcAlB0AuPcAYGDUW1GkiY2NjY1CUwS2iUDCpUuXkocPHyZ79uzZzMzX1taSpaUlRwwAZQcAhXDvAYCyg2HQowEGIIIJs7OzyalTp9L/P3r0KLl//74gAwDKDgDcewAwVOqtKIIeDTAg9+7dS4dJisw7ejIsLCw49gAoOwBw7wHA0Km3old6NMCA7Nu3Lw0yRFfmI0eOOO4AKDsAcO8BwEhQb0Wv9GiAAVpdXd0cPil6NQCAsgMA9x4AjAL1VvRCjwYYoMuXLycrKyvJgwcP0p4NAKDsAMC9BwCjQL0VvRBogAE5f/58cuHChWTPnj3J9evXk0uXLiXLy8uOPwDKDgDcewAwVOqt6JVAAwwosw4RZAizs7PJuXPn0mGUYrIdAFB2AODeA4BhUG9FEczRAH3OqKPnwqNHj9Igw8OHD9PlEVw4dOjQ5vLjx48nH3744WYgAoDqUnYAoPwAwL0HZSPQAAAAAAAAdM3QSQAAAAAAQNcEGgAAAAAAgK4JNAAAAAAAAF0TaAAAAAAAALom0AAAAAAAAHRNoAEAAAAAAOiaQAMAAAAAANA1gQYAAAAAAKBrAg0AAAAAAEDXBBoAAAAAAICuCTQAAAAAAABdE2gAAAAAAAC6JtAAAAAAAAB0TaABAAAAAADomkADAAAAAADQNYEGAAAAAACgawINAAAAAABA1wQaAAAAAACArgk0AAAAAAAAXRNoAAAAAAAAuibQAAAAAAAAdE2gAQAAAAAA6JpAAwAAAAAA0DWBBgAAAAAAoGsCDQAAAAAAQNcEGgAAAAAAgK4JNMCQ3Lt3z7EfQaurq8PeBYC+Uv50R/kADIt8uzn5MwCMBoEGGNLF8NLSkmM/os6fPz/sXQDoC+VPb5QPwKDJt9sjfwaA4ZvY2NjYGPZO0LtTp06lLV2y1hwnT57c8vqjR4+Sffv2JQsLC8mePXvG6pDH584+f3zGlZWVZJTFdxH7e+3atZbrXrp0Kdm7d2/y4MGDze+w/rstq/g8Fy9eTJ/ju4vPeOHChY4+X5bG9PR0cv/+/fT8P3HiRDI/P79t3Vg+OzubnDlzJn2ObUawJ9KoD/osLi6mv5Nz584V8lmh7GXKw4cPW5Ydy8vL6e8sHD9+PP2dRZkzjuVLs3In8pRDhw4lc3NzuXlRWcqfcSk3ohxdW1vb/P/MzMy2vF35AOOtnfukKOMiP4myq0z5drx+9uzZ9Pq2iHuEKL+Kupfq5Dq9lXbyctfvADACItDA+IivdHZ2Nve1c+fOpa9fu3ZtYxwdP3684Wcftf1cW1trus7Dhw83FhYW0uda8flWVlY2yi4+V5yPtZ8vPlecnydPnuwojfple/bsyU0jjl2kX/uI76KRWL/V9wTjLn5L8VuJ/KiVpaWl9PdXhny4qPIlb3nkG53kZaNW/oxLuRHrxeeN87J2/X379m0rR5UPUA3Nrv2inIvXa/OMbszPz28MIt+OfC/yxOz+7sqVKz1vK0urCJ1epzdLp928PLh+B4DhMnRShUSr7Wh5Ga1iooXJuInPNuqixW87+xqtf6K1T30L4ujdcPPmzaTs4vPV967JWj9fvXp18zg189FHH6Wtm2rXjfSiJXWkUT+WbdajJ1p7xbGNlmHNWvVGqzZdsKm6+N3Eb6rVUG9RpkT+FLLncdIoz85bHsuiPuvKlStJGcufcSk3ooVvfNbaFq/RAyIe9WWr8gGIa8PIU6I3Wi/zIRQ5l0KzfDvKmCibi7pWjd4GRd5jdHqd3kgneXlw/Q4AwyXQUDFxcRcVQu1U5FK8dm4I4uI7uj/Xy763+A7LLj5jdM2ul322diroDh8+nHuDkVVy1gfTYnncREbaUTHV6jhGQCKO9zgG5aATWaVLs3IjKhTGZVi3cVVkhdSolxtRYRbr13/eWDeGAauvtFM+ALX5Sbf3SZHvFHndOMh8+/Lly7n3H93q9Do9T6d5eXD9DgDDJdAAA9ROoODGjRu548NGi56oIC9ra9Ra8RmiJVK97IYk77V6cYziJqP+eMYxjvSLGGM30o7WWFBlcdMev00T2Jdb2QPVnZQb0fshlvezvFQ+wPjJWtl3k3fEe+NavYz5dsxtED0BilTEdXq3ebn8GQCGR6ChYrIWOvUXfXGBmT2i1Uht5Wr8Ha0IJyYmkldeeSVdJ5vMK1uWtTSJ5fH/WB6TdGWtVeI5WsVm6cfwTdmEbNl+xTaytOK1WDe2ES1ZehEX/pFmtu28lkFFfMZW4nO0c6EcaYY4RvGI/YpHt5OnjaIYsihuPupl58SRI0e6SjeOcVQ2NesREduI77v2/GskjnlZJ02FIp0+fbphS834LbVTEdJOXhyvZXlxPGrXibwxy3ejPBlE2dSp+Iyx/dhWXuv7do5BbD8rD7OhDmP9+Lzx/25at7Zb/sR2suOc9xhmD69Oyo1sku7s2OVd2+RRPkB1Rb4ReWWUZ/VlWquyorbVfQw/lF3D1+fX7ZQBtWkOonFRfI7YTl7vg6K1c51eq9u83PU7AAzRkOeIYICTQWeTnNVPBh3r108gFhOZ1U7UFRNxxXvrJziL99ZP9BUTltVOtBb/j4m/aifsylsWYmKvSC+beDT+3+6EYfG+vM9ev3/xGfLW6+UztiPSbTVBXHw32TGJv2Pis3hu573jII5pnBf1k2A3E99FnC9xnsT7G703Xo9jmJ3/8dxqwrhsolGoqizfySY4zpsks3ZZ/H4b5Y3t5sXZurG92t9z/M7zJpYcRNnUrHxpd5LoTo9BpBGP+om4Y//amZy7fjutypBse3E9EPljlh/HsYjHqE4inVduZBO+1h+nbPLUesoHqO5k0JF3RL4X9xx5+WQnZUWWjzbSSRnQ7rV/Vj53Oxl07Taye8UitXudnqfTvDzj+h0AhkcN2piJC7K48M0uTrNHdkFWX1EQyxpdEEc6tRetkUb9xXC8P9ZrVFkeIv28YEF2wdmswqCTi9G8Cp/Yl/rgSlYxlXdB3u1nbEek3apyqNnrsQ/tVvTEdx43TJ0+ur1JKUL2XXW7D/G9xvGLz5H33eRVkGbrN0uzvrITqqS2AiLyxrwKkdrfbKNAQ6d5cba9LK1Yt1GFyyDKpk4DDUWVR/Fafb7fKIjRS/nTqHIsr3HCKJU9eeVGdkzzjl1WIVf/mZQPUA3x+488J/LD2kfk2fHIy4s6vY9pdF/VTRnQTlC5l0BDpF8f0O9XA5tW1+l563eal9e/1/U7AAze1DB7U9Af0cU0xvJvR3Q9bTTudnQbju6p2QSfMUFYdAPOurFmojtrdLvNxtqMv7MhfrIJjPP2J7q15nUZru2+22s33kgrPkft/mZpZuOw1urmM7Yr0ujl88T+RJfjdrYbx7vdc2BUxHGPc7HbCWXj2MaxWVtbS4cdiefa7/DcuXMNz/E4rnnbrR37exBdymGUxbAR2cTQ2W+r0W+n17w4xNAKMYxPlFPxe26Upw2qbOpVp8cgJszMG5u60VwF3ZY/cRzikTcsUWyr3aGxhlH25JUb2eeMfc87diHeU/uZlA9QHZEP5F1LRz4Z148xL8DKykrhZUWnZUCv9w2jNGRSO9fpeet3mpfXv9f1OwAMnkBDhWVji0ZlRp5YXnvhm93IZ5Xd8f6scicu9PICFjFOaYiLybzxNPMu3PP2p1FFVKuxS+P1bIz9uGCPdLK07t+/v239bj5juyKNRse6HbH/efs8DFH510lFV3zu+B4afV9xjGMSurzKnk5FWnGuxc1fqzFgs/2JfWtWWTrMcclhVMTvMwINkW9neWG7+UCneXH2nthObLPZb3lQZVOvujkGvZQZ7ZY/WYOCvMqm+uDNKGlVbjSqPIvlja4paikfoFoib4jyIgIIUe7E30WWFZ2WAb3eN7Ry+fLloTRK6uQ6vde83PU7AAyeQEOFddoiMkRlRFx4R2VOtPCJ52xis1gey2pblmQXyHHR3m5L9bwLykYTdH7TC7q5qHyKfYtWrdHytVGFSrefsV2xzWbHvLblbZ7Yh24nSS5a3HAVJY5rfK5Oe4hEa6gsgFErO+dqz5es0jFrodbpb2FUK9pg0CJ//Oijj9K8MH5Thw8fbvu9nebFIdKP319t77phlU1F6PQYFNXStFn5E3llVNjnLQ/dlHf91qrcyCYPzVO/XPkAZLL8Luvx3WtZUd8roZMyoNV9Qy9iPyKPj4BKrSywki2PQES35VAn1+mNdJKXN3o/ADBYAg0VllUQNWoNEhe39RdoceEZF9vxnqz1zenTp9PlUZETF5O1rWOyivN2Wg8206iCuJWsxUy8v92LzU4/Y7uaXSyHSLtRy8z4HHGh326FT3zmSK9TUdk0yEq3+FzT09PbKotieaveDVmX73rZTVlt0KbVuo0CONn3Zdgk+Ebkg1FBEY/4/bTbC6mbvDhrcRnvef311zdbmOYZVtnUiW6PQREalT/Z588LcmeVYZ3kf4Moe9opN6KsjIBYI7UBMuUD0EivZUVU3GdD/3RaBrS6b+hF5Ld5eW6UmZEn9tKDu5vr9EY6yctruX4HgOGZHOK2GbKoPIiLzKjIyRMVSfVjj8YFX7wvq9TJ0onlcVEaN//1ojKgWeVQP8X+R2VT/YVu7YX74uJiz5+xHdGCp9nQRzdu3GjYcikqyKJ7cbsVPrF+3Mh0+hhkkCEqo+J7aDRObivxfeQFoLKWU9FSLBOfq75FVXaOh0aVpXFjKcgAW3932ZBG/c6L4z1RAR2/wevXr6eVNI0qsctQNnVzDIrSqPxpVNkV+xSVO+0MazHIsqfdciOOdTY0Sa2sBW3t96t8AOqvC2vzpE7KivprxsiHshb8nZYBre4bRl0n1+mNdJKX13L9DgDDI9Awhjpp/fLhhx+m69ePOxoXbnGBmFcBGxfJ0UKntnV9DD3QaFLQrNttfQVKrB/v62X/W723fp6JbLux78220+lnbLcLcbOWntGqvv5GJvYxm2xuFIev6FZ8J3GOxRBM8dnikU00G5+3/iYs5oSIR604VmfPnt3yPUa6cZ7FeVv7PUVlZf3NSNykxIR/zQI4cQ50MjQMjJP4PeUNbZD1Eoh8sl78HvMCpp3mxVmQO/ttRsvH+F1nw9zkGVbZ1O7yTo9BHMdOlndb/sRxzSrXsv0+duxYen0wSsNOdFJuZJO91uf7UWbE8tqWtMoHIER5F3lE5B+R/3VTVkRelA0/lOVbWTnWaRnQ6r4hk7232TBLedfReZql8corr7SVRjfX6Y32sZO8vJbrdwAYnomNdga5Z+TFBVhcVNWOqRwXX+0O8ZNV6kSrz7iJj4vbRq28sy619ZXicXGd12o8bxu1+1hb6Vvbkqjdcfuzcbizlo7x3ixQEv/PLnKzVq7ZhW2sE9vPWswW8Rlb7WdcQOf95LIJ4eImJHqYZBfa8V3k7V/Zxedr1g09WkDV3jw0Guc1vtf47uP4ZC2e4nvNq1SsXTfbdhzbZl23I61mvwUYV1F5kuWpcaMfv4UsP87y1drW7vF6/K6ySpH6PLzdvDjKgQh8Z+NaP3z4MF0n0o7fYrY88veoCKrNG/tRNjUqX7J5I9pZnuUz7R6DqOhpJ+3Yz0b5XSflT3YM4vPH8YxyJ9temcuNEOdSfPdRfsRxzcZEr6d8gPGWBamz+6TIB7KeBpE3ZHlL5BGN7j2alRX168V2Ip+P17Pgd6f3JO3k29lnygIasa34XPXlYKPr6Pq8MrsPi32Jhja16UQa2Xxn7ZYPnVynN9vHdvPyjOt3ABgegQYYoLiIjuBPfe+ErEXTuAUUyi5u8Iqc+Bpg1MofuqN8AKqWb2eBmnbmVxgm+TMADI+hk2CAoqVS3jivxhIdPRH8GfUbKYBeyx86p3wAqphvZ5NbjzL5MwAMl0ADDFB0882b1IzREzd27Q49BjDqlD/FUT4AVcy3s+ELR5n8GQCGayopmZhA6v79+5tzCcQYl51MzhsXSDEGdHbRFuM8RmuRRmn0uj3IuwCOsUOzMUjLcNFeNTEWbN6k1ADjVP7QOeUDUMV8O2/C6lEjfwaA4SvVHA1xkRVjLtZO0hWVgTHBVzuTtWYTUkUr5axiN8aajPEvI3hQO7FlEduDRmIyz7hhiHMx/o7JzQzTMxoiT4gJufVmAMa9/KEzygegqvl2NL5rNFH2KJA/A8BoKE2gIQsI1O9uo+V5IsgQvRfqW4/HhVO8Fhdw2WRbRWwPWp3Tggujx/cCjDv5nOMGlIt82/EBgDIoTaAhehfEBFQrKyvbXpuYmEh7I7Qa0ih6J4QYAikveBC9FLIJt4rYHgAAAAAAjLupMnUZbdT6O3ooRG+EVhX/Md56BA/y3h9ivoZetvfixYvkL/7iL5Jdu3alwQiAbkUM+PHjx8nv/u7vJpOTkw7kmFN+AEVQdgxep/O/5YmhWeO+48yZM+lzpBONnyLNrBFUHmUHUARlBwCVCzTEBXc2rFG9GN8+L4BQr9EkWtGjIRw5cqSn7UWQ4W/+zb/Zcj8A2vXnf/7nye/93u85YGNO+QEUSdkxGJ3O/9ZIBCdiKNd4ZOI+pNUEwMoOoEjKDgAqEWiIi/hm4sK+1TrNREuhSCOb4Lnb7UVPhvD7v//721ogf/TRR13vH+UWrc0ePnyYvPLKK1qmk3tu/ON//I9zz5v/7X/73zbzFcZb9j3HDd7u3bu3vf78+fN02L8YAnDHjh1D2EPq+U5Gj+8kSdbX19NGL8qOwYieDLVBhhA9EmJZBCCih3Sjhkv1va6jN8ONGzfSv6OHQzvva1V21PMbwblBHmUHAJUKNPRT3ADEI1oc1U8S3alsuKRf/OIXbV3sUw1xU3f37t1k//79KgjJPTf+l//lf9l2bsQFf/SeMgxbNWTfc5QdjQIN3/nOd9LXBBpGg+9k9PhOvqXsGIyrV6+mj/r537IgQdxftBMwiPJ+fn6+8LKjnt8Izg3ayVMAYKwDDa0CAL30Zjh16lTao6F2HNVetxeVQCqCqBU9XJwXdHJuyEMAYLR1Mv8bAACMu1IEGlqJi/i40O8myBCTtWVDJvV7ewAAwHjoZP63dsT7InBx+PDhdAimdh09enRbA4Vbt27l9miIoRnjGZwb1XTw4MFty+QJAFQu0BAtgxq1CooeBnFB3okYNzUu/ht1Uy56ewAAwPirn/+tlbjniHuTmJsh3hPDusaE0jH0UjuNm54+fbptHrA7d+5sWy+CDLGtGLaxfn2qzblRHU+ePMn9/gGgUoGG06dP53ZNzsSFebsuXbqUTE9PbwsyxPLshqDI7QEAAOOvm/nfsgBDJuZ1iMmhY3n9/A95VlZW2p6jwbxhODeq7fbt29uWZXPDAUBlAg0xzFEEAqI3Qe1Fe1zIh3YmWsvWjzTyejLUzr1Q1PYAAIBqyJv/rZW8ng9xrxG9HGKy6VZpdTIPmHnDcG6Ql4cAQBFK02c2LrbjIvvixYtbli8sLOS2GJqZmUkfte7du5fMzc2lLYPiwj0e8f94RIuh2q7JnW4PAACorm7nf8uT3Zc0mgcCAABGTWl6NISo4F9cXEwDBDH0UQQMIkiQ18onAgH13f8imBDBhuipkCeCCN1uDwAAqKZW8781C07E/UkMf5Sn0ZxxAAAwakoVaAjtXrznXay3M8Zpt9sDAACqp5353xpZXV3NnfA5CzBE8AIAAMqgNEMnAQAAjJJ2539rJHpK5w2PFHMzhCKGYQIAgEEoXY8GAACAYcvmf8smbq4NLmSv1crmj6vtZR1zOsR6MYF0bS+HmCfOvHAAAJSJQMOQ/eiDT5KvHj/Jfe3VXS8lH7/z9sD3CQBGnfITGLZO53/Lm0MulsV6EWyIvyO9cP369WR2djYZZfJhAABqCTQMWQQZvlj/eti7AQClovwEhq3T+d8aTfgcAYbaHg1lIR8GAKCWQMOImJxIktd27Uz//vLx18mLjWHvEQCUs/yM57feu75lPb0EAQaXDwMAUD0CDSMiLs4/ffdY+ndUjujlAADdlZ9RyVVfjsb/BR8Aiuc+BgCAINAAAIyF6LVQrzbgIIgPAAAA/SHQAACMhY/febutyUoN7QEAAADFEmgAACoVfDBEIQAAABRrsuD0AAAAAACAChFoAAAAAAAAuibQAAAAAAAAdE2gAQAAAAAA6JrJoAGAkfajDz5Jvnr8ZMuyLx9/PbT9AQAAALYSaAAARloEGb5YF1gAKIsIBr/13vUty17d9VLy8TtvD22fAADoL4EGAKAUJieS5LVdO7dVXAEwWl5sJALEAAAVI9AAAJRCBBk+fffYsHcDgAbygr/RuyECDwAAjDeBBgAAAHqWNzRSDKFk+DsAgPE3OewdAAAAAAAAykuPBgAAAMbKjz74JPnq8ZPc10xMDQBQPIEGAAAAxiqoYLgmAIDBEmgAAACgtCLI0Cyw8L3dO9NnE1MDAPSPQAMAAAClNzmRJK/t+iaokDdMkompAQD6R6ABAACA0s+pEEGGT9891oe9AwCgFYEGAKCSYgiNaN1aywShAKM5/BEAAKNNoAEAqKQXGyYLBRiV3gsR/K0f/iibU0FgGABg9Ak0AACVEr0W6pkgFGA0ei/UDn+UzamQFxiO/2e90rIgBQAAwyPQAABUSt443+1OEPptK9yN5NmzZ8nU1K+TJJlIXzPsEkAxkzfn/Z2pzasNtQQAMDoEGgAYGYuLi8n9+/eT6enpZG1tLTlx4kRy8uTJvqXR7rqPHj1KLl68mD7fu3cvefDgQXLhwoWe06Wz4TRGsxXu8yHuDUC5tZq8OS8w3GrSaAAAhkOgAYCRMDc3l8zMzCQLCwuby6KCPir1z507V3ga7a4bwYXz58+n6+3Zsyddtrq6mhw6dCgNHly5cqXwz1FlZZkMNFrh7v0bO5Kpqanky8dP0mE9AOi/vOADAADDJ9AAwNBFxf2lS5eSjY2ttbVRWR8V+u1U0HeSRifrRk+G2iBDmJ2dTZdFAGJ5eTk5fvx4YZ+D9obTGLbXdr2U/Hc//r3kwIEDyd9b+B9KERwBAACAfpnsW8oA0KalpaW08r5etuzq1auFptHJuvF3BAnqZcGF2h4NRXwOtg6nUfvQihUAAABGk0ADAEMXvQL27duX+1r0JLh27VqhaXSybqwXwx7lrRdqXyvicwAAAACUjaGTABi6mGA56yFQb+/evcnNmzcLTaOTdRsFB2KYpHDkyJHCPsfRo0eTHTt2bFl269at5Pnz58mLFy/S5/GXDTu1MeDPu7E58fRb7y1veeW733kp+df/5O9uWS/+/fY7GdY+U6tav5MkOXjw4LZlVfnsAADA6BFoAGCoYrLlZqInQKt1OkmjiO1lwyTFutm8C0Wk+/Tp02Rycmtnwzt37qSVp9Fz4u7du9teHzfPnj3bfI7PPujtxqTOX6w/2fZati+1+5d9J8PaZ7aq0u8kPHmy9TzNjgEAAMAwCDQAQIdiiKR4xPwMtZNE92plZSXZvXt3bivlqDzdv3//th4P42Zq6tfxiZOpqal0ouVB+d6er5Kpqa0Vt18+fpIGHmr3pXb/opdKfCdTU38xlH2mur+TcPv27W3L1tfX0/MSAABg0AQaABiqVhX17fQu6CSNIrZ36tSptEfDyZMnC003KkcbVZBGC+1mr4+Pic3nQX7Wf/NP//62ZW+9dz35Yv3run2Z2Pw3+06Gtc9U+XeSr6qfGwAAGL7x71cOQKnFUCi99hroJI1W60aQ4cKFC5tDJvVjHwAAAADKRKABgKGLCvioiG/UE+Dw4cOFptHt9s6fP59O/jw/P9+3zwEAAABQNgINAAzd6dOnk3v37jV8/cSJE4Wm0c32Ll26lExPT28LMsTyXtIFAAAAKDuBBgCGLoYjWl1d3TaPQUy4HI4fP15oGp1uL5bHunk9GWrTKOJzAAAAAJSNQAMAQxcV8DGx8sWLF7csX1hYSK5cubJtboOZmZn00W0anawbPRTm5uaStbW1dOikeMT/4xE9FPbt29f15wAAAAAYB1PD3gEACFERv7i4mFbkxxBFUbEflflRcV8vKuz37t3bUxrtrhvBhAg21A6RVB9E6HYfAAAAAMaBQAMAI6PRJMv1VlZWek6j3XUjUNCpTvYBAKrgy8dfJ2+9d33zbwAAxotAAwAAAH31YiNJvlgXYAAAGFcCDQAAAF149OhROi9PPMcwew8ePEguXLjQ8XB5MeTe/fv3N4fci2H7xmXIvVd3vdTVawAAlItAAwAAQIciuBDz8cRcPTF3UFhdXU0OHTqUBglizp52xDw+MzMzW+b8iUBDBC3OnTtX+u/l43feHvYuAAAwAJOD2AgAAMA4iZ4MtUGGMDs7my67evVqsry83DKNCExcunRp29w+kUYEIAAAoCwEGgAAADoUwYTovVDv+PHj6XM7PRqWlpbS4ES9bFlsAwAAykCgAQAAoEP79u1Lhzeql/VwyHutXvR6iHTyRDrXrl3zvQAAUArmaAAAAOhQoyBADIcUjhw50jKNmEA66wFRb+/evcnNmzdbpnH06NFkx44dW5bdunVr23rPnz9PXrx4kT4XY2Pzubg0+62M+9x/xZ8bjKqDBw9uW+Z7B6AoAg0AAAAFieGQojdCq4mcYzLpZiKNVuuEp0+fJpOTWzuq37lzZ9t6UZEcvSzu3r27bf1uPHv2bPM5b3ujqIz7PAhFnxuMridPnuR+/wBQBIEGAACAAsRQSPGI+RlqJ4nup5WVlWT37t0t14tWy1GRvH///m09ILoxNfXrSDWZmppKDhw4kJRBGfd5EIo+Nxhdt2/f3rZsfX097UEFAL0SaAAAACjAqVOn0h4NJ0+ebLluq0BEO70ZQlQMt1s5HK3VO1k/86MPPkm+ery1JfSXm/+fKFHl9EQJ93kwuj03KD/fOQBFEWgAAAAoIMhw4cKFlkMmtSuGsmk0UfSgRZDhi/Wvh70bAACMMIEGAIAmvnz8dfLWe9c3/waod/78+XTy5/n5+Y4OTvRqiIBCox4Nhw8fHqmDPTmRJK/t2rll2au7Xhra/gAAMDpKF2hYXFxM7t+/n0xPTydra2vJiRMn2uqanHfhfvbs2eTMmTMN3x9pz87OpuvE871799Ku0PHeeAYAxt+LjURLXqChS5cupfcm9UGGWN6qd8Pp06eTmzdvNnw97kdGSQQZPn332LB3AwCAEVSqQMPc3FwyMzOTLCwsbLn4jlZA7XZRji7N2URHV69eTYMIjUS6EdiIR+b48ePJtWvXevocAMDoa9ZK97vf+Z2B7gswmmLi52iElNeToZ05FuLeJAISsW7tnA2RbnbvAQAAZVCaQMPq6mp6Eb6xsbFleQQdDh061Hag4cqVK+lz9E6I9JqJMVEjEHHjxo307whquNgHgGLkTy46OkMTffzO2w1fe/78eXLnzp2B7g8wWuJ+IhpCxf1BDJ1UG1zIXqsVDaZC9MrOxHujd/XFixe3NKaKv+O+pdWE0QAAMCpKE2iIoYpi+KJ62bLondDNEErNRM+HTsdZBQDaY3JRoMyiEVKzxku1gYMQQYOsZ3WtCChED+oIVmTDw0aQouh7G4oPjtf2gGsWnAYAqILSBBqi+3BeoCG7aI/hjFyMA0D5mFwUKKPangntWFlZafiaxk2jH1T4Yn10etwBAIyi0gQaorVQo2GLomVQs0nUihi2KdI/fPhww2BHraNHjyY7duzYsuzWrVsN1s6GgtpIh2FovIyyiu/wxYsXvksanht/+2//7WRiYmLba1AVJhcFoEw97r63e+fmcH8vto7sCwBQWaUINLSaSC16NLQz2VqnYjLo6MIc3aJjDojoVRHzQUT35pizoZGnT58mk5OTW5Y1Gsf52bNnm8/ZOnnLKK+oSI5z6e7du9vOC6otOzciz8h7DQCA0etxlw2T9NZ71/V0AAAoU6BhWLIAQyZ6VMTk0LG8WVfp6Ba9e/futrYxNfXraLucTE1NJQcOHGi4jPKKlukRZNi/f/+2ni5UW3Zu/OpXv9p2bqyvr+eO4wwAQG+iJ0IECVrNs6DHHQDAmAUaosdCM/3ozRBqgwy1wYbo5dBs8umoMGy/QjkbLmWi5j15yyiz6MnQ2XlB1c8N5woAQH/EcEfmXAAAKNZYjOMSQ4+0CkYUJRsyKSafBgAAoByi10LMr1D7iOGRAACoSI+GEIGECCg06tEQEzUX6dSpU+kE1DEMUp5G+wIAAMDoqR8aKZhnAQCgYj0aTp8+nVb8NxLzJhRpdXU1d3z0LMBw5MiRQrcHAAAAAABlVJoeDdHD4NKlS2nvhdphkpaXlzfnTihSzL+wsLCwbXnMzdBo/gYAIEl+9MEnyVePn+QeirzJNgEAAIByK02PhggkROX/xYsXtyyPYMCVK1e2zdEwMzOTPlpNIN1oCKQLFy4kc3Nz23o5xPbztgcAfCOCDDHJZt6jUQACAAAAKK/S9GgIUcG/uLiYnD9/Ppmenk7W1tbSYEAEIOpFICBv6KN4bwzBFEGD7P8xsXOsu7S0tOX9EcSI9OPvbNim69evJ7Ozs339nAAwDmKCzdd27Uz//vLx18mLjWHvEQAAAJBUPdAQ5ufn21qv0STOecMhNRIBhtrgAwDQvggyfPrusfRvk20CAADA+CrN0EkAAAAAAMDoKV2PBgAAAChSDPEXve+yvwEA6IxAAwAAAJUW8wh9sS7AAADQLYEGAAAAKunVXS919RoAAFsJNAAAAFBJH7/z9rB3AQBgLJgMGgAAAAAA6JpAAwAAAAAA0DWBBgAAAAAAoGvmaAAAACD1ow8+Sb56/GTL0fjy8deODgAATQk0AAADE5VVb713ffNvAEZLBBm+WJc/AwDQGYEGAGBgXmwkKrAASmByIkle27Vzy7JXd700tP0BAGC0CTQAAH3XrHJKxRXA6Ikgw6fvHhv2bgAAUBICDQBA3338ztuOMgAAAIypyWHvAAAAAAAAUF4CDQAAAAAAQNcMnQQAAFBBP/rgk+Srx0+2LPvy8ddD2x8AAMpLoAEAAKCCIsjwxbrAAgAAvRNoAGBkLC4uJvfv30+mp6eTtbW15MSJE8nJkyf7lkan23v06FFy9uzZ5MyZMw3XizRmZ2fTdeL53r17ydLSUvreeAaAUTM5kSSv7dq5Zdmru14a2v4AAFA+Ag0AjIS5ublkZmYmWVhY2FJp/+DBg+TcuXOFp9HJuqdOnUr27t2b/n316tU0iNBIvD8CGPHIHD9+PLl27VpbnwEABi2CDJ++e8yBBwCgawINAAzd6upqcunSpWRjY2PL8ggCHDp0qK1AQydpdLq9K1eupM/ROyHe18y+ffvSQMSNGzfSvyN4EYEGqiVv3PPaVsIfv/P2wPcJAAAA+kWgAYChiyGFYpihetmy6EXQagilTtIoYnuNRM+H+fn5rt5LucUEqm+9dz3925jnAAAAVMnksHcAAJaXl9PW/3n27NnT1rBDnaRRxPag3ouNbwIM9UGG7+3emT5iDHQAAAAYR3o0ADB0MSRRo+GFoofAzZs3C02jiO21EsMzRTqHDx/O7T2R5+jRo8mOHTu2LLt161by/Pnz5MWLF+lzOWRDUm2UaJ87U/udfPc7v1Pzmbf67ndeSv71P/m76d9/b+HfJl+sPxnr4zJM5fud9ObgwYPbllXlswMAAKNHoAGAoXr06FHT16OHQat1OkmjiO01E5NBnz9/Pp2bIeZ6iN4TMe9DzPPQqBdF5unTp8nk5NbOhnfu3EkrTyPdu3fvbnt9FD179mzzOfZ/HNV+J+//g1ebrpsdgyocl2Eq2++kV0+ePMk9BsBwh87LmI8HAKgagQYAKFAWYMhEz4mYHDqWr62tNX3vyspKsnv37txWylF5un///m09HkbR1NSvY6+Tqamp5MCBA8k46uY7qcJxGaay/U56dfv27W3L1tfX015ZwHCGzgMAqDKBBgCGKnoQNNNO74JO0ihie83UBhlqgw3Ry6HVJNNROdqogjRaaDd7fbRkkxFMlGR/u9P5d1KN4zJM5fqdFK+qnxuGJXot5PVuiMADAEDVCDQAMNJiKJRWQw4VmUYR26uXpReTTDcLNAAA5fHxO29vWxZDKOndAABU0fgPYAvAyIteBlHB36iHQUyoXGQaRWwvz6lTp9L5GBpptE0AAACAMutLoOGHP/xhP5IFYAT0I48/ffp0cu/evYavx/wGRaZRxPbyrK6u5o6PngUYjhw50lW6APTG/QkAAJQw0BBDQ/yrf/Wv+pE0AEPWjzw+egJEJX39/AjLy8ubcxwUmUYR28sTwyLF8akXczM0mr8BgP7r5/1JlCVRrmR5fScisB1z+ESZFCIIHv+fm5vrw54CAEAJh066e/du8tOf/jR5//33k/X19X5tBoAhKDqPj4r9qKS/ePHiluULCwvJlStXtk3gPDMzkz66TaPT7WWywESjIZAuXLiwrXIoKo9iO83SBaB8ZVcEFyLPj8BAN0GGrDxZXFxMh92bmJhIy7YoN5aWlnrePwAAKP1k0FGZ8pOf/CT9+7e//W1y+fLltHVODBnx4x//uB+bBGBA+pXHR7pR2RIVNtPT08na2lpagZM3eXJU2OcNUdRJGp2sG+vEZ8xanMb/o3Vs7ENtZVDsVwQrIp34Oxue6fr168ns7GzXxwaA0Su7Is0Q6Vy6dKmrNPbt25ecOXMmuXHjRvp39HDotlcdAACMXaAhu4gPL7/8cnL27Nn0788++yxtQXT//v30gvqNN97ox+YB6KN+5vHz8/NtrbeystJzGp2sG8GDdkWAQUtUgNEyqvcnEbDupNwCAIBKBRryfP7552lLn6h8iaEnontxtBqN7sExKefu3bsHtSsAFEweD0DZKLsAAGDEAw0xRnU27vUf/dEfpcGFGG5iY2MjnQgzhpR48803t3RdjmErDKsEMPrk8QCUzaiXXbEvN2/eTA4fPmyoPQAASqkvgYbswn15eTm9eI9xRj/66KMtXZbruy7HBX10W/7n//yf92OXACiIPB6AshnVsismg455f2Juhgh4xP7FxNAx/0PM2dCOo0ePJjt27Niy7NatW9vWe/78efLixYv0+Vsbm89bl9O98h3T/HODcXTw4MFty3zvAIx0oCGGRopJNX/2s5+lF8xxsd5KO+sAMHzyeADKZlTLrizAkIkASMwVEctjf9vx9OnTZHJycsuyO3fubFsvKpIjsHH37t3N9Z89e7b5nPceOlfGY5p3bjCenjx5kvv9A8DIBhpmZ2fTrr9546DGhGf18zH89Kc/Tf7kT/4knbMBgNEmjwegbEa17KoNMtQGG6KXQzanXSsrKyttzXcXrZajInn//v2bPSCmpn4dryRTU1PJgQMHuvwU1CrjMc07NxhPt2/f3rZsfX09zQcBYCQDDd/97ne3LYuux9EqJy6EP/vss/Ri/wc/+EH6WrTaycZOBWC0yeMBKJsylV3ZkEnXrl1rK9AQFcPtVg5Ha/Wt609sPqtgLko5j+n2c4Oq8J0DMNKBhrwWNdH1+NixY5v/j/FOswv5mHgtm3wNgNEmjwegbEax7Dp16lRy7969NNCRJ4ayAQCASgcaJiayVhyNRQsdEz8DlI88HoCyGcWyKyanzpvwOQswHDlyZGD7QvG+fPx18tZ717cse3XXS8nH77zd8D0/+uCT5KvH28fQb+e9AABjEWiIMUyjNU50Pw4x/ukPf/jDhuvH63njkQIweuTxAJRNGcquGBZpYWFh2/KYmyG4Xyq3FxtJ8sX61x0FFZqtDwBQiUDDz372s82L4rggjhZDGxsbuevu2bMnXf/s2bNFbBqAPpPHA1A2gyy7Hj161HKoo5mZmfQ55oTIxPwPc3NzydLS0pZeDhcvXkyuXLmS7hflEz0P8no3ROChVgQZmgUWvrd7Z8P3AgCM/dBJ0SonJlGLFkQfffRRkUkDMGTyeADKpp9l1/nz59NeExEcyP4fwy/t3bt3S/AgRNAgltcvix4NEWyIvyOtcP369XSfKae84Y1iCKVGQYXJiSR5bdc3QYW8YZKavRcAYKznaIhxRs+cOVN0sgCMAHk8AGXTr7Irb9ijRhpN+BwBhvqgBNUSQYZP3/12UnIAgLKa7EeiP/nJT1qu88tf/rIfmwagz+TxAJSNsgsAAErWo6FdMfbo5cuXh7V5APpIHg/tyZsMNG/oDKD/lF0AADDEQMP777+f3Lx5M/njP/7jzWVHjhxpOWFaNgYpAKNLHg/Fi4k9Y8ztYNxtKJ6yCwAAShho+Jf/8l8mn3/++ZZAw9raWnL48OF0PNQ8GxsbyYMHD3rdNAB9Jo+H4r3YyA8wfG/3zs1ARKwDdEfZBQAAJQw0rK6ubgsaRJDh5z//edP3PXz4sNdNA9Bn8ngoTgyF1Oy1bJik6O2gpwN0T9kFAAAlDDS8/PLL6aPW0tJSy/ctLCz0umkA+kweTzdzDESLfLYz3wIMhrILAADGZDLo119/vZB1ABg98nhqRZBB63tg1Cm7GLU5egTlAYBx05dAwy9/+cvNv9944430+c/+7M+SK1eupMMsnTlzJvnxj3/cj00D0GfyePJMTiTJa7u+mWOgnaGCAAZJ2cUoz9EDADAO+hJo+MM//MPks88+S+bm5tJAw/Xr15N/8A/+QXLy5Mnk3LlzybVr19L1BBsAykceT54IMnz67jEHBxhJyi5GeY4eAIBx0JdAw8zMTHoxnzl//nyyb9++5PLly+n/jx07lrz//vv92DQAfSaPB6BslF0Mkzl6AIAqmBzERlZXV9PeDACMH3k8AGWj7AIAgBL0aHjllVc2/45hkyYmJpITJ05sWWdjY6OrtBcXF5P79+8n09PTydraWppuN0GMR48eJWfPnk3ni2j2/qK2BzAu+pnHA0A/KLsAAKCEgYa7d+9u/r2wsJDs2bMn+cEPfrC57PPPP99ysd+umPMhuj1Hmpmo3IoJpmPuh3acOnUq2bt3b/r31atX00BDP7cHMG76lccDQL8ouwAAoIRDJ0UF/eHDh5MDBw4ky8vLyYcffrjZ8vUP/uAP0sr76FHQaffmS5cuJfPz81uWRyVXbK9dV65cSZaWltJ5IwaxPYBx0488HgD6SdkFAAAl7NHw+uuvJzdv3kx+8YtfpJNAv/zyy1t6FMQjhtroRAQHZmdnty3PlkXvhCKHNBr09gDKoh95PAD0k7ILAABKGGjIvPnmm1v+f+zYsc2///RP/7SjtKLVbF7Ff4hhO65du1Zoxf+gtwdQNkXm8QAwCMouAAAoYaChVY+BH//4x22vf+/eveT48eO5r8WcC9G6tki9bO/o0aPJjh07tiy7detWg7WzCVM3kufPnzdZRlnFd/jixQvfJQ3Pjb/9t//2th4AZf/td5rHA8CwKbsAAGAEAw3vv/9+cvny5dzXYuzuqMhvV6uxvqOHQZHjgfe6vadPnyaTk1unv7hz507uus+ePdt8ztbJW0Z5RUVyTCAekxDWnxdUW3ZuRJ6R99ooKzKPB4BBUHYBAEDJAg0//elP04mUY7LQGL+73v3799PKtXG1srKS7N69u611p6Z+HW2Xk6mpqXRi1UbLKK9omR5Bhv3792/r6UK1ZefGr371q23nxvr6etp7ahRVPY8HoHyUXQAAUMJAQ7RkbVXJdPr06bbTix4EzRTZm6GI7UWFYfsVytlwKRM178lbRplFT4bOzguqfm6M8rlSdB4PAP2m7AIAgBIGGk6cONFynYWFhcK2FxVeea1q+2VQ2/vy8dfJW+9d37Ls1V0vJR+/83bftw0wKnk8APRK2QUAAP3VlwHj2+lh8Nlnn3Xcy6BRC9rYXgzhUaRBby/Pi40k+WL96y2Prx4/6ft2AQadxwNAPym7AACghIGGc+fOpZOtff755w3XWVpa6ijNGIaj2eSi7bRSGuXt1fda+N7unVsek9loSgBD1o88HgD6SdkFAAAlHDopLuSj1dD58+fTngExoWntvAfxWrNK/DynTp1KJx+N99amtby8nD4fP368wE8w+O3VyhsaKYZQih4NAMPWjzwe2M4QilAcZRcAAJQw0HDt2rV0aKGf/OQnaQVUvfv377ecSLReVOyfPHkyuXjx4paxv+PvK1eubJvAeWZmJn1eW1tr2n260X50uj2AquhHHg80HkIR6J2yK0n+4b/4d8kXj/5jMjX16yRJJjYDmgAAMLKBhpgo+ec//3nLoYk6FRX8i4uLaSva6enpNIgwNzeXBgTqZa1s68V7o6Xt6urq5v/jxiPWrR/qo5PtAVRFv/J44NshFOtFZWAEHoDuKLuS5Dd/+ST5zV89T5IkHgAAUIJAw4cffthyndpeAp2Yn59va72VlZVCttvu9gCqop95PGAIRegHZde3Yu6313btbBngZLQYTg8AqGSg4c0332y5zuuvv96PTQPQZ/J4AMpG2fWt13a9lHz67rEhfht0w3B6AMCom+xXwr/85S+TH/7wh+mQQ3/0R3+0ufwP/uAPkj/7sz/r12YBGAB5fDX96INPkrfeu77lYXxvoCyUXZRR9Db53u6dWx7RKwUAoBKBhl/84hfJD37wg+Tll19Ofvazn2157Q//8A+Thw8fphf6AJSPPL66vnr8JJ2cuPZh3gCgDJRdlNXH77yd9kCpfdQPfQUAMLZDJ0VwIeZIyIZHqu3REH7yk58k77//fvLGG2/0Y/MA9JE8HuN7A2Wj7GIcmbcBABj7QEMEGMzBADCe5PFES0rjewNlouxiHJm3AQAY+0DDd7/73S3/39jY2LbO/fv3+7FpAPpMHg9A2Si7GLd5G/J6NxjOEAAYu0DD3bt3k3//7/998v3vfz/9/8TE1tmqYtgkAMqpn3n84uJiGoienp5O1tbWkhMnTiQnT57sWxqdbu/Ro0fJ2bNnkzNnzjRdr4jPAUBx3J8wbvM21Hvrvevp3EkAAGM3R8O+ffvSipUjR46klSx79+5N7t27lywtLSV79uxJbty40Y9NA9Bn/crj5+bmkpmZmWRhYWFzWWzjwYMHyblz5wpPo5N1T506lX7GcPXq1TTQ0M/PAUCx3J8AAEAJAw1RyXTz5s20smV+fj5dFpVPIf4fF/oAlFM/8vjV1dXk0qVL24bai8r6Q4cOtVVB30kanW7vypUr6XMEU+J9/fwcABTP/QkAAJQw0BCiteu1a9eS3/72t2mFVLQEffPNN/u1OQAGqOg8PgIVs7Oz25Zny6IXQauhhzpJo4jt9etzANAf7k8AAKB/JvuYdurll19Ojh07JsgAMIaKyuOXl5fTCqBGrVAjqFFkGkVsr9d9AGA43J8AAECJejT88pe/TMejjlauUemye/fufm0KgAErOo+PIYmOHz+e+1psI3pNFJlGEdvrdR/yHD16NNmxY8eWZbdu3UqeP3+evHjxIn0ermxIqI0R2JfhGvx34ti3Mjq/k8E4ePDgtmVV+ezdcn8CAAAlCTSsr68nZ8+eTYeGqBcTYcaQEr//+79f5CYBGJB+5fGPHj1q+nr0BGi1TidpFLG9XvehkadPnyaTk1s7G965cyetPI3Azt27d7e9PkjPnj3bfI79qrJBfyeO/eh9J8P25MmT3GPAVu5PAACgZIGG999/Pzl//nw6AWa0bo0KlWi9GTd8UbHy85//PF0e67z33ntFbRaAAZDHD8bKykpu75BopRyVp/v379/W42GQpqZ+HXuTTE1NJQcOHEiqbNDfiWM/et/JsN2+fTu3Uj2uv/mGsgsAAEoWaPjwww/T4MHPfvaz5Ny5c+m4p/ViwtDLly8nP/3pT5Pp6enkn/2zf1bEpgHos37n8RGYbqad3gWdpFHE9nrdh0aicrRRBWm00G72+mBMbD5XoSK3lcF+J459O0bjdzI8Vf3cw74/ifw9evydOXMmOXnyZMfvX1xcTO7fv5/uw9raWtpLsJt0AACg1IGGX/ziF8nCwkLaCvP1119vuF5c3MdF/unTp9OeDTGG9fe///1eNw9AH41CHh894xpNsNyPNIrY3iDTBWA4ZdepU6c2e5DEsIIRaOjU3NxcMjMzk+5vJgINUWbEvgEAQFn0PIBttBKKcbmbXcTXt/j86KOPDJ8EUAKDyuPjfVGp0qil6OHDhwtNo4jt9boPAJS77Lpy5Uq6nRgathurq6vJpUuXkvn5+S3LI+gQAQgAAKhUoCG69x47dqyj90RroYcPH/a6aQD6bFB5fLQmvXfvXsPXo3VnkWkUsb1e9wGA/ijL/UkEKWZnZ7ctz5ZFLwkAAKhMoKHbYSBajWUNwPANKo+P4SeiZWf9PAbLy8ubFUBFplHE9vr1OQCoxv1JlA2N9jX25dq1awPdHwAAGGqgYWIim5ywM9l4pgCMrkHl8VEBHxNfXrx4cdvwETE0RX3lT4xnHY9u0+h0e5ksgNBoeKRu0wWgOGW5P4kecI22Gctv3rw50P0BAIChTga9sbEx0BsAAAZnkHl8VMQvLi6mY11PT0+nQ1/EGNVRcV8vKuzzKmc6SaOTdWOdqBCK3grZ/6OlaexDDH3RbboAFK8M9yf1Pd/yyrlW62SOHj2a7NixY8uyW7dubVtvo+b5+fPnHewt5fDtN9zp9xvrv3jxwnlRAQcPHty2TH4AwMgEGrq9IG/UGhSA0THoPL5+QsxGVlZWek6jk3WjR0InOtkHAIpVtfuTp0+fJpOTWzuq37lzZ9t6z54923zOe51y6+X7jSBDnP93797ddi4xXp48eZL7/QPASAQaouXmZ5991nHLoaxVKACjSx4PQNmUoexqNZReu70ZsuD77t27W643NfXraLucTE1NJQcOHGg7fcqhl+83WrRHkGH//v3bescwXm7fvr1t2fr6uqGtARiNQEPoZvxQQycBlIM8vpp+9MEnyVePt7Z6+/Lx10PbH4AqlV3RurzdSa2jYridyuGJmmeVyePo22+4m+83ejK0ey4xXnznABSl536Rs7OzaVe7Th9vvvlmMZ8AgL6Rx1dXBBm+WP96y+NFd8OeAwxUWcqu6NXQaLim6NFw+PDhge4PAAAMtUdDuy1tinofAIMjj2dyIkle27Vzy4F4dddLDgwwsspSdp0+fbppz4sTJ04MdH+obo/FmEA65nb43p6vkn/zT//+kPYMAEiqHmj46KOPBvo+AAZHHk8EGT5995gDAZRGWcquU6dOJZcuXUp7L9TO2bC8vJw+Hz9+fKD7Q7V6LOaZmto+UTAAwMCGTgIAAKiqbOLmRsMghZmZmfRRKwIJJ0+eTC5evLhl+cLCQjqhdasJo6HXHovf270zfcTfAAAjMRk0AABAlZw/fz65d+9esrq6uvn/a9euJXv37k2Wlpa2rBtBg1heLwIKi4uL6Xunp6eTtbW1ZG5uLg1AQDe+fPx18tZ717cNefjxO2837LH41nvLyRfrejMAAL0RaAAAGLMKJKD/oudBu1ZWVhq+Nj8/X9AeQZK82EgaDo3UypePnyhjAICuCTQAAFSgAgmA8RVB57zgdJQb7VLGAAC9EGgAAKhABRIA4yuvZ1v0gGsnOP3d77yUPHv2LJmaiuqBbyZsUMYAAJ0SaAAAGOMKJABo5l//k7+b3LlzJzlw4ECyY8cOZQwA0JXJ7t4GAAAAAAAg0AAAAAAAAPTA0EkAAAAwpmK+hRhuL/sbAKAfBBoAAABgTL3YSMzpAwD0nUADAAAAjJlXd73U1WsAAN0QaAAAAIAx8/E7bw97FwCAChFoAAAYk7G3a1uqqmACAABgUAQaAABKzNjbAAAADJtAAwBACeWNrx29GyLwAAAAAIMk0AAAUEJ5QyPFEEpfrH89lP0BAACguiaHvQMAAAAAAEB5CTQAAAAAAABdM3QSAAAAkDv3TwzLVz9HUN7wfQBAtQk0AAAAANu82EjM/QMAtEWgAQAAANjSayGvd0MEHgAA8gg0AAAAAJvyhkaKIZS+WP/aUQIAcpkMGgAAAAAA6JoeDSVjMi4AAABG1Y8++CT56vGT3NdMJA0A40ugoWRMxgUAAMCoiiCDIZYAoHoEGkrCZFwAAACUxeREkry2a2f6t4mkAWD8CTSUhMm4AAAAKIsIMnz67rH0bxNJA8D4E2gAABgz5nQCAABgkAQaAADGjDmdAAAAGKTSBRoWFxeT+/fvJ9PT08na2lpy4sSJ5OTJk31JI5bPzs4mZ86cSZ/v3buXLC0tJY8ePUqfAWAc/OiDT9KJG+tbxFM+5nQCAABgGEoVaJibm0tmZmaShYWFLcGABw8eJOfOnSs8jVgWQYl4ZI4fP55cu3atkM8DAKMgggxfrAssjANzOgEwysPzGdoPAMZXaQINq6uryaVLl5KNjY0tyyNgcOjQobYCDZ2msW/fvrQ3w40bN9K/IyARgQYAGEeTE99M3NiqhTwAUF29DM9naD8AGF+lCTTEUEUxfFG9bNnVq1dbDqHUaRp79+5N5ufnC9h7ABh9EWT49N1jw94NAGDMhucztB8AjL/SBBqWl///7d1PbFzHeQDwIcXGThvJhBgLRhAUrWhfYsQHSjTsJj2ZvhTwjZR7aHsoYOpQJC16ECEfeilgg7z6EtK3AgEa0Jc0PRQwhVwMxIUpHnwIYFjUxUHgWKDIUEEiqZS2+FZdlkvucv+9/fPe+/0AgtTbt7OrndmZ3ffNzLfRMEgQJicnq9sZtQo0ZFEGAAAAlE0v2/PZ2g8Aii83gYZIxNxs26JYebC5udm3MmLLpbjt8uXLTQMVR7388svpzJkzdcc+/fTTJmfXpn9U0qNHj1qWnd19GZSom8ePH6sjmraN7373u2lsbOzEbQAAAACQB7kINOzt7Z16e6xGaHVON2VEMuilpaVqbobI3xArIiKXw/r6ejVnQzMPHz5M4+Pjdcc+//zzhuceHBwc/m52TjO93JfBiQvJ0ZZu3bp1ol1QbrW2EX1Go9sAAAAAIA9yEWgYllqAoSZWQ0Ry6Di+vb3d9H43b95M586da+sxJiZ+FXOX08TERHrhhRc6en693JfBiZnpEWR4/vnnT6x0odxqbeOXv/zlibaxv79fXWkFAAAAAKMuF4GGWG3Qy2qFbss4GmQ4GmyIVQ6nJZ+OC4btX1CubZcy1sVF6F7uyyDFSobO2gVlbxvaCgAAAAB5UYh9XGLrkVaBhKzKqG2ZFImjAQAAAACg7HITaIggQAQDmq1GiETNWZaxsLBQzcfQTLNyAAAAAACgTHKxdVK4cuVK2tzcbHp75E3Isoytra2GCZ9rAYbZ2dk2njUAnVhZWUk7Oztpamqqmgsn+uVm29RlUUY/zo3jMzMz1Zw+8fv27dtpdXW1GtCO3wAAAABFk5tAQ6wwWFtbq16oObrF0cbGxmHuhCzLiItHy8vLJ8qI3AzN8jcA0L2rV6+m6enpur43LtpHgLfdPreTMvp1bhyLoET81MT4Yss9AAAAoKhyE2iIizRx8f/dd9+tu9ATf6+vr5/IrxAXhELMOu2mjOvXr1cvLB2dfRqrHOK+jR4PgO5F/xqB4EqlUnc8+ufYxq6dQEMnZfTr3BCr4WI1wyeffFL9OwIS7QTDAQAAAPIqN4GGEBf4Y4bo0tLS4dYVEQxotHVFBALOnz/fdRlx/7iIFLfF37H1Rbhx40Z1KwwAshNB3UZ9a+1YrCZrtYVSJ2X069wQY8+1a9da/I8BAIrjq3v30yvv3Dj8GwAon1wFGkK7F29u3rzZcxkRYLCfNkD/xRZ2zYK40RfHtkOtAg2dlNGvcyEvF4Fqnj37VPrZD74/tOcERdFrjiH5fci7x5WUvtwXYACAMstdoAGA4olVY822F4oVApubm5mW0a9zj4otl+K2y5cvt70S7uWXX05nzpypO/bpp5+mR48epcePH1d/90dtW6hKHx+jWPpfJ9nXb+OLQMWp83zVSe9eeumlE8fK8n8vYo4h+X3IqwhYd3MbAFA8Ag0ADNXe3t6pt8eqgVbndFJGv849erEotueLi0xxgSlWREQuh9i6L3I2nObhw4dpfHy87tjnn39evXga5d66devE7Vk4ODg4/B2PR2v9rpMsnf2jlA7+uD6AdfcPj6qBhyLVeZ7qJAsPHjxo+BqQvxxDQX4f8sqqOACgRqABADJUCzDUxGqISA4dx2M7jVbb/p07d67hLOW4ePr888+fWPGQhYmJX8WjpImJifTCCy9kXn4R9btOsvRf/3yyTr+3/PP05f6DQtV5nuokC5999tmJY/v7+w1zlDHaOYaC/D4AAOSdQAMAQxWrAk7TalVBp2X069yaRrNXI9gQqxxaXXCKi6PNLpDGDO3Tbu/N2OHvMlygzUp/66Tfilnn+a6T3pX1/z1M8vgAAMATAg0AjLTYCqXVlkNZltGPc2vnSBwNUCxZ5BjqV36f5tl45PNgVHLcPGmVX927n155Z6Pulm9+46n003/4iwE/n+KT3weAfhJoAGDoYuVAXLRvtmogLrhkWUa/zl1YWKhedIotkBppVs4gvfHeR+nOvfq93eMLPgCDzzHUz/w+x8nHw6jluKm1ycgZFNv5Hb+tKDmERon8PgD0k0ADAEN35cqVU2d9xoWXLMvo17kxE7XRBaFagGF2djYNWwQZvtwXWAAoen6f4+TjYdRy3Dw3eSdNTByf/PCgGngoUg6hUSK/DwD9JNBQAE+Wmt6oO/bs2afSz37w/aE9J4BOxEqAtbW16szPo3kRYkZnaLYtRbdl9OvcyL+wvLx84rlFboZm+RuGZXwspQtnnz4xdgAwuBxD/c7v0zgzi3wejEaOm//84V+eOBbfa59MiChWDqFR5nUGICuDWxdJ3zxZanq/7uf4thgAoywupsRFlHfffbfueFy0j20jjl/MmZ6erv50W0a/zr1+/Xq6evVq3XmxyiHu2+j/MUwRZPj47dfqfgSoAbITq9l66feP5vcBAIBRZ0VDjjWaeRqrGyLwAJA3cSF+ZWWlOntzamqqulVEXLRvNIszLtxEks1eyujHufG8IgARt8Xfka8h3Lhxo+2kngCUK8dQHvL7AABAKwINOdZo5un/LzUFyJ9r1661dV6zizGdlNGvc+Oi0+rqatvlAlDuHEN5yO8DAACt2DoJAACgC7EaIQIFx/MxdJJjKFbHNdoeaRTz+wAAQDMCDQAAAEPKMZSn/D4AANCMrZMAAACGlGNIfh8AAIpAoAEAAGCIOYbk9wEAIO9snQQAAAAAAHTNigYAgJL66t799Mo7N+qOPXv2qfSzH3x/aM8JAACA/BFoAAAoqceVlL7cv193LP4t+AAAAEAnBBoAAEomVi0cdzTgIPgAAABAJwQaAABKptHWSG+891G6c+9B28EHAAAAqBFoAACg7eBD5HWILZcAAACgRqChoCR3BAD6EXyI/A1WNwAAAHCUQEOJkjsCAAAAAEDWBBpKkNzRFgcAAAAAAPSLQEPB2OIAAAAAAIBBGh/oowEAAAAAAIViRQMAAAAwMmL731feuXFim+BGK/gBgNEg0AAAAACMjMeVlL7cv9/09jfe+yjdufeg4W0CEgAwHAINAAAAwNBFkKDR6oYIPBwVQYbTAhEAwOAJNAAAAABD12hrpNhCqVlQYXwspQtnn24akAAABkcyaAAAACB3Isjw8duvVX9qAQcAYDisaAAAoGMSdQIAAFAj0FAiLggAAINK1AkAw+B7LwAMh0BDibggAAAMKlEnAPQrgBB/N+N7LwAMh0BDCbggAAAMK1EnAGShVQChl++9b7z3Ubpz70HTchuNfQBAPYGGEnBBAAAAgDxqFEBodFsv33sjyCBgDgC9EWgAAAAARtIgVxOMj6V04ezT1b9tCwgAnRnv8HwAAACAwokgw8dvv1b9qQUcAID2WNFQckcTatXYgxKgGBrtN3xa8kQAAACAbgg0lFyrhFoA5Jf9hhkmiTUBAADKQ6ChpBol1LIHJUAxHd1vuJ3EipDFSkkTGQAAAMpDoKGkGiXUigsDLgoAFG+bpNp+wzCslZLPnZNYEwAAoMgEGgCgAGyTxDCdtkLmaO4nkxoAAACKSaABAArENkmMykpJABjFrf2a3d7JfY8G0QGAJwQaAKBAbJMEANDe1n5Z3LfR9pU1AhIAlIlAAwAAAFCqrf1and/ovrG6IQIPR1c5yHMIAE8INNDxstIaszMAAAAYVb1sb9TovrVcQ81WOTx37um6gAQAlIlAA5kuKwUAAICyrZA4OhGvFpAAgDIRaKDjZaVmZwAAvZBYE4A8kgAaAJoTaKDjD01mZwAAvbB6EgAAoFgEGuia2YgAQCdOS6x51BvvfZTu3HvQtAwzSgEAAEaLQAOZzkaMfx9PJO2CAADQKrHm0QkM9rUGAADIF4EGMpmNePSCgIsDAP3VaLZ3XKSFIm6n9Ny5p6u/5YgCAAAYXQINZDIbsdlFr+NbIQDQu+hvBXUp8gSGo7fVPnfIEQVAntkWEICiE2hgoFsh1NhOCaB342MpXTj79In+FfKkl3wL9RdtKung4CBNTPwqpTTmswYAQ2dbQADKRKCBkcnl0IiABEBzEWT4+O3XvESUVuPVPY+G9GwAoJ5tAQEok9wFGlZWVtLOzk6amppK29vb6fXXX0/z8/N9KyOLxyurLHI5FCEg8dJLL6UHDx6kzz77bNhPhRGjbeSrjx/m2KGtjB51Mnqreyq/30uVSkrjfzJp60YKP3514ze/+U1KX5/MrDyKw5iWPdsCAlBGuQo0XL16NU1PT6fl5eXDY/EB/O7du2lxcTHzMrJ4vDJrN5dDI70kl7b3JeTTKPfxxg4YzW0oaknQY3XPvf94pxrY/+bf/yh9ud/6swbkdfwCRt+oToIDgH4aq1Ri7tfo29raSpcuXUrHn26z472W0c3j7e/vp2eeeSb99re/TefOnWvr/1XLY/DcOdtfdBOQqCWcPrpPeavARLzW3Ti6aqKTYMaLL754uKLhzJkzXT02xXRa2+imP8mzUe7j+z12tFPfx9uKsWP49O2Dd1oy6Bjb7/34h3WBhjJ+tirb2FHG8auX+v7zf/xxqnx9Mj137qn08dtzLc+nPIxpwx/bmuXdGkTAwtgBQOlWNKyurqaZmZkTx2vHPvjgg5ZLizspI4vHo3vtfqCqfTBrtfdlLyskjt6vNoPytDKOb/f0+7m3q1s5fG/559XklHndAgrK2scbOyAf21DcO2Xlw9HzjLPkefwCikleQwCKIDeBho2NjYYfwMPk5GT68MMPW34A76SMLB6P0dn7st0VEo20ClK0DGb831647Wzj0G5Oin6uwuj0tRqFiza2y8q/Ue7jjR0wGlqNNS/+KLuLNlmNgVmMvaMwzjI64xdQLP3Ka2jsAGAYchNouH37dpqba7zE9/z582lzczPTMrp5vNqy5ljifHwblF/84hcNyzq4/7v0+MHDdHD/IO3u7rb8P1Dv3/7mxVNfktpr2uq80/z1+/+ddn7X+Mv/1DeeSv/+1syp59356k6qpEq6cOFC08f46t7Dw79/fef3KSu/vpPS5X/56YnH6Pa8Vvcdlnb/b8P0t6/8afq7V/+s4W3f+c530thY/WqXR48eVX/nZHe7Qvfx/R47Ohk/ok+L240do6NWJ4yOZyb+Jx187WDg42wjWYy9B/e/1vAz4quvvnriWNnGjjKOX71893j84A+pMv413ztoypg2eI2+p572/bPdsc3YAcAw5CLQsLe3d+rtMdOn1TmdlNHt492792TR/q1btxp+STjNF3HOv556CiPq/FJ750UdD0O7j9vL8xvW/y1Pz+9mSumfurhf9Cux/3KRjXIfP4ixo5Px49lnn627zdgxfMfrhOHbXvqrgT5ev8fPbt7nZRg7yjh+ZfXdw9hBM8a04jB2ADAMuQg05MW3vvWt9MUXX6SzZ8+emJ0M0ImYpRgXEKJfofiMH0AWjB3lYuwAsmDsAKBUgYaYydPLTKBOy+j28cbHx9O3v/3tls8FoB1lmY06yn38IMaOYPwAslKWsaOM49dxxg4gK8YOALIwngrg7t27LT+kZ1lGFo8HQP77eGMHAFmMEf0sAwAABiE3gYb4gB0ftJvN9Ll8+XKmZWTxeADkv483dgAwKuMXAACMqtwEGq5cuZJu377d9PbXX3890zKyeDwA8t/HGzsAyGLs6WcZAAAwbGOVyPyTAxsbG9UP2bu7u3XLh5sd77WMLB6P4ou2MDMzk958883q7/iSuLq6Wp19Fr+PW1lZSTs7O2lqaiptb29X7z8/Pz+U5062os7feuutals4rU47aQNlai+j3McPe+woUzsoWj+u7vLRf6on8jR+9asNex8UU6/12ukYCQCUXCVH5ufnK9euXas7Njc3V1lfXz9x7sWLF6s/vZTRybmU08zMTATq6n6ijTSyuLhYWV5erjsW566urg7o2dIP0U9E3cZP1P9p/UMnbaCM7WWU+/hhjR1lbAdF6cfVXT76T/VEHsevrNuw90ExZVGvnYyRAAC5WdHQ6ayMS5cupfPnz6cPP/yw6zI6PZfyWVhYSLOzs+mTTz5JFy9erLaPubm5E+dtbW1V2+Txt1uz4+RPzPCanp5O6+vrDfuITtpAmdvLKPfxgx47ytwO8t6Pq7t89J/qiTyPX1m1Ye+DYsqqXtsdIwEAqsRaoLeZQu2eFzOCGmk1i5N82N7ePrUuO2kD2gvaQb77ce/hfPSf6om8y6INex8UU1b12u4YCQAQcpMMGvIs9tiNWUCNxJ67jWa1Ud42oL2gHYwe7+Hivfb6WvIuizbsfVBM6hUAGAaBBshALENeW1ur/m62LUQsk28kjm9ubqqHguukDWgvaAf57se9h7PVr9dePZF3WbRh74NiyrpeW42RAABBoAF6cPfu3bS0tFT9vbi4WP0d+57Gh/uavb29U8uIGWetziHfOmkD2gvaQb77ce/hbPXrtVdP5F0Wbdj7oJiyrNd2xkgAgJqJw7+AjkVCtPjQXRPJ0d58883q8UjiB8Bo048DgDESAOidFQ3Qg6NBhqPBhpjl88EHHxzOGjqN1QzF10kb0F7QDvLdj3sPZ6tfr716Iu+yaMPeB8WUZb22M0YCANQINEDGakn52k3wHEuQW30hoNg6aQPaC9pBvvtx7+Fs9eu1V0/kXRZt2PugmHqt107HSACgPAQaoEsLCwvVPUpP+xBfEx/mj/77+Kyiy5cvq4eC66QNaC9oB/nux72Hs9Wv1149kXdZtGHvg2LKol47GSMBAIJAA3Rpa2srnT9/vumH7tnZ2cNjV65cOTVpWuwRTrF10ga0F7SDfPfj3sPZ6tdrr57IuyzasPdBMWVRr52MkQAAQaABujQ/P99wyXBtv9Kje5rGjKD4sH58T9SNjY3DvU4ptk7agPaCdpDvftx7OFv9eu3VE3mXRRv2PiimLOq1kzESAKCqAnRld3e3sri4WHfs5s2blcnJycr6+vqJ8+fn5yvXrl2rOzY3N9fwXPIn6j661NXV1abndNIGtBe0g3z3497D+eg/1RN510kbvnjxYvWnlzIoT9vodIwEABiLl0DMBboTs4SWlpaq+6DWlidfv349zczMNDx/ZWUl7ezspKmpqbS9vV1dthyzhcivqP+o+5g1Fr+jLcQssVhqvrq62lMb0F7QDvLdj3sP56P/VE/kXbttOPbbj/dXo1nq3gfF1Gvb6HSMBADKTaABAAAAAADomhwNAAAAAABA1wQaAAAAAACArgk0AAAAAAAAXRNoAAAAAAAAuibQAAAAAAAAdE2gAQAAAAAA6JpAAwAAAAAA0DWBBgAAAAAAoGsCDQAAAAAAQNcEGgAAAAAAgK4JNAAAAAAAAF0TaAAAAAAAALom0AAAAAAAAHRNoAEAAAAAAOiaQAMAAAAAANA1gQboo5WVlXTp0qU0NjZW/VlaWjq8bXp6unosfl+9elU9AGDsAMB3DwAgl8YqlUpl2E8Cii4CCWtra2l3dzdNTk4eBiG2t7fT6urqsJ8eACPI2AGA8QMAyAsrGmAAIpgwMzOTFhYWqv/e29tLOzs7ggwAGDsA8N0DAMg9KxpgQG7fvl3dJimCDrGSYXl52WsPgLEDAN89AIDcs6IBBuTixYvVIENshTE7O+t1B8DYAYDvHgBAIVjRAAO0tbV1uH1SrGoAAGMHAL57AAB5Z0UDDNBPfvKTdPPmzXT37t3qygYAMHYA4LsHAJB3Ag0wIEtLS+n69etpcnIy3bhxI62traWNjQ2vPwDGDgB89wAAck2gAQYUZAgRZAgzMzNpcXGxuo1SJIkGAGMHAL57AAB5JUcD9DnAECsX9vb2qkGG3d3d6vEILly6dOnw+NzcXHr//fcPAxEAlJexAwDjBwCQNwINAAAAAABA12ydBAAAAAAAdE2gAQAAAAAA6JpAAwAAAAAA0DWBBgAAAAAAoGsCDQAAAAAAQNcEGgAAAAAAgK4JNAAAAAAAAKlb/wtGTY8DqsjgvgAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "samples = [uniform, logunif, gauss, lognorm, powerlaw, maxwell, beta_samp]\n", "titles = ['Uniform(0, 10)', 'LogUniform(1, 1000)', 'Gaussian(30, 5)',\n", " 'LogNormal(1.0, 0.5)', r'Power-law ($\\beta=-2.35$)',\n", " r'Maxwellian ($\\sigma=265$)', 'Beta(1.4, 3.6)']\n", "xlabels = ['x', 'x', 'x', 'x', 'x', 'x', 'x']\n", "\n", "fig, axes = plt.subplots(2, 4, figsize=(16, 8))\n", "\n", "for idx, ax in enumerate(axes.flat[:7]):\n", " ax.hist(samples[idx], bins=50, density=True, histtype='step', lw=2)\n", " ax.set_title(titles[idx])\n", " ax.set_xlabel(xlabels[idx])\n", " ax.set_ylabel('Density')\n", "\n", "axes[1, 3].set_visible(False)\n", "\n", "plt.tight_layout()\n", "plt.show()" ] }, { "cell_type": "markdown", "id": "ecee49fa", "metadata": {}, "source": [ "## JSD between distributions" ] }, { "cell_type": "code", "execution_count": 4, "id": "7f063aae", "metadata": { "execution": { "iopub.execute_input": "2026-06-13T10:26:26.874877Z", "iopub.status.busy": "2026-06-13T10:26:26.874780Z", "iopub.status.idle": "2026-06-13T10:26:26.876632Z", "shell.execute_reply": "2026-06-13T10:26:26.876242Z" } }, "outputs": [], "source": [ "from gwGenealogy.utils import compute_jensen_shannon_divergence" ] }, { "cell_type": "code", "execution_count": 5, "id": "bb68c47f", "metadata": { "execution": { "iopub.execute_input": "2026-06-13T10:26:26.877628Z", "iopub.status.busy": "2026-06-13T10:26:26.877563Z", "iopub.status.idle": "2026-06-13T10:26:26.880331Z", "shell.execute_reply": "2026-06-13T10:26:26.879953Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "JSD(Power-law beta=-2.35 vs beta=-1.5) = 0.0481\n" ] } ], "source": [ "# Compare two power-law samples with different exponents\n", "pl_a = sample_powerlaw_1d(5000, beta=-2.35, xmin=5, xmax=100, seed=10)\n", "pl_b = sample_powerlaw_1d(5000, beta=-1.5, xmin=5, xmax=100, seed=11)\n", "\n", "jsd_pl = compute_jensen_shannon_divergence(pl_a, pl_b)\n", "print(f\"JSD(Power-law beta=-2.35 vs beta=-1.5) = {jsd_pl:.4f}\")" ] }, { "cell_type": "code", "execution_count": 6, "id": "9fa053dd", "metadata": { "execution": { "iopub.execute_input": "2026-06-13T10:26:26.881293Z", "iopub.status.busy": "2026-06-13T10:26:26.881226Z", "iopub.status.idle": "2026-06-13T10:26:26.883623Z", "shell.execute_reply": "2026-06-13T10:26:26.883304Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "JSD(Beta(1.4,3.6) vs Beta(2.0,5.0)) = 0.0136\n" ] } ], "source": [ "# Compare two beta distributions with different parameters\n", "beta_a = sample_beta_1d(5000, a=1.4, b=3.6, seed=20)\n", "beta_b = sample_beta_1d(5000, a=2.0, b=5.0, seed=21)\n", "\n", "jsd_beta = compute_jensen_shannon_divergence(beta_a, beta_b)\n", "print(f\"JSD(Beta(1.4,3.6) vs Beta(2.0,5.0)) = {jsd_beta:.4f}\")" ] }, { "cell_type": "code", "execution_count": null, "id": "398dcbf3-9079-4518-863c-fff0a8956724", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "gwgen_env", "language": "python", "name": "gwgen_env" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.11.15" } }, "nbformat": 4, "nbformat_minor": 5 }