diff --git a/ENTREGA.ipynb b/ENTREGA.ipynb new file mode 100644 index 0000000000000000000000000000000000000000..aa9f051f3ca12095332fd9fd182197dd72adf227 --- /dev/null +++ b/ENTREGA.ipynb @@ -0,0 +1,201 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Nombre: Jennifer Ortega \n", + "### Ejercicio No. 1 - Diagrama de Hertzsprung-Russell\n", + "\n", + "\n", + "El diagrama de Hertzsprung-Russell es un gráfico de disperción que relaciona las temperaturas efectivas de las estrellas con su luminosidad. El gráfico tiene 5 regiones principales, cada una para un grupo diferente de estrellas.\n", + "\n", + "- **Grupo 1**: Contiene estrellas promedio, donde la luminocidad esta directamente relacionada con la Temperatura, son estrellas que estan al principio o en la mitad de sus vidas.\n", + "- **Grupo 2 y 3**: Gigantes rojas y las supergigantes, son estrellas viejas de tamaño pequeño y mediano que se han hinchado, por lo tanto son muy luminosas aunque no esten calientes.\n", + "- **Grupo 4**: Enanas blancas, son super calientes pero no son muy luminosas, porque son muy pequeñas.\n", + "- **Grupo 5**: Gigantes azules son enormes y no tienen un tiempo de vida largo. \n" + ] + }, + { + "cell_type": "code", + "execution_count": 1, + "metadata": {}, + "outputs": [], + "source": [ + "import numpy as np\n", + "import matplotlib.pyplot as plt\n", + "import math \n", + "\n", + "%matplotlib inline" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "lum temp radius\n", + "304.2285727480961 3654.60109945958\n" + ] + } + ], + "source": [ + "f = open(\"data/giants.txt\", \"r\")\n", + "print(f.read(50))" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": {}, + "outputs": [], + "source": [ + "filestxt = [\"data/giants.txt\",\"data/supergiants.txt\"]\n", + "filescsv = [\"data/dwarfs.csv\", \"data/ms.csv\"]\n" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "metadata": {}, + "outputs": [], + "source": [ + "giants = np.loadtxt('data/giants.txt', delimiter=' ',skiprows=1)\n", + "supergiants = np.loadtxt('data/supergiants.txt', delimiter=' ',skiprows=1)\n", + "dwarfs = np.loadtxt('data/dwarfs.csv', delimiter=',',skiprows=1)\n", + "ms = np.loadtxt('data/ms.csv', delimiter=',',skiprows=1)\n", + "data_final = np.concatenate((giants,supergiants,dwarfs,ms))\n", + "lum = data_final[0:, 0]\n", + "temp = data_final[0:, 1]\n", + "rad = data_final[0:, 2]\n", + "lumlog = math.log10(3)\n", + " \n" + ] + }, + { + "cell_type": "code", + "execution_count": 20, + "metadata": {}, + "outputs": [], + "source": [ + "lum_giants = giants[0:, 0]\n", + "lum_supergiants = supergiants[0:, 0]\n", + "lum_dwarf = dwarfs[0:, 0]\n", + "lum_ms = ms[0:, 0]\n", + "\n", + "temp_giants = giants[0:, 1]\n", + "temp_supergiants = supergiants[0:, 1]\n", + "temp_dwarf = dwarfs[0:, 1]\n", + "temp_ms = ms[0:, 1]\n", + "\n", + "rad_giants = giants[0:, 2]\n", + "rad_supergiants = supergiants[0:, 2]\n", + "rad_dwarf = dwarfs[0:, 2]\n", + "rad_ms = ms[0:, 2]\n" + ] + }, + { + "cell_type": "code", + "execution_count": 34, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "['Solarize_Light2', '_classic_test_patch', 'bmh', 'classic', 'dark_background', 'fast', 'fivethirtyeight', 'ggplot', 'grayscale', 'seaborn', 'seaborn-bright', 'seaborn-colorblind', 'seaborn-dark', 'seaborn-dark-palette', 'seaborn-darkgrid', 'seaborn-deep', 'seaborn-muted', 'seaborn-notebook', 'seaborn-paper', 'seaborn-pastel', 'seaborn-poster', 'seaborn-talk', 'seaborn-ticks', 'seaborn-white', 'seaborn-whitegrid', 'tableau-colorblind10']\n" + ] + } + ], + "source": [ + "print(plt.style.available)\n", + "\n", + "temp_color=(-temp+max(temp)) # para cambiar la escala y que los colores coincidan con los pedidos en la tarea" + ] + }, + { + "cell_type": "code", + "execution_count": 43, + "metadata": { + "scrolled": true + }, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAyEAAAHFCAYAAADlpPZOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAxOAAAMTgF/d4wjAACQD0lEQVR4nOzdd3hTZR/G8W+6d9ktW6CyZAkyRGWqgKKIC1ARFATc4AIZigNEXIACIoKA8qKICLJEEVDZioLIpuxRZmlLV9ok7x+xhdB0pM3ouD/X1YvmjOf8QlfunGcY1qxZY0FERERERMRNvDxdgIiIiIiIlCwKISIiIiIi4lYKISIiIiIi4lYKISIiIiIi4lYKISIiIiIi4lYKISIiIiIi4lYKISIiIiIi4lYKIVIi/Pjjj7Rv3z7zo0uXLvTs2ZNRo0axZs0aLBbb5XJiYmJo3749P/74o4cqLtwy/n+WLVtmd/+zzz7L4MGDnXrNdevWMX/+fKe2Wdxd+T3foUMHunXrxogRIzh06JCnS8uRvZ+/cePG0bNnTw9WJSIizuTj6QJE3Gn06NGUL18eo9HImTNn2LRpE2+//TZLly5l7Nix+Pv7A1CmTBkmT55MpUqVPFyxZFi3bh1bt27lwQcf9HQpRUrnzp256667MJlMHDx4kC+++IKhQ4cyc+ZMQkJCPF2eiIiUUAohUqJERUVRuXLlzMe33347v/76K2+88QbTpk3jueeeA8DPz4/69eu7vT6j0Yifn5/br1uYFfb/k8JeX7ly5TK/lxs2bEhQUBBjx45ly5YtdOjQwcPViYhISaUQIiVe27Ztuemmm1i6dCkDBgwgICCAmJgYevXqxdChQ+ncuTMAe/bsYd68eezatYv4+HgiIiK45ZZbePTRRzPvoACYTCZmzZrF8uXLSUxMpH79+gwePJg+ffrQp08f+vbtC8CsWbOYPXs2M2fOZMqUKfz77780bdqUMWPG8Mcff/Ddd9+xf/9+EhMTqVixIl26dOG+++7D29s781o9e/akYcOGNG/enK+++orTp09Tp04dXnnlFcqWLcuUKVP47bff8Pb25tZbb2XgwIGZ5xuNRqZPn86ff/5JTEwMgYGB1K1bl0GDBlGtWjWn/z+npKQwe/Zs1q5dy7lz5yhXrhx33nknDz30EF5e1p6h27ZtY8iQIbzxxhts2bKFdevWkZ6ezs0338zKlSsBaxcjgIiICL7++uvMx1fL2A+watUqvvnmG44fP46XlxcRERHcc8893H333YC1q8/WrVt5/fXX+eSTTzh48CBlypThwQcf5N57781s88cff+Tdd99lwoQJfP/99/z5559ERkby+eefZ3Y/mzBhgk0dPXv2pEmTJgwbNsymjcmTJ7Nw4UI2btxIYGAgbdq0YdCgQTaB5uTJk0yaNIlt27YRGBjIbbfdRtWqVfnwww+ZN28ekZGRDn8dateuDcDp06ezrTFD+/btbb5njx07xmeffca///5LYmIipUuXpl69eowaNQpvb2+Sk5OZPn0669evJzY2luDgYGrWrMnzzz+f+T1lMpn4+uuvWblyJTExMYSFhdGhQwf69+9fqMOciIg4l0KICNCyZUvWrVvH3r17ady4sd1jTp8+TVRUFJ07dyYwMJDDhw8zZ84cTp06xWuvvZZ53KxZs5g7dy49evSgWbNm7Nu3jxEjRmR77ZEjR9KlSxd69eqV+WL81KlTNG3alO7du+Pn58fevXuZNWsWFy9eZMCAATbn//PPP5w8eZIBAwaQlpbG5MmTef3116lYsSKVK1dm1KhR/PPPP3z55ZdUqlSJe+65B7CGkKSkJHr37k2ZMmVISEhg8eLFPP3008yePZsyZcrk+v9mNpsxmUy5HmcymXjllVc4cuQIvXv3pkaNGuzevZs5c+YQHx/PU089ZXP8xx9/TIsWLRg+fDhGo5EaNWpw8eJF9u7dy5gxYwDw9fUFYPLkyTbnnj9/njFjxlC9enUAduzYwdixY7n33nsZNGgQZrOZo0ePkpiYaHNeUlISb775Jr169aJy5cqsXr2ajz/+mKCgoMwgmmHMmDF06NCBN954I0/P356xY8fSoUMH3nzzTXbu3Mns2bMJDQ3lscceAyAtLY2XX34Zo9HI4MGDKVWqFMuXL+e3337L1/UyxMTEAOSrq+Grr75KaGgogwcPJjw8nHPnzrF582bMZjPe3t5MnjyZDRs20K9fP6pUqUJ8fDz//vsvly5dymxjzJgxbNy4kZ49e9KgQQOOHDnCzJkziYmJ4c033yzQcxMRkaJDIUQEqFChAmB9AZudtm3bZn5usVho2LAhwcHBvPPOOzz//POEh4eTkJDAggULuOuuuxg4cCAAN9xwAz4+PkydOtVuu927d+f++++32ZbxDn3GtRo1akR6ejrffPMN/fv3zwwrAMnJybz77ruZ/fsvXLjAJ598Qt26dXnyyScza9i0aRO//vprZggJCQnh5ZdfzmzHZDLRvHlz7r33Xn755RceeOCBXP/fPvzwQz788EO7+64Mc7/88gs7duxgwoQJmdubNWsGwOzZs+nVqxelS5fOPL5u3bo2tQGUKlUKHx+fLN3krnycmprKxx9/TLly5Rg+fDgAu3btIiQkhGeeeSbzuObNm2epNykpiRdffDGzi1KLFi04d+4cs2bNolOnThgMhsxj27Zty6BBg3L4n8ldx44dMwNHs2bN2L17N6tXr87c9uOPP3Ly5EmmTJlCvXr1AGtY7t+/v0PXsVgsmEymzDEh06ZNo379+tx0000OtRMXF8eJEyd4++23bc699dZbMz/fuXMnHTt25M4778zcdsstt2R+/s8//7BmzRqGDRtGp06dMp97aGgoY8eO5cCBA0RFRTlUl4iIFE0KISJXuPKF5tUSExP56quv+O233zhz5gzp6emZ+06cOEF4eDgHDx4kJSWFdu3a2Zzbtm3bbEPIlS/SMpw/f55Zs2bxxx9/cO7cOZt32y9evGhzl6J+/fo2A4wzur1c/UK7atWq7Nmzx2bbmjVrmD9/PseOHbO5M3Ds2LHs/hts9O7d2+6L2ffee8/m8ZYtW4iIiKBBgwY2z+WGG25gxowZ7Nq1y6adm2++OU/Xv5LFYuGdd97h5MmTfPLJJ4SHhwNQp04dEhISMu9eNGzY0O6AbC8vL9q0aWOzrX379rz//vucO3eO8uXLF6i+q7Vq1crmcc2aNdm6dWvm4127dhEREZEZQMD6/dmmTRsOHjyYuc1sNtvM7mYwGGxC6ty5c5k7d27m48jISD788EN8fBz79R8WFkalSpWYPn06sbGxNGnShCpVqtgcU7duXVauXEl4eDjNmzcnKirKpvvgli1b8PX1pW3btjbfBxnfq9u3b1cIEREpIRRCRIAzZ84AULZs2WyPeffdd/nrr7/o27cvUVFRBAYGsnv3biZOnIjRaAQu30kpVaqUzblXvst/tauvaTabGTFiBOfPn6dPnz5Uq1YNf39/1q1bx1dffZV5rQyhoaE2jzO6KdnbfuW5GzZs4M0336RTp0706dOH8PBwvLy8GDZsWJZrZCciIoI6depk2R4YGGjz+OLFi5w+fdrmXfMrxcfH2zzO6euQnS+++IL169fz3nvvUbVq1cztTZo0YfTo0SxcuJBRo0YB1rs0Tz31FLVq1co8LjQ0NMsL84ywd/bsWZsQkp/6rhYWFmbz2NfXl7S0tMzHFy5cyPJ9BFm/l8aPH585XgagU6dONmM7unTpQrdu3TAajfz111/MmTOHN998kw8++CDH0H01g8HAe++9x6xZs5g+fTrx8fFUrFiRHj160K1bNwCee+45ypQpw4oVK5gxYwZhYWHcfvvt9OvXj4CAAC5evEhaWhpdunSxe42rvw9ERKT4UggRATZt2oSfn1/moN2rGY1G1q9fT9++fW26Tl35jjRcfnF68eJFm+2xsbHZXvvqF4InT55k7969DB8+nNtuuy1z+4YNG/L0XPJq9erVVK5c2eYFa3p6ukteCIaFhVGxYkVef/11u/uvHmDtyItjsA48/+qrr3jllVdo0qRJlv1t27albdu2JCcns23bNqZNm8bQoUOZP39+5l2DhIQE0tPTbYLIhQsXAGwCSHb1+fn5kZSUlGV7QkKCQ88lQ5kyZThy5EiW7Vd/L/Xt25fu3btnPs64A5ShbNmymUGxYcOGWCwWZs+eza+//pp5x87Pz88mAIG1+9XVKlWqxPDhw7FYLERHR/P9998zYcIEIiMjadmyJYGBgTzxxBM88cQTxMTE8Ntvv/HZZ5/h4+PDwIEDCQsLw8/Pj0mTJtl9zs4IdyIiUjRosUIp8X799Vc2bNjA3XffTUBAgN1jjEYjZrM5yzvlVy9mWLNmTQICAli7dm2Wa+RVSkoKgE03lvT0dFatWpXnNvIiNTXV5hoAP/30E2az2anXAev4ijNnzhAYGEidOnWyfFz9wtmeq+/kZNi5cyfjx4+nV69eWQaQXy0wMJAbb7yRu+66i/Pnz9sELrPZnGXQ95o1a4iIiKBcuXK51hcREcGxY8dsXsxv377dbjDJi/r163P69Gl2796duc1isWSpMTIy0ub/MrcZs3r16kW5cuWYM2dOZjeuiIgIDh8+bHPcpk2bsm3DYDAQFRWVOaGAvcUPIyMjefDBB6lZs2bm/hYtWmA0Grl06ZLd74O8/D+LiEjxoDshUqIcOHCAuLg40tLSOHPmDBs3buTXX3+lWbNmPPHEE9meFxISQv369Zk/fz5lypQhPDycFStWcO7cOZvjQkNDuf/++5k7dy5BQUE0a9aM/fv3s3z5ciBv7/BXr16diIgIZsyYgbe3N97e3ixYsKBgT9yO5s2bs27dOiZPnkyrVq3Yt28fCxcudMkCdrfeeisrVqzgxRdf5IEHHqBWrVqkp6dz8uRJNmzYwFtvvZVtAMxQvXp14uPjWbx4MXXq1MHPz4+IiAhGjRpFtWrVaN26Nbt27co83tfXl2uvvZaZM2cSGxvL9ddfT9myZTl79iwLFy4kKirKprtTUFAQ06ZNIy4uLnN2rK1btzJ06NA8fd06dOjA0qVLGT9+PJ07d+bUqVN8++23BAcH5+v/rHPnznz99de89tpr9OvXL3N2rIyZphy9W5TB39+fhx56iEmTJvHbb7/Rtm1bOnTowPjx4zO/F6Kjo226eAFER0fzySef0L59eypXrozJZGLlypV4e3tz/fXXA/D000/TunVratasSWBgINu3byc6Oprbb78dsHaN69ChA6NHj+aBBx6gbt26eHl5ERMTw+bNmxkwYIBNVzoRESm+FEKkRBk9ejRg7X5SunRprr32WkaNGkXbtm1zfVE3atQoPvroIyZOnIi/vz/t2rXjmWeeyZyFKUPfvn2xWCwsX76chQsXUq9ePYYOHcqzzz6bpxekvr6+vP3220ycOJF33nmH0NBQunTpQkREBO+//36+n/vVunbtytmzZ1mxYgVLliyhTp06jB07NnPchDP5+Pjw3nvv8b///Y+lS5cSExNDQEAAlSpVolWrVpnjWHJy5513snv3bj7//HMuXbpEREQEEyZMIDY2ltjYWJvZr+DyOiH16tVj4cKFTJ48mYSEBEqVKsUNN9zA448/bnN8UFAQo0aNslkn5Jlnnsn17kqG66+/niFDhjB//nx+++03rr32WoYPH55tF7Tc+Pr6Mn78eD7++GM++ugjAgMD6dixI/Xq1eOzzz7Ld7gB69f+m2++4auvvqJNmzZ06tSJM2fOZH4vNGzYkDfffJNHHnkk85wyZcpQoUIFvv32W86ePYufnx81atRg7Nixmd29GjVqxNq1a5k3bx4mk4mKFSvy1FNPcd9992W2M2LECBYuXMiKFSv46quvMsNk8+bNcxw7JSIixYthzZo1ltwP84zVq1ezaNEioqOjSUpKYtWqVTbdR9q3b4+fn5/NTDCTJ0+mZs2anihXJFu//voro0ePZuLEiTRq1MjT5chVMhYr/Pbbbz1dSq5effVVjh49ajPjlYiISFFTqO+EhISE0K1bN1JTU7NM+Zlh7NixmesNiBQGu3btYtOmTdSrVw8/Pz/27dvHvHnzqF+/Pg0bNvR0eVKEzJ8/n8DAQKpUqUJSUhK//vormzZtYsiQIZ4uTUREpEAKdQhp0aIFANu2bfNsISIOCAwM5J9//mHx4sUkJiZSunRp2rVrxxNPPJHvfvxSMvn6+rJgwQLOnDmDyWSiatWqvPTSSzaLAYqIiBRFhbo7VoZt27YxZMgQu92xSpcujclkIiIigrvvvpuuXbt6sFIREREREclNob4Tkpv333+fBg0a4OXlxdatWxkzZgwmkylz4awrmc1mzp8/T2BgoN6NFhERKSIsFgvJycmULVvWZgyoiBRtRTqEXDkWpFWrVtx33338/PPPdkPI+fPnefDBB91ZnoiIiDjJ/PnzsywcKiJFV5EOIVczGAyZi29dLTAwEIBjx44RFhbmzrJcavjw4YwdO9bTZUgu9HUq/PQ1Kvwc+hqlpEDf+2DTOtcWlZ0Ot8O0/0Eepp8u6tISE9n7zTf8O2MGc/79lzuc3H5Io0Y8/88/mX/HRaR4KNQhxGQyYTKZMlcgNhqNeHt74+Pjw4EDB7BYLNSsWRMvLy/++usvvvvuO/r06WO3rYwuWGFhYcUqhPj5+RWr51Nc6etU+OlrVPg59DV693XYsg481Xtn7U/w2UR49U0PFWBlMZm4tGcPF7duJW7rVuK2bSPtwgXMKSkAeAUE4BMeTnjjxoQ3a0apG24gpH59vHxyf3mQnprK5rfe4u+PP8YYHw+AN5DzsqOOu/TPP0D+F+gUkcKpUIeQn3/+mXfffTfz8R13WN9f+eijj0hKSmLatGmcOXMGb29vIiIi6N+/P3fffbenyhURkcJg0zqYPsnTVcDEd+DO7tDoerdf+tyaNRyeOpUzy5ZhSkrK9fjY9eszP/cKCKD87bdzzZNPUr5TJ7sv/k9t2cLKvn25sHu3U+sWkZKjUIeQzp0757hacevWrd1YTeHUqVMnT5cgeaCvU+GX5WuUehziV0PiVkg9DKSDdzgENYbg5hDWFgze9poSF8nTz1FyMjz7GJjNri8oN+np8GxfWPWnW7plpSckcGz2bA5PmcKlAoQDc0oKp3/4gdM//EBwVBTVBw2i6uOP41e6NOmpqWx47TW2fvABFpMpy7l1CvIERKREKRJT9DpDYmIiXbt2JS4uTl0uRCR7Cevg5Hi4uBzI+iIrk18VqDAQIoeAd7DbypNcjHsd3vdsF6gs3ngfnn7RpZc4tWgROwYNIvX0aZe071umDNVfeYXNs2e7/e5HCjACWLp0KcHB+lkTKS40152ICIApEQ49DbvawMUl5BhAAIzH4fgo2NEQ4te6o0LJTVoazJnm6SqymjnZZXdmjOfPs7VXL/7s3t1lAQQg7cIFDgwbhv/u3ej+n4g4g0KIiEjaWdh1C5yZAjh4czj1EOzuCGemu6Q0ccCS7+CM616I59uRQ7BqhdObPbNyJWuvu46TX3/t9LazEwJUB3Q/QkQKSiFEREo2UxLs6QRJfxegETMcGgjn/ue0siQfZk7xdAXZc3JtJ+bNY0vXri69+5Edb6AioI7NIlIQCiEiUrIdG1rAAJLBAocGQeoRJ7QlDjt3Fjb97ukqsrf6R0hMdEpTx778kr8eeQRLerpT2ssPAxABhHusAhEp6hRCRBxgxMIvXGIusfyPi6zmEmmOdt+RwiNhA5ye7Lz2zAnWOyLiftu3erqCnJnN8O+2AjdzatEitj9WSGb/AsoDoZ4uQkSKpEI9Ra9IYXGGdGYTyzdc5NxVA5bL4U1PSvEopamgH6mi5dT7ODwGJDdxKyFpBwQ1dG67krPCHkIAtv0JLW/K9+mJ0dH8/fDDdqfG9ZSMOyKpgNHDtYhI0aI7ISK52EEKd3CIyZzPEkAAzmHiE85zJ4f4lxQPVCj5YjwFsT+4pu3Tn7qmXcnetj89XUHuChCULBYL2/v1y9PCg+6WEURERByhECKSg2hS6c1Ru+Hjamcx8SjHOKT3A4uGhF/JdRre/Ipf7Zp2JXtHD3m6gtwVoMbDn3zC+V9/dWIxzhUAlPF0ESJSpCiEiORgNGe4SN77Xl/AxBsUwilCJatEF3bfSdkHpkuua1+ySk72dAW5y2eNiQcPsvvVV51cjPOVBvw8XYSIFBkKISLZOISRdTg+m81vJHJUd0MKP5fOYmWG1KMubF+yKCQDtXNkzt+dt/1vv43JSTNruZIXUNbTRYhIkaEQIpKNucTma8iyGZjLRSdXI05ncfXgXs9Nn1oiBQR4uoLcBQQ6fIoxNpYTblyMsKCCQSuqi0ieKISIZGMPqR45V9zEp5Rr2/d2cftiq2x5T1eQu3zUeOyLLzAXha5m/zGgtUNEJG8UQkSykVKAqVuTHRhHIh4S1Nh1bfuUAf9qrmtfsmrU1NMV5M7BGi0WC0c+LXozrSmEiEheKISIZCOsAD8eYeqQUPgFt3Bh281d17bY17iZpyvInYM1XtyyhcT9+11UjOv4AEGeLkJECj2FEJFs3EJwvs9tW4BzxU1CWoJ/lGvaLveIa9qV7BXDEBK7ebOLCnG9IjBCR0Q8TCFEJBv3EU4QBofPC8GL7oS5oCJxKoMBIp5yfrs+5aHMA85vV3JW61qoWt3TVWSvTn2IrOjQKXFbi8Aq8Nnw93QBIlLoKYSIZCMMbx7IR+/mBwgnRN2xioaIpyCgnnPbrPYBeOklmNsZDPDoAE9Xkb2+gxw+RSFERIozhRCRHLxKBZqR92k1WxDIMIrALD1i5eUPtWaBwcc57ZW6G8r3dk5b4rhH+oNfIVwuLzgYevZx6JT0pCQu7dnjooJczxdN1SsiOVMIEclBAF7MoSq3EpLrsbcRwiyq4q8fq6IlpAXUnE2BXzIFt4Sor5xSkuRT+QrQ9T5PV5HVfQ9DqGNdNFNPnsRicvVaNq7lpGgvIsWUXi2J5CIYLz6nCoupzv2E43/FOJEADDxAOD9QnelUIUg/UkVTuYfg2u+sU+vmR+l7oN7P4B3q1LIkH4a+AYGOLwroMqFh8NIoh08zpaS4oBj3cnxEnYiUJHqjQiSPGhNIYwJ5h0guYvpvUS5vfPWntngo0w1CdsLhZyB2IeRlnRjfCOsYkHIPu7w8yaNa18LwMTDqBU9XYvXWh1CpisOnWdLSXFCMe+k3o4jkRG/bijjIFwPl8aEcPgogxY1fJNReAI33QcWXIOh6MPjaHuMbAaW6QtQ8aHJMAaQwGvg8tLzJ01VAh07wSL98nerlX/SHdud/uVcRKQl0J0RE5GoBUVDtPevn5lQwngBLOviEW0OIFG5eXvDJbOjcCs6f80wNkZVgwuf5Pt07qOgv92f2dAEiUqjpToiISE68/CGgJgTWVgApSmrUgm9+dHhAuFOUKQvf/ZyvblgZAqtWLdJBxAIU/Q5lIuJKCiEiIlI8NWkGC3+xhgJ3iagIi9daFycsAIO3N2GNGzunJg8wou5YIpIzhRARESm+rr8Blq2HRk1df63mN8Ly9VCvgVOaC2/WzCnteEKqpwsQkUJPIURERIq3a+vAT5vh1bdcs5hhQAC88T4sWwfVazit2VI33OC0ttxNIUREcqMQIiIixZ+PD7w4Elb9Ce1uA4MTZrbz8oLb7oQ12+DpF62PnajMzTc7tT13SvZ0ASJS6CmEiIhIyVG/ISz4CTbugYGDIbyU422ULQfPvgJ/HIB5S613WlwguFYtyrZv75K2XSkF3QkRkdxpil4RESl5omrDmI9g5Fj4YyNs33r54/RJyFixPCAAKlaBxs2gyQ3Wf29oBW5ax+Oap57i/Jo1brmWs8R5ugARKRIUQkREpOQKDIQ2HawfhVDkPfcQUKkSKSdPerqUPDEBCZ4uQkSKBHXHEhERKaS8fHyoPnCgp8vIs3g0Na+I5I1CiIiISCFW88UXCapZ09Nl5CoduODpIkSkyFAIERERKcR8goNpPHOmc2b0cqEzgNnTRYhIkaEQIiIiUsiVa9uWa55+2tNlZCseSPR0ESJSpCiEiIiIFAH1xo0j+NprPV1GFmnAWU8XISJFjkKIiIhIEeATHEzLFSvwr1jR06VkMgEnUTcsEXGcQoiIiEgREVyrFjf+/DN+FSp4uhRMwAnA6OlCRKRIUggREREpQkKvu46bfv+dwGrVPFZDOnAcrYwuIvmnECIiIlLEhNSuzc2bNlHhzjvdfu0k4Bi6AyIiBaMQIiIiUgQFVKxIy6VLaTJrFr6lSrn8emas0/CewHonRESkIBRCREREirCqffrQbudOIu+5x2VriSQCR4A4l7QuIiWRQoiIiEgRF1CpEs2//572e/ZQY/Bgp9wZ8QkNpfqgQdCpEyfR3Q8RcS6FEBERkWIipHZtGnz0EbeeOEHjzz8nols3AqpUyfP5/hUrEtG1Kw0nT+a2EydoNHUqdy5dSt2HHnJh1SJSEvl4ugARERFxLp+gIKr160e1fv0ASD1zhotbtxL/99+kxcZiSkkBwDsgAJ/wcMIaN6ZUs2YEVKqUpS0vHx+6fPkl/uHhbJ861a3PQ0SKL4UQERGRYs6/QgUiunQhokuXfJ1v8PKi45QpVGzdmjXPPUdqbKyTKxSRkkbdsURERCRP6j/yCH137aLW3Xd7uhQRKeIUQkRERCTPgiMj6bZ4MZ2//BL/0qU9XY6IFFEKISIiIuKw+o88wmN793LT228TWrWqp8sRkSJGIURERETyJah8eVqOGEG/Q4e4e9Eiqt9+u9PWKjF4eVHzrrvoumCBU9oTkcJFA9NFRESkQLy8vYnq1o2obt2IPXCAQ8uXc2brVk5v3cqFPXuwmEy5t+HjQ5n69Ylo1oyIZs2o2bUrYdWrEx8f74ZnICLuphAiIiIiTlM6KorSzz2X+TgtKYkzf//N2e3bMcbFkZacjNloxNvfH5+AAPxLl6Z8kyaUb9wY38BAD1YuIu6kECIiIiIu4xsUROWbbqLyTTd5uhQRKUQKdQhZvXo1ixYtIjo6mqSkJFatWoW3t3fm/ujoaCZNmsS+ffsIDg6ma9eu9OnTB4OT+qOKiIiIiIjzFeqB6SEhIXTr1o2nn346y76kpCReeeUVGjRowKJFixg/fjzLli1jgQawiYiIiIgUaoU6hLRo0YKOHTtSqVKlLPt+++03zGYzjz/+OP7+/tSsWZMePXqwaNEi9xcqIiIiIiJ5VqhDSE6io6OJioqy6Z5Vt25dTp48SWJiogcrExERERGRnBTZEJKYmEhISIjNttDQUMDaVUtERERERAqnQj0wPSfBwcGcPXvWZltCQgIAQUFB2Z43fPhw/Pz8AOjUqROdOnVyXZEiIiLisJUrV7Jy5UoAjEajh6sREVcosiGkVq1arFq1CpPJlNkla+/evVSqVIng4OBszxs7dixhYWHuKlNEREQcdOWbhPHx8UyePNnDFYmIsxXq7lgmkwmj0UhaWhpgfTfEaDRiNptp06YNXl5efPHFF6SmpnLo0CHmz59Pt27dPFy1iIiIiIjkpFDfCfn555959913Mx/fcccdAHz00Uc0adKE8ePHM3HiRLp160ZQUBB33303DzzwgKfKFRERERGRPCjUIaRz58507tw52/21atVi0qRJbqxIREREREQKqlB3xxIRERERkeJHIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNxKIURERERERNzKx9MFiIiISMlmSkvjzI4dxB0+THpyMmaTCZ+AAILKlyeoVi1PlyciLqAQIiIiIm5lMZuJXrmSfYsXE7N1K2d27MCUmmr32BQ31yYi7qEQIiIiIm6RdP4822fM4K9p07h48KCnyxERD1IIEREREZcym0xs/uADfhs9mvTkZE+XIyKFgEKIiIiIuMy5PXtY0rcvJzdv9nQpIlKIKISIiIiIS2ybMYOVzzxDeopGdoiILU3RKyIiIk63cfx4lvXvrwAiInYphIiIiIhTbf7oI1YPHerpMkSkEFMIEREREafZvWABq1580dNliEghpxAiIiIiTpF45gw/PvkkWCyeLkVECjmFEBEREXGKFU8+SdK5c54uQ0SKAIUQERERKbBd8+ezd+FCT5chIkWEQoiIiIgU2O+jR3u6BBEpQhRCREREpEAOr1nDud27PV2GiBQhCiEiIiJSIFunTPF0CSJSxCiEiIiISL4lnj3LvkWLPF2GiBQxCiEiIiKSb8c3bMCcnu7pMkSkiPHxdAEixdEpkvmek+wmnkC8aU8FOlIBP+V+ESlmYrZu9XQJIlIEKYSIONmvnOV1dmHEnLltHef5luNMoDFh+HqwOhER5zqlECIi+aC3ZUWc6ALGLAEkwx4S+Ij9HqhKRMR1Tv/9t6dLEJEiKM93QiwWC3v37mXPnj0kJCQQGhpKnTp1qFu3LgaDwZU1ihQZSzllN4BkWM0ZnieKUvi5sSopNk7/BdunwKnNYEqBoAio2wvqPwp+oZ6uTkoorZAuIvmRawgxmUwsXLiQ7777jgsXLlClShWCg4NJTEzk+PHjlClThvvuu4/u3bvj46PeXVKyHSQxx/1pWDhGskKIOCb5PCzrCUdX2W6/eABOroffh0G7CdCwn3Oul5YKp/eCMRECwyGiLnjpxrlkZTGbMaeleboMESmCck0Njz/+ONWrV2fIkCFcf/31+PldfvFkNBr5+++/WbZsGUuXLmX27NkuLVaksAvNw83FEA3FEkekxsG37eHcjuyPSbsEP/cHsxEaP5n/a50/DL9NhQ0z4dIV726Xqwm3DITW/SCkbP7bl2LH4OUFBgNYLJ4uRUSKmFxfDQ0fPpw6derY3efn50fLli1p2bIl+/btc3pxIkXN7UTwHSey3V+bEGoQ7MaKpMhbNyLnAHKl1c/CNV0g/BrHr/PXApjVG9JSsu47dxC+Hwq/fAhPLYXqNzjevhRbPgEBpCcne7oMESlicr2/nl0AuVrt2rULXIxIUdeQcG6lgt19Phh4hlpurkiKNGMC7J6T9+MtJvhnmuPX2bEMZvS0H0CuFH8aJt4Gp3Y7fg0ptkrVqOHpEkSkCHKok29sbCybNm3i559/ZtOmTcTGxrqqLhEAVkdDty+h8STo9x3sOePpinL3OvV5ghqUvWLcR1NKMYkm3EAZD1YmRc6+b61BxBE7Zzp2vCkN5j4BZlPejk++CN8869g1pFir2KyZp0sQkSIoT53Tf//9d+bPn8/OnTsJCAjIHJiemppKvXr16NGjB7fccoura81i1qxZfPnllzbjVFq3bs2oUaPcXos437zt8PD8y12N/4mBb/+F9QOhYaRna8uJNwYe4xp6U40zpBKAN2U0EF3yI+6w4+cknYG0JPANytvx276HuFOOXWPfaojZC5F5u1MuxVtks2bs+PJLT5chIkVMriHkhRdeIDExkTvuuIPhw4dTsWLFzH0xMTFs3ryZr776iu+//54PP/zQpcXaU79+fT7++GO3X1dcy2yG4T9lHeuYkApvr4FvenmmLkf44EUlAj1dhhRlhnzOSGXwzvux62c43r7FAhtmwL3jHT9Xih3dCRGR/Mg1hNxzzz20adPG7r7IyEi6detGt27d+P33351enLhOXFwKmzYd55prSlGnTjlPl5PFyQQ4nE1vv/VH3FuLiMeUref4OWHXgI9/3o8/d9DxaxTkPCl2KrVsSXBkJIkxMTke5wUY/vsAsPz3kf3KSiJSnOX6Nlt2AeRqnuiOBXDgwAHuueceevbsyVtvvcWpUw52KyiB5s3bQeXKH9K581zq1p1Mz54LSE8vXH8GSgdCoK/9fZXC3FuLiMdEdYfA8o6d02iggxfJ59SqlsL1O0M8x9vXlyb9+9vfB/gBAYD/f5/7/vfh99+2wP/+deD+nYgUAw4vWJCcnMyRI0dISkqy2d60aVOnFZVXbdu2pXPnzkRERHDu3DmmTZvGSy+9xOeff05goLrB2HPhQjL9+v1AcnJ65rZvvtlJ+/bXMHBg4Zl2M9gP+lwPn27Juu/pVu6vR8QjvP2g0SDY/FbejvcNhgYOLlhYqjKcjXa8tlJVHD9HipT0xETO//EHF7ZuJeHgQSxpafgEBxPeoAFlmjWjdKNG1nVCgKYDB7LhnXewmKwTHHhhDRp57VDohTWUWIA0II/TJIhIEeZQCFm3bh3jxo3LEkAMBgO//PKLUwvLixpXTAtYvnx5XnnlFbp27cq///5L8+bN7Z4zfPjwzIHsnTp1olOnTm6ptbBYt+6oTQDJ8NNPB50eQmLOwsrfoGwp6NwWfByMvB/dCSnpMHc7pJkgPACGtoE+7s+7Ip7TahTEbIIjP+d8nJcv3Pk1BDl456Tlo7D/t3zU1cfxc6RIuLhzJ3s//pjDX31FemJitseF1KhB1MCBRA0YQFiVKtTu1o29CxfiSz7e4fyPAWsY2QVkTAStQCJSPDn0e+LTTz+lT58+3HXXXQQEBLiqpnwzGAwYDAYsOazcOnbsWMLCSm5/nipV7D/3qlWd+38y7wfo+woYjdbH9aLgl6+gov0lNOwK8IUv7ofxneFEPFxbznqHRKRE8faFbktg7WDY+QWYUrMeU+pa6DgFqt/qePvNH4LvXrJOvZtX17SA6hqMXNyYUlP5Z/Rodr/3XuYdjZxcOnSIbcOGseejj2gxdSo3jxzJkUWLMJgL3lWvPlAXMALJwOYCtygihY1DISQ2NpYHHnjAVbU4bM2aNTRt2pTw8HAuXLjAp59+SunSpWnQoIGnSyu0mjatSJcuUaxYcSBzW3i4P88808Jp10hMgkEjLwcQgN0HYPREmDbG8fbKh1g/SqJ9HGcFf3CY05QihDY0pB2NMGQO7ZQSwccfbp0Krd+Ef2dCzGZIT4HgCKjTC6rfBoZ8fk/4BcL9H8CXeezG5eMP97t/JkRxreSYGNZ06ULstm0On5ty+jS/3XsvREQ4JYBkyOiipbXYRYonh0JIw4YNOXDgAFFRUa6qxyGrVq1i4sSJpKSkEBoaSqNGjfjggw8ICsrj/Pgl1MKFPZg69Q9WrTrENdeEM3hwK6KinLeI3l//QvylrNvXbHTaJUqEnRzhExZj/m/gcBKpfM1aznCRHrT1cHXiEUHlocVQ57fb+nFIjoPvXsw6L/aV/IKg/3yodZPzaxCPSTl7llXt2xO/Z0/BGjp9mjAg3ilVWWUEEREpfhwKIfXr12fkyJF07dqVMmVsX7TecccdTi0sL8aMycfb6kJAgA9DhtzIkCE3uqT96pXBy8u61seValR1yeWKLDMWtnKSg8QSQQitqYrfFfPDLGZDZgC50lr+4XaaUZoQUknlL7ZxkIMEEEBTmlCDGlnOEclVxyFwTUtYM9G6gKEp7fK+gFBo8Qi0f14LFBZDGx59tOAB5D+BWAeWO/PuhWbNEimeHAohK1aswGAwsGzZMpvtBoPBIyFECqdqleGx+2HG/MvbfH1h+FOeq6mwSSKNN1jLPs5nbpvHDt6kPRGEYCSdI5yxe64ZM9Gc5DqqMp2ZnOZ05r6/+JuOdKAD7Vz9FDLFsYME9uBLacpxC944sEaFi1jSN0H6v+BVBnzvwGAofGPYCqVara0fcTFw5A9ITYTAcIi62RpEpNiJnjmTUz/+6NQ2Q4BUtP6HiOTMoRAyb948V9Uhxcy0MdDqevhhFZQrDU89Ajc08nRVhccCdtkEEIAzJPI5fzGCNvjghT++pJJm9/xgAvid9TYBJMNq1nA9TShNKVeUnimNOHYwjDi2Z27zIYzreIMytHTptbNjMe2DSz3B9PfljYYyWII+wODf1yM1FUnhkdDoLk9XIS5mTktj24gRTm/XC2sQcWa3LBEpfvI6hbeIQ7y9oX8P+GE6zByvAHK1jRyzu/0vTpFKOl540Yq6do8pSyh1qMJu7HefsGBhTzb7nGkPY20CCEA68ezgVVKvCljuYDFfhPgOtgEEwHIBEh/HYlzi9ppECrOj331HSi6rnOdXAGj6DBHJkUN3QsaOHZvtvuHDhxe4GJGSIqfZrTL2decmTnGBfZzI3BdGEAO5Ey+88tSGq6QQwznW291nJoVTLOUa3LyOhHEWWE5ks9MCyWPAT+/ui2Q4PHeuy9o2YF0FPcVlVxCRos6hEOLtbTs87Ny5c/zzzz+0bauZekQccTNV+ZZdWbY3p3Lm4PRA/HmR+9nHcY5whlIE04Ra+P73Y9uA+na7Y3nhRX3qubT+ZE6QU4/v5Gzu9DgilV9IYhrp7MGLygTxKAH0zD5gGZfn3KBpMxbzBQxezpsJTqSoSL90ibg//iB+2zbS4+LwCgjg4u+/YwA70184hy8KISKSPYdCyNChWaeGXLt2Lf/884/TChIp7JYRz3fEYQDuI5w7cHyhx/uoz27O8e8Vg88rE0p/si4HX5sq1KZKlu030Zp97OcYxzO3GTBwB50Jy0dNjgigEuTw8iWQygVqP5FPSWBY5mMTR4ljI0Y2Ec7EbM7Ky0spDZWVkiUxOpqD48Zx6n//w5SUZLPPH+v0t2lYw0K6k6+d31XTRaRkKPDviDZt2vDRRx/x3HPPOaMekULtcy7w9hXB4RcSeZ10HsOxd9f98eEtOrCD0xwklkhCuIFKeDswTMsff56gH/+yk2gOEkgA19OESCIdqiU/AqlIWVpz3k6XLC/8qUjXfLdt5gIJjLa7L5nZBPIofthZrdv3NkhflX3D3tdj8CqX77pEiprDkyaxd9gwzMnZT5hrwBpE/LAGkaRsj3ScBp2KSE4KHEJ27tyZpZuWiKddTINX98PP5yEqCN6MghbhBWvTgoVP7Qy4nsJ5h0NIhoZE0JCIfNfkjTeNaURj3D/yvy7D2cErxLMzc5sPIdRnNP6Uz3e7KSwnp04cKSy0H0L8+0HKB2CxP7UxAcPsb5ecxeyHNVPg7x+sU/aWrwlt+kHLXuCvhWELq13PP8+RSZMcOicA65ocCS6pSETElkMh5LnnnsNguNwfOzk5mUOHDvHwww87vTCRgnh4Byw/Z/08Ohk2XIQDN0OFAixhYcTCOUxZtp/FRBoWfEvYXDB+lKYZ07nI3ySwF19KUZ62eBNYwJZz7kVuyWYZNINXWSyhP8Gl+8F84Io9wRA0FoP/gwWsqwRaNwtmPQGmKzrqxJ+G6I2w8gN48Scok7WroHjW4UmTHA4gGXyxTq97yQl1qPOjiOTEoRDSrJntu4+BgYHUqVOHxo0bO7UokYI4a7wcQDIkmGDhGRhUgFXb/fGiBYFsuepF8I0ElbgAcqVSXE8prndae37ckuN+f7KfCMPg0xhL+D5I+wlMO/9brLA7Bq8C3gYrifb8CjP7gSWbl5Ind8OErjD6L/BSx5vCIjE6mn2vvlqgNjK6ZxkLWIv9VY5ERKwcCiF9+rh5yk2RfPAzgK8B0q4apxzshF6DY4ikL8c48d8Qzsr48FYBulOVdBZMpLIDA374Ux8AH+oQwL2ksDDL8T40wJ87c2zTYDCAXyegkytKLjmWv5t9AMlwbDvsWAGNc/6aiPscHDcuywD0/AhEIUREXMuhELJjxw7Kly9PZGQkFy9eZMqUKXh7ezNo0CDCw/VOoxQO4b7weGWYdnnSKKoHwL0VCt72tfizllps/m/4ZiuC8CnBd0EKIpFfOMurpHMKAD9qU4EJBNCQcKbiRSTJfImFBMCPALoRxjgMmnPH9eLPwr8r83bs+jkKIYVE+qVLnPrf/5zSljfWFwj5nTHLDKQ6pRIRKa4cuoc+YcIE0tOtv5I+++wzzp07x4ULF5iUz76nIq4yuR5MrQfdK8DL18DGlhDspNeuvhi4mWBuJlgBJJ/SOEwMAzMDCICRfZyiN2aSMOBPGGOpwF7K8QcV2EcppuNFWQ9WXYJcOpf7XZAM8VnXqhHPiNuyxSl3QTL4FuDc7OfjEhGxcuhl2enTp6lSxToIcdOmTUyfPp2AgAB69+7tkuJE8svbYB3/UZAxIOI68XyDxc77pCbOc4llhPEAAAaC8OFa5xdgsYChBAXIPUth81Q4vQP8Q6HBg9DySQjJ5vZgcOm8tx2sxR8Li/ht25zaXn57sJqARGcWIiJ2WSwW5s2bx4YNGzAajdx66610794dX1/rWwgxMTHMmzePgIAAAgIC6NixI9WqVfNw1Zc5vGJ6eno6x48fJyQkhLJly2KxWEhJ0ZqoIpJ36XZWes9gIpspdgvKdAHiRkDSXLAkQsCtEPYW+LfIX3snZ8KJKZB8CIJqQ5VnIfIh59ZcUMc2wYohcHST7VqOq0fDH9Og32ooXzfreeGRUKcN7P0t92u07OmsaqWA0uPinNpefmN6PM5dhd3b1xfSNMJE5GqLFi2iYcOGPPTQQ5w6dYrRo0ezZs0axo8fz7Fjx1i4cCGDBw8mJCQEgB9++IEyZcpkPvY0h7pj1atXj4kTJzJ9+nRatWoFwKlTpwgLc+3qzCLieanEYXbSmsoB3JDDvqyrxheYJR3O3gaJn4IlATBDyk9wtj0Ydzje3r7nYU8/SNgK6RcgfhPsehgOvub00vMl/gRMvxE+vxFObLK+pe2N7W/8hFPwTa/s2+j8Uu7XKV8TmnYvYLHiLF4BAU5tLz9BIoGCD2i/WvMCzvYlUhxZLBa8vLxo2LAhABUrVuS9994jISGB4cOHM2fOHIYOHWoTOO688042btzoqZKzcCiEvPDCC6SkpBAWFpbZBWvPnj107NjRJcWJiHOYMZGWz2Gi59jBzwxiGb1YzsPstzNrlaNC6Y4ftbNsD6IdgdxY4PazSF4MaX9l3W5JgoR3HWsraT8c/9j+viPvQOop+/vcxZQGc26H45tstxuw/sa/8rf+qW1wZIP9dprcBfePy/46pSvD4KXgU5CRA+JMIddd59T2sq6KlD0vPz/KdO7s1BXXAeo+9BDNhgxxcqsiRV9sbCwVK1a02RYWFsbbb7/NgQMHiIyMzOyWlcHb2xuTyZGfbNdyqDtWhQoVGDFihM22du3aqTuWFHlmi4U3EtOYm5JOZS8D74T40drPCXP6FgJb+Ik/+IlUkqlOXW7nEULJW59/IwlsYDTp/w0zNZLADj4nmEgq0TrX882kcZFfSOYQwVxHODdjwAsvAqnMAmKZQiI/YcCPELpRmicK9FyzfyI5vPOT6uC7QmcXke17xJZ0OLcEKg9wrE1n2v09nN2V/f6r+9ic3ArVs/la3jkUat8Mqz6GbT+AMRnKVoc2/aHdQAgr77SypeBKtWhhXbPF7JxlAvN637NMs2bcOGsWpRo0IPz119n05ptOuX693r3pNHMml5w42F6kuAgNDeX8+fNZtp8/f57atWuzePFirrvuOm677bbMfUajEaPR2fcq86/A8wWlp6fzwQcfcMcddzijHhGP+DQ5nbcTrX2OD5ks3H0xhSPlgwgu4oOX97ONdSzOfHyEPaxgNg8yOE/nn2JTZgC50lFW5xpCTCSxl34k8m/mtnBu4VomYcAHb0pTjhGUY0QOrTiJVw7zM3tHgPEwXJgGxoMQ0BjKPAE+2b3Azu0FXjb7zUY4v9k6IL5MK/By0VTD+5bmvN/w30dGjvLJpQvPtTdZP8D64lYLExZa/pGRlO/UibMrVhS4La/gYIKqVydul/1Aa/DxoeLtt1P7qaeo1KULhv++L1q/8QaRzZvz84ABJJ7K311Bv7Aw2n7wAQ379893/SLFna+vL5cuXSI+Pj5zWMTOnTtZsmQJ77//Pl9//TXvvfcePj4+tG/fHpPJxMyZM7n33ns9XPllTvkraLE4cwiaiPutMtrenoy1wNY0M22K+N2Q/WzLsu04+0nmEoHkZWBadiEs93B2hvk2AQQgjt+5wCrK0jkP13ai4N4Q9xp2Vy7wbQMHGoD5v/l84ufDhY+hxm/gb2dmrrJdIXpYNhfygrJ23pA58iX88xKk/jfoPqASNJkAVR7Ix5PJhcmBAbxePlCnqwPHK4AUdjVeecUpIaTmyy9z7euvk3LuHLF//UXCwYOYjUZ8Q0IIb9CA0o0a4Z3NGJSaXbvSZ+dONrz+OrtmzcKYkJCna3r7+1P7gQe4eexYQqtqakOR3HTr1o1vv/2W1NRU4uLiKFWqFMOHD8fX15fevXtTo0YN5s6dy+zZswkPD6dv375UqOCERdOcxCkhxFDE3y0WaerjxeLUy0HEH6jvU/RfcAUQlGWbD3744Jen8yvSCl+CSbtqws3q3JrruZfYmu12t4cQ74pQ7hs43/u/gekAXhDyNJz/8XIAyZAeA6eHQ7Vvs7YVch1UfAxOfZF1X5VnIeCq6Q/PrIE/+mDThSvlJGzuBYHVoGzLgjyzrKrfAv9+nfMxGaVUvRGC3Lj2yon9cGw3BIZCg1vA28E/QTFHYecW692khjdC+UquqbMIK9uuHVWfeIJj06fnu43Qhg2p9d9g8IBy5ah4++1UzOWcqwWULk2HSZO4eexYdn/1FYeWL+f01q0knjxpe1zZskQ0a0a1jh257rHHCCqvLn4ieeXn58fDDz+c7f6bb76Zm2++2Y0VOUZLD4sALwT78m+6me9STZT3MvBhqB/lvIp+uG5CW3ayyWZQehPa4JvHEOJHCK15k21MIY5o/ClNXXpSkdxfOPtj/53M7La7XGA3qHQSkn8AyyXrFL2WADiRzSDzhBy6NdX9HIIbwImpkHzw8hS9lZ/Meuz+j7A7hsRiggMToOy8/Dyb7DXuDWteh6Rz9vdbrvg4+DvM7ASPrQBf586sZOPYHpj6NGxffXlbmYrw4HC465ncz485Cu89C+uWXh7v4O0D7bvDSx9D2QjX1F1E1f3oIxJ27uTihmwmHciBX4UKXL9gAV5+efsdARB/+DB7vviChIMHSUtMxDc0lNL16lH3sccIioig8aBBNB40CICkM2dIiY3FYjbjFxZGaOXKDtcoIsWDYc2aNbn2pRo7dmy2+ywWC6tXr+aXX35xamHOlpiYSNeuXYmLi9OUwpKtFIsFP8CrGN3dO8dJ/mYtSSQQRSPq0wpDPlYASCMJHwIw5HFSvRSOsYsemLjcFcOPSK5jAT6EO3x9l0g/D3vKYzckeJeGehcKfo0VNSHxkP19YQ3g9nxMEZybE3/A3Dsh8aztdjPW0cZXP93O70C77LqYFVDMIXihJcSdtb+/99vQM4dxQWdOwOM3wulj9vdXjYKZG6FUuYLXWoykX7rEP48+yunvv8/zOSENGtD0u+8Irp115jp7jq9axT8TJnB0xQosdgbDe/n5UfPee2k0ZAgRLfK5Hg8QHx9PeHg4S5cuJTg4ON/tiEjhkqdXE97e3tl++Pj4cPvtt7u6ThG3CDAYilUAAShHJW7jIboxkOu4MV8BBMCXoDwHEIAAqlKPuZTjXkJoQgUeph5fFZ4AAuBTFkKy6RoW/ohzrhGYw+q0QXb2JZ+B0xshKSb/16zcHJ4/CF0/hUYPQ2QzSMP6Ye9tpz9m5P9aufn6rewDSMb++KwzvGT67PXsAwjAsQPwxTv5r6+Y8gkJoenChTT+3/9yDRW+ZcsSNXo0N23dmucAsuW111hy220cWbbMbgABMBuNHPj6a75v3Zqd06Y5/BxEpHjL052Q4kB3QkTELuMxOHI7pO65vC3oZqi+DLyd8Lvi2DewOZtVxVsvgUr/DQxPT4GNT8OBL8GcBgYfqNkTbp4GPlnH9gBw8SBsHg+Hf7LOclW3B9zwAviH2h636GnYNCX7Gn0C4O2ss6AVmDEFepaB1FzaHjgJ7n426/ZL8dClIqTkMkVreBlYfhL8/PNfazFmsVi4sHYtF379lfht20i/eBGvgABCrruOUq1aUeGuu7IdZG7P5pEj+WvMGIfraPPpp1w3cKDD5+lOiJRUL730Elu32h9f6WzDhw+3mc7XHTQmRERKNr+qEPUvJCwHYzQENIGQds5rv2oPuLgd9o23jgMBa8Co9xoYE+H3nmDwgqTzcOyny+dZ0iH6K+u+trOztnthH8xtDclX3EVYPxoOLIGHfgPfK4JL2aica8xtf35dis09gABcOGl/+4mDuQcQgLgL1m5bVWo6Vl8JYTAYKNu+PWXbty9wW4eXLMlXAAH4/emnqXDDDZRv1qzAdYiUBPfee69NCOnQoQND8rh4Z3p6euZHUlIS8fHxxMfHc+TIEaKjo9m5cyfnzl0eO7hw4cLCF0KGDx9O//79qVkz+1/u0dHRzJgxI8exIyIihZbBG8Lucl37DcdCrSfh1DJrqIi8E7a+AH+9ZnucD1lXiIv+H9R9GoIqQugVg/rXv2EbQDKc3gr/zoHrB13e1vRR+GkEpGUTCFoNsr+9oEJKg39g7kHk7HH7233zPjhad0HcY/sHH+T7XIvJxD8TJ9JxzhwnViRSfLVu3Zp7772XhQsXArB69Wpuvvlm2jvhDQWz2cxff/3Fjz/+yJo1a9izZw+7du2ifv36BW47r3LtjrVq1SqmTZtG+fLlad68OTVq1CA4OJjExEQOHz7MH3/8wZkzZxg0aBAdO3Z0V90OU3csEXELUzIcngJnloNPGFR7AiKuWjvk5E+wupP9841kHbeR+t+2yBuh3WdQtgF8FAppl+y3UaMzPHDVWhH/fg/zeoLpqtVyG/WAnnPBy0Vr4nz0GKyalf1+C4AXfPI31Gxku89shu5RcDKbgf0Zrm0E/9tewEIlNxd27uSbBg0K1Ia3vz+9jx8nsFzeJxJQdywpyYxGI8888wz79+8HIDg4mOnTp1OxoqMTZ2cvOjqaN954g2uvvZZRo0Y5rd3c5DrK9NZbb2Xu3Lncdddd7N+/n48//pgRI0bw8ccfs3fvXu68807mzp1bqAOIiIhbWMyw+Q7Y9RKcWw0xi2DLnXBwku1xJ5dn38bVv5UzptMFiNkIi28FY4L1jkp2DHYCRYPu8OIeaDsUaneGJg9Zp+btNc91AQSg50gIy2EtEgvWsPGDnamSvbygh52xIlfr8Vy+y5O821WAtUcymFJT2TvbTvdCEbHLz8+PUaNGERgYCFjfVH/rrbdIT7/6tnn+1apVi08++YQDBw7YdNFytTyNCfHz86NLly506dLF1fWIiAcc5g/2sQZ/QmjE3ZSmiqdLKprO/Ajn12bdvm80VH8CvK1/RPBx4N1c01WPk0/D/nlQ+174d5b9c2rfa397mRrQZVzer+0MFWvBoE/g3V7YTMx25XolALs32j+/5/Ow609Y+T/7++8dCN36Oa9eydbFPXtyPygv7ezd65R2REqKqlWr8vzzzzNunPX39+7du/n8888ZNMh5XWnDwsJ45ZVX+OGHH3j88ced1m5Oiv6S0CJSICf5l9VM4DjbiWY9K3ibNFI8XVbRFJfNLCZpsZB0RZeiax4Gu1MlG8Dw3yxFBm/r+BB7b3bFRcNNoyHEzorhVdpA/exX0HW7/X/C8o+tb3n5YP2rY8G6ZsmV3c6O7oH3Hofj+2zP9/KCt76CMV/D9W2s40R8/aDFrfDe9/Dqp+56JiVeWkJC7gflpZ34eKe0I1KSdOrUyWbg+Pz589myZYtTr3Hddde5ddmNXO+EvPvuu3lqaOjQoQUuRkTc7zB/2DxOIYEY9lCVJp4pqCgLvtb+dq9ACLji7lKp+tBiMvz5vHU6XgAvf2g5Dap2h0tHIfo72DLafntlG0F4dXj0T9g6EQ6tBJ9A6xS9TQaCj4sHae9dB7tXg38wtHgQyla1f9z+P2FEW0hNupy5vP77SMMaRDKkm2DlF7D+e/hgLdRqfHmfwQC397B+OMqYCquXwOkTEFkV2ncFB1YDl8t8Q0NzP8iN7YiUNEOGDGH37t0cP34ci8XCuHHjmD59OmXL5tDl1UFVqrivJ0SuIcRkutwXIC0tjXXr1lGzZk0qVqxITEwM0dHR3HLLLS4tUkRcJ5QKedomeVDxXgipC5eu6rZS42nwvWpCjNpPQtV74fgS6/iOKndDwH+Ddcs0gMBI+PfTrIsWlqoNUfdbPw+pCG3HWT/cwZgCk7rDPz9e3vbNUOj1AXR6Puvx/xtlDSD2+GAdhA+23bIuXYSRdwMBkG6E69tBjxegVkPH6/1+Dox7EWKv6ONctgKMmAhds1m7RbIVVquWU9oJj3LRlNAixVxgYCCjRo3imWeeIS0tjdjYWMaOHct7772Hl1fR69zk0GKF77//PnXr1qVr166Z25YtW8bu3bt56aWXXFKgs2h2LBH70kjhFz7iJP/ihTeN6EZT7vN0WUVX6mnY8xqcWQY+4daxIDWet76b76iL+2D9y3BkmbV7Vq374ab3ITiHWVFMRki+AEHlnT/g/JuhsGy8/X2vb4ZaLS4/Tk2CHiFgyeFPjBHr3ZCrx70ANj0C/QLgnUXQMpsZxez5eRE8c6/963t5wbSl0FbjHB1x9u+/WdC0aYHa8PL1pffRowRFRub5HM2OJWJrwYIFTJ48OfNxv379eOSRRzxYUf44FELuvvtuvv/+e7y9L/9hM5lMdO/enR9++MElBTqLQohIzhI4iy8BBKCuEoWOOd16tySnGbEAfh8Dmz+0hpCwqtDuTWjc1zk1WCzwVDlIvGB/f9v+0O+K2ZOS4qFXeM5tJmPbJetKaVweOwJQoQp8exi88xis7moMe//Jfn/jlvDtpry1JZkW3ngjpzfl//+t1oMPcvs33zh0jkKISFavvvoqm/77WfTy8mLChAk0bJiPO8Ye5NC9m6CgIPZeNavF3r17M6cNE5GiK5TyCiCFlZdP7gFkyyRYO9IaQADij8EPj8GBH3M+L6/SUrIPIACxJ+D4HviwJzwcAk9eA6VyWAviyoBhb58B61+ojMxx5jj8ucr+8WYzrF8JcyfBT9/C3h05BxCA7ZvhxJGcj5EsGuVxtWa7DAYaPW+n256IOGzYsGGU+2+9HbPZzJgxY0hw0uQR7pKnKXoz3HPPPbz66qt06dKFyMhIYmJiWLFiBT17qm+tiIhH/fGJ/e1/Toaozrmfn5oIK9+Gv74Bswma3AedR0FQaet+v0CIuBZO77d/ftlqMOImuJQRVBKt/3hjv7tVaGTW8S4Zrrw/b/jvwwL8ughqNIAKlS/v37sdXrwfjh64ou1S1gCTXcjJkJTNYo+SragHHyRm/Xp2TJqU+8FXuXH8eCJbt3ZBVSIlT3h4OCNGjODFF1/EbDZz+vRpxo8fz1tvveXp0vLMoTshPXv25NlnnyU6OppFixYRHR3NM888Q69evVxVn4iI5MWlU45t3/w5TGoF79aFxUNg8u3w8zg4fwhij8Kaj+CTW8F0xRzBd7xsvy0ff9i26vJdmCsZgOArusAGhMADI2D8OihnZxYWM1nDQ8Zwmu8/hXuvgfeeApMJ4mNhwG22AQQg4SL4Yn8W5AzBoVD5mhwOkOzc9NFHNHTkjobBQKt336VJIR87KlLUNGnShIcfvjwl+7p16/j+++89WJFjHLoTAtYV1G+99VZX1CIiIvlVvS3sX2Zne7us234cBavevvz47F77dw2O/QU7foAm98KF47B9MXj/d1siozuVwQdSUiEl2voXxZvL4zky1LsR7h1hne3q2hYQ9F+3v7eWwIDrL4cFy1XnXSljuykdFk6F0NIQVBZiz9o/3oC1nrRs2rvnUQjS+IL8MHh5cfOECVRq04Z/Jkzg1O+/Z3tc1c6dafzCC1Tp2NHNVYqUDH369GHbtm3s2LEDgKlTp9KwYUOiisAsdA6HEIDk5GQuXryI5YpZRypVsrNoloiIuEf7sXBsPaRcvLytVA248ap3n1Pi4dcPs56fsX7H1WHkxDZo1A3evw1OXTH1sAHw9YHkq1ZTtPfif88GKFsZImvaHnttE7i2Kez7K+fnZm/8yMIpUK9Nzudldyfkuqbwwpicz5Vc1bz3Xmreey/n//2XPTNmEH/wIGmXLuEXFkapevWo378/YTVr5t6QiOSbt7c3I0eO5IknniA+Pp60tDS++OILxowp/L/jHAohMTExvP322+zevTvLvl9++cVpRYmI5CbNNJs007uAN37eo/HxLuHTCkc0ggH/wF/T4MJ+qNgMru8PgWVsj7twCNKyWbvD3ov2crVgx4+2ASSDOd3+mI+rO/qmJMCk3jB2fdY2XvgUXrwNEuPs12TB/qrxCRdtu4rl1V0PwVuf6S6IE5Vt0ICbPvrI02WIlFgVKlTg5ZdfZtSoUZQqVYqnn37a0yXliUMh5OOPPyY8PJxPP/2UIUOGMGHCBGbMmEG7du1cVJ6ISFYm8yaM6f3I6KOTmv4QXoa/8PK6zrOFeVp4VWj/ds7HlKkBvkH2g8jVXaFKV4PrH4TfZhS8tr0b4OR+qHTVqvJ1m8P0v+CbD2D113Ap1jodsMFgXQX+UnaByQAdusNvy7O/pr0uZt16OyWAWFJSMAQEAJC27ndSJn5E+q9rwWDA9/ZOBDz/Aj433FDg64iI5MXFixfx8/Pj7bffdkrvpISEBGbPno3JZOLEiRO0a9eOW2+9lU8//RSLxUJCQgKPPPII11xzTb6v4dDA9F27djF06FCuvdb6RyQqKooXX3yRBQsW5LsAERFHmS1/Y/uKOR2zJZcpWcUqIAza2JlmNbg8VG9lXeDQ4AUN7oLGPWBEFHw7JPu/FvbGcFiwf3zSRfttVKoJQybDkvPwwwVrKFlwAnoOzf55tLgd7u4LTW/Jvq6rb5RUqAit8z+m0ZKaStr4MaTUrUZKuUCSq5QmsXNHEjq0JW3R91hiY7FcuIDx63nE33IjxgXf5vtaIiJ59ccffzBx4kReffVVrruu4G/GpaWlMWHCBHr27Mnzzz/PCy+8wPvvv88bb7zBAw88QOvWrVm7di2LFy8u0HUcCiEmkylzob+AgABSUlIoV64cJ0+eLFARIiKO8DI0xrbvkA9ehqK1SJNHdXkbHpgOVVtAQLj1jkHcWTi6Cao1hTExUPVG+PE9uHjS2u3qyjU7MlS4FkIibLdZsHbPyphaN0NYeajeKPfaQkvBtddD2Yrw8EvQ6Kasx5SvDC9NBh8fGDEZOna/fHfD2wfMXpB61Tne3jBykvWcfLCkpWG8vyvpb47EcvyYdVvsRVJ/WW1/Vfb0dC7164v54sV8XU9EJC8OHjzIG2+8weOPP+60nkk//PAD3bp1y1yHxM/PD4vFQmRkJBUrVsRsNlOlShU6dOhQoOs4FEIqVqzIwYMHAahevTo//PADy5cv1wrkIuJW3l6t8fP5DIMhCoOhDv4+X+Hl1cA1F0vcDkcGw/HXwZjNuhZFUcv+cOd4uBRn223p2B+w4VP4ZWLWcwxAqQgoWx1uex5GboDnvrGGDjPWf01XHZ/h4XfA19+xGgOC4ONfYOQX0PoOuKEjPPkOfLkdgsPhmTvgwUaw9ntIT4Xb74dVJ2DuBuh4tzVweHnBLZ1g1irofL9j17+C6Zu5mNfYLpaYbm/9kyslJWGc+1W+rykikpPz588zfPhw2rVr59TlMsLCwmjU6PKbRhkLlbdo0QKAli1bMmvWLLsrtA8YMIAVK1bk6ToOvSX00EMPcf78eWrWrEnv3r0ZMWIERqORlzT3t4i4ma/34/h6P+7aiyTtgJ0twfLf2+rnvoRGu8ArwLXXdZf9a+xv37caEs7Y33f7EOhyRTepms0gIBSS7azUGx4BUS2hy7PQOJ/doPz84c6+1o8rPdEe/lx7+bEpHX5ZAGXKw/ApMHWx9Q6FxWINIgVk+mpWlm32boBkOW9P1olcREQKKjk5meHDh1O1alWGDLHTxbYAbrvtNpvH27Ztw9vb227ouFqPHj24IY/j4RwKIW3bts38vEmTJixevJi0tDQCAwMdaUZEpGg4N+dyAAFIPQRxP0PpuzxXkzOVrmp/e5lqULcj7F6Vdd91t9s+DgiBbsPg6xG228tWgXf/snbDcrblc20DyJV++AKeHWvt1mUwWD+cwHLmdL7OM5Qq5ZTri4hkMJvNvP3226SnpzN69Gi8va/uK+tcf//9N9deey1BQUG5HtvRgTWBCvT2kI+PjwKIiBRf3qWzbvMpk3VbUXV9Dyhbw3abbwC0HQI9PoLQCrb7bh0C1a6//NiYAsf+hc7PwLNfQb1boFIduP0peGtD7gHk9BE4ujtvtxQyfDMZRvbOfn9qCsQcy3t7eWSolXXhr7z83ffr+ZDTaxGRkm3y5Mns3buXsWPHEhzs3OnGL141ji0+Pp7o6GiaNGlis33ZMjuL4zoo1zshvXr1wpCHd5L+97//FbgYEZFCpcIAODcbUvZZH5e+F0LtDJQuqvxD4Nnf4ecxcGiDdU2QjkOhShPr/jH74Y9vrF2z6t8ONZpfPnfdVzBnMFw6D/5B0G04vPFb3q57IQbefQi2/9cdrHJteOVLqNMi67HnY6whpVxFuHgePnwx59DiHwCR2dzhKQCfxwdiXLHUZpu3l7Wnl9neVMCA3yOP4uOEmWpERDIsXLiQ5cuXM3HiRCIiInI/wQHHjh1jwYIF7Nu3j2bNmtG/f3/++OMPzGYzdevWzTzu4sWL/Pvvv9x5552Z25YsWcK5c+cwm83069cvT9fLNYQ88sgjmZ+fPXuWxYsXc/vttxMZGUlMTAw///wz3bp1c+Q5iogUDb7loOG/EL8GfMIhpKWnK3K+UpXhgSn29wWGQZsnsm6POQCf9gXzfyOzU5Ng/ki4pik06ZL7Nd/rfTmAAJzYB691hTlHwP+/u+unDsPYx+DvtdbHDW+CW+4D49XTXl2lax9rVywn8+7SFe8BT2P6bLLNdj8fSA8tTXpiCqQkWzeGhhIwYBCBb491eh0iUnJt3LiRqVOnMnr0aGrXru3Utk0mE++//z6tWrXihx9+oGXLlqSmprJmzRrKlStHcrL191tycjIff/wxAwYMyDx3165dlCtXjkqVKjF16lTnhZArU86LL77I22+/bTMHcdu2bZkxwwkLWYmIFEZevlDq9tyPK0n+WHg5gFxpy4LcQ8jZ4/C3nbEmcWdhyzK45X7rnY6X74TDuy7v37Eejh/Ivl2DAbo+Ci9PyNNTyA+/Dz/BdFtn0j+fimn971gSEjBbgLhYvIOC8X3iOXw634Hvja0xhIa6rA4RKXn27dvHm2++yaBBg7jpJufekU9ISGDMmDHExcVxzz33cPz4cWJjY/nggw944oknSExM5PPPP2f79u2kp6fz0EMP2dyFMZlMtGrVig8++IDWrVvn+boODUzfs2ePze0YgDp16rBnzx5HmhERkaIsIMT+dv889E1Oy+FOhjHF+u+232wDSIbY0xASDJcSs+57bhz0fSXna584AgvnwIWz0LKtdeHCbVsgMBiatcrTLFreXbqStupn0pfbrtRuSEok/dOP8bv3PgUQEXGqs2fPMmLECO644w7uu+8+p7Wb0aPpu+++Iy4ujoEDBxIYGMjLL7+c5dgPP/ww23YaNmxIcnIya9eu5bPPPsvz9R0KIZGRkfz44482d0dWrlzp9D5pIiJSiLXqYe1+lRh7eZu3L7TLwy34SrWgRiM4dNUK934B0Py/uyhx57I///6B8PXUy12fAG69H3q/mPN1V/0Azz4ARqP18ZyPweANyf/d0aleE6Z8DU2aZ98GYElIwDg7m7v/FgvGyZPwublNzrWIiORRUlISw4YNo3bt2jz99NMOnWuxWDCZTCQnJxMfH09CQgLHjh0jOjqa3bt3s2PHDiz/jbHz8vLi9tvzf9d/zZo11K5dm0qVKrFjx448TefrUAgZOHAgI0eOZMmSJURGRnL69GkOHjzIW2+9le+iRUScZfTo0axatYp169blu43Dhw9To0YN9u/fT1RUFHPnzuXNN9/MXKypb9++pKen89VXJXgRutCyMOIX+N8rcGATVKoHD7wF1Rvn7fyX58Brd8K5E9bH/oHwwhcQVtb6uElb6/ogV4//8PaBh16A3sNg5dcQfwFa3gpNcumaYDTCq/0vB5AMFpN1FXgTcOQgPHonbD4MgdlPQ2k+fAgS7dyJydi/c0fOtYiI5JHJZOKNN97g4MGDHD16lM6dO+f5XLPZjMmU24qql7Vo0YIyZfI/++PGjRtp3749R48eJSHBzrpRdjgUQlq0aMHs2bP55ZdfOHv2LFFRUbz22mtUrFgxXwU7g8ViYdasWSxbtozExERq167N4MGDqVGjRu4ni0ihNGvWLEaOHMnx48ddep28hJaHH36Yhx9+2KV1FEnXXA/Df879uJiDEHsKopqDr591W83G8MVB+OsnSEmEprdD6BXTIZcqB0+Oh4nP27bV/y0oX9n6ea9n817rn+usXbDsyQghAOfPwrLv4P7spwA2VIjIcUosQ6XKea9LRCQHX3zxBVu2bAEgPT3dpddyJODY06pVKw4dOkRaWhr33ntvns5xKIQAVKxY0WbGLE/75ptvWLFiBePHj6dy5crMmTOHV155hTlz5mgNExERTzGZYNJj8OtX1oHmpSNh2EKoe6N1v68ftOya/fkPPAdN2sCqr8FihvYPQL2cu0ply5EV08/G5NxURAQ+d9xF+tLFdvf79cnbrDAiIrnp378//fv393QZeXLlUI28yvU3808//ZT5+fLly7P98JTFixfz4IMPUrNmTfz9/Xn88cdJS0vj999/91hNIiVNu3btePbZZ+nRowdhYWFUqVKFr7/+mh07dnDjjTcSGhpKixYtMrs0ffvttzRr1ozSpUtTrlw57r77bg4dOgTA77//zqBBgzh58iQhISGEhIQwd+5cAAwGAx9++CHNmzcnJCSEFi1a8Oeff2Zb18WLFxkwYABVqlShXLlydOnSJbOGuXPnMnbsWDZu3Jh5HXu/N2bNmkWVKlVsthmNRvr370+pUqWoVq0a48ePd8r/Y7Hy61xY++XlNT1iY2BiH8cWJry2CTw5Dp4an/8AAnDDzRBRyf6+q3srNLsx1+YCJk7Bq07dLNt9+/bDR4sTiojkSa4hZN68eZmff/nll3Y/PNU3+tKlS8TExFCvXr3Mbd7e3lx77bXs37/fIzWJlFRffvklTz/9NBcvXuTll1+mX79+DBs2jHnz5nH+/HmqVq3Kc889B0BoaCgzZ87k3Llz7NmzB4vFwkMPWV+83XLLLXz66adUqlSJS5cucenSJZvuUFOmTGH27NmcP3+eO+64gy5duhAXF2e3pt69e7N//37+/PNPjh49Sp06dbj11lsz2xw+fDg33nhj5nVuueWWPD3X77//nubNm3P27Fnmz5/PuHHjMoOS/OfftVm3ndwPF066vRR8fOC9ORB81axeJmxDSIcu0OLmXJvzqlSJ4M3bCJgxB99HH8N30DMErV5P4NTP87S4r4iI5KE71hdffJH5+ZWBpDBISkoCICTE9g9LSEhI5j4RcY/77ruPNm2sswL17duXwYMH88gjj3DNNdcA8NBDD2XeVr6y72m5cuV48803adq0KQkJCYTmMr3p888/T/369QF47bXXmDZtGj/88AO9e9v24z916hRLly5l27ZtREZGAvDee+/x1VdfsXTpUnr27Jnv59qoUSMGDhwIWPvBPvHEE8ycOVNjR65UuU7WbaFlIay8+2sBuKkjrD0ES+ZZx4fUqg/r1sDq5RAUDN0fhqdymeL3CgZ/f/we6g0PZT9+REREsufwmJArnThxAi8vL48NTA8Kss5gcunSJZvtly5doly5cnbPGT58OH5+1sGRnTp1olOnTq4tUqSEuPL3QHBwsN1tGTNm/Prrr7z55pvs2rWLxCtmGjpz5kyuIeTKSSe8vLyoXr06x44dy3JcxrZatWplbvP19aV69eocPXrUkaeWYw0Zj7///vsCtVnsdBoIq2fD8d3WxwYDPPrO5cHprnD2FHw4DNathAqVoN8rcMcVYbNMOehzxYD2u/MWRC3/rRRs0DhDt1m5ciUrV64ErN0fRUqyAQMG0L17d7p0yWUx2DxKTU1l2rRpREZGcuHCBUqVKlWgN+byy6EQ8u6773LHHXfQsGFDVq9ezZgxYwDrC/uOHTu6pMCchISEEBkZyZ49ezJXcTeZTBw4cIDbbrvN7jljx44lLCzMnWWKyBWMRiNdu3bltddeY9GiRYSGhvL333/TtGlTm/nKs3P48OHMz81mM0ePHs0yZgOgatWqAERHR9O4sXXq2PT0dI4ePUq1atVyvU5Orqwh47G9Gkq0kFLw4VZY/621C9YNd8I1uc8bn28WCwzoDHv/W3/k/Gl4qRcEh0JbxwdMAljOnMb00tNYli0CwND5Lrzfn4yhYjbjS8RprnyTMD4+nsmTJ3u4IhHP6dGjBzfccIPT2vvss88ICwvjwQcfBGDw4MFERETQvn17p10jLxz6C7xlyxauvfZawDqwdNSoUbzzzjv873//c0lxedGtWzfmz5/PoUOHSE1N5YsvvsDHxyfPfbtFxL2MRiPJycmULl2a0NBQTp48yciRI22OiYyM5Ny5c5w/fz7L+RMnTmT37t0YjUbGjBmD0Wjk7rvvznJcxYoVueOOO3jxxRc5ffo0ycnJDB06FD8/v8xZPCIjIzl69CgpKSkOPYft27fz+eefk56ezpYtW5g+fTqPPfaYQ22UCP6B0OFRuH+YawMIwNbfLweQK82bku8mTT3vwvLDd9aZvkwmLMsWkf7gnZlhWUTEHTp27Eh4eLhT2kpJSWH58uU2r5Nvvvlmm4mo3MWhEJKamkpAQABJSUkcP36cNm3a0KJFC06fPu2q+nLVo0cPOnXqxIsvvki3bt3YsWMH7777rqbnFSmkQkJC+Pzzz3n77bcJCQmhS5cuPPDAAzbHdOjQgW7dulG7dm1KlSpl80bHk08+Se/evSlTpgw//PADy5cvp1SpUnav9eWXX3LNNdfQtGlTqlSpws6dO1m1alVml68ePXpQp04dKlWqRKlSpfK8yGH37t3ZtGkT5cqV47777uOll14qVFOXl0iJ2SyOdcn+pAW5Mf+xCctff2TdsWMblo05f59YjEZSlywhee5cTB78+ygicrVjx46RkpJCRERE5rbIyEi2b9/u9jdYDGvWrMnzFR9++GHGjh3LoUOHWLFiBe+++y7Jyck8+OCDLFmyxJV1FlhiYiJdu3YlLi5O3bFEiiiDwcDPP//Mrbfe6ulSpLBJToL2VSA+1nb70A+hzxCHmzMvW4zp4Xvs7vOeNR+vex6wuy/t77+5eNddmE/8txq8ry8h775L8BDHaxCr+Ph4wsPDWbp0aeZ4M5GSYMmSJZw7dw6z2Uy/fs5Zg2jLli0MHTqUVatW4e3tDcCff/7Jyy+/7PafMYfGhNx///0MGjQIgFGjRgHwzz//ZM5+IyIi4hGBQTDxO3j1UYg5Dt7e0K0PPOzAyupXMLRsDf7+kJpqu8PXF0PrNnbPsVgsxPXocTmAAKSlcemFF/Bp2hS/1q0x+Prmqx4RKVl27dpFuXLlqFSpElOnTs0SQlJTU5k0aRIm09WLHdkyGAwMHDgws8dAfHw8Xl5emQEEwMfHGgeSkpIKbwjp3r07LVq0wNvbO3PKy0qVKvHCCy+4pDgREZE8a9kefj4MB3ZC2QgoF5HrKdkxlCuP16ixmEe+aLPd64XhGCrYbzf9zz8xZbNG1YV27bCEhhL8xBOEjRtnE0YsaWmYLl7Eu2xZDPmcLEFEiheTyUSrVq344IMPaN26dZb9/v7+vPzyyw63ay9kJP83A2DG7LHu4vBvu8qVK+Pv78/evXuJjY2latWqWaarFBFxBYvFoq5YkjNvb6jTqEABJLOpZ17Ae/UWDL36QDnr+ibmD8aQ/mRfLFffIQEs6ek5tmdJSODShx8S998LB4vZzPk33uBgZCSHKlTgcM2axM+eXeC6RaToa9iwISkpKaxdu9Zmba2CKlOmDGazOTN4gDWEBAQEuH24gkN3QhITExk3bhzr168HrLd4WrduzdChQ7MsGCgiIlKkzJwMX0yG5GS49yF4eTSGJs2wbNkA585aj0lPxzJvNuaKlfB+bazN6b4tWuBVpQrm48ezNG2+4vOk6dMJf+cdLrz/PhdGj87cnn7kCKf79sWrdGlC7Mz4JiIly5o1a6hduzaVKlVix44dNGx4eZbB/HbHioqKIjg4mJiYmMybCMePH6dOnToYDAaXPRd7HAoh06ZNIyEhgc8++4zKlStz4sQJpkyZwmeffaYuWSIiUnTNnQGvPnP58cSxkJyEpftDEJ21i5V5wbwsIcTg7U34nDlc7NYNS8Ll2brMV51rSUrCdOECcZMm2S3l4oQJCiEiwsaNG2nfvj1Hjx7NXOw3Q367Y3l7e9OmTRt+++03atSogcVi4ffff6dXr17OKjvPHOqOtXnzZl577TWioqIIDAwkKiqKkSNHsmnTJlfVJyIi4npzP7e7zZDdCu/+/nY3+7VvT7kjRwidNg2/++8nnawhxKd2bQxhYZjOnbPbRtpVi2GKSMnUqlUrjhw5wp9//ml3XEh+PfXUU5w4cYJ58+YxefJk2rdvT4cOHZzWfl45dCckNTU1S7erkJAQUu30jRURESky7M2Pb7FgaNgYw/U3YPn7T5tdXr2zny7Tq3RpggYMIPDRR0k/cQLjxo2Xd/r6EvbBB3iHhuJXrx7G3buznB/QokW+n4aIFB8ZC+s6W0hICMOHD3dJ245w6E5I3bp1mTlzJmaz9X0ds9nMrFmzqFOnjkuKExERcYtedla872nd5v2/xRjuvMc66L10GbyGvo7XMy9mPf4qhoAAyv3yC6U+/ZTABx4g+NlnqfDXXwR27QpAmbfegqtmwzKEhFCmELw4EBFxNYcWKzx69CgvvvgiJpOJiIgITp8+jbe3N++//z7Vq1d3ZZ0FpsUKRUQkR59+ZB2YnpIM9z0MQ9+y6XZlMZudPoVu0po1XPzoI9Kio/Fv1ozSQ4fif911Tr1GUafFCkWKJ4dCCFin8dqwYQNnz56lQoUKtGrViqCgIFfV5zQKISIiIkWPQohI8eTQmBCAwMBAOnbs6IpaREREir30Y8dIP3QIv2bN8MrDi+r08+eJW7YMg78/pe66C68i8MafiEhuHAohZrOZn376iT179pCUlGSzrzAMcBERESmsLBYLF55+moRp08BsxhAWRrkZMwi+//5sz7m4aBGHHnoIy38Li/mUL0/UihUENWvmrrJFRFzCoc6tEyZMYOrUqVy8eBFvb2+bDxEREcle0oIFJEydCv9N7mKJj+dcnz6YLl60e7w5OZkjjz+eGUAA0s+e5Uj//u4oV0TEpRy6E/Lrr78yefJkqlSp4qp6REREiqXkH3/Mss2SlETqb78RZGdxwsRNmzDFxmZtZ9s20mJi8I2MdEmdIiLu4NCdEF9fXypWrOiqWkRERIot76pVHdruU66c3e0Gf3+8rlqzS0SkqHEohNx1110sXLjQVbWIiIgUW6EDBuAVEWGzLfDOO/G//nq7xwc2bEhImzZZtpft0wdvhRARKeIc6o61detW9uzZw6JFiyh31Ts0EydOdGphIiIixYlPpUpU+uMP4idNIv3AAQI6dCB0wACbYywmE/GLF5P8998EtWhBjYULOfHSS1xcsACDvz9lH3uMSm+/7dB1LRYLFzdswHj+PGXatsU3PNyZT0tEJF8cCiHNmjWjmWbkEBERyRefqlUp89572e4/ct99xC9enPm4VO/eXDNnDnzxRb6ul3LqFH/ecQcJ27YB4B0cTIPp06nUq1e+2hMRcRaHQkifPn1cVYeIiEiJlrh+vU0AAbj45ZeUf/llAhs2zNyWcvAgJ8eMIemffwhu0YLKI0fil814zd3PPZcZQABMiYn807cvZTt2xL9CBZc8DxGRvMg1hFgsFgwGA2BdJyQ7Xl4ODS8RERGRKxijo7PdnhFC0uPj2XXzzaSdOgVA4p9/Er9qFQ3//RcvX1+b8ywWC6evCjUAFqORs8uWUblvX2KWL+f08uX4lS1L9cceI7hGDSc/KxER+3INIV27dmXZsmUA3HrrrZmB5Gq//PKLcysTEREpQYLbtwdvbzCZMrcZAgIIvuWWzMcXvv02M4BkSNm3j7iVKyndtavNdoPBgHdgIOlpaVmu5R0czN+DBnH4s88yt+3/4ANu+vFHyl1xPRERV8k1hIwbNy7z8w8//NClxYiIiJRUflWrUmX6dE4OHow5Ph7vMmWoPHUqPmXLZh5jSkiwe64pPh6ApL17ubRtG6E33EBgrVpU6d+fw1f97faLiMCvenWbAAJgSkpix8sv037TJic/MxGRrHINIQ2v6IfapEkTV9YiIiJSopV57DFKPfggqQcP4h8VhVdgoO3+7t059sorWK64u+EVGkqpO+4geuhQjo0fb91oMHDN6NHUHjsWU1ISJ2bNwpySQnjz5jT47DPObd5s9/qxmzfbdMMWEXEVhwamA5w9e5b9+/eTlJRks/3WW291WlEiIiIllVdwsM1A9Cv5V6/OtQsWcOSFF0iNjibwuuu4ZvJkUo4cuRxAACwWDo8eTYWHHqLB1KnU++ADTElJ+P03vX7K+fN22w+uVUsBRETcwqEQsmTJEiZNmoS/vz/+/v6Z2w0Gg0KIiIiIG5S++25K3XUX5qQkvIODATg5fXrWAy0WEv74g6CoKLyDgvAOCsrcVb5DB8q0bs2FDRtsTqk7cqRLaxcRyeBQCJkzZw4jR46kbdu2rqpHREREcmEwGDIDCEBINt2ls9tuMBi46ccf2Td+fObsWDWffppK3bq5oFoRkawcCiFGo5E2bdq4qhYRERHJh7Dmzak4YACnrhhsXmXIEILr1cv2HN/QUK576y2ue+std5QoImLDoRDSrl071q9fz8033+yqekRERCQf6kybRsV+/bj099+EtmhB6PXXe7okEZFsORRCBg0axFNPPcWSJUsoU6aMzb6hQ4c6tTARERFxTFiLFoS1aOHpMkREcuXQMueTJk3iwoUL+Pv7YzKZbD5ERERERETywqE7Ib/++iszZ84kMjLSVfWIiIiIiEgx59CdkLCwMMpesXKriIiIlEzGhAQST5zwdBkiUkQ5FEJ69erFjBkz1P1KRESkhLJYLGx++WXmVqjAvCpV+L5pUy7u2ePpskSkiHGoO9bcuXO5cOECixYtIjw83GbfN99849TCREREpPA58OWX7Hj//czH5//+m9U9e9Jtyxb+GD2aA/Pm4e3vT/2BA2k0eLBWYBcRuxwKIY8//rir6hAREZEi4PCiRVm2Xdi+nbWPP86+uXMzt2144QUwGGg8eLD7ihORIsOhENK5c2dX1SEiIiJFQGD58lm2efn6Er1gQZbtuz79VCFEROxyaEzIyZMns/0QERGR4q/+s8/iExRks632Y49hTk/PcqwpJcVdZYlIEePQnZBHHnkEg8GAxWIBsOnn+csvvzi3MhERESl0yjRowF0bN7Jz4kSSY2Kodvfd1H3iCRIvXODgVXdDoh56yENVikhh51AI+d///mfz+Ny5c8yZM4fbb7/dqUWJiIhI4VW2USPazJhhs63dZ58BcOj77/Hy9aXuY4/RfPRoD1QnIkWBQyHk6kUKIyMjGTZsGK+88gq33nqrUwsTERGRosO/dGk6ffst6cnJGLy98fbz83RJIlKIORRC7AkJCeHUqVPOqEVERESKOJ/AQE+XICJFgEMh5K+//rJ5nJKSwo8//kiNGjWcWpSIiIiIiBRfDoWQl156yeZxYGAgderU4eWXX3ZqUSIiIlI8mNPTiV6yhEsnTlCza1fCr7nG0yWJSCHgUAhZvXq13e1ms9kpxYiIiEjxYUpLY0HHjpz4/XcAfnvxRe5evJgaWndMpMRzaJ0Qe4xGI7fddpszahEREZFi5OCSJZkBBMBkNLJ+xAgPViQihUWBQwiQuW6IiIiISIZLJ05k3Xb8uAcqEZHCxikh5MpFC0VEREQAanbtipevr822qHvv9VA1IlKYOCWEiIiIiFwtvEYN7v7+e8o3aUJg+fI0GjiQtu+/7+myRKQQyNPA9Kun5r1Senq604oRERGR4qXmnXdS8847PV2GiBQyeQohV0/NezV1xxIRERERkbzKUwjJbmpeERERERERR2lMiIiIiIiIuFWuIeTXX3/NU0N5PU5EREREREq2XLtjLV26lLlz59K1a1eaN29OxYoVM/fFxMSwZcsWli5dSqlSpWjbtq1LixURERERkaIv1xDy3nvvsX79er799lsmTpyIn58fwcHBJCYmYjQaadCgAY8++ig333yzO+rNtG3bNoYMGUJAQEDmtpCQEL799lu31iEiIiIiIo7J08D0m266iZtuuom4uDj27t1LQkICoaGh1K5dm1KlSrm4xJwtXboUb29vj9YgIiIiIiJ5l6cQkiE8PJwWLVq4qhYREREp5lLj4zm4dCn+YWHUuPNOTfMvUkI5FEK2b99ud7ufnx8RERGUKVPGKUU5olevXqSnp3PNNdfw6KOP0qRJE7fXICIiIrlLjY/n6xtv5PyuXQBc17cvnb74wsNViYgnOBRChgwZgsFgwGKxZG7LeGwwGLj++usZPnx4gcLIuHHjWLlyZbb7GzduzIQJE6hWrRrTp0+nRo0apKamsmTJEl555RWmTJlCVFRUvq8vIiIirnHk558zAwjAzlmzaDdhAv7h4R6sSkQ8waEQMnz4cNauXUu/fv2IjIwkJiaGL774gtatWxMVFcWUKVOYPHkyo0aNyndBzz//PE8++WS2+zPGf5QpUyYz7AQFBdGjRw82btzImjVrcgwhw4cPx8/PD4BOnTrRqVOnfNcqIiIieRdw1ThS74AAvK+YYCbDypUrM9+QNBqN7ihNRNzMsGbNGkvuh1k9/PDDTJs2jZCQkMxtly5dYuDAgcydO5ezZ88yaNAgvvvuO5cUm5sXXniBevXq8cQTT2TZl5iYSNeuXYmLiyMsLMwD1YmIiMjq555j2yef4BMQQKdZs6jz4IM5Hh8fH094eDhLly4lODjYTVWKiKs5tGJ6fHy83QFkcXFxAJQrV47k5GTnVJaLLVu2cOrUKcxmMykpKSxYsIB///2XNm3auOX6IiIi4rgOkybxXFISzyQk5BpARKT4cqg7VuPGjRk7diwDBw4kIiKC06dPM3369MzB4AcOHKBcuXKuqDOLPXv28MEHHxAfH4+fnx81a9Zk3Lhx1KlTxy3XFxERkfzxsdMFS0RKFoe6Y8XGxvLWW2+xbdu2zDsijRs3ZuTIkZQpU4ZDhw6RmJhIgwYNXFZwfqk7loiISNGj7lgixZNDd0JKly7Nhx9+yNmzZzl37hzlypWjfPnymftr1Kjh9AJFRERERKR4cWhMSAYfHx+8vLzw8XEow4iIiIiIiDh2JyQxMZFx48axfv16wLpGSOvWrRk6dKjNjFkiIiIiIiLZcehOyLRp00hISOCzzz5j2bJlTJs2jcTERD777DNX1SciIiIiIsWMQyFk8+bNvPbaa0RFRREYGEhUVBQjR45k06ZNrqpPRERERESKGYdCSGpqapZuVyEhIaSmpjq1KBERESk5LsXE8GW7dnxcrRpbp071dDki4gYOhZC6desyc+ZMzGYzAGazmVmzZmltDhEREcm330aP5uivvxJ/7Bgrn3mGhFOnPF2SiLiYQwPTn3zySV566SV++umnzMUKvb29ef/9911Vn4iIiBR3FssVn1psHotI8eRQCKlevTpz5sxh48aNnDlzhgoVKtCqVSuCgoJcVZ+IiIgUc7e8/jrndu8mNjqa1q++SmilSp4uSURcLNcQMnPmzGz3HT58mMOHDwPw+OOPO60oERERKTlCK1Xi0d9+83QZIuJGuYaQHTt25NqIwWBwSjEiIiIiIlL85RpCPvroI3fUISIiIiIiJYRDs2OJiIiIiIgUlEKIiIiIiIi4lUKIiIiIiIi4lUKIiIiIiIi4lUKIiIiIiIi4lUKIiIiIiIi4lUKIiIiIiIi4lUKIiIiIiIi4lUKIiIiIiIi4lUKIiIiIFErpqal89/DDni5DRFzAx9MFiIjkVTrbMfE33jTAhxs8XY6IuNjBn39m39Klni5DRFxAIURECj0LFpJ4GiOzMrf5cj/BzMKgG7oixVaFhg0JCA+HuDhPlyIiTqa/3iJS6KWzyiaAAKSxgDQWeaQeEXGPUtWr88Qff3i6DBFxAYUQESn00vgtm+1r3VuIMQbifoe0c+69rkgJFhIR4ekSRMQFFEJEpNDzpmo226u7pwCLBQ4+B39WhX/bwB+V4cgI91xbRESkGFIIEZFCz4+eeHGNzTYDFfHjUfcUcPpzOPUxWNKtjy1GOD4Wzn3rnuuLiIgUMwohIlLoGQgjlNX4Mxgf2uLP04TxK16Ud08BZ//n2HYRERHJkWbHEpEiwYtIghjroatbHNwuIiIiOdGdEBGR3JTrkc32nu6tQ0REpJhQCBERyU3kQIgYQOavTIMPVHoByiuEiIiI5Ie6Y4mI5MbgBVHToOoISN4LgdeBfyX31mBOB4M3GAzuva6IiIgLKISIiOSVfzXrhzvFH4Hfn4Mjy8A7AOr1hdbvgU+ge+sQERFxIoUQEZHCymyCH26DuP3Wx+mJsGMypF2CjrNce+2Dm+D8IajSBCrWc+21RESkxFEIEREprI6uvBxArrTvf3DThxBQxvnXjD8NU7vB4c2XtzV9APp+Cb7+zr+eiIiUSBqYLiJSWCWftr/dnAapsa655twBtgEE4K9vYcXbrrmeiIiUSAohIiKFVZWO1sHoVwuPgrCazr9e/BnYsdT+vg0znX89EREpsRRCREQKq9Bq0PKqOxA+gdB2qmtmyUqJA4vZ/r4kF915ERGREkljQkRECrNmw+CaO+HQYvAJhmt7QnBF11yrXC0oew2cP5x1X92OrrmmiIiUSAohIiKFXdmG1g9X8/KCe9+HGT2sM3NlCAiDuzUmREREnEchRERELmt6H5ReD2s/hnOHoOr10GEwVIjydGUiIlKMKISIiIitGi2tH4XJz/Phy/fg4E6odA08+Czc/6SnqxIRkXzSwHQRESncvvsUhveA3X9CajIc2g3vPgWfvOrpyvIsee1azjz4IMdr1+bk9ddzccwYTOfOebosERGPUQgREZHCKz0NPhttf9+8j+Bi4X8hH/vaa5xu356kb78lff9+jNu2cXHkSE42aULavn2eLk9ExCMUQkREpPA6vAcuZLNoozEV/tno3noclLxqFXFvvWV3n+nECc727OnmikRECgeFEBERKbxCSuW8P6y0W8rIr/hPPslxv/Hvv0nZWLiDlIiIKyiEiIhI4RVZFW5ob39f1ShofJN763GQcetWpxwjIlLcKISIiEjhNnIGVKphu61UOXh7nmtWjnciQ1BQ7scEBrqhEhGRwkVT9IqISOFWuQZ8uxtWfwfR/0LlmnB7TwgK8XRluQrq3p34d9/N/gBfX4K6dnVfQSIihYRCiIiIFH5+/tD5IU9X4bCwZ57h0vTpmC9csLs/dMAAvCMi3FyViIjnqTuWiIiIi/hUqULEzz/jE3XVivM+PoQ++SRlJkzwSF0iIp6mOyEiIiIu5N+0KZX37SPlp58wbtuGITiYoO7d8alc2dOliYh4jEKIiIiIixkMBgI7dSKwUydPlyIiUigU2hBy4MABpk+fzv79+4mNjeX999+nWbNmNsckJCQwceJENm3aBECrVq0YPHgwISGFf7CiiIiIiEhJVWjHhPj6+nLLLbfwzjvvZHvMmDFjiI2NZe7cucydO5fY2NgcjxcREREREc8rtHdCqlevTvXq1bPdHxMTw+bNm5k+fTrh4eEAPPnkkzzxxBOcPn2aCM02IiIiIiJSKBXaOyG5OXDgAL6+vkRdMeNIVFQUvr6+HDhwwIOViYiIiIhITtx+J2TcuHGsXLky2/2NGzdmQh6mLExKSrI79iMkJISkpKSClCgiIiIiIi7k9hDy/PPP8+STT2a739vbO0/tBAUFcenSpSzbL126RFBQULbnDR8+HD8/PwA6depEJ81UIiIiUqisXLky8w1Lo9Ho4WpExBXcHkICAwMJDAwscDtRUVGkpaURHR1NrVq1AIiOjiYtLc2mi9bVxo4dS1hYWIGvLyIiIq5x5ZuE8fHxTJ482cMViYizFdoxIRaLBaPRmPkOSHp6OkajEZPJBEBkZCQtW7bk008/JS4ujri4OD799FNat26tQekiIiIiIoVYoQ0hp0+ftnknZNiwYXTq1Ikvv/wy85jhw4cTFhbGww8/zMMPP0x4eDivvvqqp0oWEREREZE8KLRT9EZGRrJmzZocjwkLC2PUqFFuqkhERERERJyh0N4JERERERGR4kkhRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3EohRERERERE3MrH0wVk58CBA0yfPp39+/cTGxvL+++/T7NmzWyOGTx4MDt37sTH5/LTGDhwIPfcc4+bqxURERERkbwqtHdCfH19ueWWW3jnnXdyPK5nz56sWLEi86OkBZCVK1d6ugTJA32dCj99jQo/fY0KP32NRCSvCm0IqV69Ol27dqVOnTqeLqVQ0y/8okFfp8JPX6PCT1+jwk9fIxHJq0LbHSuvfvjhBxYtWkTp0qW5+eab6d27N4GBgZ4uS0REREREsuH2EDJu3Lgc3ylp3LgxEyZMyFNb/fv3p1q1aoSEhHDw4EHeffddTp06xeuvv57lWIvFAkB8fHy+6i6sjEZjsXtOxZG+ToWfvkaFn75GhZ8rvkYZ7WX8HReR4sGwZs0at/5UJycnYzQas93v7e1NSEiIzbb27dvbHZh+tW3btvHiiy+yfPly/P39bfadPXuWBx98MP+Fi4iIiMfMnz+f8uXLe7oMEXESt98JCQwMdFl3KYPBANh/t6Rs2bLMnz+fwMDAzONERESkcLNYLCQnJ1O2bFlPlyIiTlRox4RYLBbS0tIyH6enp2M0GvH29sbb25sLFy5w4MABGjZsSEBAAIcPH2bKlCncdNNNBAQEZGnPy8tL76CIiIgUQVf3kBCRoq/QhpDTp0/Tq1evzMfDhg0DoE+fPvTt2xej0cgXX3zBsWPHMJlMlClThltuuYVHH33UUyWLiIiIiEgeuH1MiIiIiIiIlGyFdp0QEREREREpngptd6zibPXq1SxatIjo6GiSkpJYtWoV3t7eAOzatYuvvvqKPXv2kJqaSkREBA888ABdunSxaeP777/nm2++4eLFi1SrVo2nn36axo0bZ+4/ffo0EyZMYPv27fj6+tKhQweeeuopfH1989xGSZbT1wggOjqaSZMmsW/fPoKDg+natSt9+vSxmRxh1qxZLFu2jMTERGrXrs3gwYOpUaOGU9sQWxcuXGDy5Mn8/fffGI1GqlevzhNPPEGTJk0A6wx6U6ZM4ejRo5QuXZqePXvSrVu3zPONRiNTpkxhzZo1pKWl0bhxY4YMGUKFChUyj8mtDcndzp07mTFjBnv27MHLy4trrrmGSZMm4eXl5ZafLcm7UaNGsW7dOpsZKvVzJCLOoDshHhASEkK3bt14+umns+yLj4+nTZs2zJgxg6VLl/Lss8/yySefsG7dusxj1q5dy4wZMxg2bBhLliyhS5cuDBs2jDNnzgBgNpsZPnw4oaGhfPvtt0ybNo1//vmHTz/9NM9tlHQ5fY2SkpJ45ZVXaNCgAYsWLWL8+PEsW7aMBQsWZB7zzTffsGLFCsaPH8+iRYto0KABr7zyCsnJyU5rQ7KaMGECZ8+eZebMmSxevJi2bdsyfPhw4uPjiYmJ4dVXX6VLly4sWbKEoUOHMn36dH7//ffM86dMmcKOHTuYNm0a3377LaGhoYwYMQKz2QyQpzYkZzt37mTo0KH8v707j4nqev84/h4G2QWhpSKxUCgIWqUY08YYJdZqAQ1S3ButVjQ2LlGitolVu2ijYqzVGLVuQa22cYllWhe0pi7RGKvGKhJFsbiLVbQwDPvA9w/DTacjYn+lo7/280pI4Jxznzn3DGdmnrnn3puYmMiOHTuwWCxMmDABk8nksrklT2bv3r1UVlY6lGkeiUhzURLyFLz++uu8+eabhIaGOtV17dqVpKQkAgMDMZlMdO7cmc6dO3P69GmjjcViITk5mfj4eFq0aEFaWhpt27YlJycHgLNnz3L16lUmTpyIr68vISEhjB49mt27dxv3aGkqxn/d456jw4cPU1dXR3p6Op6enkRGRjJ06FCys7ONNhaLhSFDhhAZGYmnpyfp6enU1NQYb7LNEUOc3bx5k4SEBFq1aoXZbCYlJYWKigpu3LjB3r17adu2LWlpabRo0YL4+HiSk5P57rvvgIff3ubk5DB69GhCQkLw9fVl4sSJFBYWcu7cOYAmY0jTVq1aRd++fUlMTMTLywuz2UyHDh0wmUwum1vStIZkfvr06Q7lmkci0lyUhDzjbDYb58+fJzo62igrKCggNjbWoV1MTAwFBQVGfWhoKAEBAUZ9bGwslZWVXL9+/YliSOMuX75MVFSUw/Ks2NhYbt26hc1mo6ysjKKiItq3b2/Um81moqOjuXTpUrPFEGfvvPMOR48epbi4mNraWrKzswkNDeXll19u8n/+2rVrVFVVOYx5QEAAbdq0McZc8+bvqaysJC8vDzc3N8aPH09qairjxo3j0KFDgGvmljStvr6ehQsXMmLECFq3bu1Qp3kkIs1F54Q8w2pqapgzZw5hYWH06dPHKC8vL3e6ZnrLli25ffu2Ue/r6+tU31D3JDGkcTab7ZFjBw/HteFmmX9u4+fnZ4x/c8QQZx07duTHH39k0KBBuLm54e/vz5w5c/D09MRms9G2bVuH9i1btjQ+mDaMa1PP2+NiyONZrVbq6urYt28f8+bNIzo6mqNHjzJ37lyef/55l8ytP782ijOLxUJ9fT0pKSlOdZpHItJcdCTkGVVZWcnMmTOpqalh3rx5Dt/q+fj4UFZW5tDearUab64+Pj5OL+ZWq9Woe5IY0jhfX99Hjh08HNeGMf5zm7KyMqOuOWKIo7q6OqZNm0ZQUBAWi4V9+/Yxbdo0ZsyYQUFBQaNj/sd5A/+3503z5sl4e3sDkJiYSGxsLGazmYSEBOLj4zly5IhL5pY83s2bN/n666+dlmE10DwSkeaiJOQZZLVamT59OmazmQULFhhv3A2ioqK4cOGCQ9nFixeJiooy6m/fvk1JSYlRn5+fj5eXFy+++OITxZDGNSztsdvtRll+fj6hoaH4+vri5+dHSEiIw/ja7XYKCgqMZXXNEUMcWa1Wbt26xYABA/D398dsNtO9e3dCQ0M5ceIEUVFR5OfnO2yTn59v/M+HhYXh6enpMOYlJSUUFRUZY95UDHk8Pz8/QkNDG71KlSvmljxebm4upaWlvP/++6SmphpXrPrkk09YtGiR5pGINBslIU+B3W6nurqampoa4OGJfNXV1dTV1XH//n0yMjIIDg5m7ty5eHh4OG2fmprKnj17OHv2LDU1NVgsFq5fv05SUhIAcXFxhIWFsXLlSsrLy7lz5w5ZWVkkJycb8ZqK8V/3uOcoISEBNzc3srKyqKqqorCwkK1btzpcXjI1NZWtW7dSWFhIVVUVWVlZuLu706NHD4BmiSGOAgICCA8PJzs7G5vNRl1dHceOHePKlSu0a9eOxMRErl27hsVioaamhrNnz7Jnzx7efvttADw8PEhKSiIrK4s7d+5QXl7OihUrCA8Pp2PHjgBNxpCmpaWlkZOTQ0FBAXV1dRw9epQzZ87Qo0cPl80taVzPnj355ptvWLt2rfEDMHXqVMaNG6d5JCLNRndMfwpycnLIzMx0Kv/yyy85c+YM69evx8vLy6EuLi7OYZuGe3w8ePCA8PBwJkyYYNwLAR5eArHhPiEeHh706tWL8ePHOyQ1TcX4L3vccxQfH8/ly5dZunQpFy9exMfHh/79+zvdyyArK4udO3dSXl5OTEwMU6ZMITIy0ojVHDHE0Y0bN/jqq6/Iy8ujurqa4OBgBg4caKxt/+WXX1i+fLnDvQn++MHnj/c3qK6u5tVXX2Xq1KlO9zd4XAxp2ubNm7FYLJSVldG2bVtGjhxJ9+7dgeaZF03FkL/mjTfecLpPiOaRiPxdSkJERERERMSltBxLRERERERcSkmIiIiIiIi4lJIQERERERFxKSUhIiIiIiLiUkpCRERERETEpZSEiIiIiIiISykJERERERERl1ISIiIiIiIiLqUkRETkH7Bp0yaWLFnyxO3v3r3L4MGDKSkp+ec6JSIi8oxwf9odEJGnIzk52fi9traWuro6PDw8jLLMzEzi4uKeRteaXUZGBp06dWLMmDEuebz79++zZcsWNm7c2GgfqquryczM5MKFC8yfP5+wsDB69uzJhg0bmDx5skv6KSIi8rQoCRH5j9qzZ4/x+7p168jNzf1L39w/C2pqamjRooXLHq+2thZ396ZfNi0WC126dCEwMPCR9aWlpcyaNQu73c7y5ctp1aoVAP369WP8+PGkp6fj5+fXnF0XERF5pigJEREn1dXVbNiwgYMHD1JaWkpkZCSTJk0iOjoagPXr13Pq1Ck6d+7Mrl27qKysJCUlheHDh7N48WJOnDiBn58fkyZNonv37gDk5OSwbt06hg4dypYtW6iurqZbt25MnjwZb29vAMrKylizZg3Hjx+noqKC9u3bM3nyZEJDQwFYsGABlZWVBAQEcOjQIWJiYsjMzGTRokWcPHmSkpISWrVqxVtvvcWoUaNwc3Nj8eLF5ObmkpeXx/bt24GHCVjDPixbtszY7wULFmC325k5cyYAw4YNIzExkfPnz3Pu3DlGjBhBYmIiixYtIj8/n8rKStq0acO7775Lz549jTiHDx9m0KBBjxzbmzdvMmPGDCIiIvjoo4/w9PQ06l566SX8/f35+eef6dWrVzM9myIiIs8enRMiIk4WL17MhQsXWLJkCdnZ2fTs2ZMPP/yQsrIyo8358+fx9/dny5YtfPHFF2zfvp3p06eTlpbG999/z4ABA1i4cCGVlZXGNvfv3+fy5cts3LiRtWvXUlhYyPLlywGor69n9uzZ2Gw2Vq9ezbZt24wP6rW1tUaMI0eOEBMTw/bt2/nss88AaN++PStWrGD37t3Mnj2b7Oxsdu7cCcDUqVPp1KkTw4YNY8+ePQ5HgJ7EDz/8wKhRo9i1axcDBgygrq6O5ORkNm/ebOzn559/TmFhIQBVVVVcuXKFiIgIp1gXL15k0qRJdO3alU8//dQhAWkQGRlJfn7+X+qjiIjI/zdKQkTEQUlJCXv37mXKlCkEBwdjNptJS0vD19eXY8eOGe1at27NoEGDcHd3JzY2loiICNq1a0dcXBxubm4kJiZitVq5ceOGQ/yJEyfi7e1NcHAwo0ePZu/evdjtdi5dusS5c+eYNm0a/v7+eHh4MHbsWG7fvs358+eN7WNiYujbty/u7u54eXkBD5cxBQUFYTKZ6NChA7179+bkyZPNMh7Jycm88sormEwmvLy8CA4OJiEhAW9vb9zd3enXrx/h4eGcPn0aAKvVCvDI5VR5eXnU1NSQnJyMyWR65OP5+vpSWlraLH0XERF5Vmk5log4uHnzJgATJkxwKK+treXu3bvG30FBQQ71Xl5ePPfcc8bfDUusKioqjDI/Pz+HD+dt2rShtraWBw8ecOPGDex2O4MHD3bq02+//Wb8HhIS4lBXX1/P5s2b2b9/P8XFxdTX11NdXU1sbOwT7/Pj/PnxrFYrq1at4tSpU5SWlmIymaioqOD3338HoGXLlgAOR40apKamYrPZyMjIYP78+XTo0MGpjc1mIywsrFn6LiIi8qxSEiIiDhqSi3Xr1tG6detmjV1WVkZZWZmRiBQVFeHu7k5gYCBBQUG4u7uTnZ392JO/3dwcD+D+9NNPbN++nczMTKKiojCbzSxbtoyLFy82ug2Aj4+Pw1IxgOLiYuMk8ca2Xb16NdeuXWPJkiW88MILmEwmxowZQ319PQCenp6Eh4dz5coVpyTDzc2NjIwM/P39mTZtGnPmzOG1115zaFNYWEifPn0a3X8REZF/Ay3HEhEHISEhdO/enaVLl1JUVARAeXk5x48fp7i4+G/HX7lyJRUVFdy7d4+srCz69OmD2WymU6dORERE8OWXX/LgwQPg4VGHQ4cOOSULf2Sz2TCbzQQGBmIymTh9+jT79+93aBMYGMj169cdytq1a8evv/5Kbm4udrudgwcPcubMmSb7b7PZ8PT0xN/fn9raWnbs2MGVK1cc2iQkJHDixIlGY6SnpzN27FhmzpzJgQMHjPKrV69SUlLilJiIiIj82+hIiIg4mTVrFt9++y0ffPABxcXF+Pj40KFDh799/4qgoCAiIiIYOXIkVVVVdOvWjUmTJgFgNptZtGgRWVlZTJgwgZKSElq2bElcXBxdu3ZtNGZSUhK5ubmkp6fj5uZGly5d6N27t8ORkCFDhrBw4UJSUlKor69n586dxMfHM3z4cD7++GPsdju9evUiISGhyX0YM2YMCxcuZMCAAfj6+pKcnEzHjh0d2qSmpvLee+9x//59p2VrDQYOHIi/vz/z58/HarXSv39/du7cSVJSkrGkS0RE5N/KdODAgfqn3QkR+fdruETvtm3bnnZXXGLTpk3cu3ePjIyMJ2p/9+5dJk6cyJo1awgICPhnOyciIvKU6UiIiMg/YMSIEX+pfXBwMFu3bv2HeiMiIvJs0TkhIiIiIiLiUlqOJSIiIiIiLqUjISIiIiIi4lJKQkRERERExKWUhIiIiIiIiEspCREREREREZdSEiIiIiIiIi6lJERERERERFzqf1b5339snmGNAAAAAElFTkSuQmCC\n", + "text/plain": [ + "<Figure size 640x480 with 1 Axes>" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "from matplotlib import cm\n", + "import matplotlib.style\n", + "import matplotlib as mpl\n", + "mpl.style.use('classic') # seleccion del estilo del gráfico\n", + "fig, ax = plt.subplots()\n", + "plt.scatter(temp, np.log(lum), s=10*rad, c=temp_color, norm=matplotlib.colors.LogNorm(), edgecolors='none')\n", + "ax.set_xlim(max(temp), min(temp)) # cambiar el orden del eje temp(X)\n", + "plt.title(\"Diagrama Hertzsprung-Russell\")\n", + "plt.xlabel(\"Temperatura(K)\")\n", + "plt.ylabel(\"log(Luminosidad(Lsun))\")\n", + "plt.text(0.5, 0.6, r'$\\sum_{i=0}^\\infty x_i$', fontsize=20)\n", + "plt.text(0.5, 0.5, 'matplotlib', horizontalalignment='center',verticalalignment='center', transform=ax.transAxes)\n", + "\n", + "ax.grid(False)\n", + "plt.show()\n" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": {}, + "outputs": [], + "source": [ + "c=np.log" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + } + ], + "metadata": { + "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.7.3" + } + }, + "nbformat": 4, + "nbformat_minor": 4 +}