Appendix-E-Ensemble-Kalman-Filters.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "[Table of Contents](http://nbviewer.ipython.org/github/rlabbe/Kalman-and-Bayesian-Filters-in-Python/blob/master/table_of_contents.ipynb)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Ensemble Kalman Filters"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>\n",
       "@import url('http://fonts.googleapis.com/css?family=Source+Code+Pro');\n",
       "@import url('http://fonts.googleapis.com/css?family=Lora');\n",
       "\n",
       "//@import url('http://fonts.googleapis.com/css?family=Open+Sans');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Vollkorn');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Karla');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Poppins');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Arimo');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Roboto');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Lato');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Domine');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Chivo');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Cardo');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Arvo');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Crimson+Text');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Ubuntu');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Fontin');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Raleway');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Merriweather');\n",
       "\n",
       "\n",
       ".CodeMirror pre {\n",
       "    font-family: 'Source Code Pro', Consolas, monocco, monospace;\n",
       "}\n",
       "    div.cell{\n",
       "        //width: 950px;\n",
       "        margin-left: 0% !important;\n",
       "        margin-right: auto;\n",
       "    }\n",
       "    div.text_cell_render{\n",
       "        font-family: 'Lora';\n",
       "        //font-family: 'Open Sans';\n",
       "        //font-family: 'Karla',verdana,arial,sans-serif;\n",
       "        //font-family: 'Roboto',verdana,arial,sans-serif;\n",
       "        //font-family: 'Lato',verdana,arial,sans-serif;\n",
       "        //font-family: 'Domine',verdana,arial,sans-serif;\n",
       "        //font-family: 'Chivo',verdana,arial,sans-serif;\n",
       "        //font-family: 'Cardo',verdana,arial,sans-serif;\n",
       "        //font-family: 'Arvo',verdana,arial,sans-serif;\n",
       "        //font-family: 'Poppins',verdana,arial,sans-serif;\n",
       "        //font-family: 'Ubuntu',verdana,arial,sans-serif;\n",
       "        //font-family: 'Fontin',verdana,arial,sans-serif;\n",
       "        //font-family: 'Raleway',verdana,arial,sans-serif;\n",
       "        //font-family: 'Merriweather',verdana,arial,sans-serif;\n",
       "        //font-family: 'Crimson Text', verdana,arial,sans-serif;\n",
       "        //font-family: verdana,arial,sans-serif;\n",
       "        //font-family: arial,sans-serif;\n",
       "        line-height: 125%;\n",
       "        font-size: 130%;\n",
       "        text-align: justify;\n",
       "        text-justify:inter-word;\n",
       "    }\n",
       "    div.text_cell code {\n",
       "        background: transparent;\n",
       "        color: #000000;\n",
       "        font-weight: 400;\n",
       "        font-size: 12pt;\n",
       "        //font-style: bold;\n",
       "        font-family:  'Source Code Pro', Consolas, monocco, monospace;\n",
       "   }\n",
       "    h1 {\n",
       "        font-family: 'Open sans',verdana,arial,sans-serif;\n",
       "\t}\n",
       "\n",
       "    div.input_area {\n",
       "        background: #F6F6F9;\n",
       "        border: 1px solid #586e75;\n",
       "    }\n",
       "\n",
       "    .text_cell_render h1 {\n",
       "        font-weight: 200;\n",
       "        font-size: 30pt;\n",
       "        line-height: 100%;\n",
       "        color:#c76c0c;\n",
       "        margin-bottom: 0.5em;\n",
       "        margin-top: 1em;\n",
       "        display: block;\n",
       "        white-space: wrap;\n",
       "        text-align: left;\n",
       "    }\n",
       "    h2 {\n",
       "        font-family: 'Open sans',verdana,arial,sans-serif;\n",
       "        text-align: left;\n",
       "    }\n",
       "    .text_cell_render h2 {\n",
       "        font-weight: 200;\n",
       "        font-size: 16pt;\n",
       "        font-style: italic;\n",
       "        line-height: 100%;\n",
       "        color:#c76c0c;\n",
       "        margin-bottom: 0.5em;\n",
       "        margin-top: 1.5em;\n",
       "        display: block;\n",
       "        white-space: wrap;\n",
       "        text-align: left;\n",
       "    }\n",
       "    h3 {\n",
       "        font-family: 'Open sans',verdana,arial,sans-serif;\n",
       "    }\n",
       "    .text_cell_render h3 {\n",
       "        font-weight: 200;\n",
       "        font-size: 14pt;\n",
       "        line-height: 100%;\n",
       "        color:#d77c0c;\n",
       "        margin-bottom: 0.5em;\n",
       "        margin-top: 2em;\n",
       "        display: block;\n",
       "        white-space: wrap;\n",
       "        text-align: left;\n",
       "    }\n",
       "    h4 {\n",
       "        font-family: 'Open sans',verdana,arial,sans-serif;\n",
       "    }\n",
       "    .text_cell_render h4 {\n",
       "        font-weight: 100;\n",
       "        font-size: 14pt;\n",
       "        color:#d77c0c;\n",
       "        margin-bottom: 0.5em;\n",
       "        margin-top: 0.5em;\n",
       "        display: block;\n",
       "        white-space: nowrap;\n",
       "    }\n",
       "    h5 {\n",
       "        font-family: 'Open sans',verdana,arial,sans-serif;\n",
       "    }\n",
       "\n",
       "    .text_cell_render h5 {\n",
       "        font-weight: 200;\n",
       "        font-style: normal;\n",
       "        color: #1d3b84;\n",
       "        font-size: 16pt;\n",
       "        margin-bottom: 0em;\n",
       "        margin-top: 0.5em;\n",
       "        display: block;\n",
       "        white-space: nowrap;\n",
       "    }\n",
       "    div.output_subarea.output_text.output_pyout {\n",
       "        overflow-x: auto;\n",
       "        overflow-y: scroll;\n",
       "        max-height: 50000px;\n",
       "    }\n",
       "    div.output_subarea.output_stream.output_stdout.output_text {\n",
       "        overflow-x: auto;\n",
       "        overflow-y: scroll;\n",
       "        max-height: 50000px;\n",
       "    }\n",
       "    div.output_wrapper{\n",
       "        margin-top:0.2em;\n",
       "        margin-bottom:0.2em;\n",
       "}\n",
       "\n",
       "    code{\n",
       "        font-size: 6pt;\n",
       "\n",
       "    }\n",
       "    .rendered_html code{\n",
       "    background-color: transparent;\n",
       "    }\n",
       "    ul{\n",
       "        margin: 2em;\n",
       "    }\n",
       "    ul li{\n",
       "        padding-left: 0.5em;\n",
       "        margin-bottom: 0.5em;\n",
       "        margin-top: 0.5em;\n",
       "    }\n",
       "    ul li li{\n",
       "        padding-left: 0.2em;\n",
       "        margin-bottom: 0.2em;\n",
       "        margin-top: 0.2em;\n",
       "    }\n",
       "    ol{\n",
       "        margin: 2em;\n",
       "    }\n",
       "    ol li{\n",
       "        padding-left: 0.5em;\n",
       "        margin-bottom: 0.5em;\n",
       "        margin-top: 0.5em;\n",
       "    }\n",
       "    ul li{\n",
       "        padding-left: 0.5em;\n",
       "        margin-bottom: 0.5em;\n",
       "        margin-top: 0.2em;\n",
       "    }\n",
       "    a:link{\n",
       "       color:#447adb;\n",
       "    }\n",
       "    a:visited{\n",
       "       color: #1d3b84;\n",
       "    }\n",
       "    a:hover{\n",
       "       color: #1d3b84;\n",
       "    }\n",
       "    a:focus{\n",
       "       color:#447adb;\n",
       "    }\n",
       "    a:active{\n",
       "       font-weight: bold;\n",
       "       color:#447adb;\n",
       "    }\n",
       "    .rendered_html :link {\n",
       "       text-decoration: underline;\n",
       "    }\n",
       "    .rendered_html :hover {\n",
       "       text-decoration: none;\n",
       "    }\n",
       "    .rendered_html :visited {\n",
       "      text-decoration: none;\n",
       "    }\n",
       "    .rendered_html :focus {\n",
       "      text-decoration: none;\n",
       "    }\n",
       "    .rendered_html :active {\n",
       "      text-decoration: none;\n",
       "    }\n",
       "    .warning{\n",
       "        color: rgb( 240, 20, 20 )\n",
       "    }\n",
       "    hr {\n",
       "      color: #f3f3f3;\n",
       "      background-color: #f3f3f3;\n",
       "      height: 1px;\n",
       "    }\n",
       "    blockquote{\n",
       "      display:block;\n",
       "      background: #fcfcfc;\n",
       "      border-left: 5px solid #c76c0c;\n",
       "      font-family: 'Open sans',verdana,arial,sans-serif;\n",
       "      width:680px;\n",
       "      padding: 10px 10px 10px 10px;\n",
       "      text-align:justify;\n",
       "      text-justify:inter-word;\n",
       "      }\n",
       "      blockquote p {\n",
       "        margin-bottom: 0;\n",
       "        line-height: 125%;\n",
       "        font-size: 100%;\n",
       "      }\n",
       "</style>\n",
       "<script>\n",
       "    MathJax.Hub.Config({\n",
       "                        TeX: {\n",
       "                           extensions: [\"AMSmath.js\"],\n",
       "                           equationNumbers: { autoNumber: \"AMS\", useLabelIds: true}\n",
       "                           },\n",
       "                tex2jax: {\n",
       "                    inlineMath: [ ['$','$'], [\"\\\\(\",\"\\\\)\"] ],\n",
       "                    displayMath: [ ['$$','$$'], [\"\\\\[\",\"\\\\]\"] ]\n",
       "                },\n",
       "                displayAlign: 'center', // Change this to 'center' to center equations.\n",
       "                \"HTML-CSS\": {\n",
       "                    scale:95,\n",
       "                        availableFonts: [],\n",
       "                        preferredFont:null,\n",
       "                        webFont: \"TeX\",\n",
       "                    styles: {'.MathJax_Display': {\"margin\": 4}}\n",
       "                }\n",
       "        });\n",
       "</script>\n"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#format the book\n",
    "%matplotlib inline\n",
    "from __future__ import division, print_function\n",
    "from book_format import load_style\n",
    "load_style()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "> I am not well versed with Ensemble filters. I have implemented one for this book, and made it work, but I have not used one in real life. Different sources use slightly different forms of these equations. If I implement the equations given in the sources the filter does not work. It is possible that I am doing something wrong. However, in various places on the web I have seen comments by people stating that they do the kinds of things I have done in my filter to make it work. In short, I do not understand this topic well, but choose to present my lack of knowledge rather than to gloss over the subject. I hope to master this topic in the future and to author a more definitive chapter. At the end of the chapter I document my current confusion and questions. In any case if I got confused by the sources perhaps you also will, so documenting my confusion can help you avoid the same.\n",
    "\n",
    "\n",
    "The ensemble Kalman filter (EnKF) is very similar to the unscented Kalman filter (UKF) of the last chapter. If you recall, the UKF uses a set of deterministically chosen weighted sigma points passed through nonlinear state and measurement functions. After the sigma points are passed through the function, we find the mean and covariance of the points and use this as the filter's new mean and covariance. It is only an approximation of the true value, and thus suboptimal, but in practice the filter is highly accurate. It has the advantage of often producing more accurate estimates than the EKF does, and also does not require you to analytically derive the linearization of the state and measurement equations. \n",
    "\n",
    "The ensemble Kalman filter works in a similar way, except it uses a *Monte Carlo* method to choose a large numbers of sigma points. It came about from the geophysical sciences as an answer for the very large states and systems needed to model things such as the ocean and atmosphere. There is an interesting article on it's development in weather modeling in *SIAM News* [1]. The filter starts by randomly generating a large number of points distributed about the filter's initial state. This distribution is proportional to the filter's covariance $\\mathbf{P}$. In other words 68% of the points will be within one standard deviation of the mean, 95% percent within two standard deviations, and so on. Let's look at this in two dimensions. We will use `numpy.random.multivariate_normal()` function to randomly create points from a multivariate normal distribution drawn from the mean (5, 3) with the covariance\n",
    "\n",
    "$$\\begin{bmatrix}\n",
    "32 & 15 \\\\ 15 & 40\n",
    "\\end{bmatrix}$$\n",
    "\n",
    "I've drawn the covariance ellipse representing two standard deviations to illustrate how the points are distributed."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhwAAAECCAYAAAC1yg4KAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xl8W+d54Psf9n0nCBLcCZGiJGqzJNpyLNlObMtuNqWN\nbztux+l8uqTzyafJxJ02Tsep02naJL2d3nRa5zadtDd1mjhpZmrZTRzL8RJrMR1qo0SJkkiCO0GQ\nIAAS+0bi/kGBIiVRpjbLlp7vP5ahg4NzHvLoffAuz6soFotFhBBCCCFuIOXNvgAhhBBC3Pok4RBC\nCCHEDScJhxBCCCFuOEk4hBBCCHHDScIhhBBCiBtOEg4hhBBC3HDqq33j3NwcTz31FAMDAyiVSv70\nT/8UrVbLk08+iVKppKmpiaeffvp6XqsQQggh3qeuOuF4/fXXUSgUPPfcc3R0dPDXf/3XFItFnnji\nCbZu3crTTz/Nq6++ygMPPHA9r1cIIYQQ70NXPaTywAMP8Gd/9mcABAIBbDYb3d3dbN26FYCdO3fS\n3t5+fa5SCCGEEO9r1zSHQ6lU8uSTT/KVr3yFj3zkIywuWmoymYjH49d8gUIIIYR4/7vqIZWSr33t\na4TDYT75yU+SzWYXXk8mk1it1ms9vRBCCCFuAVfdw/HCCy/wD//wDwDodDqUSiWtra10dHQAsG/f\nPrZs2XLZc8g2LkIIIcTtQXG1m7el02m++MUvMjU1RaFQ4NOf/jSNjY089dRT5PN5fD4fX/nKV1Ao\nFJc9Tygkwy4r4XZbJFYrIHFaOYnVykicVkbitHK3eqzcbsslX7/qIRWDwcA3vvGNi17/7ne/e7Wn\nFEIIIcQtSgp/CSGEEOKGk4RDCCGEEDecJBxCCCGEuOEk4RBCCCHEDScJhxBCCCFuOEk4hBBCCHHD\nScIhhBBCiBtOEg4hhBBC3HCScAghhBDihpOEQwghhBA3nCQcQgghhLjhJOEQQgghxA0nCYcQQggh\nbjhJOIQQQghxw0nCIYQQQogbThIOIcS7LhCIsn9/iEAgerMvRQjxLpGEQwjxrvP7C6RSFfj9hZt9\nKUKId4kkHEKId53Pp8ZoDOLzqW/2pQgh3iXytAsh3nVerwOv92ZfhRDi3SQ9HEIIIYS44SThEEII\nIcQNJwmHEEIIIW44STiEEEIIccNJwiGEEEKIG04SDiHELUcKiwnx3iMJhxDiliOFxYR475GEQwhx\ny7lZhcWkZ0WI5UnCIYS4pQQCUfz+Aj6fGq/X8a5+tvSsCLE8STiEELeUm9noS8l2IZZ3VU9FoVDg\nj//4jxkbGyOfz/N7v/d7rFq1iieffBKlUklTUxNPP/309b5WIYR4Rz6fGr//5jT6UrJdiOVd1RP5\n4osv4nA4+Mu//EtisRgf//jHaWlp4YknnmDr1q08/fTTvPrqqzzwwAPX+3qFEOKypNEX4r3pqoZU\nHnnkET73uc8BMDs7i0qloru7m61btwKwc+dO2tvbr99VCiGEEOJ97aoSDoPBgNFoJJFI8LnPfY7P\nf/7zFIvFhb83mUzE4/HrdpFCCLHY1awGeTdWkMgqFSGWd9WTRsfHx/nUpz7FJz7xCT784Q+jVJ4/\nVTKZxGq1XpcLFEKIC5UmhnZ0vHMDX0oCOjqi1zSZdCXJhKxSEWJ5VzWHY2pqit/6rd/iT/7kT7jr\nrrsAWLNmDYcOHWLbtm3s27dv4fV34nZbruYSbksSq5WROK3clcRqdDRKb2+WpiYd1dXv7nLTC7W1\nFTh4cJSenhRqtY9IJMHGjZZLXuPJk0k0mlqs1rOUl8doanJc8e+I221ZOE8kEmTjxku/v62tQG/v\n1X3GreB2vOerdTvG6qoSjm9961vEYjG++c1v8swzz6BQKPhv/+2/8ZWvfIV8Po/P5+Phhx9e0blC\nIRl6WQm32yKxWgGJ08pdaaxKPQTRaBCd7savALlcPQ2dTo1Go8blqmV4+Cx33OEmFIpf8hqdzhzR\naD8tLQa8XhNwZf/ulOJUOo/TqV72/TqdmtZW9RV/xq1Anr2Vu9VjtVwypSgunnxxE9zKQb+ebvVf\n0OtF4rRyVxqrdyqoda0Fty58//79IVKpCozGIDt2uFf0eTei6Jf8Tq2MxGnlbvVYLZdwSOEvIQSw\n/ByF0usAO3a4l23Ir3X+woXvv5oiWl6v41wdjoJM3BTiPUYSDiEEsHzCsNJE4lqrbF74fq/XccUJ\nTiAQZc+eECMjOpm4KcR7jNTfFUIAy1foXGnlzsUFt65maONKCnYFAlHC4TgQY/1658Lrfn8Bh6OW\naHSY++67eBhGCHHzSMIhhACWb/CvpnKn319gZERHV1eQ3bu57puo+f0FDIYmjMbgknPPJ0Vx7r13\n+Z4RIcTNIUMqQtwi3ktFp3w+NdPTQRyO2hsytLHc8M2FwzDvpZgIcbuTHg4hbhHn5zQEL9sjURru\naGsr3LDlrV6vg927we+PX/WcjmPHBmlvz7F9u5bNm+svOv+l7vHCoZyVxmSlbsQqGCFuF9LDIcQt\nYqWTNkuNcG9v9oZez+IVI8eODV5xT0N7e47R0Xq+//3Iit93PVa6XMn5hRArJz0cQtwiVjrXojQJ\ntKnpxn9D9/sLjI5aOHp0hDvuqADiS65xuR6DQCCKwZBmevoIra0N+P2FFU1INZuTdHWdZvt2LeBe\nUUyupNdipRNohRAXk6dGiNtMqRF2uy0cPz58zUMEl2uwfT41XV3D1Ne7iUaHuffepStHFg95wPnz\nzCcYd/DQQz24XCp8PvXC5xw50supUwbKyzN85jMtSz4zkTDh8/lIJIIrvsbz+7L04HJdPhZXM4FW\nCDFPhlSEuI1djyGC5ephlIqF7d7tZsMGPbt3XzyZM5MZx+8/jdmcXHKe0lBIW5tzYRJo6e97eiAS\nqWBysvqi6/b51KTTPYTD8SXDMJe7z9JngUKGS4S4gaSHQ4jb2PUYIrjwHKXiW3Z7BZA9lzCw8Hd+\nf4FwOEY6XcXRoyPU1xtob0/T1BTD74+yfbsWr7eeUo8HRBfmg3R09NDcrMRu76OpyYbPV7PkWrxe\nBy7XxRNFL3efpV6LQODiYwKBKCdPJnE6czJJVIhrJAmHELexKx0iuNTQxPx/zycHpeJbQ0NnMRrN\nBALnj+3oiDAwYMNuT5JOD+N0mmlvH+au7Y2c7SnStGoNQ0M9JBIh+voCTE9XEQ5HaGvj3PkVrFnz\nAbZsufT+KnDp5OJS93nhvVzqGL+/gEZTSzTaz+IhH0k+hLhyknAIIVbsUstMSz0aDkct4fAYoMBg\niNHSYsZgaMbvD2Jz5un0B/l5zwBnRmaJFcaZyWZJpLPk54q88uYvADC2v4JRq8GoMUGqSK0jStFi\nxnLWRCFXBcQuuerk4uTh6u7lQj6fmkgkiNOpvu5LbIW43UjCIYRY4sLGe/H/X6r3YL7nwcrRo92s\nXavG692MTh8go5nmpbdfoSc0wfFnJijMXX5jaoUCUrk8qVweSAEwHJziQBD0WjVtDTX82v1N3P2B\nBlTKpdPPSj0n4XCE3buXrna53ITWdxpO8nodbNxY2tkzKitUhLgG8uQIIZYoLWU9cKCbNWuiQHGh\np2J+GOPiuRVdXVE2b17FVKaXHx75d37RP0Qoll44pwJwqF04NHY+cIeNbWudFFNGyp0FZnNW1q7W\n4/U6mElliMYznPZP8NbxKSbSETr7JxieirPv7AD7zg7w9HMH+PimVj79sTXU1bgAiEaTDA6asduT\nF93Lcr0SVzqcJCtUhLg2knAIIZYoLWXNZp3097tpbJxYMoxxYSPuLLNQKD/L37x+kIGp8ytDastt\nfGhTPTs31GIpFDn4ZoKKCjWbNlnnl6+uW9qDolROAzDmL7DZ5+HD97QA8NWvHuVYwsyk8ihxY5yR\nqTjffrODf+s4yRd/9U4ef2QTDoeJ+noVDofpontZrlciEIjS0REBFLS1OWRehhA3mCQcQtyirrYM\nd6ks+XxjPHFRY1xqxD1Vc/ztC4f41k+OMjk937NgN+n4lR1reHTnGjb7KlAoFAvv23n3/H/37w8t\nSVjO18HoJRCYn3BaKhAWCERJJHLkYgq2NW3iox91MpwK8OfP/YJwOskffucNXukc4MOr1tLQEKOt\nzclil+uV8PsLDAzYKBYNuFyphc+7mphJyXMh3pkkHELcopYbTljJXirzScfFDWcgEOXUmTQvdXXy\n/Lf85+ZbwLo6N5/52Fa21JYxOgwVRvWSZGOxC6uBlhKYdLqI3V6xpECY319g3bq1WCz96PUpDh5M\n8fDDTp79/Qf55vN9vD50hp91DnKif5ovfnTHFTX2Pp+acDgCxPD5HJeN2TuRCaVCvDPVl7/85S/f\nzAtIpXI38+PfN0wmncRqBSRO52k0ORKJGD6fGovFsPB6Z2eCVKqCfD5Febn2sucIBKJ0dibQaHKY\nzXqe+VEXX/3xfjr8AfKzc6yr8vC3v/8Qv/uh9eSjBgb640AdiUSMujrTkveXrqGvL4/F0kCxmKOu\nzoTFYqCuzoTVqmByMkyhkCYSmcVqVZDPx+nrm+G++0yMj+vIZu8gFBrH5bSx2ruZnWtc9E2FGJ6K\ncrBvkAfuaKDcbrro2hfff4nFYqClxYnVWlr+msPpVF4Us5X8Ti0X69uJPHsrd6vHymTSXfJ1STje\nJ271X9DrReJ0Xqkhv7ABLDWOGzaYUasvX2y4lJwc6OzmC//8Gv/eeZp0PkdDmYVPNG/hDz65lnqP\nnRdemCKXs6PRJCkUwqTTuYWGPJWqWEhAFn++z6cmHs8sJAVer4NQaBa/30N3d5S33w4RiUTJ55tQ\nKhO0tqoJhcbZvl1LTY2RoaFRVEUFn/74OkYiCU6PTPHCwbNE+9RYDVkmJtQXffbl7jGRiNHa6roo\nZiv5nVou1rcTefZW7laP1XIJh5Q2F+I2sbjc+I4dbqqrLz38UDouEIjiqijwvY4X+MrLb9AdmMSq\n1/HXn36Qr//qx/jglntIJEyLCn31Mb8eBQyGpiUlypdbSnqp3V0bGydIJkNoNK309EA8HqWjY4yh\noTm2b9eSSMwnDy6XFYOhicCogm8/8RHu3VDHdCrLd44d4ZWfx86VSz9NJjN+2Z1qr/eOskKIS5Mn\nTIjbxErnGZSO+84LB/mn9iPEMllUSgX3NTdwT0015rgJW12ahHJxIx3HYLCTSpVx6lQ3ZvMMW7eq\n6egwEY0mCYdNtLXNH1kqEgbxi1aRlOaOmEwJXnnlJDt3QiSSYnKyhv5+N4HAGHa7k66uINu3awmH\ne0inFUSm1Hznv36MD3/xe3SPRdg30cX2+IP4fD4OHz6EXq9YUrFUJncK8e6ThEOI28Q7FboqTSY1\nmhL870M/4wdvd1IE1lV5+H+f2EVosMjhwwr6++dXdSwuLV7ai+TZZweYmKhEoynS25sjm7UxOGim\nvl6FyzXfi7F4YuiFZdFLy2R7e/XU1xuIRNJs365laGgWmMBkyvD886cwmcrxeucABf39HmCC3bsd\nPPelX+HeP3iWQ/0T/OTYYeZGK7BYZtDpDEDskknXha/JihMhbgwZUhHiNuH1OhZ2Xr0Uv79AKGrn\nSz88xHPnko1Ht23g2T94mJaaMnw+NXb7GOFwF0eO9LJnj39hmKLUSHs8CmpqVORyIxgMGRSKs5jN\nA9jtY/h8ajKZcY4d66OpKbNwHaUGf+/eYb72tU6efvogb72V59ChHhyOWhIJ08Jy1yNH5jCZqsnl\n5oAiUEShSJ/7M1Q6LXzxk9sB+G57D2rzGjQaG9u2pWhrcywMs5jN5wuEXTikcj120BVCXEx6OIS4\nRVzrN/O8boanX3yeiVgCu0nHn/36Dry6SqZCcQYH5vD51Kxa5SUaNXLq1BinTsU4cybO44/PN9In\nTmQ4dSrOtm3gcFRhMDQxMXEUrdbLxESAiYkZXnmlgFa7iSNHutDrQ0vKpff1pTlxooGhoTB2u447\n7lBQXV0adinQ3+8hmwWtdhSnU0ddnfPcfUeoqzu/2maVtYa2htV0DJylM/IT/uI3t7J58/lltj6f\nj0QiuHD8hbU6LtUTJL0eQlw76eEQ4hax0m/mpUmho6PnJ1H+4Oen+M2/+QkTsQQbGsr52dd/A6+u\nklSqgvb23MJ5fT41DoefeNxPf/8MJ0+a2bt3hL6+cV5//Qz5fB3d3QVMpgR+/2n0+gSRyDThsJnv\nfz+C02kml+vE4+Fcsa/IwnlXrTJQKJxErY5hNI6xcWPZwrBNOBzH4fCzcWORbds83HHHNhIJ03zF\nUt8ahoZmFyaG+nxqfuMDq1EqFJyYHKWs+nwxsEv1cFzoUj1Bl4vt4km2QojlybLY94lbfRnV9XI7\nxykSCXHixAQ+3xyVlfZlj1tch8NqV/KFf3yNr//wLQpzRX79g6380x98lDKrcdHy1TmKxRyZzDg/\n/WkEvz9KJlNGNpvC5bKSyYQZHq7AaJwjmQyzceMqwuE8Pt8azOYcDQ0KUqkEBoOJkZEpGhoK6PU6\n4vFxTp6cYWjIwuxsgnvvrWB0dAattoHa2mkqKoy8+aaf55+folBwoFRCRUWat96KcurUGbTaOWpr\nlRSLOdLpHOGwlc7OSVpaTGzfWkV/cJqTgyHSuTy7tviApTVANJrcZet0lJhMOnK5+LJ1NhYvq73c\n8ttb3e387F2pWz1WUofjfe5W/wW9Xm7nOF1YUKvkwgJYpUSizqfh1776PC8f9qPTqPir332AP3z0\nbtSq+Y7P0rETE/PDHm+9laOvr4qRETU2m4KKimmcziKpVBat1kVlZZZHHtETDudxOqcIBFK0tqq5\n//5GGht1DA6mUCisnDypoqtLy+TkFDpdHclkmoYGBdu2eYnHcxgMdvT6LIVCNb/4BRSLHoLBIaqr\ny/jxj0Mkk6uZmVGj1doYHIxiMhWpq1Nx9GiQWKyc2dkZWlocNFe7+Ke9nfgDUT794TvQqFRLaoB0\ndEQ5dcpBKjVFS8vywyQmkw6VSrlsnQ0p+jXvdn72rtStHqvlEg6ZwyHE+8Q7zSNYbhXKhaswvF4H\n8fQEH/+TH+GfjFBuNfCHj9zPWoee/ftDS86/+L3bt2uZmRmguTlHfb2baLSSY8fs6HTTNDSk2b27\nDr+/gN2uY9++KUwmA7FYGBgkkTDR1JRhbCzJ5OQYUI/ZbMXlitDSYl6YFNrW5sTlKmA2OxkamiEe\nHyEYNLJzp4LBwQhWq510+jh33aVnYGCYoSEX4XCWwcEohw6l8HhmARsATVVONvk8dPonePXoAB+9\nq/mC+RqRhQmn8xu5ze+M29bmvOK9Z6ScuRDvTBIOId4n3qmOxnIN34V7l/T0T/AfvvoTRmPT2HQW\nHlvbRoV5He3tp/H5fEvOn8mM8+qrIzz0kJodOzaweXM9MJ/8PPPMOJGIirq6GLt31ywsce3qCmIy\nVTA8PIPZbKW9PXfuvFG0Wi/l5VbU6gCbNpmw2cxEown27k0CIzgcpRUpNhKJAhs2NHLnnU2k070U\ni9OUl9vYtKmOujoVBw70UywaSCYj9PQo0esfIBZ7i7a2hoV7/8QHWuj0T/APL5xgS+3SeRml5Mbn\ncy5MSlUo0rhchYX7DwSinDyZxOnMyWRRIa7RNSUcx48f56/+6q/47ne/y/DwME8++SRKpZKmpiae\nfvrp63WNQgjeeav15Xo/5idWzq/M6O2f5D/8xUvnkg09X9h1F3WeIr29p2lqyqDXBzGbk+zZEweK\ndHTMkEy28vOfnyKZ9FPayn1+75FqDAYlFkt+yed5vWpyuSBarRKbbQaDwUI63cv27Vr27h2lslKJ\ny+Wmvt5Kf7+H7m4/iYQNiyXL2rU2XK4CfX0BOjuNNDRMsmWLhUAgQS5Xi14/SltbHc8+O4Ra3QKc\n4Xd+p5GzZwP8+Mdv8ZGPKJbc/+67m3n62Tc5OhjgTE+WUs2PUpxKicXExCC53BgejxKfr2bh/fP3\nWUs02i+9GEJco6tepfLtb3+bp556inx+/h+br371qzzxxBP8y7/8C3Nzc7z66qvX7SKFEJevo3Hh\nKorFKydKdSYqq4v89t+8wkgsik1n4KW/+FV++zfWo9dX4vOtQa+vZMcON4mEiePHrbz+uoZUag6N\nJgeo6exU8cYbJjo65s+5adMs5eUBslnrua3s568jnfYyNqantXUTmYyZ4eFKzpyJ4/HY2LJFyfh4\nnHhcyenT0yiVXRQK4zidM5SXh87V+YjT15dmZsZLT48SszlJMFigWEzQ0mLB63Xg8SixWBKsXq3D\n47FRVdXEf/yPO6mqaloSl0qnhfpyG/nZAqHUOHv2hBgZ0V202iSRMLF58xZWrapcEl+fT43ZPJ+E\nyUoUIa7NVSccdXV1PPPMMwv/f+rUKbZu3QrAzp07aW9vv/arE+I2dDXLLH0+Nel0L+FwbKG3o5SA\neL0ONm+z8cR33uD0+CQus5Hv/eGHufuOeuDipaI+nxqdbhSXq8DGjXbuvz/BqlVGcrkRnM480Wj8\n3HwHWLXKiMXiARQEAlHC4RgnT3ZSLBoZGjqLXp/g9OkRhobyPPvsEM8/H0evr+fo0SlsNi9arZOH\nHtpObe0cbjccOpSkuztNsZjHaDzBunXVtLfnsNvrmJwcIRpNEghE2bWrhs2bYd269Qs9Ful0D+Fw\n/KK4tbVUAXB8IInDUcv09Hwv0aWSsgt7j7xeB/ffX0EiYZJiYEJco6tOOB588EFUKtXC/xeLxYU/\nm0wm4vH4tV2ZELepd6qncamExOt14HJZMBiaL9o0zT8wycP/9TnePj1GmUXPf//lh6gusyyca75B\nr1jYFM3rdfD443Vs3FjE4bAACjKZKsbHc+Ryo8zMpHnjDRMHD+YJBgs0Nk5QV6fk618/yne+E2Rw\ncJB8XoXHA/X1laxZU0Mmk8TvdwImMplBtmzxMDAQoKkpQ01NlpYWKwMD5aRSGzl6dBKdbg3btpkw\nGqOcPTvCm2++TSxmZGiobCEuuVyEAweOk8mMn7t/68KmcaUYHTs2iE1pBqBrKIDBEGD37vleoguT\nsstVYZUN3oS4dtft6VEqz+cuyWQSq9W6ove53ZbrdQm3PInVyrzf49TWVqC3N0ZTk2PJvYyORunt\nzTI5mcZoXE0kEmTjRssl31dd7WDjRsjkCmz/9P9Hz2QYi87Ib2/dQU3ZRs6c6SESCTI5maG2djUD\nA91UVNjJZgtUV89/biSiJZGoYGTkEPv2nWZmxo7Z7KWiIkx19RzB4DRVVS3E4zGiUTXj41VEoy5i\nsQSbNpVTKPQzMpJl7do8lZUmTpxIUVmZ4Y/+aDX796dwOFZTWZnl/vsrGB2NMjl5nDNnZli3TkUs\nNkNXVwaXS0E2u5GpqbPU1RlxuWK0ta2itzdLMNjM7OwsL7/sR6Uap7FRQyIxf/8vvtjH4cNGLJYw\nzevWAx2Mx3LE42pcLgtut2XZOF8Y76amAhs31rJx47vx039/e78/e++m2zFW1y3hWLt2LYcOHWLb\ntm3s27ePu+66a0XvC4WkJ2Ql3G6LxGoFboU46XRqWlvnH83F99LRESWVqiCdnkCj6cfpVC/8/eJJ\nozqdiVAoTrFY5LPf3EvnYACbQc8ffOgu7lxvIpHoJxZLYTQ2E4sFcbkmqanRks9X89Of9uByzQ8x\nOJ0Qjfaj0Wiprq4jGu0nkejB7dYQiWjZtm2Offu6SKdNFItR6utTpNN+TCYjR468RSIxR2Pj3ahU\nIbxeJ42NHhobtdTXl6HVRtm7t5s9e/r5H//DwWOPOUilsoyNWVGrxwkEVOj1qwgEOkkkItTUVNHU\nlONXfqUWnU6N0xln9epR3nxzkmTSzf79GjSaIjt2mAgEorzxRpRksgy1OoNNPw1AODZNf3+Bn/50\niN271YTDcaLRAuFwDp3u4n8KS/GGGDrd+/t36t1wKzx775ZbPVbLJVPXLeH4whe+wJe+9CXy+Tw+\nn4+HH374ep1aCMH5VSrr1zsu6vrv6IgwMGCjr2+YVavmE4+Xjg/ywze7MWjVPP+nj9Jaf35310BA\nTSQSpK3NycTEDAcPJtDr95HJmKmtdQNZzOYkBw/OoNcnWLvWgtNpob5+LZHIMA5HLXv3jtPXlyEU\nMlBWZsRoNPLLv2yis9PD0FAIs1lHJNJNU5OFZNJEY+MEbW3z1+31OggGBzh2zMrc3FoMhh7S6TyT\nk/Ukk2M88kgtL7/cSTJpxOlU0Nyc4/HHGxbu2+t1sGsXBIMFwmEFuVwv4XA5gcD8vivr1zfT3z/E\nJz9ZxfoNNSi+A6lCDqOpDJhaErNwOMLu3Y5FsZlP3ubntARpanLIXipCXAfXlHBUVVXxgx/8AID6\n+nq++93vXpeLEkIstVyDV3o9Gk1SLFYwMQFebwUv7jvOn+75OQDf+M8PLUk2YL7B3rjRwvHjw/yv\n/zVCKlWH0RijtdXDsWN9uN1qXnyxQCi0GpstwCc/aWfXLjUdHWNkMhnefns/p06lmJqapVjUEAxO\notW6aG4ucv/9Zk6cmKC/v0Bjo5bDhwtYrVoUijQTEzMLjXkikcNkSjE3N0RzswLQkkoNsm6dig0b\n9AQCJo4fd1EoKFmzRsPExAxf//pRYjF47DE3en0lNpuXkRE/hUKO48dtQJS6uvly54895mTz5noC\ngSg2g4HpVJqmdRHaNpbip6BYnN+2frHS3A4IsmOHG7fbstDbsVwNFCHEO5MZUEK8DyxX9Kv0usMR\nY9WqFGazjYFxP//PK/uYnSvymY9t5RMfaFk4fnFFzUceqaOjI0os5iQcjrB69RwzMwEaG+t4/vke\nwIhCcZiGBjXhsAqfzwEoiEZ9xGIR7PZmEolx1Oo8anU9RmMZmUyC3bsbcbks+P3d/Nu/OaivD6FW\nT9PQ0MTY2CgaTSWDg3243Y3s2FFg2zbnQs/Hli0FMpk5urqirF+vpFicQqHQUFc3v4nc0aNeslkr\nzz8/yFNP1XDgwAj5fDkTE4MEg2Fqa7O4XJVLdoT1+wuYtCamU2ma15kWEra2NgcuV+rcfZ13qXon\nl6uBIoRYGXl6hHgfWK7BOz/MMl+OO19w8sf/+1+JJjLct7GOpx67Z8nxiytqHjwYoqNjjHS6iEqV\nZvXqrVRC7Ya1AAAgAElEQVRWzhAIhCkWNUSjs+zaZcfhMHP8uJbTpweoqFCjUFhpaCgAAYzGGSor\nK9Dp0ng8MXbtcnLs2CAvvxzk2LExwMPp00EcDgvj4300NEQJh+cnpMbjQ7S06IEip06NcOTIHB6P\nkr6+FNGoD4cjTHW1B4ejlkQizvbtnJu86qK5eb7A15o1UXI5K4XCMPX1JiBLOBwD4qxf71iIkUp1\n8Yqf5SqzXup1KV8uxLWThEOI97grmT/wtR++xZHecbwuM3//2V9CtWj1WKlOhsMxjsNhAgzkcjbC\n4QIul5KZmVHuussGqBkby1JR0YzDMT+xbXR0DlBhMMSZmurD78+RTHpJp/VEIlOsW2diYmKab34z\nwPi4BrW6BZMpw8zMGfR6mJ724XAUUKvnzk3AVHL6dJhw2MqRI+PMzU2g1W6goiKLx1NArQ4DcxSL\nBo4e7UahyBCJlLF7t42BATU2m55AIHqul6KA2exjaCjH6dMF9Poqqqvnr7u0N4zJqILwxXGdL1g2\nXz1V5mYIcWNJwiHEe9ylhlNKSUhf3zjRqI9weAJbwwx/+8IhlAoFf//ZX8J5wc6lpUmSDQ2c24sk\nRj4forq6DoNBRUuLmUTChMFQwbZt4HKB2azkzTfHGB6Okc1qOHs2T3l5LVNTSYaH/RiNVjQaE+3t\neoaGrMzNxaiosDA724nBYAIyTE+r0ekGqKtT8Ou/XsPoaJHvfW+S2dkGBgaGsVjMOJ0WjEbIZM4S\nCLjYuHGUe++tor09jVZr4e/+LkRLSx0ORwKvt5Zo1IDfnzpXOwPATSIRoq6ujJMnjzA0ZKKvT4nX\nuxm/P7gQg9BkjP1DpZ6hAgMDNopFAy5XSnowhLjBJOEQ4j3uUsMppSRkYmIMnS7NTDrNH//dfgD+\n6P/azl1rqheOPXZskPb2HFNTI4TDrdjtSQIBPcXiNiorM2zalKK3N0Y0qqCuTgUEF4Yj9uzJsX//\nLOPjFaRSI1RUNBOL9TA7q8TtbkahmKVYPMPEhJJCIYNaPYpW66C+XkMqVUt39yyFQgGLpcDatUb+\n9V/HCAaLrFsX5u23c9hsetTqLLW1cVpbU3R2mkkk7qGnp4P/8l/q8XiifPaznSgUDQwOnuRXf7WS\nZHIGiOHzOZbMSamrU5HNBjh7NofX20CxGMTh6CEQSBKdyQLzPTVlhoqFeIbDkSXnkpUoQtw4knAI\n8R53qfkDpSRk1y47iUSR//l6N1MzKXa01vDonU1Ltplvb88xM9NKMBiiqUmFw2GiujrJgQPt3HOP\nEr2+kulpL93dMbq7Izz2GHi99XzjG2+xd6+a/v7TZLMVmM15HI5+fD4bqdQshw/34PEUMBo96PU1\naLUDuN1lrFmzAZNplnw+yORklkKhnExmlgMHokxO1pHLuUmnT/Cbv7mF1147Tnm5FVBRV7eakyf7\nGR19A4slRSAQxet1cM89Otrbc2zfbuHBBzcQCETZu3eEZ5+N4fEoiUZ953Z5LZJO56iuvpOpqWP8\np/9URSJhor+/glj6AADrW4zMhILnyrmblmxFv39/aNmVKJKMCHHtJOEQ4n3ofBLi5tWj/bx5agSz\nQcuXHr2bf38xjMNRC8TxeqGpKcMrr7zN+vVpZmbGiEYBLDzyyEby+X4ymXFOnJhibCxOS8tO9u4d\nIpEI0d4eJRhcQyZTj83WiMt1kpoaO2azjfb2XqCeSCREoTBJVZUBs3mOtrYtxGJBWlosRKMGslkf\nmcxZ9PoNWCyQzZ6mUAhyxx1ZWlvnuOuuRtrbc5w9m+If/3E/mUyITMZLT4+Db3+7E63WiV5f5OGH\n53A4PAsN/1tv6YlG7VRVdVFZWcTjAZ+vFrNZC4zyqU9VLSyJDYWCZH4x38PRsqoCbYuK/fu5KLm4\n3EqU5VYJCSFWThIOId7D3umbda4wy1P//HMAfvehjbz1RpJi0UA0Osy9987X3tDrK3nggQr8/tNo\ntVVMT6dxOGYwm4M4nWqefbZIPL4Nna6dQOAY6XSKbLaK2dkkmYwfkymEVhuhWNQyM2NmYGCCigoN\n/f0qikUF4bCJXK6H8nILJ0/6Wb8+zaFDGfz+KDU1rTzwgINgcJLmZgXr1rXQ26tn+3YtHo+ajo5Z\nDIb5YyORViKRMMXiBKmUjddfD7NmzQ602gEqKuJMT9sYHDxNOq0nnR4iGrUyNRWjUKgB0gBs3lzP\n5s0XxG6NnrliEZfFgFY9v//T4uRicYx37HBfGOKLjhdCXB15eoS4yS6XVJS+WXd09J5bjZEkkTAt\n/Hf/0Gn6x6dZ5XXQVtXEhNLO0NBZWlrMC+coNZbbt2sZGpqgNN+hxOMBg6GPmZkYNtsm0ukEIyPj\nzMxosNvrUKmUlJUZ0etXMTrajtXqIZ2OsmOHkmPHkkQiLaTTehQKC2vXtvLaax1MTjaSy/XjdI5S\nVaVjeLiKwcFBQqEUuZyJsbEgZnOEbNZJMukgl+snmTyLWm0gn49it9eg14+gUByioWGOWCxNZ6cW\nCOPx3InfP0wmo8ZorKO3N4Tb3UBHRxSXq7CQFOzZE8JuryATnK8sWuE8H5P5OM/HPRyOn9v07eLe\ni2PHBjlwIIHZPMeuXTUynCLENZCEQ4ib7HLd9aVkIZ0ukkpV0NV1Gp/PR1fXaVyVHv7+p50AfOU3\n76OlXIdGFcdgMJ/bNXb+fOcbVxNtbecb49paM9HoJJs22ejujqBUupiZieF2B9BoHFRVeYlGI9TX\nz9HQMEcmE6Ciwozf72V2dppUahKvN0skcoJi0YHF0gsUKC8P0dPzOqmUkpmZOQYG3AwOxikWjWi1\nc0QivWQyJvT6Ahs2BNBoXLS0tJFOH0ahWMPc3DAwyYYNZXzoQ/cQCBzjwIEcKpWTysowU1O9qFQN\n5HLT6HRD3HOPmQ0bMkBxIY4ADkct0egwSu/8Ellf5aWTOYgtuxNse3uOgYEGZmcTrFpVkOEUIa6B\n6stf/vKXb+YFpFK5m/nx7xsmk05itQLvxzhpNDkSiRg+nxrLBUtZLRYDdXUmrFbFuWPmKBZz+Hxz\n/O3LP2cwPM2962r4wq99YMmxQ0OjBAIROjuDDA2lOHlyCr/fQyo1RSajIZezEwj0otEUGRyMkU67\nOXOmk0TCjMuVZdOmzUSjIZqbk9jtbj70oXIefbQFj0dBX98YWq2dSCTL2JgLjSZLZeVGTKYcu3Zt\nJh7PMjioRam8G7NZQXW1gkKhwNRUL8mkiunpfqJRDdlsBQZDko9+VMfAQBifb46GhjwaTYHW1g9g\ns2Xw+VScPh0DWpiZOc7OnWXs3KkjEJhCpTKybZuXhx/2sWtX9aIYqXE6lczOZrj7bguvd4/O1yZR\nl7GqTEtlpX0h7kNDo4CC9etNl+y9UKvjxOOT1NYm2LbNcdHPR5z3fnz2bpZbPVYmk+6Sr0sPhxA3\n2XJVLBdvIpZImM4NudQDcLgnT/vACGqligcbVy9ZleL1zhfD6u+vYGgoQF2dh6GhAYJBJUpljLY2\nB+FwiP7+HMPDNqamwhgMaiYnXZhMTQwPHyUW20+h4GR4OMWGDU28/PJZOjtj9PWlqKhI4XAomJxM\nYLNtIZHYw9TUcQoFDd///kEqK5W0tpoJh4/i9UYJBMopK0uRyzUzNlbD7GwUg8FJoTCKy6VlYECP\nRmPBZLLzwQ+C2Zxk795ecrkIBw8qqa9PEY+fwmrVMjnpIJ1O89nP+hgamgOKCz0TF8ax9Ofu4fkh\nFYuylfb27MIcj1KcLrcyJZEw8bu/W33J3WSFEFdG+c6HCCFuhlKXf3t7jpERHXv2hAgEooyORvj9\n//kzAO6qbkBTKA0lnC/f7fOpcTj8FIv9hMNdwCxg4cCBDBMTM4CC7u45XnttjHBYw+TkBGvX1pLL\nHUGrdTE2ZmF83EcsZqen5xAGQ45jx+wcP17D668XmZx0cc89Ru64w09lpREop78/SSRSYGTETTod\nwWZLcPCgkr6+WmIxLU1NWaxWPxbLNIXCGLW105hMOgyGHA6HCZ0ugM+nZvPmejweBQcOWHn7bStH\njuixWMxUVm6lq2sYu72CoaE5otEEZ87Ez93PfIKwf/98jEp/Pnp0gOO9EwA0lkdpasosHFOK03LD\nKaX49/Zmb+SPWYjbhqTtQrxHmc1JurpO09SUobc3iELhYs+eEP3pMfyTYQwqLR/beD92S+iiRtPr\ndbBqVYHpaS/FooFVq/y89tooanUZf/7nA9jtM4yNmVAq58jnDTQ2KvB4CjQ36zl61IzdHiaX6ySd\nThGJlHHq1DQwSSJRJB7PceKEFYUiyy/9kh2FwsLZs2qMxhpSqR6UykpmZlRMTnrI5RxkMt1s3qzh\n/vsrGR+fYnCwDpPJx+TkCYpFHZBHrw/h8cxf+7Fjgzz33Cjj4waUymHWrKmhvz9MIBBk7do8RuMU\nUOTYMTuFgov29jE8nih79oQWlgPD/LLXNzsOE8tkMeu0fPHzGzlwYGpJj0apV6SUoCyeuFuaP9PU\nJBNFhbgeJOEQ4iZYSSGp+WEUH3p9kN271ezZE8JsdfN8+ysA3Omq59SJ3nPbsLsvOrfZnKShYRaI\n0dZWQ3X1CF//+gSzs2tJpQbx+aopFvvwerVUVzcxOtqL319ApUpQXq4gl9PR369kYkLL1JQehyNJ\nJqMkl4ujVJrp6UlQVWXB4SjnjjuKJJMZikUjc3MeAoE+VKoYKpUalytKobCRb30rjFptRq/PEIud\noKrKw/i4gUwmgVpdQyYzh99foKsrB7SiUEzR1KSlUNCQStWRTtvp7x+iWJzmoYfUbN4MExMRtm+3\n4/cXsNsrliwH3rv3GEf6BwG4q6UShUKx7PLWS03cLSUjbreFUCh+XX/+QtyOJOEQ4iZ4p0JS8xut\nxYHYwk6wu3fDd17qJZ7NYVFaUYTUmNeWMzQ0Py9h8f4q09NeGhpmaWtz4vcX2LfvND/6UZJ8Pkah\ncByvN8bmzS5aWtyMjhY5dOgsw8NxAgE3Hk8Si0VJOu0llzuB2WxndjbHzIwetTqDRmOkUOgiGp3l\n5ZcTPP64iUcfbSQaVfPaaxG6u7uBFHa7HqXSyuxsgfb2OJWVTeRyp6mpqUarBYtlAoVivr7H6GiY\nUGiA6WkHTmeBtWvVrFkzx6ZNZVRXK/je94I4HHFmZoIMD6/l298+w+bNDezaZV8o8AVZ7rvPjdc7\nX6a8u7vAaGp+1PiRO5su+/O4kjobUnVUiKsjCYcQN8GlGrjFDZnfX8BgaMJoDC40al6vg9d7zgJQ\np16FQZ3HYnEA83MUSpuzTU0lcbkMQGwhsfnxj3uYmWklEplg40YNra0Ftm3bxhtvvEl/v4JoVEk4\nnMJqNWKzTWC3Gzl7tp9YTEdlpZ/GRiuZjJKJiTxW62qiUSVDQyZyOS/PPXeKX/7lIKGQknC4mkwm\ngU63E602TCZTIJ3WUiyOMzg4Q1tbDqsVlEo9xaISna6Kw4cPEgyWA0mi0Vaamws8/rh7oQjXnj39\nuN0O4AwOh5FQKE82W8bQ0Dra2/vZvPniCaN+f4GaWheB4fn9Ze7fNJ+UlGpzhMNTCzU7ShNtV7rk\nVaqOCnF1JOEQ4ia4VAO3uCG7VEJyvH+CE4Mh9GotGzwm7tqWJxIZo65Oe+4Ixbn5GiZWrUrh8znO\nnTfIRz6i4Ec/OklVlZaWFgu7djmJRs+SSORIpy1EowUMBjtWay96vQGPp5Jc7hSFwgYikR7uusvC\n0FAFSuUI1dUTVFXNkUwGGR8fJRAw8uyzcRobvRiNERyOEfL5QfJ5A1ABjJLJaHA6a+nrg4qKGTye\nSvr7R1Ao8oTDVlSqNnK5oxgM4xiNOcJhHT/72Tg//3kKv3+K2dkqcjkbzc2V+HwR9HoFmcwptm+3\nLcRnccLm86npHouTL87SVOWkuszK/v2hhdocRqOKkREdXV1Bdu/mkj0VpfO1tRWWrFKRqqNCXB15\nYoR4j/D51HR09JBOK/D5HEvKbAcCUZ7+h/lv69uqK7jvbte59zTR2XmMvXs70etTKJVB+vqUQACz\n2bawnHbHjrvZufN8gwxw8mQcr7cKsznE7GwEv9/N8PAELS3lFAoh7rxzjrffHsdsdpFMJkkkMszN\nGZiZybBzZznp9BixWJ5kcn7YxGKJ0NxsIRK5l+PHh8nlLCQScXQ6K1rtWvL5k1gsLhIJ6O09QT6f\nQq3W4XbPkUi0s2NHloce2gSAwdDMK68cor+/mVzOjE43QXOzG612lJaWckwmBb29eqamYguTPRcn\nbDt2uInSDcD9G+sW4gvxhTkee/YEcThq8fvjl+ypOL9KJUZr69IJudKzIcSVk4RDiPeIS9WFKH3L\nHhwL8Xb/KAAVczYMhibS6R6MxiDBYIGRkdWo1WF0uikSiQ10dg6QTudwOCro6hpe+BZfOueePSEs\nliri8V7WrLFjMGQ4cSIE1NLdPcC6dasoK3OyY0eWVMpCfX2e6elOTp+2MDtb5ODBMoaHi5hMVjSa\nMDU1Wh57rIp162p45pmzFItxQqFRqqstZDLTWCx5GhsNbN5s4yc/mcJiuYNAwI9WO00qNcsdd1RQ\nW6uhrk5FZ2eMM2eO0NgYYWRknNnZAhs2mNmwwQSYMRiaeOWVQzidrbzyytvn9om5uFdo72E/ABuq\nyxbiuzhR2L0b/P74sj0VF65SkbkbQlwbSTiEuIGutJG6sNEsfct+9WQ7c8zhohq7zoLRGFyYTFoq\nlKXXpwAIhQ6zapWBTZtstLfP163o6AgtzFkoreiYnQ3T0mJhbExDT4+CNWtmGB6uRKWykU5rSCar\nKS8PUyw6iMWmOHrURrH4QQYG9qLVTpDL6ckok+jtUdauvYfR0RTJZBQoYDJZcLvryWZj1NTksFqn\nuPtuC3q9mcbGQQYHj6LTzaFQ1AAmTp9OMz2t5KWXTqDVuti8uRYAi0VBJmMkkVAQCMye2w+ml6qq\nDJlMBw89pCaZXNorFAhE+eY/H+PU8BRalRqnsuySsX6nnooLV6nI3A0hro0kHELcQIsbqdJmYRcm\nHxcmJYsbM59PTU9vgF+MjgGwqdzBffcZl3xGaYfU/ftDpFIV1Nf34nJZ8HjU577FZwmHFZw8qeTF\nFwfYsiXD8LAJn09HXZ2ePXt6mJpqorIyzNatc0QiRYaGhs8tU1XjcpVx4gTMzUWYnv53HI4sKpWW\ncDqH1manMJfjwNEI+/efJh5vYtWqNOvWGejvH2ZiwsjUVBlNTV66uibx+1XMzNioqGgmkRhEpUoR\njU4xO1vG4OAE2aybVEpFofAWBoOScFiJwaDh8GEVdns9Q0NqDh0KMDnp4e670zz44Ab27w8xOmph\nz57hc/db4MXDAQDqDeW4HNencJfM3RDi2siTI8QNtLiRWu4b8jt9cz4yOEYklabCZuH/fnILAwOz\nFw277N07Ql9fEre7l3Raj91uo6srxO7dbnw+NeEw/OIXnRSLG/nBD/YxPe3GaJzG4zESjUZRKNTE\n4yqy2SomJ8NYrTWMjg5QKFQQCJzAaMwzPa3DZGpFqcygUERpXbeW4amjpNMJopFhomMmtNoyOjtP\nYTDUUlOTY3razvR0gcOHe9FqQ4RCWozGYTZu1KHTpTh1Kk0iYcVmW4XH08nQ0DRmcw16fTnJ5Dp0\nuiAazQStrR/gxIn9HD2qJZOJo9OtZWJieCHGXV3DFIsGnnnmLNlCmq5QLwCfaNtOImG8OKhXQeZu\nCHFtJOEQ4gZa2khFL/kNeXFScmFvh99f4KXOQQB+4/41DAzMYjYngaXDLm+9pScYLEOnO80HP9hI\nV1cP5eV17N07wsTEHMPDGgYH4ygUHcTjCubmXIyNKZia8mKxFFi9eoCKCjXd3SMkErOo1Xn0+jTh\ncJ5QyEU2W8ThKJJMHqeszIzLNUMipqbOrSZfWMfhvmmKxTCZzFEslnL8/m2Ew+0YjSpUqihO5xYi\nkWOUlzexenWOz3xmNf/8zwNoNHbc7ix6/QBNTWZ8vjXk88O43VoGB09jNCrZsMFEJDLK9LSSQuFD\nZLN72bZthF277IvirObMmTCDgw76s3Hyyjm2NVXQts4oPRJCvEfIkyjEu2S5b8iLXy8Ni5SGYMaC\nEbrHh1AqQBmyM6qwUF3NQo8JRPH51JSXh+jvHyGRyDM62su2bVZOnZri+9/vIZ1WMjOjQ6erRacL\n4XA4CYfPUFenJZOZxW5XU1trxuOxEAyqCASy+Hxgs9kYHj5DMqlBoVBjMmXweHTk8xsZHDyKxeIm\nMaakomKampoEQ0N2dDoHhcIZdLojlJcnyGTSqFQzzM6+id2uRqXq4L773Pj9BXK5cgyGNBpNgLq6\nZrLZfsrKtDQ0uNm923dRnJqaMrzyymE+9SkPDz64YeH1Us2SlpYeYrFJ3urvhzn4zx/bxo47l670\n2bt3hGCwwMMPO9m8uf76/oCFEJclCYcQ7yEXDsEMTeuYnSvS4HKQiJoZTHXT0uJkz54cdnsFkGXH\nDjef+cxaotEzpFIbsFr72LTJxv/5P2cZGTGRzxfR6VSYTFl0uhTl5WWsWmXCaCxgs6UJBvXE40aG\nh0NEIlEcjkr6+gapq/Oi1aZRqSrR6dI0NCSJxx2MjZ0gl4vhcIQoLy9n06Y4995bzmuv+enuzlFW\ntp65uSlAy9SUEYfjTsJhP1ZrEwpFimRSw7p1ao4c6eXMmQSpVISyMg9r11pxOOY3ljt2bHDRDrnz\nlUP1+ko+//mLJ9+WYrZ+vZOCNc23/z6Jw2hgg9e55Di/v0Bnp4d8Xkt7e2hh11ghxLtDEg4h3kMu\nHIL5xzdPAbCushxF0YRCMcPQ0NxCAav77nMvvO93fqeS9vYxtm+fr79ht/tQqQaYm1NSXW0DxrFY\nrLhcYSorPXi9q3j77dOAiZdeOsLsbB12exSl0guUMz09SbEYp6ZmDo1Gh0pVhsUyhUpVhs22nrm5\nX+B0Wjl8eASr1cvWrR7WrrXzwx+eJJstY3Iyi9HYz/BwLw0Ndnp7j1JVZSUadeH1NqLVlpHPryeV\nOsP4eJjHHqsikTCRSlWwd+9RslkTBw8O8/jj5+e5dHT0LqkQemHMXvjumwCsd7Vy9MgMtTWuhdj6\nfGo2bZrv4di+fWkyIoS48SThEOI9yut1cDY0CcDaSjeFcIDa2gYgRHX1fAGr0rf/jo4o0Wgcr9eC\nx2PD44GtWyeZns5hNhsoFGYoFp1ks/Xk88Ns2aLgwIFRFAoXo6Mxpqdz5PNWYrEMDQ2T6PUZJicV\n1NRsQaEYxevNotVq0Wr1OBwhzpyJ0tJiQKOp5vjxDAaDj8nJM1RVjZHJ5InHM+j1BqzWSioqLExM\nDGOzmQgECnR2BgkEas5tFf82+XwYjaaZl1+O8PDDEA73EI/HOXPmBMmklljsJI8+WgMESaeLl5xg\ne+zYIC+9Ps0rR/pRK5Vs8tQDuYviuWvXfPLi8Vz+n77Fc2ncbst1/bkKcbuShEOIG+xqC0YNT87g\nH49i1GpYW30nOXcfLlcWs1lFIgETEzP4/QXC4Rj9/RUMDc1S9/+3d++BbV71wce/knW15Isky/It\ndmzFiZMmca4uubXJtpJ2BZqOjHaM0ZfxlkspLx1klNKtpO9oMzboy2DtVmCwUtaNUUqAUdoCLW0u\nbu2kcZo0seNLfJVly7rYulkX63n/UKzESZw4qd3Yye/zV6Ln0aOjny0/P53zO+dUOLDZImzaZCcv\nz0x+fiGBQCtOZx7t7UE0mtcxmRZz4MAoa9cup77+AB0d7SSTCZJJDzpdLsGgEb3eg1a7jJMn32bT\nJjuJhB6Vyk1OziIsFhU33FBGX99bRKPNOBwdHD7cgqLYOHo0QFZWEUZjMaWlSQoLuzl2zAtoSCS0\nFBYuIRLR0dDgx+UyUFOzArPZRTCYTyym0NQ0QmNjCJXKQTB4nNHRSrq7hwiFTKfW2EgPn5jNYfbs\nIRPT+vo4e0+GUICbV1Vy4zoDTqf5nJheaEbQ2XvZjJ9XW/vOfv5CiDT1dF5MURS+8pWvcOedd/LR\nj36Unp6e6by8EHPS6ZtX8pKe98rhLgDW1ZSQYx6krs7Kpk12QiETPT16nnnGR0+PHlChVh/G7W6h\nu7sBszmMy+WnsXGYQEBHfr6Nvr4oBQVrsFgsFBdXUFSkYf78GHp9iJMnc1EUB0ZjjKysBAZDGZGI\nhlAogdmcRV9fjPz81QwMxDl+fIzBwRB+v5HeXhvB4GLC4ULicSuBwHJCIQcwxMqVLaxaNczoqBqN\nZgXRaB1mcx1G4yAWSzsvveSmp8fP8eMNeDydaDQHSST6aGsL4/VaOXp0AKNxjPLyFBUVY6cSDA+Q\nTjLq6+P09OgzMV1wXZK3hg4CcMemxZPG1OnUkJ19/rU0zvw5Xeg8IcTlmdZP029/+1vi8Tj/9V//\nxeHDh9m1axdPPPHEdL6EEHPO5S4Y9dqRdMJx8/XVE/ZVSa874aaqqoJAwM3mzXaOHw8wNracUGiE\nUMhEKJTkuuvm4/MdQ6NRU1UVY3Q0TFWVgfnzBwATZnOUxsY4UIzBMEJZ2QhgYmCgH73eiKKMoFKB\nxTJCOPwqPl8Ki2WY7GwtKlUzIyN9JBIaioqSGAwaksmTqNUDrFxZh0rlQ6stwmDwYLUep7RUR15e\nHsuXl/Hmm26SyXkcOnSAVKoYtdqKoviw2ZahVh+mpaUfu72AykorlZUptm2rObUTrgqv14fNlntO\nDctzhztJpMZ43/XVZCctp+o9Tlyw3uNCP6f0+ekeD5vNP2HzNiHE5ZnWT9HBgwfZtGkTALW1tRw9\nenQ6Ly/ErDbZ0MmZN7nxc8zm8IRZGOfT1OYGoNI6cWigpMSSWVHT6UzXcRQV9WAwxBkbG8RsNuNw\n5OH1xjEa9SjKQrTaLt7/fiOjo2b+4z/aGRpSoygBkslcbLYEa9dmUVtbxdNPj5BMDpNKBTCZVIyN\nOdBoVHi9CQoLq/F6O1i3zsjJk060Wj3hcJiFC7NZs8bEiRMK5eVWDhzoR6PR0dzcQXl5BYsWLWDt\nWrkzPJwAACAASURBVDsmU5if/cxHPK5mZMSDVqvC6x0hK2sAszlKPG7E4TDzJ3+ylO7uk9TUZFNX\nZ83c/BXFCIxM2IStpMTC6819PLe3GYM2i51/cQOaZIr2djfRqOqSliI/3xb359u8TQhxeab1UxQK\nhcjJOV1gpdFoSKVSqNXTOnIjxKw0lb02xs85cuQ4Tqdz0nND0Tg9Q0Gy1GrGgqYJx86X2GzdOo9o\n1EN+/hpCoRgrV6Y3goM8WlpOUlUVZdMmO0884WdoaAHd3WqMxh40GiMLF4a4664KHI48XnqpEbc7\nl1DIQywWIZUaxOWKUlo6n56eVoqKinn99W5ycg4DfnJyNhCL9fKhDy2ivT1JW5sLgyFIb6+KJUsK\nyc4OEYvl88orKWCAkydNDA/HMZtdFBaCTpcPFKPTDVJRsYjR0RMEAm62bs3D4cjLrDVSUaHG5eqj\nokI3ofdhLOXly99/GYDP3LaW8sL0dvXpFVjf2VLkZ2/eJoR4Z6Y14TCbzYTD4cz/p5JsSAX41Ems\npuZKxamuLklr6wjV1ZZJ2zB+zqJFVkKh85/b2+vnV79P1z/Nt1tYv8424ZyjR8MMD1v56U+Pcd11\no2zYYKe2thybLYfW1hhms4qjR8Po9UGGh31YrZWo1X7s9hzWrNHw+uvNJJMxkskIdvtCliy5Dq3W\niM2mp6JCz9GjAZJJK4pSiqJoyc0dJBbLI5UapK/PTTxeiN1ezPXXn+TYsU5aWnz84Ac6yssXceRI\nEp0uB4tFjd/vZ8mSOPv2RSguLiMcNmIwVODzNVNVtZF589pwOJL4/WpKS8tIJgcxm3UkEqX4/cNo\ntTq02nJ8PjdarZ1Vq4rQat3Y7TkcPRpGqy3nR797maOdHuYV5vJ/P/4H+IZCtLbGqK7WU1tb/o4K\nPu32HCkYvUTyN2rqrsVYTWvCsWrVKl555RVuvvlmmpqaWLhw4UWf4/EEp7MJV63xHSvFhV3JOOn1\nmkzX+3gbzu6N0Os1WK1B2tsVnM44er3pnPY2NPh580R659dViwrR6zUTzrFa47z2Wgs+n57XXlP4\n3e8aMZv1rFmjwWAoprGxn0CghFhsBKezmJ6eBGZzjO9+t4nm5iCFhaV4PGHGxkYYGuokGg3Q2FjI\nL3+ZZGzMQXa2Gb+/i1TqBFlZIYLBAvr7T6LRFKMoLoqLzajVzWg02YTDA3R3W05Ne/VQXm6hqCiC\nwaAhP381r7/ehN1eyeDgSTZsGOX3v2/CbA7i9++hpsZKZeUy3n77bbKz9SxdqqerS0coBCMjo9TU\nGPH7O0gkwjQ1jeB2t7BmjYbnnotgNofxBf384PcNAHzlz28gHBylocFPJFKE3++etroL+exNjcRp\n6q72WE2WTE1rwnHTTTexb98+7rzzTgB27do1nZcXYs5pb0/S06PnyBE327aR2R/lQlMz29pcHG49\nAcCiMts51xyv4Who8NHc7KGvrwy/30xPz1Gqq0vo7PQADiorI2zYkAOo6OyE//zPKAUFJhRlEJUK\nfL4FZGXFOXlyELPZSTzuR6frwmRqR1HUqFQ2kkkIBKqJRNzk5QVYsWIR110HRmM23d0laDTp4lCv\nN0Y8nkN+vp6FC5NEIlE6Ol7DYPBw7FiE66+fz/BwmFgsTDxei1bbz+iolsHBo7hcCQoKTHR1jVFX\nZwUGAAWATZvs7NkDfr8dgyFKa6sPp7MIcLOn5wiBcIzrSh2sLk9vQS87ugoxe03rp1KlUvHwww9P\n5yWFmNPGZ5RYLOW0twcpKbnwTbGhwc+hQxaGR9NDkaUFk3e72my5bN2axauv9nHiRIKioiQjIwpd\nXSGgm8pKhYqKLF58cZhDh4YYHbXQ3d3B4sUxtNoYGk0Pen0ukUgSi6WdqioTfn8Ov/51LolEATBE\nNJpPIjHEwoVaFi1SUV09isViwu+H3FzYuFHPW2/14fFko9fryc0dYXQ0h0BgEW53F2DE4ViI398F\npMjKUjCbPWg0AcrLr2fv3qMMDS2nv/8oNptCXZ0Vmy1nQkKW3u02nYRUVOgIhdykjCP8+0uHUatU\nfLjuZjo6xigtlR1dhZjN5GuAEDPo9IySYCbBuPBNUcFmSxINBgDQjCns2eM5Z1bLeC8JuFm92sni\nxUVEoydobu5jdNSMVluMx9PFd7/bQ3d3BXr9MFqtC4PBTmOji0ikhooKNWr1AOXl1VgsCUymEN/7\nXj+BgA+1ugzIoqAgSWVlgGg0m3nzVAwMKHR351JeHqaycoRXXlExOhpkbCxFVlYErbaIvLwhQiEv\nen2S4eER3O4RKiogFpuPRtPDqlUqli0z0dnZwchIP8mkQn7+CBUVa2lo6MPvD9Pa2kx1dV4mZjZb\nTua9p1IKH/y7VxhLKdyxsYaastQ76tG43IXZhBCXRhIOIWaRioosXK44Sm8CgL7OGIai9KwWMPCL\nX/Tw4Q8Pn1o5NL3iZldXChihoiKL5mYVTuc8wuEWPJ44x4+r8fvbyMnxsXKlikSinHg8hFoNJpOb\n1auNHDgwwu7dLmIxhd5eJ6OjXej1HZSWmtmyJQ+/P4WirOEXv6insDBAV5cfrfZtRkYq0WjiBIMF\nmEyVjI31kUhcx5Ejb1JZaef48U6s1qU4HFYGB48SjarRaqvR6bJ57bWTJBImoJKFC3NYscJMWVkQ\nr1eF3+/E4+kiHLZhsfgzPR4vvniIaNSDK6uLfW/3Ys0x8H//cgsWs/EdxXwqs4uEEO+cJBxCvAMX\n+nY8fszrDWI0Vk96QzvzGuleDCfD+14CwKTTk53tZt06Hc8840KrXUp9fSf33JPuJdmzB6LRHPz+\nbvz+MKOjdvLzu8nPz+bNN60EAqMEAm7U6hp6evrYsqUEaMPjCVFamgXA4GAMo9HC8HALoVAKs/k6\n7HYzmzdHsdtTDA5G6Ov7OTk5Tvr7+9Dp1tHc7EOl2ohK9TxLl8bQaA5jNgeJRiPY7RG6u8ux2fQE\nAs0MDemZNy9KMmklmewjHjejKDo8nhgVFX62bMmjrq4cgIMH3bz11km02igWSzagYDaHOXLkOC0t\nHtyRGl4dSa/v89gn3zsh2bjcngqp+xDi3SGfMCHegQt9Oz497DFyzjLZk+3bMX7zU2cpkIC6NVYW\nV9uB9Iqa9fWdpzY98+B0ajCbw7z5Zg9Wq5mGhn6GhgYYHh7G6XRgs8HgoJucnBJCoWMkEikslnbm\nzdPR1WUjEIgyMhKitNRAd3cXDkcRipJNONzFkiV5LFhQyqFD+Wi1Nszm11CUIrKyOhkbex29fohY\nrJHKShV//MfzKCnJ4uTJPBTFSDzeismUxcCABp1uERrNQqABm60Ei8XCihVjNDZ6SCZzWbEixbZt\nTgB27+7g9ddLMRhUVFUNsXGjKRMfp9PJye4BGgd/SUpR+N83r+CP6xZMiHd6NdI8vF4f27ZNPeGQ\nug8h3h2ScAjxDlzo2/H4sWXLrOd84z5fknHmNdRqFQBFxXmZx1aunM/KlbB7dztHj55e5nvVqiL2\n7GkkmSzA7VaTn38dLtcBLJZhKivHGB52EQyaCIcX0dbmwWy2YLFAc3MPixen60BWraqgu7sNi6WI\ngoISIIe2tjBVVSFGR30UFWk5dixCIlHE2JiBRYvsqFRh/vAPy6iuHuXgQQWDwU1lZTEVFfmEQnpG\nR/P41a/cBAL93HyzEYMhDKQLQzs73bS3DzE8HD0jKgrl5QlCoX5uvnkeK1eOL+fup62tnxOqbkJj\nEZyOPB76yA3nxNvvD9PZaSY/P3zOMSHElScJhxDvwIW+HV/Kvh3j5+3Z4yESKSKZTE8LdfcPY6k+\nu0ZBNWGZb6+3D4tFi6IolJX1k58/SF6ekb4+J7GYherq4/h8w4RCWUCKNWs07N9/EpWqAI8ngtms\nYDYXkpvbgccTIZkcAIqIx/PQ6ZJ86EN2Bgby8fm6aGoaICfHiUYTwGSah8cT5le/8tLbu4aysgPM\nn6/w6qsuRkezMRgi9PYWMzQUxe2Oc999VbhcfhoafBw+nMRgcDI62pl5V+P1K+vWzWPlyvkT4vj8\n4UPUtw+g12i576bNGHTn/umyWEzMn5+FxWI655gQ4sqThEOIGTSV/VXONJ6IjPvV80Pk/akeIHOd\nujoLNlsEp9NCSUl6CfOlS0t5881j3HHHH1BWFsRsDvPv/95Ke3sKi6UYjSaI1eqmrCzKSy8ZATuJ\nRBYuVye1tVXE480kkznk5y+htHSMtWuzaGwc5OTJUh5/vJmcnGzC4Sp0OjN+f4DRUR+KosLrzUev\nHyQc7iAUcuF2l+DxgNFoA/oIBFKkUrmcOBHOvIe9e8N0do6h1f6eFSvsHDrUSShkwutN4XQuJhQ6\n/f5dLj+/eK2Ph3f/HoB7/vB6Nq2xnjfW6Sm1SczmrMyQk8w6EWL2kIRDiBl0vhqPqRQ35hi1ROIJ\ndNn2zBbs49dJ7xzrz+wzMt7LsWSJBqPRdSoRmZ/Zxv7QoTYSiRx8vhzeeCNMSckiYC9Wazbh8DyS\nySJycvysWWOls7OFzZuzMRhyWbAgzKFDCv39KRQlSjh8GBhFrS4mlVrC6KgHjWaM7GwrhYXZBIMF\ngJpw2IVOZ6Cy0sGCBWHi8WFuuCErM703Egmi1ZaQSIwxPJzHiy96icXMxOMu1q6FZctOx+T1pmH+\n8fm9jKUUPvOBNTzwkbrMsbPjOJ7E7dnDJc06kWmxQrw7JOEQYgadrz7jzCRkPHEYv9k1NPjp6HBg\n1hoYIMKo0onXa6eiQg2cvs6Z19i0yY7NlsRoXEh2tjtz00yfG2PxYivf/a4Lo7ESrbaXRKKFP//z\nHDyeXLq6UhQU+Ni6NZ9QyMD69VZeeMFHLBamokJh5coAr702SiJRgNlcQ0FBglCok0BAYXRUhccz\nzIIFFSxcqFBdXcxLLwVZvnwxanWCmpo8TCYdra0GotFoZt2Qu+8u5Z/+qYVQyEws5qO8PJdDhyLY\nbCXYbDmZ9kfjCb71u1cJjkZYXGTnLzdfN2mx7ZmJxaXOOpFpsUK8OyThEGIGnW/o5Mwb4rk3OwWV\nKkp+tgGAsSwVRmM1odB4z8bpazQ0nCAaVeFyaTJTR6urR9m9O0hnp5vR0Wy2bs1n5cr53H031Nf3\nEY/noijlRKN9XHedjXnzImdsAZ/uHYjHNfj9WioqRohGDdTVraC3t5m+viF0ujyqqsxUV+fy05+G\nUauX0Nt7krvuqiIUsvDe9/Zz4EAMg2EUyOXAgSQDA2X09jZw661vsnVrOSUl89m2LUVHh4OqqgHq\n6ixYLH5AySQJsUSSD+38Ccd6vFiNBv7PH91B58lhIEUkUkRDQ+upSIywbNnEIZZLnXUi02KFeHfI\nJ0yId9nEG6J/ws1uvA6hS5vPwW4XwUSA9vbjrFunA9L1DvX1cdat02Gz5WZuvi5XEoulnNbWbkZH\nHbzxRgSbrZL6+j5WrjxzhouGjg4jAwNgt5fgcrVMaJvTqaGtzUU8HqWpKYxOV4nN1svatUWUlOh4\n5ZUTQA4DAynq6rI4dizOe95jIBQy0dubw549bSQSNiCHVMpBUVGIt98+hMUyn2PHfGzdyqn3eboO\nBSauJBpPjnH3//sVb7S5MWQZeN/8Wgotw6cWORsDgoByTo/O9Pw8hBAzJWvnzp07r2QDIpH4lXz5\nOcNk0kuspmC2xcnl8tPUFEKrjZOTYzznsZISCxUVpsyxnBwjFRUm3IEwv3nzJGadiT9acSNDQ0N4\nPGO8/PIAsdgqPJ5+nM4Ub701QDIZxmyeRyDQR3X1KL29PvT6LoaHB7DbE5SXmzPXj0aHcbsDLF06\nxtBQjOzsAgYHfXg8SXw+DwMDGhQFPB4bLS0h1Oos1q83UVdnpbGxH5crj66uCBZLMWvX5nHjjflY\nrXF2727n1Vc70WhUZGWpKSz0U1qa3vHWZosTCIRxOEx0do5SWqqe8L6bmkJEIkWEQiOUlBq4+//9\nDy8caMds0PIXy6/njluLuP76UtraEkAFeXkJli0zEQqlZ+mMv7eZMtt+p2YridPUXe2xMpn0531c\nejiEmEHnqw+YSs3A6upiAE56h07tkRKioqKUoiIN0ehR1q3TZVYljUZbsdlibN5sp6HBj07nwGpV\nKC4uI5Uy0t4eYbxWxOsdIz+/iNZWN2vWJDlwwE0sFiUSKeDQoR6s1jivv/4mXV1gNBZTUTFARYWT\n9vYkDocKRRkkEhmjtfUgS5ZUUVdnYfduA319TlQqCyZTM+99by4mUxYvvZQE7JSUlPDe9w7jco1l\nNrE7s3ZlfEijfL6aT3zzV/y6sZ08k56f/u12llc5MjGZbCqxFH0KMTdIwiHEDDpffcBUagYWlxdg\n1GnoGQqizdZRXl6J39/NunW5hEImHA4NDgenFhaznHGj9REMDhCPjxKPN9HRAWq1mba2bPx+JxZL\nP9Fo+sb/xhudJBLziMeP0NfXRkGBnSNHuvH5FhCJ6EilDPh8x/nGN/qoqCgEIgSDI2g0iwkGe2lq\ncmCx+Fm3Tkdvbw+RSA+1tfnU1VnZvXuMWMyMy3WC0lId4DiVJAXPqV3ZtMmOvXCMT33reX7V0EaO\nUceD7/tDCgy6CTE5e+jj9NLxIxiNC6XoU4hZThIOIWbQ+eoDplIzoMlSs8LpoP54H96xISrnGdi8\n2X7Ojfrs61RUZPHzn3eQSJSTTPpJJhdx9GgWTqcXn6+HtrZhbr8dwmEXIyMBEgkD8XgCRcmhq6uT\n970vi/37vbz1VpCCgkICgSygml/+8nWqqx3k5BRgNMZQqbyEwwP4/QrbttXicOSxe7cHvb6I9vYY\n1dWj7N/vYe1aJzpdFKNxIaHQ6UJZsznM+Kyb5FiKe779a375eiu52XoefN8fUGy+7qIJxHgs/P5+\nXK7TdS5CiNlJajjmiKt9zG+6zKU4HTrUyc9+NohGE6S4OH9CbUcwOErLyRBv9fRjMOj47J0ryckx\notXGCYVGMJvDtLUlJtSGALS1JXjrLQiFslGrB8nLs+BwDPKBD9g5cmQEk2kVyWSIZDLMiRO5aDQu\nqqvzaGvLI5Uykp0dwWi0sHSpgRtuKGbJkhSHDrnIyclFrVZRXh5i8WITCxZkU1BQQ0lJgpoaC01N\nIRKJQgKBPtavzyEYNOFwlKAoAdat06Eo8VMza3wcO2ZApYqydWsZxmw9n/n2r9ldf4Ico47//psP\nssKZf8H3OG48FoqikJXloK1tmNJS9bTXdMyl36krSeI0dVd7rKSGQ4hZpr4+zvDwUurrj+Jw+Nm9\n24PFUk56FgasKlvL07zJC43ttHYMUl1VeM7iVg0Nrdhs6fqFgYFhXnihH5XKy6JFNjZvLiMcVgOF\nOBx5fPjD6c3f1q3T8eyzPrzeSvLyfGzdOo+BgS5isVxOnFCAahKJPqLRYYxGNWVlETyeMBUVelas\nKKKkZCXR6IkJs0zSw0NBbrzRfmp4xw+k60om1lX4M8uyh6Jx7n38BZ5vaEOXpeELN9UScWvACZs2\n2S+6gNd4LFwuDbt3uzP1ITKsIsTsJD0cc8TVnhFPl7kUJ40mSFtbKxZLCp9vjKwsB93drZhMUFiY\nxKhV8ebJPjyhMNmYuXF1Wea5Pp/n1AyVCAZDNaHQCAcPhunurkavz+W22xwsWmThd7/z0Nubx9hY\niC1bKlm71kZxcT7d3T4SCS0LFyZYsCCPWCxFUZGK5ctVDAz0Egh4UattHDtmIRIpIBrNJpk0YLen\nE49ly0xYrWra25OZ2TZabZyGBj/NzT5KS40sXWo7p7dhfJZMUj/IJ5/4LW+296PP0nB71Z2EXKPM\nm7eIUGiEigpTpgfjYjNRcnKMlJaqGRsbnZFZK3Ppd+pKkjhN3dUeq8l6OCThmCOu9l/Q6TKb4nS+\nKbFnKi7OBzQYDNUkkz6KirIwmVREInba2oZZvz4HrU7Fq0d7aHN5eU+5LdNb0NaWICenEo0mSDLp\nIxqNUVQUJRLxsmBBjPJyNb/9bZC2Ni2trR7icT+LF+dm2uFw6Fi4MJtly8ynlkhPTzfdsqUKnc6A\n3V5OPD5AZeUwBkOA0dFhFKWKaHSEbdvSPS1nTmetqDDR1BTi7bctBAJZGI0pKipM58Sivz9GezDF\n1194leHoKOW2fL7+kY3o4gnq6lIYjapM0jA+RXgqCcSlnHupZtPv1GwmcZq6qz1WMqQixAxK74Ka\nXi1zfOXOqUx/Pb2FfXqmydnDA3d/YA3f/Nmb+EYj/NsvOikpzqehwYffH8ZiCVJXZ828jsHg5ktf\nShdN7tnjAQwMDBzHbLaj06ULTsfbMb7pW0NDBLM5jNd7etXSM5dE7+pKMX++wpYtYV56yY3Fksvu\n3R62bTt3to3ZHCYed+FwgNNZPmG66osvdnPwkIEO1QEO9A8AsKywlK/+xWoq5tnI08m0ViGudtLD\nMUdc7RnxdHk343RmD0Z7e/Kcb/dTGRI4+5v52cMDkXCM3q4QxweHCKZCrCpeQPPxbKLRQhYsULN0\nqS0zvJKXF8DlUqPVxrFa1TQ3+5k/30F+foTS0hiKArm5ZF6ruTmKz2dDUeJEo3E6OhxEIkOsXVuC\nz+fhmWd8NDcn6O7WkJ2dYOHCPIaGwhiNxTQ3D1FTY5owbNLWlqCgoBqHQ8XSpbYJPSANhz08332M\n9hE32iw1f3nDWj78nvehyUri96cm9JTMNvLZmxqJ09Rd7bGarIdD/S63Q4irxukejPS386qq9BDE\n+Df+khLLqamrU//WfvYiVu3tSd639hYKzSZ6fSF+dugAlZXDVFUNZF4nvQDYYlpbDZn2lJRY2LbN\nzvLlBj760XIWLCghGi1h924PLpcfl8vP4OAw0Wjrqeuk93ABBUgXtGq1Sxkc9GO15jMwAEZjNYsX\n56MoXvLzizK72I5zOjVEoyfweoO4XOldbLOz3fhSHn7ScQBvsp+CHCM/f/hD3Hv7dZhMA5mFv7Kz\nZS8TIa520sMxR1ztGfF0eTfjdGYPRkmJhZoaCzU11ndUR3B2XYRWGycSDmE3Gth/soe3e7z8yRYn\nFmM+Vqt6wlRZpzOVmX56dg2EVhunqamf/PwixsYS+P0pDIZq9PoIS5fayM1VYTSmKCxM0taWIC8v\ngM/nwWDw4PdHqa1NYLOli0VrakyMjSXO6bnJyTHi8YwBFXR19RIIwkutR/m7H+8nmkiyrKyI7993\nCysWlhAMjuL3p7BaJy5zDhevfXm3yWdvaiROU3e1x0pqOISYZhdbwOtiS26f7/jZdRHjr+F0amjz\nB3juwFHu+87v2Jz7x6xfHeRjH7NMmB4aCiUnfY1t2+y0t8cy1/b53FitmlP1Jz5AhderYDQWkZ0N\nGzbAT35iAWz4fH3ceac9c83zXR/A6w2iKMM0drn4aVMz/kgIgD9ZvZTblm/FPzAICy+8vPvFal9k\nKXMh5iYZUhFihpw55DKV45PdSMcff+DOVayuLmZ4NMZvPb9lzxtDE27+53u98ccaGnyZFT7Hj2/Z\nUkRJiYWGBh+vvJLiN79J8tJLblyuQ5mhjpUrA5SXt55axTOtocFPY2M2DQ3+Ca/Z3p5kMJLHN199\nne/tP4A/EqK6KJ+nPncrNy9aQCzWlklMxodRzOYwe/Z4JryPiw2xXCyuQojZSRIOIWbIxW6cZx+f\n7EY6foM/9OYIzzywjcrCPMJKkMZYJwePBCdc78waikOHOtm3rw+X6xCgoqdHzzPP+Ojp0Z/1Gios\nFhMDAy7i8VJ+97sgAwPDlJRY+NjHlvOlL60A4IknTnDoUCegEAwO0NycXg00O9uN1ZHkv95s5KGf\n/5DjrkFyDXo+vnErX9h8AycPq4hGS7DZcoHxGTTpxb1CIdMlJw9S8yHE3CSfWCFmyMWGXM4+fvZw\nynjPht8fQqXKBRQsZiO//Oqd3Prgf9LlCfDlnz5PqfM2Vi0ozkx17e3N4Yc/bKa5eQhYiKIMsXWr\nhd273VRVVRAIuNm8+fSeI3V16edt2GDlqaf6sNkWU18/hMORfv3R0X6efHKIrKxKRkYCfPSjFbhc\nHiyWGjr73bx28iDP7mshlhwjS63iE7es5CMbF+Ppz8LrDZKfX4jf301NjY7du+OZ1VTHh4rO3sju\nYkMqU9mLRggx+0jCIcQscfaNdPzGa7GMsGCBgtNpBaAw38T/PPJn/OU3fklji4v3/e2P+eTWWhab\nywgHIwwMKMRiuSSTWpJJBaMxvQKo0RjFYhmmri49c6a3N712yOkt4ou56y5obR1i3TodDQ0+mpqy\nOHJkgLExJyMjfhyO9HDP0uv9PPGLF3mtpYNkKj2zZWVFCd/87B+xuLwAgJoF6bqS9vZYZuO5/Hwr\nfn83N95oP+97drn8eL1BYIRly6zvStyFEO8OSTiEmKVOLwpmPac40pFv4rmHtvN3/7GH7zx/iMef\nP4RZ08KS7CXcvj4Pm1WhvHwMi2UMsNPR4SAY1BGNRhgYGKa9PYnb7aOnpwyvdwCbLYdIJF0ses89\n6WSgq6sDv9+Iw+EA3Cyp1RE0h7nlwWc42OoGQK1SsbK0irqiUj79FwvOaefEhGKy/VVOa29PYjRW\nk53tloJQIa4y7yjh+M1vfsMLL7zAN77xDQAOHz7MI488gkajYf369dx7773T0kghrmaTFYtebOhA\nr9Xw1f+1hfeudvLlf3uZEy4fDSMH6NiTzUdvXsaHbl5CVbHlVEHmAMePB8jPX0B9vRun04nbPcTI\niIfjxwPcfLOa8e3ix61dm09voIvjbg+uiJ+f/76fxFgKgByjjltWV7HcMo/CnGwqKrJOrbTqy6y0\nerapDIWcb4hFCHF1UCmKolzOEx955BH27dvH4sWLMwnHtm3b+Od//mfKysr4xCc+wec//3lqamou\neB2PJ3jB4yLNbs+RWE3BXIzTnj2eU70LbjZtsk84NtUpoKmUwr/9z0H+6RcHGRwJZx5fUlHA6upi\nVi0oothsJOLX4bAlGY3mMG+eiuef95GXX4at0M+8ah2tfT5O9Hpp6fXy+vE+PMORzLXUKhU3sssl\nuwAADxxJREFULi/ng5sW877rqznYEMi0G+Dll8HnC7Bli5pt25zTHKWLm6npsnPxd+pKkDhN3dUe\nK7s957yPX/bXiFWrVnHTTTfx4x//GIBQKEQikaCsLL2j5caNG9m/f/9FEw4hrnWTFYuazWHq6ycW\nWU5GrVZx9wfWcPOqSn7xWh9v9nXy8pEujnUNcaxriKd/e2TC+TlGHQa9hmAkzmh88hkiRRYTNQ4H\niwprqKvO5v1byydt97593VgsdmD0smPxTkxl7xohxJVz0YTj2Wef5amnnprw2K5du7jllltoaGjI\nPBYOhzGbzZn/m0wment7p7GpQlydJisWPXLkOPn5RROKLC+m82SKBbaVLJ9XzLf/zy0cbh/gULub\npvYB3u7y4B2JEAjFCEbjBKOnVzo0GbQUWcxUl1pZWGZlQYmVlQuKWFhqpampi/r6OGWFqXN6Ec5s\n90c/yqljxiuyOJcMxwgxu130k7l9+3a2b99+0QuZTCZCoVDm/+FwmNzc3Is+b7KuF3EuidXUzPU4\n1dUlaW0dYdEiK6FQFtXVVZSVTe2mPf7c6moLZWUWKsqsfODGxRPOSaUUAqFROru9uPsUli42UV4+\n+YwQrdbOqlVFaLVufD7Qasvx+dzU1k6Ms92eQ21t+t+vvOKe9LyZcubrz8S1xcVJnKbuWozVtH0V\nMJvN6HQ6enp6KCsrY+/evVMqGr2ax7Gm09U+5jddroY46fUali7VAKd3Tp3qezr93NPPOXPp8ro6\nS6bHYdijwefJ5UDYjdGonfSaVmscv78DqzV93fF/X6hNZz5nrv88robfqXeDxGnqrvZYTXsNx/k8\n/PDD7Nixg1QqxYYNG1i+fPl0Xl6Iq86lDD1c7jBFe3uSkyfzUBQjNlskMwxSXa3H7z+3duRis2Wm\nUh8hi3MJIc72jhKOuro66urqMv9fvnx5pohUCHFxUyl0HE8EvN4gRmP1JRdFOp0avF4fMILTeTqR\nKCuzoNefu8JnQ0MrNtuFE5uzN22bSiIkm64JcW2T6iohrqDJCh3PvDmPJwIwQjTaSjSq4HJN/aZd\nUmJh27aLnzveFpcrREeHA693YNLnnZkoAVOaHSKzSIS4tknCIcQVNNnQw5k35zNXHJ3Jm/Z4W7ze\nEQKBKDD5Ej1nJ0pTmR0is0iEuLbJJ1+IWWj85pzeTt50xjCE/7Jv2lMd0qirs54aUpl85orUdQgh\nLpVsTy/ELFRSYjnv9u3jjw8MDJ+xXfzUnO4dufBW8OOvIXUWQojpJAmHELOY06khO/vcHo36+jjD\nw0upr49P8swLX8vl8vPKK+5T+6ycn8vlZ88ezwXPEUKIqZKEQ4hZ6kJDIOvW6VCpGonHPeze3TGl\npODMnov29iSh0IV7O6baIyKEEFMhCYcQs9SFbvgrV85nw4YSFGURHR2OS04KnE4NZvOFa0Em610R\nQojLIX9JhJilnE4NDQ0niEZV550GO9n6GlNRUmKhtvbCqx1KkacQYjpJwiHELFVSYsFmm3wa7FTX\n1xBCiNlAhlSEmMVkWEMIcbWQv2JCzGIyrCGEuFpID4cQQgghZpwkHEIIIYSYcZJwCDELzYZFt2ZD\nG4QQVw+p4RBiFprqJm0ul5+GBh+goq7OMqXt4c3mMKGQibq65ITt6S+3DUIIMRXSwyHELDTV2Snt\n7UlOnsyb0uJf4wlEfX2cSKSI1tbYtLRBCCGmQv6SCDELTXV2yqUs/jW+A+26dTpCITfV1Rc+X2bI\nCCGmkyQcQsxSU9lO/lIW/zqdQNgBsNsvvNKoEEJMJxlSEWKWks3ThBBXE0k4hJilpIZCCHE1kb9k\nQsxSUkMhhLiaSA+HEHOArIkhhJjrpIdDiDngzDUx4OLFpEIIMdtID4cQc8CZ9RxSTCqEmIukh0OI\nOWBiPYef9nYpJhVCzC3yF0uIOWY2FJNOZY0QIYQ4kwypCCEumQzrCCEulSQcQohLJmuECCEulfy1\nEEJcstkwrCOEmFsuK+EIhULs2LGDcDhMIpHggQceoLa2lqamJh599FE0Gg3r16/n3nvvne72CiGE\nEGIOuqwhlR/84AesX7+ep59+ml27dvHwww8DsHPnTh577DGeeeYZ3nrrLZqbm6e1sUIIIYSYmy6r\nh+NjH/sYOp0OgGQyiV6vJxQKkUgkKCsrA2Djxo3s37+fmpqa6WutEEIIIeakiyYczz77LE899dSE\nx3bt2sXSpUvxeDx88Ytf5MEHHyQcDmM2mzPnmEwment7p7/FQgghhJhzLppwbN++ne3bt5/zeEtL\nCzt27OD+++9nzZo1hEIhQqFQ5ng4HCY3N3d6WyuEEEKIOemyhlTa2tq47777+OY3v8miRYsAMJvN\n6HQ6enp6KCsrY+/evVMqGrXbcy6nCdckidXUSJymTmI1NRKnqZE4Td21GCuVoijKpT7pnnvuoaWl\nhdLSUhRFITc3l8cff5zDhw/z6KOPkkql2LBhA/fdd99Fr+XxBC+r4dcauz1HYjUFEqepk1hNjcRp\naiROU3e1x2qyZOqyejieeOKJ8z5eW1vLj3/848u5pBBCCCGuYrLSqBBCCCFmnCQcQgghhJhxknAI\nIYQQYsZJwiGEEEKIGScJhxBCCCFmnCQcQgghhJhxknAIIYQQYsZJwiGEEEKIGScJhxBCCCFmnCQc\nQgghhJhxknAIIYQQYsZJwiGEEEKIGScJhxBCCCFmnCQcQgghhJhxknAIIYQQYsZJwiGEEEKIGScJ\nhxBCCCFmnCQcQgghhJhxknAIIYQQYsZJwiGEEEKIGScJhxBCCCFmnCQcQgghhJhxknAIIYQQYsZJ\nwiGEEEKIGScJhxBCCCFmnCQcQgghhJhxknAIIYQQYsZJwiGEEEKIGScJhxBCCCFmnOZynhSNRvnC\nF77AyMgIOp2Ov//7v6ewsJCmpiYeffRRNBoN69ev5957753u9gohhBBiDrqsHo7//u//ZunSpfzo\nRz/i/e9/P9/73vcA2LlzJ4899hjPPPMMb731Fs3NzdPaWCGEEELMTZfVw3HXXXehKAoALpeL3Nxc\nQqEQiUSCsrIyADZu3Mj+/fupqamZvtYKIYQQYk66aMLx7LPP8tRTT014bNeuXSxdupS77rqL1tZW\nvv/97xMOhzGbzZlzTCYTvb29099iIYQQQsw5F004tm/fzvbt28977KmnnqKjo4NPfvKT7N69m1Ao\nlDkWDofJzc2dvpYKIYQQYs66rCGV73znOzgcDm677Tays7PJysrCZDKh0+no6emhrKyMvXv3Tqlo\n1G7PuZwmXJMkVlMjcZo6idXUSJymRuI0dddirFTKeDHGJfB6vdx///3EYjEURWHHjh2sWLGCw4cP\n8+ijj5JKpdiwYQP33XffTLRZCCGEEHPMZSUcQgghhBCXQhb+EkIIIcSMk4RDCCGEEDNOEg4hhBBC\nzDhJOIQQQggx4y5rWuw7kUql2LVrF2+//TbxeJzPfvaz3HjjjbIPywW0t7dzxx13sH//fnQ6ncTq\nLKFQiB07dhAOh0kkEjzwwAPU1tZKnM5DURR27txJS0sLOp2ORx55hHnz5l3pZs0ayWSSL3/5y/T1\n9ZFIJPjUpz7FggUL+NKXvoRaraa6upqvfOUrV7qZs4bX6+WDH/wgP/jBD8jKypI4ncd3vvMdXn75\nZRKJBB/+8IdZu3bttRsn5V323HPPKQ8//LCiKIridruVp556SlEURbntttuUnp4eRVEU5e6771aO\nHz/+bjdtVgoGg8onPvEJZf369UosFlMURWJ1tm9961uZ36OOjg7l9ttvVxRF4nQ+L730kvKlL31J\nURRFaWpqUj796U9f4RbNLj/96U+VRx99VFEURRkeHlY2b96sfOpTn1IaGxsVRVGUhx56SPnNb35z\nJZs4ayQSCeUzn/mMsnXrVqWjo0PidB5vvPGG8qlPfUpRFEUJh8PKt7/97Ws6Tu/6kMrevXspLCzk\nk5/8JA899BBbtmyZdB8WAQ899BCf//znMRgMABKr8/jYxz7GnXfeCaS/oer1eonTJA4ePMimTZsA\nqK2t5ejRo1e4RbPLLbfcwuc+9zkAxsbGyMrK4tixY6xZswaAG264gfr6+ivZxFnja1/7Gn/2Z39G\nYWEhiqJInM5j7969LFy4kHvuuYdPf/rTbN68+ZqO04wOqZxvHxar1Yper+fJJ5+ksbGRBx54gG98\n4xvX/D4s54tVSUkJt956K4sWLcpslnet71lzob19PB4PX/ziF3nwwQev+ThNJhQKkZNzeoVDjUZD\nKpVCrZZyLgCj0Qik4/S5z32Ov/qrv+JrX/ta5rjJZCIYDF6p5s0azz33HDabjQ0bNvCv//qvQHq4\nfJzEKc3v9+NyuXjyySfp6enh05/+9DUdpxlNOM63D8vnP/95tmzZAsDatWvp7OzEbDZf8/uwnC9W\nW7du5dlnn+UnP/kJQ0NDfPzjH+df/uVfrulYTba3T0tLCzt27OD+++9nzZo1hEKhazpOkzGbzYTD\n4cz/Jdk4V39/P/feey8f+chHuPXWW/nHf/zHzDH5PUp77rnnUKlU7Nu3j5aWFu6//378fn/muMQp\nLT8/H6fTiUajobKyEr1ez8DAQOb4tRand/0vzerVq3n11VcBaG5upqSkZMI+LIqisHfvXlavXv1u\nN23WefHFF/nhD3/I008/TUFBAd///vcxm80Sq7O0tbVx33338fWvf52NGzcCSJwmsWrVqsznr6mp\niYULF17hFs0u44n9X//1X3P77bcDsHjxYhobGwF47bXX5PcI+NGPfsTTTz/N008/TU1NDf/wD//A\npk2bJE5nWb16NXv27AFgYGCAaDTKe97zHhoaGoBrL07v+iyVP/3TP2Xnzp3ccccdADz88MMA7Ny5\nkx07dmT2YVm+fPm73bRZTaVSZYZVHn74YYnVGR577DHi8TiPPPIIiqKQm5vL448/Lr9T53HTTTex\nb9++TM3Lrl27rnCLZpcnn3ySkZERnnjiCR5//HFUKhUPPvggX/3qV0kkEjidTm6++eYr3cxZ6f77\n7+dv//ZvJU5n2Lx5MwcOHGD79u2ZGWKlpaX8zd/8zTUZJ9lLRQghhBAzTgZvhRBCCDHjJOEQQggh\nxIyThEMIIYQQM04SDiGEEELMOEk4hBBCCDHjJOEQQgghxIyThEMIIYQQM04SDiGEEELMuP8P7l9R\nI6qdqZgAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1962e39f320>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "from numpy.random import multivariate_normal\n",
    "from filterpy.stats import (covariance_ellipse, \n",
    "                            plot_covariance_ellipse)\n",
    "\n",
    "mean = (5, 3)\n",
    "P = np.array([[32, 15],\n",
    "              [15., 40.]])\n",
    "\n",
    "x,y = multivariate_normal(mean=mean, cov=P, size=2500).T\n",
    "plt.scatter(x, y, alpha=0.3, marker='.')\n",
    "plt.axis('equal')\n",
    "\n",
    "plot_covariance_ellipse(mean=mean, cov=P,\n",
    "                        variance=2.**2,\n",
    "                        facecolor='none')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## The Algorithm"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "As I already stated, when the filter is initialized a large number of sigma points are drawn from the initial state ($\\mathbf{x}$) and covariance ($\\mathbf{P}$). From there the algorithm proceeds very similarly to the UKF. During the prediction step the sigma points are passed through the state transition function, and then perturbed by adding a bit of noise to account for the process noise. During the update step the sigma points are translated into measurement space by passing them through the measurement function, they are perturbed by a small amount to account for the measurement noise. The Kalman gain is computed from the \n",
    "\n",
    "We already mentioned the main difference between the UKF and EnKF - the UKF choses the sigma points deterministically. There is another difference, implied by the algorithm above. With the UKF we generate new sigma points during each predict step, and after passing the points through the nonlinear function we reconstitute them into a mean and covariance by using the *unscented transform*. The EnKF keeps propagating the originally created sigma points; we only need to compute a mean and covariance as outputs for the filter! \n",
    "\n",
    "Let's look at the equations for the filter. As usual, I will leave out the typical subscripts and superscripts; I am expressing an algorithm, not mathematical functions. Here $N$ is the number of sigma points, $\\chi$ is the set of sigma points."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Initialize Step"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "$$\\boldsymbol\\chi \\sim \\mathcal{N}(\\mathbf{x}_0, \\mathbf{P}_0)$$\n",
    "\n",
    "This says to select the sigma points from the filter's initial mean and covariance. In code this might look like\n",
    "\n",
    "```python\n",
    "N = 1000\n",
    "sigmas = multivariate_normal(mean=x, cov=P, size=N)\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Predict Step"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "$$\n",
    "\\begin{aligned}\n",
    "\\boldsymbol\\chi &= f(\\boldsymbol\\chi, \\mathbf{u}) + v_Q \\\\\n",
    "\\mathbf{x} &= \\frac{1}{N} \\sum_1^N \\boldsymbol\\chi\n",
    "\\end{aligned}\n",
    "$$\n",
    "\n",
    "That is short and sweet, but perhaps not entirely clear. The first line passes all of the sigma points through a use supplied state transition function and then adds some noise distributed according to the $\\mathbf{Q}$ matrix. In Python we might write\n",
    "\n",
    "```python\n",
    "for i, s in enumerate(sigmas):\n",
    "    sigmas[i] = fx(x=s, dt=0.1, u=0.)\n",
    "\n",
    "sigmas += multivariate_normal(x, Q, N)\n",
    "```\n",
    "\n",
    "The second line computes the mean from the sigmas. In Python we will take advantage of `numpy.mean` to do this very concisely and quickly.\n",
    "\n",
    "```python\n",
    "x = np.mean(sigmas, axis=0)\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can now optionally compute the covariance of the mean. The algorithm does not need to compute this value, but it is often useful for analysis. The equation is\n",
    "\n",
    "$$\\mathbf{P} = \\frac{1}{N-1}\\sum_1^N[\\boldsymbol\\chi-\\mathbf{x}^-][\\boldsymbol\\chi-\\mathbf{x}^-]^\\mathsf{T}$$\n",
    "\n",
    "$\\boldsymbol\\chi-\\mathbf{x}^-$ is a one dimensional vector, so we will use `numpy.outer` to compute the $[\\boldsymbol\\chi-\\mathbf{x}^-][\\boldsymbol\\chi-\\mathbf{x}^-]^\\mathsf{T}$ term. In Python we might write\n",
    "\n",
    "```python\n",
    "    P = 0\n",
    "    for s in sigmas:\n",
    "        P += outer(s-x, s-x)\n",
    "    P = P / (N-1)\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Update Step"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In the update step we pass the sigma points through the measurement function, compute the mean and covariance of the sigma points, compute the Kalman gain from the covariance, and then update the Kalman state by scaling the residual by the Kalman gain. The equations are\n",
    "\n",
    "$$\n",
    "\\begin{aligned}\n",
    "\\boldsymbol\\chi_h &= h(\\boldsymbol\\chi, u)\\\\\n",
    "\\mathbf{z}_{mean} &= \\frac{1}{N}\\sum_1^N \\boldsymbol\\chi_h \\\\ \\\\\n",
    "\\mathbf{P}_{zz} &= \\frac{1}{N-1}\\sum_1^N [\\boldsymbol\\chi_h - \\mathbf{z}_{mean}][\\boldsymbol\\chi_h - \\mathbf{z}_{mean}]^\\mathsf{T} + \\mathbf{R} \\\\\n",
    "\\mathbf{P}_{xz} &= \\frac{1}{N-1}\\sum_1^N [\\boldsymbol\\chi - \\mathbf{x}^-][\\boldsymbol\\chi_h - \\mathbf{z}_{mean}]^\\mathsf{T} \\\\\n",
    "\\\\\n",
    "\\mathbf{K} &= \\mathbf{P}_{xz} \\mathbf{P}_{zz}^{-1}\\\\ \n",
    "\\boldsymbol\\chi & = \\boldsymbol\\chi + \\mathbf{K}[\\mathbf{z} -\\boldsymbol\\chi_h + \\mathbf{v}_R] \\\\ \\\\\n",
    "\\mathbf{x} &= \\frac{1}{N} \\sum_1^N \\boldsymbol\\chi \\\\\n",
    "\\mathbf{P} &= \\mathbf{P} - \\mathbf{KP}_{zz}\\mathbf{K}^\\mathsf{T}\n",
    "\\end{aligned}\n",
    "$$"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This is very similar to the linear KF and the UKF. Let's just go line by line.\n",
    "\n",
    "The first line,\n",
    "\n",
    "$$\\boldsymbol\\chi_h = h(\\boldsymbol\\chi, u),$$\n",
    "\n",
    "just passes the sigma points through the measurement function $h$. We name the resulting points $\\chi_h$ to distinguish them from the sigma points. In Python we could write this as\n",
    "\n",
    "```python\n",
    "sigmas_h = h(sigmas, u)\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The next line computes the mean of the measurement sigmas.\n",
    "\n",
    "$$\\mathbf{z}_{mean} = \\frac{1}{N}\\sum_1^N \\boldsymbol\\chi_h$$\n",
    "\n",
    "In Python we write\n",
    "\n",
    "```python\n",
    "z_mean = np.mean(sigmas_h, axis=0)\n",
    "```\n",
    "    \n",
    "Now that we have the mean of the measurement sigmas we can compute the covariance for every measurement sigma point, and the *cross variance* for the measurement sigma points vs the sigma points. That is expressed by these two equations\n",
    "\n",
    "$$\n",
    "\\begin{aligned}\n",
    "\\mathbf{P}_{zz} &= \\frac{1}{N-1}\\sum_1^N [\\boldsymbol\\chi_h - \\mathbf{z}_{mean}][\\boldsymbol\\chi_h - \\mathbf{z}_{mean}]^\\mathsf{T} + \\mathbf{R} \\\\\n",
    "\\mathbf{P}_{xz} &= \\frac{1}{N-1}\\sum_1^N [\\boldsymbol\\chi - \\mathbf{x}^-][\\boldsymbol\\chi_h - \\mathbf{z}_{mean}]^\\mathsf{T}\n",
    "\\end{aligned}$$\n",
    "\n",
    "We can express this in Python with\n",
    "\n",
    "```python\n",
    "P_zz = 0\n",
    "for sigma in sigmas_h:\n",
    "    s = sigma - z_mean\n",
    "    P_zz += outer(s, s)\n",
    "P_zz = P_zz / (N-1) + R\n",
    "\n",
    "P_xz = 0\n",
    "for i in range(N):\n",
    "    P_xz += outer(self.sigmas[i] - self.x, sigmas_h[i] - z_mean)\n",
    "P_xz /= N-1\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Computation of the Kalman gain is straightforward $\\mathbf{K} = \\mathbf{P}_{xz} \\mathbf{P}_{zz}^{-1}$.\n",
    "\n",
    "In Python this is\n",
    "\n",
    "```python\n",
    "K = np.dot(P_xz, inv(P_zz))```\n",
    "\n",
    "Next, we update the sigma points with\n",
    "\n",
    "$$\\boldsymbol\\chi  = \\boldsymbol\\chi + \\mathbf{K}[\\mathbf{z} -\\boldsymbol\\chi_h + \\mathbf{v}_R]$$ \n",
    "\n",
    "Here $\\mathbf{v}_R$ is the perturbation that we add to the sigmas. In Python we can implement this with\n",
    "\n",
    "```python\n",
    "v_r = multivariate_normal([0]*dim_z, R, N)\n",
    "for i in range(N):\n",
    "    sigmas[i] += dot(K, z + v_r[i] - sigmas_h[i])\n",
    "```\n",
    "\n",
    "\n",
    "Our final step is recompute the filter's mean and covariance.\n",
    "\n",
    "```python\n",
    "    x = np.mean(sigmas, axis=0)\n",
    "    P = self.P - dot3(K, P_zz, K.T)\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Implementation and Example"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "I have implemented an EnKF in the `FilterPy` library. It is in many ways a toy. Filtering with a large number of sigma points gives us very slow performance. Furthermore, there are many minor variations on the algorithm in the literature. I wrote this mostly because I was interested in learning a bit about the filter. I have not used it for a real world problem, and I can give no advice on using the filter for the large problems for which it is suited. Therefore I will refine my comments to implementing a very simple filter. I will use it to track an object in one dimension, and compare the output to a linear Kalman filter. This is a filter we have designed many times already in this book, so I will design it with little comment. Our state vector will be\n",
    "\n",
    "$$\\mathbf{x} = \\begin{bmatrix}x\\\\ \\dot{x}\\end{bmatrix}$$\n",
    "\n",
    "The state transition function is\n",
    "\n",
    "$$\\mathbf{F} = \\begin{bmatrix}1&1\\\\0&1\\end{bmatrix}$$\n",
    "\n",
    "and the measurement function is\n",
    "\n",
    "$$\\mathbf{H} = \\begin{bmatrix}1&0\\end{bmatrix}$$\n",
    "\n",
    "The EnKF is designed for nonlinear problems, so instead of using matrices to implement the state transition and measurement functions you will need to supply Python functions. For this problem they can be written as:\n",
    "\n",
    "```python\n",
    "def hx(x):\n",
    "    return np.array([x[0]])\n",
    "\n",
    "def fx(x, dt):\n",
    "    return np.dot(F, x)\n",
    "```\n",
    "\n",
    "One final thing: the EnKF code, like the UKF code, uses a single dimension for $\\mathbf{x}$, not a two dimensional column matrix as used by the linear kalman filter code.\n",
    "\n",
    "Without further ado, here is the code."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false,
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhwAAAD+CAYAAABx7dzcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xd0VMUewPHv1pRNr6SRQCq9C0gREBSkWQDposITREAQ\npEhHQRFF7CiCSgcRUBDpHaQZCC0hpJDee9123x+RRR4thGQDvPmc49G9e3fu7O/E3d/OnfmNTJIk\nCUEQBEEQhCokr+4OCIIgCILw+BMJhyAIgiAIVU4kHIIgCIIgVDmRcAiCIAiCUOVEwiEIgiAIQpUT\nCYcgCIIgCFWu3AnHuXPnGDJkCACXL19m0KBBDB06lOHDh5OVlQXAhg0beOmll+jfvz8HDhyokg4L\ngiAIgvDoUZbnpGXLlrF161Y0Gg0A8+fPZ+bMmQQHB7N+/Xq+//57Xn/9dVauXMnmzZspKSlhwIAB\ntGnTBpVKdcd209PzK+dd3IWjozXZ2UVVfh1BxNrcRLzNR8TavES8zacqYu3qanvb4+Ua4fD19eWr\nr74yPV68eDHBwcEA6PV61Go1YWFhNGvWDKVSiY2NDX5+fkRERFRC1x+MUqmo7i783xCxNi8Rb/MR\nsTYvEW/zMWesy5VwdOnSBYXiRqdcXFwA+Pvvv1mzZg3Dhg2joKAAW9sbWY21tTX5+VU/giEIgiAI\nwsOvXLdUbuePP/5g6dKlfPfddzg6OmJjY0NBQYHp+cLCQuzs7O7ahqOjtVmyqzsN7wiVT8TavES8\nzUfE2rxEvM3HXLGuUMKxdetWNmzYwMqVK01JRcOGDfnss8/QarWUlpYSHR1NYGDgXdsxxz06V1db\ns8wVEUSszU3E23xErM1LxNt8qiLWd0pg7jvhMBqNzJ8/H09PT0aPHo1MJuOJJ57grbfeYsiQIQwc\nOBBJkpgwYQJqtfqBOy4IgiAIwqNPVp27xZojgxWZsvmIWJuXiLf5iFibl4i3+TzUIxyCIAiCIFQ+\ng8FAVNRVAPz9A25arPE4EJVGBUEQBKGaxcRE07Hjk7Rt24K2bVvw1FOtiI6+Wt3dqlRihOM2QkPP\nMHPmVGrVqs31O06Ojk7Mnbvgtuf37v0sW7fuBODatVimTJnApEnTkMlkN7Ujk8no0qUrPXs+b7b3\nIgiCIDzcDAYDQ4f2JyIiHCcnJwCuXIlg8OCXOXToBErl4/FV/Xi8iyrQrFkLZs/+oJxnywCIjo5i\nxozJzJgxl7p16xMaeuY+2xEEQRD+30RGXiEiIhxnZ2dOnQoDZLRs2ZirVyMJD79M/foNqruLleKh\nTjg27LvKqfC0B2pDoZBhMNyYF9sixI1+nQLu+brbzaUdM+YNAgODiI6OoqioiHnzPsTdvQYgcfVq\nJDNmTOaDDxZSu3bAXdsRBEEQhOuMRiMAMpkcuVyBTCZDLi+b8fA4fYc81AlHdfr779OMHTvSdCuk\nVas2yGQy6tatz9ix7/Ddd1+zZ89OBg16haKiIubPn4NSqbyluur/trNkyTfIZLJqeleCIAjCwyYo\nKJjatf2Jjo6iZcvGyGQyUlNT8POrRXBwSHV3r9I81AlHv04B5RqNuJuKLvm53a2Q48ePEBRUtoeM\nm5s72dlZ/zwj48MPPyEnJ5v33pvM99//hIODwx3bEQRBEITrlEolP/+8jsGD+xEbGwOAr68fP/+8\n7rGqZ/VQJxzV6c7DWLeOTlhbW+Hm5o6bmzsvvdSXuXOn8+mnX1ZtBwVBEITHRlBQMEeOnCI8/BKS\nJFGnTr3HKtkAkXDcUWjoGcaOHWl6LEkSWq32DmffSEL69x/MyZMn+PHHZTRq1KSKeykIgiA8LtRq\nNQ0bNq7ublQZUWlUqDQi1uYl4m0+ItbmJeJtPuasNCoKfwmCIAiCUOVEwiEIgiAIQpUTCYcgCIIg\nCFVOJByCIAiCIFQ5kXAIgiAIwmNu/vy5REZeqdbKpSLhEARBEISHlMFgIC0tDZ1Od9fzYmNjSE9P\nNz2eM2cG+/btITc3i6Sky4SEuFJQcJa4uLNV3eU7EnU4biM09Axbtmxizpz5AOzfv4cVK76nTp16\nhIWdxdXVzVSq/M03xxISUreaeywIgiA8bn7/fQsTJ44jOzsbS0tL3n33PUaPHotMJuOnn5bj7e1N\np05dKCkp4aefvic4uDZt2zZFkorp3DkIe/t0tNrD2Nqq6dLFD4C8vOp7PyLhuIPr+53s3v0n69ev\nYcmSb9m8eSP9+w+md+8Xq7l3giAIwuPs9OmT/Oc/r2IwGNBobCgsLGDu3BloNJY8+2w7PD0VWFun\nce3a76jVBsaOfRKlUgGkAuDg4FK9b+A2HuqEY+BAK/bsqYwu3ihC0rmznjVriu/5CkmS2LnzD375\nZT1LlnyNRmNjOi4IgiAIVWnNmpUYDAYGDRrK1KmT+P77+WzbtpaLF5fTp09NWrVy+9fZqmrr5/14\nqBOO6hQWdpbMzAzy8/PR6/Wm4+vXr2bfvt0A1K4dwNtvT6yuLgqCIAiPqfz8fJRK8PXVo1CE8eST\nnpw8CZKkv/eLbyMq6iqXLl0kM9OJ5593xsenZiX3+N4e6oSjPCMR91LRsq0uLi4sXvwVv/22mTlz\nZvDJJ58DiFsqgiAIQpXYvv13mjVrjkxWQMeOPqSmwuHD68nKukRERAQA9evXv+92d+36k5UrVyBJ\nEBkJixf/xJo1v9CqVevKfgt3Ve5VKufOnWPIkCEAxMXFMXDgQAYPHsycOXNM52zYsIGXXnqJ/v37\nc+DAgUrvrDl5efmgUql46aV+qNUqfvrph+rukiAIgvAYS0i4xvnzO1GpwunRozXduj2HXm/g3Lmz\nlJQU07RpM3r3fum+2kxNTWHVqh+RJGjXrj1NmjSloCCfN98cbvYpAuUa4Vi2bBlbt25Fo9EAsGDB\nAiZMmEDz5s2ZNWsWe/bsoXHjxqxcuZLNmzdTUlLCgAEDaNOmDSrVo3Fv6W6mTJnJa68NwsvLm86d\nn63u7giCIAiPAaPRSHj4ZerWrUdycgTPPeeLg8ONr+XBg1+hffsOpKQk4+joiL9/IHL5/VWzuHbt\nGkajRL169XnjjdFkZ7tx9Gh3EhLiycjIwM3NrrLf1h2VK+Hw9fXlq6++4t133wXg4sWLNG/eHID2\n7dtz9OhR5HI5zZo1Q6lUYmNjg5+fHxERERUa/qluTZo0o0mTZqbHDg4O/Prr9mrskSAIgvA4SElJ\nxs3NHblcjlarZcKEt1i8eCru7qXY2Nz6lVyzpi81a/pW+Ho2NmUDBQkJ8SQkxHP5cgxFRUUoFAps\nbW+/q2tVKVfC0aVLFxITE02P/z0Mo9FoKCgooLCw8KbOW1tbk58vthcWBEEQHn9Go7Fcow+vvDKQ\nBQtm4+qqxmDIZuXKt7G0NFBVUyqDgkLw9/cnKiqKKVMmEhlZdvy110ZgaWlZJde8kwq9w38HtbCw\nEDs7O2xsbCgoKLjl+N04Olr/s264arm6mjeL+38mYm1eIt7mI2JtXo9KvM+dO8frr79OaGgoTk5O\nTJo0iUmTJplqOUmSREpKCu7u7ly7FsbGjVOwszMik5UC1mbp49y5c1i6dCmhoWG4uKh44403mD17\nNkplWQpgrlhXKOGoW7cup06dokWLFhw6dIhWrVrRoEEDFi9ejFarpbS0lOjoaAIDA+/aTnZ2UYU6\nfT8qukpFuH8i1uYl4m0+Itbm9ajEOyUlmU6dOpGVlQVARkYGkydPRqs18sYbowH4669jLF36JTNn\nvoaTUylyuZz8/BKz9dFohF37vNh3dBV//+VFvcbpjBunITu7bBVoVcT6TglMhRKOyZMnM2PGDHQ6\nHf7+/nTt2hWZTMaQIUMYOHAgkiQxYcIE1Gr1A3VaEARBECpLQUEBarW60r6btmzZRFZWFi1btmbN\nmo1s2fIr77wzli+++IxBg4ai0djg4+PA3Ln9cHDQYa7tyxISEli1JorTYV25cqkDhbk3EgB3j2JA\nY5Z+/K9yJxxeXl6sW7cOAD8/P1auXHnLOX379qVv376V1ztBEARBuIvExARmzpzGhQthODu7MGbM\neLp1637TOZGRVxg1ajhhYWdRKBT07z+IBQsWPfAchuvTCGrW9EWr1dGq1ZMAFBbmcvHiftzdLbC3\n16FUWjzQde5FkiA11YojR13ZtNWBc6FdKS26MaXByjaPZ5+LZ9TwWNxcbIHqKXv+UBf+EgRBEIQ7\nyc7OolevrsTHxwEQExPNK68MYMWK1XTv3hOA/Pw8+vV7nsTEBBQKBQaDgdWrf0Yul5sKOlZUixYt\nAfjll/WkpESi1yfSqBG0a1efunVVgBF48HmK+fkqEhOtiY2zIjpWQ1yCNcnJVqSlWpOVYU1OphV6\n/c3XsbAuJDAkGkvVQpSy1fTuOpKA2h3Iy6u+uTEi4RAEQRAeSRs3riM+Po66devz+edfs2nTRr75\n5gsWL/7YlHAcOLCfxMQEgoND2LFjH1FRkXTp8hRr167io48+NU2cvB+7d//JqVMnGTHiVcaO7cfB\ngxvQ6c4gk0FAgCfDhr1eKe8vPt6a2e/X5+BeXyRJdtdzVRZa3GtmY2O7F0vFD3y5uDceHjXYsKGA\n334rm7ha3UTCIQiCIDyS0tPTAXjmma40bNgYa2sN33zzBenpaaZziovLFid4enphY2NDrVq1kclk\n6PV6dDpduROOpKREPD290Ov1eHra06lTbRSKs4wb9xLPPdeA2NgYbG1tady4ablv1RQUFJCTk42z\nszNWVjdWrGRkWLDo07ps+sUfo0GOTG7EwSUfB6diXFyLqVGjCE/PIry9i/D1KcLPtxAXRwOWapg1\nayZRUVH8/ruBFi1a8tdfxwBwdnYuV5+qkkg47uLixQt8++0XfPHF0uruiiAIgvA/6tSpC8CKFctQ\nq9Xs3bv7puMATzzRCpVKxf79e+nTpzcJCXFIkkSzZi2wsrIynVdQUIBGozEtZ01PT8fJyQmFQkFB\nQQH9+7/Ijz9+jkKRjo+PHJnsxhd4YGAQgYFB5e63JEls3Liebdu2YDRKqFRKBg16hY4dn+XLL+vw\n7dIQtKVKQKJOszjmzrhA0wYF92wXoG/fAXz88XwOHTrIoUMHAfDxqUmbNu3K3b+qIpOqcb91cyx7\nquiSnzVrfmbnzj+wsrLm22+XV0HPHj+PylK2x4WIt/mIWJtXeeNtNBp57bUh/PHH76ZjLi6ubNu2\nk9q1A0zHNmxYy/jxb6HT6QDw9fWjceOmDB36KoGBXkiSnqFDRzB//sc0bdoChUJB587tWLhwPm5u\nGozGAmSybBwcKmd1y969u1mxYtk/79WV9PR09HpbSnWnCTtblrh4BSUzfvx5nu+Sg+zud1NuERER\nzo4d2ykoyMfX15cXX+yLRmMDQF5eDby8mpjOfeiXxf4/8PLyYf78RcybN9N07NixI2zZsonw8Ev4\n+dWiU6fOPP98n2rspSAIwv8vuVzODz/8zMaN67h48TwuLq707z8Id/caN53Xr98A7OzsOXculDp1\n6lCvni+pqVG4u6dhZVWAXC7nl1+mUFycRELC70iSgp9/noCNTT4y2fWRhcor83Dy5HEAhg17nc6d\nn+GLL/bw3bJpFBX4obLUMvStk0x4PQnLCl4yODiE4OCQSutvZXnoE46FC+cD8O670yr0ePbs2RQW\nlpoel9dTT3UkJSXZ9Dg1NYWjRw+xcOFiDh8+gNEo8dRTHSv6tgRBEIRKcH2Z671YWFjg4WFNs2Zy\n7OzScHO7+Ve4XC5Ho7FCYypRUXVfj0ajEQBLS0sOH3Hjm6VLKC22wtY5lw8XHaNr+8dzNO2hTzj+\nN1G438ezZ8+ulOGiHTu20adPfwByc3Px9PR64DYFQRCEqrNhw1peeKEP6elR+PuX0Lz5w/Grv169\nVhw7Xp933+tOWmJ7JEmOq28M701eRdf2wdXdvSrz0CccD4uCggJq1PAA4OLF83Tt2v0erxAEQRDM\nafDgfsyYMRcfn5oUF+dz8uRhAgPV+PlZVOq+XdenPsruY3JFeroF23a58ccOL86feQGdVmV6zq/e\nHt4Yspyezz5fof7Ex8cRGxuDRmND/foNHtoq3yLhKKeePZ9nz56dAPTp079Ca7cFQRCEypOVlUlW\nVhZ+frVITj7PqFHdkMvPk5t7GQsLJbNm9azU6+l0OtavX8PBg/vR6XTUrVuP4cPfIDo6mNWr/TEa\nZbi6luDqWoKLSzGSJOPocRf+OuFKWvLNt3BcvXJo2iKadm1P06WdAReX8icbBQUF6HRa7O0d+OOP\nbaxfv5rryz98fHx49933cHR0vOV11bhGBBCrVIRKJGJtXiLe5iNibV53irderyczMxN3d3cAfv11\nI1FRF+nb9wmcnO5vxKEili9fxr59u02P8/KeJD3jI9JS297ztQqlHnuXKOztduHm9AMzpj9/0/Ld\n8iguLuL777/l5MkTADg7O5GZmYVMBk2aNCcuLpaMjAyaNWvB+PETb3ptePg1PvpoNYsXf46Pj1il\nIgiCIPyf0+v1pv9OTU1l9Oj/sHr1BoqL84mPj2f8+LGsWbMGg6GUpk0d6dixOUpl1SYaUDa6ceDA\nXoxGFc89t5wtWztw8aI3AEq1jlpNotE4FlBaZEFpoSWlhRYYjXJ8a6ehLfoBC+NKnu38JKmpqZw9\ne46lSxNZvPjL+0qSli//npMnTyCXy1CrLcjMLNul9skn2zJq1BhSUpKYOHE84eGXAMjIyMHJqWxP\nFWdnb6ZPn4+nZ8NKjkz5iYRDEARBeCgUFxfTqlUTzpw5TmpqInp9Dm+88Rzp6X9gaanA01PGxo2T\nkcnCkMtl2NjIqIy9Su4lN1fFn3+6ER65mZyszpw4UVYVVKnW4dfoIs/1OMHgXo442oBWD6U60OpA\nqYDL54/w9ddf0aRJU4YMeRWDwcDrrw8lIyODwsJCbGxsytUHnU7HiRNly2k//PATXFxceeedsWRn\nZxMTE41WqyUuLh4oW5EDMH78F7z11ivUrfsUfn5ByOXm2a32TkTCIQiCIFSb114bwvTps7G2ltDp\n0lm9ei6SdBo7u7KRjlatat53mzqdjkuXLlBYWICfX+1yrSosLlaQkKAhIUHDtWs2xMTYEHnVlphY\nW9JTrW/ay8TWJZeagRew08xBUXqAAc98gvM/m7NaqrmpfoadXdkT0dFRREVFkpycjF6vR6VS3lTp\n9F6MRiNGo4RMVtamSqXC2dmF7OxskpKSGD58KDqdRGkpvPhiJzIzjXz55TI8PUOqPdG4TiQcgiAI\ngtnMmTODZ57pSlBQACUlqQwa1AlJ+hsrK2s0GjkODrao1SpKSvT3buw2cnJy+OijD0w7yMrlMgYO\nHGJaWShJkJCg4e+/nQkNdeb8eUfi4mzIyrrz/icyuREX7wzqNooF41fo835BMhQiL5Yz7NXhuLm5\n3/G1ISF1CQgI5OrVSGbNmm463r17TxSK8o/OWFhY4OcXyMmTkUyb9i4ODg5ERESTmyunfn03UlJS\nKC5WotHUomPHsXh61n/oFjc8XL0RBEEQHiuJiQlkZ2dTr1590tPjaNs2CBubeCQpEzs7FU88UePe\njdyHFSuWER8fh7OzMz4+Ppw7d5ZVq37Gyqo1W7d24+hRNzIybh1ZkMmNWNkVYW1fhMa+EHePAmrV\nyickKJ/6dQqo61e2OVpubhvOnrVCq9USElIHH5+7j8AolUomTZrC6tUruXz5Imq1mvbtO9Kt271L\nKxgMRt588xOmTHkNR8da9O49liNHphMWlgVkoVSqadCgDdOnv0tpqRa12h5393qmWyoPG7FKRag0\nItbmJeJtPiLW5ZeRkcGVK+E8+WTZyo3t238jMjKM559vjp2dDpXq3r9z7eysyMsrrtD1//OfVykq\nKmLRos+oUcODxYt/ZtOmp0lNHYvRWDaioLIsxckzC0fPbHxqZxEUkI+/XxGeLhI1HCXcHCUsVPe4\nkBlIksSRIzG0bdsfO7uyzeJKS0s5ffokxcVFNGrUFFdX1we6hlilch8kSSI6+irp6RkEBQXh5PTg\nW/CGhp5h7NiRzJ49n6ef7mI6/sor/QkOrsO0abMe+BoPg+joq+Tn59OoUZN7nywIgnAb167F8t13\nXzN79gdkZ6dw5cpF/vxzB76+1hiNpTRoIKNdu6bIZBL3+sopLi4iOjoaW1sr3N29K/RLXa1WU1RU\nRFxcKtu3t+G779eiLbUFJHzqx9L0qas0bZiHk2Uql89sJj87BptCb9rV6Vcp3x8PqqREy8GDZ+nQ\noQkFBW48++yzN916sbCweCh2fq2IRzrhKCkpYdSo4Wzf/htQ9oe2cOFiBg4c8sBt+/r6sXfvLlPC\nER19lZKSkgdu92Fy4MA+nJycRcIhCP+H/v77NF988Rnp6WmEhNRlypTpuLi43HJecXExa9eu4rXX\nRgBlv7Dffns033yzjMLCAgyGLBo29CA5+Q/s7JQ0aiSnUaOuQOo/LZTvayYy8gqLF39MXl4eAC4u\nLkycOAVvb59yvyeDQYan9zhCz9dlyCvd0GvLfmk7+6TzzEt/M/SFPAK9JGJjo5k3byZabdnusVev\nXuXChTDmzfsQe3v7cl+vKhQUFHH48EXatBmMn59/tfalsj3SCcf8+XPZvv03LC0tqV07gEuXLjBh\nwhiCg0No1qzFA7Xt7x9IfHwcRUWFWFtr2LlzB888043U1BT279/D+vVrUCgUNGzYmDfeGE16ehqL\nFi1Ap9ORmZnBiBGjaNv2KZYu/YqzZ89gMBjp0KETAwcOZcyYN5g0aRo1a/qyZcsmsrOz6NatB+++\n+zYODo60atWGVq1a89lniwCws7Nn2rSZRESEs3Llj6jVKtLS0ujd+0X+/vsUUVFX6dOnP88//xKh\noWf4/vtvUCgUeHl5M3HiVHbv/pPjx49SUlJCUlIigwYNpUWLluzYsQ2VSkVISB0OHtx/Sz8FQXg8\nnTx5ghdf7I5Wq/3n8V8cPnyAXbsOYG/vQHT0VXx8fFGpVBgMei5ePE1cXH3kcgmDQUfz5jWJjd2J\nWl2Mq6sFzz33YHuUlJSU8Nlni8jLy8PT0xODwUBqaiqffbaIjz769J6TK6/GWvDBwrqcPuZDUX5f\n03Fbl1zqNPudscPTad30xkqVzZs3odXqaN68BZ07P8v69WuIiYlm9+6d9OnT76a24+KukZKSjLOz\nM7VrB1R5cTGZzIEPPliKjc2tlUIfdY90wvHnn9sBWLVqA+3bd2D8+LdYvfpn9uzZ9cAJB0CHDp04\neHA/3br14PLliwwePIzIyAiWL/+OZctWYmFhwbx5Mzl9+iQAAwYMoXHjply4EMby5d/Rtu1T7Nmz\niy++WIqzszM7dmy76/Wys7NZsaIskXnjjVeZNm0Wvr5+bNu2lVWrfqJFi5ZkZKTx449ruXz5EjNn\nTmHDhq2kpaXy3nuTeP75l1i48AO++WY5Dg4OLFv2LTt2bEOpVFJYWMgnn3xOQkI8kyePp1u3HnTr\n1gNnZxdCQuoyY8bUcvdTEIRH20cffYBWq+WFF17i5ZcHMWvWNCIiwlm16mdGjx7LjBlTGT36Nby8\nLFEqc5k9+3kgx/T6fv1aAEagciYnXrsWS25uLm5u7ixYsAiNRsWrr75KSkoK6emp1KjhedvXZeTC\ngs8C2L6+PrrSsrWo9s4FNH8ynm7PXaN101RquFgBNy+LzcrKBKBr1+6EhNQhKekpYmKiycrKMJ0j\nSRI//7yC3bt3mo61aNGS0aPHVvrqj9jYFD78cBVz507Gy6sdlpblXy77KKlQ1PR6PZMnTyYxMRGl\nUsm8efNQKBRMmTIFuVxOYGAgs2ZV/TyH62uLdTotkiRRWloKVE55W5lMRpcuXfn44wV4eHjSqFET\nJEnCYDCQk5PDpEnjkCSJ4uJiEhMTaNiwMT/99APbtm0FblTLmzlzLt988znZ2Vm0avXkba50Y86u\nh4enKZO/di2GTz750NTW9WHFWrX8kcvl2Nra4OXljUKhwNbWDq1WS3Z2NpmZmcycOQVJktBqtbRo\n0RIvL28CA4MAcHNzNw0j/tu9+ykIwqNKkiS2bv2V06dPYmNjQ2xsNABvvTWeBg0asmHDGiIiwomJ\nCScu7hgffzwYW1uAEiorqbib6xurabWllJaWIpfr0enKPkPl8ltHN9KyYd12R1YubUpWUtm8i0bN\nUpg+5TxNmmRz4yvg9l/cHh4exMbGsH79Gtq1e8r049Xd3Z3Y2BhkMhmxsTHs3r0TpVJBvXoNCA+/\nzKlTJ9i2bSvPP//SLW2WfR+UYm195+W1d+Ll5cKwYcPw83v6vpbKPmoqlHAcPHgQo9HIunXrOHbs\nGIsXL0an0zFhwgSaN2/OrFmz2LNnD507d67s/t6kZ8/nWbLkEwYPfhkvLx/i4mJRKpXlWm5UHh4e\nnpSUFPPLL+sZOfItEhMTkMlkuLm5s3jxVygUCnbs2EZgYDDLln1Dr14v0rJla/7443d27NiGXq9n\n//49zJkzH4DBg/vy9NPPYGFhQWZmBjVr+hIREYGbmxtwc6JUs6Yf06fPwc3NnfPnz5ky8n+f878L\njBwdHXFzc+fDDz/B2lrDkSOHsLa2JjU15X+SsLLXyeVyJMmITqe7bT/d3St3uZogCOYnSRITJ77N\nypUrTMeu/1ibPfs9evV6gRMn9uHtDc2by3B0zAUqd7fR9PQ0li1bSkREOFZWljz9dBdeeKGv6cvV\n17cWnp6eJCUlMXbsKIxGA1qtjoCAQFxd3dDqIDJRzolzanbtdSA51o/Ey95IRjn2DkXMnn2Onj3i\nKe9vzb59+3Px4gUiI68QGXkFKEtCDh48wMaNGwBMRbn69OlPjx69+OuvY3z55RIuXjxvSjhKS7VY\nWJTF6sqVeMaP/5Jt2z4sV6Gtdev2IkkSPXq0RafzoXPnHlV+u6a6VSjh8PPzw2AwIEkS+fn5KJVK\nzp07R/PmzQFo3749x44dq/KEY9KkqaSkJLN+/Rri4mLRaGxYvPgLGjRoVGnXePrpLuzcuQNvbx8S\nExNwcHCkc+euvPXWCAwGIx4ennTq1IWOHTvz5ZeLWblyBW5u7uTm5qBUKrGzs+c//xmGpaUlLVu2\nxt29Bi/Xn6dNAAAgAElEQVS99DKffPIh7u4eNy1p+vcf2zvvTGHevJkYDAbkcjlTpswgPT3tpr7d\n7o9z3LgJTJw4DkkyotHYMH36XFJTU/7nrLLXBQeH8PXXn+PrW+u2/RQE4dF37NgRVq5cgZWVFaNH\njyM09Ax79+5GqYSoqIOsWHEQHx8ICPCnS5dnKv36xcXFLFjwPmlpZZNI8/ML2LJlMwaDkZdfHgiU\n1ap4553JfP75YuKScpBbheDg2RmtzdMMm2BPUoId6XEu5Gf8e0KnEfcaX9O61Td06zoFmaz861jd\n3Nx5//0P2bXrT7KysqhRowaHDx8kPT0dBwcHJEkiNzcXKEuWJEkiPT0dAIWi7GtTkiRefPE9Fi6c\ngJdXfQIDn2DnziHk5SWi1aZw6tRBcnJy6devEwC//HKAS5diGT++P1qtFQ0btsXKyhV7+8ZoNJoH\nDfMjoUJ1OFJSUnjzzTcpLCwkJyeHb7/9lnHjxnHo0CEA/vrrL3799VcWLlx413Yqa+1vSkoyGRkZ\n+PnVuqUuvVg/bz4i1uYl4m0+j3Ksf/hhKVOnTsLFxZU//9xJcXEMo0e/hEYDXbs+R25uDt7ePnTt\n2h1Jkrh2LRaFQo6vby2USiWXLl0kLS0VV1c36tatd99lsq+PDHh4eDBt2iyio6NYvPhjLC0t+f77\nHzFKMmJTZJwIs+KPbb5cPF2TvAw7kG79QSWXF2Nre4Q+fRQ880w6P//8JllZWUybNoO6detXOEYx\nMdHMmDEVR0dHFi1agtFoZNy4UWRlFaNQlNUFKS4uJjISBg58jiFDBgI26PWO1KxZH5Xq1mTn6tVI\nIiJO0aRJA2QyFSUlBlQqDU5OblhbW1e4r5Xtoa/D8eOPP9KuXTvGjx9PamoqQ4YMQae7MS+gsLDQ\nVD/+bhwdrU337h5E2ZsLusfzgjmIWJuXiLf5PIqxliSJgAA/ZDJwdTWQlbWDqKhwNBpwcnLkzTdH\nms69ciWCDz74gOzsssmhnp4euLu7Exp61nROy5ZPMGXKFNOv/PIpm4vh7deAqEwPkrQuaPzHIlfb\nM/9nGy6fdyXqXE3SY9yQpLJkRiaXcHMvwtevkAD/Ivx8i6hbN44vvuiFlRW8//5alEprdu1yJysr\nC5nMgJ1dxSda2tiUzbsoLDSQmpqJv78nMpmC5GSwtlaTnl6MwaBmxIg3eOedGTg6Ot8z8XJ1bUrr\n1k0r3CdzMtffdoUSDnt7e9MsXVtbW/R6PXXr1uXkyZM88cQTHDp0iFatWt2znezsoopc/r48yr9M\nHjUi1uYl4m0+j2Ksf/zxB5KSEunbtxO9e/uRmhrLrFnvmZ7v1etFUzXP4uJi5s1737RSRKfTkpSU\nTFJSMlZWVjRu3JRz50I5ceIk69ZtpGfP5+95fUmCxAwZsflNsG/4CRGlwRxe4k5uqgP5Ge3Iy7Cj\ntPDGBEuFwshTHRLp2yeGHj1y0WoL/6c9iY0bHUlLS2XixEk4ODhw6dJllEolbm7eFa5MCuDo6Ia7\nuzthYakMGzYNg0FNUlIBAQHB7Nt3lOzsLJycnFGpVBiNkJlZeO9GHxEP/QjHK6+8wrRp0xg0aBB6\nvZ6JEydSr149pk+fjk6nw9/fn65duz5QhwVBEP7fhYdfZuHC+cTFXcPPrxZTpkwnICDwtucWFxdz\n9uzftG7dhry8TFq08MPCwh5393zmzJnBunWriIy8gpWVFU8//Qzt2j1lem1MTDS5ubl4eHjw4Yef\noNWWMnLk6xgMRnr3fpEePXpx+PBBli79moiIy7dNOLLy4FqanPh0GfFpchLSZeQXKEiN8SU5wovU\naDeMhptvPVhb6wgOzqV793h69ozDxaVspaGlpRX/lAgxkclkjBkzjoULFxAdXbbKRqlUMGLESJyd\ny18hVJIk0/y3pKQMvvlmC1OnDuWVV2bx7rsfER19DSgmODiEn35ai1qtFnPaKkmFEg5ra2s+++yz\nW46vXLnygTskCIIglFXefO65zhQUlP36DAs7y8GD+9m16wC1atU2nZefn8cXX3zGuXOhREdHsnjx\nFAICbKlZUw2UTUa0sbFh+PCRt7sMcGPVil6vR68vW5J6fXZfYmI8RqORa9diAbCwKLt1YTBCTIqM\ni7FyLsbKSUxVkZ9pS36mHQWZtpTk2JKR4IS29EaSERgYTv36V2jf3pLGjXV4exdyP1NCatXy56OP\nPuXcuVD0eh3BwXXx9Lx9jY7bycrKo2/fGWze/CE6nSUyWQ3ateuKpWUb2rRx4eDB/qZVK0FBwY/1\nEtXqIDZvEyqNiLV5iXibT3XEety4N1m7dhWdOnXm7bcn8fHH8zl8+CDDhr3OwoWLAYiPj2PAgD5E\nRobj5QWenqDRWDJz5lxq1vQt97VKS0uZMuUd0tPTsbHRoNcbbtnKQaa0QWnjT6duIyiVeXA+wpqk\nGGeyEp3JTnIiL90e6TYTPRs0yKJnzzieey4eD4/y3fZ4kM3b/k2SJMaOXcLs2a8ik9kAThQVWeDr\nW6dKJ25KkkRsbAzZ2VkEBARiZ1e95dLv5qG/pSIIgiBUraSkRACGDHmVVq1aM3jwKxw+fJDk5CTT\nOVOnTiQ2NpyOHV159dWXOXbsCGfPnmX9+jVMmjS13NeysLDgnXemsGTJIpKTkwHwD6hLcJMXOB6W\nT15xawrzA8hLt+erT+3Iz7CjpODmSZpKpRF//1yCgnIJDMwjKCiXunVz8PK6/Vw9SZKIj48jOzsb\nT0+v2+56KkkSly5dJDU1BWdnF+rXb3DPUYf9+0MJDPTCw8OF3FwDvXu/gEzWBG/vALOMWBQVFTFq\n1HBTxWZraw1LlnxF794vVvm1H3Yi4RAEQXgIBQUFc/DgfhYu/ICUlGSWL/8OAKPRCEBRUSG+vgoa\nN4Zevbrw5JPt8PT05uzZs6aaF/fD29ubOfM+4cCZYq4kWxMZa8eZLTWJv+BLQdatv1g1Gh1Nm2bg\n53eZmJgv0Gp3otHA00+/QI8eve5axEqn0/HNN19w8uQJAGSysmJcvXq9YDrHYDDw9ddfcOLEcdOx\nhg0bMX78pFuWoZaW6rCwUGE0GomOTsPKyhsnpzp4e/vi62ve2yJz5kxnx45tWFpa4u3tw9Wrkbz5\n5giCgkKoU6euWfvysBEJhyAIwkNo/Ph32bXrTy5fvsTUqRMBqF27FnZ2lsTGHkClyqNhQycuXIBD\nh/bj7l6DI0fKaiHdbrTgbjJy4fB5BccvqYm/4kfceV/SomtgNJZNsHByKqFZswxCQnIJCckhODiX\nmjULSU9P5L33pmAwlKJUQkkJrF+/BktLS7p0efaO19uyZRMnT55ArVbj51eLyMgINmxYh6+vn2n3\n6r17d3HixHEsLa1o2rQZYWGhhIWd4/fft/Diizc2aNu06SBhYVcZN240FhbejB79zG3rYpjLH3+U\njWz88svvtGjxBCNGDOO33zaze/dOkXBUdwcEQRCEW7m4uLBr1wE6dWpL/foBNGniSa9e7XF3dwSK\nARUdO3bi4MF9JCYm8vnnnwKgVqvp27f/PduXJIhKlnHgrIIzlyyJu+BL/Hk/CnPKJprK5UY6dUqi\nb99oOnZMRqW6dbrf8ePHKC0tpXHjJrz99kQOHtzPihXL2Lt3N5aWlhQXF+HvH4C//80ra86dCwXg\nrbfepmnTZqxbt5pt237j3LlQU8IREREOQP/+A+nc+RlTAbHLl8M5fPgc7do1Ij9fR5s2XenePQhX\n1/JPHq1K10d2rk+PvP7vx71seXmIhEMQhP8Ler2epKRE7O3tsbd3qO7u3JHBYCAtLRUPD08MhkI2\nbPgYFxfptkUSraysmTFjLlu3/kpCQjz29g4891yPu04YNRohLEbOvlAF5y/aE3UqkJRITwz6sva9\nvQt4+eUYXnopBnf3kju2A5gmlnp7+6BUKk3XTU5OZOnSr03n9e378k1zGK7PpSgoyEeSJAoLy+pa\n/HujNguLshodsbExGAwGoqOjAFCp1Kxbt4+goPZ4eDTE0/Pm6tLVrXv3nvzww3f06dMTT08vYmKi\nUalUdx3x+X8hEg5BEB57hw8fZNSo4aa5DQMGDObjjz9Dra7cTcoqQpIkMjIyTLdB9u/fw9atv/LO\nOwOwscmjRo273x6wsbFh0KCh97yOwQgnLsvZG6ogIcmKK8dCiL/giyTJkMkkOnRIYvDgKNq3T0Gh\nuPPixXPnznL8+FEMBj12dmWJ259/bictLdW0pNRgMOLj44OPjy9//XWMjRvXU69efQICyipCt27d\nhqioKL777hvWr19Lbm4OcrmMli1bm67ToUMnjhw5xIED+9i+fR8yGSiV0KTJC0ye/Co2NveuZl0d\nZs6cR0pKCtu3/0ZMTDQajQ2ff/41ISF1qrtr1U4kHIIgPNauXYtlyJD+FBUV4uTkRHZ2NmvXrsLe\n3oG5c+ebvT/Xa1pcr6URERHOpElvs3TpZxgM2fj765k5szuWlmW3TSpDZKKMTYeVxKeoiP3bn+jT\ngWhLlSgURgYMiOL11yOoWfPe1TP379/LDz98d9OxWrVqExMTfdMEUEmCd999D0dHR+RyOUeOHOLq\n1UhTwvHMM93Iycnhjz9+Jzc3B0tLK1577XUCA8ue37HjL5KTsxg7djw//riCxMQs5HINY8fO5cUX\nR1RKTKqKlZUVK1asIi7uGtnZWfj7B2Bj8+iVxa8KIuEQBOGxtnPnHxQVFdKx49OsXbuJY8eO8OKL\nPdi0aUOFE45/V6uEsv2jrK2tkclkSJJEePhl0wRBvV7PkiWf8M47k4GyQl2vvTaYVau+x2jMwcIi\nh2XLRmFtHftPa/J//nlwWfnw2zElZ64oSbjgQ8zpEPJzypazdu6cyOTJYdSuXb4aDFqtllWrfgSg\nR49e2NnZsX79GmJiohk7djzFxcXY29vz3XffkJeXR3j4JerVa0Bc3DVKS+HSpRSuF6BOSEjHwsKH\nb7/94Z9zkzl6NJxGjUopKbHE07MB1tZ6goLq8/nnY1Ao1LdszPmwq1nT975qofw/EAmHIAiPNZ2u\nbPOw67+2nZ1dANDrdXd72U2uXbuGpaUDcrkcg8FAo0Yh7Nt3GLVahUwmo3Hjepw5E4ZcLkerLWb4\n8CGsX78ahUKG0VhKbm4isbEHkcv1GI1FbNw4GQuLuH9alwMW5e5LVlYmYWHnMBgM1KtXjxo1bkyW\n1OnL9i+JS5MRlyYn9Kqc+MueRJ2oQ0562a/s+vWzmDr1HK1apZf7mgC5uTmUlmqxs7Olf/9BAJw9\nG8qlSxdRqVQ88UTZ/lldujzL2rUbmTLlc7y9wWAAudweg8GJ3FwXZDJL8vL0FBZao9W6o1Y74Otb\nA1fXJ7CxaUWNGraiqN1jSiQcgiA8EEmS2L9/D1evRuLl5cMzz3St1mWJ/+uppzqiUCj49ddfiIuL\nIyoqEoBOnbrc8TVnz/5NYGAwGk3Zio0+ffrw0UfzcHS0wGDIY/36ecBf6PVloxynTn2L0XgMoxFU\nKgU7dsxDJosxtTdpUjfgegGsin/shodf5tNPF1JUVNaWSqXkjZFjUDo+yeEwBfHpMgxGGcX5lmQm\nuBB3thZZyU4A1KyZz4QJF+jePf6+yolfZ2dnj0qlJC8vn99/34q9vT2RkREAODg4snbtHvr3f5ru\n3XuTnq7jp5/+5PJlHUFBTfn8829Mt1MAvLwa0Lp1zwrHQXg0idLmQqURsTavhyHeBoOB//znVX7/\nfYvpWMuWrVm37lfTl/XDYO3aVUycOA6drmxUo2XL1qxcuQ4HB0egbNRAqVRiZ2eP0WhkzJj/0L//\n8/j6OmA05uLmpkSvN1Tr0ka9Xs/bb48mJyeH4OAQLK2suZxgh8anH5LSm9SrHhSnuZB6zZms9Btz\nBlSqFIYMOcakSUbUauMD9WH79t9Zu3YVALm5YGMDHTq0Z/Dg4Xz00QZGjRqJg4MfLi7eyOXyW249\nldfD8Lf9/0KUNhcE4ZHw00/L+f33Ldja2tGr1/Ps3r2TEyeOs2jRh8yaNa+6u2cyYMBgOnZ8mkuX\nLuDg4IibmztFRUWmhGPhwvk0blyPNm3qYDBkMndub9RqBZAPyNFoLCplb4//VVhYQH5+Ac7Ozvcc\nFcrISCcnJwc7p5q07TmXw+eVWMrVRJ2tSfw5f4oKbuwNolDk4eBwHB+fv1CrPyYqSktW1iJq1Kj4\nrqeSJNGlSzccHR04evQw+/bF07x5B15++W2srGqyePGtO8iK2hPCv4mEQxCECjt9+iQAU6dOZ/jw\nkfz55x8MHdrfdNwcLl26yIwZU4mMjMDT05N3332PTp0633KejY0NnTp1wWg08vXXS1CpjHTv3gGj\nMZ8332yPvb0SuTwHqPpS2AaDgdWrf2L37p1IEmg0GkaMGEnz5k/c9vzMPDgT44xd/Q9R2Yawbqcl\n187WIibUB72ubBJovXrZ9Ox5jZ07R2JhcZpPP/0UN7cQPv44gLCwc1y+fOGBEo4FC1bi7V2DHj16\n0qDBEKZO9cPKyureLxSEf4iEQxCECrO1LRs6PXXqBAMHDuWvv47ddLyqxcbG0KtXV/LycgFISUlm\n4MA+/PLLb7Rt29503r59u1mzZiUzZ47BaEyhX79aWFlZAGn/nGHeehy//76FXbt2AmBnZ0deXh5f\nfvkZc+cuMK1skCS4miRj798KLscpKMp1JCWyEUkRDuSkenF9JYu7+3EWLCihfft0QOL48VC0WgOl\npaUYjUZKS0uBm4tqlcemTQe5ciWekSP7IUk2vPbaWGrUqIOt7cNbNE14uImEQxCEChs8+BVWr/6Z\nzZs3sXnzJtPx114zT62En35aTl5eLm3bPkVCQhzNmjVi06YtfPXVEjZsWMvChYvR60uoVcuamTNf\nwNExg7KPver96Dt+/CgAY8eOp0WLlnz99RKOHz/OyZMn8Pb25Vx0WSXQa6lyksK9iD8XQEbSjS96\nmawEZ+ftdOiwjffea4Od3fUiWDJatmzF4cOHmDlzKhqNDbm5uVhZWdGgQcO79ik6OolffjnAW2/1\noaTEjlatetCxoxuenrVRKsVXhfDgxF+RIAgV1qBBI9au3cSUKe8QFXUVLy9vpk+fTefO5injnJmZ\nAUCnTu3p3r0eSUnnuXZtC5aWV2jQoB4pKcdQq/NxcVEDlmbpU3kYDAYAFGoHwqIVFKifQuPfnLOZ\nDTizXE1RqYyMOGdiT9QjJb5snolGo6NDh2SefTaB9u2TsbHRI5N1u6XtoUNfo6iomDNnTpGbm4u9\nvQNvvTUOJyfnW84tLdViYVE2umNpaUO9ek9gZdUWd/dbzxWEByUSDkEQHki7dk9x9OjpCq9IqKj4\n+DjCws6i0cCBA0vQaJ7mxIm/UKmgYUN/+vdvDZRi7tsl95JbCB51BlHkZM9PR+ohkymAVqidZaRn\nW6DUWRN1OojIC+4A1KhRxLhxF+nd+xoWFv9eZXL7WFtZWTF+/ESys7MpLi7G1dX1thNSJUliwIA5\nzJnzFjVrNsfXtyPBwQ9PUiY8fkTCIQhCpTBHsqHVaoGyHVHd3NyZMOFVjh5dQXj4RbZs2QyAh4cH\n/foNqPK+lFdOAcSmyIlNlRGdLCcuVYZRakdBkQOJJ63ITHCitMgBndYJSbpRIMPGRsvIkeEMGxaJ\nlZXhvq/r6OiIo6PjLccNBiMKhZzCwrIKqP7+T2JhUf7CY4JQUSLhEAThkTFy5Gv07t2Lpk1rIUlp\ndOzoQ/v273H06GGSk5NwdHSibdv2WFtb37uxByBJZVU9S7RQopNRooXiUsgtlJGdLyO7oOyflCwZ\nOQVliZjRICMvzZ7iZA9iL3qRmX5znRKZTMLVtRh392Jatkxn5MjLODlpK7Xff/99haVLtzJ//nQc\nHJrj4eFaqe0Lwt2Iwl9CpRGxNq/7ibckSRw8uJ/z58NwcXGhR49e2No+nLtt/tuePTvJzMzk2Wef\nIT8/htLSZBwdZajV5vmtJEkSsbEx5BUUUGioQZHRg5hkGbGpcgqK7zyiYzTIKMqzRpdrhy7HgcxE\nJ+JjHNBqb/Tb3b2Irl0T6Nw5CT+/fFxdS1CpqvbjOCNDR3GxO3XqPPFQ18gQnyXm80gU/vruu+/Y\nt28fOp2OgQMH0qJFC6ZMmYJcLicwMJBZs2ZVuLOCIFQeSZKYPHkCP/74g+nYZ58tYsuWP/Dw8LzL\nK81Pr9eTkBCPn18tjEYjGo0Svb4QrfYoDg5qKmv31HuRJEjJMrJs3V8k5zojt25KQbYTuan25KY6\nUJBhj6RXYaXRYq3RYWOjxcraQEGuJenJtqSlWmMw3Fo/3NIyAgeHXbRufY65c9tjbV31cybef/8n\nWrasR9OmnfHwaIylpaidIVSPCiUcJ0+eJDQ0lHXr1lFUVMTy5ctZsGABEyZMoHnz5syaNYs9e/bQ\nufOtxXcEQTCvnTt38OOPP2BlZUW/fgM5fvwIV65EMG3au6xYsapa+6bT6Thz5jStWrUGyiaCvv32\nm3z11Ufo9YnUqSOhUHibpS+5hXDyfCbnrxaTVliD9DR30mJeJy3WnawEJ4yG8n9cymQSXl6F+Pvn\nERAQx6lT72NhcRQPDwW5uTmkpOhZsyaa4cNHPnC/4+JSuXo1gU6dmmE0Gtm9+xTnz8cwcWJ/jEYj\nPXp0JSDgKTw8aj3wtQThQVQo4Thy5AhBQUG8+eabFBYWMmnSJDZu3Ejz5s0BaN++PceOHRMJhyA8\nBEJDTwPw6qsjmD37fcLDL9O+fUtCQ8+YvS+lpaUMGPASv/zyG1qtlpycDD76aA5LlrwPFCCX57F8\n+WjU6iTKVmFU3bB/XiGEx8u5kiAnOllOaoaSzPiapMW6kxbjSnHezduhu7qmYDTup169TEaMqI+r\nawm5uSry89Xk5qooKFDh5laMn18Bvr4FphUlu3fv5PLlbTRs2JBJk6YRGxvDjBlTOXnyRIUSDp1O\nT2RkAnXr+iFJEpmZBVy9mkvz5j4olbY0auRDQEARhYWugIwnn+yKQlH11VMF4V4qlHBkZ2eTlJTE\n0qVLiY+PZ9SoURiNN5ZraTQa8vPF/TdBeBjY25etVDh16gSpqans2bPrn+P2Zrn+hAljmDJlBjKZ\nltLSZF555RliY3/DwkLC0lLJzz+PA65vlV61t0xiU2Scj5ETHifnWoqSvDR7MhNcyIxzJSPBEcl4\n4yNRqczAweFPnn02m+nTg9i48Qt27NhOvXpNSUz0JSFBolGjJjRvHnLXa16fJqdQKP9pV3nT8f9V\nXFyMQqFArb6xnFevN6BUliUNubkFLFiwmiVL5qNSudGiRUfat79xm8TdvQKBEQQzqNCk0U8++QRn\nZ2eGDRsGQO/evYmLiyM0NBSAvXv3cvz4caZPn37Xdv79P5EgCFUjKSmJRo0akZGRcdPx77//nuHD\nh1f69fbt24efnx9+fn6kpyewc+dm6td3x8PDApWqehbGJWdKrN0tsXOnB5nxLuSmOpKXbo/xpnkW\nBtzcLvDKK5Z06JDMggXPYjCUrRLRaDQUFhYCIJfLMBrLPjZlMhg16k26du16x2snJSUyZswYdDo9\nPj7epKdnUFJSQocOHZgwYYLpvPT0NBYtWsTly+HIZPDkk20YM2YMlpaWtG79JmvWzMPFpTZqtTvu\n7v43JSSC8Cio0P/9zZo1Y+XKlQwbNozU1FSKi4tp1aoVJ0+e5IknnuDQoUO0atXqnu1kZxdV5PL3\nRcx2Nh8Ra/Mqb7xVKlt+/XU7EyeO48KFMJydXRgzZjy9e79MSkoOpaWlD7yMVKvVolarkSSJw4eP\nEB9/kcxMV+zsjLRt6wFAcbGO4mLdfbVrMBgICztLZmYmNWp44O8fwJo1K//ZHE6iadPmDB487I6b\niBWWwJ+nFPyx35GzOxuTk+Jkek4mkwgKyqFx4yxq1TrF3r2jcHWVeOGFOWRmZmAwaJHJQKVSU1hY\niKWlJQaDHp1OT9u27VGr1ezbt4fvvltKvXqNTSNGBoOBzZs3snfvHkpLSwgMDGHw4GGsXbuK+PgE\nAOrVa8DAgUNNO9AWFhYxc+YMUlOTkctlpKRI7Np1lNzcYkaMmMHPP6/C27uhKcnIzS2lrKjZ40l8\nlpjPQ79KpUOHDpw+fZo+ffogSRKzZ8/Gy8uL6dOno9Pp8Pf3v2vGLwiCeYWE1GHbtl2mxwaDgXnz\nZvH9999QUlJCcHAIX365lEaNmtx325s3/8LhwweYMGE4en0yL78c8K9lq7eu1Civ0tJSFi36kMuX\nL5mOOTk5kZWVZXp88OABcnNzePU/UzgVoaSoBPQGGToD6A1wIUpJ2KEQok4HIBnluLsXMXhwFI0b\nZ2JldYKDB38lPz8PCwtfate2JT4+jnfeGWdqv0OHpxk69FVAx3/ZO++4qso/jr/PXYzLvWwQBAQE\nBQQVxL1XapqWK9McaWp7D0vLbGiWo2W2zNLUtpp75FZUnCAqIkM2yL6syx3n98fN688UJ2jWeb9e\nvornPus8B8753Of5jqysPKZNm4KbmxuTJz+BIAikpJwlLS2NnJwsq+D49defWLNmtbWPhIR4Cgry\nmT17HkVFBdjZqfHy8qKiohqw7PS++uoCEhJyCAx0Zvr0mWzbdoQFC74mPv4IL70UibPzRaEkIXG3\nctP7my+99NJlZUuXLr2lyUhI3EkqKir48MNZHDy4H41GwyOPTKRv33vv9LTqhXnzPuDTT+cDIJPJ\nSEw8zfDh97Nr1wE8Pa+ewry6upqtWzczYMBAiovPEx7uRmRkDxwdc/6qUTfHJitX/sqpUyfRajU0\nbx7JwYP7KSoqQi6X8+abM5DJZLz3wWcklbZj2heOpCf4YTL8dUQrWI48sk/7UFHigCCIjB6dxIsv\nxqPRGElIOMEHH7yHyWSxPUtMPI23d0Natozk1KmTKJUKOnXqyoMPjkSpVKLVaqmosByvFBUVcuTI\nYVQqJdnZ2QBotZbEaqIosnnzRgCef/5lAgMDee+9GeTm5pKcnESbNpad359/3sbp0xk89dQT2Nr6\nMfPkYooAACAASURBVHz4czzyyMM4OTUiLGwkjRsPYd68FYhiuTXbq4TE3Y4UaVRCAot75ogRgzlw\nIMZatn37nyxY8BXDho24gzOrH3744XsAFi9eRs+evRk6dCAHD+5n3bo118z0Kooia9f+TtOmNjg6\nVtOggQqo+9gOSUmJAEyYMJlWrVqj1WpZv34tSqUCR7dg/jwixyZ4AcmHQ0g5HITJcOXHWZMmJcyc\neYjIyIs7I7/8sgKTyUyXLl1p3botS5d+T3Z2Fn369OOll6ZcsR8nJye6d+/J9u1/Mn/+h9by6Og2\neHlZjo2MRqNVIISFhWFnZ4+Pjx/nzuXy3Xd/EhXVmtJSgR49hjF8eFM0GsuuSNu2DqjVDsTFHWPg\nwL6UlpZQUVFO48ZBeHhIVqAS/w4kwSEhAWzdupkDB2Lw8PBk/vxPOXBgP598Mo/33pvxrxQclZUW\n+6nQ0DBsbW1p3DiIgwf3W8v/zvfff0toaDMaNXKnqiqF994bgkplpj4To10IUJWcfJYWLSLJLyzH\nxqM3KveuzFhiS8YJPxL3NqGmymJ/0qtXFuHhxYgiiKKAKIKXVyWDB6ddFsHzggHtwIEP0KCBF2fP\nJrF69Upr9tnaGDduAk5OThw8uB+z2UxkZCuGDn3QGrVTqVQSEBBISkoKo0e/SPfuoRw6dJDKSggP\nb4fBEIGfny8y2aVHTW5ubnz55SImTXrEKnobNPDim2+WXFZXQuJuRRIcEhJAdrbFmK9nz9707t2X\nzp278ckn88jJycZsNv/rHvrt23dk48Z1DBs2iGbNwtm0aQOANQDX77//glKpom/feykpycPBQUQU\nk1CpzmNvL+d2PDp69OjN8fiTbNyTy7bTx1A5P4Ho5EpKfEMyT3pQpbO4+zZufI6ZM5OJjr66WPh/\nGjTwoqSkhGXLlhAV1YqdO7dby6+GXC5nyJDhDBkynP37E4iMDLYacj7xxFxmzpzEmDGPMmPGbFJS\nisjM3EtZGUyZ8i5PPvnMVfu+555+xMQcITb2AAqFkk6dOqPV3h7XZQmJ24EkOCQkgMaNgwFYs2Y1\nYWHNOHw4FoDAwMb/OrEB8OGHH5Gamkxi4mkyMtIB6NChE/7+vmRnn8bJyYRcXkROzjo0GgX33HN7\nQqAbDAYyMnPIKtaSWtgOz/Zt0ZU6kHHKh+zTnpQVuFnrurrm8/TTsTz8cAU3mhZk5MjRvPfeDI4e\nPcLRo0cAaNo0hPbtO9baZu3afbRvH46rq5bq6hqWL9+Oi4sfHh6uCIINI0aMw2QKo0mTvixaNILN\nmzdQUVFO69Ztadky6rrm5eXlzcCBD9zYxUhI3CVIydsk6oy7ea1FUeSxx8azcuVv1jKVSsXy5b/S\npUu3Ozexq3Ar67116yaWLVtK//7dKC/PwsfHhcDABjg52d72WBk1RkjPEzgQX8yBuEKwDcRkUpN7\n1pvcRB/y0twRRYui0Gpr6Ns3k/vuS6dt2/PI5Tf/+MrOzmLLlo3odDoaNfKnT597L4ltkZ9fjEwm\n4ObmhEqlYM6cVXTv3gN//zDUak+0Wud/dAK0u5m7+Vlyt3E73WIlwSFRZ9zta202m1m69DtiYw+g\n0WgYOXIMERHN7/S0auVG13v//n1ER7chPz+JysosBKEYNzf7276Do6+B1FyB5GwZydky0vIETGYB\nfaWK82ke5CdryEtphMlkA4BKZaJXr2wGDTpH58651pDhdY1OV0llZTWenhYX1I8//gUvLx/uuWcQ\nAQHhGI029TKuxOXc7c+Su4l/fBwOCYl/IzKZjLFjxzN27Pg7PZU6x2w28803XyAIWYSE2KPRyACH\na7a7tTHhfKlAbrFAfrFAXrFATpFAdoGAWRSoqVJSmudETZEr2UnO5GV5XNJeo9mDu/sSli3rjodH\n3UckLi7WkZV1nvDwQERRZN26GDIyipgw4SFkMg2TJ8/A1bUhMpkMZ2fpBSghcatIgkNC4i4jMfE0\nO3b8iVZrT7t2XQkICLxmm4qKEqZNewhX1/qdW5EOTqfLSEgVOZMpUGNSIopQpbOjvMiByiINNSVO\nFOU4U5B3qeARhCoaNIjlkUccaNMmiQ8/fAizWUSj6Uhdud3+vwFwUlImW7fG4ePTBoXCmQcf7ICD\ng/ZfabMjIfFPQBIcEhJ3EatW/caTT07CYLCECLezs2Px4h/o0aP3ZXVramp4/PFHef31V7CxScPV\ntX7yFokixKXIWH9ATm6xjJoqJQUZ7pxPcaA4R0NlmTtm0+XHETY2Rvz80vDxSaJ9+0xiYl6jurqQ\nhARPYmPLMZtFmjVrVmvY8htFp6tk1KgZLF78NjKZBxERQ+je/TlJYEhI3CYkGw6JOkNa6/qluLiI\nli1Dqaqqon//gZjNBjZs2ICzszPHjydia2sLWITGxx/P5ccfl2MyFTN4cASTJo1Hra77I5RiHfy6\nS8GR07akHQlEl+1BTroTfw9prlTmEhFhomnTasLCivH3T+Hnn5+noCDbWqdt2/acPHnCmmk6MDCQ\n559/BWdn5zqZa3W1gTNnjERH98fG5sbsMa71u52cnMSZM2do0KABLVtGScakt4j0LLl9SDYcEhIS\nl3H69Cmqqqrw9w/k22+X4ubmgKdnAwoKzpOWlkpISCgLF37GmjUrSUmJxcsLNBo4fHgvc+acZ9q0\nt5DL62aXw2yG3fFy1h2QU5Cr4fi6tpQUWASNINSg1e5j3DgtHToU88cf00hJ2c/TTz9H27aWOB/z\n5s2moCCbhg0bEhzchD17dnPgQAzPP/8Sjo6OqFQ2+PhcHiDrRtm/P4GTJ9O4//7+qNWRdOzoe8vX\n/ndmz36PuXNnW3/u06cfX3/9vVUASkhIWJD2EiX+k4iiyBdffEbr1s0JCwtkzJgR5ORkX7vhbcZs\nNluTl10IAnX+fDq7d69g7dpPsLEpxNkZamoyyMo6Qbt2gRiNsTRtCm+++SLvvjsLjUZDUtIZa6jw\nm8VkgqRMgVV75by3XMnvexTkpXhy4OculBQ4EBZWzBdfbKNNG2fCwrozaNBWwsLSqahIBbjE5fTs\n2SQAnn/+JR599DG6d+8JQHr6OYKCmuDn1+imxEZRURk//LAJsNxjNzdngoKi8fbuhZtb3YuNDRvW\nMXfubGQyGZ07d8PBQcOmTRuYM+f9Oh9LQuJuR9rhkPhP8tlnH/POO29af964cT1JSWfYsmUXDg4X\njx4SEk6wZ89OJk9+st7nVF5ezvvvv8P+/TGo1WrGjZtAhw6dmDr1ZebPfx17+0LGjg3jzJmTzJ37\nGAB+fpajiIAAE5CByVSCnR3Y2toSGdkKuVyOn18jEhJOoNOV39B8RBHyigXOZgskZco4nSGjuuav\nEN5yEX1aY/aubIYoCtx7bwazZx/E3t5EcnI0e/bsYurUVxEESz8eHh7k5+ezaNFXaDQaVCobQEdc\nXBz29g6kpiYDYG9vf0NzrKioZu7cH3njjbFUVtag08kpLITSUk8UCi3BwZ2IjNTeUJ83ws6d2wB4\n9tkXeO21N9myZSOjRg1nx45tTJv2Vr2NKyFxNyIJDon/JJ9//jEA8+Z9Svv2HRg1ajjJyWf5/fef\nWbv2D+bNm40g6KmsPIdWqyctbTN2dkF4el7bI+RmMBgMPPTQkEuSxx0+vJc333ychQsnYW9fBih5\n6aWXWbJkMceOHUYQZLRp046RI0db2zg6OlpTuM+b9wEuLq4kJJxAJhPw9fW76hz0NZBZIJBxXkZq\njsDZbBnlVRdtEZzV4KXSUJxhw6EYJ06fbgbA88/H8+STp6zRPh955FFkMhl79uzCbDbTuHEQtra2\nLF36nbUvlUoJwJIli1myZDEAWq2Wtm07XHWOoiiyYsVWhgzpavF+qVLRqFEIpaWNcHb2pUEDLeHh\nQ4mPP87PP39DZWUVHTp0ZPDgYfViV6FUWnZtcnNzEUWR3Nzcv8qVdT6WhMTdjmQ0KlFn3C1rbTab\n8fZ2wWw2c/ZsBlqtI088MZFff/2JmTM/xN/flagom8sibqanFzJt2lJ++WUdCkXdavUL34zd3NwZ\nPPg+lMoiYmNX4eSkYeHCb67YRqu1o6ys6rLykydPMHfuB5ekNX/44TH07dv/knqiCCfSZBxNkpFZ\nYImVIWJ5KZtNAlQ5YFOjQV+qJTfDmVMJLlRVXbxumayc0NBHePvtxlcM3W00GjGZTBw7doRPP/0I\ntVrNwIEPEB9/nBMn4vHw8MBoNFJerqNRowAmTJiMj4/PZf0YDEaMRhN2djaYzWbmz1/JoEHDCQiI\nwtHR9TIh8eefmxkz5iGrJw9YRNDs2fOuuI7XQ22/20eOHKJ//96YTCacnJwoKSkBYPbseTzyyKM3\nPd5/nbvlWfJvQDIalZCoR2QyGc2bt+DYsaOMGDGEiIjmrFplCWnu6+tMy5aqK4b39vV14fnn7yM/\n/yze3iF1Np+9e3fz9ddfYGsL/fo15YUXeiIIMG7cKnQ6HQaD4Ya+MYeFhfPeex8QG3sAo9FAs2YR\nNGnS1Pq5+S831s2xctKyVRRmumEot0OoscdQYY+u2J6cTA0Gw+UGpnZ2STg47CEsLANf3z0kJ29h\n4UI1H3/8+WVGkgqFAoVCQWZmBgBdu3anf//7aNu2Hc899xRVVZUsXLjosjFqaoxUV+vRatUAzJix\nmGbNGtOnT39sbHx4440vazV+FUWRF198FoPBwODBwwgJCWXu3NksXvwNDz44kqio6Otex+shKiqa\nr75azIsvPkNJSQlKpZJnnnmBceMm1Ok4EhL/BiTBIfGf5KOPPmfIkAEcOnSQQ4cOAjB06CBatrTD\nxubKfxaCIBAZ2Zjq6jNkZZlRKFzw9GxwU+NXVlZib2+PKIq4u2vp0yeCqqptpKcfZMMGH86dSwPA\nw8PzprbnGzRowH33DbqkTBTh2FkZGw/JySmUkXPGm8TdEVSUXdmbws9PR2hoKU2bFhEQkIuX1wk+\n/vhFvLy8+OCD+UBTXnjhGOfPnyc3Nwd//4Ar9uPiYgkVfvToYaKiWnHokCUxnpOTpdxoNFFWVoGL\ni8XWYtGiNZhMAiNHDkMud2bKlA9xdfW9xOi0NioqysnOzkIul/P5518jk8mIj49jzZpVnDmTWOeC\nA+C+++6nX78B5Ofn4ezsUmdxQyQk/m1IgkPiP0lYWDN27Ijhjz9WotPpcHS0p3dvH+ztr/0nYWur\nJCNjL6++upjVq7de0dDx8OFYPDw8rXYT+/btITCwMQ0aeGE0GunZsxP9+/dk//5fsLOrIDg4mPbt\n23PgQAw///wjAEqlgvHj62ZbvqIaVmxTEJ8qp6rUjnP7W3A2wROAZs2KaNmyCB+fCnx9K2jYsIKA\nAB0ODgY2bFjLb7/9TGxsDWq1ZcehqKiIjIx0zGYzpaWWI4SrGXu2a9eRDRvWkZ2dzbvvzrCWDx48\nDID162M4ePA0L7zwGHK5M2PGTMHFxee6BMbfUasd0Gi06HRlfP75p4SGhrJ37y7Akom1vlAoFHh7\nN6y3/iUk/g1INhwSdcbdttZVVVUolUqqqsooK4tBq70xN0xRFCktNQKenDlTRE2NmX79+iOKIp98\n8iEBAR60aOGLTKbiu+9W0rFjJ8LCWmA2G/jxxzns2PHzJWnV3dzcGDjwflJTU7G1taVz5674+TWq\ndfzabDj+ztksgaVblRSXCZScCeLQn03R6xU4OuqZMiWOoUNTuZIH6u7dO/nyy88Bi0trTU3NFfuP\njIzixRdfrXV8s9nM/v3xJCTsIzU1FVBx9qyZRYumYzA4oVB44uLiX2c7A9988wWvv/7KJWWdO3fl\n559X3XQckrvtd/tuR1rv24dkwyEhcRuYOfNt/P196NvX/4bFBliOWJyclEARgpBIebmMjIz9mEyF\njBoVgo2NEqgBanjuOUvocVG0xMI4fnwzggATJz5GeHgEc+a8T0ZGBgqFkgkTJtXJ9ZnMsClWzubD\ncgzVSs7tjuZUnGVX4/7703jtteO4uelrbb9r1w4ARo4cTb9+/fn66y/YtWsHXl5elJSUIAgC0dFt\nGD163CXtamqM9O79HKtXz8ZksqWmxsTcub/z1VfTARGzWcG5c1U4Ovaql+injz76GI6OTixbtoTK\nygrat+/Eq69OrbOgZxISEjeHJDgk/rOMGzcce/t0HBxuPf5d69YXjDKLscTTu3KfFzwqKisrAGjR\nIhInJycCAhqTkZFBRcWNxcq4EiYzHD0rY+thOTlFMmTVauLXtCUrQ4OLSzVz5x6gS5e8K7bNzs4m\nKSkRGxsb9PpqADw9GyAIAp6eFrHStGkojz46+ZJ2c+b8SM+erfD3b45MpmXp0uU4OYXi4KBBEAQ2\nb77UniSwfryLrQwbNoJhw0bU7yASEhI3xC0JjsLCQoYMGcLixYuRy+VMmTIFmUxGcHAw06dPr6s5\nSkjUGX/8sZKWLaMwm7Nxciq4KTuBuiAwMIjTp08xZ84s/P0DrXYGgYFB12xbrIPtx+XUGM14OMrx\ncTfT0E1EpYADp2VsO6qgsExAJog41biz8vvWVFQoCQ0t5ssv99KwYeUV+92xYxuLF3+NyWQGsHqd\nfPnlApo2DSE+/jhgsX/5O9HRzfD27oqfXwsAruDdKiEh8R/npgWH0Whk+vTp1ofSrFmzeOGFF4iO\njmb69Ols3bqVXr161dlEJf7ZlJaW4OqqviT99z8BvV5PWVkZ7u7uiKJIVtY51OpSoqPdkcnu3Abf\nhAmTmDnzbdLS0khLSwNg0KAHaNq0dnfbsgrYckTO3hNyTOYLxh8Xr0GpEDEYBRRykVaBAueOBPP9\n4tDLIoFeidzcHKvYiIqK4vz582RkZGBnZ0dlZSVHjx4BoG/fe2nfviMAOTmFNGjgQnGxQJcu43F0\ndLv1hZGQkPjXctNP3NmzZ/PQQw/x5ZdfIooiJ0+eJDra4nLWpUsX9u3bJwmO/wiJiad56aVnWLz4\nTQoKwNMzDGdnlzs9LQBWrFhKZuY5xo0biMlUwPDhjVGp7vxJopeXNzNnfsChQwepqKgkKCiIkJCw\nK9YtLIO9J+TsipdjMAo42oN9WUO0Nna4++aiciohp0hGYRkEeQqkHw9i/utNKCuz7N78PRLolTh7\nNgmTyUyLFi154YVXqaqqYtKkcVRVVTFlylTKy8uRy+0JCPC3HgtNmfIFTz45ng4dRmBnd2MhySUk\nJP573NST9/fff8fV1ZWOHTvyxRdfABZL9Auo1WprimmJu5cNG9bRp0+/a+5YeHt7MHfuU7i7G1m+\n/DdKS9W88cadSV6l05Xx00/LeeSRieTlJdG1qzdarScqVcFfNW6P2Dh+/Bj79u3BbDYSHt6CLl26\nXRYRU6PR0r375aL8Qg6T4yky4pJlZBZY1l9rb8bfxo2VS1uQlnbBCrwZGk0Nbduex8engvfn+FuF\nRtu2+Tz//Alaty64bIy/c2GnMj8/j4qKcrKysigvtwRJCw1thskk8vnnqzl7Vkffvt0RBFvmz/8I\nH59IyRhTQkLiurhpwSEIAnv37iUxMZFXX32V4uJi6+cVFRVotddOmOTsbI9CUf8Pq9pcdCQuIooi\nkyZNYtasWRgMOior81m+/Gu6dWuGTKaiutpA//5DOHnyJHK5nJKSEjZt2kSPHh2AozRtarnfTz/9\nABUVRmSyclxdvW5qLkuWLEGtVjNkyJDrqn/gwAEiIiKwt7fHwUFBQUEaxcU78PISkMlu/zfvzZs3\n8dlnC6w/x8TEkJubycSJV/c+MRhFYhJgc6xI1nlLmYAZY+lRClPT2R83koKCTgAEBZXTpk0R+/a5\nkpamZuvWizEgOnYs4OWXz9CxY+FfJdd2N+3QoQ1ff+3KqVM5TJ48AVEEnQ5CQiJRqUJQqxswd+6w\n/3SOEOk5cnuR1vv2cbvW+pbjcIwZM4YZM2bwwQcfMH78eFq3bs306dNp164d/fr1u2pbKQ7HncVo\nNKLX67G3t6eoKI+tW38lPNwTb2+7y4Sg2WwmJSUXZ2cvZDItBQV6Vq5cy5NP9kCttujW/48LUVJi\nRqWKwtPz2qLDaDSSlHSG0NAwKisrOH58OzpdHs2atcPZ2Z+4uDgCAxvj4eEBwDvvTGfUqDH4+wdQ\nUaHjxRef4tFHR+Pt7YzJVICTU93akRhNUGMEe5tr1zUYDDz22AT0ej0DBgzE0dGRH39chslkZvbs\nuTRseLk1ZaX+ryOTODlllQIymUhwAxm5iSXs3JiFrqw1ZWVdEUUlcnkpEyfG8txzxbi62lJWVkVW\nlj0xMR4kJWnp3j2Hdu3OX9d1iaLI6dPnCA31p7q6hs2bD7No0Y/k5uZSU6Ni8ODRvPPO+9jYXMeF\n/8uRniO3F2m9bx93ZRyOV199lTfeeAODwUDjxo3p27dvXXUtUU988cUCKit1DB/eFgeHSu69t3Gt\ndWUyGUFBFyI16nB2hldeuafW+mq1mQEDerNy5dZrhv/OyspkxoypzJ79MkplAeHhSsAbSKek5Cxr\n166nR497CQlpitlcSWioO2VlR8nIOIGtrZk5cx5CEEzAhaODWxMbVXpIyZGRkiOQmiPjXL6AyQQD\nO5joEXllo8sLlJSUoNfr0Wo1jBgxCoBjx46QkJBAXl7eJYKjSAc7j8uJOSlHbxCwUYr4KFyI2RzK\n2jgXzOaL1yEIIs2br0GlmoC3dxuUyjHWzxo2rGTo0LQbvs7z50uYNWsZ8+bNw8UliIceuo9Ro96m\npqYGhULxjzL+lZCQuPu5ZcGxZMkS6/8vXbr0VruTuAZGo/GSTKWFhYXY2Njg4GAJoDRlyouMHTuB\n0FCLAeLixd/Qq9c91hDboigiCAJGo5E+fVqjVp9HozECdeseqlQqWLXqXSorT1FcLKemxoxcLsfN\nzeLJ8NlnHzNixChKS3MoKTnNxx+Pw96+lOPHT5CUlIRabU+7dh1xdnbm1Vfv/6vXVAD69KldGN0q\n5/IEFq5RUqUXMBlllBdqECodKS1R8VtNGmWV1QzsYEJWiwGmo6MjSqWCsjIda9aswtHRiTNnLMG+\nLlx7er7A9mNyjp2VYRYFtPYiEZ52/LkqnN92WUSdQmHG2fkYKtU2Hn7Yg6FDNezYsYRVq85jNBqu\nPPh1IIoi1dU12NqqkMtd+Pbb33B29rykzp1yFZaQkPh3c+fN9SWuyr59e9i7dzcPPzwYs9nArFkf\n0bp1B8aOHQ/AG29MoU+fHrRqFYzJVEy3boGYzac5dy4DQZBRXp5KQcEJ5HIdSqUjY8aMZfbsmWi1\nRXh6ighC/b1c7OxssLOrpro6hpUrD5OfX83UqTMRBAFB0DNt2nAyMw8hl4O7uwdBQcHExOy1tl+z\nZhWvvfamVSzVNzmFAp/9bsux7aFUn3cjL0uLyXTxW35uiicI+ymrrGFkDyNXMj9SqVQMGzaC5ct/\n4KefVgAyZLZeNG9zL8cyA/k9VkZqrqVPL1czUX5ydqwL4+2fAzCbZTg4GHj88VOMHn2WP/74jnXr\n1nDkiD1lZcGcOBEPQLNmETd9jT/9tI2UlBwef/wFfH1bSgafEhIStw0pl8o/nPT0ZE6f3k7r1g1Q\nKOSYTCZKS6sQRTtAiVJZhYOD8rq2v41GE6mpefj6umNrW/fGf9fK7WE2mykpAVGU8csvn7Nnz05k\nMgEbG1uqqiztFAoFPXv2JiUlmaSkMwQFBfPWW+/W+Vz/TkEpzP/Vhu0/tiEvxWJ3IggijRuXERiY\nw5EjXhQUONIw8DzN+x+gWaCB8X0N2FxBr1XpYeWWFI6nKNHL/BD/JupCfM208FPw55qmrFgRiF6v\nQC4389BDKTz9dII13Lher+eTT+Zx/Pgxa9uBA+9n+PCHgOvPpXIhNorZbCY3V8TZuSUeHlJkrhvh\nbn+O3G1I6337uCttOCTqjiVLFtO5c1cUinIUijTat7/ogSCXy3FxuZB/wgxcv0GfQiEnOPhixszs\n7GxiYvag1+sJDQ0jMrJVHV3BlZHJZLi4WF6ABw5YdjLeeed9vL0b8uKLz1BUVERwcDCjR4+jrKyM\nJ56YyLlzqfU6J4CScliwWsXe1S3JS/HCyUnPJ5/EEBGRxw8/fMa+fXvw8QmirGwnWSneKDa3hnsO\nMm2xgKezSAMXES8XEY29yIk0GQlpMoymUJCDl4sZXw8TDd1EfNzMyI22LP0+lEenWIQGQO/embz0\nUjxBQZf+0dvY2PDii6+SnJxEcXEx3t4++NxgCE9RFBk16m1ef/1RGjXqQHBwqGSbISEhcUeQBMc/\nEBsbG7Ky9hAR4YpMVj9uiKdOnWTOnFno9ZYMoOvXr2XAgIFWQ8f6RBRFa/hsrdYRpVKJWu1AUVER\n2dk5ZGZmkpAQB4CDQ/26a5VXwcI/lOzbEEZGQiNsbY18/fUeoqIK+fXXn9i3bw8KhYImTWqAHiQk\n7OBcYgPsbVrTYcgBcoog4/ylL3BPZzP+TjakxfmRfMyZ2CIbiv76V15+8X726ZPJU0+dJCyspNb5\nWVIFNK3189rQ62uwsVFRU2PkzTdfo1mznqjVkpuhhITEnUMSHP8wzGYzrVt74e5efUPtSkpKKC4u\nwt3d/bpe0l99tRC9voaoqCi8vX3YuHEda9f+Qdu27QkIqN/MWnK5nGbNwomPj+ONN6bg6upKRkY6\ngmAJkT5lyovWugMGDLypMSorK0lMPIXJZCI4uCmOjo6X1amotoiNmD+DOHuwCXK5mU8/jSEqyhK/\n4vDhwwA888wLREW14vvvv0UUe5KUtI9TcQ0wVnWjQ4c8AoKL8PAtRlBVk53kxdb1ASyK8bxsPACV\nykTPntk89dRJQkJKb+raAB57bA5vvTWOBg0shqh6veGv7LSweXMsmzcf5LXXXsDBIYI2baTjEwkJ\niTvPf05wbNu2ha5de/zjjOVKS0vYunUzrVs3xtW1kut17RRFkRUrfmDDhrWIIsjlMkaMGEW/fgNq\nbVNdXc358/nIZALPPvsScrmcvLw8YmMPkJGRXu+CA2DSpMf58MP3SU8/R0lJCSqVilGjRhMXbtXN\nVwAAIABJREFUF0dychL29vb07t2XXr1qd72tjbS0VObMeZ+SEsvOgZ2dHc888zwRES2sdSqq4fPV\nSg7s9uXkznAAZs48RI8eOdY6FyKDGo1G63/t7U8yfvwbLF/+AUlJjiQlXRQyMpnZ6spqY2OkX79M\nevfOwt29GhcXPa6uejQaw1VDjF8v48cPx8+vD4WF5ykvz+aBBx5j+/ZPEASB5s3b06HDGLy8/G99\nIAkJCYk64j8nONavX4NMJtKt242/yOqTkpISEhPj6NLF9oaOUbZu3cz69WsRBPDy8iInJ4dly5bS\noIFXrTYZKpUKW1tbqqur2b79T/z8GlldN6+0E1AfODu78PbbMzl7Ngm9vho/P3+cnZ3p2fPW7ovJ\nZOLTT+dTUlJCw4YNsbW1JTk5mc8++5h58z5BrXawio3dWxpzckcEoijwyitxl8WyaN26Denp5/js\ns/m4u3uQl5eHTCbQv787Eyeu59gxV44fd+H4cReOHXOlsNCWiIgihg9PZcCAdLTaS91Xq6urycjI\nQ6PR4uzsfMPXFht7iqioJhQXq+nYcQw+Ph7Y2HgAzTh48CQVFQXI5TYEB1897omEhITEneA/JThE\nUeTxx+9DLjdjMpn+UbscTk4axo1ri63tjd2SQ4cOAvDIIxPp0aMXv/zyI6tXryQ29mCtgkMmkzFo\n0GB++mk53323yFoeFBR8Sy6XN4pCoSAkJLRO+ywsLCAvLw+1Ws27785GoVDw+usvk5GRQVpaGv6N\nw1mwSsmWleGkHAoG4KWX4pg06fRlfQ0c+ABFRUVs3/4neXl5qFRKxo17lCZNmgJGOnXKo1OnPMCS\n/6SiQoGDg/GK84qJ2cuiRV9RXW05KuvWrQfjxk24JKbK1TAaTXz77XqcnMJp2bLrZYafarUatVp9\nvcskISEhcdv5TwgOs9nM118vpFevdri61iAIkJV1DF/fqMsSat1ORFHk888/ZdSoMZSUHMTV9WZu\nh8Wr+YJ4uvgiurq384ABA7Gzs2PHjj/R6/U0bRrKyJEPX/cL8P8xmUwkJJygsrIUrdaF0NBmd2xd\nlUqLC2pNjZ6SkhLkCjt0NU6oXHw5ntGQVUdVbFzRiuzTvigUZmbNimXw4HNX7EsulzNhwiQGDx5G\naWkJHh6e2NtfOTeLIFCr2EhNTWbhws8wm824u7tTWFjIjh3bcHJyYujQB2u9lq+++gMbGxUPPNAd\ng8GGuXMX4ukZfIMrIiEhIfHP4D8hOPR6PUVFBSiVGcjlFjfSjRtXUFCwmtdff/uSumfPJtGokf9t\nS1JVVVVKaupmgoJuzoMgKiqahIQEFi36ik2bNpCefs5afjUEQaBXr3tuykbi/zEYDMyf/yFxccet\nZe3atefxx5++IztIzs7OhEdEciZHw1vfVqJw8AG/9zAXOfDbGi0ZJ/wozHRHrTawYME+OnfOu64+\nb+YI5AKHD8diNpvp1KkLkyc/wZEjh5k//0MOHtx/ieDYvfs4+/cn8OijwxAEN+65ZyQajTtubgFS\nPhMJCYm7nv+E4LCxsWHEiPY4O1/Mg3HPPa0pKRGoqam5JJTznDnvM3hwb5o0CcDfv22dz2X79j/J\nysqkb98elJcnMn58S5TKm78NvXv3JTc3ly1bNpGefg6ZTGDo0AeJjm5Th7OunTVrVhEXdxy1Wk10\ndCsOHDjI/v0xNGnSlHvuuXryvitRUaHg11/9kctFIiKKCQkpwcbGfF1ty6tgX4Icnds09EWepOxr\nSFGGA5VlzojiRfGjVOYyYsTbdOjQGah/UWQ2XxhXiSAIqFTKv8rNxMenEBERSE2NEW/vQO69NxJP\nz27Y2Njg61vvU5OQkJC4bfyrBYfRaCQ5OZnS0mxcXAz8v+eHq6sWFxeR7OzDuLs3R61WU1ycx5Qp\nw3F2NlJZmcnTT3/BrFmfWvOU1AWenp5UVSUDR3F1VXKrt0AmkzF27HgGDBhIYWEhHh6eODk53XR/\np07Z8McfTjg5CURHGwgJKUWtvvJRAUBi4mlEUU54+CecP9+djh2/588/Z5CYmHjDgmPLFm/eeiuK\n3NyLxxZyuRn/gFIaNykmtFkh4RGFBATqsFFBWQVkF8rIKRLIyJNx5owjGad8yDntQ1X5xZTsMpkZ\nW9vT2NmdwN39LA4On3HyZBZr1xYxaNDgG1+kG6Rly0jWrFnJ9u1/kpKSTE5OFgBhYS2YNesHPvro\nfbTaQFq08LmjR3wSEhIS9cm/WnCcOBHHF198xPTpQ6/o+SEIAmVlpxk79jkWL/4Ee/u8v0SADLXa\njlatGlBWlouDQ1CdzcnevoTOnev+q6urqxuurm431bagwIa1a/1YvtyN5ORL5yYIIo0alRMWVkJ4\neBEREcVERBSj0RjIy7Pl+PGnOHZsAPv3N/yr/ht4ednQocOO6x4/O9uOt9+OZMsWS7wIR88SHFx0\nlOY5UV7kQPJZZ5LPOrN5vcVdV6Ey4ORVjFxhorrclupyO/QVNsDFl7WPTzkDB6bTu3cWx49/yaZN\nv9GtW3ceffQxDhwYxqeffsSBA/tvi+Bo0qQpjzwykSVLvmX37jRcXSEysjX33fciEyaE4eCgrfc5\nSEhISNxp/tWCo2XLKD744HkEIbfWOg0buvH118/g5laIIFwUJYIgMGRIF8rKTlFa6oSj4829zC/w\nzTdf0LZtJF5eRcCV7UNEUeT48aOkpqbg4KChXbv2aDT19zIyGgXefbcly5c3tiYpk8tL8fJaj8Gg\nR6drQXV1OGlpGtLSNKxff1GM+PmVk5Vlb21n75CFU0MD2Yl+ZGe/xvLlk2nRIo6OHfNrHV8UYenS\nIObMiaCiQolSZaBJh9OEt0+mZZAZlVLEZFCQl+lIerILyWdcSE1yobjQnoJzHpf0JZOZ8fSspEeP\nHAYNOkdUVKE13sWJE5bcJBeOzi7YQ9xIGqGEhFScnBxo2ND9uttcvE6RHj16ERkZzcqVBwgN7USP\nHvfflIGuhISExN3Kv/6Jd60tapVKgZeXS62fa7UKdu9eTmmpCw899PBNzUEURQwGA3J5OjY2Vz6e\nEUWRxYsXsW3bFmvZmjWrmDp1Op6edR9XQa+X8eyz7diyxQe53Ex4+DH0+lm0727AM+xh9NUG4mPf\nQWHI5okJMzhz2pkTJ5yJj3fh5Ekn0tMdkMvN9O2bga3XYQxOAoIAAZHpxG9tQU6OC2PGdGPQoHM8\n80wC/v7ll4yfn2/LK6+0ZvduS6I07ybZhHaNp0NkBcO6GlHbXqhpAvL/+mchN9eO48ddMBoFvLyq\nKCw8wMaNCykszKW62gOFYiyC0NJav2XLKNavX8vmzZs4c+YMWVkZAERGRtW6PklJmXz55WqmT38U\ng0FLbGwJTk7FNyw49u9P4OeftzN16tPY2rbi6acflHKZSEhI/CeRv/XWW2/dqcErK2vqfYyamlyg\n/Jr1robRWE1+fh5BQZFWt8sbQRAE3N0NNLiKbjh+/CjLli1BpVLRq9c9mExGcnJyyM3NpmPHLtcc\nY+1aX555pj0nTli8Kby8KlEqr/wNvqJCwaRJndi50xuttoalS3cS3PQb8sy2GF2GUVAmo7jCFpVT\nSxQunTmeYUBQl9Gpcw5PTkzhqSdPc++9mUycfJJSdRb5NTJ8Pcy8MtKMWlOByT0dQS5SkuPCqVPO\nLFvWmPR0NU2alOHkVMOWLd5MmNCZxERnbO31tOhzhJbdTjO2XzX92phQXUMGOzgYCQrS0aRJGaWl\nR1mw4E3Ky8swm83odDoOHNhHeHhzXF1dAXB398DBwYH4+DiKi4sxm820bdue0aPHXdGTxmw2U1Oj\nxs+vBYGBvfDwCKFt286EhrYjNzcPW1s9MtmVhazZbCY29hQNG7pjMBhRKl2Jju5DcHA3NBq3OrXR\nUKttLvkbMplMkpipJ/6+1hL1i7Tet4/6WGu1+spedf/69PQ63TEg55r1rofSUhMyWRAnT6bTs2fv\n63p5nD2bhFarIj9/E6tW/UJubg6uri7cf/8Qzp5NYufOHdTU6LG3V5OamkKfPv0YPXoceXm5vPji\ns2g0DixcuOiqY+za5cnEiZ0xGi++bGxtLYGpOnfOpVmzYkJCSrGzM1FSomLChM4cO+aKm1sV33+/\nC6ONjh+3CeiqFZiqMrHX/UpxfhImlT+NQnqicW9JRr6AWRRQykWaB5qJDDKzIVZOVoGMsEYmxt1j\nxN3Nki49p1Dg550K4k85kBLbhIwEX8xmGTKZmcjIQg4ftuwSeDfOI6znMZqHVPJwLwOONxG3asGC\nj4iJiaF795489NDD/PDDd+zatZNOnbrw2GNPXlK3uLiY3NwcHB2d8Pb2vqyvY8eSOHcul44d+9Co\nUecrvrxFUWTu3DewsytnzJg+ABw9eobgYF8cHOwwGk1MnjyHt956HTe3pri7+9WbIeiFtNJbt27i\ntdde5ty5NBo29OHNN9/mgQeG1suY/1WkdOm3F2m9bx9Sevp/KI6OcvLyDrN8+S+0bh2No2PtRzEA\n5eU6nn32caZNG8GCBTPQ6Sw3NSMjnbi4OMzmy909ExNPkZeXS2ysJYLotRKxxcU58+STHTAaZYwe\nnYS7ezVbt3oTF+fK1q0N2brVYswpk5kJDNSh18vJyHCgYcMKFn+3k/2p1Rw8rUQuEwn3SmX/+qmU\nVFcCEBHhyjPjArCzM1BWAbGJcvafknE4Sc7hJMvOQMdmJoZ0MSL/v3ezl6vI0w8YOBBSymrPowS2\nPkNefDCJR/w4fNgdpdJEsy4nadg8he4tjQzsYLqk/Y1QUWGZa2hoGPb29oSEhLFr104qKysuq3ut\neBpqtR2C4IqfX6dadwoEQWDs2GfJyjpFSYklyduKFTsZMWIgAQGuyGSOfPfdOuztb0/Uz8OHYxkz\n5iFrvpesrEwmTx6Po6MTPXr0ui1zkJCQkLgepB2Om+7XiNnsw9q1uygrq+CVV14H4Mcfl1kTj5WW\nZlFVlc7mzd+wfv1amjYNYcSIUaxbt8YaknzSpMdp0MCL+fM/QKe7/OhnwoSJdO9+5RfHuXNqhg3r\nSWGhLfffn8aHHx7kwnsyN9eObdu8OHLEjZMnnTh7Vms18AwKKuWLr3ax9oiRs9mW45BRPYx4uYqU\nl5eTmZmJg4MDDRs2vOzbuShCSo7AoTNyvFzMdI4wW40ztVrLDsf/U1YBv+1WcCxZTnWZHeb8hpid\nc3F01zG8q5F2YdcXY6M2fv/9F37//Ve0Wi3R0W2Jjd2PTqdj+PARDBz4wHX3I4oiBQX2NGrU+boD\nllVXV6NQKO6Y8ae7u4aHHx7LsmVLGDFiFHPmfMw770znyy8X0K/fAL7/fvkdmde/Eekb9+1FWu/b\nh7TDcReg0SiAXAYO9CEvr4Jz5/5ELnfGy0uJXK6jqGgTGo0KtVqguLgYgI4dOxMc3IQOHTpx6NBB\nBMFSJpfLCQpqwtGjR/D19aOkpAi1WkvHjiNxcOjHb7+pKStT4edXTkBAOb6+5ZSWqhg3rguFhbZ0\n7pzLe+8dIqvA8uZXKkBpX8V996fw4IgU5DKLkWhioiO5uXY0aZbPsh0C+SUyWjQ28XAvo9VuwsHB\ngZCQkFqvWxCgsbdIY+/aY3P8P1o1PNLXSHyKmV92iZRqz+JgJzKhn4FAr1vXugMGDOL06VOcPJlg\nNbiNiGh+1Wy5/09KSjbffruOZ5996obEBoCtre21K9Uz5eUWkRoVFY1KpaJ16zZ8+eUCysulh7WE\nhMQ/i5sSHEajkddff52srCwMBgOPPfYYQUFBTJkyBZlMRnBwMNOnT6/ruf4j0Wjs0WjsgRogj6go\nDXCpuvPxscSXWLNmFXp9Nbt37wIsuwUffzwPT09Pjh8/BsDAgW8wb9797NzpxMqVV749crkZe3sj\nOp2KZs2KmD1nH4s2yTmTefkxgEoh0shTpLG3mUCvEgJcSvh6k5LyKoEeLY3c18FELfaPdUpEoJlg\nnxoOnZER7m/GqY5iqalUKl59dSqHDx+ioCAfDw9PIiNbXbdwaNjQna5du+Pn1+WmQ7Hv3r2TN998\nnbS0VPz8GjF9+ju37Tijdes2rF79O++8M50jRw6xZcvGv8rrPkquhISExK1wU0cqv//+O4mJibz2\n2muUlZUxaNAgQkJCmDBhAtHR0UyfPp3OnTvTq9fVH7p385HKjVBdXc27775FWlqqtUyr1VJdXUVN\nzcUU5t26Teb77+eSnW05/3dxqcbPrwJf33K0WgPp6Q6kpGjIzrZHFAX8/MqZt2Abq2LNlJQLNPU1\n4+VixmAUqDFCjRHyiwVyii4VIjJBZEgXI53Cb+044+9c6Ujln05hITRo0A1bW7trV74CR48eZsCA\nezAYLt5HuVzO6tUbadOmfl/67u4acnKKmThxHOvW/WEt79y5G8uW/fyP2IH5tyBt8d9epPW+ffzj\nj1T69etH3759Aaxp3k+ePEl0tCVhWJcuXdi3b981Bcd/BVtbW954YwZ//rmZnJwcXF1d6dGjN+Xl\nOvbv30dNTQ1eXm14552nyc5WExlZwNdf78HRqYbcIoGUHIFincBAVxFfdxGNrYzMDDUZpdUs3WnJ\n1dG/rZFera68W1FRDak5MpJzBM6XCHQMNxHqd8dMd+44oigya9ZSBg7sRUTE4JsWGwCLF3+DwWBg\n6NAHefvtWbzzzpusWPED3377Zb0LDgCFQsGiRUvYuXM7KSnJ+Pr60qNHbymomISExD+Om3oq2dlZ\nHtDl5eU8++yzPP/888yePdv6uVqttnpkSFiwsbHh3nvvu6RMq9X+lfpcyahR3UhN1RISWszkl/fy\n6z4zKTkqKvWXKwgbpYiLtoqcQhn2NiJj7zEQchUBobaF8AAz4QE3Pm+DwcCJE3HodDr8/Brh738T\nndwClZUVxMTso6ysFD+/RkRFRd+ym6kgCERFheHn1xUHh5vPAgtQUmKxz+natTtubm50796TFSt+\nsNrt3A5kMhndu/eke/eet21MCQkJiRvlpr8G5eTk8NRTT/Hwww/Tv39/PvzwQ+tnFRUVaLVSfojr\noaJCwYQJnTl1ypmAgDImv7yXtYcA5LhoRBq5lVCaG4uxKgeNWyiOHi3ILlKQWyTQyNPMuHsMuNTT\nUut0ZcyePfOSo6D77x9Cq1bRnDqVgEKhpFWr6JvO4XItiouLeOed6eTnX4wy2qFDJ0aOHE1y8lkU\nCjlNmoRYBfDfOXkyDY3GHl9fSxj0s2czcXbW4uBgT8eOw/D2vvUcOZGRrdi4cT0zZrxBXNwx/vhj\nlbVcQkJCQuIiN2XDUVBQwJgxY3jzzTdp164dAI8//jjjx4+ndevWTJ8+nXbt2tGv39WzhRqNJhSK\n+k0PnpKyG622pF7HuFmMRoGRI9uwfbsHPj6VfPXdHr7aVIWdDbw+WqAo9yTTpr1hjbEAEBoaynvv\nvQuCAoW8fq09582bx44dO3B1daFJk6YcOHAAs9mMTCZgNlt+bezt7Zg+fTqhoWF1Pv7s2bPZu3cv\nvr6+REZGsnnzZqqrq1GpVNTUWCLjubu78dZbb+Hr6wdYjksu7IB88MEK2rRpTWRkM0TRzPz539O2\nbTM6duxD48at62SO1dXV3HvvvWzfvt1a1qlTJzZv3lyrEJKQkJD4L3JTguO9995jw4YNBAYGWh/w\nU6dO5d1338VgMNC4cWPefffda259/1eMRmtj5swWLFrUFBeXapat2MbvsXpyimRM7G8g3N/M9Omv\nk5ycTLt27WnZMoqfflpOcXExkyc/QefOXet9fs899yQFBQW8884sAgIC+eqrBezaZfGw6dChEyUl\nxZw8mYCbmxvz5n2Kk5O6To1GX3jhaXJz85k58338/QNYvPhbVq7chFYLgYGBVFVVkZ2dg49PQ2bN\nmsPGjQc4ePAUTz/9CHK5J46O/mi1rnU2n9owGo2sW/eH1Uulf/+B1kRx9YlkWHf7kNb69iKt9+3j\nH280OnXqVKZOnXpZ+dKlS2+mu7uesrIy/vxzM0VFRXh7e9OjR29rRtLaWLmyEYsWNUWhMPP55/uI\nz6kip0hBp3AT4f4W75H8/DwAHnroYVxd3cjMzGDt2j/Iy8u74Tmmpqbwww/fkZ2dhbOzK0OGDKdV\nq+irtrlwDRkZ6fj4+JKWdg4ADw8PnnjiaUwmE5MmPUJBQQGlpaU4OdVtdE07O3uOHYMNG/bQs6cv\nK1fGkZ4OLVs6Mm3a2+zbF8/LL39IRUUWKSnlNG/ejdatR+DnF3Jbc4ooFIrbkuZeQkJC4m5GMmW/\nRYqKCnnrrWkUFRVZy2Ji9jJ16lu1io64OGdef93ysn/zzaOoPYrYtVaJp7OZQR0uHp+4u3ug05Xz\n888riIxsxb59ewAQRTN//LESlUpFdHQb3Nwuz2B6/vx5TpyIA0Q8PRswb94cqqstuw86XTkfffQh\nL7/8Os2bt6j12rp168GyZUv56quFfP31Qi7shVl2FrIoLCxAr7ckMrO3t7+hdauNiopqTCYTNjZK\n2rW7nzVrPmLOnLXMmbPWWicpyZ6KilDc3W3R6UTi48HPr5c1WZuEhISExD8PSXDcIr/88iNFRUX4\n+wfQsWMnNm3aQEpKCps3b+S++wZdVr+gwIbHH+9ITY2cBx9M5r77k/ngJxVymciY3kZUyot1R4wY\nxQcfzGTv3j3s3WsRG25urqxevdJa57fffuHll1+jSZOm1rL4+ON89NFc9Ho9AHK5DJPJTPPmLRg/\nfiLr1q1hy5ZNbNy49qqCo2/f/uj1etasWY1er6dBgwbI5XKysrJ45ZUXrPW6d+91zR2d6+Xnn7eh\n0xmYNGkqI0cOQqVqxkcfzeH8+XxCQsLIyMggMzOdjh07YDAYMJtFevfuI4kNCQkJiX84kuC4RbKy\nMgEYOXI0YWHNkMsVLFmymOzszMvq1tTIeOKJDuTm2hMVVcCTzx/lq/VKdFUCgzoY8XG/1JwmLCyc\nN96YwaZNGykv1+Hm5s727VsRBIsb5vnz50lIOMHChZ8xb94nCIKAXq9nwYKP0ev1NGvWDJlMTnx8\nHGBJcObm5k50dBu2bNlkDYtdG4IgcP/9Qxg48AEMBgM2NjaUlZXx/feLOHEiHoVCQefOXRk69ME6\nWk24776+uLm1R6OxJMYbOvTBS/pPTU1h8uRHOHbsKAB9+97LJ58srLPxJSQkJCTqB0lw3CJOTpY4\nDtu3b0GpVBITsxcAZ+fLM8l+8UUIhw+74+lZybhnYvjodwU1RoE2ISa6tTRdsf/GjYN54olgAPbv\n38e2bVtp3rwFjz76GEajkYkTx3L+fD4VFeU4OGjIzc2hvLwCNzc3Xn11GoIg8OSTEykr07FmzSp0\nOh2HD8cC0KjR9cXUkMlk1h0MrVbL008/f2OLdA1SU3MoLCwlICAUD48uqNW1Z8gNCAhk06YdFBcX\noVQq0Wgk92sJCQmJuwFJcNwiAwfeT1zccWJiYoiJiQFAq9XQq9c9l9TLzrbjyy8tSdG6DzvE5jgz\ndjYwtoeBqODrCzHu4GBJQJKefo7MzAzy8nIxGIwoFHJsbCxhrC+4YpaVlZGdnYVMJkevt7iQVlZW\nsX69xRbCx8eH4cNH3OLV1w25uYUkJhbRvv1j1xX1UxAEXFykIxQJCQmJuwlJcPwfRUWFrFixlLS0\nNP7X3p3HRVmvjR//zMLAwACOLCKrqCiIiom7oVhZanl+pj4dc2vx9Ej+LFNJc0tbFJcys0fP0Xwy\nw0r99ah1OlmPHStFTYMURRPzqKmIBggyM2yz/f6gxjiu2QyDw/X+b773Pfd9cTmvmcv7u/n7B/Dg\ng3+66UyO1q3bMGvWXLZu/R/HLJVHHnn0qh/ExYs7UlWlJiq+AKOmlLZRNkbeY/5dm5i1bZtAbGxL\nTp06yQsvpDva+/cfiJdX7eCPkJBQkpI6kZt7sM457dol0q/fvRQWnkevb0rPnr2vu9eG1WpDpbr9\nWR5VVTWcPl1IfHzMLZ0fF3cXqal96mUqqRBCCPe4rXU4nKUhrcNhNBqZM2cGRUU/12mfNGkqXbt2\n+0Mx5OQE8cgj96JSW+n72D8ZOchA36Tb26W1vLyc9957h6NH8/Dy0tCnTypDhgyrs9NpRUUF69b9\nN999tx+A5ORkHn/8L/j53by6+fjjXeTlnSYt7Ul8fEzs3XuIbt0SCAi4/pRXm83GX/+6lSeffIjK\nSm8KCi6yfPl6Fi58FrXaRmCghoqKCj777BPOnTuHXq8nNbU/WVnHuP/++wgPv+cP7WfSWMlaBfVH\ncl2/JN/1p8Gvw+GJ9uzJoqjoZ6KionjqqafZu3cP27Z9yscf/88fKjhsNnjllbsAiE0+wX09jfS7\nzniNWxEQEMDEic/d8BxfX1+efvoZ0tJqa8nfs/dISkpnkpP/Dy1adMdgKOWrr7YQG9sGrdaCl9e1\nPy4KhQKt1heNpiOhobHExECvXo9ht9uprKxg9uw0Llz4FpPpyvohX3+9i+jovowa1VOKDSGEaATq\nb3WkBu7y5drlzzt16kzLlq24997ajbDKym5tWXSDQc0//hFJWVndboEtW2I4fLgp3n6VdL/nOEN6\nW65zBedTKBS3VGxUVlZz9uzPGAwWfH27kpBQu8upv7+e5cvfpWPH/8Bi6cj33xcxf/57jvfZbLVj\nTwwGO2PHziAuriNQ+7TomWfSiI1tTnx8LDt3HuDw4YuEhoYyfvwEEhM7YDRWUlxc45iNIoQQwrPJ\nE45fREZGAbBjx5d4e3tz+HBunfYbsVoVjB9/N/v2haLVWhg27DSPP36c0NAqliyp/RFOSDnKEw9V\n4d0Ahyns2ZPH7t1HmTlzKXp92FXHlUoloaExaDSDKS/XUFICgYEWRo16hRdfnEj79oPr7Lo6YcJT\nfP75Pxyvz5w5A0C7dkMxGLzp23cAH354mMrKC67/44QQQjQIUnD8onv3nnz77V5ycr7jo482AeDv\n78/YsU/c9L0rViSwb18o3t4WKivVrF/fmvffb0VsrIGiIi1Nwkp5fPRpYsPcNlzmhjq+o3j0AAAV\nBUlEQVR37ky/fmkEBt5419cmTfQMHjwCq9XKhQvHmDMnnaSkh9Bqr6wyev58AZ9//g80Gg1ffPE1\nvr6+9OnTg+rqKj788AsSElpTUVGMyYRjwzUhhBCeTwqOXyiVSp59djJZWbvIO1mD3SuM6BZt+f6M\nL3tPgLcGeida0f/bWJh9+0J46612KBR2hj61l6ZB1eR83YrsPdGcPFm7RkS/h3MZ1O32x224Qk2N\nhW+/zSMpqSN6fc/f1bWhUqmIiEgkIiLxqmO/Liam1zclPj4BpVJJdHQ0P/54nNOnTzm2utdqtUyb\nNtM5f4wQQogGr1EVHIWF57l48QLBwSFEREReNb5BpVKhbHoPx/Nqp5j+mFv3/TsOqOiVaOO+zhaa\n6KCkRMPEZ3pgsylp3e04Bs0lDAYITj5Ev4RjnM2Lxte/iufHFfGbSSQNQlFRKdu2ZdOr15NOHUcR\nE9OC4OAQLl68wLBhg/Hx8eHHH4/j6+vH0KHDOX36FM2ahTFx4nMkJrZ32n2FEEI0bI2i4LDb7Xzw\nwXq2bbuyAVifPn0ZN258nemkpwoVbPynD8ezEmkTbaFV68u0aVNGmzblXDLZ+CJbza7DKvYeVdI9\n3saq13pwqUSLPryEkY8fYUA3CyqlnfIKBeUVZsof+IGIYDvN9PXTlWK325kxYxUPPNCNfv06A7XT\nXDt2bEVsbHidcwMD9bz00pvo9c2cGoO3tzerV69l7NhHHZvN+fj4sHLl2wwa9JBT7yWEEOLO0SgK\njqysnWzb9ilKpZK2bdty4sSP7Nz5DZGRUQwaNBiA8gpY+4UXx3bHc+pALP86cOX9CoWdqCgTsbHl\naANNFFlNrPnaj9M/hOGtreHtv+3lrkSz4/zgQDtQf+M1bDYbSqWSkhIfJk9eiFqtwGTyAqwYjYcx\nGv25dMlMYKCKTZt28MADPdFquxAUFOmSeO6+uw9ZWfv56qt/YrPZuPvuPsTGtnTJvYQQQtwZGkXB\ncehQbd/In/88kgcfHMzOnV+zevVfOXQol0GDBmO1wbovvDh/TseZ3JYoFHbS0n7g3Dk/jh0L5ORJ\nf86c0XHmzNULZy1bup+7Eiuuaq+urubUqZNA7f4fztpN9Vrmzn2Hrl07MXTo1DoDOAHS0qYBYLFY\nKCkp4MyZz6ioiCYmpoXL4gEID49g1KixLr2HEEKIO0ejKDjU6tpuk8vlRrLzlfz0sw8Kla+jO+Xv\ne1WcOK/k/PcJWK1Khg8/RXp6HpcvX2bx4gU0bVpAVVUc1dVtiYh4iGPHNJSXR6PXf0Z29qf06jUZ\nne7KaNLTp0+xdOliLl26BEBQUBBTp04nOvrWlvr+vZ56agh6fe+rio26OVDTrFkMr766wiUxCCGE\nEDfSKBb+6tGjNwDbvutE+rSubD8whKbdPqCsyTRWfarmq4NqbGV68nMj0GotTJ6cB8DatWv46afT\nBAf707OnBr1+CxUV44iOHkNKyhBatlzG0aNHWL36b457mc1mli17jUuXLhEaGkpISCglJSW88cZr\nmM3ma8Z3PUeOnOJf/yq44TlWq42AgCTCw6XLQgghRMPVKAqOpKRO3D1oIYd2TeFsXgx7/18PKs5p\nsCm0HP1JhUZt59S3tTMm/vKXfMLCKgH44YcjALzwwhyef34GCQkJADRt2pRly1awcOHrABw4kO0o\nJs6dO0txcTFNmzZl0aKlLF68lMDAJhQV/Uxh4fmbxpqXd5KSEiOlpc3Izb3EoUM/U1JixmSq5qOP\nviI7+xgAly8bmTx5OUVF3kRGdnBuwoQQQggnaxRdKicKFHzwwUiMl/zx8TFTVaVh5+YHWLRoPyn3\nnGH7F1FsPtqUkJBKnnoq3/E+jUaDyWTiwoVCQkJCKC8vB2qnzyoUCtTqK+n7dQ+8X7tpLBYzZnMN\nNpsdi6W2GFEqbz43dsuWnaSmDmLQoM6MH9/ZcW2TyURwcBVeXhouXdJgtaq5774+REX1/l17pQgh\nhBDu4PEFR5nRxuyFrSn4IQpvHwtbt/6TzZtbsHp1PFOn9mDmTB/WrYsDYMqUPPz8rux10rdvP7Zu\n3czSpYtRq1VYLLWLdxUVFTFt2mQMhtod9jp0SHJsrR4REUlUVBRnz57l2Wf/L2CnqqqK6OgYmjdv\nftN4J0x4jPDw++q0KRQKdDodgwePcLRZLBaaNx8gG58JIYS4Izi1S8VutzN37lxGjBjB2LFjOXv2\nrDMvf8tqt0v/Lx4e+ieem3+BnC9q9zNZMD+buLhypk8/xLRptTNXFizoREGBH23bljFs2Ok613n4\n4eEMHPgQarUai8VKWFgYY8c+gU7nx4ULFzCZTLRu3Zr//M8JjveoVCqmTJlGixYtqKqqpKqqitjY\nlkyZ8nydNT/+ndVqo7rajK9voqN4uRG1Wo1WK8WGu9hsNsdTLSGEEDensDvxW3P79u3s2LGDjIwM\ncnNzWbVqFStXrrzu+UVFBmfduo709Od47713iE+ZyJnDGVSU6ejS5X/ZuLHuzq8bN8Yye3YyNpuS\nd9/9hpSUi9e8nsViobq6Cj+/2mmxlZWVnDt3Fm9vbyIjo1Aqr67bbDYbJSUlKBQQFBR8w24Pm83G\n8OFzePXVF+nefcgf+MvdKyTE32X/pg1FcXExU6Y8w44d21GpVAwZMowFC5bg5+dX77E0hnw3FJLr\n+iX5rj+uyHVIiP81253apZKTk0NKSgoASUlJ5OXlOfPyt+Snn07z3nvvoFarqb6UTkWZDl+/XNTq\nIRgMy/H3D3Cc++c/nyIurpySEu/rFhtQ+zRBrb6yBodWqyUurs0N41AqlYSEhNxSzEqlkkWLniMh\n4f5bOl+4h9VqZcyYR8jJyXa0ffjhegwGA++8k+nGyIQQouFzapeK0WjE3/9KZaNWq7HZbM68xU0V\nFJwDIDGxA4aLUeh0Fnr2GI9CUUlpaanjPKvVyqZNG1i79j/YvHkgb721DIOhvF5j/ZXZbKF58+43\nXEdDuF9u7gFycrIJCgoiJyePHTt2o1ar+fTTjykqKnJ3eEII0aA5teDQ6XSYTCbH61+X3K5PLVu2\nQqFQkJt7gCFDljJ69BQqKvahVqsJDr6y/fr776/jk0+2UFpaisFgZN++vSxZkoHFYrnB1Z0rO/sY\nc+f+N2VlgYSF3fiJiXC/X3fCDQ+PJDIyioSEdjRp0uSXY/L4VwghbsSpXSqdO3fmq6++YsCAARw8\neJA2bW78I6rX+zpWAXWWkBB/Xn75ZebMmcOaNc8TEQHR0TBu3JOEhdUWHDU11Xz55XYA5s6dS2ho\nCLNmzebkyZOcP3+a9u3rZ12L5OQ2+Pu3pmvXgfVemLnK9fruPEHfvj3x9fXl8OFcnnjiUcrLyyku\nLiY8PJy77mqHl5dXvcfkyfluaCTX9UvyXX/qK9dOLTj69+/P7t27GTGidvpmRkbGDc8vLb16DxJn\nGD9+EhERsWzf/jk63SUGDuxIUtJdlJfXLuh1+fJlbDYb3t4aWrZsi0qlIiIigrKyMi5cKCE6utIl\ncQFYLFaKispo0iQQs7klXbokUlJiuvkb7wCeP9DLm7feWsWECX/hs88+AyAwsAl/+9taysqqgKp6\njcbz891wSK7rl+S7/tTnoFGnzlL5verjA2UwHAQK67TZ7XYmT55IcXExiYmJ6PVBZGXtRKlUsmTJ\nGzRrFuayeHbvPsSWLbt5+eVlNGvmWcuRN5YviVOnTpKVtRO1Ws0999zn0s/LjTSWfDcEkuv6Jfmu\nP3fsLJU7hUKhIC3tGV5/fSFHjhz5pQ1Gj37M5T8eHTu2o0uX0QQHR7j0PsJ1YmNbEhvrWcWiEEK4\nWqMsOADi4+PJyFhCdvZ+zGYL7dol0qpVa5fdz2CowMdHA8RJsSGEEKLRabQFB0BwcAgDBjzo8vvY\n7XbGjVvI9OmTSU2Nd/n9hBBCiIamURcc9UWhUPDmm7OJiRng7lCEEEIIt/CMuZgNhMlkYc6ctezb\nd9TRZrfbMRrNhIR0x9vb243RCSGEEO4jBYcT7N2bR37+GWy2OMaPn0dERE9KSiA//wzjxy/BbI5B\nr3fPTAYhhBCiIZAuFScoLCzBZAqgS5d4wsNr22rX+ThJWlowUVEd3RugEEII4WZScDhB3779iIy8\nt06bUqmkefPWNG/uupkvQgghxJ1CulRuUXW1mdzcE3z77RFHW27uCS5fthAc3A21Wmo3IYQQ4nqk\n4LiJ778/zoEDP1JZ2YKysjAKC6GsrIbychPLlm3CZotDp9O7O0whhBCiQZP/lt9EaamBigpfUlI6\nOMZiVFdXU1Jyijff/CsREbKuhhBCCHEzUnDcRKdOyTRvfg8KhcLR5u3tTXi4FBpCCCHErZIulevI\nzj5GebkZvb6rrJ8hhBBC/EFScFxDdXUNb7/9d0ymCAIDg90djhBCCHHHky6Va9BovFiwYAEtWiS7\nOxQhhBDCI8gTjn9jt9spK7MTHt7Z3aEIIYQQHkMKjt8wGiv5059ewGIJQ6PRuDscIYQQwmNIwfEb\nOp2WjIypREcnuTsUIYQQwqNIwfEbNpuNmJjusmqoEEII4WRScPwiK+sQFy5YCQtr4+5QhBBCCI8j\nBQe1Tza2bt0FhKNUSkqEEEIIZ7utvgOj0Uh6ejomkwmz2cyMGTNISkri4MGDLFiwALVaTa9evZg4\ncaKz43UJpVLJjBnPEhMj02CFEEIIV7it/86vXbuWXr16kZmZSUZGBi+99BIA8+bNY+nSpXzwwQcc\nOnSIY8eOOTVYV7BYrNTUWNBq29ZZvlwIIYQQznNbBccTTzzBiBEjALBYLHh7e2M0GjGbzURGRgJw\n9913s2fPHudF6gJms4WhQ2dx+rSFZs1i3B2OEEII4bFu2qXy0UcfsW7dujptGRkZtG/fnqKiIqZN\nm8asWbMwmUzodDrHOX5+fpw7d875ETuRl5eaJUtm0L79IHeHIoQQQni0mxYcw4cPZ/jw4Ve15+fn\nk56ezvTp0+nSpQtGoxGj0eg4bjKZCAgIuOG19Xpf1GrVbYR96wwGCAjQ1mkzGivx8/OhtFRFaupD\naLW+Lo2hMQkJ8Xd3CI2K5Lv+SK7rl+S7/tRXrm9r0OiJEyd47rnnWLZsGW3btgVAp9Oh0Wg4e/Ys\nkZGRZGVl3XTQaGlpxe3c/ncrL6+s83rq1P8iNbU3Q4Y8i9FoxWg01Escni4kxJ+iIsllfZF81x/J\ndf2SfNcfV+T6egXMbRUcS5cupaamhvnz52O32wkICGDFihXMmzeP9PR0bDYbvXv3pmPHjn8oaFeZ\nOnUc4eH34Ouru/nJQgghhPjDbqvgWLly5TXbk5KS2Lhx4x8KyFWqq2uoqqoBfGnevC/+/k3cHZIQ\nQgjRaDSaVa7+/vc9vPXWx/j7dycgIMjd4QghhBCNSqPZNOT++3szeHAygYHB7g5FCCGEaHQ8+gmH\nzWbj+PFTGAxmfHw6Exwc4e6QhBBCiEbJowuOkyf/xeuvv41a3ZGgoEh3hyOEEEI0Wh7dpdK6dRyf\nfPK/GI0Wd4cihBBCNGoe/YQDQKvV3vwkIYQQQriUxxccQgghhHA/KTiEEEII4XJScAghhBDC5aTg\nEEIIIYTLScEhhBBCCJeTgkMIIYQQLicFhxBCCCFcTgoOIYQQQricwm63290dhBBCCCE8mzzhEEII\nIYTLScEhhBBCCJeTgkMIIYQQLicFhxBCCCFcTgoOIYQQQricFBxCCCGEcDm1uwNwFbvdzrx588jP\nz0ej0TB//nyioqLcHZbHsFgszJw5k4KCAsxmM2lpabRu3ZoXXngBpVJJXFwcc+fOdXeYHqWkpIRh\nw4axdu1aVCqV5NqFVq9ezY4dOzCbzYwcOZKuXbtKvl3EYrEwffp0CgoKUKvVvPLKK/L5doHc3Fxe\ne+01MjMzOXPmzDXzu2nTJjZu3IiXlxdpaWmkpqY6NQaPfcLx5ZdfUlNTw4YNG5g6dSoZGRnuDsmj\nfPLJJ+j1et5//33WrFnDK6+8QkZGBlOmTGH9+vXYbDa+/PJLd4fpMSwWC3PnzsXHxwdAcu1C+/fv\n58CBA2zYsIHMzEwKCwsl3y70zTffYLPZ2LBhAxMmTOCNN96QfDvZmjVrmD17NmazGbj290dxcTGZ\nmZls3LiRNWvW8PrrrzvOdxaPLThycnJISUkBICkpiby8PDdH5FkGDhzIpEmTALBarahUKo4ePUqX\nLl0A6NOnD3v37nVniB5l0aJFPProo4SGhmK32yXXLpSVlUWbNm2YMGECTz/9NKmpqZJvF2rRogVW\nqxW73Y7BYECtVku+nSwmJoYVK1Y4Xh85cqROfvfs2cOhQ4dITk5GrVaj0+lo0aIF+fn5To3DYwsO\no9GIv7+/47VarcZms7kxIs+i1Wrx9fXFaDQyadIkJk+ezG8XrfXz88NgMLgxQs+xefNmgoKC6N27\ntyPHv/0sS66dq7S0lLy8PJYvX868efNIT0+XfLuQn58f586dY8CAAbz44ouMGTNGvkucrH///qhU\nKsfrf8+v0WjEZDLV+c309fV1et49dgyHTqfDZDI5XttsNpRKj62v3KKwsJCJEycyevRoHnzwQZYs\nWeI4ZjKZCAgIcGN0nmPz5s0oFAp2795Nfn4+06dPp7S01HFccu1cTZo0oVWrVqjVamJjY/H29ubi\nxYuO45Jv53r33XdJSUlh8uTJXLx4kTFjxtR5lC/5dr7f/hb+ml+dTofRaLyq3an3derVGpDOnTvz\nzTffAHDw4EHatGnj5og8S3FxMePGjeP555/n4YcfBiAhIYHvvvsOgJ07d5KcnOzOED3G+vXryczM\nJDMzk/j4eBYvXkxKSork2kWSk5PZtWsXABcvXqSyspIePXqwf/9+QPLtbIGBgeh0OgD8/f2xWCy0\na9dO8u1C7dq1u+r7o0OHDuTk5FBTU4PBYODkyZPExcU59b4e+4Sjf//+7N69mxEjRgDIoFEnW7Vq\nFeXl5axcuZIVK1agUCiYNWsWr776KmazmVatWjFgwAB3h+mxpk+fzpw5cyTXLpCamkp2djbDhw93\nzHaLiIhwDLqTfDvXY489xsyZMxk1ahQWi4X09HQSExMl3y50re8PhULBmDFjGDlyJHa7nSlTpqDR\naJx6X9ktVgghhBAu57FdKkIIIYRoOKTgEEIIIYTLScEhhBBCCJeTgkMIIYQQLicFhxBCCCFcTgoO\nIYQQQricFBxCCCGEcDkpOIQQQgjhcv8foCC1PdXtEssAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1962e7dc6d8>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "import code.book_plots as bp\n",
    "from numpy.random import randn\n",
    "from filterpy.kalman import EnsembleKalmanFilter as EnKF\n",
    "from filterpy.kalman import KalmanFilter\n",
    "from filterpy.common import Q_discrete_white_noise\n",
    "\n",
    "np.random.seed(1234)\n",
    "\n",
    "def hx(x):\n",
    "    return np.array([x[0]])\n",
    "\n",
    "def fx(x, dt):\n",
    "    return np.dot(F, x)\n",
    "    \n",
    "F = np.array([[1., 1.],[0., 1.]])\n",
    "\n",
    "x = np.array([0., 1.])\n",
    "P = np.eye(2) * 100.\n",
    "enf = EnKF(x=x, P=P, dim_z=1, dt=1., N=20, hx=hx, fx=fx)\n",
    "\n",
    "std_noise = 10.\n",
    "enf.R *= std_noise**2\n",
    "enf.Q = Q_discrete_white_noise(2, 1., .001)\n",
    "\n",
    "kf = KalmanFilter(dim_x=2, dim_z=1)\n",
    "kf.x = np.array([x]).T\n",
    "kf.F = F.copy()\n",
    "kf.P = P.copy()\n",
    "kf.R = enf.R.copy()\n",
    "kf.Q = enf.Q.copy()\n",
    "kf.H = np.array([[1., 0.]])\n",
    "\n",
    "measurements = []\n",
    "results = []\n",
    "ps = []\n",
    "kf_results = []\n",
    "\n",
    "zs = []\n",
    "for t in range (0,100):\n",
    "    # create measurement = t plus white noise\n",
    "    z = t + randn()*std_noise\n",
    "    zs.append(z)\n",
    "\n",
    "    enf.predict()\n",
    "    enf.update(np.asarray([z]))\n",
    "    \n",
    "    kf.predict()\n",
    "    kf.update(np.asarray([[z]]))\n",
    "\n",
    "    # save data\n",
    "    results.append (enf.x[0])\n",
    "    kf_results.append (kf.x[0,0])\n",
    "    measurements.append(z)\n",
    "    ps.append(3*(enf.P[0,0]**.5))\n",
    "\n",
    "results = np.asarray(results)\n",
    "ps = np.asarray(ps)\n",
    "\n",
    "plt.plot(results, label='EnKF')\n",
    "plt.plot(kf_results, label='KF', c='b', lw=2)\n",
    "bp.plot_measurements(measurements)\n",
    "plt.plot (results - ps, c='k',linestyle=':', lw=1, label='1$\\sigma$')\n",
    "plt.plot(results + ps, c='k', linestyle=':', lw=1)\n",
    "plt.fill_between(range(100), results - ps, results + ps, facecolor='y', alpha=.3)\n",
    "plt.legend(loc='best');"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "It can be a bit difficult to see, but the KF and EnKF start off slightly different, but soon converge to producing nearly the same values. The EnKF is a suboptimal filter, so it will not produce the optimal solution that the KF produces. However, I deliberately chose $N$ to be quite small (20) to guarantee that the EnKF output is quite suboptimal. If I chose a more reasonable number such as 2000 you would be unable to see the difference between the two filter outputs on this graph."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Outstanding Questions"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "All of this should be considered as *my* questions, not lingering questions in the literature. However, I am copying equations directly from well known sources in the field, and they do not address the discrepancies.\n",
    "\n",
    "First, in Brown [2] we have all sums multiplied by $\\frac{1}{N}$, as in \n",
    "\n",
    "$$ \\hat{x} = \\frac{1}{N}\\sum_{i=1}^N\\chi_k^{(i)}$$\n",
    "\n",
    "The same equation in Crassidis [3] reads (I'll use the same notation as in Brown, although Crassidis' is different)\n",
    "\n",
    "$$ \\hat{x} = \\frac{1}{N-1}\\sum_{i=1}^N\\chi_k^{(i)}$$\n",
    "\n",
    "The same is true in both sources for the sums in the computation for the covariances. Crassidis, in the context of talking about the filter's covariance, states that $N-1$ is used to ensure an unbiased estimate. Given the following standard equations for the mean and standard deviation (p.2 of Crassidis), this makes sense for the covariance.\n",
    "\n",
    "$$\n",
    "\\begin{aligned}\n",
    "\\mu &= \\frac{1}{N}\\sum_{i=1}^N[\\tilde{z}(t_i) - \\hat{z}(t_i)] \\\\\n",
    " \\sigma^2 &= \\frac{1}{N-1}\\sum_{i=1}^N\\{[\\tilde{z}(t_i) - \\hat{z}(t_i)] - \\mu\\}^2\n",
    "\\end{aligned}\n",
    "$$\n",
    "\n",
    "However, I see no justification or reason to use $N-1$ to compute the mean. If I use $N-1$ in the filter for the mean the filter does not converge and the state essentially follows the measurements without any filtering. However, I do see a reason to use it for the covariance as in Crassidis, in contrast to Brown. Again, I support my decision empirically - $N-1$ works in the implementation of the filter, $N$ does not.\n",
    "\n",
    "My second question relates to the use of the $\\mathbf{R}$ matrix. In Brown $\\mathbf{R}$ is added to $\\mathbf{P}_{zz}$ whereas it isn't in Crassidis and other sources. I have read on the web notes by other implementers that adding R helps the filter, and it certainly seems reasonable and necessary to me, so this is what I do. \n",
    "\n",
    "My third question relates to the computation of the covariance $\\mathbf{P}$. Again, we have different equations in Crassidis and Brown. I have chosen the implementation given in Brown as it seems to give me the  behavior that I expect (convergence of $\\mathbf{P}$ over time) and it closely compares to the form in the linear KF. In contrast I find the equations in Crassidis do not seem to converge much.\n",
    "\n",
    "My fourth question relates to the state estimate update. In Brown we have \n",
    "\n",
    "$$\\boldsymbol\\chi  = \\boldsymbol\\chi + \\mathbf{K}[\\mathbf{z} -\\mathbf{z}_{mean} + \\mathbf{v}_R]$$ \n",
    "\n",
    "whereas in Crassidis we have\n",
    "\n",
    "$$\\boldsymbol\\chi  = \\boldsymbol\\chi + \\mathbf{K}[\\mathbf{z} -\\boldsymbol\\chi_h + \\mathbf{v}_R]$$ \n",
    "\n",
    "To me the Crassidis equation seems logical, and it produces a filter that performs like the linear KF for linear problems, so that is the formulation that I have chosen. \n",
    "\n",
    "I am not comfortable saying either book is wrong; it is quite possible that I missed some point that makes each set of equations work. I can say that when I implemented them as written I did not get a filter that worked. I define \"work\" as performs essentially the same as the linear KF for linear problems. Between reading implementation notes on the web and reasoning about various issues I have chosen the implementation in this chapter, which does in fact seem to work correctly. I have yet to explore the significant amount of original literature that will likely definitively explain the discrepancies. I would like to leave this here in some form even if I do find an explanation that reconciles the various differences, as if I got confused by these books than probably others will as well."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## References"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "- [1] Mackenzie, Dana. *Ensemble Kalman Filters Bring Weather Models Up to Date* Siam News,  Volume 36, Number 8, October 2003. http://www.siam.org/pdf/news/362.pdf\n",
    "\n",
    "- [2]  Brown, Robert Grover, and Patrick Y.C. Hwang. *Introduction to Random Signals and Applied Kalman Filtering, With MATLAB® excercises and solutions.* Wiley, 2012.\n",
    "\n",
    "- [3] Crassidis, John L., and John L. Junkins. *Optimal estimation of dynamic systems*. CRC press, 2011."
   ]
  }
 ],
 "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.5.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}