{
"cells": [
{
"cell_type": "markdown",
"id": "3d54b332",
"metadata": {},
"source": [
"# 6. Pathway Enrichment of LR loadings\n",
"\n",
"Cell-cell communication inference from single-cell data provides unprecedented insights into the molecular mechanisms underlying cell-cell communication. In particular, it enables the capture of CCC interactions at a systems level, which is not possible with traditional approaches. However, as the number of inferred interactions increases, the interpretation of the inferred cell-cell communication networks becomes more challenging. To this end, as done in other omics studies, we can perform pathway enrichment analysis to identify the general biological processes that are enriched in the inferred interactions. Pathway enrichment thus serves two purposes; it reduces the dimensionality of the inferred interactions, and also provides a biological summary of the inferred interactions.\n",
"\n",
"In this tutorial, we will show how to use the `cell2cell` and `liana` packages to perform classical gene set enrichment analysis using `GSEA` with `KEGG Pathways`, and also we will use perform a footprint enrichment analysis using a `multi-variate linear model` with the PROGENy pathway resource. We will use the [`GSEAPY`](https://gseapy.readthedocs.io/en/latest/) and [`decoupler-py`](https://decoupler-py.readthedocs.io/en/latest) packages to perform these analyses."
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "2ac22be8",
"metadata": {},
"outputs": [],
"source": [
"import cell2cell as c2c\n",
"import liana as li\n",
"\n",
"import pandas as pd\n",
"\n",
"import decoupler as dc\n",
"\n",
"import matplotlib.pyplot as plt\n",
"import plotnine as p9\n",
"import seaborn as sns\n",
"\n",
"%matplotlib inline"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "907c4925",
"metadata": {},
"outputs": [],
"source": [
"import logging\n",
"logging.disable(logging.WARNING) # This is to avoid too many logs because of GSEApy\n",
"\n",
"import warnings\n",
"warnings.filterwarnings('ignore')"
]
},
{
"cell_type": "markdown",
"id": "eded8fe1",
"metadata": {},
"source": [
"## Directories"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "c0a73062",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"../../data/tc2c-outputs/ already exists.\n"
]
}
],
"source": [
"output_folder = '../../data/tc2c-outputs/'\n",
"c2c.io.directories.create_directory(output_folder)"
]
},
{
"cell_type": "markdown",
"id": "2f5d0d38",
"metadata": {},
"source": [
"## Load Data"
]
},
{
"cell_type": "markdown",
"id": "973d1b0b",
"metadata": {},
"source": [
"**Open the loadings obtained from the tensor factorization**"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "1657732b",
"metadata": {},
"outputs": [],
"source": [
"factors = c2c.io.load_tensor_factors(output_folder + 'Loadings.xlsx')"
]
},
{
"cell_type": "markdown",
"id": "24b18655",
"metadata": {},
"source": [
"**Load list of ligand-receptor pairs used as reference to [run LIANA](./02-Infer-Communication-Scores.ipynb)**"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "e96e13f7",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
\n",
"\n",
"
\n",
" \n",
"
\n",
"
\n",
"
ligand
\n",
"
receptor
\n",
"
\n",
" \n",
" \n",
"
\n",
"
0
\n",
"
LGALS9
\n",
"
PTPRC
\n",
"
\n",
"
\n",
"
1
\n",
"
LGALS9
\n",
"
MET
\n",
"
\n",
"
\n",
"
2
\n",
"
LGALS9
\n",
"
CD44
\n",
"
\n",
"
\n",
"
3
\n",
"
LGALS9
\n",
"
LRP1
\n",
"
\n",
"
\n",
"
4
\n",
"
LGALS9
\n",
"
CD47
\n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" ligand receptor\n",
"0 LGALS9 PTPRC\n",
"1 LGALS9 MET\n",
"2 LGALS9 CD44\n",
"3 LGALS9 LRP1\n",
"4 LGALS9 CD47"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"lr_pairs = li.resource.select_resource('consensus')\n",
"lr_pairs.head()"
]
},
{
"cell_type": "markdown",
"id": "43090ea6",
"metadata": {},
"source": [
"**Generate a list with the LR pair names**\n",
"\n",
"KEGG was originally designed to work with list of genes. Here we are using ligand-receptor pairs, so instead of using Gene Sets we need to build our LR-gene sets. This list is a first step to build the LR-gene sets."
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "0550bdd3",
"metadata": {},
"outputs": [],
"source": [
"lr_list = ['^'.join(row) for idx, row in lr_pairs.iterrows()] "
]
},
{
"cell_type": "markdown",
"id": "3572fb9d",
"metadata": {},
"source": [
"## Gene Set Enrichment Analysis\n",
"\n",
"In this case, we will run the PreRank method of GSEA since our LR pairs are pre-ranked based on their loadings in each factor obtained by Tensor-cell2cell."
]
},
{
"cell_type": "markdown",
"id": "fa6385ac",
"metadata": {},
"source": [
"### Generate LR-gene set for running GSEA\n",
"\n",
"`cell2cell` includes some resources associated with the Gene sets employed in GSEA. It includes annotations of GO Terms (BP), KEGG, and Reactome, for human and mouse."
]
},
{
"cell_type": "markdown",
"id": "97f8d4a6",
"metadata": {},
"source": [
"**First, specify which organism and annotation DB to use**"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "c11ec2fd",
"metadata": {},
"outputs": [],
"source": [
"organism = 'human' # For the COVID-19 data analyzed with LIANA+Tensor-cell2cell\n",
"pathwaydb = 'KEGG' # KEGG Pathways"
]
},
{
"cell_type": "markdown",
"id": "77835f05",
"metadata": {},
"source": [
"**Generate the LR-gene set that will be used for running GSEA**"
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "2b9e896f",
"metadata": {},
"outputs": [],
"source": [
"lr_set = c2c.external.generate_lr_geneset(lr_list,\n",
" complex_sep='_',\n",
" lr_sep='^',\n",
" organism=organism,\n",
" pathwaydb=pathwaydb,\n",
" readable_name=True,\n",
" output_folder=output_folder\n",
" )"
]
},
{
"cell_type": "markdown",
"id": "b6ab0812",
"metadata": {},
"source": [
"### Extract LR loadings for each factor"
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "f43cd70c",
"metadata": {},
"outputs": [],
"source": [
"lr_loadings = factors['Ligand-Receptor Pairs']"
]
},
{
"cell_type": "markdown",
"id": "0a8b4d7d",
"metadata": {},
"source": [
"### Run GSEA\n",
"\n",
"This external function implemented in `cell2cell` builds upon the function `gseapy.prerank()`.\n",
"\n",
"All outputs will be saved in the `output_folder`. Make sure that this path exist before running this function."
]
},
{
"cell_type": "code",
"execution_count": 10,
"id": "40220386",
"metadata": {
"scrolled": true
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"../../data/tc2c-outputs/ already exists.\n"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"100%|███████████████████████████████████████████| 10/10 [00:08<00:00, 1.11it/s]"
]
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"CPU times: user 9.86 s, sys: 132 ms, total: 9.99 s\n",
"Wall time: 9.05 s\n"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"\n"
]
}
],
"source": [
"%%time\n",
"pvals, scores, gsea_df = c2c.external.run_gsea(loadings=lr_loadings, \n",
" lr_set=lr_set,\n",
" output_folder=output_folder,\n",
" weight=1,\n",
" min_size=15,\n",
" permutations=999,\n",
" processes=6,\n",
" random_state=6,\n",
" significance_threshold=0.05,\n",
" )"
]
},
{
"cell_type": "markdown",
"id": "c4654706",
"metadata": {},
"source": [
"### Visualization\n",
"\n",
"We can use the `pvals` and `scores`outputs to visualize the results as a Dot Plot:"
]
},