217 lines
101 KiB
Plaintext
217 lines
101 KiB
Plaintext
{
|
|
"cells": [
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 1,
|
|
"id": "d859e6a3-9237-4c6c-adad-8a88917ec849",
|
|
"metadata": {},
|
|
"outputs": [],
|
|
"source": [
|
|
"import json\n",
|
|
"import numpy as np\n",
|
|
"import matplotlib.pyplot as plt\n",
|
|
"\n",
|
|
"from aligner import *"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": null,
|
|
"id": "09cd94af-994f-4c36-a88a-530000ecce40",
|
|
"metadata": {},
|
|
"outputs": [],
|
|
"source": []
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 2,
|
|
"id": "573034ed-a259-4ce9-96aa-5f07a440ada4",
|
|
"metadata": {},
|
|
"outputs": [
|
|
{
|
|
"name": "stdout",
|
|
"output_type": "stream",
|
|
"text": [
|
|
"Titan's Call\n",
|
|
"280 280 280\n",
|
|
"a) total 282 selected 280\n",
|
|
"b) total 281 selected 280\n"
|
|
]
|
|
}
|
|
],
|
|
"source": [
|
|
"j = 0 # track id (0-based)\n",
|
|
"i_a, i_b = 0, 1 # annotation run ids (0-based)\n",
|
|
"\n",
|
|
"annot_a = load_track_annot(fns_all[i_a][j])\n",
|
|
"annot_b = load_track_annot(fns_all[i_b][j])\n",
|
|
"assert annot_a['contentId'] == annot_b['contentId']\n",
|
|
"\n",
|
|
"print(annot_a['title'])\n",
|
|
"\n",
|
|
"ts_a, ts_b, tsa_a, tsa_b = align(annot_a, annot_b)\n",
|
|
"\n",
|
|
"print('a) total', len(ts_a), 'selected', len(tsa_a))\n",
|
|
"print('b) total', len(ts_b), 'selected', len(tsa_b))"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 3,
|
|
"id": "c650bfa4-5c38-4b0a-81f0-0544ae711ec9",
|
|
"metadata": {},
|
|
"outputs": [
|
|
{
|
|
"data": {
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAqsAAADcCAYAAABXo+ihAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABgxUlEQVR4nO2deXQUVfr3v52QHZIQEpKwSEIISwyCRgKIAYUgEQR3AZcRxgM/Fd5hUUdxRiFuiDKIw66j4KAg4r4gGlYRIkEWEQIRmAAqCRFCFghJIH3fP2IV1d213Oqu6q4Oz+ccziHV1VW3696697nPamOMMRAEQRAEQRCEBQnwdQMIgiAIgiAIQgkSVgmCIAiCIAjLQsIqQRAEQRAEYVlIWCUIgiAIgiAsCwmrBEEQBEEQhGUhYZUgCIIgCIKwLCSsEgRBEARBEJaFhFWCIAiCIAjCspCwShAEQRAEQVgWElYJAsCMGTNgs9lw6tQpXzeF8COSkpIwZswYXzdDF/7YZnc4evQobDYbli1b5tb3ly1bBpvNhqNHjxraLiljxoxBUlKSodf0RrsJwtuQsEoQJrFixQrMnTvX1824rKmpqcGMGTOwadMmt6+xbds2zJgxAxUVFYa1y19ZuHCh28IfYSwvvfQSPv30U183gyC8AgmrBGESJKz6npqaGuTm5nosrObm5soKq0VFRXjzzTfdb6AP8KTNJKway5tvvomioiK3vqskrD7wwAM4f/48OnTo4GHrCMI6NPN1AwiCIPyVkJAQXzdBN1Zr88WLF2G32xEcHOzrpnidoKAgw68ZGBiIwMBAw69LEL6ENKsEIeHUqVO45557EBkZiVatWmHSpEmora11Oe/dd99FRkYGwsLCEBMTg1GjRuHXX38VP7/hhhvw1Vdf4dixY7DZbLDZbEhKSgJjDLGxsZg6dap4rt1uR3R0NAIDAx20d7NmzUKzZs1w9uxZ8djBgwdx1113ISYmBqGhobj22mvx+eefu7SvoqICkydPRvv27RESEoJOnTph1qxZsNvt4jmCT9/s2bPxxhtvICUlBSEhIejVqxd27Nih+azKy8vx+OOPo3v37mjevDkiIyNx880346effnI4b9OmTbDZbPjggw/w4osvol27dggNDcWgQYNw+PBhh3NvuOEGpKeno7CwEDfeeCPCw8PRtm1bvPLKKy73Lysrw0MPPYT4+HiEhoaiR48eeOeddxx+X1xcHAAgNzdX7IcZM2YAAPbu3YsxY8agY8eOCA0NRUJCAv7617/i9OnT4jVmzJiBJ554AgCQnJwsXkPwB5Tz//zf//6Hu+++GzExMQgPD0efPn3w1Vdfuf1MDh06hDvvvBMJCQkIDQ1Fu3btMGrUKFRWVornnDp1CgcPHkRNTY1Sd4k4t1nwcdy6dSumTp2KuLg4RERE4Pbbb8cff/zh8L39+/dj8+bN4nO44YYbxM/1jrm5c+eKY2737t1o1qwZcnNzXdpbVFQEm82G+fPnA+Afd3rYv38/Bg4ciLCwMLRr1w4vvPCCQ7ulfP3118jKykJERARatGiBYcOGYf/+/eLns2fPhs1mw7Fjx1y+O23aNAQHB+PMmTMA5H1WZ8+ejeuuuw6tWrVCWFgYMjIy8OGHHzqcY7PZcO7cObzzzjtiXwh9quSzunDhQlx55ZUICQlBmzZtMGHCBBdrgZ73b968ebjyyisRHh6Oli1b4tprr8WKFStknxlBeAwjCIJNnz6dAWDdu3dnw4cPZ/Pnz2f3338/A8AeeOABh3NfeOEFZrPZ2MiRI9nChQtZbm4ui42NZUlJSezMmTOMMca+/fZb1rNnTxYbG8uWL1/Oli9fzj755BPGGGMjRoxgGRkZ4vV2797NALCAgAD25ZdfiseHDRvGrr32WvHvffv2saioKJaWlsZmzZrF5s+fz/r3789sNhv7+OOPxfPOnTvHrrrqKtaqVSv29NNPs8WLF7O//OUvzGazsUmTJonnFRcXMwDs6quvZp06dWKzZs1ir7zyCouNjWXt2rVj9fX1qs9sx44dLCUlhT311FNsyZIl7LnnnmNt27ZlUVFR7PfffxfP27hxo3ifjIwM9tprr7EZM2aw8PBwlpmZ6XDNAQMGsDZt2rD27duzSZMmsYULF7KBAwcyAGzNmjXieTU1Naxbt24sKCiITZkyhf373/9mWVlZDACbO3cuY4yxs2fPskWLFjEA7Pbbbxf74aeffmKMMTZ79myWlZXFnnvuOfbGG2+wSZMmsbCwMJaZmcnsdjtjjLGffvqJjR49mgFgr732mniNs2fPMsYY69ChA3vwwQfFdpWWlrL4+HjWokUL9o9//IPNmTOH9ejRgwUEBDj0Ee8zqaurY8nJyaxNmzbshRdeYP/5z39Ybm4u69WrFzt69Kh4njB+N27cqNpncm1eunSp2JaBAweyefPmsccee4wFBgaye+65Rzzvk08+Ye3atWNdu3YVn8O3337LGNM/5tLS0ljHjh3Zyy+/zF577TV27NgxNnDgQJaWlubS3tzcXBYYGMhKS0sZY/zjTrjX0qVLVZ9HSUkJi4uLYy1btmQzZsxgr776KktNTWVXXXUVA8CKi4vFc//73/8ym83GcnJy2Lx589isWbNYUlISi46OFs87duwYs9ls7JVXXnG5V8eOHdmwYcPEvx988EHWoUMHh3PatWvHHn30UTZ//nw2Z84clpmZyQA4zA3Lly9nISEhLCsrS+yLbdu2OfSntN3C+MjOzmbz5s1jEydOZIGBgaxXr14O7znv+/fGG28wAOyuu+5iS5YsYa+//jp76KGH2N/+9jfVZ00Q7kLCKkGwS5P5iBEjHI4/+uijDIAo4Bw9epQFBgayF1980eG8n3/+mTVr1szh+LBhw1wWIsYYe/XVV1lgYCCrqqpijDH273//m3Xo0IFlZmayJ598kjHGWENDA4uOjmZTpkwRvzdo0CDWvXt3VltbKx6z2+3suuuuY6mpqeKx559/nkVERLBffvnF4b5PPfUUCwwMZMePH2eMXVrMW7VqxcrLy8XzPvvsMwaAffHFF6rPrLa2ljU0NDgcKy4uZiEhIey5554TjwmCWbdu3VhdXZ14/PXXX2cA2M8//yweGzBgAAPA/vvf/4rH6urqWEJCArvzzjvFY3PnzmUA2Lvvviseq6+vZ3379mXNmzcXn+0ff/zBALDp06e7tL+mpsbl2MqVKxkA9t1334nHXn31VZfFX8BZ8Js8eTIDwLZs2SIeq66uZsnJySwpKUl8XrzPRNjIrF692uXeUowQVrOzs0UhnTHGpkyZwgIDA1lFRYV47Morr2QDBgxwuabeMRcZGcnKysoczl2yZInLeGCMsbS0NDZw4EDxb95xxyusCn22fft28VhZWRmLiopy6Pfq6moWHR3Nxo0b5/D90tJSFhUV5XC8b9++DhtSxhgrKChwGdtywqrzuKyvr2fp6ekOz4AxxiIiIhz6UcBZWC0rK2PBwcHspptucnhu8+fPZwDY22+/LR7jff9uvfVWduWVV7rcmyDMgtwACELChAkTHP7+f//v/wEA1qxZAwD4+OOPYbfbcc899+DUqVPiv4SEBKSmpmLjxo2a98jKykJDQwO2bdsGANiyZQuysrKQlZWFLVu2AAD27duHiooKZGVlAWg0fW7YsAH33HMPqqurxfuePn0aQ4YMwaFDh/D7778DAFavXo2srCy0bNnSoY3Z2dloaGjAd99959CekSNHomXLlg7tAxrN2WqEhIQgIKBxCmloaMDp06fRvHlzdOnSBbt27XI5f+zYsQ5+iUr3ad68Oe6//37x7+DgYGRmZjqct2bNGiQkJGD06NHisaCgIPztb3/D2bNnsXnzZtW2A0BYWJj4/9raWpw6dQp9+vQBANn287BmzRpkZmbi+uuvd/g948ePx9GjR1FYWOhwvtYziYqKAgB88803qib+GTNmgDHmYJbXy/jx42Gz2Rza0tDQIGvOdkbvmLvzzjtFFw2BO+64A82aNcOqVavEY/v27UNhYSFGjhwpHtM77rRYs2YN+vTpg8zMTPFYXFwc7rvvPofz8vLyUFFRgdGjRzv8xsDAQPTu3dvh3R85ciR27tyJI0eOiMdWrVqFkJAQ3HrrrartkY7LM2fOoLKyEllZWW6PyXXr1qG+vh6TJ08WnxsAjBs3DpGRkS4uKjzvX3R0NH777TcudyGCMAISVglCQmpqqsPfKSkpCAgIEP2/Dh06BMYYUlNTERcX5/DvwIEDKCsr07zHNddcg/DwcFEwFYTV/v3748cff0Rtba34mSD0HD58GIwxPPPMMy73nT59OgCI9z506BDWrl3rcl52drbDeQJXXHGFw9+C4Cr41Slht9vx2muvITU1FSEhIYiNjUVcXBz27t3r4E+p9z7t2rVzEJqEc6XnHTt2DKmpqQ6LLwB069ZN/FyL8vJyTJo0CfHx8QgLC0NcXBySk5MBQLb9PBw7dgxdunRxOa7ULq1nkpycjKlTp+I///kPYmNjMWTIECxYsMDt9qnh7jgA9I854TlLiY2NxaBBg/DBBx+Ix1atWoVmzZrhjjvuEI/pHXdaCGPJGed+PHToEABg4MCBLr/z22+/dfiNd999NwICAkTBmzGG1atX4+abb0ZkZKRqe7788kv06dMHoaGhiImJQVxcHBYtWuTRmJT7PcHBwejYsaPLmOR5/5588kk0b94cmZmZSE1NxYQJE7B161a32kcQPFA2AIJQwXnSttvtsNls+Prrr2Ujbps3b655zaCgIPTu3RvfffcdDh8+jNLSUmRlZSE+Ph4XLlzA9u3bsWXLFnTt2lXUPgnBHo8//jiGDBkie91OnTqJ5w4ePBh///vfZc/r3Lmzw99KkcOMMdXf8dJLL+GZZ57BX//6Vzz//POIiYlBQEAAJk+eLBucwnsfd9ujl3vuuQfbtm3DE088gZ49e6J58+aw2+3IyclRDK4xGp7f+q9//QtjxozBZ599hm+//RZ/+9vfMHPmTPzwww9o166dV9uihN4xJ9UeShk1ahTGjh2LPXv2oGfPnvjggw8waNAgxMbGiufoHXdGIVx7+fLlSEhIcPm8WbNLy2mbNm2QlZWFDz74AE8//TR++OEHHD9+HLNmzVK9x5YtWzBixAj0798fCxcuRGJiIoKCgrB06VKvBS/xjINu3bqhqKgIX375JdauXYuPPvoICxcuxLPPPisbJEcQnkLCKkFIOHTokIPW5/Dhw7Db7WLEbkpKChhjSE5OdlmAnXEWdKVkZWVh1qxZWLduHWJjY9G1a1fYbDZceeWV2LJlC7Zs2YJbbrlFPL9jx44AGgVdQVulREpKCs6ePat5nqd8+OGHuPHGG/HWW285HK+oqHAQLsygQ4cO2Lt3L+x2u4N29eDBg+LngHIfnDlzBuvXr0dubi6effZZ8bigPZOi1o9y7ZLLm+ncLr10794d3bt3xz//+U9s27YN/fr1w+LFi/HCCy+4dT13UXoWRo252267Df/3f/8naiR/+eUXTJs2zeEco8ddhw4dZPvduR9TUlIAAK1bt+b6nSNHjsSjjz6KoqIirFq1CuHh4Rg+fLjqdz766COEhobim2++cUgxtnTpUpdzecelMOaKiorEeQQA6uvrUVxc7HafRUREYOTIkRg5ciTq6+txxx134MUXX8S0adMQGhrq1jUJQglyAyAICQsWLHD4e968eQCAm2++GUCjX11gYCByc3NdNE6MMYe0RxEREYqmu6ysLNTV1WHu3Lm4/vrrxYUnKysLy5cvx4kTJ0T/RaBxgbzhhhuwZMkSlJSUuFxPmmLonnvuQX5+Pr755huX8yoqKnDx4kXVZ8BLYGCgyzNYvXq16DtrJkOHDkVpaamDf+PFixcxb948NG/eHAMGDAAAhIeHA4BLih5Be+TcfrkiDhEREbLXUGpXQUEB8vPzxWPnzp3DG2+8gaSkJKSlpWleQ0pVVZVLf3Xv3h0BAQGoq6sTj+lJXeUJERERss/BqDEXHR2NIUOG4IMPPsD777+P4OBg3HbbbQ7nGD3uhg4dih9++AEFBQXisT/++APvvfeew3lDhgxBZGQkXnrpJVy4cMHlOtJ3EGj0yw0MDMTKlSuxevVq3HLLLeJYUiIwMBA2mw0NDQ3isaNHj8om/1fqC2eys7MRHByMf//73w7P7a233kJlZSWGDRumeQ1npPMc0OhSkJaWBsaY7LMhCE8hzSpBSCguLsaIESOQk5OD/Px8vPvuu7j33nvRo0cPAI3alRdeeAHTpk3D0aNHcdttt6FFixYoLi7GJ598gvHjx+Pxxx8HAGRkZGDVqlWYOnUqevXqhebNm4ualb59+6JZs2YoKirC+PHjxfv3798fixYtAgAHYRVoFKSvv/56dO/eHePGjUPHjh1x8uRJ5Ofn47fffhPzTD7xxBP4/PPPccstt2DMmDHIyMjAuXPn8PPPP+PDDz/E0aNHDdF83nLLLXjuuecwduxYXHfddfj555/x3nvvOWhvzGL8+PFYsmQJxowZg507dyIpKQkffvghtm7dirlz56JFixYAGs3NaWlpWLVqFTp37oyYmBikp6cjPT0d/fv3xyuvvIILFy6gbdu2+Pbbb1FcXOxyr4yMDADAP/7xD4waNQpBQUEYPny4rODx1FNPYeXKlbj55pvxt7/9DTExMXjnnXdQXFyMjz76yMXHVosNGzZg4sSJuPvuu9G5c2dcvHgRy5cvR2BgIO68807xvPnz5yM3NxcbN270KMhKi4yMDCxatAgvvPACOnXqhNatW2PgwIGGjrmRI0fi/vvvx8KFCzFkyBBER0c7fG70uPv73/+O5cuXIycnB5MmTUJERATeeOMNUXsvEBkZiUWLFuGBBx7ANddcg1GjRiEuLg7Hjx/HV199hX79+om5YIHGDeaNN96IOXPmoLq62iFITIlhw4Zhzpw5yMnJwb333ouysjIsWLAAnTp1cmgL0NgX69atw5w5c9CmTRskJyejd+/eLteMi4vDtGnTkJubi5ycHIwYMQJFRUVYuHAhevXq5RBMxctNN92EhIQE9OvXD/Hx8Thw4ADmz5+PYcOGie8eQRiKt9MPEIQVEVL/FBYWsrvuuou1aNGCtWzZkk2cOJGdP3/e5fyPPvqIXX/99SwiIoJFRESwrl27sgkTJrCioiLxnLNnz7J7772XRUdHMwAuKWp69erlkjLnt99+YwBY+/btZdt55MgR9pe//IUlJCSwoKAg1rZtW3bLLbewDz/80OG86upqNm3aNNapUycWHBzMYmNj2XXXXcdmz54t5lUUUvu8+uqrLveBQronKbW1teyxxx5jiYmJLCwsjPXr14/l5+ezAQMGOKQ3EtI0OadfkkstNGDAANmUOHIpfk6ePMnGjh3LYmNjWXBwMOvevbtsmqJt27axjIwMFhwc7PC7fvvtN3b77bez6OhoFhUVxe6++2524sQJ2d/+/PPPs7Zt27KAgACHtEDOaaAYa+yju+66i0VHR7PQ0FCWmZnpkCNTzzP53//+x/7617+ylJQUFhoaymJiYtiNN97I1q1b5/A9I1JX7dixQ7aN0muWlpayYcOGsRYtWjAADv3s6ZgTqKqqYmFhYS6pyQR4xx1v6irGGNu7dy8bMGAACw0NZW3btmXPP/88e+utt2RTlm3cuJENGTKERUVFsdDQUJaSksLGjBnDfvzxR5frvvnmmwwAa9Gihew8Ijeu33rrLZaamspCQkJY165d2dKlS8X+lXLw4EHWv39/8VkJfSqXZ5WxxlRVXbt2ZUFBQSw+Pp498sgjYl5oAd73b8mSJax///6sVatWLCQkhKWkpLAnnniCVVZWunyXIIzAxpjBUQsEQRAEQRAEYRDks0oQBEEQBEFYFhJWCYIgCIIgCMtCwipBEARBEARhWUhYJQiCIAiCICwLCasEQRAEQRCEZSFhlSAIgiAIgrAsTa4ogN1ux4kTJ9CiRQtdZRIJgiAIgiAI78AYQ3V1Ndq0aaNZMKXJCasnTpxA+/btfd0MgiAIgiAIQoNff/0V7dq1Uz2nyQmrQqm3X3/9FZGRkT5uDUEQBEEQBOFMVVUV2rdvz1Wit8kJq4LpPzIy0i+F1QY7Q0FxOcqqa9G6RSgyk2MQGEDuDARBEARBND14XDabnLDqz6zdV4LcLwpRUlkrHkuMCsX04WnISU/0YcsIgiAIgiB8A2UDsAhr95XgkXd3OQiqAFBaWYtH3t2FtftKfNQygiAIgiAI30HCqgVosDPkflEIJvOZcCz3i0I02OXOIAiCIAiCaLqQsGoBCorLXTSqUhiAkspaFBSXe69RBEEQBEEQFoCEVQtQVq0sqLpzHkEQBEEQRFOBhFUL0LpFqKHnEQRBEARBNBVIWLUAmckxSIwKhVLyBhsaswJkJsd4s1kEQRAEQRA+h4RVCxAYYMP04WkA4CKwCn9PH55G+VYJgiAIgrjsIGHVIuSkJ2LR/degdWSIw/GEqFAsuv8ayrNKEARBEMRlCRUFsBA56Yno1ykW3Wd8CwBYNrYXslLjSKNKEARBEMRlC2lWLYZUMKVSqwRBEARBXO6QsEoQBEEQBEFYFhJWCYIgCIIgCMtCwipBEARBEARhWUhYJQiCIAiCICwLCasEQRAEQRCEZSFhlSAIgiAIgrAspgurCxYsQFJSEkJDQ9G7d28UFBQonrt//37ceeedSEpKgs1mw9y5c81uHkEQBEEQBGFhTBVWV61ahalTp2L69OnYtWsXevTogSFDhqCsrEz2/JqaGnTs2BEvv/wyEhISzGwaQRAEQRAE4QeYKqzOmTMH48aNw9ixY5GWlobFixcjPDwcb7/9tuz5vXr1wquvvopRo0YhJCRE9hyCIAiCIAji8sE0YbW+vh47d+5Ednb2pZsFBCA7Oxv5+fmG3aeurg5VVVUO/wiCIAiCIIimgWnC6qlTp9DQ0ID4+HiH4/Hx8SgtLTXsPjNnzkRUVJT4r3379oZdmyAIgiAIgvAtfp8NYNq0aaisrBT//frrr75uEkEQBEEQBGEQzcy6cGxsLAIDA3Hy5EmH4ydPnjQ0eCokJIT8WwmCIAiCIJoopmlWg4ODkZGRgfXr14vH7HY71q9fj759+5p1W4IgCIIgCKIJYZpmFQCmTp2KBx98ENdeey0yMzMxd+5cnDt3DmPHjgUA/OUvf0Hbtm0xc+ZMAI1BWYWFheL/f//9d+zZswfNmzdHp06dzGwqQRAEQRAEYUFMFVZHjhyJP/74A88++yxKS0vRs2dPrF27Vgy6On78OAICLil3T5w4gauvvlr8e/bs2Zg9ezYGDBiATZs2mdlUgiAIgiAIwoLYGGPM140wkqqqKkRFRaGyshKRkZG+bo5uauovIu3ZbwAAhc8NQXiwqfsJgiAIgiAIr6NHXiNJ6DKhwc5QUFyOsupatG4RiszkGAQG2HzdLIIgCIIgCFVIWL0MWLuvBLlfFKKkslY8lhgViunD05CTnujDlhEEQRAEQajj93lWCXXW7ivBI+/uchBUAaC0shaPvLsLa/eV+KhlBEEQBKFNg50h/8hpfLbnd+QfOY0Ge5PyXiQ4IM1qE6bBzpD7RSHkXmsGwAYg94tCDE5LIJcAgiAIwnKQZZAASLPapCkoLnfRqEphAEoqa1FQXO69RhEEQRAEB2QZJARIs9qEKatWFlTdOY8g/AkKKiQI/8WblkGaK6wPCatNmNYtQg09z1vQxEHoQW685BWWkumQIPwYPZbBvimt3L4PuRn4BySsNmEyk2OQGBWK0spa2d2pDUBCVOPibhVo4iD0IDdeosODUFFzweVcwXS46P5raCwRhMXxhmVQcDNwXh9prrAe5LPqI7wR3RgYYMP04WkAGgVTKcLf04enWUZrSf5JTRczxrvSeJETVAGIC1LuF4UUTUwQFsdsy6CWmwHQtOcKf8uwQJpVH6CmPezfOc7Qe+WkJ2LR/ddg+uf7cbKqTjyeYDFtJWUuaLqYoS1XGy9qGGU69FfIxYbwF8y2DHrLzcCK+KMFkzSrXkZLe5hXeNLwe+akJ2Ld1AHi38vG9sL3Tw601KCkzAVNE7O05VrjRYvLMahw7b4SXD9rA0a/+QMmvb8Ho9/8AdfP2kAWC8KSmG0ZvFwDkP3VgknCqhfhMTvMXHPQlHtLX2gralMu14nDX3DHZGSmmc3TceCLoEJfmt3MXqD8zaRI+AeCZbB1ZIjD8YSoUI/9Sf01ANkT/Nn1gdwAvAiP9rC06vIUxqw4cZDJtBF3TUZmmtncHQe+Cio00uymd1ya7WLjTyZFX7zTNI94Rk56Ivp1ikX3Gd8CaLQMZqXGefwM/TEA2VP82fWBhFUvQlpBZaw2cfjTAmwmnkTLmqkt1xovcvgqqNDIiGN3xqWZC5Q/RVP74p2mecQYzLAMCm4Gj7y7CzbAYQxbMQDZCPzZgkluAF6kKZkTjMZKmQv81afHaDw1GZmpLecZL9FhQQ7HjTAd6sVIs5u749KsBcqfTIq+eKdpHrE+ZroZ8OJNFxorWjB5IWHViwjaICVxywYgIdJ6g8RbWGXi8JcF2Gw8DXrjGe+JHmjL1cbL4vuvwZYnbxSP+Sqo0KjAQU/GpVkLlJFBkWYu2L54p608j5B/sSO+DED2dtCj2XOymZCw6kV4tEHThnb1apushq8zF1BWgkt4qpHzhrZcbbxYIajQKK2mJ+PSrAXKqN9m9oLti3faqvMIZYSQxxdzhS8071ayYOqFhFUvo6U9HJwW76OWWQelicMbGgE9C3BT11AYoZHzhrbcCkKpEkZpNT0RDM1aoIz4bd5YsH3hp2dF30ByS7AOvtS8W8GC6Q4UYOUD1KIba+ov+rh11sRbgQq8C/DRUzW4ftaGJh04YVTQm1nRvP4ATyBYgA04c65e9TqeCoZmFAfxdHy4k6XAnch6X/jpWc03kIquWAtfR+X745xMmlUfYWVtkNXwpkaAx2QaHR6Euet+afIaCiM1cpfreJc+QyXsDJiwQn3cGGHKN9rFxtPxoddU7q4J2xd+elbzDbSqW8LlihU07/42J5OwSlgab5tLtBZg4S5WDJwwA381GVmJnPRELLj3amitBWrjxqiNg9ELlCfjQ8+CrXfDKnXRKSguxzPD1N/pUb3a48u9Jwxz57Gab6AVhCMemrprlYDVNO/+ALkBEJZm57EzXjeXqJlMR/Vqj9fWHfJqe3wNj8mIEp+r0zIiBGrrLs+4McOUbwTumhR5F+LYiBA8/uFP3CZsJZeh8f2T8emeEw7PLiq8Mb2Z9J02yp3HSv3lD8KRmqtX/85xPmsXL3rmQKvlFfcHSFglLM0f1XXaJ8F4jYDSAvzl3hM+aY+vUdPIUeJzbYzSbFnV18wdjS3vgg0buDeslefrFYsUvPFdMebc0wNTPvgJADDxxhQs2HiEu6CBOxsyq/SXt4QjdzetWsUl5o7q6VG7zEbvHHg5FiTwFBJWCUsT1yJE+ySYoxGQW4D9QUPhTfypgpEvMXLc+JuvmRK8C/aps3wb1tLK83jlmyJVDeystUXisdU7f/NYWysnjDgLbOltI8XPfNVf3hCO3N208gR/zVxzUFdbvGnpcXcOHJyWgMnZnfH21mJUnr8gHve1pcSqkLBKWJqMDi29Zi5RW2QEyHxzCT0Rxpc7mckxiI8McTAHS/H1uFFa3M1e9HlM5flHTnNdq/xcvaYGtrTq0udKfSGcy6OtdRZG5AS2+Ei+DbfZmOmW4MmmlSf4S9pvPG3xlqXH3SwLcm0EGrX9UwZ38dsNqJmQsEpYGjM1AtKF+OipGqwsOO4wKcotMmS+uYSeCOMe7aO81zALEhhgw9NDu2HS+3tcPvP1uFFa3Ef0SMTnP5WYvuhrmcp5N4gxzY0XCrccKsOK7b9yCSN5haWyAluZilBsBmobDDPcEjxNi2Wky5RRlh7eTZo7KaiU2ggACzYeQXrbKI/dT5oiJKwSlscMjYDSzlaK0iLjzcAJK09U+vwwL29hFYBiwQ9fmv2UFs6Sylos+a7Y5Xyz3DvUXBt4N4hRYcGGtUdg4ab/qX4uCCM/HDmtmbUEgOnR7TxaRaPdSDzNGWqUy5RRuWTzCk/i5a8Pcm3S9Pqiq7VRwF33k6YOpa4i/AIjc0QqpcFxRm2RcW7Pkzld8PecrogKCzZsQbJ6aUTy3/Ucb5cTlsKzcDrjq/RsPCmyeHKbJkReGovxkSGK5+ol/3+nNOcToDG7iVn4qkKVp8GDevtNCaNyyU5+fw/3M9Q7B+ppI1Ucc4SEVcJvMEIj4M4CDcgvMtL7z1pbhCmrjBMo/WGislric3/El5pyrbRwSvgqgbzahlWwQNycniBq0aQIf08b2lU89vTQbg6feQbfVf6orvM4l6jc931ZvtPTTStPTlppvylhVMYNPc9Q7xzI28bSyvM+60+rQm4AFqbRuf+C5cy//ozWzlYJuRRaeYUnZc/11FTqL6URyX/Xv+FNC6eEN9OzCcLob2dqxGPCnChnKrXZACYZkIKrhTRf5+C0eFl3Hj0I/rJ9U1ph/sbDmucfO33OozLNSmbhUb3a+8x/3IigUy3XKp48q2ZbeuTcGfTOgbz35gkWbGr5vLUgYdVkPPE5HLN0h/j/y9VPxWjcXWCdU2g12BleWnNA9lxPBUp3fMB85dtqpcTnhD5408Ip4S33DiX/8rzCkwhpFiDrcytVOEmDiGrqLzqc5xxw9PCAjli8Wd1PVUAqjPTp2EpVYBPQk9fVGbXgIbVCJVI89R9XmmeM2LSqBX8595sc3srU4ryG6JkDjQ4WzCssNURY9QfFGAmrJqKnIocwCXytYN71dd5KKwf6aCFt+yk3tUkZHVo6/F1QXM6d+sZ5MtF6lnrNWb52wrdK4nOj4R3zVnw3eNqklRZOCW+m2VKLnJ70/h5Ehwdptj2jQ0vV/pB+1qdjK25h1VkYURPYtMo0a21uecz8PHiywZCbZ6LDgjC2XxImDkw1ZNPqiauXtyw9cs+Qdw5Ua6OAnmDBt7ceRWZyjMfzvD8oxkhYNQk9FTnkog+d8aX519fCkCfItT3ABtXSlwLSycT5ebvrH8XzLPWYs6ySlF/vImNFAU8K75i34rvBu0nmWTid8aZ7B49/eUXNBZVPG9l57Axu6NKa6568AvzbD16LAV1aOzwDJQ1bfGSoZp5QLbOuu+5LAtINRt3FBt3fV5pnKs5fwGvrDmHptqN4+Y7uWDd1gE83rUZYepTeBa1NGu8cqNRGAJg7qqfog50YFarZ557IBGa5sZkFBViZAM8uWFqRQy76UA5fBDZoBfqs2XvCo2ABMxACEJ77Yj8elmk7bxPjVSJQ3fGP4g2a4nXaz+jQ0i+d8PMKT1o6ywHvmFcaX74MgtNqu/MCpRRlnxgViv/rn+ySa1gafW82ngpoAnp8c3mCfQCgT0orWeFALggsb2p/7vsrbYL1uC8ptduTfNQ8m4ZH3t2FDQfLHI5/+ee74u3MEZ5mjjH6GWq1UUBIbycdh2q4KxNoubEB1ls7SLNqAnorcugdDt4KbOARuieu3O0g/CVEhmB05hVIio3wicaMJ3+qgLOGNb5FCE7+ubAtG9sLGR1ailoCZ7QqEgnXP3OuHoD+oCkec5ZWNLdVnfAnv79HtybYW5pYd8a83Hm+sIK4W7ZSzYQ5cWCqePzJnC5IiAoT07OZ/buMmuf0+uZ6oiEFXDVselDaBPNujqdkd8aKgmOG+o/zbhoYHMeXL03LnrgTzB3VEy+tOWC6Dz6PD+9D/ZLw1tajmtfS+6544sbmK0hYNQGzhUlvBTbwTFLOi3ZpVZ2Ds783Jyk1/zY5nIMwpMKp1iKjVpFIev0JK3ZhUcA1iAoL1iVY8pizPtvzO9fv9GbUNg96/fa8aWp3Z8zL4YvJ3pOylUqLu3N6NgFvvNdGzXPO/uY8yAnwaptXXtw1MfMG5kwc2Al/vT7JUFP8usJS7nOVxpdVTctyDE6Lx8CurS3hg5+dlsAlrOp9V4xK8+VNyA3ABMwSJs3IW6mW88+Igeotk6i7+VMF3NHUKVUkcib3i0Lu2tbSZ65lzuIdZ6f+zO1odQQBb9nWYrG93s43a/Tk7M3J3ox7afm1mfle87jDRIcHwQZ1k727QobRlZ4E3DEx87gnCN83st1r95VwCUtaeMu0LKxnX+31bFya1fd6ESx4SrgrE/hjQRcSVg2mwc5gtzNEhwWpnhev0zTF6zMjvKyf7PoNb235Hz7ZrexLKlchKeP5PLy+7hc02JkhA1VukvI0KbYcRvm3GY0ggJWf5fObO3TyrMMzUZs0tRZzgee/OuDQr1bn+a8O4PpZG7Bm7wlVszYD8PQnP6P+ot2wexs9OXtzsjf6Xp74tRnxjqsJaAIv39Fd1udWzd/cl8wd1VO1Cpcag9MSMDm7MyKd1haz/IgFBYBR6PWvLCgu1zV+pOvZEx/uFY8rbbj0YMaaxYNgwZPDEz9as4RgMyE3AAPR4y9ZXaedN04Kj8+M2v2dByZPdOdLt6W7ldrGGekkVXm+3hSTrl6tkg38PmhGEBMRzPUs5288jPkbDyumOJOiJ5rbOWrX6qa40spaPLpit+Z55ecu4MbZmwy7r7AB8HTj4830TgI8pmI9Y16PX5s00bye2upa8EROAzDFZG8G7pqYleb2iTemYMrgLqZo/sxSAPDmBpX6vKoJVoC6C9jk9/cgpFkA99grKC53CMozcjy7g5IFzxM/WjU3Nm9m/NADaVYNgrfevEBNPX/6EJ5oRq37Syd63ujOCSt2Y0SPxnsaMWQXbz5iWvS0Hq2S8Ft4SvgZRUJUmKaWSEpJZS0efncXXv7aNSBGilI0txJC1K7Ss9ajmXcXnt+v525nONIX8cIbhasFA/DMsG6GTvZaJk6jylYKuOvXpqe2Og9akdOAdcy2zhQUl7v0l1JblbR3a/aWyM6bQGORgTwdPqV68CQDgRpvbz3qMg6E364035WpbJp41jM97gdjlu5w0MxOMng8G8GTOV3w95yuYrCjO6gJwVb0LSbNqgF46i+phdbkq/f+O3REd37+UwkW3HsNcr90vyShwOZf/lC8jw3AjM/3o0VoEE5UnNd9bS2tkhQ9Jfw8RaphCwyw6S7v+N/8Y+L/8wpP4taebV3OcQ4G0YIBeOqjn9EiNAh9Ol5KwaNHM+8pvHk99WKEUK0nCleN5786gIAAmyGTvloVJ+mYUNJERv2ZvH1gV758o4D7fm3uJr5XwyrCpxpyJWGl2kGgsb/khASlIMJbrkrAW98fVb2v2jP1pDKRHgVAQlQonhnWDblfFmrObc7jgMciKR1Tzu+4OxX/nNHrKuDLvOeAecGOVi7oQsKqAfjaX1Lv/QuO8udkK6msRcuIYIdkz2bQGKlch/v+s93huJKA5gyvSZy3hJ+zKcgd5MwpegVLKWrmLL2TS8X5C7jvP9vFVGMVNfVYuu2Y4vmeblSkyKWGMQo9CeDV4I3CVcOoCGi9Jk65MSa4gawoOM59X630bHoTzVsxHY5R8LqATX5/j0NBGOG7cv1bUlmLN7ccVb2e1jOVSx/Fu0nnVQBICyRcsDPVDCnOba48X68rgwvg+o57Gtmu5puthlXGs5GZFqxklXCG3AAMwNfpHXTfX6fyqay61mcDePL7e7hNLUom8QRJsAXvy+hsCnLHSV/JnOLJszQ6mlZINaYmqDrj6f0Hp8XLmnWNwNMNhoBWAIIUpe40IgLaXROnc3J2ATVzqjNmBXf4er40Gr0uYNJcpEZZ5aTPVCuDA+9cxhPgBjgWSODNkNLYnvNu/Xbnd9zTyHYt32wtth7+w6eBq1ZN4m80JKwagK/TO+i9f0BA4y7brOsbjdxLKI0Urb9oF329osKC8c3kSxVj9FSQUZvE9QjNwn31Vk3RwhcVzOTYeeyMx9eQCjhqC2F0WBAmD0pFy3D17BoCehPAK6EmqEl5MqeLZoEAT/pMj4lTgCeKXzhPCyXhIyosCJOzUzE4LUHzGs74ej4xEr3CpnOuW63CHrwIz5Sn7+UKQyih1ydeD+Xn6t367c7vOG/FP6VgR083T/M3HsH1szYYknXAXayyNpgJuQEYwJlzddz15s1AbwTzok3/w/j+yVjyXbHqeUpmPrXE1q0lVaCMQCniWGrecqlEJZlYeaOxeUxBevyT1DS4nk5qvtZMbThYZoipXYrSmKo4fwFz1x9CXHM+YdWdBPBKKAYgSCLqY5vzLeKNfRaleZ789/Sdx6sp8sRlQnAreH/Hry6mZXcT3/sjnrqAGWEJkApiPBkc9GZA8cR1SQ5hHMRwvjvOOL/jPC5gahYAIzZPpZW1mKzh/uANfL02mAlpVj1k7b4STFihXn7RDKSRowXF5XhmWJquiMzP9pzA5EGpCA8OlP1cy8yndK+nh2lro9yhrLpWUchzfvZ6TJ0CO4+d4U7T4wnu+kdJOVVdp5rv7/VRPREVZt4+9L/5xwyNgpXLPenMqbPKEf96E8B7kjPRWVPPq8l1d0F0x8TJu2AZISgpmZbleoEBGJqegILictVnrje/pi/xVDgwwhIgnaPNElaMLogwfXgajp8+p/t7Sm1R0wBLU5zJocflRwkh77OvaUpWC2dIs+oBPCYgMzSuSpGj4/sn49M9JzS1KkIw09z1hxTPUcvhphQgM3dUTwdtVHxkCMqq6gx5iWMjQjB19R6uc51NnTwTLe/CzbsYSKNwMzq0xM5jZ1BWXYtT1XUe+UcF2BqjzAXkgiacczkajdFRsDztlfZp6xYhKJP0F2/u0AY7w/wNh7F0azEqzl8SfvVE0zprBTM6tDQ0CEnufjxlNqXt4l2wjBCUhKhoqWlZbn4Q5sG3th7FW1uPqj5zT4KCvI1e4cA5121Gh5Ye5bJ+7Z4eDs/Q6sKKsK4AcCjLrQXPO66kAdbyo+UpnR0dFuQwZ+hByBJRVl3rVlYGHnisFs7tSG8b6XKOJ9kjzIaEVQ8wqo64HpQiR0sra/HGd8WYc08PTPngJ4/uoZW+Qkm4cJ4Unh7aDZPf36Manf/2g9fi6U/24WSV+mIMm3sR6TuPneEyOxqtIVNzU/AE5+sImi3nCGMzJxkzomD1tHfG8DSxYABvAvi8wpOY/vl+VMjkZPUkmtbs5NrumDi1ovgFjHKZcDYty80PSuNWeOZaQUHO49sqZCbHOLiF8DBtaFdxvOgp7CHHkHRHn2GeDA7eLIYiRVhXAOD6WRt0fY+3yIO76bvUBNrXR/XUvfkXstgoKZaM3IDxzDNy7YiUsb7JbRStkm/VK24ACxYsQFJSEkJDQ9G7d28UFBSonr969Wp07doVoaGh6N69O9asWeONZurG2/4hDXaGGZ8rl58EHPOvuYuwiKmZSnkW38Fp8ZrR+X1SWmHGCPWI0+nD03CKs1ypM7waU0FDpoQn5efcEVR5A4r0BE3YYExxBwFf+UeV19SL/5dbfOSS5096f4+soAp4Hk1rdnJtNRPnhBtTUPdngKG0RK9WFL9wnlloXVv6zOsv2g0NCvImeYWlqNWhMXe2PgHGBjHx9D1vYQhnVzNPEd5VvX6+PBo+tYIZY5bukC0nzsvgtHjd78rMNQexZq98lgi9WRm00Ap2VMpWUXVevYqmrwsfOGO6ZnXVqlWYOnUqFi9ejN69e2Pu3LkYMmQIioqK0Lq1q3P/tm3bMHr0aMycORO33HILVqxYgdtuuw27du1Cenq62c3VhbdNLvM3HFLdEbvjPC/Hks1H8NGu3x0Gt7s+Pc6mGeddckFxOeou2jE5uzNWFBxTcC1IwLKt6sFgSpw6W6dY8UcKjynIW+XnEiJDsXZyFno+l8d1vnO/K+WInXBjCpb/cByVbpqznPGVyVEtqEmuNCIPvNpi3vy7RifXVjJxzt94RPy/VBOiJED7Sqsmh/DMl+cfNTwoyBuo5b9VwrlfpFq/byb3537n5b4vCHVa5TmdNXpyZuq8wlLu0uF6MXqTq6fMubtlp/UKlqVVtfjnZ/tUC2QobcB4zPVSxN+09SjG9kvCxIGpDlXR3E2N5uvCB86YLqzOmTMH48aNw9ixYwEAixcvxldffYW3334bTz31lMv5r7/+OnJycvDEE08AAJ5//nnk5eVh/vz5WLx4sdnN1YXRtbjVWLuvRJePjydIF0ABT/wsnU2UUqRmh9YtgmW/f/2sDW5Pms6aZrUiA2qmIC0nfSOZNrQrgpu5b/Rwrpoj4NyvN6fHY9fxCt19a5R/lNr11SbXq6+IVvxs8vt7PPKR1lpIeSsSmeHvpXU9qWldzsSox5yqB0/nuWPlNdonWQyj8qNKx5NUISD0tNb1ef17lYqhyG3uosODFK0QRmDkJtedDQNwqew0j+XD3aDY8nP1ip8pbcDk+oNXUSQniHuarcIqhQ8Ak4XV+vp67Ny5E9OmTROPBQQEIDs7G/n5+bLfyc/Px9SpUx2ODRkyBJ9++qmZTXWLwAAbpt/SDZPf+QGA68RiA/D0wC4OyeXd4cK5c5j58W6EXDS+4o87fPtjMbLT4l3aY69pXHSE4/aaGtgvNoO9/qLDsQ0Hy2R/S9WZOji/lk++ux0McDnuLk++ux2hFxu1tc7tkrbdmUFJLcRzBKTfNxLhXmb398Y9x/HyHd3x5Mc/6/qeDcCMwd1gqz0P+5/HpM9i3a6j+Ne3vzhMxq1bhCg+a+e/5RZqqQAbUFfL3W96aR3ExGvy9O+T726H7a6rZN8F+0XH6VXuenLnNdgZfjx2Bn9U1yKuRSiu7dASgQE2rvbYAMz8eDf6/r9+LudeGx8K1J53eT+d26Z0XOl+gOM8p/aOy5EcEeBW/6m1H3Cdh+TOk5ufeK65s+gEyk9X6p6X1J5NZfml+e/Vu67Cq98UOQQSanHmdB0mL8vHKzLj8dr4UPF9lf4uufn1fJXrPKz1m3jGi3DetfGh6BBuU4xTkPse4NqfDXaGmR/vRrAH7/7Mj3djUFIk7BcbFNu9s+gEKsqrDFuD1JDrD+m44OF8VeM4eH1UT9RdtLv1btUFBgO2S5tjK6TEsjHGTMu4cOLECbRt2xbbtm1D3759xeN///vfsXnzZmzfvt3lO8HBwXjnnXcwevRo8djChQuRm5uLkyddVfF1dXWoq7vUGVVVVWjfvj0qKysRGcmvzXEXe00Niq7JMP0+BEEQBEEQZnPbLS+irtklEXnluD6maFarqqoQFRXFJa/5fZ7VmTNnIioqSvzXvn17r95frb48QRAEQRCEP+JJYLHRmOoGEBsbi8DAQBeN6MmTJ5GQIB+5lpCQoOv8adOmObgNCJpVb2ELC8Ntt7wIANj5z2yEBzc+0pr6i8h4YZ3Lcb1Ir8PLkgcy0C8lVtbHTel6CZGhuDujHeZtPKzrPv3/TEXi/HsBuPx+Pb/l3t5XYMX249xtcZdX77oKA7u2Vu0rtXZLz88rLBUDtLTMFTY0pkRxjuCU3qvgH4MwYt5WbnOZu7x4WzruuKYd17nVdReQ+eJ62c+WPJCB/1u+k/u+cuNE629p30ifld57G430GfKOF3febV6k72ZeYSle+OqAg0nZOVet0DYAivPD00O7iuNVru1q3+dFLc2b1rPTGj/O11Gbs6TnNjAmjvklD2Tg6iuiFd8BAcGPd93UAeI87E5/O4/5L/ee4HIr0xqPL96Wjn98uk9XW6QsG5uJ3skxiv2gZw5d8kAGqs5fcHB5kfueEeNLDa1286QWS4gMxWM3deZ2/RPuyduvD/fviMXf/Y/r2p5QFxhsSOo9IzFVWA0ODkZGRgbWr1+P2267DQBgt9uxfv16TJw4UfY7ffv2xfr16zF58mTxWF5enoMbgZSQkBCEhHjDm0Qem80mqssDwsMR8OdgD2h2Ufa4XqTX4WXMyn2KOdKk11s2tpdDJCkALP+pjDvg5o8LAQgID3e5rnDM+fezgAvcv6VdQgzqmplfa3nSp0WYOypMta/U+kB6/pBrO2JuaBhXZOrro3piiEygl/ReQRERmHbH1Xjk3V0AzKuQEhsXLfaZGmv3lWD65/sVn8UPJ2p0jVW5caL1t7RvpM8qo0sbRMcc0hUsFhnaDFW1xlhGpM+Qd7y4827zcqYhEAHh4Y0BKB8eaBw7knv9dh5gTvcOCA9HXuFJ2TYdr2F4+MMDWHR/GHLSE2Xb7txfnnB/nyvw7g+Om1WtZ6c1fpyvozZnCaw/Wo2X1hwQj49ZuQ/xkSGqv1FMEXXH1QhqHnHpvm70t/OYj4tryXUNrfEYGxftdj8lRoUis1tbBATYFPtBzxya2a2tojKHp4+MQqvdakFvD/VLQnZaglj8o+5TvhSSwj15+zWzWzt8sP+0oVkaBCHcObBOrTCQLzDdDWDq1Kl488038c477+DAgQN45JFHcO7cOTE7wF/+8heHAKxJkyZh7dq1+Ne//oWDBw9ixowZ+PHHHxWFW0IenhxpmckxuLVnW/RNaYXAAJtqjj459FTAWbuvBNlzNmueJ5gdHuibpJn3VJqv1ZN9n5E5HHPSE/H9kwOxclwfvHrXVeLxZWN7OfytVVVFer1F91/TWBhBgvS3C7j7DHiSwwtRt2rC4OLN5u/4leDNLbpsbC+8PqonpmR3RphCqWF34E2wb1RuRS2mfvAT1uw9oRixLndMLeqZNxetUeVR5awqWs9Oq4yrO0x6f4/LmNcq58yTX1fpXXWe15wRMtBoveta41Err7QaVtG0qSG0TqmcuDvMHdXTZR5OjArF4vuvwTPDrxTXUTmc539ntEq+Cutin5RWYgUwd5h4YwqiwxxzeCf8+Rt2/nMwVo7rg9dH9cTKcX3w/ZMDLSOoAl5IXTVy5Ej88ccfePbZZ1FaWoqePXti7dq1iI9vXKyPHz+OgIBLMvN1112HFStW4J///CeefvpppKam4tNPP7VcjlWroGQ2czdHGq8QBehboPWkFJo+PA3BzQI0KwNJK8EIx/UuVWbkcAwMsKFvSivU1EeJph1Bc+1OZojGvJkJLqmg5FIQ8T4D6Xk8CbeNSNNjdvUcntyimckx+O6XPzB33S+G/B65Z6gmVE1+fw9CmgV4ZRH452f7UH6OPwXRks1HNPOdCmlserSPcvk8r/CkWyl+5JCb07Se3ZilOzyu8c6D3LhxtlLxzLdK76rzvCaFt+KV1v158krLMSW7s6UEGCUEreDgtATM33AYS747gpp6/eWOpQxOi8ctV7Vxq3Sq1vyvJ893TnoiFt9/DZ76+GfdKcYevbETpgzuovgbfJ2eSg2vlFudOHGiomZ006ZNLsfuvvtu3H333Sa3yn95MqcLEqLCUH62zqFOvDNm5EjTI+QIvLTmALdgIM1nqie59dxRPV3y0xlZ5tTXCAKwgFxgn9wzSIwKxYgeifj8pxKn3H38QqOnufqkqC3ERiIVHqSCvVGCNyD/DHlyMgobSDO1rAzQJagC8vmV5WhMY+MqrHqa45YHrWenpfVUw5P+yEyO0RWXIPeuCse1FAaCtcXTpP16FBMAkBAZgokDO7l9P7NR2jBMHNgJKwuOawqrPFp553mYl4Lick3ljp4834IC48kP9+LDXb/paou7v8HXeEVYJYxl1toiJEaFYmi6fNCZM+7mSHtdZkJ1RzOmx49Q6YWVm4ikApvcrjejQ0sszz+qKtA3JdR2/n/P6calmZXDqBx7cguxdIEoKC5Hz/bRDn+7W79eKjxIx8nOY2fcWuAbtcIh+Nc9PXHqbJ3iMywoLufSTv5w5LRhWkhvo5TUXW2pjw4Lwn19rsACToFY6fpaz07ahgY7495Qu5v4XaCguFxXxTLpu/rbmRpR46bHPUiwtki/bxY2ADNGXGlp87/ShqGguJxrzdp57Axu6OJaVdMI9Gr9l43thT+q61THRWCADVmdY3ULq/4KCasW55t9pbLHSytr8dbWo1zXcLdiiJzwo0fIMRIezYXcjnFMv2T85/ti06qMqVXE8gVKu2YezawS7o6fh/t3RErr5ooTrrPZeMzSHZCuhWaYdXlKpTojNGnGiCvRr1OseFzuGfIK9vn/O+V2VTg97i4xEcE4c67eEI2ntHJZ3UV9JtWK8xcQEy5foU4vvM9u57Ez3Cl3dh4741GVvjFLdygGtSoh5y4EuG7glIRgpe87f89TP94AGzB/tHalJ6vC+066MzfIobRm6xlfvG5jvip57Qv8Ps9qU2bN3hN4bPVPsp+J1XxUNrpG5EgTJkRpIJY/Ifh4AcpBDdOGdlW9hpYPoloQW1OAN6jDmT4prTDsKscFTrpwygWvOK+r0s+NCKDRExQowBMwI3D01DnOq7r/HiVEhWLhvVer9onw7r9wa7qHd7t0PcCz4JqYiGCDNh9899cjfBghqPAEtWqRV3jSIRB1zNIduH7WBt3XFL43c02hbGCrHpcHOwNaRhiz0TATpfmBV6BzZ25wRm3NlsNTVyB352Z/hIRVi7J2XwkeXbFb0+dS6XNv50hz1gZI/46PDHE7+tUIlCLqAW0fMV4fRKOjkK0Ej8Avh7MJ33kh1ou7C7dzmxJlxoESzwzrxh0Vu3ZfCV5bd0j1HEGI1OszJmQxEKJ0h17VRrNPpg9Pw9CrlMe+HvQI7MrXCNPMNmLTGGABNqAZ53wW1yKE+700QlDhzZigxmSZDZyWEKwk8JRU1mLJd8WyGr3J7+9R1ADKYYVym1oozQ9GZVDQgnfNluKpskOYm3lv6c/rFAmrFkQIAvEEIxYXXpzTUo1ZusPhb2GBcp4spFH9ZqOUUkrLR4zXB7GguNyoploSNYFfirSPnTdJcguxXjzVXkkFbx5iW4Rwbfb0vLPTh6ehT8dWulKzOaeZA7Q3YcK7L4z99x7q7ZK2Ru2ezgKy1lyipeXNTI5RfN+iw4IwJTsVC0ZfoypU2Bkwd/0h1d8hcOZcPbdW8cy5elUrFS+ezgdqacbkhGB3fW0ZoEsD6C/mZrn5gXez7YlSx5M121NlR056IqZkp3Kdu/PYGbfv42tIWLUgRkRfPzOsm9cEVbn8m9Ko3MFp8fK5Qv8UqPVGpbqL4NLgbJpWg1ejYJTmQU1D7Wty0hPxzDB1QS8qXFmIMOKXGKG90jO58y7SvIFbE25MQU56omZOYwbgiZs6i38rjQXeTVhggA39UmPx8p3dYQOfS4ycgKyF0oZUzsIjFYZ3PjMYk7I7Y+hViVhw79Wa7k08n0394CdFraLzZmfqBz9pasQe6peEiTfyRcMbrYlUEoK1NtNq8Lw+RpXb9NY8pjQ/8G623cXdNdsoZUdSbIT2STDOL9cXkLBqQTyd6GwAnv/qgOkThFoaIOeoXOmiatWkw3LwCitGaB7kNNSemr2NpMHO8PxX6tqDkEDzpxRhgvdESzBxYKqq+4neRZp3EejQ6tKiorVJe+KjS8EVamNBzybME5cYLeSSpjtbeJznpFuuauMiDLeMCFEVpBgaA7bkiJf0qdrs5yzMqJ0bYAMW3tuY+F0aYKeGWZpI57XBTPO8Ua5kvAVhjEJJAFTa2BmBp/3g6fe96ZfrKygbgAXxdKIzI7+qHLzaJCEliD/mdxMqiyhpL6QR0p4gaKidF03BrOUtlw41eLQHJ724c/dESxAYYMOMEWmqpWz1LNK8i4CexcJZYDNqLCilPfLUwqGVNF0o1yugFEGvd+FWyqmrhN7NjjTA6My5OtX8ze7MB3qUCs5rg5nmeSPKbSrNa95AbhypZVDwBE/74eipGo++L/jlKmW9EfDUL9eXkGbVghgV4efObk2PGZpXWPBn0wNPCVpPNQ88GmorBHHpHU9muzF4qiVQ0zJOuDEFg9P48hgDlwK3tEaB3c7cfiZGjgV3XGL0XNfZfUDJXUjOx1Dvwu+OuwKgb14qq67F2n0lmMARQKNnPtBbhtpZCNYq0+kuE29M4bJ8aZXe5SnAYdYc4U0/W0/X7LnrfvEoM4A3/HJ9DQmrFsTd6Gtn9L6ses3QZmiTrIieyiLuoKWxtEoQl97x5BxoxzOWeeZSYeF21hK44+8rmAanZHdGlCRoZ/7GI7pcMHjf2b++86NHrh1WGQt60LsZ01r4jcogomdeio0I0RS8AmzAgnv5td5KArwzauZ4ns20O/TrpF3gQEnQFoQuXj9Oo4N+eFx4jBaQjVizZ6456FEbzPbL9TUkrFoUpYGXqCPHoh5TlB7NhwCvNsmfTQ/OLBvbS1c2AR68HcTlLu5oD5zLXyp996F+SVg5rg8OPn8zVo7rg7/2S5I9X2nh9iQ/ZV5hKeau+wWVTn6QejMP8C4WwnU90aT4eizoQe9mTG3h15NBRGt+lM5LWufCBk3BS08+Uj1lf7UyuyjNQYlRoRh+Fb91AOBfO9QEbSGAjXeMbjhYpviZ3g0oj5+tlpDtLmprtlY/MMDtwjTObTDLL9fXkM+qhZH6ljn7gAUE2PDIu7tcqtm44xSvpfmw4VJNbinCoiLXDik7j53RVYrQyvBWFtGDN4O4PIG3v6VIz5lzTw+88k2Rw6Iv57PYN6UV+qa0QmZyjEv9c6kfnbSClFxdekEonDuqp2L79Ix9nvErvLM/HDmNCSt2yQYCCdf1RJPi67GgB3c2Y8LCr9T//TvHcV1TaZzKzY9qc+mps3wuA7y/ldff/5lh3TCmXzL33OlclvrLvSfwxV6+fKq8awePoJ37RSFm392D677/zT+G61JauQjjcj7OWi4PWn62aj60k9/fg5BmAYb5g0vXbD394Clm+eX6GhJWLY5SUJLWZK7nhdOj+ejRPoqrHdIgBHdKETZFlMooajnH8wZt6K1P7g5K/c1Dy4hgfP/kQMUAHLl7KW3WnFETNtWEQj1j3/k9VOpPYTOpFLEuXFeqSeEV/tXGAm+ZTm/j7mZMrf95ygXPHdUTL3990KV/Bdcd6TXkzpXOpflHTrv1G5Tg9ZflzfMr4FyWWs+mhnft4H1nwBo3ozzzxIzP9ztsCJWESmdLDeAqoCs9L14hm3djqoTcmu1Pm0urQsKqH6NnMVdDn+YjyuW4tB15haV4e+tRXZHMVl1kjSSv8KRD8m5nAV5JY6lHU+6tTYHzuDt08izmbzys+b0/qut0Z4TwNIOElnnNXRcMrch2d8z0WgKr2ljgjbT3BZ5sxjzpf2l2Aq2sB1qZDIzaUArw+st6KuTwRIlHhwVhwX3XoE9HvgA13rF96lwdpg9Pw8N/ZtxQo7SqDvM3HMak7FTulIgCzgK6Ep5sTD2FZ/zER4Ya4grQVCGfVT9HKfpWD0aYoQUN4dcKJfyUIpmtnlvUKLTKKCr5Oyn5qyn5VxlRn5wH6bjjzT1pxUA7d8Y+j3+3XiFDLkep86usNBbc8Tf3Jjw+qGaVhdaT9UBtLjX6N2j5+xuVjF+r3TYAL9/ZHf06xXK3Xc87k5OeiIf+9D/X4rV1v2DtvhJDiuLI4cvYAHf8sK1cIMYXkLBKcEXf8kycegMptBZZTx3eeTFrUtBKOu4swPMWTlArseiLVFfeqr1tBM79q3fs80a28wgj8RLhPSosCJufuNGh74VgM62x4Iu0Z3qvp3czZkWM/A3eFOCNfvZ635lsHSngcr8oNE276OvYAK1+kGr8PQkYbaqQGwChGjijZ+LUs3PlCWzxNJUHD1rmeXdxNssq4Wx64jF7apVY9FZRCAGt8cMk55mFmgldzX9a79jn3ZDtPHZG85nUXbSLx6TturVnW/E4z1jwtmlTaWznFZ50aLszRrktAb5zHTLyNxgZd+DNdmsFWjI0BoY5u1DwaEtLKmtRzhnMJsDb/0a7crgDrx+2WsCov2zujIY0qwQAY3bfenauPIusWTts6UI3ScM87w68+ROl6DE9WTHVldL4iTcgHyYvSkuVkv+00L96xr6eZ6903ajwxnyuzgFY7ow7b48FnpRFahjhtmSU65C7FhUjfoOAN8tQG91utTRtz391QOwPqRaZh5iIYF0p8nj7nycXqlnuKM7t0OoHqxeI8QWkWSVEPN1969m5frn3hKFt54VH46mWrksLPfkTpegxPfnanKWE3PjhKX9pBHLR3EplMeVSUvGOfb3P3vm6sREheGz1TwCUU1rpiUj25ljwVjS1GlplidXSlDlfxyoBaf5YhhpoHNt2O/DoCtcAKmctYE56IqZkp+K1dYc0r5sQFaY7RR6v1lEtm4kRBV7MxttWMytBmlXCAU9233r8sHyRykOPxtPdSkF6gwPcCaTQKrFoVHCGOxipvdHD4LR4Bw3VM8O6qZbFlOtfnra7498tvW5AgE3VYqB33Bnlb86Dryut8fjn8rgO5RWetHRAmr/QYGd4/qtC2c/ktIATB6aqVh6TjlUlza3SdKJH6yjVZs+6s7t4PCosyG80lv5UFMQoSFglDIXXpOqtkooC7mo89U4Kes53N5BCrcSi2dHVVkYqFMZyZh7Q27+eBsYYbbb3ZqCOr91PjHIdemnNgSZlZvVV1Lg7lclmjEgTsxBIkRurzi4S7mxAlQgMsKHyfD3m5P0iHvOnIKbLMW8rCauE4fD4YRlVUpEXd9Oh6J0U9CbhdtdZfmDX1oZf0wx8tZCaaR73xL/bjHZ5K9Le1+4nRgnBvMGJ/oAvU/95UpmMd6yatQG1ero3b1hK/A3yWSVMgccPSysaVpq/09OIX70LndS/tu5iA/f3zEjC7Yyc3210WBDG9kvCxIGpltGo+tIv0OzIX3f9u81ql7vtcd5MqKUX4xnbCZEhpi2k3tQm+YOZVct/1+xNqxmVycy4nzPulhf3Np5k5mmKkGbVJCihLx9KWlgAhmoM9Cx0nkwKZiThlqKkEag8fwFz1x1CXqF8UQZv42vNhTfM4+7455rZLr3tkdPKSf/W03aB2ot208agN12HrG5m9VV+XSme+Eu78+6cOVen6LOqdT8pvva95kGuSIjVrGbehoRVE2jKVZnMEMKdJ668wlLDBR3exPWA66Sg9zebZZa1wgLFg7faqdUvVk1Eb4V2KW0m5GqvSxHaLqTfcqay5oJpmxFPXIekY6NleJDfm1mtIHB501967b4STFixW9Vnlfd+vva95sE5YNTMlGb+ArkBGIyvTTNm4g2zrlkmGq1E1gDwUL8kZKclOJik5H6zWiS+gJFJuAV8WdtaD95oJ+9YNKMfjMCX7eKtva60mRicloAZnxfCqPRbenDHdSivsNRhrJypcW034F9mVqsIXN4obMATHBtgA+aPNj4fuC8xI6WZNwtoGA0JqwaiR9Dyt8HiLSFcj6DTo32UrmsrTaxKArfSb9bSPgkYPdlYZYEy6v7utlPvWLRqHktftYs32HDnsTO4oYtrIF9BcTl3+i0zfp+SoJ9XWOpi0YoOD0KFgnDqjBmVo8zCSgKX2RsvnvFqZ0DLiGCu6+nxG9cTr+AP+DKfsKeQsGogO4+d8QvNl168KYTrE3T0CasA/8TqqfbJDKy0QBlxf3fa2ZQ3hN6C9x37o1p+U2aFTZOzoK+0gdESVGMigvDMLVciIdIaGnderFA6VIqZGy+z0r3xlFj2Zllfo++VV3hS9ri/WnnJZ9VAlCZ3Z3yt+dKLN/2j9Ag6ZpZL1KN98hZGJoA3MwDQzET1VvDV83d437E4hVRBVts0uZtDGQDKz11AQmSoVwtYGIE3/UV9ja/SvXkz9sToezXYGV5ac0D2MyvFN+iBhFUDUZrcnfG15ksv3tSk8Ao6Z87VmzqReKp9MgOjFiizJ2EzF1IraPX8Hd5gQ6U0Vt6smsWDuzmUBfx1rFghUM8bmDXe1PKBezObiRn3mr/hUJPKJwyQsGooGR1aWmoSNwpvalJ4BJ0RPRIxYYW5Ewnvbzl1ts6ru1NPFyhvTcJmLaRW0+r5IzzvmHCeu9/3plbPU2HTn8cKTwEWf8fb6d68mXXFjHut3VeC19Yd4jrXnzZqJKwaiNUmcaPwtiZFTdBZcO/V+PynEtMnEl7t06y1RV5PS+buAuXt1FdmLKRW0+r5K0rvWDxnnlIrafXcFTabylhxJ2epv+HN8eZNVyOj7yXM8bz400aNAqwMxhupPLyNHod0o1AKhPJW+iaeVFcCvnBYdyegwRepr4wOvPDFWGyqyL1j6W0j0X3Gt25/3xdBSjyVtZyhseJ/eGu8edPVyOh76XGJ8beNGgmrJmCVSdxIfCGEywk63pxIlH6zM/4Shd5U/D2b4obQVzi/YzX1Fz36vi/Q2sAwwCWFFY0V/8Qb482brkZG30vP3O1vGzUSVk3CCpO40VhBCPe2z6Lwm5dtLcbzX8lHVwL+kZasKfl7WmEsEtZBawNDY4XgxZtpwYy+F+/cPSW7s99t1EhYJXThayHcF/kFAwNsiOXM9GBlraTVcjN6iq/HImEttDYwNFYIHrzpamT0vXhcYhIiQzBxYCcPWu0bKMCK8Ct8FcTWFLSSTTUAkCAELodgI8J8vBnQZeS9tOZ4G4AZI670y/eCNKuE3+ELn8WmopUkf0+CIAhtvOlqZOS9muocT8Iq4Zd422exKUWhk78nQRCENt50NTLyXk1xjidhlfBbvO2z2JR2rOTvSRAE0XRpanM8CasEoYOmuGMlCAFpQYiC4nJkpcbR2CYIwueQsEoQOmlqO1aCABrLNE7/fL/495ilO5Doh1YDgiCaHpQNwEOcNRHerBNPEARhBGv3leCRd3fhZFWdw3GhOps3ywkTBEE4Q8KqB6zdV4LsOZvFv8cs3eH1OvF6IeGaIAgpQj1xuZlAOJb7RSHNFQRB+AwSVt3EHzUR/ihcEwRhLlr1xKXV2QiCIHwBCatu4I+aCH8UrgmCMB/eqmtWrs5GEETThoRVN/A3TYQ/CtcEQXiHplCdjSCIpg0Jq27gb5oIfxOuCYLwHkJ1NqUEVTYAiX5QnY0giKYLCatu4G+aCH8TrgmC8B5a9cQB/6nORhBE08Q0YbW8vBz33XcfIiMjER0djYceeghnz55V/c4bb7yBG264AZGRkbDZbKioqDCreR7hb5oIfxOuCYLwLkJ1toQoxzkgISoUi+6/hvKsEgThU0wrCnDfffehpKQEeXl5uHDhAsaOHYvx48djxYoVit+pqalBTk4OcnJyMG3aNLOa5jH+VideEK5LK2tl/VZtaFyUrCJcEwThfag6G0EQVsUUYfXAgQNYu3YtduzYgWuvvRYAMG/ePAwdOhSzZ89GmzZtZL83efJkAMCmTZvMaJah+FOdeH8TrgmC8A1UnY0gCCtiirCan5+P6OhoUVAFgOzsbAQEBGD79u24/fbbDbtXXV0d6uoupWOqqqoy7Npa+JMmwp+Ea4IgiKaGc0GWrNQ4S64VBGFFTBFWS0tL0bp1a8cbNWuGmJgYlJaWGnqvmTNnIjc319Br6sGfNBH+JFwTBEE0FdbuK8H0z/eLf49ZugOJpCggCG50BVg99dRTsNlsqv8OHjxoVltlmTZtGiorK8V/v/76q1fv728IwvWtPduib0orElQJgiBMhAqyEITn6NKsPvbYYxgzZozqOR07dkRCQgLKysocjl+8eBHl5eVISEjQ3Ug1QkJCEBISYug1CYIgCMJTtAqy2NBYkGVwWgIpDghCBV3CalxcHOLi4jTP69u3LyoqKrBz505kZGQAADZs2AC73Y7evXu711KCIAiC8CP0FGTxF3cygvAFpuRZ7datG3JycjBu3DgUFBRg69atmDhxIkaNGiVmAvj999/RtWtXFBQUiN8rLS3Fnj17cPjwYQDAzz//jD179qC8nCorEQRBEP4FFWQhCGMwrSjAe++9h65du2LQoEEYOnQorr/+erzxxhvi5xcuXEBRURFqamrEY4sXL8bVV1+NcePGAQD69++Pq6++Gp9//rlZzSQIgiAIU6CCLARhDDbGmJw7jd9SVVWFqKgoVFZWIjIy0tfNIQiCIC5TGuwM18/aoFmQ5fsnB5LPKnHZoUdeM02zShAEQRCXM0JBFgAu5bmpIAtB8EPCKkEQBEGYhFCQJSHK0dSfEBWKRfdfQ3lWCYIDU4oCEARBEATRCBVkIQjPaHLCquCC682yqwRBEAShxZVxQbgyLggAcO5stY9bQxC+RZDTeEKnmpywWl3dOAG0b9/exy0hCIIgCIIg1KiurkZUVJTqOU0uG4DdbseJEyfQokUL2GzeMbFUVVWhffv2+PXXXykDgUWgPrEe1CfWhPrFelCfWA/qE+NhjKG6uhpt2rRBQIB6CFWT06wGBASgXbt2Prl3ZGQkDWKLQX1iPahPrAn1i/WgPrEe1CfGoqVRFaBsAARBEARBEIRlIWGVIAiCIAiCsCwkrBpASEgIpk+fjpCQEF83hfgT6hPrQX1iTahfrAf1ifWgPvEtTS7AiiAIgiAIgmg6kGaVIAiCIAiCsCwkrBIEQRAEQRCWhYRVgiAIgiAIwrKQsEoQBEEQBEFYFhJWPWTBggVISkpCaGgoevfujYKCAl836bJhxowZsNlsDv+6du0qfl5bW4sJEyagVatWaN68Oe68806cPHnShy1umnz33XcYPnw42rRpA5vNhk8//dThc8YYnn32WSQmJiIsLAzZ2dk4dOiQwznl5eW47777EBkZiejoaDz00EM4e/asF39F00KrT8aMGePy7uTk5DicQ31iLDNnzkSvXr3QokULtG7dGrfddhuKiooczuGZs44fP45hw4YhPDwcrVu3xhNPPIGLFy9686c0GXj65IYbbnB5Vx5++GGHc6hPzIeEVQ9YtWoVpk6diunTp2PXrl3o0aMHhgwZgrKyMl837bLhyiuvRElJifjv+++/Fz+bMmUKvvjiC6xevRqbN2/GiRMncMcdd/iwtU2Tc+fOoUePHliwYIHs56+88gr+/e9/Y/Hixdi+fTsiIiIwZMgQ1NbWiufcd9992L9/P/Ly8vDll1/iu+++w/jx4731E5ocWn0CADk5OQ7vzsqVKx0+pz4xls2bN2PChAn44YcfkJeXhwsXLuCmm27CuXPnxHO05qyGhgYMGzYM9fX12LZtG9555x0sW7YMzz77rC9+kt/D0ycAMG7cOId35ZVXXhE/oz7xEoxwm8zMTDZhwgTx74aGBtamTRs2c+ZMH7bq8mH69OmsR48esp9VVFSwoKAgtnr1avHYgQMHGACWn5/vpRZefgBgn3zyifi33W5nCQkJ7NVXXxWPVVRUsJCQELZy5UrGGGOFhYUMANuxY4d4ztdff81sNhv7/fffvdb2popznzDG2IMPPshuvfVWxe9Qn5hPWVkZA8A2b97MGOObs9asWcMCAgJYaWmpeM6iRYtYZGQkq6ur8+4PaII49wljjA0YMIBNmjRJ8TvUJ96BNKtuUl9fj507dyI7O1s8FhAQgOzsbOTn5/uwZZcXhw4dQps2bdCxY0fcd999OH78OABg586duHDhgkP/dO3aFVdccQX1jxcpLi5GaWmpQz9ERUWhd+/eYj/k5+cjOjoa1157rXhOdnY2AgICsH37dq+3+XJh06ZNaN26Nbp06YJHHnkEp0+fFj+jPjGfyspKAEBMTAwAvjkrPz8f3bt3R3x8vHjOkCFDUFVVhf3793ux9U0T5z4ReO+99xAbG4v09HRMmzYNNTU14mfUJ96hma8b4K+cOnUKDQ0NDgMUAOLj43Hw4EEfteryonfv3li2bBm6dOmCkpIS5ObmIisrC/v27UNpaSmCg4MRHR3t8J34+HiUlpb6psGXIcKzlntPhM9KS0vRunVrh8+bNWuGmJgY6iuTyMnJwR133IHk5GQcOXIETz/9NG6++Wbk5+cjMDCQ+sRk7HY7Jk+ejH79+iE9PR0AuOas0tJS2XdJ+IxwH7k+AYB7770XHTp0QJs2bbB37148+eSTKCoqwscffwyA+sRbkLBK+C0333yz+P+rrroKvXv3RocOHfDBBx8gLCzMhy0jCGszatQo8f/du3fHVVddhZSUFGzatAmDBg3yYcsuDyZMmIB9+/Y5+NgTvkWpT6R+2t27d0diYiIGDRqEI0eOICUlxdvNvGwhNwA3iY2NRWBgoEuk5smTJ5GQkOCjVl3eREdHo3Pnzjh8+DASEhJQX1+PiooKh3Oof7yL8KzV3pOEhASXoMSLFy+ivLyc+spLdOzYEbGxsTh8+DAA6hMzmThxIr788kts3LgR7dq1E4/zzFkJCQmy75LwGeEeSn0iR+/evQHA4V2hPjEfElbdJDg4GBkZGVi/fr14zG63Y/369ejbt68PW3b5cvbsWRw5cgSJiYnIyMhAUFCQQ/8UFRXh+PHj1D9eJDk5GQkJCQ79UFVVhe3bt4v90LdvX1RUVGDnzp3iORs2bIDdbhcXBsJcfvvtN5w+fRqJiYkAqE/MgDGGiRMn4pNPPsGGDRuQnJzs8DnPnNW3b1/8/PPPDhuJvLw8REZGIi0tzTs/pAmh1Sdy7NmzBwAc3hXqEy/g6wgvf+b9999nISEhbNmyZaywsJCNHz+eRUdHO0QFEubx2GOPsU2bNrHi4mK2detWlp2dzWJjY1lZWRljjLGHH36YXXHFFWzDhg3sxx9/ZH379mV9+/b1caubHtXV1Wz37t1s9+7dDACbM2cO2717Nzt27BhjjLGXX36ZRUdHs88++4zt3buX3XrrrSw5OZmdP39evEZOTg67+uqr2fbt29n333/PUlNT2ejRo331k/wetT6prq5mjz/+OMvPz2fFxcVs3bp17JprrmGpqamstrZWvAb1ibE88sgjLCoqim3atImVlJSI/2pqasRztOasixcvsvT0dHbTTTexPXv2sLVr17K4uDg2bdo0X/wkv0erTw4fPsyee+459uOPP7Li4mL22WefsY4dO7L+/fuL16A+8Q4krHrIvHnz2BVXXMGCg4NZZmYm++GHH3zdpMuGkSNHssTERBYcHMzatm3LRo4cyQ4fPix+fv78efboo4+yli1bsvDwcHb77bezkpISH7a4abJx40YGwOXfgw8+yBhrTF/1zDPPsPj4eBYSEsIGDRrEioqKHK5x+vRpNnr0aNa8eXMWGRnJxo4dy6qrq33wa5oGan1SU1PDbrrpJhYXF8eCgoJYhw4d2Lhx41w22dQnxiLXHwDY0qVLxXN45qyjR4+ym2++mYWFhbHY2Fj22GOPsQsXLnj51zQNtPrk+PHjrH///iwmJoaFhISwTp06sSeeeIJVVlY6XIf6xHxsjDHmPT0uQRAEQRAEQfBDPqsEQRAEQRCEZSFhlSAIgiAIgrAsJKwSBEEQBEEQloWEVYIgCIIgCMKykLBKEARBEARBWBYSVgmCIAiCIAjLQsIqQRAEQRAEYVlIWCUIgiAIgiAsCwmrBEEQBEEQhGUhYZUgCIIgCIKwLCSsEgRBEARBEJaFhFWCIAiCIAjCsvx/nVZurGpoUb8AAAAASUVORK5CYII=",
|
|
"text/plain": [
|
|
"<Figure size 800x200 with 1 Axes>"
|
|
]
|
|
},
|
|
"metadata": {},
|
|
"output_type": "display_data"
|
|
}
|
|
],
|
|
"source": [
|
|
"#mn = np.min(np.diff(ts_a[idxs_a]))\n",
|
|
"plt.figure(figsize=(8,2))\n",
|
|
"plt.stem(np.diff(tsa_a) - np.diff(tsa_b))\n",
|
|
"plt.title('between annotations: interval deviations'); None"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 4,
|
|
"id": "1e19feb7-b841-4163-8276-8ccfe40cac2f",
|
|
"metadata": {},
|
|
"outputs": [
|
|
{
|
|
"data": {
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAArMAAADcCAYAAABj2im9AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABj80lEQVR4nO2dd3wU1fr/P7sh2SSQAqRsQg0hAiEUDYK5KCgEE0DEK0pRvwIiKBClKpf7kyYqduFSxHIFlSqiIIpIqAqEDgIGEDCAQhIgIQVIIdnz+yN31i0zszO7M7Mlz/v14qWZPbNzZk97znOeomOMMRAEQRAEQRCEF6J3dwUIgiAIgiAIwllImCUIgiAIgiC8FhJmCYIgCIIgCK+FhFmCIAiCIAjCayFhliAIgiAIgvBaSJglCIIgCIIgvBYSZgmCIAiCIAivhYRZgiAIgiAIwmshYZYgCIIgCILwWkiYJTRl6dKl0Ol0OH/+vPna/fffj/vvv99tdSLkM2zYMDRv3tzd1fBobty4gaioKCxfvlxSeZ1Oh5kzZ6pbKQ+Dmw8OHjzosKztPHH+/HnodDosXbrUfG3mzJnQ6XQq1NS7cGV86nQ6ZGRkKFaXHTt2QKfTYceOHaLluLa7du2aYs92lYKCAtStWxcbN250d1UIB5AwSxAEoQLz5s1DSEgIBg8ebL62ceNGnxJY33jjDaxbt87d1RDFG+pIuBehcdmwYUM8++yzmDZtmvaVImRBwizhdjZv3ozNmze7uxqEDD755BOcPn3a3dXwWG7fvo158+bh2WefhZ+fn/n6xo0bMWvWLN57ysrK8Morr2hVRUXQUlCUMk+88sorKCsrs7pGwizhCLFx+fzzz+Pw4cPYtm2bxrUi5EDCLOF2AgICEBAQ4O5qCHLz5k1Z5cvLy2EymVSqjTUmkwnl5eWaPMsSf39/GAwGzZ/rLXz//fe4evUqBg4cKPmewMBA1KlTR8VaeTdS5ok6deogMDBQoxoRtYE2bdogKSnJypyF8DxImCXcjq0tHGdj9dVXX+H1119H48aNERgYiJ49e+Ls2bN29+/btw/p6ekICwtDcHAwunfvjt27d1uVuXDhAsaMGYNWrVohKCgIDRs2xOOPP25luwv8bcO3c+dOjBkzBlFRUWjcuLFg3bm6rlq1Cq+88goaNWqE4OBglJSUSK4b9z2dOnVCYGAg4uPj8dFHH/Ha/3H2bMuXL0fbtm1hMBiwadMmAMClS5fwzDPPIDo6GgaDAW3btsVnn31m96z58+ejbdu2CA4ORv369dGpUyesWLHC/HlpaSnGjx+P5s2bw2AwICoqCr169cLhw4fNZfhs8m7evIlJkyahSZMmMBgMaNWqFd59910wxnjfYd26dUhKSjLXlXsPR1y5cgUjRoxAdHQ0AgMD0aFDB3z++edWZTh7ynfffRcff/wx4uPjYTAYcPfdd+PAgQMOn1FYWIjJkyejXbt2qFevHkJDQ9G7d2/8+uuvkuq4bt06NG/eHPHx8eZrw4YNw8KFC82/AfePw9Zmlmv/33//HU899RTCwsIQGRmJadOmgTGGP//8E/3790doaCiMRiPee+89u3pUVFRgxowZaNmyJQwGA5o0aYKXX34ZFRUVDt/hzJkzGDBgAIxGIwIDA9G4cWMMHjwYxcXF5vrevHkTn3/+ufldhg0bBkD6eOO4desWnnvuOTRs2BChoaF4+umncf36dasyUmzrbceMUB23b98OnU6Hb7/91u47VqxYAZ1Oh6ysLN5nHDx4EDqdzq7PAcBPP/0EnU6H77//HoC0sSSHCxcu4OGHH0bdunURFRWFCRMmmJ/pyCZV6vjkWL58OVq1aoXAwEAkJyfj559/tquLnDaWyrVr1zBw4ECEhoaiYcOGGDduHO+GfdmyZUhOTkZQUBAaNGiAwYMH488//7Qq88svv+Dxxx9H06ZNzf1/woQJVtp7R+MSAHr16oUNGzYI/laE+yE1AOGxvPnmm9Dr9Zg8eTKKi4vx9ttv48knn8S+ffvMZbZt24bevXsjOTkZM2bMgF6vx5IlS9CjRw/88ssv6Ny5MwDgwIED2LNnDwYPHozGjRvj/Pnz+PDDD3H//fcjOzsbwcHBVs8eM2YMIiMjMX36dEma2dmzZyMgIACTJ09GRUUFAgICJNftyJEjSE9PR0xMDGbNmoXq6mq8+uqriIyM5H3Wtm3b8NVXXyEjIwMRERFo3rw58vPzcc8995gFxcjISPz4448YMWIESkpKMH78eAA15gEvvvgiHnvsMfMicezYMezbtw9PPPEEgJpjta+//hoZGRlITExEQUEBdu3ahZMnT+Kuu+7irRNjDA8//DC2b9+OESNGoGPHjvjpp5/w0ksv4dKlS/jggw+syu/atQvffPMNxowZg5CQEPznP//BgAEDcPHiRTRs2FDwdy4rK8P999+Ps2fPIiMjA3FxcVizZg2GDRuGoqIijBs3zqr8ihUrUFpaiueeew46nQ5vv/02Hn30Ufzxxx/w9/cXfM4ff/yBdevW4fHHH0dcXBzy8/Px0UcfoXv37sjOzkZsbKzgvQCwZ88eu9/queeew+XLl5GZmYkvv/xS9H5LBg0ahDZt2uDNN9/EDz/8gNdeew0NGjTARx99hB49euCtt97C8uXLMXnyZNx9993o1q0bgBqt/cMPP4xdu3Zh1KhRaNOmDY4fP44PPvgAv//+u+jRe2VlJdLS0lBRUYEXXngBRqMRly5dwvfff4+ioiKEhYXhyy+/xLPPPovOnTtj1KhRAGAW3uWOt4yMDISHh2PmzJk4ffo0PvzwQ1y4cMG8WXQWoTrec889aNKkCZYvX45//vOfVvcsX74c8fHxSElJ4f3OTp06oUWLFvjqq68wdOhQq89Wr16N+vXrIy0tDYBzY0mImzdvokePHsjNzcW4ceNgNBqxYsUKbN++3eG9csfnzp07sXr1arz44oswGAxYtGgR0tPTsX//fiQlJQGQ38ZSGThwIJo3b445c+Zg7969+M9//oPr16/jiy++MJd5/fXXMW3aNAwcOBDPPvssrl69ivnz56Nbt244cuQIwsPDAQBr1qzBrVu3MHr0aDRs2BD79+/H/Pnz8ddff2HNmjUApI3L5ORkfPDBB/jtt9/M7094GIwgNGTJkiUMAMvJyTFf6969O+vevbv57+3btzMArE2bNqyiosJ8fd68eQwAO378OGOMMZPJxBISElhaWhozmUzmcrdu3WJxcXGsV69eVtdsycrKYgDYF198YVe/e++9l1VVVTl8H66uLVq0sHqGnLr169ePBQcHs0uXLpmvnTlzhtWpU4fZDlEATK/Xs99++83q+ogRI1hMTAy7du2a1fXBgwezsLAwc9369+/P2rZtK/pOYWFhbOzYsaJlhg4dypo1a2b+e926dQwAe+2116zKPfbYY0yn07GzZ89avUNAQIDVtV9//ZUBYPPnzxd97ty5cxkAtmzZMvO1yspKlpKSwurVq8dKSkoYY4zl5OQwAKxhw4assLDQXHb9+vUMANuwYYPoc8rLy1l1dbXVtZycHGYwGNirr74qeu/t27eZTqdjkyZNsvts7Nixdm3KAYDNmDHD/PeMGTMYADZq1CjztaqqKta4cWOm0+nYm2++ab5+/fp1FhQUxIYOHWq+9uWXXzK9Xs9++eUXq+csXryYAWC7d+8WfIcjR44wAGzNmjWi71q3bl2rZ3LIHW/JycmssrLSfP3tt99mANj69evN12znCa6NlyxZYr7G/WZS6jh16lRmMBhYUVGR+dqVK1dYnTp1rNqBj6lTpzJ/f3+rvlVRUcHCw8PZM888Y74mZSxJ5b333mMA2Lp168zXysrKWOvWrRkAtn37dvN1V8cnAHbw4EHztQsXLrDAwED2z3/+03xNahtzc6Rl/fjg2u7hhx+2uj5mzBgGgP3666+MMcbOnz/P/Pz82Ouvv25V7vjx46xOnTpW1/nqOGfOHKbT6diFCxfM18TGJWOM7dmzhwFgq1evFn0Hwn2QmQHhsQwfPtzKRu6+++4DUKM1A4CjR4/izJkzeOKJJ1BQUIBr167h2rVruHnzJnr27Imff/7ZbLsaFBRk/p7bt2+joKAALVu2RHh4OO+R38iRI60cdxwxdOhQq2dIrVt1dTW2bNmCRx55xErb17JlS/Tu3Zv3Wd27d0diYqL5b8YY1q5di379+oExZn7WtWvXkJaWhuLiYvM7hoeH46+//hI9ag8PD8e+fftw+fJlye+/ceNG+Pn54cUXX7S6PmnSJDDG8OOPP1pdT01NtTqCb9++PUJDQ81tK/Yco9GIIUOGmK/5+/vjxRdfxI0bN7Bz506r8oMGDUL9+vXNf9v2ISEMBgP0+prpsbq6GgUFBahXrx5atWrl8Ii4sLAQjDGr57rCs88+a/5/Pz8/dOrUCYwxjBgxwnw9PDwcrVq1snqvNWvWoE2bNmjdurVVn+jRowcAiGr0wsLCANQcm9+6dUt2neWOt1GjRllpykePHo06deqoGhLp6aefRkVFBb7++mvztdWrV6OqqgpPPfWU6L2DBg3C7du38c0335ivbd68GUVFRRg0aJD5mjNjSYhNmzahUaNGePjhh83XAgMDMXLkSIf3yh2fKSkpSE5ONv/dtGlT9O/fHz/99BOqq6sByG9jqYwdO9bq7xdeeMH8DgDwzTffwGQyYeDAgVb92mg0IiEhwapfW9bx5s2buHbtGv7xj3+AMYYjR45IrhM3lj0pbBhhDQmzhMfStGlTq7+5CYWzpTtz5gyAGkEyMjLS6t+nn36KiooKs31fWVkZpk+fbrYXi4iIQGRkJIqKisxlLImLi5NVV9vyUut25coVlJWVoWXLlnbfyXeN71lXr15FUVERPv74Y7tnDR8+HECNnSkATJkyBfXq1UPnzp2RkJCAsWPH2tnwvv322zhx4gSaNGmCzp07Y+bMmQ6FvwsXLiA2NhYhISFW19u0aWP+3BLbtgVq2tfWTpLvOQkJCWZBU+5zbPuQECaTCR988AESEhKs+suxY8d4+wsfTCH7Ott3CAsLQ2BgICIiIuyuW77XmTNn8Ntvv9n1iTvuuAPA332Cj7i4OEycOBGffvopIiIikJaWhoULF0p+d7njLSEhwervevXqISYmxmX7SzFat26Nu+++2yoO8PLly3HPPfcIjj2ODh06oHXr1li9erX52urVqxEREWHeLADOjSUhLly4gPj4eDuzC0d15e6VMz5t2wMA7rjjDty6dQtXr14FIL+NpWL77Pj4eOj1enNfOHPmDBhjSEhIsOvbJ0+etOrXFy9exLBhw9CgQQPUq1cPkZGR6N69OwDIqiM3limGsedCNrOExyKkGeUmFk7r+s4776Bjx468ZevVqwegZne/ZMkSjB8/HikpKQgLC4NOp8PgwYN5Iw9Y7uilYFteat2ciUQg9KynnnrKzoaPo3379gBqFq/Tp0/j+++/x6ZNm7B27VosWrQI06dPN4emGThwIO677z58++232Lx5M9555x289dZb+OabbwS1xXJx1LZK4exz3njjDUybNg3PPPMMZs+ejQYNGkCv12P8+PEOI1U0aNAAOp3OocAsFb53kPJeJpMJ7dq1w/vvv89btkmTJqLPfe+99zBs2DCsX78emzdvxosvvmi2YxRzigTkjzd38fTTT2PcuHH466+/UFFRgb1792LBggWS7h00aBBef/11XLt2DSEhIfjuu+8wZMgQq4gUWowld6FVG9sKkCaTCTqdDj/++CPvOODm/OrqavTq1QuFhYWYMmUKWrdujbp16+LSpUsYNmyYrDpyY9l2A0l4DiTMEl4Ld0wdGhqK1NRU0bJff/01hg4dauXxXV5ejqKiIrfWLSoqCoGBgbxRGviu8REZGYmQkBBUV1c7/B0AoG7duhg0aBAGDRqEyspKPProo3j99dcxdepUc1ijmJgYjBkzBmPGjMGVK1dw11134fXXXxdcgJs1a4YtW7agtLTUSvtz6tQp8+dK0KxZMxw7dgwmk8lKO6v0c77++ms88MAD+O9//2t1vaioyOGCVqdOHcTHxyMnJ8fuMy01O/Hx8fj111/Rs2dPp5/brl07tGvXDq+88gr27NmDrl27YvHixXjttdcACL+P3PF25swZPPDAA+a/b9y4gdzcXPTp08epelsi9u6DBw/GxIkTsXLlSpSVlcHf39/KTECMQYMGYdasWVi7di2io6NRUlJilSCDQ+5YEqJZs2bIzs4GY8zqnaTME3LHJ3eyZMnvv/+O4OBgs2OqWnPqmTNnrE6fzp49C5PJZI6eEh8fD8YY4uLizKcMfBw/fhy///47Pv/8czz99NPm65mZmXZlHY0PbixzmmzC8yAzA8JrSU5ORnx8PN59913cuHHD7nPuOAyo0WTZauPmz59vtv9yV938/PyQmpqKdevWWdnVnT171s6OTQg/Pz8MGDAAa9euxYkTJwSfBdSkZ7QkICAAiYmJYIzh9u3bqK6utjt+i4qKQmxsrGg4pz59+qC6utpOq/XBBx9Ap9MppoXq06cP8vLyrI53q6qqMH/+fNSrV898hOgqfP1lzZo1uHTpkqT7U1JSeFO01q1bFwBU20RZMnDgQFy6dAmffPKJ3WdlZWWiUTpKSkpQVVVlda1du3bQ6/VW/aBu3bq87yJ3vH388ce4ffu2+e8PP/wQVVVVivQboToCNZq23r17Y9myZVi+fDnS09Mla9/atGmDdu3aYfXq1Vi9ejViYmLMkSQASB5L165dw6lTpxzaJqelpeHSpUv47rvvzNfKy8t529cWueMzKyvLyu71zz//xPr16/Hggw+ataFqzalcmCzL7wRgruOjjz4KPz8/zJo1y+75jDHzHMfV07IMYwzz5s2ze6ajcXno0CGEhYWhbdu2TrwRoQWkmSW8Fr1ej08//RS9e/dG27ZtMXz4cDRq1AiXLl3C9u3bERoaig0bNgAAHnroIXz55ZcICwtDYmIisrKysGXLFtEwUFrVbebMmdi8eTO6du2K0aNHmxedpKQkHD16VNLz3nzzTWzfvh1dunTByJEjkZiYiMLCQhw+fBhbtmxBYWEhAODBBx+E0WhE165dER0djZMnT2LBggXo27cvQkJCUFRUhMaNG+Oxxx5Dhw4dUK9ePWzZsgUHDhzgjWPK0a9fPzzwwAP4f//v/+H8+fPo0KEDNm/ejPXr12P8+PFWzl6uMGrUKHz00UcYNmwYDh06hObNm+Prr7/G7t27MXfuXDubQGd56KGH8Oqrr2L48OH4xz/+gePHj2P58uVo0aKFpPv79++PL7/8Er///ruV9ohzqnnxxReRlpYGPz8/Xm2eEvzf//0fvvrqKzz//PPYvn07unbtiurqapw6dQpfffUVfvrpJ3Tq1In33m3btiEjIwOPP/447rjjDlRVVeHLL780b5ws32fLli14//33ERsbi7i4OHTp0kX2eKusrETPnj0xcOBAnD59GosWLcK9995r5ezkLEJ15Hj66afx2GOPAagJsSeHQYMGYfr06QgMDMSIESOsTgtKS0sljaUFCxZg1qxZ2L59u2gc3eeeew4LFizAkCFDMG7cOMTExGD58uXm0xQx7aLc8ZmUlIS0tDSr0FwArLJkqTWn5uTk4OGHH0Z6ejqysrKwbNkyPPHEE+jQoQOAGs3sa6+9hqlTp+L8+fN45JFHEBISgpycHHz77bcYNWoUJk+ejNatWyM+Ph6TJ0/GpUuXEBoairVr1/Ka/zgal5mZmejXrx/ZzHoymsVNIAgmLzSXbVggvlA8jNWEEXr00UdZw4YNmcFgYM2aNWMDBw5kW7duNZe5fv06Gz58OIuIiGD16tVjaWlp7NSpU6xZs2ZWYXu4+h04cEDS+wjVVU7dGGNs69at7M4772QBAQEsPj6effrpp2zSpEksMDDQqhwAwVA/+fn5bOzYsaxJkybM39+fGY1G1rNnT/bxxx+by3z00UesW7du5vrEx8ezl156iRUXFzPGasILvfTSS6xDhw4sJCSE1a1bl3Xo0IEtWrTI6lm2oX8YY6y0tJRNmDCBxcbGMn9/f5aQkMDeeecdq9BkYu9g2xZC5Ofnm9syICCAtWvXzq5PcH3lnXfesbsfNiGw+CgvL2eTJk1iMTExLCgoiHXt2pVlZWXZ9VUhKioqWEREBJs9e7bV9aqqKvbCCy+wyMhIptPprMIB2daLC1V09epVq+8YOnQoq1u3rt0zu3fvbhd2rbKykr311lusbdu2zGAwsPr167Pk5GQ2a9Ysc5vz8ccff7BnnnmGxcfHs8DAQNagQQP2wAMPsC1btliVO3XqFOvWrRsLCgpiAMztJ3e87dy5k40aNYrVr1+f1atXjz355JOsoKDA7v2cCc0lVEeOiooKVr9+fRYWFsbKysoEfxM+zpw5Yw5ltWvXLrvvlTKWuDo7Cl3FWE279O3blwUFBbHIyEg2adIktnbtWgaA7d2711xOifG5bNkylpCQwAwGA7vzzjvt6ie1jeWG5srOzmaPPfYYCwkJYfXr12cZGRm87bJ27Vp27733srp167K6deuy1q1bs7Fjx7LTp0+by2RnZ7PU1FRWr149FhERwUaOHGkOA2jZb8TG5cmTJxkAu75PeBY6xiilBUF4Io888gh+++03Xvs1wvOZPXs2lixZgjNnzsgK80ZoS1VVFWJjY9GvXz87G2lvYO7cuZgwYQL++usvNGrUyN3V8TnGjx+Pn3/+GYcOHSLNrAdDNrME4QFYplcEapwgNm7c6DB9J+G5TJgwATdu3MCqVavcXRVChHXr1uHq1atWTkKeiu08UV5ejo8++ggJCQkkyKpAQUEBPv30U7z22mskyHo4pJklCA8gJiYGw4YNQ4sWLXDhwgV8+OGHqKiowJEjR3hjPhIE4Rr79u3DsWPHMHv2bERERLgU6F8revfujaZNm6Jjx44oLi7GsmXL8Ntvv2H58uXmdNQEURshBzCC8ADS09OxcuVK5OXlwWAwICUlBW+88QYJsgShEh9++CGWLVuGjh07YunSpe6ujiTS0tLw6aefYvny5aiurkZiYiJWrVolOZwYQfgqpJklCIIgCIIgvBaymSUIgiAIgiC8FhJmCYIgCIIgCK+lVtrMmkwmXL58GSEhIeShSBAEQRAE4YEwxlBaWorY2FirpCS21Eph9vLly2jSpIm7q0EQBEEQBEE44M8//0Tjxo0FP6+VwiyX8vLPP/9EaGiom2tDEARBEARB2FJSUoImTZo4TFVeK4VZzrQgNDSUhFkeqk0M+3MKcaW0HFEhgegc1wB+ejLHIAiCIAhCexyZhNZKYZYQZtOJXMzakI3c4nLztZiwQMzol4j0pBg31owgCIIgCMIeimZAmNl0Ihejlx22EmQBIK+4HKOXHcamE7luqhlBEARBEAQ/JMwSAGpMC2ZtyAZfBg3u2qwN2ag2UY4NgiAIgiA8BxJmCQDA/pxCO42sJQxAbnE59ucUalcpglCZahND1rkCrD96CVnnCmizRhAE4YWQzSwBALhSKizIOlOOIDwdsg8nCILwDUgzSwAAokICFS1HEJ4M2YcTBEH4DiTMEgCAznENEBMWCKHgFzrUaK06xzXQsloEoThkH04QBOFbkDBLAAD89DrM6JcIAHYCLff3jH6JFG+W8HrIPpwgCMK3IGGWMJOeFIMPn7oLUaEGq+vGsEB8+NRdZEdI+ARkH04QBOFbkAMYYUV6Ugy6toxAu5mbAQBLh9+N+xIiSSNL+AxkH04QBOFbkGaWsMNScKVUtoSvQfbhBEEQvgUJswRB1CrIPpwgCMK3IGGWIIhaB9mHEwRB+A5kM0sQRK2E7MMJgiB8A9LMEgRRayH7cIIgCO+HhFmCIAiCIAjCayEzA4IgCA+k2sSwP6cQV0rLERUSSJpjgiAIAUiYJQiC8DA2ncjFrA3ZVpnKYsICMaNfIjmnEQRB2EBmBgRBEB7EphO5GL3ssF3K3bzicoxedhibTuS6qWYEQRCeCQmzKlNtYsg6V4D1Ry8h61wBqk3M3VUiCMJDqTYxzNqQDb5Zgrs2a0M2zSMEQRAWkJmBitBRIUEQctifU2inkbWEAcgtLsf+nEKkxDfUrmIEQRAeDAmzKsEdFdrqT7ijQgrMThCeiTsdr66UCguyzpRTEnJIIzwZ6p+1GxJmVcDRUaEONUeFvRKNNNgIwoNw92lKVEigouWUwt2/C6EOviIAUv8kSJhVAToqJAjvY9H2s1iw/ZzddS1PUzrHNUBMWCDyist5N8M61KTc7RzXQNV6WEKnTL6JrwiAvt4/fWXDoTbkAKYCnnxUSBAEP3yCLKCt45WfXocZ/RIB1AiulnB/z+iXqNliRg5pvokvRMyoNjHsPnMN/1p73Gf756YTubj3rW0Y8slejFt1FEM+2Yt739rmFe2jNSTMqoCnHhUSBGFNZna+pHKWpylqk54Ugw+fugtRoQar68awQM21THJOmQjvwBc2KJyQ9+R/96Go7LZgOW/un76w4dASEmZVgDsqFNKd6FBznKPlUSFBENZUmxje2HhS1j1ipylKhuFLT4rBlondzX8vHX43dk3poflxKZ0y+R7evkEREvLEyMzOU7FGyuMLGw6tIZtZFeCOCkcvOwwdYNUh3XFUqBRku0P4EvtzCpFfUiHrHqHTFDXsDy3HlrvGGp0y+R7evEERE/LE+Gz3eXSOa+ARtrNS1lHyu5GPJprZhQsXonnz5ggMDESXLl2wf/9+0fJr1qxB69atERgYiHbt2mHjxo1Wnw8bNgw6nc7qX3p6upqvIBtPOipUArLdIXwNOYu12GmKLx8H0imT7+HNGxRHQp4QXAQhd2sypa6jrm44uFOibw//hf/+8ge+PeL7SZtUF2ZXr16NiRMnYsaMGTh8+DA6dOiAtLQ0XLlyhbf8nj17MGTIEIwYMQJHjhzBI488gkceeQQnTpywKpeeno7c3Fzzv5UrV6r9KrLxlKNCV/HlxZqovchdrPlOU3z9ONDTHNII1/HmDYqz2mJPMJ2Qs466suGwFJgnfPUrZv9wEhNW+74CSnVh9v3338fIkSMxfPhwJCYmYvHixQgODsZnn33GW37evHlIT0/HSy+9hDZt2mD27Nm46667sGDBAqtyBoMBRqPR/K9+/fpqv4pTeMJRoSv4+mJN1F46xzVAtM3JCR8xIqcp3m5/KAVfO2Wq7XjzBsVVbbG7TCfkrqPObjgc2RPn+rACSlVhtrKyEocOHUJqaurfD9TrkZqaiqysLN57srKyrMoDQFpaml35HTt2ICoqCq1atcLo0aNRUFCg/AsQOHThus8v1kTtxE+vw7/7tBEtk/FAvOhpijfbH8rBV06ZiBq8dYPiSMhzhLtMJ+Ruep3ZcMixJ/ZFBZSqwuy1a9dQXV2N6Ohoq+vR0dHIy+P3LszLy3NYPj09HV988QW2bt2Kt956Czt37kTv3r1RXV3N+50VFRUoKSmx+kdI42qpNAcZb1+sidpJr8Ro0c/HPNBSVEPlzfaHcvH2UybCGndvUJyJ/iEm5InhbtMJZza9cjccUu2JfVUB5ZXRDAYPHmz+/3bt2qF9+/aIj4/Hjh070LNnT7vyc+bMwaxZs7Ssos8QGeL4GBYAzl+7pXJNCEIbFj91F55fdlhSWU/M2EUQUnHXBsWV6B+ckDfju9+sopGEB/nzxpz1BNMJZze96Ukx6NoyAu1mbgZQs+G4LyGS9z3kKpR8TQGlqmY2IiICfn5+yM+3Dkyen58Po9HIe4/RaJRVHgBatGiBiIgInD17lvfzqVOnori42Pzvzz//lPkmnoOSsSylkNysvqRjnblbfvdJOxyi9tGpuXT7e2+2PyQId6CEQzGfVvmXKQ/wlvUE0wlXnO6kbjjknv74mgJKVWE2ICAAycnJ2Lp1q/mayWTC1q1bkZKSwntPSkqKVXkAyMzMFCwPAH/99RcKCgoQE8PfWQ0GA0JDQ63+eSPuCI/FLda11Q6H8Gykbu7U3ARqZX+o9UaWIJRGSYdiKUKep9h2a7Hp5QRmqfiaAkp1M4OJEydi6NCh6NSpEzp37oy5c+fi5s2bGD58OADg6aefRqNGjTBnzhwAwLhx49C9e3e899576Nu3L1atWoWDBw/i448/BgDcuHEDs2bNwoABA2A0GnHu3Dm8/PLLaNmyJdLS0tR+HbfB7WZthzi3m1Vz55meFIMJqQn4YMsZwTIUxJnQGqlHlWLlut0RqUhd5BwHOkNmdj7e/PGUokkZCEJrtE4G4Em23ULmEUaFxrFlsiap29xZG7LRK9HoMb+RK6gemmvQoEF49913MX36dHTs2BFHjx7Fpk2bzE5eFy9eRG7u37uDf/zjH1ixYgU+/vhjdOjQAV9//TXWrVuHpKQkAICfnx+OHTuGhx9+GHfccQdGjBiB5ORk/PLLLzAYpNl3ehtahcfiND8/HLPfrTWPqCvpO3zNDofwTKQeVToql5ltbdLkCmraH45fdZTiPBNeT22J/iGE2k53nMAsJeSgrzmCaeIAlpGRgYyMDN7PduzYYXft8ccfx+OPP85bPigoCD/99JOS1fN4tNjN8mmvgBqNUP+OjWqV1zbh2Tja3HHZfnq0jnZYbs7GU2pWVTEcvauvaFcI30bq+iDXnnN/TqHkyDvuRm2nO9tTIkf4ysbBK6MZ1DbU3s0KmTAANRohQx09eiUaER1qEMxlT17bhFZI3dx9mXXeYbm8EvExsz+nEMVltwVzqLsbMu8hvAlH0T845m75Ha2M9SRrLIctOaBMBX0EOfOUryigVDczIFxHTa2olEDLszZkA4BggHny2ia0ROqm7UKh6966w5YccMrZcn9OoaaOWr6iXSF8GyUdipU0EaqNuDv2rtKQMOsFdI5rAGOosKDqSqeUY8IgFGDeE0KfELUHqZu2Zg2Cnfr+baeu8F6XY6PqrBDsLL6iXSF8H86hWAxH9pzVJoY3Np50+CyK+MGPLyqgSJj1AjKz81BexZ/dzNVO6aoJg6eEPiFqD1JjNv5fSnOH5fg2ie/+9DtveUfOlkKaIiUctZyJT6k1FDqMkIqrDsX7cwoFTd4sOXThuqx61RZ8UQFFNrMejpg9KwCEB/tjzqPtnO6UrpoweKIdIeHbWIag0cHaOcpycxdQRy9YDv/7e0p6K0z46ler61dEHEmEbFTFNEVKOWqJvauSY7DaxLA/pxBXSssl2wq7ktGJ8E3E+pGr645UJYy3OIVpidJhAz0FEmY9GCn2rJxzlrPIScdZIaAdJrwHqYKKMwKNlkiN2ShUjuPNTc5FM7BdTB1pilx11Jo7uCPe2HhSlfiUljgTk9edMbAJz8TR5qZzXAOXHIqlCsNS07HXJjxtLlcKEmY9GEf2rACQV1LhkiezmJaLw5fsamozfAuMMdSAIZ2bonlEXbPQmpmd5xVaNqmJCtKTYlBWWW2ngQWAKxKOKvmwXUzVjjjSKzEaPVpHqZaUAXAslM4d3NHuHqlh0ih0WO1BqB/lFpfj+WWHseiJO9GnfSz+3acNxq06ane/lBMHR8IwR3Iz6ampCe+GhFkPRqsA02Laq7mDO3qUAEM4h6CgUlJhldktPNgfRbdu293vqVo2KTEbq01MUAPLJ4RFhRhETQ3Cg/xhYgzVJqbYsakU1IxPKUUo5YvJq3VGJ8KzkXKamLHyCBZAJ+pQ7Gjj7KfXCQrDtuU8GU8P/edNkDDrwai9QHJHyXnFZSi8WYnxqQmY+s0JqzJCEw7hPUhZYDj4BFnAfVo2IXMH7vpf1x2H35LqLMIxOe0OvPz1ccHPi8pu48lP9yl6bOpupAilfDF5a3tGJ8IaKaeJJgaMWXEY83g0/VPSW8EYFoSwoACrzSIfvrA2WcbH9cTTL2+ChFkPRkqAaT4tkRSEMn4RvoeUBUYKWmvZhOzuHu4Qg+9+zRXMVmeLXEGqR+soSeVstdWuHJu6ghL2zc4Km+7ODOjptt21DTn9iE/T/9am0+b/91XhzlHUE087/fIWSJj1YKTYs/JpiRzhKEKCHDzpmIQWNn6U1oppoWUTs7v76Occ3nu4bHW2Y0AtQcpWW+3KsakjhNJ1uhJFwHK8XHPS61uOA6kSWNb5/LVbWLn/opXG2FcFIG9BzlhzlH3PF4U7LaKe1FZImPVwHHljc0gd+HKOnLnyYnjKMQmFBhJGaWFO7QD9cvuoJXwLgVRnEWew1FZ3aBJm97lSjlq26Tozs/NhqKN3OooA33jR62qOgPnQAYgODbQTQKSGSVNiYZZymuSLApA3wW1ulDoJ0kK401IJonbUk9oMJU3wAtKTYrBlYnfRMo4CunPIPXKWE3RaieDwzsBp8Wzfy1318QQsA9ibTAzGUOHkAVLRKkC/s2YRQlmDOGcRqTgT7D8zO4/3uisLo1i6znGrjuJf3xwXdNgChOcCofEiJsgCwNQ+rXk/5zbcUaHWYZCUDMwuVGdbpM6DhDpwmxulcJQJzFU2ncjFvW9tw5BP9srO2OdMkhC1bcy5Ov1wrPateaSZ9RKkLIhiuzpu9/mjTMFOTtBpdxyTUGgge/g0WMEBfubfw9klnkH9MG3VJobdZ6+59B18C4EcZ5EjF4tkP/Oz3efRoUm47PuEkJKuU8hZDxBP7uBI622roTU6iDMLSA+T5gxyNfWk3XIv6UkxWPTEnchYeURwgwQAYUH+KC4T7sOWqGHa5Ep8ZGdPAtW0Ma/tfjAkzPogtgPflU4uN+i01gsJhQayRmiCvlXpesKLZ7o2d1nLJnakp9Rk7KoZxLUb8s0RhEJXOYvcCAxC8CV3kOJtzmEplN6qrBK9T63QYc5q6imCwt9o7U/Qp30s3qsy8cZ25pAqyALKmza5ogRxVgiuNjGYTExUiHfWxlxJPxhvhYRZH8Ry4LvayZ0NOq3VQkKhgf5GigbLlcmOL9OcnEVSTJsBwOXJWClno4h68rMGCYWuchal+quzyR0s+f7YZUSFBCKpUagidZKL2pEWfB1nsrq5AjcnVClg5qFWSDtnlSDOCsFSNurO2pi74mPgS5Aw60PYDnwlOrmzu3e+hUQN7YC7QwO5C77fUqkQXLYILShSM4r56XWi2oznlx1GeLC/IpOxEmYQdzYNV6AmruFqfxVqM7nfa+l4FhrknuVCbp09Paavlmw8losxKw7bXRfL6uYKSh51qxnSzlkliDNCsFSFUliQP4Z3bS6anp4vepBa8763QcKsj8A38N3RyWu8ng0wMYb1Ry+pniJVSixevQ64frPS6Wd4GkKalj5JwpOgq9guKFIzisWEBWJa3zaY/cNJUWclMftPPvg875XKVucJttVKRGDgEwKkhNIKC/JHEc8xaEnZ32YGQjF91UDKGOdQO6avN7Hx2GVkrDzC+5lYVjdnUfqoW4mQdkI4qwSRKwTLSlhTdhsfbDmDVQf+xIx+ieiVaLRLDMMXPaiiyiSpTr4OCbM+At/A1/ponXMuKq8y4clP95mvq5ki1TI0kBAmBoxdcRgf6r0/XI+YhvO/u8+r8kxbIbHaxDDzO2kTdF5xOcas4F9QXYHvBNPSyctWc632EblQ6CpnkZquU4jxqXegV6IRWecK7E5CxGJXMwCV1Y4XR6GYvmogJd42h5oCkBhyT53UtmHddCLX4biTYhojtZ5KH3Ur6UDIh7PxkeUKwc4olCxPq8Q2+dz6OT41Qdb3+yokzPoAlgPfMkXt4YvSw2opQVCAH25VVtsNQLVTpKYnxWChBO9Zb49q4MheCxCPFeostpEAFmw7I1lo09KOizuC4wumHx0q3w5WKpahq5wVPi3hxnClCxqX4rJK3PvWNsGTEL7Y1VzfkeIsyKDteJISbzvjgXhM6NVK8/Et17Nd7ZjY3DzhKnLqqfQpoJLCvZBA7kx8ZLlCsDMKJamnVdz6uXL/RRhDA5Ff4vjkAvCsREdKQsKsD2BplyjHXsmVME18OOMxr1S0gfp1DaJCnC9ENZDriS4XKdrFTSdyrcwIPAnbxAKWXFEhYQKHlNBVfPAtsnzmOM7wGY+W3vYkxDKUFiC/72g9ntKTYlBWWS3oIb9w+zkkNQrTVCsr17PdlXBQUlFCsJRbT091sM3MzsebP57iFch7JRoxPvUOfLY7xyq6gJh2X64QrLavRo12vQITUu/A3C2/S1rTPSXRkdKQMOvhcAuepd0MH3LtlTxtL+bqZOgNUQ1cPVqUWvdnujbH2sN/obhMPJQSH2LaRaU0Pu5ALQ3xlPRWGNUtXlLoKkv4Np6OjhWlIqSdVyPuspbjqdrE8OYmcRtPT45xLeVk5V9rjyMk0B/3tGjo9DvIaRMjz+bVGY99OUKb0koUMcavOirqdGo73qRo94VOCfiEYDn23q7QPCJYUqZQW3wpYx5lAPNgLLOTvPT1McFyztgrGcMCFfdkdQVXd7DujGogJROMK5lmOKTWvVF4EHa+9ADvZ8ZQ4e+YO7ijaHIBLRwKve3E66l7mskWOoSyWSkhyALiGlalMyppGSVETipQW5zJ1iSlPlI926WUB2qcgJ78dJ/sucESOW3Cl9VN7nsBfwttjkbCBwM72GWKUxO5TqcLt58TzOZniW1WzqXD78auKT3sBELLjGhqTm1RIYGSMoXa4ksZ80iY9VCkpm8EalLOyhEypqS3wq4pPWRlRVILpVKkct7faj/HFilCqlLpdqUuGLN/OIm0uT/bXV86/G5kTuwmeF9YkL9oGkRntXDcbz++p2NHBcv59LOhnVS1dVUCuYKs2jEhn05pJqmcqxpVrVIbW+Ls6YsSG0kl6iPnN3clFbeUeUKvAxY9cRfvGuDM7yxVaEtLMsoWuGxRO2WrVMFOapIQoXTPSmA7Dp3R5qudMlgrSJj1QOQueHJSzgI1QeE9wehbyTA6nPe32s+xRIqQKuVoUc7kOa1voqR+wWcj6sisYdiSA6InAM5q4RiAPklG2UHU74lvKNim3ora2u0eraMklVNCo6p1+CtnTl+U2kgqUR85v7krGjMpguWCIXeiT3v+Y2VnT7mkCm2u9Bmpp5XOooZgV21iCAsKwMRedyj2nYDy65oUjbQnQ8KshuzPKZR0zCV3wZObflNuilpX4YZZeJC/1XVjWKCitjpCmmalnwNIiywwa0M29v5RIPvITohNJ3Ix+wdpNqtqaP6kanz4/v7v7vNYsP2s7Gd6wumBkqhtY5rcrL5oGwlpVOUuhUrF9JWD3NMXJTeSQvWR81tLPVmxrKOzgpWYYDlvcEf0aR8reK/c97J9rquaVyHknFa6itxxKrS2WwrfU9Yel/Rd4cE166SjfqL0uvbZ7vMun1a4E3IA0xCpXoRyB9Jbm07LKu9silpn4QzjLb2n1Y4jyKHEc5zJtsUtRFnnCiQ9w1GbuzP39t5zBSitqEJUSCCm9U3E2BX88T51AN57vIOVt7kcOUHpWK2eiNo2plLiyE7r24Z3PMhxzHHHJkNK7F05SWNcjXAi17NdTrxcS5zdAPFFrACE284yZJPQOJeiDVRjTtc6ZasrmfJcSdHNrVeZ2XmizlxqrJ9KO4dqDQmzKpOZnc97XcyLUIsFTyssB52ltze3q+cL6q4kYt8pJbqA69m2pE1lYm0uJ0mBGjzz+UHz/8eEBWJUtzisO3rZbqI1O5AJhE4SQ+lYrVoiNeIIoI13s6OYrLN/OAm9XmcVSmzu4I54Y+NJlzKOaYGYEG2rLdYiwokcz3ax8mK4sh7ImU9tBTK+ce6upBRaZbMUSpbAh9ja7myKbm4NEtqI2JZTEm8PX0nCrIpUmxje2HiS9zOxMDlqLHhhQf7mWHqcAb2UxdcZLMO9CA06sfh/WkyUUgKCK5FtK6VFBNYeviQ704wlcpIUqE1ecTk+/jkH7w/sYBfv0xVtnbOxWqUSHWJAvkzbcinw9WMxHGnzGGrMcfjSycpBLCYrt5G2jGbSKzEaPVpHCS6e3oBt/1MywonYxtdW8ODTmtne/3JaK0xaI27zKUewchYxgcx2nGt1msaHliHgpNigOlrbAdejk4jVQc3EB54aM9gRmtjMLly4EM2bN0dgYCC6dOmC/fv3i5Zfs2YNWrdujcDAQLRr1w4bN260+pwxhunTpyMmJgZBQUFITU3FmTOeF8jd2VAyaoTzsAwKfd9b21UzoHfkMc8xftVRVZwypKCU45ZeJ9w+nF3ZPfENBdtSypGdpyUpYP/7N1PBeLNCYW2UQgfg333VcSTj68eOELJnNIYFYvFTd+GXKfxh1eQgFpOV679zNlp/ruSCKNU/QE1csf20REo0BDHPdr77pc69cp175IQgkyKQWZqwuTNblFYh4KTagjta29Vm2JIDVn1x47HLikV40DLcnpKorpldvXo1Jk6ciMWLF6NLly6YO3cu0tLScPr0aURF2Xve7tmzB0OGDMGcOXPw0EMPYcWKFXjkkUdw+PBhJCUlAQDefvtt/Oc//8Hnn3+OuLg4TJs2DWlpacjOzkZgoOc0hCvHXM4cR0nFVa2PGFK1CHICciuJ1IDgIYH+krNtObIrk3sUaVtXT6SYpw/tzymUHVkDqLHhVnORdBQ/1xUciWlCwoSYNk9O8gUhpGyk1dT2e0KWIVfshzlczdgldL8U+X586h2yfjO+0yY+hzlOq3ettELzPuLsHKFV8gHLeUJMG+9J2svc4nKMWXHE5e/R4iRATVQXZt9//32MHDkSw4cPBwAsXrwYP/zwAz777DP861//sis/b948pKen46WXXgIAzJ49G5mZmViwYAEWL14Mxhjmzp2LV155Bf379wcAfPHFF4iOjsa6deswePBgtV9JMq4eczmym/FF5Njt2E42SY1C7crYHsco7bj1TNfm+OF4rkMhVcpRJF/dtbATUwqxdLJiHLpwHfe3khZSyhncGRGB7934bGx94ajQGf8AR0gZ42JItR/mq5ecTFhCdXfFaal5RLDkskJCM987OztObXEmo6Gzz3bWeU4u3ObTkRmat2ovhVArfKWWqCrMVlZW4tChQ5g6dar5ml6vR2pqKrKysnjvycrKwsSJE62upaWlYd26dQCAnJwc5OXlITU11fx5WFgYunTpgqysLI8SZjvHNUB0SACKrpcKltHrgKKCYphigng/11VWwVDlvuMM061buF3ph4MXruPS9VsO62K6VbNAc+VMt27BVFXTzUwy3uXq1euCvwlQExPvjY2nrDQHUSEGu+9/7pNd5v83hgYirW20pDr4VZRJKtdAb8Lm0Xej8+tbAQAf/V8yusZHwE+vM/8WHJZt2Sk6ELryMphEvvvq1euKtr1t23gK164WwdSkHgDrPrL/5CVcK61wub5y35uVSWt7KVi+G2Dfbzm92ZbD5/HQ/8IlyRkntnDjLcqfyfoOub+R5bgGaoSAd9cfFbxfB2DON0fQs3mo5MVS6hgXqhPHgy3CUH5/M0z5xj400vWCCoxfmoV5gzvaCaX7cwpRWFAMsUCGhQUV2H/yEto1CrWb86TcL0aUP7ObQ/ioNjHM+eYIAlQe15a/L1/bGEMD8e8+rdG1ZYQiY9a2LR9sEYbFj7XBaz+cxBUV7N8B4NDpyyirqMb4r44CgFXbWfaVHq2j0SQIgvXQAQgL9kfx/+xmpQrfzq6XrsK1Xa8WYZL6HADogoKg03mO4KtjjKmmtb98+TIaNWqEPXv2ICUlxXz95Zdfxs6dO7Fv3z67ewICAvD5559jyJAh5muLFi3CrFmzkJ+fjz179qBr1664fPkyYmL+3k0PHDgQOp0Oq1evtvvOiooKVFT83SlKSkrQpEkTFBcXIzRU3k5fLt/tPYuEYf1UfQZBEARBEIRWtDp8CPpg6acHzlJSUoKwsDCH8lqtSJowZ84chIWFmf81adJEs2c/JBKcmiAIgiAIwttQwrZfSVQ1M4iIiICfnx/y861tqfLz82E08tsZGY1G0fLcf/Pz8600s/n5+ejYsSPvd06dOtXKdIHTzGqBLigIRWszMWyJeAQHMQ69korgAP6mqjYxswnA/1t3AgCw///1xMPzdyO/RLqxPBewfsvE7th2Kh8vOoj1aVmeOza8VVmF5Ne2mOsMwPy35btsyc7nPfLzZj76v2R0albf7n1tsf1ddk15APe+tV3288xhnIL9rULARIcEIv9/9pKW/ca2bYT6kyV8R4lK4UwYKr4+5whn3luM749dluyJ7kqILWfqavmuHO881h7vbf5dtA2jQgzm41Ipz+V7jjtYOrwzuvzPWcUT6iQ058khqp4BV2RmdFQbZ97ro/9LRreESMnj5Z3H2psVP1LbkpsB+ExEgL/Xxqul5WhY14B/fXNMcfOEQ6+k4pffr5nNEsTKOTOuuPuqTQyp7++UtaZzWI4TKc+UyuEgYTNAd6CqMBsQEIDk5GRs3boVjzzyCADAZDJh69atyMjI4L0nJSUFW7duxfjx483XMjMzzWYKcXFxMBqN2Lp1q1l4LSkpwb59+zB69Gje7zQYDDAYtE3hyqHT6XDltg4VdZx/vj44GHqBgaAHkNK2Lm5VVmHy9zUhnPzr1sXUR+/E6GWHATi21zEHrH/0TuiDgzEzM0dSfS/cYjiYX2521NLXqTLfxx0/2H7P1vOlmLPjgku/hydy9bYe+uBgh+9l+7vogoKc+i04Z4ReiUY7BxnOycyy39i2jVB/siStUwuk3hWHvecKMGbFYd4IBs7CCeH5twHIeH/bPucIZ95bjMjI+pLbS+67WbL1fCn6d2wk6x7Ld+VIvas5+naOx/6cQvx4IhdfZF2wu++vMoDJ+I34nsOhpnOOLVdu68zjSaxOWiE050nlsbsa4Zez11Dhpujvj93VCJknr9iN863nS9GjdZSs9+LmQ6njJTKyvlNtqQMwMzMHqXfF2W1wubWRY3K1XvGELCwwCK9udbxessAg6AP9RcvwvTc3HvWArDUdsIhO0KYR9AKbf1fGjSfZywIamBlMnDgRn3zyCT7//HOcPHkSo0ePxs2bN83RDZ5++mkrB7Fx48Zh06ZNeO+993Dq1CnMnDkTBw8eNAu/Op0O48ePx2uvvYbvvvsOx48fx9NPP43Y2FizwOxpuMPzkfPiNYZZPzs82N+c+5nDMsezXA96uV7T41cd9agsQ0uH3413Hmvv8vdEhjg3IciNwRke5I8JqQnmuKx+eh1S4huif8dGSIlvKJrtjGN/TqHk5/rpdeiaEIFX+7eVVU8xxtzfAoF1/Jy+351hcRzFLgVqnDpdZfyqo7LjLQu1qZ9eh85xDfDTb3m8n1vepURMWK2WOF/zKP/68CW3zY3hQXWw9vAl3g3ruFVHcZ/M06MLBTcBOBfrV04fFIrVzocaUU0OXbguab08dOG6y88SWtP58IXoBHJRfQ84aNAgXL16FdOnT0deXh46duyITZs2ITq6pmNdvHgRev3fMvU//vEPrFixAq+88gr+/e9/IyEhAevWrTPHmAVqHMhu3ryJUaNGoaioCPfeey82bdrkUTFmLdEqRp4t6Ukxdpo7btJQKn6e5YJiKzAlN6tvV16L99frAMakPYv7PRwdhTnSOPG9qxT6zd8tq3xx2W3M3XIGrYwhksMcbTqRixnf/Wb+e9iSA7JjfwotBJyGVY5GrmFdg0umC+4UYqRk8FIqR4CceMu2bcyRmZ2P/h0bSQ7y7mqYtLmDO8rKiOYM7oqHuXT43bhaWqFKshl3U1xeJTp+5ZrLLNx+DkmNwpCeFCM6XgBrgUuoHztC7rq1dPjdKC67jTP5N7Bg+1nZz+OQGjPXmdi6fNiu6WFB/rzhztyVdtidaHKgkZGRIWhWsGPHDrtrjz/+OB5//HHB79PpdHj11Vfx6quvKlVFVdEqRp7Qs/mOZIWOaeUICpY7aj6BiS9gt9roAIy8Lw4f/5yj+G8t9n1Sd7+2Qr5cGy65iSVcDfrOB7cQcBuhzOw8u5iMfHBCSIN6Th5rQb4QY7vBUiIdJ6chsX1nY1gg+iQZJac6FkNOvGWhNgZqNLyGOnockKC5AlxfdHslRuOh9rHmOLquCn5CY84dGiepG18htJ77+Z4dHOCHW5XVdp+rEdOIm6PExouU9OFSkLvB7RzXAMEBdZB1rsAlYVbqiZyzJ3d8WK7plk5YtvOy2uNDqflUKdxknVP7EBrQUpDSaZRatKVqkXX4e0ERmoSuuOHIbO7gjujfsRHubFrfqd9a7HuV0DgNW3IAUSEBLn2HVEFHTtB3OX2FWwg4LLUFmdl5+Gz3eVEtTFiQ/Pd35thMCY20EEKnHvtzChURZjkcaZykBOb/1zfHJeeJv3ajAtUmJnmu4YNbbG9VhrkszEaFGuw0ylJTjnoiQtp8JZnHM1cZwwIxrW8bvPr9SV5hVmls5yih8cL1M2cTTLiqpXf11DS5WX1J9zt7cicH23lZbZScT5WgVoTm8hTSk2Kwa0oPrBx5jyw7zWFLDtjlArdk04lcpL6/U3J5MTgtMiBs+1Y/2N+s0XMkMKnNZ0M7Wf2W3HF4elIMpvVNVOw5vRKjnWo7Pq6UVipSJ0eCjiN7Ljn2Zo7gBJjp/dpiMY9dl6VdthS7U1ss75cCt8GyFYQ4jbQzY8MWPntlZ95NDEcaJykZ7aQKsgDw1qbTsuYatdkysbvdmHNnRjdXGN61Oe+4mJCaoOhzLOeqeYM7YuXIe7BrSg/Ud9G8xxks5ygx+35nsh0qYRcqZb2z9TGRez9XzhFK2KtrjZLzqauQMKsx3IDu217eTkao06ixaAs6j/3P+ejgK73MQoXcSUjpAwm9Xsf7W1abGGb/kO3wfjkTiLNtpxaOBB2pR8ZS7M3kOJBZbtosF1Ouz1guAI4Y0bW53f1S6upogzVrQ7Yqi4fUxQ0QdxTjc4zhQw1nOLlzjSVCKW2dReqYU/q5ziJWjyW7z2Na3zZ24yKjRwJiJDj1yIFPcHSH46TU439n6iZ3gyuEo/Xul5cfcOp+OQhtEj2lXwuh9nwqBzIz8BL4joXVOkYGhI9Rbb/HmUlIyaM2IYFNqpB96MJ1zR1JXEXq0ZpUOy1HC44zx/VCttoc6UkxWPjEnchYeUTQYUoHYOOJPPy7rzzNixRtpVR7VGcQMynKeCAeCdEhiAoJxPWblRi7QjjUjhSNkxrOcHLnGks4G10tjx2rTQxvbDyp2fNcqcfsH05i15Qedu06o18inv9f2CW10NJxUu7xv9S6TUlvhdjwIMXtQsXWOynJASzvl2srLsXm3ROO8YVQez6VCmlmvQjbY2E5i7YzSAn7JHeCnDu4o90ONixIPP6eGEICm1QhWykvU62Qc7TG2XPJCYtji5rH9fXrGkQ9/53tv1LbXk1NlZBJ0ZgHWprHU5/2whodqXahSps1cMidayzvE9PSKBUKzxKpkRrU5tCF66L1EOvP6Ukxipsb2KJWXxFCzvG/1Lo9dU8zh2EInUVqmENLFm0/a+7rzpzcSdkkeoLWUwruDJkIkDDrlXCdxhMWbbkTpK091/IRXWDwd74bmkyMd6BLFbKV9DLVAjlHa2JH3lKEYrWP69Xqv1LbXm1NlZTFTUjolWoXKtXmz9llX+5cA4hvQDrHNVDcTMfdiyjH3j8KJJUTqm9GjwQYQ9Xrk1LmA6WQ66SnpO2pEM7G2hY76l+w/RySX8t0elOvtkJKS9wd95mEWS+E6zSesGg7M0Fa7oAPXnBNq/LM5wd5bY2kCtnOeJkquUuWEr5s6fC7eW1PpSBkzyVFKFZ7olWr/zoTqN2duGqLLWazN29wR7z5aDsAzgkscucaDjEBU2ktk9S6/bNjrKLPtWXxzj8klROqr59eh5kPJ0IH9RJPiM0Hi564UzHNrTNOekrYngrhrJO0FNORolu3nT6l8gSFlBJ4wnxKNrNehK0dkqOwIloFFxeyEYwODbTznrUMG7bpRC4+2HLG5efzhQCTGttX7k5fLDi9VCzrs2Vid5y4VCJqZyUn5ApfiDap9s+2qD3RqtV/HSU2AHwvM46QzV6vxGgEB9ThHZ96nXCCB6G5Rqqzp5DAlpmdr7h9a+e4BojmCeHFwb3LY52a4Nujl51+TmZ2vkuRFKT0Z1dCOEpFbD7Q63VOxUPXgX++d6VufHOiM2EnXYm1LdWEhTOv6ZVolFwvwDMUUq5iGabTnZBm1kvgW4RdPUZWEj4P9syJ3ezKcTvijccuY9YGx9EGpCCUjlPqTl+qpkjMm3vcqqOYsV5a5ppoi6NEJSMkiGkfnLEHU3uiVbP/uqKR9lbE+hLf+Fww5C5eLaCjuUYMR1pvNdJZ++l1+HefNqJlZvRLxD0tGrqkeRy/6qjL3uVS+rNtW03rK/5uziA0HwiNm5iwQDzXLU60v0zt01rRuvH1Y7lhJ101lVLKvEYIbztFssUyTKe7Ic2slyCUnk5qdhUtsPVgX3/0Em+5vOJyjFlxRJU62KbjdLTTl6opkmKoLyXl47S+bTDw7iZoN3Ozw7JyUCPTlxaafzX7r7MaaaVRIwuZM/BFmPhQL/23T0+KweKn7nKYhIET2PgEBLXcWMQ0ppb2m65mYpyz8ZRzFYQ8O1LLtqo2MXz8yx+aObmJjRu+ZDRcf+l2R6Si9RDaOMiZ01yNbOKMeU21KVRyeSkniO7WevIl4QgP8sfwrs2R0SPB7RpZDhJmPRip6ek8ZdG2RMzWSE2/TL7oBGIZicavOiqpPo4SEEglIsSgeLuoFaJNq+N6NfuvoxBhaqNmFjIlkPvbc+UXbDuLJbtz7DZwnMAmZI6jBUuH342rpRVWphYcrhzjM8ClY/QeraMcF+KB0zyPW3XU6Wc780y+cSPWX0rL/+4Lrm7aHK0fUuc0V02lHJmw2HL+2i1eH45F289iQq9WvHUV65OekO3OMkW1p8gXfJAw68HIsZV096Jti7vC5ciNTiBVsFYqhJcatk9qxlXVSvPvaf1XCdTQlquB3N/eT6/DuNQEZPRoif05hbhYeBNT1h4HUBNmb+OxXIxdwR83Uwu4kwIhG3RLgcyy7mpje2okB0/KesbXX5TetDlaP6TOaa6aSkndSOgAhAX7Y+6W33n7/YLt57Bs30W8+Wg73t9DqE+GBfmLppfW6tTHG+ZnEmYJVVDK+/KzoZ1QWlGFiLoGTFzj2N5OrRzYrobwsjySr6hSNje62o5anqj593TUTGjiKfjpdSguq8T7mb+brw1bcgB6nfKnL0ov2nx1V5u95woUETYyHojHl3svoliCWZMWqLFpU2pOU8JUytFGwvLUStQM7X9RD4R+D6HxJLQp8PRTH60hBzBCFeRoIMWmdr1eh4fax6JrQoRDJw+gRvuhRoBpk4nBGCpuqB/uIPkDny2hlHiHjpxOtPCIdcaBrDbjS/EjhRByiFR6+DkbVsnRdzpKzas0i3/+w+V6A8DC7ec8RpBVKxa1UnOa0k6mGQ/E283zxrBATEhNELUj5xBLKiInQY2ayWy8FRJm3YQ3ZPRwBc7WSAjOS3PRE/aes5bziuXCJbRDFiqvJM98fhDlVdVmrRofs/q3Fbzf0paQb2F+96fTgveOX3VU9H283SPWF3Fn/Ehng8PLfYaU9LaukpmdL2nRlvPOWtWdDyWEDUf1/ulEntPfLRe1Nm1S1w8pc5qSkU3GPNASh6b1sooKsmtKDzSPqCv5O/h+DzmbArWT2XgrJMy6AVuBRgi1FiItEAuXY7kj7tP+7zA0z3RtDsBes8MtAJYayqXD73ZYXsqCIUe/WPy/nXdYsL0Gdu7gjlbCtm3azl6J0YK76dzicny2+7zos8UmJ08K0UbU4K74kWpoMfmQmt7WVd7YeNLhor3xmP07i82vatQ9OMAPgOP5RAthY8JXv2qmmVNr0yZ1/ZA6p/GFppObgMaybranVK4mFZGzKVByA+Gt8gUfJMxqjJzjLbUWIq0Q0qTa7oj99Dp0jmuAHwU0CtxwswyLk9ysvsPyUhcMqSIep5UNrKPH8me7iKYftU3b6Yo2SMrkVBvjqnoy7tCWa3n06KxGWQfAGGrAiz1aShp3UpyAxqywf2e+RCrmz1TQht+qrMaE1DskZa/SwsREK82cmps2qeuHVNQ0leLGu1Rsfw85mwKlNhBSlWreAjmAaYgzAo2neT7LwXIynZLeCsawIBhD+Z2HpOw2LcPiOAqVxeftyje5z+WJoSdGTT0qoNfp0Ld9jKDHtC1KhPZyNDmRo5bnoHUWMq0dzlzRKPfvGIv/bDvrch3EsPwd9p4rQPdWUeb3ViubUvOIYOya0gMfZP6OBdsdv5+aKUqdjV4iF62yUEoNU+kuuPH+/LLDouWEfg81NgViZYWc9rwZ0sxqiDPHW95qA2O763tr02m8vekUissqeSciuRO71FBZ3PeK7UK54ydLTavU75WKEqG9pExk5KjlOWipLdfa4cyR5hmwtmXneH9gB3z3q7YnTc98ftDqhEtK3Z0hKiQQfnoduraMkFxeTdQUljm0MnHqHNfA4+c0LqlIOI8ZGiD+e8g5yXH11MedNuNqQsKshjg7uXib57Mzx51yJ3apobKiQgIdpqH96UQuUuIbIj1Jel5tvvqKOaG4GtrLGGogBy4vRElbPTG0djhzJMToACwYcqfdJrF+3QBZG/r6wf6KCJ2Wc49Y3S2R+lxb4cFTHDLPX7ul6vdzkInT36QnxeDQK70wIfUO3qgHYmG5pG4KXN1AaGXvrjUkzGqIqztxLXbaruKsp6WUBcAY+vfvl9ysvqQFI7lZfYe70IyVR/D6D9mS7IeEFqLM7Hw7J5Se7+0w/+0otJcjyqtMyMzWzkuZUA4ttOXucDhzJMT0aR+LlPiGVrbjck8opj8kvmhLxXbuEaq7LVKfYyk8qK2t/GBgB16tty1zt/yumb+FVps2dyA3OgiXVIQv6oHY7yFnU+DKBsIb5AhnIJtZDXFkX+QItY+llMDZbFRSbAyn9mltzsQi1SZRiq2qiQGf/JLj8N3EFiK+tLhXSivN///M5wcRHuxvtl+U2/7FDgJuE7UbrWwXbZFrpy33hCItyYgPA/x4M9BN69sGs384KXk+tZ17LOv+1/VbdvbvfPb0ep199BS+lKNqZs5LSzLCUMcPY1aI22cC2ibm8IYsUXJxJTGBM7+HnPHkrI+EN8gRzkDCrIaICWBiqLUQqYErx52OFoBud0TKKp+eFIP1Ry85+Sb2iC1EUtrSMrSXlADbtt/vC1mjCHXQ2uHM9tlSF23uREXOhl5s0dbrdbLnU8u5h6v7rcowO2GWLyd9crP6OHThupXwK+R1r6ZDZp/2MZhwJQEfbDkjWMaVNNaE+9JRyxlPzgjMrirVPBUSZjVGSAATwtvihLp63Cm2ANyqrJJVXk59HDGtbxsM6xrnUhvYhvb66bc8fJF1Qdb9tDgRQqipDVQKZzf0Qou23PkUkDcn8D1XSPiVer9SSA3U76vHymriy+monR2Dng4Js27A0fGWJZ60EElBieNOuQuAWHmuPq4avEeEGBSZtCxDe/VOipElzHLQ4kQI4Q3h2QSF7tBAq/B7cr6vV6IRe88VYOyKwygSSPXqTSdcUnBXYo7agLPmct6CM5tAT4eEWTchdrzl6TH1xHDncadYfRzF/3OE0gvCldJyPNQ+1qnjHlqcfBNbR5P7EiKdGifeYLvIJ3QnNQpFu5mbnfo+P70OXRMi8OaAdhj9v7Hu7rlHbdxlJ10bcGc6aq2Qo1TzBiiagQfiDTH1xPC0UC3pSTFY9MSdkjyAbZEaRkfuV3MxKaWECJJbF8L70CoNrSehRpQHT5t71MTX01jLjSKgJLVF682NQcuII94KaWYJVfC0484+7WOxADpJHsC2SF0QpNgf2WpLpB73+MLiRPDjLkcTT0MpYUWLuUcpLbqrz/MGO2lncCWKgBKQ1tv7IGGWUA1PO+7s0z4Gi/X2E39MWCAe7hCD737NtRMo+cLu8CElLa6QQGq7+J6/dgsr91+0sh/09sWJ4MeXHU3kYCu8cGRm56N/x0ayv0/NuUdrQcvR8zxNceAqnrC58zRzOcIxqgmzhYWFeOGFF7Bhwwbo9XoMGDAA8+bNQ7169QTvKS8vx6RJk7Bq1SpUVFQgLS0NixYtQnT036FPdDr7zrNy5UoMHjxYlfcgfAuxif/l9DZ29kNCYXdssQ3jI1cgtV18M3q09JnFiRDG1x1NpCCWJ378qqMw1NF7zCZOa0FL6vM8TXHgLJ60ufNVrbevopow++STTyI3NxeZmZm4ffs2hg8fjlGjRmHFihWC90yYMAE//PAD1qxZg7CwMGRkZODRRx/F7t27rcotWbIE6enp5r/Dw8PVeg3CBxGa+LnrpeWhZmFWzhGikgKpryxOhDi1wdFEDCl54rUQXqSYOGgtaHmSYKcVnra58zWttxBa2iOrhSrC7MmTJ7Fp0yYcOHAAnTp1AgDMnz8fffr0wbvvvovY2Fi7e4qLi/Hf//4XK1asQI8ePQDUCK1t2rTB3r17cc8995jLhoeHw2g0qlF1woPRwk5NySNEEkgJR9QWRxMhPEF4kWrioHVdPeG30RpP3Nz5+jwu1P+9DVWiGWRlZSE8PNwsyAJAamoq9Ho99u3bx3vPoUOHcPv2baSmppqvtW7dGk2bNkVWVpZV2bFjxyIiIgKdO3fGZ599Bsa8f1dBiKOFtzd3pJdfYp0/njvS82XPcsI9cI4mQlsyX49gIVUoyczOU+X5QmMeqDFxsBzzWgtanijYqU1t39xpjVj/9zZUEWbz8vIQFRVlda1OnTpo0KAB8vL4J6W8vDwEBATYmQxER0db3fPqq6/iq6++QmZmJgYMGIAxY8Zg/vz5ovWpqKhASUmJ1T/Ce9BCyHR0pAfUHOn5wnEM4Tn4englR0gVSj7bfV7xzaRUEwduzGstaNVGwa62b+60REr/58p5A7KE2X/961/Q6XSi/06dOqVWXQEA06ZNQ9euXXHnnXdiypQpePnll/HOO++I3jNnzhyEhYWZ/zVp0kTVOhLKoZWQKedIjyCUpDbFRrWFE14cwdmHKrmwyh3zWgtatVGwq+2bOy1x1P85Dl24rkFtXEeWMDtp0iScPHlS9F+LFi1gNBpx5coVq3urqqpQWFgoaOtqNBpRWVmJoqIiq+v5+fmi9rFdunTBX3/9hYoKYTX51KlTUVxcbP73559/Sn9pwq1oJWTWxiM9wnNIT4rBrik9sHLkPZg3uCNWjrwHu6b08GlBFrAWXsRQYzMpd8xrLWjVVsGuNm/utERq/79a6h0mCLIcwCIjIxEZGemwXEpKCoqKinDo0CEkJycDALZt2waTyYQuXbrw3pOcnAx/f39s3boVAwYMAACcPn0aFy9eREpKiuCzjh49ivr168NgMAiWMRgMop8TnotWQmZtPNIjPAtfdzQRIj0pBiO6Nsd/d593WFbJzaQzY17rcE21NTxUbYki4E6k9v/IEO+QnVSJZtCmTRukp6dj5MiRWLx4MW7fvo2MjAwMHjzYHMng0qVL6NmzJ7744gt07twZYWFhGDFiBCZOnIgGDRogNDQUL7zwAlJSUsyRDDZs2ID8/Hzcc889CAwMRGZmJt544w1MnjxZjdcgPACthEzK+EIQ6iAlCklqolGSMKvkZtLZMa+1oFVbBbvaurnTCkf9nyO5WX3N6uQKqjiAAcDy5cvRunVr9OzZE3369MG9996Ljz/+2Pz57du3cfr0ady6dct87YMPPsBDDz2EAQMGoFu3bjAajfjmm2/Mn/v7+2PhwoVISUlBx44d8dFHH+H999/HjBkz1HoNws1oZTdWW4/0CEJNpEYhcYd9qNiY5xAa85yg1b9jI6TEN1R9XtD6eYTvI6X/c+W8AR2rhXGtSkpKEBYWhuLiYoSGhrq1Lrcqq5A4/Sera9mvpiE4gDINc3DRDAD+tIJK2lFtOpHLm+6W70jPsu2ozQjCGqHsVULjVstxbltP2zEPAPMGd3QqlS5BeBNC/Z/D3eNAqrxGqy/h8WhpN1Zbj/QIQkmcyV7lLvtQyzHvTCprgvBm0pNiYDIBY1Yc5v183Kqj8Nfr0Ke9fbIrT4KEWQ9ErexW3oyWQibZahGEazibvcpdm0luzN+qDDMLswRRG6g2Mcz+IVu0TMbKI1gAHfq091yHQxJmPRBXUqj6MiRkEoR34EoUEhrnBKEdUuLNmliN5nax3nNDo6nmAEZIIzM7n/c6pVAlCMJboVB3BOEdyAl358lZMEmYdSPVJoY3Np7k/YxSqHo+tiGHqJ0IoobamL2KILwRORtKT86CScKsG9mfU4j8EuHsGpRC1XORGnKIIGojFOqOILwDqSmlOTw1CyYJs26EUqh6J1wIIduNCJmGEMTfUFpSgvB8pKaU5vBU0yByAHMjZFfmfTgTcoggaisU6o4gPJ/0pBgseuJOZKw8AiFrOU/PgknCrBvpHNcA0aEGQVMDT+88tRFnQw4RRG2FohMQhOfTp30sFkDHG2/WG0yDyMzAjfjpdfh3nza8n3lD56mNkGkIQRAE4Yv0aR+DxU/dZWdD6w2mQaSZdTNCWWbUznpDOAeZhhAEQRC+ireaBpEw60EsHX43istue03nqY1wnp95xeW8drNkGkIQBEF4M95oGkTCrAfROa4BggOoSTwZzvNz9LLD0AFWAi2ZhhAEQRCE9pDNLEHIhEIOEQRBEITnQGpAgnACb7UrIghCGNusfvclRNKYJggvgIRZgnASb7QrIgiCn00ncjHju9/Mfw9bcgAx5IhLEF4BmRkQBEEQtRrK6kcQ3g0JswRBEEStxVFWP6Amq1+1UGokgiDcDgmzBEEQRK1FTlY/giA8ExJmCYIgiFoLZfUjCO+HhFk3Y+s9S0dZBEEQ2kFZ/QjC+yFh1o1sOpGL1Pd3mv8etuQA7n1rGzkbEARBaASX1U8oAJcOQAxl9SMIj4aEWTdB3rMEQRDuh8vqB8BOoKWsfgThHZAw6wbIe5YgCMJzoKx+BOHdUNIENyDHe5aC8hMEQagPZfUjCO+FhFk3QN6zBEEQngdl9SMI74TMDNwAec8SBEEQBEEoAwmzboC8ZwmCIAiCIJSBhFk3QN6zBEEQBEEQykDCrJsg71mCIAiCIAjXIQcwN0LeswRBEARBEK5RK4VZxmrit5aUlLi5JjW0jfRH20h/AMDNG6Vurg1BEARBEIT74eQ0Tm4TolYKs6WlNQJjkyZN3FwTgiAIgiAIQozS0lKEhYUJfq5jjsRdH8RkMuHy5csICQmBTqf+kX5JSQmaNGmCP//8E6Ghoao/j3AMtYnnQW3ieVCbeB7UJp4HtYl6MMZQWlqK2NhY6PXCbl61UjOr1+vRuHFjzZ8bGhpKHd3DoDbxPKhNPA9qE8+D2sTzoDZRBzGNLAdFMyAIgiAIgiC8FhJmCYIgCIIgCK+FhFkNMBgMmDFjBgwGg7urQvwPahPPg9rE86A28TyoTTwPahP3UysdwAiCIAiCIAjfgDSzBEEQBEEQhNdCwixBEARBEAThtZAwSxAEQRAEQXgtJMwSBEEQBEEQXgsJsxqwcOFCNG/eHIGBgejSpQv279/v7irVCmbOnAmdTmf1r3Xr1ubPy8vLMXbsWDRs2BD16tXDgAEDkJ+f78Ya+yY///wz+vXrh9jYWOh0Oqxbt87qc8YYpk+fjpiYGAQFBSE1NRVnzpyxKlNYWIgnn3wSoaGhCA8Px4gRI3Djxg0N38K3cNQmw4YNsxs76enpVmWoTZRjzpw5uPvuuxESEoKoqCg88sgjOH36tFUZKfPVxYsX0bdvXwQHByMqKgovvfQSqqqqtHwVn0FKm9x///124+T555+3KkNtog0kzKrM6tWrMXHiRMyYMQOHDx9Ghw4dkJaWhitXrri7arWCtm3bIjc31/xv165d5s8mTJiADRs2YM2aNdi5cycuX76MRx991I219U1u3ryJDh06YOHChbyfv/322/jPf/6DxYsXY9++fahbty7S0tJQXl5uLvPkk0/it99+Q2ZmJr7//nv8/PPPGDVqlFav4HM4ahMASE9Ptxo7K1eutPqc2kQ5du7cibFjx2Lv3r3IzMzE7du38eCDD+LmzZvmMo7mq+rqavTt2xeVlZXYs2cPPv/8cyxduhTTp093xyt5PVLaBABGjhxpNU7efvtt82fUJhrCCFXp3LkzGzt2rPnv6upqFhsby+bMmePGWtUOZsyYwTp06MD7WVFREfP392dr1qwxXzt58iQDwLKysjSqYe0DAPv222/Nf5tMJmY0Gtk777xjvlZUVMQMBgNbuXIlY4yx7OxsBoAdOHDAXObHH39kOp2OXbp0SbO6+yq2bcIYY0OHDmX9+/cXvIfaRF2uXLnCALCdO3cyxqTNVxs3bmR6vZ7l5eWZy3z44YcsNDSUVVRUaPsCPohtmzDGWPfu3dm4ceME76E20Q7SzKpIZWUlDh06hNTUVPM1vV6P1NRUZGVlubFmtYczZ84gNjYWLVq0wJNPPomLFy8CAA4dOoTbt29btU3r1q3RtGlTahsNycnJQV5enlU7hIWFoUuXLuZ2yMrKQnh4ODp16mQuk5qaCr1ej3379mle59rCjh07EBUVhVatWmH06NEoKCgwf0Ztoi7FxcUAgAYNGgCQNl9lZWWhXbt2iI6ONpdJS0tDSUkJfvvtNw1r75vYtgnH8uXLERERgaSkJEydOhW3bt0yf0Ztoh113F0BX+batWuorq626sgAEB0djVOnTrmpVrWHLl26YOnSpWjVqhVyc3Mxa9Ys3HfffThx4gTy8vIQEBCA8PBwq3uio6ORl5fnngrXQrjfmm+McJ/l5eUhKirK6vM6deqgQYMG1FYqkZ6ejkcffRRxcXE4d+4c/v3vf6N3797IysqCn58ftYmKmEwmjB8/Hl27dkVSUhIASJqv8vLyeMcR9xnhPHxtAgBPPPEEmjVrhtjYWBw7dgxTpkzB6dOn8c033wCgNtESEmYJn6V3797m/2/fvj26dOmCZs2a4auvvkJQUJAba0YQns3gwYPN/9+uXTu0b98e8fHx2LFjB3r27OnGmvk+Y8eOxYkTJ6zs+wn3ItQmljbi7dq1Q0xMDHr27Ilz584hPj5e62rWasjMQEUiIiLg5+dn53Gan58Po9HoplrVXsLDw3HHHXfg7NmzMBqNqKysRFFRkVUZahtt4X5rsTFiNBrtHCarqqpQWFhIbaURLVq0QEREBM6ePQuA2kQtMjIy8P3332P79u1o3Lix+bqU+cpoNPKOI+4zwjmE2oSPLl26AIDVOKE20QYSZlUkICAAycnJ2Lp1q/mayWTC1q1bkZKS4saa1U5u3LiBc+fOISYmBsnJyfD397dqm9OnT+PixYvUNhoSFxcHo9Fo1Q4lJSXYt2+fuR1SUlJQVFSEQ4cOmcts27YNJpPJvHgQ6vLXX3+hoKAAMTExAKhNlIYxhoyMDHz77bfYtm0b4uLirD6XMl+lpKTg+PHjVpuMzMxMhIaGIjExUZsX8SEctQkfR48eBQCrcUJtohHu9kDzdVatWsUMBgNbunQpy87OZqNGjWLh4eFW3o2EOkyaNInt2LGD5eTksN27d7PU1FQWERHBrly5whhj7Pnnn2dNmzZl27ZtYwcPHmQpKSksJSXFzbX2PUpLS9mRI0fYkSNHGAD2/vvvsyNHjrALFy4wxhh78803WXh4OFu/fj07duwY69+/P4uLi2NlZWXm70hPT2d33nkn27dvH9u1axdLSEhgQ4YMcdcreT1ibVJaWsomT57MsrKyWE5ODtuyZQu76667WEJCAisvLzd/B7WJcowePZqFhYWxHTt2sNzcXPO/W7dumcs4mq+qqqpYUlISe/DBB9nRo0fZpk2bWGRkJJs6dao7XsnrcdQmZ8+eZa+++io7ePAgy8nJYevXr2ctWrRg3bp1M38HtYl2kDCrAfPnz2dNmzZlAQEBrHPnzmzv3r3urlKtYNCgQSwmJoYFBASwRo0asUGDBrGzZ8+aPy8rK2Njxoxh9evXZ8HBweyf//wny83NdWONfZPt27czAHb/hg4dyhirCc81bdo0Fh0dzQwGA+vZsyc7ffq01XcUFBSwIUOGsHr16rHQ0FA2fPhwVlpa6oa38Q3E2uTWrVvswQcfZJGRkczf3581a9aMjRw50m4DTm2iHHxtAYAtWbLEXEbKfHX+/HnWu3dvFhQUxCIiItikSZPY7du3NX4b38BRm1y8eJF169aNNWjQgBkMBtayZUv20ksvseLiYqvvoTbRBh1jjGmnByYIgiAIgiAI5SCbWYIgCIIgCMJrIWGWIAiCIAiC8FpImCUIgiAIgiC8FhJmCYIgCIIgCK+FhFmCIAiCIAjCayFhliAIgiAIgvBaSJglCIIgCIIgvBYSZgmCIAiCIAivhYRZgiAIgiAIwmshYZYgCIIgCILwWkiYJQiCIAiCILwWEmYJgiAIgiAIr+X/A/sg/49UxzBBAAAAAElFTkSuQmCC",
|
|
"text/plain": [
|
|
"<Figure size 800x200 with 1 Axes>"
|
|
]
|
|
},
|
|
"metadata": {},
|
|
"output_type": "display_data"
|
|
}
|
|
],
|
|
"source": [
|
|
"# Simple linear regression\n",
|
|
"from scipy import stats\n",
|
|
"slope, intercept, r_value, p_value, std_err = stats.linregress(np.arange(tsa_a.shape[0]), tsa_a)\n",
|
|
"\n",
|
|
"plt.figure(figsize=(8,2))\n",
|
|
"plt.stem(tsa_a - np.arange(tsa_a.shape[0]) * slope - intercept)\n",
|
|
"plt.title('linear regression on a (time stability vs. global beat)'); None"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 6,
|
|
"id": "db46ce9e-e237-4501-a90a-191e32f1d0a3",
|
|
"metadata": {},
|
|
"outputs": [
|
|
{
|
|
"data": {
|
|
"text/plain": [
|
|
"(\"primary:Music/Rhythmology/01 Titan's Call.mp3\",\n",
|
|
" \"Titan's Call\",\n",
|
|
" 'Follow The Cipher')"
|
|
]
|
|
},
|
|
"execution_count": 6,
|
|
"metadata": {},
|
|
"output_type": "execute_result"
|
|
}
|
|
],
|
|
"source": [
|
|
"annot_a['contentId'], annot_a['title'], annot_a['artist']"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 7,
|
|
"id": "db35c8f5-4258-47e1-ab19-64d5c0667930",
|
|
"metadata": {},
|
|
"outputs": [],
|
|
"source": [
|
|
"#import os\n",
|
|
"#os.makedirs('ref')"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 9,
|
|
"id": "d9147484-6a1d-4639-9da9-8376e73b479f",
|
|
"metadata": {},
|
|
"outputs": [],
|
|
"source": [
|
|
"with open('ref/Rhythmology_001.json', 'w') as fout:\n",
|
|
" json.dump({\n",
|
|
" 'contentId': annot_a['contentId'],\n",
|
|
" 'title': annot_a['title'],\n",
|
|
" 'artist': annot_a['artist'],\n",
|
|
" 'beatTimesSec': list(tsa_a)\n",
|
|
" }, fout, indent=2)"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": null,
|
|
"id": "85bb17d3-7c61-43c1-8700-135b8b1bd1cf",
|
|
"metadata": {},
|
|
"outputs": [],
|
|
"source": []
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 5,
|
|
"id": "1cdb8c06-406d-4861-a1d9-2a6809d92834",
|
|
"metadata": {},
|
|
"outputs": [
|
|
{
|
|
"data": {
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAqsAAADcCAYAAABXo+ihAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABaQklEQVR4nO3deXxM5/4H8M9MJJNENhHZ0CSWIkIQkqa2lpCgSjdLKXJ7tURaa6/q/SFor7UttXalt6i1pVpSia2XprZQJSgaS8lCIgvZM8/vDz2nM5NZzpn1zOT7fr282pw5c85z1vme5zzP95ExxhgIIYQQQgiRILmtC0AIIYQQQoguFKwSQgghhBDJomCVEEIIIYRIFgWrhBBCCCFEsihYJYQQQgghkkXBKiGEEEIIkSwKVgkhhBBCiGRRsEoIIYQQQiSLglVCCCGEECJZFKwSs9mwYQNkMhmuX7/OT3vqqafw1FNP2axMRLxx48YhNDTU1sWwGO48PXXqlEnLSUpKQr9+/QTN6+j7VBeZTIbk5GSD8wm5d1y/fh0ymQwbNmwwf0HtyOHDhyGTyXD48GHR3x03bhw8PDzMWp7Q0FCMGzdO7zzcsVu2bJlZ122qESNGYNiwYbYuBhGAglVCCBEpOzsbn332Gd555x1+2p07d5CSkoKzZ8/armBmtHfvXqSkpNi6GHrZQxmJbem7LmfOnImdO3fi119/tX7BiCgUrBKL2r9/P/bv32/rYhARPv30U1y+fNnWxZC0FStWICwsDE8//TQ/7c6dO5g3b57WH0V73Kd79+7FvHnzrLKuV155BeXl5QgJCdE5T0hICMrLy/HKK6/w06xZRmKf9F2XnTt3RteuXfH+++9bv2BEFApWiUW5uLjAxcXF1sXQ6eHDh6Lmr6iogFKptFBp1CmVSlRUVFhlXaqcnZ2hUCisvl57UV1djU2bNol6fUj7VD8nJye4urpCJpPpnEcmk8HV1RVOTk5WLBlxdMOGDcM333yDBw8e2LooRA8KVolFabY749pbbdu2De+99x6aNWsGV1dX9O3bF1evXq3z/ePHjyMhIQHe3t5wd3dH7969cezYMbV5bty4gaSkJLRp0wZubm5o3LgxXnrpJbX2b8Df7eKOHDmCpKQk+Pv7o1mzZjrLzpV1y5Yt+L//+z80bdoU7u7uKCkpEVw2bjldu3aFq6srWrZsiY8//hgpKSl1fpi59n2bNm1C+/btoVAokJqaCgC4ffs2/vGPfyAgIAAKhQLt27fHF198UWddK1euRPv27eHu7o5GjRqha9eu2Lx5M/95aWkppkyZgtDQUCgUCvj7+6Nfv37IzMzk59HWvvLhw4eYPn06mjdvDoVCgTZt2mDZsmVgjGndhl27diEiIoIvK7cdhuTn5+PVV19FQEAAXF1dERkZiS+//FJtHtX2b5988glatmwJhUKBbt264eTJk4LWAwBlZWV4/fXX0bhxY3h5eWHMmDG4f/++we8dPXoU9+7dQ1xcHD/t8OHD6NatGwAgMTERMplMrX2l5j5V3YbVq1ejRYsWcHd3R//+/XHr1i0wxrBgwQI0a9YMbm5uGDJkCAoLC+uUZd++fejZsycaNmwIT09PDBo0CBcuXDC4DdXV1Zg3bx5at24NV1dXNG7cGD169EBaWhpf3tWrVwMAvy2q5+uyZcvw5JNPonHjxnBzc0NUVBR27Nihc32bNm1CmzZt4OrqiqioKPz0009qn2trs6pJs82qrjIyxhAaGoohQ4bUWUZFRQW8vb3x+uuv61xPRESEWo05R6lUomnTpnjxxRf5aVu2bEFUVBQ8PT3h5eWFDh06YMWKFTqXrY9SqURKSgqCg4Ph7u6Op59+GllZWYLahALA9u3bERUVBTc3N/j5+WH06NG4ffu21nn/+OMPxMfHo2HDhggODsb8+fPrXMtij7FQH374IUJCQuDm5obevXvj/Pnzdea5dOkSXnzxRfj6+sLV1RVdu3bFd999pzZPYWEhZsyYgQ4dOsDDwwNeXl4YMGCA2it9Q9clAPTr1w8PHz7kz30iTQ1sXQBSPy1atAhyuRwzZsxAcXExlixZglGjRuH48eP8PAcPHsSAAQMQFRWFuXPnQi6XY/369ejTpw/+97//ITo6GgBw8uRJ/PzzzxgxYgSaNWuG69evY+3atXjqqaeQlZUFd3d3tXUnJSWhSZMmmDNnjqCa1QULFsDFxQUzZsxAZWUlXFxcBJftzJkzSEhIQFBQEObNm4fa2lrMnz8fTZo00bqugwcPYtu2bUhOToafnx9CQ0ORl5eHJ554gg8EmzRpgn379uHVV19FSUkJpkyZAuDRq+Y333wTL774IiZPnoyKigqcO3cOx48fx8svvwwAmDBhAnbs2IHk5GSEh4ejoKAAR48excWLF9GlSxetZWKM4dlnn8WhQ4fw6quvolOnTvjxxx/x1ltv4fbt2/jwww/V5j969Ci++eYbJCUlwdPTEx999BFeeOEF3Lx5E40bN9a5n8vLy/HUU0/h6tWrSE5ORlhYGLZv345x48ahqKgIkydPVpt/8+bNKC0txeuvvw6ZTIYlS5bg+eefxx9//AFnZ2eDxzU5ORk+Pj5ISUnB5cuXsXbtWty4cYN/SNHl559/hkwmQ+fOnflp7dq1w/z58zFnzhy89tpr6NmzJwDgySef1FuGTZs2oaqqCm+88QYKCwuxZMkSDBs2DH369MHhw4cxc+ZMXL16FStXrsSMGTPUHlC++uorjB07FvHx8Vi8eDHKysqwdu1a9OjRA2fOnNHboSslJQULFy7EP//5T0RHR6OkpASnTp1CZmYm+vXrh9dffx137txBWloavvrqqzrfX7FiBZ599lmMGjUKVVVV2LJlC1566SV8//33GDRokNq8R44cwdatW/Hmm29CoVBgzZo1SEhIwIkTJxAREaF3/+ijq4wymQyjR4/GkiVLUFhYCF9fX/6zPXv2oKSkBKNHj9a53OHDhyMlJQW5ubkIDAzkpx89ehR37tzBiBEjAABpaWkYOXIk+vbti8WLFwMALl68iGPHjtU5V4WYNWsWlixZgsGDByM+Ph6//vor4uPjBb1d2bBhAxITE9GtWzcsXLgQeXl5WLFiBY4dO4YzZ87Ax8eHn7e2thYJCQl44oknsGTJEqSmpmLu3LmoqanB/Pnz+fnEHGOh/vvf/6K0tBSTJk1CRUUFVqxYgT59+uC3335DQEAAAODChQvo3r07mjZtirfffhsNGzbEtm3bMHToUOzcuRPPPfccgEcB965du/DSSy8hLCwMeXl5+Pjjj9G7d29kZWUhODhY0HUZHh4ONzc3HDt2jF82kSBGiJmsX7+eAWDZ2dn8tN69e7PevXvzfx86dIgBYO3atWOVlZX89BUrVjAA7LfffmOMMaZUKlnr1q1ZfHw8UyqV/HxlZWUsLCyM9evXT22apoyMDAaA/fe//61Tvh49erCamhqD28OVtUWLFmrrEFO2wYMHM3d3d3b79m1+2pUrV1iDBg2Y5uUHgMnlcnbhwgW16a+++ioLCgpi9+7dU5s+YsQI5u3tzZdtyJAhrH379nq3ydvbm02aNEnvPGPHjmUhISH837t27WIA2Lvvvqs234svvshkMhm7evWq2ja4uLioTfv1118ZALZy5Uq9612+fDkDwDZu3MhPq6qqYrGxsczDw4OVlJQwxhjLzs5mAFjjxo1ZYWEhP+/u3bsZALZnzx696+HOg6ioKFZVVcVPX7JkCQPAdu/erff7o0ePZo0bN64z/eTJkwwAW79+fZ3PNPcptw1NmjRhRUVF/PRZs2YxACwyMpJVV1fz00eOHMlcXFxYRUUFY4yx0tJS5uPjw8aPH6+2ntzcXObt7V1nuqbIyEg2aNAgvfNMmjSpzjnK0bzmqqqqWEREBOvTp4/adAAMADt16hQ/7caNG8zV1ZU999xz/DQh9w5un6nuX11lvHz5MgPA1q5dqzb92WefZaGhoWrXra7vap6vSUlJzMPDg9/2yZMnMy8vL0H3EkNyc3NZgwYN2NChQ9Wmp6SkMABs7Nix/DTuvnTo0CHG2KN97+/vzyIiIlh5eTk/3/fff88AsDlz5vDTxo4dywCwN954g5+mVCrZoEGDmIuLC7t79y4/XegxDgkJUSufNtyxc3NzY3/++Sc//fjx4wwAmzp1Kj+tb9++rEOHDvy5zpXxySefZK1bt+anVVRUsNra2jrrUSgUbP78+fw0fdcl5/HHH2cDBgzQuw3EtqgZALGJxMREtbas3BPvH3/8AQA4e/Ysrly5gpdffhkFBQW4d+8e7t27h4cPH6Jv37746aef+Lajbm5u/HKqq6tRUFCAVq1awcfHR+31Nmf8+PGi2r2NHTtWbR1Cy1ZbW4v09HQMHToUwcHB/PdbtWqFAQMGaF1X7969ER4ezv/NGMPOnTsxePBgMMb4dd27dw/x8fEoLi7mt9HHxwd//vmn3lfhPj4+OH78OO7cuSN4+/fu3QsnJye8+eabatOnT58Oxhj27dunNj0uLg4tW7bk/+7YsSO8vLz4Y6tvPYGBgRg5ciQ/zdnZGW+++SYePHiAI0eOqM0/fPhwNGrUiP9b8xwy5LXXXlOrgZ04cSIaNGiAvXv36v1eQUGB2npN8dJLL8Hb25v/OyYmBgAwevRoNGjQQG16VVUV/1o3LS0NRUVFGDlypNo54eTkhJiYGBw6dEjven18fHDhwgVcuXLFqHKrXg/3799HcXExevbsqfV6i42NRVRUFP/3Y489hiFDhuDHH39EbW2tUes35PHHH0dMTAw2bdrETyssLMS+ffswatQovTXnjz/+ODp16oStW7fy02pra7Fjxw4MHjyY33YfHx+zvT4+cOAAampqkJSUpDb9jTfeMPjdU6dOIT8/H0lJSXB1deWnDxo0CG3btsUPP/xQ5zuq6cS4NzZVVVVIT0/np4s5xkINHToUTZs25f+Ojo5GTEwMf80VFhbi4MGDGDZsGEpLS/nzuqCgAPHx8bhy5Qp/DSgUCsjlj0KY2tpaFBQUwMPDA23atBFdxkaNGuHevXtGbxexPGoGQGziscceU/ub+/Hn2gxyP6Jjx47VuYzi4mI0atQI5eXlWLhwIdavX4/bt2+rtb0qLi6u872wsDBRZdWcX2jZKioqUF5ejlatWtX5XNs0beu6e/cuioqK8Mknn+CTTz7R+p38/HwAj9KwpKenIzo6Gq1atUL//v3x8ssvo3v37vy8S5YswdixY9G8eXNERUVh4MCBGDNmDFq0aKFzW27cuIHg4GB4enqqTW/Xrh3/uSrNYws8Or6G2oPeuHEDrVu35n+AxK5H8xwypHXr1mp/e3h4ICgoSG+7SQ7TaN9nLM1t4ALX5s2ba52ueX306dNH63K9vLz0rnf+/PkYMmQIHn/8cURERCAhIQGvvPIKOnbsKKjc33//Pd59912cPXsWlZWV/HRtQaDmfgYeBYRlZWW4e/eu2qt2cxozZgySk5Nx48YNhISEYPv27aiurlbLJqDL8OHD8c477+D27dto2rQpDh8+jPz8fAwfPpyfJykpCdu2bcOAAQPQtGlT9O/fH8OGDUNCQoLosnLntuZ9wdfX1+CDEffdNm3a1Pmsbdu2OHr0qNo0uVxe53p//PHHAUDt3BdzjIXSdS5s27YNAHD16lUwxjB79mzMnj1b6zLy8/PRtGlTKJVKrFixAmvWrEF2drbag4++5kbaMMZM2i5ieRSsEpvQVbPJBQFcrenSpUvRqVMnrfNyya3feOMNrF+/HlOmTEFsbCy8vb0hk8kwYsQIrT33VWsMhNCcX2jZjOnJr2tdo0eP1hkccwFGu3btcPnyZXz//fdITU3Fzp07sWbNGsyZM4dP7zNs2DD07NkT3377Lfbv34+lS5di8eLF+Oabb3TW9opl6Niai7XWo6lx48aCA2JDdG2D0Ovjq6++0hrsqdbKatOrVy9cu3YNu3fvxv79+/HZZ5/hww8/xLp16/DPf/5T73f/97//4dlnn0WvXr2wZs0aBAUFwdnZGevXr1frzGdrI0aMwNSpU7Fp0ya888472LhxI7p27ao1qNM0fPhwzJo1C9u3b8eUKVOwbds2eHt7qwWi/v7+OHv2LH788Ufs27cP+/btw/r16zFmzJg6nQLtja2OMXdez5gxA/Hx8Vrn4QL6//znP5g9ezb+8Y9/YMGCBfD19YVcLseUKVNEZ2y5f/++1kCaSAcFq0SSuNfIXl5ear2utdmxYwfGjh2rliuvoqICRUVFNi2bv78/XF1dtWY50DZNmyZNmsDT0xO1tbUG9wMANGzYEMOHD8fw4cNRVVWF559/Hu+99x5mzZrFvyIMCgpCUlISkpKSkJ+fjy5duuC9997TGayGhIQgPT0dpaWlarWrly5d4j83h5CQEJw7dw5KpVKtdtXc6+FcuXJFrdf3gwcPkJOTg4EDB+r9Xtu2bbFp0yYUFxervcK3Zs0Mdw76+/sLOi+08fX1RWJiIhITE/HgwQP06tULKSkpfLCqa3t27twJV1dX/Pjjj2rpuNavX691fm1NDX7//Xe4u7vr7GgolL597uvri0GDBmHTpk0YNWoUjh07huXLlwtablhYGKKjo7F161YkJyfjm2++wdChQ+ukH3NxccHgwYMxePBgKJVKJCUl4eOPP8bs2bN1vj3Rhju3r169qvZ2paCgwOCDEffdy5cv16lpv3z5cp3rRqlU4o8//uBrU4FHxwMA3ylP7DEWSte5wK2Xq/F1dnYWdN9/+umn8fnnn6tNLyoqgp+fH/+3oeuypqYGt27dwrPPPitkE4iNUJtVIklRUVFo2bIlli1bpjX/3d27d/n/d3JyqlObtnLlSou1hxNaNicnJ8TFxWHXrl1qbUSvXr1ap52nLk5OTnjhhRewc+dOrSleVPdDQUGB2mcuLi4IDw8HYwzV1dWora2t0yzC398fwcHBaq/5NA0cOBC1tbVYtWqV2vQPP/wQMpnMbDWyAwcORG5urlpbwZqaGqxcuRIeHh7o3bu3WdbD+eSTT1BdXc3/vXbtWtTU1BjcntjYWDDGcPr0abXpDRs2BACLPSSpio+Ph5eXF/7zn/+obQNH9bzQRvNc8fDwQKtWrdTOA13b4+TkBJlMpnZ9Xb9+Hbt27dK6royMDLU2hLdu3cLu3bvRv39/k3OmGtrnr7zyCrKysvDWW2/BycmJ78kvxPDhw/HLL7/giy++wL1799SaAAB196FcLuffcnD7sbq6GpcuXUJOTo7edfXt2xcNGjTA2rVr1aZrXnPadO3aFf7+/li3bp3a8du3bx8uXryotee+6nIZY1i1ahWcnZ3Rt29fAOKPsVC7du1SS6d14sQJHD9+nL/m/P398dRTT+Hjjz/Wus8M3fe3b99eJ12XoXMkKysLFRUVBjN3ENuimlUiSXK5HJ999hkGDBiA9u3bIzExEU2bNsXt27dx6NAheHl5Yc+ePQCAZ555Bl999RW8vb0RHh6OjIwMpKeni263ZImypaSkYP/+/ejevTsmTpzIB30RERGCh+VctGgRDh06hJiYGIwfPx7h4eEoLCxEZmYm0tPT+fyb/fv3R2BgILp3746AgABcvHgRq1atwqBBg+Dp6YmioiI0a9YML774IiIjI+Hh4YH09HScPHlS7wgugwcPxtNPP41///vfuH79OiIjI7F//37s3r0bU6ZMUetMZYrXXnsNH3/8McaNG4fTp08jNDQUO3bs4GvENNvMmqqqqgp9+/bFsGHDcPnyZaxZswY9evQwWMPSo0cPNG7cGOnp6Wo1WS1btoSPjw/WrVsHT09PNGzYEDExMaLbSAvh5eWFtWvX4pVXXkGXLl0wYsQINGnSBDdv3sQPP/yA7t276w10wsPD8dRTTyEqKgq+vr44deoUn9KMw3WKevPNNxEfH88He4MGDcIHH3yAhIQEvPzyy8jPz8fq1avRqlUrnDt3rs66IiIiEB8fr5a6CoBZRp7SVUbOoEGD0LhxY2zfvh0DBgyAv7+/4GUPGzYMM2bMwIwZM+Dr61unpu+f//wnCgsL0adPHzRr1gw3btzAypUr0alTJ76d9e3bt9GuXTuMHTtWLbenpoCAAEyePBnvv/8+nn32WSQkJODXX3/Fvn374Ofnp7d20NnZGYsXL0ZiYiJ69+6NkSNH8qmrQkNDMXXqVLX5XV1dkZqairFjxyImJgb79u3DDz/8gHfeeYev6RZ7jIVq1aoVevTogYkTJ6KyshLLly9H48aN8a9//YufZ/Xq1ejRowc6dOiA8ePHo0WLFsjLy0NGRgb+/PNPPo/qM888g/nz5yMxMRFPPvkkfvvtN2zatKlOe1xD12VaWhrc3d3Rr18/o7eLWIH1ExAQRyUmddX27dvVvqstLQ1jjJ05c4Y9//zzrHHjxkyhULCQkBA2bNgwduDAAX6e+/fvs8TERObn58c8PDxYfHw8u3TpUp2UKlz5Tp48KWh7dJVVTNkYY+zAgQOsc+fOzMXFhbVs2ZJ99tlnbPr06czV1VVtPgA600rl5eWxSZMmsebNmzNnZ2cWGBjI+vbtyz755BN+no8//pj16tWLL0/Lli3ZW2+9xYqLixljjFVWVrK33nqLRUZGMk9PT9awYUMWGRnJ1qxZo7YuzTRLjD1KlTR16lQWHBzMnJ2dWevWrdnSpUvrpADStQ1C0ttw28kdSxcXF9ahQ4c65wR3rixdurTO9wGwuXPn6l0Hdx4cOXKEvfbaa6xRo0bMw8ODjRo1ihUUFBgsI2OMvfnmm6xVq1Z1pu/evZuFh4fzqcm4sutKXaW5DbrOOV3n7qFDh1h8fDzz9vZmrq6urGXLlmzcuHFqqaK0effdd1l0dDTz8fFhbm5urG3btuy9995TS+VVU1PD3njjDdakSRMmk8nUUkR9/vnnrHXr1kyhULC2bduy9evXs7lz52pNxzZp0iS2ceNGfv7OnTvzaZc0t09s6ip9ZeQkJSUxAGzz5s1694k23bt3ZwDYP//5zzqf7dixg/Xv35/5+/szFxcX9thjj7HXX3+d5eTk1CmzkHO/pqaGzZ49mwUGBjI3NzfWp08fdvHiRda4cWM2YcIEfj7N1FWcrVu3ss6dOzOFQsF8fX3ZqFGj1NJEMfboPGzYsCG7du0a69+/P3N3d2cBAQFs7ty5ddJACT3GYlJXLV26lL3//vusefPmTKFQsJ49e7Jff/21zvzXrl1jY8aMYYGBgczZ2Zk1bdqUPfPMM2zHjh38PBUVFWz69OksKCiIubm5se7du7OMjIw65w1juq9LxhiLiYlho0eP1lt+YnsyxizcG4EQUsfQoUNNSh1EbOuPP/5A27ZtsW/fPv7VKZGmqVOn4vPPP0dubm6dAUKkrqioCI0aNcK7776Lf//737YujsM5e/YsunTpgszMTJ2dZYk0UJtVQiysvLxc7e8rV65g7969asPQEvvSokULvPrqq1i0aJGti0L0qKiowMaNG/HCCy9IPlDVvE8A4DuE0b3CMhYtWoQXX3yRAlU7QDWrhFhYUFAQxo0bhxYtWuDGjRtYu3YtKisrcebMGUqXQogF5OfnIz09HTt27MCuXbvsouZsw4YN2LBhAwYOHAgPDw8cPXoUX3/9Nfr3748ff/zR1sUjxKaogxUhFpaQkICvv/4aubm5UCgUiI2NxX/+8x8KVAmxkKysLIwaNQr+/v746KOPJB+oAo/yJTdo0ABLlixBSUkJ3+nq3XfftXXRCLE5qlklhBBCCCGSRW1WCSGEEEKIZFGwSgghhBBCJMvh2qwqlUrcuXMHnp6eVh3+kBBCCCGECMMYQ2lpKYKDg9WG2NbG4YLVO3fuoHnz5rYuBiGEEEIIMeDWrVto1qyZ3nkcLljlhmS8desWvLy8bFwaQgghhBCiqaSkBM2bNxc0lLbDBavcq38vLy+rBqu1SoYT2YXIL62Av6crosN84SSnZgiEEEIIIboIabLpcMGqLaSez8G8PVnIKa7gpwV5u2Lu4HAkRATZsGSEEEIIIfaNsgGYKPV8DiZuzFQLVAEgt7gCEzdmIvV8jo1KRgghhBBi/yhYNUGtkmHenixoG1WBmzZvTxZqlTTuAiGEEEKIMShYNcGJ7MI6NaqqGICc4gqcyC60XqEIIYQQQhwIBasmyC/VHagaMx8hhBBCCFFHwaoJ/D1dzTofIYQQQghRR8GqCaLDfBHk7Qp9SRcCvRSIDvO1WpkIIYQQQhwJBasmcJLLMHdwOADoDFgrapRIy8q1XqEIIYQQQhwIBasmSogIwtrRXeDlpj1lbXFZNaWwIoQQQggxEgWrZtAvPBCuzk5aP6MUVoQQQgghxqNg1QxOZBcir6RS5+eUwooQQgghxDgUrJoBpbAihBBCCLEMClbNgFJYEUIIIYRYBgWrZhAd5osAL4XOz2UAgrxdKYUVIYQQQohIFKyaqFbJcCK7EPHtA7V+zqW0mjs4HE5yfRlZCSGEEEKIJqsEq6tXr0ZoaChcXV0RExODEydO6Jz3008/Rc+ePdGoUSM0atQIcXFxeue3pdTzOeix+CBGfvoL/ptxQ+s8gd6uWDu6CxIigqxcOkIIIYQQ+2fxYHXr1q2YNm0a5s6di8zMTERGRiI+Ph75+fla5z98+DBGjhyJQ4cOISMjA82bN0f//v1x+/ZtSxdVlNTzOZi4MRM5xbo7TW1I7IajM/tQoEoIIYQQYiQZY8yiyT9jYmLQrVs3rFq1CgCgVCrRvHlzvPHGG3j77bcNfr+2thaNGjXCqlWrMGbMGIPzl5SUwNvbG8XFxfDy8jK5/FrLpGTosfig3kAVAH5L6Q9PV2eLlIEQQgghxF6JidcsWrNaVVWF06dPIy4u7u8VyuWIi4tDRkaGoGWUlZWhuroavr7aOydVVlaipKRE7Z+lncguNBioAsDpG/d1flarZMi4VoDdZ28j41oBDRhACCGEEKKF9jFCzeTevXuora1FQECA2vSAgABcunRJ0DJmzpyJ4OBgtYBX1cKFCzFv3jyTyyqG0Hypd0u1DxSQej4H8/ZkqQW8Qd6umDs4nJoMEEIIIYSokHQ2gEWLFmHLli349ttv4eqqPUfprFmzUFxczP+7deuWxcslNF9qE8+66az2nsvBBC1tXXOLKzBxYyZSz+eYpYyEEEIIIY7AosGqn58fnJyckJeXpzY9Ly8PgYHaUz1xli1bhkWLFmH//v3o2LGjzvkUCgW8vLzU/lladJgvgrxdYSgRVVRII7W/9567g+SvM7XOyzUCmLcni5oEEEIIIYT8xaLBqouLC6KionDgwAF+mlKpxIEDBxAbG6vze0uWLMGCBQuQmpqKrl27WrKIRnGSyzB3cDgA6A1YVfOqpp7PQdLmM9AXhzIAOcUVOJFdaJ6CEkIIIYTYOYs3A5g2bRo+/fRTfPnll7h48SImTpyIhw8fIjExEQAwZswYzJo1i59/8eLFmD17Nr744guEhoYiNzcXubm5ePDggaWLKkpCRBDWju4Cfz0jV3FqlQzz9mQJXrbQNrGEEEIIIY7Ooh2sAGD48OG4e/cu5syZg9zcXHTq1Ampqal8p6ubN29CLv87Zl67di2qqqrw4osvqi1n7ty5SElJsXRxRUmICEL3Vn7okLJf73xCswdwhLaJJYQQQghxdBYPVgEgOTkZycnJWj87fPiw2t/Xr1+3fIHMSN8QqieyC3H/YRV++v2u4OUFebsiOkx7mi5CCCGEkPrGKsFqfTVu/UnR35k7OFxvAEwIIYQQUp9IOnVVfSKXAWte7kJ5VgkhhBBCVFCwKhGrRnbGwI4UqBJCCCGEqKJgVQJWjOiEgR2DbV0MQgghhBDJoWBVAvq09bd1EQghhBBCJImCVQk4feO+rYtACCGEECJJFKxKwN3SSlsXgRBCCCFEkihYlYAmnoZHwSKEEEIIqY8oz6oERIU0Qq2S4UR2IfJLK+Dv+WhgAMq3SgghhJD6joJVCTh4KR+L9l1SG5I1yNsVcweHU95VQgghhNRr1AxAAiZvOasWqAJAbnEFJm7MROr5HBuVihBCCCHE9ihYlSj213/n7clCrZLpnZcQQgghxFFRsCphDEBOcQVOZBfauiiEEEIIITZBwaodyC+tMDwTIYQQQogDomDVDvh7utq6CIQQQgghNkHBqsQFeT9KY0UIIYQQUh9RsCpxswe1o3yrhBBCCKm3KFiVuAU/XKT0VYQQQgiptyhYlTjKt0oIIYSQ+oyCVYmjfKuEEEIIqc8oWLUDlG+VEEIIIfUVBat2hPKtEkIIIaS+oWDVjlC+VUIIIYTUNw1sXQBimAxAIOVbJYQQQkg9RDWrEsdlWJ07OJzyrRJCCCGk3qGaVYkL9HbF3MHhSIgIsnVRCCGEEEKsjoJVCduQ2A09WzehGlVCCCGE1FvUDEDCosN8KVAlhBBCSL1GwSohhBBCCJEsClYJIYQQQohkUbBKCCGEEEIki4JVQgghhBAiWRYPVlevXo3Q0FC4uroiJiYGJ06c0DnvhQsX8MILLyA0NBQymQzLly+3dPEIIYQQQoiEWTRY3bp1K6ZNm4a5c+ciMzMTkZGRiI+PR35+vtb5y8rK0KJFCyxatAiBgYGWLBohhBBCCLEDFg1WP/jgA4wfPx6JiYkIDw/HunXr4O7uji+++ELr/N26dcPSpUsxYsQIKBQKSxaNEEIIIYTYAYsFq1VVVTh9+jTi4uL+Xplcjri4OGRkZJhtPZWVlSgpKVH7RwghhBBCHIPFgtV79+6htrYWAQEBatMDAgKQm5trtvUsXLgQ3t7e/L/mzZubbdmEEEIIIcS27D4bwKxZs1BcXMz/u3Xrlq2LRAghhBBCzKSBpRbs5+cHJycn5OXlqU3Py8sza+cphUJB7VsJIYQQQhyUxWpWXVxcEBUVhQMHDvDTlEolDhw4gNjYWEutlhBCCCGEAKhVMmRcK8Dus7eRca0AtUpm6yIZxWI1qwAwbdo0jB07Fl27dkV0dDSWL1+Ohw8fIjExEQAwZswYNG3aFAsXLgTwqFNWVlYW//+3b9/G2bNn4eHhgVatWlmyqIQQQgghDiP1fA7m7clCTnEFPy3I2xVzB4cjISLIhiUTz6LB6vDhw3H37l3MmTMHubm56NSpE1JTU/lOVzdv3oRc/nfl7p07d9C5c2f+72XLlmHZsmXo3bs3Dh8+bMmiEkIIIYQ4hNTzOZi4MROa9ai5xRWYuDETa0d3sauAVcYYs886YR1KSkrg7e2N4uJieHl5WXx9ZVU1CJ/zo0WWnTU/Hu4uFn2eIIQQQogDqVUy9Fh8UK1GVZUMQKC3K47O7AMnucy6hVMhJl6z+2wAhBBCCCHkkRPZhToDVQBgAHKKK3Aiu9B6hTIRBauEEEIIIQ4iv1R3oGrMfFJAwSohhBBCiIPw93QVNN//fr9nN9kBKFglhBBCCHEQ0WG+CPJ2haHWqDsy/0TUu2lIPZ9jlXKZgoJVQgghhBAH4SSXYe7g8DqZALQpKqvGxI2Zkg9YKVglhBBCiN1xlIT3lpAQEYSpca0FzcsAzNuTJen9R3mRCCGEEGJXHCnhvaWE+jUUPC+XHSC2ZWMLlsh4VLNKCCGEELvBJbzXTM/EJbyX+ittaxHa0Yoj5ewAFKwSQgghxC7UKhnm7cnS2h6Tmyb1V9rWEh3miwAvheD5xQa31kTBKiGEEELsgjkT3jt6m1cnuQzvDGwnaF65DLj/sMrCJTIetVklxEHUKhlOZBciv7QC/p6uiA7ztelQeoQQYm7mSnhfX9q89gsPEDSfkgGTNmdirbyLJLefglVCHIAxN14KbgkxjK4TaRH6qlrffFybV816VK7N69rR0gzYrGHeniz0Cw+U3DlOwSohds6YG299qVUg0mUPQSBdJ9LDJbzPLa7Q2m5VBiDQ+9H5pI2hNq8ySDdgszTVJhRSywpAbVYJsWOGbrwMQMp3F9TaYlFPWmJrqedz0GPxQYz89BdM3nIWIz/9BT0WH5TUuUfXiTRxCe8B1Bmhift77uBwnYGmOdu8OiopZgWgYJUQO2boxgsAuSWVWHXwKgDqSWtvHLEDiFSCQH37Vsx1IuVjJOWyCaGr/AkRQVg7ugv8NXq6B3q7GnyFb642r1LH7bsfzom/nqSYFYCaARBix9KzcgXN92H672gT6AFvNxfBtQpSew1U3zjiK2ipvII1tG+F1r6tOngVW07elOQxsvfzx1D5EyKC0L2VHzqk7AcAbEjshp6tmxg8b8zR5lXqtO07IQw1obAlqlklxE6lns/B58euC55/3p4s5JbUj1oFeyeV2kdzE/MK1lK1gkL2rdDz/8P03yV5jOz9/NFV/pziCkzYmIkFey4g41qB2mdC2zzff1gJfbPJ8CgolmLAJoSufWeIkCYUtkQ1q4TYIa6GSoyc4goUPqgUNK891yrYO6nUPlqC0CAwLSsX07adNXutoNB9u+ylSKPXYetjZO/nj77ycz4/dh2fH7sOLzf1EMZQp73U8zmYtPmM3mUD6gGbPXQE5AjZd7oESrzWnYJVYhfs6YZhDULaqmrj29DFpJ60xPLE1D7aW1MNoQ9BX2h5Y2COtEJC9y0Y9F4nhlj7GHH3x9zicmTevG/188ec92cx97aS8hr+/6dvO4uMa4UoKq/mp6k+4AgJ5OQyYNXIv88vWzalMGafGvu7ILQJhS1RsEokz97bXlmC0Laqmm4WlmPu4HBM3JgJGaB245b6a6D6wpE7gBhKOwQ8Chi0vfE3R62g0H1272Gl3utEaABrjWNkbPtEfWUTEyiZ+/5s7D7bdz6vzjTVBxxD7fWBR+ddo4YuAGybi9XYfWrsvrOHyh9qs0okzZptr8S2kbNUmzpDyxXbVlXV8vTfAcDonrSOTgq9px25A4i+tEMcfbtcSFohfcdQzL7V1+N8alxrwcuxJGPbJwLAvdJKree3mLRilrg/m3OfqWZuENNe35SsKcbeQ7jvzd9zAROM3Kf2eE8QimpWiWRZs+2V2CdZS9X2GlquMW1VNc3bk4WjM/sY1ZPWkUmlBt/UpOfmYqmmN1wQOPe7C8grEdaGWpOuGiRDx1BIzW6gl4Lft7p6nAPAlpO3BC/HEkxpnwgAC364iM+OZqud32JqE815f1Y91/waKhDo5Yq8EuOaYGgri9j2+sY2xTH2HiKkdpzbpynfXYCnqzPySypQ+LAKvh6P9ld0mK+g89teUbBKJMtabffEvu6x1OshIcsV8ipLH9V9Ftncm5/O/ahmXCtAbnF5nZugowexUhp+kat9tGVTDUsH7gkRQSivqsXUbb8a9X1tNUhCj6GufcupqFEiLSuX307V/ax6LYhdjrkZ2z5Rleq+6RceaLA2MeW7C3zwaY77c62SYdXBq1h/LFutramPWwOzB1ti2ut/f+6OoGWqPjTpOv+4DAZrXu6MgR2D6yxD1/e0YXiUN3vUZ8frfObj5ozE7qGYPSgckzbrPi+1OZFdiPsPqyR936dmAESyrNF2T+zrHksl1Re6XKGvsgzR3GdpWXn8q7+p237Fgh8uYupWaY4sZG5SHCjBlKTnprJG05taJcOi1Euiv6ctrVCtkuHYlXt4e+dvgo4ht281e5JzisuqBW2nuZZjLHO0h+VGuXt752/44mi2qAFGTL0/p57PQdS7afgw/Xe1QBUAilQ6TplLoLebwSYo3EOg2KY4Qmq5k78+g70aCfpNrR1XVVRejQ/Tr+CdXb/htV5hde4d+oxbf1Ly930KVolkWaPtntih9yw1VJ/Q5Qp9lWWI5j6bvOWszvXn2EluRmNJdfjFhIggpE/rzf+9IbEbjs7sY9FA1VqB+4nsQtFNALTVKnPtK0d9frxOwKNK8xj2Cw+Eq7OTznm1DVOsjaHlAJZ70DFn+8Si8mq8t/eioHk/TP8dqedzTLo/p57PwYSNmSgq033MzIl7wNH1EAgAy0d04q8t7nW6rqBW86FJSC23kgFJmzOx99wdZFwrwLeZfyLlu/Mm145rKiqrxic/ZePthLYmLUdq930KVh2EZqPuqhqlzTuKmErsDcMYYmsHLFXbK3R+7lWWKS9nfNycoWRM9DnhqMOwCt33x67etfr1pPka7vu/fugstX5rBe7G1Apq1iob07mIW6+QYFm1FlEXQ8ux5IMOd3+0hXl7shAV0giBXrrXr+v+bEq7e2PuezKoP+BoPgRy+oUH8P+vryOgtocmMedz8tdn+DdYX/1yU/D3xGAAUkzs28CRyn2f2qw6AG3tyzTTv9hjqidrtN0TWztgqdpeofNzr7L0tZUzpKi8GqM+Ow5vHa8vtbHn3J6GCN33qw5d4//fFtfTuPUnLb5+oSnRVH+cjemIJfb6mJnQBq/1aqmWqN2Y16fcesWMUPWYr5vWz2qVDMeu3hO0HK6HuTk7rKneH60dSuQUV2Dt4WuoqKnV+rm++7MpbW3FbqePWwMkdg9D5V+VN9w+F7LfdXUE1JY8X8z5bK24r1jPmwahpHTfp2DVzulqnK15Qdiio4g5iLlhGENsz2tL9dQWs1wnuczk3tQAUGxEuzB7zO1piDE9aK11PaVl1c0daan1i0mJxqU9SsvKNaojVnSYLwK8FILP39FPhKgFGGIDHs3rUkxwMX173U5gYnObXr9Xhh6LD5q9w5qhzAqBXq58O3djH251+fCvNHjauDnLMbBDEMqr1YNEwHr3kAERAThzsxgfpl/hp3H7vNfjTQQtQ1c2CM1gl7uHmPuVvlRI4b5PzQDsmJjaBVt1FDEHS7bdE/u6R+z8liqHrldZlnb9XpnV12lpQnJ/arL09cR1Gpqz+4JV1i/21eyCHy4i6t00o/NBOslleGdgO8HrO3gpX+1vMT+e2q4fLlgWQnP37j2Xo3W7da3bx90Zy9N/t1iHNV33gg2J3ZA2rRf/tzXv+mXVSuzIvK21s461coGmns+r0yGV2+e6HgK1EdIUR/Ue4oikkL+VglU7JrZ2wVYdRcxBV/oYcxDb89pSPbXFLtcWaUWW/9W5QixTE2Wb2lbU0HIM9ezWxlLXk2qnIX2v8sy5fmNezerqHMN1UJq58xyOXb2n85ipthE0ZMqWs2rnnZgfT23Xj9hgmbPvXA6Sv84UNK9qTaalO6xpuxdIKfWQamcdQ30RdJEBetvIatK3zxfuFZ+JAnjUFEdzoATu3lJZo8Tkvq0hkV1uFuboF2Iu1AzAjhlbNa/te5ZKAG5p5iq30Nc9uuafmdAGgd5u8HZzQa2SiSqD5jb8OKUXOs1PE1QOWxA7EIM5E2Ub+p6280Hbq2ouJ2Fyn9b8dvQLD8Tc7y6Ibh5hzldkYnIumnP9lnjNV1xeg1GfHTdb+1rV805o040vxnZF7zb+Ws9VMcEyZ/qOc4LnDfR2xYhuzdVeQ2uydJtAMTWIQpjalIA7hrMHhSNps7CgX9WsgW0xectZE0rA5SoVfr7ra4ozYWMmfNyd1R7cvFwboKTC/Km3bEUqw29bJVhdvXo1li5ditzcXERGRmLlypWIjo7WOf/27dsxe/ZsXL9+Ha1bt8bixYsxcOBAaxTVrhhbNc+1N1NN/yKFkXvEMne5xdbeqn6+OPWy6DLoSoit+npSag8NYn9cjU22b8z3tJ0Pmj8kHC4n4fqfr2PR8x2QEBFkVDolwHyvyEztNGQKS77myzFD+1ruvNtwLBvjuocZ7HzJ/f1Ey8Y2uX4GdwxCn3YBuJb/QND8mg8L5ngIr1Uy/EdgOiqhhATgunDH8B8bTuDsrWKd8w2ICEDmzaI61+LyEZ2MesAQiqshVd3nAHTuQ+4c07y/lDpIoCq1GMDiwerWrVsxbdo0rFu3DjExMVi+fDni4+Nx+fJl+Pv715n/559/xsiRI7Fw4UI888wz2Lx5M4YOHYrMzExERERYurh2xdhG3arD7AGQzMg9YkhhxCFjO79wQerHP11DWVXd3rT5KjfpE9mFKC6vllxtt5CaOGOHYxSS61N1JB1A9/lgKI9j0V+J21e/3BlZOaUGt0mTXAbcf1gl+nvamNppSJXYYCc6zFetM44lmGNoZM0hQrV1LgrQsR3cPuFGaPNwtdzP355zOdhzTnhzGdWHBXM9hJ++cd+kDpiaJvRqgbf+yt1paKhZfY78rj+DwvvDOqFWyfi3VhxLBqoAEPfBEbX9FfRXYC52H9pXjxDtpPhGz+LB6gcffIDx48cjMTERALBu3Tr88MMP+OKLL/D222/XmX/FihVISEjAW2+9BQBYsGAB0tLSsGrVKqxbt87SxbUrXO3ChI3iX6eovsIQG0xYkpAfWXOOSW1KOfU9cesqQ+r5HLz9zW96gyjV7bJGuiJjCKmJM3Y4RiFBG5cDM7lPK/xyrUDn6EVCMDzKfWhMs0ElAyZtzsRauekPR0LTRgG6O/Ppqq3Xd+5w3ympsFyCdnO+7tZ8GNRsvhMV0qhOoCO29761aD5wmPMh/G6p+QJV4FEtNfDo+hwQEYgvjl03e4YBAKKbUZmLZlCaW1xhVA2yI5BSxQjHosFqVVUVTp8+jVmzZvHT5HI54uLikJGRofU7GRkZmDZtmtq0+Ph47Nq1y5JFNRpjDIoa894UOMqyMihr9B+i/i28MaNnc6w8pD95tS7lJZXQ1ye2sKASJy7eRowVGlinZeXiP3svqdWKBHq54p2BbdG9lR+/n09fvoPCgmKzl1tZVcOvw9C+P5FdiKLCElFlSMvK5dtbCR8I72/3CyoxZUMGVozohH7hgWrlNWRwxyBRNT26yPCo5qprgCuUZfozA9y9e19Q+e7evQ9l0N+5LA+dvS7oe2tSf8Pmny7xgb8x+9QcZAAWfnMGfUO9jL7B/3g+FxuPXBa8Ddx10a+FN38c0rJyMee7C1r3h+a5w9H3HUtQPdZizl9NqvtcVlPLL6drgCtQUa52Hf947h6mbDkLBtudI9pwZ0pKv3aQVZSjWsmw8JszcNGxT/SdZ9r2ZRNnpVl/m4oKitFn6yn+/qxA3Xze5nD68h10DWlUp+zceW6p31vyN8akVz8sYxYs1Z07d9C0aVP8/PPPiI2N5af/61//wpEjR3D8+PE633FxccGXX36JkSNH8tPWrFmDefPmIS+v7mvXyspKVFb+ffKWlJSgefPmKC4uhpeXl5m3qK4HRSW49USMxddDCCGEEGJpzX85Dg8fy8dPJSUl8Pb2FhSv2X3qqoULF8Lb25v/17x5c6uu392FEioQQgghxDFIMa6xaIn8/Pzg5ORUp0Y0Ly8PgYGBWr8TGBgoav5Zs2apNRvgalatRebmhjaZp822PO5VsTHV3Utf7Ij39/9u0U4StvTxK1HoGtIIUe+mW2wdE3q1wD96hiH6vQOC5rdEmy0xTv9fnNYbS62SIe6DI4LOBf6Vcrj6NVarZDh14z5u3y/Dv3edF1WuozOfRo/Fh0R9h+g3JjYEk/u2tuj5b8jHr0ThyO93sfm4ecY013b+1ioZNv5yA4tSxeXC9PdUIP+vdpqn/y8OAETvq8XPd8CH6VeQVyK881CglyteimomqCnWx69EoVfrJiirqjH7cVTdl6rL1/a7EODpijyB6com9GqBdT/9YXC+DYnR6BrSCKdu3EdecTnO/lmEr0/cMmJL/sY1dljyYke89VfaMF33vOPZhRi3/oRJ69NlcMcgnMi+L3ifyQCsGNEJp2/cx5cZN0Svj2tulT6tt9oQw6du3Mfd0go08XTl97WQbd6QGC2qSZzMTfsQw7Zk0WDVxcUFUVFROHDgAIYOHQoAUCqVOHDgAJKTk7V+JzY2FgcOHMCUKVP4aWlpaWrNCFQpFAooFLZriSSTySBzdzfLsmqVDClp2ahoYNz2MFc3HPy/AdhwLBsLfjBvyhIpuFsth9zdHZVG7h8hVvx8G1+eyRe8DlsHq3J3d8i13bivFeBGGQMMbMfsQe34VEB1lg0gtn1DlFZUY8b34joayNzcLHqc6qM+nUItfv4bEtUmGDN/uGK2Mmg7f+UAxvRph09O5ojqcf5nOcD+Kpf8r3uy2HIGBPhi1vOdMfGvTqu6UmJtSOymlqUDADadu2uwvFFtgiF3dYa8QY3Zj6PqvmTyan75Xo28kPZOAk7fuM93XI1o6lWnE5outQph13J+tQzOHg0R274hAMDl3B3890y+yW1aZQDeO3CdL8OpvArtQ562c4NvY2+jsxTosyOrEIFeCvTrEoj/Xbmnd7AO4FGgGt+pKbwaFeCTk8I7TgJ/B+iznu8MZ4+G/HTufqwqup273m3mO/C1awq5xDpMiWXxZgDTpk3Dp59+ii+//BIXL17ExIkT8fDhQz47wJgxY9Q6YE2ePBmpqal4//33cenSJaSkpODUqVM6g1tHYswoMqqmbDmLtKxcjOseZtQIIVLXxNM6P9JFBm5Eqj4YFinJEUuEJnn381To7RSUej4HcR8cMUuZNiR2w9IXO5plWfWNLUeRUR01SGg6pFeeeAyLX+hgcL4T2YVaR28yZvhK1aWIHRFKdaQeLiVWoLd6xosAlf0QHeaLIZ2aIvavPK5Ch+y1ZA9rbps1r9lx60+i99JDKC6vUiszR1eJuH0iNHuD6lDMqedzMGmzcRk2NDEAeSqZDcatP6k2fCvHmGGTVRn6Tm5JJb4/l6MzUB0TG8L/P5dmS8hoXZqnhJiREC01/LcUWTxYHT58OJYtW4Y5c+agU6dOOHv2LFJTUxEQ8Ohg3rx5Ezk5f590Tz75JDZv3oxPPvkEkZGR2LFjB3bt2lUvcqyaYxQZbnxvYy7aCb1bmLx+S4oKaaR1+obEblgxohM2JHazcomARg1dzN4b1hyEJnnXNx+XRsdcuRqjw3wxqKPt024B6j8s9sBWPzgyPBo1iCM0HVLXUF8Mjgw2OJ+uwAN4NErc6pc7G/UwuPGXG6IDVtV9nBARhKMz++Dr8U9gxYhO+Hr8E0ib1kvv93UFudZy+sZ9ndcsl/pK11DJ+gKd4rIqQceAG4pZzOAWPm7OGBChvYmfPrq2x5hjIAPweq8wk46bDMD+C3U7gBsKJmUAVo3srHaeHZ3ZR1QaPF3bbOrw31JjlVa0ycnJOmtGDx8+XGfaSy+9hJdeesnCpZIeU0eRUc1lyJ3AmrkFvd2cdT4ZNpD405euH+voMF+4uzRAWZX1Rw755Y8Cq69T1YnsQu2vxAwMR6kvoTxg/GhK9kCGR6mi7IFcBqwa+fcPjqljyIvRyN0ZC5/vgF6PN+GnCX27IeZepi+HaKOGCqMeBhenXsaGn68LmldXDlonuUytVlHI/SUhIgj9wgNxIrsQf94v49tZarLEccwrqcDy9CuC8k+rWj6iExbtu6T2OxH41z4BgEmbzwi+D8zbkwVPV2dBbwi5JkgnsguxT+T1qC+XNXcMhDaHmxL3OCbHtUZks0ZGDQHLlUdX/wBdv8WBZsybrXre2duQ6UJJr8tXPSZ0vGtDuBpabSewvrZKqw5dM2Gtlnciu1Bn7aqtrDtiuOOBNlybt3ullSa1Lx63/qTWH1tDw1EC+mvrTG2SoouxP9KBXgqMjH4Mgd6umLnzN5PKoPlakZP8dEu0DvDElbwHWGVk3mJzWzWyMwb+VRudej4Hc7+7IOr7bs5OKK+uO0qaPj5uzkjsHorkPq3hJJepBWlRIY0Q4KXQWduu+hD03dnbgtanL/Aw5W1TvpYybkjshvsPq1D4sAq+HgoEepn/R50LcksrvLQGq6nnczBzp/Yg1hRnbxYJHoQjsrk3P71feACe6RhcJ9ABgB6LDwr+LeKWn3FN2AM81wTJ2N89fQNNOMll8BP4YBXq545aJcOCH7JErF0/zUoEawSTmg9XjoaCVQnRF2CIoVqrYUztgBS4OzuhTONHdtz6k/D3dNH5HWvWOpmKqw2uVTJ8djRbbw2oruEjObpqpkx5oje1SUr6xbqvxNKy8owaq5yrgUnLyhUdrImx+tA1rB3dBd1b+dk8WNV8ANE1spE2qp1/yqtq8I8vTxn8zsyENgj2cdP6I6p6XZ2+cR9vJ7TF1G2/1lmG6kPQj+dzMH173Xl00RV4mPK2Sdu+4q47S9P1YLHsx8v44th1i6zz65PCet4/ura91aZpC3QyrhUY+cAq7D7MHVtTf/d03avENIUy98O5tkoERw8mLY2CVYnRFWAIFeil0NsRw14COs1AlZNfWncc9rSsPCgayC1SW2FpQmpAZw1sy498pY2QV2Jin+hNbZLy9s66qa70bYM2XC0dF6gKDdaMxQCkfHcBP/2rj8EmFIYeIIyh2cNcNWWNmCYZqgFZqcAhVEc/EQJPV+c60zWDLn0PjKqvjpM2nxFYWnWagYe53jZZk74HC0sFqmIIvbaNfWCNbeGHnZm3RTVBMuV3T9f2iGkK9f25O6LWKYQxQ+QS3ex+UABHpNq4X2zv6YoaJdJ0jDFuzp7dUjJ5y1lM2JiJ4nLz1Bpbu5WPoQbyfdr6G1yGas2UJu6JXltPYF2E9GK1NIa/Owpaq/1sbkkl1h6+ZrCHrWqnI3PR7GHOMaXWR0ivb835OLo662g+MPq4OWNqXGscndkH/cID+U6extAMPEzt4a1JV/YBIVS/p2s55m7r7eNe9wHCWKoZD4Rsi9gHVm75T7RsbFQPdc1ObZtejUGgl+57kOr2aCOmp7ypD+facHt03p4su6kkkjIKViWKCzDE9p4uLqvW2kvS3D27HUXy0y3h46b+g+BvQoqs5KdbCnrA0PyB0Nb7+OjMPgAg6gHDHBklAPMHCcYYExuChIggi7Wf1eXD9N8BQO8DBJeaxhrMdUz9vYSf12KCruLyaixPv4K0rFyjj5W+wENfD++hnQxnHFClL/uAPtrSQWlbjtnPVTONhq4anKVl5QraFjEPrJrBn7E91FUfrLu39kPKs6alZRJaDqEppp5u00TPHHXpq0Qg4lAzADv0xdiumLrtV629+rW9Enbknt2mSnq6Fab2a8O/Jr9+rwybT4gfcYTTvVUTRDb31tkLmCOkTZOYdoocc9YQmNokxVRcjbIpwZpcBqN6k8/bk8XXFGprQmGJtt+6MjuY65imT+uN87dL9PZSVy2L0GOues/5V4L4GmchgYeuXvYLhkbgeHahqGYCOSJfz+q6DrW95jXXgwWnyExvi1SbaAjdFjFtSbW1gzdHpyJz9KQXUg4h27pqZGc0aqjAoct3BZefY+7zoj6iYNUOyeUyvSNoaHZWsHbNlL3hgsTU8zlYnv67UUG9avunyhphva/1tWkS+4BhKBWVsYSm4rEEpZKhVsmMCta4nyHuBya/tAJ+DRWYvv1XQW1NVa8fa3WK0JXZwVztNrnzvKzK8MOU2B9X7p5T+ED8mxuhgYe28pvSOUdbG29N+q5DbRUDlnidbKzZg9rBz1MhqHe/rnbvOoPFv7JzhPo11BuEmqNTkTmCXiHl0LWtqtdkrZIZdS1K6bywVxSs2iGhibm5HxxTnuoCvRToGuoraIg5e8S9ijel9tnYkUL0dYwS84Bh6ZFKxAQ55vSPL08hyNsVswe1E/0DoSsASnk2HBM2CsulaIvaEEM1XNZi7I+rb0MXQcdqaGQwerf1N0vqKGPeAOhLe6TK0HWouRyhDxamZHsxRLVjoup+NdS7X9c+kUIOT2v1pDe0rWIfjixViVAfUZtVOyQ2MbexPzyzB7XDnGfa4wc9Q8zZu9M37gMwra2ZKSOF6GrTJCZQcrSRSlTlFldg0uYzeDZS2La92j1U7ygwCRFBmBrXWtCybFEboqtThpjRnMzRmcPYDnaB3m6C2jrv+vUOlqReQnF5lVmCHq7Nd/LTrUR9z9B1JvQ65OYz1NbbHKMl6aPvwVXstqgyppOmvTK0rUJHybJ0JUJ9Q8GqHYoKaaT3h0SzswL3wyOWb0MXLPjBNm1dVccjF8qY2wFXSy22Fo0b4tWY4fG00Vy/0EBp9qB2Zlm/VHHn3ne/5mDlCP3BmgzA3vO5Bmt9kvu01nt+GeplbGm6HmCEjubEPYCZQmwHO9V9JvTH3NAQoGI5yWXo3spP1HcMXWfGDFusb/uXj+iEWQPD+c6U/+geCkB3ByIfd2dR9zV9D67mGIKZPKLZIXZq3ON17imOXIlgC9QMwA6JHZ1IdX4xgWfhwyqzt3V1d5ajrFqp8/NXu4ciLjxQbaQtLv/k9Xtl+PrEzTptDleM6ARFA7lRHYG4WmqxN2hzJxfXXL/QHIGar/qkwlDHJjGvQbngLa+0Qu8yhb7adZLLkPLs36/VxY7uBTzK7Wtpmg8wQh+oNh+/qbWjllhCX69r22fc69RfrhVg0uZMFAnsDGoq7roxdB8Q+nrW2GGLdbX15rJIcLV3XNMBXR2IAAh65aw56pg5t4Vop9k0IblPK4ce7tTWKFi1U2J7SYpp18XdtHw9jE/hpG2ZAPDB8E4AoLcRO6A+0pZqYMjdEDR/ANxdGqi1NeI60+SV6G87xg3fasvk49pq8YS0jZLa6yXVchqqAWT4uwOI0OFNbxSWCSqHkKDOlF7GtUpm1EhcYmk+wAh9oNqflYeod9Ow6PkOajWNXLYBMTTb8Gl7YNS1z5zkMsjlMq2BKkfoA4ZQYh7MhVw/pgxbLLStt6F2kkLu21z6sDaBnnrTQpkyBDPRj0aosiwKVu2Y2IbvqvPvv5CD9T/XTdGketPydtM9tKku3E3Qx90ZRWV//0hp/qAZ22Bf27jbqil/VG8WXO2ZvloJYxvOm5OuHwh9DxjLR3SS3Oslf08F8gR2/gMejQ0+pFNTZFwrEBSshvi6CyyHsKDO2I4jJ7ILBeUrTmjvj/1Z+UalztL2ACO01hAAisqqMWFjploOYS7bwNsDxKWXMqUGyZR2ksYy9GCuLeOCMcsTkz7JEH2Bjjlrqa2xLYRYAgWrdkg1OBP7NKf6+immRWO9Ny0uTYeYV+vc94XktdNXbs0RVlRfa2obAlLbD5DYXsLWzisqlwGrRupv02TodaKtqTbREJuflgsqhb6efCU2FJ8dzTbra0xjakOEBlapF/IxvmcYPv1ftqjlA9ofYLgHKqHZDADUCWxyiyswReSwt5rE7DNbtZNUvW5yi8tR+LAKvh4Ko7MP2LpHvDlrqW29LYQYg4JVO6QrOBNLTJoOQ5VDXFtT1e8b+0pEXzAKCE9qrbmNQvKEWjOv6KqRnTFQwAhltkodpUlbL/PoMF/89PtdUflpNYNKoa8nXRrIJfEaU0xg9f25HKwa0Rlvbj0jqIbV0ANMQkQQXu0eis+NHGOeq4GzFlu2kzT3a1lbv+Y1Zy21rbeFELEoG4DE/Xg+V+t0c/WkFZqmQ1c2gSBvV6wb3QWzB7c3S0oTXcPC5hZXYMLGTLz9zW86k1oD2sdhFjt0rZD5TRljnNtnAzuKGyrSljSHm+T8eD5X9OAFQN2gUuiwiMYO42hO0WG+CBA4dGlOcQUaeyqwamQXQfMLeYCJCw8UtCxdtB0rU85nfcSMz070o978pD6jmlUJ23vuDqZv/1XrZ5boSauLuV+p6WJotBgAau1gtc1jzs4a+uir3TbUU3z2oHZGBVX6mkZYkr5hX6du035+6qKvbZzQ15O2fo3pJJfhnYHtMFng6/T80goM6dQU6+TmaUcppu2qUOZ6W6MNtZM0D+rNT+ozClYlKvV8DpI2n9E7jzWDM2u8NjLXsLDWGnlIW9MDQz3FZQAW/HAR8RFBooIroe10zc2Ukb1UjYkNwYCIIINBpdDzzNavMfu09Rc87/V7j7IYmOuhz5i2q0Lk6Bn+11S2fsBwBNSbn9Rn1AxAgrgAQShbDAtpCebaDmu9BtPW9MBQT3FdCd/10dc0wpxJ1bUx1wPEgIgghxn5RleTCF2Wp//OHyMuyH6uSzO82rMFnuts3IhACRFBWDe6C3zcnQ3PLJK2pjTmUJ9GQbIUKTSDIcQWqGZVgsQGCI7SRsnU7bDFazDN2m1zp+ox1DTC0k1BTH2AcLRXk/qaROhjiWPE1VauOngV649lq/UU51LHiU3DZs23NcQ4VEtN6iMKViVITIBgy2EhzU1Imyxvd2cU/9VuVUqvwbhjZu5OEIYeXCwdXJjyAGHrY2JuxjaJsOQxcpLLMDmutdbcp2lZuUanYXOUtzWOytbNYAixNgpWJUhMgOAogQAgrE3Wouc7AKg7ApbQkYcsRTVnaICXQmdTALE1jbZIqq5KTGcezSFWHa0DjalNIiwZAGoLXlRr4I5dvYtVh64JXp6jvK0hhDgGClYlSMjQn0ISytsjoT2Hxb4G0+ygxEnLysOQTk21fkfI+O+awWdaVi4qqmt1zguIe8CwdboaMZ15lOzvIVQd8dWkqcGmLQJALoiNDvPFzszbgodadpS3NVJgqywehDgSClYlSMjQn0ITytsjIW2yxLwG09fOcPKWs3CWy+rkPBUz/jsXfBpqz+jj7oyFz3cQ9YAhhXQ1YhLRc0OoOiJjg00pBIBiBvhwpLc1tmarLB6EOBrKBiBRunp92mNCeWOYq+ewkHaGyV+fwd5z6j3qhY7/PiXucX5oWkPrUTSQo5/IhO5SSaouNBG9I78+5h4cxOxpKbXbFTLAB/UoNx9bZvEgxNFQzaqEUa9P0wlpZ6hkQNLmTKyT//1DLfSVb6ifu+D15JZUGtXJRgpJ1aVQw2trhtpUM/zdC58jtXa71hrgw5Ls4bW6rbN4EOJoKFiVOOr1aRox7QxVfzzEthW1dEcoWz+4UELyRww9ONjDw6U931Ps5bW6rbN4EOJoKFglDk3Ma2nVHw+xvfqt0RHK1kGGFGp4pcDQgwMFH5ahq024tpHkbM3WWTwIcTQUrBKHJnYcde7Hw0kuw9sJbTF126915tFWk1hfXpPbuoZXKmz94FDf2NtrdVtn8SDE0VAHK+LQVDsoCcH9eKSez8Gi1Eta59E2tKFUOkJZAw2bSaxNzGt1KTDUGU8GxxrQhRBLo2CVOLyEiCCsebkz9MVUqj8eunrxcmYPaqf1dSON202IZdjba/X69PBKiDVQMwBSLwzsGIxVkCFpc93k9qo/HgD0pqCSAVjww0XERwRp/aGh1+SEmJ89vlanNt6EmI/FgtXCwkK88cYb2LNnD+RyOV544QWsWLECHh4eOr/zySefYPPmzcjMzERpaSnu378PHx8fSxWR1DMDOwZhnVz/j0fGtQKTe/FSe0ZCzMte24TTwysh5mGxYHXUqFHIyclBWloaqqurkZiYiNdeew2bN2/W+Z2ysjIkJCQgISEBs2bNslTRSD1m6MfD3l43ElIf2HPqNHp4JcR0FglWL168iNTUVJw8eRJdu3YFAKxcuRIDBw7EsmXLEBysffSlKVOmAAAOHz5siWIRAkD/j4c9vm4kpD6g1+qE1F8WCVYzMjLg4+PDB6oAEBcXB7lcjuPHj+O5554z27oqKytRWfl3R5iSkhKzLZvUP/b6upGQ+oBeqxNSP1kkG0Bubi78/f3VpjVo0AC+vr7Izc0167oWLlwIb29v/l/z5s3NunxSv1AvXkKkjVKnEVL/iApW3377bchkMr3/Ll3SnpvSUmbNmoXi4mL+361bt6y6fuJ4KAUVIYQQIh2imgFMnz4d48aN0ztPixYtEBgYiPz8fLXpNTU1KCwsRGBgoOhC6qNQKKBQKMy6TELodSMhhBAiDaKC1SZNmqBJkyYG54uNjUVRURFOnz6NqKgoAMDBgwehVCoRExNjXEkFYuxRS0Nqu0rMoX0TZ7Rv4gwAePig1MalIYQQQhwDF6dxcZs+Fulg1a5dOyQkJGD8+PFYt24dqqurkZycjBEjRvCZAG7fvo2+ffviv//9L6KjowE8auuam5uLq1evAgB+++03eHp64rHHHoOvr7AOLaWljwIKartKCCGEECJtpaWl8Pb21juPxfKsbtq0CcnJyejbty8/KMBHH33Ef15dXY3Lly+jrKyMn7Zu3TrMmzeP/7tXr14AgPXr1xtsfsAJDg7GrVu34OnpCZnMOq9sS0pK0Lx5c9y6dQteXl5WWSfRj46J9NAxkR46JtJDx0R66JhYBmMMpaWlOtOZqpIxIfWvRK+SkhJ4e3ujuLiYTmSJoGMiPXRMpIeOifTQMZEeOia2Z5HUVYQQQgghhJgDBauEEEIIIUSyKFg1A4VCgblz51IKLQmhYyI9dEykh46J9NAxkR46JrZHbVYJIYQQQohkUc0qIYQQQgiRLApWCSGEEEKIZFGwSgghhBBCJIuCVUIIIYQQIlkUrJpo9erVCA0NhaurK2JiYnDixAlbF6neSElJgUwmU/vXtm1b/vOKigpMmjQJjRs3hoeHB1544QXk5eXZsMSO56effsLgwYMRHBwMmUyGXbt2qX3OGMOcOXMQFBQENzc3xMXF4cqVK2rzFBYWYtSoUfDy8oKPjw9effVVPHjwwIpb4VgMHZNx48bVuW4SEhLU5qFjYj4LFy5Et27d4OnpCX9/fwwdOhSXL19Wm0fIvermzZsYNGgQ3N3d4e/vj7feegs1NTXW3BSHIeSYPPXUU3WukwkTJqjNQ8fEeihYNcHWrVsxbdo0zJ07F5mZmYiMjER8fDzy8/NtXbR6o3379sjJyeH/HT16lP9s6tSp2LNnD7Zv344jR47gzp07eP75521YWsfz8OFDREZGYvXq1Vo/X7JkCT766COsW7cOx48fR8OGDREfH4+Kigp+nlGjRuHChQtIS0vD999/j59++gmvvfaatTbB4Rg6JgCQkJCgdt18/fXXap/TMTGfI0eOYNKkSfjll1+QlpaG6upq9O/fHw8fPuTnMXSvqq2txaBBg1BVVYWff/4ZX375JTZs2IA5c+bYYpPsnpBjAgDjx49Xu06WLFnCf0bHxMoYMVp0dDSbNGkS/3dtbS0LDg5mCxcutGGp6o+5c+eyyMhIrZ8VFRUxZ2dntn37dn7axYsXGQCWkZFhpRLWLwDYt99+y/+tVCpZYGAgW7p0KT+tqKiIKRQK9vXXXzPGGMvKymIA2MmTJ/l59u3bx2QyGbt9+7bVyu6oNI8JY4yNHTuWDRkyROd36JhYVn5+PgPAjhw5whgTdq/au3cvk8vlLDc3l59n7dq1zMvLi1VWVlp3AxyQ5jFhjLHevXuzyZMn6/wOHRProppVI1VVVeH06dOIi4vjp8nlcsTFxSEjI8OGJatfrly5guDgYLRo0QKjRo3CzZs3AQCnT59GdXW12vFp27YtHnvsMTo+VpKdnY3c3Fy1Y+Dt7Y2YmBj+GGRkZMDHxwddu3bl54mLi4NcLsfx48etXub64vDhw/D390ebNm0wceJEFBQU8J/RMbGs4uJiAICvry8AYfeqjIwMdOjQAQEBAfw88fHxKCkpwYULF6xYesekeUw4mzZtgp+fHyIiIjBr1iyUlZXxn9Exsa4Gti6Avbp37x5qa2vVTlQACAgIwKVLl2xUqvolJiYGGzZsQJs2bZCTk4N58+ahZ8+eOH/+PHJzc+Hi4gIfHx+17wQEBCA3N9c2Ba5nuP2s7RrhPsvNzYW/v7/a5w0aNICvry8dJwtJSEjA888/j7CwMFy7dg3vvPMOBgwYgIyMDDg5OdExsSClUokpU6age/fuiIiIAABB96rc3Fyt1xH3GTGetmMCAC+//DJCQkIQHByMc+fOYebMmbh8+TK++eYbAHRMrI2CVWK3BgwYwP9/x44dERMTg5CQEGzbtg1ubm42LBkh0jVixAj+/zt06ICOHTuiZcuWOHz4MPr27WvDkjm+SZMm4fz582pt64lt6Tomqm20O3TogKCgIPTt2xfXrl1Dy5YtrV3Meo+aARjJz88PTk5OdXps5uXlITAw0Ealqt98fHzw+OOP4+rVqwgMDERVVRWKiorU5qHjYz3cftZ3jQQGBtbpkFhTU4PCwkI6TlbSokUL+Pn54erVqwDomFhKcnIyvv/+exw6dAjNmjXjpwu5VwUGBmq9jrjPiHF0HRNtYmJiAEDtOqFjYj0UrBrJxcUFUVFROHDgAD9NqVTiwIEDiI2NtWHJ6q8HDx7g2rVrCAoKQlRUFJydndWOz+XLl3Hz5k06PlYSFhaGwMBAtWNQUlKC48eP88cgNjYWRUVFOH36ND/PwYMHoVQq+R8HYll//vknCgoKEBQUBICOibkxxpCcnIxvv/0WBw8eRFhYmNrnQu5VsbGx+O2339QeItLS0uDl5YXw8HDrbIgDMXRMtDl79iwAqF0ndEysyNY9vOzZli1bmEKhYBs2bGBZWVnstddeYz4+Pmq9A4nlTJ8+nR0+fJhlZ2ezY8eOsbi4OObn58fy8/MZY4xNmDCBPfbYY+zgwYPs1KlTLDY2lsXGxtq41I6ltLSUnTlzhp05c4YBYB988AE7c+YMu3HjBmOMsUWLFjEfHx+2e/dudu7cOTZkyBAWFhbGysvL+WUkJCSwzp07s+PHj7OjR4+y1q1bs5EjR9pqk+yevmNSWlrKZsyYwTIyMlh2djZLT09nXbp0Ya1bt2YVFRX8MuiYmM/EiROZt7c3O3z4MMvJyeH/lZWV8fMYulfV1NSwiIgI1r9/f3b27FmWmprKmjRpwmbNmmWLTbJ7ho7J1atX2fz589mpU6dYdnY22717N2vRogXr1asXvww6JtZFwaqJVq5cyR577DHm4uLCoqOj2S+//GLrItUbw4cPZ0FBQczFxYU1bdqUDR8+nF29epX/vLy8nCUlJbFGjRoxd3d39txzz7GcnBwbltjxHDp0iAGo82/s2LGMsUfpq2bPns0CAgKYQqFgffv2ZZcvX1ZbRkFBARs5ciTz8PBgXl5eLDExkZWWltpgaxyDvmNSVlbG+vfvz5o0acKcnZ1ZSEgIGz9+fJ0HbDom5qPtWABg69ev5+cRcq+6fv06GzBgAHNzc2N+fn5s+vTprLq62spb4xgMHZObN2+yXr16MV9fX6ZQKFirVq3YW2+9xYqLi9WWQ8fEemSMMWa9elxCCCGEEEKEozarhBBCCCFEsihYJYQQQgghkkXBKiGEEEIIkSwKVgkhhBBCiGRRsEoIIYQQQiSLglVCCCGEECJZFKwSQgghhBDJomCVEEIIIYRIFgWrhBBCCCFEsihYJYQQQgghkkXBKiGEEEIIkSwKVgkhhBBCiGT9PyYpvIY5GlMOAAAAAElFTkSuQmCC",
|
|
"text/plain": [
|
|
"<Figure size 800x200 with 1 Axes>"
|
|
]
|
|
},
|
|
"metadata": {},
|
|
"output_type": "display_data"
|
|
}
|
|
],
|
|
"source": [
|
|
"# Simple linear regression\n",
|
|
"from scipy import stats\n",
|
|
"slope, intercept, r_value, p_value, std_err = stats.linregress(np.arange(tsa_b.shape[0]), tsa_b)\n",
|
|
"\n",
|
|
"plt.figure(figsize=(8,2))\n",
|
|
"plt.stem(tsa_b - np.arange(tsa_b.shape[0]) * slope - intercept)\n",
|
|
"plt.title('linear regression on b (time stability vs. global beat)'); None"
|
|
]
|
|
}
|
|
],
|
|
"metadata": {
|
|
"kernelspec": {
|
|
"display_name": "Python 3 (ipykernel)",
|
|
"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.11.4"
|
|
}
|
|
},
|
|
"nbformat": 4,
|
|
"nbformat_minor": 5
|
|
}
|