"
],
"text/plain": [
"\n",
"Dimensions: (chain: 4, draw: 2000)\n",
"Coordinates:\n",
" * chain (chain) int64 0 1 2 3\n",
" * draw (draw) int64 0 1 2 3 4 5 ... 1995 1996 1997 1998 1999\n",
"Data variables: (12/13)\n",
" lp (chain, draw) float64 -17.41 -11.12 ... -13.76 -12.35\n",
" perf_counter_diff (chain, draw) float64 0.0009173 0.0009097 ... 0.0006041\n",
" acceptance_rate (chain, draw) float64 0.8478 1.0 ... 0.8888 0.8954\n",
" energy_error (chain, draw) float64 0.3484 -1.357 ... -0.2306 -0.2559\n",
" energy (chain, draw) float64 21.75 18.45 16.03 ... 19.25 16.51\n",
" tree_depth (chain, draw) int64 2 2 2 2 2 2 2 2 ... 2 2 3 2 2 2 2 2\n",
" ... ...\n",
" diverging (chain, draw) bool False False False ... False False\n",
" step_size (chain, draw) float64 0.8831 0.8831 ... 0.848 0.848\n",
" n_steps (chain, draw) float64 3.0 3.0 3.0 3.0 ... 3.0 3.0 3.0\n",
" perf_counter_start (chain, draw) float64 2.591e+05 2.591e+05 ... 2.591e+05\n",
" process_time_diff (chain, draw) float64 0.0009183 0.0009112 ... 0.0006032\n",
" max_energy_error (chain, draw) float64 0.3896 -1.357 ... 0.2427 0.303\n",
"Attributes:\n",
" created_at: 2022-05-31T19:50:21.571347\n",
" arviz_version: 0.12.1\n",
" inference_library: pymc\n",
" inference_library_version: 4.0.0b6\n",
" sampling_time: 6.993547439575195\n",
" tuning_steps: 1000"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"idata.sample_stats"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"The sample statistics variables are defined as follows:\n",
"\n",
"- `process_time_diff`: The time it took to draw the sample, as defined by the python standard library time.process_time. This counts all the CPU time, including worker processes in BLAS and OpenMP.\n",
"\n",
"- `step_size`: The current integration step size.\n",
"\n",
"- `diverging`: (boolean) Indicates the presence of leapfrog transitions with large energy deviation from starting and subsequent termination of the trajectory. “large” is defined as `max_energy_error` going over a threshold.\n",
"\n",
"- `lp`: The joint log posterior density for the model (up to an additive constant).\n",
"\n",
"- `energy`: The value of the Hamiltonian energy for the accepted proposal (up to an additive constant).\n",
"\n",
"- `energy_error`: The difference in the Hamiltonian energy between the initial point and the accepted proposal.\n",
"\n",
"- `perf_counter_diff`: The time it took to draw the sample, as defined by the python standard library time.perf_counter (wall time).\n",
"\n",
"- `perf_counter_start`: The value of time.perf_counter at the beginning of the computation of the draw.\n",
"\n",
"- `n_steps`: The number of leapfrog steps computed. It is related to `tree_depth` with `n_steps <= 2^tree_dept`.\n",
"\n",
"- `max_energy_error`: The maximum absolute difference in Hamiltonian energy between the initial point and all possible samples in the proposed tree.\n",
"\n",
"- `acceptance_rate`: The average acceptance probabilities of all possible samples in the proposed tree.\n",
"\n",
"- `step_size_bar`: The current best known step-size. After the tuning samples, the step size is set to this value. This should converge during tuning.\n",
"\n",
"- `tree_depth`: The number of tree doublings in the balanced binary tree."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Some points to `Note`:\n",
"- Some of the sample statistics used by NUTS are renamed when converting to `InferenceData` to follow {ref}`ArviZ's naming convention `, while some are specific to PyMC3 and keep their internal PyMC3 name in the resulting InferenceData object.\n",
"- `InferenceData` also stores additional info like the date, versions used, sampling time and tuning steps as attributes."
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAABQcAAAEiCAYAAABEE2SgAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAAB0b0lEQVR4nO3deZxU1Z3///ettbt6b2hoQBaNgIiIooDrEBWjUeMkRI2Ja4yaxIh+jdHRJBNj9BdMHCdj1NGYSVwZBzsjmmA0bhMRUXYxuOFGK/vSe1d3VXXd+/vj9K2tq5umaeilXs/Hwwd21a17zz3L55x76ta5luM4jgAAAAAAAADkHE9fJwAAAAAAAABA32ByEAAAAAAAAMhRTA4CAAAAAAAAOYrJQQAAAAAAACBHMTkIAAAAAAAA5CgmBwEAAAAAAIAcxeQgAAAAAAAAkKOYHAQAAAAAAAByFJODAAAAAAAAQI5ichADzk033aSJEydq2bJlg+I4+8tTTz2lc845R0ceeaRmzJihK664QqtXr+7rZAE5iTi2Z7Zs2aL58+frpptu0pe//GUdcsghg+K8gIGKGNZ9tm1r5cqV+vWvf605c+boyCOP1GGHHabZs2frZz/7mT7//PO+TiKQk4hje+bll1/Wv/zLv+grX/mKZs6cqcmTJ+vYY4/VFVdcof/7v//r6+ShF/j6OgEA9r3/7//7//Too48qLy9Pxx9/vCKRiJYuXarXX39dv/3tbzV79uy+TiIAdOpvf/ub5s2b19fJAIA99vnnn+uCCy6QJFVUVOiYY46R1+vV22+/rQULFmjRokV68MEHdfTRR/dxSgGgc88884xeeOEFjR8/XlOnTlVBQYE2bdqkxYsXa/Hixfrud7+rH/7wh32dTOwFJgeBTvzwhz/UFVdcoZEjR/Z1UvbK0qVL9eijj6q0tFQLFizQuHHjJElr1qzRRRddpJtvvlkzZsxQcXFx3yYUQK8bLHFs9OjRuuSSSzRlyhRNmTJFt912m5YsWdLXyQKwjw2GGGZZlo4//nhdccUVOuaYY2RZliQpGo3qlltu0VNPPaUbbrhBL7zwgvx+fx+nFkBvGwxxTJK+973v6dZbb1VZWVna62vXrtWll16qBx98UGeeeaYmTpzYRynE3uJnxUAnhg0bpi984QvKz8/v66TslYceekiS9P3vfz8xMShJRx55pM4//3w1NDToT3/6Ux+lDsC+NFji2CmnnKIf//jH+spXvpIWxwAMboMhho0ZM0Z//OMfdeyxxyYmBiUpEAjolltuUVFRkTZv3qw1a9b0YSoB7CuDIY5J0qGHHtphYlCSpk6dqjPOOEOO4wz4n07nOiYH0W9s2bJFt99+u0477TQdfvjhmjFjhubMmaN7771XTU1NWT+zYsUKXXzxxTryyCM1bdo0XXnllfroo486bNfQ0KDHHntM3/nOd3TSSSfpsMMO08yZM/Wd73xHr7/+etZ9d7Y+xMknn5z4RqSqqkpf+cpXdPjhh+v444/Xz372MzU0NOxlTvSe1tZWvfnmm5Kk008/vcP77musEwH0DuIYgIGMGLZ/5eXlJb7w2L59e98mBhgkiGP7n89nfpDK3c8DG5OD6BdWrlyps88+W4899phisZhOOukkTZs2TU1NTbrnnnuyLtb8f//3f7rkkkvU2tqqWbNmqaKiQq+++qouuOAC7dixI23bt956S7fffrs+/fRTHXjggTr11FN14IEH6vXXX9d3vvOdHt059+tf/1q33nqrhg0bpn/6p3+S4zhasGCBvv/978txnB7nRW/69NNPFY1GVV5ersrKyg7vH3rooZKkDz74YH8nDRh0iGMABjJi2P5n27Y2b94sSRo6dGgfpwYY+Ihj+98HH3yg5557Tn6/X8cff3xfJwd7gTUH0efq6uo0d+5cNTQ06MYbb9S3v/1teTzJees1a9Zo2LBhHT73yCOP6J577kk8TCMej+u6667T3/72N/33f/+3rr322sS2Bx54oBYsWKAjjjgibR/vvvuuLrnkEs2bN09f/vKXVVBQ0O10//nPf9af//xnHXTQQZKkmpoanX/++Vq5cqXefPNNHXvssd3az0UXXaTly5d3+7iSNG/ePM2ZM2e327kDzmwTg5IUCoVUXFys+vp6NTU1qbCwcI/SAcAgju27OAZg3yOG9U0MW7RokXbt2qXy8nJNmzZtr/YF5Dri2P6JY6+88opeeOEFxWIxbdmyRWvWrJHP59Ntt92mMWPG7NG+0L8wOYg+V1VVpZqaGp144on6zne+0+H9I488MuvnzjzzzLSn7Hq9Xl155ZX629/+ppUrV6ZtO3r0aI0ePbrDPg499FB961vf0gMPPKBly5bp5JNP7na6r7322kQQl6Ty8nKdf/75+tWvfqWVK1d2O5CfeOKJGjVqVLePK6nbgTccDksyP1vpTH5+vhoaGtTc3MzkINBDxLF9F8cA7HvEsP0fw7Zs2aJf/vKXkqRrrrlGgUBgr/YH5Dri2P6JY++//74WLlyY+DsvL08/+clP9M///M97vC/0L0wOos+98cYbkqTzzz9/jz53wgkndHitq3Vb4vG43njjDa1Zs0Y7duxQNBqVJG3YsEGSVF1dvUfHz3bbtHv8zFvQu3LllVfu0XEB9D/EMeIYMJARw/ZvDAuHw7r66qtVW1ur2bNn65vf/OZ+PT4wGBHH9k8cu+qqq3TVVVcpEomourpaTzzxhP71X/9Vr7zyin7729/yRccAxuQg+tyWLVskKeu3MF0ZPnx4h9fcO99isVja61u3btV3v/tdvf/++53ur7m5eY+On+2nuu4t5G4n0ddCoZAk82CSzrS0tEjSHt3+DiAdcQzAQEYM239isZiuvfZarVu3TkcddZTuuuuuvk4SMCgQx/avYDCoCRMm6JZbbpHX69Vjjz2mxx9/XJdddllfJw09xOQgBqzUNSR25yc/+Ynef/99nXbaabr88st14IEHqqCgQB6PRwsWLNDPfvazPV7wdU+O35UHH3xQn3zyyR595pxzztHRRx+92+1GjhwpyXRk2YTDYTU0NKikpISfFAN9gDi2+zgGoP8ihu1ZDLNtWzfddJMWL16sSZMm6YEHHuhy6RcA+x5xbO/HYu5DYF5++WUmBwcwJgfR50aMGKFPPvlEn3/+eeJx7r0pHA5r6dKlGjp0qH7zm9/I6/WmvZ/tqVX702uvvbbHi8fOmDGjW4H8wAMPVCAQUE1NjbZt29bhm7F3331XkvZJvgO5hDi27+IYgH2PGLZ/Ythtt92mRYsWady4cfrDH/6g4uLiPfo8gM4Rx/puLFZWVibJPEwFAxeTg+hzxx57rF5//XUtWLAgbTHY3tLY2CjbtlVRUdEhiMdiMb344ou9fsw98dhjj+2zfefl5emYY47R4sWL9dxzz+nSSy9Ne//555+XJJ100kn7LA1ALiCO7bs4BmDfI4bt+xj2m9/8Rv/93/+tkSNH6qGHHtKQIUP2+TGBXEIc67ux2IoVKyTxsLmBrnfuYQX2wrnnnquysjItXrxYDz/8cIdbsd966y3t2rWrx/sfMmSIioqK9OGHH2rVqlWJ1+PxuP7t3/4tsXjsYPXtb39bknT//fenneuaNWu0YMECFRcX65xzzumj1AGDA3EMwEBGDNu3Hn74YT3wwAOqqKjQQw89lFj2BUDvIY7tOzU1NXryyScTa9Wnev3113XnnXdKkubMmbO/k4ZexJ2D6HOlpaW6++679f3vf1/z5s3TY489pilTpqi1tVWffPKJqqur9fTTT/f4G1afz6fLL79cv/nNb3TRRRfpmGOOUUlJidauXatdu3bpggsu0Pz583v5rPqP4447ThdffLEeffRRffWrX9Vxxx2nWCympUuXynEczZs3j5+1AHuJOLZvbd++XVdffXXib3dNnVtvvTWxXuqsWbP0gx/8oE/SBwx0xLB957333tMdd9whSTrggAP0wAMPZN2OdViBvUMc23fC4bD+9V//Vb/85S81efJkVVZWKhwOa8OGDYkx2aWXXqrTTjutj1OKvcHkIPqFmTNn6plnntF//dd/6bXXXtNLL72kgoICjRo1Stdcc80eP3Uq0/e+9z1VVlbqkUce0erVqxUMBnXUUUfpmmuuSay7N5j95Cc/0aRJk/T4449r6dKl8vv9OvbYY3XVVVdp2rRpfZ08YFAgju070WhUa9eu7fD6xx9/nPj/gw46aH8mCRh0iGH7RkNDQ+IOpjVr1mjNmjVZt2MdVmDvEcf2jSFDhuiGG27Q8uXL9dFHH2ndunWJn1ifeeaZ+sY3vqGZM2f2dTKxlyxnTx+nAwAAAAAAAGBQYM1BAAAAAAAAIEcxOQgAAAAAAADkKCYHAQAAAAAAgBzF5CAAAAAAAACQo5gcBAAAAAAAAHIUk4MAAAAAAABAjmJyEAAAAAAAAMhRvr5OwEBQW1u7R9uXlJSovr5+H6UG/RXlnnv2tMzLysr2YWo6t6cxTKI+5yLKPDcN1jhGfc5NlHvu6UmZ90UcYyyG7qDMc1N/GYtx5+A+4PGQrbmIcs89g7nMB/O5ITvKPDcN1nIfrOeFrlHuuWcwl/lgPjdkR5nnpv5S7v0jFQAAAAAAAAD2u377s+JIJKJ///d/17p161RdXa36+noVFxdr9OjROvfcc3X22WfL7/d3a1+2bWv+/Pl68sknVV1drVAopOOOO07XXXedRo8evY/PBAAAAAAAAOif+u2dg83NzXriiSdkWZa++MUv6tvf/rZmz56t7du368c//rG+973vybbtbu3rZz/7mW6//XY5jqOLLrpIJ554ol544QWdc8452rBhw749EQAAAAAAAKCf6rd3DpaWlmrlypUKBAJpr7e1tenb3/62lixZosWLF+uLX/xil/t58803VVVVpenTp+uPf/xjYn9nnXWWrrzySt122236wx/+sK9OAwAAAAAAAOi3+u3koMfj6TAxKEk+n0+nnnqqli9frurq6t3up6qqSpJ07bXXpu1v1qxZmjFjhpYsWaLNmzdr5MiRvZd4SZGIo083ONqxXTr4YGnECI8iEUctLZLXa/6NRpPbBwKS1ys1NUvlZVIgYKm+3lFTk9TYJPn9UlmpFI1JhQWS12sl9iNJJSWWgkEr7fgtLVI87qimVorHJZ9XchxpV400doxUXm6OIUn5+VI8bnVIWyBg9i1JO3Y4am2Vhg+Xioo82rLF1kcfm/QUFUllZSZdmcd0t3fTlXrMlhYlzlGShpRLFRXmeJl5FQgk05mfb7Z391VSYikaNcd18yc/XwoGrbRjxuPS1m1SKGTyM/W843FLLS2Otm6TSkukkSOT+wz4Tdnk55nP1NRKXo9UWGjKraZWKi+LKR63O5RFY6OtrdvM9m4V3N25uOluaVHifTf/y8rMeUSj5r/GJikWM9v4/VJRYTJdTc0d8yMzXfl5yTx3j1FQYD4bj0u2LRUXSSUlHdPb1CS1Rsz7qeWWeqxIxNGmTY527pQqKky+pn7etpPn1Fld3rHD0a4ac36jD0jWp9S6np+vDnXALddsedhVG3KPmZnfLre91tW1yeu109IzGDQ22olzCwSsRFtMbStunqbW0cw6nGnXLltbtkojKqUhQ9LzLLWddvb5xka7va0lY2Rn22fuz01ravqzHaOrY2bWu8z6lG2fnaUjHncS8T41PnY3jan1PrXeph6rqcnpkN+Zn3M/E41KkWhccpwO+6uplRrqpVCBlBdMxvRs+0rtJ9y6476/a5etzz438XXoUKW1vdR8SY2vqXEyMy6mpjO1v+isTuyJ1Hxx+8HOjt/ZZyUTOz7faPJt1KiOfVNX7Sdb/nb2txv7stWnrurI3uRRfxeJOKqrcxLtyW1znfWJmZ/dXZ+8u3rQ3TRm20+2/jmzjmeOIbpTL7O1fbcfzwsmxyjRaHrfvru631l+dVbXMmNFZh/aWSzMlobUMVY0Jvl8EZWV2mkxL7Ut5+dL9fVSQ2PHsU1X55U5Ts48r87qh9R539hV/5Kt38gWt7PFp+6kMRi0Esfvqk10Vm/cPMk8zu76sWxp+HyjeW/0AenXP7ad3pdkG394vY4kW5GIMyjjmRvHWlocffKp1BaTxoyRPB7ps8/Tr8HicSX664oK089mu7aTkm2hqz6zszFXZ21BMuPshkYTT9w+3L3WS40p2fqs1Lac+tnUei6lx8KaWqklbPIkc1yZOuYsLLQ6tI3U66T6emnLFvOaz2+uE93rjWxxILX+urHajSVer6P6eqVdQ6TmcWefcf+VbDU22h32ka383G0KCjrOEbhxLjW2u5918zLzWqylxcwTSOZ6UjLlFo+belVUlH0ss6dtNrPeuWU1pFzKz+94Pd7VGMa9zvzsM6mtzZzH8OGmLHfV7H4M3FXf2tTk6LPPk/MGbr59usHR5k3SsGFSaWl6HHbPpahQkmXmEDLLJtuYJBRy1B/028nBzti2rddee02SNGHChN1uv2zZMoVCIU2bNq3DeyeeeKKWL1+u5cuX66tf/WqvpXHzlrj+4z8cLV1mKmlxkXTWmbaGVUibt0iffWYaa12ded/vl/LyzARPcbGpJCUlpjJ+ukGKtEoerxQMShVDJX9AGjXSUaTVVLb8PGn8eEez/kmqHG5p6zZHa9c6Wv+htO5dqbY2GfxbW6VgwFzkTZzgKBCQWlukQFAqL3dUsyuZNssyjevAcabBv79eirdJBxwgVVTY+vurUk2NCQQlxaaBVFaaxuwes7hYmjBemvM1WxVDLS1e7OiDD80x447U1Cht2mwCu8drOrJpRzgaMUIKh5N5FQ6bhllcLI0Z46gg37y+ZatkSSosdFTfkOx8Jox3NOFgaeQoR++9K33wobRtq9k+GjOTpCNHSqNGmvNubZGamh29974JvsGg2YfPJ+3aZSYU/T4TaBxHitvJSbNYmxSJSFKjRlRKM6Yny+If79h66ilTjk3NJniE8k1n0Nm5jB/v6NBDTdBpajZ5vnmLyf/WFsnySAUhqaHBdAotLVJbe8fj9ZoOtWKoKb9gUHKUzI+pU9PTtXGT5POZuiCZYzQ2mHoimbz0+Ux6R42QJkxMpvfd98zAxI5LZeXSlMmm3Hw+U4enTjX7WLDA0WtLk2U47UhHwyrM5zdvaZ80DZoBQllpx7r87LOO/v6aqWvBoDRxvPStb9maMtmTqOtNzVJdrbRxc7IOjBrpyOuRhlaYNKfmYTAoHXhg9jYkSc8+62jF6vT8bmnp2F6HDGlWWal0+ukmPYPBP96x9fzzUmskLDsujRrlyOuVdu6U8vJNW3Hz1JEp73hbxzrslqHr1dds/ekpUz6FhdI5c2zNOtHk2dZtTiI2dPZ5N10NjeZLh5ISR41N2bfP3F9lpaOiIlMHd+406R9ZmWwPu8sL9yLSLWe3Xrp1JBiUDjssvY11lY7t26X1H5kB27AKs9+KoZbWrnW0ecvu05ha71PbWuqxPB5H1Z+1D0ja83viBCvtc258fOttE7MrhzfpkEOS7WDxYkevviZtqHZjnDRkiHToJOmLs+z2uJnc1/Jl0orVpj4MGSodfJCjsjLzfmvE0V+fN+cmScOHmQFcfl56vqx7V2psNPF16BBzgTNmjMmLkaOScTH1vN18++yz7H3innLLzc2X8nLp8Cnpcdk9frayccth21YTY1sjJh4PHWr6B7dvSo3xXfUBbv529rcb+9zxRmp96qqOuHlUVrbHWdTvbd3m6MOPIvroI/PFVNx2tGmzOe9sfWLmZzPzKbNPdmNfZ/Wgu2nMLKPO+ufUMVE0Zk4iEEyOIXaXnsxjuW3/zeWmH49FzRjMvYBrbR8zlpVLx850NHNm53W/s/xyx16ZMTo1frpjyjlfS/ahblozY2G29r9jp5MYY9XUmn45P69JFRXSt75pYl5qWw6FzPittkYKt18gumObzHibel6Z4+Rs4wY3f1PrR1d9Y1f9S7Z+wy3fzLidGZ86i4WpaSwsMP358uXS9h3p4+bMPMhWb1LH8V5P8jipacjWj2XuLy/f0XPPSx9/YtJ3wChp7Bhz/bN5i+kjK0eYvsQdh6SOP9w+ffToVpWVOj1ui/2VG8deX+rojTeTfaPfbyYH29rMNUlJsVRaJskx5RmPm7o9bqx00EHp13buOLakxFxPppZfap/Z2Zirs7ZQWWmuFd9eZ47j8UiVw6XCInOtt3Vb8nrhCwc6irWl91mSEm25NWL6/8IiE5vCYVPPDxyXHMt99pm5ftyy1Vz7DB9m2rw7rkwdcwYC0qRDHJUUJ9uGO34JBk3scK9b29pMWysqkg49xKRt21alxQF3TJKX7+jvfzex2rbNeG7kKEebNpl62dZmriFSr28qK7N/prUlOcYuLQ3r40/S95F5HVJSYmJZU2Py9bKy5BxBXZ2Z7GxsMnUkP99cx/r9yWt7Kf1azPK033TSaD7j9Zr/2mKmTCTz+WHD0scybl/S3TabGZ/y8hy99IopA9uWhg1ztHlz8nr82GMcnX+eOX62vuaVl6WX/27qnW2b8nZjZn6+yZfOxsCZMTO1b22ok976h0lXW5s0bpx0xOGOqqul1W+ZfLcsM6E5fryJw+657NhhzrWk2JRTcVGybLye7GOSI6ZGNP5gp89jWL+fHIxGo/rd734nx3FUV1enN954Q5988onmzJmjY489tsvPhsNh7dixQxMmTJDX6+3w/tixYyWpW3cgdlck4uiFF1u09E0zaRQMSnX1UtX/Smd+2UzobPjMVJi2tvbJuqBpAD6fmbjbvsNMltlOsiNwHCncbAY+hYVSzS7zzUZBSApWSh+slwpCjmb9k6O1a6VtO6T1H0pbt5pjuN84OI4J2LW10pvLpKOPMulrbTXHbWpKbm9ZJlBs2WpeLyoykyLrP5ReX5pMV1ubmZmPtZlA5Mgcy7JM+jd8Kv35z2aA+/56SY60s8Y0nEhr8hsJn0/aucPse8J40yls+Cz5bWltnQngwYAJ6Dt3mgm+WNycSzDPzN43NEgff2Qa/pq3pJZWM0Bz893vM2ncsMHsu3SXee3jT8y55+WZf99cZgKg45hJOL/PNPBIxAww3YlDj2XSFAiahv7OOlMWxx7j6NlnzXEjrVJzs9lPcaEJ1qEs52J5TaD+9FNpzFjTGS5ZIn30kZQfMnlVX2/y1rbN+aQuvekOFpqazCRkcZEkj8mP0hIzMPQfZdL12UYpPyiFW6Wlb5jP54fMQKCp/U5OyzL7jMWkcJP5FiTcaupVY/ukrVuf3nhDmjhBOu54qb5RWrXKUWOjtPTN5HbNzaZ8K4aaPI5ETXpjNaZ+Fxam1+Vly6Q3lpn67tarDz809Wlkpa21a82xQiFpxSqzv+Ji00HU1UiVIyVZ0vpa6dNPzPlFolJ9g0m/P5DehoIBR21tZl/unRP19SbtwaAp+9T2Wlxk7mh8/nlp3JiBfwdhY6MZDNbWSweMtLRqjRl8feEg81rLZtOZy5I2NJlBw7QjpRX/kLZtS9ZhtwxPnZ28Y+xPT5m2WTFE2rFL+tNT0mGH2iostLRqlYkNXo/kWB0/n5quoeXSP94xMfKQQ0wcTN1eUtr+2mxp5SozMVNcbOKd12fKce1aR2Wzst9pk3rM4RUmpj7/vKl3y5Yl60gkatrMh+uTbaxsVufpKCs37auhQfKXm9j57LPSIRMdhVtNW6rtIo2RiLm4qm80A5JdNeY4sZgSxwq3Sm+vNfs48EAzsf7kn6STZzmKO+Zz27ab+OhOLsXapF01tt59z7QDSVrztlT9WXqcqa2VPvrYxO+ZM0yZb9tuzq26/csCf8DErNpd0tlnS5u2SC++ZNp8fp6J5es/lA4cZz7v5ktDg6lHtm3qVvVnUnN7LPN6pbf/IY06wFwEuOftOFJDk4mhn200MTW1T3TrUHdFIo5WrXL0znvJfKmplf6xLhmX3eNnKxu3zKMR6dNqE/N8XhO/Nm+W3nvP9E3uuZSXm9iZ2n5S+4Dhw0z+pp576t/FxaYuxtqkQw81fZJbn6JtndeR1HZ20IH949vq3uK2kZZWW7W1Jt5s3tx+seMxZZHaJ3ZWhm4+ZfbJy5eb48yYbvqfruLI7tKY2o4765/fXGZicGGRed3fPnKPtZn63tJq2kdn6ck8ltv269u/nG1uH1vatrkzRO2nkZdn+tJly0z9PPCgjnVfUtb8+ugjc6EV8KfXNbdfX7HK9P3BPDM5+eyzpg8NBMzF267a9FiYrf0vW+Zo3Tvt48SwiR+SuTDctl1a8KQ060RH775v2nIkaj7X2GjGNqGQaS/hJlMXUuNt6nnJkWrqkuPkcNjEsdRxg89nyqq4OFk/uuobjz3G6bR/WbVKHfqN4cOlY4+RdtV2jNup8WnkqOyxMDWNQ8pNua9Yac4j3pYcN2f2YdnqTUurSdvOGpN/Q4eY83PbychR2fuxgmOdxJjNTcMbb5ovhbxeU//e/8BcpHt9pn5IZsy3a4f5FVZ9oxLjD8ncLVVXL4VCcdlOz9pif+W229q6Nr211rRFxzH55OaNm287d5nrirY2M171ek2Zvv+BmdAKt5r661535uWZcUs8bsbjI0em95nRqJN1zDV+fPa24Kswcaq+waQn2j6239Bi4lWszVzTejxmPLJtqznuYYebOvTnP5u0bPjM9J2RiPlSMuAzEzUFheYcd+w0k2kFBdJHn5i25fGY/7Ztl6r+ZMaVkhJjziFlpi+urZWmTDZtw+sx70WiZjJr67bkL7Akk/7mZtMOa2vNl5R+vzn3aNSMSRxHWrW6PQ+Cpp+pqzfltH2nuTb0eNuvldqvb/6xTnp1sWlDqZ9pakpenwQD0qbNUX2+MbmPhob065C8PJPXzc3JmxWi0fbr0wJp107z5Ucsavq71lZz40cgaOKe1Z5ftm3SFomafTQ3J8df7sSzlPzXknk/dSzj9iWbtyTvbOyqzeblpcent96WqjeYNPm8Jv+2bjX/nx8y+3zjTWnEcEfBoDr0Ne++a75UdicG3fLbtav97sr2fizbGPjtf0jDK5MxM7VvtWTqfVP7rzfb2sy8QV1t+0R3+xxILGbaQnGxtDouff6ZOd9IxNSN5vZfHbaEzfnU1UhDh5kxdOaYZERlXOHmvo9h/X5yMBaL6d577038bVmWLrvsMl1//fW7/WxjY6MkqbCwMOv77uvudp0pKSmRx9O9i/7aWluffd6cmCUOBD2yZKulVbKsgBw5KiiIq7XFVl7QBMOA3zTgwgJJllehfFsNDY4smQrp2O0/lfC6t/1a2rbNUUmJR6GQpaFDvGppdRSN+WTbAdlOTKWlUltbq/LzbUXag1AkYo7lkQkQLa2Sz+tVQciS1xNXJGoplO8oGnEUzLfk9Zrbgpua7Pa7OLzKz/eouSmmWMwEx1hM8rRPlvi9Zp+FheZ8igrMZFpenqXGRo/q630qKHCUn2dpV02bPJ64bNucl8cyafK0T3pZ8ire5lFBQVyxqKOyUkuxqKPCAkttca/8fkcer62hQ3xqjUi2E1Uo35LP49GIEeaCsrDAr0/CbQqF2huYE5PfZwJtKM90pCYAelVa4igatZWfb/IoEDDfxng85lwKCtrPVeYzHo8pr+YmyeuXPD6ppNhSJCr5fF5FY341h31qDreooMCWbTvKizpmcjIgeTxW1nPJy/No67Y2hcNS5bCAAgFLPl9EttOmslKPbNuR32cu8tyJOzkmsMgxr7l3NwaDki2PRlRIkail8rKAbMdSc9ir5nCLiosdFRd5FGq1VVsbl0dScZGl+jrHdOax9jsP27898vmkaNQnv1+yHfPVeH5+e9kFTP3y+XwqKgyqrNTSxk1x7dwVV1s8mrh9u7XFUXPY1BePZam0xFFtrWT5TYAsKw0obifrcmNTRLYTkT/gKC/PtBevJTU1+xRuyZftxDV2tEd1dbY83ogKQrZ8Ho+GDXO0caOj4RU+MyCxbNlOXGWlHtXXS16PaZOlpeltqDnsVUuLI4+3TUMrpG3bbPl9jmmngfZvpFLaq+1YOqDSr4ZGyXFCKivrnyG1uzGsrq5NrZGwxo72yPJYKi72qrnJVrzNpxHDLX30cVwjKj1mQJVn8iTg9ysYjMnjjSfqcG1dXNGYT3l5eSot9aj6s6haW5o1YoSlgpBXvkBcu3Y6amkp0JAhPrVGWlVQYMpHUofPZ6aroCCq5iZbJcV+lZR407Z3HKXtr7XV0bZtMfn9HlnyaOxoS60RR5XDfO3nEVRpace8ST1mKORRXp6t7TtshVvy1NgUk8cb09AKaddORwG/rbY2j8rLgrIdS3l5wU7T4bEstcWkESMkx7Y0otLSzl2OwmGfhld6tWNHrMs01tbasp2Ixo72KBi0VFjgaOOmuFpbHRUU2Cor9Wjb9rjidkwlBVJRgU/BoKNtW23V1ft12ORA+6DD1nvvReT3m8FkeZ7U2uooL8+v5rD5Ms1STHLiCvjVHtuSMSYcNuc7pNwnydbata2y5Kiy0iPblurr2yRZCvjzFQo5am1tUUFICoU8amm1zUWM11JJiU87dsblsSxFI+YcWsLmy6jaWtN/yPKrrNSr6uo2VQ7zq7zcmzhvy5KGD/No48aoiottc7dwSp/o1qHuqq211RppVTDYplCBk8gXv99KxGX3+DW1doeyccvc73fk2DF5PO0xX+aLrEjEo8KCgDZuiqtymF+BgKVgMJrWflL7gPJyryRbH38cTZx76t+xmOTxRlSc76i40K/CkLRxU1zhsE9f+IKv0zqS2s5aWhyVDYDbB7sbx9w2UlRkyesLaGi50153pNo6S8Mz+sTOytDNp8w+uaDA3F5QVBhQWanVoR50R7Z2XFNrZ+2fNzTGE2Oi4mJbdtyRJUv5+Y7a4t5ELO4sPZnHctu+x3LksWzl5TmybTOecC+uvF5z0VwQsmTbjpqafVnrfmacc/Orod5RSalHw4eZWOLWNbdf93hjqhxqyeezVF8fVzjskeMUKBj0yHYiGj7MSouF2dr/hx+1qbExpoICR62ttjweE5987XdANjdb2rnLp2BQChU4Kix21NRofmbh95mxqnuXUDzuV+UwfyLepp5Xfp6lhsZ4YpxcWGDGM6njhvx8S2NHexWLOYn60VXf2Bw2/V62/qU10tah38jL86qoMCC/X1njthufykr92rq1rUMsTE1jMGipORxXa0urRo/yqKZGGjHCUXOT0tqE4yhrvQmFLJUUe7Srpk2BQFx5eV4NHeJVXb2dSENtbbxDP2bbftlOLLG/5nBczU1h+XxmLBZvc9QSdhJlMqTcav+lj6mf8Ta/xo72qL7BVl6eeT29Tw912af3B3t6PWk7EckyP793P+b3mf7YlRc0f1tW+9g/YK5z4u0TiOFWS6F8KRpxEtedecFk3XfLL7XPbGmxs465CkI+8xPMjLZQUuLR5xtj5o7FPHNd0BY3E31qT1deMHknmpmgtFRcGFBhSPrk0zbF41JBgWN+9ZRna8cOM2Hk8VgaOsS0hdZWye/3KN5mKRiMSY65acHjseT3OwqHPWppKVDcVmLM6fd5VFTUpuZmR+EWr4qLpWjElscyk6HuZHeblbzxwh3neDxSY7OlkSO9GjrUq4bGuPKCtmT5VBDyKNwcUWGhpcJCM64Jh22FW5I3+9i2pYDfSVzf1NXH1dxsq7g4/TORqEejR3v0+UZHQ4dY+vCjeNo+Mq9D8tqv6S3LXJPY7fnuzhFs2WLu2rPbv1hxl4jyeqVQvqVAwFyPxeNO4lps1672iSrHzD1YVvqEqWUlrzs9nuRY5pNwmzyWI5/XUVG5s9s2G2+z0uLTxk1tammxNWa0pXCLR0VFcTU2muvuYEAqKrTU2ORo5y6/Sku9aWOYrdvatHVLXLGYk2gfjpOczDT1zaOKodKuGqfDGPjjj6MqK/UlYmZq3xpusRSPx+X3mTzOyzMTqe4NVIFA8liWJEseObbU0mJr5AgpErFUVGTu4CwoNJP0ZWWWdu6Uigs92rzZ7jAmyc+zZDuhPo9h/fNKNkVBQYE++OAD2bat7du365VXXtFvfvMbvfXWW/r973/f6cRfb6qvr+/2tpGIozGjA+ZbiUYpGLUTdwc6TlSW2mfmHRM0YzHTyBy56w3GFW7/psV2kt/munfo5eVLtbWOggGpOWxGcTt3mcFOwB+Tx9Mij2VmyX0+c2ecI6mlPYhEY+ZOG/dnTm3xuJrDprKHQu3rBThSuMUEaJ/fSQTInbviysuLK9Z+S3Mkmpw1l8yFTyDQ/rqkxvaf0ba2OjrggLhKSuLatMVMqLVGlBiMRmOmoTlqn7z0SI7i8vriam42HUxtnaO2uPnp7yivyVM7Lu3cFVVb3HRE4RZHgWBcO7dIZSVSU3Pc3IbdfuegLPNtgL89bbLcOxbNcQKB9vxyzDl42idjPV6pua69TOxkkG1y99Em2R6pvsGRz28emhPwt6kgZO5K27hJirZ/a2fbZmLJF3CynovHa94PhaSt29tUXt7+TaAl1dbZikTav9FoH8Tbtikvl+Mk78KMRKRgwNa2Heabs5raFvONW3u6Nm+RYpG4wq3Jb0saGh3FbXMMWe354JE8cZOOQKBN4VaTHlkmvzweMwns95pzb2xqU0OjVJBvvlH2ec1EbCTiJNLm1ou6evPZWKz9pz51UTMp3V6XiwrNsWLR9glFj5QXkAoL2hTKb5THMncrhUImD1uiktcXV81O843jth0xDato/8ayPQ9bW5N3MGW2oYJQzHTGcXMXq5vfHo85fmZ7HWE52rK1TUPKJcuqV21t18G8ry7AuxvDvF7zpUX159IBIwNqaIjL65G8vpi2bDNtYcvWuIYMMXkYi0nRWEyRSHodtuOmDFtbW1Vba6mwwFZevrmzuGJIXDt2mW/Z8vMbFYlYygs6am42nadjdfx8arqGlpsY6vFI9Q0xNbfE0raXlLa/NtvsLxaz5eS372OItHV7VEPKpNbWFtXWdvyGLvWY7p0dZSVSKD+mosJkHWlpNfWkpNhOtLHW9t/kZUuH7ZgYvGWL+fZyy1Zzt0Yo1KZt2008qd7eeRojEUcey1H158k7jgryzaDOPVa4fbDb1Cw1NreppsZcCJeWRFX9eTTxrWlByGwTbm6/+yLfo9bWmApCJqg7kmQll2Jw6388LhW2n29evtlXcZG5I3DrVlv+QPs3yAWOorGwwmEzsGoOm6VBWlvb74CIO6qvjyXyJRA0dwbYtuk/zTfejuREVVvX/g3/9pgcJc/bcczxLZlv14OB9D7RrUPdFYk4ygs65lvf9nxxlxIoLDBx2T1+SVHHsnHLPBpJxmn323tJCgZtNTW3Js6lvFwd2k9qH+DInF/quaf+XVxsPtvQKDU0RbRrl7k7PRRq67KOpLaz/Px81dbWdjuP+nscc9tIY2O+4m1R7awxZbizRvJ6nA59Ymdl6OZTZp/c3GyO09gUU0Njx3qwJ2lMLaOSouz9c+qYqKEh/e6G4cPMF12xWOfpyTyW2/brG027cy92zBezkqz2O4wikixH5aWmz81W9yVlza+iIqm1Na4tW2Jpdc3t1+24uTsk2P5LjeHDbFlWvSIRSx7L6RALs7X/okJznG07JLstObFpxorSyBGOhg6JmTucmqXWqBJ3Rcbaf2ERi0k+j+T1xrR1eywRb1PPq7lJ6ePkZpNfqeMGn8/0FcXFyfrRVd9YEFKn/UteUB36jdZWW41NMdXUdozbqfGptq4tayxMTeOQcnONkJcv7dxlqy1m7kgsK1GHPixbvWlpNWXcGmn/hUWr+SI4HjftpLauLWs/5vEoMWZz01BQKDVvb58Iax/P+tp/3bOrxix9lLgbzhdT9efm55nu0jfuOY0aEdDW7eEu+/RMfRHH9vR60mM5khNUXp75ssxxpJidvp37c0+n/QaBSDR5R2EgIIXyHIVb0q87WyPtk0ZtyfKTkn2m1+tkHXM1h9uytoX6+rg87RNqre13/rW1mfG/e+NCa8Tsw99+M0Brq5Pss4rMfnbVSJEWc80qy9wcY9uOdu6KJ+pFLGYrUNC+1Illri+CQUeRqHRAqa38fHPDjzvmHFJmJpv8fimUH9f2Heb4tmPu5soLmuuP1Ekl91rKtqWiAkexWJt27mxL3Nww0omqOWwm6JuaHcXbbDW133UXyjft045LHq+ZLAsGzfWNmQDt+JlgwNbOXbY8XnOHZEGBpZra5D4yr0PcvHSc9mvo9rs18/PMHEFe0Nw5GG//6XYsZtpePG6ulYN5jiyP5JUUccy1mO0k8yEW63jnYGISzHKXrTJjGbcvaWsvv9212by89PjkyHyZs7PGkc9rysrrbV9SzSPF6kwcGDokpmAw1qGvKSszv05xr4lTxeNSPG6ug4uLO46B84JSbV0sETNT+9axpSYdLS2mnkWb3cl083oslqwnXq/kyJblaf+lY72pt43tNxk1t999WFvrKOCTGpriWcckLa2Oigpb+jyG9fvJQZfH41FlZaW+9a1vqaysTP/v//0/3X///brhhhs6/UxRUZEkqcn9jWQG93V3u94QDFr60qn5evediJYuM421tEQ660wl1hwcNyZl7YeC9DXMgnnSQUPMWgKffW4qXfY1B5VYL81jmVu9jz7aUlGRpalTza3oE8ab9911FPLzTeP0+9w1B00jdu/aKy9XhzUHh1WYn365aw5GI+Ynv8cfp8Sagz5f6pqDSltzsKhIGneg+WlZxVDzDfQHH5rOpqysszUHlVhfJxpJrjlYXmYa9/BK6aCUdfoCPumYmUqsOVhcLH3hYGnkCOnoo5VYH2XcmM7WHDSBpqhIaWsOHj4l+dPhtrbsaw4Oq0hdc9DSiEpHkw8zZTFkiKUzz7QVecqsoVEQ77jmYOa5yDbreR16qMzaOk3S4YeZAdb7601HmV+552sOfuFgk9apU5Ppctc0Cgakw9t/of/++vY7Xv3m78w1B0eMMheZ48ZmrDlYZm7XHzHCDCJKisyxJGnbtuSagwUF5qc2qWsOhvKTaw76fel1eeZMRzt3JtccDATM+2efbRaXdet6U7M0/ajkmoOlpaZsvR5Tp8aONgNwNw+Li8zPLTPb0DHHmDTX15v1kOyU/M625qA/YLWvOdhxQfWBqKjIo9NPNz+n3VXjaMwB7sLD7ROzI01bcfPUkan7kyeZhZvdOjxxgilD97b4IUM8OmeO+WnxrlpT98+Zk1wg+KijpObm9nVEnI6fT01XbZ1pyyUl7Q8yynK81P15LbN8grtOjSVTX4eUmTra2a37qcfcuctcOJ1+ukmzWy/dOlJcJI2fkGxju0uHu+ag12sukk4/XYk14iLtF+qdpTEYtDR1qvn5VF1deluLRs2x8oPmJ2jVn5lY7ua3u95KXZ3Z95TDlLZ21ZByjw45xE60g2jUUUN9+pqDZWXSwV+QvjjL1IXEvs5WYs3Btpj5OcnBB5k+Z9QI6fzzlFhzMBRKrjno96XnS1tbcjmNUSPbF1WvMO33yCNMXMw8bzffxoQ79ol7+tOMYNBKlNtb7f1neXteuXE59fiZZZNa5geO7bjm4KRJyb7JjfGZ7Se1D0gtq2x/p8a+utpkPXXrU2d1JLOdhcN7lE39mttGPvzIo9pa057Ky5S2vk9qn9hVGVpOxz75kAlmH6l9XU/qWbZ2nLV/Pix9TJS6LtLwSiXWcuosPZnHSm37brzpbM3B0lJp5kzzX2d1v7P8Sqw5mFLXUvv1FavN+YwZLZ15ZrIPddOaGgs7a/8HfSG55mBxsemXvR4TN75xnol5lpVsy2UHpK85OGSIiU/DKjvG29TzKi9NHydnGzesXeuk1Y+u+sYhQ6xO+5ejjnKy9hvNzdnjdmZ8ird1jIWpaayrM/l5zhwl1hx0x82ZbaKzepM6jve2T/pNnJSehsx+LPX6xE3D+ecpseagZUmHTDQPTExdc3DcONOXuOOQ1PGH26cPGeJVWWnP2mJ/lYxjPh0xNaKWlj1bc7Cw0F1zMP3azr3uTF1zMLPPDAatTsdc2dqCv/06LHXNwby8rtYcNNdOqX2WlFxz0Jbp/zuuOZgcy0UjZjJsy9b2NfKHSeeekxxXumPO2npzfTrpEJNPdvuv8YoKk2sOen0d1xwsKDBt+fTTlVhz0D33oRXSyErpC+cosX5gcZG7fqBUsim5XmAolLy+OXSSuU7O9hmz5qDa1xwMKC8vkraPPVtzUGlrDgaD6WsOutf2Uvq1WGlpcs1B98GLna05mDqWcfuS7rZZtx/zWNIRh0vHzFBizcHiYjN+dNccDATMeHbWrI5jGLevKX+58zUHZZnz6mwMnBkz3b519BhTX9w1B30+d81Bpa05GAiYcfzw4Wbe4LhjlFhzMB7PtuagaXPDh3Uckwwf7tX4g/s+hlmO484LDxyNjY06+uijdfjhhyeeRtyZE044Qc3NzVq5cmWHdQcffPBB3XXXXfrVr37V5QNJ9uQbdcnM5G7dWqMN1Y527pIOGpd8WrG5NdxJrLfmCrQH6HDYfeqNpYYGs/hrU7OZACspSQYJn89K7MeR+Ulr5tPAzPoSphG0td+xZcdNJR892jytuKHBMbP2eebW5cy0BQJm35K0c6ej1ogJKu7Tij/5xATQomIzCerzWR2O6W7vpiv1mK2tSpyj5ZjgP3SoOV5mXrlrB9i2pbw8kz53XyXF5slXdXXJ/MnLSz6JyN0u3mYuQPPyk4vhu+dt25bCYUfbt0vFJdLIEcl9+nxK3P3irsPj8ZhvPb0e83dZWbHa4g0dyqKx0db27Ur8vMzN167OxU23uzZIav6XlpjgF2m/m66pWWprLy+/3wwI3HS5DwFJzY/MdOXlJfPcPUYov31dzPY7eooLTcDOTG9jU/s6f4Xp5ZZ6rEjE0ebNZm2+IUNMvqZ+3nGS59RZXd6506xBFPCZyarMp/q5x8ysA265ZsvDrtqQe8zM/M5sr35/iSyrvtsTg311x82exrDGRluOY84tELASbTG1rbh5mlpHM+twpl27bG3bZjrRbE8r3t3nzVOU02NkZ9tn7s9Na2r6u/u0YveYmfUusz5l22dn6WhrcxLxPjU+djeNqfU+td6mHqupyemQ35mfcz8TiUrDKkrlOPUd9rer1gxA8/LMANON6dn2ldpPuHXHfX/XLlsbN5r4OrT97tNs+ZIaX1PjZGZcTE1nan/RVR3qrtR8cfvBzo7f2WfdPmfTJimQZyYhMvumrtpPtvzt7G839mWrT13VkWDQUllZ2YC4c3BP41goVKotW2oT7cltc531iamy5VNmee2uHnRHZ/UpW/+cWcczxxDdqZfZ2r7bj+cFkmMUdz1Vt2/fXd3vLL86i9GZsSLb04qzxcJsaUgdY0Vikt9fqNLiprSYl9qW8/PMl6sNTR3HNl2dV+Y4Odu4IVv9kDrvG7vqX7L1G9nidrb41J00uuv5Zhs3d6fepI7js7WTzvqxbGnYuMlcJI8alX794zjpfUm28YfH46i4uEytrXV71Bb7Io7taQyTknEsHHa0YYMUbTMT6h5L2rgx/RrMtpXor4cMMf1stms7KdkWuuozOxtzddYWJDPObmgy8cTtw91rvdSYkq3PSm3LqZ9NredSeizcVWsmpg4Ylf1pxe4YqLDQ6tA2Uq+TGhrMBKblMdcaoVDyeiNbHEitv26sdmOJx+OooUFp1xCpedzZZ9x/i4vL1NBQ22Ef2crP3SaU33GOwI1zqbHd/aybl5nXYm6eWo65nnRkjhtv/+K2oDD7WGZP22xmvXPLqqxMCoU6Xo93NYZxrzM/32gm+MrKpMphyTXxdzcG7qpvbWpytHFjct7AzbcN1Y62bDXrqJeUpsdh91wKCkwafL6OZZNtTDJiRJnC4bpux4Z9FcMG5OTgxx9/rDPOOEPTpk3TE0880eW2P/zhD/Xss8/q8ccf1/Tp09Peu+iii7R8+XK98sorGjVqVKf76MnkYE86AAxslHvuGawX1RL1ORdR5rlpsMYx6nNuotxzT0/KfKBMDlKfcw9lnpv6y1is3/4O7qOPPlJLS0uH11taWjRv3jxJ0qxZsxKv19TU6OOPP1ZNTU3a9uedd54k6e6771Y0mpxyf/XVV7V8+XKdcMIJXU4MAgAAAAAAAINVv11z8LnnntNDDz2ko446SqNGjVJhYaG2bdumxYsXq66uTkcffbQuvfTSxPbz58/Xvffeq6uvvlpz585NvH7MMcfo3HPPVVVVlebMmaNZs2Zpx44d+utf/6rS0lL99Kc/7YOzAwAAAAAAAPpev50c/OIXv6jt27drzZo1euuttxQOh1VYWKiJEyfqzDPP1Ne//nX5fN1L/i9+8QtNmDBBTz75pB599FGFQiGdeuqpuu666zRmzJh9fCYAAAAAAABA/9RvJwenTJmiKVOmdHv7uXPnpt0xmMrj8ejiiy/WxRdf3FvJAwAAAAAAAAa8frvmIAAAAAAAAIB9i8lBAAAAAAAAIEcxOQgAAAAAAADkKCYHAQAAAAAAgBzF5CAAAAAAAACQo5gcBAAAAAAAAHIUk4MAAAAAAABAjmJyEAAAAAAAAMhRTA4CAAAAAAAAOYrJQQAAAAAAACBHMTkIAAAAAAAA5CgmBwEAAAAAAIAcxeQgAAAAAAAAkKOYHAQAAAAAAAByFJODAAAAAAAAQI5ichAAAAAAAADIUUwOAgAAAAAAADmKyUEAAAAAAAAgRzE5CAAAAAAAAOQoJgcBAAAAAACAHMXkIAAAAAAAAJCjmBwEAAAAAAAAchSTgwAAAAAAAECOYnIQAAAAAAAAyFFMDgIAAAAAAAA5islBAAAAAAAAIEcxOQgAAAAAAADkKN/efLilpUXr1q3Tjh07FI1GO93uq1/96t4cBgAAAAAAAMA+0OPJwbvvvlsPP/ywWltbO93GcRxZlsXkIAAAAAAAANAP9Why8Pe//73uv/9+eb1ezZo1SwceeKAKCgp6NWHbtm3Tc889p8WLF+uTTz7Rzp07VVJSomnTpunyyy/X1KlTu7WfZcuW6eKLL+70/Xnz5mnOnDm9lWwAAAAAAABgwOjR5GBVVZXy8vI0f/58TZ48ubfTJEl67LHH9Pvf/15jxozR8ccfr/LyclVXV+ull17SSy+9pLvuuktnnHFGt/c3Y8YMzZgxo8PrkyZN6s1kAwAAAAAAAANGjyYHt2zZomOOOWafTQxK0uGHH67HHnusw4TeypUrdemll+rnP/+5Zs+erUAg0K39zZgxQ3Pnzt0XSQUAAAAAAAAGpB49rbiiokL5+fm9nZY0X/rSl7Le6Xf00Udr5syZqq+v1wcffLBP0wAAAAAAAAAMZj26c/CMM87Qn/70J4XDYYVCod5O0275fL60f7tjw4YNevjhhxWJRDR8+HAde+yxGj58+L5KIgAAAAAAANDv9WhycO7cuVqzZo2+//3v6xe/+IXGjh3b2+nq1ObNm7V06VJVVFRowoQJ3f7cokWLtGjRosTfPp9PF154oW688UZ5vd59kVQAAAAAAACgX+vW5GC2p/3atq3ly5frjDPO0MiRI1VZWSnLsjpsZ1mWHnnkkb1PqaRYLKYbb7xR0WhUP/rRj7o1qVdeXq7rr79eJ510kkaNGqWWlhatWbNGd911lx5++GFZlqWbbrqpV9IHAAAAAAAADCSW4zjO7jY65JBDen4Ay9J7773X48+7bNvWDTfcoEWLFum8887Tbbfdtlf727Fjh84++2w1NDRo8eLFGjJkSJfH9nh6tDwjAPQ5YhiAgY44BmAgI4YB6O+6defgyy+/vK/T0SXbtvXjH/9YixYt0tlnn61bb711r/dZUVGhU045RVVVVVq7dq1OPvnkTretr6/fo32XlZWptrZ2b5OIAYZyzz17WuZlZWX7MDWd29MYJlGfcxFlnpsGaxyjPucmyj339KTM+yKOMRZDd1Dmuam/jMW6NTk4atSofXLw7rBtWzfffLOefvppnXXWWbrjjjt67VsXN1NbWlp6ZX8AAAAAAADAQNKjWbZ77723W3cTvvLKK7r33nt7cghJ6RODZ5xxhn7961/36sND1q5dK6lvJz8BAAAAAACAvtLjycGXXnppt9u98soruu+++3pyiMRPiZ9++mmdfvrpuvPOO7ucGKypqdHHH3+smpqatNfXrVuXdftHHnlEy5Yt07hx4zRlypQepREAAAAAAAAYyLr1s+KeisfjPf4J8H333aeFCxcqFApp3Lhxuv/++ztsM3v2bE2aNEmSNH/+fN177726+uqrNXfu3MQ211xzjXw+nw477DANHz5cLS0tWrt2rd59910VFxfvdtIRAAAAAAAAGKz26eTgRx99pOLi4h59dtOmTZKkcDisBx54IOs2o0aNSkwOdub888/XkiVLtGLFCtXV1cnj8WjkyJG65JJLdNlll6mysrJH6QMAAAAAAAAGOstxHKc7G958882J/1+4cKHGjh2radOmZd02Ho/r008/1bp16zR79mzdc889vZPaPtKTJ2DxlKHcQ7nnnv7yZKnd6Um9pD7nHso8Nw3WOEZ9zk2Ue+4ZKE8rZiyG7qDMc1N/GYt1+87BhQsXJv7fsixVV1erurq6y89MnDhRN954Y89TBwAAAAAAAGCf6fbk4KOPPipJchxHl1xyiU488URdccUVWbf1+/0aNmwYTwEGAAAAAAAA+rFuTw7OmDEj8f9f+9rXdNRRR6W9BgAAAAAAAGBg6dEDSebNm9fb6QAAAAAAAACwn+3V04qj0aj+9re/aeXKldq+fbskadiwYTrqqKN02mmnKRgM9koiAQAAAAAAAPS+Hk8OLl26VDfddJN27NihzAceP/nkk7rzzjt1xx136Pjjj9/rRAIAAAAAAADofT2aHFy7dq2++93vKhaLaerUqTrzzDMTDx/ZvHmznn32Wb311lv63ve+p8cff1xTp07t1UQDAAAAAAAA2Hs9mhy8++671dbWpp///Oc6//zzO7x/0UUXacGCBbrlllv029/+Vn/4wx/2OqEAAAAAAAAAepenJx9au3atDjvssKwTg65vfOMbmjJlit56662epg0AAAAAAADAPtSjyUGPx6MxY8bsdrsxY8bIsqyeHAIAAAAAAADAPtajycHDDz9c69ev3+1269ev1+GHH96TQwAAAAAAAADYx3o0OXjttddqw4YN+u1vfyvbtju87ziOfvvb32rDhg269tpr9zqRAAAAAAAAAHpfjx5I8sknn+hrX/ua7r//fj3zzDP60pe+lHha8aZNm/Tiiy9q06ZNOvfcc/Xpp5/q008/Tfv8V7/61b1OOAAAAAAAAIC906PJwZtuukmWZclxHG3atEkPPfRQYm1Bx3ES21VVVamqqirxt+M4siyLyUEAAAAAAACgH+jR5OAPfvADHjQCAAAAAAAADHA9mhycO3dub6cDAAAAAAAAwH7WoweSAAAAAAAAABj4enTnYKr3339fb7/9tmpra3XwwQfrlFNOkSRFo1FFo1EVFhbudSIBAAAAAAAA9L4e3zn4ySef6Pzzz9fXvvY13XLLLfqP//gPvfTSS4n3//KXv2j69OlavHhxryQUAAAAAAAAQO/q0eTgli1bdOGFF+qtt97SSSedpBtuuCHtKcWS9OUvf1l+v18vvPBCryQUAAAAAAAAQO/q0c+K77vvPtXW1ur222/XOeecI0n69a9/nbZNKBTSpEmTtHbt2r1PJQAAAAAAAIBe16M7B1977TVNnDgxMTHYmVGjRmn79u09ShgAAAAAAACAfatHk4O7du3SgQceuNvt2tra1NLS0pNDAAAAAAAAANjHejQ5WFpaqi1btux2u08//VQVFRU9OQQAAAAAAACAfaxHk4PTpk3TP/7xD7333nudbrN8+XJ9+OGHmjFjRo8TBwAAAAAAAGDf6dHk4He+8x05jqOrrrpKr776quLxeNr7b7zxhm688Ub5fD5dcsklvZJQAAAAAAAAAL2rR08rnjp1qn7yk5/ol7/8pb73ve8pLy9PlmXphRde0EsvvaSmpiZZlqVbbrlFhxxySG+nGQAAAAAAAEAv6NGdg5J0wQUXaP78+TrppJNkWZYcx1Fzc7Oi0ahOOOEEPfbYY/rGN77Rm2kFAAAAAAAA0It6dOeg64gjjtB//ud/ynEc1dbWyrZtlZWVyev19lb6AAAAAAAAAOwjezU56LIsS+Xl5b2xKwAAAAAAAAD7Sa9MDu4L27Zt03PPPafFixfrk08+0c6dO1VSUqJp06bp8ssv19SpU7u9L9u2NX/+fD355JOqrq5WKBTScccdp+uuu06jR4/eh2cBAAAAAAAA9F/dmhw85ZRTenwAy7L00ksv7fHnHnvsMf3+97/XmDFjdPzxx6u8vFzV1dV66aWX9NJLL+muu+7SGWec0a19/exnP1NVVZXGjx+viy66SNu3b9dzzz2n119/XQsWLNC4ceP2OH0AAAAAAADAQNetycFNmzbt8Y7dh5T01OGHH67HHntMM2bMSHt95cqVuvTSS/Xzn/9cs2fPViAQ6HI/b775pqqqqjR9+nT98Y9/TGx/1lln6corr9Rtt92mP/zhDz1OJwAAAAAAADBQdWty8P333+/w2u23366nnnpKF1xwgc4880wdcMABksxE4rPPPqv58+fra1/7mn7605/2KGFf+tKXsr5+9NFHa+bMmVqyZIk++OADTZkypcv9VFVVSZKuvfbatInEWbNmacaMGVqyZIk2b96skSNH9iidmbZssbXw6TqtXm3L75dCRVLNTikSkTwe858tyeeRfH7J75faolK0TXIkBbxSJCo5juTzmW3y8iS7TWpplbweKdZm3nf348i8JjuZDvdY8kj5AcmyzNuWIxUUSrW1UiwmBQOS1y9FWtP3abU/x9q2zb49Sm7bEpbicZOuIUOkaERqiST3kXlMd3s3XT6fSXPcNscK5kmFhVI8JjU0Sq1Rc7zUvLI8Ji1x27zX1v6vu69Ym9ku4E/mT5tt8iT1mJZHKi6UYlEp3Jrcnzwm7yUpVCC1tkhNLcl9xu3kvttsc362nSy3/IAUCO5SQ2PHsvD5pPw8s32bnczXzHPJzF+3DDPfj7SavLU8ktdrysHXXl7xuNTaml6fMvMjM12xtmSeh/KlYL5Jg9fTfhyvOWa0rWPeu3U40ppebqnH8njMfkOFUrjJ5Gvq5z2e5Dm59cI9jlt+wYCUX2DaQWNTsj6l1vU2u2MdsJU9j1Pre9CXnu9u3hQWmP9vDifz22XbUigknfzFsKYebquoqMcPfu+XqqttzX+iVv9YZyvelmyLHnXM02x1NDUvUwUCUkGB1NRs2mCq1Haa2YZcPp8p39QYmbl9Zpx130+NKan1a3cy20nqcTLrU2qdTc2X1HS66XBjTiSaHh+7m8bUY6SmKTW+hoKmPTc0JPM783PuZyyP5PfvUlNzx/gbzJPy85N9VbYYldpWU+NUavkEAlJxsYmv4Uj28nH7Dze+ZosH2c67qzLvblln5q+bL24/2NnxO/us0172FUOllhaprq779STbsbr62419u6tPmccpK5bOOqtFkw5xFAxae5ZJ/Vhjo61ly5v18v/Zqq1L7+c76xNTpdblzvrkrupBd+pJV9uljhtS20FqHc8cQ3SnXmZr+6ljUXeM4vUm+3a37+8qLmWL313F9Gz5my2tmf1OtjSk5pXHJxUW7dLOHekxL7Utx23T7wfyOo5tujqv1DFVZmzLzN/MvrGzfEjt0zLHNdn6jdR8SB0bZ8aMzFiYLY2yk8fvqk1k+4zbJ0bbOh5nd/1Ytv25463GpvTrH4/H1A/HyX7+7rj+sEMbdcJxtsaO9Wgwqa62taCqTu++Yysal4aWm3yoq5XijlReln4NZlnJ/jrcZPpZN7/cOu+2c7ctdHYdKXU+5soWj1KvFd3juLHLveZ1r2cbm0xaU/fdWf+fWc+zjU8CgfRxjit1zOmOZd1juNdzbpsOBqSiEnNS8bgUjSavNzLHVqn1OrM9pOZ3d8dCmTEu6N+lSKzjPjKvQ9xY5vWbc8+cI8gs88zPunnr9SpxLebzmustyzHXk45l+gKPpLyQFG011+rZrrO7arOpcTFb/Ar4TFmFW9rLymeu193r8Wzj3NSYE8qXSstNHrSETRptSQX56dccXcXAbNehAZ8pK3fewK2DRYVSUakUbpSaW9LjsFvvIhFJjtlfZtmkjknitmnLM2aENflQW0OG9G0c69Gag48++qj+53/+R0888USHybmJEydq4sSJOvXUU/XNb35TBxxwgC699NLeSGuCz+dL+7cry5YtUygU0rRp0zq8d+KJJ2r58uVavny5vvrVr+51up75i6177pNaW7NcFQMYdF54sUWjRko//YmtKZMHx6D0sfm2/vCQ1JY6kgAwaD3/YliTD5VuvUWqHD7wJwj/8Y6teb+SPvusta+TAmA/eO21qP7rj9Jll9q66ILBMxb748NSLMY1JZALnv5LiyoqpP93ja1ZJ/ZdHOvRkRcsWKCZM2d2edfelClTNHPmzMSde71l8+bNWrp0qSoqKjRhwoQutw2Hw9qxY4cOOOAAeb3eDu+PHTtWklRdXb3X6dqyxdYjj5lZdgC5Y9Nm6fH55k6Vga662tZj86W2tr5OCYD96Z13pT//xVEk0vPlYPqDxkZbj8+XPvusr1MCYH+KxaQnFphxzEBXXW3r8f825wQgd+zYIT02X9q1q+/iWI/uHPz88881ceLE3W5XUlKilStX9uQQWcViMd14442KRqP60Y9+lHXCL1VjY6MkqbCwMOv77uvudp0pKSmRx9P1POq6d1rV2Njc5TYABqft2z1ynCKVlfXPB8B3J4ZJ0pvLWhSNhvdDigD0N5997lNeXpFKS/vnnTfdiWN1dW3avr1+P6UIQH/S2ipt316gI47I6+ukZMVYDMDu1NVJLS1FKivr+rka+0qPrmRLSkq0YsUKRSIRBYPBrNtEIhGtWLFCxcXFe5VAl23buummm7RixQqdd955vfIz4O6qr9/9QHNIua2iIrOuEIDcMmyYLcuqV21t14O+srKy/ZSidN2JYZI0otJWIMCdg0AuGjO6Ta2tdaqt7fqnxf05jnm9toYNkz78aD8kCEC/kpcnDRvWrNra3V+M9UUc29OxGHcOArmntFTKz2/ss2vKHn09fOqpp2rHjh265pprtHHjxg7vb9y4Uddee6127typU089da8Tadu2fvzjH2vRokU6++yzdeutt3brc0VFRZKkpqamrO+7r7vb7Y0RIzy65CLTMQHIHaNGShdeoEHxUJKxYz266AKzQC+A3DH5UOnsr1gD/qEkRUUeXXiBNGZMX6cEwP7k90vf/IYGxUNJxo716MJvmXMCkDsqKqSLLlCfPpSkR5eA1157rd588029+uqrWrJkiQ477LDE0343b96sd955R21tbTrooIN07bXX7lUCbdvWzTffrKefflpnnXWW7rjjjm7dki1JoVBIFRUV2rhxo+LxeIefIbtrDbprD+6tf/6KRzOOtrV8hU9r1rTJ55cKQlJNjXmSlMdjnlZjO+ZpQH6/5PeZu3QiMfMEm4DfPN3Gdsx7fr+UFzRPsmltkTxe802S4yT34yj5mss9lizzFCWPZZ6kY1lSfkiqrzPHzAuafbRG0vfpZnHcNvu2rOS24RapLW72W15u3m9tSe4j85ju9m66/O4Tn+Jmf3l55qk+bXGpqcHklWWl55XHY9ISj5v9tsXNNu6+YjHztKVgIJkXbfHk05Lc7TyWVFRktgm3JPen9qdmWTL509pqnm7k7tOOS/6A1BaTYu3n7tjJcsvPkwLB9qdlZZSFz2eeouTYyXzIdi6Z+euWYeb7re3542l/wlReUPK2t+R43DxtKrU+ZeZHZrraYsk8LwhJwaA5f4/HpM+yTJ2MxDrmvVuHI5H0cks9lsdj9hsKSeGwydfUz3us5Dm59cI9jlt+eUGT1rgtNTYm8zG1rrfFO9YB28mex6n1PRhIz3c3b4raVyNobk7mtytum3P64hfzdfhhLYNiYtB10QUe/dMJtpa+4dU/3omrrS3ZFj1WxzzNVkdT8zJVIGCeStjUbJ4Elyq1nWa2IZfPZ8orNUZmbp8ZZ933U2NKav3ancx2knqczPqUWmdT8yU1nW463JgTiaTHx+6mMfUYqWlKja+hfNOeGxqT+Z35OfczHo9pj+5TBFPzMC/P/NfW1v7EtiwxKrWtpsap1PIJBKSSYhNfwi3Zy8ftP9z4mi0eZDvvrsq8u2Wdmb9uvrj9YGfH7+yzjsy/Q4aYc048rbgb9STbsbr62419u6tPmccpK5XOPCOkQya2DPiJQdeUyR797j5b77+fp1f+3qqa2vR+vrM+MVVqXe6sT+6qHnSnnnS1Xeq4IbUdpNbxzDFEd+pltrafOhZ1647Xm+zb3b6/q7iULX53FdOz5W+2tGb2O9nSkJpXHq9UWCjt2Jke81Lbcjxu2orbXlLHNl2dV+qYKjO2ZeZvZt/YWT6k9mmZ45ps/UZqPqSOjTNjRmfXB5l1wD1+V20i22fcPjES63ic3fVj2fbnjrcaG9OvfzyWqR+2k/383XH95MkBHXdMdFBMDLrcsdjyFT69+16bYjHTl3g8pi+xbXOHUeo1mGUl++tw2OSdm19unXfbudsWOruOlDofc2WLR6nXiu5x3NjlXvO617ONjZKs9H131v9n1vNs4xO/P32c40odc7pjWfcY7vWc26aDQXN96Mhc80WjyeuNzLFVar3ObA+p+d3dsVBmjMsLms9l7iPzOsSNZb72uYHMOYLMMs/8rJu3Xm/yWsw9H1nmetIdx3g85snQkajU2JD9OrurNpsaF7PFL7/flFVLS3tZeU15uNfj2ca5qTGnIGTag7e9T2lqMMcJhdKvObqKgdmuQ/1+U1buvIG7XXGR+a85bNpaahx26120/Uncdrxj2aSOSey4VFYuzZgR0qSJ4T5/WrHlOHs6ZDYaGhp011136ZlnnlFrxlM48vLydPbZZ+v6669XSUlJjxOXOjF4xhln6N/+7d92u85gph/+8Id69tln9fjjj2v69Olp71100UVavny5XnnlFY0aNarTfdTW1u7RMcvKyvb4Mxj4KPfcs6dl3lc/x+tJvaQ+5x7KPDcN1jhGfc5NlHvu6UmZ90UcYyyG7qDMc1N/GYv1+MdjxcXFuvXWW3XTTTfpnXfe0fbt2yVJFRUVmjx5skKh0F4lzP0p8dNPP63TTz9dd955Z5cTgzU1NaqtrVVZWZnKy8sTr5933nl69tlndffdd+uPf/yjAgGzuOOrr76q5cuX64QTTuhyYhAAAAAAAAAYrPZ6Zan8/HwdffTR3d6+qqpKq1ev1rx587rc7r777tPChQsVCoU0btw43X///R22mT17tiZNmiRJmj9/vu69915dffXVmjt3bmKbY445Rueee66qqqo0Z84czZo1Szt27NBf//pXlZaW6qc//Wm30w4AAAAAAAAMJvt92fnVq1fr6aef3u3k4KZNmyRJ4XBYDzzwQNZtRo0alZgc7MovfvELTZgwQU8++aQeffRRhUIhnXrqqbruuus0hlWrAQAAAAAAkKP67TMp77jjDt1xxx3d3n7u3Llpdwym8ng8uvjii3XxxRf3VvIAAAAAAACAAW/wPNYJAAAAAAAAwB5hchAAAAAAAADIUUwOAgAAAAAAADmKyUEAAAAAAAAgRzE5CAAAAAAAAOQoJgcBAAAAAACAHMXkIAAAAAAAAJCj9vvkoOM4+/uQAAAAAAAAALLw7e8DXnnllZozZ87+PiwAAAAAAACADHs1OfjRRx/pySef1Ntvv63a2lqdcsopuvHGGyVJq1ev1rp163T22WertLQ08ZmDDjpIBx100F4lGgAAAAAAAMDe6/Hk4EMPPaS77rpLbW1tkiTLslRbW5u2zbx58xQIBHT++efvXSoBAAAAAAAA9LoerTn497//Xb/61a9UWVmpe++9V0uXLu2wluC0adNUXl6ul19+uVcSCgAAAAAAAKB39ejOwYceekj5+fl66KGHNHr06E63O+SQQ/Tpp5/2OHEAAAAAAAAA9p0e3Tn4zjvv6IgjjuhyYlCSysrKtHPnzh4lDAAAAAAAAMC+1aPJwVgspoKCgt1uV1NTI6/X25NDAAAAAAAAANjHejQ5eMABB+j999/vcptoNKoPPvhA48aN68khAAAAAAAAAOxjPZocPPnkk7Vp0yY99NBDnW7zX//1X6qpqdGXvvSlHicOAAAAAAAAwL7ToweSXH755frLX/6iX//611q7dq1OPfVUSdKuXbv04osv6sUXX9Rf/vIXHXDAAbrgggt6NcEAAAAAAAAAekePJgdLSkr00EMP6ZprrtHzzz+vv/3tb5Kk1157Ta+99pocx9HBBx+s++67T4WFhb2aYAAAAAAAAAC9o0eTg5J04IEH6plnntErr7yi119/XZs2bZJt26qsrNRxxx2n0047jYeRAAAAAAAAAP1YjycHJcnj8Wj27NmaPXt2b6UHAAAAAAAAwH7SoweSAAAAAAAAABj49mpycMmSJfrBD36gE088UYcddphuvvnmxHuvvfaa5s2bp23btu11IgEAAAAAAAD0vh7/rPj222/X/Pnz5TiOQqGQ2tra0t6vqKjQI488ohEjRujSSy/d23QCAAAAAAAA6GU9unPw6aef1uOPP67Jkydr4cKFWr16dYdtDjnkEI0YMUKvvPLKXicSAAAAAAAAQO/r0Z2DTzzxhIqLi/Xggw+qvLy80+0mTpyo9evX9zhxAAAAAAAAAPadHt05uH79eh155JFdTgxKUmFhoXbu3NmjhAEAAAAAAADYt3r8QBLLsna7zfbt25WXl9fTQwAAAAAAAADYh3o0OThu3Di98847isVinW7T1NSk999/XwcffHCPEwcAAAAAAABg3+nR5ODpp5+uHTt26K677up0m3//939XY2OjzjzzzB4nDgAAAAAAAMC+06MHklxyySV69tln9cgjj2jNmjU65ZRTJEmff/65Hn74Yb344otatWqVDj30UJ177rk9TtwzzzyjVatWad26dVq/fr1isZjmzZunOXPmdHsfy5Yt08UXX9zp+3u6PwAAAAAAAGCw6NHkYF5enh5++GHddNNNWrx4sd5++21J0sqVK7Vy5UpJ0vHHH68777xTgUCgx4m7++67tWnTJpWVlWnYsGHatGlTj/c1Y8YMzZgxo8PrkyZN6vE+AQAAAAAAgIGsR5ODklReXq4HH3xQ77//vpYsWaJNmzbJtm1VVlbq+OOP1+GHH77Xibv99ts1duxYjRo1Sg8++GCXP2PenRkzZmju3Ll7nSYAAAAAAABgsOjR5ODVV1+tiooK3XLLLTrkkEN0yCGH9Ha6JEnHHXfcPtkvAAAAAAAAgB5ODr766quaPXt2b6dln9qwYYMefvhhRSIRDR8+XMcee6yGDx/e18kCAAAAAAAA+kyPJgcPOOAAtbS09HZa9qlFixZp0aJFib99Pp8uvPBC3XjjjfJ6vX2YMgAAAAAAAKBv9Ghy8Mwzz9Qf//hH7dixQxUVFb2dpl5VXl6u66+/XieddJJGjRqllpYWrVmzRnfddZcefvhhWZalm266qct9lJSUyOPx7NFxy8rK9ibZGKAo99wzEMq8JzFMGhjnht5FmeemgVDujMXQXZR77hkIZc5YDN1Fmeem/lDuPZoc/O53v6u3335bF154oX70ox/pi1/8ovx+f2+nrVeMHz9e48ePT/wdCoU0e/ZsTZ06VWeffbYee+wxXXHFFRoyZEin+6ivr9+jY5aVlam2trbHacbARLnnnj0t874K+nsawyTqcy6izHPTYI1j1OfcRLnnnp6UeV/EMcZi6A7KPDf1l7FYjyYHTz/9dDmOoy1btuiaa66RZVkqLy9XMBjssK1lWXrppZf2OqG9raKiQqeccoqqqqq0du1anXzyyX2dJAAAAAAAAGC/6tHk4KZNm9L+dhxHO3fu7JUE7U/ujOtAWz8RAAAAAAAA6A09mhx8//33ezsdfWLt2rWSpFGjRvVxSgAAAAAAAID9b89XRZV0880363//9393u93ChQv14x//uCeH2GM1NTX6+OOPVVNTk/b6unXrsm7/yCOPaNmyZRo3bpymTJmyP5IIAAAAAAAA9Cs9unNw4cKFkqSvf/3rXW63atUqLVy4UL/85S97chhVVVVp1apVkqT169cnXlu+fLkk6aijjtK5554rSZo/f77uvfdeXX311Zo7d25iH9dcc418Pp8OO+wwDR8+XC0tLVq7dq3effddFRcX684775TX6+1R+gAAAAAAAICBrEeTg90Vi8X2auLNnVxMtXr1aq1evTrxtzs52Jnzzz9fS5Ys0YoVK1RXVyePx6ORI0fqkksu0WWXXabKysoepw8AAAAAAAAYyPbZ5KDjOHr33Xf36jHLd9xxh+64445ubTt37ty0OwZdV155pa688soepwEAAAAAAAAYrLo9OXjxxRen/f3aa691eM0Vj8f12WefaefOnfrnf/7nvUshAAAAAAAAgH2i25OD7jp/kmRZlnbu3KmdO3d2vmOfT1/84hd144037l0KAQAAAAAAAOwT3Z4cfPnllyWZnwvPnj1bp512WqcTf36/X2VlZfL7/b2TSgAAAAAAAAC9rtuTg6NGjUr8/9VXX61JkyalvQYAAAAAAABgYOnRA0muvvrq3k4HAAAAAAAAgP3M09cJAAAAAAAAANA3mBwEAAAAAAAAchSTgwAAAAAAAECOYnIQAAAAAAAAyFFMDgIAAAAAAAA5islBAAAAAAAAIEcxOQgAAAAAAADkKCYHAQAAAAAAgBzF5CAAAAAAAACQo5gcBAAAAAAAAHIUk4MAAAAAAABAjmJyEAAAAAAAAMhRTA4CAAAAAAAAOYrJQQAAAAAAACBHMTkIAAAAAAAA5CgmBwEAAAAAAIAcxeQgAAAAAAAAkKOYHAQAAAAAAAByFJODAAAAAAAAQI5ichAAAAAAAADIUUwOAgAAAAAAADmKyUEAAAAAAAAgRzE5CAAAAAAAAOQoJgcBAAAAAACAHOXr6wR05ZlnntGqVau0bt06rV+/XrFYTPPmzdOcOXP2aD+2bWv+/Pl68sknVV1drVAopOOOO07XXXedRo8evY9SDwAAAAAAAPRv/Xpy8O6779amTZtUVlamYcOGadOmTT3az89+9jNVVVVp/Pjxuuiii7R9+3Y999xzev3117VgwQKNGzeudxMOAAAAAAAADAD9+mfFt99+u1555RW9+eabOv/883u0jzfffFNVVVWaPn26nnrqKd1www268847dd9996murk633XZbL6caAAAAAAAAGBj69Z2Dxx133F7vo6qqSpJ07bXXKhAIJF6fNWuWZsyYoSVLlmjz5s0aOXLkXh8LAAAAAAAAGEj69Z2DvWHZsmUKhUKaNm1ah/dOPPFESdLy5cv3d7IA9BORiKO6OkeRiNPXSUE/lWt1JNfOF+httCEA6B7iJdB/9Os7B/dWOBzWjh07NGHCBHm93g7vjx07VpJUXV29v5MGoB/Yus3R2rWOmpqlwgJp6lSpcrjV18lCP5JrdSTXzhfobbQhAOge4iXQvwzqOwcbGxslSYWFhVnfd193twOQOyIRMyCpb5RKS6X6RmntWr65RFKu1ZFcO1+gt9GGAKB7iJdA/zOo7xzsLSUlJfJ49mwetaysbB+lBv0Z5T5w1Nbasp2Ixo72KBi0VFjgqKbWVl5eUKWl3W/vA6HMexLDpIFxbvtSb9WRgcKcbyhnzhdJA6GtD4SxWK7FjP5qINRn9K6BUOaMxdIRLzs3WMscXesP5T6oJweLiookSU1NTVnfd193t+tMfX39Hh23rKxMtbW1e/QZDHyU+8ASiTjyWI6qP5eGlEu7aqSSIqm1tUW1td37ScOelnlfBf09jWES9VnqnToykIRCpfJY4Zw5XxiDNY71RQzLtZjRH9F35Z6elHlfxDHGYumIl9kN5jJH5/rLWGxQT8uHQiFVVFRo48aNisfjHd531xp01x4EkDuCQUtTp1oqKZLq6syAZOpUS8Fg7g5IkC7X6kiunS/Q22hDANA9xEug/xnUdw5K0owZM/Tss89q9erVmj59etp7r732miR1eB1AbqgcbqlsltTaKuXliQEJOsi1OpJr5wv0NtoQAHQP8RLoXwbNnYM1NTX6+OOPVVNTk/b6eeedJ0m6++67FY1GE6+/+uqrWr58uU444QSNGjVqv6YVQP8RDFoqKeGbSnQu1+pIrp0v0NtoQwDQPcRLoP/o13cOVlVVadWqVZKk9evXJ15bvny5JOmoo47SueeeK0maP3++7r33Xl199dWaO3duYh/HHHOMzj33XFVVVWnOnDmaNWuWduzYob/+9a8qLS3VT3/60/18VgAAAAAAAED/0K8nB1etWqWFCxemvbZ69WqtXr068bc7OdiVX/ziF5owYYKefPJJPfroowqFQjr11FN13XXXacyYMb2ebgAAAAAAAGAgsBzHcfo6Ef1dT56AxVOGcg/lnnv6y5Oldqcn9ZL6nHso89w0WOMY9Tk3Ue65Z6A8rZixGLqDMs9N/WUsNmjWHAQAAAAAAACwZ7hzEAAAAAAAAMhR3DkIAAAAAAAA5CgmBwEAAAAAAIAcxeQgAAAAAAAAkKOYHAQAAAAAAABylK+vEzBYvP3227rnnnu0Zs0atbW1acKECbr00kt1xhln9HXSkOKZZ57RqlWrtG7dOq1fv16xWEzz5s3TnDlzsm7f1NSke+65Ry+88IJ27NihYcOG6bTTTtPVV1+tgoKCDtvbtq358+frySefVHV1tUKhkI477jhdd911Gj16dNZjvPbaa/rd736nd955R5ZlafLkybrqqqt07LHH9uq556pt27bpueee0+LFi/XJJ59o586dKikp0bRp03T55Zdr6tSpHT6Tq+VOHBsYiGO5hRjWfcSwgYEYlnuIY91HHBsYiGO5JxfiGE8r7gVvvvmmLr/8cgUCAZ155pkqKCjQCy+8oE2bNulf/uVfdNlll/V1EtHu5JNP1qZNm1RWVqZQKKRNmzZ1GsjD4bC+9a1v6b333tMJJ5ygSZMm6b333tOSJUs0ZcoUzZ8/X8FgMO0zP/3pT1VVVaXx48dr1qxZ2r59u5577jkVFBRowYIFGjduXNr2zzzzjG688UaVl5cnOv2//vWvqq2t1X/8x3/o9NNP32d5kSv+7d/+Tb///e81ZswYzZgxQ+Xl5aqurtZLL70kx3F01113pQ24crXciWMDB3EstxDDuocYNnAQw3IPcax7iGMDB3Es9+REHHOwV2KxmDN79mznsMMOc959993E6w0NDc6XvvQlZ/Lkyc7GjRv7MIVI9frrryfK43e/+50zYcIE53//93+zbnv33Xc7EyZMcO6888601++8805nwoQJzgMPPJD2+htvvOFMmDDBueCCC5xIJJJ4/e9//7szYcIE57LLLkvbvq6uzjn66KOdmTNnOlu2bEm8vmXLFmfmzJnOzJkzncbGxr06XzjO3/72N2fZsmUdXl+xYoUzefJkZ/r06WnllYvlThwbWIhjuYUYtnvEsIGFGJZ7iGO7RxwbWIhjuScX4hhrDu6lN998U5999pnOOussTZo0KfF6UVGRvve97ykWi2nhwoV9mEKkOu644zRq1Kjdbuc4jqqqqhQKhXTVVVelvXfVVVcpFAqpqqoq7XX372uvvVaBQCDx+qxZszRjxgwtWbJEmzdvTrz+/PPPq6GhQRdeeKEqKysTr1dWVurCCy9UbW2tXnrppR6dJ5K+9KUvacaMGR1eP/roozVz5kzV19frgw8+kJS75U4cG1iIY7mFGLZ7xLCBhRiWe4hju0ccG1iIY7knF+IYk4N7afny5ZKkE044ocN77msrVqzYr2nC3tuwYYO2b9+uadOmKRQKpb0XCoU0bdo0ff7559qyZUvi9WXLliXey3TiiSdKStaX1P/vqu6kbo/e5/P50v7N1XInjg1OuVqfcwkxrPtpIoYNPLlan3MNcaz7aSKODTy5Wp9zzWCJY0wO7qUNGzZIksaOHdvhvYqKCoVCIVVXV+/nVGFvuWWW+bt+l/u6W/7hcFg7duzQAQccIK/X22F7t36k1oWu6k627dG7Nm/erKVLl6qiokITJkyQlLvlThwbnHK1PucKYlgSMWxwytX6nEuIY0nEscEpV+tzLhlMcYzJwb3U1NQkydzynU1hYaEaGxv3Z5LQC9wyKywszPq++7pb/t3dPrUudFV3sm2P3hOLxXTjjTcqGo3qRz/6USII52q5E8cGp1ytz7mAGJaOGDY45Wp9zhXEsXTEscEpV+tzrhhscYzJQQA5xbZt3XTTTVqxYoXOO+88ffWrX+3rJAFAtxHDAAx0xDEAA91gjGNMDu6l3c3KNjU1dfoNEPovt8zcGflM7utu+Xd3+9S60FXd2d23h+gZ27b14x//WIsWLdLZZ5+tW2+9Ne39XC134tjglKv1eTAjhmVHDBuccrU+D3bEseyIY4NTrtbnwW6wxjEmB/eS+3vxbL/n3rFjh8LhcNbfgaN/c8vM/S1/Jvd1t/xDoZAqKiq0ceNGxePxDtu79SO1LnRVd7Jtj71j27ZuvvlmLVy4UGeddZbuuOMOeTzpITBXy504Njjlan0erIhhnSOGDU65Wp8HM+JY54hjg1Ou1ufBbDDHMSYH99L06dMlSUuWLOnwnvuauw0GjnHjxmnYsGFavXq1wuFw2nvhcFirV6/WAQccoBEjRiRenzFjRuK9TK+99pqk9LrQnbqT7XHp2HNuEH/66ad1xhln6Ne//nXWxV5ztdyJY4NTrtbnwYgY1jVi2OCUq/V5sCKOdY04Njjlan0erAZ7HGNycC8de+yxGj16tBYtWqT33nsv8XpjY6MeeOAB+f3+QfH781xjWZbOPfdchcNh/ed//mfae//5n/+pcDis8847L+119++7775b0Wg08fqrr76q5cuX64QTTtCoUaMSr3/5y19WUVGRHn/8cW3dujXx+tatW/X444+rrKxMs2fP3henl1Pc276ffvppnX766brzzjuzBnEpd8udODY45Wp9HmyIYbtHDBuccrU+D0bEsd0jjg1OuVqfB6NciGOW4zjOHn0CHbz55pu6/PLLFQgEdOaZZ6qgoEAvvPCCNm3apH/5l3/RZZdd1tdJRLuqqiqtWrVKkrR+/Xq98847mjZtWuKW26OOOkrnnnuuJDOr/81vflPvv/++TjjhBB166KF69913tWTJEk2ZMkWPP/648vLy0vb/05/+VFVVVRo/frxmzZqlHTt26K9//asKCgr0P//zPzrwwAPTtn/mmWd04403qry8XGeccYYk6a9//atqa2v1m9/8Rl/+8pf3dZYMevfcc4/uvfdehUIhXXzxxfL5fB22mT17tiZNmiQpd8udODZwEMdyCzGse4hhAwcxLPcQx7qHODZwEMdyTy7EMSYHe8nbb7+t3/72t1qzZo3a2to0YcIEffvb304UEvqHm266SQsXLuz0/a997Wu64447En83Njbqnnvu0QsvvKCdO3eqoqJCp59+un7wgx9kfdS4bdt6/PHH9eSTT6q6ulqhUEjHHXecrrvuOo0ZMybrMRcvXqzf/e53evfddyVJhx12mL7//e/ruOOO28uzhbT7MpekefPmac6cOYm/c7XciWMDA3EstxDDuo8YNjAQw3IPcaz7iGMDA3Es9+RCHGNyEAAAAAAAAMhRrDkIAAAAAAAA5CgmBwEAAAAAAIAcxeQgAAAAAAAAkKOYHAQAAAAAAAByFJODAAAAAAAAQI5ichAAAAAAAADIUUwOAgAAAAAAADmKyUEAAAAAAAAgRzE5CAAAAAAAAOQoJgcBSffcc48mTpyop556qq+TAgB7jBgGYKAjjgEY6IhjGMiYHAQAAAAAAAByFJODAAAAAAAAQI5ichAAAAAAAADIUUwOIqe8/PLL+sY3vqGpU6dq5syZmjt3rj799NOs25588smaOHGiHMfRY489prPPPltTp07VP//zP0uSHMfRokWLdN111+m0007TEUccoSOPPFLnnHOO5s+fL9u20/b31FNPaeLEibrnnnvSXq+rq9MhhxyiiRMnqqqqKu299evXa+LEifrud7/bi7kAYKAihgEY6IhjAAY64hgGI19fJwDYX5544gn9/Oc/l2VZOvroo1VRUaG1a9fq3HPP1UknndTp52655RY99dRTmj59ur7whS8oFotJkqLRqK6//nqVlpbq4IMP1qGHHqq6ujqtWbNGv/jFL/SPf/xDd9xxR2I/06dPlyQtX748bf8rVqyQ4ziSpGXLluncc89NvLds2TJJ0owZM3onEwAMWMQwAAMdcQzAQEccw2DF5CBywqZNmzRv3jz5/X7df//9OvHEEyVJsVhMN998s/785z93+tkXXnhBCxcu1Pjx49Ne93q9uu+++zRr1iz5/f7E6zU1Nbriiiu0cOFCff3rX08E8NGjR2vUqFF66623FIlEFAwGJSWD9fjx47VixYq0Y7hBn0AO5DZiGICBjjgGYKAjjmEw42fFyAn/+7//q0gkojPPPDMRxCXJ7/frJz/5ifLz8zv97BVXXNEhiEuSz+fT7Nmz04K4JJWXl+v666+XZG45TzV9+nRFo1G99dZbideWL1+uL3zhCzrttNO0detWVVdXSzK3mK9YsUKFhYU69NBD9/icAQwexDAAAx1xDMBARxzDYMadg8gJK1eulCSdccYZHd4rKyvT8ccfr5deeinrZ08++eQu9/3ee+9pyZIl2rx5s1pbW+U4jpqbmyVJGzZsSNt2+vTpevrpp7V8+XLNnDlTdXV1Wr9+vc4///zENznLly/X2LFjtX79etXW1mrWrFnyer17esoABhFiGICBjjgGYKAjjmEwY3IQOWH79u2SpFGjRmV9v7PXJWnkyJFZX49Go7r55pu1aNGiTj/rBnTXzJkzJSVv7XbXhpg5c6aOOOIIBQKBxBoR3P4NwEUMAzDQEccADHTEMQxmTA4Cu+Gu45Dp4Ycf1qJFizRhwgTdcMMNmjx5soqLi+X3+/Xpp5/q9NNP7/CZ0aNHa8SIEYk1Ity1IWbOnKlgMKgjjjgisUYEgRxAbyCGARjoiGMABjriGPo71hxETqioqJBkFpHNZvPmzXu8zxdffFGS9O///u/6p3/6Jw0ZMiSxVsTnn3/e6edS14hYvny5Dj74YJWXl0syQXvr1q3asGGDVqxYoYKCAk2ePHmP0wZgcCGGARjoiGMABjriGAYzJgeRE44++mhJ0vPPP9/hvbq6Or3++ut7vM+GhgZJUmVlZYf3nnvuuU4/535r8+KLL2r9+vWJ28JT33v88cdVW1uro446irUhABDDAAx4xDEAAx1xDIMZk4PICXPmzFEgENBf/vIXLV26NPF6LBbTvHnzFA6H93if48aNkyQ98cQTaa8///zzeuaZZzr9nBu4n3zySTmOk3aLt7tGxJNPPilJiUfWA8htxDAAAx1xDMBARxzDYMbkIHLC6NGjddNNNykWi+k73/mOLrroIv3whz/Uaaedppdffllf+cpX9nifl19+ubxer+666y7NmTNH119/vb7+9a/r2muv1SWXXNLp58aMGaPKykpFIhFZlpUWyN01IiKRiCSlfQMEIHcRwwAMdMQxAAMdcQyDGZODyBkXXHCB7rvvPk2ZMkVvv/22lixZokMOOUQLFizQ2LFj93h/06dP13//93/rmGOO0caNG/V///d/8vv9uueee3TBBRfs9rOS0taGcLmBPRQKsTYEgARiGICBjjgGYKAjjmGwshzHcfo6EQAAAAAAAAD2P+4cBAAAAAAAAHIUk4MAAAAAAABAjmJyEAAAAAAAAMhRTA4CAAAAAAAAOYrJQQAAAAAAACBHMTkIAAAAAAAA5CgmBwEAAAAAAIAcxeQgAAAAAAAAkKOYHAQAAAAAAAByFJODAAAAAAAAQI5ichAAAAAAAADIUUwOAgAAAAAAADmKyUEAAAAAAAAgR/3/z8fmyQst/o4AAAAASUVORK5CYII=\n",
"text/plain": [
"