diff options
Diffstat (limited to 'PVCM/cama/fr/ma42 Surrelaxation pour Gauss-Seidel -- Exercice.ipynb')
| -rw-r--r-- | PVCM/cama/fr/ma42 Surrelaxation pour Gauss-Seidel -- Exercice.ipynb | 236 |
1 files changed, 189 insertions, 47 deletions
diff --git a/PVCM/cama/fr/ma42 Surrelaxation pour Gauss-Seidel -- Exercice.ipynb b/PVCM/cama/fr/ma42 Surrelaxation pour Gauss-Seidel -- Exercice.ipynb index bb2a484..5fad716 100644 --- a/PVCM/cama/fr/ma42 Surrelaxation pour Gauss-Seidel -- Exercice.ipynb +++ b/PVCM/cama/fr/ma42 Surrelaxation pour Gauss-Seidel -- Exercice.ipynb @@ -1,25 +1,4 @@ { - "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.8.10" - } - }, - "nbformat": 4, - "nbformat_minor": 2, "cells": [ { "cell_type": "markdown", @@ -32,7 +11,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 2, "metadata": {}, "outputs": [], "source": [ @@ -182,21 +161,58 @@ "cell_type": "code", "execution_count": null, "metadata": {}, - "outputs": [], - "source": [] + "outputs": [ + { + "data": { + "text/plain": [ + "[np.float64(7.0), np.float64(-2.0), np.float64(-1.5)]" + ] + }, + "execution_count": 7, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "def solve_triangular(L,b):\n", + " x = [b[0]/L[0,0]]\n", + " for i in range(1,len(L)):\n", + " t = L[i,:i]\n", + " x.append((b[i] - t @ np.array(x))/L[i,i])\n", + " return x" + ] }, { "cell_type": "code", - "execution_count": null, + "execution_count": 12, "metadata": {}, "outputs": [], - "source": [] + "source": [ + "def gauss_seidel_r(A, b, x0, w, n):\n", + " D = np.diag(np.diag(A))\n", + " L = np.tril(A,-1)\n", + " U = np.triu(A,1)\n", + " xi = [x0]\n", + " l = D/w + L\n", + " for k in range(n):\n", + " B = (((1 - w)/w) * D - U) @ xi[k] + b\n", + " xi.append(solve_triangular(l, B))\n", + " return xi" + ] }, { "cell_type": "code", - "execution_count": null, + "execution_count": 13, "metadata": {}, - "outputs": [], + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[np.float64(0.9999852128485041), np.float64(0.9999904902977096), np.float64(1.0000020873614475), np.float64(1.000040295309542)]\n" + ] + } + ], "source": [ "np.random.seed(123)\n", "A = np.random.randint(10, size=(4,4))\n", @@ -209,7 +225,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 15, "metadata": {}, "outputs": [], "source": [ @@ -228,9 +244,27 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 16, "metadata": {}, - "outputs": [], + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Itération du minimum : 100 100\n" + ] + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAABHwAAAF0CAYAAABR4y7lAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjMsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvZiW1igAAAAlwSFlzAAAPYQAAD2EBqD+naQAAe2ZJREFUeJzt3XdclfX7x/H3OYeNgAIKooC4B+49cGSpuDIbZuUqLW05mlbf1reyr5WZmVpW2s4sszTLtFScqbjFvcAF4mDKPPfvD5NfpJYocB/g9Xw8zh/c5+Y+17k/xufquj/DYhiGIQAAAAAAAJQaVrMDAAAAAAAAQOGi4AMAAAAAAFDKUPABAAAAAAAoZSj4AAAAAAAAlDIUfAAAAAAAAEoZCj4AAAAAAAClDAUfAAAAAACAUoaCDwAAAAAAQClDwQcA8I+OHTumSpUq6bHHHjM7FAAAAIdH7gRHQcEHZdbs2bNlsViu+Fq+fLnZIRabzp07Kzw8vFg/8/Dhw7JYLJo9e3axfm5pNnToUFWrVi3fsWrVqmno0KHXfM2cnBzdeeeduuGGG/Tmm29eX4AAgFKBHOrfXU+ec/z4cb344ovasmXLJe+9+OKLslgs1x8gJJE7ofRzMjsAwGyzZs1S3bp1Lzlev359E6IBCtf3338vb2/va/798ePHy8nJSZ9++ikJJgAgH3KoonH8+HG99NJLqlatmpo0aZLvveHDh6tHjx7mBFZGkDuhNKHggzIvPDxcLVq0KNDvGIahjIwMubu7X/Le+fPn5ebmVix/4HNzc5WTkyNXV9ci/ywU3Pnz5y/7b6Q4NW3a9Lp+/4033iikSAAApU1ZyqH+Ke7iVLVqVVWtWtXUGIoSuRNQuJjSBVwFi8Wihx9+WDNmzFC9evXk6uqqTz75JG9I86+//qp7771XFStWlIeHhzIzMyVJc+bMUdu2beXp6aly5cqpe/fu2rx5c75rd+7cWZ07d77kM/8+xPTi0OCJEyfqlVdeUVhYmFxdXbVs2bIrxv3ee++pY8eOqlSpkjw9PdWwYUNNnDhR2dnZlz1/5cqVatOmjdzd3VWlShX95z//UW5ubr5zpk+frsaNG6tcuXLy8vJS3bp19cwzz+Q7Z8eOHbr55ptVoUIFubm5qUmTJvrkk0/+6RZf9jtfdLnhy4ZhaNq0aWrSpInc3d1VoUIF3XbbbTp48OC/fs7F6+3cuVMDBw6Uj4+PAgICdO+99yopKSnfuRkZGRo/frzCwsLk4uKiKlWq6KGHHtK5c+fynVetWjX17t1b8+bNU9OmTeXm5qaXXnpJy5cvl8Vi0ZdffqmnnnpKlStXVrly5dSnTx/Fx8crJSVF999/v/z9/eXv769hw4YpNTU137UL2o5/j+uvw5LtdrteeeUV1alTR+7u7ipfvrwaNWqkd955J9/v7du3T3fddZcqVaokV1dX1atXT++9994l109OTtbjjz+e7/6MGTNGaWlp/xobAKD0K6k51JXilq6+j/y7/fv3a9iwYapVq5Y8PDxUpUoV9enTR9u3b887Z/ny5WrZsqUkadiwYXnT5F588UVJl+ZE/fr1U2hoqOx2+yWf17p1azVr1izvZ3IncieUPYzwQZl38QnPX1ksFtlstnzH5s+fr5UrV+r5559XYGCgKlWqpA0bNkiS7r33XvXq1UufffaZ0tLS5OzsrNdee03PPfechg0bpueee05ZWVl64403FBERofXr11/zcOcpU6aodu3aevPNN+Xt7a1atWpd8dwDBw7orrvuyutQtm7dqldffVW7d+/Wxx9/nO/ckydP6s4779TTTz+tl19+WT/99JNeeeUVnT17VlOnTpUkff3113rwwQf1yCOP6M0335TVatX+/fsVExOTd509e/aoXbt2qlSpkqZMmSI/Pz99/vnnGjp0qOLj4/Xkk09e0/f+uwceeECzZ8/Wo48+qv/97386c+aMXn75ZbVr105bt25VQEDAv17j1ltv1YABA3Tfffdp+/btGj9+vCTl3RvDMNSvXz/99ttvGj9+vCIiIrRt2za98MILWrt2rdauXZvvyeCmTZu0a9cuPffccwoLC5Onp2de5/3MM8+oS5cumj17tg4fPqzHH39cAwcOlJOTkxo3bqyvvvpKmzdv1jPPPCMvLy9NmTIl77oFacd/M3HiRL344ot67rnn1LFjR2VnZ2v37t35krCYmBi1a9dOISEheuuttxQYGKjFixfr0UcfVWJiol544QVJUnp6ujp16qSjR4/qmWeeUaNGjbRz5049//zz2r59u5YuXcpQZgAoxUpzDnWluK+2j7yc48ePy8/PT6+//roqVqyoM2fO6JNPPlHr1q21efNm1alTR82aNdOsWbPyvnuvXr0k6Yqjeu69917dfPPN+v3333XjjTfmHd+9e7fWr1+fL58gdyJ3QhlkAGXUrFmzDEmXfdlstnznSjJ8fHyMM2fOXPYagwcPznc8NjbWcHJyMh555JF8x1NSUozAwEDjjjvuyDvWqVMno1OnTpfEN2TIECM0NDTv50OHDhmSjBo1ahhZWVkF/r65ublGdna28emnnxo2my3fd+nUqZMhyfjhhx/y/c6IESMMq9VqHDlyxDAMw3j44YeN8uXL/+Pn3HnnnYarq6sRGxub73hkZKTh4eFhnDt3Lt/3mTVr1hW/80UvvPCC8dc/V2vXrjUkGW+99Va+8+Li4gx3d3fjySef/McYL15v4sSJ+Y4/+OCDhpubm2G32w3DMIxffvnlsufNmTPHkGR88MEHecdCQ0MNm81m7NmzJ9+5y5YtMyQZffr0yXd8zJgxhiTj0UcfzXe8X79+hq+v7xVj/6d2vNz9Cw0NNYYMGZL3c+/evY0mTZpc8fqGYRjdu3c3qlataiQlJeU7/vDDDxtubm55nzlhwgTDarUaGzZsyHfet99+a0gyFi1a9I+fAwAomcpCDnWluK+2j7xcnvN3OTk5RlZWllGrVi1j7Nixecc3bNhwxd/9e06UnZ1tBAQEGHfddVe+85588knDxcXFSExMNAyD3IncCWUVU7pQ5n366afasGFDvtcff/xxyXk33HCDKlSocNlr3Hrrrfl+Xrx4sXJycjR48GDl5OTkvdzc3NSpU6fr2r2ib9++cnZ2vqpzN2/erL59+8rPz082m03Ozs4aPHiwcnNztXfv3nznenl5qW/fvvmO3XXXXbLb7YqKipIktWrVSufOndPAgQP1ww8/KDEx8ZLP/P3339W1a1cFBwfnOz506FClp6dr7dq1Bfm6l7Vw4UJZLBbdc889+e5vYGCgGjdufNX39+/ft1GjRsrIyFBCQkLed7kY+1/dfvvt8vT01G+//XbJ79euXfuyn9W7d+98P9erV0+S8p7c/fX4mTNn8g1NLkg7/ptWrVpp69atevDBB7V48WIlJyfnez8jI0O//fabbrnlFnl4eOS7vz179lRGRobWrVsn6UI7hIeHq0mTJvnO6969O7u0AEAZUJpzqMvFXZA+8nJycnL02muvqX79+nJxcZGTk5NcXFy0b98+7dq165q+k5OTk+655x7Nmzcvb2pVbm6uPvvsM918883y8/OTRO5E7oSyiildKPPq1at3VQsOVq5c+arfi4+Pl6S8Odh/Z7Vee631n+L4q9jYWEVERKhOnTp65513VK1aNbm5uWn9+vV66KGHdP78+XznX24Yb2BgoCTp9OnTkqRBgwYpJydHM2fO1K233iq73a6WLVvqlVde0U033ZR37uViDAoKynet6xEfHy/DMK449Lh69epXdZ2LSdBFF4cYX7w3p0+flpOTkypWrJjvPIvFosDAwEu+yz+1ja+vb76fXVxc/vF4RkaGypUrV+B2/Dfjx4+Xp6enPv/8c82YMUM2m00dO3bU//73P7Vo0UKnT59WTk6O3n33Xb377ruXvcbFQl98fLz2799/xeT5cgVBAEDpUVpzqCudX5A+8nLGjRun9957T0899ZQ6deqkChUqyGq1avjw4QXuz//q3nvv1VtvvaWvv/5aDzzwgBYvXqwTJ05o2LBheeeQO5E7oWyi4ANcpX+aT/v39/z9/SVJ3377rUJDQ//xum5ubpcsdidd+Q/+1c7rnT9/vtLS0jRv3rx8MWzZsuWy519MsP7q5MmTkvJ37sOGDdOwYcOUlpamqKgovfDCC+rdu7f27t2r0NBQ+fn56cSJE5dc6/jx45L+/95cjpubW95ijX/193vh7+8vi8WilStXXnZ3jcLatczPz085OTk6depUvsTFMAydPHnykmS0KOZcF7Qd/42Tk5PGjRuncePG6dy5c1q6dKmeeeYZde/eXXFxcapQoYJsNpsGDRqkhx566LLXCAsLk3ShHdzd3a84F/6f2hoAUHaUtBzqSucXpI+8nM8//1yDBw/Wa6+9dkm85cuXL1Bsf1W/fn21atVKs2bN0gMPPKBZs2YpKChI3bp1yzuH3IncCWUTBR+gCHTv3l1OTk46cODAJUOV/65atWqaO3euMjMz8zrb06dPa82aNfL29r7mGC52oH/twA3D0MyZMy97fkpKin788cd8Q3W//PJLWa1WdezY8ZLzPT09FRkZqaysLPXr1087d+5UaGiounbtqu+//17Hjx/PG9UjXRj27eHhoTZt2lwx5mrVqikhIUHx8fF5T6CysrK0ePHifOf17t1br7/+uo4dO6Y77rjjKu7GtenatasmTpyozz//XGPHjs07/t133yktLU1du3Ytss++qKDtWBDly5fXbbfdpmPHjmnMmDE6fPiw6tevry5dumjz5s1q1KhR3lOzy+ndu7dee+01+fn5/WOCCwDA1XKEHOpKPDw8rrqPvByLxXJJYeWnn37SsWPHVLNmzbxjfx81czWGDRumUaNGadWqVVqwYIHGjRuXb/FscidyJ5RNFHxQ5u3YseOSHSYkqUaNGpcMR71a1apV08svv6xnn31WBw8eVI8ePVShQgXFx8dr/fr18vT01EsvvSTpwjSp999/X/fcc49GjBih06dPa+LEidedqNx0001ycXHRwIED9eSTTyojI0PTp0/X2bNnL3u+n5+fRo0apdjYWNWuXVuLFi3SzJkzNWrUKIWEhEiSRowYIXd3d7Vv316VK1fWyZMnNWHCBPn4+OQ9sXnhhRe0cOFCdenSRc8//7x8fX31xRdf6KefftLEiRPl4+NzxZgHDBig559/XnfeeaeeeOIJZWRkaMqUKZdsDd++fXvdf//9GjZsmDZu3KiOHTvK09NTJ06c0KpVq9SwYUONGjXquu7fxXvYvXt3PfXUU0pOTlb79u3zdppo2rSpBg0adN2fcTUxFKQd/02fPn0UHh6uFi1aqGLFijpy5IgmT56s0NDQvN1K3nnnHXXo0EEREREaNWqUqlWrppSUFO3fv18LFizIm58/ZswYfffdd+rYsaPGjh2rRo0ayW63KzY2Vr/++qsee+wxtW7dutDuBQDAsZTWHOqfXG0feTm9e/fW7NmzVbduXTVq1EjR0dF64403LtmBq0aNGnJ3d9cXX3yhevXqqVy5cgoKCsr3IO3vBg4cqHHjxmngwIHKzMy8ZA0dcidyJ5RRJi4YDZjqn3aYkGTMnDkz71xJxkMPPXTFa/x9pf2L5s+fb3Tp0sXw9vY2XF1djdDQUOO2224zli5dmu+8Tz75xKhXr57h5uZm1K9f35gzZ84Vd5h44403rvo7LliwwGjcuLHh5uZmVKlSxXjiiSeMn3/+2ZBkLFu2LO+8Tp06GQ0aNDCWL19utGjRwnB1dTUqV65sPPPMM0Z2dna+OLt06WIEBAQYLi4uRlBQkHHHHXcY27Zty/e527dvN/r06WP4+PgYLi4uRuPGjS/ZaeJKu1csWrTIaNKkieHu7m5Ur17dmDp16iU7Ulz08ccfG61btzY8PT0Nd3d3o0aNGsbgwYONjRs3/uN9uXi9U6dO5Tt+sT0PHTqUd+z8+fPGU089ZYSGhhrOzs5G5cqVjVGjRhlnz57N97uhoaFGr169LvmsiztNzJ0797Kf9fd/O5eL7Wrb8Wp2mnjrrbeMdu3aGf7+/oaLi4sREhJi3Hfffcbhw4fz/d6hQ4eMe++916hSpYrh7OxsVKxY0WjXrp3xyiuv5DsvNTXVeO6554w6deoYLi4uho+Pj9GwYUNj7NixxsmTJy+5HwCAkq8s5FBXivvi9f6tj7xcnnP27FnjvvvuMypVqmR4eHgYHTp0MFauXHnZ3ca++uoro27duoazs7MhyXjhhRcMw7h0l66/uuuuuwxJRvv27a/4vcidyJ1QtlgMwzCKrpwEAAAAAACA4sa27AAAAAAAAKUMBR8AAAAAAIBShoIPAAAAAABAKUPBBwAAAAAAoJRxyILPwoULVadOHdWqVUsffvih2eEAAAAAAACUKA63S1dOTo7q16+vZcuWydvbW82aNdMff/whX19fs0MDAAAAAAAoEZzMDuDv1q9frwYNGqhKlSqSpJ49e2rx4sUaOHDgVf2+3W7X8ePH5eXlJYvFUpShAgCA62AYhlJSUhQUFCSr1SEHHZcZ5E8AAJQMBcmfCr3gExUVpTfeeEPR0dE6ceKEvv/+e/Xr1y/fOdOmTdMbb7yhEydOqEGDBpo8ebIiIiIkScePH88r9khS1apVdezYsav+/OPHjys4OLhQvgsAACh6cXFxqlq1qtlhlHgLFy7UY489JrvdrqeeekrDhw+/6t8lfwIAoGS5mvyp0As+aWlpaty4sYYNG6Zbb731kvfnzJmjMWPGaNq0aWrfvr3ef/99RUZGKiYmRiEhIbrcDLOCPGny8vKSdOHLe3t7X/sXAQAARSo5OVnBwcF5fTeuXU5OjsaNG5dvSnz//v2veko8+RMAACVDQfKnQi/4REZGKjIy8orvT5o0Sffdd1/eU6fJkydr8eLFmj59uiZMmKAqVarkG9Fz9OhRtW7d+orXy8zMVGZmZt7PKSkpkiRvb28SFgAASgCmEF2/650Sf7ENyJ8AACgZriZ/KtYJ81lZWYqOjla3bt3yHe/WrZvWrFkjSWrVqpV27NihY8eOKSUlRYsWLVL37t2veM0JEybIx8cn78VwZAAAUNJERUWpT58+CgoKksVi0fz58y85Z9q0aQoLC5Obm5uaN2+ulStX5r13vVPiAQBA6VOsBZ/ExETl5uYqICAg3/GAgACdPHlSkuTk5KS33npLXbp0UdOmTfXEE0/Iz8/vitccP368kpKS8l5xcXFF+h0AAAAK28Up8VOnTr3s+xenxD/77LPavHmzIiIiFBkZqdjYWEm67inxAACg9DFll66/JyCGYeQ71rdvX/Xt2/eqruXq6ipXV9dCjQ8AAKA4mT0lPjk5uRC+BQAAcCTFOsLH399fNpstbzTPRQkJCZeM+gEAAABT4gEAwLUp1oKPi4uLmjdvriVLluQ7vmTJErVr1644QwEAACgRmBIPAACuRaFP6UpNTdX+/fvzfj506JC2bNkiX19fhYSEaNy4cRo0aJBatGihtm3b6oMPPlBsbKxGjhxZ2KEAAACUGkyJBwAABVHoBZ+NGzeqS5cueT+PGzdOkjRkyBDNnj1bAwYM0OnTp/Xyyy/rxIkTCg8P16JFixQaGlrYoQAAAJR4TIkHAADXotALPp07d77sThF/9eCDD+rBBx8s7I8GAAAodf46Jf6WW27JO75kyRLdfPPNJkYGAAAcmSm7dAEAAOD/MSUeAAAUNgo+BWS3G7JaLf9+IgAAwFViSjwAAChsFuPf5l+VMMnJyfLx8VFSUpK8vb0L7bq/747X1N/364a6lfTwDbUK7boAAJRVRdVno+CKqi1mrz6kbzYe1cjONdS3cVChXRcAgLKqIH12sW7LXpKdTcvWpthz+m7TsX9dowgAAADS8aQMxZxI1qp9p8wOBQCAMoeCz1XqER4od2ebDiWmaXPcObPDAQAAcHjta/pLklbtS+SBGQAAxYyCz1XydHVSZHigJGnepqMmRwMAAOD4WlarIBebVceTMnT4dLrZ4QAAUKZQ8CmA/s2qSpIWbD2hzJxck6MBAABwbB4uTmoWWl6StGp/ornBAABQxlDwKYC2NfxU2cdNSeez9fuuBLPDAQAAcHgd/pzWtXofBR8AAIoTBZ8CsFkt6te0iiTpu03HTI4GAADA8V1cx2fNgUTl2lnHBwCA4kLBp4D6/1nwWb4nQadTM02OBgAAwLE1rOIjLzcnJWfkaMexJLPDAQCgzHAyO4CSplaAlxpV9dG2o0n6cetxDWsfdtW/a7cb2nE8SWsPnFZmjl0WSRaLZLFYVMHDRU2Cy6tOoJdsVkvRfQEAAIBi5GSzqm11P/0aE69V+xPVOLi82SEBAFAmUPC5Brc2q6ptR5M0b9Oxfy34ZObkas2B01oaE6+lu+IVn/zPo4I8XWxqVLW8modWULcGAWpYxUcWCwUgAABQcnWo5a9fY+K1en+iHupS0+xwAAAoEyj4XIM+jYP034Ux2n4sSXvjU1Q7wOuScwzD0OKdJ/XyghgdT8rIO+7hYlOHmv7yK+cqyZBhSIYhHTt3Xlvizik1M0drD57W2oOnNXXZftWqVE79m1XVLU2rKNDHrRi/JQAAQOG4uI7PxsNndT4rV+4uNpMjAgCg9KPgcw18PV3UpW4lLYmJ13ebjmp8ZL187x85naYXftyp5XtOSZIqermqe4MA3VgvQG2q+8nN+fJJTq7d0L6EFG2OPadV+xO1NCZe+xJS9b9fdmvi4t3qVLuiHuhYQ22q+zLqBwAAlBjV/T1V2cdNJ5IytPHIGUXUqmh2SAAAlHoUfK7Rrc2qaElMvL7deFQWWeThYpOHi02nUjM1a/VhZeXY5WKz6oFO1fVg55pX9STLZrWobqC36gZ6a2CrECVnZOvn7Sf0XfQxrT98Rsv3nNLyPafUJLi8RnaqoW71A2RlvR8AAODgLBaL2tf017fRR7VqfyIFHwAAigEFn2vUpW4lVfBw1um0LM1YceCS9zvU9NfLNzdQ9YrlrvkzvN2cNaBliAa0DNHhxDR9vPqQ5myI05a4cxr5ebRqVPTUo11rqU+jIAo/AADAoXX4s+Czen+i2aEAAFAmWAzDMMwOojAlJyfLx8dHSUlJ8vb2LtLP2hx7Vr/tSlB6Vq7OZ+coPStXObmGejasrJ4NA4tk2lViaqZmrz6sT9YeVkpGjiSpQZC3nuxRVx1r+TPVCwBQYhRnn41/VhxtkZCSoVav/iaLRYp+7ib5eroUyecAAFCaFaTPZoTPdWgaUkFNQyoU62f6l3PV493r6IFO1TV79WG9H3VQO48na8jH69W2up+ejqzLdqcAAMDhVPJyU50AL+2JT9HaA6fVq1Fls0MCAKBUo+BTQnm5OeuRrrV0d5tQvbdsvz5be0RrD57Wze+tVv9mVfRUj7oK8GZXLwAA4Dja1/TXnvgUrdqfeEnBJzUzR/sTUrUvPkX7E1J1KiVTKZk5SvvzlZFtl5uzVe4uNnm4OMndxSY/TxdVKe+uKhXcFVTeXcEVPORfzoURzwAAiCldpcbRs+matGSv5m06JunC9u8Pdq6h4RHVr7grGAAAZiqrfbYjKq62+H13vO6dvVEuNqsqernKYpGsFosyc3IVn5xZKJ/h6+miepW9VC/QW3Ure6txVR/VrFSOIhAAoFQoSJ9NwaeU2Rp3Ti8t2KlNseckSVXKu+s/veupe4OiWVMIAIBrVdb7bEdSXG2RlpmjjhOX6XRa1mXfr+jlqtoB5VSrkpcq+7ipnJuTyrleeLk525SRnXth7cSsXKVl5ehUSqaOnzuvY+fO69jZ8zqZnCH7ZTLbCh7OalHNV62q+aplmK8aVvGRjQ0vAAAlEAWfMp48GoahH7ce1+s/79aJpAxJUkQtf73Yt4FqXMeuYQAAFCb6bMdRnG2RlJ6tI2fSZBiS3TBkNySb1aJqfh4q73F9CzlnZOdqb3yKdp9IUcyJZMWcSNa2o+eUkW3Pd14FD2d1rF1RnetUVMdaFeVXzvW6PhcAgOJCwYfkUZKUnpWj6csP6P0VB5WVa5ezzaL7OlTXIzfUlKcryzcBAMxFn+04SnNbZOXYteN4kjYcOqMNh8/oj0Nn8nY6lSSLRWoaXF49G1ZWZMPKqlLe3cRoAQD4ZxR8SmnCcq0OJ6bppQU7tWzPKUlSoLeb/tO7fpFtHQ8AwNWgz3YcZaktcnLt2hR7Tsv3JGj5nlOKOZGc7/3GweXVq2Gg+jauokAfNsAAADgWCj5lJGEpqKUx8Xpp4U7FnTkviWleAABz0Wc7jrLcFieTMrR450kt2n5C6w+f0cXM2GKR2tfwV/9mVdS9QSCjowEADoGCTxlNWK5GRnaupi8/oOkrDigr58I0rxER1fXwDTXl4UIiAwAoPvTZjoO2uCAhJUOLd8ZrwdbjWn/oTN5xDxebejasrIGtQtQspDwjpAEApqHgQ8Lyr46cTtOLP/7/NK8gHzc917u+IsOZ5gUAKB702Y6DtrhU3Jl0fb/5mOZtOqrDp9PzjtcJ8NLAVsG6pWlV+Xg4mxghAKAsouBDwnJVDMPQrzHxenlBjI6duzDNq0NNf73Yt75qVvIyOToAQGlHn+04aIsrMwxD0UfO6usNcVq47Xjejl9uzlb1a1JFQ9pVU73K3DMAQPGg4EPCUiDns3I1fcUBzfhzmpeT1aJ7O4TpkRtqysuNJ1cAgKJBn+04aIurk3Q+Wz9sOaYv/4jV7pMpecdbhflqaLtq6lY/QE42q4kRAgBKOwo+JCzX5MjpNP13YYyW7kqQJPmXc9WTPerotmZVZbUyzQsAULjosx0HbVEwhmFo45Gzmr36sH7ZeVK59gvpdJXy7rq3Q5gGtAxWORZ5BgAUAQo+JCzXZdnuBL28MEaHEtMkSY2q+uiFPg3UPLSCyZEBAEoT+mzHQVtcuxNJ5/X5uiP6an2czqRlSZK83Jx0d+tQDW1Xja3dAQCFioIPCct1y8qxa/aaQ5ry236lZuZIkvo2DtKTPeqoagUPk6MDAJQG9NmOg7a4fhnZufpu01F9tPKQDv750MzZZtEtTavogU41VKNiOZMjBACUBhR8SFgKTUJKht74ZY++3XRUhiG5OFk1vEOYRnWuwfo+AIDrQp/tOGiLwmO3G/p9d4I+WHkwb2t3i0WKDA/UqE411bCqj8kRAgBKMgo+JCyFbsexJL3yU4zWHbyQuPiXc9GYG2trQMtgObM4IQDgGtBnOw7aomhEHzmr6cv3562PKEkda1fU6K411TzU18TIAAAlFQUfEpYiYRiGlu5K0IRFu/KGKof5e+qJ7nUUGR4oi4WFnQEAV48+23HQFkVr98lkzVh+QAu2nchb4Ll9TT89ekMtta7uZ3J0AICShIIPCUuRys6168s/YjXlt306/efihI2Dy+upHnXUroa/ydEBAEoK+mzHQVsUj9jT6Zq2fL++jT6qnD8LP63DfDX2ptpqQ+EHAHAVKPiQsBSL1MwczYw6qJkrDyo9K1eSFFHLX491q6MmweXNDQ4A4PDosx0HbVG8jp5N1/TlBzR341Fl5dolXRjxM/bG2mpRjaleAIAro+BDwlKsTqVk6t3f9+mr9bHKzr3wz+nGegF6rFtt1atMGwAALo8+23HQFuY4kXRe7y3brzkb4vJyqI61K2rcTbV5eAYAuCwKPiQspog7k653ftuneZuO6s9RyurVqLJGd62l2gFe5gYHAHA49NmOg7Yw19Gz6Xpv2X7N3fj/U7261Q/QY93qqE4gORQA4P9R8CFhMdWBU6l6e8leLdx2QtKFrUh7NaysRyn8AAD+gj7bcdAWjiH29IWHZ99vvvDwzGKR+jWpojE31lKon6fZ4QEAHAAFHxIWh7DrRLKm/LZPP+84KelC0tKzYWU9ckNN1Q2kbQCgrKPPdhy0hWPZn5CiSUv2atH2CzmUk9Wiu1qH6JEbaqmil6vJ0QEAzETBh4TFofy98CNdGKb8yA211LCqj4mRAQDMRJ/tOGgLx7T9aJLe+HWPovaekiR5uNg0PKK67u9YXeVcnUyODgBgBgo+JCwOadeJZE39fb8W7Tihi//qutSpqIdvqKXmoRXMDQ4AUOzosx0HbeHY1hxI1P9+3q2tR5MkSX6eLnq0ay3d1TpEzjarydEBAIoTBR8SFoe2PyFF7y07oB+2HMtb3LltdT89fENNtavhJ4vFYm6AAIBiQZ/tOGgLx2cYhn7ecVJvLt6jg4lpkqQwf0891aOuujcIIH8CgDKCgg8JS4lwODFN05cf0LzNR/O2Im0SXF4Pd6mprvUqkbgAQClHn+04aIuSIzvXrq83xOmdpXuVmJolSWoRWkHP9KqnZiGMmAaA0o6CDwlLiXL83Hl9EHVQX62PVWaOXZJUN9BLD3WpqZ4NK8tmpfADAKURfbbjoC1KntTMHL2/4oBmrjyojOwL+VPfxkF6KrKuqpR3Nzk6AEBRoeBDwlIinUrJ1EerDumztYeVlpUr6cJQ5VGdaqhf0ypycWKOOgCUJvTZjoO2KLlOJmXozV/36LtNR2UYkquTVSMiqmtU5xryZGFnACh1KPiQsJRoSenZmr3msGatOaRz6dmSpCrl3TWyU3Xd3iJYbs42kyMEABQG+mzHQVuUfDuOJem/C2P0x6EzkqSKXq56onsd3dasqqyMlgaAUoOCDwlLqZCamaMv/ziimSsP6VRKpqQLycsDHavrrtYh8nDhqRUAlGT02Y6DtigdDMPQrzHxem3RLh05nS5JaljFRy/2ra/mob4mRwcAKAwluuATFxenQYMGKSEhQU5OTvrPf/6j22+//ap/n4Sl9MnIztU3G+M0Y/kBHU/KkCT5erpoeESYBretpnIMVwaAEok+23HQFqVLVo5ds9cc0pTf9is1M0eS1K/JhfV9Kvuwvg8AlGQluuBz4sQJxcfHq0mTJkpISFCzZs20Z88eeXp6XtXvk7CUXlk5dn2/+aimLT+Q99SqvIezhncI0+B21eTt5mxyhACAgqDPdhy0Rel0KiVTby7eo2+i42QYkruzTY90ran7OoTJ1Ykp8gBQEhWkz3a4VXArV66sJk2aSJIqVaokX19fnTlzxtyg4BBcnKwa0DJEv43rpEl3NFZ1f0+dS8/Wm7/uVYfXf9eU3/YpJSPb7DABACh2cXFx6ty5s+rXr69GjRpp7ty5ZocEB1DRy1X/u62Rfnyog1qEVtD57FxN/GWPekxeqeV7EswODwBQxApc8ImKilKfPn0UFBQki8Wi+fPnX3LOtGnTFBYWJjc3NzVv3lwrV668puA2btwou92u4ODga/p9lE5ONqv6N6uqJeM66Z07m6hGRU8lZ+Ro0pK9ipi4TO8t26+0P4cvAwBQFjg5OWny5MmKiYnR0qVLNXbsWKWlpZkdFhxEw6o+mjuyrd4e0FgVvVx1KDFNQ2dt0IhPNyruTLrZ4QEAikiBCz5paWlq3Lixpk6detn358yZozFjxujZZ5/V5s2bFRERocjISMXGxuad07x5c4WHh1/yOn78eN45p0+f1uDBg/XBBx9cw9dCWWCzWnRzkyr6dWwnTRnYVDUqXhjx88biPYqYuEzvrzig839u7w4AQGnGCGn8G4vFoluaVtXvj3XSiIgwOVktWhITrxsnrdDU3/cpM4ecCQBKm+taw8disej7779Xv3798o61bt1azZo10/Tp0/OO1atXT/369dOECROu6rqZmZm66aabNGLECA0aNOhfz83MzMz7OTk5WcHBwcxBL4Ny7YYWbD2ud37bp0OJF55qVvJy1SNda2lAi2C5ODncDEYAKNPK0roxUVFReuONNxQdHa0TJ05ckj9JF0ZIv/HGGzpx4oQaNGigyZMnKyIiosCftXHjRg0dOlQ7duy46t8pS22BC/bFp+j5H3Zq7cHTkqTq/p566eYGiqhV0eTIAAD/xLQ1fLKyshQdHa1u3brlO96tWzetWbPmqq5hGIaGDh2qG2644V+LPZI0YcIE+fj45L2Y/lV22awW9WtaRUvGdtSbtzdW1QruSkjJ1H/m79CNk1bo+81HlWt3qDXKAQBlBCOk4WhqBXjpyxGt9c6dTeRfzlUHE9M06KP1evjLTUpIzjA7PABAISjUET7Hjx9XlSpVtHr1arVr1y7vvNdee02ffPKJ9uzZ86/XXLVqlTp27KhGjRrlHfvss8/UsGHDy57PCB9cSVaOXXM2xGrK7/t1KuXCv5E6AV56OrKuOtepKIvFYnKEAFC2ldVRJYyQhqNJzsjWpF/36tO1h2U3JC9XJz3Zo47uah0qm5V8CQAciem7dP39f6QNw7jq/7nu0KGD7Ha7tmzZkve6UrFHklxdXeXt7Z3vBUgXdvUa1LaaVjzRWU/1qCtvNyftiU/RsNkbdNfMP7Tt6DmzQwQAgBHSMJ23m7Ne7NtAPz7cQY2q+iglM0f/+WGn+k9fo53Hk8wODwBwjQq14OPv7y+bzaaTJ0/mO56QkKCAgIDC/Cjgqnm4OGlU5xpa+eQNeqBjdbk4WbX24Gn1nbpaD3+5id0pAACmSkxMVG5u7iW5UkBAwCU51ZWsXr1ac+bM0fz589WkSRM1adJE27dvv+L548ePV1JSUt4rLi7uur4DSofwKj76/sH2eqlvA5VzddLWuHPqO3W1Xlu0i40wAKAEKtSCj4uLi5o3b64lS5bkO75kyZJ8U7wAM/h4OGt8z3r6/bFO6t+0iiwWaeG2E+o6aYX+98tupWRkmx0iAKAMY4Q0HIHNatGQdtW0dFwn9WwYqFy7oQ+iDqrb5BVaue+U2eEBAAqgwAWf1NTUvERCkg4dOqQtW7bkLSo4btw4ffjhh/r444+1a9cujR07VrGxsRo5cmShBg5cq6oVPDRpQBMtfKSD2lb3U1aOXdOXH1CXN5frq/WxLOwMAChWjJCGIwr0cdO0u5vroyEtVNnHTXFnzmvQR+s17pstOpOWZXZ4AICrUOCCz8aNG9W0aVM1bdpU0oUCT9OmTfX8889LkgYMGKDJkyfr5ZdfVpMmTRQVFaVFixYpNDS0cCMHrlODIB99OaK1PhjUXNX8PJSYmqXx87arz7urtPHwGbPDAwCUEYyQhiPrWi9AS8Z10tB21WSxSPM2HdONk1box63HdR17vwAAisF17dLliMrqjh+4Plk5dn269rDe+W2fUjJyJEm3NK2i8ZF1VcnbzeToAKB0Kkt9dmpqqvbv3y9Jatq0qSZNmqQuXbrI19dXISEhmjNnjgYNGqQZM2aobdu2+uCDDzRz5kzt3LmzWB6alaW2wLXbFHtW47/brj3xKZKkG+sF6NVbwhVArgQAxaYgfTYFH+AvTqdm6o3FezRnY5wMQ/J0sWn0jbU0rH2YnG1FsqkdAJRZZanPXr58ubp06XLJ8SFDhmj27NmSpGnTpmnixIk6ceKEwsPD9fbbb6tjx47FEl9Zagtcn6wcu6Yt36/3lu1Xdq4hLzcnPdernu5oEXzVa04BAK4dBR8SFlynrXHn9MKPO7Ul7pwkqU6Al169JVwtqvmaGxgAlCL02Y6DtkBB7TmZoie/3aqtRy9s2x5Ry1+v39pIVcq7mxwZAJRuBemzGbIAXEbj4PKaN6qd3ritkXw9XbQnPkW3zVirp77dprMsVAgAAMq4OoFemvdgez3bs55cnaxauS9R3d+O0lfrY1nbBwAcBAUf4AqsVotubxGs38Z10sBWwZKkORvjdMNbyzV3YxzJDAAAKNNsVotGdKyun0dHqHloBaVm5mj8vO0a/PF6HTt33uzwAKDMo+AD/IsKni6a0L+RvhvVVnUDvXQ2PVtPfLtNgz5ar9jT6WaHBwAAYKrqFcvpmwfa6rle+Uf7fLOBB2QAYCYKPsBVah7qqwWPdNDTkXXl6mTVqv2J6jZ5hWZGHVROrt3s8AAAAExjs1o0PCL/aJ8nv9um4Z9sVEJyhtnhAUCZRMEHKABnm1UjO9XQL2M6qk11X2Vk2/Xqol3qP32N9pxMMTs8AAAAU10c7TM+sq5cbFb9tjtB3SZHacHW42aHBgBlDgUf4BqE+XvqqxFt9L9bG8rLzUnbjiap97sr9d6y/Yz2AQAAZZrNatEDnWpo4aMdFF7FW+fSs/XIV5v1yFebdS6dzS8AoLhQ8AGukcVi0YCWIfptXCfdWK+SsnMNvbF4j/pPX6O98Yz2AQAAZVvtAC99/2B7je5aSzarRQu2HlePySu1ct8ps0MDgDKBgg9wnSp5u2nm4BaadEdjeV8c7TNllWasOKBcOwsVAgCAssvZZtXYm2pr3qh2qu7vqZPJGRr00Xq9+ONOnc/KNTs8ACjVKPgAhcBisah/s6paMq6TutatpKxcu17/ebcGfrBOcWfYyQsAAJRtjYPL66dHIzSkbagkafaaw+r97krtOJZkcmQAUHpR8AEKUYC3mz4c0kITb2skTxeb1h8+o8h3Vurb6KNsSwoAAMo0dxebXro5XJ/c20qVvFx14FSabpm2WjNWHJCdUdEAUOgo+ACFzGKx6I4Wwfp5dEe1+HNb0sfnbtWozzexUCEAACjzOtWuqMVjOqp7gwBl5xp6/efduuvDdTp+7rzZoQFAqULBBygiIX4emvNAWz3Zo46cbRb9svOkIt9ZqXUHT5sdGgAAgKkqeLpoxj3N9b9bG8rd2aZ1B8+ox+Qo/bTthNmhAUCpQcEHKEI2q0UPdq6p7x9sr+r+njqRlKG7Zq7TpF/3sH07AAAo0y7ueLpodIQaV/VRckaOHvpyk576dpvSs3LMDg8ASjwKPkAxCK/iowWPdNDtzavKbkhTft+vAR+s09GzLOgMAADKtjB/T307qp0e6lJDFos0Z2Ocek9ZxYLOAHCdKPgAxcTT1Ulv3N5YUwY2lZerk6KPnFWvKau0NCbe7NAAAABM5Wyz6onudfXF8NYK9HbTwcQLCzp/uPIgCzoDwDWi4AMUs76Ngy4MXQ4ur6Tz2Rr+6Ua9tmiXspniBQAAyrh2Nfz18+gIdat/YUHnV37apXs/2aDTqZlmhwYAJQ4FH8AEwb4emvtAW93bPkyS9EHUQQ14fy27UwAAgDKvgqeL3h/UXK/0C5erk1XL95xS5DsrtWZ/otmhAUCJQsEHMImLk1XP96mvGfc0l5ebkzbFnlPPKSsVtfeU2aEBAACYymKx6J42ofrh4faqWamcElIydfdHf+jNxWx8AQBXi4IPYLIe4YH66ZEINazio3Pp2Roya72m/r6P+eoAAKDMqxvorQUPd9DAVsEyDGnqsv0aOHMdo6IB4CpQ8AEcQIifh+aObKs7W15IZt78da/u/yxaSeezzQ4NAADAVO4uNk3o30hT77qw8cWGw2fVc8pK/baLjS8A4J9Q8AEchJuzTa/f2kiv928oFyerlu6K181TV2nPyRSzQwMAADBd70ZBWvhoh7xR0fd9slGvLIxRVg5TvADgcij4AA7mzlYh+nZkW1Up767Dp9N1y7TV+mXHCbPDAgAAMF2on6e+HfX/G198uOqQbn9/reLOpJscGQA4Hgo+gANqVLW8FjzSQe1r+ik9K1cjP9+kt5fsZV0fAABQ5rk62fR8n/r6YFBz+bg7a2vcOfWaslK/7jxpdmgA4FAo+AAOytfTRZ8Ma5X3BOud3/Zp5OfRSs3MMTkyAAAA83VrEKifHu2gJsHllZyRo/s/i9Z/meIFAHko+AAOzMl2Yev2N25rJBebVb/GxKv/tNUMWwYAAJBUtYKHvnmgrUZEXHhA9tGfU7yOniVXAgAKPkAJcHuLYM15oI0qeblqb3yqbn5vtTYcPmN2WAAAAKZzcbLq2V71NXNwi79M8VrFLl4AyjwKPkAJ0TSkgn58uIPCq3jrTFqW7p75h76LPmp2WAAAAA7hpvoBeVO8ks5f2MVrwqJdys5liheAsomCD1CCBPq46ZsH2qpHg0Bl5dr12Nyt+t8vu1nMGQAAQP8/xWtY+2qSpPejDmrgB+t0Ium8uYEBgAko+AAljIeLk6bd3UwPdakhSZq+/IAe+nKTMrJzTY4MAADAfC5OVr3Qp4Gm391MXq5O2njkrHpNWaWovafMDg0AihUFH6AEsloteqJ7XU26o7GcbRb9vOOkBs5cp9OpmWaHBgAA4BAiG1bWgkc6qH7lC9Phh8xar0lL9iqXkdEAyggKPkAJ1r9ZVX12X2v5uDtrc+w53TJtjQ6eSjU7LAAAAIdQzd9T8x5sp4GtQmQY0pTf9mnIx+uVyEMyAGUABR+ghGtT3U/fjWqnYF93xZ5JV//pa7T+EDt4AQAASJKbs00T+jfU2wMay93ZplX7E9Vrykp2PAVQ6lHwAUqBmpXK6fsH26txcHmdS8/WPR/+oUXbT5gdFgAAgMO4pWlV/fBwe9WsVE7xyZm684N1en/FARkGU7wAlE4UfIBSwr+cq74e0Ubd6gcoK9euh77cpE/WHDY7LAAAAIdRO8BLPzzUXv2aBCnXbmjCz7t1/2fRSjqfbXZoAFDoKPgApYi7i03T72mue9pcmKf+wo87NfGX3Ty5AgAA+JOnq5PeHtBEr/QLl4vNqiUx8er97kptP5pkdmgAUKgo+ACljM1q0X9vDtfj3WpLkqYtP6DH525Tdq7d5MgAAAAcg8Vi0T1tQvPWQYw7c163Tl+jz9cd4UEZgFKDgg9QClksFj18Qy1NvLWRbFaLvtt0VCM+3ajzWblmhwYAAOAwGlb10cKHI3TTn1Pin5u/Q2PmbFFaZo7ZoQHAdaPgA5Rid7QM1szBzeXmbNXyPad0z0d/KCmdOeoAAAAX+Xg464NBzfVsz3qyWS36Yctx9Z26SntOppgdGgBcFwo+QCl3Q90AfX5fa3m7OSn6yFnd8f5axSdnmB0WAACAw7BYLBrRsbrm3N9Ggd5uOnAqTTe/t0rfRh81OzQAuGYUfIAyoEU1X30zsq0qeblqT3yKbp2+RocT08wOCwAAwKG0qOarnx7toIha/srItuvxuVv15LdbmRYPoESi4AOUEXUDvfXdqHYK9fPQ0bPndduMNYo5nmx2WAAAAA7Fr5yrZg9rpXE31ZbVIn2z8aj6vbdaB06lmh0aABQIBR+gDAn29dC3I9upXmVvJaZm6c4P1mpT7FmzwwIAAHAoNqtFj3atpc/vay3/chdGSPd5d5V+2HLM7NAA4KpR8AHKmIpervr6/jZqHlpByRk5uufDP7Rmf6LZYQEAADicdjX9tWh0B7Wp7qv0rFyN/nqLxs/bpoxspngBcHwUfIAyyMfdWZ/d10oRtfyVnpWrobM3aElMvNlhAQAAOJxKXm76YngbPXpDTVks0lfr43Tz1NXan8AULwCOjYIPUEZ5uDjpwyEt1L1BgLJy7Br5eTTDlAEAAC7DZrVoXLc6+uze/FO8vmMXLwAOjIIPUIa5Otn03l3N1L9ZFeXaDY2Zs0VzNsSaHRYAAIBD6lDrwhSvdjX8dD47V4/N3arHvtmqtMwcs0MDgEs4bMEnPT1doaGhevzxx80OBSjVnGxWvXlbY93TJkSGIT313XbNXn3I7LAAAAAcUiUvN312X2uNvfHCLl7fbTqqPlNXaefxJLNDA4B8HLbg8+qrr6p169ZmhwGUCVarRf+9OVwjIsIkSS8uiNG05ftNjgoAUFA8MAOKh81q0egba+mrEW0U6O2mg6fSdMu0Nfp07WEZhmF2eAAgyUELPvv27dPu3bvVs2dPs0MBygyLxaJnetbT6K61JEkTf9mjt37dQ9ICACUID8yA4tW6up9+Hh2hG+tVUlaOXc//sFP3fxats2lZZocGAAUv+ERFRalPnz4KCgqSxWLR/PnzLzln2rRpCgsLk5ubm5o3b66VK1cW6DMef/xxTZgwoaChAbhOFotFY2+qracj60qS3v19v15btIuiDwCUADwwA8xRwdNFMwe30At96svFZtWSmHj1eCdKa/Ynmh0agDKuwAWftLQ0NW7cWFOnTr3s+3PmzNGYMWP07LPPavPmzYqIiFBkZKRiY/9/IdjmzZsrPDz8ktfx48f1ww8/qHbt2qpdu/a1fysA12Vkpxp6qW8DSdLMlYf0wo87ZbdT9AGAa8UDM6B0s1gsGtY+TPMebKfqFT0Vn5ypuz/6Q//7Zbeyc+1mhwegjHIq6C9ERkYqMjLyiu9PmjRJ9913n4YPHy5Jmjx5shYvXqzp06fnJSHR0dFX/P1169bp66+/1ty5c5Wamqrs7Gx5e3vr+eefv+z5mZmZyszMzPs5OTm5oF8JwGUMaVdNrk5Wjf9+uz5de0SZ2Xa91r+hbFaL2aEBQIlz8YHZsGHDdOutt17y/sUHZtOmTVP79u31/vvvKzIyUjExMQoJCZF04YHZX3Oei3799Vdt2LAh74HZmjVrivz7ALi88Co+WvhIB/13YYy+Wh+n6csPaM3+RE2+s6nC/D3NDg9AGWMxrmOuhsVi0ffff69+/fpJkrKysuTh4aG5c+fqlltuyTtv9OjR2rJli1asWFGg68+ePVs7duzQm2++ecVzXnzxRb300kuXHE9KSpK3t3eBPg/Apb7ffFSPfbNVdkPq1yRIb97eWE42h1z+C0AJk5ycLB8fnzLXZ/89f5Kk1q1bq1mzZpo+fXresXr16qlfv35XNWpn/Pjx+vzzz2Wz2fIemD322GMFemAWHBxc5toCKEo/bz+hp+dtV9L5bLk72/R8n/q6s2WwLBYengG4dgXJnwr1/9oSExOVm5urgICAfMcDAgJ08uTJwvyoPOPHj1dSUlLeKy4urkg+ByirbmlaVe8ObCYnq0XztxzXI19tVlYOQ5MBoLBkZWUpOjpa3bp1y3e8W7duVz1aZ8KECYqLi9Phw4f15ptvasSIEVcs9lw838fHJ+8VHBx8Xd8BwKUiG1bWz6Mj1La6n85n52r8vO0a8Wm0TqdeOlIPAIpCkTym/3vV2jCMa6pkDx069B9H90iSq6urvL29870AFK5ejSpr+j3N5WKz6ucdJzXq82hlZOeaHRYAlAo8MANKr6Dy7vpieGs927OeXGxWLd0Vr+6TV+r33fFmhwagDCjUgo+/v79sNtslyUlCQsIlSQyAkuWm+gGaOaSFXJ2s+m13gkZ8ulHnsyj6AEBh4YEZUDpZrRaN6Fhd8x9qr9oB5ZSYmql7Z2/U099tU2pmjtnhASjFCrXg4+LioubNm2vJkiX5ji9ZskTt2rUrzI8CYIJOtStq9rBW8nCxaeW+RA2ZtZ5EBQCuEw/MgLKhfpC3fny4g4Z3CJPFIn29IU49Jkdp3cHTZocGoJQqcMEnNTVVW7Zs0ZYtWyRJhw4d0pYtW/K2XR83bpw+/PBDffzxx9q1a5fGjh2r2NhYjRw5slADB2COtjX89Nl9reTl6qT1h87ong//UFJ6ttlhAUCJxQMzoOxwc7bpud719dWINqpawV1Hz57XwJnr9N+FMUyXB1DoClzw2bhxo5o2baqmTZtKulDgadq0ad7CgAMGDNDkyZP18ssvq0mTJoqKitKiRYsUGhpauJEDME3zUF99OaKNyns4a0vcOd05c50SWYAQAK6IB2YA/qpNdT/9Mqaj7mwZLMOQPlp1SD3fWamNh8+YHRqAUuS6tmV3RGV1i1fADHtOpujuD/9QYmqmqlf01BfDW6uyj7vZYQEoIcpSn718+XJ16dLlkuNDhgzR7NmzJUnTpk3TxIkTdeLECYWHh+vtt99Wx44diyW+stQWgKNZtjtBT8/bpvjkTFks0r3tw/R4tzpyd7GZHRoAB1SQPpuCD4DrcigxTfd8+IeOnTuvKuXd9eWI1gr18zQ7LAAlAH2246AtAHMlnc/WfxfG6Nvoo5Kkan4emnhbY7UK8zU5MgCOpiB9dpFsyw6g7Ajz99Q3I9sqzN9Tx86d1+0z1mrPyRSzwwIAACgxfNyd9ebtjTVrWEsFervp8Ol03fH+Wv1n/g6lZLBWIoBrQ8EHwHWrUt5dcx5oozoBXkpIydQd769V9BHmoAMAABRElzqV9Ou4C2v7SNJn646o29tR+n13vMmRASiJKPgAKBSVvNw054E2ahZSXknns3X3h39o2e4Es8MCAAAoUbzdnPX6rY305fDWCvH10ImkDN07e6Me/Wozm2QAKBAKPgAKTXkPF30+vLU616mojGy7Rny6Ud9vPmp2WAAAACVOu5r+Wjymo+7vWF1Wi/Tj1uPq+tYKzdkQq1K2DCuAIkLBB0Ch8nBx0szBLXRL0yrKsRsaO2erPlx50OywAAAAShx3F5ue6VlP8x9qrwZB3ko6n62nvtuuAR+s0/6EVLPDA+DgKPgAKHTONqveur2x7m0fJkl65addevHHncq18zQKAACgoBpVLa8fHmqv53rVk7uzTesPnVHkO1Ga9OseZWTnmh0eAAdFwQdAkbBaLfpP73oaH1lXkjR7zWE98NlGpWflmBwZAABAyeNks2p4RHUtGddRN9StpOxcQ1N+369ub0dp2R7WTQRwKQo+AIqMxWLRA51qaNrdzeTqZNXSXQka8P46JSRnmB0aAABAiVS1goc+GtJC0+9upkBvN8WeSdewWRs06vNonUg6b3Z4ABwIBR8ARa5nw8r6ckQb+Xq6aPuxJPV7b7VijiebHRYAAECJZLFYFNmwspY+1knDO4TJZrXo5x0n1fWtFXp/xQFl5djNDhGAA6DgA6BYNA+toO8fbKfq/p46npSh/tNX64ctx8wOCwAAoMQq5+qk53rX18JHOqh5aAWlZ+Vqws+7FflOlFbvTzQ7PAAmo+ADoNiE+nnq+wfbq2PtC9u2j/56i15ZGKOcXJ5CAQAAXKt6lb0194G2euO2RvLzdNGBU2m6+8M/9NAXm5jmBZRhFHwAFCsfD2fNGtpSD3WpIUn6cNUhDfpovU6nZpocGQAAQMlltVp0e4tg/f54Zw1pGyqrRfpp+wnd8OYKvbdsvzJz2M0LKGso+AAodjarRU90r6sZ9zSTp4tNaw+eVs8pK7WGoccAAADXxcfdWS/dHK4Fj3RQi9AKOp+dqzcW71G3t6P02654s8MDUIwo+AAwTY/wypr/UHtVr+ip+ORM3f3RH5rw8y4WGgQAALhODYJ8NHdkW00e0ESVvFx15HS67vtko4bNWq+Dp1LNDg9AMaDgA8BUtQK8tPCRDhrYKliGIb2/4qBunb6GRAQAAOA6WSwW9WtaRb8/3lkPdKouZ5tFy/acUvfJUXpt0S6lZGSbHSKAIkTBB4DpPFycNKF/I824p5nKezhr+7Ek9ZqySh+tOqRcu2F2eAAAACVaOVcnjY+sp8VjOqpLnYrKzjX0QdRBdXlzhb7ZGCc7+RZQKlHwAeAweoRX1i+jO6pdDT+dz87VfxfG6JZpq7XjWJLZoQEAAJR41SuW06xhrTRraEuF+XsqMTVTT367Tf2mrVb0kbNmhwegkFHwAeBQAn3c9Pl9rfXaLQ3l5eakbUeTdPN7q/XqTzFKz8oxOzwAAIASr0vdSlo8pqOe6VlX5Vwv5Fu3Tl+jsXO26GRShtnhASgkFsMwStX4veTkZPn4+CgpKUne3t5mhwPgOiSkZOjlBTFauO2EJCnQ202Pdaut/s2qyma1mBwdgOtFn+04aAug7EpIydCbi/dobvRRGYbk7mzTQ11qaHhEdbk528wOD8DfFKTPpuADwOEt252g//ywQ0fPnpck1Q300jM966lj7YomRwbgetBnOw7aAsC2o+f00oKYvKldVSu467le9dS9QaAsFh60AY6Cgg8JC1DqZGTn6tO1hzX19/1KzrgwtSuilr/G3lRbzUIqmBwdgGtBn+04aAsAkmQYhn7celwTFu3WyeQLU7vaVvfT833qq15l/jYAjoCCDwkLUGqdTcvS1GX79enaw8rOvfDnq31NPz1yQy21DvPlCRRQgtBnOw7aAsBfpWflaMbyA3o/6qAyc+yyWqS7WofosZvqqIKni9nhAWUaBR8SFqDUiz2drqnL9mnepmPK+XMr0ZbVKmhU5xrqXLuSrKzxAzg8+mzHQVsAuJy4M+l6/efd+mn7hfUUfdyd9Vi32rqrVYicbOz/A5iBgg8JC1BmHD2brvdXHNScjXHKyrFLkqr7e2pY+2rq36yqPF2dTI4QwJXQZzsO2gLAP1l74LReWrBTu0+mSJLqBHjphT711a6mv8mRAWUPBR8SFqDMiU/O0EerDumr9bFK+XONHy83J93ZMlj3tAlVqJ+nyREC+Dv6bMdBWwD4Nzm5dn21IU5v/bpH59KzJUk9GwbqmZ71VLWCh8nRAWUHBR8SFqDMSsvM0XebjmrW6sM6lJiWdzyilr/uaROqrnUrMQQZcBD02Y6DtgBwtc6lZ+ntJXv12bojshuSq5NVozrX0MhONdjGHSgGFHxIWIAyz243tGxPgj5de0RR+07p4l+6QG833dEyWHe0qMrTKMBk9NmOg7YAUFC7TybrxR93at3BM5KkKuXd9Z/ebOMOFDUKPiQsAP4i9nS6vlwfq7kb43Q6LUuSZLFIHWtV1MBWwepaL0DOjPoBih19tuOgLQBcC8MwtGj7Sb36U4yOJ13Yxj2ilr9e7NtANSqWMzk6oHSi4EPCAuAyMnNytXhnvL5eH6s1B07nHfcv56rbmlfVnS2DVc2ftX6A4kKf7ThoCwDX43xWrqYt36/3VxxUVq5dzjaL7u0QpkduqKVybKABFCoKPiQsAP7F4cQ0zdkYp7kbjyoxNTPveNvqfrqzVbB6hAfK1Yl56EBRos92HLQFgMJwODFNLy+M0e+7EyRJAd6ueq5XffVuVJlpXkAhoeBDwgLgKmXn2vXbrgR9vSFWK/b+/1o/vp4uur15VQ1sFcKoH6CI0Gc7DtoCQGH6bVe8XloQo9gz6ZKk9jX99FLfcNWsxDQv4HpR8CFhAXANjp07r282xOmbjXE68ec8dEnqUNNfd7cO0U31A9jhCyhE9NmOg7YAUNgysnP1/oqDmrZ8vzJzLkzzuq9DdT3ataY8XJjmBVwrCj4kLACuQ06uXcv3nNIXfxzR8r+M+qns46Z72oRqQMtg+ZdzNTdIoBSgz3YctAWAohJ7Ol0vLdip3/6c5hXk46YX+zZQtwaBJkcGlEwUfEhYABSSuDPp+npDrOZsiFNi6oUdvlxsVvVuXFnD2oWpYVUfkyMESi76bMdBWwAoakti4vXijzt17Nx5SdKN9SrphT4NFOzrYXJkQMlCwYeEBUAhy8zJ1aLtJzR7zRFtjTuXd7x1mK+GR1RX17qVZLWyGCFQEPTZjoO2AFAczmfl6t3f9+mDqIPKsRtyc7ZqdNfaGh4RJmemzQNXhYIPCQuAIrQl7pxmrz6khdtOKMd+4U9oNT8P3dchTLe3CJabM7t7AVeDPttx0BYAitO++BQ9N3+H/jh0RpJUJ8BLr/VvqOahFUyODHB8FHxIWAAUgxNJ5/XJmiP68o8jSs7IkST5l3PRsPZhGtQ2VN5uziZHCDg2+mzHQVsAKG6GYei7Tcf06k8xOpueLYtFuqtViJ7sUVc+7uRQwJVQ8CFhAVCM0jJzNHdjnGauPJQ3L93L1Ul3twnV8IgwFngGroA+23HQFgDMciYtSxMW7dLc6KOSJP9yrnqxb331alhZFgvT5YG/o+BDwgLABNm5di3YelzTlx/QvoRUSZK7s02D2obq/o7VKfwAf0Of7ThoCwBmW3vgtJ6dv10HT6VJurCo88s3hyuovLvJkQGOhYIPCQsAE9nthpbuitd7y/Zr69EkSRcKP4PbhmoEhR8gD32246AtADiCzJxcTVt2QNOW71d2riFPF5ueiqyre1qHsjkG8CcKPiQsAByAYRhavueUJi/dm1f48XCx6b4OYRrRsTpr/KDMo892HLQFAEeyNz5FT3+3TZtiz0mSmoWU18TbGqlmJS9zAwMcAAUfEhYADuRi4eftpXu17c/CT3kPZz3YuYYGt63Grl4os+izHQdtAcDR2O2GvvjjiP73yx6lZubIxWbVo11r6oFONdjCHWUaBR8SFgAOyDAMLd4Zrzd/3aP9f67xE+jtpnHdauvWZlVlY6gyyhj6bMdBWwBwVMfPndez32/Xsj2nJEn1K3tr4m2NFF7Fx+TIAHNQ8CFhAeDAcnLtmrf5mCYv2avjSRmSLiQvz/Wqp3Y1/U2ODig+9NmOg7YA4MgMw9D8Lcf00oIYnUvPls1q0chO1fVo11pydWKkNMoWCj4kLABKgIzsXH269rDe/X2/UjJyJF3YkeLpyHqqWamcydEBRY8+23HQFgBKgsTUTL3w4079tO2EJKl2QDm9eXtjNapa3tzAgGJUkD6byY8AYBI3Z5vu71hDK57ooiFtQ2WzWrR0V4J6TI7Sqz/FKCUj2+wQAQAAHIZ/OVe9d1czTb+7mfzLuWhvfKpumbZG//tltzKyc80OD3A4FHwAwGS+ni566eZwLR7TUV3rVlKO3dDMlYd0w1sr9P3moyplAzEBAACuS2TDyvp1bCf1aRykXLuh6csPqM+7q7Tt6DmzQwMcCgUfAHAQNSuV00dDW2rWsJYK8/fUqZRMjZ2zVbfPWKuY48lmhwfAwR06dEhdunRR/fr11bBhQ6WlpZkdEgAUGV9PF707sKlm3NNc/uVctS/hwmifSb/uUVaO3ezwAIfgkAUfEhYAZVmXOpX0y5gIPdmjjtydbdp45Kz6TF2l1xbtUnpWjtnhAXBQQ4cO1csvv6yYmBitWLFCrq6uZocEAEWuR3igfh3bUb0bVVau3dCU3/er33urtfskD8sAhyz4kLAAKOtcnWx6sHNN/f54J/VsGKhcu6EPog7qpklR+m1XvNnhAXAwO3fulLOzsyIiIiRJvr6+cnJyMjkqACgevp4umnpXM029q6kqeDgr5kSy+ry7StOW71eunanxKLscruBDwgIA/6+yj7um3d1cHw9toSrl3XXs3Hnd98lGjfo8WgkpGWaHB+AqRUVFqU+fPgoKCpLFYtH8+fMvOWfatGkKCwuTm5ubmjdvrpUrV1719fft26dy5cqpb9++atasmV577bVCjB4ASobejYK0eGxH3VgvQNm5hib+skcD3l+rI6eZMYKyqcAFHxIWACh+N9QN0JJxHfVAp+qyWS36ecdJ3TQpSt9Fs6gzUBKkpaWpcePGmjp16mXfnzNnjsaMGaNnn31WmzdvVkREhCIjIxUbG5t3TvPmzRUeHn7J6/jx48rOztbKlSv13nvvae3atVqyZImWLFlSXF8PABxGJS83zRzcXG/c1kjlXJ208chZRb6zUl+tjyVnQplT4KEzFxOWYcOG6dZbb73k/YsJy7Rp09S+fXu9//77ioyMVExMjEJCQiRdSFgyMzMv+d1ff/01L2HZsmWLKlWqpB49eqhly5a66aabruHrAUDp4eHipPGR9XRz4yp68rut2nEsWY/N3aoF247rtVsaKqi8u9khAriCyMhIRUZGXvH9SZMm6b777tPw4cMlSZMnT9bixYs1ffp0TZgwQZIUHR19xd+vWrWqWrZsqeDgYElSz549tWXLlivmT5mZmflyseRk1roAUHpYLBbd3iJYbar76bG5W7X+0BmNn7ddS2Li9b9bG6miF0uGoGwo8AifyMhIvfLKK+rfv/9l3/9rwlKvXj1NnjxZwcHBmj59et450dHR2rFjxyWvoKCgfAmLq6trXsJyJZmZmUpOTs73AoDSrH6Qt+Y/2F5P9qgjFyerlu85pW5vR+lrnlwBJVJWVpaio6PVrVu3fMe7deumNWvWXNU1WrZsqfj4eJ09e1Z2u11RUVGqV6/eFc+fMGGCfHx88l4XC0UAUJoE+3ro6xFt9GzPenKxWfX77gT1mBylpTGsh4iyoVDX8CFhAYDi4WSz6sHONbXo0Q5qFlJeqZk5enredt33yUbW9gFKmMTEROXm5iogICDf8YCAAJ08efKqruHk5KTXXntNHTt2VKNGjVSrVi317t37iuePHz9eSUlJea+4uLjr+g4A4KisVotGdKyuHx9pr7qBXjqdlqXhn27UM99vZ/dTlHqFWvAhYQGA4lWzkpfmjmynZ3rWzXty1f3tKC3afsLs0AAUkMViyfezYRiXHPsnkZGR2r59u3bs2KFJkyb947murq7y9vbO9wKA0qxuoLfmP9ReIyLCJElf/hGrXlNWaWvcOXMDA4pQkezSRcICAMXHZrXo/o41tOCRDmoQ5K2z6dl68ItNGv31ZiWdzzY7PAD/wt/fXzab7ZKHYwkJCZc8RAMAXDs3Z5ue7VVfXw5vrUBvNx1KTNOt09ewfTtKrUIt+JCwAIB56gR66fsH2+uRG2rKZrXohy3H1fOdldp4+IzZoQH4By4uLmrevPklu2otWbJE7dq1MykqACi92tX01+IxHdWrYWXl2C9s3373h+t0Ium82aEBhapQCz4kLABgLhcnqx7rVkdzR7ZViK+Hjp07rzveX6u3l+xVTq7d7PCAMis1NVVbtmzJ24ji0KFD2rJlS9626+PGjdOHH36ojz/+WLt27dLYsWMVGxurkSNHmhg1AJRePh7OmnpXU028rZE8XGxad/CMekxeqZ+ZFo9SpMDbsqempmr//v15P19MWHx9fRUSEqJx48Zp0KBBatGihdq2basPPviAhAUAilmzkAr66dEOeuGHnZq3+Zje+W2fVu1P1Dt3NlHVCh5mhweUORs3blSXLl3yfh43bpwkaciQIZo9e7YGDBig06dP6+WXX9aJEycUHh6uRYsWKTQ01KyQAaDUs1gsuqNFsFpW89WYrzdr69Ekjfpikwa2CtHzvevL3cVmdojAdbEYBdzDd/ny5fkSlosuJiySNG3aNE2cODEvYXn77bfVsWPHQgn43yQnJ8vHx0dJSUms5wMAkuZvPqbn5u9QamaOvN2c9ObtjdWtQaDZYQH02Q6EtgBQ1mXn2jVpyV7NWHFAhiHVrFRO7w5sqnqV+ZsIx1KQPrvABR9HR8ICAJeKO5Ouh7/anLcTxX0dwvRUj7pycSqStfuBq0Kf7ThoCwC4YPX+RI2ds0UJKZlycbLqmci6GtKuWoE2IQKKUkH6bDJ9ACgDgn09NPeBtrqvw4WtSD9adUi3v79WcWfSTY4MAADAcbSv6a+fR0fohrqVlJVj14sLYjTi02idTcsyOzSgwCj4AEAZ4eJk1X9619cHg5rL281JW+POqdeUlVq2O8Hs0AAAAByGXzlXfTSkhV7sU18uNquW7opXzykr9cfB02aHBhQIBR8AKGO6NQjUT49GqHFweSVn5GjY7A2atGSvcu2laoYvAADANbNYLBraPkzzHmyn6v6eOpGUoYEz1+mdpfvImVBiUPABgDIo2NdD3zzQRoPaXNgBaMpv+3Tv7A06l85wZQAAgIvCq/howSMddGuzqrIb0ttL9+ruD9cpPjnD7NCAf0XBBwDKKFcnm/7bL1yT7mgsN2erVuw9pV5TVmnHsSSzQwMAAHAYnq5OeuuOxnp7QGN5uti07uAZRb6zUsv3MC0ejo2CDwCUcf2bVdW8Ue0V4uuhY+fO67YZa/Tj1uNmhwUAAOBQbmlaVQse6aD6lb11Ji1LQ2dt0IRFu5Sdazc7NOCyKPgAAFQ/yFsLHu6gTrUrKiPbrke/2qzXf97NHHUAAIC/qF6xnOY92E6D216YFv9+1EHd8f5aHT3LzqdwPBR8AACSJB8PZ308tKUe6FRdkjRjxQHd98kGJZ3PNjkyAAAAx+HmbNPLN4drxj3N5OXmpM2x59TznZX6dedJs0MD8qHgAwDIY7NaND6ynt65s4lcnaxavueUbnlvtQ6eSjU7NAAAAIfSI7yyFv1l59P7P4vWSwt2KiuHKV5wDBR8AACXuLlJFX03qp2CfNx0MDFNt0xbozUHEs0OCwAAwKEE+3po7gNtNbxDmCRp1urDum3GGsWeZooXzEfBBwBwWeFVfDT/4fZqGlJeSeezNfij9fpqfazZYQEAADgUFyernutdXx8ObiEfd2dtO5qkXlNW6uftJ8wODWUcBR8AwBVV8nLTVyPaqG/jIOXYDY2ft13/XRjDYs4AAAB/c2P9AC0aHaFmIeWVkpmjUV9s0gs/7FBmTq7ZoaGMouADAPhHbs42vXNnE429sbYk6aNVh/TAZ9FKz8oxOTIAAADHUqW8u+Y80DZvE4xP1h7RbdPX6sjpNJMjQ1lEwQcA8K8sFotG31hL7w5sKhcnq5buitedH6xTQkqG2aEBAAA4FGebVeMj6+njoS1UwcNZ248lqfeUVfppG1O8ULwo+AAArlqfxkH6akRrVfC4MD/9lvfWaH9CitlhAQAAOJwb6l6Y4tUitIJSMnP00Jeb9J/5O5SRzRQvFA8KPgCAAmke6qvvH2yvan4eOnbuvPpPW6O1B06bHRYAAIDDqezjrq/ub6NRnWtIkj5bd0S3Tl+jw4lM8ULRo+ADACiwav6emvdge7UIraDkjBwN+Xi9Fmw9bnZYAAAADsfZZtVTPepq1rCWquDhrJ3Hk9X73VXkTihyFHwAANfE19NFnw9vrV4NKysr165Hvtqsj1YdMjssAAAAh9SlTiUtGh2hltUqKDUzR498tVnPfL+dKV4oMhR8AADXzM3ZpncHNtXQdtUkSf9dGKPXFu2SnW3bAQAALlHZx11fjWijh7rUkMUifflHrPq9t1oHTqWaHRpKIQo+AIDrYrVa9EKf+nqqR11J0gdRBzXumy3KyrGbHBkAAIDjcbJZ9UT3uvpkWCv5ebpo98kU9Xl3lb7ffNTs0FDKUPABAFw3i8WiUZ1raNIdjeVktWj+luMa/ulGpWflmB0aAACAQ+pYu6J+Hh2httX9lJ6Vq7FzturxuVvJn1BoKPgAAApN/2ZV9dHQlnJ3tilq7ynd/eEfOpeeZXZYAAAADqmSt5s+H95aY26sJatF+jb6qPq8u0q7TiSbHRpKAQo+AIBC1al2RX0xorV83J21Ofacbp+xVieTMswOCwAAwCHZrBaNubG2vhjeRgHerjpwKk03v7dan687IsNgXURcOwo+AIBC1yykguaObKsAb1ftS0jVrdPX6CCLEQIAAFxR2xp+WvRohLrUqaisHLuem79DD325SUnp2WaHhhKKgg8AoEjUDvDStyPbKczfU8fOndftM9Zq5/Eks8MCAABwWH7lXPXRkJZ6tmc9OVktWrT9pHpOWanoI2fMDg0lEAUfAECRCfb10NyRbdUgyFun07J05wfrtPEwCQsAAMCVWK0WjehYXd+NaqdQPw8dO3ded7y/Tu/+tk+5dqZ44epR8AEAFCn/cq766v42almtglIycjToo/WK2nvK7LAAAAAcWuPg8lr4SAf1axKkXLuht5bs1d0frtOJpPNmh4YSgoIPAKDIebs569N7W6tT7Yo6n52r+z7ZoJ+3nzA7LAAAAIfm5easyXc21aQ7GsvDxaZ1B8+ox+SV+mUHeRT+HQUfAECxcHexaebgFurVsLKycw099OUmzd0YZ3ZYAAAADq9/s6r66dEINarqo6Tz2Rr5+SY9/d02pWflmB0aHBgFHwBAsXFxsmrKwKa6o0VV2Q3piW+36bO1h80OCwAAwOGF+Xvq25HtNKpzDVks0tcb4tR7yiptP8qmGLg8Cj4AgGJls1r0ev9GGtqumiTpPz/s1IwVB8wNCgAAoARwcbLqqR519cXw1gr0dtPBxDTdMm213lu2nwWdcQkKPgCAYme1WvRCn/p6uEtNSdLrP+/WpF/3yDBIVAAAAP5Nuxr++mVMhHo2DFSO3dAbi/dowPtrFXcm3ezQ4EAo+AAATGGxWPR49zp6skcdSdKU3/frlZ92UfQBAAC4CuU9XPTeXc301u2NVc7VSRuPnFXkOys1d2Mc+RQkUfABAJjswc419VLfBpKkj1Yd0jPf75CdIckAAAD/ymKx6NbmVfXz6Ai1rFZBqZk5euLbbXrgs2glpmaaHR5MRsEHAGC6Ie2qaeJtjWS1SF+tj9Vjc7cqJ9dudlgAAAAlQrCvh76+v62e7FFHzjaLfo2JV7e3o9i+vYyj4AMAcAh3tAjWO3c2lZPVou83H9PDX25WVg5FHwAAgKths1r0YOea+uGhDqob6KUzaVka+fkmjZ2zRUnns80ODyag4AMAcBh9Ggdp+j3N5WKz6pedJ3X/ZxuVkZ1rdlgAAAAlRv0gb/3wcHs92LmGrBbp+83H1O3tFVq2O8Hs0FDMKPgAABzKTfUD9NHQFnJztmr5nlMaOmu9UjNzzA4LAACgxHB1sunJHnU1d2Q7Vff3VHxypobN3qDHvtmqpHRG+5QVFHwAAA4nolZFfXpva5VzddK6g2d094d/6Fx6ltlhAQAAlCjNQyto0egIDe8QJotF+m7TUXWbvEK/7Yo3OzQUAwo+AACH1CrMV1+OaK3yHs7aGndOd36wTqdS2G0CAACgINycbXqud319O7Jt3mif+z7ZqDFfb9aZNB6olWYUfAAADqtR1fKac39bVfRy1e6TKRrw/lodP3fe7LAAAABKnOahvlo0OkL3d6wuq0Wav+W4bpy0Qj9sOSbDMMwOD0WAgg8AwKHVCfTS3Afaqkp5dx1MTNPtM9bqwKlUs8MCAAAocdycbXqmZz19/2D7vJ28Rn+9RcM/2chDtVKIgg8AwOFV8/fU3D+HIR87d153zFirHceSzA4LAACgRGocXF4/PtxB426qLWebRb/tTtBNk1Zo1upDyrUz2qe0oOADACgRgsq7a+7Itgqv4q3TaVm684N1WnfwtNlhAQAAlEguTlY92rWWfno0Qs1CyistK1cvLYhR/2mrFXM82ezwUAgo+AAASgy/cq76akQbtanuq9TMHA3+eL2WxLDLBAAAwLWqHeClb0e20yv9wuXl6qStR5PUZ+oqTVi0S+lZOWaHh+tAwQcAUKJ4uTlr9rBW6lY/QFk5do38PFpzNsSaHRYAAECJZbVadE+bUP32WCf1alhZuXZD70cd1E2TovTrzpNmh4drRMEHAFDiuDnbNO3uZrq9eVXl2g099d12TVqylx0mUKa9/fbbatCggerXr69HH32U/x4AAAVWydtN793dTB8PbaEq5d117Nx53f9ZtIZ/skFHz6abHR4KiIIPAKBEcrJZNfG2Rnr0hpqSpCm/7dOT325Tdq7d5MiA4nfq1ClNnTpV0dHR2r59u6Kjo7Vu3TqzwwIAlFA31A3Q0nGd9GDnGnKyWrR0V4JunLRC7y3br6wccq2SwiELPjyhAgBcDYvFonHd6ui1WxrKapHmRh/VfZ9sVGom881R9uTk5CgjI0PZ2dnKzs5WpUqVzA4JAFCCubvY9GSPuvp5dIRahfkqI9uuNxbvUY93orRy3ymzw8NVcLiCD0+oAAAFdVfrEM0c3ELuzjZF7T2l26avYdgxHEpUVJT69OmjoKAgWSwWzZ8//5Jzpk2bprCwMLm5ual58+ZauXLlVV+/YsWKevzxxxUSEqKgoCDdeOONqlGjRiF+AwBAWVUrwEtz7m+jSXc0ln85Vx08laZBH63Xg19E6/i582aHh3/gcAUfiSdUAICC61ovQF/f30b+5Vy0+2SKbp66WhsOnzE7LECSlJaWpsaNG2vq1KmXfX/OnDkaM2aMnn32WW3evFkRERGKjIxUbOz/L0jevHlzhYeHX/I6fvy4zp49q4ULF+rw4cM6duyY1qxZo6ioqOL6egCAUs5isah/s6r6/fFOGta+mqwWadH2k+r61oVpXpk5uWaHiMsocMGHJ1QAAEfVOLi8fni4g+pX9tbptCzdNXMdO3jBIURGRuqVV15R//79L/v+pEmTdN9992n48OGqV6+eJk+erODgYE2fPj3vnOjoaO3YseOSV1BQkJYuXaqaNWvK19dX7u7u6tWrFyOkAQCFztvNWS/0aaCFj0SoZbUKOp+dqzcW71H3t6P0++54s8PD3xS44MMTKgCAI6tS3l3fjmqrng0DlZ17YQevlxfEKIfFnOGgsrKyFB0drW7duuU73q1bN61Zs+aqrhEcHKw1a9YoIyNDubm5Wr58uerUqXPF8zMzM5WcnJzvBQDA1aof5K1vHmiryQOaqJKXqw6fTte9szfq3tkbdDgxzezw8Cengv5CZGSkIiMjr/j+X59QSdLkyZO1ePFiTZ8+XRMmTJB04QnVlcydOzfvCZWkvCdUHTt2vOz5mZmZyszMzPuZhAUA4OHipKkDm2lKwD5NXrpPH68+pO3HzmnKwKaq7ONudnhAPomJicrNzVVAQEC+4wEBATp58uRVXaNNmzbq2bOnmjZtKqvVqq5du6pv375XPH/ChAl66aWXrituAEDZZrFY1K9pFd1YP0Dv/rZPH606pN93J2jVvkTd2yFMD99QU+VcC1xyQCEq1DV8zHhCNWHCBPn4+OS9goODr+s7AABKB6vVojE31taMe5qpnKuTNhw+q15TVmn5ngSzQwMuy2Kx5PvZMIxLjv2TV199Vbt27dLOnTs1ZcqUf/zd8ePHKykpKe8VFxd3zXEDAMq2cq5OGt+znn4Z01Eda1dUVq5dM1Yc0A1vLtd30Udlt7PrtlkKteBT2E+oGjVqpBo1avzjEyoSFgDAP+kRXlkLH+mgBkHeOpOWpaGzNmjiL7uZ4gWH4e/vL5vNdkmulJCQcElOVVhcXV3l7e2d7wUAwPWoWamcPhnWUh8ObqFQPw8lpGTqsblb1X/6Gm2JO2d2eGVSkezSVZxPqEhYAAD/ppq/p74b1U6D24ZKkqYtP6Db31+rA6dSTY4MkFxcXNS8eXMtWbIk3/ElS5aoXbt2JkUFAEDBWSwW3Vg/QL+O7ainetSVh4tNW+LOqd97q/XYN1uVkJxhdohlSqEWfMx4QgUAwNVwc7bp5ZvD9d5dzeTl6qTNsefU852Vmhl1ULkMNUYRS01N1ZYtW7RlyxZJ0qFDh7Rly5a8TS3GjRunDz/8UB9//LF27dqlsWPHKjY2ViNHjjQxagAAro2rk02jOtfQssc7q3+zKpKk7zYdVZc3l2vacrZxLy6FWvDhCRUAwNH1alRZi8d2VEQtf2Xm2PXqol264/21OshoHxShjRs3qmnTpmratKmkCwWepk2b6vnnn5ckDRgwQJMnT9bLL7+sJk2aKCoqSosWLVJoaKiZYQMAcF0CvN006Y4m+v7BdmoSXF5pWbma+Mse3TQpSot3npRh8NCtKFmMAt7h1NRU7d+/X5LUtGlTTZo0SV26dJGvr69CQkI0Z84cDRo0SDNmzFDbtm31wQcfaObMmdq5c2exJC3Jycny8fFRUlIS07sAAFdkGIbmbIjTKz/tUmpmjlydrBrVuYZGdqohN2eb2eGVCfTZjoO2AAAUNbvd0Pwtx/T6z7uVkHJhp+32Nf30fO8GqhPoZXJ0JUdB+uwCF3yWL1+uLl26XHJ8yJAhmj17tiRp2rRpmjhxok6cOKHw8HC9/fbbV9xWvbCRsAAACuLYufN6+rttWrkvUZJUpby7/tO7nro3CCzQ+nMoOPpsx0FbAACKS1pmjqYt36+ZKw8pK8cuq0W6p02oxt5YWxU8XcwOz+EVacHH0ZGwAAAKyjAMLdp+Uq/+FKPjSRcWE2xf00//6V1fdQPpS4oKfbbjoC0AAMUt7ky6Xlu0Sz/vuLAGsI+7s8bdVFt3tw6Rk61I9pcqFSj4kLAAAK5BelaOZiw/oBlRB5WVY5fFIvVpFKQxN9ZS9YrlzA6v1KHPdhy0BQDALGsOJOrlBTHafTJFklQnwEsv9KmvdjX9TY7MMVHwIWEBAFyHuDPpev3n3fpp+wlJktUi3dqsqh7tWkvBvh4mR1d60Gc7DtoCAGCmnFy7vtoQp7d+3aNz6dmSpO4NAvRcr/rkXn9DwYeEBQBQCHYcS9LbS/bqt90JkiQnq0V9Gwfp/k7VmepVCOizHQdtAQBwBOfSs/T2kr36/I9Y5doNuThZ9UDH6hrVuYY8XJzMDs8hUPAhYQEAFKJNsWc16de9WrU/Me9Y5zoV9UDHGmpT3ZfFna8RfbbjoC0AAI5kz8kUvbxwp1bvPy1JquzjpvE966lPo8plPu+i4EPCAgAoAlvjzumDqIP6eccJ2f/sPRsEeWtI22rq2ySI7dwLiD7bcdAWAABHYxiGFu88qVd+2qWjZ89LklqF+eqlvg1Ur3LZ7aso+JCwAACK0JHTafpw5SF9szFOmTl2SVJ5D2cNaBGse9qEMtf8KtFnOw7aAgDgqDKyc/VB1EFNW75fGdkXtnEf1CZU426qIx8PZ7PDK3YUfEhYAADF4Gxalr7ZGKfP1h3Je/JksUidalfUwFYh6lq3EtuK/gP6bMdBWwAAHN2xc+f12k+78jbVqODhrKd61NUdLYJltZadaV4UfEhYAADFKNduaNnuBH2y9rBW7vv/dX4CvF01oEWwBrQKUZXy7iZG6Jjosx0HbQEAKCnW7E/Uiwt2am98qiSpcXB5/ffmBmpUtby5gRUTCj4kLAAAkxxOTNNXG2L17cajOp2WJen/R/3c2TJEXetVkjOjfiTRZzsS2gIAUJJk59r16dojenvJXqVm5shikQa2CtET3eqogqeL2eEVKQo+JCwAAJNl5uRqSUy8vvwjVmsOnM47XtHLVbc3r6o7W4YoxK9sr/VDn+04aAsAQEmUkJyhCT/v1vebj0m6MM3r6ci6ur156Z3mRcGHhAUA4EAOJ6ZpzsY4zd14VImpmXnHI2r5665WIbqxfkCZHPVDn+04aAsAQEm2/tAZPf/DDu0+mSJJahZSXv/tF64GQT4mR1b4KPiQsAAAHFB2rl2/7YrXl+vjtHLfKV3sgf3LuerOlsG6q3WIgsrQWj/02Y6DtgAAlHQ5uXbNXnNYby/Zq7SsXFkt0uC21fRYt9rycis9u3lR8CFhAQA4uLgz6fpqfay++cuoH6tFurFegAa3rab2Nf1ksZTOocgX0Wc7DtoCAFBanEzK0H9/itFP2y7s5lXJy1Uv9Gmgng0DS0VuRcGHhAUAUEJk59r16854fbbusNYdPJN3vHpFTw1rV039m1WVp6uTiREWHfpsx0FbAABKm6i9p/T8Dzt0+HS6JKlznYr6783hCvYt2WsoUvAhYQEAlED74lP02bojmrfpmFIzcyRJXm5OGtgqRIPbhqpqhZKdoPwdfbbjoC0AAKVRRnaupi0/oBnLDygr1y43Z6se7VpLIyKql9j1Eyn4kLAAAEqw1MwcfbsxTrPXHM57KmW1SJHhlfVAp+pqVLW8uQEWEvpsx0FbAABKs/0JqXpu/va80dR1A700oX9DNQ2pYHJkBUfBh4QFAFAK2O2Glu1J0KzVh7Vqf2Le8bbV/fRAp+rqVLtiiZ6LTp/tOGgLAEBpZxiG5m06pld+itHZ9GxZLNLgNqF6vHudErWoMwUfEhYAQCmz60SyZkYd1I9bjyvHfqHrrhvopYdvqKnI8MqyWUte4Yc+23HQFgCAsuJMWpZe+SlG8zYdkyQFervpv/3CdVP9AJMjuzoUfEhYAACl1LFz5/XxqkP6an2s0rNyJUk1Knrq4Rtqqk+jIDmVoPno9NmOg7YAAJQ1q/Yl6tn523Xkz+nzvRpV1ot9Gqiil6vJkf0zCj4kLACAUu5cepZmrT6sWasPKTnjwgLPoX4eeuSGWurXpGQUfuizHQdtAQAoizKyc/X20r36cOUh5doN+bg76z+96+vWZlUcdto8BR8SFgBAGZGSka1P1x7RR6sO6UxalqQLI37G3lRbPcMry+rAU73osx0HbQEAKMt2HEvSk99uU8yJZElSRC1/Tejf0CF3SKXgQ8ICAChj0rNy9OnaI5qx4oDOpWdLurDGz+Pd6qhrvUoO+ZSKPttx0BYAgLIuO9euD1ce0uSle5WZY5eni01P96ynu1uFONQDNAo+JCwAgDIqJSNbH686rA9XHlRK5oWpXq2q+Wp8z7oOt/UofbbjoC0AALjg4KlUPfntNm08claS1Ka6r/53ayOF+nmaHNkFFHxIWAAAZdy59CzNWHFQs1YfUmaOXZLUq2FlPdG9jqr5l7yEBUWLtgAA4P/Z7YY+WXtYE3/Zo/PZuXJ3tumpHnU0uG0100f7UPAhYQEAQJJ0Ium8Jv26V99uOirDkJysFg1uW02ju9aSj4ezqbHRZzsO2gIAgEvFnk7XU99t09qDpyVJrcN89cZtjRXiZ97aPgXpsx1/Cw8AAHDNKvu4643bG+vn0RHqUqeicuyGPl59SJ3fXKbP1x1Rrr1UPfcBAAAoNCF+HvpieGv9t1+4PFxs+uPQGfV4J0qfrT0sewnIoSj4AABQBtQN9NasYa302X2tVKtSOZ1Nz9Zz83eo15SVWnMg0ezwAAAAHJLVatGgNqH6ZXRHtQ7zVXpWrv7zw07d89EfOnbuvNnh/SMKPgAAlCERtSrq59EReqlvA/m4O2v3yRTdNfMPjf56sxKSM8wODwAAwCGF+HnoqxFt9GKf+nJztmrNgdPq8XaUvo0+KkddKYeCDwAAZYyTzaoh7app+eOdNahNqKwW6Yctx9X1rRWatfqQcnLtZocIAADgcKxWi4a2D9PPozuqWUh5pWTm6PG5W3X/Z9FKTM00O7xLUPABAKCMquDpov/2C9cPD3VQ4+ALSctLC2LUd+pqbY49a3Z4AAAADinM31NzR7bTkz3qyNlm0ZKYeHV7O0qLd540O7R8KPgAAFDGNazqo+9HtdNrtzSUj7uzYk4kq//0NXppwU6lZeaYHR4AAIDDsVkterBzTf34cAfVDfTSmbQsPfBZtJ78dqtSHSR/ouADAABktVp0V+sQ/f5YJ/VvWkWGIc1afVjd3o7Ssj0JZocHAADgkOpV9tYPD7fXyE41ZLFI32w8qsh3orTx8BmzQ6PgAwAA/p9fOVdNGtBEn9zbSlUruOvYufMaNmuDRn+9WacdcG46AACA2VydbHo6sq7m3N9WVcq7K+7Med3x/lpN/GW3snLMWxuRgg8AALhEp9oV9evYjhreIUxWi/TLjpMOMzwZAADAEbUK89UvYyJ0W/OqshvS/M3HdD4717R4nEz7ZAAA4NA8XJz0XO/66tskSPsTUhXq52l2SAAAAA7Ny81Zb97eWDfWq6TyHi7ycXc2LRYKPgAA4B81qlpejaqWNzsMAACAEqNHeGWzQ2BKFwAAAAAAQGlDwQcAAAAAAKCUoeADAAAAAABQylDwAQAAAAAAKGUo+AAAAAAAAJQyFHwAAAAAAABKGQo+AAAAAAAApQwFHwAAAAAAgFKGgg8AAAAAAEApQ8EHAAAAAACglHEyO4DCZhiGJCk5OdnkSAAAwD+52Fdf7LthHvInAABKhoLkT6Wu4JOSkiJJCg4ONjkSAABwNVJSUuTj42N2GGUa+RMAACXL1eRPFqOUPVaz2+06fvy4vLy8ZLFYCvXaycnJCg4OVlxcnLy9vQv12riAe1w8uM/Fg/tc9LjHxaOo7rNhGEpJSVFQUJCsVmaZm6mo8if+Gy0e3OfiwX0uetzj4sF9Lh6OkD+VuhE+VqtVVatWLdLP8Pb25j+MIsY9Lh7c5+LBfS563OPiURT3mZE9jqGo8yf+Gy0e3OfiwX0uetzj4sF9Lh5m5k88TgMAAAAAAChlKPgAAAAAAACUMhR8CsDV1VUvvPCCXF1dzQ6l1OIeFw/uc/HgPhc97nHx4D7jWvFvp3hwn4sH97nocY+LB/e5eDjCfS51izYDAAAAAACUdYzwAQAAAAAAKGUo+AAAAAAAAJQyFHwAAAAAAABKGQo+AAAAAAAApQwFn6s0bdo0hYWFyc3NTc2bN9fKlSvNDqnEmjBhglq2bCkvLy9VqlRJ/fr10549e/KdYxiGXnzxRQUFBcnd3V2dO3fWzp07TYq4dJgwYYIsFovGjBmTd4z7XDiOHTume+65R35+fvLw8FCTJk0UHR2d9z73+frk5OToueeeU1hYmNzd3VW9enW9/PLLstvteedwjwsuKipKffr0UVBQkCwWi+bPn5/v/au5p5mZmXrkkUfk7+8vT09P9e3bV0ePHi3GbwFHR/5UeMifzEH+VHTIn4oW+VPRKHH5k4F/9fXXXxvOzs7GzJkzjZiYGGP06NGGp6enceTIEbNDK5G6d+9uzJo1y9ixY4exZcsWo1evXkZISIiRmpqad87rr79ueHl5Gd99952xfft2Y8CAAUblypWN5ORkEyMvudavX29Uq1bNaNSokTF69Oi849zn63fmzBkjNDTUGDp0qPHHH38Yhw4dMpYuXWrs378/7xzu8/V55ZVXDD8/P2PhwoXGoUOHjLlz5xrlypUzJk+enHcO97jgFi1aZDz77LPGd999Z0gyvv/++3zvX809HTlypFGlShVjyZIlxqZNm4wuXboYjRs3NnJycor528ARkT8VLvKn4kf+VHTIn4oe+VPRKGn5EwWfq9CqVStj5MiR+Y7VrVvXePrpp02KqHRJSEgwJBkrVqwwDMMw7Ha7ERgYaLz++ut552RkZBg+Pj7GjBkzzAqzxEpJSTFq1aplLFmyxOjUqVNewsJ9LhxPPfWU0aFDhyu+z32+fr169TLuvffefMf69+9v3HPPPYZhcI8Lw98Tlqu5p+fOnTOcnZ2Nr7/+Ou+cY8eOGVar1fjll1+KLXY4LvKnokX+VLTIn4oW+VPRI38qeiUhf2JK17/IyspSdHS0unXrlu94t27dtGbNGpOiKl2SkpIkSb6+vpKkQ4cO6eTJk/nuuaurqzp16sQ9vwYPPfSQevXqpRtvvDHfce5z4fjxxx/VokUL3X777apUqZKaNm2qmTNn5r3Pfb5+HTp00G+//aa9e/dKkrZu3apVq1apZ8+ekrjHReFq7ml0dLSys7PznRMUFKTw8HDuO8ifigH5U9Eifypa5E9Fj/yp+Dli/uRU6FcsZRITE5Wbm6uAgIB8xwMCAnTy5EmToio9DMPQuHHj1KFDB4WHh0tS3n293D0/cuRIscdYkn399dfatGmTNmzYcMl73OfCcfDgQU2fPl3jxo3TM888o/Xr1+vRRx+Vq6urBg8ezH0uBE899ZSSkpJUt25d2Ww25ebm6tVXX9XAgQMl8W+5KFzNPT158qRcXFxUoUKFS86hfwT5U9Eifypa5E9Fj/yp6JE/FT9HzJ8o+Fwli8WS72fDMC45hoJ7+OGHtW3bNq1ateqS97jn1ycuLk6jR4/Wr7/+Kjc3tyuex32+Pna7XS1atNBrr70mSWratKl27typ6dOna/DgwXnncZ+v3Zw5c/T555/ryy+/VIMGDbRlyxaNGTNGQUFBGjJkSN553OPCdy33lPuOv+K/y6JB/lR0yJ+KB/lT0SN/Mo8j5U9M6foX/v7+stlsl1TbEhISLqncoWAeeeQR/fjjj1q2bJmqVq2adzwwMFCSuOfXKTo6WgkJCWrevLmcnJzk5OSkFStWaMqUKXJycsq7l9zn61O5cmXVr18/37F69eopNjZWEv+eC8MTTzyhp59+WnfeeacaNmyoQYMGaezYsZowYYIk7nFRuJp7GhgYqKysLJ09e/aK56DsIn8qOuRPRYv8qXiQPxU98qfi54j5EwWff+Hi4qLmzZtryZIl+Y4vWbJE7dq1Mymqks0wDD388MOaN2+efv/9d4WFheV7PywsTIGBgfnueVZWllasWME9L4CuXbtq+/bt2rJlS96rRYsWuvvuu7VlyxZVr16d+1wI2rdvf8m2uHv37lVoaKgk/j0XhvT0dFmt+bsrm82Wt60o97jwXc09bd68uZydnfOdc+LECe3YsYP7DvKnIkD+VDzIn4oH+VPRI38qfg6ZPxX6MtCl0MVtRT/66CMjJibGGDNmjOHp6WkcPnzY7NBKpFGjRhk+Pj7G8uXLjRMnTuS90tPT8855/fXXDR8fH2PevHnG9u3bjYEDB7JFYCH46y4ThsF9Lgzr1683nJycjFdffdXYt2+f8cUXXxgeHh7G559/nncO9/n6DBkyxKhSpUretqLz5s0z/P39jSeffDLvHO5xwaWkpBibN282Nm/ebEgyJk2aZGzevDlvy+yruacjR440qlataixdutTYtGmTccMNN7AtO/KQPxUu8ifzkD8VPvKnokf+VDRKWv5Ewecqvffee0ZoaKjh4uJiNGvWLG8LTBScpMu+Zs2alXeO3W43XnjhBSMwMNBwdXU1OnbsaGzfvt28oEuJvycs3OfCsWDBAiM8PNxwdXU16tata3zwwQf53uc+X5/k5GRj9OjRRkhIiOHm5mZUr17dePbZZ43MzMy8c7jHBbds2bLL/i0eMmSIYRhXd0/Pnz9vPPzww4avr6/h7u5u9O7d24iNjTXh28BRkT8VHvIn85A/FQ3yp6JF/lQ0Slr+ZDEMwyj8cUMAAAAAAAAwC2v4AAAAAAAAlDIUfAAAAAAAAEoZCj4AAAAAAAClDAUfAAAAAACAUoaCDwAAAAAAQClDwQcAAAAAAKCUoeADAAAAAABQylDwAQAAAAAAKGUo+AAAAAAAAJQyFHwAAAAAAABKGQo+AAAAAAAApQwFHwAAAAAAgFLm/wBXtfgSVPWMPAAAAABJRU5ErkJggg==", + "text/plain": [ + "<Figure size 1400x400 with 2 Axes>" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], "source": [ "plot_convergences(res, np.ones(4))" ] @@ -246,10 +280,21 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 17, "metadata": {}, - "outputs": [], - "source": [] + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[np.float64(1.0), np.float64(1.0), np.float64(1.0), np.float64(1.0)]\n" + ] + } + ], + "source": [ + "res = gauss_seidel_r(A, b, x0, w=1, n=100)\n", + "print(res[-1])" + ] }, { "cell_type": "markdown", @@ -265,10 +310,29 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 54, "metadata": {}, - "outputs": [], - "source": [] + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[np.float64(1.000000000000001), np.float64(0.9999999999999989), np.float64(1.5754107389499148), np.float64(0.4245892610500852)]\n", + "[np.float64(5.670568200880977e+92), np.float64(-1.8901894002936588e+92), np.float64(3.1679574348921706e+93), np.float64(-3.7350142549802686e+93)]\n" + ] + } + ], + "source": [ + "np.random.seed(308)\n", + "A = np.random.randint(10, size=(4,4))\n", + "b = A.sum(axis=1)\n", + "x0 = np.random.random(4)\n", + "\n", + "r1 = gauss_seidel_r(A, b, x0, w=0.2, n=1000)\n", + "r2 = gauss_seidel_r(A, b, x0, w=1, n=1000)\n", + "print(r1[-1])\n", + "print(r2[-1])" + ] }, { "cell_type": "markdown", @@ -283,15 +347,55 @@ "cell_type": "code", "execution_count": null, "metadata": {}, - "outputs": [], - "source": [] + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Itération du minimum : 4 5\n" + ] + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAABJMAAAF0CAYAAACXCngZAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjMsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvZiW1igAAAAlwSFlzAAAPYQAAD2EBqD+naQAAUh1JREFUeJzt3Xt8VNW9///3zCSTGyFAkGAIBFCrYLjYgBaUS6CCqeCltV6ogKi1ShQhtVjRUy0/NR61PtQj4tGfBan1QK1AbaUiiCgKKoRQUYqIRYNc5ZoQyHXW949kJplkkswkc9khr+fjkYfOnpU9a3Zs9+fx2Z/1WTZjjBEAAAAAAADgB3ukJwAAAAAAAIC2g2QSAAAAAAAA/EYyCQAAAAAAAH4jmQQAAAAAAAC/kUwCAAAAAACA30gmAQAAAAAAwG8kkwAAAAAAAOA3kkkAAAAAAADwG8kkAEDE7NmzR926ddOvf/3rSE8FAADA8oidYBUkk4AQWLhwoWw2W6M/a9eujfQUw2b06NHKyMgI62d+8803stlsWrhwYVg/93R20003qXfv3l7HevfurZtuuqnF56ysrNT111+vMWPG6Mknn2zdBAEApwViqOa1Js7Zu3evHnroIW3ZsqXBew899JBsNlvrJwhJxE44/UVFegLA6WzBggU677zzGhzv379/BGYDBNeyZcvUsWPHFv/+fffdp6ioKC1atIjgFQDghRgqNPbu3avf//736t27twYPHuz13q233qrLLrssMhNrJ4idcDohmQSEUEZGhoYMGRLQ7xhjVFpaqri4uAbvnTp1SrGxsWG5eVRVVamyslIxMTEh/ywE7tSpUz7/GwmnCy64oFW//8QTTwRpJgCA0017iqGamnc4paWlKS0tLaJzCCViJyC4WOYGRJjNZtOdd96pF154Qf369VNMTIxeeeUVT5n3O++8o5tvvllnnHGG4uPjVVZWJklasmSJhg0bpoSEBHXo0EHjx49XQUGB17lHjx6t0aNHN/jM+mW37nLpxx9/XA8//LD69OmjmJgYvffee43Oe968eRo5cqS6deumhIQEDRgwQI8//rgqKip8jl+3bp1+9KMfKS4uTj169NB//dd/qaqqymvM/PnzNWjQIHXo0EGJiYk677zzNGfOHK8xn3/+ua688kp17txZsbGxGjx4sF555ZWmLrHP7+zmq6TbGKPnn39egwcPVlxcnDp37qxrrrlG//nPf5r9HPf5vvjiC91www1KSkpSSkqKbr75Zh0/ftxrbGlpqe677z716dNHTqdTPXr0UE5Ojo4dO+Y1rnfv3powYYKWLl2qCy64QLGxsfr973+vtWvXymaz6bXXXtO9996rM888Ux06dNDEiRN14MABFRcX67bbblPXrl3VtWtXTZs2TSdOnPA6d6B/x/rzqluq7XK59PDDD+vcc89VXFycOnXqpIEDB+qZZ57x+r2vvvpKkyZNUrdu3RQTE6N+/fpp3rx5Dc5fVFSke+65x+v6zJw5UyUlJc3ODQBw+murMVRj85b8v0fWt3PnTk2bNk3nnHOO4uPj1aNHD02cOFFbt271jFm7dq2GDh0qSZo2bZpn6eBDDz0kqWFMdNVVVyk9PV0ul6vB51100UX64Q9/6HlN7ETshPaHyiQghNxPpuqy2WxyOBxex5YvX65169bpd7/7nbp3765u3bpp48aNkqSbb75Zl19+uf70pz+ppKRE0dHRevTRR/XAAw9o2rRpeuCBB1ReXq4nnnhCI0aM0KefftriEvBnn31WP/jBD/Tkk0+qY8eOOueccxod+/XXX2vSpEmem9W//vUvPfLII9q+fbv++Mc/eo3dv3+/rr/+ev32t7/V3Llz9dZbb+nhhx/W0aNH9dxzz0mSFi9erOnTp+uuu+7Sk08+Kbvdrp07d2rbtm2e83z55ZcaPny4unXrpmeffVbJycl69dVXddNNN+nAgQOaPXt2i753fb/61a+0cOFCzZgxQ//93/+tI0eOaO7cuRo+fLj+9a9/KSUlpdlz/OxnP9N1112nW265RVu3btV9990nSZ5rY4zRVVddpXfffVf33XefRowYoc8++0wPPvigNmzYoA0bNng90dy8ebP+/e9/64EHHlCfPn2UkJDgCQzmzJmjrKwsLVy4UN98843uuece3XDDDYqKitKgQYP0f//3fyooKNCcOXOUmJioZ5991nPeQP6OzXn88cf10EMP6YEHHtDIkSNVUVGh7du3ewV427Zt0/Dhw9WrVy/94Q9/UPfu3bVy5UrNmDFDhw4d0oMPPihJOnnypEaNGqXvvvtOc+bM0cCBA/XFF1/od7/7nbZu3arVq1dT3g0Ap7HTOYZqbN7+3iN92bt3r5KTk/XYY4/pjDPO0JEjR/TKK6/ooosuUkFBgc4991z98Ic/1IIFCzzf/fLLL5ekRquRbr75Zl155ZVas2aNfvzjH3uOb9++XZ9++qlXPEHsROyEdsgACLoFCxYYST5/HA6H11hJJikpyRw5csTnOaZMmeJ1vLCw0ERFRZm77rrL63hxcbHp3r27ufbaaz3HRo0aZUaNGtVgflOnTjXp6eme17t27TKSzFlnnWXKy8sD/r5VVVWmoqLCLFq0yDgcDq/vMmrUKCPJ/O1vf/P6nV/+8pfGbrebb7/91hhjzJ133mk6derU5Odcf/31JiYmxhQWFnodz87ONvHx8ebYsWNe32fBggWNfme3Bx980NT9v8INGzYYSeYPf/iD17jdu3ebuLg4M3v27Cbn6D7f448/7nV8+vTpJjY21rhcLmOMMW+//bbPcUuWLDGSzIsvvug5lp6ebhwOh/nyyy+9xr733ntGkpk4caLX8ZkzZxpJZsaMGV7Hr7rqKtOlS5dG597U39HX9UtPTzdTp071vJ4wYYIZPHhwo+c3xpjx48ebtLQ0c/z4ca/jd955p4mNjfV8Zl5enrHb7Wbjxo1e4/76178aSWbFihVNfg4AoG1qDzFUY/P29x7pK86pr7Ky0pSXl5tzzjnHzJo1y3N848aNjf5u/ZiooqLCpKSkmEmTJnmNmz17tnE6nebQoUPGGGInYie0VyxzA0Jo0aJF2rhxo9fPJ5980mDcmDFj1LlzZ5/n+NnPfub1euXKlaqsrNSUKVNUWVnp+YmNjdWoUaNatcvJFVdcoejoaL/GFhQU6IorrlBycrIcDoeio6M1ZcoUVVVVaceOHV5jExMTdcUVV3gdmzRpklwulz744ANJ0oUXXqhjx47phhtu0N/+9jcdOnSowWeuWbNGY8eOVc+ePb2O33TTTTp58qQ2bNgQyNf16R//+IdsNptuvPFGr+vbvXt3DRo0yO/rW//7Dhw4UKWlpTp48KDnu7jnXtfPf/5zJSQk6N13323w+z/4wQ98ftaECRO8Xvfr10+SPE8c6x4/cuSIV7l2IH/H5lx44YX617/+penTp2vlypUqKiryer+0tFTvvvuurr76asXHx3td35/85CcqLS3Vxx9/LKn675CRkaHBgwd7jRs/fjy7+QBAO3A6x1C+5h3IPdKXyspKPfroo+rfv7+cTqeioqLkdDr11Vdf6d///neLvlNUVJRuvPFGLV261LPcrKqqSn/605905ZVXKjk5WRKxE7ET2iuWuQEh1K9fP7+aR5555pl+v3fgwAFJ8qx5r89ub3mOuKl51FVYWKgRI0bo3HPP1TPPPKPevXsrNjZWn376qXJycnTq1Cmv8b5Km7t37y5JOnz4sCRp8uTJqqys1EsvvaSf/exncrlcGjp0qB5++GFdeumlnrG+5piamup1rtY4cOCAjDGNlmP37dvXr/O4Ayw3d9m1+9ocPnxYUVFROuOMM7zG2Ww2de/evcF3aepv06VLF6/XTqezyeOlpaXq0KFDwH/H5tx3331KSEjQq6++qhdeeEEOh0MjR47Uf//3f2vIkCE6fPiwKisr9T//8z/6n//5H5/ncCcRDxw4oJ07dzYamPtKNgIATh+nawzV2PhA7pG+5Obmat68ebr33ns1atQode7cWXa7XbfeemvA9/O6br75Zv3hD3/Q4sWL9atf/UorV67Uvn37NG3aNM8YYidiJ7RPJJMAC2hq/XL997p27SpJ+utf/6r09PQmzxsbG9ugcaHU+M3E33XUy5cvV0lJiZYuXeo1hy1btvgc7w7e6tq/f78k78Bh2rRpmjZtmkpKSvTBBx/owQcf1IQJE7Rjxw6lp6crOTlZ+/bta3CuvXv3Sqq9Nr7ExsZ6Gm/WVf9adO3aVTabTevWrfO5C0uwdrdLTk5WZWWlvv/+e6+gyBij/fv3Nwh0Q7HGPdC/Y3OioqKUm5ur3NxcHTt2TKtXr9acOXM0fvx47d69W507d5bD4dDkyZOVk5Pj8xx9+vSRVP13iIuLa7T3QFN/awBA+9HWYqjGxgdyj/Tl1Vdf1ZQpU/Too482mG+nTp0Cmltd/fv314UXXqgFCxboV7/6lRYsWKDU1FSNGzfOM4bYidgJ7RPJJKCNGT9+vKKiovT11183KN+ur3fv3nr99ddVVlbmuZEfPnxY69evV8eOHVs8B/fNuW5wYIzRSy+95HN8cXGx3nzzTa/y5ddee012u10jR45sMD4hIUHZ2dkqLy/XVVddpS+++ELp6ekaO3asli1bpr1793qqkaTqUvj4+Hj96Ec/anTOvXv31sGDB3XgwAHPk7Py8nKtXLnSa9yECRP02GOPac+ePbr22mv9uBotM3bsWD3++ON69dVXNWvWLM/xN954QyUlJRo7dmzIPtst0L9jIDp16qRrrrlGe/bs0cyZM/XNN9+of//+ysrKUkFBgQYOHOh52ufLhAkT9Oijjyo5ObnJ4BkAAH9ZIYZqTHx8vN/3SF9sNluDpM1bb72lPXv26Oyzz/Ycq1/t449p06bpjjvu0Icffqi///3vys3N9WqETuxE7IT2iWQSEEKff/55g51IJOmss85qUKLrr969e2vu3Lm6//779Z///EeXXXaZOnfurAMHDujTTz9VQkKCfv/730uqXjr2v//7v7rxxhv1y1/+UocPH9bjjz/e6iDo0ksvldPp1A033KDZs2ertLRU8+fP19GjR32OT05O1h133KHCwkL94Ac/0IoVK/TSSy/pjjvuUK9evSRJv/zlLxUXF6eLL75YZ555pvbv36+8vDwlJSV5njQ9+OCD+sc//qGsrCz97ne/U5cuXfTnP/9Zb731lh5//HElJSU1OufrrrtOv/vd73T99dfrN7/5jUpLS/Xss8+qqqrKa9zFF1+s2267TdOmTdOmTZs0cuRIJSQkaN++ffrwww81YMAA3XHHHa26fu5rOH78eN17770qKirSxRdf7NmR5IILLtDkyZNb/Rn+zCGQv2NzJk6cqIyMDA0ZMkRnnHGGvv32Wz399NNKT0/37GrzzDPP6JJLLtGIESN0xx13qHfv3iouLtbOnTv197//3dMPYebMmXrjjTc0cuRIzZo1SwMHDpTL5VJhYaHeeecd/frXv9ZFF10UtGsBALCW0zWGaoq/90hfJkyYoIULF+q8887TwIEDlZ+fryeeeKLBTm1nnXWW4uLi9Oc//1n9+vVThw4dlJqa6vWQrr4bbrhBubm5uuGGG1RWVtagZxGxE7ET2qkINv8GTltN7UQiybz00kuesZJMTk5Oo+eovyOD2/Lly01WVpbp2LGjiYmJMenp6eaaa64xq1ev9hr3yiuvmH79+pnY2FjTv39/s2TJkkZ3InniiSf8/o5///vfzaBBg0xsbKzp0aOH+c1vfmP++c9/Gknmvffe84wbNWqUOf/8883atWvNkCFDTExMjDnzzDPNnDlzTEVFhdc8s7KyTEpKinE6nSY1NdVce+215rPPPvP63K1bt5qJEyeapKQk43Q6zaBBgxrsSNLYLicrVqwwgwcPNnFxcaZv377mueeea7Bzidsf//hHc9FFF5mEhAQTFxdnzjrrLDNlyhSzadOmJq+L+3zff/+913H333PXrl2eY6dOnTL33nuvSU9PN9HR0ebMM880d9xxhzl69KjX76anp5vLL7+8wWe5dyR5/fXXfX5W/f92fM3N37+jPzuS/OEPfzDDhw83Xbt2NU6n0/Tq1cvccsst5ptvvvH6vV27dpmbb77Z9OjRw0RHR5szzjjDDB8+3Dz88MNe406cOGEeeOABc+655xqn02mSkpLMgAEDzKxZs8z+/fsbXA8AQNvXHmKoxubtPl9z90hfcc7Ro0fNLbfcYrp162bi4+PNJZdcYtatW+dzV7r/+7//M+edd56Jjo42ksyDDz5ojGm4m1tdkyZNMpLMxRdf3Oj3InYidkL7YjPGmNClqgAAAAAAAHA6afmWBQAAAAAAAGh3SCYBAAAAAADAbySTAAAAAAAA4DeSSQAAAAAAAPAbySQAAAAAAAD4jWQSAAAAAAAA/BYV6Qm0NS6XS3v37lViYqJsNlukpwMAABphjFFxcbFSU1Nlt/P8LFKInQAAaDv8jZ9IJgVo79696tmzZ6SnAQAA/LR7926lpaVFehqnjauvvlpr167V2LFj9de//rXZ8cROAAC0Pc3FTySTApSYmCip+sJ27NgxwrMBAACNKSoqUs+ePT33bgTHjBkzdPPNN+uVV17xazyxEwAAbYe/8RPJpAC5y7M7duxIQAQAQBvA0qrgysrK0tq1a/0eT+wEAEDb01z8RAMBAACAMMrLy5PNZtPMmTODet4PPvhAEydOVGpqqmw2m5YvX+5z3PPPP68+ffooNjZWmZmZWrduXVDnAQAATn8kkwAAAMJk48aNevHFFzVw4MAmx3300UeqqKhocHz79u3av3+/z98pKSnRoEGD9NxzzzV63iVLlmjmzJm6//77VVBQoBEjRig7O1uFhYWeMZmZmcrIyGjws3fvXj+/JQAAON2RTAIAAAiDEydO6Be/+IVeeuklde7cudFxLpdLOTk5mjRpkqqqqjzHd+zYoaysLC1atMjn72VnZ+vhhx/WT3/600bP/dRTT+mWW27Rrbfeqn79+unpp59Wz549NX/+fM+Y/Px8ff755w1+UlNTW/CtAQDA6YhkEgAAQBjk5OTo8ssv149//OMmx9ntdq1YsUIFBQWaMmWKXC6Xvv76a40ZM0ZXXHGFZs+e3aLPLy8vV35+vsaNG+d1fNy4cVq/fn2LzgkAANonGnADAACE2OLFi7V582Zt3LjRr/Gpqalas2aNRo4cqUmTJmnDhg0aO3asXnjhhRbP4dChQ6qqqlJKSorX8ZSUlEaXzvkyfvx4bd68WSUlJUpLS9OyZcs0dOjQBuPmzZunefPmeVVXAQCA0wPJJAAAgBDavXu37r77br3zzjuKjY31+/d69eqlRYsWadSoUerbt69efvnloOxMV/8cxpiAzrty5Uq/xuXk5CgnJ0dFRUVKSkoKaI4AAMDaWOYGAAAQQvn5+Tp48KAyMzMVFRWlqKgovf/++3r22WcVFRXVaOXOgQMHdNttt2nixIk6efKkZs2a1ap5dO3aVQ6Ho0EV0sGDBxtUKwEAADSFyiQAAIAQGjt2rLZu3ep1bNq0aTrvvPN07733yuFwNPidQ4cOaezYserXr59ef/11ffXVVxo9erRiYmL05JNPtmgeTqdTmZmZWrVqla6++mrP8VWrVunKK69s0TkBAED7RDIJAAAghBITE5WRkeF1LCEhQcnJyQ2OS9W7uV122WVKT0/XkiVLFBUVpX79+mn16tXKyspSjx49fFYpnThxQjt37vS83rVrl7Zs2aIuXbqoV69ekqTc3FxNnjxZQ4YM0bBhw/Tiiy+qsLBQt99+e5C/NQAAOJ2RTAIAALAQu92uvLw8jRgxQk6n03N8wIABWr16tZKTk33+3qZNm5SVleV5nZubK0maOnWqFi5cKEm67rrrdPjwYc2dO1f79u1TRkaGVqxYofT09NB9IQAAcNqxGWNMpCfRlribSB4/flwdO3YMyjlPlVdp8sufSJL+/MuLFBPVsNwdAAAEJhT3bAQuVH+Hf+8r0gPLP9eZSbF6btIPg3ZeAADaM3/v21QmWUCVMdr07VFJkssV4ckAAAC0ASVllcr/9qh6J8dHeioAALQ77OZmAXU34zWiUAwAAKA5tpoAisgJAIDwI5lkATZb82MAAABQV3UARcMGAADCj2SSxRAQAQAANK+2MongCQCAcCOZZAG2OgvdCIcAAACa546eeBAHAED4kUyygLrL3NhcDwAAoHk2G8vcAACIFJJJFkM8BAAA0DxaTgIAEDkkkyyABtwAAACB8fRMojQJAICwI5lkMcRDAAAAAADAykgmWUDdBtyscwMAAGieO34idAIAIPxIJlmAVwNuQiIAAIBm1S5zi+w8AABoj0gmWQAtkwAAAFqGB3EAAIQfySQLsNUpTeLpGgAAQPOoTAIAIHJIJlkM8RAAAEDz6JkEAEDkkEyygLrL3NjeFgAAoHlUJgEAEDkkkyzAxmZuAAAAAamNn4ieAAAIN5JJFlC3ZxIAAACa51nmRi4JAICwI5lkMQREAAAAzfMsc4vsNAAAaJdIJlkM29sCAAD4j36TAACEH8kki/CsdCMeAgAAaBZNAgAAiBySSRZBLgkAAMB/LHMDACBy2nUy6eqrr1bnzp11zTXXRHoqNOEGAAAICA24AQCIlHadTJoxY4YWLVoU6Wl4ISACAABonqcyieAJAICwa9fJpKysLCUmJkZ6GpLqLnMjIAIAAGgOLQIAAIicgJNJ8+fP18CBA9WxY0d17NhRw4YN0z//+c+gTuqDDz7QxIkTlZqaKpvNpuXLl/sc9/zzz6tPnz6KjY1VZmam1q1bF9R5hFPt07XIzgMAAKAtsNE0CQCAiAk4mZSWlqbHHntMmzZt0qZNmzRmzBhdeeWV+uKLL3yO/+ijj1RRUdHg+Pbt27V//36fv1NSUqJBgwbpueeea3QeS5Ys0cyZM3X//feroKBAI0aMUHZ2tgoLCz1jMjMzlZGR0eBn7969AX7r0LOxJwkAAIDfqEwCACByogL9hYkTJ3q9fuSRRzR//nx9/PHHOv/8873ec7lcysnJ0TnnnKPFixfL4XBIknbs2KGsrCzNmjVLs2fPbvAZ2dnZys7ObnIeTz31lG655RbdeuutkqSnn35aK1eu1Pz585WXlydJys/PD/TrRQ4P1wAAAPxGzyQAACKnVT2TqqqqtHjxYpWUlGjYsGENT263a8WKFSooKNCUKVPkcrn09ddfa8yYMbriiit8JpL8UV5ervz8fI0bN87r+Lhx47R+/foWnbM58+bNU//+/TV06NCQnN+NgAgAAKB57qpuIicAAMIv4MokSdq6dauGDRum0tJSdejQQcuWLVP//v19jk1NTdWaNWs0cuRITZo0SRs2bNDYsWP1wgsvtHjShw4dUlVVlVJSUryOp6SkNLp0zpfx48dr8+bNKikpUVpampYtW9ZosignJ0c5OTkqKipSUlJSi+feGE+pNhERAABAs+g3CQBA5LQomXTuuedqy5YtOnbsmN544w1NnTpV77//fqMJpV69emnRokUaNWqU+vbtq5dffrm2aWIr1D+HMSag865cubLVcwiWIFwOAACAdoedcAEACL8WLXNzOp06++yzNWTIEOXl5WnQoEF65plnGh1/4MAB3XbbbZo4caJOnjypWbNmtXjCktS1a1c5HI4GVUgHDx5sUK3UVtCAGwAAIHBUJgEAEH6t6pnkZoxRWVmZz/cOHTqksWPHql+/flq6dKnWrFmjv/zlL7rnnnta/HlOp1OZmZlatWqV1/FVq1Zp+PDhLT5vJFGqDQAA4D+qugEAiJyAl7nNmTNH2dnZ6tmzp4qLi7V48WKtXbtWb7/9doOxLpdLl112mdLT07VkyRJFRUWpX79+Wr16tbKystSjRw+fVUonTpzQzp07Pa937dqlLVu2qEuXLurVq5ckKTc3V5MnT9aQIUM0bNgwvfjiiyosLNTtt98e6FeyFEq1AQAAmudubUDkBABA+AWcTDpw4IAmT56sffv2KSkpSQMHDtTbb7+tSy+9tMFYu92uvLw8jRgxQk6n03N8wIABWr16tZKTk31+xqZNm5SVleV5nZubK0maOnWqFi5cKEm67rrrdPjwYc2dO1f79u1TRkaGVqxYofT09EC/kiXQgBsAAMB/nsIkYicAAMIu4GTSyy+/HNB4X0kmSRo8eHCjvzN69GgZP7Iq06dP1/Tp0wOaj1UFoyE5AABAe+FpEUA2CQCAsAtKzyS0nqcyKaKzAAAAaBvcm5dQ1Q0AQPiRTLIYfyqyAAAA2rvayiQAABBuJJOsgoAIAADAb7X9JomeAAAIN5JJFkEDbgAAgADwIA4AgIghmWQRNOAGAAAIHA/iAAAIP5JJFmFjf1sAAAC/2cSDOAAAIoVkksXwdA0AAKB5dYu66ZsEAEB4kUyyCE/PpIjOAgAAoG2gLgkAgMghmWQR9EwCAADwX93YicIkAADCi2SSRbCbGwAAgP/qPoYjfAIAILxIJlmEzbO9LeEQAABAc+iZBABA5JBMshhiIQAAgObV3c2N8AkAgPAimWQZ1QERySQAAAA/eFUmRW4aAAC0RySTLIL+2wAAAP7zWuZGbRIAAGFFMskiPA24CYYAAACa5dWAm/AJAICwIplkMQRDAAAAAADAykgmWQTL3AAAAPxnqxM88TAOAIDwIplkETYacAMAAPjNa5kbbQIAAAgrkkkWQWUSAACA/2zs5gYAQMRERXoCqEYDbgAAAP/Z1LafxBljVOUychnJZYxcdV+7al4bI2MaT5Y1Fjc2lVxr7C3TyC+RqAMA6+rWMUYxUY6IfDbJJIvhhg0AANA8r8qkyE1DFVUu7T5yUt8eOanCwye1v6hUx06W60hJuY6erFBxaaXKKqpUWlGlskqXSiuqVFrpUpWLoA8A0Dpv3nmxBqZ1ishnk0yyCHcTScIKAACAwDRWVRMKFVUurf/6sNbt+F4Fu4/p8z3HVVbpCsln2WyS3WaTvZkirGartJp4u7n6LloxAIB12SP4f9IkkwAAANDmhLsyaefBE/rjR7v01mf7dPxUhdd78U6HenWJV68u8UrtFKcuCU51jo9Wp3inOsZFKzbKrthoR81P9b9H2W1y2G2y2ar/6bDZZLNJDrvNk0CykckBAFgUySSLcMcK4XyyBgAA0FbVrcYJZfh0oKhUj7z1b735r72eY107xOjS/t00tHcXXdCrs3onx5P4AQC0KySTLMKTTIrsNAAAANoEr9xNiAKopZu/0+/+9oVOlFXKZpN+3C9F04b31kV9k+Vobu0ZAACnMZJJFkNhEgAAQPO8c0nBDaBcLqOH3/q3/vjRLknSoJ6d9MhVGcrokRTUzwEAoK0imWQRtaXaZJMAAAACEcyHcS6X0QN/+1yvfVIoSZox9hzdPfYcKpEAAKiDZJJFsMweAADAf3V7FAXzUdy893bqtU8KZbNJT14zSD/LTAvi2QEAOD3YIz0BVPPUJVGYBAAA0CyvZW5BCqDe235QT63eIUl65KoBJJIAAGgEySSLcD9dI5cEAADQvLpV3cGIn46frNDsNz6TMdKNP+qlSRf1CsJZAQA4PZFMshgqkwAAAJpnC3KPgLx//lvfF5epb9cEPXB5/6CeGwCA0w3JJIugZRIAAEDLtPZh3Jf7i7Vk025JUt5PByg22hGEWQEAcPoimWQVNdmkYK35BwAAON25i5NMKxe6PfPuDhkjZWd010V9k4MwMwAATm8kkyzC04A7orMAAABoOzyV3a0IoL7cX6wVW/fLZpNm/vgHwZgWAACnPZJJFkNhEgAAgH+CsYHJnz7+RpJ02fnddW73xNZPCgCAdoBkkkXUBkNkkwAAAPzhqexuYfh0oqxSyzbvkSRNHpYenEkBANAOkEyyCBpwAwAAtExLH8a9uWWvSsqr1LdrgobRKwkAAL+RTLIIG02TAAAAAuJpwN3C+Gn5luqqpOsv7OmpEgcAAM0jmWQRNrV+zT8AAEB70pr46dCJMm385ogk6ScDzgzirAAAOP2RTLIYGnADAAD4yVOZFHgAtXrbARkjDeiRpLTO8UGeGAAApzeSSRZBZTUAAEBgWtOA++0v9kuSxp+fErwJAQDQTpBMshh2cwMAAPBPSx/GnSqv0vqdhyVJ48/vHsQZAQDQPpBMsgh300eWuQEAAPjH1sL9cAsKj6q8yqXuHWN1drcOQZ4VAACnP5JJFsFmbgAAAIFp6W5uH/+nuirpor5d2MUNAIAWIJlkMS1pIAkAANAe1T6MCyx++nhX9S5uP+qbHOQZAQDQPpBMsggeigEAAASmJW0CSiuqtKXwmCSSSQAAtBTJJIvwlGlHdhoAAABtRkvaBBQUHlN5lUspHWPUOzk+FNMCAOC0RzLJIjwNJMkmAQAABCSQNgGffXdMkpSZ3pl+SQAAtBDJJIsJdM0/AABAu9WCZ3Gf7y2SJJ2fmhT8+QAA0E6QTLIIHowBAAAExrPMLYBs0hd7j0uSzk/tGPwJAQDQTpBMsoiWBEMAAADtWe0yNf8CqJKySu06VCKJyiQAAFqDZJJVtGA3EgAAgPbMs4GJn/HT9v1FMkbqlhijMxJjQjcxAABOcySTLKIlu5EAAAC0Z4HGT5/vqe6XlNGDqiQAAFqDZJLFBLIbCQAAQHsW6G5s9EsCACA4SCZZBA24AQAAAhNoz8mvDp6QJJ3bPTE0EwIAoJ0gmWQRLHMDAAAIjKdnkh8RlDFG//m+uvn2WWd0COW0AAA47ZFMsggbDbgBAAAC5H/8dLikXMdPVchmk/p0TQjxvAAAOL2162TS1Vdfrc6dO+uaa66J9FRUu8qNbBIAAEAg/EkmuauSenSKU2y0I8QzAgDg9Nauk0kzZszQokWLIj0NL1QmAQAA+CeQZW5ff1/dL6kvS9wAAGi1dp1MysrKUmKiNRow0oAbAAAgMIE04N51qLoyqS9L3AAAaLWAk0l5eXkaOnSoEhMT1a1bN1111VX68ssvgzqpDz74QBMnTlRqaqpsNpuWL1/uc9zzzz+vPn36KDY2VpmZmVq3bl1Q5xFONvea/wjPAwAAoK0I5GHct4erk0m9k+NDNBsAANqPgJNJ77//vnJycvTxxx9r1apVqqys1Lhx41RSUuJz/EcffaSKiooGx7dv3679+/f7/J2SkhINGjRIzz33XKPzWLJkiWbOnKn7779fBQUFGjFihLKzs1VYWOgZk5mZqYyMjAY/e/fuDfBbh4G7TJtsEgAAgF9sATTg3n3klCSpZxeSSQAAtFZUoL/w9ttve71esGCBunXrpvz8fI0cOdLrPZfLpZycHJ1zzjlavHixHI7qZoc7duxQVlaWZs2apdmzZzf4jOzsbGVnZzc5j6eeekq33HKLbr31VknS008/rZUrV2r+/PnKy8uTJOXn5wf69SLGU6ZNbRIAAIBfAumZ9N3Rk5JIJgEAEAyt7pl0/PhxSVKXLl0antxu14oVK1RQUKApU6bI5XLp66+/1pgxY3TFFVf4TCT5o7y8XPn5+Ro3bpzX8XHjxmn9+vUtOicAAADaFn97Jh0/VaGi0kpJUlrnuNBOCgCAdiDgyqS6jDHKzc3VJZdcooyMDJ9jUlNTtWbNGo0cOVKTJk3Shg0bNHbsWL3wwgst/txDhw6pqqpKKSkpXsdTUlIaXTrny/jx47V582aVlJQoLS1Ny5Yt09ChQ32OnTdvnubNm6eqqqoWz7spNpa5AQAABMRm86/n5O4j1VVJyQlOxTtbFf4CAAC1Mpl055136rPPPtOHH37Y5LhevXpp0aJFGjVqlPr27auXX37Zc/NvjfrnMMYEdN6VK1f6PTYnJ0c5OTkqKipSUlKS37/nLxpwAwAAhIZ7iVsaS9wAAAiKFi9zu+uuu/Tmm2/qvffeU1paWpNjDxw4oNtuu00TJ07UyZMnNWvWrJZ+rCSpa9eucjgcDaqQDh482KBaqa2orUwinQQAAOAPf+On747WNN9miRsAAEERcDLJGKM777xTS5cu1Zo1a9SnT58mxx86dEhjx45Vv379PL/zl7/8Rffcc0+LJ+10OpWZmalVq1Z5HV+1apWGDx/e4vMCAACg7ahtwN20vcdKJUmpnUgmAQAQDAEvc8vJydFrr72mv/3tb0pMTPRUByUlJSkuzvsG7XK5dNlllyk9PV1LlixRVFSU+vXrp9WrVysrK0s9evTwWaV04sQJ7dy50/N6165d2rJli7p06aJevXpJknJzczV58mQNGTJEw4YN04svvqjCwkLdfvvtgX4lSwjCqj8AAIB2qbnC7gNF1cmklI6xYZgNAACnv4CTSfPnz5ckjR492uv4ggULdNNNN3kds9vtysvL04gRI+R0Oj3HBwwYoNWrVys5OdnnZ2zatElZWVme17m5uZKkqVOnauHChZKk6667TocPH9bcuXO1b98+ZWRkaMWKFUpPTw/0K1mCp2cSq9wAAAD8Yqvdz63JcftrkkndSSYBABAUASeTAu3pc+mll/o8Pnjw4EZ/Z/To0X59zvTp0zV9+vSA5mNVtWXaZJMAAAD84e9uuPuP1ySTkmJCPCMAANqHFjfgRmhQmQQAAOAff+qSXC6jg8UscwMAIJhIJgEAAKBNstmabxNw5GS5KqqqB3RLJJkEAEAwkEyyCH+CIQAAANTyVCY1EUC5l7h17eCUM4rQFwCAYOCOahH+tY8EAACAh6fnZOPcO7l1T6IqCQCAYCGZZBG1DSRJJwEAAPijtjKp8TEHi8skscQNAIBgIplkEVQmAQAABMbdJqAph09UJ5O6dnCGejoAALQbJJMAAADQJtU+jGv8cdyhE+WSpOQOMWGYEQAA7QPJJIvwPFmjNAkAACAwTcRPh0tqkkkJVCYBABAsJJMswp8nawAAAKjlz7O42mVuVCYBABAsJJMsorYBd2TnAQAA0FbYah7HNRU/Ha5Z5kYyCQCA4CGZZBnNN5AEAABArdrKpMazSYdLqiuTkmnADQBA0JBMshgKkwAAAALTWGVSlcvoiLtnEskkAACChmSSRbDMDQAAIDDuDUwaC5+OnSyXq+bNLvEkkwAACBaSSRZBA24AAIDAeOKnRp7GuXdy6xwfrSgHYS8AAMHCXdUiqEwCAAAITHO7uR2q2cmtSwJVSQAABBPJJAAAALRJttrSbp+OnayQRDIJAIBgI5lkEZ6tbSM8DwAAgLaiNn7yHUG5k0lJcdFhmxMAAO0BySSLqH2yRjoJAADAH574qRHHT7mTSVQmAQAQTCSTLKK5Nf8AAADwrbFnccdOVTfg7hRPZRIAAMFEMskiPGXaZJMAAAD80lxh9/GaZW6dWOYGAEBQkUwCAABAm2SrKe12NZJNcvdMojIJAIDgIplkFe5lbpQmAQAA+KW5NgHuZW5J8fRMAgAgmEgmWUQzO9sCAACgHrut6TYBx1jmBgBASJBMsghbM8EQAAAAvNX2TPIdQLl3c2OZGwAAwRUV6QmgWjM727YZxhi5TJ1/ysiY6iSZCUPdVbiSceT8ACC4YqLsinbwjAuBsXt6Jvl+312ZlERlEgAAQUUyyWLClaSoqHJp95GT2nPslPYcPaV9x0t1/FSFikorVHSqQkWnKnWyolLlla7anyqXyipdqqhyeRJGxlQ3vTSiqgoA0HKP/2ygrh3aM9LTQBtT2zOpYRBSWlGlUxVVkqROcfRMAgAgmEgmWYQtxA24i0sr9OFXh/ThzkP6fM9x/Xt/scorXSH5LAAAgHBwx0++KpOKapa42WxSYiwhLwAAwcSd1SJCsczNGKMPdx7Sqx9/qzXbD6qiyjvSiot2KK1znHp0jtOZSXHqFB+tpLjanzinQzEOu5xRtT/RDrucDrscdptsturycpuqez55v27kmEUX9NmsOS0AaDei7PwfMQJX24C7YTbJ3S+pY2y07Pz3BQBAUJFMsohgN+DO//aoHn5rmwoKj3mO9emaoNHnnqEh6V2U0aOjenaOJ7gCAABtVlO7uRWVVkqiKgkAgFDg7nqaKa906YmV2/XSul2SqquPrh2SpkkXpevc7okRnh0AAEDw1C5za5hNOlHmTibRfBsAgGAjmWQRnq1tW9GC+0RZpX75yiZt+M9hSdLPM9P0m/HnqlvH2CDMEAAAwFqaquw+4a5MiiHcBQAg2Li7WoWnAXfLfv1UeZVu/P8/0Zbdx9QhJkpPXTtI487vHrz5AQAAWIz7YZyvyqTi0uqeSSxzAwAg+Li7WoS7MXVLckkul9GvX9+iLbuPqXN8tBbdfJEGpCUFd4IAAAAW42796Ct+ci9z60AyCQCAoLNHegKo1prdxF7Z8I1WbN0vp8OuF6cMIZEEAADahaZ2c6MBNwAAoUMyyWICXea2+8hJPf72l5KkByb009DeXUIwKwAAgGpXX321OnfurGuuuSbSU6nTgLvhe+6eSR1iaMANAECwkUyyiJY24P7//rFNpyqqdGGfLrrxovTgTwwAAKCOGTNmaNGiRZGehqSmG3DTMwkAgNAhmWQRthY04P5i73G9s+2AbDbp0aszZLe3Yq0cAACAH7KyspSYmBjpaUhqugG3u2cSySQAAIKPZJJF2BR4Iui5NTslSRMGpursbtYI6gAAQEPz58/XwIED1bFjR3Xs2FHDhg3TP//5z6B+xgcffKCJEycqNTVVNptNy5cv9znu+eefV58+fRQbG6vMzEytW7cuqPMIJ0/PJB/veRpwx5BMAgAg2EgmWUSgDbh3Hzmpf36+X5J015izQzAjAAAQLGlpaXrssce0adMmbdq0SWPGjNGVV16pL774wuf4jz76SBUVFQ2Ob9++Xfv37/f5OyUlJRo0aJCee+65RuexZMkSzZw5U/fff78KCgo0YsQIZWdnq7Cw0DMmMzNTGRkZDX727t0b4LcOPXtNJNt0A256JgEAEGw8qrEYX8GQL3/ZtFuSdMnZXfWDFKqSAACwsokTJ3q9fuSRRzR//nx9/PHHOv/8873ec7lcysnJ0TnnnKPFixfL4XBIknbs2KGsrCzNmjVLs2fPbvAZ2dnZys7ObnIeTz31lG655RbdeuutkqSnn35aK1eu1Pz585WXlydJys/Pb/H3DDd3ZbfLRwfuEzU9k6hMAgAg+KhMsohAeia5XEZv5H8nSbp2aM8QzgoAAARbVVWVFi9erJKSEg0bNqzB+3a7XStWrFBBQYGmTJkil8ulr7/+WmPGjNEVV1zhM5Hkj/LycuXn52vcuHFex8eNG6f169e36JxNmTdvnvr376+hQ4cG/dxunvjJx3vFpfRMAgAgVLi7WoStiTX/9X2257j2Hi9VgtOhcf1TQjsxAAAQFFu3btWwYcNUWlqqDh06aNmyZerfv7/PsampqVqzZo1GjhypSZMmacOGDRo7dqxeeOGFFn/+oUOHVFVVpZQU79ghJSWl0aVzvowfP16bN29WSUmJ0tLStGzZMp8Jo5ycHOXk5KioqEhJSUktnndT3PGTj8IkGnADABBC3F0toqndSOpb+UV1wDf63G6KjXaEcFYAACBYzj33XG3ZskXHjh3TG2+8oalTp+r9999vNKHUq1cvLVq0SKNGjVLfvn318ssve5InrVH/HMaYgM67cuXKVs8hWOyeym7v+KnKZXSyvEoSy9wAAAgFlrlZRCDL3N799wFJ0rjzqUoCAKCtcDqdOvvsszVkyBDl5eVp0KBBeuaZZxodf+DAAd12222aOHGiTp48qVmzZrXq87t27SqHw9GgCungwYMNqpXaCs9ubvXiJ3dVkiR1oDIJAICgI5lkEU1tbVvX98Vl2nHghGw2aeQ5Z4R+YgAAICSMMSorK/P53qFDhzR27Fj169dPS5cu1Zo1a/SXv/xF99xzT4s/z+l0KjMzU6tWrfI6vmrVKg0fPrzF542kxiq7T5ZXJ5Oi7DbFRFHFDQBAsPGoxiLcwVBzu7lt+M9hSVK/7h3VOcEZ4lkBAIBgmDNnjrKzs9WzZ08VFxdr8eLFWrt2rd5+++0GY10uly677DKlp6dryZIlioqKUr9+/bR69WplZWWpR48ePquUTpw4oZ07d3pe79q1S1u2bFGXLl3Uq1cvSVJubq4mT56sIUOGaNiwYXrxxRdVWFio22+/PXRfPoQa6znpXuIW7ySRBABAKJBMsghbI2Xa9W34ujqZNOys5FBPCQAABMmBAwc0efJk7du3T0lJSRo4cKDefvttXXrppQ3G2u125eXlacSIEXI6ax8cDRgwQKtXr1Zysu8YYNOmTcrKyvK8zs3NlSRNnTpVCxculCRdd911Onz4sObOnat9+/YpIyNDK1asUHp6ehC/bfi42wQ0qEwqcyeTCHUBAAgF7rAWUbu1bdPZpILCo5KkC/t0CfWUAABAkLz88ssBjfeVZJKkwYMHN/o7o0ePbrbCWZKmT5+u6dOnBzQfq7I30nPSvcwtPobKJAAAQoGeSRZhU+Nb27qdKq/SjgPFkqTBPTuFYVYAAADWVduAu37PJJa5AQAQSiSTLMKf3dy27Tsul5G6JcYopWNseCYGAABgUbXL3LyP1yaTKMIHACAUSCZZhN2PZW6ffXdckjQwLSkcUwIAALC0xnpOlriXuVGZBABASJBMsgh/GnB/vqdIkpTRg2QSAACAezfchg24q5NJCVQmAQAQEiSTLMIdDDXVOHPnwep+SeemJIZhRgAAANbm6ZlU7/jJiuplbnFUJgEAEBIkkyyiucokY4y+/r5EknR2tw7hmhYAAIBl1e7mVr8yqTqZlEAyCQCAkCCZZBGNNZB0O1BUphNllXLYbUpPTgjfxAAAACzK/TCuwTK3cndlEsvcAAAIBZJJFuFZ5tZIA+6dB09IktK7xMsZxZ8NAACgsd1wT5a7eyZRmQQAQCiQlbAIezPL3Nz9ks5iiRsAAIAkySZ3ZZL3cXdlUnwMlUkAAIQCySSLsDWy5t/tm8MnJUl9z2CJGwAAgFSnZ5LqL3OrrkyKpzIJAICQIJlkEbXL3Hz77ugpSVLPzvFhmQ8AAIDV2e2+K7s9lUkkkwAACAmSSRbR3G5u3x2trkzq0TkuXFMCAACwNPfDOFe9dW4lnmQSy9wAAAgFkkkWUbubm+9s0p5j7sokkkkAAABSnYdx9Y6fogE3AAAhRTLJItwNJH2lko6fqlBxaXVQ1KMTy9wAAACkxh/GlZRVVybFkUwCACAkSCZZhL2RrW2l2iVuyQlOgiIAAIAajcVPpyqqk0kJ7OYGAEBIkEyyiKZ2c9tT03w7jSVuAAAAHnZPz0nfu7nFRfMQDgCAUCCZZBFNNeB290tK7UQyCQAAwM3TgLtO/GSMUWmFS5IUSzIJAICQIJlkEZ7KJB9dk74vLpMkpXSMDeeUAAAALK22AXdt/FRW6fL8e2w0oS4AAKHAHdYi3A24XT4qk9zJpDMSY8I5JQAAAEurbcBde6y0pl+SRGUSAAChQjLJImxNNOA+SDIJAACgAV89k9zNt6PsNkU7CHUBAAgF7rAWYfdjmRvJJAAAgFq+dnNz90ui+TYAAKFDMski3MvcfFUmfX+iJpnUgWQSAACAm7tnkqtOAOVe5hZDMgkAgJAhmWQRtcvcvLNJVS6jwzXJpG4dSSYBAAC4+WoT4F7mRvNtAABCh7usRdTuRuLtSEm5XKa6jDs5gWQSAACAm93WcAMTd2USy9wAAAgdkkkWUfNgrcFubu5+SV0SYuRwNwYAAACAJ36qW9ldVtMziZ3cAAAIHZJJFtHYMjdPvySabwMAAHix+6jsZpkbAAChx13WInwFQ5J0pKQ6mZSc4AzzjAAAAKzN/TDOVwNuKpMAAAgdkkkW0Vhl0tGSCklSp/jocE8JAADA0jw9J716JrHMDQCAUCOZZBG1a/69jx87VZ1M6hxPZRIAAEBddh+VSaeoTAIAIORIJlmErydrknT8ZLkkKpMAAADq8/UwrnY3N8JcAABChbusRfha8y9JR09WVyYlxZFMAgAAqMtud/ecpGcSAADhRDLJIhprwM0yNwAAAN/clUkuV+2x2sokkkkAAIQKySSLqC3T9k4nHWOZGwAAgE+eNgFelUnVmaUYkkkAAIQMySSLqN3Nzfv4sZPs5gYAAOCLu7LbVSd+qm3ATZgLAECocJe1CFtjy9w8lUkscwMAAKir9mFcw55JLHMDACB0SCZZhGfNf51gqLLKpaLSSklSJxpwAwAAeLH7qOx2L3OjATcAAKFDMskiPJVJdYIhdyJJYjc3AACA+mxyL3PztZsbYS4AAKHCXdYiPE/W6hxzL3FLjI1SlIM/FQAAQF02H/ETy9wAAAg9MhQW4WvN/7FT1c23qUoCAABoyFcD7tLK6mQSu7kBABA6JJMswl2mXXeZ24maZW6JsSSTAAAA6vP1MO5Uec0ytyiSSQAAhArJJIuoLdOuDYZOlNUkk2KiIjElAAAAS7P76DlZVuluwE2YCwBAqHCXtQh3A26Xq/aYuzKpQyzJJAAAgPrcD+PqNuAur0kmxVCZBABAyJBMsoiaWMirMqm4pjKpA5VJAAAADfjaDdddmeSMIswFACBUuMtahK8ybSqTAAAAGufeDbeqTgBV5m7ATTIJAICQ4S5rEbUNJGuPnSir3s2NnkkAAAAN2WubTlb/w5g6y9wIcwEACBXushbha5nbCZa5AQAANKp+ZVKly8hVE0rRMwkAgNAhmWQRvtb8F7PMDQAAoFHuyiR3A253vyRJimE3NwAAQoa7rEX42o2EyiQAAIDGeZJJNeVI5XWSSU4HYS4AAKHCXdYiape51XI34E6kMgkAAKABh91dmVT92t18O9phk929Bg4AAAQdySSL8Lmbm6cyKToSUwIAALC0+pXdtc236ZcEAEAokUyyiNrd3GqzSfRMAgAAaJy7MqnK5d0zyclObgAAhBR3Wouot7OtJHomAQAANKV+ZXdZhbsyiRAXAIBQ4k5rEfV3czPGeJJJ9EwCAABoyP0wrsq9zK2qumcSlUkAAIQWd1qLcLeIdK/5L61weUq2qUwCAABoyOHeza0mfqIyCQCA8OBOaxH1K5OKyypqjkvxTppIAgAA1Ofesc1Vr2cSDbgBAAgtkkkWYa/XM+lUeXWZdny0w5NoAgAAQC27pzKp+jUNuAEACA/utBZhk7syqToaOlVRnUyKoyoJAADAJ/fDuNrd3KrjJ5a5AQAQWtxpLcKzm1vNkzV3ZVJsNMkkAAAAXxx274dx5ZX0TAIAIBy401qEJ5mkepVJJJMAAAB8ci9zc+/mxjI3AADCgzutRbiXubnX/JeyzA0AAKBJjfVMogE3AAChRTLJImqXudVUJpVXB0MscwMAAPDNXhPJundzK6cyCQCAsGiXd9qrr75anTt31jXXXBPpqXi4n6y5d3NzVyaRTAIAAPDN4alMogE3AADh1C7vtDNmzNCiRYsiPQ0vDRpwe3omtcs/EQAAQLNsLHMDACAi2mWmIisrS4mJiZGehpeaXJJnmVspDbgBAACa5N7NjWVuAACEl+XutB988IEmTpyo1NRU2Ww2LV++vMGY559/Xn369FFsbKwyMzO1bt268E80yGz1lrmdKqcBNwAAQFNqckkscwMAIMwsd6ctKSnRoEGD9Nxzz/l8f8mSJZo5c6buv/9+FRQUaMSIEcrOzlZhYaFnTGZmpjIyMhr87N27N1xfI2C2esHQKXomAQAANMndc7LKeFcmxdAmAACAkIqK9ATqy87OVnZ2dqPvP/XUU7rlllt06623SpKefvpprVy5UvPnz1deXp4kKT8/P2jzKSsrU1lZmed1UVFR0M5dV+0yt+p/nmKZGwAAQJPsdt89k5wOkkkAAIRSm7rTlpeXKz8/X+PGjfM6Pm7cOK1fvz4kn5mXl6ekpCTPT8+ePUPyOZ7d3GqCodKK6mCIyiQAAADfPLu51WSTyirclUnETwAAhFKbSiYdOnRIVVVVSklJ8TqekpKi/fv3+32e8ePH6+c//7lWrFihtLQ0bdy4sdGx9913n44fP+752b17d4vn35Ta3dxowA0AAOCP+j2TyqtqkklUJgEAEFKWW+bmD3ezajdjTINjTVm5cqXfY2NiYhQTE+P3+JayyXcD7lgacAMAAPhUd5mbMaa2ATc9kwAACKk2daft2rWrHA5HgyqkgwcPNqhWamtqK5Oq/0nPJAAAgKbZ6zxMNKbOMjd2cwMAIKTa1J3W6XQqMzNTq1at8jq+atUqDR8+PEKzCo7GdnMjmQQAAOCbo04yqcoYzzI3J8kkAABCynLL3E6cOKGdO3d6Xu/atUtbtmxRly5d1KtXL+Xm5mry5MkaMmSIhg0bphdffFGFhYW6/fbbIzjr1qu/zM3TM8lJMAQAAOCLrU6Y5DKmTmUSD+MAAAglyyWTNm3apKysLM/r3NxcSdLUqVO1cOFCXXfddTp8+LDmzp2rffv2KSMjQytWrFB6enqkphwU9ppgqHY3t5qeSQRDAAAAPjnqLXPzNOCmMgkAgJCyXDJp9OjRnh3NGjN9+nRNnz49TDMKD09lUr1lbjTgBgAA8K1uz6Qql1FZTfzEMjcAAEKLO61FeBpw17w+VV79ZI2eSQAAAL7Z6y9zq2SZGwAA4UAyySLsnt3cqtNJpTTgBgAAaFLdyiSXSyqvpAE3AADhwJ3WMqqDIZepTih5dnNjmRsAAIBPdXsmeVcmEeICABBK3Gktwl2Z5DJGFVVGVa7qCqVYKpMAAAB8qpNLUqXL0IAbAIAw4U5rEY6abJLLZVRaWeU5HhvNnwgAAMAXm83mSSi5WwRILHMDACDUuNNahHvNf5UxKi2vqjkmOR38iQAAABrjXupWN5lEA24AAEKLTIVF1FYmqbZfUrRDtrr12wAAAPDifiDnjp9sNinaQfwEAEAokUyyCHcyqYrm2wAAAH6z10SzJ2squ50OOw/jAAAIMZJJFuFZ5uYyOlUTDNF8GwAAoGn1K5Novg0AQOhxt7UId2WSJE8yKY5kEgAAQJM8PZPclUn0SwIAIORIJlmEo0459omySklUJgEAADTHHUJRmQQAQPhwt/XTvHnz1L9/fw0dOjQk57fX+UuUlFcnk6hMAgAAaJq7utuTTIomvAUAINS42/opJydH27Zt08aNG0Ny/rrL3E6U1fRMogE3AABAkzw9k8rdlUnETwAAhBrJJIuw11nmVlLmrkzizwMAANAUe80DuVKWuQEAEDbcbS2ibmVSbTKJJ2sAAABNcYdQJ8tJJgEAEC7cbS3C4VWZVLObG8vcAAAAmuSOoU56dnMjvAUAINS421qE3UdlEmv+AQAAmhblqA5nT5YTPwEAEC4kkyzEvdTthHuZG5VJAAAATYqqiZ9KytnNDQCAcOFuayENkkn0TAIAAGiSO3466ansJrwFACDUuNtaiHvNPw24AQAA/ONe5uapTGKZGwAAIUcyyULqVybFsswNAACgSZ5lblQmAQAQNtxtLaT+1raxBEMAAABNinJ47+ZGMgkAgNDjbmshjnpP1mjADQAA0LSoepXdJJMAAAg97rYWQgNuAACAwDRowE38BABAyJFMshB7TQPuskqXJJJJAAAAzYmuacB9soJlbgAAhAt3WwtxP1lzowE3AABA09zxkzHVr50kkwAACDnuthbirkxyozIJAACgaVF273CWyiQAAEKPu62FNKhMIpkEAADQpKh68VNMFPETAAChRjLJQuonk6hMAgAAaJrDUT+ZRHgLAECocbe1kHq5JJJJAAAAzYiuX5kUTXgLAECocbe1kIYNuPnzAAAANMVRr2eS08HDOAAAQo1shYXUbcBtt0lOB38eAACApjTomURlEgAAIcfd1k/z5s1T//79NXTo0JB9Rt3KpLhoh2z1dncDAACAtyh6JgEAEHbcbf2Uk5Ojbdu2aePGjSH7jLrJJHZyAwAAaB67uQEAEH4kkyyk7jI3kkkAAADNi6rXFoDKJAAAQo+7rYV4LXNzkkwCAABoTv3KJCfJJAAAQo67rYU4bN49kwAAANC0+rvhUpkEAEDocbe1kPoNuAEAANC0BsvciKEAAAg5kkkWUjeZxLa2AAAAzWvYgJsYCgCAUONuayH2OsFQh5ioCM4EAACgbYhy1MZPdlvD5BIAAAg+kkkWUicWUryTZBIAAEBz6iaPEpxRstlIJgEAEGokkyyk7u4j8ezmBgAA0Kwoe238xG64AACEB8kkC6lbjRQfQzAEAADQnLoP4BJoEwAAQFiQTLKQuk/TEljmBgAA0Kz4OgkkKrsBAAgPkkkWEl9nK1uCIQAAgOYl8DAOAICwI5lkIXWfrFGmDQAA0Ly6ld20CQAAIDxIJllI3WokKpMAAACaV7caifgJAIDwIJlkIfGUaQMAAAQkIabuwzjiJwAAwoFkkoXE1emZ1DnBGcGZAAAAtA1xdRJICVQmAQAQFiSTLCTKYfP8e3pyfARnAgAA0DbUTSAlxUVHcCYAALQfJJMs5PCJcs+/J1OZBAAA0Ky6S9vO6tYhgjMBAKD9IJlkIT/qmyxJ6pLglM1ma2Y0AAAAnFF22WvCpuFndY3sZAAAaCfoUmghGT2S9NaMS5SaFBfpqQAAALQZ6+4do7KKKp2RGBPpqQAA0C6QTLKY81OTIj0FAACANqVHJx7EAQAQTixzAwAAAAAAgN9IJvlp3rx56t+/v4YOHRrpqQAAAAAAAEQMySQ/5eTkaNu2bdq4cWOkpwIAAAAAABAxJJMAAAAAAADgN5JJAAAAAAAA8BvJJAAAAAAAAPiNZBIAAAAAAAD8RjIJAAAAAAAAfiOZBAAAAAAAAL+RTAIAAAAAAIDfoiI9gbbGGCNJKioqivBMAABAU9z3ave9G5FB7AQAQNvhb/xEMilAxcXFkqSePXtGeCYAAMAfxcXFSkpKivQ02i1iJwAA2p7m4ieb4XFdQFwul/bu3avExETZbLagnbeoqEg9e/bU7t271bFjx6CdF964zqHHNQ4PrnN4cJ1DL5TX2Bij4uJipaamym5nZX+khCp2kvjfaDhwjcOD6xweXOfQ4xqHhxXiJyqTAmS325WWlhay83fs2JH/0YUB1zn0uMbhwXUOD65z6IXqGlORFHmhjp0k/jcaDlzj8OA6hwfXOfS4xuERyfiJx3QAAAAAAADwG8kkAAAAAAAA+I1kkkXExMTowQcfVExMTKSnclrjOoce1zg8uM7hwXUOPa4xWoP/fkKPaxweXOfw4DqHHtc4PKxwnWnADQAAAAAAAL9RmQQAAAAAAAC/kUwCAAAAAACA30gmAQAAAAAAwG8kkwAAAAAAAOA3kkkW8Pzzz6tPnz6KjY1VZmam1q1bF+kptRl5eXkaOnSoEhMT1a1bN1111VX68ssvvcYYY/TQQw8pNTVVcXFxGj16tL744guvMWVlZbrrrrvUtWtXJSQk6IorrtB3330Xzq/SpuTl5clms2nmzJmeY1zn4NizZ49uvPFGJScnKz4+XoMHD1Z+fr7nfa5z61RWVuqBBx5Qnz59FBcXp759+2ru3LlyuVyeMVzjwH3wwQeaOHGiUlNTZbPZtHz5cq/3g3VNjx49qsmTJyspKUlJSUmaPHmyjh07FuJvB6sifmo54qfIIH4KHeKn0CJ+Co02Hz8ZRNTixYtNdHS0eemll8y2bdvM3XffbRISEsy3334b6am1CePHjzcLFiwwn3/+udmyZYu5/PLLTa9evcyJEyc8Yx577DGTmJho3njjDbN161Zz3XXXmTPPPNMUFRV5xtx+++2mR48eZtWqVWbz5s0mKyvLDBo0yFRWVkbia1nap59+anr37m0GDhxo7r77bs9xrnPrHTlyxKSnp5ubbrrJfPLJJ2bXrl1m9erVZufOnZ4xXOfWefjhh01ycrL5xz/+YXbt2mVef/1106FDB/P00097xnCNA7dixQpz//33mzfeeMNIMsuWLfN6P1jX9LLLLjMZGRlm/fr1Zv369SYjI8NMmDAhXF8TFkL81DrET+FH/BQ6xE+hR/wUGm09fiKZFGEXXnihuf32272OnXfeeea3v/1thGbUth08eNBIMu+//74xxhiXy2W6d+9uHnvsMc+Y0tJSk5SUZF544QVjjDHHjh0z0dHRZvHixZ4xe/bsMXa73bz99tvh/QIWV1xcbM455xyzatUqM2rUKE8wxHUOjnvvvddccskljb7PdW69yy+/3Nx8881ex37605+aG2+80RjDNQ6G+sFQsK7ptm3bjCTz8ccfe8Zs2LDBSDLbt28P8beC1RA/BRfxU2gRP4UW8VPoET+FXluMn1jmFkHl5eXKz8/XuHHjvI6PGzdO69evj9Cs2rbjx49Lkrp06SJJ2rVrl/bv3+91jWNiYjRq1CjPNc7Pz1dFRYXXmNTUVGVkZPB3qCcnJ0eXX365fvzjH3sd5zoHx5tvvqkhQ4bo5z//ubp166YLLrhAL730kud9rnPrXXLJJXr33Xe1Y8cOSdK//vUvffjhh/rJT34iiWscCsG6phs2bFBSUpIuuugiz5gf/ehHSkpK4rq3M8RPwUf8FFrET6FF/BR6xE/h1xbip6hW/TZa5dChQ6qqqlJKSorX8ZSUFO3fvz9Cs2q7jDHKzc3VJZdcooyMDEnyXEdf1/jbb7/1jHE6nercuXODMfwdai1evFibN2/Wxo0bG7zHdQ6O//znP5o/f75yc3M1Z84cffrpp5oxY4ZiYmI0ZcoUrnMQ3HvvvTp+/LjOO+88ORwOVVVV6ZFHHtENN9wgif+WQyFY13T//v3q1q1bg/N369aN697OED8FF/FTaBE/hR7xU+gRP4VfW4ifSCZZgM1m83ptjGlwDM2788479dlnn+nDDz9s8F5LrjF/h1q7d+/W3XffrXfeeUexsbGNjuM6t47L5dKQIUP06KOPSpIuuOACffHFF5o/f76mTJniGcd1brklS5bo1Vdf1Wuvvabzzz9fW7Zs0cyZM5WamqqpU6d6xnGNgy8Y19TXeK57+0X8FBzET6FD/BQexE+hR/wUOVaOn1jmFkFdu3aVw+FokBE8ePBggwwkmnbXXXfpzTff1Hvvvae0tDTP8e7du0tSk9e4e/fuKi8v19GjRxsd097l5+fr4MGDyszMVFRUlKKiovT+++/r2WefVVRUlOc6cZ1b58wzz1T//v29jvXr10+FhYWS+O85GH7zm9/ot7/9ra6//noNGDBAkydP1qxZs5SXlyeJaxwKwbqm3bt314EDBxqc//vvv+e6tzPET8FD/BRaxE/hQfwUesRP4dcW4ieSSRHkdDqVmZmpVatWeR1ftWqVhg8fHqFZtS3GGN15551aunSp1qxZoz59+ni936dPH3Xv3t3rGpeXl+v999/3XOPMzExFR0d7jdm3b58+//xz/g41xo4dq61bt2rLli2enyFDhugXv/iFtmzZor59+3Kdg+Diiy9usDXzjh07lJ6eLon/noPh5MmTstu9b30Oh8OztS3XOPiCdU2HDRum48eP69NPP/WM+eSTT3T8+HGueztD/NR6xE/hQfwUHsRPoUf8FH5tIn5qVftutJp7a9uXX37ZbNu2zcycOdMkJCSYb775JtJTaxPuuOMOk5SUZNauXWv27dvn+Tl58qRnzGOPPWaSkpLM0qVLzdatW80NN9zgc0vFtLQ0s3r1arN582YzZsyYdr1NpT/q7kZiDNc5GD799FMTFRVlHnnkEfPVV1+ZP//5zyY+Pt68+uqrnjFc59aZOnWq6dGjh2dr26VLl5quXbua2bNne8ZwjQNXXFxsCgoKTEFBgZFknnrqKVNQUODZpj1Y1/Syyy4zAwcONBs2bDAbNmwwAwYMCMrWtmh7iJ9ah/gpcoifgo/4KfSIn0KjrcdPJJMsYN68eSY9Pd04nU7zwx/+0LMtK5onyefPggULPGNcLpd58MEHTffu3U1MTIwZOXKk2bp1q9d5Tp06Ze68807TpUsXExcXZyZMmGAKCwvD/G3alvrBENc5OP7+97+bjIwMExMTY8477zzz4osver3PdW6doqIic/fdd5tevXqZ2NhY07dvX3P//febsrIyzxiuceDee+89n/9fPHXqVGNM8K7p4cOHzS9+8QuTmJhoEhMTzS9+8Qtz9OjRMH1LWA3xU8sRP0UO8VNoED+FFvFTaLT1+MlmjDGtq20CAAAAAABAe0HPJAAAAAAAAPiNZBIAAAAAAAD8RjIJAAAAAAAAfiOZBAAAAAAAAL+RTAIAAAAAAIDfSCYBAAAAAADAbySTAAAAAAAA4DeSSQAAAAAAAPAbySQAAAAAAAD4jWQSAAAAAAAA/EYyCQAAAAAAAH4jmQQAAAAAAAC//T+xocceYJNgTAAAAABJRU5ErkJggg==", + "text/plain": [ + "<Figure size 1400x400 with 2 Axes>" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "plot_convergences(r1, np.ones(4))" + ] }, { "cell_type": "code", - "execution_count": null, + "execution_count": 52, "metadata": {}, - "outputs": [], - "source": [] + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Itération du minimum : 0 0\n" + ] + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAABIAAAAF0CAYAAABbk0LyAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjMsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvZiW1igAAAAlwSFlzAAAPYQAAD2EBqD+naQAAcn1JREFUeJzt3XlcVXX+x/H3ZRdEFFEUEURzQ1QSrNwlC0NFpZqxmtRcpkzLlBobtdVUmqYca1wamxmtaZr8NaWWQxFWiqWNuRVormG4IYoLi7Lde35/lHci0EQv97C8no8Hv9+c7/2ecz7ne4375XO/i8UwDEMAAAAAAACos1zMDgAAAAAAAADViwQQAAAAAABAHUcCCAAAAAAAoI4jAQQAAAAAAFDHkQACAAAAAACo40gAAQAAAAAA1HEkgAAAAAAAAOo4EkAAAAAAAAB1HAkgAECVHD16VM2bN9ejjz5qdigAAAA1Hn0n1BQkgIAfrVixQhaL5ZI/69evNztEpxk4cKAiIiKces9Dhw7JYrFoxYoVTr1vXXbfffepTZs25cratGmj++6776qvWVZWprvuuks333yzXnzxxWsLEABQJ9CH+mXX0s85duyYnnnmGe3cubPCa88884wsFsu1BwhJ9J1Q97mZHQBQ0yxfvlydOnWqUB4eHm5CNIBjrVq1So0aNbrq82fOnCk3Nze98cYbdDgBAOXQh6oex44d07PPPqs2bdooMjKy3GsTJ07UbbfdZk5g9QR9J9QlJICAn4mIiFB0dHSVzjEMQ0VFRWrQoEGF1y5cuCAvLy+n/MK3Wq0qKyuTp6dntd8LVXfhwoVK/4040/XXX39N5//xj390UCQAgLqmPvWhLhe3MwUHBys4ONjUGKoTfSfAsZgCBlwFi8Wihx56SK+++qo6d+4sT09Pvf766/Yh0B9//LHGjx+vZs2aydvbW8XFxZKklStXqlevXvLx8VHDhg01ePBg7dixo9y1Bw4cqIEDB1a458+HpF4cSvzCCy9o7ty5CgsLk6enpz777LNLxr148WL1799fzZs3l4+Pj7p27aoXXnhBpaWlldbfuHGjbrrpJjVo0ECtWrXSk08+KavVWq7O0qVL1b17dzVs2FC+vr7q1KmTZs2aVa5ORkaGRowYoSZNmsjLy0uRkZF6/fXXL9fElT7zRZUNdzYMQ0uWLFFkZKQaNGigJk2a6M4779R33333i/e5eL1du3bp7rvvlp+fnwIDAzV+/HidO3euXN2ioiLNnDlTYWFh8vDwUKtWrTRlyhSdPXu2XL02bdpo2LBheu+993T99dfLy8tLzz77rNavXy+LxaK33npLjz/+uFq2bKmGDRsqPj5eJ06cUH5+vu6//34FBAQoICBA48aNU0FBQblrV/V9/HlcPx3GbLPZNHfuXHXs2FENGjRQ48aN1a1bN7388svlztu/f7/uueceNW/eXJ6enurcubMWL15c4fp5eXl67LHHyrXPtGnTVFhY+IuxAQDqvtrah7pU3NKVf0b+3IEDBzRu3Di1b99e3t7eatWqleLj45Wenm6vs379evXs2VOSNG7cOPu0umeeeUZSxT7RyJEjFRoaKpvNVuF+N954o3r06GE/pu9E3wn1DyOAgJ+5+A3QT1ksFrm6upYrW716tTZu3KinnnpKLVq0UPPmzfXVV19JksaPH6+hQ4fqH//4hwoLC+Xu7q758+friSee0Lhx4/TEE0+opKREf/zjH9WvXz9t2bLlqodHv/LKK+rQoYNefPFFNWrUSO3bt79k3YMHD+qee+6xf8B8/fXXmjdvnvbs2aO///3v5epmZ2frrrvu0u9//3vNmTNH//nPfzR37lydOXNGixYtkiS9/fbbmjx5sh5++GG9+OKLcnFx0YEDB7R79277dfbu3avevXurefPmeuWVV9S0aVO9+eabuu+++3TixAnNmDHjqp775x544AGtWLFCU6dO1R/+8AedPn1ac+bMUe/evfX1118rMDDwF69xxx13aNSoUZowYYLS09M1c+ZMSbK3jWEYGjlypD755BPNnDlT/fr10zfffKOnn35amzdv1ubNm8t9c7h9+3Z9++23euKJJxQWFiYfHx/7h/msWbMUExOjFStW6NChQ3rsscd09913y83NTd27d9e//vUv7dixQ7NmzZKvr69eeeUV+3Wr8j7+khdeeEHPPPOMnnjiCfXv31+lpaXas2dPuU7Z7t271bt3b4WEhOill15SixYtlJKSoqlTp+rUqVN6+umnJUnnz5/XgAEDdOTIEc2aNUvdunXTrl279NRTTyk9PV3r1q1j6DMA1GF1uQ91qbiv9DOyMseOHVPTpk31/PPPq1mzZjp9+rRef/113XjjjdqxY4c6duyoHj16aPny5fZnHzp0qCRdctTP+PHjNWLECH366ae65ZZb7OV79uzRli1byvUn6DvRd0I9ZAAwDMMwli9fbkiq9MfV1bVcXUmGn5+fcfr06UqvMWbMmHLlWVlZhpubm/Hwww+XK8/PzzdatGhh/PrXv7aXDRgwwBgwYECF+MaOHWuEhobajzMzMw1JRrt27YySkpIqP6/VajVKS0uNN954w3B1dS33LAMGDDAkGWvWrCl3zm9/+1vDxcXF+P777w3DMIyHHnrIaNy48WXvc9dddxmenp5GVlZWufK4uDjD29vbOHv2bLnnWb58+SWf+aKnn37a+Omvr82bNxuSjJdeeqlcvcOHDxsNGjQwZsyYcdkYL17vhRdeKFc+efJkw8vLy7DZbIZhGMZHH31Uab2VK1cakoxly5bZy0JDQw1XV1dj79695ep+9tlnhiQjPj6+XPm0adMMScbUqVPLlY8cOdLw9/e/ZOyXex8ra7/Q0FBj7Nix9uNhw4YZkZGRl7y+YRjG4MGDjeDgYOPcuXPlyh966CHDy8vLfs+kpCTDxcXF+Oqrr8rV+/e//21IMpKTky97HwBA7VQf+lCXivtKPyMr6+f8XFlZmVFSUmK0b9/emD59ur38q6++uuS5P+8TlZaWGoGBgcY999xTrt6MGTMMDw8P49SpU4Zh0Hei74T6iilgwM+88cYb+uqrr8r9/Pe//61Q7+abb1aTJk0qvcYdd9xR7jglJUVlZWUaM2aMysrK7D9eXl4aMGDANe2OMXz4cLm7u19R3R07dmj48OFq2rSpXF1d5e7urjFjxshqtWrfvn3l6vr6+mr48OHlyu655x7ZbDalpaVJkm644QadPXtWd999t9asWaNTp05VuOenn36qQYMGqXXr1uXK77vvPp0/f16bN2+uyuNWau3atbJYLLr33nvLtW+LFi3UvXv3K27fnz9vt27dVFRUpJycHPuzXIz9p371q1/Jx8dHn3zySYXzO3ToUOm9hg0bVu64c+fOkmT/Zu+n5adPny43lLkq7+MvueGGG/T1119r8uTJSklJUV5eXrnXi4qK9MknnyghIUHe3t7l2nfIkCEqKirSl19+KemH9yEiIkKRkZHl6g0ePJhdYACgHqjLfajK4q7KZ2RlysrKNH/+fIWHh8vDw0Nubm7y8PDQ/v379e23317VM7m5uenee+/Ve++9Z5+KZbVa9Y9//EMjRoxQ06ZNJdF3ou+E+oopYMDPdO7c+YoWMGzZsuUVv3bixAlJss/h/jkXl6vPxV4ujp/KyspSv3791LFjR7388stq06aNvLy8tGXLFk2ZMkUXLlwoV7+yYb8tWrSQJOXm5kqSRo8erbKyMr322mu64447ZLPZ1LNnT82dO1e33nqrvW5lMQYFBZW71rU4ceKEDMO45FDltm3bXtF1LnaKLro4JPli2+Tm5srNzU3NmjUrV89isahFixYVnuVy742/v3+5Yw8Pj8uWFxUVqWHDhlV+H3/JzJkz5ePjozfffFOvvvqqXF1d1b9/f/3hD39QdHS0cnNzVVZWpj//+c/685//XOk1Lib+Tpw4oQMHDlyyM11ZghAAUHfU1T7UpepX5TOyMomJiVq8eLEef/xxDRgwQE2aNJGLi4smTpxY5c/znxo/frxeeuklvf3223rggQeUkpKi48ePa9y4cfY69J3oO6F+IgEEXKXLzcf9+WsBAQGSpH//+98KDQ297HW9vLwqLJ4nXfoD4ErnBa9evVqFhYV67733ysWwc+fOSutf7HD9VHZ2tqTyH/bjxo3TuHHjVFhYqLS0ND399NMaNmyY9u3bp9DQUDVt2lTHjx+vcK1jx45J+l/bVMbLy8u++ONP/bwtAgICZLFYtHHjxkp373DUrmhNmzZVWVmZTp48Wa4jYxiGsrOzK3ROq2POdlXfx1/i5uamxMREJSYm6uzZs1q3bp1mzZqlwYMH6/Dhw2rSpIlcXV01evRoTZkypdJrhIWFSfrhfWjQoMEl59Jf7r0GANQfta0Pdan6VfmMrMybb76pMWPGaP78+RXibdy4cZVi+6nw8HDdcMMNWr58uR544AEtX75cQUFBio2Ntdeh70TfCfUTCSDACQYPHiw3NzcdPHiwwtDmn2vTpo3eeecdFRcX2z98c3NztWnTJjVq1OiqY7j4gfrTD3TDMPTaa69VWj8/P1/vv/9+uaG9b731llxcXNS/f/8K9X18fBQXF6eSkhKNHDlSu3btUmhoqAYNGqRVq1bp2LFj9lE/0g/DxL29vXXTTTddMuY2bdooJydHJ06csH9DVVJSopSUlHL1hg0bpueff15Hjx7Vr3/96ytojaszaNAgvfDCC3rzzTc1ffp0e/m7776rwsJCDRo0qNrufVFV38eqaNy4se68804dPXpU06ZN06FDhxQeHq6YmBjt2LFD3bp1s3+rVplhw4Zp/vz5atq06WU7vAAAXKma0Ie6FG9v7yv+jKyMxWKpkGj5z3/+o6NHj+q6666zl/18VM2VGDdunB588EF9/vnn+uCDD5SYmFhuMW76TvSdUD+RAAJ+JiMjo8IOFpLUrl27CsNXr1SbNm00Z84czZ49W999951uu+02NWnSRCdOnNCWLVvk4+OjZ599VtIP06r+8pe/6N5779Vvf/tb5ebm6oUXXrjmjsutt94qDw8P3X333ZoxY4aKioq0dOlSnTlzptL6TZs21YMPPqisrCx16NBBycnJeu211/Tggw8qJCREkvTb3/5WDRo0UJ8+fdSyZUtlZ2crKSlJfn5+9m90nn76aa1du1YxMTF66qmn5O/vr3/+85/6z3/+oxdeeEF+fn6XjHnUqFF66qmndNddd+l3v/udioqK9Morr1TYir5Pnz66//77NW7cOG3dulX9+/eXj4+Pjh8/rs8//1xdu3bVgw8+eE3td7ENBw8erMcff1x5eXnq06ePfSeL66+/XqNHj77me1xJDFV5H39JfHy8IiIiFB0drWbNmun777/XwoULFRoaat8N5eWXX1bfvn3Vr18/Pfjgg2rTpo3y8/N14MABffDBB/b5/dOmTdO7776r/v37a/r06erWrZtsNpuysrL08ccf69FHH9WNN97osLYAANQsdbUPdTlX+hlZmWHDhmnFihXq1KmTunXrpm3btumPf/xjhR2+2rVrpwYNGuif//ynOnfurIYNGyooKKjcF2s/d/fddysxMVF33323iouLK6zBQ9+JvhPqKRMXoAZqlMvtYCHJeO211+x1JRlTpky55DV+vpL/RatXrzZiYmKMRo0aGZ6enkZoaKhx5513GuvWrStX7/XXXzc6d+5seHl5GeHh4cbKlSsvuYPFH//4xyt+xg8++MDo3r274eXlZbRq1cr43e9+Z3z44YeGJOOzzz6z1xswYIDRpUsXY/369UZ0dLTh6elptGzZ0pg1a5ZRWlpaLs6YmBgjMDDQ8PDwMIKCgoxf//rXxjfffFPuvunp6UZ8fLzh5+dneHh4GN27d6+wk8WldsdITk42IiMjjQYNGhht27Y1Fi1aVGHHi4v+/ve/GzfeeKPh4+NjNGjQwGjXrp0xZswYY+vWrZdtl4vXO3nyZLnyi+9nZmamvezChQvG448/boSGhhru7u5Gy5YtjQcffNA4c+ZMuXNDQ0ONoUOHVrjXxZ0s3nnnnUrv9fN/O5XFdqXv45XsZPHSSy8ZvXv3NgICAgwPDw8jJCTEmDBhgnHo0KFy52VmZhrjx483WrVqZbi7uxvNmjUzevfubcydO7dcvYKCAuOJJ54wOnbsaHh4eBh+fn5G165djenTpxvZ2dkV2gMAUPvVhz7UpeK+eL1f+oysrJ9z5swZY8KECUbz5s0Nb29vo2/fvsbGjRsr3c3sX//6l9GpUyfD3d3dkGQ8/fTThmFU3AXsp+655x5DktGnT59LPhd9J/pOqF8shmEY1ZdeAgAAAAAAgNnYBh4AAAAAAKCOIwEEAAAAAABQx5EAAgAAAAAAqONIAAEAAAAAANRxJIAAAAAAAADqOBJAAAAAAAAAdZyb2QFUN5vNpmPHjsnX11cWi8XscAAAwGUYhqH8/HwFBQXJxYXvqcxC/wkAgNqhKn2nOp8AOnbsmFq3bm12GAAAoAoOHz6s4OBgs8Oot+g/AQBQu1xJ36nOJ4B8fX0l/dAYjRo1MjkaAABwOXl5eWrdurX98xvmoP8EAEDtUJW+U51PAF0cttyoUSM6MAAA1BJMOzIX/ScAAGqXK+k7MbkeAAAAAACgjiMBBAAAAAAAUMeRAAIAAAAAAKjjSAABAAAAAADUcSSAAAAAAAAA6jgSQAAAAAAAAHUcCSAAAAAAAIA6jgQQAAAAAABAHUcCCAAAAAAAoI4jAQQAAKqsuMxqdggAAACoAjezAwAAALVHUalVr244qHe2HlHy1H7y83Y3OyTUclaboa2HTuubI+eUdfq8zl0oVVGpVUVlNpVZbTIMyZAhST/+b/34f34oNwzTQgcAoMpeGxOtJj4eptybBBAAALgin+8/pSfXZCjzVKEk6d3tRzS+b5jJUaE225udr0fe3qE92flmhwIAgFOU2mym3ZsEEAAAuKyT+cWa95/dWr3zmCSpma+nno4P19CuLU2ODLXZmcIS/eav/9WpgmL5ermpX/sAtWnqI38fDzXwcJWXm6vcXC2yWCyy/HiOxSJdPPrhf/+vHACA2qCRl3mjp0kAAQCAStlshv71VZb+8OEe5RWVyWKRxtwUqkcHdzS184K6YfmmQzpVUKx2zXy08oFeCmjoaXZIAADUaSSAAABABbuP5Wn26nTtyDorSeoS1EjzE7qqe+vGpsaFuuOjjOOSpIdvbk/yBwAAJyABBAAA7AqLy7Rw3T79/YtDstoMNfR006OxHTT6plC5ubJ5KBzjVEGx9p0okMUiDezYzOxwAACoF0gAAQAASVLKrmw98/4uHT9XJEka0rWFnhrWRS38vEyODHXNwZwCSVLrJt5q7G3OTigAANQ3JIAAAKjnjpw5r2fe3611356QJAU3aaDnRkQoplNzkyNDXXVxJ7k2AT4mRwIAQP1BAggAgHqq1GrT8i8y9afU/bpQapWbi0X392+rh29urwYermaHhzrs+9PnJUltmnqbHAkAAPUHCSAAAOqhbd+f1uxVGdqTnS9J6tmmieYldFWHQF+TI0N9cKawRJLUjMWfAQBwGhJAAADUI2fPl+gPH+3Vv7ZkSZKaeLtr5pDOurNHsFxcLCZHh/ri3IVSSZKft7vJkQAAUH+QAAIAoB4wDEOrdx7V3LXfKvfH0Rd3RgVr1pDO8vdhEV44lz0B1IAEEAAAzkICCACAOu7gyQI9sSpDm7/LlSRd17yh5o6M0E1tm5ocGeqriwmgRiSAAABwGhJAAADUUUWlVi1Zf1Cvrj+oEqtNnm4umjqovX7br6083FzMDg/1GCOAAABwPlN6fwkJCWrSpInuvPNOe9nevXsVGRlp/2nQoIFWr15d7rzz588rNDRUjz32mJMjBgCgdtm4/6RuW5imVz7ZrxKrTQM6NFPq9AGaEnMdyR+YjgQQAADOZ8oIoKlTp2r8+PF6/fXX7WUdO3bUzp07JUkFBQVq06aNbr311nLnzZs3TzfeeKMzQwUAoFbJyS/S3LXf6v2vj0mSmvt66un4LhrStYUsFhZ5hvkMw1BBcZkkydeLwegAADiLKZ+6MTExWr9+/SVff//99zVo0CD5+PjYy/bv3689e/YoPj5eGRkZTogSAIDaw2oz9NaWLL3w0R7lF5XJxSKN6dVGj8Z2kK8XoyxQc5RaDRnGD//b083V3GAAAKhHqjwGPC0tTfHx8QoKCpLFYqkwTUuSlixZorCwMHl5eSkqKkobN26s0j3+7//+T6NGjSpX9thjjykpKamq4QIAUOdlHD2n25du0pOrM5RfVKaurfy0ekofPTO8C8kf1DilVpv9f3u4Mh0RAABnqfKnbmFhobp3765FixZV+vrKlSs1bdo0zZ49Wzt27FC/fv0UFxenrKysK7p+Xl6evvjiCw0ZMsRetmbNGnXo0EEdOnSoargAANRZBcVlem7tbg1f9Lm+PnxWDT3d9OzwLlo9pY+6BTc2OzygUj9NALm7Mi0RAABnqfIUsLi4OMXFxV3y9QULFmjChAmaOHGiJGnhwoVKSUnR0qVLr2gEz5o1azR48GB5eXnZy7788ku9/fbbeuedd1RQUKDS0lI1atRITz31VIXzi4uLVVxcbD/Oy8uryuMBAFDjGYahlF0n9OwHu3T8XJEkaWjXlnoqPlyBjbx+4WzAXCU/JoAsFsnVhQQQAADO4tA1gEpKSrRt2zb9/ve/L1ceGxurTZs2XdE1/u///k/3339/ubKkpCR78mjFihXKyMioNPlzse6zzz57FdEDAFDzHTlzXk+v2aVP9uRIklr7N9CcERGK6djc5MiAK1Nq/WEBIA9XFxYmBwDAiRyaADp16pSsVqsCAwPLlQcGBio7O9t+PHjwYG3fvl2FhYUKDg7WqlWr1LNnT507d05btmzRu+++e9UxzJw5U4mJifbjvLw8tW7d+qqvBwBATVBqtelvn2fq5XX7daHUKndXi+7v31YPxbRXAw8W0kXtUVL2wwgg1v8BAMC5qmUXsJ9/m2MYRrmylJSUSs/z8/PTiRMnLnvt++6777Kve3p6ytPT88oCBQCgFth66LRmr8rQ3hP5kqQbwvw1b2SE2gf6mhwZUHUX1wBydyMBBACAMzk0ARQQECBXV9dyo30kKScnp8KoIAAAcHlnz5fo+Q/36O2vDkuSmni7a9aQzrozKpipM6i1Lo4AYgFoAACcy6FfvXh4eCgqKkqpqanlylNTU9W7d29H3goAgDrLMAy9u+2Ibn5pgz358+voYH366ED9Kro1yR/UahdHAHkwAggAAKeq8giggoICHThwwH6cmZmpnTt3yt/fXyEhIUpMTNTo0aMVHR2tXr16admyZcrKytKkSZMcGjgAAHXRgZwCPbE6XV9+d1qS1L55Q81L6KobwvxNjgxwjIuLQLuzBhAAAE5V5QTQ1q1bFRMTYz++uODy2LFjtWLFCo0aNUq5ubmaM2eOjh8/roiICCUnJys0NNRxUQMAUMcUlVq1+LMDenXDQZVaDXm5u2jqoPaa2LctIyVwzdauXatHH31UNptNjz/+uCZOnGhaLCwCDQCAOaqcABo4cKAMw7hsncmTJ2vy5MlXHRQAAPVJ2r6TenJNhr7PPS9JiunYTHNGRKi1v7fJkaEuKCsrU2Jioj777DM1atRIPXr00O233y5/f3NGldkXgSYBBACAU1XLLmAAAOCX5eQVac7a3Vr7zXFJUmAjTz0T30W3RbRgnR84zJYtW9SlSxe1atVKkjRkyBClpKTo7rvvNiWeEiuLQAMAYAa+egEAwMmsNkNvbD6kQS9t0NpvjsvFIo3r00brEgcormtLkj+1RFJSknr27ClfX181b95cI0eO1N69ex16j7S0NMXHxysoKEgWi0WrV6+utN6SJUsUFhYmLy8vRUVFaePGjfbXjh07Zk/+SFJwcLCOHj3q0DirgkWgAQAwB5+8AAA4UcbRc7p9yRd6as0u5ReXqVuwn95/qK+eju8iXy93s8NDFWzYsEFTpkzRl19+qdTUVJWVlSk2NlaFhYWV1v/iiy9UWlpaoXzPnj3Kzs6u9JzCwkJ1795dixYtumQcK1eu1LRp0zR79mzt2LFD/fr1U1xcnLKysiSp0qn7ZiYZmQIGAIA5mAIGAIATFBSXacHH+7RiU6ZshuTr6abf3dZRv7kxVK4ujPipjT766KNyx8uXL1fz5s21bds29e/fv9xrNptNU6ZMUfv27fX222/L1dVVkrRv3z7FxMRo+vTpmjFjRoV7xMXFKS4u7rJxLFiwQBMmTLAv7Lxw4UKlpKRo6dKlSkpKUqtWrcqN+Dly5IhuvPHGq3pmR2ARaAAAzMEnLwAA1cgwDH2Yfly3vLRBf//ih+TPsG4t9cmjAzSmVxuSP3XIuXPnJKnSxZVdXFyUnJysHTt2aMyYMbLZbDp48KBuvvlmDR8+vNLkz5UoKSnRtm3bFBsbW648NjZWmzZtkiTdcMMNysjI0NGjR5Wfn6/k5GQNHjy40ustXrxY4eHh6tmz51XFc0Uxsw08AACmYAQQAADV5PDp83pqTYY+23tSkhTi763nRkZoQIdmJkcGRzMMQ4mJierbt68iIiIqrRMUFKRPP/1U/fv31z333KPNmzdr0KBBevXVV6/6vqdOnZLValVgYGC58sDAQPu0Mjc3N7300kuKiYmRzWbTjBkz1LRp00qvN2XKFE2ZMkV5eXny8/O76rgup/THEUDurAEEAIBTkQACAMDBSq02vbbxO73yyX4Vldrk7mrRpAHtNCXmOnm5u5odHqrBQw89pG+++Uaff/75ZeuFhITojTfe0IABA9S2bVv97W9/c8h6PD+/hmEY5cqGDx+u4cOHX/N9HMFq+2EEkBuj3wAAcCq+egEAwIG+OnRaQ1/ZqBc+2quiUptuauuvDx/pr0djO5L8qaMefvhhvf/++/rss88UHBx82bonTpzQ/fffr/j4eJ0/f17Tp0+/pnsHBATI1dW1wiLSOTk5FUYF1RSGfkgAsdkdAADOxQggAAAc4ExhiZ7/cI9Wbj0sSfL38dDsIZ11e49WbOteRxmGoYcfflirVq3S+vXrFRYWdtn6p06d0qBBg9S5c2e988472r9/vwYOHChPT0+9+OKLVxWDh4eHoqKilJqaqoSEBHt5amqqRowYcVXXrG4/DgCSC/9dAADgVCSAAAC4BoZh6N/bjmh+8rc6c/6HLb7v6tlaj9/WSU18PEyODtVpypQpeuutt7RmzRr5+vraR+H4+fmpQYMG5erabDbddtttCg0N1cqVK+Xm5qbOnTtr3bp1iomJUatWrSodDVRQUKADBw7YjzMzM7Vz5075+/srJCREkpSYmKjRo0crOjpavXr10rJly5SVlaVJkyZV49NfPduP29KT/gEAwLlIAAEAcJUO5ORr1qoMbck8LUnqGOireQkRim5TcRco1D1Lly6VJA0cOLBc+fLly3XfffeVK3NxcVFSUpL69esnD4//JQa7du2qdevWXXJR5q1btyomJsZ+nJiYKEkaO3asVqxYIUkaNWqUcnNzNWfOHB0/flwRERFKTk5WaGjoNT5h9TAYAQQAgClIAAEAUEVFpVYt+vSA/pJ2UKVWQ17uLpp2SwdN6BvG1tb1iHExk3GFbr311krLIyMjL3nOwIEDr+g+kydP1uTJk6sUj1kuPg/5HwAAnIsEEAAAVbB+b46eWrNLWafPS5IGdWquZ4Z3UWt/b5MjA2qHi2sAsTYWAADORQIIAIArcCKvSHPW7tZ/vjkuSWrRyEvPDO+iwV0C+UMWqIL/TQEzNw4AAOobEkAAAFyG1WboH5sP6cWP96mguEwuFmlcnzBNv7WDGnryMQpUlY0pYAAAmIKeKwAAl5B+5JxmrUpX+tFzkqTurRtr3sgIRbTyMzkyoPa6uAYQi0ADAOBcJIAAAPiZ/KJSvfTxPr2x+ZBshuTr5aYZt3XSPTeEyJV5K8A1ubikNQkgAACciwQQAAA/MgxDyenZevaDXcrJL5YkDe8epCeGdVZzXy+TowPqBlsVd08DAACOQQIIAABJWbnn9dT7GVq/96QkqU1Tbz03MkL92jczOTKgbrHZF4FmBBAAAM5EAggAUK+VlNn02sbv9Mon+1VcZpOHq4smDWynyQPbycvd1ezwgDqHXcAAADAHCSAAQL313+9yNXt1hg7kFEiSerVtqrkJEWrXrKHJkQF1l8EuYAAAmMLFjJsmJCSoSZMmuvPOO8uVu7m5KTIyUpGRkZo4ceIVnQMAQFWdLizR7975WqOWfakDOQVq6uOhP43qrrd+eyPJH6Ca2dgFDAAAU5gyAmjq1KkaP368Xn/99XLljRs31s6dO6t0DgAAV8owDL2z9Yjmf/itzp4vlSTdfUOIHr+toxp7e5gcHVA/XJwCZiEBBACAU5mSAIqJidH69eur/RwAAC7adyJfT6zK0JZDpyVJnVr4al5ChKJC/U2ODKhfbPYEkLlxAABQ31R5ClhaWpri4+MVFBQki8Wi1atXV6izZMkShYWFycvLS1FRUdq4ceMVXTsvL09RUVHq27evNmzYUNXQAACo4EKJVS98tEdDXt6oLYdOq4G7q2YN6aQPHu5L8gcwwf+mgJkcCAAA9UyVRwAVFhaqe/fuGjdunO64444Kr69cuVLTpk3TkiVL1KdPH/3lL39RXFycdu/erZCQkMte+9ChQwoKClJGRoaGDh2q9PR0NWrUqKohAgAgSfpsT46eXJOhI2cuSJJu6dxczwzvouAm3iZHBoA1gAAAcK4qJ4Di4uIUFxd3ydcXLFigCRMm2BdxXrhwoVJSUrR06VIlJSVd9tpBQUGSpIiICIWHh2vfvn2Kjo6uUnzFxcUqLi62H+fl5VXpfABA7Zd9rkhz1u5Scnq2JCnIz0vPDO+i2C4tTI4MwMURQKR/AABwLofuAlZSUqJt27YpNja2XHlsbKw2bdp02XPPnDljT9wcOXJEu3fvVtu2bascQ1JSkvz8/Ow/rVu3rvI1AAC1k9VmaPkXmbplwQYlp2fL1cWi3/YLU2riAJI/QA1hTwAxAggAAKdy6CLQp06dktVqVWBgYLnywMBAZWdn248HDx6s7du3q7CwUMHBwVq1apVKS0v1wAMPyMXFRRaLRS+//LL8/f0ve07Pnj0rxDBz5kwlJibaj/Py8kgCAUA98M2Rs5q1Kl0ZR38Y+RnZurHmJ3RVeBBTiYGa5OIuYEwBAwDAuaplF7Cff6NjGEa5spSUlErPS09Pv+Q1L3XOz3l6esrT0/OK6gIAar+8olK9lLJXb3z5vQxD8vVy0+O3ddI9N4TIhVVmgRqHXcAAADCHQxNAAQEBcnV1LTfaR5JycnIqjAoCAOBaGIah/6Qf15wPdisn/4cpxCMjgzR7aLia+fJFAFBTGewCBgCAKRyaAPLw8FBUVJRSU1OVkJBgL09NTdWIESMceSsAQD32fW6hnlqzSxv2nZQkhQX4aO7ICPW5LsDkyAD8EsM+AogMEAAAzlTlBFBBQYEOHDhgP87MzNTOnTvl7++vkJAQJSYmavTo0YqOjlavXr20bNkyZWVladKkSQ4NHABQ/xSXWfVa2nf686cHVFxmk4erix4c2E4PDmwnL3dXs8MDcAX+twi0yYEAAFDPVDkBtHXrVsXExNiPLy64PHbsWK1YsUKjRo1Sbm6u5syZo+PHjysiIkLJyckKDQ11XNQAgHrny+9yNXtVug6eLJQk9bmuqZ4bEaG2zRqaHBmAqrCxCDQAAKaocgJo4MCB9rnblzJ58mRNnjz5qoMCAOCi3IJizU/eo3e3H5EkBTT00BNDwzUiMogpJEAtZOjHEUAmxwEAQH1TLbuAAQBwrWw2Q+9sO6ykD/fo7PlSWSzSPTeEaMbgTvLzdjc7PABXiW3gAQAwBwkgAECNszc7X0+sTtdXh85Ikjq18NX827uqR0gTkyMDcK1YAwgAAHOQAAIA1BjnS8r0yicH9NeN36nMZsjbw1XTb+mgcX3ayM3VxezwADgAu4ABAGAOEkAAgBrh0z0n9OTqXTp69oIkKTY8UE8P76JWjRuYHBkAR7o4AsiF/A8AAE5FAggAYKrj5y7o2fd366Nd2ZKkID8vPTsiQreGB5ocGYDqwBpAAACYgwQQAMAUZVabXt/8vRZ8vFeFJVa5ulg0oW+YHhnUXj6efDwBdZV9FzDyPwAAOBU9bACA0+08fFazV6Vr17E8SVKPkMaal9BVnVs2MjkyANXNZvvh/7MGEAAAzkUCCADgNHlFpfrjR3v15n+/l2FIjbzc9Pu4zrqrZ2u5sCAIUC+wBhAAAOYgAQQAqHaGYeiDb47rubW7dTK/WJJ0+/WtNGtoZwU09DQ5OgDO9OMSQLKIDBAAAM5EAggAUK0OnSrUk2sytHH/KUlS2wAfzR0Zod7XBZgcGQAzGIwAAgDAFCSAAADVorjMqr9s+E6LPjugkjKbPNxcNGXgdZo0sK083VzNDg+ASWzsAgYAgClIAAEAHG7TwVN6YnWGvjtZKEnqe12AnhsZobAAH5MjA2C2iyOAmAEGAIBzkQACADjMqYJizU/+Vu9tPypJCmjoqSeHddbw7kHs+ANAEiOAAAAwCwkgAMA1s9kMrdx6WM9/uEfnLpTKYpHuvTFUjw3uKL8G7maHB6AGYRcwAADMQQIIAHBN9mTnafaqDG37/owkKbxlI81LiND1IU1MjgxATcYAIAAAnIsEEADgqpwvKdPL6/brr59nymoz5O3hqsRbO+i+3m3k5upidngAaqj/jQAiAwQAgDORAAIAVNm63Sf09Pu7dPTsBUnSbV1a6Kn4cAU1bmByZABqOpvN7AgAAKifSAABAK7YsbMX9OwHu5Sy64QkqVXjBpozoosGdQ40OTIAtYUhRgABAGAGEkAAgF9UZrVpxaZDWpC6T+dLrHJzsWhCvzA9Mqi9vD34KAFw5dgFDAAAc9BrBwBc1o6sM5q1KkPfHs+TJEWFNtG8hAh1atHI5MgA1EbGj2sAkf8BAMC5TFmlMyEhQU2aNNGdd95pLzt8+LAGDhyo8PBwdevWTe+88065c/70pz+pS5cuCg8P19SpU+2dBwBA9Th3oVRPrE7X7Us36dvjefJr4K7nb++qdx7oRfIHwFUz7COAzI0DAID6xpQRQFOnTtX48eP1+uuv/y8QNzctXLhQkZGRysnJUY8ePTRkyBD5+Pjo5MmTWrRokXbt2iV3d3f1799fX375pXr16mVG+ABQpxmGofe/Pqbn1n6rUwXFkqQ7egRr1pBOatrQ0+ToANR2NvsIIDJAAAA4kykJoJiYGK1fv75cWcuWLdWyZUtJUvPmzeXv76/Tp0/Lx8dHklRWVqaioiJJUmlpqZo3b+7UmAGgPsg8VagnV2fo8wOnJEltm/lo3siu6tWuqcmRAagrLq4BRPoHAADnqvIUsLS0NMXHxysoKEgWi0WrV6+uUGfJkiUKCwuTl5eXoqKitHHjxirdY+vWrbLZbGrdurUkqVmzZnrssccUEhKioKAg3XLLLWrXrl1VQwcAXEJxmVUL1+3T4IVp+vzAKXm4uejRWzvow0f6kfwB4FAXJ/GzCDQAAM5V5QRQYWGhunfvrkWLFlX6+sqVKzVt2jTNnj1bO3bsUL9+/RQXF6esrKwrun5ubq7GjBmjZcuW2cvOnDmjtWvX6tChQzp69Kg2bdqktLS0qoYOAKjEpgOnFLdwoxau26+SMpv6tQ/Qx9P66+FB7eXp5mp2eADqmIvrOLqYshIlAAD1V5WngMXFxSkuLu6Sry9YsEATJkzQxIkTJUkLFy5USkqKli5dqqSkpMteu7i4WAkJCZo5c6Z69+5tL1+3bp2uu+46+fv7S5KGDh2qL7/8Uv3796/0GsXFxfbjvLy8Kj0fANQXJ/OLNT/5W63acVSS1MzXU08NC9ewbi1ZmwNAtTHsU8D4PQMAgDM59LuXkpISbdu2TbGxseXKY2NjtWnTpsueaxiG7rvvPt18880aPXp0uddat26tTZs2qaioSFarVevXr1fHjh0rvU5SUpL8/PzsPxenkQEAfmCzGfrnf7/XoJfWa9WOo7JYpDG9QvXJowMU3z2I5A+AamVjG3gAAEzh0EWgT506JavVqsDAwHLlgYGBys7Oth8PHjxY27dvV2FhoYKDg7Vq1SoVFxdr5cqV6tatm31doX/84x/q2rWrbrrpJg0ZMkTXX3+9XFxcNGjQIA0fPrzSGGbOnKnExET7cV5eHkkgAPjRt8fzNGtVunZknZUkdQlqpPkJXdW9dWNT4wJQf9js28CTAQIAwJmqZRewn397bBhGubKUlJRKz7PZbJe85rx58zRv3rxfvLenp6c8PdmmGAB+qrC4TC9/sl9/+zxTVpshHw9XPRrbUWN6hcrNlYU4ADiPwQggAABM4dAEUEBAgFxdXcuN9pGknJycCqOCAADO8fGubD3z/i4dO1ckSYqLaKGn47uohZ+XyZEBqI8MRgABAGAKhyaAPDw8FBUVpdTUVCUkJNjLU1NTNWLECEfeCgDwC46evaBn3t+l1N0nJEnBTRpozoguurkTCXkA5mENIAAAzFHlBFBBQYEOHDhgP87MzNTOnTvl7++vkJAQJSYmavTo0YqOjlavXr20bNkyZWVladKkSQ4NHABQuVKrTcu/yNSfUvfrQqlVbi4W/bZ/W029ub0aeLCtOwBz/TgAiF3AAABwsiongLZu3aqYmBj78cUFl8eOHasVK1Zo1KhRys3N1Zw5c3T8+HFFREQoOTlZoaGhjosaAFCp7VlnNOu9dO3Jzpck9WzTRPMSuqpDoK/JkQHADy6OAHIh/wMAgFNVOQE0cOBA++J9lzJ58mRNnjz5qoMCAFTNufOl+kPKHv1rS5YMQ2rs7a5ZcZ11Z1SwXPgrC0ANYl8DiN9NAAA4VbXsAgYAcA7DMLRm5zHN/c9unSookSTdGRWsWUM6y9/Hw+ToAKAi+y5gJscBAEB9QwIIAGqp704W6Mk1GfriQK4k6brmDTV3ZIRuatvU5MgA4NJsP44AsrAKNAAATkUCCABqmaJSq5auP6il6w+qxGqTp5uLpg5qr9/2aysPNxezwwOAy2IXMAAAzEECCABqkc/3n9KTazKUeapQkjSgQzM9NyJCIU29TY4MAK6MfQ0gMkAAADgVCSAAqAVO5hdr7n92a83OY5Kk5r6eejq+i4Z0bcE0CgC1isEuYAAAmIIEEADUYDabobe2ZOkPH+1RflGZLBZpbK82SoztoEZe7maHBwBVZl8DiGWgAQBwKhJAAFBD7T6Wp1mr0rXz8FlJUkSrRpqf0FXdghubGhcAXAtDrAEEAIAZSAABQA1TWFymP6Xu0/JNh2S1GWro6abHYjtodK82cmXOBIBa7uIaQAAAwLlIAAFADZKyK1vPvL9Lx88VSZKGdm2pJ4eFq4Wfl8mRAYBjMQIIAADnIgEEADXAkTPn9cz7u7Xu2xOSpNb+DTRneIRiOjU3OTIAcCwGAAEAYA4SQABgolKrTX//PFML1+3XhVKr3F0tur9/Wz0U014NPFzNDg8AHM5gEWgAAExBAggATLLt+9OavSpDe7LzJUk3tPHXvIQItQ/0NTkyAKh+TAEDAMC5SAABgJOdPV+iP3y0R//acliS1MTbXTOHdNavooJl4S8iAHUek8AAADADCSAAcBLDMLRqx1HN+8+3yi0skST9OjpYv4/rLH8fD5OjAwDnsE8BI98NAIBTkQACACc4eLJAT6zK0ObvciVJ7Zs31NyREbqxbVOTIwMAc7AGEAAAzkUCCACqUVGpVUs+O6BXN3ynEqtNXu4umjqovSb2bSsPNxezwwMAp2MCGAAA5iABBADVZOP+k3pidYa+zz0vSRrYsZmeGxGh1v7eJkcGAOYxfpwDxhQwAACciwQQADhYTn6R5q79Vu9/fUySFNjIU0/Hd1FcRAsWeQaAH/HbEAAA5yIBBAAOYrUZeuu/3+uFlL3KLyqTi0Ua27uNEm/tIF8vd7PDA4AagSlgAACYgwQQADhAxtFzmr06Q18fPitJ6hbsp3kju6prsJ+5gQFADcWASAAAnMuUFUgTEhLUpEkT3XnnneXK165dq44dO6p9+/b661//Wu41Nzc3RUZGKjIyUhMnTnRmuABwSQXFZZrzwW4NX/S5vj58Vg093fTs8C5aNbkPyR8AqIRhHwJEBggAAGcyZQTQ1KlTNX78eL3++uv2srKyMiUmJuqzzz5To0aN1KNHD91+++3y9/eXJDVu3Fg7d+40I1wAqMAwDKXsytYz7+9Wdl6RJGlYt5Z6cli4Aht5mRwdANRchsEkMAAAzGDKCKCYmBj5+vqWK9uyZYu6dOmiVq1aydfXV0OGDFFKSooZ4QHAZR0+fV4TX9+qSW9uV3ZekUL8vbViXE8tuqcHyR8AuEJMAQMAwLmqnABKS0tTfHy8goKCZLFYtHr16gp1lixZorCwMHl5eSkqKkobN278xeseO3ZMrVq1sh8HBwfr6NGj9uO8vDxFRUWpb9++2rBhQ1XDBoBrVmq1aen6g7r1Txv0yZ4cubta9PDN1+nj6f01sGNzs8MDgFqB8T8AAJijylPACgsL1b17d40bN0533HFHhddXrlypadOmacmSJerTp4/+8pe/KC4uTrt371ZISMglr1vZcOCfbpd86NAhBQUFKSMjQ0OHDlV6eroaNWpU1fAB4KpsPXRas1ala9+JAknSjWH+mpcQoeua+/7CmQCAcn7s8jEACAAA56pyAiguLk5xcXGXfH3BggWaMGGCfaHmhQsXKiUlRUuXLlVSUtIlz2vVqlW5ET9HjhzRjTfeaD8OCgqSJEVERCg8PFz79u1TdHR0VcMHgCo5U1iiP3y0R29/dViS5O/joVlDOuuOHq3KJakBAFXD71AAAJzLoYtAl5SUaNu2bfr9739frjw2NlabNm267Lk33HCDMjIydPToUTVq1EjJycl66qmnJElnzpyRt7e3PD09deTIEe3evVtt27at9DrFxcUqLi62H+fl5V3jUwGojwzD0Lvbj2p+8rc6XVgiSbqrZ2s9flsnNfHxMDk6AKi9mAIGAIA5HJoAOnXqlKxWqwIDA8uVBwYGKjs72348ePBgbd++XYWFhQoODtaqVavUs2dPvfTSS4qJiZHNZtOMGTPUtGlTSdK3336rBx54QC4uLrJYLHr55Zftu4P9XFJSkp599llHPhaAeuZATr5mr8rQfzNPS5I6BDbUvISu6tmm8t87AIArd3HaP+N/AABwrmrZBv7nQ3oNwyhXdqndvYYPH67hw4dXKO/du7fS09Ov6N4zZ85UYmKi/TgvL0+tW7e+onMB1G9FpVYt+vSA/pJ2UKVWQ17uLnpkUAdN7Bcmd1dTNk0EgDqLGWAAADiXQxNAAQEBcnV1LTfaR5JycnIqjAqqLp6envL09HTKvQDUHRv2ndSTqzOUdfq8JOnmTs317PAuau3vbXJkAFC3MAUMAABzODQB5OHhoaioKKWmpiohIcFenpqaqhEjRjjyVgDgEDl5RZqzdrfWfnNcktSikZeeGR6uwV1asEApAFQDw74LGL9jAQBwpiongAoKCnTgwAH7cWZmpnbu3Cl/f3+FhIQoMTFRo0ePVnR0tHr16qVly5YpKytLkyZNcmjgAHAtrDZDb375vV5M2av84jK5WKRxfcI0/dYOauhZLbNjAQA/QY4dAADnqvJfOVu3blVMTIz9+OJ6O2PHjtWKFSs0atQo5ebmas6cOTp+/LgiIiKUnJys0NBQx0UNANcg4+g5zVqVrm+OnJMkdQ/207yEropo5WdyZABQ9xlMAgMAwBRVTgANHDjQvnvDpUyePFmTJ0++6qAAoDrkF5VqQeo+vb7pkGyG5Ovpphm3ddQ9N4bK1YWvogEAAADUXcxzAFDnGYahDzOy9ewHu3Qir1iSFN89SE8O7azmjbxMjg4A6hf7GkDk3QEAcCr2NQZQpx0+fV7jV3ylyf/crhN5xQpt6q03xt+gP999PckfAHXa2rVr1bFjR7Vv315//etfzQ7HjglgAACYgxFAAOqkkjKb/vr5d3rlk/0qKrXJ3dWiBwe00+SY6+Tl7mp2eABQrcrKypSYmKjPPvtMjRo1Uo8ePXT77bfL39/f7NDs2GkRAADnIgEEoM7ZknlaT6xO174TBZKkm9r6a+7IrrqueUOTIwMA59iyZYu6dOmiVq1aSZKGDBmilJQU3X333SZHJvsQINI/AAA4F1PAANQZpwtLNOPfX+vXf9msfScK1NTHQwt+3V3/+u1NJH8A1CppaWmKj49XUFCQLBaLVq9eXaHOkiVLFBYWJi8vL0VFRWnjxo32144dO2ZP/khScHCwjh496ozQfxG7gAEAYA4SQABqPcMw9M7Wwxr00nr939YjkqS7b2itTx4doNt7BDPNAECtU1hYqO7du2vRokWVvr5y5UpNmzZNs2fP1o4dO9SvXz/FxcUpKytLkirdsbWm/S6sYeEAAFDnMQUMQK22/0S+Zq/O0JbM05KkTi18NS8hQlGhNWedCwCoqri4OMXFxV3y9QULFmjChAmaOHGiJGnhwoVKSUnR0qVLlZSUpFatWpUb8XPkyBHdeOONl7xecXGxiouL7cd5eXkOeIrKVZKbAgAATsAIIAC10oUSq/6YskdDXtmoLZmn1cDdVTPjOumDh/uS/AFQp5WUlGjbtm2KjY0tVx4bG6tNmzZJkm644QZlZGTo6NGjys/PV3JysgYPHnzJayYlJcnPz8/+07p162qL/2L+x8IqQAAAOBUjgADUOuv35ujJNRk6fPqCJOmWzs31zPAuCm7ibXJkAFD9Tp06JavVqsDAwHLlgYGBys7OliS5ubnppZdeUkxMjGw2m2bMmKGmTZte8pozZ85UYmKi/TgvL69ak0ASU8AAAHA2EkAAao0TeUWa88Fu/Sf9uCSppZ+XnhneRbHhgTVubQsAqG4//71nGEa5suHDh2v48OFXdC1PT095eno6NL5LqWx9IgAAUP1IAAGo8aw2Q//YfEgvfrxPBcVlcnWxaFzvNpp+awf5ePJrDED9EhAQIFdXV/ton4tycnIqjAqqyUjbAwDgXPzlBKBGSz9yTrNWpSv96DlJUmTrxpqXEKEuQX4mRwYA5vDw8FBUVJRSU1OVkJBgL09NTdWIESNMjOzK2Mf/kAECAMCpSAABqJHyikq14ON9emPzIdkMydfLTY/f1kl33xAiVxf+agBQtxUUFOjAgQP248zMTO3cuVP+/v4KCQlRYmKiRo8erejoaPXq1UvLli1TVlaWJk2aZGLUV4YZYAAAmIMEEIAaxTAMJadn69kPdikn/4ctiUdEBmn20M5q7utlcnQA4Bxbt25VTEyM/fjiAs1jx47VihUrNGrUKOXm5mrOnDk6fvy4IiIilJycrNDQULNCrjJ2AQMAwLlIAAGoMbJyz+vJNRnasO+kJCkswEfPjYhQ3/YBJkcGAM41cODAX1wsefLkyZo8ebKTInI81u4HAMC5SAABMF1JmU2vbfxOr3yyX8VlNnm4uujBge304MB28nJ3NTs8AICDsAMYAADmIQEEwFT//S5Xs1dn6EBOgSSpd7umem5khNo1a2hyZACA6sQAIAAAnIsEEABTnC4s0fzkb/XvbUckSQENPfTE0HCNiAyShXkBAFAn/XQAEL/rAQBwLhJAAJzKZjP0721HNP/Db3X2fKkk6Z4bQ/T44E7y83Y3OToAQHViAhgAAOYhAQTAafadyNcTqzK05dBpSVKnFr6al9BVUaFNTI4MAOBsjP8BAMC5SAABqHYXSqx65dP9ei3tO5XZDDVwd1XirR10X582cnd1MTs8AICT/HQRaGaAAQDgXDXuL68//elP6tKli8LDwzV16lQZhqHDhw9r4MCBCg8PV7du3fTOO++YHSaAK/TZnhzd+qcNWrr+oMpshm4ND9S6Rwfot/3bkvwBgHqGKWAAAJinRo0AOnnypBYtWqRdu3bJ3d1d/fv315dffqk2bdpo4cKFioyMVE5Ojnr06KEhQ4bIx8fH7JABXEL2uSI9+8EufZiRLUkK8vPSM8O7KLZLC5MjAwDUBBYmgQEA4FQ1KgEkSWVlZSoqKpIklZaWqnnz5mrZsqVatmwpSWrevLn8/f11+vRpEkBADVRmtemNzd/rpY/3qrDEKlcXiyb0DdMjg9rLx7PG/coBADiRwRAgAABM49D5F2lpaYqPj1dQ0A/bOK9evbpCnSVLligsLExeXl6KiorSxo0b7a81a9ZMjz32mEJCQhQUFKRbbrlF7dq1K3f+1q1bZbPZ1Lp1a0eGDsABvj58ViOXfKE5a3ersMSq60Ma64OH+mrWkM4kfwCgFli8eLHCw8PVs2fP6r8ZA4AAAHAqhyaACgsL1b17dy1atKjS11euXKlp06Zp9uzZ2rFjh/r166e4uDhlZWVJks6cOaO1a9fq0KFDOnr0qDZt2qS0tDT7+bm5uRozZoyWLVvmyLABXKO8olI9tSZDI5d8oYyjeWrk5ab5CV317qTeCg9qZHZ4AIArNGXKFO3evVtfffVVtVzfEItAAwBgFod+JR8XF6e4uLhLvr5gwQJNmDBBEydOlCQtXLhQKSkpWrp0qZKSkrRu3Tpdd9118vf3lyQNHTpUX375pfr376/i4mIlJCRo5syZ6t279yXvUVxcrOLiYvtxXl6eg54OwM8ZhqG13xzXnLW7dTL/h//uEq5vpVlDOquZr6fJ0QEAahqmgAEAYB6nbcFTUlKibdu2KTY2tlx5bGysNm3aJElq3bq1Nm3apKKiIlmtVq1fv14dO3aUYRi67777dPPNN2v06NGXvU9SUpL8/PzsP0wVA6rH97mFGrv8Kz38rx06mV+ssAAf/XPijfrTqEiSPwCAX8QAIAAAnMtpCaBTp07JarUqMDCwXHlgYKCys3/YJeimm27SkCFDdP3116tbt25q166dhg8fri+++EIrV67U6tWrFRkZqcjISKWnp1d6n5kzZ+rcuXP2n8OHD1f7swH1SXGZVX/+ZL9i/5SmtH0n5eHmoum3dNCHj/RTn+sCzA4PAFBLWJgDBgCAUzl9Vdaff9gbhlGubN68eZo3b165On379pXNZrui63t6esrTk9EHQHXYfDBXT6xO18GThZKkvtcF6LmREQoLYEc+AMAvYwoYAADmcVoCKCAgQK6urvbRPhfl5ORUGBUEoGbJLSjWvORv9d72o5KkgIaeenJYZw3vHsQ3uACAq8KnBwAAzuW0KWAeHh6KiopSampqufLU1NTLLuoMwDw2m6G3t2Tp5pc26L3tR2WxSPfeFKJPHh2gEZGtSP4AAKqEXcAAADCPQ0cAFRQU6MCBA/bjzMxM7dy5U/7+/goJCVFiYqJGjx6t6Oho9erVS8uWLVNWVpYmTZrkyDAAOMDe7HzNXpWurd+fkSR1btlI8xIi1COkicmRAQBqK6aAAQBgHocmgLZu3aqYmBj7cWJioiRp7NixWrFihUaNGqXc3FzNmTNHx48fV0REhJKTkxUaGurIMABcg/MlZXr5k/3628ZMldkMeXu4KvHWDrqvdxu5uTpt0CAAoI6zMAkMAACncmgCaODAgTJ+4audyZMna/LkyY68LQAH+eTbE3pqzS4dPXtBkjS4S6Ceju+ioMYNTI4MAFAX/LSXyBQwAACcy+m7gAGoeY6fu6Bn39+tj3b9sEh7q8YN9OzwLrolnAXaAQCO80tfFAIAgOpDAgiox8qsNq3YdEh/St2nwhKrXF0smtg3TI/c0l7eHvx6AAAAAIC6gr/wgHpq5+GzmvVeunYfz5MkRYU20byECHVq0cjkyAAAdRVTwAAAMA8JIKCeOXehVH9M2aN//jdLhiH5NXDXzLhO+nV0a7m40BsHAAAAgLqIBBBQTxiGoQ++Oa7n1u7WyfxiSdLt17fSrKGdFdDQ0+ToAAD1wU+XAGIXMAAAnIsEEFAPHDpVqCfXZGjj/lOSpLbNfDR3ZIR6twswOTIAQL3CGtAAAJiGBBBQhxWXWfXq+u+0eP0BlZTZ5OHmoodirtMDA9rK083V7PAAAPUYawABAOBcJICAOmrTwVN6YlWGvjtVKEnq1z5Az42IUJsAH5MjAwDUVIsXL9bixYtltVqr5frGT4YAkf8BAMC5SAABdcypgmLN/8+3em/HUUlSQENPPRUfrvhuLWXh61YAwGVMmTJFU6ZMUV5envz8/Bx+fYMpYAAAmIYEEFBH2GyGVm49rOc/3KNzF0plsUj33hiqxwZ3lF8Dd7PDAwCgHL6UAADAuUgAAXXAt8fzNHtVurZnnZUkhbdspPm3d1Vk68amxgUAwE/9dAAQ6R8AAJyLBBBQi50vKdPL6/brr59nymoz5OPhqsTYjhrbK1Ruri5mhwcAQDkGc8AAADANCSCgllq3+4Sefn+Xjp69IEmKi2ihp+LD1dKvgcmRAQDwy5gBBgCAc5EAAmqZY2cv6Jn3d+nj3SckSa0aN9BzI7vo5k6BJkcGAMDllZsCRgYIAACnIgEE1BJlVptWbDqkBan7dL7EKjcXiyb2a6upg66Ttwf/KQMAaj5mgAEAYB7+agRqge1ZZzR7VYa+PZ4nSYoObaJ5CV3VsYWvyZEBAAAAAGoDEkBADXbufKleSNmjt7ZkyTCkxt7umhnXSb+Kai0XF4bOAwBqF+PHSWDM/gIAwPlIAAE1kGEYev/rY3pu7W6dKiiRJN3RI1izhnRS04aeJkcHAAAAAKhtSAABNcx3Jwv05JoMfXEgV5LUrpmP5o7sql7tmpocGQAA1+jHNYAYAAQAgPORAAJqiKJSq17dcFBLPjuoEqtNnm4uevjm63R//3bycHMxOzwAAK7ZxTWg2QEMAADnq3F/Vbq5uSkyMlKRkZGaOHGivfzFF19Uly5dFBERoTfffNPECAHH++LAKQ15eaMWrtuvEqtN/Ts008fT++uhm9uT/AEAAAAAXLMaNwKocePG2rlzZ7my9PR0vfXWW9q2bZskadCgQRo2bJgaN27s/AABBzqZX6x5/9mt1TuPSZKa+Xrq6fhwDe3akm9HAQB1jsEUMAAATFPjEkCV+fbbb9W7d295eXlJkiIjI/XRRx/prrvuMjky4OrYbIb+9VWW/vDhHuUVlclikcbcFKpHB3dUIy93s8MDAKBaGPZJYAAAwNkcOrckLS1N8fHxCgoKksVi0erVqyvUWbJkicLCwuTl5aWoqCht3Lix3Ot5eXmKiopS3759tWHDBklSRESEPvvsM509e1Znz57Vp59+qqNHjzoydMBpdh/L0x2vbtLsVRnKKypTRKtGWjOlj54dEUHyBwBQLzDIFQAA53PoCKDCwkJ1795d48aN0x133FHh9ZUrV2ratGlasmSJ+vTpo7/85S+Ki4vT7t27FRISIkk6dOiQgoKClJGRoaFDhyo9PV3h4eGaOnWqbr75Zvn5+alnz55yc6sVg5cAu8LiMi1ct09//+KQrDZDDT3d9GhsB43p1UauLvSEAQB13/+mgPG5BwCAszk0ixIXF6e4uLhLvr5gwQJNmDDBvrjzwoULlZKSoqVLlyopKUmSFBQUJOmHUT/h4eHat2+foqOj9cADD+iBBx6QJE2cOFHXXXddpfcoLi5WcXGx/TgvL88hzwZci493ZeuZ93fp2LkiSdKQri301LAuauHnZXJkAAA4DxPAAAAwj9O2FyopKdG2bdsUGxtbrjw2NlabNm2SJJ05c8aevDly5Ih2796ttm3bSpJycnIkSXv37tWWLVs0ePDgSu+TlJQkPz8/+0/r1q2r65GAX3TkzHlNfH2r7v/HNh07V6TgJg20/L6eWvKbKJI/AID6iwFAAAA4ndPmUZ06dUpWq1WBgYHlygMDA5WdnS3ph8WeH3jgAbm4uMhisejll1+Wv7+/JGnkyJE6e/asfHx8tHz58ktOAZs5c6YSExPtx3l5eSSB4HSlVpuWf5GpP6Xu14VSq9xcLLq/f1s9fHN7NfBwNTs8AAAqtXjxYi1evFhWq7Varm/8OAeM/A8AAM7n9IV0fr61tWEY9rLevXsrPT290vMujhL6JZ6envL09Ly2IIFrsO37M5q9Kl17svMlSTe08dfchAh1CPQ1OTIAAC5vypQpmjJlivLy8uTn52d2OAAAwIGclgAKCAiQq6urfbTPRTk5ORVGBQG10bnzpXr+oz3615YsSVITb3fNHNJZv4oKrpD4BACgPrIvAs3HIgAATue0NYA8PDwUFRWl1NTUcuWpqanq3bu3s8IAHM4wDK3acUQ3v7Tenvz5VVSwPnl0oH4d3ZrkDwAAP8MuYAAAOJ9DRwAVFBTowIED9uPMzEzt3LlT/v7+CgkJUWJiokaPHq3o6Gj16tVLy5YtU1ZWliZNmuTIMACnOXiyQE+uztCmg7mSpOuaN9S8kRG6sW1TkyMDAAAAAOB/HJoA2rp1q2JiYuzHFxdjHjt2rFasWKFRo0YpNzdXc+bM0fHjxxUREaHk5GSFhoY6Mgyg2hWVWrVk/UG9uv6gSqw2ebq5aOqg9vptv7bycHPawDoAAGoVpoABAGAehyaABg4caN/d4VImT56syZMnO/K2gFNt3H9ST67O0KHc85KkgR2bac7wCIU09TY5MgAAajZD7AIGAIBZnL4LGFBb5eQXae7ab/X+18ckSYGNPPV0fBfFRbRgnR8AAAAAQI1GAgj4BTaboX9uydILH+1RflGZXCzSmF5t9GhsB/l6uZsdHgAAtcb/poDxxQkAAM5GAgi4jF3Hzmn2qgztPHxWktQt2E/zRnZV12A/cwMDAKAWurhQAOkfAACcjwQQUImC4jL9KXWfln+RKZshNfR00+8Gd9S9N4XK1YVuKwAAAACgdiEBBPyEYRhK2XVCz36wS8fPFUmShnZrqaeGhSuwkZfJ0QEAULvZNwvhuxQAAJyOBBDwoyNnzuvpNbv0yZ4cSVKIv7fmjOiigR2bmxwZAAB1w+X3igUAANWJBBDqvVKrTX/7PFMvr9uvC6VWubta9ED/dnro5uvk5e5qdngAANQ5DAACAMD5SAChXtt66LRmr8rQ3hP5kqQbwvw1PyFC1zX3NTkyAADqHnYBAwDAPCSAUC+dPV+i5z/co7e/OixJauLtrtlDw3VHj1Z0SgEAAAAAdQ4JINQrhmHove1HNS/5W50uLJEkjYpurd/HdVITHw+TowMAoK77YQgQ37UAAOB8JIBQbxzIKdATq9P15XenJUkdAhtqXkJX9Wzjb3JkAADUD2wCBgCAeVzMDgCobkWlVr308V7FvZymL787LS93Fz1+WyetfbgfyR8AQL2QkJCgJk2a6M477zQ7FAAAYBJGAKFOS9t3Uk+uydD3ueclSTEdm2nOiAi19vc2OTIAAJxn6tSpGj9+vF5//XVT47i4DTzr7QEA4HwkgFAn5eQVac7a3Vr7zXFJUmAjTz0T30W3RbSg0wkAqHdiYmK0fv16s8NgChgAACZiChjqFKvN0D82H9KglzZo7TfH5WKRxvcJ0yePDlRc15YkfwAAVXL06FHde++9atq0qby9vRUZGalt27Y57PppaWmKj49XUFCQLBaLVq9eXWm9JUuWKCwsTF5eXoqKitLGjRsdFgMAAKgfGAGEOiPj6DnNXpWur4+ckyR1C/bT/ISuimjlZ3JkAIDa6MyZM+rTp49iYmL04Ycfqnnz5jp48KAaN25caf0vvvhCN9xwg9zd3cuV79mzR40bN1aLFi0qnFNYWKju3btr3LhxuuOOOyq97sqVKzVt2jQtWbJEffr00V/+8hfFxcVp9+7dCgkJkSRFRUWpuLi4wrkff/yxgoKCqvjk1cdgFzAAAExDAgi1XkFxmRZ8vE8rNmXKZki+nm763W0d9ZsbQ+XqQg8TAHB1/vCHP6h169Zavny5vaxNmzaV1rXZbJoyZYrat2+vt99+W66urpKkffv2KSYmRtOnT9eMGTMqnBcXF6e4uLjLxrFgwQJNmDBBEydOlCQtXLhQKSkpWrp0qZKSkiTJoaOSqtPFKWBMAgMAwPmYAoZayzAMfZh+XLe8tEF//+KH5E989yB98ugAjenVhuQPAOCavP/++4qOjtavfvUrNW/eXNdff71ee+21Suu6uLgoOTlZO3bs0JgxY2Sz2XTw4EHdfPPNGj58eKXJnytRUlKibdu2KTY2tlx5bGysNm3adFXXvJzFixcrPDxcPXv2dPi1AQCAuRgBhFrp8Onzevr9Xfp0T44kKbSpt+aMiNCADs1MjgwAUFd89913Wrp0qRITEzVr1ixt2bJFU6dOlaenp8aMGVOhflBQkD799FP1799f99xzjzZv3qxBgwbp1VdfveoYTp06JavVqsDAwHLlgYGBys7OvuLrDB48WNu3b1dhYaGCg4O1atWqSpM8U6ZM0ZQpU5SXlyc/P8dPobYvAs13NAAAOB0JINQqpVab/roxUy9/sk9FpTa5u1o0aUA7TYm5Tl7urmaHBwCoQ2w2m6KjozV//nxJ0vXXX69du3Zp6dKllSaAJCkkJERvvPGGBgwYoLZt2+pvf/ubQzYg+Pk1DMOo0nVTUlKuOQZHsK8BZHIcAADUR0wBQ63x1aHTGvrKRv3hoz0qKrXpprb++vCR/no0tiPJHwCAw7Vs2VLh4eHlyjp37qysrKxLnnPixAndf//9io+P1/nz5zV9+vRriiEgIECurq4VRvvk5ORUGBUEAABwOTUmAZSfn6+ePXsqMjJSXbt2LTfHPjMzUzExMQoPD1fXrl1VWFhoYqRwtjOFJXr839/oV69u1r4TBfL38dBLv+quf/32Jl3XvKHZ4QEA6qg+ffpo79695cr27dun0NDQSuufOnVKgwYNUufOnfXee+/p008/1f/93//pscceu+oYPDw8FBUVpdTU1HLlqamp6t2791Vf1yxMAQMAwDw1ZgqYt7e3NmzYIG9vb50/f14RERG6/fbb1bRpU913332aO3eu+vXrp9OnT8vT09PscOEEhmHo3e1HNT/5W50uLJEk3X1Daz1+Wyc19vYwOToAQF03ffp09e7dW/Pnz9evf/1rbdmyRcuWLdOyZcsq1LXZbLrtttsUGhqqlStXys3NTZ07d9a6desUExOjVq1aVToaqKCgQAcOHLAfZ2ZmaufOnfL397dv8Z6YmKjRo0crOjpavXr10rJly5SVlaVJkyZV38MDAIA6p8YkgFxdXeXt7S1JKioqktVqlWEY2rVrl9zd3dWvXz9Jkr+/v5lhwkkO5ORr9qoM/TfztCSpY6Cv5iVEKLoN7z8AwDl69uypVatWaebMmZozZ47CwsK0cOFC/eY3v6lQ18XFRUlJSerXr588PP73JUXXrl21bt06NW3atNJ7bN26VTExMfbjxMRESdLYsWO1YsUKSdKoUaOUm5urOXPm6Pjx44qIiFBycvIlRyLVBhZWAQIAwOkcNgUsLS1N8fHxCgoKksVi0erVqyvUWbJkicLCwuTl5aWoqCht3Lix3Otnz55V9+7dFRwcrBkzZiggIED79+9Xw4YNNXz4cPXo0cO+ECPqpqJSq15M2au4lzfqv5mn5eXuot/HddLaqX1J/gAAnG7YsGFKT09XUVGRvv32W/32t7+9ZN1bb71VXl5eFcojIyPVunXrSs8ZOHCgDMOo8HMx+XPR5MmTdejQIRUXF2vbtm3q37//NT2XWZgCBgCAeRyWACosLFT37t21aNGiSl9fuXKlpk2bptmzZ2vHjh3q16+f4uLiyi2k2LhxY3399dfKzMzUW2+9pRMnTqi0tFQbN27U4sWLtXnzZqWmplaYB4+6Yf3eHMX+KU2LPjugUquhQZ2aK3X6AE0a0E7urjVmuSoAAAAAAGodh00Bi4uLU1xc3CVfX7BggSZMmKCJEydKkhYuXKiUlBQtXbpUSUlJ5eoGBgaqW7duSktLU3BwsHr27Gn/5mzIkCHauXOnbr311krvU1xcrOLiYvtxXl7etT4aqtmJvCLNWbtb//nmuCSppZ+Xno7vosFdAh2ydS4AAKgZ2AYeAADzOGVYRUlJibZt26bY2Nhy5bGxsdq0aZOkH7ZNvZisycvLU1pamjp27KiePXvqxIkTOnPmjGw2m9LS0tS5c+dL3ispKUl+fn72n0sNuYb5rDZDr286pFte2qD/fHNcLhZpQt8wpSYO0G0RLUj+AABQx/xvChif8QAAOJtTFoE+deqUrFarAgMDy5UHBgYqOztbknTkyBFNmDDBPvf9oYceUrdu3SRJ8+fPV//+/WUYhmJjYzVs2LBL3mvmzJn2BRSlH5JJJIFqnvQj5zRrVbrSj56TJHVv3VjzRkYoopWfyZEBAAAAAFD3OHUXsJ9/22MYhr0sKipKO3furPS8X5pe9lOenp5sE1+D5ReV6qWP9+mNzYdkMyRfLzfNuK2T7rkhRK4ufBsIAEBdZpgdAAAA9ZhTEkABAQFydXW1j/a5KCcnp8KoINRNhmHow4xsPfvBLp3I+2GNphGRQZo9tLOa+1bcMQUAANQ9xo9zwJgBBgCA8zklAeTh4aGoqCilpqYqISHBXp6amqoRI0Y4IwSYKCv3vJ56P0Pr956UJLVp6q3nRkaoX/tmJkcGAAAAAED94LAEUEFBgQ4cOGA/zszM1M6dO+Xv76+QkBAlJiZq9OjRio6OVq9evbRs2TJlZWVp0qRJjgoBNUxJmU2vbfxOr3yyX8VlNnm4umjSwHaaPLCdvNxdzQ4PAAA42cUpYIwAAgDA+RyWANq6datiYmLsxxcXYh47dqxWrFihUaNGKTc3V3PmzNHx48cVERGh5ORkhYaGOioE1CBbMk9r9qp07c8pkCT1btdUz42MULtmDU2ODAAAmMW+CxgbwQMA4HQOSwANHDjQPq/7UiZPnqzJkyc76paogU4Xligp+Vu9s+2IJKmpj4eeGNZZIyNbseUrAAAAAAAmceouYKi7DMPQO9uOKCn5W505XypJuufGED0+uJP8vN1Njg4AANQMLAINAIBZSADhmu07ka8nVmVoy6HTkqROLXw1L6GrokKbmBwZAACoicj/AADgfCSAcNUulFj150/3a1nadyqzGWrg7qrpt7bXuD5hcnd1MTs8AABQw/zCagEAAKAakQDCVflsb46eWpOhw6cvSJJu6RyoZ0d0UavGDUyODAAA1FT/2wWMMUAAADgbCSBUSfa5Is1Zu0vJ6dmSpCA/Lz0zvItiu7QwOTIAAAAAAHApJIBwRaw2Q29sPqSXPt6nguIyubpYNL5PG027pYN8PPlnBAAAftn/toEHAADOxl/u+EXfHDmrWavSlXE0T5J0fUhjzRvZVeFBjUyODAAA1CYGGSAAAExDAgiXlFdUqpdS9uqNL7+XYUiNvNz0eFwn3d0zRC4u9NwAAAAAAKgtSAChAsMw9J/045rzwW7l5BdLkkZGBmn20HA18/U0OToAAFBb2ReBNjUKAADqJxJAKOf73EI9tWaXNuw7KUkKC/DR3JER6nNdgMmRAQCA2s4+A4xdwAAAcDoSQJAklZTZtCztoP786QEVl9nk4eqiyTHtNGlAO3m5u5odHgAAcILFixdr8eLFslqtZocCAAAcjAQQ9OV3uXpidYYO5BRIkvpc11TPjYhQ22YNTY4MAAA405QpUzRlyhTl5eXJz8/P4dc3fpwExvgfAACcjwRQPZZbUKz5yXv07vYjkqSAhh56cli4hncPYmg2AACoNnQzAABwPhJA9ZDNZuidbYeV9OEenT1fKotFuueGEM0Y3El+3u5mhwcAAOoq45erAACA6kECqJ7ZdyJfs1el66tDZyRJnVs20ryECPUIaWJyZAAAoK773y5gDAECAMDZSADVExdKrHrl0/16Le07ldkMeXu4KvHWDrqvdxu5ubqYHR4AAKhHmAIGAIDzkQCqBz7dc0JPrdmlI2cuSJJuDQ/UM8O7qFXjBiZHBgAA6hODKWAAAJiGBFAddvzcBc35YLc+zMiWJLVq3EDPDO+iW8MDTY4MAADURwaLAAEAYBoSQHVQmdWmNzZ/r5c+3qvCEqtcXSya2DdMUwe1l48nbzkAADAXu40CAOB8ZAPqmK8Pn9Xs1enKOJonSeoR0ljzErqqc8tGJkcGAADqO6aAAQBgnhq1+m9CQoKaNGmiO++8016Wn5+vnj17KjIyUl27dtVrr71mYoQ1V15RqZ5ak6GRS75QxtE8NfJy0/yErvr3pN4kfwAAQI3wv13AAACAs9WoEUBTp07V+PHj9frrr9vLvL29tWHDBnl7e+v8+fOKiIjQ7bffrqZNm5oYac1hGIbWfnNcc9bu1sn8YklSwvWtNHtoZwU09DQ5OgAAAAAAUBPUqARQTEyM1q9fX67M1dVV3t7ekqSioiJZrVYZjB+WJH2fW6gn1+xS2r6TkqS2AT6aOzJCva8LMDkyAACAii724VgCCAAA53PYFLC0tDTFx8crKChIFotFq1evrlBnyZIlCgsLk5eXl6KiorRx48YruvbZs2fVvXt3BQcHa8aMGQoIqN8JjuIyqxZ9ul+xf0pT2r6T8nBz0fRbOujDaf1I/gAAgBrLPgWMBBAAAE7nsARQYWGhunfvrkWLFlX6+sqVKzVt2jTNnj1bO3bsUL9+/RQXF6esrKxfvHbjxo319ddfKzMzU2+99ZZOnDjhqLBrnS+/y9WQlzfqxY/3qbjMpr7XBShlWn89ckt7ebq5mh0eAAAAAACogRw2BSwuLk5xcXGXfH3BggWaMGGCJk6cKElauHChUlJStHTpUiUlJV3RPQIDA9WtWzelpaXpV7/6VaV1iouLVVxcbD/Oy8urwlPUXLkFxZqfvEfvbj8iSQpo6KEnh4VrePcgtlIFAAC1w49DgCwsAw0AgNM5ZRewkpISbdu2TbGxseXKY2NjtWnTpsuee+LECXsSJy8vT2lpaerYseMl6yclJcnPz8/+07p162t/ABPZbIZWfpWlQQs26N3tR2SxSL+5MUSfJA7UiMhWJH8AAECtQ/cFAADnc8oi0KdOnZLValVgYGC58sDAQGVnZ9uPBw8erO3bt6uwsFDBwcFatWqVXFxcNGHCBBmGIcMw9NBDD6lbt26XvNfMmTOVmJhoP87Ly6u1SaB9J/I1e1W6vjp0RpLUuWUjzUuIUI+QJiZHBgAAUHWG2MgDAACzOHUXsJ+PVjEMo1xZSkpKpeft3Lnziu/h6ekpT8/avf35hRKrXvl0v15L+05lNkPeHq5KvLWD7uvdRm6uThm0BQAA4HCGfQoYAABwNqckgAICAuTq6lputI8k5eTkVBgVVN99uueEnlqzS0fOXJAkxYYH6unhXdSqcQOTIwMAAHAQ5oABAOB0TkkAeXh4KCoqSqmpqUpISLCXp6amasSIEc4IocY7fu6C5nywWx9m/JAkC/Lz0rMjInRrOAkyAABQNxjMAAMAwDQOSwAVFBTowIED9uPMzEzt3LlT/v7+CgkJUWJiokaPHq3o6Gj16tVLy5YtU1ZWliZNmuSoEGqlMqtNb2z+Xi99vFeFJVa5ulg0sW+Ypg5qLx9Pp87QAwAAqFYX8z+M/wEAwPkclmHYunWrYmJi7McXF2IeO3asVqxYoVGjRik3N1dz5szR8ePHFRERoeTkZIWGhjoqhFrn68NnNXt1ujKO/rDLWY+QxpqX0FWdWzYyOTIAAIDqwwwwAACcz2EJoIEDB8r4hXG9kydP1uTJkx11y1orr6hUL6bs1T++/F6GITXyctPv4zrrrp6t5eJCjwgAANRNv9RXBAAA1Yc5Rk5kGIbWfnNcc9bu1sn8YklSwvWtNHtoZwU0rN07lwEAAPwSpoABAGAeEkBO8n1uoZ5cs0tp+05KktoG+GjuyAj1vi7A5MgAAACcy8IcMAAAnI4EUDUrKbPptY3f6ZVP9qu4zCYPNxdNGXidJg1sK083V7PDAwAAcBpmgAEAYB4SQNVoR9YZ/f7ddO09kS9J6ntdgJ4bGaGwAB+TIwMAADDDDxkgxv8AAOB8LmYHUBcVFJfpmfd36falm7T3RL78fTz0p1Hd9Y8JN5D8AQAATpeQkKAmTZrozjvvNDsUAABgEhJADrZh30nFLtigFZsOyTCk23u00rrEAUq4Ppj57gAAwBRTp07VG2+8YXYY9ilgdIkAAHA+poA5SFGpVS98tFd//yJTktTav4Hmjeyq/h2amRwZAACo72JiYrR+/Xqzw7CzMAkMAACnYwSQA+TkFemOpZvsyZ+xvUKVMq0/yR8AAOqIpKQkWSwWTZs2zaHXTUtLU3x8vIKCgmSxWLR69epK6y1ZskRhYWHy8vJSVFSUNm7c6NA4nIU1oAEAMA8JoGt09OwF3b50k3Ydy5O/j4f+fl+0nh0RIW8PBlcBAFAXfPXVV1q2bJm6det22XpffPGFSktLK5Tv2bNH2dnZlZ5TWFio7t27a9GiRZe87sqVKzVt2jTNnj1bO3bsUL9+/RQXF6esrCx7naioKEVERFT4OXbs2BU+pXPYdwFjABAAAE5HluIaFJdZ9eCb23TkzAW1aeqtN8bfqJCm3maHBQAAHKSgoEC/+c1v9Nprr2nu3LmXrGez2TRlyhS1b99eb7/9tlxdXSVJ+/btU0xMjKZPn64ZM2ZUOC8uLk5xcXGXjWHBggWaMGGCJk6cKElauHChUlJStHTpUiUlJUmStm3bdrWPaAryPwAAOB8jgK7B/311WN8cOafG3u56cyLJHwAA6popU6Zo6NChuuWWWy5bz8XFRcnJydqxY4fGjBkjm82mgwcP6uabb9bw4cMrTf5ciZKSEm3btk2xsbHlymNjY7Vp06aruublLF68WOHh4erZs6fDry1JBpPAAAAwDSOArlKp1aZXN3wnSUq8tYOCm5D8AQCgLnn77be1fft2ffXVV1dUPygoSJ9++qn69++ve+65R5s3b9agQYP06quvXnUMp06dktVqVWBgYLnywMDAS04rq8zgwYO1fft2FRYWKjg4WKtWrao0yTNlyhRNmTJFeXl58vPzu+q4L4VdwAAAMA8JoKv02Z4cHT17Qc18PfXr6NZmhwMAABzo8OHDeuSRR/Txxx/Ly8vris8LCQnRG2+8oQEDBqht27b629/+JosDsh0/v4ZhGFW6bkpKyjXH4EjsAgYAgPMxBewq3RoeqLfvv0lzhneRl7ur2eEAAAAH2rZtm3JychQVFSU3Nze5ublpw4YNeuWVV+Tm5iar1VrpeSdOnND999+v+Ph4nT9/XtOnT7+mOAICAuTq6lphtE9OTk6FUUG1QXSbJnr13ihNv7WD2aEAAFDvMALoKlksFt3UtqnZYQAAgGowaNAgpaenlysbN26cOnXqpMcff9y+yPNPnTp1SoMGDVLnzp31zjvvaP/+/Ro4cKA8PT314osvXlUcHh4eioqKUmpqqhISEuzlqampGjFixFVd00wt/RqopV8Ds8MAAKBeIgEEAADwM76+voqIiChX5uPjo6ZNm1Yol37YBey2225TaGioVq5cKTc3N3Xu3Fnr1q1TTEyMWrVqVelooIKCAh04cMB+nJmZqZ07d8rf318hISGSpMTERI0ePVrR0dHq1auXli1bpqysLE2aNMnBTw0AAOoyEkAAAADXyMXFRUlJSerXr588PDzs5V27dtW6devUtGnlo4a3bt2qmJgY+3FiYqIkaezYsVqxYoUkadSoUcrNzdWcOXN0/PhxRUREKDk5WaGhodX3QAAAoM6xGIZRp/fjvLiLxblz59SoUSOzwwEAAJfB53bNwPsAAEDtUJXPbBaBBgAAAAAAqONIAAEAAAAAANRxJIAAAAAAAADqOBJAAAAAAAAAdRwJIAAAAAAAgDqOBBAAAAAAAEAdRwIIAAAAAACgjnMzO4DqZhiGJCkvL8/kSAAAwC+5+Hl98fMb5qD/BABA7VCVvlOdTwDl5+dLklq3bm1yJAAA4Erl5+fLz8/P7DDqLfpPAADULlfSd7IYdfwrNpvNpmPHjsnX11cWi8Wh187Ly1Pr1q11+PBhNWrUyKHXxv/QztWPNnYO2tk5aOfqV51tbBiG8vPzFRQUJBcXZqqbhf5T7UYbOwft7By0c/WjjZ2jutq5Kn2nOj8CyMXFRcHBwdV6j0aNGvEfihPQztWPNnYO2tk5aOfqV11tzMgf89F/qhtoY+egnZ2Ddq5+tLFzVEc7X2nfia/WAAAAAAAA6jgSQAAAAAAAAHUcCaBr4Onpqaefflqenp5mh1Kn0c7VjzZ2DtrZOWjn6kcb41rw76f60cbOQTs7B+1c/Whj56gJ7VznF4EGAAAAAACo7xgBBAAAAAAAUMeRAAIAAAAAAKjjSAABAAAAAADUcSSAAAAAAAAA6jgSQFdpyZIlCgsLk5eXl6KiorRx40azQ6o1kpKS1LNnT/n6+qp58+YaOXKk9u7dW66OYRh65plnFBQUpAYNGmjgwIHatWtXuTrFxcV6+OGHFRAQIB8fHw0fPlxHjhxx5qPUKklJSbJYLJo2bZq9jHZ2jKNHj+ree+9V06ZN5e3trcjISG3bts3+Ou18bcrKyvTEE08oLCxMDRo0UNu2bTVnzhzZbDZ7Hdq46tLS0hQfH6+goCBZLBatXr263OuOatMzZ85o9OjR8vPzk5+fn0aPHq2zZ89W89OhpqL/dPXoPzkffafqQ9+p+tF/qh61vv9koMrefvttw93d3XjttdeM3bt3G4888ojh4+NjfP/992aHVisMHjzYWL58uZGRkWHs3LnTGDp0qBESEmIUFBTY6zz//POGr6+v8e677xrp6enGqFGjjJYtWxp5eXn2OpMmTTJatWplpKamGtu3bzdiYmKM7t27G2VlZWY8Vo22ZcsWo02bNka3bt2MRx55xF5OO1+706dPG6GhocZ9991n/Pe//zUyMzONdevWGQcOHLDXoZ2vzdy5c42mTZsaa9euNTIzM4133nnHaNiwobFw4UJ7Hdq46pKTk43Zs2cb7777riHJWLVqVbnXHdWmt912mxEREWFs2rTJ2LRpkxEREWEMGzbMWY+JGoT+07Wh/+Rc9J2qD30n56D/VD1qe/+JBNBVuOGGG4xJkyaVK+vUqZPx+9//3qSIarecnBxDkrFhwwbDMAzDZrMZLVq0MJ5//nl7naKiIsPPz8949dVXDcMwjLNnzxru7u7G22+/ba9z9OhRw8XFxfjoo4+c+wA1XH5+vtG+fXsjNTXVGDBggL0TQzs7xuOPP2707dv3kq/Tztdu6NChxvjx48uV3X777ca9995rGAZt7Ag/78A4qk13795tSDK+/PJLe53Nmzcbkow9e/ZU81OhpqH/5Fj0n6oPfafqRd/JOeg/Vb/a2H9iClgVlZSUaNu2bYqNjS1XHhsbq02bNpkUVe127tw5SZK/v78kKTMzU9nZ2eXa2NPTUwMGDLC38bZt21RaWlquTlBQkCIiIngffmbKlCkaOnSobrnllnLltLNjvP/++4qOjtavfvUrNW/eXNdff71ee+01++u087Xr27evPvnkE+3bt0+S9PXXX+vzzz/XkCFDJNHG1cFRbbp582b5+fnpxhtvtNe56aab5OfnR7vXM/SfHI/+U/Wh71S96Ds5B/0n56sN/Se3azq7Hjp16pSsVqsCAwPLlQcGBio7O9ukqGovwzCUmJiovn37KiIiQpLs7VhZG3///ff2Oh4eHmrSpEmFOrwP//P2229r+/bt+uqrryq8Rjs7xnfffaelS5cqMTFRs2bN0pYtWzR16lR5enpqzJgxtLMDPP744zp37pw6deokV1dXWa1WzZs3T3fffbck/i1XB0e1aXZ2tpo3b17h+s2bN6fd6xn6T45F/6n60HeqfvSdnIP+k/PVhv4TCaCrZLFYyh0bhlGhDL/soYce0jfffKPPP/+8wmtX08a8D/9z+PBhPfLII/r444/l5eV1yXq087Wx2WyKjo7W/PnzJUnXX3+9du3apaVLl2rMmDH2erTz1Vu5cqXefPNNvfXWW+rSpYt27typadOmKSgoSGPHjrXXo40dzxFtWll92r3+ov/kGPSfqgd9J+eg7+Qc9J/MU5P7T0wBq6KAgAC5urpWyLzl5ORUyPTh8h5++GG9//77+uyzzxQcHGwvb9GihSRdto1btGihkpISnTlz5pJ16rtt27YpJydHUVFRcnNzk5ubmzZs2KBXXnlFbm5u9naina9Ny5YtFR4eXq6sc+fOysrKksS/Z0f43e9+p9///ve666671LVrV40ePVrTp09XUlKSJNq4OjiqTVu0aKETJ05UuP7Jkydp93qG/pPj0H+qPvSdnIO+k3PQf3K+2tB/IgFURR4eHoqKilJqamq58tTUVPXu3dukqGoXwzD00EMP6b333tOnn36qsLCwcq+HhYWpRYsW5dq4pKREGzZssLdxVFSU3N3dy9U5fvy4MjIyeB9+NGjQIKWnp2vnzp32n+joaP3mN7/Rzp071bZtW9rZAfr06VNhG959+/YpNDRUEv+eHeH8+fNycSn/ceXq6mrfxpQ2djxHtWmvXr107tw5bdmyxV7nv//9r86dO0e71zP0n64d/afqR9/JOeg7OQf9J+erFf2na1pCup66uI3p3/72N2P37t3GtGnTDB8fH+PQoUNmh1YrPPjgg4afn5+xfv164/jx4/af8+fP2+s8//zzhp+fn/Hee+8Z6enpxt13313p9nnBwcHGunXrjO3btxs333xzvd6S8Er8dCcLw6CdHWHLli2Gm5ubMW/ePGP//v3GP//5T8Pb29t488037XVo52szduxYo1WrVvZtTN977z0jICDAmDFjhr0ObVx1+fn5xo4dO4wdO3YYkowFCxYYO3bssG/J7ag2ve2224xu3boZmzdvNjZv3mx07dqVbeDrKfpP14b+kznoOzkefSfnoP9UPWp7/4kE0FVavHixERoaanh4eBg9evSwb8GJXyap0p/ly5fb69hsNuPpp582WrRoYXh6ehr9+/c30tPTy13nwoULxkMPPWT4+/sbDRo0MIYNG2ZkZWU5+Wlql593Ymhnx/jggw+MiIgIw9PT0+jUqZOxbNmycq/TztcmLy/PeOSRR4yQkBDDy8vLaNu2rTF79myjuLjYXoc2rrrPPvus0t/FY8eONQzDcW2am5tr/OY3vzF8fX0NX19f4ze/+Y1x5swZJz0lahr6T1eP/pM56DtVD/pO1Y/+U/Wo7f0ni2EYxrWNIQIAAAAAAEBNxhpAAAAAAAAAdRwJIAAAAAAAgDqOBBAAAAAAAEAdRwIIAAAAAACgjiMBBAAAAAAAUMeRAAIAAAAAAKjjSAABAAAAAADUcSSAAAAAAAAA6jgSQAAAAAAAAHUcCSAAAAAAAIA6jgQQAAAAAABAHUcCCAAAAAAAoI77f0sjkE9Y3VndAAAAAElFTkSuQmCC", + "text/plain": [ + "<Figure size 1400x400 with 2 Axes>" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "plot_convergences(r2, np.ones(4))" + ] }, { "cell_type": "code", @@ -320,10 +424,27 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 81, "metadata": {}, - "outputs": [], - "source": [] + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[np.float64(0.9960440949231335), np.float64(1.0021427377274625), np.float64(1.1881434006902023), np.float64(0.8133186744545533)]\n", + "[np.float64(0.9958710879907714), np.float64(1.0022285003471707), np.float64(1.1883127324043288), np.float64(0.8132176353399935)]\n", + "[np.float64(0.9951033518959111), np.float64(1.0026010928526303), np.float64(1.1891375637452173), np.float64(0.8127002814019945)]\n" + ] + } + ], + "source": [ + "r1 = gauss_seidel_r(A, b, x0, w=0.01, n=1000)\n", + "print(r1[-1])\n", + "r1 = gauss_seidel_r(A, b, x0, w=0.0099, n=1000)\n", + "print(r1[-1])\n", + "r1 = gauss_seidel_r(A, b, x0, w=0.0095, n=1000)\n", + "print(r1[-1])" + ] }, { "cell_type": "code", @@ -332,5 +453,26 @@ "outputs": [], "source": [] } - ] -}
\ No newline at end of file + ], + "metadata": { + "kernelspec": { + "display_name": "Python 3", + "language": "python", + "name": "python3" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 3 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython3", + "version": "3.13.3" + } + }, + "nbformat": 4, + "nbformat_minor": 2 +} |
