Skip to content

Instantly share code, notes, and snippets.

@vijayvd
Created May 7, 2015 03:20
Show Gist options
  • Save vijayvd/d5b10b27ba4f9b9da9b0 to your computer and use it in GitHub Desktop.
Save vijayvd/d5b10b27ba4f9b9da9b0 to your computer and use it in GitHub Desktop.
Problem with pykalman
{"nbformat_minor": 0, "cells": [{"execution_count": 1, "cell_type": "code", "source": "import pandas as pd\nimport numpy as np\nimport pykalman as kf\nimport scipy as sp\n%pylab inline", "outputs": [{"output_type": "stream", "name": "stdout", "text": "Populating the interactive namespace from numpy and matplotlib\n"}], "metadata": {"collapsed": false, "trusted": true}}, {"source": "### Convention", "cell_type": "markdown", "metadata": {}}, {"source": "$$\n\\begin{split}\nx_{t+1} \n&= \nAx_t + b_t + \\epsilon_t,\\ \\epsilon_t\\sim \\mathscr{N}(0, Q_t) \n\\\\\nz_t\n&=\nC_tx_t + d_t + \\delta_t,\\ \\delta_t\\sim\\mathscr{N}(0, R_t)\n\\end{split}\n$$", "cell_type": "markdown", "metadata": {}}, {"execution_count": 2, "cell_type": "code", "source": "n_obs, n_lat = 1, 2\nT = 1000", "outputs": [], "metadata": {"collapsed": true, "trusted": true}}, {"execution_count": 3, "cell_type": "code", "source": "x = np.zeros((T,n_lat))\nz = np.zeros((T,n_obs))", "outputs": [], "metadata": {"collapsed": true, "trusted": true}}, {"source": "### Initial conditions", "cell_type": "markdown", "metadata": {}}, {"execution_count": 4, "cell_type": "code", "source": "x[0] = [1, 2]\nz[0] = [3]", "outputs": [], "metadata": {"collapsed": false, "trusted": true}}, {"source": "### Dynamics of the system", "cell_type": "markdown", "metadata": {}}, {"execution_count": 5, "cell_type": "code", "source": "A = np.array([[0.5, -0.25], [-0.25, 0.5]])\nC = np.array([[1., 1.]])", "outputs": [], "metadata": {"collapsed": true, "trusted": true}}, {"source": "### System noise and Observation noise", "cell_type": "markdown", "metadata": {"collapsed": false}}, {"execution_count": 6, "cell_type": "code", "source": "Q = np.array([[1, 0.0], [0.0, 1]]); Qhalf = sp.linalg.sqrtm(Q)\nR = np.array([[0.25]]); Rhalf = sp.linalg.sqrtm(R)", "outputs": [], "metadata": {"collapsed": true, "trusted": true}}, {"source": "### Simulate the system", "cell_type": "markdown", "metadata": {}}, {"execution_count": 7, "cell_type": "code", "source": "for t in range(1, T):\n eps = Qhalf.dot(np.random.randn(n_lat))\n delta = Rhalf.dot(np.random.randn(n_obs))\n x[t, :] = A.dot(x[t-1,:]) + eps\n z[t,:] = C.dot(x[t,:]) + delta ", "outputs": [], "metadata": {"collapsed": true, "trusted": true}}, {"execution_count": 8, "cell_type": "code", "source": "plot(range(T), z)", "outputs": [{"execution_count": 8, "output_type": "execute_result", "data": {"text/plain": "[<matplotlib.lines.Line2D at 0x7f667cace3d0>]"}, "metadata": {}}, {"output_type": "display_data", "data": {"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAEACAYAAAC9Gb03AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztfXm8LUV19drv3vsew0NwBJWngFOMRhESIEYUZzSKJtEv\nGofEmC8maiTOgtOLQ4wxMQ6JGYwQpxi/YFRMxIgoYlQcgSiDigr4VEAQQRR40/7+6N737LPPrqGn\nc849t9bvd3/nnu7q6jrV1atXr9pVRcyMgoKCgoLFxYZZF6CgoKCgYFgUoi8oKChYcBSiLygoKFhw\nFKIvKCgoWHAUoi8oKChYcBSiLygoKFhwdCZ6ItqPiE4hoguJ6AIiOqqPghUUFBQU9IPlHvJ4M4CP\nMvPjiGgZwN495FlQUFBQ0BOoy4ApItoXwDnMfEh/RSooKCgo6BNdrZuDAfyIiE4moq8S0duJaK8+\nClZQUFBQ0A+6Ev0ygMMAvI2ZDwPwMwAv6VyqgoKCgoLe0NWj3wZgGzN/qf5+CgzRE1GZTKegoKCg\nBZiZ+sinE9Ez8+VE9D0iuiszfxPAQwCc76TrpbBrHUS0lZm3zroc84BSFyOUuhih1MUIfYrkPqJu\n/gTAe4loI4BvA3haD3kWFBQUFPSEzkTPzOcB+JUeylJQUFBQMADKyNjp4sxZF2COcOasCzBHOHPW\nBZgjnDnrAiwiOsXRZ52AiItHX1BQUNAMfXJnUfQFBQUFC45C9AUFBQULjkL0BQUFBQuOQvQFBQUF\nC45C9AUFBQULjkL0BQUFBQuOQvQFBQUFC45C9AUFBQULjkL0BQUFBQuOQvQFBQUFC45C9AUFBQUL\njkL0BQUFBQuOQvQFBQUFC45C9AUFBQULjkL0BQUFBQuOQvQFBQUFC45C9AUFBQULjkL0BQUFBQuO\nQvQFBQUFC45C9AUFBQULjkL0BQUFBQuOQvQFBQUFC45C9AUFBQULjkL0BYOCCIcS4bdmXY6CgvUM\nYuZhT0DEzEyDnqRgbkGEzwM4ihmlDRQUNECf3FkUfUFBQcGCoxeiJ6IlIjqHiD7SR34FC4VhXxkL\nCgqS6EvRHw/gApSbumCBQISPEeHfZ12OgoKu6Ez0RHQggEcC+Geg+LAFC4WHo2rbBQVrGn0o+r8B\n8EIAu3vIq6Bg3lDadcGaRyeiJ6JHAbiSmc9BUfMFi4liRxaseSx3PP6+AI4jokcC2APAzYjoXcz8\nVJ2IiLaqr2cy85kdz1tQMC0URV8wFRDRMQCOGSTvvuLoiegBAF7AzI8220sc/ToGET4H4FfXYhw9\nERjANcy4xazLUrB2QISbAdjJjJ93y2d+4+jLa26PIMJbibBx1uVY5yhtuqApLgNw6qwLoVFGxs4x\nakV5EDMunXVZ2mIBFP1VzLj1rMtSsHZQt5ttzNjSLZ/5VfStQYTziPC2WZejYH2CCPcmwnecXcWj\nL2iDuRI2c0P0AO4F4OhZF6Jg3eJIAAc72wvRtwARHk6ElVmXo6DCPBE9ANw06wLMAkQ4kAibZ12O\ngbBWPO5dAECEFSIcqbYXom+HjwE4btaFKKgwb0R/YyoBEV5NhBdMozBTxPdQjSz2MFevgBpEuAMR\nfmkOyrGJCId0zGZn/fkEAGer7YXo22Np1gWYIVbvWyLct/btZ4Z5I/ocRf8yAC8fuiAzwK1mXYAW\n+DiA/511IQC8BMC3O+YhRL/JbC9E3x7zxi/ThBZoh86sFDXm7UIkFX2NuVW5A2CerQ9LirPCzXvI\nY1f9advWPNf/vGM93acW+rfPXCysVaJfRIQIZZ5vlkV6Nd8V2D7zm3QNI4tfiPDLRHhQ08yJQER4\ndDrlTFCIPoJC9GsLOUS/Vn5XSNHP/CZdw5gQKbVfbaNxPgLgDDcDwhOJ8MpA/lswZwOTFArRR5Ab\ndTPPKrct1gohasxb++kC8egL0fcHr318FlWHt8ZOJ51ga/3nYZ7vGd2OZl7OebtRp6boibAXEV4b\n2PdKInx0WmWpMfPG0AKdrZv6Otyuj8J0RIhsCtG3R4hf7GSKOyJ5xO6Leb42RdFHME3r5jAAJwb2\nPQHAI6ZYlgkQ+W8ttS8Z7Xwkwn5E+P1hSjaGPtrP3wL4fg/5dEWxbvpH7pt3jOhjWCvXZublnDei\nn+aAqVgjnIW6tuck8yl4OoAfJ/J6MoB39FGoBProjD3A21g/rO6dmUcfVp7cjCXqpj/k8kvMuonV\nPwNhUTRHmHkbmjein6ain7fGESJ6e41yJkqa1m8bMurmbwCc23emRNiTKDrVhu0onLkaW8MItUO7\nva11I/nMY/RXsW4siFYvVO4r3LyRdB/IVfQ5v30tEX2orHv0kLeHPwJwVqQchej7w9CKXvLvuoDS\nEChE70Bu6mkS+Lw9LGyDlmvThuinhSHDK4d63Q2VWerbzv/f+iYlwolEU7HQ5hUhfrHXtq2il/yL\nok9gXoh+z/ozl8T6ILt5IkwPIetmnhR9H+1n2v5lamCaJfou5Xs6MJVO8XlFsB0S4atEuHv9tW1n\nbFH0mZgXoh/qNT2GtdoZm3PNFsG6GeoapIjeWjdtSQgAbuhw7CIg1FYJwH1QrTkNdLdustohEfYm\nwkE5aXtAiaN30FTRj1UcEQ4iwh0bnnPmir6OLBEFGbJuFl3Rh9Dk5ujj98pv6ZPoG60ZSoSjiLBP\nh/PNG1LXRUJau3bGJhU9EY4A8BYA302l7QlF0Tto6tHbdBdiPmZRbIprALy1/j+p6IlwJ1QzNaYw\n84fYtECEJ6NZO25q3cTUZgpNF4f+PICXdjjfvCF1XYTo29ZxlqKvxdQX0M/kd21QFH2NporeYg8A\nezU8ZmpkSIQPEuG3ArvvVH/mWDdP6rVgARBhJxEOz0wemgws+3Qdj383gP0bpE8Rvb0npkn0wOQb\nxQRq5X/LFnlngQhMhNv3kFWqM7arol9THv0s4/3nhegXPermsQBOCVxo2ZZj3WzPPF/X37YE4GFE\n+GpG2qFeS5uooD4UU0gdTs26MeWI4fMA3tAi7ya4RQ955Fo3Q3v0oYd4b6gfjrdJJJsZ384L0csT\nuUvUTVPCmUVnrFffsQ4r/QlMj+gB4B6oOsxSGErRD0X0WYpePZQnfh8Rnpq59GObztjcazePIYUW\nqbYtdbtnIF0KuR794ERfQ8/ZpK/jzMNA1yrRe5i5D5YB7/elboZZKHogvz53EeEVU4xmWIUi4z6I\n3r5BxSyUd6J6S0uhzUNwkfpXkoqeCPcC8MhImj4V/RIAEOFDRHhi4pg2oMT/hejrzy6KvinRz+KG\nihF9zoCpXKLvA7lvSLsA/BmAp/R8/pzr6Q4qq1+j922Yr70ZN5ntFrbTNla+JnCPIcJpRHiw3tQi\n72kjpzO2yxKauR69fYg/BsDzQ4mJ8E6izv1hnqKfWV/CWiV6D2vBumni0Xe2boiC83jnILcOZunR\nL5lPjaYRFpbohchD94hL9ER4PhH+0uTZuBxEuIlo7OF5LID7tchvlthAhGcS4U/s9vqTka6jPjpj\nvbfjgyLpn4pqqoymCCn6Yt3UaEr03sXPIiYibCbCfvZc9fS/R2Wevy2aKPo+rJtnZ6b3kEvgoVkf\npwFpN94NFOrgy7VuNpnvFiFr50QAL6z/b030qB4kRwBjFtV3ggcRDolEds0KhCp2/S1mu36Ipjho\nVp2xOW9sofOE/l/3RN+0ArpYN58A8A2MlJOEqd0TVTTDkGjSGevZErkRIHJsF7Wde6zUe9tO1S4P\niDFFT4QPE+FP621NFxKxZJBS9KGF0XVbbvPb9PnkeIlK+1nkuNcBOCX3JERYToRQ9jUIzbv+UrdL\n6MZBXTpjU6IpGeYaOY/9v5F1Q4S7EeEeLc4fxLwQfSvrhgj/QIQ31V+TxESEDQCOBMbCoK6qP6fx\ntG2j6D9GhOfV/zdV9F0sqKaKPoReJjUjwtFEONCkkWsmdXgcgN9NnDfl0d+1/kwp+pDi84i6CTyy\n8B76Nu+m/u8LAGxreEwURHgzEV6nN8FvH0KiS0jXUax9dVH0qbUvuip6b3sux/wngK+3OH8Qa5ro\nATwDwB/X/+cQil5p3p5rGlE7TTx6mYv95gB+vf7fVakRZdZF0efWR0rRZ4EIhyWSnIXRKGKBtBvd\njuUGbVoeyeNIk09T62ZD4P9cxMRALL+mRH8bACDCI4hwz4bHhvAcYPWNCggrek30fVg3uZ2xmmhT\nb8d9WjdNPfpLWpw7is5ET0RbiOhTRHQ+EX2diJ7TIptpRd3cOpIHB/7vE00UvV55SVaUCtXPNvOq\nN01F37WupKzPI/KvgVqv4DpzrNcZm4qWSSl6gZBRlqKvPXJGd0XvHR+aslqjrf35UVThohZtr6sl\nOK8dyTXqat10UfSpt+OQNZeLLh79VekkzdCHot8B4LnMfA8ARwF4FhHdPXGMxRjRE+HTRHhbwzxy\niEkPG58XRR8ieo0JoifCNSaNHsDTB9HLMm2pPLLfGoiwh3QsEmGTydu2Rb3v+PrzWpPGI/qUos8l\n+qXAdnsewSFOWfom+j4VvYZ3Dfvw6HOsm2koevktR6ptKaLPVvSqszxUZ03DK2N9Ma3QmeiZ+XJm\nPrf+/3pUE4zdLn7UBOxNdX8Aj46kb6voYzMDTkPRN+mM1fAU/X6R9JJu3qybG1C93gMBQg5ME7F3\n/XmtSrcZwC/XX/u0buz3XI9e6rpPj96WJebRd+ljGrK9e3lrJd6F6HOVstdp2xvRw3+z6GLd9G6p\n95ohER2Eatj8F8a342aJQ6c1Mlb7qjFFr5XzezL84xROC5wTyFP0NxJhbzSPoz6QCC9ueIxgKOvm\nbvWnXHPrhXt1dGn9qa2brQA+XP+fVPT1IKrbRsobIs6oR0+EVxHhH1W+0yJ6i6aKXue1eq2bTLxV\nr7/7iEi+IUUfJXoi3DbjTVLnk+Ixb/80iT75QCLC44lWo6u62kYT6I3oiWgzqvCu42tlr7D1WiLa\nWv8d4xzuEX2swXmNIIeY9M2Qa908CcD/ycg7htirXY7NshvAizCKz64OpLGFwkN1lzOtsYemir4p\n5KGrQ+0A/8Enc8boV1p9I+p2HPPoD8DIkvpFs88j+l1IK/onAvhD+ETfdWTstKwbT+TkEP6TUHn8\nIbRV9Dq6SvfVbCDCW1W7z3kAhvan2m0bog/xS4518xrgiU8koq3AM+4DvKLLrKnBAnYCEa0A+ACA\n9zDzhyZTbAUzb63/znSyaEr0eh+bzxiWUd28Zzj5x6ybrHoiSo7GbKvodwN4hbP9sozztLVvmir6\no+u54XNhiT5GZt4+Xb5cj54wKq8lKHveJVT9Tymil7eNGKE1QW54pUWXsShHqbfWXJUMpNuIrm8N\nfY4mJHxLVIMAf8HJJ1WOUBlCCI183qkdinpBEwnJbaToiXBrNY5nCXjfecy8FfjHbwOv6nXwWx9R\nNwTgHQAuYOY3pdIHsIwqdLCLdbN/HSefOs8n6/9X09bHxYg2eRMRYU+MvPSJ3eZTI4focx5inoLI\nPXaUCeEZDY+TdA9ENTd8LizR2zryVKa+Dppk9O+N2YChB4WXPkX0Uv4fmu/AaDKzvqwbL2rEpus6\nj4osYu6dKwRv0jZdrj0wbgvJ/EMpRR8SXfba5r59uG2BCCcT4VHOvtgDbAnj43C+AOC/1D4PIY/+\nbIwWTFrG6F7YhHScfyP0oeh/DcCTATyQiM6p/45tmIcQvUauotdIzSi4gmqO8I0Y/+2rF6D2KC3J\n5ail2I3WlehzoPMOEWIOms7T4uZfj+57SOQ4mVTKEr1XH00UvT3Gbgu9ATZV9HJOabenO2XrGnVz\nr7pvZmiPXn9vouhT9sJ+GK/nn5i8m3bG2nETXTz6JQC/B+BpiWM92LrWUUSCHI9+D4wCV5ZVPvNH\n9Mz8P8y8gZkPZeb71H8fa5jNEiYVfZub5NZ2AxGeQITfr78uo/J7VzCpAmNKpuuo2VjeOb8z5zqF\niL7pA6SpGg09SE7GOPlZyDw8dgRqrI5CDzDv+qSsm5Si34CqTW6oO3LtkHS52T2ySxI9EW6WWIgG\nqMKVXw6/btqIkRia2ESCHKJv3BmLcJttS/QxgZV6K/Fg69ojeu9c9gHxPZPn/BJ9T1hGpZ66Er03\nNe1JGL2WLqNS9E2Jvms9dVX0Tc/fRdF7C2007RgHRqsrpR40tvM0VkchRZ/7AI0pes+60eLjbma/\n3LQeUeQ8LK+Frya9MM8hFH0I7n1AhJOIcKJJ6xG9LuPN0a4zNmTd2A77Lp2xMaJPwRL6LQPb7bns\nfj0uZO6tm0ao1ZENt2xK9KF9XhinvgFi1k0s1jVbLSVUmrdP8u5K9HdRUyG09ujhq9HY+UP55w76\nCHXGtvXo7TEa+k3TPgBT1k3IE2+k6IlwNtHqwB1v6gp7zA2BMngPJjnHrYlwfyfv2HlS1s3TgNU3\nY0GKJPtQ9H1YN7H7rg3Rhx6qY9uJ8EZzfmsR38Ecu3CK3i7k7RF9LnRD8Ihek0Nb62YPovF5Vojw\nz0Rj4wVyOobaKvoctXYygE8552lK9JLe1lsqvUVovVRbByHrJujRU7U2p441zrVuboZmij5G9DGi\n8AZPCY4E8OZAGbzz3Ig8QtPX6PUAPj2RMWGFCLejai3g0PJ9sXPZpRFT1s1mpIm+yT1vib5TZ2z9\n2SYqLST87HYJy5ZzfVKtxPZYjKKH5FgRPXuguu69YapET7T6Q+yEQk0VfYhcvJGvltBD1o1u4Db/\n2wF4thq+fysAT0c9X7g5T0xdxtR+jJBzp0yVB+g8E71FzoApq4iuAFZn9NTHjh1Tx13f0mzP9eil\nTUre/06EQ9X+GNGnrBtR9DESFGiiz7VuQuR1AoDvoxrUaOvFnt87140AQITDifADpK2bJYzGLezv\nnKNpZ2zIumnTGdvFujmSCB90tofeznWd3LH+tIuvL5SiFyL2birbGRsrW6jBh560cjFXUKkSz7oJ\n9YzrbXtF0uSQeVv/P3fwhleGpoqlE9ET4S71v7nWTeh1PBV1oxG6Hn+AyQmicqNuLNEDwN3NfqBl\nZ6xJpxGzbjZE0mmRE1LasalJbNvZUC/GcygR/r7eJirzKGA0ylhNOmexjNFvvEht72rdhAg/hL6t\nm8fBj/CzRO/d9zudbcDAHv201zCUibdCEQ5trBubjwep3GVUQ5/JpE1ZNzLXys1REVjsBm2q6HOs\nmy5E3zZscznwv4Wti2/W5xdFvwV5yKm/3Fdmyee2zvamil7bjCtmv5eP3pZqzzkRHzci3n4EerBe\niOhjE67Z+l8C8CyMTw1t7QRJuwL/t8j9BozPzZSybpp69E2tG/0AzyZ61f+WaoeW6PX5d5i0+tiF\nUfQhoidUN0fXztjQ9t1EOBxxoteNT5SKhNTJzS6NNTXaL1SuPoj+25F0Nk+g/cjYrtaNDB47xGy3\n3wVNo268/TY/r2y5Hr28Zeq8NdHHZreMefReOo2YdZN6y/1p/X9orvWch7ec448wOf+/ePSW8HT7\n1PWxP8Y7HAUPV+da/U1Eq28lE0RPhDsg3J+TmubDXqPt6tgmM3fKMaG6S0XdAKNro3+31MPCRN2E\nrJsNqJ6suUQfQuiYPQF8GeOv4yHrRm+XVV6E6CWSQc+/8S9EuB/iRJXt0RPhfCK82qTRN9JFCKMP\nRd+XRx+6Fg8LbG8adaMRsm68mzik6L1XafuW6Sn6mCLtw7rZrcqWyk+mrh5T9HWk214Yr6dQv4+c\n40hnX0jR27n5U28gstaCtW6+W396bepSYHWRcVsfvxQ4jy2n4Ca0s27kfKG6C3n0KetGjlupSX8F\n+avJZWHaRC8WSEjRt4E762QAUoExRR+zbv7OOc/voprcKVd1WViS+kVgYli2N7w+ln8jj56q9UMf\najZrctdhYbeq1VUq/6Ztyz5ovevah6IXNO2MBcYJLUb0qfKGyuDlp9uqK4QUscrsnp51sxFxRW/L\n7FluNurGI3pC+j7Ux+v2f8va7z8OWP1d+hreUh2nP1MgAD9Q329Sx3r3U6jdpBS9bI9ZN55okeNW\nUNXlduZ+p46eNtGHXnc3YFjrRrCM0UUOefSeOrThoM813xn+RbXl8upbzqdv+r1NGn0jNSX6nAbz\nCAAfD5QLGG/YHwdwacZ0tk3blpxvMxEuwWhBbJ1XzquxIEX0XnrBjahuOkv0nqJvEzQgyLFuNNHn\nPuhCHcQ5ij729iB2gp2p0+aVS/SEyWt6OIDXqv2xTnltf6T4Qitkbd00GRk7RvS1HaxhFf0BdaSW\nvj7etdQWWO+2jT3ZNBAiw1WPXs2P0qZsOUTfxLoRWKK3IwT1cU09em+EpR1Io2+knDefpkTvpQ8R\nvXT6peyE1PX7BoAvOekfgioETc+9n1LITaybnPIyfEWf69HndJ4is3y6rep9tyTCfUxZJN0q0SsC\nlN8ksB38OW8h1vvvSvRW0QM+sdvv3r2WeuDqsmdZN0S4U+D88nvtDJPe2+h/Ybw+lswnMK7oF4Lo\nQ+QgN+UfAmg7A6aXr0WOdbPaGWv2Vycg/LmTr1b0TYleGo0+p32w6JsyRvRtFb0X5RB6zbcPgrZE\nvxujmR91PtJhfysnr1YevVF7Y+UlwllE+BdM1tmsiP6OauSs5PHK+n9dlocC+KopywTRm3LG+l1y\nymzHB3jXRaybHO/bI3pdJ3qAmz6PxyMTbyhEuD8RTsCoD1CgrZvY/fQO891aN/Zh6RE9w38gWf6R\n/BaC6EONSV+IewTSjKXPmBDKg0TdxKybVJ38amB7287YnKXv2lg3Z9f/57wBeA+DECnY1/Ymvrnd\nr2PtbYz0JgAgwmMxeoVvoug3YNIrtf8DwNEAft3kvRtpos9ZFS1J9ER4IhFer7bdA6NrJ3k8Wv3v\nwT50NNFri8lOB+Ihdt121Ss/3cGk9eym3DdJa1Pq464G8GCnbDGy1Hg5gD/H5Ntdyrqx50M9Cv7l\n9Vc7xbYgJET0dfMejqudsVgQopfzPUq9dgL+a7bubHoZEf6YaGzIsO081PmHsALfo9cDqDRBeAgN\nTW5r3dg0Hpoqek2i2YrePDxDRL/b7G+r6JcwTvRy88jQfPnNH6zz2onwjRSqcy+e3bsWy+a7KPpY\neGVOZ2yK6BnAi1GtHhZCjj3RRtG3Ivr60xL9Uj0thSCX6DdgciLC3PatPwH/98g2W57cqBud/xEA\nHmm2d1H0+lr8g8pvI3qOuAGmOGCqDhuSC/X3qDxavVJMrMJtuCEAPKX+bBJ1sxGjASi68vfAKKIg\nNf9GiOjdm5sID8boLaAt0beJupF0TRS9Tptr3bRV9EsAZMnJszAaEm5VE1D9Hku6qXPpm3uD2W4R\nI/qQPZD6/al9OfsB/yFlkUv0OYo+x7qx9boEQK9BMRTRxxR9E6LPVfS/RoT/y4y319/tG2Jf1o2M\nsl2BvzZHZ0xT0dubaR8ibFYx6KlwNwuZwEyn2w8AiPBmotWnpMZGjBS9bhh7YFIthOARvb6Y9vhn\nq//tvp2RfRpBRR+xsOxNCSJsIlodrJKCVhz3cs7fujO2LrMm+u2Ij/5NKfpUZ2yKLG089274UTdN\nwysJCF4jOW+KENsQvfRNbDD7Yoo+p/1Le32CSbsB4+2yiXXThujfVltIE9YN0dhkd5rodb6fRx7R\nA6P+ESBM9J+tP712qPvvgMp2XgJWl27UWKnzaDMtQxTTJnp9vs2owqg+A1/Rp8rmvY08sJ6X49nA\n6pJ4wGg4vlb0luhjnbEaNpZYELpRQqMGgXGiT1k3IYUesiUmiB7VbHreojApj/7t6v8+FL28UYl1\nkyL64xFfqDtF9G0VvSX6vcwxkt+pgXLtEdguyCF6ff5NgYeGtdHkcwVhRS+Ta51mjoldt9BAoCW0\nI3pP0ceO8/r2BCv1yFotxFZUOp3vF5E/YOr2RKt9JFY4SH1KUEGOor8bqvv+Mc65NmIBiP7XUC0g\nLtgHWJ0N0FP0KYRGof0RJhuEEIpM/zmEog/5svo855t9uUS/otJKI/gcqpGQIdvCI/qQVZciei9t\nF49e1KUm+juGk2MLxknLwvtdmmxSHffWrgsRvZ4GWxO9HYMg8GYT1Wiq6LcCeEEgH2DyrZIQVvSC\nbeZ7EyvqJWp7W0VvZ5xNrdeqoR8SelIwQci6+aHKK4d3JNTXEr3cY+Kpe5xkFb3Xt6jLu6aJ/loA\nj8ekshJfdguaWzehm371WCIcXH+XGNo9EbZuunTGPgJYjZyIKXqLJopefoPU0/sx2YhSHn3oek+b\n6IWAhOjtYDkPMUUv2220Sq51YzvnGcCfohqqr7drUtK/P9R2Y0R/OdJtzTv2nmM7K6vCEr18f3x9\nHsnHeyDuVPtDZQ2VRVbdaqvo/xCYiFWPHWevv57IbQUjy0q/0Xjl+SFGnf45RC9EbtuT1J3MMeQp\netl+AoAzEO4bAUaKvu0sAUFMqzPWzmMjkG33w6SlkLrxc+Zo/079KQ1/L1RETxgnYK3oYw8QwJ8w\n6k4YNVjbGGPl1E/uFNFLo5JGIGozRPSSPjYAReDZOSn1b+eRt8hR9OLR5yiYXUhfmzEvFPHOWNkn\n+epjpY4faM6pFaQm1xA5rRI9EbajEgSC89DcugEmr8uNAA7G+INQPt+l0ml1r2E7/lIPaA+WnH4V\n+aOS7bz4sfvAll8T/TLGr/cuhD36H6nzrP5eIhwP/361I4Kton8TqplSjybCd+Er+t0Y9f2EIOtk\nrFlFn6PYmir60Ao5HnQssafoNyHcmWpxcGJ/zLqxaNIZa4k7RvS6saQW0Q4hRfSbEF8VLEaA1qPP\niTLQk3uFYOsiFkd/HyKcj0lFC3XcLTB+/casm1o5Hlyn17HPnqJfAfArKo10LjdV9N41XMK4zeTV\n0wZgYsi+lMOmyy2L4AEY/x23QT7R554DmLznraKXsuuRpsB4W/wrhN+k/xLAG5zzyrW1il7OtwPA\nFwAcA+AgTL7F713nkSL6Ne/RhxS9TaORSm9Hj8agK1cUfcij/2Iir8cn9g9h3eiQqxjR70OE30TY\no29yvVcAvMrZLo19E8YXxViFmnY1NF2utW4s2XgkwUg/qKxyT3XG/qI6t7VugCqEVz/w9OCeZVSD\nrR5apx+y9jgCAAAgAElEQVQbRVnXwSaE235bovcewEsYDwUN9dt4kN9/l3qyulgbCcX7vwbtrBsP\nsfPbe17Pca+J3g5osg99W1ZBaLGcENHr6Qw0f+lOeEbVzi7EaCBeCGveo88heqvoV4jw5Ej63OX1\ngPHKlRvsrmqbJvquaKvoU0RviVviyq1a/QAyPfp6BaE7B84pi0nYG1a+y9iDUBTIGNET4TWmHJro\nbcP2Rgb2regFl9Sfeqk7qTPP0tL5rM6fzoxnAWPLy0nHf+jmFqsl5cf2pehD7Us/jD8USZfCoeb7\nUIpeDyZaUefRUX1W0ceI3obVekgp+qsjxzKqt76LkWfdrGmi96wbWzHej/uLSJ52StAY9M0h6R+k\ntunO2Nw8Q+jbo38Jxok+5tHrvHIU/UMBfCtwXnmLsNcpqejhED2Al5ryaY/eKvpQZFNbRa89VVtf\nV9efuqNV6sw+cGwnnB1opueC2Quj9Yl/rLYLchV9yqMH8oj+a+p/rVxDE5U1xSvU/4xhiH4vjDo+\ngfEHZcy6GXvDMFMA22viwdaRJvrfYcbVCJOzjuDSRP9DJ+1CWjd21GGTmQaBZh3Jq2kD8zynRsM2\nQRPrJvQKqXERxtfejFk3Oi+P6C1RxuwvCemU6RFeZvLbhPFl7jQ8orfli3n0qUFpIVjlbq2bHzjl\nlRvPG/0bGgUdspr0DbqnOn4TJjFG9ER4b+BcNiqlraLXHcn6uujfshvAVwLlaIK/xWR8fC5S97we\nx6IflFrRr5hP/UCw+eu6ChGspLHzUulO6JiiF/LWb3eXOGnXvKIPEb1NYxErn1xEOymShxwl2Fdd\nNLFubJSIB5lH3BK9WDepOHqt3mIzBVqIdSNpZBoKTfQ3OHkCVYx/jOjFupGBbLYc3lwfmuhzFjsh\njOp32WzbiZFaXzafwKRyX91OhHdhciCUJXpCeii7RPvIsccF0tm1D9p69Bo7A//3Fda3XzpJEKmy\n67ahBZAOqFhW2yTqxYvAAoCtRHh+/X/o90sb/1H9mQp6sJB0OhLI47uFUPQW9iYYUtGnoEkB6Gbd\n2HCxmKIfU6B1B57gvaiiPMT2sgo9V9Ffp7bb6x1qULuQZ91cB1+t3gujKQQ8WOvG/obYcUB40if7\n4LQ3vtTXbox+u0f0K0Dw7e8p6n8bN63r07th9TlsuGjuAvBtFb1GSNF3afd9IfVmrcuuQyr3weT1\nBiqLJNRfI7hL/Rm6HzwrCBh/U4hZNyLUtHXj1fWeWONE73n0OUSfo+j7QIgw2+Cz5nuuohfyE+xi\nxk8xuug5nbECHT66w2zXCKkQUR6r1o2CVvTXIBzmuj/ChL0EAMy4kXniIQv4Slgr+tA0rrmKXr+p\nSD3pug+JiFAEjCX6fVHZYvaG1VFMcv3suIQUvLLJgLpcoh9a0XdB6j7U+zXR74tJ6wYY2XUxopdB\nZSlFb4m+i3XjEf2tMa9ET0THEtFFRPQtInpxIJln3dgf04Tov4h2RB96lRbPuBcYZR7L1xLTAeq7\nbUAe0YfyJ4wW8ND7l+ryvcfkaSHKYyzqpp6MSb5vQdUxFrLFbosw0evOZSmvhnccA9hAhP0RrlP7\n4LQ3fkjRyxuMLp+HUMeop+gfivgNm9sZa+GV7X/gD5gKIebRzxoporcTs0mZb6aOXVaTm12LyUgb\niyvqzxA/yBKknqJPWTes0mlF76XfDF8gdEYnciOiJVQdL8eiihV9IhHd3Uma49E3sW4+hXwFpBF6\n5e9T0QPjD6HcASgE4DL1PRTRoUmF4Xeo6tA/rx/gSSZvC23daCLaUx3zSFSNPzTJ220QV/Sxjugg\n0aNSX3Z+FEFI0UsYnqfofwHV72i6KMcOdYz3+n4t4uSyE8DvYHwlrRzokEK7vaui751gWiBVdvvm\n5Sn6ZQCPq/+XB35M0Ut7C/3+JVRtxI4G1+04Vne51s2VqATSfBE9qsn4L2bmS5h5B4B/gz8rW07U\nTZPO2NRsh/bcgtgrf2hQTRtsJMLe9QCUWF6x+G5pQHaeDf0A2A3g0kC+YiFtIMIt6reMV5l0MUXv\nefRabVyFONHL3NqhfTFFn+qMDSGm6KUNirrTv8sSfeg8evt2hBX9+RifcdTDTlQPQ08YxRCKlGoy\nMV9Ixa8Fj94SvZTfWjfy1irXPEb0S0S4MnLuropeWzcxRX8Nqgf/3BH97QF8T33fhsmFrYH2Hn2o\n4pssteXF0HvnSUXANHmtXQHwVlQknEv0If/8RvNdk4pE5FhsRDVA4zcB/AaqWPGX6QT1q20s0kfO\nYRW9+JkyICikSJdR2Tuhfa2sm0B+Ar3/HQCeac6nrRtdb1eq7290yuPlrxW9JfqdmHyYWbS9me8R\nKJ9dJS0GXb+6HGvNutGKfjP8zli55qGoG6CyTG+Nkd3pnXM72nfGelE3HhfdiIGsm66RK5kK4LkH\nANdvrJ4Bx9R/Ez/G+3F9EL1GqOPQWjchos99MG7EyG+PNdzYW4SN5fasm9CNuQcmB60cZNLcAmF/\nXROWPodEBcj/seuwjHDoqyVBW6+pztgQdB3qh4woP9lvRyheW6d5DjPeSoTnBPLX9eURvb7pU0Qf\nC738FkaRILloMl9TSNHPA9E3tW52q/89otdvcfLdQsYqhCzBkKJv0hmbYd2cugfwofsCANHJWwP5\ntUJXov8+xm+oLZic3xrA31yGyn/WU6xald9E0bddUzFG9LkdqDlYwYjkJn4DEf4IwJfNvpB1EyJ6\naz9obErsB6pY59D1l+NiRL8H4kS/gvircFOPHhi/LndDtRxlaL8ti7Vu5HecjapONyKu+uz2rore\nEv1PMIo/b0O480L0N8EPuc1FW0Uvg/SAyWUotUfv4cD6M0T0oug3qe9y/pw4ei/qxkl/3BXAcT8G\ncC3zSVuJ6JWRPBuhK6F9GcBdiOggItoI4Lfhr7ZDwNjiwUDeFAih8oVGLLaFtW68xtbkJtgTwP0j\nef09RqvTh87Zt6K3+cugJQ05x63Ud53HbTG6JjmK/vOBfTaErI1100QYWOtGR9m8FKMbMKb6gDDR\n29f3nZi0pyws0V+j/veOe2MkL6B6UGmlGUOoA7YPj/5/Ox7fB9HHrBsv/02RfZKfVvSa6HM6Y8Gc\n1Rl7E+Yx6oaZd6Jatu+/AVwA4P3MfKGT9EjUfo1CaCCORt/WTagRWkUfmj0xF7rDLPQbNqGboo9F\ndciIQL3/kSaNR/SWYK2iPxyT1s3L4GMJ/pSvQLozNjZgKvaqHKrrWyOs6L+GEfHHyADI74zNsW70\nvmchHf2Smsp5BYCeOTSGIcMru5JUiuj1b9NEGyL6nM7Y1BuIRN1Y60aHd8ZGa+s0KaLfG/NG9ADA\nzKcx892Y+c7M/LoGh+YQ/YZ6EWCLHKK3c4f8FOPzgWvkKPomRK8tohjR6/o/WP2/G3FFvy8zvov4\nAA+r6O0blUf01hKz4ZU3x6R1cznCCJXPhpA1UfS2c1oj1J7/HZOKfgnAfzDjRxgRf4rodP56sfKQ\ndZMKrxRcgfHf7B2Xc/OLbdKE6PW5+iCY3IfF8wPbU2UPzV6piV5bMJbovfxTEXxW0evtKesmFJlj\n0zOqtn0Mmq0bkYXeBgm1QI5108WjP91838kcjaNP1UUTtaMVfShfq+h1/4VW0pbYdjOvTmsQa1xW\n0XtpUkRvw2KlnnZjpOhj5wjtOxXtid4uqagRU4PSJ+S9FYh102Qmyd2oFrGQsgHjRJ8TXinYhbSy\nTpHwNlShrnsi3Zb1ufSDOvfBFEOuIPp6YHtK0f8Gxucp0opejj1EpddtNpR/avCljbrRx6WsmxDR\nxwTs45x9nTBPRN9EoeUoektasd9qo25yPfoLAvkFFT3R6qAoS/Qa2n+3pKcbSIroYzcdIW3d2DqT\nY3ZipOjbEL3d1xfRx66xtm68/g5L9F7dhfoIunbGWqJvat1cxYwtqESBnXI7du5fx/ikd7H2kmuX\n5gqiULoU0W9DNU4BCFs3esZPq+g9/ELinCFFr69x7Pd4DwNvahERdaEwz9aYJdEfivEK1hUlizw3\n9eh15eUS/XXIi6P38InA9ph1I1FKmui9mTwZGJtYy/P2mlg3FjkevX2lFXW0E90UPdDOo9+A0XX1\nflvsldfG0QOTij6luHS78OLPm3j0ltibWjcvdNI3VfQ2MitU70D6Lfqb9WeO/bMF7Yle5x/qjI1Z\nN6n8PTyzPpedtTTUgX+y+l8HHhyotnvWjfBal6glF7Mkegv9w19Sf8ZGxnrQlb3LfIbyuhqT1o0X\ndtjkdTpnPVvtpdolzLyIGm8ypBjR2/BVC4/o7QPUNjg5ZidG89H3QfRN4uitotf1EWvP2rqxKixX\n0YcG33lx9Kmom13m/5R1Y+vkfU76XKLXQ/77Ino5Z9K6Yca2SLpXB7YLrJ0o5ddx9NY67Ur0QNWh\nb7FRnV/fK1Y0yrX+TbU9NGBqEMwT0euGL8TXR9RNKmTuGkxaN0KqqSlcQ6/TuUQv5/SIPnQj6O2x\nKR1SA4w8orfl8BS9EL2cvw+it4hNgWCJPneBdRt1I3lKXnaCMa/+U+M+bHRFE0Wf6hTV275u0seI\n/nNOXroOc62WXKL/dCLdVercHu6bON7WjafoLenqa57ivC8Ftj/V2ab7YUJEr8/tQY7/EtpHEyYx\nT0SvKyNF9PLkuyqwXyM2iOT3ALwZo0EVAiF63bi9Gz90I+csXK4HFP3Y7PPmgZe0uhw/hQ8h+piV\nkUP0Xuz9Mkb1Mg8efezNQKMP60Zv86wn/abQlOhTUxHo9Mcin+ivjOS125w3NoAyRPQfrz83APgP\njJZmdMG8qozbxuzvwng79Tx6bbHYAVNyrR4Lf8yPV/c3Afi2s30PlV4Tvf5fK3qPr8TyegBGbftI\nJ10nzBPR6waXeoWRlYlCg4z0Pvn0fut/oiLLkKLXT9ic12lBzoRrMpvi/VENNNPwFL1H9NfDh6jT\npkR/rfn+b84xK+iH6D2C/u/6M0b0dpI3fQ1SHn0T60aX3Qvp9B5Uluhzo1gs0acUvaw/astyI8Yn\n90qd2yr6WPRJyNZ5NKr+ArELe48BN9iF8aiUFNFb6wYAwIwPA/iuk38omsvbrmdz3WS26/NLnRyh\n86zXYvjfujzbMboulzjn6oR5InpdkSH1IBUmN54XFWIRI3p5ZbadsUIYKUXfhejlnOdicgZKz6Nv\nQvQ51s1nMB7SCVRTTq+C2Y3CWcHoAdh3Z6ysZyD1auecIXSzbqyil7p8MKq3sJB1I6o4pfj1gyMU\nXvlGlUaXLfQQ0WkEIaK/AZVKDdkIAqkze96Yog8RuIy1kDY3ONEz4woAMjAz17rxPPoQJ1jo8EyN\nkHWjHzR6fIZ+E5FyeIPX5m/AVI/IIXqpCKkca3l4sMpeQxq67YwVZZPy6EMXRHfSxCA3h53q1yN6\nrxxNrBvvNfXwsYwZpyXOR2hO9IcjXXf2IebND57j0ccUfawz9mbmuy4LMForVOefo+i99iFtO9b5\nmnp7DBG9pNFCw2v3IY8+RvTuda6H9muiz423bwsbZOEper2Aj32490n0+vxXqG17mv8lT8+D94i+\n9zqcJ6LXFZwiekl7BoDXB9La6Q5CF1UaqW4AHtE3tW5yOlYIANc3i803pzO2iXVznZMupy/B3hj6\ntyWJnhlfDZzbI0rZ1pboU4oeqObktx69Pof3/3ec/D2bUJNPKOpGjrPWjRcxhsC2ENFL3qk3yjbW\nTew6t1H0baNfLMFrRS95aqIPefShMoSIPnQ/Svr3orLNgHGi36zSaF4rih7hH2oJYAkjtWXxHfM9\npuj3xrhtITeKvgjeSL6Yok8RvXj0IbXbxbqR10U7ra5FTnSQRlvrJmVF2N+23Uljid4jzFh79tS4\nN2gFqOpVFNpBgPum4z2oNIGGFH1borerQsWIPjXSM6To5bjLAfyaOSZF9E09+r6I/skAXoCRovfu\nm9AUCF57yXmTnkjPDFYj1vV9pe8Rj+j1tnWh6HOI3j7xdKgfoC5GPdAo1bMvhHpns91bWPvpzvGh\nRvGnGE05G4IN+7Ll6ju80qvT0HzxIbTtjG1K9FLvP1dpGNXQdqlXIe5c60aTn5THLpoi2++KOvKB\nGZfCvxYpRR8ieu9NRPoPbBqYNKjLlPL0U9aNvo+8qJufM0+EZbZV9H2rU0v030AV+SNx9N5AyZCg\nylX0rAYuWnj1IlbsW1F1+O6qMnHzXihFn7OIQhuiX8Z4RYWicEKNVAjVEqZHIqHY7i4IvRLmKvrQ\n7/KsG69Ob5dRRg2r6HMHTOUSvSVBvbaB/O7DUSls6Z/Q1z/WnjXRy7kPMGlkNPIP68nOBF75PaLX\nvzel6G1nbK6i/3enTH1YN29BfNBTLtHnQv/eP8g9SBGurutdGCl6++ZqFb0OfMj16OVch2K0alko\n/UGoppYAqs7XkIW3eNYNMy7OSJYKLwMmrZvU6MMo0deNRm7K2Pl0Xho5Ha4x2DeSCwE8A/mKPvTb\nvc7Y3MbzOPhhZ5JvG0Uv5/6Qsw8IK3q9TKWe8vivVb+G98DwoDsavSgIW6bU9tjDS9qUTfMF5Fk3\nsbxf66TJtW4kD8+6uRDpzuwYQh2WORaNnWk2B02I/s6oFknaD8DLE2ULPtSZcR4mF1YaS8+MS5VI\nYEyOz7CYWNbR6bPrjHm1bkI/1D7xNiLuZ8UUpYARVkCpqJsuF0Q6iXQelzDjn5A/YCp0fiH6y9S2\nENGPvc0w4wMIj2MQotdeYkzF2fp/WqA8G0w6yf+HKs0ZTr46bQr6ppZO6FhnbGp7zLpxwyuZcRQm\nf6OkTxE9m89cRe/l64VX7kZc0XvhpnqftLkYsT8rsL3Nm7El+o2oyN4LCT4MwNnMuJY52aeTus/t\ntQndg6egmvSwiaIfDPNE9NqDDF14S/T7In6j25vDw260J/qcePkQllFNndzUusmxuJZQ+YqfUOlD\naf/W2Ra6WQlqyb2MfhBL9KGyiw2zSmL1YJLQOAb9f2i6Wwv9m2Ruftv+Q7/lGmfbtDtjLZm3sW4k\nnafoU0Tf1rohky61PReW6O8I4A8xabESqhj3n2MSXjvXZZG+u8ZEz4zHo1oesgnRd7WCg5gnoo81\npC/Wn7YjZl9UKwSFkPLoZV9otrjUjIJfdLblYi+E45ZzrZuYR29v/hDRe+o9RPR2KTVbHotcon8m\nqjVg7fUK/VbZvjcqqysH+jftX3/a9h+qz9NQEYlGbGzAbqTDK2MevZTj42rbmGgxHXu2LCnrxvPo\nd6OZdaPfziaIvn5QA3lE3+bN2BI9UM2G61k3K852IK3oLd/Y/4H4HECedXOWSVOIvsaNAH5Xp1Eq\neN86TvsJgWNziF6sm7MxOZ9F9CIw4yxU83v0jdwBUyHy1p1PS4m0uUpnF6oH4g6ME0lbol/9v36l\n/qZNZ952JoieGT9H/ox/+jfJyMUsRV+Hzl1mNnuKPqcz1nvDChG9F4njldGudJVS9N5cN+Jzh2Dr\n344pyFmhK4Q+FD1Q/S59z/4BRnajR/QpRa/ryZ5XELNeJtoBMx5gzr1w1s1jE/tjpBVq5KkQxvNR\nzRuRUvQ2ekeQsm6kfDm4LcY7F2N5eYre8+hzOqJTaT2i1O1Cl0mIXr+J9KHobV6pzs+uHv1H689c\n68ZDStGnrBt7HT3rxhug5Z1Xrsc/1p8pj96bpliHeOpzyKRb9prYsuVE3Qxl3eg+I30ffxLNFb2G\nXAeP/AUpRb8xcG4vjn7tK/p6EqEYQmSsh6zringJgK0qjYcHAPglAP+3/ouddwcmK/qcjPLl1iEj\nn5S8zlidT6hMD0MV2eGlS0Uyaej6lM7a7RgRfe5rd1ui12X/rHOstgjaqEiZvC3XuvHgKfpY1I1d\nMMU+uLw6XQJWbZCYol+u052JaqBgSNFb60aizuScuZ2xQHeiz9mezKt+6wspeqnXENF/J1EWr93a\nfGJiS2zh2D22cIo+BVthp6BSwG4jZ8brmfG2WIbM+BkzrmfG+5nxz6Fk9aclsIMxPntjG0WvZ9lj\n5F/UVh49M06H/1C0A2M0PHLTv+l+9acm+tBUAKG8Pa8zZmt46Vyibwkhwr4UfU5nrF3nN0fRe30h\nMUUvx6asG91G2nbG2jZA6jOEoDUWOSYELxzbEr08gEJE/ypMBiN4nf7BN8hE2cWj9wTeuu2MtTf+\ny1DZCnqaz74q4mj1v/Ut9fYc6yZUhy9ENUOkPj5X0ecOmPIQIpIQ0UcjJZhXp0zdjsrb3o781+4s\nj97Jq2mHWA50ntLHINfuXU6aJvl5it6GVwrR51g3nkefS/TykPFgr4Ml+phHbwnw3RitmBSy3GL1\neUVkXw605Sh1J4t4C6KKvu7Qfr7d7PyvtzVR4NYq05AyfwmjN8x1QfS2MuQihawbjabzZnhEtQPj\nE6HtQjXIQpcnlZdeWOArmCT2aRB9iEiaKHqvXWxHFZZ4FeJE/3aMOqjtQyfVv6CvhU0XU/TnoApl\n8+BNXWGJXqZH7mrd6Dl6bNRNG0Wf2xmrif5u9XF/Ycqm800peg233dVvyx80+2zZ7lH/TYAZXwPw\nV96+TKwSPTN+DOBjqCLZmih6mQd+C4DnebvrT12HfRH9hvr85zMn+xo7Y56I3pKgXCRt3YRuxB80\nPNfZAJ5i8tyBKnrn2bKdGdcCuI8qj4fQfDLWF8+xbmKdsTqfGKy6Oh7VvO6hcLyoolfYjmrdzCsQ\nJ/ptGC3HNkZO5jU3Zt3YGHO9b+KczDgMk1M92+Njij5GoiF4aT+Dql/Is27sWAFbF147yrVuPCV+\ncv3pvSl4b0ehzlhL9C8H8ApzLrdszLiQGRc46QRd/GkbRPA5VJOJSZ6HYmQlBYm+Luc2jCYI9Or3\n5up/LxInmHX96Z3bu2brQtF7RJ+l6JnxaeRNuSvpdzDjPfXXVXKpn+6ivCzB5LxN2GmNrfUj3x+U\nKGJuZ2zMSxVyfQszTsI46VySyCNE9ECa6HUZYjdDrqL36j/nhji+/vTeovog+glFz4zddcit9cKB\nSevGeszHA3iNybuNotd5hsrs+dta0evlAMeInhmvYcbbTb6a6M93zmvThb43gSX6mzBO9Bcg3Rmr\nkWsP22myY5A0nmJft0Qfsm7kfyBSEcy4IbY/Aksu9okd8h8FIUVvPXn9/fJEmWLhlSEL5o7qPPpz\nrGx1BMd59baHqP2poeESfXM10p2xngq3iHn0KUUfOxYAvopq5S6bl8B2xuY8mGLnsw9Gj1CvNfu+\ngZFltKueLfJkk6aNRy+Qc0etG/WWxRiRz2PqzzdhtCJWqtMRqN6Cz1WDpULpmsLOLwP4RL8PRm92\neuWrtkTvvSW2UfT7O/u8e+ybzrZeME9E7yl6aXi5N6I0sKManNf6aHaahdgNps+pj5F8Q9ZNLPZW\n8gl59G4dqAE9UaLXeTHjDJXubs65NLarz5iiT4UCeuUReJ5mbtRNyBYYyrrxPHq7T6exin43Kl/Z\nK6Nn3fSp6EMPyiUAP6g9bzDjuanINiePJshN7z2sLdHfiCoq7EZg1SZkVHbjvkgTvY1gI+aJSe+A\nyXYVg4y49x4YE4o+8ZDshHkiensh5GksFwzIbBjMq7HkObC+cEjRN/XouZ5ASYaKa0WfQ/Sh86Ua\nWDbR63TMYzHFMaLPCa/0ojssYtZNI4/eYBuAz2Pyeuo8ZqHovagbq7BtOZ4H4KnOcRZelI1HjvZ8\ndlR0LBY+R9Hn9h81hfdb7FoK8sapAyKkHvdAvqLPLUvOb7ll/enZNLEIp97RieiJ6A1EdCERnUdE\n/0FE+6aPCsKSn4R7MdoprlyErJvQzWeh57vxSOVjalsXRa/3xRCqqxDRe/D271CfuYq+qUefq+iD\n1g0ztjDj2XAUPfPqzdW3orfHeUR/pdm2KmCMfaKP/woz3m329aXoNzOPkSLQfE55QcretOma4pSM\nvITor1GqWKdp69HbtZabePTS3rzrM1WR3fVkHwdwD2a+Nyp/6YQOednFP1or+oZoo+g/gHp+HWac\nCOBO9XZr3dh8chX95wBcZLZFrRuFVoo+cC6NkKKX8shqTR7RewN4Gnn0JmInpx3ElmWzk7O1aV86\n7fvMPqucT0G12tDqtsDMn9a60ftjD6PYAhpe1M2ugC0Rm6+mD0VvEQs3/tZq5owTMJq2QvA8jHdy\nipVjVyUT5Fo3qQ7jJouDvL/+XNuKnplPZ2apoC8AOLBDdpb8PKJv8mqdC5t3jkd/LvPqRdTHeiGM\n+jNrTUhm/BPz2GCrD2O0OMNEAzO+XqiuuhK9PIhdRc88FqlhVbi3ZGFu1M3VkXSxbaG6fh9GA6R6\nsW6YJ0SKbQ/n1hFdqXPELK+mZOrVb6qTfNrWTew+eJL5PtYmmXFTHf4skGug/XD9O1NjWEKK3taV\nffMPQs0wuraJ3uD3MfnUbYJZWTf25sqxbiwRWiWm0+uboMlr3yhzxmMxGlySOjZUV2fWf0A76ybX\no/cUvRf6mhtH/yJMrumb0w7czlhm/A4zXq3Kqs/b1rrJOrfZ5p0vJmialjE3vFKjq3XTp6K3eaW4\nSojeVfQZ0yzkEr2U+Q0AHpzIU+CR+lStm9B86KsgotMxubYmAJzIzB+p07wUwHZm/tdAHlur/14J\n4Jj6bwKedSOvklL5l2FydfqusI3ftW6YwRSmR89y8BS9PtdHADza5JMzT0iMKLzzVl+qztYH1l9T\nnamx8EovOkqDMPlAyyX6CUVfT0X8bZMuFV6pyxkjExvV1da6sRClmSL6UJ7ecX0o+s2JMnRV9E09\n+pid0pbotaJvcj1DdmqoXV3HjE9m5p2l6InoGATIsSuSRM/MD43tJ6LfA/BIRJ5uzLy1SotXRrKS\nin6BHIZJ6+ajaLCQcCZsIxUPtc1rfErRa+/2x82K2fhmijXyptbNGaj6Y56MeGfsKajWhb27Kate\n2OVSVDH/uR69xosA/GWgzCESiRF9F+smltbGzOcq+ph10/RhJOlvRNV39jpUE/XpfYL/RjXdR9s5\n5bsVfyMAABJzSURBVIewbixSRC8efci6SUEeFFGPnhm7a8H3U+TD41knvJLPxOitG0QU48tG6Bp1\ncyyqybsew8y5C0CEIDfmafWnKHqthHfXCq9PpEIBNan93Nmmj31HIl+97cOoOl2zoF4921o3Go2s\nG2Y8BKMRj0HrhhmPZ8ZX4FtUP6/THIRqUe7cqBuN2OjFNp1mEnUytKLX0HVio17aKPqPALiLc56b\nUE3R/XQArzf7xsrFjGNRXR+53zxMuzPWIlfRX6e2tVH0KetG0ESozb91k8BbUUVUnE7VY+7zzPzM\nSPpbYbJxC2zj3o3RzJXT6IxNjuZkxt5EUZJ5H6r6+Dv4jV8/sD4IrE4K1QR9KPrUwiHeg0CTcOqN\nwBL9DowvurIL+XH0Xhly2oFV9Lact0dNDglbLoQcRR+1buoO7JyOdL1P5/UcAJ9mxsUm7Rm1MFhd\nT9f8vlD+MtWwh2l3xlrkEr2e3K4JX+RaN4Jcon82gB8526faGduJ6JnZUxKR9Lg6ckN5c8tswPi8\nL0MQfUrR5xCh58N71s2q39+qpOP5hpBz0+2TyMO7qTRppzx+S/S/hMk503OjbrwytFH09hW86UR4\nqfNptPXobX9R1LphXg3ZDOUTQqzuh/Tov2y+D2Hd6IXc+1D0oTzs0pIumPF3gV1rStF3xRkYefte\nOOM0wistMdtz5Gg9z4e3N05Oo8s5V651E0u3ObIPqDqJbQeqjojR5fQIcyyUkBnfMPtDA8K8qBtv\nf05dZofBKTTR9bH6lZkQYyGOHvrqjE3FeMeIPvS7fhjJL6dshzFPjA0ZQtFros8d7Qo0I/plszh7\nG6zZ8MomkMb0F2qbZ93YqJuulesh27qJIOTD2319zGOR+3qcq+gn0jHjLOaxEb1AQNEz4ypnfo5U\nGOk0Ff1QRB+rX1GXMWUey9ObG6XJQy539LR37lD+/4rRoLhGx9bzxpzj7Orz2shvXu0kZcYOjNZG\nSCFk3UzUVQ8kD6wTote+qSwMEVPCbTt7cpCybjyErJvdTj59lz3Xuok1xpSij53XevQehKhDv7mt\nRx97GNtzdXlw56AJ0WvkkITMvNlW0X8/sb+VdWMGxblpEvst+lT0VwP4PSdaLvfh3dS66Yp1Yd1o\nUt/gbNOf+6rOsiEq3cbNdw0v8/oa9GdXpPJJRa4A40SfWy6ttlON1Ju2wO7PHTDl7U+WWbWZ3JW5\ndNocxNrJTSZNrkf/EwBPxMh+iHr0ARwAjI0Y9dBG0ffh0Vt4bfQ/MZquWSPa5uq29s6G59fIVvQ9\nYV0o+gmid57EOaqtz7KkRg1qhMIrY52xOchJm6voY2rpCRgtXJ5bviaKPscj9l6JdwM4KPKAaGLd\nNN3fFJLf5yd2TA591yHBwbphBjPj3+CTbpZgYMYVzBNT+ALjEV6x+gu1rxyi70PRn8Xs2i1Dc1Vo\nioqhFP3aibrpAG3d2AsYW4l+CNgJqM5HNQAshph1M3RnbC6hxZZO+0DGeSz6JPqdoTTMuDSjDE1U\nVhOVnkp7OSZHid8vkn4ZVcSR7ozOKXuMiNveF3bRbIuUMJkW0YfyaEv0udc/NGBqKEXvTSo3GObB\no7dPthjRD67oa1V1WiS9h5h103dnbG7UTa7/2VTR29hvDzmKvk2HVhtF36TNpH7XoQCeq/NNjKBe\nYcbX605BQc7v9oizzRuihrc8Y+qc3v7Yvnkk+lw0jaPvgoMBHDFAvkHMo6KPTcc7DesmB66ir31h\nO0Cnb49eE4XnxSYVfSB9CvJA/i5mR/SxkbEejsZoebbO9c+MK4iSy0BqeAuCtCX6ru1It4c2bwxD\nePSfBXBi5nmGVvQhYdQ75zCPrdk8FcyDR28vxLStmybzSzfJbxqdsW/H5PD3HI++Da4F8JN6St6s\nztjE/q5hrBYT25jxPw0Hp/VhnWn0SfRdFX3KukmRdewBJ8emOoHHD2L8jBmvq7/K+sWhyJ5BuUq1\nE3t9psVDg2IeiH6urBsHlwB4rdkWIwQ36qbjaFiN1XIyY6cz/F3Wg80l0qxyMeNHzLi5PkcEQ1s3\nbd72+o66yUHugisWQ3j0fw7gKW3yr+PgY30nckxoepMk6vWLDwRWV9SySIV2dkY9HsQKpKE8+qli\nHjz6UGfsXFg3zNjBjJc1yM+zbl6ReWwfinIavflPAvDYyP558uiHQs75rnO2zUTRM+MyZrzH5JU6\nZ3b29WdrogcAZnw/IlB+A8AdWmTb9OEdCrRY05iFR/8TAF8BcAjmIOqm5YRWsSPGrJtayb86nLwx\nUgpj8GvKPBlSaPBt1EstBjAE0feBPq2bO8AfuDQrj97LP3XOpvl1IvroCRjXYHx6g2mhEH1LHIBq\nfvLHY5LotwK4sP5/2qqty2u7Llesw3AaUyA0VfRDdDbtBsaWWrRo69E3maZY4xOYnFDLQ29Ezzw2\nW6dGDtFrkWDP28f1auPRxyDHXh9NNRt0VfQLYd1MneiZcRPR2KLMet6UP9NJW2R/EaY8EAEAmHEd\n0erqTbYztm+k8p050WdgCI8+CGZEF8+ZMtp69H2+zbSJo48h1JG5FlGIvkdodeL2EwQslVQj7BKb\n2mlgDfPqyjBNQwCblqFvop8FggOmEojNodMHAabmiOkDFwB4fovj1oJ1s4hEf76bao1hVp2xulE1\nIaZoI2TGT5kbLfGl0VddFEWfRltFL1MJ9E70dWTJN9Mpu52v7tx/Y4tD+1T0bSY1a5vvrNHVLv0n\n+GGyawqzno+eAZwK4LAZl6MpsjtjB8Dce/QZaOvRC9HP8nV6GvU1RHil4AHM7tzyfXj080j0nVD3\nk/Q9JmXqmLl1w4xzEQ/V844bAn3FUHfpjM35fYtg3cydom+AvgfA5aIXRc+Ms0K7OuTf5SExNJre\n119PJ1l7mAfrZl7Q1+RXs1b0TR/es7gGp6DdDSUTQc1Tu5kWhr5nikdf4a2YvdPRO2au6FseNwT6\nIvq11hk7ddJkxkktD50HRT8rTIvo25D1EERvl59si0aKvrZqFuGBNYa1pujnhehjiCn6acxeuRas\nm7YQRT+tUdMeZuXRDz1YTPL9eTRV/NjeCJIZ16OKUOoKb1nGdYe1RvRDYi1YN09GWulsGuC884LY\nTZs7W2dfmLYfPS1FPxdEX6MPYfQMAPfqIZ81jVl7UYuo6LtYN1Ew470ZyW7WNNs2ZZkFmLGrHlvh\nRUGE5hPvGz+pP6cdiTH0g0XyXyiir9e5HXxCtHnHWlP0a4GUhu6MTaEp0f8knWR+UMe7e6Q+LaL/\nFKp5bKb9BrEWFP08Rt0UYO0R/ZDo27oZqjM2haZE/ykAW3o476wxFaKvVyD7HoZV9MHO5h6nuw6d\nc6EUfUGFmRC9aqxNG+2QveF91UVs2uNpNNyNaEB6NXFtG7A808K0Ffa0rZuhRZHk37rzssEaCAVT\nxqw9+ia4P6qlx+YByfDKAZVXCt8CgjMnLjKmZd0Ipq3oh0YXRT8rZ6AgE2uG6JnxmYFP0ZfanrWq\nuTcWMA44A4tE9LNAF6IfKqS3WDc9ofOTmIieT0S7iegWLQ5fqx59DDMlWWbcEOisXHRM+zcPaRXN\nYkCY5N+mHouin3N0UvREtAXAQ4HoepJrAX8A4MMN0v8otIMZNxC1nqStKJj2uGnK5xtS0d84YN4h\nCNG3eYAN1W7L/dATulo3bwTwIjQjybkDM97RIPntAVyZyO+cwK4/B/A/Dc5VkI9pd8YOdj5mXEOE\n/ezmoc5n8m/zACuKfs7RmuiJ6DEAtjHz/1KLRVfXKpjxgw7Hno/4QgbzZGWtNUzbuhnUomPGtWbT\n9QD+a8BTSt9Sm981FNHP4s1mIREleiI6HdUarxYvBXACgIfp5D2Wq6CgKd4ABN+khsBUO2OZsQvA\no4Y8Rf05T4r+EQD2GSjvdYUo0TOzu9YmEd0TwMEAzqvV/IEAvkJERzDzhK1BRFvV1zOZ+cy2BS4o\n8FCPBfiXKZ5y2lbRoFBLd86NomeeytKOcwMiOgbAMUPk3cq6YeavA9hfvhPRdwEczsw/DqTfGsqq\nzfkXGOWtaO1g0cIrBXND9OsNtQA+U74T0Sv7yruvC1QIu2C9YaEUvUIboi8CZc7RC9Ez8yEhNd8R\njwZw+gD5FhR0RVH0IxRFP+eY9QWKvgkw4z+Zxzp8CwrmBR8F8JhZF2IAtBnZPWseKUhgzUyBsE5Q\nXoHXCJixA8Cpsy7HACiKfgEx6wtUvP2CgvlCIfoFRLlABQUFGsW6WUCUC1RQUKBRFP0ColyggoIC\njTaK/qMAXt93QQr6Q+mMnR/8DoALZ12IgnWPxoqeGVcCeMkAZSnoCYXo5wTMeN+sy1BQgPW5aM3C\nY9bWTYm6KSiYLxSiX0DMmugLCgrmC4XoFxCF6AsKCjRmveZxwQCYNdHfMOPzFxQUjKMo+gXELDtj\nb8McXnu1oKBgJiiKfgExM0VfSL6gYC5x/awLUNA/SnhlQUGBYP86Jr5gwUDMw0Y4EhEzc5mVsaCg\noKAB+uTOWXfGFhQUFBQMjEL0BQUFBQuOQvQFBQUFC45C9AUFBQULjkL0BQUFBQuOQvQFBQUFC45C\n9AUFBQULjkL0BQUFBQuOQvQFBQUFC45C9AUFBQULjkL0BQUFBQuOQvQFBQUFC45ORE9Ef0JEFxLR\n14no9X0VqqCgoKCgP7QmeiJ6IIDjANyLme8J4K96K9WCgoiOmXUZ5gWlLkYodTFCqYth0EXR/zGA\n1zHzDgBg5rKQSBrHzLoAc4RjZl2AOcIxsy7AHOGYWRdgEdGF6O8C4P5EdDYRnUlEv9xXoQoKCgoK\n+kN0hSkiOh3AAc6ul9bH3pyZjyKiXwHw/wAc0n8RCwoKCgq6oPUKU0R0GoC/YOZP198vBnAkM19t\n0g27hFVBQUHBgqKvFaa6rBn7IQAPAvBpIrorgI2W5IH+ClpQUFBQ0A5diP4kACcR0dcAbAfw1H6K\nVFBQUFDQJwZfHLygoKCgYLYYdGQsER1LRBcR0beI6MVDnmvWIKItRPQpIjq/HkD2nHr7LYjodCL6\nJhF9nIj2U8ecUNfNRUT0sNmVfhgQ0RIRnUNEH6m/r8u6IKL9iOiUenDhBUR05DquixPqe+RrRPSv\nRLRpvdQFEZ1ERFfULohsa/zbiejwuv6+RURvzjo5Mw/yB2AJwMUADgKwAuBcAHcf6nyz/kMVnXRo\n/f9mAN8AcHcAfwngRfX2F6PqwAaAX6zrZKWuo4sBbJj17+i5Tp4H4L0ATq2/r8u6APBOAL9f/78M\nYN/1WBf17/kOgE319/cD+N31UhcAjgZwHwBfU9ua/HZxYL4I4Ij6/48CODZ17iEV/REALmbmS7ga\nVPVvAB4z4PlmCma+nJnPrf+/HsCFAG6PavTwO+tk7wTw2Pr/xwB4HzPvYOZLUF3II6Za6AFBRAcC\neCSAfwYgHfLrri6IaF8ARzPzSQDAzDuZ+Vqsw7oAcB2AHQD2IqJlAHsB+AHWSV0w82cAXGM2N/nt\nRxLRbQHsw8xfrNO9Sx0TxJBEf3sA31Pft9XbFh5EdBCqJ/cXAOzPzFfUu64AsH/9/+1Q1Ylg0ern\nbwC8EMButW091sXBAH5ERCcT0VeJ6O1EtDfWYV0w848B/DWAy1AR/E+Y+XSsw7pQaPrb7fbvI6NO\nhiT6ddnLS0SbAXwAwPHM/FO9j6t3rVi9LESdEdGjAFzJzOdgpObHsF7qApVVcxiAtzHzYQB+BuAl\nOsF6qQsiuhOAP0VlRdwOwGYierJOs17qwkPGb2+NIYn++wC2qO9bMP4kWjgQ0Qoqkn83M3+o3nwF\nER1Q778tgCvr7bZ+Dqy3LQLuC+A4IvougPcBeBARvRvrsy62AdjGzF+qv5+CivgvX4d18csAPsfM\nVzPzTgD/AeBXsT7rQtDknthWbz/QbE/WyZBE/2UAdyGig4hoI4DfBnDqgOebKYiIALwDwAXM/Ca1\n61RUHU6oPz+ktj+BiDYS0cGo5g76IhYAzHwiM29h5oMBPAHAJ5n5KVifdXE5gO/VgwoB4CEAzgfw\nEayzugBwEYCjiGjP+n55CIALsD7rQtDonqjb03V15BYBeIo6JoyBe5kfgSr65GIAJ8y613vg33o/\nVH70uQDOqf+OBXALAJ8A8E0AHwewnzrmxLpuLgLw8Fn/hoHq5QEYRd2sy7oAcG8AXwJwHioVu+86\nrosXoXrQfQ1V5+PKeqkLVG+3P0A1wPR7AJ7W5rcDOLyuv4sBvCXn3GXAVEFBQcGCoywlWFBQULDg\nKERfUFBQsOAoRF9QUFCw4ChEX1BQULDgKERfUFBQsOAoRF9QUFCw4ChEX1BQULDgKERfUFBQsOD4\n/1/ik5hW/pvFAAAAAElFTkSuQmCC\n", "text/plain": "<matplotlib.figure.Figure at 0x7f66a006f310>"}, "metadata": {}}], "metadata": {"collapsed": false, "trusted": true}}, {"execution_count": 9, "cell_type": "code", "source": "plot(range(T), x)", "outputs": [{"execution_count": 9, "output_type": "execute_result", "data": {"text/plain": "[<matplotlib.lines.Line2D at 0x7f667c9b63d0>,\n <matplotlib.lines.Line2D at 0x7f667c9b6650>]"}, "metadata": {}}, {"output_type": "display_data", "data": {"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAEACAYAAAC9Gb03AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXe4JEd57n9f98xJm1fSKqAcDUgICUkgogQCgwkyYAzY\ngMHGJtnXRGOSV0OwfS+2cMAXsE0woEu4wAWDMcGAwAIDkhAGSSSBUNxdpc0nzEx33T+qqit098yc\nPbNng/p7nvOcmZ7u6urqqrfeeuurr0QpRWONNdZYYwevJfs6A4011lhjje1da4C+scYaa+wgtwbo\nG2usscYOcmuAvrHGGmvsILcG6BtrrLHGDnJrgL6xxhpr7CC3JQO9iKwVkU+IyI9E5HoRecg4MtZY\nY4011th4rDWGNP4W+LxS6jdEpAWsGEOajTXWWGONjclkKQumRGQNcI1S6sTxZamxxhprrLFx2lKl\nmxOAO0Xk/SLyPRH5JxGZGUfGGmusscYaG48tFehbwNnA/1ZKnQ3sBv50yblqrLHGGmtsbLZUjf5W\n4Fal1JXm+yeIgF5EmmA6jTXWWGN7YEopGUc6SwJ6pdRmEblFRE5VSv0UuAi4ruK8sWT2QDcRuUQp\ndcm+zsf+YE1ZOGvKwllTFs7GSZLH4XXzR8BlIjIB/Bx4wRjSbKyxxhprbEy2ZKBXSv03cO4Y8tJY\nY4011thesGZl7PLa5fs6A/uRXb6vM7Af2eX7OgP7kV2+rzNwMNqS/OhHuoGIajT6xhprrLHF2Tix\ns2H0jTXWWGMHuTVA31hjjTV2kFsD9I011lhjB7k1QN9YY401dpBbA/SNNdZYYwe5NUDfWGONNXaQ\nWwP0jTXWWGMHuTVA31hjjTV2kFsD9I011lhjB7k1QN9YY401dpBbA/SNNdZYYwe5NUDfWGONNXaQ\nWwP0jTXWWGMHuTVA31hjjTW2zCYdWSUdeeFy3a8B+sZKJsL/FOGIfZ2Pxho7iO3XgX9arps1QN9Y\nlf0J8MR9nYnGGjuIrbecN2uAvrE627s70jTW2L3bGqBvbL+wBugba2zvWQP0je0X1gB9Y43tPesv\n580aoG+sscYaW35rGH1j+4U1jL6xxvaeZQDSkWXB4AboG6uzBugba2zvmZj/k8txswboG6uzBugb\na2zvWWr+N0DfWGONNXaQWgP0je0X1jD6xhrbe2axd2o5b9ZYY7E1QN9YY3vPLKO/aDluNhagF5FU\nRK4Rkc+OI73G9gtrgL6xxvaeWew9ZjlvtlT7Y+B6GnA4mOyAeZciPFWE2/d1PhprbBFmGX17OW62\nZKAXkaOBXwP+Gecy1Fhjy2mPAo7c15lorLFFmMXe1nLebCn2DuA1QD6GtBrbf+yAYfQ0BKOxA88O\nHEYvIk8C7lBKXUPT2A42a4C+scb2ni0r0C912PBQ4Cki8mtoN6HVIvJBpdTz/JNE5BLv6+VKqcuX\neN/G9r41QN9YY3vPStKNiFwAXLA3brYkoFdKvR54PYCIPAp4dQzy5rxLlnKfxhobYg3QN3agWYnR\nGwJ8uf0uIhvHdbNx+9EfSCywscF2r3qXIhwqwiv3dT4au9dYgm5j+79G75tS6utKqaeMK73G9rkd\nSEA/DkZ/IvCcMaTTWGOjWAoscAB53TR2cNqBBPQvHUMawjI1usYaQwP9PAcao2+ssQPcEpxu2lhj\ne9sSNKNvgL6xfWoHEqMfhzWMvrHlNCvdNEDf2D61exvQj8TopSMnSkdOWob83GtMOpJIR87d1/lY\nZrOMvtHoG1uaSUfeLB15+aKukWJi84AAei+/S06K0RrdtcB1Y7pnY9qeCnx3X2dima3R6Bsbm70J\neOMir5Ho//5u46rDowL9NMvUOO9F9pB9nYF9YI1009hYrb/I85Po//5u45pAXcxk7GLLtLHBNrGv\nM7APrJFuGhur7SnQHyiMflxAv5jJ2N6Y7tmYtgOlro1k0pH/lo58cMhpjXSzHCbCYSJcua/zsQyW\nLfL8Aw3ox8WIFsPoV0hHlmULuMYOSHsAw2PWNO6Vy2RnA+fs60wsgzWMfjRbrHvlM8Z038bundas\njF0mW7OvM7BM1mj0o9liF0w1jL6xQTaMKDWMfpmsAfpqGyujl45MSkfWjiOtGttXjH5yTPdt7MAZ\nPY7TGo1+mWz1vs7AMtm+lm7eC2wdU1pV9pUxpbNYRt8AfWNLsca9cpmsYfTVNm6gP2FM6dTZ6WNK\nR4B0EQuwGulmfHZALM5bpA2rR21glkaj3+t2b2mo+xroD5S9hO3zjsrqx15/pCPvk448aNzpHgBW\nW9ekI6ulI4cvZ2bGZIOe6VzgVWigbxj9XrZ7y7PvqXvluMpnrwH9GMMfgHveUYG+kG6kIyukI+OQ\nAl8APGsM6RxM9mlg877OxJjtvuZ/A/TLYPeWZz+YGf04wwrb5x11KO0z+q8BPxpjXg4ak44cKR05\naglJHDm2zOw/ZsnXHPc26UaEz4vwtmW85X7z7HvZDmagH2cj2WNGD5wKLAXMfDvY9OprgOuXcP2B\nWh6D2o8F+nslo38C8BvLeL/96dn3pi0J6AetAJWOtKQjwyrqgQL0i2X0PtCPE4wOVGCrs/UMd3wY\nBIoHyhxPrYnwoEhm9Bl9Wzqy191L9zewW86Xem/x3d0zoD/qqrXSkS66MtbZ/wO+MyS9vflOx8mG\nFsvo/fP2W3CWjly0j7Ow1Pe/35btIuwq4KHed1smC+bzXsfhgwboRbhJhFMWcUnts4uwQoRD9jQv\n4zAR3iwyFiDbs8nYY765geFA+iSGSxYHK6P3zyvASDpyvnTkcUvIx7iB7ct7ecHaMFvq8xyojD4m\nkv4I0LbJvvnb6/LNAQf0IrRFKhvjscCDF3GvQc/+ceCuRaRVa9KRU6Uj75aOPH2Rl74JKO1kJB05\nTzpy8pB7PkY6YhvYnjH6+bWjRmi8YcjvBwrQL5bRVwI98CjgV8eSoyWaJwksqp1LR9QYXRrvrYw+\nBnq/Xlmg75m/Bugr7Grg32t+W7GIew169mMWkc4w+wnwIuD/7MG1VUzsO8BX/QMirIzOuZ/3ebFA\nryvk/JpRgX7LkN8HvlPpyAnSkcePeK/Y9iWj9xuu/4wt9p9Nxm0d3xMgWXagl468cinX7+fm4829\nHuhH6b3PoJ65zyziXoOefW9sLLEnmyvUTWLFjXCnCKd63/1nW+yzmCGmDKwb0pFV5uPskPSGNdR/\noL7jHmZ7g9EvSbph6UA/TgZr87gndW9cc1iLeZ6/XuL1YzER/rcIvz3uZL3Ptk1Y6Wavu1jub0A/\nau9dVwmnF3Ev61VypnQkbgglcJSOfF068shFpD8OW1dxbBfVDddn/0sBel2Gkhdp1HgF2BWcw0Bt\n2DtdSh2sbCBGejh0kWktdmXsgQD0linuS6AfpU0vyutGZKwj7ip7CfCXS0wjfqak4rNwL2X0SwX6\nPZFuvg88P/qtSrZ4JPBri0h/HFbF6OvmDqTm88hAL0KCXQgkuZ9GFaA+wORlGBsZ9k73CFCkI2t5\n2nPODY4J/0ME62Wy2JWqgXQjHTlBXr/qCpHaWD11QJ8y5nYlQirCP+3BpUth9OOyUdr0oPIKrhfh\nPsDNS8rRcOsBR485zbTmcwP0A6wu3wOlGxE+IcIDKtKIwaYOHJdlGOn53FZVgG7dZd7nRTN6Ec5C\na4dmVKT8NKqAwgZlWiqjHwr00pG2dGRDdPgtPOCyy8zn7cBu4G9xTGyx7yqejH04k7sehiYCVbac\njH4V8MI9SGt/kG4q3790ZO2Ik8VxeSyHW/TCGNII8/nY15wlHTnffEu9c3rcC6UbNWL8kj1l9E8H\nXms++89e6MymI6iTaJZrYsjmraqB1uWhDuhHzfMR5n9JuqnJx6g75NTeX4RJ+pOjgNAllCd9Y3e1\nOEbPYoE+noy1+a4bGdT50ZeA3sTC+cAi8+NbBnsU22efSTciPFiElPr3vxV48gj3it+jLYu9yYLH\n35k87K86wGfMN1+6uVe6V54H/MsI59W9iFFihB9bkYY/ofh/B1y7XBNDFijGAfSjVlrrCWCAfiij\nXzLQA3/PzQ9/1Ah5O6LiWOyuFoduqH1X0pGTK+KvxIx+WAc5EtCbYGcnAL8zJL2qtIpkKu45io2N\n0ZvOajEx+L8NPIXB5WidCkaWbrx8LcbxYrE2DqA/TDryW9Ex+ywxoy+AXoSzjYQ6VltygiJyjIh8\nTUSuE5FrReR/LDHJh49y25rjpYpoFj+t8g7ZhVApUtQhf6g2qGIuF9AvldHXfR5kFuitRj8M6Fss\nHeiPH7FIqxZ9LRroDVg9GvgZ8MX4Z/PfPs+whWZDpRvpyBloWWnWfB+lvVUVyJ66SY5To98FvGaR\n16xgtBe8GKC3730x83GLtUUB/YAQBn8RfbfPYp/Xulf6delq4OLF3H8UG0fP0QNeoZS6P/AQ4GUi\nct8h1wyyUSpGXcFWVejLgWu97+vN/xYTu+yxUaWO5ZJubGWuYlB1ALRURm+1fKvR+9dV5cMy+kqW\naRZ2rWRwmYnX2Y6St/j+1jKK58wrpRvpyHXA63E7UsWgOZTRRw160GSsTePY6Pd4vcOoVgC9dORh\ni/Aoss+4J7thVdWbxa6wnWA04rQY6caWxd5k9Iu126Ujr6s4Hj+7z+h/jl4PUyXdjH3yfMlAr5Ta\nrJT6vvm8Cx2udVGR/IyWV3wd5ZKa41UV+r64BgdOc02Z2GmP1TXa2PY6o5eO3MTrV77ZfK164YsF\n+lHfcSzd+OntiXTzHfTq3kENPUXGzOiTzOY7riP3I1xpHE+4jcLo/bIcRbqxXlP2uj2NWe8z+iuA\nO6Ujo6wCXzSj91ZUV7Wxxdb/AOhFuFSkciHWYiZjC6AXQUR4yqiZEeFVIiOteB+Z0ZvO/wiqO54S\n0EtHDkGTzW+rjSrnQPS6EZHjgbMYHugqtt9f7K1qjldV6Phce05KUrTlgQGqpCMzdb/tBTuWdMHO\nzi+ndGPT1c8quV8me6rRr2Yw0CcjFukojN4Afd8er6rb/rEY6PVv0m9JR06gOt/+Pa0b5kMJiYQP\n9Kuj60bZvnKQdOOX9UtHSGsp0k3VSG2PgV6ENwOvIAwPYdPbU+nmWNwE5yj2MBg9htWIE772nVYF\n/4vznqGJ8N97v1UB/dhxZmxALyIrgU8Af2yYvf/bJd7fBRWXr684NvB2NcerGH0p5oTxXmghBdD7\nDaiqgX/Q/F/SCxg9HGlxWq10Ix05re4i9ozRR6xzKKMfRaOfYqh0MzZGb1IsGH3Vcz/D+1zN6J/6\n/AuAX1Cd71bF51jr39uM3toodWlfA30bV45W3/fT2BOg96WbcMJYmBThxkWkVWc23Z+PcO59zP9R\nZNYcOCz6rQ+0pCOvlHXyBLgEph78XDlPfjpiXkeysQC9iLSBTwIfVkp9Ov5dKXWJ93d5RRJ76goX\nWy3QG3lIAfPoSuJPxg6Tbo4z/4fFbblg0O+MWt4OYwcx+ivi29fcZ1RGH4KTnoy9Ej06q2r0AzV6\nY9OMh9EPA3qFTSjpj7oVYjXQT+y04Kzved476+6ZBueFx+29Y6AfBXBHnYwd5b3a85++BzHPq8pv\nKdKN/T/Iq6jK6qSbKoKxFjh+WFrSkWnpSGrDeIhQR5pGWYFrgbsKe+o0egjj3bSBv+blAJfARQ/8\nGU9cVCTeoTYOrxsB3gtcr5T6mz1Mxn+ZI2n0IpWjgEENaQK92GgnegFKHdBXAZPNn5KO/Ll0pMQa\nzMTj14bke8TyHgno4wm5pTJ6WwZrgTnjXqnQFXIQ0A9i9DMMA/rRGH2VdBN3zvo+yUBGH9678nvh\nGaPTe0wwx+aXgwXROG9VjD6N/g+yUaWbUcye/zQW76VSlde6xU91m9NMetdYYPPTeK905LEsTqMf\nVJbDsMOmNQu8BdghHUmBH4ssKnyKb7aMRwX6XvSbL91sAyBdGHtQvHEw+ocBzwEuFJFrzN9ioxHu\niSRS5YY5SLqZRANTFdDXRSKMjylz3+MrzhmlLEdl9LY8Rqk81sbF6FejG0Ji7rUUoH8yQ4Fe/yxC\nIlLylDnUrCYchdFboI/dLOvsl9F3fb6ScMFUHjye/8U2zjhvg6SbPV0BuSj3SrNO4PkMGKmKcLII\nrxqQzEjSjQHqus1pfEZvyylO42RGkG6kI18w3weFkx4V6AHOGZDewHTMJLAlYYsFeruy3Wf0dt5A\nk4bWwtgnZ8fhdXOFUipRSj1QKXWW+fvC8CsD2xO3xSoWMYjRx0DfWoR0Y489kvoQBABIR26Vjjy3\n5udxMPoCWAwbKV3EIoFehE/j5Kk1wKyRbiyjr1qwNKof/aBO3JduWiYvfn4vBb7F8MlYJ91IWbqp\nkS3qgk6l4e/BaVUT1IMYvdXkB4HTKLZY6eaNwPuHnP9HwF8NSGNUjX6Qh92Ed80ggjJKUDM7iTuo\nLBMo1SHf/PzrtTVvmbftuRbopSNT0pE/8tZBPBj4VHRdJdBHeclxkqEfwfK95rN+X2l3/wP6Mdli\nGL0992MVvw1n9L/2sg284pi3Ekk3Zv/T59Tkxa9sdXEw7H3uA5xbc05Q3maZ+G+Zz48U4ZOAr9FP\neucq45rmM0g//ksV0N8R37PGLgYeYz5b6UbQz30G8NGKa0bR6IeZL91UDckHuTrGo7CSdCMdeYR0\n5Ac1eYwbkwV2e09dbsGcdCWjX4x0MwqjH6bRj9JW7Dn+/eJ6MGy+YFSgH7RvgS/d2P/vKeq5y1ep\njkpHOtKRT3rX2Q7bryd2/i12ja2rk35no4E+qXTIiDuKRwF/h3MaWYlbUzCM0cdAX8Xowzy0FsYe\n+2Z/AfrFTBSVh49SuLdVFVAI9Gd85DDW3PoEyhr9/YEP4Ves33xGYlizX0FGCSq2o+acuLy/DVxm\nVu4+E62l+knFYYrXET6/77sr0pHjpSP/gKvotzB62dqKGjP6QbFeBjF66688qI75jN7m+ZqK80Zh\n9EYDCjT6R6M7qjiPs5SBLmaKVfkeRaP3F0xVMnoRzhRZ1CImX/qxADEKCw7ASzryj9IRG/t/2P39\nMA62fVWx8r455zcrfmtRlm5WUdRzoAbo0RFln0bY3p5JWJbx6G0Y0JcZfbXnXWy2rlgdP8G9/2GL\nG+PFmF3vM4RAr9NMugct0AeFJMIzRLig5twqVnG9+T8c6J3+3Y40+rKGeN9P/SNwU3TPOgbjl+Wo\nQG/NSQ9gWaQiZOwAq9l9qA/+Mx6bEbQm/lJc5VsYcM/Y7DtYBcyaTIT5Cm0Y0N9k/lcOQ43GeXYF\noz/dP838Dxi9dOS+OG8HqNboE1xjihv+Tpsv6cgl0pEXF/dSMgjo/WcdRbqxE6BxJ/J9oC5UyDBG\nP0o0UvvcbewEn07jYsDOnw1j9BuMgwG4dwmAdOTj0hFbD21+3lP8LjzQfBzm5AC63AeVtd/ejiIs\ny5QVW6AM9IPamTUD9JXybdyJ2t9mvN99jb5HveNEHdD77pXhfdL+gQn0IkP9h2O9/ePAB2rOLV6W\nB3K2QQ0HevcOp6JhW9krQNTDcX6y1uoY/Z4D/Rum/56Hvv3EKMubgMMjje9V3HHGg7zv07gH8kce\nPtAvltGvJJyMrWughUYvHfmqdCSWq4ZNIOrj1Q0NEc7ltnPub77G0s316BGYNc+9shLo43qxy8vX\nRnR0zIFgIcJLufXBfsyktpESBk3GWgZYJd3ULeEfBvSVK6NFuE2kiLrqSzdfQUeKFMK0hwH9e4Fv\nxLcx/58BRbhvC1S+pKlHZUd9dx1uNfJiNfqqUBQ78aWbs957GK85AoYweumISEeeBVzoHY6lm0G6\nuE3Xvk8hZPSz1Es3MdDHGn2Z0R/AGv12EY41OvS3Kn73C8m+9FG0yLiDqCqgOkY/FYFMlVeATe8R\n3rFRgH7XCOd4d5l/Pqd/9FeK75rR287Cd4tre0PNWTRYVAGbPRYwehN/pu6d24YvwCHcc8oL8Jly\n2XyN/kLKm7IM8x03jbHE6K29kntOPsN8vtTkf9B8gAYcCbxubGFVMXo/Xx7YiF+evv0DV770Au+7\n9UaqYvT2Wls/F+3ZEdlA6cZEOzwK+Lo55HdwPXQhWykO6cgpPPbV9/Ouf4cIz62YtD4r+u7/bl+c\nBaquicPjQjO88Pyn4cqg7nnrpBsLrv51uzjpC1YnT1m5ZWXxObwmft9t4CNI7k8c63NHY/QWB6qk\nmxZ6L4RRNPoM9/7qNfrkwAV60AzgIuD8it+mgHglWB3Q2xdyB2X9eNCQZxLowjCgz30JZTEbf/hl\nWZf3+vLurvQYkWDyswBMBPKMy/NOFs/ov0O4BN03v6Iebx4hZ3TpJm5cowF9WbqxVrVqts5fO8cC\nbtL3NXrbmOJ6UUg3xb0K6WaARp/0/Dx20c82iNHb/Nrvl0lHzqM1B7/z6CfVPEuVDWP09vmuNP+d\nPKnLJQacP+Rhf+2PiF4OvIrhnU/V77aj66IX8X29+CUJdimr66SHAb3/2wzPfcLnSLs6vaSv69b0\nXfHosbouphVNt1qjr5NufEbvQqloeWwVZYs1+nNxI6GydGPxJu0f0EB/JPUeK62K34a5XG7nVUdd\nKh3xNxQeJN1MUGL0WVnDTfrmmIIkr1oZV6eRjrJIqb68eytcLdSh4H0fdsdYHNBfx+f/7vWEE36j\nSDd1k3D+8QVzn1EY/UhAX8EWBzP6kz+/gdNLjlV1i1oUBaPPfUY+SKOPGX1iPtUDfRoAvV3oMmgy\nNmb0AM9k3S/ghK/5Etyw8BjDJmPt/bab/z6j7+MYPd7x2OJzbL580KkiMzY/tv4OYu7RpcXxqmuq\n6pWeFzjiGn28vUuX7/Ffnx5wjfvemi/fxbWnMrie8nmkI2dSBvqY0d9JdRiXWLqJfwOf0ectI2ce\n2ECfoMMPIB05Tjrycu+3lPIk5zDpZgurNv0W8DI/nQof2kEave1c3HDbVoa1v6y7bwn4RFC894rn\ne4cGAr0IR4gEeiEsrHIdnSoYvQX6Fq15OOaK9aZiPhTYxPy6hwNW269j9HFe6phVBPQKBk/Gxn70\nxZJyY365+9/96/3DzstDEM788IWUbZCureuP5PZ+o2r0AIeSLuhy2nS23fVIf7/Dmxt2AdNgCKM3\nwN2iDKBHhLs0FibRf6Qjh0lHfCnIH6XE9wTXecXSjWX09l3WLTaqytj/iM4J7arftyODvPYcgAs2\nuhDNYWjqYZOxT/WOaaCf3AGQMrFb14f1N4wG9I9+o/6WtV0Gqt0rtf32EwG+QDWjj4Hed5LwV7/W\nYUFZusna+v1tPuuimmv22PYJ0AMvBt4hHXmWdKRHqHWOqtFvMv/9WhPu1tKR+3NJ4JO+4KU7RdL3\ngV6f2FrQLyCrnauqDhPcXeV7iyTm/nEitsJ81fw584GeAOingBbnXwq/94iHmkaSo2Se9qyfn1En\nY0cB+nQR0o1N7wTCSWgLtFWs1uXD+amHPurVYDgI6HX9Oe+dFlB8ULw5Oj9m9HDGR7QLoQMhnYHt\nx7lzQre3OkZvpRu7KjQnLPMiTG80XxIt1AK0PPknLA3oBzJ66cjvclwx51pV6P4m2WWN/rpn/l3w\nHdpM311O5YI3u9FYEjxCHdBX1dO13vWtAuin7xkN6M/7B3NHT1Jy0k15r4C5taAXCw7yukmBu4G1\n0pHfNscsNinCuSLfyoxepboNVklMS7TlA/oHfPg4/uCcZwKgivt+BNcwBjJ6M9HjH9tUcV6bUALy\nWilnAT1PuhGkbxupz+j1RGo10JTy5TIYrJW3F/9ldJYFnLIklLddZVBiFytlwI3A02jvtvcByOmu\nyGjNgQPoKqmiS/kd21HFb0fxPSaDz4uXbo6suM98lD/fDOCVGT0wES1UsjYc6M9+ny+91MXujxk9\nLKyy0o89ovObeDjuPHpgONDPoCfpYg267TVkv8zrJoGPjH6zGfRHT/aZLfjYOnoevkafpxKdD/Be\nnlCQ9qpKP1i6qdo45rj/rEjGMwley8noubvYquSltd71Ka15XR8md0xF11QDvc1vAPRF/mMPI8iK\n1xNv4BJLN/aFPks6chKO0E6bc6tcsssafZ7q99cax97koS0f0J/1vvM46mrtvZK3Yw3KB3oV/bd2\navT9nprzfPNr4euAZ3rAMh/50euyePzLNTMfbecjRHTMOSQY1ttyrdpkoS7P/jENUis2YwD+/sXE\npAX63oqEtAehC18M9LupZ/QfJvSU8RnupLnfYhh9zIgS83vJ60KEKfSWaT6j98tvcpGM3k3GOhuk\neZcZfXelLjvHRhOyFhz9X0hHdFyUpOeDTxc9inlUcUMts/hAP0uZsaqiIf/NL94kUrgf1gG94jee\n+Qjvty56K0Tf7Tdm9Pad/Q6+182uw20oC38i1gfeYUBfZvSunaj4p1pL+r5X2mI2KdLySGKAPjVA\n3561dczmtV6uTBcgb1VJN2VzdTDuQGL3yj7wWfQakBvM/b8E/AT3zmKLpZvdKAP06YEM9H5h9yfj\nCUGrJYKtTKtvTqQj/gRHXFg2kNIgRNY17oSv2sITj9HPkxSTsa4x3u9TJ+gjtcnqZddu2P0Oc37V\nZOzwMZi9z92nXuBlW4P2a448mqc+T6dnzysY/cqUJAB6f1JrAvgbtAdEDHh+Pv2H9N/JHR6jH6TR\nW7YIZa8DAeZZWGUDNiUifF+Ei4Fp/izVIYDTrmU/EaPfQ+nGWUK9F1aZ0dv39xCrRPBB8jas+yVg\nXILLjD6YUEW/CzsZ648uw1W8tiH3p14H/F50Tvzgqzn945ea12Df8fMJgT5m9P57rfK6Cdc86FFL\nnYQSAL1X7xMRr14uxiTztZ1he/OC7jCvpcToF3TdT3ujSTeg2bKv0Xv5L23KXQZ63+XTesPZCfLr\ncBFlE7SHmx3NVTH6XvT/fajkIGD0ARCW/KHL/sivPO4U4MvekRg0Z83/4Yz+dx4DD7jMHnNA75iM\nv5wa4GLShRtr0oxfurYkr2L0dUDv8mzzEKhS4mSH1beZ/JnfNQPJ6a6MGf1fMr/GAuE0Gpxi9y6i\n7/5NfdfFR5ifhkk3Gfa99ScPBeOj7e4zT3eFX/nPRI8iJklyOP7ycmNqzcMj3np0DdDXhdpdLNCX\nQyDoHbUn1IVmAAAgAElEQVTCZ81b4eRlyOh7lMt2GsfoffIS5sM2ZD3xH7oR7tqwWoQHSUfeFRzX\n8zGJSXsLYZjqOo0eAo2+ZpDj5KlBk6L+/TH/p7x6e0txlih4UeyCD2A2BUl793jHRgH6Xwe+h+3o\nLaNPjLaSLowu3bTmIZvwN6rxzwvrRD2jT7zvFrv8+cHEHLuI5110JtXPaOuG81pScjAweuUKWyVx\nA6zS6CGUPixo/hI96WeBfhCliHvu74eMvvCZj3XUHklWVwHb0X+b/mIYvQNYO3SUQCVxEzi6srUj\n6SajP4lh9Bbcj+L2B1kPnGlzfczk/LwBnCSC3TjXfyd3IJkF+dBFQu/iozdu8YFeJTaAl10PoYE+\n7fl+7fY+uoFObfOzp+vHRa+Fx7zpuzWgVMHoiw5pMUA/R/z+9IYlYRpuQt7o9XnIzMt1zwK9vbdl\n09WMvj0HIXDCTy7+deAqtMMCRT4nd9pzCinAY9e2/QwC+nopqx7od1H2L/dBb8ZjxD/yH5Ejv191\nJ90e0gUH9JseeGZtvpzZ+Sqzn3Gu72/nxWYP2yDCyYwC9OkC5JVeNxDXF/e6C3Inwou8720co/c3\n+Rbz/SiOv/xfGBAjCNfo+2D96A/kyVhRfuFXvYhBUfDADYNvBp7CaEDvMdZcAed61HnegBmEsgdA\n3/Oxj82+zNdIR7YyaVyXq4F+WNf8bcfoM6I5A/1FVza3uGvDdQA52aSKGD1QePlYoK/yjfbL/kQq\nvA2UIkcQVKV082m09mi9mHSFbc9VdSjzhdsiJDzibXDSFw7DAv3kdjvx7PI1Y0b1o0o3uqEaoC/5\nZw9i9FVAH9bBLJpKkiwG+rhsYkZvQbYVXGcZ/XMeDyXiUBqgGqDfAY5RW+DzR5dVgdrAulf2J6Q0\nwb3akHAHdlVAH2v0dUDvrq2XHhbM71uLI/Nrj6472bMI6A2jR+m8Xfub70bPWwxb0zFQuildVy3d\nvBs7X/jada8C3gz0+Y8/73j38UmDgrwKS+I5yV7B6A9o6SZg9JXSTajR2y86ToXtITHxxnOcRj9I\nuvEkkjxTKgCDeaSW0Wckeczo/fghAH8IrOXRb7QgvSfSTVKwKck9dq8sW7aVbTJ6zJz+pIo0elDB\nUu06Ru/nsy4mD0ieoZICzLwFPeejXe6m0e+gbuSjJ2PTni0L4TFvhOc+4clYQJra7hrTzJ36WN9M\nFYwK9GKkLOhX+GfXLTzxQfGXwFaSrEw28ujy0I++aqLaB3p/VWoI9JbRr7kF4DARFFf/nvZaSrIw\nTftOJ3aCq6eWQfrA5j+Tn0/d2cyvLZfFK4+1z1U8YXhvmWXHUc6b6vazT2bzGXYuJwZ6V8+qFiZp\nM0A/5zT6tDtKKAgL9KF0I0o/f9qfRkdLHczotx5fZvQygNG7OhivuNXkaHLHi8z3LNqgxgG9kpyk\nX4VRttBtXmYRpaXJA5vRB4w3ZlotjrimLpb8VeggZ0arnGuhC2cgo5eOHI4v/VifefGB3mMyX3+D\nP7GWISXpxl9WDnZXmAe/E86/FH+cByMz+jSQbqoYve4fJ6NwADn9KQyj98CvYPSTeIzexLK3nhbx\nhGK1SZaZih5LL1Ynt0Bft1JYT8a2FuLr4Y2Teu5lcgdgOpA/2aB98qxcMnNnVZoVjL73M+DtQD8a\nhltGfxluYwdtNzxuLSGTnjMEosect+4la4cNVILJ2LJ0c9nn3orTbYcz+pseDjZC6RWv+765RwT0\nxre6zOh9oG8RylEx0Od0V6VDNPryZOzuw2bYfNbDiu+3PPSJ3PxwC24DGH3dhlN0eerz4Nhvup7A\n1Q8ARDgh3mkM12FqLd4yejFDLt1xbqYccoLge29Gs+XZQ5zjQMToRTij+Oa4aQvtMmnLWwO9Siwi\n9yNS4M/P5KS9KqCPGf1Wkr5ZJzBKgNLF2b6RblQS3/fxvPjsp3Loj6DMks5Gu7Dp69u7Y6AvFaII\nr0fH/fhIcTDJ9EtxPrShRv/TJ+swq73pHOh6HjlFrs3/MjNac/NipRs3yecz+qLjqWD0YYXUQB8z\n+jz1vQ8s0NvntRU4FbfMP3yWc94FD/5b80Um+MnFF1MGel3Zs/YqdOUv1UoRVpi1Ej6184Cgq938\nEouD9sLc+S7f71NUWBno19/wMqX4DLpz9n+xQH8tej3D3xa//Ocb3g+0vI0sumZE1gvS8BfNPfEl\n8OB3vpBQWw3r3twhv4YqTcYqqjR6JYq7T6X4zYJ/zOjzlgavyQpGf8VrDzHAVDB6EV7GjRee56XQ\nQwlc//Rj46IrrM69Mm/1TR0z5yk4711vNd/MXEswr6StXQv0C5z5ITjjMjdKKLP/eC9k0GXupJuk\nkG7MBP4CaEZve+kUQDryu9KR/yBrpeSpJkateeh7fgchOUiBHxTfQulmnpjRuwB4WSTzTYOd+5Kc\npDdIurG/bUUyA/TDVOzF2zIBvSIai1drp7oyV9GOw8BsLjJRAP0g6eZtxLHcJdMNNG/b1htq9BZI\nb3rkTmC+4uXcZv6XgT5v1QF93Rvz7tvXNS3JqqUbq9HH0k1vxjJ6V2tVywK9DX7lpBtbiZysEF4L\n8KSXwuNfYfKQE2m6EQikK6ln9LvoT6+gDuiLIwWTdN/rVyRDFdCv3Gw3OMlrGH1fbVQ3AK8uftGN\n0npJ6AUv2nUyHBX4LO2M/wNpb4rqTa61iVmxGTL6CummB3krM3VOP5Mdrks/TDc373T1rZDO29GC\nZvRXvvh/oYHJLtpRwDu561ce6qXQZ/dh6/j+C/6dOquTbvJWz9SxIuuefYiXPOAPqxn9PMyvgS/+\nte9dA1bKXLXJHU+7cEsQ57CqPvXwPcgklm4WQIditkCvz7v75FcBj2F+7WRBjNIFfyGUz+hvY/Bk\nrL840M5r6bo7t64VSTczuAi2qobRx5OxW0lMG00PVKDXHiWuJObXrin/TpU25SPNuwE49psr0RU9\n7hF9Kw+rRdkGap851Oidn7peCFSejH02ZS8EHXWrHugBvliRP2spSd+tyCwqnaqQbmoZvau1lv3p\nPFoN2XiM9Nd5vz3GfHajgUkj18/b/qDkqBEOh5O+lQuqNXqV2AVVxRWlc5I+qKQdfO9XKXiFzYgU\nHa62Y75tQUMqGL2/ctX9qBmd9VIxjD5LARUxepfnuWJJhx/OOiQlSR+kkFfqpBt9Xt7KTKdidFlT\nVKIi6SbRBfLEl8Fv/sYT8Rl90rcdn3Xx21ak76yPSpOKjjvMj/kUHM/b3SCtsA6u4tAfPynS6PUX\nyXSbmF8Tv3PdwFdu2VkcSRdgZ7BmKovuAw7o9UgoyeClpz+T9TfoUapm9IqI0XPPyToUc3+6XUid\nrYWwjrn33Seu4y4jmtG7ucXVJi+6QL/zh88rzefYuq8kr9Hc44WhW5FM4+IBy+iTDC8qINx6frjk\n2Wp6k/Vzg4U95G+O5Nd/50GUe0RnK2/PiSMdTu6MT58PKrgD2cT8Fr+d7cDthIxeN6y8Bc/71Ud4\nx/3Z9+2UzZ+MtcG4fJ9650evT401+oz+tJie30dGS4fLjD7pzXi/fc581mV09H/B60zfO7/G5oEa\njV7nSQN9pXRjzhFCRl9GmaQfgocGQP15oSrqKyuIV1Kef6ldMS01jL4HoDYaAL397O+x+YEAPS7/\ns8NQBuiTXAOzn4Y/HHeMzZ9EC4HBNehBjD41z5mbZ18ZXBtP3FmgBzjt356IXpiWA4dy6r/Z3sdK\ndVtNGn4KJuyHl+yt5/7IP6Ee6NNeSboJ85ZHjP46dh/2laIu96difDHzZF4Zp12JRnF54cnmzEo3\nOk+Sw4brjmXV5vubNECX+Toe1YE3Tv0/edPUWznCuHhmrVYto3fv2/diMs/n7UYHC2w/xs75mcpp\ngL4/Mx0w+u4K4bPveok5R7mIuGHqxfNq20GSmbpwwAJ9n2Jm4xKB6WhEN2HiuJRfcNmOvGYFD7js\n7Wz4gW1k5UJ8wQXxZAxc0NG6sKuYC6F0Y/3UlWaiaXd3lIKd+feBXr+RcNim09NW5zbqa/SO0ReV\nzmP0p3wBjv3Guqg/U/SmxDRCr5UUo6Y2H/vEn/LJD724iCvkGrPfMWjwX+FNfC5YoC8YYBnoWwug\nkkxtVHZSsGy6PAdLN6FLqW6INp9l9qmoWjA1ufNL0pEjKDN66woY5m/zA280j9TlqpdcS39qJbou\naLnMT8P3CD7kBvvpM+gtG8tMvbUACoWSFvOrp+iuaJEnEp3XMozX3mtFcS1Afzoc7apKl0kNSk94\n+UPMdzty0cRDp3ut+c2L4WPK9K773jdITde78mRs3uoGoBMz7aydeuA/CZyodX0jQ9YBvdvAXYN0\nqG+nTJT27elhO/jeTPye7WgoAVZx7BXQWthAa+ENrNps8jkxUTD69Td8LehY3DOVGb32OgPdxnbz\no6c933tWe0+YPWQ+eIaFVSm7jtRhMzSjr5JuYkY/C8rMQRyok7FHXbnJsFRtU1vD39tmXnVqO4TM\nr871SnHx71n5IZJossFDn8DrJgdlQokW0k2mGX1rIQb6jHKvr2lYBdBLR56JnlcY1D2npD27M5LP\n6J1GD/DElx0fNbJWJaOX3HmS7LzPE9h+7JONhAJpsarTH+noz/70yYINb6+IXBxDr4q8ZZ+rplbm\nMXCU61pawegLrTpqG1krY+cRYeA0p5ZNUM/o40VQbul50p8xk/NdfW7eJ/HL2cvD7mKOcE5tVJ+j\nPMkKG64FlWYoSfnuH/4qtzz0PLqrVvDDZz2aYnLOMPqsbSWK0KXOrSS25VC3F6lvb0bXF42QGij+\n3vy2DZWooDyzWnksAX4O/Km+SyuUbuJ2tbB6VfH+lJyBkgnm1m1yjH6yWrrxy7W1IEb2eDzWs6XM\naL1QvhPt0qhLd5J6FFWWUEClLbKJ3BCJJABl8Rj9Hzzo33ng+70Lvb0rYLdXFv5cGMyt7wb37c34\nJEGRVGr0MdDPIbmWXnU5D1oftGhbHqB/9BtvDOlRZDYy4+hLf3MQ53Vz8e+6+PEbWzY+SbU5YDFe\nN5KRJ6kDesPoWwux60CGrkxrvWN1jL4NfBS9sGsQo09qNHpPukEXXQh8kxrouzkh0DtGn6c6GQv0\ndvn+1hN874sn6vS9/tSyHccAv1nk1VprHpRpjXnVa1UgSgg38I59tL372DP6bthaYo+TLbYdf3J4\nrMDAFsM1epuO8W+m640gNKOXvF+rY88W+OuvPtVlev3TYevxc1z0OkiyFkkOkk2QtfQipTtOfwJu\nk+7VhXSjAWuKR77Vxj4vP3d1u4nnRR6IL2/odG393QyoQKPvT1Fjdo2KrmwqCb1uYg8Zlbj8zh5y\nGJvOatObmS9Gp/2pOO9dL38uTf0ev4EdhZWJmgsvMLkTkp5LoLvC4oaeqC63RchbbfqTirSrgd4v\nUp/RH/W9EzjLA3rH6KeAXV6+jZun2ftgYfVCcN/ejNcuRN+3PxmDfR79ny2AvrWQM9w1e1G2PEC/\ncvMh0crY0Cyj1wU5yko5RdZ2s9pnvR8e/HcDL/AvNdZFMoWSjLzV1q59xaTqPK25WfKAIubAfcH3\ns61l9H5M9CrGazNxQlFpywumXIvPE8Hf4vBtuz5AbzohXYiB3jF6ZVbLK+MCZhn9Ly56NvFCLp+5\nFwCs7Oeno0HKnZT2NAgA3H3q8aWnk9wCfV2IW9cQfGBLe57OXdKD8bYJNGkUsWj0sw3Q6AvrT+rv\n2aTohp8L0EPyFmkvD96lXy5Of7eI7xj9dc9weQA90mh1W6g0B7GdiR3GPoAzP+Rr9HDmB919YqAX\nUrLWbHDszd2XhicV6yYM0PfByWZ3mJfplwEVdjqfee+j8ed1VNIP2PUgoBfVQiXWMcF4UJWCF5Y1\n+iSz0o3eyOWIa1ZXMPo+toNtLcDELqf9dldGjD6Je0rD6M1KcsmS4B2HGr3DIvA1esfoX3IGxc5W\nViq9+5TdwSihPw0qMWUo2o9+23FdLlG+ZBYz+nkkb5P0/tWoH2MV6pdrMlaCVhNX5ond1ccHSTd5\naofCx+sjIz5K4Uefm6Gz9Augz9sgGUrRpzU3i0p9tCl7l+QtM0FUAnqrJa9l0Av7jWfBWR9YafKl\nG0DWspqJz+il8Ov/wW9Db8VD6E0ntBYynv3ks5na+gvg1R6jl4LRYxm92V9TT3z/MCoU9zHpueG1\nSjA6fBjASzJdgQH6U+Wxsms8vqb+wuAcm1V/pJL09bC1P1mWbrQfdFgf+lM+0Pdr/Oj70TXGzbaV\n0ZrLTWfWJ8lTkn4ejFB8buJ8w+2Izk3G6qKQ4LrWfEreUqaDgjz9IVa+WX2rXp1p8+s7IZT3yU3I\nJkMZMW//bXTOBH4kUZ2ufe47sIyeIYx+0wNfgU8wlPTdnAk7y0CvHZX0PXPNlPO2HlklWdWkonE8\niNp13kZtVBkgvPjsL9Uwel1/r34hTO5wQN9b4Wv0KapV0U7TNllby1dpNw3ecajRO3Uh/NEB/eHX\nwqrbZ8zPLd62C7YfvzvAgP6kNxoQRdpVhpTEC+7AuG2rjSpHpT2m77arnutWnO+RLZd7ZajXxhWg\n2FTDPNsh8T7hJctJspipjZYXOxEkeUYSM/o2Rg6ByR07PcYIVQWvUiMDlAYrltFPM0i6Of1jcOYH\nV5n8mAnQNNMunl58DJUKYr7bStqfEtKFnNM+dyTHfGuCpJex+tb13jWakaskJV2AM/6PDr+sAetQ\n9NJ/c67X7o75r28W+XHtUc9jWEsyUEldJ2feoxLCODqvD84pgN7r3E/88gtJu9CbhpjR5y2QLGb0\n9mLtFVHeuagC6Kf1+1BJRnuuZ8pT+wMm/X7wLn1AcEzPTpa6yVj9DElwXbqgGb0SPRrozewAtGdW\n1tarbu2z28nHPC2Tnbn1h7LriIGLC4zFjD5TG5WojWoWJcoboVVr9E99Lkzssm652vrTu4r5tGzi\nw7WMXolC8rRg9LZzKC+Gil2ctXWjcEvVHYQmFLs3QHvWyardFT6jTysZPdJCJYqsBa2FkNFL4eUU\nM/rzILcZ0dJNMYdmJBtRqdHms0Cjz1sho096dvQZLpwD1Eb1VawXT97qsvIOS6wOQI1esiSIBRPP\nKtvCtZX84heUpt0jU7Rnw4IYtIFAmBf9AtKuBbMMlbZpzd5G3nbubYf9+MbAj7qS0ac+o/PNZ7KD\nh2DpvBT5TzIzpFdJsDI3T6Xw+bcjl2yyVbhtrdiS8uQ/OI8N1zvZSzc6LZ+c/jF49J89AbCMfi1u\nQQcBo291zZxBMBkbuhLq96fvXbXASUs3CfCi8o82RdPYZu50Lj9Peuk/k/b00DcuU5WW643zOZ9E\nS3H+r9UafW/Kzi30ac/2zagt02UeMfqgXAps8YHee//e5t4qgdZCizzVjL49Zzc2sYv2IJvwgN4Q\nnWwiLz235LDzqLo4/ObptwN5RmK8O3Q5+Q7wKki3itGf+WE4/IdT+NLN7KF3M22mFvQoJbzGAX2m\nGX2i5c+0pzutXaV9d6zuHEasjIF+EKNXKaTeYsZeoNG3UGkFIUvaqESRtwyj96Wbvj0/AHqluLKS\n0QPaFVdpuUqn1afnvaK85REg0ZKRrmdVjB61Uem2mLV7rNii2B+BXkQeLyI/FpGfichrq0/KJVhQ\nFDdYW9EtWE/fUzVL7ZuiNR+eE7tc1ZnVedOuchp92iLtzZFNUCxBP/LqLRHQVxW8vxGCb8OA3uW9\n1bMjDIzbXQ4qCfypVZoU+bKMM2+1itW7K+9IWP/z0FsjT6E3rUBNBfnTjH4Nevcpk35QDfR9fvUV\nLwOzyCxPFN965ZXSkfO173vmHro3U35XScHory395s7REtH0PeHiCcvoY1ZXpdFTAP0EenVrcAcq\nNXqf0e/ukycG6POUpJ8FjN4ybSXQKjR6X7rxGX0o3VhGj2iA1GDmRkHZhF7489jXeHmb1HKi714o\nGcyvqfe3kxxetxae9NJT+LNJvbF5Ekg3OD96k8VeuMSkYPozd4eMfsfRLvhYnko90Cc5KC2J5G3t\n9punsOlBcEnQbGyde3SQTgno+75HlX2WieKe7VlX1gurLDhrRl8N9C1Uqkxn1QrecZLZctL/JwKV\nzNavafzJWCtNKVGm7fS556SwXLYdr9ObX9Mi7Skz+vRjN+8ktmyiv18CvYikwDvRrlH3A54tIvct\nn5gnpD1XunHQI5/RJz049Cd2tUy1Rq8kZ2KXA5hvvcpW7gcPzXTStw031x2O5KikRdrVvrC2E2rP\nqkiWKFcgUUcWaYY2OtAXaeWW0WeISoOgarkn3RSMfiIlNZ3BzJ0J0/eEK4xUqkFY1FQA5K7sHZqE\nniaGFWaqOJ63Ux566ZHAtzjrfbaSG6CPQAN0QxAFb+7Vew7Y6yVyc7GxSETBV97qjmtJQEx5wCcv\n81l/HaOv8qN3K2lP/ErLzKdZoA8Zvd1WMJQFv2OfAMvobz+X4JXmKaTdFnmas3ILHPMtWFjlIrBm\nHqM/8SvuOsvon/YcdyztQXdFfaNPF94JwDnvcZvTG+nGOytk9BZYdVA173nv0kA/e4impzuPcmCk\nRynh+SoxcwqiCkaft/UkdwGmweuNXZa1xYvj0l4gLxn9XmsjeQrt3e4lza/To5FTPr8WSAtJMcgn\nLY/RJxz9X58tfnOLmfR1fmfmymwlAaPvGxwp7iWeV5Z+5u3HKT79frj1/PuRdiE3832X3vQ4dHjw\nbxJb3s6Y3rr/AT168+EblFK/VEr10C6FF5fOkkyYudONF+MFEb5Gf/Y/V01IWbPgkAeMfmGVnUT9\nbumKbjTqTfqiF08sQNLrkac5Km2TdrtkE07aaS2EGu2nPhhvfg07jtIRLMt1axjQlzswq9F3V06Q\npxNBfAyVOq8bm6dsIi0YvWTQng2F1zzVbl6STwe3c5OKzr2gaiJbMnNffbfi+MQuW8l1fnoVqoJl\n9AMnyM1LrgR603nEermdX5lfCz9/rC/dTFDW6EtArzYq4bYHay1i5p7785C/W1Fo9JK3SHqa0e8+\nVNN3GxvfSh3v/3ofTWrAAv1dp13HnfcL66xKIZtYrxk9eg+B3nSGbbyW0UtfBWVkgX7dL7yn6MHu\nI+Bnj6fSWvN3lI5JxOiRPNDouwZY7eI4Wz+m79HeXt961ZsAmD3ULd5SaVJa6Oikm1xr9Kl+Nivd\nlK0a6Kukm/Kcm5Nupra7NmWjjZ762XVAqyAQYYJtVJIX0s3am251PxXSjZMsrTmgXwHsKjyv0q7O\noyrmAyYrOanBJU+jhx3HzqqN6sZipbZvWVsxuWO/BPr74G8hBrcS7mVp7pInrLjLibmT0ajFZ/RP\nij3HAnO+dz6j70+VwbZiCb10JCHJta6cdhVpr2cWuLRJe3rRQ5KJdORPWHvjvwSN8CcX/7CU4Huu\neSzXPL/M6FUwCRn8KMI5wBHEZhdMdVe0jTzSKhqmSqiWbjILlqo0UakKoA9F9CP+G57+bAhWyFYy\nevE2BVHBz5KBSnIRLmB+bVnwlUx7VtT5pIOdoAc4gTx1lXr1rQ5YFZnJ5ze4/hmuE9adkJ+tKkb/\nEuACyu6tYQPT5ZkHjP7td36H9/0nfFBHUy7mIfpTLS4pGqiWbqrcRPMUtp748OC5sqk+hRtoH7JJ\nRWthgfU/c+EIssmcmbvgcK+qpV3ozszwzddUBylv7y4DfczoVZQ/C6xxG5ncEbr1bjrbBYLLU1i5\nJTxf8h6H/PRyVKrnlbRGr6Wbam/q0YA+YvTG/pqfPeGbrNz0fY7+9k+KozY2kx5ZpjXhBlrkhtG3\n5lMkdx2FmH0nnvSiE0tXuVhYK8nTuWKdT9o1nVnLA3rP7IpyO6+ULuRFpz9ozi5vK73r2v4H9MO0\ndG1f37WGrwFfQ+8aGQ8BJ3ZD1pobwb3Si/HcDYFeMn9zDPirTfq/HYplLd0ws5YdWiuSXp/5tStR\nSUrS7ZK1LRv6XYS1HqPV1/zbO8NeqLsKtp1QpdGv9EYS8Yu90vwPg50V0k2qg10l/ZlCFlHiAb2V\nbtqtQsfXfstlP/X+NEg2USrX0z8G+m0AvCF6jSZNj9FL5vQZUdbrJgfeU7Bv3861250KfOyTsLPc\nr0V5dRk47EdODlJJrjaqrtqoHsVVL6LYnFvrwj77mgD+jYmdP/NSPcf8Hwz0kzu0K6uohLSXG4BS\n3PxwrTHrfOj/WcuBtWX0esI7TNZOXPobXGhGb/YE7Sr6E7qzc5ujQ38yD6QcMGsWWkKSbabKJneW\n44bEGr31o7flZUdhxSpoY+25BEVenOeDtd8WimNi57m0102eQj4hRqpw5+046jHsPBLqXAZ7UWSL\npFea5Fcb1Qe57PNXMLVtU+EqDDB3iE1De0GlFSGBlbRRqZFueimSu07TMvpz/vGk0nV+u5lfk7m9\nfhdClh4DvTU9CsXT6KEG6EXo8JOVR/HTz57Bl1jNl1lddd6e2lKB/jbgGO/7MWhWH9qFU7u5ELgQ\nOKEiFb2kPh56V5n18c5JvNmabNJWble7bEWxoyvNIPUSaWWAvrvqaGYPWas1esPo9XnWNcrdWaWw\n/ZjyC81bVdLNYewuPA7qevDbg2/OvVI3tLSXumfIKJZR2waUT7S8kAnKix/imL9KQSX9Ugwh3ZBt\nb9uLOgLD6HMBQQQh7a8PftbSjfYiqZJuLui4zz96GuyqAPqklwL/ZfIagm/B6FPPZdHT6HWH6PIK\nk0rxBV5w4d2ULS7/8F560ldr9DN3XUMlGNm5ipabFIScPG2z6wgdg8Uf8ueplsh8+aE/nWH3VZ3c\nOceuI+9jNhXx3BmnyiClmSOgtpZ+A2jvLk/UljV6Xafse7YdqZVw7MR3e9bIbco+r0vhB8/5/eKz\nDXy34foJHvzOC01wM83ob37YM0i7KugkPvHR03RHUzOn7MppW5GfanfpFUg2D15YCOvZo+thSrpQ\nMf+VTaDSnMkdcNRVRyBZl7/YBtf+pmP0Vea/000PWl2QRr3C1i+fKvdXMS7BGD96e7yOqf8ZGw6D\nxzxJWEUAACAASURBVKy/g8fxHR5LHOJ5SbZUoL8KOEVEjheRCeCZwL+Wzkq7tWuudS76oJJuBaOP\nzYJT6KpkCzRYnWg+upWedqVmnzyFtTf+iKntLR1+VlLSXi+YjIVQu85T2L2hAuir3P6QIvTq7PqU\na5/xmtJ18epU53WTOD9sc//DfwCrNukHKqSbtFXcd+aOlrf70UKQ97zVY+bOqs7GZrpX08HaPWPD\nOiKFdKOoA/rYqrT6tNvCzhP4Egc4jT73R1TtQdKNbWgPoWyDGT3A3accztz6+/PIt6xgzS23lM4p\nQgdMO39uUGw98RgWVq8tJasSf2m/tt50pjaqjBseN4dkLTMfBCgP6EvL5LXlKSBhffn4/9WAPbGr\njJ5ljV4VRAJcRxpLJu3ZtHYdgS2D+dVv4N3Rxt8qUaD0tf2p1aQLoSPD3af8PUlWvyGJu+ehevRa\nKd0ArCDJ5hGPQdvw0XlbgLQygNjKzadx12nzTO6ElVtWkvT7LKzBeNnVg454q9G3H9MKQrWkPT8Y\nm653H/jqy7yr2550o7xz6zB3nrwNa28+CriO/Um6UUr10XunfhG4HviYUupHpRPTeEY0/r0H5ZWN\nMEi68X/LUzj2iq8QLLk3P285A96+xfbObdpzO1i5+Tae/dTvMHNnqtmmcWfUL95dHIQGSGBhdbkS\nuU4mtJ1m7vaXF57CJz7+He8X22Ad+OaJ50ffToyW2Ka1MMvW43ew4k446uowAFnWdhr96f/3EKa3\nnhCkrxupkLd7zNxVRaUc0FcxehBIykCfLuxm+p4fGM2xX+l1E1sV0OsRmW454cI0j9En3jtuGelG\nYWKH+EPrAZMBIwD97g1rkAxO+spL0SutI6A3+e/NOH9uMOBW8Wwq1aPUbMLlq7syK55Jcj1a0/73\nPqOvBno9yanf68c/ro/d/DDNrCd31jH6snulvZUF0W4kmUzsEvKWFEzWgvWWM2DT2fY5vse24+P8\nKRLD6FWqRyzh6tOrkRxe+JAXVz6fz3bFe79lW0Ham8dn0FtP0hOyeesI4JjK2O9TO07gjtNdL1Os\n9k2o2DLUzzdkZgJx55FZ4URipRv//QL88kJ/8nEqkG6ytj23LhSMdgZpzU2gCdD+A/QASql/V0qd\nppQ6WSn1F9V3ySbJKlZQFr/3QNQojN5sB6hCAFIprLvR+TX7NrFLr6bTcSsmgT5iApStuLNlZIUE\n6auANdp0i88J7Dyqgj2pee7z3a+Vjlu5ojfdIoz5YrVWl1Z3lZNu8pYOX9BaaDGxeyfIq/GtYOrt\nduBrn2S6M+1NTQV5z9p9Zu6qKthhQJ9UM/o8Y2rb5oLRW9D46KfqJ5nqvW/M6pSI0RcafUvELpfP\nWzC5c4ojr9EdQRjobTFAX5E/iXenqgP6BXzpRvK0uHWVdOOz0t6MA/rEAH2S4W2IA/2Z8Dm+/ccm\nvRbF/qS3PsQdW7UZLnrdU0rPU+deWTB6k694fmViVxKOokxzuuYFFCs/k35ZUijkRBN2Y2rbLlbf\nelPxe9r9L9qz82y4/rGla737FJ4oF78Qf6WpuJAJK0h68/iaeHclfO+FeBOl7kV89t0mgWyG7kpX\nP9OufqEqgfbuARKJgmxC6zU77+PCJ1vpJlxnA8UCOgUw6Uk3kBWRPOuA3ioKdi3D/gX0I1t31QBf\n4B74C15+9NRLAb1hQLhwwmr0GTGjT3uPBcINTcBb9JIo4NfRWwxqoJ+5I2HXEZiVqHp45TP6u09V\nfPXNJiGBhbXll5QnfSZ3lsekdqJLD9l8oLcrQV15dFfqRji5HeOLDO3ZNpCz7sZQm7VD6N6KsypH\nEnnb+RtrS7n/J8tx3O39T/3cIRzyk4qfEVM9wjqS9PsmqJlm9BYMbvjVAT7ztdVsMKOX/k6sF5e9\nz5P/wGnLweij1gZr9ABZO4kWA0XnWLfEFQv40g0k7tniydg58OYMHaM3QFower8ezEShAcxr0wDa\nLdIGB4Qbrvt97wrdqSUZvOuaD4qYQGp24lpy+NolV7HDLKAOd3aCyZ0JubchuiM6u4vnTHvluQL7\nm46WqudSxOtoeise5nlZOfvaRj15nrd8MHfP7MzGRJkh6UdAP6PoTxqW3YdjvuXcoK9+kU5H1Ep6\nK1x+iphQAr/ymU8Rb2ZTnJdpogSw6/AJB/RWupkM9EPg/YDtaO900k1PGWkJ6jE3MV5/E1Ds9zw2\nWz6gr94xSNvxXwPJNaPPE/j6m/T+lgtr1gasOiu8GBaIGb22cJ9Y8Fc3KnQkRrBAf+hPtxnmLYV0\nI4F0o7jhCUFqpfRVS5jYsbJ0vNCZTRiCnUX9E/Rm514DX6Eb4W89BdbcpIFuYrfdJSoCqsDlsh7g\nlJFupre5sMphp6kr0rN+/c08/pXl43WM/sQvT/LkFz3O6Oou2qNKBgzZarNpNfoI6E3ZTey+HdC7\nCFmNU681iLJa5PGjFfeolm6+cCnMrted0y0Pu4A1N5XPAdh6guIOsxYpnwilG/Gkm2GM/p6TdUW0\ncpQDeg8QI6C3YJ63IE99SS6cKHWWFIt5uqseRLGKV9yI8b+fe0fxPrYdD7ee97ji6ondQtbOSpOx\n2URaEIzWnO6cP/ER6E/o95dkbv5IJZD2kuC5sslzqPK1vfs0ExKiFB4gnvey4anbxjXSKQS9mS7Z\npAbf1px2Irjj/u7KrA3za46nN+25VNICzjWrbHtKsalclECSu/q9+/BJ2t5kbNrtmTDMNvy0Gxlo\nDNlSxP1Ju9AvZJ46Rt82Xn+a4B2wjD52ofJt1WZA6U0Okhzm1jl/W79nv/1sHSJY6BIzem3l5zF1\n0TD67wGXYIF+avsUu47QWmuSKRPrhiJtJVHskwo3qum7E9bf+MhSPbaTlNqNvc07bsZscjKBDlnr\n67duWL3ulxro27vsC6+WHqoCYAHsOiIrfi9dE4CDvjgpaWoeSxbtQujbhus03dYT7N4WalKx44NN\ncYh0EzN6K91M7rgVuL+IH2dcQXcmHsoMlW5EOEaEl2Of79uvgB89/UYAbnkYnPZZ/xqbn5fwt7+Y\n85ivPxmb05tZ7YDeexcFo590+fr5YzUQ5Im+wE3GumfvrQyJRNGJpqDSkNFnlcU9T2veTup7+TGe\nNNqrSydw+UbozpzA0d8NZ1bztqtvNo3etFsQ1J63ISSgZebekm6rOKadE0Kg16y6/I7yVmp+D+fc\nwNu7OLAUlAZ6O2JWrV7B6CWH3nQv3AKyDe3ZCbor/DqjnTKUlB0pfJN+Tm+Froxzh7hEWwuQLtj7\nWKB3mLXh2muBlN6MXiOUjCTdtI2cZxn9gQr0FZN2//oeJ3lIvqB9UxOYPcyLw+Iz+qLhdKlm9GEh\nfvUt8A0bNFEUOmZFF12QM8A08+v6YHaYiqUblYRuYlWMft3Pq5msBXpdGVbryk6CjctSAnqvHvan\nYGK2htF7zxwPh//yHth6Ylb8Hlu4+46qiC6oj2szyx5J+PljP8r1ZjBkQwVPbz0UuIGpbf9t7ldX\ngUcB+mqvm5m7bkMzehdVMu31OeRnVwE+i06i/77ZfL0aeIf3fNCb1ox+9lBYfVvpQqXMhvRugi8H\nRITPMb96im3H3dcDek/ysF43Hi+Y3ZAHWcwmyx31wuqwftmFhFqj13kdzOjnNKvtl+eXijmgNAV6\nXH4J3P0rC8TupHkaAjRo33j3DvvBb+D2OrAafdJL8Out7miqgN5c19L1zLebHln1fC1E9YAp7jkp\n4xuv72MnMdMFs75D8qDPyNp67kFPpFvTbcuPp+/bul8gHVlNkuUsrNZSxOx6V+B6ZWzPvIMy0Lfn\n9HqL7kq9ODTt2pE9pef086Tduy3Be3jNeXtkywj0M+V7bTnTdaeSd01EyS69FW6DggAkihcYVgyf\n0femtvLJD+tv33gjfPVt+rOOnz6FA/rDgbvJkz6Im4wNenhRJUZ/9ynRMzwgj/Jgn1f/153TsW7S\nLqsBelPh/mqTYfS77aRM6Bvur5aNETRvQ99s9lDpDeI1gPMvPYU3DvSYce6VH/rSv7HjPhoA7ERW\n0psEWmy49qu197PzEcM0en/iDpxGP7FrG3ajD5uGnoC2epzLq7bwJXzhHZ/Ardx+ANFFXPHan/Kp\nD8XvTgXnEKzwtfXuiew86lDjFlxK1vOj1xd+65VgHQVyszbAkgpRrveNXVX9zj83q3WcVEaFzdGe\nM95bEaOXgNHbhRVliSBv9bj6RfCZ93rSTRvau24wZ/RL90+7aXFMyyHT3j1sRzOA0Uex2ncfmvOz\nX6t6PhvSYpL+tOKrbxOgp6WbIhqtCup53ob2vASTsT7QVzH6Pz4J4F0k/bx4zvm1rteevhtO/Epi\niJN9Ts/jRiVAysIqLR2nPeV56JQIkQgpkHDilz9upJtMbVQjbKA9ui0f0FeFRvUriwV6/bZ0pU57\nBA3IhS9tAcIOE23Br9TZxC5++NtV98pxjN6C/gJ6FUda7HQTMHopM/q//ync4rlrf+tPdE3ZciZc\n8/w5cpPJYvXhyj5wnM5nkpH0LNC7h097W4oKt+twA/Sztmf/cflhsMPhkNb5Q1ZdtkYPNnMbvo68\n7ufhJtRlM+hgwv3aHaVsGknWBiZAGb/9SvVkt5cXbe+65if893MtaGjK2uqGOyituVmvGrZ6rG95\nmmO3kSx73YT1+dsvP82LKXJBdC7sOmqOHzwnHv3EQJ9w5csAfoBPMHREUfds4kkw2YRu4Hqyzp5j\nPseTsZ4ffQz0RWyVBafRu2BhL6dsWrrRG524o5KlnPavn9EhHlptQqCPGH27z+4Nfa75XU8mmsAw\nabBMPYgAmTt3ZJVCe36SAOiHSTctO1Jy7y+KdyXC3cBZpk5MmEV2KdB10o1l9ME99P/uqjLQI4Pm\nuSaCIHf5hGNFZ34IHvXWVcaH3zJ5F8BLR+pdSXcVTFhGX0zGrhPh+dG92kCP9b+41WzKMlbZBvY5\n0At803gPCr1ilWKxt2Q/XHXqet82kHDbuXCJusqrdKpUQ+C17reA0WstTEkPjNdN3rKAZRqzKLYd\nDz/4LVtJNOj87Imw6awiVwDcfQp85v138L3f/yfANVodh3wtsBOV5qS9MqM/6ctf0EHUEn1YA30L\n3bPv5kv/y3scAbiSP12/prQ9Y96uePzi2UOrd2W1deIkth+bAeuA9xFLSJK1gEkkX/DyFZsBeu+3\nLQ+cZfOZtlFYgA9plUp2mnvoYbq17/zRZtpzQnk/zTrpJs7UTdGxqtGYivIj3H4OaqM6Ex/o8zQJ\nR5xeEfdm9K5Rds9U/ZpSERIvIqgBpvRW4PNA2a/dX3Kfm5lAl9f3UTYt3cQa/Ybrj+Oc91yswxWk\n/gRiRgz02stEh462aWjmGobzrVs/4I5f6I4PlW4coag3vTLKdv5uXqdXTMbqOq2CV1wAfQ2jlyxs\nG7PFIvCdJFnuPY8DeitD677X6p+eDK0S4MhirUd7NvH86J+N9c4J89PFtbGhIQIWa/ua0Stvkrav\nw5t6jD7ph4DkgH4GEO776U8BPwwaatKPV0z8i/4nPqO/Bz2Tn5lEtXSjCqnG1I4kpz8Fn7rMZkID\n/TfeCO/5nk3fvHUBmGdhtX5JBdCvspV4nt5Mm6S3hhjoYZZ0wbGkULqBb5UW1orOe7RuwF++Hhyv\nAOH27rpdi+zJZ7LprC52kwhRIRjrBU8TSOZA94o/idOynib2+5FARm/a3sOyoXi3MOMmmIdA35tq\nMbFLKBpX8awXSUf+E13Or+b7z7MurFW9Txnoy/H434KLwhqfb4C+JYF0IxHQ62X8TrvWeyFnxYlZ\n2+rmKTb0cYnRW6Cfp1gw5fr2qiBn81q66VOp4WuvthS9ot0+a8SA233s/rbO6wYc+Nj1AOX0rXSj\nzT2MEkrbB+r0jStwy07GDgJ68wxKSzcO6LvBZKxK8qC+2x2ywrUamygmYyNGv7tw3HsUkzuOCZ5n\nyxk6bMm6X5q0t4EDejfhVey9IdoVeGpr22P0VXNZ70VHyDwIgL5qSbOdpdfWM9KMnmABo2FWAv16\n8Hasd42sxcTuWeDJ3l0sEFm5pgd8GT3ZkaESxeTO1ZzzT0/VmnwLbIVzmwNPRv/rbIH5tTqTttfX\nwK/zmrXhoX/d5uIXTEKgWc4FYZH7U9Cea1EZA7+QlqoXiIWMvoJuGpveWnYJDa9Zyfw6b/Y/r3JP\nmECUA5z/+J/m1KTHX2yDsnQzB2T0p+0By+hdBMbrnwY2zrdm9I5J9afatHcnWPbknvVJ6PeZAD/h\ntnOtXtoS4dia56N4tnBglCvFt5UqQnlUA30o3fS5/UFu6G4nk+1+uvqcN+nP5l3kViZULWwdjT3T\njvvPrwNWuokZfdUk/RwTO+/xJ2NFPPAUJcYN9lbgV9HMvorRbwvulbWxIbGdFFZBKPLUz583tyR+\nW/LOb9kRjy3XEYA+17Knc8nV7pVHfm+7CbEQMvqpHTC3zpeoLgfegGP0YfoO6E8sInJqmy42lrd2\ny/lQBfR+hZpfCzN3Tw1ZGfsM89+2sQNZuqnCSPGZR4+Tv3S5YXGOJfov4tpnwaffD1pOsPRVeedY\nXdt3J7EacoZj9HrHeeiTTRzOlJETldh3ZBdUxJWzHuiVYfRza3Vl6BdAn6Nfru5EHvnncNYHtjC1\n3a/UcyVG35oLd/t5z5V6UnHNzZhnDxeNuYyU8xvvrwswsbsu/pCtE1P0pxzQx4xeWxqs7LSWtxZM\nrHMN9G4FppYKnKuZBfqXACfpj+IBfb6Az+izqTYTsymomqAppkO98qVbeMsCwGlouca3UaQb3/wN\nIkLpxjH6Hh/64k/45qtfwYe+6E3ET7gY6u72+p4rtnxFM3qVYutrzOin79bP2ZqHItaNkcNVJRjM\n8bTn/F40GesBfS7k7TPQndmXlKpg9NmEA/pizUZKKSpk1SRmyOjjTYBy9Pve7J1vpZtVwOlBXuvW\nC0lmt4K0J/wVKzd9mrU3reERbyszeoCdR/qTzneojUrLtyqB1be/RTpycdEBx8qDz+hzt80Ytz1o\nB19+O7gO12P03nBn64mw5pY1ZO1h7pUQj5rGaMsH9FVDSbvriza7GCV0KUxyuPU8/XlhDXz/+aDB\nezX25SUloPcLKmb0djLWaPSBV08o3cQTO1Bebu6u1feaWx9KNzo+TllTmd7qA/0srXmF3czDMXp3\n/03nmA4rtTcrVQalENq7/XgbafFLKbu1dckBvQ6tq6+VvMoXMzEx40tZMf810LvolRro+yWg3642\nql+4q1ML9JF0M2ORMJZurKVADsn5lXvZUloAVsXo40Sv8j6H0o3T6HuoVpcvv/37/PxxmA7ShTSI\n65i+3u4P6zF6D+gvURR7CWg32Op49KHN05pvR5Oxvp+lGLnEr9cx0GcE+wmb/Mdx3n2gv9oszg01\n+jCssmb0LeAa76hdBQe6Qz1Rf1WEOzYF9zWTwYVL7mU85Q/eAcCp//Zt414UXpNNKnTdexaazYOd\njNX2vOJzHg0qXN2YIW97uFSEH/kxcAX4i668F37PybBq01rywutmEObatnQAM/o6oFfBUHQC6CuF\n4v99AO40uxJ+8D/s+fbce4BD0I0y8+SdNtFkWsF8VKDRZ/gTMkUe27bjqWP09WYZ/exhRroZAvRT\n23wkmmX9DbeQLmjgm18DK+4sT5T95T3w2feAk27KNnAy1vutyt1N29nSkYcA0ya0ri6LufVVm0ak\nIOXRQloglgaMEtBP2UK3abpnUQJ5Sz+33lfWSTdudexc6Xm02ZFOnQkhClTpzXGim6PfqqUbDdb6\nOU78j48ATr7LE8/HvthlqmcWTKXFdbZz6k1rLdgu1b/6RbDj6J1cvnHAowF6F6TV0cburuG5/BZl\nVNrpKJss1zsdUTI8z4Y2zlpb+MFz9Ofjvv7p+jkEseRKp7Nrw2coSzXOd3z7cbDldEpmNw0xHZlp\n37qjXHHnLr0he1S1uyvbQKI2qo+pjcp6fPltv1vkKyYIPqPP2u6Z3A5wP1aKR2Df4ZUvhoXV/1ic\nZ6WgrPCjH7CC/KBg9DXyW8joz8HGt7/9XOcmVgTyKirRVjTQ6+FgWbopF9TUzsOAUyl53XhFkLVV\nIN1UD49rrAD6kNHPr3FA74Pw5LZ2Edcb5nToVfMgR3/7A5z8xaNxm2fo3MyvszJIUrkVGRAs3AH4\nzh9uZ/tx7y6+X/HazwI/G7Ai8FB0nPipYLOM3RvcSOHnj72aK1/yr2jpJkyoN7XgeQPpmbAy0Nvf\ntagaP4tdtJP0FnBAnxd1yI4uyl2VnbcZZGVGH0o38fWxdCPmGl+6sZ5iOt+rbtd7MhSSlQf0+YSW\n07KJTDsbqLTIk62Lt51nQzno9O44HfL2ApdfMuTRuIP2/D8X3lva3MO5DqC+jPoT5VWZKoG0W83o\nxdsycnL7Dq+ftPVCH/j/7Z15vCRFle+/cZfeWBpFZR82UVx5Awroc2kVHFw+wsxze47LjIw6742j\n80ZFwaVgnDduOKAwOvNUfIrI01FGFmVplMVB9k2k2VpAulmarRvo5fa9VXXeHxGRGRkZuVTdqlu3\n6sb38+m+VVm5REZG/vLkiRMnUh+93s+J647JuDj8stry+ozNmFQQGT2xAjyhU0x4DUO7Y97s7anp\nPBBSAfdHHKdFPJC2I/Tpg88aOgIg531T8aX1pyTr2RG89Sz6EfDRh0dhujeZH2KHmZeR1FeYseif\nQSL0Sb1MYK38YopdN75FL8kNeippMrIwiUW/oy6M9fVpi95Yba7QP7UoI/Tan663Xbp+NXmyg3iK\nyN4cD3P+KZv52t0nJUsu/uJNwCZ2ve5uf1OPxcai10K/9mXpPHKr3rqKn3/jOnNe2eu2aadH3W/A\nbmze0X7XD2FP6L0TSH30OoLK8dHnLHpfCeoIve+j/0KZ60aEK0WSbbI++mfdipnUxVqVthNZv6nM\nLHHfJjXLHr0XgKnlOkJGd8bezUVfnnKSg+n9XHTiJ7n0c+80W9Zx3ehyjmWqwJmjIbmXiu8PbdEH\nhN6z6BNDwRH6RZs22+XOw9scK6kDR/CDb2TO2QSKOTFtxqksdTe27lk9taNv0Yci/rL3fjpgxx9R\n7urD9R9M3XitpGM2I/T2iObvhiTHV2rQZg6glIlqs2XK/u0Zcyf02wZmQhOvM9ZFT/eXtejTi+C6\nbjZX+Oh9shb9c895nXNM66M3vtUkJOt/VZwdiUW/ceejuPy4VJSmt1WELPolGxalkzMbobeCIGPh\naeO8g4VJLPpPoDvVlpIVVH0jjzdL9gHAjDl/LfQPvOTRZMRxc8mMOacxlBeN8/2VZ/LT099nCwO0\n7Xmam7/tpHfNDpRKTsFa9M0pMkJlU+W2/agbSzcW/ecqXDeZkiXbP/lHzwXstIxtUgMCxpq6d18/\nKElCbgEmNz9hlk2b8MnlwBP85hPNtJ0bhbvhA49x6Qm/M1vWEfrQILhUWNLzLK6j1qJ8B602uop8\n9O7I4Wbi0knxBT6Nlslfv7/JfAtNiDO+VdfDIy9w210q9DIesOiDMRSu0C+hbQr+7PPPzqzlto3b\n/uw+YBE/OB/OPPtXZmlI6G3Bf+RlsYV8Z6yba2gEfPThp2qx0O//H2806XDdyravBVnXzaP7298X\nmWVNs7+vBY7oCn2TFY1V6U/eSNh0lveCCBcHMZ2xMn4gv/rf+qH0tdUNGLOv5pJpC3/6F29iaRKB\nthk9+Mjk22+F3h6yozWLSG+eGXSjWkJW6K0/siqUbYr0zUefmxXamUTo8xb94/tNccu77TSJ+qEy\nldGfzcwstT6Np8i6pzCRT7rBT0xlB0alN4tJkdH2e1yDQq+UO1lIzqJvVbhu/N+ybpZFm0CfZ1rW\nsbZ2WTWXtDj1Nrj+g+552PDKlk59wc7ovoyZVDCV+wpv23wdof9cYFn6oEzbdvE56jTF2d+3bp+3\n6GeW6Mbrum6UtHj3G/yABd91kwp9/pm6W+bbuf8GZ30/u8b4jNaJR58TCriYcAY9puhspwd6x3I7\nY/dietsH2bKDuy9Ntm3MSENmWH0EbNjHujLLhH5r6rpJLkNRWLO73dBa9Dtx3jdDy/3O2JT9z743\nmU3eXrjFT65GD2W3Qi/AJta9GOD7ZF0360QIDRPfgmvRL3vcuYEkOwWaUU0ThlZhARuL3mX9vhvR\nYqjDK30L5oKT3TKNMTE9BbyIHVdXCf1vioshzo2EHr2abbx3kRXwIqzQW3PoyURo9ZR6Vujd6/Zx\n4ERSoTAWfWae4w20lljlF2nI9bkjtxOhd8VNHL+sTcfgd2ylHZtZbI5s+3aVHEkk0wavAD4a2D5Z\nP9k+KybWos++qcj4NI/ur0U9OQtl+wXECNAuJELvWfT6mts6qJxARRqylrxIOG9EyXmWuG4ms2/E\nJ9/9HR57bj7q5pZ33WULm94aYoMcMsUCXAPErry1cqKhe18Dv31PdtnkZl0f2UlTbBtclLPov3UV\nXPRV8COJshb9zsxss44vrXf3ZUqfkciZwOec0Iskn2dSoU+qxeZcCjHcFr005OHc/JSgL36RRT8x\nNZOz6JuLF4lwGVqElqIr5B6zhW1ktqEWWUCupZoN5fSTmGXzpBcJvX6Fl4DQY9Ihhyx6gHteaz9Z\ni7slwu9wM+HleSPwgcJf05tnBp1jR4nQZuNOq7jndaAnqHYF/Fe5fWim0PW5A9pV9q+pRb9sBvgM\n8GzPdXOeCE+QNlRt0WcnuFhPS3dIFnYoW4t+crP7gFpbY4LloqgbG6s3TijqJr3mJ0tDvlewb3vc\nrEV/wUl2P+4cCf4gGrdMaXilnkBFkRP6MVfo3ekz6+Cvlz7Qi1036UT1et7a9PcNe+tz8C16xsyA\nPrUhY9Hno0r8a2zrP+S6MSXe6Iuy+5u+r7PRMbadbEd7QmUewvcfYueV9d2ELWc9wV6rtEz/pH8p\nFPqmtyw8L7MdBBfONgrZazH0Fr1fYclSpwKyFsvip6ZzPvrmYvu0sKGYIsJKtBg1gfcCz0QPWCzU\n3gAAIABJREFUkjnZ2du7aU2ebz5nLfrM3JpkLXoZc1tike/WZEc0rpssOl1pKLwS3JmSbCNsed9d\n7PZNx2JwsbHlfjSA5sSHjuXMc0CLimvRFzWqaWA/9IjOK3AtencSB8Q9jj3/bIPdulzHhWs2lDR6\nzCA6vd2SDe6D8wpHkNuELe8iH/0uzu/uGATbJ2LXq7rBUqG3enXD0XY7N61FsdDbNiXjbcf4eYqw\n0LdJsyLWvfn9trE0Sf4nhRZ96h9pL7KDi/bO7M8Pr9Ttei+mdnhFKpjSJJPFEUjrxA8v3MpDB8Af\nXpl/e53YWvz2Yn9zcv1LQzagc8hsz7YPPcDyNfcEtswYTyIIz77wB/Zr2o6T09Tx/tnOWff+tse3\n7f8U3Pw+mmbisg7PH5A5ICMh9GGD2LWgs1NQLX1sOrHo7bbadQNaQK0/HmNF2srZX4RNIiQhTiKc\nwfjMqebrFG4cvSuICt9HH7oIPvpmrrbo852HYzM2M9oWZ31IG+WKwPGClp00ZKl3DN/CsPu2nYbW\noi+yFJve5y3pRBTbuG9L7md3cJr966v6+uCkKEmpl0F7XB978VN2f48CR7tCJcLXgZ+SjXOvI/Su\n1e8PmKov9FbcdGqPLWiBs/fTFcZVWWzRy1jbhuCakZrTgQeOmLZcPG5CcxZpR6y/3jJOu3wjsNF7\nUIYxnbEi3JtZnp90e1Iasla++PiDqUXfbgEXQC7thEtqYGzZEb57+Q2BdapHj+YHxD0IbM+SJzbw\ntHvuC2yXN552u9bqyVjaP5YEM5gBihmJdB8W7pSSmOt0qXeEmUToi42bkMYMp+sGyHeQaJRTAba3\nVovessfzg1ne+FHbuWrT14Zee4qworGFrOvGsbbIum6yQh+q/PtT67pU6PPhlYBMb2N73H1L2DbK\nB9zVS8qRknXduNgf7MCeKou+CZznfE47Y91MgGd/5zYu/Mp3zDd7/onrRsTJY6N5smDUKvz+sP/K\nBSelrpuxli3bIyJM+a4HachbAbfzp0uLXl1iyxsuWEIq9PaBqi21bwBfsL9JQ7Zy3V9D/gGedkq2\nJ9s89hx339OORZ+NQXfLGmZaGvJkwXrLkDGT2jepP99iXpt8ykfd6DJc9tkbgbc6y1PlSnP0TElD\nRBqyxlnPt+gnIePHDmlQWaCAGVyWi6TZSOq6zb3xipR2Zi9mzCTnS/u4/NxCkBV6OyiwPC4+GTSX\nPSVnjtwRc924Fr0kDW3MsabOAPaVhtgKNEIfjHG2Fn0odrUIK0426mYMaGX9xMpz3YyXPkhE2D35\noidc9huTznWvj3VqYeRRsdDnozWy5Qj48hOLpMiit0K/2FsO2aiNFmDDJLW7yL5+bt0uPc/mshZX\nftzOLu6fh62/W5z9TjH1NPj5KV/MFf30lbezdXnqukn3Y/KWBy1v9xp1Z9F/domdVL6+0CdzI+gH\nvAibyb+2+m84JEHu7fE2l/wDpFFHM854kZBlVyb0RfUBsAwZ34ruqLTL/Hvlm9zwlwfpXxb7Qq8b\n7cZdt0hDfuosT2+UpLNdbSCPPSn7PRQp5VM9ejSdyMPijrK2B1uKO5drHnvsJaiWadOeRX/3YdBc\nbKcfc98K3MixIlLXzZL1v/R+sxfDvU9HwKLPYMLHRDlzgTKdyXeSCL02kswyV8T8nB1VN6mCpAPQ\nrus1ePFcN6re/q/70F78+jhbLhfbOTUGXMbmZ4SiaQ4gbYwtU8ZQB1zIos/fWMWuG7udO2DM398v\nnM9NUgvdTB9nLfrt/MbpPkTcfYYart7ntR8uHhSmE2u522c7QMP1AuGom/WkQm87Y/0RiKGHaIhU\n6Cf87hgA7iQ7sK68M/apXXGijt7GPhe/Q3/MRN24xy7Cbcd2mxegI6zSZFwFFr00pM3vX6/L3coJ\nvVXx4mgeGz4rQaG3WGH2hT70ql8t9O58vJqNzu/WnTIFnFmjTEsYa/oWvf5++1Hwj1N2rI1v0U+I\nUDbmRZLooP9+1LHoqSwtVmjOcpaNgkXvIDZOOCP0fkO2lg/kX8XKYleLcM81XKEq1xlbT+jP+9e1\ntCfODKzj+ujbXHDSjSZKIz1EQ35LeuO521+NP42gplzoZdzemNbqfq+3RhNdbyGLvu2tlxWq1qR2\nPTSXui4q+3aESE7Yi4U+LFxWBMNCP7nltkCZ/fEFLXQnsj3vtVRY9M5bXdlgKbu+LkveZ40ID4jw\nLGeRdd2kCa8S100mogsRfs+b/kaPvOzcdZO7hiKsQuc4Pxo732zqisy//W5dbrKu5lw3tu/Mb9up\nLyJ5U1WhIAJfkH23TEiDynSpjkXvTsNXkB0tU7YljM/YqJus60avYq9jRuglnNDP5RKai37slNs9\nrwlnOd7nobbov+B8nuHCEx/n8WePl8T2hix6WwGhnBBVrpt8elr/mHkffS2hF6ElwrvIpzHO+ujv\nOWwDV+VD+6UhvmWJNORQaWTS8Yaszvy8kpd9boLzvwbpw9DWy7gpq72Rs30i2gLMnK+TCtcM+5+0\nHaXu62qoXvwGW1fozcM9cd2kLj6AP3/TBwPb5lw3IqwGrjHLVgMvIhVp16Lv1HJKhX5ys/tWUxR6\nOwXsDPxfZ5mNugkd255/qH2XlbXoYb0r8NpkGsJt1tk+n/y9sumZ+nzyFr0VenfZHYAz/sGcvh9r\nn/kx66N3yGpQPmOsjxH6nI/eNUqOxoY9w7eArxbsyx57wuRVwolacztj7VPdfbqXTrgsghLhWpEx\n85ZmxtKkjHt/bdndvz1jzoRehOPSb2qGKz+2mfbkhNNJ4V9g10cfct3421RVzp3YKdIKKzQ3Mrau\nRWV5nve9iRbUp3nHCkUaQPn1CLlurNB/O1ly+1Fw9UcgH9/rvg630Q31dMAO+b6X4gen3rY1CZIT\n+jZwlVdW/0Hq1uNUYBmZZXqWI3f7Me97kevG9dHb874KHXI7Rt6iL+w0KcAKfYsd77ijamW0734d\nmRBeYzzYeXzz+4fmEvvw7caiz9erjPt9J3mh37KjXqYFNCT07n5fBPxZvhjBCa21wKcuxXIffda4\nClHkurEi3JaGrJeG3Akgwi9F+HjBvtJ9JDPTJYvcuZD9NvVD9L3TCRNkHxRlQj/UFr2LHZ4/7syT\nGXbdhLPtdey6kYaslYYs99bNbrP08bXYya7B73GvHJmIFk0333YLncN+F7KdRDcXbF9kGX4EsPP0\nuXUxDSANCQ2g8i1615KyFv1VZN+Oijq3zQQRk5gIDvcVvSXCf5JaUG4ZO7XojY9+0r8+tp2G3sSK\nOmNdK8yN8HAt+jcEj1+MFfoxxqfr3Dsl4ZXjxec/syz0MOy2MxYnva7/RpyyYe/UPZeudxepUZKI\nlAgzIpl9fICT/gBnf+fykjIW+eh9oa/TTwLticy9Io3ABDjVpMdWNIFjGJ+25+W6bkDnsr8PQIQ/\nF+EbHR7LF3preI2WRZ9BVNqZmk62EHbdtIMWfTeuG5e8YLQmZtjnl4ez5+WXpuXM9AlVVr4IF5KO\nWt2V/A1oy18k6EGhEeEUs2+37FB+zr7Ihiz6aW+9ojckve2yR17I+PSJZBtsy5TxzsCxu/PRt3IW\n/bj3vROLvumsY/tL7H5CfSBltLEdumPNLoVeWYt+WX51U/at24XqqHuhbyc9x03vb/7YMubG7O8P\n/JX5HOx9NghP/BGsenvVgxLyrptdMt+keEYcg2kLhZZ/J5rm3odNEb7ivHm4rhslwo8KZvWqw4XA\nKrIBEr5F7waJzC+LXin1FaXUbUqpm5VSZymlQtnzQtg4+CYffv6eZlnm5FK/debVyXfddNIZ65J/\nhV3z8mvk1FUPm6kMTYlqJ7sK7btN9oZyt6/KHFmXOkJv6ygk9Fud332hd/c9AyD3v/RWOaH5eQJC\n7+FbJuuc32oI/aJuLXo36saWcTPkLHq7/fvoDMH6Wsdnygb1WAIjY+05Lt45sL4R+u17a9HLuF/n\nxUKvgxFsh27biXcvF3pwY+NdlLsOWYv+SvRgOGdPlRZ9qBPTpc518cvm7i/NxQMwPnU65SGalUhD\njpCGPEb2vnlQKZ5P2rZdoZ93Fv1FwAtE5AC0D/zYmttNY4U+dQOUNeQ6Fn03Qu9Y9IvzHXTtcbch\nlDX00L7bXpncz0X7qmMR1R07UCX0S3BnRUrF3+KW1w+FXF+wnn9s+/e9pKMly3z09uFuQ+WynbH1\nLHq7jr2pNlIs9AWTJBTSRtdjm7GmKyi/LVg/L/TLHtUDyNrjuwTW1+e1ecfZ+OjzOf5tZ2waURLK\ny2Is+vGiUbh127+Pb9Q8Rnbk96u81asMqkwocoBOhN7VP/9e0uf7iV0+LBKsr27wQ7W2Iy3vOuar\nRS8iK0WSiJGrwRlAVMbGnU9EP9knSS96mVDX6YwNhaQVkRf6Pa642tuP77p5K3BQjX27D6EiS+uT\n+Ol569Ot0P8neni6/W0RWf+1oAc22Rwydt+fIzv61IbuHWM+lwh24tZ5UsTmBKph0SM2bNR33YR8\n/lWum01khf6zzn79G7iOj16/NajWBCTRFfkMnBpf6Ldjzyt0iKiS9wOHB/YPT+22JfM9u48Q7m+B\nhGAmvHKsZftRyiz6IqHv9KFo8YX+5cCzAUSYNrHuKfeuWAmcVrCvNaaj1X1Y+3Rr0ftjC3wjoxfY\n9mY7rZ8DHAl8CPhj5rFF7/J+sgNuivgnTlpjR9hNkjagAneGJP+RF3o3aqKTp2HeBbBocz55lGPR\ni7BGJOmYugf4x4J9F1n0SbmkIRuC6Xk7p7bQi3CvSK7zMbHSzND1JvDvZtF+ZrvPi3A/edYEllnK\nroV9gwu1PX19VU7o/f3W6YwNCb2lzFddhrXohcVP1bHyfKFPyi3NxbeJcHFg/7B+7yrXjX8+bn2E\nRmqGRmz7GKEf67VFn0Ea8og0pHiQ0Rm/+Lk05OjALw8Beznfi8SwE00rs+gt3T7gQlyFbo/Wxf1q\n8/d+E51VZIDMmrIRaAAopVaiY4F9jhORc806nwamReSHBfs4HtCntSMH4Db4hog6QUF5qFuR6+Zl\nzjqdCL1foV8nHUHn+FMLjYNzRfhsxb59H32LaouxU9dNJ29BLq81f7eS5uzQKzfkQXM9qvpbyhpj\nyJeu9y+0lH587hDYzkSkjPlCb8WrtkUvgpjjhASqW8sptejHWnUsx9+bv+5bXnX/zKP7z8ZH363Q\n28FqRWkkZiv0dfulitZrep2hfbXoAcybQ88Q4VpgW6WS62pDDvW5nM0BbA9cy/vV8aosVXnHVAq9\niPivlxmUUn+BzpH+uqJ1ROR4AHWCagDCWcGLtD6wDJu5wO7KXUh2pFonE+va428FkIa4KW/TC16W\nZbF632Wum9lQ13UTSppkcfP+hKI/oHoGqq6E3mHHwDK93fi09cPa8/PftupY9BbXPZWhwO1URuqj\nr564ZXuyw/Lt9mXiYSYOXxJ6SNcVetd1cwjapboVuIzUgsyLmiDqVZ8X7nuFO5+Ay2yE/nBU65yS\n319M2s9RZJH7Ze670M8BNk+1Ppcj0SOjX8OnpSHTSqlGrw5UKfRlKKWOQM9N+moRqfOK8xBwO95F\nKn1y6uHVvr/56ebvbF03oYx2jutmoqhcdW66boS+U+EpE1t7bUP7fJKw0LpUCX1ZNsA61+LpgWW6\nrNvd361FH8p1ExKo0HX9a8pm7kqPW0voRYKjh6uE/hH022VocFxRn4bvrz6OdBCcDXmdJnu9wqL2\n68+0gX3Rk9P4VEbdFP4oXKxOKP39FuUMXQqssoliN55Przpj54KsRW+MLpO2uqfM1kd/CvqptFIp\ndaNSqmoQwa7SkE/QyStzVugtdsDRNs4y/2FQhh+ZkTli8qnYdVN2jDrhlbOhrusmky/bw95W2wd+\n87cvoqwxVln0ZwA/8xeaV/N95da3zwDXOvvxLeMyi9797evoiCG/DnJiKw35N2lUGiudWPQutYRe\nGjJj3i5rWfTSkDG0KLuu0DulIXbUpjUEtlJH6PUx3gns4y1/GDLjJHyKzukt5l/ZOkX7Oh6423y+\nknyZi+6nruPoO9iuV1j9sudSdj/OillZ9CKyX0frm+RRjv+0Gp3POWPhSEPuNn5k1+3QjdCX+W97\n7brpVU96XdeNvbahRFP2ZrgGOLRg+164btaGfhTh3UUbiuibWxpyMIA6nrcBv/P2W+ajT76L6Agi\nxydq6db36gr949Rra0kZRZAyy9Y7TrKdwW0/vkAVtS13dHDV2Ae7fBKvv0yEncoKW4QI53axmTLb\nnqAUr0c/dELnWHQOnfi2B23RZ1032r32sX4caFApEOqjLfqijkVX6KtGnbqUWfRVrpu/BU4NLPe3\nD7luetEZ6xJqnJ8yfyfQaQlC7gibsvlRiutrNkJvwwNvr9hHJSL8RCTZT1Uc/URFuSy9EPpDyaZ9\nKMMtY51jl7lu/IylLYpFyo0uyk60HqYTY6lfhDI8Nqnvo3+yYHmIQVv0GaGXhjwlDfnnfhxoVhb9\nnNCebKFD+Z5BXphDHYl1bqQyH316wZ9292V4eWlESkXe3T4k9L8B3lyjfGWUWvQifEkpngXc5PmJ\nXdybqai+ql5n/Jmj3DJsKtnvbAi5hFwxnKAsoqQ3x58A2tKQ0NwCRXQr9CHXzUNkI5YKLXrnzVlR\nT+jTbI6Dw60fV+jrWvRlE4GUHasn4aMdYiPb+nGvZJj/Fv2+F16KDgc8CD1JteW/Ae9yvndi0Vsx\nKHfdrPiHs0Q4qm5Rve1z4ZUi/ItIadhilSDdC5mY9qD1KsLHSkQeqq97pVVsXCylqVr7QB2Lvo7Q\ndxsbncbRd0anQh+yrO0yP1NomevGsi31hH6J93cQuG3TvlV2YtF3IvTusex2c/k2Y330fX+wzn+L\nfnJqJxE2gGRS+0pDzvLW7OQC2YtabtF351d347h7Gl4pwt7eok4atUuV0Depdt3YGXzmkm7CKyHf\nNq4DXtrF8dM4+s5wj19bRL3wT3vM/wn8vbO8RXU7XU42kqbKjz0fhd4/x7+C3P3wKdJ5juvgPnQ3\nBpb1m1AGy74eaL7yEQiOypwt1o8Xslxn22la1IHW8/wV6NQEl3axXZXr5u/oYwTALKjqjG1m5wAO\nrgPlo3qrjr+I2V3LjgIYvGMDbBXJ5LOpY9FvD/wz8E3jVqtikELvtsdzIcnnnw3JFm4kmxIcaciX\nOjxWyKIfBAvbopeGnNLJ6uZv5RNZGtIyUTuhEDm3QYVGb1bhzlPZaXhlR6+NJt62KLd9GaUWvTTk\n/3Sxz7mgyqKvErzV6Dwr3Y46dDtjO93Ocj3VYxhC2HMLpT8oO+8Z4C6TDbbqvG9Gz2HcD6HveGSs\niM6npBTfpT+dpSGLfhD0XYcH6aP/UI/314mP3hLyMbsN6pkdF0KndrV9B/226Ltl/vfNhKm06Cu2\ntyMNu/XDdiv0achnQ46Xhry8i2MXpQKusui3AT5T8xi/Nn87iVyB3vm1ryc8c9MMfUj0RVYrBmnR\n5zOO9phB3vDfojuLuYhuGlvIcrENah2YIcndY/f15ZrD7eeqI2hYhb6qM7ZIDOw69sHe7fnPWuj7\nQKnQm9mg6pbX9lm9ddal6o7TRIKTwYQ6Y3uB2w5C403mgpPpzv3aEQNz3RjhC80x2fUuzd+6Fv3F\nQGjqM9ugXmAmC5gN9ga8ZZb76TWq4PN8pyq8skjQbGeXfbAPo9AX7aNOZ2xd9NSUUpR3albMpp3V\n6YfoBlumN0hDBuW6uaeLnEsdM6999B3SqX+7KFlbL3NCd5sOt9/UiaOfj1RZ9EUCbKM37LleSL2U\n2j4z6IfFfHLD9VIELwL+tEf76obCzJX00aKXhlxQtWIfmZO2NEpCb5mtcHWSHK3uvuqGIUbXTTnd\nWvS20/3bwLVm4pQ3dXH8B9CZWju9Tr3wwRad2w+Am3qwf0S4BD1P7HyjXxb9j+mt+7gbotB3SDed\nsSHyUwp2j93Hgz3YVy8ZSqGXhrRNtJQrtHVT+bZN3P/VJetUcR+wJx1GaIjwW6V43iyOiwhPKMWf\n5Jbnx5PMVw6he73pi0UvDTkHKEufPBdEoe+QXlnDneS1r7uvB3qwr17iCv116EyRw0pdoe/FrD0P\nAruSH51aiZOvpw6PoOPe/X1c1Olx5wvSkDpvHWWum35Y9CGK5v/tF1Hou2S2Fr1Ni9BLi7546rTB\nkAi9NGQNcPAAyzJb6rhuoDytcl1s27ixdK1ZIkKTPmUx7BP97ufpl48+xKnAv8zRsSAKfcf0KilR\nz4RehLZSLOrhLPK9YihdN4ZjyMY8132T64XQW7FZ14N9ReozZxa9GVk9l/luotB3iM193SuLvicX\noEORj52xFUhDvuItcq9T2bXvxcPWnYIxktLvdjuXFv1cMydCP7Q3vI80xD7xZyv0M2Z/g8zJHalP\n3cyQvRDnGe9vRPMfwIf7uP+59NHPNdGiHxCDvInn6uESmoNzWPGzVxbRS4s+Cr2Dzi7bE7/2fOiM\nnWui0A+IUX1FdHkBg51FqJfUseg/hA6NnC3RddNf/NmzLOeh540dRaLQD4iRt9akIX8YdBl6SOUD\nS4ReZeO0RsDIt5EBsCfFcwzfTTpR+KgRhb4L2vRuwNQgGBUrey6p2xnbC6Lrpk+I9OSNaxiJQt8F\nz6dHnbGRoWHO8s6I0DJzsI6qvzgy90Sh7xRpyB092M2gLPob6S7R1kKnk4lHekV884r0iij0A2Ig\nQi8NOXAQxx0B3BtlIXSkR0aLGEc/IOZbuoJIOa2Cz/0kWvSRXhEt+kEgDXmA4crRvtCJrpvIMBMt\n+kikBoMQ+kikV0Shj0RqEH30kWFmOIReKfUxpVRbKfX0XhQoEumQ6LqJDDPzX+iVUnsAhwOjNNIy\nMlwMojM2EukV81/o0bPgHNOLgkQiXTLXrpu3A2fPwXEiC4P5HXWjlDoSWCsiv1UqBqlEBoZ7o1ze\n74OJ8O/9PkZkQTF4oVdKrQR2Dvz0aeBY4PXu6j0sVyRSF3ujvE4a8quBliQSqcdP0G7GdzAfhF5E\nDg8tV0q9ENgbuNlY87sD1yulDhaRhwPrH+98vVRELu22wJGIR9v7G4nMa0R4m1IswxN6pdQKYEU/\njtmV60ZEfgfsZL8rpe4BDhKRxwvWP76r0kUi1dgO2Cj0kWHCJk9M2q0xgC+135VSjV4drFdx9DHc\nLDIookUfGUasgTIn2tmTFAgisk8v9hOJdEEU+sjQIUJ7LmNY4sjYyLAThT4yzMyJBkehjww70Ucf\nGWai0EciNYgWfWSYiUIfidTACnwMCIgMI1HoI5EaRIs+MsxEoY9EahCFPjLMjM/FQaLQR4ad2Bkb\nGWaiRR+J1CBa9JFhZk6i6eOcsZFhJwp9ZFh5PXDJXBwoCn1k2IlCHxlKRFg5V8eKrpvIsBOFPhKp\nIAp9ZNiJnbGRSAVR6CPDTrToI5EKotBHhp0o9JFIBVHoI8NOFPpIpIIo9JFhJwp9JFJBFPrIsBOT\nmkUiFUShjww10pAo8JFIBVHoI5FIZMSJQh+JRCIjThT6yKgwh1MtRyLDRRT6SCQSGXGi0EdGhZlB\nFyASma9EoY+MAntKQx4edCEikfmKkj5HpymlRESi/zQSiUQ6oJfaGS36SCQSGXGi0EcikciIE4U+\nEolERpxZCb1S6m+VUrcppX6nlPpSrwoViUQikd7RtdArpV4DvAV4sYi8EDixZ6UaUZRSKwZdhvlC\nrIuUWBcpsS76w2ws+v8BfEFEZgBE5JHeFGmkWTHoAswjVgy6APOIFYMuwDxixaALMIrMRuj3A16l\nlLpKKXWpUuolvSpUJBKJRHrHRNmPSqmVwM6Bnz5ttn2aiByqlHop8GNgn94XMRKJRCKzoesBU0qp\n84Evishl5vtq4BARecxbL+YLj0QikS7o1YCpUou+gp8BrwUuU0o9B1jkizz0rqCRSCQS6Y7ZCP1p\nwGlKqVuAaeC9vSlSJBKJRHpJ33PdRCKRSGSw9HVkrFLqCKXU7Uqpu5RSn+znsQaNUmoPpdQlSqlb\nzQCyj5jlT1dKrVRK3amUukgptYOzzbGmbm5XSr1+cKXvD0qpcaXUjUqpc833BVkXSqkdlFI/MYML\nVymlDlnAdXGsuUduUUr9UCm1eKHUhVLqNKXUOuMFscs6Pnel1EGm/u5SSn2t1sFFpC//gHFgNbAX\nMAncBDyvX8cb9D90dNJ/MZ+3Be4Angd8GTjGLP8kugMb4PmmTiZNHa0GxgZ9Hj2uk78HzgDOMd8X\nZF0A3wPebz5PAMsXYl2Y87kbWGy+/wh430KpC+CVwB8DtzjLOjl364G5BjjYfP4FcETVsftp0R8M\nrBaRe0UPqvp/wJF9PN5AEZGHROQm83kjcBuwG3r08PfMat8DjjKfjwTOFJEZEbkXfSEPntNC9xGl\n1O7AG4Fvk07zt+DqQim1HHiliJwGICJNEXmCBVgXwJPoCWKWKaUmgGXAAyyQuhCRXwPrvcWdnPsh\nSqldgO1E5Bqz3vedbQrpp9DvBqxxvq81y0YepdRe6Cf31cBOIrLO/LQO2Ml83hVdJ5ZRq5+TgE8A\nbWfZQqyLvYFHlFLfVUrdoJT6llJqGxZgXYjI48BXgfvQAr9BRFayAOvCodNz95ffT4066afQL8he\nXqXUtsBPgY+KyFPub6LftcrqZSTqTCn1ZuBhEbmRgkm7F0pdoF01BwLfEJEDgU3Ap9wVFkpdKKX2\nBf4O7YrYFdhWKfVud52FUhchapx71/RT6O8H9nC+70H2STRyKKUm0SJ/uoj8zCxep5Ta2fy+C2Cn\nvPPrZ3ezbBR4OfAWpdQ9wJnAa5VSp7Mw62ItsFZErjXff4IW/ocWYF28BPiNiDwmIk3gLOBlLMy6\nsHRyT6w1y3f3llfWST+F/jpgP6XUXkqpRcA7gHP6eLyBopRSwHeAVSJysvPTOegOJ8zfnznL36mU\nWqSU2hudO+gaRgAROU5E9hCRvYF3Ar8SkfewMOviIWCNGVQIcBhwK3AuC6wugNuBQ5USpJPuAAAA\n0UlEQVRSS839chiwioVZF5aO7gnTnp40kVsKeI+zTTF97mV+Azr6ZDVw7KB7vft8rq9A+6NvAm40\n/44Ang5cDNwJXATs4GxznKmb24E/GfQ59KleXk0adbMg6wI4ALgWuBltxS5fwHVxDPpBdwu683Fy\nodQF+u32AfQA0zXAX3Zz7sBBpv5WA1+vc+w4YCoSiURGnDiVYCQSiYw4UegjkUhkxIlCH4lEIiNO\nFPpIJBIZcaLQRyKRyIgThT4SiURGnCj0kUgkMuJEoY9EIpER5/8DgJtaoi8xohcAAAAASUVORK5C\nYII=\n", "text/plain": "<matplotlib.figure.Figure at 0x7f667cb0d610>"}, "metadata": {}}], "metadata": {"collapsed": false, "trusted": true}}, {"source": "### Use Kalman Filter to predict the unobserved states", "cell_type": "markdown", "metadata": {}}, {"execution_count": 10, "cell_type": "code", "source": "k = kf.KalmanFilter(initial_state_mean=x[0], n_dim_state=n_lat, n_dim_obs=n_obs)", "outputs": [], "metadata": {"collapsed": true, "trusted": true}}, {"execution_count": 11, "cell_type": "code", "source": "k1 = k.em(z)", "outputs": [], "metadata": {"scrolled": true, "collapsed": false, "trusted": true}}, {"execution_count": 12, "cell_type": "code", "source": "k1.smooth(z)[0]", "outputs": [{"execution_count": 12, "output_type": "execute_result", "data": {"text/plain": "array([[ 1.92111207, 2. ],\n [ 1.50039026, 2. ],\n [ 1.29593564, 2. ],\n ..., \n [ 1.66067946, 2. ],\n [ 1.82033731, 2. ],\n [ 1.85784542, 2. ]])"}, "metadata": {}}], "metadata": {"collapsed": false, "trusted": true}}, {"execution_count": 13, "cell_type": "code", "source": "x", "outputs": [{"execution_count": 13, "output_type": "execute_result", "data": {"text/plain": "array([[ 1. , 2. ],\n [ 0.78215862, 0.51617268],\n [ 0.5260967 , 0.50095401],\n ..., \n [ 1.0281408 , 1.42488997],\n [ 0.78545079, 1.95924174],\n [ 1.59494923, -0.39612278]])"}, "metadata": {}}], "metadata": {"collapsed": false, "trusted": true}}, {"source": "### Note the imputed value of latent variable is very different from actual value of $x$.", "cell_type": "markdown", "metadata": {}}], "nbformat": 4, "metadata": {"kernelspec": {"display_name": "Python 2", "name": "python2", "language": "python"}, "language_info": {"mimetype": "text/x-python", "nbconvert_exporter": "python", "version": "2.7.9", "name": "python", "file_extension": ".py", "pygments_lexer": "ipython2", "codemirror_mode": {"version": 2, "name": "ipython"}}}}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment