{ "cells": [ { "cell_type": "markdown", "id": "collected-allah", "metadata": {}, "source": [ "# Simulation and illustration of several results of \"Points and lines configurations for perpendicular bisectors of convex cyclic polygons''\n", "\n", "In this notebook we illustrate several of our results, by means of simulation and explicit computation. We follow the subsections, references and order of the revised version of the paper." ] }, { "cell_type": "code", "execution_count": 1, "id": "exceptional-substitute", "metadata": {}, "outputs": [], "source": [ "# We make the simulations repeatable by setting the seed\n", "set_random_seed(0)" ] }, { "cell_type": "markdown", "id": "bright-serbia", "metadata": {}, "source": [ "---\n", "\n", "## 1.1 Deterministic results\n", "We begin with **Theorem 1**, which states that a word is realizable *iff* its signature is interlacing." ] }, { "cell_type": "code", "execution_count": 2, "id": "diagnostic-possibility", "metadata": {}, "outputs": [], "source": [ "def sim_data(n):\n", " # simulates a configuration for n iid points uniform on the circle\n", " # the data is returned as a sorted list of elements of type (pos,index)\n", " # where 'pos' is the position in [0,1] in increasing order, and 'index' is either:\n", " # - 0 for a point\n", " # - 1 for the diametrically opposite point\n", " # - 2 for a perpendicular bisector\n", " # - 3 for the diametrically opposite part of the bisector\n", " P = [random() for _ in range(n)] #n iid points\n", " PP = [(P[i] + 1/2) % 1 for i in range(n)] #diametrically opposite points\n", " P.sort()\n", " L = [ (P[i] + P[i+1]) / 2 for i in range(n-1)] + [ (P[n-1]-1+P[0])/2 %1 ] #bisectors\n", " LL = [ (L[i] + 1/2) % 1 for i in range(n)] #diametrically opposite side of bisectors\n", " # the next list, U will contain, the data described before\n", " U = [(x,0) for x in P] + [(x,1) for x in PP] + [(x,2) for x in L] + [(x,3) for x in LL]\n", " U.sort() #default sorting is w.r.t. the first entries\n", " return U" ] }, { "cell_type": "code", "execution_count": 3, "id": "impressed-dominican", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[0, 1, 0, 1, 1, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 0, 1, 0, 1, 0]" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "def sim_word(n):\n", " # simulates data and deduces the occupancy word\n", " U = sim_data(n)\n", " V = [] # will be the occupancy word\n", " current_number = 0\n", " for (_,x) in U:\n", " if x==2 or x==3:\n", " V.append(current_number)\n", " current_number=0\n", " if x==0:\n", " current_number+=1\n", " V[0] += current_number\n", " return V\n", "\n", "# an occupancy word for n=10 random points:\n", "sim_word(10)" ] }, { "cell_type": "code", "execution_count": 4, "id": "included-governor", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[1, 0, 2, 1, 0, 1, 1, 1, 2, 1]" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "def word_to_sig(V):\n", " # takes a word and returns its signature\n", " n = len(V)/2\n", " return [V[i] + V[i+n] for i in range(n)]\n", "def sim_sig(n):\n", " # simulates the signature of the occupancy word\n", " V = sim_word(n)\n", " return word_to_sig(V)\n", "\n", "# a signature of a word for n=10 random points:\n", "sim_sig(10)" ] }, { "cell_type": "code", "execution_count": 5, "id": "healthy-dispute", "metadata": {}, "outputs": [], "source": [ "# We now test if signatures are always interlacing.\n", "def is_interlacing(S):\n", " # Takes a signature S (list of 0, 1 and 2) and says if it is interlacing\n", " n0 = S.count(0)\n", " n2 = S.count(2)\n", " # If there are no 0, no 2, or not the same number of 0 and 2, fail:\n", " if n0==0 or n0!=n2:\n", " return False\n", " # Otherwise, we run through S and each time we see a 0/2, we check if the previous was different.\n", " # Initially, 'previous' is the one *different* from the first 0/2 in S.\n", " i0 = S.index(0) #first index for 0\n", " i2 = S.index(2) #first index for 2\n", " if i0= 14 or so.\n", "n = 13\n", "\n", "l = all_words(n)\n", "v = l[ int( len(l)*random() ) ] #one word taken uniformly at random\n", "s = word_to_sig(v)\n", "\n", "# Lists for the three functions of 0s, 1s, 2s in the signature at each position\n", "l0 = []\n", "l1 = []\n", "l2 = []\n", "for i in range(n+1):\n", " l0.append( (s[0:i].count(0) - float(i/6)) * 2 / sqrt(n) )\n", " l1.append( (s[0:i].count(1) - float(2*i/3)) * 2 / sqrt(n) )\n", " l2.append( (s[0:i].count(2) - float(i/6)) * 2 / sqrt(n) )" ] }, { "cell_type": "code", "execution_count": 32, "id": "interim-profile", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA2YAAAGFCAYAAACWmuD7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAA9hAAAPYQGoP6dpAACc0ElEQVR4nOzdd3gUxRvA8e+k0RN6E6QoLfQSaQroD0FEpBdFBUVQgShFERQVKypFBIKogIiCHL0JAkpHiiA1AQIIgtJbQg0kmd8fc8EQ0nN3e5e8n+e5B7K3t/vebX1nZmeU1hohhBBCCCGEENbxsjoAIYQQQgghhMjqJDETQgghhBBCCItJYiaEEEIIIYQQFpPETAghhBBCCCEsJomZEEIIIYQQQlhMEjMhhBBCCCGEsJgkZkIIIYQQQghhMUnMhBBCCCGEEMJikpgJIYQQQgghhMWyRGKmlGqilNJKqbz2v7srpS45eZ1TlVILnLmO9HLF908rpVQbpdQhpVSMUmqMi9dd2r5/1HDletNKKVVUKbVSKXXVyu2nlDqqlOpn1fozKuH5ICtyxDa04rzqKdxlH1NK5VRKzVVKRbpDPK7mKed2R1JKDVNK7XTCcj32t0zs3KSU6qWUOq6UirXiemb/Ldu4er3x1r8mpXstF90rp3gtctY+7a7SlJjZkw2tlBqcYHobpZR2bGhOZQPKWx2EuMPXwBygJPCOi9d9HCgG7HXxetOqPybOGnjY/puai4Cn8fAENQj4xsHLtPS8qoxhSqkTSqnr9n2uslXxuIluwENAA8y5I8LacITV3KXQwMXuODcppfyB8cBnwD04/lzoCdoR717LwuvZHdciqxPW5Cil2imlliulziVVSKGU+lopddh+DTqrlFqolKqYlvWkp8bsBvCmUipfOj6bJKWUnyOXlxyt9XWt9RlXrS+rUEr5pvNzuYHCwHKt9Qmt9WVnri8hrXWM1vqU1jraEctzovuA7Vrrg7L/iozQWp/VWl9z8DKtPq8OAgYAfTEX+1PASqVUHgtjstp9wD6t9V77OS7NBahKKW+llFNa19iTaR9nLDuzc9T1LytI5Nx0L+AL/Ky1Ppnec6EnbwOt9YXU3ms5OQ6HX4ucKBewERiczDzbgeeBSkBzQAErlFLeqV6L1jrVL2AqsBjYB3web3obs6g75m0PhAJRwFFgYIL3jwJD7cuMAL4HugOXgCeAA8A1TC1KLkzJ31HgIjAO8I63rGeAbcBlzMV4BlA43vtNAA3ktf/dHbiUIBad8BXv/XswJS4XgfPAQqB0vPe9gdH22M8Dn9u/z4Jkfsu479rc/nteAX4BisWbZw0wJsHnFgBTE/kdp9mX8TfQGihkj/MKsAeok8i62wDhmGR7JVAywbpaYXayG8BfwHuAT7z3NfCyfT1XgfeT+K757PFdtG/TZUC5BNsm/qtJEstJdH2piLMisMH+fhjQ1L6sNvb3S9v/rhHvM42BrZj99yTwaYJlrgHG2rf1Bcx+NyxBvMOAY/ZlnADGpnB8vQIcBm5i9v9nk9lHpyaznOft+9QNYD/QO8H7n9m3+zX77/Uh4Jtgnicxx9QN4BwwL0EsbwFTMMfcMaBXCueNhNu4DHAIeD3BvFWAWOC+eNv8Ffs+cx04AnRM8Jlkj89E4onb51oCu+zfcQtQNcF8DYB19vUet2/vXPG2f8LvpICzQPt4y9gJnIn3d33gFpDb/ncAprTwDBAJrAKqp+M4fBGYb9+mB4EnU9jXjgL90rIM4HH7fnMdWI05jyR5Xk3FfuSHOX7+xRzPW4h37AOlMNebi/b3Q4HHk/g+CnOcvhlvWjbMee6lZH6HIMy57xzmOrQWqJXIeSdDv00S685r3/an7b/PXuCJ+L8lyVwfUnpx9z66xj49yfNxgnU/gTlfRgNlElm+NzAZc0xex5yzXkshpib2WJrb94ubwMOYBHKh/be4AvwBNE1kn032vAM8AOyw/57bgLak79w+Dhhj/41OA70w9yHf2dd9GGiRzPcMBvbE+7uNPY4+8aYtB4bH+zvJ83+8/TCx699ge4yX7dvjU2BnEnGV5u7z1lT7e49hrpOXMOfRJdjPwwk+W8P+txfwLWa/L+XEc9Xta3W8aZeA7gniaoc59q5hzuv1E+7T8f6f8Dcond5tgLnO7wRewOyTV4CvMMfHIMy9wRng7aS+F+a8Pz7B+wUw++gjqTze5wLj4v09xr6Oyva/fez7SPN4+/mYJM4VOr3nIfv2Hxjv7wWYc4i//e+i9nVUiHdc94v3//hxHLVPj/uNn7XPEwHMBPIkE0eaY0/ti0TuGZOZt5p93vtSvfw0BjPV/iO3xZyIS8Q/6cSbrzYQg6kmLW//ga5hP5DibYAI4HXgfvurO+aAWAHUBBphLpjLMTdegZiLRRTQOd6yXgBaAGWBesAmYGm895uQfGJWyL6zFMXc5G0C1tnfy4k58UwGqmKy4OmYG14/+zyD7DtAe/v7kzA3WQtS2GluYm4K6gC1MBfB6fHmWUPqErPzwEtAOWCCfd3LgI7233++fdkqwbr/wNwo1sbcFG2Mt9zm9u3Tzf67Poq5AL+X4MRy2v77l8V+ck7kuy60r/8hoDrm4DiIKbHys8cYd2ItGve7JrKcu9aXUpyYi8d+zD5VHXjQ/l3jnxRLc+cF5x7MSTcEk9S1wdxwD0uwbSIwF55ywHOYhOJR+/sd7O+3wJTOPQD0TGZ/aGvfJr3tv8cAzMns4Xj76DLMcVAUCEhiOT0xSWA7TPLTzr5/dIs3z1BM0lEacyE9BQyK935L+7rfx+zP1YG3EtnnemOO28GY471iEjEFAL9jbkLjjjNvzE1WaIJ5RwNrE2zzc5iLeXlMEhkNVErt8ZlIPE3syw2z7y9VMQnAEewJqn3aZaCfffs2AP4EvrO/nx+TrL0T953s029fHDE3wDcxN3eB9mlDgM32/yvMjdAizDmgHDDS/n3zp/E4PA48Zd8eX9pjz5/M/naUuxOzJJeBaWJ8A3OxrwB0te83yZ1XU9qPpmNKHx/C3Ji/bl9HXKHNEsxxW9X+3Z8AGiXxfcraY6mZyLnn+2R+h0cwBXuV+O/cfYp4F3xH/DaJrNcLc53Za9+mcd+vRWqvDym9MPvoN5hjr2i8eJM8HydY90bMfl8Be4FEguX72rdtEOZc0xVz3uyUTExN7L/LLvv3vg9zA1odcw2rijkOPsLcY9yb2vMOJnE6g7lhq2z/PQ+TvnN7JOY8Wc7+bwywFHN+jbvOngNyJvE9q2KuBwXtf39hX88s+99xN8mPpeb8H28/THj964S5H3rRvp0+sse+M4m4vDHXBG1fz+1rCeb+pZ39+9XAnJd2A17290vH/ZaY6/YcTBJc2MnnqtvX6njTLnF3YrYPc84pD8zG7C8+Cc9NQA7gf/bPBPHf9Si922CY/TvMxtyjtrJvk18whXkVMIWlGqiX2PcCnsYU8GaL9/6r9t9PpfJ4T1gYsAOzz/W2/52wUHAN/yVmSV3PupPG8xAwClhs/7/CHLNnsReq2bf9ycSuRZj7HG1fb1GgkH163G88F1N4+xCmUOXjZOJIMXbMOetKCq+uiSy7NKlIzDDnpC8whRSJ3o8k+rnUzmhfyVTsyQbmojLZ/v823JmYTQdWJPjs58S7CbNvjPmJ/JCaO0tpJmJOpLnjTfsFmJhMnEH25cTtgE1IoWQ33me/tMcWt0O8gLnJU/Hm8cMkms3sf5/gzpJaH8xOviCFnSbhd+0NnIr39xpSl5j9EO/vuNKID+JNq2efVjTBuuvGm6eifdoD9r/XAUMSrPsZ4ES8vzXwRQr7TDn7fA3iTStg//062v/OSzI1ZcmtL6U4MaWAt+K+u31asjVmwMeJbPPemBND3EVqDbA+wXq3Ap/a/z8AU+Lmm9x3ivfZjcA3CabNwjS1SHTbJ7GcY8BTCaYNBX5P5jNvANvi/f078GMy8yfc5xTmYvVyMp9JbF8uhrnwxe1zvpgbq24JtvlXCT63GZiQ2uMzkVia2Jcbv3Anv/0znex/TwO+TvC5BzE3aNnj/Q79Esxz++KIqbn+A3Mxibs4Lo+3jzyCuZHJlmAZh7DXBJD64/DDeH/nwtwUPpbCNuyX2mUAnxCvcMc+7VOST8yS3I8wN+SxQPEE038FPrH/fzfxbupS2O8b2GNJuLxvME2kU1yGfX5vzI3tE478bRJZTzP7vlQ+ife7k8L1IZXfZwz2mjL736k5H8etu3pa1mX/bAgwJ5n3m9iX3ToVywoF+ibYZ5M872Bqtc4TL1nC1G5k6Nxu3yeuANPiTYu7ztZLIvY7as8xN8mDgdP2vxPeJKfm/J/Y9e93Ej8/7kzFNkh034w3X9xNchX736Xtfz+IudndQLwCQpx3rrp9rY437RJ3J2Y94r0faJ8Wl7R3585zUw3i1ZRlcBsMw9yjxi/M+QWTVHnFm7YfGJzY98LU7p8nXqGGfZ95Lw3H3u3CAP4rFHyb/woDbhcKxtvPxyQ4vvolWGZ30ngewiSmlzCFT9Ux1/TR2FvZYfoTmJnUepPY3on9xp/H/z6JxJFi7EAe/qsYSup1V60cKSRm9vVc4b8Cg1TXlmmtM9Qr45tAN6VUYCLvVcLs5PFtBMolaGe5LZHPXtNaH47392lMdeaVBNMKx/2hlKppf8Dub6XUZcwOB6amItWUUr2AHpiLxln75NqYjXNZKXVFKXUFU7KRHbhPKRWAucHcFLccbZ5VSuy7JZTwu56M/73SYHe8/5+2/7snkWnxl31HjFrr/ZiDqZJ9Um3g3bjvbP/e3wLFlFI54y0npe9Zyb6uLfHWdR6TtFRK6kPJSLi+lOKsABzXWp+K95mtqYh5k7YfYXYbgdxAiXjTdnOn+NtvNqZk7i+l1LdKqbYpPEuR1DGT6t9IKVUIU3o/OcHvMRRzIxw3Xwel1Aal1Cn7+x9y57FSA/gthdXd/u723+kUadx3tdYngZ8xyRWYEu7smN8uvk2J/B1/P03y+EwhhPjH7AXu3CdrA90T/I7LMRebMskscw1QWSlVENNkao391di+/RtgmsvFrSM3cD7BesrEiz21x2H87XEVc6OZ1nNJcsuohLkIxj8mEm6XhGqQ9H5UC3PzGp7guzXmv+8+FhiqlNqolHpfKVUtFd9BJ/hbJTLtvzeVKqyUmqiUCldKRWAS5dzcfe1wxm/zj9Y6PJl5HHV9iC+15+Ob3H1+u4tS6mWl1Db7Q+5XMDVKqbnu3nEeV0rlUkp9rpQKU0pdsi+rYiLLSu68UwnYpe98XiXhdkjzuV1rHYO5aU7pmnqbffnrgCb2TjYqYwqZvZVSlTDJ0Z/x7mtSe/5PeP2rROLnxzRTSt2nlJqhlPpLKRWJSSzg7m3wE+b3aqa1jt+ZjCvPVYmJv7+etP+bluWmdxuAuUeN/7zWaSBMax2bYFpS+0sU8CP2a6G9U4nqmMqQ1NqL2U8bY2qUdmFqPRvb32/Cf9eetEjreWgdJuGpaV/3WkwT04zGkfA3Ts35MNnYtdaXtdaHUnil5zm86fz3/Q8Cs5RS2VP74XQ/dKu1XqeUWo4pLZya4O3ELoYqkcVcTWTarYSrSmKaF5gTOqa5ywpM6cxZzIlkOabkPFWUUk0w7cqf0lrviveWF6bNbNdEPnY2kWlpkdj3iv87xXL375bYw6a3l6O11kqphMuO2xYJE/HEbljiz/seMC+ReW7E+39i2zC+xLZ73PQkb5iSkXB9KcWZnvUkt//Gn57kfqm1Pq6UqoBpztEU0/TlDaVUY611ws/F/3xKcSQnbvv2JN6Nl10MgFKqHqapz3uYYyQC6AIMjDfv9VSsK8nvnkaTgB+UUv0xzT1sOnUPAsffTx15fMZf7teY5CChY8l8Pv7FsTHwLqb2/G1MTX4OTElz3DpOYi5UCV2KN09qjkNHbI/klpHUcZyc5PYjL8w+GdfsPb4rAFrrSfZrTEtMDdMQpdRArfW4RJYXV/BSlP9uysBchE/fPfttUzE1A/0wz+dGYW5sE147XPnbJLfO9KwrvtSej68nSF7u/oBSnTDNdAZifrPLmNr3uqmII+F5fASmKdzrmBrj65imco7eDhk5tyd2nU3uGFuDqcV7CJMwXlJKrcOcF5rwXwFy/HWkFGtK19uMWIw5V8U1h/fCnM8SboOlmHutephno+I461yV2H6f7H0QSd/zpCS92yBN961JmATsVEqVwCRov2mt/04+3HgLN/vkOsy+dROzf+3FFAZUxRQKjknt8uJJ03lIax2hTNf2TezrXAWsB2oopcphmomucVAcKW3fZGNXSnXFXOeT85LWenqqIoxbiSmwiAAOKqU2Yx5naIsp1EhRRntDGox5IC9hqV8Ypro7vgZAuL30yZEqYqpuB2utjwMopeqkZQFKqfsxzY0+0VonPKn8CXTGPMAfmcTnT2JOUuvsf/tgbjj+TEsciTiLqY2LW483pn3t6gwuF8y2r4O99sieROTFVLeDib2C1vpQBtcTZl9XXUyzC5RSBTAH574MLhtSiFMptR+4VylVRGsdd4MWlIqY2yulVLybkwaYG49/UxuY1vo6psRqkVIqBPPbViXx/WIf5piZFm9aA9LwG2mtTyul/gXKJnMiaQj8rbX+OG6CUqpUgnl2Y9rgf5fadafCTUyToISWYi50r2Cex2uUyDz1uPN3qYdp5gGpOD6TUQ97kqVML7PluXP/r5zC/n/Xd4p3cWyNOVbXY/YbX0yzqj/jlcD9iUkkorXWR5NYh6OOw4wKwzRZj69eCp9Jbj/agfntCmut1ye1APs5fSIwUSk1HHPTmFhidgSTnD1qX3ZcT7+NMa07kvIQppnpUvtnSmKuJ2mR3t+mhFKqfAq1Zo7myPPxQ5gm0hPiJiilUqqlTm5ZU7XW8+3LyY1pLpQWYcCzSqkc9nMv3L0dHHJuT6U1mEcjOvDfjehaTEFdA/t7cdJ7/t9H4ufH5Ny0/3v73GXfByphbkLX26clvIeL8xXmhn+RUqql1jqu9sNZ56qE90HlMM8WO1qGr8EZobXeo5TahjnHPY1pFp9WazCFATeBd+3Xo/WYAo8c3F0jGF9S1+j0WIPp1KeuPY5LSqkwTOudM1rr5H7TWw6MIyWLuLsQO6HkCvZSS2Gaq6ZKhhIz+440nbt3oFHAH0qpdzCdFdTHdGHcOyPrS8IxzA4VrJSaiLkZSvU4WEqpHJiSop3AN0qponHv2Zu/TceUAi5USr0L/IOpkWsHjNBa/4M5wQ5WSh3EHMQDMElORq0CRiulWmIeYu7voOWC2fnHKaVetf9/PKY5Tlwzvw+AJUqp45imZbGY3mWqaq2HpnYlWuuDSqmFwLdKqZcwF8BPMRfBhQ74HinFuRLz232vlBqEqWKPS0qSKhGegClBH6eUGo9pDvk+MDpB04QkKaW6Y04uWzDPbzyLKQVOqgRsBKa6+09M869WmH2saWrWF88wYKy9KcoyzMmgDpBPaz0aUxp9r1KqC+b5p5aYkpz43gd+U0odxtSu+WA6Jfg8jbHEdxSoq5QqjakRuaC1jtVaxyilpgLDgUNa68Sa4XS0X7A2YGrGHsA0OYbUHZ9JeVcpdR5z4v0Y8zD/Avt7nwGb7Qn1t5jksRKmc5e4891RoJFSaiYQpbU+Z5++BlOTsCMuWbQna10xbe3j/IqpaViglHoT05ysOKaHvwVa62046Dh0gInAQKXUaEwJY21MG/7kJLkfaa3D7deOaUqpgZhkqiDmubs9Wuulyox7twxT8JfP/l6iF3T7DcgY4C37efggpnOZa5heepNyCHMzvw3wxxyHqanNii/Nv43Weq19n5irlBpgj6Oi/av8ksb1p5qDz8eHgOeUUs0xifGzmEKvI8l+KulltVNKLcb+HBJpr+2YgTmOJyulPsIkdq8nmCfD5/Y0iKs974opqAFzbhhl//+GePOm9/z/JebaFv/8WBnT2UBS/sb8xk8opZZi9veL9lh72Qua78XsF4nSWo+zFxQvUUq10FpvwHnnqlVAX3vNgxfm3JxUq5OMcNQ1OCMmYe7FrmE6bUurNZh9IhpTKBg3bRSmUDC5wsujJH49S481wGuYxwrC4k0LJvEa1YRx/E8ptdEex8UMxJEseyFpqpsqKqXyY46N4vZJFew156e01qeUUmUxBcUrMAUK92AKBq9jCqFTJSPPmMV5hwTVmlrrPzG9BXXBnJw+wGTNUx2wvjto8yxYd0wPhGGYWryEJ+PkFMFcEB/BVN+fjPfC3qyqESYBnIe5MZiCKX2I28lHYUpZpvJfk470HFQJTcF0uz8NU9J2BMfUloE58D/DXMw2YXacLnFvaq2XY575eRRzA78Zk3Cmumo9nucxzc2W2NelMD30ZPjkmlKc9hraNpg28X9gTnwf2T9+g0Rorf/F3Bw/gGmnPRHT699Hic2fhEuYkq+N/Fdz0Eqb5zkSW+cCzInsDcxD7y8Bz2ut16RhnWitJ2F66OqOeSZirf3/R+zvL8QkDeMxhRENMDdB8ZexBnM8PWmfZxWpa56UnJGYJmth/NfcOM5kTHOZKUl89j3Mvrkb0+NXV611mD3W1ByfSRmMuYhtx5TIPqm1vmlf7m5MbUs5zAVuB+Z3it9M7l3Mzd9h7mw2uRqTlK+JN22tfdrttvX2EvvHMTXtUzAJyEz7Mk/b53HkcZhuWutjmF7bWmGOiZcxiU9yn1lD8vvR85hz2yhMUrrI/v5x+/vemM4k9mEepj9A8oV7n2Oa6kzAPAdyD+Y5mOQuvC9gkr4dwA+YpqtpGostPb+NXXvMNv0Jc1x8ThpKipW6PVBw6bTEi+POxxMxx5wNUwBVAPPbp0d/THLwO6agdDlpbHGizfNarTAdP+zAJGlvJpjHEef21Maj+e94j7tJ3o1p4nS70MY+7wLScf7XWtsw91efYbZpKUyNVnKf+RdzTv0Uc54Zb09Ku2AKFfZirhFvpLCcMfblLFVKNXDiuWog5pywDnO/MhJz/+JQjroGZ9BPmKRqhtb6jvsTpdQwpdTRFD6/F1PAuCve/nXXtScJSV3P0mNd3Lrj1UynNo6BmH3oOP+1jHEXT2Ji+tn+90z73y/b/76Bqf1fiilsmoUp1G2g0zDGp9LJNyMXItNRSjXElC7er+98MFRYwL491mCG3zid4D0NtLVfNIUQdvZa+bcxQzE4owZBCOFC9qbUR4EgewVH/PemAmitu7s8MOFSGX3GTAi3p5Rqi2k+d5D/xk7ZKEmZtZRS2TC9SH6I6dLXEW25hcgqHsOMCydJmRAeTCnli2m18SnmkZLEaosbk/gz2CKTkcRMZAV5MM2ESmKq+X/lzl4IhTWewjQj2ol5PkUIkUpa6y4pzyWE8AANMU3gwzGdxdxFa53cMC0iE5GmjEIIIYQQQghhMUd0/iGEEEIIIYQQIgMkMRNCCCGEEEIIi0liJrIEZfgr+6ATQgghhBBCuBPp/ENkBik+KBkREUFAQAARERGuiEcIIYQQ1pACWOGxpMZMCCGEEEIIISwmiZkQQgghhBBCWEwSMyGEEEIIIYSwmCRmQgghhBBCCGExScyEEEIIIYQQwmLSK6PI1EJCQggJCSEmJsbqUISwxG9//YaX8uLhMg9bHYoQIhO4efMm0dHRlsbg4+ODn5+fpTEI4QxK6xR7GhfC3aW4E0dGRt7uLt/f398VMQlhuVsxtyg6qigXrl+gbcW2jGo2ijL5ylgdlhDCQ928eZPQ0FBiY2MtjcPLy4vKlSsnlZxJd/nCY0mNmRBCZFK//vUrF65f4MOHP2TitolUCqnEoIaDGPzgYHL65rQ6PCGEh4mOjiY2NpbSpUuTI0cOS2K4fv06R48eJTo6WmrNRKYjiZkQQmRStlAbFQpU4O2H3qZ/vf4M3zCczzZ+xtSdUxnVbBQdAjuglBQuCyHSJkeOHOTMaW3hTkREBDdu3CBbtmzkypXL0liEcBTp/EMIITKhG9E3mL9/Pp0rd0YpRS6/XHz0yEeE9Q6jZrGadJrTiUemPcKe03usDlUIIdJsyZIlTJs2DZvNxtWrV60ORwiHkMRMCCEyoeWHlhMZFUnnKp3vmH5f/vtY2GUhy7ou48TlE9T8uiavLnuVi9cvWhSpEEKkXa5cufDz8yMyMpKoqCirwxHCISQxE0KITMgWaqNq4aoEFgpM9P3H7n+MPa/s4dOmn/Ldzu8oN64c32z/hphY6cFUCOH+smXLRrZs2awOQwiHksRMCCEymWu3rrHowCI6V+6c7Hx+3n683uB1wvuG07J8S15a8hIPTHqA34//7qJIhRBCCBFHEjMhhMhklh5cytVbV+9qxpiUYnmK8X2b7/n9hd9RKBpOaciz85/lxOUTTo5UCCGEEHEkMROWUEr1VkodUUrdUEptV0o9lML82ZRSHyul/lZKRSmlDiulXnBVvEJ4kpl7Z1KrWC3uz39/mj5Xv2R9tvbcyqRWk1h+aDkVxlfg842fExUtz28IIYQQziaJmXA5pVRnYAzwMVATWA8sU0rdm8zHZgH/A3oAFYCngP3OjVQIz3M56jI/H/yZLpW7pOvzXsqLHrV6EB4cTo+aPXjrt7eo+lVVlh1c5uBIhRBCCBGfJGbCCgOAyVrrSVrrfVrrfsBx4JXEZlZKPQY0Bh7XWv+qtT6qtd6qtZYHYYRIYHH4Ym5E36BT5U4ZWk7e7HkZ89gYdr28i5IBJXl8xuO0+qkVhy4cclCkQgghhIhPEjPhUkopP6A2sCLBWyuABkl87ElgGzBIKfWvUipcKTVSKZUjqfVERUURGRl5x0uIrMAWaqNeiXqUylvKIcurXLgyvz77K3M6zmH36d1UnlCZIb8O4crNKw5ZvhAi65owYQJlypQhe/bs1K5dm/Xr11sdkhCWksRMuFpBwBs4nWD6aaBoEp8pCzwIVAHaAv2ADkBIUisZPnw4AQEBt18lS5bMaNxCuL1LNy7xy6FfUuyNMa2UUrQPbM++PvsY8uAQxmwZQ4XxFZixZwZaa4euSwiRNdhsNvr168fbb7/Njh07eOihh2jRogXHjh2zOjQhLCOJmbBKwrs5lci0OF7297ramzAuxTSH7J5UrdmQIUOIiIi4/Tp+/LjDAhfCXS3Yv4BbMbfoGNjRKcvP6ZuTYU2Gsa/PPuqXqE/XeV156LuH2HFyh1PWJ4TIvEaPHk2PHj148cUXqVSpEmPGjKFkyZJ89dVXVocmhGV8rA5AZDnngBjurh0rzN21aHFOAv9qrSPiTduHSeZKJPYBGXhSZEW2UBsPlXqIe/zvcep6SuctzZxOc/j1r1957ZfXqP1NbXrV7sVHj3xEwZwFnbpuIYSbuXYN9qfQF1fFipAz5+0/b968yfbt2xk8ePAdszVr1ozff5fHx0XWJTVmwqW01jeB7cCjCd56FEjqbLwRKK6Uyh1vWnkgFvjH4UEK4YHOXzvPr3/96vBmjMlpWrYpO1/ayRfNv2Dm3pmUH1eekK0hRMdGuywGIYTF9u+H2rWTfyVI3M6dO0dMTAxFihS5Y3qRIkU4deqUK6MXwq1IjZmwwmjgB6XUNmAT0Au4F5gIoJQaDtyjtX7OPv8M4B3gO6XUe5jn1EYAU7TW110dvBDuaN6+ecTqWNpXau/S9fp6+/Javdd4qupTvPXbWwQvC+abP79h7GNjaVy6sUtjEUJYoGJF2L495XkSoZS642+t9V3ThMhKJDETLqe1timlCgDvAsWAvZiu8P+2z1IMk6jFzX9FKfUoMA7TO+N5zLhmQ10auBBuzBZq4+HSD1Mkd5GUZ3aCwrkKM+nJSbxU+yWClwXT5PsmdK7cmRGPjqBkgHS+I0SmlTMn1KqVpo8ULFgQb2/vu2rHzpw5c1ctmhBZiTRlFJbQWk/QWpfWWmfTWtfWWq+L9153rXWTBPPv11o/qrXOqbUuqbUeKLVlQhinr5xm9dHVLm3GmJSge4L4vcfvTG09lTVH11AxpCIfr/uYG9E3rA5NCOEm/Pz8qF27NitXrrxj+sqVK2nQIKmRc4TI/CQxE0IIDzcnbA5eyot2ldpZHQoAXsqLbjW6ER4cTu86vRm2dhiBIYEs3L9QutcXQgAwYMAAJk2axJQpU9i3bx/9+/fn2LFjvPzyy1aHJoRlJDETQggPZwu18WjZRymQs4DVodzBP5s/I5qNYM8reyhXoBxtbG1oMb0F+8+l0IObECLT69y5M2PGjOGDDz6gRo0arFu3jqVLl1KqVCmrQxPCMpKYCSGEB/s38l82HNvgFs0Yk1KxYEV+6foLC7ss5OCFg1T9qiqvr3idyKhIq0MTQliod+/eHD16lKioKLZv306jRo2sDkkIS0liJoQQHmx22Gx8vX1pU7GN1aEkSynFkxWeJLR3KMMaD+OrbV9Rflx5vt/5PbE61urwhBBCCMtJYiaEEB5s5t6ZtLi/BQHZA6wOJVWy+2Tn7UZvs7/PfpqUbkL3hd1pMLkBf/z7h9WhCSGEEJaSxExkaiEhIQQGBhIUFGR1KEI43NFLR9ny7xa3bsaYlJIBJZnZYSZruq3h2q1r1J1UlxcXvciZq2esDk0IIYSwhCRmIlPr06cPYWFh/PGHlMaLzGdW6Cxy+OSgVYVWVoeSbo1LN+bPl/5kXItxzNs3j/LjyjNm8xhuxdyyOjQhhJOsW7eOVq1aUbx4cZRSLFiwwOqQhHALkpgJIYSHsoXaaFm+Jbn9clsdSob4ePnQ54E+hAeH81SVpxiwfAA1vq7Bb3/9ZnVoQggnuHr1KtWrV2f8+PFWhyKEW5HETAghPNDB8wf58+SfHtmMMSkFcxbkqye+Ynuv7eTPkZ+mPzSl/az2HL101OrQhBAO1KJFCz766CPatXOPsReFcBc+VgcghBAi7WyhNnL75ebxco9bHYrD1SxWk3Xd1/HT3p94Y+UbVAqpxJsN32RQw0Hk9M1pdXhCiHiuXYP9KQxNWLEi5JRDV4gUSWImhBAeyBZq48kKT2baREUpxdNVn+bJCk/yyfpPGL5hON/t/I7RzUbTrlI7lFJWhyiEwCRltWsnP8/27VCrlmviEcKTSWImhBAeJuxsGHvP7OXjRz62OhSny+2Xm0/+9wkv1HyB/sv702F2Bx4p8whjHxtL5cKVrQ5PiCyvYkWTeKU0jxAiZZKYCSGEh7HttRGQLYDm9zW3OhSXuT///Sx+ajFLDy6l3y/9qD6xOn2C+vD+w++TN3teq8MTIsvKmVNqw4RwFOn8QwghPIjWmpmhM2lTsQ3ZfLJZHY7LPV7ucfb23ssn//uEKTunUG5cOSb9OYmY2BirQxNCCCEyRBIzIYTwILtO7yL8fDhdqnSxOhTL+Hn7MajhIA70PUCL+1vQc3FP6k6qy6bjm6wOTQiRCleuXGHnzp3s3LkTgCNHjrBz506OHTtmbWBCWEwSMyGE8CC2vTYK5CjA/8r8z+pQLFc8T3GmtZ3Ghuc3oNE0mNKAbgu6cfLySatDE0IkY9u2bdSsWZOaNWsCMGDAAGrWrMm7775rcWRCWEsSMyGE8BBaa2yhNtpVaoevt6/V4biNhvc2ZOuLW/nmiW9YenApFcZXYOTvI7kZc9Pq0IQQiWjSpAla67teU6dOtTo0ISwliZnI1EJCQggMDCQoKMjqUITIsG0ntnHk0pFMNai0o3h7edOzdk/C+4bTvUZ3Bv86mKpfVeWXQ79YHZoQQgiRKpKYiUytT58+hIWF8ccff1gdihAZNnPvTArnKkzj0o2tDsVt5cuRj7EtxrLjpR0Uz1OcFtNb0Hpmaw5fOGx1aEIIIUSyJDETQggPEKtjmRU2i46BHfHxkpFOUlK1SFVWPbeKWR1msePkDgInBPL2b29z9eZVq0MTQgghEiWJmRBCeIBNxzfxT+Q/0owxDZRSdKzckf199zO44WBGbRpFhfEVmLl3Jlprq8MTIssaPnw4QUFB5MmTh8KFC9OmTRsOHDhgdVhCWE4SMyGE8AC2UBv35LmHhvc2tDoUj5PTNyfvP/w++/rs44F7HuCpuU/ReGpjdp3aZXVoQmRJa9eupU+fPmzevJmVK1cSHR1Ns2bNuHpVarRF1iaJmRBCuLmY2Bhmh82mY2BHvJScttOrTL4yzOs8jxXPrODctXPU+qYWfX7uw/lr560OTYgs5ZdffqF79+5UrlyZ6tWr891333Hs2DG2b99udWhCWEoeVBBCCDe37u91nLpyis5VpBmjIzx636PsenkX47eOZ9jaYcwMnclHD39Er9q98Pbytjo8ITzKtVvX2H9uf7LzVCxYkZy+OZN8PyIiAoD8+fM7NDYhPI0kZkII4eZsoTZKBZSi7j11rQ4l0/D19qV//f48XfVp3vrtLXov7c3X279mXItxPFTqIavDE8Jj7D+3n9rf1E52nu29tlOrWK1E39NaM2DAAB588EGqVKnijBCF8BiSmAkhhBuLjo1m7r65vFDjBZRSVoeT6RTJXYTJrSfzUp2XCF4WTKOpjXiqylOMeHQE9/jfY3V4Qri9igUrsr1X8k0QKxasmOR7ffv2Zffu3WzYsMHRoQnhcSQxE0IIN7bqyCrOXTsnzRid7IF7HmBTj01M2zWNN399kwrjK/D2Q28zoP4Asvlkszo8IdxWTt+cSdaGpSQ4OJhFixaxbt06SpQo4eDIhPA88hS5EEK4MdteG/fnv5+aRWtaHUqm56W86F6jO+F9w3mp9ku8u+ZdKk+ozJLwJdK9vhAOpLWmb9++zJs3j1WrVlGmTBmrQxLCLUhiJoQQbupmzE3m7Z9H58qdpRmjCwVkD2BU81Hsfnk3ZfOVpdVPrWg5oyXh58OtDk2ITKFPnz78+OOPzJgxgzx58nDq1ClOnTrF9evXrQ5NCEtJYiYytZCQEAIDAwkKCrI6FCHSbMXhFVy6cYkuVbpYHUqWVKlQJZY/s5z5neez79w+qkyowpsr3+Ry1GWrQxPCo3311VdERETQpEkTihUrdvtls9msDk0IS0liJjK1Pn36EBYWxh9//GF1KEKkmS3URmChQKoUlp7KrKKUok3FNoT1DuPdxu8ybus4yo8vzw+7fiBWx1odnhAeSWud6Kt79+5WhyaEpSQxE0IIN3Qj+gYL9y+kc2Xp9MMd5PDNwdBGQ9nfdz+NSjXiuQXP8eCUB9l+QgbEFUII4RiSmAkhhBtadnAZl29elsTMzdwbcC+2DjZWd1vNlZtXCPo2iF6Le3H26lmrQxNCCOHhJDETQgg3NDN0JjWK1qBCwQpWhyIS0aR0E/586U/GthjL7LDZlBtXjrFbxhIdG211aEIIITyUJGZCCOFmrt68ypLwJVJb5uZ8vHzo+0BfDgYfpEuVLvT7pR81JtZg1ZFVVocmhBDCA0liJoQQbmZJ+BKu3bpGp8qdrA5FpELBnAWZ+MREtvXaRkD2AP437X90nN2Rvy/9bXVoQgghPIgkZkII4WZsoTaCigdRNl9Zq0MRaVCrWC02PL+BH9v+yMZjG6kUUokP1n7A9VsyNpMQQoiUSWImhBBuJDIqkqUHl0ozRg+llKJrta4c6HuAV+u+ykfrPqJSSCXm7ZuH1trq8IQQQrgxScyEJZRSvZVSR5RSN5RS25VSD6Xycw2VUtFKqZ1ODlEISyzcv5ComChpxujh8mTLw6dNPyW0dyiVC1em/az2NPuxGWFnw6wOTQghhJuSxEy4nFKqMzAG+BioCawHliml7k3hcwHANOA3Z8cohFVsoTYalmxIyYCSVociHKBcgXL8/PTPLHlqCUcvHaXaV9Xo/0t/Lt24ZHVoQggh3IyP1QGILGkAMFlrPcn+dz+lVHPgFWBIMp/7GpgBxABtkltBVFQUUVFRt/+OjIzMSLxCuMTF6xdZcXgFo5qNsjoU4WAty7ekadmmjNk8hg/Xfcj0PdMZ/r/hPF/zebyUlJEKz3L9unXPTVq5biGcTRIz4VJKKT+gNvBpgrdWAA2S+dzzwH3AM8DQlNYzfPhw3n///QxEKoTrzd8/n+jYaDoEdrA6FOEE2Xyy8eaDb/JMtWd489c3eXHxi0zcPpFxLcZRr0Q9q8MTIkU+Pj54eXlx9OhRS+PQWhMTE2NpDEI4g5KHkYUrKaWKA/8CDbXWv8eb/hbQTWt912i6SqlywAbgIa11uFJqGNBGa13DPstdO3FiNWYlS5YkIiICf39/R34lIRym2Q/NuBV7i9XdVlsdinCBDcc2ELwsmJ2ndtKtejc+bfopRXMXtTosIZJ18+ZNoqOjiYiIYMmSJfj6+uLn5+fSGKKjo7l16xZRUVHcvHmT5557jvz588e9rVwajBAOJDVmwioJkymVyDSUUt6Y5ovvaa3DU7vwbNmykS1btoxFKIQLnb16llVHVhHyeIjVoQgXefDeB9nWcxuT/pzE26veZt6+ebzX+D2C6wbj5+3aG10hUsvPzw8/Pz+01nh7e3Pp0iVL4/H395frvcg0pMZMuJS9KeM1oKPWen686V8CNbTWjRPMnxe4iHmuLI4XJpGLAZpprVPsDCQyMpKAgACpMRNua+K2ifRd2pdTr5+iYM6CVocjXOzC9Qu8u/pdvtr2FeXyl+PLx76k+f3NrQ5LiGRdvXr1jtYpVsiWLRu5cuWKP0lqzITHksRMuJxSaguwXWvdO960MGCh1npIgnm9gMAEi+gNPAJ0AI5ora+ktE5JzIS7e/j7h/Hz9mP5M8utDkVYaPfp3by67FXW/r2W1hVaM7r5aBloXIi0kcRMeCzpCkpYYTTwolLqBaVUJaXUF8C9wEQApdRwpdQ0AK11rNZ6b/wXcAa4Yf/7qmXfQggHOXn5JGuPrpVBpQXVilRjdbfV2DrY2H5yO4EhgQxdNZSrN+VUJ4QQmZ0kZsLltNY2oB/wLrATaAQ8rrX+2z5LMUyiJiyiNURHWx1F1jE7bDY+Xj60rdjW6lCEG1BK0alyJ/b32c8bDd5g5O8jqRhSEduu6ejYWKvDE0II4STSlFFkBinuxNKUMfVu3oQnn4Rbt+A3GcrbJRpOaUj+HPlZ/NRiq0MRbuivi38xcFk/FhxcTKOLAYxrP4lqD8mQCkIkQZoyCo8lNWZCiNu0ht69YflyWLUKtm61OqLM73jEcX4//rs0YxRJKpuvLPOzd2f5D3BGXaPmrx3pO7gaF04ctjo0IYQQDiSJmRDithEjYPJkmDIFypSBceOsjijzmxU6i2ze2XiywpNWhyLcmc1Gs9zV2fXJBUbkeJJpXnsoP7YcX3/RlZhbN62OTgghhANIYiaEAGDePHjzTXj7bXj+eejTB2w2OH3a6sgyN1uojcfLPY5/NmliK5Jw9SosWQKdO+OXIzcDBi8k/MVdtIq+j5cjZ1BnUF42LJlgdZRCCCEySBIzIQR//AHPPAOdOsEHH5hpL7wAvr7wzTfWxpaZHb5wmD9O/CHNGEXyFi+Ga9eg83/7SdGy1fhu5EE215uEL148tL0Pzwwsw7/h2ywMVAghREZIYiYytZCQEAIDAwkKCrI6FLd17Jjp7KNaNZg6FbzsZ4V8+UyyNnGi6QhEON6s0Fnk9M3JE+WfsDoU4c5sNggKgrJ3j2dWt3kPNo+8xOR83Vnh8zcVvg/i048eI+pqpAWBCiGEyAjplVFkBtIrYzpFRsKDD5p/t2yBIkXufH/PHpOwzZx5R2G9cJAaE2tQsWBFZnaYaXUowl1FRkLhwvDxxzBwYLKzXjr9Nx982Y6xvn9S5oovY2q9Rcuuw1wTpxDuQ3plFB5LasyEyKKio6FLF/j7b/j557uTMoCqVaFJE+kExBkOnDvArtO7pBmjSN7ChRAVZdoZpyBvkVKM/mQ7u1sspFRMbp449D4t+xcmfNtyFwQqhBAioyQxEyKLGjAAVqyAWbOgcuWk5wsOho0bYccO18WWFdhCbeTxy0OLci2sDkW4M5sNGjSAkiVT/ZHAek+ycvQ55pUaRJj3RaoseozBQ+tx+fwJJwYqhBAioyQxEyILGjfuv1fz5snP++ST5p5Qas0cR2vNzL0zaV2xNdl9slsdjnBXFy6YQQXT0Y5YeXnRtvtnhL17mqHeD/MlW6jwWUl+/OoVdGysE4IVQgiRUZKYCZHFLF0K/fpB//7wyispz+/jYwadnjEDzp1zenhZwt4ze9l3bh9dKnexOhThzubPh5gY6Ngx3YvI4Z+fd99bxf6nNtAwuhjPnpnIgwPy8ueq6Q4MVAghhCNIYiZEFrJ7tyl8b9nSDCadWi++aP6dNMk5cWU1tlAb+bLn49H7HrU6FOHObDZo3BiKFcvwokpVbsjs0f/wW9WRRHjdpM66Z3hpUCXOHT/ggECFEEI4giRmQmQRJ0/CE0/A/feb2i9v79R/tmBBePppmDDBdBoi0k9rjS3URtuKbfHz9rM6HOGuzp6FVasc3h3qI+0GsvPTS3yZqz2zvA9QbkIlxo/oSPTNGw5djxBCiLSTxEyILODaNfOsWEyMGas2d+60LyM4GI4fh0WLHB9fVrLj1A4OXThE5yrSG6NIxty55t/27R2+aB+/7AS/MYfwV0LpGFOBV6/OodbgfKxZMMbh6xJCCJF6kpgJkcnFxsKzz0JYmEnKSpRI33Jq1oSGDaUTkIyauXcmBXMW5JEyj1gdinBnM2fC//4HhQo5bRWF7q3EN5/v44+HfiB3rA8P7+pP5wH3cixsk9PWKYQQImmSmAmRyb31lulDYMYMqFUrY8sKDoY1a8zA0yLttNbMCp1Fh0od8PHysToc4a5OnIB161w2qnvt/z3DhlEXmVboJdb5/EvFGQ348IP/cePKJZesXwghhCGJmRCZ2OTJ8NlnMHIktG6d8eW1awfFi8P48RlfVla05d8t/B3xtzRjFMmbM8d0h9q2rctW6eXtw7O9JxL+5nH66iA+jF5F4LBCLPh+iHSvL4QQLiKJmRCZ1KpV8PLL8NJLpmt8R/D1Ncv88Ue4eNExy8xKbHttFM1dlIfufcjqUIQ7s9mgWTPIl8/lq85ToDiff7yVPa2WUiEmL22PfspjAwqzb8sSl8cihBBZjSRmIlMLCQkhMDCQoKAgq0Nxqf37TZ8BDz9snglTynHL7tXL9Mw4ZYrjlpkVxOpYZoXNomNgR7y90tAlpshajh2D3393WTPGpFQIasHSUadZVHYoh3wiqPZzKwa+VYeIM8csjUsIITIzpbW2OgYhMirFnTgyMpKAgAAiIiLw9/d3RUyWOXcO6taF7Nlh40bIm9fx63j2WbPsgwfT1u1+Vrbu73U0ntqYjS9spEHJBlaHI9zVyJEwdCicOQNucq66ceUSX4zpzEfXV5An2otPSz7Pc69MxMtbnpMUbsmBRZFCuJbUmAmRiURFQZs2cPkyLFninKQMTCcgR47A0qXOWX5mZNtro6R/SeqVqGd1KMKd2Wzw+ONuk5QBZM+dlyFDl3Pgua08El2S589Ppv7redm64jurQxNCiExFEjMhMgmtoUcP2LYNFi6EMmWct64HHjAv6To/daJjo5mzbw6dKnfCS8lpVyTh8GFzAFvcjDEpJSoEMWPUUdbVGkcUMdTd9AI93ijP6SN7rQ5NCCEyBblDECKT+PBDmD4dvv8e6td3/vqCg2HlSvM8m6ucuHyCt397m38j/3XdSh1g7dG1nLl6hs6V3fOGW7iJWbMgZ0544gmrI0nWQ636sv3zCCbk7sIC70OU/7YqX3zahls3rlkdmhBCeDRJzITIBGbMgPfeg48+cl1he8eOULiwa7vO/3Lzl3yy4RMqjK/A8PXDiYqOct3KM2Dm3pmUzVeWOsXrWB2KcGc2m0nKcuWyOpIUefv68crAnwjve4CuMZV5/fpCqr+Vj1/nfG51aEII4bEkMRPCw23cCM8/D889ZwaTdpVs2UxX/N9/D5GRzl+f1ppZYbPoWrUrvWr34t0171J5QmWWhC/BnTsxuhVzi3n759G5cmeUI7vHFJnL/v2waxd06WJ1JGlSoEQ5Jny2l+2PzKRgbHYeDX2Tdv2Lc2T3OqtDE0IIjyOJmRAe7K+/TGcf9erBN984tlv81Hj5ZbhxA6ZOdf66tv67laOXjtKjZg9GNx/Nrpd3USZfGVr91IqWM1oSfj7c+UGkw69//cqF6xekGaNIns0GefJAixZWR5IuNRp3Zu3oi/xULJitPqepNLsx777XiGsR56wOTQghPIYkZkJ4qEuXoGVL0/PivHmmBsvVihc346WNHw+xsc5dly3UDM7cqFQjAAILBbLimRXM6zSPfef2UWVCFd5c+SaXoy47N5A0soXaqFCgAtWKVLM6FOGutDaJWevWZpwLD6W8vOjSayz7h/zLQNWAz2LXU/GjosyePADt7BOEEEJkApKYCeGBbt2CDh3g9Gn4+WcoUMC6WIKDzXhmK1Y4bx2xOpZZobPoUKnDHYMzK6VoW6ktYb3DeKfRO4zbOo7y48vzw64fiNXW3whGRUcxf/98acYokrd3L+zb57a9MaZV7vxF+fiDjYS1/40a0QXp9M8XPDKgAHs2zLU6NCGEcGuSmAnhYbSG3r1h3TqYPx/Kl7c2ngYNoGZN53adv/HYRv69/C9dqiT+/E0O3xy80/gd9vfdT6NSjXhuwXM8OOVBtp/Y7rygUuGXQ78QGRVJ5yqZ44ZbOMnMmabqu1kzqyNxqPtqPMKiL06xtPz7nPC+Ss2VHXh1cHUunjxidWhCCOGWJDETwsOMGgWTJplnyho3tjoa81xbcDAsWwaHDjlnHbZQGyX8S1C/ZPLjANwbcC+2DjZWPbeKyzcvE/RtEL0W9+Ls1bPOCSwFtlAbVQtXJbBQoCXrFx4grhlju3bg52d1NE7R4ql32fPRBYZne5zvvHdT7sv7+GbMM8Tcuml1aEII4VYkMROZWkhICIGBgQQFBVkdikPMnw+DBpneF7t3tzqa/3TpAvnzQ0iI45cdExvDnLA5dApM/eDMD5d5mB0v7WBsi7HMDptN+fHlGbdlHNGx0Y4PMAnXbl1j0YFF0umHSN6ff5qBpTNJM8ak+OXIzRtv/Uz4CztoGV2GlyKm88CgvPy+9GurQxNCCLchiZnI1Pr06UNYWBh//PGH1aFk2Pbt0LWrebbsww+tjuZOOXJAz54wZQpcueLYZa/9ey2nr55Oc3NAHy8f+j7Ql/C+4XQK7MRrv7xGza9rsvrIascGmISlB5dy9dZVacYokmezQcGC8MgjVkfiEsXuq8H3Iw/z+wNfo1A0/ONlnh1YlhMH/7Q6NCGEsJwkZkJ4gOPHoVUrqFbNjBvm5YZH7iuvmKTshx8cu9yZe2dSJm8Zgoqnr9azUK5CfN3qa7b12oZ/Nn8emfYInWZ34ljEMccGmoAt1EatYrW4P//9Tl2P8GBxzRjbtwcfH6ujcan6LXqxdWQEk/J1Y7nPUSpMrc3nHz9O1FUXDIoohBBuyg1v74QQ8V2+bJIyPz9YuNDUTrmje+81Y6qNH2/uNx3hVswt5u6b65BeDWsVq8WG5zfwQ9sf2HBsAxXHV+SDtR9w/dZ1xwQbz+WoyywJXyLNGEXyNm+GY8cyfTPGpHh5+9Dj1amEDzhCj9gavHVzGVXfLciynz6wOjQhhLCEJGZCuLGYGHjqKTOQ9JIlUKSI1RElLzgYwsJg1SrHLO+3I7+ZwZkd1BxQKcUz1Z7hQN8DvFr3VT5a9xGBEwKZv28+2lHZJLA4fDE3om/QqXInhy1TZEI2GxQtCo0aWR2JpfIWKcWY4TvY1XwBJWNy83j4e7TqX4RDO36zOjQhhHApScyEcGMDB8Ivv8Ds2VClitXRpKxxYxOno7rOt4XaKF+gPNWLVHfMAu3yZMvDp00/ZW/vvVQqWIl2s9rR7MdmhJ0Nc8jybaE26pWoR+m8pR2yPJEJxcaaA7tjR/D2Tnn+LKByg9b8Ovocc0q+zm6f81Se15Qh79TnyoVTVocmhBAuIYmZEG4qJAS+/NIkOc2bWx1N6sR1nb94MRw9mrFlRUVHMX+fcwdnLl+gPEu7LmXJU0s4euko1b6qRv9f+hNxIyLdy7x04xK/HPpFmjGK5G3YACdOZNlmjElRXl60f2EE+4aeYoh3Y8bozVT49B5mTOyDjrV+0HghhHAmScyEcEO//AKvvgr9+plONTxJ167g7w8TJmRsOcsPLyciKiLJQaUdqWX5lux9ZS8fPfIR3/75LeXHl2fKjinE6rTfCC7cv5BbMbfoGNjRCZGKTMNmgxIloH7yY/NlVTkDCjJs2Br2dVlP/ehidD09gUYD8rFzrc3q0IQQwmkkMROWUEr1VkodUUrdUEptV0o9lMy87ZRSK5VSZ5VSkUqpTUopD6lDSrs9e6BTJ3j8cRg50upo0i5XLujRwwyCfe1a+pdjC7VRpXAVlw3OnM0nG4MfHMyBvgd4tOyj9FjUg3qT6rHlny1pWs7M0Jk8eO+D3ON/j5MiFR4vOto0Y+zUyT27WHUjpas8yJzR/7Cy8mdc8Iqi9qouvPJmZc7/c9Dq0IQQwuHkiiBcTinVGRgDfAzUBNYDy5RS9ybxkUbASuBxoDawGlislKrp/Ghd69QpeOIJuO8++Oknz330pHdvuHQJZsxI3+ev37pu2eDM9/jfw4/tfmT98+u5FXuLepPr8fzC5zl1JeXnXM5fO8+vf/3qklo+4cHWrIGzZ83I7CJVmnYYxM5PLjA6V1t+8gqjXEgFJozqQvTNG1aHJoQQDiOJmbDCAGCy1nqS1nqf1rofcBxItNGe1rqf1vpzrfUfWuuDWuu3gINAK9eF7HzXrkHr1nDrlnlGK3duqyNKv7JlTYI5blz6us5fenApV25esfQ5rQfvfZBtPbcxseVEFh9YTPlx5Rn1+yhuxtxM8jPz9s0jVsfSvlJ7F0YqPI7NZg6SOnWsjsSj+GbPyWuD5hH+0l7aR5en72Ubtd/Mx9qFX1odmhBCOIQkZsKllFJ+mFqvFQneWgE0SOUyvIA8wIWk5omKiiIyMvKOlzuLjYVu3WDvXpOUlShhdUQZFxwMu3fD+vVp/+zM0JnULFqTcgXKOT6wNPD28ualOi8RHhzOc9WfY9Cvg6g+sTorDifcfQ1bqI2HSz9MkdxuPq6BsM6tWzBvnmnG6KRObTK7wqUr8+2I/WxpOJWc+NBkZz+eGliK4/vS1uxYCCHcjSRmwtUKAt7A6QTTTwNFU7mMgUAuYFZSMwwfPpyAgIDbr5IlS6YrWFcZOhTmzoXp06F2baujcYymTaFixbR3nX/l5hV+Dv/ZrZoD5s+Rn/GPj2fHSzsokqsIzX9sTpuZbfjr4l+35zl95TSrj66W3hhF8n79FS5ckN4YHSDo0W5sHHmR7wv2ZLXPP1T8sR4ff/goN65csjo0IYRIF0nMhFUSNnBTiUy7i1LqKWAY0FlrfSap+YYMGUJERMTt1/HjxzMUrDN99x0MHw4jRkCbNlZH4zhKQd++MH8+pOXnX3xgMdejr7vl4MzVilRjdbfVzGw/k+0ntxMYEsg7q97h6s2rzAmbg5fyol2ldlaHKdzZzJlQvjxUd+zYfFmVl7cPz/X5hvDX/6Y3dRh261cqDyvEoh/elu71hRAeR+n0PAAiRDrZmzJeAzpqrefHm/4lUENr3TiZz3YGvrN/9ud4b6W4E0dGRhIQEEBERAT+/v7p/wIOtno1NGsGL7wAEydmvpZNly/DPfeYZo0ff5y6z7SZ2YZTV06x+cXNzg0ug67evMqnGz5lxO8jKJSrEDm8snF/QTMumhCJunEDihQx42C8/77V0WRK+7cupd9Pz7E873maXyrAmKemUfGBx60OS2TQ+fNw5QqUKpWq2TPZlVRkJVJjJlxKa30T2A48muCtR4Hfk/qcvaZsKvB0gqTMYx04AO3bQ5MmMH585kvKAPLkgeefh2++MfekKYm4EcGyQ8s8ojlgLr9cfPjIh4T1CaN2vsocvHSYrouOmvEOhEjM8uUQGSnNGJ2o4gOPs2zUGRaWfYuD3hFUXdKSN94OIvLsP1aHJtLp5k1o1w5atjTPYwuRmUliJqwwGnhRKfWCUqqSUuoL4F5gIoBSarhSalrczPakbBrm2bLNSqmi9leAFcE7wrlz5iJTtKgZzsjX1+qInKdPH/N9bakYF3bB/gXcjLlJx8qeMzhz2XxlWXCxOYcm+PL0rhioUcNUEV5Ism8akVXZbFClCgS6Zmy+rEp5efHksx8TOuws7/s+ygS2UX5kKb4P6UlsTLTV4Yk00Bp69YLNm+Hrr2XYP5H5yS4uXE5rbQP6Ae8COzHjlD2utf7bPksxTKIW5yXABwgBTsZ7eWQfyVFRpvQvMhJ+/hny5rU6IucqXx4eeyx1XefbQm08dO9DlPD3sG4pbTbuq9cCtXsPfP45fP+9+eLffAMxMVZHJ9zBtWuwaJHUlrlQ9tx5eeudFex/ZjMPR5eg+7lJNHw9H9t+nZbyh4VbGD7cnE6/+w4aNrQ6GiGcTxIzYQmt9QStdWmtdTatdW2t9bp473XXWjeJ93cTrbVK5NXditgzQmvo2RO2boUFC6BMGasjco3gYNi+3ZR6JuX8tfOs/GulRzRjvMPRo7Bli7nh9vODgQNNO9UnnoCXXoKgINi40eoohdV+/hmuXpXEzAIlK9Xlp1F/s7bmWK4RzQMbuvHiGxU4czTU6tBEMmbNgrffhmHD4OmnrY5GCNeQxEwIF/r4Y/jhB1P61yBVo7ZlDo89Bvffn3zX+fP3zydWx9IhsIPrAnOEWbMge3ZoFW+882LFYOpU2LQJvL3hwQfh2WfhxAnLwhQWs9mgVi0oZ+3YfFlZoyeD2f7ZRcbn6cw8n4OU/7oKX37ejls3rlkdmkhg82Z47jmTkL37rtXRCOE6kpgJ4SIzZ8I778AHH8BTT1kdjWt5eZlnzWbPhpMnE5/HFmqjSekmnjc488yZ5oHBPHnufq9ePVObNmmS6fihfHn47DPTnlVkHZcvmxozqS2znI9fdnoPnMnBPgd4KjaQ/tfmU/Ot/Pw2d4TVoQm7o0ehdWuoUwcmT86cHWMJkRRJzIRwgU2boHt3eOYZM5h0VvT885Atm3mAO6HTV06z6sgqz2vGePAg7NgBXZIZDNvLC3r0gPBwePFF0zanalVYKt3qZxmLF5tuSTu539h8WVWBEuX46rNQtjeZQb7YbDTdO4gOA0pwdO8Gq0PL0iIiTCvw3LnNGJjZs1sdkRCuJYmZEE525Igp/QsKMhUnWbX0LyDANE2ZONF0fxzf3H1zPXNwZpsNcuWCx1MxTlLevDBmDOzaBSVLmlq2J56AQ4ecHaWwms0GdetC6dJWRyISqNnkKdaNvsiMon3Z5HOSSjMf4v33H+Z6pPSq6mrR0abs4p9/TAVzoUJWRySE60liJoQTXbpk7r/9/U3pX7ZsVkdkrb594fRp06QxPluojaZlm1IwZ0FrAksvmw2efBJy5kz9ZypXhl9/hblzYe9e8/eQIWb0VJH5XLoEy5ZJM0Y3pry8eOqlcRwY/C/9VX0+iV5DpQ8KM3fKG2gZOMsltDadRK1aZU6NFStaHZEQ1pDETGRqISEhBAYGEhQU5PJ137plSv9OnjSlfwU9LOdwhsBA+N//7uwE5MTlE6z/e73nNWMMDTWJVXpuuJUyYyaEhcFbb5matAoVYMaMlMcUEJ5lwYL/qgKEW8udvyiffPg7oe1+pWpMATocH8mjAwoS+vtCq0PL9L780rSm+Oorc40QIquSxExkan369CEsLIw//vjDpeuNK/1bvRrmzTP33MIIDjb9YcRtktmhs/H19qVNxTaWxpVmNptpn/nYY+lfRs6c8N57sH+/6aaza1d46CHz3JrIHGw20yvnPfdYHYlIpftr/o/FX5zm53LDOOZ9herL29BvSE0unf475Q+LNFu8GAYMgEGDzGO4QmRlkpgJ4QRffGE6ufjmG3j4YaujcS9PPGEetYmrNZsZOpPm9zUnb/a8VoaVNlqbG+42bRzTPrVUKdO+87ffTNO32rXh5Zfh3LmML1tY5/x502xVmjF6pMeffo+9H13gE78WTPbaSfnRZZg8tjuxMdFWh5Zp7Nhheilu29YMJi1EVieJmRAOtnAhvP46DB5seiIUd/L2ht69TV6z/fDfbP5nM12qJNOroTvatcv0sujoG+5HHjF3KmPGmG74y5eHkBDTFE54nnnzIDYWOnjY2HziNr8cuRn09lIOdN/OYzFlePHi99R9PYDNv0yyOjSP9++/ZvjHSpXM+J5eckcqhCRmQjjSn3+aATHbtTODSYvE9ehhErS3p88iu092WpVvlfKH3InNBvnzQ9Omjl+2ry+8+qrpir99e9P2s3ZtWLvW8esSzjVzJjRpAkU8bGw+cZfi5WoxbeRhNtSZQCya+lt60v31+zn1126rQ/NIV66YpMzLCxYtSlv/SUJkZpKYCeEg//xjLjSVK8O0aVL6l5z8+c3jVKvO2GhxX0vyZEtkcGZ3FdeMsV07k0Q5S6FC8O23sHWruWtp0sSMl3b8uPPWKRzn9GlYsyb5Me6Ex2nY8hW2fn6JbwKe5Wefvyg/qTojhz/BzevSq2pqxcSY8//Bg7BkCRQrZnVEQrgPuXUUwgHiSv98fKT0L7VaP3+IW4W2U+qKhz1/88cfZnA6V91w16kDGzfC99+bWrOKFU117I0brlm/SJ85c0zpTDsPG5tPpMjb14+e/aYR/tphusdWY/CNn6k2tADLbdJMIjXefNMkZDYbVKtmdTRCuBdJzITIoJgY03zx8GFzsSla1OqIPMOuaBte0bnY8kNLq0NJG5sNCheGxo1dt04vLzM694ED0KcPvP++GXtg4ULpXt9d2WymqWuBAlZHIpwkX7EyjP10FzsenUOxmJw8tn8orfsX469da6wOzW19/TWMGmW6x3/8caujEcL9SGImRAa98YYZp2zWLKha1epoPIct1EaDAk+yaV1Odu60OppUio01G7pDB1M96mr+/vD557Bnj+kYpE0baNHCdLcv3Mc//8CGDdIbYxZR9cH2rBp9nlkl+rPD5yyBcx5m6LsPcvXiGatDcysrV5pypeBg6NvX6miEcE+SmAmRAV99ZbrGHzcuY8NZZTX7zu5jz5k99H+0MyVL3jngtFvbtMncdFt9w12hAixbZtrNHjxoSgRefx0iI62NSxizZ5vnD9u0sToS4SLKy4uOPUazf+gp3vRuxEi9kYqfFGfmN6+iY2OtDs9yYWGmPKt5cxg92upohHBfkpgJkU7Ll5uSv9deM92/i9Szhdrwz+ZPywqP8corMGOGGfLJ7dlsULy4GTDYakqZBxtDQ03Txq++MrVo339vavaEdWw2U1KTN6/VkQgXyxlQkPeHrWVfx7UERRfmqZPjaDIgP7vWzbI6NMucOQMtW5rhGmfOtKaxgRCeQhIzIdJh717o2NHce40aZXU0nkVrjS3URpuKbcjmk42ePc1jUpPcfVigmBhTE9Kpk3t1uZk9O7z1lnn+7JFHoHt3aNDAdFIiXO/oUdiyxfpaVWGpMtUaMe+LE6wIHM5Zr+vU+q0zfd6syoUTh60OzaWuX4fWrU1fRUuWQB4P6oBXCCu40d2FEI4XEhJCYGAgQUFBDlvmqVOm9K9sWfjpJzMel0i93ad3s//cfrpUNr0aFiwITz0FEya4+TjK69aZje+uN9wlSpiqx7Vrzd1Q3brw4oumuFq4zqxZJllu5WFj8wmneLTjYHZ9cpGROVrzo/deyo0tx8TRTxNz66bVoTldbCw8/zzs2mVaXd97r9URCeH+JDETmVqfPn0ICwvjDwfVHsSV/t26BYsXS+lfethCbeTPkZ+mZf8bnDk4GI4dM7+p27LZTFucunWtjiR5jRrB9u0QEgLz55vmjWPGmJ1WOJ/NZkpu5OQg7Hyz56T/4AWE99xD65j7eeXyT9QZlJf1i8dbHZpTDRtmDocffwQHlo0KkalJYiZEKsXGQrduphnj4sVQsqTVEXmeuGaM7Sq2w9f7v8GZa9Uyre/cthOQW7dg7lzTjFEpq6NJmY8PvPIKhIebsRwGDoQaNeC336yOLHM7eBD+/FMGlRaJKlKmClNGhLOl/hT88KLRn8F0HViaf8O3WR2aw02bBh9+CJ99JkP5CZEWkpgJkUrvvGPGjJ0+HWrXtjoaz7T95Hb+uvgXnavc3RwwOBhWrzaJr9tZtQrOnfO8G+4CBUwb0e3bzf+bNoX27c1zUMLxbDbIlUsGaBLJeqDZ82waeYkp+Z/nV59jVPg+iOEfNSfqauboVXXdOtOKukcPM5yMECL1JDETIhWmToVPPjFDSEkP2Ok3c+9MCuUsRJPSTe56r317KFYMxrtj6x6bDe6/H2rWtDqS9KlRwzx7NmOG6ZiiUiV47z24ds3qyDIXmw2efBJy5rQ6EuHmvLx9eD54CuEDj/KSrs27N1dQ+d0CLPnxXY/uXv/QIWjbFh56yJQJeUIDAyHciSRmQqRgzRro1Qt69jQtwkT6xOpYZoXOomNgR3y87u4v2dcXXn4ZfvgBLl60IMCk3LxpntXq3Nmz7zKUMr2s7N8P/fvDp5+aBG3OHNMtpsiYsDBT3euuncMItxRQ+F5GfbKN3S0XUzbGn1aHP6TlwKKEb1tudWhpduGCebyyUCFzWvHzszoiITyPJGZCJCM83LSPb9TI9KXgyfflVtv8z2aORx5PtBljnF69zONc333nwsBSsmIFXLqUeW64c+c21b+hoVC9uhn3oWlTN21D6kFsNvD3l5HmRbpUqvsEy0efZX7pwezzvkiVRY/x5tC6XD5/wurQUuXmTdPq4fx5+PlnyJfP6oiE8EySmAmRhPPnTelfkSKm9M/XN+XPiKTZ9toonqc4D96b9ODMRYua/jVCQsywYW7BZjM1S1WqWB2JY91/v+nDeulS+Ocf09zxtdfcrLrSQ2htRs5t2xayZbM6GuGhlJcXbboNJ+zd07zj8wjj2Er5z0ryw4SXiY1x37FEtDatHX7/HRYsgPvuszoiITyXJGZCJCIqytxjXbpkSv/y5rU6Is8WExvD7LDZdAzsiJdK/rQTHAx//QXLlrkouORcv27uNLp0ybzVpS1awJ49MHw4TJliutefNMmNMmMPsGuXqV7PLLWqwlI5/PPzzru/sf/p32kUfQ/Pnf2aBwfmY/tvP1odWqI++8y0cpgyBR5MutxNCJEKkpgJkYDWpknd1q2wcKEZSFpkzPpj6zl55SSdK6d841q3rhnzxi26zl+2DK5cyfw33H5+pvu08HCTqPXsaTbEpk1WR+YZbDbIn980CRXCQe4NrI9t9DFWVRvFZa9bBK1/ll6DKnH22D6rQ7ttzhwYMsT0JdS1q9XRCOH5JDETIoFPPjFjsHz3nRlbS2Scba+NUgGlqFeiXqrmDw42j3YdOODkwFJis5nnsCpUsDgQFylWzOz8GzeaEooGDeC55+DkSasjc19am/2kXTtp7yyc4uG2A9jx6SXG5urAbO8DlP+qMuNGdCD65g1L49q6FZ591vQp9N57loYiRKYhiZkQ8dhsMHQovP++udiIjIuOjWbuvrl0qtwJlcrmgJ06QeHCFnedf/UqLFmS+WvLEtOggbnr+uYbU2tYvjyMGGGe8Bd32rYNjhzJmvuJcBkfv+z0fWM24a+E0immAq9dnUvNwXlZNW+UJfH8/bcZGaJWLdOEMbO29BbC1SQxE5laSEgIgYGBBAUFpTjvpk3QrRs884wZTFo4xuojqzl77WyqmjHGyZbNNCedOhUirRpzdckSM85XVr3h9vY2TRrDw+GFF0x7papV3eThPzcyc6YpRWjSxOpIRBZQ6N5KfP35PrY1+hH/WD/+t+d1Og4owd+hG10WQ0QEPPGEGa5vwQLInt1lq3YqpZSPUqqnUmq1Uuq8UuqGUuovpdRcpVRrq+MTWYPSMn6N8Hwp7sSRkZEEBAQQERGBv7//Xe8fOWIeqalQAX79VTpWc6QXF73ImqNrOBh8MNU1ZgD//gulS8Po0aZpo8u1bWuC2LrVgpW7ob174dVXYfVqc1f2xRemZ8esLDYWSpUyVQchIVZHI7IYHRvLjG/68sZfX3PJL5bBvg/zRv855PDP77R1Rkebw3/zZlOYWamS01aVEWmuv1NK5QMWAw3tk/4GzgL3AoWBv7XWpR0VoBBJkRozkeXFlf7lyWPGEZakzHFuxtxk7r65dK7cOU1JGcA995jHdsaPN/e/LhUZaWqGsmptWWKqVIHffoPZs2H3bqhcGd56y3SOklVt2mSGGpD9RFhAeXnR9eUJHHjzOK9Sl49iVlPpg8LM+24Q2gknTa1N2cxvv8HcuW6blKXXFExSdhiop7UurbUO0loXAcoBUvIiXEISM5Gl3bplxtc9ccJ0i1+woNURZS4rD6/k0o1LdKnSJV2fDw42LelWrnRwYClZuNCMmdCpk4tX7OaUgg4dYN8+07Txiy+gYkX46Sdz15bV2GxQvLj0ES4sladAcT79aDOhrVdQOSY/7Y+NoNmAQoRtXuTQ9YwdC199BRMmwP/+59BFW0opFQS0AaKAFlrrLfHf11of0lqPsCI2kfVIYiayLK3Njf/q1ab0r2JFqyPKfGyhNioVrESVwukbnLlhQzPuscu7zrfZTAcYJUu6eMUeImdOGDbMJGh168LTT0OjRrBzp9WRuU5MjKk97NgRvORSKqxXrvaj/PzFGZbc/x5HvS9TbVlr+g+pRcSZYxle9uLF0L+/GVWjZ08HBOte4p4fm6+1PmhpJCLLk6uJyLLGjIGvv4aJE+GRR6yOJvO5EX2DBfsXpKsZYxylTPK8dCkcPuzgAJNy4YLpq79L+mr5spTSpU2pxsqVcP481K4NvXub/2d269bBqVOynwi307LrMPZ+cI6P/R7jW68dlB9Vmiljnyc2Jjpdy9u50/RS3Lo1fPqpY2N1E3GNMjdbGoUQSGImsqhFi2DgQHjzTejRw+poMqdfDv3C5ZuX6VwlY8/fPPWUGbvXZX0rzJ9vnnDv0MFFK8wEmjaFXbtg1CiYPh3KlTPtnWJirI7MeWw20/FH3bpWRyLEXbLl8ufNt5dxoNsfPBpdih4Xp1Lv9bxsWT45Tcs5ccI8g12xIvz4Y6atHI7rEeySlUEIAZKYiSzozz/NzX7btmYwaeEcM/fOpFqRalQsmLE2ojlywIsvmrFyXNLPhM0GjRubwZZF6vn6Qr9+cPCg6bWlb19Tg7ZundWROV50tKkp7NRJBnASbu2e8nX4cdQR1tcO4Rax1Nv8Is+/Xo5Tf+1O8bNXr0KrVmYXX7QIcuVyQcDWuGz/N6+VQQgBkpiJLObECXOhCQyEH37ItKV/lrt68yqLwxfTpbJjmnm98gpcvmxKbJ3q7FlYtUp62cuIwoVh0iTYssUMcNS4sSkJOX7c6sgcZ9UqOHdO9hPhMR58ojfbPr/ERP+nWexzmPKTqjNqeCtuXk+8tCsmBrp2NZ0vLVli+rjJxELt/9azNAohkMRMZDGdO5txcxctMv0XCOf4+eDPXLt1LcPNGOOUKmWebxg/3smd/82da/5t396JK8kigoLg99/NKOGrV5u2UB9/DDduWB1Zxs2cCffdB7VqWR2JEKnm7evHS/2nE/7qQZ6LrcqgG0uoNrQAy20f3zXv4MGmw4+ZM6F6dQuCda0F9n/bKKXuszIQISQxE1lC3KMuhw+b0j9ppeZctlAbdYrXoWy+sg5bZnAwhIaae3ynsdlMP9CFCjlxJVmIlxd06wYHDphqz2HDzPhnixZ5bvf6N2+a5xC7dJFmjMIj5S9+H+M/3c2OprMpGpODx/YPpU3/Yvy1aw0A33wDI0ea0TBatrQ2VlfQWm8H5gPZgWX27vNvU0rdr5R63ZLgRJajtKdeHIUAlFIqIiLirpE0o6KiiIqKuv33i6/OZ/n8V/n+++O0aeOfcHbhQJFRkdw/9n6GNhrKq3VfddhytYb69aFsWZgxw2GL/c/Jk6ZWZ/x4ePZZJ6xAEB5uetxZtcokwJ9+CuXLWx1V2vzyi6l6//13k2QK4cF0bCzzpg1m6KFvOZc9ltaRzZkzdQovPpebESM8s+whICAgALis03CDq5TKB/wM1LdPOgqcA0oCRYC/tdalHRupEHeTxEx4NKWUPxBhdRxCCCGEcBsBWuvItHxAKeUL9ASeBqoA2YCTwHbge631EodHKUQCkpgJj5baGrMvRn7B6C9Ho/pBLyoz5JUZ5Cta2qGxBAUF8ccffzh0mZ64jk6zO3HxxkVWPrsy0fcjIyMpWbIkx48fx98/bbWXV69CpUrQvTt88EHy86b5ezz6KOTNawYNTiVP2B4pycj2SK1Ev8ONG2bk8JEjIU8es0G7dEl3jzyu2BYP1q7NhlOnTI+TQ4Y4fPmu2BaQOfZbV6zDsmPDgnWcPWsqsX0CwrmvwZOsyHuShpf8+bz1OKo0aOOQdWREWrZFemrMhHAXPlYHIERGpPbE2//1/oz+cjRDczzCF2oVc6bU4JOiT9OjzxS8ff0cEou3t7dTb6Y8YR0Xr19k1YlVjGw2MsVl+Pv7p3k9/v5m3LmpU2H4cNOVflLS9D2OH4etW01XnWmIyd23R1qkZ3ukVqLfwd8fPvwQevWCN94wz6B9/71J1urUccw6HKzJ9ev4X7linptz4rqcuS0g8+y3mfbYcPE6btwwrbdv3IAta+tQqtQJlts+5rWtH/DQpm68sqYqH7w6n/zFk+4Xw522RVpryoRwJ9L5h8hSevX+nvAXdtAyugwvRUzngUF5+X3p1w5Zdp8+fRyyHE9ex4L9C4iOjaZDoPMGZ+7TBy5eTPk5szR9j1mzIFs2ePLJNMbi3tvDXST7HUqWNF2/rVljqkQfeMAMXHfmjOPW4SDBRYqYLuoqZmxsPqtllv020x8bLliH1vD887Bjh+mTp1QpM71557fZ/dF5Ps/+BNO89lB+bDm+/qIrMbdupnkdQojUk6aMIjNIcSf+559/bjeDKFGiBACbln1D8Ir+bM97jWciy/DZy3MoXk66v86Ix358jBvRN1jTfU2S80RGRhIQEEBERES6S1hbtTKVXDt2OOjh9KAgkyDMm+eAhXkWR2wPh4mOhq+/hnfegdhYeP996N3bDF5ttatXzRhtQ4c6pRkjuNm2EFlie7z3nmlFPHs2dEiiPO3UX7sZMqE9U/McosalHIz730gefKK3S+NM47bwwC5LhDCkxkxkCdmyZbvjX4D6LXqx5fOLfJv3OX7xPUqFqbX5/OPHiboqrSDS4+zVs/z61690qZL8oNLZsmXjvffeu2NbpFVwMOzaBRs2pHsR/zl8GLZty7KDBTtieziMj4+pEg0PN8+b9e8PNWvCb79ZHZkZZ+PaNafuJ261LUSm3x4//miSsuHDk07KAIqWrcZ3Iw+yqe63+KB4aHsfnhlYhn/Dt7ks1sy+LYSIIzVmIjNIcSdOqbTt4skjDBvblhC/XZS94suXdYbS4ql3nRJsZvX1tq/ps7QPJweepFAu544DFhsLgYFQrZpphZghw4fDRx+ZpnO5cjkkPuEgO3aYLHzjRjPo98iRULq0NbG0a2eqaZ3cUYMQrrB+PTRtCs88A5Mmpb7lQWxMNFNDejL43++55qMZmqM5/fvPIlsut6pRlBoz4bGkxkwIIF+xMnw5fCc7m82jZExuHg9/j1b9i3BohxuU1HsIW6iNR8o84vSkDEzHfX37mpaH//yTwYXZbKZtpCRl7qdmTXMH+eOPZtywSpVM88br110bR2QkLF2aZWtVReZy6BC0bQsNG8JXX6WtObiXtw8vvPod4QOO0DO2JkNvLafKuwX5efowp8UrRFYiiZkQ8VRp2JZfR59jTsnX2e1znsrzmjLknfpcuXDK6tDc2snLJ1lzdA2dK7vuxrVbN8iZEyZOzMBC9u83bSLlhtt9KQVdu8KBA9CvH3z8sUnQ5s41PRe4wsKFEBUFnTq5Zn1COMnFi9CyJRQsaA4hv3R2Spy3SCm+GP4nu1sspFRMbp449D4t+xcmfNtyxwYsRBYjiZkQCSgvL9q/MIJ9Q08xxLsxY/RmKnx6DzMm9kHH3jVkmgDmhM3B28ubtpXaumydefKY8cy++cZ085wuNptZUIsWjgxNOEOePKbZaWgoVKliHop59FHzt7PZbNCgAdx7r/PXJYST3LxpWgSfP28emcyXL+PLDKz3JCtHn2PuvW8Q6n2BKoseY/DQelw+fyLjCxciC5LETIgk5AwoyLBha9jXZT31o4vR9fQEGg3Ix861NqtDczu2UBvN72tO/hz5Xbrevn3NwKjpes5Ma3PD3bo1ZM/u8NiEk5QrZ+4qlyyBv/823df36weXLjlnfRcvwooVUqsqPJrWZqjAjRth/ny4/37HLVt5edHu+c/Z9+4Zhno/zJdsocJnJfnxq1ekMFOINJLETIgUlK7yIHNG/8PKyp9xwSuK2qu68MqblTl3/IDVobmF4xHH2Xh8o0ubMcYpXx6aNzdjEqe5VdvevbBvn+n9T3ieli3NNvz4Y9N7QfnyMHmy6RnGkebPN934J9dtnRBu7vPPYcoUc4g89JBz1pHDPz/vvreK/U9toGF0MZ49M5EHB+Tlz1XTnbNCITIhScyESKWmHQax85MLjM7Vlp+8wig/oRIhIzsRfTO97egyh9lhs8nmnY3WFVvfnjZhwgTKlClD9uzZqV27NuvXr0/y8/PmzePRRx+lUKFC+Pv7U79+fZYvT/1zCsHBprf7LVvSGLjNZtryPPpoGj/oWdKyLeLbuHEjPj4+1KhRw7kBZkS2bPDmm6Z7/ebNzcDUdevC5s2OW8fMmdCoERQv7pDFpXV7REVF8fbbb1OqVCmyZcvGfffdx5QpUxwSS1aX1m0xffp0qlevTs6cOSlWrBjPP/8858+fd1G06TdnDgweDO++a3phdLZSlRsye/Q//FZ1JBFeN6mz7hleGlQp2cLMdevW0apVK4oXL45SigULFqS4nrVr11K7dm2yZ89O2bJlmZihB46FcBNaa3nJy9NfSRo/fryuVKmSLl++vAZ0REREcrOn2ukje3WP18tr9R66Wr/ses2CMQ5Zrid64NsHdJuZbW7/PXPmTO3r66u//fZbHRYWpl977TWdK1cu/ffffyf6+ddee01/9tlneuvWrTo8PFwPGTJE+/r66j///DNV64+J0fq++7R++uk0BB0baz70wgtp+JDnSeu2iHPp0iVdtmxZ3axZM129enXXBOsIGzZoXbOm1qB1t25anzyZseWdOaO1t7fWX33lkPDSsz2efPJJXbduXb1y5Up95MgRvWXLFr1x40aHxJOVpXVbrF+/Xnt5eekvv/xS//XXX3r9+vW6cuXKuk2bNonO7y62bNE6e3atn3rKnPZc7eb1q/rLz9rpgMHovIOVHvd5B30r6vpd8y1dulS//fbbeu7cuRrQ8+fPT3a5f/31l86ZM6d+7bXXdFhYmP7222+1r6+vnjNnjtbW35PIS17pflkegLzk5YBXiiIiIhyamMXZumKqrtsvl2YYunP/kvpY2GaHLt/dHb5wWDMM/dOen25Pe+CBB/TLL798x3wVK1bUgwcPTvVyAwMD9fvvv5/q+UeP1trXNw334du2aQ1aL1+e6nV4ovRui86dO+uhQ4fq9957z7MSM621jo7W+uuvtS5QQOs8ebQeOVLrqKj0Leurr0xiduaMQ0JL6/ZYtmyZDggI0OfPn3fI+sV/0rotRowYocuWLXvHtLFjx+oSJUo4LcaM+vtvrYsU0bp+fa2v350LudSZv8P0i69X0Oo9dNX+2fXq+V8kOW9qErNBgwbpihUr3jHtpZde0vXq1dPa+nsSeckr3S9pyihEBgQ92o3fR15iasEXWeP7LxV/rMfHHz7KjSuXrA7NJWaFziKnb05alW8FwM2bN9m+fTvNmjW7Y75mzZrx+++/p2qZsbGxXL58mfz5U9+RyPPPg68vfP11Kj9gs5n+oh95JNXr8DTp3Rbfffcdhw8f5r333nN2iM7h7Q29epnmjd26waBBZiTyNDSPvc1mM/tIoYyPzZee7bFo0SLq1KnD559/zj333EP58uV5/fXXue7qcdwymfRsiwYNGvDPP/+wdOlStNacPn2aOXPm0LJlS1eEnGaRkfDEE2ZIkYULre/fqNC9lfh2xH62Pvg9ubQPD+/qT+cB93IsbFO6lrdp06a7tl/z5s3Ztm0bSilfR8QshBUkMRMig7y8fejW51vCX/+b3tRh2K1fCRxWiIXT3sr0PVLZQm08Uf4JcvmZwZnPnTtHTEwMRYoUuWO+IkWKcOpU6saCGzVqFFevXqVTGsaMypsXnnvOjGl282YKM2ttunHs0AF8fFK9Dk+Tnm1x8OBBBg8ezPTp0/Hx9N8mf37TK8yOHVC0KDz2GLRpA3/9lbrPnzwJa9c6rDfG9GyPv/76iw0bNrB3717mz5/PmDFjmDNnDn369HFITFlVerZFgwYNmD59Op07d8bPz4+iRYuSN29exo0b54qQ0yQ62vRpdOyY6bzUAeUKDlOn6XNsHHmR7wv2ZJ3Pv1Sc0YAPP/hfmgszT506lej2i46OBijosICFcDFJzIRwEP9CJRjx8R/seeJnykUH0ObIcFoMLMz+rUutDs0pws+Hs/PUzkR7Y1RK3fG31vquaYn56aefGDZsGDabjcKFC6cpnr594dQpM2hqsrZsMd2sZ5Huz1O7LWJiYnj66ad5//33KV++vKvCc75q1WD1alP7tX07BAbC0KFw9Wryn5s929S+tXXs2HxpOTZiY2NRSjF9+nQeeOABHn/8cUaPHs3UqVOl1swB0rItwsLCePXVV3n33XfZvn07v/zyC0eOHOHll192Rahp0r+/GeFhzhyzu7sbL28fnuvzDQfe+Ju+OogPo1cROKwQC74fkqbCzMS2X9x/HRetEK4liZkQDlbxgcf5ZfQZFpZ9i4PeEVRd0pLX365D5Nl/rA7NoWx7beT2y02L+/8bnLlgwYJ4e3vfVep85syZu0o371qezUaPHj2YNWsWTZs2TXM8lSubVmcpFmDPnGlqUJzVZ7SbSOu2uHz5Mtu2baNv3774+Pjg4+PDBx98wK5du/Dx8WHVqlWuCt3xlIJOnWD/fnjjDRg5EipWNMmaTuIezmYzPT2moUltctJzbBQrVox77rmHgICA29MqVaqE1pp//slc5xNXSs+2GD58OA0bNuSNN96gWrVqNG/enAkTJjBlyhROnjzpirBTZdw4GD8evvoK0nEadSn/QiX4/OOt7Gm1lAoxeWl79FMeG1CYMqmo7ypatGii289e0+/+XWUKkQRJzIRwAuXlxZPPfkzosLMM823KV2yn/MhSfB/Sk9iYaKvDc4iZoTNpXaE1OXxz3J7m5+dH7dq1Wbly5R3zrly5kgYNGiS5rJ9++onu3bszY8aMDD2zERwMmzaZipFExcaampCOHU1tSCaW1m3h7+/Pnj172Llz5+3Xyy+/TIUKFdi5cyd169Z1VejOkysXfPghhIVBnTqmvdfDD8Pu3XfOd/w4/P67Q2tV03NsNGzYkBMnTnDlypXb08LDw/Hy8qJEiRIOiy2rSc+2uHbtGl5ed94yedvPITqp5N7Ffv7ZjLX++uvQs6fV0aRehaAWLB11mkVlh3LIJ4Ljr8DiFQOJOHMsyc/Ur1//ru23YsUK6tSpg9b6lrNjFsJprO59RF7ycsArRc7qlTG1joVt1p37l9QMQ9ftl0tvXTHVkjgcZc/pPZph6MUHFt/1Xlw31JMnT9ZhYWG6X79+OleuXPro0aNaa60HDx6sn3322dvzz5gxQ/v4+OiQkBB98uTJ269Lly6lOa7oaK1LlTI9pSdq7VqtwXSrngWkdVsk5JG9MqbF8uVaV6igtZeX1n36aB3X++HIkVpny6a1g88Xad0ely9f1iVKlNAdOnTQoaGheu3atbpcuXL6xRdfdGhcWVFat8V3332nfXx89IQJE/Thw4f1hg0bdJ06dfQDDzxg1Ve4w86dWufOrXWbNuY86GkuX76sd+zYobdsXKsbP4jO+Ra6yBtKjxveWcdE37pre8R1l9+/f38dFhamJ0+eLN3lyytTvCwPQF7ycsArRVYnZnHWLBijq/bPrtV76B6vl9enj+y1NJ70GvrbUJ3307w6KjrxbshDQkJ0qVKltJ+fn65Vq5Zeu3bt7fe6deumGzdufPvvxo0ba8wzAXe8uiWZXSXvs8/MPXWiPZz37q11yZJm8LMsIi3bIqFMn5hpbbrSHzXKdK1foIDpIr9OHXOH6wRp3R779u3TTZs21Tly5NAlSpTQAwYM0NeuXXNKbFlNWrfF2LFjdWBgoM6RI4cuVqyY7tq1q/7nn39cHPXd/v1X6xIltK5VS+srV6yOJn1Wr159x/m/qD+6fns0w9AP9MulX+zY4K7tsWbNGl2zZk3t5+enS5curb/6b7xBq+9J5CWvdL+U1u5RBS9EBqS4E0dGRhIQEEBERAT+/v6uiClJ0Tdv8PXY53jn/Bxi0QzL24Y+r03HN3tOS+NKLa01FcZX4MF7H2RK6ylWh3OX8+ehRAl45x146614b0RHwz33wLPPmmeMhIjv1CkYMgSmTjV///STaeoohBu7ehUaNza779atULy41RE51vrF4wle9Qa78t6gx5XyfPvZPpRXik/hpNzTlBBuSp4xE8LFfPyy0+f1WYT33keXmEoMuL6AGm/l57e5I6wOLVV2ntrJwQsHE+2N0R0UKABdu5qH36PjP863Zg2cOZNlemMUaVS0KHz3HWzebDoIad3a6oiESFZsLDzzjOnTZsmSzJeUATzUqi/bP49gQu4uFMqePzVJmRAeTWrMRGbgUTVmCe1Y8xOvLniZDfkiaRdRnFEv2Chd5UGrw0rSmyvfZPKOyZwceBJfb/ccx3PXLqhRw/Tz0aGDfWLPnrBqFRw6ZHrpE0IIDzZoEIwaZQaQfuIJq6NxK3KCFx5Lih5EphYSEkJgYCBBQUFWh5Kkmk2eYt3oi0wv0pvNPqeoNPMhhg1rwrWIc1aHdhetNbPCZtEhsIPbJmUA1aub3vBvd51/6xbMm2dqyyQpE0J4uG+/hREjYPRoScqEyEykxkxkBh5dYxbflQun+PiLtoxmM0VveDO6Un/adf/MbZpvbPlnC/Um12PVc6t4uMzDVoeTrNmzzdBVu3ZBtX+XweOPw86dJmsTQggP9euv0KIF9OplxiyTsqa7yC8iPJZ73O0JIQDInb8owz/cxN42K6gWXYAOx0fSdEBBQn9faHVoANhCbRTNXZRGpRpZHUqK2rQxfX2MG4cZLLhCBahWzeqwhBAi3fbtM82zmzaFL7+UpEyIzEYSMyHcULnaj7L4i9P8XG4Yx72vUH15G14bUoNLp/+2LKZYHcus0Fl0qNQBby/3H5zZ1xdeeQWmT9dcmLdGmjEKITza2bPQsiWULGnKmnx8rI5ICOFokpgJ4cYef/o99n50gU/8WjDFaxflvijDpC+7EXPrpstj2XhsI/9e/pcuVTynC/GePSHmViyTL3eU3hiFEB7rxg3TCuDqVdMDoxu3yBdCZIAkZkK4Ob8cuRn09lIOdN9Oi+gy9Lw0jbqD8rFp2TcujcMWaqOEfwnql6zv0vVmROHC0KXERib49iOmQqDV4QghRJppDS+8AH/+CYsWQalSVkckhHAWScyE8BDFy9Vi2sjDbKgzAY2mwdaXeO71+zh5eKfT1x0TG8OcsDl0CuyEl/Kg08a1awSfepujt+5hyRKrgxFCiLR7/30z3vm0aVC3rtXRCCGcyYPusIQQAA1bvsLWzy/xTcCzLPM5QvkpNRnxyRPcvH7Faetc+/daTl89TecqHtYccOlS6tzYQL0aN/7rOl8IITzE9OkmMfvkE+jY0epohBDOJomZEB7I29ePnv2mEf7aYZ6PqcaQqJ+pOjQ/v8z8yCnrs+21USZvGYKKu+94cImaORNq1SL4jez89huEhVkdkBBCpM6GDaYJY/fuMHiw1dEIIVxBEjMhPFi+YmUY++kudjw6h+IxuWhx4B2e7F+UQzt/c9g6bsXcYs6+OXSu3BnlSb0aXr4MP/8MnTvToQMULWrG/BFCCHd3+LDp7KN+ffj6a+lQVoisQhIzITKBqg+2Z9Xo88wq0Z+dPueoPLcpb7/bkCsXT2d42b8d+Y0L1y94XjPGxYtNV2adOuHnBy+9ZJ7RiIiwOjAhhEjaxYumW/z8+WHePPDzszoiIYSrSGImRCahvLzo2GM0+4eeYrB3I0bp36n4SXF++iYYHRub7uXaQm2UL1Ce6kWqOzBaF7DZoF49KF0aMIlZVBR89521YQkhRFJu3oT27c2YZT//bJIzIUTWIYmZyNRCQkIIDAwkKMjDno3KgJwBBXl/2Fr2dVxL3ZiiPH1yPI0H5GPXullpXlZUdBTz9833vGaMly7BsmV3jF1WrJh5eD4kBDKQpwohhFNoDb17m2fL5s+HcuWsjkgI4WqSmIlMrU+fPoSFhfHHH39YHYrLlanWiLmj/2VF4HDOed2g1m+d6f1mFc7/czDVy1h+eDkRURF0ruxhzRgXLIDo6Lu6MQsOhkOH4JdfrAlLCCGSMmIETJ4MkyZBo0ZWRyOEsILSWlsdgxAZleJOHBkZSUBAABEREfj7+7siJrdy68Y1Qr7synsRC/CJVXxUqDO9Xv0eb9/kH17oOq8ru0/vZs8re1wUqYO0aAHXrsHatXdM1hqCgqBQIVOhJoQQ7mDePNOEcehQ+PBDq6PxeB7UvEOIO0mNmRBZgG/2nPR7cz7hPffQJuZ+el+ZSe1BAaxblPTgXtdvXWfRgUWeV1t2/jz8+usdzRjjKGVqzX75BcLDLYhNCCES+OMPeOYZc8p6/32roxFCWEkSMyGykCJlqjB5RDhb6k8hG9403vEqTw8szT8H7m7qufTgUq7cvOJ5idm8eeYhsvbtE327c2coWNA8ayaEEFY6dgyefBKqVzcdE3nJXZkQWZqcAoTIgh5o9jybRl7iuwI9+M3nGBWmPcAnHzbjxpVLt+exhdqoWbQm5Qp42BPoM2fCww9DkSKJvp09O/TqZW6CLl92cWxCCGEXGQlPPGHOSQsWQI4cVkckhLCaJGZCZFFe3j507zuJ8IFHeVnX5r1bK6k8rBCLfhjKlajLLAlf4nm1ZadPw5o1iTZjjO+VV8wjaNOmuSYsIYSILzoaunSBv/+GJUuSLEcSQmQxkpgJkcUFFL6XUZ9sY3fLxdwX7U/rvz4m6L1iXI++7nmDSs+ZY9oCtWuX7GwlSkDbtjB+vOkQRAghXGnAAFixAmbPhsqVrY5GCOEuJDETLqWUyqeU+kEpFWF//aCUypvM/L5Kqc+UUnuUUleVUieUUtOUUsVdGHaWUKnuEywffZYFZYYQFR1Fk6OK0p9MMO1tPIXNBo8+CgUKpDhrcDDs32/6CRFCCFcZN868QkKgWTOroxFCuBPpLl+4lFJqGVAC6GWf9A1wVGvdKon5A4A5wLfALiAfMAbw0VrXsc8m3eU7WMy1q0SPGkG24Z9DQAB8/jl07ereT6b/8w/ce695eKxbtxRn1xpq1IBSpWDRIueHJ4QQS5dCq1bQrx+MGmV1NJmWdJcvPJYkZsJllFKVgDCgntZ6i31aPWATUFFrfSCVywkCtgKltNbHkMTMeY4fhzfeMDVR9eubYt7ata2OKnFffAGDB8OZMyaZTIVJk0xHIIcOQdmyTo5PCJGl7d4NDRvCI4+YzmO9va2OKNOSxEx4LDcu/haZUH0gIi4pA9BabwYigAZpWE4AJhm7lNQMUVFRREZG3vES6VCypOnlcM0auHLFjM7cqxecPWt1ZHez2eCxx1KdlAE8/TTkzQsTJjgvLCGEOHnS9MBYrhxMny5JmRAicZKYCVcqCpxJZPoZ+3spUkplBz4FZmitk8y2hg8fTkBAwO1XyZIl0xWwsGvcGP7809SYzZlj7i7GjjVdi7mDo0dhyxbTzVka5MwJL74IkyfD1avOCU0IkbVdu2bGKouJgcWLIXduqyMSQrgrScxEhimlhimldAqv5J4HU0lMT7geX2AmZr/tndy8Q4YMISIi4vbr+PHjaf1aIiEfH+jTB8LDTQLUr595SGvVKqsjg1mzzCBArRJ9VDFZvXub/k2mT3dCXEKILC02Fp59FsLCTLf499xjdURCCHcmiZlwhPFApRRee4FTQGKjtRQCTie3AntSNgsoAzyaXG0ZQLZs2fD397/jJRykYEGYOBG2bzftAP/3P+jY0QzIYxWbDVq2TFdRdOnSJp8bN066zhdCONZbb8H8+fDTT1CzptXRCCHcnSRmIsO01ue01vtTeN3AdPIRoJR6IO6zSqm6mGfGfk9q+fGSsnJAU631eSd/JZEaNWvC+vXw44+wcSNUqgQffADXr7s2joMHTTPLFAaVTk5wMOzdC2vXOjAuIUSWNnkyfPaZ6X3xySetjkYI4QkkMRMuo7XeB/wCfKuUqmfvkfFbYEn8HhmVUvuVUm3t//fBdJdfB+gKeCulitpffq7/FuIOSplu9A8cgNdeg48+MgnavHmuq36y2UxN2eOPp3sRjzwCgYGm1kwIITJq1Sp4+WXz6tfP6miEEJ5CEjPhal2BPcAK+2s38GyCeSpgatHAjHn2pP3fncDJeK+09OQonClPHhg+HEJDoUoVaN/ejJwaFub8ddtspjg6Z850L0Ip6NsXFiyAY8ccF5oQIuvZv9+cAh95xPSRpKTzdiFEKsk4ZiIzkHHM3M3PP5ti4iNHTDvBYcPS1I19qoWFQeXKsHBhhtsKXbkCJUrAK6+YHFMIIdLq7FmoVw+yZ4fff3fOaU+kSFJh4bGkxkwI4XgtW5qHtj7+GL79FsqXhylTTBdljmSzmTuf5s0zvKjcueH55024rn5MTgjh+W7cgLZtTSHPzz9LUiaESDtJzIQQzpEtG7z5pulev1kz6NHDFCVv2ZLyZ1NDazP4dZs2Zl0O0KcPXLhgFiuEEKmltTnFbdtmKvBLl7Y6IiGEJ5LETAjhXMWLww8/wIYNZkDqevVM1dSpUxlb7q5d/42p5iD33w8tWkjX+UKItPngA5gxA77/3pzihBAiPSQxE0K4RsOG8Mcf8PXXsHixad44ahTcvJm+5dlsUKCAGUfNgYKDYccO83yIEEKkZPp08xjtRx9laNQOIYSQxEwI4ULe3tCrl6np6tYNBg2C6tVhxYq0LUdrk5i1awe+vg4NsVkzKFdOus4XQqRs40Z44QV47jkzmLQQQmSEJGYiUwsJCSEwMJCgoCCrQxHx5c9vMp8dO6BIEdN5R5s28Ndfqfv8tm2mx0cnFE97eZmu8+fOhRMnHL54IUQmcfiwOW3VqwfffCPd4gshMk66yxeZgXSX78m0htmzYeBA09f0G2/A4MGQK1fSnxk40LQf+vdfUwvnYJGRcM890L+/eXZECCHiu3gRGjQwj81u3mxaVQu3ISmy8FhSYyaEsJZS0KmTGZX1jTdgxAioWNE0VUys4Cg2FmbNgg4dnJKUAfj7m5aWX38NUVFOWYUQwkPdumVOP6dPm27xJSkTQjiKJGZCCPeQKxd8+KEZNLpOHdPb4sMPw+7dd863aRP884/Tn7Lv2xfOnDGVeUIIAaas6JVXYP16mD/f9GEkhBCOIomZEMK9lC1r7niWLzdd6tesabKkCxfM+zabaWfYsKFTw6hYER59VDoBEUL8Z+RImDzZDETfuLHV0QghMhtJzIQQ7qlZM1NbNmIETJtmiqYnTjRVWB07ml46nCw4GLZuNS8hRNY2bx68+Sa8/bZp6iyEEI4mnX+IzEA6/8jsTp2CIUNg6lTz9+bNULeu01cbE2O6zm/Y0IyRLYTImrZtg0aNoFUr+Oknl5QLifSTzj+Ex5LETGQGkphlFVu2mIc7Bg50Wd/Uo0aZnPD4cdOzvxAiazl+HB54AEqVgtWrIUcOqyMSKZDETHgsScxEZiCJmXCaixehRAnTg/8771gdjRDClS5fhgcfhIgIUy4khTMeQRIz4bGkMl4IIZKRLx8884x5vO3WLaujEUK4SnS06Rz26FHTLb4kZUIIZ5PETAghUtC3L5w4YR7+F0JkDQMHms5hZ82CypWtjkYIkRVIYiaEECmoWhWaNJGu84XIKsaPh7FjzTHfvLnV0QghsgpJzESmFhISQmBgIEFBQVaHIjxccDBs3Ag7dlgdiRDCmZYuhddeg/79zWDSQgjhKtL5h8gMpPMP4XTR0Wbs66ZNYcoUq6MRQjjD7t1meIyHHzbj3Ht7Wx2RSAfp/EN4LKkxE0KIVPDxgd69YcYMOHfO6miEEI526hQ88QTcf785ziUpE0K4miRmQgiRSi++aP6dNMnaOIQQjnXtGjz5pBlUfvFiyJ3b6oiEEFmRJGZCCJFKBQvC00/DhAmmaaMQwvPFxsJzz0FoqEnKSpSwOiIhRFYliZkQQqRBcDAcPw6LFlkdiRDCEd5+2wyFMWMG1KpldTRCiKxMEjMhhEiDmjVN5wDSdb4Qnm/KFPj0Uxg5Elq3tjoaIURWJ4mZEEKkUXAwrFkDe/ZYHYkQIr1Wr4aXXjKv/v2tjkYIISQxE0KINGvXDooXN4PQCiE8z4ED5jh++GFT+62kg3UhhBuQxEwIIdLI1xdefhl+/BEuXrQ6GiFEWpw7By1bmsKVWbPM8SyEEO5AEjMhhEiHXr1Mz4wy2LQQniMqCtq2hchIWLIE8ua1OiIhhPiPJGZCCJEORYpAp04QEmLGPhJCuDetzViEf/wBCxdCmTJWRySEEHeSxEwIIdIpOBiOHIGlS62ORAiRko8+Ms2Pv/8e6te3OhohhLibJGYiUwsJCSEwMJCgoCCrQxGZ0AMPmJd0nS+Ee/vpJ3j3XfjwQ+jc2epohBAicUprbXUMQmRUijtxZGQkAQEBRERE4O/v74qYRBbx44/w7LOwbx9UrGh1NEKIhH7/HR55xDQ9/v576YExC5AtLDyW1JgJIUQGdOwIhQtL1/lCuKO//jIDRz/wAHz7rSRlQgj3JomZEEJkQLZsZoDa7783Pb0JzxYTA7t2WR2FcIRLl0y3+Hnzwvz55lgVQgh3JomZEEJk0Msvw40bMHWq1ZGIjOrXD2rUgFWrrI5EZMStW9ChA5w6ZbrFL1DA6oiEECJlkpgJIUQGFS8O7dub5oyxsVZHI9Jr3DizDfPmlQ5dPJnW0Ls3rF1rasoqVLA6IiGESB1JzIQQwgGCg+HgQVixwupIRHr8/LOpLRs4ED79FBYtgr//tjoqkR6jRsGkSfDNN9CkidXRCCFE6kmvjCIzkF4ZheW0htq1oVgxc5MvPMeuXfDgg/C//8HcuaZZaokS0KsXfPaZ1dGJtJg/39ReDx4Mn3xidTTCItLFi/BYUmMmXEoplU8p9YNSKsL++kEplTcNn/9aKaWVUv2cF6UQaaeUqTVbtgwOHbI6GpFaJ0/CE09AuXIwfTp4e0OuXPDCC6bW5fp1qyMUqbV9O3TtahKzjz6yOhohhEg7ScyEq80AagCP2V81gB9S80GlVBugLnDCOaEJkTFdukD+/BASYnUkIjWuXoVWrUxt5+LFJiGL06cPXLwIM2ZYF59IvePHzbasWhWmTQMvubsRQnggOXUJl1FKVcIkYy9qrTdprTcBPYEnlFLJPp6tlLoHGA90BW45PVgh0iFHDujZE6ZMgStXrI5GJCc21gwMvn+/ScruuefO98uWNV2tjx1rEjfhvi5fNkmZry8sXGiOQyGE8ESSmAlXqg9EaK23xE3QWm8GIoAGSX1IKeWFqVUbobUOTc2KoqKiiIyMvOMlhCu88opJyqZNszoSkZwhQ2DBAvjpJ6hZM/F5goNh925Yv96loYk0iImBp54yA0kvWQJFi1odkRBCpJ8kZsKVigJnEpl+xv5eUt4EooGxqV3R8OHDCQgIuP0qWbJk2iIVIp3uvRfatDHdrktNi3uaNAk+/xxGjzY1LUlp2tR0tS5d57uvgQPNc52zZplmjEII4ckkMRMZppQaZu+QI7lXHfvsid2qqiSmo5SqDbwGdNdp6EJ0yJAhRERE3H4dP348zd9LiPQKDoZ9++C336yORCT022+mVvOVV+C115Kf18sL+vY1Pf3JKcT9hITAl1+axPmxx6yORgghMk66yxcZppQqCBRMYbajwNPAaK113gSfvwT011p/l8iy+wGjgfjD9nrb/z6utS6NdJcv3IzWUK2aeU5p4UKroxFx9u2D+vWhbl0zpIGPT8qfuXzZPH8WHAwff+z8GEXq/PKLeQYwOBjGjLE6GuFmpLt84bEkMRMuY+/8Iwyoq7Xeap9WF9gMVNRaH0jkMwWAYgkmL8c8c/ad/TOSmAm388038PLLcPgwlCljdTTi7FmoV890DLFxIwQEpP6zr75qnkU7fhyyZ3dejCJ19uyBhg2hUSNT8OHtbXVEws1IYiY8ljRlFC6jtd4H/AJ8q5Sqp5SqB3wLLImflCml9iul2to/c15rvTf+C9Mr46nEEjkh3EXXrubmf8IEqyMRN25A27amU5YlS9KWlIFpznjuHNhszolPpN6pU2bcubJlTbIsSZkQIjORxEy4WldgD7DC/toNPJtgngpAGm+dhHAvuXJBjx4weTJcu2Z1NFmX1mY7bN9ualdKl077MsqXh+bNzbNM0sjEOteuQevWcOuWGeIgTx6rIxJCCMeSxEy4lNb6gtb6Ga21v/31jNb6UoJ5lNZ6ajLLKK21HuPkUIXIsN694dIlmD7d6kiyrg8+MINEf/+9acqYXsHBJrnbvNlxsYnUi42Fbt1g716TlElHu0KIzEgSMyGEcJKyZU2zK6lpscb06TBsmOm0o1OnjC2rRQu47z7pOt8qQ4fC3Llmm9aubXU0QgjhHJKYCSGEEwUHm84K1q2zOpKsZeNGeOEFU8syZEjGl+flBX36wOzZcPJkxpcnUu+772D4cDP2XJs2VkcjhBDOI4mZEEI4UdOmUKmS1LS40uHD5ga+fn3TO6ZyUB9tzz8Pfn7w9deOWZ5I2Zo10KsX9OxpBpMWQojMTBIzIYRwIqVMr34LFsggxa5w8aJpPpovn2n65ufnuGXnzQvPPWcSs5s3HbdckbjwcGjXDho3NoNJOyrBFkIIdyWJmRBCONlzz5leGr/6yupIMrdbt6BDBzhzxgwgXaCA49fRt6/psn3OHMcvW/zn/HkzgHSRIua39vW1OiIhhHA+ScxEphYSEkJgYCBBQUFWhyKysNy5TTO4b781Y2oJx9MaXnkF1q+HefOgXDnnrKdyZXjkEWma6kxRUWbcuUuXTIKdN6/VEQkhhGsoLV2FCc+X4k4cGRlJQEAAERER+Pv7uyImIe5w8KAZD+u776B7d6ujyXxGjIBBg2DqVNPhhzMtWGAShz/+gDp1nLuurEZrc3zYbLBqFTRoYHVEwgNJo1fhsaTGTAghXKBcOdPlunSd73jz5sGbb8Lbbzs/KQNo1QpKlZJaM2f45BOYNs0UYEhSJoTIaiQxE0IIFwkOhj//hE2brI4k89i2DZ55Bjp2NINJu4K3txk8fOZM8zybcAybzYxX9v778NRTVkcjhBCuJ4mZEEK4SPPmpuZMaloc4/hxU3tVrZppwujlwitajx5mfd9+67p1ZmabNpnazmeegXfesToaIYSwhiRmQgjhInGDFM+ZAydOWB2NZ7t82XSLny0bLFwIOXK4dv0FCkDXrqanzVu3XLvuzObIEWjdGoKCYNIk6RZfCJF1SWImhBAu1L27SSZkkOL0i46GLl3g6FHTa1+RItbEERwM//5rOgMR6RMRYRLsPHlg/nxzbAghRFYliZkQQrhQQIBpsiWDFKffwIGwfDnMmmW6r7dK9erw0EPSNDW9bt0yzwaeOGES7IIFrY5ICCGsJYmZEEK4WN++cPo0zJ5tdSSeZ/x4GDvWJEPNm1sdjak1W78edu2yOhLPorX57VavhrlzoWJFqyMSQgjrSWImhBAuVqkSNG0qNS1ptXQpvPYa9O9vBpN2B23awD33yLZMqzFjTK3xxIlmwG4hhBCSmAkhhCWCg2HLFjNIsUjZ7t3QuTO0bGkGk3YXvr4mSZw+Hc6ftzoaz7BokWmO+uabpndLIYQQhiRmQghhgZYtoXRpqWlJjVOnTAcR998PM2aYccTcSc+eEBsLkydbHYn7+/NPM0ZZ27ZmMGkhhBD/kcRMZGohISEEBgYSFBRkdShC3MHb23Sdb7PJIMXJuXYNnnwSYmJg8WLIndvqiO5WuLDpJXLCBBOnSNy//5px5wID4YcfXDvunBBCeAKltbY6BiEyKsWdODIykoCAACIiIvD393dFTEKk6MIFKFEC3noLhg61Ohr3ExsLnTrBsmWmg41atayOKGnbtplxuObPN8+diTtduQKNGsG5c6YJb7FiVkckMjEZCU94LCmvEkIIi+TPD888I4MUJ+Xtt2HePNN80Z2TMoA6daBePWmampiYGDMY98GDsGSJJGVCCJEUScyEEMJCwcFmHKf5862OxL1MmQKffgojR0Lr1lZHkzrBwbBqFYSGWh2Jexk0yCRkNhtUq2Z1NEII4b6kKaPIDKQpo/BoTZqYWoX1662OxD2sXg3NmsGLL5rntpSHNEy6eRNKlTJNGb/6yupo3MPEiabXynHjzPh9QriAh5wxhLib1JgJIYTFgoNhwwbYudPqSKx34AC0a2fGtho71nOSMgA/P3jpJZg2DS5dsjoa661YYZKx4GBJyoQQIjUkMRNCCIu1bg0lS8rzSefOmWEEiheHWbPMGGGe5qX/t3fvwVFW9x/HP4eLoEICFhSLAR0sl2iBUSJI63hrvQygohYmggIK/ICyjG3Btg4Vp9YC2kqlJEpAUSGQcgfBAW2H/hyKUVoFAbVYWi20IvXGgoQEwtM/vlAQE5OQ3T37PPt+zewEQ/Z5vpidJJ+cc77f/7OVszlzfFfi17Zt0ve+J11/vfTYY76rAYBwIJgBgGeNGtl2r/nzM3dIcXm5zbbat09avVrKzvZd0ak591wLJAUF1lUyE334oQXs9u2lkhJ7fQMAakYwA4A0MGKEFATS7Nm+K0m9ILDzZBs3SitW2ODtMIvFpB07rM1/pikrszN25eXW8KN5c98VAUB4EMwAIA20aiXl51uzi8OHfVeTWr/4hTRvnvTss9ZyPux69ZIuvTTztqYeOSINGyZt3mzDwNu1810RAIQLwQwA0kQsJv3zn/ZDbaZYsEB64AELZwMH+q4mMZyzz+XatdbMJFNMmmRnA4uLba4bAKBuaJePKKBdPiLjW9+SmjSxeVhRt2GDdV8cOFB65plwdWCsycGD1tAlP9+6S0bdc89JQ4ZIjzwiTZjguxpkuAh9JUGmYcUMkVZQUKDc3Fzl5eX5LgWolVjM5nht3eq7kuT6+9/tLFLPnlJRUbRCmSQ1bSqNHGmBc98+39Uk18sv2xnB4cOl8eN9VwMA4cWKGaKAFTNExqFD1s3upptsOG8UffaZ1Lu3/VtLS6Wvfc13Rcmxa5c1MvnNb6I7x+vdd+1MXffu0po14RxxgMiJ2K95kElYMQOANNK4sTRqlDR3rvTpp76rSbxDh6Tbb5d277a2+FENZZJ03nk2AmDGjGi2zv/4Y2uLf/bZ0uLFhDIAqC+CGQCkmZEjLcBEbUhxEEhjxtjWt2XLpI4dfVeUfLGYNQD5/e99V5JYFRXSrbfaLw9WrZJatvRdEQCEH8EMANJMmzbSgAE2pLiy0nc1ifPrX9uctqIi6corfVeTGldcIXXtGq3W+UFgvzwoLZWWL5c6dPBdEQBEA8EMANJQLGYNMqIypHjZMum++6T775eGDvVdTeoca52/erV9PqNg8mSbOTdnjnURBQAkBs0/EAU0/0AkXXaZbRFbu9Z3JfXzl7/YylHfvlJJidQgw34leOCAnTcbNsxWDcNs4UIbb/Dggza3DEhDNP9AaGXYt0cACI9YTHrxxXAPKd65U+rXz7bzPfts5oUySTrjDGsl//TT0uef+67m1JWWSnfdJd1xhw0FBwAkVgZ+iwSAcBgwwDrezZjhu5JTs2+fhbLTTpNWrJBOP913Rf6MGSPF49K8eb4rOTXvvSfdfLPUo4f01FPRmzsHAOmAYIaUcs61dM7Ndc7tPfqY65xrUYvndXHOrTz6nH3OuVLnXLsUlAx406TJ8SHF8bjvauqmslLKz7dzVatWSeec47siv84/30Lqb39rzTPCZO9e24barJmdFWza1HdFABBNBDOk2nxJ3SXdcPTRXdLcr3qCc66DpPWS3pF0laRukh6SdDB5ZQLpYdQo6eBB2wYYJj/6kQ0cXrRIuvhi39Wkh1hM2rZN+uMffVdSe4cP28rtrl3WwKR1a98VAUB00fwDKeOc6yLpLUm9giB49ej7ekl6RVLnIAiqPEnjnCuRdCgIgjurufSXXsTl5eUqLy//33/H43Hl5OTQ/AOhNHCgtGmT9Pbb4TijVVAgjR0rFRZKo0f7riZ9BIGF1E6dpKVLfVdTs2Nz52bPtpB97bW+KwJqhY22CK0QfItHhFwuae+xUCZJQRCUStorqXdVT3DONZDUR9J259xa59we59yrzrlbvupGkydPVnZ29v8eOTk5iftXACkWi0nbt0svveS7kpqtWSONGyfdey+h7GTOWWBdsUJ6/33f1dTs8celJ5+UnniCUAYAqcCKGVLGOXe/pKFBEHQ86f3bJc0JgmByFc9pI+kDSQckTZS0TrYF8peSrg6C4P/FihkiLgikSy6R2ra181rpassWm2t15ZU2eLhhQ98VpZ/9+611/qhR0pQpvqup3vPPW7OPCROkqVN9VwPUCStmCC1WzFBvzrkHnXNBDY8eRz+8qt8EuGreLx1/ja4IgmBaEASbgiCYImmVpFHV1dSkSRNlZWV94QGE1bEhxS+8IO3Y4buaqu3ebQ0iOnSQFiwglFWnWTObZzZrllRW5ruaqr3xhjVu6d/fhkkDAFKDYIZEmCGpSw2PrZJ2S6qqN1trSR9Wc+2PJB2WnU070duS6MqIjJGfL511lp3fSjcHDtjqyqFDttLSrJnvitLb978vffqpBdh0869/WffILl2kuXPDcaYRAKKCL7motyAIPgqC4J0aHgdlTT6ynXOXHXuuc66npGxJG6q5doWkjZI6nfRXHSWF4JQGkBinn358SPH+/b6rOe7IEWnIEGnrVgtl553nu6L0d+GF0o03pl/r/P37LZQ5J61caYOxAQCpQzBDygRB8LakNZJmOed6He3IOEvSqhM7Mjrn3nHO9T/hqY9KGuicG+Gcu9A5N1ZSP0mFqawf8G30aBvanE5DiidOlJYskYqLpUsv9V1NeMRi1mnzT3/yXYmprJQGDZLefdfOMZ57ru+KACDzEMyQaoMkbZH04tHHm5JOboPfSbaKJkkKgmCZ7DzZfUefO1zSbUEQrE9FwUC6aN/etgzOmJEeKy1z5tgZpEcflW65xXc14XLdddI3vmGrZungxz+2QFZSInXr5rsaAMhMdGVEFNT4Io7H48rOzqYrI0Jv3TrpmmukP/zB3vqs47rrpLvvtpbqjj5odTZ9uvTDH1rr/LZt/dUxc6Z1iZw+3VbygJDjqxFCi2CGKCCYIWMEgdS1q51TWrbMTw1//at0+eW2dfGFF6TGjf3UEXbxuAWye++VHnrITw0vvWTn3UaPTp/VO6CeCGYILbYyAkCIHBtSvHKlnyHFH30k9ekjtWkjLVpEKKuPrCxrnFJUJJ0wdjFl3npLuv12W/mcNi319wcAfBHBDABCZvBg+6G+MMXtb8rLpVtvtZWe1aulFi1Se/8oGjtW2rNHWrgwtffds8cCdvv2dq6sUaPU3h8A8GUEMwAImTPPtLNds2enbkhxEEgjRkivvSYtXy5dcEFq7ht1nTtL3/1uarcRlpVZE5myMmv4we5uAEgPBDNEWkFBgXJzc5WXl+e7FCChjg0pnj8/Nfd7+GEbODxnjtS7d2rumSliMWnjRunVV5N/ryNHpGHDpM2bbe5cu3bJvycAoHZo/oEooPkHMlK/ftLOndIbbyS3K2JJiZSfL/3859LPfpa8+2Sqykprnd+7d/Jn1D3wgDUaWbLEtqUCEUTzD4QWK2YAEFKxmK18rE/iRL9XXpGGDrVzbRMnJu8+maxhQ1sBXbhQ2r07efd57jkLZVOnEsoAIB0RzAAgpL7zHalTp+SdT/rHP+wsUl6enWdjVlny3H23dbgsKkrO9V9+WRo+XLrnHmnChOTcAwBQPwQzAAipBg2sq9/SpdKuXYm99mefWde+rCybl9akSWKvjy9q2dJWJZ98UqqoSOy1//Y3qX9/6YorrJMnARsA0hPBDABCbMgQ6Ywz7Af6RDl0SBowQPrgA2uL36pV4q6N6o0da//Ply5N3DU/+cQCduvW0uLF0mmnJe7aAIDEIpgBQIg1b25nwIqKpIMH63+9ILCza+vWWUDo1Kn+10TtfPOb0lVXJW5rakWFdNtt0scfW8Bu2TIx1wUAJAfBDABCbuxY6T//ScyQ4mnTpJkzLehdfXX9r4e6icWkDRuk11+v33WCQBo1yq61fLnUoUNCygMAJBHBDABCrmNH6frrbaWlPhNQVqyQxo+XfvITm3WF1LvpJiknp/6rZlOn2sy5p5+Wvv3txNQGAEgughkAREAsJv35z6c+pPj116U77rA26g8/nNjaUHuNGkljxkgLFtgq6KlYvFj66U+lSZOkQYMSWx8AIHkIZgAQATfeaNvVTmWlZdcuG1Z90UU266oB3xm8Gj7c3s6eXffnvvaadOedNhB80qTE1gUASC6+/QJABDRoYEOKFy2q25Di/fstlDVqJK1caR0e4VerVrbS9cQT0uHDtX/e++/bVshLLrEtjLTFB4BwIZgBQEQMG2ZDimfOrN3HV1ba9sUdO6RVq6Q2bZJbH2ovFpN27rRzf7URj0t9+1qwXr5cato0qeUBAJKAYIZIKygoUG5urvLy8nyXAiRdixbSXXfVfkjx+PHWRv13v7NW7Ugf3btb047p02v+2MOHpYEDLcitXm0zywAA4eOC+rTwAtJDjS/ieDyu7Oxs7d27V1lZWamoCfBi2zbp4oul4mJbDatOYaFtfZwxw94i/SxcaIFr82apa9eqPyYI7PM3a5a0Zo107bWprRFIQ2ziRWixYgYAEXLRRdI113x1E5A1a6Rx4+xBKEtf/ftLbdt+9edy+nQ7i1ZYSCgDgLAjmAFAxIwbJ5WWWvv8k23dKg0YIN1wg/TYY6mvDbXXuLENiS4ulj755Mt///zz0g9+IE2YII0Ykfr6AACJRTADgIjp21c6//wvr7Ts3i316SNdcIHNyWrY0Et5qIORI61Jy1NPffH9mzZZS/ybb5amTPFSGgAgwQhmABAxDRvakOKSEmnPHntfWZn9EF9RYR0Ymzf3WyNq5+yz7ZxZYaEFNEn6978tfHfuLM2bx9w5AIgKvpwDQATdc48FtFmzpCNHpCFDpC1bbPtbTo7v6lAXsZj03nsWqD//3ObOOWdz584803d1AIBEaeS7AABA4p111vEhxfG4tHixtGSJ1KOH78pQV3l5Us+e0uOPS888I23fLq1fL339674rAwAkEu3yEQW0yweq8OabUrdu9udHHrEmEQin4mJp8GDbtrhypZ0VBFAl2uUjtAhmiAKCGVCNwYOlc86RfvUr2/6GcKqosDEIgwdbp0YA1eIrHUKLYIYoIJgBAACJYIYQo/kHAAAAAHhGMEOkFRQUKDc3V3l5eb5LAQAAAKrFVkZEAVsZAQCAxFZGhBgrZgAAAADgGcEMAAAAADxjKyMygnMuS9JeSdlBEMR91wMAAACciGCGjOCcc5KaS9oX8KIHAABAmiGYAQAAAIBnnDEDAAAAAM8IZgAAAADgGcEMAAAAADwjmAEAAACAZwQzAAAAAPCMYAYAAAAAnhHMAAAAAMCz/wKZaAaAin5c1gAAAABJRU5ErkJggg==\n", "text/plain": [ "Graphics object consisting of 3 graphics primitives" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plot0 = list_plot([(i/n,l0[i]) for i in range(n+1)],\n", " plotjoined=True,color='red',legend_label='0',\n", " title='Normalized number of regions of each type between indices 0 and cn, for a random word taken uniformly, with with n='+str(n),\n", " axes_labels=['c',''])\n", "plot1 = list_plot([(i/n,l1[i]) for i in range(n+1)],\n", " plotjoined=True,color='blue',legend_label='1')\n", "plot2 = list_plot([(i/n,l2[i]) for i in range(n+1)],\n", " plotjoined=True,color='green',legend_label='2')\n", "(plot0+plot1+plot2).show()\n", "\n", "# In the limit, the green and red curves should look like the same brownian motion W_c,\n", "# while the blue one looks like -2*W_c." ] } ], "metadata": { "kernelspec": { "display_name": "SageMath 9.2", "language": "sage", "name": "sagemath" }, "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.9.2" } }, "nbformat": 4, "nbformat_minor": 5 }