{ "cells": [ { "cell_type": "markdown", "metadata": { "Collapsed": "false" }, "source": [ "# Timber Harvesting -- using 2 nodes\n", "\n", "**Randall Romero Aguilar, PhD**\n", "\n", "This demo is based on the original Matlab demo accompanying the Computational Economics and Finance 2001 textbook by Mario Miranda and Paul Fackler.\n", "\n", "Original (Matlab) CompEcon file: **demdp01a.m**\n", "\n", "Running this file requires the Python version of CompEcon. This can be installed with pip by running\n", "\n", " !pip install compecon --upgrade\n", "\n", "Last updated: 2022-Oct-22\n", "
" ] }, { "cell_type": "markdown", "metadata": { "collapsed": false }, "source": [ "## About\n", "\n", "Profit maximizing owner of a commercial tree stand must decide when to clearcut the stand.\n" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "Collapsed": "false" }, "outputs": [], "source": [ "import numpy as np\n", "from compecon import NLP\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "Collapsed": "false" }, "outputs": [], "source": [ "price = 1.0 # price of biomass\n", "κ = 0.2 # clearcut-replant cost\n", "smax = 0.5 # stand carrying capacity\n", "ɣ = 0.1 # biomass growth parameter\n", "δ = 0.9 # discount factor" ] }, { "cell_type": "markdown", "metadata": { "Collapsed": "false" }, "source": [ "### Code the growth function" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "Collapsed": "false" }, "outputs": [], "source": [ "snodes = np.array([0.2, 0.4])" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "Collapsed": "false" }, "outputs": [], "source": [ "def h(s): return s + ɣ*(smax - s)" ] }, { "cell_type": "markdown", "metadata": { "Collapsed": "false" }, "source": [ "## SOLUTION\n", "\n", "### Code the approximant and the residual" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "Collapsed": "false" }, "outputs": [], "source": [ "def vhat(c, s): return c[0] + c[1]*s\n", "def vhat1(c,s): return price*s - κ + δ * vhat(c,h(0))\n", "def vhat0(c,s): return δ * vhat(c, h(s))\n", "\n", "def resid(c,s=snodes): return vhat(c,s) - np.maximum(vhat0(c,s), vhat1(c,s))" ] }, { "cell_type": "markdown", "metadata": { "Collapsed": "false" }, "source": [ "### Solve collocation equation" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "Collapsed": "false" }, "outputs": [], "source": [ "cc = NLP(resid).broyden(np.zeros(2))" ] }, { "cell_type": "markdown", "metadata": { "Collapsed": "false" }, "source": [ "### Compute critical biomass" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "Collapsed": "false" }, "outputs": [ { "data": { "text/plain": [ "0.36199999999813565" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "scrit = NLP(lambda s: vhat0(cc,s)-vhat1(cc,s)).broyden(0.0)\n", "scrit" ] }, { "cell_type": "markdown", "metadata": { "Collapsed": "false" }, "source": [ "## ANALYSIS" ] }, { "cell_type": "markdown", "metadata": { "Collapsed": "false" }, "source": [ "### Compute refined state grid" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "Collapsed": "false" }, "outputs": [], "source": [ "ss = np.linspace(0,smax,1000)" ] }, { "cell_type": "markdown", "metadata": { "Collapsed": "false" }, "source": [ "### Plot Conditional Value Functions" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "Collapsed": "false" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Optimal Biomass Harvest Level = 0.3620\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAuwAAAGZCAYAAADSJ2QTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAACLwElEQVR4nOzdd3hUZdoG8PtMn0kBQkJC6AJJSEKaUkWpAiJNERSRjgjRtX4WdtXtll07SBDBQpOigguIFFEUURASWhodAiG9TqbPnO+PmBPGZEJImwm5f9fltctzzpl5JiflmXee930FURRFEBERERGRR5K5OwEiIiIiInKNBTsRERERkQdjwU5ERERE5MFYsBMREREReTAW7EREREREHowFOxHRDeDCWuRp+D1JdPNjwU5EDcpgMOCLL77ArFmzMHjwYERGRmLgwIF49NFH8d1337k7vSrGjh2L0NBQXL58WYqFhoYiNjbW6byzZ89i5syZKCwslGIHDx5EaGgoHn300SbL90YtXrwYoaGhWLlyZY3nzZ07F6GhoUhISLjuY5aUlCAyMhK9e/dGUVHRDeUzbNgwhIaGoqCg4IauawwvvvgiQkNDr/vfhAkT3J1qtYqLi/G3v/0Ne/bscYpX9/1LRM2bwt0JENHN4+TJk3jyySdx+fJl+Pj4ICQkBFFRUbhy5Qr27duHH374Affeey9effVVyGTNa7xg/vz5TkX9zebee+/F/v37sX37dixcuLDGc3fs2AGr1YrRo0ejdevWTZNgI4qNjUXHjh1dHu/QoUMTZlN7r7/+Or766iv079/f3akQUSNjwU5EDeLMmTOYNm0aTCYTFi5ciHnz5sHb21s6npqaiscffxybN2+Gt7c3XnrpJTdmW7NvvvmmyhuK6toOoqKi8M0338DLy6upUms0d911F3x8fHD69GmcOnUKISEhLs/93//+BwC47777miq9RjVlypRm+VpctcJU9/1LRM0bf6KJqN5EUcRzzz0Hk8mEp556Ck899ZRTsQ4AvXr1wvLly6FQKLBu3TqPHq3u3r07unXrdt3ztFotunfvjqCgoCbIqnGp1WqMGTMGALB9+3aX5125cgVHjhxBu3btMGjQoKZKj25Abb9/iaj5YMFORPX222+/ISUlBUFBQZg7d67L87p3745JkyZh8ODByM7OdjqWl5eHf/3rXxg2bJjU9/7kk08iLS2tyuOEhoZi+vTpyM/Px1/+8hfcfvvt6N27N8aPH48NGzZU+9ynTp3CE088gYEDByI2NhYLFizAuXPnqj332h7gij71K1euAAAGDBiA0NBQp2PV9bDv2rUL06dPR1xcHKKiojBhwgR8/PHHsFgsTudV9JgfPHgQ33zzDSZNmoTo6Gj069cPzz77LDIyMqo8dllZGZYtW4ZJkybhtttuQ2RkJO644w48++yzOH36dLWvqTYmTpwIoHyE1pWtW7dCFEVMmDABcrkcQPnXdtGiRRg+fDh69+6NmJgYjBs3DkuWLIHJZLru89bUc13dHAMAyM/Pd/p+GTRoEBYtWiTdp8ZSUw/+o48+Kt3La88fNmwYysrK8Prrr2PIkCGIjIzEqFGjsHz5cthstiqPYzQasWzZMowbNw4xMTEYPHgwnnjiCZw5c0Y6JzQ0FJs3bwYAPPnkk07P6+rrmZKSgieffBIDBgxAZGQkhg8fjn//+9/Iz893Oq/i+3rx4sVITk7GI488gttuuw2xsbGYMWMGDh06VOWxr1y5gkWLFmHkyJHo3bs3Bg4ciIULF+Lw4cO1/MoSUU3YEkNE9bZjxw4AwIgRI6BSqWo89x//+EeV2Llz5zB9+nTk5eWhS5cuGDZsGDIzM/Htt9/iu+++w9tvv42RI0c6XVNUVIQHHngAhYWFiIuLg8lkwm+//YZXXnkFFosF06dPl85NTEzE3LlzYTAYEB0djaCgIPz222+YOnWqVHS64u/vj3HjxuG7776DwWDAqFGjrvsaX331VXz22WdQqVTo06cPtFotfvvtN7zxxhv47rvvsHLlSmg0GqdrPvvsM3z33XcIDQ3FnXfeiaSkJGzbtg2//fYbdu/eDbVaDaC8WH/wwQdx6tQpdOjQAf369YPZbMbJkyexbds2fP/999i2bRuCg4NrzLE6cXFx6Nq1Ky5cuIDjx48jKiqqyjlbt24FUN7zDgA///wzFixYAKvVipiYGERERCA3NxdHjx7FqVOnkJycXKuJrDfi4sWLmD59OrKzs9G5c2cMGTIEmZmZ+Oqrr7B37158+umn6NWrV4M+Z31YLBbMmjULp06dQmxsLLp3746DBw/irbfeQkFBAV588UXp3OLiYsycOROpqakICAjAHXfcgZycHOzcuRM//vgj1q5di4iICIwbNw5Hjx5FRkYGbrvtNrRv3x7+/v4uc9ixYweee+45WK1W6WcgJSUFq1atws6dO/HZZ59VGZU/evQoli9fjrZt26Jfv364dOkSDh48iMTERGzYsAEREREAgIKCAsycORMZGRno1asXhg4dipycHOzduxf79u3DsmXLcOeddzbOF5eopRCJiOpp5syZYkhIiLh58+YbvtbhcIjjx48XQ0JCxLffflu02+3SsW+//VYMDw8XY2JixKysLCkeEhIihoSEiJMnTxZzc3Ol+JYtW8SQkBBxxIgRUsxms4ljxowRQ0JCxI0bN0rx0tJS8eGHH5YeKyMjw+nxY2JinPIcOnSoGBISIubn50uxX3/9VQwJCRHnz58vxXbu3CmGhISIQ4YMEc+dOyfFS0pKpOf797//LcXff/99MSQkRAwNDRW/+OILKa7X68UJEyZU+bomJCSIISEh4qJFi5y+Vnq9Xnr8pUuXVnn8FStWVPflr6Li8V999dUqx06ePCmGhISIDzzwgBQbOXKkGBYWJh44cMDp3OPHj4uRkZFiSEiIePXqVSle3dexuq93hXvuuafK/bn//vvFkJAQccmSJU5fg82bN4shISHiqFGjRJvNdt3X+sILL4ghISHil19+ed1za8q/wvz588WQkBDx119/rXL+XXfdJV64cEGKHzp0SAwLCxNjY2PFsrIyKf7KK6+IISEh4p/+9CfRZDJJ8S+//FIMCQkR77333ir579ixwymPP349s7KyxN69e4sRERHid999J8VtNpv4+uuviyEhIeLEiRNFh8MhimLl93VISIj48ssvi2azWbrmL3/5i/T9V2HJkiViSEiI+O677zrlUXE/HnrooRq+okRUG2yJIaJ6y8nJAQC0bdv2hq89ePAg0tLSEBERgaeeesppstyoUaMwdepUGAwGrF+/vsq1zz//vNOo4vjx4+Hr64tLly5Jyy8eOXIEZ86cwYABAzB58mTpXG9vb7z66qsQBOGGc67JqlWrAAAvv/yy04ilj48P3n77bahUKmzYsAEGg8Hpuv79+2PSpEnSv728vKQWlRMnTkhxrVaLO++8E08++aTT18rLywvjx48HAFy9erXO+U+cOBEymQw7duyAw+FwOlYxul4xQbO0tBQxMTGYNm0aBgwY4HRu79690bt373rn80dHjhzB8ePHcdttt+Gxxx5z+hpMnDgRw4cPx/nz5/HTTz/V+jEXLVrkcknH+Pj4Bsn7scceQ5cuXaR/9+nTByEhISgrK5Nas8xmM77++mvodDq8+uqr0qcqQPnXfODAgdBoNCgtLb2h596wYQPMZjMefvhhDBs2TIrL5XI8//zz6NWrF1JSUvDrr786XafT6fDnP//Z6ROladOmAXD+nszNzQUABAQEOF0/YcIE/OUvf8H8+fNvKF8iqootMURUbwpF+a+S6vpxr6eix3XkyJHVFs9jxozB6tWrq/TCCoJQpWVDEAS0bdsWJSUlMBqNaNOmjXRddRMkO3XqhJCQEKSnp99w3tWx2Ww4evQoNBpNtS0AAQEB6NOnD37++WecOHEC/fr1k45FR0dXez4Ap+J+5syZmDlzptN5JSUlSEtLw2+//QYAsFqtdX4NQUFB6N+/Pw4cOIDDhw+jb9++AACHw4Ht27dDo9FIk1N9fHzwxhtvOF0viiKuXLmCEydOIC8vr975/FHFa7z2a3etQYMG4bvvvsOhQ4cwZMiQWj1mTcs6RkZG1inPP6ru/la82TQajQDKi2Cj0Yg77rijyqRtAPjkk0/q9NzX/oz9kSAIGD16NFJTU3H48GGnN16hoaFVWrcqcr72e7JPnz74/PPP8e9//xupqakYPnw4+vXrB61WixkzZtQpZyJyxoKdiOotICAA6enpddoMp2J0zlXPdcUa2BXnVdBoNNX2kle8eagYHa64LjAw0OXjN1TBXlRUBKvViuDgYCmP6p4PgFTMVvDx8alybkV/vfiH5fuys7OxZs0a/Pbbb7hw4YL0aULFG54/nn+j7r33Xhw4cADbt2+XCvZffvkFOTk5GDduXJVi8uDBg9i4cSPS09Nx6dIlmM3mBs3nWllZWQCADz74AB988IHL8/44qbkmTbGso6+vb5WYq+/V9u3bN+hzVzyuq/XkXf2MVZdzdd+T99xzD44fP45Vq1Zh48aN2LhxI1QqFQYMGICJEyfi7rvvbvBPsohaGhbsRFRv4eHh2L9/P44fP+7U1lGdpKQkHD58GHfccQfCwsKkP/yu/qBXHP9jcd5QBYCrwrourvdagMrirK6v55dffsGCBQtgMpnQsWNH9OvXDz179kTv3r2rTGCsq5EjR+Lvf/87du7ciZdffhkKhUJae/2P9/eVV17Bhg0boFAoEBERgQkTJiAkJAS33nor3n33Xezbt69eudjtdqd/V3z94uLiatzQqKFGxm/EH3O9URWvraGL2+u9Yarv9yRQ3lY0ffp0fPvtt9i3bx+SkpKwb98+7Nu3D9988w2WLFly44kTkYQFOxHV29ChQ7F8+XLs27cPFoulxlVU1qxZg23btuHixYv417/+hXbt2gGAy3XZK+J16Y8HKkfWMzMzqz3+x1HF+mjdujWUSiWysrJgs9mqfTNQ8Xr8/Pxu+PFFUcTLL78Mk8mEN998E+PGjXM6vnHjxrol/gcajQZ33303Nm3ahF9++QV9+vTB7t27ERwc7LSr5qFDh7BhwwZ07doVK1eurNJWUlJSUqvnEwTBZbH7x37tijahESNG1LiEaGOpKGL/2N8PVM31RlW0m1R8ivBHBw8eRG5uLgYMGHBDPw/t2rXD+fPncfny5Wo/aarvz1iFjh07Yt68eZg3bx4MBgN2796Nv//979i9ezeSkpJcLt1JRNfHSadEVG+xsbGIiYnB1atXa+yzPX78OHbt2gUAeOCBBwAAt956KwBg9+7d1Y4EfvvttwCA2267rU65VRSY3333XZVj+fn5OHnyZJ0etzpKpRLR0dEwmUz48ccfqxzPy8vDkSNHoNPppCXxbkRBQQEyMjLQsWPHKsU6UL7EIlB9MXmjKia87t69G/v370dZWRkmTpzoNOp67NgxAMC4ceOqFOsFBQVITU2tVT46nQ5msxl6vd4pfuXKlSpvqCq+X1xNKn3nnXdw3333SRNkG5pOpwOAKmuXG43GerdWRUREQKVS4ciRI1UmJQPA+++/j2effVZ6I1TbEfBrf8b+SBRFKd6nT5865f3nP/8Z/fv3d5pcrNPpMGHCBAwfPhxAw048JmqJWLATUb0JgoBXXnkFSqUS77zzDt5//32UlZU5nXP48GE89thjsFgseOCBB6QVRPr374+QkBAkJyfj3XffdSru9uzZg88//xw6nU4qIG9UdHQ0oqKicOzYMSxbtkyKm81m/OUvf6n1hMiKFTv+WFT+UcX67//85z9x4cIFKa7X6/F///d/sFqtuO+++6pM5quN1q1bQ6PRICsry2lDKYfDgU8++UR6c1PRQ14ft912G7p06YK9e/di9+7dEAShSp93xQ6v+/fvd9oQqqCgAE899ZS0adIfN4v6o5CQEADA6tWrpZjRaMTf//73KucOGDAAPXv2xC+//IKlS5c6fb/89NNP+Pjjj5Gamip9fzW06nK12Wz417/+JU0erStvb2+MHz8eer0ef//7352+N7ds2YLDhw8jIiJCWn2o4pOs643sT5kyBRqNBmvWrMHevXuluMPhwFtvvYXU1FSEhoYiLi6uTnkHBASgsLAQb775plPOubm5OHToEGQymVtalIhuJmyJIaIGERERgY8++giPPfYYPvjgA6xatQoRERFo3bo1zp8/L40+jh8/Hi+//LJ0nSAIePvttzFz5kwsW7YMO3bsQK9evXD16lUcO3YMarUar776ao39yjURBAGvv/46ZsyYgXfeeQfffPMNunbtiqNHj6K0tBQ9e/as1e6gXbp0wblz5/Doo4+iZ8+eeP3116s9b/To0Xj44YexZs0ajBs3Dn379pU2TioqKsJtt92G//u//6vTa5HL5Zg2bRpWrlyJyZMno2/fvtBoNEhOTsbVq1fRo0cPnDlzpsqE1rqaOHEi3nvvPWzbtg19+vRBp06dnI4PGzYMHTt2RFJSkrTDpV6vx5EjR2C1WnHLLbfg3Llz1207mj17NpKSkvDuu+/i+++/R0BAAI4cOQKFQoEBAwbgl19+kc4VBAFvvfUWZs2ahffeew9ffPEFevXqhby8PBw9ehRAeV99165dG+Rr8EfTp0/Hzp07sWnTJhw7dgxdu3bFiRMnUFxcjOHDh1f7Sc6NeP7553H8+HFs2bIFBw8eRFRUFDIzM3HixAl4eXnhv//9r3RuxWt866238NNPP2H27NnVtp20b98er776Kp5//nksXLgQ0dHRaN++PVJSUnDp0iUEBQXhnXfecVoi80bMmzcP3377LbZt2ya9qbBYLNInBXPnzkXnzp3r9NhEVI4j7ETUYAYMGIBvvvkG8+bNQ8eOHXH8+HHs3r0b+fn5GD58OJYvX47//ve/UCqVTtf17NkTmzdvxkMPPQSr1YrvvvsOV69exYQJE7Bp0ybcfffd9cqre/fu2LhxI+677z4UFBRg37596Ny5Mz777LNaFxIvvvgi4uLicOXKFRw8eBAZGRkuz3355ZfxzjvvICoqCklJSfj555/RoUMHvPTSS/jss8+g1Wrr/FqeeeYZvPjii+jSpQuOHDmCo0ePIjAwEIsWLcKWLVvg7++Po0eP1mnFnj+qWJPdZrNJO5tey8vLC6tXr5bWf9+3bx8uX76MQYMGYfXq1dIE2O+//77G5xk1ahQ++OADxMTEID09HYcPH8btt9+OL774QhrFv1ZoaCi2bNkirQm+b98+ZGZmYvDgwVi1apUUbwwxMTH49NNPMWDAAFy+fBkHDhxAWFgYNm7c2CC7q7Zq1Qqff/45FixYAJVKhb179yIjI0OaU9C9e3fp3AcffBBjx46FxWLBTz/9hFOnTrl83HvuuQfr16/HyJEjcenSJezduxcymQzz5s3D5s2bnR73Rvn4+GDNmjXSzsE//vgjkpKSEB4ejjfffBPPP/98nR+biMoJYkOut0VERERERA2KI+xERERERB6MBTsRERERkQdjwU5ERERE5MFYsBMREREReTAW7EREREREHozrsF+Hw+GA3c6FdIiIiIio8SiVcpfHWLBfh90uoqio6hbRREREREQNJSDAx+UxtsQQEREREXkwFuxERERERB6MBTsRERERkQdjwU5ERERE5MFYsBMREREReTAW7EREREREHowFOxERERGRB2PBTkRERETkwViwExERERF5MBbsREREREQejAU7EREREZEHU7g7ASIiIiIitxJFqC59D3X6l7D7hcFw25/cnZETFuxERERE1DLZrVCf2QpdUgIU+am/B7+GIXYBIFe6NbVrsWAnIiIiopbFaoA25XNoj30EeellKexQesHQ7zmPKtYBFuxUBw6HA99/vwfbt2/FhQvnUFhYAF/fVoiOjsWUKVMRGRnl7hSJiIiIqhCMBdCe+ATaE59CZiqU4g5tAAzRc2GKnA5R3cqNGVZPEEVRdHcSnsxqtaOoyODuNDyGXq/HX//6Zxw8eAAxMXEYOHAQfH1bISvrKnbu/AZXr2bi6aefw6RJD7g7VSIiIiIAgKwkA7qjH0KTuh6CzSTFba26whizAKaw+wGFxo0ZAgEBPi6PsWC/Dhbszv761z9j797d+POf/4q77x7rdMxsNuOFF57G0aOJWLVqAzp37uKmLImIiIgAeV4KdIlLoT6zFYJol+LWdtEwxC6E5Za7AZncjRlWqqlg57KOVGvHjx/Fd9/twqhRY6oU6wCgVqvxf/+3CDabDd98s9UNGRIREVGLJ4pQXv4ZrbZOg9+GkdCc3iIV65bOg1E0YQOK7t8GS4+xHlOsXw972KnWdu36FgAwffpsl+d07NgJ772XgPDwSCl2//3j0KdPPzgcDuze/S1atWqNjz9eizZt2uDYsSR8/PFHSEk5AQDo1SsCc+bMR0xMHABg1qyHAACffrpOerwvvliPd999E48//hQefPBhKT5r1kPw9/fHm2++33AvmoiIiJoHhx2q899Cl7gUypxjUlgUZDD3GAdDbDzsARFuTLDuOMJOtXb06BG0beuPLl261njerbf2gVardYrt2bMTp0+fwhNPPItx4yaiTZs22L9/H/70p0eRnZ2FmTPnYubMucjOzsKTTy7E/v37AAADBtyOs2dPo6ioSHqspKQjAIBjx45Ksfz8PJw9exoDBgxqkNdKREREzYTNBE3yGrRZNwStvn1UKtZFhQbG3rNQ8PB+lI78oNkW6wBH2OkG5OTkoEuXqn3pJpMJJpPJKSaTyeDr6yv922w245//fB0dOnQEANhsNrz99n/g7x+AlStXwcvLGwAwYcIkzJjxAN566w307387+vcfiNWrP0FS0mEMHToCoiji6NFEBAS0w4kTRyGKIgRBwKFDv0IURQwcyIKdiIioJRDMxdCcXA3dsZWQGXOluEPdGsbes2CMmg1R29aNGTYcFuyNKPlqCVb8egkGi/36JzcRnUqOef07I6K97/VP/gNRdKC6KcorVizD+vVrnGJBQe3xxReVfewdOnSUinUAOHUqDTk52Vi48E9SsQ4APj4+uO++KfjwwyVIS0tBRERveHt748iR8oL9zJnTKC4uRnz8k1i69D1cuHAe3brdgoMHf0HXrregffvgG35dRERE1HzI9FehPbYCmuS1kFn1Utzu3QHGmEdg7DUVUHm5McOGx4K9EX2eeAX7zxW4O40qvFRy/OueGy/Y27YNQEFBfpX4hAn3oV+/AdK/lyx5F2Vleqdz2rTxc/r31auZAFDtSjJdu3YFAGRlXUVkZBRuu60fEhN/AwAkJR1GmzZ+uOeecUhIeB/HjiWia9duOHz4IMaMGXfDr4mIiIiaB3nhGWiTEqBJ/wqCwyrFbW3DYIhdCHOP8R634VFDYcHeiKbGdUCZxe5xI+xTb+14/ROr0bt3FHbs2IbLlzPQsWMnKd6pU2d06tRZ+rePj0+Vgl0mc54uUdNqog5H+TGlsvyHrn//gfjhh++Ql5eLxMTDiImJQ6tWrXHLLT1w7NhRhIWFo6ioCP37316n10VERESeS5F1pHxpxvM7neKW4P4wxi6EpcswQBDclF3TYMHeiCLa++KdeyOvf2IzMXLk3dixYxs2bFiHZ599oV6PFRRU3rpy8eIF3HGH87FLly4CANq1CwQA9O9/OwRBwG+/HcTx48cwd+6jAIDY2Dj8+OMP6Nq1G7y9vREVFVOvnIiIiMhDiA6oLu6FNjEBqqsHK8MQYLllFAyx8bAFxbkxwabFVWKo1vr06Yfhw+/Cli1f4KuvNlV7zq5dO5CennrdxwoL64W2bf2xefMXTqPxZWV6bN68CW3b+iM0tBcAwN/fHz169MRXX21ESUkxYmPLf0BjYuKQk5ON7dv/hz59+kOh4PtPIiKiZs1uhTrtC7RZfxdabZ8lFeuiTAVj+FQUPvQDSu5e0aKKdYAj7HSDXnjhJdjtdrz99hvYvv1/uOOOwfD390dOTg727fseZ8+ehp9fW/zpT0/X+DgKhQJPP/0cXnllEebOnYFx4yYAALZu/Rr5+Xn45z/fcGqj6d//dqxe/QlatWqFbt26AwBiYm6FIAjIzLyC2bMfabwXTURERI3LUgZtyjpoj30EuT5TCjtUPjBFPAxj9Fw4vILcmKB7sWCnG6LTeeFf//oPfv75J2zf/j9s3boFBQX58PLyQo8eIXjmmRcwZsw4aDSa6z7WkCHD8fbbS/DppyvwyScfQaFQIDw8EosWvYzo6FincyuWd4yKioXwe59a69at0bVrN1y4cB79+w9slNdLREREjUcw5EF74hNoT3wKmblYitt17WCMngtTxHSI6htfKONmI4g1zf4jWK12FBUZ3J0GERER0U1DVnwRuqMfQpO6AYLdLMVtrW+BMXYBTKGTALnajRk2vYAAH5fHOMJORERERE1CkXsS2sSlUJ/dBkF0SHFrYCwMcfGwdB0JyORuzNAzsWAnIiIiosYjilBe3g9dUgJUGT86HTJ3HgpjXDyswf1v+qUZ64MFOxERERE1PIcd6rPfQJuUAGXucSksCnKYe06AIXYB7P7hbkyw+WDBTkREREQNx2aEJm0TdEkfQl5yUQqLCi2M4VNhjJ4Ph2/dNnFsqViwExEREVG9CaYiaE+ugvb4x5AZ86S4Q9MGxqg5MPaeBVHTxo0ZNl8s2ImIiIiozmSlmdAe+wja5LUQbJUr69l9OsEQMx+mXg8CSq0bM2z+WLATERER0Q2TF5yCLikB6lObIThsUtzWNhyGuHiYe4wFZCw1GwK/ikRERERUa4rMQ+WF+oXdTnFLh4EwxMXD2mkwV3xpYCzYiYiIiKhmogOqC3ugS1wKZdbhyjAEWLqPgSF2IWyBMe7L7ybHgp2IiIiIqme3QH1qM3RJy6AoPC2FRbkaptD7YYx9FPbWt7gxwZaBBTsREREROREspdAkr4P22EeQl2VJcYfKF6bIGTBEzYHo1c6NGbYszaZgz8jIwBtvvIFDhw4BAIYMGYIXX3wRfn5+NV73yy+/4P3330daWhq8vb0xevRoPPXUU/Dy8mqKtImIiIiaDcGQC92xldCcXAWZpUSK270CYYx+BKaIaRBVPm7MsGUSRFEU3Z3E9RQWFmLSpEmwWCyYMWMG7HY7Vq5ciQ4dOmDTpk1QqVTVXvfrr79i9uzZiIiIwL333ourV69i1apViIiIwNq1ayGTya773FarHUVFhuue19IYDGX4+uvN2LNnJy5fvgS73Y5u3W7B2LETMW7cROlre//94xAU1B5Llix3c8bOrFYrduzYhh07tuLSpYswGo3o1Kkz7rprNKZMecjl99T1GAxlMJstaNOG68wSEVHzISs6D93RD6FJ2wTBbpbitjY9YIhdCHPIRECudl+CLUBAgOs3Qs1ihP3TTz9FVlYWtm7diu7duwMAoqOjMXv2bGzZsgVTpkyp9rr//Oc/aN++PdasWQONRgMAaN++Pf7xj3/gp59+wuDBg5vsNdxMLl26gBdeeAZXr2Zi5Mi7MWbMOFgsFuzfvw///e+rOHYsCS+//A8IHjpDPD8/D4sW/R9SU5MxePAwDB06AjKZDIcOHcSyZUtw6NCvePPN92+4aE9LS8WLLz6DV175J9q0ua2RsiciImo4ipxj0CUuhersNxBQOYZrDboVhth4WLrdBQjXH+CkxtUsCvbt27ejb9++UrEOAAMHDkS3bt2wffv2agt2s9mMNm3aYOTIkVKxDgB9+/YFAKSnp7NgrwOz2YwXX3wWxcXFWLFiNXr06Ckdmzr1Ybz11hvYvHkTevWKwOTJD7ox0+qJooi//e0vOHPmNN5+ewn69OknHbv//gexdu1nSEhYjA8/XII//emZG3rsc+fOIC8vt6FTJiIialiiCGXGj+WF+pWfnQ6Zu46AITYetuC+bkqOquPxb5mKi4uRkZGBiIiIKsciIiJw8uTJaq9Tq9VYuXIlFixY4BRPTU0FAAQHBzd8si3A5s2bcOnSRfzpT087FesVHn/8Sfj4+OLrr790Q3bX9+OP3yMp6QimT5/lVKxXmDZtJrp06YodO7bDbDa5IUMiIqJG4rBBfWoLWm8cjdZbp0nFuihTwBQ2GQUP7kHJPZ+yWPdAHj/Cnp2dDQAIDAysciwgIAB6vR6lpaXw8al5AsSVK1dw8OBBvPHGGwgJCcFdd93VKPne7L77bhe0Wh3uumt0tcfVag2WL/8UQUHtXT7GyZPHsWLFMiQnl7/ZiozsjUceWYjw8EjpHFEU8fXXX2L79v/hwoULsNttCApqjzFjxmHatJlSu839949Dnz794HA4sHv3t2jVqjU+/nityx7yPXt2AQDGjbvXZX5vvvk+2rTxg1pd/snM44/PR1bWVXzxxVan866Nr1z5IT755CMAwBNPLEBQUPsq5xMREbmF1QhN6nroji6HvDRDCosKHYwR02CMfgQOHw5kejKPL9jLysoAAFqttsoxtbp88oPBYKixYC8qKsKwYcOkx3nppZeka6n2RFHEqVPp6N07GgqF62+dTp06uzz222+/4rnnnkLPniF45JEFsFgs+OabrXj88fl4550PEB0dCwD46KMErFr1Me6+eyzGjbsXBkMZvv32GyxbtgRt2/rj7rvHSo+5Z89OdOrUBU888SwKCvJrnPCZnp6KoKD28Pf3d3lO+/Y3/ktr8OBhyM/Pw//+txnTp89Gr15VPxEiIiJqSoKpENoTn0J7/BPITAVS3KFtC2PUHBgjZ0DUcJGE5sDjC3aHw3Hdc6632osgCHjnnXdgsViwevVqzJ49G2+//TZGj65+lJiqV1RUBLvdjrZtXRe7NXE4HPjvf19Dr14RWLJkOeRyOQBg0qQHMHv2Q3j33f/ik0/WwWaz4csvN2D48JH4y1/+Jl0/btxEjBs3Ej/88J1TwW42m/HPf76ODh06XjeHgoJ8dO9etZWnvnr06InIyCj873+b0adPP8TFcdIpERG5h6z0CrRHl0Obsg6CzSjF7b5dYIh9FKawyYCi6kAoeS6PL9gr1ks3m81VjlXErremeqtWrTBmzBgAwOjRozF27Fi8/vrrjV6wK7KToDv8HgSLvlGf50aIKm8YbnsStsDYG75WLi9/Y1SbN1HVOXUqHZmZVzBx4v0oLS11Onb77Xdgw4Z1yMnJRrt2gfjf/3bBbrc5nVNUVASdzgtGo9Ep3qFDx1oV6wAgk8nhcNjrlD8REZEnk+enQpeYAPXpryGIlX/rrP6RMMbFw9x9DCDz+NKPquHxd61icmhubtXVN3JycuDr6wudTlfrx9NoNBgyZAhWr16NgoKC6268VB/aYyugvrCn0R6/rkSlN0pHLrnh63x8fKFUKlFYWHD9k6tx5cplAMDSpe9h6dL3qj2nomBXKpX45Zf9+Omnfbh06SIuX85AaWn5Bg5/fMPQpo3zPTSbTdDrnd8kabU66HQ6+Pm1RWFhYZ3yJyIi8jiiCGXmr9AmJUB9ca/TIUvHO2CIi4e14yDAQ5daptrx+ILd19cXHTt2RHJycpVjKSkpiIyMrOYq4OzZs3jkkUcwd+5cTJs2zelYWVkZBEGo8+Y4tWWMngfBWuZxI+zG6Hl1ulYQBERE9EZ6ehpsNpvLPvbly5fiypXLeOIJ52URK0a2581bgIiI3tVe27lzV4iiiH/842Xs2bMTUVEx6N07ChMm3IeYmDg88cSCKtf8sSXqu+9249VX/+4Umz37Ecyd+yh6947Cjh3bkJeX57KP/csvN+Lw4UN49NHH0LVrt+q/GKj7Jw1ERET1JjqgOr8TusSlUGYnVYYFGczdx8IYtxC2gOr/1lLz4/EFOwCMHDkSq1atwtmzZ6W12A8cOIDz589j7ty51V7TpUsXlJaWYv369Zg8ebJUnF+5cgW7du1Cnz594O3t3ah52wJjUXLPp436HE1t8OBhOHo0EXv37sbIkXdXOW42m7Bt29dwOOxo1aq107GKyZw6na7KkoqpqckoKSmBWq3GsWNJ2LNnJ2bNmod58yoLdJvNhpKSYgQHd6gxx759B+Cddz5wilVcM3jwUOzYsQ3bt3+NmTOrfu84HA5s27YF58+fw3PPLQIAyOVyWK2WKucWFOTXmAcREVGDs5uhSf8S2qRlUBSdk8KiXA1TrwdgiJkPR6uu7suPGkWzKNgfeeQRfP3115g1axbmzJkDs9mMFStWICIiAhMmTAAAZGRkIDExEXFxcejUqRMUCgVeeuklPP/885g+fTrGjx+PwsJCrF27FoIg4OWXX3bzq2qexo+/Fxs2rMWSJe+iR4+euOWWHtIxu92ON998HQUF+YiPf7LKCHxYWDjatvXHpk0bMHbsRKmVqaxMj1deWQSLxYIvv9yG4uJiAKgyur116xaYTCbY7TX3oPv7+7scPR80aDAiInpjzZpP0bt3dJXJoR9/vBynT5/CffdNhp9fWwCAn19bJCUdQV5eLvz9AwCU72p6+XKG0/KVFSP9oiiCiIioIQnmEmiS10B7bCXkhmwp7lC3gjFyJoxRcyDq6rYoBHm+ZlGw+/n5Yc2aNXjttdfw/vvvQ6PRYMSIEXj++eelkfPffvsNixYtwmuvvYZOnToBACZMmAClUokVK1bgtddeg06nQ//+/fH000+jWzfXrQ7kmlqtxquv/hfPPPM45s2biZEjRyMsLBwlJcX4/vs9OH36FIYOHYEHH5xW5VqFQoGnn34Or7yyCHPmPIxx4yZApVJj69bNyMq6ilde+ScUCgV6946Cl5cXFi9+G9nZWfD29kFi4mHs3bsbKpUaBoOhXq/h739/FU8+uRBPP/0YhgwZhqioWJjNJhw4sB9HjyYiOjoWCxc+IZ0/YsQo7N79LZ599gnce+8kFBQU4MsvN6Bjx86w2azSea1bly+NtXnzF8jPz8fIkVyFiIiI6kdWlg3t8ZXQnFwNmaVywQa7d3sYo+fDFD4VoqpxOwbI/QSRw4E1slrtKCqqX4F4M8rNzcHGjZ/jl19+RnZ2FkTRge7de2L8+HsxZsw4p42NgoLaY8mS5dK1hw8fwqpVHyM1NQUymYBbbumOhx+ejdtvv0M65/jxo0hIWIwzZ05DpVKiU6cumDx5KlJSTuKLL9Zj8+Zv4OfXttrHr42SkmJ89dUm/PDDXmRlZcJisaJLly4YNWoM7r//wSqfDnz11SZs3Pg5srOvolOnzpg16xEcPXoEBw7slzZIstls+Mc/XsbPP/8IlUqNLVt2cL1/IiKqE3nROWiTEqBJ+xKCo7It0+YXCkPsQph7jgfkjTsXr6UoNdmQkl2KlKxSnM0rQ2zHVpgU3fQbSQUEuN5TiAX7dbBgJyIioqaiyEqELikBqnPfQkBliWZt3xeGuHhYugwDhJr3nyHXTFY70nP0SMnWI/lqCVKz9bhU6LxctEwAfnpiEFSKpv0611SwN4uWGCIiIqKblihCdel7aBOXQpX5q9Mhc7dRMMQuhK09N+S7UTa7A2fzDUjOKh89T8kqxbm8MthrGKoO8lHjvuj2TV6sXw8LdiIiIiJ3sFuhPvM/6JISoMhPk8KiTAlTyH0wxi6A3a/hd+e+GTlEEZcKjVJhnpKlx6lcPcw210swt9EqER7kg/Ag79//1wd+Os9sM2LBTkRERNSUrAZoUz6H9uhyyPVXpLBD6QVTxMMwRs+Fw7vpe6ibC1EUkV1qRkq2XirQU7NLoTe7XkXOSyVHWKA3wgPLC/OI9j4I8lFLc+48HQt2IiIioiYgGAugPf4xtCc+hcxcJMUd2gAYoufCFDkdorqV+xL0UEVGq1SYV7S3FBisLs9XygWEBJSPmkf8PnLexU8LWTMpzqvDgp2IiIioEclKMqA7+iE0qesh2ExS3NaqK4yxC2AKvR9QaNyYoecwWOxIza5sa0nJLkVmscnl+TIBuKWtl1NbSw9/LyjlntWDXl8s2ImIiIgagTwvBbrEpVCf2QpBrGzXsLaLhiF2ISy33A3I5G7M0L0sNgdO55Vd03deigsFBjhqmBTasbVGamsJD/JBWKA3tMqb/2vIgp2IiIiooYgilFcOQJe0FKpL+5wOWToPhiE2HtYOA4Fm3J5RF3aHiAsFBqfWljN5ZbDWsGSLv5fKaVJor0AftNYqmzBrz8GCnYiIiKi+HHaozu2ALikBypxjUlgUZDD3GAdDbDzsARFuTLDpiKKIK8Ump7aW9Gw9DFbXk0J91IrKtpbfR9Db+XDzwQos2ImIiIjqymaCJv0LaJM+hKL4vBQWFRqYej0IQ8x8OHw7uzHBxpdXZnFqa0nJKkWxyebyfLVChrB2lT3n4UE+6NRa02xWbHEHFuxEREREN0gwF0NzcjV0x1ZCZsyV4g51axh7z4IxajZEbVs3Ztg49GZbZWH++26hOXqLy/PlMgE9/H+fFBpYvpxit7ZeUMhYnN8IFuxEREREtSTTX4X22ApoktdCZtVLcbt3Bxhj5sMYPhVQ6tyYYcMxWe04lVvmtFPopUJjjdd0aaNFRPvKtpaeAV7QtIBJoY2NBTsRERHRdcgLTkObtAyaU19BcFSuAW5rGwZD7EKYe4wH5M13QqTNIeLsH1ZsOZtvgL2GJVuCfNTXtLV4o1egD7zVLC0bA7+qRERERC4oso6UL814fqdT3BLcH8a4eFg6D212K744RBEZhUakZP8+KTSrFOk5ephtDpfXtNYqpbaWiiK9rZeqCbNu2ViwExEREV1LdEB1cS90iUuhvHqoMgwBlltGwRAbD1tQnBsTrD1RFJGjtzgtp5iaXQq92fWKLTqlHGGBzjuFtvdVc1KoG7FgJyIiIgIAuwXq019Dl7QMioJ0KSzKVDCFTYIxZgHsbbq7McHrKzJanVdsydYjv8z1pFClXEDPAO/fC/PyIr1LGx3knBTqUViwExERUctmKYM2ZR20x5ZDrr8qhR0qH5giHoYxei4cXkFuTLB6BosdaTmVbS0pWaW4Umxyeb5MALq11Tm1tfQM8IJSLmvCrKkuWLATERFRiyQY8qA98Qm0Jz6FzFwsxe26QBij58IU8TBEta8bM6xktTtwOrdMamtJySrFhQIDapgTig6tNE6TQsPa+UCn4ootzRELdiIiImpRZMUXoTv6ITSpGyDYzVLc1voWGGMXwBQ6CZC7b5dNu0PExUIDkq9WtrWcztXDanddnbf1Ujm1tfQK9EFrbfNdtYacsWAnIiKiFkGRewLaxASoz26DIFauiGINjIUhLh6WbqMAoWnbQ0RRRGaJyamtJS1bD4PV9aRQH7UCvQKddwpt563ipNCbGAt2IiIiunmJIpSX90OXuBSqyz85HTJ3GQZjXDys7fs12dKM+WWWayaElvefFxmtLs9XK2QIbecttbWEB/qgUxstZCzOWxQW7ERERHTzcdigPvsNtEkJUOaekMKiIIe55wQY4hbC3rZXo6agN9uQ+ntRXtF3nl1qdnm+XAC6+3s5Lad4i78XFFyxpcVjwU5EREQ3D5sRmrRN0CV9CHnJRSksKrQwhk+FMXo+HL4dG/xpTVY7TuU67xR6sdBY4zVd2mid2lpCArygUXJSKFXFgp2IiIiahFqtgJeXGjKZAIdDRFmZGWazrUEeWzAVQXtyFbTHP4bMmCfFHRo/GKNmw9h7FkRNmwZ5LptDxLm8Mqe2ljN5ZbDXsGRLoI+6vDAPrJwU6qNhGUa1w+8UIiIianRqtQI+PhppYqRcLsDHRwPAVK+iXVaaCe2xj6BNXgvBZpDidp9OMMQ+ClPYA4BSW+fHF0URGUUmp51C03P0MNscLq9ppVE4jZyHB/nA30tV5xyIWLATERFRo/Pyqrq1vSAI8PJS16lgl+enQ3d0GdSnNkNwVF5vaxsOQ1w8zD3GArIbL3NySs1Sv3lKVilSs/UorSE/rVKGsMDKnvPwIG8E+2q4Ygs1KBbsRERE1OhkLiZOuoq7osg8BF3SUqgv7HGKWzrcDkNcPKyd7qz1ii/FRuvvLS2Vu4XmlVlcP7dMQEg7b6mtJTzIB139dJBzUig1MhbsRERE1OgcDhFyedXC1lHTVp0VRAdU53dDl5QAZdbhyjAEWLqPgSF2IWyBMTU+hNFqR1q23qm15UqxyeX5AoBubXVObS09/b2gUjTtOu1EAAt2IiIiagJlZWanHnagvD+8rMz1MoewW6A+tRm6pGVQFJ6uvE6uhilsMowx82FvfUuVy6x2B87klV2zU2gpzucbUNN7g+BWGoQH+iCifXlbS1g7H+hUXLGFPAMLdiIiImp05X3qJmg0CiiVClitNphMtmr71wVLKTTJ66A99hHkZVlS3KHyhSlyBgxRcyB6tSuPiSIuFBic2lpO5ephtbuuzv10ymt6zn0QHuiD1jplg79moobCgp2IiIiahNlcfYFeQTDkQndsJTQnV0FmKZHidq8gGKMfgTF8KjJNKqRcLkVK1jmkZJciLVuPMovd5WN6q+XoFXhtce6NQJ+qE2CJPBkLdiIiInIredE5aJM+hCb9Cwj2yhYZc6vuONl5BnbJ7sSJ8yak/pKMQqPV5eOoFTKEBHgjPMi7vLUl0Aed2mghY3FOzRwLdiIiInILRc4x6BKXQnX2GwiobGE5rQrHcvt4fJEdCTFbBuBqlWvlAnCLv5dTa0v3tjoo5JwUSjcfFuxERETUZPb/9D3yrpxGO8MpTDCudzq2xx6LZbZxOGwKq3Jd5zZap7aW0Hbe0Cg5KZRaBhbsRERE1GhsDhHn88uQerUQqtPboc++iBwEQCmW951bRTm+dtyO5bZ7cErsBABo561y3ik00Ac+GpYs1HIJoijWYgHUlstqtaOoyHD9E4mIiFo4URRxucjktFPoxZx8jBf34hH5N+gky8WnmIyLQid0FK9AtJVgo2Ic/IK6ITyocrdQfy+Vu18KUZMLCPBxeYxvV4mIiKhOckrN0jrnKVmlSM3Wo8RUvgpMa5Rihnw3Zip2oq1QKl1jE8pLj4LWURhxzxSMb6Xhii1E18GCnYiIiK6r2GhFanblWufJWaXIK7NUOa8DcjFP8Q0ekP8AnVC54ovJqxPMcQsgnlUDVzPR2kuLjq21TfgKiJovFuxERETkxGi1Iz1bj5Ts0vLdQrNLcbnI5PJ8AcDQNjlYIN+G2/TfQ4bKddGtAb1hjI2HufsYQCYHzm1sgldAdHNhwU5ERNSCWe0OnMkrk3rOU7L0OJdfBkcNM9yCW2kQHuiD8EAv3K5MR++Mz6DL+MHpHEvHO2CIi4e14yCALS9E9cKCnYiIqIVwiCIuFhgri/PsUpzK0cNid12d++mUf1ixxRtttAqozn0LXVIClNlJ0rmiIIO5+1gY4xbCFtC7KV4SUYvQbAr2jIwMvPHGGzh06BAAYMiQIXjxxRfh5+dX43U//fQTEhISkJycDJlMhujoaDz11FOIiYlpgqyJiIjcQxRFZFVMCs2qnBRaZrG7vMZLJUev35dRjAjyRniQDwJ91JWTQu1maNK/gDZpGRRF5yqfS66GqdcDMMTMh6NV10Z+ZUQtT7NY1rGwsBCTJk2CxWLBjBkzYLfbsXLlSnTo0AGbNm2CSlX98k+HDh3CjBkz0LNnT0yaNAk2mw3r1q1DTk4O1q1bh6ioqOs+N5d1JCKi5qDQYEFKlh7JWSXSxNBCo9Xl+Sq5gNB23k6j553baCGrpn1FMJdAk7wG2mMrITdkS3GHuhWMvWfB2Hs2RJ1/rfJMS0tGSUkxfH1bISws4sZfKNFNqqZlHZtFwf7OO+/go48+wtatW9G9e3cAwIEDBzB79mz885//xJQpU6q9buLEiSguLsY333wDrbZ8JnpeXh7GjBmDiIgIfPLJJ9d9bhbsRETkafRmG9Ky9U5LKl4tMbs8Xy4At/h7SYV5RKAPuvvroJDLanweWVk2tMdWQJO8BjJL5dKMdu/2MEbPhzH8IUDl1WCvi6gla/brsG/fvh19+/aVinUAGDhwILp164bt27dXW7AXFxcjLS0Ns2fPlop1APD390efPn3w888/N0nuRERE9WG2OXA6V+80KfRCgQE1jbZ1bqNFr0BvaTOi0Hbe0CjltX5OeeFZaI8ugybtSwiOyqUbbX6hMMQuhLnnBECurMerIqIb4fEFe3FxMTIyMjBq1KgqxyIiIvDDDz9Ue523tze+/fZbp2K9QmFhIeTy2v/iIiIiago2h4gL+QantpYzeWWw1bBkSztvlVNbS69Ab/hq6lZMK7ISoUtaCtW5nRCueUtgbd8Xhrh4WLoMA4SaR+WJqOF5fMGenV3eKxcYGFjlWEBAAPR6PUpLS+Hj4/wxglwuR9euXatck5aWhsTERAwaNKhR8iUiIqoNURRxucjk1NaSlq2HyeZweU0rjaJ8Uug1E0P9vdX1TQSqi3uhTUqAKvNXp0PmbqNgiF0IW/vb6vcc12APO9GN8/iCvaysDACqHSlXq8t/SRkMhioFu6vHeuGFFwAA8+fPb8AsiYiIaparNzu1taRml6LYZHN5vkYhQ69Ab/T6va0lPMgHHVppKldsqS+7Feoz/4MuKQGK/DQpLMqUMIXcB2PsAtj9ejbMc10jLS0ZmZmXERzckQU7US15fMHucLgeaaggk13/4zmj0YiFCxciLS0Njz76KPr27dsQ6REREVVRYrIiNct5p9BcvcXl+QqZgJ4BXtLIeXh7H3Tz00Eua4QNh6wGaFM+h/bocsj1V6SwQ+kNU8Q0GKPnweHdvuGfl4jqzOMLdi+v8tnnZnPV2e8VsYpzXCkpKcGjjz6KxMRETJo0CU8//XTDJ0pERC2SyWpHeo4eydesd55RZHJ5vgCgq58O4e0r21p6BHhDrWjc3nDBWADt8Y+hPfEpZOYiKe7QBsAQPRemyOkQ1a0aNQciqhuXBfuiRYtu+MEEQcCrr75ar4T+KDg4GACQm5tb5VhOTg58fX2h0+lcXp+fn4+5c+ciNTUVDzzwAP7+97833MeJRETUotjsDpzJK5PaWlKyS3Eurww1bBSKYF+106TQ0Hbe8FY33XiZrOQSdEeXQ5O6HoKt8o2ErVVXGGMXwBR6P6DQNFk+RHTjXP7G2Lx5c5WYIAioWLa9bdu2cDgcKCwsBFDeY96qVcO/M/f19UXHjh2RnJxc5VhKSgoiIyNdXqvX66VifdasWXV6E0JERC2TQxRxqcAoTQhNySpFeo4elhqqcz+dsrKtJcgH4UHeaKOrfnO/xibPTYYuaSnUZ7ZBECt3N7W2iy5f8aXbaEDGFdOImgOXBXtaWprTv0+dOoVp06Zh5syZmD59ulSc6/V6rF27FsuWLcPrr7/eKEmOHDkSq1atwtmzZ502Tjp//jzmzp3r8rp//OMfSE1NxYwZM1isExGRS6IoIru0fFJoRWtLarYeZRa7y2u8VHKntc7Dg3wQ6KN276e4ogjllQPQJS6FKmOf0yFL58EwxMbD2mEgwE+aiZqVWu90OmXKFHTv3h2vvfZatcf/+te/4tixY9iyZUtD5gcAKCgowNixYyGXyzFnzhyYzWasWLECnTt3xvr166FSqZCRkYHExETExcWhU6dOOHv2LMaMGQMfHx/8+c9/rnbd9QkTJlz3ubnTKRHRzafQYJHWOa8YQS8wWF2er5ILCG3nXdnaEuiDzn5ayDyl8HXYoTq3A7qkBChzjklhUZDD3GMcDLELYQ/wjBVZtmzZKK0SM3Fi9TuVE7VEDbLTaXp6OiZOnOjyeEhISKMU6wDg5+eHNWvW4LXXXsP7778PjUaDESNG4Pnnn4dKVf5R42+//YZFixbhtddeQ6dOnXDo0CEAQGlpqcvR9doU7ERE1LyVWWxIy752p9BSZJZUXciggkwAuvt7/d7WUl6kd/f3glLugRsG2UzQpH8BbdKHUBSfl8KiQgNTrwdhiJkPh29nNyZIRA2h1iPso0ePRpcuXfDhhx9WOWa32zFz5kyUlpbi66+/bvAk3Ykj7EREzYfF5sDpXD2Ss/RI+X230AsFBtT0h65Ta400ch7x+6RQjdKze7sFczE0J1dDd2wlZMbKRRkc6tYwRs2GsfcsiNq2bszQtf37v0deXi78/QMwaNBQd6dD5DEaZIR9+vTp+Oc//4knnngCkydPRufOnWE2m3HhwgV89tlnSExMxJtvvtkgCRMREV2P3SHifL7Bqa3ldG4ZbA7X5XmAt0rqNw8P9EGvIG/4apRNmHX9yPRXoT22AprkNZBZy6S43bsDjDHzYQyfCihdr5zmCVikE924Wo+wA8DixYvx0UcfwWqt7PMTRRE6nQ7PPfccpk6d2ihJuhNH2ImI3E8URVwpNkmTQlOzSpGWo4fR6npzPV+NwqmtJTzIBwHe6ibMuuHIC05Dm7QMmlNfQXBU/g22tQ2DITYe5h7jAHnzeeNBRFXVNMJ+QwU7ABQWFuLXX39FZmYmBEFAp06dMHDgwOtuXtRcsWAnImp6eXpzeVvL7yPnqVmlKDbZXJ6vUcgQ9vuKLeGBPoho74MOrTTNft8NxdXD0CUuhfrCLqe4Jbg/jHHxsHQeyhVfiG4SDVqwtzQs2ImIGlepyea01nlKVily9BaX58tlAnr6eyGifeV6513b6qCQ3SSFq+iA6uJe6BKXQnn1UGUYAiy3jIYhdiFsQXFuTLB+8vJyYDaboVar4e/fzt3pEHmMBulhB4AzZ85g27ZtyMvLg91edW3axtjplIiIbh4mqx3pOXqntc4vFRpdni8A6Oqnc2pr6RngDbXCA1dsqS+7BerTX0OXtAyKgnQpLMpUMIVNgjFmAextursxwYaxf/8PXNaR6AbVumD/9ttv8cwzz8DhcN0vyIKdiIgq2OwOnM0zIPma0fNzeWWoYaNQtPdVO+0UGhboDW/1DY0tNT+WMmhT1kF7bDnk+qtS2KHygSlyOoxRc+HwCnRjgkTkbrX+LfjBBx8gODgYb7/9NsLCwqT1z4mIiByiiEuFRqe2llO5ZTDbXA/ytNEqfx8190ZEkC96BXnDT9dy/rYIhjxoj38M7cnPIDMXS3G7LhDG6LkwRTwMUe3rxgyJyFPUumC/cOECXnjhBURFRTVmPkRE5OFEUUR2qfn3FVvKJ4amZpWizFK1VbKCl0qOsEDvyiUVg3wQ5KNu9pNC60JWfAG6o8uhSd0AwV65gZOtdXcYYxfAFHofIG+eq9kQUeOodcEeFBQEk8nUmLkQEZEHKjJYndpaUrJKUWCwujxfJRcQ0s5bamsJD/JBFz8tZC2wOL+WIvcEtIlLoT67HYJY+cmDNTAWhrh4WLqNAoSbsDefiOqt1gX7tGnT8Nlnn+G+++6Dn59fY+ZERERuUmaxIS1b71ScZ5aYXZ4vE4Bb2npJk0IjgnzQ3d8LSjkLTwCAKEJ5+SfoEhOguvyT0yFzl2EwxsXD2r4fl2YkohrVumC3Wq0QBAEjRozAbbfdBj8/vyofZXLSKRFR82GxOXA6rwzJVyt3Cr2Qb0BNa/12bK1x2ik0NNAbWqW8yXJuNhw2qM9+A21SApS5J6SwKMhh7jkBhriFsLft5cYEiag5qXXB/tZbb0n//8cff6z2HBbsRESeye4Qcb7A4DRyfjq3DDaH6/I8wFt1TVuLN3oF+qCVlrtp1shmhCZtE3RJH0JeclEKiwotjOEPwRj9CBy+Hd2YIBE1R7Uu2NPS0hozDyIiaiCiKOJKsen3wrx8UmhadimMVtcrtvioFZVrnf9epLfz4cTH2hJMhdCeXAXt8Y8hM+ZLcYfGD8ao2TD2ngVR08aNGRJRc9agO53a7XbI5TfXR6Pc6ZSIPF1emUUaNU/OKl+xpdhkc3m+WiFDWLvKnvPwIB90bK1pkSu21JesNBPaYx9Bm7wWgq3yb4XdpxMMsY/CFPYAoNS6MUMiai4abKfTw4cPY+fOnTAYDE4bKNntdpSVleHIkSP49ddf654pERHVqNRkk/rNK/7L0Vtcni+XCejp7yW1tYQH+aBbWy8oZCzO60Oenw5dUgLUp7dAcFS+ObK1DYchLh7mHmMB2U2+4RMRNZla/zb58ssv8dJLL6FiQF4QBFw7OK9SqTBkyJAGT5CIqKUyWe1Iz9Ej5ZpVWy4VGmu8pquf1qmtJaSdN9QKrtjSIEQRyquHoE1KgPrCHqdDlg63wxAXD2unO7niCxE1uFoX7KtWrUKnTp3w4Ycfwmq1Yvz48di3bx8EQcBnn32GTz75BA899FBj5kpEdNOy2R04m29wam05l1cGew1Ni0E+ammd84pJod5qjuo2ONEB1fnd0CUthTLrSGUYAizdx8AQuxC2wBj35UdEN70b2un08ccfR7du3QAAXl5e+O233zB27Fg899xzOHXqFJYvX47+/fs3WrJERDcDhyjiUqHxmrYWPU7l6mG2uZ4U2lqr/L3fvLytpVegD9p6qZow6xbIboH61GbokpZBUXhaCotyNUxhk2GIeRSO1t3cmGDztGXLRmRmXkZwcEdMnDjF3ekQNQu1LthlMhnatKmc4d6lSxekpaVh7NixAIChQ4di6dKlDZ8hEVEzJooiskvNTm0tqdml0JvtLq/xUskRFui8U2h7XzUnhTYRwVIKTfJaaI99BHlZthR3qFvBGDkDxt6zIXq1c2OGRNTS1Lpg79y5M06dOiX9u2vXrkhNTZX+bbPZUFZW1rDZERE1M0VGq9OE0OSsUhQYrC7PV8oFhAR4S20tEUG+6OKnhYzFeZMTynKgO/4xNCdXQWYpkeJ2ryAYox+BKWIaRJW3GzMkopaq1gX7mDFjsHjxYmi1WsTHx+P222/HX//6V2zevBndu3fH2rVrpXYZIqKWwGCxIy2ntHyn0N/XO88sNrk8XyYA3drqEB7og4j25SPnPfy9oJRzUqg7yYvOQZv0ITTpX0Cwm6W4rU0PGGIXwhxyLyBn+xERuU+tC/a5c+fi9OnT+Oijj7BgwQJMmDABW7ZswaJFiyAIAuRyORYvXtyYuRIRuY3F5sDpvDKn0fMLBQbUsFEoOrbWOLW1hLbzhk51c+1V0Zwpso9Cl7QUqrM7IKDyRlqDboMhLh6WriMAgW+miMj9bnjjpLy8PPj7+wMArFYrtm/fjqKiItx+++3o2bNnoyTpTtw4iajlsTtEXCioXLElJVuP07l6WGtYssXfS+W01nmvQB+01iqbMGuqFVGEMmMfdIlLobpywOmQuetd5Su+BPd1U3ItAyedElWvQTZOWrJkCUaOHImQkBApplQqMXHiRADA8ePH8de//hV///vf654pEVETE0URmSWm8paW33vO07P1MFhdTwr1USvQK9B5p9B2PuomzJpumMMG9Zlt0CYlQJmXLIVFmQLmkHthiFkAe9tQNyZIROTaDRXsXbt2dSrYr5WYmIivvvqKBTsRebT8MotUmJev2KJHkdH1pFC1QobQdt5SYR4e5INOrTVcsaW5sBqhSV0P3dHlkJdmSGGH0gum8GkwRs+DwyfYjQkSEV2fy4I9IyMDc+fOhd1eOcr06quv4p133qlyriiKyMnJQdeuXRslSSKiutCbbU5tLSlZpcguNbs8Xy4T0MPfq7yt5ffe81v8vaCQsThvbgRTIbTHP4H2xCeQmQqluEPbFsaouTBGzoCoae2+BImIboDLgr1Tp06YOHEifvnlFwDAlStX0Lp1a7Rt27bKuXK5HDExMZg3b17jZUpEVAOT1Y5TuWVOyyleKjTWeE2XNtprdgr1QUiAFzRKTgptzmQll6E9thzalM8h2Crvv923Cwyxj8IUNhlQaN2YIQ0aNARmsxlqNdvIiGqr1pNOhw0bhr/85S8YPnx4Y+fkUTjplMjz2BwizuWVSW0tKVmlOJtvgL2GJVsCfdTXtLV4o1egD7zVte4KJA8nz0uBLikB6tP/gyBWfjJsDegNY2w8zN3HADK+GSMiz1XTpNMbXiWmpWHBTuReDlFERqERKdml0sTQ9Bw9zDaHy2taa5VObS3hQT5o68V1tG86oghl5q/QJi6F+tL3Tocsne6EIXYhrB0HAZxvQETNQL1WiSkuLkZKSgoGDBggxUwmE9asWYMjR45Aq9VixIgRGDNmTMNkS0QtliiKyNFbnNY6T8kuhd7sesUWnVKOsEDva1pbvBHsy0mhNzXRAdW5b6FLSoAyO6kyLMhg7j4WxriFsAX0dmOCREQNq8YR9rVr1+LNN9+ExWJBcnL5MlgWiwXTpk3DyZMnIZfL4ePjg6KiItx11114//33myzxpsIRdqLGU2S0IjX7957zq+UTQ/PLLC7PV8oF9AzwRnigt7RTaJc2Osg5KbRlsJuhSf8S2qRlUBSdk8KiXA1TrwdhiJkPR6subkyQamP//u+Rl5cLf/8ADBo01N3pEHmMOo2wHzhwAP/85z8RHh6OqVOnQhRFCIKAjz/+GCdOnEC3bt3w6aefIjAwELt27cLTTz+NDRs24IEHHmiUF0FEzZvBYkdaTmVbS0pWKa4Um1yeLxOAbm11Tm0tPfy9oFJw58mWRjCXQJO8GtpjKyE35Ehxh7oVjL1nwdh7NkSdvxszpBuRl5eLzMzL7k6DqFlxWbCvXr0a3bt3x8aNG6FQVJ72xRdfQBAEPPXUUwgMDAQAjBw5EqNHj8ZXX33Fgp2IYLU7cPqaFVtSsktxPt+AGuaEokMrjVNbS1g7H+hUnCTYksnKsqE9tgKa5DWQWUqluN07GMaY+TD2mgqovNyYIRFR03BZsB8/fhwzZsxwKtbPnz+Py5cvQ61WY+hQ54+x+vTpgx9++KHREiUiz2R3iLhYaPi9OC8fPT+Vq4fV7ro6b+ulQvi1feeBPmitUzZh1uTJ5IVnoU1KgCb9KwiOyhYpm18oDLELYe45AZDz+4WIWg6XBXtJSUmVNdd//fVXAMCtt94Klcp5xQWZTAaHw/WqDUTU/ImiiKslZqflFNOy9TBYXU8K9VbL0SvQx2mn0HbeKk4KpSoUWYnQJS2F6txOCKh8w2dp3w/GuHhYugwFBLZEEVHL47Jg9/PzQ3Z2tlPsxx9/hCAIuOOOO6qcn5aWBn9/9hAS3UzyyyxObS0pWXoUGa0uz1crZAht5y21tYQH+qBTGy1kLM7JFVGE6uJeaJMSoMr81emQudsoGOLiYQu61U3JERF5BpcF+4ABA7BlyxbMmjULXl5euHDhAn766SfIZDKMHj3a6dzs7Gxs3boVd911V6MnTESNQ2+2/b5ii17aKTS71OzyfLkAdPf3ctoptHtbHRRyjoBSLditUJ/5H3RJCVDkp0lhUaaEKfQ+GGMXwt6mhxsTJCLyHC4L9oULF2LixIkYO3YsoqKicPDgQdhsNkyfPh3t27cHAOTm5uLnn3/GkiVLYDKZMGfOnCZLnIjqzmxz4FSO3qm15WKhscZrOrfRIjyosrUlJMALGiUnhdINshqgTVkH7dGPINdfkcIOpTdMEdNgjJ4Hh3d7NyZIROR5XBbsXbp0wbp16/Df//4XP/74I7y8vLBgwQL86U9/ks5ZuXIlPv30U3h5eeHtt99Gjx4cDSHyNDaHiHN5ZU5tLWfyymCvYcmWQB/175NBy9tbegX6wEdz3X3WiFwSjPnQHv8E2hOfQmYukuIObQAM0XNhipwOUd3KfQkSEXmwGjdOup7ExERkZGRg2LBh8PFxvdh7c8aNk6g5EUURGUUmp51C03L0MNtcTwhvpVE4tbWEB/nA30vl8nyiGyEruQTd0Q+hSd0AwVa57r6tVVcYYxfAFHo/oNC4MUNqamlpySgpKYavbyuEhUW4Ox0ij1HTxkn1KthbAhbs5MlySs3SyHny1VKkZutRara5PF+rlCEs0Of3zYjKdwsN9tVwxRZqcPLcZOiSlkJ9ZhsEsXIVIWu7aBji4mHpNhqQsaWKiKhCnXY6JSLPUmy0IjW79Pe+8/KJoXllFpfnK2QCegZ4OS2n2NVPB7mMxTk1ElGE8soB6BKXQpWxz+mQpfNgGGLjYe0wEOAbRCKiG9JsCvaMjAy88cYbOHToEABgyJAhePHFF+Hn51frx3jppZdw8eJFrF69urHSJGoQRqsdadl6pyUVLxeZXJ4vAOjWVufU1tLT3wsqBVdsoSbgsEN1bgd0SQlQ5hyTwqIgh7nHOBji4mH3D3djgkREzVuzKNgLCwsxc+ZMWCwWzJs3D3a7HStXrkR6ejo2bdpUZROn6mzatAmbNm1C3759myBjotqz2h04UzEp9PfR83P5ZahhTiiCW2mktpbwIB+EBXrDS9UsfpzpZmIzQZP2BbRHl0FRfEEKiwoNjL2mwhgzHw7fTu7LjzwSe9iJbpzLv/BFRUVo3bp1E6bi2qeffoqsrCxs3boV3bt3BwBER0dj9uzZ2LJlC6ZMmeLyWrvdjoSEBCxZsqSp0iVyySGKuFhgRHJWidTWcjpXD4vddXXup1NKo+YRQeX956113Jad3EcwF0NzcjV0x1ZCZsyV4g51axijZsPYezZEbe0//aSWJS0tGZmZlxEc3JEFO1EtuSzY7777bjzzzDOYPHkyAGDJkiUYOXIkQkJCmiy5Ctu3b0ffvn2lYh0ABg4ciG7dumH79u0uC3az2YzJkycjPT0dEydOxC+//NJUKRNBFEVcLTE7tbWkZetRZrG7vMZbLUevwGtWbAn0RqCPmpNCySPI9FehPbYCmuQ1kFnLpLjduwOMMfNhDJ8KKHVuzJCI6ObksmA3GAzIz8+X/r1kyRJ07dq1yQv24uJiZGRkYNSoUVWORURE4IcffnB5rdlshl6vxzvvvIMxY8Zg2LBhjZgptXQFBotTW0tKVikKjVaX56sVMoQEeEttLeFBPujcRgsZi3PyMPKC09AlJUB9ajMER+X3tK1tGAyx8TD3GAfI+akPEVFjcVmwR0VFYfHixdi9eze8vLwAAAkJCdi4caPLBxMEAZ999lmDJpidnQ0ACAwMrHIsICAAer0epaWl1a4D7+3tjV27dkGhYG8vNSy92SZNCq3YLTSr1OzyfLkA3OLv5bRTaPe2OijknBRKnktx9TB0iUuhvrDLKW7pMADG2IWwdB7KFV+IiJqAy0r2jTfewOuvv4709HQUFhZCEAQUFBTAaKx5+/KGVlZW/rGrVqutckytVgMo/zSguoJdJpNBJmNBRPVjtjlwKkd/zU6hpbhYYERNGxh0bqN1amsJbecNjZJrTlMzIDqgurgXusSlUF49VBmGAMsto2GIXQhbUJwbEyQianlcFuzBwcF4//33pX+HhYXhz3/+M8aNG9ckiVVwOFzv0FiBRTk1FJtDxPn8Mqe2ltN5ZbDXsGRLO2+V806hgT7w0fBTHWpm7BaoT38NXdIyKArSpbAoU8EUdj+MMY/C3qZ7DQ9ARESNpdZVxapVq5wmfTaVinYcs7lqu0FFrOIcohshiiIuF5mcRs7TsvUw2Vy/SWylUaDXNYV5RJA3/L3VTZg1UcMSLHpoUtZBe+wjyPVXpbhD5QNT5HQYo+bC4VW1JZGIiJpOrQv2vn37Qq/X46233sKePXuQmZkJpVKJ9u3bY+jQoZg/fz68vb0bPMHg4GAAQG5ubpVjOTk58PX1hU7HVQno+nL1Zqee89RsPUpMNpfna5UyhLXzRniQrzQxtEMrDVdsoZuCYMiD9vjH0J78DDJzsRS36wJhjJ4LU8TDENW+bsyQiIgq1LpgLywsxEMPPYTz58+jW7duuOOOO+BwOHD+/HksX74cO3fuxKZNm+Dr27C/4H19fdGxY0ckJydXOZaSkoLIyMgGfT66ORQbrUjNrmxrSckuRa7e4vJ8hUxAzwAvp9aWbn46yGUszunmIiu+AF3Sh9CkbYRgr/zk0ta6O4yxC2AKvQ+Q81Mjajz+/gFO/0tE11frgv29995DRkYG3n33XYwePdrp2J49e/D0009j8eLF+Mtf/tLgSY4cORKrVq3C2bNnpbacAwcO4Pz585g7d26DPx81L0arHenZeqmtJSWrFBlFJpfnCwC6ttU5tbX0DPCGSsG5EHTzUuSegDZxKdRnt0MQK9u+rIGxMMTFw9JtFCDwZ4Aa36BBQ92dAlGzI4iiWNNiF5I77rgDo0ePdlmQ//vf/8aePXvw/fffN2iCAFBQUICxY8dCLpdjzpw5MJvNWLFiBTp37oz169dDpVIhIyMDiYmJiIuLQ6dO1W+FPWzYMHTo0AGrV6+u9XNbrXYUFRka6qVQPdnsDpzJK7umtUWPc/llqGFOKIJ91U4j52GB3vBScVIotQCiCOXln6BLTIDq8k9Oh8xdhsEYFw9r+35cmpGIyAMEBFRd8bBCrauWoqIi3HLLLS6Pd+vWzWmjpYbk5+eHNWvW4LXXXsP7778PjUaDESNG4Pnnn4dKpQIA/Pbbb1i0aBFee+01lwU7NS8OUcTFAqPTTqGncvSw2F1X53465R9WbPFGG52qCbMm8gAOG9Rnv4E2KQHK3BNSWJQpYO45AYbYBbC37eXGBImI6EbUeoR99OjRCA0NxXvvvVft8SeeeALp6enYuXNngybobhxhbxqiKCKr1HzNTqHlk0LLLHaX13ip5OUrtvze1hIe5INAHzUnhVLLZTNCk7oRuqPLIS+5KIVFhRbG8IdgjJkPh08HNyZIBOTl5cBsNkOtVsPfv5270yHyGA0ywn7//ffjzTffxOuvv4758+fDz88PQHm7yocffojdu3fjySefrH+21CIUGizShNCKVVsKjVaX56vkAkLbeTuNnnduo4WMxTkRBFMhtCdXQXv8Y8iMlZ90OjR+MEbNgbH3TIiaNm7MkKjS/v0/IDPzMoKDO2LixCnuToeoWah1wT5nzhwcP34cn376KT777DNpNZiSkhKIoojhw4fjkUceabREqfkqs9iQlq1H8tXK9c6vllRdV7+CXABu8a9csSUi0Afd/XVQyDkhjuhastJMaI8thzZ5HQRb5SeBdp9OMMQ+ClPYA4Cy6i7RRETUvNS6YJfJZHj//ffx448/Yu/evbhy5QpEUUSHDh0wbNgwDB48uDHzpGbCbHPgdK7+mtYWPS4UGFBT31XnNlr0CiwfPY8I8kFoO29olPImy5mouZHnp0OXlAD16S0QHJV7CVj9I2CMi4e5+z2AjBOriYhuFjf8G/3OO+/EnXfe2Ri5UDNjd4g4n29w2in0dG4ZbDUs2dLOW3XNhFAf9Aryhq9G2YRZEzVTogjl1UPQJiVAfWGP0yFLh9thiIuHtdOdXPGFiOgmxCEYqhVRFHGl2OTUc56WrYfJ5nB5ja9GgfBAn993CfVFRJA3/L25IQvRDREdUJ3fDV3SUiizjlSGBRnMt4yBMW4hbO2i3ZggERE1NhbsVK1cvdmprSU1uxTFJpvL8zUKGXoFeqPX720t4UE+6NBKwxVbiOrKboYmfTO0R5dBUXhGCotyNUxhk2GIeRSO1t3cmCARETUVFuyEEpMVqVnOO4Xm6C0uz1fIBPQM8JLaWsKDfNC1rQ4KGYtzovoSLKXQJK+F9thHkJdlS3GHuhWMkTNgjJoDUcct3YmIWhIW7C2MyWpHeo5eamtJySpFRpHJ5fkCgK5+ut/bWsqL854B3lAruGILUUMSynKgO74SmpOrIbOUSHG7VxCM0Y/AFDENosrbjRkSEZG7sGC/idnsDpzJK5PaWlKyS3Eurww1bBSK9r5qqaUl/PcVW7zV/DYhaizyonPQJn0ITfoXEOyVy53a2vSEIXYhzCETATl36yUiaslqvdMpAFgsFmzYsAE//PADMjMz8eqrr0Kj0WDbtm2YO3eutJnSzaS57HTqEEVcKjA6tbWcyi2DuYZJoX46pVNbS3iQN9roWBgQNQVF9lHokpZCdXYHhGsWPrUG3QZDXDwsXUcAAj/JIiJqKRpkp1O9Xo9Zs2bh5MmT8Pf3R35+PkwmE3Jzc7Fy5Up8++23WLt2LYKCghokaXJNFEVkl5p/X7GlfOQ8NasUZRa7y2u8VHJprfOK9c4DfdScFErUlEQRyox90CUuherKAadD5q53wRAXD1v7Pm5KjoiIPFWtC/b33nsP6enp+PjjjxEWFoaBAwcCAEaOHImlS5fi2WefxXvvvYfXXnut0ZJtqYoMViRnlyLlmp1CCwxWl+er5AJC2nkjPNAHEe3LR9A7+2khY3FO5B4OG9RntkKXmABFfooUFmUKmEPugyHmUdjbhroxQSIi8mS1Lth37tyJhx56CAMHDkRhYaHTsWHDhmHatGnYtm1bgyfYEl0tMWFPeq7U2pJZYnZ5rkwAuvt7XbPeuQ+6+3tBKedH6URuZzVCk7oeuqPLIS/NkMIOpRdM4dNgjJkHh3ewGxMkIqLmoNYFe2FhIbp37+7yeMeOHVFQUNAgSbVkdoeIh1cnosTFmuedWmucdgoNDfSGVilv4iyJqCaCqRDa459Ae+ITyEyVAxwObVsYo+bCGDkDoqa1+xIkcqMtWzYiM/MygoM7YuLEKe5Oh6hZqHXB3rFjR5w4cQJTplT/w3XgwAF06NChwRJrqURRRDtvNUpMNgR4q6QJoRFBPugV5A1fjdLdKRKRC7KSy9AeWw5tyucQbEYpbvftAkPsApjC7gcUWjdmSEREzVGtC/bJkyfjrbfeQlRUFO68804AgCAI0Ov1SEhIwO7du/HUU081Vp4thkIuw9oZcSgx2dBay+KcqDmQ56VAl5QA9en/QRArJ39bA3rDGBsPc/cxgIyfhBERUd3UumCfPXs2Tp8+jZdffhkyWXl/9OOPP46ysjKIoojhw4dj3rx5jZZoSyITBBbrRJ5OFKHM/BXaxKVQX/re6ZCl050wxMbD2vF2gJO9iYionmpdsAuCgNdeew0TJ07Erl27kJGRAbvdjg4dOmD48OEYPHhwY+ZJROQZHHaozu+ELnEplDlHpbAoyGDuMQ7G2AWwBfR2X35ERHTTueEtLPv164d+/fo1Ri5ERJ7LboYm7Qtoj34IRdE5KSzK1TD1ehCGmPlwtOrixgSJiOhmVeuCfcuWLbU6b+LEiXVMhYjI8wjmEmiSV0N7bCXkhhwp7lC3grH3LBh7z4ao83djhkREdLOrdcH+4osvQhAEiKJY5di1u2WyYCeim4GsLAvaYyugObkGMqteitu9g2GMmQ9jr6mAysuNGRIRUUtR64J91apVVWJ2ux15eXnYvn07Ll68iISEhAZNjoioqckLz0KblABN+lcQHBYpbvMLhSF2Icw9JwByTgonIqKmI4jVDZnXwSOPPAJ/f3+89tprDfFwHsNqtaOoyODuNIiokSmyjkCXlADVuZ0QUPlr0dK+H4xx8bB0GcYVX4gaQF5eDsxmM9RqNfz927k7HSKPERDg4/LYDU86deWuu+7CW2+91VAPR0TU+EQRqot7oU1aClXmQadD5m6jYIiLhy3oVjclR3RzYpFOdOMarGA/f/48bDZbQz0cEVHjsVuhPvM1dIkJUBSkS2FRpoQp9D4YYxfC3qaHGxMkIiKqVO9VYiwWC1JTU7Fx40YMHz68ofIiImp4ljJoUz+H9uhHkOuvSGGH0humyIdhjJ4Hh1eQGxMkIiKqqtY97GFhYS5XiQGAiIgILF68GMHBwQ2aoLuxh52o+ROM+dAe/xjaE59BZi6S4g5tAAzRc2GKnA5R3cp9CRK1IPv3f4+8vFz4+wdg0KCh7k6HyGM0SA97davEAIBMJkNAQAC6dOGGIUTkWWQll6A7+iE0qRsg2ExS3NaqG4yxC2AKnQQoNG7MkKjlycvLRWbmZXenQdSs1Lpg79u3b2PmQUTUYOS5ydAlLYX6zDYIol2KW9tFwxAXD0u30YBM7sYMiYiIas9lwV7bnU3/iBsnEZFbiCKUVw5Al7gUqox9TocsnYfAEBcPa/AALs1IRETNjsuCvaadTV0RBIEFOxE1LYcdqnM7oEtKgDLnmBQWBTnMPcbBEBcPu3+4GxMkIiKqH5cFu6uedSIij2AzQZP2BbRHl0FRfEEKiwoNjL2mwhgzHw7fTu7Lj4iIqIG4LNjZs05EnkgwFUF7cjW0x1dCZsyT4g51axijZsPYezZErZ8bMyQiImpYN7RxksFgwIEDB1BWVubUKmOz2VBWVoZff/0VCQkJDZ4kEZFMfxXaYyugSV4DmbVMitt9OsIQMx+mXg8CSp0bMyQiImoctS7YExMT8eijj0Kv10uxiqJd+H0SV5s2bRo4PSJq6eQFp6FLSoD61GYIDqsUt7XtBUPsQph7jAPkSjdmSERE1LhqXbC/9957cDgc+Nvf/gZRFPG3v/0NH3zwAQwGA9avX4/k5GSsX7++MXMlohZEcfUwdIlLob6wyylu6TAAhth4WDsP4YovRM1QWFgEgoM7wteXm5UR1Vatdzq97bbbMHXqVDz77LOwWq2IjY3F4sWLMXToUFgsFkyaNAkhISF46623GjvnJsWdTomakOiA6sJ30CUthfLqb5VhCLB0vxuG2IWwBca6MUEiIqLG0SA7nZrNZnTt2hUAoFQq0aVLF6SmpmLo0KFQqVSYOHEi1q5dW+9kiagFslugPv01dIkJUBSeksKiTAVT2P0wxi6AvfUtbkyQiIjIfWpdsAcGBiI7O1v6d6dOnZCeni7928fHB/n5+Q2bHRHd1ASLHpqUddAe+why/VUp7lD5wBQ5HcaouXB4BboxQyIiIverdcF+5513Yu3atYiJicHAgQMRGxuLjz76CBkZGWjfvj2+/fZbBAbyDysRXZ9gyIX2+CfQnvwMMnOxFLfrAmGMngdT5MMQVa4/GiSi5istLRklJcXw9W2FsLAId6dD1CzUuoc9Ly8PDz30EDIyMnDgwAEIgoCxY8eisLAQWq0WZWVlePrppzF//vzGzrlJsYedqOHIii9Al/QhNGkbIdjNUtzWujuMsQtgCr0PkKvdmCERNbYtWzYiM/MygoM7YuLEKe5Oh8hj1KmH/YknnsCECRMwePBgKBQK+Pv7Y+vWrfjuu++k5Rs3btyIxYsXo7i4GHfeeScefPDBhs/+dxkZGXjjjTdw6NAhAMCQIUPw4osvws+v5g1S6nodETUcRc5xaJMSoD67HYLokOLWwDgY4uJh6TYSEGRuzJCIiMhzuRxhj4yMhN1uh6+vL+6++26MHz8ecXFxTZ0fAKCwsBCTJk2CxWLBjBkzYLfbsXLlSnTo0AGbNm2CSqVq0OuuxRF2ojoSRSgv/wRd4lKoLu93OmTuMhzGuIWwtu/HpRmJWhiOsBNVr04j7AcOHMDOnTvxzTffYNOmTdiwYQOCg4Mxfvx4jBs3Drfc0nQrNnz66afIysrC1q1b0b17dwBAdHQ0Zs+ejS1btmDKlOp/4Ot6HRHVg8MG9dnt0CYmQJl3UgqLMgXMPSfAELsA9ra93JggERFR81KrHvaCggLs2LEDO3bswJEjRwAA4eHhmDBhAu655x60bdu2UZMcMWIEOnbsiE8//dQpPnr0aAQGBuKzzz5r0OuuxRF2olqyGaFJ3Qjd0eWQl1yUwqJCC2P4QzDGzIfDp4MbEyQiT8ARdqLq1Xsddj8/P0ybNg3Tpk1DdnY2duzYgW+++Qavvvoq/vOf/6B///4YP3487rrrLmi12gZLHACKi4uRkZGBUaNGVTkWERGBH374oUGvI6IbI5gKoT3xGbQnPoHMWLm0q0PjB2PUHBh7z4SoaePGDImIiJq3Wi/rWCEwMBCzZs3CrFmzcPnyZezatQvff/89/vznP+Nvf/sbEhMTGzTBirXfq1syMiAgAHq9HqWlpfDx8WmQ64iodmSlmdAeWw5t8joItspPoey+nWGImQ9T2AOAsmHfwBMREbVEN1ywX8vX1xd+fn5o27Yt1Go1jEZjQ+UlKSsrA4BqR+7V6vLl3wwGQ5XCu67XEVHN5Plp0CUtg/r0FggOmxS3+kfAGBcPc/d7AFm9frUQERHRNW74r2pRURH27NmDb7/9Fr/++ivsdjtCQkIQHx+PsWPHNniCDofjuufIZFWXg6vrdURUDVGE8uohaBOXQn3xO6dDlo6DYIiLh7XjHVzxhYiuy98/wOl/iej6alWwFxQUSEX6oUOHYLPZEBwcjNmzZ2P8+PHo2bNnoyXo5eUFADCbzVWOVcQqzmmI64joGqIDqvO7oEtKgDLrSGVYkMF8yxgY4xbC1i7ajQkSUXMzaNBQd6dA1Oy4LNjz8/Oxa9cu7Ny5E4cPH4bNZkOrVq1w3333Ydy4cejTp0+TJBgcHAwAyM3NrXIsJycHvr6+0Ol0DXYdEQGwm6FJ3wzt0WVQFJ6RwqJcDVPYFBhi5sPRupsbEyQiImo5XBbsd9xxB0RRhFKpxLBhwzB+/HgMHjwYSqWyKfODr68vOnbsiOTk5CrHUlJSEBkZ2aDXEbVkgqUUmpNroD2+AvKybCnuULeCMXIGjFFzIOr4MTYREVFTclmw33rrrZgwYQJGjx4Nb2/vpsypipEjR2LVqlU4e/astAHSgQMHcP78ecydO7fBryNqaYSyHOiOr4Tm5GrILCVS3O4VBGPMfJjCH4Kocu/vASK6OeTl5cBsNkOtVsPfv5270yFqFmq1cZK7FRQUYOzYsZDL5ZgzZw7MZjNWrFiBzp07Y/369VCpVMjIyEBiYiLi4uLQqVOnWl93Pdw4iW5m8qJz0CZ9CE36FxDslfM9bG16whC7EOaQiYD8+j8nRES1xY2TiKpX742T3M3Pzw9r1qzBa6+9hvfffx8ajQYjRozA888/LxXdv/32GxYtWoTXXntNKthrcx1RS6TIPgpd0lKozu6AgMr37Nb2fWCIjYel63BA4CpKREREnqBZjLC7E0fY6aYhilBe+qG8UL/yi9Mhc9eRMMQthK1900wmJ6KWiyPsRNVr9iPsRFQPDhvUZ7ZCl5gARX6KFBZlCphD7oMhdgHsfiFuTJCIiIhqwoKd6GZlNUCTuh66o8shL70shR1KL5jCp8EYMw8O72A3JkhERES1wYKd6CYjmAqhPf4JtCc+gcxUKMUdWn8Yo+bCGDkdoqa1+xIkIiKiG8KCnegmISu5DO3RD6FNXQ/BZpTidt8uMMQugCnsfkChdWOGREREVBcs2ImaOXleCnRJCVCf/h8E0S7FrQFRMMTFw3LL3YBM7sYMiYiIqD5YsBM1R6IIZeYv0CUuherSD06HLJ0GwxC7ENaOtwOC4J78iIiIqMFwWcfr4LKO5FEcdqjO74QucSmUOUelsCjIYO4xDsbYhbAFRLovPyIiIqoTLutI1NzZTNCkfwnt0Q+hKDonhUW5GqZeD8IQMx+OVl3cmCARERE1FhbsRB5MMJdAk7wa2mMrITfkSHGHuhWMvWfBGDUHoratGzMkIiKixsaCncgDycqyoD22ApqTayCz6qW43TsYxpj5MPaaCqi83JghERERNRUW7EQeRF54FtqkBGjSv4LgsEhxm18oDHELYe4xAZAr3ZghEVH9bNmyEZmZlxEc3BETJ05xdzpEzQILdiIPoMg6Ur7iy/ldEFA5D9wS3A/G2HhYugzjii9EREQtFAt2IncRRagu7oU2aSlUmQedDpm7jYIhLh62oFvdlBwRERF5ChbsRE3NboX6zNfQJSZAUZAuhUWZEqbQ+2CMXQh7mx5uTJCIiIg8CQt2oqZiKYM29XNojy6HXJ8phR1Kb5giH4Yxeh4cXkFuTJCIiIg8EQt2okYmGPOhPf4xtCc+g8xcJMUd2gAYoufCFDkdorqV+xIkIiIij8aCnaiRyIovQndsOTSpGyDYTFLc1qobjLELYAqdBCg0bsyQiIiImgMW7EQNTJF7EtqkBKjPbIUgOqS4tV0MDHHxsHQbBcjkbsyQiIiImhMW7EQNQRShvPwzdEkJUGXsczpk6TwEhrh4WIMHcGlGIiIiumEs2Inqw2GH6twO6BKXQpl7XAqLghzmnuNhiF0Iu3+4GxMkIvIsgwYNgdlshlqtdncqRM0GC3aiurCZoEnbBF3SMshLLkphUaGBsddUGGPmw+HbyY0JEhF5Jn//du5OgajZYcFOdAMEUxG0J1dDe3wlZMY8Ke7QtIGx92wYe8+CqPVzY4ZERER0s2HBTlQLMn0mtMdWQpO8BjJrmRS3+3SEIWY+TL0eBJQ6N2ZIRERENysW7EQ1kBecgi5pGdSnNkNwWKW4rW0vGOLiYe4+FpAr3ZghEVHzsn//98jLy4W/fwAGDRrq7nSImgUW7ETVUFz9DbrEpVBf2O0Ut3QYCEPsQlg7D+GKL0REdZCXl4vMzMvuToOoWWHBTlRBdEB14TvokpZCefW3yjAEWLrfDUPsQtgCY92YIBEREbVELNiJ7BaoT22BLmkZFIWnpLAoU8EUdj+MsQtgb32LGxMkIiKilowFO7VYgkUPTco6aI99BLn+qhR3qHxgipwBY9QcOLwC3ZghEREREQt2aoEEQy60xz+G9uQqyMzFUtyuC4Qxeh5MkQ9DVPm4MUMiIiKiSizYqcWQFV+ALulDaNI2QrCbpbitdXcYYxfCFHovIOfOe0RERORZWLDTTU+RcxzaxKVQn/sGguiQ4tbAOBji4mHpNhIQZG7MkIiIiMg1Fux0cxJFKC//BF3iUqgu73c6ZO4yHMa4eFjb9+XSjEREROTxWLDTzcVhg/rsdmgTE6DMOymFRZkC5p4TYIhdAHvbXm5MkIioZQsLi0BwcEf4+rZydypEzYYgiqLo7iQ8mdVqR1GRwd1p0PXYjNCkboTu6IeQl1ySwqJCB2PEQzBGPwKHTwc3JkhERETkWkCA6wUvOMJOzZpgKoT2xGfQnvgEMmO+FHdo/GCMmgNj75kQNW3cmCERERFR/bBgp2ZJVnoF2mMfQZu8DoKt8hMQu29nGGIehSlsCqDUujFDIiIioobBgp2aFXl+GnRJCVCf/hqCwybFrf6RMMYthLn7PYCM39ZERJ4qLS0ZJSXF8PVthbCwCHenQ9QssLIhzyeKUF49CG1iAtQXv3M6ZOk4CIa4eFg73sEVX4iImoG0tGRkZl5GcHBHFuxEtcSCnTyX6IDq/C7oEpdCmZ1YGRZkMN8yBsa4hbC1i3ZjgkRERESNjwU7eR67GZr0r6BNWgZF0VkpLMrVMIVNgSFmPhytu7kxQSIiIqKmw4KdPIZgKYXm5Bpoj6+AvCxbijvUrWCMnAFj1ByIugA3ZkhERETU9Fiwk9vJyrKhPf4xNCdXQ2YpkeJ2ryAYY+bDFP4QRJW3GzMkIiIich+ZuxOorXXr1mH06NGIiorCuHHjsH379ht+jEuXLiEqKgoHDx5shAzpRsmLzsH7++fht2oAdIkfSMW6rU0ISoa/g4LpB2CMmc9inYiIiFq0ZjHCvnLlSvznP//B6NGjMWvWLOzevRvPPPMMBEHAmDFjavUYJSUleOyxx2A2mxs5W7oeRXYSdEkJUJ3dAQGVG+1a2/eBITYelq7DAaHZvJckIiIialQeX7CXlJRgyZIlGDt2LN566y0AwJQpUzB9+nT85z//wahRoyCXy2t8jLNnz+Lxxx/HuXPnmiJlqo4oQnnpB+iSlkJ15RenQ+auI2GIWwhb+z5uSo6IiIjIc3n8MObevXthMBgwdepUKSaTyfDQQw/h6tWrSEpKqvH6zZs3Y8KECSgqKsLkyZMbO136I4cN6vSv0GbDSLTeNl0q1kWZEqawKSiYuhcl93zMYp2IqIXw9w9AcHBH+PtzEQGi2vL4EfaTJ08CACIinDdXCA8Pl47fdtttLq9PT0/HPffcg+eeew4//vgjNm3a1HjJUiWrAZrU9dAdXQ556WUp7FB6wRTxMIzRc+HwDnZjgkRE5A6DBg11dwpEzY7HF+w5OTlo1aoVtFqtUzwgoPydeWZmZo3XP/PMM1CpVI2WHzkTjAXQnvgU2hOfQGYqlOIOrT+MUXNhjJwOUdPafQkSERERNTNuK9hzc3NrPK7T6eDl5YWysjJoNJoqxytiRqOxxsdhsd40ZCUZ0B5dDm3qegi2ynti9+0CQ+xCmMImAQptDY9ARERERNVxW8E+aNCgGo8vWLAATz/9NBwOBwRBqHK8IlbdMWo68rwU6BKXQn1mKwTRLsWtAVEwxMXDcsvdgKzmScFERNRy5OXlwGw2Q61Ww9+/nbvTIWoW3Faw/+tf/6rxeK9evQAAXl5eMJlMVY5XjKx7eXk1fHJUM1GEMvMX6BKXQnXpB6dDlk6DYYiLh7XDQIBvpoiI6A/27/8BmZmXERzcERMnTnF3OkTNgtsK9tqu2NK+fXsUFxfDYrE4tbfk5OQAAAIDAxslP6qGww7V+W+hS1wKZc4xKSwKMph7jIMxdiFsAZFuTJCIiIjo5uPxk04jIiIgiiJSU1MRHR0txVNTUwEAvXv3dldqLYfNBE36l9Ae/RCKosq17EW5GqbwB2GIng9Hqy5uTJCIiIjo5uXxBfvgwYOhVquxevVqqWB3OBxYt24dOnTogJiYGPcmeBMTzCXQnFwF7fGPITfkSHGHuhWMvWfBGDUHoratGzMkIiIiuvl5fMHepk0bzJ8/H4sXL4Yoiujfvz927tyJw4cP491333Xa5XTPnj0AgBEjRrgr3ZuCrCwL2mMroDm5BjKrXorbvYNhjJkPY6+pgIpzB4iIiIiagscX7ADw2GOPQavVYu3atdi9eze6du2K9957D6NHj3Y679VXXwXAgr2u5IVnoE1KgCb9KwgOqxS3+YXCELcQ5h4TALnSjRkSERERtTyCKIqiu5PwZFarHUVFBnen0agUWUfKV3w5vwsCKr8dLMH9YIyNh6XLMK74QkREDWLLlo1cJYaoGgEBPi6PNYsRdmoEogjVxb3QJi6F6urByjAEWLqNhCEuHragW92YIBEREREBLNhbHrsV6jNfQ5eYAEVBuhQWZUqYQifBGLsA9jY93JggEREREV2LLTHXcdO0xFjKoE39HNqjyyHXZ0phh9IbpsiHYYyeB4dXkBsTJCIiImq52BLTggnGfGiPfwztiU8hMxdLcbuuHYzRc2GKmA5R7evGDImIiIioJizYb1Ky4ovQHV0OTep6CHazFLe16gZj7AKYQicBCo0bMyQiIiKi2mDBfpNR5J6ENikB6jNbIYgOKW5tFwNDXDws3UYBMnkNj0BEREREnoQF+81AFKG8/DN0SUuhyvjR6ZCl8xAY4uJhDR7ApRmJiMjtuKwj0Y1jwd6cOexQn/0G2qQEKHOPS2FRkMPcczwMsQth9w93Y4JEREREVF8s2JsjmxGatC+gS1oGeclFKSwqtDCGT4Ux+hE4fDu5MUEiIiIiaigs2JsRwVQE7cnV0B5fCZkxT4o7NG1g7D0bxt6zIGr93JghERERETU0FuzNgEyfCe3RFdCkrIXMWibF7T6dYIiZD1OvBwClzo0ZEhEREVFjYcHuweQFp6BLWgb1qc0QHFYpbmvbC4a4eJh7jANkvIVERERENzNWex5InpcCr4P/hfrCbqe4pcNAGGIXwtp5CFd8ISKiFikx8TCSko5g7txH3Z0KUZNhwe5p7Fa03jJZ2pVUhABL97thiF0IW2Csm5MjIiJyj23btkCpVCEgoB0AIDn5JPbt+w7x8U+6OTOixseC3dPIFLC36gohPw2m0Ekwxi6AvfUt7s6KiIjIre6+exy+/vorfPXVBygtLYHBUIaZM+e6Oy2iJsGC3dMIAoru3wpAYNsLERHddAYNGgKz2Qy1Wn3D1wrX/F2UyeQQBFlDpkbksQRRFEV3J+HJrFY7iooM7k6DiIioRdu2bQsUCiXatQtEUtIR9O9/O1ti6KYSEODj8hhH2ImIiMjjjR07EUD5pFMAiIiIREREpBszImo6HGG/Do6wExERNZ6yMj2++GID9u7djaysqxBFoF27doiLuw3PPPOCu9MjajIcYSciIiKPsH//98jLy4W/fwD69r0d8fGPIDv7KsaMGYeuXW+ByWTCuXNnkJFxyd2pEnkMFuxERETUZPLycpGZeRkAsH//jzh79jTefnsJ+vbt7+bMiDwXp1cTERGRW5SWlgAAUlJOwuFwuDkbIs/Fgp2IiIjcYujQ4ejRIwQrVizDxIl34403/oWff/6JxTvRH7AlhoiIiNzC17cVVq5cjcOHD+GXX37G/v37sHXrFkRGRmHx4g+hVCrdnSKRR+AIOxEREbmNXC5Hv34D8NRT/4cNG7Zg1Ki7cfLkcZw5c8rdqRF5DBbsRERE1OQMBgP+uLK0XC7/fQdTAQEB7dyUGZHnYUsMERERNbkffvgeW7Zsxh13DEaHDp0gig4cPPgrDhz4CQ89NB3+/gHuTpHIY7BgJyIioibXpUsXtGrVGnv37kFRUSF8fX3RtesteO21N3HHHUPcnR6RR2HBTkRERE0mLCwCwcEd4evbCmFhEe5Oh6hZEMQ/NpCRE6vVjqIig7vTICIiIqKbWECAj8tjnHRKREREROTBWLATEREREXkwFuxERERERB6MBTsRERERkQdjwU5ERERE5MFYsBMREREReTAW7EREREREHowFOxERERGRB2PBTkRERETkwViwExERERF5MBbsREREREQejAU7EREREZEHY8FOREREROTBBFEURXcnQURERERE1eMIOxERERGRB2PBTkRERETkwViwExERERF5MBbsREREREQejAU7EREREZEHY8FOREREROTBWLA3sYyMDDz++OPo27cv+vbti+effx4FBQWNdh25R0Pcr5deegnTp09vpAypIdT1Pv/000946KGHEB0djdjYWMyaNQtHjx5t/ISpTup6n3/55RdMnToVsbGxuOOOO/Dvf/8bZWVlTZAx3aiG+J2dlpaGyMhILF68uJGypPqq632+//77ERoaWuW/J554ogmyLqdosmciFBYWYubMmbBYLJg3bx7sdjtWrlyJ9PR0bNq0CSqVqkGvI/doiPu1adMmbNq0CX379m2CjKku6nqfDx06hEceeQQ9e/bE008/DZvNhnXr1uHhhx/GunXrEBUV1cSvhGpS1/v866+/Ys6cOYiIiMD//d//4erVq1i1ahVOnjyJtWvXQibjeJmnaIjf2TabDYsWLYLVam2CjKku6nqfRVHE2bNnMWLECIwcOdLpWIcOHZoidSkRaiJvv/222KtXL/HMmTNS7OeffxZDQkLEDRs2NPh15B71uV82m01cvHixGBoaKoaEhIgPP/xwY6dLdVTX+zxhwgRxyJAhosFgkGK5ublinz59xFmzZjVqznTj6nqf7733XnHo0KGi0WiUYmvWrBFDQkLEH374oVFzphvTEH9jlyxZIkZERIghISHi+++/31ipUj3U9T5funRJDAkJEb/88sumSNMlvsVvQtu3b0ffvn3RvXt3KTZw4EB069YN27dvb/DryD3qer/MZjPuvfdeLF68GBMmTEBgYGBTpEt1VJf7XFxcjLS0NIwePRparVaK+/v7o0+fPkhKSmr0vOnG1OU+m81mtGnTBlOmTIFGo5HiFZ+YpaenN27SdEPq+zc2PT0dCQkJiI+Pb8w0qZ7qep/PnDkDAE7XuQML9iZSXFyMjIwMREREVDkWERGBkydPNuh15B71uV9msxl6vR7vvPMO3njjDSgU7FjzVHW9z97e3vj2228xa9asKscKCwshl8sbOlWqh7reZ7VajZUrV2LBggVO8dTUVABAcHBwwydLdVLfv7EVrTADBw7E+PHjGytNqqf63OfTp08DqCzYDQZD4yR5HSzYm0h2djYAVDtqGhAQAL1ej9LS0ga7jtyjPvfL29sbu3btwpgxYxo1R6q/ut5nuVyOrl27VrkuLS0NiYmJiI2NbZyEqU4a6vfvlStX8NVXX+Hf//43QkJCcNdddzV4rlQ39b3HH330ES5evIh//OMfjZYj1V997vPp06fh5eWF1157DbGxsYiNjcWIESOavMOBQ3hNpGJlgGs/Bq+gVqsBlL9r8/HxaZDryD3qc79kMhknojUTDflzWVZWhhdeeAEAMH/+/AbMkuqrIe5zUVERhg0bJj3OSy+9JF1L7lefe3z69Gl88MEHeOWVVxAUFITLly83brJUZ/W5z2fOnEFZWRlKS0vxn//8ByUlJVi1ahWeeeYZWK1WTJw4sVFzr8CCvYk4HI7rnlNdsVbX68g9eL9ahoa6z0ajEQsXLkRaWhoeffRRrgrkYRriPguCgHfeeQcWiwWrV6/G7Nmz8fbbb2P06NENlSbVQ13vsd1ux6JFi3DrrbdiypQpjZEaNaD6/CxPmTIFDocD06ZNk2L33HMPxo4di//+978YN25ck7QzsnJoIl5eXgDK+5T/qCJWcU5DXEfuwfvVMjTEfS4pKcGcOXNw8OBBTJo0CU8//XTDJ0r10hD3uVWrVhgzZgwmTpyItWvXIjg4GK+//nrDJ0t1Utd7vHLlSqSlpeHZZ59FQUEBCgoKUFJSAqD8jXhBQUGtikRqGvX5WZ46dapTsQ4AGo0GEyZMQF5enjQptbGxYG8iFZOMcnNzqxzLycmBr68vdDpdg11H7sH71TLU9z7n5+djxowZSExMxAMPPIB///vfEASh0fKlumnon2eNRoMhQ4bg6tWr3PjOQ9T1Hv/000+wWq2YPHkyBgwYgAEDBuDee+8FUF7MDxgwAJmZmY2bPNVaY/xt9vPzA9B0k1DZEtNEfH190bFjRyQnJ1c5lpKSgsjIyAa9jtyD96tlqM991uv1mDt3LlJTUzFr1iwsWrSoMVOleqjrfT579iweeeQRzJ07t8rIXFlZGQRB4IZ3HqKu9/iFF16QRtQr5OXl4bnnnsOECRMwceJEBAQENErOdOPqep+zs7MxZ84c3H333Xj88cedjp0/fx4A0LFjx4ZPuBocYW9CI0eOxC+//IKzZ89KsQMHDuD8+fM1rgxS1+vIPXi/Woa63ud//OMfSE1NxYwZM1isNwN1uc9dunRBaWkp1q9fD4vFIsWvXLmCXbt2oU+fPvD29m703Kl26nKPIyMjMXDgQKf/4uLiAACdOnXCwIEDObnYw9TlPgcGBqK0tBSbNm2CXq+X4levXsVXX32Ffv36NdkbM0EURbFJnolQUFCAsWPHQi6XY86cOTCbzVixYgU6d+6M9evXQ6VSISMjA4mJiYiLi0OnTp1qfR15jrre5z8aNmwYOnTogNWrVzfxK6DaqMt9Pnv2LMaMGQMfHx/8+c9/rnai0oQJE9zwasiVuv48f/3113j++ecRExOD8ePHo7CwEGvXroXVasW6desQEhLi5ldGFRrqd/bly5cxfPhwPP744/jTn/7UxK+Crqeu93nPnj147LHH0LNnT0yePBllZWXSz/Lnn3/edBsquXWf1Rbo7Nmz4rx588SYmBixf//+4gsvvCDm5+dLx7/88stqt8C93nXkWep6n681dOhQ8eGHH26KdKmObvQ+r1u3TgwJCanxP/I8df153r59u3jvvfeKERERYp8+fcQ//elP4rlz55o6faqFhvidnZGRIYaEhIjvv/9+U6RMdVDX+7x7927x/vvvFyMjI8XbbrtNjI+PF8+cOdOkuXOEnYiIiIjIg7GHnYiIiIjIg7FgJyIiIiLyYCzYiYiIiIg8GAt2IiIiIiIPxoKdiIiIiMiDsWAnIiIiIvJgLNiJiIiIiDyY4v/budOQKrc9juNf9WRpg1kWFipFRUlaWFFt0wZ3EWYaYqU4BL0JsQGNChOFBm0wjIjmyEqyUgspU7SiOdAXEZRhg5pRUpCVJebsvi+kTfsa555O5173td/n1bPX818+i/Xqt9f+P/b0AkRE5J+TmJhIfn6+xVifPn0YPnw48+bNY926dTg5OVnUPnv2rCeWKiIif5ECu4hIL7R582acnZ0BaGlpobKykpycHB4/fsy5c+ews7MjPDwcg8HQwysVEZH/RIFdRKQXmj9/Pm5ubhZjo0aNYuvWrdy5c4d58+bh4+ODj49PD61QRET+KvWwi4j8JmbMmAFAZWVlD69ERER+hgK7iMhv4t27dwC4u7sDXT3s48ePt6ipra1l48aNzJw5E29vb0JCQsjNzbWoSUxMZPHixTx48IDw8HAmTZqE0WgkPz+ftrY2MjIymDVrFtOnTyc+Pp5Pnz5ZzC8uLiY6OpqpU6fi5eVFQEAA6enptLa2mmtaW1tJS0vDaDTi5eXFnDlz2Lp1K58/fzbXmEwmDhw4wMKFC/H29sbX15eNGzfy9u3bf3TfRER6mlpiRER6oS9fvvDx40cA2traqKqqIjU1lYkTJxIQEPDDOa9fv2b58uW0tLQQHR3NsGHDuHr1KikpKdTU1LBp0yZz7fv374mNjWXZsmWEhISQlZVFUlISBQUFNDQ0EBcXR1VVFdnZ2Tg4OLBz504A8vLySE5OJiAggA0bNtDW1sa1a9c4ceIEjo6OrFmzBoBt27Zx5coVVqxYgbu7Oy9evCA7O5tXr16RmZkJwJEjRzh48CBRUVGMHz+eN2/ekJWVRXl5OVeuXMHOzu6/ucUiIv8zCuwiIr1QaGhot7F+/fqRlZWFvb39D+fs3buX+vp6Lly4wMSJEwGIiooiLi6OzMxMQkNDGTduHAD19fWkpKQQHR0NgJubG6tWraKmpobi4mLzMyoqKrh37575GZmZmfj4+HDo0CFsbGwAiIyMxGg0UlJSYg7sBQUFhIWFsX79evNcR0dH7t69S2NjI/3796egoIDZs2eTnJxsrhkxYgTnzp2jtrYWDw+Pv71/IiLWRIFdRKQX2rNnDy4uLkDXCXttbS3Z2dlERUVx9OhRZs2aZVHf0dHBrVu38PPzM4d1AFtbW2JjY7l58yY3btwwB3aABQsWmK9HjRoFgL+/v8UXAjc3Nx4+fGj+fPnyZZqamsxhHeDDhw8MGjSIr1+/msdcXV0pKirCy8uL+fPnM2jQIOLj44mPj7eoKSsr4/Tp0wQFBeHi4kJERAQRERF/c9dERKyTethFRHqhKVOm4Ovri6+vL3PmzCEyMpIzZ87g4OBAampqt/pPnz7x9etXRo8e3e3emDFjgK7+9u8NHTrUfP2t/eT7sW/jJpPJ/LlPnz6Ul5eTlJREREQEvr6+zJ49m+fPn1vUbdmyBZPJxObNmzEYDERFRXHq1CkaGhrMNZs2bcLZ2ZkdO3bg5+dHWFgYBw8e5P379z+zVSIiVk+BXUTkN+Hs7MyMGTOorq62eHkTsAjL/66zsxOgWyvNH390/5H2+5PzH8nIyGDlypVUVFTg6enJ2rVruXTpEtOmTbOoMxgM3Lx5k7179xIYGEh1dTU7d+4kODjY3Js/YcIESkpKOHToEEuXLqWuro79+/ezaNEiqqqq/nQdIiL/T9QSIyLyG/kWvm1tLc9rhgwZgqOjI9XV1d3mvHz5EuhqQfkVtbW1HDt2jCVLlpCenm5xr66uznzd2tpKRUUFrq6uBAUFERQURGdnJydPniQ9PZ3CwkIiIyN5+vQpAwYMwGg0YjQaASgqKiIhIYG8vDwSExN/ab0iItZCJ+wiIr+Juro6SktL8fT0ZODAgRb37Ozs8Pf35/79+zx58sQ8bjKZOH78ODY2NsydO/eXnv/tVH/s2LEW47dv36ampob29nag64XW8PBwjh49aq6xtbXF29vbfN3R0cGKFSvYsWOHxd+aPHmyuUZEpLfQCbuISC90/fp1nJ2dga7Q/e7dO3Jzc2lqaiIhIeGHczZs2EBZWRkxMTHExMQwbNgwrl27RmlpKStXruwWtH/W2LFjGTlyJEeOHKGlpQVXV1cePXpEfn4+ffv2pbGxEYDhw4cTHBzM2bNnaWpqwsfHh/r6es6cOYOLiwuBgYHY29sTExPD4cOHWb16Nf7+/jQ3N5OTk4ODgwNhYWG/tFYREWuiwC4i0gt9+7/n0HV67uTkhLe3N2lpaRgMhh/O8fDwIDc3l3379nH+/Hmam5sZM2YMaWlpLF269JfXZG9vz7Fjx9i1axdZWVmYTCY8PDxISkqivb2dtLQ0ysvL8fLyYvv27bi7u1NYWEhhYSEODg4YDAYSEhIYMmQIAOvWrWPw4MFcvHiR3bt3Y2dnx5QpU9izZ4/5RVkRkd7AxvRnbxqJiIiIiEiPUpOfiIiIiIgVU2AXEREREbFiCuwiIiIiIlZMgV1ERERExIopsIuIiIiIWDEFdhERERERK6bALiIiIiJixRTYRURERESsmAK7iIiIiIgVU2AXEREREbFi/wJAXfVrumCYkwAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig1, ax = plt.subplots()\n", "\n", "ax.plot(ss,vhat0(cc,ss),label='Grow')\n", "ax.plot(ss,vhat1(cc,ss),label='Clear-Cut')\n", "\n", "ax.set(title='Conditional Value Functions',\n", " xlabel='Biomass',\n", " ylabel='Value of Stand')\n", "\n", "ax.legend()\n", "\n", "vcrit = vhat(cc,scrit)\n", "ymin = ax.get_ylim()[0]\n", "ax.vlines(scrit, ymin,vcrit,'grey',linestyles='--')\n", "ax.annotate('$s^*$', [scrit,ymin])\n", "ax.plot(scrit,vcrit, 'wo')\n", "\n", "print(f'Optimal Biomass Harvest Level = {scrit:.4f}') " ] }, { "cell_type": "markdown", "metadata": { "Collapsed": "false" }, "source": [ "### Plot Value Function Residual" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "Collapsed": "false" }, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAucAAAGZCAYAAAAqx5/kAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABzvUlEQVR4nO3dd1xV9f8H8Ndl772XuC4gIODGEeYk90hz5ygztaEN0+pnQ8ss0zS/lWU50koyNTVz79wbxQGobAHZ697Lvef3x5WbxPAC93Iv8Ho+Hj6Ezzn33ve9h/Hicz9DJAiCACIiIiIi0jkDXRdARERERERKDOdERERERHqC4ZyIiIiISE8wnBMRERER6QmGcyIiIiIiPcFwTkREWsVFwYiI1Gek6wKIqGH7448/MH/+/ArthoaGsLOzQ5s2bTB27Fj07t27To9z5swZTJo0CT179sR3331X7rHHjx+P//u//6vT/euLpKSkGr1W586dg42NjRYrqr3c3FwsX74c3bp1Q9++fVXtfn5+sLCwwKVLl3RYnVKvXr2QnJxc6TFTU1NYW1vDz88Pw4YNw5AhQ+qtrokTJ+Ls2bP4/fffERwc/MTzBw0ahDt37uDgwYPw8vKqhwr16zoSNSYM50SkEd7e3ggNDVV9LpFIkJGRgdOnT+P48eN466238MILL+iuwAbGwsJCrZBubGxcD9XUzpIlS/DHH3+gS5cuui7lifr06QNzc/NybQUFBYiLi8PJkydx8uRJxMXFYc6cOTqqkIiaCoZzItKIDh06YMmSJRXa79y5gzFjxmD58uUYNGgQ3NzcdFBdw2Nvb48vvvhC12XUSVXDWf766y8YGOjXqMr58+dX2uMsCAI2bdqEjz/+GN9//z1GjRpVLz3Tn332GYqLi+utF5yI9Id+/XQkokandevWGDRoEEpLS3HkyBFdl0N6oGXLlmjevLmuy1CLSCTChAkT0LZtW8jlcpw8ebJeHtfDwwMtW7aEqalpvTweEekPhnMi0jpXV1cAQFFRUYVjN2/exKuvvorw8HAEBQWhf//++Oqrryo9Vx2rVq2Cn58fzpw5g23btmHIkCFo27YtIiIisHz5csjlchQWFmLRokXo3r072rVrhzFjxuD8+fMV7quwsBDffvstRo4ciQ4dOiAoKAg9evTAG2+8gTt37lT5uH/99RdGjhyJkJAQdO7cGW+88QYSExNr9Xxq6s6dO3jttdcQHh6Odu3aYfr06YiPj8c777wDPz8/JCUlqc718/NDWFhYpfczaNCgCucDwIULF/D6668jIiICQUFBaNeuHUaMGIF169ZBLpeXu+9t27YBAF577TXVa1Pd4964cUNVe1BQEHr37o3Fixfj4cOH5c47c+YM/Pz8sGrVKly/fh0vvvgiOnTogLCwMEyaNAlnz56t3YtXDXd3dwDKcfSPUygUiIqKwqhRoxAWFoZ27dphwoQJOHDgQKX3c+DAATz//PPo1q0b2rZti/79++PTTz9FVlZWufMmTpwIPz8/XLt2rVz77du38eqrr6Jr164ICwvDjBkzEB8fX+ljafP6EpH2cFgLEWldWcBo27ZtufYDBw7g9ddfh0wmQ2BgIDp06IDr16/jf//7H44cOYINGzbA2tq6Vo+5Zs0anDhxAmFhYejSpQvOnDmDb7/9FsXFxbh48SLu37+Pdu3aISMjA5cuXcLkyZOxfft2tGrVCoAymI8ZMwa3b9+Gp6cnOnfuDIlEgujoaOzatQuHDx/Grl274OHhUe5x169fj4MHD8LPzw9PPfUULl26hF27duHcuXPYv3+/VntCL168iBdeeAGFhYUIDg6Gu7s7zp49i7Fjx6J169Z1vv8//vgDCxYsgJGREcLCwhASEoLk5GRER0fj+vXrSEpKwnvvvQcAGDx4MC5fvozExER06NAB7u7ucHJyqvK+9+zZg7feegsymQwhISFwc3PDjRs3sGHDBuzduxfr16+v0Nt++fJlrFmzBo6OjujcuTMSEhJw5swZXLx4Eb/99hsCAwPr/JwB5R+VZX+8lX19AMohL2+++SZ2794NKysrtGvXDiKRCOfOncOsWbPw6quvYtasWarz9+3bh1dffRUmJibo0KEDLC0tER0djXXr1uHIkSPYvn17hXHvj7t48SKmTZuGoqIi1Wt07tw5jB07FoaGhnV+njW5vkSkRQIRUR1s3bpVEIvFwrx588q1S6VSISUlRfjyyy8FsVgsvPTSS+WOp6enC+3atROCg4OFY8eOqdplMpnw/vvvC2KxWHj33XdV7adPnxbEYrEwffr0Co/94YcfqtpWrlwpiMViQSwWC3v27FG1HzlyRNX+zDPPCOnp6apjb7zxhiAWi4XPP/9c1fbNN98IYrFYmD9/viCXy1XtBQUFwoQJEwSxWCz873//q/C4fn5+wu+//17u/KFDhwpisVjYtm3bE1/PxMREQSwWC08//fQTz32cTCYT+vfvL4jFYuHnn39Wtefk5AijRo1SPffExETVMbFYLISGhlZ6fwMHDix3flFRkdC+fXshNDRUuHXrVrlzDx48KIjFYiEkJESQSqWq9nnz5lW4DpU9blpamhAcHCwEBgYKBw8eVLWXlpYKS5YsEcRisTBs2DBBoVAIgvDv14JYLBbef/99QSKRqG7z7rvvqq6bOp5++ukKr4sgCIJCoRDy8/OFM2fOCGPHjhXEYrEwZMgQobS0VHXOpk2bBLFYLIwePbrc11NiYqLQp08fwc/PTzh37pyqvXfv3kKbNm2EuLg4VZtUKhWmTp0qiMViISoqStVe9jV29epV1WsxYMAAQSwWC1u2bFGdl5+frzq3vq9vdfdPRLXHYS1EpBHbtm2Dn5+f6l9QUBB69uyJb7/9FsHBwVi+fHm583///XcUFBRg2rRp6NGjh6rdyMgI7777LpycnLB9+/YKwwjU1alTJ0RGRqo+j4iIgIWFBQBg9uzZcHZ2Vh3r168fACAhIUHVZm5ujqeeegqvvfZaucmLlpaWqiX1UlNTKzxuly5dMHLkyHLnDxs2DAAqDFGoTnJycrnXs7J/q1atUp1/7tw53L17F506dcL48eNV7ba2tli6dClEIpHaj12ZzMxMRERE4KWXXoJYLC53rFevXnBzc0NxcTGys7NrfN+//fYbJBIJJkyYgF69eqnaDQ0N8fbbbyMgIAA3btzA6dOny93OwsICCxYsgImJiaqt7LnX5LUGgN69e5d7bf39/dG+fXtMnDgRFy5cQJcuXbBmzZpyPdTr1q2DSCTC0qVLy309eXl5Yd68eRAEARs2bFC1Z2RkwMjICA4ODqo2Y2NjzJs3Dx9//HGVQ1AA5XCT2NhYhIeHY9SoUap2KysrfPLJJ3p9fYmoZjishYg04r9LKcrlcuTk5CA6OhrXrl3D+PHjsWbNGtXQhnPnzgEAOnfuXOG+TE1N0bFjR+zZsweXL19GREREjev57xAaALCzs0NRURH8/f3LtZetEy6RSFRtzz//PJ5//vly5+Xl5eHmzZuq2mUyWYXHCAkJqdBWFtxqMo5enaUU/fz8VB9fuHABAPDUU09VOM/X1xdt2rTB9evX1X78//L29sayZcvKtSkUCty/fx9XrlxRvXaVvSZPUjZkpOyPpMeJRCJERkYiJiYG58+fR3h4uOqYn58fzMzMyp1f9vVV0zkLjy+lmJ2djVOnTkEul+Ppp5/GG2+8UWFY0IMHD3D//n14eHigWbNmFe6va9euMDAwUH2tAMoVjU6cOIERI0Zg1KhR6NmzJwICAiAWiysE4v8qe426d+9e4Zi3tzfEYjFu3bpVo+f83/vQ1vUlopphOCcijahqKcWioiLMmzcP+/btw7x587B27VoAQFpaGgBUCMD/9eDBg1rVY2trW6GtrHfRzs6u0vbKHvvnn3/GuXPncO/ePVWvYdn5QiVLBVY2Rr6st7Wy86tS06UUMzIyAKDKpSq9vb3rFM4BZf0HDx7Ejh07cPv2bSQnJ6vCWnWvyZOU1e7p6Vnp8bL2svPKVLb5Um1ea6DiUopxcXGYMmUKDh8+jDZt2lQazgEgJSWl3B9J/5WVlQWZTAZjY2N8/PHHmDVrFm7cuIEVK1ZgxYoVcHFxQZ8+fTB+/Phy49n/q+y5l02u/i9PT886hXNAe9eXiGqG4ZyItMrCwgIffPABDh48iBMnTuDBgwdwdXVVrfzQr1+/aidJlq2SUVN1nSB36tQpzJgxAyUlJfDy8kLnzp3RunVrBAcHIysrC++8806lt6vr8ILaetLjGhnV7Mf9f1fmkMvlmDFjBo4dOwZTU1MEBweja9eu8Pf3R8eOHfHKK68gNja2xnUDTw58CoUCAMoNXwG0+1q3bNkSK1euxLhx47B69Wr4+PiohicB/74+zs7OT9xkqbS0FMbGxvDw8MAff/yBM2fO4MCBAzh+/Dju3buHzZs3Y8uWLfjyyy/Rv3//WtWrz9eXiGqG4ZyItM7R0RH29vbIzMxEamoqXF1d4eLignv37mH69OlqbU9enwRBwPvvv4+SkhJ88cUXGDx4cLnjW7Zs0VFlVSv7I6aqrejT09MrtIlEoiqXx8vPzy/3+Z9//oljx44hLCwM33zzDezt7csdz8vLq03ZAAAXFxfcvXsXSUlJlfYMly335+joWOvHqI3Q0FBMnToV33//PRYtWoTw8HBVfWVDlezs7Gr0DodIJEKXLl1UgT4pKQlr1qzBb7/9Vm04L3vclJSUSo//912FssfSh+tLRDXDCaFEpHW5ubnIyckB8O+wi3bt2gEAjh07Vultpk6dirFjx+L27dv1UuPjsrKykJiYCC8vrwrBHIBqI5qyHl19UDZ2v7L1tXNycnDlypUK7RYWFpBIJCgoKCjXnpycXCHsld1+zJgxFYLbnTt3VOH/8V5wdXu227dvDwDYv39/hWOCIKjaO3bsqNb9adLs2bPh5eWF/Px8fPbZZ6p2Ly8vuLm54d69e5WuYR8TE4PIyEi8/fbbAJSv6dChQ/HSSy+VO8/Lywvvv/8+DAwMVEO9KlMW5g8ePFjh2MOHDxEdHV2hXdvXl4i0g+GciLRKIpFg8eLFKC0tRfv27VXhfPTo0TA1NcWaNWvKBXRBEPD111/j5MmTePDgAVq2bFnvNdvZ2cHMzAxpaWm4efOmql2hUOCnn37C33//DaD8BFJdCw0NRVhYGK5du4bVq1er2qVSKd5///1Kay2bhLhx40ZVW3FxMT788MMK55ZdtyNHjpQLaImJiZg7d67q88cfp2wYyn97af9r9OjRMDMzw88//4xDhw6p2hUKBZYtW4aYmBj4+fmp/qCrT2ZmZliwYAEAYPfu3Th16pTq2IQJEyCTyfDWW2+Ve2ciOzsb7777Lu7evataB9/DwwM5OTk4duxYuecIAH/99RcUCkW17yCFhISgbdu2uHLlCr799ltVu0QiwbvvvlvpRE1tX18i0g4OayEijTh//jzefPPNcm3FxcU4d+4ccnNzYWdnVy4UeHp6YtGiRXjnnXfw4osvok2bNvDy8sLt27dx7949WFhYYMWKFRrZXKWmDA0NMX78eKxduxajRo1Cp06dYGZmhuvXryM1NRWtWrVCbGwsMjMztVZDdnZ2hdezMkOHDlUtRbl06VKMHz8eK1euxN9//43mzZvj6tWrqte/7N2LMlOmTMGlS5ewYsUKHD58GM7Ozrhw4QKMjIwQHh5eLogOGzYMP/74I/bs2aMKy1lZWbh48SIMDAzg4+ODhIQEZGRkqCY2+vr6AgCWLVuG48ePY8qUKZUuF+ju7o5PPvkEb7/9Nl5++WWEhITA3d0dN27cQEJCAtzc3LB8+fJyS1rWp969eyMiIgJHjx7Fxx9/jB07dsDY2BhTp07FhQsXcPjwYURGRqJt27YwNTXFuXPnUFhYiI4dO+Lll18GoHwX4aOPPsKMGTPw8ssvo23btnB3d0dqaiquXr0Kc3NzzJs3r8oaRCIRlixZgkmTJmH58uX466+/4Ovri8uXLyM/Px+tW7eusGuttq8vEWkHe86JSCMSExOxc+dO1b9du3bhn3/+gaurK6ZMmYI///yzwooXQ4YMwa+//op+/fohLS0Nhw8fhkKhwIgRI7B9+/ZKl0OsL3PnzsU777yDZs2a4cKFC7h8+TJcXV0xf/58bN++HU5OTrh8+XKFbdc1paioqNzrWdW/uLg41W18fHywdetWPPvss3j48CGOHj0KHx8fbNiwodw63GX69++P1atXIzQ0FLdu3cL58+fRrVs3/P777xVWfXFzc8OmTZvQu3dvFBYW4tixY8jIyMAzzzyD33//HVOmTAEAHD58WHWbMWPGYNCgQZBKpTh+/Hi1Q5QGDhyo+lpISEjAoUOHYGBggBdeeAHbtm3TyTsoj3v33XdhYmKCuLg4rFu3DoDyj7jVq1dj4cKFaNmyJa5cuYLz58+jWbNmmD9/PtauXVtusnNERATWrl2L7t27IyEhAQcPHkRaWhqGDBmCP/7444lzL1q2bIktW7ZgxIgRyMrKUl3f9evXw8fHp8L52r6+RKQdIqGBDiB77733cP/+/XJv1wHKgPDZZ5/h7NmzAICePXvinXfeKbfpAxFRUzNo0CDcuXMHBw8eLLdkIBER6ZcGOawlKioKUVFR6NSpU7n27OxsPP/885BKpXjhhRcgl8uxdu1a3Lp1C1FRURWW4SIiIiIi0icNKpzL5XJ88803+Prrrys9vm7dOqSlpWHnzp2qt0BDQkIwZcoUbN++HaNHj67PcomIiIiIaqTBjDmXSCQYPnw4Vq1ahaFDh1a6Fu7u3bvRqVOncmMTu3btiubNm2P37t31WS4RERERUY01mJ7zsrValy9fjgEDBqBXr17ljufm5iIxMbHSDRwCAwNx5MiReqqUiEj/7Nq1S9clEBGRGhpMOLeyssK+ffuq3KL4wYMHAFBpj7qzszMKCgqQn58Pa2trrdZJRERERFRbDSacGxgYVLvGbWFhIQDA3Ny8wrGypayKioqqDecKhQJyeYNcvIaIiIiIGghj46r38Ggw4fxJ1NlG+0kbWMjlAnJyijRVEhERERFRBc7OVXcWN5gJoU9iaWkJoPKthcvays4hIiIiItJHjSace3h4AAAyMjIqHEtPT4eNjQ0sLCzquywiIiIiIrU1mnBuY2MDLy8vXL9+vcKxGzduICgoSAdVERERERGpr9GEcwDo168fTp06hbi4OFXbP//8g7t372LAgAE6rIyIiIiI6MkazYRQAHjxxRexY8cOTJ48GVOnToVEIsEPP/yAwMBADB06VNflERERERFVq1H1nDs4OODnn3+Gv78/Vq5cifXr16NPnz744YcfYGJiouvyiIiIiIiqJRIEgQt7PyKTybmUIhERERFpVZNYSpGIiIiIqKFjOCciIiIi0hMM50REREREeoLhnIiIiIhITzSqpRSJiIiofikEAQ8LpUjLk+BBvvJfeoEE6fkS5BTLEOZli0kdvWFmbKjrUokaBK7W8hiu1kJERPQvQRBQIJEjLb8ED/IlqgCe9iiEP8grwYMCKeSK6qOEn4sVlg5pAw9bs3qqnEi/VbdaC8P5YxjOiYioKZGWKpBeoAzdlQXw9HwJCqXyGt+vqZEBXK1NUSiV42GhFABga2aExYMC0LmZvaafBlGDw3CuJoZzIiJqTEpkcqTlSZCcV4LU3BKk5pUgJVeC1Dzlx1lFshrfp6EIcLYyhZuNKVytTeFqbQZX68c/N4WtmRFEIhFK5QqsPHYXv1xMBgAYiIBZ3ZtjYkcviEQiTT9dogaD4VxNDOdERNSQSEoVqqCdmluC5MeCd0pu7cK3nbkx3B6F7McDt/JzMzhZmsDQoGbBek/MAyzedweSUgUAoI/YCe/394OFCcehU9PEcK4mhnMiItInpXIFUvMkSM4tRkqe5LHe7xKk5ElUQ0bUZWQggpuNKdxtzOBuYwo3azO4PgrgZYFcWxM3b6UX4O0d15GSJwEAtHC0wBdDA+Ftb66VxyPSZwznamI4JyKi+lYklSMppxhJuSVIzilGUk6J6vMHeSWQ1+C3tKGBCK7WpvCwNYPHoxCu/NgM7ra16/XWpJxiGd7ffROn72cDAKxMDfHxAH90b+Gos5qIdIHhXE0M50REpGmCICCrSIaknGIk5z4K3jklSMopQXJucY2GnhiKAFdrU7jbmimD96Pw7W5rCg8bMzhZmcJIh+FbHXKFgG9O3sP6s4kAABGAF7s2w7QuPjDgOHRqIhjO1cRwTkREtSEIAh4WyZCQXYSErGIkZBcj8VEYT84pQZFM/RVPrE2N4GVnBk9bc3jZmak+9rA1g4u1/odvdR26nYEP/76tem16tHDAh8/4w9qMW7BQ48dwriaGcyIiqk6BpBQJ2cWP/hU99nFxjZYcdLEygaedObxszeBlpwzhZZ/bmhtr8Rnol/iHhXhrxw0kZBcDADxtzbBkcAD8XasOLkSNAcO5mhjOiYhIWqpAUm6xqge8LIjfz1Z/CIqhgQietspeby9bc3ja/RvCPWzMuFvmYwokpfhgzy0cjXsIADAxFOGNXq0wPNiNyy1So8VwriaGcyKipiOvRIa7D4twL6sIdx8W425WIe5lFSMtrwRP2PBSxdXaFD725qp/zewt4GNvDndbs0Yz/KQ+CIKAn88nYfXxu6oJsAPauOCdPq1hzj9kqBFiOFcTwzkRUeMiCAIyC6WPhfAi3H30v7q94HbmxuUCeNk/bztz9oBr2KWkXCzYFYPMR0tEtnC0wGdD2sDXwULHlRFpFsO5mhjOiYgaJoUgICW3pEIIv5dVhALJk8eCGxmI4GNvjuaOFmjmYIFmjwXwpjQGXB88LJTivb9u4nxCDgDAwtgQ7/ZrjX7+LrotjEiDGM7VxHBORKTfylZFicsoRNzDQsRmFCLuYRHiMwtR8mj3yeqYGxvA18ECzR0t4OtggRaP/ve0M+cwFD0iVwhYc+o+fjydoGobHeqB1yJawMTIQIeVEWkGw7maGM6JiPRHobQUcZlFiM0sRHxmIWIzCxGXWYSc4icPR7EzN0ZzB3P4PgrfzR0t0NzBAq7Wppxk2ICcjM/Cwj03kVtSCgAIdLPGp4MD4G5jpuPKiOqG4VxNDOdERPWvVK7AvexixGaUBXDlv9RH27xXx8bMCC2dLNHKyRItnf4N4fYWJvVQOdWH1LwSzN8Zg+tp+QCU1/z9fmL0bO2k48qIao/hXE0M50RE2lUoLcWd9ELczijA7Uf/x2UWQvqEPepNjQzQwtECLZ0sH4Vx5cdOlibsCW8CZHIFvjoaj98upajangvzwKtPcZgLNUwM52piOCci0gxBEJBRIC0Xwm+nFyAxp6Ta2xmIAB9788dCuPJ/T1szGHJMeJN38HYGFu27rZrk6+dihcUD/dGMq7lQA8NwriaGcyKimlMIAhKzixHzQBnAywJ59hPGhluZGkLsbAU/FyuIXSzR2skKvo4WMGVPKFUjObcY7+2+iehU5TAXC2NDzOvTCgPauOq4MiL1MZyrieGciKh6giAgObcEMQ8KEJOWj5gH+Yh5UPDErevdbUwhdlaGcOX/VnC34eRMqp1SuQLfnLyHDeeSVG0DA13xdq9WsDDh2vOk/xjO1cRwTkT0L0EQkJYvQUxaPm48CuM30wuQ92jljMoYGojQwtECYmdLiF2UveKtnS1hY8a1wknz/rmbhQ/23FK9S9PM3hyfDg5Aa2crHVdGVD2GczUxnBNRU5ZVJEV0aj6up+U/6hUvqHbZQkMR0NLJEgFu1mjjaoUAN2u0dLTkBD2qVxkFEvzfnluqTYtMDEWY07MlRoa4850Z0lsM52piOCeipkImV+B2RiGiU/JwLTUP0an5SM6terKmgQho7miBAFdrBLhao42bFVo5WXL7etILcoWAdWcTsOaf+1A8SjW9Wjvh3X6t+a4N6SWGczUxnBNRY/UgX4Lo1DxcTVEG8ZsP8qtcvlAEoJmDuTKIP+oVF7tYwZxBnPTcpaRcvLc7BukFUgCAi5UJPhrgj/bedrotjOg/GM7VxHBORI2BtFSBmAf5qiAenZqnCiuVsTM3RpC7NYLdbRDkbo02btawMjWqx4qJNCenWIaP997GsbiHAJR/bE7s6I0Z3ZrB2JBDrkg/MJyrieGciBqi/JJSXE3Nw+WkXFxJzsX1tKp7xQ0NRBA7WyqDuIcykHvamnFsLjUqgiBg29VUfHkkHpJSBQDA38UKHw/0hy/XRCc9wHCuJoZzImoIHuRLcCU5F5eScnElJQ+xGYWo6ge5i5UJgj1sEORug2B3a/i5WHGcODUZ9x4W4b2/buJWegEA5U6zc3u2wPC2nCxKusVwriaGcyLSN4Ig4F5WMS4l5eBSch6uJOciNU9S6bkiAK2cLRHmaYsQTxuEeNrC1dq0fgsm0jMyuQLfnryHjeeSVH/EPtXSEe/1aw17CxOd1kZNF8O5mhjOiUjXBEFAQnYxLiTm4EJiLi4k5eJhYeXjxU0MRQh0t0HYoyDe1sOGY8WJqnA+IQcL99xUzb9wtDTBwkgxwn0ddFwZNUVNKpwnJibis88+w9mzZwEAPXv2xDvvvAMHhyd/8zGcE1F9K9tx83xCDs4n5uBiUi4yqpi8aWNmhBAPG4R62iLUyxb+LlZcU5yoBnKLZfj0wB0cvJ2pahvTzhOzezSHKb+XqB41mXCenZ2NkSNHQiqVYtKkSZDL5Vi7di08PT0RFRUFE5Pq375iOCei+pCWV4KzCTmq3vEH+ZUPU7E1M0I7bzt08LZFO287tHC0gAHHyRLViSAI2HX9AT4/FItimXKyaAtHC3z4jB/8XasOTESaVF04b1Tvf65btw5paWnYuXMnWrZsCQAICQnBlClTsH37dowePVrHFRJRU1QgKcWFxFycvZ+N0/ezkZBdXOl51qZGaOdli/Y+dmjvZYtWzpYM40QaJhKJMDjIDaGetvi/PTcRnZqP+IdFmLz5Ml7o4oPJnX1gZMDvO9KdRtVz3qdPH3h5eWHdunXl2iMjI+Hq6or169dXe3v2nBNRbZlu3QLLxR/CIDkJck8v3Jr5Nv4M7Ikz97NxLTUfckXFH7WWJoYI87JFB287tPe2RWtnKxgyFBDVm1KFgHVnEvDD6QTV9+jMpNN47fA6mKQmQ+HphcJ3F0Iykp17pFlNouc8NzcXiYmJ6N+/f4VjgYGBOHLkSP0XRURNgunWLbCa8woMSpQ94kZJiWj+f28iNXI2Lgc+rTrPyECEth426OJrj04+dvBztWYPHZEOGRmI8EJ4M3Rv4YD/23MLwcd2Y/bfX8O0VDnUzDApEdZzXwEABnSqN40mnD948AAA4OrqWuGYs7MzCgoKkJ+fD2tr/RtPtjM6DbuuP6j2HLGLFd54uqXq81vpBfjycNwT7/u750LKff7Sb1eeeJu5T7eEn4uV6vNlh+Nw+9EasVUZFOiKwUFuqs/5nJT4nP7V2J6TtFSB5zddQl6JDNs+XwCbkvJDVSxKJXj72Aac7/oMnmrpiM7N7NDOyw7fnLyH0/eycfpett49J6DxXSeAz6kMn5NSVc/J1swI849vhEVp+TkgouJimHy0EJNL/Z5Yn749p8c1luv0OE09J33TaMJ5YWEhAMDc3LzCMVNT5Tq/RUVFehnOU/NKcDEpt0a3KZCU1vg2ANS6TYGktNznt9MLnni79t625T7nc/r3cfmclBrDc0rPl+Dk3Sz8czcLZ+/noEgmBwC452VUer5HXib+r78Y7b3tVG369pz+qzFcp//ic1Lic1Kq7jm55lb+vWyaltJgn1N1j8vnpJ8aTThXKBRPPMfAQD+XSXK3MUM7L9tqzxE/9hcoAFg9mjhWU+rc5r/rJP/3sSvjbmNW4XM+Jz6nxzXE5yQAKJSUIrekFHklMgxcc6bS+0uzdYZHbnqF9kx7F717Tk/SEK/Tk/A5Vf7YlWnqzynT3gUu2RV7a5OtnWBrZgQfe3MYG6qfJfThOVV1vw35OlV1v7V9Tvqm0UwIvXnzJoYOHYr3338fEyZMKHdsyZIl+Omnn3Dp0iVYWFhUeR+cEEpERVI5Tt/PxtHYTJyMz0JuSWmFc+zMjdG1uT26NXdAF197OO/eBuu5r0BU/O/QFsHcHPlfruI4VaIGxHTrlgrfy8XGppjXfzb+DHwatmZGeLNXK/T3d4aIKylRHTSJCaEeHh4AgIyMim9Jpaenw8bGptpgTkRNV1aRFMfjHuJI7EOcvZ8Nqbxin0WAqxW6NXdA9xYOCHCzLrfEYVkAL1uthSs8EDVMlX0v5771HiQ2ocDNDOSWlOL9v25i/60MzO/TCk5WprotmBqlRtNzDgC9e/dG8+bN8cMPP5Rrj4yMhLu7O3766adqb8+ec6KmIyG7GEdjM3E09iGupuThvz8IzYwMEN7cAT1aOKBrcwc4Wla/iRkRNW4Hb2fgswOxyC6WAVDuSzD36RYY2MaVvehUY02i5xwA+vXrhw0bNiAuLk61CdE///yDu3fvYtq0aTqujoh0SRAE3E4vxME7GTgS+xB3H1b8Q9ze3BhPtXRERCtHdPSxg5mxoQ4qJSJ91FvsjPZedvjicCz23sxAvqQUH/59G/tvZWBBXzFcrdmLTprRqHrOs7KyMGjQIBgaGmLq1KmQSCT44Ycf4OPjg19//RUmJtX3fLHnnKhxKQvkB25n4ODtDCTmlFQ4x9vODBGtnNCzlSOC3G24CRARPdHR2Ex8eiAWDwulAJQbir0W0QLDgt3Yi05qqa7nvFGFcwCIj4/Hp59+ivPnz8PMzAwRERF4++234eDg8MTbMpwTNXyCIOB2RiEO3Ko6kAe4WuHp1k6IaOWI5g4W/GVKRDWWVyLDl0fisfuxtbg7+djh3X5ieNiaVXNLoiYWzuuC4ZyoYSoL5AdvZ+DArcoDeZC7NXqLndFb7FRhyS4ioto6GZ+FT/bfRnqBshfdwtgQs59qjpEh7uUmjhM9juFcTQznRA1LUk4x9t5Mx54b6bifXVzhOAM5EdWHAkkpVhyNx45raaq2dl62WNC3NZo5cKU4qojhXE0M50T6L6dIhv23M/B3TDqupuRVOM5ATkS6cuZeNhbtu420fAkAwMRQhGldmmFiR68abV5EjR/DuZoYzon0U4lMjmNxD7EnJh2n7mVDrij/Y6u1syUi/V3Q19+ZgZyIdKpQWor/Hb+HqMspqiVaWzpZ4N2+YgR72Oi0NtIfDOdqYjgn0h8KQcDFxFzsuvEAh29nokgmL3fc1doU/f1d8EyAC1o5W+qoSiKiyl1NycPifbcR/2jZVhGAZ0M9MLO7b4Xt7anpYThXE8M5ke6l5ZVg1/UH2Hn9AVJyy0/stDY1Qm+xEyIDXBDmZcvJVkSk12RyBTacS8Ta0wmQPdp52MXKBG/3bo2IVo46ro50ieFcTQznRLohKVXgaGwm/oxOw9n7OeV26zQ2FKF7C0c8E+CCbs0dYGLEcZtE1LDcyyrCp/vv4GJSrqqtV2snvNWrJZysuHlRU8RwriaGc6L6IwgCbqYX4M9raard9h4ndrbEkCA39A9wgZ25sY6qJCLSDIUg4M9raVh57K7q552VqSFeeUq5eRHfCWxaGM7VxHBOpH1FUjn23kzH1iupuJVeUO6YrZkRIgNcMDjQDX6uVjqqkIhIezILpVh2KA4Hbmeo2kI9bTCvT2u0cuL8maaC4VxNDOdE2hObUYitV1KwJyYdhdJ/J3caiIDOzewxJMgNT7V05LAVImoSjsc9xGcHY/Hg0bKLhgYijGvniRfCm8HCxFDH1ZG2MZyrieGcSLMkpQocupOBrZdTceU/a5K7WJlgWFt3DAlyg6s1x1wSUdNTKC3Fdyfv47dLyShbIdbV2hRvPN0SPVs5QsShLo0Ww7maGM6JNCMppxh/XEnFn9FpyC35dyy5CEB4c3uMaOuBbi0cYGTAXzxERLfSC/DZgTu4lpqvauvewgFv9moJT1tzHVZG2sJwriaGc6LaEwQBFxJz8evFZByLe1huxRV7c2MMCXbD8LZu/EVDRFQJhSBgx7U0fH38LvIedWqYGhlgWhcfTOjAHUYbG4ZzNTGcE9WcpFSBvTHp+PVSMu5kFJY71s7LFiND3PF0ayf+YiEiUkN2kRSrjt3FzusPVG2+DuaY17s1OvjY6a4w0iiGczUxnBOpL6NAgt+vpGLblVRkF8tU7aZGBngmwAXPtfPkygNERLV0OSkXSw7eQVzmv7kkMsAFr0e0gKOliQ4rI01gOFcTwznRk91OL8DP55Ow/1YGShX//vhwtjLBqFAPDA92h50F1yUnIqqrUrkCv1xMxven7qNYpgCgXBv9pa6+eDbUg/N2GjCGczUxnBNVrmw8+fpziTh9L7vcsSB3a4xt54lerZ1gxKErREQal5ZXgmWH43Ak9qGqraWTBd7q1Qrtve10VxjVGsO5mhjOicqTKwQcjc3E+nNJuJH27yoCBiKgt9gZY9t5ItjDRocVEhE1HSfiH2LZ4Tgk5ZSo2vqInfFaRHO42ZjpsDKqKYZzNTGcEylJShX468YD/Hw+CQnZxap2UyMDDAlyw/gOnlx1hYhIBySlCmy+kIQfTyegpFQ51MXMyABTOvtgfAcvmHIjtwaB4VxNDOfU1JXI5Pjjaio2nEvCw0Kpqt3GzAijQj3wXJgH7C04EYmISNfS8kqw8thd7L+VoWrztDXD3KdbokcLB25gpOcYztXEcE5NVbFMjt8vp+Dn80nIKvp35RVXa1OM7+CFoUFu3E6aiEgPXUjMwReH4hCb+e9Stl2b22Nuz5Zo5mChw8qoOgznamI4p6amSCpH1KNQnvPYcojN7M0xpbMP+vs7c5InEZGeK1UI2Ho5Bd/9cx/5EuUGRkYGIoxr74mpXXxgaWKk4wrpvxjO1cRwTk1FgaQUUZdTsOl8EnIf7UQHAM0dLDCtiw/6+DnDkEt0ERE1KNlFUvzvxD3suJam2qXZydIEM7v7YmCgKww41EVvMJyrieGcGrsSmbKnfP3ZxHKhvIWjBV4Ib4ZerZ0YyomIGrgbafn44lAsrqX+u8pWgKsV5vRsiTAvWx1WRmUYztXEcE6NValcgR3RafjhVAIyH5vo2drZEi908UHP1k7sUSEiakQUgoA9N9Kx+sRdZBT8+3O/V2snvPJUc3jZccUtXWI4VxPDOTU2CkHA3pvp+O7kfSTn/rsubnNHC8zo5ouerRwZyomIGrFimRwbziZi4/kkSB4tvWhsKMLYdp6Y0tkHVqYcj64LDOdqYjinxkIQBByPz8I3J+6Vm8HvYWOK6V19ERngwuErRERNSFpeCf534h72xKSr2uzNjTGjuy+GBrnxd0I9YzhXE8M5NQY30vKx4mg8LiXlqtocLIwxrUszDG/rBmOuvkJE1GRdT83DssPxuJaap2pr5WSJ13u2QOdm9jqsrGlhOFcTwzk1ZJX1ilibGmFiRy+MaecJc2OuU05ERMp3V/ffysCqY3eRli9RtXdv4YDXnmoBX0euj65tDOdqYjinhqhQWooNZxOx6UKyajyhkYEIo8M8MK2LD2zMjHVcIRER6aMSmRy/XEzGujOJKJLJAQCGImBYW3e8EN4MTpbcEVpbGM7VxHBODUmpQsCf0Wn47uS9crt69hY7YXYPzsQnIiL1ZBZK8c2Ju9gZ/UC1Prq5sQEmdvDG+A5e3CFaCxjO1cRwTg3FpaRcfH4oFncy/p3s2cbNGnMiWiCUa9gSEVEt3MkowKpjd3HqXraqzcHCGC91bYYhwe4w4qRRjWE4VxPDOem7zAIJVh67W25cuZu1KWb1aI5+/s5cFpGIiOrszP1srDp2F7fSC1Rtvg7mmN2jOZ5q6QgRf9fUGcO5mhjOSV+VyhX47VIKvj91H4VS5bhAUyMDPN/JGxM7eMGMkz2JiEiDyvbJ+ObEPaTm/TtpNNTTBq8+1QLBHjY6rK7ha1ThPCEhAYMGDcL333+Pzp07Vzi+efNmbNiwASkpKWjWrBlmzJiBgQMHqnXfDOekj84lZOPzg3G4m/Xv12bPVo6Y07MlPGzNdFgZERE1dpJSBaIup+CnMwnIKylVtfcWO2Fm9+bwsef8ptqoLpw3qG2h8vLyMGvWLEgkkkqPr127FkuXLkVkZCQmT56M/fv3Y+7cuRCJRBgwYEA9V0tUN1lFUiw/Eo+/HxvC4mNvjjd7tUS4r4MOKyMioqbC1MgAEzp4YUiQK9adScRvl5IhlQs4eDsTR2IfYkRbd0zt4sOVXTSowfScx8XFYfbs2YiPjwcAbNiwoVzPeV5eHiIiItCrVy8sW7YMAKBQKDBx4kQkJyfj4MGDMDSs/q1/9pyTPhAEATuvP8DKo/HIfdRLYW5sgGldmmFce09uIkRERDqTmleCb/6zp4aZkQHGtPPExI5eXL5XTdX1nDeI3/Lbtm3D0KFDkZOTg1GjRlV6zqFDh1BUVISxY8eq2gwMDDBu3Dikpqbi0qVL9VUuUa0lZBdjZtRVfLz3tiqY92zliKgpHfF8J28GcyIi0il3GzN8NMAfP09ohy6PdhQtKVVg3dlEDPvhHNadSUDxozXTqXYaxG/6W7duYeDAgdi5cyfatWtX6TnR0dEAgMDAwHLtbdq0KXecSB/J5Ar8eDoBY9efx/nEXACAs5UJlg5pg8+HBsLV2lTHFRIREf3Lz9UKq54Nxrej2yLYXdkLnC8pxeoT9zB87TlsuZQCmVyh4yobpgYx5nzu3LkwMal+LFN6ejpsbW1hbl5+YoKzszMAICUlRWv1EdXFjbR8fLT3FuIylUOqRABGhrhjVo/msDJtEN+iRETURLX3tsPasaE4Hp+Fb07cQ2xmIR4WSvH5oVhsupCEl7o2Q39/FxhyjXS1Vfmb39/fv8brWIpEIty4cUOtczMyMqo9bmFhAUtLSwB4YjAHgMLCQpiZVVy5oqytuLhYrbqI6otMrsAPpxOw/kwC5I9mfrRwtMCCvq0R4smNhIiIqGEQiUR4qqUjujV3wL5b6fju5H0k55YgJbcEC/fcwvqziXi5my8iWnGNdHVUGc6HDRum1Rewe/fu1R6fMWMG5syZo/b9KRSKSusta+MXA+mTW+kF+PDvW6odPo0MRJjaxQeTOa6ciIgaKEMDEZ4JcEUfsTN2XEvD2tMJyCyUIv5hEd768waC3K3xcjdfdPSxYy6rRpXhfMmSJVp94EWLFlV7PCAgoEb3Z2lpiZKSkgrtZT3mZb3wRLpUKlfgp7OJWHs6AXKFsru8tbMlPoj0g9jFSsfVERER1Z2xoQGeDfXAoEBXbLmUgvXnEpFXUoro1HzM+v0a2nnZYnrXZmjvbafrUvWSRga0KhQKFBYW4vTp0+jbt69at6lq1ZXacnd3R25uLqRSablhMOnpyqV+XF1dNfp4RDUVm1mID/bcUm2HbCgCpnT2wdQuPuwtJyKiRsfM2BCTOnljeFt3/Hw+Eb9cTEaxTIGLSbmYseUqOvjYYUbXZhzK+R9qh/P8/Hx8+umn2LdvH4qKilDV8ugxMTEaK64mAgMDIQgCYmJiEBISUqGe4OBgndRFJAgCfruUglXH4iF9NLi8pZMFFkb6IcC16nVOiYiIGgNrMyO83L05nmvniQ1nk/D7lRRIShU4n5CDFxJy0KWZPaZ3bYZgDxtdl6oX1O6u+/zzz/HHH3+gVatWCA8PhyAIGDRoEMLDw2FkZARTU1OsWrVKm7VWKyIiAqampti4caOqTaFQYPPmzfD09ERoaKjOaqOm62GhFK9vi8ayw3GQygUYiIDJnbyxYXw7BnMiImpSHCxM8HrPFtj+QieMbecJE0PluPPT97Mx9ZfLeO2Pa7ielq/jKnVP7Z7zI0eOoG/fvli1ahWys7MRHh6OiRMnom3btoiJicH48eNVu3fqgr29PaZPn45Vq1ZBEAR06dIFe/fuxfnz57FixYon7g5KpGkn72bho79vIatIBgBwtzHFxwP8+fYdERE1aU6WJpj7dEtM7OiFdWcSse1aKmRyAf/czcY/d7PRo4UDXurqCz/XpjkXS+1wnpWVhW7dugFQBmFXV1dcvXoVbdu2RUBAAJ599lns3LkT06dP11qxTzJr1iyYm5tj06ZN2L9/P3x9ffHVV18hMjJSZzVR0yMpVeDr43fx68VkVVt/f2fM690a1mZct5yIiAgAnK1M8VbvVsqQfjYRO66loVQh4Hh8Fo7HZ6FnK0dM79oMrZ2bVkhXOylYWlpCofh3pycfHx/cvn1b9blYLMbWrVs1W10lRowYgREjRlR6TCQSYdq0aZg2bZrW6yCqzP2sIszfFaNaItHC2BDz+rTCMwEuXDaKiIioEm42ZninT2s838kbP55OwM7rDyBXCDgS+xBHYh/i6dZOmNbFB35NZFUztcect23bFnv27IFcLgcAtGrVCufPn1dNDL17965amwURNVYHbmXg+U2XVME8yN0amya1w4A2rgzmRERET+BuY4Z3+4mxdWoHDAlyxaMh6Th8JxMTNl7E3G3RuJ6ap9si64FIqGrZlf84deoUpk2bBg8PD2zduhX379/H6NGjER4ejmbNmmHr1q3o3bs3VqxYoeWStUcmkyMnp0jXZVADIy1V4Kuj8dhyOUXVNqmjF17u5gsjLpFIRERUK4nZxVh7JgF/33ig2kkbALo0s8e0Lj4I9Wq4c7icnateFELtcA4AJ06cwPr167FmzRqIRCL8+uuv+Oyzz1BcXIyQkBCsXLmyQa8nznBONZWSW4L5u2Jw49HschszI3wQ6YceLR11XBkREVHjkJRTjPVnE7Hr+gOUKv6Nre29bTGtiw86eDe8HUc1Fs4rI5VKUVJSAhubhr82JcM51cTxuIf44O9byCspBQC0cbPGp4MC4GFrpuPKiIiIGp+0vBJsOJeEHddSVfuGAEBbDxtM6+KDcF/7BhPStRrOGxOGc1KHQhCw9nQC1vxzX9X2XJgHXn2qBUyMOIyFiIhImzILJNh4Pgl/XElFSem/i5W0cbPG1M4+eKqlg96HdI2Ec39/f7WeqK52CNUEhnN6kkJpKT7YcwtHYh8CUK7G8l5/Mfr6Oeu4MiIioqYlu0iKTReSEXUpBUUyuaq9tbMlpnT2Qa/WTjA00M+QrpFw/s4771QI53K5HJmZmbh06RLc3NwwatQoTJ06tW7V6hDDOVUnKacYb2y/jviHyq8RH3tzfD60DVo4Wuq4MiIioqYrt1iGXy8m49dLySiQ/BvSfezNMbGDFwa0cdW7d7a1PqwlLS0NY8eOxUsvvYQxY8bU9e50huGcqnLmXjYW7I5RjS8P97XH4oEB3FSIiIhITxRISrHlUgo2X0hC7qPf1wDgbGWCce29MLytGyxN9OP3dr2MOf/+++/x+++/Y+/evZq4O51gOKf/EgQBmy8kY+WxeJRNEJ/U0Rszu/vq7VtlRERETVmxTI7t19Lw87lEpBdIVe02ZkYYHeqB58I8YWdhrMMKqw/nGvvzwczMDKmpqZq6OyKdK5UrsORgLHZcSwMAmBoZ4P1+YvQPcNFxZURERFQVc2NDjG3niWdD3PF3TDo2nEvEvaxi5JWU4ofTCfj5fBKGtXXH+PaecLPRvxXW6txzLpPJEBMTgzlz5sDKygo7duzQVG31jj3nVKZAUop5f97A2YQcAICrtSmWDQ2En2vT2DqYiIiosVAIAo7EPsS6MwmIeVCgajcyEGFyJ2+81M233muqt9VavvrqK/Tr169m1ekRhnMCgNS8Erz+R7Rq4meAqxW+HBYIJytTHVdGREREtSUIAs4l5GDd2USce9T5BgCnXu9e7zt6a2RYy7BhwyoN5wYGBnBycsKAAQPg5+dXuwqJ9MT1tHzM3RaNrCIZACCipSM+HugPc2NDHVdGREREdSESidCpmT06NbPH9bR87IpOQwsny3oP5k/CTYgew57zpu1obCbe3X0TkkcbGoxr74lXn2rBiZ9ERESkUfUyIZSoIdt+NRWfHrgDhQAYiIA3e7XCqFAPXZdFRERETUyV4VzdHUH/qyHvEEpNjyAIWHc2Ef87cQ8AYGZkgE8HB6B7C0fdFkZERERNUpXhvLIx5gcOHIBEIkH37t3RokULKBQKJCYm4ujRo7CyssKoUaO0XjCRpigEAcuPxOPXi8kAAFszIywfHoRgDxsdV0ZERERNVZXhfMmSJeU+37hxIw4fPowdO3agefPm5Y4lJSVh3LhxteppJ9IFmVyBD/++hb03MwAALlYmWPVsMFo4Wuq4MiIiImrK1J6e+sMPP2Dy5MkVgjkAeHl5YcKECYiKitJocUTaUCKT443t11XB3NfBHGvHhjKYExERkc6pPSE0Pz8fJiYmVR5XKBSQSqVVHifSB0VSOeZsi8bFpFwAQJC7NZYPC9L5Nr5EREREQA16zkNDQ7Fx40Y8ePCgwrHY2FisW7cOnTp10mhxRJpUICnFK1uvqYJ5Jx87rH62LYM5ERER6Q211zmPjo7GxIkTYWBggIiICHh7e0MqleLu3bs4ceIErK2t8euvv6JZs2barllruM5545VbLMMrW6+ptu3t1twBnw1pA1Mj/dp4gIiIiBq/6tY5r9EmRLGxsVi5ciVOnDiBoiJliLWyskLfvn3x2muvwc3Nre7V6hDDeeOUXSTF7N+v4XZGIQCgZytHfDIoAMZ6tiMYERERNQ0aC+dlBEFAdnY2RCIR7O3t61ScPmE4b3wyC6WYFXUV8Q+V17WP2BkfD/DTu616iYiIqOnQ+A6hIpEIDg4OtS6IqD5kFUkxc8tV3M1SBvMBbVzwfn8/GBlwyU8iIiLST1WG84CAACxduhSDBw8GoN6OoSKRCDdu3NBshUS1kFMsw6yoa6pgPjTIDfP7toYhgzkRERHpsWp3CPXx8Sn3OTcZooYgv6QUr/x+DbGZyjHmgwJdsaBfaxjw65eIiIj0XK3GnDdWHHPe8BVKSzH792uITs0HAPT3d8aHz/izx5yIiIj0hsbHnJeRyWQ4efIkDAwM0LVrVxgZ1enuiOqkWCbH639Eq4J5r9ZO+IDBnIiIiBoQtdO0VCrFokWLkJSUhB9//BFSqRTPPfccbt68CQBo2bIl1q9fD0dHR60VS1QVmVyBt3fcwOXkPABAjxYOWDTQn5M/iYiIqEFRez25r7/+Glu2bIG7uzsAYPv27YiJicHEiRPxySefICMjA1999ZXWCiWqikIQ8MGeWzh9PxsA0MXXHksGt+E65kRERNTgqN1zvmfPHjz77LNYtGgRAGDv3r2wtrbG22+/DSMjIyQmJiIqKkprhRJVRhAEfHEoDvtuZQAAgt1tsHRIG5hw508iIiJqgNROMGlpaQgNDQUAFBcX49y5cwgPD1eNM3d3d0deXp5WiiSqyg+nExB1OQUA0MLRAsuHB8Lc2FDHVRERERHVjtrh3MnJCZmZmQCA48ePQyqVomfPnqrjt27dgouLi8YLBICCggIsWrQITz31FIKCgtCrVy8sX74cUqm0wrmbN29GZGQk2rZti8GDB2P37t1aqYl07/fLKVjzz30AgJu1KVaNDIatubGOqyIiIiKqPbWHtXTu3Bnr16+HqakpNm3aBHNzc/Tp0wd5eXnYunUrtmzZgjFjxmi8QEEQMHv2bJw9exajR4+Gn58fLl++jO+++w6xsbFYvXq16ty1a9di6dKliIyMxOTJk7F//37MnTsXIpEIAwYM0HhtpDuHbmdg6cFYAICduTFWPRsMF2tTHVdFREREVDdqr3Oel5eH1157DadOnYKFhQU++ugjDBo0CBcuXMD48ePRsWNHrF69GjY2Nhot8MCBA5g1axb+7//+D+PHj1e1L1++HN9++y02b96M9u3bIy8vDxEREejVqxeWLVsGAFAoFJg4cSKSk5Nx8OBBGBpWP9yB65w3DNdS8vBy1FVIShWwMDbEN6Pboo1b1euFEhEREemT6tY5V3tYi42NDX766Sf8888/OH36NAYNGgQACAgIwG+//YaNGzdqPJgDwNmzZwEodyh93DPPPAMAuHTpEgDg0KFDKCoqwtixY1XnGBgYYNy4cUhNTVWdRw1bUk4x3th+HZJSBQxFwJIhAQzmRERE1GjUeEkLBwcHZGdn48qVK8jPz4eRkRGCg4O1URsAYObMmdi+fTssLS3LtWdnK5fNK5uQGh0dDQAIDAwsd16bNm3KHaeGK7dYhtf/iEZ2sQwA8Haf1gj3ddBxVURERESaU6NwfuHCBYwYMQI9e/bEmDFjEB0djbNnz6Jnz57466+/tFKgnZ0dAgICKrT/8ssvAICwsDAAQHp6OmxtbWFubl7uPGdnZwBASkqKVuqj+iEtVeDtP2/gfnYxAGBSR2+MaOuu46qIiIiINEvtCaFXr17FlClT4O7ujueffx7r168HANja2sLIyAhvvvkmLC0tERERodb9ZWRkVHvcwsKiQm95me3bt2Pv3r3o3LkzQkJCAACFhYUwMzOrcG5ZW3FxsVp1kf4RBAGL9t3GxaRcAEAfsTNm9fDVbVFEREREWqB2OP/qq6/g5eWFP/74A0VFRVi3bh0AIDg4GH/++SfGjh2L7777Tu1w3r1792qPz5gxA3PmzKnQfvDgQbz33ntwdnbGkiVLVO0KhQIiUcWt2svaKjtGDcO6s4nYE5MOQLnJ0MJIMQx4PYmIiKgRUjucX7p0CTNnzoSZmVmFXmgrKyuMHj0aK1euVPuBy3YarUplQ1l27dqFd955B5aWlvj+++/h4eGhOmZpaYmSkpIKtymrtapeeNJvx+Ie4psT9wAAHrZmWDasDcy4yRARERE1UmqHcwAwMTGp8phEIoFCoVD7vkaNGlWTh8avv/6KDz/8ELa2tli3bh38/f3LHXd3d0dubi6kUmm5OtPTlT2urq6uNXo80r27D4vwf3/dhADA3NgAy4YFwt6i6q9BIiIiooZO7QmhISEh2LVrV6XHioqKEBUVpbVVW7Zv346FCxfCyckJmzZtqhDMAeUqLYIgICYmplx72efaXFGGNC+vRIY3d1xHoVQOAPjwGX+0cuK7H0RERNS4qR3OX331Vdy4cQMTJkzA9u3bIRKJcPXqVWzYsAFDhw5FUlISZsyYofEC4+Li8P7778PBwQEbN25Ey5YtKz0vIiICpqam2Lhxo6pNoVBg8+bN8PT0RGhoqMZrI+2QKwS8u/smEh6tzDI9vBmebu2k46qIiIiItE/tYS1hYWH47rvvsHDhQnz22WcAlLt0AsrlCpcvX44uXbpovMBVq1ZBKpWiR48euHLlCq5cuVLuuJ+fH/z9/WFvb4/p06dj1apVEAQBXbp0wd69e3H+/HmsWLHiibuDkv5YffwuTt9TrmPfs5UjpoX76LgiIiIiovohEgRBqMkNBEHAjRs3kJCQAIVCAU9PTwQFBak2A9K0bt26ITMzs8rjs2fPxiuvvKKq7ccff8SmTZuQmZkJX19fzJw5E5GRkWo9lkwmR05OkUbqpto5cCsD83cphyK1dLLA2rGhsDTRztcWERERkS44O1e9u3mNw3lV5HI5vvvuO8ycOVMTd6cTDOe6dT+rCM9vuoRCqRzWpkbYMCEMXnbmT74hERERUQNSXTivtksyLy8PUVFRuHLlCgRBQJs2bTB+/HjY2NiUO+/KlSt47733EBsb26DDOelOiUyO+btiHpsA6sdgTkRERE1OleE8ISEBEydORHp6Oso61/fv349NmzZhy5Yt8PDwgFQqxRdffIFNmzZBLpdj8ODB9VY4NS6fH4rFnYxCAMCkjt7o0dJRxxURERER1b8qV2v58ssvkZ6ejjfeeAMnT57ExYsXsWLFCigUCnz00UfIy8vDhAkTsGHDBnh4eGDt2rX4/PPP67N2aiR2Rqfhz+gHAIAwL1u83N1XtwURERER6UiVPecXLlzA8OHD8cILL6jaIiMjIZFI8O677+Ktt97CtWvXMGnSJLzxxhswNTWtl4KpcYnNKMRnB2MBAA4WxvhkoD+MDEQ6roqIiIhIN6oM59nZ2QgLC6vQ3qFDB5SWluLMmTNYvXo1evXqpdUCqfEqksoxb+cNSEoVMBABiwcGwMmKf+QRERFR01XlsJbS0lKYmZlVaLewsAAATJs2jcGc6uSLQ7GqjYZe6uqLDj52ui2IiIiISMfU3iH0v9q1a6fJOqiJ2X8rAzuvK8eZd/Sxw+TO3jquiIiIiEj3ah3ODQxqfVNq4lLzSvDJ/tsAAFszI3wQ6QcDEceZExEREVW7znl8fDzOnTtXri0/Px8AcOvWrUp3Be3YsaMGy6PGRq4QsPCvmyiQKNczf6+fGC7WHGdOREREBFSzQ6i/vz9EVfRmCoJQ5bGYmBjNVVfPuEOo9q09fR/fnrwPABjR1h3z+7bWcUVERERE9atWO4TOnj1bK8VQ0xWdmofv/1EGc18Hc8zp2ULHFRERERHplyp7zpsi9pxrT4lMjgkbL+J+djGMDERYNy4Mfq5Wui6LiIiIqN5V13POWZ1UL745eQ/3VcsmNmMwJyIiIqoEwzlp3cWkHPxyIRkAEORujQkduWwiERERUWUYzkmriqRyfPj3bQgATI0MsDDSD0YGXDaRiIiIqDIM56RVK4/FIyW3BAAws7svfB0sdFwRERERkf5iOCetOXM/G1uvpAIAwrxsMaadp44rIiIiItJvaofzSZMm4dSpU1UeP3ToEAYOHKiRoqjhK5bJ8ck+5S6g5sYG+L/+Yu4CSkRERPQEVa5zXlxcjOzsbNXnZ8+eRd++fdGsWbMK5yoUChw7dgxJSUnaqZIanDX/3EdKngQAMKt7c3jZmeu4IiIiIiL9V+U651lZWYiMjER+fr5adyQIArp164a1a9dqtMD6xHXONePmg3w8v+kSFIJydZYfxoTCkJNAiYiIiADUcodQBwcHfP7557h27RoEQcDq1avRt29f+Pn5VTjXwMAADg4OHNZCKFUIWLTvDhQCYGggwrt9xQzmRERERGqqMpwDQEREBCIiIgAAKSkpGDNmDEJCQuqlMGqYfrmQhFvpBQCA5zt6oZWzpY4rIiIiImo4qhzW0hRxWEvdJOcW47l1FyApVcDH3hybJ7WHqREXBCIiIiJ6XK2GtVTm2LFj2LlzJzIzMyGXyyscF4lEWL9+fc0rpAZPEAQsORALSakCALCgb2sGcyIiIqIaUjucb9q0CYsWLQIAODo6wsTERGtFUcNz+E4mTt9Tru4zNNgN7b3tdFsQERERUQOkdjjfsGED/P398f3338PJyUmbNVEDUyyT48sj8QAAWzMjvNKjuY4rIiIiImqY1B53kJqaiueee47BnCr46UwCHuQ/WtO8R3PYmhvruCIiIiKihkntcO7j44PMzExt1kINUEJ2MX4+r9x8qo2bNYYGu+m4IiIiIqKGS+1wPn36dGzcuBF37tzRZj3UgAiCgC8OxUImFyAC8HbvVjAQcU1zIiIiotpSe8z5hQsXYGlpiaFDh6J58+ZwcHCA6D9BjKu1NC1HYx/i1GOTQAPdql4WiIiIiIieTO1wfvz4cQCAm5sbiouLkZycrLWiSP+VyOT48kgcAOUk0FndOQmUiIiIqK7UDueHDh3SZh3UwPxyMRmpecpJoC9394WdBSeBEhEREdVVrXaJefDgAa5cuYL8/HxIpVIoFApN10V6LLNQinVnEgEArZwsMSzYXccVERERETUONQrnFy5cwIgRI9CzZ0+MGTMG0dHROHv2LHr27Im//vpLWzWioKAAn3zyCXr16oXQ0FCMGTMGJ06cqPTczZs3IzIyEm3btsXgwYOxe/durdXVVK355x6KZModYl+PaAFDA04CJSIiItIEtcP51atXMWXKFBQWFuL5559Xtdva2sLIyAhvvvkmjh49qvECFQoFZs+ejc2bN6N///546623AAAvvPACjh07Vu7ctWvX4sMPP4Sfnx8WLFgAFxcXzJ07V6t/ODQ1sZmF2HEtDQDQtbk9Ovva67giIiIiosZDJAiCoM6J06ZNQ2pqKv744w8UFRWha9eu+OmnnxAeHo6CggKMHTsW1tbW2Lx5s0YLPHToEF5++WUsXLgQ48aNAwBIJBJERkbCyckJUVFRAIC8vDxERESgV69eWLZsGQBlsJ84cSKSk5Nx8OBBGBoaVvtYMpkcOTlFGq2/sXll6zWcvpcNQxGw+fn2aOFoqeuSiIiIiBoUZ+eqV7hTu+f80qVLGDFiBMzMzCosoWhlZYXRo0drZQ30vLw8+Pv7Y8iQIao2U1NTBAcH4/bt26q2Q4cOoaioCGPHjlW1GRgYYNy4cUhNTcWlS5c0XltT88/dLJx+tHTisLbuDOZEREREGlajMecmJiZVHpNIJFqZGDps2DDs2LEDVlZWqja5XI47d+7Aze3f3Sijo6MBAIGBgeVu36ZNm3LHqXZKFQK+OhoPALA0McT0rs10XBERERFR46N2OA8JCcGuXbsqPVZUVISoqCgEBwdrrLCqHufatWt4/fXXER8fj5kzZ6qOpaenw9bWFubm5uVu4+zsDABISUnRam2N3a7oNMQ/VA75mdzJGw4WVf+hRkRERES1o/Y656+++iomTpyICRMmoHfv3hCJRLh69Sru3LmDjRs3IiUlBR9++KHaD5yRkVHtcQsLC1halh828eWXX2Ljxo0AgP79+6Nfv36qY4WFhTAzM6twP2VtxcXFatdG5UlKFfj+1H0AgKu1Kca089RxRURERESNk9rhPCwsDN999x0WLlyIzz77DACwfPlyAMre6eXLl6NLly5qP3D37t2rPT5jxgzMmTOnXFvv3r3RpUsXXLlyBT/99BMmTZqETZs2wcTEBAqFosJYeACqtsqOkXq2XklBeoEUADA9vBnMjKufWEtEREREtaN2OAeAbt26Yf/+/bhx4wYSEhKgUCjg6emJoKAgGBnV6K6waNGiao8HBARUaAsPDwcA9OnTB56enli4cCH+/PNPPPvss7C0tERJSUmF25T1mP+3F57UUygtxU+PNhzysTfHgEBXHVdERERE1HjVKFGnpKRg8+bNePHFF1UTL9esWYO9e/fihRdegKOjo9r3NWrUqJpV+h/PPPMMFi5ciBs3bgAA3N3dkZubC6lUWm7ianp6OgDA1ZWhsjZ+uZCMnGIZAGBGN18YccMhIiIiIq1Re0Lo7du3MXz4cPz0009ITU1Vtefl5WHTpk0YNmwYEhMTNV7gxx9/jL59+1ZYCaawsBCAcllFQLlKiyAIiImJKXde2efanqzaGOUUy/Dz+SQAgNjZEr3FTjquiIiIiKhxUzucL1u2DJaWlti9ezf8/f1V7W+++SZ2794NY2NjfPHFFxov0MPDAwkJCdizZ0+59h9//BEA8PTTTwMAIiIiYGpqqpowCig3Idq8eTM8PT0RGhqq8doauw1nE1EolQMAZvZoDgOO2yciIiLSKrWHtVy+fBkzZ86Er69vhWPe3t6YMGECfvjhB03WBgCYMGECtm/fjgULFiA6Ohre3t44ceIEDh48iGeffRadOnUCANjb22P69OlYtWoVBEFAly5dsHfvXpw/fx4rVqx44u6gVF5GgQRbLiuXnwz1tEFXX3sdV0RERETU+KkdzgVBgEQiqfZ4ZRMy68rU1BTr16/HF198ge3btyM/Px8+Pj549913MXHixHLnzpo1C+bm5ti0aRP2798PX19ffPXVV4iMjNR4XY3dj6cTIClVDiWa2b05V7shIiIiqgciQRAEdU588cUXER8fj23btsHGxqbcscLCQowcORLu7u746aeftFJofZDJ5MjJKdJ1GTr3IF+C4WvPQiYX0MXXHqtGcrw+ERERkaY4O1tXeUztnvPZs2djwoQJGDRoEAYPHoxmzZpBJBIhISEBu3fvRkZGBj799FONFEy6tfFcImRy5d9sL4Y303E1RERERE2H2j3nAHD+/Hl89tlniI6OxuM38/f3x/z589G5c2etFFlf2HMOZBZIMGztOUhKFejczA5fP9tW1yURERERNSoa6TnPzs5Ghw4dEBUVhaysLCQnJ0OhUMDd3R0uLi4aKZR0b+P5JNVY8xe6sNeciIiIqD6pHc6HDx+OUaNGYdasWXBwcICDg4M26yIdyCqSYusV5Rr2HbxtEeplq+OKiIiIiJoWtdc5z8rKgrOzszZrIR3b9HivOceaExEREdU7tcP54MGD8dtvvyEpKUmb9ZCO5BTJEPVoXfMwTxu0Y685ERERUb1Te1iLgYEB4uPj0b9/f/j4+MDR0REGBuWzvUgkwvr16zVeJGnfb5eSUSxT9ppPC2/Gdc2JiIiIdEDtcH7y5EnY2yt3iZRIJEhJSdFaUVS/imVyVa95GzdrdPKx021BRERERE2U2uH80KFD2qyDdGjHtTTklpQCAJ7v6MVecyIiIiIdUXvM+eMePHiAK1euID8/H1KpFAqFQtN1UT0plSuw6bxyHoGPvTkiWjnpuCIiIiKipqtG4fzChQsYMWIEevbsiTFjxiA6Ohpnz55Fz5498ddff2mrRtKi/bczkJYvAQCM7+AFQwP2mhMRERHpitrh/OrVq5gyZQoKCwvx/PPPq9ptbW1hZGSEN998E0ePHtVKkaQdgiBg4zllr7mDhTEGtnHVcUVERERETZva4fyrr76Cl5cXduzYgenTp0MQBABAcHAw/vzzT7Rs2RLfffed1golzTt1Lxt3MgoBAGPaecLUqFajnIiIiIhIQ9ROY5cuXcKIESNgZmZWYcKglZUVRo8ejTt37mi8QNKeDecSAQAWxoYYGeKu42qIiIiIqEZdpSYmJlUek0gknBjagNx6UIALibkAgGFt3WBjZqzjioiIiIhI7XAeEhKCXbt2VXqsqKgIUVFRCA4O1lhhpF2/XEoGABiKgLHtPHVcDREREREBNQjnr776Km7cuIEJEyZg+/btEIlEuHr1KjZs2IChQ4ciKSkJM2bM0GatpCEPC6XYdzMdANCztRPcbMx0XBERERERAYBIKJvZqYaTJ09i4cKFSEpKKtfu7OyM999/H/369dN4gfVJJpMjJ6dI12Vo3ff/3MeaU/eVHz8XglAvWx1XRERERNR0ODtbV3msRuEcUC6/d+PGDSQkJEChUMDT0xNBQUEwMlJ7s1G91RTCubRUgcHfn0FWkQz+LlbYMCGMO4ISERER1aPqwvkTE7VMJkNsbCxKS0vRqlUrmJubIzAwEIGBgRotkurHgdsZyCqSAQDGtvdkMCciIiLSI9WG83Xr1mH16tUoKCgAoFytZdy4cXjjjTcaRU95UyMIAn65oJwI6mBhjD5iZx1XRERERESPqzJhb9++HUuWLIGnpyeGDh0KAwMDnDlzBuvWrYNcLseCBQvqs07SgCvJebiZrvxD69kQD5hw0yEiIiIivVJlON+8eTNCQ0Oxfv16mJqaAlD2vM6ZMwe//fYb3nzzzWrXPSf989ulFACAsaEII7jpEBEREZHeqbLrNC4uDoMHD1YFcwAQiUSYPHkypFIp4uPj66VA0ozMQikOx2YCAPr6OcPRkn9YEREREembKsN5cXExrK0rziT18vKCIAjIy8vTamGkWTuj0yBXKBfmGRnioeNqiIiIiKgyVYZzhUJR6UoehoaGAAC5XK69qkij5AoB266mAgBaO1si2L3q5XuIiIiISHc4I7AJOH0vG6l5EgDAiLbuXD6RiIiISE9Vux5iTk4OUlJSyrXl5uYCALKysiocAwAPDw6Z0Ddbryivk7mxASIDXHRcDRERERFVpcodQv39/avsYRUEodJjIpEIN27c0GyF9agx7hCalleCoT+chUIAhrd1w4K+Yl2XRERERNSk1WqH0OHDh2ulGKpf26+l4dE8UIxsy3c1iIiIiPRZleH8008/rc86SAtK5QrsuJYGAAhyt4afq5WOKyIiIiKi6nBCaCN2PD4LmYVSAMqJoERERESk3xjOG7E/o5W95pYmhujr56zjaoiIiIjoSRjOG6nMAgn+uZsFAOjv7wIzY0MdV0RERERET9Igw/nrr78OPz+/So9t3rwZkZGRaNu2LQYPHozdu3fXc3X64a8b6aqJoIODXHVbDBERERGppcGF83379mHPnj2VHlu7di0+/PBD+Pn5YcGCBXBxccHcuXPx119/1XOVuiUIgmpIS3NHCwS6cUdQIiIiooagynXO9VF2djYGDhyIvLw8yGQy3Lp1S3UsLy8PERER6NWrF5YtWwYAUCgUmDhxIpKTk3Hw4EEYGlY/tKOxrHN+NSUP0365DAB49anmmNjRW7cFEREREZFKdeucN6ie848//hj29vbo169fhWOHDh1CUVERxo4dq2ozMDDAuHHjkJqaikuXLtVnqTpV1mtuKAKeacMhLUREREQNRYMJ5wcOHMCePXuwePFimJiYVDgeHR0NAAgMDCzX3qZNm3LHG7timRwHbmUAALq1cISTZcXXioiIiIj0U4MI5zk5Ofjggw/w/PPPIzQ0tNJz0tPTYWtrC3Nz83Ltzs7KJQRTUlK0XaZeOHQ7E4VSOQBgcCB7zYmIiIgakip3CNW2jIyMao9bWFjA0tISgHI4i7m5OV5//fUqzy8sLISZmVmF9rK24uLi2hfbgOy8rhzS4mBhjO4tHHRcDRERERHVhM7Ceffu3as9PmPGDMyZMwcHDx7E7t27sWHDhkrDdxmFQgGRSFShvaytsmONTWpeCS4k5gIAIgNcYGTYIN4YISIiIqJHdBbOFy1aVO3xgIAA5ObmYuHChRg0aBBatWqFrCzlpjpSqXJL+qysLJiamsLS0hKWlpYoKSmpcD9lPeZlvfCN2d6YdNXHAzgRlIiIiKjB0Vk4HzVq1BPPOXPmDDIyMrBz507s3LmzwvHw8HAMHz4cS5Ysgbu7O3JzcyGVSstNGE1PVwZWV9fGH1b/vql8rs0dLSB2bvx/jBARERE1NjoL5+rw9/fHTz/9VKF97dq1OHHiBH766Se4uLgAUK7SIggCYmJiEBISojo3JiYGABAcHFw/RevInYwCxGUq12h/JsClSQzjISIiImps9Dqc29raomvXrhXa//zzTwAodywiIgKmpqbYuHGjKpwrFAps3rwZnp6eVa7y0lj8/diQlv7+LjqshIiIiIhqS6/DeU3Y29tj+vTpWLVqFQRBQJcuXbB3716cP38eK1aseOLuoA2ZQhCw96Zy9ZsQDxt42FY9cZaIiIiI9FejCecAMGvWLJibm2PTpk3Yv38/fH198dVXXyEyMlLXpWnVpaRcPMiXAFCu0kJEREREDZNIEARB10XoC5lMjpycIl2XUWOL993G9mtpMDQQ4e+XusDOwljXJRERERFRFZydras8xoWwGzhpqQIHb2cCAMJ97RnMiYiIiBowhvMG7tS9LORLSgEoV2khIiIiooaL4byB+ztGORHU3NgAT7V01HE1RERERFQXDOcNWIlMjhPxDwEAT7V0hJlx412RhoiIiKgpYDhvwP65l42SUgUAoI/YWcfVEBEREVFdMZw3YAdvKYe0WBgboouvvY6rISIiIqK6YjhvoJRDWrIAAN1bOHBICxEREVEjwHDeQJ2+l40imRwA0FvspONqiIiIiEgTGM4bqIN3lGubmxkZoGtzBx1XQ0RERESawHDeAElKFTgep1ylpXsLrtJCRERE1FgwnDdAp+9lo1DKIS1EREREjQ3DeQN06I5ylRZTIwN0a8EhLURERESNBcN5AyMtVeBorHJIS7fmDjDnkBYiIiKiRoPhvIE5l5DDIS1EREREjRTDeQNzNE65SouxoYhDWoiIiIgaGYbzBkQhCDgWp9x4qKOPHSxNjHRcERERERFpEsN5AxKTlo+HhVIAwFMtHXVcDRERERFpGsN5A3L00drmANCjBcM5ERERUWPDcN6AHHsUzgNcreBibarjaoiIiIhI0xjOG4iknGLEZRYBACJasdeciIiIqDFiOG8gjj02pIXjzYmIiIgaJ4bzBqIsnHvYmKKVk6WOqyEiIiIibWA4bwByi2W4nJQLAOjR0hEikUjHFRERERGRNjCcNwAn72ZBLig/5nhzIiIiosaL4bwBKBvSYm1qhDBPWx1XQ0RERETawnCu52RyBU7fywYAdG1uDyNDXjIiIiKixopJT89dTs5FoVQOAOjOjYeIiIiIGjWGcz33z11lr7mBCOjia6/jaoiIiIhImxjO9dw/d7MAAIFu1rAzN9ZxNURERESkTQzneiwtrwTxD5W7goY3d9BxNURERESkbQzneuyfRxNBAaAbwzkRERFRo8dwrsf+iVcOaXGwMIa/q5WOqyEiIiIibWM411MyuQLnEnIAKCeCGnBXUCIiIqJGr0GE86NHj8LPz6/Sf7dv3y537ubNmxEZGYm2bdti8ODB2L17t46qrpvLybkokimXUOzqyyEtRERERE2Bka4LUMedO3cgEonw6aefwsCg/N8T7u7uqo/Xrl2LpUuXIjIyEpMnT8b+/fsxd+5ciEQiDBgwoL7LrpPHl1DszCUUiYiIiJqEBhPOPTw8MHz48CrPycvLw9dff41BgwZh2bJlAIDRo0dj4sSJWLp0Kfr37w9DQ8P6KrnO/l1C0YZLKBIRERE1EQ1iWMudO3fQokWLas85dOgQioqKMHbsWFWbgYEBxo0bh9TUVFy6dEnbZWrMg3zJY0sosteciIiIqKnQ+3AuCALi4+PRqlUrAIBEIkFpaWmF86KjowEAgYGB5drbtGlT7nhDcC7h3yUUuzRjOCciIiJqKvQ+nCcmJqK4uBgJCQkYPHgwQkJCEBoaijlz5iArK0t1Xnp6OmxtbWFubl7u9s7OzgCAlJSUeq27Ls7czwEAWJkaIsDNWrfFEBEREVG90dmY84yMjGqPW1hYwNLSEnfu3AEAXLlyBS+++CK8vLxw4cIFbNiwAbGxsYiKioKZmRkKCwthZmZW4X7K2oqLizX/JLRAEAScva/sOe/gbQcjAy6hSERERNRU6Cycd+/evdrjM2bMwJw5c+Dt7Y1Zs2Zh8ODBaN68OQCgT58+aNasGRYuXIjff/8dEyZMgEKhgKiStcDL2io7po/iMouQVSQDAHT04ZAWIiIioqZEZ+F80aJF1R4PCAgAAIjFYojF4grHR44ciY8//hhnzpzBhAkTYGlpiZKSkgrnlfWYW1paaqBq7Tv72Hjzzs3sdFcIEREREdU7nYXzUaNG1en2xsbGsLGxQVGRclUTd3d35ObmQiqVwsTERHVeeno6AMDV1bVOj1dfzj4ab+5qbQofe/PqTyYiIiKiRkXvJ4SuWLECvXv3RkFBQbn2nJwcZGVlwdPTE4BylRZBEBATE1PuvLLPg4OD66fgOpDJFbiYlANA2WveUIbiEBEREZFm6H049/DwQFJSEn7//fdy7atXrwYADB48GAAQEREBU1NTbNy4UXWOQqHA5s2b4enpidDQ0HqrubaupeahWKYAAHTieHMiIiKiJkfvdwgdPnw4oqKi8Pnnn+PevXsQi8U4deoU9u3bh+eeew4dO3YEANjb22P69OlYtWoVBEFAly5dsHfvXpw/fx4rVqxoELuDlg1pAYCOHG9ORERE1OSIBEEQdF3Ek+Tk5ODLL7/EwYMHkZubC29vbzz33HOYNGkSDAz+7fwXBAE//vgjNm3ahMzMTPj6+mLmzJmIjIxU63FkMjlycoq09TSeaOrmS7iWmo/WzpbYPKm9zuogIiIiIu1xdq56H5sGEc7riy7DeYGkFL1X/wOFAIxv74XXe7bQSR1EREREpF3VhXO9H3PeVFxOzoXi0Z9JHX3sdFoLEREREekGw7meuJCYCwAwFAEhnjY6roaIiIiIdIHhXE9cTFKGc39Xa1iZ6v08XSIiIiLSAoZzPVAgKcXNB/kAgHZetjquhoiIiIh0heFcD1xJzlONN2/vbafTWoiIiIhIdxjO9cCFxBwAgAHHmxMRERE1aQzneuACx5sTERERERjOde7x8ebtOd6ciIiIqEljONcxjjcnIiIiojIM5zrG8eZEREREVIbhXMc43pyIiIiIyjCc65CkVMHx5kRERESkwnCuQyaGIrjbmMHIQIT+AS66LoeIiIiIdEwkCIKg6yL0hUwmR05OUb0+ZoGkFDK5AvYWJvX6uERERESkG87O1lUe4yBnHeM4cyIiIiIqw2EtRERERER6guGciIiIiEhPMJwTEREREekJhnMiIiIiIj3BcE5EREREpCcYzomIiIiI9ATDORERERGRnmA4JyIiIiLSEwznRERERER6guGciIiIiEhPMJwTEREREekJhnMiIiIiIj0hEgRB0HURRERERETEnnMiIiIiIr3BcE5EREREpCcYzomIiIiI9ATDORERERGRnmA4JyIiIiLSEwznRERERER6guFcixITEzF79mx06tQJnTp1wttvv42srCyt3Y50QxPX67333sPEiRO1VCFpQm2v8/HjxzFu3DiEhIQgLCwMkydPxuXLl7VfMNVKba/zqVOnMHbsWISFhaFHjx5YvHgxCgsL66FiqilN/My+efMmgoKCsGrVKi1VSXVV2+v87LPPws/Pr8K/V199tR6qVjKqt0dqYrKzs/H8889DKpXihRdegFwux9q1a3Hr1i1ERUXBxMREo7cj3dDE9YqKikJUVBQ6depUDxVTbdT2Op89exYvvvgiWrdujTlz5qC0tBSbN2/GhAkTsHnzZrRt27aenwlVp7bX+fTp05g6dSoCAwPx5ptvIjU1FRs2bEB0dDQ2bdoEAwP2g+kLTfzMLi0txfz58yGTyeqhYqqN2l5nQRAQFxeHPn36oF+/fuWOeXp61kfpqkJIC7788kshICBAiI2NVbWdPHlSEIvFwm+//abx25Fu1OV6lZaWCqtWrRL8/PwEsVgsTJgwQdvlUi3V9joPHTpU6Nmzp1BUVKRqy8jIEDp27ChMnjxZqzVTzdX2Og8fPlx4+umnheLiYlXbzz//LIjFYuHIkSNarZlqRhO/Y7/++mshMDBQEIvFwsqVK7VVKtVBba9zQkKCIBaLha1bt9ZHmVXin/Nasnv3bnTq1AktW7ZUtXXt2hXNmzfH7t27NX470o3aXi+JRILhw4dj1apVGDp0KFxdXeujXKql2lzn3Nxc3Lx5E5GRkTA3N1e1Ozk5oWPHjrh06ZLW66aaqc11lkgksLe3x+jRo2FmZqZqL3sn7NatW9otmmqkrr9jb926hW+++QYzZ87UZplUR7W9zrGxsQBQ7na6wHCuBbm5uUhMTERgYGCFY4GBgYiOjtbo7Ug36nK9JBIJCgoKsHz5cnz22WcwMuIIM31V2+tsZWWFv//+G5MnT65wLDs7G4aGhpouleqgttfZ1NQUa9euxYwZM8q1x8TEAAA8PDw0XyzVSl1/x5YNZ+natSuGDBmirTKpjupyne/cuQPg33BeVFSknSKfgOFcCx48eAAAlfaGOjs7o6CgAPn5+Rq7HelGXa6XlZUV9u3bhwEDBmi1Rqq72l5nQ0ND+Pr6VrjdzZs3cfHiRYSFhWmnYKoVTf38TU5Oxh9//IHFixdDLBajb9++Gq+Vaqeu1/j777/H/fv38dFHH2mtRqq7ulznO3fuwNLSEp9++inCwsIQFhaGPn361PvIBXbXaUHZDP3H38ouY2pqCkD515i1tbVGbke6UZfrZWBgwEliDYQmvy8LCwsxb948AMD06dM1WCXVlSauc05ODnr16qW6n/fee091W9K9ulzjO3fuYPXq1fi///s/uLm5ISkpSbvFUq3V5TrHxsaisLAQ+fn5WLp0KfLy8rBhwwbMnTsXMpkMw4YN02rtZRjOtUChUDzxnMqCWW1vR7rB69U0aOo6FxcX4+WXX8bNmzfx0ksvcXUePaOJ6ywSibB8+XJIpVJs3LgRU6ZMwZdffonIyEhNlUl1UNtrLJfLMX/+fLRv3x6jR4/WRmmkQXX5Xh49ejQUCgXGjx+vahs4cCAGDRqEzz//HIMHD66XIYlMDlpgaWkJQDmu+L/K2srO0cTtSDd4vZoGTVznvLw8TJ06FWfOnMHIkSMxZ84czRdKdaKJ62xra4sBAwZg2LBh2LRpEzw8PLBkyRLNF0u1UttrvHbtWty8eRNvvPEGsrKykJWVhby8PADKP7qzsrLUCoRUP+ryvTx27NhywRwAzMzMMHToUGRmZqomjGobw7kWlE0AysjIqHAsPT0dNjY2sLCw0NjtSDd4vZqGul7nhw8fYtKkSbh48SKee+45LF68GCKRSGv1Uu1o+vvZzMwMPXv2RGpqKjeR0xO1vcbHjx+HTCbDqFGjEB4ejvDwcAwfPhyAMriHh4cjJSVFu8WT2rTxu9nBwQFA/U0Q5bAWLbCxsYGXlxeuX79e4diNGzcQFBSk0duRbvB6NQ11uc4FBQWYNm0aYmJiMHnyZMyfP1+bpVId1PY6x8XF4cUXX8S0adMq9LgVFhZCJBJx8zg9UdtrPG/ePFVPeZnMzEy89dZbGDp0KIYNGwZnZ2et1Ew1V9vr/ODBA0ydOhXPPPMMZs+eXe7Y3bt3AQBeXl6aL7gS7DnXkn79+uHUqVOIi4tTtf3zzz+4e/dutSt01PZ2pBu8Xk1Dba/zRx99hJiYGEyaNInBvAGozXVu1qwZ8vPz8euvv0Iqlarak5OTsW/fPnTs2BFWVlZar53UU5trHBQUhK5du5b7165dOwCAt7c3unbtyom/eqY219nV1RX5+fmIiopCQUGBqj01NRV//PEHOnfuXG9/hIkEQRDq5ZGamKysLAwaNAiGhoaYOnUqJBIJfvjhB/j4+ODXX3+FiYkJEhMTcfHiRbRr1w7e3t5q3470R22v83/16tULnp6e2LhxYz0/A1JHba5zXFwcBgwYAGtrayxYsKDSSURDhw7VwbOhqtT2+3nHjh14++23ERoaiiFDhiA7OxubNm2CTCbD5s2bIRaLdfzMqIymfmYnJSWhd+/emD17Nl555ZV6fhb0JLW9zgcOHMCsWbPQunVrjBo1CoWFharv5V9++aX+NifS6f6kjVxcXJzwwgsvCKGhoUKXLl2EefPmCQ8fPlQd37p1a6XbxD7pdqRfanudH/f0008LEyZMqI9yqZZqep03b94siMXiav+R/qnt9/Pu3buF4cOHC4GBgULHjh2FV155RYiPj6/v8kkNmviZnZiYKIjFYmHlypX1UTLVQm2v8/79+4Vnn31WCAoKEjp06CDMnDlTiI2Nrdfa2XNORERERKQnOOaciIiIiEhPMJwTEREREekJhnMiIiIiIj3BcE5EREREpCcYzomIiIiI9ATDORERERGRnmA4JyIiIiLSE0a6LoCIiGrnnXfewbZt28q1GRsbw8XFBU8//TReffVV2Nraljv31q1buiiViIjUxHBORNTAzZ8/H/b29gAAiUSC2NhY/Pbbb7h27Rp++eUXGBoa4rnnnkN4eLiOKyUioidhOCciauD69OkDLy+vcm2+vr748MMPcezYMTz99NMICwtDWFiYjiokIiJ1ccw5EVEj1LlzZwBAbGysjishIqKaYDgnImqE0tLSAADe3t4AlGPO/fz8yp2TnJyMt956C126dEFwcDCGDBmCLVu2lDvnnXfewaBBg3DhwgU899xzaNu2LXr37o1t27ZBJpNh2bJl6NatGzp16oTXX38d2dnZ5W7/999/Y8KECWjfvj2CgoLQq1cvLF26FFKpVHWOVCrF4sWL0bt3bwQFBSEiIgIffvghcnNzVecIgoCvv/4a/fv3R3BwMLp27Yq33noLqampGn3diIh0jcNaiIgauLy8PGRlZQEAZDIZ4uLisGjRIgQGBqJXr16V3iYxMRGjR4+GRCLBhAkT4OzsjH379uH999/HvXv38Pbbb6vOzcjIwIwZMzBq1CgMGTIEGzZswIIFC7Bz507k5+dj5syZiIuLw6ZNm2Bubo5PP/0UABAVFYX33nsPvXr1wptvvgmZTIb9+/dj7dq1sLCwwOzZswEAH330EXbt2oVJkybB29sbd+7cwaZNm3D//n38+OOPAIBvv/0Wq1evxvjx4+Hn54ekpCRs2LAB0dHR2LVrFwwNDbX5EhMR1RuGcyKiBm748OEV2szMzLBhwwaYmJhUepsvv/wSOTk5+P333xEYGAgAGD9+PGbOnIkff/wRw4cPR+vWrQEAOTk5eP/99zFhwgQAgJeXF6ZPn4579+7h77//Vj1GTEwMTpw4oXqMH3/8EWFhYfjf//4HkUgEABg3bhx69+6NvXv3qsL5zp07MXLkSMydO1d1WwsLCxw/fhyFhYWwtLTEzp078dRTT+G9995TnePu7o5ffvkFycnJ8PHxqfXrR0SkTxjOiYgauM8//xxOTk4AlD3nycnJ2LRpE8aPH4/vvvsO3bp1K3e+XC7HkSNH0L17d1UwBwADAwPMmDEDhw8fxqFDh1ThHAD69u2r+tjX1xcA0KNHj3Lh38vLC5cuXVJ9/ueff6K4uFgVzAHg4cOHsLGxQVFRkarNzc0Nf/31F4KCgtCnTx/Y2Njg9ddfx+uvv17unDNnzmD9+vUYOHAgnJycMGbMGIwZM6aWrxoRkX7imHMiogauXbt26Nq1K7p27YqIiAiMGzcOP//8M8zNzbFo0aIK52dnZ6OoqAjNmzevcKxly5YAlOPRH+fo6Kj6uGwIyeNtZe2CIKg+NzY2RnR0NBYsWIAxY8aga9eueOqpp3D79u1y533wwQcQBAHz589HeHg4xo8fj3Xr1iE/P191zttvvw17e3t88skn6N69O0aOHInVq1cjIyOjJi8VEZHeYzgnImqE7O3t0blzZ8THx5ebWAmgXDD+L4VCAQAVhsMYGVV8o/XxHvHKLFu2DFOmTEFMTAwCAgLwyiuvYMeOHejQoUO588LDw3H48GF8+eWXeOaZZxAfH49PP/0UgwcPVo2l9/f3x969e/G///0Pzz77LDIzM7Fy5UoMGDAAcXFx1dZBRNSQcFgLEVEjVRa0DQzK98M4ODjAwsIC8fHxFW5z9+5dAMphJHWRnJyMNWvWYOjQoVi6dGm5Y5mZmaqPpVIpYmJi4ObmhoEDB2LgwIFQKBT46aefsHTpUuzevRvjxo3DzZs3YWVlhd69e6N3794AgL/++gtz5sxBVFQU3nnnnTrVS0SkL9hzTkTUCGVmZuL06dMICAiAtbV1uWOGhobo0aMHTp48ievXr6vaBUHA999/D5FIhJ49e9bp8ct661u1alWu/ejRo7h37x5KS0sBKCebPvfcc/juu+9U5xgYGCA4OFj1sVwux6RJk/DJJ5+Uu6+QkBDVOUREjQV7zomIGrgDBw7A3t4egDJgp6WlYcuWLSguLsacOXMqvc2bb76JM2fOYOLEiZg4cSKcnZ2xf/9+nD59GlOmTKkQqmuqVatW8PDwwLfffguJRAI3NzdcvXoV27Ztg6mpKQoLCwEALi4uGDx4MDZv3ozi4mKEhYUhJycHP//8M5ycnPDMM8/AxMQEEydOxDfffINZs2ahR48eKCkpwW+//QZzc3OMHDmyTrUSEekThnMiogaubF1xQNkrbmtri+DgYCxevBjh4eGV3sbHxwdbtmzBihUr8Ouvv6KkpAQtW7bE4sWL8eyzz9a5JhMTE6xZswZLlizBhg0bIAgCfHx8sGDBApSWlmLx4sWIjo5GUFAQPv74Y3h7e2P37t3YvXs3zM3NER4ejjlz5sDBwQEA8Oqrr8LOzg5bt27FZ599BkNDQ7Rr1w6ff/65ahIrEVFjIBKqmxlERERERET1hgP1iIiIiIj0BMM5EREREZGeYDgnIiIiItITDOdERERERHqC4ZyIiIiISE8wnBMRERER6QmGcyIiIiIiPcFwTkRERESkJxjOiYiIiIj0BMM5EREREZGe+H84DonmIlB0zQAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig2, ax = plt.subplots()\n", "\n", "ax.set(title='Bellman Equation Residual', \n", " xlabel='Biomass',\n", " ylabel='Percent Residual')\n", "\n", "ax.plot(ss, 100*resid(cc,ss) / vhat(cc,ss))\n", "ax.hlines(0,0,smax,linestyles='--')\n", "ax.plot(snodes,resid(cc),'ro')" ] }, { "cell_type": "markdown", "metadata": { "Collapsed": "false" }, "source": [ "### Compute ergodic mean annual harvest " ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "Collapsed": "false" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Ergodic Mean Annual Harvest = 0.0311 after 13 iterations\n" ] } ], "source": [ "s = h(0)\n", "for n in range(100):\n", " if s > scrit: break\n", " s = h(s)\n", " \n", "print(f'Ergodic Mean Annual Harvest = {s/n:.4f} after {n+1} iterations') " ] } ], "metadata": { "kernelspec": { "display_name": "Python 3.9.7 ('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.9.7" }, "vscode": { "interpreter": { "hash": "ad2bdc8ecc057115af97d19610ffacc2b4e99fae6737bb82f5d7fb13d2f2c186" } } }, "nbformat": 4, "nbformat_minor": 4 }