Skip to content

Instantly share code, notes, and snippets.

@HajimeKawahara
Created March 30, 2023 22:45
Show Gist options
  • Save HajimeKawahara/3638737db26f47301b9926ae666306a2 to your computer and use it in GitHub Desktop.
Save HajimeKawahara/3638737db26f47301b9926ae666306a2 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"#sampling"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"from jax.config import config\n",
"config.update('jax_enable_x64', True)"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"def RBF(t,tau):\n",
" Dt = t - np.array([t]).T\n",
" K=np.exp(-(Dt)**2/2/(tau**2))\n",
" return K"
]
},
{
"cell_type": "code",
"execution_count": 15,
"metadata": {},
"outputs": [],
"source": [
"from jax import random\n",
"from scipy.stats import multivariate_normal as smn\n",
"from jax.random import multivariate_normal as jmn\n",
"import numpyro.distributions as dist\n",
"import numpyro\n",
"\n",
"key = random.PRNGKey(4)\n",
"N = 1200\n",
"t = np.linspace(-5,5,N)\n",
"ave = np.zeros(N)\n",
"\n",
"def samplegp(tau):\n",
" tip=0.00001\n",
" cov = RBF(t,tau)+tip*np.eye(N)\n",
" #scipy\n",
" #ds = smn(mean=ave ,cov=cov , allow_singular =True).rvs(1).T\n",
" #jax\n",
" dj=jmn(key,ave,cov)\n",
" #numpyro\n",
" mn=dist.MultivariateNormal(loc=ave, covariance_matrix=cov)\n",
" dn = numpyro.sample('a',mn,rng_key=random.PRNGKey(20))\n",
" return dn"
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"0.0\n",
"1.0000049999875\n"
]
}
],
"source": [
"from scipy.linalg import det\n",
"print(det(cov))\n",
"print(np.sqrt(np.trace(cov)/N))\n"
]
},
{
"cell_type": "code",
"execution_count": 25,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAD4CAYAAADxeG0DAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAClaklEQVR4nOydd5wmRZ3/3x2e/DyT0+bZzLLkJEkUBMUERuQ8050Bs955nJjDqdzPU+/MZ8asCAbMBwoCSs5hI5tmZienJ4d+un5/dHi6n6efNDuzu7POh9cyM93VVdXVVd/61jdKQgiWsIQlLGEJixfyke7AEpawhCUs4dCwRMiXsIQlLGGRY4mQL2EJS1jCIscSIV/CEpawhEWOJUK+hCUsYQmLHOqRaLSrq0v09/cfiaaXsIQlLGHR4oEHHpgQQnSXXz8ihLy/v5/777//SDS9hCUsYQmLFpIk7fe6viRaWcISlrCERY4lQr6EJSxhCYscS4R8CUtYwhIWOZYI+RKWsIQlLHIsEfIlLGEJS1jkWCLkS1jCEpawyLFEyJewhCUsYZFjiZAfQcymCySyhSPdjSUsYQmLHIdMyCVJCkqSdK8kSY9IkvSEJEkfm4+O/T3gvn1T3LNn6kh3YwlLWMIix3x4duaAi4QQSUmSfMCdkiT9Xghx9zzUvYQlLGEJS6iDQybkwkgxlDT/9Jn/ltIOLWEJS1jCYcK8yMglSVIkSXoYGANuFkLc41HmTZIk3S9J0v3j4+Pz0ewSjmEMz2b48/ZRdH2JJ1jCEuphXgi5EKIohDgFWAmcJUnSCR5lvi6EOEMIcUZ3d0XwriUswYXdY0l0HfJF/Uh3ZQlLOOoxr1YrQogZ4Fbg0vmsdwl/f1BlY2oWlgj5EpZQF/NhtdItSVKb+XsIuATYfqj1LuHvGz5FAqC4JFpZwhLqYj448mXArZIkPQrchyEj/8081LuEQ4SuC3aOJshri4+rVWSDkBeKS4R8sWAmneeWJ0fJFopHuit/d5gPq5VHgVPnoS9LmGeMJ3McmEyT13ROWNF6pLvTFGTJIOS6WCLkiwWD0xkAplJ5lreFjnBv/r6w5Nn5d4DFSAstjnxJtHJ4MZPOox2iXsLcg5dwGLFEyP8OIBahWb9FDJYI+eFDoahz/75pHhuaPdJdWUKTWCLkxzAWM2NUIVrJpyC9FM5gIWFtmomsdoR7Mv+YSOaOadn9EiE/lmFS8sUsWrEZ8h//A3x6LRSyR65TxzisU9AinC61oeuIH17BxM/fe6R7smBYIuR/B1jMC7Oo61DUYO9fjAuD9x3ZDh3DkBb1Ga4GRh6le/g2Vm77BiTHjnRvFgRLhPwYxrGwMIUApvaULkzsPGJ9+XuBWIxHuFqY2FX6/eBDR64fC4glQv53gMW4MK0uC4D4YOlGYuRIdOfvAotRKd4IxPTe0h9DDx65jiwg5iOM7RLmgMNBXBezGZhFVITATbwTw0emQ38HcG2exxKm95EN9aD524gePDYJ+RJHfoRwOJhki44v+oUZP2j87Nq8xJHPIw5MpplI5ipvHOKEOepEejMHyIZXkmjdDOPHZvSQJUK+hKMSJe5QGMQ72Aod65Y48nnEztEEDx+YOdLdWHikJ8kHO8hEV8PsIGgem9cixxIhP0I4nFzyIhSRlwi5AJKjEO2DSBekJ49ov/4ecMzJytNTFPxtpGP9IHSY3n+kezTvWCLkRwiLQQGZymnc8uQo44kjx8EIAeTiBkcebIXsktfhQmERTMnmIQRkpij420lH1xjXnFZQxwiWCPkRwuFdM3NrbTZTAGA0fvidcFxcYS4JgSgE26CQBi1/WPowncovig13vuBSMB8ryKeQinkKgTZDtAIw9dSR7dMCYImQ/x1grgvTeuxIWL+4ZOS5BPijEGozLh4GrnxgKs0D+6cZ91IGHuOY83w5GjeAjBHWoeBvo+BvN051Sxz5EhYT6q4rXW9o9R1JKwQhgHwSAi3GIgTIzix4u5Opw8P1H6sQQizoaWYimeO2HWP1g6qZ8XkKgTaDI+lYB5NLHPkS5orBB2B8/r0SLfFHLXhO9eQYfPVc+N5lhgu813PmQjzi9ui5JGkpyEPj5pscBo7cfvejzZSuCSSyhaZC0s4n3b1txzh37p6YvwrLsHssiVYUpPN1AnyZm37BZzIBnRvdnp7HCJYcgg4HsrPwzYuM3z88DfL87J8DU2l2jCQ4dXUbndFAxf2aC/Pmj8D4NuPfzj/AlhfYt7SijiJL9vMLTcgLeoHvPP4dEvkEQ8khJjOTPDj2IMe3ns2rNr+RE/IJpgoBNF+L8UBmZmE7BEhHfPc6dNyzZ4q2sI8z+jsaKj+f/HNRFwsagli3mYw63ymXNPrjixp/9xwHj10P2TgEWxasf4cbS4T8cODAPaXfJ3Yak8mBVE4jEmj+UyRzBjeSqRKes6oZmZaDJ38Jp70GdvwBnvi5i5DftmOc3pYgbWFf031qFj/f9XM+ftfHKQrjHbpCXeSLhljjydm7ef+9d/O7nk7+WYWC31x4h0G0stiVnFb/Z9L1T2zzjcOyB5qfR67XVi4BgGYR8u4txs/x7bDqrIXp2xHAEiE/HBjfVvp97MkKQj6bKcyJkFtomuYM3GNYf2x+nsHdDj9SUWQ0nqU1ZBDyhRAv7JndwzW3X8O2qW30Rfp45XGv5BmrnsHalrVIksQ9+w4yEk9x18i3uWXoJu6Z/R1nyglOlyBwGAi5hcVqUz2XfWgxbV4Ws193buYNjlxTI8bf1tob23ZMEfIlGfnhwMwBUIMgybaLsHPNqHXZinmGFQp29dnQc7yhxS9kqhafbw7rYPIgr/j1K9g2tY0XrX8Jv3nxb/inE/6Jda3r7KNyUAkTUVt47fKX8/OhEc4JbeRv07fxtt4eRpIH57dDHljsopXFQ5LnBmuDrfuZTI68aBHytn7whQ2G6hjCEiE/HJjeD12boG21p6JFVRboM1RbzcOPQHs/hNoNDkXonuFh7cUyj12ayk5x1c1XkS1mec26D/D2k95HQCnJ96dTeWbSJYsRpZBitabx1u6XcPmqN3NPKMiVgzcxEB8gldO4d+8UhUPMMVkTi5QizkVRfaivejhPLw2L33MJhKSgW3NMlqHvpGMuCuISIZ8HCCEYjWerH01n9kP7GmhZWQoA5cBcGfJ6i7TqXB953JjMAB3rjZ9TeyuKzbeyczY3y2t//1qGU8Ncc8rnOLXzwgoRwAP7p7l/37Tdd7lQOho/s+9lfG8qT14UedXvX8Ud+54gnikwtWQqWIHFtP9MJHNNK0attVZXGpRPQiDmnsQrzzCYmcPkWHY4sETI5wFDMxkeG5xlaKaKeCI5ZsQKaVkOiUpCflhFk0XN2Fg6Nxh/t5tuy9P7KoqWunXolLygF3jv7e9lMDnIFy/6Iid31ZZP2hxl3lRWKcbR+CQ1xg/9myiKIh974G0MpBbWlOxoJ4gHJtPktEpl91zCKhwJEXkyp/HwgRm2j8Sbeq7hruaSCL+h6LRp+cozoZiDkceaavNoxhIhnwfkNONon9c8jvi6bpgfhtqgZZkRyU8I1zF0odaP58KMD4KuQcda4+9gK4Q6vAn5PK7srz78Vf568K988Gkf5Jzl55TaqPL2JY48BUDBtDoQgRbWFvJ87eKvUdQ1PvfkW3hg7B7POo51JLIFdo4meOKgmwjquuBJ81pTp6l5+tzNNFkw10wm31xi5IbnZi5ueAU7sfps4+eBu5pq82jGIVutSJK0Cvge0IsxFb4uhPj8odbbDIQQDM9mWdYaPPqUVLk4IIw4IbIKxbwRwS/QfmT6Y4lQ2teWrrWvMbh0vBfIoQ7pg6MP8o3HvsHz1z2fl256qetevfWomKIVS1kl/C2QjbO1aysfOO3zfPC+N/CZBz/OOStvoD04/2N6NBtyWF0rlDEQ+tHc6TJUCwOR13Qy+SKtVUxgXSEcaiFf4shtxPqM+X/gLjj37bX7VzaWuWKOoBokq2VRZIU9M3uYyEwQ9UeZSE8Qz8c5mDpIUS+S0TIMJYfY1L6JjFY6rb9000tZ17qudr+bxHyYH2rAe4QQD0qSFAMekCTpZiHEYVMLH5zNsu1gHK0oWN0ZPlzNNgbLVC7UVuIMEsMuQm5NFl0X7J9Ks6YjjDwPliyek9xKe9XeX7oW7TPiNFc8b+BQeqLpGtfeey094R4+dPaH7OuNZmyX85ZDh0nIAy2QMPq6tmUzr9/4H3zvqf/g7X96O9+59Dv4Ff8h9HZxoVrikLn64RwJU8sSoXTPsvv3TZHOF7n4+N4qzzXYQD6F8Ecqr686C/bezkhqhBt33UjUF2U4NcxDYw8RUAIURZF4Ls6++D77kZAachHkRnHrwK2okmrPzfNWnHf0EXIhxDAwbP6ekCRpG7ACOGyE3BJp5IvNHc/mG55zy/JCDLZBuNP4PTkGXcdXFB2ayfDUWBJdCNZ3Ryvul2NO9t3T+0DxG/J6C9Fu8EiBVVJ2zp2Uf+uxb7F9ajufPP+TRHylBWX1vd4R2eDIJYpqGHSTkOdKooQT2s7hfWd8nI/few1v/L83ct2l183rqWwx2JGXD6FwaTcOfSxu3zlOf2ekLpPkmSruoR/C5C644N/BX/l8NY483aSopSoKGUSoy/4zX8zzpwN/4tu5nQy1KyRuuMRVXJVU+lv7aQ20siK2guHUMAW9wIldJ9IebGff7D5WxVbRHe4m5ouxP76fk3tOZlP7JkJqiEQ+QX9LPyE1hKZrFEWR1kArbYG2BZUWzKtDkCRJ/cCpQIXQUpKkNwFvAli9evV8NutsZYHqPQRYcUFCbRDtMX5PjrmKWJPZWgjz5drsLSM/aBBxWSldi/RAagJ0HTGPYzieHue7T36XM/vO5LL1l7nuec1pJ1G3fpXzKQjEzEUgDEJeFmvlGSuexau2vIofbPsBX3joC7zrtHfN2zssBIQQTS3qyWSOjoi/4hnr7/LNZr6jF+Y1nZ2jiaZPu/7MOPz6rcYfkgzP+nDVNhds5RbS4Asxnh3ihv3/w7vvfci+1SN0nrPsGbR3HcfFay5mQ9uGRXuimzdCLklSFLgReLcQokIFLYT4OvB1gDPOOOPoZ3PmC5ZoJdgGUfOYmBx1LZoKjqrJ0WmqfGoCwl3ua9EeEEXITBuKz3nCFx76AulCmvee+d6Ke7ZopcY4gMmRO2SceqAFtGyF6di/nfFvZLQM33zsm5zcfTLPXPXM+XiFeZeR379vipl0oarIoBxjiSyPDsyyqTdWQUirET+XjPwQ7cjnovC2Hln51I8MAt66Cnb83puQ2449cyPl9bqXKWR4f2Eftzz2WgCWRZbxjlPfwXPat+L7wqlIJ50Op73erEuQzGlED8HL+khhXqxWJEnyYRDxHwohfj4fdR4zsEUrrUZyBF+kgiO3UJIbzxNH7nUxPWGkTHMi0m38THn3ay7YPb2bXz/1a67YfAWbOzbXLf/EwVnPSI5yIengyE1lJ7jEKwCKrPDBsz9If0s/77/z/Tw1c3SGKm029kmuYFp1eMTTkaoIyeczVtWh1NU68SD0ngAnvMRwOPOy267DkTezkfxl5zgP7J+y//7b0N94YbvCLdoEqyObeeOm/+CPL/0jL1z/Qvxt/UhqyBWb/MBUmrufmmwooujRhkMm5JKxwr4FbBNCfO7Qu9Q8juoYEU5lJxjcb3LUVWS+jsb1MJ3KoyU9OHKLsKfKw47OrSNCCK6991rCvjCv2/o6zzIlGTkUijrDM1keGpiuaFsupIwN0IRuRazLzlaMkyqrfOYZnyFTyPCvt/0rs7lDD3d7FM8sb5k0c7da8VpHc7aAETqt048Zzjc9Ww2TV4/MPNVk5KU+Nd5kQdOZThUo6kW+8OAXuOqWq0hKcHlwDe/a8kVOaD+nxPnLshmbfLf9fCJrBqGbL/n8YcR8cOTnAa8GLpIk6WHz3/Pmod5jA5kZkJSSeCDai5YY9Uyf1uzpsq5nZ9kqeGDfFHJ6AiKd7oLBNuNndmZeCNddw3dx78i9vO2Ut7E8utyzjLPv1q+6hxm+JVqxmc9AiZB7YXPHZr588Zc5ED/AB+/8IIXi4uOunGiEkFWI5hzjeKiyZ1tX0QSlEAjCiT2ohSSsOAO6Nho3JndXzEm7/rmKVsr+TmlxrvztlXzjsW9wyeqL+fOBIT7YcS6y5PECne4kE1YXFpP5poVDJuRCiDuFEJIQ4iQhxCnmv9/NR+cWGmOJLBMLncorO2Nw49YsifaQnR5m12iyVKZJGblW1LnlyVEOTKab6oqipZD1QiVHbp0W5inO93WPX0dPqIeXb3p51TIlqUADViuBWGmR+WsTcoBzl5/Lu057F7cN3sYn7/nkIZ3YjubTnq0kLxvDOXPkHtfsuN9VtoSiLjwd4VonHzV+WXG6bSGlzx7kT9vG2DeRcrRZu69NqX+0Wb6642p2TO3gPae/h8+e8zHCQiB8Ie8HYstdp2NrM/m7JOSLGY8OzPLwgZmFbSQbL6UoA4j24M96Z06pZoVQjpyXB6kHyuejP2fID0W4nCM3+zcP4WF3TO3gruG7uPK4K2taAHgpO72gFAyrFQtF6/dcbZfu153wOl5/wuu5cdeNXL/j+rr93j2W5JYnR+uWOxpRPoYLQYiqMcz37Jnk9p3jFddjM9spKiEjWFy4C2QfhZkhAPZPlRiQel1tdCOdzc3y1R3/zmjmAF+9+Ku87oTXIWnGqdci5BWbUbTHmEd5oz+NzsmjEccUIT/anDoBw/zJ6ZAQ7cWfn0EqlhQ/5S4R8z2RrHp9JiHXQ2WE3B8zrAsyM3UXTjqveYciMPE/D/4PLf4Wrth8RUO9ElSzljB+lkQrZcrObP3YHO849R1csPICrr33Wu4YvKNmWSeX6NWPoxF24KgaZZqxBvF613qZeLzsvYWAYGqATGSlIZORZYj1ISWGAXeQuLoy8rp9FqQKKd7+p7czkjnA6zd+nFXBU8kWisbaA4RaxWzSsiIzlfwWR340f/NqOKYI+ZFGJl+stAEvZIz4xxZMW3J/btK+dKhRBhs9nlocebGcI5dlgytvgCP/2+5J7txdyYGBYaly59CdvHbra2kNtHqWqehbrVUjhKnsLHHknqKVKuOmyAr/7+n/j03tm3jPX97DY+OLN0hSrbkxX+Ifr3lUSuDQyPMlhFKDBiG3EFuGnBw+pP55QdM13n/H+3ls4jFete59HNd6JjtHEzx4YNqOsW9x5EVdcNAZ2M42B7YIufHnYhStLD6DSQ8cLcM+MpslWyi6cyQWMkZSCQvm5Alkx8mFl7meb9QLr9H3rRCtZE1CHvSwFQ+2GRx5A/U6lZK6Lrhj9wTH9cX4wbYfEFACNWXjFhpx0Zf1PLLQIBB1yMhjgFRXtGIh6o/ylYu/wqt+9yre9qe3ccNlN9AT7mno2cWKRr5ho8R/TnHNhSCUGmS652mli9EepIndRvTNYpHk7bfjW7mS7M9vomUyjiTpzJxxMsXZWeRggNiOAZSJMcZXdMHMDNrkJNnHH6fzTW8CoPuXv0Hr6uHm8ACdmT18/JzL6bpvBvx/QBkfRVJVkse3IY34yTx+gGD+Noqd3ex+skjPM05B8vmQw12Gt/DkAPKqswAJCvmmCLnQdZCkklhU19HTaSS/HwnQpqYozsbRxsYQhTzpe+6l6x3vQIl6hA04BBwThPxwQ9cFd+2ZZGNPlJ6WoOtehZ1wIW0kcLBgceQOOflCmx+Wi1YKwQ4q1D+htjnJyHOaTkHTeXBwkF8/9Wsu33D5vAWvsgJm4S9x5EiywaE3IFqx0BXq4ksXfYkrf3slb73lrVx36XVEywMpmWjW6/JIQpT9nLcKHbA48oasSgoF1IMDaA/vYOYRidT923nqh89H7e1Bmt5PZu8sK750MQADjsesrzv885/a16y0yOUq7bFPfxoAP+Df+SQnAicC/PWXFd0x2uiC276J8wzqDny8HK7/IEgfgnCYFakUuVgLu6IR0HW0UUNvonR2onZ0kNu1C39/P8VEAingRztonDKU9nb0bBaRz0OtUCGKQvhpZxG76KLqZeaAY4KQW1Ns73iKsF9h+3CCZ2zqnpfAU17QdEEmX+SJ4Tg9LcHafHQhA06tucmR+zOV4on5dwhy1+PPTVFUAuiKh8zQ5MibhWay53eO3kRez/Oq41/V0HPOMfPauASgaqbc2iFaEQhTDNScjfiG9g18/sLP8/Y/vZ23/eltfPXirxL2ecT+EG7uc1GcsufgGdzoe9kcedk1bWSE1D330HL3I/j27GSfXyJw/wP0AmkgTQswSB7QEwn0ZBxZ0Ulc+Hxkn0q3niF86mlM7drD6IYTaQ8qrOlfhtLaitLawh3DOZAkzmyDSF+PQSAVBUmSKIyN8Y0Hf8Ptu7/NGfJ63njCPxPesJG7B5PIiThyMgF6kVOCQ2R/+WmSbZcwsupk1JEhKGqs3tyPNjJKcXoKac8fKSirUNdsJaMGSE1OE1q5gohfJpvJof3+t8b7+03FvSwjhI5/9WqkYACRzRHcsgXfypUG4R8fB1VBaWkFIRCFAuHTT8Pf348ca0Ht7sLX25hXbzM4Jgi5c07uGEkYJlFFnaAznsg8wlroxWJ9ZROFDOM5hcR4knXdUduLMuAlI2+w/ZrH4nzaMzgRgD8/RSHQgeb1eKgNZgc8btRGURdoep7bhn/F+SvObzqqW61XUS2O3OEQJARgBs5qlsaet+I8rr3gWq65/RrecstbPIm5LgTyURSzp9am3vQmM/kUbP8NnPN2vNRj3kpngTq4H1mRmLz916QfeIDkrbfa8jVrixUnn4S+9SSy3cvoiqXoyf2Cx879PGe87BXIqgp//Tzc/GFuffHbUUItnLbJWAepyTTZ0QS0BYkuL+lVhGkW6F/bheJzr+Mpf4Hr7vkxXZvWc+Wzv05Hm0EYNdxWR1FxF9EdSZTXvovdmZI/wxnO8AifXg9bToMXfp59EykOjCVZ0xlmeW/MsGJ649UNh1M4kjgmCLkT9bTg89JGMwuokCaLnz3jKYOQqwHy/jb82RJHblfXoPlT+W27/MGH4evPgJd9B054iYf54TT5QAda0SOmhCUj9+KOa3RI0wUPTt1KvDDFq49/de2OV3mbaspTxeLI/VH3SwcrA2c1ikv7LwWoSswrxvao0cDMA274Zxh+GLq3wIZLqpcrFIj/8f9I33svMzfeSG/WMOMbA5BlfCtX0vqC54Oq8vjmp1HsXc7FJyzj4YEZphM5Ngx9l9BfC6S2nI2kmETYjOHjy89QDLUwMJVmWWvQkRe2Qf2QELzvzvehiyL/sPZq2r30PfZ7GIpNXa1iRw6GV3PaYKqqe5ca+Q46In6CvoVhDg8VxxwhPxrWXSZfJOQ3P7iWpai45ej5YGdVW3JoxOSqyo0nfmH8vOdrRnyLMvhyU+QDHeyZSJLOFTlzbQetITNwv2W10iSbVyzq/GXkRlZG1nHOsnPqP1CGnKZ7enSCQ0YeaEEyHWF1ISj4YvhSI023ZaGcmH/pWV+y76VyGsmcNue6DyuqfCrPzUeIkjv68COIMkKeHxwk+fkv0XfbrSiJOEPmdSkcJnfCKWgXPIuTX/Y81GXLXDqEosP23trwldn95APtFH0OPUTYJOS5abLaSnaMJJhJF4gFa5Mg53RM5TR+/tRPuWf4Hi5f9WZWRjaWwhR4zVvT/JBq5odgbDDp6Zp9mM0UePJgnOVtIY5f3lKz7JHCMUHIXRH0zEl8OGWb5Rv5wHSaTb0xoxOFNHoFIe8m4OTIK9yW65gTVrt/4G7j58GHQMtVLGdfbppUbC3pnKGMyRaKJUIeiBnxMIqVnq61evPwxP0czOzhDce9tykloVTHZlcIp4w8CiYh3zWWRCkE6XUoO+cSc/vS/kuRkHjv7e/lTf/3Jl6x6mNEfa3cv89Y1F2xQM3+HSnouiCZ12gJ+po7LSRGjETEAOPbEbpO5A+/InzbH9kxPICeMHKj0tZB+rwL2fiyFxI5/3zicoAH9k8T9iv4lndVr98BdfaA2/QQXBy53aVsA4Tc8Y4/f+x+PvfkZ7lg5QU8o7eBTFMWR+4LUnUWhztcgbPwKGkpfNP5o3eDPyYIuROHY+E1vIBMz7Ki4j7a5YNdtEw9WqP+Ou1XI37TexFKALmYM1K6yatc9/0mR+4JM4aJyMWxnXVsbqfUcM/AH5jsPQ8jsx/csPtHRNU2zu2tcVSfI2wZucPCpFgUaL7onEUrTjyn/zkIBO+/4/18Nv5mXrDqDZze+SyjnfkMITgPsLaqHaMJhqYznLehRFRry9HNe9N7EQLSE2HiP32A+KfPpS1ubIaBc88lfPwW8pe+kO1yGwAtllw4lTfbaBxK8iDZyHqzfVNkYQY7s78pUNAb34pShRTffeo/CClRPn7ux3lor5uoesr3C2kkQKghDBWsB8IdMHhfzbaP5jANFo5ZQn4kx97mD20ZnZsjzwW73OaHjr6q+Vmgu2b9nq+WTyElRxlfcQk9QzfD5C5El4MrKmRQtDQFByF3jZGlUMwlKamwyjD8CCfd9U4m+p4Op/yGgfgAd4/cybOXvxq/EqjZ57lAcVmtlDxhNX8LUi5uv8ChyLEv7b+UrmAXb//Tu/nBnmuZzo/xrL4rbVvihZpGQghymt60zNUKsapVk0dBRaczTzzB1Oe/SPKuPvSCjCRniVz8DAY2nkr21DPRO7q4+Phew1nGTNpsmWHO5f2V1Bj5zrPdF83NWC24PWgtIjmZylHUBUqZpZk1R7/52DcZzR7g9Rv/g85QJ1iKzRqilVwmRUBWQfbO+wmYopUpEKLuye5otko9Rjw7Kz/iUaGkMmV0lRx5N6qWNuKIOBC6/3955i/PJDpcO7t3hSgGYNpInjzRd4FxcWKXewTMELVVOXK/k5A76nVi6AEAYtNGFr+f7foZsqRwTvfza/YXjDR2tzw5WuHeX120Ijw5cgDNF0MSRWStuaBh1XBG3xm878TvsCF2Cr8d/BYffOgljGeaj7uiFRuLgQOwbzLNnbsmKo7rQoiGA6LVYlZ8O56g/b//g54XX8S+l76M+O0PE2wv0Hv58ax77hjLPvP/SD/ruegdbnFJz8Dv2Hr3e+zTZDUksgVP4ilrWeR8nFzI8JfYNZZk52jCNiG1N2frHcyfuYLOkwcrfQMEMBAf4Mfbf8z62ElsbT27ooyzHifGJmcoyoGasYlEuBP0AuQS1ctUvWO2k8iSyNaOsjmbKTDbZCz6ZnCMEPJKlM+xsUSWeJ3Bnmvd5bB37oIlWimXkRuLx+mmjxBE7v5vALp2/JjReBa9yvHe86oZxyLVugEt2Gnk5nQibRDyQlXRikks80mPm2aLZoJmWc+TzCe5fsf1nLfsGbT568tOB81ASVnNkM/boWtrDKaipYxcnWUxVDWfSRTy1Rdfs4ioLbxl86d5es+LSBcTfOC+f+R7T32CRL4xx6OxeJbbdow3PMemUoYuwkocYUEzv/lTE8Z3aPZkmbz9dmbfdhU973sb4Tv+BIpC97/8Cxs/9SLWPGuGjpc9D3+siPBIbiJ0jZPuejfLDvwabv5I1Tbi2QL37Jlir0d8mkDWqDcXNAj5wFTa2JQsjryckDuVmR4yaCEEn7r3U8jIXNn/b1X1MF7jJBcNQ4Ph2eqb0njR8LDMxmsYH9T5Bo8OzHLPnqmaZe7bO8V9+2qXORQcu4S87O9HB2a5t8pgL5gMzOTIK0UrhujEMkEUAogPIedmAAiOPshjg7PsGE14WlCUEz8hhC0z1nytFGIrYeaAO/FAytg0nBy5qx7Te1LyJOQmTELuKyT4y/4/kSqkuHxdveBY3hBlP72gFlKGPLwMFiGXLTd9q5L9d8EfP2BvoM1ClmResubtfODE77G59TQemrqN1/7p+Vx7z7XE6xD0sYRBmJPZxhRi1eLrNONTYI9dLsfkt77FzrPPYeBNV6Hddy+5LScy+rlvMXbDzXRd9SZUJW54FZshZaVk5YnDP2Sk2s0FOuG+b8Lg/Z5rI2tmK0pkK+XU1pzOhcrEg6ofofgrTqH1+N3f7f01dw7dyVUnX0VXsDK2vXXyLnichuRiFl0JotSgctPCmEta0lvU6WzjqMwJbOKYIOSeu/H+O+G374HMjOdHrvf8oUAXBnGdmjWIa6VoxXAYDjhNEMd3ADDdfSah9EGkYp6h6Qx3PzVJBSomGjYhL/hjFKIrK517bI685D7vIuQ2R17iciv0DSYhB/jJ9h/R39LPiV2nVPavCVTbRAWgaEmKamVMihJH7iCuhSx8/0Vw15dKZphzRFdwOW8+7lO8e8sXifla+NH2H/GcG57DW255C3/Y94eamXTKZbzVYA9ppuBSrNaLNuhsOnf3XXR+4r30vuHljP3XZ0BV6bzqKtr/fAcTn/wiWv/60lRJmSn+HHljy+EbfwKAB5/5fcNB7NuXErnlmoYWiDUmAdNjOR+s1PMIf7RCtFILRV3jB9u/y8b2jbxqS22P4ZQHw6MUc+hKZdJqJ/L+VgDkbHPc8mQyd1RlEjomCLkXAv/3XoOruO8bdS0Q5jva2YHJNPfvn2Z4wgxS5WF+CCU3fYGACSMCxPjyi5CETih1oGr9Wtn7CIEdJ0XztZCPrTCIrpkqRgg8ZeQufZl59JVcMvKycYkPoalhdvt8PDz1JC/b9DJkDs1Botqn0YVALaQo+kxC7liLBTMCokXIBcCBu0py3YMPHlKfrDrXRLfwrYtu4ivP+gpPW/Y07hy6k6v/cjXPvvHZfOmhL7FrepdNwKw51mimG0tstms0yfaR0oZkTcWq+0E+R/DeOxl7/T8x87Y349/2GPnjTmDVN77BpjvvoOdf3g0RR2o8a4Czs4bTl53Wr5JBUGf2UPDFSLWsR7z5LlhzDsGHv0372N2Vg1MGq9+BjFu04irji6Jq7hNfraV3x9gv2Rffw9tPeTtSHVLlNY9sjrzGNyn4DcZGSnswTGV9dFbz0IEZ7tpTXRxzuHFMWK2Uf0NfdhJlYrvxx5M3Ic7+16aenw/MpgsEcgb3Ue5Zlve3IyTZ7RQ0sQM90MZM1xkAhBP7Sbds8KzbU3mTnUXIPnQlSD6yArQsSmYCiJDMFdCS48iSiuYrOTQUPTnyStFK6QgfJx1bx8+lQVRJ4YXrX0i+wQRL5WNc4mqraTsNGbm3aMW4JhccMvLhR4yfHeuMRL+HClvEIfP0lU/n6Sufzv74fm7efzN3HbyLrz/6db726NfoCnVx0aqLWBE4nV7f8ShyYwHDnG/tFMdUzcijF5m+/npaPvc/yDPTFHt6CL/hTew57Rnoy1YQ3eLtRm4PczYO7f22PTeZKSg77CizB8hEVoEkIaI98Mrr0T+7hRV7fspTq86vqNuLPvqzYwhZpRBoq+yLP+ISrUhUX3uz+Qn+MPRdzll2PheuurBuEmgvZkzWcxSVYM1TUtbXav5S3Zy1GiNYy3DocOOYIOTlaB+/1/hl9TmG23odjnuh4g9LtrKzzDRPVsgHOm2nICGAiV1oHRtIR9cAEE7uq1pvTis/0hkycj3QApJEIbYCACU+COpmdB3GR4foCLS7Vp+XjNxJyCuGJZ9mpmsFv5ameVbLBjqCHQznMjSD8jqrc+SGzXEusLLinmYuPtVMriwEMPo4omU5ia5TiY3cfei5Kj1IzJqWNbzhxDfwhhPfwHh6nN/t/R2/eupXXL/zeuB6fHKAp49dwNNXncvl6y/Hp1Q3e6s250qiFfPv0RE6P/kJtIfvZUTXobOLqfd8hK2vuIyirFA8MNPYu+bihi23PwKKH7LT7vcVAiU+QCKyFjDk4Km8RMuG59K587c85UgEavV8LJ6ruBbIjqOFuo0olWUQ/mhNZacTNw18naLQePcpVxuikWoiOKceqAxKMYuuBGouf82Kb5+Z9vzmq7d/i/AAsOr11Ss5CnBMEvJI3EyoetwL4MBdiNQotV51oXSdctEg5LpSGesh77AlFwATO9HWXIQWaCPvbyOc2Fe13nJOwBCtzCLMZA4506vOlxiE9s0AqJkpCkF3QgmXVYyighqqboZVNLw+bwvIzBQVXhLxPi3URbl8v4ZoxbBaqZSR5wMdCElGdUaQHHmcdPsWxvyraYnfWDN4WEPdtPvl3cHucDev3fpaXrv1tYynx/neI3/k3tHbuOPgbfx58GY+dtfHOLPnfF699QqeueqZFXJa53tLrs3V+CnPzjB67TfJ3fhzgskERCIs/9jHeGLTGWQ0kFSfEQvbq+9efc7FDacvSYJQO1LaLRPWdYOQZ7oNzvvefVMUi4KnrTyb4OM/IjyzC6gfxz2QGacY8S4n/FGUstAUXn29b+Q+Hpz6M89e/mqWR425PBdvZ7mYo+BrqWmKLCkqBV8MKVPa2KzyocQ+Nj36/wBoDZ0G0ad51nE04Jgk5IHMCHq4G7lrEwDyzAGgelS+Zgl5o+UVzeBWy5WdYDkFGYRIzs1CchStwyCO6Vg/4eT+6u2Xxy8Hg5Cb3nO5qMGRq4lBME/6/uyErWS1UM4N6/4owkHIbacYAZhH4j8ywYqCxtlKrHSvAZQcbKyf7ute5dUqohVkBRHuwpc2CbmWhYmd5E69mIzP5OBnB6B7c2Od84DVrUYcPLvD3ZzX+1xObr2ETctUfnfgZ9w3uJsHx//MfbfeyfGdx3PFpit40YYXud7PiXTeiPEiDw3S/pnPEXrwbqayGaR165l4y7/R+9xLaF3eDk9NgFZEIBqft0IYG7Q5Pwh1oKfKCHl6Cl8xSzZsWIZYkT3zK434OS1j9wDnucbG3YRx0Z8dR2tb7d0NfxS14J7Xrg0NKOgFPnXPp+jw9/GsZVeWyjXwiuWwZOS1xklCouBvRclOV5TrGrnd/r1t4n6S/bUJ+eNDs/R3RdzB6A4TjgllZ/kHCKaHSQZ6od0QU0gz+2o/75gmhsfd/GijFZMjL5aZH4IZb8W0I1enDEWn1rERgEx0NaFkdWWn58R0cORFfwsEWvElS1YmvtwkhYCbkDtlf0IIsnKI2RkHZ2IfWwXk0wyoKo8WR7k8nUeuZaboASu+S8XmUePIrBaSaL5KjhxARHpsjlyZ3AGiSL7reHLhPqOAw8LmUNAosbTKxfwx3nTiVfzD2qv58Mk/4nVbX8dgYpCP3vVRLrnhEnbMPoAu9Ip6737wKfZ/9BOMv+iFhP92K9rWk1j3m18T+N5PyZ7zDCR17idKLRM3FN9WXPdQO8kZd8RJkfZ2GNNbVpMPtBOZ3l77/c2fvvw0WpWIhLovUukQ5Oi7AG7ceSO7Z3bzotVvwS/XFos42/UqZ1itBCrmmNPxSJZA87d5xuJvnXyYbKiXQssaV0gNJ/c/5EgdNzKb5YmhQw8dMRccE4S8HMHMKNlwH5icgTRTnbsFN3HZMZrgjp0TTXnpVYNcNF30HVYrlsNINryMQGYUuZhj5oBh9mUT8sgqgpmRUoLm9BT85B9hr3cCYUu0opuEXAgBbatQE0N2mUB2klzQ7bjjnOBCYIgxcpUycl0HCml+EYsgIfH8nGpn6GlWKmUnDLY2iWq6Ti2HrBc8RSsAeqQbn2X1M/I4ANmO48lb5pWZaYQQDEylqzpWNQJrjNJ5jZ2jifpHfEpjEvO1854z3sOfXv4nPvC0D5DX8/zvzvfyyUdfw2DKCNQkpRL4r/8+vW9/DZHf/4LAsy9l9IvfI/0fnyWwwSm+cotmdFFiQMLxPUaIWjP4U3kXC2mTuJjxdAh3uIJXAQhT1KKZ5nj2dQnS0X5CDlFfVVGFEPhz1Qm5YX5Y3WM1V8zy9Ue/zmk9p3FC27lmW3bVNeHVJ8MhKFBx5+BMhgf2mycSCQr+VuTsdEW56MwOEu1byfScQuvkI56B3rZ5eKMeCcwLIZck6duSJI1JkvT4fNTXLMo/YiA9TC7UZ2TmifQg1UmY4Fyc46ZjR7mJX632qkExJ5JT8WM5JaVa1iMJnXBiD5H4boqyHy1mbDyZyCokoRNMm8lqH/6hkRDgR68AUXmkFqayUwRKFim0rsKXMN5b1tIoxQyFsgVmEampVJ7tIwk0NeJSRlnvqQuBlp3ll9EIpwTW0aHEIFfJeaRyGrc8Oepp01tqs/xv77G0rBss0Uq5FYce6cWXNmyhswOPUlSCFFr7bXMy0lOMxnPsGEmwZ6K504MTVvceGZjlwGSa3CM3wvdfDDv/r8Yz7ncKqkGuPO5Krn/B9ZzYdh5T+RF+9OerSF77Tyx77eUEvvO/FPrXM/bZb7LnzdegrfAWTVRrY90TX4DHb4Q/f8KzrGYRclu00o5aQciNk1iFtYmAdHQNwUSJGfJaGkKArKWQ9QKaR6o/IQTCb84vJwPhWEu/H/gB45lx3jUyaM9De+Ovs+a8RSu5qqKV6ZTBUFmiFcmR5lAI439GAunVZDpPIJgZQbGU67X6UbOXC4f5EuZcB3wJ+N481TdnKIUkvkKCbMg8Ykd7atqIgnti1gucM5XK8+D+6ZplLMhatiKErYVki8F9R2Z3E5vZQbp1I5gZjTJRI2phKDVIJraGwvY/4ANDTj1zAF24uabBqQybHRw5AG2rUffeAULYjkf5MtGKNcGt9+nxRVwmkc5xuWPkHsZVlVfGzkLzzXrmzByJG6KkwekMiWyB45e3EPa7p1ij+UmtWORVOfKWVQTTo0jFPNHZ7SRbN4GslDjK9KRtXlkr1kY9ODezQHqYwO/ebMTm2HcnvOMB+9RXep/qJGd5dDlvCryO0E3jRLZvQynu5bEtMZa98hpyW89rtENlfws6Rv9m/L7rZmOjLyuiW2IDa36E2vHlZozBN7lMS/lZ8LdVNJeJ9RPY/0tbgex1KhEIVNMSRvMIAyEECF8USeiG7No0ybWqms1PcPPQz3huMsVp4/ezres3DK2/0vW853DYhL4SltVKLfIqScY7S2VWPKQmUIoZMpGViJZ+AALxvcCaBbNyOxTMC0cuhLgdWLhAAk0gYAY7yloZ6kPtSJnaXfOcmFW+1cBU44GalGKmwhnIQjq2Fl1SiMZ3G4SofYs93axYzqHUAEohiTJwN/ry04yb03srpqVczCFpWcP80ELbKpRCErUQt5Mu58tEK+XQVLedrzVhdSG4ceg2ujWNE6InGVxyrhQlrxwTyRwz6QJ7xlOk8waXbkO4f6m2JCyOrJqMvNi6BglBKD1EdGY7ibbjjPpklWKgFTJT9pZ8KOtueCbLTDqPENAz+H9IegFe8ysjdvv933a8lltMVY7C6Bij115L77v/iZYnnqRw8pn8/BOv4doXa3xUfJ3RjFsnUs+kUDdlOMH0EP78DPRsNb5JwiPhhqXADpZEK4qet62qANscsYKQC0jF+o0/qohuLChmHcUqybeF+S29xCu/H/ouRV3jHdMG12vJpBv9dBXEVS8aorm6yk5DnCRnZxC6QzdmimMz0ZXkW/sBCMzuNfpUo74jReMPm4xckqQ3SZJ0vyRJ94+Pe6f2miucgxdMGxPZVnqFO+sScifqJUAuNvGlDNGKd5opofjJRFbRPnYXgewE6baShUUu2IMu+wgl99M69Qiy0MifZLooT++vIJ6q6RgjnBx5q8HVB1ND+LNmnJUy0Ur5mxTLRStmgYnsOHdMb+NFyRT4YoaLfI0s9kGfMa3S+SLTZRHfyo/l1WTOllKsqHpnuy+2Gpxw2/j9+PMzJFuPK7URbIf0VMNelvVw/z7Dxrh76Bb0ruNg3TNh06Xw8I8rVq5BX0vX9FSK0WuvZfdFFzH13e+hdfcyec0nmfnwp3na8f/MP2/4GOlCks9veye74w833Cerjdi0oV/hpJcbPz2coSSLkAdKohVwJ3kgM40uKXb4A2c7GdO3gem91fsjQDXd3D05ckD3VxJyIWD77H3cM/F7XhjczCpNg74T7QibzSqbLci6ISLVlWDNzUCSTNGK0N1xhsygc5nwCnLR1QgkAvH95rscoxx5IxBCfF0IcYYQ4ozu7trxtg8FgYxByLOhPnRdGGEq000Q8rK/hRAumW8zijNZy1QEzHJiqvdc2iYfBiDVscXxoEKyZSOxmW3EpkxF6KbngaTAzP6KaaSaruoWRy5JErSZ4pn0QdvxKB9wc+TlRFTzuWNhWPf/cvC36AhenDCiEbo58urvrxX1ivHyNJ30gBXCtmCKVsqfK7YZjiu9A78HsDlyAD3YAenJklPNPLBJcnaG9on7jO8AsPl5kByB8UprDiFAnpmi9Wv/zY6zz2Hqu98jfPrp9N9wA6Nf+ynZs87Dmmlb2s7iI6d/jRZfB/+78xrum6gue69oBwgnTf3Ppuca1yZ2Vc7RvKXsLFmtAPjyJT2HlJkynGM87N0tk8R6lkCqefIrVJORm7lRlWLJ0qMoivxp+CeElRhvLHYZm826C4nGd4FumFmmchr7JuvFaHG/s2L7cATqrlnrFKKYQesAmyPPRlZSVAIUAu22TqYmR36EiPyitVq5e88kDx4wjnJujtxQEOZCvdy7b4q9qYAh/xLV5aSj8Zwd0c2CVedjQ7Pc9dSkHUe7mcwxSjFXlSMHOLj2ZQDokkK86zQXYY13nEDL9BO0TD9OOrISIt3QusKMaujNkbtl5AYXFUwNEUwPo0uKHXWxGgyOPA3mEXMmXUAXOrcP/5anRVaxStMoKqGqHHn5Jqjpoup4laxW6nHk1axWeigEOugcvZOiEiTRcWLpPYLthmil0WTWDRD6yNQTSEKnuMZ0VV/7dOPn/r+6yhUTSUbf9U763vhyon/8Fb7ly1h93XdY873vEjphq2fdXcFlvGvLF1gXPZEf7f00vzrwNYd1TzXhsPEjmB6i4GuBrk0gKYwf3MuOEbdTl5wtE614EHIy0xViFauZgr/NmMf1CLktI/cWrehm7kwnR37f2F/YnXiE5698PS3JA9CxFjo3IOsFghljLT+wf5rBKW8PYmt0Kjhy0+KrKFdarTghmVYrgK3wFAKYtvKORhAYEUtdDmhzxEJFWl20hDyZ1ZhK5iuuBzIj5AKdCMVPMqtRCLQjCd3mWr0wMJW2czXamezNe5Ybckle3Hgf5WLGVLZ4I9FxAo+e8wXuefZNSKq7XKL9BHz5WXoH/0Ci/QTjYrQXkmMVk1YtmNyxk5CHO9GVIMH0EMHUELnwMoRcrnh0QzOJppNj2hV/iInsiO3JaXDkMSO3p+YdaMXqX1EX3gG+GoCVSaZQRbSCJDG14cUAjKx+oWucLdGKpbiux5E30qfwrJG4uNhpOJnRtsYIQjXymNGdRILoL3/C1PMuIf3nP5PffAKjn/sWG/74RyJneydDcCKkRnnTpk9xYtv53Db6M67f+5Wai96IsAnB1EGykeVGzPZYH9rMwcrCubhhOWUGRhPBNgCX5YqcmS5Z/DggTIVoLrIcZg6Y1yqbEAKU7BRCkktu7877UOLITUIu0Pn+zi/QE1zN2d3PJZAcgPa1RrwcMJziRGPizIpTqjA5crW2jBxKhNzNkZfyjgohyAe78ZkBwZo94ekuf42mHm0Y82K1IknSj4FnAl2SJA0CHxFCfGs+6q6HsYQ79nQwPULOUnRSCtvqz0+jeQTysWBx5Nbir0yIbPxs5iMqxSw5X+WkdmJslZHRvdyZPN5xkv37+PKLjPhGkR6Y3lcxaX22aCVmJyk2Yq6sJJQ6iD83QTayotJssezvooNjsjKg3z3+O8JKjGeoXWaZUMnbMpfw6HlpURV14emMMZXK41cNHqJa4CFLtGJHP/To+9BpV5OQWxle8yLXPT3U7hKn1ftijXzRcPwpCr4WMv5uomCwcn0nIoYfZfoHP6TlS19BnplCPelkWq+4gkc3P80Ie1C1zVKr1slBlX3804aP8NN9n+PmoZ/x1Ue6uHTFa2r2NZQeIh1ZZSTni/XZ0QedkPJxCMQYTeSIBlRCwXYkyjnyKQp+j9CzZoO5yHLCDXDkBX8booqqVtjzy2AU7hz9NdP5CV697v3IkmL4BcSW2YQ8lDgwZ2WnUnTIyGusWYmSaEXOTEPQ/DYz+8lEDL2VwIivHhuvrey1C5fBuREtlOBlXgi5EOIf5qOeueDRgVmWt5XEF8HMCOlIKemw9ZF8uWmIra1bn1TGkVtwmqE1CkXLeHp1VoOz6kTbFg5seDX+3BSjq57HWoBoNwze68GRG0fnor+1RMiBfHQloeQB/NlJppY9vW77xTKOKVmY4bGZv3Jez2X4tJyRxFaSSwqx7CzCi5A7+uc1XiOzWVZ31o6DUk+0IhAUlQD7try54p4ebIdCCmGGtZ0X0crsblIt63lkcJbzN/oISILpHSEm/7gLLf0JxIpVTL7pX9j40ucRDflhd22T12qQJIkr+t+Nqgi++shX0bUIx4UvrewzBkseTA0x1WNy/LFlBIa9lJ1xRKCFxwZnURSJC9Z4yMizMxRMhzQv5CLLYfjPVe8LBGpumnyg3ZNYCQG6vyQjz+s5/jD4I9bGjuPUjguRtYyxeUd7ILaMohIgnNxvWEhWbdUponNft72qPRyCLFicsmU7r1gxyfUizAyQ7XmWXXc+2GWIVmqYl1ZDuQf1QiSoOCZirTiHNpAeYbr7LPtvy9Ov3JOtbp3zEKLSiPVQ2mRquv6Xf1tJYudpH3Jfi3RDehK96Ha4scRGhkOQZge8z3SfTM/gbQAk27ZQiTIFoskxqSYhv3/yFopC4+zu5yEGfmREzgObWyeXgMAyyuEkjHMN9aloKcN0TFar5q+sZpWiW9mOcnHAX7etRhZmJP4U48svhFyO6R//mMSXv0RxdhZfRKfv/Vez7ZzLyRZ0JFk+5OOzLCm8bvPVqL4M33jyv3nlWh9ruy9zldEFKLlZVC1VUkbGlhHYU+n9K+WStrNYsSgQaghd9rnWhJyZskUMTljvko2sgNS4nVDcC2rOEM94iTwFwiVa+dPwT5jJj/OW4z+EJEm2Qp5oL8gymchqQsn96FhBxZoTj6kOq5VqzFe6UESSSmGlFTN2j5oaBb1ghPTFmM+5YA+yXjA8huUqycmr9FI/DBz5opWRe0EppPAV4iVnIMo48gZg2x436LhSsz9lduS18vrVc0QCDNGK0G0zLwtqIW7EIjedLFI5jZHZLKnl59pl4h6ZfCpEK4qx0GQtjRCCu8Z/R3/keJaF+kkn43Y0QUuWXh673CsTSzX+pd54qoWk3c7O0URFLU8MxavGmdZNGa0zScahwJebxpeZJHHnDCv+4TnMfOqT+FasYMU1r2f9C8Zov+A411Gu+alS+R6KpPLx8z5Om7+dH+y5lgfHKxNyy2ae1my4DyEEergLX37WsHV31p6P25ubT5URkkTB34ZqcuRSMY9USHvKyC3krM0iftDzmxrmh9OeccgtWMrOg7lhbhn+Ead1PpPNbacApRRxxIy46unYGsLJAxSKOoUGHLrK+6SaopWiEqyq1xqZNTclWaEY6kJJG2IpyyPalpHj8MHwyKxUD0490UI5Ex0ThNwaG9v00CEjt8QAdkb2OvCKpwBz20ktF2EL5Rnk3e02UKGZ3aWQcGvPffm4HYvcQjxbINV3JvG240lF+0m0HV+3eidHvjf5OGPZA5zdbZjbFbJJ8FkcuUnIc0nvRd3Aq9SDqjmyA3n1VRfMpCs5PzB1BWCbSNYb25oOI9kMLb/+IU/9tgdxh2Hb3POFL9F/4w20XP4yo+6pPa53ni/LhI5gB194+g8IKVE+/eAHeWrmKZe1j5IyiEou2MOjg7PsSBgK33KmRc7FEYESFymEoeCzRCsWZ+5ttWI0mIuYhHzGO5jb3okUam7Gsw6rTUtG/uvsoyiSysvWvhVrtthpD81UdOnoGkKpAwxM1l63djTNcqsVmyOvDJplIRYsxYsvhrtRTY7cnzDeMRtZYdedC5mheRMjzYtWiguv7FyUhLyaXWgwU+YMREkMoBaqxNkug7Xod466yze9OIUwZeTVzQ+daMh5JWw4Wvic2nUMjtxleoi5cCSV+y6+gXue/St0jyQHlVYrpaPv3eO/JyCHOaXjGaiKRJisLVrRLEuSKhEQ52OyKoVUqZ0qddbjyEVZ1pdmfACUiTHavvJfLH/lcwn+6AZknyD5jrcwdOOthJ9pxhePWVzqcKmfzM9GZr1ZzN/KO7d8Hp/s4123votUoWS/L5uEPB/qYTyRo2A64vhz7hOblEuQlR3p34TwJuRlcwiMVHQAOZOoMTvo+S1GZrMo+eqEHADFz+OBILdp+zm763m0+bvsuqy0hxYhz0TXoBRz+NONccCV5ocmIZcDddauMdIZfxda3KAfvvgAIJEJr7BLlTjySmVyPaQLNUSq84RFR8izhSKPVQkVGbC8OkOltFdCVhG+cMOE3FJMeGUIbwaynkdCVHXRL0eNJCil9s00Xb68m+NSCwmKgVbvUAOyWpFqzr5XxWolm5/m4em/cFrnRQSUEIosIRcyYMo4bbf5XLwBRWK196n9oKola3LkUENGbnKfB0dKRGD/ZIo/bx/zPBU5+6KMDRP70Tfp+dfXE7nlt+jhCPKlp7L22eMkz3++kQbNKq/6IdxpOAY565tPrksIVklR3nvatQwlh/jOrk9R1DWTkFv5MQ1rEyvevD/nVrRKuTgzuhXbxOic5m+zCbjFGNQiwvlwLyBVJvU2IWtZlGLOU85uQUPnE10dtODjuSte67oXyI4jJMUYTyAdNTx3ncG6aqF8Psm2sjNYVU/j/E5pf6d9KvAnByC2DKH4zXKi5IORHK35fTP5YoVOx+mjsiRaMbFzNGFHKLRgDY7lnu+UkQMIf6xCtOJF9B4emKFQJXxts+Mva1YI28Y48nr8uBCiKkfuy8dtLtQub/5Xt04HLAuRuzKPU9BznN39XLNvEpKWtpVVtiVJFRn0fHi3GRx5bUJedZMwx8L5zYdnjYVdTeHs27WNzo/8K31v/gdabvgBhdXrmPjQpxn+/m8In9GOFoiVlLxORPsgUdowRBlL7oxXXQ2SBIriPQOi26/nGTedw5nJWa458xqemLmHH+39tC1a0XxRe6O2FftO0YoQyPmEK0FHNdGKVoOQC9lnmAbODnp+3VId3oRcCPjL4M084Vd5A6sImactqy5/doJCoJ0iMrc8OWqHBahLyKtYrZSyczXGSOWCXcYGKATyzAGmHUp8gTHndSXYkIy88jTv/ft8YtFZrXhyVLaMfNh2BnLdD7R4EPLKuicSOZa3hTg4k6EtbIgi6nHKRhkPZVWNpBJzhs2Rz9iX/KpsilZWu/o5lcq7Ig82Yj1SVMPowB9z21gRXs+qsOH8IkkgF1IlqxWTcyef9CbZDUzWuspOLUXaSXw8Kq2aXchv6UUqT2HlT+T27GX8uuvouf5649loC9Nv/TeyT3u6LWcLpodLliHlfY/1QmLYvVgdrdQK6euELEkUPd4xsu8Wo5ntP+UVr/o+9w3u4Y9DP+I3+2/gn1KjLm9dW7TiUIbLxRyy0OxNSJj9K/hbbWWnRfjzgXaCPqXCy9l+p7ZVVTnyWnJ2gFwxx/e2f40NGlxCmCeteoX1/CxaoB3dZKSyoT6KSojwzA5Y41llWf/cULSSi34jyAe7kfUCan6WQHI/Uz0lQwHDYlBCC3bgz0wfkuxsoaxWFh8h9+CYrcGxnIHKia8IxFC0xpSdltzV+mnXVUvs4UFQmuUIJKkBTtYfpqgEPAh5gmwZR57OFe2sPNWQLhRdpxtdCXBfKMSAHucfe99qb1AShiWLRcgNL8FYDY780KE4rFaqoZyQW3+VE3IJhzWSAKFppO+9l6kf/ojkn/5kmLud8wziV7wWbU1lSsBgeshNyJ03o30wvsNxzx0vvl6IgnrwTxgh/v2jDwHwwtX/zIHkbr755Oc5OxNjbbCUH7Pgb0NIsh3tEkpj4AyGJYTBOataGqmYx2+K6gr+dmI+2ZOQA9C6EoYe9LxlcffVrFau3/kTRtJD/E8miM9nnFIKRZ1cwVjPan4WLdCGzxoXWSHefjzRyUc96ytH+VyQHVYr1SAcGg1LmRlJ7CWYGSPt4XOiBdvx1wmJXa0l+7cFYskXFSEXQngSJ2twAplRMpGVlYTc34KanGmsDXPQ4/lZrnv8Jr7/xE+JqK3IO3Nki2mmMrP0BldzUsfTObXjQqDXsx7Lzb1xGXlt4Yr1OlqgHb9DtCJh2JG7Qtg2iGJR8MhAqS4kiZujMVQkTmg/11VWLmQQfgdhDZQCZy0E1CqJl52ouiQUH0Ul5OLIJUlCSqdIXPcrhr79LaREHCkcpuV5zyV21Vu4t1h9/ILpYWa6Ti+165xcsV7juO2RZR4al4l6LvBCBl/csKDwJQ8yODGLLMn847pr+OL2t/PvhYN8TnX4B0gyBX+bS0ZeQciF0T9Llu3Lz+LLTaP7wuhqsKaIT4utQIn/2h3u1WrH5sgrRSuJwjTffuIbnNV7HqcndqBbmbMcPJkvP0OhZTWqY/TiHSezavcPUApJb7EWpbH2Unbqss+O8V8Pls1454hhh5+OlY4BFk3QAkYwtkMRHR5CsqqaWFQy8gNVYoGXOPJhcuG+ymwygVjFMbvaeO6J7+BrO67hXX97EZ994DOMZQfQhcayyHL6W/vJFJPsSz3JTQNf42OPXMkH7/wgU2WZwaF0tGvUaqUerInqVFIBSMUcip6n6G89ZPlbrpjhN+EAT6eDoOL0vBTIWspWdgJG3I580lvBeoiiFUkvoBSzFXLdyjqqV6L5zG8uBPJD9xP6+DUsf9XzSXz+vyn09JF/67+y6Y7bWfG5z+FfVz0xt1JI4svPujhyF6J9oGuopnhCCHe/qsrx6ymJgdT4PgAme85BFkUy5t8RtYX3nXYtM7LgQ8oAaa00t/OBTpfVSjkhN04MopKQmyGOa/ETu3PtSMU8SrpyvtsceUU8c8EP9lxLVsvy+i3vNGLee8Qj9+Vn0crm8Oiq5yLrefq3/W/1TpmYTrlNUeViFl2u7QzmbCtjKle7hm8FjPR2pYLGj2KwvalIql7tLJRsZVFx5PGMt7xRCKcz0LIKzWEjys4dsw/wiwNfYTS7H78c5Pze53NK+8VE1Ta6gss5bU07HRE/d+waZzQ1yo7ZB7hh/+f51VO/4raB27lk2as5p/sFyGZaN698nbUg0RgBLATaXFYrPts9v3mOvBwPTv2ZlCxxmebmqmS9gCR0NyEPRBtSds5lc7G+VXls7Ip2atSdzUbI37KLnu+9Bt+QIddNn3shLa/8B4aWb6Ij6keO1Ob4oRRN0+mb4GrWdGDxZcYh2lZx3ylasTxuPd/F49qO7U9yGjDdczadY3eRGN6FtsLox9pAH/9vbIJ/65X56b7P8br1H0aSJPKBDpdoxadZY+k0PywRXDU/gz83RTFkZY+qTsmtzUxJDIKy2XWvRMjdc+feiT+yM/4g/3La1ayK9aMrIRStMvmFLz9LMdDmOsHEO04i2bGVtdu/Tia6hoPrXu7Zr0S2UHHNyNdZf+2VGKQW8oF2WqafoKgESbVWhivQzPDIh4KFCnO7qDjy6hC2M5DBkbthcORlhNz8mdYSfP+pT/K/O9/LZO4gl/W/ko+c/GP+ccO/0h89nq5gJSfW5u/mad2X8l9n/J7vPOc7+JUAN+z/Al/f+T6yRSNGyHFdxh5ZK4ytE/VEK1aHtUC7y2rFcs8/VEKuiyK3Dl/PJk3mpLL90nLZF37H8TYQcyk7BfMXQcJTrutRrsJpK5Mh9Jebybz6FRz8SY7sfeMIn5/iW9/Nwe/exPS/fQTppFOhrK+1NoRg2khgnXXYFLvKm+ZyFvEUZfctwhTPFvjr7kpOttoLCmGEIAaY6j3HaCqx33YukVMjXJzO8PLIuTw6fQf3TxpK0Xyww82Ra2UcuSkWdnHk+WmKDXDkGZOQq8mhinu+/Ay67HNZaaW1BDcNfI3+yPG8ZP0VgHFCdUbXBEMMohQzaIFWt1Jektj+zK8x3XUGWx74kD0e5ePktXaM8Bi1CXn5Z0+YYSwS7ccbVjpl5bRAG2RnEMX6CmzXqcxx3bKemm8sKo68fDeLBlWSWcOm1jY9DC+riM2g+1uMzDcO2Z6uC36w5z95wFwA5/VcxgtWvoF1nZ0cnMl4uOh7r/Yz+s7ga8+8kU/c8yEemrqVTz/+Jq7adC2qYipbGrRaqR9NwoAWcItWFJOQa/4Wz92+NewjV9ApClHT1fnR6b8ynhviHYUW1PKFZh2FXaKVmB18vxwNiVaqvK0klfJ1alXkoq6+JWZRD+wl9Lfb0H//SzowxtHfHyJynMrjb/oWslwSYTerbAqmKjlyF8osiTJ5DSVQIgIWQ15L8VxNbBTMjCCQSLRvpagECaVKFiOTwwfoB57d+RzukWb52b7/pjPQx+ZAp50RCsDnMZaG1UobYDgP+XLTaF0GB1prM7Y4cjUxBGWicDU/a9TpIKq/GvhfMsUUL1r9FhRZAXSKarhCtKI6zB/LdQr5yDKePPNTnPf7Z9M5fDtDGyrj83n1WdZyDVusWNi/+Y20Tj7C0LpXuK7bXLtpFSRlpqlHOo0NpvL6RDKHEKI+49YkFhUhL4flECKAYMoIsZmJrKTc4U83OUkrjVmqkOKTd19rE/F/2vBRTmo/3/1QOYdUox8SMq9e/37O77mcr+96H1/c/m62rLqM02hctCJLEloVG3aAVF5jYDpNixUjw5wpTo5c8uikLEm0hnyMxmtzAg9N/ZlWXyfnZqIo2rTrns2R+0rclvBHkHIJT6uMQzk8SpJTtFKFkBc1/Nseo+Un1xF48hHXw1P/8iGWnXsWkbvfT3TWiATo5PKs7joXUq3jbjA9hC77Si7axgMlmLb9qnlK2jeRdi3ghrPAlzMOGGaBBX8rQvaRDS+zxTyAHWSqEOrlnzZ8hC9uezff3PUhTvadz6pCHKmYRyj+CjGVIcMvOc0FMiP4c9N2GsBaBKboj6H5oqiJoQqDAl9+1iVW2Zt4gvsmbuYZvS9hTXSLfVIxRCtuQm6JZbRAJSEHw8szG+qhfeI+T0LuBbmYa0y04vh9qu88bnvR/RUKUuvbzMoxVoEpXvE2cvCqVwhQFQmtKMgVdCZTebqizW0y9bCoCbkzA0woPYQuqeSCPZTr3p3maDumdvDG/3sj07lpzuu5jBevfhuKVPpwcw3yBLA2tpV3HPc/fGvXh3n7/p/x1YC/YdEKwP4qUf7AiOMN4FdbkUURtZBA87c4OPIYlU74xhjV2/zHsoM8Pv03Luh9CaS3oaTdyQnso7A/as/QwbTKilzS9mJzK/icx8rmxlOSJE/RCtNTBB97mNDf/kLw0QeQkyWLmdSFl5L4h39mxeZ+MpNp6I6YaesqZfjWxtMoPxRMD5MN9Roml1Z7eY1W08/A5sgdTjiNmB/WgxCGB68VT7+CkJtJxnOhXmK+KK/b8BE+v+2dfLrwID/BeDYf6sVnuvS7wh0AQvGTC/YQTuxH1VIUTK/QauNivUcmvMJbtJKbsQl5Spvlu099nHZ/Dxcte4X9PmD4ICjFrHFEMsfUTcjLBwKQJGa6zqBt/D7PieNF/BsSrXhNwhpWLjmf8S0ys+Pgq03IdSFQqozm/AexXcSEXC7mkEzSpQsjLnM2vAxkpYJwFf0tCODPY7/kph2/wyeFeM2G93Fq+7Mq6rW+beV8amxBLg+v423HfYavP3EVb+zr4S3Z/awKdtZ9Tmsw3qtlp+vLz6D5W+zsQEV/i+fHbGTS/N/Q91FlPxctewXFkU/bHDjA6s4w+qzFkYfBNA4o+iIVsVbmQ40jYWy4WlZGvvcxYvGHUA/sIfjI/cgpo738puNJPudy8ltOoNi3wrMezRe1sww5YS16N9dcvT/B9MFS5D8TTx6M0xMLoCoy+ILgC9sceTlsHUI9K5WKE6DAn3MS8uV0Dd9u3w9kxtDUiG2Wtzy8jiv738P39nySa7o7eVlmgr6Va1ELcQq+mItAWQQsG15G24yR4UgrS8xdDq1YeqYlUZlgwpefsaMF/m7wOyQK0/zL8V8m5nNHVCw6kktYIRicnqXVRF8z3WfSN/A7U7xU8tyu5r8sFbNoDYhWmrGwsr5FMT1RIVqqVW95D+dbrAKLjJBbg7Ppwf9g5Z6fsPPSHzMTNTLphFJDdrSy8mEq+MK8v7uT30z8gqf1nc/z+95VMcHK22hOlOou3BHo47s9z+Y1Az/jK099lPds/V9PpakT5SnRqsEKNerLTZOJrrY58mIV12hJkmraMqe1BA9P38653c8n5mtHK5NhbuqNsWtnJSHX1KixmZpH+KlUnkkzDvVcCLo8PQmKQvgvf0S97/fseryPKF833i3WQv6k00if9XTyW06k2NPnWYeLE1ajJb2Ig4g1G+silD7IVPfTKq5nCkViismlhzurhkme65JN54rGNzaJYy7Uhz87bo93IDPqFvcAp3ZeSHzqXn7JzfjHfsYnzn0Gk7l4hdu8NdUykeW0DhiiqVzMsKOuR2Oy4eV0DD7sKVpJtJ/Ao1N38Lfx33Bez2WsjDgsP4RB0CxzXKXoIOS56hy59edM12kAtE4+DJzpLmMWWtYWZHjGECFKWhY9OHerJ68+WGtPyUxXJeSylmXFnusR3S+B7vX2dafebokjB3zZKVbv/j4Ayx/5PIPnfYtkViOYHmSy7wKjkGOkdsUf4quDX2VnNMLzo+dw0coPIUvVj0/WQq/Y55ukTn1IfGdshpevauerO/6dtx/3OdoDPVXLNxJzGRzx1U0uRs3H0WU/uuzNfUi4NwlZdsuM7xr/LUVR4Myu5wCmVUGZDNOKG4M/AiaTaznrqFqKguKvaitbbbHIOx4l8sBj+Hc8gXpwEP+eUmYbuS2M8OvEn3slqWdeirZmHYoiucKBemEimbObtzhVRUtT9Lu9GsEd/71arWGfhD8z5grCZqHg7Euo3ZMjjwbVhqeNc5y6YgEmEjl8+WniZlLpbHg5EoJAZpRsdJVJyCv79Zzel6EP/IKbuJO1j3+LlxVmjeTMzrYQdp0WstF+EPXj4mfDy1FyMyiFFLpSMt/05WfZ5VP4zlMfY3loHZetusqzTZuQO+aYaotWWis2WuvvVMsGirKf6My2ij5ZXHx3LMDwTJbNfTH8IkemrmgF4plK08Vq9VscuZypbku+atd32fjYZ8nF74HX3GC348QCMOSLj5C3TBtHwdmOE4kN321bOAQzY7Z3ljUZ7xz7FTfu/yIAH5mYZHPv0xitQcSheqLXWgvS+Ug4oJDOFZG0DKslH1dt/k++sv3f+PG+/+KqTf/pksc70ShHLkzlWisJJjGsVjR/DB1vblOScGU7cnIG+WKW20Zu4LjWM1llck9FNWLYjZucH5TSrglHNEIrAqKipewJbvfR9qAowvAQ0T/8Af+eneihMIHtjyOnEkjxWdqs8ooxJqlnPY/8i1/BKblfEXzga/z5ZW9taEwsOO20LUWpqiXdhFzLccafriCABlfdbCfL8EKXmkEWxVIIUwdcMX/CHfhmKzlypyVSM3yABCDKRSuG1UwwM2IT8pnuMyue1YJdfGJiij1t6/niQ18koEd5YblpqtkZZ27bXLgPUlp9jtyMSx5ID1OIGQm5jZNZhi/pu5CQuKL/X/F7MBZClJKXOAm5Lz+DLqkIX7RiDluyeSH7SLVuIjb9pOu+pbwFUGWZi4/vNV8xT1d7bflHeXCretCVAJoarkjs4oRxYgAx7A4tsBBcuBOLj5BPPYZAYu/xb+eUO6+iY/Sv9kJLth0HgCby/Gbg2/xp5CcA/Nepn+LSn7+Kbfn6H862MihnyBtciSeuaEWSJOTdRtjX1ZHNPHfF6/jlwFf55YGv8NI17/B8rtEjf948Joc1Q6Si5OO2DkDz4FilMsNGxRGc6a7x35LUZrh42Svt+84EzJpJyG2O3BcGDAJW9EUpFiR8g3thdj9yIo4yOYYyMYak66h7d+PfuwtwJ1sTPh8s7yIgTdC6Ns22532R7KnngEnMgz4F6fFEhcVKs2aDNiEvJHHGygyOP0KbudjY/1fYeEnVulXTjC/voeNw6TRCHaim16UTiix5xgaqhKj4S9FSyHqBvGnyZhPy9EEQOoHMmKGELe+XrwUhqXxEOYmPth3kcxNP0Ku1uxa61Vq83eD2M5GV6A5lbi1Y/QilR0iahFzNz/LZjja265Nc2f9vrIkeV/mGDmUnlBIwg+me728FSaoI7ub8NIm2LXQP3Uy+LBaMVcRJLCUth+KfJ69qx+8Ff3spt6cHIom9AAQzo5CahEhnhSliQ9nAmsTiI+TTj5NqWcdk3/lo/ha6D95GvH0rYHzo6dwo39v2cfYld7Cl9Sz+ecNH2dplfNBGYpJb665RsrF7LOmyNlFkyYg6WMiAaa73jL6XMp0f4y+jN7IivN7OugNw8qo2do0myFUTrVgsh2wstILagi4klNw0CIGSmTEIuRBGCF5NA1W1n5VSCTCTSkjpFL5ChmK2gKZleXzHj3iafxNbDsoInkCZnqSYiTOxLUo4/Buk4QkOJCfJjuxiINVB/sE30Dc2jpyIQzDAztQyWm78UGWXVZX85q3k121C6esj3buCzHkXUVi7ASSJdY9/gXVPPgCAuqHLJuIAsgRyLlHh1dkMHRfCkJFDZWaoyPjDpT/23QkbL6laj5IxHHgsYuqEyxqlioy8gVSTdn/LYTn1VHDk6WF8uWlkoXmKVpAk8sEOYvlZvvysL/OPP7qAD/qGeH38ITa2GM5QFtMw23kK6ed8liflzUgeSmAvZENGP3zJIdsC7y+jP+f61haeHT6Fs0wRned74hCtFJ0c+SyFQBuS5GZoFEVyJQNJtB/Pir0/47Fd24HSKUl49b2QhXkKj+H8hoVAO3K2ij5E1wilBki0HUdsZjtM7oKIh6HDkmgFojM7mOk6DSH7SK58Jp1Dt6EWEqSCPdyeeohfD34TXWi8cu17Ob3zWciSjK4G0SUFtYEIiLaMvGx1VbMq2TdhCo0LedAFhYEBMukkxSdH0KcUQnf+GWVsmFclob+wnMJfPktk6hf4Vm4icstvmQBae/rQVRUhSaijI0haAa2rB4RAnRxHD0dA15GzGbTeZeyc7EP89Gf0ddzKwakJJL+E/2cvRR4dZgUgVB/C77efCQMhRUEqujmZTwEwCbzNdX2cFmKPfBOwReIkCQI7oa2DQv8G5KiP8J6HSZ/5dOInXkSxp5diZzfFjm5joirG1OqKBYiXxY8PJUuORMH0sH2SAkAyEiHUIuR+Va6ZNg9KHHm5CWJo8nGyoV5EuJvQ6OM161AzBkde8CDkLlFYqN2w53eY1IEhxipZrdSn6P1dYZa1htg1liyFljUVbLoaIu9vI5geJmR6m+aqOCkVTKegzlAn141O8oaVa/nWrg/zho2fYEPLyRQ0S1EgUTj5tWSGZgh66A68kA/1ICTZNoW8ffQX/GL8Bi5OpfmHVZczU2UnsN6+GkduKWSdhFyVJXKOU6bFsKkjD0PfxXa9Xr4BaBnwBRsKQ10PTn1ZPtBOuAoh92fGkPUC091nGYQ8PsTQTIaDMxkCPue8OLT+eGFxEXItSzB9kHTsJQDMbriMtj03kRv9M29ZuYGd+z5Lb3AN7z7xkwQdJkpIUimIUh1YHIAoaqhDQ8izM8jpJJmgxMTUKK233Iau6/j2PYWcTqEHg+jRVtQJw673YFl9Hb//uP378+3fdsF2Q+yALCMlZlEzGbS+5eixFpTpSfS2DoSiIGfSSJk02rKV6EWNwpr1RNURRCDGTMdGepdPUfB1kAisRxcCPdaKtnwlQlWRikVC7a1oio+spiOCIXx9vaSKOjcNfZMOLcwzVrwCva0DZWIcqagR0ofZOPgNnrjof0gsO5ULz9zA2C+upmPXjRSu3sffnjKIW2R2Jyf+8QU8es4FZFZV52q9EE7uI9u+ieD0TkNU4ICEhJSvJOROtIV9jMVzVe+DW7TiRDC+n1RsHWq0i9CkQcirrXOLkOcDlVxVsYyQSwjTtr8kl1UkqSkiosoykYCKEMJ2+XduIrnwMgLpYcJpM3aMGeipHHbgrEKWvnyGDwcv4D3afXxlx9W8ev37keVnuspLSA2nwROySi7USzAzzIOTf+YXB77MOt8yrh2/l0dOrm7CaI2X5Vchuwj5rOGRjTs6oCrL5Cht2Im24ykqQSIH77IJOZQIrU0g9SIU8zCfuQBMFPxtKFXyllrhgBPtZn7c+EG2mf4fzg3yqLVakSTpUuDzgAJ8Uwjxn/NRbzkCiQNICDtW8FDv6fx22Ra+7U+SlzKc3f08XrTqLfRFWphOlbTRQlh2xWY0vGQcNI3IH28ifMctZM44DyU+gzq4H3lmiraZqQruVQPGocLpRo+2gqJQWL6K7OnnsGrdcoLtrSgPfQU5oPLgKR+AYhERCqNHouwIjHDdY//OstbjeMuWz3LGhl52jyUb0p5bWPenK/AFY+w571pO+uWZzKy/iIHTPmzHdnYi1BYil9eIp43628I+/rL/Jv7QkuRNG99Puu0sV3l17F7Ct30ZuSuK3t6J5PejFrOG6aFjBtpWKx622k5MJCoJbjA9Qmr1hfhn91YScgmkXBLNX8UlnkZkjKJq0u1A4gAzKy4h0NoLe38LWp5qS0vJWhx5pamqSx/hyNzkJOSyMy1czd4acHJqfpMjd7adDS8nlDxAxOHF7IV8oINQch9kZwCIBft49+ov8Z3dH+F7T32C8ewglyz7R/PEIJClkpK/EW4xE+7jZ4WdfH/PA6wKb+IjwfMIintqporTdB0hhEsHY0HNz1IwT2XOk0u5h7ZQ/Mx0nUHH2N1lN4wfdnEz8ijqPHHkZaIVtQpHbm2+6egaQykaP4it0XfgqLQjlyRJAb4MXAIMAvdJknSTEOLJ2k82j8DsHgDGQl386sDXuG30ZxCEjbFTeNGKN7NCXoY8OYNy522EcxrqyEECj9xPrrebkQMSubEHWOZ7PnLGbV4X+/X1CFVFqCrF3uVktpyI3tVDMRxFTiXInnEuHcu66F+3nIfGcuSEBH6/NQCuuto3dhH0KTD7VYh2oW890bXoN7CMFxz3r/x032f529TvOUO81vNda03AYqCNYG4C9CJKPo4eavck4h7doyg0/jT8Y1aFN3Fca6XVg5d5mKylK8LxVhNd1IXQ8eWmKIR7yIb6XN6KYCaByMfRIpuqVlE3vpgoeTI6T2FKIYkvN0Umuhq5bbkhCpkdQIS9OVs1M4kWaEfIlcvEpZwOmbb9+RkyjnQ2zm9Yi5hUOgOVZORO+XyqZR2dI7fTktxLLtCJrnpb3NiBs8zk05q/hVZ/J2/e/Gk+v+2d/H7oOobST/HyNe8G2l3mNfXGNl/M8h+RIjfLCdZEtnDVpmtp2fl9s6/VHd+sNWDNI1Urk5HbohXj2rruCJNloWkBpnrPZeOjnyaUPGCHnrWGz87hWnCYy84DXITc346SjyPpBVdgLXBuvh3kQr2o8UoPWDh6OfKzgN1CiD0AkiT9BLgcmHdCPrv7bm54spftB97HCbNw2WyIUOdaWiYn8A2+0VXWafyk791FsVPF1+Vj5oQL0SNREILcSacjx2fJnXQ6emubW+lWZm+dCfu4e6YA/gaPa6ay06fIaA7uXpLgrK7ncO/EH7jxwBfpbVM4o/3FFY/XIuSx9l6kA7vxFWaREHYsaS+U13Pr8K+YyB3kDRs/4ckZeHFMspZGV8OuCdgoR14OX34WWRTRgl3kg5VKQlmWkD1k5OXvVAu6KL2HkyO3oudlIisIRkzRW2IYTEJePlZKZgLNDO9qBWgrteFFyGfd7+KQkTcC50nDl5uiKPtdyTWSrZuQ9QId+35Ddnmlk5KFfKDTIJRx47RjEcmgEubft36D20Z+xm+GvsX22ft4OPVCntn1z6hSsKIPTggheGjqNm4a+F9m5Rlekkhz3un/jSyr+HNTaGoEvYYoI6fpSFJJtGKFfZCKeVQtheYvKTv9qsy67ihTqUrrkOE1L2T9Y59j5e4fseuUaxBCOGTk1gCYc9JX3bR0rrDzouZnK8xSfY7NNx/sIpIqRbtcqPC1FuaDkK8AnIn8BoGKWSZJ0puANwGsXu3NAdXDvj/cztZHFbbag5JGy06gt7WTPeVMRDBEbsuJhPtXE8/pFJevpNjaxpo1vbT+8rUE00Pse/bVDbVVTkRrxZL2hJl13mtZyJLMmzZ9ig8/fAVff+xLrDrjTIJ1gvC4no90IGenS+FsQ95eqkZbEj4zsW+yMMsv9n6bTS2nc3yrNyGwCIfiINCylqboC7tN2GSVohKwbcwbhRWZTwt2ovlbXZH6AMN6IpeomhHGeqdaEAiQFTOJQYmQW1nSc8Ee9IiR63LH7l2M5o3Y2l6EvFglBolLpGyOv1pOyOXGlJxeZfy5KUM+7njXZKuZQ7WYJ9RfnZAXLM7YVObmHXk9JUniwmVXcMHqp/GNbZ/mxl0/40/7b+WSFVeyLnIqHdoKLLKgC53p/Ch3jf2Wp5KPsS/5BBG1hbeHn84b9v6Qv+amyId68OcmPW3trTEI+VRmM3nawn6E4keXfXZEzfI45roQ9vf1+sz5UC/jKy5h+b4beeqEdwGxSievgsmE+MMV5reHioIjwXX5O/tzU+iSguaLGSep1D7POha1slMI8XUwfK7POOOMOY3sWR/7Ng/fczuB45+PUH2sXt7OkEegqVAsQM4lm5VMGXnjYgCv424jsD9SIW2bH5bfFwKCSoT3nfAdPv3EP/OVJz/CWzZ9xs4sbvS4xgQMdyAXUnYEPD1UnSOXJdi6vJXbd47zh6HryBYzvGj1W6rK6WxO1sGRK4U0OTVGqkyOX1QjdkTJRmGlISuEuij4WwnHn3LdV4opJIeM2wvOrpdzyuA4SZXFoQ9kxwCDsBUjhpetlBi1LWDKwx+r2SkybRsd90oo6oKDMxn8qkyXQ7TihM2Ra3k6bn4Xy6KnMrz2pXXfSwiDuys3e0y2bEJXQ4aicN0zIVNZh/F+5nPDpvt9mSs/wHEdW/nx83/MvSP38qE7P87P9n7JuPE4dAaW4ZeDzOTHyRSN8YuoLVzR/y88reu5dI/cgcIPCSUPkA/14MtO1hSrBHwyRb2Uy7SohGyrlbAwRF8GITciBJbLxp1QZImBja+id/D3rH3yy7DqY5U6hrwj7HJjDtMNw0nIy+GzHLgkmYK/ncLE/Z51HK125ENgRHc0sdK8Nu/o3HAifn2ZpwOAE+XXZ9J5Ag1arVRDszE60LKehNyJ9kAPnzz3M1x9x7v4/p5P8caNnywRjBrfWpiEO2w6H0jhWoRcwq/KTOT38Lfx33LJyhezLNRftbwzFoZdh5am4O/hiSF3jk7NF7U9a52QdA1/dtzTPM5vcsVaqIuCvx1/fgZFluiI+BlP5FDy9WOROzny1pCvgpBbx1jhj7o2GqvtfLALPdAGSsC+BpVDrmQmKC47x7MP06m8nV5sQ1eQfipFK4aHJrD3L0S2Xc9Wrmd8xSVoZZ6WXjPL5sidkBX01/4WeXY/9J8P28Y8+2bHTh+4F10J2DFCXFVJEj7Fx3krzuNTZ/6AbVNPcPvITUhKhgPxAwSVCGujPYTVFs7reSGrI5vt0BbpqKEHCKcOMNt9Bv7cFJnoKjqifqaSlXJtMLL4zJoKd2dyCb8jRVx5ujavRSBJRgCt4TWXsWbndxh++jsRZnKPEhPlEK3UNm5qGs447uXw56bssc4HOoz0f2Umqa5+ziPmg5DfB2yUJGktBgG/Enhl7UfmBqVsBKpxleWXZ9IFOn0xg+hUi/heDw3ScQmT5bY4csnjvqOys/rO4dmrXs4fBn7C3RO/45zu51MXJiG34m2LWB9U2aNkSaKoF/nBU58lqrbyz8e/mWFvpbtRl+yjKPtdIhNFy6B5KNY0Neq5Oa5//L/p3/4NHnzGd5jqPc91z+bIg50U/K2ohSQKGqoZqc4zhG0ZnEPqtb/aCkZ/OUc+bkYMjCCQINprc+lQZokidOTsjJGnkdpTZvdElpW+aCUhl4xNJbn7LqxtKTqznZmesyorKYMvN+XOG2nVueI0WHV65QMOWKEqmNlPIbbKs/OKg+1VZIX+2PGsCm+hvytS8o2ogmx4ObqkEEoYZnj+3CSp7lM5bXU7tzw5Wtln3PomZ3KJainiwHvMZfPUNLDh1SzbfxPBA7eT3vxisx3zgXxJtIL3vjJnWKekahy5db8QaEcSuuns1H7IljP1cMip3oQQGvB24I/ANuB6IcQTh1qvF1RF5twNpSNctSOY19FF88WQRdEz8WsjaJQjlyQMG1ah1+XILbxi/VVsajmN6/f9N7eNGIF2am01Fgcend2JQEKKVg/GJUsSP9j2A/Ynd/Di1W9lZWsnJ69qo7+ruiKoPIuLXEjZIhcnCv5WfHk3ly4V86zeeR0Aq3Z+r+IZX3YKIclogTZ78aqOOhSbkFfnyOslhbA+lR6IuZzAAtlxW665fzJNIdyF3xRPlUMpGCIeK4Veva0/73fnUjV7iq5DevAxI5QsEPMI+lR6r9JbGXFWKk9atcQOFoq+KDnzPQsRbzNOJyG3Tw40xuMIxU82vJxwcj8IHX9uCjnmPQe91qJzfpWyA1UScs/6zOri7ScYDlJ7by7JyJ1iTQDf/FitOGGJVqpy5IESRw648qcuJOYlZ6cQ4ndCiE1CiPVCiE/OR53VEPaXDhHVlF5el0sOInMTrzS1oxYcMroGoEgK/7Tho3QFlvOrgf9lb+KJmramwpTJRmZ3ooW6UFSvlBIGds1s4/MPfp6TO87jlI5ngmREiVPk6p++WJbpXNYyFeaHYCWCnnFdi87uQNYL5ILddIzdjaS7xR7+3IQxySXZEcmxxMmqOTMRQlnEPiecQ9MSrHz3kmglRtQhSPZnxsk5FH9JXxeBjDchtxIy6IGWija9oPnbPDlygEh8F9M955APtBOJ76rsb5lnoqTlULWUbSHhrrOx02QuZFjlZLpO9rxfXo0dp72h2o2s8225QcLFOJLQ6eoxgmmdvqZ+n4tKCL9u2Hpbm7gnR+7Rrr3mZYWJ5RcSPPAXhBkVrmR+6FR2zi+E7EPzxTw2bTOZh4MjhxLn7iQfCyFaWdTJl5sj5N4OIvMNCUp2rF4cuUffJEkiqIR58+b/R1AJ84Xt72J3/NHKglZ5kyP352fQIr0VIicLicI0/3H/1XSGOnntpquRJMkem1pzyRXKVjciOXpx5JqVds6B1imj34Pr/wGlmLFT8FnwOxRjmp0AeMbm3JSCtbAbSya9qiNMX6vb7M2yKBGBmC1zN9oed1lwFELd+HPeCZGthB263Y/SiEWDlRLJgr+1YiwkAL1IKHmAVMs6sqFlRjClOrCDdZlEQVWaX/n7N78egq3MbHyJ533n2nFa6zS6UaSjawjG93NupyGGkdsMcU57xF/rMcCYX2qxJFoRSDVFaU44uxdvPxElO4WUNHwR7Fv5EkfuFUjuUFFwJLg+eVUbJ61sBb1oZEkyCbh1mvI7CH5p7c0/JV/UhLzanKsmWoG5c+RNwSbk4YY+mlWiM7CMN2w0DjRf2/4BBlI7vctHe21nBC2yzHPxCSH4yd7PEs/P8oULv0DM1+auo0a3XKIVk7uxwo86YYhWZl3HlXB8D5oacWR+3+t6xjBV60TgyHbkiOWtmkkGjlu7ik293ou7vOuxMsJqR7D0xyBX+t6B7Di5UImQawGr/5WmDdaG7+TIN/ZGOX1NO2s6a4yFs5+SRCA7jiyKZMPLyYV6CWQqFZS2i7nVthXjpYE8mtUwuvr5cM0B9vvWe96XavzVCDLR1Ybn6JAR/IyOtd4Fpcrai2rYdtFXc0amK6dC0BL7eL23k2lJtBtZ7/3jj5vlzRu2snOegmaVQQt02CITVTZe0PLpKIlWTI7cjJToNDFd4sjLIFcRGNbmyOOVN+cRkiSVCHkDsR4kJFd/18dO5P0nXodPDvCVHVfzxMzdlc8oPrSoIfvMtm+ukJsWRZGf7vssT87ezVtPfidbOrc0tVSLarhk7WFyN9Vk5LLQXIrRUGqATHS1HUYhktjnesbJkVuiFScna/3e0dnrCjTkRPkCL98sbRtvfwxycUI+mc5gEbWQdNn+av42JKG7bOZL/TAzLwVKJ4M1nRHaI37Pk2DBXylmUiTJ9lzNhvvIhXrsXJsulDGNPsu6JuBtmz0fKOfI7d8bfD7Zamb/2fZr42d7f8NtF5UQqmm1oubjFWIVSw/m1RenbD/ZehwCyUHIncpOacEIeSHQYXtxWk36yjxxLesVb458/rG4CXk1jrwmIa/PkfvVuQ+LW7QSbmj3LSdE3cGV/OsJ/0Obv5tv7vog1+3+eIXTSGrDCwFIrrnY9Xxez3Hd7o9yz8QfeNOJb+b1J77GaMPsiFeGnHIYMnLzHUwi5ykjL8tWBBBOHiAdXW3Y0fpbCSf3ASXxgKEQsgi5JVopEXIlFzd0C6q/Yb1E5RgbD+rBFkBw7kofp7QZdmhOGblm9b8wW16BPU+Eh7LTc375K7l7RZHsWDLZ0DJyoV78uSmkotuUotwO2m8mN85WiW44H3ARcsf1RrnFRJsZGGrvXyC2vCl3+OU9XcjmvFKtWOQm/KpMQK2e/MXZv6Ivgta2Dv/4E7hUPoW0MYfmyPr2tNTO9ZkPttuEWzL/K4+No6tBirLfpci31txCxFpZ5IS8CkfuKVrxjobnBZ8io+bjrNj9Y1qsJATNwFZ2ejgEeRT3eo3+1jW847j/pt3fwyPTt/Od3R/lYHqPXX7q7Gu44/m3kek93a50x+wDvO/By3h85i5euvodvPlkD8efBqwTXHk78xYh95KRlxFioRNMDRrmb5JEOtpvi1bOWd+JrGVQtRS5YKcZyCyGkGTXRqDm4xBsq9456nM0tozcNNOUMlNISYMTdsrItUDlRmLBVnaa7+gcR2+OvLWCu1ckiUBmBDAiF1qOOYEyS5nyDcufOogu++zyh7Lsq+m0nabNkouoN9ZaIdgJraaH9hpvW3ujPo+55gsjFRyiFV+JkFc7KbhrLCHftZXAxOPufudTFZmfumP1EzFbOGllW837BX+HocQUwhAdSQ6O3GGzr5WJ26qdMOcDiyuMbRmqysgPkSMP6glO+tMVRBJ7KCpB7njhnRVOHDX75FR2ztEhoSMSYG3XakLKD/n14Ne5deRnPDbzV87vuRw9fAEBsZxMqItEYhc7kyPcsPt3PDJtZFl/7foPc0rHBZ6cViMxH1wy8lx1c8CCQ1kJRmZ3Rc/bdszpWD/tY/cY91TFdr6xzeokmYKvBTU3Y/dPzs9CsDFTtGqw7cjDpmgiNQkJk6A6ZOTFsv47oZrZpPRADDKFuuTNlvfnZ+zUcrIMgcw4RSVouG2bm4g/O2EnCgenjNw8tSSHyIb6KHckmQvKvVXt69V+b2bXOO01cOsnYNOlzXXKH0EuZkEvouZnSbaWxqKeeWV5/wrdJxDZfRNqYRYwTSDN8BhO+JT6Y/m0dR0NlcsHO1D0PIqWQpKM06UzYJbdN39L6WQHnLa6nZl0wSUemi8sakJenSOvRFENIySlIULes/OnRBJ72Lf5jfTv+AYrd/+Qfce/pfGOOVOjNUDIvR0foCPiR5IkLlt1FWd2Ppsf7/0v/jp2E3eO/coso6ALIwaMXw6ypfUsvvzs/+SRfZXKuw3dUR4bmiUaqP/Ji2oYXzHNcctiMJ4wr1UenS05oCXXDqWMkDtWVLpUbB3L9v/KjqsRMJ2BDBtnYdbR5nKuUHKzEGqr2b/y8bI3AdmYE3ZyEDPgFekJMDnyXLAU00Yzia+ar9SbqIUEuhoysysV6sqRnWKmrOnorEgS/uyE8b6SZKeM82cnXMmkyznyQGK/i9AfCpyp/ZyoF6+mIZz/L7D2AljldnCqm/jDJLJKMYOamSTf43Dvd41zZR/Lr+Q7jfgz4fhewJTbF1IVop5GXjfmYcrqBad3p8RqJEpmhs6wCprPbckU9Cn0tdbOGTxXLGpC7vw2y9tCHJwxCKinElSSKPq8PRHL0bH758x0nsbuk68mOrudlU/9kH1b3tzQbHApO30NREqUqk1Y97Vl4bX869avkCzMkg8+zI6JAfbNDHNKzymcvuw44jN9qLKP7nA3YBAt52Jtj/i5YJMzeFLpZzkhKaphpEKala1BGGyAIzetTkJJi5CXOHLAcBxhLYFcyUXemniavxXF6RCUn4WY4QLebMQ4Q3FcSuYgRSyOfAISI+iSanPORv+NU5Z1/N3YG2XXqCF6Uwtx9ECr/RVciQGqiFacdYGhmHOaPOYcHLkTtm724D0gbSQ8tY3JDa9y3T9/Y5c7mcUhoqqysxn6rqiwujJ419PWdTAwlbE9RCvmd8CYS/7cFGohQaSjpAuot8GUj32h1VCqh5OOZA/5yjhH87JxWdU7vDut/lgRIK2E5QCaP1b61gvs2bmoCbnz46ztitiEXK1ydCk6jjpeUBUJNT5IaHo7Aye9F4DR1S+g6947iM08ScJMNVUXDoeg8r5EAipaUZDKlRxlnHPMyal5Iepr5eLNV7CvM8XusSQr2kOs6Qzzt/hkRdlac7eWwqWohpEQxskiV4MjD7QZVgMmpx1KDSAk2VbS2YTctFyxOHInIS8E2gjmSn1viCMvN54rexfb/DDs5sjzwS6XuEIrs5oJ+UrcklpIoPtL5o/1OfJKMY0sS6jZCXscLGudQHaiLMyAoOvgn+m4883Gc8Bs58mu+0Hf3Di5ajPJ/T61N6lmEVAV2sM+9lWrz/wukfhuADr7VtEZ9TOZzNcV85Rf0lpWIZBMZsGEpex0oJ4047hljdmxA+RMhWaHnCDiV8hruume3+4qV/BVBoVbKCxuZaf5dTb2Rl0fvdruWwxU2vo6EVAVOocNOfPE8guNn33PQCDRdfDPjXfMISM/YYVb3isBp6xqc19zdLevJei6Vv681zPVHaMaOUFUXrMVm/lUTRm5kH3kg522bXQoOUA2vNy2cbeDK1mE3HZ06bS55oK/1SUjV3IlZWfDVivl/bKu+6NGAt7UBCSGXfJxAKEGKSoBe044x0stJNADLZ7j0xqqPIJrHl6qhmil5E0qFD8Ff6unE1Lf/l87+i8x1XtutddtCtXC6M6ZC58DKuq3CPms6eUa6bHHvppJcbW6hBKkEF1ux30BTGVnuWiler0tIR8r2xuPXW4pNNcE0y6OvDykguZvsUNYLHQ88kVNyCUJLj6+lzWd7o9WzRNOC7S7nE/KEfDJtE0+QCHcY9tBF4IdJNq20DbxQOMdc4TRDPqUuuZMTo6ofN11VPGUWwjvMAsuQp5PIpBcVittYQchi/YRNANPhVIDrvRjuhomG+qzLVcCuUkK/lbX8bPgb7fjbaAXjVgrDcrIw/5KLlXCETteAiJdkJ6ExCj5YGU8EMORx2jfqYRS83Hb9LAcsiyxtts95wo+t5gGQBEF/PkZl+16LtjlKVppnXyI7KbL4K1388BltzQUe+Ts9Z2s7JibrXS1+VPrNDhvMJXQUYuQR0sb7Fz0gPmWNQSdHHkuDgH3t6tVb7MbWd5vWkOlDcZEwnDPL08JaASFS3g6nM03FjUhr4ZqLutFj9ggFtZ0hlEkidbJR0h1n+r6uvGOE2mZerxxFrGQMpyBZIPQlD9WTVnnda0hO/R5pulaOUfuj7hEEs6NU4v02k4u5YQcDPGKZUvuz4y67LjBmuxJKBZskz+bI6/Rx2du7ubsdWXhS8sgSxjcX2oCEgcr4nILDIWUxTU5542vkEAPGDGya7Vh16X40dRIaVPCCIMLuAh5PtBVYX7oy0wSSh+ksOx06NlCoiz1XLVxiAZUuqK1mQRnbKJqcL6bps8P0al5GjQ58mjc9FyO9JR0EXXMD702oFxLP6GEg5BnZyssn2r1p9nlU1Qj6LLPJuRgcOTl8eM1XwtWUu6FxqIm5NU+QLXjmRbsqErIe1qCqLlpwsn9ZHpOcd2Ld5yIrxAn5FSo1ELeLaOrRZAk6dDshCVpfhU54JCH51Oe3I3LKSPSRyA9ipJPEMhO2BYrFtKxtbZoJZg6aFtjWGNS8u6Ml6xH6nDkYETCtL6z60Tj6qcE0R4Y3wGZ6Yqs84Wi7oqRojhOcmohgQiW3rtCLu/Rp0JZ4CwlZRBs5+aVD3ZWcORBc6PTOzYAsKxt/rK/n7KqjZPLRHlQXUbuFUdmLqg5I+2gb7vNRksbbL2Z7OUWkW9ZY9hxZ2YMrikbryDkFklY3Rnmgk3drm/dtF5AkgyibXHkkplUoizuu5cCfKGwqAl5NVTlyIPtqPk4LYHK15YkCI8bAZ/S3ae47sXbTwCgZfqxxjpQSLtkdPXSfTm5b0ssZDFG1abYXNyqXc/XuFcSrSQglwS/Wz7ufFZrXYM/N0n7+L0AJNq2uMqmo/348zOQniKYHiITdpvVWVYkSna65GFZzyGoRuddMS0AOjfCrLEBp2PrXGVzBd0tWrGN7YWp7HQQ8rI2vSxICv5Wl+hOTlkZiZwceadhtuao0JLv6m39AGzujbkIzaGIV/2q3LAzTFvY1xAHf8hQVIrBdmShGf4ZaqBkQjoH56SsqUxmeq/BfIgiBKv7ffhV2VXzXOy6C4GOEkeeT6EUs5UcuTl/1Hz86I9HfiRR7RhWjUPVAm1ICE+XbAkITzyMQCJbFvoz1bqRohKgZapBQp5PNsyRl1p3c9fNHHHn2+XXDiGbnYVcAhFwa/Sd7WmdxwHQd+AmoJKQp1oMXQMHH0TJTuPvWotflR3KzjYAlOxMyebW5MgjHjJwL9R0DOstWRqlygg5mN53pkjH8oKUizlkveBSdpa34eU4UhE4y7Zdd0RcDHaiFpJIxZKDQShxwNBDtFpJoCXXHK6lKJuPL1+yw69NOmtEPq5apxDCW6FuxkkvhPs8n2sGuRZDqc7UXmPOQlNOZcf1NW6xYiEfaLc5cjljeHV6ycjBOwTEfGNRE/KqqDIZNDOanBWRrBzhsYdJtWxABNwcqJB9pFo22hl56iKfdrkI19qNJSBkEixdL5lOWhxfQ5YnjfWqYdixLzLThoy8nJA7y1qEfOD35IJd5Mvk0OmoSch3/AEArXW1iyxZhFzOTZcSMZuJkdvC/obiWzv/colWkGDZSfbf5fJ7o/2SaMVnUqrKELaVWN0R5sy1bg5MjZaJ7uywAE6OvDLhQDB9kFyoG0ktcc6NJG0+FDjHzIsb9nymCSpbj5vWWg1fg1zbBlf5ei76Xteylshsak9VQl7eHyu8bVcsMCfTToMjN0Rk1skrV5aM2VKAezmczTcWNSF3TUZqTwCAQtjw6vN7hBJFCELjDxPvOKnyHpBqWUekUZvQQrosO4l7UZb3z2m/3BE1LDo66yixnKhnslULXoy/TcjT08bpwl9pemihGFtJ3iTGlgjKiWxkBbqkwqPXG393nejZlupByOHQZLYGR34CnPkGeNm3beWzu/02VC2NVMwjyxIXH99LQDPdqoMtng5BYIx5uRmiGul0xyRPDJMPtLusdCxC7ncwE4HMqJ0IYj4Q8Mm0eJhIVkOJiNYrNzd4EfVCj7HOUj2nVX+uAWWnEIKiGjYYiKm9hk4HKvQ6kYDx7ctDHqfz7sQnjSLvL3HkUsrcsMuYmIpYRAuIxU3Im5xZVtorvxla1Al5Zi9qdsrliOFEKraeYGbEM9lwBcqC9tRjrpwLqCXo4+Lje22zw/nmti3UiikhFD/CHy1x5MFy0Yrjd1lidNXzEEgMrv+HyrpklXjHCZCbhVAHWmwVQgib47SOo0pmCn9uAl0Juk4AXs5dFVeqseRWZ5//WTjBO3u9Jcf0OcIb2/ExAtVl5F7Qg212fPbzN3ZBfJhcmcljiZCXLB4sQu5OYed4pVonOsczy9tCKIrE0zd2c1bZaaHWc7b4CKnmezalVK9TNHHyG3nyjE8xddwrXX1otA3rFLtnPIUuhCFeqcGRd0YDnL2+k2WtzZlrukxtHdCCHUi5OGh55FSlCA1KnsNLVitzRLWpYHPkqZGKe8qgoayb6fJObJtqMQL0lydK8G7IbbWydXkrEUeMk/n2pJsLLA6lKxbwVoaF2iEzZVgABGrLEHee+gHueOEdTJpOVOUYW2kGVdr0nAqKWFQjFJUgSmYcf3aSYqjTVcZrfCosF0xCV9SFS57cCFEoDzMATkLeWnJUaZCQW/HZVVkynZB6XWUKZrwVl2glM0I21Ouet3OQrKzrjnDhZu/cmSs7QlW5dPvUUfaO5XldmyHkdc01AzEOrnsZosxxxymHL+e+V3eGHZtOCdlCkVxLv0HIk+ZpO+IWcwANxRkqR7Uxy4es8A9jyMlRBFKFQ5CuhNAltSJz1ELgmCHkjcyx9pYW8oF2wtnK4P7K0H0U/S2kWtZ71pVqMRRlDYlX8m6rlZBfYcsyb3dv4+/qnV8oOh/2q5y0qpWty1tY1VHp1SaC7RA/aHDSkV6PGkr9E7Kv4ljpxMDGV8M//xGe9xmjPA46JUkUwj2o6TH8uUmK4e6q9VRDNa6pkbGzCPlx7SUZkyXT1AMtdgCuRoiYlUvVl581vmliuMJ23ZaRZyc5ZXUbSiGFWkiSC1cXrdTTsTSC4/paXB7FXjLycmzoqX4Sq4d69Ve3/fcWkV58fC+bemOla457uoB8Sz8kR2DKXJ/RBkVVdU/L3h3NR4wcpcwMIKXGjKxXctlGIUku786FxKIm5M16N65sD6G2r7HtwZft/Tln//EFhONPoe78HZkV51AtdGgmugZdUhsj5IVUw4mXJQ6NWB+Kh2dPLFhVxCJC7TBmZnyPuQl5IycK52Uh+2D12XawpHIUQj2oqTGC6RGKNTYNr/aBqokIGhkZy0KnXS4lm3Zy5LreOCHXgxYhnwFdg+SYiyOXZaM9XVLx56boigbok42oeblQr3vMFjrKkgtGw4ZPw+FlKhqxz/eyeXc+pwthEHKA/XcZp8lGAtY1gGonsWTQDPQ1O4CcHHXFuW+PlBgLw+FtiSOviepmZ9WJi9y3FXViG+HEXrbedw3R2Z2c+4fnIqXHmTnx9UY5j+kkZB+Z6GoiiUY58sZjN9RaH3MRvZy8qq1qjJaGEVsGpjZeilXnbqr1ruq3gQouSAt3o6RHCaUG0Fr7G67LC07utXzsFI/QDSULnRn7mmW1IgIxO0lFQxx5oETIpeQYIFyEXJIk4wQS6MBnBgrzmfqachFMo2iKS5a8f7ev1dn6mhOt1C5bfY06mIQG2xJCkLdMEAfuNubuPKHaO2fDFkd+ACk15rJYURzyIc1ncOQLLUFd1ITciUbGSZIk6N2Kkh7j5DuuQlMjPHHmf7J/0z+Rf86nya48r+bzqZb19TlyvQjFXJnVSvW+Gmu7+a9c67jdHQtUZJavBa/WRfua0v0aC6M6wS7dOKO/vbJ+R/8LoV4CM0+hFLNoLasrylaKoqp2x03Iy+6dsLy1IkuLnTAkO2NfU/NJw9LGF7bNQBuxoRam/buan4WESaCDzjgiRo/ywQ58WYMT95uKsmyor2bMnUYwV2LhtJWfL2WnVVJU6Vf5pZHZLACzmULV/pX31YKuQ8HJAETntil6odo762rYCDUwO4BUpgtxcvEFf8thkZEv6jC21VBzuq020lJFkvt4/KxPM9L/IgB6NnQhzRpRC6vN11TLeroO/hmpmHeZlLlgpkabL468GizZbQMJTeYE0VYi5MT6YKpY+tvjqFsBx+W2cGmsDBm5m0rlYqvs3zXTu9FdlUQj2r+ATyZXKMm6y80yu2MBumPd3PJkSUdiZY4q58g1XxQkyc7f2oh/lrBEK7lZmDHEd1lLlkppgecDHXboX1/KyOmZC/e5xmwhzMjdDETznO98cpVVTYSLes0y1bogAi2mgn4aujYeegdr9AGMJBG0rYbRJ5FTo3akT6OPpYc0f6s7xO4C4ZDIgCRJL5ck6QlJknRJks6Yr07NsS+NFVxxGvFLPsvuE/6FkTWXu+uo82iqZT2yKLqD2JfDEYu8kb7WM/mqBouQz4fVi2cVToJqxfVu9Fma25wyHSVv0ELnlor7FRy5Rx0XHtfDeetLx9tGRUtC9hl28i6OPI7mb0FCor8zzKqOcEMnHGGGFvDnJpFmjMVrpb0z+m1y5IFOfFmLkBu25rrithzauqLB1IJInr83A2sO1ZtKzdTuJdd212VcK9/U60vMa9zpNJyLWNN4COB6+2W1MVnTGYa2NWBau1nRUstR8LUsCoegx4GXALfPQ1+axlxpWO6kVxup28o05M5M82f0t9NpOudYxdLRfgA7mp8nbI688azic1mAzchu5wKx/FRoWQknXlEx0O5FCluWt1Q4WlTfYCoXcHLZ2aSWn8NE3wXo0Qbkmx5VK7Lk4sCbGpZgWxlHnrCVoKois7kvZnPmNeELkA31EEoNwcx+CLVTdMRxt7qXC/XgT4+CEPhTw2RDyyr6vKw1xEmrmtNzzFm0UuPe6s4SQzKfnp3NWcBUftfyvqiKBC/9Fjz7E3DcCxquu9rJx1JYVntnIYA1JVFsNUJuh4BY4FC2hyRaEUJsg/nhCg8V5XLnmmUb4CLbwn56WopMJvPIkkRRCDvTiytkZjmqcOSWHXl/V6SpvlaD5VCzYCOvBuEdD4Bax8NUghVtIVa0hVwii1p21xWLR5LZ/7wfc3A2x2aP54zNagGtOEJtbo68kPBMpNEIMpFVhFIDSHLE4NgcsNZJLtSHrOchPYUvNUzSzKhU/urzHXO+EQujcmzqjSEB+yfThyyDd12rUrbL4dPg3Z73k6osQfsaOPcdh9w3gFNWtVMo6kyl8p73V7aHYOMl8HsQkuyK4+Oss+BvQRI6SiHVVL+axWFTdkqS9CZJku6XJOn+8fHx+g80Uudcj5LVrps3Ko575nXN30re31aVIz9+eUtVjtynyFx8fG/d+NGNYqE5cknCMOHycsip8sw56zur2nSXo4KYy0rVVTWXV2zqkXKOPB+n4G9tul0JySDkyQMw+rgrYBeU3sNWjMWH8KUOkq3int9Q+/Pw+RtdR3MSrdTZf615cGZ/B+t7oqzrKq2bWnO7/NZc10E181tFlgj6lOr27rIEHWvhFT8k+8Y7PcM/gBHvHhbeTb8uIZck6RZJkh73+Hd5vWedEEJ8XQhxhhDijO7u5p0+5hNVORPr2C+svyvLG4kSvDny5W0hwxMSmoq+digy8mYi0lVt/xCogfPZSEC1w6BWq1OSzJDRhWLV+/X611B/m3klT468+Yh4kgTplnUEMyOQGoflp3p2KWs5/0zuRs3N1nQGaqr9Q7wP3uN/OCzaW8M+1nZFPJWwTsaq2lpR57gQvLJMOVF3rm15AXQdV/V2yU1/YeXkdUUrQoiLF7QHh4BqdrF15XMN1OeEU0yQia6hfeye6pVXCdpTq525EFKLA5nrBF5o1HNpF8LweM3ki3WdopoxP7TLNDOmHhy55muh8bBTVpsYuTYf+5xxYdVZMF1ZzrZB3nMbUIrKWM5VNkt45yrirPfYoVrQNCEhcRexTsg1TEotzHUZbK4TwtZrbNZ1N67/csYkX0gcnVRggeHJ9UnVCYZzgaWj/QQzI8ha1rtym5A3ztHNZf2t746wvidKb518oI11oMnikmSH/qwmhmqkqUabnQt5ampMnRy5lkcpZkr25U0i13MSk73nGQq3PnckTYvQ5oPdRrCwJ38JlOKkV86/w6t7ksz/ymF942bo+Xx0vfyEbNTrXXGtIHC1UO+58tb6WoOs63brT2ptIl65XBcCh2p++GJJkgaBc4DfSpL0x/npVoPtu/vi+L3ec15yX6mGyKWENRsNuWcoNehduZl1vlaGkvK+zmXOq4rsOoqesKKVE1ceojenA/X6VC28bDWrAu+yHkTDi1pUnGDqo2kZeSENWt6Onlfwxea0wZ25tgv9Vb+AK39YMRHtPyWJTPtxkJ1FINkWD3M5mTWzMc7dqqWSoDb6THm7p65uc5WrlwaxUcwly09DaEAWX34qdhaxQtke1Ry5EOIXQoiVQoiAEKJXCPGc+erYfMHTnrjKN68ucindkTuNKIihakb+2ThISsOxVsrrnyv6WoP0tsxfrsd6mCvRcI1lg6+9UApdG1aO0OyMTcgbyWJfDgnDUcQrmuSm3hgBhwljqsvg1vPRleiqEVq1EXv5hYbXUFvfqZlkF9U+meUc1sim1Yy1S7n563yhMg59ZZlam4gdb3+BZeSLWrTSGNfX+LVqsmvXd+owjsFVnYKsjDoNEh/vw+zhhfd4VO+V644ovydVlqmCRhNilNOPuXL7VWHlCM3Oujjy+fwy5Q5FE5teAd1bGNvymqrPNG01U+8kOmcZuvGcV57ShUQjIYzBiH5ZLXDaoffB/XezTEVRjaBLyoKLVo5JF30nvIa92qeoxnC4dtxQOwVfS3WO3CPrfD0cBWb484Z6ohXnVeewzvcYNFWdxZFnZmwdx1xk5M28Q751Lbztbsb2T0Gq4Pl8Y1zrwk8eq4liAxx5oVBgcHCQbDZLp2ZYJqXHJDrNTWDH9ikkydgUOos60ozEtmRlohcolQHYts2I364VdTp1gZSQ6LR8KRIS27ZVhqYG7D5Yz9e7XqsPANlxiW1TlfyvVR9Acar0vgDbn/NTimqIbdu21WzLiWAwyMqVK/H5GlO5L2pC3tAUbpLbdJczfpZ7DKajq6tz5Nl4Tfm4VxtH2qGq2da9Ti9eZerXM0cOcU5P1YBlKpqdsXUcmq+leY64Rs8q6/LW0xwK5jqe9axCLC60EcnK4OAgsViMNWvWkMwZxM2vyuQ1gxjGgiqSJFEo6mTyRVRZIlwl4UNB08mYZqpWgoecViRX0FFkyT4hKLLkStziRNwMwlWeIKLa9XJoRZ10vkSk/arsmeMzr+kUijpFXeBTZFfMmOishOQPI3WsrdmWBSEEk5OTDA4OsnZtY88satFKI/BWbHqVc/xeVkApu5CJriGUqkbIZ5qyIbdw8qo2TilTBB1JNEoSKpw06zzvHMq5yr4b2ySaqNASrWRmSjLyOdiRzwVui4yym4285zz0IWdyk0GfNzmw+JhGRCvZbJbOzs66m4ptPusRWthGDRn54RXy1IdflavKyiVZRRLefhOe5SWJzs5OstkqlnEeWNSEvLpCrfS719hWe86aHBUOKGXlM9HVBFNDSHplyE3SkzWDTFVDdywwb16fC4GTGrSIKYlW6pd1iVaOpKbApeycAQwlVdM9aoImNaLIO1wHNSso2Io273yWVoq0le2N5bts5GSgyBLRoIq/adm2pXltrLQiSwtn0dIghKwY4a2bQLOnq0VNyL3gmX+yDJ5cupdixetZSSIdXYMsigRTQ5UFmiTki0U83uOwiJEohXet5Kpr8+TOsW+UI7fsmBfM9ynUDkiQmoDMNLrsM5JAN4lm1l6Pxzydi0PQfCDoUzhhRSuqInuug85ogHPWd7KuO8p5GypzYdZEudxfav77e6MxSh4JqN5ilzk2Pac9YQ6EvOkmFrT2BUb5pHvm5m5OXNFal1hU5cgbNK+yEjFHZ3e6b+g6pKfmxJEfSTRqHeDExp4Yxy2LVd0457pGG80YXw+1OPxTVre5rUgUn/HNkqOQmjSS6ErV/QrcfXK2WR3OebihJ8rKdsM89VBFBPPNtSuyxLK2YEUyEIsYNtqeJEkEVJmIX53XDale+zMzM3zlK1+pW080oFb1haiGkF/Bp8hce+21bNiwgc2bN/PHP3q4zkiwb99eLrrgPE45YQtXvv6d5HOGmOT222/ntNNOQ1VVbrjhhqbar4VFTcjLoSpyhUlbwxOvzjPnbujkgk1GjJhk6yaEJBOb2WbEJbaQmwVRbI4jP0pNVuolhFZkySZG5fegQfNDs/DhHoOuaMD93cDIKpMcg/QE+eDcNuJaWdpluWRoWlWWWoN7PZzYurzVlQxkrgj4lEMSa9R6stoG2CghlyWp6ROBT5HZtm0bP/nJT3jiiSf4wx/+wFvf+laKxUpu+yMffD9vfcc7efjxbbS3t/GtH94AQrB69Wquu+46XvnKVzbVdj0saquVRuAtRvEuW4sbtIJBAehqiFRsHbHpJ5lxFkqbpkyLjSNvsnytxWkLVqqZHzr1F02yEU31s65VTVmBaI/BkSMMjpzGjtHWnDl7fX0FX63nvfp0NG7xzbziztEEiWzBtOgwXjQSaFwmXiwKsqa1i/VcyKd4MhAWrrnmGp566ilOOeUULrzwQh599FGmp6cpFAp84hOf4PLLL2ffvn284AUv4PHHHwfgM5/5DMlkko9+9KN1+/SrX/2KK6+8kkAgwNq1a9mwYQP33nsv55xzTkkRKwS3/+U2vnXd9/EpEq/7xyv56H98grf8e5H+/n4A5HmWEx6ThNw12TyVSrWJeyOTNd5xIp3DtzPjXIkpMzzvMSgjd6LcimeucHJEh5v5rGgu2guTT4EkkW87xSzTeKeqEf1lbUGyhbklFWjIOmdRzqCFw3/+53/y+OOP8/DDD6NpGul0mpaWFiYmJjj77LO57LLLaj7/X//1X/zwhz+0/xaA0AXnnv90/vcrX2JoaIizzz7bvr9y5UqGhty6ssnJSVpb21BVlaBPYeXq1QyNjINeAGVhSO4xScidaNSzE7DtQyP++sMS7ziJ5ft+gZoYAEwF0KwZf6V1xRx6unhQyyOzGa50wV3va6Ci6VgvJEdAVin0Ghux1ATTVO1dti6vbe3TjNv70YBmNo5NvYYJZ65QJGdy1vXstp0o6jqpnNuOvNyuG6qfEIUQvP/97+f2229HlmWGhoYYHfV2HLJw9dVXc/XVV3v2YS6QJAkk852LBfA1ZvnTLI59Qt7oNUmiOxbg9DXttEfqywenu4wUpeGDf4Nuc5ePmztzy8IQ8rXdETrmQXZZjvmkp7VtVtxoVH56wopW9o6nyGm67chRtx/1RCvlPWxdBcU8FPPkQt1mmSOLucYjWUgc3vaqn5yd+1+gShq+H/7wh4yPj/PAAw/g8/no7+8nm82iqiq6I5u20167nCMH0B0c+YoVKxgYGLDvDQ4OsmKFe713dnYyOzuDpmmAj8HhEVb0dUM+CYnhBaEPx5Sy04JbcVk5GWpxgo0QcYBU6ya0cA8dI38rXZwdhEBr056djWJFW6jh/h2tcL6v6kHIy8PiArQEfZy8qm1+N5zyutr77V+T7UYC6GZODHM9XcxX9L/DhcPZpWoGrBVXHAMVi8VIJAzv3NnZWXp6evD5fNx6663s378fgN7eXsbGxpicnCSXy/Gb3/zGfv7qq6/m4Ycftv898OCD3HnP/Xz6s/8NwGWXXcZPfvITcrkce/fuZdeuXZx11lkV/Xn6Bc/gl7+4EYDv/uDHXP7sZxrZw6xUkPOMY5KQO+HF9DUarKkmJAllw0Wo+24txSafHTzmxSr10AzxcX6HFW0hVEWatwiOTX/h1efYvybajFDFTb1LE2UbNVk8ZjDXl6zzXMivVChPOzs7Oe+88zjhhBN4+OGHuf/++znxxBP53ve+x3HHGZl8fD4fH/7whznrrLO45JJL7OuNYOvWrVxxxRUcf/zxXHrppXz5y19GUYw+XH7ZCxg+eBCAj33iU3z5C59nw4YNTE5N8fpXvhTyKe57+AlWrj+On/3sZ1x11VVs3bq1VnMNY1GKVpa3hTg4k6l6350u6tCXSjUxpjjlVciP/oQVe37CwKbXwcjjsOK0puo+Ws0P5wo7+mGV1/JK5QWGnfIzN/csYM/cqOCgA1F4+nsgn0YLtIEuDgtHbo1HeZzuoxWHc77WE4tWSwrxox/9qG7d73znO3nnO985p1584AMf4AMf+EDF9V/d9BtbF7Bl00bu+OtdpTgy4zuhkOLMU7YyOHCgOQVMA1iUhPz45S1GouMGsJDzTlp7PsU159O/7WtMLLsQZg/A2W9euPYW6F2OlOVDMzbn8w1PDvpZHwZAbK+tEPPCoX6buThlHQkc8S4d8Q40Br8q297PgJHIvJACJTDvRBz+LkQrC/vltQs/jD83xXm/v8S44DiizzcWg6lZfSWj8/eFe59GgzYdaplG26tWttZjf+/KzkV3Wq3VXZ+Z5zPQeL7PZnBMEnK3snNh2xIrzmDnqeYx64SXVWROXwxYiDE62jedRvQk8z0utQJleVtSzW/784EjTVyPwiFpDOEOQ6HesnJBql+UopVmUM15ZVNvDJ8q8cTQ3FMwWZN6YONrGFx3Jc86cdWc62qsvQWtvi5WdoSIBWvbAdfr41yT5DaLQxkqSydyOOzcD1kks3hJ26LAvI2uJJnB2RYGxzwhr7ZQVpuxNg6FkDvrF8rCmwUe6SV7XN/cMss74cytOFcidriIVyOthPwKmfyhRLYzWmkg57T3000Oxfkbu9APsxPSfH6tI30iOFpxTBJyt7v9wn74w8kRLYZJPB+y6YV4thksazWsohpp7qy1Ha5sMI3Aq1ovD8+Goi821TKe2W2WsPhxTMrIl3D04lBsqLcub2Ftd4TWcH0370Oh+VuWxXjm5u6GCKlPkV0B1ZpFbWVnI88f/Zv7QqBaZqFGox8eCqqFsXX26Etf+hIbNmxAkiQmJiYWtD9wjBLyeZvcDVRjm9AdhvW0GJZsvT7Kh0DJgz6F9d3RpvtUsz8eK0CSJNTDLMufq7BjMcyJ+UYsqBKqcrJYaEL+5JNPNhTG9rzzzuOWW25hzZo1C9YXJw5JtCJJ0n8BLwTywFPw/9u72+CoyiuA4/+TbMwSAoIJkDQvBEpGYIDE6QqOUgVDOkApjNOOWtpOLDqIU6a+dQRk2jrD1CHq2JaxjsPAB7TMtIwt0jLVglo/+EEoWmylSStToPKmNEiptoCR0w+7G5a4m+xm7937suf3BXb3Zvc8d3fPPve5zz0P31bVMw7EZdJwbR65g8872HOlW96t0GO2SbObqz0p3JVu6C/dLgh8Z/vF1XDyL0RUqeg7g5zbD2RFsiZK8u9qpiML1mfc3vEytv3eg4HK2Ka65prCzl7Ld4x8N7BGVXtFpBNYA6zKPyxnlZeVuLYeZiG/a0E4jJa+k3fpk/PlS32R2Na9OAbih/HiSz3yIY6R+/8jUVBOl7GFy4tmZVPG1gt5JXJV3ZVy8w3ga/mF444vNo/xOoRAGZXFGPRgsruYJXNvtOgMcR/4+sc90XPu7b3IuU8uL0ebrf8mql3m+nfgTBnbi6p8dK4359cuNCdnrSwDfpnpQRFZDiwHaGxsdPBl89NYVTG0BVUTfP1FylJq0r3u81VEM5QFzek5szm/kPjXjZrcQXlbJo4ZzsfnewNX1TLtgsY+41UZWy8M+m6IyMtATZqH1qrqjsQ2a4FeYGua7QBQ1Y3ARoBYLOabPliy+L2JG2jdyWykDhGMiEb43KjMhfTFxaGVoBgRLeP6QVamb6p257LuoZrdXJ22BLEf5FrGtrKykp07dzJ//nxg8B754sWLWbp0KQ888ADHjx9PW8bWC4N+a1V13kCPi8gdwCKgTYO23IkD+i6x9ufnOitOxp78BIjArIkDL3lnQyuDmzd1nNchfIYfzi1kklrG9tprr6W7u5vp06cTi8XSlrGtq6sbsIxt/69GahnbSCRyWRnbxV9ZxE+feobGxno2bNjAY489xsmTJ5kxYwYLFy5k06ZNbjU771kr84GHgJtU1Z2K6SE27Ar/fiGGKpecnPwB8WrWigknZ8vYflamMra/+e3OvvVZ83n+och3oOspoBzYnehdvaGq7tVx9bFhZbntypsnF672diElD8oGOtnZ0jCK8rKSvql/xZbHrT6KcVq+s1YmORVIUJWUCC0Noxg5LLdd6cgqRT6UOrSSyZgR8amgH5+Pjz1mmqoYNskjjyAPwxWaX8fi/cb/p54DIJmYgsqNr0o2ySq5+LIbY65+TJZ9P3LehhEoFQGYHeMHtpeMK7K9IGdGw5WMrgjW1Luh6jvuKLJMXmTN9YQlcuPoXPhshlZSjR3hzGLL/flxHDqb8wehFuBm+32ExxL5AJKLUpT6/V30kWIZ7x6K5J7x47CPSU9EGHZFacYFavwilNUPnVJVWU7zuEqurrGLhrJl48CZFe2+KWCD3ah+WFZa0jc5oaenh7lz51JZWcnKlSsz/s3p06dpb2+nubmZ9vZ2PvzwQ0dj6s8S+SDGVw0v2PJkXnHye5ZcObw84u0ceX92oOKZ3IuKi35QiFa7XcY2Go2ybt06nnjiiQG3W79+PW1tbbz77ru0tbWxfn3mio1OCPXQShgvuPG72iujlJYIYwM+k8cNuZ4/CLrOvZ10n+5G9dLUy3yHKSdfNZlVMzMXWHW8jG0/w4cPZ/bs2Rw8eHDA7Xbs2MFrr70GQEdHB3PmzKGzszPrduYqtIn8+klVoe9JO8XZeuTCuJHunMAMur4xchf6pjmW+Q4tN8rYAtx4441s2LBhkFdPXuCmvP/++9TW1gJQU1MzaNXFfIU2keez/JbfNFZV8O9EOU83hKGCY39+bJFeyuSOmt1c7csT8sme8yefXuR/Fz5FJF4krFCcKGObi+TFS/0vYhIR179j4cl2IWYVGsPBrSs7/VzEyktOlbGF7HrkJSXCiGgkflQ6bhwnTpygtraWEydOMHasuyU5LJGbUPLjUUZleYSPzvVSZuMgrnG6jG2ukp+7xYsXs2XLFlavXs2WLVtYsmRJni0bmCVyYwpkcs0Iaq+M2kl4FzldxjadpqYmzp49y4ULF3jhhRfYtWsXU6dO5a677mLFihXEYjFWr17NrbfeyubNmxk/fjzbtm1zo7l9LJGbUPJffxwipSVUubR2bBAU6j1xu4zt4cOH096fWm+8qqqKV155ZUjPPxR2jGeMMQFnidwYYwLOErkxxgScjZEPQcUVpYFYRbyY+XDSStGSNP8zzrJsNASDrXpujEnD8rhrbGjFhJIf55Eb4xZL5MaY0HC7+uHevXtpbW2ltbWVlpYWtm/fnna7Q4cOMWvWLCZNmsRtt93GhQsXXIsJLJGbkJk58SomjhnudRgmjTCUsZ02bRr79u1j//79vPTSS9x999309vZ+ZrtVq1Zx//33c/DgQUaPHs3mzZtdiwlsjNyEzMhoGSMLWJjJZHby0Uc539WNAp9eVETIe6Wd8imTqXn44YyPu13GtqKiou//586dSzuEp6q8+uqrfRcmdXR08Mgjj3DPPffk2NrsWSI3xriqkIv/FaKM7Z49e1i2bBlHjhzhueeeIxK5PI329PQwatSovvvr6+s5duyYQy1ML69ELiLrgCXAReAD4A5VPe5EYMaYYEv2nJNlbCOlUtDy0m6VsZ01axYHDhygq6uLjo4OFixYQDTqbQ3+fPfq46r6fQAR+S7wA2BF3lEZY8LDo7VK3S5jO2XKFCorK3nnnXeIxWJ991dVVXHmzBl6e3uJRCIcPXqUuro6F1p4SV6JXFXPptwcTmGPooyDZk28yqbsGVdcSgruf77cLmN76NAhGhoaiEQiHDlyhO7ubpqami7bRkSYO3cuzz//PLfffntBytjmPWtFRH4kIu8B3yDeIzcBNCJaRqVdrWpc4c6CGumklrHdv38/+/btY/r06Tz77LNpy9i2t7fnVMb29ddfp6WlhdbWVm655RaefvppqqvjFwguXLiQ48fjI8udnZ08+eSTTJo0iZ6eHu68807nG5tCVAfuRIvIy0BNmofWquqOlO3WAFFV/WGG51kOLAdobGz8QvLX0RgTHl1dXUyZMuWy+1SV870XKY+U2FFfDtLtSxF5U1Vj/bcdtAumqvOyfN2twO+AtIlcVTcCGwFisZgNwRhTJETElqNzWV5DKyLSnHJzCdCdXzjGGGNyle+g6HoRuZr49MMj2IwVY4qeqtoQSp4GG/LuL99ZK1/N5++NMeESjUbp6emhqqrKkvkQqSo9PT05zU23aQrGGMfU19dz9OhRTp065XUogRaNRqmvr896e0vkxhjHlJWVMWHCBK/DKDpW/dAYYwLOErkxxgScJXJjjAm4Qa/sdOVFRU4Rn64YNNXAv7wOooCKrb1gbS4WQW3zeFUd0/9OTxJ5UInIvnSXx4ZVsbUXrM3FImxttqEVY4wJOEvkxhgTcJbIc7PR6wAKrNjaC9bmYhGqNtsYuTHGBJz1yI0xJuAskRtjTMBZIh8CEXlQRFREqr2OxW0i8riIdIvIn0Vku4iM8jomt4jIfBH5m4gcFJHVXsfjNhFpEJE/iMhfReSAiNzrdUyFICKlIvInEdnpdSxOsUSeIxFpAL4E/NPrWApkNzBNVWcAfwfWeByPK0SkFPgZsACYCnxdRKZ6G5XreoEHVXUqcB3wnSJoM8C9QJfXQTjJEnnufgw8ROri4CGmqrtUtTdx8w0g+9qawTITOKiq/1DVC8AviK96FVqqekJV30r8/z/Ek1udt1G5S0TqgS8Dm7yOxUmWyHMgIkuAY6r6ttexeGQZ8KLXQbikDngv5fZRQp7UUolIE3ANsMfjUNz2E+IdsYsex+Eoq0fej4i8DNSkeWgt8DDxYZVQGajNqrojsc1a4ofiWwsZm3GfiFQCvwLuU9WzXsfjFhFZBHygqm+KyByPw3GUJfJ+VHVeuvtFZDowAXg7sYRVPfCWiMxU1ZMFDNFxmdqcJCJ3AIuANg3vhQfHgIaU2/WJ+0JNRMqIJ/Gtqvprr+Nx2Q3AYhFZCESBkSLyc1X9psdx5c0uCBoiETkMxFQ1iBXUsiYi84EngZtUNbTrd4lIhPjJ3DbiCfyPwFJVPeBpYC6SeI9kC3BaVe/zOJyCSvTIv6eqizwOxRE2Rm4G8xQwAtgtIvtF5BmvA3JD4oTuSuD3xE/6bQtzEk+4AfgWcHPivd2f6K2agLEeuTHGBJz1yI0xJuAskRtjTMBZIjfGmICzRG6MMQFnidwYYwLOErkxxgScJXJjjAm4/wNmkBlJjrmY4QAAAABJRU5ErkJggg==",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"import matplotlib.pyplot as plt\n",
"\n",
"plt.plot(t,samplegp(0.01),alpha=0.3, label=\"tau=0.01\")\n",
"plt.plot(t,samplegp(0.1), label=\"tau=0.1\")\n",
"plt.plot(t,samplegp(1.0), label=\"tau=1.0\")\n",
"plt.plot(t,samplegp(3.0), label=\"tau=3.0\")\n",
"plt.legend()\n",
"plt.savefig(\"gp.png\")"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3.8.8 ('base')",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.8.8"
},
"orig_nbformat": 4,
"vscode": {
"interpreter": {
"hash": "72bc7f8b1808a6f5ada3c6a20601509b8b1843160436d276d47f2ba819b3753b"
}
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment