L1.5.ipynb

{
 "metadata": {
  "name": "",
  "signature": "sha256:9716ad7e9679a9b787f0cdcb8f673365e681c2ea6e22759a6d7681b353da6702"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "heading",
     "level": 1,
     "metadata": {},
     "source": [
      "J\u0119zyk Python - Wyk\u0142ad 1.5"
     ]
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "TYPOWANIE"
     ]
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Wynik operacji \n",
      "\n",
      "    '3' + 3\n",
      "    \n",
      "w r\u00f3\u017cnych j\u0119zykach"
     ]
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Javascript: `'3' + 3 == '33'`"
     ]
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "PHP: `'3' + 3 == 6`"
     ]
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "LUA: `'3' + 3 == 6`"
     ]
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Ruby: `'3' + 3` **[error]**"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "'3' + 3"
     ],
     "language": "python",
     "metadata": {},
     "outputs": []
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "AUTOMATYCZNE RZUTOWANIE"
     ]
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Any object can be tested for truth value, for use in an if or while condition or as operand of the Boolean operations below. The following values are considered false:\n",
      "* None\n",
      "* False\n",
      "* zero of any numeric type, for example, 0, 0L, 0.0, 0j.\n",
      "* any empty sequence, for example, '', (), [].\n",
      "* any empty mapping, for example, {}.\n",
      "* instances of user-defined classes, if the class defines a \\_\\_nonzero\\_\\_() or \\_\\_len\\_\\_() method, when that method returns the integer zero or bool value False. \n",
      "\n",
      "All other values are considered true \u2014 so objects of many types are always true"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "doc_tests = [None, False, 0, 0L, 0., 0j, '', \"\", (), [], {}, list(), dict()]\n",
      "\n",
      "for test in doc_tests:\n",
      "    if test:\n",
      "        print(test, 'evaluates to True')\n",
      "    else:\n",
      "        print(test, 'evaluates to False')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": []
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "class Foo1(object):\n",
      "    \n",
      "    pass\n",
      "\n",
      "class Foo2(object):\n",
      "    \n",
      "    def __nonzero__(self):\n",
      "        return False\n",
      "\n",
      "class Foo3(object):\n",
      "    \n",
      "    def __len__(self):\n",
      "        return 0\n",
      "    \n",
      "class Foo4(object):\n",
      "    \n",
      "    def __nonzero__(self):\n",
      "        return False\n",
      "    \n",
      "    def __len__(self):\n",
      "        return 1\n",
      "    \n",
      "def bar():\n",
      "    pass\n",
      "\n",
      "others = [set([1]), 1, 123, -1, \"ala ma kota\", \"ala\"[0], \"ala\"[:0], Foo1(), Foo2(), Foo3(), Foo4(), str, bar, Foo1]\n",
      "\n",
      "for test in others:\n",
      "    if test:\n",
      "        print(test, 'evaluates to True')\n",
      "    else:\n",
      "        print(test, 'evaluates to False')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": []
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "WYSZUKIWANIE STRINGA"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import string\n",
      "import time\n",
      "from math import sqrt\n",
      "import random\n",
      "\n",
      "def random_word(n):\n",
      "    return ''.join(string.lowercase[random.randint(0,25)] for _ in range(n))\n",
      "    \n",
      "def random_bits(n):\n",
      "    return ''.join(str(random.randint(0,1)) for _ in range(n))\n",
      "\n",
      "def test_string_search(n, m=1000, warmstart=100, gen=random_word):\n",
      "    \"\"\"Zwraca czas wykonywania m wyszukiwan losowego slowa dlugosci sqrt(n) w losowym slowie dlugosci n\"\"\"\n",
      "    dummy = 0\n",
      "    for i in range(m+warmstart):\n",
      "\n",
      "        if i == warmstart:\n",
      "            start = time.time() #warm start\n",
      "\n",
      "        haystack = gen(n)\n",
      "        needle = gen(int(sqrt(n)))\n",
      "        \n",
      "        if needle in haystack: #lets hope python does not optimize too much\n",
      "            dummy += 1\n",
      "        \n",
      "    return time.time() - start, dummy\n",
      "\n",
      "for gen in [random_word, random_bits]:\n",
      "\n",
      "    lengths = []\n",
      "    times = []\n",
      "    for length in range(20, 500, 10):\n",
      "        lengths.append(length)\n",
      "        times.append(test_string_search(length, gen=gen)[0])\n",
      "\n",
      "    %matplotlib inline\n",
      "\n",
      "    from matplotlib import pyplot as plt\n",
      "\n",
      "    plt.figure()\n",
      "    plt.title(gen.__name__)\n",
      "    plt.plot(lengths, times, linewidth=5)\n",
      "    plt.show()"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAAEKCAYAAAAYd05sAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAH+1JREFUeJzt3Xmc1WP/x/HXNO0pd1maqLRI8ZOkZMvdKEvTTbit4Q5R\nWuyELDX4IdxkKcIvla0s3ZJUsg1F5W7BKKVStClKiyamab6/Pz7fMaczZ5uZs37P+/l4nMec8/1e\n58zlq8f7XHNd1/e6QEREREREREREREREREREREREREREREREROIlF3g50ZWIkWKgRaIrIamrSqIr\nIFJJTqIrIJKsFPASb1UTXYEkpGsiMaGAl3hYDdwGfAP8DtwFrAC2A4uBc3zKXgHMBh4FtgA/AN19\nzjcHPnXfOxPY3+939XQ/8zfgE6CNXz1udeuxAxgDNASmA9uAD4C/hflvGQ/c7D4/GOtGGei+bgls\n9inbF1juHnsHaORzruR9y4Fl7rHBwHpgLdAnTD1ERJLCamAhFog1gfOBLPfchVjoN3RfXwEUAlcB\nGUB/YJ3PZ80B/g1UA07Ggv4l99xh7md1AzKxwFxOaQt5FfAFcABwELDRrVc7oAbwETA0zH/LlcAU\n9/kl2BfVRPd1H+Bt93lX4BfgaKA68BT2xVSiGHgf+0KpgX2J/QwcAdQGXkN98CKSAlZhwR3MIqzl\njVtuuc+52ljQHQg0BXYDtXzOv0ppwN9DadiCfUGsBf7uU49ePuffAkb5vL6W0oAOpiX2l0UG8CzQ\nD1jjnhsP3Og+HwMM93lfHeyLq6n7uhjI9jn/IvCgz+tWKOClktRFI/Gyxud5byzUf3MfRwL7+Zz/\n2ed5gftzH6zV/Ruwy+f8jz7PDwJ+8nntuL/3YJ9jG32e7/J7/Yf7e0JZCezEWuYnA1OxbpXDsC+S\nklZ6I7+67cS6anzr4ntNGvm99v3vEKkQBbzES8lsl0OA54FBQAOgPvAt1iIOZ4NbvrbPsUN8nq/z\ne50BNGHvLh5/kfxef58CF2DdROvd11e4dfvKLbMeaObznjrYl5hvXXxnAG2gtHWP33ORClHAS7zV\nwYLtV+zf35VYCz4SPwLzgXuxcO0MnOlz/k3gH1j/dzXgFqxV/kU0Ku7jU6w75zP3dZ77ehaloT0B\n+28r6d9/EJhL8Jb5G9iXxOHYF9iwKNdZ0lC4gG+CzURYjLWyrg9QJhubgbDIfdwdxfqJ9ywBHsMG\nS3/Gwn22z3mHsnPbfV9fAhyH9YMPxfq9SywDLgOexgY4/wGcBRSFqI/j9zySefWfYV05JQH/OTYu\n8JlPmY+wMYFJWGu+OXBxkN8LMAN4AvgY+N59v+b4S0xlYX2NYP+gl2EtDF/ZlM4qEBGRJBGuBf8z\npX2KvwPfYQNZ/irSjykiIkmiGdYH6j/LoAs2O+BrYBo2j1cklV2K3Qjl/8hPZKVEYmUfbHDrnADn\n6lI6qyEH6z8UEZEEi6RrpRo213c6NggUziqgAzYI9peWLVs6K1euLHcFRUTS3Erg0Iq8MVwffAZ2\nR94Sgod7Q0q/KDq5z7f4F1q5ciWO4+jhOAwbNizhdUiWh66FroWuRegHdvd0hYRbxe4kbNrZN9gU\nSIA7Kb0J4zlsXZEB2FS0AvaeCiYiIgkSLuBnE76VP4q91/MQEZEkoDtZEyA7OzvRVUgauhaldC1K\n6VpERzznrztuf5KIiEQoIyMDKpjVasGLiHiUAl5ExKMU8CIiHqWAFxHxKAW8iIhHKeBFRDxKAS8i\n4lEKeBERj1LAi4h4lAJeRMSjFPAiIh6lgBcR8SgFvIiIRyngRUQ8SgEvIuJRCngREY9SwIuIeJQC\nXkTEoxTwIiIepYAXEfEoBbyIiEcp4EVEPEoBLyLiUQp4ERGPUsCLiHiUAl5ExKMU8CIiHqWAFxHx\nKAW8iIhHKeBFRDxKAS8i4lEKeBERj1LAi4h4lAJeRCSBVq+GP/6IzWeHC/gmwCfAYuBb4Pog5Z4C\nlgNfA+2jVjsREQ9bvRpOPhl69oSdO6P/+eECfjdwE/A/wPHAIOBwvzI9gEOBVkA/4Nko11FExHPW\nr4du3WDtWvjgAzj9dNi6Nbq/I1zA/wx85T7/HfgOOMivTE9gvPt8HvA3oGG0Kigi4jW//gqnnQY/\n/FB67IsvoGtX+OWX6P2e8vTBN8O6X+b5HT8YWOPzei3QuHLVEhFJLYWF8PXX4fvTt22D7t1hyZKy\n5zZsgO3bo1enqhGW2wd4C7gBa8n7y/B77QT6kNzc3L+eZ2dnk52dHeGvFxFJXv/9L+TkwObNUKcO\n3HYbDB4MtWrtXa6gAM48ExYsKPsZDRrAhx/CmjV5vPxyXlTq5R/MgVQDpgLTgScCnB8N5AET3ddL\ngS7ARr9yjuMEzH0RkZS1fTscfrj1qftq3hyeeALOOgsyMuDPP20wdebMsp9Rty58/DF07Fj2XEZG\nBkSW1WWE66LJAMYASwgc7gBTgN7u8+OBrZQNdxERT7rzzrLhDrBqFZx9NvToYd0xvXoFDvdateC9\n9wKHe2WF+1boDHwGfENpt8udQFP3+XPuz5FAd2AncCWwMMBnqQUvIp4yZw6cdBJUNNqqV4cpU+CM\nM4KXqUwLvkJvqiAFvIh4RmEhHHMMLF5csfdnZsKbb8K554YuF8suGhERCeDRRwOHe6tWkb1/7Njw\n4V5ZCngRkXL6/nu4//6yxzt2tP72SZPgkEOCv/+ZZ+Bf/4pd/Uoo4EVEysFx4JprbFaMr8xMeP55\nqFoV/vlPC/qhQ6FGjb3LPfwwDBgQn7qqD15EpBzGjoU+fcoeHzwYHnmk7PFVq+DZZ22mzUUX2bTJ\n8tAgq4hIHGzaBG3awG+/7X28eXPIz7ebnKJNg6wiInFw441lwx2shR6LcK8sBbyISBh//mnhPmFC\n2XOXXhp6HnsiqYtGRCSEpUvh4ottITF/DRrY+QMOiN3vVxeNiEiUOQ6MGQMdOgQOd4DHH49tuFdW\npKtJioikja1boX9/eP314GV697ZHMlMLXkTEx0cfQfv2wcO9WjUYMQLGjbNVIpOZWvAi4llbt8In\nn0C9etC5c9mbjnzl58Ptt8P06cHLHHYYTJxoXwCpQC14EfGkefPgiCPsrtJTT4UDD4TLL4dp02yh\nsBJr19qNS+3ahQ73Pn1so45UCXfQLBoR8aBVq+C444Lvb1q/vgV//fowcmTobfbq1YPnnrOZNImg\nO1lFRFxbt8KJJ8J331X+s048EV55xe5UTRRNkxQRAXbvhgsvrHy4N2sGr70Gs2YlNtwrS4OsIuIJ\njgPXXQcffFDxz6hfH+65BwYODD0gmyrUghcRTxgxwvrK/dWtC19+Cf/5j63mWLt22TI1asBtt8HK\nlXDTTd4Id1AfvIh4wDvv2O5I/hFTpQpMnQo5OaXHdu60mTRvvQXLl8MJJ9j0yKZNSUoaZBWRtLVw\nIZx8MhQUlD03ciQMGhT/OkWTBllFJO0UF8Po0dC1a+Bwv+661A/3ylILXkRSTn6+bZs3Z07g8z16\nWLdNVQ9MI1ELXkTSQkEBDBkCxxwTPNzbtrXlBLwQ7pWlSyAiKeH9922z6lWrgpc55BAbVK1bN371\nSmZqwYtI0nIcWyysWzfo3j10uF95pa0Vk6yzYRJBLXgRSTqOAzNmwP/+L3zxReiyrVvbYGt2dlyq\nllIU8CKSNIqLYcoUC/YFC0KXrV4d7rwT7rjDOzcmRZsCXkSSwqJFcNVV9jOc7GxrtbduHfNqpTT1\nwYtIQu3eDffeC506hQ/3Nm3g1Vfh448V7pFQC15EEiY/3zbhCBfsRx0Fd99ta7hnZsanbl6gFryI\nxF1RETz0EHToEDrcO3a0G5YWLYILLlC4l5da8CISV/n5cPXVtsJjMG3bwqOPwumnJ//G1slMLXgR\niYvFi2253nbtgod7ZibcdRfMnw9nnKFwryy14EUkphYvhvvugzffLLucr6/DD4fx4+HYY+NXN69T\nC15EYmLxYtuoum1beOON4OFepYpttrFwocI92iIJ+BeBjUB+kPPZwDZgkfu4Oyo1E5GUtHo19O5t\nwf7666Fb7a1awezZ8PDDULNm3KqYNiLpohkLPA28FKLMp0DPqNRIRFLS5s3wwAMwahQUFoYum5Vl\nd6Bec42CPZYiCfhZQLMwZTQUIuJRRUW2Hszu3XDQQfaoV690ALSgAJ58EoYPh+3bQ39WSbD36we1\nasW+7ukuGoOsDnAi8DWwDrgVWBKFzxWRBJs1y5YPWL587+O1a5eG/YoVsH596M9RsCdGpC3vZsC7\nQNsA5+oCe4ACIAd4EjgsQDnt6CSSIgoKbLrik0+G7kMPZ//9bUGw/v0V7BVVmR2dotGC3+HzfDrw\nDNAA2OJfMDc396/n2dnZZGt9T5Gk88UXcMUVZVvt5VG7Ntx8MwwebN05Erm8vDzy8vKi8lnRaME3\nBDZhXTWdgDcI3GevFrxIEtu1C+65Bx5/vOKt9sxM6NsXhg6FRo2iW790FesW/ASgC7A/sAYYBlRz\nzz0HnA8MAIqwbpqLK1IREUmcL7+0qY3LlgUv07KlBf/69fDHH2XPn3eezaLRKo/JI56zX9SCF4ki\nx7FHlUrcrlhcDI88Yi33oqLAZWrUsOC+8UZroTsObN1qQb9+vbX8jz5aW+XFSmVa8Ap4kRQ0frwN\nXhYU2BK6ubnQpEn5PmP9evjXv2xt9WCOOw7GjbN12CUxFPAiaeT9920Dal81asD118OQIVC/fvjP\nePdd26R68+bA52vUsPVjbrlFS/QmmgJeJE389pstAbBuXeDz9etby/7aawPfIbprl81sGTUq+O/o\n2NH+QjjiiOjUWSpHAS+SJnr3hpdfDl+uaVPbKWnHDti0CTZutJ9r19qXRCAZGbbo13332YbWkhwU\n8CJp4O23rb89Fho1si+Obt1i8/lScQp4EY/btAmOPBJ++WXv43XqwJ49gactRuof/4CxY+GAAypX\nR4mNygS81oMXSXKOAwMGlA13sL705cuhT5/yT5esUQOeesoGXBXu3qQWvEiSe/VVuOyyssd79oTJ\nk0tXdfz2W5tFM3Vq6M+rUsX2Oh0+3LbPk+SmLhoRj1q3zrpmtm7d+/h++1mgZ2WVfc/cufDhhzaY\n2rChPQ48sPRxwAFagz2VKOBFPMhxoEcPmDGj7Lk33oALLoh/nST+1Acv4jGOAw8+GDjcL75Y4S6R\nUQteJMmsW2eDpjNnlj2XlWVdM/vtF/96SWKoBS/iERMn2p2qgcId4IUXFO4SOQW8SBLYsgV69bJH\nsDtN+/aFM8+Mb70ktUVjRycRqYSZM23hr1D7mg4YYNvniZSHWvAiCbJ6NZx/PpxxRvBwz8qCadPg\nmWegWrXAZUSCUcCLxNnOnbbBRps2MGlS8HLnnw/5+ZCTE7+6ibeoi0YkThwHJkywFRuDLfcLsO++\nMHIkXHpp6V2qIhWhgBeJg2++gYED4fPPQ5fr2tV2UCrv7kwigaiLRiTGxo6FTp1Ch3vDhlbugw8U\n7hI9asGLxEhhIdxwA4weHbxM9epw0022C1O9evGrm6QHBbxIDKxfb4Okc+YEL9OzJzz2GBx6aPzq\nJelFAS8SZbNn21oxP/8c+HybNvDEEzY9UiSW1AcvEiWOY7NfTjkleLj37g0LFyrcJT7UgheJgq1b\n7W7TiRMDn69a1VrtAwdq6qPEjwJepJJmzbIdl376KfD5rCx4803o3Dm+9RJRF41IBe3eDXffDdnZ\nwcP9hBNgwQKFuySGWvAiFbBihd1p+uWXwcsMGGDdMtWrx69eIr7Ughcph+JiGDMG2rcPHu716tlG\n2c88o3CXxFILXiQCjgNvvw3DhtmOSsGcdBK88go0axa3qokEpRa8SAiOA1OnQocOcN55wcM9MxPu\nuw/y8hTukjzUghcJwHFsXZihQ2HevNBlW7SwLpnjj49P3UQipYAX8bNkCVx3HXz8cfiyvXvD009r\nHRlJTuqiEXHt2AG33grt2oUP95NOsjLjxyvcJXmpBS9pz3HsDtRbboENG0KXPfZYuP9+OP103ZEq\nyU8BL2lt8WIYNAg+/TR0uaOPtkHUM89UsEvqiKSL5kVgI5AfosxTwHLga6B9FOolElObNlk/e7t2\nocO9RQt46y27G/WssxTukloiCfixQPcQ53sAhwKtgH7As1Gol0hMFBTAAw/YGuwjR8KePYHL1axp\nLfbFi216ZBWNVkkKiqSLZhbQLMT5nsB49/k84G9AQ6zVL5IU9uyxAdF77rHNOEI5+2wYMQKaN49P\n3URiJRp98AcDa3xerwUao4CXJOA4MG0a3HFH6DtQAVq2tCmPOTnxqZtIrEVrkNW/Z9KJ0ueKVMiu\nXbZkwBNP2Lz2UPbZx74AbrnFumZEvCIaAb8O8N0HvrF7rIzc3Ny/nmdnZ5OdnR2FXy9SasMGGDXK\nNrrevDl02cxMuOYaW1/mwAPjUz+RcPLy8sjLy4vKZ0U6J6AZ8C7QNsC5HsC17s/jgSfcn/4cx1HD\nXqLv999t043XXoPXX7d12sP55z/hwQehdevY10+kMjJs6laF5m9F0oKfAHQB9sf62ocB1dxzzwHT\nsHBfAewErqxIRUQiVVgIc+fCRx/ZY948KCqK7L0nnACPPmp3oop4XTxn9aoFLxVWWAiTJ8O4cTZv\nvaCgfO/v0gVuvllz2SX1xLoFL5IwP/4Izz9vm2xsLOe8rGrVoFcvuPFG26BDJN0o4CXp7NkDM2bA\ns8/aFMfy/uG3//62Xd6AAdCoUWzqKJIKFPCSNPbssXXV770XfvihfO9t0gS6dYPTToNzz4VatWJT\nR5FUooCXhHMcmDIF7rrLlgaIxH77wSmnWKh362ZLD6hvXWRvCnhJqE8/tZuM5s4NX7ZqVZve2L+/\nDZpqfRiR0BTwkhALFsDdd1tfezhNm9oNSX36QFZW7Osm4hUKeImb4mIL9H//Gz75JHz500+3JX1z\ncuyuUxEpHwW8xNyff9rg6WOPhV8XBuwmpIcegpNPjn3dRLxMAS8xs22brQvz1FORzWFv29aCvUcP\nDZiKRIMCXmJi8mQYODD8Hqdg667ff7/dlKSBU5HoUcBLVJVshffGG+HLtmljS/T27g3Vq8e+biLp\nRgEvUeE4tprjDTeEX6Y3OxtuvdUGT9ViF4kdBbxU2tq1Njf9vfeCl8nMhAsusBZ7x47xq5tIOlPA\nS4UVF9tCYLffDtu3By+Xk2MbXLdoEb+6iYgCXipo6VLo2xdmzw5epkEDePJJuPRSzYoRSQT1gEq5\nFBbajJd27UKH+4UXwnffwWWXKdxFEkUteInY3Llw9dWhFwTLyrJlfs85J371EpHA1IKXkIqLraV+\n1VVw4omhw71PH7tTVeEukhzUgpeA8vNt2uOECbarUiiHHmqDraecEp+6iUhkFPDyly1b4LnnLNi/\n/TZ8+cxMGDwYhg7VBhsiyUgBLwBMnw6XXw6//BJZ+Q4d4P/+D44+Orb1EpGKUx98misstJuPevSI\nLNzr1rXlfufOVbiLJDu14NPYihW2wNf8+aHLZWbaXqeXXGIDqHXrxqd+IlI5Cvg09dprtkvS778H\nL3PCCRbqF14IBx4Yv7qJSHQo4NPMzz/DkCEwblzwMmecYeu4t2wZt2qJSAwo4NPAnj0wcya88AK8\n+y4UFQUuV7Wqbbhx881a5VHECxTwKcpxYMQIC+3t26FVK2jd2tZYb93aHlWrWkv9xRfhp59Cf16L\nFjbnvVOnuFRfROIgnquEOI7jxPHXeduIEdbSjoZevWD0aKhXLzqfJyLRk2GLOVUoqxXwKej7722x\nrz/+qNzn1K5ty/hecYUWBBNJVpUJeHXRpJg9e2zNl8qEe6NG9hn9+0PjxtGrm4gkFwV8ihk5Ej7/\nvPzvq1LFbmbq29d+VtX/eRHPUxdNClmxAo46Cnbt2vt4ixZ2d+n339tGHMuW2c9t22yq42WXWYtd\nrXWR1KM++DRQXGybVc+aVfbcJ5/YOV+OY+/JzIxH7UQkVioT8JrtnCJGjQoc7oMGlQ13sEFThbtI\nelMLPgX88AO0bQsFBXsfb9bM1m3fZ5+EVEtE4kAteA8rLrbdlPzDHWDMGIW7iAQXScB3B5YCy4Hb\nA5zPBrYBi9zH3dGqXLras8c23Bg7Fi69FPLyypbp3x+6do171UQkhYRr9mcCy4BTgXXAf4FewHc+\nZbKBm4GeYT5LXTRBOA589hlMnmxL9y5cGLjFXqJpU/sC0LK9It4XyxudOgErgNXu64nA2ewd8BX+\n5QKbN9uyvZMmRf6eMWMU7iISXrgumoOBNT6v17rHfDnAicDXwDTgiKjVzuNmzLDB0/KEe79+cOqp\nsauTiHhHuBZ8JH0qC4EmQAGQA0wGDgtUMDc396/n2dnZZAea35cGCgrgttts6mN5nHOOLTImIt6V\nl5dHXqCBtwoI17VyPJCLDbQCDAGKgYdDvGcV0AHY4ndcffBYH/tll9ndpqHUrw8dO9rj2GPtoTtR\nRdJPLPvg5wOtgGbAeuAibJDVV0NgE9ba7+RWxD/c097OnfDoo/DAA8E33Khf3zbcOPVUW35AKzyK\nSGWEC/gi4FrgfWxGzRhsgPUa9/xzwPnAALdsAXBxTGqaooqKbMON3FzYsCF4udNOs2mRB/uPcIiI\nVJDuZI0Rx7Fpj0OGhO6OqVkTHnnElhzQNnki4k/rwSeZ2bNtEHXOnNDl2reHV16BIzTvSERiQAEf\nJZs3w8SJ8NJL8OWXoctmZtoXQG4uVK8el+qJSBpSwFdCYSFMnw7jx8PUqbB7d/j3XHCBDbS2ahX7\n+olIelPAl9NPP9myvZ99Bv/5D/z6a2Tv69LF+to7dYpt/URESijgQ3AcWLzY+tRnzbLHmjXh3+er\nbVsYPhxycjTtUUTiSwEfxLp1dufo/PkVe3/79nDjjbYapDbeEJFE0DTJALZvh86dbTON8sjKsrtU\ne/e2lruISGVpmmQUFRXBRRdFHu41a1pL//LL7Q7UqrqiIpIkFEc+HAduuMFWeQymalXo0MFa+J07\nwymnwL77xq+OIiKRUsD7ePJJeOaZssdr14bBg+Hvf4fjjoM6deJfNxGR8lIfvOudd+Dcc60V7ysj\nw5Yc6BluvyoRkRjQptuVtGABXHJJ2XAHePxxhbuIpKa0b8GvWWPdLoFWehw0CJ5+WvPXRSRxKtOC\nT8uALyiwgdRJk+Ddd2HHjrJlcnJgyhTNihGRxFLAR+D33+G99+Ctt2DaNAv5YI46yu5e1cbWIpJo\nCviQvxRGj7Z12bdtC1++USOYNw+aNIl93UREwtGNTkHs2AH9+tkyvpGoW9e6bBTuIuIFnp1Fk59v\nG1ZHEu716tnyAt9+azcxiYh4gSdb8GPHwsCB8Mcfwcs0aGBLDJx3HnTrBjVqxK9+IiLx4KmALyiw\nqY3jxgUvc8klcOWVtj57tWpxq5qISNx5IuCLi+Htt+Gee+C77wKXqVMHXngBevWKb91ERBIlpQO+\nqMj62B96CJYsCV7uyCPhzTehTZv41U1EJNFSMuD//NM2tx4+HH74IXTZyy+3BcRq145P3UREkkVK\nBbzjwIsvwrBhtuNSKDVrwqhR1t+upQZEJB2lTMD/+CNcfTV8+GH4sp06wfPPQ7t2sa+XiEiySvp5\n8I5jYd22bfhw79IFZs6EuXMV7iIiSd2Cj7TV3r073HWX7bAkIiImKVvwkbbazz0X5s+H6dMV7iIi\n/pKuBR9Jq71xY5vT3r17/OolIpJqkmo1yS1boGVL2Lo1eJk+fWyXJW10LSLpwDNb9jVoYK33QBo3\ntq6YMWMU7iIikUiqFjzArl3Qvj0sW1Z6TK12EUlXntvwY84cGzQ96CD1tYtIevNcwINtrXfaaWq1\ni0h682TAi4hI7AdZuwNLgeXA7UHKPOWe/xpoX5GKiIhIdIUL+ExgJBbyRwC9gMP9yvQADgVaAf2A\nZ6NcR8/Jy8tLdBWShq5FKV2LUroW0REu4DsBK4DVwG5gInC2X5mewHj3+Tzgb0DD6FXRe/SPt5Su\nRSldi1K6FtERLuAPBtb4vF7rHgtXpnHlqyYiIpURLuAjHRX1HwDQaKqISIKFG5k9HsjF+uABhgDF\nwMM+ZUYDeVj3DdiAbBdgo99nrQBaVryqIiJpaSU2zhl1Vd0PbwZUB74i8CDrNPf58cDcWFRERESi\nLwdYhrXAh7jHrnEfJUa6578Gjolr7UREREREJLoiuVHKS17Exh/yfY41AD4AvgdmYlNJSwzBrs1S\n4PQ41TFemgCfAIuBb4Hr3ePpeD1qYtOIvwKWAA+5x9PxWpTIBBYB77qv0/VarAa+wa7Fl+6xlLgW\nmVjXTTOgGoH78L3mZOxuXt+AfwS4zX1+OzDcfX4Edk2qYddoBUm2hHMlZQFHu8/3wbr6Did9r0dt\n92dVbKyqM+l7LQBuBl4Fpriv0/VarMIC3VdKXIsTgBk+r+9wH17XjL0DfimlN39lua/Bvol9/6qZ\ngQ1Ue9Vk4FR0PWoD/wX+h/S9Fo2BD4FTKG3Bp+u1WAXs53csKtci1skfyY1S6aAhpdNGN1L6P+4g\n7JqU8PL1aYb9ZTOP9L0eVbDW10ZKu67S9VqMAAZj065LpOu1cLAvu/lAX/dYVK5FrPdk1Q1PZTmE\nvi5evGb7AJOAG4AdfufS6XoUY11W+wLvY61XX+lyLc4ENmF9ztlByqTLtQA4CdgAHID1uy/1O1/h\naxHrFvw6bKCtRBP2/vZJFxuxP7MAGmH/uKHs9WnsHvOSali4v4x10UB6Xw+AbcB7QAfS81qciK1h\ntQqYAHTF/n2k47UAC3eAX4C3sTXAUuJaRHKjlBc1o+wga0m/2R2UHTCpDjTHrlU81+iPtQzgJezP\ncV/peD32p3QmRC3gM6Ab6XktfHWhtA8+Ha9FbaCu+7wO8Dk2MyZlrkWgG6W8bAKwHijExh+uxEbI\nPyTwlKc7sWuzFDgjrjWNvc5Yt8RX2J/ji7Bps+l4PdoCC7Fr8Q3W/wzpeS18daF0Fk06Xovm2L+J\nr7CpxCUZmY7XQkRERERERERERERERERERERERERERERERES86P8BL2rCtbsLGuQAAAAASUVORK5C\nYII=\n",
       "text": [
        "<matplotlib.figure.Figure at 0x7f71e8852f90>"
       ]
      },
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAAEKCAYAAAAYd05sAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAH25JREFUeJzt3X2cjPX+x/HXtEuF7qxwQil3iVLk5pSydaRScU5FN7px\nH926617ZTjcofjmioiKlQ8UhxamjtKlOORFO5D7KzZKSmw5i7fz++My2Y3eu2dndmWturvfz8ZhH\nM3N9d+frar3343t9r+8XREREREREREREREREREREREREREREREREJF6ygNdd/LzaQB5whMPxB4GX\nXOuNSJD0eHdAJMr88e5AIUODntcGvsP+3uXFpTfiKU5Vh4hbvFhk+OLdAfEGBbzEwwbgPuC/wK/A\nw8BaYDewHPhzUNuuwGfAM8AOrAK+LOj4qcAnga/9F1Cl0Gd1CHzPX4CPgdML9WNQoB97gFeAasA/\ngV3AXOD4CP9MPYDNwBZgYND7WRQMGc0P/Hdn4PNaAnUD/d8JbAemRvh5IiIJaQPwNVADOAq4Fqge\nONYZC/1qgdddgQNYgPqAPliQ5vsCGAGUAy7Agv61wLH6ge/1JyANuBdYQ8G/GtYD/wZOBE4CtgX6\n1QQ4EvgIeLSYP0ttbLjlDeBooDHwY+AzAYZQEPCnUHS8fgo2Tg9QHjivmM8TEUlo67HgdrIYq7wJ\ntFsTdKwCFpJVgZOBg1iw5nuDgoB/hMMrYh+wCbgwqB83BB2fBowNen0nMCPcH4SCgK8f9N5w4OXA\n8ywKAj6/bXDATwLGYb/sRKJKQzQSLxuDnt+ChfovgUdjICPo+Nag53sD/62EVd2/APuCjn8f9Pwk\n4Ieg1/7A5waH6bag5/sKvd4f+JxIBP95fgh8diTuw37x/AdYBnSL8OtEiqWAl3jJn+1yCjAeuAOo\nDJyABV0kFyJzAu0rBL13StDzzYVe+4BaHD7EU1hpL4CeXOh5qM8INcNnG9Ab+6VzG/A8cFop+yBy\nGAW8xFtFLPh+wn4eu2EVfCS+BxYCj2Fj8K2BK4OOvw1cAVwcOD4Qq8r/HY2OFzIYGypqhA0rvRmi\nzXZsiKZO0HudgJqB5zuxc6EplBIVxQV8LWzmwXKsqrrbod1obJx0KXBO1HonXvAtMBK7WLoVC/fP\ngo77KVr5Br++EZuNsgO7IDop6Ngq4CbgOSxcrwCuAnLD9Mdf6Hkk8+r92EyYtcCH2IyfD0N8j73A\nk8Dngf62BM4FvsRm1byD/R3bEMFnipRZdeDswPNK2F+YhoXatAfmBJ63xH5YRUQkycykYPpXvheB\n64Jer6RgipuIiMRJScbga2PDLwsKvV+Dw2cQbKJgTFEkFXTBhlAKP76JZ6dEoqUSdjHrzyGOvQuc\nH/T6Q6CpG50SERFnkawDUg6YDkzGhmgK24xdjM1XkxBTxOrUqeNft25dafooIuJl67AlLUqsuCEa\nH7Y+x7fAKIc2s7AbVQBaYVO9thVutG7dOvx+vx5+P0OGDIl7HxLloXOhc6FzEf7B4dNqS6S4Cv58\nbJrZf7E7DQEeouCmjnHYDJr22BSx/6E78UREEkJxAf8ZkV2IvTMKfRERkSjSnaxxkJmZGe8uJAyd\niwI6FwV0LqLDzY0H/IHxJBERCfLjj1C1auhjPp8PSpnVquBFROJk/34YOBDq1IE1a4pvX1Kq4EVE\n4mDxYrj5Zli+3F63bAmffQbpha6MqoIXEUkSubnw1FMW6PnhDrBgAQwd6vx1paEKXkTEJWvXwi23\nwBdfhD6eng5LlkCjRgXvqYIXEUlgfj+MGwdNmjiHe/nyVsGffnro46URyVIFIiJSSjk50LMnzJnj\n3KZJE3j9dTjzzOh+tip4EZEYeestaNzYOdx9Prj/fht/j3a4gyp4EZGo+/lnuPNOmDrVuc2pp8Jr\nr0Hr1rHrhyp4EZEomjPHqvZw4d6zJyxdGttwBwW8iEhU7NkDvXvDFVfA1q2h21StCrNmwUsvwTHH\nxL5PGqIRESmj2bNtSGbDBuc211wDL7wAJ57oWrdUwYuIlNZ330GHDnDllc7hfvzxMHkyvP22u+EO\nquBFREps3z4YPhyGDYPffnNu164dvPIK1IzTLtUKeBGRCPn98O670K8frF/v3K5CBRgxAvr0samQ\n8aKAFxGJwPr1cNddNt4eTmamXUStW6pdVKNLY/AiImEcOGBLCDRqFD7cTzoJpkyBefMSI9xBFbyI\niKP5822YZcUK5zbp6TBgADzyCFSq5F7fIqGAFxEp5Kef4L77YOLE8O3atoXnnovuAmHRpIAXEQnw\n+y3U770XduxwblejBowaZXPb43kRtTgKeBERbBimTx8blnGSlgb33ANZWe7ciVpWCngR8bR9+2yH\npeHD4eBB53YtWtia7mef7V7fykoBLyKe9eGH0Lev7bTk5Ljj7IamXr2sgk8mmiYpIp6TkwNdusAl\nl4QP9xtvhFWrbOgm2cIdVMGLiIccPAhjxsCQIbb6o5O6dW1hsLZt3etbLCjgRcQTPvnEVnxctsy5\nTblytsPSQw/B0Ue717dYUcCLSErbsgUGDbK7TMO58EJ48UVo2NCdfrlBY/AikpIOHoSRI6FBg/Dh\nXrkyTJgA2dmpFe6gCl5EUlB2NtxxB3z7rXMbn892YHryScjIcK1rrlLAi0jKyMmxu1DfeCN8uxYt\nYOxYOPdcd/oVLxqiEZGkl5sLf/ubrQkTLtyrVIGXX4Yvvkj9cAdV8CKSpH78ERYtssdbb8E33zi3\n9flsLvsTT9iYu1co4EUkYR06BNu2waZN9lixwgJ94ULYuDGy79GyJTz/PDRtGtu+JiIFvIgkjJkz\nrRrfsMECfcsWC/nSyMiw9WW6dYMjPDoYrYAXkYQwahT071/27+OF2TGRcnMlY7/f73fx40QkWSxY\nAK1b28XSsmje3GbHNG8enX4lAp8tOF+qrFbAi0hc7doF55xjm1qXRFoanHGGzYZp1gxatbJx9kTe\ngKM0yhLwGqIRkbjx+212S7hwP+EEqFmz4HHWWRboTZpAhQru9TUZKeBFJG5efRWmTi36/rnnwuTJ\nFugVK7rerZShIRoRiYuVK60S37v38PcrVYLFi23JXinbEI1HJw+JSDzt3w/XX1803MFWdFS4R4cC\nXkRcd999sHRp0fdvvdV2WpLo0BCNiLgmL89uZLrhhqLH6tWDr7+2IRopoFk0IpKwNm6EuXPt8dFH\nsH170TblytnFVoV7dCngRSTqtmyBZ5+Fd9+1TauL8/TT3lwrJtY0RCMiUXPwIDz3nG1q/euvkX3N\nFVfYL4JUu0EpWmI9i2YCsA1wWowzE9gFLA48BpemIyKS3ObPtyp84MDIw71OHZg4UeEeK5EM0UwE\nngNeC9PmE6BDVHokIkll61bbRWny5MjaH3WUbXB96aW20uMJJ8S2f14WScB/CtQupo1+/4p4TF6e\nLew1eDDs3u3czuezyv6SS6BtWzj/fAt5ib1oXGT1A+cBS4HNwCAgzFa3IpLsdu6Em2+G995zbpO/\ni1JWFlSt6lrXJEg0Av5roBawF7gcmAnUD9UwKyvr9+eZmZlkZmZG4eNFxE3LlsFf/gJr1zq3ad7c\ndlHywr6n0ZadnU12dnZUvlekQyu1gXeBMyNoux5oBuwo9L5m0YgkuTffhO7dQy8xALbf6dCh0LOn\nd3dRirZ4r0VTLejDWwSeFw53EUliubkwaJDz+jFgob5qle2mpHBPDJEM0UwB2gBVgI3AEKBc4Ng4\n4FqgL5CLDdNcH/1uiki8bN8O110HH38c+vgxx8Drr0PHju72S4qnG51ExNGMGXDHHZCTE/p4w4bW\npkEDd/vlJfEeohGRFLNlC1x9tT2cwv2aa2wvVYV74lLAi8jv8vJsPfb8yjyUI46A4cPh7bdteEYS\nlxYbExHAdljq1Qs++8y5TUaGrfrYtq17/ZLSUwUv4nH798Njj9km1uHC/YorbCs9hXvyUAUv4mFz\n58Ltt4e/aalqVRg9Gjp31qJgyUYVvIgH5eTYrkrt2oUP927dYMUKmyapcE8+quBFPOTQIVtCoLgF\nwurUgfHj4eKL3eubRJ8CXsQjVq6Em26CRYuc26Sn2x2rjz4KRx/tXt8kNhTwIh7wj39A166wZ49z\nmwsugBdegEaNXOuWxJjG4EVSWG4uPPCA3ZTkFO5VqtiuSp98onBPNargRVLU9u12IfWjj5zb9Opl\nqz9mZLjXL3GPAl4kBX31lVXtGzeGPl6vHrz6Kpx3nqvdEpdpiEYkhfj98NJL0Lq1c7h37Gi/ABTu\nqU8BL5Iifv7Zbkbq3RsOHCh63OeDJ5+0C67HHed+/8R9GqIRSQFz59osmS1bQh+vXBmmTLEbm8Q7\nVMGLJLH9+6F/fwtup3Bv1szmvivcvUcVvEiSWroUunSB5cud23TvDmPHwlFHudcvSRyq4EWSzKFD\n8Mwz0KKFc7ifcIJtkP3KKwp3L1MFL5JEVq60BcC+/NK5zcUXw6RJULOme/2SxKQKXiQJ5FftZ5/t\nHO7ly8PIkXbBVeEuoApeJOFFUrU3bgxvvAFnneVevyTxqYIXSVAHDxZftft8Novmq68U7lKUKniR\nBOP3w7Rp8PDDsGaNc7u6dW2RsNat3eubJBdV8CIJZN48mx3TubNzuPt80K+fTZNUuEs4quBFEsDi\nxbas77/+Fb6dqnYpCVXwInH05ZdWrTdtGj7cVbVLaaiCF3FZbi5Mnw7PPgsLFhTfvm1bGDbMlhwQ\nKQkFvEiU/PCDBfHq1XDiiVCr1uGPjAx46y0YM8Z5Kd9gTZva97vkktj3XVKTAl6kjPx+ePllGDgw\n/J6nkapTx5b17dQJjtAgqpSBAl6kDDZuhJ49i784Gonq1WHwYNtGr3z5sn8/EQW8SCn4/TBhAgwY\nALt3l+17nXOO3ax03XUKdokuBbxICW3aZFX2+++X/nv4fLZ1Xv/+cMEF9lok2hTwIiWwbBm0aQM7\ndji36datYE/U/MemTbB5s118bd8e7r7bxtpFYsnNusHv9/td/DiR6MrNtXVhnNZgP+kk2/C6fXt3\n+yWpzWf/vCtVVusavUiEJk1yDvdbb7XqXuEuiUQVvEgE9u6FevWK7ntavTqMHw9XXRWffknqK0sF\nrzF4kQiMGhV6U+vp0+G889zvj0gkVMGLFGP7drsgWvgmpquvtoAXiSWNwYvE0OOPFw33tDQYOjQ+\n/RGJlAJeJIy1a+GFF4q+37s31K/vfn9ESkJDNCJhXHedLRAWrGJFWLcOqlWLT5/EWzREIxIDCxYU\nDXeA++5TuEtyUAUvEoLfD5mZMH/+4e9Xr25b6VWqFJduiQepgheJsvfeKxruAFlZCndJHqrgxdNW\nrbIgz8mxee45OfZYubLozJkGDexu1XTdPSIu0o1OIiW0e7et4/7225F/zfDhCndJLqrgxXOWL7eb\nlFavjvxrzj8fPv1Uy/qK+2I9Bj8B2AZ8E6bNaGANsBQ4pzQdEXHD1KnQokXJwj0tDUaMULhL8okk\n4CcCl4U53h6oC9QDegMhbgsRia+DB6FfP7jhBls4LFJ168LkydCqVez6JhIrkYwofgrUDnO8AzAp\n8HwBcDxQDav6ReIuJ8c2sP7889DHK1WyO1Nr1bI13f/wh4JHxYru9lUkmqJxyagGsDHo9SagJgp4\nSQCzZ0OPHrDN4aexYUP4xz/g9NPd7ZeIG6I1J6Dw6GTIq6lZWVm/P8/MzCQzMzNKHy9yuB074J57\nbHjFSadO8MorcMwx7vVLpDjZ2dlkZ2dH5XtFetmoNvAucGaIYy8C2cDUwOuVQBuKVvCaRSOumDED\n+vZ1rtrT0uCZZ2xMXhdOJdHF+07WWcAtgeetgJ1oeEbiYPt2Wxzs6qudw71aNZg3D/r3V7hL6otk\niGYKVpFXwcbahwDlAsfGAXOwmTRrgf8B3aLfTZHwpk2zqv2nn5zbtGsHEyfahVQRL9CNTpL0Ro2y\nitzJccfBs89C166q2iX5aKkC8ayRI2HQIOfjV10FL76oql28SatJStJ6+mnncK9cGd54A955R+Eu\n3qUKXpLS0KHw0EOhj3XoAOPHa1MOEVXwknSeeMI53Lt2tRuXFO4iCnhJMn/9KzzySOhj3bvbjUtp\nae72SSRRaYhGksLOnTB4MIwdG/p4r152MfUIlSwiv9NfB0loubnwwgtQr55zuN92m8JdJBRV8JKw\n5s61+e3Llzu36dsXxoxRuIuEor8WknBWrbL56+3ahQ/3O++0ql7hLhKa/mpIwjh0CJ56Cho3hvfe\nc2539NG2P+ro0bozVSQcDdFIQvj+e7j5Ztv3NJybbrI58DVrutMvkWSmCl7ibupUaNIkfLi3agVf\nfgmvv65wF4mUAl7iZvduuPVW2yd1167QbWrVgilT4N//hpYt3e2fSLLTEI3ExRdfQJcusH596OM+\nH9x7LwwZAhUquNs3kVShgBdXHThgd6MOHQp5eaHb1KxpQzHa0VGkbBTw4pply+xC6pIlzm06dYJx\n4+CEE9zrl0iq0hi8xNyhQzBiBDRr5hzulSrZbktvvqlwF4kWVfASU999Zys8hpsh07IlTJ4Mdeu6\n1i0RT1AFLzGRl2frw4Sb/pieDo89ZscV7iLRpwpeom7pUujTx+atO2nYEF57Dc49171+iXiNKniJ\nml9/tS30mjULH+79+8OiRQp3kVhTBS9RMWuWLf61caNzm5NPhldfhYsucq1bIp6mCl7K5Pvv4c9/\nho4dw4d79+7wzTcKdxE3qYKXUtm3D55+GoYNg/37ndudfrpdbG3Txr2+iYhRwEuJ+P22qfXAgVa9\nOznqKNtib9AgOPJI9/onIgUU8BKx5cvh7rth3rzw7dq1g+efhzp13OmXiISmMXgp1po1cPvtNqc9\nXLhXr25L/77/vsJdJBGogpeQ/H7bE3X0aJgzx147KVcO+vWzIZljj3WvjyISngJeDvO//9lKjqNH\nw4oVxbe/9FL429+gQYPY901ESkYB71EHD9o6MatW2WP1avvvkiWwZ0/xX3/aaTBqFFx5pfZFFUlU\nCniPWbjQhlI+/NBWeSypihXhoYdgwACbKSMiiUsB7yGzZkHnzvDbbyX/2tq17U7V7t21nK9IslDA\ne8SkSdCjR8mr9sxMuOceuOoqSEuLSddEJEYU8B7wf/9nNyZFqmJFuP56uOsumxopIslJAZ/C/H54\n+GHb/zSUjAxo1Ajq17dZMA0a2PPTTrOpjyKS3BTwKerQIbs5afz40MfbtYPp022rPBFJTQr4FPTb\nb3DTTTBtWujjnTvbXPfy5d3tl4i4S0sVpJjcXBs/dwr3Pn3g739XuIt4gQI+heTlQa9eMHNm6OOD\nB9siYJoNI+INGqJJEX6/Lc376quhjz/7rK0XIyLeoYBPEU89ZSFemM8HEyZA166ud0lE4szNVUT8\n/nBLEkqpPf883HFH6GNjxjgfE5HE57PFnkqV1RqDT3J//7stIRDK448r3EW8TBV8Eps92za8zs0t\neqxfP7uDVSs9iiQ3VfAes2+f3aHqFO633gojRyrcRbxOFXyS+eADu0P1u+9CH+/Y0ebAp+vyuUhK\nUAXvAVu3wo03wmWXOYf7RRfZnqgKdxGByAL+MmAlsAa4P8TxTGAXsDjwGBytzondvDR+PDRsCFOm\nOLf74x/hnXe0CYeIFCiu1ksDxgBtgc3AV8AsoPBunZ8AHaLeO4/butXWlPnoI+c25cvDgw/CAw8o\n3EXkcMUFfAtgLbAh8Hoq0JGiAa/LeVE2b54NyWzb5twmMxNefFEbXotIaMUN0dQANga93hR4L5gf\nOA9YCswBzoha7zzo0CEYMgTatnUO94wMW5Jg3jyFu4g4K66Cj2Tay9dALWAvcDkwE6gfqmFWVtbv\nzzMzM8nMzIykj56RkwNdusDHHzu36doVnnkGqlRxrVsi4qLs7Gyys7Oj8r2KG1ppBWRhF1oBHgTy\ngOFhvmY90AzYUeh9TZMMY+5cG2//8cfQx6tWtTXc27Vzt18iEl+xnCa5EKgH1AbKA9dhF1mDVQv6\n8BaB54XDXRz89BP07g2XXuoc7hddBEuWKNxFpGSKG6LJBe4EPsBm1LyCXWC9LXB8HHAt0DfQdi9w\nfUx6mmIOHYJx42yN9l9+Cd3G54NHH4VHHtEa7iJScrqTNQ4+/9wWCFuyxLlNtWq2kNjFF7vXLxFJ\nPLqTNUnk5MAtt0Dr1uHD/U9/suMKdxEpCwW8C3bssJuR6ta1C6VOjj0WRo2y9WaqV3evfyKSmrRq\nSQzt2WOBPWIE7N4dvm23bjB0qA3NiIhEgwI+Bvbts12Whg2zWTLhNG0KY8dCq1bu9E1EvEMBH2XT\np8Pdd8OWLeHbVa5sFXuPHpohIyKxoTH4KMnLsymN114bPtyPPBL694fVq23+u8JdRGJFFXwU7N1r\nuyhNm+bcJj3dqvXBg6FmTff6JiLepYAvo82bbRelRYtCH/f5bAmCIUOgTh13+yYi3qaAL4NFi6BD\nB+chmcsvt4XBGjVyt18iIqAx+FKbNg0uuMA53B99FN57T+EuIvGjCr6Edu6Ehx+2aZChHHkkTJwI\nN9zgbr9ERApTwEfI77e1YQYOdN6Io3p1mDkTWrZ0t28iIqEo4COwahXcfrvtoOTk7LNh1iyoVcu9\nfomIhKMx+DD27bOles86K3y4/+Uv8NlnCncRSSyq4B18/rnNbV+3zrlNxYrw+ONwzz1whH5VikiC\nUSwVcuCArfx44YXhw/2aa2DlSrsrVeEuIolIFXyQZcvspqSlS53bnHoqjBkD7du71y8RkdJQ7Ylt\nnzdiBDRr5hzu5crZ9MhlyxTuIpIcPF/Bb9hgY+3z5zu3adXK5raffrpr3RIRKTPPVvC5uVa1N2rk\nHO7p6fDkk/Dppwp3EUk+nqzg//MfuO228PuinnGGba/XtKl7/RIRiSZPVfC7d9tmHK1aOYe7zwcD\nBthCYgp3EUlmnqjg/X6YMQPuuiv8ZhwnnwyTJkFmpmtdExGJmZSv4Ldvh06dbN66U7j7fBb+33yj\ncBeR1JHSFfz06dC3r4W8kyZNYPx4aNHCvX6JiLghJSv4n3+GG2+0/VGdwr1CBduMY+FChbuIpKaU\nq+BnzbLNrJ2W9AXbaen556F2bde6JSLiupSp4HftshuWOnZ0DvfKlW1N99mzFe4ikvpSooKfPx9u\nvhl++MG5TYcOMG6cbcohIuIFSV3B56/8mJnpHO7HHw+vvWY7LSncRcRLkraCX7ECunSBxYud21x+\nObz0EtSo4V6/REQSRdJV8H4/jB1rd5k6hfsxx8DLL9tYu8JdRLwqqSr4LVugZ0/45z+d25x/vq0h\nc+qp7vVLRCQRJUUF7/fb7JfGjZ3DPT0dnngCsrMV7iIikAQV/Pbtdjfq9OnOberXh8mToXlz9/ol\nIpLoErqCnzHD1msPF+59+sDXXyvcRUQKS8gK/pdfbFnfyZOd21SrZhdSr7zSvX6JiCSThAv499+H\nHj3CL+t7/fW28XVGhnv9EhFJNj4XP8vv9/vDNtixA045BX79NfTxjAxbQ6Zz5xj0TkQkAfl8Pihl\nVifUGHzlyjByZOhjV10Fy5Yp3EVEIpVQFbw1sjtQP/jAXh97LIweDbfcYhtziIh4SVkq+IQLeIBN\nm2zOe/PmMGEC1KoV456JiCSolAt4gNWroV49Ve0i4m0pGfAiIpJCF1lFRCR6FPAiIikqkoC/DFgJ\nrAHud2gzOnB8KXBOdLomIiJlUVzApwFjsJA/A7gBaFioTXugLlAP6A28EOU+ppzs7Ox4dyFh6FwU\n0LkooHMRHcUFfAtgLbABOAhMBToWatMBmBR4vgA4HqgWvS6mHv3wFtC5KKBzUUDnIjqKC/gawMag\n15sC7xXXpmbZuyYiImVRXMBHOq+x8BQezYcUEYmz4uZWtgKysDF4gAeBPGB4UJsXgWxs+Absgmwb\nYFuh77UWqFP6roqIeNI67Dpn1KUHvnltoDywhNAXWecEnrcCvoxFR0REJPouB1ZhFfiDgfduCzzy\njQkcXwo0dbV3IiIiIiISXZHcKJVKJmDXH74Jeq8yMBdYDfwLm0qa70Hs3KwE2rnUR7fUAj4GlgPL\ngLsD73vxfByFTSNeAnwLDA2878VzkS8NWAy8G3jt1XOxAfgvdi7+E3gvKc5FGjZ0UxsoR+gx/FRz\nAXY3b3DAPw3cF3h+PzAs8PwM7JyUw87RWlJr+YjqwNmB55Wwob6GePd8VAj8Nx27VtUa754LgAHA\nG8CswGuvnov1WKAHS4pz8Ufg/aDXDwQeqa42hwf8Sgpu/qoeeA32mzj4XzXvYxeqU9VMoC06HxWA\nr4BGePdc1AQ+BC6ioIL36rlYDxTeYToq5yLWyR/JjVJeUI2CaaPbKPgfdxJ2TvKl8vmpjf3LZgHe\nPR9HYNXXNgqGrrx6Lp4F7sWmXefz6rnwY7/sFgK9Au9F5VykR7WbRemGp6L8hD8vqXjOKgHTgXuA\nPYWOeel85GFDVscBH2DVazCvnIsrgR+xMedMhzZeORcA5wM5wInYuPvKQsdLfS5iXcFvxi605avF\n4b99vGIb9s8sgD9gP9xQ9PzUDLyXSsph4f46NkQD3j4fALuA2UAzvHkuzsPWsFoPTAEuxn4+vHgu\nwMIdYDswA1sDLCnORSQ3SqWi2hS9yJo/bvYARS+YlAdOxc5VKm1S6ANew/45HsyL56MKBTMhjgbm\nA3/Cm+ciWBsKxuC9eC4qAMcEnlcEPsdmxiTNuQh1o1QqmwJsAQ5g1x+6YVfIPyT0lKeHsHOzErjU\n1Z7GXmtsWGIJ9s/xxdi0WS+ejzOBr7Fz8V9s/Bm8eS6CtaFgFo0Xz8Wp2M/EEmwqcX5GevFciIiI\niIiIiIiIiIiIiIiIiIiIiIiIiIiIiEgq+n/KE1dZUTxCmAAAAABJRU5ErkJggg==\n",
       "text": [
        "<matplotlib.figure.Figure at 0x7f71f804d110>"
       ]
      }
     ],
     "prompt_number": 73
    },
    {
     "cell_type": "raw",
     "metadata": {},
     "source": [
      "1149\tstatic int\n",
      "1150\tstring_contains(PyObject *str_obj, PyObject *sub_obj)\n",
      "1151\t{\n",
      "1152\t    if (!PyString_CheckExact(sub_obj)) {\n",
      "1153\t#ifdef Py_USING_UNICODE\n",
      "1154\t        if (PyUnicode_Check(sub_obj))\n",
      "1155\t            return PyUnicode_Contains(str_obj, sub_obj);\n",
      "1156\t#endif\n",
      "1157\t        if (!PyString_Check(sub_obj)) {\n",
      "1158\t            PyErr_Format(PyExc_TypeError,\n",
      "1159\t                \"'in <string>' requires string as left operand, \"\n",
      "1160\t                \"not %.200s\", Py_TYPE(sub_obj)->tp_name);\n",
      "1161\t            return -1;\n",
      "1162\t        }\n",
      "1163\t    }\n",
      "1164\t\n",
      "1165\t    return stringlib_contains_obj(str_obj, sub_obj);\n",
      "1166\t}"
     ]
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Fast search/count implementation, based on a mix between Boyer-Moore and Horspool, with a few more bells and whistles on the top. For some more background, see: http://effbot.org/zone/stringlib.htm.\n",
      "\n",
      "[...]\n",
      "\n",
      "When designing the new algorithm, I used the following constraints:\n",
      "\n",
      "* should be faster than the current brute-force algorithm for all test cases (based on real-life code), including Jim Hugunin\u2019s worst-case test\n",
      "* small setup overhead; no dynamic allocation in the fast path (`O(m)` for speed, `O(1)` for storage)\n",
      "* sublinear search behaviour in good cases (`O(n/m)`)\n",
      "* no worse than the current algorithm in worst case (`O(nm)`)\n",
      "* should work well for both 8-bit strings and 16-bit or 32-bit Unicode strings (no `O(\u03c3)` dependencies)\n",
      "* many real-life searches should be good, very few should be worst case reasonably simple implementation"
     ]
    }
   ],
   "metadata": {}
  }
 ]
}