{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [], "source": [ "from urllib2 import urlopen\n", "stxt=urlopen('http://www.gutenberg.org/cache/epub/100/pg100.txt').read()" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "904061" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "#rough total numer of words\n", "len(stxt.split())" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**starts with:**
\n", "1609\n", "CYMBELINE\n", "by William Shakespeare\n", "\n", "**ends with:**
\n", "End of the Project Gutenberg EBook of The Complete Works of William\n", "Shakespeare, by William Shakespeare" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "(7598, 5568865)" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "#find beginning and end\n", "stxt.index('1609'),stxt.index('End of the Project')" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "\"1609\\r\\n\\r\\nTHE SONNETS\\r\\n\\r\\nby William Shakespeare\\r\\n\\r\\n\\r\\n\\r\\n 1\\r\\n From fairest creatures we desire increase,\\r\\n That thereby beauty's rose might never die,\\r\\n But as the riper should by t\"" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "#check them\n", "stxt[7598:7598+200]" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "'End of the Project Gutenberg EBook of The Complete Works of William\\r\\nShakespeare, by William Shakespeare\\r\\n\\r\\n*** END OF THIS PROJECT GUTENBERG EBOOK COMPLETE WORKS--WILLIAM SHAKESPEARE ***\\r\\n\\r\\n***** Thi'" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "stxt[5568865:5568865+200]" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "903810" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "import re\n", "swords=re.findall(\"[a-z']+\",stxt[7598:5568865].lower())\n", "len(swords)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[('the', 27554),\n", " ('and', 26725),\n", " ('i', 20296),\n", " ('to', 19637),\n", " ('of', 18127),\n", " ('a', 14545),\n", " ('you', 13611),\n", " ('my', 12474),\n", " ('that', 11122),\n", " ('in', 10989),\n", " ('is', 9567),\n", " ('not', 8723),\n", " ('for', 8219),\n", " ('with', 7992),\n", " ('me', 7772),\n", " ('it', 7685),\n", " ('be', 7074),\n", " ('your', 6866),\n", " ('his', 6853),\n", " ('this', 6805)]" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from collections import Counter\n", "Counter(swords).most_common(20)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "collapsed": false }, "outputs": [], "source": [ "wlist,nlist=zip(*Counter(swords).most_common())" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "('the', 'and', 'i', 'to', 'of', 'a', 'you', 'my', 'that', 'in', 'is', 'not', 'for', 'with', 'me', 'it', 'be', 'your', 'this', 'his')\n" ] } ], "source": [ "print wlist[:20]" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "brown_top10: ('the', 'of', 'and', 'to', 'a', 'in', 'that', 'is', 'was', 'he', 'for', 'it', 'with', 'as', 'his', 'on', 'be', 'at', 'by', 'i')\n", "\n", "24:00 min point here:\n", "http://www.radiolab.org/story/91728-words-that-change-the-world/\n", "\n", "tiny subset of 'un-' words potentially first coined by Shakespeare:
\n", "unnerved, unaware, uncomfortable, unearthly, unhand, undress, uneducated, ungoverned, unmitigated, unwillingness, unpublished, unsolicited, unswayed, unclogged, unappeased, unchanging, unreal, ...\n", "\n", "or things like:
\n", "madcap, ladybird, eye-drops, Eyesore, eyeball, fainthearted\n", "\n", "or all these expressions first documented use is Shakespeare:
\n", "So truth will out. What's done is done.\n", "Crack of doom. Dead as a doornail. A dish fit for the gods. \n", "A dog will have his day. , fool's paradise, forever and a day,\n", "foregone conclusion, the game is afoot, the game is up. \n", "Greek to me, I’m in a pickle, in my heart of hearts, in my mind’s eye, kill with kindness. knock, knock, who's there?\n", "sorry sight. love is blind, what the Dickens, all’s well that ends well.\n", "Something wicked this way comes. " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Shakespeare Zipf" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "(-1.030628927624901,\n", " 11.770585933080547,\n", " -1.4677534901237208,\n", " 14.753251384499743)" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from scipy.stats import linregress\n", "import numpy as np\n", "\n", "xz=np.arange(len(wlist))+1\n", "slope1,intercept1=linregress(np.log(xz)[:1000],np.log(nlist)[:1000])[:2]\n", "slope2,intercept2=linregress(np.log(xz),np.log(nlist))[:2]\n", "slope1,intercept1,slope2,intercept2" ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZAAAAFoCAYAAAB5fq24AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAMTQAADE0B0s6tTgAAIABJREFUeJzt3Xd4FOX2wPHvm0ACiYQiXboUQXqTjogQEQJeQJogxXtt\nWMAKQa4FBbH8BIUrqDQFRIpIEQmKiIAIKCJFijQpgtISWkJIcn5/TBITCEl22d3Zcj7PMw/Zmdl3\nz0DIyduNiKCUUko5KsjuAJRSSvkmTSBKKaWcoglEKaWUUzSBKKWUcoomEKWUUk7RBKKUUsopmkCU\nUko5RROIUkopp2gCUUop5RRNIEoppZyiCUQppZRTNIEopZRyiiYQpZRSTtEEopRSyimaQJRSSjlF\nE4hSSimnaAJRSinlFE0gSimlnKIJRCmllFM0gSillHKKJhCllFJO0QSilFLKKZpAlFJKOUUTiFJK\nKadoAlFKKeUUTSBKKaWcoglEKaWUU7w+gRhjyhtj/jTGfGuMmW53PEoppSx57A4gl5aKyIN2B6GU\nUuofXl8DSXWXMWa1MaaP3YEopZSyGBGxO4ZsGWPyAsGpL2OAe0TkjI0hKaWUwgdqICJyWUQSRCQB\nWAPcbHdMSimlPJxAjDHjjTEHjDEpxpjaV1yrbIxZZ4zZbYzZYIypnno+PPVPAzQEDnsyZqWUUlnz\ndA1kHtAcOJjFtcnAJBGpBrwBzEg938wYswlYC8SIyF+eCFQppVT2bOkDMcYcALqIyNbU18WA34Ei\nIpKSeu4Y0FxE9ueiPAOUBs65L2qllPK4AsCf4qWd1d4yjLcscCwteaQ6BJQDckwgWMnjiDsCU0op\nm5UBjtodRFa8vhM9l84BHD58mLi4uOs+Bg8e7JJ7r3Utq/NXnsvu9bW+1ufW59bn9p/nPnz4cKaf\nb97IW2ogh4FSxpigDLWQcli1kFyLiIggIiLimtf/vvA32/7aRttKbbMtJyoqKttycnvvta5ldf7K\nc9m9zvh1SEhIrmPNiT739d2rz53zeX3unF87Ep/tRMTjB3AAqH3FuW+B/qlfdwc2OlBeBCCDBw+W\n5cuXy7Us27NMwl4Lk0+3fXrNe3zN0KFD7Q7BFvrcgSXQnnv58uUyePBgAQSIEBt+Tufm8HTimIRV\n20gEjgF7MlyrCvwA7AY2Arc6UG4EIHFxcTn/w/y+XG4YfYP8b+P/crzXF2SXMP2ZPndgCcTnjouL\n8/oE4vUz0XPDGBMBxA0ePJioqCgiIyOzvf+Hwz/QcXZHnmn6DNEto7EGcfmwTz6BggWhc2e7I1FK\nuUBMTAxLlixh4sSJAAVF5KzdMWXFrxJIXFxcrtsOt/61lfaftKdPrT683f5t304ic+fCQw/Bzz9D\npUp2R6MckJCQQGJiot1hKBuFhISQL1++q86fPXuWggULghcnEG/pRPe42iVqs27QOhbuWujbyQOg\nRw9Ys8b6c906CA21OyKVCwkJCVSsWJHjx4/bHYqyUcmSJTlw4ECWScTbBWwNxO9cugQtWkCTJvDe\ne3ZHo3Ih7TfMw4cPB+73bYA7e/YsZcuWJaufXVoD8bDo6Ohc9YH4pdBQqymrfn1o2dKqjSifkNPw\ncxV40vpAvJ3WQPzNokXwwguwZQsEB+d8v7JN2m+Y+n0buLL7HvCFGoi/zER3qZ0ndnI4zkcX/e3S\nBTZt0uShlHI7v0og0dHRxMTEXHc5n2z9hBbTWrDn1B4XRGUDH+yMU74pOTmZoKAgvv/+e7tD8Ssx\nMTFER0fbHUaO/CqBjB492iX9H6/d8Rp9avahxdQWbD622QWRKaUAgoKC+Pbbb+0Ow+tFRkYyevRo\nu8PIkV91oruKMYYxd46hSP4i3DHjDhb3Xkyr8q3sDksppbyKX9VAXO3Z5s/yVvu3uHvW3Szds9Tu\ncJyXkgIJCXZHofzAiRMn6Nq1K4ULF6Zy5cp89tln6deOHz9OVFQUJUuWJCIigrp16zJ//vz06zVr\n1sQYk75YYMeOHQGYP38+DRs2pEiRIhQvXpwuXbpw8OBBTz+acoLWQHLw7/r/pnC+wpy6eMruUJz3\n3//CoUMwYwb4+qTJACMC51y8mHeBAs5/G9x3333kyZOHgwcPIiL069cvfSJucnIyDzzwAHPnziUk\nJISPP/6YPn36cOutt1K9enW2b99OUFAQS5cupU2bNullRkREMG3aNGrVqsXp06cZOHAg9913H+vW\nrXPF4yp3snsxLlcc5HI13oD1558ixYuLfPSR3ZGoDNIWy8tuEdC4OBErjbjuyMWao1k6evSoGGNk\nx44d6ee2bdsmxhhZvXp1lu+pU6eOTJgwIf21MUZWrlyZ7eds3rxZgoKC5Pz5884F6kOu9T3gK6vx\n+lUNZPTo0TqePiulSsHs2dYQ34YNoU4duyNSuVSgAMTFub5MZxw5cgRjDBUqVEg/V7FixfSvY2Nj\nefbZZ1m5ciWnT5/GGMOFCxf4+++/sy139erVjBo1it9++42LFy+m/VLI33//nan8QBIZGUnTpk3T\nFlP0WtoHEijatoXnnoN774WzXjknSWXBGIiIcO3hbPNVmTJlADL1Txw4cCC9CWvYsGHs2bOHNWvW\nEBsby5kzZ6hRo0Z6QrCeJ/OHX758maioKO6++2727t1LbGwsq1evBsj0PuWdNIFchy3Ht3A5+bLd\nYeTeiBFQoQL8+99Wa4ZSDihdujRt27blueeeS08QI0aMSL8eFxdHWFgYhQsXJjExkffee48dO3Zk\nKqNUqVLs3r07/XViYiIJCQkUKlSIsLAw/vzzT1544QWPPZO6PppAnCQiPPrlo3Sb2434y/F2h5M7\nwcEwcyaULQtJSXZHo3zQzJkzCQkJoUKFCjRs2JBevXqlX3vttde4ePEiJUqUoFKlSpw4cYIWLVpk\nev+YMWN4/fXXKVKkCJ07dyY8PJwPP/yQUaNGpY/M6qHruPkMv1oLK7cbSrnKmfgzdPq0E3mC8rC4\n12IK5ivokc9V/kHXwlLX+h7QDaU8yM7FFC8kXqD7vO78df4vlvddTvHw4h79fOW7NIEoXUwxwIWH\nhLOo1yKq3liVltNacijukN0hKaWUR2gCcYGQ4BBmdZ1Fh8od+CP2D7vDUUopj/CreSB2Cg4KZtxd\n4+wOwzkHD8Lp09ZmVEoplUtaA1GwdKk1yfDkSbsjUUr5EE0gCgYPhkaN4P77rYUXlVIqF/wqgbhq\nQylX8/qNqYyBqVNh924YO9buaJQKeL6yoZRfDeNdvDiOYsUiCAvjqiNvXnsWoj158SSV363MS7e/\nxJAmQzwfgCM2b4aWLWHZMmjd2u5o/J4O41W+PozXrzrRhwyBS5fg4kXruHTpn2vBwRAefnViufKI\niID//Adq13ZNTEXDivJ1v6/pMKsDpy6e4pU2r1y1HpDXqF8f3noL+vSB33+3/kKUUuoa/CqB/PKL\nlQDSJCVBfPw/CSXjceFC1ucPHYImTeDhh+Hll51fuTSjRjc1Ys3ANbT7pB2n40/z3t3vEWS8tPXw\n4YehcWNNHsohFStWZOTIkQwaNMjuUJQn2b2evCsOUvcDyW5fBUfs2CHSurVI6dIic+eKpKS4pFg5\neOagVHm3ivRZ0EcSkxJdU6jyWbnZD8RXVKhQQaZMmWJ3GLkWHx8v9957r1SpUkWCgoJk5MiRuXrf\nvHnz5JZbbpGwsDCpUaOGfP7555mud+7cWW666SaJiIiQ0qVLy6BBg+TUqVPXLC+774G0a3jxfiBe\n+muwvWrUgFWr4PXX4bHHoEMH2Lv3+sstX6g8awet5Uz8GY6dP3b9BSqlnGKMoXnz5nz44Yfcdttt\nuXrPhg0b6NevH2PGjOHs2bOMGjWK++67j82bN6ff8+qrr7Jv3z7i4uLS9zd58MEH3fUYttMEcg3G\nQL9+sGsX3Hyz1Sfy0kvXv7V48fDiLLtvGeUKlnNJnEp50oQJE6hcuTIFCxakVKlS12yyWr9+Pa1a\ntaJIkSLcfPPNDB8+nMTExPTrFStW5KWXXqJt27YUKFCA2rVrs2LFikxlLFu2jCZNmlCkSBGqVavG\ne++957LnCA0N5cknn6R169aEhobm6j2TJ0/m7rvv5p577iE4OJiuXbty11138f7776ffU6tWrfTy\nUlJSCAoKYs8eLx+FeT3srgK54sDFTVhZ2bhRpEEDkZtvFvnqK7d9jAoguW3CSricIHEJcVkeSclJ\nWb4nKTkpy/sTLic4He/vv/8uYWFh8ttvv4mIyIULF2TNmjUikrkJ648//pDw8HB599135fLly7J3\n71659dZbZejQoellVahQQYoXLy7r16+X5ORkmTJlioSGhsrBgwdFROTbb7+VQoUKyapVq0REZMeO\nHVKuXDmZPXt2ehmFChWSwoULS6FChdKPtNeFCxeWTz/9NFfPdfvtt+eqCatevXry+uuvZzo3evRo\nadCgQaZzw4cPlwIFCogxRsLDw2XBggXXLFObsAJEo0awYQM89RT06mVt7HfkiN1Recjnn1vVL2WL\nMWvHUPD1glkeu07uyvI9u07uyvL+MWvHOB1HnjzWmJvt27dz7tw5wsLCrtrvA2D27NlUq1aNxx9/\nnDx58nDzzTfz6quv8sEHH2S6r3///jRp0oSgoCAGDRpE3bp1mTlzJgDjxo3j0Ucf5fbbbwegRo0a\nPPTQQ0ydOjX9/WfOnOH06dOcOXMm/Uh7ffr06Ux7lbjC2bNnKVSoUKZzhQsX5uwVO3yOHj2as2fP\nsnfvXp566ikqV67s0ji8iV+NwnK34GB49FHo1g2eeQaqVoXKlaF0abjppqv/rFgRihTJXdkiwvHz\nxylVoJR7H8IZ1atD//5WFu3Y0e5oAs7wFsN5qulTWV4Lzxue5flbit5C3LCrN1MPDc5dc01WKlSo\nwJw5c/jf//7Hgw8+SNWqVXnqqafo2bNnpvsOHz7MzTffnOlc5cqViY+P5+TJkxQtWjS9vIwqVqzI\n4cOHAfj9999ZuXJlevOQiJCSkkL58uUdjvuRRx5h5syZGGMoX74827Ztc7gMgIiICGJjYzOdO3Pm\nzDXn8FSqVCl9f6IjR44QHBzs1Od6M79KINHR0R7ZUKpECfjkE2uH2L174ehR+PNP68+ff/7n9enT\n0KwZdO1qHdl973938DvunXcvy+5bRuObGrs1fodVrw6TJllLnWzenP2DKJcLzRNKaB7HfvAHBwUT\nEer6yYlRUVFERUWRkpLCggUL6NWrF/WvWISzbNmybNiwIdO5vXv3kj9//vTkAZn3Vk973alTJwBK\nlixJnz59st3etkCBAlnOqRIRjDFMnjyZ3r178/7772fqp3BW3bp12bRpU6ZzmzZtol69etd8T2Ji\nIn///TdxcXEUye1vk/yzoZTXs7sNzRUHHugDccaff4r8738id94pkiePSMOGIqNHi+zbl/X949aP\nkwKjC8g3+77xbKC59eCDIrfdJnLpkt2R+AVfG8a7e/duWbZsmZw/f15ERJYvXy7BwcGyf//+q/pA\nbrjhBpkwYYIkJibK3r17pVatWjJkyJD0sipUqCAlSpSQH3/8UZKSkmTatGmSL18+2b9/v4iIfPHF\nF1K8eHFZuXKlJCUlSVJSkmzfvl2+//57lz3PpUuXJD4+Xlq1aiXR0dGSkJAgiYnXHl7/448/Sv78\n+eWLL76Qy5cvy/z58yUsLEx+/vlnERHZs2ePfP7553L27FkREdm1a5c0a9ZMmjRpcs0yfb0PxPYA\nXPIQXppAMjp5UmTqVJFOnUTCw0Wu9f9gxpYZEv5auCz47dodb7aJjxepW1ckQ2eocp6vJZBt27ZJ\n8+bNpVChQlKwYEGpVatWeqd2xYoVM80D+eGHH6RFixZSuHBhqVChggwbNkwuZfjFo0KFCvLiiy/K\nHXfcITfccIPUrFlTli9fnunzYmJipHnz5lKkSBG58cYbpWnTprJw4UKXPU+FChUkKCgo09GmTZv0\n6w8//LDcfffdmd4zf/58ueWWWyR//vxSvXr1TPHs3r07/ZkLFCggFStWlEceeUT++uuva8agCcQL\nDl9IIBlNnixSsKDIli1ZX1+0a5GEvxYuUzZ74cSs338XqVVLJDbW7kh8nq8lEFfytYmH7uLrCURH\nYdngwQfhuefgrrtg376rr3eu1pmlfZYyfct0EpMTr77BTpUrw5YtYC3yppQKYH7Vie5Lhg+HEyeg\nfXtYuxZKXTH46vYKt7N6wGrvXHgxSH/vUNfHK7+vlcM0gdjEGHj7bTh1yqqJrF4NVwwx1/9kym/t\n37/f7hCUC+ivkjYKCoIpU6BsWWjRAj76CM6dszsqpZTKHU0gNsubF+bPh0cegf/9z2rKGjTImvWe\nlcTkRC4kXvBskDlJSbHa45RSAcUnEogxZogx5mu743CXfPmsbck3b4Y1a6ytONq1gx49rP1JMnrr\nh7doP7M9Z+LP2BNsVqZNg7ZtrQ1VlFIBw+sTiDEmD1AHazib36tXDyZMgD17rB0Uq1eHV16xNsYC\nGNpkKEXDitJ6emuOnfOSJeHvv9/ayeuJJ+yORCnlQV6fQIB+wGy7g/C0kiWtX+xXroTFi6F5czh2\nDPLnzc+CHguoV6oeLaa1YP8ZL+iMzJsX5syBRYtgxgy7o1FKeYhXJxBjDUNqLyJfAwE5JKlJE/jh\nB6hVy/r6t98gT1AepnWZRlTVKFpMbcG2v5xbHM6lypSBWbOsHbh27LA7GuVhFStWzLRSrgoMXp1A\ngK7AYruDsFtICEyfDgMGWDWRjz+GzT8H8Vztd3i44SOMXTfW7hAt7dvD0KHWWvfnz9sdjVLXlJCQ\nQI8ePahatSrBwcH897//dej948ePJygoKNP7Dh8+TIECBYiIiEg/8uXLR968eTl9+rSrH8EreGwe\niDFmPNAZKA/UFZGtGa5VBmYARYFYYICI7ASqAa2NMfcDdYwxg0QkIH/NMQZeftlaIv7tt63Vfk+e\nNHTrNpI5nyXbHd4/XnzR6g/RyYbKi6VtaTt48GCGDx/u0Ht3797Nu+++S+3atTOdL1u2LOeuGIff\nrVs3kpKSHFqJ15d48n/5PKA5cDCLa5OBSSJSDXgDK5kgIqNFJFJEOgBbAjV5ZDRgAPz6qzVq9u+/\nYdMmGPeOF+0zEBxsbZYSFmZ3JMoNAnlLW7C2qe3fvz/vvPMOhQsXzvbeo0ePsmTJEgYPHny94Xov\nTy++BRwAamd4XQyr1hGU4dwxoJIDZfrUYoqutGGDSFiYSOrOosqH+NpiioG+pa2IyKhRo6Rv3765\net/IkSOlcuXK2Zbn64spesNSJmWBYyKSkuHcIaAc4NAQo+joaEJCQgCIjIx0+8ZS3qBxYxg7Fnr2\ntDrb0/Z6ir8cT/68+e0NTrnGpUvWkZXwcKvWd6XkZLiQxYTT0FDrcELGLW3LlClDgQIFctzSFkjf\n0rZv37783//9X/p9aVvaAgwaNIgPPviAmTNnMmLEiGy3tO3duzdg7QboSVu2bOGjjz5iy5YtOd6b\nlJTElClTePrppx36jJiYGGJiYgAy1di8lTckEJcZPXr0NbeX9GeDB8Pu3dCypTXst0oV6Dm/J9WL\nVuf1O1/XNbV83ZgxVgdYVrZvh1tvvfr8rl1Qs+bV51980en97QN5S9ukpCT69+/PuHHjrtoXPSuf\nf/45sbGxDBw40KHPyfiL79mzZ5k4caLDsXqUp6s8uLEJa/DgwVdtShMoUlJEnntOpEQJkZ9/Ftlz\nco+Uf6e8/HvRvyUpOcne4A4cEJk7194YvFCum7ASEkTi4rI+kq7xb5uUlPX9CQkuiT05OVnmzp0r\nQUFBsmfPnkxNWGPGjJEGDRpkun/hwoUSHh6e/rpChQry7LPPZrqnSZMm8uqrr4qISJs2bWTUqFHZ\nxnDDDTdIgQIFrjrSzmds7spObpqwDh48KEFBQVKsWDEpWrSoFC1aVPLmzSthYWFSs2bNLMscOHBg\njp99re+B5cuXy+DBg72+Ccv2BJJ67lugf+rX3YGNDpYZsH0gGaWkiIwdK5I/v8hLL4nsP3lYqk+o\nLt3ndpeEy675weGUtWutjpoNG+yLwQv5Wh9IIG9pm5KSIkePHs10NG3aVIYMGSLHjx/PdO+OHTvE\nGCObNm3KMQZf7wPxZOKYBBwGElNrGHsyXKsK/ADsBjYCtzpYtiaQDH76SaRmTZGyZUUee+6E3Dqu\nkbT7uJ2cu3TOvqDefFOkfHmRU6fsi8HL+FoC0S1tM2vTpk2WNZfHH39cGjVqlKsYfD2BGBHfX2LK\nGBMBxA0ePJioqKiA6DzPyeXLsGIFzJ4NC5edw/S+h4Y3dKZvlSfp189awNGjRKBLF+vPRYt0nghW\nG3fBggWJi4sLuL67ihUrMnLkyGsOAw4U1/oeiImJYcmSJWl9IAVF5KxtQWbDrxJIIP5HzI0LF2DB\nF5fY/FNeln8VROPG1pJVHu9bP30a6te3ev2ffdbDH+59NIFoAsnueyDtGl6cQPTXwAAQHg733xfK\nuHeCWLXKGqn12mvWSE+PKlIE5s61RhRduU69Cig6MtA/+NUw3ujoaG3CykGpUvDFF9Ctm7WmVrNm\nUKAA3HwzPPigByaQN24MO3da2zCqgKVb2mYvrQnL22kTVoBKTLR2QtyzN4nzZ/OwbJm1J5QLV4tQ\nOQjkJixl0SYs5ZNCQqBrjwTmF61L24e+YvFia/+RMWNg3Tq7o1NK+QJNIAEsX558jGw1ku7zuvNT\nwhymTrWSR2SktXx8SkqORSilAphfJZDo6Oj0dWRU7vSs2ZMFPRbwnyX/4XSlSSxdCgsXWtt63HST\ntSKGR/hBU6pSrhITE0N0dLTdYeRI+0AUAOsOraPTp514ttmzDG8xnORkw9ChsGULrF7t5mkb330H\nb75p7d2b1cKAfiqtjfvw4cP6fRugzp49S9myZX22D0QTiEq39a+t3DXzLj7v+TlNyjThwgWoWxda\nt4bx463hwG5x/jw0agTdu8OoUW76EO+TkJBAxYoVOX78uN2hKBuVLFmSAwcOkO+K2b2aQDxEE4jr\nxCXEUTBfwfTXR49aE8i3bYMPP4T773fTB2/fbm36vmCB1QkTIBISEnxi2W7lPiEhIVclD9AE4jGa\nQNwvJga6doW1a6FePTd9yIwZ1m6Gv/wCZcq46UOU8g2+kEC0E13lSmQk3HcffPSRGz+kf3/o3Bl6\n9bIW81IqQGknugdpDcT9RIQ1awxdu8KRI25cjPHiRejRA959FypVctOHKOX9tAai/MLuk7tpOqUp\nFWofoVw5a3jvgQNu+rCwMFi6VJOHUj5AE4jKUeUilalVvBYtpzdndswebr8dJk+2OyqllN00gagc\nBQcF80HUB/Su2ZtW01vQrt8vTJkCuh6eUoHNrxKIdqK7jzGG1+98nWeaPcOwXW1oM+B76taFQYPg\ns8+0z1spV9JOdA/STnTP+vDnDxm2chizGu9jdUwhPvnE6vv+8UeoWtUNH5iSArt3Q/XqbihcKe/k\nC53omkCUU46fP07JG0oC1s/3oUOtJU+++QaKFnXxh33/PURFwc8/Q+XKLi5cKe/kCwnEr5qwlOek\nJQ+w1sl6800oVw4ef9wNH9aqldVWdu+9kJDghg9QSjlDE4hyiZAQazOqhQth40Y3fMDYsRAaCkOG\nuKFwpZQzNIEolylfHkaP/mftrPPnXVh4SIjVWz93Lsya5cKClVLO0gSiXGbpnqVsqzSQnn0uU7++\ntdf6u++68APKl7c2cn/4YatTXSllK79KIDqM114NSzdk8/HN7G/YjZOx8cyaBe+8AxcuuPBDOnWy\ntkssW9aFhSrlXXQYrwfpKCzvEZsQS8fZHckblJf53RbTo0sEu3bByy/Df/5jd3RK+Q4dhaUCTqF8\nhVjRdwX58+Yn8tM2fLroBBMnWqOzDh2yOzqllCtpAlEuFx4SzqJei6hSpArd53flnnuEe+6BJ56A\ns175e5RSyhnahKXcJjklmT/i/qBS4Urs3An33AMi1uaDISF2R6eUd9MmLBXQgoOCqVTYWpa9enXY\nudOatf7GGy6eD3j4MLzyipWdlFIeowlEeUxQEEyZYh3Fi8PEidYqJSkp11lwSAi8/76uMa+Uh2kC\nUR7VujXs3Quvv27tG9Wli7XM1alT11FoiRIwZw48/TRs3uyyWJVS2dMEojzu/358kxuaf8xXX1kz\n1kWgQgUYPvw6WqFat4YRI6z1suLiXBmuUuoaNIEoj6tfqj6PLXuM8T+Op0wZ+PJLaxXfjz+GTz65\njoKHDbPWkx80SPtDlPIAv0ogOhPdN7St1JZv7v+GUd+P4r+r/gsIt91mzVp//nk4csTJgoOCrAz0\n00/www+uDFkpj9KZ6B6kw3h9028nfqP9J+2555Z7eLfDuxiC6Nnzn9pIp05OFnzunLUQl1I+TIfx\nKpWNGsVqsG7QOlbsW8EL376AMdZiuxMnWh3rDz0EsbFOFKzJQymP0BqIst1f5/8iKSWJmyJuSj/3\n668wYADceCN8/TUYY198StlBayBK5UKJG0pkSh4AderAqlXWKK277oING2wKTil1TZpAlNcqVAg2\nbYJbboEmTWDJkuso7PJll8WllLJoAlFerVw5GD/e2kNqxgwnC9m1y1pL5e+/XRqbUoFOE4jySimS\nwjMrnuG3E78B8Oyz8N13MGGCE4VVqwYNG0LfvpCc7NI4lQpkXp9AjDGNjDFrU49X7I5HeYbBEJ43\nnFbTWrHx6EYqVYLFi2HoUKuD3bHCDHzwARw8aG3arpRyCa8fhWWMCRaR5NSvvwHuEZHzV9yjo7D8\n1PgfxzNy1UgW9lxI20pt6dIF/vwT1q1zYkn4X3+F5s1h0SJo29Yt8SrlKjoKywUyJI9g4Chw0d6I\nlCc92eRJJtw9gS5zurBw50KmTLG6MpYudaKwOnWsDpU+feDYMZfHqlSg8foEAmCM6Q38BsSKyPUu\n/q18zP117md2t9n0W9iPTWe+ondveOwx+OMPJwobNMiaYKId6kpdN69vwsrIGDMfeFFEdlxxXpuw\nAsCGIxuoWbwmoUHh3HcffPuttX5W3752R6aU62kT1hWMMeONMQeMMSnGmNpXXKtsjFlnjNltjNlg\njKmeej6L6TdqAAAgAElEQVRjS/c5wJV72SkfcluZ2wgPCSdPHmv7j2eegccf14V3lbKLp5uw5gHN\ngYNZXJsMTBKRasAbQNqo/87GmFXGmO+AwyKyzxOBKu9mDDz5JFy6ZC11opTyvDye/DARWQtgTOaV\njYwxxYAGQLvU+xYYYyYYYyqJyHxgvifjVL4hXz548UXo1cuaI1K7do5vUUq5kLd0opcFjl3RQX4I\nKGdTPMpHDB5ykWqPDqdz9wvOzxFMSYGVK10al1KBwKM1EHeLjo4mJHVyQGRkJJGRkTZHpNwthSRC\nbl7Pn23b8/70pTz2QGHHCzlyBLp2tTYh6dLF9UEqlUsxMTHpm+IlJibaHE3Ocj0KyxgT5KohtMaY\nA0AXEdma+roY8DtQJO0zjDHHgOYisj8X5ekorAAWfzmeFuN78suBgyzuEUOn1qUcL2TuXGsDkp9/\nhkqVXB+kUg7yt1FYB40x0ak/7F1KRE4Am4F+AMaY7lgd5jkmj4x0S9vAlD9vfn4cuoBbCtWl25ct\n+H6bQ982lh49rPHAPXpYPfNK2cjvtrQ1xlQDHgF6AyuACSLi0C4NxphJQEegBHAKOCciVVOvVQWm\nAzcCccDAK+d7ZFOu1kAUf/2dQotXh7IvdB4PBW/gnZfKki+fAwVcugQtWlhrx7/3ntviVCo3fKEG\n4vBEQmNMGNAXGAH8DbwDfCo2zkjUBKLSiAgvfvY580f9i1trBDFnDgQHO1DAgQNQv77VHxIV5bY4\nlcqJLyQQh0ZhpQ6/bQ/0BM4Dn2LVSD53fWiO0yYsZYzhlV7d+H51EFu3wn33wfHjDhRQsaI1sUQX\nW1Q28scmrOHAg8AO4F0RWZHh2u8iUsU9IeYqNq2BqKscPgz/+Y+1GaGO0lW+xhdqII4M4y0NRIrI\nniyu9XJRPEq5TNmy1pInhQtb+6u3aWN3REr5F0easN7EmtwHgDEmvzGmLICI/OzqwJyhTVjqSoUK\nWZtQDXziGBM3TLI7HKVyxR+bsDYCrUQkIfV1fuA7EbnNjfHlijZhqexcvgzVW2/j0O138tBtA3i3\n8+tcsZqOUl7HF5qwHKmBhKQlDwARiQdCXR+SUq6VNy/8+nUt2hxYy+R1n9Fx8oMkpziw7snJk9Cv\nH5w/n/O9SgUQRxKIGGOKp70wxpQE9Nc45RPCw2Hx9Co8mm8tX21fR7N3enEpKZcTBgsVsnrkH31U\n145XKgNHEsi7wHpjzEvGmJeAtcBbbonKSdoHorITGgrjXinDO7W/Z/O+P2g+MYqExKSc35gnD3z6\nKcTEwJQp7g9UBTy/6wMBMMbcDtyd+nKJiKxxR1CO0j4Q5QgRGDrsHFN/+ILysf3YtInczVhfuRI6\nd4YffrD2V1fKjXyhD8SntrS9Fk0gyhmJidCgAdx4IyxcaA33zdErr8DMmfDTT6Dfa8qNfCGB5LoJ\nyxhTwRjzvjFmhTHm27TDncEp5U4hIbBmjZVIeveGhNxsljxihDVbfdkyt8enlLdzdBjvSmA9kD6E\nRUS+dE9ouZdWAxk8eDBRUVG6D4hyyMmTcOedULIkLFlijdrKVlKS1S+ilJvExMSwZMkSJk6cCF5c\nA3EkgWwVEa/cNFSbsNT1OnXKWv6qZN0tPDb0Ap3qNLc7JBXg/KoJC9hujNEtZpVfuvFG+Oor2H9x\nC53n3cXUNcvtDkkpr+dIAikG/GqMWWaM+TztcFdgSnlaqVKwY9YA2pz7kAdiujF0yhyd9qFUNhxp\nyJ2Zeijlt/LmhZXjezF8akFe33cvX/aNZfVbD1Mqp11yY2OtCYdKBRBnNpQKFRGv2vNT+0CUO8Ts\nXEen2Z2ofPx5dn447No3njkDVavC4sXQtKnnAlR+za/6QIwxtYwx24F9qa8bGGPecFtkTtCZ6MqV\nIqs3Z9m9q9n1QxXWrcvmxsKFreG9PXpYQ7qUuk5+NxPdGPMd8ALwnojUS92dcLuI3OrG+HJFayDK\nnYYNsyYafv89lChxjZtEoFs3azLJ0qUQ5NBmn0pdxa9qIMANIrI27UXqHuiJrg9JKe8yapQ1d7BG\nDfjuu2vcZAxMnQq7d8PYsZ4MTynbOJJAkowxeQEBSN1MyoE1sZXyTXnzwpdfwgsvWEth7d59jRsL\nFYJ58+DVV2H1ao/GqJQdHEkgE4AvgGLGmFeBNYBX9YEo5S7BwdbOhv37Q/v2MHv9t5y8mEV/R/36\n8MEHuVydUSnf5uhqvM2ALlj7gCzO2KRlJ+0DUZ5y+TIMHAify/2Uve0nVg5YQZmIMnaHpfyQL/SB\n6Gq8SjkoKQnu7pjMhuIPE1F3BSsHfE3VG6vaHZbyM76QQBwZxrsq4yq83rgarw7jVZ6QJw98tSyY\nyMQPSPq1Fy2ntmTL8S12h6X8iD8O4+2Y4WU+oA+wR0SGuyMwR2gNRNkhPh7atYPTNcbyZ6UxLOm9\nhJblW9odlvITvlADyfVSJlcu226MWQR4VQ1EKU/Knx/mz4fy5Z+n5ZNFib9wjXXgU1Lg44+hb19d\nBl75leuZ7RQMlHZVIEr5opIlYf16YPMD3F27Ce+/n8VNCQnw1lswcqSnw1PKrXL965AxZiGpc0Cw\nkkdtQLdlUwGvfn345htYtAh69oRLl+DJJ625hQCEhVnzQxo3hhYtoGPHbMtTylc40gfSP8PLJGCv\niGxwS1QO0j4Q5S2WLoVevaBLF/jwQyt3pJs1C554AjZvhvLlbYtR+QZf6APRYbxKudiff1oJxBgY\n/OYq+jRvTd48qa3FDz0Ev/5qLawVEmJvoMqr+VUCMcZMze66iAxySURO0ASivE1CAjz02EU+Dq9N\nwbPN+OXlKVQsl9e60LSpVU15/nm7w1RezBcSiCOd6IlAE2A/1pLujYFLwM+ph1IqVb58MOOjMHY/\nv4aU4luoMqIbsz6Lty4sWQKPP253iEpdN0dqIN8DndIyYepv/UtFpJUb48sVrYEob3b64hka/F9H\nDu4LoeulxcycEkH+/HZHpbydv9VAimV8iNSvi7k+JOfpTHTljYqEFWb7c1/Tokk+lhZtQ7d+J0hK\nsjsq5c38cSb6HCABmJJ6aiAQJiK93BRbrmkNRPmCxOREHlrwDF/992nqVyrPokXWUvFKZcUXaiCO\nJJAbgBeBO1JPfQOMEpHzboot1zSBKF9y9Cg0bGhto75okbWNiFJX8oUEkusmLBE5LyLPikiD1ON5\nb0geSvmam26CTZus7dPvvReS4xOtTUauuVOVUt7JkdV4yxpjlhpjtqS+rmuMGeq+0JTyX2XKwKpV\nsHUrdOoaQlLtelY2uXjR7tCUyjVHOtEnA3OwNpMC2A7YNvdDKV9XvDisWwc7dwlNtnUjpUBBa6a6\nUj7CkQRSXERmAikAIpKEtaSJUspJlSvD9GXb2dzoDloX7k7KF4tgxgy7w1IqVxxJIEnGpC8PhzGm\nMP/URpRSTrq9ei3mdVvED3VH0uPm7qQMfgy2b7c7LKVy5EgCmYfVjBVhjPk38DXwkVuiysAY09wY\ns94Ys1b7XJS/6lavLd8NWskXbefxVtn6JK1YaXdISuXIkVFYbwOrgJ+A9sD/icgEdwWWwT6gpYi0\nAKKMMfk88JlKeVzLSo1Y/+D3RHfaS7WDv5MiKXaHpFS2cpVAjDHBxpixIvKpiPQUkR4iMtvdwQGI\nyPHU/haw+lz0f5XyW40q1CCm1zoOHD/JsJHxdoejVLZytaGUiCQbY9q4O5jsGGPaAftEJNHOOJRy\nt7YNKvDto3No0wZKFoGnnrI7IqWy5kgfyDJjzAhjTGljTETa4bbIMjDGlAGeB/S/kgoIt99uzVJ/\n+mlYvdruaJTKWo4JxBhza+qX/wVGAUeA2NTjjCMfZowZb4w5YIxJMcbUvuJaZWPMOmPMbmPMBmNM\njdTzIcA04BER0Tq9ChidO0N0NHToAFu+OAB+sPmb8i+5qYF8kvrnehEJuuIIdvDz5gHNgYNZXJsM\nTBKRasAbwPTU832A6sBkY8y3xphSDn6mUj7r1Veh173JhP0rkkPD3+HI2SN2h6RUuhwXUzTG/Aa8\nDIwGhnDF3A8RWezwhxpzAOgiIltTXxcDfgeKiFhDT4wxx4DmIrI/F+XpYorKr43ptI4nvrqDfz1U\nkHEvf0eNYjXsDkm5mS8sppibTvRhwMNAca7ugxDA4QSShbLAsbTkkeoQUA5rB8RciY6OJiR1n+nI\nyEgiIyNdEJpS9hu2pDnL7niVD6a9xu1FWjD3keU0vqmx3WEpF4uJiUnf0ygx0fvHCzmynPt4EXnS\nJR96dQ2kPjBLRKpnuGcD8LyIfJeL8rQGovyfCNsrd+Fg7H56Df2Dxfct4o6Kd+T8PuWTfKEG4shE\nQpckj2s4DJQyxmSMpxxWLUQpBWAMt26cTuOE8zzyv450nNmZL3Z9YXdUKoA5MozXbUTkBLAZ6Adg\njOkOHM5N/0dGuqWt8nfmxiIU+24+JZoNIWHWbKatX2h3SMoN/G5LW5d8mDGTgI5ACeAUcE5EqqZe\nq4o18upGIA4YKCI7clmuNmGpgCIC//oX/PIL7NsHeXI1JVj5El9owvJoAnEXTSAqEMXHQ9mycOed\nMGeO3dEoV/OFBOIVTViuok1YKpDkzw9ffw2ffQZPurOHUnmcNmF5kNZAVCBbtQruuMNKIm+/DbGX\nTlEgtAAhwSF2h6aug9ZAlFJu1+Z24deB7zBn/HEaN4ZnVjxDlzlduHhZ91dX7uVXCUSbsFSgqp34\nM/sa92bL5mRi54wn/nI87T5px5l4h5arU15Cm7A8SJuwVMA7fx4aNWLnrd2psWAUI16KZ2vVnvwR\n9wcxfWMoeUNJuyNUDtImLKWUZ9xwA8yfT/Xl7zD3geW89lJ++udfQJ0SdWgxtQUHzhywO0LlhzSB\nKOUvbr0VJk7k3i/6MuDOI3T/V14G3zSdjlU6Mme7jvNVrudXTViDBw8mKipKF1FUge2BB0hKgtqb\nprBzJ+zeLVSpAsaYnN+rvEJMTAxLlixh4sSJ4MVNWH6VQLQPRCmsGYYiJIeGcdNN1qz1o0d1trqv\n0T4QpZTn5c8PYWEEB8OaNfD339C1q91BKX+kCUQpP1alCsydC0uWwBNPZL4WlxBnT1DKb/hVAtF5\nIEpd7d57Ydw4eO89eP5561xySjLNpjbjjXVv2BucypLOA/Eg7QNRKgfJyYx5I5joaBg7Fp57Dn45\n9guRMyMZVG8QY9qO0U52L6N9IEop+4lAx44MrzKfJ56waiHLl0O9UvVYO2gtn27/lIeWPkRySrLd\nkSofowlEKX9nDDz4IDzwAOMf30vTptChA/z+O1S9sSrrBq1jzaE19FrQi0tJl+yOVvkQTSBKBYKu\nXWHQILj3XtZ8nUChQlC1Kpw5A2UiyrBm4BoOxR1i1cFVdkeqfIj2gSgVKBIToVUrqFuXM2MmUaSI\ndfrkSbjxRkhKSSJPkE4W8RbaB+JhOgpLqWyEhFi7T82dS+Flszh0yDpdvz6kpKDJw4voKCwP0hqI\nUg748kv46Sd48UW2b4dateCuu2DZMqu7RHkHrYEopbxPx47w4osA1KwJX3xhjcp6882sb09ISvBg\ncMqXaAJRKsB16UL68N7x4zNf+/Pcn1R+tzI/HP7BnuCUV9MEopRi3Djo1w+GDIENG/45X7pAaV5o\n9QKRMyNZvne5fQEqr6QJRCmFMTBlirWlSJMmsGfPP9cebvgwH0Z9SLe53fhs+2f2Bam8jg67UErB\n9u3kLV+e1asLUL48VKtm7ZIbHm5d7lWzFwVDC9Jjfg9iE2J5qOFD9sarvIJf1UB0GK9STho+HB58\nkBuLCNu3W6eaNYPTp/+5pUOVDiy/bzmvfP8KR88etSfOAKHDeD1Ih/EqdZ3++gvq1rVGZz38MOvW\nQYsW0KABfPstZPxvFX85nvx589sXa4DQYbxKKd9QogTMmQNPPw2bN9O8uTUv5OefoXdvaz3GNJo8\nVBpNIEopS+vW8MIL1gYisbF06AAzZ1qJ5L337A5OeSNtwlJK/SMlBTp1ggIFrGVPgF69rC/XrbP6\nRa4lOSWZ4KBgDwXq/7QJSynlW4KC4OOP4amn0k9NngyVKsHjj8PevVm/bcW+FTSf2pxTF095KFDl\nDTSBKKUyK1oUbrst/WXBgtZyJ7/9ZtVGErJY2aR1+daUiShDy2ktOXL2iAeDVXbSBKKUylGtWjB3\nrtWpPm7c1ddD84TyWffPaFa2GS2mtuD3U797PkjlcZpAlFK5EhUFAwdaU0ZOnLj6enBQMB9GfUjP\nW3vSYloLthzf4vkglUf5VQLRiYRKudf770Pp0tCyZdbXjTGMbTeWp5s+Te8FvXWfdSfpREIP0lFY\nSrnZW29Z/SItW7J3L1SpYtVIFi++9ltiE2IplK+Q52L0MzoKSynlH/Lnt3rQ//6bypVh4UL45hsY\nNerab9Hk4f80gSilcvboo9baJn37QnIyd98NTz4JL71kjc5SgUkTiFIqZ8bAhx/CwYMwejQhIVby\nqFbN2tUwLi53xYgI/tBsriyaQJRSuRMRAfPmwdixsHIloaGwYweEhFi75Kat4pudN394k6ExQ0mR\nFPfHq9xOE4hSKvfq1LEmgowZA1gVkyVLrBrIwoVw4UL2b+9eoztL9ixh4KKBJKUkeSBg5U6aQJRS\njnngAfjyy/SX7dpBhw7w2mvw8svZv7VS4UqsHbiWX479Qre53UhIymJau/IZmkCUUo4xBkJDM516\n4w1rId8VK6xV4bNTqkApVg9YzYkLJ+gwqwNnL3nlCFWVC16fQIwxhY0xPxlj9LtMKS8WFQUVKsCw\nYTnfWzh/Yb7u9zUhwSE8/tXjbo9NuYfXTyQ0xgQDBYC5ItL+GvfoREKlvMD27VC/vtWpPmNG5p0M\ns5KYnMjZS2cpGlbUMwH6EJ1I6AIikiwisYCxOxal1DUkJsLly1SvDh99ZE0y3L8/57eFBIdo8vBh\nXp9AMvDuqpJSgeyRRyA6muBguP9+a4fcfv2s/dSV//KlBKKU8lbDhsEHH8CiRQB8+inceCOsX29z\nXMqtPJpAjDHjjTEHjDEpxpjaV1yrbIxZZ4zZbYzZYIypfuXbPRiqUsoRVapYM9UHDID9+2nUyJqh\nPnMmvPqq48U9sOgBFu/OZqVG5RU8XQOZBzQHDmZxbTIwSUSqAW8AM9IuGGO+BuoaY1YYY2p4IlCl\nlIN69LDWyurRAy5d4tFHoX17mDTJ8aI6Vu1I7wW9+eTXT1wfp3IZW0ZhGWMOAF1EZGvq62LA70AR\nEWuNA2PMMaC5iOTYFaejsJTyEpcuWYsu3nYbTJjAtm3QpInVutW0qbW3em6t3L+Sf332L0a1GcWT\nTZ50X8xeyhdGYeWxO4BUZYFjackj1SGgHJCLsRyW6OhoQkJCAIiMjCQyMtKlQSqlchAaau19O306\niFC+vKFlSxgxwqqNfPBB7otqW6kt39z/DR1mdeB0/Gleuv0ljPHvluyYmJj0TfESExNtjiZn3lID\nqQ/MEpHqGe7ZADwvIt/lojytgSjlxd54A375xepcd9RvJ36j/SftebH1i/ynwX9cH5yX8oUaiLeM\nwjoMlDLGZIynHFYtRCnl4yIi4KuvoGFDWL3asffWKFaD9Q+sp2/tvu4JTjnNKxKIiJwANgP9AIwx\n3YHDuen/yEj3RFfKO/XpAx9/bC39/ssvjr+/bMGy5M+b3/WBeSndEz2rDzNmEtARKAGcAs6JSNXU\na1WB6cCNQBwwUER25LJcbcJSygf06wfVq4MP/Gy0nTZhXUFEHhaRsiISIiKl0pJH6rU9ItJMRKqJ\nSOPcJg+llA/Ytg127qRwYWsnwwIFYOpUu4NS18srmrBcRZuwlPJS8+ZBt268PvICW7ZAp06wd+/1\nFXkh8QL3fX4ff8T+4ZoYvYg2YXmQNmEp5eWSkuDOO6FcOZgxg6eeNojAO+84X2SKpPDEV0/wxa4v\nWNFvBTWK+dccY23CUkopgDx5rDG8MTEwZQphYfD553D33c7NVAcIMkG81+E9BtUbRKtprdh0dJNr\nY1Y58qsEok1YSnmxUqWsJDJkCINb/MqIEdaCi6nrLzrFGMMrbV7hhVYv0Pbjtnx7wD+W/9UmLA/S\nJiylfMgrr8Bnn8HWrXw8K5jp012z7PvHv37Mo18+yrpB66hTss71F2gzX2jC8palTJRSgWLECPjX\nvyA4mNBQSEhwTbH317mfW4vdSq0StVxToMqRXzVhKaV8QHAw1LJ+yBcqZO0ZEhJiHW+/fX1FNyjd\ngCCjP9Y8xa/+prUPRCnf0r49HDwIu3ZB//7wh/+NyHWK9oF4kPaBKOX7RoyAkydh8mS7I/EOvtAH\n4lc1EKWU7woNtbYTcbXDcYcZsnwIicnevzy6r9EEopSy3+TJNN80jh9/hKFDrePzz11TdGieUL7/\n43u6zOnCxcsXXVOoAvwsgWgfiFI+qm5d2qwcwSN112MM7NwJEye6puji4cVZ1X8VFy9fpN0n7TgT\nf8Y1BbuR9oF4kPaBKOUHxo2zhmH98guzYooyeTJ8/73rio+/HE/P+T35I+4PYvrGUPKGkq4r3A20\nD0QppXLrySehUSO4/35C8qRw+bJri8+fNz8LeiygTok6tJzWkguJF1z7AQFIE4hSyjsYY63xvns3\nNZaOxR1bgucNzsv0e6YzqeMkwkPCXf8BAUYTiFLKexQqBPPmUXrTIi6dS2T3btKPiy7q/w4yQbSt\n1NY1hQU4v0og2omulB+oX5+Ds9dz6HgI9epBvXpQsyY89ZTdgXmOdqJ7kHaiK+XfxoyB336DTz6x\nOxLP0U50pZRygTx5rD2p3Gnj0Y1M/UX32XWEJhCllNfLkweXj8q60qWkS3z5+5ekSIp7P8iP6HLu\nSimvF0Ii+S+cBYq67TNalm9Jy/It3Va+P9IaiFLK6zVc9QbPb+zq/nYs5RBNIEopr7czcihh8af5\n+6GRbN1K+nHunN2RBTa/SiA6jFcp/3RT1XD65J1H/qkTeK3Zl7RuDY0bgw+MdHWKDuP1IB3Gq1SA\nmD0bHn8cfvmFl6aW4/BhmDLF7qDcQ4fxKqWUK/XpA/feCz17klcSSU62O6DApglEKeVbxo2Dnj0J\nyhOkCcRmmkCUUr4lXz4YMgSTN48mEJtpAlFK+aTgYDSB2EwTiFLKJ2kCsZ8mEKWUT9IEYj9NIEop\nn5QpgWzdCkuW2BpPINK1sJRSPik4GA4ftqaG3PTzQZr9ry/LR23iXOlqgLU7bpUqNgfp5/wqgURH\nRxMVFUVkZKTdoSil3KxePYiIgEmTADrzcNGHue3Fe3m47o/8fjSMVq3g44/tjtI5MTExLPGBGpXO\nRFdK+YfLl6FNG6hWjderTGHbNpg1y+6gnKcz0ZVSylPy5oU5c2DxYmptnoEf/G7s9TSBKKX8R5ky\nMGsWdywdSmhCnN3R+D1NIEop/9K+PTOe2c6FPAXtjsTvaQJRSvmd+MKltQnLAzSBKKX8jjGQolub\nu50mEKWU3wkK0gTiCZpAlFJ+JyiIq5uwtE3L5bw+gRhj3jXGrDbGDLM7FqWUb7iqCWvePGszKk0i\nLuXVCcQY0wC4LCKtgQbGmGJ2x6SU8n5XNWHdfjusWQMTJtgVkl/y6gQC3AZ8m/r1aqCBjbEopXzE\nVU1YxYrBZ5/BsGGwcaNtcfkbb08ghYC0KfznUl8rpVS2suxEb94cXn4ZevSA06dticvfeHsCiQXS\nFrcqkPpaZRATE2N3CLbQ5w4sjj73NYfxPv001K4NAwZof4gLeCyBGGPGG2MOGGNSjDG1r7hW2Riz\nzhiz2xizwRhTPfXSRqBN6tetgJ89Fa+v0B8ogUWfO3eyHIUFVmaZMcO6ePKka4ILYJ6sgcwDmgMH\ns7g2GZgkItWAN4AZACLyE5DPGLMa2CIiJzwRqCPfrNnde61rWZ2/8lx2r931Q0Sf+/ru1efO+byn\nnju7eSAxGzdam08Vy3lMjq89t6d5bD8QEVkLYIwxGc+njqxqALRLvW+BMWaCMaaSiOwXkcdy+xln\nz7pmxeMlS5bQtGnT6773WteyOn/luexeZ/w6MTFRn/s66XNf373e+NyXLsHOnTB48NXX1q1bwuLF\n137uBx6AypWzf7bsrrnquV317+xOHt8PxBhzAOgiIltTX9cHZolI9Qz3bACeF5HvclnmTcARN4Sr\nlFJ2KyMiR+0OIiv+siPhn0AZrJFaSinlLwpg/XzzSt6QQA4DpYwxQSKS1mpZDjiU2wLEqkZ5ZYZW\nSqnr4NXtWLYP403tGN8M9AMwxnQHDovIflsDU0oplS2P9YEYYyYBHYESwCngnIhUTb1WFZgO3AjE\nAQNFZIdHAlNKKeUUj3eie5Ix5l2gDvCViLxudzyeYIwpDHwNVBWRiJzu9yfGmObAW0AysEBE3rE5\nJLczxjQC0p7zWxH5r53xeJoxZgjQUUTa2R2LJxhjygPrgV3AIREZYGc8tjdhuUsAL8R4FrgT+NHu\nQGywD2gpIi2AKGNMPrsD8oDNItIi9ZmbGWNusDsgTzHG5MH6BdF/fwvO2lIRucPu5AF+nEAI0IUY\nRSRZRGIBk+PNfkZEjotIUurLJMDvtxQSkWQAY0ww1kCSi/ZG5FH9gNl2B2GDu1K3uOhjdyD+nEAC\nfSHGQPutLJ0xph2wT0QS7Y7FE4wxvYHfgNgMIxn9WuqE5PYi8jWB9cvSn0BVIBJ4KLXJ2jb+nEB0\nIcYAZIwpAzwPPGV3LJ4iIp+mLgN0kzHmVrvj8ZCuwGK7g/A0EbksIgkikgCsAW62Mx6fSCCBuhCj\nk8+dfosHQ3U5B5+9Rur5EGAa8IiIxNsR9/Vw5t879ZnTnAMSPBmzKzj5fV4NGGCM+QqoY4wZ5PHA\nr5OT/97hqX8aoCHWPDr7iIjXH0ALoDSwH6h9xbWVQL/Ur7sBGzNcm4DV/xFt9zN4+Lm/Bk4AK4Aa\ndj+Hp54dGIC1pM23qUcpu5/DA8/cHVgFfAe8YvczeOq5r7hnhd3P4MF/73bAJmAdMNT2Z7A7AAf/\nwtHOTjAAAAJNSURBVA9k/IsGimE1TQVlOHcMqGR3rPrc+uz6zPrc/v7cPtGElY2ywDHJ3HF4CGsp\nFH8WqM8NgfnsgfjMoM/t9c/t6wlEKaWUTXw9gaQvxJjhnEMLMfqoQH1uCMxnD8RnBn1ur39un04g\nEqALMQbqc0NgPnsgPjPoc+MDz+0Ta2EF6kKMgfrcEJjPHojPDPrc+PBz+0QCUUop5X18uglLKaWU\nfTSBKKWUcoomEKWUUk7RBKKUUsopmkCUUko5RROIUkopp2gCUUop5RRNIEoppZyiCUQppZRTNIEo\n5SLGmP7GmIV2x6GUp2gCUSoDY0zwdRahawOpgKEJRAW81D2pXzLGbARGG2NqGmPWGGN+MsZsN8ZE\nZ7j3RWPMHGPMYmPMDmPMN8aYQlmUWdoYs9EYM8CTz6KUJ+WxOwClvMRlEWkMYIwJB+4QkcvGmHzA\nD8aYb0RkY+q9jYH6IhJrjPkUeAgYm1aQMaYmMAd4UkRWevYxlPIcTSBKWaZl+DoMeN8YUxdIAcoA\ndYG0BLJcRGJTv14P1Mzw3prAYqCLiGxzb8hK2UubsJSy+i3OZ3g9GjgB1BGRusBqIF+G6wkZvk4m\n8y9ifwLHgbbuCVUp76EJRCkwV7wuDBwRETHGVAPaOVDW6dT77zHGjHRVgEp5I00gSl09cupVYJAx\nZgtWbcShfgwRuQDcBTQzxozN6X6lfJXuSKiUUsopWgNRSinlFE0gSimlnKIJRCmllFM0gSillHKK\nJhCllFJO0QSilFLKKZpAlFJKOUUTiFJKKaf8P0nttc3JPRC7AAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%matplotlib inline\n", "import matplotlib.pyplot as plt\n", "\n", "plt.figure(figsize=(5,4.5))\n", "plt.loglog(xz,nlist)\n", "plt.loglog(xz,np.exp(intercept1)*xz**slope1,'--')\n", "plt.loglog(xz,np.exp(intercept2)*xz**slope2,'--')\n", "plt.ylim(1,1e5)\n", "plt.xlabel('rank')\n", "plt.ylabel('frequency')\n", "plt.legend(['data','slope=-1.03','slope=-1.47']);" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Shakespeare Heap\n", "\n", "Recall `swords[]` contains the tokenized words, consider reading it in blocks of 10000 words\n", "(`swords[m:m+10000]` for m from 0 until greater than len(swords) ).\n", "\n", "Appending each block of new words to a set will accumulate the distinct words (the 'vocabulary'), and the length of that set will give its size.\n", "\n", "Must grow much more slowly than linear, since starts at vocab size 2243 for first 10000 words, and makes it to only 26962 for the full 903810 words:" ] }, { "cell_type": "code", "execution_count": 12, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "(2243, 26962)" ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" } ], "source": [ "#recall len(swords) was 903810\n", "len(set(swords[:10000])),len(set(swords))" ] }, { "cell_type": "code", "execution_count": 13, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAqgAAAFACAYAAABwR44uAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAMTQAADE0B0s6tTgAAIABJREFUeJzs3XdcVfX/wPHXhyUOQMiBuEfL0rJlZaaWSs6fozQ1La0c\nWanV1xQVwZllw8xRmbnLzG0qmqWpmWXcvlbfHKTm3gouZL1/f5zLFRRS4MK9wPv5ePCQcz7nnvO+\ncvmcN591jIiglFJKKaWUu/BwdQBKKaWUUkqlpQmqUkoppZRyK5qgKqWUUkopt6IJqlJKKaWUciua\noCqllFJKKbeiCapSSimllHIrmqAqpZRSSim3ogmqUkoppZRyK5qgKqWUUkopt6IJqlJKKaWUciua\noCqllFJKKbeiCapSSimllHIrmqAqpZS6hjGmsjHmsDHmO2PMDFfHo5QqXLxcHYBSSim3tUJEero6\nCKVU4aMtqEoppTLzhDFmgzGms6sDUUoVLi5JUI0xUcaY34wxNnvld7d9f2ljzCpjzC5jzHZjTP00\nrylqjJlnjNltjNlhjGmfpswYYyYaY2Lsr+171fWG2st2G2NG5d07VUqpfOswcAsQCvQyxgS6OB6l\nVCHiqi7+p0QkDsAY0waYAdwNjAO2iEgzY8x9wGJjTBURSQbeAOJF5GZjTBVgqzHmOxE5A3QFbhOR\nGvZK1GYv+8sY8yjQEbgTSAE2G2M2i8iqvH3LSimVf4hIIpAIYIzZCFQHtrk0KKVUoeGSFtTU5NSu\nJJBs//4pYKr9mG3AIaCBvaxjmrJ9wHqgrb2sA/CpvewMMB/olKZstojEi0gCMD1NmVJKFXjGmAnG\nmL3GmBRjTO2rymoYYzYbY3YaY7YaY2637y9u/9cA9wEH8j5ypVRh5bIxqMaYmcaY/UAk0M0YEwR4\nicjxNIf9A1Syf1/Jvp1qnxPKlFKqMFgA1MOq/672MTBVRG4F3gZm2vc/bIz5BdgERInIsbwIVCml\nwIWz+EXkWQBjTFesSrErYFwVjz0WA4QA51wZh1LKJfyAwyIirg7E2URkEzjqOAdjTGngXqCJ/biF\nxpiPjDHVRGQtsPZ659Z6U6lCLdfqTZcvMyUis40xU+2bicaYMmlaUasA++3f/wNUBo6lKYuyf7/f\nXrY1g9ellpFB2dVCgIPZeBtKqYKhAtbQosKiInBERFLS7NuP1cu05wbPofWmUoVbrtSbeZ6gGmMC\ngGIicsS+3QY4JSKnjTELgD5ApDHmfqyKb4P9pV8DvYGfjTFVscam9rGXLQBeNMZ8jTWmtSPQIk3Z\nR8aYiViTpHoAwzMJ7xzAgQMH8Pf3d9ZbdpmwsDDGjBnj6jCcQt+Leyoo7yUuLo6KFSuCtgJmR4b1\nprt+NnIrrpyeN6uvv9Hjr3fcv5VntUx/5rn7enf7med2vemKFtQAYIExxhcQ4DjQ0l42CJhtjNkF\nXAa62GfwA7wDTDfGxABJQF8ROW0vm401iH83VhI6XkT+BBCRDcaY+cAf9ut9KSIr/y1Af3//ApGg\n+vj4FIj3Afpe3FVBei+F1AGgnDHGI00raiUy72XK1NX1prt+Nlq1apUrceX0vFl9/Y0ef73j/q08\nq2X6M8/d17vjzzw35XmCKiL7gbqZlB3HWnMvo7KLwNOZlKUAr9i/MiofBej6p0oplYaInDDGRGPN\nAZhpjHkSOCAiN9q9n++EhmZ4i3H5ebP6+hs9/nrH/Vt5dsvcjf7Mb7zcnX7m+iSpAiw/VSDXo+/F\nPRWk91KQGWOmGmMOAOWBKHsvVareWAvx7wQGAt2dcU39bBQ++jNXzmQK4ITVbDPG+AOxsbGx1zRj\nx8fHk5CQ4JrAlFvz8fHB19fX1WGoHIiLiyMgIAAg4Kp1mtV1/Fu9qW6M3l+Uu8voPpfb9abLZ/Hn\nB/Hx8VStWpWjR4+6OhTlhoKDg9m7d68mqapQCwsLo1WrVtqKlkV6f1H5wdX3uaioKJYvX56r19QW\n1DQyawlI/SuhoMzuV86TOotRW4/yN21BzT5tQc0Zvb8od5fZfU5bUN1IQZndr5RSyr3o/UWp9HSS\nlFJKKaWUciuaoCqllFJKKbeiCapSSimllHIrmqAWUhs2bMDDw4OUlJTrH5yJyMhI6tev78SolFL5\nVVhYGFFRUa4OQ7mY3lsKh6ioKMLCwnL1GpqgFmD79u2jc+fOlC9fHn9/fypUqEDLli05duwYAMaY\nHF/DGedQSuV/Y8aM0SWmCgm9t6jQ0FDGjBmTq9fQBLUAa968OQEBAfzvf/8jLi4Om81Gx44d9Rc/\nE7pQtlJKXZ/eW7JG7y3ZowmqE9lsEBlp/evq85w+fZodO3bQq1ev1HXKKF26NF27dqVMmTKO4xYv\nXsytt95KQEAATZs25ciRI46yyZMnU6tWLQICAggJCaFbt26cOnUq02uuXLmScuXKsXDhQgBSUlJ4\n9913qVmzJiVLluT+++/nu+++cxy/fft2GjZsSGBgIEFBQdx///3s3r0buNLFExYWRtmyZQkJCWHg\nwIEkJyc7Xn/48GG6dOlChQoVCA4OpkuXLpw8efKG40+9Rnh4OOXLl+eee+4BIDY2lj59+lClShVK\nly5Ny5Yt2bt3b7Z+DsrJEhNdHYFSec5Z9xZnnEvvLXpvySuaoDqJzQYtW0JEBLRoAZs2QVxc1r82\nbYLmza3ztGyZ/UokKCiI2rVr06tXL2bMmMHvv//O1Q9lEBGWLFnCr7/+ysGDB7l48WK6MSXlypVj\n8eLFxMbGsnXrVnbv3s2rr76a4fU++OADevfuzYoVK2jfvj0AI0aMYO7cuSxbtoyzZ88ydOhQWrdu\n7fiFfOmll2jcuDGnT5/m5MmTfPbZZ5QsWdJxzq1bt+Lj48PBgwdZv349Cxcu5O233wasv0gff/xx\nKlSoQExMDHv27MHLy4vOnTtnKf6tW7dSpEgR9u7dy7Zt2wBo06YNFy5c4L///S+HDx+mVq1atGzZ\nMl0FpvLY5cswcSJUrgz//a+ro1Eqzzjr3uKs+4veW/TekmdERL/sX4A/ILGxsZJWbGysZLQ/rYgI\nEXD+V2Rkppe8rtOnT8vw4cOlbt26UrRoUQkKCpI33nhDEhISZP369eLh4SEHDx50HD9p0iSpWbNm\npudbvHixlCpVKs17jpC6detKz5495a677kp3LhGRgIAAWbNmTbp9TZo0kdGjR4uISKNGjeTFF1+U\nv//++5prRURESHBwsKSkpDj2TZkyRWrUqCEiIgsXLpQKFSqke83BgwfFGCOHDh264fivPkd0dLQU\nKVJELly44NiXlJQkRYsWlc2bN19zzhv5bKgcSEwU+ewzkUqVrvxSPPmk0y+T+nME/MUN6qL89JVZ\nvaluzPXqkNy6t+Tk/qL3luvHn9N7izvJ7DOa2/WmPknKSVq3hk8+gcOHoVw5+OorqF076+fZvh06\ndIAjRyAkBFq1yn5MgYGBREREEBERQWJiIqtWraJbt274+fnRoEEDwPpLMFXx4sU5d+6cY3vRokW8\n++67xMTEcPnyZZKTk7l48SIi4hhrFBMTw2+//cbs2bMpX76847XHjx8nLi6Op556Cg8Pq6FeREhK\nSuLmm28GYObMmYwcOZLHHnuMlJQU2rdvz+jRoylWrBgAFStWTDemqWrVqhw4cMBx3aNHjxIUFOQo\nFxGKFi3K/v37CQkJuaH4K1eunO7/bPfu3SQmJlKhQoV05zXGOK6t8kBKCixYAOHhsGuXtc/HB3r3\nhlyeOaqUO3HWvQWcd3/Re4veW/KCJqhOUqcOrFgBy5dbv/R16mTvPI88At98k/PzXM3b25vWrVvT\nuHFjoqOjHZVIZg4dOkSHDh2YN28ebdu2xdvbmyVLltC+fft0v4Q1a9Zk2LBhdOzYkaSkJDp16gRA\nyZIlKVq0KCtWrOCRRx7J8BoVK1bkk08+AWDPnj20bt2aEiVKMHLkSAAOHDiQ7lp79+51/HIHBwdT\npUoVx7ii7MafWsGlCg4OpkiRIpw4cQJPT8/r/r8qJxOxfgGGDr3Sle/pCd27w7BhUKmSa+NTmQoL\nC6NVq1Y6k9/JnHVvgdy5v+i9pXDeW6Kioli+fHmuXkPHoDpRnTpWg09Of+mdcZ6zZ88yePBgfv/9\ndxISEkhJSWHdunV8//33jgpERDJ9/fnz5xERbrrpJry9vdm9ezdjx47N8NgmTZqwatUq+vXrx6RJ\nkwDw8fGhd+/eDBw4kB07dgBw6dIlNm7cSExMDGD9lXvo0CEASpQogZeXF97e3o7znj59mhEjRpCQ\nkMDOnTsZP348PXr0AKBdu3YkJiYSHh5OXFwcYP1l/dVXX2U5/rQeeeQR7rzzTnr37s2JEycAOHPm\nDIsWLSI+Pv66r1c58P33UK+ededMTU47dYK//oJPP9Xk1M3pMlO5x1n3FmecS+8tem8BXWZK5YCP\njw8nT56kQ4cOlC5dmlKlSjFgwAAGDx7MgAEDgH9fZ+7WW29l7NixdO3alYCAALp3707Xrl0zPb5u\n3bps2LCBcePGMWLECADGjx9Pp06deOqppwgMDKRatWq89dZbJNpnYn///fc88MAD+Pn5UadOHerV\nq8fAgQPTnTMhIYEKFSrQsGFD2rVrx5tvvglYlc6WLVvYv38/tWrVomTJkjzyyCNs3LgxW/Gn8vDw\nYO3atRQrVoy6desSEBBAnTp1WLJkiS6hklu2boXGjeGxx2DLFmtf69ZWkjpvHti77ZRSrqf3Fr23\n5BXzb3/pFDbGGH8gNjY2Fn9/f8f+uLg4AgICuHq/yj2RkZGsW7eOH374wdWh/Cv9bOTA9u1Wt/2y\nZVf2Pf44jBoFDz6Yp6Gk/hyBABGJy9OL53OZ1Zvqxmgdkrfyy73FnWT2Gc3telNbUJVSeWv3bujc\nGe6++0py+uCDsG4dfPtthsmpM9eBVEop5f40QVVK5Y0DB+DFF+H22+GLL6wJUbVrWzM2fvzR6uLP\nwBdfwMMP53xtYKWUUvmHJqjKLQ0fPly7YAqKY8egf3+oUQOmTYPkZGtc6ZdfXlmF/KoxWKktpu+9\nB926QeocgsOHrXxWKaWyQ+8t+YcuM6WUyh1nzsD48fDBB3DxorWvYkWrKbRbN/C6tvqx2WDyZJg/\nH9Ism+jg55eztYGVUkrlD5qgKqWc6/x5mDDBSk7PnrX2lS0LQ4ZAz55QpEi6w202ayhqSAi8/vq1\niamfn7WvRAn4+GPnrQ2slFLKfWmCqpRyjvh4mDoVxowB+zp/BAbCwIHwyitQvHi6w2026/CvvrqS\nx14tJMTKc3fvdu6DK5Tz6UL9ShUeebFQvyaoSqmcSUqCGTOsQaMHD1r7iheHAQOsJtGSJR2H2myw\nZAkkJlqNrKk9/6lSW0tLlYInn7QaXDUpzR/GjBmjyyQpVUiEhoby0EMPOR6gkBs0QVVKZY8IrFwJ\n//mP9bQnsLrvX3oJBg2CMmWAK1345cpZ+er58xmfTltLlVJKpdIEVSmVdTYbvPEGfPedte3pCT16\nWM9QtD/TOrULf+FCOHXq2lNoa6lSSqnM6DJTBdTRo0fp0aMHHh4evPDCC/z999+OssWLF+Pv70/3\n7t05fvy4U663bds2unTpkm5f/fr16dGjB6NGjaJDhw4EBAQQGRnJ8OHDue+++/joo48AmDNnDmvW\nrKFv375OiUXlogMH4Nln4d57rySnLVtaT4X65JN0yWnz5tautMmpn5/1b0iINeEpMhLWrIEpUzQ5\nVUopdYU+6jSNgvao0/j4eEqUKMHOnTupXr26Y//u3bvZunUrzzzzjFOuk5KSQuvWrTlz5gybN28G\nYNeuXfzyyy+OpHXo0KHExMTw5ZdfArB582bOnz/PhQsX2L9/P/3793dKLHktv342siwuDt56C95/\n/8qipPfcY/XJN2rkOMxms1pM162Dn35Kfwp37sLXR51mnz7qNGcKTR3iAnPmzKFMmTIsXbo0V8dK\n5rVGjRpRv359RowYkSfX00edKqfz9fWlbNmyxMTEpNu/aNEipyWnAEuWLKFbt24cOXLEsW/z5s10\n7tzZsf3999/TKE0iEx8fz2233caGDRto1qwZAOczG5yoXCcx0VqYtEYNGDvWSk4rVoTZs+GXXxzJ\nqc0GvXpBgwYwerRVZNVbVhd+796wYgV06mSNAnCn5FQpdePyS+/cokWLOHnyJE2bNi1QyWlhomNQ\nnSUpCY4ede45g4MzXMw8K6pWrcqePXsc2zNnzqRr166O7WPHjvHBBx9gjCFta3rqtjGGRx99lCee\neCLD8x84cIAqVarg5eWVLkHt3r274/sLFy7wyy+/MGPGDMe+xx9/HIAuXbqwZcsW9u3bh6+vLw0a\nNMjR+1VOImI9smngQNi509rn7w+DB0O/flC0qONQmw2aNoWTJ6+8PDkZnn7aajV1t9ZSpVT2BQcH\nM3nyZGbNmsXgwYPT9c7deeedTJ482am9cxEREZw5c8axb9euXfTu3Ttd75yHhwfDhw8HoGnTppw7\nd45Vq1bx0ksvAVbjR4kSJZwSU15LSEjAx8fH1WG4hojk6RdQBFgM7ABsQBRQzV62HtgDRNu/+qV5\nXVFgHrDb/tr2acoMMBGIAXYBfa+65lB72W5g1L/E5g9IbGyspBUbGysZ7U/nwAER67buvK8DBzK/\n3g3q2rWrvPHGGyIismfPHpk7d26Oz5nW2rVrRUTk+PHj4uHhIadOnbrmmNWrV0uFChWcel13cUOf\njfzml19EGjS48jn08hJ5+WWR48fTHRYdLfL88yKVKl370Q0Jscrzi9SfI+AveVwn5vevzOpNdWPy\nax0SEhIiq1evTrfvrbfecuo1Fi5cKPPnz5eqVas69k2fPl1SUlIc2w8//LBMnTrVsf3tt9/KP//8\nI1u3bpXPP/9cVq9eLevXr8/0Gg0bNpSXX35Z2rVrJ35+fnLzzTfLrFmzHOVnz56VF198USpUqCCl\nS5eW5s2by86dO0VEZPny5RISEuI4dsaMGWKMkWXLlomISFJSkvj7+8umTZtERCQ+Pl4GDx4s1atX\nl6CgIGnQoIHYbLYM43n66aclMDBQ+vTpk2HMw4YNc2xPmjRJ7rzzTvH395dy5cpJ165dHffirMaY\nkcw+o7ldb7qqBfVjEVkNYIzpC0wDHrO/0X4iktHqr28A8SJyszGmCrDVGPOdiJwBugK3iUgNY0wg\nYLOX/WWMeRToCNwJpACbjTGbRWRVbr9Jd1CtWjW2b98OwNy5cxk6dKjTzr1u3Tp+/PFHoqOjSUlJ\nwdPTkyNHjhAUFJTuuKu795Wb+ucfCAuDefOu7GvTxhp7euutjl0//WQ1pG7aZHUcFC0KN91kTYbS\nGfmFly7Un0tyo3cOctxD5+69c5UqVeKBBx64ofcyffp0FixYwIIFC4iKiqJNmzbUqFGDhx56iGee\neYaLFy8SHR2Nn58fgwcPpnHjxuzYsYOGDRty6tQp/vzzT+644w7WrFnDLbfcwpo1a2jVqhVbtmzB\nGMNDDz0EQK9evTh8+DAbN26kbNmyTJ06ldDQUHbv3p1ubOeMGTP4+uuv+eKLL4hPHfP/L8qVK8fi\nxYupUaMGBw4c4KmnnuKVV15h7ty5WY7xRhXIhfpF5DKwOs2un4DX02xnNi62I9DDfo59xpj1QFtg\nOtAB+NRedsYYMx/oBITby2aLSDyAMWa6vcy5CWpwsDXD2dnnzKGqVauyePFi5syZc804HkhfiVzt\n3yqRs2fP4uXlRXh4uGPfpEmTOHz4MHfccUe6Y9evX0/v3r1z/F5ULjl71hpfOmECXL5s7bv/fnj3\nXahfH7C68RcsgD17rCc/pbnfcOkSdOum3fmFnS7Un0uOHrXGfTvbgQOOVTeyo1q1ao4Ede/evXh7\nexMSEuIoL1u2LGPHjs32+Xfu3Enjxo05ceIECQkJnD59+prGj02bNlG2bFluvvnmbF8HoHnz5jRv\n3hyAZs2a0bZtW6ZPn07VqlX55ptv2L59O6VLlwbgnXfeYc6cOaxYsYIOHTpQt25d1qxZwx133MG3\n337Lxx9/zMCBAwFYu3YtjRo1wsPDg1OnTjFr1ix27txJuXLlAHjppZf44IMPWLFiRbo5G61bt3b8\noefr63vd+Nu2bev4vmLFigwaNIgXX3wRgBIlStxwjFlRWBbq7wcsSbM9zhgzAvgfECYie+37KwH/\npDlun31fZmV105RtvKqsoxPiTs/LK0e/7LmlWrVqxMTEkJCQQNWqVa8pz24lMmPGjGtm3pcpUybd\nX7pgzfL79ddfeeyxx7J8DZXLEhKurPWUuhZU5cpWi2mHDuDhgc1mTdxfuNB66pOHR/rkFKzEtFcv\nTUyVKkwKUu/c1ffGqlWrYrPZOHDgAMYYqlWr5ijz8vKicuXK7N+/H7DGvK5Zs4aGDRsSGBhImzZt\n6NOnD3v27GHt2rWOVuXUyWR169Z1nEtESExM5GDqE/gyied6Fi1axLvvvktMTAyXL18mOTmZixcv\nOhqZbjRGd+PSBNUYEwZUB3radz0jIofsZX2BFcAdmbxc3YDq1atz77330qNHD6ecb+PGjYwYMYK9\ne/fSoEED6tizkmHDhrF3714+//xzgoODueWWW5g9eza//fYbAFOnTuXuu++mQ4cOTolD5YCI9bzR\nN9+01nwC63GkQ4bAyy9j+8uXZSOtrvuICKuFNFVKii6wr1Seyo3eudTz5kBB6p3bt2/fNdsVKlSg\nYsWKiAh///03tWrVAiA5OZn9+/dTqZLVPtakSRPGjBnD0qVLHa2eTZo0Yf78+Wzbto1Zs2YB1uQy\nYwzbt2+nwnUas7LSmnno0CE6dOjAvHnzaNu2Ld7e3ixZsoT27ds7/p9vNEa3kxsDW2/kC2tM6c+A\n378ccwkItH//O/BAmrL5QA/79yuADmnKxgEj7N9/BAxMU9YHmJXJ9fwB6du3rwwYMEAGDBggq1ev\nzreD2FXuy3efjZ9+EqlXTxyzmby9Rfr3l/9+d1IiIkTmzRMpVUqumfSUdvLTvHkikZH5axJURlav\nXu34Pe/bt69Oksp+Xa6TpHIg39Uhdj/88IMUK1ZMPvvsM6ee9/33379m3z333CMzZ85Mty82Nla8\nvLzkn3/+ydH1GjZsKMWLF5eVK1dKcnKyrFq1Snx9fWXz5s0iItKiRQtp0qSJHDt2TC5evCivvfaa\nVKhQQc6fPy8iIikpKRIUFCSBgYGyYsUKERGZNWuWlCxZUqpVq5buWu3atZPWrVs7Yo6Li5NVq1bJ\n0aNH08WTdgJUZjGnHrNjxw7x8PCQb7/9VkREdu3aJQ888IB4eHhIcnJylmPMiKsmSbmqQnsN2Ia1\nuGvqPk+gTJrt9sDeNNvDgen276sCR4Eg+/azwFqs8atBWN34d9jLGtiT26JYKwj8AjTPJK7sz+JX\nhVK++Wzs2SPSsaOkyzaffFJk926JjrYSTxDx9Ex/iJ+f9W+pUiK9e+f/pDQzOos/R/W5Jqg5kG/q\nkKscOnRI6tev77Tz/fDDD9K4cWOpXr26RKepaIYOHSqBgYHSsGFDiYqKkr1798qIESOkXbt24uXl\nJYMHD5b58+dn+7oNGzaUV155xTGLv0aNGumS4TNnzsiLL74o5cuXl1KlSkmzZs1kx44d6c7x1FNP\nia+vr1y4cEFERI4ePSqenp7Su3fvdMfFx8dLZGSk3HLLLeLv7y/ly5eX9u3by5EjRxzHNGrU6LoJ\n6tXHjBs3TsqVKyf+/v5Sr149mThxYroENSsxZsRVCWqeP0nKGFMeOAD8DZzDWiIqHngc2AD42N/w\nCeA1Efnd/rpiWBOi7gOSgCEistBe5gFMAJpjzdSfICIfpbnmUKC7/bxfikiGg2UK2pOkVO5z+8/G\nmTPWyvkTJ1pjTgEefNCaAPXww0RHW08u/eOPKy8pUQLOn3fvJz85mz5JKvv0SVI54/Z1SAGX109l\nyo9c9SQpV8ziP0TmM/Xv/5fXXQSezqQsBXjF/pVR+ShgVNYiVSofS06G6dOt9aBSJ0BVqwbjxmGr\n1p6lywyxC6y8NTkZjLHaTAtTUqqUUsp9ucMsfqWUM23dCi+/DNu2WduBgRAezm8P9WHK9CLMnw+x\nselfIgING8J772lSqpQqPDKaxKXcgyaoShUUx49bLabTp1vbxlhT7EeNwnagFE88YR2SVuqM/JAQ\nTU6VUoXPd9995+oQVCY0QVUqv0tKgilTYNiwK02jdevCpEnYPO5l6hD49ttrk1PtzldKKeWuNEFV\nKj/74QerO//3363t0qXh7bex1erGu+97sGCBNTfK1xeCguD0aV2/VCmllPvTBDUL4uJ0cq9Kz2Wf\nicOHYeBAmDvX2vb05HiHl5leKYKyUpJ+jayu+1Tx8dZsfX0cqcotYWFhtGrVyrEQuFKq4IqKimL5\n8uW5eo08X2bKnWW2XEp8fDxVq1bl6NGjrgtOua3g4GD27t17Q89MzrGEBPjwQ+vxpOfPA3Dunkd5\nv+pEPlxf2zFh/2ohIbBihSammdFlprJPl5nKGV1mSrm7QrPMVH7k6+vL3r17SUhdR1KpNHx8fPIm\nOV27Fl59FXbsACDOL4Qf246n09KnORudfiaqPo5UKaVUfqYJ6g3y9fXNmyREqav98w+8/josXAhA\nsocXEzxeY/i5oZyf5XfN4Tr5SSmlVH6nCapS7io+Ht55h5QxY/GIvwTAGprwasqH7Ey5zXGYtpYq\nlf/pHAflrlz12dQEVSk39PcHywkc0Z+gM3vwAP6hEv35gCW0AUy69Uu1tVSp/MvHx4fg4GAqVqzo\n6lCUylRwcDA+Pj55ek1NUJVyI38ujSH51f7U3v8NAPEU4W0G8haDuEQxQJNSpQoSneOg8oM8m2uR\nhiaoSrmYzQarvr5A0+ix1Fr9DkWwblTLaMUA3ueEX3UuaRe+UgWWznFQ6lqaoCrlQrZoYdLjCwk/\n+xqVOABADNXpxwRW0oKQEPhYW0uVUkoVMpqgKuUify38H/HdXmXaxXUAXKQo7/gMZWzCa/iV8qW3\ntpYqpZQqpDRBVSoP2WzwzcJ4Ht0wkoc2vY03SQB8xVO8U2Y8r31QiTBtLVVKKVXIaYKqVB6w2eD9\n9+Hogo1h4DKsAAAgAElEQVRMjH+BW9kFwP+4nVeYSErDx/nkPU1KlVJKKdAEValcY7PBsmUgApPG\nxjEi4U36MBWAy/jwts8wRiYMpHSIDys0OVVKKaUcNEFVKhfYbBAaCidOQEuWY6MPFTgEwCbqMaT0\np/SecDtDtTtfKaWUuoYmqEo5kc0GS5bA6tXAieN8was8zXwAzlGCkcXHcf6Z3nzQy0OTUqWUUioT\nmqAq5SQ2GzzxBBw/LnRjNqvMAILkNACH6rRgQaMpdHqmoiamSiml1HVogqqUE/zyC3TsCEWP72M1\nvQhlDQhcKFaKE0M/pMqgp+lvjKvDVCrXhIWF0apVK0JDQ10dilIql0VFRbF8+fJcvYYRkVy9QH5i\njPEHYmNjY/H393d1OMrNpU6C8vCAUZHJ9EmeyGiGUJyLAJxu0ZWgGe9Zj4BSbi0uLo6AgACAABGJ\nc3U8+YnWm0oVTrldb2oLqlLZYLNB8+Zw9CjcwR9s4AUeZCsAR4tU4vz4j6nx8hMujlIppZRygpQU\nWLPGmv2bR72BHnlyFaUKkA0boE0bOH30MhEMJ5p7eJCtpGCYVvxVjq37U5NTpZRS+Z8ILF4Md98N\nzZrBypV5dmltQVXqBthsMGcO7Nxp/X4+KD+yiheoyV8AxJavycJm07j3pYe4SydBKaWUys9EYPly\niIiwboBACoY/p/9CrRYt8iQETVCV+hc2G0yeDPPmwcWLEOBxjgkSRl8m4YGQaLw52XMI5SYMokeR\nIq4OVymllMo+EVi1CoYPh23bHLu/4ikiGc6OJXcw5wvo1Cn3Q9EEValMpF1sH+AJVjE1pTeV2Q/A\nr94PUmzeNG5/8g4XRqmUUkrlkAisXQvh4bB1q2P3ItoSQQS/U9vakQKffKIJqlIuYbPBwoXWH5En\nTsBNnOQD+vMMcwFILFKcdY+NoeyIvtx+n6eLo1VKKaWySQS++85qMd282bF7Ga2IIAIb96Q73MMD\nevbMm9A0QVXKzmaDqVPhiy/g3Dnw9BBeKP4FYy70ozQnAYh7+An8503licqVXRytUkoplQMbNnDu\ntXD8on9w7FpJM4YTyTbuv+bwqlVh9Oi8aT0FTVCVAiA6Gh5/HM6etbYrsp8pKX1occGasZhU8ia8\nJn6Af5cuebbEhlJKKeV0mzdzbkA4fr98h599VxRNGU4kW3nwmsPzOjFNpQmqKtRsNvjyS+vLSk6F\n3kzlbQbix3nroE6d8JowAUqXdmWoSimlVPb99BNxA4bj/9MaR2K6jscYTiSbeSTdocbAY4/BO+/g\nssdza4KqCi2bDZo2hZMnwdsbbg86xvjT3WnOKgASgiviM20K5NGSGkoppZSz7Zi7jQuvD+feYytJ\nfdbbD9QnnBFsoGG6Y90hMU2lCaoqdFIfUfrTT1ZyCtA0cQXzLvbAH/uU/Z498Rk/Hvz8Mj+RUkop\n5WZS51P4/22jzX8jqHdymaNsMw8Tzgi+4zHgynC1KlXg3nthyBDXJ6ap8jxBNcYUAb4EbgcuAceB\nl0Tkb2NMaWAWUB2IB/qKyEb764oCnwH3A8nAEBFZaC8zwIdAMyAFmCAik9JccyjwHCDAfBEZmgdv\nVbkhmw1atoTDh63tMiUuMvz8G7zEFOsTV6oUfPYZtG7t0jiVUkqprPriCxj/7O+EJUbQnkWO/Vt5\ngHBGsIampE1MixSBzz/P+/GlN8JVLagfi8hqAGNMX2Aa0AgYB2wRkWbGmPuAxcaYKiKSDLwBxIvI\nzcaYKsBWY8x3InIG6ArcJiI1jDGBgM1e9pcx5lGgI3AnVvK62RizWURW5fF7Vm5g6dIryend2Fjp\n1Zly7LB2hIbCjBkQHOyy+JRSSqmsSG0xTfztT574OZJfWeAo28a9DCeSlTQnNTH19ob69eGWW6wl\no9ylxfRqHnl9QRG5nJqc2v0EpK7Z8xQw1X7cNuAQ0MBe1jFN2T5gPdDWXtYB+NRedgaYD3RKUzZb\nROJFJAGYnqZMFSK//GJNhvIghdcZz1bqUu7sDutPyAkTrGeYanKqlFIqH7DZoH17eLbuDhp80plp\nP9eigz05/Y27aM1S7ucXVtICLy9DlSrW8Vu3wrp1MGWK+yan4B5jUPsBS4wxQYCXiBxPU/YPUMn+\nfSX7dqp91ymrm6Zs41VlHZ0Qt8pHfvwRGjWC0gkH+d7rWR5N+s4qqFUL5s61/lVKKaXcnM0Go0bB\njuW7GZQ4gs7Mw5MUAH7nToYTyRLa4OnlwWOPun9LaWZcmqAaY8Kwxpv2BIq5MhZVMNls8NVX1hib\nVglf8wk9CUo6YxX27w9jx4Kvr2uDVEoppa4jNTH9fekeBiePZD6z8SIZgP9xOxFE8DVP4lPEg3Yt\n3WvCU3a4LEE1xrwBtAEeF5F4IN4Yk2SMKZOmFbUK2B98brWQVgaOpSmLsn+/3162NYPXpZaRQVmG\nwsLC8PHxASA0NJTQ0NCsvTnlFmw2aNYMLhw7x0TTj+f4HIBjHsGcmziTGi81dXGEypWioqKIirKq\nkISEBBdHo5RSGbPZYORI2L5sH28mj+ZLZuBNEgA7uYVIhrPEpyMPPeJJr3zaWpoRIyJ5f1FjXgM6\nYyWnsWn2Twf+EZFIY8z9wCKgiogkG2OGA5VFpIcxpiqwBagpIqeNMc8CzwChQEkgGmghIn8aYxoA\nHwEPYE2S2gQMF5GVGcTlD8TGxsbi7+9/dbHKR379FZ5+GoJitjKXLtTgbwB23PZ/JEyaRu3HSrk4\nQuVO4uLiCAgIAAgQkThXx5OfaL2pVO6ZNw+GPnuAgUmj6cF0fEgEIIbqRDKcRT6daNbKyyWtpbld\nb2ZrkpQxppW9UsIY84Yx5mtjzJ03+NrywHggAPjeGGMzxmyxFw8CHjbG7MKazNTFPoMf4B2gmDEm\nBliFtQTVaXvZbGAHsBurFXW8iPwJICIbsCZN/QH8CURllJyqgmPtWnjkwSSejhnJZupRg7+5aIqx\nf8jH3Pa/xZqcqnwlJ/WtUir/sdmgVy94ptEhznR5mb+SatCbj/Ehkb1UoQefcbfPX1xq35VNP3nx\n9dcFo8X0atlqQTXGbBeR2saYu4CZwBTgGRGp7+wA85K2BORvNhtMmwbrP9/Lx5e68gibAdhZ4l6S\nZ82lZttbXRyhclfu3ILq7vWt1ptKOc8XX0BYj6P0i3+L3kzFl8sA7KcioxjKFz7PEdrKxy3Gl7pl\nCyrYBz9AU+ATEfkYKO6ckJTKOpsNHn9MiJ08hy2X7uIRNpOCYWKJQVxa96Mmpyo/yxf1bVhYmGNM\nr1Iqa2w26N7iOIc7v8Gf8dXozwR8ucxBytOHydT23c3p9j354Scft2gxjYqKIiwsLFevkd0W1P9i\nzbyfAHS3L4j/u4jk67V6tCUgf9q4EZ5vf5bIE33oxJcAnC5RkeVPzqL2qw1d/ous3J+bt6C6dX2r\n9aZS2ZO6wP6Ff05y97rx9EmaSHEuAnCEYMYQxoKAF2nbydctJz7ldr2Z3Vn8Q4GPgXX2yvJWYJfz\nwlLq39ls1tJR+/fD2WU/sFa6Utm+OMNS345UXTGFZxsEujhKpZxC61ulCgibDZYtA2Ng0sjTvJL0\nHv2YgB/nAThGGd5iENM8exPapiir3KAr31VcMovfXWlLgHtL/cUuUgQiIiD5ciIjTARvylg8EOLw\n46PbJtFs7jPUucdc93xKpXLnFlR3p/WmUv8u9d7l7Q2jR4P3xbP05wMG8D4BWNXNCUoxjjeZQh8o\nVpxp06CTmz/z0q1aUI0x4f9WLiIjchaOUtdK7Qb5+ms4bV+34WZ2MZcu3C/bAPiRh3i9zBwmz6tW\naP/aVAWL1rdK5V+pSamPj5WUXrgAfsTxGhN4jfcI5CwApwjiHf7DnICXeahJCbqVKjjrmOZUVrv4\n/ez/VgAeB5YBArQG1jkxLqUA65e8RQs4ciR1j/AC0/iA/hTnIkl48lf7Yay7YwiT23jpL7UqSLS+\nVSofSZuUjhkD561ee0pwjkF8xBuM5yasVpYzlORdXudT31ep38Kf5YW4Kz8zWUpQReQ/AMaYNcDd\nInLYvh0OzHB6dKrQSv1F//nnK8lpEKeYxgu0ZQkApwOrceK9OdR67iHcYraIUk6k9a1S+ccXX1gt\nn6lJKUAxLvASkxnI25TmJACx+DPJewC+g/rj41WS1a00Mc1MdidJhaRWlgAicsS+AL9SOZJRd36J\nElDt/H9Z5tGGyin7rJ3PPkvQxIkE+fllei6lCgitb5VyUzYbTJkCM2dC6hOTi3KR3kzlTcZRFuvJ\n7ecowfZG/dn8wACadQzSpPQGZDdBPWiMiQSm2befBw46JyRVWNls0Lw5HD2afv/7db/k2R964J14\nCfz84NNPoWNH1wSpVN7T+lYpN2OzwcSJMH8+XLRWhqII8fTkEwYzlnJYN7IEn+L8/MAr+Ee+Tr3H\nSlHPhTHnN9lNUJ8DPgR+wxoT9a19n1JZZrPBokXw7bfpk1NPkviweBgvrHvH2nHLLbBkCdx+u2sC\nVco1nkPrW6Xcgs0GH3xg9fKlJqY+XOYFphHGGMpjdXYkehfl9NN9KTv+PzxSpowLI86/spygGmM8\ngXoi0iEX4lGFjM0GjRtb3fkeHhAQALGxUCPoNEuLPU3Ng2utA1u0gLlzrQOUKiS0vlXKPaQOP5s7\n15qRD+BNAs8xg6GMohIHAEgp4otHn954v/kmZYODXRhx/pflBFVEko0xQ4CFuRCPKkSio+HJJ6+M\nNU1Jgaefhrs9ttN9WRuKHNxrFQwbZi186pHdJ/MqlT9pfauUa2U0L8KLRLoxi2GMpAr/AJDi7YNH\n7154DBoEISEujLjgyG4Xf7Qx5hER2eTUaFSBd/Vi+5cvW3lnSor1O/1m1a+oOqK71XdSogTMmgVt\n27o6bKVcSetbpVzg6mUOPUmiC3MJZwTV2QNAipc3Hi88j0dYGFSs6MJoC57sJqgPAs8ZY/YAjkUV\nROQep0SlCiSbDUJD4cSJ9PtTUuCxBsnMqzqEsoPGWTtvvtkab1qzZt4HqpR70fpWqTxms8Err1jJ\nqQfJPM2XDCeSW9gNgHh5Ybp3x2PIEKhc2cXRFkzZTVD7OjUKVeCtXw9duqRPTv384Nw5qBl8msUJ\nnfGfEWUVNG9uDfQpWdIlsSrlZrS+VSqPpHbpf/UVxJ5NoQMLiCCC29kBgHh6Yrp1wwwdCtWquTja\ngi1bCaqIbAAwxoTYtw//+ytUYWSzWY2gBw/C55+DSPru/PHjIXbT79Z40y1WdwlDhkBkJHh6ujZ4\npdyE1rdK5b60Y03PnE6hHYuIIII7+RMA8fDAdOmCGTbM6uFTuS5bCaox5nbgayDEvn0QeEpEdjgx\nNpWPpe3ON8ZKTsFKThs2hPfegzoxC2Bmd2tKZPHi1njTdu1cGrdS7kbrW6Vyl80GzZrBsWPC/7GU\nSIZzF9sBSMFwtlkngt4Ph1tvdXGkhUt2p0VPBkaLSKCIBAKjgSnOC0vlVzab1QA6dOiV7nwRqzsf\nrJbT995Jps5Xg6FDBys5rVEDtm7V5FSpjGl9q1QuWbcOWrUU7ju2gm3cxxLaOpLTbdU7sPPrPwha\nOVeTUxfI7hjUQBGZl7ohIl8aYwY5KSaVT9ls0LIlHLZ3QKaOMU3tzt+9G9o0OEPtYZ1h9WrroGbN\nrPGmgYGuC1wp96b1rVJOZLPBl1/Cjr+EhOWrWcRwHuCXK+VV2+E7djj3daztwihVdhPUZGNMTRH5\nH4AxpiaQ7LywVH5js0H//leSU4DOna3ktFUrrOcO//EHtGkDf/9tHTB4MIwcqeNNlfp3Wt8q5QQ2\nG0yeDPPmCg9f+pYRhPMQPznKN9/UmpsmRlCnUx0XRqlSZTdBDQN+MMZst2/XAro4JySV39hs8MQT\ncPz4lX0hIdCrlz0xBVi4EJ599sp40xkzrFX6lVLXo/WtUjlks0GTJlDr1PesJpz6XFlW+Bua81Gp\nSMasuY/bNDd1G9mdxR9lH7hf177rJxE56bywVH7y0Ufpk1PHJKg6QHIyhIfDmDFWYfXq1tT+O+90\nRahK5Tta3yqVfTYbLFgAB+ZtZMGpcBqx3lH2fZFQEodE8qvUZUyrNA0qyi1kKUE1xkQC3wFbROQE\nsCJXolL5gs0GY8da+WZAAMTG2idBpSanZ89a/fyrVlkvCA2FL77Q8aZK3QCtb5XKGZsNwhr+yGtx\nw2nCt479G7wf5+fmkTQeXo86daCpC2NUmctqC2oZ4GOgojFmC/A9VgX6s4jomKhCJDoaGjWCuDgI\nCrJaUXfvTjPe9M8/rfGmMTHWCwYNglGjdLypUjdO61ulssFmg83v/0zthcNZdXG1Y/+uco+y5YkR\n1H6lAf/R1lK3l6UEVUT6ABhjygEN7V8zgWBjzEYRaeHsAJX7+fVX+L//s5JTgNOnreQ0PNx+wKJF\n1njT8+ehWDFrvOlTT7kqXKXyJa1vlcq6HXN/5Xj34byc+I1j3ybq8eFNIxi8ohHP3mNcGJ3KiuyO\nQT1ijFkIHAGOAp2Au50ZmHJP0dFQvz5cupT+qVCtWmFthIfD6NHWwdWqweLFUFuX6lAqu7S+Ver6\ndnz5G0nDIrgzZim32ff9RF22Nh9B7ANNGNza6BjTfCarY1AfxforvhFQHvgJ+AFoISK7nR6dcjth\nYVZyClc9FarqWWjVBVautAqbNrXGmwYFuSxWpfIzrW+VugF//MGZ/hHctm6hY5fN8z6GJEfy33LN\nWDFKE9P8KqstqOuxKskRIrL6OseqAsRmg0mTrKdulCoFJ0+mmRAVuA/qhsKuXdbBAwdas/Z1vKlS\nObEerW+Vythff3GmfyQBa74iEOtZ2jbuJpwRlH++JQ+WN4zWmfn5WlYT1NS/6N8wxnwE/IxVia4X\nkV3ODU25C5sNmjeHo0fB3x8+/DDNhKiiO+CRxnDokDXedPp06NjR1SErVRBofavU1XbtghEjkHnz\nCBQrMf3d1GJ8iUhmnWtDSIhhRG9NTAuCrE6S2gRsAkYZY3yw1uVrBCwzxpQQkQq5EKNysUWLrOQU\nrIlRjglRNhvUb2o1pwYGWo8vfeABl8aqVEHhLvWtMaY/1rCCJnlxPaUy9PffMHIkMns2JiUFA/xJ\nTSKIYKG0p2dnDyLTPrlQ5XvZmiRljAnBqigbAo9hLYey6d9eo/Kn6GiYPx+8vCApKc2EqM2boUUL\na/HT4GBYu1YX31cqF7iyvjXGeAF3gb0PVam8tm+ftUThjBmQnIwBdnAr43wjWFXiKY6d9Lz2yYWq\nQMjqJKlPgQZACJC6Lt8zwC8ikuT88JQr2WzWJKhz56y5Th06QM+eUOfEGmjbFi5ehMqV4dtvoUYN\nV4erVIHiJvVtV2Ae8J88up5Slv37YfRoZPp0TJL1cY+hBhEM5ws6kRLvSa9nrzSaaHJa8GS1BfUA\n8DzWo/YScyEe5UamTLGSU7DWOi1XDursWwxPPw0JCXDrrVbLacWKrg1UqYLJpfWtMcYATUWkkzFm\nYF5fXxVShw5xYsAYghZ+imdKIgbYQ1XG+YRz/4Rn+H6kFymH0VbTQsAjKweLyAgR2ZjTytIYM8EY\ns9cYk2KMqZ1m/3pjzB5jTLT9q1+asqLGmHnGmN3GmB3GmPZpyowxZqIxJsYYs8sY0/eq6w21l+02\nxozKSeyFxeXLVsOon5+1HRICz5pZ1oL7CQlw993www+anCqVS5xV3+ZAO2CZi66tCpsjR6BfP1Kq\nVaf0gsl4piSyj8o8zzRuZSefJDzH4eNerFgBkZGwYoUmpwVdtsagOsECYBzXjqMSoJ+ILM/gNW8A\n8SJyszGmCrDVGPOdiJzB6oa6TURqGGMCAZu97C/7WoIdgTuBFGCzMWaziKzKnbeW/9ls0L8/+PhY\ny0qtWgXdL06iYvjL1gEPPwzffAMlS7o2UKXUDTHGTABaA5WBu0Vke5qyGlhPqCoFnAWeE5G/gFuB\nBsaYbsBdxpgeIjI976NXBdqxY/D22zB5MsTH4wEcoAKjGMrndMfXz4ekc+m78jUxLRyy1ILqLCKy\nSUQOAxk9cyyzmDoCU+2v34e13Epbe1kH4FN72RlgPtbTVlLLZotIvIgkANPTlKmr2GzWGvs//GB1\n63t5Qbj3WCqOsyenTZrAmjWanCqVvywA6gH7Mij7GJgqIrcCb2Mlq4jIGBEJFZFmwG+anCqnOnkS\n3nzTeuLge+9BfDyHCWHCzR/xcOkYPqEXpUN8+PhjbTEtrLI7i7922r/AM9uXTeOMMSOA/wFhIrLX\nvr8S8E+a4/bZ92VWVjdN2carynShzkzMnm3VGwAnTggXXx0Mm8ZZO9q2tZ4OVaSI6wJUqpBxRn1r\nX7IqdVxp2vOUBu4FmtiPW2iM+cgYU01E9qR5fdOcvAelHE6dgnfftRbUvnABgKOUZSyDWXRTT5bN\nL8qjwPLlOvmpsMtuC+qMG9yXVc+IyG0ichdW9/8KJ5xT3aANG6x19n19wZDCzGIvUS81Oe3aFb76\nSpNTpfLejBvclx0VgSMikpJm336u/PGvlHOcOQPh4SRXrgpjx8KFC5wwpRkVOJ5q7OFD+nHwVFGW\nL7eS0vBwTU4Lu6wuM1UGCAaKGmNqcaWLPgAontNgRORQmu8nGWPGG2MC7d32/2CNnzpmP6QKEGX/\nfr+9bGuasv1XlZFBWYbCwsLw8fEBIDQ0lNDQ0Oy9oXxkyxaraz8hAYJvSmRphed4IGaeVfjSSzBx\nIni4ZESIUrkiKiqKqCirCklISHBxNNfK7frW2QpjvaluQGwsTJhgdePHxuIJnOQm3vMcyO0f9aVF\n3eJMaQmXDqdZZ1u5rTytN0Xkhr+AfsBe4LL939Sv34Dns3Iu+/n2ArXt33sCZdKUtQf2ptkeDky3\nf18VOAoE2befBdZitQgHYXXj32EvawD8DhQFigC/AM0ziccfkNjYWCksoqNFwsJEypUTAZEiXJIl\ntLY2QGTwYJGUFFeHqVSuio2NFaxJmv6SxXost76cXd/KVXWufbs01sQojzT7jgDVsnDOQldvqhsQ\nFycyerRIYKDjfnLaBMpgRksJ4gREIiOtQ6Ojre+jo10bssqa3K43s/qo0wnABGPMMBEZmZXXpmWM\nmQq0AMoCUcaYc1hPK/nG/kg/AU5gzTpN9Q4w3RgTAyQBfUXktL1sNnAfsBtrpv54EfnTHvMGY8x8\n4A/7eb8UkZXZjb0gsdmgeXPrMaZeXlD5pvNMO9WGxqyzDnjrLWsQu1Iqzzmrvr3ONU4YY6KxVkKZ\naYx5EjggacafKpUlFy7ApEnWzPxTpwCI8whggudrJL3Uj5kLAjh/VWupzsxXGTEiWX+CnTGmFbBR\nRM7atwOBh0XkGyfHl6eMMf5AbGxsLP7+/q4OJ9cNHQqjR1vfl+QMv5ZpTrXjPyHGYCZNgj59XBug\nUnkkLi6OgIAAgAARiXN1PGk5o769qlHgFHBORG6xl92CNab1JiAW6J76B/4NnrtQ1ZsqExcvWk93\nGTcOTpywdnn5MT6pP+8zAN/gQFbam4Z0AlTBkNv1ZnYT1N9E5O402wb4VUTucWZwea2wVLQ2Gyxa\nBEuWwK5dUDLhGN97NaVm0nbw9LSeefzMM64OU6k84+YJqlvXt4Wl3lSZiI+Hjz8mceRYvE9ZU0Qu\neRZngvRj+S2v8+OOIMehkZHW5CdVMOR2vemUhfpFRIwxns44l8pdNhu0aGE9tMPbG+aM2c9jbzWh\n1Kld1sr8X30F//d/rg5TKZUJd61vw8LCaNWqlU6OKiwuX4Zp02DMGDh8GG/gAsWYxMv80+4N3hxf\nmtBT0LIlHNYJUAVOVFQUy5dn9Ewl58lugnrOGPOwiPwIYIypB5xzXlgqtyxebCWnAFUSdxE6qjEB\nsQegWDFYuhQaN3ZtgEqpq+WL+nbMmDHagloYJCRY6xGOHg0HDwJwCV8m8xJvM5DjlCXyTqhUyfpa\nsUK79Aui0NBQHnroISZNmpRr18hugjoQWGyM2WHfvpkrT3VSburXX2HePKsXv2bydtZ5NCEg9rj1\nVKiVK+Ghh1wdolLqWlrfKtdLTISZM2HUKPjHei7OZVOEmUV6caDLID5bWY7jR65tKdUJUCq7spWg\nisgWY8ztQGpG82PqAH7lnqKj4dFHrXHsTf1/YnF8M4olnIUyZaxHl951l6tDVEplQOtb5VJJSTBn\nDowcCXusxR0S8OYTejKx2GCmR5WnZz1oZ9OWUuVc2R6DKtbi+bpcUz4RHm4lp434joVxrSnGBahQ\nAb79Fm691dXhKaX+hda3Ks8lJ1uPto6MhJgYAJKMF597PM/I5DAOUAkuwLp1UK+etpQq58vWo4GM\nMU8YY3YYYxKMMcnGmBRjTLKzg1M5Z7PB889DVBR08VvGSppTggvEV6wBmzZpcqqUm9P6VuWp5GT4\n8kvia9xpPeI6JoYkPJlT5Hk+/c9uam6YSnKI9SRcnfikclN2W1A/BF4BtgBaUbopmw2aNYNjx+D5\novP45EI3PEjmUo1aFN24BoKDXR2iUur6tL5VuS8lxVp/MCIC/vwTXyAZD2bTlf3dhvHapOqUKGEd\nqhOfVF7IboIaJyJRTo1EOd3SpVZy2pOPmXKpDx4I1K1L0ZUrISjo+idQSrmDfFHf6jJT+ZSItSh2\nRARs3w5ACoZ5dGYE4ezmFiKr40hOQbvzVd4sM5XdhfojgN9EZInTI3Khgrbg9JtvwuV3PuQD6QfA\nufsa4ffdUvDzc3FkSrkXN1+oPwI3rm8LWr1ZaIhYTaHDh1vdbViJ6SKvjvzeNpxPN93OEfus/BUr\nNCFV13LXhfr7AQHGmEvAZcBgrR+tzXJuYuFC+Gv8NyyR/gDE/n979x0fVZX+cfzz0IsUQVHwByIq\nFlQWWFFw7QUbiItSRGxY1ro2LLEgqOiq64plBV0RZQVFXaS4GgsgYkGQQURBYKkqVnonyfP741zI\niAGSkMnMZL7v12temTnn3rnPZcLNM+eeckx7amUPhypVkhyZiBSRrrdSctzh7bfDyNkpU7YU/6f8\nOUnGVDUAACAASURBVMzv0ZuLHjmEc+pCR43KlyQrboL6hx1vIskyaRL07fIVH+Z1oxzOqlbHUuvd\n18JKUSKSbnS9lZ3nHmZtuftu+PTTLcUj6MjfqtxD35HNuemU/M11G1+SrbjzoC4s6UCk5PS++hdG\n5LanJquYxz68fuJr9FJyKpKWdL2VnTZuXEhMJ07cUvSmncndfg9TaQXrQ856yinbeQ+RUlakaabM\nbFj0M2ZmU7d+JCZEKaxYDC7oupHbp55DE+azil3ouftoTuq6W7JDE5Ei0vVWdtqECXDccXDCCVuS\n0w+qn0a3JpP4bsBofmjQCtB0UZKaitqC+kj08/qSDkR2TiwGp53q9P3pGo7lA9yMMd2G8ejNzXSb\nRiQ96XorxfPxx6HF9P33txRN3e1kbt/Qhw4PtGHIFVChAhx+uPqZSuoqUoLq7p9HT/d190HxdWZ2\nCfBBSQUmRTNyJHT+6Qku51kA3jv5Ibq9dGaSoxKR4kq3662mmUoBkyaFUfnZ+bOSfV7reHqt6cOB\n5x7N0Huhbt38zdXPVIorlaeZmuruLbcq+9zdW5VYZEmQztOlDD7vHXoMO43y5DG86oXsP/F5WrS0\nZIclkhZSfJqplL7epvN1s8z4/POQmL755paixY2P5qKFfRjrx7P77iFnVTIqJSmlppkys9ZAG2B3\nM7surqoWULkkA5PCG/f0LDoO60x58ljcsC1NXx3IH5SciqQ1XW9lh6ZNC4npqFFbita1aEOfcn15\ncuaJrPHwd+Dnn8OtfCWokk6KNEgKqE+Y8qQa0CLusRtwUYlGJoUybexSGl7Vntqs4NvyjVg+6D/8\n4Qj97RIpA3S9lYJ9+SV06hQyzig5zTu8NUN7vMVusz5i1REn8d+3jAYNwuYaBCXpqKh9UEcCI83s\nNHd/K0ExSWFt2kTuOZ3Zj7msoRpn5I6i08d7cOhJyQ5MRHaWrrfyO19/DX36wPDh+WUtW/JFp750\nHnw61b40xo2H1q1D1ZgxGgQl6auoLaibNTKzWgBm9qSZTTGzY0owLtmBWAzeb34DrZaFUZrn829+\nadBc35JFyh5dbzPdN99A9+5wyCH5yWnz5iwf/AY9DprC0Q+ewVVXG5Mn5yenEJLSu+9WcirpqbgJ\n6tXuvsLMjgIOBe4gf0oUSbBYDF4+7mlOnPkUAF90vo8Wfc7WeskiZZOut5lq7ly48EI4+GAYOjSs\nBtWsGfMeeo0O/zeVhtecxbr1xsyZ8Ne/hqmjRMqK4v4650Q/TwBedPdsM3ughGKSHZj+2FjuX3kt\nAEPpxtyDs7j77iQHJSKJouttppk3D+67D158EXJzQ9mBB8I99zAs51wuvLgcmzZBnTpwxx2w117J\nDVckEYrbgppnZl2ALsB7UZnW0iwNc+dy7vBzqEAun3E4d9V/jvYdNGJfpAxLi+ttVlYW2XHzb0ox\nLFwIl18OBxwAzz8fktP994eXXmL5xBlcO7ELF1wUklOApUtDH1OR0padnU1WVlZCj1HcBPUaoBvw\nrLsvNLOmwNiSC0sKtGIFK45pT7X1y1hZYy8m3vQGr71ZVbf1Rcq2tLje9uvXT5P0F9e338JVV4Vk\n9NlnIScHmjSBF17Av/qaIbnnccDB5Zk1K3RB1eh8SbZ27drRr1+/hB6jWLf43f1ToGPc69nAddve\nQ3Zabi6/ntKVuktmsZaqdK48kge6N1ByKlLG6Xpbhn3/PTz4IAwcCBs3hrLGjeGuu6BHD2Z8U5Gr\nTwpdUZ94As49F8zCJhqdL2VdUSfqv8nd/25mjxZU7+43lkxYsrW8m3pR97O3AbiIwWT/0oq2mnhZ\npMzS9bYM++EH+NvfYMAAWL8+lDVsCHfeCRddxKoNlehzOzz1FFx9dZguqkaN/N21RKlkgqK2oK6O\nfq4o6UBk2xb2fo69+/8DgEd3uZtXV3fWrR2Rsk/X27Lm55/hoYdC5rluXSjbay/IyoKePfFKlXn1\nVbjxxnCHf/LkMLOUSCYyd092DCkjFdeUnv3chzS+9EQqsYnRlTux+rnhzPlfOd3aESlBiV5TuixL\nxetmyvn1V3jkkXCffs2aULbnniExvewyYjOr8Pzz8NlnMH9+2PT888PtfJFUlejrZpH7oJrZ4cDN\nQLOoaAbwiLtPKcnABJg/n/rX/JlKbGIqLei64QVu/V85TSklkiF0vU1zy5bBo4/CY4/B6qhBvF49\nuO02+MtfoGpVPvoITj01VFerBm+9BcdoGQaRoo3iN7M2wDvAPOBO4C5gPvCOmR1R8uFlri8mrmJJ\n6w7UWP8LP9oenMVIajeortv6IhlC19s0tmJFWJK0ceMwn+nq1VC3bri9P28e3HADmypU5emnoV27\n/Nx17VoYPz6ZgYukjqK2oN4CXOLuI+LKRpjZp8DtxI00leKLTcnl55O603zDDNZTmWm93+Aya6jb\n+iKZRdfbdLNqFTz+eLhHv3x5KNt1V+jVC665BmrUwB1GjYRbbw0LQ/XtC3//exjQr7EFIvmKOg9q\ns60ulgC4+0jg4MK+iZn1N7P5ZpZnZofFle9uZm+Z2Wwzm25mR8fVVTWzoWY2x8xmmVmnuDozsyfM\nbG6079VbHe/OqG6Omd1XxHMudWtuuJNTNoTZly/lX0yyI7WeskjmKZHrrZSC1avDdFGNG4eR+MuX\nQ+3acO+9sGAB3H471KjBpElw7LFw2WVw3XUwY0YYEDVmTGhw1XLVIvmK2oK6djt1a4rwPq8CfwMm\nblX+IPCJu59mZn8ktBY0dvdcQj+s9e6+v5k1BiaZ2Vh3Xwb0AA509/3MbFcgFtXNNLNjCCuwHALk\nAR+Z2Ufu/lYR4i09Q4bwp4kPAvAAtzGuwfncpG/UIpmopK63pSIrK4v27dtn1mT9a9fCP/8Zpoz6\n5ZdQVrMm3HADXH99SFKB//0vjIcaMyY/IY0fT6ZpoyTdZGdnMzrBy5gVNUGtbGaHAgWNLaxS2Ddx\n94kQWj63quoM7BttM8XMvgOOJaya0gW4JKpbYGbjgbOBQdF+z0Z1y8zsFcLKK3dHdUPcfX10zEFR\nXeolqJ9+Sl7PSykHTGvUnk0X38+Ys3ThEslQJXK9LS39+vXLnFH869aFyfUffBB+/DGU7bIL/PWv\nIQOtUwcIg/fvvTds2r07zJ4dZpUSSXft2rWjTZs2PPXUUwk7RlET1KrAqG3U7dR8VWZWB6jg7j/F\nFS8EGkXPG0WvN1uwg7oj4uo+3Kquy87EmhCLF7PpjI5U3LSR6RxK540v8cpZ5ZScimSuhF1vpZjW\nrw9LkT7wACxZEsqqVYNrr4Wbb4bddtuy2eOPQ79+0LZtmD7q0EOTGLdIGipSgurujRMUR2ZbswY6\ndKDi0h/5md3owCgW/lCD0VopSiRj6XqbQjZsgEGD4P774bvvQlmVKmGZp1tuCVNHAXl5MHQo3HFH\nGLT/+utw4olJjFskjRV5HtREcfelZpZjZvXiWlEbA4ui5wuBvYEf4+qyo+eLorpJBey3uY4C6gqU\nlZVFpUqVgNCMndA+VXl5LDvrQnadNo2NVKRn7f+wcHljjeYUSbDs7Gyys8MlZOPmddBF4m3aBIMH\nh6miFkV/NipXDnOY3nor1K8PQCwWpjqdNCnc/b//fjjvPChX1GHIIrJFUleSMrP5wFnuPj16PQhY\n6O59ogmq/wM0dvdcM+sN7O3ul5jZPsAnwMFRYnshcD7QDqgNTAXOcPevzOxY4EmgNWGQ1ESgt7v/\nt4B4Sn1FlCWX96b+s30BuLb6c7R99hLmzEFTSomUIq0kVXxlciWpnBwYMiR0IJ0/P5RVqgSXXx4m\n2Y/rSDpmDJxzTmhkrVkT3nkHjtAstZIBUm4lqZJgZgOAM4A9gGwzW+XuTYHbgCFmNhvYAHSPRvAD\nPAwMMrO5QA5wtbsvjeqGAH8E5hCS0Efc/SsAd/8gGjQ1g9Bv6+WCktOkeOWVLcnp37mRJ9dcwu5z\n0EpRIiLJkJMDw4aFyUnnzg1lFStCz55hGH7Dhls23bABHn44TA+VkxPKVq6E7GwlqCIlIaktqKmm\nVFsCZs6Eli1h/Xre5lTOYAx7NiivefBEkkAtqMVXJlpQc3Nh+PCQbX7zTSgrXx4uvjh0KG3c+Deb\nv/NOmHd/88D9rKz8ifZ1DZdMUSZbUIUwvHP9euba/rzX82V6Nyyv2/oiIqUpLw9eew3uuSc0GkDo\nOHrBBXDXXdCkyW82//bbMItUdnboZ3rllSGPPewwGD1aXbNESpIS1GRYtAh/+WUM6O29Gf/fWvrW\nLSJSWvLy4I03oHfvsJwThMS0e/eQmO6//28237QJ+vcPDawdO4ZG1j33zK/XRPsiJU8JajI89hiW\nk8NCGjGczuR8j6aUEhFJNPdwse3dG6ZNC2Vm0LVr6Px/4IG/22XCBLjqqvxdjzuudEMWyVSaBKOU\nTZ+wnPVPPgvAwKo3kENFTSklIpJI7vDmm3D44XDWWfnJ6bnnwpdfhslLt0pOf/wx3Ok//fTwMxZT\ncipSmpSglqJYDEafMYAqm1azjNq0eOpS+vRRp3oRSX9ZWVlb5pVNGe6hw+iRR8KZZ8Lnn4fys8+G\nL74IA6OaNfvNLrm58NRTcMABsHZt6Jp6yy1hlikRCbKzs8nKykroMTSKP06iR6Ped9cGet7XmPr8\nwP1kkdvnfk0pJZICNIq/+FJyFL87jB0bbtt//HF+efv2YUBUy5YF7jZpUridv2IFPPkknHpq6YQr\nko4Sfd1UC2opOmr+v6nPD2ygEq/uca1u64uIlLQPPgj34k86KT85Pf10+OwzGDWqwOT011/hiivg\nhBNCD4AZM5SciiSbEtRSMnFCHg2GPQLAVy0v4Pm39tRtfRGRkjJxYlj4/rjjwsgmgFNOgU8+ye9/\nupW8PHjuuXA7/9tvYfr00OhapUrphi4iv6dR/KUgFoPHTn6T1/JmkYdR5c6bOVjJqYjIzvv00zAq\n/5138stOOCHMCfWnP21zt2nTwu38776Df/0rtJyalUK8IlIoakEtBf37w183PgTAKDrw2pcHJDki\nEZE0N3lyuHXfpk1+cnrMMTBuHLz/foHJaSwGt98O3bqF3Y49Fr7+OsxtquRUJLWoBTXB1q2Dpf/9\nlKOZCMDzdXtxj/qeiogUTywWWkxHj84va9sW+vYNLafbyDSnTg3dUpctCyPyhw2DP/+5lGIWkSJT\nC2oCxWKhC9S16x8GYHHDttzz7lHqeyoiUlTTp4eMsmXL/OT0iCPg7bfz+59uIzn9/ns477yQnAJs\n3Ji/gJSIpCYlqAkSi0G7dvDjxNmcuGoEAA0f76XkVESkKL76Cjp3hubNYUS4ltKqVRj49Mkn4UK7\njcTUHQYNClOd7rNP/vKkWhxFJPXpFn+CjBwJP/8MfXmUcji/1G3Kbh06JDssEZH0MGtWGOj0yish\n04SQpPbtG7LLHXQaXbAALr885Lcvvhh2icVC42v79locRSTVKUFNkCpVYE/7kYt8MABrr7wZyqnB\nWkRku+bMCUno0KFhHiiAQw4JyWrHjju8jublwT//GQZDde0aFouqXTvUtWihxFQkXShBTYB16+Dp\np+GVo5+kyoQNbKq7B43u6JHssEREUte8eXDvvTBkSFhvFOCgg0Ji2qlTob7gz54NPXuGqaPeeCN0\nSxWR9KQmvQS48UaouHENbb/4JwAVb7hWMz+LiBRkwQK47LIwW/7gwSE5bdo0tKB++SWce+4Ok9Oc\nHHjoodA62qpV2E3JqUh6UwtqCXv/fRgwAK5hEBVYSm7V6pS/8spkhyUikloWL4Z+/cJSTps2hbJ9\n9w1TSHXrBhUK9+dp+nS45BJYvRrefTfMOCUi6U8tqCXsnnugPDncyKMATG5+KdSpk9ygRERSxXff\nwTXXwH77hW/zmzZB48ZhuP2sWdCjR6GS040bw/X2yCPDdH7Tpik5FSlL1IJagpYtg88/h4trvMY+\nqxaQQ3l2ufOGZIclIpJ8P/wADz4YktING0JZo0Zw551w0UVQsWKh32ry5NBqWr58mAK1ZcvEhCwi\nyaMW1BL05JNwRGun/15hYv6Vp3bmkDP2TnJUIiKJl5WVRXZ29u8rfvoJbr4ZmjQJ6z5v2AB77RWG\n2s+eHfqfFjI5XbcOevUKS5R27RoSVSWnIqUvOzubrKyshB7DfPP8coKZ1QRWrFixgpo1axZp348+\nCsvovXLFWDr0j3rnT52qOU1E0sDKlSupVasWQC13X5nseNLJNq+bv/wCjzwCTzwBa9eGsvr1ISsL\nLr20yANHJ0wII/Tr1g29AQ4+uOTOQUSKLtHXTd3iLwGxGJx+OqxfDzUGPBQKTzpJyamIZJ5Nm8LU\nUP37h5FLAPXqhYlJr7gCqlYt0tutWgW33QYvvBBmobruunBrX0TKNiWoJeCNN2DlSjiU6Ry/IbrF\n1atXcoMSEUmGChVg3LiQnO62G9x6K1x5JVSvXuS3ys4Oq0Htuy988UX4KSKZQQlqCahVK3yjvzn3\nEQDWNm1OtZNPTnJUIiJJYAYPPAAffxxG6++yS5HfYunSMJ/0iBHw8MOhR4AW4hPJLEpQS8Cbb8Jt\n3RfT/d/DIA+q3d1rh+tEi4iUWcccEx7FMGIEXHVVGPw0YwY0bFjCsYlIWlCCupO++io0FIy5+DHK\n5+WEq2nnzskOS0Qkrfz4I1x7bVjspH9/6N5d3/NFMplumuyEWCxM39fxuOVUHfJMKLzxxiLN5yci\nksnc4aWX8kflf/01nH++klORTKcW1GLaPHL/hx+gQ/UBsGY11K4dOkuJiMgOffst/OUvMGUKPPss\n/PnPyY5IRFKFWlCLadSokJxWYgOXrukfCq+8slgDAkREMol7SEibNQsD/b/+WsmpiPyWEtRiat8+\nzKbSnZeozw/kVawUJugTEZFtmjcvTBN9773w8ssweDDUqZPsqEQk1ShBLaaNG6FyxTwe3C1MLVXu\nwgtgzz2THJWISGrKzYXHHoPDDoMDDggj9E87LdlRiUiqUh/UYnrmGXj4+Dep99+ZoeCmm5IbkIhI\nipo5MyxT+vPPYVq+Y49NdkQikupSrgXVzBaY2Uwzi5nZVDM7Nyrf3czeMrPZZjbdzI6O26eqmQ01\nszlmNsvMOsXVmZk9YWZzo32v3tkYP/wQhgyBTvMeDgVnnQUHHrizbysiUqZs2gT9+kGrVtC2bVgN\nSsmpiBRGKrag5gGd3f3LrcofBD5x99PM7I/ACDNr7O65wM3Aenff38waA5PMbKy7LwN6AAe6+35m\ntisQi+pmFie4WCz0P22V8yn1Zn0YCrWsqYjIFrEYDBwYVjytWBHGj4fWrZMdlYikk5RrQQUsemyt\nMzAAwN2nAN8Bm7+Ld4mrWwCMB86O2+/ZqG4Z8ArQrbjBjRwJK1ZAL0Lr6eKGbeCoo4r7diIiZUos\nBu3ahQR1yRJ47jklpyJSdKmYoAIMMbMvzOxZM6trZnWACu7+U9w2C4FG0fNG0evNFhSyrsiaNIH9\nmcPZjABg0/W3FPetRETKnFGjQl9TgFWrIDs7ufGISHpKxQT1aHdvDrQEfgVeiMpTYl2RCRPgySaP\nUg5n/d5NaXJ9h2SHJCKSMjp0gAYNwvMGDUKXKBGRokq5Pqju/m30M9fMHgO+cfelZpZjZvXiWlEb\nA4ui5wuBvYEf4+o2f29fFNVNKmC/AmVlZVGpUiUA2rVrR7t27QBYuRLeG/oTz+Q+D0CVrJugXCrm\n+CKyI9nZ2WRHzXsbN25McjRlR4sWMGYMjB4dktMWLZIdkYikI3P3ZMewhZlVAyq6+4ro9Y1AB3c/\nzswGAQvdvY+ZHQ78B2gcJbK9gb3d/RIz2wf4BDg4SmwvBM4H2gG1ganAGe7+VQHHrwmsWLFiBTVr\n1vxdfAMHQrk+d3PZknuhXj1YuBCqVEnIv4WIlJ6VK1dSq1YtgFruvjLZ8aSTHV03RaRsSvR1M9Va\nUPcAXjezcoRb+vOAC6K62wh9U2cDG4Du0Qh+gIeBQWY2F8gBrnb3pVHdEOCPwBzCDAGPFJSc7kgs\nBn+7ew0zVj0VCq67TsmpiIiISAKkVILq7vMJfU8LqvuJ0ApaUN1aoOs26vKAa6NHsWweldrl50FU\nYym5VatT/sori/t2IiIiIrId6kBZCKNGwdKfc7iRRwGYclhPLR4tIhInKytrS59eESnbsrOzycrK\nSugxUqoParJtqy9VLAYP//EVhuZ1JYfyfDNmLs3OaJy0OEWkZKkPavGpD6pIZkr0dVMtqIVQYxfn\npryHAFh5amclpyIiIiIJpAS1EKb+fRytmApAnX5a1lREREQkkZSgFkKj4WFZU046SZP6iYiIiCSY\nEtQd+HXcdI5c9nZ40UutpyIiIiKJpgR1BxZf/wgAa5s2h5NPTnI0IiIiImWfEtTtmPHWYppNHwbA\nLT/1IjbNkhyRiIiISNmnBHU7VvTtT0VyWERDBi7vzOjRyY5IREREpOxTgroty5fTetpAAP7BDdRr\nUJH27ZMck4iIiEgGUIK6LQMHUnH9alaWr80eWZcyZowG8IuIiIiUhgrJDiAlbdgA/fsDMONPV3Lb\n/TWSHJCIiIhI5lALakGGD4clS9hAJcpff22yoxERERHJKEpQC/L44wC8VO4C/nBa/SQHIyIiIpJZ\nlKAWJDcXgHEtb6Jy5STHIiIiIpJh1Ae1IFOm0KfbLPZudWCyIxERERHJOGpBLUi5cgyedSRHH53s\nQEREREQyjxLUArz7LixYANWrJzsSERERkcyjBLUAl14afnbpArFYcmMRERERyTRKUAuwfHn4+f33\naHlTERERkVKmBLUAe+4ZfjZogJY3FRERESllGsVfgOHDYdy4kJxqeVMRERGR0qUEtQDNm6MR/CIi\nIiJJolv8IiIiIpJSlKCKiIiISEpRgioiIiIiKUUJqoiIiIikFCWoIiIiIpJSlKCKiIiISEpRgioi\nIr9jZoeb2cTo0TfZ8YhIZtE8qCIiUpCp7v4nADN7z8x2cffVyQ5KRDKDWlBFROR33D0XwMzKA98B\na5MbkYhkkoxJUM1sPzP7yMy+MbNJZnZQsmMSEUllZtYN+BpY7u55yY5HRDJHxiSowEBggLsfADwE\nvJDkeBIuOzs72SGUGJ1LaipL5yK/5+7DomvmXmbWrCj76ncj8+gzl5KUEQmqme0OtAJeAnD314GG\nZtYkqYElWFm6WOhcUlNZOpeyzMz6m9l8M8szs8O2qivw7pKZVYrbbBWwvijH1O9G5tFnLiUpIxJU\noCGwZKtbVIuARkmKR0SkNL0KHAUsKKBuW3eXOpjZODMbDyx29/+VRqCJlqgkamfft6j7F3b7HW23\nvfri1qUafeaFr0+lz1yj+AuwcuXKZIdQIjZu3KhzSUE6l9RTFs5he9x9IoCZWXx53N2lk6PtXjez\nJ82sibu/BrxW2GNs/W+Yqr8bo0ePpk2bNin3vkXdv7Db72i77dUXtU6feWL3T7XPPNGftbl7Qg+Q\nCqKL8BygzuZWVDNbAhzl7vPittsL+DY5UYpICvg/d/8u2UEkipnNB85y9+nR65bAS+5+UNw2k4Bb\n3X18Id9T102RzJaQ62ZGtKC6+89mNhXoAbxgZucQblnN22rT74H/I/S3EpHMUoNwDZCi0XVTJHMl\n7LqZEQlq5C/AYDPLAlYAF2+9gYfm5DLbeiIi25V69yYTbzFQ38zKxfXRb0Too18oum6KZLSEXTcz\nJkF199lA22THISKSKopwd0lEpFRlRB9UEZFMZmYDgDOAPYBfgVXu3jSqawoMBuoS3V1y96+SFKqI\nCJA500ztULJXmjKzymY2wsxmmVnMzLLNbN+obncze8vMZpvZdDM7Om6/qmY21MzmRPt2iqszM3vC\nzOZG+1691THvjOrmmNl9W9X1jPaZY2YDo+UOi3NeF0dzL3ZI13Mxs0rRsWeb2Rdm9mK6nkv0Hqeb\n2efR79l0M7sgXc7HtjGfZ6rFbmZnmtlMC9eT18xsl+2dV6K5+1/cvaG7V3L3+puT06hutru3dfcD\n3L11aSSnZna9mb2b6ONI8pnZ3mb2vZmNNbPByY5HSoeZXW5m70Wfe7HyB9xdj9CK/D7QI3reCfis\nlI9fGTg17vXVwLjo+SDg7uj5Hwn9xspHr+8CBkXPGwM/ArtGry8A3o2e70qYA/Gg6PUxwJdAFaAS\nMBk4Larbh9CnbPfo9UjgymKc097AR9GjQ1T2XLqdC/APoH/c63rpei7Rfr8CzeI+o3VA9XT4PQP+\nBDQA5gGHxZWnzGcR/Vv+AOwfvX4CeKg0ryep/CB0LXseeCfZsehRKp/33sAzyY5Dj1L9zBsCj+/s\n+6gFldRYacrdN7j723FFnxL+YwOcCwyItptC+MN4bFTXJa5uATAeODuq6ww8G9UtA14BusXVDXH3\n9e6+kZCcbK7rBIx095+j1wPi6grFzAz4F3ANsDGuqnM6nYuZVQMuAe7YXObuP6XjucTJIyRjALWA\nXwifUcr/nrn7RHf/HrCtqlLpszgNmOruc6LX/9zReWWYHsDQZAchpepUM/vAzM5LdiBSKk4BdjGz\n982sd3HfRAlqkIorTf0VeMPM6gAV4pIigIXkx9Yoer3ZggTXFdaNwIfuHttckKbnsi+wFLjDzCZH\nF9kT0vRcNusKjDCzBcAE4ELCVCFpeT4p+FkUVLenmWX89Tb64nqKu7/L779kSNn0PdAUaAdcYWa7\n7mB7SX/1gDx3PxHYx7ZaXrmwMv6CmYosTIW1L5CV7FiKw8yaEVqY7k92LCWgAqEle4a7H0744vBy\nVJ52f2CjvkB3Ah3dvTFwEvBv0vR8JO38GRiV7CCk9Lj7puguxHrgQ8LfNinbVgAfRM8nAAcU502U\noAZb5gKMKyvSXIAlxcxuBjoS+qOud/elQI6Z1YvbrHFcbAvJ7wqwdd2iBNQVxtHR/nMsrFxzJPAM\n4ZZpup3LIiCX6Jaku08jtIgdCmxKs3MB+ANQ390/gi23w78FDiM9z4cU/D+yKHq92T78/g5NFOMB\n+wAABiJJREFU2tvWgLWobluDTg8ALjKzt4DmZnZJqQcuxVacz9zMqkc/jfz+4ZImivn//BPC30ii\nnwuKdfBkd6ZNlQcwFrgwen4OpTxIKjrujcAUoNZW5YOA3tHzw/ntAJDe5A8A2YcwOKNO9PpC4F3C\nF5E60S/J5sExxxIGgFQlDNCaDJwe9z7fEprpjTAA5KqdOK9xQPt0PRfgbX47OOYnoH6anks9wrfb\nA6PX+xH6oP5fOp0PMJ/fDpJKmdiBXaJjNI1el8lBUmxjwFpUt8NBp2iQVNo9ivOZAydH/3c+Am5I\n9jnokfjPPHr9D8Lf/gHFPnayTz5VHoQ+Mh8D3wCfbf4jVYrH34sweGUOMBWIAZ9EdfWAbGB29Afz\nmLj9qhFuOc8FZgGd4urKRX8c/xe97zVbHfPOqG4ucN9WdT2j8jmE1s/yO3FuY8kfxZ9250JIRsYC\n06PPpWO6nkv0Hl3izuULoEu6nA9hMNJiwqCuJcDsVIwdOBOYGcXzH6BGaV5PSvPB778s7A4sB8rF\nlS0BmiQ7Vj30meuRPp+5JuoXEZFii7rxnOXu06PXLYGX3P2guG0mAbe6+/jkRCklSZ955knGZ64+\nqCIiIiKSUpSgiohISUqZQadSavSZZ56Ef+ZKUEVEpMR4WLxgKmFCfszsHGCxu89LamCSMPrMM09p\nfObqgyoiIkVmZgOAM4A9CMvnrnL3plFdU2AwUJcwa8TF7v5VkkKVEqLPPPMk8zNXgioiIiIiKUW3\n+EVEREQkpShBlZRgZjEzm2pmX5lZTvQ8ZmbDdrDf8WZ2ciGPcaKZTS6gvImZrYzv7B2tnPFM3Ouj\nzGyn+9aYWXMzm7Oz7yMiIlKWVUh2ACIA7t4CwMz2BmLu3rKQu54AVCGsBlSoQxVw7HlmtpSwDN9n\nZtYYWAq0jdvseMKqGYVmZuXdPbcwMYiIiEg+taBKWjCz28zsSzObbmYvmNkuZtYKuBToEbW43m5m\nlczsbTP7LNr2RTOrUohDjAeOi54fB/wXWG5mDeLKxkWx7GFmI8zsi+ixZT1xM1tsZg9EExY/F5Xd\na2ZzzOwz4Ny4beuZ2TvRe0yLb7EVERHJZEpQJeWZ2ZnAecCR7n4YkAPc7+6fA/8CXnT3lu7+gLtv\nBLq6e+to23XAVYU4zDhCKynRz3HAB8DxZlYJaEN+C+pTwHR3bw6cAvSNVtXYrLa7H+HuF5lZB6A9\n0NzdWwP7xm13AfCNuzd39z8AtxTtX0ZERKRsUoIq6eAk4GV3XxO9fhoosN9p1I/0lqj/6nSgHfCH\nQhxjHHCUmVUg3Nr/GJhASFaPABa5+49x8QwEiMreiMo2ez7u+YnAK+6+Nno9MK7uE6C9mT0UJeFr\nEBERESWokpa214ezB3AUcFTUgtqf0Ed1+2/ovgj4CegKLHH39cBHwNGEJHXsdo6/9evVhYnd3T8C\nWgCTgc7AZzuKU0REJBMoQZVUZFu9fg/oYmbVo9dXANnR85VArbhtdwV+cfe1ZlaTcBu9sMYDd0U/\ncffVwHKgO1H/07h4LoPQjxToyLYHab0HdDaz6lHr7mWbK8xsH2C1u78KXAccaGbVihCviIhImaQE\nVVLRb1ok3X0M8BLwaXTbvgohkQR4HWi9eZAU4fZ6LTP7GhhN6EdaWOOA/YgS1MgHhH6j8QnqNUAL\nM/uCkJje7e6xbcQ+GhgFxIBPgflx1ScAn5tZLDrO9XFdAURERDKWVpISERFJQ9GXWwcqAwcA0wl3\noGa5e7ft7Hc8UMHddzg9n5mdCDzo7odvVd4EmEYYFJoXlc0H3nX3y6PXRwFD3L1Jcc4v7ljNgdfc\nff+deR9JL5oHVUREJA1p/mgpy3SLX0REpAzS/NGSzpSgioiIlDGaP1rSnRJUERGRskfzR0taU4Iq\nIiJS9mn+aEkrSlBFRETSn+aPljJFCaqIiEj60/zRUqZoHlQRERERSSlqQRURERGRlKIEVURERERS\nihJUEREREUkpSlBFREREJKUoQRURERGRlKIEVURERERSihJUEREREUkpSlBFREREJKUoQRURERGR\nlKIEVURERERSihJUEREREUkpSlBFREREJKUoQRURERGRlPL/S5AyPEv8YgoAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "#same data on linear-linear, then log-log, can read off exponent from slope on right\n", "plt.figure(figsize=(8.5,4))\n", "..." ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.13" } }, "nbformat": 4, "nbformat_minor": 0 }