
{
"cell_type": "code",
"execution_count": 11,
"id": "587f2ede",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA1wAAAQaCAYAAABXQ3k7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/OQEPoAAAACXBIWXMAAA9hAAAPYQGoP6dpAADkf0lEQVR4nOzdd3QU1cPG8WcTktBCCgkdpWgCCAEJvUiRJl1AkW6jKtgVeO29N4oIiCA/qdISkN4EpIMUQZBO6C1AEiBl5/0j7pIlm7o7hPL9nONxd+bOnTuzs2GenTt3LIZhGAIAAAAAuJ1HTjcAAAAAAO5UBC4AAAAAMAmBCwAAAABMQuACAAAAAJMQuAAAAADAJAQuAAAAADAJgQsAAAAATELgAgAAAACTELgAAAAAwCQELgAAoJkzZyo0NFShoaGKiorK6ebcltauXavQ0FDVqVNHV65cyenm5Jj169fbj6X169fndHPSdPLkSVWsWFEVK1bUwYMHc7QtPXr0UGhoqHr06JGj7YA5cuV0AwAAALIiNjZWu3bt0vbt27V9+3bt2LFDx44dkyQVL15cy5Yty1J9e/fu1f/+9z/9+eefOn36tPLmzasyZcqoTZs2euyxx5QrV8anS1arVR9//LEk6emnn1aePHlSlWncuLG9nSnlypVL+fPnV5kyZVSnTh09/vjjKly4cJa2AVlXpEgRdejQQVOnTtVnn32mUaNG5XSTcIcicAEAgNtKv379tGHDBrfUNW3aNL3//vtKSEiwT7t27Zo2b96szZs3a+bMmfrxxx8VGBiYbj3z5s3T3r17FRAQoK5du2apDYmJiYqOjtaWLVu0ZcsWjR8/Xp9++qmaNm2arW1C5vXt21czZ87U8uXLtX37doWFheV0k3AHInABAIDblr+/vypWrKgtW7YoLi4uS8uuXLlS77zzjqxWq4KCgtSvXz9VrlxZ0dHRmj59uhYtWqTt27fr+eef18SJE+Xp6ZlmXbarI507d1bevHnTXW+hQoX0008/2d8nJibq2LFjmjFjhpYvX66YmBi99NJLmjZtmipUqJClbULWFC9eXM2aNdO8efM0cuTIHLvKNXHixBxZL24O7uECAAC3ldatW+urr77SokWLtH79ev30008KCAjIUh0JCQn64IMPZLValT9/fk2ePFk9evRQWFiYHnroIQ0bNsx+pWrz5s2aM2dOmnWtWbNG+/btkyS1bds2w3V7eXkpJCTE/l+FChXUtGlTjRo1Sk899ZS9fT/88EOWtgnZ07p1a0nJAfzo0aM53BrciQhcAADgttK5c2e1bt1a9957b7brWLx4sf3kum/fvrrnnntSlXn99dfl5+cnSQ5XpG7022+/SZIeeOABlS1bNtttkqSBAwcqd+7ckqTVq1fLarW6VB8yVr9+ffn7+8tqtWrGjBk53RzcgehSCAAAMu38+fOaMGGCVq5cqaioKF27dk3BwcEKDw9X586dVa1atQzr2LRpkyZMmKAtW7bo0qVLKlSokOrWratnnnlG9957r3r06KENGzaoRo0apnW1Wrp0qf31o48+6rRMnjx51KJFC02dOlX79u3TwYMHVbp0aYcy165dsw/S0axZM5fblS9fPpUtW1Z///234uLiFB0dneH9Y5K0ceNGde/eXZL0wQcf6PHHH0+3/I8//qivv/5aUvL9Z/fdd5993tGjR7V48WJt2LBBe/fu1dmzZyVJBQsWVOXKldWhQwc99NBD2d3ETH++w4YN0/DhwyVJe/bsSbPc5cuXNWnSJC1fvlyHDh1STEyMvatp+/bt1bx5c1ksljSX9/LyUuPGjTVz5kz9/vvvevHFF7O9bZcuXdKvv/6qFStW6MCBA4qLi5Ovr68CAwNVunRp1a1bV82aNVNQUJDDcuntk9DQ0EyvP71BY/7++29NnTpV69ev1+nTp2UYhgoXLqxatWrpySefTHVsu2O7kIzABQAAMmX16tV64YUXFBMT4zD92LFjOnbsmCIiItStWze9+eab8vBw3olm9OjR+vrrr2UYhn1aVFSUpk6dqrlz5+r77783dRtsNm/eLEkqXbq0goOD0yxXo0YNTZ06VZK0ZcuWVCel27Zt09WrVyVJlStXdkvbvLy87K8ze4WrWrVqKlasmI4fP67IyMgMA9fcuXMlSeXLl08Vtpo0aeJ0mePHj+v48eOaP3++2rZtq08++SRTIziaae3atXrxxRcVHR3tMP3MmTNavny5li9frgYNGuibb75Rvnz50qyncuXKmjlzpg4fPqxDhw6pVKlSWW7L/v379eSTT+r06dMO0y9cuKALFy5o//79WrJkiaxWqz0c3wxWq1WfffaZJkyY4PC9k6RDhw7p0KFD+u233/T222+rc+fOqZa/VbfrdkLgAgAAGdq9e7f69eunhIQEeXl5qVu3bmrcuLHy5Mmj3bt3a/To0YqKitKvv/6qPHny6LXXXktVx++//66vvvpKUvJgF88++6z9itimTZs0ZswYvfzyy1m+HyurYmNjdeLECUlSmTJl0i2bcv7+/ftTzd+0aZMkyWKxqGLFii63LTExUQcOHJCUHLz8/f0ztZzFYlGrVq00ZswYbdq0SadOnUpzaPl//vlHe/fulSS1adPGYZ7VapWXl5fq1aununXr6r777pOfn58uXryogwcPatKkSfr3338VERGhkiVLatCgQdnfWBdt3rxZvXv3VkJCgoKCgtS9e3eVK1dOhQoV0unTp/X7778rIiJCK1eu1ODBgzVs2LA060o5OuHGjRuzFbhee+01nT59Wl5eXnrsscf00EMPKSgoSIZh6OTJk/rrr7+0ZMmSLNcbGRmZ7vzNmzfr3XfflSQVK1Ys1fwPPvhAkyZNkiRVr15djz76qEqWLKncuXNrz549mjBhgv7991+9/fbbCgoK0sMPP3xTtutuQuACAAAZeuutt5SQkCBPT0+NGjVK9erVs88LCwtTixYt1LVrV+3bt0/jxo1T+/btdf/999vLxMfH66OPPpIkBQQEaOrUqQ73YD344INq0qSJnnjiCR06dMjUbTl58qT9dZEiRdItm3J+yuVstm7dKkkqWbKkfH19XW7b1KlTdenSJUlSeHh4lq4gtWnTRmPGjJHVatXcuXP1zDPPOC1nu7rl4eFhHzDCJjg4WMuWLVOhQoVSLVe7dm116dJFQ4cO1cyZM/Xzzz/rqaeecst2Z1VCQoJee+01JSQkqH79+ho2bJjDs88eeOABNWrUSNWrV9dbb72lRYsWac2aNapbt67T+kJCQuTl5aWEhATt3LlTjz32WJbac/ToUf3999+SpMGDB6e60hMWFqZmzZrptddes3++mRUSEpLmvKioKPtV4aCgIH3xxRcO89esWWMPWx9++GGq7QoLC1Pbtm3Vp08frVu3Th999JEaNGhgP+7M3K67CYNmAACAdNkeLixJjz32mEPYsvHz89P7778vKfkqie0kz2bJkiX2e4EGDhzodMCL0qVL67nnnnN381OJjY21v85oCPeU850NO28LYQULFsx2exITE3X48GF988039lAqKc3AlJbQ0FD7yXlaV0UMw9C8efMkJV/tuPEqWN68eZ2GLRuLxaI33nhDnp6eiouL059//pmlNrrLvHnzdOzYMfn4+Ojzzz93+qBpSXr88cftV69mzpyZZn25cuWyD5ASFRWV5facOXPG/jq9+xgtFot9Pa6KjY1V//79df78eXl7e2vEiBEqWrSoQ5nRo0dLkpo3b55miPTx8dHbb78tKbl78Pr16+3zcmK77kQELgAAkK6UJ9WdOnVKs1x4eLh9lL4bT8Rt7z08PFJ1Y0upbdu26Q5w4A7x8fH21ynvl3LG29vb/tp2r1ZKFy5ckCQVKFAg0+s/duyYQkND7f898MADatasmUaNGqWkpCRZLBa9+OKL2RqYwrZvd+/e7bQL5ObNm3X8+HGHsulJSEjQyZMntX//fu3du1d79+7V6dOn7V0d//nnnyy30R1sA0NUr149w0FFbEHhr7/+SrecLTCkDBmZlfI+wFmzZmV5+ayyWq165ZVX7F1DP/jgA1WpUsWhTExMjP0B4c2bN0+3vrJly9q78tqu2ko3f7vuVHQpBADgDnXq1CldvHjR6Tw/P7807/G50b///ispOZyUL18+3bKVK1fW/v37dfjwYcXHx9sDi62OkiVLphtO/P39VbJkSR05ciRTbcuOlCEqISEh3bIpw5ltuPaUbIM1ZCVwpSV//vyqVauWnnrqKYerCQkJCTp48GCay5UuXdoeHFu3bm0flCQiIkIvvfSSQ1nblS9vb+80T8ITEhI0bdo0zZkzR7t27Up3H9kC5822c+dOSckDuWR2FD/bFda02D7DK1euZLk9JUuWVLVq1bRp0yaNHz9eq1evVrNmzVSjRg1VqVIlzStw2fXVV19p+fLlkqRnn31W7du3T1Vm165d9kFXXn75Zb388suZqjvlfrrZ23WnInABAHCH+uabb9L8VfrRRx/Vp59+mql6bKHC398/w3uKbMNCG4ahS5cu2d/bgl9mhjgPDAw0NXClHK3OWTfBlFLOd9b90MfHRwkJCbp27Vqm11+oUCGH53p5enrK19dXwcHBTq/unTp1Kt2rUUuXLlWJEiUkJQ+aUK1aNW3cuFFz5851CFwJCQlasGCBJKlhw4ZOQ2J0dLSefvpp+307GcnKdrvT+fPns7yMsyuUKdm2JbsjL3799dd64YUXtHXrVu3bt0/79u3TyJEj5eXlpcqVK6t169bq0KGDfHx8slW/zezZszV27FhJUqNGjfTKK684LXfu3Lls1X/jfrpZ23UnI3ABAIC7Ssore84GwkgpowE2AgMDFRMTk2pY8vR4eXmlOxCCq9q0aaONGzcqKipKW7du1YMPPigp+WqQrZ1pBbiPPvrIHraaNGmijh07KjQ0VAULFpSPj489EDZs2FAnTpxINcz4zZKUlCRJeuihh5yOiJkdth8Fsnu1snDhwpoyZYrWrl2rRYsWaePGjdq3b58SEhK0adMmbdq0SePGjdPo0aPTfeZVerZu3aq33npLUvJgGl9++WWaj2BI+UiB999/334cZOTGe7Fuxnbd6QhcAADcoT799NNMX8VKj+1+nejoaCUmJqZ7BcDWHclisTicuNpO4jJzZSI7Vy+yIn/+/CpatKhOnDhhH4I9LSnn2+5PS8l2Nc7MEdpKlCiR7oN/b9SiRQt98MEHSkhIUGRkpP1E2zY6oa+vrxo2bJhquZiYGM2fP19SciD78ssv01xHWl1VM8MW2jJ6xlh6Vx/9/f11+vRpJSQkuC282rbpxoEnsqp27dqqXbu2pOQul2vXrtXUqVO1bt06HTlyRC+99JJmz56d5XqPHz+u559/XvHx8QoICNDIkSOVP3/+NMunfKRA7ty5Xd5PZm3X3YBBMwAAQLpsw7snJCRo9+7d6Zbdvn27JOnee+91uFfK9nDdo0ePpnuyHh0draNHj7ra5AyFh4dLkg4ePJjuIAm2QQckqWrVqqnm205ijxw5kumHFJvNz8/PPuDG/PnzlZiYqLi4OC1dulRS8gAKKT8bm0OHDtnv12rZsmWa9e/fvz/DrpjpsXXpzCikpvd4gAoVKkhKvpcr5X122XX27Fn7A71TPgjaVQEBAWrZsqUmTJigxo0bS0oe0CSrjz6Ii4vTgAEDdPbsWXl5eWnYsGEqWbJkusuUL1/eHm63bNmSrfanxV3bdbcgcAEAgHTVqVPH/nrGjBlplrPd43HjMpLsv4zbnhGVloiIiJvSTS3lw13Tus/typUr9nue7rvvPqfdpWzBLS4uzumogDnF1mXw/PnzWrNmjZYsWWIfDCKt7oS2bnpS+gNHTJkyxaW22e43O3jwoD3k3Oj8+fPpDjlvO8m/fPlyusO9Z5btsQdS8sAvZrB9B6SsDTZiGIZef/11+48db7/9tqpXr57hcoGBgfaRC+fOnWvalePsbtfdhMAFAADSFRYWpooVK0qSpk+frrVr16Yqc/nyZfuzfDw8PNSlSxeH+U2bNrU/q2rYsGFOB8U4dOiQRowY4e7mO9W0aVP7FYIff/zRaXs+//xz+9W4tJ6JlXI0QdvVvVtB48aN7d3NIiMj7SG3cOHCqlGjhtNl7rnnHvsVkVmzZjkNvsuWLdOvv/7qUtts609ISND//ve/VPMTEhL05ptvpjvIxaOPPmrv+vfZZ59p48aN6a5z06ZNDlcrb2T77Hx8fDIVZm60e/fudK/+GoZhD5AWi0XFixfPdN3ffvutFi9eLEnq0aOHHn/88Uwv279/f0nJ3UUHDRqU7lXF+Ph4/frrrw4DoZi5XXcT7uECAAAZ+vDDD/XYY48pISFBvXv3Vo8ePdSoUSPlyZNHu3fv1ujRo+1dAZ9++ulU94v4+Pho6NCheuWVV3ThwgU99thj6t27tz2wbNy4UWPHjpXValWpUqXS7Zp0+PBhbd682WGa7WHGsbGxqa541K9f3+F5QlLywBVvvfWW+vXrp5iYGHXp0kX9+/dXWFiYLl68qOnTp2vhwoWSkq9itWvXzmlbSpQoodDQUO3Zs0dr165Vx44dM9iTN4ePj4+aNWummTNnaunSpfaugq1atUpzkIWAgAA1aNBAK1as0KpVq/T000+rS5cuKlasmM6dO6dFixZp1qxZKlmypC5dupTtKyYNGjRQ8eLFdezYMX333Xe6cOGCmjZtKh8fH+3bt08TJ07Url27VKVKlTSfneXt7a1vv/1WPXr0UFxcnHr16qWWLVuqSZMmKlGihKxWq86cOaO///5bixcv1t69e/XWW2+lGTZtPyLUq1fP6fD/Gdm9e7eGDBmiSpUqqVGjRnrggQcUFBSkxMRERUVFaebMmVqzZo2k5DCc3sOlb2zXqFGjJCV37e3UqZP92VvOeHl5OVyJbdCggXr27KlffvlFGzduVMuWLfXEE08oPDxc/v7+iouL05EjR7Rp0yYtXrxYFy9eVPv27e0jDpq1XXcbAhcAAMhQ+fLlNWrUKL3wwguKiYnRuHHjNG7cuFTlunXrluYw1a1bt9bRo0f13XffKTo6Wl988YXD/Dx58ui7777T6NGjdejQoTSHmd68ebOGDBnidF50dHSqeb/88kuqwCUln4y+9957ev/993X27Fl98MEHqcqEhYVp+PDh8vT0dLo+SXr88cf1wQcfaOnSpbpy5cot82yiNm3aaObMmQ73W7Vt2zbdZd5991117dpVx48f159//pmqW1+xYsU0YsQI9enTJ9vt8vb21ueff67evXsrLi5O48eP1/jx4+3zPT09NXToUF28eDHdhxVXqVJFEydO1IsvvqgTJ04oMjLS/pwxZ9IaYCIqKsq+noz2T0Z27Njh0D3xRg8++KA++uijTNd34sQJ++t///03zeBvU7x4cftDoW2GDh0qPz8//fDDDzpz5oyGDRuW5vJ58+Z1eqy7e7vuNgQuAACQKfXq1dPixYs1YcIErVy5UkePHlV8fLyCgoIUHh6uJ554wqGLnTP9+/dX9erV9fPPP2vr1q26fPmygoODVatWLT3zzDMqW7asvv76a0nJo+mZ7fHHH7efuK9du1anT59Wnjx5VLZsWbVp00aPPfZYhs9latu2rb744gvFxcVp2bJlatWqlentzoxatWopODjYPijIfffdl+GDq4sWLaqZM2dqzJgxWrp0qY4fPy4fHx8VL15cTZo0Uc+ePVMNG54d1apV04wZMzRq1CitXbtWFy5ckL+/v6pWraonn3xSVatWTTcY2FSpUkWLFi3SzJkztXz5cu3atUsXLlyQh4eHAgMDVbZsWVWvXl3NmjVTmTJlnNYxb948GYahQoUKOdzblxWtW7dWwYIF9eeff2rHjh06deqUzp07p8TERBUsWFAVKlRQy5Yt073CaBaLxaLnn39e7dq105QpU7Ru3TpFRUXp8uXLyp07t4oWLary5curXr16atKkicMVvlt5u24nFiOnHqAAAABwg4SEBFWrVk1Xr15V//799eKLL+Z0kzLl3Xff1eTJk1WnTh39/PPPOd0cZJLVatUjjzyiQ4cO6ZVXXnHpyh2QFqIoAAC4ZSxZssQ+WIJthLXbwXPPPae8efPqzz//TLcbHG4tv//+uw4dOqSAgAB169Ytp5uDOxSBCwAA3DSHDx9Oc15UVJT9Qc1BQUGqV6/ezWqWy4KDg9WrVy9JumkjLcI1hmHYB6QYOHCg/flggLtxDxcAALhpHnnkET300ENq1KiR7rvvPuXNm1fnzp3T+vXrNWXKFPuw1a+//nqG907dap599ll7m2+lwTPg3OnTp9W8eXO1bt1aTzzxRE43B3cw7uECAAA3TWhoaLrzPTw89OKLL6pv3743qUUAYC4CFwAAuGmWL1+uP/74Q1u3btXZs2cVHR0tb29v+wN5u3XrluoZXgBwOyNwAQAAAIBJGDQDAAAAAExC4AIAAAAAkxC4AAAAAMAkBC4AAAAAMAmBCwAAAABMQuACAAAAAJMQuAAAAADAJAQuAAAAADAJgQsAAAAATELgAgAAAACTELgAAAAAwCQELgAAAAAwCYELAAAAAExC4AIAAAAAkxC4AAAAAMAkBC4AAAAAMAmBCwAAAABMQuACAAAAAJMQuAAAAADAJAQuAAAAADAJgQsAAAAATELgAgAAyMD69esVGhqqBQsW5HRTJJnTnmHDhik0NDRTZUNDQzVs2DC3rRu4k+XK6QYAAHA7iI+P13fffac5c+bo0qVLCg0N1Ysvvqi6detmuOypU6f08ccfa82aNbJarapZs6aGDh2qkiVL3oSWIy2ZDRe//PKLyS25u8TGxuqnn37Stm3btGPHDl28eFGffPKJOnTokKnlN27cqJ9++km7d+/W+fPnVaBAAZUrV04DBgxQeHh4qvJbtmzRF198oV27dil//vx65JFH9NJLLylfvnxuaxNuPVk9Tpxx199uAhcAAJkwePBgLVy4UD179lSpUqU0a9Ys9enTRxMmTFC1atXSXC42NlY9e/bU5cuX1bdvX3l5eWn8+PHq3r27Zs+erYCAgJu4FUjp888/d3g/Z84crVmzJtX0smXLav/+/TezaXe0CxcuaMSIESpWrJhCQ0O1YcOGLC1/6NAheXh46IknnlBQUJAuXbqkiIgIde/eXT/++KMeeughe9ndu3frySefVNmyZTV48GCdPHlS48aN06FDhzR27Fi3tQm3nqwcJ8649W+3AQAA0rVt2zYjJCTEGDt2rH3a1atXjSZNmhidO3dOd9nRo0cbISEhxrZt2+zT9u3bZ5QvX9746quvTGszsu69994zQkJCnM5bt26dERISYsyfPz9bdcfFxbnSNLe3x5nvv/8+ze2/UUhIiPH9999naz3Xrl0zTp8+bRiGYWzfvt0ICQkxZsyYka26bOLi4ow6deoYTz/9tMP0Z5991qhbt65x+fJl+7Rp06YZISEhxqpVq0xtE1I7d+6ccerUqRxbf1rHiTPu/NvNPVwAAGRgwYIF8vT0VOfOne3TfHx81KlTJ23dulUnTpxIc9mFCxeqUqVKCgsLs08rW7asateurfnz55vabrif1WrVDz/8oIceekiVKlVSr169dPjwYYcyPXr0UOvWrbVz505169ZNlStX1tdffy0puWvq999/r6ZNm6pixYpq0KCBPv/8c8XHxzvUsWbNGnXp0kXVqlXTgw8+qObNm9vryGp7JGn+/Pnq0KGDwsLCVLNmTb366qs6depUhtsbHx+vjz/+WLVq1dKDDz6ofv366eTJk07L7t+/X8ePH8+wTm9vbwUHB2dYLivy5MmjwMBAXb582T4tJiZGf/75p9q2bav8+fPbp7dr10558+Z1+P6Z0Sak9u+//6pRo0bq37+/li5dqsTExJu6fmfHSVrc+bebLoUAAGRg9+7dKlWqlMNJmyT7P8S7d+9W0aJFUy1ntVq1Z88edezYMdW8SpUqafXq1YqJiUlVL25dY8aMkcVi0dNPP62YmBiNHTtWr776qqZPn+5QLjo6Wr1791arVq3Utm1bFSxYUFarVf3799fmzZv1+OOPq2zZstq7d68mTJigQ4cOaeTIkZKST0r79u2r0NBQDRo0SN7e3jp8+LC2bNmSrfbMnDlTQ4YMUaVKlfTyyy/r3Llz+uWXX7RlyxbNnj1bBQoUSHN7/+///k8RERFq3bq1qlatqnXr1qlPnz5Oy7Zs2VI1atTQxIkTs7NrsywmJkbx8fG6cOGC5syZo71796pfv372+Xv27FFiYqIqVqzosJy3t7fKly+v3bt335R24rry5curf//+mjVrlgYMGKDg4GA9+uij6tixo0qVKuV0mYSEhEwFJEny9/eXh4fj9aSMjhNn3P23m8AFAEAGzpw54/TXb9u006dPO10uOjpa8fHxGS5L4Lp9XLt2TbNnz5a3t7ckqUCBAvroo4+0d+9ehYSE2MudOXNG7733np544gn7tDlz5ujPP//UxIkTHe77u//++/XOO+9oy5Ytqlq1qtasWaOEhASNGTNGgYGBLrUnISFBX375pUJCQvTrr7/Kx8dHkhQeHq6+fftq/PjxGjRokNO6//nnH0VERKhr16565513JEndunXTK6+8oj179mRj77nXCy+8oNWrV0uSvLy81LlzZw0YMMA+/8yZM5KkQoUKpVo2ODhYmzdvvjkNhV2BAgX0/PPP67nnntP69es1Y8YM/fLLLxo9erSqV6+uTp06qUWLFsqdO7d9mS1btqhnz56Zqn/p0qUqUaKEw7SMjhNn3P23m8AFAEAGrl69aj+hTcl28nr16lWny127dk2S0l3WVga3hw4dOjh8nrbgdPToUYfA5e3tnWqEuwULFqhs2bIqU6aMzp8/b59eq1YtSclDvVetWtV+xWnp0qXq2LFjql/ss9KenTt36ty5c3r++eftx5wkNWzYUGXKlNGKFSvSDFwrV66UlNxFMqVevXpp7ty5qcrf7BD26quv6umnn9aJEyc0e/ZsJSQkKDExMdX3Mq3vX1rfW5jPYrGoVq1aqlWrlt555x3NnTtXM2fO1BtvvKEPP/xQrVq10ssvvyw/Pz+VK1dOP//8c6bqdRaQMjpOnHH3324CFwAAGcidO3eqe2yk6//gpvw1NiXbP8zpLZveP/q49RQrVszhvS0cXbp0yWF64cKFU52sHT58WPv371ft2rWd1n3u3DlJyV3zpk+frjfffFNfffWVateuraZNm6pFixapwldG7bHdU1W6dOlU6ytTpky6V3mOHTsmDw8P3XPPPamWuxWUL1/e/rpt27bq0KGDhgwZou+//17S9e9lWt+/tL63uLny58+vJ554Qo8++qh++OEHjRo1SlOmTNETTzwhPz8/+fn5qU6dOtmuP6PjxBl3/+0mcAEAkIHg4GCnAwyk12VJSr6fwNvb214uK8vi1pTW1SbDMBzeOzuZt1qtCgkJ0ZAhQ5zWUaRIEfuyv/76q9avX68VK1Zo1apV+v333zV16lSNGzdOnp6eWW7Pnc7b21uNGzfW6NGjdfXqVeXOnTvdLr9nzpzhu3eL2L59u2bMmKHff/9dly5dUuXKldWpUyeVLVtWUnLouXjxYqbqCgwMdPh+3MjZceKMu/92E7gAAMhAuXLltH79+lQ3SW/btk2S4y+oKXl4eNi7dd1o+/btKlmyJPdv3UXuuece/fPPP6pdu7YsFku6ZT08PFS7dm3Vrl1bQ4YM0ahRo/TNN99o/fr1Wfq133YF7ODBg6murB08eDDVFbKUihcvLqvVqiNHjjhc1Tpw4ECm138zXb16VYZhKDY2Vrlz51ZISIhy5cqlnTt3qmXLlvZy8fHx2r17tx555JEcbO3d7dy5c5ozZ45mzpypf//9V/7+/nr00UfVqVMnh665krR161aX7uG60Y3HiTPu/tvNsPAAAGSgRYsWSkpK0tSpU+3T4uPjNXPmTFWuXNk+QuHx48dTPSC3efPm2rFjh3bs2GGfduDAAa1bt04tWrS4ORuAW8IjjzyiU6dOadq0aanmXb16VXFxcZKSb9i/kS3UO+vilJ6KFSuqYMGCmjJlisOyK1eu1P79+9WwYcM0l7U9GPbGUQcnTJjgtHxmh4XPitOnT2v//v1KSEiwT7N1vUzp0qVLWrRokYoWLaqCBQtKknx9fVW7dm1FREQoJibGXnbOnDmKi4vj+5cDTpw4oQEDBuihhx7S559/ruDgYH3zzTdatWqVhg4dmipsSbLfw5WZ/1Lew5XZ40Qy/283V7gAAMhA5cqV1aJFC3399dc6d+6c7r33Xs2aNUvHjh3TRx99ZC/3xhtvaMOGDQ6DB3Tt2lXTp09X37599fTTTytXrlwaP368ChYsqKeffjonNgc5pF27dpo/f77eeecd+wAZSUlJOnDggBYsWKCxY8eqUqVKGjFihDZt2qQGDRqoePHiOnfunCZNmqQiRYooPDw8S+v08vLSq6++qiFDhqh79+5q1aqVfVj44sWL68knn0xz2fLly6t169aaNGmSLl++rAcffFDr1q1z+pwvKWvDwv/vf//TpUuX7N39li9fbn++V48ePeTr6ytJ+vrrrzVr1iyHKxe9e/dW4cKFVblyZRUsWFDHjx/XzJkzdfr0aX3zzTcO63nppZf0xBNPqEePHnr88cd18uRJ/fzzz6pXr549UGa1Tci+I0eOaNeuXerTp486duyY4dUoSdm+hysrx4nZf7sJXAAAZMLnn3+ub7/9VhEREbp48aJCQ0M1atQoVa9ePd3l8ufPr4kTJ+rjjz/WDz/8IKvVqpo1a2rIkCEZDvmNO4uHh4dGjBih8ePHa86cOVq8eLHy5MmjEiVKqEePHvaBLRo3bqxjx45pxowZunDhggICAlSjRg0NHDgwWyf9HTp0UO7cuTVmzBh9+eWXyps3r5o0aaLXXnst3WdwSdLHH3+sgIAARUZGaunSpapZs6ZGjx6tBg0aZGsf2IwbN07Hjh2zv1+0aJEWLVokKXlgg/S2s2PHjpo3b57Gjx+vy5cvq0CBAqpcubK++uorh+H2JemBBx7Qzz//rC+//FKffPKJ8uXLp06dOunll192a5uQOZUrV9ayZcvSHXnTXbJynDjjzr/dFuNuu6sSAAAAAG4S7uECAAAAAJMQuAAAAADAJAQuAAAAADAJgQsAAAAATELgAgAAAACTELgAAAAAwCQELgAAAAAwCYELAAAAAExC4AIAAAAAkxC4AAAAAMAkBC4AAAAAMAmBCwAAAABMQuACAAAAAJMQuAAAAADAJAQuAAAAADBJrpxuAAAAAABkx4EDB7R27Vpt3rxZGzdu1MmTJxUfHy9vb28VKVJE1atXV3h4uGrXrq0yZcrkSBsthmEYObJmAAAAAMii+Ph4zZ49W8OGDdPq1aslSV5eXkpISEhVNuX0evXqaeDAgWrfvr28vb1vWnsJXAAAAABuCxEREerTp49OnTolT09PJSUlZXpZW/nChQtr9OjRatu2rYktvY57uAAAAADc0s6fP6/u3burXbt2OnPmjCRlKWylLH/mzBm1a9dO3bt31/nz593e1htxhQsAAADALWv//v1q2LChTpw4keWQlR5PT08VLVpUK1asUNmyZd1W740IXAAAAABuSfv27VOdOnV04cIFJSYmur3+XLlyKSAgQGvXrjUtdBG4AAAAANxyzp8/r8qVK+vkyZOmhC2bXLlyqWjRovrrr78UGBjo9vq5hwsAAADALWfQoEE6ceKEqWFLkhITE3X8+HENGjTIlPq5wgUAAADglhIREaF27drd9PXOmTPH7aMXErgAAACgPXv2aOPGjbp8+bJ8fX1Vo0YNhYSE5HSzcBeKj4/XPffcozNnzshqtd609Xp4eCg4OFhHjx6Vl5eX2+rN5baaAAAAcFsxDEOzZ8/Wd999p5UrV6aa37BhQw0aNEjt27eXxWLJgRbibjR79mydOnXqpq/XarXq1KlTmj17th577DG31csVLgAAgLtQQkKCnnnmGU2cODHNB8japvfo0UM//fSTW3/1B9JSv359rV271q1DwGeWp6en6tSpoz/++MNtdRK4AAAA7jKGYahXr1763//+p8ycClosFnXv3l0TJkzgShdMdeDAAVOfiZVZ+/fvV5kyZdxSF6MUAgAA3GVmz56tiRMnZipsSckBbeLEiZo9e7a5DcNdb+3atTndBEnSunXr3FYXgQsAAOAu8/3338vT0zNLy3h6emrYsGEmtQhItnnz5hzvuurl5aXNmze7rT4CFwAAwF1kz549WrFiRZbvj0lKStLy5cu1d+9ek1oGSBs3blRCQkKOtiEhIUEbN250W30ELgAA3Ix7XHArc/VEcsOGDW5qCZDayZMnc7oJktzbDgIXAABuZAtbhC7cqi5fvuzS8Xn58mU3tgZwFB8fn9NNkCRdu3bNbXURuAAAcCPbIAQMAoxbla+vr0vHp6+vrxtbAzjy9vbO6SZIknx8fNxWF4ELAAA3I2zhVlajRo0cXR5IT5EiRXK6CZLc2w4CFwAAwF0kJCREDRs2zNYohY0aNVJISIhJLQOk6tWr3xKjFFavXt1t9RG4AAAA7jKDBg3K1iiFAwcONKlFQLLw8PBbYpTC8PBwt9VH4AIAALjLtG/fXj169Mj04BkWi0U9e/ZU+/btzW0Y7nq1a9fO6SZIkmrVquW2ughcAAAAdxmLxaKffvpJ3bt3l6Q0uxfapvfo0UNjx45l9E2YrkyZMqpXr16Wu7y6i6enp+rXr68yZcq4rU4CFwAAwF3Iy8tLEyZM0MyZM1W/fn2nZR566CHNnDlT48ePz/H7anD3GDhwYJa7vLqLGV1nLQZDKQEAANz19u7dqw0bNujy5cvy9fVVjRo1GCADOSI+Pl733HOPzpw5I6vVetPW6+HhoeDgYB09etStPzAQuAAAAADcUiIiItSuXbscWW+bNm3cWiddCgEAAADcUtq2batu3brdtHu5PD091a1bN7eHLYkrXAAAAABuQefPn1flypV18uRJJSYmmrYeT09PFStWTH/99ZcCAwPdXj9XuAAAAADccgIDA7VixQoFBAQoV65cpqzD09NTgYGBWr58uSlhSyJwAQAAALhFlS1bVmvXrlXRokXd3r3QdmVr7dq1Klu2rFvrTonABQAAAOCWVbZsWf3111964oknJCWPJugK2/JPPPGE/vrrL1PDlsQ9XAAAAABuExEREerTp49OnTolT0/PLD2vy1a+cOHCGj16tNq2bWtiS6/jChcAAACA20Lbtm119OhRTZs2TXXq1LFPT+u5WSmn16lTR9OmTdPRo0dvWtiSuMIFAAAA4DZ14MABrVu3Tps3b9bGjRt18uRJXbt2TT4+PipSpIiqV6+u8PBw1apVS2XKlMmRNhK4AAAAAMAkdCkEAAAAAJMQuAAAAADAJAQuAAAAADAJgQsAAAAATELgAgAAAACTELgAAAAAwCQELgAAAAAwCYELAAAAAExC4AIAAAAAkxC4AAAAAMAkBC4AAAAAMAmBCwAAAABMQuACAAAAAJMQuAAAAADAJAQuAAAAADAJgQsAAAAATELgAgAAAACTELgAAAAAwCQELgAAAAAwCYELAAAAAExC4AIAAAAAkxC4AAAAAMAkBC4AAAAAMAmBCwAAAABMQuACAAAAAJPkyukGAAAAIHPWr1+vnj172t8/8sgj+vbbb9NdZvDgwZo1a5Ykac+ePQ7zhg0bpuHDh2epDSNGjFCTJk3SbN/s2bO1bds2nTx5UlevXlXu3LlVqFAhlSpVSmFhYapTp47CwsLk4cHv/rg7ELgAAABuUwsWLNCePXsUGhqao+2IjY3V66+/riVLljidd/DgQR08eFDLly/Xd999pzFjxuihhx7KgZYCNx+BCwAA4DZlGEa2rlI58/HHH6tSpUoZlitWrFiqaYMGDdLq1aslSffee68ee+wxVapUSX5+foqLi9Phw4e1ZcsWLVu2TOfOnXO5rcDthMAFAABwGwoICNCFCxe0ePFi7dq1SxUqVHCpvhIlSigkJCTLy61YscIeturVq6cffvhB3t7eDmXCw8PVoUMHJSUlacmSJSpatKhLbQVuJ3SeBQAAuA316NHDHmy+//77HGvH0qVL7a8HDx6cKmyl5OnpqebNm+v++++/GU0DbgkELgAAgNtQ0aJF1blzZ0nS8uXLtX379hxpx/Hjx+2v77333hxpA3ArI3ABAADcpvr06aPcuXNLkr777rscaYOXl5f99f79+3OkDcCtjMAFAABwmypUqJC6dOkiSVq9erU2bdp009vwwAMP2F+///77On/+/E1vA3ArY9AMAACA21jv3r01depUxcXF6fvvv9cvv/ySrXqioqIUEBCQbhkvLy+VLl3aYVqnTp30008/6cqVK9qyZYsaNWqkBg0aqEaNGgoLC1O5cuXSva8LuNMRuAAAAG5jBQsWVLdu3TRmzBitX79e69atU61atbJcz9ChQzMsU7x4cS1btsxhWtGiRfXNN9/o5ZdfVlxcnK5evaqFCxdq4cKFkiRvb2+FhYWpadOmat++vfz9/bPcNuB2RpdCAACA29wzzzyjfPnyScqZe7kaNWqk+fPnq0ePHqmuksXHx2vTpk365JNP1LRpU82ePfumtw/ISVzhAgAAuM0FBASoV69eGjlypLZs2aJVq1apfv36Warjl19+Uc2aNbPdhiJFiujNN9/U0KFDtXv3bv3111/6+++/tXnzZh06dEiSdOnSJb3xxhtKSkpSx44ds70u4HbCFS4AAIA7wFNPPaUCBQpIkoYNG5Zj7fDw8NADDzygbt266eOPP9bChQs1Y8YMhYeH28t89tlniomJybE2AjcTgQsAAOAOUKBAAT355JOSpG3btmn58uU526AUKlasqLFjx9qf03Xx4kWtXbs2h1sF3BwELgAAgDtEr1697INS5ORVLmfy5s2rVq1a2d8fPnw4B1sD3DwELgAAgDtE/vz59cwzz0iS/v77by1evDiHW+SoUKFC9tcWiyUHWwLcPAQuAACAO0j37t1VsGBBSdL3338vwzBMXV9W6t+5c6f9dcmSJc1oDnDLIXABAADcQfLmzavevXtLkvbu3auVK1eaur533nlHo0aNUnR0dLrl1qxZYx8SPm/evKpTp46p7QJuFQwLDwAAcIfp0qWLfvrpJ505c0YXLlzI1DJRUVGpnqHlTEBAgIKDg+3vo6OjNXXqVA0fPlwNGjRQ9erVFRISIn9/fyUmJurIkSNatmyZ5s+fL6vVKkl64YUXlD9//uxtHHCbIXABAADcYXLnzq1+/frpgw8+yPQyQ4cOzVS5nj176v/+7//s7233ZSUkJGjJkiVasmRJmsv6+Pho0KBB9tEUgbsBgQsAAOAO9Pjjj2vs2LE6ceKEqet588039fTTT2vVqlXatGmT/v33Xx0/flyxsbHKlSuXChQooPvuu081a9ZUu3btVLRoUVPbA9xqLIbZd1ICAAAAwF2KQTMAAAAAwCQELgAAAAAwCYELAAAAAExC4AIAAAAAkxC4AAAAAMAkBC4AAAAAMAmBCwAAAABMQuACAAAAAJMQuAAAAADAJLlyugEAANxuTp8+rV9++UXbtm3Tzp07FRcXp19++UU1a9bM1PLDhg3T8OHDU0339vbWjh073N1cl7m6vZJ06tQpffzxx1qzZo2sVqtq1qypoUOHqmTJkia2PHsuXbqkL774QosXL9bVq1dVqVIlDR48WA888ECGyw4ePFizZs1KNb106dJasGCBGc3NlPj4eH333XeaM2eOLl26pNDQUL344ouqW7duhsveTp+dTXa393b7btrExsbqp59+0rZt27Rjxw5dvHhRn3zyiTp06JCp5V055m82V7Z15syZGjJkiNN5q1evVnBwsLubK4nABQBAlh08eFBjxoxRqVKlFBoaqq1bt2arnnfffVd58+a1v/f09HRXE93K1e2NjY1Vz549dfnyZfXt21deXl4aP368unfvrtmzZysgIMCklmed1WpVnz59tGfPHj3zzDMKCAjQpEmT1KNHD82cOVOlSpXKsA5vb299+OGHDtN8fX1NanHmDB48WAsXLlTPnj1VqlQpzZo1S3369NGECRNUrVq1NJe7nT67lLK7vTa3y3fT5sKFCxoxYoSKFSum0NBQbdiwIdPLuuOYv5lc2VabQYMGqUSJEg7TChQo4K4mpmYAAIAsuXz5snHhwgXDMAxj/vz5RkhIiLFu3bpML//9998bISEhxrlz50xqoXu5ur2jR482QkJCjG3bttmn7du3zyhfvrzx1Vdfubu5Lpk3b54REhJizJ8/3z7t3LlzRrVq1YyXX345w+XfeOMNo0qVKmY2Mcu2bdtmhISEGGPHjrVPu3r1qtGkSROjc+fO6S57O312Nq5s7+323bS5du2acfr0acMwDGP79u1GSEiIMWPGjEwt6+oxf7O5sq0zZswwQkJCjO3bt5vZxFS4hwsAgCzKnz+//P393VJXTEyMDMNwS11mcXV7Fy5cqEqVKiksLMw+rWzZsqpdu7bmz5/vhha6z8KFCxUUFKRmzZrZpwUGBuqRRx7R0qVLFR8fn6l6kpKSFBMTY1Yzs2TBggXy9PRU586d7dN8fHzUqVMnbd26VSdOnEhz2dvps7NxZXtTuh2+mzbe3t7Z7g7nrmP+ZnFlW1OKiYlRUlKSG1qUMQIXAAA55OGHH1Z4eLiqVq2qV199VWfPns3pJrmd1WrVnj17VLFixVTzKlWqpCNHjtwywUSSdu/erQoVKsjDw/EUqVKlSrpy5YoOHjyYYR1XrlxReHi4wsPDVaNGDb333nuKjY01q8kZ2r17t0qVKqX8+fM7TLeFqN27dztd7nb77Gyyu70p3Q3fTRt3HPO3m549eyo8PFyVK1dWv379dOjQIVPXxz1cAADcZAUKFFD37t1VpUoVeXt7a9OmTZo0aZJ27NihGTNmpDpRvJ1FR0crPj7e6S/StmmnT5++Zbb5zJkzTu/xKVSokKTktoaGhqa5fHBwsJ599llVqFBBhmFo1apVmjRpkv755x9NnDhRuXLd/FOvM2fOZLj/nbndPjub7G6vdHd9N21cPeZvJ7lz51aHDh1Us2ZN5c+fXzt37tT48eP1xBNPaNasWSpatKgp6yVwAQDualarVQkJCZkq6+3tLYvF4vI6e/Xq5fC+efPmCgsL06uvvqpJkyapT58+Lq8jLTd7e69du2av60Y+Pj4OZdwtO9t69epVp221Tcuora+88orD+1atWqlUqVL65ptvtHDhQrVq1SqTrXeftLbJtv+vXr3qdLmc/Oxckd3tlXL2u5lTXD3mbyctW7ZUy5Yt7e+bNGmievXqqXv37vrhhx/0/vvvm7JeAhcA4K62ceNG9ezZM1Nlf//9d5UtW9aUdrRp00afffaZ/vzzT1NP6m729tpOcp3dB2I7kbOVcbfsbGvu3LmdttU2LTttffLJJ/Xdd9/pzz//zJHAldY22fZ/7ty5nS6Xk5+dK7K7vWm5Wd/NnGLGMX87qVatmipXrqy1a9eatg4CFwDgrlamTBl98sknmSpr62JjliJFiujixYumruNmb6+/v7+8vb115syZVPNs08zar9nZ1uDgYKdttXVDy05bc+fOLX9/f9M/27QEBwfr1KlTqaZntP9z8rNzRXa3Nz0347uZU8w45m83RYoUMfVeNQIXAOCuFhwcnOmHg5rJMAwdO3ZMFSpUMHU9N3t7PTw8FBISop07d6aat337dpUsWdK0+2Kys63lypXT5s2bZbVaHQYR2L59u/LkyaPSpUtnuR0xMTG6cOGCAgMDs7ysO5QrV07r169XTEyMw77etm2bJKl8+fJOl8vJz84V2d3etNys72ZOMeOYv90cPXrU1GfKMUohAAAmOn78uPbv3+8w7fz586nKTZo0SefPn1f9+vVvVtNM4Wx7mzdvrh07dmjHjh32aQcOHNC6devUokWLm93EdLVo0UJnz57VokWL7NPOnz+vBQsWqFGjRg73uhw5ckRHjhyxv7927ZrTUftGjhwpwzBy7LNt0aKFkpKSNHXqVPu0+Ph4zZw5U5UrV7YPFHC7f3Y2rmzvnfzdlJKvWu3fv9/h3sasHPO3E2fb6uzzXblypf7++29TP1+Lcbs8YAAAgFvIyJEjJUn79u3TvHnz1LFjR5UoUUKSNGDAAHu5Hj16aMOGDdqzZ499WuXKldWyZUuFhITI29tbW7Zs0bx581SuXDlNnjxZefLkubkbkwmubG9MTIweffRRxcbG6umnn1auXLk0fvx4JSUlac6cOTl25ceZpKQkde3aVXv37tUzzzyjgIAATZ48WcePH9dvv/2mMmXK2Ms2btxYkrRs2TJJUlRUlB599FG1atXKXm716tVauXKl6tevr9GjR6caevtmeeGFF7RkyRL16tVL9957r2bNmqUdO3Zo/Pjxql69uqTb/7NLKbvbezt+N23+97//6dKlSzp9+rQmT56sZs2a2a/m9ejRQ76+vho8eLBmzZqlpUuX2r+/WTnmbxXZ3VZbuYoVK8rX11e7du3SjBkzFBwcrN9++01BQUGmtJfABQBANqQ3THLKEzhnJ3Vvvvmm/QGs8fHxKlasmJo1a6Z+/frdkl20JNe2V5JOnjypjz/+WGvWrJHValXNmjU1ZMgQ3Xvvvaa1ObsuXryozz//XEuWLNG1a9dUqVIlvf7666pUqZJDuRsD16VLl/TBBx9o27ZtOn36tJKSknTvvfeqTZs2evrpp+Xl5XXTt8Xm2rVr+vbbbxUZGamLFy8qNDRUL7zwgsOv+nfCZ2eT3e29Hb+bNo0bN9axY8eczrOFDmchRMr8MX+ryO62fvPNN1q5cqWioqJ09epVBQcHq0GDBnr++edNC1sSgQsAAAAATMM9XAAAAABgEgIXAAAAAJiEwAUAAAAAJiFwAQAAAIBJCFwAAAAAYBICFwAAAACYhMAFAAAAACYhcAEA4KL4+HgNGzZM8fHxOd2Um4LtvXPdTdsqsb13ultlewlcAAC4KD4+XsOHD8/xf9RvFrb3znU3bavE9t7pbpXtJXABAAAAgEkIXAAAAABgEgIXAABukC9fvpxuwk3F9t657qZtldjeO92tsL0WwzCMnG4EAAC3istnT8g3qGhONwMAcIcgcAEAcINvOtyna7GX5GGRPCzJ0zw8Ury2SB4WI8Xr5OkWS3I5KbkLiadD+dSvHcpbDKdlbmxD8jSLZJ9mSTHfIst/bzwsye+vtz35tcVD18voehnLDfXY2+hhcVivxUP2emzdZJLX6/nfdE/7aw8PT1k8U0+3WDwdyl+fn1yjxdNxvsWS67/XHo51pyzj6aQei6csHrn+a7Cn9F/55Mb/N93DM3me/bVtusf16RZPe3mLxfP6B5FiukP9N9YpSUq5/lwObTEsKepL8dpQyumW69MdyifXkzzNVqfF/tqwpNwOjxTl067TMJKnG4Yh21liytfWlNOtxn/zJdsJZbrLWZ2XsTopb1jlWI+tfqshq229VkNWw3A+PeVr2RuRdhl7PXKYfr3t/722GvaGWa2O063219Y0Xqcon5RGGcMqI8l5PdeXNRzKO9TjZFkjyZBhtV5fNun6a9v868tdn28kGTJSlLW3K2Xbk1K0N2U9KdaZsn4jybGNhpF0fbv/O0CSp9teJ9lfG1arrNbk8oZhldVWv2GV1VbGVl86dTiuJynV67z5c2tx1AS5Sy631QQAwB3iWuwlxcddznLg8nAIUJI1owDlQuCynR8bHhb9d34sI0XgMizJ75OnS0aKwOVhT1OyV27xuB7ilOK1xcNyPVlZrocv3bBeh7Dx32vDw1PytE3P5bSMUgYv2zTrDfNThqO06jDSmp4iENnKWDyl/06sZKQIR4anZHEy3cNTkvX6svZgY5U9Zlis18vIKlns8eC/+ddP0pMn2157Xi+Tsnyq1/YPIcU0i+N7iyVF2ZTrT+t1yrpcdz0cuet1WkEsG4HLhTJS6mBlZCZwpRGaUgaVtAJaWoErM8s6Xe8Nwcl54Eo9/8bpDoErxWvHwGVNY50ZBK40ApLj6xsDl7PwlbX6nAUuW13uwj1cAAAAAGASAhcAAAAAmITABQAAAAAmIXABAAAAgEkIXAAAAABgEgIXAAAAAJiEwAUAAAAAJiFwAQAAAIBJCFwAAAAAYBICFwAAAACYhMAFAAAAACYhcAEAAACASQhcAAAAAGASAhcAAAAAmITABQAAAAAmIXABAAAAgEkIXAAAAABgklw53QAAAG41PvkKSJI8LMn/SZKHR4rXFsnDYqQqY7Ekl5OSf9H0dCif+rVDeYvhtMyNbUieZpHs0ywp5ltk+e+NhyX5/fW2J7+2eOh6GV0vY7mhHnsbPSwO67V46Ho9ul7G4uH532tPx9eeTqZbHMskV+55fQM9PZPf26ZbcqX4EFKWT/Ha00k9Fk/JI1eK157XG2+b7uGZPM/ZumzTU9WTQf031pm8USnqy5XitYeu71RPx9dKOd3ipPwNr+2/o1uU4tNJ57V7WTKo0pLOap1OtljsW2RIMizXXyvFdNuEVNOdvJblejuNG1/b6rFYZPnvjfHf/y0ehgzrf6+thmQkv7Zar0+3Wg1Z7a8lq0eK1/bpluvlPSwO0+2vDYuMJCfTUy6b5JGivDVFGausSSle/zfdSDJkWK3/LWvImnT9tW3+9eWuzzeSDBkpytrblZRinUnWFNNT1JNinSnrN5Ic22gYSfbtMAxrium210n214bVKqs1ubxhWGW11W9YZbWVsdWXTh2O60lK9Tpv/txyJ4thGEbGxQAAAAAAWUWXQgAAAAAwCYELAAAAAExC4AIAAAAAkxC4AAAAAMAkBC4AAAAAMAmBCwAAAABMQuACAAAAAJMQuAAAAADAJAQuAAAAADAJgQsAAAAATELgAgAAAACTELgAAAAAwCQELgAAAAAwCYELAAAAAExC4AIAAAAAkxC4AAAAAMAkBC4AAAAAMAmBCwAAAABMQuACAAAAAJMQuAAAAADAJAQuAAAAADBJrpxuAAAAt5pz5y7LMLK3rMUiFSzo61IdSBv711zsX3Oxf83njn1sq8NdCFwAANzAMOTyyZA76kDa2L/mYv+ai/1rvltpH9OlEAAAAABMQuACAAAAAJMQuAAAAADAJAQuAAAAADAJgQsAAAAATELgAgAAAACTELgAAAAAwCQELgAAAAAwCYELAAAAAExC4AIAAAAAkxC4AAAAAMAkBC4AAAAAMAmBCwAAAABMQuACAAAAAJMQuAAAAADAJAQuAAAAADAJgQsAAAAATELgAgAAAACT5MrpBgAAcKcwDEPRZy/pWkyc4q4kyK9gAVkslpxu1h3DMAxdvhCnpCuJir1yTfn987J/3cgwDF08d1nxsVcUdyVBBQJ92b9ulHz8XpH1apJi4q4pv38e9q+bGYahS+djlRgXr9gr8fINyHdL7GMCFwAALkhKTNKfCzZr/uQV+nvDXsVcjLXPy++XTw/UCFHLro1Uu3lVeebyzMGW3p6SEq3a+sderZqzXf9ui1Lc5av2eXl9c+v+yiX0ULvKqvLQ/fLMRcedrEpKTNK6xdu0eOpq7d68XzEX4+zz8vvlVfnwsmr2RD3VbFKZ4zcbkhKt2rb6gNZE/q39O44r7vI1+7y8vj4qW6mY6rV5QGH1ynD8ZlNSYpI2r/hHK2dt0d6/jij20vW/EfkK5FZIlXvUsEO4qjYIzbFj2GIYhpEjawYA4BZ19uxlZeZfx7/W7NLnA3/Qqaiz8vD0kDXJmqqMbXrhEkF6Y3h/Va5TwYQW35l2bzqsse9G6tyJS7J4WGRYU38otukFixZQ7/faqFz4vTnQ0tvT9rV79M3LP+v0sXPy8LTImpR6/9qmFypeUC9/85Qq1QrNgZbenvZsOaqfP1ik8ycvZ3j8Bhbx1VNvNVdo1RI50NLb166NBzXqzZk6ezxaHh4WWZ3sY9v0oGL+6v9hB5WvXjrDei0WKSjI123tJHABAHCDjAKXYRia+OUMTfhiRpr/yN/Iw8NDVqtVvV7rqB6vdrwlurncqgzD0JwxqzX7x1VpnqjeyFaufd/6ate7Hvs3HYZhaPK3kfr1m8gsHL/J5bq91EZdXmzD/k2HYRiaO269Iseuy/Lx2+bZWmr9dE32bwYMw9DMUcs1Y+TyTO9j2zHccUAjdejXKN197O7ARZdCk/366696//33FRYWpunTpzstExqa/GvRG2+8oaefftph3syZMzVkyBD99ttvCggI0MMPP5yp9S5dulTHjh1Tz5499d1336lFixYO848ePaqff/5Za9as0cmTJyVJxYsXV82aNdW5c2eVK1fOXnbYsGEaPnx4mutavXq1goODFRUVZW/f999/r+bNmzuUs9Wzdu1aBQYGplmfbZttvL29VaxYMdWtW1cDBgxQUFCQJGn9+vXq2bOnvZyHh4f8/f1VvXp1vfDCCypbtqzT+pcvX67//e9/2rlzp+Li4lSsWDE1btxYffr0UUBAgNNl1q9fr4kTJ2rr1q26ePGifH19VblyZXXo0EHNmjWTJIftd+aVV15Rnz59JElWq1URERH69ddfdfjwYSUkJKhQoUKqXLmyunbtqipVqtiXi4qK0ogRI7Rx40adOnVKBQoUUKlSpVSzZk0NGjTIXq5Hjx66cOGC5s6da5/WuHFjHTt2TN27d9dbb72VapvccXxkZOXKlerTp4+Cg4P1xx9/yMMjdZeJrLbT2THi5+en0NBQNWjQQB06dFD+/PkzbNuNx1CuXLlUpEgRhYeHa+DAgSpZsqSk1J+txWJRgQIFFBYWpueee04PPvigQ72Z/c5I0vnz5zVy5EitXr1ax48fV758+ez7esCAAcqXL58kafDgwZo1a5bT+ry9vbVjx45U2zRjxgxVrFjRoezgwYO1cOFCbd26NdV+TEvx4sW1bNky+/tOnTppx44deuedd9S1a9dU683Inj17HP62VapUyWH+5s2bNW7cOG3dulWXLl1SoUKFVK9ePfXr10/FihVzKGvb1wULFtTSpUuVJ08eh/mNGzfW/fffrx9//DFTbcssW9iSlKmT1eRyyVe/bMv1fK2TW9t0J7GFLUmZOpFKWc62XPs+9c1p3B3AFrakrBy/yeVsy3V9qa05jbsD2MKWlPXj17Zcm2dqmdO4O4QtbEmZ38e2Y9i2XMf+jc1pnBMELpNFRkaqePHi2r59uw4fPqx77027q8NPP/2kLl26pDphsAkMDNTnn3/uMO3nn3/WyZMnU500BQYG6tixY07rWb58uV566SV5enqqTZs2KleunDw8PHTgwAEtWrRIkydP1tKlS1W8eHGH5d59913lzZs3VX0FChRINW3EiBFq1qyZS7/QDBo0SCVKlFB8fLw2b96syZMna+XKlZo7d67DPurRo4cqVaqkxMRE7dmzR1OmTNH69es1d+5c+0mtzWeffaZx48apXLlyevbZZ+Xv76+///5b//vf/zRv3jyNHz9eZcqUcVjm+++/14gRI1SqVCl17txZxYoVU3R0tFauXKmBAwfqyy+/VJs2bezlW7durYceeijV9lSocL0b0Ycffqhff/1VDz/8sNq0aSNPT08dPHhQq1atUsmSJe2B6/Dhw+rUqZN8fHzUsWNHlShRQqdPn9auXbs0ZswYh8CVnmnTpqlPnz4qXLhwhmWze3ykJSIiQsWLF9exY8e0bt061alTxy3tlK4fI4mJiTp79qw2bNigjz/+WOPHj9fIkSMzHQxTHkO7du3S1KlTtXLlSkVERDi0xfbZWq1WHTp0SJMmTVLPnj3122+/2X84SSmj70x0dLQ6duyomJgYdezYUWXKlFF0dLT27NmjyZMnq0uXLvbAJSUHqw8//DBVfZ6ezvukDx8+XKNGjUpzu6tXr57qb8qbb76psLAwPf744/ZpKdtw6NAh7dixQ8WLF1dkZKQ9cJUtWzZVXV9//bXy5s2rfv36pdmGG02cOFEfffSRSpYsqe7duys4OFgHDhzQ9OnT9fvvv2v06NGqWrVqquXOnTunyZMnp/rRygx/rdllD03ZNeGLGapctwLdC53YvemwPTRl1+wfV6lc+D10L3Ri+9o99tCUXb9+E6lKtUPpXujEni1H7aEpuyLHrlPIgyXoXpiGXRsP2kNTds0YuVwVqpXOVPdCdyBwmejo0aPaunWrhg8frrfffluRkZF6/vnnnZYtX768du/erSlTpuipp55yWiZv3rxq166dw7Tff/9dly5dSjU9LUeOHNHLL7+sYsWKafz48SpUqJDD/FdffVWTJk1yehWiefPm6V6ZunFbFi9ebL/6kx0PPfSQ/Zfvxx57TP7+/vr555+1dOlStW7d2l6uWrVqDldoSpcurXfffVezZ89W79697dPnzp2rcePGqWXLlvryyy/tJ6mPPfaYOnTooJ49e+qFF17QrFmzlCtX8ldjwYIFGjFihJo3b66vvvpKXl5e9vqeffZZrVq1SomJiQ7trlChQrqfx9mzZzVp0iQ9/vjj+uCDDxzmGYah8+fP29+PHz9ecXFxmj17dqqAc+7cufR34H/uv/9+HTx4UGPGjNGbb76ZbllXjg9n4uLitGzZMr388suaOXOmIiMj0wxcWWmnTcpjRJL69u2rtWvXql+/fhowYIB+//135c6dO8N6Uh5DHTt2VKlSpfThhx9q9uzZ6tu3r73cjZ9teHi4evfurcmTJ+vdd99NVW9G35nffvtNx48f1+TJk1OFiJiYGIfjTUq+ApfZ73r58uW1fPly/f3333rggQeclilZsqT9Kp7Nu+++q5IlS6a5noiICBUsWFCDBw/WoEGDFBUVpRIlSigoKCjVMmPGjFFAQECm27x582Z9/PHHCg8P19ixYx1+WOnSpYu6dOmiQYMGad68efLz80u1vT/99JO6du2aqc88u5ISk/T5wB8y3Q0rLR4eHvrs+R80ccO3DESQQlKiVWPfjcx0F6G0WDwsGvNOpD6fPYCBCFJISkzSNy//7Ibj16KvX/pZY1d9xPGbQlKiVT9/sMgtx+/PHyzUR9Of4vi9QVJikka9OdPlfezhYdEPb87UN/NevCnHMJ+iiSIjI+Xn56cGDRqoefPmioxM+xelqlWrqlatWho7dqyuXr2aZjlXjR07VnFxcfrkk09SnUxLySd0PXv2VNGiRbO9jpYtW6pUqVIaMWKE3HmLYK1ayZfXo6Ki0i1XrVo1ScmBN6Xhw4fLz89PH3zwQaorAmFhYXr22We1d+9eLVy40D79u+++k7+/vz7++ONUJ7+SVL9+fTVq1ChL2xEVFSXDMJz+Sm+xWFSwYEH7+yNHjqhw4cJOryalLJee4sWLq127dpo2bZpOnTqVbll3Hx+LFy/W1atX1aJFC7Vs2VKLFi3StWvXnJbNSjvTU7t2bQ0YMEDHjh1TREREtupw9VjLrCNHjsjT09OhC6lN/vz55ePjk616Jal79+7y8/PTsGHDsl2HM3PnzlXz5s3VsGFD+fr6OnRhddXIkSNlsVj06aefprrSf8899+i1117TmTNnNHXq1FTLPvfcczp79qwmT57stvY48+eCzToVddalk1UpuXvhqaizWrtwi5tadmfY+sdenTtxyaUTKSm5i9G5E5f01x//uqlld4Z1i7fp9LFzbjh+DZ0+dk7rl2xzU8vuDNtWH9D5k5fdcvyeP3lZ21cfcFPL7hybV/yjs8ejXd7HVquhs8ejtWXlHje1LH0ELhNFRkaqadOm8vb2VuvWrXXo0CFt3749zfIDBw40/YRh+fLluvfee1W5cuUsL3vx4kWdP3/e4b9Lly6lKufp6an+/fvrn3/+0eLFi93RbEnJJ6eS5O/vn245W1fKlF0dDx06pIMHD+rhhx9O896e9u3bS0reR7ZlDhw4kO4yzly5ciXVfjp//rz9SpjtHpQFCxboypUr6dZVvHhxnTx5UmvXrs30+p3p37+/kpKSNGbMmHTLuXJ8OBMZGamaNWsqODhYrVq1UmxsrMO9QNltZ0ZsV1RWr16dreVdOdZSyug7U7x4cSUlJWnOnDmZbpuzYysmJiZVufz586tXr172q1zusG3bNh0+fFitWrWSt7e3mjZtmu4PSVlx5coVrVu3TuHh4amuutm0bNlS3t7e9u9oSuHh4TflR6v5k1fIw9M9/3R6eHpo/uQVbqnrTrFqznZZPNwzWICHh0WrIggEKS2euloenm7av54eWjx1jVvqulOsifzbrcfvmrnu+dt9J1k5a4s83LiPV8y6OT96EbhMsnPnTh04cECtWrWSlHwyUKRIkXRPTqpVq6aaNWvqp59+MuWEISYmRqdPn9b999+fat6lS5ccTuCcrb9FixaqXbu2w38p7/NIqU2bNi5f5YqJidH58+d18uRJ/f777xoxYoRy586d6opSbGyszp8/r9OnT2vVqlX6+OOPZbFYHLoz7tu3T5Kc3mdjU6JECeXPn18HDiT/orR//35JUkhISJbaPWzYsFT7qXbt2tq5c6ckqVChQmrfvr1WrFihBg0a6Pnnn9e4cePs60upR48e8vLy0pNPPqn27dvro48+0pIlSzIMajcqWbKk2rZtq2nTpun06dNOy7h6fNzo3LlzWrt2rf07UKxYMVWpUiXd70Bm2pkZRYoUka+vb6avPKU8hlauXKmPPvoo1TEkXQ/TZ8+e1aZNm+z3Tt44QIxNRt+Zjh07KjAwUIMHD9Yjjzyid955R3PnztXly5ed1hcXF+f02HrhhReclu/Zs6f8/PzSHcAjKyIiIlS0aFGFh4dLklq1aqV9+/Zp9+7dLtd9+PBhJSYmpnvfnbe3t0qXLm3/jt7o+eef19mzZzVlyhSX2+OMYRj6e8Nep0O/Z4c1yaq/N+xxa0+A25lhGPp3W5TLv1zbWK3/1cf+lZS8f3dv3u906PfssCZZtWvzfvbvfwzD0P4dx916/O7fcZz9m4JhGNr71xGXr9DaWK3J9d2Mfcw9XCaJjIxUUFCQatasKSm5q1jLli0VERGhwYMHp3mT+8CBA9W9e3dNmTJFTz75pFvbZPsV3NlN/D169NA///xjf//666/rmWeecSgzbNiwVFd60hrgw3aV64033tCSJUvUtGnTLLf3xu0vXry4vvzyy1QDKgwdOtThvW1wkbCwMPu02NjkB5GmvPnfmXz58tn3k+3/GS1zo86dO6ca9U+S7rvvPvvrTz75RGFhYZoxY4YWL16sxYsX67PPPlOtWrX0+eef27fx/vvv1+zZszVy5EitWLFCu3fv1i+//KK8efNqyJAhaQZeZwYMGKCIiAiNHj3a6T1Srh4fN5o3b16q0NK6dWt9+umnunjxYqp7cDLbzszKmzev/XPPiLNj6NNPP001et6wYcMcuujlzZtXgwcPdvp528qn950JCgrSnDlzNGLECC1ZskRTpkzRlClT5OXlpf79+2vAgAEOA8/4+Pg4HQQjrdE1fX191bNnTw0bNky7du1yGLglqxITE/X777+rffv29jbVqlVLBQsWVEREhMqXL5/tuqXMf99SfkdvVL16ddWsWVNjx47VE088kal7ueLj4xUfH+8wLX/+/HI23k/02UsODzV2h8vRsbp0/rL8g5xfJb2bXL4Q5/BQY3eIvXRVMRfjVCAga3/H70QXz112eKixO8REx+ryhRj5FXTf8Nm3q8sXrjg81NgdYi9dU+zFK/INSP3v8t3o0vlYh4cau0PsxSuKiY5TgUDHvxHuHpWfwGWCpKQkzZs3TzVr1nS4ByQsLEzjxo3T2rVrVa9ePafL3njC4E62E5m4uNR/cN9//33Fxsbq7Nmzeu2115wuX61atUwNmmHTpk0bjRw5UiNGjFCTJk2y3N63335bpUuXlqenp4KCglS6dGmngzU899xzqlatmuLi4rR48WLNmzcvVTnbtmd0Ah4bG2u/N8p2opzZk3abe++9N92R+KTkG+a7deumbt266cKFC9qyZYumTJmiP/74Qy+99JImTZpkL1u6dGl98cUXSkpK0r59+7RixQqNHTtWb731lkqUKJHhumxSXj2yDU+fkqvHx40iIiIUFham6OhoRUdHS0oe2CAhIUELFixQ586ds9XOzIqLi8v0fW62Y8jDw0MBAQEqW7asfeCUlGxh+tq1a1q3bp0mTpyopKSkNOvNzHemUKFCeu+99/Tuu+/q0KFDWr16tcaMGaPvv/9ehQoV0mOPPWYv6+npmenP26ZXr16aMGGChg0bph9++CFLy6a0Zs0anT9/XmFhYTp8+LB9es2aNTVv3jy99tprmR5MxZnMft9iY2PTDWVZ/dHqxx9/dLgCmC9fPm3ZskUFnZxAXotx78mqTZ7cudz6vJfbVdKVxIwLZUO+3D7sX0nxsVnrGZFZHL/JrFfT/rfAFfnycPzaJMbFZ1woG/Lm9jJ9HxO4TLBu3TqdOXNG8+bN07x581LNj4yMTDNwScndYnr06KEpU6akeW9Idvj6+io4OFj//pv6JmLbPTsZDRKQFbarXIMHD9bSpUuzvHxYWFiqKwzOhISE2E9CmzRpoitXruitt95SeHi4fXAH2zO59uxJ++bIY8eOKSYmxl7WNjz83r17s9z2rLA9X+3hhx9Wjx49tGHDBh07dizVQBmenp4KDQ1VaGioqlSpop49e6Y76p8z/fv3V0REhMaMGZMqBLvz+LANHS7J6UiVkZGRaQaujNqZGSdPntTly5d1zz33ZKp8ymMoPSnDdKNGjeTh4aGvvvpKNWvWzNSxmh6LxaLSpUurdOnSatiwoZo1a6aIiAiHwJUdvr6+6tWrl/0qV3bZBiB58cUXnc7fsGGDfbCR7LjnnnuUK1eudL+j8fHxOnjwYKpni6VUvXp11ahRI9M/WvXt29fpyLDnzqV+8HHclYQM68uOK1cTdfas826kd5PYK+69OmCv9+o19q84fs0WE2fS8XuF49cm9oo5gSvuakKqfWyxyOkPb9lF4DJBZGSkChYsqLfffjvVPFv3sffeey/N7i41atSwnzAMGDDArW1r2LChpk+fru3btzt0uTNL27Zt9cMPP2j48OFq3PjmPGDu1Vdf1ZIlS/TDDz/o/fffl5R8lahUqVJaunSpYmJinA6CMXv2bEmy3yNmO/ldunRphr+qu0vFihW1YcMGnTlzJt3nXNlOOLN6n9M999yjtm3baurUqU4HxnDX8REZGSkvLy99/vnnqa56bN68WRMnTtTx48dTPcQ2s+3MiG0QivR+2HCH/v37a/r06fr222/1008/ua3ekiVLqkCBAjpz5oxb6rNd5Ro+fHi2fsSxDe/fsmVLp/erffjhh4qMjHQpcOXNm1c1a9bUunXrnP7gICU/BiM+Pj7DkUEHDhxo/9EqI97e3vL29k413TCUKnD5FSyg/H753Nqt0Nc/nwoE+qZa190ov39e5fXN7dZuhfkK5FZ+v7zsX0kFAn2V3y+vW7sV5vfPJ9+A/OxfSfn98yivr49buxXmK+CjfH552L//8Q3Ip3wFcru1W2E+vzzK72/+3wgGzXCzq1evatGiRWrYsKFatGiR6r9u3bplOFKblHzCcObMGU2bNs2t7Xv22WeVJ08eDR06VGfPnk013903Dtqucu3evTvDbXaXe+65R82aNdOsWbMcTlife+45Xbx4Ue+8806qbmA7d+7U2LFjFRIS4nBFZtCgQYqOjtabb76Z6nlbUvIoeM5GTEvPmTNn7IN4pBQfH6+1a9fKw8PDfmVm06ZNSkhI/avkypUrJSWHwqzq37+/EhMTNXbs2FTz3HV8REZGKjw8XC1btkz1HXj22WclKcPhxNNrZ3rWrl2rkSNHqkSJEmrbtm2Wls2qAgUKqHPnzlq9enW2Bo7Ytm2b0y6c27dvV3R0dLY+X2dsV7mWLl2arXYuXrxYcXFx6tatm9O/a40aNdKiRYtS3QuVVf3795dhGBo8eHCqgVmOHj2qL7/8UsHBweleHZUcf7RK6zEE2WGxWPRAjRC3jlL4QI1Qlx4QfyexWCy6v3IJt47ydn/lEuzf/1gsFpUPL+vWUQorhJdl//7HYrGobKVibj1+y1Yqxv5NwWKxKKTKPW4dpTCkyj03ZR9zhcvNli1bptjY2DSv5lSpUkWBgYGKiIhQy5Yt06zHdsKwYcMGt7avVKlS+vLLL/XKK6+oRYsWatOmjcqVKyfDMBQVFaW5c+fKw8NDRYoUSbXswoULnQ6oULduXQUFBaW5Ttu9XO4YySyznnnmGc2fP18TJkzQq6++Kin5atuOHTv0yy+/aP/+/WrTpo0KFCigXbt2acaMGfL399d3333n8Lytli1bas+ePRo1apR27dql1q1bq1ixYoqOjtaqVau0du1affXVVw7r3rVrl9Nhvu+55x49+OCDOnnypB577DHVqlVLtWvXVlBQkM6dO6d58+bpn3/+Ua9evez3/YwZM0Z///23mjZtah9hcdeuXZo9e7b8/f3Vq1evLO8b29WjWbNmpZrnyvFhYxs6vFu3bk7nFy5cWBUqVFBkZGS692il106bP/74QwcOHFBSUpLOnj2r9evXa82aNSpWrJh++OEHl55jlVk9e/bUhAkTNHr0aH3zzTcO8zL6zsyZM0eRkZFq0qSJKlasKC8vL+3fv18zZsyQj4+P+vXr57BcYmJimkPIN23a1Om6UrZz/Pjx+ueff9It50xkZKT8/f314IMPOp3fuHFjTZs2TStWrHDpYefVq1fXG2+8oU8++URt27bVo48+quDgYB04cEDTp0+X1WrV6NGj0xxwJaXnn39ePXv2zHZb0vJIl4Zav3irW+qyJln1SJeGbqnrTlG/XZi2rU79g1R2WK2G6rd1zyMu7hRNO9fTxmU73FKXNcmqpp3ruqWuO0XdNg9ox58H3VKX1WqobmvnD62/mzV4tKq2/uGeWz2sVkMNH039TFQzELjcLCIiQj4+Pqpb1/kfIQ8PDzVs2FCRkZG6cOFCmqOLSeadMDRp0kSRkZEaN26c1qxZoxkzZshisahYsWJq0KCBunTp4nRo5nfffddpfb/88ku6gStXrlzq37+/fQjtm6FSpUqqUaOGJk+erL59+8rXN7kf7v/93/+pZs2amjRpkn788UdduXJFRYsWVbdu3dS7d2+nAxy89NJLqlWrliZOnKjJkyfr4sWLKlCggCpXrqyRI0fq4Ycfdig/d+5cp1dvHn30UT344IMqXbq0hg4dqpUrV2rSpEk6d+6cvL29FRISog8//FCdOnWyL9O3b1/NnTtXGzduVGRkpK5evWp/ptWAAQPSfF5RRmz3SDkb8CG7x4eNbdj39LqQNm7cWMOGDdM///yTbl3ptVOSvv/+e0mSl5eX/P39FRISoqFDh6pDhw5ZenaaKwoXLqw2bdpozpw5OnLkiMN9Yxl9Zzp37qzcuXNr3bp1WrZsmWJiYhQQEKC6deuqb9++qUYVjI+P1+uvv+60zqVLl6YbpAoUKKBevXpleYj4lMP7pzW6au3atZUnTx5FRES4FLik5NFJK1asqHHjxmnChAmKiYlRcHCwWrRooX79+qXb1TalmjVrmvKjVZ0W4SpcIkhnjrv28FgPDw8FFwtU7eY35x/728WDD4WoYNECOn/KtYfHWjwsCizsqyoPpX7Mxd2sVtPKKlS8oM6eOO/i8WtRUNFA1WxCoE2pcr0yCiziqwunY1w+fgMK5VdYvTJubN2dIbxhOQUV89e5kxdd2sceHhYFFvFT1QZpPy7InSwGA/wDAODg7NnUg2bY/LVml1559AOX1/H17LdUuU72h+q/U+3edFif9f3V5XoGj+6mcuH3uqFFd5bta/doSOcvXa7n02mvqlKtm3OyejvZs+Wovnpuhsv1vDKik0KrlnBDi+48uzYe1IdPj3O5nrfGPa3y1Z133bdY5NaRC7mHCwCALKhSt4J6vdbRpTp6vdaRsJWG8tXuVfu+9V2qo33f+oStNITVDlW3l9q4VEe3l9oQttIQWrWk2jyb/QGEJKnNs7UIW+moUL20Og5If/CkjHQc0CjNsGUGAhcAAFnU49WO9tCV2eeP2cr1eq2jerzqWmC707XrXc8eujI7CIGtXPu+9dWut7kjlN7uurzYxh66MjsAga1ct5faqMuLrgW2O13rp2vaQ1dWj982z9ZS66drmta2O0WHfo3soSurx3DHAY3UoZ9rgS2r6FIIAMAN0utSmNK2P3fps+d/0Kmos/Lw9JA1yZqqjG164RJBemN4f65sZcE/mw9rzDuROnfikjw8LE7vO7JNL1i0gHq/14YrW1mwY90eff3Szzp97FyGx2+h4gX18jdPcWUrC/ZsidLPHyzU+ZOXMzx+A4v46qm3mnNlK4t2bzyoH96cqbPHozPcx0HF/NX/ww6ZurLl7i6FBC4AAG6Q2cAlSUmJSVq7cIt+n7Rcuzbu1eXo68/p8vXPpwdqhOqRrg1Vu1lVeeZyPvAI0paUaNVff/yrP+Zs077tUQ7P4MlXILfur1xC9dtVVpX698szFx13siopMUnrl2zToimrtXvLAcWkOH7z++dThfCyavZEPdV4OIzjNxuSEq3avvqAVkf+rQM7jyv20vVHVeQr4KOylYqpbpsHFFa3DMdvNiUlJmnLyj1aMXOz9m47qtiLV+zz8vnlUUiVe9Tw0aqq2iA008cwgQsAAJNlJXClZBiGLp2/rDy5c+nK1UQVCPTlOTpuZBiGYi7GKV9uH8Vevab8fnnZv25kGIYuX4ixH7++AfnZv25kGIZiL15Rvjw+ir1yTfn88rB/3cwwDMVExylvbi/FXU1Qfv/s/Y0gcAEAYLLsBi7p+j/UrtSBtLF/zcX+NRf713zu2MeMUggAAAAAtwkCFwAAAACYhMAFAAAAACYhcAEAAACASQhcAAAAAGASAhcAAAAAmITABQAAAAAmIXABAAAAgEkIXAAAAABgEgIXAAAAAJiEwAUAAAAAJiFwAQAAAIBJCFwAAAAAYBICFwAAAACYhMAFAAAAACYhcAEAAACASQhcAAAAAGASAhcAAAAAmITABQAAAAAmIXABAAAAgEkIXAAAAABgEgIXAAAAAJiEwAUAAAAAJiFwAQAAAIBJCFwAAAAAYBICFwAAAACYhMAFAAAAACYhcAEAAACASQhcAAAAAGASAhcAAAAAmITABQAAAAAmIXABAAAAgEkIXAAAAABgEgIXAAAAAJiEwAUAAAAAJiFwAQAAAIBJCFwAAAAAYBICFwAAAACYhMAFAAAAACYhcAEAAACASQhcAAAAAGASAhcAAAAAmITABQAAAAAmIXABAAAAgEkIXAAAAABgEgIXAAAAAJiEwAUAAAAAJiFwAQAAAIBJCFwAAAAAYBICFwAAAACYhMAFAAAAACYhcAEAAACASQhcAAAAAGASAhcAAAAAmITABQAAAAAmIXABAAAAgEkIXAAAAABgEgIXAAAAAJiEwAUAAAAAJiFwAQAAAIBJCFwAAAAAYBICFwAAAACYhMAFAAAAACYhcAEAAACASQhcAAAAAGASAhcAAAAAmITABQAAAAAmIXABAAAAgEkIXAAAAABgEgIXAAAAAJiEwAUAAAAAJiFwAQAAAIBJCFwAAAAAYBICFwAAAACYhMAFAAAAACYhcAEAAACASQhcAAAAAGASAhcAAAAAmITABQAAAAAmIXABAAAAgEkIXAAAAABgEgIXAAAAAJiEwAUAAAAAJiFwAQAAAIBJCFwAAAAAYBICFwAAAACYhMAFAAAAACYhcAEAAACASQhcAAAAAGASAhcAAAAAmITABQAAAAAmIXABAAAAgEkIXAAAAABgEgIXAAAAAJiEwAUAAAAAJiFwAQAAAIBJCFwAAAAAYBICFwAAAACYhMAFAAAAACYhcAEAAACASQhcAAAAAGASAhcAAAAAmITABQAAAAAmIXABAAAAgEkIXAAAAABgEgIXAAAAAJiEwAUAAAAAJiFwAQAAAIBJCFwAAAAAYBICFwAAAACYhMAFAAAAACYhcAEAAACASQhcAAAAAGASAhcAAAAAmITABQAAAAAmIXABAAAAgEkIXAAAAABgEgIXAAAAAJiEwAUAAAAAJiFwAQAAAIBJCFwAAAAAYBICFwAAAACYhMAFAAAAACYhcAEAAACASQhcAAAAAGASAhcAAAAAmITABQAAAAAmIXABAAAAgEkIXAAAAABgEgIXAAAAAJiEwAUAAAAAJiFwAQAAAIBJCFwAAAAAYBICFwAAAACYhMAFAAAAACYhcAEAAACASQhcAAAAAGASAhcAAAAAmITABQAAAAAmIXABAAAAgEkIXAAAAABgEgIXAAAAAJiEwAUAAAAAJiFwAQAAAIBJCFwAAAAAYBICFwAAAACYhMAFAAAAACYhcAEAAACASQhcAAAAAGASAhcAAAAAmITABQAAAAAmIXABAAAAgEkIXAAAAABgEgIXAAAAAJiEwAUAAAAAJiFwAQAAAIBJCFwAAAAAYBICFwAAAACYhMAFAAAAACYhcAEAAACASQhcAAAAAGASAhcAAAAAmITABQAAAAAmIXABAAAAgEkIXAAAAABgEgIXAAAAAJiEwAUAAAAAJiFwAQAAAIBJCFwAAAAAYBICFwAAAACYhMAFAAAAACYhcAEAAACASQhcAAAAAGASAhcAAAAAmITABQAAAAAmIXABAAAAgEkIXAAAAABgEgIXAAAAAJiEwAUAAAAAJiFwAQAAAIBJCFwAAAAAYBICFwAAAACYhMAFAAAAACYhcAEAAACASQhcAAAAAGASAhcAAAAAmITABQAAAAAmIXABAAAAgElymb2CZcuWaf78+bpw4YJKlCihxx57TA888IDZqwUAAACAHOdS4Fq3bp1eeukl+fj4KCIiQgUKFHCY/+233+rHH390mPbbb7/pww8/VPv27V1ZNQAAAADc8lzqUrhy5UpduHBBlSpVShW2/vnnH/34448yDEOGYahAgQIyDEOJiYl65513FBUV5VLDAQAAAOBW51Lg2rJliywWi2rXrp1q3uTJk2UYhvz8/DRz5kytX79e06dPl5+fn+Lj4zVlyhRXVg0AAAAAtzyXAteZM2ckSffff3+qeStWrJDFYlG3bt1UoUIFSVKlSpXUvXt3GYahP//805VVAwAAAMAtz6V7uM6fPy9J8vX1dZh+5MgRnTp1ShaLRU2bNnWYV61aNUnS0aNHXVk1AAAA3MgwDCXGnFWccU4JMVZ55guSxWLJ6WYBtz2XApdhGJKkmJgYh+mbNm2SlBzEypcv7zDP399fknTlyhVXVg0AAAAXGUmJiv57gc6un6TYQxuVdCXaPs8zj7/ylaquoJrd5P9Ac1k8TR/cGrgjufTNCQoK0vHjx7V//377lStJWr16tSSpatWqqZaJi4uTJPn5+bmyagAAALjg8r41OjR5oOIvHJUsnpKR5DA/6Uq0Lv2zTJd2L5Z3QEmV6jJMvvfVzaHWArcvlwJXlSpVdOzYMU2ePFlt27ZVnjx5dPToUS1dulQWi0V16tRJtcyhQ4ckJYc1AID066+/6v3331dYWJimT5+ean5oaKj9tcViUVBQkEJCQtS3b1/VrFnToWxCQoImT56sOXPm6MCBA5KkMmXKqF27durSpYu8vLwcyjdu3FjHjh2zvw8MDFTp0qX11FNPqWnTppo5c6aGDBmS4TYUL15cy5Yts7/fvHmzxo0bp61bt+rSpUsqVKiQ6tWrp379+qlYsWKplt+0aZNGjRqlPXv2KDo6WgULFlS5cuXUqlUrtWnTxmFfdOvWTW+//bZ92vnz5zVy5EitXr1ax48fV758+VS8eHHVrFlTAwYMUL58+TJsP3A3MQxDJxZ9qRMLP5cs/93Of0PYul44eXp89DHtHdleRZu/rqLNXqWrIZAFLgWuTp06ad68edqzZ49at26tBx54QBs3btS1a9eUJ08etW7dOtUytu6GpUqVcmXVAHDHiIyMVPHixbV9+3YdPnxY9957b6oydevWVbt27WQYhqKiojR58mT16tVLP/74oxo0aCApuQdB3759tWHDBjVq1EgdOnSQxWLRqlWr9NFHH2nx4sX68ccflTdvXoe6y5cvr6eeekqSdPr0aU2dOlXPP/+83n33XdWrV0+ff/65Q/k333xTYWFhevzxx+3TUoaaiRMn6qOPPlLJkiXVvXt3BQcH68CBA5o+fbp+//13jR492qEHxPz58/XSSy+pfPny6tmzp/z8/BQVFaWNGzdq2rRpDoHrRtHR0erYsaNiYmLUsWNHlSlTRtHR0dqzZ48mT56sLl26ELiAG9jDliQZ1swt9F8523LFmr9mRtOAO5JLgat27drq2bOnfvnlFx07dkzHjx+339f12muvKTAw0KH8tWvX7Fe/qlev7sqqAeCOcPToUW3dulXDhw/X22+/rcjISD3//POpypUqVUrt2rWzv2/atKnatm2rX375xR64Pv30U23YsEFvvfWWunfvbi/btWtX+1W0zz77TO+9955D3YULF3aou3379mrWrJnGjx+vLl26qGTJkg7l3333XZUsWdJhGZvNmzfr448/Vnh4uMaOHas8efLY53Xp0kVdunTRoEGDNG/ePHvX8uHDh+u+++7T1KlT5e3t7VDfuXPn0t1/v/32m44fP67Jkyen6sYeExOT6ooecLe7vG/N9bCVTScWfi7fsnXoXpgN164m6Pzpy7oWnyBvH/4+3S1cvvtx6NChql27thYsWKCzZ88qODhY7dq1c/psrqVLlyp//vzKnz+/GjVq5OqqAeC2FxkZKT8/PzVo0EDNmzdPM3DdKDQ0VAEBAfaHyJ88eVK//fabatWq5RC2bLp166aFCxfqt99+U//+/VWkSJE06w4ODlaZMmW0Z8+eLG/PyJEjZbFY9OmnnzqELUm655579Nprr+mNN97Q1KlT1adPH0nJI9u2bt06VdiSpIIFC6a7viNHjsjT01NVqlRJNS9//vxZbr+rLpyJ0cbFu5VwJVFeeXKpRrPy8g+6+e0AnDGSEnVo8sDkboSZvbLljMVDhyYPVMWhGxhII5P2bo3S4slbtH31ARmGZLFIYfXKqGmXqgp5sERONw8mc8u3pFGjRpkKUC1btlTLli3dsUoAuCNERkaqadOm8vb2VuvWrTV58mRt375dYWFh6S538eJFXbp0yd798I8//lBSUpLat2+f5jLt27fX+vXrtWrVKj322GNplktISNDJkyfto8pm1pUrV7Ru3TqFh4enuipm07JlS7311ltavny5PXAVK1ZMa9eu1cmTJ9MNgs4UL15cSUlJmjNnjh599NEsLetO1iSrpg1boaVTtsiQ5OFhkdVqaPr3K/XwE1X1+MCG8vB06dGXgMui/16QPECGqwyr4i8cVfTfCxUQ1sr1+u5wK2Zs06Qvl8vD06L/OoLJMKQdfx7UtlUH1PXVRmrYsXLONhKm4q8/AOSQnTt36sCBA2rVKvmEJTw8XEWKFFFkZGSqsteuXdP58+d1/vx5bd++XS+88IKSkpLUokULSdK+ffskSeXKlUtzfbZ5+/fvd5iemJhor/uff/7RG2+8obNnz9rrzqzDhw8rMTEx3TZ4e3urdOnS9gE9JKl37946ceKEmjRpop49e+q7777Tpk2bZLVm/At8x44dFRgYqMGDB+uRRx7RO++8o7lz5+ry5ctZarurpg1bocWTNstqNWRYDSUlWmVYDVmthhZP2qzpw1be1PYAzpxdPyl5NEJ3sHjq3IZJ7qnrDrZ3a5QmfblckmRNMhzm2d5P+nK5/v3rWKplcefgOjAA5JDIyEgFBQXZRxq0WCxq2bKlIiIiNHjwYHl6Xj8x+u233/Tbb7/Z3/v4+Oipp55Sr169JEmxsbGSlO4AEbZ5Nz47cfXq1Q7dwD09PdWuXTu9+uqrWdoeW70ZDVKRL18+hzZ06tRJhQsX1vjx47V+/XqtX79eI0eOVMmSJfX55587fcSITVBQkObMmaMRI0ZoyZIlmjJliqZMmSIvLy/1799fAwYMMH00tQtnYrR0ypZ0yyyZslnNu1eneyFyjGEYij20Me3RCLNcYZJiDm2QYRiMWJiOxZO3yMPTkipspeThadGiyZt1f5XiN7FluJncGrhiYmIUFRWl2NjYTP0yycAZAO5WSUlJmjdvnmrWrGm/D0uSwsLCNG7cOK1du1b16tWzT3/44YfVvXt3WSwW5cuXT/fdd5/DaIO2kGMLXs6kFcoqV66sF198URaLRblz51bZsmVVoECBLG+T7Z6p9Npgm39jG+rXr6/69evrypUr+vvvv/X7779rypQp6tevn+bPn5/uvVyFChXSe++9p3fffVeHDh3S6tWrNWbMGH3//fcqVKhQut0n4+PjFR8fn2o7snL+uHHxbqV9KpXMkLRx8T9q1rVaBiWREYtFmjtjh2o3LJXTTbmtJMacdXiosTskxUUrKe6cvPLzqB9nrl1NsN+zlR5rkqHtqw4o/lqCfHIzkIarLBbpo8Gz1e+1h12qw53cErimTZumSZMmae/evfZRCjNisVi0a9cud6weAG4769at05kzZzRv3jzNmzcv1fzIyEiHwFWkSBGnzza0KVu2rCRpz549Kl++vNMytkEw7rvvPofpAQEB6dadWffcc49y5cqV7mAb8fHxOnjwoCpWrOh0fp48eVStWjVVq1ZNAQEBGj58uP74449M3Z9lsVhUunRplS5dWg0bNlSzZs0UERGRbuD68ccfNXz4cPv7fPnyacuWLSpY0DfD9dkkXEmUh4dFSda0//3z9LQoPi5BQUGZrxfOzZ2xQ5K0dsUhte5YKYdbc/uIM9If8TO7/PJZlJfj2qnzpy9nGLZsDEPK4+2lQPalyz4aPFuSNOqLpfq/T9vnaFtsXApcSUlJGjhwoJYvT+6bmtmwBQB3u8jISBUsWNDhAb42ixcv1uLFi/Xee+8pd+7cmarvoYcekqenp+bMmZPmwBmzZ89Wrly5VL9+fVeanqa8efOqZs2aWrdunY4dO6bixVN3j/n9998VHx+fqYGWbKHszJkzWW5LyZIlVaBAgQyX7du3r/0ZZCmdO5f5EyWvPLlkTSdsSVJSkiHvvF46e/bm3lt2J6rdsJTWrjik2g1LsT+zICHGhVEJ03Ex1lCchc/BmWvxCbJYlKm/JRaLdCU+gWPaDfq99rBGfbFU/V57ONv702JRln54y4hLgWvKlClatmyZpOR+9B06dNADDzwgPz8/eXgwHgcAOHP16lUtWrRILVq0cDowRaFChTR37lwtW7Ys0yO7Fi1aVB06dND06dM1adIkde3a1WH+5MmTtW7dOnXu3DnLIwFmRf/+/bV27VoNHjxYY8aMcQiMR48e1Zdffqng4GB17tzZPn3t2rVOHyWycmXyQBOlS5dOc33btm3T/fffn+phztu3b1d0dHS6939JyYN4OBuO3jAyd5IkSdWbltf071em263QIql603KZrhPpa92xks6ezXwohuSZL0ieefzd2q3QM6+/PPMW5HNIg7ePl8LqldGOPw9meA9Xpbql5e3jxb50k//7tP0t9TfCpcA1e/ZsScndU3799Vf7QywBAGlbtmyZYmNj1bhxY6fzq1SposDAQEVERGTpURpDhgzRgQMH9N5772nVqlX2K1mrV6/W0qVLVaNGDQ0ePNgt25CW6tWr64033tAnn3yitm3b6tFHH1VwcLAOHDig6dOny2q1avTo0Q7/XgwYMEAlSpRQo0aNVLJkSV25ckV//vmnli9frkqVKqV7NWzOnDmKjIxUkyZNVLFiRXl5eWn//v2aMWOGfHx81K9fP1O3V5ICgvPr4SeqavGkzWmWadIlnAEzkKMsFovylaquS/8sc8/AGRZP5S9VgwEzMtC0S1VtW3Ug3TLWJEPNuoTfpBYhJ7gUuPbv3y+LxaIBAwYQtgAgkyIiIuTj46O6des6ne/h4aGGDRsqMjJSFy5cyHS9+fLl0/jx4zVp0iRFREToiy++kGEYKlOmjIYOHaquXbvKy8v8G7KffPJJVaxYUePGjdOECRMUExOj4OBgtWjRQv369UvV1fDDDz/U0qVLNX/+fJ0+fVqGYahkyZLq16+fevfurVy50v6nqnPnzsqdO7fWrVunZcuWKSYmRgEBAapbt6769u2rChUqmL25kqTHBzaURRYtmbJZhpLv2UpKMmRRcth67PkGN6UdQHqCanbVpd2L3VOZkaSCNbpmXO4uF/JgCXV9tZH9OVwpr3TZ3nd9tREjFN7hLIYLN15VrVpVV65c0cyZM9O8SRsAgNtNdruiRJ+N0cbF/yg+LkHeeb1UvWk5rmy5mcUiBQX53lLdhW4XRlKidn5cQ/HRxyTDhXu6LB7y9i+uikM3yOLJE4Yy49+/jmnR5M3avip51EKLRQqrX0bNuoQTttzMHX8jbHW4i0vfklKlSmn37t2Kjo52U3MAALh9+QflV7Ou1QgEuCVZPHOpVJdh2juyvWsVGVaV6jKMsJUF91cprvurFFf8tQTl8fbSlfgEefswBPzdwqWRLVq2bCnDMLRixQo3NQcAAABm8b2vroo2f92lOoo2f12+9znvEo30+eT2UmAhX563dZdxKXD17NlT5cqV0+TJk7Vp0yZ3tQkAAAAmKdrs1euhy5LJU8H/yhVt/rqKNnvVpJYBdyaXApe3t7d++uknPfDAA3rqqaf0+eefa/fu3bp27Zq72gcAAAA3slgsKtb8NYUMmC1v///uH7J4plE4ebq3f3GFDJitYs1fY2RCIItcGjQj5UAZhmFk6QtosVi0a9eu7K4aAADTuONma+7hMgf7172MpERF/71QZ9f/qtjDG5UUF22f55nXX/lL1VBQzW7yq9CMe7bcgOPXfHfcoBk3ZjUXshsAAABuMotnLgWEtVJAWCsZhqGkuHPyy2fRxVhDnnkLcjULcAOXAtfzzz/vrnYAAAAgB1ksFnnlD1LeIF/FWbgCA7gLgQsAAAAATOLSoBkAAAAAgLQRuAAAAADAJG4dbiYhIUG7du3S3r17dfHiRUmSn5+fQkJCVKFCBXl58ZA3AAAAAHcPtwSuK1euaOTIkZo2bZouXbrktEyBAgXUuXNn9e/fX3ny5HHHagEAAADgluZyl8Ljx4+rffv2Gjt2rC5evCjDMJz+d/HiRY0ZM0bt27fXyZMn3dF2AAAAALiluXSFKyEhQb1799bhw4clSWXKlFGHDh1UuXJlBQUFSZLOnj2r7du3a9asWdq3b58OHz6s3r17a9asWcqViwfoAQAAALhzuZR4pk+frv3798tisahv374aOHCgPD09HcqULl1a1atX11NPPaVhw4bphx9+0L59+zR9+nR16dLFpcYDAAAAwK3MpS6FCxYskMViUZMmTfTiiy+mClsOK/Lw0AsvvKCmTZvKMAwtWLDAlVUDAAAAwC3PpcD177//SpI6duyY6WU6deokSdq7d68rqwYAAACAW55Lgevy5cuSpEKFCmV6meDgYElSTEyMK6sGAAAAgFueS4HL399fkhQVFZXpZWxl/fz8XFk1AAAAANzyXApcFSpUkGEY+vXXXzO9zKRJk2SxWFShQgVXVg0AAAAAtzyXAlerVq0kSRs2bNCQIUMUFxeXZtkrV67o//7v/7Ru3TpJUuvWrV1ZNQAAAADc8lwaFr5t27aaMmWKtm7dqtmzZ+uPP/7QI488orCwMBUsWFAWi8X+HK758+fr/PnzkqSqVauqbdu2btkAAAAAALhVWQzDMFyp4OLFi+rbt6/++uuv5AotFqflbKupUqWKfvzxR+7hAgDcss6evazs/utosUhBQb4u1YG0sX/Nxf41F/vXfO7Yx7Y63MWlK1xS8uAXkyZN0uTJkzVp0iTt37/fabmyZcuqa9eu6tKlizw8XOrJCAAAAAC3BZcDl5T8UONu3bqpW7duOn36tP79919dvHhRUnIgu//++7M0dDwAAAAA3AncErhSKlSoEOEKAAAAAOTiKIUAAAAAgLQRuAAAAADAJJnqUjh8+HD76+eff97p9OxIWRcAAAAA3GkyNSx8uXLl7MO979692+n07EhZFwAAtwqGhb91sX/Nxf41F/vXfLf1sPBp5TIXH+MFAAAAAHesTAWuf/75J0vTAQAAAAAMmgEAAAAApiFwAQAAAIBJXApcPXv2VK9evXTs2LFML3Pq1Cn7cgAAAABwJ8v0oBnObNiwQRaLRVeuXMn0MteuXbMvBwAAAAB3MroUAgAAAIBJbnrgiouLkyT5+Pjc7FUDAAAAwE110wPXH3/8IUkqXLjwzV41AAAAANxUWbqHa8iQIU6nf/vtt/L1Tf9pzPHx8Tp69Kh27Nghi8WiGjVqZGXVAAAAAHDbyVLgmjVrVqrBLgzD0NKlSzO1vGEYkiQ/Pz/16dMnK6sGAAAAgNtOlgJXsWLFHN4fP35cFotFwcHBypUr7aosFot8fHwUHBysBx98UF26dKFLIQAAAIA7XpYC17JlyxzelytXTpI0btw43Xfffe5rFQAAAADcAVx6Dlf16tUlSXny5HFLYwAAAADgTuJS4Jo4caK72gEAAAAAdxwefAwAAAAAJiFwAQAAAIBJXOpSaBMfH6/IyEgtWbJEe/bs0YULF3T16tV0l7FYLNq1a5c7Vg8AAAAAtySXA9fBgwf13HPP6eDBg/bnbAEAAAAAXAxccXFx6t27t6KiouTh4aGHH35YgYGBmjZtmiwWi/r376+LFy9q586d2rZtmywWi6pUqaK6deu6q/0AAAAAcMtyKXBNmTJFUVFR8vT01NixY1W7dm39+++/mjZtmiRp0KBB9rK7du3S66+/rm3btqlVq1bq3r27ay0HAAAAgFucS4NmLF++XBaLRS1atFDt2rXTLVuhQgX98ssvCgwM1KeffqqdO3e6smoAAAAAuOW5FLj27dsnSWratKnT+Var1eF9YGCgnnrqKSUmJurXX391ZdUAAAAAcMtzKXBdvnxZklSsWDH7NG9vb/vruLi4VMtUrVpVkrRhwwZXVg0AAAAAtzyXAlfu3LklJQ/xbuPr62t/feLEiTSXPXv2rCurBgAAAIBbnkuBq0SJEpKk06dP26cFBgbKz89PkrRly5ZUy/z999+SJC8vL1dWDQAAAAC3PJcCV8WKFSUp1QAYtWvXlmEY+umnnxQdHW2ffvToUY0ePVoWi0Xly5d3ZdUAAAAAcMtzKXDVqVNHhmFo2bJlDtN79OghKTlgNW/eXIMGDVLv3r3Vrl07+9Wwxx9/3JVVAwAAAMAtz6XA1ahRI1WvXl358uXTkSNH7NPDw8P13HPPyTAMXbx4UYsXL9bq1avtg2h06NBBbdq0ca3lAAAAAHCLc+nBx3ny5NHEiROdzhs4cKCqVaum6dOn699//9X/t3ff4VHUaxvH701ISKMHQarUhBY6gnAoASmho9ioCgoiggUV1BcVURA9cKRYUBFQpKgQklAE6SCgVOm9g0gILb3t+wcnewjZJLvZHZKN38915dLMzs48+2QS9t75zW9SUlJUsWJF9ejRQx06dHBktwAAAADgEhwKXNlp1qxZtjdEBgAAAID8yqEhhQAAAACAzDkUuIKDg9WuXTudOXPG5udcvHhRbdu2Vbt27RzZNQAAAADkeQ4NKbx48aJMJpOSkpJsfk5ycrIuXLiQ7mbJAAAAAJAfMaQQAAAAAAxyzwPXrVu3JEleXl73etcAAAAAcE/d88AVFhYmSSpbtuy93jUAAAAA3FN2XcPVv39/q8vHjBkjb2/vLJ+bmJio8+fP6+rVqzKZTGrevLk9uwYAAAAAl2NX4Pr9999lMplkNpsty8xms/bt22fXTsuXL68hQ4bY9RwAAAAAcDV2Ba7GjRun+/6PP/6QyWRSrVq1sjzDZTKZVLBgQZUsWVL169dX586d5ePjk7OKAQAAAMBF2BW4vvvuu3TfBwYGSpImTpyoqlWrOq8qAAAAAMgHHLoPV48ePWQymVS4cGFn1QMAAAAA+YZDgWvixInOqgMAAAAA8h2HAldmkpOTdfPmTUlS4cKFVaCAIbsBAAAAgDzNaUnoxIkT+uGHH/Tbb7/pzJkzlpkMTSaTKlasqObNm+uJJ57gWi8AAAAA/xhOCVz//ve/NWvWLKWmpqabMl66PW38qVOndPr0ac2fP1+DBg3SK6+84ozdAgAAAECe5nDgev/99/XDDz9YglaVKlUUFBSkkiVLSpKuXLmiffv26fjx40pJSdFXX32l2NhYvf32247uGgAAAADyNIcC186dOzVv3jyZTCZVrVpV48aNU4MGDayuu3v3br3zzjs6evSo5s2bp5CQkEzXBQAAAID8wM2RJy9cuFCSVK5cOc2fPz/LAFW/fn3NmzdP5cuXlyQtWLDAkV0DAAAAQJ7nUODasWOHTCaTnnvuORUqVCjb9QsVKqRnn31WZrNZO3bscGTXAAAAAJDnORS4IiMjJUk1atSw+Tk1a9aUJF29etWRXQMAAABAnudQ4CpYsKAkKS4uzubnpK3r6enpyK4BAAAAIM9zKHCVK1dOkrR27Vqbn5O2btq1XAAAAACQXzkUuFq2bCmz2azvv/9eW7duzXb9bdu26fvvv5fJZFLLli0d2TUAAAAA5HkOBa4BAwbIz89PycnJevbZZzVu3DgdOHBAqamplnVSU1N14MABjRs3ToMHD1ZSUpL8/Pw0YMAAh4sHAAAAgLzMoftwFS9eXP/5z3/0/PPPKykpSfPnz9f8+fPl4eGhIkWKyGQy6fr160pKSpIkmc1meXh46NNPP1WxYsWc8gIAAAAAIK9y6AyXJLVo0UKLFi1S7dq1ZTabZTablZiYqCtXrujvv/9WYmKiZXnt2rX1448/6qGHHnJG7QAAAACQpzl0hitNjRo19NNPP+nPP//U1q1bdfToUd24cUOSVKRIEVWvXl3NmjVTUFCQM3YHAAAAAC7BKYErTVBQEKEKAAAAAP7L4SGFAAAAAADrCFwAAAAAYBCnDSlMTU3V8ePHde7cOcXExKSbGj4zPXr0cNbuAQAAACDPcThwxcfH6/PPP9eiRYt0/fp1m59nMpkIXAAAAADyNYcCV3x8vAYMGKA///xTZrPZWTUBAAAAQL7gUOCaPXu29u7dK0mqVq2a+vbtq1q1aqlIkSJyc+PyMAAAAAD/bA4FruXLl0uS6tevrzlz5sjT09MpRQEAAABAfuDQaahz587JZDJp8ODBhC0AAAAAuItDgcvDw0OSVKZMGacUAwAAAAD5iUOBq3LlypKkK1euOKUYAAAAAMhPHApcPXv2lNls1sqVK51VDwAAAADkGw4Frscee0xNmzbV0qVLFRER4ayaAAAAACBfcGiWwkuXLun//u//9Pbbb+u1117T6tWr1bVrV1WqVEne3t7ZPp9rvwAAAADkZw4FruDgYJlMJkmS2WzWqlWrtGrVKpueazKZdPDgQUd2DwAAAAB5mkOBS7odtKz9PwAAAAD80zkUuCZMmOCsOgAAAAAg33EocPXs2dNZdQAAAABAvuPQLIUAAAAAgMwRuAAAAADAIAQuAAAAADCIw7MUSlJycrLWr1+vnTt36ty5c4qJiVFKSkqWzzGZTJozZ44zdg8AAAAAeZLDgWvHjh16/fXXdenSJcuyrKaHN5lMMpvNlvt3AQAAAEB+5VDgOnHihJ599lnFx8fLbDbLw8NDDzzwgIoUKUKgAgAAAPCP51Dg+vLLLxUXFyd3d3e9+OKL6tevn3x9fZ1VGwAAAAC4NIcC17Zt22QymdS/f38NHTrUWTUBAAAAQL7g0CyF165dkyS1a9fOKcUAAAAAQH7iUOAqXry4JMnLy8spxQAAAABAfuJQ4GrYsKEk6dixY04pBgAAAADyE4cC18CBA+Xu7q65c+cqOTnZWTUBAAAAQL7gUOAKCgrSmDFjdPjwYQ0fPlxRUVHOqgsAAAAAXJ5DsxROnz5d0u3gtX79egUHB+uhhx5S5cqVbbqua/jw4Y7sHgAAAADyNJPZbDbn9MmBgYHpbnBsNpvtuuHxoUOHcrprAAAMExl5Szn919Fkkvz9Czm0DWSO/hqL/hqL/hrPGT1O24azOHSGS7odsrL6HgAAAAD+qRwKXIcPH3ZWHQAAAACQ7zg0aQYAAAAAIHO5FriuXbuWW7sGAAAAgHvingaupKQk/fLLL3r++efVsmXLe7lrAAAAALjnHJ40wxa7d+9WaGioVq5cqZs3b96LXQIAAABArjMscF24cEGhoaEKCwvT2bNnJf1vBkM3Nzc1bNjQqF0DAAAAQJ7g1MAVHR2tFStWaOnSpdq1a5fMZrMlZBUoUECNGzdW+/bt1b59e5UoUcKZuwYAAACAPMfhwJWamqpNmzYpNDRU69atU0JCgqT/nc0ymUwaMGCAhgwZomLFijm6OwAAAABwGTkOXIcOHVJoaKiWLVumq1evSvpfyKpWrZq6d++uTz75RJIUFBRE2AIAAADwj2NX4Lpy5YrCwsK0dOlSHTt2TNL/QtZ9992nLl26qFu3bgoMDJQkS+ACAAAAgH8imwPXoEGDtG3bNqWmplpClq+vr9q3b69u3bqpadOmMplMhhUKAAAAAK7G5sC1ZcuW208oUEAtWrRQ165d1a5dOxUsWNCw4gAAAADAldl142OTySR3d3cVLFhQBQsWlLu7u1F1AQAAAIDLs/kMV+PGjbVjxw4lJCRo1apVWrVqlQoXLqyOHTuqW7du3FcLAAAAAO5i8xmu7777TmvWrNGIESNUsWJFmc1m3bhxQ4sWLVLfvn3Vrl07ffrppzp58qSR9QIAAACAy7BrSGGZMmU0bNgwrVy5UgsXLtQTTzyhwoULy2w26/z58/riiy/UuXNnPfroo5ozZ45RNQMAAACAS8jxfbjq1q2runXr6q233tL69esVGhqqDRs2KDk5Wfv379eBAwcssxbu3r1bTZs2VfHixZ1WOAAAAADkdSZz2hzvTnD9+nVFREQoNDRU+/fvv72D/4Yud3d3NWzYUO3bt1f79u1VsmRJZ+0WAACnioy8pZz+62gySf7+hRzaBjJHf41Ff41Ff43njB6nbcNpNTkzcN3pxIkTCg0NVUREhC5dunR7Z/8NXyaTSfXr19e8efOM2DUAAA4hcOVd9NdY9NdY9Nd4eTFw2XUNlz2qVKmiV199VWvXrtW3336rHj16yNvbW2azWampqdq1a5dRuwYAAACAPMGwwJXGZDKpWbNmmjhxorZs2aKPPvpIzZo1s5ztAgAAAID8KseTZuSEt7e3unfvru7du+vy5cv3ctcAAAAAcM8ZfoYrM6VKlcqtXQMAAADAPZFrgQsAAAAA8jsCFwAAAAAYhMAFAAAAAAYhcAEAAACAQQhcAAAAAGAQAhcAAAAAGITABQAAAAAGIXABAAAAgEEIXAAAAABgkALO3NitW7e0Zs2aLNfx8vJSx44dnblbAAAAAMiT7Apc69at08cffyxJGjZsmLp06ZLu8b/++kujR4+WyWTKcjteXl5q3bq1fZUCAAAAgIuxa0jhlClTdPLkSRUvXjxD2Lqb2WzO9Gvy5MkOFQ0AAAAArsDmM1z79u3T0aNHZTKZ9Morr2S7/oQJEzIsu3HjhiZOnKhjx45p7969qlu3rn3VAgAAAIALsTlwrV69WpJUu3ZtNWjQINv1e/bsaXX5+vXrtX37dq1Zs4bABQAAACBfs3lI4f79+2UymdSqVSuHdti2bVuZzWb9+eefDm0HAAAAAPI6mwPXiRMnJN0+w+WIqlWrSpJOnz7t0HYAAAAAIK+zeUjhzZs3JUklSpTIdB0vLy8FBARkOUth0aJFJd2+ngsAAAAA8jObA1dKSoqk27MPZqZ8+fJaunRplttJTk5Otz0AAAAAyK9sHlKYdmbq2rVrDu0wKioq3fYAAAAAIL+yOXCVLl1akrRnzx6Hdpg2WUba9gAAAAAgv7I5cDVu3Fhms1m//PJLjndmNpu1cuVKmUwmNWrUKMfbAQAAAABXYHPgatu2rSTp1KlTWrBgQY529uOPP1pmO0zbHgAAAADkVzYHrgYNGljOcn3wwQd2n+n69ddfNW7cOMvZrYYNG9pdLAAAAAC4EpsDlyS988478vX1VXJysl566SW99NJL2rlzZ5bP2bVrl1555RW9+OKLSk5Olo+Pj8aOHetQ0QAAAADgCmyeFl66fdPiKVOmaOTIkYqLi9Mvv/yiX375Rb6+vgoMDFTx4sXl7e2tuLg4RUVF6fDhw4qJiZF0+/otLy8vTZ48WdWqVTPkxQAAAABAXmJX4JKkli1bauHChXr11Vd17NgxSVJ0dLTVM1133rOrWrVq+vjjjxUYGOhAuQAAAADgOuwOXJJUvXp1hYeH69dff9XSpUv1xx9/6Pr16xnWK1q0qBo3bqzu3burXbt2jtYKAAAAAC4lR4ErTbt27SxB6vLly7p27ZpiYmLk6+urYsWKqVSpUk4pEgAAAABckUOB606lSpUiYAEAAADAHeyapRAAAAAAYDsCFwAAAAAYxK4hhZcuXVJERIQkKSgoSA8++GC6x48dO6Zu3bpluQ0vLy8tX75c999/v52lAgAAAIBrsesM15QpUzR58mR98803qlSpktV1zGZzll9xcXH697//7ZTiAQAAACAvs/kMV2RkpJYtWyZJGjlypO67775M1zWZTGrUqFGG5bGxsTpw4IBWrlypMWPGqESJEjkoGQAAAABcg82B65dfflFKSor8/f31yCOPZLv+d999l2FZSkqK2rdvr4sXL2r58uXq16+ffdUCAAAAgAuxeUjhjh07ZDKZFBwcLE9PzxztzN3dXZ07d5bZbNauXbtytA0AAAAAcBU2B67Dhw9LUoaJMuwVFBQkSTp06JBD2wEAAACAvM7mwHX16lVJcvjmxmnXfl27ds2h7QAAAABAXmfzNVyxsbGSJF9f30zXqVixokJDQ7PcTtpwxJiYGFt3DQAAAAAuyebA5evrq5s3b+rGjRuZruPp6anAwMAst5P2/KyCGwAAAADkBzYPKUybwv3kyZMO7fDUqVOSpOLFizu0HQAAAADI62wOXLVr15bZbNbmzZsd2uHGjRtlMplUu3Zth7YDAAAAAHmdzYHroYcekiStX79eR48ezdHOjh07pvXr10uSmjdvnqNtAAAAAICrsDlwderUSSVLlpTZbNZLL72kmzdv2rWjmzdv6qWXXlJqaqpKlCihkJAQu4sFAAAAAFdic+AqWLCgRowYIbPZrFOnTqlHjx7asGGDTc/dtGmTevbsqRMnTshkMmnkyJE5vnkyAAAAALgKm2cplKTevXtr3759WrRokS5duqShQ4eqQoUKatGihWrWrKnixYvL29tbcXFxioqK0qFDh7Rp0yadPXvWso1HH31UvXv3dvoLAQAAAIC8xq7AJUnvvfeeSpQooS+++EJms1lnz57VDz/8kOVzzGazTCaThgwZopEjR+a4WAAAAABwJTYPKUyTNiRw4cKFatu2rUwmk8xmc6ZfJpNJbdu21cKFC/Xyyy/Lzc3uXQIAAACAS7L7DFeaoKAgzZgxQ7du3dKOHTt05MgRXbt2TTExMfL19VXRokUVGBioRo0aqVChQs6sGQAAAABcQo4DV5pChQqpTZs2atOmjTPqAQAAAIB8g/F9AAAAAGCQXAlc0dHRio6Ozo1dAwAAAMA9Y3Pg6tmzp3r16qVz585ZffzWrVsKDQ1VaGholts5duyYGjVqpCZNmthVKAAAAAC4Gpuv4Tp06JBMJpMSEhKsPv7XX39p9OjRcnNzU48ePbLdntlstrlIAAAAAHBFDk+acTeCFADgn8psNutmVIySYxMVE5eoQsV8ZTKZcrssAHmA2WxW9PV4pSakKiY2Qb5FvPj78A/h9MAFAMA/SUpyinauP6wNS3bp6J6zirkZb3nMt7CXqteroNa9GqpBqwC5F3DPxUoB3GspKanat+W0tq04rFMHLiv21v9GivkUKqhKtUqpWUgN1X6ootzdmcsuvyJwAQCQQwf/OKUv3l6syIvX5eZmUmpq+lEeMTfjtXfzMe3eeFT+ZYrq+fG9VKNxpVyqFrAuPjZRO9Yd17aVRxT11y0lJqbI09NdxUsXUtOOAWrUpqq8fDxzu0yXc2zPBX0/cZ2iLkfL5GaS+a6/D7G3EnTw93M6sO2sipfyU9/RwapWr0wuVZs/xMUkKuryLV09d1NxiUkqXqqQvH1z/9glcAEAYCez2azFX6zTz5+tk8nt9pCgu8NWmrTlUX/d0PvPzNIjw9qo19A2DCVCrkuIS9KyOTu0ffVRJSempH8sPlmXzlzTki+2KXzWH3rw4erqPKCRCnp75FK1rsNsNmvF3J1aMXuH5e/D3WHLsu5/l1+7EqOpL4ep08BG6tS/IX8f7HT26BVtDj+onWuPKyU51bLcvYCbGgZXVYuuNVWheslcq49zl7gnzp49q7Fjx6pt27aqU6eOGjRooCeeeEJz5sxRfPz/ht8EBwcrICDA6tegQYMs602bNk0BAQGKiorKdJ/bt29XQECAVq5c6VDtjz76qAICAvTDDz9YfXzx4sXp6qxTp45atGihQYMGae7cuVZvgZBd/cHBwRoyZEi6ZZn1JSAgQGPHjk237tq1a9W3b181a9ZMdevWVdu2bTVy5Eht3LjRss758+cVEBCgb775JsP+L168qLFjxyo4OFi1a9dWs2bNNGzYMO3cuTPDuml9DggI0P79+zM8Pnr0aNWvX9/q67TWk8DAQF26dCnD49HR0QoKClJAQIDGjRuX4fFr167po48+UocOHVSnTh01adJEgwYN0rp16zKsm9Vrv5Otx2OaQ4cOadSoUWrVqpVq166tJk2aaODAgfr555+VkpJi076/+eYbBQQE6Pz585Zl/fr1U5cuXbKsNa1/mX1duXLFsm5UVJTGjx+vjh07KigoSM2aNdOjjz6qjz/+WDExMVnuJ42tv9NpUlJS1KJFCwUEBGjDhg1ZvoacHAPOeE32SAtbUuZvpO6WFrx+/mydFn+R8bgE7qVb1+I07fUI/bbicIawZfHfQzs5MUW/rTis6a8v061rcfeuSBeVFrYk2/8+pK23YvYOrZib8d9aWBcfm6gv316pyS+GaseaY+nCliSlJKdqx5pjmvxiqGb+30rFxybmSp2c4YLh1q9fr5EjR8rT01Pdu3dX9erVlZSUpJ07d+rjjz/W8ePH9f7771vWr1Gjhp5++ukM27nvvvvuZdmSpNOnT2vfvn0qW7aswsPD9dRTT2W67ogRI1SuXDklJycrMjJSv//+uz788EPNnj1bn332mQIDAx2up3nz5urevXuG5ZUq/W+I0jfffKNJkyapSZMmGjJkiLy8vHTmzBlt3bpVy5cvV8uWLbPcx86dO/Xcc89Jknr37q0qVaooMjJSS5YsUZ8+ffTWW2+pX79+Vp87ffp0ffHFFw68QsnT01MRERF69tln0y1ftWpVps85efKkBg4cqKioKPXq1Ut16tTRzZs3FR4erqFDh+qZZ57RG2+8kaN6bD0ef/zxR73zzjsqUaKEunfvrooVKyomJkbbtm3TW2+9pStXrmjo0KE5qsEe7777rnx8fDIsL1y4sCTp+vXreuSRRxQdHa1HHnlElStX1vXr13XkyBHNnz9fTz75pHx9fbPch72/05K0bds2XblyxfK71KpVq0y3b+8x4IzXZI+Df5yyhK2c+vmzdarZqBLDC5ErEuKS9OXYlbp89rpdgeCvs9f05diVenFSF850ZeLYnguWsJVTK2bvULW6ZRhemI342ERNGxWhi6duf3idmpLJKIP/Lj+047ymjYrQi590uedDZAlcMNS5c+f08ssvq0yZMpozZ066N6l9+vTRmTNntH79+nTPKVWqlNVQkRvCwsJUokQJjR49WiNGjND58+dVrlw5q+u2bNlSderUsXw/ZMgQbd26VUOHDtWwYcO0fPlyeXl5OVTPAw88kGVvkpOT9dlnn6l58+aaNWtWhsevXr2a5fZv3LihkSNHysvLS/Pnz1eFChUsjz399NMaNGiQPvzwQ9WqVUsNGjRI99waNWpo3bp1OnDggGrVqmXnK/ufVq1aadmyZRnebEdERKh169b65Zdf0i1PSkrSyJEjdfPmTc2bN09169a1PDZw4ECNGjVKs2bNUp06dRQSEmJ3PbYcj3v27NE777yjevXqaebMmfLz80tXw759+3Ts2DG7950THTp0UPHixTN9/KefftLFixc1f/78DD/D6OhoeXhk/SYqJ7/T0u3fpVq1aqlHjx6aMmWKYmNjrQZDyf5jwNHXZI+U5BR98fZiq9dj2MPNzaTP316sKcteYiIN3HPL5uzQX3aErTSpqWb9dfa6ls3ZoV5DmxlUnetKSUnV9xPXOfz3weRm0vcT12rsvKeYSCMLcz5cq4unouz60ODiqSjNnbBWz73f0eDq0uOnCEN9/fXXio2N1QcffGD1DFXFihU1YMCAXKjMNhEREerQoYNat26tQoUKKSIiwq7npw3Fu3DhgsLCwgyq8n+uXbum6OjoDG8605QoUSLL5y9cuFBXrlzRa6+9li5sSZKXl5cmTpwok8mkGTNmZHhu3759VaRIEU2bNi3nL0BSly5ddOjQIZ04ccKy7MqVK9q2bZvVYXWrVq3S0aNH9eyzz6YLW5Lk7u6ucePGqXDhwg7XlZXp06fLZDLpk08+SRe20tSpU0e9evUybP/2OHv2rNzd3VWvXr0Mj/n5+algwYJZPj8nv9Px8fFavXq1QkJC1KlTJ8XHx2vNmjWZ7sPeY8DR12SPnesPK/Ki/W9U75aaalbkxevateGIkyoDbBMfm6jtq4/m+Bg2p5r1++qjuTY0Ky/bt+W0oi5HO/z3wZxqVtTlaO3/7YyTKst/zh65okN/nLO71+ZUsw7+fk5nj17JfmUnsvsM15gxY+Tt7Z1heVzc/8b09u/fP9Pn37ke8r9169apfPnymQYAa5KTk61e2+Tj4+PwGSJ77N27V2fOnNGHH34oT09PPfzww5Yhavbo3r27Jk+erM2bN+uxxx5L99iNGzesPic1NdXq8oSEBKu98fPzk6enp0qUKCEvLy/LNVxFixa1q9a1a9eqYMGCmZ4JKl++vBo2bKjt27crPj4+3c/Dz89PAwYM0NSpUx06y9W4cWOVLl1aERERGjlypCRp+fLl8vHxUevWra3WLCnTG64XKlRIbdu21ZIlS3TmzBlVrFjRrnqyOx7j4uK0bds2NWrUSGXK5P7wD2vHVIECBSxDCsuWLauUlBQtXbpUPXv2tHv7OfmdXrt2rWJjY9W5c2eVLFlSTZo0UXh4uLp27Wp1fXuPAUdfkz02LNlldTbCnHBzM2n9kl1q3LamEyoDbLNj3fHMr9myUVJiinauO6HmnWs4qar8YduKww6f3UpjcjNp24rDqvsvhh1bsznioNzcTZkOI8yKm7tJW8IPqsKrmQ9tdza7A5e1i+LTpM2o8scff+S8IuQb0dHRunz5stq2bWvX8zZv3qxmzTIOVXj11Vct1xbdC2FhYbr//vvVsGFDSVLnzp31888/69ChQ6pRw/Z/ZEqXLq1ChQrp3LlzGR7r2DHzU9oBAQEZlv3000/66aefMiyfPHmyOnfuLDc3Nw0aNEgzZsxQmzZt1KhRIzVs2FD/+te/bApAJ06cUKVKleTpmfnY5oCAAP3+++86c+ZMhhr79++vOXPmaPr06fr888+z3V9mQkJCtGzZMsub7fDwcD388MNW6zpx4oQKFSqksmXLZrq9tOvnTpw4YXfgyu54PHPmjJKSklS9enW7tmsUa8dUpUqVLJPHPPLII5o9e7ZGjx6tmTNnqkmTJmrcuLFatWqlQoUKZbntnP5Oh4WFqX79+rr//vsl3f5deu+99xQVFZXp8Ed7jgFHXpM9zGazju4565SwJd0+y3V0z1mZzWZmJLNTQnxSbpfgsratdMJZVZO0deVhAtcdzGazTh247JSwJd0+E3Ny/1/8fbAiLiZRO9cez1HYkm5f07Vj7XH1HNpMXvdoyni7ApfZ7JyDCP8MabPz2Xuxet26dfXSSy9lWG7vG2VHJCcna/ny5erRo4flD13Tpk1VokQJhYWF2RW4pNtnQ6zNlDZt2jSrQ9Bee+01q9tp27at+vbtm2H5nW/2R4wYocqVK+uHH37Q5s2btXHjRk2ZMkU1a9bUJ598oipVqmRaZ0xMTLY/r7THrc2+WKhQIfXv31/Tpk3TwYMHVbNmzj6579q1q2bNmqU///xTRYoU0b59+/TKK68YUnN2sjsec3qcG8XaMXXnqAR/f38tXbpUM2bM0K+//qoFCxZowYIF8vDw0PPPP69hw4Zl+o97Tl7rtWvXtHnzZo0ZM8ayrH379ho3bpxWrFihPn36WH2ePceAI68pMTFRiYnph0b5+fnJ2uo3o2LS3dTYGWJuxCn6eqwKF88bx48rWDJ7uw7uPK+aDcup58AHc7sclxP11y3HN2LW7XtLkQMsoq/Hp7upsTPE3kpQzM14FSqacWTZP1nU5VsZZiO0V0pyqqIu31LZKtYvtXD2sW1z4Dp8+LBz94x8L+1Nn71TMhcrVkwPPfSQESXZbMuWLYqKilJQUJDOnPnfGOoHH3xQy5Yt02uvvSY3N9svgYyNjbV6/VSjRo2sfsKf2TUnpUuXtqk3Xbp0UZcuXRQdHa29e/dq8eLFioiI0NChQxUREZHp9n19fbP9eaU9ntmb7gEDBmjOnDmaNm1ajs9y1axZU5UrV1ZERIQKFy6skiVLqmnTppnWfO3aNYdqzkp2x2NOj/Ps5PQTzcyOqTvdd999eu+99/Tuu+/q9OnT2rx5s7766itNnTpV9913n3r37m31eTl5rcuXL1dSUpJq1KiR7ncpKChI4eHhmQYue44BR17Tl19+qenTp1u+9/X11a5du1SiRMYzY8kGXbPi4+Uhf3/nnYnLzxLik3Rw5+1bJhzceV79h7dSQS9my7NHooPDCS3bSUjmuL1DaoJjASAzvt4F6fNdrp676ZTteHneu7+9zFIIw/j5+em+++67Z7OzOVPaBBfWzmxI0u+//57lm787/fXXX7p161aGSSjuBT8/PzVv3lzNmzeXh4eHlixZor1796pJkyZW169SpYoOHjyoxMTETIcVHjlyRB4eHnrggQesPl6oUCENGDDAcpYrp7p06aL58+fL19dXnTp1yjTgVqlSRYcOHdLFixczvYbqyJHbQ2iqVq2a43oyU7FiRRUoUEBHjx61af20sGvtXlXS/65zdeZED5kxmUyqVKmSKlWqpNatW6t9+/YKCwvLMnDZ+zsdHh4uSXryySetPn7u3DmVL1/e6mO2HgN3svc1DRkyxOq0/1ev3tLdgzpi4owJXLHxSYqMdMJZh3+Img3LWc5w3YqO161o5551zO88Pd2VEJ/s+HYKFuC4vUNMrHPPblm2G5dAn+8Sl+icIcXxiZn/7TWZZPWDt5zKtVkKV6xYkVu7xj3Upk0bnT17Vrt3787tUmwWGxurtWvXKiQkRJ9++mmGr5IlS1reRNpi6dKlkqQWLVoYVbJNateuLUnpboB7t9atWyshISHT38/z589r586devDBB7OcwGTAgAEqXLhwujMH9uratauuXLmi06dPZzq5QlrNkhQaGmr18ejoaK1Zs0aVK1c2ZFiqt7e3mjZtqh07dli9We/dihcvLm9vb506dcrq46dOnZK3t7eKFSvm7FKzVL58eRUuXDjL40Oy73f63Llz2r17t/r27Zvh92jKlCny8PDI8nfJ1mPAkdfk6ekpPz+/dF+SZDZn/CpUzFe+hZ07cY9vEW/5FfWxuj++rH/1HPigJn7bRz0HPpjrtbjiV/HSTngTaZKKl/LL9deSl758i3jJp5BzPyjzKVRQvoW9cv215bWv4qUKyb2AYxHGvYCbipcqlOV+nMnmat977z2n7XTJkiUaNWqU07aHvGvw4MHy8fHR22+/rcjIyAyPnz17VnPmzMmFyjK3evVqxcbGqk+fPurYsWOGrzZt2mjVqlUZrvuwZuvWrfrss89Urlw5devWzfDa4+LiMn0jvHHjRknpb5J8t8cff1wlSpTQxx9/nGGSj4SEBI0ZM0Zms1kvvPBClnWkneVas2aNDh06ZOeruK1ChQp688039eqrryooKCjT9Tp06KCqVavqq6++0r59+9I9lpqaqnfeeUc3btzQ8OHDc1SHLV544QWZzWa9/vrrVofb7d+/X0uWLJF0e6r65s2ba926dbp48WK69S5evKh169apefPmcnc35t5Me/fuVWxsbIblf/75p65fv57l8SHZ9zudFqYGDx6c4fcoJCTEMlthZmw9Bhx9TbYymUyqXq+C3NycM7jfze329rgg3n4MI8y5ph0zTshkN7PUrGOg49vJR0wmkyrVKiWTk/4+mNxMqly7NH8frPD29VTD4Kpyc89Zb9zcTWoUXPWeTZgh2TGkcP78+SpYsKBGjx7t0A4XLFjg1PCGvK1ChQr65JNP9PLLLyskJETdu3dX9erVlZiYqN27d2vlypUZ7k90+fJly1mhO/n6+qpdu3bpls2ePTvDmRY3N7d0U7evWrVKJ0+ezLC9nj17WmZNu1N4eLiKFi2q+vXrW31NwcHBWrRokdavX6/27dtblm/cuFEnT55USkqKIiMjtX37dm3ZskVlypTR559/7pQhYqdPn7baG39/fzVv3lxxcXF64oknVK9ePf3rX/9S6dKldevWLf3666/asWOH2rVrl+VEFsWKFdPUqVP13HPPqWfPnurdu7eqVKmiyMhIy7Tqb731lk1Tgvfv31+zZ8/W4cOHM73BbXZsuUebp6enpk6dqgEDBuipp55Sr169VLt2bd26dUsRERE6cOCAnnnmGXXu3DnDc7du3aqEhIzDQNq1a2eZiMSW47FBgwYaO3as3nvvPXXq1Endu3dXxYoVFRMTo99//11r165NNzz1lVde0WOPPaaePXvq8ccfV9myZXXhwgUtXLhQJpPJ6uQQUVFR+uyzzzIsvzvM//LLL1b73bx5c8vkEuHh4WrXrp1q164tDw8PnThxQj///LMKFiyY7W0P7PmdDg8PV40aNaz+nkm3f5fef//9LG8jYMsx4Ohrskerng20e6Ntw0ezk5pqVuuetk+vDzhDozZVFT7rD4emhvfwdFfDNplPwPRP1bRToA5sO+uUbZlTzWraiVCbmRZda+r3VTn7W5yaYlaLbvf2dhx2XcM1Z84ceXh46NVXX83Rzr799ltNmjRJZrM5y2mnkb+0bdtWYWFh+uabb7RmzRrNnz9fnp6eCggI0OjRozPcm+rQoUN6/fXXM2ynbNmyGQLXl19+mWE9d3f3dG+wli1bZrWuJk2aZHgjePXqVW3dulWdO3fO9AxDs2bN5O3trbCwsHSBa+rUqZIkDw8PFS1aVNWrV9ebb76pXr16WZ2JMCe2bNmiLVu2WH0tzZs3V+HChTV+/HitX79eixcv1pUrV+Tu7q5KlSrp9ddfV79+/bLdR6NGjRQWFqYvv/xSK1eu1JUrV+Tn56f69evrgw8+UKNGjWyqtXDhwhowYIBDwwptVaVKFYWFhWnmzJlau3atFi9eLC8vL9WuXVuff/65goODrT5v06ZN2rRpU4blZcuWtQQuW4/HJ554QnXq1NGsWbMUGhqqa9euycfHRzVr1tSECRPShaIqVapo0aJFmj59un766SfduHFDRYoUUfPmzfXCCy9YnUny6tWr+vTTTzMsb9asWbptv/vuu1Zf69y5c+Xv76/HH39cXl5e2rZtm9auXavo6GgVK1ZMzZs315AhQ2yaWdKW3+kDBw7o5MmTGjZsWKbbadOmjd5//32FhYXl+L5tkpzymmzVsHWg/MsU1dW/bjg0/bObm0nFSxdRg1ZOONsA2MHLx1MPPlxdv604nKNj2ORmUpOHq8vLh/dxd6vT/AEVL+Wna1diHPr7YHIzqVhJX9V+qKITq8tfKlQvqRqNy+vwzvN29drkZlKNRuVUvlpJA6uzsl+zjXO9h4SE6OTJkzKZTHr++ec1YsQIu3b02Wefadq0aTKbzfLy8tLUqVPVsmXLHBUNAICRIiMzTpqR5uAfpzT+mVkO7+P/Zj2jGo25qam9TCbJ379Qlj8jZC0hLknTXo/Q5bPX7bqvnJubSaUrFNPwSZ1V0JthndYc23NBU1+2/TrvzIyY0k3V6lmfCAq3xccmatqoCF08FWVT6DK5mVSmUnG9+EmXbD8wSPs74yw2X8M1Z84cVaxYUWazWZ9//rld0z1PnjzZErZ8fHz0xRdfELYAAC6pZuNKemRYG4e28ciwNoQt5JqC3h4aMq6jSlUoavM1R6b/hq3nxnUgbGWhWr2y6jTQtpEgmek0sBFhywZePp568ZMuqtGonCRlek1X2vIajcppxL+75srZWZvPcEm3r2Xo27evzp07J5PJpFdffVWDBw/O8jkffPCBvv/+e5nNZhUqVEhffvmlTdd/AACQW7I7e2I2m7X4i3X6+bN1cnMz2XSWIG29R4a1Ua+hbbgYPoc4w+U8CXFJWjZnh35ffVRJiSmSSdKdPf3v9x6e7mrycHV1HtCIsGUDs9msFXN3asXsHTK5mWw++2JONavTwEbq1L8hfx/sdPboFW0JP6gda4+nuymyewE3NQququZda6pCdduHETr7DJddgUu6PYtW3759dfHiRZlMJo0ZM0b9+/e3uu7YsWP1448/ymw2q2jRovr6668tU1MDAJBX2fpm/tAfp/T524sVefF6psErbbl/maJ6fnwvzmw5iMDlfPGxidq57oS2rjysqMvRSkxIlmfBAipeyk/NOgaqYZsqXLOVA8f2XNT3E9cq6nJ0psErbXnxUn7qOzqYM1sOiotJ1LXLt+Tl6aH4xCQVL1UoR7MR5nrgkm7fW6V///66dOmSTCaTxo4dm+6mlqmpqRo9erTCw8NlNpvl7++vb775RgEBXBwMAMj77Hkzn5Kcol0bjmj94p06uvecYm7EWR7zLeKt6vUqqHXPBmrQKkDuBYyZ7v+fhMBlLPrrXCkpqdr/2xltXX5Ipw5cVuyt/82M61OooCrXLq2mIYGq3ayi3N1z7fa4+YozjuE8Ebgk6cyZM+rXr5/+/vtvubm5ady4cXr00UeVlJSkV199VatXr5bZbFapUqX07bffqnLlyk4rGgAAI+X0H2qz2azo67Hy8fJQbHyS/Ir6MDTIyQgExqK/xjGbzYq5GS9f74KKiUuQb2Ev/j4YIC8GLrumhb9TxYoVNXv2bPXv31+RkZEaO3asUlJStG7dOq1fv17S7WmTZ8+erfLlyzurXgAA8iyTyaTCxX15wwogA5PJpEJFvfn78A/k0LnLypUra/bs2SpWrJhSU1P17rvvasOGDZJuB7J58+YRtgAAAAD8Yzk8WLRq1aqaPXu2ihYtKrPZLLPZrGrVqmnevHkqXbq0M2oEAAAAAJdk85DC6dOnZ/l4o0aN9Ouvv8rNzU1NmjTRggULst3m8OHDbd09AAAAALgcmyfNCAwMdPqFfYcOHXLq9gAAcAZnXGzNNRrGoL/Gor/Gor/Gc/lJM3I4oaFVzMoCAAAAIL+zOXDNnTvXyDoAAAAAIN+xOXA1adLEyDoAAAAAIN/hltYAAAAAYBACFwAAAAAYhMAFAAAAAAYhcAEAAACAQQhcAAAAAGAQAhcAAAAAGITABQAAAAAGIXABAAAAgEEIXAAAAABgEAIXAAAAABiEwAUAAAAABiFwAQAAAIBBCFwAAAAAYBACFwAAAAAYhMAFAAAAAAYhcAEAAACAQQhcAAAAAGAQAhcAAAAAGITABQAAAAAGIXABAAAAgEEIXAAAAABgEAIXAAAAABiEwAUAAAAABiFwAQAAAIBBCFwAAAAAYBACFwAAAAAYhMAFAAAAAAYhcAEAAACAQQhcAAAAAGAQAhcAAAAAGITABQAAAAAGIXABAAAAgEEIXAAAAABgEAIXAAAAABiEwAUAAAAABiFwAQAAAIBBCFwAAAAAYBACFwAAAAAYhMAFAAAAAAYhcAEAAACAQQhcAAAAAGAQAhcAAAAAGITABQAAAAAGIXABAAAAgEEIXAAAAABgEAIXAAAAABiEwAUAAAAABiFwAQAAAIBBCFwAAAAAYBACFwAAAAAYhMAFAAAAAAYhcAEAAACAQQhcAAAAAGAQAhcAAAAAGITABQAAAAAGIXABAAAAgEEIXAAAAABgEAIXAAAAABiEwAUAAAAABiFwAQAAAIBBCFwAAAAAYBACFwAAAAAYhMAFAAAAAAYhcAEAAACAQQhcAAAAAGAQAhcAAAAAGITABQAAAAAGIXABAAAAgEEIXAAAAABgEAIXAAAAABiEwAUAAAAABiFwAQAAAIBBCFwAAAAAYBACFwAAAAAYhMAFAAAAAAYhcAEAAACAQQhcAAAAAGAQAhcAAAAAGITABQAAAAAGIXABAAAAgEEIXAAAAABgEAIXAAAAABiEwAUAAAAABiFwAQAAAIBBCFwAAAAAYBACFwAAAAAYhMAFAAAAAAYhcAEAAACAQQhcAAAAAGAQAhcAAAAAGITABQAAAAAGIXABAAAAgEEIXAAAAABgEAIXAAAAABiEwAUAAAAABiFwAQAAAIBBCFwAAAAAYBACFwAAAAAYhMAFAAAAAAYhcAEAAACAQQhcAAAAAGAQAhcAAAAAGITABQAAAAAGIXABAAAAgEEIXAAAAABgEAIXAAAAABiEwAUAAAAABiFwAQAAAIBBCFwAAAAAYBACFwAAAAAYhMAFAAAAAAYhcAEAAACAQQhcAAAAAGAQAhcAAAAAGITABQAAAAAGIXABAAAAgEEIXAAAAABgEAIXAAAAABiEwAUAAAAABiFwAQAAAIBBCFwAAAAAYBACFwAAAAAYhMAFAAAAAAYhcAEAAACAQQhcAAAAAGAQAhcAAAAAGITABQAAAAAGIXABAAAAgEEIXAAAAABgEAIXAAAAABiEwAUAAAAABiFwAQAAAIBBCFwAAAAAYBACFwAAAAAYhMAFAAAAAAYhcAEAAACAQQhcAAAAAGAQAhcAAAAAGITABQAAAAAGIXABAAAAgEEIXAAAAABgEAIXAAAAABiEwAUAAAAABiFwAQAAAIBBCFwAAAAAYBACFwAAAAAYhMAFAAAAAAYhcAEAAACAQQhcAAAAAGAQAhcAAAAAGITABQAAAAAGIXABAAAAgEEIXAAAAABgEAIXAAAAABiEwAUAAAAABiFwAQAAAIBBCFwAAAAAYBACFwAAAAAYhMAFAAAAAAYhcAEAAACAQQhcAAAAAGAQAhcAAAAAGITABQAAAAAGIXABAAAAgEEIXAAAAABgEAIXAAAAABiEwAUAAAAABiFwAQAAAIBBCFwAAAAAYBACFwAAAAAYhMAFAAAAAAYhcAEAAACAQQhcAAAAAGAQAhcAAAAAGITABQAAAAAGIXABAAAAgEEIXAAAAABgEAIXAAAAABiEwAUAAAAABiFwAQAAAIBBCFwAAAAAYBACFwAAAAAYhMAFAAAAAAYhcAEAAACAQQhcAAAAAGAQAhcAAAAAGITABQAAAAAGKZDbBQAAACD3xEUnaMvyfVq/eI8iL1xXQnySCnp5yL9sUbXuVU/NQ+rI269gbpcJuCyT2Ww253YRAADkJZGRt5TTfx1NJsnfv5BD20Dm6K/zxMcm6sdp67Rx6V4lJSZLd/fTdPs/Hp4F1LJ7XfV+sY28fDzveZ35Ccev8ZzR47RtOItLDilcvHixAgICVKdOHV2+fDnD4/369VOXLl2sPjclJUUtWrRQQECANmzYkGGb2X0FBwdLkkaPHq369etnWmP9+vU1evRoy/fbt29Pt50aNWqoWbNmGjFihE6cOJHpdjZs2KCAgAC1aNFCqampVtcJDg7WkCFDMt1GZtJe8759+9Itv3Xrlh599FHVqVNHGzdulCRNmzZNAQEBioqKsqyXXQ+y2sedzp8/n2XPZ86caVm3X79+6R4LCgpS165dNXv27Ez7Y01KSop+/vln9evXT02aNFHt2rUVHBysMWPGpKs1u+Niz549lnUDAgI0bty4LPeb1bGZlTuPn6VLl1pd54knnlBAQECG7QcHB2da/6BBgyzrpf2M077q1q2r1q1ba+jQofr555+VmJiYYZ/ZHQOPPvqoAgIC9MMPP1h9Ldl9Sfb/DNK+atasqSZNmqhXr14aP368jh8/blOv02R3jNy8eVMtWrTQE088IWufW+3Zs0eBgYEaMWKEza/3/PnzludfvHhRY8eOVXBwsGrXrq1mzZpp2LBh2rlzZ7r9ZPXzvfNr8eLFlufExsZqxowZ6tq1q+rWrauGDRvqqaeeUmhoaLrXsmXLFgUEBGj69OkZXt+5c+dUt25djRgxwrIsq+Nh9erVGjx4sB588EHVrl1bLVq00MiRI7V161bLOmnHxsqVKyXJ5r7NnDlTAQEB2rRpk9V9P/vss2rYsKHVfy+Af6obV6P1wTNzte6nXUpKsBK2pNvLzFJSQrLW/bRLHzwzVzeuRt/rUgGX59JDChMTEzVz5kz93//9n83P2bZtm65cuaKyZcsqPDxcrVq1kiQ1btxYkyZNSrfu22+/raCgID322GOWZb6+vg7V3K9fP9WpU0fJyck6cuSIFixYoO3btysiIkIlS5bMsH5YWJjKli2rCxcuaNu2bXrooYcc2n92oqOj9cwzz+jIkSOaPn26WrZsaej+7tSlSxer+6tZs2a670uXLq1XXnlFknTt2jVFRERowoQJunbtml5++eVs9xMfH6/hw4dr06ZNaty4sYYMGaIiRYrowoULWrFihZYsWaL169erdOnSlueMGDFC5cqVy7CtChUq2PsyHVKwYEFFRESoe/fu6ZafP39eu3fvVsGC1od81KhRQ08//XSG5ffdd1+GZe+++658fHyUmJioy5cva/PmzXrzzTc1Z84cffnll7r//vttqvX06dPat2+f5XftqaeekiRVqVIlw+/a5MmT5ePjo6FDh2a6PVt/Bs2bN1f37t1lNpsVHR2tw4cPKzQ0VPPnz9eoUaOs9uFuth4jb775pl5++WUtWrRIjz/+uOX5ycnJeuedd1SmTBlNnDhRbdu2Tbf9b7/9Vn/99ZfGjBmTbnnx4sUlSTt37tRzzz0nSerdu7eqVKmiyMhILVmyRH369NFbb72lfv36SZLefPNNxcTEWLaxceNGRUREaMyYMSpWrJhleYMGDSRJkZGRGjhwoE6cOKGQkBD17dtXCQkJWrVqld544w1t2LBBn3zyidzd3dW8eXN16dJFX375pTp37qxKlSpZtvfee++pQIECeuutt7Lspdls1ptvvqnFixerZs2aevrpp+Xv768rV65o9erVGjhwoObPn2+p7053HydLly7Vli1bMixv3ry5wsPD9d577ykiIkJeXl6Wx1asWKGNGzdq7NixKlWqVJa1Av8U8bGJ+njYfF08FanUVNtOAaSmmnXh5BV98sICvTWrP2e6kOfERSdo24rD2r3huBJik1TQx0P1W1VV006BuT4k1qUDV40aNbRo0SI999xzNv9DGhYWplq1aqlHjx6aMmWKYmNj5ePjo/Lly6t8+fLp1n333XdVvnz5DG9uHdGoUSN17NjR8n2lSpX07rvvKjQ0VM8++2y6dWNjY7V27Vq98sorWrx4scLDww0NXNHR0Ro0aJAOHTqk6dOnW8LovVKzZk2bel2oUKF06z355JPq1KmTvvvuO40YMULu7u5ZPn/SpEnatGmTxowZo4EDB6Z7bPjw4Zo9e3aG57Rs2VJ16tSx6XUYqVWrVlq7dq2ioqIsb84lKSIiQv7+/qpYsaJu3ryZ4XmlSpWy+Tju0KFDum0PHz5cYWFheuONNzRy5EgtWrTIpu2EhYWpRIkSGj16tEaMGKHz58+rXLly8vf3z1DLV199pWLFimVZo60/gwceeCDDdl599VU9//zzmjhxoipXrpztsW3rMRISEqIlS5bo3//+t9q2bSt/f39J0ty5c3X48GHNnDlTPj4+GepZvny5bt68afX13rhxQyNHjpSXl5fmz5+fLlA+/fTTGjRokD788EPVqlVLDRo0ULt27dI9PzIyUhEREWrXrp3VgPrGG2/oxIkTmj59erog2L9/f3300UeaNWuWatSoYQl8b775pjZt2qR33nlHc+fOlSQtW7ZMmzZt0ttvv53t395Zs2Zp8eLFGjBggMaMGSOTyWR57Pnnn1doaKgKFLD+T9Hd/dm7d6+2bNlitW/vv/++nnzySX322WeWD2Sio6P14Ycfql69enryySezrNNZrl+J1o41R5Ucn6ICXu5q3C5ARfwd+6AOcLYfp63TxZO2h600qSlmXThxRT9OW6d+b3QwqDrAfttWHNK8j9cqKTEl3dnao3suaPHnm9XntWA17VQj1+pzySGFaYYMGaLU1FR99dVXNq0fHx+v1atXKyQkRJ06dVJ8fLzWrFljcJVZa9SokaTbw3Putnr1asXHx6tjx44KCQnRqlWrlJCQYEgdMTExGjx4sA4cOKBp06apdevWhuzHCAULFlTt2rUVExOjq1evZrnuX3/9pYULF6p58+YZ3khLkru7uwYNGpTu7FZe0rZtW3l6elqGXKWJiIhQp06dsg2bOdWtWzf17t3b8obXFhEREerQoYNat26tQoUKKSIiwpDabFGsWDFNnjxZBQoU0Oeff57luvYeI++8844SExM1YcIESdKlS5c0bdo0hYSE5OhDi4ULF+rKlSt67bXXMpy98/Ly0sSJE2UymTRjxgy7t71nzx5t3rxZPXv2zHDWTbodTB944AF9/fXXio+PlySVKFFCo0aN0vbt27VkyRLdvHlTEyZMUJ06ddSnT58s9xcfH6+ZM2eqcuXKeuONN9KFrTQ9evRQUFCQ3a/lbvXq1dMTTzyhWbNmWYaP/uc//1FUVJTGjRsnNzdj/7lLTUnVj1M36s1e3+rnGVsU/u12/Txji8b0nKUfp25UaortQ54BI8VFJ2jj0r12h600qalmbVq6V3HRxrwfAey1bcUhzR6/WkkJKRmHxpqlpIQUzR6/WttWHMqV+iQXD1zlypVT9+7dtWjRIpvG5q9du1axsbHq3LmzSpYsqSZNmig8PPweVJq5CxcuSJIKFy6c4bHw8HA9+OCDKlmypDp37qyYmBitXbvW6TXExcXp2Wef1f79+/Xpp5+qTZs2Tt+HrXVERUVl+EpOTs72uRcuXJDJZLLaxztt3LhRycnJ6tatm121RUdHZ6jr2rVrdm3DGby8vBQcHKxly5ZZlh0+fFjHjh3L8tqw5ORkq71Ne1Nti7Sebd68Odt19+7dqzNnzqhz587y9PTUww8/7PDvmqM/gzJlyqhx48bau3evoqMzvwbB3mOkXLlyevHFFxUREaEtW7Zo/PjxKlCggN58802ba7vT2rVrVbBgQYWEhFh9vHz58mrYsKG2b99u189PktatWyfpdsixpkCBAurSpYtu3LihXbt2WZb37t1bDRo00EcffaT33nvP5hCzc+dOXb9+XV26dDHsw4A7vfrqqypWrJjGjh2r/fv364cfftAzzzxjuR7QSD/P2Kw1C/coNdUsc6pZKcmpMqealZpq1pqFe7T4M9s+qACMtmX5vtsTZDggMTFZvy3f76SKgJyLi07QvI9te2887+O1ufZBgUsHLun2kJSUlBSbznKFhYWpfv36lmtQOnfurC1btqSbCMJoMTExioqK0t9//61Nmzbpww8/lMlkUvv27dOtd/XqVW3dulWdO3eWdPvNYr169QwJiKNHj9aff/6p//znP1Y/9b5Xpk2bpmbNmmX42r8//R/1lJQUyxvukydPatKkSdq/f79atWqV7toNa9ImKLH3DdjAgQMz1HUvr2+7U9euXbVz505dunRJ0u3junz58qpXr16mz9m8ebPV3qYNEbNF9erVJVk/G3u3sLAw3X///WrYsKGk279rx48f16FDOf90yRk/g2rVqik1NTXd5BR3y8kxMmDAANWoUUOvvPKKfv31V40aNcrqNZm2OHHihCpVqiRPz8yvjwgICFBSUpLOnDlj17bTzvwEBgZmuk7aY3dO5mMymTRu3DhFR0crIiJCffv2zXBtpTU5/X3LKT8/P7399tvauXOnBg0apDJlyuiFF14wfL/Xr0Rr3Y97s1xn7aI9uhEZk+U6wL2w/ufd1ifIsNO6xbsd3wjgoG0rDt8eRmiDpMQUbVt52OCKrHPpa7ik25/2duvWzXItl7VJAKTbkyts3rw53UXq7du317hx47RixYpsh8Y4y92fehcvXlyTJk3KMKRm2bJlGYJYly5dNHHiRN24cUNFihRxWk2RkZHy9PS0eTIEozz++OPprm9LU7Vq1XTfnzx5Us2aNUu3LDg4WB988EG2+0g7s2Hv5Cdjx45NN2GAJMOHKGWmefPmKlKkiJYtW6ZBgwZp+fLl2V6fVbduXb300ksZllesWNHm/fr4+EhSugkarElOTtby5cvVo0cPyxCypk2bqkSJEgoLC1ONGjkbQ+2Mn4EtryEnx0iBAgU0btw49e7dW/Xq1Us30Y69YmJist132uNZnanLbNt3Pt+ebfv5+cnDw0NJSUlq3ry5TfvL6e+bIzp06KBWrVppw4YN+vjjj7P9ECYxMTHDDJx+fn6yMvoxUzvWHM32/av5v+u1eyLrmV2RvVU/7tH2Ncf0YNtqat+7Xm6X43IiL95wfCNmKfLCdbt+T3CbySS9/cJCvTzO+igG2Gf3huO2f4BglnZvOKHg3nWzXdXZx7bLBy5JGjZsmMLCwjRz5ky9/fbbVtdZvny5kpKSVKNGjXSfCgcFBSk8PNyQwGXteoUXXnhBjRo1UmxsrFavXq1ly5ZZfdMYFhamoKAgXb9+XdevX5d0e5KQpKQkrVy5Mt2MaI4aN26cJkyYoMGDB2vevHmqXLmy07Ztj4oVK9o0KUjZsmU1fvx4paam6uzZs/riiy907dq1TGfou5Ofn5+k7EPD3YKCgvLEpBmS5OHhoY4dOyoiIkJBQUG6dOmSunbtmuVzihUr5vCEK7GxsZKyf/OcdtY4KCgo3e/agw8+qGXLlum1117LUVh1xs/AltfgyDEiSbVq1bL6u28rX1/fbPdtS3DKbNtpz89s+G1m2x43bpxMJpPKli2rjz76SA899JA8PDyy3F9Oe+moOnXqaMOGDapdu3a263755Zfppr339fXVrl27VKKE7fdfSY5PkZubSSlZXBPj5mZSUlyyU+/r8k+1fc0xpSSnavuaY3rq+X/ldjkuJyE+ySnbSYxP4njOgbdfWChJmjJ2ucbPcN57uX+qhFj7jueEmMRcOW7zReC6+yyXNWlD8TKbqercuXMZZinMiqenpxITE2U2mzO8uTKbzUpISLA6JKh69eqWN77t2rVTXFyc/u///k8NGza0nGFKm05bUoahhmmvxZmBq0qVKvrqq680YMAAPfPMM5o/f36un+3Kio+PT7rw0KBBA/Xq1UtTpkzJNHCnSQuTR44cyfGZlryga9euWrBggaZNm6bAwMAMZwGNcPToUUnZT4UfFhYmSVbPqEnS77//rqZNmzq1NlsdO3ZM7u7uVmfvS5Pbx0iVKlV08OBBJSYmZjqs8MiRI/Lw8NADDzxg97Z//fVXHTlyRI0bN85021L6M8urVq3S2rVr9eabb+qBBx7Qc889p2+++SbLafyl9L28ezbFvGLIkCFWbxVw9artN8ws4OWe7QQEqSlmeXgXUGTkrZyUiTs82Laa5QwX/bRfQS8PxcdmvK+ivTy9POh/Drw8LkRTxi7Xy+NC6J8TFPTJ+oO/DOv7etrUd5NJdn3wlp18Ebik29dyhYWFWb2W69y5c9q9e7f69u2b4U1GamqqXn/9dYWHh2vYsGE2769s2bJKTk7W2bNnMwzLOnPmjFJSUlS2bNlstzNq1Cj9+uuv+vzzzy03zQ0PD5eHh4cmTZqU4UzAzp079d133+nixYsqU6aMzfVmJygoSJ999pmee+45Pf300/rhhx/STQ2elwUGBqpbt25asGCBnnnmmSz70rJlS7m7uys8PDzTiQNcQcOGDVWmTBn9/vvvGjVq1D3ZZ1qQ+te/Mv9EOe1WBiEhIerQIeOUwePHj1d4eHiuBK6LFy/qjz/+UL169SxnXqzJ7WOkdevW2r17t1asWGF1qOj58+e1c+dONWvWLNvhcta2/eWXXyo0NNRq4EpJSVF4eLiKFCliuS9WdHS0xo8fr1q1aqlv375yd3dXhw4d9Pnnn6tz585ZflDVsGFDy/DXoUOH3pOJM+zl6elpNdiazbI5cDVqW12LZ2zJclSLyXR7PVu3icy1711PTz3/L0VG2h6K8T/+ZYro/PErjm3EJPmXLUr/c2j8jMc5fp2kfquqOrrngm3DCk1S/VZVcqXvLj9pRpoKFSqoW7dulimV75R2dmvw4MHq2LFjuq+QkJAczVaYdrH+999/n+GxefPmpVsnu7rbt2+vJUuWWOoODw9Xw4YNFRISkqHewYMHS5IhU2w3a9ZMkydP1tmzZzV48GC7rw/JTYMHD1ZycrK+/fbbLNe7//771bt3b23evFnfffddhsdTU1M1a9Ys/fXXX0aV6hQmk0lvvfWWhg8f7tT7xGUmPDxcP/74o+rXr5/h+rk7rV69WrGxserTp0+GY7djx45q06aNVq1aleGaGaNdv35dr7zyilJSUrI9K5Pbx8jjjz+uEiVK6OOPP84wQUlCQoLGjBkjs9mco8kgGjRooIceekiLFy+2zFh4pylTpuj06dMaPHiwJcz95z//0ZUrV/Tee+9ZAtNbb70lNzc3vf/++1nuz9vbW4MHD9aJEyf0ySefyGzlX7mlS5fqzz//tPu15CVFS/qpTTbXBAQ/Xo/7cSFPaP1IfckJ16e06cX1iMh9TTsFysPTtg/zPDzd1bRj5pNGGSnfnOGSpKFDh2rp0qU6deqUqlWrZlkeHh6uGjVqZDpMLjg4WO+//74OHDigWrVq2bSvGjVqqHfv3po7d67OnDljGeL222+/acOGDerdu3eWM4HdadCgQVqxYoXmzJmjhx9+WGfOnMn0mrJSpUqpZs2aCg8PTzd88syZM/rss88yrF+zZk277qn18MMP6/3339ebb76p559/Xl9//XWW10YlJSVZ3W+RIkXSvYaff/5ZmzZtyrBe//79Lf9/8OBBLV26NMM6FSpUUP36Wf9hr1q1qlq1aqWffvpJw4YNU7FixTJdd/To0Tp37pzGjx+vVatWqU2bNipcuLAuXbqklStX6uTJk5bZIdNs3LhRJ0+ezLCtBg0apPuEf//+/Vb70aRJE8s916KioqyuU65cObumq2/Xrp3Nw7QuX75stbe+vr4ZtvHLL7/Ix8dHSUlJunz5sjZv3qxdu3YpMDBQn376aZb7CQ8PV9GiRTP9eQUHB2vRokVav3691eGyWbH1Z3D69GktXbpUZrNZMTExOnz4sFauXKnY2FiNHj3apg9CcnKMOEuxYsU0depUPffcc+rZs6d69+6tKlWqKDIyUkuWLNGZM2f01ltvWc5A2eujjz7SwIEDNWzYMHXp0kWNGjVSYmKiVq1apd9//10hISEaNGiQJFmmVu/Tp0+66+dKlSqlkSNHasKECfrll1+sns1MM3jwYB0/flyzZs3S9u3b1aFDB/n7+ysyMlK//vqr/vzzTy1YsCBHryUveeSFFjKZTFq7aI/Mun3NVmqKWSbT7bDV63nbJhoBjNY8pI4W/metkhJyPjW8p2cBPRSS/TWSgNG8/Qqqz2vBmj1+dbbr9nktWN5+2V/vb4R8FbgqVqyobt26acmSJZZlBw4c0MmTJ7McLtimTRu9//77CgsLszlwSbcvIq9evbp+/vlnTZ48WZJUqVIlvf3223ZNwlGnTh01adJE8+fPt9y4Nzg4ONP1g4ODNW3aNB0+fNgS6k6dOmX1zfCjjz5q902MH3nkEd24cUMfffSRRo4cme6C8rslJSVZ3W+FChXS9WD+/PlWn9+rVy/L/0dERFg9c9ezZ89sA5d0O7iuX79e33//vV588cVM1/P29tZXX32lxYsXKzQ0VJ999pni4+N133336cEHH9Qnn3yiUqVKpXvO1KlTrW5rwoQJ6d7s7927V3v3ZpweeuTIkZbAdfXqVas9a9asmd33B7PVoUOH9Prrr2dYXrZs2QyB691335V0+4bSxYoVU40aNfThhx+qa9euWU5VfuetDDIbOtasWTN5e3srLCzM7sBl689gy5Yt2rJli9zc3OTn56dy5cqpR48eevzxx22+1i0nx4gzNWrUSGFhYfryyy+1cuVKXblyRX5+fqpfv74++OADy7GUE/fdd59+/PFHffvtt1q5cqVWrVold3d3BQQEaOLEiZbZJVNSUjR27FiVKFHC6vV4/fr1U2hoqD788EO1aNEi0wk83NzcNGnSJLVt21aLFi3SrFmzFB0drWLFiqlx48Z67bXXbPr9zuvc3N306Iv/0sNPNtCONUeVFJcsD+8CatS2Ome2kKd4+xVUy+51te6nXTm6+bGbm0n/6l431964Andr2un29dbzPl57e4r4Ow9r0+0zW31eC7aslxtMZmtjPAAA+Adz5PoKk0ny9y/ENRoGob+Oi49N1AfPzNWFk1eUmmJ7E93cTSpbuaTemtVfXj6ZfwCHzHH8GicuOkHbVh7W7g0nlBCTqIK+nqrfqoqadgy0+wOCtJ+TsxC4AAC4C4Er76K/znHjarQ+eWGBLpy4YtOZLjc3k8pWKalRM55QkRKZTzyErHH8Gs8ZPXZ24Mo3k2YAAADANkVK+OmtWf3V5tEG8ixY4PZEGndPpvHfZZ4FC6jNow301qz+hC0gB/LVNVwAAACwjZePp/q90UGPvtBavy3fr3WLdyvywnUlxifJ08tD/mWLqk2v+noopDbXbAEOIHABAAD8g3n7FVTbxxqq7WMNGfIGGIAhhQAAAABgEAIXAAAAABiEwAUAAAAABiFwAQAAAIBBCFwAAAAAYBACFwAAAAAYhMAFAAAAAAYhcAEAAACAQQhcAAAAAGAQAhcAAAAAGITABQAAAAAGIXABAAAAgEEIXAAAAABgEAIXAAAAABiEwAUAAAAABiFwAQAAAIBBCFwAAAAAYBACFwAAAAAYhMAFAAAAAAYhcAEAAACAQQhcAAAAAGCQArldAAAAeY3J5PhzHdkGMkd/jUV/jUV/jeeMHjv752Mym81m524SAAAAACAxpBAAAKeKjo5WgwYNFB0dndul5Ev011j011j013h5sccELgAAnCwmJia3S8jX6K+x6K+x6K/x8lqPCVwAAAAAYBACFwAAAAAYhMAFAIATeXp6avjw4fL09MztUvIl+mss+mss+mu8vNhjZikEAAAAAINwhgsAAAAADELgAgAAAACDELgAAAAAwCAELgAAAAAwCIELAAAAAAxC4AIAAC7j9OnT2rdvn+Li4nK7lHzpxo0bmjBhgo4dO5bbpeQ7ZrNZMTExuV1GvnHr1i3Fx8enW3bjxg1t3bpVmzdv1pUrV3KpsowIXAAA2On48eN6/fXX9cgjj2jw4MFasmSJrN1lJSwsTDVq1MiFCl3bjz/+qJCQELVo0UKvv/66bt26patXr+qxxx5Tp06d9Nhjj+mhhx7SnDlzcrvUfCc6Olpz587VuXPncrsUl/Tnn3/q+vXr6ZYdPnxYzz77rOrWratGjRqpXr16Gj58uE6ePJk7Rbq4+Ph4vfDCC2rSpIkaNmyojz76SJI0b948tW7dWs8884yeffZZtW7dWmPHjlVqamouV8x9uAAAsMvp06fVs2dPpaSkqGrVqrp69aouX76sBg0a6NNPP1XJkiUt64aFhemNN97QoUOHcrFi17Ju3To9//zzCgwMVOnSpbVp0yY9/PDDSk1NVUxMjDp06KD4+HiFhobq0KFD+vzzz9W6devcLttldO3aNcvHk5OTderUKZUpU0a+vr4ymUwKCwu7R9W5vho1amjSpEmWPu/fv199+/aVJLVt21alS5fWuXPntHbtWvn6+uqnn35S+fLlc7NklzNjxgxNmzZNPXr0kL+/vxYsWKB+/frpiy++UI8ePdS2bVslJSUpIiJCa9eu1csvv6znnnsuV2sukKt7BwDAxfznP/+Rr6+v5s2bp4oVK0qSli5dqvfff1+PP/64vv76a1WuXDmXq3Rd33zzjRo3bqy5c+fKZDJp9uzZmjRpklq2bKlvvvnGst5TTz2lbt266bvvviNw2eHYsWPy8fFRrVq1rD6emJgoSfL19VXRokXvYWX5w93nMSZNmiRvb28tXLhQFSpUsCw/fPiwnnzySc2YMUMTJ06812W6tIiICPXq1UsffvihJKlmzZp69dVX9cgjj2j8+PGW9Tp16qRnn31WS5YsyfXAxZBCAADssHfvXvXt29cStiSpe/fuWrhwodzc3PTUU0/pzz//zMUKXdvx48fVvn17mUwmSbfPCqSmpqpTp07p1itQoIC6du2q/fv350aZLmvkyJEym80qUKCAxo4dq++++y7d17///W9J0ksvvWRZhpxJSUnRzp079cwzz6QLW5IUGBio3r1767fffsul6lzXpUuXVL9+fcv3DRo0kNlstvrBS5s2bXT+/Pl7WJ11BC4AAOxw/fp1+fv7Z1hepUoVLViwQKVLl9aAAQO0adOmXKjO9cXHx8vLy8vyvZ+fnyTpvvvuy7Cuv78/kxDY6fnnn9fKlStVtGhR9ezZU++99166a47Sgi4cFx8fbxl6bE21atUUFRV1j6tyfUWKFEl3zKb9/93XzqUtywtnaglcAADYoWzZsjpy5IjVx/z9/fX999+rZs2alje2sI+/v7/+/vtvy/deXl56/PHHVbp06QzrXr58OU+8mXI1pUqV0pQpU/Ttt99q165dat++vWbPnq3k5OTcLi1f2L9/v1atWqUtW7bI19dX165ds7re1atXLR8owHYNGjTQggULdOLECV2/fl1Tp06Vl5eXfvnll3R/O86cOaN58+apZs2auVjtbVzDBQCAHZo0aaKVK1fqjTfeUIECGf8Z9fPz07fffquRI0dq7dq1nDGwU61atbRnzx7L997e3nrvvfesrrtjxw4FBATco8ryn8aNG2vJkiWaP3++pk6dqvnz56t///4csw6aM2dOuhk0N2zYoF69emVYb/fu3emGJsM2L730kh599FF16dJF0u3r5kaOHKmyZcuqffv2qlmzplJTU3Xw4EGlpqZq+PDhuVwxgQsAALv07NlTkZGR2r9/v+rVq2d1HU9PT82YMUMTJkzQ4cOH722BLu7FF1/UhQsXsl0vKipKfn5+CgkJuQdV5V9ubm7q06ePOnfurClTpuiDDz7I7ZJc2ty5czMs8/DwyLAsKipK8fHx2c4aiYwqVqyo8PBwhYaGKiYmRk2aNNG//vUvSZK7u7vmz5+vqKgoPfTQQxo8eLDq1KmTyxUzLTwAAAD+68yZM7p8+bKqVaumYsWK5XY5QL5A4AIAAAAAgzBpBgAAAAAYhMAFAAAAAAYhcAEAAACAQQhcAAAAAGAQAhcAADkQFxen4cOHKywsLLdLyZfor7Hor7Hor/FcqccELgAAcsDb21u//fab4uPjc7uUfIn+Gov+Gov+Gs+VekzgAgAghxo2bKjdu3fndhn5Fv01Fv01Fv01nqv0mMAFAEAOjR07Vjt37tSUKVP0119/5XY5+Q79NRb9NRb9NZ6r9JgbHwMAkEP169dXSkqKkpKSJEnu7u7y9PRMt47JZNLOnTtzozyXR3+NRX+NRX+N5yo9LpCrewcAwIV16NBBJpMpt8vIt+ivseivseiv8Vylx5zhAgAAAACDcA0XAAAAABiEIYUAADggOjpas2fP1vr163Xx4kVJUpkyZdS6dWsNHDhQfn5+uVyha6O/xqK/xqK/xnOFHjOkEACAHLp8+bL69Omj8+fPq3LlyqpcubIk6dSpUzpx4oTKly+vefPm6b777svlSl0T/TUW/TUW/TWey/TYDAAAcmTUqFHmunXrmtevX5/hsfXr15vr1atnfv3113OhsvyB/hqL/hqL/hrPVXrMNVwAAOTQpk2bNGDAALVq1SrDY61atVK/fv20YcOGXKgsf6C/xqK/xqK/xnOVHhO4AADIobi4OJUoUSLTx/39/RUXF3cPK8pf6K+x6K+x6K/xXKXHBC4AAHKoSpUqWrZsmRITEzM8lpSUpGXLlqlKlSq5UFn+QH+NRX+NRX+N5yo9ZtIMAAByaMWKFXr55ZcVEBCgp556Sg888ICk2xdsL1iwQEeOHNGUKVPUsWPH3C3URdFfY9FfY9Ff47lKjwlcAAA4YPHixfr3v/+tq1evymQySZLMZrNKlCihUaNGqWfPnrlcoWujv8aiv8aiv8ZzhR4TuAAAcFBycrL279+f7h4wtWvXVoEC3O7SGeivseivseiv8fJ6jwlcAADkUGhoqBo1aqRy5cpZffz8+fPasWOHevTocW8Lyyfor7Hor7Hor/FcpcdMmgEAQA6NGTNGu3fvzvTxP//8U2PGjLmHFeUv9NdY9NdY9Nd4rtJjAhcAADmU3SCR2NhYubu736Nq8h/6ayz6ayz6azxX6XHeGNgIAICLOHz4sA4fPmz5fseOHUpJScmw3s2bN7VgwQJVqlTpXpbn8uivseivseiv8Vyxx1zDBQCAHaZPn67p06dLkkwmU5afsBYuXFgfffSR2rRpc6/Kc3n011j011j013iu2GMCFwAAdvj777/1999/y2w2q3fv3hoxYoRatmyZbh2TySRvb29VqFAhz8yS5Sror7Hor7Hor/FcsccELgAAcuj3339X1apVVbx48dwuJV+iv8aiv8aiv8ZzlR4zaQYAADlUvXp1/f3335k+fuTIEd24ceMeVpS/0F9j0V9j0V/juUqPOcMFAEAOvfHGGzp16pQWLVpk9fEnnnhClStX1ocffniPK8sf6K+x6K+x6K/xXKXHnOECACCHtm3bpuDg4Ewfb9OmjbZu3XoPK8pf6K+x6K+x6K/xXKXHBC4AAHIoKipKxYoVy/TxokWL6urVq/ewovyF/hqL/hqL/hrPVXpM4AIAIIdKliypgwcPZvr4gQMH8vzF3HkZ/TUW/TUW/TWeq/SYwAUAQA61a9dOP//8s9asWZPhsV9//VWLFy9Wu3btcqGy/IH+Gov+Gov+Gs9VesykGQAA5NCtW7f01FNP6fjx4woMDFS1atUkSceOHdPhw4dVpUoV/fDDDypcuHAuV+qa6K+x6K+x6K/xXKXHBC4AABwQGxurr7/+WqtXr9bZs2clSRUqVFD79u01aNAg+fj45HKFro3+Gov+Gov+Gs8VekzgAgAAAACDcA0XAAAAABikQG4XAACAK0tISNAvv/yigwcP6tatW0pNTU33uMlkyvWbbroy+mss+mss+ms8V+gxgQsAgBy6cOGC+vfvrwsXLqhw4cK6deuWihQpolu3biklJUXFihXLE9cPuCr6ayz6ayz6azxX6TFDCgEAyKFJkyYpOjpaixYt0sqVK2U2mzVlyhTt3r1bo0aNkpeXl7755pvcLtNl0V9j0V9j0V/juUqPCVwAAOTQtm3b9OSTTyooKEhubv/7J9XT01ODBw9W06ZNc30oiyujv8aiv8aiv8ZzlR4TuAAAyKH4+HiVLVtWkuTn5yeTyaRbt25ZHq9fv7527tyZW+W5PPprLPprLPprPFfpMYELAIAcuv/++3X58mVJUoECBVSqVCnt2bPH8vjx48dVsGDBXKrO9dFfY9FfY9Ff47lKj5k0AwCAHGratKnWrFmj4cOHS5J69uypmTNn6ubNm0pNTVVYWJi6d++ey1W6LvprLPprLPprPFfpMTc+BgAghy5evKh9+/apTZs28vT0VEJCgsaNG6dVq1bJzc1Nbdq00dtvvy0/P7/cLtUl0V9j0V9j0V/juUqPCVwAANhh8uTJCgkJUWBgYG6Xki/RX2PRX2PRX+O5Yo+5hgsAADvMnDlTx44ds3x/7do11ahRQ1u3bs3FqvIP+mss+mss+ms8V+wxgQsAAAcxWMRY9NdY9NdY9Nd4eb3HBC4AAAAAMAiBCwAAAAAMwrTwAADY6cKFCzpw4IAkWW6yeebMGRUuXNjq+rVq1bpnteUH9NdY9NdY9Nd4rtZjZikEAMAOgYGBMplM6ZaZzeYMy+5cfujQoXtVnsujv8aiv8aiv8ZzxR5zhgsAADtMmDAht0vI1+ivseivseiv8Vyxx5zhAgAAAACDMGkGAAAAABiEwAUAAAAABiFwAQAAAIBBCFwAAAAAYBACFwAAAAAYhMAFAAAAAAYhcAEAAACAQQhcAAAAAGCQ/wcfAcMizp2DxwAAAABJRU5ErkJggg==",
"text/plain": [
"
"
],
"text/plain": [
" Factor Term NES P-value Adj. P-value\n",
"41 Factor 4 APOPTOSIS -1.260816 0.001 0.016266"
]
},
"execution_count": 13,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"gsea_df.loc[(gsea_df['Adj. P-value'] < 0.05) & (gsea_df['NES'] < 0.)]"
]
},