yelp_analysis.ipynb

{
 "metadata": {
  "name": "",
  "signature": "sha256:443772b49ca2108d210f189d78f582b349664dfe5603a2887e549cce5e4f57b3"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import pandas as pd\n",
      "from matplotlib import pyplot as plt\n",
      "%matplotlib inline\n",
      "from matplotlib import rcParams\n",
      "rcParams['figure.figsize'] = (10, 6)\n",
      "rcParams['figure.dpi'] = 150"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 50
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "yelp = pd.read_csv('yelp.csv')\n",
      "cats = pd.read_csv('categories.csv')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 3
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "yelp.head(1)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "html": [
        "<div style=\"max-height:1000px;max-width:1500px;overflow:auto;\">\n",
        "<table border=\"1\" class=\"dataframe\">\n",
        "  <thead>\n",
        "    <tr style=\"text-align: right;\">\n",
        "      <th></th>\n",
        "      <th>Unnamed: 0</th>\n",
        "      <th>address</th>\n",
        "      <th>categories</th>\n",
        "      <th>city</th>\n",
        "      <th>display_address</th>\n",
        "      <th>id</th>\n",
        "      <th>image_url</th>\n",
        "      <th>name</th>\n",
        "      <th>phone</th>\n",
        "      <th>postal_code</th>\n",
        "      <th>rating</th>\n",
        "      <th>rating_img_url_large</th>\n",
        "      <th>review_count</th>\n",
        "      <th>url</th>\n",
        "      <th>complete_lookup_address</th>\n",
        "      <th>latitude</th>\n",
        "      <th>longitude</th>\n",
        "      <th>latlng_confidence</th>\n",
        "    </tr>\n",
        "  </thead>\n",
        "  <tbody>\n",
        "    <tr>\n",
        "      <th>0</th>\n",
        "      <td> 0</td>\n",
        "      <td> [u'3201 California St']</td>\n",
        "      <td> [[u'Banks &amp; Credit Unions', u'banks'], [u'Inve...</td>\n",
        "      <td> San Francisco</td>\n",
        "      <td> 3201 California St, Laurel Heights, San Franci...</td>\n",
        "      <td> sf-fire-credit-union-san-francisco</td>\n",
        "      <td> http://s3-media2.ak.yelpcdn.com/bphoto/t31F2Pk...</td>\n",
        "      <td> SF Fire Credit Union</td>\n",
        "      <td> 4156744800</td>\n",
        "      <td> 94118</td>\n",
        "      <td> 4.5</td>\n",
        "      <td> http://s3-media4.ak.yelpcdn.com/assets/2/www/i...</td>\n",
        "      <td> 184</td>\n",
        "      <td> http://www.yelp.com/biz/sf-fire-credit-union-s...</td>\n",
        "      <td> 3201 California St, 94118</td>\n",
        "      <td> 37.787171</td>\n",
        "      <td>-122.447225</td>\n",
        "      <td> P1AAA</td>\n",
        "    </tr>\n",
        "  </tbody>\n",
        "</table>\n",
        "<p>1 rows \u00d7 18 columns</p>\n",
        "</div>"
       ],
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 4,
       "text": [
        "   Unnamed: 0                  address  \\\n",
        "0           0  [u'3201 California St']   \n",
        "\n",
        "                                          categories           city  \\\n",
        "0  [[u'Banks & Credit Unions', u'banks'], [u'Inve...  San Francisco   \n",
        "\n",
        "                                     display_address  \\\n",
        "0  3201 California St, Laurel Heights, San Franci...   \n",
        "\n",
        "                                   id  \\\n",
        "0  sf-fire-credit-union-san-francisco   \n",
        "\n",
        "                                           image_url                  name  \\\n",
        "0  http://s3-media2.ak.yelpcdn.com/bphoto/t31F2Pk...  SF Fire Credit Union   \n",
        "\n",
        "        phone  postal_code  rating  \\\n",
        "0  4156744800        94118     4.5   \n",
        "\n",
        "                                rating_img_url_large  review_count  \\\n",
        "0  http://s3-media4.ak.yelpcdn.com/assets/2/www/i...           184   \n",
        "\n",
        "                                                 url  \\\n",
        "0  http://www.yelp.com/biz/sf-fire-credit-union-s...   \n",
        "\n",
        "     complete_lookup_address   latitude   longitude latlng_confidence  \n",
        "0  3201 California St, 94118  37.787171 -122.447225             P1AAA  \n",
        "\n",
        "[1 rows x 18 columns]"
       ]
      }
     ],
     "prompt_number": 4
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "cats.head(1)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "html": [
        "<div style=\"max-height:1000px;max-width:1500px;overflow:auto;\">\n",
        "<table border=\"1\" class=\"dataframe\">\n",
        "  <thead>\n",
        "    <tr style=\"text-align: right;\">\n",
        "      <th></th>\n",
        "      <th>Unnamed: 0</th>\n",
        "      <th>id</th>\n",
        "      <th>short_category</th>\n",
        "      <th>full_category</th>\n",
        "    </tr>\n",
        "  </thead>\n",
        "  <tbody>\n",
        "    <tr>\n",
        "      <th>0</th>\n",
        "      <td> 0</td>\n",
        "      <td> sf-fire-credit-union-san-francisco</td>\n",
        "      <td> banks</td>\n",
        "      <td> Banks &amp; Credit Unions</td>\n",
        "    </tr>\n",
        "  </tbody>\n",
        "</table>\n",
        "<p>1 rows \u00d7 4 columns</p>\n",
        "</div>"
       ],
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 5,
       "text": [
        "   Unnamed: 0                                  id short_category  \\\n",
        "0           0  sf-fire-credit-union-san-francisco          banks   \n",
        "\n",
        "           full_category  \n",
        "0  Banks & Credit Unions  \n",
        "\n",
        "[1 rows x 4 columns]"
       ]
      }
     ],
     "prompt_number": 5
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "merged = pd.merge(yelp, cats, on='id')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 8
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "'''most popular categories'''\n",
      "grouped = merged.groupby('short_category').id.count()\n",
      "grouped.sort(ascending=False)\n",
      "print grouped.head(20)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "short_category\n",
        "hair                  2856\n",
        "contractors           2162\n",
        "coffee                1831\n",
        "realestateagents      1675\n",
        "othersalons           1548\n",
        "grocery               1546\n",
        "massage               1528\n",
        "drycleaninglaundry    1427\n",
        "skincare              1326\n",
        "sandwiches            1316\n",
        "mexican               1306\n",
        "womenscloth           1294\n",
        "apartments            1147\n",
        "banks                 1075\n",
        "chinese               1060\n",
        "jewelry               1032\n",
        "realestatesvcs         971\n",
        "pizza                  880\n",
        "bakeries               854\n",
        "beer_and_wine          833\n",
        "dtype: int64\n"
       ]
      }
     ],
     "prompt_number": 23
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# '''Making a map for coffee'''\n",
      "category = 'coffee'\n",
      "mean_lat = all_locations.latitude.mean()\n",
      "mean_long = all_locations.longitude.mean()\n",
      "\n",
      "all_locations = merged[merged.short_category == category]\n",
      "lower_perc, upper_perc = all_locations.rating.describe()['25%'], all_locations.rating.describe()['75%']\n",
      "\n",
      "good = all_locations[all_locations.rating >= upper_perc]\n",
      "bad = all_locations[all_locations.rating <= lower_perc]\n",
      "ax = plt.subplot()\n",
      "plt.plot(good.longitude - mean_long, good.latitude - mean_lat, '*', color = 'red', alpha = .5)\n",
      "plt.plot(bad.longitude - mean_long, bad.latitude - mean_lat, '*', color = 'blue', alpha = 0.5)\n",
      "ax.set_aspect(1)\n",
      "plt.ylim(0.02,0.1)\n",
      "plt.xlim(-0.15, -0.09)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 158,
       "text": [
        "(-0.15, -0.09)"
       ]
      },
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAATEAAAFwCAYAAADQYUyhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXlYVPX+x18zw7CquLEo4AqyiaYCllmipoilXNNKrVzS\n8lqJdqtbt9ui1TXNvNfM6qb9Sropdq9ZbgiulLu4ooMouwiyyjrAMMyc3x84E+CwgzDjeT2PzyMz\nZ/nOmTnv8/l+vp9FIgiCgIiIiIiRIm3vAYiIiIi0BFHEREREjBpRxERERIwaUcRERESMGlHERERE\njBpRxERERIyaBkUsIiICDw8P3NzcWL169V3vx8XF8dBDD2FpacnatWubtK+IiIhIS5HUFyem0Whw\nd3fn4MGDODk54efnR1hYGJ6envptcnJySE1N5ddff6Vbt268/vrrjd5XREREpKXUa4mdOXMGV1dX\n+vXrh1wuZ+bMmezcubPGNnZ2dvj6+iKXy5u8r4iIiEhLqVfE0tPTcXFx0f/t7OxMenp6ow7ckn1F\nREREGku9IiaRSJp94JbsKyIiItJYzOp708nJibS0NP3faWlpODs7N+rAjd3X1dWVxMTExo5XRETE\nSBk4cCAJCQmtftx6LTFfX1/i4+NJSUmhoqKCn376ialTpxrctvb6QGP3TUxMRBCEDv/vgw8+aPcx\niOMUx2nM42wrY6VeS8zMzIwNGzYQGBiIRqNhwYIFeHp68s033wCwaNEiMjMz8fPzo6ioCKlUyuef\nf05sbCydOnUyuK+IiIhIa1KviAEEBQURFBRU47VFixbp/+/o6Fhj2tjQviIiIiKtiRix30gCAgLa\newiNQhxn6yKOs+NTb7DrPRmAREI7D0FEROQe0Fb3umiJiYiIGDWiiImIiBg1ooiJiIgYNaKIiYiI\nGDWiiImIiBg1ooiJiIgYNaKIiYiIGDWiiImIiBg1ooiJiIgYNaKIiYiIGDWiiImIiBg1ooiJiIgY\nNaKIiYiIGDWiiImIiBg1ooiJiIgYNaKIiYiIGDWiiImIiBg1ooiJiIgYNaKIiYiIGDWiiImIiBg1\nooiJiIgYNaKIiYiIGDWiiImIiBg1ooiJiIgYNaKIiYiIGDWiiImIiBg1ooiJiIgYNaKIiYiIGDWi\niImIiBg1oojdJwiCQETEcQRBaO+hiIi0KqKIGSlarZa1H69Dq9U2avvY2ER++aWY2NjENh6ZiMi9\npUERi4iIwMPDAzc3N1avXm1wm5CQENzc3Bg6dCgXLlzQv/7555/j4+PD4MGD+fzzz1tv1CJE7o7k\n0JZsIndH1rtddHQMK1ZsJSwsFUfHQMLCUlmxYivR0TH3aKQiIm1LvSKm0Wh49dVXiYiIIDY2lrCw\nMK5evVpjm/DwcBISEoiPj2fjxo0sXrwYgCtXrvDtt98SHR3NpUuX2LNnD4mJohXQUrZuDmP62Bf5\ncc1xfO0e48c1x5k+9kW2bg4zuL2vrw/Tp/uhUglIJBJUKoEZM/zx9fW5xyMXEWkb6hWxM2fO4Orq\nSr9+/ZDL5cycOZOdO3fW2GbXrl3MnTsXgJEjR1JQUEBmZiZXr15l5MiRWFpaIpPJGDNmDDt27Gi7\nT3KfMHPOM7wQMg21CqRSKWoVLFz2JDPnPGNwe4lEgkQioaioEkGIpKioUv+6iIgpUK+Ipaen4+Li\nov/b2dmZ9PT0BrfJyMjAx8eHo0ePcvv2bUpLS9m7dy83b95s5eHff0ilUqRSCSUlWuKFCEpKtPrX\n6yInp5CQEFeWL59ISIgrublFLR6HIAgcj4gQFwpE2h2z+t5s7NPa0A/Zw8ODt956i4kTJ2JjY8Ow\nYcPqvdFE7kar1fKvlet57Z2QGtfuRnIGS1Y9QuCUQCJ3R3IzOaPe4wQEjND/39vbtVXGlhgbS/Ev\nv5Do4oKrt3erHFNEpDnUK2JOTk6kpaXp/05LS8PZ2bnebW7evImTkxMAL7zwAi+88AIA77zzDn36\n9DF4nuXLl+v/HxAQQEBAQJM+hKmic957+UQSFBykf33RsgX6/1d/vSUIgkBk5AkCA0fV+/CKiY7m\nSng4Dmo1gY6OHA4L44xczuDJkxni59dm5xUxPqKiooiKimr7Ewn1oFarhQEDBgjJycmCSqUShg4d\nKsTGxtbYZu/evUJQUJAgCIJw8uRJYeTIkfr3srKyBEEQhNTUVMHDw0MoLCy86xwNDOG+ZMv3W4Un\nAxYKsx/+u/DeI4eE2Q//XXgyYKGw5futbXbOK1fihZde2idcuRJf73ZarVa4fvmycOCNNwThgw+E\nA2+8IcRfuSJotdo2Pa+I8dNW93q9lpiZmRkbNmwgMDAQjUbDggUL8PT05JtvvgFg0aJFTJ48mfDw\ncFxdXbGxseH777/X7z9jxgzy8vKQy+V89dVXdOnSpS312GSYOecZbG1tCV15AqmNFHUZLHzzSQKn\nBLb6uaKjYwgPv4Ja7XAnBOMwcvkZJk8ejJ/fkLu21y0UVBYVEWltjVBUpH+9Lc8rIlIXkjsK2X4D\nkEhE57AB9u3cxxdvH6WbvZb8bClLVj3SalPH6giCgEKRQGhoKjY2j6FUHmTevH54eQ2sU5jORUVh\na2fHQC8vEmNjKczJIadc3qQpYXPOK2LctNW9brKedq1Wy8cff9PoiPaOhs55/58jK1my6pEGnffN\npTkhGCMCAnD19kYikeDq7Y2lnXOTswHE0A+R1sJkLbGdOw/x9tuJrFo1kODg8a1+fFMiKuocdna2\neHkNJDY2kdzcIsaMGd7gftWnhDLZODSaw8jlWY2eEjb3vCLGSVvd6yYnYps3/0xoaDRqtStS6Qto\ntd8hlycwd64f8+ZNb7XziIhTQpGm0VYiVq9j3xiZM2catra2rFyZTKWqADMLgTffnMCUKWPbe2gm\nR/UpobV1JEVFgv51EZF7hcmJmC6i/XZeOdZlKyiyctW/LtL66LIBqk8JRUTuJSY3nYyJjmbNe1/g\n1tmMv3s48Y+4dOKLNbz50avNCsY0JQQxsFSkHRFXJxuJj68v73/2V0b164FMJmNUvx58sPav+Pj6\ntvfQ2p2m1hQTTCg/0thXq0XqxuSmk60VjGlKNDew1JTyI3fvPsKWLRJ8fI6Iq9UmhslNJ+HuYMyi\n3FyGjxnTqudoKoIgcCIyklGBgfdcUJu6ilg9P3KcTMZhjYasFuRHtifianXHQZxONoHawZjtLWAA\nCQoF19evJ0GhuOfnbmpgqY+vL37TpyOoVFU/PJUK/xkz8PH1Nbop5pw50wgJeQzVnfprKhUsWzaB\nOXOmtffQRFoJk5tOdjR0Vo08PZ3RmZlcXL+eaCene27VNGUVsb4peYJCYVRTzOr116ytv6SkRK5/\nXcQ0EEWsjREEgaysLLrExODauTO/x8RQbGZ2zy2ZptYUK8zJwTUkRD8lP3f0KGe2b2+1Ejz3kuTk\nLFatGsgTTwTw4Ydfk5SU3d5DEmlFxMdRGzPEz48hEycilJWxHxDKyvAJDOzwN37tKfnTixYZnGJ6\nDx/ONx9/3KFX/ZYtm01w8Hji4pLJzHRl4sQ/rr0gtrIzekRL7B7w++8X6OvmhuDlhezqVZT5+Ua3\nWlrXFDNqzx4kW7ZwxMeH8cHB7TxKw9S3Omttbc0vvxTj4pLYalVvRe4tJrk62ZFQKBJ47+3feXFh\nLyZNndRhVkubQ/VV369XrkTxv/8xrFMnXpBK+U6rJd7MDPmIMXz02fIOJdKGVmdHjFARF1dAZaVj\ns5LXRZqOuDppZFTv9+gzfD7Ho8358MMw8ks1RilgUHOK+ee//Y1pH3wAKlWVk1ylYtD0WWQXj2xx\ng97WnuIZWp0dPHgQM2b4i63sTIAOIWK6H61Wq+W3vXtZMG4clZWV7T2sFmHq/R6lUikSqRRtSQl/\nLSpnV4oZh3/X0KtXUJMa9BoK2WiLbuW1Oz7l5RWL9cxMhA7hE9P9aAVVJDHL30W4dInP33uP1z/5\npL2H1mzuhwoPWcnJDFy1igVPPMH3X37Lzl/z9YI9a5Y/Xl4DGzxG9ayA/FJNm5WsNrQ6GxV1Tkxe\nNwE6hE9s4fwfOBOuxCxnO17CDZyEG2SYacixsGDIrFms3rSpPYfYbO6non8KRQLr1yfg6CghM1Mg\nJMS1Xke5oayATDMzLN2HcOp8lzarTyYmwbcfJl1PzLa7I56jHEjaf4BxylScJGoiNKDt04en77R8\nM0baot9jR6UxwbTVBcTH1xcrKytSQ0OR2NggqFSMnDWLcsGc/VGJbWa96qx+cTXSdOgQIlZcrKGn\n/UUSzLpwEvCXSJBotXTWaOjSuXN7D0+kETRGsGNjE9mxo4ic1M0899I8gyEbublFbTLFE7srmS4d\nYjp55Uo8ZdlpfBDyJj3zlfR2sKYkKYkuXbow7vnnjTb5uDm0Z6J4W1FdQArzBpNzajUOo/owwseW\nUY/435NEfbGUdvtj0tNJb29X8Hblwy8+o0vPngiCwOl//hNvQeC2SoX/rFkM9PJq72HeE0yp/I0O\nX18fbiQk8N0Xx+hhJsHZqguelUnIMrtTWjpYnxXQltwPCy33Kx1CxHSMCAgAqio+2Jmbk+PoiJCZ\nCZj+j626o9vYchMbQiKR4O7jQ2encoT0vRSruiFUVOA/YwaZN25wbN8+Hp40qc2/Y7GUtmnSoURM\nR+3k46LcXP17prq6ZMjRbUoWaG5uEXOf6wZ7i7kh03D7ZhE3EhNJ/c9/yMzKwuGO5dmWU+n7aaHl\nfqJDipjOIgPummaY6upS9dzECCsr4hQK+gqCyQh1QMCIqrSlpUvprVQStno1e997j2FmZgzMyODH\nuXPR9O+Pp1arFzQRkcbQIRz7jRlCdeewVDqW1OubGODZicmTfUxmdUmXmygIAuc//BDLoCCC589v\n72G1OoIgsGfbNo6uXo1bbi5jKir4QSJB6+zM08OGkefoeF8t5twv3Pe5k9XTeApycymJjUeqzWHE\niMHtPbRWQ25jw5nt27mxbRtPe3rSKTGRrStWEBMd3d5Da1UkEgkePj70tLenSK0mXyajs5UVvV1c\nGOrkVKOSrIhIQxiNiEkkEq5fucL5w6couLwWSaXAwf/msebV11t8kwuCwLF9+9i371i7VtSoryy0\nqVGUm4vto4/i/OijJAYHI7ezQ5WfT4RWS3jUWbRabZtOpbVaLWs/Xteh66CJNA6jETGAbg4uTJnZ\nAzS36G5lhrnEl9RKf34Jv96oZOO6SIyNRfHdFkK/v9GqScdNpYZfTKvllCIVwYT8YtUZERCA7+jR\njFi+nFlff83Al1/Gbf58hBGjOJ83jMhdkfpt26Kuf+TuSA5tySZyd2SD24qC17ExGp+YjvgrVzi0\n/BMOZrhyO2cwA0ZX8tobfs0KWoyJjmbnxlASrpfirHXmBmMolvxObw8b5r84qV18bTq/WLlgzscf\nRjM9CJ6aP/uej+Nes3VzGD+HHsZS7cBA6TgStYcpl2cxfe44/P2GEP/55xSPGMFTL73ULFHXBRGn\nZN5mR+gRg+eZPW+WwX337dzHF28fZcmqRwgKDjLZFfK2xqSDXZtCUW4ufZ9/DsmPhaTcVCBN7gE0\nL47Mx9cXS0tLflu7gevXZUjQ0NPdk5dDRjRp5bM1f9Ram+5s2X4JtdoBd89ZXEo8TOyKrSafHjNz\nzjPY2toSuvIEUhsp6jIYH+yDNjmO1AQFrjIZv33zDesuXqSnpyfPLVnSpGutCyL2e/VVutp2JXTl\nCSTWEjKzC3j90yeZNDXwrn2qC6uv3WP8uOYw367bgf/4wSSluZvcCrmxYlTTSai6yX/6JQapkMWD\nPYbRo6eWNWt+btZ0UiKRIJVK0SqVKFRSNKpdlCi1+mldY2lu/StDxf90Cxjl5RqOHj1HebnGpOqQ\n1UX1rkTxQgQlJVr6ubnR1dOTE0eOcEOhYKyFBfKcHPJDQwn/6adGHTcmOpqtK1aQGhZGoKMjaT/9\nxPGffuJWjpKLyl9QZyu5lXbDYPejmXOe4YWQaajvtHvLul2OlUMvCku97uRfNr5umkjb0aCIRURE\n4OHhgZubG6tXrza4TUhICG5ubgwdOpQLFy7oX//kk0/w9vbGx8eH2bNno1KpWjxgc6EcD+tMhJu3\ncevSBeHmbTytMzEXypt1vMKcHDpN/hOfhs5keeifmR5EoyO5q1dvbc6P2pD46QT08uUbxMREcvny\nDf3rps6N5AyWrHqE/xxZySsrR/PfjZtx8/LCuVcvjhUUEJ6aivbaNRwdXTGPiWnUyq2hxZJi806M\nHpnHqxM7sXJiJUm79xs8Vm1hNRfkPDyqP2q1xCQLXRor9U4nNRoNr776KgcPHsTJyQk/Pz+mTp2K\np6enfpvw8HASEhKIj4/n9OnTLF68mFOnTpGSksKmTZu4evUqFhYWPPPMM2zbto25c+e2aMBD/PxI\nz8jj2PFTYLGf0jJzhgRObHY80YiAAEYE/PG32+DGh2z4+vpgZWV1J6m48cUA66uooFDEs25dJBLJ\nADw930YR8zrPPnuGZcsCTb5j9aJlC/T/t5SqGZUcQ9SuXbg+/zwOe/ZwvbCQWyfOcV7jymj7PALe\nWNJgRoOhBichf10CQGpoKKO9BlKuVNJvxgyDx9IJa+CUQCJ3R3L6ZAxFRb3E/MsORL0idubMGVxd\nXenXrx8AM2fOZOfOnTVEbNeuXXphGjlyJAUFBVV9Frt0QS6XU1paikwmo7S0FCcnpxYPWCKRkJ+v\nZMqgRPp6DqD31URu5w9plx9Sc5OK6xO/ESMGY2try8qVKRTl52OWL/Dk/H73TcfqnzdvJjo0FFe1\nmoV2dnwXEcH27dvp8uCjlNt5kda5J5VmIzh48RiXvzrAc/M0DfoKDaWxdenZ02Bz4NpUF9ag4CCs\nbO3vKnQp0r7UK2Lp6em4uLjo/3Z2dub06dMNbpOens7w4cN5/fXX6dOnD1ZWVgQGBvLYY4+1yqA9\n+vfEdvnfDOZW3muak1Rcn/hJJBJOHjlMUqyEHrK3sZI6kLnjF945eQC/uXOZPm9ei8bb0Uv9TJsz\nB1tbW1JWrkRqYwNlZcxYsYKAJ55g+w8/caC0F059ppF+o4LRvtpGTeVGBATc9bnPRUXVmZ9bH2L+\nZcejXhFr7I/c0LJpYmIi69atIyUlBVtbW5566im2bNnCs88+27yRVqO+3Mp7TXN/1PWJn2M/b0KW\npNDrwAUcyuVsLdIyfdkyxk6Z0uLxdvRSP9UbkHxpbY28pAQAmUzGYH9/Dp9KAPYjM++Kl79ro3+j\ntT93fb+hji70IjWpV8ScnJxIS0vT/52Wloazs3O929y8eRMnJyeioqIYNWoUPXpUhUA8+eSTnDhx\nwqCILV++XP//gIAAAqr9wEyV+sTvL395lkM7d5K480du2tszXlsIYHAFrbG0Rqmf5t7cTQ1B0TUg\nGfP447w49w06JVYtfDS1BLZEImnwcxsaW0cXemMhKiqKqKioNj9PvSLm6+tLfHw8KSkp9O7dm59+\n+omwsLAa20ydOpUNGzYwc+ZMTp06RdeuXXFwcMDd3Z2PPvqIsrIyLC0tOXjwIP7+/gbPU13ERKrQ\n3chjp0zhyO7d5CQnt+h4rVHqp7k3d1Mrj8xetgyAnTsPceq8F396umqhpLrwe3oOYNPKlTzyyAM1\nxL32uWp/bm15OWYuLgweMeKu7TWl+S0W+srKSp4NnMGWyO2YmRldGGarUtsgWbFiRZucp96rbGZm\nxoYNGwgMDESj0bBgwQI8PT355ptvAFi0aBGTJ08mPDwcV1dXbGxs+P777wF44IEHmDNnDr6+vkil\nUoYPH85LL73UJh/C1BAEgb4eHnqLZ3xwcIuPaWiVTvd6QzTXimtuXfvNm38mNDQalWoAqlJnVq9O\nZN26/YwabsnHd7qLH9m9G8mWLRzx8WF8cLDBc5mZnaZ7dy1B44bpP3dGairWN27wa1d7LseV1Rpb\nJoPdfRDOnWy20K95/2NuHLVgzfsf87eVyxu9n0jzMbq0o/uBBIWChPXrcQ0JadXpjC6lqak17QVB\nIEGhIDU0lMdsbDioVNJv3jwGennpRVCr1fKvlet57Z0QvWXU3Lr2Go2G5cu/YltYHnkZHvToHcfr\nS/vSJ+YEKf36cePgQVzVal6QSvlOqyVBLmfEnDmUqSAmwZlOnSagVB5k9GgJe/dW8NhD+dhYyyiI\ni8M8PZ2HVSpinZ05m19KXPFw+g18Xj82c6GcxC++QHKnqnBjv4PXXnyFI/87TU+tE6MYwwl+I1ea\nztinRvKvTV828hsybcS0o/uAti5R3dwFkcZYcbqEai+fSIKCg2rs19QQlM8+2cDGDdEoS/tja9aF\njJRKvvwwlLdH2eNmpyTD3p7Sy5eR2ttDWRkT3nwTlwEDOPLham5KJtPZ9hsuXMiksNCLXr1mEJP4\nh5XVpbiY5MxMBJWKoYETSdrLnQ7gVWMrys3FNSSEAZ6ebN4Uhl1OTqOu0ZqvP2d1txUc+PIKUrmU\nSrWUp195nKGjx5lsEn9HQRSxDkRHLlFdV8nwuvILdQnVTQlB0U0JVSoveveRkpZ4HdfuanJUWfSx\nkzNrqAeHy8rwGD0a5aVLvJWXR+/iYn7bupUxbm4Msu6KmeUpLBy74eXmyOnzKnr31sXhjcRcKOfo\n9q1cupjA0AfKsfBWsnTpAzWbG98ReoUigdPnezBytHOd49UhCAIHD55GIpOiUpvzm1kkKnUXsnKL\n2LWrjH79xBzLtkQUsQ5ES/xWbU1dVpyhxO2Fbz5J4JSqhOqmhKBUDwLu3as3CXEF0CmSThUyHuxl\nwQGJBKGoiLybN+nxyiuMuXyZ24MHQ0YGgkrF2P490SiV9Hv6acoFc+ISE/VW1rXLl0mK3IUEC8ps\np5NPHMLve/EZ3BeJRIKn5wD+tXI9lpZSIiJim+TH0y0O5CUomPrGYMZOncq7b37N0d81TJkp9rhs\na0QR62DU1ySlI1Ijv9A6gpISqf71piKRSFAornPkyFkK8oqxdxqLo1cat05sJc5qAO8uX074f/+L\nOiKCXnI5Yx0d+eHECaQ9epCQkYFmyJA6m/AePXqONOsh5MRl49ZpNOfSS7Eb7kC5YA78MR32HHyL\n6dP96k0l04VldOtmw/r1u9BoujFo0CvY2cmpkGchlcpY9/VfCQ1N1edYNiYdTaR5iCLWwTBk8XT0\n+lW18wtvJmc0+1j29g68++5oftmeQfo1GVmxyXwVPJzM7jaEffgh3kFBDHr9dVJDQylQKlFfvYr5\nU08xbvFiXL299cL/6KOPsGnlSjw938Hbu0rM9u2M4LtjGSSaQ0VpJZMmuhMfq2D1Wxv00+Etnx0m\nW7UVrY0no8doDPrxdJZXUFA5OTndSE0V6N49CRubP8QqNjZR7HF5jxBFzAjo6B2eaucXtoSxY31R\nKBKwsK7AY0Q2KSctkUgqoKIC/9mzGejlRfhPP5F0+DDx1tb0MzMjLy2N6J9/pqysTL8Acmjnzhoh\nGBKJhJL8fMzchzHAQ8X1uGGU3M5n1tyZdLXtyuaVxzmVm4VtT4HxwQ/yePAEBg92q+HH0/nsUlNL\nuXQpj337slCrh1JRMYi9e/+DtXURvr6P4e3tKva4vIeIIRYdmOqxTzLZODSaw8jlWSbvW4mKOoed\nnS3mQjmHV6yiyLorPpZqfbjD2SNHKCwsRHvsGP1LSzmcl8e4999noJcXO0JD9Qnk1UMw/ObOpUc/\nn7uSt8eMGc6+nft4/5UIsm5X4tDdjA+/nGRQjHUhI5s3p6BU9uX333egUvWnU6fBDBumIDi4B127\ndq3hB2wM90ua033f7eh+pHqHJ4C4uDSmT/cz+fpVAQFVlXWLcnMZv/xvvLH5C1xDQvT+Qd+xY+nr\n5oamuJhEOzv6dO2KRqNh6dxXmfrsszwWEgIqVZVfTqViwrJl/On55ykvL9fHqHl7uzJmzHC2bg7j\nvde/pEunUuaOeJIunUp57/Uv2bo57K5x6RZeios1wAlKS7WYm+cgkx1BqRRwc+vbZAGDapkQsbEt\nvXT3JeJ0sgNTPc6qpGQbCkUFiYk3GDzYrb2Hdk+oL66t9gLIj998z8Vdpfzw7+/o18fxrgTyuLhk\ng1Py6qurMpmMHp0deOPdJ5n4xEQiIo7f5YfMySnk1VcH8N0/w3nhBW+cnR3x8/MiOjq2wSljbYur\nreMC7xdMxhIzVOrZFDh69BwyWRzdu9sREDCW7duTjLYkcmt+RyMCAnD19ubT5SuZ9/hSon/JIqDL\nI+z49BjL5v+Dc0N8WXzkCJULFvPvrw/WWX1XKpUikUjIziwiXltVFhv+EL3aJccDAkaQnXQVz1MH\nGDXEkfnzgxk82I3584MZM2a4wbFqtVo+/vgb4q9coXDHDj4JCUGr1d5XLfraEpMRsebWue/oLFr0\nNC+//DgVFeDg4IajY3+jLYncku+orrZtr7/3FqMDvahUV/VLqFRLWLZiAV9t3YxUKmXx60t4f80S\nkuKSAAyWlD5/8hyPeSbzwZfPErjAha/rEL3//t//MdhhOAmrVrHQzo7kNWt4e+xYft68ud6xr1+7\nka8/z2Xbh//CvbAQeVgY/3j2WS6fPftHXKAgUNmB4gKNCaOfTjY3ydhYaG7qTkeiNb6jBIWC6PXr\nsXd2rlFCPOXaNdQxMZQo+3DaIpKyckugqv4YVF2nm0mJFCiSuW33IxVaO/3ruuncQ3I14x4azuGf\nfqKnmRlz5z7MqfOCPk7sqaeG895fljNmciApJc9wOG0Li1y761Oe6qrzpktkr6gYiKX1bDbtXMUu\nSQYB1pb0jItjy5tv4vDgg0w1orjAjojRi1hz69wbE7rlek/PAWzatJ2cnML2HlKTaMl3pBMabUYG\n7llZnP7iC6J798aib19Uqak4qNVYVkoY9VA+Qx4cRrKyktRrKTX2vZWQyacBFmR2j+Nc/nlOHc0z\nWKZHUKnwnzmTk8eiKSy00j80Fsx/jfMXOhF1Jh4rq3fYn5lDj9zLPGRzg9eompJqtVo2rVzJi++8\now/0nTNnmr7UeA87O7LzO+MhS2WcRSUxFRVMfO01xk6Zot9erF3WPIxexEzBUmkI3YqXQpHA+fNd\nGD3atp3pgBZNAAAgAElEQVRH1DRa8h0JgkBedjaWMTFM6tSJb2NiKDcz45FJk7Dx9SU1NJSPJz/K\nKkUqDz79NM95e+uPq+sremLFCrwHeJNZWsr0VxYywNOT4xERegd79TSvG4mJ3Nq7n+AXnmflR5+R\nGK8F9UD6yYeTWnILZUkaZnIzPl39HJ4uNvo6b7VLA0HNbAYr88+psOpJZ5klZ3vbo83O1m8j0jJM\n4grqLJXlyycSEuJqcoGFutZwW7emUFLiydatKXc59+vyGXUUmvsdDfHzw2vCBCgtJRGgtBSviRMZ\n6u+vF6CN2YUcVfQmITGthjBKJBJuJCRwKuocWzMz9T6npKtX9SENulVOx8cfJ97MjFvbt/PX4YOI\n27KF0lQXRgzvjKXMBolEioA5nTtvxMraHqlUwsQnn0TetStvjx1L8po1LLSzI/qdd3irmp8sOTmL\nVasGcvbC33jm8VtIg6bz7pkzPPLFF+QkJ9f43kx1caqtMXpLDEy/eYNuOrZ27UXi4tSYm+fzxhs1\np2MdvaRyc78jXaR9F3d34t3dSYo6h11eHhKJhBNHz3DOzBMbW2/MBqj43/8yOHOxqlu6BSquhIdz\n/tAxblYGkJlbzFVJBQf++lcmDxumD2nIqhbSYG1tzbY33+WVo8e5oXbg4X7T+C0+lJxSC7SSXdjQ\njaWL3enpbM+1a7eAmo1NCgB5djY+f/4z0+bMAWDZstn6z/Kv7z7X/398cDBarZZPli7lgbIyEl1c\nUGHRoTMzOiomIWKmztmzl/n++8NcutSJzp0TuHSpkC+/3MP8+eP0N6spxxo59++P7QcfUC6Yk3zN\nk+F3jK3nF81juCKBNWtOcP5CMeMndGbGjCpx3xoaxn/2xCHNsMddGEJE1HFkA+RMCPZHKMzR+8B8\nn3mGI7/+qi9X7YQKuZ8jiVECUqmUYqWKaU/m8l3YDlYt/4zEq6msWP2ufmxSqZTUxEQyU1L42soK\nV0Eg57ff2Hb7dr3fQUx0NNs/+4xuhw6R8oA/P728lkLBBcdBLxIWdu2uhY/7Jaq/OYgi1s40Jrnb\n19eHy5eTKSsrok+fQG7c+BF/f3t9mEBr1CDryDeJ1qY7W7Zfory8J7FxafQfOIzYFVvp29eC6OgM\njh0tRpnvz7GjZ+jUaR/z5o1h1pyZlBUWseW9Q8QrS8ijlBmP+TMlOIikDRuItLZGU1jIP99+m2FJ\nSRwZMoSS27eRdOpEz0ceQXXuGokVu/Hq7cOzcx7B3Nyc91e+Y3B8ljIZ/d57D/vUVLrfvElyTg6T\n6mjGC1W9Nc+EhuKUn4+jYEPB1UuUaBSUuC3EwcmJ0tK4uxY+Orql3Z6YhE+sOsbmVzAUO1X7M0gk\nEkaO9Mbc3A7Yj7m5HX5+XnqHeWvEGt2L1Jfm+n906VenT5/k+nUbTp8+yYwZ/phpyzl98AhCSSlD\n7AcilJRy+uARrl+5wi8//MDWVavIUQqUSA8i11oQtX496xcvRjNuHMX9+7P18GG6HzjAw2VlnHn1\nVfauWEHnigqu7o5kpF8u/9g4nyWrHmmwKsezr73Gg+PHoy0pIc/DgwfcqjIq6voOps2Zw4SQEDJL\n1ISXPcz1Cgn9+venh+MAYD9FRZX6/S+dOcMnzz1HytatBDo6khoWxtYVK4iJjm7axTdhTM4SUygS\n+Hq9Amdn+w6dnlNf7JS1tfVdvhGdY9zDoz8fLnuL7Oye6B7ILalB1tLUl7osOEOvV7cmmuL/Wb58\nHZs2naWoyBNz85mcOfNvJkz4gIULhzN/8eN89EEqWEVSqpLy5itTmDV3JgC7t0fS5bcDBGnLuCVY\nk2xty6AhQ7imUJD5+++4mpkxwtwc1e3b5Go0aAcN4smhQ+laVqbvIfBAIx8GTfkOzp27wn/+e57M\nHFe0tmPIVqsp1loT5JnB35bPJjY2kZycQo5HRGDn5IRjeTn5GRlI+vbtUNV+OwomI2I6UchKl9Ap\n05Gv1p/FwSm6wwa9Goqd8vGxJTz8Mmq1413CpnOMH9q5E6eDe9E+9ghQlZ7SkmbCLS2JXdc0p/rr\npaWleqHsrjXjo5fXUiC4ILEOYuvWFMzN6w98ff/9EHJzP+HH/2jp2q0HSqUZfn7deP/9pbz52irQ\nFiNYx8BtC+IT0vVhC88tfIbI87+huK0mhTEMkp/DOzCQgCee4PDOnYS+8hZDNRp2SyR069QJ+x49\n9NVjtVot//jHRt5558W7wiAMxYRV/w4GenlxIjKyztr6vr4+eLr8j7Kxk/Ac+hJXL2kZ636LN1a8\noU9O37ttG0lr1tDN15fRdnb8dvo0nyUl0atHD/phWiFELcVkRMxcKMciK5rcmC54dg7gasw5+poV\nYS50zFUeQ7FT3t6D8PbGYFDoz5s360vMLLSz47s1a3h73Tr85s5l+rx5LR5HU0ti12XBVQ9C1b2e\nKZdj6+6OcO4cfr17MlqdyYECR+RpOWRlKHmt1kqrVqtl5cpNegGRyWT0sLWgslyFWv0+KpUtSUlS\nQkI2cPFsDnZ9HKi07UlPaTdOnzrLihVVK5Qnfz/FeZexdLHtQqVkLBfyzEj//ihJqbeQIEFR4IdV\nbzm+ci0Jw4fj0qcPE5cvJzE2lp2/HGDLFmt8fI4QHDy+xmc3FBNWHZ2AJzg7k3iz8C5/p0QiYcrc\nuSSvr+pm3s2hP5PnTKiRRWBfUcEsf39+uXiRHWVlaJ2dmbFqFYAY1V8LkxGxIX5+pGfkcfTYKQqV\n21BKHBkSOLFDr9AZKpzXs2cXg0Gh1ZfyJdbWXMkq5Yk17zNu6tRGn6+uqV9zpqN1WXADPD1JjI2t\n8frIWbMQBIHEI0fYWFDApeiz5Gu785B9KucvWbPny9No5s/Wf1e7dx9hyxYJPj5H6N+7C1fCw7kY\neZFFwzRky7txMUOFtnIQlZW9wayYQR5DuHEjk0fHjMfauiszZvTHy2sgw4f/g4zps1F7PYJTn6dI\nv6Glq+U5fgs9Qpb5QAZPfQ+7gTc5kfw7bo5qQlat0qcKqdWu2Nm9wJo137Fu3dvMnetHZ4rrfZDU\nFvbvvtjI1rP2lBZk8+TMaQ1+97Wvq1QigdJSHEaNwkGjQSKRiE59Axi9iFW/MfPzlYztdQGLnJs4\n2zlzO39Ihza7DcVORUWd0/u+3lq2Qu/7kkqlSKRStCUlvCuYcz5nKD5pWU2K+K5r6tec6WhdFpyu\nKkTt13Wt0CZ6evKvlV/gfOwww3r0wkWej7X/RHx8fWsJyELWrPkOuTyeSeN7ETLOlfHW1nybmoml\nWxBOfZ7hxo3tODk506NHDrGxWm7d2Ihd76H68UXt2YPDpUucHRSAmex7kmMymDWinLn+vfjogpzS\nrCwyzDP567tz9ZZg9VQhGxspZWXw5psTmDJlLID+QSK1sbkrd1InQJGf/JOtlwsos3bH0/9Ffv31\nVy5f28LkyT76KXNdcXPVr+vO8nK6uLsz/JVXqqx20QIziNGLmO7GDC8s5OZv+3Dq14vZIx9ga1YW\nN6N2E+Pdp0NbY7XR/bjDf93H2V8LGDf+j76HJ347ScIDMxAEJ/oOHMOeXb9w8/bWBv1+LXHe1xcC\nUpcFZ+j16kIZPO0xDl2K5os4J4IHXGOknwcSiaQOAZmI5wB7kjZsYL+1NVqlkhKJFthPbm45Cxf2\noUsXWx7wvkxh5CWGPD6Un//zE/85XdVg99HOPZBkh6Ist+DBUX4UZpch7dGZ0jLo98ANzp3NRKvV\n6j9b9VQha+svKSmR618H9A+S6rXK9O/dESBrVREulUmcynfCLD+f0qsJ+M8a1+jKI7rrF9jEJsf3\nK0Zbnrr6jTlOJuNQZSWX8/NxKC5m1oABbFAomLR8Oa7VcumMAV0fR6GgM71vDSKj13UkXYuZPncc\ns+bOZN/OCP7vg0NUuszBzSOL+fP7N9hRuzEdvHXb1RYshSKB9esTCAlxbZUo8ujoGL79ajcy+QB6\nOj5DbuZPaNVJzP/z4xw4cBpv74G8804S9vZqsrPlrFo1EGdbmb5z+dbNP2NOBTPmzSI2NpFTR09j\nlZWAfUUF6pQ85P16kGlmRpm1NZLt21loY8O3SiX9330XlwEDSNqwgZhyObalBcgeGU3ov67wxAIv\n3nr3Zf0Y163bSv/+DkyZMpbdu4+QnJyNn0c3RgUGEvb55zj078/YKVM4sns3OcnJzFy2DKj6Te7/\n6iv6yuVIcov5x1Ebukpyce7uQcCYUqx69zCpIOSmInYAr0VtnwwVFQwJDKRi715Cs7OxUCi4kZhY\no2yLMeDlOYCL9vuIiVPTw8KChMzbDB2kRa4tI+zDDynLyEOqtiX27DYyk2SM8i1pUFwa67yv3pCk\ntLS01Uoc6WLCAAIDR2H1etWqrFQqxcbWjnnz/IiPT2HLFgkPPXSMVaseqSYgOQQvm6k/1rPzZ+j/\nr+tilKBQ8NvaDUTHOeNnfpOAN5aQGh9PYi2LqTgvD9eQEPKiL7Hlqz3I42IZ6/IkMRGHmX7oRZ6c\nO5Z+jt1ZunSW/toEB48nQaEgYf16El1cmH1HsAC9U18n/hMnPoTVnU5MMpU5Tz2cjWvXrnhZlRN5\nu5QxIYvE0Ig2wGhFzNCNef3cOSRyOZ5duzIsIIB927ez9eJFo3r6DfX359zDFzl94AwJ5vspU1kz\naPRDTJ83j1/++yub/7eV/lZSRvabhNr5GrsiruEyMKZBYanPeW8oZs3MLBN3d1vOV6urZah8js4n\n+dDEiezff9LgtDM2NpHvvotHIrGmT5+kGquy0dGJREXtwMJiOHZ2C7l+/TvWrTtAfn4B8+ZNb/B6\nnT17mR83HybrojWyzqM4eHEPl786QE+LTEatWlXDYhp/R4AGeHrS1bZrVcNf6R8NfwcOcCZpwwa9\nz7Cx0/A/xD8Jizu/SU0/R3rFZiGUlpLety/e2jJADI1oC4xWxODuG9MuJ4cuPXuSGhpKDwcHHDt1\nol896R8dEYlEwq30HCb4ZjNi5AOcO32RzPRcpFIpQwYP4kkPCw5mdMVSpaZzZwdee8O3UXW56nPe\nG673NRJBEIiKSsTaOpLCQi3njx2961w6n2SkSuDXcEmN4NXo6Bg2btzH9esVaLUPAf1ZvHgPXbpk\n8Nxzj/L000Fs3JhPTs5Qdu3CoCPdENUXc3x9fUi6coUDpY449RnLzdQsOlkk8PfPVup9VbXDIKoW\nHyA58xbZttuRlsiJ2hqGuVufGmLlHRSE3/TpdcbQGRL/nPRTBE6axLSn/8SuO5UsJs6bJxY8bEOM\nWsQM3ZgJCgWVRUVEWFlxS6Ggbx0Bhx2N6haNbY/uTPrnJ7h6e+NW7cdflJtL3+efo/tesJDlknFT\nDbT86W4oZk0QBI7sPcCSJRPw9nYlYlcEl789QOJo/xpWSkFqJjlJGq4rjlMmG8+6tSdx7ls17fT1\n9cHS0pK1a09x/TqABHf33rz++mS8vAayZVMo58874uNTSUlJqUFHuiFqr7IO9vfn8KmqmKsiZTnF\npYO4ejWp3mn26ZOXsezjQw/HLjw2RkJeRjqCqhCsrUmLi+PhO/7UxNjYOqfhhsQ/5I1n9D7K4Pnz\n9ecTQyPaDqMWMUPorDNBELh19SpXz541Cr9YdYvmSuIAxkks74oLGhEQQFTUOZYutcXLa1KLmrLW\njhmrHbd0MfoSbkmXSL3SnZift+KgVvPmcLe7rJSUzZt51FLDh9Hg7OeNuXWlPk5LIqmqe69UalCp\nkoE4lEo7FIrrbPp6L1knbmD34GR27drGAw8UsWzZFDIz80hOztGPsfpCQ13TuzIHV8aN60pcXAHd\nunVCpXqYsLCrBn140dExbNp0kLg4S5A+QXZ2Mtt+PomjnZIumhvkOHanopo/tb5p+P1QkNMYMDkR\nk9vYcGb7dhzUauZ5eXE4IYGtK1Z0WL9YYyya6jdha9VOq23N6I6rG4+jWs04R0cOXb5MXkEBysJC\nxru61phSJcbGoiku5rjUipISLYIQSXGxnf4cERHHsbCwYPLk7vj5VU3Bdvx3DzHhifhYSJBb9SBb\nocC8soCnFgzB339oDQGoWhm9rs+DNRRg6zdzJllpaYwKDCQ2NvGOVdQbpTLWoA9PZx3+858XuXat\n6lyurv3w6J3M2f1JPOcxAN+AAKK2byfrjj9V9yAxZE2Jnb7bH5MTsZbmAt5rdOOtz6JpTRpyVhta\n9fWaOJGKvXuJFIQaUyqdlWKRXcaSSTnY21phaedMbm6R3tkdEuJAQMAIvUXlvWKZfjWxWGWGVtiL\nUt4Xt2qhMDpfU3q6nMzM0axffxGnO3mw3axld5WTVoeHk9SnDxKJRYNW0R/WoYBKFQ5YUFJiRk8H\nB/zz8oiJjcWmc2f6OTrSvxH+1KY+VDpyySNjpcFw74iICDw8PHBzc2P16tUGtwkJCcHNzY2hQ4dy\n4cIFAK5du8awYcP0/2xtbVm/fn3rjt4ArVWa5l6hG6+muJjjglDNotHo328OdZWrbkyvwwvHj6Mu\nLNRfv5Lbt3FbupSJy5fX6MSt6/04dqwvQcFBjAgIoLS0lKiouBotz159dQPffHOZ2NhE/efNyy1k\nuH8uo50ysbLScu5cXI2xZ2VlERNTROfOrsTEFJGVlY0gCAbLSetK1ISt+ZzAsWUGS2BXvx45OYUE\nBVkSGhrIs9MKUMZsIm5TKPO7dSPv2DF2/forX+48fFcCd2NKgDe0jdjtu/Wp1xLTaDS8+uqrHDx4\nECcnJ/z8/Jg6dSqenp76bcLDw0lISCA+Pp7Tp0+zePFiTp06hbu7u17QtFotTk5OTJs2ra5TtSot\nKU3THtRn0TSmaKIh6koxaihmLDE2FtWePfRYuJBxU6fqr5+rt3dV0GxaAYGBj9Z53urO7oKCbBQK\nBT4+3fHxmU1Y2BHk8jP0d1Dh9tzTXI4rJUsjxcXmQRISkvSJ235+Q8jIyOX48WQsLPZTViYQGPhH\nyk5ExHE6CUUE/vnP3PjhB70Yz3vzz/oA3tpWUfXrERAwQm9BeXm9zhBvV0688w7/zcnHVipl6Ojx\nnEv14pt3P2TO39/UuyEaU5iwrm3Ebt9tR70idubMGVxdXenXrx8AM2fOZOfOnTVEbNeuXcydOxeA\nkSNHUlBQQFZWFg4ODvptDh48yMCBA3FxcWmDj3A3LSlN0x7oxutqYKgKRUKT6q435mYxJPLV93tu\n+HAOR0cTdscnpEt5qR4MW9dYqju7e/W6QEGBAxUVFUilUsrKasaaDVIk3PFh9UepTNS/J5FIyM8v\nYdAgGZ6eAlevysjPVyKRSLhyJZ7Vn5xnpv1pKp5/usEA3oauhy4n1UEQONnNgdQULWlpzvTqGkCp\n/RB+Cb/OFUUCpF6v95o2dZre0d0cxkS9Ipaenl5DeJydnTl9+nSD29y8ebOGiG3bto3Zs2fT2uhM\nd4CHJ01q0ynjvfZlNLbhbO1xNeZmMSTygiDUuV9Tm99Wd3bv2nWYb7+9Tp8+d/up6lvZ69/fmeXL\nvfUO86NHz/PyS58Sdy6P9PQgzivTCX/vB1zdLJi7YAE2NjYGLe7GXI+s5GRcV68moH9/tv15KUcy\n1DhbVlChljBjxkg8PQfcVZnDf9Ys+nt48PHH3/DOOy82eJ7mljwSaZh6fWKNvcC15//V96uoqGD3\n7t089dRTzRhe/STGxhL/3Xekf/99i30MWq2WtR+vQ61W8+KcN9BoNHed6176MnQlmVWqKr+MSiUw\nY4b/XUnEtcfVXJ9gffs1diw6AgJG4O3tikQiwda2K6tWTeCDDyYwfHiRvvGvIAgc3rOfJUsGGvRh\nVT+Gt7crw4e7g5kZuereVFZmcynTEmQyHn5qBkP8/HD19q6RJK37PnV+reqfSyf8ut/t7GXLGB8c\nTMnt2/i/9VcGTZxGN/8C8nNLgVqVOapdmz17otiyRcLu3Ucadd11FnBt36JIy6jXEnNyciItLU3/\nd1paGs7OzvVuc/PmTZycnPR/79u3jxEjRmBnZ0ddLF++XP//gIAAAqpZCoaIiY5m38aNVFy/zkNa\nLf2BPYsXo/Hw4LEXX2zQx6CrzLnwb3/j1IEDjAoMJHJ3JIe2ZJOU+g4Ruzqz6asf+POS+e3my2go\nBqm+camVymb5BOvyJbYkHqquxr+JsbG4JV3CUjLGoA+rNjo/WfjeYxSVVlAhKJlkp2HI4EH6cVT3\nH+q+Ty+fSOxtrWp8rtjoaMxPnjRYkqiuOLzq1+bTlet599n1WHby05cMWrfubcaM7MRz9Vx3Y3Nz\ntJSoqCiioqLa/Dz1VrGorKzE3d2dQ4cO0bt3b/z9/QkLC7vLsb9hwwbCw8M5deoUy5Yt49SpU/r3\nZ86cSVBQkN5vdtcAmpHZLggC8VeucGrtWnpfv84A4Jy7Ow+8/nqjqlYc2rmTxLffxvrll0n6eQ8n\nlN3JyLQmI3sMGuEKVhJ7ZDbX6NyliJlzhzPrqccbrADRFkRFncPOzrZGDNKYMcP116AxlSnuxVjq\no/pUVCYbR0baj9xO/p3x/c1Y7GLPYY2mRu/H+o7zj3/8j4vnzbGzm0JRyUG06qO8+2Ywcxe/AMCV\nK/E8O3sTPc3TcLQcyEDpOBK1hymXZzF97jgGe7vWqHzS2HPXRqvVsnPn4TslgxaiVH7Lu+/254kn\nAvQPRXGaeDftUsXCzMyMDRs2EBgYiEajYcGCBXh6evLNN98AsGjRIiZPnkx4eDiurq7Y2Njw/fff\n6/dXKpUcPHiQTZs2teqgdbE+GqWSZJWKOMBOqdRbDHWhK/Hc8/ZtJpSWcuDTT7GRSulv50xRpRue\nPeFSthkOtg+glBYS8vooXnx5Dslxce3iy6gvBqmlPpam+vh0Y9FoNKxe/R3ff/9Ro85TOzXHtrsj\n06Y+idnRg0QmZiFz7KSvCBsRcZzAwFFoNBqeDZzBlsjtmJmZ6Y8THJxEt25q+vQZyvGj4WRfHYBd\nr153hPIy50+lIjMLQN79AhevJGLdKxW1tCq5O3BK1eesqxqt7tyNuRZ11RxLjosT26q1Aw0GuwYF\nBREUFFTjtUWLFtX4e8OGDQb3tbGxIbeN5v2FOTl0nzwZrztP0NjoaIPTpupTDF2J5+SVK3G2t0eS\nlsaYMWMQVAJXb3enVPIblRoX8glDo+oDgEwm67AhGy0ZV3P7GH71VRi7dvXkq6/CWLLkuQa3rz0V\nLS4WqpKvb2SxM8mdYNU1+gNXryaxY0cRqan/Iy9FwY2jFqx5/2P+tnK5/jj+/oP5MfS/nDryApXl\nfRjfewJr/r4DjfltcpVaNAXu2PZ5lJQUBaWlKm5oTyArrXJjVC9cWFv4r15NanLn7eTkLFatGsiU\nKWNZv3Yj27/+gU6+LmL4RDtgtEURG0vton66qWS+jQ3yhARkEyeSmVaA6+zpXLl6iwq1Eh8vdzCT\nkXDtJp9teK/NxnYvqCv/0L6iAsuUFMr79SPb3LzBG+6DD/5FaOglNJrBSCR/obhwObbdbjB37lBW\nrHit3jFUn4r+9FM4WzYfwKXPSHo6PoMi5jNSUxLxHvowYVsUSLUqesrimCQfSopZHLnSdMY+NZJ/\nbfqSqKhz9OjRmdT4eL54/zes1B6cy07Gz+siqRmeZOaq6NHtCdLzjzBjehlf/d9n7N+zn5vJGby4\nbAEA56Ki9AUWf/nvr+yPiKWn04PIZOPQaA4jl2c1uWbavZ7aGytiUcQmUldYgLQwhlGrVtGlc2eu\nKRTINBoenzCBotxcFixZ1PCBjYzasV26UIBTa9diFxdHjrk5D73xRoPxSu+/H0K3bltYtSoHQVuO\npkTN40/b8f77IQ2Oofq0+JlnJjN48CB9UcT+rsN48OHurFr1O2rNQ1iap6Cp7EmCRkmxtoxn/vIE\nb374bo3jHI48RFJWLhrZZUrKH+XyLXOK8pSUlFdA0XoGeQ7nwTEjkclkBAXXnEVUd65Pe/pPDPIe\nbLC7VFMQwyfaF5MVMcM1svzx8vqjaqffuHHtPMq2oy4Rd+srIS/6GOYXLzKwc2dSL15k31dfMWbe\nvHotMZlMRkz0GQpy7bGRvolU4kDugV/5+2OPNaltXO3p5ZYf15OTpULLg0jworwij0w0lFacxdOs\nakqv84vpsJDJ+ct7o0m84czJ45ZcvlxOJ2l/bDtnINXa0r17PiqVqsljaUkVio7qcrgfMFkRq/4D\ntbTcR1TUVebM6XPfPB3rEnFPzwHsEspRl5bi2qcPF2/cwNzXV587WZ/D39LWmRfm5jL0wkkqCmXs\nKIKXli3Td/tpLNWDYf/0JxfmP/8WlxWlyLBFzS06Wak59PsmDuz4mdSrSXft/+fXFqJQJHBxfQI+\nQ29yPcEClfYEcvkoutgOxcLsNFlZaqKjG65421AVisYugNxv4RMdicb3+zJCdD/QESPk5OaWs3Pn\nb+09pBbRmARkHdVFXKuN4LoiCUEQkEqlDPb3p7O5OfuBzubmePv718idrCuo96uv3mb61EeQlJUh\ndbZmVs+qkstNaRsHNQNZhw3zwsPdCbAGyRdI6Ib7oApKSir428rl/PuXHwweQ/fdPvPMIN57rz/e\n3m44OJShFWTY2vVmxoyRjeouVDuotnbYSGJsLEU7dvC/jRvb1Hcr0nxM1hIDSElJYcWKqj6GLi5v\nExHxHYcOVTVCbUz99o5GU1cUdTe6uVCO09UfuXS2B4MHu9WY+iQoFBzbuxelUoli374Gg3qzkpMZ\nWKt2fUuJv36Lp56EH8LWEzz5Bc6eH4idnW29++j8Y97eVYLdt29v9u6tQCZL4ObNSuAPR3JzEuir\nBxP7y2Rs//e/2RAT0+C0W+TeY9Krk3UFJU6ZMrbJ1kN7Urs9XWODNBuzn66Tz8AlSwAa1datLXJI\nawfFVlYeIjP1OItC/oS//9AG968rGLe5LecEQWDvtm0kffEFT5iZkVhZSZaPD5JevfB5/PEOI2TG\nVJ9MXJ1sBg01QjUWmlsBob797kpb2raN2Px85BkZaIYMqXOFrbnxZQ1R24eXlZGLX/l1ks+fxM+v\n4WsY/AUAACAASURBVE7utQODo6NjWLFia7NbzkkkEtx9fFA6ORGZnk5/lQqhooIHn3qq1SpPNNdK\nrE5bfR/GhEmLGNQMStT1MTQ2qi/hN6UBSn1L/4YEro+vL153kqlrr7C1dQ6p7rOcPXGBrmaHUBdb\n0cdVSsbGjXweE8O4+fMblRO7cuUm3nnnxXpWpxsfPlGUm0vP556jy969ZMpklN28WWOsLaUxpY3q\nQqxP9gcmL2LLlv1RAig4eHw7jqRlNLcBSkNJ3dUFztvfv8568veiHtal6Es80+siMv+hXAr7mbg4\nM4JtLDhfWcnl8HCAem/Q3buPsGWLBB+fIwQHj29x+MSIgICq4NilS5nUiqETTS1tZAixPtkfmLyI\n6TAm34EhmtsApb6l/6bENrVlQGf15iSzvV35X3w8clUxWjMzhAo5qNWMrKfe/ebNPxMaWrWAU72q\nxMiR9oSETG1RE4+2CJ1oDStRDLD9g/tGxIzdd9AWT96m3qBtFdBZ+7NVFBczYvFiul++TLRMRmUD\n07g5c6Zha2t7ZwGnZgNenf+zJZ2hGkNDD8na77dGkK0YYFuFyYuYqfgOdD98dWEhn2dkMMjGRv96\nYxAEgWP79nH63HX+8veQZi1utFVAZ22roqetLapOnRi0dGmjbtC2WsBpivXe0EOy9vut0epNDLCt\nwuRFzJR8B4U5OZg//jiD9u7F/PHHa9zYDd1wibGxnP783+yJd8V7SORdOYX3EkOrcoasCl1zEl1f\nyfpo6QKOoevXGOu9oYdkfe/rrMO2thJNHZMXMVPxHcRER3Ptt99wUKuZ1KsXh6OjiZPLMbO2Zoif\nX71ddr77x2dcu5DLAK0TLnizask2/vnhZuYv+ROz582655/F0KpcXVZFY90ALV3AqX6e0tLSRlvv\nDT0kTekh2lExeRED0/Ad1HUzKJVKtq5YUW+XncUf/p3d733MtWgZEhl06tKDV96bwKSpVdaNISuk\nLRZCmrIqd6/cAIbOkymXY+vujnDuXIPC09BD0lQeoh2Z+0LETMF3UNfNMMTPD2tra4NPep0Q2Ts7\nQ3k5seoudBPCUZb2RSL5w2dkyNppyAJqTqBmU1blWmLBNEWADZ1n5KxZCIJA4pEjjRKehh6SpvAQ\n7cjcFyJmKhi6Gep70icoFBT/8gu5Dz5I+QP+/H3xCFwHOhO5K5KbyRkGrZBd2dmoBQFfB4d6LaDm\nBGo2ZVWuJRZMU1ai6zpPUW5uo4WnoYekKTxEOzImnTt5v1C9WmlibCznjh5Fk5XVYK6loYqkfefO\nRRAEbvzwA4/Z2HCgpIR8NzdmvPgiUqn0rhzHplZDbUrDkdqfqyg3t0Zbtto0N8e0qecRaR5tda+L\nImaCNKVcsi4BXOLoiJCZiWtIVaVW3WsZsbFYSySM+OAD/WqhQt+5+zGUyoPMm9dP37m7PRHLRHds\nxATwJmDs0fktpSlTMUNTVEEQ0IwbR0FcHDZdu2KRlMS+jz6iu4cHPo8/jsy6W6tUQ21tRCf6/YlJ\nipixR+e3Bo11Jtflr9FZNafWrkW4dQtlWhpes2bh4+vLb7+db3GgZlshOtHvP0xqOtlcn4jI3cRE\nR/Pzp5/CsWNMMTMjSa3mopMTtsOHE/TSS+L1FGky4nSyEYiBha2Hj68viUFB3Far6ZadjVVpKT3t\n7ZmydOl9a92KdEyMqzpgA9TwiQgClaJPpNlIJBJ8Ro6kk7k5+1QqMmUyunXqpL/GIjURBIGIiOON\nsjQa0yuhKce73zEpEYM/fCITly/HNSTEKHwiHeEHa2gMhTk5WAYFERgaSkBoKD1q5Ws297imiC5u\nLjY2scFt62vG0pzj3e+YlE/MWGluHXhjGENH+GxtSVPi5hrjszXUayA19SQhIVMb1WugIyPGiZkg\ndTXHeGlJMPn5yhbVXm/uGJoavHqvj9vRaErcXGPi2GofLzX1PyiV5rz//nAGD3Zrj4/YarTVvW5y\n00ljwtfXh+nT/VCpqurlZ2Xk4q+KJ+/WrXs2lag9BpVKYMYM/0b1bGyP43Y0qqdSCUIkhYVqjh07\nX++29flsddskJCQQFbWCixdL6dHjMbZtu8GKFVuJjo65J5/LmDCp1UljQ/eDTUnMQJP3N/6/vXOP\ni6rO//9rhosIBmriDMygKAPKTVBBqrUNNEXxK2by/S7bRWut9eG3Qtvq+7XaLdhKcak2i/qu2Sr0\nS8mtNLwglJdx0wJGRVHQBAQEFBxQ7veZz+8PmmkGh2Eu58ycmfk8Hw8eDxjOnPOeM+e8z/vzvt5q\n7MQYT+Df755Gj9NCfPDeTxBPNa73uqkyMJ28ytZ+uYhmg8MDB47js89aMH++7ppSQ/LY5PI2vPXW\nQhBC8P77F9DT029SC2tHgS4nrYxUehaTJnnClfTi5HuZcL7ThyM3xBBGp8DV/QKefnoa6yU9xtQz\nMrlfe6qs0OUaqKn5EfPm+SMl5UmTPp/KnygU8tDYSGzer2jXy0lHiWDpCq3Hxs5FWFgg+Hw+xK4E\nPw+2o6NzEDXXs9He/uska6aOr+s8x8bORWioBDweD6GhEkYUmCH7NSRKZysMXz7fuHETzc0TcOrU\nmLvcAoakWAC/WnipqYuRkiLhVGUEl+DEctKc+Xu2AiEEX336KcadOYNdtS14+o+PaSknVevpCTu+\nxP0rPXG+dCKWLXNl9MJl+jybOvzVXuYeaKLpy7p+/Vs0NLhAJFqA7u4QrF9/CDNnKvDssw8jOnqW\nwWVxwwcCU3QzqiWWn5+PmTNnIjAwEFu3btW5TUpKCgIDAxEREYGSkhL1662trUhKSkJwcDBCQkJQ\nWFio8/05ObUQCuOxZ08NnnjiLRQXXzDx43CTUpkM2557Du3btwNtvdi3vQp/e+5PKJXJAAAXiovx\n3Wef4cq3R9Db6YGf/90NdEzHodwbkEqvmOXMJYRgx4f/QFrabvV5zsmp1ekkNtRCUGFqLlN4VBSi\nV60C6esbWmL09WFeUhLCo6KM2g/XUPmydu16ATExAejrU0Ch4CEgwBcpKcvhouzBlieeQM2ePYgX\nClGbk4M9aWnq64BiGnotMYVCgeeffx5Hjx6FSCRCdHQ0EhMTERwcrN4mLy8PlZWVqKioQFFREdav\nX69WVhs2bEBCQgK+/vprDA4OoqurS+dx+voIPDyGTPD+/iC4u7sz+BGtS6lMhot5eei+0YLqW144\nIeehk/jiXON49OZdRR/GYLy7O4S9vfBwdsZDIRPwVxkgjg6Fq/sgkpKmmeXMrSovh/jSOYgSVuLE\naaK3o6rKQqgUi1FV3zaihWXu8Fd77TahspzKyirh6uqBMWPk6OvLQ1eXN3g8Hjw8PCDs7cWdGzfA\nmzqVlsUxhF4lVlxcDIlEAn9/fwBAcnIycnNztZTYgQMHsGbNGgBATEwMWltb0dTUBDc3N/zwww/I\nzs4eOpCzM7y8vHQep6KiAjXXDoLnFIL77kvAl1+eh4uLjPWcIrYcy5rLLFU9p0t1NSI87yCjSQgv\ngSfu8fJA+Aw3XD28D8LBQcz39sbJoiJk8ceipWUeBKQAHR3eAEy7uTWXbEt8fLDzyBGcOzMZgbOv\ngzj5au13+PJu50efYs+ZyTh38nu8uvnNu47PxPBXe+42IZe3ISFhIqKjh5TTvn8dQk7GKcSJvdTf\n87vXrsHn3nvhD9tX3tZGrxJraGiAn5+f+m+xWIyioqJRt6mvr4eTkxO8vb3x9NNP48KFC5g7dy62\nbdum08r641NTUfv/fsBRBc+i4WS2WvYM9z3xeDz0tbTg4pSpUN5wA3/qFdxpFiMoLApjwgJRm52N\n9p4ejImIQOSDCVgwcRwE48fCzVtssk9seDF8a2s3Xn5jLuKXx+Py5Wta+1Vt+13633GqgY+feyfB\nV7IKpfk7sLF9G554aoHWw4SJ9Al7btkcGzt3mD9rozrJVfU9r3jpJQCwK+VtLfQqMWMGsw5/3+Dg\nIM6dO4fMzExER0dj48aNSE9Px1//+te73v/l++kY192Pi/U3MW7sQUyPTALgz9oTii3H8kjLLOee\nn+Hh54fx7QSbFnei3vUq7njeQdGpMXjowUgMtrejeepUeDc2QhITYvJNPdzRrrlkCx3vBIlkCvh8\n/l1OYtW2092d0DSxAzWn+iDqr8TUsZ4IHryGq4dbMQZ9Wufm1q1WzJnTjmefTbpLKVK00fwu5FOm\nYPDMGQBAYFiYlSVjF6lUCqlUyvpx9CoxkUiEuro69d91dXUQi8V6t6mvr4dIJAIhBGKxGNG/XPhJ\nSUlIT0/XeZzVMTFwlvciQNyECc5KCJ+YzOpNwVbLnpGWWcHByagqL/+l3ESgVW5y7uRJ9bKqsqwM\npw4fNrmd8nAL0JAlm0rxTRrTj8ANGyBWuuDcy8cwWH8AHZgK0t+PmP/8z7vOzeTJ4/Hll564fPna\niJEze8oDMxfVd0EIwc3Ll3H5zBm7V2KxsbGI1bC409LSWDmO3uhkVFTUkL+qpgb9/f3Yu3cvEhMT\ntbZJTEzE559/DgAoLCzE+PHjIRAIIBQK4efnh6tXrwIAjh49itARLIymqiqcqS6F1/VKzJ80Cdf2\nfo76EwdYi9qw1bJneAmKKs+Lz+ePeLy5sbGQhIaq3yusqjI6b6qo6DwWLHgBu3dXa0UflR4T1fuW\nhIbqHH6hUnxjJw8tqVtaOvDcf/vjj7/hYc68Ztxu1j43Mlkp0tL2jBrpBOwrD8xcXDw8UPz117j+\n5Zd4KiQEHpWVNDLJEHotMWdnZ2RmZiI+Ph4KhQJr165FcHAwtm/fDgBYt24dEhISkJeXB4lEAg8P\nD+zatUv9/o8++giPP/44+vv7ERAQoPU/Teb99a8ghKDo/fdx5eZNuNXUQJiczGrInS3HsmYJimbr\nZn3HM3d5e/NmC2pqAuDp+TMiIxcZ5FPUF2G8B2PhtWEDluiQ1RCnviXzwEzNVbM0tGEne3Cm7Ojw\n3r04k5EBn7FjEdfVhfNRURgQiWw6AdJQTJnSQwjByy+n49y5VgwMBKKn5xE0NGTBza0a/v5T8PHH\nq/QmSJoztWi0chilUomvP/0UEyorsWjcOFanDtlSqx9dk6XsLaihD7tvTy2YPBnz//xnKE+dwr3d\n3ei7fRv3bdjgEE8qU/Kmysur0N4egUceacEXX/RBoejHpEnOePvtRxEQ4DeqT9GcCONI1qaKa5cv\ngxw/jnpCUBASwkoemCm5atb20dlzWok14YwSi4qLG3pSHTkCmVCISQoFAMfJoTH0Ate8eX19l+KH\nHz5GTU0r/Px2QaHwBo8HhIcHGXTM0ZTRSIxUDqO5jJw/dix+HDMGzXw+xi9YwPgNa0qumrWnYNlz\nWok14YwSAxznSaXLIjD0Ah9+87a0dOCNN/zx3//9exw6JEV1tdxgOZiuzdP0+4imTcOEri74JyWx\nsow0xpK0x1pNyq9wSok5ypPKHItg+M0bFDQbCxZI4OTkhBUrFrIksXGyWaqcyFBLkjrV7RtOKTF7\nhymLwJCb11r+H0ta04ZakvZaq0kZgjPRSUfAlCikqagiYY4WARuJs1IpvLy9tZSrrrw5CnvYdXSS\nEOIQT0VLWATU/6M7d8xRXBWOCCc6u2577jmHyVxmey4mm726jO03Zi3ozEbHghOWmGhwEBfz8gDA\n7q0Fti0CNq09a6cojIa5fc4otgknLDEMDCBmmLVg6ac+F6wMpmRg2torlcmwJy0NtTk5d3Uk5cJ5\nU+EoY+Io2nBCifX+cpNpWguWLh7mQrEyUzJoFpWPVPhtDPqWqFw4bypGKsB3BH+rI8OJ6GTFpUvq\naJEho96ZxNLH46oMozG87k+xYAHarlzhnMxsjZ+jmA9b0UlOKDFNEYxJQ2Cig4El0x50Hbug4Ecs\nXny/Rr8xy8pgKMNTFNrkcnhOmmSyzNauY6RYHrueO6mJMb2+TIlCDffhsNVbzBBU8l++fE0tQ75S\niUtlZZxLOxm+RJ0bG2vWeePSMpRi23AiOjmc0bK+zYlC6YqwWbpmU5f88oZCxC9ZgvDQQLRduYJy\nmYzznT9NOW+W7jVGrT37h3PLSUMwpRcWl/xOuuS/f24Heq5cgHBwkDH5uHgTj7R8nx4cjO+++4nR\n5oa0aoFb2HXGvrGY0guLS0XAuuQPCgtTTz5iSj4u5nWNlMd2+fI1xqaT06oFx8ImlRhgfC8srhUB\n65JfNGkMI/Jx/SbWXIbu/9e3ePfdbzFJdB9jCapcemBR2MdmlZgpvbC41K9Ml/xnpVJG5LP0TWzs\nslWzamHlfz2CoNAwdX+03l4l/PycMXeu6f5Arj2wKOxis0rMFLheBKySjxCCpro6PBAfb9J+zLmJ\nTfGjMdkf7fr1VtTWdmP+/JFHwRkClx5YFHaxSce+vcOEQ9rU1jPGHJupYIlUehaNjY24cqUVN28S\n9PYGw8+vFS4uTbTu0Y5wmGRXR8aaEVRTjs1korA505dU7+daJJaijcMku9oahBDk559m5Mths40O\nG8dmMlHY3LpHmjzruDiUT4xpCCH49NOvcPbsPYykBljTIW3qsZn0PZkyfYnrkVgmsZVBwZbG5paT\nXFk2yGSl2LXrOAoLxyEm5iHce+91Rnw41myjfFYqheekSWiqq4PAzw8dLS2cb+FszdpXS2NLg4J1\nQZNdf4ELCZxDZUMXMTjoizFj5qKsrAYKxY9ISZljdu8qa0ZQ58bGorKsDJ3ffgthSgrnFRjwqwU5\n0NaGd250YI6HUv26vUCbPerHZpQYl5YNqtmP7713Hv39BDzeNYhEYoSGBjFy81jD2hzt/CqVSvx9\n84d48bUU8PnccqW2yeXAspWoPQzMWQa7S6cwZVCwI8Gtq1EP1nR6D0f19G9u7sW8eTKEhYnxxBNT\n0NLSwcj+reGkHu38FhwswLHdt1BwsMBiMhmCTFaKQydv4LTMFT4+S3Ba5oqD0gbIZKXWFo0xaLNH\n/diMJca1LGy5vA3p6fcz2nzPmtbmSOc3J/tLfJN9HG4DAkR5P4wvMo7jsw/2YdWaBXjsqd+zKpMh\nOIqVYkrQw1GwGSUGWCcLe6SlnSllT6Nh7Zo/Xec3efXv4OXlhezNP4LvwcdAD/DMK48ifrlp1QRM\nY0ozAC5grMuAjevNXrApJWYNp7clAwnWtjZHOr98Pg+dnUpUuOejs5P/y2vc8UTYopXChQCVvWBz\nKRaGYm5OjbWy51UpFtODg/H1jh0ImDFDS7lYg+0f/BNTpvkifnk8Cg4WoL76Bp7duNaqMtkqXOpr\nZ2mslrGfn5+PmTNnIjAwEFu3btW5TUpKCgIDAxEREYGSkhL16/7+/pg1axZmz56NefPmMSe1AZg7\nQFXT0Q0A5VeuIHrVKtYDCao20NcuX4bnuXPw8vZm9XiGsG7jWixdsRR8Ph9LVyylCswMuBSgshf0\nLicVCgWef/55HD16FCKRCNHR0UhMTERwcLB6m7y8PFRWVqKiogJFRUVYv349CgsLAQxpXqlUiokT\nJ7L7KTRgKqdGc2mX3dmJMWVluF5VxXrLaGOc+zSD2/awtsvAHtGrxIqLiyGRSODv7w8ASE5ORm5u\nrpYSO3DgANasWQMAiImJQWtrK5qamiAQCADA4sXdTEarzv7wA245OyN4/HjMjo3Fka+/xp7z51k1\n/Y1x7qusTSZKniiWg7YJYha9SqyhoQF+fn7qv8ViMYqKikbdpqGhAQKBADweDw8//DCcnJywbt06\nPPvsswyLfzdMRqv+a906dUnLvQIBhOPGwT8pidVooSFPaprBrRtbsUy53tfO1tCrxAy9EEaytk6d\nOgVfX1/I5XIsWrQIM2fOxIMPPmi8lEbCVLTKWqb/aE9qW86NYlPRUMvUMdGrxEQiEerq6tR/19XV\nQSwW692mvr4eIpEIAODr6wsA8Pb2xsqVK1FcXKxTiaWmpqp/j42NRayZ0Tgmc2qsYfqP9qS21dwo\ngB1FQy1TbiKVSiGVStk/ENHDwMAAmT59OqmuriZ9fX0kIiKClJeXa21z+PBhsnTpUkIIIT/99BOJ\niYkhhBDS1dVF2tvbCSGEdHZ2kgceeIAUFBTcdYxRRGANpVJJjhw5RZRKpVWOby4nTpwhly5VEKVS\nSS5dqiBS6Vn1/7j42YqLL5DU1N3k9dePkjfeUJLXXz9KUlN3k+LiC2bvW6lUkosXr5KXX/6evPkm\nIS+//L363FgLa30HSqWSnDpyhFPfvQq27nW9KRbOzs7IzMxEfHw8QkJC8Lvf/Q7BwcHYvn07tm/f\nDgBISEjA9OnTIZFIsG7dOnzyyScAgMbGRjz44IOIjIxETEwM/uM//gOLFy9mWyePCvmliWFZWaVZ\nKRhsQYZNKCcjNF2MjZ2L0FAJeDweQkMlWiVP5qaXsEFUVDhWrYpGX9/QZPO+PoKkpHlmd/0AuFlb\naK3vwBGbQ9ptsutI7N17GBkZZzBnTgh8fJKgUBznVC93zR73ASEh6qaLGzYEjrr80lxWOTkt4Nxn\nU/XDEgp5aGwkjPbFkkrPwtvbi9FaVlOw1ndgC0m0tMe+magurv7+ybh9ewpKSn6Am9s8+PldxKZN\n89S93ImVmi4Ovwj/r+4W9l1sQ0P7fYhdsNSgpovEzD71bMMVRcMm1voOiA00h6RNEc1EM6LH4/HR\n00Mwe/YlKBQTAPy69LBWTZtmftiZ1h4UXunCJP+H0FW/EGVltQY1XeS6w98Ripit9R04chKtwygx\nzYurt/d7BAW54bnnon7pC9Zu9aaLmhdhi48AD9y5g+8HeOjvh1FNF22xGNrekMvb8MILAaira4Kf\nn8Bi34GjJtE6zHIS0L+c4YI5rtlfP/9APrZsvYJ7J/bB23cWli0bg/HjJ9jd8stesfV++GxAfWIW\nQOVU5wmFII2NZg2vNRep9Cx62m4Bh/cDy1bCfbyAKjAbgOvBFWtC504ayPAUBWNQmeOLU1MhSUmx\nmjleKpPhxslDcJWdxhIfH7jKTqNBehClMplV5DEVpVKJ7W+/DaVSaW1RLAabqSQU3didEjMnT0bV\nBofH40ESGmrStJ+R8rqMwR7atRBC8HFaGni7d+PEwYPWFsdicDFnzd6xGyVWKpNhT1oaanNyEC8U\nojYnB3vS0gy2Xgyx4AxRUEwkOTI5WdsafJOVhQ0xMej5xz/wuLs7qjMysCkuDt9kZVlbNIugCq6k\npi5GSoqEBldYxm6ik+b2pzcktUJf3R/T9Xu2GmkqlcnQU12NqX5+8KyvR39zM1oGBhD44otYuXq1\ntcUzCUNzB4mO4nbq1Gcfm1Ziwy8aU/JkDEmtMERBMd1ZwlbbtageJgcuXAAZGMDHTk4Y7+YG/8BA\nTvXlNwZDcwdpFw3rYNNKrKysEv/3YRnE4skICws0yXoxxIIzREFxPdHUFEypXlCdh9amJpC4OMzx\n8sLtWbMgr6lhV1gWMDR3kHbRsC42qcRUF01TAw/jGoX45MMzEIhkSEgIUz8pDbVeDLHgDFVQlk40\nZbtEytTqhTa5HGs++0zrYWJKkMTaGOqisOX+bvaATSoxV9KLMU0yNJfeg3EDEWguvYypzh1wJaaZ\n8IZYcIYoKEPKaphUPGyVSJlSvaD5uWx1KTwcQ10U9miF2xI2qcRmRUej4UYLCo5J4dK+B/J+CWbF\nLza5PMiQm86Yuj9dDl4VTCgetkukTAmS2OscRUNdFLTcy3rYpBK7eOYMjmTtwW95pXhc4ond8koc\n2XUDIt9JJt3ESqUSOzZvxrOvvcaI81mXg5dJxcP2pHBjgiTWrjllG0OtSkcobucqNqHEhi/BwqOi\nsHDFAgxMcIbEzw8DP12EeNoktMnlIISobzZDldOJgweHkjLDw7FwxQqT5dTn4I1iUPFYomOBoRYI\n2wqVQhkNm1Biw5cqPB4PYfPmobKwEDvk7fixahriukpQWV8Pn6lT1U/M0ZTTN1lZkGVnQzIwgGe8\nvbEzIwObPvgA0WvWYNVTTxktpz4HL9OKh+08MkMtEEduAUPhBpwuANfXrfL82UsoOFaBpspeKG8G\nIUBxHH1jWzB2Uhf4nm4Y7OhAzNix+AOfj51KJSpdXO5STkqlEsdzc1GzeTOe8fDAZ11dmPbnPyNu\n+XKTl5X6updqdqmwRtROn6/OHKz9uSjasPU9m4tDNkXUt1QJj4rC7N9UYNv7Z9BG3IFbHvjtxFq4\nzgpGzJ/+hOuVlajdsgV8Dw+gpweLXnkFccuXa+2fz+eDx+dD2dmJj93d4dLZqX7dVPQ5eE2N2hly\nURoS9WQrGdNeopG2juoa8BIHOFTSLaeV2GhLFT6fj85OoPP2l7jRPw6i/n5M7eqCk5MT+E5OBimn\npupqBKSnI275cpw4eBDy6mqzZGbDwWuI8tEXHaTJmMzBVSsHAPb/61vsy9gLzzkrIPRJdpjvmdPL\nSUD/UkUqPYuLp47jzncH4Zv0OMbecy/u4XVCPG0arpSUQDBtmpZySt640VIfixEM6U1lyIAIrvfe\ntyW42OxQdQ1M7u+H4HYX3irxgsItDvf41eGVTfdz5nt2yOUkMPJSZajnVh5CBgawIC4Ox2/VoelO\nI/x/uXk1fTLmRBytiSGZ4IZEB2kypvlw2ZrVvAb4PB66e4Bps+vQp5gEwP6/Z84rsZFQfXE1WVmo\na2yEUii0ydC+Pl+WIcrH0OigpZIxubzcMgculxZpXgOHe12wPKgKC577T/Tz3Bwi6dZmlZjqi7t9\n/ToOVFXBu78f02B7T53RMt1VQyda6yox3k93bypD0i0slYxpr50cuG7Nqq6BeAeMEHPeJzYSpTIZ\njnz8MURjxuAxgQB7mprQ0NeHpc89ZxOZ4sYMO9UcqMvV6J8j9JZ3hLmZbEIHhQyDC9OJzMEQ+Q11\n2ltj2K+uz0ODBxR90EEhw7D1Fs6GyG9Ir31zZgowCe0tT7EWNusTA2y3hbOK0eRXKYaBtja8c6MD\nczyU6te5WHhNOzlQrIHNLicdhbNSKW619WD/YWDlMkAw3h1zHnrI5pfTFMeDLicdEJmsFIdOBET0\nGQAAF6FJREFU3sBpmSt8fJbgtMwVB6UNkMlKbX45TTEec2aq2jNUiXEYzUGsAHDlSjlWrYpWD2Id\nPuy3TS5n5CJn42ahN+DdGHtOuOL/5BpUiXEYTWd5U1MWZEUDqKy8rra2hg/79fL2ZuQiZ+NmoTfg\n3Rh6TsydqWrvUJ8Yx/nHP/bi4sVbUA5MwvUfXOEsuYk5URONrp80BKb2w/Y+bR1jz4kx/k+upNzo\nwmo+sfz8fMycOROBgYHYunWrzm1SUlIQGBiIiIgIlJSUaP1PoVBg9uzZWD6sDQ7FMO6fMw1+5Bpa\nSq8gavJEjJHL4dYkgyvpVW9jSCqGITC1H7b3aesYe06M8X86osWrN8VCoVDg+eefx9GjRyESiRAd\nHY3ExEQEBwert8nLy0NlZSUqKipQVFSE9evXo7CwUP3/bdu2ISQkBB0dHex9CjtGNRTlh1OFuN71\nJbp4wruGojDVXZWNLq208+vdmHJORkvH4WLKjaXQq8SKi4shkUjg7+8PAEhOTkZubq6WEjtw4ADW\nrFkDAIiJiUFrayuampogEAhQX1+PvLw8vP7663j//ffZ+xR2DI/Hw507XbhfeA6Tmxtwa5IIt+/M\nuuuCZypnjo3cO1vP52MDY8/JaI0nHXnWgV4l1tDQAD8/P/XfYrEYRUVFo27T0NAAgUCAF198ERkZ\nGWhvp0mPplIqk+Hmvw9jhr8PVt03G980NqL25CGUhk7ResIy1V2VjS6ttPPr3TB9TjQTo7fduIEg\nDw/16/aOXp+YoSdguLOOEIJDhw5h8uTJmD17NnXcm0F4VBRWvPQSxnt6ov7aNYz39MQjL788ov+E\npjI4Lm1yOVyXLUOQhwdcly1zGItXryUmEolQV1en/ruurg5isVjvNvX19RCJRPjmm29w4MAB5OXl\nobe3F+3t7Vi9ejU+//zzu46Tmpqq/j02NhaxGk8pR0dfyyFdkSh7HWJL0U+pTIafT56EYGAAS3x8\ncFwmwxUXFzi7u1vNJyaVSiGVSlk/jt4Ui8HBQcyYMQPHjh2Dr68v5s2bh5ycnLsc+5mZmcjLy0Nh\nYSE2btyo5dgHgJMnT+Ldd9/FwYMH7xaApljoRV/LIXd3d3WLnu7ubprK4MDYQhmaVdpTOzs7IzMz\nE/Hx8VAoFFi7di2Cg4Oxfft2AMC6deuQkJCAvLw8SCQSeHh4YNeuXSN+AIrxhEdFYezLLw+1Hubz\nIfT0xISoKFzMy4NQIxLV6OICrxkzQM6edTjHLsWxo8A02dUGUDVF5AmFII2NCHjhBQC466lLCEHV\nRx+pt+NyE0UK83B9/qfDDgqh6A7He06adNdTt725GQEvvICmujoI/PwcxrFLGcJRo8DUErNRRnrq\n2kIra4pjQttTU/RiTD0el+vrKPYL7SdG0Ysx9XiOWF9HsV+oT8xOMCQ65cj1dRT7hSoxO0IVAJge\nHIysHTnwlsu1/u/I9XUU+4UuJ+0IVZPEy5evoejcvXDz1q6u4HJLa1ouxRyOdi6pJWZHaA6wFQrj\nkZNzHC4uxVoNFK3VUUKpVOLvmz/Ei6+lgM+/+9lJy6WYw9HOJY1O2hFMDbAlhCA//zTGkXbweDw8\nEB+P7777CfHxD2jtx5go55HcI/ho0w94If1BLF2xVP067fzKHFw/lzQ6STGI06dL0Nraj4qKHWht\n7cP27f/SeeHoW3KUl1che9d1nPzg/9CwaxcKDhZg//4OlJdXaW1nSJRzT1YOVsU9iy8yTiPK+2F8\nkXEaq+KexZ6sHAC08yuTOOq5pMtJjsBE7lZ5eRUOHerD4sX9uHhRjEmTLmP37jFYuPAEVqxYqLWt\nriWHTFaKXTvycaX4JpS3w3GnzxOHeuUgJz+Fl+hpfPDeTxBPLUbgVB5Qe9WgKGfy6t/By8sL2Zt/\nBN+Dj4Ee4JlXHkX88ngAjl3zxzSOei6pEuMIhvoxCCEoKPhRa2mn8oXV1gLt7RH48MNcyOWFGDNm\nLNzcnsAbb3yPDz74HmvWRGNO6JQR0yyioqLg5uaGbe+fQUuJO1yaXLBS2ILDzqHwv+8+jHG/gKSk\naQgOno6q8nKDopx8Ph98Pg+dnUpUuOejs5Ovfl0F7fzKHI54Lh1eiY1kAWkqC0IIdmzejGdfe02n\nU9ocjM3dKi+vwv79HfDzq0JoqATA0HzKsWPHIiurBm5uInR3C+HpKUJ9/W1MnuwHwAsbN87G8uVx\n4PF4I6ZZ8Hg88Pl8dHYBiv6DaBl0xtmBQfQM8AFSgI4ObwBDCsiYJ/716ht4If1BxC+PR8HBAtRX\n39D6v6PW/LGBI55Lh1RimgpqJAtIU1k0VpaBt3s3ToSHY+GKFYzKYmju1miRRx6Ph44OBfj879DS\n0oe+vjJ0dweAkB3o6poA4FfrR58Cksvb8OhSwHXJb+A60QenfwrGrHuUeGC+H9y8xWhuHtremCf+\nuo1r1b9rOvUpFCbgRHRSqVRatJavrKwSqX/5EVHCEkRN9NCK5GBqECpqCQYGBPj54kTUn98Gb1KH\n1+/txEUPF1S6uCB6zRqseuopxuQZ3mpHV/H2aJFHqfQsGhsbceTIRXR3u6K6uh8TJ3bCx8cDYWGT\n4eLigY0bkwFwv2ULxT6x61Y8lspr0bRmgsOeQE2jC8pO5MMz5B6Q8WMx7/e/x/Tg4KHoXHYtQiIj\nwePNQ2hFKeZ53oOL3d1Y9MoriGN4hqYhVo3mNHB39wK0txP16wAQGzsXhBCEhQUhO7sWoaEjp1g4\n4pKDYr9wQompxrOzXcun8h0NWTN8DAzw8BtfBVp8JgJNTQB+9fcMKYvv0N3DA+nuxie4By6dnept\nmMRQpSKXtyElRYKQkACUl1epl3YqRlN0FIo9wgklRvr6LFLLN/wmv93cjSnPPI74xCVaFpCmsvjb\npjfh8uB6rH/tBZw4eBDy6mpWZDOE2Ni56t9VTv3hjKboKBR7gxM+sSN//KPFWipLpWfh7e2ldZM/\n9NAc1o5HoVCGsOumiBWXLlEnM4Vi59i1EqO1kxSK/UNrJykUCkUHVIlRKBSbhioxO8XRGuNRHBeq\nxOwUOgyE4ihQx74elEola4XfbMH1xngUx4U69q3AiYMHhwq/Dx60tigG46iN8SiOC1ViOvgmKwub\n4uJQnZGBZ7y9UZ2RgU1xcfgmK8vaoo2KqiphoK0N2yoqMNDWpn6dQrFHqBLTwcrVq/FwSgrQ1ze0\njOzrw6KNG7Fy9Wpri2YQbXI5XJctQ5CHB1yXLXOIxngUx4UTtZNcg8/ng8fnQ9nZiY/d3Vkr/GaD\nUpkMP588CcHAAJb4+OC4TIYrLi5wdnenPjGKXUKV2Ag0VVcjID0dccuXW73w2xjogFyKo0GV2Ag8\ntnGj+nemu7myharV9mSx2OGGRVAcF6rEbBRdswFUuWHN992HkBdeQFNdHQR+ftQnRrFrRnXy5Ofn\nY+bMmQgMDMTWrVt1bpOSkoLAwEBERESgpKQEANDb24uYmBhERkYiJCQEr776KrOSOziayaylMhn2\npKWpm0uOq6rC4U8+Qe2OHeDxeLQrCMWu0ZvsqlAoMGPGDBw9ehQikQjR0dHIyclBcHCwepu8vDxk\nZmYiLy8PRUVF2LBhAwoLCwEA3d3dcHd3x+DgIObPn493330X8+fP1xaAw8muXERXMmujiwu8ZsyA\n29mzmNzaitzycswMDcUqoRAnlEqa7ErhBFZJdi0uLoZEIoG/vz9cXFyQnJyM3NxcrW0OHDiANWvW\nAABiYmLQ2tqKpl9aPbu7uwMA+vv7oVAoMHHiRMY/gKMxPJlV2dsLF6EQQaGhGGxvxw0fH4h9fICB\nAfD5fJrsSrF79CqxhoYG+Pn5qf8Wi8VoaGgYdZv6+noAQ5ZcZGQkBAIB4uLiEGLHETJCCPLzT7Nu\nVWpNeSYEt69fR9fRo7h85gwkKSmIT0vDlCefRG9zMwoIwSB17FPsHL1KzNALf/iNq3qfk5MTzp8/\nj/r6evz73/+GVCo1TUobQDWnsry8yux9KZVKbH/7bSiVSp3/b5PLoYiLQzOPh4EJEzDL3R3jqqpQ\n/PXXuHjmDMZ7eeH+9HQsTk2FJCWFOvYpdo3e6KRIJEJdXZ3677q6OojFYr3b1NfXQyQSaW3j5eWF\nZcuW4cyZM4jVmOyjIjU1Vf17bGyszm24ymhDbU1BXbM5wrDeubGxIISgsqwMtdnZmOvjg6NdXVqT\nvFXQkWwUayGVSi1iuOh17A8ODmLGjBk4duwYfH19MW/ePL2O/cLCQmzcuBGFhYVobm6Gs7Mzxo8f\nj56eHsTHx+PNN9/EwoULtQWwccf+aENtjeGbrCzIsrMhGRjAH/h87FQq9Q7rNWToLkU/ulJVKOxg\nleG5zs7OyMzMRHx8PBQKBdauXYvg4GBs374dALBu3TokJCQgLy8PEokEHh4e2LVrFwDg5s2bWLNm\nDZRKJZRKJZ588sm7FJg9wOSsx5WrV8PLyws1mzeD7+EB9PToHdZryNBdin4sNbiZwh60nxgDMDkG\n7lhuLqo2bcLA5MlwuXULAenpNlMxYEvQvmuWxyqWGMUwDBlqq2K05Yut1mzaGrTG1H6gSszCjLZ8\nscWaTVtEK1WF1pjaNJxYTqam7jYrmmcLGLN8oc5my3BWKoWXtzcd3Gwh7Hp47qVLFSZF82wJzZSI\nhz08cLSrC/5PPXVXSgTwa9SRRhsp9oTd+8TsWYEBhi1fNK21eKEQx3NyUEydzRQrYSsrAk4osebm\ndmuLYBFGS4mgzmYKl7CV9BNOLCdtPcWCSWgCK8XasJV+YvfLScoQNIGVYm1sbUVAlRjHmKtRN0ot\nMIo1sLX0E6rEKBTKXdjSioD6xCgUO4IQgoKCHxEf/wDnLCerdHalGI6lmiJSKCNBCMEXn2Zh3772\nEfva2eN1SpUYQzDVFHG0hogUii5kslJsfH4bDv/jIsY4zUFOTi3S0vZAJivV2o7J5p1cgS4nzUSz\nKaKT0wIoFMfh4tJkchmVqosF7V5BMZRSmQwXDx9Gz40WXLl0L245/xYtgz9gbcosrPzdI+DxeIxf\np6ZAl5McJSoqHKtWRaOvj4DH46GvjyApaR6iosKN2s83WVnYFBeH6owMPOPtjeqMDGyKi8M3WVns\nCO6AEEJwOj/fph+augiPisK8pCRgYAAdfc4gg4dxjzgQQWFhar8YU9cpF6HRSTNhqimisQ0RKcZj\nKxnoxqK6Blua2zBnnhumKLqBZRPQ0tJx1zZMNO/kGlSJMYBc3oaUFIlWU0Rj4fP54PH5UHZ24mN3\nd7h0dqpfp5iHI9SktsnlWJX+52EdObQbczJxnXIR6hPjEHs++ACCadO0GiIma/QXo5iGMR1EKOxB\ny44cANoQkR1sLQOdYhxUiVEcAlvKQKcYB11OUigUi0BTLCgUB8deU0TMhSoxO4de+PaDOkWkvNza\nonAKqsTsHGMvfE2lRxUgNyiVybAnLQ21OTmIFwpRm5ODPWlpKJXJrC0aJ6COfTvF1NwozYRQAHaZ\nHGpr2FqTQktDlZidYuyFr6n0IhUK7F2/Hu08HhKDglBrh8mhtgRNEdEPVWJ2irEXvqbSE4hE8FEo\nMBHALF9fHOvupk9+K0NTREaGKjE7xpgLX1PpfefuDtLVhUFC8B1An/wcgLYtHxmaJ0ZRozkR+0BW\nFng8HpavWUOnY1MYwa4ngFMlRqHYPzTZlUKhUHRAlRiFQrFpqBKjUCg2jUFKLD8/HzNnzkRgYCC2\nbt2qc5uUlBQEBgYiIiICJSUlAIC6ujrExcUhNDQUYWFh+PDDD5mTnEKhUACAjMLg4CAJCAgg1dXV\npL+/n0RERJDy8nKtbQ4fPkyWLl1KCCGksLCQxMTEEEIIuXnzJikpKSGEENLR0UGCgoLueq8BInCC\nEydOWFsEg6ByMguVkznYutdHtcSKi4shkUjg7+8PFxcXJCcnIzc3V2ubAwcOYM2aNQCAmJgYtLa2\noqmpCUKhEJGRkQCAcePGITg4GDdu3GBcEVsCqVTKyH4Iy/WITMnJNlRO49F37XBJTkszqhJraGiA\n3y91dAAgFovR0NAw6jb19fVa29TU1KCkpAQxMTHmymzT0E4EFFOh145uRs3YNzRLe/jTQfN9nZ2d\nSEpKwrZt2zBu3DgjRbQPHGFYBYUd6LUzCqOtN3/66ScSHx+v/nvz5s0kPT1da5t169aRnJwc9d8z\nZswgjY2NhBBC+vv7yeLFi8nf//53nfsPCAggAOgP/aE/dv4TEBBgrLvLIEa1xKKiolBRUYGamhr4\n+vpi7969yMnJ0domMTERmZmZSE5ORmFhIcaPHw+BQABCCNauXYuQkBBsHGFqT2Vl5WgiUCgUyoiM\nqsScnZ2RmZmJ+Ph4KBQKrF27FsHBwdi+fTsAYN26dUhISEBeXh4kEgk8PDywa9cuAMDp06fxxRdf\nYNasWZg9ezYAYMuWLViyZAmLH4lCoTgSVq+dpFAoFHOwSMb+7du3sWjRIgQFBWHx4sVobW3Vud0f\n/vAHCAQChIeHa72empoKsViM2bNnY/bs2cjPz+eknCree+898Pl83L59m5Ny/uUvf0FERAQiIyOx\ncOFC1NXVcVLOV155BcHBwYiIiMCjjz6KtrY2zsn41VdfITQ0FE5OTjh37hzj8jElp6Hvt5ScIyXQ\nX7hwAffffz9mzZqFxMREdHR0jHpMiyix9PR0LFq0CFevXsXChQuRnp6uc7unn35ap4Li8Xj405/+\nhJKSEpSUlLC2HDVXTmCoSuH777/H1KlTWZGRCTn/53/+BxcuXMD58+fxyCOPIC0tjZNyLl68GGVl\nZbhw4QKCgoKwZcsWzskYHh6O/fv347e//S3jsjEpp6Hvt4ScCoUCzz//PPLz81FeXo6cnBxcvnwZ\nAPDMM8/gb3/7G0pLS7Fy5UpkZGSMflBWwgXD0IxW3rx5k8yYMWPEbaurq0lYWJjWa6mpqeTdd99l\nVUZCzJeTEEKSkpLIhQsXiL+/P2lpaeGsnCo2b95M/vd//5dxGQlhVs59+/aRxx9/nLMyxsbGkrNn\nzzIunwpz5TTm/WzL+eOPP2plPGzZsoVs2bKFEEKIl5eX+vXr16+TkJCQUY9pEUusqakJAoEAACAQ\nCNDU1GT0Pj766CNERERg7dq1rJnC5sqZm5sLsViMWbNmsSGeGibO5+uvv44pU6YgOzsbmzZtYlpE\nAMzIqWLnzp1ISEhgSjQ1TMrIJubKaanPachx9CXQh4aGqiuCvvrqK4NcHYy1p160aBEaGxvvev2d\nd97R+lvVBtkY1q9fjzfeeAPAkD/npZdewj//+U9Oydnd3Y3Nmzfj+++/V79GzIiZsHk+Vft55513\nkJ6ejhdffFEdUeaanKp9ubq64rHHHuOsjExgKTnNfb+5cuo79s6dO5GSkoK33noLiYmJcHV1HVUe\nxpSY5s07HIFAgMbGRgiFQty8eROTJ082at+a2z/zzDNYvnw55+SsqqpCTU0NIiIiAAD19fWYO3cu\niouLjf68bMo5nMcee8wsC4dtObOyspCXl4djx45xVkamYFNOJj+nuXKKRCItC6uurg5isRgAMGPG\nDBQUFAAArl69isOHD48qj0WWk4mJicjOzgYAZGdn45FHHjHq/Tdv3lT/vn///hGjguZijpzh4eFo\nampCdXU1qqurIRaLce7cOVZuCnPPZ0VFhfr33NxcdQ4f05grZ35+PjIyMpCbmws3Nzc2RDRbRk3M\nsbxHw1w5mfyc5h5HM4G+v78fe/fuRWJiIgBALpcDAJRKJd5++22sX79+9IOa58YzjJaWFrJw4UIS\nGBhIFi1aRO7cuUMIIaShoYEkJCSot0tOTiY+Pj7E1dWViMVisnPnTkIIIU8++SQJDw8ns2bNIitW\nrFA7DrkmpybTpk1jzbFvrpyrVq0iYWFhJCIigjz66KOkqamJk3JKJBIyZcoUEhkZSSIjI8n69es5\nJ+O+ffuIWCwmbm5uRCAQkCVLljAuIxNyjvR+a8mZl5dHgoKCSEBAANm8ebP69W3btpGgoCASFBRE\nXn31VYOOSZNdKRSKTUPbU1MoFJuGKjEKhWLTUCVGoVBsGqrEKBSKTUOVGIVCsWmoEqNQKDYNVWIU\nCsWmoUqMQqHYNP8fS0mL+ZSVs/MAAAAASUVORK5CYII=\n",
       "text": [
        "<matplotlib.figure.Figure at 0x11c9f8050>"
       ]
      }
     ],
     "prompt_number": 158
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "area = good[pd.np.logical_and(-0.14 <= (good.longitude - mean_long) , (good.longitude - mean_long) <= -0.12 )]\n",
      "area = area[pd.np.logical_and(0.02 <= (area.latitude - mean_lat) , (area.latitude - mean_lat) <= 0.035 )]\n",
      "print area"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "       Unnamed: 0_x                   address  \\\n",
        "14833         10194      [u'1551 Dolores St']   \n",
        "15240         10450      [u'3139 Mission St']   \n",
        "15314         10489           [u'78 29th St']   \n",
        "15406         10542      [u'3459 Mission St']   \n",
        "15415         10548      [u'3438 Mission St']   \n",
        "15528         10621     [u'1600 Guerrero St']   \n",
        "15557         10642  [u'3017 Mission Street']   \n",
        "15596         10667       [u'35 Precita Ave']   \n",
        "15648         10703      [u'3100 Mission St']   \n",
        "15831         10829     [u'833 Cortland Ave']   \n",
        "16044         10980     [u'745 Cortland Ave']   \n",
        "16050         10983     [u'400 Cortland Ave']   \n",
        "16066         10990     [u'1600 Guerrero St']   \n",
        "16081         10997      [u'3188 Mission St']   \n",
        "16120         11024        [u'1580 Valencia']   \n",
        "\n",
        "                                              categories           city  \\\n",
        "14833  [[u'Coffee & Tea', u'coffee'], [u'Wine Bars', ...  San Francisco   \n",
        "15240                     [[u'Coffee & Tea', u'coffee']]  San Francisco   \n",
        "15314                     [[u'Coffee & Tea', u'coffee']]  San Francisco   \n",
        "15406  [[u'Coffee & Tea', u'coffee'], [u'Mexican', u'...  San Francisco   \n",
        "15415  [[u'Coffee & Tea', u'coffee'], [u'Breakfast & ...  San Francisco   \n",
        "15528                     [[u'Coffee & Tea', u'coffee']]  San Francisco   \n",
        "15557                     [[u'Coffee & Tea', u'coffee']]  San Francisco   \n",
        "15596  [[u'Grocery', u'grocery'], [u'Coffee & Tea', u...  San Francisco   \n",
        "15648  [[u'Coffee & Tea', u'coffee'], [u'Ice Cream & ...  San Francisco   \n",
        "15831  [[u'Bakeries', u'bakeries'], [u'Coffee & Tea',...  San Francisco   \n",
        "16044                     [[u'Coffee & Tea', u'coffee']]  San Francisco   \n",
        "16050  [[u'Coffee & Tea', u'coffee'], [u'Sandwiches',...  San Francisco   \n",
        "16066  [[u'Coffee & Tea', u'coffee'], [u'Cafes', u'ca...  San Francisco   \n",
        "16081                     [[u'Coffee & Tea', u'coffee']]  San Francisco   \n",
        "16120  [[u'Coffee & Tea', u'coffee'], [u'Food Stands'...  San Francisco   \n",
        "\n",
        "                                         display_address  \\\n",
        "14833  1551 Dolores St, Bernal Heights, San Francisco...   \n",
        "15240  3139 Mission St, Bernal Heights, San Francisco...   \n",
        "15314  78 29th St, Bernal Heights, San Francisco, CA ...   \n",
        "15406  3459 Mission St, Bernal Heights, San Francisco...   \n",
        "15415  3438 Mission St, Bernal Heights, San Francisco...   \n",
        "15528  1600 Guerrero St, Bernal Heights, San Francisc...   \n",
        "15557  3017 Mission Street, Mission, San Francisco, C...   \n",
        "15596  35 Precita Ave, Bernal Heights, San Francisco,...   \n",
        "15648  3100 Mission St, Bernal Heights, San Francisco...   \n",
        "15831  833 Cortland Ave, Bernal Heights, San Francisc...   \n",
        "16044  745 Cortland Ave, Bernal Heights, San Francisc...   \n",
        "16050  400 Cortland Ave, Bernal Heights, San Francisc...   \n",
        "16066  1600 Guerrero St, Bernal Heights, San Francisc...   \n",
        "16081  3188 Mission St, Bernal Heights, San Francisco...   \n",
        "16120  1580 Valencia, Bernal Heights, San Francisco, ...   \n",
        "\n",
        "                                                id  \\\n",
        "14833                        noeteca-san-francisco   \n",
        "15240                     coffeeshop-san-francisco   \n",
        "15314                  cafe-seventy8-san-francisco   \n",
        "15406               nenas-restaurant-san-francisco   \n",
        "15415                  cafe-st-jorge-san-francisco   \n",
        "15528            the-other-alex-cafe-san-francisco   \n",
        "15557                  paisanos-cafe-san-francisco   \n",
        "15596              precita-food-mart-san-francisco   \n",
        "15648           sweet-chinito-coffee-san-francisco   \n",
        "15831                 sandbox-bakery-san-francisco   \n",
        "16044     martha-and-brothers-coffee-san-francisco   \n",
        "16050            progressive-grounds-san-francisco   \n",
        "16066        new-alternatives-cafe-san-francisco-2   \n",
        "16081              tierra-mia-coffee-san-francisco   \n",
        "16120  st-lukes-hospital-coffee-cart-san-francisco   \n",
        "\n",
        "                                               image_url  \\\n",
        "14833  http://s3-media2.ak.yelpcdn.com/bphoto/wxOWWii...   \n",
        "15240  http://s3-media4.ak.yelpcdn.com/bphoto/ABoU3IE...   \n",
        "15314  http://s3-media2.ak.yelpcdn.com/bphoto/Me5vvLj...   \n",
        "15406  http://s3-media4.ak.yelpcdn.com/bphoto/sAMw4ix...   \n",
        "15415  http://s3-media4.ak.yelpcdn.com/bphoto/lQk_hLz...   \n",
        "15528  http://s3-media3.ak.yelpcdn.com/bphoto/K4qD7mp...   \n",
        "15557  http://s3-media2.ak.yelpcdn.com/bphoto/qnz5hry...   \n",
        "15596  http://s3-media4.ak.yelpcdn.com/bphoto/xAfS-IU...   \n",
        "15648  http://s3-media4.ak.yelpcdn.com/bphoto/4Mt626B...   \n",
        "15831  http://s3-media2.ak.yelpcdn.com/bphoto/zFG2M6t...   \n",
        "16044  http://s3-media4.ak.yelpcdn.com/bphoto/pNX602h...   \n",
        "16050  http://s3-media4.ak.yelpcdn.com/bphoto/c5yPKDV...   \n",
        "16066  http://s3-media2.ak.yelpcdn.com/bphoto/xmDyPex...   \n",
        "16081  http://s3-media4.ak.yelpcdn.com/bphoto/iye6oMl...   \n",
        "16120                                                NaN   \n",
        "\n",
        "                                  name       phone  postal_code  rating  \\\n",
        "14833                          Noeteca  4158245524        94110     4.5   \n",
        "15240                       CoffeeShop  4153683802        94110     5.0   \n",
        "15314                    Cafe Seventy8  4159702233        94110     4.0   \n",
        "15406                Nena's Restaurant  4152858985        94110     4.0   \n",
        "15415                   Cafe St. Jorge  4158142028        94110     4.5   \n",
        "15528              The Other Alex Cafe         NaN        94110     4.0   \n",
        "15557                    Paisanos Cafe  4158261865        94110     4.5   \n",
        "15596                Precita Food Mart  4158212113        94110     4.5   \n",
        "15648             Sweet Chinito Coffee  4158213388        94110     4.5   \n",
        "15831                   Sandbox Bakery  4156428580        94110     4.0   \n",
        "16044         Martha & Brothers Coffee  4156427585        94110     4.0   \n",
        "16050              Progressive Grounds  4152826233        94110     4.0   \n",
        "16066            New Alternatives Cafe  4158265282        94110     4.5   \n",
        "16081                Tierra Mia Coffee  4155253142        94110     4.0   \n",
        "16120  St. Luke's Hospital Coffee Cart         NaN        94110     4.0   \n",
        "\n",
        "                                    rating_img_url_large  review_count  \\\n",
        "14833  http://s3-media4.ak.yelpcdn.com/assets/2/www/i...           278   \n",
        "15240  http://s3-media3.ak.yelpcdn.com/assets/2/www/i...            86   \n",
        "15314  http://s3-media2.ak.yelpcdn.com/assets/2/www/i...            95   \n",
        "15406  http://s3-media2.ak.yelpcdn.com/assets/2/www/i...            25   \n",
        "15415  http://s3-media4.ak.yelpcdn.com/assets/2/www/i...           106   \n",
        "15528  http://s3-media2.ak.yelpcdn.com/assets/2/www/i...            14   \n",
        "15557  http://s3-media4.ak.yelpcdn.com/assets/2/www/i...             7   \n",
        "15596  http://s3-media4.ak.yelpcdn.com/assets/2/www/i...             3   \n",
        "15648  http://s3-media4.ak.yelpcdn.com/assets/2/www/i...            47   \n",
        "15831  http://s3-media2.ak.yelpcdn.com/assets/2/www/i...           259   \n",
        "16044  http://s3-media2.ak.yelpcdn.com/assets/2/www/i...           123   \n",
        "16050  http://s3-media2.ak.yelpcdn.com/assets/2/www/i...           222   \n",
        "16066  http://s3-media4.ak.yelpcdn.com/assets/2/www/i...            25   \n",
        "16081  http://s3-media2.ak.yelpcdn.com/assets/2/www/i...            44   \n",
        "16120  http://s3-media2.ak.yelpcdn.com/assets/2/www/i...             1   \n",
        "\n",
        "                                                     url  \\\n",
        "14833      http://www.yelp.com/biz/noeteca-san-francisco   \n",
        "15240   http://www.yelp.com/biz/coffeeshop-san-francisco   \n",
        "15314  http://www.yelp.com/biz/cafe-seventy8-san-fran...   \n",
        "15406  http://www.yelp.com/biz/nenas-restaurant-san-f...   \n",
        "15415  http://www.yelp.com/biz/cafe-st-jorge-san-fran...   \n",
        "15528  http://www.yelp.com/biz/the-other-alex-cafe-sa...   \n",
        "15557  http://www.yelp.com/biz/paisanos-cafe-san-fran...   \n",
        "15596  http://www.yelp.com/biz/precita-food-mart-san-...   \n",
        "15648  http://www.yelp.com/biz/sweet-chinito-coffee-s...   \n",
        "15831  http://www.yelp.com/biz/sandbox-bakery-san-fra...   \n",
        "16044  http://www.yelp.com/biz/martha-and-brothers-co...   \n",
        "16050  http://www.yelp.com/biz/progressive-grounds-sa...   \n",
        "16066  http://www.yelp.com/biz/new-alternatives-cafe-...   \n",
        "16081  http://www.yelp.com/biz/tierra-mia-coffee-san-...   \n",
        "16120  http://www.yelp.com/biz/st-lukes-hospital-coff...   \n",
        "\n",
        "          complete_lookup_address   latitude   longitude latlng_confidence  \\\n",
        "14833      1551 Dolores St, 94110  37.744485 -122.424377             P1AAA   \n",
        "15240      3139 Mission St, 94110  37.747092 -122.418914             P1AAA   \n",
        "15314           78 29th St, 94110  37.743976 -122.422201             L1AAA   \n",
        "15406      3459 Mission St, 94110  37.741485 -122.422537             L1AAA   \n",
        "15415      3438 Mission St, 94110  37.741585 -122.422473             P1AAA   \n",
        "15528     1600 Guerrero St, 94110  37.745592 -122.422286             L1AAA   \n",
        "15557  3017 Mission Street, 94110  37.748917 -122.418127             L1AAA   \n",
        "15596       35 Precita Ave, 94110  37.746720 -122.418262             L1AAA   \n",
        "15648      3100 Mission St, 94110  37.748101 -122.418269             L1AAA   \n",
        "15831     833 Cortland Ave, 94110  37.738920 -122.414090             P1AAA   \n",
        "16044     745 Cortland Ave, 94110  37.738907 -122.414763             P1AAA   \n",
        "16050     400 Cortland Ave, 94110  37.739178 -122.417611             P1AAA   \n",
        "16066     1600 Guerrero St, 94110  37.745592 -122.422286             L1AAA   \n",
        "16081      3188 Mission St, 94110  37.745968 -122.419653             L1AAA   \n",
        "16120        1580 Valencia, 94110  37.746832 -122.420213             L1AAA   \n",
        "\n",
        "       Unnamed: 0_y short_category      \n",
        "14833         14833         coffee ...  \n",
        "15240         15240         coffee ...  \n",
        "15314         15314         coffee ...  \n",
        "15406         15406         coffee ...  \n",
        "15415         15415         coffee ...  \n",
        "15528         15528         coffee ...  \n",
        "15557         15557         coffee ...  \n",
        "15596         15596         coffee ...  \n",
        "15648         15648         coffee ...  \n",
        "15831         15831         coffee ...  \n",
        "16044         16044         coffee ...  \n",
        "16050         16050         coffee ...  \n",
        "16066         16066         coffee ...  \n",
        "16081         16081         coffee ...  \n",
        "16120         16120         coffee ...  \n",
        "\n",
        "[15 rows x 21 columns]\n"
       ]
      }
     ],
     "prompt_number": 159
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "def get_num_and_mean(x):\n",
      "    return [x.rating.mean(), x.shape[0]]\n",
      "\n",
      "\n",
      "rating_by_cat = merged.groupby('short_category').apply(get_num_and_mean)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 184
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "ratings = rating_by_cat.copy()\n",
      "# ratings.sort()\n",
      "print ratings"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "short_category\n",
        "accessories             (3.8757183908, 696)\n",
        "accountants             (4.26470588235, 17)\n",
        "active                  (4.10714285714, 14)\n",
        "acupuncture             (4.5873015873, 126)\n",
        "adult                   (3.76086956522, 46)\n",
        "adultedu                (4.44230769231, 26)\n",
        "adultentertainment      (3.22580645161, 31)\n",
        "advertising             (4.63333333333, 15)\n",
        "afghani                 (3.82954545455, 44)\n",
        "african                           (3.8, 10)\n",
        "amateursportsteams       (3.83333333333, 3)\n",
        "amusementparks                     (5.0, 1)\n",
        "animalshelters                     (4.3, 5)\n",
        "antiques               (4.11344537815, 357)\n",
        "apartments            (2.84350479512, 1147)\n",
        "...\n",
        "watches                     (3.97983870968, 124)\n",
        "waterdelivery                           (1.5, 1)\n",
        "waxing                                (4.38, 25)\n",
        "web_design                   (4.68085106383, 47)\n",
        "wedding_planning             (4.40178571429, 56)\n",
        "weightlosscenters            (3.65789473684, 19)\n",
        "wholesale_stores            (3.98557692308, 104)\n",
        "wigs                         (3.88461538462, 13)\n",
        "windowsinstallation         (4.01225490196, 204)\n",
        "windowwashing                (4.52777777778, 72)\n",
        "windshieldinstallrepair       (4.16666666667, 3)\n",
        "wine_bars                   (3.92068965517, 145)\n",
        "wineries                             (4.28, 100)\n",
        "womenscloth                (3.80873261206, 1294)\n",
        "yoga                          (4.6320754717, 53)\n",
        "Length: 607, dtype: object\n"
       ]
      }
     ],
     "prompt_number": 183
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [],
     "language": "python",
     "metadata": {},
     "outputs": []
    }
   ],
   "metadata": {}
  }
 ]
}