examples/Two_Moons_Starter.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "009b6adf",
   "metadata": {},
   "source": [
    "# Two Moons: Tackling Bimodal Posteriors\n",
    "\n",
    "_Authors: Lars Kühmichel, Marvin Schmitt, Valentin Pratz, Stefan T. Radev_"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "d5f88a59",
   "metadata": {
    "ExecuteTime": {
     "end_time": "2025-04-11T19:54:01.403328Z",
     "start_time": "2025-04-11T19:53:24.823026Z"
    }
   },
   "outputs": [],
   "source": [
    "import os\n",
    "# Set to your favorite backend\n",
    "os.environ[\"KERAS_BACKEND\"] = \"jax\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "0551e46f",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING:2025-04-21 12:20:32,803:jax._src.xla_bridge:967: An NVIDIA GPU may be present on this machine, but a CUDA-enabled jaxlib is not installed. Falling back to cpu.\n"
     ]
    }
   ],
   "source": [
    "import numpy as np\n",
    "\n",
    "import matplotlib.pyplot as plt\n",
    "import seaborn as sns\n",
    "\n",
    "import bayesflow as bf"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c63b26ba",
   "metadata": {},
   "source": [
    "## Simulator<a class=\"anchor\" id=\"simulator\"></a>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "9525ffd7",
   "metadata": {},
   "source": [
    "This example will demonstrate amortized estimation of a somewhat strange Bayesian model, whose posterior evaluated at the origin $x = (0, 0)$ of the \"data\" will resemble two crescent moons. The forward process is a noisy non-linear transformation on a 2D plane:\n",
    "\n",
    "$$\n",
    "\\begin{align}\n",
    "x_1 &= -|\\theta_1 + \\theta_2|/\\sqrt{2} + r \\cos(\\alpha) + 0.25\\\\\n",
    "x_2 &= (-\\theta_1 + \\theta_2)/\\sqrt{2} + r\\sin{\\alpha}\n",
    "\\end{align}\n",
    "$$\n",
    "\n",
    "with $x = (x_1, x_2)$ playing the role of \"observables\" (data to be learned from), $\\alpha \\sim \\text{Uniform}(-\\pi/2, \\pi/2)$, and $r \\sim \\text{Normal}(0.1, 0.01)$ being latent variables creating noise in the data, and $\\theta = (\\theta_1, \\theta_2)$ being the parameters that we will later seek to infer from new $x$. We set their priors to\n",
    "\n",
    "$$\n",
    "\\begin{align}\n",
    "\\theta_1, \\theta_2 \\sim \\text{Uniform}(-1, 1).\n",
    "\\end{align}\n",
    "$$\n",
    "\n",
    "This model is typically used for benchmarking simulation-based inference (SBI) methods (see https://arxiv.org/pdf/2101.04653) and any method for amortized Bayesian inference should be capable of recovering the two moons posterior *without* using a gazillion of simulations. Note, that this is a considerably harder task than modeling the common unconditional two moons data set used often in the context of normalizing flows. However, free-form models (e.g., flow matching, diffusion) will tend to outperform normalizing flows on multimodal data sets."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "21bf228e706a010",
   "metadata": {},
   "source": [
    "BayesFlow offers many ways to define your data generating process. Here, we use sequential functions to build a simulator object for online training. Within this composite simulator, each function has access to the outputs of the previous functions. This effectively allows you to define any generative graph."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "f761b142a0e1da66",
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-10-24T08:36:22.160542Z",
     "start_time": "2024-10-24T08:36:22.155240Z"
    }
   },
   "outputs": [],
   "source": [
    "def theta_prior():\n",
    "    theta = np.random.uniform(-1, 1, 2)\n",
    "    return dict(theta=theta)\n",
    "\n",
    "def forward_model(theta):\n",
    "    alpha = np.random.uniform(-np.pi / 2, np.pi / 2)\n",
    "    r = np.random.normal(0.1, 0.01)\n",
    "    x1 = -np.abs(theta[0] + theta[1]) / np.sqrt(2) + r * np.cos(alpha) + 0.25\n",
    "    x2 = (-theta[0] + theta[1]) / np.sqrt(2) + r * np.sin(alpha)\n",
    "    return dict(x=np.array([x1, x2]))"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "722cb773",
   "metadata": {},
   "source": [
    "Within the composite simulator, every simulator has access to the outputs of the previous simulators in the list. For example, the last simulator `forward_model` has access to the outputs of the three other simulators."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "4b89c861527c13b8",
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-10-24T08:36:22.305265Z",
     "start_time": "2024-10-24T08:36:22.301546Z"
    }
   },
   "outputs": [],
   "source": [
    "simulator = bf.make_simulator([theta_prior, forward_model])"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f6e1eb5777c59eba",
   "metadata": {},
   "source": [
    "Let's generate some data to see what the simulator does:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "e6218e61d529e357",
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-10-24T08:36:22.350483Z",
     "start_time": "2024-10-24T08:36:22.345161Z"
    }
   },
   "outputs": [],
   "source": [
    "# generate 3 random draws from the joint distribution p(r, alpha, theta, x)\n",
    "sample_data = simulator.sample(3)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "46174ccb0167026c",
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-10-24T08:36:22.470435Z",
     "start_time": "2024-10-24T08:36:22.464836Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Type of sample_data:\n",
      "\t <class 'dict'>\n",
      "Keys of sample_data:\n",
      "\t dict_keys(['theta', 'x'])\n",
      "Types of sample_data values:\n",
      "\t {'theta': <class 'numpy.ndarray'>, 'x': <class 'numpy.ndarray'>}\n",
      "Shapes of sample_data values:\n",
      "\t {'theta': (3, 2), 'x': (3, 2)}\n"
     ]
    }
   ],
   "source": [
    "print(\"Type of sample_data:\\n\\t\", type(sample_data))\n",
    "print(\"Keys of sample_data:\\n\\t\", sample_data.keys())\n",
    "print(\"Types of sample_data values:\\n\\t\", {k: type(v) for k, v in sample_data.items()})\n",
    "print(\"Shapes of sample_data values:\\n\\t\", {k: v.shape for k, v in sample_data.items()})"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "17f158bd2d7abf75",
   "metadata": {},
   "source": [
    "BayesFlow also provides this simulator and a collection of others in the `bayesflow.benchmarks` module."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f714c3a178b5a375",
   "metadata": {},
   "source": [
    "## Adapter\n",
    "\n",
    "The next step is to tell BayesFlow how to deal with all the simulated variables. You may also think of this as informing BayesFlow about the data flow, i.e., which variables go into which network and what transformations needs to be performed prior to passing the simulator outputs into the networks. This is done via an adapter layer, which is implemented as a sequence of fixed, pseudo-invertible data transforms.\n",
    "\n",
    "Below, we define the data adapter by specifying the input and output keys and the transformations to be applied. This allows us full control over the data flow."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "5c9c2dc70f53d103",
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-10-24T08:36:26.618926Z",
     "start_time": "2024-10-24T08:36:26.614443Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Adapter([0: ToArray -> 1: ConvertDType -> 2: Standardize(exclude=['theta']) -> 3: Rename('theta' -> 'inference_variables') -> 4: Rename('x' -> 'inference_conditions')])"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "adapter = (\n",
    "    bf.adapters.Adapter()\n",
    "    \n",
    "    # convert any non-arrays to numpy arrays\n",
    "    .to_array()\n",
    "    \n",
    "    # convert from numpy's default float64 to deep learning friendly float32\n",
    "    .convert_dtype(\"float64\", \"float32\")\n",
    "    \n",
    "    # standardize target variables to zero mean and unit variance \n",
    "    .standardize(exclude=\"theta\")\n",
    "    \n",
    "    # rename the variables to match the required approximator inputs\n",
    "    .rename(\"theta\", \"inference_variables\")\n",
    "    .rename(\"x\", \"inference_conditions\")\n",
    ")\n",
    "adapter"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "254e287b2bccdad",
   "metadata": {},
   "source": [
    "## Dataset\n",
    "\n",
    "For this example, we will sample our training data ahead of time and use offline training with a very small number of epochs. In actual applications, you usually want to train much longer in order to max our performance."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "39cb5a1c9824246f",
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-09-23T14:39:46.950573Z",
     "start_time": "2024-09-23T14:39:46.948624Z"
    }
   },
   "outputs": [],
   "source": [
    "num_training_batches = 512\n",
    "num_validation_sets = 300\n",
    "batch_size = 64\n",
    "epochs = 50"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "9dee7252ef99affa",
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-09-23T14:39:53.268860Z",
     "start_time": "2024-09-23T14:39:46.994697Z"
    }
   },
   "outputs": [],
   "source": [
    "training_data = simulator.sample(num_training_batches * batch_size)\n",
    "validation_data = simulator.sample(num_validation_sets)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "2d4c6eb0",
   "metadata": {},
   "source": [
    "## Training a neural network to approximate all posteriors\n",
    "\n",
    "The next step is to set up the neural network that will approximate the posterior $p(\\theta\\,|\\,x)$.\n",
    "\n",
    "We choose **Flow Matching** [1, 2] as the backbone architecture for this example, as it can deal well with the multimodal nature of the posteriors that some observables imply.\n",
    "\n",
    "* [1] Lipman, Y., Chen, R. T., Ben-Hamu, H., Nickel, M., & Le, M. Flow Matching for Generative Modeling. In *The Eleventh International Conference on Learning Representations*.\n",
    "\n",
    "* [2] Wildberger, J. B., Dax, M., Buchholz, S., Green, S. R., Macke, J. H., & Schölkopf, B. Flow Matching for Scalable Simulation-Based Inference. In *Thirty-seventh Conference on Neural Information Processing Systems*."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "09206e6f",
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-09-23T14:39:53.339590Z",
     "start_time": "2024-09-23T14:39:53.319852Z"
    }
   },
   "outputs": [],
   "source": [
    "flow_matching = bf.networks.FlowMatching(\n",
    "    subnet=\"mlp\", \n",
    "    subnet_kwargs={\"dropout\": 0.0, \"widths\": (256,)*6} # override default dropout = 0.05 and widths = (256,)*5\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "851e522f",
   "metadata": {},
   "source": [
    "This inference network is just a general Flow Matching backbone, not yet adapted to the specific inference task at hand (i.e., posterior appproximation). To achieve this adaptation, we combine the network with our data adapter, which together form an `approximator`. In this case, we need a `ContinuousApproximator` since the target we want to approximate is the posterior of the *continuous* parameter vector $\\theta$."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "76722c33",
   "metadata": {},
   "source": [
    "### Basic Workflow\n",
    "We can hide many of the traditional deep learning steps (e.g., specifying a learning rate and an optimizer) within a `Workflow` object. This object just wraps everything together and includes some nice utility functions for training and *in silico* validation."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "96ca6ffa",
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-09-23T14:39:53.371691Z",
     "start_time": "2024-09-23T14:39:53.369375Z"
    }
   },
   "outputs": [],
   "source": [
    "flow_matching_workflow = bf.BasicWorkflow(\n",
    "    simulator=simulator,\n",
    "    adapter=adapter,\n",
    "    inference_network=flow_matching,\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "708b1303",
   "metadata": {},
   "source": [
    "### Training\n",
    "\n",
    "We are ready to train our deep posterior approximator on the two moons example. We use the utility function `fit_offline`, which wraps the approximator's super flexible `fit` method."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "0f496bda",
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-09-23T14:42:36.067393Z",
     "start_time": "2024-09-23T14:39:53.513436Z"
    }
   },
   "outputs": [],
   "source": [
    "history = flow_matching_workflow.fit_offline(\n",
    "    training_data, \n",
    "    epochs=epochs, \n",
    "    batch_size=batch_size, \n",
    "    validation_data=validation_data\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "66248a2f",
   "metadata": {},
   "source": [
    "## Swapping Inference Networks <a class=\"anchor\" id=\"swapping\"></a>\n",
    "\n",
    "Using BayesFlow, it is easy to switch to a different backbone architecture for the inference network. For instance, the code below demonstrates the use of a **Consistency Model**, which can allow for faster sampling during inference."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "80a97232",
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-09-23T14:42:38.639240Z",
     "start_time": "2024-09-23T14:42:38.637439Z"
    }
   },
   "source": [
    "### Consistency Models: Background\n",
    "\n",
    "Consistency Models (CM; [1]) leverage the nice properties of score-based diffusion to enable few-step sampling. Score-based diffusion initially relied on a stochastic differential equation (SDE) for sampling, but there is also a ordinary (non-stochastic) differential equation (ODE)that has the same _marginal_ distribution at each time step $t$ [2]. This means that even though SDE and ODE produce different paths from the noise distribution to the target distribution, the resulting distributions when looking at many paths at time $t$ is the same. The ODE is also called Probability Flow ODE.\n",
    "\n",
    "CMs leverage the fact that there is no randomness in the ODE formulation. That means, if you start at a certain point in the latent space, you will always take the same path and end up at the same point in the target $\\theta$-space. The same is true for every point on the path: if you integrate to get to time $t=0$, you will end up at the same point as well. In short: for each path, there is exactly one corresponding point in latent space (at $t=T$) and one corresponding point in data space (at $t=0$). \n",
    "\n",
    "The goal of CMs is the following: each point at a time point $t$ belongs to exactly one path, and we want to predict where this path will end up at $t=0$. The function that does this is called the _consistency function_ $f$. If we have the correct function for all $t\\in(0,T]$, we can just sample from the latent distribution ($t=T$) and use $f$ to directly map to the corresponding point at $t=0$, which is in the target distribution. So for sampling from the target distribution, we avoid any integration and only need one evaluation of the consistency function. In practice, the one-step sampling does not work very well. Instead, we leverage a multi-step sampling method where we call $f$ multiple times. Please check out the [1] for more background on this sampling procedure.\n",
    "\n",
    "When reading the above, you might wonder why we also learn the mapping to $t=0$ of all intermediate time steps $t\\in[0, T]$, and not only for $t=T$. The main answer is that for efficient training, we do not want to actually compute the two associated points explicitly. Doing so would require to do a precise integration at training time, which is often not feasible as it is too computationally costly. Learning all time steps opens up the possibility for a different training approach where we can avoid this. The details of this become a bit more complicated, and we advise you to take a look at [1] if you are interested in a more thorough and mathematical discussion. Below, we will give a rough description of the underlying concepts.\n",
    "\n",
    "\n",
    "**Training** First, we know that at $t=0$, it holds that $f(\\theta,t=0)=\\theta$, as $\\theta$ is part of the path that ends at $\\theta$. This _boundary condition_ serves as an \"anchor\" for our training, this is the information that the network knows at the start of the training procedure (we encode it with a time-dependent skip-connection, so the network is forced to be the identity function at $t=0$). For training, we now somehow have to propagate this information to the rest of the part. The basic idea for this is simple. We just take a point $\\theta_1$ closer to the data distribution (smaller time $t_1$) and integrate for a small time step $dt$ to a point $\\theta_2$ on the same path that is closer to the latent distribution (larger time $t_2=t_1+dt$). As we know that for $t=0$ our network provides the correct output for our path, we want to propagate the information from smaller times to larger times. Our training goal is to move the output of $f(\\theta_2, t=t_2)$ towards the output of $f(\\theta_1, t=t_1)$. How to choose $\\theta_1$, $t_1$ and $dt$ is an empirical question, see the [1] for some thoughts on what works well.\n",
    "\n",
    "**Distilling inference** In the case of _distillation_, we start with a trained score-based diffusion model. We can use it to integrate the Probability Flow ODE to get from $\\theta_1$ to $\\theta_2$. If we do not have such a model, it seems as if we were stuck. We do not know which points lie on the same path, so we do not know which outputs to make similar. Fortunately, it turns out that there is an _unbiased approximator_ that, if averaged over many samples (check out the paper for the exact description), will also give us the correct score. If we use this approximator instead of the score model, and use only a single Euler step to move along the path, we get an algorithm similar to the one described for distillation. It is called Consistency Training (CT) and allows us to train a consistency model using only _samples_ from the data distribution. The algorithm for this was improved a lot in [3], and we have incorporated those improvements into our implementation.\n",
    "\n",
    "**Improving consistency training** We have made several improvements to get to a standalone *consistency training* algorithm. As a consequence, the introduced hyperparameters and their choice unfortunately becomes somewhat unintuitive. We have to rely on empirical observations and heuristics to see what works. This was done in [4], we encourage you to use the values provided there as starting points. If you happen to find hyperparameters that work significantly better, please let us know (e.g., by opening an issue or sending an email). This will help others to find the correct region in the hyperparameter space.\n",
    "\n",
    "[1] Song, Y., Dhariwal, P., Chen, M., & Sutskever, I. (2023). Consistency Models. _arXiv preprint_. [https://doi.org/10.48550/arXiv.2303.01469](https://doi.org/10.48550/arXiv.2303.01469)\n",
    "\n",
    "[2] Song, Y., Sohl-Dickstein, J., Kingma, D. P., Kumar, A., Ermon, S., & Poole, B. (2021). Score-Based Generative Modeling through Stochastic Differential Equations. In _International Conference on Learning Representations_. [https://openreview.net/forum?id=PxTIG12RRHS](https://openreview.net/forum?id=PxTIG12RRHS)\n",
    "\n",
    "[3] Song, Y., & Dhariwal, P. (2023). Improved Techniques for Training Consistency Models. _arXiv preprint_. [https://doi.org/10.48550/arXiv.2310.14189](https://doi.org/10.48550/arXiv.2310.14189)\n",
    "\n",
    "[4] Schmitt, M., Pratz, V., Köthe, U., Bürkner, P.-C., & Radev, S. T. (2024). Consistency Models for Scalable and Fast Simulation-Based Inference. _arXiv preprint_. [https://doi.org/10.48550/arXiv.2312.05440](https://doi.org/10.48550/arXiv.2312.05440)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "24b69117",
   "metadata": {},
   "source": [
    "### Consistency Models: Specification"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "3016f08f",
   "metadata": {},
   "source": [
    "We can now go ahead and define our new inference network backbone. Apart from the usual parameters like learning rate and batch size, CMs come with a number of different hyperparameters. Unfortunately, they can heavily interact, so they can be hard to tune. The main hyperparameters are:\n",
    "\n",
    "- Maximum time `max_time`: This also serves as the standard deviation of the latent distribution. You can experiment with this, values from 10-200 seem to work well. In any case, it should be larger than the standard deviation of the target distribution.\n",
    "- Minimum/maximum number of discretization steps during training `s0`/`s1`: The effect of those is hard to grasp. 10 works well for `s0`. Intuitively, increasing `s1` along with the number of epochs should lead to better result, but in practice we sometimes observe a breakdown for high values of `s1`. This seems to be problem-dependent, so just try it out.\n",
    "- `sigma2` modifies the time-dependency of the skip connection. Its effect on the training is unclear, we recommend leaving it at 1.0 or setting it to the approximate variance of the target distribution.\n",
    "- Smallest time value `eps` ($t=\\epsilon$ is used instead of $t=0$ for numerical reasons): No large effect in our experiments, as long as it is kept small enough. Probably not worth tuning.\n",
    "\n",
    "You may find that different hyperparameter values work better for your tasks.\n",
    "\n",
    "A short note on dropout: in our experiments, dropout usually lead to worse performance, so generally we recommend setting the droput rate to $0.0$. Consistency training takes advantage of a noisy estimator of the score, so probably the training is already sufficiently noisy and extra dropout for regularization is not necessary."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "d53a41b8",
   "metadata": {},
   "outputs": [],
   "source": [
    "consistency_model = bf.networks.ConsistencyModel(\n",
    "    total_steps=num_training_batches*epochs,\n",
    "    subnet_kwargs={\"dropout\": 0.0, \"widths\": (256,)*6},\n",
    "    max_time=10, # this probably needs to be tuned for a novel application\n",
    "    sigma2=1.0,  # the data adapter standardizes our parameters, so set to 1.0\n",
    ")\n",
    "\n",
    "# Workflow for consistency model\n",
    "consistency_model_workflow = bf.BasicWorkflow(\n",
    "    simulator=simulator,\n",
    "    adapter=adapter,\n",
    "    inference_network=consistency_model,\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "9fbcca16",
   "metadata": {},
   "source": [
    "### Consistency Training"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "c3c1a812",
   "metadata": {},
   "outputs": [],
   "source": [
    "history = consistency_model_workflow.fit_offline(\n",
    "    training_data, \n",
    "    epochs=epochs, \n",
    "    batch_size=batch_size, \n",
    "    validation_data=validation_data\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "a94a43f6",
   "metadata": {},
   "source": [
    "## Good 'ol Coupling Flows\n",
    "\n",
    "Of course, BayesFlow also supports established coupling flow models with a variety of parameters, including the timeless *affine* and *spline* flows."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "69bec2e2",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Use a deeper affine coupling flow (default depth is 6)\n",
    "affine_flow = bf.networks.CouplingFlow(subnet=\"mlp\", depth=8, subnet_kwargs={\"widths\": (64, 64)})\n",
    "\n",
    "affine_flow_workflow = bf.BasicWorkflow(\n",
    "    simulator=simulator,\n",
    "    adapter=adapter,\n",
    "    inference_network=affine_flow,\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "id": "27b83a8f",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Use a shallower spline coupling flow (default depth is 6)\n",
    "spline_flow = bf.networks.CouplingFlow(subnet=\"mlp\", transform=\"spline\", depth=4)\n",
    "\n",
    "spline_flow_workflow = bf.BasicWorkflow(\n",
    "    simulator=simulator,\n",
    "    adapter=adapter,\n",
    "    inference_network=spline_flow,\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "8aecf471",
   "metadata": {},
   "source": [
    "### Coupling Flow Training\n",
    "\n",
    "First, we train the classic coupling flow:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "f52e8e49",
   "metadata": {},
   "outputs": [],
   "source": [
    "history = affine_flow_workflow.fit_offline(\n",
    "    training_data,\n",
    "    epochs=30,\n",
    "    batch_size=batch_size,\n",
    "    validation_data=validation_data\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "7fc745f3",
   "metadata": {},
   "source": [
    "Next, we train a spline flow for fewer epochs. Spline flows will generally outperform affine flows on multi-modal, low-dimensional problems."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "afa9839f",
   "metadata": {},
   "outputs": [],
   "source": [
    "history = spline_flow_workflow.fit_offline(\n",
    "    training_data,\n",
    "    epochs=30,\n",
    "    batch_size=batch_size,\n",
    "    validation_data=validation_data\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f6ffbb96",
   "metadata": {},
   "source": [
    "## Validation"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c0ca73c8",
   "metadata": {},
   "source": [
    "### Two Moons Posterior\n",
    "\n",
    "The two moons posterior at point $x = (0, 0)$ should resemble two crescent shapes. Below, we plot the corresponding posterior samples and posterior density. \n",
    "\n",
    "These results suggest that these generative networks can approximate the true posterior well. You can achieve an even better fit if you use online training, more epochs, or better optimizer hyperparameters."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "id": "073bcd0b",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABdIAAAF9CAYAAAAA85qPAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjEsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvc2/+5QAAAAlwSFlzAAAPYQAAD2EBqD+naQABAABJREFUeJzs3Xl8VNX5+PHPnSX7vgEhISxBFgNhEQMCsoNVjIBWFhdEW1txR/tDW1trq3UtVdtql68itqKiogatsgYV0KAsgQhBEiAhhCSTyT7ZZjm/P8KMGTJZWLIQnvfrxQu425wzy3PvfebMczSllEIIIYQQQgghhBBCCCGEEB7pOrsBQgghhBBCCCGEEEIIIURXJol0IYQQQgghhBBCCCGEEKIFkkgXQgghhBBCCCGEEEIIIVogiXQhhBBCCCGEEEIIIYQQogWSSBdCCCGEEEIIIYQQQgghWiCJdCGEEEIIIYQQQgghhBCiBZJIF0IIIYQQQgghhBBCCCFaIIl0IYQQQgghhBBCCCGEEKIFkkgXQgghhBBCCCGEEEIIIVogiXRxXvTt2xdN01r88+KLL7q2P3bsGJqm0bdv305r89lo3M/777+/xW2ff/5517YGg6GDWnjmtm7diqZpTJ48+Yz3dfZPCNHxNm7cyJIlS7jkkksICgrC29ubXr16MWPGDP7yl79gMpk6u4lnROJJ+3Ceb51/Pvzwwxa3v+aaa1zbTp8+vUPa+MYbb6BpGrfddtt5Od5tt92Gpmm88cYb5+V4QpyLxMRENE3D29sbs9nc4rZlZWXcfffdxMXF4eXl1eT67MSJE9xyyy1ER0djMBjcPjcX2vv++PHjPPbYY4wdO5bIyEiMRiMhISGMGjWK+++/n2+//bazm3hWWrrHcd5HHDt2rMPb1VbO+4LW/pSVlbn2+f3vf4+mafz+97/vtHYLcT4dPnyYe+65h6FDh+Lv74+Pjw8xMTGMGTOGe+65hw8++OC8PVZzsft8Xxu1N4vFwssvv8xVV11FdHQ03t7eBAQEMGjQIG6++WY+/vhjHA5HZzfzrDQXuy+U825bYvpHH33k2v5c8kMXi66b3RMXpPHjxxMfH+9x3dChQzu4Ne3rrbfe4vnnn8fLy8vj+tdff/28P+axY8fo168fcXFxXfoiXAjRfoqLi1m4cCGbNm0CGi7upkyZgr+/PwUFBezYsYNNmzbxu9/9jk2bNpGUlNTJLe44EiNb9/rrrzN37lyP606cOMH69es7uEVCdF/ffvst+/btA6C+vp7//ve/LQ7EuPPOO3nvvffo27cv8+bNw8fHh8GDBwOglGLevHns3LmToUOHMmXKFIxGIxMmTOiQvpxPzz33HL/97W+pr68nICCApKQkoqKiqKysZP/+/bz88su8/PLL/OpXv+K5557r7OZetBYvXtzsuubuf4S40K1du5ZFixZRV1dHeHg448ePJzIyktLSUvbu3cvf//533nnnHa6//vrObmqXsWHDBm6++WZMJhMGg4HRo0czceJEbDYb2dnZvPXWW7z11luMGTOGnTt3dnZzL1qzZs2iZ8+eHtf16dOng1tzYZNEujivfvazn10w35qei8suu4zvvvuOjz/+mJ/+9KdN1u/YsYPMzEzGjBlzwY6oaYuDBw92dhOEuKiUl5czYcIEDh06xODBg/nXv/7FxIkT3bapq6tj1apVPP7445w8ebKTWnrmJJ60L71ez7Bhw/j8888pKCjweCG9atUq7HZ7tz93CdFRXnvtNQB69+7NiRMneO2115pNpFutVj788EN8fHxIT08nKCjIbX1OTg47d+6kT58+pKenN/m149NPP80jjzxCr1692qcz58kjjzzCs88+i9Fo5IUXXuCee+7B29vbbZtvvvmG3/zmN/zwww+d1Mr2sXnzZqxWK7179+7sprRJVx9lKcT5VlhYyOLFi6mrq+Ohhx7iySefxMfHx22bXbt28f7777d7W+bOncvYsWMJDg5u98c6F59++inXXXcddrud22+/naeffpqoqCi3bXJzc/nTn/7EmjVrOqmV7eNCOe86PfLIIzLK/DyR0i5CnIXbb78daH7UufPGyblddzV48GDXSCkhRPu79957OXToEH379mX79u1NkugA3t7e3Hnnnezdu5chQ4Z0QivPjsST9nf77bdjs9lYtWqVx/UrV67Ex8eHRYsWdXDLhOh+qqurefvttwH4z3/+Q0BAAPv372/2S6qTJ09is9no0aNHkyQ6NCQiAPr16+exZGCvXr0YPHhwl066bN68mWeffRaAd999l4ceeqhJEh1g7NixbNq0iYceeqijm9iuBgwYwODBgzEajZ3dFCGEB5988glVVVVER0fzwgsvNEmiA4wePZqnn3663dsSHBzM4MGDu3SS1mw2c/PNN2O327nvvvt47bXXmiTRoWG08z/+8Q+38iHdwYVw3hXtQxLpokvKy8vj3nvvZeDAgfj4+BAcHMz48eP55z//id1ud9s2JSUFTdNITk5ucpylS5eiaRpGo5GKigq3dV9++SWapnHllVeecfuGDRvGZZddxoYNGzhx4oTbuqqqKtasWUNMTAwzZ85s9hgHDhzg8ccfZ/z48fTu3RsvLy/Cw8OZPn26x29rb7vtNvr16wc0jEo6va7V6Xbt2sXixYvp168fPj4+hIWFkZiYyK9+9StycnI8tslqtfLss89y6aWX4uvrS3h4OPPmzWt2pGhzj924jlhqaiozZ84kNDQUX19fRo0axZtvvtns82I2m7nvvvvo06cP3t7exMXF8cADD1BWVnbB1CEToj0cOXKE1atXA7BixQrCwsJa3L5Hjx4MGjSoyfJ33nmHadOmERYW5vqM3X777c2O/Dvbz3N5eTmPPfYYw4YNw9/fH29vb6Kjoxk/fjy/+93vsFqtbts3F09OnjzJ/fffzyWXXIKPjw9+fn7ExsYybdo0XnjhBdd2ZxMjb7rpJlesCQsLY9asWfzvf/87r88DNJRkWLt2LbNnz6Znz554eXnRs2dPJkyYwLPPPktNTQ3Q8DN6TdNavEFbs2YNmqZx+eWXN7tNc2666Sa8vb1ZuXJlk3VffPEFWVlZzJ07l5CQkBaPk5mZyZIlS4iLi3M9d9OmTWtxpJHNZuPFF19k2LBh+Pj4EBkZyfXXX8/+/ftbbfcPP/zAL37xCwYMGOC6Jrjyyiv573//2+q+QnSW9957j4qKChISEpgyZQrz588Hfhxs0ZimacTFxQFN45ezTu6kSZOAhs9q4/XOMlbNXSM1rl9tMpm4++67iY2NxcvLi9jYWO699163etenO5+fvyeffBKA5OTkZktMNX5OPH1ZfKbxp7Vrx+bqEDdebjabufvuu92uTR988EFKS0tb73QjzdXZnTx5MpqmsXXrVvbu3cu8efOIiIjA29uboUOH8uc//xmllMdjWiwWfvvb3zJw4EDXefb222/nxIkTXaZ2+fr165k9ezZRUVF4eXkRHR3N/Pnz+e6775psO2/ePDRNY+3atW7LbTYbwcHBaJrGjTfe2GS/22+/HU3T2qW0prh4FBYWAhAZGXnG+zb+fH/44YdMmDCBoKAgAgMDmTx5crPXl81pLjY1rl99pvfuAKWlpTz++OOMGDGCwMBA/Pz8GDZsGE8++STV1dVn1Ma//e1vlJWVERUV1aYyXJ7yLiUlJfz617/m0ksvxc/Pj8DAQEaPHs1zzz3nuj5urLXa8c3NU9F4uc1m47nnnnM9bxEREdx4441kZma2qd9O7XHeVUrx+uuvc9lll+Hn50d4eDg/+clP2LFjR5epXX4m5+GXX34ZTdO47777mqy7+uqr0TSNnj17NjnHvfnmm2iaxq233tpu/TgnSojzIC4uTgFq5cqVbdr+6NGjClBxcXFN1u3cuVOFhYUpQPXp00fNnz9fXXXVVcrHx0cBatasWaqurs61fXl5uTIYDCooKEhZrVa3Yw0cOFABClAff/yx27rf/va3ClBPPPHEGffzq6++Uq+88ooC1JNPPum2zWuvvaYA9Zvf/MbVT71e3+RYd9xxhwLU4MGD1axZs9T8+fPVuHHjlE6nU4B68MEH3bb/97//ra6//noFKH9/f7V48WK3P40999xzruNccskl6sYbb1TXXnutGjJkSJPXKTU1VQHqiiuuUNOnT1d+fn7qqquuUtdff72KjY1VgAoJCVFHjx5t0gfnc9vc8/Tb3/5WaZqmRo8erRYsWKDGjh3r2ucvf/lLk/3y8/PVgAEDFKDCwsLUvHnz1Jw5c1RoaKgaNGiQmjNnzhm9z4ToTl566SXX59Fms53x/g6HQ916660KUAaDQU2dOlUtWLBAXXLJJQpQfn5+6rPPPmuy39l8ni0Wi0pISFCAioyMVNdee61asGCBmjx5surZs6cCVGlpqds+nuLJyZMnVXR0tOt8cN1116n58+eriRMnqrCwMBUcHOza9kxi5IsvvuiKkSNGjFA33HCDmjBhgvLy8mr2vHC2ca2+vl7NmzdPAUqn06mxY8eqhQsXqhkzZqjevXsrwBVfd+3a5eprc6/xlVdeqQC1atUqj+tPd/p56MYbb1SA2rZtm9t2t9xyiwLUxo0b1cqVKxWgpk2b1uR4n3zyiet8PGjQILVgwQI1depUpdfrFaBuv/32JvvY7XZX/Pby8lIzZ85U8+fPV3379lU+Pj5q6dKlCmjyOiml1Jo1a1yPN3jwYDV37lw1depU5e/vrwC1ZMmSJvssXrxYzhWi002cOFEBasWKFUoppbZv364AFRwcrKqrq922Xbx4cbPx66uvvlKLFy9Ws2bNUoDq0aOH23qTyeQ6hqf3/eOPP+76bMbExKgePXqoefPmqauvvloFBwcrQI0ZM0bV19c36cPZfP6aU1pa6oq7H3zwwRk+mw3OJv60Fg+c8e70+ONcnpycrAYMGKBCQkLUnDlz1Ny5c1VoaKirDUVFRW77tXSP4zyPnH5NPWnSJAWoRx55RHl5eakhQ4aoBQsWqEmTJrn6dv/99zc5XlVVlRozZowCVEBAgJo9e7b66U9/qnr16qWioqLUbbfdpgD1+OOPt/CsunPeF5xJmsD5HvP0OI899pgClKZpavz48WrhwoVqxIgRrvPSa6+95rb93//+dwWou+66y2258/MDqPDwcOVwONzW9+nTRwEqJyenze0W4nT/+c9/XO/NTZs2ndG+zs/3gw8+qAB12WWXqYULF6rLL7/c9d59+eWXm+zXXIxqLjady737999/79qmV69e6qqrrlLXXnut6tGjh+uauKysrM19HjlypALUvffe2+Z9GsvOznY9b5GRker6669XycnJKjAwUAFq1KhRqqSkxG2f5p4Xp+ZicOPl8+bNU0ajUU2fPl0tWLBA9e/f3xVHd+zY0eSYzcXu9jjv3nXXXa57hkmTJqkFCxaoSy+9VOn1evXQQw8pQE2aNKm1p9aN8/2Xmprapu2d7zFPj3Om5+Hvv/9eAWrIkCFuy+vr613XEoBKT093W++8L2nr/U5Hk0S6OC/OVyK9trbWdaxf/vKXbsElOztb9e3bVwHq17/+tdt+48aNU4Davn27a1lOTo4C1PDhwz0GeE/7tLWfX331lSorK1O+vr4qPj7ebZvx48crTdNUdnZ2i4n0rVu3quzs7CbLMzMzVUxMjAJUWlqa27qWLs6dPv74YwUoHx8f9e677zZZ//3336sDBw64/t/4gnnkyJHq5MmTrnU1NTWuG7c777yzybFaS6QbjUa1bt06t3XOk5+nG8m5c+cqQE2ePFmVl5e7lpeWlqoJEya4Hk+SI+Ji5LygmDp16lnt/+qrrypARUREqD179riWOxwO1wVfSEhIk6TA2XyeV61apQD1k5/8pMlFot1uV1u3bnX7QlQpz/HkiSeecMWf02+a6+vrm9zktCVGfv7550rTNBUREaG++OILt3X79u1zxd+tW7ee8/OglFLLli1TgOrbt6/au3ev2zqHw6E2bdrkdtMyfvx4Bai1a9c2afv+/ftdNxu1tbXN9rGx089D69evb3KhW15ervz8/FTfvn2Vw+FoNpFeUFDgugF48skn3V6Tb7/91pVc+te//uW239/+9jdXArDx+cdqtbpuGDzdFO3bt095e3srHx+fJom3Y8eOqWHDhnm8yJZEuuhshw4dcsWLxjF18ODBClBvvvlmk31ai18t3dgq1foNPaBuu+02t9iRm5vr+kJv9erVbvud7eevOZs3b3a1Izc3t037NHa28edcE+mAGjt2rDKbza51paWl6oorrlCAWrBggdt+55JIB9Q//vEPt3WbN29WmqYpvV6vjh8/7rbOmbAbOnSoys/Pdy2vqalRN9xwg+uYnZVI/+yzz1z3JBs2bHBb93//93+uz0hGRoZrufOzM3DgQLftndcDzvu6Xbt2tbqPEGeqsrLSFRM1TVOTJ09Wf/zjH9Wnn37a5Pr4dM7Pt6Zp6r///a/bunfeeUdpmqYMBoPav3+/27qzTaSf6b17dXW1a8DaY4895nYdbrFY1MKFCxW0/QtSq9Xq+nLU0zmtLZKSkhQ0fGFZVVXlWl5UVKRGjRqlALVo0SK3fc41ke68F2qcuLXZbOree+917Xf6NfbZJtLP9LzrzOMEBAQ0yVH9+c9/dh2zsxLpZ3sedg6KOnHihGvZF1984RbT//znP7e6T1ciiXRxXjiDS3N/Tv8QNhfgnN8CR0dHe0wSvP/++wpQgYGBqqamxrXcObr897//vWuZc2T466+/rqKiotTgwYNd61oaxd6Wfn711VdKKaVuuukmt4RLZmamKxHcuJ+eEukt+ec//6kA9atf/cpteVuSRM5RHqcHo+Y4A6WmaU2SPEop9c033yhA9e/fv8m61hLpy5Yt8/iYzhvJL7/80rXs2LFjStM0pdPp1MGDB5vss3//fqVpmiRHxEXrqquu8njT3lbOi2dPo2EcDofrQuapp55yW3c2n+fnnntOwY8jMdvCUzxxjlT2lFT2pC0x0nnR/v7773tcv2bNGgWo66+/3m352TwPhYWFrlHu3333XZv64Hx8T6PBf/GLXyhAPfroo206llJNz0N2u1316dNHBQQEuG5a/vGPf7idQ5tLpP/xj39UgBo9erTHx3rhhRc8JjPi4+MVoF599dUm+9TU1Lh+pXD6TdH8+fMVoF544QWPj7dz506P7ZFEuuhsy5cv9xhHnLHR001weyfSY2JilMViabLfM888o6DpKLKz/fw155133nHF+bZ+EdjY2caf85FIb/zls9O+fftc162NE9znkkifN2+exzY6z/+Nk1XV1dUqICBAAWr9+vVN9ikqKlJ+fn7nlEhv7k9z77HTH2fatGktnjdnz56tAPXzn//cbblzxGzj0eUTJ050fakDqGeeeca1rrlR7EKcjczMTNe14ul/RowYoV599VWPvxp0fr7nzJnj8bjOXx2d/n4/20T6md67OwfUzJ4922P7KisrVVRUlDIYDE1GgXtSUFDgel4+//zzVrc/3VdffaWg4RexBQUFTdZ/9913CmgSY89HIv3FF19ssl9tba0rwf3WW2+5rTvbRPqZnnenTp3a4nW+8xdIZ5tIb+5Pc++x0x/nbM/DzsFgb7zxhmuZM3/38ccfK4PBoK666irXuuZGsXclUiNdnFfjx49n8eLFTf5cddVVbdp/69atACxYsMDj5EPz5s0jNDSUyspKdu3a5Vo+ffp0ADZt2uRa5vz3zJkzmTZtGpmZma565lu3bsVmszFp0iSPEza11emTjjr/busko1VVVbz33nv8+te/5s477+S2227jtttu44MPPgDg0KFDZ9SegoIC9u7di06n44477jijffv06UNiYmKT5c7JCk+vBd8W1157rcflno751VdfoZRi1KhRHiccTEhIYPjw4WfcBiFEw7wT2dnZQEMd7tNpmsaSJUsASE1N9XiMM/k8jxkzBoDnnnuON998k5KSkrNqt7MO+COPPMLatWupqqo6q+M4FRcXs3PnTnx9fZvtj7Pu4I4dOzyuP5PnITU1lfr6ekaPHs3o0aPb1Ma5c+cSGxvL5s2b3Wo1lpeX89///he9Xs9dd93VpmN5otPpWLx4sWs+D2g4d+l0umbrTTo5z9Ge3kOA67xz+PBh8vPzgYbnIysrC4Cbb765yT4+Pj4e6906HA4+++wzAFdt6dNddtllBAQEsGfPHmpra1tsuxAdpfGEvqdfD956660YDAa+/PJLV0zuKNOmTcPPz6/Jck+xqyt+/s4m/pwPiYmJjBgxosnyYcOGMXLkSBwOB19++eV5eawzOb/s2rWLqqoqIiIiPM7JFBkZyYwZM86pPZ7u6RYvXkx8fHyr+9psNrZv3w7Q7LnF+Zqdft3hvK/buHEj0FAH/ptvvmHChAnMmjULo9Ho8Z7PuZ8Q52LQoEF88803pKWl8bvf/Y5Zs2a5aqbv3buXu+66i6uuuor6+nqP+zcXo5zLnbHsXJ3pvfunn34KNB/TAwICuOyyy7DZbM1Oin0+OZ+Hq666ih49ejRZP3r0aBITE3E4HHzxxRfn9bE9vUbe3t6u5+Z8vUZnct612Wyue4+bbrrJ4/EWLVp0Tu2ZNWuWx5g+YcKENu1/tufh5nJ1fn5+XHXVVYwZM4avvvrK9Zm6EGL62WcQhfDgZz/7Was34i1xBhPnhHGn0zSNfv36UVpa6hZ4xo0bh7+/P2lpaVRVVeHv78+WLVsYMmQIvXv3Zvr06bz99tts2rSJxYsXn7cP55QpU+jXrx/vv/8+L774Im+++SZBQUHccMMNre67bt06lixZgtlsbnab0ydIbU1ubi7QMIP0mc4e3adPH4/Lg4KCAKirqzuj47XlmI1vvPLy8gCaTAzSWN++fUlPTz/jdgjRHTgv4ouKis54X2e8DA8Pd33+TjdgwAC3bU93Jp/nyZMns3z5cp5//nnXBJoDBw5k/PjxXHfddVx77bXodK1/l3/LLbewceNG3nrrLa6//nr0ej1Dhw5lwoQJ3HDDDUydOrXVYzR29OhRlFLU1NR4/LK2MZPJ5HH5mTwPzomdPX052ByDwcDSpUt59NFH+dvf/sbf/vY3AFatWoXFYnEl2s/FkiVLePLJJ3n99de5/PLL2blzJ9OnT3dNdtic1s7RISEhhIWFUVJSQl5eHtHR0a7YHhERQUBAgMf9PB3PbDa7zoFt6a/ZbKZ3796tbidEe/v0008pKCigd+/ezJo1y21djx49uPrqq0lJSeH111/nqaee6rB2nUnsao/PX+PJ+4qKis44jp1N/Dkfmns857rdu3e74ty5ao/r5nPR3AStbWE2m13tbe45bO66Y/r06axcuZJNmzZxxx138MUXX2C1WpkxYwb+/v6MHTuWbdu2UVtbi5eXF6mpqeh0ujO+JhCiJZdffrlrQIdSij179vD888/zzjvvsGnTJl566SV+9atfNdmvufe7c3lHxYvT792PHDkCNFxb33LLLS0eu7lr4MbCw8PR6XQ4HI5zujdpKcYOGDCA9PT0sxrQ15yQkBBCQkI8ruvo16hxTC8uLnb9v7nYfa4x/ZFHHjmniUrP9jzszLlt3rwZaMhxffvtt8yYMQMvLy+mT5/O119/zddff82kSZMkkS5ERzEajVx55ZV89tlnbN26lT59+lBYWOj6VrHxyIbzmUh3zhj9+OOPs3jxYgoKCrjzzjvx9fVtcb8TJ04wf/58ampq+H//7/9x00030bdvXwICAtDpdGzYsIFZs2Y1mb24PbUlqdURx9Q07azWCdHdjR49mv/85z/s3r0bu92OXq/v0Mc/08/zM888wy9/+UvWrVvHtm3b2L59OytXrmTlypWMGTOG1NRU/P39W33M//73v/z617/m008/Zfv27Wzfvp1XX32VV199lWuvvZYPP/ywzc+Fw+EAGkbdXH/99WfUn8Ztam8///nP+cMf/sCbb77J008/TUBAAK+88goA99xzzzkfv1+/fkyePJnU1FQeffRRoO2/pOooztcKmh/50lhrX4wI0VFee+01oOEGedKkSU3WO29E33jjDf7whz90WCw/k9jVHp+/kSNHupIu33777Tl/IXi+NO7r2Tpf1+ty3dxg2rRpaJrG5s2bUUq57tucI+ynT5/OV199xbZt2wgKCqKsrIwxY8Y0mxwT4lxpmsaoUaN4++23qa6uJiUlhY8++shjIr01nRUvnLGuuRHgjbU2sAIaBn4MHz6cvXv38u2337aanO8oF3pMb8mFGtOjo6MZMmQIBw8eJCMjgyNHjmCz2dxi+h//+Ec2btzI+PHj+eKLLzAYDOeU9G9vkkgXXYpzNIvzG1NPjh496rat0/Tp0/nss8/YtGmT69s/Z6K8T58+DBw4kM2bN3PixAkOHjxIdHQ0Q4cOPec233bbbTzxxBOsW7cOaFsyYt26ddTU1DB37lyeffbZJusPHz58Vm1x9vvkyZOUl5ef8aj0zuR8PY8dO9bsNi2tE6K7mz17NsuWLaOsrIyUlBTmzp3b5n2dny/nKENPo9Kdcfd8jurt27cv9957L/feey8A3377LTfffDPffvstzz33HE888USbjjN06FCGDh3Kr371K5RSbNmyhUWLFrFu3TrefPNNV1ma1jgTN5qmucqZtCdnTG5coqUtwsPDuemmm/i///s/3nzzTS655BIOHTrE0KFDz9uIu9tvv53U1FTWrVtHaGhom95PvXv3JjMzs9lzdHl5uauMj/N95Py7uLiYqqoqj6PSPcX2iIgIfH19qamp4YUXXiAiIqKtXROi05w8eZL//e9/QEO8dZa28CQ/P5/PP/+ca665pqOa12bt8fkLDQ1l4sSJfPHFF6xatYp58+ad0f5nE38AvLy8AKisrPS4n/OXQ81x3nd44oxdMTExLR6jPXT16+bw8HC8vb2pq6vjyJEjHsszNnfd0aNHDxISEti/fz/p6els2rSJiIgIV4md6dOn8/jjj7Np0ybX9UxXHrkoupeZM2eSkpJCcXGxx/VHjx71WHKlM+MFNFwDZ2Zmcscdd7Tp1/Ntcd1117F3717effddnn/++TMa1NCWvI+nGHGuMb2srIyysjKPX7x15mvUOGbm5OR4zFN1di7kbM/D0BCjDx48yKZNm1z7O+O2s7rEpk2buPrqq6moqGDcuHHN/oq6K5Aa6aJLcX7r9O6773qst/jhhx9SWlpKYGBgk3qzjUedb9q0qcm3WNOnT6egoIAXX3wRaBjtcD706dOH6667jvDwcMaOHUtSUlKr+zgDjKdve5VSrF692uN+zhOHzWbzuL5nz56uWmLOeu0XiokTJ6JpGrt27eKHH35osv7AgQNS1kVc1AYMGMDChQsBeOihh1qtO15UVOSaZyEmJsb1E2pPP9VWSrmWT5ky5fw1+jRjxoxh6dKlQEOdybOhaRrTpk1z1QlsfJzWYmR0dDTDhw+nsrKSzz///Kwe/0xMnToVLy8vdu3axe7du89o3/vuuw+Av//9767yLnffffd5a9v1119PXFwc4eHhLFmyBB8fn1b3cZ5TnfWfT+c87wwcONB1AR0TE0P//v0BPJ7b6urqeO+995os1+v1rpEqzlruQnR1b7zxBna7naSkJJRSzf75f//v/wE/jl7vatrr8/eb3/wGgJSUFD788MMWt1VKsW3bNtf/zyb+wI838wcPHvT4GM5a8M3Zt28f+/bta7L8+++/Z/fu3eh0Oq688soWj9EeRo8ejZ+fHyaTya3urFNxcbGrxnhnMBgMrrq7zZWIcb5mnq47nPd1b731FhkZGa5R6tBQciMoKMh1z9d4eyHORVtGIjtLqTaXbP3Pf/7jcfmbb74J0GmjbH/yk58A5zem33vvvQQHB1NUVMTy5ctb3f6rr75y/dv5PHz++ecUFhY22XbPnj2uud8ax1hnTG9ukIqzFnxLPL1G9fX1vPvuu25t60hGo5Fx48YBnq+XAd5+++2ObFITZ3sehqa5up49ezJs2DDgx+oS3333He+//77b9l2VJNJFl/LTn/6UPn36kJ+fz7Jly9ySIUePHuWhhx4CGoL26Tf9w4YNIyoqigMHDpCamsrYsWMJDAx0rXd+GJ0JifP54Vy7di3FxcV8/fXXbdreOcHE+++/z8mTJ13L7XY7v/vd75qd5C4yMhIvLy8KCgqaTaI9/vjjQMPNinPS0sYOHDjg8Wais/Xt25drr70Wh8PBXXfd5fYtc3l5OXfddVeHlroRoiv661//Snx8PEePHmXChAluSQan+vp6Xn/9dUaOHOn2WX/44YcB+OMf/+j2pZRSiieffJK9e/cSEhLCz3/+83Nu54cffsiXX37Z5OeVVqvVlcBuy89G33zzTbeJpZ0qKytdE940Pk5bYuSTTz4JNNQJd/6SqDGlFGlpaWzYsKHV9rUmKirKNTHoT3/6UzIyMpo81pYtWygvL2+y77Bhw5g6dSoHDx4kJSWFoKAgbr311nNuk5Ovry/Hjh2juLiYP//5z23a5+c//zlBQUHs3r2bP/3pT24xec+ePa7n9vSfOj/wwAMA/P73v3e78bHb7Tz88MPNTgz4+OOP4+Xlxa9+9StWrVrl8ee6GRkZrF27tk3tF6K9OW8iWyuH4vwsf/LJJ22qRdsZ2uPzN2PGDNe1/IIFC1ixYoXHOXh27drFrFmzeOGFF1zLzjb+OK/3//Of/3DgwAHXcqvVyvLly1udVE8pxV133UVpaalrWePr0uuvv75TytT4+fnxs5/9DIAHH3zQLRFVV1fHPffcg8Vi6fB2NeZ8rV999VVXbVynN954g5SUFIxGI/fff3+TfRvftyml3CZONRgMTJo0ib1797J9+3Z8fX0ZP358O/ZEXCxeeeUVFi9e7PFeXCnF2rVrXbmEBQsWeDzGhx9+yDvvvOO27P333+eDDz7AYDC4fqXZ0e68807i4uJ47733WL58uccR3QUFBfz73/9u8zHDw8N588030el0vPTSS/zsZz/zWC/9xIkT3HPPPcyZM8e1bMKECSQlJVFTU8MvfvELqqurXeuKi4v5xS9+ATQ8z41jrPOLtAMHDjRJiL/33nu8/PLLrbb7j3/8o9s1ucPhYPny5eTl5REbG3vW5R/PlXMQzcsvv8w333zjtu6ll14iLS2tM5rlcrbnYWhIwhsMBrZs2cLBgweb5OKmT5+O3W7n1Vdfdf2/S1NCnAdxcXEKUCtXrmzT9kePHlWAiouLa7Ju586dKiwszLV+/vz56uqrr1Y+Pj4KULNmzVJ1dXUej7tw4UIFKEA98cQTbutKS0uVTqdzrT9x4sSZdtPVz6+++uqM+qnX692WW61WNXr0aAWogIAAdc0116gbb7xRxcXFKaPRqJYvX64ANWnSpCbHvOGGGxSgYmNj1cKFC9Udd9yh7rjjDrdtnnrqKaVpmgLU4MGD1fz581VycrIaOnRok9cpNTW12cdycj5nbV3ufJ6OHj3q8XiLFy/2+H45ceKE6tu3rwJUeHi4mjdvnpo7d64KCwtTAwcOVMnJyQpQb731VrNtFaK7KywsVJMnT3Z9/vr166euu+46tXDhQjV16lQVEBCgABUUFKTS0tJc+zkcDnXLLbcoQBkMBjVt2jS1cOFCNWjQIAUoX19f9b///a/J453N5/n+++9XgIqIiFAzZsxQN910k0pOTlZRUVEKUL1791bHjx93O46neHLdddcpQEVHR6urr75a3XTTTerqq69WwcHBClAJCQmqoqLCbZ+2xMiXXnpJGQwGBaj4+Hh1zTXXqEWLFqkZM2a42rh8+fJzfh6UUqqurs4Vu3Q6nRo3bpxatGiRmjlzpurdu3eLx/zoo49cz8u9997rcZvWNHceasnKlSsVoKZNm9Zk3bp161zn48GDB6uFCxeqadOmuZ7PJUuWNNnHbrera6+9VgHKy8tLzZo1Sy1YsED169dP+fj4qLvuuksBavHixU32XbNmjfLz81OAiomJUTNnzlQ33XST+slPfqJiYmIUoObPn++2T3OvhRDtaevWrQpQ3t7eqqSkpNXtR40apQD1wgsvKKVavjZWqvXrtebe948//rgC1OOPP37Gxz2bz19bPPXUU8poNCpABQYGqunTp6tFixap5ORk13Wgpzh8NvFHqR/PJb6+vmrGjBkqOTlZxcTEqKCgINf56vT444yDycnJqn///iokJETNnTtXzZs3z3WfMnDgQFVYWOi2X0uvY3PnkUmTJilApaamemx/c69hZWWl2/1EcnKyuvHGG1V0dLSKiIhwvSeeeuopj8f1xPl+OJM0QUvvsccee0wBStM0NWHCBLVo0SLXe1+v16vXXnvN4zGrqqpc7xFA5eTkuK1/6aWXXOtmzJjR5rYK0ZK//OUvrvdVZGSkmjlzplq0aJG6+uqr3WLTzTffrOx2u9u+zs/3Aw88oAA1ZswYtWjRIpWUlOTab8WKFU0es7nY7YxBp8emc7l3z8jIcPUjJCREXXnllWrRokVqzpw5aujQoUrTNNWjR482P19O//vf/1RERITrHmPs2LFq/vz56vrrr1cjRoxw5SXGjh3rtl92drbreYuKilI33HCDuu6661RQUJAC1KhRozyeTxu/TuPGjVM33HCDuvTSS5Wmaeq3v/2txxjsjM19+vRRc+fOVUajUc2YMUMtWLBADRgwQAHK39/fY56nudjdHufdO++80xUfJ0+erBYuXKgSEhKUXq9XDz744FnFPOdz1dw55kzad7bnYaWUGjdunKstq1atcluXnp7uWufv76/q6+vPpIsdThLp4rw4n4l0pZTKzc1Vd999t+rfv7/y8vJSgYGBaty4cerVV19VVqu12eO+9tprrg/g9u3bm6wfM2aMAtSQIUPa1M7Tna9EulINF7+//vWv1aBBg5SPj4+KiopSc+bMUd99912LwctsNqtf/OIXqk+fPm4XmKf7+uuv1cKFC1Xv3r2V0WhUYWFhKjExUf2///f/3C5Gu1IiXSmlioqK1N13361iYmKUl5eXio2NVXfffbcym81q6tSpClDr169vtq1CXCw+++wzdeutt6r4+HgVEBCgjEaj6tmzp5oxY4Z68cUXldls9rjf6tWr1eTJk1VISIgyGo0qNjZW3XbbbSozM9Pj9mfzed6zZ4965JFH1IQJE1Tv3r2Vl5eXioyMVKNHj1Z/+tOfVHFxcZPjeIonX375pXrggQfU5Zdfrnr27Km8vLxUz5491bhx49Rf//pXVVVV1eQ4bY2R+/fvV3feeacaOHCg8vHxUX5+fqp///5q1qxZ6uWXX27yZeu5xDWHw6FWr16tZs6cqcLDw12v1cSJE9Xzzz+vampqPB6zsrJS6fV6pWlas69Pa853Il0ppQ4cOKAWL16sYmJilNFoVCEhIWrKlCnqnXfeafaYVqtV/fnPf1ZDhw5V3t7eKjw8XF133XVq7969zd4sNu7Dgw8+qBISEpS/v7/y8fFRcXFxavLkyeqZZ55RWVlZbttLIl10BucXlTfccEObtn/xxRfdrkm7YiLd2a4z+fy11bFjx9Sjjz6qxowZo8LDw5XBYFDBwcFq5MiR6v7771e7d+/2uN/ZxJ/a2lr12GOPqf79+yuj0aiioqLUwoULVVZWVrPxp/HyoqIi9Ytf/MLt2vS+++7zeJ7tyES6Uj/eTzjvmXr27KluueUWlZOTo26//XYFqH/+85/NPjenO9+JdKUarleuvvpq1+vcs2dP9dOf/tTty35PJk6c6PrC4nTff/+9q53PPvtsm9sqREsqKirURx99pO699151+eWXu+KMr6+vGjBggFq4cKH67LPPPO7b+PO9Zs0aNW7cOBUQEKD8/f3VxIkT1bp16zzu15GJdGcfn3vuOTVu3DjXvUCvXr3UmDFj1K9+9Su1Y8eOZo/bksrKSvWXv/xFzZgxw3XN7ufnpy655BJ18803q08++UQ5HI4m+5nNZvXoo4+qIUOGuK7HR44cqZ555hlVXV3d7OOtWrVKjRo1Svn4+KigoCA1depUtXHjxmZjcOPlVqtVPfXUU2rw4MHK29tbhYWFqeuvv159//33Hh+rIxPpDodD/fvf/3b1LSQkRM2cOVN9+eWX6s0331SAWrhwYbPPiyfnM5Gu1Nmdh5VSri85oOmgVofD4RrM9JOf/KRN7exMmlJSK0EI0fWVlZXRv39/ysvLKSwslInnhBDd3v/93//x85//nJkzZ7J+/frObo4QQlw03njjDZYsWcLixYubrfHdlVmtVhISEvjhhx/YtWsXo0aN6uwmCdGt9e3bl5ycHI4ePUrfvn07uzniNMeOHaNfv37ExcV1+qSdZ+v2229n5cqV/PnPf2bZsmWd3ZyLmtRIF0J0KTt37myyzGQysXjxYkpLS5k9e7Yk0YUQ3Z7FYuHpp58GfqwzK4QQQjS2a9euJvXrq6qquOeee/jhhx8YPny4JNGFEOIC8f333zeZ38LhcPDvf/+bN954Ax8fHxYuXNhJrRNOhs5ugBBCNJaUlERMTAxDhgwhPDycEydOsGfPHqqqqujTp49rghchhOiOnn/+eTIyMti2bRtHjhzhqquuYubMmZ3dLCGEEF3Q9ddfT3V1NcOGDSMqKoqioiL27t1LSUkJYWFhF+RoeiGEuFg9//zzrFmzhpEjR9K7d28sFgsHDhzg2LFj6PV6XnnlFXr16tXZzbzoSSJdCNGlPPbYY2zevJn09HRKS0vx8vJiwIABzJ49m2XLlhEeHt7ZTRRCiHbz6aef8sUXXxAREcFtt93GihUrOrtJQgghuqhly5bx4YcfcuDAAbZv345OpyMuLo6bb76Zhx9+mNjY2M5uohBCiDaaP38+FRUV7Nq1i71792Kz2YiKimL+/Pk88MADjB07trObKACpkS6EEEIIIYQQQgghhBBCtEBqpAshhBBCCCGEEEIIIYQQLZBEuhBCCCGEEEIIIYQQQgjRAkmknwdKKSoqKpAqOUIIcWGTeC6EEBc+ieVCCHHhk1guhOiKJJF+HlRWVhIcHExlZWWHPabD4aCgoACHw9Fhj9nRunsfu3v/QPooLjwdHc8vhveP9PHC1937BxdHHy8mEsvPP+lj99Dd+9jd+3exkTzL+dfd+wfSx+6iK/dREulCCCGEEEIIIYQQQgghRAskkS6EEEIIIYQQQgghhBBCtEAS6UIIIYQQQgghhBBCCCFECySRLoQQQgghhBBCCCGEEEK0QBLpQgghhBBCCCGEEEIIIUQLJJEuhBBCCCGEEEIIIYQQQrRAEulCCCGEEEIIIYQQQgghRAskkS6EEEIIIYQQQgghhBBCtEAS6UIIIYQQQgghhBBCCCFECySRLoQQQgghhBBCCCGEEEK0QBLpQgghhBBCCCGEEEIIIUQLJJEuhBBCCCGEEEIIIYQQQrRAEulCCCGEEEIIIYQQQgghRAskkS6EEEIIIYQQQgghhBBCtEAS6UIIIYQQQgghhBBCCCFECySRLoQQQgghhBBCCCGEEEK0QBLpQgghhBBCCCGEEEIIIUQLJJEuhBBCCCGEEEIIIYQQQrRAEulCCCGEEEIIIYQQQgghRAskkS6EEEIIIYQQQgghhBBCtEAS6UIIIYQQQgghhBBCCCFECySRLoQQQgghhBBCCCGEEEK0QBLpQgghhBBCCCGEEEIIIUQLum0i/e9//zt9+/bFx8eHpKQkdu7c2ab93nnnHTRNY86cOe3bQCGEEK2SWC6EEBc+ieVCCHHhk1guhBDdNJH+7rvvsmzZMh5//HF2795NYmIis2bNoqioqMX9jh07xsMPP8zEiRM7qKVCCCGaI7FcCCEufBLLhRDiwiexXAghGnTLRPqKFSv4+c9/zpIlSxg6dCj/+Mc/8PPz4/XXX292H7vdzk033cQTTzxB//79O7C1QgghPJFYLoQQFz6J5UIIceGTWC6EEA26XSK9vr6eXbt2MX36dNcynU7H9OnT+frrr5vd7w9/+ANRUVHccccdHdFMIYQQLZBYLoQQFz6J5UIIceGTWC6EED8ydHYDzrfi4mLsdjs9evRwW96jRw8yMzM97rNt2zZee+019u7d26bHqKuro66uzvX/iooKABwOBw6H4+wafoYcDgdKqQ57vM7Q3fvY3fsH0sf2otN1u+9Am+iIWA6dH8/lM9I9dPc+dvf+Qef1sbvHc4nl3Yf0sXvo7n2UWN4+LpZY7nws+Yxc2KSP3UNXzrN0u0T6maqsrOSWW27h3//+NxEREW3a5+mnn+aJJ55ostxkMlFbW3u+m+iRw+GgvLwcpVS3PXF39z529/6B9LG99OzZs0Me50JyNrEcOj+ey2eke+jufezu/YPO66PEc3cSy7su6WP30N37KLG8a7hQYznIZ6Q7kD52D105z9LtEukRERHo9XoKCwvdlhcWFnp8UrKzszl27BjXXnuta5nzGw+DwcChQ4cYMGCA2z6PPvooy5Ytc/2/oqKC2NhYIiMjCQoKOp/daZbD4UDTNCIjI7v1B6c797G79w+kj+LsdUQsh86P5xfD+0f6eOHr7v2Di6OPnUFiefchfeweunsfu3v/OsvFEsud7ezO76Hu3j+QPnYXXbmP3S6R7uXlxejRo9m8eTNz5swBGl6AzZs3c8899zTZfvDgwezfv99t2WOPPUZlZSUvvfQSsbGxTfbx9vbG29u7yXKdTtehL7CmaR3+mB2tu/exu/cPpI/i7HRELIeuEc8vhveP9PHC1937BxdHHzuaxPLuRfrYPXT3Pnb3/nWGiymWQ/d/D3X3/oH0sbvoqn3sdol0gGXLlrF48WIuu+wyLr/8cl588UUsFgtLliwB4NZbb6V37948/fTT+Pj4kJCQ4LZ/SEgIQJPlQgghOo7EciGEuPBJLBdCiAufxHIhhGjQLRPp8+fPx2Qy8bvf/Y6CggJGjBjB559/7pocIzc3t8t9oyGEEMKdxHIhhLjwSSwXQogLn8RyIYRooCmlVGc34kJXUVFBcHAw5eXlHVq7q6ioiKioqG57wurufezu/QPpo7jwdHQ8vxjeP9LHC1937x9cHH28mEgsP/+kj91Dd+9jd+/fxUbyLOdfd+8fSB+7i67cx67VGiGEEEIIIYQQQgghhBCii5FEuhBCCCGEEEIIIYQQQgjRAkmkCyGEEEIIIYQQQgghhBAtkES6EEIIIYQQQgghhBBCCNECSaQLIYQQQgghhBBCCCGEEC2QRLoQQgghhBBCCCGEEEII0QJJpAshhBBCCCGEEEIIIYQQLZBEuhBCCCGEEEIIIYQQQgjRAkmkCyGEEEIIIYQQQgghhBAtkES6EEIIIYQQQgghhBBCCNECSaQLIYQQQgghhBBCCCGEEC2QRLoQQgghhBBCCCGEEEII0QJJpAshhBBCCCGEEEIIIYQQLZBEuhBCCCGEEEIIIYQQQgjRAkmkCyGEEEIIIYQQQgghhBAtkES6EEIIIYQQQgghhBBCCNECSaQLIYQQQgghhBBCCCGEEC2QRLoQQgghhBBCCCGEEEII0QJJpAshzlpWrokVq1LJyjW1uF328eI2bSeEEKL9mLKySF2xAlNWVmc3RQghhBBCCNECU2EWqRtWYCqUa/euxNDZDRBCXJiyck0sX5FCVm4x5jIL4SH+JE9JIL5PZJNt123NYPWnuwFYtnhKRzdVCCEEkJGSwq7VqwGYsmxZJ7dGCCGEEEII0ZyM9BR2fXPq2n2mXLt3FZJIF0K4ZOWaSEnNaDYh3lhKagZZucXE94kAYPWnuwDPifJrJycAGslTEs57m4UQQrRNQnKy299CCCGEEEKIrikhMdntb9E1SCJdCOGSkprRYkLcKSvXhLnMwuxJl7JkbpJrubnM4irfkpKawbWTLyXQGwbERrBs8RRXKRhnQr2tSXshhBDnLjI+/oxHopuystifkkLvKVOIiopqp5YJIYQ437JMR0nJ2EhywgziI/t1dnOEEEK0gakwi4z0FBISk4nsEd/qSPTTtxftTxLpQggXZ4K7uZHjzhHr2ceL2bA9kxuvGulKgoeH+LPq453sPpgHQH5ROaBYNGuoa//GiXpoeRS7EEKI9mfKyiIjJYWE5GQi45tefGekpLD77bex+fkxMDHxjPcXQgjROVIyNrJ611oAlk25s5NbI4QQoi0al3NJSExuNUku5V86niTShRAu8X0iW0xqOxPhkWEBbsudI9TDgv3Y/0M+Br2OIQN6cOS4mf+kfMucWUkMjIvymKiXci9CCHF+NJfUPrR5M1teeIGI/v2Z/OCDrnWmrCxSli+n+NTkownJyU32T0hORgG9k5Iozs7m+3Xr3NanrVzJnjVrsJjNzH7qqY7tsBBCiCacI9ETo4cA80hOmNHZTRJCCNFGjcu5tCVJLuVfOp4k0oUQbeZMeicOiib9UL7r/ys/TGPN53uIOpVg7xcTTv+YCNZtzSChXzBpB4t45kH3wN5a0l4IIUTLTk+cOycTtZjN+IeHu5Zve+UVTuzeTUFGBhEDBjBl2TJMWVm8f++9FB48SFi/fhRnZ/PBvfdSnp8PuCfVo4cPZ8c//0nZ/v1UnjyJxWwmackSMlJSqCkr69wnQQghLmJZpqOsTFsDwJKkG4mP7NdoJPo8t5Homw99xVMbXqavb28WTZjH9MFXtnhcKQsjhBDnV1vKsDQu5+JMjkfHJJK6YQXRMYnk56Vz6fBrgcAm24uOIYl0IYRHm785xCvvbOOKEX3Z9M1h+seE8+Ctk121zrekHWblh2lMTRrI7oN5WG12yqtq8DIaGDqgJ2n7jlFXb8Nqs5O+9yh/eXMrA2IjXOVckqckuNVIP5OJToUQQuBKnANMWbaMhORkLGYzebt3uxLiU5Yto+8VV3B89258goPJ27uXT37zGwAKDx7EWlND4YEDFH7/PQA+wcFkrFtH3t69HPv6ayxmM4WHDlFSXk7ZkSPodTrgx5HofceNI3r4cAZOndoJz4AQQlzcUjI28tZ3a6mz1bH18A4mD7yCqQOvwGyZgtlSSpbpqGu7dRkbOVhwGF0IrNj6b6YPvtItYe7cLjlhhpSFEUKIM9SWJHlGego7t6/iSNZ2kq9/1rWdc19notw/IJL9ez9iwuSlRMck8tGah6ixlNGj9xAsFSaUgqGJizqye6IRSaQLcZE5PWHdXAL7lXe2sWPPMfb9kE9FVS0Zh08SEuhLeIg/5jILb/9vN3VWG59++T1Wq53IsACqa+oZOaQ3IYG+5BWUAeBwgM2u2HMwj5BAX5KGx2Eus7DywzTWb88EGpLqy1ekkJVbDEjNdCGEaIuE5GS3v0tycsjcsIHaigrC+/enODub/9x6K1mpqdjq6rBaLBwqKADANzQUL39/rNXVbsesLS+noLwcU1YWDquVw6mpXHHXXXz9/vugFEY/PwZOncrhLVsAsJhMVJlM7Fu7ltQXXiD8VPkYQGqnCyFEO3GORC+rqSDAy49qaw1ZxcfILTtBWU0FBwoO8UPREVIP7wAgtzSPYJ8g1/7lNZX8Z+f7vPzl61TWVbHx0JeUV1eQW3aC7OIcHpz8MwASo4ewIvVfMjJdCCFOc3rivLUyLKbCLCxVZvz8wzh+bBcfrL6XKbMeZt/utRzM2EBdbQWgoWkaBoMXVmsN1VUlVFeXUlGWj6bpiIjoz+ChM+jVezj79nxEnz7xZB1Kpaa6DF+/EJLGL5EJRzuAJNKFuMg4y7BkHy8mJNCX3QfzTk0M6p7AnjN1GLknS5kx7hL2ZObTPyacg0cK+HrvMaLCAwCF1WqnoLiSiBB/hvTvAYCptIqwYF/8fL2orql3Ha+k3ML7G/Zitdux2RzE9AghMsyfxEHRpKRmkJVbTHyfCKmZLoQQHniqfx4ZH8+UZT9eqG975RVKc3MxeHkREBlJ5oYN1FVWohwO1zbOf1ebzWh6fbOPZ6+rA8B89Ci5aWk47HagIdH+9u234x8ZyeCZMxk+bx756elkrFtHwfffk7d3b8N+R45QcPAghzZuZMrDD5Ofni5JdSGEOAuNa55vObzjVKL8B34oykbTNJL6jKQuv57Kegt1tnpSMtZTZ2u4Bs8qPuo6jsliRkfDr4qOFefy8Md/RKHwN/qyP/8gNrsNB4oj5hziI/uxbMqd/OaTZ3l790dsPPQlf73+j5JMF0KIU05PnLdWqzwjPYXv0z/B7rBjs9ZyPGcv7666E6u1BofDdmorhVJgtdaiaRq1tZVUVRQBEBDck+LiI8T1H8vWDSuordNzaL+FooJDOBxWdHojpSXHKSs5js1WR0yfkUye8aAk1tuBJNKFuEgdyC4gK7cYvU5j9KWxrgR2Vq6Jv7y5le17jlBlqWPj1z/wzIPXAnDLI//F7lCcNFW6juNwKIpKqvh8WyYaoIDdB/JAgab9+HhV1fX4+hipstSjAVm5xeg0jd/+9X+MHBLD7EmXsmRukmuU/MoP0wBcy4QQ4mKWkZLCzlWrOLRxIzGjRpG0ZEmTpPSEpUvJ27OHeouF47t3U2exoE47js5obEiKOxzoDAbspxLknmg6HQZvb46lpUFwsGu5taaGstxczJGRDJo2DYBtr74KNJwDnEl0W00NBQcPNiT4c3Ka1G4XQgjROmeZle1HYtlzIoOa+hocyuGK7zuOfYdd/fiFqTOJ3hI7DtSpI1RbawHw0htx4GBk9KWu5H1ZTQX1tnr25R9kZdoanpq9/Lz3TwghLkTRMYkcidhOdEwi8GOtclNhFqkbVrgS6s6SLcVF2dgddmosZtSpmF1XV9nM0RsS6uWleae21bBUFlNZlk/esd2gaYSEJVBrK8XhsALgsNv44cAmV1K+pPgovn4hzJ73VLs+DxcjSaQLcZFZMjeJ8BB/so8Xk5NfSlx0KIP79SAnv4SU1Ayyjxfz/vq9rovzCksdT/5zAz7eBmx2R4vHdu6jTv1DO21dSIAPkaEBFBRXUFVdj0Mpso8Xc7ygjFuTx7gS5impGaz5fA8A4SH+UupFCHHRS0hO5sj27eTt2UPhoUMUZmaS/OyzRMbHu41W7z9hAj9s3EiN2YzykCR3WK2uf9sb/dsjnQ5lt1Nx8iRBjRLpTicPHOA/t95K7s6d1JY3/LLJ6OPDwOnTKT95kmqHg/B+/Qjq2ZOgnj05vns35iNH2PX221z7zDOuJLwQQojmJUYPYV1GIDkleXjrjVQpi9v6xkn0s+FMqNfZGxLwb373Pl9kf81R83GCfQPpExpDQWXROT2GEEJ0N/l56ZQW55Cfl86goT9e0zpHqluqzBQWZFJcmIWffxhFBZk0/DC0+UEsp1Ou+K5w2Otd/3ZmWiorTjbeutHIdg3Qsfe7DziatYOkCbdjqTK1WL9dtJ0k0oXo5pyThi5dMIFpYwe5ls+bPtxV2uWTL75nx94jHM4pxtfb2GQE46FjhXgbzi1cWK12ck6WYdA31P1y8vEyuo1ch4aa6eYyC2WVNZjLLGz+5hDph/JlIlIhxEUrMj6e5GefJW3lSg5v3cqxr7/mX7Nn0+fyyyk7fpzirCy++Otf8QkKAk1DqdMjuQeOVr4ctdmw2WxwaoLR09lqasj8/HP3Q1qt7PjHP1BK4RMURKXJROXWrYT26UPhwYM4bDZqy8vZ9sorkkgXQlz0Gk/26Syb4qx/DjB14BU8ueFlDhT8gEJh0JovyXW+1NutHC4+BoC5uowwv1DumbjENSGpEEII3Eq5NJ4s1FJlpkf0UNJ3rXXVLi8rzcPhsOM+1LB9KWWnprqEmuoSvtz8smu0o6f67eLMSCJdiG7OOWkowLSxg0hJzWDVxzuJ7xOBv68X+3/Ip0+vEE6aKqm32qm3Nv2G1GZT2GytjFxsI5td4Ry7rtdr1NZZCQrwoayyhqxcE/F9IonvE8lT98/mNy99wprP97D7YB6mkipAJiIVQly8IuPjmf3UU/x1yhQcNhvVJSUcWr8e9HqUzQanJhQ9E5rRiGptZHob6b280BuNDaPTNY26igoAQuPiCIiMJH/vXjS9npDYWCYsXXpeHlMIIS40jZPnzrIte/O+52DhYe678nYyi7JZtXMNSik+/X4TRZXFrlHj/l5+WOqrsam2j2g8V9nFOWzL/hazpZQlSTdKnXQhhMC9lEvKB8spLswiokc8pcU5WK21VFvMAK6/G7RhoMt5odwG1ZSX5RMU3IviomwOHdjM4cwtADI56VmSRLoQ3dzSBROorrXSMyKIrFwTiYOi0TTY9f1xfLyN1FvtHM0rweHoqKD+I7u94THLq2pZl5rBkTwzo4bEMLhfFB9t2Y+/rxcA/WPCmTFuEImDolmxKlVGpgshLlqmrCxs9Q0/7dR0OhSg07Szviw/X0l0aFR/HVyjXnRGIz2GDMEnMBCjry9RgwZx/V//KjXShRAXrZSMjazauYZ1GRtBKZLiRvFF1jcUVRXzu89eoG9YDHblwO6wU1hhws/Ll8p6Cxoaer0eTdPQKx12zq2kS1s5cLDtaBppObvYfXw/f71BJh0VQginjPQUCk9m4uXlR2TUIErNubTlh6EdyWG3UlaSy97SPLIOfYG1vhpNp8c/IFxGqJ8FSaQL0c3FRYfh52Nk687DHMkzU1tnpchchcPhwGprSHjYOyGJDmDQaziUwuGAunobu74/zv4f8vHy0mOpric0yI+Yng11eZOnJJCSmsHqT3cBP45Mb650jRBCdEcZKSmUnzjx42zODgeOVkq0dBSrxYLB29v1f53BgKZpZKWmkpCczMR77iEhueFnsKkrVrj+7azvLsl1IcTFIDlhBtuPfMs3Obux2q2U1pTjY/ABoNZaR3ZxjmskoR0H1fU1QEMt85LqMvyMvlQ7apocV6fp0Gs6rK4aueeX1WFjX/4Bt0lHPZWmEUKIi4GznIt/QCReXn7U1lSQvus96mor8fULQdN0gA6lziUma5zXUezKQVVlEX7+YQQF93JNlCrOjCTSheimsnJNpKRmYC6zkJVbjNGoZ/8P+TgcjlPlVcBh67jki6bR5JtZm13h622ktt7qWmez/VhepqKqFktNPUfzSti5P4f7brqSRdeMJnlKgusYp5euEUKI7qLxJKLOJHNCcjKHNm7kRHq6a6ohdXoi3VPAbatT++qMxoYE/akR5prBgN7Li+DevSk9dszjRKYAtro6178dNhteAQHY6uo4eeAAkx98kMj4eFJXrGDX6tUAWMxm9qxZg8VsZvZTT51dm4UQ4gLTMyiSuNAYKmoqMFeXUm+rx6AzMCiiP8fKjrtdo58+8rzGWote07kmGXWmWRzKQXuXDbA57KQe3s7d7/2aEN9gymrK2ZD5BWZLqSu5LoQQF4O07SvZ8+0aInrEU19X7ZoYNCQslpCwWI5l7TjHJDq0JaZrmh51BuW+NE1HXW0VZY68JhOliraRRLoQ3ZSzFrq/b8NkniGBPhQWVxLg543VZqem1tokLOv1mqvcyvmk10FEWCAnTZUAGA06rKduEGx2O34+Rmx2B3X17icAf18vIsP8KS2vIfdkKa9/mMa1kxPctlm6YILb30II0V1kpKSwc9UqjmzfTvKzz7qWTXn4YQ5v2cLhrVspPny4YWNNQ3PWSj+HJLpOp8Nht+OwWvELC6O2spKwvn2JHTsWFR9PcXY2FUYjtmYS6aerr2qY36Lw4EFenjiRhLlz8Q8NZfCsWUQnJpL6wgs/loMRQoiLQErGRjZkfoFdOegZGIm5uhQF2Bw2skpyXMkYgLiQ3gR4+xMbGs2+/IPkVxQCEODtT0VtFQpFqG8ItbY6qq01hPuFoqFRZDFj1BlQKGwO9xiruf7W8DF4UWez4mhjmRiFIrs4h6PmXHy9fBkUOeB8PCVCCHHBcI5Er6kuAyAioj9eXn4UnsykR6/BTJi8lC0bXnCL5edK0xkIC+9LbW0FoaGx5OWlu9adSRK9YXuFUg7Co+JdE6aKMyOJdCG6qeQpCWz8OpM9B0/gcCjKK2twKIgMCyA40Jd9h064ktlOUWEBFJc2TFR3+rrGzvQHRnYHFJqrXP8PCfKluNSCpmlomoal5scavc4qMxpQXVvP8QIbsT2DqappqAncuLRLVq6J9EP5PLssWWqmCyG6DedI9OjERCLi4yk+9X/n6O2RN97I7Kee4tWf/MS1j9HPr2Gi0VZGo+uNRrwDA9EZDFjM5obR7EqBphHUqxdVJhPY7Wg6HUG9euGw24mfPBnvnj1Je/ttlNVKr0svxeDry8n9+7HW1GA/VbO9xcd2NKRp9n/wAb6hoXgHBlJTVkZ5fj4xI0cycOpUV7kXKfEihOhuGpdAUQ4HDuUgyMu/oYzLqatqDQ2lHIT4BFFaU46maVzedxQDIuJITpgBwA2v30l+RSFVdQ3X62G+IUwbNIG16f8DoNZWx++veojX094BNJRDcaDoBww6vSuh7ozSXnoDaBpeBgO1tvo290WhUGjEhfbm4am/YMvhHQBsPvQV6fkHpcyLEKLbajyx6KWJs5kw5S6iYxI5nLmFHj0HE9VzMOs+eITq6lJ8fIOprang3H4ppNEzeigzZz/WcNyqYkLDYwkJ6d36npoOpTTg9ES7wuGwcTJv3zm06+ImiXQhuqn4PpH0j4lg94E8NA0G94sCTYeppJI9B/NwOBQaDXkPZ/K6rt6Gzd4ww3N0VBDDL4km/dAJSspq0Os1qmsbEt4+3kZq6qxogE6nuWqsaxoY9DqPSfhAP2/Kq2oBqKisI8DPG18vA8Vl1W7b6fU6Egf1orKqDnN5NbV1VkYOieWGmSNJHBTNlrTDp8rVmDzWTBdCiAtdRkqKq/RJ8rPPusq7bP3LX7DV1lJTVgbA1IcfZu3991NdWkrUoEEUZGRgr69vmITUQ7kX74AAxi9dik7TiE5MZOeqVRz9+ms0oNewYQT36sXhL77AUlREYM+ezHzsMfLT07n02mspraqi7oYb0ICkJUtIWb6c2ooK9EYjEQMHUm+x4BMSQtGBAy32TQH1NTVUl5RQX1PD8OuuI2nJErc+T1kmkx4JIS5sp9cO/8vW/yMlYz3ZxTnszNlD5alEuE7TsJ/6AtKoN6CAS6LiKagsoqKmgpPlhWzI3IrZUsqSpBsZHj0Uk8WM1d5QLqBXcBSgERUQSbHFjMOh+Gj/enyMPuw98T09AiII8PLHz8uXitpKHPYfzw02hwO73YrONUb9TChKLGW8kPpP+ofHkZazi8zCLHJKjwOwbMqd5/L0CSFEl5SRnkJxYRYRPeJJGr+EyB7xpG5YQWbGekaPXcTutNWUmnMxevnhsFtdg0w0TYdOZ8But9I4sa43eOPtHYDe4EVl+UnXcp3OgNHbH7utoWTihk+epKqiiMCgngQERGI6eZjAkB6A1qS0S0TUQGpryqmqLAI09AYv7B6+LHXYrWzd+Bcmz3iQjPQUEhKTiewhg1naQhLpQlwENDTKq2oJDvTjpKmiIYmugZ+PEUuNFU2DyFB/Rg2N5cvvsgkJ8uXPv5rDtLGDuP6B19hhPkbCJb0pr6zmaF4J/WJCKausxd/HC3NZNXaHA0tNPXq9Dp2meUykV1rqMBp0GPR6Avy9MJVYqKm14jht9KLd7uBwTjF19TYMeh2jL41l3vThpB/KJy46jPAQf1Z/uovwEH9XrfTGNdOFEOJC55yE0zk625lY9g0JweDjg29ICACDpk0jftIk9n/8MeV5eWg6HZqmgU6HTq/H6OPjSqD7hYZSXVKCTtNcx9v2yivUV1Sg9/JC0zRy0tIYOGkSFQUFTFi6lEHTpjFo2jSKDh8mJy2NpNtuQ9M00lauxMvfn55Dh9Jz6FCGz5vH4S1bOLh+feudUwpbTcMkedUlJRRmZlKSk4PFbGbwrFmuvgshxIUsJWMjq3etxWwppaymnE8yNlFnq2fnsd14GYwEewcS4huEUW+k5FR5AJvDhkFnoLSmjKMluTgcdix51VgdNspqKkjJ2MjXx74D1TCxaIR/KEN7DmJD5lZqbXV4G70J8g5gz4kMJvQbw4jelxLpH47JYuZIcS5hfiEMCO9LMP707hlNWu4eSqrLCPEL5vI+I0nN2k5dG0emO5TiZGURBZUmABaNnkekfxgf7V9PYvSQ9npahRCiUzlLoTiTzqbCLCxVZkLC+7D9i39ht9ejaTr0egN2Wz1eXv7U11WhTv2mPyQ8lprqcpTDgZ9/KD4+gVRbSvDzD6OyvABnkv2SodO56trfkfLBck7k7MFmq0PTdPQfOB6z6QgKCI/qT31tPkqzU1nWsG9gSC+unvMH9u1eS0b6J9isNRiMvthtNvBQwisvdw8pHyynMD+TI1nbSb7+WUmmt4Ek0oXoxkICfTEa9NjtDnJPluFVXMUlfSMxlVioqq5j4ugBmEqr6B8TwYO3TgYaaqsnT0lwlUppXIN8S9phTCV70DQdlup6xo/sT0igL59++T0VVbX07hHMgp+M5D/rvqNnRCAZPzR8q6oBdqVw2BTDB/Wmts6KqcSCze5A00DTNFdy32jQU1dvw253nKrfpdiSdpiPNu/j7f/tcptwNL5PpIxEF0J0O42T540lLVmCf3h4k2SzstupLi3F6OND7GWXEdyrF74hIQycOpX89HTX9s5yMc4SKsPmzKE4O5vo4cO5fPFi17anl1ZJe+MNsjMyqDt+nIDwcPasWQPAhLvuYsqyZaSuWMGeNWuoq6riTCiHg+O7dvHRQw+BUly+eLGUdRFCdAvJCTMwW0rZeng7WafKt2hoFFtKqLHVYdQbOV5+kiFR8fQL70N5dQX5FYUM6TGQ/hFxHDXnojfqsNptaJpGiG+Q65hlNeWE+AazJOlGAEJ8g1zLymoq2JC5ldjQ3oyISWDVzjWE+YXia/SmrKaCjJMHeTDpZ9w5/Vbu/eC3pGSs5/I+I4gNjcbfy6/NiXQnH4M3vgYfsotz2HjoS/LLC0jPP8i0QRPb42kVQohOFdkjnikzf7xGz0hPITNjPTU15dTWlAPg5R1AfV01mgZh4X0wm44CCrvdhq2+lrqaCgwGL8Ii4ijMz8TL24+QsFiKCn7A4bACOtexJ0xeyuHMLdRUl+HrFwJAeVk+MXEjCA2NJSfLgsNhJSi4F5YqE7WWctZ98Ah2Wz0Rkf0pKsjEYbcSGNKDyrKTnK687CQl5hz8/EIpzM/kg9X3Eh7ZH1+/ENeIe9GUJNKFuAA5y5o0Tnh7Wj81aSAAZZUNo/9CAn2ZmjSQLWkNk9MtmZvUZP/TE9PTxg5i2thBrv9nHi3E39eLvIIyQgJ9CQ/xp7bWRlCgD/fddCW3JF/OQ7dN4zcvfcL+H/LRNPD386LC0nBh7uttYPrYgRzOMVFvtZ8qp6uICgugrt5KpaXh21a7Q2F3KLbtPkqAnzf1Njum/FJ+/8rnXDVhCDn5JaSkZpA4KJr0Q/nNPhdCCNEdOOumn57o9g0JwejrC5pGvaVh7onJDz5IRkoKYXFxDJo2zW3fxiVUAPQGAzEjRrhGn7ektrwcDRg8cyYAFrOZQ5s3N4wmnzmT2spKju/eTbXZjDo1gajHMjONWKursdfXE9y7N9GJiefwDAkhRNeyO28/WcXHUIBRZ8CgN1Bvt2LQ6THqDPQMjOSxWfcTFxbDyrSGLyidyfEBEXEkRg9x1R9fknQj8ZH9eGr28iaP03hZlumoW0317Ue+Zc+JDACG9BzIqN7DSIobCcCDk3/GgIg4zJZS1uxJaUjao7lqtrdFja2Wnbl70DQNo97IyN4JrscWQojuLiExGUuVmdycXRTmHyQgKJKevYZy9PA2onoOYsqsh9m29RVyj+xEKRvBYTEEBDbkLPr2v4LiwmyqKk3E9bucpAm3cTgzlVJzLmUlx8k/nk5Ej3jXKHFTYRZp21dyaeJsLr/iNnJzj3AiZwD+gRH4+ARSW1tJ/vF9rqT70OGzMRp9CI/sz6GDmz2231rfUGZMKYXNWsvxnD3k5e7Fy8sX/4Bwty8NxI8kkS7EBSYr18TyFSlk5RYDnmuDN64d/tT9s137rfwwjRfe2EJ+UQWLr7v8jBPP6YfyyckvZdb4wYxYEOMqqbJ9zxGycov5Zl8OH23Z7xrF7uttZECfCG6dO5Tfv7qeKksdBcWV/GPNDmw2O95GPfVWOz0jAgnw9eZw7qnRjKfNcJ2dW0xNbUMJmkpLHWs37eNInhlTSRUbvz5EVo4Jc5mFp+6f3eqXDEIIcaExZWWRsnw5xVlZgHsNcecodf/ISHavXk1Qz56krVxJ5qkyK1OWLXMlzy1mM4CrhEpJTg5Htm9vNYGddNtteG/eTHlGhmuyU//wcHatXk1hZialOTkMnjWLiAEDqKuq4tiOHYTExREYEcGJffuaT6SfWu6w26ksLOTwli1NEv8yQl0IcSFambaGfScONJrY04jFWoOmafjqfRjUYwCjYoYRFxZDSsZG1memsmj0PNcknc4a42c6sjs+sp9bffKlE2511TF/cPLP6B8eR1FRkdu2mw99xe68/UT6h5N+4ntOVhad0WNa7TZiQ6KptdXi7+V3RvsKIcSFLLJHPP4B4ZiLsjEYvLgsaREJicmkbV/p2iYoqCc6nQGbrY6ykjzCwmIpL8vn2JEd1NZWoBx2fP1CmD3vKT5Z+xv2fLuGntFDqagooLgwi4z0FKbMXOYa/T567CIiogaQtuMjTEXZKE1RmH+A0Ig4QOHrF0JtbQU7vvgHmqbD6O1HeHg/8iwlzfbDUmXGWfpFKR2h4XGuMjaiKUmkC3GBSUnNICu3mPg+Ec3WBvdUOzwlNYM1n+/B7lCMHNL7rOqKNz5u4yT1s8uSSUnNYOPXh9h78ATQUArm0LFCJo0ewMkyC79fehUfbdlPSXk1FVW1BAX4MC6xL1/vPUZ0VAh7Dh73+Jg6Dczl1dRbG0Y36nUaRr2O/jHhzBg3iOzjxRw6WsSOvUf4zUufALB+eyYgE5AKIbqHjJQUirOyiIiPb7aGeL9x47CYTOxctYrg6Gi3euPOvy1mM5nr1zN60SIi4+PJSEmhNCeH/PT0FkejRwwYwPDAQL49/mOcdh4zOjGR/PR0irOz2bNmDTq9HuVwYKupYcrDD7P5+ec5sXu3az/NYEDZbO4PoBS22lpKTx1fJh4VQlzIskxH2Z23H6UayrnoNR19w2Kps9ejASNihhHiG8T6zFTC/UNdI7jbYyR3ev5BTFXFzBh0JfGR/XB4+GLTuY2f0ZdaWz2apqGUItDbH0t9DY5TA1yaG62uUOSWNVz/bzm8DUt9Nc8m/9r1pYAQQnRnzlHpANExiWzd+BeOHN6O3W6lsCCT/OP7UCh8fIOpr62isOAQMXEjmTB5Kft2r8VsOsLAwVMBSBq/BP+AcFcS2zkJqPNxnI/x6Ue/JT8vF2t9DUajLwGBkXh5+XNp4mwGDp7Kug8eodSci15vIOfITjRNB5oeH58AbNY6bLba03rREOf1Bi8Mei/6xV8hZV1aIIl0IS4wzSWz27KfuazhpzueSrq0xek1yRuP/l62eAqJg6J54Y1UekYEsSXtMDn5pXyzL4cdGYXMvGIw40f2JzLUn4+27GfO1GFkHi0i9qpQpiYNZPGv36LeYXd7PJ1OIzTQl/KqGjStYZlS6tTI9FqSpySQk1/Czv05HM0rIa+gnBuvGumqoS6EEBcaT6OxT598tLHGSeeE5GSObN9OcVYWg2bMIDI+3u14gFuN9cZ/t2UUeNJttxFwan9nHXdTVhb56emubfzCwqgtL6euqop9a9dSkJHx4wE0zVXyxZO83bv55De/YeDUqW7tE0KIC0GW6SgpGRvJLj7GwYLDeBu9qamvwabsZJmPcXmfkeSUHneVXnEm0U8fRX4+tSVJ71yXGD2E+9c+jlIKvU5Pz8AoegREsO3YtwCtlnzRoRHg5U9mYRYpGRvbrU9CCNGVRPaIZ/a8pwBI3bCC/Xs+xm6vJyg42lXjHGDg4Kmufzvrj+fnpZNzJI38vHQGDZ3WpAZ7438716VuWMHeb9/HL6AfmqZhq6/BVJSFqSiLmD4jGDR0GvAM6z54hMqKQmzWOkA1zLcRGkNB/sFm+2I0+jJyzI0kjV/ittx0amS8c5LVi50k0oW4wLRlgs3GpV2c28b3iXSVeTlfGj9O8pQE0g/l0z8mnA3bM5k5fjALrx7FpX2DiOvTG3NZNas/3cWs8Q0J9cyjRazfnsmia0YzbewgAvy8KClvqOXeMKd1w+hzm92B3f7jpbsCaupsbPr6EJaaegb360G91c6QAT0YNST2rL8kEEKIrsDTaOzTJx89tHkz2155hQlLlzZJsic/+6xb4vz04zU+TuPjpq5Y0eoo8IgBA1zrnIl35yj3wbNmMeGuu/CPjGTX6tUEREbyw+bN2OvdJ67z8ven/tSkpHqjEbvV2rBC06guKWH322/jHx4uI9GFEBeclWlreOu7tTiUw1ULPcI/jCDfQEbGJDBv+E9Izz/Y7snzxtryOM5tskxHOTV5ET56L46U5GKpr27zYzlQVNVbGNZriNRJF0JclBISk9mfvo7igsP0HzieQUOnnUpsN2j8b+f2jf9u62NUVZmpqzPi7T2aSxol6J3HGTR0GmHhb5O2fSWlJcfJy9lNXV0Vvv5h6HR6j79OArDW15CX0/BL0oGDp5Kfl05CYjIZ6Sns+ubUfYLUTZdEuhDdkafSLu39OM6kemRYANAwsekDt0xm41e72bH3KNdNGUZ4iD/mMguvffA1RqOe8SP7kzwlgaxcE5cPi+OrXdnU1ttRDgUoDHod5VWn/+yogV6nJ/NIIdW1VmZPulQS6EKIbuH00eKebHvlFY7taJiA7o4PPmg2Od7W453Jdk7OBP3gWbMYPGuWa9+MlBQq8vOxFBdTW1HhvpNSDUn0Uz8xcjgcoGnodLqGL0s1jZC4OCxmM6asLKmPLoS4IDhHopfVlFNnq8OuHOg0HSE+QVTbarlu4FWuCUHPtOZ5R8gyHWVl2hp25+3HpuxEB/Wgf3gcNbZalENxsrIITYFOp8PmaP5XRdBQL73WVtdBLRdCiK4lskc8C2/9t1tJlta2P9PEdGSPeK6Z80eKioqIiopCp9O5Jegbjx53jpQ/dGAz27a+gsVSgsNhbfbYDoeNvNy9FOQfJPP7DdhtDQNizibh353pOrsB7eXvf/87ffv2xcfHh6SkJHbu3Nnstv/+97+ZOHEioaGhhIaGMn369Ba3F6Krc45ab+/EcuPHSZ6SwKJrRvPwbVO4a8EElsxNAuB/Xx7k6z3H+GjLfpYtnsKSuUn4+XpRUFxJQXEF8X0iSUnN4EB2IYuuuYx504cTGuQLNIw8b05YiC9Go56D2QWutojuR2K5uNg4E+EtJZEnLF1K3yuuYMLSpefleGeynVNCcjKjFy1yTXaauX69ayR8RHw89dXVhPXti85obLrzqRGPym4HpfAODsbg7U3UJZegAd+99Rbv33svplOTq4oLn8Ry0Z2lZGxk9a61hPgGM7jHQAw6AwPC+xDkG0SdrZ6ymnKyTEdZkfqvhlHfXUxKxkbW7EnhUGE2l/YcxDWXTie3LI8Zg65E02mgoGdQFAZd62PwFIoDBYf4w+cvtn/DRYeTWC5E65zJ8c4qgeIcPZ6RnuJadjhzC/nH9xEaGktIWBw6nZGGOgDulHKglAObrY6qShOgUVyUDdCpfepqumUi/d1332XZsmU8/vjj7N69m8TERGbNmuWaofx0W7duZeHChaSmpvL1118TGxvLzJkzOXHiRAe3XIgLlzOpPm3sIFdyPft4MT5eBnpEBDBn6jDXds88eC3jR/Zj6YIJQMOI9lnjB1NWWcPeg3mYy6udeZYWmUot2Oyef5YkLnwSy4XwbNC0adzxwQctThDa3hon3p1Jdedo9h6DBzN8zhxiRo7Ecaomut7Ly/0AjYJ8TUkJ9vp6aisrKc7OxlZTQ+HBg2SkpCAufBLLRXeXnDCDWYMbSik+NvM+liTNJ8g3iPzyAnSaRohvsCvZnpKxsZNb21RywgxuHJnMwtFzeDb51yxJupFZg6dgtpSyaNQc4iPiqLXWYWxDIh0ayjB+kf11l/zSQJw9ieVCXBgSEpMZPXaR2+jx0pLj1Nc1lFYckjATTXMW0/VEoZQda30N1ZYSMvamkLZ9Zfs3/AKiKdWWdNWFJSkpiTFjxvC3v/0NaPjpcGxsLPfeey+PPPJIq/vb7XZCQ0P529/+xq233trq9hUVFQQHB1NeXk5QUNA5t78tHA6H2085uqPu3se29K/xZJ5tGXF9ptufb1m5JlZ+mAY0TGiakrqfL77eR2ZuJb+YP8FjbXdnm81lFtZ8vocKS22LSXRNa8i/BAf4YHcoBvWL5K+/vqHTRqR39/dpZ+roWA4dH88vhveP9PHC19b+mbKySFm+nOKsLC6dPZujO3ZQ8P33AGh6fYsTjTZspOEfGUl9VRU9hwzh+r/+tcPKu3T317AzSSzvHqSPnjWeYPSzA6kM6jGAUTHD+OT7jUQH92RUzDCWJN0INIz8dtZH7yxt7eOK1H+xetdaFo2ex9u7PiKnNA+DpsemWonjjfx0xLX8/adPnY9mt9nF8D7tLBdDLIfu/x7q7v2D7t9HU2EW+/em0DtuCqEhAXy/b12rk4C++uJPOJGzm6jooWgKCk8eRKnWByT6+YdTV1vJsJHX8dOb/34+u9Gqrvw6drsa6fX19ezatYtHH33UtUyn0zF9+nS+/vrrNh2juroaq9VKWFiYx/V1dXXU1f1Y+63iVA1Qh8PRbNH+883hcKCU6rDH6wzdvY9t6V9K6n7e/t9uQPHALZNbPeaZbn++paTu5/31ewAID/Fj2MBeHM7OoU9Mb66dfGmTvmYfL+bRF9eRnVvM1VcO5carRrD6k+9aLOkCzmS6g6snDiU40Iev9xzhkb+k8MsbxzM16ZJ2658nnfE+7WonkvbQEbEcOj+ed/c4B9LH7qCt/dufkkLxkSNEDBxITXk5puxsMBgavv3U6Wj1p0aaRuzo0dRbLIz/5S8J79+/219Xdfd4LrG8+5A+erbymzW8v3cdvYN7YbNb+T4/kwHhfVg85kauTZjOgIi+rm0fmPQz1+N0lrb28dpLp2OuKsVcVcqCEcm8uv1Nqq216M7gB+2bMr/gcNERt+egvUksbx8XSyx3PlZ3jnXdvX/Q/fu4f28Ku9PewWr3I09XzZ6db6MUTJ7xAMVF2WSkryMh8Voioga49hk1ZhEVZYXYrfWUluShAE1nQKfpsbdQM93LOwCHw4GPX0iHP59dOc/S7RLpxcXF2O12evTo4ba8R48eZGZmtukYy5cvJzo6munTp3tc//TTT/PEE080WW4ymait9Twx4vnmcDgoLy9HKdVtT9zdvY9t6d+UUb3xM9hIGt672Z/Nncv259uUUb2x141w/TttXw5+XoqB/YIJ9HY0adPmbfvQ2Wu5ckQ0N0xtSIDbaqs4kF3QajIdoNpSTlFRERmZ2dTUWnn9vVTCAzR6RQaf9741pzPepz179uyQx+lMHRHLofPjeXePcyB9vBCVnzxJTloacUlJBPfq1ab+lZ88icVuZ+iiRVwyfTo/bNpEaGIigZGRhPTpQ3VJCcWHD2O3Nn+xHtKnDxVVVdiVIjcri5CE9p0wu7HOeg27ezyXWN59SB/dnSwv5OP96zl08hDx/nGMiU7kkoC+nKgoIM6nN4uGJoODTrkeb0lb+xiIH7HePdhyeBu9gnpw86Vz+CZnD3aHrdliAKfTgM17vyRwuN95aXtbSCxvHxdLLIfuH+u6e/+g+/exd9wUbHY/wqKGEeDvi1350TsuiaKiIvbt2Uz2oW+xOfwYPjLQtU95hQVv31hQDkIjwrBZa2iI0i1HdL3Bm54Rw+gVc2WHn8+6cp6l2yXSz9UzzzzDO++8w9atW/Hx8fG4zaOPPsqyZT/OrFtRUUFsbCyRkZEd+pMjTdOIjIzslsEBun8f29K/qKgoEi8d2OZjnun259vpj+8fGAIopo5PJCqqaemVaRNGUG0zcO3kBAbERvDif7ayYWcuJeXV1NW3nkjPL6mnb0w4R4+X4lCKkyVWUnef4IFbOu456O7v0wtVW2I5dH48vxjeP9LHC8+B1avZ9/bbGKqrGfjAA23q34HVqzm0Zg2jFi5kYGIiymymKC2N8YsWERYXx8r586nIz2/xcauPHcNaXU1ITAwjpk0jIiqqPbrnUXd7DbsLieVdh/TR3V+/XcXKve/iUA6ig6IYENufnbv30793LHOSfkJURMfFrzNxJn2cNuJKNh3fwZbjO4gMiAAvjWPm41Rb25bQHBDRl2kjruzQ5+JieJ9eiC6UWA7d/z3U3fsH3b+PUVFRDIgfhslkIjIykksGjXCtGzFyGgZdNQmJ7tfRI0ZOo67mOAC+vmF8tfllbPa60w99moZEe7nZmwEDEoi6bMb570wLuvLr2O0S6REREej1egoLC92WFxYWtvrtwgsvvMAzzzzDpk2bGD58eLPbeXt74+3t3WS5Tqfr0BdY07QOf8yO1t372Fr/Orvm+bmK7xNJkE9DEt1THwfGRZE8ZRgpqRkkDorGXFZNSJAfplILaBoOR8vfkBaXVTO4fw9+/tPxJA6KZkvaYcxl1RzJM3fo89Xd36edoSNiOXSNeH4xvH+kjxeWYcnJaEBCcrKrP6317/R9Tu7bR+mxY5zct4+T+/ZRVVgIrfw0U9ls+IWFccn06ay97z78IyMJjY0lacmSDqmV3p1ew65CYnn3In1svKGz1KCi3mHj9Z3vklmUha+XDwOj+ndMY89SW/s4MKo/zyQ/yvKUP7HnRAYAtfZ6HLT+M3u9puMPVz/UKc/FxfA+7WgXUyyH7v8e6u79gwu/j6bCLDLSU1qsfe6pj1E9BzK157Im20b1HEiPnoP5cvPLhITFAg50mobD0fK8F37+4VwyZCrDRiR3ynPZVV/HrtWa88DLy4vRo0ezefNm1zKHw8HmzZsZN25cs/s999xz/PGPf+Tzzz/nsssu64imCtGqlNQMVn+6i5TUjM5uSrtx9vGVd7bx0eZ95Jws4UzmQN65P5fIUH/SDzWMdFy/PdPt+crKNbFiVSpZuabz3nbRfiSWi4uNKSuL1BUrMGVldXZTiIyPZ8qyZWeUvD59n4TkZAbPmoXFbKayuBiHrfVfGVmrq6krL+eHjRs5vns3P2zcyJ41a8hISTnrvojOJbFcdFdLkm7ktqT59AqKothSAiiu6HsZSye0bRLFri7LdJQVqf8C4NnkX/OTIVPxM/qizqBWbXr+wfZqnuhgEsuF6FgZ6Sns+mY1Genndg1sKswidcMKDh3YzPp1f6DUnEP+8XQcDlurSXQATdMzecaDLU5kejHqdiPSAZYtW8bixYu57LLLuPzyy3nxxRexWCwsWbIEgFtvvZXevXvz9NNPA/Dss8/yu9/9jtWrV9O3b18KCgoACAgIICAgoNP6IUTylAS3v7uLzd8c4pV3trF0wQRX3xIHRfPIX9ZRXGqhV1QQwy+JZvM3P2C1NX/BrgH1Vju/+9tnWG12piYNZNE1o92eL2eiHmDZ4int2i9xfkksFxeTjJQUdq1eDcCUZU1HklxoIuPj8Q8PZ9fq1dSUl7c+0SjgGxbGoGnTiBs7lt2rV7tGpCckJ3dAi0V7kVguuovNh77ilW1vsnTCrUwbNJGnZi9n6sAr3JZ1FykZG1m9ay0AidFD2H5kJyWWUuxtGI0OYNAZWJexgcToId3qebmYSSwXouMkJCa7/X22nAn5IxHb0TQdPr7BXDHpl+TnpZN16ItTtdKbV1NTStr2lcye99Q5taO76ZaJ9Pnz52Mymfjd735HQUEBI0aM4PPPP3dNjpGbm+v204BXX32V+vp6brjhBrfjPP744/z+97/vyKYL4Sa+T2S3TP6+8s42duw5BsCzyxpODnHRYTzz4LW88MYW/LyNmEqr8PYyYLXVt3gsX28j1TX1KGDn/hxMpRYSB0W7Srt01y8jLgYSy8XFxJks7k5J44TkZCxmMwUHD5KbltbiRKMAg6ZN46d//zumrCyKTk1e1lFlXUT7kVguuoMs01EeWfc0x8safgHpTA5PGzSxWyaKkxMaauEmRg/hkXVPk1/RUNJDh4ajDdON1tutfF/wA4+se5q3w2KIj+zXru0V7U9iuRAdJ7JHPFNmntnAGk/lYJyJ+OiYRA5nbgFg+Mjr0Ok0jhzeRksTjnp5B2C31VFTXXa23ei2umUiHeCee+7hnnvu8bhu69atbv8/duxY+zdIiFZc6PXQz8TSBRNcfzceMZ48JQE/Hy/S9uVgtTkIDfIFmk+kK8DbqMdmt+NlNFBpqWP3geO88s42po0dBHTfLyMuFhLLxcXCWRqlO3EmwE/s3Yvey6vVRHrmhg28d/fd+IaEsGfNGgD8w8O73fNyMZJYLi50KRkbqaitJNA7gDnDZnV2c9pdfGQ/lk25kxWp/6K6vhovvZF6u7VNSXQvnZF6hxUNjcJKEyvT1vDU7OUd0GrR3iSWC9F1OUefA64kfOOE/OHMLez5dg2lJccpKzmOXu+FplU3W1bXLyCcitJ8CvIPYCrMkvIujXTbRLoQXV328WI2b9vHtAkjGBgXdVGVIJk2dpAr0R0XHQY0JNFXfpjGt/tzATAYNAw6jfg+4ZSUV1NS7vlnR0ajHm+HkbjoEGJ7hmIqtbgS9UIIITqHKSuLjJQUasrKcFitKKXQdLoW6+vWlpez9/33GTxrFoNnzsR85AjRiYkd2GohhGhq86GvWJexEaPOgEJhspR0dpPOqyzTUVIyNpKcMKPJyPHkhBlkF+ewN28/9TYr9XYrJTVl1LXwi9F6R8OXpgqFzWFj9/H9ZJmOyqh0IYRoR20tB5Oft4+KspM0NxLdqcycA2iYTUfJSE854xHy3Vm3m2xUiAvFuq0ZfPFdNuu2NkyMmTwloUl974uBc8S4cxS+1Wan3mrHanVQVGrhaF4JFVV1ze5fUl6Nn4+RInMVgf4++PkYO6rpQgghmuGs+e4bEkJCcjLeAQFtm0haKQ6npgJQZTKRn57ezi0VQoiWvbLtTQ4U/EBpTTl+Xn4kRg/p7CadV8566CkZG5usi4/sxxFzDlnFOdRYa6m21hAbHI3WxmP7GHzIryjweGwhhBAdJ2n8EiZMuYvLxt5CUHAvjF7+bdhLERzaG0uVGVNhVru38UIhI9KF6CTXTk7Az2Bj2oSGxLmUIIElc5PYsfcIB7MLUQq8vAxYrTbsjuaTL3aHoqS8mt49gtl7MI+s48VU11pdI96FEEJ0HOdIdOdI8oTkZDJSUtj34YdtmnAUwF5Xx5Ht27n0mmu6Vc14IcSFo/HEonOGzSK7+Bh+Xr5Y6qtJzz/YreqiO+uhO/8+XaR/GDqdjjpbPeV1leg1PS3V1XXyNXozof8YYkN7N3tsIYQQ54eztIulyox/QLhbrXRoKPOSkJjMB6vvpa62AqPRG2u9xeOxDEZfABx2K3qjF5kZ6/EPCJdR6adIIl2ITjIgNoLAacOJioro7KZ0KVeM6M/QAb0ICfRl14Fcdh840eo+/r5e1FvtBPh542U00D8mvANaKoQQ4nQZKSnsXLWKoOhoIvr3J23lSgZOnUpQ796U5eS0vLPWMMbR4OOD3WrFPzxcJhoVQnSKV7a9yfaj33KoMAuD3oDVbmXKwOmE+4d2u6Swsx56c2JDexPg5YelvhqAstoKVBtqpfcLjSM9/wCWes/lGYUQQpw/zpIulipzk1rpTmnbV5Kftw+lFHaHHZ3eiE5nwGZ1j9M2aw2azoBy2CkpOsKAQZNaLRlzMZFEuhCi0zknWjWXWVi/PdNV4ua/675t0/52h4PoqCAWXT0aU6nloiuPI4QQnck5Cj0hOZmE5GSObN/OiT17OLl/P5qmUXr8ODUlbagprBRGX18m3ncfOk2T0ehCiA6TZTpKyv6NXBo0gC+/XYW/ly89AyMxVZmxOez0Ce3NkqQbL7o635sPfUXq4R14672ICAknu+RYm8u65JQeJyIgnIyTmSxP+RPPJv/6onv+hBCiozgnFjUVZrlGpDuZCrPISE+htOQ4AHq9Eau1Bm+fQCJ6DOREzu4mx1MOO6CwWmsoPHlQJhttRGqkCyE6VVauieUrUlj18U7KKmuIiw4lMtSf5StSqLfa2nQMS3U9WTnFmEotbvXWhRBCtL+0lSvZ9uqrpK1cSWR8PMnPPsuohQuJHDgQTa/nyLZt1FVWoulav+y01tRweNMmpixrGEGTumIFpiypySiEaF8pGRt5e/eH/O/gFt7fu46vj+3imkunM2f4TxgVO4xnrn30oksCZ5mO8si6p8kqPkqRxUyNrYYwvxDsqvlJoxsbFTOM6voaHEqRWZglddKFEKIDOBPqjRPfzrIvVVUmjF5++AdFAhp1tZXk5zY3H5ECTYe3TyBXTruvQ9p+oZAR6UKIdnXSVM7q9QdInjLMY4I7JTWDrNxi4vtEEBLoS9q+HD7asp+s3GL6x0RwJM/cYo10AL1eh7eXnsjQtkyYIYQQoj1Fxscz+6mnMGVlkbJ8OSczMgiIiqKmvBx7XfOTRwOgaYT37w/8OGEp4EqsCyFEe0hOmAEKLg0aQEBIEGhclCPQG0vJ2Eh1fTVhfiHUWeuosdZi1LmnD3SaDk0BmoZd2V3LvfVG5gy/itfT3uVQUTYxUQO6XUkcIYS4UDhHpzscCktlMSFhsZSZcwFQjWJ3EwqiY0fQb8C4jmjmBUMS6UKIdpW2L4e3/7cP0DxOpuosw+L8OzzEn8RB0WxJO0zqzsOoNkxOZ7M7MJVa+GjLfm5Jvvy8tl8IIUTLkpYswT88vEkplsj4eIbNmUNpbi46o5GqoqIWj+MbFsagadOY/OCDmLKysJjNDJ41S0q8CCHaXXxkPx6Y/DOKioqYMWoKujb8gqa7cya+s4tz+Hj/51jrbUwZcAU7jn1HdX0NvUN6cbwsHweqybyj9XYbH+1fj4/RG7vDTml1Wcd3QAghLkLOMi6NJxt1jlJ/7ZXrKS/Nw2arR9N0qFZ/YeQgJ3sHGekpMtFoI3KFIIRoV0nD41h49ahm65bH94l0lWNx/jsuOozMo4UczTPTymB0oGF+On8/L64Y0ZcVq1LJyjWd514IIYRoTmR8PFOWLXObGNSUlcV7d9/N/x57jNKcHEqOHm31OPb6emorK4GG0eiZ69fLhKNCCNFJnJOQhvgGYbXbqLdb2ZHzHUoDnU7HyYpCj5OOeuu9GNpzIEsn3MrDU35Bn9DeWO1WKe0ihBAdwFnGJSM9pcm6CZOX0jf+Cmqry9qQRG8QN+AKmWj0NDIiXYguwjnhZvKUhG5V47tXZDAP3DLwjEb2pKRmkHmkEJ1Oa7Wsi0Gvw9/XiJfRwI69x8jJLwXwOPpdCCFEx8hISWH/xx9jr68HQNlb+NnoKfVVVRzauJF6i4UJS5cCyGh0IUS7yjIdJSVjI8kJM+gfHtfZzemSBkcNwMfgRZ29nh4BkeSU5mF1/DiPkYb7gHSHcuBj9CEuLIb4yH68HRbjeo6FEEK0r4TEZCxVZixVZkyFWW610sPC4+gfP56ysnxKirJbPZbB6EPPXoObjHC/2MmIdCG6iJTUDFZ/uouU1IzObkqnS56SwOD+PQjw88bLqG92O53WUNal0lJHn+gwli6YwKJrRjc7+l0IIUTHiE5MxCckBE2vR+/l1badNA3vgACO79rF4S1bSEhOJiMlRSYbFUK0m5Vpa/jbVyu59/3fkl18rLOb0yWt3v0RtbZ6egZGUW2taZJEBwjw8sega7hmtznsHCrKJiVjo9sXFRdzvXkhhOgokT3i8Q8IJzNjvWtUuqkwi9QNK0jbvpJd36ymurK4TccKjxoA0OwI94uVjEgXoos4vVb4xSy+TyTPLktm5YdprN9+kNyTZR638/E2UlNnRafTGDUkhmljBzFt7CCg+47wF0KIC0F+ejrKZsPLzw+Djw/VZjOaXo/Dam12n55Dh1Kak+MavS6TjQohOoLVbuWQKZt1GZtYNFR+BXO6SP9wdJqGn5cvkweM483v3qfe3hDLnSPRq+otru0HRMQxZeB4khNmkJKxkdW71gKwbMqdHd10IYS4KDlLsUTHJJK6YQWWKjP7dn+EXm8kOnY4WZmprR9E09F/wBUkjV+Cf0C4lHdpRBLpQnQRzvrgooEz+X2isKzZbYwGPVabndAgPwb3i2LFqlRX4tw5wh+kzIsQQnQkU1YWxdnZWOvqsFos1J2qe64czddiNPr7Y6urI3rkSKzV1QycOpWwuIYyC1LeRQhxvjlHSk8deIVr2bUJ06FtJWMvGlmmo5gsZnSajtzSE3ydswsvvRcGnYFqa02T7X2NPsRH9AXg66O7WLd/A4HegSRGD+nglgshxMXLOblo6oYV7PpmNYMTZqHXG6koz6eyohCtDWV3dZoOh0NJWRcPJJEuhOhSGo8kB5qdbNRo0GOpqXeVdvloy363+ugywl8IITpHRkoKmRs2YK+vR9Pp0On12FsYiQ5gtVgozs6msrAQlOLwli3MfuopGYkuhGgXf9n6b1IyNpCcMJO///RPADgcDoqKijq5ZV1LSsZGckvy8DF6Y7XbOFlRRL29nqkDJ2CqMrMrb5/b9jXWWjYc+hI/L1/C/ULJLTuBl95Iev5Bpg2a2Em9EEKIi5NzFHlCYjKlJcepKM9HKTsaOnz9QqmpLm12X4fDRvqu9/D1DQZgyky5JneSRLoQoktxjiQ3lzX8RDQyxJ+iUkuT7ay2HyeuCw32ZemCCWxJO4y5zEJWrklG+AshRCdJSE6mODubvD17qKmooLasrPmNNQ1UwzemBh+fhoS7anmSaSGEOBebD33FZwdSqbdZOVKc29nN6dISo4egaTpq6msJ8w+lqs6CpmnEhkaz+PIbuPWtB7DaT/uiVCn8DD4sGJnMph+20T+ij0w0KoQQncA5Mh0gNCyWhmkyHTjs1haT6AB6vRdXTPolOp0mZV1OI4l0IUSX4hxBbi6z8P6GvZRX1ba4vV4HDy2ewrSxg0g/lM/qT3cRHuIvSXQhhOgEpqws0lauxHzkCOX5+ThsNrwDAqguKfG8Q6OkeUT//vS74gpqTiXeTVlZRMbLz0iFEOeHs5zLxswvsdRXE+QTyMNTf9HZzerS0vMPUlZTjk3ZqbPW4Wv0xkvvxYaDW/l4/+dNk+iAQW+g2laLptPx1xv+SErGRnJK8mTSUSGE6ERJ45ewb8/HWCrb8ssrjWEjr2PaVQ+5lpgKs6TMyymSSBdCdCnOkeT/SdlJeVUtjuZqu5yi0+nIPNpwMpByLkII0bkyUlLYs2YNdqsVTdNQSuEXFtZsIl1nMOAfFYV/aCgzH3uMQdOmkbpiBbtWr8Y/PFxKuwghzoss01GWp/yJrOKjTOh/OX5eviydcKuUG2lFcsIMsouPsf3It5iqzFjrbG7rjToDNocN0FAoovzDiQ3rTf/wOLfJRrcfiSWn9Dggk44KIURHcybBQ0NjW0yk6w3egKJn7wQmz3jQLXmekZ7Crm9WA1LmRRLpQoguafX/dreYRDcadESE+lNeWcuOvUf5zUufsGRukoxEF0KITpSQnIzFbCZv925Ofv89DrudKpMJTa9H2e1Ntlc01EdPuPVWBk2b5nYMi9kso9KFEOds86GveGTd01TUVpLQazAPTv65jIpuo/jIfjw4+eeE+AZzvDSf3Xn7MVWZUTRco1sdNnSahkFnIMI/jFDfYPLLC5gx6EriI/u5SrokRg8hPf+glHgRQogOZirMIuWD5RQXZmGz1bW4rcNuRSkHlWUFpG1fSV7ObsrL8gH3eusXO0mkCyG6JF/vlsOT1eYgOiqY+D6R7DmYR15BmZR0EUKIThYZH8/sp57ClJXFB/feS96ePdSWlzfUQtfpwOFA0+kaUjCaht5oJCAykujERLdj+IeHy6h0IcQ5yzId5ZF1T3O8LJ/YkGieTf61JNHPUErGRt7f+ykO5WDGJRP5+PsNbiVdHEoxMKIfPl4+5JbkMbjHQFfCPD6yn2sEuoz+F0KIjpeRnkJhfiZ6gxGHo+mglsaCQntTV1NBYEhPvt3xJmgagUE9iI5JdKu3frHTdXYDhBDCk5q6pjUXdRr4ehvR6zSiwvx5+LapPLssmYVXj+bGq0ZKSRchhOgiIuPjuf6vfyU4JqZhgVJoOl1DEl0pUAodYKutpez4cfLT0932T0hOZvSiRSQky6gXIcSZyTIdZUXqv1w10UuryzHqjSwYmSxJ9DZo/PxBQ3kXh3JQXlvBRxnrm9RFD/MLITY0mv35B4kKjJQvK4QQogtJSEymR/Rgqi0l1NVWAprH7XQ6A0GBUfj6BqNDQ9P06HVG6uuqyc9L97jPxUpGpAshupzN3xyits6G0ajDanWg0zR0Og1/Xy8qLXWgwaQxA5k2dhAAT90/u5NbLIQQojHnpKO22oYJozWdjh6DGmJ2cXY2OoMBTafDOyCA/uPHN0mYR8bHy0h0IcRZcdblhoaSIpqmUW+rZ8exXTzUyr4Xu8a15KGhnnl8ZD9+mngN/9m1tuELUEdDnXSdphHhH86k+HEcOHkIq91GaXVZ5zVeCCFEE5E94hk2Yg452Wko5cDHN5jamnLXep3eSFh4X3x9gxiVtAhLlYnomES2bX2FwvxMekQPlnIup5FEuhCiy3nlnW0cPFKIr7cBh12haaDX6/D20lNepdBpGiGBvp3dTCGEEM1wTjpqq6vDy9+foF69mPnYYwBse+UV+l5xBcd27GDC0qWu2uhCCHGuskxHMVtKSYobjdlSypbDO/DSG+kT2pulE27t7OZ1eSkZG8kqPkp8RD+3euaHTEewOWz0C+uDqaqYytoqlAIfgzcbMrcyru9oquqrqa6vISVjo0woKoQQXcjutNXY7fXoDV5cMemX7Nu9luKiLHx8g+g7YBxVVSbKS/OxVJlc5VvCwuNcE41G9pD5ihqTRLoQostZumAC2ceLMZVU4WXUExUWwOXD+zJ2eByr/7eb/jHhLJmb1NnNFEII4YEpK4u8vXvR9Hp0BgP11dWU5+e7yreU5uS4/s5PT5dEuhDivGg8mjo+oh9pObuYNXgKd4xbSHLCDCk30gaNJwdNydjoet6WTriVamsN/cPjOGLOYdfxfQCYLGY0NEyWEu678nY+2r+exOghndkFIYQQpwmP7M+J43vRaQZ0Oo2b71jlmoC0qsqE6eQhInsNIjomkdQNK1zJc6mJ7pkk0oUQXUpWrom1m/Zhqa7DanNgtTkoMFcxIDaCW5Iv55bkyzu7iUIIIVqQkZJCVmoqtrq6hklGlcI3JAT/yEh2rV5Nj6FD8QkMpMfgwVIDXQhxzpy10LOLj7HnRAaDogawdMKtpOcflAT6GXJODroi9V9u5XFe2fYm/cPj2HYkjTC/EEJ8AimvraTeZiXcP5RDhdms3v0Rpqpi0vMPysSiQgjRhUye8SC+fiEAriT5hMlL2bb1FZRS2Gz12OprOZy5hcyM9QCuJLqpMEtGpp9GEulCiC5l5YdprN2Yjt2hgIbJRWdPvrTJRKJZuSZSUjNInpJAfJ/IzmiqEEIID6ITE4mIj8dSWorVYsEnOJi+SUnsXr2a/L17GyYgVYrRixYBkLpihSuhnpGSQkJyMpHxcqEuhGgbZ010L50XNdZaIv3DmTZooiRzz0Fi9BC2H4klMXoIT254me8LDjEwoi/RQT05ZMomLqQ3FmstGhAZEEZe2Un6h8cxY9CVbiVhhBBCdL7IHvHMnvcU0JAYT92wAkuVmdLiHLx9AlHKjrn4KP3ir2D02EX/v717j4+yPvP//5pDzgdymgCBgEAwgIEoWOMBWwIottV42nWVfqtSq+26davY/WkP28O6tnXXYle3ttsTxbboUk872lakEmwBDcghEkkCE2IChCSTczI5zeH+/RFmSCAJAZJMMnk/Hw8fIffc98znYiYX8Zprrk9gJrqzxhHoXAfUoX6SCukiMuaYzSa8PgOz2cSNyy7hJ9/6+zPOsecXsfGPewBYe0/uaC9RREQGUFVYSHtDA+72drpaW4lKSKCioICZOTmERUez8JZbcDmdZOXlUWS3s2fjxsC1/j9ro1EROZt3Sv/G89tf4OqLljAzMZ3qlhp8Ph9OV32wlzbubT28kw+rDrL18E663F0AVDYeZ0aiCYAFUy/m6tmfAGD53KvV/S8iMk4UFdrZtWMD0TFJRETG0dnZCiYzSbZZ5Fyzpk/XeVGhnboaBymTM7ThaC8qpIvImLI8Zy57i49hS4whfUrigLPQ/R3qp3eqn06d6yIioysrL48jO3ZQ/dFHxE2ZwuT580lMTydnzZozOs39nei9R7xo3IuInI3DWc6jr/8bJ1pq2H+8iNiIGJbOzmFK/GRtKjqMmjpaiAiLwISJLq+bsvoKLpuexSPL7u9TNFf3v4jI+JCVnccRxw6OVezD53VjGAZmk5nZc64G6DMj3V8811iXvlRIF5FzNtzF6d73V1hahbOhjeuuyhy00zxjhm1InejqXBcRuXDNJ05wcONGFg5h7IotI4O8p56iyG7HVV9PyebNTL/00n6vs2Vk9Ok+Vye6iAzmndK/8fTW/6HT00VDexOYTLi9HjJSZvHIsi+qI3qYrMm5g+SYROpdjVQ1VxMdFonVYsXt9TA7eSb2oi1kp81n6+GdgfP1dy8iMvbZJmeQd/tTFOxYT2PDUZoajpKQlE5HexOvbHyI5qYqoGeMizYc7Z8K6SJyzoa7ON37/obaaT5Uw31/IiITUUVBAR+++CImhlbs9hfInQ4HMcnJ6jIXkQvmcJbz+Bs/oLLxOGEWKwvT5jM7eSYJUfEq5A4z/6ajDmc5b5e8y9GmKi6bMpfrMj9JvauRDbs2EW4Jp97VCCYoqXHwVN439ByIiIwD/pnp+W+vo+poIe3tjbS11AIQG59K2vTsIK9wbFMhXUTO2UgUu+ubXNQ3uYDh7Rwfaue6iIgMbGZODtb29nMuiJ/ecS4icr7sRVto6WwlJjyapbM/wbdveESF2xHgcJZjL9oSmHn+w5u+zvPbX+DBpXezIvNaHM5ySmoclNQ4yJw8h053J/uOF7G+YBNP3vjYkO9XRESCyz/m5ejHezAMA4CO9ia2b3uepOSZGucyAHOwFyAi44+/OD1cM8f99/Pin/bw0PdfwVHpHJb7FRGR4TFp6lSWPfzwWce6iIiMlLys68iaOo+4yBgunZ6lYuwIsRdtYeOeV7EXbRnwnHmTM7hl0Q08d/sTXD3rE1hMZy8rDOV+RURk9PjHvFx+5efIXnIb2UtuY3LafOpqHBQV2oO9vDFLHekiMma4PT5Ky2ux5xepi1xEZAw5lxnpIiLDweEsZ33BJpo6WgLjW57K+0agq1lGRl7WddS7Gql3NeJwlvP89hfY+fEHAMxMms5j9u/jqCvnnit6xun456mf7Tnx367nTkQkuJwnC+X+TURvvO3Jfm+T/qmQLiJjwppbcwJ/1jxzEZGxpfeM9Ky8PIrsdrJUVBeREeAvoO89doDSmjK6vd2YTCaaOpqZk3KRRoOMsAzbLJJjEtm451WSYxJ5cOndADy49G7sRVtw1JWTkTIrUBD3z1Mfyv0O5TwRERlZRYV29ry/EegZ79K7qK4NRs9OhXQRGRMyZth48qs3BnsZIiLSj94z0ovsdvZs7PnlW/PPRWS42Yu2sGmfHbfXQ0JUPNHhUVS3OjlSV0lBxV4AFWRHiH+OeXbafOC2wJsWKzKvBXo60gG9mSEiMo75u839RXR/Ub13Af30rnU5RYV0ERERERnUpKlTmfvww5jN5sCGo+e68aiIyNm8U/o3tpT8lasuWoLT1UBVczW5c68hOSaR7LT5FFYVazTICPLPMYfbWJv7AA5nOd9886k+43X0JoaIyPjWu+u8d1G9t4EK7KJCuoiIiIicA1tGhjrRRWREPL31f9h77ABT41P50S3fDhTO/d3P/s5oGRmnzzH3fzqg09NFpDWC5JhEFdJFZEIItY7sgeIZaJTLQAV2USFdRERERERExoDZKTPYf7wIZ1s9Ww/v5MkbHwv2kiaU0+eY+zce9Xek69MAIjJRhFpH9rnGo1npA1MhXUREREREREaVf1NRgDU5d5Bhm8Ujy+7nSH0lpTVlQV6d+Oelr8m5A+jpTu99XHPSRSSUhVpHdqjFE0wqpEvIc1Q6secXkZebRcYM21mPi4iIiIjIyHpm2y959cM/EWEJD4wMybDN4rnbnwgUaiV4Ts1L79HfnzXmRURCVah1ZA8lnlAbZzNSVEiXkGfPL2LjH/cAsPae3LMel+FTdrSOd7Z/yIqllzJ3ZmqwlyMiIiIiQdK7A31e6hzeKt6K1+clIbbvyJDTx4tIcJw+L73e1Ui9q5Hlc68GbtMbHSIi49jpRXNnjQP7K49RV+MAQmOczUhRIV1CXl5uVp+vZzsuw+eNbUXs3l9Gu8fK2nuWB3s5IiIiIhIEDmc5j9m/z75jRVjMZmLCo2ntchEXEcOPbvm2RoSMQae/oZEck8jGPa9qw1ERkRBw+sz0okI7dTUOUiZnaPzLWaiQLiEvY4at347zgY7L8LlpWRbRVg8rlurNChEREZGJyl60BUddOZmT5zA7eSb7jx0g3BLGDfNzWZF5bbCXJ4Pwf5KgqaOFVfNy1YkuIhICTp+Z3vurxroMzjycd3bgwAG+/OUvs2rVKv7pn/6JgwcPnnHO/v37mT179nA+rMiQOCqdrNuQj6PSGeylTBhz0lO4ZcUi5qSnBHspcg6Uy0VExj/lchkrHM5y6l2NLJ2dw+LpCwE40eJkYdp8Hll2f5BXJ2ezvmATG3Zt4v8OvAWgTw+MMuVyERkJ/pnp/qL56d/LwIatkL5jxw4+8YlPsH37dhITE3nrrbe49NJLeeKJJ/qc19XVRUVFxXA9rMiQ+Wei2/OLgr0UkTFLuVxEZPxTLpexwOEsZ13+z/m3t57hhd1/4GB1KZtL8jlSX4HFbGbx9IUqyo4jRrAXMAEpl4uIjD3DNtrlG9/4Brfddhu///3vMZlM+Hw+nnvuOb7+9a9z8OBBfvOb3xARETFcDydyzjQTPThOOJvZuPkgebkLyZhhC/Zy5CyUy0VExj/lcgk2/0z0kprDdHq66PZ042yr5+aFN7B87tUUVhVrRMgY5XCWYy/aQl7WdWTYZrEm547Abb3/LCNPuVxEZOwZto70wsJCvvCFL2AymXru2Gzmq1/9Kvn5+Wzbto3ly5dTV1c3XA8ncs78M9FVzB1dBR9W8OKf9uqTAOOEcrmIyPinXC7B1LuIHh0ejdUcRnxkPIZhkByTyIrMa1mb+4C60ccoe9EWNu55FXvRFoBAMT05JjHIK5t4lMtFRMaeYSukR0dH09bWdsbxnJwc3nvvPZqamrjyyispLi4erocc1E9+8hMuuugiIiMjycnJYdeuXYOe/4c//IF58+YRGRnJwoUL+dOf/jQq6xQJdTmLZnLXZxbrkwDjhHK5iMj4p1wuweIvon9YdZBOTxfpCdP4u0s/y8/u+AH3XbVaXejjQF7Wdaxecluf5+r04rqMDuVyEZGxZ9gK6YsXL+b//u//+r3toosuYufOncyYMYMvfvGLw/WQA/rf//1f1q5dy3e+8x327t1LdnY2q1atora2tt/zd+7cyV133cV9993Hvn37uOWWW7jlllsoKlIHrciFmmqbxMOfX6ZPAowTyuUiIuOfcrkEi71oC466csBEa5eL3Uf3qQs9BPRXXJeRp1wuIiPJWeMg/+11OGscwV7KuDJshfS7776bQ4cO0djY2O/tkyZNYvPmzdx7773MmDFjuB62X+vWreP+++9nzZo1LFiwgJ/97GdER0fz61//ut/z/+u//osbbriBf/mXf2H+/Pk88cQTLF68mP/+7/8e0XWKiIw1yuUiIuOfcrkES17WddxzxR18e9XDLElfRF7WKhVfxxl1n48dyuUiMpKKCu3seX8jRYX2YC9lXBm2zUbvuOMO7rhj8M1HwsLC+OUvfzlcD9mv7u5u9uzZw9e//vXAMbPZzMqVK3nvvff6vea9995j7dq1fY6tWrWK119/vd/zu7q66OrqCnzf0tICgM/nw+fzXWAEQ+Pz+TAMY9QeLxhCPcZQjw8U40gxm4ftPdAzTKRcDsHP5/oZCQ2hHmOoxwfBi3Gk8rlyuXL5cBsoxq2HtvOzHb/j5qzrcboauClrJQ9/qqc79nOX39bn+rFuIj+Pvd10yUower76z7Mf2MKLe18DAx5eNvLdz+dLuXxkTJRc7n+sUM4DoR4fKMZzccmimzCMnq9j7e9rLNdZzrmQ3tDQwLPPPstbb71FRUUFAImJicydO5fLL7+cZcuWsXTp0sCGGKOtrq4Or9fL5MmT+xyfPHkyJSUl/V5TXV3d7/nV1dX9nv+DH/yA733ve2ccdzqddHZ2nufKz43P56O5uRnDMEa0qBZMoR5jqMcHinGkTJky5YLvQ7m8R7DzuX5GQkOoxxjq8UHwYrzQfK5c3kO5fOT1jvFAVTGvHXiLKXGpVDYdp7mxkd9vf4WYsCiiPeHELYoO9nLPy0R7HgeKMY5oVi/IAx+BsR2503KI9oSTM+2yAUd5jAXK5SNjouRyCP08EOrxgWI8N3EsyF6ND8Zcbh/LdZZzKqQfPnyYT37yk9TW1mIYRuB4TU0NpaWlvPnmm0BPgrzvvvtYu3YtiYmht7v317/+9T7vrra0tJCeno7NZiM+Pn5U1uDz+TCZTNhstpBODqEcY6jHB4pxrFIuPyXY+Xw8vn7O1WAx1pWVUfTGG2TddBMpc+YEaYUXLtSfx1CPD8ZnjMrlpyiXj7zeMf7iz99lx5HdhFksJETG43Q1YmAwI3EaKy79JKkpqcFe7nmZaM/jucSYmppK9tyFI7iy4TEen0Pl8lOCncthfL6GzkWoxweKcTTV1ZZRVPgGWdk3kZI6vP8vN1Zi7M85FdK/9rWvUVNTw5e+9CW++tWvkpqaSnt7Ow899BB2u527776brVu3cvToUb7//e/z05/+lJ/+9Kf8/d///Uit/wwpKSlYLBZqamr6HK+pqRnw3YUpU6ac0/kRERFERESccdxsNo/qE2wymUb9MUdbqMcY6vGBYhyLlMtPGQv5fLy9fs7HQDF+9MYb7N24EROQe9rHf8ebUH8eQz0+GH8xKpefolw+Oqpbanmx+A2uumgxjrpyosOjaO5oYd7kOUSGR/K13C8xN3V2sJd5QSbC8xjqMY63+JTLTxkLuRzG32voXIV6fKAYR8tHH77B3oKNmEyQe/3w/7/cWIixP+e0mr/+9a9cfvnl/PSnP2XevHkkJSUxffp0EhISAFi/fj0VFRV88MEH3H///bS1tXHnnXfy3HPPjcTa+xUeHs6SJUt45513Asd8Ph/vvPMOV111Vb/XXHXVVX3OB9iyZcuA54uIjGfK5TJWZOXlsWT1arLy8oK9FJFxR7lcRovDWc6//vE/+Z+dv2P9+y/y0j47Bga5c6/hvqtW84u7/pM/f/l3rMi8NthLlRHgcJazLv/nOJzlwV5KSFIuFxkaZ42D/LfX4axxBHspclJWdh5LrlxNVvbE+n+5cyqk+3w+MjIyznre4sWL+dnPfkZhYSGZmZmsXbuW3bt3n/ciz9XatWv5xS9+wYYNGyguLuYf//EfcblcrFmzBujZ/br3Rhlf/epXeeutt/jRj35ESUkJ3/3ud/nggw/4yle+MmprFhEZLcrlMlbYMjLIXbsW2xBejyLSl3K5jJb1BZv47e4/UF5fSaeni5bOVjJSZrEm5w7W5j5Ahm1WsJcoI8hetIWNe17FXrQl2EsJScrlIkNTVGhnz/sbKSq0B3spcpJtcga516/FNnli/b/cOY12yc7O5oMPPhjy+ZmZmbz99tvMnz+fH/7wh7zyyivnvMDz8Q//8A84nU6+/e1vU11dzaWXXspbb70V2OyisrKyz0cDrr76ajZu3Mi3vvUtvvGNbzB37lxef/11srKyRmW9IiKjSblcRGT8Uy6XkfZO6d94fvsLYBh4fD4Mw8Dr87Fo2iU8lfcNFdAniLys6/p8leGlXC4yNP6u54nW/Sxjj8novaPFWWzatIk777yTJ554gm9+85uB42vWrOGFF17A6/X2e91tt93G9u3bx9wusMOlpaWFSZMm0dzcPKqbYNTW1pKamjrm5gUNl1CPMdTjA8U4VimXD2y08/l4fP2cK8U4/oV6fDA+Y1QuH5hy+YVzOMu5a8M/cbSpipToRFo6W/nUtBxmpM1gzZV3hGQRPRSfx9NdaIwOZzn2oi3kZV03Jl8D4/E5VC4fmOoswy/U4wPFGCrGcozntJo77riDz33uc3z729/mC1/4AtXV1UO6rqmpiY6OjvNaoIiIDC/lchGR8U+5XEbSM9t+QVVzNSkxiVycOgeTycSc5Jk88dl/GZMFVBkdGvMy/JTLRUTGl3Ma7QKwYcMGkpOTefbZZ9m4cSM33XQTlZWV/Z7b1dXFM888w7Zt27SphIjIGKJcLiIy/imXy0hwOMvZcWQ3bp+HDncnYDp5i2mwy2QC0JiXkaFcLsPNWeOgqNBOVnbehJtfLTLSzrmQbjab+fGPf8ztt9/Ov/7rv/aZyRUXF8fUqVOJjY2lu7ubI0eO0NXVhcVi4YknnhjWhYuIyPlTLhcRGf+Uy2UkrC/YRJ2rAYDWLheFVUWEmcOIiYgO8spktPUe5QKM6bEu45lyuQw3/8acALnXrw3yamQkNDed4GDhRhZeqjdLRts5F9L9rr32WrZt24bD4eDll19m27Zt7N27F4fDcerOrVZWrlzJd7/7Xb1bKiIyBimXi4iMf8rlMtzcXg8AZpMZqzmMhVPnsfLipUFelYw2/ygXgHpXI5v22al3NfLkjY8FeWWhSblchos25gx9FeUFfLj7RUwmvVky2s67kO6XkZHB448/zuOPPw5Ac3MzjY2NmEwmpk2bhtV6wQ8hIc5R6cSeX0RebhYZM2zBXs6EpedhYlMuFxEZ/5TL5UI5nOUAxEfG0tzZyvRJU1l9+a3cdMlK4nzqSJ8o/J3o2WnzgdvIy7qO9QWbgr2sCUO5XC6UbXKGiqshbuasHKzmdr1ZEgTDnoEnTZrEpEmThvtuJYTZ84vY+Mc91De5SE6IUSE3SPzPA8Dae3KDvBoJNuVyEZHxT7lchsrhLGd9wSb2HjtAVXM1189bRnVLLbcsXIXz5JgXmThOdaLfxtrcBwBYk3MHyTGJmo8eBMrlInK6SQlTmXvxw5jN5mG5P83VHzq9lSlBl5ebBUB9k0uF3CDyPw/+ryIiIiIS+hzOch6zf599x4oAuGx6Fo8s+yIZtlmsy/95T0HVgNUL8gLna1Z2aOtvU9EM26xAUV1EREKL5uoPnQrpEnQZM2ysvScXR6WT5IQYsjPTWLchn5uWXUJcRLBXN3H4nwcRERERCX29u9CP1FUSHRZJVFgkDe1NVDQcI8M2K1BIvemSleDrua733GwVVkOTiuYiIhOLf0RM2vRs8t9ep870QQzPZwBEhoG/kFtYWsXGP+7hjW1FwV6SiIiIiEhIshdtYdM+O6W1ZcRHxtLu7uRYSzUHqw/x/PYXgFMF1TkpFwWuy8u6jtVLbtOIDxERkRDhn6tfdayQPe9vpKjQHuwljVnqSJeg8G9smZ2ZRmFpVZ+56P7RIjctu4RA64uIiIiIiFyQ3mNZ8rKuo97VyNHG4zjqPmZqfCqT42x0eDp5cOndA96HupVFRERCk78zXZuYDkyFdAkK/8aWO/YdwVFZx459R3hqbR4ZM2yBznSfz0dtbW2wlyoiIiIiEhL8Y1nK6j7m4IlDdHndtHS0UOuqJ8Iazu2XflZFchERkQnK35kuA9NoFxl1jkon9U0uVl0zjwfvXErGjBQclXXY8zXKRURERERkpPjHshypq+SjmkM46spxuhpIikogN+Nq6l2NOJzlwV6miIiIyJikQrqMOnt+EZt3lABQWFrFg3cu5Z6brwiMdJHR5ah0sm5DPo5KZ7CXIiIhqK6sjA9ff526srJgL0VEZMLzj2X52vIvccnki4kKi8TAYFbyDNITp7Fpn531BZuCvUwRERGRMUmjXWTU+Qvm9U0uNv5xDwBr78kN5pImNP+YHdDzICLDr+iNNyjbvRtrezvL1+pjgiIiweSfkZ6dNp+bFl6P4fPx0j47q5fcQkmt3vAUERERGYwK6TLq/DPQHZVOkhNi1IkeZP6/fz0PIjLcnA4Hrvp6pi9eTNZNNwV7OSIiE55/RvqW0hQcznKmJ6TR1NHM+xV7SYiaxB2X5bEm545gL1NERERkTNJoFxl1/lEi0NMBnTHDFuQVTWz+Nzb0PIjIcCuy2yl9+20i4uJImTMn2MsREZnw/DPSZyfPPHnEwO318FZxPq9/+GeSYxLJsM0K6hpFROT8OWsc5L+9DmeNI9hLEQlJ6kiXUadRIiIiE0NWXh4GMC0nJ9hLERERemak52Vdx/qCTVx10RIcdR9jYNDa5SIuIpa8rOuCvUQRETkPzhoHRYV2XG31lBRtBiD3eo1VFBluKqTLiHFUOln/WgEAa27NCXQ8a5SIiMj45XQ4KLLbycrLw5aRMei5towMlj38MLW1taO0OhEROZtntv2C/zuwGZPJRJenGwCTycSitAXqRhcRGaP8hfKs7Dxsk8/8Hbyo0M6uHRuYlJDGvKxVZGXnBWGVEsrO9hqcKFRIl2HhqHRizy8KFMfXv1bA3uKjlJY7sZhNJCfEBOai+8/TKBERkfHDX0B31ddTsvlkl4s2DxURGTccznLWF2zi3cPv0e11B45HhUViMVtIT0wL4upERGQwRYV29ry/EYCs7DwKdqyno72JqOgEcq5ZQ1Z2HkccO6ircZB5yXUTutApI6P3a3Aif9pBhXQZFr3HtQBsemsfXp9B5iwbi+enk52ZxroN+dQ3udi8owTQWBcRkfGkyG5nz8aNzFu1iiWrV5OWnU3+unVD6kwXEZHgsxdt4cW9r9Pp7iQqLJKEyHiSYhL5Qs4/4HQ1aKyLiMgYlpWdh6utHldbPdu2PEPRfjuG4SMsPJqY2GRyr19L3u1PBTqGRYab/3U10V9fKqTLsDh9XEt9k4um1g4S4qJYc2tOoNC+6pp5rP7sEo11EREZZ7Lyen5hSsvOpqqwkMNbt6ozXURkHMlOm098RCwd3Z10ebqpa2/kxqzruGrWEuxFW4K9PBERGYRtcgYd7U0U7bcTO2kyhmGQmHwR02dciqutntKD71B1rHDCj92QkWObnDGhO9H9VEiX89Z7Bvq8Walsea+EHXvLaO9y87V7l1NYWsWvXnmPt3eW8M+f+2SggK6RLiIiY9vpc9BL33mH7c8/z8JbbmH7889TU1JC0owZzMzJwVVfj9PhUFe6iMgY5XCWYy/aQlldBU0dLYRbw+jydIFhAD2d6hv3vArA2twHgrlUEREZRL3zCB5PNy2NJzBbLIDBkcM7cHs6KPnobbwn971QsVNk5KiQLufNnl/Eprf2AZAQH0XliUYMA0wmePo3W5k9PYXObg91jS5e33qAV358X5BXLCIig3E6HBSsX0/5zp00VlTgqq/nxiefZPvzz/Pxzp00Vlbi7e4mPDqa5qoqwqKjqSgoICY5WV3pIiJjhL9wnpd1HRm2Wawv2MTvPngFt9eDx+cBwGqysChtAWty7gCg3tVIvasRh7NcG46KiIwhzhpHYB66NTwai8WK1+vGMLzU1ZYBBhZLGO7uDqamXTLhx26IjDQV0uW85eVmUd/kAno60p/61V+ob2pnUlwknV0e/vy3g/h8BlNT44mJCueb//Umy3PmUlhapc50EZExqMhuZ9+mTXS7XGAYNB49yh/+6Z9oOXGC1HnzWHDjjRz+y1+IsdlITE8ndd48Drz+OmnZ2cFeuoiInGQv2sKGXZvYcWQ3Dy69m53lu+lwd/Y5Z1JUHIvTF1LRcIzCqmIANpfkkxyTqK50EZExpGDHenbvfAGv14theAATAIbhO3mGCa/XgwWYPnOxxrqIjDAV0uW8Zcyw8eRXb8RR6cSeX8TiBelsLTiM1WKhtr6VmWlJgEF3t5d33j+EyWRib/ExnA1tgDYbFREZa7Ly8qgrK+PQ1q10NDdzbO9eOhob8brdhMfEsGv9erra2rA4HCSmp7N340ZqSks5vHUrmStWBHv5IiIC5GVdx5bSv7LvWBH//vazHHaW95RdTCbMhgkfBonRCWwuyaekxkFF41FWzctl9ZLbtOGoiMgY4axxUFRop7HhKB5PN2CcvMXoc15Ph7qHiMg45s5bPurrFJloVEiXC+bfSDQizILb46OptYO5M200trTT0NROt8eLYRhYLGZmT0/muqsytdmoiMgYZMvIIGXOHD764x8xvF7aamuxRkWBx0O3y0W3y0V0UhIXr1hB+c6d1B7qeZNURESCzz/SJTttfuBYa2cbHp8Xs8mC1/DixcBiMnPZ9IXMSZlJdtp8CquKA2NgRERkbCgqtLNrxwZ6OtCNAc/zet0AdHS0sPXtpzlcspWca9aoM11khKiQLufN34menZlG2dE6/m/rAQzDICLMwglnCw3N7ZjNYPh65qanJsVy28pFrLgyE0elk3Ub8jXiRURkjPBvMBpjsxERG4u7vR0AT0dHzwkmExgGlogIKnbtorWmBpPJxNSFC8lZsyaIKxcRETg10sWEiYb2JuIjYmnubMEX+Ph/j9jwGG5bdAMrMq8FCHwVEZGxIys7j6L9b1BbcwiTyXpyrMtATERExHDi6AGcJ0qJiU3WhqMiI0SFdDln/gJ6fZOLN9/9iIwZKVQ7W+h2ewFoa+/GoGe3aN/J39ujI8MwDCgsrWLFlZmBLnY4c8SL//5vWnYJcRGjF5eIyERWsH49H/z+95gtFnyevr+om8xmpl12Gcf376f1xAkAYlJTWXTzzeSsWYMtQx0vIiLBlpd1HTuO7GZ35X66vN00djTjM3q6GMNMVtwnizDt7g6e3/4CM5OmqwtdRGQMa2k+ge9kx/ngDGLjUmlpPkHqlExtOCoyglRIl3PmL4KvumYeaamT2Fd8DJ/v1EeNTv/QkdlswueDOenJ1De5cFQ6A6Nd+hvxcqrIbrB61YIRjERERHrzdnXh9noJj4kBsznwbqjh89HR3Hzq3VEgffFibnzyyWAtVURETlPRcIyG9iYirOF4fF4uSkqnob2R+vYmvHgJM1uJDo/C4/XwUXUp9qIt2lhURGQMctY4sL/yGO2uhkHO8o9X7KnApKTOYe68ZSO9NJEJT4V0OWe9i+BNrR0cOFSFx+sb8Hyfz6Cjy82BwyeoPNFEckIMa+/JHXCzUf/937TsEmDg+xURkfPnH+WSlZeHLSODnDVr6Ghq4tA779De2Nhz0slxLgD1DkcQVysiIv3xz0XPy7qOp/P/h4+qSzFjIswaxqXTs9h2eCcACZGTWJF5LTuO7OJESy1WS1ifWeoiIjI2+Ivox48WDnqe2WzFMIzAyJejFXtpajhKa0stNdUl5N3+lOaki4wAFdJlyByVTta/VgDA8py5rH+tgH3Fx/B4fX060gdy3VUXkzEj9awbjWbMsLH2nlx8Ph+1tbXDsnYREemrYP169m3ahKu+nhuffJKGigqqDx6k2+UKFM8H0+Z0jsIqRURkMPaiLWzc8yoAs5Nn8uHxg8RFxrLi4mtp7WzD6aoHoKWzlYPVh2hobyI2IhqzyaTxLiIiY1BRoZ2aqhK87q5Bz/P5PJjNZgwDLNYIXK1OOtubiE+YSl2Ng6JCu+aki4wAc7AXIOOHPb+ITW/tY9Nb+3j+pe1semsfHx+vH1IRHaD8eANr78nts7mof9NRR6UKMiIiwfSnf/1Xqj/6CE/X4L+0A0ROmkTK7Nk4HQ5K33mHX91+O6XvvDMKqxQRkd7ysq5j1bxcyuoqAFg0bQEWk4V9x4r4y6G/Bs7zGT4OVh/C6/NyzawriAqLZM/RD1lfsClYSxcRkX5kZeeRlDwDn8/LqfEtZzKbrYH/pkzLIj5hKnPnL+em23/IJdk34mqrx1mjT5SKDDd1pMuQ5eVmUd/koqm1A4ApKfG8vaOE5rbOQa8zmQADbImxZ9w22KajIiIycnLWrCEmOZm07Gzy162jvWGwGYy9mEx4PR4++uMfiUpIoKakhI939owOyFyxYgRXLCIifr1HuiTHJLJpnx2vz8fMxGmEWayU1X2M0WvnIhMmDHzYYpOJi4ylztVwxr5GIiISXM6TneTJttkcrdjLGTvQmcxg9Iy/9fncge2LXK1OurtcJCalk7lgBVXHCtnz/kZiYpPVlS4yzFRIlyHLmGHjya/eyLoN+Wz84x5yFs3E4/We9TrD6Cmmp09JPOO2wTYdFRGRkWPLyCB37Vry161j14YNRCclAeD1eOhqaRnwOpPJhNvlwmQ2c2zvXhavXg3A0gcfHJV1i4hI35EueVnXUe9qZGf5bg7VHgEgKiySdncHFpMFr+HFMEF8RDyP5n6Jktoywq3hZKbOYU3OHcEMQ0REeinYsZ59uzeRnDqn/xMM/x5yJ7sVT2prdWL4vHS0NwE9Xe29v4rI8NFoFzln2ZlpzExLZF/xMVwd7kE+bHRKYnwUy3PmjvjaRETk3GTl5TEpLY3m48cxW614u7sHPtlsJjIhAZPZTFRCAs1VVbicTu575RV1o4uIjKK8rOtYveQ28rKuI8M2izU5d1DbWo/b58Ht89Du7vkEqXGy6BJuCaPL08WP8v+HpKhJXDYti6/lfknz0UVExqC66sOc0Y3eizUsEqs1ksioSYSFxYDhw2Tqqczkv70OgNzr12qzUZERoEK6nBNHpZPnX9pOyZEaWto6TnsfdGAR4WEUlladcdw/2sWeXzTsaxURkbOzZWQwffFiPN3dtNXW4ukcZFyXz0dnUxNmq5WLV6zginvuIStv4na6OB0O8teto66sLNhLEZEJqN7VyDPbfsk333yKf3vrmcDGor0ZQGx4DJdMuZguTzdVLTW8tM9OReNRCquKR3/RMuY5nOWsy/85Dmd5sJciMuHMnbccW2oGUdFnfpq/t9j4VKJjk/B63ZjMJlJSM5g6redT/nve30hRoX00lisyIWm0i5yT9a8VsK/4ONGRYbS1d2G1mnF7fINeE2Y1s3Tx7H7Ht2i0i4jI6HI6HBTZ7WTl5WHL6OlSyVmzhuLNm2mqrDzr9SazmZS5c1n2yCOB6yeqIrudPRs3YgALTo64EREZCb1nomfYZmEv2sKmfXY6PV1EWiNwez39Xmdg4Opup/TkyJdJkXHMnzyX9MQ08rKuG80QZJzoPTZobe4DQV6NyMRSdawQZ62Dzo6BxywCNDccA0xExyQxd96niIpOoKRoM9NnLmbJlas10kVkBKmQLufkaHUjro4uOjq7sVjMpCTGUF3XijFIW7rb4+Pd3Q4euXvZGbdlzLBpk1ERkVHkL/4C5K49tflQd1vboNdZwsIwAJ/HM/j4lwnE341/yU03MfhbyiIiF8Zf3Kx3NZIck4gtJokM2yxsMcmkJ6bx6/df7HN+bFgMAG1uFwYGbd0uACwmCwdrSrl0+iUa6yJ9+N+syU6bD9ymN1pEgiBtejaxcal0dbkwfP2/QQpgGF5MJjMedzspqXPIys4jJjaZrOw8jXMRGWEa7SLn5Gh1E4YBXl9P5dzZ4Bq0iO5X19g24PgWR6WTdRvycVQ6h3OpIiLSj6y8PJasXt1nJEvB+vW0NzUNep3X7e75g2HQ+PHHFNntgdEmTocjcF5/x0KVf8PWlDkDbAglIjJM/DPRATbueZXXD2zG2VZHemIayTGJxEfG9Tm/ze2ize0KfG8xWwCIi4hhZmL6yWKpyCn+N2sKq4pZm/uA3mgRGWXOGgfbtz1PQ135oEV0gKjoJFLT5jN/4afVfS4yytSRLudkwZwpFJdVg8nElOQ4KqubhnRdRHgYeblZOCqd2POLyMvNImOGDTg1Jx1Qd7qIyAjzF3/h1JiXyj17wHf2nmqf2018Whqzr7mGrLy8Pt3t/u9d9fWUbN4M9O14FxGR85dhm8Xa3AdwOMtJjkkkO20+hVXF1Lsa2bjnVeIiYmnsaO732piwKO5afAvvVezF2VpHSc1hCquKWZF57ShHIWOZvwNdnegiwVFUaKeuxoHFGoGve/BCeldXK23NFpi6gFc2PoTb3Um7qwHo2WTUWeOgYMd6AHKuWaMudZFhpEK6nJNH7l7Gjn1HqKpt4VhN/7+s9+dTn5hDxgwb6zbkn1E015x0EZHgKFi/nn2bNuEZ4qgWS3g4l3z2s+SsWYMtI4OsvDxc9fW46uspWL+eks2bmbdq1Rkd7yIiMjz8BXWAFZnXBgrrL+9/s9/zrSYL8VFxpMQlkxSdQHHNYeIiYslOm99nlEdhVXFg/rpMTL1fWyIy+vyd5a0tdRRs/+Wg5/q8bsIjoql3HuFoxV4sljBmzskJ3EdRoZ19uzcBEBObTO71am4RGS4qpMuQOSqdrH+tgOjIcEwm8J2c6WIxg3eQRkazCZwNrXzzv95kec5coG/RXHPSRUSCKzoxkZaOjsFPMpmYccUVlGzeTExyMrlr12LLyCAmOZk9Gzf2KaBP9E1IRURGwzulf+P57S9wy8JVmAY4Jzo8mqToRPYfKyImPJopcTbcXjeFVcUUVhWzcc+r7DiSTkXjUUCbS4qIjDZnjYOiQjtZ2XlkZefx65/+3VmviYiI45Mr/plJCWnkb36aZNtsll33SKDzPCs7D1dbPR3tTbja6nHWONSVLjJMVEiXIbPnF7HprX14fT5mTEmgpqGNyUmxdHu8nHC2Dnidz4A9B49TfKSWkvIanlqbFxjrIiIiwZOzZg0xyckc27+flqqqwU82DOqOHGHS1KmkZWcHDvs7z1VAFxEZXU/n/w97jx5g79EP6fa6MQGnb13U0tWKo64cR105kdYI7risJ2fXuxpZPvdq4LY+HekiIjK6igrt7Hl/I662eo5V7KWtpRb6zeg9YuNScbs7yX/7R6RNX0SybTZR0Qk01FdQVGgnbXo2VccKyblmTeC+1ZUuMnxUSJcBnT7PPC83i/qmnk2Ldu4/QsWJJipONGEe4pa1YVYzjso67PlFfTrQT38c//c3LbuEuIiRiExEZOLyz0X3F75jbDbK3n13SNe2VlXhqqnhw1dfJXPFCqDvzHURERk9s5Nnsv/4R3R5ujGbzURaI+jwdAFgwUxUeBRJ0QnMSJzG0cbjtHa7KK4+TGNHMxWNxwB48sbHADQvXUQkSPzd40cr91JTVYzJZMJkMmEY/RTSTRZ8Pg9eTxctTa20tdRgNluxhkVSU11CY10FR1J20FhXEbjv3l9F5MKpkC4DOn0T0IwZNtbcmoM9v4j0KYkcLKvBMIa0Px0Arg43118z/4xZ6Kc/zqnvDVavWjCcIYmITHi9NwjNXbuWvz77LO6zjXUBTBYLZosFb3c31QcPjvQyRUTkLG5bdANH6iuwxSQTFxlLwcf7qGzqKZB78dHp6aK2rQ5bXDInWmvp9rrZUb4bq9mCeaidMCIiMqL8I1ecJ0qJiUvB4+4kIjKeemcZp3elmzBhDYvEMAyiY5JJS19EZGRc4PbJU+Yxd95yqo4VkpWdh21yhjrRRYaZCukyoP42AV3/WgGb3trHVZdexII5U6isaqC1vZuIMAtdbu+g9+fx+jhyrD7QfQ49RfTszLR+H++mZZcAQ6zSi4jIkPQexQJw2Z13kv+f/4lxtndFT3bFmC0WpizoeZPz9O52EREZPYVVxTjb6pidPJNdFfuobq3pc7vH58Hj8xBljSQhMp46VwNmk5mLbbO5evYnmJc6h9t/dT8PLr1bHekiIkFmMlvweT20tTppa60lPCKW7q62U7ebzBiGF8MwSLHNJnfV18hc0PMJ0Tdf/Sb7dm/isk/cQeaCFYHjIjL8Qq4VoaGhgc997nPEx8eTkJDAfffdR1tb26DnP/TQQ2RmZhIVFcWMGTP453/+Z5qbm0dx1WNP73ErAOs25OOodNLU2kFHl5u9B49yrLqJ1vZugEGL6BaziTCrhbTUeGyJMfz0pe0888I2HltnZ8P/7aKwtCrQ8Q6nNh+dk54y8oGKyJikXD5y/KNY/IXvjoaG/j86epI5LKzn46U+H9aICBbddhvLHnkEONXdXmS3j8raRWR8US4fWXlZ17F6yW3sP1ZEReMxur3uM84xY6bD3YnFYmH+lLlkTp7DgqkXsybnDl4/sJmdH3/A89tfCMLqZSxxOMtZl/9zHM7yYC9FxiDl8pGXc80alub+I7mrHiUiMhaA7q42zJawwDmG4QMMWpurqa0uZfu253HWOIK0YpGJK+Q60j/3uc9x4sQJtmzZgtvtZs2aNTzwwANsPPkx9tNVVVVRVVXF008/zYIFC6ioqODLX/4yVVVVvPzyy6O8+rGj97gVoM+fPV4fzkYXliF+JNTrM4iOsvLoPbn8+rUC2rvcHCyrpqG5nYwZKWeMehERUS4fZb0L6WZzn5ldPrc7cNzrdhOVkBAowp/e3S4i0pty+cjKsM0iL+s6fr7zdwOeE2a1ctm0S4gOj2JKfCpvl2zjWNMJ5qRcxINL7wYIfJWJy160hY17XgVgbe4DQV6NjDXK5SPLWeOgqNAeGMUyKSGNP77+rzQ3HMVzct+L3kwmE253F8crCwPXdbQ3kZA4nY72Jpw1jsC4GBEZfiFVSC8uLuatt95i9+7dXH755QA899xzfOYzn+Hpp58mLS3tjGuysrJ45ZVXAt/PmTOHJ598kv/3//4fHo8HqzWk/oqGrL+xLnm5Wax/rQCMvjUXv4H3lYbOLjcb/7SXiqpGzCYTC+ZMYU56SmCD0aE4fVNSEQlNyuWja+7y5RRs2ICvu+cTRgNufOHz4enqonjzZgBy1qzRRqMiMiDl8tFhL9pCp/vMQotfl6eb3+99DY/Xw9T4KUxPSGPBlIvJy7oOgGtmf4KZSdNHa7kyRvlfD/6vIn7K5SOvqNDOnvd73pTIys6j6lghttQM6mvL+j3fMHwYhg+vx0xMrA37K49xvGIfHk83dbUOoqITuPG2J0czBJEJJaQy2HvvvUdCQkIgwQOsXLkSs9lMQUEBt95665Dup7m5mfj4+AETfFdXF11dp35hbWlpAcDn8+Eb6s6bF8jn82EYxog93uzpyTz8+U8Fvvf/+d5bruCN/APUNrRhGD4iwy14PF58Jyvopn7uK8xqwTAMWtrauXbxLI7WNGKiZwa6YRh887/eOHnfOX3GuZweoz3/AC/+aS9g8PDnl41A1KNrpJ/DsUAxjoxQ3yBstHI5BD+fj4WfkUNbt/a8CTrE11VzVRW7f/c7DOCzTzxx1vPHQowjLdRjDPX4IHgxhnI+Vy4fWWV1H/NG0V9YlDaPGy9ZwdbSHTR3toFh4D1tj6Eud88bpceaqrCaLDS3N3PrwlV8WFXCi3tfAwMeXvbFQR9PeSA0nB6j/3V0U9ZKHv7UFwPnjFfK5cNvIuVy/2ON9mvokkU3YRg9Xw/st7N752/p7GwFk4n+Kyw9TCYre3a9SEvjCWxT5+Hp7qSxvhKDgX+OJ2KeC0WKcWQMNZeHVCG9urqa1NTUPsesVitJSUlUV1cP6T7q6up44okneOCBgT/S9oMf/IDvfe97Zxx3Op10dnae26LPk8/no7m5GcMwRvQf7hPOZgo+rCBn0Uym2iYRFwH/+ch1/OrVAhqa2od0HyZTTwd7mBXcXW3EWH1UVB7nR7/6I91uL1W1zYRZLLyTEkHcikWB606PMXfxNKKtHnIWTaO2tnakQh41o/UcBpNiHBlTpkwZlccJltHK5RD8fB7Mn5HmEyeoKCigzeNh0snNQzFOfuTIn7j7ERYT09O1npAwpFysPDD+hXp8ELwYQzmfK5ePrHc+/Cu7y/YQ7QlnZuQ05sTOICElnob2Jrq83QNeZ8KEgcHP3/kt9+XcyX2L/oGcaZedNZ8rD4SG02Ps/TqKWxQd7OVdMOXy4TeRcjkE6zUUx4Ls1TQ2naChoYWo2DlYw9vwGQYYvfei6zUDwNSzNq8vltS0HBZeditx8alUlBcwc1YOhw8VBv48KWFqkOMbXYoxNIzlOsu4KKQ//vjjPPXUU4OeU1xcfMGP09LSwmc/+1kWLFjAd7/73QHP+/rXv87aXh9lb2lpIT09HZvNRnx8/AWvYyh8Ph8mkwmbzTYiL6qyo3W8sa2I+iYXb+8spd1j5aZlifzm9QIAHrhzBd/9yZ9xVNad9b789ZioCCtuTwPxsVFMscXz9vsVdHZ7sJhNLMxMY8XSS0lN7duR3jvG1NRUsi+ZO+yxBstIP4djgWKU3sZaLofg5/Ngvn4ObtzIvt/+Fp/XS3tjI1GJibjOVhg3mZh3/fV0u1zMveyyM/7Hqj8T4Wck1GMM9fhgYsQ4XJTLzxSM18+KSz9Ju7WbmJhYnt/6W+pdTZhM4PF5zzjXXzwPt1hJiU6kuq2eSFc4V7Z8godX9nQh9+5MnpNy0Rn3MRF+RiZijP7X0YqsT5KacvZ/08e6ifAcDhfl8v4F9Xfzwo0Uf/gi3V3tGHhO9rSc3thihl6fOmqqNxEZGU+KLZmWmGQuvewmUlLnsG3Lj/lw94tYze3MvfjhwPkT4WdEMYaGsRzjuCikP/roo9x7772DnjN79mymTJlyRjeFx+OhoaHhrO8stLa2csMNNxAXF8drr71GWFjYgOdGREQQERFxxnGz2TyqT7DJZBqxx3xj20ds/ONeVl0zj7s+s4S83Czs+UW8+Ke9uD0+wET2vOmUH2/A5zPw+gaajg4mA6xWC2Ci22PQ0NzB5JR4Zk5LpvxYPfPnTOHZr9/e79zzkYxxtAw22z0U4jsbxSh+Yy2Xw9jI58F6/SzMy+PQli0c//BDTCYT6ZdeSslbb531uq7WVurLytjx05+SfNFFgY1HBzMRfkZCPcZQjw8mRozDQbm8f6P9+pmbOpu1yx/g9l/dT3Wbs+fgwL+OA9Dp7SYyPIo1V/4DTR3N1Lc3cqS+ggzbLN746C9s3PsqmAbebHIi/IxMtBj9r6NQMhGew+GgXD6woP1ufmke9c4yDhW/Q0dH88lPh54+1uLU9yaTmUmJ05mfdT0d7U3s372Jdlc9N972JNPSsykv28G09Owz4pgIPyOKMTSM1RjHRSHdZrNhs519c8mrrrqKpqYm9uzZw5IlSwDYunUrPp+PnJycAa9raWlh1apVREREYLfbiYyMHLa1j1e9Nxv1F3/zcrN4Y1sRhyucNLV2AD3/qKUkRnHC2TrwnZnA7fHi8fZ0yBgYlFXWkZQQzScvn0NcTCTPvLCNhLgo1tyaE3Ibidrzi9j4xz0ArL0nN8irEQke5fKxxZaRwfTFi6ktLcWWmUlievqgI138ag4eZMoll1DncFBkt5OVlxf4OlBRvfnECQ5u3MjCQc4ZSc5eaw3G44uEEuXyseXBpXdTVvcxNW11ePvpRj9dReMxjjZW4XTVU9XcM5YhOSaR7LT5wG3abFJkglAuH3tskzOIik6gq7N1gCJ6X2Hh0czPup4bb3uSN1/9Zp/bDpdsperohxwu2UrmghWB43W1ZXy47x0uvWwFqVNC59P+IqNpbJX1L9D8+fO54YYbuP/++9m1axc7duzgK1/5CnfeeWdgN+njx48zb948du3aBfQk+Ouvvx6Xy8WvfvUrWlpaqK6uprq6Gq/37L+MhqqMGTbW3pPbp6idMcPG1ZfOIioijIS4KBLioogMtxITFYHVYu7ZC+M0ZhOB26wWC1aLmYz0FJISojnhbCF/l4M/vvsR9vwiNr21D3t+EdDTxf3j327jhLN5yGt2VDpZtyEfR6XzQsMfVnm5Waz+7JLAmxMiMjjl8tGTs2YNi++6i/TFi4lKSiI8JmZI1y198EGuuOeeQBF9z8aNFNntA55fUVDA3hdfHPSckTSUNYrI8FIuHx0rMq/l5S/8nC/k3MkN85YxI2EaEdbwAc93+zy8XbKNPUc/pL27k53lu9mwaxOFVcWszX2ADNusUVy9iIx1yuWjq6O9CTARFn72PQuiY5IBeOetH/HRh3/kojlXkXPNmkGvKSp8g7LSdykqfGMYVisyMY2LjvRz8fvf/56vfOUrrFixArPZzO23386zzz4buN3tdlNaWkp7e89GmXv37qWgoGfud8ZpXWrl5eVcdNFFo7b28WDNrTkkJ8SQl5tFRVUDJeU1xESFU1HV0G8To88An6fnnVS3p+cfzckpcfzTXdfywHf/l1ZXF6lJMXzqE3NJiIsKFJv9Y2SirZ7AXPTBRqT4rxmLnd/+NyVEZOiUy0fPsb17OfHRR/g8Hnwez1nP72hu5vVHH+WWH/0IW0YGWXl51JWVUbplC2nZ2WSuWHHGNTNzcrC2t5OVlzcSIZxVVl4ervp6XPX1OB0OdaWLjBLl8tGRYZvFkzc+hsNZzmP27wc6zU9nxUJSTAJOVwMAzZ0tuOvdLElfRHbafNbl/5y8rOtUTBeRPpTLR56zxkFRYU/Dh8UajjGETxg1Nx7lo8I36e5up7OjpwGxqNBOVnYec+ctp6a6hLnzlve5Jiv7Jjy+aLKyV/R53KzsPGyT9fuxyFCEXCE9KSmJjRs3Dnj7RRddhNGr4rts2bI+38vg/EVhR6WTp3+TT2l5LUsXzyI8zHJydvrZHTh0gplpScyYmsDhijo+9Ym5PHL3skA3OvhHyxjkLJoWOHa2QnnvcTQDOVsxXkTGBuXy0VFkt1NbWoqnq+usI138DK+Xlqoq7I89RmdzM1d/+cu0VFdzfP9+tj//fL+F9ElTpzL34YeDMt/OP9YFoGTzZmKSk8nttZGViIwc5fLR4XCWYy/aQr2rkZIaBxFh4fjcXRiGgdFrcLphMmjtcvU5Nis5nafyvoG9aAsb97wKDDwjXUQmJuXykVdUaGfP+xuZl7WK6TMvo7J8N2BisM0vDMNHdEwSyakZ1DvLSLbN5m/v/DclB7eQkjK739EuKalzWHRZHCmpqX0eFyD3ev1+LDIUIVdIl9Gx/rUCDhyqwgCOVjfS1d3TxRgdGUZntwdfP5uPmk09HepWi5n1rxVQdrSe8HALs6cn8dg6O47KOoDASJmHP7+sz8YmZyuUD6Xze6S61lWgF5HxKC07m9TMTJpOnKD1xImeg0P4H5+wqChaTpzA53az/b//mzt/+UugZ+TLWOF0OChYv57D+fk0VlYyd/lylqxeHbSueBGRkeIvgq+al8u8yRkUnSjB4/PR7enqc57X8NHh6exz7FhTNf/21o9xdbnweL3YYpJGc+mD8r9B4O+S93+fnTafwqriwPF3Sv/G89tf4MGld7Mi89pgL1tE5JxlZff8fhoTa+No5V7SL7qc2upDdLY3YWDg87r7va6p8RjtrgauuOYe6mrLKO/awfGKfXS0n31ErrPGgautnnlZqwKPf77U2S4TiQrpclb9FYn9m41ePNPGgjlTOPSxE4sZpqVOChTET2cAFrOJ6KgwfvfGB3R2eaALfvzbvxIRZiFxUjTZmWmBx1t08VT2HThMUwesufXKYRmRMpSu9aE4/e9krI6VEREZTFVhIW1OJzGJibhqajBZLETExtLR1IThG/hTRu6ODjjZXW4JD6eqsJClDz5IVWEhSTNnjonRKQXr17P7hRfwut0YPh8up1Od6CISkvwbhGanzaepo5nSGgc+nxcTZoyzbFbX3NnC5pJtgS71jXtf5/NX/B3Qt5A9O3nmyAbRj95d8nlZ1/GY/fs46srJSJmFo66cHUd281TeN3h++wvs/PgDAGYmTe+32C4iMpbZJmeQe/1afvX87VRV7mdS4nQsZgsxsSm0NFcNeF1nRzMedxc73/05iSkXAT2d6iZMpKUvOmO0S29FhXZKijaz5MrVF1z8Vme7TCQqpMtZ9VckToiLIjzMSmSElUmxEZhMJswWE81tHYSFWeh29zPTywCvYXC0uu+7o13dHqwWE91uL4WlVRSWVrHxj3vYuT8Rd0crJZWtJCfEDkuB+nyK8f29kXD638lwFehFREaTvzu7rqyM+vJywqKi8Hk8gxbR/fwfNrWEhbFn40ZKtmyhzuHAVV9Pzpo1FNntZOXlkTx79pDX4x/DkpWXhy0j44zv+zvndKXvvMP2558nPCYGk8VCVEwM8VOnkvu1rw15HSIi40mGbRZrcx9gXf7PebvkXTrcnfgMH76zFNEBTJj6fH+gqpirn7mZf7zm82zc+zqltWXUuxp54jP/MlLLB87sPoeeNwZ2HEnHFpPEQy//K8U1h5k/eS4PLr2b57e/QEnNYR565V+xxSRxadolPLj0btYXbGLTPjsZKbOoaqkOFNtVTBeR8WDpsp5Pdy689BZcbU7KDm+npfkEYGAyW/qdne7xdOLxdOL1uAkLiyIqOoHpMy6l4kgBVccKA6NdnDUODuy3M21mLqmpqaRNz+ZIyg7Spmef0xr76z73d7RfaGe7yHigQroMylHppL7Jxapr5vUpEq+5NYe9xUcpLXdSfKSGbreXMKuZtvZuYqMjaGju2WgkNjqctvZuoO90L7PZ1Gf8S3xsJBkzUsjO7Nn5e8e+I9y8LIvqmlqy5vcUqEdyfMpg9+0vmtc3uQIbrZ5eONeGoiIyHtkyMshduxanw0H9kSMc3bt3yLPSTRYLEXFxWKOimHnFFQDUORx0NDVhf+wxakpKOLJjBzf98IcQFzek+yyy29lzcgZn7tq1Z3zf3zmn2/zv/07NRx9hjYoiefZsOhoayLrppn5nt4uIhJK8rOuodzXS1NHCwepSPqo+dNZrDIw+xfRurxtHXTnP/vXXONvq6fJ209TRAsCJ5ho2HrSTt3D4u7xP7z63F22hrK6CD6sO0t7dQXHNYbq83dhik9h6eCcYBu3dHew/9hEWs5l7rriDFZnX9twG2GKTqHM1UFJzGHvRFs19H0R/b2KISHBkLljRZ6b5vg9exl9JsVojcXe7+r3OZLKQeckq2tpqWLrsQZKSZwaK3X5FhXZ27/wtx49XYHhvYse7P6WuxhEotg91PEvBjvXs270JV1s9N972JHCqo95Z4yD/7XUD3odGwEgoUCFdBmXPL2LzjhJWf3ZJnwJzxgwbs6enUHS4muRJ0XR1e0mzxdPU2onV0vPLuAkCRfTTRUeG0dHlwev1YTL1dKXvKz7O1oLDJCfE4Kis48W39nLFPBsQEVhL7y7w4Sysr3+tgE1v7aO+ycWTX70RR6WT9a/17DS+PGcuAPVNrj6Pr8K5iIQKW0YG0xcv5nhhIT6PZ0jX+NxuOpub6Wxupjo2luu/9S2qDx7k8Lvv4nO7CYuK4vi+fRT85jd84qGHzrjeP8McYO7y5VQVFpKW3dMR4++UP/3rQMf899V49CgNR44A4OnowARccc89ZOXlBTrVlz74IJkrVpy1s11EZDxKjklkTc4dAKz66Wpau/ovuvRmnLaZXZQ1ghkJaYRbrFQ0HgcMyuo+5je7NvHXql1g6il2ry/YBMCanDvOuwDrcJbzzLZfcrD6EDkzlwSK6Bt2bcKECa/hY3bKDDo9XRx2HqGsroJ8x066Pe7Aug2vj53lu3E4y1mTcwfJMYnUuxoprDrIvMlzA6NvpH/aaFZkbHDWOCjY0fO7cc41a7BNziDM2lMLCQuLIiY2maaGAQrpZjMfl+2ko72J/K6nuX31c2cUtrOy8zji2ElL0wl2vPsz6mocpEzOCBTbe49nycrOO6+C99nuQyNgJBSokC6DGmhkiaPSyZFjdYRZLSzKTMPV0c0tyxdSUl5L8ZFqwEl7Z/eAhXRXezfxsZE0t/VseNTU0gEmE9t2HebS+dNJmhTFgdIqPJ1tHDnRDphoau3AlhSLYfi4/eFfMSUlnoIPK4ChzyXvXXwHAn9uau2gs9sTmP1uzy9i01v78PoM9hYfZfH8dObNSqWkvCbQNS8iMlYNViQeaHzK3OXLObp3L8f37RvSaJfeGj7+mA9ffZXa0lJ8Hg/WqChiUlLoamkJnFNXVsZHb7wRKIDbH3uMY/v2YbZYqCkpoc7hICUjg7ynngqsy19oP11dWRkv3X8/UxYsYOaVV/LXZ5+lzenE29WFz3vyI68mEwnp6YGu9Vceeoije/fibm8nc8WKs3a2i4iMN70Lotlp85kcZ8Pj89Lh7jzLlX11eLr4W/kuwsxhgEFrp4uvvvodaPdiwkR22nzsRVvYtM8O9BTvz6cA63CWc/9L/x8fVZee7Ik3uPt3D9PQ3ojFZKa+vYnZyTO4cuZidhzZjdlkpr69Ea/Ph9Vixe11E2YOw2d4Ka8/yvqCTSTHJAYK5/4/q8t6cP6/L73hIBJcRYV29u3ueYMyJjaZ3OvXcv2N32L7tueJj59C8YE/c2q44ikmk5moqITALPVjlft5acP9TElbQL3zCM1NPcezsvNInZIJJDA38zJOHP9wwPEsgxW8585bTk11yRnz1/2bl86cnYOrrZ5tW57h4IE/U/rRFm5f/RzAWTc3Vce6jAcqpMugBhpZYs8voqq2hcvmTyN9SiKbd5RQUl5LSXkN+0t6EnV0ZBjtnW4MwzhjUoABuDq6MZl6pggYAIaBo7KOw5V1hIVZ4OTolxlpSQC8vaOErm4P+4uP4fUZXJIxhdWfXUJ2ZhrrNuT36UwfqFu9d1d7fZOLF/+0hy3vlTJ7ejJhVgtHjtXjqHSSl5tFfZOLnfvLOXDoBKXltVw2fzoVVY0Ullax4srM4f2LFhEZRoMViQcbn7L8a1/j7X//d1wNDbjq6vC53QQS9SC6XS72/+EPmE5uQOrp6KChvJy07Gxy7r0XH1Dwm9+wf9Mm6srKaKmupqakhMmZmUxfvJi5y5ez/fnnqTtZ1A+s68UXe9YAgXPCY2IoffttDJ+P6o8+4sPXX8fndmMJDydx1izqjxwBnw8Mg6ajR9n1299y4PXXsUZFYQ0PxxoVxa9uv52Ft9zCktWr+3S2i4iMZ70Loo/Zv095w1Emx6ZgMVtwe910efpvcBmI2+fGBHxYdZCm9mbmxl5EbVs9Ww/vZE3OHdS7Gvs8bm+9x4UAZ4wOcTjLeeiVf6X45PiZ6LAoWjtdVDYd73M/zrZ6nv3rr6ludRIXEUt7dwc+w0dceAwdwMzEaZxodZKZOoemjhY27bNT72rkyRsfU3f1EPln7ItIcGVl5+Fqqw/8GSApeSazM64hbXo2UdEJlDt2Ul11kN7FdMPw4Wpz9vreS231IZw1hzEMg9h4GzGxNuyvPEZdzREuu/JLXDx/OfMuWQmcKl73npXe38xz/3mutnoa6yoCI2F27fwtf33nWRKS0jn68Qf4vB4MYFJCGu7uDo4f/ZBtW56hru4IzhOlLM65a8AiuTrWZTxQIV3OS3ZmGhkzUnjwzqXMTEsiOSGG+iYXjso6Mmel0tnVTfmxBubPnkz6lATe3/8xTW19u2E83jM7Hv3/HLjdXqzmk8cMH69sKcTr8+H2+PCeLLC3ujp5+e39/PrV9+nscvOCfReLLk4jfUoiAG+++xFb3isJdJP/+rUCuro9LJgzmfomF02tHbg9PkrLa5k9PZn42AiKy6pZ/1oBa27NCazJajGTOWsyD965lMLSqn4L9yIiY0l/408Guq331yK7na7WVmKSkmg9caLngiHOTAf6dLIbhkH64sUAbP/Zz/jof/+X7pYWit9+G2t4OEkzZhBjs1G+cycdTU0sffBBPnz1VUq3bCEtO5u07GzCo6Npb2yko6mJ7c8/T/mOHVijovo8jr/Q7u3upr68vKeIflJjRQV/ffZZmo8dY9qll7LskUco3bKF4/v393SmX6fuOxEJHRm2WYHRKLcsXAXA1RctYefHe3A4yznRWnvO92nCRG1bPXHh0QB4DR//9+FbzEudM2DHt8NZzmP27+OoKw8c6z065J3Sv/H4Gz+gqqUG38nf/uMj4/F43Wc8vsVkoaWzlaiwSCZFxhFhTaKi8Thur4e5ttlgGHx6fi6PLLs/MGpGRGQ8sk3OCMwc9wsUlq+EG297ktKD7/D2m/9OvfNj3G4X1rBIPP186igmLgWvp4uO9iY62ps4sP/1k6Nc5pKUPIttW37MtPRsqo4V4mqrp6RoM0dSdtBY1/OJf38BvbzsPeyvPMbSZQ9yuGQrH7z3eyIiY8nI/BRp07P5w+/+iQP7/g+vt5uWpiq8vfJ4U0Ml0FPY/+jDP2EY3rP+f4U2LZXxQIV0Oav+ursLS6v6dGf7Z5YnJ8SQnZnG07/J73MfZou5nw8h9Rio2dFn9NRDistr8Z12e2x0OMeqmwLHLWYTLa4uTjhbiI+J5Ppr5mEywYelVRSX1WCxmGl1dQHQ0tbJ9r1HmJmWRO4VGTgbXbS6Oqmtb8Pr8/H2jmL2Fh+jtLwWMFhySTpPrc0jY4aNFVdmsm5Dfp9Z6SIiY41/E9Gh3Nb7e39Rff/LL/fcaDIRY7PR0dgIJhO+7m6s0dHg9eLp6ur/wS0Wpsybx5QFC8hZs4b316+nfPduutvaAOhqbWXaNdcwed48dm3YgLe7m5riYg5t3Yq7sxNPRwf5Tz+Nu7OTttqeos/BP/+ZiPh4rFFReLsH6aj0j3UBopKSuPT220mdN4+9GzeSPHs2adnZ1JWVYQDJs2drtIuIhBz/eJfVS27jlft+wTfffIoPqw6SEBl/TvdjNpnxGT58GPh8Hlo7XZgie94krXXV8923foTX5+ONA29z9exPsHzu1RRWFWOLSeLZv/6als5WsqbOIy/rOioajvFGURwv73+TLcXvUniiGI+v734c1a21RFrD+xyzmi1clDSD6tYa2rpcHG2uYkHq3J4COgZgUOIsIzI8EnvRFpbPvbrPaBcREb/xMjLk9HX6C8oxsTZ+9fztxMdPoauzlSlp86ipLh3wflqbT+AfAxMeHktLUzXRMUnMnHUV7/31F9RWFZA6JYOuzlYmpy3A6/NgGEZg7Ip/Q1GrNYK21lqOfryX8Iho3O523O526uqOkL/5aY5W7KUnH5vw9vo9/HTubhdR0Umk2GYxd97yPl3wVccKA/H6Ny0VGctUSJezOn2TT+h/drp/DMy6DflUVjUQFmahrLKOg2XVg77xeA7NjgGdXe4+xXV/l7phQFt7F3/b46C6rqdoYzKZaO889c5oZISVuiYXxWU1nHA2097pxmw24fP5MAyoONFElbOFiy9K5epLZ7Hm1pw+necDzY0XERnv/EX1urIy6srKAIiz2ZicmUnLiRPUlZURGRtLe0PDgPdhAlqdToyDB/n13/0dcVOn9rk9PDqa8JgYOpqaiJsyhaajRzF8Ptrrez7KarZaaaqqOtURD3S3tQUK8UN11f33s+jmm9n2zDM0nzhBQ2UlLdXVNFZUBEa6+GfFi4iEioHmXcdGxhJpjaDb68ZnnH0fjNPP8fVqh7GarbR3d+A1fJQ6j3Cs+QQlNQ4qGo/i8Xo50VKLLTaJmPBo7n/p/6O1s42jTcf7bajprcPT9w1ar89HUXUxBmCLSaK120VidAIHThTjNXxMibUxOTaFKGsUP92+gTsuy+PJGx87a2wiMvGMl5Ehp6/TX1j+1fO387FjJ9PSL2XJlauJibXx13eepanh2CD31pN1XW21uNp6mlNqqg8xKXE+PsNNvbOcOZmfpOrYh7Q0naC1uYbW5moOl2yjufE4bnfPXnXQUwh3d5/a6LTB+THu7nZOtUqaMZvN+AbYZ8lkMhMRFUvVsQO8/ea/Myvj6jO64Mfy8yLSmwrpclbnWjjOzkwjOiqcruZ2rGFm3B4TPgxMgMViwe059U7lQF3qZ+PxDnyV12cEiugA3e5Tj2cxm3C7vfhOFt4bWzoChfykSVFYLRYaml2YTCauvnQWT371xjPuf6C58SIioWLRbbdRffAgBpCQns7H773H1IULaW9ooL2x8dSGnv0wvF5ctbW4TnaTt9TUkLBwYeD27rY2St56C4AYm63nHdBeH03yeTx9iujna8fzz3PwzTdxHj6M1+0mfupU4qdMYfK8eYGNVtWJLiKh5vR512ty7iA5JpHstPk8v/0FdlcW0unpDHScnwsTJiwmM1PjJwdmmcdHxBIXEcOHVQeZdHI8i4FBU0cLm0u2YZzlN/3EqAQaO5r6PS8mPBqf4aPD3Um7uwNXdweFVR9xzaxP4GxroNRZhsVkpsPTcU5xiMjEM15Ghgy0zqXLHgx8zVywgjdf/SbtbfVERU+i3TVwg8vpDMMLJjNgwu1up+LILsDAZDJjGF7qah2nX9Hv/bS76omOSQ7cT1hYBCazme6uM0d0mc1WpqVn0+3ppMlbSc2JYqakLQi8IXBg/+t95rOLjHUqpMtZ9Vc47q9L3W9rwWHqm9qJj42go9NNfEpkYHb58py5rP2P1zjhbMVqMQ1aEB8JXp/BibrWwPcmkwnjZPFmxZWZ3LZyEU//Jp/Z05P7zEkXEZlIqgoL6WptZcnq1bjq6/n4vfdoranp04luMptPzSofbENSw+gzt7w3l9N56pwLYTJhjYxk+pIlHN+zB3d3N12trdQUF5M8Zw6R8fGkzJ5NRUEBS1avxpYxdj/SKyIynHoX1mcmTSf3v+8AejrOTZjOWujuzYsPr+ELFNGtJgsLplzM9vJdADR3nvodu7u/eeeYsVgseL1evPgwm8w0dTSfcZ7VbAEDOtw9Xe/Q889ETHgUrV0uth/Zzfc+/SgltT2fnPKPldFIFxEZSDBHhpzLWJmB1pm5YAWZC1b0OWYyW5gx6wocpe/icZ/DG4qGD3+BvKO9AbMlrKfAfo462hsxTubonu71/tmmZnL76udoqK/g9U2P0tHeRFR0ArnXryX/7XV9Ni4VGQ9USJfzMpQu9ejInjmHbreXS+dND8xRN5l6Ph7kHeUien98PoOoCCs3LsvikbuXYc8vwtnQxnVXZWojURGZsE7fiDQmOZmiN94I3B4WFYW794z0Cy2EXyjDwNPZSUdjI8mzZ+M8dAivYWAJC2PusmXc+OSTOB0OjXIRkQktwzaLhz95Hz/ZvoH4yFga2puIj4yjvbuDtm4XJkyY6DvGZTAew8vO8t1DfnwvPrxeHxHWcLyebkxG/49lNlnA1FPit2DB7XWzKG0+Uyal8vL+P9LW7eL1A5t55b5fBK5ZkXntkNchIjKaRmKszNx5y6mpLiEyMo6IyFhMJjMed0egsH0ufP288TkUAz+W6eR/Pbe7WpwU7FhPzjVr+MI/vhx4UwHGzycFRHpTIV3Oy2DjTdbcmkNJeQ37io/j8/mwJcViS4xh3YZ86ptctLm6CLNaSIiLpKPLTVv7IBvHjZCemeg9v7hPn5zAT77198CpNwayM9NYtyG/zwarIiITkS0jg6y8vMDM9CkLFpA0ezZ/e+453O0Dd5+MKLP5zC53w6C2uBiT2YzP07OJXdyUKeSsWQMMvgGriMhE8eiKL/Poii/jcJbz0Mv/SqmzjCnxqZTVf4xhGFgtYedUVBlq0b23bo8bi8mMxWzB6+3J5bHhMXR5unD7PLh9bkwGxEXG4fa6MZlMzJsylzU5Pd30R+oqeXDp3ef8uCIiwTASxeKqY4XU1Tjo7m5n0WW3UH2imI8dO8FkgfPoLh9eBv6Od7MljK7ONvbt3kRMbDK516/t82aCNheV8UiFdBl2GTNsPLU2j/WvFbC3+ChVtS28vvUAFVWNrLpmHomTomk90UiLq4u83Cze3e3A2dAGpp5RKz7fhXU2WswApsAGpP0Js1iIiQun1dXJpfOn91m7f8PUgUbXnAtHpRN7fpEK8iIyrhTZ7ezZeLJzZu1aiux2KgoKmJmTQ0t1Nc6ysj5F9D5jXkaYNSqKWJuN7ra2wKgZS1hYz6/rhkFkfDwdLS1EJyVx0w9/qDEuIiL9yLDNYnH6Qhx15Vw2/RIyUi5iR/kuDAO8Pm+/89NNQ7jfoYyLyUiZSbg1nJKassCxTncnk6Li6HR3YTKZSY5J5J8/+YXA+JY1OXeQYZvFT/7+++cUp4hIsI1EsTgrO48jjh3UVJUA0FR/FMMwejb8DP4H/wOslnCWXLkaUNe5hA4V0mVEZMyw8eRXb+Sd90t5/qXt3LJ8Ic5GF3m5WSzPmcvTv8nHlhjDwbJq2tq7iI4Kw9XhJiYqDFd793ltQOrX09hy5j2YAKvVfHImukG328uV2RfxyN3L+pznqHRS3+Ri1TXzhrzB6kAGmyUvIjJWnT7axf+1dMsWju/fj8li6TnRZMJksTApLY2mo0dHZcSLp7OTpqNHmbJgAeGxsbjq65m1dCnVBw6QtmgRkXFxFP/5zySmpwOQv25dYHNREZGJyOEsx160hbys68iwzQocA7jjsjzW5NzB+oJNGOUQHxlLl6frjEK62WTu89u1xWTBZ/jOKJr3/t6MKdCxbjVbuShpOhHWCG5csILffvAKYRYrsxLSibCG42xroM7VwGXTs1g8vWeD6qtmLeHzV/zdcP91iIiMe7bJGeTd/hT2Vx7j6Md78HrcRETG4vP5MBthPWMPPd34x6uMLhNmsxWLxcrS5V9hxQ2PBmENIiPHHOwFSOhwVDpZtyEfR6UzcKywtIqKqkacjS7W3pNLxgwbK67M5M8/+zKXzptORVUDHq+PiPAwLGYTbo+PmOjwEVmf1drzcrdaLCQnxvDpa+fz1Nq8MzrF7flFbN5RQnJCzAV3keflZrH6s0suuCAvIjKa/GNQ/MVn//e5X/saF119NZfffTeRkyYRERfH9OxsEmfMwBoZ2bPpqJ/5zF8xzNae9+9NFgvRycnntzjDwGSxsODGG7l89Wpik5NxlpTQWl1NTXExi267jWmXXUZLVRXbn3+ePRs3UmS3n99jiYiEAHvRFjbueRV70ZY+xzaX5JMckxgorltMZq6Z/QkWpS0g3BJGhCWcCGt4T0H8tMJ6ZFg4c1JmEmYJ6/cxrWYLV190OWFmKyZMrMi4hghrBIedR1i/63850VJLhDWCJz7zL+Q/9DL/ddv3uGbWJ/ha7pdIjklkc0l+n/WKiEhf/mL65LT5hIVFMD/rBtIvWkx0dALJk+dwqrmwv88Tne0zRkP5DNKZrNYI0mcuxjC8xMansuiym3HWOMh/ex3OGsd53afIWKOOdDkv/Y0s6a/7erCZ43m5WZQdrWPHviN0dLq5dP50TjibqaptCTyOCbBazHR7en55N5s4p48qWS1mIsIsdHR7MIye2ehJCdEYBsxJT+m3UD6UjVSHarBZ8iIi403mihVkrlgBwFVf+AJFdjuu+np2v/ACnu5uwqKjsVitdLe3c/GKFdQ5HNSVl4PJRHRSEl0tLVgjIrBGRREeE0NXWxve7u4+newmq5XMlSup3LWr53aPB2tYGJ5em5saHg8f/Pa3zL7mGuatWkXqvHn89dln6W5vp6qwkLynnqLIbictO5uqwkJtMCoiE1pe1nV9vvZ3bE3OHSTHJAa+f8z+ffYc/RCPz8PslJm0dbad3Iq0Z3yLx+flWNMJwsxWws1WOjxdGIaByWQiJSYJi9nMkYZK3L6ePStKnGVUtdTg9rqJj4wnKiyK9u4OCquKWZF5beA/gJlJ089Yr4iInMk2OYO/u+u5Pht4FhXaqasto8FZjs/n6dnc2WTCao3EZxiAiTBrFG5326k7Mpl6/T5uYUraPDyeLupqy+jv0/5+Fks40bEp+HxuurvauXb5V4iLT8VZ66Crs5Wiwp5mluHebFUkmFRIl/MyWNG8dwF6sJnjGTNszElPYfveI2TNncpTa/OoqGrgi99+EVeHG7PZhMkMSZOicDa2YzJB5qzJdHd7OHKsHqvVjNdr4PFvUhQVRltHz+ZIURFhfOoTc7juqkw2/mkvtsQY4mIiSYiLYnnOXApLqwYslKv4LSJydv4udafDQUdTE0d27KCjqQmf18u0RYtITE/ninvu4dDWrXSFhXH0L3+h4+RMc09XF90uFxddeSULb7mF93/9a9obG/F2djLjiitwOZ20NzZiMpkIi4rC8HgwW62BTUTDYmLoaGqi5O23WfqP/8gVn/88s666iiK7PTDGxb+xqL/wLyIyUWXYZrE294EBj/U3+uWpvG9w/0v/H4edR7hsehYPf+qL/P6vr1Bf1kRVay1dnm5iw2P49IJcwMSJ5hoO1ZaxOH0h91zx9xRWFbP/2Ec42+pJjJrEP3/yC7xfsY+D1YdIT5hKXGQcCVHx/RbL+1uviIj07/QZ7LnXr8VZ4yAldQ4xsTb2FmwkJs5GYlI6jQ3HcNbWEhYeidvdRnRMEtExyRgYdLgaaXfVEx4eSUJSOl1dLtpdDYRHxBITl0JzwzGsYVG0Np8ADAzDOHldPVOnZdHcVMXHR3Yyeco8wsOjSZmc0WcuumakS6hQIV3Oy2BF86Ge759FfuOnLmHNrTlkzLCRMcPGL//tLp5/aTtXXzqTyqMnmNoO735QRmx0BOlTEgBwdXbjdntZckk6WwsO0+32Eh5uJdzjIzY6nG63j/Qpiby+9QD7io8RFxPBz759ByuuzAQIfBURkQtjy8jg73/yE5wOBwXr1wPQ0dTE3hdf5Njevdz27LPsf+cdmo8fxzAMwmJi8HZ2EpWYyNIHHyRp5kxqS0o4tncvDZWVuJxOYmw2Ymw2OhobezrWfT1Tdk1mM5Pnz+f6b32Lw1u3Aqfmt/cunouIyND5R78AgQJ2hm0WV8+6nIrGYxypqwTg5oXX447wUVzj4JCzjMXTF5IQNYnNJfmYMFPf0YSru4Oq5hpe3PM6d16Wx6XTLwkUy52uBhKi4tm0z47X8HHZNI0+FBEZCb2L61dc/fnA8drqw+zf9w6T4mMoKnydpcsepOpYIXve38jF85dTUb6L1pYaDh3ccrJB3ceUaZeQd/tTFBXaSZuezRuvPE5Tw1HiJk0hOjqRhvoKkm2zCYuIpq7GweQp87jimnvIys7DNrlnTKQ60SWUqJAu52U4urb9s8hXf3ZJnxErK67MZMWVmfh8Pmpra2ntMtP2zBvsKz5O/q6euVphVguZs1JJn5LII3d/ipf+vI/0KQkcOHSCqbZJHKtuAuDBO5fy4aEqWto6efyZN3jxP5MueO65iMhE53Q4+nR/Q08h+8YnnwTgzW9+E6/bTXVxMW98/etk338/UYmJuF0ukmfOJDopiTqHg6rCQqoKC9m3aROG10t0cjLVxcUYRUVk5eURlZBAR1MTAJ2trbicTnK/9rU+I2ZEROTC9Df6BXrGvew9doDSmjJ+U/AH0iMm83bpNu5acitL53yCjXteJT0xjdVLbsPw+Vhf8L9Ut9Ty1Ds/obatnp/t/B1//vJvybDNYl3+z9mwaxNp8VO4ft6nOFJXiaOupxNe3eciIqMjJXUOiy6LIzU1lZxregrsSckzgVMd469sfIjqqmImJU4jMjqepcsepKjQzp73N+LKqmfmrCuIjU0hd9XXSEqeecZYmd4FdJFQpEK6jIqhjoLpz5z0FJ5am8czL2zjYNkJkuKjae9yM3t6Mpt3lGBLiqWppYMrFs5k6eI5ZGemBUa3ZMyw8bNv38Hjz7xBe0c39vwijW0REblARXY7ezaenHXYTxd4zpo1ABzdu5e6sjIaysu5+T/+gx0//WmgC91fiAeoKyuj/sgRFq9ezZ6NG3GWlhKVkBAozAPkr1vHno0bqSosDBTRS995h+3PP8/SBx9UYV1E5DwNNEolwzaLxdMX4nCWAzArKZ2Zielkp83vM8c8wzaLb775FPUdTdS66slImUWXx43FbA4UyvOyrmPHkd046sq5bt4neWTZ/YFxMiIiEjynj4a5ffVzZxTE/cV2V1s9FUcKWHLlajIX9PzuffpYGZFQp0K6jIqhjILpbwPT3ufOSU+h4MMKkibF4GxwsXh+Oqs/u4Syo3U4KpwkxEUF7q/36JYVV2by4n8mBe5bREQujL8APtAmnv7udKfDwQG7nWk5OczNzmbeypWBc/zz1YvsdqISEmhzOnE5nfzdc8/1KbIP9pjbn3+ej3fuBDQLXURkJPg3Ib3pkpW8s/+vVDQeDWwQenrxPcISTkJsPE985mvMTJrep1CeYZvFU3nf6DOLXZ3oIiLB56xx9Cmcn15Yh1PFdmeNg5jYZM07lwlNhXQZFUMZBdNf13pv/iL46R3njkonc9JTBi2SawNREZHhM9R55LaMDJY9/DC1tbX93u7vbJ+3ahVLVq8+Y6PQsz3m0gcf7PNVRETOX3+bjvoL3j6fj7aZrbRbu8nLuu6Mc/0F997XDrbBqYiIjA3+sS1w9o7y/orsIhONCukyZpxt1EvvYnjvjvNzKZIP1vUuIiKjq3eXuX/Wem/9zWLvTbPSRUSGT3+bjvY2ddJkHp77RcxmM+vyf97nXBXJRUTGJ393ubrMRYZGhXQZM0aja/xsXe8iIjJ6ztbZfrZZ7CIiMnwG2nT0dA5nOfWuRlbNy9WMcxGRcU5d5iLnRoV0mVCGusGpiIgE39lmsYuIyPAZale5vWgLm0vyWb3ktsAYFxEREZGJQIV0mVA0K11EZPwY6ix2EREZHWV1H6sbXURERCYsc7AXIHI2ZUfrWLchH0elM3DMUek845iIiIiIiIycN4r+wuaSfJJjEtWNLiIiIhOOOtJlzHtjWxEb/7gXODXXXLPORURERERG101ZK8F09jnqIiIiIqFIhXQZ825algWY+sw116xzEREREZHRNSfloiHNURcREREJRSqky5g3Jz3ljK5zzToXERERERERERGR0aIZ6TImOSqd/Pi32zjhbA72UkRERERERERERGSCUyFdxiR7fhEv/mkvBR9WBHspIiIiIiIiIiIiMsFptIuMST2zzw1yFk0L9lJERERERERERERkglNHuoxJGTNsPPz5ZUy1TQr2UkRERERERERERGSCUyFdRERERERERERERGQQKqSLiIiIiIiIiIiIiAxChXQRERERERERERERkUGokC4iIiIiIiIiIiIiMggV0kVEREREREREREREBqFCuoiIiIiIiIiIiIjIIFRIFxEREREREREREREZhArpIiIiIiIiIiIiIiKDUCFdRERERERERERERGQQIVdIb2ho4HOf+xzx8fEkJCRw33330dbWNqRrDcPg05/+NCaTiddff31kFyoiIgNSLhcRGf+Uy0VExj/lchGRU0KukP65z32Ojz76iC1btvDmm2/y17/+lQceeGBI1/74xz/GZDKN8ApFRORslMtFRMY/5XIRkfFPuVxE5BRrsBcwnIqLi3nrrbfYvXs3l19+OQDPPfccn/nMZ3j66adJS0sb8Nr9+/fzox/9iA8++ICpU6eO1pJFROQ0yuUiIuOfcrmIyPinXC4i0ldIdaS/9957JCQkBBI8wMqVKzGbzRQUFAx4XXt7O6tXr+YnP/kJU6ZMGY2liojIAJTLRUTGP+VyEZHxT7lcRKSvkOpIr66uJjU1tc8xq9VKUlIS1dXVA173yCOPcPXVV3PzzTcP6XG6urro6uoKfN/S0gKAz+fD5/Odx8rPnc/nwzCMUXu8YAj1GEM9PlCMI8VsDqn3QM8wWrkcgp/P9TMSGkI9xlCPD4IXYyjnc+Xy0KIYQ0Oox6hcPvwmUi73P5Z+RsY3xRgaxnKdZVwU0h9//HGeeuqpQc8pLi4+r/u22+1s3bqVffv2DfmaH/zgB3zve98747jT6aSzs/O81nGufD4fzc3NGIYRsv9wh3qMoR4fKMaRMl67OsZaLofg53P9jISGUI8x1OOD4MU4HvO5cvmZ9DMSGhTj+KdcPnTK5f3Tz8j4pxhDw1ius4yLQvqjjz7KvffeO+g5s2fPZsqUKdTW1vY57vF4aGhoGPAvZOvWrZSVlZGQkNDn+O233861117Ltm3bzrjm61//OmvXrg1839LSQnp6Ojabjfj4+CHFdKF8Ph8mkwmbzRbSPzihHGOoxweKUfoaa7kcgp/PJ8LrRzGOf6EeH0yMGIeLcvmZJsLrRzGGhlCPMdTjG07K5f0L9ddQqMcHijFUjOUYx0Uh3WazYbPZznreVVddRVNTE3v27GHJkiVATxL3+Xzk5OT0e83jjz/OF7/4xT7HFi5cyDPPPMNNN93U7zURERFERESccdxsNo/qE2wymUb9MUdbqMcY6vGBYpRTxlouh7GRzyfC60cxjn+hHh9MjBiHg3J5/ybC60cxhoZQjzHU4xsuyuUDC/XXUKjHB4oxVIzVGMdFIX2o5s+fzw033MD999/Pz372M9xuN1/5yle48847A7tJHz9+nBUrVvDCCy9wxRVXMGXKlH7fSZ0xYwazZs0a7RBERCY85XIRkfFPuVxEZPxTLhcR6WtslfWHwe9//3vmzZvHihUr+MxnPsPSpUv5+c9/Hrjd7XZTWlpKe3t7EFcpIiKDUS4XERn/lMtFRMY/5XIRkVNCqiMdICkpiY0bNw54+0UXXYRhGIPex9luFxGRkaVcLiIy/imXi4iMf8rlIiKnhFxHuoiIiIiIiIiIiIjIcFIhXURERERERERERERkECqki4iIiIiIiIiIiIgMQoV0EREREREREREREZFBqJAuIiIiIiIiIiIiIjIIFdJFRERERERERERERAahQrqIiIiIiIiIiIiIyCBUSBcRERERERERERERGYQK6SIiIiIiIiIiIiIig1AhXURERERERERERERkECqki4iIiIiIiIiIiIgMQoV0EREREREREREREZFBqJAuIiIiIiIiIiIiIjIIFdJFRERERERERERERAahQrqIiIiIiIiIiIiIyCBUSBcRERERERERERERGYQK6SIiIiIiIiIiIiIig1AhXURERERERERERERkECqki4iIiIiIiIiIiIgMQoV0EREREREREREREZFBqJAuIiIiIiIiIiIiIjIIa7AXEAoMwwCgpaVl1B7T5/PR2tpKZGQkZnNovh8S6jGGenygGEdSXFwcJpNp1B5vohjtfK6fkdAQ6jGGenwQ3BiVz4efcvnwU4yhIdRjVC4PLaqzDL9Qjw8UY6gYy3UWFdKHQWtrKwDp6elBXomITBTNzc3Ex8cHexkhR/lcREab8vnwUy4XkdGmXD78lMtFZLQNJZebDP/bfHLefD4fVVVVo/oudEtLC+np6Rw9ejRk/8EO9RhDPT5QjCNJXS8jY7TzuX5GQkOoxxjq8UFwY1Q+H37K5cNPMYaGUI9RuTy0qM4y/EI9PlCMoWIs11nUkT4MzGYz06dPD8pjx8fHh+wPjl+oxxjq8YFilPEjWPl8Irx+FOP4F+rxwcSIcSJQLh85ijE0hHqMoR7fRKE6y8gJ9fhAMYaKsRhjaA7TEREREREREREREREZJiqki4iIiIiIiIiIiIgMQoX0cSoiIoLvfOc7REREBHspIybUYwz1+EAxipzNRHj9KMbxL9Tjg4kRo4ycifD6UYyhIdRjDPX4ZOSF+mso1OMDxRgqxnKM2mxURERERERERERERGQQ6kgXERERERERERERERmECukiIiIiIiIiIiIiIoNQIV1EREREREREREREZBAqpIuIiIiIiIiIiIiIDEKF9HGioaGBz33uc8THx5OQkMB9991HW1vbkK41DINPf/rTmEwmXn/99ZFd6AU41xgbGhp46KGHyMzMJCoqihkzZvDP//zPNDc3j+KqB/eTn/yEiy66iMjISHJycti1a9eg5//hD39g3rx5REZGsnDhQv70pz+N0krP37nE+Itf/IJrr72WxMREEhMTWbly5Vn/TsaCc30e/V566SVMJhO33HLLyC5QxhXl8/7PVz4PLuXygSmXS3+Uy/s/X7k8+EI9nyuXy3BSLu//fOXy4Av1XA7jOJ8bMi7ccMMNRnZ2tvH+++8bf/vb34yMjAzjrrvuGtK169atMz796U8bgPHaa6+N7EIvwLnGeODAAeO2224z7Ha74XA4jHfeeceYO3eucfvtt4/iqgf20ksvGeHh4cavf/1r46OPPjLuv/9+IyEhwaipqen3/B07dhgWi8X4j//4D+PgwYPGt771LSMsLMw4cODAKK986M41xtWrVxs/+clPjH379hnFxcXGvffea0yaNMk4duzYKK986M41Rr/y8nJj2rRpxrXXXmvcfPPNo7NYGReUz8+kfB5cyuUDUy6XgSiXn0m5PPhCPZ8rl8twUy4/k3J58IV6LjeM8Z3PVUgfBw4ePGgAxu7duwPH/vznPxsmk8k4fvz4oNfu27fPmDZtmnHixIkxneAvJMbeNm3aZISHhxtut3sklnlOrrjiCuOf/umfAt97vV4jLS3N+MEPftDv+XfccYfx2c9+ts+xnJwc40tf+tKIrvNCnGuMp/N4PEZcXJyxYcOGkVriBTufGD0ej3H11Vcbv/zlL4177rlHv7BLgPK58vlYpFzeP+VyGYhyuXL5WBXq+Vy5XIaTcrly+VgV6rncMMZ3Ptdol3HgvffeIyEhgcsvvzxwbOXKlZjNZgoKCga8rr29ndWrV/OTn/yEKVOmjMZSz9v5xni65uZm4uPjsVqtI7HMIevu7mbPnj2sXLkycMxsNrNy5Uree++9fq957733+pwPsGrVqgHPD7bzifF07e3tuN1ukpKSRmqZF+R8Y/y3f/s3UlNTue+++0ZjmTKOKJ8rn481yuXK5XLulMuVy8eiUM/nyuUy3JTLlcvHolDP5TD+83lwfwpkSKqrq0lNTe1zzGq1kpSURHV19YDXPfLII1x99dXcfPPNI73EC3a+MfZWV1fHE088wQMPPDASSzwndXV1eL1eJk+e3Of45MmTKSkp6fea6urqfs8favyj7XxiPN1jjz1GWlraGf+wjRXnE+P27dv51a9+xf79+0dhhTLeKJ8rn481yuXK5XLulMuVy8eiUM/nyuUy3JTLlcvHolDP5TD+87k60oPo8ccfx2QyDfrfUH9QTme329m6dSs//vGPh3fR52gkY+ytpaWFz372syxYsIDvfve7F75wGXE//OEPeemll3jttdeIjIwM9nKGRWtrK5///Of5xS9+QUpKSrCXI6NI+Vz5fKJSLpdQolyuXD6RhVo+Vy6fuJTLlcsnslDL5TD28rk60oPo0Ucf5d577x30nNmzZzNlyhRqa2v7HPd4PDQ0NAz4UaKtW7dSVlZGQkJCn+O333471157Ldu2bbuAlQ/dSMbo19rayg033EBcXByvvfYaYWFhF7rsC5aSkoLFYqGmpqbP8ZqamgHjmTJlyjmdH2znE6Pf008/zQ9/+EP+8pe/sGjRopFc5gU51xjLysr4+OOPuemmmwLHfD4f0PPOf2lpKXPmzBnZRUtQKJ/3UD7vMZ7yuXK5crmcolzeQ7m8x3jK5RD6+Vy5XIZKubyHcnkP5fKxZ9zn86BMZpdz4t8g4oMPPggc27x586AbRJw4ccI4cOBAn/8A47/+67+MI0eOjNbSh+x8YjQMw2hubjauvPJK41Of+pThcrlGY6lDdsUVVxhf+cpXAt97vV5j2rRpg26CceONN/Y5dtVVV435TTDOJUbDMIynnnrKiI+PN957773RWOIFO5cYOzo6zvi5u/nmm43ly5cbBw4cMLq6ukZz6TIGKZ8rn49FyuV9KZfL2SiXK5ePVaGez5XLZTgplyuXj1WhnssNY3zncxXSx4kbbrjBuOyyy4yCggJj+/btxty5c4277rorcPuxY8eMzMxMo6CgYMD7YAzvJm0Y5x5jc3OzkZOTYyxcuNBwOBzGiRMnAv95PJ5ghRHw0ksvGREREcZvfvMb4+DBg8YDDzxgJCQkGNXV1YZhGMbnP/954/HHHw+cv2PHDsNqtRpPP/20UVxcbHznO98xwsLCjAMHDgQrhLM61xh/+MMfGuHh4cbLL7/c5/lqbW0NVghnda4xni6Yu0nL2KR8rnw+1iiXK5fLuVMuVy4fi0I9nyuXy3BTLlcuH4tCPZcbxvjO5yqkjxP19fXGXXfdZcTGxhrx8fHGmjVr+vxQlJeXG4CRn58/4H2M9QR/rjHm5+cbQL//lZeXByeI0zz33HPGjBkzjPDwcOOKK64w3n///cBtn/rUp4x77rmnz/mbNm0yLr74YiM8PNy45JJLjD/+8Y+jvOJzdy4xzpw5s9/n6zvf+c7oL/wcnOvz2Jt+YZfTKZ8rn49FyuXK5XJulMuVy8eqUM/nyuUynJTLlcvHqlDP5YYxfvO5yTAM4/wHw4iIiIiIiIiIiIiIhDZzsBcgIiIiIiIiIiIiIjKWqZAuIiIiIiIiIiIiIjIIFdJFRERERERERERERAahQrqIiIiIiIiIiIiIyCBUSBcRERERERERERERGYQK6SIiIiIiIiIiIiIig1AhXURERERERERERERkECqki4iIiIiIiIiIiIgMQoV0EREREREREREREZFBqJAuMko2bdrEJz/5SeLj40lISCAvL49Dhw4Fe1kiInIOlMtFRMY/5XIRkdCgfC6jzWQYhhHsRYiEMq/Xy7333svvfvc75syZw+WXX47D4WDPnj2kpqby0UcfkZKSEuxliojIIJTLRUTGP+VyEZHQoHwuwaKOdJER9uijj/K73/2O733vexw6dIiXXnqJDz74gC996UvU1tby7LPPBnuJIiJyFsrlIiLjn3K5iEhoUD6XYFFHusgI+tvf/sYnP/lJ7r77bjZs2NDntuLiYhYsWMAnPvEJdu3aFaQViojI2SiXi4iMf8rlIiKhQflcgkmFdJERtGzZMnbu3El5eTnTpk3rc1t9fT0pKSmkpqZSU1MTpBWKiMjZKJeLiIx/yuUiIqFB+VyCSaNdREbI4cOHeffdd7npppvOSO4A7e3tQViViIicC+VyEZHxT7lcRCQ0KJ9LsFmDvQCRUPXyyy8DcPToUe69994zbm9qagIgMTFxFFclIiLnQrlcRGT8Uy4XEQkNyucSbCqki4yQbdu2AbB7925279494Hlz584F4IMPPuDZZ59l586dlJWV8c1vfpN///d/H42liojIAJTLRUTGP+VyEZHQoHwuwabRLiIjZO/evURHR2MYRr//3XPPPQAsWbIEgB07dvD++++zdOlSJk2aFMyli4jIScrlIiLjn3K5iEhoUD6XYFMhXWQENDU1UVdXx9SpUwc8Jz8/H4Dc3FwAHnroIQ4dOsRvfvMbEhISRmOZIiIyCOVyEZHxT7lcRCQ0KJ/LWKBCusgIaGxsBCA2Nrbf2/fs2UNlZSWTJ09m6dKlAJjN+nEUERlLlMtFRMY/5XIRkdCgfC5jgV5RIiPAZDIB0N3d3e/tv/jFLwD4whe+gMViGbV1iYjI0CmXi4iMf8rlIiKhQflcxgIV0kVGwLRp07BarRw5coSOjo4+t5WWlrJ+/XoSEhJ49NFHg7RCERE5G+VyEZHxT7lcRCQ0KJ/LWKBCusgICAsLIzc3l66uLp555pnA8WPHjnHrrbfS3d3NT3/6U5KTk4O4ShERGYxyuYjI+KdcLiISGpTPZSxQIV1khHz3u9/FYrHwzW9+k9zcXG699VbmzZtHaWkpzzzzDHfeeWewlygiImehXC4iMv4pl4uIhAblcwk2FdJFRsjVV1/NW2+9xZVXXsmuXbv429/+xsqVK9mxYwcPP/xwsJcnIiJDoFwuIjL+KZeLiIQG5XMJNmuwFyASylauXMnKlSuDvQwREbkAyuUiIuOfcrmISGhQPpdgUiFdZIxwOp28++67ALS3t1NSUsLLL79MTEwMn/70p4O8OhERGQrlchGR8U+5XEQkNCify3AzGYZhBHsRIgLbtm0jNzf3jOMzZ87k448/Hv0FiYjIOVMuFxEZ/5TLRURCg/K5DDcV0kVEREREREREREREBqHNRkVEREREREREREREBqFCuoiIiIiIiIiIiIjIIFRIFxEREREREREREREZhArpIiIiIiIiIiIiIiKDUFAz5fcAAAB7SURBVCFdRERERERERERERGQQKqSLiIiIiIiIiIiIiAxChXQRERERERERERERkUGokC4iIiIiIiIiIiIiMggV0kVEREREREREREREBqFCuoiIiIiIiIiIiIjIIFRIFxEREREREREREREZhArpIiIiIiIiIiIiIiKD+P8B951xJvGoTFYAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 1500x600 with 4 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Set the number of posterior draws you want to get\n",
    "num_samples = 3000\n",
    "\n",
    "# Obtain samples from amortized posterior\n",
    "conditions = {\"x\": np.array([[0.0, 0.0]]).astype(\"float32\")}\n",
    "\n",
    "# Prepare figure\n",
    "f, axes = plt.subplots(1, 4, figsize=(15, 6))\n",
    "\n",
    "# Obtain samples from the approximators (can also use the workflows' methods)\n",
    "nets = [\n",
    "    flow_matching_workflow.approximator, \n",
    "    consistency_model_workflow.approximator,\n",
    "    affine_flow_workflow.approximator,\n",
    "    spline_flow_workflow.approximator\n",
    "]\n",
    "names = [\"Flow Matching\", \"Consistency Model\", \"Affine Coupling Flow\", \"Spline Coupling Flow\"]\n",
    "colors = [\"#153c7a\", \"#7a1515\", \"#157a2d\", \"#7a6f15\"]\n",
    "\n",
    "for ax, net, name, color in zip(axes, nets, names, colors):\n",
    "\n",
    "    # Obtain samples\n",
    "    samples = net.sample(conditions=conditions, num_samples=num_samples)[\"theta\"]\n",
    "\n",
    "    # Plot samples\n",
    "    ax.scatter(samples[0, :, 0], samples[0, :, 1], color=color, alpha=0.75, s=0.5)\n",
    "    sns.despine(ax=ax)\n",
    "    ax.set_title(f\"{name}\", fontsize=16)\n",
    "    ax.grid(alpha=0.3)\n",
    "    ax.set_aspect(\"equal\", adjustable=\"box\")\n",
    "    ax.set_xlim([-0.5, 0.5])\n",
    "    ax.set_ylim([-0.5, 0.5])\n",
    "    ax.set_xlabel(r\"$\\theta_1$\", fontsize=15)\n",
    "    ax.set_ylabel(r\"$\\theta_2$\", fontsize=15)\n",
    "\n",
    "f.tight_layout()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "895b03ae",
   "metadata": {},
   "source": [
    "The posterior looks as we have expected in this case. However, in general, we do not know how the posterior is supposed to look like for any specific dataset. As such, we need diagnostics that validate the correctness of the inferred posterior. One such diagnostic is simulation-based calibration(SBC), which we can apply for free due to amortization. For more details on SBC and diagnostic plots, see:\n",
    "\n",
    "1. Talts, S., Betancourt, M., Simpson, D., Vehtari, A., & Gelman, A. (2018). Validating Bayesian inference algorithms with simulation-based calibration. *arXiv preprint*.\n",
    "2. Säilynoja, T., Bürkner, P. C., & Vehtari, A. (2022). Graphical test for discrete uniformity and its applications in goodness-of-fit evaluation and multiple sample comparison. *Statistics and Computing*.\n",
    "3. The practical SBC interpretation guide by Martin Modrák: https://hyunjimoon.github.io/SBC/articles/rank_visualizations.html\n",
    "\n",
    "Check out the next tutorial for a detailed walkthrough of the workflow's functionality."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "886bce5c",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "bf",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.11.11"
  },
  "toc": {
   "base_numbering": 1,
   "nav_menu": {},
   "number_sections": true,
   "sideBar": true,
   "skip_h1_title": true,
   "title_cell": "Table of Contents",
   "title_sidebar": "Contents",
   "toc_cell": true,
   "toc_position": {
    "height": "calc(100% - 180px)",
    "left": "10px",
    "top": "150px",
    "width": "165px"
   },
   "toc_section_display": true,
   "toc_window_display": true
  },
  "widgets": {
   "application/vnd.jupyter.widget-state+json": {
    "state": {},
    "version_major": 2,
    "version_minor": 0
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}