
{
"cell_type": "markdown",
"id": "56fd5694",
"metadata": {},
"source": [
"## Footprint Enrichment"
]
},
{
"cell_type": "markdown",
"id": "ad5bb159",
"metadata": {},
"source": [
"Footprint enrichment analysis build upon classic geneset enrichment analysis, as instead of considering the genes involved in a biological activity, they consider the genes affected by the activity, or in other words the genes that change downstream of said activity [(Dugourd and Saez-Rodriguez, 2019)](https://www.sciencedirect.com/science/article/pii/S2452310019300149). This is particularly useful when working with transcriptomics data, as it allows to identify the genes that are regulated by a biological process, instead of the genes that are part of the process. In this case, we will use the PROGENy pathway resource to perform footprint enrichment analysis. PROGENy was built in a data-driven manner using perturbation and cancer lineage data [(Schubert et al, 2019)](https://www.nature.com/articles/s41467-017-02391-6#Sec8), as a consequence it also assigns different importances or weights to each gene in its pathway genesets.\n",
"\n",
"To this end, we need an enrichment method that can take weights into account, and here we will use multi-variate linear models from the `decoupler-py` package to perform this analysis [(Badia-i-Mompel et al., 2022)](https://academic.oup.com/bioinformaticsadvances/article/2/1/vbac016/6544613)."
]
},
{
"cell_type": "markdown",
"id": "ca5209d4",
"metadata": {},
"source": [
"### load the PROGENy pathway resource:"
]
},
{
"cell_type": "code",
"execution_count": 14,
"id": "68972dfc",
"metadata": {},
"outputs": [],
"source": [
"net = dc.get_progeny(organism='human', top=5000)"
]
},
{
"cell_type": "markdown",
"id": "bac3fcc0",
"metadata": {},
"source": [
"Now we will use the PROGENy genesets to assign pathways and weights to the ligand-receptor pairs in liana's consensus resource. Specifically, we will use the weighted bipartite networks from PROGENy, where the weight represents the importance of the genes to a given geneset, and assign a weight to each ligand-receptor interaction, based on the mean of the weights of the proteins in that ligand-receptor interaction. We keep only ligand-receptor weights only if all the protein in the ligand-receptor interaction are present for a given pathway, and are sign-coherent."
]
},
{
"cell_type": "code",
"execution_count": 15,
"id": "25731683",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"