{ "cells": [ { "cell_type": "markdown", "id": "cb5649fb-fbec-4281-aa8b-5f5c802d8e13", "metadata": {}, "source": [ "# Lifetime regressions" ] }, { "cell_type": "markdown", "id": "d27b45aa-cfc2-462f-ba23-5b6d2e0e7a4e", "metadata": { "ExecuteTime": { "end_time": "2025-09-02T12:31:51.641448Z", "start_time": "2025-09-02T12:31:51.638613Z" } }, "source": [ "A built-in dataset of insulator string lifetimes and covariates." ] }, { "cell_type": "code", "execution_count": 1, "id": "613e4f41-4c7b-4998-92ba-c32b1266356c", "metadata": { "ExecuteTime": { "end_time": "2025-09-02T12:31:52.179542Z", "start_time": "2025-09-02T12:31:52.173357Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "('time', 'event', 'entry', 'pHCl', 'pH2SO4', 'HNO3')\n" ] } ], "source": [ "from assetlife.datasets import load_insulator_string\n", "\n", "dataset = load_insulator_string()\n", "print(dataset.dtype.names)" ] }, { "cell_type": "code", "execution_count": 2, "id": "0e83fb9f-d750-4109-a1c6-d7a2639d47d3", "metadata": { "ExecuteTime": { "end_time": "2025-09-02T12:31:52.567088Z", "start_time": "2025-09-02T12:31:52.563513Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[0.49 0.76 0.43 ... 1.12 1.19 0.35]\n" ] } ], "source": [ "print(dataset[\"pHCl\"])" ] }, { "cell_type": "markdown", "id": "211ec5be-cfc0-4897-9987-44b6b8fbf50b", "metadata": {}, "source": [ "Covariates must be grouped in one ``np.ndarray``." ] }, { "cell_type": "code", "execution_count": 3, "id": "7eae8b87-3cf6-4850-90ab-f445bfbfc7a1", "metadata": {}, "outputs": [], "source": [ "from assetlife.lifetime_models import Gompertz, ParametricProportionalHazard" ] }, { "cell_type": "code", "execution_count": 4, "id": "2826f516-aa7a-4bad-9eb7-86e296de1f21", "metadata": { "ExecuteTime": { "end_time": "2025-09-02T12:31:54.407376Z", "start_time": "2025-09-02T12:31:53.549002Z" } }, "outputs": [ { "data": { "text/plain": [ "ParametricProportionalHazard(Gompertz(shape=0.21810962041584048, rate=0.02945802131524193), [4.108879842427079, -2.660506112183457, 3.234117694929174])" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "regression = ParametricProportionalHazard(Gompertz())\n", "z1, z2, z3 = dataset[\"pHCl\"], dataset[\"pH2SO4\"], dataset[\"HNO3\"]\n", "regression.fit(\n", " dataset[\"time\"], covar=(z1, z2, z3), event=dataset[\"event\"], entry=dataset[\"entry\"]\n", ")" ] }, { "cell_type": "code", "execution_count": 5, "id": "704e8ca0-4165-40a2-a473-560233bb1c57", "metadata": { "ExecuteTime": { "end_time": "2025-09-02T12:31:55.538299Z", "start_time": "2025-09-02T12:31:55.535144Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[ 4.10887984 -2.66050611 3.23411769 0.21810962 0.02945802]\n" ] } ], "source": [ "print(regression.get_params())" ] }, { "cell_type": "markdown", "id": "da36f246-9f28-4252-88a9-0a641ebd9f9b", "metadata": {}, "source": [ "Notice that the regression parameters include parameters from the baseline (Gompertz) and the covariate effect." ] }, { "cell_type": "code", "execution_count": 6, "id": "e4be562b-f824-478e-9ec9-3ed867801e3c", "metadata": { "ExecuteTime": { "end_time": "2025-09-02T12:31:56.225005Z", "start_time": "2025-09-02T12:31:56.065229Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjEAAAGdCAYAAADjWSL8AAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjExLjEsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvctoD+AAAAAlwSFlzAAAPYQAAD2EBqD+naQAAjnBJREFUeJzs3Xd4W+X58PGvtuQlz9gZzt47gYQQRhghBCiQhoaWUuYbyh7lB7RsOqEttGkLpYvSsjdllwSyIQ0jAbJ3YifxHpK1xznvH8dWosRL1pFX7s916ZJ0xqPnSLZ1+1m3QVVVFSGEEEKIHsbY1RUQQgghhOgICWKEEEII0SNJECOEEEKIHkmCGCGEEEL0SBLECCGEEKJHkiBGCCGEED2SBDFCCCGE6JEkiBFCCCFEj2Tu6gq0h6IoHDx4kMzMTAwGQ1dXRwghhBDtoKoqDQ0N9OvXD6NR/3aTHhHEHDx4kOLi4q6uhhBCCCE6oLS0lAEDBuhebo8IYjIzMwHtTcjKyuri2gghhBCiPdxuN8XFxbHvcb31iCCmqQspKytLghghhBCih0nVUBAZ2CuEEEKIHkmCGCGEEEL0SBLECCGEEKJH6hFjYtpDURRCoVBXV0OITmGxWDCZTF1dDSGE6FK9IogJhULs2bMHRVG6uipCdJrs7GyKiopk7SQhxDGrxwcxqqpSVlaGyWSiuLg4JYvpCNGdqKqKz+ejsrISgL59+3ZxjYQQomv0+CAmEong8/no168faWlpXV0dITqFw+EAoLKykj59+kjXkhDimNTjmy2i0SgAVqu1i2siROdqCtrD4XAX10QIIbpGjw9imsi4AHGskZ95IcSxrtcEMT3NY489xm9/+9s2j/voo4+47LLLOOecc1izZk2HX68j5bS3jt1db7kOIYQQ8Xr8mJieasOGDUQikVaP2bt3L+eddx6PPPIIY8aMYdiwYR16rfaU86c//Qm32829996bUB07qqKigqeeeopVq1Zx1VVXcfHFF8ftVxSF1157jSVLluB2u5k0aRI33HAD2dnZrZbb2deRKs1dhxBCiHgSxHRj69atIzMzkx/96EcpL2fTpk1UV1cn9Trt9d///pdrrrmGSy+9lM8//5xZs2YddcwFF1xARkYGZ555JpmZmfz5z3/mb3/7G1988QX5+fktlt2Z15FKveU6hBAilRIOYnbv3s37779PdXU1I0eOZP78+djt9lbPaWho4MUXX2Tfvn2MGDGC733ve22ec6x45513ePvtt/F4PJx//vl8//vfB+Cvf/0rf/7zn2loaGDu3LlYrVbefvvtZstYu3YtTz31FBUVFYwYMYJbbrmFgQMHtrucZ555hvfff59QKMTcuXMBeOSRR9qsY5OdO3fy5JNPsmvXLoqLi7nqqquYOnVqi9c8bdo0du/ejcVi4bnnnmv2mH/9619xwcr5559Pv379eP7557n11lubPSfZ6zjciy++yL///W8AnE4nU6dO5eabb46bAVdXV8cTTzzB119/TV5eHpdcckksIGttX1vvWUvXMXny5BbrK4QQx6KEgphFixbx9NNPc8YZZ5CVlcVvfvMb7rvvPj755JMW16qorq7mxBNPJDc3l9mzZ/PYY4/xpz/9iZUrV5Kenq7LRRxOVVX84aju5baHw2JKaLDlO++8Q3V1NVdeeSUHDx5k4cKFqKrKpZdeyqxZs9i1axelpaXcdtttLU6hXbp0KXPnzuWmm27izDPP5MUXX2Tq1Kls2LCBvn37tqucE088kQkTJuB2u7ntttsAKC4ubrOOAJ9//jlz585l4cKFXH755WzZsoXTTz+dV199lTlz5jRb57y8vDbfmyNbW9LT00lPT8fj8bR4TjLXcaQTTjiBnJwcAGpra/nrX//Km2++yZo1a2Kf8fnnn4/dbueHP/whHo+HX/ziFyiKwumnn97qvrbes9auQwghxCEJBTFz5szh1ltvjf0R/8lPfsKAAQN47rnnuPPOO5s95xe/+AUGg4Hly5fjcDj4v//7P0aOHMnjjz/Oj3/84+Sv4Aj+cJSxD3yoe7ntsflnZ5Nmbf9bmpubyzvvvIPFYgFg69atvPLKK1x66aWMHj2aiRMnYrVaY/+NN+euu+7i6quv5ne/+x0AF198MZMmTeKXv/wljz/+eLvKGTFiBMXFxVRXVx91TGt1BLjlllu49dZbeeCBB2LnGAwGHnzwwRaDmI544YUXKC8v55xzzmnxmGSu40hDhw5l6NChsefz5s2jsLCQ1atXc8oppxAIBPjkk09Yt24dU6ZMAeDqq6/G7Xa3ug/afs9auw4hhOhJooqa0vITCmLGjh0b9zwSiRCNRlsdbPmf//yHSy+9NLY4V25uLueffz7/+c9/UhLE9CTTpk2LfakCDBkyhK+++qrd54dCIdavX89Pf/rT2DaDwcAFF1zAhx/qE8i1VseGhgbWrl2Lqqp89tlnqKqKqqpUVFSwY8cOXV4f4Msvv+Taa6/lvvvua7WbqjWJvtfhcJgXXniB5cuXU1lZSTQaRVVVdu7cySmnnILdbmfSpEncfvvt3H777ZxyyilkZ2eTlZUF0OK+znrPhBCiO2gIpHYdq4THxJSXl/Poo4/i9XpjM0uuvPLKZo8NBoPs27fvqNkww4YNa3F8R9N5wWAw9rzpP9j2cFhMbP7Z2e0+Xk8OS2Krph65QJ/BYEgo/5PL5UJRFJxOZ9x2p9NJbW1tQnXpSB3dbjeqqjJ//nwmTpwYd5xe6R+++uor5syZw1VXXcXPfvaz2PYdO3Zw8803x55fccUVXHLJJS2Wk+h7fe2117J69WpuvvlmBg4ciM1m45tvvsHr9caOWbp0KY8//jiPPPIIF198MWeeeSZ/+9vf6NevX4v7moKWVL5nQgjRXTQEUjszNOEgxmKxUFRURH19PUajkXXr1lFXV0efPn2OOtbn8wHE/jtt4nQ6474MjvTwww/HtS4kwmAwJNSl05MVFBSQkZHBjh07OPnkk2Pbt2/fzpAhQxIqqyMLp/Xp04e0tDQsFktKuj2++eYbZs+ezSWXXMIf//jHuH2FhYWx8SIAI0eOBPRZAE5RFF566SXefPNNzj5bC4iDwSB1dXVxx+Xm5vLAAw/wwAMPUFNTw9lnn81DDz3E3/72txb3PfHEE+16z2QhOyFETxeJKvjDqQ1iEv7XLy8vjzvuuINf/OIXfP7551RWVvLQQw81e2x6ejoGg4H6+vq47fX19WRmZrb4GnfffTculyt2Ky0tTbSax4wf/OAHLFq0KPYeb9myhZdeeonLL788oXIKCgooLy9P6ByLxcJVV13Fb37zG7Zt2xbbfvDgwRZnHbXXxo0bOfPMM/ne977H448/ftT+rKws5s6dG7s1jV/pyHUcyWg0Yrfb2bVrV2zbgw8+SCAQiD2vqqriX//6V6w1Jy8vj9zcXKLRaKv72vue6XEdQgjRlRoCEdTUDolJbsVem83G1KlT2bx5c7P7rVYrw4YNi/tjDdqgyjFjxrRablZWVtxNNO+Xv/wlTqeT4cOHM3PmTI4//ni+973vcdlllyVUzoIFC9iwYQPTp09n7ty57R6b8+ijjzJnzhwmT57M9OnTmThxIqecckqryTj3798fCz6qq6t5+umnmTt3Lvfdd1/smO9///u4XC527twZF6wsWrQoJddxpF//+tfcfvvtzJw5kxEjRrB06dK4GXjp6emsWLGCfv36MWvWLIYNG0ZpaSn33ntvq/va+57pdR1CCNFVUt2VBGBQ1fbHSStWrIhb66K2tpZJkybxne98h9///vcALF68mFWrVvHzn/8cgHvuuYfnn3+eb775BqfTyYEDBxg7diy/+MUv4sY0tMbtduN0OnG5XEcFNIFAgD179jBkyJAetfbMhg0bUFU1blzEnj17KC8v58QTTwSgrKyMrVu3cvrpp7dZ3ubNmykvL4/NbDlce8txuVxs2rQJt9vNtGnTOHjwYJt1bFJRUcGWLVvIz89n9OjRmM0td+k1jac6UkFBAccddxwAq1atarbLsbi4mHHjxqXsOg5XVlbGtm3byMvLY/z48XzyyScMHDgwtgYPQE1NDZs3byYnJ4dx48bFdQO1tg/afs+OvI4jp6b31J99IcSxYeMBF35vA9NHFTf7/a2HhIKY7373u+zbt49Jkybh9/v54IMPGDduHG+88Qa5ubkAPPTQQ3HdG263mzPOOAO3280pp5zCkiVLGDVqFO+99167M0/3xiBGiGTJz74QoruKKiqf763FFPGnNIhJaATsyy+/zIYNG1i7di0Wi4VbbrmF448/Pu6YOXPmxC1UlpWVxZo1a3j//fcpKSlhwYIFzJkzR2ZiCCGEEL2UpxPGw0AHZidNmDCBCRMmtLh/5syZzJw5M26bxWLhwgsvTLx2QgghhOhx3CleH6aJNIcIIYQQQlcSxAghhBCix1EUFU8nzEwCCWKEEEIIoSNPKEKKUybFSBAjhBBCCN24/Z3TlQQSxAghhBBCR52xyF0TCWKEEEIIoQtVVSWIEUIIIUTP4w1FiXbWgBg6sE6M0MeaNWvYtm0bV155Zdz2QCDAypUr2bNnD06nk/HjxzN+/Pg2y1NVlaVLl7J9+3YGDhzInDlzsFgsrZ7j8XhYunQpBw8eZNSoUZx22mlxS+Nv3ryZxYsXH3XejTfe2GbZQgghjj2dOR4GenEQs2ZXTae+3onD8to+6DDvvfceL730UlwQ88Ybb3D99ddTWFjIjBkz8Hq9PPLII+Tl5fHmm2+2uGRzJBJh3rx5rF+/ntmzZ/PYY4+Rm5vLRx991OI569ev57zzzmPgwIFMmDCBRx99lBEjRvDWW2/F0kF89tln3H///fy///f/4s5NIFOFEEKIY0hnrQ/TpNcGMT3NihUruPjii/nNb37D7bffftS+aDTa4rn/+Mc/WLlyJZs2baK4uJi6ujomTpzIww8/zMMPP9zsOddffz0nnHACb775JgA+n4/x48fz+OOPx72+0+lsM3O0EEIIAZ07qBdkTEy38cADD3D88ccfFcAAzJo1i5ycnBbPfeWVVzj33HNj2atzcnK45JJLeOWVV1o85+uvv+ass86KPU9LS+Pkk0/mxRdfjDsuGAzyr3/9i3//+998/fXXiV6WEEKIY4Q3GCES7dyWegliuoFgMMinn37K3LlzO3T+5s2bGT16dNy20aNHs3v3bgKBQLPnDBkyhC+//DL2XFEUvvrqKzZu3Bh3nMFgYMmSJbz77ruceOKJXHrppSiK0qF6CiGE6L06uxUGpDupW6ivrycSiTBgwIAOnd/Q0EB2dnbctqaWG4/Hg91uP+qcRx55hIsvvhifz8eECRP48MMPiUajBAIBIpEIZrOZGTNmsHv3bjIyMgDYtGkT06ZN44knnuDmm2/uUF2FEEL0Tp09HgakJaZbyMrKwmAwUF5e3upxlZWVLFq0KHb77LPPAK0rqKGhIe5Yt9sd29ecCy64gI0bNzJt2jTcbjd33HEHV1xxBWlpaZjNWmw7evToWAADMG7cOObMmcOSJUs6fK1CCCF6p4YuCGKkJaYbcDgcTJkyhZUrV7Z6XDAYZO/evbHn48aNA2DkyJHs2rUr7thdu3bRr1+/FoMYgOHDh8eNwfnHP/7BtGnTWq2D1WqloqKi1WOEEEIcW/yhKKFI589clZaYbuKee+5hyZIlvPbaa0ft27VrFx6Ph+Li4riWmKaBuRdeeCHvvfce9fX1gBbsvPLKK8ybNy9WxqZNm1i0aBHhsBYp79+/P/YYYN26dbz//vtx3URbtmyJq0dZWRlLlizh5JNP1uuyhRBC9AJd0ZUE0hLTbVx00UU89thjXHbZZTz77LOxdWI2btzI9u3bWbVqVVzXzuFuuukmXn75ZU499VS+/e1v89FHHxEIBLj//vtjx6xZs4Yf/ehHLFy4EIvFwq5du/j2t7/Neeedh9vt5qmnnuLWW2/loosuip3z8MMPU1lZyYwZM/D5fDzzzDOMHTuWu+66K+XvhxBCiJ6jK7qSQFpiuszMmTO56qqr4rbdfvvt7Nq1i3PPPReXy0V6ejoLFy5kw4YN5OW1vJheWloan3zyCbfeeit+v59LLrmEr776iqKiotgx48eP59Zbb40tZDdr1iyeffZZzGYzOTk5fPTRRzz66KNx5T7zzDPcc889mEwmsrOzeeqpp1i9ejVOp1PHd0IIIURP5/J3/swkAIPaA5ZfdbvdOJ1OXC7XUSvQBgIB9uzZw5AhQ5qdhSNEbyU/+0KI7iAQjrK+pL7Zfeaon+mjipv9/taDtMQIIYQQosO6ajwMSBAjhBBCiCR0xSJ3TSSIEUIIIUSHdXbm6sNJECOEEEKIDglGogTCXZeKRoIYIYQQQnRIV3YlgQQxQgghhOigruxKAglihBBCCNFB0hIjhBBCiB4nHFXwhaJdWgcJYrpIWVkZZWVl7Tq2pqaGbdu24fF4knrNRMtJpI7dWW+5DiGE6E66uhUGJIjpMnfffTd33nlnq8coisK3v/1thg0bxoUXXsjatWs79FrtKaeiooKDBw8mXMdkRCIRtm3bRm1tbYvHNDQ0JBSAdMV1pEJz1yGEEN1JV4+HAQliurXVq1fz/vvvs3fvXrZu3cqZZ56ZsnIefPBBbrnllmSr3C5VVVU8+OCDDB06lHHjxvG3v/3tqGMWL17MzJkzGTBgAFOmTKGwsJB//vOfbZbdmdeRSr3lOoQQvVdXrtTbRIKYbqChoYGKioq4bZWVlXz22WdkZWVRXl7O9u3bWy3D7XazZ88egsFgwuVUV1dTX19PQ0MDW7duZevWrfh8vjbreLhgMMjevXuPev3mrF+/HqPRyJo1a+KSVB5u5cqV/O53v6O+vp7y8nIeffRRrrnmGlauXNliuXpcR5OamppYGa21BCmKwoEDB/B6vQntg5bfs/ZchxBCdKVINxgPAxLEdKmysjJmz57NqFGjGDZsGNOnT6eyshKAp556it/97nfU1dUxb948FixY0GwZfr+fyy+/nIKCAmbOnEleXh6/+tWvYvvbU87zzz/Phx9+yJo1a5g3bx7z5s3j66+/brOOoH1R33333RQUFHDqqaeSm5vL1Vdf3eqX7pw5c3jwwQfp379/i8f84he/YMaMGRgMBgAuu+wyBgwYwNKlS1s8J5nrONKrr74aK2PSpEn07duXt956K+6Y1157jaKiIqZNm0ZxcTHz58+nurq6zX1tvWetXYcQQnQHDYEI3SJ9tNoDuFwuFVBdLtdR+/x+v7p582bV7/drGxRFVYOerrkpSruv6YorrlCNRqP61ltvqaqqqm63W50wYYL6ox/9KHbMs88+qxYWFrZazj333KMOGjRI3bNnj6qqqrp06VLVYrGo77zzTkLlXHvttepFF12UcB1//vOfq5MmTVJLSkpUVVXV6upq9YQTTlDvuOOONt4BTf/+/dWHH364zeOqq6tVu92u/v3vf0/JdbTlueeeUzMzM9WKigpVVVU1Go2qWVlZ6j/+8Q9VVVVVURT17bffVlevXt3qPlVt33vW3HUc6aiffSGE6CR7qz3qpzur27x9tq20xe9vPZi7OojSXdgHv+rXNa99z0Gwprf78FNPPZULLrgAgMzMTC644AI++eSThF7yySef5Oc//zmDBw8G4PTTT2fBggU8+eSTfOtb30qorETrqKoqf/jDH7j33nsJBoPs2LEDVVW56KKL+Mtf/sJvf/vbpF+/6XWuu+46+vTpw3e/+13dr6M1LpeLyspKjjvuONLS0lizZg0XXnghkUgEn89HYWEhAAaDgfPPPx+AUCjU4r7Oes+EECKV6n1dPx4GoPcFMT3Ikd0pGRkZNDQ0tPv82tpa6urqmDBhQtz2iRMn8tRTT6W8jhUVFVRXV7No0SL+8pe/xB2Xlpamy+sD/OhHP+Kjjz5i2bJlZGZmAhAIBNi7d2/smIKCAvLy8losI9H3+ssvv2ThwoVs27aNwsJCbDYbdXV1HDhwAACr1covf/lLvvOd73DSSSfFgsdRo0a1uq+z3jMhhEiVYCTaLcbDQG8MYixpWotIV712J3I4HBgMBvx+f9x2v9/fKV+IFosFgD/+8Y+xVg693Xnnnfz73/9myZIlTJ48ObZ9165dXHTRRbHnt912G9ddd51ur/uDH/yAs88+m7Vr12K1WgEtEFKUQ4nO7rrrLhYuXMjHH3/MBx98wKRJk3j22WdZsGBBi/vOOOMMILXvmRBCpJKrG0ytbtL7ghiDIaEunZ7M4XAwcuRIli1bxtlnnx3b/vHHHzNlypSEyrLb7YTDif1g5uXlMXr0aF577bWjvpDD4XAsyOmou+66i7///e8sWbKE448/Pm7fuHHj2Lp161HndOQ6jhSNRtm+fTtPPvlkLIDZsmVL3Lot0WgUo9FIbm4uCxYsYMGCBQSDQV599VXmz5/f4r4FCxa06z3T4zqEECIVXN2kKwl6YxBzjPn5z3/O5ZdfzoABA5g+fTovvvgiX3zxRbvWVDlc0xfrp59+Sm5uLgMHDmzXeY899hjz5s0jJyeH7373uwQCAT7++GPKysparEMoFGL37t2AtuBdVVUVW7duJSMjgwEDBgBw//33s2jRIp566ikyMzNjAUt2dnaL07KTuY7DmUwmpk+fzq9+9SvS0tKoqqrizjvvxGQyxY4pKytj3rx53HLLLUyYMIH9+/ezfPly7rjjjlb3tfc9a+46pLtJCNHVVFWVlhgBffv2JRqN71PMy8uLDdAFcDqdjBgxotVyFixYgNFo5K9//St/+9vfGDlyJKtXr447rz3lXHHFFWzZsoXbbrsNt9vN008/3a46nnvuuaxevZpFixZxww03kJ+fz1lnncW9997b4msdOHCAefPmAVpQ8t577/Hee+9xyimn8Pe//x2Azz77jKFDh/LLX/7yqOv9+c9/rvt1HOnll1/m3nvv5Yc//CF5eXncc889vPDCC+Tk5AAwYMAA/vGPf7Bo0SL+8Ic/kJOTw3333cd1112HwWBocV9737PmruPEE09ssb5CCNEZvKEo4Wh3mFutMahqt5jp3Sq3243T6cTlcpGVlRW3LxAIsGfPHoYMGYLdbu+iGgrR+eRnXwjR2fbX+Sit9bd9YCNz1M/0UcXNfn/rQRa7E0IIIUS7dJep1U0kiBFCCCFEmyJRBU+w6zNXH06CGCGEEEK0yd1dUg0cRoIYIYQQQrSp3hfq6iocRYIYIYQQQrSpO02tbiJBjBBCCCFaFQhHCYSVtg/sZAmvE+N2u1m3bh3hcJiJEyfGkty1ZPPmzXzzzTdx26xWK/Pnz0/0pYUQQgjRBbrbrKQmCQUxDzzwAE899RQjRozAZDKxZs0a7r//fu6+++4Wz3njjTdYtGgRs2fPjm1LT0+XIEYIIYToIbpjVxIkGMQUFxezfft20tO13ETvvPMOF1xwAWeccQYnnHBCi+cNHz6cl156KbmaCiGEEKLTdbdUA4dLaEzMNddcEwtgAM4//3ysVitff/11q+f5fD7ee+89Pv74Y6qqqjpWUyGEEEJ0OncgQlTpZnOrGyWVO2nlypWEQiHGjx/f6nH79u3jj3/8I/X19XzzzTf89Kc/5a677mrx+GAwSDAYjD13u92JV27PqsTPScaQUxI6/LHHHuOdd95h+fLlsW3RaJR///vfvPbaa+zZswen08n48eO5+eabmTRpUqvlffnllzz00ENs376dgQMH8uMf/ziuC+9Iqqry9NNP89xzz3Hw4EFGjRrFAw88wHHHHRc75rXXXoslLTzc5s2bJRmhEEIcI9zdtBUGkghiampquPrqq7nggguYOXNmi8edccYZ3HzzzTidTgBeffVVLr74YqZMmcJZZ53V7DkPP/wwP/3pTztatR6hrq6O/fv3x56HQiHOO+88tm/fzoMPPsiMGTPwer1s2rSJq6++mrfeeiuW4flIu3bt4vTTT+fKK6/kZz/7Ge+++y7nnHMOK1asaPGz+cUvfsHvfvc7/vznPzNx4kTee+89TjvtND777DPGjBkDgMfjIRgMsmbNmrhzHQ6HTu+CEEKI7q67DuqFDiaAdLlczJ49G5PJxJIlS8jMzEzo/EmTJjFr1iz++Mc/Nru/uZaY4uLmE0i1mASvm7fE3Hfffbz00kvs3LkTgAcffJBf//rXbNmyhSFDhsQdq6oqiqJgMpmaLevGG29k+fLlbNq0Kbbt7LPPxmKx8O677zZ7Tr9+/fjhD3/IQw89FNs2d+5cCgsL+fe//w3Av/71L+677764YEt0H5IAUgiRauGowpf76jq8Um+3SwDpdruZM2cORqORDz/8MOEABiArK6vVsTE2m42srKy4W2/31FNP8YMf/OCoAAbAYDC0GMAALF++nDlz5sRtmzt3LitWrKClGNXr9ZKdnR23LTs7m6VLl8Ztq6mpYfLkyUyZMoWrrrqK3bt3t/OKhBBC9HQuf7jbpRo4XEJBTFMAA7B48eJYF9HhNm7cyBtvvBF7XlZWFrd/3759rFu3jmnTpnWkvr1SfX09Bw4cYNy4cR06f//+/RQVFcVtKyoqwuPxtDie6Oyzz+Zvf/sbBw8eBOCLL77g3Xff5cCBA7HAx2KxcNttt/H3v/+dxx9/HJfLxeTJk9mzZ0+H6imEEKJn6c5dSZDgmJhvfetbbNq0iUcffZQPPvggtn38+PGxwb2vvfYaixYtiq0DM3/+fKZMmcKUKVOorq7miSeeYOLEiVx77bU6XkbPFoloWUE72iWgKApmc/xHabFYAG2wcHOeeOIJfvjDHzJkyBCcTidOp5Mrr7ySP//5z6iqisFg4JJLLsFoPBTnnnDCCYwbN47f/OY3PPnkkx2qqxBCiJ6ju06tbpJQEDN48GD69evHsmXL4gsxm2NBzPjx47noooti+1asWMHzzz/P2rVrSUtL49FHH+Xiiy+O+3I81uXm5pKVlcWuXbtaPW7Dhg2cf/75sed33HEHN910E/n5+dTU1MQdW11djcViaba1DKCgoIA333yTcDhMTU0NhYWF/PjHP6Z///6xz+bIz8hsNjN9+nQ2b97ckcsUQgjRg/hCEUKR7pdq4HAJBTHPPPNMm8d85zvf4Tvf+U7sudVq5aqrruKqq65KvHbHCKPRyLx583j++ed56KGHyMjIaPa4UaNGxU3JzsnJAWD69Ol8+umncceuXr2aqVOntjqWBrQWm6KiIhRF4a233uLcc89t9fg9e/ZQUFDQjqsSQgjRk3X3riSQBJDdxq9+9SsMBgMXX3xx3ODZvXv3cvvtt3Pw4EGsViuDBw+O3ZpaWa6//npWr17N66+/DmgBzBtvvMH1118fK+fll19m8ODB+Hw+AD766CP++9//oqoqgUCAH/3oR1RVVXHffffFzrn//vvZunUroHVL/e53v+OTTz6RgFQIIY4BEsSIduvfvz+fffYZubm5TJo0icLCQrKyspg9ezaDBg1qtfXjtNNO44knnmDhwoU4nU7OPvts7rrrLq644orYMQ0NDezbtw9F0ZoGJ06cyJNPPkl+fj55eXl8+eWXLF++nOLi4tg5J510Et///vfJyckhMzOTJ598kpdffpkLLrggdW+EEEKILqcoKg2B7h/EdGidmM7mdrtxOp2JrRPTzdXX1+P1eunfv/9R+6LRKFVVVWRlZSW0Mm4kEqGyspK8vDxsNlvcPo/HQ3V1NYMGDcJgMMS2N81eam0au9vtxmw2yyq93UxP/dkXQnR/9b4QW8oaki4n1evEJJV2QHRcdnb2Ueu0NDGZTEdNmW4Ps9lMv379mt2XkZHR7Fib9vxQHQvr9AghhDikJ3QlgXQnCSGEEOII3X1qdRMJYoQQQggRE4xE8YWaX2Osu5EgRgghhBAxrh7SlQQSxAghhBDiMD2lKwkkiBFCCCFEI1VVJYgRQgghRM/jCUYIR7v9yisxEsQIIYQQAuhZXUkgQYwQQgghGtV5e1YQI4vddZF///vfKIrSZh6iTZs28dJLL1FdXc11113HpEmTOvR6HSmnvXXs7nrLdQghRCoFI1E8wUhXVyMh0hLTRZYtW8bHH3/c6jEVFRWccMIJuN1uJk2aFEv4mKj2lPPCCy/wt7/9LeE6dlQwGOSFF17g+uuv58MPP2z2mM8++4xf/vKX/PjHP+aFF14gEmn7l6uzryNVmrsOIYRIJb1bYVRVTXnLjgQx3dj//vc/LBYLf/jDH7juuusYPHhwyspZuXIlixcvTq7C7bRq1SqGDRvG22+/zYsvvsj69euPOubKK6/k1ltvxefzkZOTw89+9jNmzJgRy8Ldks68jlTqLdchhOg5ar0hXcsrqfVxzTNf6FrmkaQ7qYutX7+e999/H4/Hw3nnncfJJ58MwKuvvspf/vIXgsEg1113HRaLhT/96U/NlrFv3z6ef/55KioqGDFiBJdffnks31F7ynn77bdZtWoVgUCA6667DoDbbrutzTo2qaur47nnnmPXrl0UFxfz3e9+lwEDBrR4zYMGDWL9+vUUFBS0eNydd97JuHHjYs8XLlzIwIEDefbZZ7n22mubPSfZ6zjcBx98wFtvvQWA0+lk6tSpLFiwAKPxUNwfDod5+eWX+frrr8nLy2P+/PmMHDmyzX1tvWctXcfo0aNbrK8QQiQjElVw65y1+n+7a0n1PCdpielCS5cuZeHChZhMJrxeL2eccQb//e9/Aejfvz+DBg3CYrEwefLkFsewrF+/nvHjx/PFF1/Qr18/nn32WaZOnYrL5Wp3OX379iUvLw+n08nkyZOZPHkymZmZbdYRYMeOHUyYMIHVq1czcOBAtm7dysSJE/nii5aj74EDB1JQUNDqe3N4AAOQn59PQUEB5eXlLZ6TzHUcqbCwMFaG0+nkoYce4tvf/nbcMRdddBG//vWvKSgowOfz8b3vfY+1a9e2ua+t96y16xBCiFSo84VRdYw4VFVl7Z4a/QpsQa9riVFVFX/E3yWv7TA7MBgM7T5eVVVWrFgRyy5dW1vLP/7xD+bOncvMmTPZvXs377//fuy/8ebcfvvtzJ07l1dffRXQ/mMfM2YMv/71r/nVr37VrnKmTZvG2LFjY4N+21tHgBtvvJEFCxbw+9//PnZOXl4eP/7xj3Udh7J48WJKSko47bTTWjwmmes40tSpU5k6dWrs+bXXXsuAAQP4/PPPmTZtGoFAgHfffZfPP/+c4447DoB7772Xurq6VvdB2+9Za9chhBCpUOfTvyupzBXAYmr/d2JH9Logxh/xc8ILJ3TJa6/9/lrSLGntPn7mzJmxL1WAMWPG8Pbbb7f7/EgkwurVq2MBDIDNZmPBggUsW7as3eV0tI4+n4+lS5dis9m46aabUFUVVVXZvXs369at0+X1Afbs2cNll13G1VdfzaxZszpURkfe6xUrVrB8+XIqKyuJRqPYbDa2bNnCtGnTsNvtDB48mMcee4y7776b8ePHY7PZKCoqAmhxX2e9Z0II0V6KolKvc76ktXtqAZhSnMNuXUuOJ91JXcjhcMQ9N5lMRKPtzxxaW1tLJBI5qmumoKCAioqKlNexrq6OaDTKqFGjGD9+PBMmTGDixInMmzePX//617q8fklJCWeccQYnnXQSf/3rX+O2X3fddbHbBx980OHraM7dd9/NggULcLvdjBw5ksmTJ2Oz2XC73bFjlixZgt1uZ+7cuRQWFnLTTTfF9re0rzPeMyGESITLHyaq6NeXpKoqa3drXUkzhuXqVm5zel1LjMPsYO3313bZa3emgoIC7HY7JSUlnHTSSbHt+/btY+DAgQmVlUg3WJP8/HysVisjRoxocbBtMkpLSzn99NOZPHkyL7/8MmbzoR9Xh8PB5MmTY88LCwuBjl3HkVRV5YknnuCZZ55h3rx5ACiKwh133BF33LBhw/jnP/8JwDfffMOCBQuwWCz8/ve/b3HfI4880q73TI/rEEKI9qjVuSuptM7PwcaupOMH5eha9pF6XUuMwWAgzZLWJbfO/uIxGAzMmzePJ598klBI+yEsLy/nxRdfPGoQaltycnJiYzbay2azcdFFF/Hb3/6Wqqqq2Hafz5f09OD9+/dz2mmnMXHiRF555RUsFkvc/oKCgriWmKbxKx25juYoihI3nfvxxx/H6/XGntfV1bF8+fLY84kTJzJu3Diqq6tb3dfe90yv6xBCiNaoqkq9zkFMUyvMxAHZpNtS21bS61pijjW/+c1vOOOMM2KzWJYtW8bUqVO54YYbEirnvPPO47HHHuPiiy8mNzc3bmpyax5//HG+/e1vM3bsWE4//XQCgQAbN27krrvuavGciooKHnzwQQDq6+v5z3/+w969exk3bhw333wzAAsWLKC0tJTTTjsttg3g5JNP5gc/+IHu13E4g8HAnXfeyQ9/+EPeffddqqqq2L59O/n5+bFjTCYTP/vZz7j99tuZMGEC+/fv56uvvmLJkiWt7mvve9bcdcgUayGE3hqCEUIRfbuS/tc4K+mEIantSgIJYrrMFVdcgXrEfLa5c+fGrSVywgkn8Jvf/KbVcoqLi9m4cSMff/wx5eXl3HDDDUetf9Keck466SQ2btzI559/jtvtJjMzs111zM3NZcWKFXz22Wds3ryZ/Px8TjzxRPLy8lp8LZvNFusKOrxL6PAusFtuuSU2TfzI603FdRzpwQcf5Fvf+hYbNmwgLy+PM888k3feeYfx48cDkJWVxdKlS/nqq6/YsGEDOTk5nH766aSnpwO0uq8971lz1yGEEHqr03mBu/11fg7WBzAbDRw3KAdI7Yq9BvXIv+7dkNvtxul04nK5You4NQkEAuzZs4chQ4Zgt9u7qIZCdD752RdCJGt9SR2BsKJbea9+Wcob6w4wdWAOd549CnPUz/RRxc1+f+uh142JEUIIIUTbfKGIrgEMwNrd2tTqGUNT35UEEsQIIYQQxyS9cyWV1vo4UO8/rCsp9SSIEUIIIY5BemeYbkozMHGAkzRr5wy5lSBGCCGEOMYEI1E8wYiuZf6vsSvphCEtT+zQmwQxQgghxDFG71aY/XVaV5KpE7uSQIIYIYQQ4pij93iYplxJE/s7U77A3eEkiBFCCCGOIZGogjug83iYxlV6TxjaeV1JIEGMEEIIcUyp84XRc4W4A/V+Suu0rqRU50o6kgQxQgghxDGkLkW5kiZ0clcSSNqBLrNlyxYOHDjA7Nmzj9q3fft29uzZg9PpZNSoUeTktC+y3bVrF9u3b2fgwIGMGzeuzeMVRWHjxo0cPHiQkSNHMnTo0Lj9JSUlrFu37qjzzj//fEwmU7vqJIQQovtQFJV6n75dSf/b0zQrqXMWuDtcrw1iPi//vFNfb1rRtISOf/7553nppZfYuXNnbNv//vc/rr32Wg4ePMiUKVPwer3s2rWLc889l8cff5y0tLRmy1JVlZtuuol///vfHHfccXzzzTfMmjWLV155BavV2uw5u3fv5qKLLqKmpoZRo0bx5ZdfcvHFF/Pkk0/GsnEvXbqU66+/nrPPPjvu3Llz50oQI4QQPZDLHyaq6NeXdLDeT2mtD5PBwPGDJIg5Zn399decccYZXHvttfz617+OBR/RaJTnnnsOr9fbYhDz0ksv8c9//pO1a9cyceJESktLmTp1Kr///e/58Y9/3Ow51157LU6nk7Vr12K1WqmoqGDSpElMnz6dq6++OnZcXl4e//nPf3S/XiGEEJ2vRudZSf9r7Eoa3z+LDHvnhxQyJqabuOeeexg8eDCPPvpoXOuJyWTiiiuuoKCgoMVzn332WebMmcPEiRMBLdPzd7/7XZ555pkWz/nkk0/4zne+E3utwsJCZs+ezdNPPx13XCQSYfny5axYsYLq6upkLlEIIUQXUlWVer3HwzR1JXXyrKQmEsR0A5FIhI8//pj58+d3qJvm66+/ZtKkSXHbJk+ezNatWwmFmv+BLSoqiuvKAm1MzVdffRW3ra6ujnvvvZc77riD/v37c++99yZcPyGEEF3PHYgQjurXlVRW76eksStpWhd0JYF0J3ULtbW1BINBBg0a1KHzXS7XUYN/8/LyUBSFhoYG8vKOjpDvuecebrjhBhwOBxMmTOCDDz5g3759eDweIpEIZrOZ8ePHs3PnToqLiwH46KOPmDt3LmPGjOEHP/hBh+oqhBCia9Tp3ZXU2Aozrou6kkBaYrqFprEutbW1rR5XX1/Pf/7zn9ht+/btAFitVnw+X9yxHo8HAJvN1mxZCxcuZMmSJXi9Xt5++22mT5/O7bffjs1mw2zWfhiPP/74WAADMHv2bObMmcObb77ZsQsVQgjRZVI1tXpGJ+ZKOpK0xHQDGRkZjB49ms8++6zV42pqavjXv/4Ve37ppZfGpkaXlpbGHVtaWkp+fj4ZGRktljdr1ixmzZoVe/69733vqG6pI2VnZx/1WkIIIbo3bzBCIKzoVl6Zy8++Wh9GAxw/uHMXuDucBDHdxG233cb111/P6tWrOfnkk+P2uVwubDYbw4YNa3am0DnnnMO//vUvgsEgNpsNRVF4/fXXOeecc2LH7Nu3j/Xr18fWeKmvryc7Ozu2f8+ePbz99tv84Q9/iG2rrq4mPz8/9tztdrN8+XIWLFig34ULIYRIOf1bYbSeg/H9nGTaLbqWnQgJYrqJa6+9lk2bNnH22WezcOFCZsyYgdfrZePGjbz99tusWbMGu93e7Lm33XYbzz77LN/61rf4/ve/zwcffMDu3bt5+eWXY8csWbKEa665hoaGBjIyMvj00095/PHHufjii3G73Tz66KOcd955LFy4MHbO1VdfTb9+/ZgxYwY+n48nn3wSh8PBT37yk5S/H0IIIfRT49F7PEzX5Eo6koyJ6SJjx47lrLPOitv2xz/+kRUrVpCWlsZbb73F559/zvDhw1m/fj2FhYUtlpWXl8dnn33GjBkz+OCDDxg4cCDr1q1j+PDhsWMGDx7MhRdeGBvvcu6553LHHXewevVqvv76a373u9/xyiuvxBa6A3jjjTeYOXMmq1at4quvvuK6665j06ZN9O3bV+d3QwghRKp4gxF8oahu5ZW7Auyr6fquJACDquqZBio13G43TqcTl8tFVlZW3L5AIMCePXsYMmRIiy0VQvRG8rMvhGiPkhofB+r9upX3xrr9vPrlfib0d3LPuWNaPdYc9TN9VHGz39966FBLjNvtpq6uLqFzAoEApaWlLa5bIoQQQgj9VXmCupWlqiqrdmgLn540PL+No1MvoSDm9ddf57jjjqO4uJihQ4cyfPhw3n333TbPe/DBB8nNzWXy5Mnk5+fzpz/9qcMVFkIIIUT7uPxhQhH9ZiVtr/BQ7g5gMxu7JOHjkRIKYpYtW8Y//vEP6uvrqa2t5aqrruI73/lObL2S5jz77LP89re/5aOPPqKmpobnnnuOH/3oRyxZsiTpygshhBCiZTU6tsIArNxRBcCMoXnYLV2fCDihIObxxx9nypQpGAwGDAYDP/nJT1AUhZUrV7Z4zpNPPslFF13EzJkzAbjgggs45ZRTePLJJ5OruRBCCCFapCiqrgkfg5Eoa3Zps5JOHdlyPr/OlNTspN27dxMOh+nfv3+z+xVFYd26dZx44olx208++WS++OKLZF5aCCGEEK2o94eJ6Jgr6Yu9dfjDUQoybIwuytSt3GR0eJ2YSCTCNddcw6RJk46aKtykoaGBYDB4VO6e/Pz8VjMiB4NBgsFDTWBut7vN+vSASVZC6Ep+5oUQranWuStpxXatK+nUkfkYD1uOoyt1qCVGURSuuOIKduzYwRtvvBFbe+Sowo1a8eFwOG57KBRqNVvzww8/jNPpjN0Oz99zpKZyZNaTONY05cuyWLputUwhRPcUVVRdEz7WeIJsPOAC4JQR3aMrCTrQEqMoCldddRXLli1j+fLlDB06tMVjMzMzcTqdlJeXx20vLy9nwIABLZ539913c/vtt8eeu93uFgMZs9lMWloaVVVVWCyWWOAkRG+lqio+n4/Kykqys7Nb/YdACHFsqvEGUXRsrF21sxoVGNM3k8Ks7rMuVUJBjKIoXH311SxZsoRly5YxcuTIo46pra3F5XIxZMgQAE477TQ+/PBD7rjjjtgxH3zwAaeddlqLr2Oz2VrMvnwkg8FA37592bNnD/v27UvkcoTo0bKzsykqKurqagghuqHqBv1aYVRVZWVTV1I3aoWBBIOYH/7wh7z++uu8/PLLmEwmdu7cCUBubi65udp88T/+8Y8sWrSI+vp6AO69915OPvlkHnjgAc4//3z+9a9/UVpayv/93//pdhFWq5URI0ZIl5I4ZlgsFmmBEUI0KxRRcAfCbR/YTjsqPZS5mtaG6dpcSUdKKIhZs2YNhYWF3HLLLXHbb7nllti23NzcuC6madOmsXjxYh555BFef/11RowYwYoVK+Ly+ujBaDTK0utCCCGOeTXeIHqO+29qhZk+JBeHtXv989TjcycJIYQQ4pAN+114ghFdygpFFK5//kt8oSj3nTeGcf2cCZ3fLXMnCSGEEKL7CYSjugUwAF/sq8UXipKfYWVM3+7XiCBBjBBCCNFLVDWkaG2YEQXdZm2Yw0kQI4QQQvQSei5wV+sNsaFxbZjukmbgSBLECCGEEL2AJxghENYvY/WqHVWoKowu6l5rwxxOghghhBCiF6jWsStJVdVYxuru2goDEsQIIYQQPZ6qqtR49QtidlV5OFgfwGoycsKQXN3K1ZsEMUIIIUQP5/KHCUX0WzFlxWFrw6RZO5wrOuUkiBFCCCF6uGqPfivWhyIKn+6qAWBWN+5KAglihBBCiB5NUVRqdcxY/WXj2jB56VbG9ut+a8McToIYIYQQoger9YWI6piyeuWOakAb0Nsd14Y5nAQxQgghRA9Wo2NXUq03xNf76wE4ZUS+buWmigQxQgghRA8VjirU+fQLYlbvrEZVYVRhJn2dDt3KTRUJYoQQQogeqtYb0i1jtaqqsYzV3XltmMNJECOEEEL0UHrmStpV5eVAvR+ryciMod13bZjDSRAjhBBC9ED+UJSGgH4Zq5dvqwRgWjdfG+ZwEsQIIYQQPVCFO6BbWd5ghNU7tVlJZ4zqGV1JIEGMEEII0eMoikqVjhmrV2yvIhhRKM5xMKZv914b5nASxAghhBA9TI03RCSqz4heRVVZsrkCgDnjijB087VhDidBjBBCCNHD6NmV9M1+F+XuAGlWEycP7/5rwxxOghghhBCiB/GFIroO6F28qRzQ8iTZLSbdyu0MEsQIIYQQPUiFW7+xMBXuAF+V1gNw1thC3crtLBLECCGEED2EoqhU6zigd/HmClRg0gBnj1ih90gSxAghhBA9RLU3qNuA3kA4yorGtWHmjCvSpczOJkGMEEII0UNU6tiV9MmuaryhKH0ybUwekK1buZ1JghghhBCiB/AG9RvQq6oqizc1TqseW4TR2HOmVR9OghghhBCiB6jUMU/StvIGSmp9WE1GZvWgFXqPJEGMEEII0c1FdR7Q++FmbVr1ySPyybD1jDxJzZEgRgghhOjmajz6Deit9Yb4fE8dAHN64LTqw0kQI4QQQnRzeq4N8/HWCqKqyuiiTAblpetWbleQIEYIIYToxrzBCJ6gPgN6I1GFj7c0Tqse2zOnVR9OghghhBCiG9MzT9LaPbW4/GGy0yxMG5KjW7ldRYIYIYQQopvSBvSGdCtvceOA3tljCjEbe34I0POvQAghhOilqj1Booo+A3r3VHvZXuHBZDRw5ug+upTZ1SSIEUIIIbopPVfobcpWfcKQXLLTrLqV25UkiBFCCCG6IY+OA3obAmE+2VUNwNk9NE9ScySIEUIIIbohPQf0Lt9WRTiqMjgvjRF9MnQrt6tJECOEEEJ0M5GoQo1OA3oVRWXJ5sY8SeOKMBh6Zp6k5kgQI4QQQnQzNd6QbgN615fWU+UJkmEzc9KwfF3K7C4kiBFCCCG6GT27kt7fUAbAaaMKsJp719d+77oaIYQQoodz+cJ4g1FdytpR0cDmMjcmg6FXDehtIkGMEEII0Y0cdPl1K+utrw8CWrbq/AybbuW2V156al+z5+bfFkIIIXoZXyhCvS+sS1kltT6+3FeHAbhgUj9dykxEmtVEXkZqwwxpiRFCCCG6iYP1+o2FefurAwBMH5JLv2yHbuW214AcR8pnQkkQI4QQQnQDwUiUao8+K/RWuAN8ursGgAsn99elzERk2MzkdUL3lQQxQgghRDdQ7gqg6jOrmne+PoiqwqQBTobkp+tTaAKKczun5UeCGCGEEKKLRRWVygZ9WmFqvSFWbK8CYF4XtMJkOcydlptJghghhBCii1W4A0Si+jTDvL+hjIiiMqowk9F9s3QpMxHFuWmd9loJBTHhcJhXX32VM844g+zsbF555ZU2z3n00UfJzs6Ouw0aNKjDFRZCCCF6E1VVKXPpM6C3IRDmoy1aioELJ3f+jKTsNAtZdkunvV5Cc58ef/xxPvnkE+655x7OOussQqG28zoEAgFGjBjBkiVLYtt6U94GIYQQIhnVnhChiKJLWR9uKicYURiUm8bk4mxdykxEZ7bCQIJBzG233caPfvSjhF/EZDKRnZ2d8HlCCCFEb1em0+J2/lCU/24qB7QZSZ3dYJCXYSXD1rnLzyXUndTRN2Tjxo0UFxczYsQILr30Uvbs2dOhcoQQQojepN4X0i3FwMdbK/AGoxRl2TlhSK4uZbaXwQDFOZ3bCgOdMLA3JyeH3/72t6xcuZJXXnkFt9vNtGnTqKioaPGcYDCI2+2OuwkhhBC9jV6L24UiCu81Jnq8YHI/jMbObYXJz7DhsJo69TWhE4KYG2+8keuvv54hQ4YwZcoUXnnlFUwmE3/5y19aPOfhhx/G6XTGbsXFxamuphBCCNGpvMEILr8+KQZW7qii3hcmN93KKcPzdSmzvQwGbXXertDpU6wdDgfjxo1j27ZtLR5z991343K5YrfS0tJOrKEQQgiRenqNhYkqKu80Jnr81sS+mE2d+9XeJ9OG3dL5rTDQBQkgo9EoO3fuZOLEiS0eY7PZsNk6P9umEEII0RkC4SjVnrZn+LbHmt01VDYEybSbOX1UH13KbC+jAfp3USsMpKAl5pFHHolbB2bhwoV88803RKNR6urquPHGG6moqODqq6/W+6WFEEKIHkGvFAOKqsYSPZ4zvm+nt4gUOe3YzF3TCgMJBjErVqyILVgHcO2115Kdnc0tt9wSOyYQCOByuWLP58+fzzXXXENWVhYDBgxg69atLF++vNWWGCGEEKK3ikQV3VIMrCupo7TOj8NiYs7YQl3KbC+T0dAl2bEPZ1DV9seCkUgEj8dz1HabzYbDoV1IIBAgGAzidDrjjgmHw1gsHVvFz+1243Q6cblcZGV1/hLKQgghhF4O1PspqfElXY6qqtz/1kZ2VXm5YFI/Lpk+UIfatd+AHEebi9ul+vs7oTExZrO5zUXr7HY7drv9qO0dDWCEEEKI3kJRVMp1GtC76aCbXVVeLCYD54wv0qXM9rKYDPR1Hv1d39kkAaQQQgjRSaq9QUKR5AfDqKrKm+u1sTCnj+rTaVmjm/TPcXT6LKjmdH0NhBBCiGNEmU6L22044GJzmRuz0cC3JnZuokeH1URRVte3woAEMUIIIUSnqPEE8YWSTzGgqCovfFYCwFljCynI7NwlSQbnpXWbRM4SxAghhBAppqoq++v0GQuzZlcN+2p8OCwm5k3pr0uZ7ZWTbun0rqvWSBAjhBBCpFi1J6RLK0wkqvDKF9oq9t+a2Jcse+dNmjEaYHBeeqe9XntIECOEEEKkkNYKk/yUaoCPt1ZS2RDE6bBw7oS+upTZXkVOe5elF2iJBDFCCCFEClU2BAmElaTL8YeivNE4I2n+1P6dGlBYzQYG5LS+JkxXkCBGCCGESBFF0W8szHsbynD7wxRl2TljdOfmSCrOTcNk7B6DeQ8nQYwQQgiRIhUNAUKR5FthXP4w723QMlVffPwAzMbO+/rOtJvpk9nBKdVK8tfeGglihBBCiBSIKioH6/VphfnP+gMEwgpD8tM5YWieLmW216C8JLqR3Af1q0gzJIgRQgghUqDcHdBldd4Kd4AlWyoAuGT6QIyduEZLQaaNzI7OgPLVQtCtb4WOIEGMEEIIobNIVKFMp1aYV7/cT1RRmdDfyYT+zrZP0InJaGBgGwkeW6SqULdX1/o0R4IYIYQQQmdlrgDhaPKtMHtrvHy6sxqg07NU989xYDV3MEzwVELIq2+FmiFBjBBCCKGjcFSh3K1PjqSXPitBBU4cmseQ/M5baM5uMdK3o/mRFAXqS/StUAskiBFCCCF0VFYfIKJDK8ymgy6+3u/CZDBw8fHFOtSs/QblpWPs6JTqhoMQ0SeIa4sEMUIIIYROQhF9WmFUVeXFxiSPZ4zpQ5Gz87JGZ6dZyE3vYH6kaARc+/WtUCskiBFCCCF0cqDeT1RJvhXms7217KryYjMbmd+JSR4NyeZHcpVCNKxfhdogQYwQQgihg0A4SqUOrTBRReXlz7Ukj+dN6NupWaOLsuw4rB1MZxAJQkOZvhVqgwQxQgghhA4O1PvRoRGG5dsqKXMFyLSbOW9i5yV5tFmMFHd0SjVog3mV5DN1J0KCGCGEECJJgXCUqoZg0uX4QhFe+1IbU/LtKf1Js5qTLrO9huandzw/UsgHngp9K9QOEsQIIYQQSdpf50PVoRXmtS/3U9+Y5HH2mMLkC2yngkxbct1WdXvR5Q1IkAQxQgghRBJ8oQjVnlDS5eyr8fLhpnIArpw5GIupc76irWZDcvmRAm7w1ehXoQRIECOEEEIkYU+1N+lGCEVVefqTvSgqTB+Sy6TibF3q1h6D89KTC5jq9uhXmQRJECOEEEJ0ULUniNsfSbqcVTuq2FbRgM1s5PIZg3SoWfvkplvJy7B1vABfrdYS00UkiBFCCCE6IKqolNT6ki7HE4zwwlptYbv5UwckF1QkwGwyJJfKQFW7tBUGJIgRQgghOuRgvZ9gWEm6nFe+KMUdiNA/28G544t0qFn7DMpN63iCR2hM8ph8EJcMCWKEEEKIBAXCUQ7W+5MuZ3eVh482a1OTrz5pMOZOGszrdFjo09EEj9CY5HGffhXqIAlihBBCiATtq/ElvbCdoqr885M9qMBJw/IY28+pS93aYjIaGFqQZEZs9wFthd4uJkGMEEIIkQCXL0ytN/kp1cu2VrKryovDYuLSThzMOyDHgd3SwdQCAOGAliOpG5AgRgghhGgnVVXZW+NNuhx3IMxLjfmRvnPcAHI6KT9Spt1M32QzYtfu7vT0Ai2RIEYIIYRop3J3AF8o+S/wlz4rwROMMDA3jbPHdc5gXoMBhhakYzB0MLUAaFOqu2hhu+ZIECOEEEK0QziqsL8u+cG82ysaWLatCoCrTxrS8XxFCeqf7UguF5MShZpd+lVIBxLECCGEEO1QUusjEk1uNG9U0QbzAswaWcCookw9qtamNKuJ/tmO5AqpL4FIQJ8K6USCGCGEEKINnmCESnfys3GWbK5gX42PdKuJS6YP1KFmbTMYYFifDIzJtPiEfOA+qF+ldCJBjBBCCNGGvdXJD+at94V45QttMO93pxXjdFiSLrM9inPTyLAl0Y0EULMT1OQX9tObBDFCCCFEKyobAjQEks+P9PzaEvzhKEPz0zlzdKEONWub02FJvhvJUwkBlz4V0pkEMUIIIUQLoopKqQ75kdaX1LF6ZzUG4OqThyTXtdNOFpOB4X0ykiskGoHars2P1BoJYoQQQogW7K/zEYokN5jXE4jwt1W7AZg7vohhBUkGFu00rCAjudxIoKUWiCa/sF+qSBAjhBBCNMMfilLmSn42zj8/3UO9L0y/bDvfm9Y5g3mLnHZy0pNcQC/YAA1l+lQoRSSIEUIIIY6gqiq7qjyoSeZHWrOrhjW7ajAa4PpZw5NvGWmHNKuJQblpyRWiqo2DeZN8A1JMghghhBDiCOXu5Afz1vtCsTVhLpzcP/nxKe1gNMCIwiSnU4PWAhP06FOpFJIgRgghhDhMIBylpCa5wbyqqvL3VbvxBCMMzktj/pT+OtWudYPz05NblRcgEoK6ffpUKMUkiBFCCCEaqarKzkoPSpK9KCu2V7GupB6z0cANpw3HbEr9121uupXCrCSTOwLU7QEl+SnlnUGCGCGEEKKRHt1IVQ1BnlmjtWQsOL6Y4mTHp7SD1WxkaEF68gX567V1YfSgqrDhNX3KaoEEMUIIIQRaN1JpbXIJHhVV5S8rduEPRxlZmMG3JvTVqXYtMxhgeJ8MLMm29ihRbTCvXnYshs//pl95zZAgRgghhAB2VnqIJtmPtHhTBZvL3NjMRq6fNbxTFrXrn+3QJ4VB7W4IJ5+lG4D6UvgstQEMdCCI2bVrF3fddRezZ8/m448/btc527Zt44YbbuC8887jtttuo7S0NOGKCiGEEKlS5vIn3Y1UVu/nxc9KAPj+CQMpcuowPqUNmXYzA3KSTCsA4KuFhvLkywFtcbyVv4ZoEPpO0afMFiQUxDz11FOcffbZ5Obm8vHHH1NW1vYiONu3b2f69Ol4vV6uvPJK9uzZw/Tp0ykv1+nNEkIIIZKgRzdSVFH584pdhKIKE/o7OWtM6nMjNaUVMBiSbO2JhqF6hz6VAvjin1C3F+xOmPVj/cptRkJBzIUXXsiOHTv4yU9+0u5zfvaznzF69Gj+/e9/s2DBAl577TVsNhuPPfZYwpUVQggh9LarKvlupHe+OcjOSg9pVhPXnjo0+cCiDQYDjOiTid1iSr6w6h36pRYo+R9sfVd7fPLtkJarT7ktSCiIyc/PT/iDWbx4MRdeeGHsucVi4Vvf+haLFy9OqBwhhBBCb+WuAG5/ct1I+2q8vPblfgCuOHEweRk2ParWqkF5aTjTdBgH01ABvprkywHwVsOnf9Aej/029D9On3JbkeSKOK3z+XxUVVUxYMCAuO0DBgxg7969LZ4XDAYJBoOx5263O1VVFEIIcYwKhKOUJJmhOhxV+PPyXUQVleMH5XDKiHydateygkwrfZ06jIMJB7TBvHpQorDqUS3fUt5wmHq5PuW2IaWzk0IhrXnK4Yh/s9PS0mL7mvPwww/jdDpjt+Li4lRWUwghxDFIj26kZ9bspaTWR5bdzMJTUt+NlGEzMzRfh/QFqgrV2/Vb1G7DK1CxEcwOOPVOMOnQStQOKQ1iMjIyMJvN1NbWxm2vqakhJyenxfPuvvtuXC5X7CazmYQQQuhJj26kVTuq+GhLJQbghtOG6zPNuRVWs4GRRTrkRQJwH4CAK/lyACo3w9cvao9nXA9ZnZNiAVLcnWQ2mxk/fjzr1q2L275u3TomT57c4nk2mw2bLfV9ikIIIY49enQjldb6+McqLbnj/Kn9mVScrUPNWmYwwIjCTGxmHQbyBj365UYKemDlb0FVYOjpMOwMfcptJ91bYp555hnmzZsXe37llVfy6quvsnOntgrgF198weLFi7nyyiv1fmkhhBCiVU25kZLpRvKFIvz+o+2EogoT+zuZP2VA2yclaUh+Oll2HVp6FEXrRlKV5MtSVVjzR/BWQWZfrRWmkyUUxKxfv57Zs2cze/ZsAB555BFmz57N7373u9gxu3fvZvny5bHnN910E+effz6TJk3iuOOO45RTTuHGG2/k4osv1ucKhBBCiHbaV+NLalE7VVX568rdlLkC5KVbufGM1K/K2yfLpk9iR4D6vRDy6lPW9v/Cvk/BaNbWg7GkPkfUkQyqqrY7HK2trT2qawigf//+jBkzBtCCmNLSUmbNmhV3zL59+ygpKWHYsGH069cvoUq63W6cTicul4usrKyEzhVCCCEAar0htpU3JFXG+xvKePZ/+zAZDTx0/liG98nUqXbNy7SbGds3S59AyV+vDb5t/9d+y+r2wXs/0taXOf5qGDe/2cPcEQvOkSem7Ps7oTExubm5sVaYlgwdOpShQ4cetX3QoEEMGjQosdoJIYQQOgiEo+yq8iRVxrbyBl5Yq6UVuGzGoJQHMFazkZGFmfoEMEpUW9ROjwAmEmxMKxDS1oIZOy/5MjtIEkAKIYTo1RRFZUeFh0i041/gLn+YP3y8naiqcuKwPOaMTW1aAaMBRhZmYDXr9DVdswsigeTLUVVY8zjUl4A9G066DQxdF0pIECOEEKJX21frwxPs+DgYRVH509Id1PnC9M928MNOWA9mSEE6mXoM5AVwl4GnQp+yNr0Ju5dpgcupd4Kj5eVSOoMEMUIIIXqtGk+QcldyLRCvflnKpoNubGYjP5o9Up98Ra0YkOOgT6ZOA3kDbv1W5d3/OXz5tPZ4+g+h7yR9yk2CBDFCCCF6pUA4yu7q5GbirNtXx3++OgjAD08dSv8cHZb7b0WfLBvFuTrN8omEoGqrPtOp60u19WBQYeRcGHVe8mXqQIIYIYQQvY6iqGyvaEhqHEylO8Cfl2trnM0ZW8jMYanNi5STbmFofro+hamqFsBEgm0f25ZgAyz9GYR9UDgepl+rrb7XDUgQI4QQotfZW+PFG4x2+PxQRGHRxzvwhqIM75PBZTNSO7s2025mZJ9M/cba1O3RJ62AEoUVv4aGMkjvA6fd3Wl5kdpDghghhBC9SrUnSIW74y0Qiqry5Iqd7Kn2kmEzc+uZIzCbUvd16bCaGFWk01RqAE8VuA7oU9YXT0HZV2C2wxn3g92pT7k6kSBGCCFEr+EPRdldldw4mJc+K+F/u2sxGQ3cNnsE+Rmpy+VnNRsZXZSJRa8gKeSFmh36lLVjMWx5W3t88u2QO0SfcnUkQYwQQoheQVFUdlQ2JJUXacnmCt75pgyAa08dyrh+qWt5MJsMjOmbqd9sp2gEKrdqXUDJqtgE//uz9njypTBoZvJlpoAEMUIIIXqF3dXJjYNZX1LH059qmam/c9wAThlRoFfVjqItZpdJmjWhhfNbV71dG3ybLE8lLP8VKBEYdDJM/F7yZaaIBDFCCCF6vP11PqoaOj4OZk+1lz98vANVhVkjC5g/pb+OtYtnMMDwPhk4HToOkK0vAV9N8uWEA7D059qg4NxhjSvydo+ZSM2RIEYIIUSPVu0JUlrr7/D5NZ4gv/1wK8GIwvh+WSw8eUhKV+QdnJdOnp7jbHy1WhCTLFWBT36nzWyyZ8Pp94JFp0X3UkSCGCGEED1WQyDMrsqOJ3b0hSL8+r9bqfOFGZDj4LbZI1M6E2lAjoMip46BQTigdSPpkdhx/bOw71MwmuH0eyCjT/JlppgEMUIIIXqkQDjKtvIGOjqON6IoLPpoB6V1frIdFu46ezTpNh3HqByhyGnXbzVe0AbwVm6GaDj5sja9CRte1R6feBP0GZt8mZ0gdZ+WEEIIkSKRqMLW8gbCHVyRV1VV/rl6DxsOuLCZjdw1dzQFmambSt3XaWewXqvxgtbyUrlFm1KdrB1LtPVgAKZeAcNnJ19mJ5GWGCGEED2Kqqpsr/DgD3V8JtJbXx1k2bYqDAa45YwRDNEzwDiC7gEMQPUO8NclX07JGljzJ+3xuPkw/jvJl9mJJIgRQgjRo+yu9uLyd7wL5ZOd1bz8RSkAV84czNRBOXpV7SgpCWBq94CnIvlyyr7WUgqoCgw/C467qlvPRGqOBDFCCCF6jAP1fiqTSCmw8YCLv6zYBcB5E/oyZ2yRXlU7SkoCGPdBcO1Pvpzq7bD0F9paMANP1MbB9LAABmRMjBBCiB6ixhOkpKbji7ltLnPz2w+3EVFUThiSy/dPGKhj7eL1y7YzKE/nAMZbA7W7ky+nvhQ+eggifiiaCKfeCUadVg3uZNISI4QQottrCITZmcRU6q3lbn7z362EogqTi7O58fThGFPU8pCSACbggqqtyU+l9lTCkvsh6Ia8EXDGfWCy6lPHLiBBjBBCiG4tEI6yvaLjU6l3VDTw6/9qi9lN6O/kR7NH6pdw8QgpCWBCPqjYrI1dSYa/XgtgfNXgLIbZPwWLjlO+u4AEMUIIIbqtcONU6lCkYxHMrioPD3+wlUBYYWzfLP5vzkis5h4UwESCWjJGJZJcOSEffPQguA9AegGc9XOwZ+lTxy4kY2KEEEJ0S+GowpYyd4enUu+p9vLw+1vwh6OMLsrkzrNHYTOnZuxH/2wHA/N0btWIRrQAJhJIrpxIEJb+DGp3gd0JZ/0C0vP1qWMXkyBGCCFEtxOJKmwta+hwVup9NV5+9f4WvKEoIwszuOvs0dgtqQlgBuQ49F2JF0BRtNV4k13MLhLUMlJXbNS6jmb/DJypS24Zx+KArOKUvoQEMUIIIbqVqKKytbwBT7BjXSiltT5++f4WPMEIwwrS+fHc0Tis+gcwBgMMzU+nT5bOSRJVVZsCHXAlV07Yp2WkLt8AZhuc8QDkDdOnjm2xZUCfceBLshWpDRLECCGE6DaiisqWMjcNgY4FMAfq/fzy/S00BCIMyU/n7nPGkGbV/6vOZDQwok8GOek6z+xRVW01Xm9VcuUEPfDxg1C1TWsROfMhKBynSxXb5MiBPmMap21LECOEEOIYoCgq28obOhzAlLn8/OK9zbj8YQblpnHPOWNSktDRYjIwqiiTTLtF34KbWmA8lcmVE3Bps5Bqd4MtU+tCyh+hTx3bklGovVYnLZwnQYwQQogupygq2yoaOpxOoMId4BfvbaHeF6Y4x8E9540hw67/V5zNYmRMUZb+3VOqqrWaJNsC46uBxfeDqwTs2TDnF5AzWI8ati27uPNeq5EEMUIIIbqUqqrsqPRQ7+tYALOn2ssj/92K2x+mf7aDe88bS5berSRAus3E6KIs/ado6xXAeCph8b3QUAZpeTDnl+AcoE8dW2MwQO4wyOqb+tc6ggQxQgghukxTAFPrDXXo/I0HXPxuyXb84SiDctP4yTmjcTr0D2CcDgujijIxGXXuJlFVbSVeb3Vy5bgPwOL7tEAoo1ALYDJTlxcqxmCEgtGQnpf612qGBDFCCCG6hKqq7KryUOPpWACzZlcNf16+k4iixhayS8Ug3oJMK8MKMjDoPc5DrwCmbh8suQ/8dZA1QOtC6ox1YEwWbQCv3Zn612qBBDFCCCE6nRbAeKlq6FgAs3hTOf/6dC8qcMKQXG44bXhKVuJNySJ2oAUwlVu0MSzJqNkJSx7QciHlDNZW4nXk6FLFVpnt2mwna9emLZAgRgghRKdSFJWdHWyBUVWVV7/cz5vrDwBw1thCrjxxMEadu3mMBhicn06h3mvAgLaQXdXW5AOYyi1aNuqwV0vmeNbPtNlIqWZ3ai0wJv277RIlQYwQQohOE4kqbKtowO1PfBp1VFF5avUelm3TpiAvOG4A357SX/duHqvZwIjCzJQMDtYCmC3gq02unL2rYPXvIRqCPmO1dWA6o1Ukqx/kDu20KdRtkSBGCCFEpwhFFLaWuzuUSiAUUfjT0h18sa8OgwH+30lDOHNMoe51zLSbGVGYkZocS3oEMKoKG16B9c9qzwdMg1N/DJYUtBgdzmDU1n/J6JPa10mQBDFCCCFSLhCOsqXMTSCsJHyuJxjh0Q+3sa2iAYvJwM1njGDa4Fzd69gny8aQvHTdu6YAiIa1XEgBd3JlrPkT7FqqPR9zIRx/dePKuClktkOf0Z3TVZUgCWKEEEKklCcYYVu5m1BETfjcWm+IRz7YQmmdnzSriTvmjGJM3yxd62cwwJBUjX8BCPu1bNRhf8fLCLhg2a+gcpPWKnLCdTDqXP3q2JJuNP6lORLECCGESBmXL8y2igaiSuIBzPaKBn6/ZDv1/jDZaRbuPmcMA3XOFp3S8S+gtbxUbtZaUTrKtR8+/qm2iJ0lDWb9BPpP1a+OLelm41+aI0GMEEKIlKjxBNlZ6aED8Qsfb63g6U/2ElVUinMc3Hn2KAoy9W0pSen4F9DWf6naBmriXWgxZV/D8l9ByKstYnfGA5AzSL86Nqebjn9pjgQxQgghdFfuCrC3xouaYAATiSr869O9fLxVm4F0wpBcrps1DLtF30AjpeNfQGs9qd2TXBk7FsOaJ0CNaqvinn4fOLJ1qV6LzHat+8iWkdrX0YkEMUIIIXRVWutjf13i4z/qfSEWfbSDbRUNGICLpxVz4aR+uk6hNhkNDM5Lo0+qxr+oqpY92n0wiTIU+PJfsOkN7fmQWXDSrWCy6lLFFqXnQ97wbjv+pTkSxAghhNBFVFHZXeWhugOL2O2s9PD7j7ZT6w2RZjVx8xnDmVys78qz6TYTI/pk6p+BuokS1bqPklnELuiBT34PpWu155MugUnfT+24FKNJG/vSGbmWdCZBjBBCiKQFwlG2lTfgCyW+BszybZU8tXoPEUWlf7aD/5szkr5Oh67165dtpzgnLXXdR5GQNnMo6Ol4GdXbYcWvwVMBRrPW+jL0dP3q2BxbBuSP6vL0AR0lQYwQQoik1PtC7Kz0EI4mNgAmoig8u2YfizdXAHD8oBxuOG24ri0lVrOBYQUZZKelsCsm5IWKzRAJdOx8VYWt78AX/wQlog3gnfUTbXBtKjn7Q/ZgMOqfc6qzSBAjhBCiww7U+ymt9SU8gNflD7Poo+1sLW8AtBQC86b0x6hjt0lOuoVhBRlYTCn8kvZUakkYlcRboAAIeeCTP0LJp9rzgTPhpFvAmsKBtSYrFIzsnESRKSZBjBBCiIQlM/7lm/31PLl8F/X+MA6LiRtPH85xg/T7QjUaYFBeOkXOFC7FryjaAN6Gso6XUbMTlj8CnnKt++j4q2H0+akd/5KWp7Xw9KDBu61JOIjxeDy88sor7Nu3jxEjRnDxxRdjtbbcTLds2TI+/PDDuG0Oh4MHH3ww8doKIYTocoFwlO0VDQnnQApFFF76vIQPNpYD0D/bwY/OGkn/bP3Gv6RZTYwozCDNmsL/0cMBLQdSR8e/qCpsfRe+eOqw7qMfQ/5Ifet5OINRG7yb1Td1r9EFEvqUa2pqOOmkk8jIyODMM8/kV7/6FX/84x9Zvnw5aWnNDwr65JNPeOGFF7jhhhti2+z2FCeqEkIIkRIuX5gdlQ0Jj38prfXx+LKdlNT6AJgztpBLTxiE1axPV4/BAEVZdgbmpnDwLmjJG6u3d3wF3pAXPv0T7FutPS+eASfdltp1WWyZWoDUQwfvtiahIOaXv/wl0WiUlStXkpaWxl133cXIkSN54oknuPPOO1s8r1+/fvzkJz9JurJCCCG6zsF6PyUJjn9RVZXFmyt4fu0+wlGVLIeFa08dytSB+nUfpVlNDC1IJzNVqQNAaz2p3wf1pR0vo2anNvuooQwMJq37aMwFqes+MpogZzBk9u3WqQOSkVAQ88Ybb/D9738/1uqSl5fH+eefzxtvvNFqEFNVVcUvf/lL7HY706dP55RTTkmu1kIIITpNKKKwu9pDnTex1od6X4i/rtzNV6X1AEwuzubaU4fqNlPIaIB+2Q76ZztS2/oSDUPVVvDXd+x8JQIbXoVvXtYepxdos48KRulazThpudrCdWZb6l6jG2h3EBMMBtm3bx/Dhw+P2z58+HDeeeedVs/NysrC7XZTUlLCgw8+yAUXXMDzzz/f4iqMwWCQYDAYe+52J5G6XAghRIfVeUPsrvYknIF6fUkdf1m5G7c/jMVk4NITBjFnbKFuq+9m2s0MLUhP7dgX0BI4Vm2FSLDtY5tTuwc+WQS1u7TnxTO09V9smbpVMY7JArnDIKMgNeV3M+3+9H0+rR8zKys+BbrT6cTr9bZ43ve//33uvffe2A/u9ddfz/Tp0zn33HP5wQ9+0Ow5Dz/8MD/96U/bWzUhhBA6UxSVfbU+yl2JrX0Siig8v/bQ2i/FuWncfPpwinXKPm0yGhiYm0Zhlk3XdARHUVUt/1H9PhKePw5Ht75YM+CEa2HIaanr2snoow3e7SUzj9qj3UFMeno6BoMBl8sVt72+vp6MjJYHJA0dOjTu+eTJk5k6dSorV65sMYi5++67uf3222PP3W43xcXF7a2qEEKIJHiDEXZUevAnuPruljI3/1i9m4P1WuBzzvgivjdtoG6Dd3PSLQzJT09d1ukmYb82eDfQwV6Aur2w+vfxrS8n3pi6dVnMdsgbpnUhHWPaHcRYrVaGDh3Ktm3b4rZv27aNMWPGJPSi0WiUQKDl6N5ms2Gz9e5+PCGE6G5UVeWgK8D+Wh9KAo0PnmCEFz8rYWlj5ulsh4VrZw1jcnG2LvWymAwMzk8nP6MTvhfcB7UgpCOL1ykR2PAafPNS57S+GAzaoN2cwdog3mNQQuHxggULeOWVV2JjVMrKynj33XdZsGBB7JiPPvoorito9erVcWV8+umnrF+/njPPPDOZegshhNBRMBJlc5mbkpr2BzCqqrJmVzV3vPp1LIA5c3Qffrtgki4BjMEARU47k4qzUx/AhANQvgFqdnUsgKnbC+/9H3z1nBbAFJ8AF/5Zy32UigDGkQ39pmgtMMdoAANgUNX2d/a5XC5OO+00AoEAp556Kh9++CFDhw7lgw8+iLWcPPTQQyxatIj6+noALr74Yvbv38+UKVOorq7mrbfe4vLLL+cvf/kLxnbma3C73TidTlwu11FjcoQQQiSn2hNkT7WXSAJrv1Q1BPnnJ3tiM4/6ZztYeMoQRhfp8zc6O83CoLy01A/cBWgo11bf7UjwEg3Dxtc6b+yLxQE5QyA9T/+yUyDV398JBTEAoVCId955h5KSEkaMGME555yDyXQoCly9ejVffPEFt912W2zb+vXrWbt2LWlpaRx//PGMHTs2oUpKECOEEPoLRRT21nipSSB1QFRR+e/Gcl79spRgRMFsNDBvSn8umNRPlxxFdouRwXnp5KSnMGFjk0hQW7vFV9ux8/d/AZ//TeuCAq31ZcaNqRmbYjRDdjFk9utRCRu7XRDTFSSIEUII/aiqSoU7SGmdL6HWlz3VXv6+ajd7qrUZqaOLMll4ylBd0gaYTQb6ZzsoyrKnds2XJp7Kxq6jSOLnusvg87/D/s+05/ZsmPb/UtP6YjBoaQmyB4G5EwI7naX6+1sSQAohxDHEG4ywu8qLJ9j+L29PMMIb6/bz303lqCqkW01cesIgZo0qSDrrtMEABZk2inPSdJvF1KpwQOs68tUkfm4koE2b3vgGKGFt1d0xF8CkS1KzpL89S1vzJZUpCXo4CWKEEOIYEFVUSmt9lLsD7V72JKIofLS5gtfXHYgFPTOH5XHZjEG6rLqb5TAzOC+ddFsnfBUpCrhKtbVfVCWxc1UV9n2iJWz0Vmnb+k6G6ddqXTx6szi0GUfp+fqX3clS3dkjQYwQQvRytd4Qe2u8BMPt+/JWVZUv99Xx/NoSyt3achgDchxcNmMQEwdkJ12fDJuZ4lyHbukH2uSr1Vpfwv7Ez60vgc/+CmVfa8/TC2DaQhg4U/+uI7NdC4oyCnt8riNVVan0VVLmKkvp60gQI4QQvVQwEmVvtY9ab/sH7u6u8vDc2n1sKWsAIMth4eLjBnDaqD6YkhyrkmY1UZybRm5nDNoFLWip3d2xgbvBBm29ly3vghoFowUmfAfGX6QFG3oyWRuDl6IeNWi3Ja6gi9KGUvwRP0YltdcjQYwQQvQyiqJS7g6wv85PtJ2LvtR4grz8eSmrdlYD2gJz503oywWT+uOwJrcOicNqYkCOo3MWq4Pkuo7Cftj8Fmx6A8Jauh2KZ2itL5lF+tbTZAHngB4346glgUiA0oZS6oP1nfaaEsQIIUQvUtWgzTpqb9dRIBzl7a8P8t43ZYSi2jknD8/nu9OKkw46bBYjA3IcFGSkOM/R4bw1WutLJLGcT0RDsO0D2PAKBBrT6+QMhuOugv7H6VtHoxmc/SGrf69YqC6qRDnoPUiFtwKVzp3wLEGMEEL0Ai5fmJJaX7tnHQXCUT7aUsE735Th9ocBbcr0D2YMYlhBcrNhrGYteOmT2YnBS9CjrZrrr0vsPCUKuz6Gr188NGg3sy9M+QEMPgUMOraQGE2Q1Q+yBoCpd3z9Vvur2d+wn7AS7pLX7x3vohBCHKN8oQgltT7qvO37EgmEoyzZXMG73xzEHdACnsIsG9+fPohpg3OSCjrsFiP9srWWl05Z6wW07p+6fYcCkPZSFW3G0frnwH1A25aWp02XHj5bay3Ri8mqBS+ZfXtN8NIQaqCkoQRfU5dbF+kd76YQQhxjgpEo++v8VDUE2zVlOhCOsnhTOe9uKKOhMXjpk2lj3pT+nDIiH3MSYzIy7Wb6Ou3kpls7r+UlEtRmDnkqaPeccdCOPfCFFrw0ZZm2ZcGEi2H0uVrAoRdLmtZtlN6nV4x5AfCGvexv2I871MEM3zqTIEYIIXqQqKJysN5PmSvQrkG7/lCUxZvLefebslhXU2GWjXmT+3NyksFLTrqFvk4HToelw2UkLBrWBu26yxIbtKtEtZaXja8fCl4sDhg3H8ZcqO9idbZMbcBuWl6PnyrdxBf2sd+zH1fQ1dVViSNBjBBC9ACRqEK5O0C5K0C4HakCfKEIizdV8N6GQ8FLUZadeVP6c/Lw/A5PlzYYID/DSr9sR+ckZ2yiRLXZRu6DiaUKiARh58fabCNPubbNbINR58H472ir4uolLVcbrOvI1q/MLuYL+zjoOUhdMMGxRo1SPdBXghghhOjGwlGFcleAcnegXXmOajxBFm+u4OOtFXiDWlbmvk47357Sn5nDOh68mE0GCjJs9M22YzN34owaJaplmXaVaq0w7RX0wLb3YcvbEKjXttmyYPS3tJtewYvBCBl9tDEv1nR9yuwG/BE/Bz0HqQ10LDlmKBpizcE1fLjtQ51rFk+CGCGE6IZCEYUyl58Kd7Bd3Ua7qjy8v6GMtbtriTaOEenntPPtqQOYOTSvwwNt020mCrPs5GfYkl7sLiHRsNbq0lCWWPDirdbWedn+X4g0rtCb3gfGfRuGnwUWnRaqszi0gboZhb1msC5oa70c9BykJtCB3FJAXaCOZaXLWF66HE/YQ9Qf1beCR+g977wQQvQCgXCUMleASneAtmKXqKLyxd5a3t9YxvYKT2z7mL6ZnDu+L1MH5nQoeDEaIC/DSp8sO1n2ThzvAlqCRvcBbcCuksAXYO1urdVl9/JD3U05g2HcRTDkFH1mGxkM4MiFrL7gyEm+vG7EE/JQ7i3vcLfRXtdeluxbwmflnxFVtc8tz57H7IGzuZ/79axqHAlihBCiGwiEoxyob99sI18owtKtlXy4qZxqj5ZSwGQ0MHNYHueM78uQ/I51a9gsRgqz7PTJtGExdfJsmpBXG/PirWr/bKNoGPathq3vQ9WWQ9sLx2vpAfofr8/AWpMVMgu1lhdzJ6063AlUVaUuWEe5txxv2Jvw+YqqsK5iHUv2LWFH/Y7Y9pE5Izlr0FlMLpiMJWSRIEYIIXqrel+IcneAel+4ze/u/XU+Pt5SyfLtlQQaV+TNtJs5a0whZ40t7FBCRYMBnA4LRVl2stMsnTdFuom/XgteElmkrqFc6y7asRiCjVN9DSYYeCKMvRD6jNGnbnanFrik5/eaWUagrbBb7a+mwldBMBpM+HxPyMMnBz/h45KPqfZraSpMBhPTi6Zz1qCzGOwcrHONWyZBjBBCdLKoolLVEKTcHcAfar3LJBCOsmZ3Dcu2VrKj8lCX0YAcB+eO78tJw/OxmhNvNXFYTRRk2sjPsHbuQF3Qchv5qrUxL8GGdp4ThYPrYOt7cOBLaJr1kpYPI8+GEWdrs4OSZXFoA3XT++g3fqabCEfDVPgqqPJXEUlkhhdaq8uW2i2s2r+KdRXriKja+RmWDE4rPo3Ti08nx975XWwSxAghRCdpGu9S7Qm2OtNIVVV2VXlZtq2ST3dVx1pdTAYDUwdlM3tMIRP6OxNuNbGYDORlaIFLZmePdQFtdd2Gcm28S3sH6/rrYOdHWl4jb+Wh7f2mwKhzYcD05PMPGc1aa0tGH631pZfxhDxU+iup9dcmPOW51l/L6oOrWX1gdazVBWBg5kDOGHgGM/rOwKrnAoEJkiBGCCFSrM57qMuoNQ2BMKt3VrNsWxWltYeWcy/KsnP66D6cOiI/4S4jgwFy0qzkZ1jJSbN2XjqAJqoKvlptllF7u4wiAShZC7uXwsH1hxa1s2ZoKQFGnaOtx5IMgwHs2drsorS8XrOibpOIEqHGX0O1vxpfJLHUABElwleVX7HqwCo2Vm+MBT4Os4MZfWdw6oBTGZQ1KBXVTpgEMUIIkQLBSJSqhiBVDcFYS0pzIorChv0uVu2s5vM9tUQapyRZTAZmDMnj9NF9GF2UmXCrS4bNTH6mlfyMLhikCxAJaYvLNZRrC861RVWgYiPsWqqtrBv2H9pXMApGnqMlZEx2YK0tU2t1Se8D5q5rQUiVhlADVf6qhFtdVFWltKGUNWVr+PTgpzSEDnXzjcoZxakDTuW4wuO6tNWlORLECCGEThRFpcYboqohiDvQ8kBdRVHZUu5mza4a1u6pjcs8PSgvjTNG9+GkYfmk2xL7E51pN5ObbiU33Yrd0snjXEBrdfHXad1Fvpr2zTJylcKuZbB7WXwSx4xCGHo6DDs9+VaXpsAlLb/XjXMBCCvhWKuLP+Jv+4TDlHvLWVu2ls/KP6PMWxbb7rQ6Oan/SZzS/xQK0wv1rrJuJIgRQogkuQNhqhqC1HpDLY510ca5ePh0Vw1rdtfEdS05HRZOGJLLaaP6JDw9ussDF9AG53oqtYXmoqG2j/dUwr5PYc8KqDk0NRdLOgw+WQtc+ozVVsPtqF4euKiqijvkpsZfQ20gsVaXan81n5d/ztqytZQ0lMS2m41mJuZP5KT+JzExfyKmZMcadQIJYoQQogOCkSjVHq3VpaUZRqqqUlrnZ82uaj7dVUNlw6FulXSrielDcpk5LJ8xfbMSWg03024mr3GMS5cFLuGANtDWUwXhdoy5cB/QApd9n8YHLgYj9D8Ohp2hDdLtaHeRwaCNmUnPh/SCXrWey+EaQg3U+GuoC9YlNMPIFXTxefnnfFb+GTvrd8a2mwwmxuaNZXrRdKYWTsVhduhWV4vRQpqeiTWbIUGMEEK0UyiiUOsNUeMN0hCINNtboigq2ysa+LKkji/21lHuDsT22cxGjhuUw8xh+Uwa4MTczrEqJqOBLIeZnDQr2WmWzp8S3SQa0bp8vJUQcLd+rKpC/T4taCn5FOr2HrbTAIXjYNBJ2jiXjiZMNJq0lXMdudr0alMXzLjqBL6wj5qA1uISak9LV6MKbwXrq9bzVeVX7KjbEWutMWBgVM4opvedznGFx5FpzdStrmnmNJw2J9m2bDKsGbjdbfycJEmCGCGEaEUwEqXOG6baE8QTbD5wCYSjbNjv4suSOtaV1NEQOPQfstloYHJxNjOH5TFlYE67W05sFiPZDgs5aVacDkvnzypqEg1rs4t81drCdGrLg5RRFajeAaX/04IX94FD+wwm6DtRC1yKT+j4sv0Wh3ZuWq42u6gXLUJ3uEAkQG2glhp/DYFooO0T0NZy2V2/Oxa4HD7GBWCocygn9D2B4wuP121NFwMGsmxZZNuyybZld/rAXwlihBDiCMFIVGtx8YTiApLD1flCrCup48u9dWw86CJ82FiYdJuJKcU5HDcoh4kDnKRZ2/5TazBoM4qy07TAJdFBvbqKBLWBud5qbUXc1gboBtzaNOgDX2iL0QVch/YZLdp6LoNOguLp2jiVRBkMWvbptFytxSXF3RNdyRv2Uh+spz5Q3+5p0cFIkM21m1lfuZ6vq76Om1VkMpgYlTuKyQWTmdxnMvmOfF3qaTPZyLJqgUuWLQtjMmOXkiRBjBBCoOUjqvOFqfM2H7iEowrbKxrYcMDFN/td7KmOzzXTJ9PGcYNyOH5QDiOLMjG3Y90Ru8WI02Ehy2HB6bB0zVToJiGfFrj4qiHoafk4VdGSLe7/QgtcqrfHt85Y0hoDl5nQf1rHgg5rmtbK4sjRFp/rAQNMO6JpcG59sJ66QB1hpe0FAFVV5YDnAJtqNrGpZhPbarfFnecwO5iYP5EpfaYwPn88aZbkgz6z0UymJZMsWxZOmxObqfuMN5IgRghxTFIUFZc/TJ0vRL0/TPCItVxUVeVAvZ9v9rvYcMDFljI3wUj8McMK0jl+UC7HDcphQI6jzbVcLCYDzsaAJcth6bpBuaC1rgRc2pRof60WxLTEXwfl38CBxhaXQH38/pzB2uDc/sdreYsSzRhttjUGLdnafS9cv6VJRIngCrqoD9bjCrpiGZ9b4wq62FyzmY01G9lcsxlX0BW3P9+Rz+SCyUzpM4UROSMwJ5mx24CBTGsmWdYssmxZpJnTOj+nVjtJECOEOGYEwlHqfWHq/SFcvjDKEb0k9b4Qmw66G1tb6qk7YoVdp8PCxP5OJgxwMqG/s83Vc01GA5l2cyxw6dIuItBmFPnrtFugXstH1JyAW1t4ruxrKN8ArpL4/WY79J0MA47Xgpf0gsTqYbKCPUtrZbFn9+ouIlVV8UV8uIIu3CF3XHdPS4LRIDvqdmitLdWb2O/ZH7ffarQyMmckY/PHMj5vPP0z+icVZBgwkGHNIMOSQaY1k0xrZpd2ESVCghghRK8ViSo0BCK4/GFc/jC+w6ZCq6pKtSfI1vIGtpQ1sLXcTZkrfgClxWRgTFEWEwY4mTggm+I2WlusZgOZdguZdjOZdgvpVlPX/gerKFqw4q/TBuW2NBU65G0MWr7RWlzq9sKR647kDIG+k7TApc+4xGYCWRzauJamwMWi3zTe7igUDcWCFnfI3eZUaF/Yx476HWyv3c72+u3sde09qoVmUNYgxuZpQcvw7OFYkpiJZTKYSLekk2XNIsOaQbolvccELUeSIEYI0WsoihoLWtyBcNxsIlVVOVgfYEu5m63lDWwtc1PjjZ+uagAG5qUxob8WtIwqzGw1Q7TdYiTTbiGrMWhxWLt47IaiQKhB6yYKuLQWlSNnE6mqNk26agtUbtXua3cffZxzoDabqGgiFI7XApD2MBi0cTF256HApZeu2dIkqkRpCDXgCrlwB91tziZyBV1sr9seu+1v2H/UYnU59hzG5Y1jXN44xuSNIcvazve/GRajhQxLBhnWDLKsWTjMbXd99hQSxAgheixFUfGEIrh8jUFLIBLrIvKFIuyu8rKzysOuSg/bKxpwHzFg12QwMLQgndFFmYzum8WowswWu3xMRgPpNhOZNgsZdjMZNnOrAU6nUBRt9lDA1XjfTNASDWtBStVWqNyiBS2+mqPLyuoPRRO0oKVoQvunQJtt2iJztsxDt146ELdJWAnTEGrAE/LQEG7A18pifxElwgHPAfa49rDbtZsddTuo8FUcdVxhWiEjc0bGbvmO/A4FGkaDkXRLOumWdDIsWitLd8t3pCcJYoQQPUa4sXvIE4jgDoTxBrWgJaqolNb52Fnpid0O1vuPWojdYjIwok8mY/pmMrooi+F9MpodXGswQJrVRIZNC1Yy7GYcli7uGgJtDEtTsBJwacv9Hx60qKqWcLF2p7ZeS9VWqNl5dCoAgxFyh0Gf0VAwRlviP70d02+NpsaAJUNrZbFl9vpWFtDGqDSEGmKBS0stLaqqUhOoYXf9bna7tFuJu4SQcmSLn4EBmQPighanzdmhuqWZ00izpMUClt7UytIeEsQIIbotXyhCQ6DpFiYQVogqKgfr/eyt8bK32suuKi97qr2EokcvwpafYWV4nwyGF2QyojCDofnpR62S2xSwpFnNpNtMpNvMZFjNXbe43OFCXi1QCTZowUvYz2H9Y1qixZqd2jL+NTuhZheEmpkebcuEgtHazKGCMZA3ou18QkYzWNO1gMWaoT22pPXaxeWaKKqCN+yN3TxhT7Or5KqqSl2gjn0N+yh1l7LHrbW0NDdw12F2MNQ5lCHOIQx1DmVkzsiEpz4bMJBmSYsFLU0BS08dy6IXCWKEEN1CIBzFG4zgDUbxhiJ4glqLS0mtj301XvbW+Nhb7aW0zhe3sFwTh8XEsD4ZDC/IYFifdIYXZBw1e8hsMuCwaC0saTYT6VathaVbBCzR8KFgJdigrdXSNCBUiYBrvzbgtm5vY8Cys/mAxWjWBuHmDdMCl4IxkNWv9eDDbNOCFGs6WDMbA5belzTxSKqq4o/444KW5haZiypRyn3llLhLKGkoocRdQmlDKZ7w0e+/yWCiOLP4UNCSPZTCtMKEgg2TwRQXsKSZ0465Fpb2kiBGCNHpjgxY3P4wB+r8lNb52V/nY3+dn301XspcgWZz89otRgblpjMoL42hBekML8ikb7YdY+MfeYNBC2rSrCYcja0saVZT167LcrhIUAtSQh6ttSXk0bapqjZepSlYqW+8d+0/FNAc7vCAJW+41sKSPbDlmUNGszad2ZLeeJ+mBSy9NOfQ4ZoCFl/EFxe4KId1x6mqSn2wnoOegxzwHOCg5yClDaXs9+xvdiE6o8FIv/R+DMwayKCsQQx1DmVg5sB2zxwyGozYzXYcZkcsUHGYHb16DIveJIgRQqSMoqj4wlF8oQj+UBS3P8LeGg97q7VApSlgOejyN9u6ApDtsDA4XwtYBuelMzgvnT5ZNowGAwYD2C2mIwIWU/cYvwJaUBL2HwpUmoKWSFBb0t9VqgUorv3a47q9zbeugDYtOWcIZA86FLRkD2o+ADGateMtafHByjEwfgW0lhNfxIcv7NPuIz78YX9sBpCqqrhCrligcnjQ0tJy/zaTjeLMYgZlDYrd90vv166AJRasmByxoMVh1h6L5EgQI4TQRSAcxRfSWljK3QF2VXrY29iaUuYKUFbvp9wdaDFYsZqM9M9xMCDHwYCcNAbmpjE4L43sNCsWkwG7RWtJcTQGKXaLEbu5m3QFgRasNAUsYZ9276/XkiA2BSlNAYt7vxbINMdghKwBkDNIC1qa7tML4ruEDAZt0TmLQ7uZHYcCl1684u3hFFUhEAngj/jxR/wEIgF8ER/BqPbeBiNBKnwVVPgqKPeWa4+92mNvxNtsmQYM9EnrQ/+M/vTL6EdxZjEDMwdSkFbQZpeQ1WTFbrJjN9tj99KykloSxAgh2k1RVPyNLStlrgD7arzsq9FaUw7U+ylvDFg8wZYX97KYDPTL1gKVAY1BS3FOGn2z7Y3BiQm7uTFIaQxYujSn0OGaWlYigcZAxQfeSm0Ks+sANJRpN/dBbZaQr7rlsoxmyOwHzgGHbjmDwVl8qHWlKVAx2xqDFHv8fTvyM/UGiqrEApWmYMUf8ROMBglGglT5q7SbryoWsJT7yqkL1LVYZlOw0i+jnxawpPejX0Y/+qb3bbV1pSlQsZls2M3afdPN1MunlndHEsQIIeKoqkowolDvC7OvRpv5U1Lro7TOR1l9gMqGIJUNAQLho2cDHS4v3Upfp52+2Q7t3mmnODeN4pw00m1mLUixmLCZtftuE6iA1krS1LISqIe6fVpXj3s/NFRogUtDuXY7Mo/QkSzpkD0Asoobg5XG+8xCLZAxmhuDFHtjcGI/1MJitvf62UBNFFWJBSWBaIBAJEAwGsQT8sQFKVX+Kqr91bHH7pC71XLTLekUpRdRlFZEYXohRWlF9EnrQ1F6UbMtJCaDCZvJhtVkjQUndrM9Frx0i25KESNBjBDHoHBE4aDLz74aHyW1Pg7U+TjoClDu0oKUak+w2UzOR8pNt9In00afTBt9sx0U5zpi41acDgs2ixGbyYTNYsRqMnafrh9FaWxN8YOnXAtSXKVaC4r7oDZ12VMBnkotOWJb7E7I7HvoltVXa2XJ7AtpuYcFJ7ZDQYrJ2nh/7PwZjigRQtEQgWhAu48EqA3UcsBzgHJvObWBWmoCNdT6tfu6QB11wbq4wbfNSTenk5+WT4GjIBaoNN1nWDPijrUYLVhMFmxGW1ywYjVZsZqsSSdPFJ1LPi0hehl/KMqB+sYunjo/Ze4AFa4AFe4AVZ4QNZ4gtd4QkSOzHzbDYTFpQUqWjSKnnf7ZWtfPwLw0inMdZNosWM1G7WYydp//UpWotiBcLDA5AO4yaDigtaR4KrRZQN5qaGbWyVHMNkgvhIw+kFEIGQXavXMg5A4CR552jMl6xL3tmOnygUNBSigawhv2Uu4tp8xbRpWvikp/JbWBWuqCddQH6qkP1lMXqGtziX7QWkfyHVqQ0hSsFDgKKEjT7pvWXGkKUKxGLSCJ3TfeLEbLMb+uSm8jQYwQPUA4qlDrDVHVEKTcFaDc3RiUNASpaghS7dUCk1pvCG+whczERzAYICfNSkGmjaIsO4VZdvo6bfR1ai0qxTlp5KRbsZlNXb+8fhMlquX9cR/Uxp54KrXgxFOuBSbeKi0w8ddqa620lz1bW7E2LV+7z+yrLcPfNE4lPf9QUGKyaEGKyXrMBSjBSJCaQA3l3nKq/FVU+iqp8ddQG6ilNlAbS3pYH6ynIdRwVD6glmRYMshz5JFrzyXPnkeuo/HenkueI488ex42sw2LUQtQLCbLUY8tRkv3CaJFp5EgRoguEAhHqfeFqWkMPmo8QaoaQlR5grGWklpviDpfmHpf6KicP22xmY3kpVvJz7RRkKF19xQ6tUCln9NBvxytVSXN2sV/AiIhLeDw1WgBiacKfJVaIOKp0gbGNgUlvjoIuo7ODdQao1nLAeTIbQxOirQWlMx+WpePc4AWrFjTDwUmJssxMw4lokSoDdRS5dPGmVQHqqn1awFJU1eOO+jGFXLFlt0/Mrtya0wGE9m2bLLt2WTbssmx5ZBty44FKX3S+lCYVkiGNQOL0YLZaI61phz+XIiWSBAjRAcpikpDUFuorSlrstuvPa/3h6jxhKjzacFIvT9MvS8cS1QYjCTwRdzIaIBsh5XcdO2Wn2ElL9NGnwyb1pritFOUZadfjoNMm7nz/itVVW18ib9Om1IcqNce++q0IMRXC4G6xvt6bXvTMa0kzmuVLasxOMmBtDytmye9jzZYNqNxTEpWP22b2dqrExI2LeLWEGrAHXLjCrqoDx7qrqkP1lMfqMcVcuEKHgpGPGFPi2uitCXdkk6WNQunzYnT6iTHnkOOPYdcey659lzyHfkUphWSa8+NjTNpCkrMRrN06QjdSBAjjknhqII3qC1t7w1G8QTDuBtz9Hga8/Q0BLT9Ln849rwpWGk6rn2N5c0zGsDpsJCdZiU3zUJOupW8dBt5GVbyMqwUZNjJz9S2NQUuJr0HxipRbXG1pmXum5a9j21r0BIN+usa7+vjsyY3HdPcarLtZTBqQYk9GxzZh4KT9HwtCMkoaAxQirRbep9eMxhWURV8YZ+2emzEizfkpSHcgDfsxR104w65Y8GJJ+TBE/bEApCm576wj4iaxPuPFpRkWjNxWp1k2bLItmXjtDm1VhN7Ljm2nEOtJ44+2C12LSAxdGKwLEQzOvSX4JtvvmHfvn2MGDGC0aNHp+wcISJRBX84qt1C0cY1SrTHvpC2XknTAmv+UBRvKIo/FGm8j2r5dxpz8HhDkdhS980lC+woq9lIlt1Mpt1Clt1MlsNCdpoWnOSlWclJt5KdZiEn7dC9M82SWGuJokDY2zjt16utTxL32Hdo3ZLDl7IPeQ4tbx/0aMcHD9vf0ZaQ5hhMWqJBR3ZjUOJsDEiywZ6jzdJx5GjjTjIKtBYUR44WvHTzsSVNi6o1Tf1tWqekaQn7phVhY88bA5Om5e2bgo24VWTDvnYNam0vo8EYy2Scbkknw5JBli0Lp9UZazE5vFsn16EFJzn2HFnfRPRYCQUxoVCIiy++mFWrVjF58mQ+//xzFixYwD/+8Y8W/xh35BzR/UQVlVBEIRRRCEajhx433cLRQ48jUYLhQ48DYYVAOEqgcXsgHG28KQQihz0OHwpSmh5HWljdVS9WkzGWuTjTbibDZibLbiHLYSHT3rTt0OMsuxmnRSHbqpJpiZBpVrEThmhQ61KJhCDi09YZiQQO3YJB8Pgbt/sP7W9aiyTs17aHGu/DvvjtLa3uqhejpTFjceahrMW2TC0ZoMOpBRr2bLBnHbpvClSaHtsyUz6WRFVVwkpYmwGjhGIzYYLRYOy+tceBaIBgJP5xIKqtR9K0LklTgBILWBqnA6eS0WAkzaxlJk6zpMWCkab7TGsmGZYMMm2ZZFoytaDE5iTTmkmWNYssa5YkCBTHpISCmEWLFvHJJ5/w9ddfM2DAADZt2sTxxx/PaaedxmWXXabbOb2VqqpEFZVI0y2qEFG0beGocmhfVCWiKLH7cFTbFm7aFlUIN51/2PZwVCEUPfQ43HgfiSqEYtu04CPctC2indO0PdR0H4l/3p7puAm8E5hQMBM97D6KufGx3RAho/GxxRDBhIKFCOkWSDerOEwq6WaVdItKuilKWuM2h0nBYYxiNyrYTAp2QxSbMYrdGMVmiGI1RLESwdJ4MxHFpDQGINGQFoBEQ+APQUPjtmhICyCaHqf4y6xdmpaWj9035sVpum8KQmJZiY9+HLU4iFrSiFjTiFgcRI0momqUiBIhokTiHkfUyKHHcTcfYb+LiPfQtrASbv4+GiaiavdhJRwLRJoeN7f98EAlHA0TUkLNJuHrbIcvgJZmTovdO8wO0izafbo5PfY4zXIoOGnanm5Jj3tsM9kkABGiAwyqqrb722nChAnMmjWLxx9/PLZt/vz5eDweFi9erNs5R3K73TidTv66+GusjnRUVUVRFBRVQVG0pl5UBSWqEm16rGhBAaqKokRRFBVVbTpHRVW0Y1TlsG1qNLZPbSxDbTpfVVGVaNy5NJanPdfObyqvqYwj742oGFExcOixEQVD7LH2BR+336DtNzVuMzSec+jYQ/u07Qomg9LsMaam/Yc9NqLGjj98v5kjt0UwNh1nUDEbGo8zaDejQcVENHaMkShGohhUFYMhikGNYkDBoEa1sSQGUKGxdholtk07QKHxeePf96ZjlcbzlSOObTq/qdxD52k31WA49BhQDIfKU458bmgs47A6ND1vOlYxmrWbyYxiNB16ftjjqNGEajQRNZpQDCYUo4mo0Yhi0G7adiNRgxHFaCSKAcVgIGpovG98zagBokBUjWo3JYqiKkTUiPbzftj2iBpBUZXY46gSv6+3MBvNWI3WoxYtO3wp+LhtZltsyfjDHzctH2832bGZbbHkfE0J+5py4NhMNhmUKkQCmr6/XS4XWVlZupff7paYSCTCli1buPnmm+O2T5w4kb/85S+6nQMQDAYJBg81n7tcLgA+/PJCLHZji4Mpm9t++DbVEL8tbh/x/wUdeUzs3gCYQD2iC1k9rIxmy2/nax/+OocfG3/8ofPUFso+dDMcVV7zxxqafU2t/ET+QzQ13o5FKqhhiIa1aKOHMRqMmAwmzAZt9ojZaNaeN84oibsZzJiM2rFH7jt8aqzJaMJisMTNTImbPmuyYDVoi5CZTVpA0vTYYrRoC5QZLHGLlVlMls5duEwBQhAOhQnT9S1BQvQkbreWFiKB9pKEtDuI8Xg8RKNRcnJy4rbn5eVRX1+v2zkADz/8MD/96U+P2v7G3XvbW10hhBBCdBM1NTU4nU7dy213EGOz2QDw+eJnM3g8Hux2u27nANx9993cfvvtsef19fUMGjSIkpKSlLwJ3ZXb7aa4uJjS0tKUNMN1V3Ldct3HArluue5jgcvlYuDAgeTm5qak/HYHMQ6Hg6KiIkpKSuK2l5SUMHToUN3OAS34aQqADud0Oo+pD79JVlaWXPcxRK772CLXfWw5Vq/bmKJlFBIq9ZxzzuGNN97QBrUCgUCAd955h3POOSd2zObNm3n77bcTOkcIIYQQIlEJBTEPPPAA+/fv56KLLuLvf/875513HmazOa7r55VXXuHyyy9P6BwhhBBCiEQlFMQMHjyYdevWMXr0aJYvX84pp5zC559/Tl5eXuyYsWPHcuGFFyZ0TltsNhsPPvhgs11MvZlct1z3sUCuW677WCDXnZrrTmidGCGEEEKI7kJWbRJCCCFEjyRBjBBCCCF6JAlihBBCCNEjJZQAsivs2rWLDRs20KdPH2bMmJGyueZd6cCBA3z99dc4nU6mTJlCWlpa3P6PP/6YqqqquG2DBg3ixBNP7Mxq6qq91xQMBvnkk0/weDxMnz6doqKizqymrrxeL++8806z+6ZOncrIkSOB3vN5V1ZWsnz5coYPH87UqVObPaY9v9897W9AW9cdCoX46quvqKqqYvTo0QwbNixuf0VFBcuWLTvqvPPOO4/MzMyU1TtZrV13ItfUmz7vtWvXsmfPnqPOsVqtzJ8/H+iZn3dbP8NNOuP3u1sHMffddx+///3vmTlzJlu2bKFfv358+OGHR6Ux6Klqa2u55ppr+Pzzz5kwYQKlpaVUVVXx/PPPc8YZZ8SOe/DBB6mvr2f8+PGxbTNnzuxxX2qHa8817dy5k7POOgubzUZhYSFffPEFf/rTn7j66qu7ospJ8/l8/Oc//4nbdvDgQVatWsXrr78eC2J6+uddVlbG//3f/7FixQp8Ph+XXnpps1/m7fn97kl/A9pz3S+++CL33XcfeXl5FBQUsGrVKr71rW/xzDPPYDZrf443bNjAJZdcwne/+924c0899dRu+aXWnutu7zX1ts973bp1rFixIm7b4sWLGTBgQCyI6Wmfd3t+hqETf7/VbmrZsmUqoK5cuVJVVVV1uVzqyJEj1euuu66La6afkpIS9fXXX1cVRYltu/HGG9X8/Hw1Go3Gtp100knqgw8+2AU1TJ32XNOsWbPUOXPmqJFIRFVVVX3yySdVq9Wq7tu3rxNq2DluvfVWtaCgQA2FQrFtPf3z3rZtm/rcc8+pgUBAPeGEE9Qbb7zxqGPa8/vd0/4GtOe6X3nlFfXAgQOx57t27VIzMzPVP/zhD7FtS5YsUbvxn+ajtOe623NNvfHzPlJlZaVqsVjURYsWxbb1tM+7PT/Dnfn73W3b6Z577jmmTZvGKaecAmhLNS9cuJAXXnghtvpvT1dcXMz8+fMxHJYlev78+VRXV7N///64Y/fv38+bb77J//73P/x+f2dXNSVau6b9+/ezYsUKbr31VkwmLSv2//t//4+0tDReffXVrqiu7kKhEM899xxXXnklFoslbl9P/rxHjhzJpZde2uq6EO35/e5pfwPac90LFiygX79+sedDhw5l8uTJfPXVV0cd+/HHH/P++++zd+/eFNRWP+257iatXVNv/LyP9Mwzz2A0GrnsssuO2tdTPu/2/Ax35u93tw1iNmzYENecDjBhwgTcbvdRuZh6k48++ojs7Gz69+8ft33x4sU89dRTXHrppYwYMYKlS5d2UQ3109o1bdiwASDuZ8BisTBq1KjYvp7uP//5DzU1NSxcuPCofb3x8z5ce36/j4W/ATU1NXz11VdHXafFYuGnP/0pjz76KKNHj+byyy8nHA53US310dY1HQuf9z//+U8uuuiio5Ih9uTPu7mf4c78/e62Y2JcLtdRH3TTKr/19fVdUKPUW7VqFY899hiLFi2KtT4A3H///Zx11lkYjUai0Sg33HAD3/3ud9m+fXu37Ctuj7auyeVyATT7M9BbPv+nnnqKWbNmxcbCNOmNn/eR2vP73dv/BkSjUa644gr69OnDNddcE9s+aNAgNm/ezPDhwwEtH92MGTMYMWIE999/f1dVNyntuabe/nmvWbOGzZs388QTT8Rt78mfd0s/w535+91tW2JsNhsejyduW9Nzu93eFVVKqS+//JILLriAm266ieuvvz5u39lnnx0bsW0ymXjooYeorq5mzZo1XVFVXbR1TU1NtM39DPSGz7+kpISPPvoo7he/SW/8vI/Unt/v3vw3QFEU/t//+398+eWXvP/++3EDOEeMGBH7QgMtlcsll1zS4sy2nqA919SbP2/Q/mkZMWIEs2bNitveUz/v1n6GO/P3u9sGMcOGDTuqSWnfvn2YTCYGDRrURbVKjXXr1nHWWWdx5ZVX8thjj7V5fNMPy5HTcHuyI6+pacrekT8DJSUlDB06tHMrlwJPP/00TqeTiy66qM1je+Pn3Z7f7976N0BVVRYuXMiHH37IsmXLjmqJa05WVlav+vzh6GvqrZ83aMsrvPzyyyxcuDBuDGRLuvvn3dbPcGf+fnfbIObcc89l2bJlVFdXx7a9/PLLnH766Tgcji6smb6++uorzjrrLC6//HJ+//vfH7Xf7Xbj8/nitr355psAHH/88Z1SR72155omTJhAcXFx3CDetWvXsnfvXs4777zOq2wKqKrK008/zWWXXXbUfxy98fNuTnt+v3vj3wBVVbnmmmv44IMPWLp0KaNHjz7qmLKysrjnoVCId999l2nTpnVWNXXXnmvqjZ93k5dffplgMMgVV1xx1L6e9nm352e4M3+/u20CyFAoxMknn0w4HGbhwoWsXbuW119/nZUrV3Lcccd1dfV0UVJSwtSpU+nbty/33ntv3L6zzjqLvLw8duzYwbx585g/fz5Dhgxhw4YN/OUvf+Gmm27it7/9bRfVPDn/v507ZF0dCuM4/usWkzDEIojBYp6ITQTfgE0FgwbBJCbDMMqyRdAgGDSaVoxWgwu+ARGMJhG8Qa5cuXrv4H+vf498P21sjPPw7Bx+Y4cFrWk+n6tUKqnVaikajcp1Xdm2rel0+o2j/zrP85TP5x9ubPuEfp/PZ81mM0nX/T3xeFzValXhcFiFQkFSsPlt2hoQpO52u61+v69ut3u3+FuWdfvM0Gg0dDgclMvldLlcNB6PtdvttFwulUgkXl/YXwSpO0hNn9jvnzKZjCzLul3/K9P6HeQZfuX8ftsQI11/DjYYDLRerxWJRFSr1ZRMJr97WP+M7/tyHOfhuV6vd/tOut/vNRqNtN1uZVmWisWiMT8+eyZoTavVSpPJRMfjUdlsVuVy+W7Ts4mGw6E2m41c13143vR+P3vjjMVid0EsyPw2aQ0IUrfjOPJ9/7dr0um0Op3O7XixWMjzPJ1OJ6VSKVUqFYVCof83+C8I2u8gNX1av6XrBtZ6va5msynbth/ey6R+B32GXzW/3zrEAAAAPPO2e2IAAAD+hBADAACMRIgBAABGIsQAAAAjEWIAAICRCDEAAMBIhBgAAGAkQgwAADASIQYAABiJEAMAAIxEiAEAAEYixAAAACP9AKRG34rvL2bQAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt\n", "import numpy as np\n", "\n", "i, j, k = 12, 1073, 7999\n", "\n", "timeline = np.linspace(0, 200)\n", "regression.plot(\"hf\", timeline, z1[i], z2[i], z3[i], label=f\"hf of the {i}-th asset\")\n", "regression.plot(\"hf\", timeline, z1[j], z2[j], z3[j], label=f\"hf of the {i}-th asset\")\n", "regression.plot(\"hf\", timeline, z1[k], z2[k], z3[k], label=f\"hf of the {i}-th asset\")\n", "plt.show()" ] } ], "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.0rc1" } }, "nbformat": 4, "nbformat_minor": 5 }