S74_drone_perception.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "view-in-github",
    "tags": [
     "no-tex"
    ]
   },
   "source": [
    "<a href=\"https://colab.research.google.com/github/gtbook/robotics/blob/main/S74_drone_perception.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/></a>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "id": "JoW4C_OkOMhe",
    "tags": [
     "remove-cell"
    ]
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Note: you may need to restart the kernel to use updated packages.\n"
     ]
    }
   ],
   "source": [
    "%pip install -q -U gtbook\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "id": "10-snNDwOSuC",
    "tags": [
     "remove-cell"
    ]
   },
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import plotly.graph_objects as go\n",
    "try:\n",
    "    import google.colab\n",
    "except:\n",
    "    import plotly.io as pio\n",
    "    pio.renderers.default = \"png\"\n",
    "\n",
    "import gtsam\n",
    "from gtbook.display import show\n",
    "import gtsam.utils.plot as gtsam_plot\n"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {
    "id": "nAvx4-UCNzt2"
   },
   "source": [
    "```{index} perception; visual SLAM\n",
    "```\n",
    "\n",
    "# Visual SLAM\n",
    "\n",
    "> When we do SLAM with cameras, we call it visual SLAM.\n",
    "\n",
    "<img src=\"Figures7/S74-Autonomous_camera_drone-04.jpg\" alt=\"Splash image with strange drone presumably thinking about maps\" width=\"40%\" align=center style=\"vertical-align:middle;margin:10px 0px\">"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "ca30f321",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "```{index} pair: structure from motion; SfM\n",
    "```\n",
    "Visual SLAM *almost* the same problem as the \"SLAM with Landmarks\" problem from Section 6.4, but we will now use 3D poses and 3D points, and we will use the inertial measurements from the IMU to inform us about the movement of the platform. This is a typical pipeline on drones, phones, etc., everything with an IMU. \n",
    "\n",
    "Many unmanned aerial vehicles use cameras as their main sensors in addition to inertial sensing, as they are an incredibly rich sensing modality, do not consume a whole lot of energy, and are quite small. Especially the latter two considerations are very important in the weight and energy constrained world of autonomous drones. Furthermore, in military applications there is another reason to favor cameras over active sensors like LIDAR: cameras do not *emit* energy and hence are stealthier.\n",
    "\n",
    "Below we first introduce how we can estimate the trajectory of a drone using an inertial measurement unit or IMU, which we introduced in Section 7.2. Here we integrate IMU measurements over time to yield an estimate of the trajectory of the drone over time. Below that we discuss how we can build 3D maps from images alone (!), in a process called **structure from motion** or **SfM**. We then discuss Visual SLAM, which combines both IMU integration and SfM in a real-time, incremental pipeline. Implementing visual SLAM is beyond our scope, but it should be easy to imagine how inertial sensing and 3D reconstruction done together and in real-time is an incredible asset to autonomous flight."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Integrating the IMU\n",
    "\n",
    "We start off by encoding the information provided by the inertial measurement unit (IMU) in a factor graph.\n",
    "\n",
    "As an example in code, below we define a scenario with forward velocity 2m/s, while pitching up\n",
    "with angular velocity 30 degree/sec (negative in FLU). We set up a simulation using a built-in GTSAM class called the `ScenarioRunner`. To use this class, we need to define four things:\n",
    "\n",
    "- a scenario: below we use a \"constant twist\" scenario, which is very much like the example in Section 7.2, a constant angular and linear velocity specified in the body frame\n",
    "- noise parameters for the gyroscope and accelerometer\n",
    "- the sample frequency, which we set at 10Hz\n",
    "- bias parameters, to simulate how a bias actually corrupts the IMU measurements\n",
    "\n",
    "We then create a `runner` with these:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "e0615c59",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "outputs": [],
   "source": [
    "# Set up the scenario with a constant velocity model, both angular (W) and linear (V)\n",
    "wb = np.array([0, -np.radians(30), 0])\n",
    "vb = np.array([2, 0, 0])\n",
    "scenario = gtsam.ConstantTwistScenario(wb, vb)\n",
    "\n",
    "# Create noise parameters\n",
    "noise_parameters = gtsam.PreintegrationParams.MakeSharedU(9.8)\n",
    "kGyroSigma = np.radians(0.5) / 60  # 0.5 degree ARW\n",
    "kAccelSigma = 0.1 / 60  # 10 cm VRW\n",
    "I3 = np.identity(3, float)\n",
    "noise_parameters.setGyroscopeCovariance(kGyroSigma**2 * I3)\n",
    "noise_parameters.setAccelerometerCovariance(kAccelSigma**2 * I3)\n",
    "noise_parameters.setIntegrationCovariance(0.0000001**2 * I3) # ignore this for now\n",
    "\n",
    "# Set the sample time\n",
    "dt = 0.1 # 10 Hz\n",
    "\n",
    "# Create bias parameters\n",
    "accBias = np.array([0, 0.1, 0]) # bias in accelerometer\n",
    "gyroBias = np.array([0, 0, 0]) # bias in gyroscope\n",
    "bias = gtsam.imuBias.ConstantBias(accBias, gyroBias)\n",
    "\n",
    "# Instantiate IMU simulation class as `runner`\n",
    "runner = gtsam.ScenarioRunner(scenario, noise_parameters, dt, bias)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Our approach will be to use the `runner` class to generate noisy measurements, and then use GTSAM to optimize for the optimal trajectory given those measurements.\n",
    "\n",
    "With the `scenario` class, we can create some (good) initial values by calling the `scenario.navState` at $t=0$ as well as every 5 steps, which is twice per second, given our sampling time $dt=0.1$:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Adding initial state with position [0. 0. 0.]\n",
      "Adding navState 1 at t=0.5 with position [1.  0.  0.1]\n",
      "Adding navState 2 at t=1.0 with position [1.9 0.  0.5]\n",
      "Adding navState 3 at t=1.5 with position [2.7 0.  1.1]\n",
      "Adding navState 4 at t=2.0 with position [3.3 0.  1.9]\n",
      "Adding navState 5 at t=2.5 with position [3.7 0.  2.8]\n",
      "Adding navState 6 at t=3.0 with position [3.8 0.  3.8]\n"
     ]
    }
   ],
   "source": [
    "# Define key naming scheme for GTSAM\n",
    "X = gtsam.symbol_shorthand.X # for poses\n",
    "V = gtsam.symbol_shorthand.V # for velocities\n",
    "B0 = gtsam.symbol('b', 0) # for bias\n",
    "\n",
    "# Create values and add initial state\n",
    "values = gtsam.Values()\n",
    "state = scenario.navState(0)\n",
    "print(f\"Adding initial state with position {np.round(state.pose().translation(),1)}\")\n",
    "values.insert(X(0), state.pose())\n",
    "values.insert(V(0), state.velocity())\n",
    "\n",
    "i, T = 0, 3.5 # run scenario at 10Hz for 3.5 seconds, adding state every 5th iteration\n",
    "for j, t in enumerate(np.arange(0, T, dt)):\n",
    "    if j>0 and j % 5 == 0:\n",
    "        # create IMU factor\n",
    "        i += 1\n",
    "        state = scenario.navState(t)\n",
    "        print(f\"Adding navState {i} at t={t} with position {np.round(state.pose().translation(),1)}\")\n",
    "        values.insert(X(i), state.pose().retract(np.random.normal(scale=0.1, size=(6,))))\n",
    "        values.insert(V(i), state.velocity() + np.random.normal(scale=0.1, size=(3,)))\n",
    "\n",
    "# Add bias, as we will need it below\n",
    "values.insert(B0, bias)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now we will create the most complicated factor graph of the entire book, as it involves a factor that is connected to no less than five unknowns at every step. GTSAM contains a powerful, state of the art factor called the `ImuFactor` that allows us to accumulate IMU measurements and periodically add the resulting relative pose information to the graph. \n",
    "\n",
    "The accumulation of (high-rate) IMU measurements is done in a class called `PreintegratedImuMeasurements`, abbreviated \"PIM\", which is done at 10Hz.\n",
    "Every 5 steps we add an `ImuFactor` that is connected to two successive poses with keys `X(i)` and `X(i+1)` and their corresponding velocities with keys `V(i)` and `V(i+1)`, as well as the (single) bias unknown with key `B(0)`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "ed36f3f8",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "outputs": [],
   "source": [
    "# Create graph\n",
    "imu_graph = gtsam.NonlinearFactorGraph()\n",
    "Iso = gtsam.noiseModel.Isotropic\n",
    "state = scenario.navState(0)\n",
    "imu_graph.addPriorPose3(X(0), state.pose(), Iso.Sigma(6, 0.1))\n",
    "imu_graph.addPriorVector(V(0), state.velocity(), Iso.Sigma(3, 0.1))\n",
    "imu_graph.addPriorConstantBias(B0, bias, Iso.Sigma(6, 0.1))\n",
    "\n",
    "i, T = 0, 3.5 # run scenario at 10Hz for 3.5 seconds, adding state every 5th iteration\n",
    "pim = gtsam.PreintegratedImuMeasurements(noise_parameters, bias)\n",
    "for j, t in enumerate(np.arange(0, T, dt)):\n",
    "    measuredOmega = runner.measuredAngularVelocity(t)\n",
    "    measuredAcc = runner.measuredSpecificForce(t)\n",
    "    pim.integrateMeasurement(measuredAcc, measuredOmega, dt)\n",
    "    if j>0 and j % 5 == 0:\n",
    "        # create IMU factor\n",
    "        factor = gtsam.ImuFactor(X(i), V(i), X(i + 1), V(i + 1), B0, pim)\n",
    "        imu_graph.push_back(factor)\n",
    "        i += 1\n",
    "        pim.resetIntegration()\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The graph displayed below shows (a bit messily) that the IMUFactors tie together the `Pose3` states $x_k$, the velocities $v_k$ and the unknown IMU biases $b_0$."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "78794d08",
   "metadata": {
    "slideshow": {
     "slide_type": "-"
    }
   },
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<?xml version=\"1.0\" encoding=\"UTF-8\" standalone=\"no\"?>\n",
       "<!DOCTYPE svg PUBLIC \"-//W3C//DTD SVG 1.1//EN\"\n",
       " \"http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd\">\n",
       "<!-- Generated by graphviz version 9.0.0 (0)\n",
       " -->\n",
       "<!-- Pages: 1 -->\n",
       "<svg width=\"360pt\" height=\"145pt\"\n",
       " viewBox=\"0.00 0.00 360.00 144.85\" xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\">\n",
       "<g id=\"graph0\" class=\"graph\" transform=\"scale(0.622362 0.622362) rotate(0) translate(4 228.75)\">\n",
       "<polygon fill=\"white\" stroke=\"none\" points=\"-4,4 -4,-228.75 574.44,-228.75 574.44,4 -4,4\"/>\n",
       "<!-- var7061644215716937728 -->\n",
       "<g id=\"node1\" class=\"node\">\n",
       "<title>var7061644215716937728</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"240.68\" cy=\"-79.76\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"240.68\" y=\"-74.71\" font-family=\"Times,serif\" font-size=\"14.00\">b0</text>\n",
       "</g>\n",
       "<!-- factor2 -->\n",
       "<g id=\"node18\" class=\"node\">\n",
       "<title>factor2</title>\n",
       "<ellipse fill=\"black\" stroke=\"black\" cx=\"219.69\" cy=\"-1.8\" rx=\"1.8\" ry=\"1.8\"/>\n",
       "</g>\n",
       "<!-- var7061644215716937728&#45;&#45;factor2 -->\n",
       "<g id=\"edge3\" class=\"edge\">\n",
       "<title>var7061644215716937728&#45;&#45;factor2</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M235.81,-61.68C230.48,-41.89 222.4,-11.89 220.24,-3.85\"/>\n",
       "</g>\n",
       "<!-- factor3 -->\n",
       "<g id=\"node19\" class=\"node\">\n",
       "<title>factor3</title>\n",
       "<ellipse fill=\"black\" stroke=\"black\" cx=\"132.96\" cy=\"-99.08\" rx=\"1.8\" ry=\"1.8\"/>\n",
       "</g>\n",
       "<!-- var7061644215716937728&#45;&#45;factor3 -->\n",
       "<g id=\"edge8\" class=\"edge\">\n",
       "<title>var7061644215716937728&#45;&#45;factor3</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M214.33,-84.49C186.42,-89.49 145.06,-96.91 135.15,-98.69\"/>\n",
       "</g>\n",
       "<!-- factor4 -->\n",
       "<g id=\"node20\" class=\"node\">\n",
       "<title>factor4</title>\n",
       "<ellipse fill=\"black\" stroke=\"black\" cx=\"180.96\" cy=\"-160.41\" rx=\"1.8\" ry=\"1.8\"/>\n",
       "</g>\n",
       "<!-- var7061644215716937728&#45;&#45;factor4 -->\n",
       "<g id=\"edge13\" class=\"edge\">\n",
       "<title>var7061644215716937728&#45;&#45;factor4</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M228.59,-96.08C213.33,-116.7 188.31,-150.49 182.29,-158.61\"/>\n",
       "</g>\n",
       "<!-- factor5 -->\n",
       "<g id=\"node21\" class=\"node\">\n",
       "<title>factor5</title>\n",
       "<ellipse fill=\"black\" stroke=\"black\" cx=\"231.24\" cy=\"-156.49\" rx=\"1.8\" ry=\"1.8\"/>\n",
       "</g>\n",
       "<!-- var7061644215716937728&#45;&#45;factor5 -->\n",
       "<g id=\"edge18\" class=\"edge\">\n",
       "<title>var7061644215716937728&#45;&#45;factor5</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M238.44,-97.94C236.04,-117.44 232.45,-146.67 231.48,-154.5\"/>\n",
       "</g>\n",
       "<!-- factor6 -->\n",
       "<g id=\"node22\" class=\"node\">\n",
       "<title>factor6</title>\n",
       "<ellipse fill=\"black\" stroke=\"black\" cx=\"288.07\" cy=\"-110.42\" rx=\"1.8\" ry=\"1.8\"/>\n",
       "</g>\n",
       "<!-- var7061644215716937728&#45;&#45;factor6 -->\n",
       "<g id=\"edge23\" class=\"edge\">\n",
       "<title>var7061644215716937728&#45;&#45;factor6</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M260.23,-92.41C270.83,-99.27 282.58,-106.87 286.62,-109.48\"/>\n",
       "</g>\n",
       "<!-- factor7 -->\n",
       "<g id=\"node23\" class=\"node\">\n",
       "<title>factor7</title>\n",
       "<ellipse fill=\"black\" stroke=\"black\" cx=\"325.18\" cy=\"-86.26\" rx=\"1.8\" ry=\"1.8\"/>\n",
       "</g>\n",
       "<!-- var7061644215716937728&#45;&#45;factor7 -->\n",
       "<g id=\"edge28\" class=\"edge\">\n",
       "<title>var7061644215716937728&#45;&#45;factor7</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M267.88,-81.85C288.72,-83.45 315.19,-85.49 322.95,-86.09\"/>\n",
       "</g>\n",
       "<!-- factor8 -->\n",
       "<g id=\"node24\" class=\"node\">\n",
       "<title>factor8</title>\n",
       "<ellipse fill=\"black\" stroke=\"black\" cx=\"358.17\" cy=\"-90.19\" rx=\"1.8\" ry=\"1.8\"/>\n",
       "</g>\n",
       "<!-- var7061644215716937728&#45;&#45;factor8 -->\n",
       "<g id=\"edge33\" class=\"edge\">\n",
       "<title>var7061644215716937728&#45;&#45;factor8</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M267.93,-82.18C298.84,-84.93 346.27,-89.14 356.28,-90.03\"/>\n",
       "</g>\n",
       "<!-- var8502796096475496448 -->\n",
       "<g id=\"node2\" class=\"node\">\n",
       "<title>var8502796096475496448</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"111.44\" cy=\"-206.75\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"111.44\" y=\"-201.7\" font-family=\"Times,serif\" font-size=\"14.00\">v0</text>\n",
       "</g>\n",
       "<!-- factor1 -->\n",
       "<g id=\"node17\" class=\"node\">\n",
       "<title>factor1</title>\n",
       "<ellipse fill=\"black\" stroke=\"black\" cx=\"1.8\" cy=\"-195.51\" rx=\"1.8\" ry=\"1.8\"/>\n",
       "</g>\n",
       "<!-- var8502796096475496448&#45;&#45;factor1 -->\n",
       "<g id=\"edge2\" class=\"edge\">\n",
       "<title>var8502796096475496448&#45;&#45;factor1</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M84.34,-203.97C55.93,-201.06 14.04,-196.77 4.01,-195.74\"/>\n",
       "</g>\n",
       "<!-- var8502796096475496448&#45;&#45;factor3 -->\n",
       "<g id=\"edge5\" class=\"edge\">\n",
       "<title>var8502796096475496448&#45;&#45;factor3</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M115.09,-188.49C120.56,-161.13 130.46,-111.59 132.56,-101.08\"/>\n",
       "</g>\n",
       "<!-- var8502796096475496449 -->\n",
       "<g id=\"node3\" class=\"node\">\n",
       "<title>var8502796096475496449</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"183.44\" cy=\"-206.75\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"183.44\" y=\"-201.7\" font-family=\"Times,serif\" font-size=\"14.00\">v1</text>\n",
       "</g>\n",
       "<!-- var8502796096475496449&#45;&#45;factor3 -->\n",
       "<g id=\"edge7\" class=\"edge\">\n",
       "<title>var8502796096475496449&#45;&#45;factor3</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M175.33,-189.45C162.62,-162.35 138.94,-111.84 133.92,-101.12\"/>\n",
       "</g>\n",
       "<!-- var8502796096475496449&#45;&#45;factor4 -->\n",
       "<g id=\"edge10\" class=\"edge\">\n",
       "<title>var8502796096475496449&#45;&#45;factor4</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M182.46,-188.44C181.93,-178.52 181.33,-167.36 181.08,-162.67\"/>\n",
       "</g>\n",
       "<!-- var8502796096475496450 -->\n",
       "<g id=\"node4\" class=\"node\">\n",
       "<title>var8502796096475496450</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"255.44\" cy=\"-206.75\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"255.44\" y=\"-201.7\" font-family=\"Times,serif\" font-size=\"14.00\">v2</text>\n",
       "</g>\n",
       "<!-- var8502796096475496450&#45;&#45;factor4 -->\n",
       "<g id=\"edge12\" class=\"edge\">\n",
       "<title>var8502796096475496450&#45;&#45;factor4</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M235.48,-194.33C216.52,-182.54 189.89,-165.96 182.77,-161.53\"/>\n",
       "</g>\n",
       "<!-- var8502796096475496450&#45;&#45;factor5 -->\n",
       "<g id=\"edge15\" class=\"edge\">\n",
       "<title>var8502796096475496450&#45;&#45;factor5</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M246.98,-189.17C241.41,-177.62 234.72,-163.73 232.23,-158.56\"/>\n",
       "</g>\n",
       "<!-- var8502796096475496451 -->\n",
       "<g id=\"node5\" class=\"node\">\n",
       "<title>var8502796096475496451</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"327.44\" cy=\"-206.75\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"327.44\" y=\"-201.7\" font-family=\"Times,serif\" font-size=\"14.00\">v3</text>\n",
       "</g>\n",
       "<!-- var8502796096475496451&#45;&#45;factor5 -->\n",
       "<g id=\"edge17\" class=\"edge\">\n",
       "<title>var8502796096475496451&#45;&#45;factor5</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M306.09,-195.6C280.85,-182.41 241.2,-161.69 232.82,-157.32\"/>\n",
       "</g>\n",
       "<!-- var8502796096475496451&#45;&#45;factor6 -->\n",
       "<g id=\"edge20\" class=\"edge\">\n",
       "<title>var8502796096475496451&#45;&#45;factor6</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M320.22,-189.08C310.1,-164.31 292.5,-121.27 288.77,-112.15\"/>\n",
       "</g>\n",
       "<!-- var8502796096475496452 -->\n",
       "<g id=\"node6\" class=\"node\">\n",
       "<title>var8502796096475496452</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"399.44\" cy=\"-206.75\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"399.44\" y=\"-201.7\" font-family=\"Times,serif\" font-size=\"14.00\">v4</text>\n",
       "</g>\n",
       "<!-- var8502796096475496452&#45;&#45;factor6 -->\n",
       "<g id=\"edge22\" class=\"edge\">\n",
       "<title>var8502796096475496452&#45;&#45;factor6</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M383.01,-192.54C354.93,-168.25 299.77,-120.54 289.67,-111.81\"/>\n",
       "</g>\n",
       "<!-- var8502796096475496452&#45;&#45;factor7 -->\n",
       "<g id=\"edge25\" class=\"edge\">\n",
       "<title>var8502796096475496452&#45;&#45;factor7</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M388.96,-189.74C370.41,-159.64 333.1,-99.11 326.27,-88.02\"/>\n",
       "</g>\n",
       "<!-- var8502796096475496453 -->\n",
       "<g id=\"node7\" class=\"node\">\n",
       "<title>var8502796096475496453</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"471.44\" cy=\"-206.75\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"471.44\" y=\"-201.7\" font-family=\"Times,serif\" font-size=\"14.00\">v5</text>\n",
       "</g>\n",
       "<!-- var8502796096475496453&#45;&#45;factor7 -->\n",
       "<g id=\"edge27\" class=\"edge\">\n",
       "<title>var8502796096475496453&#45;&#45;factor7</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M454.09,-192.46C418.37,-163.03 338.95,-97.61 326.76,-87.56\"/>\n",
       "</g>\n",
       "<!-- var8502796096475496453&#45;&#45;factor8 -->\n",
       "<g id=\"edge30\" class=\"edge\">\n",
       "<title>var8502796096475496453&#45;&#45;factor8</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M456.63,-191.51C428.57,-162.63 369.54,-101.9 359.6,-91.67\"/>\n",
       "</g>\n",
       "<!-- var8502796096475496454 -->\n",
       "<g id=\"node8\" class=\"node\">\n",
       "<title>var8502796096475496454</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"543.44\" cy=\"-206.75\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"543.44\" y=\"-201.7\" font-family=\"Times,serif\" font-size=\"14.00\">v6</text>\n",
       "</g>\n",
       "<!-- var8502796096475496454&#45;&#45;factor8 -->\n",
       "<g id=\"edge32\" class=\"edge\">\n",
       "<title>var8502796096475496454&#45;&#45;factor8</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M523.61,-194.28C479.18,-166.32 374.37,-100.38 359.85,-91.25\"/>\n",
       "</g>\n",
       "<!-- var8646911284551352320 -->\n",
       "<g id=\"node9\" class=\"node\">\n",
       "<title>var8646911284551352320</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"111.44\" cy=\"-134.75\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"111.44\" y=\"-129.7\" font-family=\"Times,serif\" font-size=\"14.00\">x0</text>\n",
       "</g>\n",
       "<!-- factor0 -->\n",
       "<g id=\"node16\" class=\"node\">\n",
       "<title>factor0</title>\n",
       "<ellipse fill=\"black\" stroke=\"black\" cx=\"5.77\" cy=\"-63.09\" rx=\"1.8\" ry=\"1.8\"/>\n",
       "</g>\n",
       "<!-- var8646911284551352320&#45;&#45;factor0 -->\n",
       "<g id=\"edge1\" class=\"edge\">\n",
       "<title>var8646911284551352320&#45;&#45;factor0</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M92.55,-121.94C65,-103.26 16.11,-70.1 7.19,-64.05\"/>\n",
       "</g>\n",
       "<!-- var8646911284551352320&#45;&#45;factor3 -->\n",
       "<g id=\"edge4\" class=\"edge\">\n",
       "<title>var8646911284551352320&#45;&#45;factor3</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M121.7,-117.75C125.85,-110.86 130.11,-103.81 131.98,-100.71\"/>\n",
       "</g>\n",
       "<!-- var8646911284551352321 -->\n",
       "<g id=\"node10\" class=\"node\">\n",
       "<title>var8646911284551352321</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"183.44\" cy=\"-134.75\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"183.44\" y=\"-129.7\" font-family=\"Times,serif\" font-size=\"14.00\">x1</text>\n",
       "</g>\n",
       "<!-- var8646911284551352321&#45;&#45;factor3 -->\n",
       "<g id=\"edge6\" class=\"edge\">\n",
       "<title>var8646911284551352321&#45;&#45;factor3</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M164.64,-121.47C152.98,-113.23 139.36,-103.6 134.66,-100.28\"/>\n",
       "</g>\n",
       "<!-- var8646911284551352321&#45;&#45;factor4 -->\n",
       "<g id=\"edge9\" class=\"edge\">\n",
       "<title>var8646911284551352321&#45;&#45;factor4</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M181.67,-153.13C181.47,-155.18 181.29,-157 181.17,-158.31\"/>\n",
       "</g>\n",
       "<!-- var8646911284551352322 -->\n",
       "<g id=\"node11\" class=\"node\">\n",
       "<title>var8646911284551352322</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"255.44\" cy=\"-134.75\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"255.44\" y=\"-129.7\" font-family=\"Times,serif\" font-size=\"14.00\">x2</text>\n",
       "</g>\n",
       "<!-- var8646911284551352322&#45;&#45;factor4 -->\n",
       "<g id=\"edge11\" class=\"edge\">\n",
       "<title>var8646911284551352322&#45;&#45;factor4</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M231.05,-143.15C212.73,-149.47 189.67,-157.41 182.91,-159.74\"/>\n",
       "</g>\n",
       "<!-- var8646911284551352322&#45;&#45;factor5 -->\n",
       "<g id=\"edge14\" class=\"edge\">\n",
       "<title>var8646911284551352322&#45;&#45;factor5</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M239.16,-149.37C236.43,-151.83 234.02,-154 232.59,-155.28\"/>\n",
       "</g>\n",
       "<!-- var8646911284551352323 -->\n",
       "<g id=\"node12\" class=\"node\">\n",
       "<title>var8646911284551352323</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"327.44\" cy=\"-134.75\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"327.44\" y=\"-129.7\" font-family=\"Times,serif\" font-size=\"14.00\">x3</text>\n",
       "</g>\n",
       "<!-- var8646911284551352323&#45;&#45;factor5 -->\n",
       "<g id=\"edge16\" class=\"edge\">\n",
       "<title>var8646911284551352323&#45;&#45;factor5</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M301.41,-140.63C276.7,-146.22 242.08,-154.04 233.32,-156.02\"/>\n",
       "</g>\n",
       "<!-- var8646911284551352323&#45;&#45;factor6 -->\n",
       "<g id=\"edge19\" class=\"edge\">\n",
       "<title>var8646911284551352323&#45;&#45;factor6</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M307.29,-122.3C300.03,-117.81 292.85,-113.37 289.71,-111.43\"/>\n",
       "</g>\n",
       "<!-- var8646911284551352324 -->\n",
       "<g id=\"node13\" class=\"node\">\n",
       "<title>var8646911284551352324</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"399.44\" cy=\"-134.75\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"399.44\" y=\"-129.7\" font-family=\"Times,serif\" font-size=\"14.00\">x4</text>\n",
       "</g>\n",
       "<!-- var8646911284551352324&#45;&#45;factor6 -->\n",
       "<g id=\"edge21\" class=\"edge\">\n",
       "<title>var8646911284551352324&#45;&#45;factor6</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M373.33,-129.05C344.01,-122.64 299.29,-112.87 289.85,-110.81\"/>\n",
       "</g>\n",
       "<!-- var8646911284551352324&#45;&#45;factor7 -->\n",
       "<g id=\"edge24\" class=\"edge\">\n",
       "<title>var8646911284551352324&#45;&#45;factor7</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M379.93,-122.01C361.02,-109.66 334.18,-92.13 327,-87.45\"/>\n",
       "</g>\n",
       "<!-- var8646911284551352325 -->\n",
       "<g id=\"node14\" class=\"node\">\n",
       "<title>var8646911284551352325</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"471.44\" cy=\"-134.75\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"471.44\" y=\"-129.7\" font-family=\"Times,serif\" font-size=\"14.00\">x5</text>\n",
       "</g>\n",
       "<!-- var8646911284551352325&#45;&#45;factor7 -->\n",
       "<g id=\"edge26\" class=\"edge\">\n",
       "<title>var8646911284551352325&#45;&#45;factor7</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M446.96,-126.63C409.23,-114.12 339.87,-91.13 327.19,-86.93\"/>\n",
       "</g>\n",
       "<!-- var8646911284551352325&#45;&#45;factor8 -->\n",
       "<g id=\"edge29\" class=\"edge\">\n",
       "<title>var8646911284551352325&#45;&#45;factor8</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M447.97,-125.52C418.42,-113.89 370.27,-94.96 360.1,-90.95\"/>\n",
       "</g>\n",
       "<!-- var8646911284551352326 -->\n",
       "<g id=\"node15\" class=\"node\">\n",
       "<title>var8646911284551352326</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"543.44\" cy=\"-134.75\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"543.44\" y=\"-129.7\" font-family=\"Times,serif\" font-size=\"14.00\">x6</text>\n",
       "</g>\n",
       "<!-- var8646911284551352326&#45;&#45;factor8 -->\n",
       "<g id=\"edge31\" class=\"edge\">\n",
       "<title>var8646911284551352326&#45;&#45;factor8</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M517.68,-128.56C470.75,-117.27 374.79,-94.19 360.08,-90.65\"/>\n",
       "</g>\n",
       "</g>\n",
       "</svg>\n"
      ],
      "text/plain": [
       "<gtbook.display.show at 0x117392670>"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#| caption: A (somewhat messy) factor graph for the IMU preintegration problem.\n",
    "#| label: fig:imu_graph\n",
    "position_hints = {'u': 2, 'x': 1, 'z': 0}\n",
    "position_hints.update({'v':2})\n",
    "show(imu_graph, gtsam.Values(), hints=position_hints)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f7dd5960",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "We can optimize even *without* any external measurements, but unless you pay for a very expensive IMU, your estimate will drift very rapidly."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "6d75b7d4",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "initial error: 20190512.217341486\n",
      "final error: 1.4256262246403582e-23\n"
     ]
    }
   ],
   "source": [
    "print(f\"initial error: {imu_graph.error(values)}\")\n",
    "try:\n",
    "    lm = gtsam.LevenbergMarquardtOptimizer(imu_graph, values)\n",
    "    imu_result = lm.optimize()\n",
    "except RuntimeError:\n",
    "    print(\"LM Optimization failed\")\n",
    "print(f\"final error: {imu_graph.error(imu_result)}\")\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "ed87f830",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbEAAAGgCAYAAAAkWlmPAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/NK7nSAAAACXBIWXMAAA9hAAAPYQGoP6dpAADWG0lEQVR4nOydd3wb9fnHP9qWPGTLOx5xPOI4y9mJHUZaZoCwwggtuwXaUkp/FEoZhdJCy2qBFkJLyyqlhbJ3gDIDYSSxZTveey9JXtrj7veH+V5Osrbu5HO49+uVF2Ccu5N0+n7ueb7P83kkNE3TEBERERERWYBI5/sCREREREREokUUMRERERGRBYsoYiIiIiIiCxZRxEREREREFiyiiImIiIiILFhEERMRERERWbCIIiYiIiIismARRUxEREREZMEiipiIiIiIyIJFFDER3nnqqacgkUjQ09Mz35cieH7zm99AIpF4/Wz37t146qmn5ueCREQEjihiIrxz6qmn4osvvkBubu58X8qCRBQxEZHAyOf7AkSOfDIzM5GZmTnflxE2VqsVGo1mvi9DREQkDMRITIR3fNOJ27Ztw8qVK/HFF1+guroaarUaRUVFePLJJwEAb731FtatWweNRoNVq1Zhz549Xse79NJLUVRUNOc8/lJxobj00kuRlJSEhoYGnHjiiUhOTsZxxx0HAHA6nbjzzjuxbNkyqFQqZGZm4rLLLsP4+LjXMT788ENs27YN6enpUKvVKCwsxM6dO2G1WgEAH3/8MSQSCT7++GOvv9fT0wOJRBI0yioqKkJjYyM++eQTSCQSSCQSv69dROTbihiJicwLIyMjuOyyy/DLX/4S+fn5+Mtf/oLLL78c/f39ePHFF3HzzTdDq9Xit7/9Lc4880x0dXVh0aJFvFyL0+nE6aefjquuugq/+tWv4Ha7QVEUzjjjDOzduxe//OUvUV1djd7eXtx+++3Ytm0bDhw4ALVajZ6eHpx66qk4+uij8cQTTyA1NRWDg4PYs2cPnE5nzBHdK6+8gnPOOQdarRa7d+8GAKhUKi5etojIEYEoYiLzgtFoxLvvvov169cDADZs2ICsrCzcfffd6OjoYARr0aJFWLNmDV566SVcc801vFyLy+XCbbfdhssuu4z52XPPPYc9e/bgpZdewtlnn838vLKyEhs3bsRTTz2FH//4xzh48CDsdjvuu+8+VFZWMr/3ve99j5NrW7t2LdRqNVJSUrBlyxZOjikiciQhphNF5oXc3FxGwABAp9MhKysLa9as8Yq4KioqAAC9vb28Xs/OnTu9/vvNN99EamoqduzYAbfbzfxZs2YNcnJymNTgmjVroFQqceWVV+Lpp59GV1cXr9cpIiLijShiIvOCTqeb8zOlUjnn50qlEgBgt9t5uxaNRoOUlBSvn42OjmJychJKpRIKhcLrz8jICAwGAwCgpKQE//vf/5CVlYWrr74aJSUlKCkpwUMPPcTb9YqIiBxGTCeKLDgSEhLgcDjm/JwIS6T4KwbJyMhAenr6nKISQnJyMvPvRx99NI4++mh4PB4cOHAAf/nLX/Dzn/8c2dnZ2LVrFxISEgBgzjVHe70iIiKHESMxkQVHUVERxsbGMDo6yvzM6XTi3Xff5ewcp512GoxGIzweDzZs2DDnT3l5+Zy/I5PJsHnzZjzyyCMAgJqaGuZ6AaC+vt7r919//fWwrkWlUsFms8XwakREjlzESExkwXH++efjtttuw65du3DDDTfAbrfjz3/+MzweD2fn2LVrF5599lmccsopuPbaa7Fp0yYoFAoMDAzgo48+whlnnIGzzjoLf/3rX/Hhhx/i1FNPRWFhIex2O5544gkAwPHHHw8AyMnJwfHHH48//OEPSEtLw+LFi/HBBx/g5ZdfDutaVq1aheeeew7PP/88iouLkZCQgFWrVnH2WkVEFjKiiIksOJYsWYLXXnsNN998M8455xzk5ubiuuuuw/j4OO644w5OziGTyfD666/joYcewjPPPIM//OEPkMvlyM/Px7HHHsuIyJo1a/Dee+/h9ttvx8jICJKSkrBy5Uq8/vrrOPHEE5njPfPMM7jmmmtw4403wuPxYMeOHfjPf/6DDRs2hLyWO+64A8PDw7jiiiswMzODxYsXixZeIiLfIKFpmp7vixAREREREYkGcU9MRERERGTBIqYTRY5IKIoCRVFBf0cuF29/EZGFjhiJiRyRXH755XP6u3z/iIiILHzEPTGRI5Kenp6QfVjhFFWIiIgIG1HEREREREQWLGI6UURERERkwSKKmIiIiIjIgkUUMRERERGRBYsoYiIiIiIiCxZRxEREREREFiyiiImIiIiILFhEERMRERERWbCIIiYiIiIismARRUxEREREZMEiipiIiIiIyIJFFDERERERkQWLKGIiIiIiIgsWUcRERERERBYsooiJiIiIiCxYRBETEREREVmwiCImIiIiIrJgEUVMRERERGTBIoqYiIiIiMiCRRQxEREREZEFiyhiIiIiIiILFlHEREREREQWLKKIiYiIiIgsWEQRExERERFZsIgiJiIiIiKyYBFFTERERERkwSKKmIiIiIjIgkUUMRERERGRBYsoYiIiIiIiCxZRxEREREREFiyiiImIiIiILFjk830BIt8uaJoGRVGw2+2Qy+WQyWSQyWSQSCTzfWkiIiILEAlN0/R8X4TItwOapuFyueDxeGCz2QAAEokEUqkUcrlcFDUREZGIEUVMJC5QFMUIGAC4XC5GqCiKAvs2FEVNREQkXEQRE+EVmqbh8XjQ0dEBjUaDrKwsJiLzFSaappk/FEUx/18qlUKhUEAmk0Eul0MqlYqiJiIiAkDcExPhEXb6cGJiAsBs+pA8N9E07SVGEomE+W+ZTOYlana7nfkdImokUhNFTUTk24soYiK84PF44HK5QFEUpFIppNLDhbAURWF6ehqJiYmQywPfgqKoiYiIhEIUMRFOoWkabrcbbrcbNE0zgkIiMKvVirq6OkxPT0MikUCr1SI1NRU6nQ7JycleYudLuKJGURTkcjnUarUoaiIiRzjinpgIZ1AUBbfbzRRvsEWnvr4eFEXBYDAgNzcXRUVFcDgcmJycxOTkJCYmJkBRFLRaLdLS0pCWloakpKSgouYLEbSWlhYoFAoUFRUxUSC7UEQUNRGRIwcxEhOJGVKI4XK5mH0utkh4PB5MTk7CbrejsrISmZmZcDqdSExMRFJSEvLz80HTNCwWCyYmJjAxMYHe3l7QNI3U1FQvUQsmPuS8UqmUKQIh1+ZwOGC320VRExE5whBFTCQm2MUbAOYI2MzMDOrq6uB2u5Gfn4/s7GxQFDXnOBKJBElJSUhKSkJBQQFomobZbGZErbu7GxKJBGlpaYywJSYmBhQff2JKIjWPxwOPxwOHwwGJRMIIHvmn798TERERLqKIiUQNu/fLN5qhaRoDAwNoaWlBUVER49ABICyBkEgkSE5ORnJyMgoLC0FRFGZmZjAxMQGj0YjOzk7IZDImSktNTYVGowkqQOxIjVwj2cMjJf9sUSPCJoqaiIhwEUVMJGJINON2u5nqQ/Yi73K50NjYiImJCaxbtw7p6elobGz0amiOVBSkUim0Wi20Wi2AwxWOExMTGB0dRXt7O+RyOdLS0uBwOKDRaOaU8PsSiaiRPjWSfhQREREGYmGHSESESh9OTk6irq4OiYmJWLVqFVQqFQCgqakJUqkUy5YtA03TcDqdzN/nAo/Hw4ja0NAQnE4nEhISvPbUEhISIn6tZE+NnZ4URU1ERDiIkZhI2FAUBafT6Tf6omka3d3d6OzsRGlpKYqKiuY0MvP5vMROLZIoKiMjAxMTExgcHERLSwsSEhK89tSIwAYiVKQG+LfIEkVNRCR+iCImEhKSPiTVh74C5nA4UF9fD6vVik2bNjEpPzZ8i5jvuaRSKXQ6HXQ6HQDA7XYzpfz9/f1oamqCRqPx2lNTKpUhj+tP1FwuF5xOJ/P/RVETEYkfooiJBMW398tXwAwGA+rr66HT6bB169aADhykCTle+AqmXC5HRkYGMjIyAMzu2xFR6+7uhsViQWJiopeoKRSKoOfwJ2qk2IVEar6iRqofRUREuEEUMRG/hOr9oigK7e3t6OvrQ0VFBfLy8kIWUQhp+1WhUCAzMxOZmZkAAKfTyYhaZ2cnrFYrkpOTGUFLTU0NapEFgNkvI7DfQ5vNhq6uLpSVlUGpVM6pfhQREYkOUcRE5sC2jgLmFm8Q6yiPx4OqqiokJSWFPGa804mRolQqkZWVhaysLACzKVLSo9be3g673c6IWlpaGrRarZdgBboO8js0TWNkZASlpaVe6UfSmM0uFBFFTUQkfEQRE/GC3fvFTpURRkZGcOjQISxatAjl5eUhF3I28YzEYj2XSqVCTk4OcnJyAAB2u50RtebmZjidTqSkpDCilpKSEvS9INdDhIr8jLzfbFETZ6mJiISPKGIiAGYXVIfDAYvFArVaPWfvy+PxoKWlBcPDw1i5ciWzuIeLVCoVVDoxUhISEpCbm4vc3FzGcJiI2tDQENxu9xxR81fQ4VuxyRY04HAFKHETEUVNRCQ4ooiJMBV2BoMBTU1NOOaYY/xaR8nlcmzduhVqtTrq88QDvlOXEokEarUaarUaixYtYtz5yZ7awMAAPB4Ps5em0+nCqnwEIIqaiEiEiCL2LYfd+0XGm5BFkW0dtXjxYpSWlkZdLi60wg4ukUgkSExMRGJiIvLy8rzMjCcnJ9HX18e89v7+fuh0urDMjAFvUSN/HA4H0ywuiprItx1RxL6l+Ov9Yqf8/FlHxYKviPEtavMpmP7MjE0mE+rq6jA5OYmenh5IJBIvN5FgZsbkmIFmqTkcDq9ITRwQKvJtQhSxbyG+1lHswZUURXlZR1VXV4d0tgiHeC6kQlu0SaQGAKtXr/Zy6Pc1MybCRsyMgx1TnHotIiKK2LcOEn35s46SSCTweDzYv3+/X+uoWDiS04mRIpVKkZKSgpSUFCxevNjLzHh8fBwdHR2MmTERNrVazYmo+Tr0i6ImstARRexbgm/vlz/rqKamJlAUhS1btiA1NZXT88dbxIQmmOR6/AmGVCplikAAbzPj4eFhtLa2QqlUMqIWjplxIFGjKEoUNZEjClHEvgWQXiRi++TbvEyso1JSUgCAcwEj5xRys7OQYJsZA7OiNjU1NcfMmL2nFq6ZMTk+W9S+/vprFBQUID09XZx6LbLgEEXsCIbdTOsvfUhRFDo6OtDb24uKigqkpaVh7969vFzLtz0SA6IXV5lM5tfMeHJyMmYzY3I8cm/4Tr0me2ri1GsRoSKK2BFKoOINgs1mQ11dHdxuN2MdRdJMoYZJxnJN31a4fO3BzIx7enpgNpsjMjMmDzjkD7le9vBTceq1iFARRewIhG0d5S8dRKyjcnNzsWzZMqYXid0fxvXCJKYT+SNcM2OSfvQ1M/b3eUcy9do3/SgiEk9EETuCYD85+0sfhrKOIr9L/i6XiOnE+ImrPzNjImq+Zsapqalhfd6RiJo49VoknogidoQQKn1oNpuh1+shl8tRXV0NjUYz5xjsBYpr/IkYX2lLITKfoqpSqZCdnY3s7GwA3mbGra2tcDqdaG1tRUZGRlhmxoA49VpEOIgidgQQrPeLpmkMDg6iubk5pHUUOxLjGjGdKBzYZsYA8OmnnyIjIwM2my0iM2M2oaZeA6KoifCDKGILGHbvF7GOYi/gxDrKZDJh7dq1TCFAIPiOxOLJtzmdGA1ZWVlISkoCTdOw2WxMpEbMjLVaLSNqycnJUYmaOPVahA9EEVugUBQFt9sdMH04NTUFvV4PjUaDrVu3hmUdFa9IjKIoDAwMQKFQQKfTBa2cO1IQoqgS2GldiUQCjUYDjUbjZWZM9tSImTG7Ry2UmTE5bqCp12SWGlvUxKnXIuEiitgCg/3lJ4uPb/qwp6cHHR0dKCkpwZIlSyJaCPhK+5HjktJ+kmJqbGz0mpicmpoa0aBNkdgJtjfJNjPOz8/38n2cnJxEd3d3xGbG5LihRI08mCmVSlHURAIiitgCgiwgFEVBqVTOETCHw4GGhgZYLBZs3LgxKucNvoZXSiQSuFwu7Nu3D9nZ2SgpKWF+RlJXLS0tcDqdTOpKp9OFlbrydy4hItTriqTARiKRIDk5GcnJySgsLARFUV5mxl1dXYyNFhG1UGbG5Lj+RO3rr79GcXExdDqdOHZGxC+iiC0QyFNqW1sbFAoFysvLvf6/0WhEfX090tLSUF1dHXWKjo9IjKQPHQ4HVq9ejdzcXKYQRaVSIScnBzk5OX73YyiKYhZEnU4X1lM+ILz0ndCuh0AKMKIVA39mxjMzM37NjMnnGMrMGDgsahRFMaIlDggV8YcoYgLHt/eLfJkJbOuoZcuWIT8/P6YvMhnHwhV2ux11dXWw2+1QqVTMJORA5/bdjyFP+SaTCV1dXYyvoE6nYxZEf8cRiQyu3jOpVAqtVgutVouioiJQFMX4Po6OjqKtrc3LzJg49AeCpmlGnMSp1yL+EEVMwPjr/ZJKpcx/s62jtmzZguTk5JjPyWU6kRgLZ2Zmori4GI2NjQDmGhAHwl/qanp6GiaTiXF3V6lUXqIWyjNwPhHighrMXZ8LpFJpSDNj8hn6MzOmKMqvmwgAv6ImTr3+9iGKmEAhT5m+vV9SqRQulwujo6M4dOgQcnJyvKyjYoWLdCJN0+jo6EBPTw8qKiqQn5+PiYmJmI/rO7LE7XZjamoKJpMJvb29aGxsRGJiIuRyOdN+wLZXmk+EnE4E4iew/syMiaj5mhmnpqYy1mnBYIuav6nXbFFjmxmLDv1HBsL4hoswkPQhqT7090UzmUwYGRnxax0VK7GmEx0OB+rq6uBwOLyiQz4WC7lcjvT0dKSnpwM47BnY19cHs9mMvXv3ejXtarVasbnWh3iLmC++nyExM56cnERvby8oikJDQwPS09PDMjMGxKnX3zZEERMQoXq/zGYzU+wQyDoqVmJJJxqNRtTV1SE9PR3r1q3zioK43mvzB/EMdDgcUKlUKC0tZYpEiBMFu2ouOTk5rouWEBdI9ow5IcA2M6ZpGh999BEKCwthNpvR1dUFi8WCpKQkrz21UNF2uKImDghdmIgiJgDC6f0i1lFarRZyuZwXAQOiExuaptHZ2Ynu7u6AxSXz4dihVquhVquZYhKLxcKIWk9PDyQSCbMY6nS6sKrmYrkeITLfkVgwyLVlZmYiLy8PQGAzY7ZDf7i+j4A49fpIQBSxeYZtHQXMLXpwu91obGyE0WjE2rVrYbVaMT4+ztv1RBqJORwO1NfXw2azYfPmzcx0aF9899r4XAz8HZvdtFtQUOBVCj42NoaOjg4oFAovUQvH5STW65pvhCxi/qJEf2bGRNRaW1vhcDiYFHJqaiq0Wm1MouZwOGC325miKlHUhIcoYvMIe+oygDn7NVNTU6irq4NarWaso0g6kS8iKewwmUyoq6tDWloa1q5dGzSt43vc+Y5MfEvB2VVzAwMDaG5uZgoMdDpdWHsxCxF/kb9QCPS9YJOQkMD0GQLw6jMcHh6Gy+WK2syYwB4QKk69Fh6iiM0DoeZ+BbOOkkqlvItYqOPTNI3u7m50dnaivLwcBQUFYTWvCnmeGLtqrqSkxGtaMnuwJCnlD+cJP5briRdCHodD3rNIinF8U8h8mRmTY3/yySeorq5mKiPFqdfxRxSxOBNq7pfT6URDQwNmZmawYcMGpr+GEA8RC7bgkuszm83YtGkTtFotJ8cVGr7Tkh0OB9N03dzczDzhE1ELdzEUGkIWMdIjFu31+Wuet1qtjKj5mhmnpqaGVexDrok88JEMhDj1en4QRSyOBOr9IhDrqNTUVGzdutVv+opvEQu2JzYxMYG6ujpotdqorK0W8jwxf/ZYJpPJazFkN+z62mMJVcCFLmJcLvwSiQSJiYlITEz0MjMmEXdPTw8AhG1mTBx0yDWKU6/nB1HE4kCo3i+KotDZ2Ymenp6Q1lF8RzT+0ons9GZZWRkWL14c8cJHfj9eiybf7xF5wmcvhiaTCUajEZ2dnYxfINutQogIXcT4LgAijjAFBQWgaZop9iE2Z8HMjAM1YgdKP4pTr/lBFDGeCZU+tNlsqK+vh8vlCss6Kt6RmMvlQkNDA6anp6N2xgfiL2LxhL0YEhNcUiRC7LEUCgU8Hg9GR0cFZY8l5M+D60gsFBKJJCwzYyJq4QpPIFETp15zgyhiPEKir0DpQ2IdlZ2djYqKirAKBeK5JzY5OQm9Xo/k5GRUV1fHtPDGs29svhdlX79At9uN/v5+9PX1MfZYkTbs8oUoYoEJZGY8OTmJ0dFRTE1NAQCamprCMjMm+BO1YANCxanXwRFFjAd8e798Bczj8aC1tRVDQ0NYsWIFcnNzwz52PETM4/Ggp6cH7e3tKC0tRVFRUcxfIHYkFg+EtAcll8uh1WqhUqmwadMmOJ1Oprigra3Nq7dJp9OFVQbOFUIWMZJ6Fwrsh5MlS5ZgdHQUXV1dUKlUGBoaCmlmHAiyX0ZgixpJP/74xz/Geeedh7PPPpu317dQEUWMY3x7v3yrq8xmM+rq6iCVSqOyjuJbxAAw/TX+qiOjJd4iJjTYr1upVHo17JIycJPJhMHBQVAUBa1Wy1Q+JiUl8eokIlQRm+9ILBQ0TUOpVKKkpARAYDNj9p5aONkMf6LW2trKPBSLeCOKGEeQp6fu7m5oNBrodLqA1lGFhYUoKyuL6gvKpwfh1NQUxsfHoVQqsXXrVk73bXxF7EhOJ0ZKIHssk8mE7u5upriAPUONq9codBET6rUBsxkVttj4mhm73W6m8pE9ZYE9IDScCl+JRAKr1cqb1dxCRxQxDmAXb4yNjSEjI4O5kQFv66g1a9YwvUfRwOW8LwJN0+jv70draysSExORnp7OeeHBtzmdSAhnQQ5kj2UymbyGShJBCzdlFQihi5iQIzFSYh8IuVyOjIwMZGRkAJgtkpqYmMDk5CS6u7tx6NChsPZGyYMNF/MCj0REEYsRkj4k5ba+k5fZ1lHV1dVISEiI6XxExLhafNxuNw4dOoSJiQmsX78eIyMjMR/TH2I6MbrXzS4uWLJkCTweD/N0T1JW5Ome/ImkSETIIia0PTFfwpl1xkahUCArKwtZWVkA4LU32tHRAZvN5iVqxOwbACwWCxITE3l5HQDwm9/8BnfccYfXz7Kzs3lbD7hEFLEoCWQdRfasaJpGb28v2tvbUVxcjOLiYk4WC/KlCfUUGA7T09PQ6/WMwKpUKoyOjvKSrvQnYnwJmlAXZS6QyWRz5m+RhbCzsxM2mw3JycleC2Gw+0TIKbuFHomFwndvlLjCsM2MX331VchkMiQlJfHu3blixQr873//Y/6bq0G7fCOKWBQE6/2SSqVwOp2oqakJaB0VC1yIGE3TGBgYQEtLC5YsWYKSkhLm+vnccyPnZv+T7/MICT7Ewvfp3m63MwshscciXoE6nW6OrZKQIzEhCywQeSQWCrYrDDBb8NPT04N33nkHBoMBmzZtwpYtW/Dd734XZ555JiorKzk7NzCb/uR6yG48EEUsQkL1frlcLgwPDyMjIyOgdVQssEUsGtxuN5qammAwGLBu3TqvvTtyfCLOXCKmE+PzuhMSEpCbm4vc3Fy/XoHAYVslnU4naKFYCJEYn/19arUal112Gc444wwUFRXh4MGDOHDgAD766CPs3buXcxFrb2/HokWLoFKpsHnzZvz+979HcXExp+fgA1HEwoTd+xXMOmpsbAypqalYs2YNL4sDOWY0IjYzMwO9Xg+VShVwf47PSCxeJsBCXZTjjT+vQOJAYTAY0NnZydzHQ0ND0Ol0Me/ZconQRczj8cTFecVsNgMAVq9ejXXr1uHKK6/k/BybN2/GP//5TyxduhSjo6O48847UV1djcbGxjkPukJDFLEwoCgKbrc7pHWU0+lk3LL5WkhJ31mkYkBmZBUVFaGkpCTg4sBH9SMhnk72Qoz45ltc/dkqdXZ2Ynx8HENDQ2htbUVCQoJX5eN8zlATemEHF/vS4UDK6/l8L7Zv3878+6pVq1BVVYWSkhI8/fTTuO6663g7LxeIIhYEdud8oOGBY2NjaGhoQHZ2NtavX4++vj7MzMzwel2RNDx7PB40NTVhbGwsrPJ+PoWG/d65XC44HI5vTe+LEEVVKpVCo9EgMTERlZWVTF8T6U8jJeBE1FJTU+O62b8QIrF4XJ/ZbA7qps8HiYmJWLVqFdrb2+N2zmgRRSwAvsUbvgIWyDoqHo4a4Z7DbDZDr9dDoVBg69atYaWK4pFONJlM0Ov1cDqdTGM4mZ7MxR7DfEc8Cwl21sC3r8npdDLjZki1HHuGGt/2WELerwPiF4nxXV7vD4fDgebmZhx99NFxPW80iCLmB9/eL98vUjDrKL4KI9iEI2JDQ0NobGyM2B2E75TfwMAABgYGUFZWhvT0dExPT8NkMqG9vR12u51ZJEklXbSLpBAjHyEuyMFS30qlck61HBG1gYEBUBTlZanEtT0WRVHzms4MxZEUiV1//fXYsWMHCgsLMTY2hjvvvBPT09O45JJLeDsnV4gixiJQ7xebwcFBNDU1oaCgAEuXLp1zE/s2O/NBsGjJ4/GgubkZo6OjqKysZEqvw4WvPTGypzg8PIyNGzciOTkZTqfTa3oy20OQLJKkik6n03FqtxRvhCiqQGQl9mq1Gnl5ecy+r9lsZiofiT0Wu+k61s9L6OnEeO6J8R2JDQwM4IILLoDBYEBmZia2bNmCL7/8EosXL+b1vFwgitg3hJr7xS5ND7a3NJ/pRIvFAr1ez0SI4YyF8IWPdOLMzAxqa2sBzG4ap6am+o1WfT0EfWc5KRQKRtCCmakuVKGbD6ItQmLPUCssLARFUZiensbExARjj8V2dNfpdBFX8gm9sCPekRifPPfcc7wen09EEcPsE5XT6QwYfRHrqISEhJDWUfFKJ/o+2Q8PD6OxsRH5+fl+I8Rw4TqdSNKaRUVFGBgYCHvPy7eSjm23xJ7JRUTN15lCaJGPUJuKubouYlKcmpo6xx6rr6/Pyx4r3P1PMRKbxWKxICkpiffzLFS+1SJG0oek+tBXwKKxjop3JEZRFFpaWjA0NIRVq1YxFjaxHJsLAWBfF0lrDg4ORu1i72u3RIoOTCaTlzMFmR4gNBETKnyJazB7LLL/SeyxyEOIr2AJvbAjXpGY1WoVRSwI31oRC9X75XQ60dDQELF1VDz2xIiIWa1W6PV6SCSSqGaT+YOLdKLdbkdtbS1omva6Ll9xiWWBYhcdEGcKUnRgMplA0zQaGxuZRVIITbxCXJDjFSEGsscymUxobGyE2+32mqGWnJwsRmLfMB/ViQuJb52IhdP7ZTKZUFdXh9TU1Iito+IViU1MTKCxsRGLFi3CsmXLOPuyxxqJGY1G6PV6ZGdno6KiwutLzleExHamKCgowPj4ONrb26FWq5kmXrVa7dXEy6ddkD+EGhnOV5ozkD2WyWRCb28v83sqlQparRYajUZQDwFkHYmHiJnNZjESC8K3SsTY1lHA3N4vmqbR0dGBnp4elJeXo6CgIKq0F597YiQCm5ycxOrVqzk37Iw2EqNpGl1dXejq6kJFRQXy8/P9HjtetlMymYxJAbvdbmaB7Ojo4LSUf6EjhJSdrz0WRVEwm81MJmT//v2Qy+VeDyHzHVmzMzh8Y7FYIq4y/jbxrRExiqIwPj4OmUyGxMTEOTef3W5HXV0dnE4ntmzZEvUAOj4jMZvNBr1eD7fbjSVLlvDiOB2N0LhcLmbB2bRpE7RabcBjxwPf88jlckGU8s+3WPhDiAUnUqkUKSkpUCgUWLJkCXQ6HaampjAxMYHBwUG0tLRArVZ7lfPHu5+MfMePlBL7hcwRL2Ls3q/u7m5otdo5AkWso7KysrB+/fqYUk3seWJcLg7kGnNycqBSqXj70kaaTiQzyRITE1FVVRW0jFooBRe+pfxmsxkmkyniUv5IEMLr9oeQy9jJnphMJmM+DwBMZE360w4dOuQ1Qy0e9lhsJx++Eac6B+eIFjHf3i+5XO4VJVEUhdbWVgwMDGDFihVYtGhRzOckXx6u8uUURaGtrQ39/f3MNdbV1fFqDRXusUnjt+9MsmDHFpoBMLvfKZpS/oWOECMxQqDCDt/Imj1MsqWlBU6nc84MNa6Fmny/4yViYiQWmCNWxPz1frH3qywWC+rq6gAA1dXVnN0kXE5ettlsqKurg9vtRlVVFbO5y2fKMhyhoSgKzc3NGBkZwdq1axmvPS6OzQWxLCyRlPLrdLqIrJaEKBZCFrFwo0T2MEmappl0sa89FomsubBwild5PSD2iYXiiBOxYL1fRMRCWUfFAjmWx+OJKeU3Pj6O+vp6v1V+fIpYqHSizWZDbW0tU9YfqSuIUNNqgQhWyt/T0wOpVMosjsFK+YX6uoUsYtGU2EskEmg0Gmg0mjn2WEajEZ2dnZDJZMznReyxorm2eETk5J4TI7HAHFEiFso6SiKRYHx8HMPDw2GNJYkGcs5oRYaiKHR0dKC3txfLly9HXl7enN+Zr8GVRFhzcnJQUVER1QIjtHRiJPiW8hOrJZPJFFYpvxDFQugiFuu1BbLHMplMGB4eRmtrK2OPFckeaDwjMbPZLO6JBeGIETESfQWyjpqenkZ/fz8kEknYY0miJdpIiVRIulwur/QhV8cPB38CSdM0Ojs70d3dHVBYwz220NOJkcC2WgpVyk8mgguNhVDYwSXszwyYLRKZmppi+tPIHii7SMRfoVe8IjFA3BMLxYIXMd/eL3/WUX19fWhra4NOpwNN07z3mETjn2gwGFBfX4+MjIyQFZJ8+jP6phOdTifq6+thtVqxefNmpKSkRH3shR6JhSJYKb/BYABN06ivrxeUK79QIzGapuMisHK5fM4eKBkMyh4PRESN2GPFKxIj6URxTywwC1rEiPMGiRx8m5edTicOHTqE6elpbNiwARaLBcPDw7xfVyTWU+wG64qKCuTl5c2rPyM7EpuamkJtbS1SUlJQVVUVc1m/EBdLPmGX8vf392N8fBxarZbXUv5IEbKIAfFpJmajVCq97LHYRSJDQ0Nwu91MCT8RWj7fP6fTCbfbLaYTg7AgRYxtHRUofWgymVBfX4+UlBRUV1dDqVTCbrfz7jAPhC8yDocDdXV1cDgcETVY813YAQD9/f1oaWlBSUkJlixZwskXlc+9vIWAXC7H4sWLBVXKL1QRI/fJfKc6fXsKLRYLI2hWqxV79+71GjfDdXRtNpsBQIzEgrDgRCxU8QZ7/2bp0qUoLCycU53IN+Gk+4xGI5NaWrduXUQN1nyKGDluW1sb1q1bx6RZuCCejh1C3H9iw2cpfyQIXcSEdG0SiQRJSUlISkoCRVGYmZlBQUEBJiYmMDY2xkTXbFFTqVQxndNisTAVlyL+WVAiRqIvko/2vcHZ1lH+9m/iKWKBRIbtMRitPyNfEY3VakVNTQ0AYPPmzZw//S0EceGTYJ8zV6X8kSJ0EZvvSCwQHo8HMpkMWq0WWq0WRUVF8Hg8jD3WwMAAmpubodFomM8tNTU14pQ8KeoQ4mckFBaEiLGtowKlD8OxjoqXiAXaE3M6nairq4PNZoupSIKPSGxsbAz19fXIzc2F2WyO+QnSH0dadWIkRPK6Yy3lj/S6hPh+kfJ6IV4b4L86kW2PVVJSApfLxaSMOzs7YbVakZyczHxu4aSMyVRnob4PQkDwIhYqfRiJddR8RmITExPQ6/VITU2NuUiCq8GVgHdhyYoVK5CdnY3+/n7e+qyO5OpEvghWyt/Z2QmbzRa1K79QRUzIpf/A4UgsGAqFYo49FomuScqYfG5khprvaxbL60MjaBEL1fsVqXUUn6Xpgc5D0zS6u7vR2dk5Z48uluNzEYmxI0NSWMKeGM0HR5K4RApXYsGlK78QRrH4Q+gDMT0eT8TVpCqVymuGms1mY0Str68PNE17OfNrNBpYLJa4zlL7wx/+gJtvvhnXXnstHnzwwbicM1YEKWLs3i9f6yjC0NAQGhsbkZ+fj/Ly8rBueFIWy/cXhKQTyXRos9kcdERJpHAhYpOTk9Dr9dBqtaiurmbSUeR9XsiRmBAXZT5fdyyu/EKNxIQqroRYm53Z9lj5+flen5vRaERdXR2uueYaLFu2DDRNo7e3F4sXL+bwFcxl//79eOyxx7B69Wpez8M1ghMxiqLgdrsDpg/dbjeam5sxNjaGysrKiIbFsR3m+RQxqVQKi8WCffv2MSX+XI5OiUXEaJpGf38/WltbUVpaiqKiojnWXHxOYP42R2LxIFJXfqGKhdAjMa6vz/dzq6iowO7du/Hss8+iqakJpaWlWLx4MY477jj85je/QW5uLmfnBmb33r7//e/j73//O+68805Oj803ghExdu8XeTr0Zx1VV1cHlUoVlXUUETGPx8PbeHrSSzI9PY3y8nIsXryY80Ui2upEj8eDxsZGGI1GrF+/npnPxNXxQ/Ft3hObr4gnUCk/2ZdxOBzo6uqCxWLhtZQ/UhbCnhif16dSqbB9+3a0t7dDIpHgmWeewaeffooPPviAlz2yq6++GqeeeiqOP/54UcSiwbd4w1fA2NZR4c6u8geJ6vjaFyMTjs1mM7Kzs1FUVMTLeaKJxCwWC2pra6FQKFBVVRX0AYDLwhE2vhG10WiEVqvl3KlCCIuwUPEt5d+3bx+0Wi2mpqZ4LeWPFKFHYuEUdnABcbBPTk7GqaeeilNPPZXzczz33HOoqanB/v37OT92PJh3EaMoCmNjY3A4HMjOzp6zALGto4JFD+HCV4Xi1NQU9Ho9kpKSojbIDZdIRWZ0dBQNDQ3Iz88Pa/QM35GY2WxGTU0N3G43nE4nU3acnp6OlJQUQS9eRxLkYTErKwtpaWm8lvJHilDTnIR4GQCbzWZe3Tr6+/tx7bXX4r333pu3B5ZYmTcRY/d+GY1GWCwW5OTkeP3OxMQE6urqvKyjYiUSX8NwYEeJxKKJOJjzRbiRGEVRaG9vR19fH1atWjXn/Q12fL4iMbPZjC+//BIFBQUoLCxkemmMRiMaGhq8KuvS09OjmvUEiOnEcGGn7fgs5Y+UhRCJxeP6LBYLZwVh/jh48CDGxsawfv165mcejweffvopHn74YTgcDsFPMp8XEfNNH8rlcq/oKJh1VKxIpVLG8T5WXC4XDh06hMnJSWzYsAFpaWkAuBdKX8IRMbYvY7CxLv7gY++KpmlMTU1henqaKchxOp1zpvKSCq2xsTG0t7cjISHBq7IunEhAiGIhVIKJq28pv91uZ6yxIi3ljxShi1i8IjGr1cprZue4445DQ0OD188uu+wyLFu2DDfeeKPgBQyYJxEjNzpJZ7BFzG63o76+Hna7PebRH/7gSmCmpqZQV1cHjUaDrVu3ekWJfHobhnN80lidlpYWsS8jwH060eVyob6+nom2c3Jy/B7ft0LL7XYzYzFIJMD2E0xOTl5QgiXEa40kQkxISMCiRYuiKuWP5rqELGLxjMT49E1MTk7GypUrvX6WmJiI9PT0OT8XKvOWTmSnrGQyGdxuN2MdlZmZGdXiGw6x7omxS9SLi4tRXFw8ZxHgu6maREq+CxDpJ2lvb0dZWVnUlZFcRmIzMzOora1FYmJiUDcVf8jlcmRkZCAjIwMAmOZQk8mEvr4+SCQSZtH0NVsVYjpRiESb5vRXyu9vuGS0rvxiJDaLxWIRHexDMO+FHcDsom+1WlFXVxfT5OBwiEXE3G43Dh06hImJiaAO7/FIJwLeCxD72tipzWiPz8WiOzIygoaGBhQVFaG0tBRtbW0xpXLVajXy8vKQl5fHuIgbjUYMDg6iubmZeYJUq9WCFQ2hwdVeHds3EJhbyh+pK7+QCzviYZhAmA8R+/jjj+N6vliZdxGzWCzM4nbUUUfx7hMWrYiRiEKtVqO6ujqoQS7fkRj58pAvktlshl6vh1KpDHlt4RBrOpGmabS1taG/vx+rV69GdnY2c1wiLlxYbxEH8eLiYrhcLkxMTMBoNGJkZAQURUGv1yM9PR06nS6u1j2BmO/z+4OvgpNwXPnZ+2m+lXFCjsTIdyOeJfYigZk3EZNIJBgaGkJTUxMyMzPhcDji8mFFKmI0TTNP+uH2qMUrEqMoiol2CgsLUVZWxskXP5ZIjPgx2u12bNmyxesp0jdNyeXiqVAomIm8ZrMZ+/fvR3p6OoxGIzo7O5n9mvT0dKSlpXHqoBIOQo0M41E1GcyVf3h42G8pv5BFjO0mxCfEOEFMJwZn3kSsp6cHra2tWL16NRITEzE6OhqX80ZSdOF2u9HU1ASDwYC1a9cyezNcniMayJenra0NIyMjXtEOF0S7JzY9PY3a2lokJyejqqpqzp5mPL0TJRIJCgoKUFBQ4LVf093djcbGRqY3TafTxa03TWiRmL991XjALuUH4LeUX6lUQqVSYWpqitdS/miI56wzs9kc9sT3byvzJmKLFi1CRkYGEhISYLfbQVFUXJ6+wo3ESIpOoVCguro6okZAvkXM4XAAAEwmE6qqqjiPYKNJJxJD5mDRajwnO7Nh79eUlpYyIzFMJhMaGhoY93B2qfi3Aa5Su7Hir5SfWGLV19fzWsofDYGG8vKB1WoVI7EQzJuIqVQqRrDIE3s8ylbDEbHBwUE0NTVh8eLFKC0tjfia+NwTM5lMzPiZyspKXlKwkaQTKYpCW1sbBgYGQhoyC8U70XckxszMDEwmE0ZHR9HW1saktnQ6HVJTUzmpkhViOlEoIuZLQkIC1Go1UlJSUFxczGspfzTEqzKRpBPFPbHgzHthB+BtzMv3XoVMJmMiGV88Hg+ampowNjaGNWvWME+G0ZyD60iMpmn09PSgo6MD5eXlaG1t5U3ww43EnE4n9Ho9nE5nWBGhEF3sJRIJUlJSkJKSgqKiIq/UVnt7O+x2O7RaLVMgEotBrtDEQqgiBhwu7OC7lD8a4tUjZrPZQFGUmE4MgSBETCKRcOqkEYxAkRhJH8rl8qgc8tlwnU50u91oaGjA1NQUNm7ciNTUVHR0dPCWsgxHbKamplBbWwutVht2T99CmCfmm9oiVXUmk8nLIJcUiYQbBQhNvAFhi1igZmeuS/mjIZ49YgCOiHQisefTaDTMd2t8fBwGgwFFRUUxpfDntTqRja/1FF/4S/WR/RyuKvyIiHGxaT4zMwO9Xo+EhAQv/0g+991CpRNJupV4RUbyGoWQTowE9uBCUlVnNBoxMDCA5uZmryggNTVVUAUIoQgpYhQFuN1AnNN1s6cOb3881lL+aIiXg73FYoFUKl2wxrxsxsfHceONN+KSSy7B9u3b0dPTg9NOOw1NTU2oqKjA22+/HfXQT0FEYsBh1454nIeImMfjQUtLC0ZGRiIesBnqHEDsT2xEXEmzMHux4VPEAkVMFEWhtbUVQ0NDUaVb+TIWjhfsqrqSkhI4nU6mN62pqQlut9trwfTtTRNaxBNMxJT33APVXXfBecUVcPzxj/G+tKianaMp5Y9mvzOellOJiYmCu28igTzIDw8PY+/evfjrX/8Kh8OBZ555BlKpFPX19fjNb36D3/zmN3jyySejOodgRCxekRjZr7JYLNDr9ZBKpaiurua0Io3dxxWNiFEUhZaWFgwNDQUUV75FzPfYDocDer0ebrcbVVVVUfu5CT2dGAlKpRLZ2dnIzs5mNuFNJhMMBgM6OzuhVCoZQYvX03skBI3EvmmYl5jN8bwkBi4qlcMp5Y/GlT9ePWxkDMtCFjHC5OQk5HI5UlNTMTg4iK+//ho/+tGPsHLlSpx22ml48MEHoz62YNKJfM358kUmk8Fut+OLL75AXl4eysvLOb8hyfGiKVSx2+3Q6/XweDyorq4OKBZ87i/5RkyTk5Oora1FWloa1q9fH3W1XrwLO+LZAyWRSJCUlISkpCQUFhbC4/Ew5sXd3d2wWCxQKpWQy+VMb9p8L06BJqgDAE32YWZm4nxVs/AhFJG48qelpQV0eYlnOvFIqUyUy+XQaDT46KOPGH9XIlwTExMxpUwFE4nFI51IURQGBwdhs9mwdu1aThuE2ZCqqkgjJaPRiLq6OmRmZmL58uVBvyjxisTI3k9paSmKiopiWngXQmEHV8hkMqSnpzP+mk1NTXA6nbBYLOjv7wcAL/Pi+dj3CCby9DcVcfMVicXDxT5aV/54jmERgl0aF1RUVODEE0/EZZddhsTERBx33HEoKSlh3us1a9ZEfex5FTH2osZ3OtFqtTIRDkkD8UkkIkPTNLq7u9HZ2Ylly5ahoKCA0+NHChGxxsZGjIyMBDU7jvS4vrZTC3mPLBJkMhmSkpJQWlrK9KYZjcY5ezXp6elITU2NWx9SSBE7giKxYERSyu9wOOIiLHxPdY4H5Duu0+lw9913Y/369dDpdDjppJMAzGZ58vPzccwxx0R9DkFFYnyJ2OjoKBoaGpinroMHD/JyHjbhNjy7XC40NDRgZmYGmzZtCnuKK58iRlEUxsbGoFKpON0vjPcTpVCnKbN705YsWcLs1RiNRrS2tsLpdDJl4unp6bxt7gd9f8jiOY97YvP52QUr5TcYDPB4PHA4HLyW8i9030TixKTRaOByuSCTyXDxxRcDmE3JUhSFsrIy3HTTTTGdRzAiJpfLOU8nkmq6wcFBrFy5Ejk5OTCbzXEtIAkGccbXaDSoqqqKyHmALxGbmJjAyMgIVCoVNm/ezGlE8G2KvPwRzgRlmqaZuWlGoxE9PT1eC6pOp+PMoSJoJPbN4rmQCzu4hF3KT+7j5ORkr/5Brkv5F/qe2O7du2EymXDnnXfizTffxD//+U8UFhZCo9Ew+8fJycmgaRpVVVVYvnx5VOcRTDqR60jMZrNBr9eDoigvNwkiLnw/pYcSGdJrFa4zfqTHjxT2sE+tVguNRsN5SuvbtCfmS7ivWyKRzOlNI2mtvr4+NDU1eZkXa7XaqBd7Ie+JCU3E2NA0jYSEBOTn53v1D3Jdyr/QRaygoICpDJ2YmGACiqmpKVitVsaRZHJyEnfdddfCFDE2weygIoVMiM7JycGyZcu8FmOuerhCESidSFEUmpubMTIyEpO1FZeCQFEUY7e1fv16GI1Gzj4LNvNRnbjQIU/4aWlpTG8aqahrbGyEx+OZ05sWLmGlE2dmAJoG5iEVLFQR8+0T46uUfyGnE2maxrnnngtgNhV7+eWX4/LLL+flXIIRMS4KOyiKQnt7O/r6+rBixQosWrRozu+wfRr5FDF/6UQSHdI0HfNeE1eRmN1uR21tLQCgqqoKarUaJpOJl1TltzmdyFXk7+tQwa6oa29vh0qlYnweQ0UAYUViFAVYrUCcIwIhR2KhHoC5KuU3m82cFFTNB6Q4jKZpKJVK7N+/H/n5+cjNzYXb7QZN05DL5ZBIJJiamoJGo4naN3fe04mEWEvsSX8VacYN9ATD7uHiE1+RMRgMqKurQ3Z2NioqKmIWUC5EzGQyQa/Xzynp51Nsvq3pRD7wV1EXKAJIT09HcnKy13sTtHgiMRG0RAIJTUNiNoOeBxET6ucYqWNHtKX8VqsVhYWFvLyGRx99FI8++ih6enoAACtWrMBtt92G7du3c3YOdjbqkUcegdvtxn333Yfc3Fzmd3p6erBr1y7ceOONOOuss6I6j2AisVj2xMbHx1FfX4+srKyQ/VUSiSQujdXkA6RpGp2dneju7kZFRQXy8/M5O360IkbMONva2lBeXo6CgoI5llZ8iM23PZ3I96Isk8mQkZHBDG8lBSImkwn9/f2QSCRMBJCenh48OpRIZlOKMzOzfzhuSfly8EuMWkZxxtIz5vw/MqxzoUZiwQinlL+rqwt6vR4zMzO8jZnJz8/H3XffjdLSUgDA008/jTPOOAO1tbVYsWIFZ+ch79NFF12En/3sZ7jssstw9913Y+nSpXjiiSdw++23Y9OmTdi4cWPU5xCMiEWTTqQoCh0dHejt7cXy5cuRl5cX1t+Lh4jJZDK4XC7U1NTAYrFg8+bNSElJ4ez40TRTA7NPkY2NjTAajdiwYQPS0tI4O3Yovu3pxHijVquRl5eHvLw8UBTF9KYNDQ2htbUVSqUSHo8HRqPRb28anZQEyczMbCTG4XW93fk2Ln3zUlCgkJuUi02LNnn9/3hOTo4GLr0T/ZXyUxSFDz/8EPv27cP777+Pjz76CCeccAJOOukkVFZWcnLeHTt2eP33XXfdhUcffRRffvklpyIGzN77xx13HL766iucfvrpuOCCCyCRSDAzM4OHHnoIF154YUzHF4yIRZpOtNvtqKurg9PpxJYtWyKaucPn0EqCx+NBd3c3UlNTUVVVxfmctGhG19hsNtTW1kIqlaKqqipgGTCfkRhhYmICDQ0NUKlUyMjIgE6n46wfSqhpqPlEKpVCq9VCq9WiuLgYLpcLvb29GBoaQktLC1wuF1JTU5kFNTEx8XCZPYcNz0/VP4Wff/BzUDSFk4tPxsrMlXN+h9x7QhUxPovClEoljj/+eObPzp07odFo8P777+Pzzz/Ha6+9xvk5PR4PXnjhBVgsFlRVVXF+fPJ9tNvtqKqqwnPPPYfu7m488sgjOP/882M+vmD2xCKJxIg9U0ZGRlRefnxHYgMDAzCZTNDpdFi3bh0vi2qkQmM0GqHX65GTk4OKioqgCwRfEROJ8Pr7+9HS0oIlS5YAmBW0rq4uKBQKxqop2rJkNkKL+oQkrgqFAikpKZiYmMCGDRu85qaRz+IohQJJADyTkzGfj6Zp3PPlPfj9F78HAFy88mI8ePyDkEvnfsYLIRKLl+1UaWkpzj77bPz0pz/l/PgNDQ2oqqqC3W5HUlISXnnllajL3EPxxhtv4K677sL4+Djuv/9+DA4O4qabbkJdXR3uvffemLJUgorEQgmL7/5SXl5eVAsDH5OXgdmbu7m5GWNjY8jIyJizkc4l4e6J0TTNGG6Ga2nFVzqRpmm4XC60tbVh3bp1SE5OBkVRjGHu1NQUjEYjU5TAdq1Y6G7eQhNUwNsA2HeEyeTkJFPM0XrgAKZzcpiqx5SUlIgExkN5cP2H1+Px+scBADdsvgG3Vt8a8PMk955QP+94VE6SqQh8TnUuLy+HXq/H5OQkXnrpJVxyySX45JNPOBUyco/dcsst2LBhA/785z8zRXdHHXUUzj//fGzbtg2ffPJJ1K9VcCIWaLPZ4XCgvr4eNpst4vRhoHNxCfFmlEgkqK6uRnd3N2+2UEB4IubxeHDo0CGYTCZmInS4x+Z60XU6nWhpaQFFUdi6dSvUajVcLhfz/9l7A2VlZV6uFb29vYyhLvmdYOlZoS5+QiPQd41Mr074ZgTQsrw8DOXnw2QyoaGhgSkRJ59HsFYRm8uGH7z9A7zZ+SYkkOD+796PK9ZcEfS6SGWiUD/HeM4T47NPTKlUMoUdGzZswP79+/HQQw/hb3/7G2fnIJ/hAw88gOOOOw7A4Qe6NWvW4NChQzj77LNjej8FlU4EZm8Q3zSSyWRCXV0d0tLSsHbt2pjTTFyLGKmOzM3NxbJlyyCVSnnfdwslYlarFbW1tZDL5aiurobqm/lQ4cB1JDYzM4ODBw9Co9EwIxlCiaRvUQKJ0np7e9HU1BS0dJwgtOhHaItyyN61bxZQud2O3Nxc5ObmMiXiRqMRo6OjaGtrQ0JCgpd5Mfl+TtgnsOu1Xfhi8AuoZCo8fsrjOL3s9LCuS6ipRJqm4+piH0/HDpqmeTE5AMAIGOD9PVAoFHjjjTdiOragIjFgttudfAlomkZXVxe6urr8loLHci4uBIamaXR0dKCnp2dOc7VUKoXT6Yz5HIEIJjSkJ40tqpEemysBGBkZQUNDA5YsWYLMzEwcOHAg4mOwXSuA2ajcaDTCaDQypeNkEdXpdIIbPgkIT1CB0CLGFHZMTzM/Y5eIFxUVwe12Y3JyEkajEe3t7bDb7dBqtXAmOHHNV9egdaIVWpUWz53xHLbmbw3ruoTe6AyA93uMDO7lKxK7+eabsX37dhQUFGBmZgbPPfccPv74Y+zZs4eX8wFAbW0tenp64Ha7QVEUkpOTkZSUxPi0RotgRMw3enE6naivr4fFYonI3T0cuBAxp9OJurq6gOlNvvbdCP4iMfZIl0haDvwdO9ZFly3wq1evRnZ2NmZmZjhZzFUqFdM8SnzriKARb0EAmJ6eRnp6umAiIKFcByFUxEN/s9kezD9RLpfP6U37uvtr/OCjH2DMMYZ0RToeXP8gimXFcDgcYWUEhNzoHK+iE6vVypgM88Ho6CguuugiDA8PQ6vVYvXq1dizZw9OOOEEXs73xhtv4LrrrmMeeJKSkmCxWEDTNFJSUjAZQ/GQYNKJwGFxmZiYgF6vR2pqKqqrq3kpT49FYKamplBbWwutVovq6mq/6U0+R6WQ47MFwe1249ChQ5icnIxZ9GONxNxuNxoaGjA9PT1H4LmOSNi+dcRb0GAwoKWlBU1NTQDgFaVFklY90gk3nRjJOJY6Ux0u+ugiTDomsTRtKZ464SlonBoMDg6iubmZmclFzIv9RTRCjsTIwy/f12exWACAt0js8ccf5+W4/jCZTLjppptw/vnn44orrsD69etRU1ODL7/8Evfffz8z4TlaBBOJAbM3Rn9/P4aGhrB06VIUFhby8kQWrcUV2+k91KTjeO6JWSwW1NbWQqlUorq6OuYu/1j2xKxWK2pqaqBUKueMl2GLI18b98RbsKWlBZs2bWJSj4ODg2hpaUFiYiJTxh9plV0s+BPvDz+UQa+XYetWN9aupcCTOUPQawornRimiL3V8RYue+sy2D12bMrdhOfPfB7p6lnvP9KbRsr4m5ub4XK55pgXk3tPyCJGJrfzidVqhVwuPyIeuoxGI8bHx3HnnXeis7OTefA877zzYLFYcMcdd+Ddd9+N+viCETGn0wm3243R0VHO04e+EDeNSCBOFwaDgZlOGuoc8Ugnjo+Po66uDnl5eSgvL+fkyx9tOpH0oi1atMjvtcR7FItEIvFq8GU7wDc0NICmaa8qOy5mQIVzXYQXXlDg2WcVAFRISKCxcaMHVVUebN3qwaZNHt49d0OKWATTnZ+sfxL/98H/MU3MT536FDQKb0d9hUKB7OxsZGdnMyXkpAK1s7OT8RBUKBSCTifGY8/VbDbzNgw1XpD7y2KxQKVSgaIo2O12pKamYmBgAMuXL4dGo0FbW1tM5xFEOnFycpIpT1+6dCmvAgbMCozdbg/79y0WC/R6PWQyGaqrq8Na7PiOxCQSCWN6HMixP5ZjRyLAbC/GYP6Q8/2F9HWAn5mZ8ZoBpdFoGEFLTU3lPRo49lg3zGbg889lMBik2LtXjr17Z7+ScjmNykoKW7d6sHWrG5s3exDiuSliwk4nBhExmqZx95d34w9f/AFA8CZmNhKJhBmMSPoEJycnYTKZMDo6CofDgQMHDnhVoAohOotXeT0RsYUMubc0Gg02bNiA2tpaLF68GMXFxbj//vtx6aWX4rHHHsOaNWtiOs+8ihhN0+jp6UF7eztKS0sxPj4el/NGUtgxOjqKhoaGiCMdPiMxt9uNzs5OxrGfS09GILJIjD2LLJAXI4FEYnwPJCUEew0SiQQpKSlISUlBUVERXC4XJiYmYDQa0dTUxMzpCqcXKtpr2bXLjV273KBpoL1din37ZPj8cxn27ZOhv1+KgwdlOHhQhj//eTbPuGLF4UitutqD3NzYotpQBRSh0okeyoNffPgLPFH/BIDQTczBIH2ARLD6+/uxaNEimEwmDA4OMlEzST3G+nlESzzL6xfqLDFfCgoKcO2114KmaWRkZOCaa67B5ZdfjqeffhqbNm1a2HtiHo/Ha/GbmJjg3dMQCC9KYs8mW7lypdf4gHDPwYeImc1m1NbWQiaTMdZBXBNu2s/hcKC2tpaZnh1qYYlXJBbNeRQKBbKyspCVlcWkuti9UGq12itKi2YhC3RdEgmwdCmFpUspXHrpbJq7r0+CfftkzJ+2NhkaG2f//OMfs39vyZLDkVpVlQdLltARza4MO53oR8SiaWIOF4qiIJfLvcaXEPPikZER5vNgjy+JV1tFPBudA80ZW2io1Woce+yxzH+feuqpGBkZwcjISMTrqj/mVcTkcjk2b97MLJixzhSL5LzBRMzhcKCurg4OhyPobLJg8JFOHBsbQ319PQoKCpCdnY2amhpOj08IJ51IKjTT0tKwcuXKsBYR8oVk2x3xSbT7b+xU1+LFi5lJvUajkTHLZY80CWeacqTXUlhIo7BwNlIDgPFxb1FraJCiu3v2z7/+NVu9m5tLMZHa1q0eLFtGIdh6GzIiJlWlPulEk82EXa/twpdDX0bUxBwuvoUd7Kh5yZIlXpOT29ra4HA4vMyL+bQoi+ee2JESiQHe9xr5dy4EDBBQYQfAzXTncAgWJZHy/rS0NKxbty5qdxAu04nsnisSFc7MzPCWrgwViQ0PD+PQoUMoKSnBkiVLwl4w2Dcx+59Chz2pl6ZpWK1WGI1GGAwGdHR0ICEhwWuaMh+LXGYmjTPOcOOMM2ZFbWoK+Oqrw6J28KAMw8NSvPyyFC+/PCtqaWk0qqrc2Lp1Ng1ZWUmB3a0STXXiwMwAzn75bLQYWyJuYg6XUP1rvpOT2ebFPT09XhZmOp2O05lc8YzEFvqeGBv2fcb1A8a8ixh7wYzHnK9A52EXJ5SVlWHx4sUxvdlcpRNdLhfT9M3uueKzDy3QnhhN00yKtbKyElnfeOuFi6+I8QlfT+Jss1xSkECitLa2NjidTmi1WmZ/h50S4vKatFrgxBM9OPHE2fvYZgMOHDgsal99JcPEhARvv63A22/PKldiIo1Nm2b307Zu9UCrBVSqMETM4QCcTjRPd+Lsl8/GoHkQuYm5ePnsl7Eik9vZU0DkfWIajQYajQb5+fmMRRkZBNrU1ISkpCTmIUOr1cYkQvGKxI4UEYtHu8S8ixibSKsGYzkPW8TcbjcaGxthMplCFieECxciYzabUVNTA41GM2cmGTk+H0US/gow3G436urqGDGNJtURTxGL13l8pymTKI2MNFEqldDpdHA4HLy2XKjVwNFHe3D00bP3tcsF1NWRYhE5vvhChslJCT76SI6PPpr92isUlVi+3Irjj1egqsqDzZtnhY2B1aT+ZedHOO/DK2abmHVL8fLZL6MwpZCX1xKLYwfboow0v5MorbGxER6Ph0k9pqenQ61WR3SuI8X8l2/I2hGP90pQIhavdCJbxMxmM/R6PdMozFVzIUknRisyxHNw8eLFKCsrm3MM3/wyl5Abj93nUVNTA7VazcmAz4WSRowGEhUUFBR4lY1brVa0t7djfHycidL47ANSKIANGyhs2EDhZz9zgaKA5mYpU/24b58MIyNS1NUloa5u9u9IpTRWrz68r1ZVJUFSQgJeK7LjgncvhJ1yzGli5gMun9592ypIb5rBYEBnZyfzkEH+hNo+ECOx0JB148svv4RSqcS6deu8fs7G4XBAqVTG9D2YdxFjX3y8CjuIiI2MjODQoUMoKChAWVkZp08N5FiR3vQkZdfb28t4DoY6PtdPO2yBJGbCeXl5WLp0aUznOhLSiZHALhs3m83MoE/23g17CCjX9mpspFJgxQoKK1ZQuPJKF2ga+N//utHQkIr29lx8/rkMPT1S6PWzLiKPPjr793I0eoxoPgaG/4KTqwr8NjFzDV8u9sF607q6utDY2MhMRyBz03zvo3gNxLRYLMye30KDrHmPPfYY9uzZg0cffRRnnHEGk+Fh//Oss87Czp078YMf/CDq8827iLGJVyQmkUiYWVurVq0KKBSxEI2IuVwu1NXVwWq1hqyKZEdLXEO+uD09Pejq6orJTNjfcY+kdGIkKJVKLFq0iNm7IQtod3c3s4ASUeN7CKhEAuTn21FSMoHi4tn0+dCQdwVkU5MMI9ZyoKYc27/TimdP/3XIJmYuiJftFPshoqysDHa7nUk99vf3A4BXlJaQkBC3a7NYLMzk84XK4OAgrFYrbr75ZhgMBvzgBz+YUzQ2MDAQ82gtQYlYPAo7yHBNANi8eTNvLtHkRvd4PGE9Yc/MzKCmpgZJSUlhpezYIsk15Cbr7e2NaJhmKAKlRL+NkMGTOp0OpaWlzALKHgJK/n96ejovUZpvemfRIhrnnOPGOefMZkOMRuDJ3/0b77al4ZmLr4RcKofEYIBkeBjUypWIqCktAuKVsvMlISHBqzdtenp6jqMLACQmJvIekcV7lhgfDA4O4p577oHBYMCNN96I6elp/PznP2daeMhWRazry7yLWDzTiezhmlNTU7x65RGT0HBEhpSsL1myBCUlJWEt7uR3uBYxu92O2tpaAMD69es5tQD7tqUT2YTau2QvoOwKu76+PmYIKKmw85fm4uOa0tOB6x/8Hq5n/Uz+3HNIuPlmuM47D3bSdc0xQjAAZvtuLlmyhHF06ejogMlkwt69e5kqVJ1Ox/n+5kLuEyPvg9lsRkZGBq666iqo1WrccccdmJ6exg033MA8ENhstpgNG+ZdxNjwlU6kaRq9vb1ob29HeXk58vPzMTIyEnaUFC2hKhSJK0h/f3/EJeuk8odLEZuamkJNTQ3S09MxNTXFaX8NIV4mwAsZ3wo7h8PBRGnsIaAkSov2c4pm70nxyisAAM/GjVGdMxyEIGK+EEeXkZERptGdpB67urogl8s57U07Emyn7HY7E01ed911yMrKws9+9jMYjUbcdtttyMjIgN1uP7JEjI9IjD1ni50a49ugN9Q5yFBNu90eU8k6VyI2NDSExsZGZsTM8PAwb/tt7HEsfCI0sYz29apUKuTm5iI3NxcURTEWTAMDA2hubkZycjIjaJGMl4m0slUyMADZ/v2gJRK4zzgjqtcS7nUJTcQIJI1IegULCgq89jdJ5Ew+k2h70xZydSLB5XIxmRy3240LL7wQubm5uPDCCzE8PIyHH34Ybrc75i2deRcx9peIRGJclY0Tn8GEhIQ5c7bisf8WyLVjenoatbW1SElJQVVVVdQbm1xEYjRNo7W1FQMDA1izZg1TEcWlQLKJ9zgWocDVa5ZKpX7HyxiNxjnjZdLT04O2jET6PZO/+ioAwFNVBTonJ9aXEhAhRmIEf9fG3t8EZvfdSQP8oUOHQFHUnLlpwSCtAHzt1/MNuacmJiYYP1Wyth933HH4+OOPsWPHDpx99tlM1W4szLuIsSEbpVxs7JJ9pkB9VvEQMX8iQyKe4uJiFBcXx+wKEsviSKohbTYbqqqqvJ78+BIbMZ3ILf7GyxiNRqYYITExkYnSfCOCSEWMpBLdZ53F+etgE0uzM9+EU9ChUqm8PhOz2QyTyYSxsTG0t7cjISHBy7zY30PsQk4nks/uvPPOQw7rYUcmk4GmaZSXl+Orr77Cjh07ACDmiFNQIkY+TLfbHbWIURSF1tZWDA4OBt1ninc6kX1d7Ign1uNHGy0RN5DExERs2bJlzt5grAIZiHiKmNDEku+F2dcol0xSNhqNjFsFO0qLRCzilUoEhB+JRbI2SSQSJCcnIzk5mTGTJqnHzs5O2Gw2aLVaRtRI9MVnOvEPf/gDXn75ZbS0tECtVqO6uhr33HMPysvLOT3P008/PScTQO63tLQ0fPzxx3jnnXfCMtAOhqBEjDibRysuZEikx+NBdXV10DeH78nL7HMQV3yn0xnyuiIhWhEj06ALCgqwdOlSvwsZX+lEYFZcnE4n9Ho9HA4Hs6jG6msnZOZDUH0nKZvNZq/xMuT7lpCQEHIIaCypRA/lwZB5CAUpBWH9vpBFLFbbKblc7mVTZrPZmAKRvr4+7N69Gy6XCzqdDhaLhavL9uKTTz7B1VdfjY0bN8LtduOWW27BiSeeiKamJk6FM5T7kVwuZ6KxWJh3EfN1N442zWc0GlFXV4fMzEwsX7485NNSvNKJZrMZLS0tSE1NjckVP9DxI53A3N3djc7OzpDToPmMxKxWK/R6PZKTk5GTk+Pla6fT6ZCRkQGdTheTBZgQ01HzeU3siKCoqAhutxsHDhwARVFobm5mxsuQknHfB61oU4kOtwNX7rkSXwx+gfd2vYcibVHIvyPkwg6ue9jUajXy8vKQl5cHiqJAURRee+01fPHFF6iursaKFStw4okn4pJLLsGqVas4OeeePXu8/vvJJ59EVlYWDh48iGOOOYaTc8STeRcxX+RyeUQViuyFedmyZcjPzw9rsYiHiDmdTnR2dmLp0qUoKirixag3XBEjDiUTExPYtGlTyP4vvtJ+FEXh0KFDKCoqQlFREZxO55z9nMHBQabqjkRp0fRGCS2dKCTkcjnkcjny8vKQlZXFeAqOj48z+zZE0HQWS1SpRLPTjO+//n181PcRFFIFmg3NYYnYkRyJBUMqlWLbtm0oLy/H448/jqGhIXz22Wd499130dPTw5mI+TI1NQUATGHKQkNwIhaJuLhcLjQ0NGB6ejqshTna80QKRVFoaWmB1WpFfn4+b/Yx4UZidrsdNTU1kEqlqKqqCivC4TqdSHr13G43SktLUVJSApfLxaQWyYiTpKQkLFmyxKvqrq6ujumNIqLGZ38fHxBB/e1vlSgpoXDOOW5w5DUd0zWRlCLbU5Ds25DxMvkvvAAtAOv69TAnJ0MTRkGI0WbEOa+cg4MjB5GoSMSzpz+L7y7+bljXJdTCDpqm4+ImYrFYoFQqkZ2djXPPPRfnnnsub+eiaRrXXXcdjjrqKKxcuZK38/DJvIuYv6rBcCKxmZkZ1NbWQqPRzCmfDwe+RMzhcECv18Ptdoc99Tdawkn5TUxMoLa2FllZWVi+fHnYT5FcphMpikJTUxPGx8ehUqmQmprKtFIoFAomjcL+PKRSKbKyspCTkwOKojA9PQ2j0Yi+vj40NzeH9BkU4iI4OKjAn/6kBEVJ8JvfzBrxXn65E+n8GcIHJVB1InvfhqZpJNxyCwBgeOtWNO/fD6VS6TUE1DdFPjAzgDNfOhNtpjboEnR48awXsSF3Q9jXJdRIjDzU8X1tZrOZ1wkHbH7605+ivr4en332Ge/n4ot5FzFfwnHtIGXqkdg0+cKHiE1OTqK2thY6nQ4rV65k9nn4IlQkRpphly5disLCwsgaWzlKJ5ICDpfLhc2bN2P//v2wWCxISUmBTCbzss8iT7rk38nDDNnPSUlJYRwsjEajl89geno6MjIyvBZVoaUTtVoPbr/dib/+VYHhYSl+9zsV7r9fie99z4Wf/MSJsrL4Xm84JfbSgQEoDx4ELZEg5+qrkZGZyURpHR0dsNvtXvO5Bh2DOOvlszAwM4C8pDy8uvNVlKdHVvUmdBHjOxIjIsY311xzDV5//XV8+umnyM/P5/18fCEIEQt3ujPZhB4ZGYm5TJ3rEvv+/n60tLR4TYXmuwIykIiRcv6hoSGsW7cO6VE86nORTiRl/MnJyVi7di0kEgmysrLQ3t6Ozs5OpKenIzMzEzqdjkkPsnsFSbTmG6XJZDLk5OQwPoNkUSUlyyTSI5Y2QojKaJpGcjKN//s/J66+2olXXpHjkUeU0OtlePxxJR5/XImTT3bjpz914uijPXx56865plDvjfy11wAAnupq0Dk5kAFMBAzM9jORtO97h97DHZ13YNo9jeKUYry681UUpRVFdV1CFDFyD/J9bcT8l6/7lqZpXHPNNXjllVfw8ccfL3i3fEGIGJtA6USbzQa9Xg+aplFVVRVzmo6rSIykykZHR+cIBtfehr74Oz6xs3I4HDG9T7FGYgaDAXq9HoWFhSgtLQVN0/B4PCgrK0NpaamX8DQ0NCA1NZVJYSUmJkIqlXo59ZO/z47WyHVqtVqkpqairKwMNpsNRqMRU1NTaGpqYsSSzOyaD3d0AlmUlErg/PPdOO88Nz7/XIaHH1bgnXfk2LNn9s+qVR5cfbUT55zjBg/2lQzhiBhTlXjmmX7/PxkC2u5px+2f3g6z24zl2uW4o/wOdOu7YUwxRjxeRsiRGDH25hO+Laeuvvpq/Pvf/8Zrr72G5ORkjIyMAAC0Wi3jsLGQEJyI+UsnksGM2dnZqKio4GQhkslkcDqdMR2DOL7TNI3q6uo5N4BUKuXVld83WiLjXJKTk7Fly5aYyvlj2RPr7e1FW1sbli9fjkWLFsHj8XgtAGwDWyI8BoPBa9ouETSdTsd83qGiNIlEAoVCgUWLFqG3txcVFRWgKIopUHA6nV7NvvP9hZVIgKOO8uCoozxob5fg0UeVePZZBRoaZPjRj9TMvtlll/GzbxZKxCT9/WFVJb7a9ip++M4P4fQ4sa1wG549/VkkK5OZHijf8TJkP81fcQ5N04KOxOJxXXw72D/6zcTTbdu2ef38ySefxKWXXsrbeflCECLmm04kCz9N0+jq6kJXVxcqKio4zdvGGolNTExAr9cjPT0dK1as8CusXAhlMNiR2OjoKOrr61FUVITS0tKYnxajSSeSqsyRkRFs2LABWq12joD5Q61Wo6CgAAUFBfB4PJiYmIDBYEBLSwucTifTO5aRkQG1Wh0wSmMXiJCfE8+6pUuXwmq1wmg0MmXkarWaEbRQzb6xEuqBoKyMxp/+5MCttzrw5JNKPPbY7L7Zb3+rwn338bNvFkosfFOJ/niy/kn8/H8/Bw0aZ5adib9v/ztU8tmyS98eqKmpKUbQiEkuef+Tk5O97jkhpIB9idecM75niQltrzhWBCFibORyORwOB1wuF+rr62E2m7F58+aY7fp9iVbEaJpGf38/WltbQxZMxCud2NnZia6uLqxatcrLqyzWY0dys7tcLsaBY8uWLUhISPDaQwh3UZLJZIxglZeXw2KxwGAwYHR0lBlMSP4/W3jI4kLaLuRyOdRqtVeUlpCQgPz8fKaMnEQJTU1NTKN1OMa50RLOe6DTAb/4hRPXXOPEyy/L8fDDStTXH9432759dt/sqKNi3zcLFYkFSyXSNI0/fv1H/Pbz3wIALl99Of743T9CJvW/yLPHywBginPIFGUSnbOnTAgNvgdhEo4EB/t4IjgRk8lksNls2LdvH5KSklBdXc1LT1A0IubxeJhS8fXr14dsDuRbxIDZCIyiKM6FPpJIzGKxoKamBhqNBps3b/ayDotlMWL3LxUVFTFegAaDAQ0NDaAoCjqdDpmZmUhPTwdN09Dr9ZDL5di8eTNjOEoiM98S/oyMDGRlZXlZMrGNc0nFIxfFIZE+/SqVwK5dbpx/vhuffSbDww8r8c47cuZPZaUHP/mJEzt3Rr9vFkzEgqUSKZrCLZ/cgkdqHgEAXL/5evy6+tcRvUcqlcprCCiZojw4OAgAqK2t9YrShCBq8dqrW8gDMecDQYgY++Ynrg1lZWUxu7wHI1IRI4UlAFBdXR3WVGg+G6ptNhtGR0chkUiwdetWzgdYhlvYYTQaodfrkZ+fj7KyMkY0+NgA9/UCnJ6ehsFgQH9/PxobG5mG6dLSUsjlcub87CpRdqqRXcKv0WiQmJjIiCUp4a+vrwdN08yCysXAw0iQSICjj/bg6KNtaG+XYPduJf79bwXq6mS46io1br+dwlVXze6bRWq4EEzEAqUSXR4Xrn7vajzX/BwA4A/b/oCr110d3Yv7BqlUitTUVKSmpmLRokXYt28f8vPzYTKZmPefvZfGR5QcDmIkJkwEIWLA7A3S3NyMsbExJCcno6SkhNfzRRIlmUwm6PX6qBqG+YjEyPUkJCQgOTmZl0U1nHRiX18fWltbUVFRgby8vLD2v7iCPT5eq9Wirq4O6enpkEqlqK2thVQq9VvCz95LY/8J1GhNxJJMVfa1wyJ7OeFecyyUldF44AHvfbORESnuuMN736y0NLyoL5iI+UslWl1WXPrWpdjTtQcyiQy7T9qNC5ZfENNr8oW4dZAhoL52ZC0tLUyUHO3AyWiJV2GH1WpFdnY27+c5UhCEiFmtVhw8eBASiQTl5eUYGBjg/ZzhREnEKqm9vR3Lli1DQUF4LtwEPsa9EOEoLy+H0+mE1Wrl9PiEYOlEdh/a+vXrmb6seAkYm/7+fqYSMjc3l7m+aEr4AzVaJyUlITk5GcXFxV6N1n19fUyjNVlUA1WEcrmZnp4OXH+9Ez/7mRMvvTS7b9bQIMM//qHE448rvtk3c2Hr1uD7ZoFEzF8qcdI+ifNfOx9fDH6BBFkCnt7xNLYXb+fsNRF8U3bBxsuQgZOk0jU9PT2sDEks1xavSExMJ4aPIERscHAQqampWLZsGUwmE+/GvEBoEfN4PGhsbITRaMSGDRuimj7KZbMzu9Gb7Md1dXXxtucWKJ1IBmna7XZs2bLFq3gingJG0zTa2towPDyMdevWeX0+7Em7XJbwk99hN1qTirvu7m40NjZCq9Uyosa3dZBSCVxwgRu7drmxd+/svtmePXK8/bYCb7+twJo1s/1mZ53lf98soO2UTypx1DKKs146C4cMh6BVafHfM/+LqrwqXl5TqIrJQONlRkZG0NbWBo1Gwwga1xWn8YrERBGLDEGIWFlZGbNQxMNdPtR5yKgQYpgb7dMdV+lEp9OJ2tpauN1ur340vsalBDo2iZj9FXCQ/q944PF40NDQAIvFgk2bNoVs6A5Vwp+WlobMzMywS/jZHnopKSnQarUoLS316ovq7u6GQqFARkZGxAMoI0UiAY45xoNjjrGhrU2KRx9V4N//VkCvl+GKK2b3za680oVLLz28b0b6sfxdEzuV2DXZhTNfOhM9Uz3ITszGK2e/gpWZ/BnFRlI84TtexuVyYWJiwqvilLRYcNEXKEZiwkQQIsb+IkU6iiVaAokYKVTIzc3FsmXLYnry4iKdOD09jZqaGqSmpmL9+vVe6So+qx9904kmkwm1tbVYtGgRysvLmcU83ulDMvhUJpNh06ZNEVeuxlrC79uHxu5rUiqVyM3NZfYHSUqzvb0dDocDHR0dyMnJ4dUYeulS6pt9MyeeeEKBxx5TYGhIit/8RoV771Xi+9934cc/dqK42H8/FjuVWHvMUpz53IkYs46hSFuE13a+hiWp/FoUxSL2CoUCWVlZTMWpxWKB0WjE2NgY0xfIjtIiFSQxEhMmghMxIi7hWOLEAok0yEJM0zR6enrQ0dHBWWN1rOnEkZERNDQ0oLi42G+lJt8iRiIxYiRM5rXN1/4XmVyg0+kiKrAJRKQl/CqVKmiU5lscQvqiysrK8PnnnyMlJQUGgwEdHR28N1qnp9O44YbZfbMXX5z1aTx0SIa//12Jf/xjdt9s69Z0bN3q/fmRVOInp6zAmR9ejCnHFFZlrsLLZ7+M7ET+iw24KmNnf7aLFy+G2+1morSWlha4XC6kpqZ6ubeEupfjGYmJ1YnhIwgRY0N6e/gWMRLRkIji0KFDmJycjHguWTCiFRmaptHR0YGenh6sXr06YKUSnyJGjt3S0oLBwUFm32m+BMxgMDCOJEuWLOHl3L77LTMzMxgfH2dK+FNSUpgoLSUlJWSUxi4OkUgkTKUke0Ftbm6G2+32ssPisjhBpQK+/303vvc9Nz75ZHbf7L33yL7ZUfjvf9245hoXzjrLDYViNpX45lLg3I0tsDvcqMqrwvNnPI/UhFTOrikYfPViyeVyZGZmIjMzEzRNM+4t5KFCpVKF9Nj0eDycTmb3B4kgk5OTeT3PkYTgRIzcJG63m9d+HPJFMZvNOHToEBQKRdgDIyM5R6TpRLfbzQz63LJlS9CbmevBlWxomsbY2BhkMhk2b94MjUYzLwUcgP8KRL5hV8WVlJTA6XQyxSF9fX2QSCSMoJEhnb5RGvljMpngcrkgkUjgdDq9yv/JomUwGLyKE8hx2WIZ2+sBtm3zYNs2G1pbpXj4YRn+8x8F6urk+OEP5bjtNgo/2jWODGs7rtkFeODGycUn46lTn4JGwd9MPF/i4ZtI+gkTExNRWFjI7JOaTCa0t7cz42VIxSkp0DlSbKeONAQhYuwFkdzAfBd3kPMcOHAAeXl5KC8v5/zLE2lUabVaUVNTA5VKhaqqqpAizldhh9VqxfDwMCQSCeN8MR8FHMEqEOONUqn0cpiYmpqCwWBAd3c3/v3vYYyNZePKK53IyvIu4R8ZGUFjYyOWLl0KrVbrN0pTq9UoLCz0SmkajUY0NDR4Nfqmp6dz8mBXXk7hj3+047jjvkZr6zb8/e8qDA1JcdufsgFFP/D5nbjgqj48fMLDUMjiO0F7Phzs2fukgPd4ma6uLigUCqSnp8dNXMQ9scgQhIixIXO4+CzuoGka3d3dAMAM1uQD9lN5qCc4UlBCCifC+SLzkU4kk6DVajU0Gg0TTUokkrguLpFWIMYT9n5Xbm4ZLrtMhf5+GfbuncJPfnIAhYVuZGRkgKIojIyMYPXq1V6z70I1WmdmZnqlNNmNvklJSYwdViSN1r7QNI2UFCduuMGJn//cjRdekOOWX/dgwliBbZpsPHrSTZBK4m/1xGcVZ7iQ8TJk/3dychImkwkWiwUzMzMwmUzMg0W442XCxe12w263iyIWAYITMSC86c7RQtJ1U1NTzBMWX7CjykAiRtM0+vr60NbWFnFBCdciNjg4iKamJpSXl8PlcmFycpK59oVUgRhPEhOB66/34OabpdDrtfj5z7+LX/3KhC1b9LBazZBIJOjv74fdbg9Zwk+idnaUlpiYiKSkJCxZsgROp5NptNbr9ZBIJF6N1pG8TySCl0gkUKmACy90Y9euHNz3x3/iZxecdVjAnE7wOtTMB6HNEmM3s09PTyMzMxMymWzOVHHSlxjrvWo2mwFA3BOLAEGKGF+9YhaLBbW1tVCpVKiursaXX37Ja9qSvVD5gwzUHBsbi6qhmisRI2m7/v5+rF27Funp6RgZGUFXVxc+//xzpocqPT2d9z0BrisQ+UYiAX74QzeOP96DH/9YiU8/leG229JRWbkWjz9Oo7CQwvj4eMQl/P4araVSKbKzs5Gbm8uY5rJHm6SkpDBRWqhGa5Li9m1vuenGs5j/ln38MRKuuQa2//wH1MqVMDvNSFLyGyEIdZYYMPvZqNVqZGZmzhkv09PTw3wGJEqLJlImDjxiJBY+ghAx3w+aj16xsbEx1NfXIz8/H0uXLmWehvkWsUDFFw6HA7W1taAoClVVVVE1YnIhYm63mxl5s2XLFqaAIyMjA8ceeyympqYwPj6OtrY2OByOOY3BXBKPCkS+KCqi8eqrZtxxxzj+9rclqKvTYds2Gr//vQuXXz5rLOx2u5mKOHYJPxG1UCX8vo3WycnJTOGJ3W5norTe3l7I5XKvajvfqrpQ+7TShgaov/99SGZmoPjrX/HUT7bi1k9vxZvnvImKjAre3kehRWJsfPvE2Gnl0tJS2O12Zi+tr6+PcY6JxDjaYrEgISFhXieQLzQEIWLA3MGYXIkLTdPo7OxEd3c3Vq5c6VXdxqUtVCD8Cc3U1BRqa2uRlpaGlStXRn3Dhus0HwibzYaamhoolUps2bJlTgEH276JNAazo4rExERkZGQgMzMTWq02JtGZjwpELiFTvs89NwGXXpqDq69WY98+GX72MyVefVWG3budKCiQzynhNxgMGBwcZIyFIynhZzdaKxQKptHa1zvSZrN59URpNJrgY1j6+qDeuROSmRm4jz4aD1y8FL/acxUA4F+N/8Jdx97F2/sohD2xQIRysU9ISJgzXoYImm+UFmi8j9ls5t2u7EhDMCLGhqvCDhJlzMzM+C1Xj4fFle85hoeHcejQIZSUlMQcbcQSiU1OTqKmpgbZ2dlYtmyZ16IY6EmYlCWzx5WMj48zI2qIoKWnp4fdTyOkCsRoMZvNqKmpQUZGBuPy8u67DuzeLcfttyvw4YcybNyYgHvuceLii2dNedkl/MXFxVGV8AdrtCajTcrKypieKFJtp1KpkJKSwvx9r4XZZIL67LMhHRmBe3kFfv2L9bj7s1sBAD9e+2P87pjf8fpeCjkSi6TEnv0ZlJSUwOFwMFEaMTj3N16G70bnTz/9FPfddx8OHjyI4eFhvPLKKzjTz9DThYQgRYyLwg6z2cxU2QUqV4+HiBGhYe87VVZWIisri7NjR8rQ0BBT9l1QUMBcXyTl8wqFAjk5OcjJyfEqOSeu8WlpaYyoBaosFHIFYriYTCbU1dWhsLDQy1VFKgV++lM3TjrJgx/9SIkvv5ThJz9R4ZVXPHjkESfy8rwj6GAl/IcOHYJWq2XSuKSEHwiv0VqlUiEvL8/LO3J4eBgURWHv3r2HG601Guh27YKsrQ3uvEW45pb1+Kv+QQDALdW34Jebf8l7hCB0EYv22lQqldd4GRKlkSg8KSkJb775JnQ6HecVj2wsFgsqKytx2WWXYefOnbycI94IRsS4TCeOjo6ioaEBhYWFKCsrC3hDxEvEXC4XampqYLFYsGXLFs42bUmfWLh9aDRNo729HX19fVizZg0yMjK8LL6i/eL4WixZrVYmqiCedWQBJsUMC6kCMRCkB2zZsmXIy8vz+ztlZTTee8+Bhx+W4447FHj//dmo7N57nfj+9/2PSvF9P8Nx4Q/UaO0vSiP7M5OTk1i7du2sv+DwMFJvuQXyL7+ETZuEi26uwEvt/wIA3P/d+3Hlmiu5fwP9INTCDvaw11hhz8IjUfjw8DBaWlqwb98+2Gw2nHPOOdi+fTtOOukkTizwCNu3b8f27dyP0JlPBCNibKJNJ7LtmlatWoUc1kTaQOeJh2N+U1MTkpKSUFVVxeliHUkfGtsJZPPmzUhMTOTNQkqj0aCwsBCFhYV+ixm0Wi2mp6eh0+mwcuVKQS5aoejt7UVnZ+ecHjB/yGTAtde6cfLJHlx5pRIHDshw1VWzUdnDDzuRmxt8XzOQC39raytTbEOi3kAl/L6z0hwOB9NoXZCfj6UPPQTll1/Cqpbj1F8twcfjH0AGGW5deStOyzoNDocjLhOVKYri3dopGkjGg4+CC6VSicWLF+O///0vnnzySfzzn//EmjVr8OSTT+Kqq65Cb29vwIckEYGKmFwuh8vliujvuFwu1NfXM9FOOH0WfIuYwWCA1WpFZmYm1q1bx3mKIFwRs9vtqKmpgUwmw5YtW7zStXxbSMnl3sUMfX19aG9vh0KhwNjYGA4ePOiVJhP6hjZ7D2/9+vUR+WyWl9P44AMHHnpIjjvvVGDPHhk2bEjA/fc7sWtX8AGWBLa7BPEAHB8fZypIwynht9ls6O7uhlarhcfjQcJDD0H5j39gMgE47bYSfO5oQIIsAX87/m9Yo1mDoaEhtLa2Mo3WwQoTYkWohR18ihgbu92O7Oxs3Hrrrbj11lsxNTXFmZfrkYogRUwmk8Fut4f9+6S3KDExMaJoh68Se7YjvkajQW5uLi9fTHLMYBWKU1NTqKmpQWZmJioqKsIq4OCLgYEBdHR0YMWKFcjNzYXdbofBYMD4+Dg6OzuhUqmYiCItLU1wERoxip6eno56D08uB37xCze2b5+NymprZfjhD1V49VU3HnrIiRDJAy/YHoD+Svg9Hg/TN0ZK+Elrh1arRUVFBRT/+Q/Uv/sdxhKBE3+5CHWOViQrk/Gf0/6DrQVbIZFImEZrUphQV1cHiUTiZYfFVYZBqHti8frO+FpOiQIWGsGIWLQzxci4kqKiIpSWlkYkFnzYW7EnQm/cuBHt7e28RXuhmqlJJWRZWRkKCwvnbQZYoArEhIQE5OfnM/Y+JpMJ4+PjaGxshNvtZkxyMzIyeDWDDgcy0drj8WDTpk0xX8/y5TQ+/tiBBx6Q4667FHjzTTn27ZPhT39y4pxzwovKfPGNetkl/E1NTUhMTITdbkd6ejpWrFgB+YcfQvXTn6JXC5zwszS000PI1GTixTNexMqMlXP20rKyspCTk8MUJpBKyubmZqbROlYrJqGKGNkP4/t7I/omRo5gRIxNOGk+drVfsHElsZ4nEkivEABmIjSfvWikGMP3+GRvsLe3F5WVlcjMzJy3ESrhViDKZDKvURlkAe7v72d6bEiUxmX11mefSeFwAMcdF/gzIulYtVqNtWvXcpZSksuBG26Y3Su76ioV6uqkuPTS2ajsgQeciKWA1beEf3x8HPX19VCr1TCZTNA/8QSqbroJLalunHiVGgOyCRSmFOKN895AaVppyEbrpKQkJCcnM+Xj7EZrtlWTTqeLaI9LqIUd4kBM4SJIEQsViTmdTtTV1cFut8dU7celiE1OTqK2tpZ5yiULHZ8zv/wdn4jG1NQU7wUcoSCpK5lMho0bN4YdvfguwA6Hg6nO6+npgUKhmFOdFw2NjRLs3KmC0wn8+98ObN8+93MiqWp2DxjXrFpF45NP7Lj/fjnuvluBV1+V47PPZHjgASfOPjv2+5OkF0lETnd1IeGyy1CbasVJl8hgUtiwWLMYjx/1OHIUOV5CEk6jtVwuR05Ojld7AOlJa2xshFarZfrdNBpN0HtQqJFYPAdipqam8nZ8s9mMjo4O5r+7u7uh1+uh0+lQWFjI23n5RDAi5m+6sz+mp6dRW1uL5ORkVFVVxVTJxJWIkXRNWVkZFi9ePGe0DN/WVmRBIZGgVCqNewGHLzMzM9Dr9UhLS4vZA5H0ORE3iomJCYyPj6OlpQVOp5OZvpyRkRHRQMmyMhrHHefBa6/JsWuXCk8/7cSZZx7+rAL1gPGBQgHcdJMbp5ziwZVXqnDokBQXXaTCa6+58cc/OvHNlJCIIW0AjBOKwQD12Wfjk8RxnP59KcwKD9ZmrcXurbvhmfHgq6++CquEP1ijtVarZayYbDabV6O1UqkMOnxSqIUdR0okduDAAXznO99h/vu6664DAFxyySV46qmneDsvnwhGxNgE2qsiTbrFxcWcLCqxihhN02htbcXAwADWrl3LzCPyPQdfkdiHH0rw2mtLIJFIsX79FOMasXz58nkt4CBP/osXL+bcA5EMlExPT/eywiJ9NklJSYyghaqgUyqBf/7TiR/+EHjhBTkuvliJxx6brRQMpweMDyoraezda8c99yhw331yvPiiHJ98IsOf/+zE6adHdq/29/ejvb0dlZWVs/em1QrVeefhDWk7zr8QcMgpHFt4LJ4/63kkK2erecMt4QfCa7QmTdzssSZGoxFtbW1wOp1eE63VarUYifE81Xnbtm28zCCcTwQpYr6OHRRFoa2tDQMDA1izZk3IvpxwiUVgyEa/zWZDVVVVQKsYPiOxF1+U4YknKqBSjcHj+RqlpaVYvHgxI2Dxjr6A+HogSiQSJCUlzRlVMj4+zhiwksU30N6MXA48/rgTCQnAM8/I8cMfKjE0NIw1a5oOL/5xRqkEfv1rF047zY0rrlChuVmKCy5Q4bzz3Lj/fidCTQ8i8/J6e3uxbt262fSU2w3lpZfiX/av8IPzAY8U2FG6A0/teAoJ8sPRq78SflJBGqiEP5JGayJabDus8fFxpine5XLBYrEITsziGYktROea+UQwIhYoneh0OqHX6+F0OoOKRTRIpdKoqhOJVx4p6Q+W0oz2HOGgVM4+UY2OTs57AYcQPBCVSiVj7UNMcMkCabfbAzrwy2TA7t1OqFQ0/vEPBX7960W4665EfPe78+sisnYtjc8/t+P3v1fgT3+S47//nY3K/vIXJ0491f+DEfkcRkZGsGHDhtmnepqG4rrr8LDxLfzfN5NWLlx5IR456RHIpYHvXXYJ/+LFi8Mq4Q+n0Zocm1Snkqb4iYkJNDU1YWBgAH19fYzgkWPPJ/GKxKxWqzhLLEIEI2JsSIQ0OTkJvV4PrVaLdevWcd7JH00kNj4+zuyTBLO0Yp/D6XTGcpl+8Xg8mJ42AchFSkoGMjLUgq9AjCf+HPgNBkNAB36AxuWXH4TJlIOXXy7CLbdoIZE4ce21/E0YDweVCrjjDhd27JjtK2ttleK881T43vfcuPdeJ9jPCmQ+3eTkJDZu3Mh8DrL77sVvux7HnSfP/t41G67B77f9PuLJzcFK+In/H3lPg7nwB5qVRvrNKioqIJfLYTQaMTIygra2NiQmJno1Wsc7SotXJEZc7EXCR5AiRsRq//79nLi9ByKSPTGSouns7Jwz0iUYfFQnOhwO1NTUQCIpBgB4PPNXwBFtBWKkmM0ATQPRPqSyIwriwG8wGKDX65lqPJlMhr/+tRSlpS7ce68CN9+shN0O3Hjj/AoZAGzYQGHfPjvuvFOBhx6S49//luOjj6R4+GEnTj55VhDq6+tht9uxceNGJnKR/PNpXLf/t9h97OxxfnP0b3D95utjvkeCufDP3puRufATUXO73aBpGhqNxmtiAmm0bmhoAE3TXo3W8eghjEckRtK3YiQWGYIRMfKlIvtfALBixQosWrSIt3OGK2IejweHDh3CxMQENm3aFFEXPdd7YtPT06ipqUFaWhoKC2ftHWZm7LDbpVCr1Qu2AjEYBw5IcdllSmzZQuHvf489qmU78NtsNhw4cICZn/bFF/tw2mlpsNvL8ec/Z+G3v1XCZpPg9ttdUTUgc0lCAnDnnbNR2VVXKdHeLsXOnQn4/veduOCCGhywfoBFhYtQpaoCAHj2vIWrPrga/9kESGjggRMexBVrr+Dl2sJx4SeilpSUNCdKI96eSqWSGcwKHJ5rl5mZyUSAxAF+YGAALS0tSE5OZgQtmmnK4SDuiQkXwYgYMPtUr9fr4Xa7IZVKkZKSwuv5SDoxmAs8aXSVSqWoqqqKODfPZXXi6Ogo6uvrUVxcjCVLluCdd2Z/brV6sG/fl0xlXmZmJm9fZoLBYEBdXQNGRirwwgsFeOCByLwuI8HlAnp6JOjqkuOUUzw46yxuHgr89YDZbDaMj4/j3HPbMDMzjiefXIH77lNgYsKOP/1JApls/osNNm+m8MUXdtxxhwIPPyzHs88q8fo7KzFzyp2Qlf0PyzKWYcWAAxe9egHeXklDTkvwjx1P4NyK8+Jyfb4u/Gx7MVJmzy7hB4CGhgY4HA5s3LgRcrk8aAl/cnIykpOTmQiQlPD39/czaeSMjAykpaVxaod1JFQnHokIRsSsViu++OILpKWlYcWKFfj000/jMrASmH3K8rffNjExgdraWmRlZUUdaXCRTmSnMletWoXs7Gx4PB4oFLMilZKSiWOOOYap9CLj6Ymg6XQ6Tp8iSQVicfEKXHzxEkxOSnDaaR6ceCI/rQRVVRR+8Qs37rtPgZ/9TIktW2yItfCR9ID5tgGo1WrGgX/1ajeKioy4/fZ0/OMfyRga6sUttwwjOzszbmmsQKjVwO23T2Px4g789r5iTI9mAv/ag9ITPkHGMTKc8Z/t+KzEA7VHimfPfg4nLT113q7V117Mt4RfJpNBJpOhsrLS6z0Np4RfJpMhOzubKeYhjdbd3d1MozWJ0mIxmI5HJOZyueB0OkURixDBiJharcayZcuQnZ3N3Jx8VfUR2F8SXwYGBtDc3IylS5eisLAw6ps/1l40YjprMpmwefNmJCUlMV9olWr243M651bmkYbg5uZmuFyuqHwIZ5wzsLvtyNTMtjT4q0C84AI3Hn1UgaeekuPEE7kvYCHcfLML770nQ12dFD/+sQqvvOKIOr0Xbg+YXC7H9dfLkZnpwNVXK/H224uhUCTj6qv1zAI5Xw78MzMzOHjwIOqyn8P0D14CPvg98NW1aH3/WFR+2QvXWcdA6/wYL537KqqWHhe36woFu4Tf6XSipqYGbrcbCQkJ+Prrr6HRaJh7NZoSfnajtd1uZ6K07u5uKBQKr0brSArFPB4P7xWSZrMZAETbqQgRjIhJpVKv+V9cTHcOBVl0fHvSWltbMTQ0hHXr1iE9VFNOCGKJxEjRBE3T2LJlC5RKpVcBB/lO+RY/+jYEm81mjI+PMz6EZPHNzMwMWAk1ODOInS/tRII8Ae+c/w6UUqXfCsRLL50VsbfekmFkBBG5sEeCUgk8/rgDRx2VgPffl+Ef/5Djiisif8ghc8Ai6QG75BIPVConrrxSidde00GpPAZ/+cs0pqbmx4F/cnISB2oO4EnDk3i191VACfz+HivW2e047/xxTE8thnzPH/HW64NYu/Q7IY83H7hcLtTW1kKlUmHjxo3MQ2ukJfwAAkZpCoUCubm5yMvL82q07ujoYFouyPck1D5UPPrWLBYLAIjViREiGBEDuJ3uHO75gvWkcbHBGq2IkSft1NRUrFy5EgC8Nrtn3RBmf9fhCBwBSCQSrz0E3/EnCQkJjKBptVpIpVLUjtbi3JfPxbB5GJmaTLQb2mHttfqtQFy5ksamTR58/bUM//qXHNdfz1/0XFFB43e/c+GGG5S46SYFtm3zoKwsPPcBdv9UpHPAAGDXLg8SEpy45BIlXnhBDocjBU8/neDlwG8wGHh34DcYDDigP4BHxx7F/wb/B6lEikdOegQXr7oYAIX6GilO3PlP3HlNIdauEqaAkepajUaDVatWMeIQyoU/OTmZETRyrwIIO0ojogWAabQmopaQkMAIWmpq6pz9L4/Hw/ueGCnqEFKT90JAUCLGJh7pRHIej8eDmZkZ1NTUICUlhdOetGiqE8fGxlBXV4clS5aguLiY+VKSSi0CWRsjaUNj70+43W5m/EldXR0AoMnThNvqb4PNbUNFegWePulpmDpMQSsQL7vMja+/luHpp+W47jo3+PwO/uhHbrz9tgwffSTDD36gxAcfOBBq756Mx5menvbqn4qUM8/04D//ceL731fi9dfluOACCZ591oGEhMMO/MuWLWMi34GBAU4d+EdGRrC/fj8eHHkQX45+CZVMhad3PI0dZTuY38nMSUPt5+cc/ks0DYyNAVFMeeADu92OgwcPIiUlBStWrAi4YPsr4Sd7vrW1tRGV8NM0PSdKI36cBQUFTKO10WhES0sLXC6XVwl/QkJCXAo7SI+YEL0jhYxgRSwe6URgVsQMBgO6urqimkkWzvHDjcTYwzTZBRyBGpgDpRPDRS6XIysrC1lZWaAoCn/8/I+4Y/8doEFjTfIa3LH0Doy0jCA/Px9Lly4N+L7s3OnBjTfS6OqS4pNPpPjOd/h07Qf+9jcnNm1KwMGDMtx7rwK33BK4MtLlckGv14OiKE7mgJ1yigcvvODArl0q7NkjwznnqPD88w6QDJBv5OvrwE8KbiJ14O/v78eBpgO4Z/AeNBgbkKxMxvNnPY9jC4/1/xdoGtJ334Xizjshcbth37cPvD5dhIHNZsPBgweh0+lQUVER0ffMd8830hL+cBqtyRggi8UCo9GI0dFRxmrL5XIhOTmZ17Si1WoVy+ujQFAiFu90InlS6+zsxOrVq7325Lgi3HQiRVFobGyEwWDApk2bkJKSEtKBg9hOxWoI4qbcuOGDG/CY/jEAwA8qf4CfLPkJBnoHkKBOQH9/PyYmJgKW7ycmAuef78bf/z5b4PGd7/BX4AEAeXk0HnzQiUsvVeGee+Q48UQPNm6c+x6z54CtXr2asyfp44+n8MorDuzcqcJHH8lw1lkqvPSSw28jdqwO/DRNo6urCwfaD+CuvrvQOdWJDHUGXj3nVazNWTv3hDQN6YcfQnHnnZB9/fXsjxITIWlpAb18OSevPxosFgsOHjyIrKwslJeXx/SgGKiEnzyMEqd8UpkbyoXfd1aaWq1GQUEB0xg/MTGBlpYWDA4OYnBwkInSdDodp8UeYiQWHYISMTaRTHeOBrfbjUOHDsHlcmHp0qW8CBgQXjrR6XSitrYWHo8HW7ZsgUqlYl57MAeOcPbEQjHtmMYlb1yC97rfgwQS3HXsXTgp5SSMDM7uHaWlpXmlcgKV71966ayIvfaaDOPjAEcezQE591wP3nrLjRdemDXt3bfPDvZ+ON9zwI4+msIbbzhw5pkqfP65DDt2qPDqqw6kps4KDw16jq2TPwd+g8EQ0IEfAFpbW3Gg5wDu6LoDQ5Yh5Cfn443z3sBS3dI51yTduxeK3/0Oss8/BwDQajXcV14J1//9H/8fSBDMZjMOHjyIRYsWcZ7pAEKX8BMX/oyMDK89p3BmpRETaaVSidLSUiiVShiNRsZqi91oHWpqQijEgZjRIVgR4zMSs9lsqKmpgVwuR0pKCq+lszKZjMnJ+7vByRc8JSUFq1atgkQimVPAEQgiYq4o+4z7p/ux86WdaDQ0Qi1X4x/b/4EiexGMRqNXBWI45ftZWZlYu7YItbWzlkjx8Bx84AEn9u2ToqNDiltuUeDBB2ffiEA9YFyzeTOFt9+24/TTE7B/vwynnJKAV1+z4t76X8LqsuKRkx4JeG62A39RUZFfB36ZTIbmqWb8rut3MNqNKNeV4/VzX0d+Sr7XsaRffQXFb38L2ccfAwBolQruyy+H6/rr+SsXDRPiMFNYWMjrZ0EI5cKvVquZ/clAJfy+jdZkPy3Q1ATiHiKRSJhGa51OF3GjtcViESsTo0BQIsa+weVyORwOB+fnMJlMqK2tRW5uLpYtW4aamhreh1YC/jv+SUHF4sWLUVJSwnxpfAs4AkG0N5q3qWakBue8fA5GLaPITszGv0/7NzAEuGSuoB6Iwcr3q6osqK2txN//TuMHP7AgKYnfL2Ra2uz+2GmnJeDvf1dgy3cMOHaTFU1NTXGbA7Z2LY2337bjtNMSUFcnReXRJkyf/xKQNIaLV12MLXlbwjoO+0HB5XKhpqYGX45+id+1/w42yoZlycvwxDFPIF1xuOVDWlMzG3m99x4AgFYo4L70UrhvuAF0HGegBYJMOy8uLsbixYvjfn5/LvykkClYCT/gbYTQ0tLCpBmDNVoTO6ze3l6moId8V8Ip6BFFLDoEJWJs+IjE+vr60NraimXLlqGgoIC387AhXwp2iS5N0+jt7UV7eztWrFiB3NzcqBzoFYro9sTeaH8Dl715GWxuG1ZkrMDTJz0NQ6chYg9E3yKG4mI7nnqKQne3Co8/3oYNG6xzyve55jvfoXD11S488ogCP7xSggt+92f84exfxHUO2KpVNF5504TjTqYx3V8IPPUp/vTPurAFjA0pRPnc+Dl+1/E7OCknjs47Gn/a9CfYpmz4/PPPkT0ygmX/+Q+0JPKSyeC58EK4brwR9DyIhT9MJhP0ej3KysqY79p8wy5kCubCT0r4AaCtrQ0TExPYsGEDVCpV0BL+lJQUpKamoqSkBA6Hg4nSent7IZPJGMEM1GgtphOjQ9AixtWeGEVRaG5uxujoKNavX8/4tZHz8Clivq4gZFzG2NgYNm7cCK1WG/UIlcPVieH9HZqm8dD+h3DrJ7eCBo0TlpyAB6ofQE9rDwoLC2Oelp2ZmYDzzqPw1FNS1NVtwgUXDHqV75M0Tnp6OqdjdSoveB747zrQ48ux559nYvcPUjk7djiMW8fxc/1OOC+ahOSZD0EbyvGXq8uw/W0HCgvDn6JL+qc+MH2Ae5rvAUVTOL3sdDx52pNIkCdA0twM+R//CMVrrwEAaKkUg9u2wfiTnyB57drZUnO+XmQEGAwG1NfXY9myZbwaeMdCoBJ+g8HAlPArFAq4XC6sX7+emT8X7qw0uVyOnJwcxhCZNFp3dnbCZrMhNTXVq9FaIpHERcR2796N++67D8PDw1ixYgUefPBBHH300byek28ktIBmVbvdbkZQhoeH0dvbiy1bIn+aZUOKJtxuN9atW+c1DBEAGhsboVAosHTp3I1yrnj33Xdx9NFHQy6XQ6/Xw+VyYe3atUz/SbQzwDo6JKisVCMlhcbwsC3o77o8Lvzig1/g8brHAQBXrrkSP1v6M3R1dHE6hfngQSmOOSYBKhWNjg4bdLpZ8ZyamsL4+DjGx8dhtVqZqrzMzEy/VXnh8tDXD+HmT24GhtdA8o+vQXsU+NvfHLjwQv7bMwCgd6oXZ7xwBton2pGuTsfuLW/hpss2oqtLivx8Cm+/7UBJSeivmNVqRU1NDd6aeAt/bvkzAODiVRfjLyf+BYrObih+/3vIXngBkm++ru5zzoHzppswmZOD8fFxGAwGWCwWpKamellhxZuxsTE0NDRgxYoVvBVL8Q2ZjzcxMQGlUgmbzTanhN/3u+rbaM1eVn333Ww2GxOlkSKUd955hykiefjhh3l5Xc8//zwuuugi7N69G1u3bsXf/vY3/OMf/0BTUxMKCwt5OWc8EJSIkXlCwOyXob29HVu3bo36eGRTmbhe+Hv6b25uBgBUVFREfZ5QvP/++1i9ejVaW1uRlJTEuBSwq6CiiYD6+yVYtkwNlYqGyRRYxKYcU7jo9YvwQc8HkECCu79zN45POh4jIyOorKzkdAozTQNVVQloaJDivvuc+MlP5kbTVquVEbTJycmo3PcpmsKNH9yI3bW7AQBXrbkKObUP4Y47VEhOpvHVV3YsXszvrX1o/BDOfPFMDJuHUZBSgNfPfR1LdUsxNCTBaaep0NoqRU4OhTffdKCiIvC1EHeWFydexJMdTwIArt14LX6/+Aoo774bsv/8B5JvHu7cp58O1y23gP7GxYWNzWZjihhMJpPfIgY+GR4eRnNzM1atWoXMeayGjAWaptHS0gKj0chEYOwSfqPRyJTwk0Zr371uf43WZJkl+93knx6PB3V1dXj44Yfx0UcfYWpqCieeeCJOOeUUnHLKKViyZAlnr23z5s1Yt24dHn30UeZnFRUVOPPMM/GHP/yBs/PEG8GKmNFoxKFDh3DssQGaOUMwMjKChoaGb/ZqAqfJ2tra4HQ6GWsnPvjf//4HmqaxePFilJaWej2pxbKwjIwAJSUaSCQ0ZmZsfg1xe6d6sfPlnWg2NEOj0OCJU55AgbUAFosFa9eu5aW58rHH5Pi//1Ni+XIKX39tD2rU63K5mIXXYDCE5b7vcDtw+RuX49WOVwEAdx5zJ36+6eegKAlOOkmFL76QYetWD955xwG+TBb2DezDuS+fi0nHJJZnLMdr57yGRcmHU2ejo8COHQlobJQiI4PGG2/YsXr13K8a8UF8xvQMXux+EQDwu8pf4MY3JiD/5z8h+eb74Nm+Hc5bbgG91k9vmB/YRQwGgwEURTG9U3w48A8ODqK1tRWVlZUx+43OFzRNo7m5GSaTySuFyIZ4MJL31V8Jvy++jdaBorSLLroIS5YsQX5+Pt566y24XC7s3buXk9fmdDqh0Wjwwgsv4KyzzmJ+fu2110Kv1+OTTz7h5DzzgWD3xKJ17KBpGh0dHejp6UFlZSWysrKC/j6X87780dfXB7fbjeLiYpSVlUW9/+UPsidG0xK43Zhjv3Rg+ADOefkcjFvHkZuUi3+f9m9Qg1TICsRYOf98N26+WYGmJim+/lqKzZsDv7/EpDVc9/1J+yTOefEcfDH8BeQSOf62/W/YtWIXAEAmA/7+dye2bEnA55/L8Oc/y/F//8d9qf/bHW/jojcugt1tR1VeFV44+wWkJXhHs9nZwDvvzJbf6/VSnHJKAl5/3YF16w6/F+Pj46ipq8Fjhsewp3/PrA/i5FG46vy/QPJNtY7nuOPg+vWvQW3cGNE1+hYxTE9Pw2AwoLe3l3HgJ1FarA22fX196OzsxNq1azmN6uMJW8A2bNgQMMVNCjSIUJNeP98SflLAEaqEn91oPT4+jmOOOQbXXXcdfvGLX3C6V28wGODxeJDtYz+WnZ2NkZERzs4zHwhKxNhfpGgKLtxuN+rr62E2m7Fly5aw5vLwVdhBURRaWlowMjIClUoFnU4XVgNzJLA1yOn0FrFXWl/BD9/+IexuO1ZlrsLTJz2NsY4x3qcwA4BWC5x9tgfPPivHk0/KsXlzeOWTodz3HSoHfnXoV+gydyFJnoTnzn4O31nsbXK7ZAmNe+914ic/UeGOOxQ47jiP3wgoWv7Z8E/89N2fwkN7sL1kO/6545/QKPxHs+npwFtv2XHWWSp8/bUMp546O0JmyxYKw8PDONhwEH8Z+ws+G/4MSlqKf70iw7l1nwIAPMccA9ett4KKIZ1OkEgk0Gq10Gq1KCkpmTOkMhYH/u7ubvT09GDdunURmyoLhXAFzB/+Svh9jaCDlfCTKK2urg779+/Hhg0bmGPz4dXou+4EGwi8UBCUiLEhkVi4bzLZGFepVMzYknDgQ8RImbTD4cCWLVtQU1MDh8MBmqY5EzDgcCQGzPaKJSbO3pR/+vpPuO3T2wAAJxefjD9W/xHdLd2cVCCGy+WXu/Hss3K8+KIM99wzK2yR4Fu+XzNQg52v7MSYfQw6hQ6/X/575DnyMDExMad8/+KLPXj7bTfefFOOH/xAhb177YihdgTA3Pf1wpUX4pGTHoFcGvwrlJoKvP66A+eco8Jnn8lw+ukq7N7dD3nKV7i/727UTjQg0Qm8+hyF47soeLZsmY28tm2L7YKD4OtwEWrh9Qexw+rv78eGDRsW7CDHWATMl0hL+Mm+WFNTE84++2z8+te/xq233srhqztMRkYGZDLZnKhrbGxsTnS20BDUnhhFUXB9Yz/hcrnwwQcf4Pjjjw9Zjm00GqHX67Fo0SKUl5dH9CQ5NDSE/v5+bN68OaZrJ1gsFmbMxOrVq5mbdGhoCMnJycjKyuIkhQPMFlEkJ6tB0xJ0dlqRnunCte9fi6cbngYA/Hjdj/HT0p+is6OT0wrEcK9t48YENDdL8eCDzqhmfxE+6fsE5718HswuM0q1pXjt3NegcWmY4hBgbvn+2BiwaZMa4+MS/OxnLvzhD1HammC2iOSmj27Cwwdnq8au23QdfnvMbyP6/KxW4PzzVfjwQxkUCjeyz92FgaUvQWcF3nkWWJ+9bla8TjgBUU/7jBGappno12AwYHp6GsnJycweJanKo2ka7e3tGB4exvr16xdsbxOXAhYKdgm/wWBAe3s73n33XWzatAm7d+/Gj370I/z2t5HdU5GyefNmrF+/Hrt372Z+tnz5cpxxxhliYQdX0DQN5zd7ARRF4b333sO2bdsC3lw0TaOvrw9tbW2oqKhAfn6+398LxujoKDo7O1FdXR3TtQOHxTQ/Px9lZWVem7hut5tZdI1GIzPHKysri3kqi4b0dDXsdgm+0o/ixppd+LjvY0glUtzznXtwfPLxGBoa4rwCMVweeUSOX/5SidWrKezbF7zAIxD/bf4vrnz7SrgoFzbnbMZL577ktf8UrHz/4MFF+P73kyGR0HjrLQeOPTbyvU+nx4kfvfMjPN/8PADg7u/cjWs2XBPxcWiaRn19G376kzTU6IsAmQPpZ+zAxyM9KPv5XfCccsq8iVcgSMMuuWflcjkyMjJgt9thNpuxYcOGBeu6Hk8B84WiKDQ0NODhhx/GBx98gLGxMWzduhWnnHIKduzYwVuRGSmx/+tf/4qqqio89thj+Pvf/47GxsZ5cVThCsGKGAC899572Lp1q99+F3bTcCwbysRRPNaGv/7+frS0tKCiooJxLPd4PH7Thx6PB0ajEWNjYzAYDJBIJF4VeZHkwnNz1ZielqD41yejS/YukhRJeOLUJ5BnyeO1AjEcTCagtFQNh0OCvXvtXkUN4cD0gAHYUbIDT53+FBLkwRcb3/L9xx5bj7feysOiRR7s329Hamr4QmFxWvC9176H//X8D3KpHH89+a+4YMUFEb0G4PCEgqmpKSxfsQabj/0a48YKfHD7R1j/w13zPiIlHCiKgslkQltbG6xWK+MTGMyBX6jMp4ARurq6sH37duzcuRPXXXcd9uzZg7fffhsKhQIvvPACb+fdvXs37r33XgwPD2PlypV44IEHcMwxx/B2vnggaBH74IMPsHHjRsbRm+BwOFBbWwuKorBu3bqYbkKTyYT6+npsi3IPgqIotLa2YmhoiBHTSCoQSTc/WXidTidTkZeZmRnSRDQv24NJczLw45VYVDKB53Y8B9eACzKZDJWVlbxVIIbL5Zcr8fzzclx+uQt/+Ut4KT2KpvDLD36JR2tn+1l+tOZHuO/4++a4wofC5XKht9eIU05ZhMHBBBx33CDuuWeEKWAI9rBgtBmx86Wd2D+8HxqFBv86/V84qfikiM4PgNm0dzqdWLt2LVQqFSxDQ+ixKbGiJH7WWLFCURQOHToEs9mMtWvXwuPxMMUhU1NTzH5PZmZmzG7ufCIEAevt7cXJJ5+MU089FQ8//LA4yTlGBCViALxMfz/++GOsXr3ayyZqamoKtbW1SEtLw8qVK2Ou4JmamsLBgwfx3e9+N+K/63K5UFdXB7vdzohpLA4c7D2JsbExmM1mpKWlMYLm27cyPD2EsqUq0DMFKD97Hd6ZkMKxdAPsJ5+M0uOPF8SXY+9eKU4+OQFJSTQ6O20ItX1id9tx+RuX47WOWWulu465C9duujamRfGrr6Q4/ngVKEqCu+/uwpo17X7L9wn90/0444Uz0GpqhS5Bh5d2voRNizZFfF6Xy8VYGK1Zs2buA8nMDKTNzaCVStBr1kT9+viGoijU19cz97nvgxG7189oNDLOEySzwKXFWCwIQcAGBwdx0kkn4bjjjsPf/vY3QXxHFzqCFrHPPvsM5eXlTPf/0NAQGhsbUVJSwtlYB7PZjC+++AInnHBCRH/ParXi4MGDzMBFdl9btA4cvthsNiZCm5iYYJwtsrKyZjfZ7XZ898TP0D/lwtt9v0Klq5P5u9TKlfDs2AH3qafOLpDzViwArF2bgPZ2KR55xIFLLw1cCTphn8A5L56DL4e/hFwix2PbH8P5K87n5Dp++1sF7rlHAZ2Oxldf2ZCcPMO8tzMzM9BqtcjMzIRRZsSuN3dhyDw7u+u1c1/DsvRlEZ+P+CAmJCRgdXk55J2dkDY2QtrUBGlTEySNjZD29QH4xkLq6ac5eZ1cQyJJt9uNtWvXhswMkMwCKWAge5RE1Pw1EMcDIQjYyMgITj75ZFRVVeGJJ57gpYT+24jgRMzpdDLFEF988QWWLFmC7OxstLW1ob+/H5WVlZxa2thsNnzyySc46aSTwhYeMs6FVENy5cARDPK0OzY2BqPRCIVCgaysLGRkZGDGPIPhzz5HZW8v0j79FNLPP4eE1cBN5efDc9pp8Jx2GqijjprbFc0zDz4oxy23KLF+vQeffup/bszA9AB2PL8DbZNtAXvAYsHlAr7zHRVqa2U4/ngPXn3Vweg66Zv6sP1D3FB7A8weM5YkLcHzpz2PiryK8D9TjweSnh64amow9uGHSBsaQlp/PyQdHYxtlC9Ubi48p5wC15//zNEr5Q63243a2loAwNq1a6OKqMgepcFgwMTEBDQajddkg3ikHYUgYGNjYzjllFNQWVmJZ555RjDR6ZGAoEVs//79yMrKYsxN161bx3k5r9PpxIcffogTTjghrCejgYEBNDc3Y9myZcjPzw9awMEXpLdnfHwcw8PDjKVQXl7ebD/IxARke/ZA9tZbkP3vf5BYrczfpVNT4TnxRHh27IDnhBOAOPT3jI8DZWVquFwS7NtnQ2Wl9y3XMNaAM/57BkZto8hSZ+H1817HqqxVnF9HS4sEW7cmwG6X4E9/cuKqqw6X/e/p3IMLX78QNrcNazPX4p7Ke+Camt3Dm+O+T9OQjIxAcuiQd2TV0gKJzb+HJZ2aCmr5clArVoBevnz235cvB1ipciFBUqFyuRyVlZWcRA0ul8vLCgsA0zeVnp4e8RDJcBCCgBmNRpx66qkoKyvDc889x8vr/DYjaBH7+uuvYTabkZKSgsrKSl4+fI/Hg/fffx/f/e53gxZB0DSN1tZWDA4OYs2aNdDpdJxaSEUKcdo2m80oLS3FzMwMxsbGYLfbodPpmH40pccD2UcfQfbmm5C99RYk3yweAEArlaC2bYPntNPgPuUUgMc+sosvVuKll+S48koXHnjgcIHHx70f4/yXz4fZbUZZahneOP8NFKTwN3/q0UfluP56JdRqGp9/bkd5OY1nDz2LH+/5MTy0BycVn4RndjyDRGXirF1Tby8sX30Fl16PhPZ26IaGkNTTA9nUlN/je5RKOEtLIV+zBvSKFaCWLwe9ciXo3FzBldAHwul0Hk6FftPryDXs1gi+HPiFIGATExPYsWMH8vPz8eKLL857odWRiOBEzOVygaIoGAwGHDx4EKmpqdi0aRNvIkHTNN59992g/Whutxt1dXWwWq3MOJdYCjhixeFwQK/XQyqVzqlAtFgsGBsbw/j4OKanp5m9nqysLGhUKki//npW0N58E9KODq/jejZuhOfUU+HZsQN0eTmni+5HH0lx2mkJ0GpnR7RoNLM9YFe8dQXctNtvDxgfUBRw5pkqfPCBDOvWeXDGPffi9s9ny/i/n3sy/kqfDlVzK6SNjZA0NUE6NOT3OLRUCmt+Ptzl5ZCuWQNbcTEaJRLkHXMM8hdwz43D4cDBgweRlJSElStXxq3wgO3APzExgYSEhJgc+GmaRlNTEzPQcj4EbGpqCqeffjoyMjLw6quvBnQ/EYkNwYmY0+lEV1cXOjo6oNVqkZqayuusL2C2H626utpvqpJtZ1VZWclLAUckzMzMQK/Xh+WB6HA4mEpHk8kEjUbDRGgpycmQtrVB9sYbs2nH/fu9/i5VWnp4H23TJsRqB09RwOrVCejuluJvf3NgvOx+3PrprMXO6SWn48nTnwzZA8YVQ0MSbNqgxMSUDDj2DuA7v8Ev9gH3vg9I/XwbqIKC2RTgN5EVtXw5nMXFMHxTSTo+Pg6KoqDT6bB48eKQ5ftCxWaz4eDBg8y9NV9l8mwPQvLeRuLALwQBm5mZwVlnnQWNRoM33nhj3gpavg0ITsRqamowNjaGdevWYXh4GDRN8zrrCwjcjzYxMYHa2lrk5OSgvLx8zlygeH/JDQYDGhoaovJAdLvdXg3WMpnMu8F6dBSyt9+G7I03IP3kE8ZFHQDozEx4TjkF7h07Zj39ovxC3n+/HLffrkT2sk6M7ioFAPx4zY9x7/H3RtwDFiu3/n4vHrjrJEDixg1rluPe2nbQ6emHhYrsXVVUBDV+7OvrQ3t7O4qLi5mHhmDl+0KFVNtmZGRg2bJlgunzYjvwj4+Pw2w2B3XgF4KAWSwW7Ny5ExKJBG+//fa8DCf9NiE4ERsaGkJiYiISEhLQ0dEBm82GVau43+Rn8/HHH8+xZhocHERTUxPKy8tRUFAwr/tfwKwjSFtbGyceiOyRJ2NjY/B4PEhPT0dWVtbsBrvNBtn//jcbpb37LiSs/R9ao4Hn+ONno7STT561ag+T3kEHVixLAk0pgJ+swP19o/j5yBIgNRV0WhrotLTZf//mv5GWNuffodFwkuakPR5sPfpNbCofw58vKpktssjODvvYbAPctWvXMg7u7F4/3/L9zMxMZhS9kDCbzTh48CByc3NRVlYmuOtj429AJXlYSEtLQ0tLCyYnJ7F+/fp5ETCbzYZzzz0XTqcT77zzzoI1Rl5ICE7E3G43k67r7u7G1NQU1vDcCLp3715UVFQgIyMDNE0z5fxr1qxBenr6vAoYMVvlywORuG2TfTSLxcLYCWVmZiJBJoN0797ZlOObb0I6MHD478pkoLZund1HO+000EVFQc9lmjRhXfVBjPefgitzr8LfBp+M/HqVytCil5oKWqeb8//A0Z4EKfIhGYNgFbMkOiPTlolnJikxn+9mVzJVuqCgIG4TDrjC4/EwD2PE7UYqlaK4uBi5ublx34Oy2+244IILMDk5iffee2/BjqZZaAhaxPr6+jA2NuY1Y4cP9u3bh5KSEqSnpzPzyNatWweNRjOvBRwej4ex+lmzZk1c0hKkr2dsbAxTU1PezvsaDaR1dZC/9dZs2vHQIa+/S61cOVvpeNppcxqsTSYT6urq4LBmgtaM4JTsRZBMTEAyOQlMTBz+98lJSEymw//+zf/D5CQz5ThaaLV6Npr7RgCZf09NhWfbNlDbt4c8BvFBnJ6eZop8woV4ZgZz348nU1NTqKmpQVFREZYsWRLXc3MJTdNobGzExMQEsrOzMTk56eXAn5GRgeTkZF6/v06nExdeeCGGh4fx/vvve7kMifCL4ETM4/EwwyO5HpMSiK+++grZ2dkYHByEQqFgyvnns4AjWAVivHA6nV4N1iqVipmXpNVqIe3tPRyhff65V0MvlZ8P169+Bc9ll2FkZASNjY2oqKjAokWLorsYmgYslllB+0bwGHFjCx35OfkdIpIhbnPXtdfC9fvfB/0dtg+iP/ulyF5OYPf9jIwM3gsBJiYmoNfrUVJSgsLCQl7PxSdkD8w3hUjuXV8H/mhMtkPhcrlw6aWXoqurCx988AEyMhaOJ+aRgKBFjMsxKcHYt28frFYrcnJymCISMjJc6BWI8cJfFEFK93U6HWSTk4cbrN9/HxKrFY5HH0XnMcegq6sLq1evnr8vN0UB09NBoz7PUUeBOvnkgIcgzb9SqRRr1qzhPGrydd8nFmOZmZmcRxFGoxF1dXUoLy9HXl4eZ8eNN4EEzBeyB0xEzeFwcObA73a78cMf/hCNjY346KOPkJWVFfWxRKJD0CJmNBrR2NjI66iAoaEhNDQ0ICsrC5WV/9/eeYdFca5t/F6KSmhSFuyigKCCSNFYohElKIIsaIqJGvUznmgSY9qnSTwnJucymgTPiZ5ESdFEUlAjTVCxUywx0QWkiQGlKGWXIrggy5aZ7w+/mbPgooC7O7Pw/q4rf2QBeXZc5573fe/nfnzY7UNm6qqheRwHoqGgKApNTU3sOVpbWxv7lOvo6Ih+ajVM0tLwl4MDKuVy+Pr6PuD8NCbkcjmys7NhYWEBb29vvdvnOwbqmpqa6mwVIZVKkZ+fj7Fjxxp0SKqu6aqAafs5zQeGx0ngV6vVWLt2La5cuYK0tDSjvp7GDK9FrLGxEdnZ2QgM1F2GHgNN0ygpKUF5eTmsra3h6OiIkSNHcupAvH37Nq5fv27wKcyPg7bkfVtbW6jVaiiVSvj7+xvt4ETgv32CdnZ2GDu2GzmKOkLTSfq49n1mW9fb29uoVww9FTBtMA8MzH9MAj8ThdXZipuiKKxbtw7nzp1DWloahg/XX8oM4eHwTsQoioJSeT+WSCaT4dKlS91OmH8UKpUKeXl5uHv3Lvz9/VFaWgqBQAB3d3eYmpr2OgeiIZHJZLh69SqUSiXUajUsLS3ZczRmvL2xIJPJkJWVxRvr+ePY96uqqlBUVMTttq4O0KWAdYTZYWCisDpL4KcoCu+++y5OnjyJtLQ0uDzClUvQL7wWsZ4kzD8KuVyOrKwsmJqasjOeqqurce3aNbYB2MnJCXZ2dgZ56ubCgagvWltbkZ2djSeeeALe3t5sfBhzUzA3N2evb0+ihAwJY3xwcXGBi4sL5wKmjY72/f79+7OCpnl9b926heLiYjbz01jRp4Bp4969e+2isPbs2QMHBwcolUqkp6cjIyMDrq6ueq2hM6KjoxEdHY2ysjIAwPjx4/HRRx8hpAsO294Gr0Wsuwnzj4KxFDs6OmLcuHHs7wPu/wNpbGxkz3nUajV7w3VwcNDLOQgfHIi6QiaTITs7m0186ChQTE8Pc30pimJvuFzYyx9GbW0t8vLyMGbMGAwbNozrcrqEpvGmrq6Ovb4A2H62gQMHclvkY6ApYAEBAQbvAVOpVIiNjUVMTAyysrJgYWGBsLAwhIWFISQkxOA9YSkpKTA1NYWb2/3km5iYGERFRSE7Oxvjx483aC1cwzsRo2kaiv+PPKIoCidPnkRgYOBjf2irq6uRn58PNzc39uyrsxlgmvZniUTCGheY+V26SNNvbm5mJ1TzxYHYU5gesJEjR3ZpWKnm9dVM3mdEjcug1KqqKly7dg1eXl5wdnbmrI7Hgbm+xcXFaGxshEAgMKh9X9dwLWBMDVu2bMEPP/yA06dPo7m5GUeOHMGRI0fw2Wef8WIFZG9vj6ioKKxatYrrUgwKr0UMAE6cOIEZM2b02BxA0zRu3LiB0tJSdqCmWq0GTdNdss/TNM0mwzPGBc1RJz35B2UMDsSuUl1djcLCwsfqAesseZ/JxjMUFRUVKCkpgY+PDxy6EafFNxjTUlVVFfz9/WFiYtLOvm9pacle3+648biAaWRuamriVMCioqLw9ddf4+zZs5gwYcIDX+fyGqrVahw6dAjLly9HdnY2u8vUV+C9iJ05cwaTJ0/uUQYZM3OrqamJjQd63Aip1tZWVtCamprajzrpgtAaowNRGzRNo7y8XOc9YJrnPPX19ewkYCcnJ73dcJkHndu3b7fLQTRGNCOx/P39H3gI0Kd9X9fwRcB27tyJ7du349SpU/D39zd4DZ2Rl5eHqVOnQi6Xw8rKCrGxsZg/fz7XZRkc3okYcP9GxqAtnLcrML09AoEAEydORL9+/XQeIdVx1MnDnHi9yYHI3CglEolee8A6Ju+bmJiwK2B7e3udbMHSNI2ioiLU1tbqZXK4IdEcAunv7//IbcOO9n2FQtG+34/DM1q+CNju3buxdetWHD9+XO/JQd1FoVCgoqICjY2NiI+Px549e5CRkUFWYnxAU8Q0w3m7CmPgcHBwYA859Z3AwTzhMjdczYgmKysrFBQU9AoHoqab0tfX12A9YB2T91UqFXtO2dPR9hRFIT8/HzKZrNs5iHxDM9OxJ849bfZ9Gxubdtu6htoyYwSMeS9cCdiePXvw0Ucf4dixY5g+fbrBa+guQUFBcHV1xbfffst1KQaFP5YwDQQCAWu6MDMzY5ufu0JNTQ3y8vLg6uoKFxeXhxo4dIm5uTkGDx6MwYMHs04xqVSK7OxsqNVqmJubw9PT06hvlEqlEjk5OaBpGpMmTTLok7qJiQkcHBzg4OAADw8PNnm/tLQU+fn57ZP3u3AD18xBNPR70TUURSEvLw/37t3r8apFIBDA2toa1tbW7Waj1dbW4ubNm53a93UNXwTsp59+wt///nekpKQYhYAB9+vWXAD0FXi5ElMoFKzw/Pnnnxg6dOgjM96Y+U7MGY2TkxPnM8Cam5uRlZUFS0tLWFhYoK6uziDWfX3QsQeMT3VrS95nrrG2FYS+cxANCSPGSqUSvr6+ehHjzuz7um6P4IuAxcbG4u2338bhw4cxZ84cg9fQFT788EOEhIRg+PDhkMlkOHDgAD777DMcP35c5+EQfIf3IiYWiyEUCh+atM1scd25cwd+fn6wtrbmXMC0ORCZKbVSqVRv1n198KgeMD7RWfI+s4Joa2tDVlYWL8W4u6hUKnZlzDTu6xvN9oi6ujq0tLTAzs6OFbWe7jTwQcAA4NChQ3j99ddx6NAhXtjmO2PVqlU4c+YMqqurYWtriwkTJmDjxo19TsAAnoqYUqlkz7CuXr3KbnFog7kpAYCfnx/Mzc05nQEGdM2BqA/rvj7obg8Yn+iYvE/TNGiaZv/R8/WhoSswq0kmeYYrMe6YatET+z5fBCwpKQmrV6/G/v37ER4ezkkNhO7DexHLz89H//794e7u/sD33b17lw1nHT9+PAQCAacjVB7Hgfi41n19oIseML7AmH0sLCygUqkeTN43ojMxhUKBrKws9O/fHxMmTODNarIn9n2+CNiRI0ewcuVK/PTTT1i0aBEnNRB6Bu9F7Nq1awDAzvlikEgkyM3NxejRozFq1CiDGTg6Q5cZiN2x7usDffWAcQWTgzhq1Ci4uLhoXQUPHDiQfWjgs/mm43YoX7d2tdn3mfR9oVCIfv368UbATpw4gWXLlmHPnj1YvHgxJzUQeg4vRUylUrFTlYuLi9HW1gYvLy8A92+wpaWluHHjBry9veHs7AyKoqBWqznbPtRnBuLDrPu2trY6f7+G6gEzFF3JQZTL5exDA7Mlxmzr6nusfXeQy+UQi8WwtbU1qqgy5qGBEbS7d+/CxsYGFEVBpVIhICBA72G+nZGWloYXXngBu3fvxrJly3jzd03oOrwXsZs3b0Imk7EDK/Pz81FfXw8/Pz/Y2NhwbuAwZAaipnW/traWbf7VVeq+5mrS2PumgJ7lIGpuifEpeb+1tRVisRj29vYYO3asUd9s5XI5cnNz0dzcDJqmDWbf78i5c+fw7LPPYseOHfif//kfo76mfRnei1hFRQVqa2vh7e2NrKws0DTNWom5NnDU19cjNzeXkwzEjs2/j2vd1+wBYxJOjJny8nLcuHHjscaPUBSFhoYGzpP3W1paIBaL4ezsjDFjxhj1zZamabbB3N/fH2ZmZg/Y95lztJ42sXeF33//HZGRkfjss8+wdu1ao76mfR1eipjmdOfKykqUl5dDqVTC1tYWXl5eMDExYUWOCwMHwK8MxMe17vO5B6y76CsHUdNaXltbi9bWVoMk7zODOYcOHQpXV1ejvtl2FLCO14z5HDPXWFf2/Y5cuXIF4eHh+OSTT/Dmm28a9TUlGIGIXb9+HaWlpXBzc8Po0aPbGTiMzYFoCDqz7jOrtI43DqYHTCgUwsPDw2jOWbTBZAfW1dVpDb/VJcwZj1QqZc94mHM0Xf3epqYmZGdnsyt9Y+ZRAqYNXdj3O5KTk4PQ0FB8+OGHeO+994iA9QJ4K2JKpRJlZWUoLi6Gubk5Zs2axbmBwxinMHe07jM3WycnJ7S2tiI3N5fX04u7Cpc5iPpI3m9sbER2djZGjx6NkSNH6qFqw9ETAetIR/u+iYkJK2hdTd/Pz89HSEgI3nnnHXz44YdG/Xkn/BdeiphKpcLVq1dRW1sLNzc33LhxA0899RSn51+9YQpzR+s+TdMQCoVwdXU1iHVfX6hUKuTm5uo1eqk7tXRM3mcEravJ+w0NDcjJyTGqydKdoQsB6whFUWhsbGQ/y9rs+x25du0aQkJCsGbNGnzyySdG+1knPAgvRezWrVsoKSmBn58fFAoFLl26BDc3Nzg7O3Nixe1NU5iZHrAbN25g+PDhaG1tNYh1X18oFArk5OTA1NQUPj4+vMpB7NgrpVQq2zVYazurZFoCPD09jb7BnEnW16WAdaQz+75QKIRKpYKrqytKSkoQEhKC5cuXY+vWrZz9+922bRsSEhJQVFQECwsLTJs2DZ9//jk8PDw4qae3wEsRU6vVaGtrg0AggEqlQmVlJSQSCbsd5uzsbLCmVC4diLqmsx4wfVv39YVcLmcDlr28vHhtSKFpGjKZjF09MKYF5hxtwIABkEgkyM/P71ZLAF8xhIBpg9ltkEgkWLBgAZ544gmYmpriySefRGxsLKer9Hnz5mHx4sWYNGkSVCoVNm3ahLy8PBQWFhrF0QRf4a2IyeXyBwwcmh/QO3fuwNraGk5OTnB2dtZLNBOfHIiPS1d7wHRt3dcXLS0tyMrKgr29PcaNG2d0DxfMWWVtbS0aGxsxYMAAyOVyjBkzBsOHDze696OJpoAFBARwJhyFhYVYs2YN2traIJVKoVQqMX/+fGzevFlrjJ2hqa2thZOTEzIyMjBz5kyuyzFaeCli69evR0VFBcLDwzF//nxYW1s/8D0KhYK90dbX18PS0pJdoT3udF6+OxC7S097wDSt+1KpFHK5nBep+0xm5tChQ+Hm5mbUN3wAKCsrQ0lJCWxsbCCTydjmX6bB2pjeH18ErLKyEnPnzkVQUBC++eYbAPfHOqWkpGDNmjUYPnw4J3VpUlJSAnd3d+Tl5bGJRITuw0sRy83NxYEDB5CQkIDy8nIEBQUhPDwcoaGhWs9rOkYzWVhYsCu07hoWjNGB+DB01QPWXeu+vmBS9ZkcRGOnoqKCbcq2s7N7IHkfQLsGa76shLXBFwGrqanB3LlzMX36dOzdu5eX14ymaYhEIty5cwfnzp3juhyjhpcixsA4m+Li4pCYmIjr169j1qxZiIiIQFhYGOzt7R8QKMYdJpFIUFdXh379+rErtEfZnXuDA1ETffaAPcy6r6/UfalUivz8fHh4eDxySKoxUFpairKyMvj5+WltyqZpup0Lr62tDQ4ODuxKmE+fT74ImFQqRUhICPz8/BATE8Mro48mr7/+Oo4ePYrz588bvQOVa3gtYpowpoT4+HgkJCQgNzcXM2bMQEREBBYsWAAnJ6cHBKqjYcHMzIy90XbcpulNDkTgv4YUQ/SAGSJ1n8lB9Pb2hpOTkw6q5g7NVBF/f3+t2+XafoZZCdfW1kImk/EmeZ/p0WtubuZUwOrq6hAaGgoPDw/s37+ft/Pi1q1bh6SkJGRmZmLUqFFcl2P0GI2IaULTNG7evMkK2pUrVzBt2jSIRCKEh4djyJAhD9w4mRw8iUSC2tpaCAQCdsuRoqgHpjAbM1zOAdNH6j4zFsbHx6fHOYh8gaZp/PXXX5BIJPDz8+vx+S1fkvcZAWtpaYG/vz9nAnbnzh2EhYVhxIgROHToEK9WqQw0TWPdunVITExEeno6L8wlvQGjFDFNaJpGRUUFEhISkJCQgN9//x2TJk2CSCSCSCTCiBEjtAranTt3IJVKUV1dDbVaDTs7O7i4uHS5IZWP8G0OmOZKuK6ujn1w6Kp1n6ZplJSUoLKykp1aYMzQNI2ioiI2FktX266dJe8LhUK9tkjwRcCampoQHh4OoVCIxMRE3kxF78hrr72G2NhYHD58uF1vmK2trdFPjOASoxcxTWiaRlVVFRITExEfH4/z58/Dx8eHFTTNAFW1Wo3i4mJUV1fD1dUVra2tkEgkvLWUPwq+zwHrrnWfyUFkxu4Yu8GGoigUFhaisbER/v7+ertpMTsOzHVmUuGZ66yrMyK+CJhMJkNERASsrKyQnJzMazHobHX8448/YsWKFYYtphfRq0RME5qmIZVKkZSUhPj4eKSnp2Ps2LEQiUQIDg7Gpk2b4Obmhs8++4y9QTKWcolEwsbZODo6wtnZGY6OjrwVNGObA/Yo676pqSny8vLQ0tICPz8/zgYm6grNMyNDNv5qXmddJu/zRcBaWlqwaNEimJiY4OjRo0b/oEPoGb1WxDShaRoNDQ04fPgwYmNjcf78eQwZMgSRkZF44YUXtBo5mIQFZryJ5o1WKBTyxvVk7HPAtFn3TU1NYWZmhokTJ3bJ9MBn1Go1cnNz0dbWBj8/P07/fnSRvM8XAWttbcVzzz0HhUKB1NRUo/+cEHpOnxAxhvz8fISGhmLKlCmYO3cuUlJScOLECQwdOhQRERGIiIiAj4+PVkFraWlhV2gtLS1wcHCAs7MzhEIhZy6o3jQHDLjfwC4Wi0FRFMzNzdvdaPVp3dcXarUaOTk5UKvV8PX15ZVbTjN5v6GhARYWFuwKrTMDDl8ETC6X48UXX0RTUxNOnDihs5lxBOOkT4lYUFAQZsyYgY8++oj9RyqTyXD06FEkJCQgNTUVjo6OCA8PR0REBCZNmqT1UJxZOUgkErbpl7nRGuoftmYPmKenp9E7KjVzEL29vWFiYmIQ676+UKlUyM7OhkAgwMSJE3mzctcG01vJiJq25H2+CJhCocDSpUtRXV2N06dPG32aDuHx6VMi1tbW9tBzgHv37uH48eOIj4/H0aNHYW1tjfDwcIhEIkydOlXrSodp+pVIJLh79y4GDhzI3mj1dZZjyB4wQ8DkIDo4OGDs2LFa348+rPv6QqlUIisrC+bm5vDx8TGqFTIz5oQ5R1MqlXBwcIBcLodarea0D0ypVGL58uUoLS3F2bNn4eDgwEkdBH7Rp0SsO8jlcpw6dQrx8fFITk5G//79ERYWhsjISEyfPl3r1pBcLmfPdhobG/WSuM9lD5g+6EkOYsdopu5a9/UJsyVqYWGBCRMmGG27BvBfY0hBQQFaW1tB0/QDyfuGQqVSYdWqVSgsLERaWprRN7wTdAcRsS6gUCiQlpaGuLg4HD58GDRNIzQ0FJGRkXj66ae1Ppl2lrjv5OTUIxeVZg+Yj49Pr3gK1UUOoubKgevUfWZL1MrKCl5eXkYtYMCDZ2BqtZrd3m1sbISVlRUraPrc3lWr1Vi7di2uXLmC9PR0DBo0SC+/h2CcEBHrJiqVCpmZmYiLi0NSUhJaW1sRGhqKiIgIzJ49W+vTqbbEfc2A4kfB9x6wnqCPHMTOrPuMYUGfxorW1laIxWI2toxP25s94VFnYAqFol2Dtb6S99VqNd58802cP38e6enpvSIzk6BbiIg9Bmq1GhcuXEB8fDwSExPR1NSEefPmISIiAs8884xWN113E/eNrQesK1RWVuL69evw8vLS27aQNus+sxWm69R95kyPCVru7QLWEbVajYaGBvYcDdBN8j5FUXjnnXdw6tQppKenY+TIkT36cwi9GyJiOoKiKPzxxx+soEkkEgQHB0MkEmHevHla+1jUajXq6uraJe4zgmZjYwOVSmXUPWDaKCsrQ2lpqcFzEPWVut/c3AyxWIzBgwfD3d29zwlYR3SVvE9RFN5//30kJycjLS0Nrq6uPXk7hD4AETE9QFEUsrOzERcXh4SEBFRUVCAoKAgikQjz58/X6qbrmLhvamoKtVoNKysr+Pr68tqi3RX4lIOoK+s+Y0oZPnx4rwiO1rWNXlvyvq2tLXuO1tnDA0VR+Oijj3DgwAGkp6djzJgxj1UHoXdDREzPMDPRDh06hMTERPz1118IDAxEREQEQkNDtc5Ea2pqQlZWFvr16welUtkucX/gwIFGZxjgcw6i5vZufX09uxp+lHWf+TvqLcM5mUkO9+7d01sfGJO8zzRYW1pasudoTPI+TdPYsmULfvjhB6SlpWHcuHE6r6M7ZGZmIioqCmKxGNXV1UhMTERERASnNRHaQ0TMgDAp5syQz7y8PMycOZOdiSYUCpGcnIyysjKEh4fDxcUFNE2ziftSqRQ0TUMoFMLZ2dkoEveZm6Mx5CB21brf0NCAnJwcuLu782LM/eNiCAHriFKpbDfhIDo6GlZWVrC2tkZSUhLS0tLg7e2t9zoeRWpqKi5cuAA/Pz8sWrSIiBgPISLGEcxgRGYmmlgshq+vL/Ly8vD3v/8db7/99gOrAOa8gWmu5nvivkqlwtWrV6FSqeDr62tUZ3qdWfcHDBiA8vJyeHp69gqnHBcCpq2G5ORkfP/99/j999/Rr18/NjVn7ty5vMlFFAgERMR4CBExHkBRFD788EPs2LEDfn5++PPPPzF58mR2hMzw4cO1ChqfE/cVCgWys7NhZmYGHx8foz7TY651aWkpu0Jj3HdcZmc+LnwQMOD+9d21axe2bduG1NRUmJiYICkpCYcPH8bLL7+MjRs3clJXR4iI8RMiYjxg48aN+Omnn3Ds2DFMnDgRVVVVSEhIQHx8PC5cuICJEyeygqbNQMC3xP3e1vQLADU1NSgoKICXlxcsLS0NYt3XJ3wSsO+//x6bN29Gamoqpk2b1u7rarWa8wcyBiJi/ISIGA84ffo03NzcHjAI0DQNiUTCzkTLyMjAuHHjIBKJEBERgTFjxmgVNC4T97uSg2hsVFVVoaioSOu0bH1Z9/UJnwQsJiYG77//PlJSUvD0009zUkdXISLGT4iIGQnMTLSkpCQkJCSwwicSiRAZGYmxY8c+NHFfKpVCJpPpNXG/qakJ2dnZGDZsWLsp2sbMrVu3UFxcjIkTJz6yr80YUvcZAWttbeV0vhlN04iNjcU777yDw4cPY/bs2ZzU0R2IiPETImJGCE3TaGpqQnJyMhISEnDy5EkMGzaMXaFpm4kG6Ddxn8lBHD16dK9JVmAas319fTFw4MBu/WxPrfv6hE8CdujQIbzxxhuIi4vDvHnzOKmjuxAR4ydExHoBzEy0+Ph4pKamQigUIjw8HJGRkQgICNAqaLpM3JdKpcjLy+s1yfo0TePmzZu4deuWThqz+ZC6zxcBA4CkpCSsXr0aBw4cwIIFCziroys0NzejpKQEAODr64t///vfCAwMhL29PUaMGMFxdQSAiFivo6Wlpd1MNFtbWyxYsAARERGYMmWK1kNyhULBClpDQ0O3EvcNkYNoSJhkkaqqKvj7+3cpoLk7dGbdFwqFenOV8knAjhw5gpUrV+Lnn3/GwoULOaujq6SnpyMwMPCB15cvX459+/YZviDCAxAR68W0trayM9FSUlLQv39/LFiwgJ2Jps21qFQq2REyHRP3LS0t222DcZWDqC+YaQFSqRT+/v56TxbRlrrP5AzqyoSjKWD+/v6ctgOcOHECS5cuxQ8//IAXXniBszoIvQsiYn0EhUKBs2fPsjPRBAIBOxNt5syZWp/OVSoVa1TQTNwXCoWQSCSorq7uNaNhaJpGYWEh7ty5A39/f4NPC9BH6j6fBCwtLQ0vvPACoqOjsXTpUl6YXAi9AyJifRCVSoWMjAx2JlpbWxs7Ey0wMFCryUMzcV8qlQIABg8ejGHDhsHGxsaob0oURaGgoAAymYw30ViPa93nk4CdO3cOzz77LHbu3ImVK1ca9WeFwD+IiPVx1Go1zp8/z46QuXv3LkJCQiASiR6YiaZSqVBQUIDm5ma4uLjgzp07qK2thZmZGXuD1eVAREPAl56ph9Fd6z6fBOz3339HZGQkPv/8c6xZs8aoPhv6gKKodkYemqb7/DV5XHqNiO3evRtRUVGorq7G+PHjsWPHDsyYMYPrsowKiqJw6dIlVtCkUimCg4NZU8iSJUvw3HPP4dVXX2VvjBRFsQMRpVIpJ867nqJWq3H16lUolUr4+fkZRXzUo6z7NE0jNzcXcrmccwG7fPkyRCIR/vnPf2LdunV9/matUqlgZmaG5uZm5OXlYdKkSUYdx8YXeoWIHTx4EMuWLcPu3bsxffp0fPvtt9izZw8KCwuJDbaHUBSFrKwsxMXF4bfffkNtbS2GDRuGdevWITIyUusWIuO8Y7Yc+Zy4rzlw1FjntXWcqCwQCGBiYgKBQIBJkyZxGoGVnZ2NsLAwbNq0Ce+++26fFzAmPquhoQFz5szBzJkzsWLFCvj6+pLV2GPSK0TsySefhJ+fH6Kjo9nXxo4di4iICGzbto3DyoyfiooKBAcHY+TIkfD390dycjKKi4sxe/ZsiEQihIWFwc7O7qGJ+1KpFCqVijeJ+0qlEtnZ2TA1NcXEiRN5k833OKhUKmRlZeHevXsQCASgKErv1v3OyMvLw/z58/Huu+/igw8+IDfo/4cZR8T0mw0aNKjdg13HrUZC1zB6EVMoFHjiiSdw6NAhREZGsq+vX78eOTk5yMjI4LA644amafj7+2Py5MnYtWsXTE1N2QGXzEy0goICdiZaWFgYhEJhp4n7TFoIl4n7CoUCWVlZ6N+/PyZMmNArBIyiqHZbiGZmZnq37ndGYWEh5s+fj7Vr1+Ljjz/uswJWVlYGa2trODg4sK9t2bIFJ06cwLlz5wDcz+Q8e/Ysbt++jeXLl2Pw4MFEyHqA0YtYVVUVhg4digsXLrRLwN66dStiYmJw/fp1Dqszfm7duoVhw4ZpvRkxjcHMTLTs7GxMmzYNIpEI4eHhGDx4sFZBa25uZrccW1tb2cR9R0dHvd5g29raIBaLYWlpCW9v715xs+goYB2vnz6s+53x119/ISQkBCtWrMDWrVv7rIApFAp4enrigw8+wOrVq9nXd+zYgYMHD+LXX39FcnIyLl68iLNnz2LkyJFoaGiAWCzuFf2WhqbXiNjFixcxdepU9vVPP/0UP//8M4qKijisru9A0zTKy8tZQfvjjz/w5JNPIjw8vNOZaMD9WB9mhcYk7jM3WF0Kmlwuh1gshq2tLcaNG9cnBEwb+krdv3nzJubNm4fnn38e27dv7xXX93GoqanBoEGD2Mns9vb2SEhIwLZt21BVVQVbW1u8+uqrCA0Nxa1bt7B+/XokJyc/MMmC8GiM7zS7A8x2VE1NTbvXpVIpnJ2dOaqq7yEQCODi4oJ3330X77zzDiorK9mZaJs2bYKvry8iIiIgEokwatQoVtCsrKxgZWWF0aNH4969e5BIJLh9+zauXbsGOzs7doTM46wY7t27B7FY3KvGw/REwADAwsICI0eOxMiRI9tZ90tKSnqcul9eXs72GRIBuw9z71m8eDFyc3ORkZGBhQsXYuTIkaiursbMmTNhaWkJU1NT5OXlgaIoqNVqjqs2Tox+JQbcN3b4+/tj9+7d7GvM3C1i7OAWZiZaYmIiOxNt/PjxrKBpm4kG6C5xv6WlBWKxGM7Ozp3+LmOjpwL2MBjrfm1tLerq6rqcul9ZWYng4GAEBwcjOjqacwHjutWm4xDP69evQyQSYeDAgUhISGgXkC2RSHDp0iUsX74cmzdvxttvv22wOnsTvULEGIv9N998g6lTp+K7777D999/j4KCgl4zFqQ3QNM06uvrcfjwYcTHx+PMmTNwd3dvNxNN282yp4n7MpkMYrG4V80304eAdUSbdZ9pldDs/aupqcHcuXPx1FNPYc+ePZybZPjUavPee+9h7dq1cHV1RUVFBebOnQsrKyvExcWxq7GtW7ciMzMTL730EjZu3GjQ+noTvULEgPtPYF988QWqq6vh5eWFL7/8EjNnzuS6LEInMBb8lJQUxMfH4+TJkxgxYgQ7E23ChAlan+o7Ju5bWVmxgqYZ2NvU1ISsrCy4uLhg1KhRhnxresMQAqbtd2q2ShQVFeH48eMICgrC3r174efnh59++olzAQP402pz8eJFzJ49G+Xl5ey2YmVlJebNmwdzc3MkJCTAxcUFly5dQltbG+8nWvOdXiNiBOPm7t277Ey048ePw8nJiZ2J5u/vr1XQOkvct7CwQFFREVxdXXvNSpyiKFy9ehUKhYKzdBGappGfn49du3bhxIkTqK+vR2hoKBYtWoSwsDBOnXVcttp03EKUy+Xw8fFBdHQ0Zs+ezX69pqYGCxYsgEwmQ3JyMsaMGaO3mvoS5ASWwAtsbGzw4osvIi4uDhKJBF988QUkEgkWLFiAcePGYcOGDbh48WK7w29zc3MMGTIEvr6+mDVrFlxcXNDQ0ICCggKYmppCqVTi7t27MPbnND4IGHDfvDNs2DDk5eVhypQpEIvFmDx5Mnbu3ImRI0dCLpdzUhcA1NXVQa1WP2DmcnZ2fsD0pWsYASsvLwcA9OvXDyYmJrh8+XK7rw8aNAjHjh1DW1sbTp48qdea+hJExPRIZmYmFixYgCFDhkAgECApKYnrkowCS0tLPPvss9i/fz8kEgm++uor3L17F88//zw8PDzw9ttvIzMzEyqViv0ZMzMzUBQFmUyGcePGwdPTE/fu3cOVK1dw4cIF/PXXX2hqajI6QeOLgAH3t2hFIhGGDBmC3377DT4+Pti0aRPEYjFKS0t5kf6vrS/REGehmzZtgqenJwICAjBz5kzY29vj4sWLSE5ORlVVFWQyGQBAKBSiqKgIb7zxht5r6isYvcWez7S0tMDHxwcrV67EokWLuC7HKLGwsIBIJIJIJIJCocDp06cRHx+PZcuWQSAQICwsDJGRkSgvL0dUVBROnjyJoUOHArj/FK5Wq1FfXw+pVIqsrCyYmpqyZ2h8T9znk4DJZDJERkbC3t4e8fHxD7Q8ODo6clTZf38/l602zOewtLQUly9fxuXLl5GSkoLCwkIoFArI5XIMGTIEO3fuJGf1OoaciRkIgUCAxMREREREcF1Kr0CpVLIz0fbv34/W1lYsWLAAL730EmbPnq21r8yYEvf5JGAtLS1YuHAhzMzMcOTIEb1PvO4phmq16Uo01I8//oitW7fiwoULqKqqQn5+PpRKJVauXKmzOgj3ISsxglFibm6OoKAglJaW4tdff8W2bdtQVlaGt956CzKZjJ2JFhQUxCZRmJiYwNHREY6OjvD09GQT9/Pz83mVuM8nAWttbcXzzz8PAEhOTuatgAHAO++8g2XLliEgIIBttamoqMCaNWt09js0TRzHjx+HXC6HQCCASCQCcP/hytzcHNOnT4dSqYRcLsfEiRMxceJEndVAaA9ZiRkIshLTPVevXsXMmTNx5MgRtqGVmYnGBBTX1dWxM9GYXp2O0DSNpqYmNs+Ry8R9PgmYXC7H4sWLcffuXZw4cQK2trac1dJV9Nlqo3m+9sYbb+DkyZMQCAQQCATw8fHBwYMH2e+tr6+Hu7s7du/ejcWLF+vk9xO0Q0TMQBAR0w9MRp02KIqCWCxmBe327dsICgpCREQEQkJCtN6UH5a47+DgoNe5Y3wSsLa2NixduhQSiQSnTp2CnZ0dZ7XwjY8//hh79+7FsWPH4O3tjQ8++ACff/45goKCcPToUfbvzcfHB+vWrcMrr7zCccW9GyJiBoKIGLcwjcKMoJWUlGDOnDkQiUQIDQ3tdCaaoRL3+SRgSqUSL7/8MsrLy3HmzJl240T6OiUlJdiwYQPWrFmD4OBgxMbG4rXXXsP//u//Ijo6Gt7e3khMTMSAAQNw9OhRhIaGcl1yr4eImIEgIsYfNGeiJSQkoLCwEE8//TQ7E83R0bFbiftCoRD9+vXrcT18EjCVSoVVq1bh2rVrOHv2LJycnDirhY9QFIWff/4ZYWFhKCgowJIlS/DZZ59hyZIl2LBhA7Zv3w43Nzfk5+c/1meC0HWIiOmR5uZmlJSUAAA7zTUwMBD29vYGz3IjaIeZicYIWk5ODqZPn87ORBs0aJBWQWMS96VSKWQyWY8T9/kkYGq1GmvWrEFWVhbS0tI63abt7TBnX/X19e1WoR17zrZs2YLCwkLs3bsXFhYW+Oqrr1BYWAhra2t88cUXXJTeJyEipkfS09MRGBj4wOvLly/Hvn37DF8Q4aHQNI2ysjJ2Jtqff/6JKVOmsDPROhsO2nFGV1cT9/kmYOvWrcOFCxeQnp7O9tr1VfLy8hASEoLLly9j8ODBWr9n+fLlyMnJwdWrV9HU1ITnn38e8+fPx/r16w1cbd+GiBiBoAWapnH79m0kJCQgISEBFy5cgJ+fHztCxsXFpdPEfSbPkUncd3JygrOzc7vEfT4JGEVRePvtt3HmzBmkpaX1mrzJx+HixYt44YUXkJOT0+mZ4B9//IGwsDBYWVlhwIABsLS0xJUrVwxcKYGIGIHwCGiaRk1NDTsTLTMzE15eXqygubu7axW0zhL3hUIhiouLeSNgGzduREpKCtLT0zF69GjOauETCoUC7u7uiIqKYvvktH1Pfn4+kpKSYGdnR+aBcQQRMQKhG9A0jbq6OnYm2tmzZzFmzBh2hExnM9E0E/fr6upgYmKCESNGYPDgwbC0tOQk/oqiKPzjH//Ab7/9hvT0dLi7uxu8Bj7QMYGDoigolUrMmjULkZGR2LBhA4fVER4Ff3J2CAZj27ZtmDRpEqytreHk5ISIiAhcv36d67KMAmZA5CuvvIJjx46hpqYG7733HnJzczFjxgz4+/vjk08+wdWrV0FRFPtz5ubmbIafjY0NPDw8cO/ePfzxxx+4ePEiSkpKDJq4T9M0tmzZgv379+P06dN9VsCA+0kuJSUlOHToEG7fvg2ZTIb+/fsjMjIS2dnZANBueoImZA3APWQl1geZN28eFi9ejEmTJkGlUmHTpk3Iy8tDYWEhr2OF+M7du3dx5MgRdibaoEGD2JloHh4eePHFF/Hss89iyZIl7BaiWq1GXV0dO0W5X79+rCnE1tZWLys0mqbx+eefIzo6GmfPnoW3t7fOf4exQNM07t27h/DwcGRlZUEoFKKpqQlTpkxBbm4uBgwYALFYjCeeeOKBuWEEfkBEjIDa2lo4OTkhIyODJGzriObmZqSmpiI+Ph5Hjx6FlZUVLC0t8a9//QuzZ8/WejNUq9VoaGiARCJBbW0tTE1NWVOIrhL3aZrGjh078K9//QtnzpyBr6/vY/+ZvYH6+nrY2tqisLAQV65cQUNDA86ePYuysjJ4e3tjz549sLa2JkLGQ4iIEVBSUgJ3d3fk5eXBy8uL63J6FXK5HBERESgtLYWfnx+OHz+OJ554AgsWLEBERASmTZumNcpKH4n7NE1j165d2LZtG06cOIHJkyfr4i32CrQl07e1tSE+Ph47d+7EkCFDEBMTAxsbmy6l2BMMBxGxPg5N0xCJRLhz5w7OnTvHdTm9jkWLFqGyspIN0JXL5Thz5gwSEhJw+PBhmJiYsII2c+ZMrU5FiqLQ2NjIpoUwiftMQHFXbqg0TeP777/H5s2bkZqaimnTpunj7Ro9TEMzI1QKhQIHDx7Ed999B5qmcfToUaMIQu5LEBHr47z++us4evQozp8/j2HDhnFdTq/j/Pnz8Pb21nrjUyqVSE9PR3x8PJKSkqBUKhEWFgaRSITAwECtyR89SdynaRoxMTF4//33ceTIEd5uGX/66ac4evQocnJy0K9fPzQ2NnJaDyNoarUa+/btw7Fjx7Bnzx4ShswziIj1YdatW4ekpCRkZmZi1KhRXJfTp1GpVDh//jzi4uKQlJSE5uZmzJ8/n52JptkozdCVxH2apvHrr7/i3XffRXJystYEGb6wefNmDBw4ELdv38bevXs5FzGg/cqsra1N698DgVuIiPVBaJrGunXrkJiY2Kf7g/iKWq1uNxOtvr4ec+fOZWeiaXOQdkzc3717N+rr6zFu3Dj89NNPSEhIwNy5czl4N91n3759eOutt3ghYsCDmYkEfkFErA/y2muvITY2FocPH4aHhwf7uq2tLXnS5BkUReHKlSusoFVVVbWbiWZjY6P153JycrBz506cOnUKMpkMwcHBWLRoEUQiEe9Hq/BNxAj8hohYH6Szp8off/wRK1asMGwxhC7D5C0ygnbz5s12M9E0bfhHjhzBypUr8csvv8Db2xvx8fGIj4+HjY0NTp8+zfE7eThExAjdgYgYgWCE0DSNwsJCdoTMtWvXMGvWLIhEIlhZWeG1117Djz/++EDun1wuf2iyvq75+OOP8cknnzz0ey5fvoyAgAD2/4mIEboDETECwcihaRrFxcWIi4vDb7/9hqtXr2LXrl1Yu3Yt52c5dXV1qKure+j3uLi4tBNWImKE7kBEjGcwfx1c33wIxglN07hy5QoCAgKM9jNERIzQHUjbOY9oaWmBQCAw2puProiOjsaECRNgY2MDGxsbTJ06FampqVyXZRQIBAJMmjTJKD9DFRUVyMnJQUVFBdRqNXJycpCTk4Pm5mauSyPwGLIS4xGLFy+GlZUVdu3axTa6MskBfcnmm5KSAlNTU7i5uQEAYmJiEBUVhezsbIwfP57j6gj6YsWKFYiJiXng9bS0NMyaNcvwBRGMAiJiPOL8+fMIDQ1FRUUFbG1tWeGqqanBoEGDuC6PU+zt7REVFYVVq1ZxXQqBQOARZDuRRwwZMgTDhg3DsWPHANxPQt+5cydGjx6Nb775huPquEGtVuPAgQNoaWnB1KlTuS6HQCDwjAfjswmcQFEURo8ejYEDB6KwsBCNjY14+eWX8ddff2H79u1Ys2YNgL6THpCXl4epU6dCLpfDysoKiYmJGDduHNdlEQgEnkG2E3lGbGwstmzZAplMhhEjRuDrr79mZz51NgKiNzoaFQoFKioq0NjYiPj4eOzZswcZGRlEyAgEQjvIdiIPYMbYV1dXo6CgAEVFRZg1axaOHDnSbmhhZyM3GEcjRVGdjlE3Nvr16wc3NzcEBARg27Zt8PHxwc6dO7kui0Ag8AyyncgxzKTY8vJyPPfcc+yYh4CAANjZ2UGpVGqdMQXcTz4/deoUpFIpgoKCMHToUEOWblBomkZbWxvXZRAIBJ5BVmIcY2pqiszMTMyYMQPW1tb48ccfsXr1apw9exZqtbpTAQOAGzdu4MKFC9i7dy/Gjh2LOXPm4PLlyw9838N2jPm4m/zhhx/i3LlzKCsrQ15eHjZt2oT09HQsWbKE69IIBALPICLGISqVCq+++iqWLl2KkJAQpKamYsiQIVi4cCEyMjLQ0tLy0J8fM2YM/vGPfyAzMxM3btzAgAEDEB0dDYVCAQBoampCdXU1u9XIQFFUu3O0mJgYKJVK/b3RbiKRSLBs2TJ4eHhgzpw5+OOPP3D8+HE888wzXJdGIBB4BjF2cIharUZ0dDQGDx6MyMhItqm5rKwM8+fPx+bNm7F48WKtP0tRFG7evImMjAyMHj0agYGByM/PxzPPPIO0tDR4enqiuLgYc+bMwaeffoply5Zp/XMkEgmeeuop7N27l7cTfwkEAqEziIjxDMZCP23aNPj4+CA6Olqrrf7LL7/EF198AXd3d9TU1ODOnTvw9PTEn3/+idbWVtYE8txzz8HS0hL79u0DcD+Xrrq6Gq+99hqsra3Z7+vM+UggEAh8hty1eAYjVv/85z9hamra7jUGuVyOL7/8Eq+++iri4uIgFouRlJSEyspKzJ07FzRNs9uHK1aswLFjx1BcXIzNmzdj7dq1qKysBHDf7fjpp5+iuroaJiYm7bYcVSpVu/8nEAgEPkLciTwlKCgIQUFBWr8ml8vR2NiIIUOGwMnJCQAgFArR0tKCp59+GqampqwAjRkzBmPHjkVERATMzc3x66+/YuHChQCAgoIC/OMf/0BAQAAGDx4MExMTSKVSODk5wcyMfDQIBAL/ISsxI8TS0hLr16/Hpk2b8MEHH+A///kPgoODIZfLERISAuC/PWXx8fE4d+4cAGD//v1YuHAh20v2yy+/YNy4cXjyySfR3NyM7777DgsXLsSIESOwbNky/Pnnn9y8QSNj27ZtEAgEeOutt7guhUDocxARM0LMzc3xySefYOfOnSgoKICpqSmcnJzg4uKCUaNGAQCqqqqwdu1a/Pvf/2aHI44dOxbAf7cnDxw4gODgYAwcOBApKSmIiopCcHAw9u7dC4qisGHDBqSlpXH2Po2By5cv47vvvsOECRO4LsUglJWVYdWqVRg1ahQsLCzg6uqKzZs3s45YAsHQkD0jI8XExAQvvfQSXnrpJQCAt7c3amtrYWFhgZiYGHz66acYNGgQEhMTIRQKkZKSgszMTMycORMmJia4efMmysvL2ZWbnZ0daJrGkiVL4OrqimeeeQapqakYPXo0l2+T1zQ3N2PJkiX4/vvvsWXLFq7LMQhFRUWgKArffvst3NzckJ+fj9WrV6OlpQXbt2/nujxCH4S4E3shly5dQlpaGlavXg1HR0cAgJ+fHxYtWoRNmzYBAD766CMcPHgQZ86cwbBhw1BZWYmlS5fizp07ePPNN7F8+XLWWELQzvLly2Fvb48vv/wSs2bNwsSJE7Fjxw6uyzI4UVFRiI6Oxs2bN7kuhdAHISuxXsiUKVMwZcqUdq89//zz2LNnD/72t79BKBTi4MGDCA0NZeeUDR06FMeOHcPXX3+Nb775BnV1ddiwYQOx3nfCgQMHkJWVpTUhpa/R1NQEe3t7rssg9FHI3amP8P777+O3336DUCjErVu3UFxcjFmzZsHMzAz79u1DYWEhLCws8N577+Fvf/sboqKiIBaLiYBp4datW1i/fj1++eUXDBgwgOtyOOXGjRv46quv2FFBBIKhIXeoPoSfnx/UajWGDx+OnJwcBAYGorm5GdnZ2di4cSMSEhIglUqhUChQX1/f56dJd4ZYLIZUKoW/vz/MzMxgZmaGjIwM/Oc//4GZmZlRThL4+OOP2WkInf135cqVdj9TVVWFefPm4bnnnsMrr7zCUeWEvg45EyOgrKwMX331FWJjY2FiYgJ3d3dMnjwZX3zxBdlO1IJMJkN5eXm711auXAlPT09s3LgRXl5eHFXWc+rq6lBXV/fQ73FxcWFXnlVVVQgMDMSTTz6Jffv2kc8IgTOIiBHacf36dfTv3x8uLi5cl2JU9CVjR2VlJQIDA+Hv749ffvmFGIAInEKMHYR2eHh4cF0CgcdUVVVh1qxZGDFiBLZv347a2lr2a2T7mcAFZCVGIBC6zL59+7By5UqtXyO3EgIXEBEjEAgEgtFCTmMJBAKBYLQQESMQCASC0UJEjEAgEAhGCxExAoFAIBgtRMQIBAKBYLQQESMQCASC0UJEjEAgEAhGCxExAoFAIBgtRMQIBAKBYLQQESMQCASC0UJEjEAgEAhGCxExAoFAIBgt/wezobxr8XdohgAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "#| caption: Result of the IMU preintegration optimization.\n",
    "#| label: fig:imu_result\n",
    "i = 0\n",
    "while imu_result.exists(X(i)):\n",
    "    gtsam_plot.plot_pose3(1, imu_result.atPose3(X(i)), 1)\n",
    "    i += 1\n",
    "plt.title(\"imu_result\")\n",
    "gtsam.utils.plot.set_axes_equal(1)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Incorporating Visual Measurements\n",
    "\n",
    "> Structure from Motion is 3D reconstruction from images\n",
    "\n",
    "In this section we provide an introduction to reconstructing a 3D environment from camera images over time, which allows for planning and obstacle avoidance, which we discuss in more detail in the next section. However, mapping the environment is in and of itself a useful capability in many applications. In some industrial and military applications, the map *is* the product of interest. In this section we discuss 3D reconstruction as a *batch* process, done after the drone has collected images."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "0a8e2d5d",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "![Balbianello](https://github.com/gtbook/robotics/blob/main/Figures7/Balbianello.jpg?raw=1 \"Balbianello\")\n",
    "\n",
    "Structure from Motion, or SfM, is the problem of computing the 3D structure of a scene from multiple images, either taken from different cameras, or taken at different times as the camera moves in the environment. In the figure above we show one such image sequence, five images of the \"Balbianello\" villa in Italy, situated on lake Como. "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "b9715f32",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "```{index} Bundle Adjustment\n",
    "```\n",
    "SfM is almost the same problem as Visual SLAM, but typically there is no odometry information to relate\n",
    "different camera views. The resulting optimization problem is sometimes called **Bundle Adjustment**. \n",
    "The most difficult problems associated with SfM are feature detection, data association, and initializing the optimization.\n",
    "\n",
    "Below we load the Balbianello data from a file and create a factor graph, as well as an initial estimate. A lot of work is done in the code snippet below! The particular file format, which stems from the \"Bundler\"\" library, is transformed into a data structure `sfm_data`. The method `sfmFactorGraph` then creates a factor graph with *all* the projection factors derived from the close to 1500 visual measurements. The `initialCamerasAndPointsEstimate` methods next instantiates a `Valeus` object with initial estimates for all the cameras and the points, extracted from the Bundler file. We then print the initial error and see that it has a value of $126.9$, which is normalized with respect to 1-pixel standard deviation given in the noise model."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "ff7d45ce",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "initial error: 126.9\n"
     ]
    }
   ],
   "source": [
    "filename = gtsam.findExampleDataFile(\"Balbianello.out\")\n",
    "sfm_data = gtsam.SfmData.FromBundlerFile(filename)\n",
    "model = gtsam.noiseModel.Isotropic.Sigma(2, 1.0)  # one pixel in u and v\n",
    "graph = sfm_data.sfmFactorGraph(model)\n",
    "initial = gtsam.initialCamerasAndPointsEstimate(sfm_data)\n",
    "print(f\"initial error: {graph.error(initial):.1f}\")\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The code above already did most of the heavy lifting that was so cumbersome and manual in earlier sections of this book, but it does not yet optimize! We once again use the Levenberg-Marquardt non-linear optimizer below, which is now performing *structure from motion*: we optimize for both the unknown camera poses (calibration is assumed known here) and the unknown 3D points in the scene:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "704cbe24",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "final error: 125.5\n"
     ]
    }
   ],
   "source": [
    "try:\n",
    "    lm = gtsam.LevenbergMarquardtOptimizer(graph, initial)\n",
    "    result = lm.optimize()\n",
    "except RuntimeError:\n",
    "    print(\"LM Optimization failed\")\n",
    "print(f\"final error: {graph.error(result):.1f}\")\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "24491487",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "The final error is not much smaller than the initial error, which tells us that the initial estimate was already very good. It is quite fortunate to have found that initial estimate in the file, because in many cases initializing a structure for motion problem can be quite a difficult undertaking. When we are flying a drone, we can build up an estimate of the 3D environment gradually, and hence an initial estimate is always available to use from the previous time step. However, we will not discuss this more in depth here.\n",
    "\n",
    "The result of SfM is a sparse point cloud and a set of well-determined camera poses, which can be \"densified\" using multi-view stereo techniques. This is out of scope of this book, however, and below we simply extract the 3D points and display them in 3D with `plotly`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "8fa3f485",
   "metadata": {},
   "outputs": [],
   "source": [
    "structure = gtsam.utilities.extractPoint3(result)\n",
    "valid = structure[:,2] > -5.0\n",
    "bRc = sfm_data.camera(0).pose().rotation()\n",
    "wRc = gtsam.Rot3([1,0,0],[0,0,-1],[0,1,0]) # make camera point to world y\n",
    "wRb = wRc.compose(bRc.inverse())\n",
    "structure = (wRb.matrix() @ structure.T).T"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#| caption: Structure from motion result as a sparse 3D point cloud.\n",
    "#| label: fig:sfm_result\n",
    "fig = go.Figure()\n",
    "fig.add_trace(go.Scatter3d(x=structure[valid,0],y=structure[valid,1],z=structure[valid,2],\n",
    "    mode=\"markers\", marker=dict(size=2, color=\"blue\"), showlegend= False))\n",
    "fig.update_layout(margin=dict(l=0, r=0, t=0, b=0))\n",
    "fig.update_yaxes(scaleanchor = \"x\", scaleratio = 1); fig.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Visual SLAM\n",
    "\n",
    "> Visual SLAM is based on Structure from Motion\n",
    "\n",
    "Visual SLAM combines both IMU integration and SfM in a real-time, incremental pipeline. Typically, The drone starts with no pre-built map, opens its eyes as it were, and starts building a 3D map of the environment incrementally. Initial estimates for the camera positions are given by integrating the IMU forward. New landmarks are added incrementally as the drone progresses throughout the environment. Fast, incremental optimization techniques can be used to keep the computation under control, but discussing these is out of scope for an introductory textbook. Suffice to say that, modulo a *lot* of engineering efforts, this process can be made to work incredibly well. When done well, this provides a superpower to autonomous drones, allowing them to map their environment in real time and affording it to do obstacle avoidance and planning in hitherto unknown environments.\n",
    "\n",
    "Visual SLAM from images and IMU integration compliment each other perfectly. The IMU excels at dealing with high frequency motion like vibrations, which is common on drones. However, it suffers from drift: the inertial measurement units available for consumer drones can drift away from the true trajectory quite rapidly. High quality IMUs are typically very expensive and only available in high-end applications. In fact, this technology is so sensitive that even *disclosing* the accuracy of military-grade IMUs is a crime in many countries. Visual SLAM, on the other hand, is not great at dealing with high frequency motions, but does not suffer from drift as much, especially in environments where many loop closures can be made. It is only in very long, open-ended trajectories that drift is again be a problem."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Summary\n",
    "\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In summary, above we introduced Visual SLAM as an extension of SLAM with landmarks and its critical integration with an IMU integration process. We discussed how the latter can be done with factor graphs, and while these are some of the most complex factor graphs we encountered, the principle stays the same. Structure from Motion was introduced as a method for constructing 3D maps from images, and while we did not go into too much detail about the specific factors, the same optimization pattern can be found at the center of it. Visual SLAM was then introduced, at a high level, as a combination of IMU integration and SfM."
   ]
  }
 ],
 "metadata": {
  "colab": {
   "collapsed_sections": [],
   "include_colab_link": true,
   "name": "S74_drone_perception.ipynb",
   "provenance": []
  },
  "kernelspec": {
   "display_name": "Python 3.8.12 ('gtbook')",
   "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.9.18"
  },
  "latex_metadata": {
   "affiliation": "Georgia Institute of Technology",
   "author": "Frank Dellaert and Seth Hutchinson",
   "title": "Introduction to Robotics"
  },
  "vscode": {
   "interpreter": {
    "hash": "9f7376ced4243bb13dfcffa8a3ba834e0602aa8334cd3a1d8ba8d285f4628083"
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}