Assignment+3.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---\n",
    "\n",
    "_You are currently looking at **version 1.5** of this notebook. To download notebooks and datafiles, as well as get help on Jupyter notebooks in the Coursera platform, visit the [Jupyter Notebook FAQ](https://www.coursera.org/learn/python-data-analysis/resources/0dhYG) course resource._\n",
    "\n",
    "---"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Assignment 3 - More Pandas\n",
    "This assignment requires more individual learning then the last one did - you are encouraged to check out the [pandas documentation](http://pandas.pydata.org/pandas-docs/stable/) to find functions or methods you might not have used yet, or ask questions on [Stack Overflow](http://stackoverflow.com/) and tag them as pandas and python related. And of course, the discussion forums are open for interaction with your peers and the course staff."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 1 (20%)\n",
    "Load the energy data from the file `Energy Indicators.xls`, which is a list of indicators of [energy supply and renewable electricity production](Energy%20Indicators.xls) from the [United Nations](http://unstats.un.org/unsd/environment/excel_file_tables/2013/Energy%20Indicators.xls) for the year 2013, and should be put into a DataFrame with the variable name of **energy**.\n",
    "\n",
    "Keep in mind that this is an Excel file, and not a comma separated values file. Also, make sure to exclude the footer and header information from the datafile. The first two columns are unneccessary, so you should get rid of them, and you should change the column labels so that the columns are:\n",
    "\n",
    "`['Country', 'Energy Supply', 'Energy Supply per Capita', '% Renewable']`\n",
    "\n",
    "Convert `Energy Supply` to gigajoules (there are 1,000,000 gigajoules in a petajoule). For all countries which have missing data (e.g. data with \"...\") make sure this is reflected as `np.NaN` values.\n",
    "\n",
    "Rename the following list of countries (for use in later questions):\n",
    "\n",
    "```\"Republic of Korea\": \"South Korea\",\n",
    "\"United States of America\": \"United States\",\n",
    "\"United Kingdom of Great Britain and Northern Ireland\": \"United Kingdom\",\n",
    "\"China, Hong Kong Special Administrative Region\": \"Hong Kong\"```\n",
    "\n",
    "There are also several countries with numbers and/or parenthesis in their name. Be sure to remove these, \n",
    "\n",
    "e.g. \n",
    "\n",
    "`'Bolivia (Plurinational State of)'` should be `'Bolivia'`, \n",
    "\n",
    "`'Switzerland17'` should be `'Switzerland'`.\n",
    "\n",
    "<br>\n",
    "\n",
    "Next, load the GDP data from the file `world_bank.csv`, which is a csv containing countries' GDP from 1960 to 2015 from [World Bank](http://data.worldbank.org/indicator/NY.GDP.MKTP.CD). Call this DataFrame **GDP**. \n",
    "\n",
    "Make sure to skip the header, and rename the following list of countries:\n",
    "\n",
    "```\"Korea, Rep.\": \"South Korea\", \n",
    "\"Iran, Islamic Rep.\": \"Iran\",\n",
    "\"Hong Kong SAR, China\": \"Hong Kong\"```\n",
    "\n",
    "<br>\n",
    "\n",
    "Finally, load the [Sciamgo Journal and Country Rank data for Energy Engineering and Power Technology](http://www.scimagojr.com/countryrank.php?category=2102) from the file `scimagojr-3.xlsx`, which ranks countries based on their journal contributions in the aforementioned area. Call this DataFrame **ScimEn**.\n",
    "\n",
    "Join the three datasets: GDP, Energy, and ScimEn into a new dataset (using the intersection of country names). Use only the last 10 years (2006-2015) of GDP data and only the top 15 countries by Scimagojr 'Rank' (Rank 1 through 15). \n",
    "\n",
    "The index of this DataFrame should be the name of the country, and the columns should be ['Rank', 'Documents', 'Citable documents', 'Citations', 'Self-citations',\n",
    "       'Citations per document', 'H index', 'Energy Supply',\n",
    "       'Energy Supply per Capita', '% Renewable', '2006', '2007', '2008',\n",
    "       '2009', '2010', '2011', '2012', '2013', '2014', '2015'].\n",
    "\n",
    "*This function should return a DataFrame with 20 columns and 15 entries.*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "umich_part_id": "009",
    "umich_partlist_id": "003"
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Rank</th>\n",
       "      <th>Documents</th>\n",
       "      <th>Citable documents</th>\n",
       "      <th>Citations</th>\n",
       "      <th>Self-citations</th>\n",
       "      <th>Citations per document</th>\n",
       "      <th>H index</th>\n",
       "      <th>Energy Supply</th>\n",
       "      <th>Energy Supply per Capita</th>\n",
       "      <th>% Renewable</th>\n",
       "      <th>2006</th>\n",
       "      <th>2007</th>\n",
       "      <th>2008</th>\n",
       "      <th>2009</th>\n",
       "      <th>2010</th>\n",
       "      <th>2011</th>\n",
       "      <th>2012</th>\n",
       "      <th>2013</th>\n",
       "      <th>2014</th>\n",
       "      <th>2015</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Country</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>Australia</th>\n",
       "      <td>14</td>\n",
       "      <td>8831</td>\n",
       "      <td>8725</td>\n",
       "      <td>90765</td>\n",
       "      <td>15606</td>\n",
       "      <td>10.28</td>\n",
       "      <td>107</td>\n",
       "      <td>5.386000e+09</td>\n",
       "      <td>231</td>\n",
       "      <td>11.810810</td>\n",
       "      <td>1.021939e+12</td>\n",
       "      <td>1.060340e+12</td>\n",
       "      <td>1.099644e+12</td>\n",
       "      <td>1.119654e+12</td>\n",
       "      <td>1.142251e+12</td>\n",
       "      <td>1.169431e+12</td>\n",
       "      <td>1.211913e+12</td>\n",
       "      <td>1.241484e+12</td>\n",
       "      <td>1.272520e+12</td>\n",
       "      <td>1.301251e+12</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Brazil</th>\n",
       "      <td>15</td>\n",
       "      <td>8668</td>\n",
       "      <td>8596</td>\n",
       "      <td>60702</td>\n",
       "      <td>14396</td>\n",
       "      <td>7.00</td>\n",
       "      <td>86</td>\n",
       "      <td>1.214900e+10</td>\n",
       "      <td>59</td>\n",
       "      <td>69.648030</td>\n",
       "      <td>1.845080e+12</td>\n",
       "      <td>1.957118e+12</td>\n",
       "      <td>2.056809e+12</td>\n",
       "      <td>2.054215e+12</td>\n",
       "      <td>2.208872e+12</td>\n",
       "      <td>2.295245e+12</td>\n",
       "      <td>2.339209e+12</td>\n",
       "      <td>2.409740e+12</td>\n",
       "      <td>2.412231e+12</td>\n",
       "      <td>2.319423e+12</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Canada</th>\n",
       "      <td>6</td>\n",
       "      <td>17899</td>\n",
       "      <td>17620</td>\n",
       "      <td>215003</td>\n",
       "      <td>40930</td>\n",
       "      <td>12.01</td>\n",
       "      <td>149</td>\n",
       "      <td>1.043100e+10</td>\n",
       "      <td>296</td>\n",
       "      <td>61.945430</td>\n",
       "      <td>1.564469e+12</td>\n",
       "      <td>1.596740e+12</td>\n",
       "      <td>1.612713e+12</td>\n",
       "      <td>1.565145e+12</td>\n",
       "      <td>1.613406e+12</td>\n",
       "      <td>1.664087e+12</td>\n",
       "      <td>1.693133e+12</td>\n",
       "      <td>1.730688e+12</td>\n",
       "      <td>1.773486e+12</td>\n",
       "      <td>1.792609e+12</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>China</th>\n",
       "      <td>1</td>\n",
       "      <td>127050</td>\n",
       "      <td>126767</td>\n",
       "      <td>597237</td>\n",
       "      <td>411683</td>\n",
       "      <td>4.70</td>\n",
       "      <td>138</td>\n",
       "      <td>1.271910e+11</td>\n",
       "      <td>93</td>\n",
       "      <td>19.754910</td>\n",
       "      <td>3.992331e+12</td>\n",
       "      <td>4.559041e+12</td>\n",
       "      <td>4.997775e+12</td>\n",
       "      <td>5.459247e+12</td>\n",
       "      <td>6.039659e+12</td>\n",
       "      <td>6.612490e+12</td>\n",
       "      <td>7.124978e+12</td>\n",
       "      <td>7.672448e+12</td>\n",
       "      <td>8.230121e+12</td>\n",
       "      <td>8.797999e+12</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>France</th>\n",
       "      <td>9</td>\n",
       "      <td>13153</td>\n",
       "      <td>12973</td>\n",
       "      <td>130632</td>\n",
       "      <td>28601</td>\n",
       "      <td>9.93</td>\n",
       "      <td>114</td>\n",
       "      <td>1.059700e+10</td>\n",
       "      <td>166</td>\n",
       "      <td>17.020280</td>\n",
       "      <td>2.607840e+12</td>\n",
       "      <td>2.669424e+12</td>\n",
       "      <td>2.674637e+12</td>\n",
       "      <td>2.595967e+12</td>\n",
       "      <td>2.646995e+12</td>\n",
       "      <td>2.702032e+12</td>\n",
       "      <td>2.706968e+12</td>\n",
       "      <td>2.722567e+12</td>\n",
       "      <td>2.729632e+12</td>\n",
       "      <td>2.761185e+12</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Germany</th>\n",
       "      <td>7</td>\n",
       "      <td>17027</td>\n",
       "      <td>16831</td>\n",
       "      <td>140566</td>\n",
       "      <td>27426</td>\n",
       "      <td>8.26</td>\n",
       "      <td>126</td>\n",
       "      <td>1.326100e+10</td>\n",
       "      <td>165</td>\n",
       "      <td>17.901530</td>\n",
       "      <td>3.332891e+12</td>\n",
       "      <td>3.441561e+12</td>\n",
       "      <td>3.478809e+12</td>\n",
       "      <td>3.283340e+12</td>\n",
       "      <td>3.417298e+12</td>\n",
       "      <td>3.542371e+12</td>\n",
       "      <td>3.556724e+12</td>\n",
       "      <td>3.567317e+12</td>\n",
       "      <td>3.624386e+12</td>\n",
       "      <td>3.685556e+12</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>India</th>\n",
       "      <td>8</td>\n",
       "      <td>15005</td>\n",
       "      <td>14841</td>\n",
       "      <td>128763</td>\n",
       "      <td>37209</td>\n",
       "      <td>8.58</td>\n",
       "      <td>115</td>\n",
       "      <td>3.319500e+10</td>\n",
       "      <td>26</td>\n",
       "      <td>14.969080</td>\n",
       "      <td>1.265894e+12</td>\n",
       "      <td>1.374865e+12</td>\n",
       "      <td>1.428361e+12</td>\n",
       "      <td>1.549483e+12</td>\n",
       "      <td>1.708459e+12</td>\n",
       "      <td>1.821872e+12</td>\n",
       "      <td>1.924235e+12</td>\n",
       "      <td>2.051982e+12</td>\n",
       "      <td>2.200617e+12</td>\n",
       "      <td>2.367206e+12</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Iran</th>\n",
       "      <td>13</td>\n",
       "      <td>8896</td>\n",
       "      <td>8819</td>\n",
       "      <td>57470</td>\n",
       "      <td>19125</td>\n",
       "      <td>6.46</td>\n",
       "      <td>72</td>\n",
       "      <td>9.172000e+09</td>\n",
       "      <td>119</td>\n",
       "      <td>5.707721</td>\n",
       "      <td>3.895523e+11</td>\n",
       "      <td>4.250646e+11</td>\n",
       "      <td>4.289909e+11</td>\n",
       "      <td>4.389208e+11</td>\n",
       "      <td>4.677902e+11</td>\n",
       "      <td>4.853309e+11</td>\n",
       "      <td>4.532569e+11</td>\n",
       "      <td>4.445926e+11</td>\n",
       "      <td>4.639027e+11</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Italy</th>\n",
       "      <td>11</td>\n",
       "      <td>10964</td>\n",
       "      <td>10794</td>\n",
       "      <td>111850</td>\n",
       "      <td>26661</td>\n",
       "      <td>10.20</td>\n",
       "      <td>106</td>\n",
       "      <td>6.530000e+09</td>\n",
       "      <td>109</td>\n",
       "      <td>33.667230</td>\n",
       "      <td>2.202170e+12</td>\n",
       "      <td>2.234627e+12</td>\n",
       "      <td>2.211154e+12</td>\n",
       "      <td>2.089938e+12</td>\n",
       "      <td>2.125185e+12</td>\n",
       "      <td>2.137439e+12</td>\n",
       "      <td>2.077184e+12</td>\n",
       "      <td>2.040871e+12</td>\n",
       "      <td>2.033868e+12</td>\n",
       "      <td>2.049316e+12</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Japan</th>\n",
       "      <td>3</td>\n",
       "      <td>30504</td>\n",
       "      <td>30287</td>\n",
       "      <td>223024</td>\n",
       "      <td>61554</td>\n",
       "      <td>7.31</td>\n",
       "      <td>134</td>\n",
       "      <td>1.898400e+10</td>\n",
       "      <td>149</td>\n",
       "      <td>10.232820</td>\n",
       "      <td>5.496542e+12</td>\n",
       "      <td>5.617036e+12</td>\n",
       "      <td>5.558527e+12</td>\n",
       "      <td>5.251308e+12</td>\n",
       "      <td>5.498718e+12</td>\n",
       "      <td>5.473738e+12</td>\n",
       "      <td>5.569102e+12</td>\n",
       "      <td>5.644659e+12</td>\n",
       "      <td>5.642884e+12</td>\n",
       "      <td>5.669563e+12</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Russian Federation</th>\n",
       "      <td>5</td>\n",
       "      <td>18534</td>\n",
       "      <td>18301</td>\n",
       "      <td>34266</td>\n",
       "      <td>12422</td>\n",
       "      <td>1.85</td>\n",
       "      <td>57</td>\n",
       "      <td>3.070900e+10</td>\n",
       "      <td>214</td>\n",
       "      <td>17.288680</td>\n",
       "      <td>1.385793e+12</td>\n",
       "      <td>1.504071e+12</td>\n",
       "      <td>1.583004e+12</td>\n",
       "      <td>1.459199e+12</td>\n",
       "      <td>1.524917e+12</td>\n",
       "      <td>1.589943e+12</td>\n",
       "      <td>1.645876e+12</td>\n",
       "      <td>1.666934e+12</td>\n",
       "      <td>1.678709e+12</td>\n",
       "      <td>1.616149e+12</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>South Korea</th>\n",
       "      <td>10</td>\n",
       "      <td>11983</td>\n",
       "      <td>11923</td>\n",
       "      <td>114675</td>\n",
       "      <td>22595</td>\n",
       "      <td>9.57</td>\n",
       "      <td>104</td>\n",
       "      <td>1.100700e+10</td>\n",
       "      <td>221</td>\n",
       "      <td>2.279353</td>\n",
       "      <td>9.410199e+11</td>\n",
       "      <td>9.924316e+11</td>\n",
       "      <td>1.020510e+12</td>\n",
       "      <td>1.027730e+12</td>\n",
       "      <td>1.094499e+12</td>\n",
       "      <td>1.134796e+12</td>\n",
       "      <td>1.160809e+12</td>\n",
       "      <td>1.194429e+12</td>\n",
       "      <td>1.234340e+12</td>\n",
       "      <td>1.266580e+12</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Spain</th>\n",
       "      <td>12</td>\n",
       "      <td>9428</td>\n",
       "      <td>9330</td>\n",
       "      <td>123336</td>\n",
       "      <td>23964</td>\n",
       "      <td>13.08</td>\n",
       "      <td>115</td>\n",
       "      <td>4.923000e+09</td>\n",
       "      <td>106</td>\n",
       "      <td>37.968590</td>\n",
       "      <td>1.414823e+12</td>\n",
       "      <td>1.468146e+12</td>\n",
       "      <td>1.484530e+12</td>\n",
       "      <td>1.431475e+12</td>\n",
       "      <td>1.431673e+12</td>\n",
       "      <td>1.417355e+12</td>\n",
       "      <td>1.380216e+12</td>\n",
       "      <td>1.357139e+12</td>\n",
       "      <td>1.375605e+12</td>\n",
       "      <td>1.419821e+12</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>United Kingdom</th>\n",
       "      <td>4</td>\n",
       "      <td>20944</td>\n",
       "      <td>20357</td>\n",
       "      <td>206091</td>\n",
       "      <td>37874</td>\n",
       "      <td>9.84</td>\n",
       "      <td>139</td>\n",
       "      <td>7.920000e+09</td>\n",
       "      <td>124</td>\n",
       "      <td>10.600470</td>\n",
       "      <td>2.419631e+12</td>\n",
       "      <td>2.482203e+12</td>\n",
       "      <td>2.470614e+12</td>\n",
       "      <td>2.367048e+12</td>\n",
       "      <td>2.403504e+12</td>\n",
       "      <td>2.450911e+12</td>\n",
       "      <td>2.479809e+12</td>\n",
       "      <td>2.533370e+12</td>\n",
       "      <td>2.605643e+12</td>\n",
       "      <td>2.666333e+12</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>United States</th>\n",
       "      <td>2</td>\n",
       "      <td>96661</td>\n",
       "      <td>94747</td>\n",
       "      <td>792274</td>\n",
       "      <td>265436</td>\n",
       "      <td>8.20</td>\n",
       "      <td>230</td>\n",
       "      <td>9.083800e+10</td>\n",
       "      <td>286</td>\n",
       "      <td>11.570980</td>\n",
       "      <td>1.479230e+13</td>\n",
       "      <td>1.505540e+13</td>\n",
       "      <td>1.501149e+13</td>\n",
       "      <td>1.459484e+13</td>\n",
       "      <td>1.496437e+13</td>\n",
       "      <td>1.520402e+13</td>\n",
       "      <td>1.554216e+13</td>\n",
       "      <td>1.577367e+13</td>\n",
       "      <td>1.615662e+13</td>\n",
       "      <td>1.654857e+13</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                    Rank  Documents  Citable documents  Citations  \\\n",
       "Country                                                             \n",
       "Australia             14       8831               8725      90765   \n",
       "Brazil                15       8668               8596      60702   \n",
       "Canada                 6      17899              17620     215003   \n",
       "China                  1     127050             126767     597237   \n",
       "France                 9      13153              12973     130632   \n",
       "Germany                7      17027              16831     140566   \n",
       "India                  8      15005              14841     128763   \n",
       "Iran                  13       8896               8819      57470   \n",
       "Italy                 11      10964              10794     111850   \n",
       "Japan                  3      30504              30287     223024   \n",
       "Russian Federation     5      18534              18301      34266   \n",
       "South Korea           10      11983              11923     114675   \n",
       "Spain                 12       9428               9330     123336   \n",
       "United Kingdom         4      20944              20357     206091   \n",
       "United States          2      96661              94747     792274   \n",
       "\n",
       "                    Self-citations  Citations per document  H index  \\\n",
       "Country                                                               \n",
       "Australia                    15606                   10.28      107   \n",
       "Brazil                       14396                    7.00       86   \n",
       "Canada                       40930                   12.01      149   \n",
       "China                       411683                    4.70      138   \n",
       "France                       28601                    9.93      114   \n",
       "Germany                      27426                    8.26      126   \n",
       "India                        37209                    8.58      115   \n",
       "Iran                         19125                    6.46       72   \n",
       "Italy                        26661                   10.20      106   \n",
       "Japan                        61554                    7.31      134   \n",
       "Russian Federation           12422                    1.85       57   \n",
       "South Korea                  22595                    9.57      104   \n",
       "Spain                        23964                   13.08      115   \n",
       "United Kingdom               37874                    9.84      139   \n",
       "United States               265436                    8.20      230   \n",
       "\n",
       "                    Energy Supply Energy Supply per Capita  % Renewable  \\\n",
       "Country                                                                   \n",
       "Australia            5.386000e+09                      231    11.810810   \n",
       "Brazil               1.214900e+10                       59    69.648030   \n",
       "Canada               1.043100e+10                      296    61.945430   \n",
       "China                1.271910e+11                       93    19.754910   \n",
       "France               1.059700e+10                      166    17.020280   \n",
       "Germany              1.326100e+10                      165    17.901530   \n",
       "India                3.319500e+10                       26    14.969080   \n",
       "Iran                 9.172000e+09                      119     5.707721   \n",
       "Italy                6.530000e+09                      109    33.667230   \n",
       "Japan                1.898400e+10                      149    10.232820   \n",
       "Russian Federation   3.070900e+10                      214    17.288680   \n",
       "South Korea          1.100700e+10                      221     2.279353   \n",
       "Spain                4.923000e+09                      106    37.968590   \n",
       "United Kingdom       7.920000e+09                      124    10.600470   \n",
       "United States        9.083800e+10                      286    11.570980   \n",
       "\n",
       "                            2006          2007          2008          2009  \\\n",
       "Country                                                                      \n",
       "Australia           1.021939e+12  1.060340e+12  1.099644e+12  1.119654e+12   \n",
       "Brazil              1.845080e+12  1.957118e+12  2.056809e+12  2.054215e+12   \n",
       "Canada              1.564469e+12  1.596740e+12  1.612713e+12  1.565145e+12   \n",
       "China               3.992331e+12  4.559041e+12  4.997775e+12  5.459247e+12   \n",
       "France              2.607840e+12  2.669424e+12  2.674637e+12  2.595967e+12   \n",
       "Germany             3.332891e+12  3.441561e+12  3.478809e+12  3.283340e+12   \n",
       "India               1.265894e+12  1.374865e+12  1.428361e+12  1.549483e+12   \n",
       "Iran                3.895523e+11  4.250646e+11  4.289909e+11  4.389208e+11   \n",
       "Italy               2.202170e+12  2.234627e+12  2.211154e+12  2.089938e+12   \n",
       "Japan               5.496542e+12  5.617036e+12  5.558527e+12  5.251308e+12   \n",
       "Russian Federation  1.385793e+12  1.504071e+12  1.583004e+12  1.459199e+12   \n",
       "South Korea         9.410199e+11  9.924316e+11  1.020510e+12  1.027730e+12   \n",
       "Spain               1.414823e+12  1.468146e+12  1.484530e+12  1.431475e+12   \n",
       "United Kingdom      2.419631e+12  2.482203e+12  2.470614e+12  2.367048e+12   \n",
       "United States       1.479230e+13  1.505540e+13  1.501149e+13  1.459484e+13   \n",
       "\n",
       "                            2010          2011          2012          2013  \\\n",
       "Country                                                                      \n",
       "Australia           1.142251e+12  1.169431e+12  1.211913e+12  1.241484e+12   \n",
       "Brazil              2.208872e+12  2.295245e+12  2.339209e+12  2.409740e+12   \n",
       "Canada              1.613406e+12  1.664087e+12  1.693133e+12  1.730688e+12   \n",
       "China               6.039659e+12  6.612490e+12  7.124978e+12  7.672448e+12   \n",
       "France              2.646995e+12  2.702032e+12  2.706968e+12  2.722567e+12   \n",
       "Germany             3.417298e+12  3.542371e+12  3.556724e+12  3.567317e+12   \n",
       "India               1.708459e+12  1.821872e+12  1.924235e+12  2.051982e+12   \n",
       "Iran                4.677902e+11  4.853309e+11  4.532569e+11  4.445926e+11   \n",
       "Italy               2.125185e+12  2.137439e+12  2.077184e+12  2.040871e+12   \n",
       "Japan               5.498718e+12  5.473738e+12  5.569102e+12  5.644659e+12   \n",
       "Russian Federation  1.524917e+12  1.589943e+12  1.645876e+12  1.666934e+12   \n",
       "South Korea         1.094499e+12  1.134796e+12  1.160809e+12  1.194429e+12   \n",
       "Spain               1.431673e+12  1.417355e+12  1.380216e+12  1.357139e+12   \n",
       "United Kingdom      2.403504e+12  2.450911e+12  2.479809e+12  2.533370e+12   \n",
       "United States       1.496437e+13  1.520402e+13  1.554216e+13  1.577367e+13   \n",
       "\n",
       "                            2014          2015  \n",
       "Country                                         \n",
       "Australia           1.272520e+12  1.301251e+12  \n",
       "Brazil              2.412231e+12  2.319423e+12  \n",
       "Canada              1.773486e+12  1.792609e+12  \n",
       "China               8.230121e+12  8.797999e+12  \n",
       "France              2.729632e+12  2.761185e+12  \n",
       "Germany             3.624386e+12  3.685556e+12  \n",
       "India               2.200617e+12  2.367206e+12  \n",
       "Iran                4.639027e+11           NaN  \n",
       "Italy               2.033868e+12  2.049316e+12  \n",
       "Japan               5.642884e+12  5.669563e+12  \n",
       "Russian Federation  1.678709e+12  1.616149e+12  \n",
       "South Korea         1.234340e+12  1.266580e+12  \n",
       "Spain               1.375605e+12  1.419821e+12  \n",
       "United Kingdom      2.605643e+12  2.666333e+12  \n",
       "United States       1.615662e+13  1.654857e+13  "
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def answer_one():\n",
    "    \n",
    "    import pandas as pd\n",
    "    import numpy as np\n",
    "\n",
    "    energy = pd.read_excel('Energy Indicators.xls', skiprows=17, skip_footer=38)\n",
    "    energy = energy.drop(labels=energy.columns[0:2],axis=1)\n",
    "    energy.columns = ['Country', 'Energy Supply', 'Energy Supply per Capita', '% Renewable']\n",
    "    energy['Energy Supply'] = energy['Energy Supply'].replace('...',np.nan) \n",
    "    energy['Energy Supply'] *= 1000000\n",
    "    energy = energy.replace({'Country':{\"Republic of Korea\": \"South Korea\",\n",
    "                               \"United States of America\": \"United States\",\n",
    "                               \"United Kingdom of Great Britain and Northern Ireland\": \"United Kingdom\",\n",
    "                               \"China, Hong Kong Special Administrative Region\": \"Hong Kong\"}},regex=True)\n",
    "    energy['Country'] = energy['Country'].str.replace(r\" \\(.*\\)\",\"\")\n",
    "    energy['Country'] = energy['Country'].str.replace(r\"([0-9]+)$\",\"\")\n",
    "\n",
    "    GDP = pd.read_csv('world_bank.csv',skiprows=4)\n",
    "    GDP = GDP[['Country Name','2006','2007','2008','2009','2010','2011','2012','2013','2014','2015']]\n",
    "    GDP.columns = ['Country','2006','2007','2008','2009','2010','2011','2012','2013','2014','2015']\n",
    "    GDP = GDP.replace({'Country':{\"Korea, Rep.\": \"South Korea\",\n",
    "                                  \"Iran, Islamic Rep.\": \"Iran\",\n",
    "                                  \"Hong Kong SAR, China\": \"Hong Kong\"}},regex=True)\n",
    "\n",
    "    ScimEn = pd.read_excel('scimagojr-3.xlsx')\n",
    "    ScimEn = ScimEn[0:15]\n",
    "\n",
    "    df = pd.merge(ScimEn,energy,how='inner',on='Country')\n",
    "    dF = pd.merge(df,GDP,how='inner',on='Country')\n",
    "    dF = dF.set_index('Country')\n",
    "    \n",
    "    return dF.sort_index()\n",
    "\n",
    "\n",
    "answer_one()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 2 (6.6%)\n",
    "The previous question joined three datasets then reduced this to just the top 15 entries. When you joined the datasets, but before you reduced this to the top 15 items, how many entries did you lose?\n",
    "\n",
    "*This function should return a single number.*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<svg width=\"800\" height=\"300\">\n",
       "  <circle cx=\"150\" cy=\"180\" r=\"80\" fill-opacity=\"0.2\" stroke=\"black\" stroke-width=\"2\" fill=\"blue\" />\n",
       "  <circle cx=\"200\" cy=\"100\" r=\"80\" fill-opacity=\"0.2\" stroke=\"black\" stroke-width=\"2\" fill=\"red\" />\n",
       "  <circle cx=\"100\" cy=\"100\" r=\"80\" fill-opacity=\"0.2\" stroke=\"black\" stroke-width=\"2\" fill=\"green\" />\n",
       "  <line x1=\"150\" y1=\"125\" x2=\"300\" y2=\"150\" stroke=\"black\" stroke-width=\"2\" fill=\"black\" stroke-dasharray=\"5,3\"/>\n",
       "  <text  x=\"300\" y=\"165\" font-family=\"Verdana\" font-size=\"35\">Everything but this!</text>\n",
       "</svg>"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%%HTML\n",
    "<svg width=\"800\" height=\"300\">\n",
    "  <circle cx=\"150\" cy=\"180\" r=\"80\" fill-opacity=\"0.2\" stroke=\"black\" stroke-width=\"2\" fill=\"blue\" />\n",
    "  <circle cx=\"200\" cy=\"100\" r=\"80\" fill-opacity=\"0.2\" stroke=\"black\" stroke-width=\"2\" fill=\"red\" />\n",
    "  <circle cx=\"100\" cy=\"100\" r=\"80\" fill-opacity=\"0.2\" stroke=\"black\" stroke-width=\"2\" fill=\"green\" />\n",
    "  <line x1=\"150\" y1=\"125\" x2=\"300\" y2=\"150\" stroke=\"black\" stroke-width=\"2\" fill=\"black\" stroke-dasharray=\"5,3\"/>\n",
    "  <text  x=\"300\" y=\"165\" font-family=\"Verdana\" font-size=\"35\">Everything but this!</text>\n",
    "</svg>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {
    "umich_part_id": "010",
    "umich_partlist_id": "003"
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "156"
      ]
     },
     "execution_count": 33,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def answer_two():\n",
    "    \n",
    "    import pandas as pd\n",
    "    import numpy as np\n",
    "\n",
    "    energy = pd.read_excel('Energy Indicators.xls', skiprows=17, skip_footer=38)\n",
    "    energy = energy.drop(labels=energy.columns[0:2],axis=1)\n",
    "    energy.columns = ['Country', 'Energy Supply', 'Energy Supply per Capita', '% Renewable']\n",
    "    energy['Energy Supply'] = energy['Energy Supply'].replace('...',np.nan) \n",
    "    energy['Energy Supply'] *= 1000000\n",
    "    energy = energy.replace({'Country':{\"Republic of Korea\": \"South Korea\",\n",
    "                               \"United States of America\": \"United States\",\n",
    "                               \"United Kingdom of Great Britain and Northern Ireland\": \"United Kingdom\",\n",
    "                               \"China, Hong Kong Special Administrative Region\": \"Hong Kong\"}},regex=True)\n",
    "    energy['Country'] = energy['Country'].str.replace(r\" \\(.*\\)\",\"\")\n",
    "    energy['Country'] = energy['Country'].str.replace(r\"([0-9]+)$\",\"\")\n",
    "\n",
    "    GDP = pd.read_csv('world_bank.csv',skiprows=4)\n",
    "    GDP = GDP[['Country Name','2006','2007','2008','2009','2010','2011','2012','2013','2014','2015']]\n",
    "    GDP.columns = ['Country','2006','2007','2008','2009','2010','2011','2012','2013','2014','2015']\n",
    "    GDP = GDP.replace({'Country':{\"Korea, Rep.\": \"South Korea\",\n",
    "                                  \"Iran, Islamic Rep.\": \"Iran\",\n",
    "                                  \"Hong Kong SAR, China\": \"Hong Kong\"}},regex=True)\n",
    "\n",
    "    ScimEn = pd.read_excel('scimagojr-3.xlsx')\n",
    "\n",
    "    df1 = pd.merge(ScimEn,energy,how='inner',on='Country')\n",
    "    \n",
    "    df2 = pd.merge(ScimEn,GDP,how='inner',on='Country')\n",
    "    \n",
    "    df3 = pd.merge(energy,GDP,how='inner',on='Country')\n",
    "    \n",
    "    return len(ScimEn)+len(GDP)+len(energy)-len(df1)-len(df2)-len(df3)\n",
    "\n",
    "answer_two()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Answer the following questions in the context of only the top 15 countries by Scimagojr Rank (aka the DataFrame returned by `answer_one()`)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 3 (6.6%)\n",
    "What is the average GDP over the last 10 years for each country? (exclude missing values from this calculation.)\n",
    "\n",
    "*This function should return a Series named `avgGDP` with 15 countries and their average GDP sorted in descending order.*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 65,
   "metadata": {
    "scrolled": false,
    "umich_part_id": "011",
    "umich_partlist_id": "003"
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Country\n",
       "United States         1.536434e+13\n",
       "China                 6.348609e+12\n",
       "Japan                 5.542208e+12\n",
       "Germany               3.493025e+12\n",
       "France                2.681725e+12\n",
       "United Kingdom        2.487907e+12\n",
       "Brazil                2.189794e+12\n",
       "Italy                 2.120175e+12\n",
       "India                 1.769297e+12\n",
       "Canada                1.660647e+12\n",
       "Russian Federation    1.565459e+12\n",
       "Spain                 1.418078e+12\n",
       "Australia             1.164043e+12\n",
       "South Korea           1.106715e+12\n",
       "Iran                  4.441558e+11\n",
       "Name: avgGDP, dtype: float64"
      ]
     },
     "execution_count": 65,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def answer_three():\n",
    "    \n",
    "    import numpy as np\n",
    "    import pandas as pd\n",
    "    \n",
    "    Top15 = answer_one()\n",
    "    Top15['avgGDP'] = np.nanmean(Top15[[str(x) for x in range(2006,2016)]],axis=1)\n",
    "    avgGDP = Top15['avgGDP'].sort_values(ascending=False)\n",
    "    \n",
    "    return avgGDP\n",
    "\n",
    "answer_three()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 4 (6.6%)\n",
    "By how much had the GDP changed over the 10 year span for the country with the 6th largest average GDP?\n",
    "\n",
    "*This function should return a single number.*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 70,
   "metadata": {
    "scrolled": true,
    "umich_part_id": "012",
    "umich_partlist_id": "003"
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "246702696075.3999"
      ]
     },
     "execution_count": 70,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def answer_four():\n",
    "    \n",
    "    import numpy as np\n",
    "    import pandas as pd\n",
    "    \n",
    "    Top15 = answer_one()\n",
    "    \n",
    "    Top15['avgGDP'] = np.nanmean(Top15[[str(x) for x in range(2006,2016)]],axis=1)\n",
    "    avgGDP = Top15['avgGDP'].sort_values(ascending=False)\n",
    "    \n",
    "    diff = Top15['2015'] - Top15['2006']\n",
    "    \n",
    "    return diff[avgGDP.index[5]]\n",
    "\n",
    "answer_four()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 5 (6.6%)\n",
    "What is the mean `Energy Supply per Capita`?\n",
    "\n",
    "*This function should return a single number.*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 83,
   "metadata": {
    "umich_part_id": "013",
    "umich_partlist_id": "003"
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "157.6"
      ]
     },
     "execution_count": 83,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def answer_five():\n",
    "    Top15 = answer_one()\n",
    "    return float(Top15['Energy Supply per Capita'].mean())\n",
    "\n",
    "answer_five()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 6 (6.6%)\n",
    "What country has the maximum % Renewable and what is the percentage?\n",
    "\n",
    "*This function should return a tuple with the name of the country and the percentage.*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "umich_part_id": "014",
    "umich_partlist_id": "003"
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "('Brazil', 69.648030000000006)"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def answer_six():\n",
    "    Top15 = answer_one()\n",
    "    \n",
    "    return (Top15['% Renewable'].idxmax(),Top15['% Renewable'].max())\n",
    "\n",
    "answer_six()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 7 (6.6%)\n",
    "Create a new column that is the ratio of Self-Citations to Total Citations. \n",
    "What is the maximum value for this new column, and what country has the highest ratio?\n",
    "\n",
    "*This function should return a tuple with the name of the country and the ratio.*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "umich_part_id": "015",
    "umich_partlist_id": "003"
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "('China', 0.68931261793894216)"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def answer_seven():\n",
    "    Top15 = answer_one()\n",
    "    \n",
    "    Top15['Ratio Citations'] = Top15['Self-citations']/Top15['Citations']\n",
    "    \n",
    "    return (Top15['Ratio Citations'].idxmax(),Top15['Ratio Citations'].max())\n",
    "\n",
    "answer_seven()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 8 (6.6%)\n",
    "\n",
    "Create a column that estimates the population using Energy Supply and Energy Supply per capita. \n",
    "What is the third most populous country according to this estimate?\n",
    "\n",
    "*This function should return a single string value.*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 89,
   "metadata": {
    "umich_part_id": "016",
    "umich_partlist_id": "003"
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'United States'"
      ]
     },
     "execution_count": 89,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def answer_eight():\n",
    "    Top15 = answer_one()\n",
    "    \n",
    "    Top15['Pop Estimate'] = Top15['Energy Supply']/Top15['Energy Supply per Capita']\n",
    "    return Top15['Pop Estimate'].sort_values().index[-3]\n",
    "\n",
    "answer_eight()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 9 (6.6%)\n",
    "Create a column that estimates the number of citable documents per person. \n",
    "What is the correlation between the number of citable documents per capita and the energy supply per capita? Use the `.corr()` method, (Pearson's correlation).\n",
    "\n",
    "*This function should return a single number.*\n",
    "\n",
    "*(Optional: Use the built-in function `plot9()` to visualize the relationship between Energy Supply per Capita vs. Citable docs per Capita)*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 116,
   "metadata": {
    "umich_part_id": "017",
    "umich_partlist_id": "003"
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.79400104354429435"
      ]
     },
     "execution_count": 116,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def answer_nine():\n",
    "    Top15 = answer_one()\n",
    "    Top15['Energy Supply per Capita'] = Top15['Energy Supply per Capita'].astype(float)\n",
    "    Top15['PopEst'] = Top15['Energy Supply'] / Top15['Energy Supply per Capita']\n",
    "    Top15['Citable docs per Capita'] = Top15['Citable documents'] / Top15['PopEst']\n",
    "    \n",
    "    dfCorr = Top15[['Citable docs per Capita','Energy Supply per Capita']].corr(method='pearson')\n",
    "    \n",
    "    return dfCorr.iloc[1,0]\n",
    "\n",
    "answer_nine()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 10 (6.6%)\n",
    "Create a new column with a 1 if the country's % Renewable value is at or above the median for all countries in the top 15, and a 0 if the country's % Renewable value is below the median.\n",
    "\n",
    "*This function should return a series named `HighRenew` whose index is the country name sorted in ascending order of rank.*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 120,
   "metadata": {
    "umich_part_id": "018",
    "umich_partlist_id": "003"
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Country\n",
       "Australia             0\n",
       "India                 0\n",
       "Iran                  0\n",
       "Japan                 0\n",
       "South Korea           0\n",
       "United Kingdom        0\n",
       "United States         0\n",
       "Brazil                1\n",
       "Canada                1\n",
       "China                 1\n",
       "France                1\n",
       "Germany               1\n",
       "Italy                 1\n",
       "Russian Federation    1\n",
       "Spain                 1\n",
       "Name: HighRenew, dtype: int64"
      ]
     },
     "execution_count": 120,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def answer_ten():\n",
    "    Top15 = answer_one()\n",
    "    \n",
    "    Top15['HighRenew'] = [1 if x >= Top15['% Renewable'].median() else 0 for x in Top15['% Renewable']]\n",
    "    \n",
    "    return Top15['HighRenew'].sort_values()\n",
    "\n",
    "answer_ten()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 11 (6.6%)\n",
    "Use the following dictionary to group the Countries by Continent, then create a dateframe that displays the sample size (the number of countries in each continent bin), and the sum, mean, and std deviation for the estimated population of each country.\n",
    "\n",
    "```python\n",
    "ContinentDict  = {'China':'Asia', \n",
    "                  'United States':'North America', \n",
    "                  'Japan':'Asia', \n",
    "                  'United Kingdom':'Europe', \n",
    "                  'Russian Federation':'Europe', \n",
    "                  'Canada':'North America', \n",
    "                  'Germany':'Europe', \n",
    "                  'India':'Asia',\n",
    "                  'France':'Europe', \n",
    "                  'South Korea':'Asia', \n",
    "                  'Italy':'Europe', \n",
    "                  'Spain':'Europe', \n",
    "                  'Iran':'Asia',\n",
    "                  'Australia':'Australia', \n",
    "                  'Brazil':'South America'}\n",
    "```\n",
    "\n",
    "*This function should return a DataFrame with index named Continent `['Asia', 'Australia', 'Europe', 'North America', 'South America']` and columns `['size', 'sum', 'mean', 'std']`*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "umich_part_id": "019",
    "umich_partlist_id": "003"
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>size</th>\n",
       "      <th>sum</th>\n",
       "      <th>mean</th>\n",
       "      <th>std</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Continent</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>Asia</th>\n",
       "      <td>5.0</td>\n",
       "      <td>2.898666e+09</td>\n",
       "      <td>5.797333e+08</td>\n",
       "      <td>6.790979e+08</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Australia</th>\n",
       "      <td>1.0</td>\n",
       "      <td>2.331602e+07</td>\n",
       "      <td>2.331602e+07</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Europe</th>\n",
       "      <td>6.0</td>\n",
       "      <td>4.579297e+08</td>\n",
       "      <td>7.632161e+07</td>\n",
       "      <td>3.464767e+07</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>North America</th>\n",
       "      <td>2.0</td>\n",
       "      <td>3.528552e+08</td>\n",
       "      <td>1.764276e+08</td>\n",
       "      <td>1.996696e+08</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>South America</th>\n",
       "      <td>1.0</td>\n",
       "      <td>2.059153e+08</td>\n",
       "      <td>2.059153e+08</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "               size           sum          mean           std\n",
       "Continent                                                    \n",
       "Asia            5.0  2.898666e+09  5.797333e+08  6.790979e+08\n",
       "Australia       1.0  2.331602e+07  2.331602e+07           NaN\n",
       "Europe          6.0  4.579297e+08  7.632161e+07  3.464767e+07\n",
       "North America   2.0  3.528552e+08  1.764276e+08  1.996696e+08\n",
       "South America   1.0  2.059153e+08  2.059153e+08           NaN"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def answer_eleven():\n",
    "    import numpy as np\n",
    "    import pandas as pd\n",
    "    \n",
    "    \n",
    "    Top15 = answer_one()\n",
    "    \n",
    "    Top15['PopEst'] = Top15['Energy Supply'] / Top15['Energy Supply per Capita'].astype('float')\n",
    "    \n",
    "    ContinentDict  = {'China':'Asia', \n",
    "                  'United States':'North America', \n",
    "                  'Japan':'Asia', \n",
    "                  'United Kingdom':'Europe', \n",
    "                  'Russian Federation':'Europe', \n",
    "                  'Canada':'North America', \n",
    "                  'Germany':'Europe', \n",
    "                  'India':'Asia',\n",
    "                  'France':'Europe', \n",
    "                  'South Korea':'Asia', \n",
    "                  'Italy':'Europe', \n",
    "                  'Spain':'Europe', \n",
    "                  'Iran':'Asia',\n",
    "                  'Australia':'Australia', \n",
    "                  'Brazil':'South America'}\n",
    "    \n",
    "    return Top15.groupby(ContinentDict)['PopEst'].agg({'size':len,'sum':np.sum,'mean':np.mean,'std':np.std}).rename_axis('Continent')\n",
    "\n",
    "answer_eleven()\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 12 (6.6%)\n",
    "Cut % Renewable into 5 bins. Group Top15 by the Continent, as well as these new % Renewable bins. How many countries are in each of these groups?\n",
    "\n",
    "*This function should return a __Series__ with a MultiIndex of `Continent`, then the bins for `% Renewable`. Do not include groups with no countries.*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Continent      Bins            \n",
       "Asia           (2.212, 15.753]     4\n",
       "               (15.753, 29.227]    1\n",
       "Australia      (2.212, 15.753]     1\n",
       "Europe         (2.212, 15.753]     1\n",
       "               (15.753, 29.227]    3\n",
       "               (29.227, 42.701]    2\n",
       "North America  (2.212, 15.753]     1\n",
       "               (56.174, 69.648]    1\n",
       "South America  (56.174, 69.648]    1\n",
       "dtype: int64"
      ]
     },
     "execution_count": 37,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def answer_twelve():\n",
    "    \n",
    "    import pandas as pd\n",
    "    import numpy as np\n",
    "    \n",
    "    Top15 = answer_one()\n",
    "    \n",
    "    ContinentDict  = {'China':'Asia', \n",
    "              'United States':'North America', \n",
    "              'Japan':'Asia', \n",
    "              'United Kingdom':'Europe', \n",
    "              'Russian Federation':'Europe', \n",
    "              'Canada':'North America', \n",
    "              'Germany':'Europe', \n",
    "              'India':'Asia',\n",
    "              'France':'Europe', \n",
    "              'South Korea':'Asia', \n",
    "              'Italy':'Europe', \n",
    "              'Spain':'Europe', \n",
    "              'Iran':'Asia',\n",
    "              'Australia':'Australia', \n",
    "              'Brazil':'South America'}\n",
    "    \n",
    "    Top15['Continent'] = [ContinentDict[x] for x in Top15.index]\n",
    "    Top15['Bins'] = pd.cut(Top15['% Renewable'], bins=5)\n",
    "    \n",
    "    serie = Top15.groupby([\"Continent\", \"Bins\"]).size()\n",
    "    \n",
    "    return serie\n",
    "\n",
    "answer_twelve()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 13 (6.6%)\n",
    "Convert the Population Estimate series to a string with thousands separator (using commas). Do not round the results.\n",
    "\n",
    "e.g. 317615384.61538464 -> 317,615,384.61538464\n",
    "\n",
    "*This function should return a Series `PopEst` whose index is the country name and whose values are the population estimate string.*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 69,
   "metadata": {
    "scrolled": true,
    "umich_part_id": "020",
    "umich_partlist_id": "003"
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Country\n",
       "Australia              23,316,017.316017315\n",
       "Brazil                 205,915,254.23728815\n",
       "Canada                  35,239,864.86486486\n",
       "China                 1,367,645,161.2903225\n",
       "France                  63,837,349.39759036\n",
       "Germany                 80,369,696.96969697\n",
       "India                 1,276,730,769.2307692\n",
       "Iran                    77,075,630.25210084\n",
       "Italy                  59,908,256.880733944\n",
       "Japan                  127,409,395.97315437\n",
       "Russian Federation            143,500,000.0\n",
       "South Korea            49,805,429.864253394\n",
       "Spain                    46,443,396.2264151\n",
       "United Kingdom         63,870,967.741935484\n",
       "United States          317,615,384.61538464\n",
       "Name: PopEst, dtype: object"
      ]
     },
     "execution_count": 69,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def answer_thirteen():\n",
    "    \n",
    "    Top15 = answer_one()\n",
    "    \n",
    "    Top15['PopEst'] = (Top15['Energy Supply'] / Top15['Energy Supply per Capita']).astype('float')\n",
    "    Top15['PopEst'] = Top15['PopEst'].apply(lambda x: '{:,}'.format(x))\n",
    "    \n",
    "    return Top15['PopEst']\n",
    "\n",
    "answer_thirteen()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Optional\n",
    "\n",
    "Use the built in function `plot_optional()` to see an example visualization."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "This is an example of a visualization that can be created to help understand the data. This is a bubble chart showing % Renewable vs. Rank. The size of the bubble corresponds to the countries' 2014 GDP, and the color corresponds to the continent.\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA7AAAAF3CAYAAACcz6fnAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xec3FW9//HXmbZ9N50UCKGETmhBpYh0QRFEAbmWi+hP\nrNer15/lXtu1XtTftRdELKCIIk0sNCkKKCUJIQkQSE9IL9vrlPP7YzeRkrJJdnZ2Nq+nj33szne+\nc85n4oNk33NaiDEiSZIkSdJQlyh1AZIkSZIk9YcBVpIkSZJUFgywkiRJkqSyYICVJEmSJJUFA6wk\nSZIkqSwYYCVJkiRJZcEAK0mSJEkqCwZYSZIkSVJZMMBKkiRJksqCAVaSJEmSVBZSpS6gP8aMGROn\nTJlS6jIkSZIkqSQa23tobO8hmQiEIvWRj5FAYMKIKirSgzvWOXPmzA0xxrE7uq8sAuyUKVOYMWNG\nqcuQJEmSpEEVY+TGR5Zz28wV1FelSSaKFV97tXfnCCHw+TcdyX5ja4va1wuFEJb15z6nEEuSJEnS\nEHX/02u5dZDCK0BNRYpYiHzltnlsausuen87ywArSZIkSUPQupYufvG3xdRWpAYlvG5WU5mioyfP\nT+5fSIxx0PrtDwOsJEmSJA0xMUauuncB+UIkkxr82NZQlWL2skYefm79oPe9PQZYSZIkSRpiHpy/\njmdWNlNfVZpti0IIVGeS/PSBRbR2ZktSw9YYYCVJkiRpCIkxcvPjK6hMJwlh8KYOv1RFOkl3Ls+D\nz64rWQ0vZYCVJEmSpCHk2dUtbGjtpvIlR9lc9e4TuPHzb+fGz72N3/33v7Jm4ZwB6e/+n3+FTSsX\nA/Crj7+RztamLc9VppP88YmVFApDYy1sWRyjI0mSJEl7ijueXAXwstHXZKaCS77wKwCWz3uER276\nIW/81FUvuqdQyJNIJHeqv9Mu//Q2n6tMJ2nuyDLv+SamTR65U+0WQ9ECbAjhYOC3L7i0P/A54Lq+\n61OApcAlMcbGYtUhSZIkSeUily8wY8km6iq3H9Wyne1U1NQBsHL+TGbcfg3VDWPYuPw5Lv3Kb7nj\nex+nbdNa8tkepp35Fg479UKWPPE3Hr/t6t5+erop5LO8/eu38fuvvZ8TLvkw4/Y7dKt9xRh5dNGG\n4R1gY4zPAkcDhBCSwErgVuBTwL0xxitDCJ/qe/zJYtUhSZIkSeViVVMnAUhs5dicfE83N37+7eSz\nPXQ0b+D8j/9gy3PrFj/NW750A/VjJwJw2uWfobK2gVxPFzd96XL2n346+x1zCvsdcwoAd//wv5h4\n8DH9qqkyneTZVa27/+YGwGBNIT4DWBRjXBZCuAA4te/6tcADGGAlSZIkieUb2tnW0asvnEK8ZuFc\n7r3mC7zlSzcAMG7/w7aEV4C5f7mRJbMeAKB901qa1i5nfO2RADxxxy9JZio44oyL+1VTJp1gVVMn\nPblCSY70eaHB6v1S4Ia+n/eKMa4G6Ps+bpBqkCRJkqQhbfHaNgrbSrAvMP7AI+lqa6KztXc1ZjpT\nteW5lfNn8vzTj3Hhp3/KJV+8njGTDyaf7QHg+acfZ9Hj9/Kaf+3/GGIiBBIBVjV27OS7GXhFD7Ah\nhAxwPvC7nXzdFSGEGSGEGevXD63DcyVJkiSpGDa0dZPayvThl2pcvZRYKFBZ2/Cy53o626ioqSdd\nUUnj6qWsXTQPgNYNq/nbL7/O2e//KqlM5U7VFUKgrSu3U68phsGYQnwuMCvGuLbv8doQwoQY4+oQ\nwgRgq4cKxRivBq4GmD59+tDYs1mSJEmSiqgnX9jm2a+b18ACECOnv/tzW91xePIRJ/DU/bfy28+9\njRHjJ7PXAUcAMP/hP9Hd3syd3/8EADUjxvD6j357p2ortRD7MTy9Wx2E8Bvgrhjjz/sefwPY+IJN\nnEbFGD+xvTamT58eZ8yYUdQ6JUmSJKnUvvGnp5m7vJHaynSpS3mR1s4s//e8wziqSDsRhxBmxhin\n7+i+ok4hDiFUA2cBt7zg8pXAWSGEBX3PXVnMGiRJkiSpXFSlkxSG4vzTEMgkS7uBExR5CnGMsQMY\n/ZJrG+ndlViSJEmS9AIHjKvl4eeG1h5AMUbyhcjEkVU7vrnISh+hJUmSJEkATB5TQyoxtGJarhCp\nqUjRUJ0pdSkGWEmSJEkaKiaPriFXKFDsvYp2Rnc2zwHjaktdBmCAlSRJkqQho64qzV4NlXRlS7/j\n72a5PBy174hSlwEYYCVJkiRpSDnv6En05IZGgM0XIiEBJ04dW+pSAAOsJEmSJA0pJxw0lkQikBsC\n5662duU44cAxQ2L9KxhgJUmSJGlIqalIccohY2nrypW0js3rcF87bWJJ63ghA6wkSZIkDTFvPn4y\nmXSS7my+ZDU0d+Z4xf6jOXCvobGBExhgJUmSJGnIGVVbwbtfcwAdPfmS7Ejclc1TlUly+Wv2J4Qw\n6P1viwFWkiRJkoagkw8ey1GTR9LcmR3UfvOFSGdPnveefuCQWfu6mQFWkiRJkoagEALvP3MqY+sq\naRmkEFuIkZbOLOceNZHj9x89KH3uDAOsJEmSJA1RDdUZPnvhkYyozhR9JDZfiDS1Zzn98PG84+T9\nhtTU4c0MsJIkSZI0hI2pq+CLF01jfEMVje1Z8oWBXxPb2ZOnpTPLecdM4v+cesCQDK9ggJUkSZKk\nIW9UbQVfveQoXn/0RFo7cwN2xE4hRpo6sqSSCT75hsN5+xAded0sVeoCJEm7LhYKFNaspdDeBvk8\npDMkx4wm0dBQ6tIkSdIAq0gnefvJ+/HKA0fzvbufY2NrN4kE1Fakdjp05vIFWrvyhAAnTh3DO0/Z\nn9rKdJEqHzgGWEkqM7nly+m85y9kZ84ku3ARFPKQSBLoO3A8lyUxahTpww4jc9JJVL7mFBJVVaUu\nW5IkDZCp4+v55tuOZfayRv74xEqeW91KJJJOBirSSVKJ8LJAG2OkO1egO5snEkglAmcdMZ4zjxzP\n5NE1JXonOy+U4kyhnTV9+vQ4Y8aMUpchSSUTY6Tnscdov/7XZOfNgwhUVBCqqgjJ5MvupaeH2NUF\n+Twhk6bqvPOovujNJMePL80bkCRJRbOqsYPHF2/k6ZUtLF7XSkd3nkQisDnCRiL5PIxrqOTgCXUc\nvvcIjt9/FFWZoTOeGUKYGWOcvsP7DLCSNLTlN22i9X+/Sfff/wGpFKG+fqemCcVsltjSAuk0dR/6\nIFVvOG9Ir22RJEm7LsZIc0eWls4s2XyBRAikUwnG1lVQkU7uuIES6W+AHTqRW5L0Ml3/eISWL3+Z\n2NlFGDVql4JnSKcJo0cTu7tp/da36LrvPho+9zmSo0YWoWJJklRKIQRG1GQYUZMpdSlF4S7EkjRE\ndd51N82f/jQQSOxieH2hUFFBGDWa7Ny5NH7o38ivXz8whUqSJA0SA6wkDUGdf/0rLVd+jVBTSxjA\nDZhCCCRGjiK/di2NH/kohcbGAWtbkiSp2AywkjTE5J5/npavfJVQW0uoqChKH4mRI8mvWk3zlV+n\nHPZCkCRJAgOsJA0pMZ+n5Sv/A/lC0cLrZmHUKHoee5Sue+8taj+SJEkDxQArSUNI5x/+QPaZZwgj\nRhS9rxACoaaW1v/9FvlNTiWWJElDnwFWkoaANWvWcOlb3sJhb3s7p8x8nLc+/BCLWluL1t+U224B\nIFRWEru66Przn4vWlyRJ0kAxwEpSicUYufDCCzl533159FWv4qHXnst/HXEk67u7BqX/UFNDx003\nE3O5QelPkiRpV3kOrCSV2P333086neZt+QI9id4Dxo8cMYK2XI43/+2vNGV7yBUKfOrwIzh34iSW\nt7fzLw8/yCtHj+HxjRsZX1XFdSeeRFUyyS+XLOaXSxbTUyiwX00tPzj+FVSnUixrb+d9jz1CPkZO\n32v8lr7bcjkue+xRmjo7KRx8MF/55je54IILSvVHIUmStF2OwEpSic2bN49jjzmGniefJNTWbrle\nmUjwixNO5N4zzuKWU07l83Oe3LJj8OK2Ni4/4EAePPu1NKTT/HHl8wC8fuIk7j79TB4482wOqq/n\n10uXAPCZJ5/gnfsfwN2nn8nYysqX9XH3cdP543vfx8c+9jF3JZYkSUOWI7CSNATEtjYgEJLJf14D\nvjJvLv/YsIFEgDWdnazr7gZgcnUNR/Zt9DRt5EhWtLcDML+lhf95ah4t2R7aczlO7RttfWzjRn72\nqhMBuGTyvnxp7pwX97FuHcl5c1jZ3s7atWsZP/6fo7SSJElDhQFWkkrs8MMP58ZrroFxe73o+s3L\nl7Oxp5u/nHEm6USC4+74E935PAAVyX9OoEmGQFffqOmHZzzGL044iSNGjOA3S5fy8IZ1W+4LW+l7\ncx/3nH4GqY4OXrngWbq6BmftrSRJ0s5yCrEkldjpp59Od2cXv1qxYsu1JzZtYkVHO2MqKkknEjy0\nbh0rOjp22FZbLsdelZVkCwVuWrFsy/VXjB7Nrc/3tn/T8uVbrrdks4ypqCSTTvPQhvUsW7bsZW1K\nkiQNFQZYSSqxEAI3fPSj/HXTJo6/88+8+u67+MYzT3Hm+AnMbtzEWff+hZtXLGdqXd0O2/rkYUdw\n7v33cvGDf2NqXf2W618+6hh+tmghZ9/3F1pz2S3X3zx58pY+blm7lkOmTi3Ke5QkSRoIoRw265g+\nfXqcMWNGqcuQpKLpuOlmWr/3fRJjx5ashsKmTYy97RYS9fU7vlmSJGkAhRBmxhin7+g+R2AlaQgI\n9fWQKt22BLFQgBAI1dUlq0GSJGlHDLCSNASk9t8PEqX7Kzl2dZHadzKhhCFakiRpRwywkjQEpCZP\nhkK+dyS0FDo7SU+bVpq+JUmS+qmoATaEMCKEcFMIYX4I4ZkQwgkhhFEhhHtCCAv6vo8sZg2SVA5C\nJkN66kHEvvNcB7+AQOYoA6wkSRraij0C+x3gzhjjIcBRwDPAp4B7Y4xTgXv7HkvSHq/64osglxv0\nfmM2S6jIUHHCCYPetyRJ0s4oWoANIdQDpwA/BYgx9sQYm4ALgGv7brsWeGOxapCkclJx8kmETIbY\n0zOo/caWFqrOP59QWTmo/UqSJO2sYo7A7g+sB34eQngihHBNCKEG2CvGuBqg7/u4ItYgSWUjVFRQ\nfcnFxNbWQeszZrOQSlJ9wQWD1qckSdKuKmaATQHHAj+KMR4DtLMT04VDCFeEEGaEEGasX7++WDVK\n0pBS89Z/ITlxAoWWlqL3FWMkNjdT+653kZw4oej9SZIk7a5iBtjngedjjI/2Pb6J3kC7NoQwAaDv\n+7qtvTjGeHWMcXqMcfrYsWOLWKYkDR2hooKGT38acrne0dEiik1NpA6a2rv2VpIkqQwULcDGGNcA\nK0IIB/ddOgN4GrgduKzv2mXA74tVgySVo/Rhh1J7xf8hNjcT8/mi9FFoaSHUVNPw+c959qskSSob\nxf6t5d+A60MIGWAxcDm9ofnGEMK7geXAxUWuQZLKTvVb3kKho5OO634JDQ2EdHrA2i40NxMqKxj5\nrW+RmjRpwNqVJEkqtqIG2BjjbGD6Vp46o5j9SlK5CyFQe/k7SdbX0/qjq4iJBKG+nhDCLrcZ83li\nUxPJ8Xsx4mtXkpo8eQArliRJKr5inwMrSdpFIQSqL3ozo358FclJk4ibNhI7O3e6nVgoUGhqIjY3\nUX3Rmxn9s58aXiVJUlly4ZMkDXHpAw9g9E9+TMftf6Dj1zdQaNwEIUGoqYFUaqujsrFQgO5uYmcH\nEMgcewy1l19O+vDDBv8NSJIkDRADrCSVgZBOU/PmN1H9xgvomTWLzt/fTnbuXAqNjZBMEgFiJITQ\n+3M+T3Lvval49flUv/48j8mRJEnDggFWkspISCapOP54Ko4/HoD8pkbyS5ZQaG+HfI6QzpAYO4bU\nlCmEiooSVytJkjSwDLCSVMaSo0aSHDWy1GVIkiQNCjdxkiRJkiSVBQOsJEmSJKksGGAlSZIkSWXB\nACtJkiRJKgsGWEmSJElSWTDASpIkSZLKggFWkiRJklQWDLCSJEmSpLJggJUkSZIklQUDrCRJkiSp\nLBhgJUmSJEllwQArSZIkSSoLBlhJkiRJUlkwwEqSJEmSyoIBVpIkSZJUFgywkiRJkqSyYICVJEmS\nJJUFA6wkSZIkqSwYYCVJkiRJZcEAK0mSJEkqCwZYSZIkSVJZMMBKkiRJksqCAVaSJEmSVBYMsJIk\nSZKksmCAlSRJkiSVBQOsJEmSJKksGGAlSZIkSWXBACtJkiRJKgsGWEmSJElSWTDASpIkSZLKQqqY\njYcQlgKtQB7IxRinhxBGAb8FpgBLgUtijI3FrEOSJEmSVP4GYwT2tBjj0THG6X2PPwXcG2OcCtzb\n91iSJEmSpO0qxRTiC4Br+36+FnhjCWqQJEmSJJWZYgfYCNwdQpgZQrii79peMcbVAH3fxxW5BkmS\nJEnSMFDUNbDASTHGVSGEccA9IYT5/X1hX+C9AmDy5MnFqk+SJEmSVCaKOgIbY1zV930dcCvwCmBt\nCGECQN/3ddt47dUxxukxxuljx44tZpmSJEmSpDJQtAAbQqgJIdRt/hk4G5gH3A5c1nfbZcDvi1WD\nJEmSJGn4KOYU4r2AW0MIm/v5dYzxzhDC48CNIYR3A8uBi4tYgyRJkiRpmChagI0xLgaO2sr1jcAZ\nxepXkiRJkjQ8leIYHUmSJEmSdpoBVpIkSZJUFgywkiRJkqSyYICVJEmSJJUFA6wkSZIkqSwYYCVJ\nkiRJZcEAK0mSJEkqCwZYSZIkSVJZMMBKkiRJksqCAVaSJEmSVBYMsJIkSZKksmCAlSRJkiSVBQOs\nJEmSJKksGGAlSZIkSWXBACtJkiRJKgsGWEmSJElSWTDASpIkSZLKggFWkiRJklQWDLCSJEmSpLJg\ngJUkSZIklQUDrCRJkiSpLBhgJUmSJEllwQArSZIkSSoLBlhJkiRJUlkwwEqSJEmSysIOA2wI4aAQ\nwr0hhHl9j6eFED5T/NIkSZIkSfqn/ozA/gT4TyALEGOcA1xazKIkSZIkSXqp/gTY6hjjYy+5litG\nMZIkSZIkbUt/AuyGEMIBQAQIIVwErC5qVZIkSZIkvUSqH/d8ELgaOCSEsBJYAry9qFVJkiRJkvQS\nOwywMcbFwJkhhBogEWNsLX5ZkiRJkiS92DYDbAjhP7ZxHYAY4zeLVJMkSZIkSS+zvRHYukGrQpIk\nSZKkHdhmgI0xfmEwC5EkSZKGitjdTc+sWWTnP0vP7NnkV6wgZrOERJLE6NGkpx1J+vDDyRx3HMlR\nI0tdrrTH2OEa2BDC/sB3gFfRuxPxP4CP9q2NlSRJkoaN/Nq1dNz+BzpvvY3Y0wO5HFRUECoqIJki\nxkhuzRpyS5bQ+Yc/QAhUnPxqai65mPRhh5a6fGnY688uxL8GfgBc2Pf4UuAG4JX96SCEkARmACtj\njOeFEPYDfgOMAmYB74gx9uxs4ZIkSdJAiYUCnX/4A60/+BHksoTaOhLV1Vu9N6TT0PdczOfpfvBB\nuv/2N6rOeS21H3g/idrawSxd2qP05xzYEGP8ZYwx1/f1K/rOhO2nfweeecHjrwHfijFOBRqBd+9E\nW5IkSdKAKjQ30/gfH6P1298lVFSQGDWakMn067UhmSQxciShoYHOO+9k47++k+yzzxa5YmnPtc0A\nG0IYFUIYBdwfQvhUCGFKCGHfEMIngD/1p/EQwt7A64Fr+h4H4HTgpr5brgXeuDtvQJIkSdpVhaYm\nNv37R8jOmUsYNap3qvAuCMkkiVGjKbS30/jvH6Fn3lMDXKkk2P4U4pn0jrSGvsfvfcFzEfhSP9r/\nNvAJ/rmj8WigKcaY63v8PDCp39VKkiRJAyR2d9P4iU+SX7GCxKhRA9Jmoq6OQns7TZ/4BKOu+hGp\nyZMHpF1JvbY5Ahtj3C/GuH/f95d+7b+jhkMI5wHrYowzX3h5a11t4/VXhBBmhBBmrF+/fodvRJIk\nSdoZbb+6ntyChYQRA7uLcKKmhtiTpfnLXyHmcjt+gaR+688mToQQjgAOAyo3X4sxXreDl50EnB9C\neF3f6+rpHZEdEUJI9Y3C7g2s2tqLY4xXA1cDTJ8+fWfW3EqSJEnblX3uOb761a9y6/p1JEKCRID/\nd+xxHDdq9E61c+eqVTzX2sKHDz7kRddDQwO5BQvp+N1N1PzLpQNZurRHCzFuPxuGED4PnEpvgP0z\ncC7wUIzxon53EsKpwP/t24X4d8DNMcbfhBCuAubEGH+4vddPnz49zpgxo7/dSZIkSdt111vfxqf/\n9CduO/0MKpJJNnZ3ky0UGF9VNWB9xGwWuroYc+vN29zRWFKvEMLMGOP0Hd3XnxHYi4CjgCdijJeH\nEPaib1OmXfRJ4DchhC8DTwA/3Y22JEmSpJ2SW7mKlXPnMqq6iopkEoDRfZs3HXfHn7hg7314eP06\nAH70ilexf20td61axbfmP0NPocDITIYfveKVjKus5DdLlzK7cRNXHnMs/zbjMepSaWY3NrKuq4vP\nHzmN11dX0/XAX6l+3bkle7/ScNKfY3Q6Y4wFIBdCqAfWATtcA/tCMcYHYozn9f28OMb4ihjjgTHG\ni2OM3TtftiRJkrRrOu+8k1NHjGRVZyevuusOPvHELP7+gj1X6tJp7jr9TN51wIF89snZALxyzBju\nOO107jvzLC7cZx++/9zWj8pZ29XFH089jetPOpkvzZsDFRV0/va3g/K+pD1Bf0ZgZ4QQRgA/oXdn\n4jbgsaJWJUmSJBVJzyOPUFNby1/OOItHNqznofXrec+j/+CzR0wD4MK99wHgTftM5nNzngRgVWcn\n73n0EdZ1ddJTKDC5pmarbZ87cSKJEDi4vp713d2E6mpyy1dQaGsjUVs7OG9QGsZ2GGBjjB/o+/Gq\nEMKdQH2McU5xy5IkSZIGXszlyC1ZQqirJxECJ40dx0ljx3FYfQO/XbYUgBD+eXDG5p/+a/YTvG/q\nQZwzcSIPr1/HN55+eqvtVySS/+wrRkIIxFSK3OLFZKZNK9bbkvYYO5xCHEK4LoTwnhDCITHGpYZX\nSZIklav8mjVAYFFHB4tbW7dcn9fcxN59Gy39/vkVANz2/Aqm9+1K3JLNMqFvg6ffLlu2k53myC9b\nvvvFS+rXFOJfACcD3wsh7A/MBv4WY/xOMQuTJEmSBlrs6oJEoD2X4z9nP0FLtodkSLBfbS3/e+xx\n3LNmNd35POfcdy8FIle94lUAfPyww3j3I/9gQlUVx40axfL29v73mcv39itpt+3wGB2AEEISOB44\nDXgfvRs7HbL9Vw0cj9GRJEkaemKMdOY6yRWyEAKZRIbKVGWpy9qu7MKFbHr/B0k0NGz1+ePu+BN3\nn37mll2JB0Jh/XrqPvRBqi/u9ymU0h5nwI7RCSHcC9QA/wAeBI6PMa7b/RIlSZJUTgqxwILG51jY\nuJAFTc+xvHUZPfkeQkhAhAIFatI17Fu/L1NHHMQhow5hn7rJL1pTWmqJ6hqIhcHtNJkkbGPTJ0k7\npz9TiOcAxwFHAM1AUwjhHzHGzqJWJkmSpCGhraeNx9c8xj3L76a1u5VCzJNJZqhIVlCVqtpyX4yR\nfMyzoPE5nt7wFLeHBONrJnD2vmdz1LijqUgO3KjmrkqM3wsSSWI+T0gmX/b8zHNfP/CdJpOkpkwZ\n+HalPVB/diH+KEAIoRa4HPg5MB4o/d9AkiRJKpoYI4+ueZQbn72BnnwPlclK6ivqtnl/CIFUSJFK\npCDd+/qNneu57ulrqVt4C5cdfjmHjBq0VWhbrzGRIH3gAWSXLB2UUdEYI+RyJPebUvS+pD1Bf6YQ\nfwh4Nb2jsMuAn9E7lViSJEnDVGNXI9c/80ue2fQM1akqGiq2vmZ0e0IIVKd7d/btzHXy/Se+w4kT\nT+LCqW9+0cjtYKs48USy8+fDYATYtjbSB00lUVW69ysNJ/2ZQlwFfBOYGWPMFbkeSZIkldiqtpV8\nZ9a3ac+205CpH5A1rFWpKiqSFTy86iEWNS3i34/9CPW7EIoHQuU5r6Xt5z8nFgqExA5Pldw9uSzV\nb3lLcfuQ9iA7/C82xvgNIA28AyCEMDaEsF+xC5MkSdLgW9W2im/O/H9057toqBiY8LpZIiQYUTGC\ndZ3r+N8Z36Clp2XA2t4ZyTFjqDjxJGJTU1H7id3dhMoqKk46saj9SHuSHQbYEMLngU8C/9l3KQ38\nqphFSZIkafA1dzfz3VnfIlvIUZMu3vTa+kwdm7o38f0nvks2ny1aP9tT+973QDJJ7OkpSvsxRmJr\nK7Uf+gBhAI/kkfZ0/ZkzcSFwPtAOEGNcBWx79b4kSZLKToyRX8+/nrZsO7VFDK+b1aXrWNW2kjuX\n/rnofW1Nau+9qX3ve4itLb0bLQ2w2NhI5rhjqTrnnAFvW9qT9SfA9sTe/6ojQAjBQ6wkSZKGmZlr\nZzBvw1zqM4MzThFCoDZdy91L72J5y7JB6fOlqi+8kMz06cRNmwY0xBZaWkg01FP/yU8MqTNwpeGg\nPwH2xhDCj4ERIYT3AH8BflLcsiRJkjRYunJd3DD/11SnqgY1cCUTSRIhwXVPX1uUUdAdCckkDV/8\nAulpRxI3bSQWCrvVXoyRQlMTiaoqRn77WyTHjRugSiVt1p9NnP4fcBNwM3Aw8LkY4/eKXZgkSZIG\nx6y1M+nJd5NJZl723DfO/iYATaubeeqep3fYVtPqZq7+15/2u++adA1r2tewuHlx/wseQImqKkZ+\n7UoqzziD2LiJ2NGxS+3EbJa4aSOpffZh5A9/QGrffQe4UknQvxFYYoz3xBg/HmP8vzHGe4pdlCRJ\nkgZHjJG7l91FOpne7n3Na/oXYHdWCIEA3L/ivgFvu981VFbS8JlPM+KrXyGk0xQ2baLQ0dGvUeHY\n00Nh00Ziezs1l1/OqKuvIjVp4iBULe2ZdngObAjhTcDXgHFA6PuKMcb6ItcmSZKkIlvaspSNnRuo\n28Ha1/uveoCNyzZyzeU/Z9o5R3DQKQdx+5f/SLazdxfh1370TPY+cu8Xvea6D17Paz9yJntN3av3\n8ft/xTkfO5txB754am1tppY562fT1tNK7SCtwd2aihNOYPT1v6Trvvvp+M1vyK9d1xtiQ4BMpvfM\n2BiJ2SykKeIVAAAgAElEQVTkspBMETJpqi+9lOrzziM5cULJapf2FDsMsMDXgTfEGJ8pdjGSJEka\nXEuaF1OgsMO1r6e971QeveExLvn6RQBku7K89ZtvIVWRYtOKTdz2hT/wrmsue9Frjj5vGnPumMtZ\nU/di4/JN5LL5l4VX6D0fFgLLW5dz2OjDB+y97YpEbS3V57+BqjecR27BAnILF9Ezdw65JUt7z3VN\npUiOG0f66KNJH3AA6cMP85gcaRD1J8CuNbxKkiQNTwsbF5AK/fmV8MXyuQJ3f+se1i5cS0gk2LRi\n08vuOfS0Q3j42r9z+gdOY86f5zDt3CO22V6ukGNF64qSB9jNQgikDzqI9EEHUfW6c0tdjqQ+/VkD\nOyOE8NsQwr+EEN60+avolUmSJKnolrQsoSK58yOIj934ODWjqvk/P38X7/rJZeRz+Zfdk65MM2X6\nFJ57cAHP3Defw886bJvtZZJpFjQu2Ok6JO1Z+vNxWz3QAZz9gmsRuKUoFUmSJGlQFGKB5u4mRlSM\n2OG9meoM3R09Wx53t3VTP66OkAjMuWMuMb/1DY+OfsNR/O6TN7PPUXtTVV+1zfbTiTTrO9bu/JuQ\ntEfZYYCNMV4+GIVIkiRpcOULeSD06+zXcQeMJZFMcM07f8a0c4/kuAuP4ebP3MYz989n32P2JV21\n9V2MJxw8nkxNhmmvO3K77YcQyMbcrrwNSXuQ/uxCfBDwI2CvGOMRIYRpwPkxxi8XvTpJkiQVz45z\nKx+/+z8ASKaSvO07l77oufdc+64tP5/2vtcAMGJCA1dc9+4t11s3tEIhsv/x+22/owihPwVJ2qP1\nZw3sT4D/BLIAMcY5wKXbfYUkSZKGvFRIEULo13mnu2LunfP4xXt/yWuuOIWQ2H44LVDYpbW4kvYs\n/VkDWx1jfOwlU0uc3yFJklTmQgjsVb0XLd3NVKQGPjweec4RHHnOtncefqGefA/71O0z4DVIGl76\nMwK7IYRwAL0bNxFCuAhYXdSqJEmSNCgOaDiA7nx3qcsgHwscMOLAUpchaYjrT4D9IPBj4JAQwkrg\nI8D7i1qVJEmSBsUBIw6kOBOId04qJNm7du9SlyFpiOvPLsSLgTNDCDVAIsbYWvyyJEmSNBimjjyI\nQO+ROonQn7GNgZctZEkl0uxTN7kk/UsqH9v9WyqEkAwhjAGIMbYD3SGE94QQnhmU6iRJklRUo6tG\nc8joQ2nPtpesho5sJ6/Z+zWkk1s/ikeSNttmgA0hXApsAuaEEP4aQjgNWAy8DnjbINUnSZKkIjt9\nnzOIMRZtN+LtKcQCACdNOnnQ+5ZUfrY3hfgzwHExxoUhhGOBfwCXxhhvHZzSJEmSNBgOHnUI42sn\nsrZ9DXWZ2kHtu7WnjWP3Oo7RVWMGtV9J5Wl7U4h7YowLAWKMs4AlhldJkqThJxESvPPwy4lEcoXB\nOy2xK9dFVaqSSw66ZND6lFTetjcCOy6E8B8veFz7wscxxm8WryxJkiQNpkm1k3jdfq/nj4tvpyHT\nQAihqP0VYoGufDdXHPk+ajN1Re1L0vCxvRHYnwB1L/h66WNJkiQNI2fv+1oOGXUozT0tRV0PW4gF\nWrpbOWXv1zBt7LSi9SNp+NnmCGyM8Qu703AIoRL4G1DR189NMcbPhxD2A34DjAJmAe+IMfbsTl+S\nJEnafclEkiuOfC/fn/09ljQvpj5TP+Ajsb3htYVXTngVFx90SdFHeiUNL8U87KsbOD3GeBRwNHBO\nCOFVwNeAb8UYpwKNwLuLWIMkSZJ2QkWqkg8d82EOG304zT0tZPPZAWu7K9dFS08rp+5zGm8/7F9L\ndu6spPJVtL81Yq+2vofpvq8InA7c1Hf9WuCNxapBkiRJO68iWcH7j/og/3LIW+kpZGnpad2tKcWF\nWKC5u4VUIsUHjvogFx/8FsOrpF2yvU2cdlsIIQnMBA4EfgAsAppijJu3t3semFTMGiRJkrTzQgic\nPOnVHDrqUG6YfwPzNz0DQE26mlSif79C9uR76Mx1Egi8csIredPUi6hJ1xSzbEnDXL8DbN/036/S\nu6b1GzHG23b0mhhjHjg6hDACuBU4dGu3baO/K4ArACZPntzfMiVJkjSARleN4UPH/Btr29fy4Mq/\n8fCqh+jMdVGIBUKAdCJNggSx7389+R5CSBAIVKUqecP+F/DKia9iRMWIUr8VScNA2NZ0kBDC+Bjj\nmhc8vhF4FxCAv8cYj9ypjkL4PNABfBIYH2PMhRBOAP47xvja7b12+vTpccaMGTvTnSRJkoogm8+y\nqn0lK1tXsqh5IavbVtNT6CEQqEhWsE/dPuw/4gAm1k5ifPV4kolkqUuWVAZCCDNjjNN3dN/2RmCv\nCiHMpHe0tQtoAt4KFICWfhQwFsjGGJtCCFXAmfRu4HQ/cBG9OxFfBvx+R21JkiRpaEgn0+xbP4V9\n66dw4qSTSl2OpD3MNlfPxxjfCMwG/hhCeAfwEXrDazX923hpAnB/CGEO8DhwT4zxj/SOwP5HCGEh\nMBr46e69BUmSJEnSnmC7a2BjjH8IIfwZ+ABwC/CVGOOD/Wk4xjgHOGYr1xcDr9iFWiVJkiRJe7Bt\njsCGEM4PITwE3AfMAy4FLgwh3BBCOGCwCpQkSZIkCbY/Avtl4ASgCvhzjPEV9E79nQp8hd5AK0mS\nJEnSoNhegG2mN6RWAes2X4wxLsDwKkmSJEkaZNucQgxcSO+GTTl6dx+WJEmSJKlktjkCG2PcAHxv\nEGspG7lCjjXta1jVtpJ1HevoynfRk+8mW8iRDEkyyQyVyQrqMvXsXbcPk2onUZ2uLnXZkiRJklTW\ntrsLsSDGyOr2VaxoXcGS5sUsbFrE2o7VBBLEWKBAgUAghEAgbHlNgQIQSIUkuZinIVPPlIb9mTpi\nKnvX7c2+9VPIJDOlfXOSJEmSVEYMsNvQneti1rpZ3L3sLjZ0bAAikUhFsoLadC2JsL3Z1y8WYyRX\nyPHUhrnMWT+bREiSSiQ5edIpvHrSqxlbPa54b0SSJEmShgkD7EusaV/Ng8//jYdXPUw+5kgn0tRl\nagkh7HKbIQTSyTTpZHrLtVwhx33L/8J9y+9l6sipnDH5LA4ddSjJRHIg3oYkSZIkDTsG2D6LmhZx\n28JbWNK8BIjUpGtIJYq3bjWVSNFQ0UCMkcVNi1jYtJCaVDWvnXIup+z9GoOsJEmSJL3EHh9gu3Nd\n3L7odv76/AMkQ4L6TN1ujbburBACtZlaAHryPdz03O94ZPU/eOcR72JCzYRBq0OSJEmShro9OsAu\nalrEL576KY1djdRmakmG0o56ZpIZ0ok0q9tX8T+Pfpk37H8Bp08+w9FYSZIkSWIPDbAvHHVN903l\nHSpCCNRl6sgVcty28FZmrH3c0VhJkiRJAvq/le4wsaZ9NV9+9Ev89fn7qc3UDNnzWXvXyNZvGY19\n6PkHS12SJEmSJJXUHjUCu7xlOd994ttk8z1DatR1WzaPxmYLWW549nrac+2cve9rB3WNriRJkiQN\nFXvMCOyipkV8a9b/ki/ktmyaVC7SiTR16TpuX3gbty28lRhjqUuSJEmSpEG3RwTYZS1L+d4T3yHA\nkJ0yvCPJRJK6ijr+svxu/rj4dkOsJEmSpD3OsA+wq9pW8d1ZveG1KlVV6nJ2SzIkqUvXceeSO7lv\nxb2lLkeSJEmSBtWwDrAtPS18d9a3ycVc2Y68vlQykaQ2U8MtC25m1tqZpS5HkiRJkgbNsA2wMUZ+\nM/8G2rJt1KZrSl3OgEolUlSnqrj+mV/R3N1c6nIkSZIkaVAM2wA7e91snlw/m/pMXalLKYpMMkN3\nvptfz7/e9bCSJEmS9gjDMsC29LRw/fxfUpWqHNZHztRn6pi3fg4z184odSmSJEmSVHTDLsBunjrc\nleuiIllR6nKKKoRAVbqKG+Zf71RiSZIkScPesAuws9c9wZxhPHX4pSqSFXTne5xKLEmSJGnYG1YB\nNlfIccOzv6ZymE8dfqn6TB1PbZjHs43zS12KJEmSJBVNqtQFDKR5G+bSmeugPlNf6lIGVQiBRAjc\nt/xeDhl1aKnLkaQBkS9E1jR1snxjO0vWtdHUmaUnWyASqUglqa9Ks/9etUweXcP4hkpSyWH1mawk\nSdqKYRVg71l2N8mQLHUZJVGTruHpjU+zsXMjo6tGl7ocSdol7V05/r5gPQ88s5ZlG9q3XM8VIqlE\nYPPcmkhvwE2EQAi9jyePrubVB4/j1QePo64qXYryJUlSkQ2bALuqbRXLWpbtMWtfXyoREkDkH6se\n5rwDzt/1hmKE1lXQvAI6NkDnRmh+HlpXQk8rFHJQyEMi1ftVORLqJ/V+VY2G6jEwcj+oNkRL6r+l\n69u4e+5qHnx2PflCJJMM1FSkSCb6txykUIisbuzkVw8v4dd/X8qJB43lnGkT2W9szR61pESSpOFu\n2ATYB5//KxD36F9UqlPVPPD8/bx2yrmkk/0cfehqgdWzYN08WPU4rHsK8j2QSPYG1UKu9+dECkIS\nXvjnG2NfoM1BLEAyDYTex5UNMG4aTDoexh0O44+G1PDeFVrSztvU1s019y9k9vImIFJbkdqlqcCJ\nRKC6IkU1vSOzDz27noeeXc/hkxp47xlTGVPn3z+SJA0HwyLAduY6+cfqv1OTril1KSWVTqZp6W5l\n3sa5HDPu2G3f2Pw8LH8InvsjrH0SQgLyWUhXQaoKKnZzFDtGKGRh5SOw/G+Q6AvTk0+Cqa+DvU+A\n6lG714ekshZj5O/PreeaBxbRnSvQUJUasA8gk4nAiOo0MUaeWtnMx66fyeWnHMBrDh23R3/IKUnS\ncDAsAuyT62aTK+SpSQ+Lt7NbkokEf1l2z8sDbE8HLLobnvgZNC3puzkDVaN6A+xACqG37WTmn9cK\nOVj+MCx7sDfgTjwOjn4nTD4Zkv7/Ju1JWjqz/OgvzzF7WSPVmSQjqouzXjWE3iDbkyvw4/sW8PBz\n6/nAWQcxsiaz4xdLkqQhaVgkh2c2PU3CT9UBqEpVsbx1GdlClnQiDY1LYN5v4enfQa6ndxpv1egX\nTwUeDIkUVI3s/TkWYM1suOPDUFEPR10Gh74RasYObk2SBt3Gtm6+dOs81jV3MqI6PSgjoplUgnQy\nzVMrm/ncTU/y2TceybiGyqL3K0mSBt6wCLCLmxdTkXR9E/Ru5hRIsGH1DCbMvREW/6V3xLOyoTcs\nDgUh8c8wm+2ER78Dj38fjnwbHPeefz4naVjZ2NrN526eQ1NHDyMGeRR082hsY3sPn71pDl+6aJoh\ndpC1Z9tZ1baSFa0rWNaylM5c55YPW6tS1ezfsD+T6iYxsWYS1enqUpcrSRqiyj7AduW6aOzaVJSz\nX9s2tnHPd+9l9fw1JNNJRkxo4KBXT2XBQwu55OsXvez+P115B694y/GM3W/MgNfSXzU97Zy49FHq\n593Zu561GFOEB1K6qvcrn4Unr4N5v4Hp74Npb4eMv8BIw0VzRw9fvG0uzR09NJTwiJv6qjStnVm+\neOtcvnjRNEbV+uFnMbX1tPLo6kd5YMV9NHY3kghJ8jFHMiT7PnANRCKFWGDm2sdJhiT5mGdM1VhO\nm3w6x+/1CsOsJOlFQoyx1DXs0PTp0+OMGTO2+tyipkV8d9a3qM3UDmifMUaue/+vOPKcIzj2jccA\nsHbBWhY8vJBVT6/eaoAtpUShwDGr5/KqFbOIhSzp6nGMq51Q6rJ2Xq4bult6N5I69b/hgLMHf7qz\npAEVY+RLt81j/qpmRlQPjfWnTR1Z9h9XyxfePK3fR/Wo/1a1reTupXcxa91MCjFSmaogk8j0a8p4\njJHufDfd+R6SIcHxE17JWZPPZq+avQahcklSqYQQZsYYp+/ovrIfgV3VtpJ8zA94u8tmLSeRSmwJ\nrwB7Td2LrrZuls5cxs2fuZX1SzYw4eDxnP/Z8wgh8Kt/+zVnfPA0JhwygW+c/U2Ov2g6C/++kFRF\niov+583UjqphwcMLefjav5PP5amqr+L8z72B2lG7t3vyqI5GzlnwAGM6NtKZqiSXypAtdO/uH0Fp\npCogNRayHXDXf8B+p/cGWc+V1U7IFXKs7VjLqraVNHU10VPoIZNIU5OpZVLtJCbUTCSTHBpBak9w\n39NreWZlc9E2a9oVDVUpFq5t5e65qzj3qEmlLmfYyBVy3LX0Tu5c8mcikdpMLcmQ3Kk2QghUpiqp\nTFWSL+R5ZNXfeWz1o5x/wAWcts/pJBM7154kaXgpWoANIewDXAeMBwrA1THG74QQRgG/BaYAS4FL\nYoyNu9rPwqYFJIowRXb9kvWMP2j8Vp9bu2At77nu3dSNqeO6D/yK5+euZJ9pe7/onmxnlkmHT+TU\nK07hvh/ez+w/PMnJl53I3tP25rIfv4MQArP/8CSP/PpRzvzQ6btU4+ZR1xNWzCQC7elqCIEQIz35\nHgoxlu/mVulqSFXCkgdg5evgtC86Gqvtysc88zc+w73L72VB07MkSFCIkXzMbZmmmAxJQkhQiHn2\nqZvMmfuexbQxR/X/3GTttPUtXVz7t8XUVAzcMTkDIYRAXUWKXz+8jKP3HcWEEVWlLqnsPd+6gp/P\n+xlrOtZQl64dkKCZTCRpqGggV8hx68JbmLH2cd55+LsYX7P1f58lScNfMUdgc8DHYoyzQgh1wMwQ\nwj3AO4F7Y4xXhhA+BXwK+OSudrK8Zdmgj6RMOHQC9eN619yOO3AczaubXxZgk+kkB554AADjDx7P\nkseXAtC6rpXbPv972ja2kc8WGDGhYZdqqMx28fpn72VS62o6U5UUXvCLQgiBECFbyFJRzqNMIQE1\nY/45GnvYxXDKp8GwoReIMfLomke5beEttGfbSRCoTddu94OtGCNr2lfzi3k/I53McO6U13H65DMc\n2RlgMUauuncB+UIkkxp6a/HTqQSdPXl+cM9zfPHN00g4lXiXzdswl5/MuRqINGTqB/zDilQiRUOm\nnpVtK/naY1/lA0f/G1NHTh3QPiRJ5aFov1HEGFfHGGf1/dwKPANMAi4Aru277VrgjbvTT1e+uygj\nsGOmjGHNc2u2+lzqBefNJpIJCvnCy+5JpBJb/gEPibDlnru/fQ/HvelY3nPtuzn3468l15Pb6dpG\ndTTy1jm3MrF1De3p6heF138KxPjyuspSurr36J+nfwe3vRM6Npa6Ig0RjV2NfH/29/jl09eSK2Sp\nz9RRm9l+eIXeD3mq09XUV9STSiS5bdGtfO3xK1nVtmqQKt8zzF/VwvxVLdRXDd3VKnVVKRava2XO\niqZSl1K2nlw3m6ue/BHpZIraTG3RRtpDCNRn6kiEBN974jvM3/RMUfqRJA1tg/KReAhhCnAM8Ciw\nV4xxNfSGXGDc7rSdLWRJFOFtTDluX/I9eZ64ffaWa6ueWc3y2ct3q93u9m7qxtYBMPeOuTv9+v02\nLePSub+nMttFR6Z6u1NqCwz9Dbr6LZGE6rGwdg7ceBFseLbUFanEFjct4suPfpHnNs2nIVO/y0dp\npRNpGjL1rGlfzZWPfZVZa2cOcKV7rjueXAWBITV1+KVCCCRC4E+zV5a6lLI0f9N8rpl3NVWpykE7\nzq4yVUk6keKqJ3/IkubFg9KnJGnoKHqADSHUAjcDH4kxtuzE664IIcwIIcxYv379Nu/LFXJF+eUo\nhMCbv3ohS2Ys5Ydv+TFXv+MaHvzZQ9SO3r3djl99+cnc+tnbuO6D11PVsHNHAxyx5hne8Ow95EOC\n7vSOzy8cNiOwm4UA1WOgqwlufius3PrO1Br+FjUt5LtPfJsYC9RX7P50xRACdZlaKpIZfjbvGh5f\n89gAVbrn2tTWzcwlm6irGLqjr5vVVaZ46vlm1jV3lbqUstLS3cxP5v6YTCIz6Et5KlOVBAI/fvJH\ndGQ7BrVvSVJpFfUYnRBCGvgjcFeM8Zt9154FTo0xrg4hTAAeiDEevL12tneMzkfu/zBVqcqiTCMe\nSo5eNY9Tlj5CV7qSfD/W6eULeSbWTqJmuJ6f19PWe3bseT+CfU4odTUaRGvaV/O1x68kAFWpgd94\nJ5vP0pnv5ENHf5iDRx0y4O3vKW55fAW/e3Q5FbkWHv7Nt1m7aB4VNXUkU2mOPucd7H/cqaUu8UUa\n27Ocf+wk/uXEKaUupSzEGPnxnB8xb8M8GioG/hz2/mrubuaVE17FOw67rGQ1SJIGRn+P0Sla6gu9\nQyI/BZ7ZHF773A5s/pfmMuD3u9NPKpGiMNxGGl9i2uqnOGXZI3T2M7xuNpSn7e22TC0kM/DH98PK\nx0tdjQZJvpDnF0/9nFwhV5TwCpBOpkkn0vz8qZ85srMb7ntqDVXpwJ3f/wQTDzqGt3/9Vi7+/HWc\n9d4v0964rl9tFAoDf0TattRUJLnv6bWUw9noQ8GsdTOZu2Eu9Zm6ktZRl6njkdWP8PTGp0tahyRp\n8BRzbtdJwDuAuSGEzQtJ/wu4ErgxhPBuYDlw8e50kk6kyced3wipXBy29jlOXfqPl+003B8JhnGA\nBcjUQA/wh/fCG38B46eVuiIV2X3L72VFy4qij/hUpapo7m7mlgU38fbD/rWofQ1H7d05NrZ107Z0\nNolkmsNPe9OW5+rGTODIMy+hUMj/f/buO0yq6nzg+PfcO312trOwu5SlLb2vgIiAIIiKiGCMlWJJ\n1KhBY8HYUBNNDInGaGLsHf1FBRMbioCABQRFRHpZOuyyffrMvef3x8ACUrbNFuB8nsfH3Zk7954Z\ndmbue8973pdv3nmaXWu/w4hG6D58At2GjWfn2uUs++/zuJLSKdq2nvNufZwP/zaVFh17sXfzKtJb\ndaTzGWP49v3nCJQXM+JXD9G8XTf2bv6JL2c+jhEJoVvtnHX1faRktmHt4g/IX7GIaDhIecEO2vYd\nxumX3Myahf+leOcmzrjsVgA2LH6fvdu3UHbZqyS7GzYdVkpJebicnd6d7KjYTlFgH2Ezgi40HBYn\nrTytyE7Iprm7BVat8auwR4wIb62didPiaPQLpZrQsOs2Xl/9Kg+f8UdVSVxRFOUUUG8BrJRyMRwz\nghoRr+M4dAcVkWovrT2hZJXvYcTmRQRrEbwCiJM8rRqIBbGhCvjgevjlu+DJbOwRndR0XadHjx5E\no1Hatm3La6+9RnJyclz2vWzZMl599VWefPLJo97vDXv5YPN/SbC5qzxpfv3mN/EW+bDsX395xsTT\n6XJW9dKBV370I7vX7WHk1BEs2f0Nw1oNp6WnZdUPPIr8/Hy++uorLr/8cqDq53iy2F7kw6IJSnZt\noVmbo68QWbvwv9idCVx8/8sYkTCzHrmOVt0GAlCweTW/fHgmic2yKN+3i7KCHYy68RFSs+7mnYcn\ns2HJHMbd/Sz5Kxby3Ycvc+7NfyElsw3jpj2DplvY8dNSlrz3T0b/5s8A7Nu2nl9Mfw3damXm7y+h\nx9mX0GHASP7v/pcZ+Iub0S0W1n35IXmX3s7WIl+DBbDFwWK+2rmYRTsXxVpACY2oGUXXNMT+r09T\nmggEQmhITDom5zK89dl0Se3SaMHayn0/EIwGSGzE1OFDOS1OysMVrC1eQ7f07o09HEVRFKWeNf3q\nGlVondiG7wuWN1j1w4aSGKzggrWfEtUsNUobhtjVfIlsElfqG4TdA4Fi+PBGmPAmWOsntVQBp9PJ\nihWxhIpJkybx9NNPc88998Rl33l5eeTlHXvZw7d7lmJIA4tWvY+tC+8fQ2bn2l/Q0ISGBBbt+ILL\nulxxzO2i0SgWy9HHlJ+fz5tvvlkZwFb1HE8WW/f5iJpHpuIufO0x9mz4Ac1ixZPWgqIdG9m0bB4A\n4YCXsr3b0CxWMtp1JbFZVuXjEtMzSWvZAYDUrHZkdzkNIQRp2R2o2Ld7/+N9zHv+IUoLtiMQmMbB\nzJyWXU/D7ooV4EvJbEvFvt0kpDYnu0seW39YTEpWDqYRxdOiHfmFXnq1Tqm31wZgp3cn72+ctT/t\nVeK0OEm0eaq8MGNKk02lG9lQsgGX1cU5OecwtOVZDRrISimZkz+n8n346NDHaNauWeX9Fz8yvtb9\nzetCQ/D5trkqgFUURTkFnPABbMeUjizfe3KtgbQaEcaunYPVjBKoRTAmkdh0G9rJvAb25xwpULQB\nPr8HzvnrcdsLKfFx+umns3LlSgAWLFjAjBkz+OCDDwC46aabyMvLY/LkyUybNo3//ve/WCwWRo0a\nxYwZM/jPf/7Dgw8+iK7rJCUlsXDhwsP2sXTpUqZOnUogEMDpdPLCiy/w2b5P2fjZJrZ8vZVoMELJ\nrlI6ndmR4TeeVe0xr5rzE9++uxwjYpDVNZPRt41C0zV++HAlX7/+De60BNJapaDbYh+Nmk9w36/v\nZ4b/b2hC44knnuCMM85g+vTp7Nq1i/z8fNLT03nkkUe46qqr8Pl8ADz11FMMGjSIadOmsWbNGnr3\n7s2kSZPo06dP5XMsLi7m6quvZvPmzbhcLp599ll69uzJ9OnT2bZtG5s3b2bbtm1MnTqVW265Jc7/\nevVr7a5ydE2Qmt2OzcvnV94+5Ko7CVSU8u5Dk0lIa87gK26ndfeBhz1259rlWG2Hf+5p1oMzokLT\n0K37L85ponKd7NJZ/yarcz9G3/wY5ft28d8/31j5GN1y8GKepmmVj+ly5li++/BlkjNz6Dx4DFZd\nsGZnORf2i8/r8HNRM8qn+XP4eMuHQKzydU0KEGpCI8EWC8TDRph3N7zD17u/Zkq3q8lKyK6fQf/M\nTu8Odnl3Vq59tdgtXPvSlGNub0ZNNEv9ZwO5rW7Wl6yn0F9AM1eduvMpiqIoTdwJH8BmubPQxQn/\nNA6SkhGbFpEaKMVnc9dqF6Y0cVpO0urDxyIEuNJg4xxo3hP6TG7sEZ3UDMPg888/55prrjnudsXF\nxcyaNYu1a9cihKC0tBSAhx56iDlz5pCdnV1526E6d+7MwoULsVgszJ07l9/d9Ts6/LYtVs1KwYa9\nXP3iFCxWnWeueI68Cf1IbH5kKuP7D31QmUJ8+ROX4i/xsXreGib+8wp0i84nf/2Unz5bTdu8HBa9\nuCYFXqQAACAASURBVJgpz0/GkWDnjVtm0jy3OQCf/2M+vSf04r7L7yclkMo555zDmjVrAFi+fDmL\nFy/G6XTi9/v57LPPcDgcbNiwgcsuu4xly5bxpz/96bDAfsGCBZXje+CBB+jTpw+zZ89m3rx5TJw4\nsXJ2e+3atcyfP5+Kigo6derEDTfcgNV64mRUlPjCWDRBdpc8lrz7L1bNf5fuZ00AIBqOtapp3W0g\nP81/l+zOeegWC6V7tuFOaXa83R5X2O+tfPy6xR9W6zHN23fHW1xA4dZ1/PKhNzA1jVJ/uNZjOJ5C\nfwH/XvkMu327SbC6q51JcCw23YZVs7LXt5dHlz7C2HYXcnabkfW+JvWnop9iac3HOc7Kj35k49eb\niIajRIIRfvHoBP5z93sEK4KYhsHQa4eQe2ZHSneX8fYd/0erHi3ZsWonnmYeLn50PFa7leIdJXwy\nYw7+Uj9C0xj/8IWkZKfwzZtLWDN/LdGwQachHRlyzZnAgaKFknUl61QAqyiKcpI74SO/rIRsDBlF\nStnoxSTioV3xVnKLNuOrY/sbp6XqPrEnHaGBMxm+eRzaDIHUdo09opNOIBCgd+/e5Ofn069fP0aO\nHHnc7RMTE3E4HFx77bWcf/75jBkzBoAzzjiDyZMnc8kllzB+/PgjHldWVsakSZPYsGEDQggqAhW0\nlzkA5PRrgyMhtmQgPSeNsr3lRw1gf55CvHruavas28tL170KQDQUxZXiwrZ6F637tMadEnvPdRnR\nmeLtJQDkL9vK3s0FfPfvi0i0JVFeXk5FRQUAY8eOxemMzRRGIhFuuukmVqxYga7rrF+/vsrXcvHi\nxbz77rsADB8+nKKiIsrKygA4//zzsdvt2O12MjIy2Lt3Ly1b1m4dbmMIR2MBjhCC0Tc/xldvPc6K\nj1/D6UnBYncw8Be/oX3eCMqLdvPOgxORSJyeZEbf9JdaH7P3uVcx74UHWfnpTLI7V38Ktf1pIyja\nvh67O5FgxCAcjX9V+13enTy+/G8EjSBJtrr3LT7gQP/iqBll9sb3KA2VcHHuJfX6XbihZAM2/eDF\nlGgoyvNTXgIgOTOJix+JvZ93/rSTa1++GmeiEzNqcvEjF2F32/GX+nnl+tfoODiWEl68o4RxD4zl\nvLvO5b37Z7NuwXq6n9ON/z70P06/ciCdhuQSDcW+4zcv3ULxjhImPzsRJPxn2rtsW7Gd1r1bAbEZ\n6o2lGxicfWa9PX9FURSl8Z3wAazD4iDVkUogGmjwRurx5owEGLlpIWHdVqcUWIHAdpKtCa423QYS\nmHsXTJgJ+gn/J96kHFgDW1ZWxpgxY3j66ae55ZZbsFgsmObBE/9gMDbLZrFYWLp0KZ9//jlvvfUW\nTz31FPPmzeOZZ55hyZIlfPjhh/Tu3bty5vGA++67j7POOotZs2aRn5/PaWechiS2pvJAei/sTwc1\nqhdwSAk9RnfnrOuHHnb7uoXrj3nCL6XksqcvoW16O27Lu/2w+9zugxkSjz/+OM2bN+eHH37ANE0c\njqovIB2tXcuBcdjtB9+/uq4TjZ5YldYPfWbu5HRGXv/Ho243cMKNDJxw42G3ZXfud1gAmpiexaUP\nz6z8ffg19x/1vhYdenD5o+9U3td//PUAdB48hs6Dx1Teft7UQ7u6wZ4NP9Bz1GVArOpgvLvo7PHt\n4W/L/0rUjNRbyxmLZiHRnsiC7fMRQmNCx4vrJYiVUrK1PP+wmhPHSiFum5eDMzF2gUciWfDsQrat\n2I7QBBWFXnzFsXT75MxkmneMZTxkdmpB6Z4yQv4QFfsq6DQkt/IYAFu+3cKWb7fwwtUvAxAJhCne\nUVwZwNp1O5tLN8X9eSuKoihNy0lRprZtUntCRqixh1E3UnLW5q+wGREieu1TBQ8UcLJpJ3YwXyfO\nFChcDT+81tgjOWklJSXx5JNPMmPGDCKRCG3atGH16tWEQiHKysr4/PPPAfB6vZSVlXHeeefxxBNP\nVAaqmzZtYsCAATz00EOkp6ezffv2w/ZfVlZGdnZsTd/LL7+MIY06X6DK6deGtV+sw1cSO3EOlAco\n21NGdtcstn6/DX9ZACNqsGb+usrHtD0thx/f/4k9/j0ARwTah443MzMTTdN47bXXMIzYGkuPx1M5\nY/tzQ4YM4Y033gBiqcXp6ekkJjaNqq51ZdVFk++nGvJX8ObdF6Pb7LTsehoQC15tcVyvGYwGeer7\nvxMxwrittVsSUl2a0Ei0JzJ/+zy+3LW4Xo5RHi4nEA2gi6qLRlmdB9+vP326Gn+pn6tfmMS1L03B\nneoiGo69R3TrwX0JTcQuSB3jT0dKOP3K07n2pSlc+9IUbnjr1/Qe0+vgMTUrRcFiQtFgLZ+hoiiK\nciI4KQLYrmldK2dnTlQ5pTvoULwFfx0r6JrSxGFxnFoFnH5OCHAkw5K/Q9n2qrdXaqVPnz706tWL\nt956i1atWnHJJZfQs2dPrrjiCvr06QNARUUFY8aMoWfPngwdOpTHH38cgDvuuIMePXrQvXt3hgwZ\nQq9evQ7b95133sndd9/NGWecgWEYsSUCdexr3KxtOkOvPZOZt/0fz016kZm3vo23yEtCegJnTjmD\nV69/jTenvk2L/etfAUb99mz2rtvLS1NepmvXrjzzzDNH3feNN97IK6+8wsCBA1m/fn3l7GzPnj2x\nWCz06tWr8rkfMH36dJYtW0bPnj2ZNm0ar7zySp2eX1Pislswm3gAa3d5uPzRdzjnxkcrbzOkxGWP\nX0XfWRvfozRUWll4qb5pQsNtcfHO+v9QFNgX9/2XhcpiLX5q+P0S9IVwJbvQLTr5322lbM/xW9/Z\n3XY8zTysWxhLxT+wlrZd/7as/HAl4f3rlCsKKyovSEEsg0EXOuXho180UhRFUU4OoqlfJQfIy8uT\ny5YtO+b9wWiQaYvuxK7b6lwYozFopsnEFf/BGQkSttRtlilqRsl0Z+FpoBOmJs2/D9qNhNF/q3pb\npUn789JHKfDvxdEIa7sN0yBqRpkx7PGqN1YAmPXtNt5eso3UBuqnGi8lvjAX9GnJ5Wfk1Hlf60vW\n8eR3T+CxeWpUaTgeysMVtE9qz2/73hrXVOLNpZv4+3dPkHBIgcG/jPobd3x622HbHeilfM6tsTXy\n/lI//5n2LkbUpHmHDHas2sEv/3IJAP931zv86tVYMbhvZi4hHIgw5OrBFG8v5uMZc/CXBtAtGhc9\nPI6UrGSW/mcZP3zwAwA2p42x940hJftg2yNv2Mu0AfeQ6VY9wRVFUU40QojlUsoq+w2eeNHeUTgs\nDgZlncHCHQtIsjd8/7m6yt23EU/Ii99Wt8JNUpqxK/D1nKp2wnCmwqbPoHAtNOt81E2klJiFhUQ3\nbSKyfgPRzZuRfj8yFAIpEXY7wuFAz8rC2qUzlnbt0Fu2ROgN13dRAZfVhSGNRjm2IQ3sp2JRtDpo\nm+HBop14WSC6Jmjfou4X/6SUvL76NayatcGDVwCPNYGNpRv4oXAFvTP6xG2/R7vc/fPgFaDneT3o\neV6Pyt9dyS4mPXPVUfd5IHgFGHjZgMqfU1ulcsXfLzti+/6/yKP/L45/bnMiXJhXFEVRau+kCGAB\nhrQcwsIdX5xw1YgtRpQzty2t07rXAwxpkupIO7XThw8ltFg68Zd/gQufryyMZRQWEpz7OeFly4is\nXYsMBkHTIBIFqzX284HXUEqQEhkOIyw6IECApX17bL164Th7BJYOHU6ov7kTUbvE9qwtWgON0Ekm\nbITpmJzb8Ac+gbVKc8W9GFJDkEDrtLpfAFxfso6SYDGenxVtOtps5fFs/X4bS2Yu5ZLHLmb94g3s\nyy9i0JUDq3ycEAKLZuGzrZ/GNYC1alaOuUC1iZAcGKeiKIpysjppAtgW7kzaJrVle8W2E2oGsvve\ntTgioTjMvsZOKk7EGeh65UyFnUuQu5YRKbDgf+c9Ql9/vb9aiw3hcKCl1Oy1l4ZBdOs2ouvW4//P\nO1jat8N92WXYzxiEsJ+i1Z/rWevE1uiNtDwgKqN0SOnYKMc+UaW6bditOhHDxKqfGKUWDFOiC0Hz\nxLrPts/fPi92rSuOF7ZyB3ckd3D1/w5dFhdby7eyy7uLrISsuIwh0Z6IIc0me6FYSokpjSMuHCiK\noignl5MmgAU4u81Inv/xucYeRrVppslpu1bE2ubUkSkN3FY31hNwDXB9khLMCi/hxy6jfFUOIBDJ\nyQit9ifVQtcRbje43UgpiW7bTtkf/4iwO3COuxDXxRejp6ZUvSOl2lp6WiEb6cRZFzqtE1s36DFP\ndEIIBrRP44u1e0l2nRjrYL3BKP3bp6HVMfW5PFTGT/tWHbdw09bvt7HoxcU4k5wUbtlHZqcWjL1v\nDEIINi3ZzGdPfo4ryUmL3BaVjzl0XemGLzfy5StfYUQNnIlOxt5/AQmph1+4FUIgkXy5azG/yL2k\nTs/pgCRbEk6LA0MaWETT+66JmBFSHamNslZeURRFaTgnxqXxauqW1h2XxUXwBCmh37psJ45omGhd\ne5XKWA3mFLsKmg5lBgJE1q3F2FOKxV6Anu5AS0mpU/D6c0IINI8HLSUVLBb8M9+i6KqrCMxfoNZh\nxVGKI4UOKR3xRXxVbxxHISNEgtVDu6T2DXrck8GonplA02+nAwfaj8HoXnWfqdxavhUhtCrXvu7d\nsJeRt4zg169dS+muUnb8uJNoKMpHj33CJX+awFVPX4G32HvUx7bs2ZJJ/76Ka16cQtcRXfjmzSVH\n3c6pO1hTtLrOz+kAIQStE9s02bZ1ISNE26R2jT0MRVEUpZ6dVAGsRbNwWefLCRohTGk29nCq1HfX\nSsw6tgYBiEoDtzUBp7VuacgnC2lKorv3EF23DsJhsNhAE9iTdtbrcYXNhpaWBgjKH3qIsnvvwygu\nqddjnkpGtD67wVuzBKNBzm4zslEK8ZzoctLdtEpz4Q83TvGtmghEDJon2sltUffU020V2zDMaJXb\nZXbJJDEjEaEJMjpkULa7jKJtRSRnJpHaKhUhBN1HdTvqYysKKnjrd//Hc5Ne4JuZS9m35egtc2y6\njQJ/AREjUqfndKjclE6E47i/eDKkSYdkle6vKIpysjvpzsp6NetN34y+eMNHv3LdVCQHysgu30PQ\nUrc1k+b+ysPNXc3jEAqf+A7Mupp794Buif0HSMOCI3kbiPo/mRZOJyI1jdDX31B01UQ1GxsnXVK7\nkuxIxh/xN8jxwkYYm26jf4sBVW+sHEEIwYV9WxIxmv7ffjgqGdu3ZVzS0zeUbMBWjWUhFuvBzBtN\n1zCNAxddqx7Dp098Rr/xfbnulWs4945ziIaPHjALIdCExm7f7mqNvTq6pnVFE1qT+0w70Cu6U+rR\nK84riqIoJ4+TLoAVQvDLTpfisDiadCpx7r5NCOTBare1ISWmNGnuao5FU21djLJyouvX7591tR7+\n2kodhIHVVdQgYxFCoKWmAlD+4EN4X3q5SZ7wSRlCykiTG9vR6JrO5G5TiJiRes+wkFLijwa4JPfS\n465lVI4vr10ayS4rvlDVM5KNxR+OkuCwMLBjelz2VxDYW+squGmt0yjbXUrJzljmxuq5a466XcgX\nwtMsNlv848c/VrFXSVHw6DO0tdEyoRVZCVkEooG47TMefBEfHVM6kuHKaOyhKIqiKPWs6VVhiIME\nm4cru0zk2R//jU23Ncn0v877NtW5eNOB1OEEVXERo6QEY+tW0PRYG5yjkljde4n4Gu4ERzidYLPh\nf/U1ZHk5nltujusa3OoyTS/h8PeEw6sJhZYSDv+ElGUcOtuj65nYbL2x2fKw2Tpjs/VEiKZ1YaR9\ncgeGthzGgh0LSLIl1ltBp/JwOZ1TuzAgs+qWJcqx2SwaN43qxMOzfsRplXUukBRvppSEIiZTR3fB\naYvP12HECNf6O8dit3DunaN5+853cCU5adWzJYWbjww+z5wymFn3zSahmYfsrlmU7i475j4lEK1G\nSnN1CSEY1eYcXvnp5bjtMx5MJCNaj2zsYSiK0sRJKTH37iW6ZQum1wdSIpxOLK1bobdsidCb1nmP\ncnQnZQAL0LNZL/pm9OX7gu9Jsic29nAOkxgsJylYjt/qrPU+VOrwQUZpaSx41fVY79djkIYVW8Je\n/AXdqU6aXrwIXYeUFAKzZ4MQsSC2gSrpRiIb8Hrfwu9/DzAxzSiaZkcIJ0I0rxxHbDbWTyAwl0Dg\nE4SwoGkpJCRMweUai643nQJhYzuMY0v5FraVbyPR5on7a1kR9pJsT2Zi10lNslXIiaZrdhIju2cy\nd9Uekt1Nqz9neSDKmZ0y6N0mfn/fx8tlONADtk2f1rTpc7Cy9Tm3Hgy82g9oR/s3jixE1PO8HvQ8\nrwcAuWd2JPfMxlvr2atZb+wWGyEjhF1v/NZhwWiQBGsCXVK7NPZQFEVpgqSURNetw//eLEKLFiMj\n4dj5ork/m0vTYm0rhMDWry+uCROw9evbKBMOSvWctP8yB1KJPTYP3nDDVi6tSuvS/cWEanlyLCtT\nh1uc8qnDZnkFRn7VwWuMhhAmuuPYsxX1Reg6IiWVwKzZDZJOHIlspKDgCvbuHY/P9zZCuNG0VCyW\nDDQtCSFshwVnQgiEcKDrqeh6BpqWimkGKCv7C7t3D6OkZDqmWVGvY64uu27nN71vJjshm/Jwedxe\nSykl5eEKPNYEpvb7neqpHEeXD8oh2W3FF2w6qcT+UCx1eNKZ8a1aa9WsyOOGsQ1LECtwGE9W3cql\nnS4nEA02+vIDU5qEjBBXdrkK/RT/PlQU5UiRjRsp/tWvKb7xJoJzPwebDS05BS0pCS0lJfZfUhJa\ncjIiIYHwsuWUTrubfb+8jNDSpY09fOUYTtoAFmKpxL/tOxWrbm2wwi/V0XnfJoxaftFKKTGkQTNX\nMzyn+No8GQoT3bIlduWsuil7wsTqKqjfgR3r0LqOSEnB/9rrhBYtqpdjSBmlvPw59u6dQDi8Ck1L\nQ9ebIUTNZ740zYWuN0PTEvH53mXPnvMIBr+qh1HXnNvqZmq/2+ia1p2ycDlhI1yn/UXNKGWhcrIT\nsrn9tLtId8ZnPaQS47Dp3HZeFyQQaAJViYMRg6gpue3czrgd8Q3uMlwZca36W1cCQaojNe777dc8\nj+7pPSgPN+6FrYqwl9Na9KdbevdGHYeiKE2LjEbxvvQyxb++nmj+VkRKClpqKsJ67PMhoetoyclo\nKSmYXi+ld91N2aN/xvQ27cKwp6KTOoAFaOHO5Ld9p4IQTaLohGaatPAWENZrHlAcCF7THGmnfM9X\nKSG6bWvshxqkeEjTgs0dv4ImNSV0HeFyUf7YX+LeYica3UVBwSWUlz+BpiWg66lxSYEVwoKupyNl\nkH37rqOk5AGkbPwTdKfFyQ29bmRS18kY0qA0VErErNm4DNOgLFRGIBpgXIeLuCPvLlIcp/Z7q750\naO7hzgu6EjZMgo0YxAYjBqGIye/O60LnrPjPsndMzq3zBZV4OfCdkemue3/bnxNCcHnnK7DptkYr\nmBiIBnBb3Vyce0mjHF9RlKZJhsOUTX8Q36uvIjyJaElJNT4f0txuREoKwU/nUPLbqZilpfU0WqU2\nTvoAFqCVpzW39JkKCHyRxk0nTg7G0ldlDYt8HDgRSXWkkupMq4+hnVDMon1Inw8sNZs9kaaObvMC\njdcnWDidSH+Aiscfj1v6XSSyhYKCS4lENqJpzRCibgXCjkbTEtC0NHy+dygquhnTbPwq30II+mcO\n4P7TH2Rkm1FETYPyUDkV4QrCRviI11dKScSI4A17KQuVE4gGOT3rDH4/4D5G5oxSKYj1rHvLZKZd\n0I2IKfE3QmVifzhKOGpy23md6ZMT/1lJgFaJrdHjnLJbW2EzTLqzWbXa+tRGkj2J63r8ipARbvCg\nPRQNYZgmv+p5PW6ru0GPrShK0yUNg7KH/0Bo8ZeI1DREDc8TDyU0DZGaRnRLPiW/u13NxDYhp0QA\nC9A2qS239bsdm26nohF7xKb7i49f5eMopDRjacPOZqQ500/5ok0yFMbYuTO27rXGBBLQbI2bUi6S\nkwktXhyXVOJodAeFhRMxzXJ0Pa1eCw8JoaNpzQgGF1Fc/NsmMRMLkGxPZlyH8fzpzMeY1G0KvZv1\nwWFxUh4upyLsxRv24g37KA+Xo2k6XdK68stOl/LomX/mii5X0sLdorGfwimjR6tk7r+oO1aLTqk/\ngtkAayillJT6Iuiaxj0Xdqdf2/q7CNgmsQ2mNOu91VN1BKOhei9s1CWtK9d0v5ZANEjICNXrsQ4I\nRoOEzTC/7nUD7ZPbN8gxFUU5MfhnzSa0aDEiLT7nQ0IIREoK0c2bqXjq6TiMUIkH0dgFGKojLy9P\nLlu2LC77KvAX8PT3T1IULCLBmtDgMy5nbvmaPntW47O5qt54/6wrCDJcGaqoDPtbQmzYgPT5azz7\neoDQw/j39iDszYzv4GpIBgIgIP21V9FSape2appe9u69CMMobNBKwVJKTLMQl+siUlP/0GDHralA\nNEBFuIKoGUEXFtxWt+rr2kR4gxFeWbSFxesKcFh1nLb6+SwOhg38EYMB7dO5emg7klz1Mxt5qKdX\n/IN1xWvxNGKLswNFyab1/z0tPS3r/XgrC3/g+R+fQ0C9vceklHgjPjShcUOvG+mU2rlejqMoStM1\na9Ysxo8fz5o1a+jc+fDPgOiOHRRdfQ3C4UTYjv5Z/9HOnbT3eOiUWLMOJTO3bGbFnj38a/YsXlq+\nHJfLxcSJE2v9PJSjE0Isl1LmVbXdKTMDe0CGK4N7Bt7P8NZn4434GjylOLtiL5FqpJdJaRKVBnbd\nQZvENpXBa/6+AN3vPbwq2vTZW5jx8bbj7m/ZlnJueWMDAAvWlvDVhppX4s25/Wv2VRyZJvbiwt30\nuPdbet73Ld3vXcr738XWmL68eDe7Sqq+Il/d7QCk11ur1OGf7aVRKhH/nHA6kV4f/tnv13ofZWUz\nMIzdDd7mRgiBpqXj988iGKyfglTx4LQ4yXBlkJWQTXN3cxW8NiEJDiu/GZnLXRd0w6JrlPjCBCPx\nWxsbihiU+MIITXDbuZ2ZOrpTgwSvAMNbjWj0SsSBaICWnpYNErxCrHXdHafdSZozndJQGYYZ33XO\nUTNKWbicTHcmd/W/WwWvinKKmjlzJoMHD+att9464r6Kp/8JhnnM4BXg4907WV9eftT7ouaxM2eE\n0MBioXzGX/n1r36lgtdG1jQW6jQwm25jfMcJ9Mnoy0urXqC4AWdjPSEv0eMd55BZ12bOZiQ7khFx\nSBrOa5tIXtvY1aYFa0tJsOsM6lj3Gd0dxUH++MFWvpueR5LLgjcYpbAillb68uI9dM92k5Vy/D6B\n1d0OwCgopM49XKWGZm0aVamFx4P/vfdwX3H5cT9wjyYY/Aaf7x00rX7W8lVFCA0h3BQXT6NFi4/Q\nNJUhoNRc7zYpPHFlP75Yu5f/fb+TMn8EISDBYUGrYfqXlJKKYBRTQqLDwuWDWjGsS3M8zobtP9sp\ntTOJtiQCUT9OS+37fdeWlJKwGWFkm1ENetxWntbc3f8e5uR/zCf5nwDgsSWg1bDmw6EM06icdb2w\n/ThGtD5brVVXlFOU1+vlyy+/ZP78+YwdO5bp06ezYMECZsyYwfvPP094yVJ+v3ULvcvSuDQnh4d/\nXMmc3bvQhcaw5s05PzubObt28XVhIX9bu5oXBw7i1uXLOC0tjaVF+zgnM4v2CR4eX7uGsGmSYrPx\nr/4DyHA4ABAWC+a+Iu67/nqSc3O5/fbbee6553j22WcJh8N06NCB1157DZerGlmWSp2ckgHsAW2T\n2nLvwPv5YPP/mLftczQhSLAm1NsaQiElzmgQn/Xof9gH1k05dAct3C1qVXhj2J++Z0C7ROavLaXU\nH+WFqztxZm4yC9aWMOOT7Tx1ZUeemb8LXRO8/vVe/nFlRzq3cHH9q+vZVhSbBX3i8g6c0TGJIm+E\ny55ZTWFFhP5tPUedTygoj+Bx6CQ4YicUCQ4LCQ4L73xbwLL8Cq54dg1Oq8bX9/blLx9v538righE\nDAZ1SOLfk3J5d1nhEdut3uXntrc24g0apCdYefnazmQm2/n7J1t55rPtWDRBlzQLr4+pWfrHAVIK\ndEvjV6QGEHY7ZkkxoW++wTFkSLUfZ5p+iovvQggHQjTeyZymuTGMQkpL/0xq6iONNg7lxOZ2WDiv\ndzaje2axakcpH63YxcptJeiaRtQ00YTAbtWwaKLy81lKiWFKghETU4JFE0RNk67ZSVzQtyXdWyZh\n0RsnyUgTGld2uYqnVvwDu26vUwBXGxVhLzmJOfTJ6Nugx4VYj9gx7cfSO6MPc/LnsKLwe6Q0cVgc\n2DRbtb5fpZSEjBAhI4wmNPq3GMConFG0cDfusg9FURrX7NmzGT16NLm5uaSmpvLdd99V3hf4+JNY\nZ4r9kxwl4TAf7drJV6NGI4SgLBwmyWbjnKwsRrXI4oKWB7NTyiIR3h96FgCl4TAfnzUcIQSvb9nM\nU+vX8VDPXgcHIQSR1ashNxeA8ePHc9111wFw77338sILL3DzzTfX8yuhnNIBLBw+G/vfTbPZULIB\nkLit7rg3f3dEg0gEHPIFLqXElAYSsAgL6a50kuxJdZp1jZqSpff346Mfinjw/Xzm3tG78r6cdCfX\nn5VFgl3n9nNbA3D5M6u5dVRLBucms60oyDl/XcmaR/rz4Pv5DO6YxP0X5vDhD0U8+8XuI47Vq3UC\nzRNttL3jG0Z0TWF8v3Qu6J3Oxadl8NTnO5nxy/aVM783jcjm/gtzALjq2TV88EPREdtFoiY3v76B\n92/pTrNEG28vKeCed7fw4jWd+fNH21k30YPdYaM0WJcCKRqapWGKjVSLbsH/1ts1CmADgU8xzWJ0\nvfH7lWpaKn7//0hMvAWLRRVDUmpP0wQ9W6fQs3UKwbDB9mI/24p8rNlZxsa9FfhCsSrCUkpsFo0E\nh5XebTx0zkqkTbqbVmkunLam8bXWJa0rAzMHsmT3Nw1avyBiRtCEYFK3KQ0eOB+qpacV1/S4HbfE\nuQAAIABJREFUlvJwOUt2f8P87fMoD5WhC52oNLAI/bDxmdIkKqPowoIpDVIdaVzQbjh5Lfqr1H9F\nUYBY+vDUqVMBuPTSS5k5cybnn38+AKGvvoL9M6UAHosFu65z63fLOLtFJqMyj91ObFzLVpU/7woE\nuG7JNxQEA4RNk9buw6ucC48HY8vmyi4Hq1at4t5776W0tBSv18s555wTt+erHFvT+KZvAtomteW3\nfW+lwF/A4p2LWLRzIb6IH6tmwWlxxmVW1hUJYFbOHpgY+6tUuqwuUuwpuKzuKsPWYwW2h45vfL9m\nAPTL8ZC/r+pWJ3NXl7B618G1wOWBKBWBKAvXlfLeTbHm8Of3SiPFfeSfi64JPvldT77dUsHnq0u4\ndeZGludXMH1c2yO2nb+2lMc+3oY/ZFLsi9At28UFvQ8PwNbt8bNqp4+RM34AwDAhM9mGlJLuqYJJ\nnwYY28FkbIeq042PSQqEJUysJFT913Te5vNx5VeLWTjy4IfaY6t/wm2x8JvcTgiPh8jadUQ3b8HS\n7uDrtmzZMl599VWefPJJFixYgM1mY9CgQbE0yYrnEaLq16B//6/5+ON+pKXZjnn7ypUVXHfdTzz/\nfDd27Qqxfr2Pm29uU+3nF5sBlvj9s0hMvKHy9gNpPR988EG196UoBzhsOh1beOjYwsOIbifmhZEJ\nHX/BT0U/4Yv4GqTViylNvGEfl3S6hAxXRr0frzoSbYmMbDOKkW1GURGuYKd3J9srtrG1PJ9gNEjE\njGDVrDgtTnKS2tLK04osd7YKWhVFOUxRURHz5s1j1apVCCEwDAMhBGPHjsU0DKJbtiA8iYT2r2O1\naBpzzhrBooICZu3YzoubNvLekGFH3bfLcjCT7fcrvuf6jrmMzsriy8IC/rJ69WHbCl0HKZFlsVoq\nkydPZvbs2fTq1YuXX36ZBQsW1MvzVw6nAtifyXBlML7jBM5vN4YfClbw6dZP2evfgyB2cmDTbTVO\nCZNSEjEjmIFiDGlimLE3XaojlSR7Elat+uuz0hIslPgO759Y7IvSNv3gOiu7JRaU6RpEjaoLiZhS\n8vW9fY9aBbQ6cbsQgv7tEunfLpGR3VKY8sK6IwLYYMTgxtfWs+z+frRKczB99haCkSNnUaWEbtlu\nvr738NQ3s8LL+2PdLNoj+WBTmEe+8bNicgoWrTYBaOwxQosgzYYp6nI8QgikaRJcsICEQwLYvLw8\n8vJihdgWLFhAQkICgwYNIhL5kWh0K5pW91Ygq1d7ue66n3jmma706OGhRw8P55xT81ldTfPg9b6K\nx3MtQjTsekNFaapcVhe39Pktf132F/wRP65jLB+JB1OalIXKOTP7TIa2PKvejlMXHpuHzqmd6awK\nMCmKUkPvvPMOEydO5N///nflbUOHDgVg9apVhHI7ETFNFhXsZUBaOt5olEA0ytmZmfRLS2PAJx8B\nkGCx4o0euwVgeSRCpjN2Tv321q3H3M4sKQGgoqKCzMxMIpEIb7zxBtnZ2XV+rkrVTrkqxNVl1+30\nzxzAPQPu5Z4B9zKp2xSGtRpOhqs5voiPinAF5eEKSoIllIZKKQuVUR4upzxUTlmojJJgCWWhMrxh\nL+XhchwWJ12SOpJiT6GlpyXtktqT7kyvUfAKsTWmmck2Pl8de+MUeyN88mMxg3Orn6LmcVioCB6s\nEDmqWypPfb6z8vcV2yoAGNIpmTe+3gvAxyuLjgicAXaVhPguv+KQx3ppk24/4jgHgtV0jxVvMMo7\nywqPOp5OmS4KK8J8vTF2ZSsSNflppw/D52O712BYaxuPDnFTFjLxhmtW5TP1yYPHRAoQTaOF1Lgv\nFvDw1nyG3H03ubm5LNrfG3bBggWMGTOG/Px8nnnmGR5//HF69+7N3Ll/p6gownXX/cS55y7n3HOX\ns3Rp7PUqLo5w6aU/MHLkMu68cx3H65K1caOfq69exZNPdqZPn1ia99tv7+b3v18PwNSpa7j33g1c\ncMF3DBz4DR98UACAaUruvns9w4YtZeLElVx55Uo+/LAc0/Tzv//9g86dOzN48GDee++9ymMVFxcz\nbtw4evbsycCBA1m5ciUA06dPZ9KkSYwaNYqcnBzee+897rzzTnr06MHo0aOJRJpGn1lFqa2shGx+\n2/dWNE3HW089yA3ToDxUzqCsM7ik86X12gtaURSlMcycOZOLLrrosNsmTJjAm2++yYTRoxn+zdfc\nsHQJ3ZNjXRl8kQhXfrWYoZ99yrgv5vNQr9hyunEtW/H0+nUMn/sZW7xHfibf0bUr13zzNRcsmE/q\nsYprSomMxM6JH374YQYMGMDIkSOPaOuj1B81A1sFIQQt3Jm0cGfSr3lsNswwDQoDBez07qQwUEg4\nGiJkhAibYXRhqZylTbInkeXOIishG4fFAfkLYe0nUMeqlK9e24XfvL6B3721EYAHLsyhfUb193lB\n7zQufvon3v9+H/+4siNPXtGB37y2gZ73fUvUlAzJTeKZSZ144MIcLntmNX0fWMbQzsm0TjsyZTVi\nSG5/exO7SkM4rBrNPDaemRRb2D55cAuuf3V9ZXGm64Zk0uO+b8lJd3Ba24MFmH6+3Tu/6cYtb2yk\nLBAlakimjmpJ26wAU+b4KQv7kRJu6esi2VH76y8SmkwAC2DoOp/06cu3N/+GBx98kLlz51bel5OT\nw/XXX09CQgK33347e/ZcwPXXF3Ldda0ZMCCZHTuCXH75ShYu7M/f/pZP//5J3HZbDnPnFvH660eu\nWz5gypRV/OMfXRgwIPmY2xQUhHn//T5s3Ohn8uQfGTMmg48+KmT79iDz5p3Gvn1hhg79lksvbUEg\nEObGG//A/PlL6NChA7/85S8r9/PAAw/Qp08fZs+ezbx585g4cSIrVqwAYNOmTcyfP5/Vq1dz+umn\n8+677/LYY49x0UUX8eGHHzJu3Lg4vMKK0nhaJ7bh9rw7eeaHf7IvsA9PHKve+yI+oqbBuW3P4/x2\nF6jgVVGUk9LRUnNvueUWACIbN3Ln1m1oiYdP5swZfvYRjxmQns7iUaMrf589dNhh95+blc25WUfO\nol6ak8OlOTkA3NG+PUn7z3FuuOEGbrjhhiO2V+qXCmBrQdf0yqC2RuJ0YtE12838u3of9b4F0/pU\n/pzusZE/43QAhnVOYVjn2FWp3BYuVj582mGPe/vGbkfsKy3Byqe3H6y89vhlHY7Ypk26g3nHGMuE\nvGZMyGtW+fsfJrTjDxPaVbld79YeFt7d57Btwj+uYv6lKXV+Db1hyYT3yygNGoRCC5jWtQfnZmWz\nzefj0sWL6Juayo+lpbT3JPBUXn9cFgsz1qzm0927CBoGp6WmMaNvP4QQjPtiAX1TU/mysICySIQn\n+uUxML3ZEcc81ogPvf38lq2QkQh9WrUiPz//mOOXMkQ0ms/ixRVs2LDh4PPyRvF6o3zzTSkvvBBb\nt3z22WkkJx/7LX7mmSm8+eZuhg1LRdePPsrRo9PRNEFurpvCwths6NKlZYwZ0wxNE2Rk2Bk0KBYA\nb9kiadXKTseOHQG48sorefbZZwFYvHgx7777LgDDhw+nqKiIsv3rR84991ysVis9evTAMAxGj459\nsfTo0eO4r4WinEgy3ZncO+B+PtzyP+ZunYsuNNxWd60DzogRwRf108zZjCndr6ZNYk58B6woinKC\n0NxuOE4P1/gfUEMk1H9dA+XYVApxQxI6jdzb/oQkI1EwonG5AOCwwH/GJvLtZA+zhg3igZU/VFaS\n2+it4Kq27fhi5Cg8Fisvbd4EwDXtO/Dp8LNZOPIcAqbBp7sPzmpGpWTO8LP5Q8/eRyz0PyDFbqc0\nHD7sttJwmFT7wRltu66BEMht24lGj0zVPiAS2QzomKbkf//ry9y5pzF37ml8990gEhJiwWp1X6Y/\n/jEWaN599/pjbmOzHV4xO/b/Y23tAHyV2x3qaLcdOHG3738dNE3DarVW3q5p2nFfC0U50Vh1K+M6\njOd3eXeQk9SW8nAFZaFyomb1/s6llPgiPspDFZiYnN9uDPcMuE8Fr4qinNK0Fi0QugXZUOcMhoml\n7ZHFSpWGowLYhqTbGqLo7UlHBgNxm72WEu5b7KPPSxVMWPA1ewIBCkKxljrZTicD0mMFjC5u3YYl\n+/YBsLiwgNHzPmfoZ3NYXFDAuoryyv2dvz/NpGdKCtv9Po4mwWKhucPJwoLYeuKScJh5e/cwIO3w\nYkkyGiWyaeNR9+HxeKioqCAa3Q5Ihg5N5aWXDq5bXrUqtg554MBk3nsvdpx584ooLT32h7kQ8M9/\ndmXTJj+PPbblmNv9XP/+SXz0USGmKSksDPP116UAdOiQyLZtQTZu/BGIrVc5YMiQIbzxxhtALA0o\nPT2dxMTa9fFVlBNd26S23Nrvd9w38AGGtTqLiBmlPFxBRdhLSbAUb9iLL+KrrKFwaD2FTHc2V/e4\nhkcHP8Z5bc/HqquiaYqinNqEEFhyOyIDgXo/loxEEA47WrMjM+6UhlNvKcRCiBeBMUCBlLL7/ttS\ngbeBHCAfuERKWVJfY2hynCnHm75SjkFGonGbuZ65JsQ+v8nSSW68W0bR7+OPCBmxAlI/T+UTAoKG\nwV3ff8dnw88m2+XisdU/ETQOFsCy67FrQLoQGMf5t33qtP5MW/EdD6yMtQe6vUtX2iYc3iZCWCyY\ne/Ye9fEXXHABF198MbNmvcYDDwgefrgDv//9BkaM+JZoVDJwYBJ//nMnbrsthxtvXM2oUcs4/fRk\nsrOP32rHbtd46aUejB//Pc2a2XC5qr6mdf75zVi8uISzzvqWdu2c9OmTiMdjwenU+dOfsrngggmk\npzdn8ODBrFq1CogVa5oyZQo9e/bE5XLxyiuvVHkcRTnZNXc35+LcXzCh48XsC+xjl28nW8vyKQoW\nETbC6JqOQ3fSJrEN2QnZB+spKIqiKIdxDD+LyKqf6v04sqICx+hzVL2BRiaOltoXlx0LMQTwAq8e\nEsA+BhRLKf8khJgGpEgp76pqX3l5eXLZsmX1Ms4GFSiFFweDW121qQmjqBhj+zaw1G2mIfXJQh4c\n7GZTqcETIxL4cEk3xi/6gmWjzwMg75OP+HDYcE5LS+O25cvo4PFwRU5bTv/0E5afez6GlJw3/3PG\nZLfkzq7dGPfFAqb37EnvlFSKQiFGzZvL8nPPr/X4zLIy7AP6k/yHh4+5jc/3PiUl96LrNW91E08+\nXxS320JxcYTzz1/O++/3ISPDjmEUkZk5F11vGj0oFUVRFEU5+Znl5RRO+AXC7UZY6md+TkqJLCkh\n9d//wrq/3ocSX0KI5VLKvKq2q7cZWCnlQiFEzs9uvhAYtv/nV4AFQJUB7EnDkRSb1pMm1KCP7ClP\n1n1hftSU2HXBZZ0djJ9dyumvltHFtY2OHk/lNrkeD29vzef275bTLiGBye3a47JYuKptO4Z+NodW\nLje9U1LrPJZj0jTk/nTmYxHCRlPIQ5848UfKyqJEIpKpU9uQkXFgplcixPFnfRVFURRFUeJJS0zE\nOXoUgQ8/RqTWz7maLC/HmtsRS4cji5oqDauhqxA3l1LuBpBS7hZCnFrTNEKAIxnMaGw9rFJNdQ/Y\nVu8zaJesk+7SWHRlApFAOr7dfSvv3+bzoQnBjL79jnjs3d26c3e37kfcfmjp9TS7vU6zrwcI/fhv\nSU1LpiksXX/33T5H3CalCQiEUJX5FEVRFEVpWAnXXkvwi4VIvx/hcsV13zISASlJvPNOlT7cBDT+\nmfAxCCF+JYRYJoRYVlhY2NjDiZ/EVhA9/iyb8jNa3T4onv0hwFUfljP9jFhgJTQTI9wEgyzTBMfx\nZy+t1o6AcdSqvo1NygAWSzuEUN25FEVRFEVpWFpSEknTpiEDgbhWJJZSIsvKSJg0EUs7VX24KWjo\nAHavECITYP//C461oZTyWSllnpQyr9nJVOkrsx9E6r9K2klFaNRlFvZXvZz8MCWVkTmxWW8pNYzQ\n4c2uW7vdLBx5Tl1GWXemibAfP4DVtDQ0zQNEGmZMNSClH7v9yBlsRVEURVGUhmAfdDruiVchS0vi\nEsRK00QWFWEfOgTX5ZfFYYRKPDR0APtfYNL+nycB7zfw8Rtf8x5QRZqocjhhj2/7IQEYIU+V2zU4\n08CSm3vcTYQQ2Gw9kdLfQIOqCYHd3ruxB6EoiqIoyinMPXkS7kmTkGWlSH/tz5dkOIwsLsZ+1jCS\n7r0HoetxHKVSF/UWwAohZgJfA52EEDuEENcAfwJGCiE2ACP3/35qSTt+gKIcSTgccWw/FNuPGYnv\n2oi4sFixtm9X5WYOx0ikbKBm3dUk9xfastmqLBynKIqiKIpSb4QQJEyZTPIjfwSLBbO4uEazsdI0\nMUtKkMEAnqm3kHTfvQir6rndlNRnFeJjzbOPqK9jnhCSWsdSYs0oaGomtjqEpoHNDtEoaHW75iK0\nKNGQh6ZQyfdQUkqIRrC0qzqAdblGU1b2R6SMIETT+EA1zRIcjsFYLJmNPRRFURRFURTsp59O+muv\n4H3hRYKffIIZNcBqRTidR7TakYaBDAYhGAQhsA86nYQbbsCSndVIo1eOR0VQDU3TYutgdy+PVSRW\nqkVLcGMWl8QhgDWI+JrgmupwGC01FS0xscpNNc2N2z0er/ctdL2pPBdJQsKkqjdTFEVRFEVpIFpS\nEom33UrCddcS/Hwewfnzia5bjxmNHDynNCUIsLRvj33Q6ThHj0bPOLUapZxoVADbGHLHwM4ljT2K\nE4pwJ0BxcZ33I6VGxN/0PpRkMIitd69qb+92X4bX+39IGW30qr+mWYGut8Bu79+o41AURVEURTka\nzePBNe5CXOMuREqJWVCA6fWBNBFOF3qL5mqN6wlEBbCNodWg2P+ljPWGVaoUl35ewgDT0jQLOBlR\nbKdVf/2o1dqehITL8Xpfb9RZWCkNpAySkvIHhGiyXbkURVEURVGA2BpZvXlz9OaNPRKlttQZZ2NI\nyIDkthDxNfZIThjC6UA4nGAYtd+HFiXsbU6TW/9qGKDrOIYOrdHjEhN/i65nYprl9TSyqplmMW73\nL3A4BjTaGBRFURRFUZRThwpgG0unsaofbA0IQGueAfur3dZuH4Kwr+ldbpNlZTiGD0dLSqp640No\nmpPU1MeQMoSUDd8X1jTL0PV0kpJ+1+DHVhRFURQl/qSUmF4vRnExMhCIFZlUlCZGpRA3ltzzYcnf\nYwGZSr2sFi0pCUPTapd6LQxM00LUn1o/g6slKWOFA1zjL6rV4+32PiQl3UlZ2Z/QtLQGWw9rml5A\nIy3tX2haQoMcU1EURVGU+mHs3o3/w48IzJod652qaWAY6NnZuC67FMdZw9Dc7sYepqIAaga28Xgy\nofWZEChp7JGcMISmoaWng1HzHqiaHiVYmkNT+5OXPh+W1m2w5Na+P3BCwlV4PDdhmsUNMhMbS1k2\nSE9/FputS70fT1EURVGU+EpIiF18NgMBSqc/xL4rr8L/5pug62hpaWgpKYi0NIziYioef5x94yfg\ne3OmmpFVmoSmdTZ/qumt2o7UlJaWBojYLGy1SUAQLs+up1HVjpQSQiFcl1+KqEMxLyEEiYk3kJR0\nB6ZZimlWxHGUB0kpMYx9gIVmzV7Bbu9bL8dRFEVRFKX+mRUVlNx0C6EvvkAkJSNTUhE2W+X9Qgg0\ntxstJRUcTrzPPkfF3x5XQazS6FQA25iyTgN3MwirYk7VpdntaBnNajQLK/QwYW8G0rDX48hqTpaV\nYenUCcfw4XXelxACj2cyzZq9iKZ5MIxCpKx9waufkzKIaRZitw+gefP3sdl6xG3fiqIoiqI0vNJ7\n7+eL5csZv+pHbli2lKGfzQFg4ldfcvbnn3Hmp3N4dfNmAITVSrvFC7n/ySfp2aYNAwcOZO/evY05\nfOUUpgLYxqRpkHcDRPyNPZITit6iBdjs1QxiJUJAsKRtvY+rJmQ4DEKQdM/dce07Zrf3p3nzD3C7\nf4lpluwPZGuecl05Thncv48IKSl/JD39OSyWFnEbr6IoiqIojcA0ifz4IyLRw4qSYu7u1oPFo0YD\n8Pe805g7YiSfjjib5zdtoDgUAsBvGPTLzGRubmfO7N+f5557rjGfgXIKU0WcGlvnC2HZvyBUDjZV\nDKc6hKZhyWlDdP36Kgs6xWZfm2GEalbhtz5JKZEV5STceAOWVq3ivn9Nc5OSch9u98V4vW/g938A\nRBHCjhCu4xZ6iqUFRTBNH2CgaYkkJt6C2z0eXU+P+1gVRVEURWkE0SjoOsKI0ic1lTaHFGh6buMG\nPtq1E4Cdfj+bvV5S7XZsmsY52S2RJSV013QW5ec30uCVU50KYBubboXB0+CTW8Hqrnl13VOU5nKh\nZWRgFhSAxXqMrWKzr4GiTg06tqrEUoc74xo/vl6PY7N1ITX1DyQl3Y7f/wGBwEdEImsxzShCCKQ0\nkPvbEgmhAzqxoDUJh+OM/f1dBzdYZWNFURRFUepfdNu2/2/vzsOkKs+8j3/vWnqlu7EBAQEBkUVZ\nRGwR9ZVNMGociTG4RCOaTBLjEp03TjJqjIyJMxoTs6jxijEK45jkTdTMGBMTkMSYNwMqqwgIKKIg\nyNJANzS91HLPH6dUlEa2qj5d1b/PddVF7fV7upo+5z7Pc54HT6exykrYuoWy6Afb+b9v2cwLmzfx\nh/ETKYvF+NRfn6c5HZySFDML5uwoLSU1bx7JwYc+AaXI4dCeaXtwzCSoPhbq3oKSzmGnyRvRHj1I\n19VDczPE9v5VjkRbaK7vRTrRfqZ9Tzc0YPF41ocOf5xotDMVFZdTUXE57mmSyXUkk6tIp3fg3gxE\nMCsiGu1FPD6IaLRLm+QSERGRtpd4bSVAqxNI1icSdI4XURaLsbq+ngXbavd6jpWWkm7YjTc15Tyr\nSGtUwLYHZjD2VvivqyCdgkjbFDb5ziIR4scOILFqNSQSHy5iLYUTobH22PACfoQ3NkIyQed7v5+T\nocMHwixCPN6XeLxvKJ8vIiIi4fLG3cECDa2Y2L0HM9esYdzsWRxbUcFJ1a0f1DYDT2ZvskiRg2H5\nMBV2TU2Nz58/P+wYuffnb8GKp4KZieWAeXMLidWrg0mdojHAicSaadg8jJb63mHHA4Li1Zua6Pzv\nd1I8enTYcUREJMcSqQSbdr/LrsQu0u4UR4s4sqw7FUUVYUeTDq7xD89S/73vEdlHcXog0nU76DLj\nUWK92tcShZLfzGyBu9fs73nqgW1PTv9nWPs8NO+EYm3gDpQVFxEfNJDE6tch0YIVp0k0VrebdV/T\nDQ1YKkXnu++iuOaksOOIiEiO1LfU8+LGeczbOJdNDZuIWgT4YJhmypOUxzsxtMtQxvYeT9/Kvoe1\nDrjIoYh0qQ5WwjhEnkpB2ol01mlvEg4VsO1JcQVM/i48/Y8QL9NQ4oNgRUERm3xjFSQaaXjnOPbc\naQiDu+M7tmOlpXS+57sUDRsaah4REcmNhkQDv339KV7cOA/P9LZWFHUiYh8uEtydpCd56d2XeOnd\nl+hZ3pNLh3yWYzoPCCm5dERFo0ZhpaV4czNWXHzQr/f6eorHnkGkvP3MMSIdi9aBbW/6jIGhU6Gx\nNlgiRg6YRSPEenUh2W8aqe0J0nV1oWXxpiZ8Wy3Fo0+hy8wZKl5FRArU8tpl3DF3OvM2zKVTvJyq\n4kpKYiV7Fa8QTJoTj8SpKq6ksqiCLY2buXfB93li1W9IpBIhpJeOyOJxyi68EN+166Bf+96ph2UX\n5nYlBZGPowK2Pfo/N0P3EdC4Lewk+cMdGmux4ZdQet19VP/kAaI9upOurcWTyTaM4aS3b4N0mspv\nfpOqf/sO0S6a1VdEpBD9dd3z/GTx/STTCaqKK1stWvfFzCiPl9MpXs5f1s3hx4t+RGOyMYdpRT5Q\n+slzseLiYILJg+B1dcT69yM+VAfmJTwqYNujWBGcez+UHgFN4fUi5pXdtdBjFJxxC5gRHzyYLg//\njLLLL8N37SS9bRve0pKzj/dUKviMHdspPuUUujw2k9IzJ+rcJhGRAvX/1/+NX6/6FeXxckpiJYf8\nPtFIlKqiKtbUvcFPFj+gnlhpE9Fu3aj619vxpsYDLmLTdXVEKjrR+c7vaP9GQqUCtr0q6wLn/RQ8\nDQkdkf1YjTuCmZvP/TFE4+/fbcXFVHzh83T9z8cov/wySKdIb99GeudOsjX7tjc2kt5Wi+/aScnk\nSVQ/cD9Vd6rXVUSkkK3fuZ7/lyleY5HDn07EzKgqqmRN3ev87o2ns5BQZP+KTzmFqm9/G5LJYF8m\n0frBE29sJF1bS7RLF464/z6i3bu3cVKRD1MB2551GwJn/wASDZDUYtGtaq4Peqz/4adBj3Uroj16\n0OnzV9HtySeo+tZtxPv3w7dvJ719O+naWtK7d+Pp9H4/yt3xpqYPXrdjB1ZaSqdrrqHbE7+h6uZ/\nIT5kiI5Kioi0Y3feeSdDhw5lxIgRjBw5khdffPGgXp9MJ5mx7BE2LN7ApmWb3r//d3f+nhV/eW2/\nr7/nrHvfv/763Dd48NKHqNtUj5lREa/gz+vm8GbdmweVSeRQlZx2KtWPPEzZhRdCczPpbdtIbdlC\nautW0lu2kN6+HSsvo9N111L98ENaNkfaBc1C3N71nwCT74FZNwW3D2OYUsFprgeLwJRHocvA/T7d\niospGTeOknHjSL79NomVq0gsW0bilaUk33oLNwsmLnb/YAItM8AwMzyZJNq9O/Exp1B0wgnEBgwg\nNmQwdhhT0YuISNuZO3cuzzzzDAsXLqS4uJitW7fScpCnl7y08UU2Nmxk89LN7CgrovfwQ1tv/M35\na5n1w+e49N6LqOpeCQTDiWORKL987XFuHn3rXgdEk8kksZh23SS7Yr16UXHtNZR//ipa5s4j9e67\neEMD1rmK+DEDiJ84Uvs60q7or2A+GHg2pFPw3DeCwipeGnai8DXVQSQGUx6BI48/6JfHjj6a2NFH\nUzp5EgCeTJJ6ex2pDRvwREswjCadxoqKIR4nWl1N9Jj+REr1sxcRyVcbN26ka9euFGcSdHWFAAAR\n5UlEQVSWDunatev7j82ZM4ebbrqJZDLJySefzIMPPkhxcTH9+vVj/vz5dO3alZdffpnPfvkyzvz6\nBBY9vQSLGK/OWs5ZNwbbknVL1vPSr1+mobaBCV8Zz3EThrSa4+0l63j2u3/konumckSvYPRQ3bt1\nPPPvz7J7x26KqoqY8NhETh16GldeeSXV1dUsWrSIUaNGcccdd3D99dezdOlSkskk06dPZ8qUKaxd\nu5bPfe5zNDQ0AHD//fdz2mmn5fLHKQUmUlpKycQJYccQ2S8VsPli8CeDobJ/+hqkk8GasR1V4/ag\niJ/yaDDMOgssFiN2TH9ix/TPyvuJiEj7c9ZZZ3HHHXcwaNAgJk2axMUXX8y4ceNoamriyiuvZM6c\nOQwaNIgrrriCBx98kBtvvPFDr9/QsIFEKkH3Xkdy4pSRFJXGGXPpKQAseeYVdtXu4ooHLmfrW7U8\ncfOTrRawqZYUT9z8FJf/+FK69v1gvoQ//WA2w88eyohzhjP3v17kuq9ez4I5CwBYtWoVzz33HNFo\nlFtuuYWJEyfyyCOPsGPHDkaPHs2kSZM48sgjmT17NiUlJaxevZpLL72U+fPn5/CnKSISDo0HyCcD\nJsMFM4OJihq3dbx1Yj0NDVugsjdM/U3WilcREekYOnXqxIIFC3jooYfo1q0bF198MTNmzGDlypX0\n79+fQYMGATBt2jReeOGFvV6/fud6wPc518GgMwZiEaNb/640bNvd6nMisQi9h/Vi8e9f+dD97yzb\nwNDJwYiiE885gRULVrz/2NSpU4lGowDMmjWLu+66i5EjRzJ+/Hiampp4++23SSQSfPGLX2T48OFM\nnTqV5cuXH/TPR0QkH6gHNt/0PBEufgp+dzVsfyOYrfgg1p3LW+lksFRO/4kw+W4oKg87kYiI5KFo\nNMr48eMZP348w4cPZ+bMmYwcOXKfz4/FYqQzE/2tqX3jYyfqi8aj7193Wj/IbBHjgjum8It/+hV/\n/4+5nH7FqXt/ZiSG4zQkguHA5eUfbPPcnSeffJLBgwd/6DXTp0+ne/fuLFmyhHQ6TUmJ5swQkcLU\nASqfAlTREz7zSzjmTGjYWvgzFDfvDIYN13wFzvmxilcRETkkK1euZPXq1e/fXrx4MX379mXIkCGs\nXbuW119/HYDHHnuMcePGAdCvXz8WLAiG8v7tD38jktl1Ki4romX3oa0vHi+Jc9Hdn2HZ7GUsfmYJ\nAL2H9WL5c0Gv6/LZKzhqWE82796012s/8YlPcN99972/HNyiRYsAqKuro2fPnkQiER577DFSqdQh\nZRMRae9UwOarojI4+4dw5ncg1QK7txbekOJ0KhgyXFwRTNY05nrQLHgiInKIdu3axbRp0zj++OMZ\nMWIEy5cvZ/r06ZSUlPDoo48ydepUhg8fTiQS4eqrrwbg9ttv54YbbuCMM87AzYPZ6oFjTzuWlS+s\n5uGrHuXtJesOOktpZSmXfO8i/v4fc1n1t9VMvmESrzy7lJ9Ne4Slf1rGuGvH0pzau0C+7bbbSCQS\njBgxgmHDhnHbbbcBcM011zBz5kzGjBnDqlWrPtRrKyJSSMzzoOipqalxTUTwMXa+C3/+JqybGxR7\nhTBLcXN90LM89GI47WvqdRURkdB9/YWbMIIhvrm2K9HAdSOvZ+ARg3L+WSIi7YGZLXD3mv09T91Z\nhaCiB5z/s6A3Fg96LVs5apsXEruD/KXVwSzD47+l4lVERNqFyqIqkulkm3yWu1Me1/ZPROSjVMAW\nCjM47gK4Yg6cfM0HhWAbbWgPW7I5yBuJwYQ74LJnoffosFOJiIi8b0DnAbS0wQHitKdxnO5lPXL+\nWSIi+UYFbKEp7gSjr4ErnoPhnw2G4jZsgURj2Mn25h5M0NSwNegxPvX/BrmP/zRENUG2iIi0L8dU\nHdMmn9Ocauao8p5EI9H9P1lEpINRlVCoyrrA2FvgpC/Ba/8NS2YEy9BYBEqqwl16J52Exh3BRBhV\nfWHUP8KAs4KJqURERNqp47ocD2akPU0kh9vRRCrBqT1Py9n7i4jkMxWwha68K5z0BRg5Ddb9DyyZ\nCe+8RND57lBcCdF47nMkm4Pe4EgUMBh4Doy4HLoPD4Y/i4iItHOVRZWM7HYiizYvpKq4MiefkUqn\nMItQ00On0YiItCaUAtbMzgZ+BESBh939rjBydCjRGPQbG1ya6mH9PHjjT7D2eWjZFfSKRqIQK4VY\nyeEVlZ4Ohiwn3xu2HAmGNg+9OFi79qiTIFacjVaJiIi0qcl9z2Lx5oWk0qmcDPHdldjFGb3G0qmo\nU9bfW0SkELR5AWtmUeABYDKwHnjZzJ529+VtnaXDKqmEY88KLqkEbFkOW1fChpfh3SWwcwNE4sEQ\n33QK0olgciXLbKjNMmvOevC4p4PHI9Hgfk9B5/7QcxT0PBG6HgddBqqnVURE8t7RlUczoc+Z/Hnd\nnKz3wjYmG6koquT8Yz+V1fcVESkkYfTAjgZed/c1AGb2K2AKoAI2DNE49DghuAy7KLivpQHq1kHj\nNmisDSaBql8fDAFOJYICNRILLmVdoKIXlHcLlr4p6wqVfSBWFG67REREcuS8Af/A0q2vUNu0lYqi\niqy8ZyKdoCWV4MsjvkJprADWcxcRyZEwCthewLo9bq8HTgkhh+xLUTl0GxJ2ChERkXapKFrE9aNu\n4J6X72Zny87DLmITqQQNyd1cNuRyBldr+ysi8nHCmIq2tXGkvteTzL5kZvPNbP6WLVvaIJaIiIjI\ngakuqeafT/4GXUq6UtdcRyqdOuj3cHd2tuykKdXM546fxmm9Ts9BUhGRwhJGAbse6LPH7d7Aho8+\nyd0fcvcad6/p1q1bm4UTERERORDVJdXcfMqtTOwziYZEA3XN9aR8/4Wsu7M7sZu6lnqO6tSLW065\nlTE9x7RBYhGR/BfGEOKXgYFm1h94B7gE+GwIOUREREQOS1G0iE8PupCTe57Mc2/NZtHmhTiOu1Mc\nLSYaiWIEa8e2pFpIeooIRvfyHpzV9xOc1L0mJ7MZi4gUqjYvYN09aWbXAX8iWEbnEXdf1tY5RERE\nRLKlT8XRXDXsC3ym5SJe27aCNTvWsKbuDepb6nGcsmgJg44YzLGdBzKg87H0qeiDaXZ+EZGDZu57\nnX7a7tTU1Pj8+fPDjiEiIiIiIiI5YGYL3L1mv8/LhwLWzLYAb4Wd4wB1BbaGHSLHOkIbQe0sNB2h\nnR2hjaB2FpKO0EZQOwtNR2hnR2gjqJ3tTV933+/kR3lRwOYTM5t/IEcO8llHaCOonYWmI7SzI7QR\n1M5C0hHaCGpnoekI7ewIbQS1M1+FMQuxiIiIiIiIyEFTASsiIiIiIiJ5QQVs9j0UdoA20BHaCGpn\noekI7ewIbQS1s5B0hDaC2lloOkI7O0IbQe3MSzoHVkRERERERPKCemBFREREREQkL6iAzRIze8TM\nNpvZq2FnyRUz62NmfzGzFWa2zMxuCDtTLphZiZm9ZGZLMu3817Az5YqZRc1skZk9E3aWXDGztWa2\n1MwWm1nBLihtZp3N7Akzey3zf/TUsDNlm5kNznyP713qzezGsHNlm5n9U+Zvz6tm9kszKwk7Uy6Y\n2Q2ZNi4rpO+xtf0BM6s2s9lmtjrz7xFhZsyGfbRzaub7TJtZ3s94uo823pP5O/uKmf3WzDqHmTEb\n9tHOb2fauNjMZpnZUWFmzIaP21c3s5vMzM2saxjZsmkf3+d0M3tnj+3nuWFmPFwqYLNnBnB22CFy\nLAl8zd2PA8YA15rZ8SFnyoVmYKK7nwCMBM42szEhZ8qVG4AVYYdoAxPcfWQhTSHfih8Bf3T3IcAJ\nFOD36u4rM9/jSOAkYDfw25BjZZWZ9QK+CtS4+zAgClwSbqrsM7NhwBeB0QS/r+eZ2cBwU2XNDPbe\nH/gXYI67DwTmZG7nuxns3c5XgU8DL7R5mtyYwd5tnA0Mc/cRwCrg5rYOlQMz2Lud97j7iMzf22eA\nb7V5quybQSv76mbWB5gMvN3WgXJkBq3XJD94bxvq7n9o40xZpQI2S9z9BWBb2Dlyyd03uvvCzPWd\nBDvIvcJNlX0e2JW5Gc9cCu5kcTPrDXwSeDjsLHJ4zKwSGAv8HMDdW9x9R7ipcu5M4A13fyvsIDkQ\nA0rNLAaUARtCzpMLxwHz3H23uyeBvwIXhJwpK/axPzAFmJm5PhP4VJuGyoHW2unuK9x9ZUiRsm4f\nbZyV+Z0FmAf0bvNgWbaPdtbvcbOcAtgP+ph99R8AX6cA2ggdoyZRASuHxMz6AScCL4abJDcyQ2sX\nA5uB2e5eiO38IcEf7HTYQXLMgVlmtsDMvhR2mBw5BtgCPJoZEv6wmZWHHSrHLgF+GXaIbHP3d4Dv\nEfQEbATq3H1WuKly4lVgrJl1MbMy4FygT8iZcqm7u2+E4GAwcGTIeSQ7Pg88G3aIXDGzO81sHXAZ\nhdEDuxczOx94x92XhJ2lDVyXGRb+SL6fxqACVg6amXUCngRu/MgRuoLh7qnMsJnewOjMcLeCYWbn\nAZvdfUHYWdrA6e4+CjiHYNj72LAD5UAMGAU86O4nAg0UxhDFVplZEXA+8Juws2RbZqdiCtAfOAoo\nN7PLw02Vfe6+AribYDjmH4ElBKepiOQFM7uV4Hf28bCz5Iq73+rufQjaeF3YebItc/DsVgq0OP+I\nB4EBBKfGbQS+H26cw6MCVg6KmcUJitfH3f2psPPkWmYY5vMU3vnNpwPnm9la4FfARDP7z3Aj5Ya7\nb8j8u5ngfMnR4SbKifXA+j1GCjxBUNAWqnOAhe6+KewgOTAJeNPdt7h7AngKOC3kTDnh7j9391Hu\nPpZguNvqsDPl0CYz6wmQ+XdzyHnkMJjZNOA84DLvGOtR/gK4MOwQOTCA4GDhksz+UG9goZn1CDVV\nDrj7pkznTBr4GXm+L6QCVg6YmRnBOXYr3P3esPPkipl1e29WQTMrJdihfC3cVNnl7je7e29370cw\nFPPP7l5wvTxmVm5mFe9dB84iGLpYUNz9XWCdmQ3O3HUmsDzESLl2KQU4fDjjbWCMmZVl/uaeSQFO\nyAVgZkdm/j2aYOKfQv1OAZ4GpmWuTwP+O8QschjM7GzgG8D57r477Dy58pFJ1c6nwPaDANx9qbsf\n6e79MvtD64FRmW1qQXnvAFrGBeT5vlAs7ACFwsx+CYwHuprZeuB2d/95uKmy7nTgc8DSzPmhALfk\n+0xmregJzDSzKMFBnl+7e8EuM1PgugO/DeoAYsAv3P2P4UbKmeuBxzPDa9cAV4WcJycyQ74mA18O\nO0suuPuLZvYEsJBgeOIi4KFwU+XMk2bWBUgA17r79rADZUNr+wPAXcCvzewLBAcppoaXMDv20c5t\nwH1AN+D3ZrbY3T8RXsrDs4823gwUA7Mz25Z57n51aCGzYB/tPDdzUDQNvAXkdRuhw+yr7+v7HG9m\nIwnmBVlLnm9DrWOMfBAREREREZF8pyHEIiIiIiIikhdUwIqIiIiIiEheUAErIiIiIiIieUEFrIiI\niIiIiOQFFbAiIiIiIiKSF1TAioiItAEzS5nZYjN71cx+995604f4Xs+bWU0284mIiOQDFbAiIiJt\no9HdR7r7MIL1Mq8NO5CIiEi+UQErIiLS9uYCvQDMrJOZzTGzhWa21MymZO7vZ2YrzOxnZrbMzGaZ\nWemeb2JmETObaWbfCaENIiIibU4FrIiISBsysyhwJvB05q4m4AJ3HwVMAL5vZpZ5bCDwgLsPBXYA\nF+7xVjHgcWCVu3+zTcKLiIiETAWsiIhI2yg1s8VALVANzM7cb8C/mdkrwHMEPbPdM4+96e6LM9cX\nAP32eL+fAq+6+525Di4iItJeqIAVERFpG43uPhLoCxTxwTmwlwHdgJMyj28CSjKPNe/x+hRBr+t7\n/geYYGYliIiIdBAqYEVERNqQu9cBXwVuMrM4UAVsdveEmU0gKHAPxM+BPwC/MbPY/p4sIiJSCFTA\nioiItDF3XwQsAS4hOI+1xszmE/TGvnYQ73MvsBB4zMy0TRcRkYJn7h52BhEREREREZH90tFaERER\nERERyQsqYEVERERERCQvqIAVERERERGRvKACVkRERERERPKCClgRERERERHJCypgRUREREREJC+o\ngBUREREREZG8oAJWRERERERE8sL/AnGhRi0jmKeKAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f6035133a20>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "def plot_optional():\n",
    "    import matplotlib as plt\n",
    "    %matplotlib inline\n",
    "    Top15 = answer_one()\n",
    "    ax = Top15.plot(x='Rank', y='% Renewable', kind='scatter', \n",
    "                    c=['#e41a1c','#377eb8','#e41a1c','#4daf4a','#4daf4a','#377eb8','#4daf4a','#e41a1c',\n",
    "                       '#4daf4a','#e41a1c','#4daf4a','#4daf4a','#e41a1c','#dede00','#ff7f00'], \n",
    "                    xticks=range(1,16), s=6*Top15['2014']/10**10, alpha=.75, figsize=[16,6]);\n",
    "\n",
    "    for i, txt in enumerate(Top15.index):\n",
    "        ax.annotate(txt, [Top15['Rank'][i], Top15['% Renewable'][i]], ha='center')\n",
    "\n",
    "    print(\"This is an example of a visualization that can be created to help understand the data. \\\n",
    "This is a bubble chart showing % Renewable vs. Rank. The size of the bubble corresponds to the countries' \\\n",
    "2014 GDP, and the color corresponds to the continent.\")\n",
    "    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "#plot_optional() # Be sure to comment out plot_optional() before submitting the assignment!"
   ]
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "coursera": {
   "course_slug": "python-data-analysis",
   "graded_item_id": "zAr06",
   "launcher_item_id": "KSSjT",
   "part_id": "SL3fU"
  },
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.2"
  },
  "umich": {
   "id": "Assignment 3",
   "version": "1.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}