support_vector_machines.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Support Vector Machines\n",
    "\n",
    "- [1, What are SVMs?](#what-are-svms)\n",
    "- [2. Primal approach - Hard-Margin SVM](#hard-margin-svm)\n",
    "- [(Optional) - Deriving the margin requirement](#margin-derivation)\n",
    "- [3. Primal approach - Soft-Margin SVM](#soft-margin-svm)\n",
    "- [4. Solving the primal optimization problem](#solving-primal-svm)\n",
    "   - [4.1 Hinge loss function](#hinge-loss)\n",
    "   - [4.2 Updated objective function](#updated-objective-function)\n",
    "   - [(Optional) Three parts of the objective function](#two-parts)\n",
    "   - [4.3 Subgradient descent](#subgradient-descent)\n",
    "   - [(Optional) Subgradient descent vs. gradient descent](#subgradient-descent-vs-gradient-descent)\n",
    "- [5. Implementation of primal approach](#primal-implementation)\n",
    "   - [5.1 Toy dataset](#toy-dataset)\n",
    "   - [5.2 SVM class](#primal-svm-class)\n",
    "   - [5.3 Training and testing an SVM](#train-test-svm)  \n",
    "   - [5.4 Visualizing the decision boundary](#decision-boundary)\n",
    "- [6. Dual approach](#dual-approach)\n",
    "    - [6.1 Recap Lagrange multipliers](#recap-lagrange-multipliers)\n",
    "    - [6.2 Recap Lagrangian](#recap-lagrangian)\n",
    "    - [6.3 Dual optimization problem](#dual-optimization-problem)\n",
    "- [7. Primal vs. dual approach](#primal-vs-dual)\n",
    "- [8. Kernels / non-linear SMVs](#kernel-svms)\n",
    "    - [8.1 What is a kernel?](#what-is-a-kernel)\n",
    "    - [8.2 What are kernels good for?](#what-are-kernels-good-for)\n",
    "    - [8.3 Example](#kernel-example)\n",
    "    - [8.4 Can we also use kernels in the primal SVM?](#kernel-in-primal-svm)\n",
    "- [9. Sources and further reading](#sources)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Link to interactive demo\n",
    "\n",
    "[Click here](https://mybinder.org/v2/gh/zotroneneis/machine_learning_basics/HEAD?filepath=support_vector_machines.ipynb) to run the notebook online (using Binder) without installing jupyter or downloading the code.\n",
    "\n",
    "Sometimes, the GitHub version of the Jupyter notebook does not display the math formulas correctly. Please refer to the Binder version in case you think something might be off or missing.\n",
    "\n",
    "I also wrote a [blog post containing the contents of the notebook](https://alpopkes.com/posts/machine_learning/support_vector_machines/)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 1. What are support vector machines? <a class=\"anchor\" id=\"what-are-svms\"></a>\n",
    "\n",
    "Support vector machines (short: SVMs) are supervised machine learning models. They are the most prominent member of the class of [*kernel methods*](https://en.wikipedia.org/wiki/Kernel_method). SVMs can be used both for classification and regression. The original SVM proposed in 1963 is a simple binary linear classifier. What does this mean?\n",
    "\n",
    "Assume we are given a dataset $D = \\big \\{ \\mathbf{x}_n, y_n \\big \\}_{n=1}^N$, where $\\mathbf{x}_n \\in \\mathbb{R}^D$ and labels $y_n \\in \\{-1, +1 \\}$. A linear (hard-margin) SVM separates the two classes using a ($D-1$ dimensional) hyperplane.\n",
    "\n",
    "Special to SVMs is that they use not any hyperplane but the one that maximizes the distance between itself and the two sets of datapoints. Such a hyperplane is called *maximum-margin* hyperplane:\n",
    "\n",
    "<img src=\"figures/separating_hyperplanes.png\" width=\"800\"/>\n",
    "\n",
    "In case you have never heard the term margin: the margin describes the distance between the hyperplane and the closest examples in the dataset.\n",
    "\n",
    "Two types of SVMs exist: primals SVMs and dual SVMs. Although most research in the past looked into dual SVMs both can be used to perform non-linear classification. Therefore, we will look at both approaches and compare them in the end."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 2. Primal approach - Hard-margin SVM  <a class=\"anchor\" id=\"hard-margin-svm\"></a>\n",
    "\n",
    "When training an SVM our goal is to find the hyperplane that maximizes the margin between the two sets of points. This hyperplane is fully defined by the points closest to the margin, which are also called *support vectors*.\n",
    "\n",
    "The equation of a hyperplane is given by $\\langle \\mathbf{w}, \\mathbf{x} \\rangle + b = 0$, where $\\langle \\cdot, \\cdot \\rangle$ denotes the inner product and $\\mathbf {w}$ is the normal vector to the hyperplane. If an example $\\mathbf{x}_i$ lies on the right side of the hyperplane (that is, it has a positive label) we have $\\langle \\mathbf{w}, \\mathbf{x}_i \\rangle + b \\gt 0$. If instead $\\mathbf{x}_i$ lies on the left side (= negative label) we have $\\langle \\mathbf{w}, \\mathbf{x}_i \\rangle + b \\lt 0$.\n",
    "\n",
    "The support vectors lie exactly on the margin and the optimal separating hyperplane should have the same distance from all support vectors. In this sense the maximum margin hyperplane lies between two separating hyperplanes that are determined by the support vectors:\n",
    "\n",
    "<img src=\"figures/delimiting_hyperplanes.png\" width=\"450\"/>\n",
    "\n",
    "### Goal 1:\n",
    "When deriving a formal equation for the maximum margin hyperplane we assume that the two delimiting hyperplanes are given by:   \n",
    "$$\\langle \\mathbf{w}, \\mathbf{x}_{+} \\rangle + b = +1$$ \n",
    "$$\\langle \\mathbf{w}, \\mathbf{x}_{-} \\rangle + b = -1$$\n",
    "\n",
    "In other words: we want our datapoints two lie at least a distance of 1 away from the decision hyperplane into both directions. To be more precise: for our positive examples (those with label $y_n = +1$) we want the following to hold: $\\langle \\mathbf{w}, \\mathbf{x}_n \\rangle + b \\ge +1$.\n",
    "\n",
    "For our negative examples (those with label $y_n = -1$) we want the opposite: $\\langle \\mathbf{w}, \\mathbf{x}_n \\rangle + b \\le -1$. This can be combined into a single equation: $y_n(\\langle \\mathbf{w}, \\mathbf{x}_n \\rangle + b) \\ge 1$. This is our first goal: **We want a decision boundary that classifies our training examples correctly.**\n",
    "\n",
    "\n",
    "### Goal 2:\n",
    "Our second goal is to maximize the margin of this decision boundary. The margin is given by $\\frac{1}{\\mathbf{w}}$. If you would like to understand where this value is coming from take a look at the section \"*(Optional) Deriving the margin equation*\" below.\n",
    "\n",
    "Our goal to maximize the margin can be expressed as follows:\n",
    "$$ \\max_{\\mathbf{w}, b} \\frac{1}{\\Vert \\mathbf{w} \\Vert}$$\n",
    "\n",
    "Instead of maximizing $\\frac{1}{\\Vert \\mathbf{w} \\Vert}$ we can instead minimize $\\frac{1}{2} \\Vert \\mathbf{w} \\Vert^2$. This simplifies the computation of the gradient.\n",
    "\n",
    "### Combined goal\n",
    "Combining goal one and goal two yields the following objective function:   \n",
    "$$\n",
    "\\min_{\\mathbf{w}, b} \\frac{1}{2} \\Vert \\mathbf{w} \\Vert^2\n",
    "$$\n",
    "$$\n",
    "\\text{subject to: } y_n(\\langle \\mathbf{w}, \\mathbf{x}_n \\rangle + b) \\ge 1 \\text{ for all } n = 1, ..., N\n",
    "$$\n",
    "\n",
    "In words: we want to find the values for $\\mathbf{w}$ and $b$ that maximize the margin while classifying all training examples correctly. This approach is called the *hard-margin support vector machine*. \"Hard\" because it does not allow for violations of the margin requirement (= no points are allowed to be within the margin)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## (Optional) Deriving the margin equation <a class=\"anchor\" id=\"margin-derivation\"></a>\n",
    "\n",
    "We can derive the width of the margin in several ways (see sections 12.2.1-12.2.2 of the [Mathematics for Machine Learning book](https://mml-book.com)). Personally, I found the explanation of [this MIT lecture on SVMs](https://www.youtube.com/watch?v=_PwhiWxHK8o) easiest to understand.\n",
    "\n",
    "The derivation of the margin is based on the assumptions that we have already noted above:\n",
    "$$\\langle \\mathbf{w}, \\mathbf{x}_{+} \\rangle + b = +1$$ \n",
    "$$\\langle \\mathbf{w}, \\mathbf{x}_{-} \\rangle + b = -1$$\n",
    "\n",
    "Including the label of each example we can rewrite this as \n",
    "$$y_i (\\langle \\mathbf{w}, \\mathbf{x}_{i} \\rangle + b) -1 = 0$$\n",
    "\n",
    "Let's say we have a positive example $\\mathbf{x}_{+}$ that lies on the right delimiting hyperplane and a negative example $\\mathbf{x}_{-}$ that lies on the left delimiting hyperplane. The distance between these two vectors is given by ($\\mathbf{x}_{+} - \\mathbf{x}_{-})$. We want to compute the orthogonal projection of the vector onto the line that is perpendicular to the decision hyperplane. This would give us the width between the two delimiting hyperplanes. We can compute this by multiplying the vector ($\\mathbf{x}_{+} - \\mathbf{x}_{-})$ with a vector that is perpendicular to the hyperplane. We know that the vector $\\mathbf{w}$ is perpendicular to the decision hyperplane. So we can compute the margin by multiplying $(\\mathbf{x}_{+} - \\mathbf{x}_{-})$ with the vector $\\mathbf{w}$ where the latter is divided by the scale $||\\mathbf{w}||$ to make it a unit vector.\n",
    "\n",
    "<img src=\"figures/maximum_margin_derivation.png\" width=\"500\"/>\n",
    "\n",
    "$$\n",
    "\\begin{align*}\n",
    "\\text{width} &= (\\mathbf{x}_{+} - \\mathbf{x}_{-}) \\cdot \\frac{\\mathbf{w}}{||\\mathbf{w}||} \\\\\n",
    "&= \\frac{\\mathbf{x}_{+} \\cdot \\mathbf{w}}{||\\mathbf{w}||} - \\frac{\\mathbf{x}_{-} \\cdot \\mathbf{w}}{||\\mathbf{w}||}\n",
    "\\end{align*}\n",
    "$$\n",
    "\n",
    "For the positive example $\\mathbf{x}_{+}$ we have $y_+ = +1$ and therefore $(\\langle \\mathbf{w}, \\mathbf{x}_{+} \\rangle = 1 - b$. For the negative example $\\mathbf{x}_{-}$ we have $y_- = -1$ and therefore $- (\\langle \\mathbf{w}, \\mathbf{x}_{-} \\rangle) = 1 + b$:\n",
    "\n",
    "$$\n",
    "\\begin{align*}\n",
    "\\text{width} &= (\\mathbf{x}_{+} - \\mathbf{x}_{-}) \\cdot \\frac{\\mathbf{w}}{||\\mathbf{w}||} \\\\\n",
    "&= \\frac{\\mathbf{x}_{+} \\cdot \\mathbf{w}}{||\\mathbf{w}||} - \\frac{\\mathbf{x}_{-} \\cdot \\mathbf{w}}{||\\mathbf{w}||} \\\\ \n",
    "&= \\frac{(1 - b) + (1 + b)}{||\\mathbf{w}||}\\\\ \n",
    "&= \\frac{2}{||\\mathbf{w}||}\\\\ \n",
    "\\end{align*}\n",
    "$$\n",
    "\n",
    "We conclude that the width between the two delimiting hyperplanes equals $\\frac{2}{\\mathbf{w}}$. And therefore, that the distance between the decision hyperplane and each delimiting hyperplane is $\\frac{1}{\\mathbf{w}}$."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "##  3. Primal approach - Soft-margin SVM <a class=\"anchor\" id=\"soft-margin-svm\"></a>\n",
    "\n",
    "In most real-world situations the available data is not linearly separable. Even if it is, we might prefer a solution which separates the data well while ignoring some noisy examples and outliers. This motivated an extension of the original hard-margin SVM called *soft-margin SVM*.\n",
    "\n",
    "A soft-margin SVM allows for violations of the margin requirement (= classification errors). In other words: not all training examples need to be perfectly classified. They might fall within the margin or even lie on the wrong side of the decision hyperplane. However, such violations are not for free. We pay a cost for each violation, where the value of the cost depends on how far the example is from meeting the margin requirement.\n",
    "\n",
    "To implement this we introduce so called *slack variables* $\\xi_n$. Each training example $(\\mathbf{x}_n, y_n)$ is assigned a slack variable $\\xi_n \\ge 0$. The slack variable allows this example to be within the margin or even on the wrong side of the decision hyperplane:\n",
    "\n",
    "- If $\\xi_n = 0$ the training example $(\\mathbf{x}_n, y_n)$ lies exactly on the margin\n",
    "- $0 \\lt \\xi_n \\lt 1$ the training example lies within the margin but on the correct side of the decision hyperplane\n",
    "- $\\xi_n \\ge 1$ the training example lies on the wrong side of the decision hyperplane\n",
    "    \n",
    "We extend our objective function to include the slack variables as follows:\n",
    "$$ \\min_{\\mathbf{w}, b, \\mathbf{\\xi}} \\frac{1}{2} \\Vert \\mathbf{w} \\Vert^2 + C \\sum_{n=1}^N \\xi_n $$\n",
    "\n",
    "$$ \\text{subject to:} $$\n",
    "\n",
    "$$ \\begin{equation}\n",
    "y_n(\\langle \\mathbf{w}, \\mathbf{x}_n \\rangle + b) \\ge 1 - \\xi_n\n",
    "\\end{equation}$$\n",
    "\n",
    "$$ \\xi_i \\ge 0 \\text{ for all } n = 1, ..., N $$\n",
    "\n",
    "Note: the objective function is somewhat not displayed correctly within the GitHub version of the notebook. It should look as follows:\n",
    "\n",
    "<img src=\"figures/primal_optimization_problem.png\" width=\"220\"/>\n",
    "\n",
    "The parameter $C$ is a regularization term that controls the trade-off between maximizing the margin and minimizing the training error (which in turn means classifying all training examples correctly). If the value of $C$ is small, we care more about maximizing the margin than classifying all points correctly. If the value of $C$ is large, we care more about classifying all points correctly than maximizing the margin."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 4. Solving the primal optimization problem <a class=\"anchor\" id=\"solving-primal-svm\"></a>\n",
    "\n",
    "Theoretically, the primal SVM can be solved in multiple ways. The most well known way is to use the [hinge loss function](https://en.wikipedia.org/wiki/Hinge_loss) together with [subgradient descent](https://en.wikipedia.org/wiki/Subgradient_method).\n",
    "\n",
    "### 4.1 Hinge loss function <a class=\"anchor\" id=\"hinge-loss\"></a>\n",
    "The hinge loss function given the true target $y \\in \\{-1, +1\\}$ and the prediction $f(\\mathbf{x}) = \\langle\\mathbf{w}, \\mathbf{x}\\rangle+b$ is computed as follows:\n",
    "\n",
    "$$\\ell(t)=\\max \\{0,1-t\\} \\quad \\text{where} \\quad t=y \\cdot f(\\mathbf{x})= y \\cdot \\big(\\langle\\mathbf{w}, \\mathbf{x}\\rangle+b\\big)$$\n",
    "\n",
    "Let's understand the output of this loss function with a few examples:\n",
    "- If a training example has label $y = -1$ and the prediction is on the correct side of the marghin (that is, $f(\\mathbf{x}) \\le -1$), the value of $t$ is larger or equal to $+1$. Therefore, the hinge loss will be zero ($\\ell(t) = 0$)\n",
    "- The same holds if a training example has label $y = 1$ and the prediction is on the correct side of the margin (that is, $f(\\mathbf{x}) \\ge 1$)\n",
    "- If a training example ($y = 1$) is on the correct side of the decision hyperplane but lies within the margin (that is, $0 \\lt f(\\mathbf{x}) \\lt 1$) the hinge loss will output a positive value.\n",
    "- If a training example ($y = 1$) is on the wrong side of the decision hyperplane (that is, $f(\\mathbf{x}) \\lt 0$), the hinge loss returns an even larger value. This value increases linearly with the distance from the decision hyperplane\n",
    "\n",
    "<img src=\"figures/hinge_loss.png\" width=\"600\"/>\n",
    "\n",
    "\n",
    "### 4.2 Updated objective function <a class=\"anchor\" id=\"updated-objective-function\"></a>\n",
    "\n",
    "Using the hinge loss we can reformulate the optimization problem of the primal soft-margin SVM. Given a dataset $D = \\big \\{ \\mathbf{x}_n, y_n \\big \\}_{n=1}^N$ we would like to minimize the total loss which is now given by:\n",
    "\n",
    "$$\n",
    "\\min _{\\mathbf{w}, b} \\frac{1}{2}\\|\\mathbf{w}\\|^{2} + C \\sum_{n=1}^{N} \\max \\left\\{0,1-y_{n}\\left(\\left\\langle\\mathbf{w}, \\mathbf{x}_{n}\\right\\rangle+b\\right)\\right\\}\n",
    "$$\n",
    "\n",
    "If you would like to understand why this is equivalent to our previous formulation of the soft-margin SVM please take a look at chapter 12.2.5 of the [Mathematics for Machine Learning book](https://mml-book.com)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### [Optional] Three parts of the objective function <a class=\"anchor\" id=\"two-parts\"></a>\n",
    "\n",
    "Our objective function can be divided into three distinct parts:\n",
    "\n",
    "Part 1: $\\frac{1}{2}\\|\\mathbf{w}\\|^{2}$\n",
    "\n",
    "This part is also called the *regularization term*. It expresses a preference for solutions that separate the datapoints well, thereby maximizing the margin. In theory, we could replace this term by a different regularization term that expresses a different preference.\n",
    "\n",
    "Part 2: $\\sum_{n=1}^{N} \\max \\left\\{0,1-y_{n}\\left(\\left\\langle\\mathbf{w}, \\mathbf{x}_{n}\\right\\rangle+b\\right)\\right\\}$\n",
    "\n",
    "This part is also called the *empirical loss*. In our case it's the hinge loss which penalizes solutions that make mistakes when classifying the training examples. In theory, this term could be replaced with another loss function that expresses a different preference.\n",
    "\n",
    "Part 3: The hyperparameter $C$ that controls the tradeoff between a large margin and a small hinge loss."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 4.3 Sub-gradient descent <a class=\"anchor\" id=\"subgradient-descent\"></a>\n",
    "\n",
    "The hinge loss function is not differentiable (namely at the point $t=1$). Therefore, we cannot compute the gradient right away. However, we can use a method called [subgradient descent](https://en.wikipedia.org/wiki/Subgradient_method) to solve our optimization problem. To simplify the derivation we will adapt two things:\n",
    "1. We assume that the bias $b$ is contained in our weight vector as the first entry $w_0$, that is $\\mathbf{w} = [b, w_1, ..., w_D]$\n",
    "2. We divide the hinge loss by the number of samples\n",
    "\n",
    "Our cost function is then given by \n",
    "$$\n",
    "J(\\mathbf{w}) = \\frac{1}{2}\\|\\mathbf{w}\\|^{2} + C \\frac{1}{N} \\sum_{n=1}^{N} \\max \\left\\{0,1-y_{n}\\left(\\left\\langle\\mathbf{w}, \\mathbf{x}_{n}\\right\\rangle\\right)\\right\\}\n",
    "$$\n",
    "\n",
    "We will reformulate this to simplify computing the gradient:\n",
    "$$\n",
    "J(\\mathbf{w}) = \\frac{1}{N} \\sum_{n=1}^{N} \\Big[ \\frac{1}{2}\\|\\mathbf{w}\\|^{2} + C \\max \\left\\{0,1-y_{n}\\left(\\left\\langle\\mathbf{w}, \\mathbf{x}_{n}\\right\\rangle\\right)\\right\\}\\Big]\n",
    "$$\n",
    "\n",
    "\n",
    "The gradient is given by:\n",
    "$$\n",
    "\\nabla_{w} J(\\mathbf{w}) = \\frac{1}{N} \\sum_{n=1}^N \\left\\{\\begin{array}{ll}\n",
    "\\mathbf{w} & \\text{if} \\max \\left(0,1-y_{n} \\left(\\langle \\mathbf{w}, \\mathbf{x}_{n} \\rangle \\right)\\right)=0 \\\\\n",
    "\\mathbf{w}-C  y_{n} \\mathbf{x}_{n} & \\text { otherwise }\n",
    "\\end{array}\\right.\n",
    "$$\n",
    "\n",
    "With this formula we can apply stochastic gradient descent to solve the optimization problem."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### [Optional] Difference subgradient descent and gradient descent <a class=\"anchor\" id=\"subgradient-descent=vs-gradient-descent\"></a>\n",
    "\n",
    "The subgradient method allows us to minimize a non-differentiable convex function. Although looking similar to gradient descent the method has several important differences.\n",
    "\n",
    "#### What is a subgradient?\n",
    "\n",
    "A subgradient can be described as a generalization of gradients to non-differentiable functions. Informally, a sub-tangent at a point is any line that lies below the function at the point. The subgradient is the slope of this line. Formally, the subgradient of a convex function $f$ at $w_0$ is defined as all vectors $g$ such that for any other point $w$\n",
    "\n",
    "$$ f(w) - f(w_0) \\ge g \\cdot (w - w_0) $$\n",
    "\n",
    "If $f$ is differentiable at $w_0$, the subgradient contains only one vector which equals the gradient $\\nabla f(w_0)$. If, however, $f$ is not differentiable, there may be several values for $g$ that satisfy this inequality. This is illustrated in the figure below.\n",
    "\n",
    "<img src=\"figures/gradient_vs_subgradient.png\" width=\"900\"/>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Subgradient method\n",
    "\n",
    "To minimize the objective function $f$ the subgradient method uses the following update formula for iteration $k+1$:\n",
    "\n",
    "$$ w^{(k+1)} = w^{(k)} - \\alpha_k g^{(k)}$$ \n",
    "\n",
    "Where $g^{(k)}$ is *any* subgradient of $f$ at $w^{(k)}$ and $\\alpha_k$ is the $k$-th step size. Thus, at each iteration, we make a step into the direction of the negative subgradient. When $f$ is differentiable, $g^{(k)}$ equals the gradient $\\nabla f(x^{(k)})$ and the method reduces to the standard gradient descent method.\n",
    "\n",
    "More details on the subgradient method can be found [here](https://web.stanford.edu/class/ee392o/subgrad_method.pdf)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 5. Implementation primal approach <a class=\"anchor\" id=\"primal-implementation\"></a>\n",
    "\n",
    "### 5.1 Toy dataset <a class=\"anchor\" id=\"toy-dataset\"></a>\n",
    "To implement what we have learned about primal SVMs we first have to generate a dataset. In the cell below we create a simple dataset with two features and labels +1 and -1. We further split the dataset into a test and train set."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAe4AAAGDCAYAAADtffPSAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Z1A+gAAAACXBIWXMAAAsTAAALEwEAmpwYAACPm0lEQVR4nOzddZxUZRfA8d9zpzdh6Q4REBFREQMssDAwMbBFsbsVxUKxuxADVERfMTARRbAwQEBBBUFE6Vq2puc+7x93WHZ27uzOLhssnO/n4+dlbz4DvJx56hyltUYIIYQQDYNR3w0QQgghRPokcAshhBANiARuIYQQogGRwC2EEEI0IBK4hRBCiAZEArcQQgjRgEjgFkKglNJKqS713Q4hROUkcAuxjVFKFZf5z1RKBcr8fEY9t61jPMg7t4f3CNEQyf8phNjGaK2zNv9aKfUPcIHW+ov6a5EQYlsiPW4hGgillEcp9bhSamX8v8eVUp74uflKqWPLXOtSSq1XSu2R4lk3KKVWxZ9zfrlzRyul5iilCpVS/yml7ixz+uv4/26KjwDsp5TaSSk1TSm1If7ON5RSjco87yal1AqlVJFSaqFSamD8uKGUulkptSR+79tKqbxU79ma3zshticSuIVoOG4D9gV6A7sDfYER8XPjgTPLXHsUsEprPaf8Q5RSRwLXA4cBOwOHlrukBDgbaAQcDVyilDo+fu7A+P820lpnaa1nAgq4H2gN7AK0A+6Mv6sbcDmwt9Y6GzgC+Cf+jCuA44GD4vfmA89U8B4hBBK4hWhIzgDu1lqv1VqvA+4Czoqfex04SimVE//5LOC1FM85BXhFaz1fa11CPMhuprWerrX+TWttaq1/Bd7ECq62tNaLtdZTtdaheLseLXN9DPAAPZRSLq31P1rrJfFzFwO3aa2Xa61D8XacLPPaQlRMArcQDUdrYFmZn5fFj6G1Xgl8B5wUH6YeBLxRwXP+K/ecUkqpfZRSXyml1imlCrACbNNUjVJKtVBKTYwPhxdifYloGm/XYuBqrKC8Nn5d6/itHYD3lFKblFKbgD+wAn2Lin4ThNjRSeAWouFYiRXsNmsfP7bZOKzh8iHATK31ihTPWYU1nF32OWVNACYD7bTWucDzWMPhAHblBO+LH99Na50Tb8Pm69FaT9Ba94+3XQMPxE/9BwzSWjcq85833m4pWyhEChK4hWg43gRGKKWaKaWaAndg9W43ex/YE7gKa847lbeBc5VSPZRSGcDIcuezgY1a66BSqi8wtMy5dYAJdC53fTFQoJRqA9yw+YRSqptSakB8EV0QCMTvB+sLwSilVIf4tc2UUsdV8B4hBBK4hWhI7gVmAb8CvwG/xI8BoLUOAJOATsC7qR6itf4UeByYBiyO/29ZlwJ3K6WKsL4cvF3mXj8wCvguPsS9L9Zc+55AAfBxuXd7gNHAemA10By4JX7uCaye/efxd/0A7FPBe4QQgNJaRqSE2F4ope4Aumqtz6z0YiFEgySrN4XYTsT3QA9jy0pzIcR2SIbKhdgOKKUuxFrs9anW+uvKrhdCNFwyVC6EEEI0INLjFkIIIRoQCdxCCCFEA9IgFqc1bdpUd+zYsb6bIYQQQtSJ2bNnr9daN7M71yACd8eOHZk1a1Z9N0MIIYSoE0qpZanOyVC5EEII0YBI4BZCCCEaEAncQgghRAMigVsIIYRoQCRwCyGEEA2IBG4hhBCiAZHALYQQQjQgEriFEEKIBkQCt9guxKIxTNOs72YIIUStk8AtGrTFc5dy+b63MMh7OsdknskjFzxHoDhQ380SQoha0yBSngphZ93yDVx70EgCRVagjoQifPnGN6xauoaHv7yzfhsnhBC1RHrcosH64JnPiIQiCccioQh//vgX/yz4r55aJYQQtUsCt2iwlv66jGg4mnTc4XSwfNHKemiREELUPgncosHq1rcLbq8r6Xg0HKVjz/b10CIhhKh9ErhFg3XsxYfj9rlRSpUec3td7HXY7rTduVU9tkwIIWqPBG7RYDVu0YinfrifvkftgdvrIjsvixOuOpoRb19b300TQohaI6vKRYPWdudW3PvhLfXdDCGEqDPS4xZCCCEaEAncQgghRAMigVsIIYRoQCRwCyGEEA2IBG4hhBCiAZHALYQQQjQgEriFEEKIBkQCtxBCCNGASOAWQgghGhAJ3EIIIUQDIoFbCCGEaEBqLXArpV5WSq1VSs0vc+whpdSfSqlflVLvKaUa1db7hRBCiO1Rbfa4XwWOLHdsKtBTa90LWARIdQghhBCiCmotcGutvwY2ljv2udY6Gv/xB6Btbb1fCCGE2B7V5xz3+cCnqU4qpYYrpWYppWatW7euDpslhBBCbLvqJXArpW4DosAbqa7RWo/RWvfRWvdp1qxZ3TVOCCGE2IY56/qFSqlzgWOAgVprXdfvF3Vr1d9reP+pT1n2x3J23b8bx15yOI2a5dZ3s4QQosGq08CtlDoSuBE4SGvtr8t3i7q34PuF3HzEPUTCUWKRGL99/TvvP/UJz/z8AC07Nq/v5gkhRINUm9vB3gRmAt2UUsuVUsOAp4FsYKpSaq5S6vnaer+of48Nf55gSYhYJAZAOBihOL+EsTennCERQghRiVrrcWutT7c5/FJtvU9sW0oK/SxftCrpuGlqZk2ZW/cNEkKI7YRkThO1wuVxoQxle86X7avj1gghxPZDAreoFW6Pi/4n9MXlSRzU8WS4Oe7SI2zvCQfDLP9rFSWFsvxBCCFSqfNV5WLHcfXzw1m/YiN//fI3DpeTaCjC/sf1Zcj1gxOu01rz9sOTef2ed1BANBrjsLMO4oqnh+F0yV9RIYQoS/5VFLUmMzeTx76+h6Xz/2X10rV02q297WryL9/4htfu+h8hf2jLsde/xu11cdkT59dlk4UQYpsnQ+Wi1nXq2Z79ju2TcgvYhPveTQjaAKFAmE/GfkkkHKmLJgohRIMhgVvUu/zVm2yPmzGTQFGwWs80TXMrWiSEENsuCdyi3nXbeyfb4zl5WWQ1zqzSs754fQant7+II5yncmqb4Xz68pc10UQhhNhmSOAW9e6CB87Em+lBqS3bxzwZbi55/FwMI/2/ol+8PoNHL3yB9cutonQbV+XzzJUvS/AWQmxXJHCLeteldyeenHkf/Y7vS/P2Tel1UA/umXwzB5/SL+1n/PnTXzx07jNEQolz4iF/mPEj367pJgshRL2RVeWiQmuWrSMcDNO2a+uEHnFN69SzPSMnXV+te4P+EDcfcS+maV+zZv2KjZimWaXeuxBCbKskcAtbKxav4u6TH2H5opUYDoOsRpnc/PqV7H7QrnXelvw1m3jz/vf48ZNfyGmSzcnXHMOBQ/Yr/SLx40ezK1yM1qxdUwnaQojthgRukSQWjXHdwSPZuHoTOt6LDZaEGHHM/bzy5xM0bdOkztpSuKGIi/a4gaINRUQjMVYuXs3Dw57l7/n/ct7dpwFQlF+CGbMP3A6nwfn32aXNF0KIhkm6ISLJrClz8RcFS4P2ZrFojM9e+apO2/LBM59RsqmEaLzCGFhfIt55eDKFG4sA6D2gZ1JbAVBw0rXHcugZB9ZVc4UQotZJ4BZJNqzaZNuDjYSirP13fZ22ZfbUeYSDyUlYXB4Xi+f8A0DbnVtx1AUD8WZ6Ss97Mz3sObAXw+4bavvckkI/87/9gxWLkyuYCSHEtkyGykWSHvt1BZ3cg/Vmeel98JY5btM0UUrV6qK1Fh2a8/vMRUk96nAgTKNmOaU/X/rE+ex1eG8+efELwqEIA884gAGn9y+d254y7ivGj3ybDSvzyWqcSUmBH4/PTTQcpcsenbjr/RvJbZqDEEJs65S2+Qd6W9OnTx89a9as+m5Gg1OwvpCp42ewYvFqevbrzgEn74vb40rr3lGnP8bMD2eXpiJ1e1207tKSZ2c9wMola3jy0hf57Zs/cLocDDi9P5c8fh6ZORk1/hkWzV7CtQfdQcgfTjiuDEXbrq158vtRZDWqOEnLh89P4YXrX0tKq7qZ0+Vg137deXjanTXVbCGE2CpKqdla6z625yRwb58Wz1nKdYeMJBqJEQ6E8WZ5adKqMU/9cB/ZjbMqvT8Wi/Hp2Gl89MLnhIMRDjmtHydfewyhQJjzul1FSUFJaafc5XHSZY/OPPHdvbXS+57xv5ncN/TxpOF7l9vJ0cMP47InUxci0VozpMUwCtYXVfgOt9fFuL+eqtOFd0IIkUpFgVvmuLdT95/5JP7CAOGA1VMNFgdZs2wdr9/zTlr3OxwOjrnoMJ7/5SFe/v1xzrpjCL4sH5+M/ZJwMJwwkh4JRVn62zIWzf67Nj4KfQ7vhd33gUg4yvS3v6/w3mBJkOJNJZW+w+FyULihuLpNrLZYLMabo9/ltLbDOTbnLG4/7gGWL1pZ5+0QQjQcEri3Q/lrNrHq79VJx6PhKDP+N3Ornv333H9sF4sppVi+sPoBZ+WS1fw+cyGB4kDysw0jZU/ecFTcw/dkeMhIYwjfMAzadW+dXmNr0KMXPM8b905iw8p8gsVBfvxoNpfvcwvrV2yo87YIIRoGCdzbIYfLYbe2DLDmc7dG1z474fa5k46bpqbDrm2r/LxN6wq4qt9tXNjrOm4ZNIohLS9k0uMfJVyTke1jl/26YjgS/7q6vS4OPfOgCp9vGAZn3n4SngyP7XmltuRFd7nTm/+vKetXbOCrid8lzN9rrQkFwrz7xCd12hYhRMMhgXs7lJOXTbe9d0oKdB6fm6MuGLhVzz5y2AC8GW6UsaWn6/a66N63C116d6ry8+466WEWzVpCOBDGXxgg5A/xyoiJ/PzZnITrbhp3BXktG5GR7cPpcuDL8tJ5946cecfJlb7jhCuP5oLRZ9CoeS4ATdvm0e/4vnTdqzP7H9+X0VNu54hzDqly27fWst+X4/Ymf1mIhqP88cOiOm+PEKJhkO1g26lb37iKaw68g6L8EmLRGIah6LFfN4ZcP7jC+5bM+4eZH87C5XZx4JB9adWpRcL5nLxsnv5xNM9c9TK/fPErLo+Lw885mGH3n1HlNs76fC6/f78wKcd4yB/inUc/Yu8j9yg91qJDM8YveZrvP5jF/G//oM3OrTjmosNwuqy/wuuWb8BwGDRp1TjpPUopjr98EMdfPmibylneqnMLIuFo0nGH06Djru3qoUVCiIZAVpVvx2LRGLOmzGXNsvV023snuu3dpcLrx9z4GpOf/YxIKIrhMDAMxeVPDWPQsK3rpdtZPGcpV/W7zXa+HKBjz3a8+OujCccWfL+QO098iFDA2tbldDkZdv9QJj32EWv+WYfW0L57G0a8dQ1tu9b9fHV13HrUKOZOX0CkzO+DN8PDc7882GA+gxCi5sl2MFGpP378ixsG3pm0X9rtdfH6P8/ROD7MXFNuOvwefvniV9tzTreTE648iuEPnlV6bNqEbxl99pP2qU3LadQshzf+fT7tPev1KVAS5JkrX2bahG8xYzHadmvD1c8Pp2e/7vXdNCFEPZLtYKJSM97+nnAgufdrOAx+/Gh2jb9v4c+LU55TCk65YcuQ/hdvfM1D5z+dVtAGKC70M/ODn7e6jXXBl+nl+pcu5cOi13gvfxxjf3tUgrYQokIyx92AxKIx5n/3J9FwlF37dccbXym99r/1THn1KzasyGfPQ3ej3/F9cTirtnrccCiUss10mrAQraY0ap5LSYHf9lwkHGXV32tp1CwXrTUv3vAa0XDM9lo70VCUV++YyLwZCzj8nIPp3ndnIuEIc76cT8H6QgBym+bQe0DPbaZX7nA68FXxz0wIsWOSwN1A/P7DIm4fPJpofDGTGTO54ZXLyMzNYOQJD2HGYkRCUb6c8A1vP/QBj0y/C4/PfguUnUNO68/kZ6ckDZWbMZP9jrUdrdkqp9wwmMeGv2B/UsPU8dPZZZ+d2bgqn01rC6r8/OWLVrFi8Wo+Hzedw885mOlvfU8oEC5NSOPyOHF73dz1/o31UmNcCCGqS+a4G4CgP8RpbYYn9VDdXhdurzspM5gnw82595zGydccW6X3jL/rbd568AN0zEQ5DNBw3UsXM+D0A7b6M5RnmibHZJ5JJGS/OG3QBQPpuldnnr92HKFA2PaamuDL8vLWyjH4sny19g4hhKgqmeNu4H78aDammVxmMxqJ2qbzDPnDTHvjmyq/5+yRp/DC3Ic5/76hDH/wLF77++kaD9qmafLZK19xed9byMi1D5Yuj5OdenesMGgrQzHogoEpE6tUxczJO+6XQiFEwyND5Q1AUX6JbX1sM5Z6tMTlTc5ulo62O7fi5Gur1lOvikcveI4Z/5tJsMTa0lV+Xt3lcTJw6AEs/Hkx4RS98YwcHxc+eBbHDD+Mjj3b8cptb5Y+r6pM08RfFKzWvUIIUR8kcDcAewzsmfaK6s123X/bW5m8YvEqvpr4XcLeba3B7XPT54jdade1NfsN3pse+3Xl9sGjbT+zL9vLjeMup99xfQE48cqjadOlFW/e9y7L/lhOcX7lBUXK0qZmr8N7bd0Ha+BM0+Sv2X8TCUXo1rdLnad+FUJUjQyVNwBturTi2IsPx5u5ZVjYm+khIyf1vOyHz33GPwv+q4vmobXmly9+5dHhz/P0lS+xcNaShPMlhX5eHfkW1xx4h22msHAgjDfDwwWjz2TX/buhlKLf8X1th8GjkRi79d8l4ViXPTqx+p+1REKJz/ZmePBlezn2ksPxZCSPQHh8bk68+uik7HA7ksVzl3JGh0u4YeBd3HbM/ZzcfBgzP5SpAyG2ZbI4rYHQWjNrylw+GfsFoUCEgUMPIBqO8tQVY5NWgoM1BN3nyD247+Nba71dD573DNMnfle64t0wFMdeeiSXP3k+4WCYS/a6kVVL1yZkByvL6XYw5LrBnD9qaOmxcDDM1f1H8N+fKwn6QygFbp+Hs+88hVPKpW194Jyn+OrN74hFE7eMtercgjG/PoI3w8PCWUv45MWpLPt9OdrUtOrcgkHDBrL7wdaK8qXz/2XcHW/x50+LadmpGWeMOJm9j+hdg79TlQv6Q/z86RyCJSH2PKyXbfrWmhQOhjmt7UUUbUwsZ+rJcPPib4/u0F9ohKhvFS1Ok6HyBkIpxd5H7pGQv1trze8zF/LJ2C+TrtcaFny3sNbb9ds3f/DVm98Si2wJmqap+eDpTzngpH1Yu2w9a/9dnzJogzVXv+ehvcrcb+LyuHj823v5fNwMvpk0k+wm2Qy+5AhadmrOyiWradW5RWmpzx8+mp0UtAHW/ru+dG1Atz470a3PTrbv//vXZVzV7zZC/jBaazas3MhdJz3EVc8N57CzKq4+VlPmf/sHtx1zP2gwtcaMxjhr5BBOu+mEWnvnj5/MSfhz2ywWjfH5q9M5565Ta+3dQojqq7XArZR6GTgGWKu17hk/lge8BXQE/gFO0Vrn11YbtndKKS5/ehhTX/vadltVbtOsWm/DjLe/t/3HH+CV294kp0lWpQvHzJjJbcfcz9XPXcikxz/m73nL8GR4OObiwzh/1Okcc9FhrFi8inuGPMp/C1egDEWjZrncOuEqeuzXDVeKJCpKpVfGdOwtrye1MeQP88J14xl4xgFJRUlM0+Szl79i8rOfESwJceDJ+3LKDceR1Siz0nfZCQfDjDh2NP7CxFrkr9/zDr0O2pUe+3at1nMrU7SxmJjNosdoOFatvfNCiLpRm3PcrwJHljt2M/Cl1npn4Mv4z2IruNwuDj/3INy+xODlyfBwyg3HEYvF+N8jkxna/mIG557N7YNH89/CFTX2/lQrvwFW/r2GWZ/PS+85gTAPnvsMS+b+g9aaYEmQD5+dwsPnP0s0EuXag0by96/LCAcjhPxh1ixbx81H3Ev+mk0cNWxAUo1wh8vB3oP2wF3J6vrfZy5k1mf2bQwUB2wD2GPDX+DZq19hydx/WPHXKt559EMu3+cWgv4twT8UCLFyyeqEY6n88sVv2E1ZhQMRprzyVaX3V1evg3qgbbYZerO89KnjaQIhRPpqrcettf5aKdWx3OHjgIPjvx4HTAduqq027Cgufew8SgoCfPfeT7g8TqLhKMdfMYijhx/GE5eM4YvXvyEUDyA/fvwLv379Oy/+9ijN2zWt8ruW/bGcb9/9kWg4ijfLg8emnvRm2Y0zKaniKu+yQoEw3777I3se1otgcTApuMWiMd4c/R5rl60DrLzqTrcTh8OgeYdmXDf2kgqfX7ixiJuPvNc2aAKgVFIveuWS1Uyb8E3CyvhIKMqGlRuZNuFbBg0bwLg73+KdRz6ytrqZmsGXHckFo89IWU40FAiDTRM2f4GpLW13bsUR5w9g6rjppSMOngwPO+/RiX2P3avW3iuE2Dp1PcfdQmu9Kv7r1UDK1S9KqeHAcID27dvXQdOqZvlfqxh3x0R+nfE7jVs24vSbT+CgU/avl7a4vW5um3A1+WsLWL98A627tCQzJ4P8NZv4fNyMhGF0rTXhQIRJj33EJY+em/CckkI/wZIQeS0blc4flzXhvkm8MepdouFo6dyxYShQ2AaezNzMCnvk6XC6nUx9bQb+okDSuXAwwuRnPsOM6dLga8ZiDLt/KCdccVTCZ9i4Op8Pnv6MP39eTKfd2nPCFUfx48e/2O6PB6vHvv9xeycNw//54184XA4oN2cfLAnxyxe/EiwJ8s4jH5V+UQKY/OwUMnMzOOO2k2zftcfAnkQjyavtvZkeDjx5vxS/MzXjiqeGsefA3fh4zBeEAiEGDj2Aw889GIdD8qYLsa2q1VXl8R73R2XmuDdprRuVOZ+vta506ey2tqp81d9ruHjPGwgWBzHje409mR7OHHFSrS4mqqq5X83nzhMfsi3m0X2fnXlq5n0AFOUX89C5z/DzlLkoQ9G4RSOuG3sJew7crfT6/xau4OI9bkhZP9uO2+vCjJlEy8+Bpwj0tqpybVxey0acfsuJ9D9pH5q2zmP5opVcse+thAIhIqEoTpcTl8fJQafuz2cvTUv5XpfHSZOWeYx893q69O4EwLzpC7j9uNEEyiVtMRwGrbu0ZPXStaWr68vKapzJexteTdnmj174nOevHUck/qXIm+VljwE9ufPdG1L21IUQ269tKeXpGqVUK4D4/66t4/fXiAn3TSJYEioN2gChkhCv3zMprTnNutKqcwvbRWuGw6Bjj7alP9929P38PGUu0XCUSDDC2mXrGHncAwlz4d+99xOxqH3vNJVwMJIUtA2HQbvurTEc5f7qpSpAVo3vlRtXb+LFm1/nnC6X88XrM3j+unGUFPhL93lHI1ECxUF++/oPvFnelO+NBKOs/mct1w+4s/TPdbcDd6FRs9yk9psxk+ULV9oGbYDi/BLbtLWbHXPR4Tz14/2ccNVRDBo2gNvfulaCthDCVl3/qzAZOCf+63OAD+r4/TXit2//tB1idTgNVi5eXQ8tste8fVM69eqA4UiMig6nwSFD+wPW/uW/f/0nKeCEg2Hee+KT0p9XLFlNLJZ+ac1UzJjJun83cOzFh+H2uvBmevD43LTu3DLl6vDqCAfChIMRHhv+AnOmzbedx16xeBXtu7dJWthWXrAkxPfv/wSAYRg8PG0kO+/ZCbfXhaeSezdr1711pUG4U8/2XPzwOVz74iX0HbSHBG0hhK3a3A72JtZCtKZKqeXASGA08LZSahiwDDiltt5fm1p1as6KRauSjkdCUfJaNar7BtkwTZN7TnmUZQv+S8ppbjgMRhw7moOG7McBJ++b4n7Nb9/9CVjDuNPe+LZavV87wZIQ3773ExNXjuHvucvIbZbDhPveZeXfNf+lRxkqdbs1/PvncqLhKNl5WUmJSDaLRWLkr9myurx5+2Y8/eNo1i3fwNfvzGTcyLeShs7L8vjcXPLYeZW2tWB9IW+Ofo+ZH8wiMzeDE648ikPPOtB2vYEQYsdVm6vKT09xamBtvbOunHbzCfz2zR8JGcvcXhf7HL0XjZrl1mPLLKFAiHce+ZCfPp1TWn864Xy83d9M+oHm7ZulrMC1YcVGwqEIL974OuFgzZbW3LhqE6/f/Q4XjD4Dl9tlVTmrheUWdlnlygoWh+LXVTzF0fOAXZKONWvbhA492qFSjPN7Mj3sun83zrpjCD37VZw7vqTQzyV73Uj+moLS0Y8nL3uRRbOXcNkT51d4rxBixyJjcdWw+0G7cu2LF5PbNBuPz43L46L/iftw47jLa/3ds6fO48r9b2NIy2HceNjd/PHjXwnnv5n0A0NaXMBrd79jG7TLCvnDfPHajJQBs2hjMSsXr0bXQkTVWvPR858z4pjRaK1p06XlVj0vac68isLBiNU7t9G4RW7KrGveTI/tynmPz80T397LA1NurzRoA3zy4hcUri9KmLIIloT4eMwXrF+5Mc1PIYTYEUjK02oacPoBHHTK/mxYsZGsxllkZKcu+FFTvpn0Aw+c81RpL3LOl7/x+8yFjP5sBD3778KqpWt44OynUvag7WxcnY/hNDBtFp65vS42rsonGt76uW074WCE32cu5N0nPubzcdO36lk1MZqslMKX4yXkDxOLxlBK4clw88j0u2yv//a9Hxl95pNJW7lcHheXPz2MnXbvmHA8FAjx82dzCfnD7HlYLxo33zI688sXv9r+ubk8ThbNWkLTwXlb/wGFENsFCdxbweFw0Lx9szp5l9aa5659NWnoN+QPM+bG13jy+/uYOn5GlReQVRSUw6EIdxz/ADvt3pE/f/or5XWpGA4Dt9dFKBBOWZY0Eorw4k2vp0ybCtC5Vwe67b0Tn4+fTixivzK7qive7TRv15QnZo7is5e+ZPGcf+jaZyeOumAgOU2yk641TZOnLh+bFGyVUuw9qDdHnjcg4fi8GQu447gHQFt/lrFojGGjz+DEK48GoGXH5hgOI2nRoxkzadJagrYQYgsJ3A1EKBBmw0r7tO5/z1sGQOH6wpSB2Om2MqrZBYeUtPXFYOPqfDwZ7krni8tr1q4J4/56ijlfzue9Jz9h1pS5Se9OJ+CuXraObn27pAzaNcFwGAwbPZS8Fo0Yeqt9opSy8tcU2Nb+1lqXFncpKSjh1Tve4quJ31GwvjBpSuLlWybQ68AedOndieOvGMTU12Yk/B47nAYtO7Wg616dt+7DCSG2KzLH3UC4vS68NvWpgdKV7H2P2st2X7LD6WD4g2fywtyH6LJnpyqvUt64Kp9rxlxU5TYffMr+OBwOOu3WHuVQ6X9hKMdf4OeHj2anVTBkM0+GG4fTsLKclWE4DDwZHhzOxONKKd564IPU6U/LycjxkerSRs1yiEaiXNVvBB+9MJWCdclBG6zRhs9ftXKRd+jRjtvfupZGzXLwZnpweVzssl83Hvh8hKwqF0IkkMDdQBiGwYnXHIOnXPD2Zng4Y8TJAPQ5Ynd23b8b3swt13gzPQy+9AhOuPJo/v1jBf/+vjzt4FT6bofBgSfvxwlXHYXDGf8rk0YsmTbxO4ryi7lkrxv58cPZVXpnefmrNyUlc6kong0YegBvrx7LfZ/cxsGn7s8+x+zJefeexs2vX4nT5UgqAxqLxlj+1yr+/GlxWu3xZXo56JT9cJfL1e6NF3eZ+eFs1v67PmVCFrC23JXdRrbP0XsxceUYnp31AK/9/QyPzbibvJa1W5NbCNHwyFB5A3LWHScTCUV4/6lPAY3D6eDM20/m8HMOBqzgPuqjW/hq4ndMm/ANLo+Loy4YSN+j9gTgnUc/rLTEZnkOl4Nue3fBXxjg0sfOo/+J+/DVm9+S2yyXlYtXM3PyzymfWbyxmFdvn2ht9apBLreTVju1YOWS1SmnBr5681syczO46KGzS1O3hgIhTm5xAcFi+z3XSlnpbL0ZbiY/O4X1KzbS96g9Oezsg2xHO656bjiB4iA/fTInXtwlxsnXH8thZx/EyyPeJJDiPZt5s7z0PylxH73D4aBdtzbp/DYIIXZQtZqrvKZsa7nKa8OGVfnkr9lEu26t8fjsh8Q3CwfDrFi8mpkfzWbZ/H/pfchuDDyjf1IJy5JCPzMnzyLkD9HniN7cMPAuVv29pkrtsipuOeKpUxXa1ChlDb+fctNxdNq1PY9e+LxtkMrI8dGuWxsW/pxeLxYFbo+LSCiSchi69NnZXvwVJD0Ba3rhzf9eKF1c9v7Tn/LMlS+nvt7n5ty7TmXcnW8RCcVzhmd6aNauKU//eH/KnQMbV+ezbvlG2nZtRWZOBgCfvvQlz179SsovNd4sL3sd1os7/nedZEgTQiSpKFe5BO56VlJQwqjTH2fe9AU43U5MUzNs1Okcf8VRKe9ZPHcpV+x7a8IwbHbjTF5Z+CS5TXMAa3vRyBMeRCmFaWq0adK5VwcWz1maXPTDThrFPRxOB8MfOosXb36daCh5SNjtdXH4uQfz0fNTK38f0KlXBx7+ciQnNx9W5eF8OxnZPnoe0J2/Zi8lI8eHqU1WLU79xWXvI3uz4PuF+AsTK5G5fW7OumMIp910fNrvDhQHOLPTpRRtLEn4LMpQKKXIyMng1BsGM+SGwVKJSwiRZFsqMiLKue+MJ5g7bT7hYAR/YYBgcZCxt0zgp0/npLzn5iPuTZo7LcovYdRpjwMQ9Ie488SHCJaECBQHCflDhIMR/v7tX3zZvtJ52QoXPaURN2PRGK+MmIgZtf8i0O/4vpx87bFpLa5ShqJL745orVMmQkm4XqlKr/MXBZg1ZR75azax4q9VFQbtVju1YOhtJ9qeCwfCfP3OzJT3Lp6zlEcueI5bBo3i/ac/JVASxJfl4/Fv76Xb3jvhdDlwOA0Mh4E2NWbMpDi/mNfvncRTl42t9LMKIURZErjrUf6aTcz5cj6RckE45A/x1oPvp7ynYF2h7bm5X80nEo4w+/N5tsEyEoyw7zF7ccqNx9Fj/24cOGQ/cpvlbFXWsWBJMCkX+mbzpi+gVecWnHbz8ZUuZtOmZsb/ZjLq9MfTSqbiyXCXji7YUQ4rsKezkl0pePCLO8jOy065PS2rUWbSsfy1Bdx2zH1c2udGPnt5GrOmzGXszW9w6V43UlLop123Njz1w/38b81LHHbOwUl/JiF/iM/HzyB/bUHSs4UQIhUJ3PVo07pCnG77YdJUe7YDJanndbXW3DDgLkIB+3nVzUO259x5Kk98ey8jJl7Ds7MeYM+Bu2E4DWsYN43ebrpKCv3Mm76A8+49nROuPAqX20lGjg9vpoeWnZsnbe8KB8LM/Wo+LTu3wKikHU63k7ELHuPYSw7HcFi92c1fDgyHolGz3JRJX8o7cthAWnZoTvvubWjW1j7ZicPp4M6THuLV299k/YoN5K8tYHiv6/jpkzkJ8/Ehf4i1/67ng2c+Kz2W1SiTJXP/SVrJDtaXqTvKlVAVQoiKSOCuR212bmVboMJwKLRp8uodE1n73/qEc+v+3VDhM5fM+wd/YSApDSdYW8MOPrVfwrHm7Zpy/2cj+LjkDS557LwaLa2pNfz7xwqUUlz62Hm8tvRZbhp/BQ9/dRcnXnW07aiANjUbVmykUYtGGM7kv54Ot0GTNnkceuaBfP/+Twy7/wzeWfsS1798Kde9eAkvzHuY9zeN56Ah+yXt1bbjy/LS7/i+gDX83qSNfeCe9flcvnvvJ95+eDLn73I1L970GsX59tXEwsEI3076MeFYx13bpRzZ+POnv7hi31tZv6LiP1shhAAJ3PXK7XEx7IEzkoKlGdOsXLKGtx/8gGE9rmb+t3+gtWbS4x9xy6B7K3xmsMTKh33508Pw+Nw4nA6UsoL2foP7sPeRvW3vc7qctO/eOr56vGYYhqJDj7alPzdp1Zj9B+9Ntz47kZOXnXIBWjQS5aAh+9muto6FTYo2FDH52Sk8c9XLnNxiGPec8ijrlm9g32P3ovNuHfBlejnu8kE43eV2O6rEvd8ur4u23VrT54jdAWs0Y970BfYfJt7USChKoDjIjLdnVrjILzsvK+HnU244DneqL0Xayoz3bpn650IIkYrs465Hv3z5G6+OmIjhUFYqUjM+v7o5SISjRMJRHjjnaY447xAmjn6fiM3q7bKUocjJy2LQ+QPp2a87X77xDf6iAPsP3pvdD941qZdbUuinZFMJ/qIgd570cPLwchqry1Np3aUlux+8q+25/Y/rk3LLl9Pp5IePZqdMXhIOWl8uyhZbWfDdn/zvock88f0o2ndvQ9udWzHqo1t4eNizbFiVjzY1ewzoSedeHZj+9veg4dCzD+KUGwaXfr5ZU+al/XErKnPq9rk5/opBCcc67tqOUZ/cygNnP8Xaf9cn3RMNR6uVD14IseOR7WD1ZNO6As7qfFlaCVFcXheGYVRaMxrAk+Hh4Wkj6d535wqvC5QEefSC5/ju/Z9RhkLHzKRFcltjl327ct8nt9ou6trs5RFv8uZ979qe82V7E7KKpUMp6HXQrjw87c7SY1prNq7ehDfTQ+H6It568AMW/ryYtt3boNB8//7PhIMROu/egX7H92XCfZPSqoZmOA1cblfSn4lSijNGnMQ5d51qe9/KJau5cLdrS798bOZwOjj6okO54qkLqvSZhRDbp4q2g0mPu55Mf+t7zDQXT2nTRFey1NrlcaKU4oL7h1YatAHuP+MJZk2ZW2kPviJOtwNvhpfmHZrx97x/So+369aaUR/fkhC0SwpKyF9TQPMOzYiGo3z/wc/kNsmmZadmrFu+Mak6WFWDNlhz6vOmL+DMTpdy/n2nc8hp/a1561aNWTr/X67qdxvhQIRYNMbiOUsT7l0y9x+WL1qZVnfb6Xay3+A+7LR7ByaMehfDYRCNxGjbtRW3v30t7bu3TXlv652sUYh50xckBG+Xx8VJVx9T5c8shNjxSOCuJ0Ubi4kEK59PNhwG3fp2ZcncpbbnMxtlMPSWE2nbrTW7HbAL2Y2zbK8ra8OqfGZNmVftoK2UNRw8aNhAdtm3K49e+HzC+VVL1zLqtMcYPeV2wqEIj1/0AtPf+t7KER4z0aaJ0+UkFo1hmrrGi2isWbaOB895mmULlnPevacD8MJ14yr9MhANR9ll3678NftvIuEIZkzj8jjJbpJNcX4JDqe1D7tdtzZEIzEmjHoXh8tBLBKjXfc2PDxtJD9+9Av3DX2C4vwS9jlmL84ccRKNWzRKeM8d71zPc9e+yhfjZxANR+nUqwNXPTec1ju1rNHfByHE9kmGyuvJ/G//4JZBo5KGyg2HYQ3DupygIK9lYx7+6k6mvDyNN0e/nzA06/G5GfXJrex+kP08cioLZy3hhgF3VppLuyLWULGTWNS0nYt2eV2MX/w040e+xZcTviUcqFpJ0JpgOAyGXH8sA884kCv2vSWtsqS77Lsz1429hKnjZxAoCdLv+H3YY0BP/v1zBf/89i8tO7dg1pQ5vHnfewm1uJ1uJ83bNWXj6vzSP1Ony0FO0xxe/O0RcvKSa3rHYjFiUTP1ojUhxA5Lhsq3Qbv2685eh/Vi9tRfS/+h92Z62PeYvThr5Cl8/b+ZLJq9GE+Glz9mLuLUm4/Hl+3jzdHvUbC2kHbd23Dxo+ekHbQDxQFeuH48X7z2NeFQpOIh4TRWaJlRk1A0dSB0uZysWbaOL17/ZqtWqvuyPASKq1YYZTMzZvL2Q5N5/8lP00oy44wXVOnQox0XjD4z4VyHXdrSYRdrCPzOEx9KCNpg9dZXLlmdeCwSozi/mI9emMrQW5KzsjkcDkl3KoSoMgnc9UQpxe3/u44Zb33PlFenowzFkecdwoFD9uODZz5j4gPvEQlGME3NDx/OZtd+3Rj18S2ceNXRlT7bXxRg09oCmrZtUtqbu/Wo+1j485JKg6jhNPBleyjJD1R4XWVMrWncPLfSRCqVadm5Bbv07conY79IOmc4Dbrs0YlFPy9Jeb82NaFAGJXGxkeX181J11jzzJFwhJ8/ncvG1ZvYdf+udNqtQ+l1wQqS4JQXDkaY88VvtoFbCCGqQwJ3PXI4HAwYegADhh5QeqxgfSFjb3o9YeFSsCTIgu/+5Pv3f+aAcmUgy4pGojx1+Ut88doMDKcDBZx916nsflAPFs9ZmhS0lVKJBTCUVU96a4O2y+PiwtFn0KJjMzwZ7qTeaVWs/nstlz1+nm3gNqMm+as34XQ7Kl0Jrs3kz7uZw+lgtwN34dLHz6Nlx+YsX7SSaw68nVBJmFgshjIU+xy9F7dOuAqHw0HfQXtYiwvLpVO1e77hMGi1U4tqfHIhhLAnCVi2MfOmL8DpSv4+FSwJ8fU7PyQd11rz1cTvuPqAEZza+kKmvPIV4WCEYHGQQHGQV2+fyGcvf2WfpUxrdt6zE627tCSnSRYHDtmffif03erPEAlFKN7kx+FwcNEj5+CxqWWdLm+Wlx8//iXl+XX/bcDhdJDZKKPSZ6UK2ufdezpXPTeccDBCJBzhugF3smltIYGSIOFghJA/zA8fzebTsdMAuGD0meTkZeH2WWVUHS4H3iwvbbu1Skrj6vI4OaHcnm4hhNgaErjrgWmaTH5uCud0vYKTmw/j/jOfYM2ydQB4M722BTmUocjI8ZG/ZlNCLvJnrnqZRy98jgXfLaRwQ3FSPuyQP8TPU+bYbj1z+9wcOGR/xi16iknrXuHWCVexYtGqGvmMr94xkb9/XcbhZx/MCVceVe0c6OFAiP898mGl151z16m8uuhJrnr2wpTX2H15cbodfDJ2KhfvcT03HnoXxzU6m402eeLDgTAfj7HKkzZr24T7p4ygaevGoKy59G59duL2t69ltwN74PK48GZ4aNwil9vfujZhmF0IIbaWDJXXg2euepnPX5lOML5CfPrE7/j5s7mMnf8oewzsicNmIZXT6eC7935k6vgZAAwY2p9jLzmcj1+YWml97eL8EnbeqxMLf9oyx62UwuNzc9QFA5k3YwHPXPkyS3/7F4erZhZLaVPz3hMfc+mT5/Pek5/YFvxwe13EYqb1ZSPFYriSgsqH7UP+MKuWrKFNl1YYhpFy6Nztc4PWpUP3hsMg5A+zsoJyn2VtXJ1PoCSIUopbjhxF4fpC0FZP/rdv/mDEsaMZt+gpijeVUFLgp1XnFrZpW4UQYmvIvyp1LH/NJj4dO600aAOYpiZYEuT9pz7F5XYx6pPbyG6cSUaOD1+2D6fbialNCtYXEQlFiIQiTH1tBlfsc0ulQVsp6N53Z+77+FYOP/dgvJkeHE6DPQ/djSdn3sfaf9dz29H3s/S3fwGSEqFsjdX/rGP+t3/isCkWAtCpVwcmF77GURccutXv2ryiu1XnFrTp0ipp1MKT4WHYfadz9p2n0KJDM1S8p1wVheuLGNruYv738GSC/lDCKIYZMynaWMyPH/9Co2a5pV8ihBCipkmPu479/esy3F5X0kKxSCjKb9/8AcAu++zMW6te5JcvfiNYHOTTl6cx+/N5CdebKepGl6WUwpPhZtj9Q/Fl+bj6ueFc/dzwhGvuOvnhWtljbTgN9jp8d9xeV8redPN2TXF7XLTo0HSr3zd32nyC/hDeDA8jJ13PtQeNJBwIE41Yi8v2PrI3gy87EofDwWevTEuZJ70i0UiM4k0lTHr8I4I2e+AjwQgr/qqZqQYhhEhFAncda9GxuW3CEsNh0HbnVqU/u9wu9jlqTwDG3/V2ld+T2yyHnv27c85dp9KpZ/uU1y2d/2/KKl220qzCoVAcPfxQMrJ9VvWzosQhb2+mh6MuPBTTNPn4xeQV41WllKJwfSHe9s1o160NE/59jh8//oWNqzbRY/+udOndCbC2eS3fynn8cCBsrZYvl9DF5XHRabfUv9dCCFETZCyvhsViMQIlwZTBsO3Orei+z864ypWcdHlcnHTtsbb39Niva1oJRDZzOK00nBc+cGaFQRugS++OVVo4ttPuHfD43JWuFHe6naz4axUlBX669ulcetxwGDg9TgZfdiR9Dt+d377+g8L1RWm/PxXDadC4ZaPSn11uF/1P2IfBlx5RGrTB+r3xxFeDp+Lxudnz0N1sFwmC9SUhOy87YQW50+2kZafm7HX47lv1OWZPncfdQx7mlqNGMXX8DNu66kKIHZsE7hoSi8YYc+NrHN/oHE5ofA5ndr6UmR/ap2m96/0b2XdwH1xuJy6Pi+YdmnLXezfQcdd2ttcPvfUkPBmetHN6x6IxSgr9jL3pjUqvPeO2k3B7Kw5km7Xt1ornf3mYZ2c/yMh3byC3eU7Ka6ORKNcPuJOTWwxj1pS5pce11uQ2yS6tnrVh5UZiaQz7V+askafgcleeOtQwDI65+HBrCN+Gw+Xg0ifPY+e9dkoVt9Eanpv9AAPPPJCMbB+ZuRkMGjaAx76+e6vmtV++bQJ3nvgQ30z6kVmfzeWh857mzE6Xsva/5DKgQogdl+QqryFPXPoiU8dPTxg+9WS4Gf3ZCHr238X2nkBJkGBxkEbNcysNyv/+uYKXb5vAb1//TnaTHEL+IEUbSyos9ZmRk8EHm8ZV2vbff1jEc9e8yqJZS6zetyZpW5k308ML8x6mdWerEMbiOUu55qA7bOd6K+PL8nL1Cxcx4PT+rPp7Ded1v3Krgvdeh/di9Ge3A9Y++P89Mpn1KzbS54jdOemaY2ncPJcNq/J568H3+WXqrzRp3RhlKH6d8TuxiFXoxOly4MnwsPvBuzJ76jzCwYjtSniA7Lws3l3/SrXba2ftf+s5r9uVSeU+ATJyfLy68MmkYiVCiO2X5CqvZSWFfqa8+lVSta+QP8zr97zD6Cm3297ny/Tiy/Sm9Y723dtw56Qbtjw7EOLzcTP4ZtJM5n61wDbIZOR40bry6ls99u3KUzPvY8b/ZvLw+c/YFj55+sfRpUH7k7Ff8PTlL1W7fnegOMiyBf+VbqMyY1v35fHaFy8pbdezV79S+uXp39+XM+WVrzj91hN58cbXrWFnDct+X44nw8P5o4ay+8G74nAamKbGNE2uPeCOSjO91cY89txp8zEc9n9OgeIgkx7/iAvuP9P2vBBixyJD5TVg46p8nE77/c//LVpZK+/0+Dwce/HhPDh1JAPPOMBaAFaWgvXLN3Jy8/N557EP01qA9u27PyQFbbDSoP4z39ouVryphGeufKXaQXuznz6bw9X9R/D4RS9U3LZKZgc69myPL8vL4rlLef7acQkjHpFwlE1rC3nu6letBYFlXhPyh3jtrv/RYdd2dO7VkS69O/Hjx79U+rk8Ge7SYf6alNU4k1Tl2bWp+WXqrzX+TiFEw7TD9bjz1xZQsqmEVp1b4EgRbKuqefumtpnJlFJ03XOnSu9fPHcpMz+chcfr5uBT96d5+2ZVev+Vz1zAxtWbmP/tnyilrOHzeHMKNxQz7va3iIajnHbTCRU+J6tRJoahkj+LsoZrAX6d8TtOt4Nw9SuCArD4F/v64kkq+b6x7Pf/OLXVhRhOI62yneX99+cKdtq9IwBOp9P6/Cm2srft2opLHjuPXgf2qPJ7KtPniN64XE7C2HwGZf0dE0IIqKcet1LqGqXUAqXUfKXUm0qp9MaLt0LhxiJuGXQvZ3S4hEv63MQprS5gxtvf18izPT4Pp954XNJKa7fPzVkjh1R473PXvMLV/Ufw+t3v8MrtEzmv+1V8Pn56ld7vy/LxwJTbGTPvYfLKrKzeLOgPMXH0+8RiFSdXGXTBQFw2i7YcTgd7DNwNsHqc2xJtaiLhaLWCdiQSpVHz3NKfDxyyL4bNlzm318WEf5/jlT+fpO+gPbaqvam4PS4e+nJkUq5zsFa5n3zd4Fp5rxCi4anzwK2UagNcCfTRWvcEHMBptf3ekcc/yNyvFhAJWQU4CjcU89D5z7Dw58U18vwzbz+ZSx47h5admuPN9NDroB48Mv0uOvdKnad6/rd/8PGLXxLyhzFjJtFwlHAwwhMXj6FwQ9W3SLXp0oqC9YW250KBMP7CitOHdt1rJy4YfSZur4uMbB8Z2T5ymmRz/2cjSguflBT4CaSxIC3dFfAVPsNgq8uCVmT3A3vQpFXj0p/bdGnF8Aetz795y5vb6+bK5y6kWdva7/HuvGdn3lz+Al377ITT5cCX5SWrUSZXPz+cnv261/r7hRANQ30NlTsBn1IqAmQAtTMRHLdi8Sr+mv13UuKTcCDC/x79kBFvXrPV71BKcfSFh3H0hYelfc/0t76zzVrmcDr46dM5HHrmgVVuR9uurVk0K7k+tS/LS2Zu5RW0jr98EAOG9ufXGb/jzfTS+5BdS4P2tDe/4dELn0+52nozl8fJWSNP4c373yVQVP0x9c69OrB84apqlQX1ZXsrfHdmbgY3vXZF0vHjLhtEv+P7MnPyLJRhsP9xfchr2djmCbWjUbNcnvlpNBtX51O0sZg2O7eyrRYnhNhx1XmPW2u9AngY+BdYBRRorT8vf51SarhSapZSata6deu26p0bVubjdCf/46e1Zs3StVv17K2hDKPCJB/VMez+M5ISjHgyPJxz1ykJe4w3rMrng2c+493HP2bF4lWEg2GK8ovRWpOTl03/E/ahz+G7JwSNF296Pa0hadPUDDi931bnPf/39xVc/Ng5Vu8/x2f1viv5bVFKsfuAXa296RVcGwlFeGz4C7bnmrZpwrGXHMExFx1Wp0G7rLyWjenQo50EbSFEkvoYKm8MHAd0AloDmUqppH0uWusxWus+Wus+zZpVbbFWeZ17dUjKDQ5Wz3Dz3G19GDD0ANvkJ7GYSd+jqjeXuufA3Rj57g102q09TreTFh2bccUzwzjusi01oae9+Q1n73Q5L974Gi/e/Drndb+KY7PP4pSWF3ButyuZN31B0nOjkSgbVmxMqw2GoXD7PJxy4/F4t6IWdzQaIxYxOf++oVz86Lk89eP9HHvx4VYyGkPRpmurhEQqhsPAm+lh78N7WyMZFQwMhIMRZn0+j2V/LK92+8rKX1vAWw+9zxOXjOHLN74hbPP3TQghakKdJ2BRSg0BjtRaD4v/fDawr9b60lT31EQClldHTmTSox+VbndyOB1kN85kzG+P0rjMAqW69vJtE5j0+MfomInhMNDAja9ezkFD9quV9xWsL2Ro+4ttE31s5s308MzPD9C+e5vSY1prTm5+PoUbiit9R4cebRk7/zG01tw95BG+fe9HFGA4HDgcBrFYLO2EK95MT3y7mKLHfl2598ObcXvdpfvT5341nzdGTWL132vpsX9Xzrz9ZKZN+JbX73mn0me7vC6uH3sJA4YekFZbUvnzp7+48dC7iUVjhIMRvFkemrVtylMzR5GZm7lVzxZC7Ji2tQQs/wL7KqUygAAwEKj1tGjn3HkqnXZtz/8emUzBuiL2HtSbobedVK9BG+D8UUM57OyD+OGjX3B7XRxw0j41PjwbDkWY/MxnTB0/wxoOr2SOOhyM8O7jH3H18xeVHlNKccaIk3jltokJJUnLU0px82tXAjBn2nwr3am2Or+xaMzKyFaFWYCy+8oXfPcnE+57l3PvPq10KqH3IT3pfUjPhHs69GiLL8tb6SK6SChCy07N02+MDa0195/5ZMK7gsUhVi9dw5v3v8cFoyVpihCiZtVLylOl1F3AqUAUmANcoLVOGQ0aQsrTbZVpmlw/4C4W/by4Sou82nVrzVEXHkrP/t3ptncXlFJorZn0+Ee8fs87lGzy296X16oxb60YA8DdQx7hm0k/JF+UZoUxO01a5zFxuf3c9GaRcITzul3F+hUbk1K3lvd57O0qryfQWhMKhHF7XWxYmc+5Xa+wHcFo1q4JE5Y9X6VnCyEEbHs9brTWI4GR9fHu7dnU12bwxj3vsGH1JnbavSPDHzqLYHGQxb/8XeWV2f8tXMmYG1/D7XHRZY9OdOjZFjQcclp/Jq17meG7Xce/f65IuMflcTJgaP/Sn/2F9sHdcBgp64krQ2EYKuVQeiRc+dyxy+3iqR/u49mrX+Gbd39MuUjOl+2tctD++p2ZPH/dODaszMfldrLv4D4p98evX76RhbOW0K1P5Ul4hBAiXVJkZDvxzmMf8ertExOKjngy3Aw4/QA+felL23sMw8A00y/uoZS1Qn3gmQdyzEWHcf0hdxINRwkFwviyvDRt24QbXr2MT178gn9/X44vy8v8b/9M+tLg8rpwOh0ES4KU/etnOAzMmInD5cCMmkmpUB0uB4efczDXjrk47TabpslJzc6nOL8k6dz5o4Zy+i0VZ5Mra9bn87jzxAeTVtYrBan+b5TXqhETl4+pkX3tQogdxzbX4xY1KxqJ8tpdbydVCgv5wyz4fiGeDE/SOXeGm32P2ov2u7YhWBTk63d+YO2/FZeP1Nqac/7itRkMGjaQ8Uue5ss3vmHl4tXsun83cppmc8OAu4iEIpgxE6fbgRnTuH0uwoEIhsPA5XZyxbMXsMs+XXntrrdL27d22brS4eayPWSX10UkGMGb6SGnaTbnjzq9Sr83hmHw1Mz7uKrfiC1JbRQcfs7BVQraAOPvfMt2O1xF3303rtrEv38sp0MP+5KtQghRVRK4twMbV29KObRcsL4QhzN515/b7eL6Vy4trU528Gn9uf6QkbZFRsoLByP88OEszrnrVE644qjS4+ftclXCF4RoOIZSig67tKXNzq1p1CyHoy4cSKfdrGxyt8UT39xx/AMsX5icg8eX7eWgIfujDEWPfbty8Gn9qrW9rG3X1ryz9iWWzP2H/LWb6Nl/l7SrspW1csmaKt8DsHjOPxK4hRA1RgL3diC3aXbKbl/rnVpy2RPncc8pj7JpnZUONa9lI+5457qE4JXXslGlC7k2czgdeMsFvkBxgFU2gU1rzfK/VvPsrAdTPs9lkxwHrBXqh5zenz2rsdd+4+p8pr/1PcWbStjrsN3psV9XuuzRqcrPKavTbu2ZO22+/ckUC+6UUuQ2y96q9wohRFkSuLcDHp+HYy89gg+f+zxpjvvsO0+h295deO3vZ1i+aCVKKdrs3CppzrVZ2ybsfvCu/PLFb5ixiue9DUNx8Kn7JxxzeVwYDsM2+GfGK4ulcuT5A/jp0zm2dcB7HbhLhfeWp7Vm0mMfMfaWNwBr2P1/D09mn6P34tYJVyVkj6uq8+49nesOHpmUOhfA6XQQjcaSgrfb66qVamJCiB2X1OPeTlww+gxOvPoofFleHE4HTVo35vqXL6XP4bsDVs+vXbc2tO3aOuVCqRFvXUv3fXZO+Q6H08DtdXPNixfTokNiNjuny8khp/dLyGQG1mK2E648ior0OaI3R10wELfXhdvrwpflxZft5e73b6pSys9oJMqtR43ihevHE4vESufKgyUhfvx4Nt+991Paz7LTY9+u3PfJrXgzE4frPRkeLnjgTPY5as/STHgOp4HL6+LWCVfbZscTQojqklXl25lYLEawJERGtq9aK5ljsRinthlOwdrEKmMuj5OTrz2WIdcPJrtxlu29QX+Ie055hLnT5uPyuIiEIgw880Cueu5CHI7Ka58vX7SSX774jYwcH/sftzcZ2ck99VVL17Dir9W036UNzdslVux6/+lPGXPDa7bpbQH2G9yHu9+/yfbchlX5zJ02n4wcH3sdvjtuT3J5083CoQhfvfkt30z6gey8LI69+HB67NcNrTVzps3n58/mkJOXxcAzDqhybXUhhICKV5VL4BZJFs9dyk2H30M0HEWbmmgkSp8jenPk+QPoc0TvCoMaWMF19dK1dOjRtsaywIUCIe455VHmfPkbLo+LcChC/xP24cZXLyvtlV+85w0smftPymfsf9ze3PXejUnHJ9w3idfvmYTT5UApheEwuP+z2+jeN/XogxBC1KaKArcMlYskXXp3YuLyF7h1wtXsMbAXYKUvfeDspzit9YUsmp1cNrSsVp1asMeA3Wo0devz141nzpe/EQ5GKCnwEwlG+P79n5gwalLpNRXNzbt9Lo4495Ck4/O//YMJ971HJBQhUBzEXxSgeFMJtx19H9FI8lx2fYnFYmklnxFCbP8kcAtbLrcLl8fFL1/MIxKKEiwO4i8MUJRfwm1H3Zf2CvSaEIlEmfLKtKS0oqFAmMnPbakIe+hZByWVNAVAwSGn9mO/wclfXj8Z+yXhQPIWuGgkZlslra75iwI8cO7THJt1JsdknMFlfW9m8Zyl9d0sIUQ9ksAtUvpkzFTbfd3hYIT53/1Z6+8PhyI8dflYjm90DpGQfe83UBQo/fXxlx9Jlz074c2ytqo53Q5cHic3vno51798me2cf6A4kDKBSjq1x2vapnUFLPx5MYUbrWQxI465nxlvfUckFMU0NYtmLeHag+5g3fINdd42IcS2QbaDiZQCqZKxKKx61zVMa83cr+azZO4/NG2bx8u3vWm7N7y0GQp2K7PVyu118+iMu5n9+Tx+/eYPmrbO45DT+5GTl3of9UFD9mfWlHlJX1Ci4Ri9Dqq7bVzRSJTHLx7DtDe/xeV2Eg1H2W9wHxbOWpL0pSUajvLhc1M4f9TQOmufEGLbIYF7B7P6n7Us/HkJTdvk0WO/rhWuPD/41P35dcaCpKAWi8bo2b97jbYrUBzg+oF38d8fK4iEo8RiMXQs9cJJp9uB2+vmkkfPSThuGAZ7H7kHex+5R1rvPeCkffn05Wn8/v1CgiWh0rSslzx+LlmN6q6W9qu3T2T6xO+IBCNE4lMC37//s20RtUg4yt+/Lquztgkhti0SuHcQpmny+EVj+PKNr3G6nGitadImj4e+HEnT1nm29xxyWj8+f3U6f/60mGBJEIfTwOlyctVzw/FlVZxUparG3fk2S3/9N+VWrrI8GR6OufgwTrr6GJq1bbJV73U4Hdz3ya388OFsvn3vR7IbZ3Lk+QPp3KvDVj23KrTWTH5uSlIxlohNohewkrpUtN9eCLF9k+1gO4hPxn7Js1e/kpBZzXAY7LLvzjz+zb0p74vFYsycPIvvJ/9MTpNsBg0bSIdd2tZ4+4a0GFaakrUyR5x7MNe/fFmNt6G+mKbJka7TkqqhgZU4x+V1lU5NKEOR1SiTl/94nEbNcuu6qUKIOiLVwQQfPP1pUoUwM2ayaNbfbFydn3LrlsPhoP8J+9D/hH1qtX2xStKslrbH5eC0W06s1bbUNcMw6Lhbe5baDH9369uFPQb25JMxXxD0h9jz0F5c/Mg5ErSF2IFVGLiVUjlAM631knLHe2mtf63Vloka5S+z+rosw2EQKA7WcWuSHXjyfkx59SvbPOCbGQ6Dez64ibY7t6ry8wPFAT4fN51fvviNVp2bc+wlR9CmS9WfU1uueGoYtwwaRSQYxjQ1hsPA7XVx+ZPn023vLpx/ryxEE0JYUm4HU0qdAvwJTFJKLVBK7V3m9Ku13TBRs/qd0Ne2CldWowxadW5RDy1KdP59p9O8fVN88a1c3kwPngw3jZrn4s30sPeRezDm10fSXnRWVuHGIob3up4Xb3qD7z/4mfef/oyLet/ArM/n1fTHqLbdDtiFp2aO4qBT+9Fx13YMGNqfp38aTbe9u9R304QQ25iUc9xKqbnAIK31KqVUX2A8cIvW+j2l1BytddX/Ba0mmePeeoUbi7h0r5vYtK6QkD+E0+XA4XJy13s3sNdhu9fKO7XW/PzZXD4Z+wUhf4gBpx/AIaf3S1k4JBKO8N17P/HXL3/TpksrDj6tn22+8qoac+N43n/y06TFXnmtGvHmfy9sVcUwIYSoDdXKVa6U+k1rvVuZn1sBHwHjgHO11nvWRmPtSOCuGYnDxS0YfOkRtN6pZa29b8yNr/Hhc1NKt5N5Mz1077szoz8fkVbRkZpydpfLWfV38n5wT4aH5+c8VK2hdyGEqE3VXZxWpJTaafP8drznfTDwPrBrTTdS1D5flo/jLhvEcZcNqvV3rVq6hvef/rR0TzJY5TX//HkxP378C/sP3ruCu2tWql67GTNLh+aFEKKhqGiM8BIgITuH1roIOBI4vzYbJRq+eV8twOFI/usVLA7yw0ez67Qtx11+JJ6MxBrahsNg570606RVzRVCEUKIupAycGut52mtF9scj2it36jdZomGLrNRpu3csdPlILdZTp225YjzDmHg0P64vC4ycnz4sry07tKS29+6pk7bIYQQNUH2cYtasc9Re2DY9LgdTodtec3aZBgG14y5mKG3ncSfPy1OK92rEEJsqyRwi1rh9roZ/fntjDjmfivrl7LmlK9/6dJ6WwzWokMzWnRoVi/vFkKImpJW4FZK+YD2WuuFtdwesR3p1mcnJq54gT9mLiIUCLNrv+54y801CyGEqJpKA7dS6ljgYcANdFJK9Qbu1loPruW2ie2Aw+GgZ/9d6rsZQgix3Ugn88SdQF9gE4DWei7QqdZaJIQQQoiU0gncEa11Qblj235JMSGEEGI7lM4c9wKl1FDAoZTaGbgS+L52myWEEEIIO+n0uK/AypQWAiYABcDVtdgmIYQQQqRQWVlPB/Cx1voQ4La6aZIQQgghUqmwx621jgGmUiq3jtojhBBCiAqkM8ddDPymlJoKlGw+qLW+stZaJYQQQghb6QTud+P/1RilVCNgLNATa4X6+VrrmTX5DiGEEGJ7VGng1lqPq4X3PgF8prU+WSnlBjJq4R1CCCHEdiedzGlLsdm3rbXuXJ0XxufLDwTOjT8nDISr8ywhhBBiR5POUHmfMr/2AkOAvK14ZydgHfCKUmp3YDZwlda6pOxFSqnhwHCA9u3bb8XrhBBCiO1Hpfu4tdYbyvy3Qmv9OHD0VrzTCewJPKe13gNrwdvNNu8do7Xuo7Xu06yZVHQSQgghIL2h8j3L/Ghg9cC3phzocmC51vrH+M/vYBO4hRBCCJEsnQD8SJlfR4GlwCnVfaHWerVS6j+lVLd4mdCBwO/VfZ4QQgixI0kncA/TWv9d9oBSamurg10BvBFfUf43cN5WPk8IIYTYIaQTuN/BmpMuf2yv6r40Xhq0T2XXCSGEECJRysCtlOqOVVwkVyl1YplTOViry4UQQghRxyrqcXcDjgEaAceWOV4EXFiLbRJCCCFECikDt9b6A+ADpdR+ko5UCCGE2DakM8c9Ryl1GdaweekQudb6/FprlRBCCCFsVZqABXgNaAkcAcwA2mINlwshhBCijqUTuLtorW8HSuIFR44G9qndZgkhhBDCTjqBOxL/301KqZ5ALtC89pokhBBCiFTSmeMeo5RqDNwOTAaygDtqtVVCCCGEsJVOPe6x8V/OAKpVylMIIYQQNaPSoXKlVAul1EtKqU/jP/dQSg2r/aYJIYQQorx05rhfBaYAreM/LwKurqX2CCGEEKIC6QTuplrrtwETQGsdBWK12iohhBBC2EoncJcopZoAGkAptS9QUKutEkIIIYStdFaVX4u1mnwnpdR3QDPg5FptlRBCCCFsVVQdrL3W+l+t9S9KqYOwio4oYKHWOpLqPiGEEELUnoqGyt8v8+u3tNYLtNbzJWgLIYQQ9aeiwK3K/Fr2bwshhBDbgIoCt07xayGEEELUk4oWp+2ulCrE6nn74r8m/rPWWufUeuuEEEIIkSBl4NZaO+qyIUIIIYSoXDrbwYQQQtQyHf0Pwj+BkQueA1HKXd9NEtsoCdxCCFGPtNbootHgnwDKgTUb6YS8cShXj/puntgGpZM5TQghRG0JTQf/RCAE2g+6BHQBOn84Wpv13TqxDZLALYQQ9Uj7JwIBmxMlEPmtztsjtn0yVC6EEPVJ2wRtABToYJ02xY42S9CByRD9DRxdUBknooxGNff82FqI/QOO9ihHyxp77vZMArcQQtQj5TsWHZ1nE8A1uHvXR5O2tCC2Br3hJDCLsEYFvOiSZ6HJRJSzy9Y9W0fQBbdC8FNQHtAhtPdQVO6DsjCvEjJULoQQ9cl3HDh3BZURP+AEvJBzP0p56rNl6KIHwNzAlqH8IOgidMGIrX928dMQnAKEQRdZ/xuchi56bKufvb2THrcQYrugdRR0GGVkVH7xNkQpN+SNh9A0dGgGGHko38koZ4f6bhqEvgJi5Q5qiMxF69DWfbHwvwGUnwoIQmAi5NxU/efuACRwCyEaNK1D6ML7ITAJiKId7VE5d6E8+9Z309KmlBO8h6O8h6e8RusYRBcCBji7oZRKeW3NcaU4brDVA7ban/K41rqOPl/DJEPlQogGTW+6Ph60Q0AMYkvR+RehI3/Wd9NqjA7/jF7XH71xKHrjaeh1h6AjC2r/xb7jgfK9aid4BqBUqqCeJtfu9sedPSVoV0ICtxCiwdKxNfHh3FC5MyF0yYv10aQap2Mb0PkXWnPN2m/9Z65EbzwHbabotdYQlX0NuHoCPus/lQGOjqjcuytub/ArdOQ3tE5dn0rl3B6f19+cXdsBKgOVc0dNfoTtkgyVCyEartiK+IrkcLkTJkQX10uTalzwQ9Dl55kBohD6AnyDq/VYHV2OLn4Swt+D0RiVeQF4Byf0dpXyQd4EiPwK0T/B0RHcfW17xFYGuEfA/yooN2CC0QryXrHd5qVcPaDJZHTJSxCZD65dUJkXoJwdq/V5diQSuIUQDZezI+jyvW0AB7h61XVraoWOrSN5RAHQkfiKb5tTOgqhGRD7F5y7gHufhGCrY6vRG44HXQyYYK5FF9wB0aWo7KsTnqWUAvfu1n8VCU2BwGtYq8TjX6Ri/6DzL0U1fdf2FuVsj8q9q+LniiQyVC6EaLCUkQe+IYC33AkvKvPCemlTTVOefcpsFSvLCe69k47q2Br0+sPRBdejix5Bb7oYveEktFmy5ZqSsfHFYWVTqgag5CW0WZj0zHToknE2e9FjEP0LHV1erWcKe/UWuJVSDqXUHKXUR/XVBiFEw6dyboesq8FoAcoH7v6ovLdQzvb13bSa4e4Pzp4kfjnxWRXEXD2TLtcFt0JslZUylbAVoKOL0MWPb7ko/DMQTX6XckN0SfXaaRakOGHE92mLmlKfQ+VXAX8AOfXYBiFEA6eUgco6H7LOr++m1AqlDMh7Ge3/HwTeBeVE+U61EreUo3UIwjNJ3nsdhsAHkHOb9aOjvTVnTbnFYzoMjhbVa6jnUPAvSX4mAbTRFFknXnPqpcetlGoLHA2MrY/3CyFEQ6KUGyPzDFTei+DsiC66B71mb8z8KzBD36FLF+dpkgPnZluGxVXmBSRv83KDe2+Uo3Wl7dFmCTq6zPqiUHp73xRXOyHwfqXPFOmrr6Hyx4EbSZxgSaCUGq6UmqWUmrVu3bo6a5gQQmyLtA6jNwyBwIfxYfBia0FY/vnotfuhg1NRyguu3iT/0+4E7xGlPyn37pD7EBhNsIbg3eA5BNXoyUraEMUsvAe9dl/0huPQa/tiFj9rJUzRa0laawBANN67FzWlzofKlVLHAGu11rOVUgenuk5rPQYYA9CnT5/UmwGFEGIHoAOTIbaa5LlpbeUP33QdNJ2Myr0fveHU+Gp7v7WwzWiCyr4u4S7lPRxtNIHoPHDuhuFJ1WMu86aix8H/P6za4fGDxS9YQ+GuXUjZ26/FnOs6tt5aQa+c1pcPY/uffa2POe5+wGCl1FFYX89ylFKva63PrIe2CCFErdFmCQQ/RceWoVy7gmcg4LT2RUd+BiMPPIejjKyKnxP5HQpHApEKroqiA29jZN8Izb6E4Mfo6D/We72HJ1Tc0mYhesOJEFuONfBpYDp3QTWZYO3dtmuDNiHwOsn5xQNQ8jyq2TS0ox3E/kq+OfgxOntEjeeRN0smQNH9oAzQCrgdnfsIhu+wGn3PtqbOA7fW+hbgFoB4j/t6CdpCiPqgo0sgthZc3VFG4xp+9j8JPV+tMkA1t/aeh3/ECsJuUPdC3quoFPvOtdboTVdRcdAGiFqfBawvAhmnplwQpjdeaO3xLmVCdAG64FZUo1TVuUIp9swDsfh+cmcX+8CNw1o05x1YyWdIn47+DUWjSej9AxRch/bMqPE/z22J7OMWQuxwtJmPueEU9PoT0JsuQ689ELPo0QpTdFb5HQU3gd4ExNOSaj+Y/0H4G6wymVHrnC5G519m9WjtxJZDbE0ab8xAeQ5CaxMd/Bwz/1LM/MvQwS8SPpdpRiE6x/4Rwc8qeL4XjOQMaAC4elj/qyrqC9bsunId+BD7LW0GBL+s0Xdta+o1cGutp2utj6nPNgghdjx60zVWmk2C8exhIfCPg+DHNfN8swQiv5E85xsjeasW1j7nVAu4lGHznPLc4OyE9hyOLrgBvelGKx1qaKqViKXw5jLv2ljBc+xSq8aboRQqZwSJC9AU4EVlW89XvhOsvfRJTPDsV8lnqCIdLyqTdNzENtPcdkR63EKIHYqObYDwLJJ6azqA9r9SMy9R1fmn1T44K0cbcLQldY/VAa49wdESih6B4FRKe/lg9fQDn5ZWE1NGo9TPUpmJLYrMx9x0HeaGoZjFz4J7L1TeS+De1+p9uw9GNZlorVIHcPcD73FYwd0Z/18vKvexlHPn1WWVQLVbxa7Bc3CNvmtbI7nKhRA7Fl0IymEfJ81NqW+LrUH7X7d60s7uqMxzUI5Wttcq5UO7+0D4JxJ3vW6uY13uS4PKsHKKp6AaPYneOBR0kC29SQNUC9D5EJmNNQfuTH42ABF0aAbKtStKudHeEyBokz88a8vKczPwCRTcDIStzxD5De1/E9X0A4y88Sk+t0Ll3o3OOB1CX1tfBLyDUI4mKT9bdSl3b7Tv+Pge8SDWlxE3ZF1mfdnZjqmanNOpLX369NGzZs2q72YIIbYDWsfQa/e3Al4CJ/hOx8i9fcu10aXownsh/ANWYDSwArELlBuV94ZV5cruPbGV8cVpxdawrnKD0RGMHIj+Gs8V7gVloBqPRbn7VNLuAASnomMrrdXozt5QcAnE/kvzk7sg43yrVCcmuvAuCLyD9Q3GCVlXYWRdGH9XJP57VD6NqQsyzsbIuSnNd9YurTVEZqEDn4FyoXyDU/55NDRKqdlaa9u/FNLjFkLsUJRyoHPuhoIbKO1N4gYjB5V1cel1OrYOveHk+Bz45g7O5t5zBHQEXXgXqslb9u9xtIZm09DFT0DgI6u37NoJMq9EmSvR4R9RRhPwHZPWCmilfOAbXDrIrWNrrXrkaYuAf7z1SVTM2pLm7g++k1DegSjl2nJpdCn2q9gjEJoGbBuB26pctjfKptjK9kwCtxBih2P4jkA726JLXrFWbbv7oTLPsKqNxWn/G/EFUBWMSkbmorVp5RO3of1vQslrWKvIsRa/hb6CJu9jlCufWVUahf2weEUC4B+D9U9/BPgDIj+h9R2ojJOs52qNjsyPD8vbMBpVu801xeppb94L3yS+Fz6z8hu3ExK4hRA7JOXaFdXo4dQXRH7F6pFXxEOqhV5ah6H4MUqDNgAmaD+65FlU7uiqNbi8otFUvtrctmUk9KZ1AIruQ/uOtea/ix4B/3hsM1IrHyrzvIqfHl1srQNwtAVXn4Q64DVB6yh60xUQ+p7EvfDjbKullW2XLnnJqn7m2hOVeS7KkWJ72zZOArcQQthx7VImUYodD/hOTB2YEhKclGXGF61VnzbzITiF6gVuuwdG0YFJaOdu4H+VlF9YnL3Bc4TtKSugXh1PP+qwDhotIO91lKNZuWtNMFeDykQZuVVrqv/deNDe/IUoamV9zb8cmn1l++ehQz+g84ezZaHdAnTgf9DkXZSzQ5Xevy2Q7WBCCGFDZZxpLShL4gC84O6LyrnZ5nyc0QR0iqBvVF6Bq0KxNVB2TnqrBaDwQdh4ChXUfoLonym/qOiSV6yV5ISshXfaD7F/0QXXJ14XmoFedyB63ZHotftjbrwQXcFq/iTBt0kcxdj84HyIJmdt01qjC2/HWnledo1CCbroofTfuw2RwC2EEDaUoxUq7w1w9aJ0q5FnEOQ+Gt8S9ZJVjSvV/UZj8AwguXymN2ERXFXpyML4HHRlKVBtODqmSJACUII1Z17BvLnOT51dLjCR5DzmMQjPQpuF1u2RRej8K8BcG782AuHv0fkXWT320Ex08Cu0WVxBG1KNMihsv3ToIoitsLnejNcub3hkqFwIIVJQrh6oJu+gdRRwVGm+VusIuPpYW8l0GGsLmQeyb0J5+le5LVqH0PmXWMljlANrmFyRPFyeB+TbHPehsi61qmkVP4kV5Cqbwy/H0Tn170GqxWyo0hzn2nYYPgKReeg1m7dxuQADnXMXRsaJyY/znQBFf5H0JUFlgbOrzeu9pOyjqoZZSUx63EKIBkfHVlvzlrFVdfI+pZxVDNoanX8hFD8Sz1ceD7KufTAyTqlWG3TxkxD+GStNawlb9pVngWpk9e4zhmENIycHbVy7gfdoVMbQeJWyioK23Wf1onJus2+b1uDc1f5RjtZgNLV+HV2G/VB82WMRIASFd1oL3cq3LOMUq+a42lxpzGvNlTd6IsXqfpdV/MTm85B5jn2bt3HS4xZCNBhaR9AFN1sLs5QHdAjtGYBq9HBC2cp6F54JkbkkzsWGIPIdOvIrOHeLb2eaB44WaGd3VPhbwAGeQ1GOpsnP3FwHO0HMekfG9eBsXWZvelkGeI9C5d4DaPTG01LnRd9MNYPMs6ytbLoInDuhcm6x3S+tY2vRG04Bc2Xye5UXlfvgli897n2tz5xWTz+M9r+TtI5AKTfkjbOG2MM/o4ym4Ds65V54XfwURP9OPuHeD5Vxdhrt2PZI4BZCNBi6+Ol4Lu5wfPgZCE1HFz2CyrmlXttWlg7/FM+MVv5EBB36AYo29543F8kIo3Fj9aBHoXNGYWQMLndzqsIZMSh5oILWmBBbgVJOdPAziKZa7V62nWut0YLNovPRRQ9C3viEnONaR9AbTgBznf1zcp9FuXuX/qgyz0AHJoBZQOV70E2bzG3x5ygFnn4oT7/EO8JzoWQMxFaBa2/ry4f/ZZLn3gFdkHL//bauYbZaCLFj8k8g+R/hIATss5fVFyuRi93CNTdE/ohvBwtg9Tw39z7DWJ8tBIW3WcVQEm49gGr/k21uQJeMRQenkVCApCoi86wtV2WFpleQ392ESGKqamXkoZq8B76TrK1iqgmp+48elCf9+t1mwUjYeKpVFS26AAKvwvojQaf4ghD9J+1nb2skcAshGg5dkuJ4IHU96/rgO8a+QpgyILIA2x5gAgNCUxNvzbkFVC72XwgqEVuMLnocgh+yVQOt4W8w/WVqdkf/oaJSoPjHo83EMqLK0RIj9x6M5t+gmn2cYpW7siqeeQ5JOKojizALbsfcOAyz+OXS1ec6PBsC/yN5bj9C6j3pdvPeDYMEbiFEw+Haw/64s+c2NeypjDxU4xdBNcZaJe20VjA3fhFUOl8wIpQPiMrRBtXsc8i6ykqEgqOC+w0SA7zGCmAxqp4mtZyiMhnfXDtbaw1S0QF0ycspTysjD5U3EVy7Y7VZgdEMsu9C5b2EUls+ow5+YeWOD7wD4W+g+DH0uoMwN5yJLryngs9l1QxP5EVtZcrZ+rTt/E0XQohKqJwR8dXEm3uNDisNZ+7I+myWPSMPa7X05jKeUWvxmOcwkvd2lxeNr/xOpIxcjKxhqLznSR243ZBxodVjrQ165ZbRDfcBYLQhdSiJQPDLCh+nXDtjNPkfqsVcVIv5GM2/w8g8DaW2jAxoHUUX3Io1UrH5C03IWjgX+Qmiv1fwBgMyr7CG5nFaJVkbP19pNbZtmQRuIUSDoVy7oJpMBt+pVi/NdwqqyWSUq1edt0VrjY6tQsfsF2bpgput2t+bF5VpP8RWg7kRnJ3KbGdKwWiW8pQy8iD7OpKHvRUYbVDZl4KzRSWfwAve41BZN0DmJYAb68tAGqvzY9YqbaUckPcGFX4RSaPymfUsb2KFsrKiS6nynvPNnD0xsi/EaP4NRsvfMZpORnn2r96zthGyqlwI0aAoZ/t672Hr8Dx0wXVW6lG0tZ2r0WMoZzu0WYwuuC2+7am8KAS/RLWYCaFp1nW6MPkyo2XCMLEdI/M8tGtPa3935E9rrtg3BJV5prXy23u8VafaLj0oAEFQTlTWhVYal6wr0aHpEP7F+mJR/KSVTzyJ23pffI5YRReilZEibboTVRN7pY0s0BXMpaei8lCNn9r6929jJHALIUQV6Nh6dP45idu9ovPRG89AN/0M1h1pbadKqQhC36C8R6BDsyAwrtx5NxhNMdf2B0crVOalKO8htk9S7t1ReS/Zv8a9L2ScFq/0ZRf0XBBZillwJ/hOtPaVFz1gZWXTMVJuP1NOcLTb8rO5GlItDHR0As/h9ueqQDlaoV3dITKfChfDAWScA2Y+eAaivEdsU2sfaooEbiGEqAIdeM+m92da861FT1QStAE0etPV6KxL4/m9y4tZ25nQYK5Fb7oKnTMiKeOa1qH4dqxCcO+DcrZPOK+UQjvaYw1/2wW7CER/sf4LTKA0fWqFBcecoPLQxU9YueB8J6KdPUkZ5H0nJ2Sc09qE8PcQXRQP6gdWOrJQKud+2HimVUwkZSMzMcpkd9M6hg5OQ4e+BiMP5TsR5Wyb3vu2YRK4hRCiKmLLsQ1UOgyRn9N8SAiKH031gnI/B6HoIbTvxNIFWzqyAL3xHKySlhow0RmnobJvLQ2UWgeg+EHSnxtOo0SoygBzBYSXW3eEfwZPv9T36vwtvzQL0BvPsAp+6LBVec1oBk0mxve9V9Ayswg2DY9vB6ygyEij+7bcoyPojedbtcHxAy50yVho9DjKO6Dyz7oN2/7GEIQQohZZq5Ht9h5HqjAPW8U62joUr6gV70VuPNuaG9d+rDnsEPjfhtBXZZqziIq3jFWVIz4fX7btIQh9Q8q95eFftnyEotHWIrPNedZ1CcSWowsqX6+g/RMgtpbkL0wK8ICjC6rxyxjeQVtOBT6whv9LE85EgCC64Hq0ruZCt22EBG4hhKgK7xHxilM2YkuonX9WtVVIhHjda11kc00A7S+TQc7Iq17pz5RSfSmJxP8rz2EtciOeAjbwns11UQh9mbpU6Gab09zavEM1HovR7JOk9Kc68AEpF+bZLhxsOCRwCyFEFSjlTlyclXDSA1nXk7w9Kt1ZSSdWwpbyNNo/wZoj9r+d+vbYljzkytkOXD2r8O6KVLJFzGhJcu/eBc7OmAV3WEPWtlXBAKLoNbthru2HWfyifQY8u7zvm+/VVnDWkYWYRY9hFj2KjvxO6s+tsf89bjhkjlsIIapAaw26OMXJiDV/mjkMHZkXXzjWB4KfWsUuzE2AAeYam5sV5IyC0A8Qeq/cuTAUP46O/VtBbnAgtgwdXQbRv9HFj8cLiritZytPPG93iCoP1Vc45K7ALB9YHVabix4h3UpgmOug+Cm0uRqVc3u5V2SluM9AEcIsfh6Kn8Xq0ZvokhdJ+RlVBtTDvv+aJD1uIYSoAh2YZC2wSqLAtQvK2RmlFIa7N4b3QAwjAyPjJIymH6OafQOeQ23u9UHmZSjf8RD5NsWbw/F83Db7vktF0QUPojddBdE/gBKsOV4NniOg2Vfgrk7ykVR7wQGygE0kDqXHsHrYVZ1LDoL/bbRZriqYZ2/sw5UbjQuKn2FLVjVd5v1lWWlnVeMXGvwWsYbdeiGEqGv+V0i5/Sn34Qpv1UWPQODd5BMZp6GyrgAdsDKrpWQXkMqJfEFyEZMoBN+F9QPBaFXx/VVWXHmbqkK5ILYs8VDGUJt1BS5wdUfF/iG9EQQFzaajXD1rqKH1RwK3EEJURfneYCl3Qn7t8rQOgv81bHuv0T+tbVzKCypzKxuYKohpa644WH4YfmtVddi9sseFwZG411o52qAajwNnV6xheJeVYKXxWKwZX2XzoGSqsuQtDYQEbiGEqApPf2znfI3ceCGLFGJrSRlgoksArCHczIuwz/2dXnCq3LYcvBzgGWS/r9vVHdXkA1Tzn1Et5qAaPQYqG7xpZmYzmloV2rYDsjhNCCGqQGVdhQ5+Fd+PHMbq/7hROfdUPHfqaFGanCyJs+uW52degEZDybNlVlM7rLnp8LcpHrAdMRKHxHX4J2uvd2wp4LZqdMf+i2eXc4PvJMi+DYpGYf1ZxEieWzcg67qELG4Nmap0/9w2oE+fPnrWrFn13QwhhADi+cr94yD8EzjaozLPR7l2qfQ+s/gpKB5L4nC5F5X3Gsq9e+I7tAnaj9YKZXispCFr+2EtONueuVHNv0YZeejIIvSGIVS8OM4D7n1RufdZe8LNgjKL1TZzgms3jCZvpXpI2nTkT3RgMhBBeY8E15618oVAKTVba21be1R63EIIUUXK0RSVfV3V78u8HK1yoWRMvLxnV1TOrUlBG+LD5ioLzCL0pmsh9CX1M8ztxKonXlMMKlzMpjzW1IE7L76tK8VCwFIhCP8A2o/KOA1d9BjJv09RiPyJjvyOcvWodsvN4peg+AmsHr1G+98G3/Ho7FtQSltV2eqABG4hhKgjSilU5tmQeXba9+j84fHUnTWZBa0KVFMgDHqj9UUi1R729B9orWzXQdCbSBr6L7s4LbqItFasKzfE/gFnh/g9Nr9XymGtVq9m4NaxVVD8OIlfJAIQeAsCb6MB7doNlTsa5excrXekq84Xpyml2imlvlJK/a6UWqCUuqqu2yCEEA2BjvwFkQXUW9AG0KutbWqN30I1esJKYLJVYlYSmcwLSc5x7gHPAShHfMuaqxdp5VvXYXDsVOYem8V9OgLObtVuNaEZ2C8QNLF6+DGIzENvONUqilKL6mNVeRS4TmvdA9gXuEwpVf2xCyHEDkOH52AWjMDcdAM6+JV9eswGSIemY248F3P9YMyix9Cbs6PF/rPqX1eZwvrnvabmXgNQPArc/cDVG/tV71V8XugrVOMxVnnP0naGIbYRHVkAgMq8MHVe+FJecO0BkZ/Q4XngOzX+5aJsePOAZ/+t6glrnFS+MFAD4fgceO2p86FyrfUqYFX810VKqT+ANsDvdd0WIUTDYRY/DcVj2JyyU4c+B/fBVpnGOlgtrGOr0f7XrR6wqwcq4yyUo+VWP9csHgPFT1O6mCr6NzrwPjSdDK5uVSgU4iZxNXUNf6mJzLfm3Ru/CIEP0CXjIbY4fjJepQt/+u/VfpRnH7R7HwispDQVa/QXq/xnk/dQzk6QNxFddB+E54CRDd5jILLQKqGqfIAPInPRhb9a9zu6Qt5rUPwUhOOVyzJORWVdVu2PrnU4niO+svl2rNGJ2D/Vflc66nWOWynVEdgD+NHm3HBgOED79u3LnxZC7EB0bCUUv0DCP5w6AOEZ1sIkz361+/7IX+iNp1pDsoQh/LNVajJvIspV/eFXM1YIxeUXU4XB3ID2v4GRdSnaeyQEp5CcDa0MlQV5b8CGk7GG1Wtjt5CB1lGUckHGyaiMk9E6BtGFgBNtBiH/5DSf5QXv0WgzP141rFxA1CF0yRhU7v0oVzdU3jjbp5gFd8Qz0YW3fOToH+B/BaPxU2m1RJvFaP+rEPgUDB8q4wzwHp/wZVD734bon2l+Nk+tZ2ertwQsSqksYBJwtdY6Kfmu1nqM1rqP1rpPs2bN6r6BQohtR+gbbP+50n508Itaf70uvCu+KGtzjzYMusQ6vjWKn8V+pXgYQtMBULmjIesKMJqkeIgD8t5ERZdYi7RqjYZQ4u+1Ug6UqwfK1dUqopIWD2CA/3V0/qUprolBZH7ljwpOJnnPdhgCH6XVEq1D1naz4hcg9hdEfkUX3okuHJF4YeADKvzilCCC9qSZFKaa6iVwK6VcWEH7Da21TeJeIYQoQ2WAbXITJxhbmyI0DZHZKY7/Unkt6YqEPkx9ztg8DG+g3H0h4zxw7g5s3nKkrKHijAswXN2sHN81NqdtJ4YOTqvg9JI0nqGwhtL9YK6M/77aDT8bCUlpUtKpiphE0vtzCXxstaP8SE5gMjq6pUSq9XubLhdKV5RvfuvV+VC5ssYfXgL+0Fo/WtfvF0I0QJ5DSJFyzKqoVduUN54pLfn4Vs2vmxtSn8s4G9MsgQ2nxOeSNaUlOt0HgKM5ynciyr03AFo1rqBudVmbF61VZ094BX09Z9f4sHlFnKS3Qt6Nyrwojcs2Z5MrO6+urIQsafy56PB3VqAuTzkgMgec1jStyjgVXbjA/trkp1ZQhrRm1EePux9wFjBAKTU3/t9R9dAOIUQDoYwsVKPnrAIcKiteiMMDObfX+p5ZAHxDSF5J7bHSbW6VCvpOjnaw/hhrCLf0S0sYCEF0ISpnVGnQNoueh/yzSG9hmBMyh1t5vjGsPdOOXUE1ptKQ4DvB9rCO/geRPyp5rye99jm7ofJesobfK6Fy7ojnH9+88txrle7MubPy9wA4WmH/Z6Cs3OabeY8Fz2Hx93jif/98JP+dcIGnP8rITe/91VQfq8q/pXbHc4QQ2yHl2Q+az4TQt0AE3PvX+j+Qpe/Ovg4dWwah761hUx2xenXZN1T7mTq2mtT1ql1W7W1zlf1pMx9dNBq8g9CqOZRUZfAyjMo4BZV9TWJ7tGnV8Q59ju3ohmqBcvdGh6aDNsG9D8rIxAzPg42nY59dzbA+i3NnyLoONl1ChT19ow1G0wqmD8o3ydkemk1F+ydB9Hdw7oLKOAllNErvft8p6JLXyrXdsL7UuPfF9L9rZUozV4OjDWRdhVJuMJqgPYdYuwH84621BToCrl6o3AfTbn91Sa5yIYRIk47+A9G/wdnJ6hGHpqFDX4ORh/KdbAWSdJ8V+ARdcD0pA55nEIQ+ruAJTsAFRuP4PG0VeAajcm5BORIXvJnh32HjEJKHs52Q8yAU3U5pv0vHrCQqJc/ZXI91XeYNGNkXbHl+4SjwTyTltirnrhhNa7rsaMV06Bt0wQ3xTG4mODugGj2DDv8AhaNIyivf6GFUmYpk2syHyCJwtEQ5O9RYuyRXuRBC1ADl7AjOjmgdRm882+rlaT/gRJe8is59GMOXbpnJPMCFbeB2tLfSd4YcpO6hRq3/zDT2FpcXmoxe9yU02bKdzQx8CAU3Jr9PNYHMS6DotuQ53pInK3iJTprzVtk3oqPLIDw9xT1bt/dc6xi65CWr7rkujo+K3GT9uaWgPAdAs++s/OjKi3K2txa2FZ9GcnGTILrokYTArYzG4Nlnq9pdVVKPWwghqkj737PKSpYuBosCQSi8yUrWkQ733lYN76SZQwc0egLlOxkrsFemusGuBF1wCzq2BjP0DRRcj+2XBL0Bih+0eqRV5WyV8KNSLsi+npRpTI3mVX8HYPrfxFx7CHpNLyh+FMw11mLC0DT0hpPi0xKpWdvaupYZMYmmXjgYW176Sx1ZiA68a5UercPRa+lxCyFEVQU/TLHCWEFknhWUK6GUA/LGo/MvgthqayUzQM4ojM0lQhs/id50A+gQFe8jTlXouxLR+eh1h7K52lVqaX4ZKS/4Paa5DpVxLsrVDa0jUHgHKb9shL/GLLwXlX1bylXhOrIIHfwclAPlPQJdcD9EZqRogAYdRPtfq+J6BKe1b95cn3zK0RatI+hNl0No5pZtikZLyHsd5WiafE8Nk8AthBBVlbJ8o4m1ZSvNxzg7QtPPILrY6iG6eliLnzaf9xwMzWeiS16B4odJHVwd4NobIr8DBWm/31KNofa0GBD7FQLz0YGPodET1ghF9E9Sfw4N/rfRzu6ojOQsbGbRk1AyFmtOXaGLn6by7WURK11qFSil0FlXQeF9JM1xZ1+LLhlrBW2CWz5K7F90wY2ovHQT0VSfDJULIUQVqYxT7YO3ygbXblV7llIo184od++EoL3lvBPlaE7FPeIoRP9EtfgJsu/C+vJQk0VGqsMs879BdME16IKRaew1D4L/laSjOrIoHrSDWEP6UdLbE+6wVrVXkZFxKuTcEU+Eo6xtc7mjUd4j4gvsyo+ARCH8I9rc2rKnabSt1t8ghBDbG8/A+N5uF1v+GVVWAQysohTaLKi5eU93v8qv0ZvQOoSReTqqxVxUkw8h534gr2bakLYUXxa0H0jKbm0vusyqkhZbu+X24OdUr7ypG5V5bsqzpv8DzHWHYK7eBXPd4ejg1NJzRsZJGM2/xmi5EKPZNAxfPOWIrmiUoprTClUggVsIIapIKYXKPC+eCnNzcNbgfwO9fjB6zV7otfuh1x2EGZha0aOsldDhWejQD+gUAUE5moKqrBKZhtC0ePviGcoK7wTKpt9UQC6orZiHNdqSXEc7obXVf3apMJS8hF4/CB3dnEq1quFKgaOLlczF2cn2CtM/CQpvh9gKIAaxf9CbrkMHv6z40d7DsJ1pdnRAGbX/RUkCtxBCVIMuHhvPlV22Vx2E2CKseeOolbij4Dp0+Gf7Z4Tnodf2Q+dfiN50KXrtvqmDRpM05k79b2x5dslLJPf+NFAEvuOoOPhWwFwbnw5IVZNbV+HZFQX5MOgidOHd1pW+I7FfluVIfo7nKFTzXzCafYJy226FthQ/SvKQdxBd9HDFrc66GoxmbMkb7waViWr0QIX31RQJ3EIIUR2RudgnTykviC5+Numo1gF0/nmgN1oL03SxVXFs0zXo2Iqk6w1nF8h7K56aNIXY0i2/jv6L/eptE/yvW/nfnbun0f7ywhD5FZy7pDhvgHcwVjDLIuXWLzyQfQ/kTQZvBaljwzMxS16zUttmX2vdV/a/nLuhyWRwHwLOHpB5A6rR/ahKis9oHbFfNQ4Q+9f+eJxyNEE1/QSybwbvcZB1CarpFFQV1zdUl6wqF0KI6nB2tmo/p7OPOrYs+VjwS+wXnMXQ/vdQ2ZcnnTHce6Cbf4te09P+XnMTWoetRW6efawEMbbzwiEITUM1/QBdMj5e07oq+7Q1xP6zP6W8qMyhkH01RP+wFmsV3FTu+YaVaSxjiLWC23UfOvgRKVe4F92PqXIwMs9Dew6D0JfWM7yHgVmA3jjUSjlKEKL/oAOvQ5N3kzLDJXJaSXDs9ms7Wlf6O6CMTFTm6cDplV5b06THLYQQ1aAyLyTtrV92c8q60EobmiQCelPqRykXW4ZoyzPAtLaDqYxz4sUwUtEQ+gaVcxtkDCF1zzjFe5w7Yd/30+DsjHI0RXkOwPANgtxRZYrDeMG5Myrv1dK92kqpeAGTVL+fUSi8ETMwBeVsi8o8B5V5FsrREl1wC+gitnwx8IO5Dl1ccf52pRRkXkny76UXsq5N8/ehfkjgFkKIalCu7qjGz1vpSSsLetE/kmtZu/fFvphHBspzYMXPc/VM0SivlbscUI5mqKbvg3PXFNc6QPlQyoWRczs0nwOeYyp+bykNObfEK4yV++yuvpTP+Gb4jkU1/wGVNw7V9AOMph+iHG0Sm5NzSwXD7/F3Ft6ELpPBTZvF8X3h5UUh+Hmln8LIPB1ybo3PVwNGK8gdZX3Z2IZJ4BZCiGpSnv1RTaeimv8U33qValFWCF38ROK9zs5WWVCVUeaoD1x7gLt/xe/Nvs7mXT7Iujq+ojx+naM1qvHT2C4k0xrK5Nw2DC9G40ch94UUb3UDGdazckdjuHa1kqqUXxgWno5ee4BVRaxsm5Ub5eqVtMJba40OTkPnX2qNQlS4YM2AMgv9rEV/KdYZ2OyJt31ixqkYzb9DtfgTo/kMDN+xad1Xn2SOWwghtoJSClQ2KuMkTOWIF+qwYTMnrHJGgucAtP9ta5+z7wSUbzBKVdynUu49IO9ldNFDViEPozlkXo6RMTj5YqM15IyCwhGgnFi9fI1q9IR9+cvQJ1h9unJz98oN2SNQ3gFb7gtMSr4OQK9H51+MzroBI2tYhZ9FF40G/wTSy+Cm2dyb12YhFFyd4jp3Qq10HfoeXfKilVrWvT8q60KUI3F7XWW/59sSCdxCCFFDlHcwunAUaJu0o84uydcrhWmut1ao6yKILkCbayHzopS5ukvvdfdBNXkr6bgOfYP2vwFmvjWHvnkBnWsf8B5mZWFz90MZGckPBQh/j30wDlhfLsrOm0d+t78WrOPFD2KahSjfEShXj+RHRpfHt7Clm7TECe69rF8GP0+dTE41Q2VdZrXC/1Zi6tLAMnTwQ2g6OSl4NxQN5yuGEEJs45QyIOtq7BY8qezrk643Ax9bNZ91PhC1toQVP2f1DqvBLH4qXvxiGkTmQPRXrFXlMYj8AMVPgnu/0qCtg1Mx1w/GXNMXc+N56Mj80jnyZDEoegi94dQt88yublQcRjT4n0NvOM2qxV1e5GcqHRrHa31ZUJmoxs/HF+dh/V6lKnmqV6MLR2KaxVA0msR849bvsy5ONSWw7ZPALYQQNcjIPMPaW+xoD3jAuSuq8Qsod9/ki4ufJHkbVgBKXkBrqyertbb2HFdCx9ZD8QspqpZBac7wwPvWT/630JuutxZ36U0Q/g694QzwDiL1qvUARBej/RMBUJkXkToRS1lBCLyNDv+SeFg1osLA7dgJsq9H5dyNavYtanNvG8DTn9QhLAaBT+Lb0OxEIfxtGu3eNkngFkIIG1oH0KGZ6PDc0iCaLiPjOIxmX2C0/A2j6Xsoz372F8ZWpXi5H62LMQvvR6/tjV7TE3PdUejQj6lfGpkdT8FaAR2A6EK0jkHRwyT2RLF+Dn4NmeeSemtWEIrHobVGubqh8sZiLVqrhA6iAx8lHvP0q6DSGlYilOJnwbVLUkIV5exibWNLeX8QQjPi2e1sGC0qb/M2SgK3EEKUY/onW+lHN12Gzj8Pve4gdGRhlZ6hdRCzcDTmmr0x1+yOmX8FOrYy8SLnTvY3G3lQeA/434z3oDXEFqPzL8QMfIoOfJzcHpWbRqt84OxuJR3RKRKuROdA8BPIuYeU29z0qtJV8sq9N6rp5DTerbbUrt58RLlReeMryAYXAp1vLXSzKdiiskegGqVYNQ+AYa3ST/oS4kNlXZBGm7dNEriFEKIMHV1srcDWgdI0pJhr0BvPSWvIuvQ5+cOthVe6wHpWaCp6/Ylos6j0GpV9I8nbupS1nzj4MXZ5tCm4Bl0wAr1hCOaGs9Cbh8bde8dTjKYaejbAyED5jgMjt4LrsDK9Fd5EyjlkTCh5ubSEpXK2B0fy4rtEbpQvedW7cnVDNZ8JjZ5PUfxEg7kOon8l36sUynMAeA6y/zzKA42eAndfwBOfK/dB9rVWrfMGSgK3EEKUof1vkzJNaPj79J4R+R3C80jc4mRaQ+CBSaVHlGc/VOMx4OxW9m6Izid1HnQTKAGCEJmNLhhtPUs5UHnjwGhj7Q1XmVg9ZhfgBM9BqCbvoIwslPJAxmmkNz+dgnImpHJVuaOw5sbLhxWn9Z7M81GuXvaPUgaGdwA4WqV4mUFFKVlV9tXxIfOy7/ZB9o0YjsYYeS+jmn2BynsD1fxHjMxzKvt02zTZDiaEEGWZG0jZ0zRttnnZiS4EpWy2KwUhNAtt5KBLXgWzGLwDwNE23qPcPJee7px6FIIT0eYNVkB2doZmX0J0AZgl4O6FspkD1jocz+2dTpGUFHQYjC3bqZR7D2jyDrrkBWvBm6MTOHdCGTngOSRlac0EvmOh6C+Sg7TLKiCSgnJ2sd5d9KS1tc7RGpV1Ccpz0JZrHC3A0XDntcuSwC2EEGUoz4B4aU1/4gkdiQ+5psHRMfW58FR0eBqlQdM/AeuLQqpNyZXR6JIX4tnU4glh4ilRtQ5ZC8JiK61j7v2soh6brrMWbqUcCq+MGzwDkop4KNfOqEYVl8RMaHlsLbr4aQhNByMbfGdYQ+7m39aecZyAE9XooYSMcHaUswuq8ZNV/ygNkARuIYQoy3s4+MdD5E9KV10rH2Sck37CDldvK3hHF5Hcq9XljlXU63WQVnANfATxwF36lugy9MbTrEVoOmjN9zq7onMftAJlWpnK7BjgG4zKucP2rNZRQKFU6vztWmu0uQE2HA/mJkprlxc/AO4BEPXHS5Sa1tx9HZXLbCgkcAshRBlKuSDvNbT/PWuBmJGFyjit8sIfCc9QkDfeyr8d+SnNu8qnGXWD5zAIfUalwdsmSOpN14K5kdKevPZD5A8ofs5KX6qrGrh91kpz72EYhs3we2wFetNN8aQqGo0bfENROddZc+pYQ/S66GHwv4X1pUiRMNKgAxD6uOxTITwTveFMaPpRg0pLWpskcAshRDlKuVGZp0LmqVW+V5sbITjd+sF7DETmkzTsnvxGMNqDudLqGeswuPdG5d6L3lQcTxaSKnh7wHdychuif5I8/B6C0PcpgrbDShoTW4ltb1wZKGcrlF3Q1gH0hpPL1bYOQ+BVa6Fek9fiQ/Q3QOgrtsxhpzM9ELV+X8I/gWffNK7f/kngFkKIGmL634fC27f0gHWU9OaRtRWcHB0h6zKUqyfK2c46lXsPesOp8frd8T3dwOb5X9y9UZnnp99IpcBoC7G/yx33ohqPRRc9ACGbkpgadGyN1UeOrbN6885OKOWG4KdQZptbguhciPyKdrSE0Jekn5e87LtNiP0DSOAGCdxCCFEjdGy1FbQJletIOrC2XW3uxXqtDGc6QOL8dhhii6H4cWj6WelR5WgJzaZaPdXoUrRzJ9AxlLkWXL3AtXtSQRJl5KGNZtaXgSQOiK0od0xB9kiUsx3afQCEviE5q1oJFFyHWXin1XblBhQ6++Z45bNUAdmE6O/W+c2jCSk5rOvL98SVKrdlbscmgVsIIWpC8LMUJxxWiUldAGYxyncUeI9BF9wMwY9IDFImmGsg8iu4dy89qpQbvEdYv06jKVrr1JnRzPJBG6sNwY8h43hUxmC0/+UUQ+Y6XjOb+HYyrCIpmedhZSezC8oOa7ubo1MF8+oOrDl+F8nTCm4r25urd4p7dzwy0y+EEDVBx6twJTFRjlYYjR7DyHsR5TvBWgBnbsJ+jtcoN1dcnbYE7EuLViSeTEUpH6rJJMi6KF4opbKvCgEr2YzRxP600dSqge1oCr7BJGeK80Hj8eA5BNvEN86dUI1fqbTM6Y5EArcQQtQE7yHYD2I6wDsw+bDnQJKDGNZQcooMYxXRph8dXWylIVVeK3uaLbsAaICrTA/fyMLIuhyVdXUFzynDXItq8i64Dy7zfAWufVBNJpZuDVM590DWxaDyADe4+lrn3XvFF63ZBO7YitS1w3dQMlQuhBA1QDm7oDPPhZJX2TJk7Lb2fzuT83gr38lo/3iIrdlyvfKB72yrd5omrTW6+DHrvcphLYjLGAIZF4D/+XJlPr3WF4aEOWxlLUzLuiT54c4uoCtbXOcCTz+Uowkqb4zVJrMQcCRX9FIOVNalkHVpuc9QwSK+CufEd0wSuIUQooYY2deivYehA9ZeZOU7KnV+biMTmryHLhlnreJWOajMs8FzeJXeqf3jwT8OCG4Zefe/A5nnQ+bFUPKiNbessqziGr5T0IF3oGSMtTLctQcq+0YrXWr5Nrq6o917QPgX7BO2OEFloTITK20pI6dKn0EpJ9q1p1WaNGH6wIiPTIiylF2ptFp/qVJHAk9grUgYq7UeXdH1ffr00bNmzaqTtgkhRENiru0P5trkEyoT1fwXrOImJVaATZHAROsAuuQlCEzGyow2BJV5Fkq50TqILnoUApOs3q+zG6WL1NwHoLIusvKAbyUdXRLf9hYGgtbog8pANXkXlbL4yPZLKTVba93H7lyd97iVNdnxDHAYsBz4WSk1WWv9e123RQghGjxzk/1x7Qei1kI4lboHrHUMveGMeJGTeK+6+Am0/3WrGInS4D0G1Ww6ysiu6daXUs6doNlUtP9dq0iLqyd4BoC5Dq28KCNVze4dT30sTusLLNZa/621DgMTgePqoR1CCNHwuXaxP+7oYAXtyoRmxJOxlB0KD1rbxvQ6MNeD/w30xtPic9G1RxmNMbKGoXLvh+hiWD8IvfFc9NoDMQtuq/X3NxT1EbjbAP+V+Xl5/FgCpdRwpdQspdSsdevW1VnjhBCiLmkdRgc+QRc/iw5+XuXgpLJvxVqdXna1uDdlEZCk90d+iffOKxKxkraEplWpbdWlS56DwPtYyWyKrf8NfIgufqpO3r+t22YXp2mtxwBjwJrjrufmCCFEhXRsLQQ/RpvFKE9/cPWudO+xjq1GbzgFdJEVPFUGGM2gyVtpDw1bdbAnWuUxI39Y+56zLke5e6d3v6MVGi/JNbDLN9aPjvyO8lZt8Vy1lIy3aU8Q/K9D9jW1//5tXH0E7hVAuzI/t40fE0KIBkkHp6E3XY21IjqMLhkL3kMh96EKK1rpgpFgrqN0K5QugVgIXTQalftA2u//f3v3HmNFfYZx/Pvs5bC7oCywGFuholargEVRiEpatTZeKkGNmmBaW9umUVupjdbGS/inl6St91ht0io2WlNjvRJLRUSrpq2AImgRbay1itFoo1huy7LL2z9mCMu6Fxb2zDBznk+y4ZyZPTvPb5fd98zvzPm9apyIRt22a+GbZsK663ei30cLqh+3a8cYrOhj3fNYT0TU/GIseUyVLwMOlnSApAowG5ifQw4zs90WsYn4+DKSM8TNJBVwU9JQo5+p5Ygu6HiGT75/uRPaF1Ytb0+qG4lG3wX1+5NMuQ8jecNP9+KoZG3yptOyCdU4qfftDYfVfNGGHAp3JC/gXAIsBFYD90XEqqxzmJkNiY6l9PqnNDYSmx7JPM6uUONk1PY4aluAxj6O2hZB41SSSdlGaDw8WeGsx4IqVcuz11ygme3f1zqgGe09N5Pj7+lyeY07IhYAC/I4tpnZ0Orv/KfvfVJ90omr41l2POtugGGnDFW4nSYJGrZPhWvMH4it64AY9IIqu52lMgXaHiDW/zp53b7xUDT8ItR4SKY59lR77MVpZmaFUJnex45m1HJWvw/VyB/3cnFaG9r7ykFFiM7XiQ13J+01K8egltlDUmyr+b7tAY/d8FnUen1ux9+TuXCbme0GaRi03kJ89N3kZeHoJGnlOQsqx/f/2Pp9YewT0L4Iut6EhoNh2JeQdv5Pc2x+hvjoEpIGHV3Q8Tyx8W4Y8zCq76NjlxWaC7eZ2W7SsBmwz9NJT+7YAJUZqPHQnXusKtB8+i4dN2Ir8fFV7PjWqXbY2kVsuM2vCZeUC7eZ2RBQXSu0zM72oF1r0gVKetoC7YvBhbuU3I/bzKyo6kb03XYz4wvKLDsu3GZmBaW60VDZ9rat7juaoeXruWSy6nPhNjMrMI28MbmoTc1Jz20q0HQWaj4772hWJX6N28yswFQ/BsY8DJ2roes9aJw0JP2xbc/lwm1mVnCSoHFi8mGl58JtZmZVEdEJmxcTHcuh/tOoedZOdz2zvrlwm5nZkIutG4gPz4Out9J+303E+ptg9F2o8XAiAjr+TnQsTS6ya56Z/GsDcuE2M7MhFxvugM5/k3RMA2iHgFh7GTHmT7D2QtjyYtKMhSZYfwOM+i2qTMszdiH4qnIzMxt67fPZXrS76XoPNtwJHcvTM3FIivpGYu33k3an1i8XbjMzq4J+JnQ3PwZs+uT2aIfOV6qWqCxcuM3MbOg1nws09dgoaDgQ1HP7NoHL0sD8HTIzsyGn4edDZRrQDFRAw6FuDGq9GbWcm27v+aAR0HBYxkmLxxenmZnZkJMqMOp22PISbFkB9fvCsBORKkT9eGh/CjqegdhCUtjrUOutSD6fHIgLt5mZVYUkqExJPnbYXo9G3UJ0rIQty0CjoelkVDcip6TF4sJtZma5UC9F3QbmOQkzM7MCceE2MzMrEBduMzOzAnHhNjMzKxAXbjMzswJx4TYzMysQF24zM7MCceE2MzMrEBduMzOzAnHhNjMzKxBFRN4ZBiTpA+A/g3hIG/DfKsXZ03nstadWxw21O/ZaHTfUztj3j4ixve0oROEeLEnPR8TReefIg8dee2Ov1XFD7Y69VscNtT32bTxVbmZmViAu3GZmZgVS1sL9m7wD5Mhjrz21Om6o3bHX6rihtscOlPQ1bjMzs7Iq6xm3mZlZKZW6cEuaI+lVSask/TLvPFmTdLmkkNSWd5YsSLo2/Xm/JOkhSa15Z6o2SadKek3S65KuzDtPFiSNl/SUpFfS3+1L886UNUn1kl6U9GjeWbIkqVXS/env+WpJx+adKQ+lLdySTgTOAKZExCTgupwjZUrSeOBk4K28s2RoETA5Ij4P/BO4Kuc8VSWpHrgVOA2YCJwnaWK+qTLRCVweEROBY4Dv1ci4u7sUWJ13iBzcDDwWEYcCU6jN70F5CzdwMfDziNgMEBHv55wnazcCPwJq5iKGiHg8IjrTu88B4/LMk4HpwOsR8UZEdAD3kjxZLbWIeDcilqe315H88d4v31TZkTQOOB24Pe8sWZI0EvgicAdARHRExNpcQ+WkzIX7EOALkpZIelrStLwDZUXSGcA7EbEy7yw5+hbw57xDVNl+wNvd7q+hhgoYgKQJwJHAkpyjZOkmkiflW3POkbUDgA+AO9OXCW6XNDzvUHloyDvA7pD0BLBvL7uuIRnbaJKptGnAfZIOjJJcRj/A2K8mmSYvnf7GHRGPpJ9zDcl06j1ZZrNsSRoBPAD8ICL+l3eeLEiaCbwfES9IOiHnOFlrAKYCcyJiiaSbgSuBufnGyl6hC3dEfLmvfZIuBh5MC/VSSVtJ1rj9IKt81dTX2CUdTvLMdKUkSKaLl0uaHhHvZRixKvr7mQNIugCYCZxUlidp/XgHGN/t/rh0W+lJaiQp2vdExIN558nQDGCWpK8ATcDekn4fEV/LOVcW1gBrImLb7Mr9JIW75pR5qvxh4EQASYcAFWpgYfqIeDki9omICRExgeQ/+9QyFO2BSDqVZApxVkRszDtPBpYBB0s6QFIFmA3MzzlT1Sl5RnoHsDoibsg7T5Yi4qqIGJf+bs8GnqyRok36N+xtSZ9LN50EvJJjpNwU+ox7APOAeZL+AXQA36iBM7Ba9ytgGLAonW14LiIuyjdS9UREp6RLgIVAPTAvIlblHCsLM4DzgZclrUi3XR0RC/KLZBmZA9yTPlF9A/hmznly4ZXTzMzMCqTMU+VmZmal48JtZmZWIC7cZmZmBeLCbWZmViAu3GZmZgXiwm1WApK6JK3o9jFhF77GmdVs1iHpMUlra62jldlQK/P7uM1qyaaIOGI3v8aZwKMMYlELSQ3dGrsM5FqgBbhw8NHMbBufcZuVlKSj0gY7L0haKOlT6fbvSFomaaWkByS1SDoOmAVcm56xHyTpL5KOTh/TJunN9PYFkuZLehJYLGm4pHmSlqbNH3rtUBYRi4F1mQzerMRcuM3KobnbNPlD6VretwDnRMRRJCsJ/iz93AcjYlpEbOtn/O2I+BvJcqlXRMQREfGvAY43Nf3ax5M0tnkyIqaTLDN8ba12bTLLgqfKzcphh6lySZOByWxf/rUeeDfdPVnST4FWYATJkqmDtSgiPkxvn0zS+OKH6f0m4DMkTwrMbIi5cJuVk4BVEXFsL/t+B5wZESvTbmon9PE1Otk+K9fUY9+GHsc6OyJe2+W0ZrbTPFVuVk6vAWMlHQtJG0xJk9J9ewHvptPpX+32mHXpvm3eBI5Kb5/Tz7EWAnPSrl1IOnL345tZX1y4zUooIjpIiu0vJK0EVgDHpbvnAkuAvwKvdnvYvcAV6QVmBwHXARdLepGkl31ffgI0Ai9JWpXe/wRJzwJ/BE6StEbSKbs6PrNa5u5gZmZmBeIzbjMzswJx4TYzMysQF24zM7MCceE2MzMrEBduMzOzAnHhNjMzKxAXbjMzswJx4TYzMyuQ/wPwVSqgTYIGdAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 576x432 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "from sklearn.datasets import make_blobs\n",
    "from sklearn.model_selection import train_test_split\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "data_features, data_targets = make_blobs(n_samples=600, centers=2, n_features=2, random_state=42)\n",
    "\n",
    "# The function outputs targets 0 and 1 so we need to convert targets 0 to -1\n",
    "transformed_data_targets = [-1 if t == 0 else +1 for t in data_targets]\n",
    "        \n",
    "# Plot the dataset\n",
    "plt.figure(figsize=(8, 6))\n",
    "plt.scatter(data_features[:, 0], data_features[:, 1], c = transformed_data_targets)\n",
    "plt.title(\"Toy dataset\")\n",
    "plt.ylabel(\"Feature 2\")\n",
    "plt.xlabel(\"Feature 1\")\n",
    "plt.show()\n",
    "\n",
    "# Split data into training and test set\n",
    "features_train, features_test, labels_train, labels_test = train_test_split(data_features, \n",
    "                                                                              transformed_data_targets, \n",
    "                                                                              test_size = 0.3)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 5.2 SVM class definition <a class=\"anchor\" id=\"primal-svm-class\"></a>\n",
    "\n",
    "Next, we would like to implement an SVM class. We will use the knowledge we already acquired:\n",
    "\n",
    "1. Our objective function using the hinge loss function is given by: \n",
    "$$\n",
    "J(\\mathbf{w}) = \\frac{1}{2}\\|\\mathbf{w}\\|^{2} + C \\frac{1}{N} \\sum_{n=1}^{N} \\max \\left\\{0,1-y_{n}\\left(\\left\\langle\\mathbf{w}, \\mathbf{x}_{n}\\right\\rangle\\right)\\right\\}\n",
    "$$\n",
    "2. We can minimize this function by computing the gradient: \n",
    "$$\n",
    "\\nabla_{w} J(\\mathbf{w}) = \\frac{1}{N} \\sum_{n=1}^N \\left\\{\\begin{array}{ll}\n",
    "\\mathbf{w} & \\text{if} \\max \\left(0,1-y_{n} \\left(\\langle \\mathbf{w}, \\mathbf{x}_{n} \\rangle \\right)\\right)=0 \\\\\n",
    "\\mathbf{w}-C  y_{n} \\mathbf{x}_{n} & \\text { otherwise }\n",
    "\\end{array}\\right.\n",
    "$$\n",
    "3. Given the gradient we use stochastic gradient descent to train our model\n",
    "4. After training our model we can make predictions using the [sign function](https://en.wikipedia.org/wiki/Sign_function)\n",
    "\n",
    "As mentioned previously, we will assume that the bias $b$ is contained in our weight vector as the first entry $w_0$, that is $\\mathbf{w} = [b, w_1, ..., w_D] = \\mathbf{w} = [w_0, w_1, ..., w_D]$."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "from sklearn.utils import shuffle\n",
    "\n",
    "class LinearSVM:\n",
    "    \n",
    "    def __init__(self, regularization_param):\n",
    "        \"\"\"\n",
    "        Initialize the model by setting the regularization parameter \n",
    "        and a boolean variable for our trained weights.\n",
    "        \"\"\"\n",
    "        self.regularization_param = regularization_param\n",
    "        self.trained_weights = None\n",
    "    \n",
    "    def add_bias_term(self, features):\n",
    "        \"\"\"\n",
    "        Add intercept 1 to each training example for bias b\n",
    "        \"\"\"\n",
    "        n_samples = features.shape[0]\n",
    "        ones = np.ones((n_samples, 1))\n",
    "        return np.concatenate((ones, features), axis=1)\n",
    "    \n",
    "    def compute_cost(self, weights, features, labels) -> float:\n",
    "        \"\"\"\n",
    "        Compute the value of the cost function\n",
    "        \"\"\"\n",
    "        n_samples = features.shape[0]\n",
    "        \n",
    "        # Compute hinge loss \n",
    "        predictions = np.dot(features, weights).flatten()\n",
    "        distances = 1 - labels * predictions\n",
    "        hinge_losses = np.maximum(0, distances)\n",
    "\n",
    "        # Compute sum of the individual hinge losses\n",
    "        sum_hinge_loss = np.sum(hinge_losses) / n_samples\n",
    "\n",
    "        # Compute entire cost\n",
    "        cost = (1 / 2) * np.dot(weights.T, weights) + self.regularization_param * sum_hinge_loss\n",
    "        \n",
    "        return float(cost)\n",
    "    \n",
    "    def compute_gradient(self, weights, features, labels) -> np.ndarray:\n",
    "        \"\"\"\n",
    "        Compute the gradient, needed for training\n",
    "        \"\"\"\n",
    "        predictions = np.dot(features, weights)\n",
    "        distances = 1 - labels * predictions\n",
    "        n_samples, n_feat = features.shape\n",
    "        sub_gradients = np.zeros((1, n_feat))\n",
    "\n",
    "        for idx, dist in enumerate(distances):\n",
    "            if max(0, dist) == 0:\n",
    "                sub_gradients += weights.T\n",
    "            else:\n",
    "                sub_grad = weights.T - (self.regularization_param * features[idx] * labels[idx])\n",
    "                sub_gradients += sub_grad\n",
    "                            \n",
    "        # Sum up and divide by the number of samples\n",
    "        avg_gradient = sum(sub_gradients) / len(labels)\n",
    "\n",
    "        return avg_gradient\n",
    "    \n",
    "    def train(self, train_features, train_labels, n_epochs, learning_rate=0.01, batch_size=1):\n",
    "        \"\"\"\n",
    "        Train the model with stochastic gradient descent using the\n",
    "        specified number of epochs, learning rate and batch size.\n",
    "        \"\"\"\n",
    "        # Add bias term to features\n",
    "        train_features = self.add_bias_term(train_features)\n",
    "        \n",
    "        # Initalize weight vector\n",
    "        n_samples, n_feat = train_features.shape\n",
    "        weights = np.zeros(n_feat)[:, np.newaxis]\n",
    "        \n",
    "        # Train the model for a certain number of epochs\n",
    "        for epoch in range(n_epochs):\n",
    "            features, labels = shuffle(train_features, train_labels)\n",
    "            features, labels = train_features, train_labels\n",
    "            start, end = 0, batch_size\n",
    "            while end <= len(labels): # Training loop over the dataset\n",
    "                batch = features[start:end]\n",
    "                batch_labels = labels[start:end]\n",
    "                \n",
    "                grad = self.compute_gradient(weights, batch, batch_labels)\n",
    "                update = (learning_rate * grad)[:, np.newaxis]\n",
    "                weights = weights - update\n",
    "                start, end = end, end + batch_size\n",
    "                \n",
    "            current_cost = self.compute_cost(weights, features, labels)\n",
    "            print(f\"Epoch {epoch + 1}, cost: {current_cost}\")\n",
    "                \n",
    "        # Set the trained weights to allow making predictions\n",
    "        self.trained_weights = weights\n",
    "\n",
    "    def predict(self, test_features) -> np.ndarray:\n",
    "        \"\"\"\n",
    "        Predict labels for new test features.\n",
    "        Raises ValueError if model has not been trained yet.\n",
    "        \"\"\"\n",
    "        test_features = self.add_bias_term(test_features)\n",
    "        if self.trained_weights is None:\n",
    "            raise ValueError(\"You haven't trained the SVM yet!\")\n",
    "            \n",
    "        predicted_labels = []\n",
    "        n_samples = test_features.shape[0]\n",
    "        for idx in range(n_samples):\n",
    "            prediction = np.sign(np.dot(self.trained_weights.T, test_features[idx]))\n",
    "            predicted_labels.append(prediction)\n",
    "            \n",
    "        return np.array(predicted_labels)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Compute some values to make sure the cost is computed correctly\n",
    "# I calculated the values for this example by hand first\n",
    "svm = LinearSVM(regularization_param=1)\n",
    "weights = np.array([1, 2])[:, np.newaxis]\n",
    "features = np.array([[0.5], [2.5]])\n",
    "new_features = svm.add_bias_term(features)\n",
    "\n",
    "labels = np.array([-1, +1])\n",
    "assert svm.compute_cost(weights, new_features, labels) == 4.\n",
    "gradient = svm.compute_gradient(weights, new_features, labels[:, np.newaxis])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 5.3 Training and testing an SVM <a class=\"anchor\" id=\"train-test-svm\"></a>\n",
    "\n",
    "After defining our SVM class we can train a model and test it on unseen examples."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Epoch 1, cost: 27.249462152357992\n",
      "Epoch 2, cost: 7.348722759212155\n",
      "Epoch 3, cost: 3.3972581191442885\n",
      "Epoch 4, cost: 2.27870883462286\n",
      "Epoch 5, cost: 1.8875111235175854\n",
      "Epoch 6, cost: 1.6020090021589712\n",
      "Epoch 7, cost: 1.413202164641454\n",
      "Epoch 8, cost: 1.2452165334256207\n",
      "Epoch 9, cost: 1.1141080363838292\n",
      "Epoch 10, cost: 1.0118493571787424\n"
     ]
    }
   ],
   "source": [
    "# Initialize a new SVM and train it on the given toy dataset\n",
    "regularization_param = 100\n",
    "lr = 0.000001\n",
    "svm = LinearSVM(regularization_param)\n",
    "trained_weights = svm.train(features_train, labels_train, n_epochs=10, learning_rate=lr)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Accuracy on test dataset: 1.0\n",
      "Recall on test dataset: 1.0\n",
      "Precision on test dataset: 1.0\n"
     ]
    }
   ],
   "source": [
    "# Predict lables for unknown test samples\n",
    "from sklearn.metrics import accuracy_score, recall_score, precision_score\n",
    "\n",
    "predicted_labels = svm.predict(features_test)\n",
    "predicted_labels = predicted_labels.flatten()\n",
    "\n",
    "print(\"Accuracy on test dataset: {}\".format(accuracy_score(labels_test, predicted_labels)))\n",
    "print(\"Recall on test dataset: {}\".format(recall_score(labels_test, predicted_labels)))\n",
    "print(\"Precision on test dataset: {}\".format(precision_score(labels_test, predicted_labels)))    "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 5.4 Visualizing the decision boundary <a class=\"anchor\" id=\"decision-boundary\"></a>\n",
    "\n",
    "Given our trained model we can visualize the decision boundary, as done below."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbgAAAFNCAYAAABylXoZAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Z1A+gAAAACXBIWXMAAAsTAAALEwEAmpwYAABrfElEQVR4nO2dd5xcZdWAn3NnZmf7poeQQkJJKIFQAtJb6L1JlSYC+iGiogioFBUEUcAGioAgvSNVOkqREqrU9JC66cn23Zl7vj9mZjM7e6fuvXtndt6H3/7YKfe8Z25m77nnvKeIqmIwGAwGw0DD8lsBg8FgMBi8wBg4g8FgMAxIjIEzGAwGw4DEGDiDwWAwDEiMgTMYDAbDgMQYOIPBYDAMSIyBMxhyREReFZFv9VHGpSJyaw7ve1ZETu/LWmnk3iEiv3JbbgF67C0iC/3WwzCwMQbOUPSIyO4i8qaIrBWRVSLyhojsKCI7i0iLiNQ6HPOBiHxXRMaLiIrIBymvDxORThGZ128fBFDVq1U1q5FU1YNV9c7+0MlgGKgYA2coakSkHngK+CMwBBgNXAl0qOpbwELguJRjJgNbAvclPV0dfz7BycBcD1U3eISIBP3WwVAaGANnKHYmAqjqfaoaVdU2VX1eVT+Ov34ncFrKMacBz6jqyqTn7gJOT3nPPzItLCL7i8gXcc/xT4CkvP5NEflcRFaLyHMislHSa1uJyAtxj7NRRC6NP3+FiNwd/71SRO4WkZUiskZE3hWRkfHXusOhImKJyM9EZL6ILBORf4hIQ/y1hId6uoh8JSIrROSnWc7psLhuTSLy7xS9d43rsTb+/12TXpsnIvslPU7+LBn1EJGqeHh0tYh8BuyYci4vFpHZcZ0+E5Gjk147I+613yAiK4FfxM/r1knvGSEirSIyPMtnN5QRxsAZip0ZQFRE7hSRg0VkcMrrdwF7ishYiBkDYt5ZanjvbuBEEQmIyJZALfB2ukVFZBjwKPAzYBgwG9gt6fUjgUuBY4DhwGvEPUYRqQNeBP4FbAhsCrzksMzpQAMwFhgKfBtoc3jfGfGffYCN47r/KeU9uwOTgGnAZSKyRbrPBpwC/DL+uT4E7onrPQR4GvhDXJ/rgadFZGgGWamk0+NyYJP4z4H0vNmA2Pndg9j5uBK4W0RGJb3+NWAOMDKu+/3AN5JePwl4SVWX56GrYYBjDJyhqFHVdcQumgr8DVguIk8kPB1VXQC8CpwaP2QaECZ2oU5mIfAlsB8x7+2uLEsfAnyqqg+rahdwI7A06fVvA79W1c9VNQJcDWwb94YOA5aq6u9UtV1Vm1TVyZh2ETMkm8a90/finzeVU4DrVXWOqjYDlxAz1smhuivj3u1HwEfAlAyf7WlV/Y+qdgA/BXaJ3yAcCsxU1btUNaKq9wFfAIdnOVfJpNPjeOAqVV0V/zf7Q/JBqvqQqi5WVVtVHwBmAjslvWWxqv4xrlcbsRuYk0Qk4VWfSvZ/U0OZYQycoeiJG5EzVHUMMJmYV3Rj0lvuZL2BOxW4P26UUvkHMU/oJLJfDDcEFiTpoMmPgY2A38dDi2uAVcRCmKOJeWSzc/hodwHPAfeLyGIR+Y2IhNLoMj/p8XwgSMybSZBsfFuJeXnpSP5czXHdN3RYJ7HW6CyfI5l0evQ4n6nriMhpIvJh0vmcTMzD7KVzXO+34/L3FpHNiXnJT+Shp6EMMAbOUFKo6hfAHcQugAkeBcaIyD7EQobpsg8fIealzFHVr7IstYSYoQIg7imMTXp9AXCuqg5K+qlS1Tfjr22cw2fpUtUrVXVLYFdinl/qfiLAYmIGNcE4IAI0ZlsjDcmfq5ZY8s5ih3USay2K/94CVCe9tkEea/Y4n3G5CR02IuadfxcYqqqDgE/ouefpNPbkTmJhylOBh1W1PQ99DGWAMXCGokZENheRC0VkTPzxWGIe2FuJ96hqC/Aw8HdgvqpOd5IVf9++QC61bE8DW4nIMfFQ4PfoeUH/C3CJiGwV16tBRL4ef+0pYJSIfF9EwiJSJyJfc/hs+4jI1iISANYRC1naDrrcB/xARCbEDdLVwAPx0GghHCKx0osKYvtZb8XDhs8AE0XkZBEJisgJxLJRn4of9yGx0GhIRKaSkr2ahQeJna/B8X/L85NeqyFmwJYDiMiZ9LyBScfdwNHEjFzGhCFDeWIMnKHYaSKWYPC2iLQQM2yfABemvO9OYt5Hxgudqk5X1azhQ1VdAXwduAZYCWwGvJH0+mPAtcTCi+viOh0cf60J2J/Y3tVSYvtJ+zgsswExw7wO+Bz4N86h09vjz/+HWGlDOz0NRL7cSyzpYxWwA/FkjXjW6WHEzu1K4CLgsPi5APg5sSSR1cQSQe7NY80riYUl5wLPk/Q5VfUz4HfAf4l5pVuTdK7TETfK7xMzjq/loYuhTBAz8NRgMJQqInI7sQSUn/mti6H4MAWTBoOhJBGR8cT2XLfzWRVDkWJClAaDoeQQkV8SCwtfp6qmI43BEROiNBgMBsOAxHhwBoPBYBiQGANnMBgMhgFJSSWZDBs2TMePH++3GgaDwWAoIt57770Vqtqr0XZJGbjx48czfbpjDa/BYDAYyhQRSW0xB5gQpcFgMBgGKMbAGQwGg2FAYgycwWAwGAYkJbUH50RXVxcLFy6kvd00EjcUD5WVlYwZM4ZQyGn6jcFg6A9K3sAtXLiQuro6xo8fz/rZhwaDf6gqK1euZOHChUyYMMFvdQyGsqXkQ5Tt7e0MHTrUGDdD0SAiDB061EQVDAafKXkDB5SEcYtGokS6Mo/vUlXsqI2qoqpEI1HsqNN4sNzXtO3cjldVujq7iEaimWVGo0Q6IyS3eHND12LAjtp0pXy2fEg+D258J5tWN9Pa1NZnORDTbcWilTStbgbAtm1WLV1NZ3unK/JTaW/tYN3KpoLPZTZy/V7nQq46qnai0UZUs58zVY29116bpy6daHRxTmvkpIO9CrVb4o+jOA+67zuqXWjkK9RudnjNRu0Wz74LmfA8RBkfZ3EYsExVJ6e8diHwW2B40swpz1BV1q1oYuWS1XFjIwSCFhXhEIFQgJZ1rWg06R9B6J4jHAgFGDSigZr6KsLVYaKRKCsWraK9pYPlKxq5/JrLmTlrBrZtc+ghh3LttdfS2dJFe1sHzataug3H8lXLufYPV/Pggw+xYuFK2ls6EEtixi2y/o/27AvO4ne/uoH6unoQIRgKUDeklsqaMIFgACsQuzdpXddG85oWIp0RQuEQNfVV/OKXvyAcquSsU2NzPa2gxaARDVTXVREKh+jq6GLJnEaiXTGdghXBJOMKofD6x+GqCoaPGYodVRbNXIKqcuI3v879tz/EwsUL+fCTDzjikCOxo8r/PvuYfz7zOD/70c/jJxxClSGGjR5CZ1snqxrXolEbsYS6IbUM3XAI7S3tNK9uoa25nWgkihWwEEsQEWoaqhk2ekj351yzbC3tLR2oKoFggGAoQCAUIBQOEaoIUlVXRfOaFppWNhONRrv/LavrqxgxbhihcIhoJMqaZWtpXtNKIGjRMLyeUWM24KM3P+k+H4lzEggFqAiHqKmvJhKJIgJiCXZUiXRFaG9uxwpY1A+tY+WS1XS1r794iCUsW76CC7c/ibGTNuS0K0+geVULVsDia4duT2VtJc/87UX+dfvLtDW1M2bShtQNqSEUChKuCfPc31+hvaUDBMZM3JAzf3kCn789i8b5y9n96J2oG1TLIzc+yfwvFlE/uI69T9yVo757MFW1Vb2+9x+9+im/OfPPrGlcQzRqM2rCCNauaKK9pZ1oxKZhWB0VVRWsW9FERVUF+56yOzX11Tz1l+dpWtVMVW0Vka6Y4R8xbjjhqgo62jvZbLuNOfEnRzJszFBmvjeHwRsMYviYoVx/9s28+c93ica/z+HqCsSysCPR7huImoYaKmsqWLl4NWjsO1oRDlEzqJqvHboDcz6cx4z35sQMWfw7WTeklorKCtqa21m7fB3V9VVsuctE9jp+VyZsPY7O9i6GbDCYpfMaWTRjKWM335BQRZDrz/4Li2YvpbI6zHHnj+Tk895AaANrGNjLgDZUhsS+sIFhUH0GUnUcIoJGV6DRJdDxCrT+HTQKYqE130Jqzou9R6No6z3QchfYzRCcCNE5oGsAGw1NhcqjoOtViCyA4GZIzSlIaJuYEYwuhehaaPkjRD7uvuisvxJVxP5nDQdrNFRMQioPRK2NoOla6HgBsMEaBfWXYVXuiXa+i669FKKLY7KsIWCvjOkT3BxqL4DAxtD1DgQ2gtAUJDofrPrYOkQQCceMbcs/oHM6aASC44AARBdCYDRSfRzatQBafgvaBbSnjF4fRGwQfBQIoERBqqDyUKg+HQlugoh3fpbnzZZFZE+gGfhHsoGLT2a+Fdgc2CEXAzd16lRNLfT+/PPP2WKLLbLqoRq7QLeu69sdceLim/BWVJWvn3EsJx17MscecRxRO8rPr/opDfWD+MkFF/c4NhKJEAx6v+35x1t+T3VVTbeB84q333uL2+++jb/e8DdP13EDK2AxeuKomGHvjPa4m9xuz2344D8fu77m4hWLuOqQP7ouNx1iwb4n78nxPzqCjbfZKKbD7KWcM+VHdLR2eLJmIBhALCFcVUGkKwooHW2dUAQ93IPhAJGO1IiE8s2fLuaE87JdbipjF/7oXGJ3uqnnLwTVZyFVh6DrboKu58jvQ1fE5c8H3PaiLZwHw2ejitjnTBxbS0y3TPpZxM5P5shP2mOtIUjD1Uh47wKOX4+IvKeqU1Of9/xqq6r/ic9tSuUGYhOD/+m1DgBL5y2jdV0bdUOeYNjo3xKsWEKkcxQrFv2IplVH5CxHbUWTvshvvftfwuEwxx5xHAABK8ClP/gp047ah++dewHPvvAMz7/yHK1trdjRKNdccR3f/sHZPPXAs7S1t3HxlRcxc/YMJmy0McuWL+Oyi65g6y23Zt8j9uLhfzxGa2srZ19wFjtM2YEPPn6fkSNGctNv/0plZSUPPnY/Dzz2AF2RLjYasxG/+cVvqarsfQef4OIrLqIiHOaTz/9HS0szF3//UvbZY186Ojq44prL+OTz/xEIBLj4B5ey89RdmDl7Bpf84mK6urqw1eaP1/6Z8ePGdxuE3/3pOmbPnc2RJx/O0YcdzRaTtuT2u2/j5t/9lf2O2ofH73ky5oECBxwzjXv/9gCWJVz+68tYvDR2Z3nphT9jhyk79NAz3br/96Nvs7RxCR0dHZx24hmccMyJQMxAnXjsyfznjVcZPmwEPzzvQq77w7UsXrqES3/4U6bttR8PP/4QL7z6PM3NTTQub+SIg4/ku2d/r9c5uvWuv/HsC8/Q2dXJ/nvvz/fO/T6tba18/5LvsXTZUuxolP8767sccsChOX9n+gu14aW7/8OrD77BBhsNZ7tp27BuZRMd7d4YN6A7MhHpzBx+94Pexg3GbtrBkBERVCFzFLkdol9meL0LWv+Ctt4KFPLZOyE6s4DjcqHQ8G3qzX/vcKN7a8WPtVegq78HQx9BQpv1QZYzvmRRisiRwCJV/SjbXoWInAOcAzBu3LiC1uto66BpZTN1Q55g5EaXYgVim/+h8GJGbnQpQF5GLpmZc2ay1eZb9XiutraOUSNHMX9BrHvMZ19+yhP3Ps2ghkEsXLyw+333PnQPDXUNPPPgc8yYNYOjvnG44xrzF8zj+l/dwK9+djUXXHI+z738L4485Cj23+dAjj86dpG/4ebrefifD3HqCadl1HfR4oU8fMejfLXwK077zinsutNu3PPQXSDCk/c/w+x5sznru2fw3CMvcv+j93HaiadzxMFH0tnV2WuP7cLv/riHB/f2e28BYFkW++65Hy+88jzHHnEcH33yIRtuMJphQ4dx4c9+wOknn8nUbaeyeOlizjr/TJ596LkectOte/XPr2FQwyDa29s57vSjOWDfAxk8aDCtba3sPHUXfnLBxZz34+9w483Xc/uf72T2nFn85IqLmLbXfgD879OPefL+Z6iqrOK4049mr932Yestt+5e9/W3XmP+V/N4+M5HUVW+c+G5vPv+O6xas4oRw0Zwy423AtDU3JTxHPtNtDPKoplLWTRzqd+qFBVb7NDCrx+YTSikWYxbPhSfYS892tGW25FBv3Zdcr8bOBGpBi4FDsjl/ap6C3ALxEKUhay5cvFqAIaN/m23cUtgBdoZNvq3BRu4XNhtp90Z1DCo1/PvfTSd0048HYCJm05k0qaTHI8fs+EYtpi0JQBbbT6ZRUsWATFP58a/3EBT0zpa2lrZfec9supy8P6HYFkW48eNZ+zoccyZN5v3PnqPbxwfM4ybjN+EDUeNZu5Xc9l26+34y+03sXTZUg7Y50DGjxuf82c+ZP9D+fOtf+TYI47j6eef4pD9DwHgzXfeYNacWd3va25ppqW1hZrqmu7n0q171wN38sKrLwCwpHEJ8xfMY/CgwYRCIfbcdc/YedxkEhUVFYSCISZuOolFS9bfUOz6td0YPGgwAPvvcyDvfTi9h4F7463XeePt1znqlNh3obWthXkL5jF12x259sZfc90ff8M+u+/D1O12zPk8GIqH865eRFV1EcRODb3peMUTsX54cJsAE4CE9zYGeF9EdlJVT24525pjRi1YscTx9XTP58KmEzbluZf/1eO55uYmljQuYaOxG/HZF59SVZU+bJgLFaGK7t8DVoCOaOzzXPyLn3DTdTez+cQtePTJR3jnvbezyhJ63rpm8qAPP+gIpkyewquvv8o53z+LKy/5FbvsuEtOOm+3zXZ8tXA+q1av5MV/v8h3vnkeALatPPj3hwmHw3mta1nCm++8yQO3P0RVZRWnnnsyHZ2x0FsoGOr+HJYl3efLsiyi0fVhqtTPmvpYVTnnjG9z4jEn9dLp0bv+yb/feJUbb76enXfcle+efX5O58FQHIgom2zlTkaqwQN0tSdi+71MQFX/p6ojVHW8qo4HFgLbe2XcILZvBhDpHOX4errnc2GXnXalrb2Nx59+DIil0V9z4685+tBjMu6HAWw/ZQeeffFZAGbNmcmMWTPyWrulpYXhw0bQFeniyX89kdMx/3rpWWzb5quF81mw6CsmbLQxU7edypP/im2Fzp0/lyVLF7PxRhNYsPArxo4ex2knns60Pffjy5lf9JBVU11LS4tznF5E2G/vA/j1DVezyfhNuj2n3Xfenbse/Ef3+z7/8rNexzqt29TcRENdPVWVVcyeN5sPP/kwp8+bzBtvv8GatWtob2/nxX+/wPYpe3+777IHjzzxMC2tsbTqxmVLWblqJY3LG6mqrOLIQ47irFPP5rMvP817bYO/qEJne/GXE5Uv3njW/VEmcB+wNzBMRBYCl6vqbV6vm0woHKSjtZMVi37UYw8OwI5WsmLRjwqWLSL8+Tc3ceW1l3PTbX/Ctm322m1vfnjehVmPPfm4U7j4ih9zyPEHsvFGm7DpxptRV1uX89oXfPv7fP3MYxkyaAhTJk+hpaUl6zGjNtiQ4844hpaWZq68+JeEw2FOPu4bXHHNZRx+4iEEAgF+ffm1VFSEefbFZ/jnM48TDIYYNnQY5575nR6yJm02CSsQ4IiTD+OYw47pDqMmOGT/Qznu9KO55vLfdD/30x/9nF9cewWHn3Qo0WiEqdvtxC8u+WWP45zWra6q4v5H7uPgrx/IhI0msO3kbXM+Twm22Wobzv/JeTQuW8oRBx/ZIzwJsPvOezB77mxO/ObXAaiurua6X/yO+Qvm85s/XIslFsFgkCsu/kXeaxv8Rnj3pTp2P2ydi/tvBvdIH9HpC56XCbhJoWUCKxavYlV8H66vWZRuEo1GiUQihMNhvlo4nzPOO51/Pfx8j5Ckm1x8xUXsvcc+HDTtYE/kFzOPPvkIn3z+Py676Ip+W7O/ywQMGRCbC69fwP7HrzEGrhgJTMIa/mTBh/tWJlAMJBfgNq06wjeDlkpbexunfecbRCKx4tfLL7rSM+NmMJQroQqb3zw8iy2nmj24oiXgfokAlImBi3QVUoToPbU1tTz6j8f7bb1rrvhN9jcNUI45/FiOOfxYv9Uw9DP1Q7r468tfMnh4cV4DDHG08ES/TAyIXpTZsCwTkzAYyo2qmih/eHomQ0ZETViy2Ol6D7Xdz6QsCwPXVYRdFgwGg5cov35wFqM28qa5sMEDIvNcF1kWBq4Y2wgZDAZvCARtrv/nLDbf1owrKikChZdrpaMs9uDEksJ6gRoMhpIiGLLZ7ZC1bLVjq9+qGPLBGoUENnBfrOsSi5CGYfWeyl/auITvXHguBxwzjf2O2odf/faXdHY5d+BuXN7I935yXlaZZ19wFuua1hWkzx9v+T233XVrQcf2lRdffYFZc9Y3kf39X27kzbff8EWXQnn0yUf4xW+u8FsNQ56IpVzwmwVc/Oev/FbFkC+DbvZEbFkYuMEjGzyTrap896Lz2G+v/Xn+0Zd47pEXaW1r4Yabru/13kgkwsjhI/nDtX/OKvdvv7+tuxN/KfHiqy8wa+76XpMXfPv77Pq13XzUyFAOVNVEufjP89n/+DVYZXFVG2B0vOyJ2LIIUa5bub77+7xn5/LRTR/S2thC9cgapvzftow/eELBsnuNywkU97ic2tpaPvnsfyxfuYIff++i7qJvpzExAH++9U888ew/GTJ4CKNGjmKrzSdz1qnfclz78y8/5+XXXuKdD97h5ttu4o+/+RM33fpn9t5jH6qrqnn4iYf4wzV/AnrOknv9rdf44y2/p7Ozk7FjxvHry67t0XwZ4B/338n9j95LIBBk0wmbcsPVv+fjTz/iqt/9ko6ODirDlVx92bVsPH5jHn3yEV789wu0tbUxf8E8vvmNs+jq6uKfzzxORUUFt9x4G4MaBnHquSczabMtePf9d4hGI1x92TVss9WUHuuuWr3ScbzPO++9zVW/+xUQG7ty9y33UVtTW/D3yNAXlKvvm8OmW7eZbMlSJbLIE7FlYeBa1sbi8fOencs7V79FtD22Ide6tIV3ro6NeCnUyJXauJxlK5Zx760PMGfebL5z4bkcNO3gtGNiwuFKnn/5OZ649ym6Il0cc+qRbLV5bGZturX33WOaY7eUXXfajcuu/hmtba1UV1XzzAvPcMj+h7JqzSpuvv0m/v7nf1BdVc0td/6Vv99ze69mxrfc+Vde/ucrVFSEu0O3G2+0Mffccj+VwTBvv/0Wf7zpRv76m1uwsJg5ewaP3f0EnZ0d7H/0NH50/kU8fs+TXH39r3j86cc44+QzAWhvb+Of9z7Ju++/w6W/uJinHni2x7pX/e5XjuN9br/7Vi77yRXsMGUHWlpbCFd402rIkJ2JU9qYsEU7FeHS6cpkSCE0xhOxZWHgErd1H930YbdxSxBtj/LRTR/2yYvLRjGNy9lvr/2xLItNN96MFatWAunHxLS0tDBtr2mEw2HC4TD77LFvt5x81w4Gg+yxy5688trLHLjvQfz79Vf48fkX8e777zBrzixOOusEALoinWy79Xa9jp+06SR+9PMfMm2v/dlv7/2B2Fy2i6/8CV99NR8RoSvShYUQIsDOO+wS86hqaqmrrWPfuO4TN53ElzPXD7I89MDYTcWO2+9Ec0tzr33PdON9tp+yA9fccDWHH3QEB+xzADUj3c8AM+TG2E3bKaGOgwYnAoXN+sxGWRi4usE1tK5tpbXRuRlxuudzodTG5VRUJLUCi18V0o2JuePev6eVU8jahxxwGPc8eBcN9Q1M3mJramtqUVV2+9puXH/VjRmPveXGW3n3g3d45bWX+cvfb+bJ+57m93+5kd132I1vXfd3FixewHHfPq77/VUVld2/W5ZFKP65LbGIRteXjaSGtFJH6KQb73POGd9mr9334d9vvMpJ3zqBW//4dzYZv0nWc2Bwn/kzKrHEWLjSZognUstiO7ZuSC1iCdUjaxxfT/d8LpTauBwn0o2J2X7KDrzy2st0dHTQ0trCq6+vH0qYbu2ampq0Uw122n4nPvvyUx58/AEOOeBQALbdelve/+g95i+YB0BrWytz58/tcZxt2yxpXMLOU3fhR+dfRFNzE61trTS1NDFqRMxzevCpB3sckzr3Lh3PvPAMANM/nE5dbV2vaQ7pxvt8tXA+kzadxDmnn8vWW27N3HlzclrP4D6z/lfNrE+q6Oo0G3AlS+ez2d9TAGXhwQmCqjLl/7btsQcHEKgMMOX/ti1cdomNy3Ei3ZiYbbbahn33nMYRJx/K0CHDmLjJJOpqazOufcgBh/Hzqy7lrgfu5A/X/qnHOoFAgL1334fHnnqUa6+4DoAhg4fy68t/ww9/+oPu0orvf/sHTNhofcg4akf58WUX0tzchKpy2gmnUV9Xz7dOPZtLrvwJv7/t90zbfVqPtTTH+VLhijBHnXI4kUgsySSVdON97rzvDt6e/hZiWWy28WbdE8UN/vCzb0zg/o8+w/QqL1Fs57mSfaUsxuVEI1FmfzQP1P0syr7Q3+NyCqGltYWa6hra2ts45ZyT+OWlv+pONCkGggQIE+rhsSlKB11EslT3n3ruyVx0wSW95sK5hRmX03+Mm9jGTc/PJFRROtczQxLVZ2LVX1Lw4WU9LscKWASCAaJdUcYfPME3g5ZKKYzLuezqnzFrziw6Ojs4+tCji8q4Ad1GrIIQFoKN0pmDcTOUPoGgEggogaBy6DdWEgwZ41ayBLfK/p5CxHoitcgQESqrw93lAsVCf4/LKYTf/eoGv1XISoRoQQbtrr/e64E2hv5gz8NX88PrF9DeGqBuUIRgyG+NDIUTQCr38URyWRi4jraOojNuBoOhcM766RKqapSqGtNIvfSxIboErNzzD3JlQGRRZttHXLeiKePrBoPbKIraJmTmDcrIsWYMzsBB0eY/eCK55A1cZWUlK1euzGjkbHOhMfQjitIeaWfJrOV+qzJAERoXmpjkgKJrtidiSz5EOWbMGBYuXMjy5ekvJh1tnaxZsbYftTKUM2orS2Yt58Ern/RblQHL+Qdvxm2vfUH9YNtvVQxuYHnTx9VzAycitwOHActUdXL8ueuAw4FOYDZwpqquKUR+KBRiwoTMWZGqyuVHXct/n3yvkCUMBkOR0bQ6yC/PHs81988hUPK36QYiX6DagYi7PV37I0R5B3BQynMvAJNVdRtgBlB4AUQOiAgbbzueQLDkI7IGgwFQFTrbLVqbzd/0gCHa6LpIz78dqvofYFXKc8+raiL96S3Am1bSSfzrtpeJRkw4w2AYKNTURamstomYfJMBQBSsYa5LLYbbn28C3jQiS6K9pcPrJQwGg2f0TBSrrI5y2GkrCVVgQpQDgdA2iFXtulhfDZyI/BSIAPdkeM85IjJdRKZnSiTJRjRiOlsYDKVMsMKmsiZKKGxz8Ckr2eWg2GgjM+R0AFB/lSdifbv3EZEziCWfTNMMOf6qegtwC8R6URa6nvHgio9h2sB4RhImRAddzKORFWKyXQ3ORDrhsNNWcsw5Kxg5xsQlBxQtt8Mg942cLx6ciBwEXAQcoaqetxhpbWrzeglDngzTBjZjNJVUIAiVVLAZoxmmDX6rZihKBBAev3U4Mz/u23xFQxHS/ijqwUQBzw2ciNwH/BeYJCILReQs4E9AHfCCiHwoIn/xUoePXv3US/GGAhjPSAIpX78AFuMZ6ZNGhtJA+NMlY8wE7wFHFI2ucF2q5yFKVT3J4enbvF43Gds22ZPFRhjnThTpnjeUMaKEK20sC+oHR9n/hBWsW2XRMNT8XQ8oAoNcF1kW+UeRDtOQtdjooItKeo8G6iC2t2L25wwAiIIKHW0BANrbLF55bAgnX2DaoA0sqhBxf3uiGMoEPGfpfPPHUGzMo5EoPe/Ao9jMo9HszxnWoz1TJNUWVi4N8dXMSp8UMnhCzbmIB+mwZWHgPnvzC79VMKSwQtYyk0W00xlrTkwnM1nECllr9ucMGbEC0LIu4LcaBheRgDd/22URoly73IzLKUZWyFpW0Dvs2Jf9ORPaHFgEgjbRSM+bHduGiVPMfMeBhGoLXpQzloUHt8UuE/1WwZAHiX24XJ9PYEKbAw9VIVwZa9JgWTbhKpvvXrWIikqTRjmg6PivJ2LLwoNbtWS13yoY8mAejWzG6B5hysT+XCY2YVTa0KaTp2gofuyoUD3UZsKW7Ww0qZ3DT1/JZtuYutYBR+RzT8SWhYF755n3/VZhwONmaHCFrAXFUV66dYZpA0Gc92VM6UFps3pZiJa1AX7z8GzCxnMbmAQ29ERsWRi4aNTUy3hJIjSY8J4SoUGUPhm5VK8r0zrjGYmkieJnC20aih1l5NhOY9wGMqFtPBFbFntwk/fYwm8VBjT9lfWYaZ10XpqiWUObhuLnkpvm+62CwUuiSz0RWxYe3E/+cT7HDfum32oMWPqrK0mmddIVjkexGc9IJukYuogiQJCAybAsISqrbTaZ3O63GgYv6XrPE7FlYeCWzm4kVBGkq9N0NMmVfPbU0hmXCFGm6kTXUvbTrQOwkiY2YHAPD8/GxopnUwJUJH3d3QijGrwnGLLZ68g1AKia0TgDFtubZhxlEaKc+8kCrKApDM2VfNPtnbqS2NgEsFxN2Z9HI0rvfRhBGEpdr8LxKDZWhq+4l8Xjw7SBqTqR3XQrpupEU6qQF4pl2VTVRBm9cQdn/3yxMW4DHkXtVa5LLQsPbuykDckwcs6QQqa9Lqd0e6esRwurh8eUTUYy6bzHFbKWSTrG8ZgwoV6JKbvpVlk/qxcZll4k3ZQbCkw7bhXnXbUYqyxuw8udQPzHXcriq7PlLhMZtuFgv9UoGQrZU1sha5kuM3hDPmW6zCBUYMp+Nu8xnyLwXLInvciwNK3G+oqgtsVz9w+labWJvJQFoSmIZZotF8TyhStZtsD9WUMDlVyNSKYwXKHdSDIVa0PmJs2pOL03l+P6ihkF5A7BoPLBa3V+q2HoD8IHeSK2LAzcPb96mEhn1G81SoZcjEg2TysfQ5QsM1uxdqJJcycRNOm/jRnVy9CmNnTuJEJX/Ljk5s5uU6hxN/REBMJVpoa1LGh71BOxZbEH9+ErZqJ3PmTqJJIg2z5dLjJSyVSs3UXsBmWYNrAxowgR6H5vslGspIJJjKFOq5krS9I2dPaSQluNGVIQ2GEv0yi9LIjO9URsWRi4uiE1fqtQNDglcICzIcpkGHIJw+VrXDKF8CykV/JGOgRhQ4YwXBtYzlqGUkeYEAoIeF4DV4hxN8SaKe84rYmv7beOdauCbLdnk2mqXDZ0oPYaxBrkqtSyMHBHn38o15z6B7/V8B2n7L6JjAboTqfPNeMv20TuQshU5xbAYiJjsHIcqiEIFQTZkCHdnl7iyP7IavTDcyxlgiGbq++bw6RtW6msVlMWUI7Ya8BlA1cWe3B7n7gr1fVVfqvhO05hRSv+XzK5ZPwVsseWiWHaQADLsc5tva75X/HShTxNVmNxsd/XV3UbNzDGrfwQCDiXAPWFsjBwgUCAw87Z3281fCefLL5s7800kTtfEp5liGBag5Tu+b5gshqLh32PWdNt3AxlSGAyIu4HFMsiRAkQtU02VqYQYCq5XGrcCMMN04asoUdF0xo4RbFRLCRvI1jKWY0DbXJ5Z0dZ3Gsb0hGc4I1YT6QmISK3A4cBy1R1cvy5IcADwHhgHnC8qno2lbRpdTNP/PlfXokvGZyy+9IZDyF2EXXrojlBRzGKIQgx47mEVTTRymaMTmvcFKWDLgJYhBy+qjbKDBayQtYyQUf12G9LluH0+aLYrKTJ1V6ZqXhlhAZip5Rn7x7C5J2aqaoxXlxZ0vGSJ2L747bpDiC1iu9i4CVV3Qx4Kf7YM2Z/OI+uDtNo2SmsGMG5PlAQ1/aoEsYn4WVZ8SzHbBmRiT2+IIFee3MxwxXbVxymDQylLq0hS3xeO14z104nS1nNBgx2tVdmMvn288yHgdgp5Y1nG3j50cGYjnrlSrMnUj334FT1PyIyPuXpI4G947/fCbwK/MQrHSprwl6JLjlSw4rDtIFJjHE0DmFCrng4oxw8q5ihy0wAq0eNW7KRS8hLGI50XmAAi7fk8x7PpQuLpvbK7IsHlm8/z3wYiJ1SQhU282dUmuQSg6v4tQc3UlWXxH9fCt7eeo6eOMpL8SXNClnLJjrKMQQIdO/Z5RMGSzYMidozJ7LtmTkZxcR+WzKZsi+d2otlMogJI5EpDAjZa9y8NEJelGj0L8ox5yxnxsfVzPigmtpBUY49dxnHnGPa6ZUv9Z5I9T3JRFVVRNIGJkTkHOAcgHHjxhW0xpwPS3MacH8lEsxmSU57c8keSDrdUg2D2zfk+RjLRNlCqsHNlNCSCNmm88AmMabHeukMv5dGqNQ7pUzZtYXTfrzU7LcZ1hPa2hOxfqUuNYrIKID4/5ele6Oq3qKqU1V16vDhwwtarGm1N/FdL/FyDycVp725dIQJZdTNyTC4Sa6ZkooSxaaO6h66ZqulC2AxTBvSelrikK3ptP/ldp1gMm6WaPjB/sevIlxljJshia73UXU/AuGXB/cEcDpwTfz///Rysc2/tpmX4j3Byz0cJ1L35nbSzXvNc4NYhC5Tx/9i2Qdy6mSSCxYWGzMqr5IK6B169LpdVyl3SrECama8GVLoAG0GcXesmedfMxG5D/gvMElEForIWcQM2/4iMhPYL/7YMz59/QsvxXuCn4kEsa7+zl8NC8na8b+YKKRAPEQg66idVJxCjytkLfNopIMuwoS6Mz7LnZcfHUxbi7FwhmRsFPdHI/VHFuVJaV6a5vXaCT5/Z1Z/LeUaXu7hZNrby6XwOpPRcKPjSCJhxIvuJbmS6oFl0iVd6HEg1qu5wfRX6njtmXr2PnINofg9kcmeNIhEcdsk+Z5k0h+M3nQDrICFHS2dbiZeJRJkyw7MlGHYXwiSsSel1yQSTRJGbiKjHQvIIXbDsZImxjOSSTqmxw1Df4eZSwNl3GYd7LBHM3ZEiKAEi8/xN/Q31gaIuF/OVRZxgn1P2g3LKq1bRK8SCTJddPuaIOKnx+UWNjazWdL9eDwjezWjhpgR/JKFWFhsyBDHhJuBWK/WV7bZdR1X3TOHwSMiVFYroQrjvRmA2p97IrYsDFxNQw2HffsAv9XImxWylukygzfkU6bLDFfCWpkuuv194U1M43aiv4ylokSIdt9EzEi5iUh3ToIE2IzRVDg0iE7cMJjJ3r2xowGGbtBlkkwMPel4yhOxZfE1U1X+++R0v9UoCjJddNO9lskQ9QXbxzBkgkQ4NF0iSPpzQkZvN0zI01KBUiUYsunsMC6bIYWO51Ftc11sWRi4eZ98xerG8tnzGKYNTNWJ7KZbMVUn9rhgO110FSVMiFY6HF+D/DyqXI1hIRMA8iFRC5fNQAcJpK03TGeksmndQVfJ16t5wSfv1JDbrApDeSEQXeq61LIwcF2dEbra0xcvDySyFYgnLrqJsBysL14eTC1raaGTSK/X8iFdQkZ/kyjsVgedUt+XTHLhdjojlSnMqCgraeo+3u0wcykT6bR4/LbhpqmyIYUIWCNcl1oWBi4UDpXNH1QuneZXyFqsuP+UTMLI2diuelZ9kdVX45ite0k6+al7b4lzGibExoxiJU1p6+QEYQMGm5o3B47/7jKO/tYKk1hiSEEQq8Z1qWVRJmAFrJIrEyiUXDP3Ml1f+ivZJBfD53UI00YJpFljqk5kJU2MYnCPTMoKgoxiMEtYzVDqHOvkAlhMZEzZ17wlU9sQ4Rs/bCRcWSZ3m4Y88ObvvCw8uHGbjyYQKIuPmjF0ltiPm6DppytIPKTnRGI3q4tIjzBmX3Ca85bv8X3RI1MWZyUV8Tl2vb87FhZDqWO6zEgr20I86x9aikzatpWuTuO6GRywhnoitiw8uFVL19DVWR4DT1fS5Nh/MXHBThQtZ96TiiVSpBaZJydIJObI9ZUuotjY3V1V8vUe0+335eL55bJeJjmJYzP1rTSF3cqUXVvYcdo6KqujBALGezM4YHkzMa0sDNx91zzmtwr9wjBtYAMGZ7woO3kjqSS6cTi185qgo+IDTPsePlSUOSzpEcLbTbfqk8xE2n/quB+nx32depDwlp26ziRTroXdIsqlf5nPjvs2Ea6yiXZBIAS2jamDM/QkUNgotGyUhYGb8W7p9aIsBDdG1STqtJy61U/QUXl350+HorTQ3t3iqouoa1H4hH5OE8CTH4cI5hTedJqNpyhWfLTO+pZezj08y7Wwe+cD17HjPk1U1cT2vq14JybV2I9JNDGsxxvPviwMXFV9ld8q9At98RRie2vRXh5VMqNcMm4QMzA1VHbLcxrN48YaTsYpn/ek8wQT43h6NE9WSnoQqdvse/Rqqmp7J3YZw2boRcdzqL0GsQa5KrYsDJxGyyPun8sMs0RKhpVy0V7MKubKEsdjEtMH3J/O7f2VLtc12ukkTKjbkwwS6HWenOQl77F5PQOu1OjqFBOONOSIDZGZULGjq1LLwsC1t3T4rUK/4LQXZKM99stiHotNF1GCBBzH5SRfoFfSxAYM9nRKdz5k88gKITE9AGLNlhPnI9f9wGTPuZQHkbrN8w8MYdeD1lJVUx43mIY+Ym3gvkjXJRYhW+460W8V+o3k1lRdRIgQ7WUQLCyi2N3dNYDu1l6TGNOjC8qGDMlo3ApN0Xe7u0lf9Eht1TUxntqf695Zue6xZePD12t59p6hpjTAkB0ZjgTHui62LAzcMRcc6rcKnpNo0ZXobh/r4GERyjB9e6pOZIKO6tHayykZwwk3GjDbeUzMzqZPcvZkvvKcbgAmMYZWOrLW6ZXzHls6ghU21XVRQmHl4//WYpfJFoGhDwx5yBOxZRGiHDF2WKy4awD/naVr0ZUIUaaS7KEVGvLrS6iwv2fH5RvaTLQtczL4iXNa7nts6bjn3U9pXFTBBmO7aBgazX6AocyxsEIbeiK5LAwcQLiqgo7WgdtwOV0GpVPRds/Xs1/0s9WUuU0hEwwS709n0BN7jvlMK0/vLcIb8mleupUTbz7XwMGnrDbZkoYccb8HZYKyCFHGGNh/bZnmvCV3w8+XKDaLWdWjm77XZOv+n4mEQU8mis0clrj2Dchnzy3T6KKByqTt2rCN42bIFfGuEULZeHCDRtTTOG+532oUTGp2Y2pozCmDMrVoe6pOzFpGkIyNdrfnmsv6EoJc5BTq5fV1kkGmLizjdWRenx96f4589twS+6KJf5NKKpjEGOq0Om1JRqmz4YQORk/oIFA2VxZDn5Fqz0SXjQe365FT/VahYLLNeIP0c8tSjWC6ES9OCM6d8HORY7N+2Gh6+b0NWRQ7Lx2TScxhSy1z2JhR7KZbESZUkBdb6MBSp31RQdiQIQPWk6uujRKNDuxoicFl7IXYK0/3RLSv91ki8gPgW8SiUv8DzlTVdi/WKuVauEwz3pJrrrLVYKUWIkPmUKDC+lZUDnLStaZSlKWsLiiBJZgm6zMbCcOVvGa6JJpMLbxS6aAr48SATKTfF5UB24B53heV6MCfSmVwm67/YkdWYgXdnSrgmwcnIqOB7wFTVXUyEABO9GKtaDTK83e86oXorLixB5PrjLdc1kyeMP0lC7u9EyfPJtO4lxWyNqNpGEpdQaHGQiaIJx+XS5mDIDnto/W1BCDTGgO1AXOky+LGH4+hvVWImn04Qz60P+G6SL9DlEGgSkSCQDWw2ItFPv73Z0Qj/X9bmUtoMRcyJZD0Zc1UY2c7GLnUaeC56lXsF/BECNOJQsKRTsyjMW1IdCAXh7/21CAuOm4TmtYU5o0bypTARq6L9M3Aqeoi4LfAV8ASYK2qPu/FWn55b5lCi/ngtOeVzrsodM1MHlk6Y5VOr5U0ZVyrGEgknzjp/yULmS4z+lzftkLWsphVZVYcrgzfsJP/u2oRg0wNnCEPJOx+noRve3AiMhg4EpgArAEeEpFvqOrdKe87BzgHYNy4wmYGzfpwbp90LZR8Q4vpyKeJb1/WTNesOeFtJGdyxlL51w8sTe5rGWvMnDnM6HUtXba1V9KUd3PkbJmsTsyVJTRpa5k0YFYOPGkl37tmEcHiduANRUcY7FVg1bsq1c8kk/2Auaq6HEBEHgV2BXoYOFW9BbgFYOrUqQX1Ihm64WDmfbKgb9oWQDaDkQ+5NvHty5qZSg1SU94TpqmCYLfXk7hoT9LcJn3b2DkNYE0mlwSRXEbkDKWOufHRQLmcV6eU/x6jcjJQLg2YDz1tJedcttgYN0P+SAgC7ncz8XMP7itgZxGpFhEBpgGfe7GQX70o8wktFsOaTqUGS1kdG0rKmLTdUFJDoLkY0+QO/vmQSCRRlE4i3f/viv+eayF6vl60W+HmgYqIcvqPl1JZPYD74Rk8Igi1P0YkvxrVHCX7g6q+LSIPA+8DEeAD4p6a2+x08PYEQ0EiXREvxKfFj/lgfVmzL6NywoTYTbfqbguWbYhobNZa4fdXFhY2Ed6QL3q9lkshejYjnHou3Ao3D1Qqq21q6s2em6EQIlCxmyeSfa2DU9XLgcu9XuedZ9/vd+OWwI/wVCFrOoXg8qljS7wvuY4tk5FLN+UgHzIlv6SGWpNRNKNH63QuMhWIT9BRA7YzSa60t1q0NQeoG2yMnKEA1l4Gw+5wXazfZQL9wq0X3+O3CkVPuq4bfSERSkwmYfTcSDBJ54UlQq0Roo7rr6Y5o0eb7lw4GTlBGMWQArQfWKgK99w4gvZW08XEUACRt1G7xXWxZWHgGueXbg/K/iKfUJui3WNVcyGxp2dnSf6wsbv30rLJzrSvmAgvBuKDXRO62iiLWcVnMj+j7EwdSJyfN4Dyzkt1rFlhmlAaCiEAUfez3cvi2zhi3DBfsihLiXTZl07NhhMF0LvqVlkv7smtrnbTrRzfo2ivvcJ0+2hO701mgo7qEVoNEiCKzYw8irbTnYt0mLQKCFfZnPqjZWwwbuAWsBu8RMFyP2GrLDy4b159irnNzkK67MvUUTnJ3T2WOBQxpx6f7GVl6nySWlhdSBH2MG1w3DfMN9sxUwcSp5DnElblLHugogq7HLDObzUMpYo1AgkMd11sRg9OROqB4ao6O+X5bVT1Y9e18YhdDtuBHQ/Zjnef/sBvVYqWTNmXyaNykpkrSxipgxwbJCcaLuc60icffdIVXGcqME8NO2Yq2l4ha6nT6l7GMorNWloYRC2JAfFLWFX2CSagbLt7M5XVpsuyoUDsRdh2M5ZV66rYtAZORI4HbgSWiUgIOENV342/fAewvauaeMzy+Sv8VqHoKST7chaLHTMWBWEDBtOkrT0MRz4lDE76ZCq4zrSPmOw95lK07VUHkkK6oRQ/wtKv3K9hMpQZ7S9B9ZGuiszkwV0K7KCqS0RkJ+AuEblEVR+jBAN+KxaaMJIXZBqdU8hIn0wM04aM62TaR0z2Et0aP5SrzgmDFiFKAKu7/i+fbijFTk2d8d4MfSTg/ozETAYuoKpLAFT1HRHZB3hKRMZSgvvqoUpTkOsVK2Rt2vZcfSmETjYOiQLyTCHIL1nYy5vUeOZksgHpr6LtVE8x5PDn5mRYS43K6igHnWRuIA19IYB4UOydycA1icgmif23uCe3N/A44JwOV8RU11WyeqnfWpQmuYTV3Oy7mVgz2ThkG4TaQVfOIVC3dU2Hk6foROl1Q1GsQKw9V6hC2W6PJvY/wRg4Qx8ITiW2E+ay2AyvfYeUUKSqNonIQcDxrmviIcu+Ws6iWca6FUKuTYbzSSDJhVyNQ+o6uYQV3dY1HbkarlKaDWcFlJFjOzj6rBV0tFtss0sLk7ZrRUpu08JQVETe80RsWgOnqh+leb4LKKnWIA/97km/VShZ8tmvcrPvZq7GwUazDiZ18kBnssjzZI9c6ulKbTacHYWzf7aE7fdqoqqm5HYqDEWLN60Uy6LQ+9M3vyzBXcPiIJ/9Kjf7buZiHBRlRtKYHifSeaAzWdRdgO4VTp6ijU0Em1DS/LzSSjARfnXOePY9dhXnXLaEBjPU1FDElIWBGztpNLPen4MaI5c3Xu5XZdrby6VhcmryiBO5eqBe4Mc0if5g9MYd7DStyTRWNriIN2UmORk4EakCxqnql55o4THH//gIXrnvNb/VKEm82q/KtreXahxijZPJ2/Pxe8zNQBt2OmGLNm54YhbhKhurLPogGfqFmrM8EZvVwInI4cBviZnYCSKyLfALVT3CE408YJMp4wmEgkQ6/RmZU8p45YXk4lm5YRz6K2OyXDjrZ4uprLZNUonBVaT2fE/k5uLBXQHsBLwKoKofisgET7QxFCVeeCH95Vn1V8ZkubDFDiZj0uA2FrgwHzKd5Gx0qWrq1a2kdrPamtuIRsx+QTGRqfGymyRmw6VrGG3IB6Wt2cQlDW5jox2veiI5Fw/uUxE5GQiIyGbA94A3PdHGIz5540sqq8O0Nbf7rYohTn96VgNtH8wflG13b2LoBibMb/CA6BxgH9fF5nI7dj6xziUdwL3AWuD7rmviIZXVYTraOv1Ww5CE8axKCaUibHPRHxaYxBKDN3gwCw6yj8sJAE+r6j7ATz3RoB/YdIeNsaOmGWyxYTyrYkYBAVEGDY1wwxOzjPdm8I7me6HqMNfFZjRwqhoVEVtEGhz24UoG2+y/GQx5M2HLVg48YTUHn7LKzHozeEt0uidic9mDawb+JyIvAC2JJ1X1e55o5AGfvlGS5XsGg2/UD47wp3/NJFgWrSAM/uNNam4uX99H4z+uIyKDgFuBycRiIt9U1f+6vc7D15telAZDPpz98yVYphzA0F8Ep3ojNtsbVPVOT1aO8XvgX6p6nIhUANVeLLJ66RovxBoMA5aa+iiWN6VJBkNv6i7zRGzWnCgRmSsic1J/+rqwiDQAewK3Aahqp6qu6atcJ3Y6ZHsvxBoMA5bXnh5EW4tJmTT0E61/9kRsLt/gqcCO8Z89gD8Ad7uw9gRgOfB3EflARG4VkRoX5Pbi6z86HDHtFwyGnPnPk4OY8WEV7a3m78bQD3R96onYrAZOVVcm/SxS1RuBQ11YOwhsD9ysqtsRS2C5OPVNInKOiEwXkenLly8vaKGO1g7UjBIwGHImGhEuPXljPnuv2kzhMHhPcBNvxGZ7g4gkx/csYh6dG7lVC4GFqvp2/PHDOBg4Vb0FuAVg6tSpBf2p3fHzBwrV0WAoU5QDTljFdru3mN6TBu+pvcATsbkYqt8l/R4B5gLH93VhVV0qIgtEZFJ8DM804LO+ynViVeMaL8QaDAMMBQEroHzrZ4s59pyVfitkKBu8aaOYi4E7S1V7JJW4OE3gfOCeeAblHOBMl+T2YLejduL9Fz72QrTBMCAIVdjsuG8TQ0Z0svUuzex28Dq/VTKUE633QsUOrouVbHtTIvK+qm6f8tx7quq+NlmYOnWqTp+ef8W7bdscWn2KmQdnMDiw0aQ2rnt4NsGKWM/JYAgTljT0M2Fk5AeIFLb7FbdJvYrp0koTkc2JNVluEJFjkl6qByoL0sInopEow8cOZclsMwPMYOiJcuUd86gbHDWNlA0+EgHtgAINXDoySZsEHAYMAg5Per4JONtVLTzmlfveYMVCs59gMKSy0cQOBg2LGONm8BdrCGK5XyWW1sCp6j+Bf4rILl60z+pP/v3Qf+nqMOFJgyGVRXMreOOZevb7+hq/VTGUMzLYE7G5+IMfiMh5xMKV3aFJVf2mJxp5QO0gT+rHDYaSJ9Jl8fufjKW2wWbnA9YnlqiafThDP2IN8UZsDu+5C9gAOBD4NzCGWJiyZJh2yu5+q2AwFC2d7RZ/v3aDHs8Z42boVyoP9ERsLgZuU1X9OdASb7x8KPA1T7TxiM128KZK3mAYKCydX+G3CoZyJrraE7G5GLiu+P/XiMhkoAEY4Yk2HrFqiTcnz2AYKIzeuMNvFQzljN3n/v2O5LIHd4uIDAZ+DjwB1ALezDbwiHt+9bDfKhgMRUswZHP6RUv9VsNQzkS92fXKZR7crfFf/w1s7IkWHtM4f4XfKhgMRUttfZRtd2v2Ww1DORN53xOxuTRbHglcDWyoqgeLyJbALqp6mycaecAex+3MjOmz/VbDYCgKdt5/DSecv4zNt29j5dIg//7nIEJhMzLA4CPqjQeXyx7cHcBzwIbxxzOA73uijUdM2skkmRgMAHsevppLbv6KLae2YVkwfMMIx357hSn0NviOaqfrMnP5Wg9T1QcBO6aERoCo65p4yFM3P++3CgZDUXDmJUuprO7prZmSAENREJnnushcDFyLiAwFFEBEdgbWuq6Jh0S7bL9VMBiKgpFj3L9LNhhcwYNi71wM3A+JZU9uIiJvAP8gNuamZDjsOwf4rYLBUAQojQvLrN6tch0ybC4yciYybC5UmjFAxUktEhjmutS0Bk5ExgGo6vvAXsCuwLnAVqpaUsPVtt17K79VMBh8RgHh9qs3oL21Z0zSHqgBjsp1SP0yJBhBhNj/65cZI1eMBDb3RGwmD+7xpN8fUNVPVfUTVe1Kd0CxYgXMDrrBAPDaU4P53Q/GsmR+BbYNyxaFmPlxSU2/yhmpXYlYKfuNliK1ZrJI0RH9GLXdL1XJVCaQfJtXkvVvCRZ8schvFQyGouE/Tw7mP0/GureHKmz+/t/PfNbIIwJpJoike94tKtfFjGggAtEg2jwU2uv7X0ZJEYDoXLC2dlVqJgOnaX4vOT55/Qu/VTAYfKZ3qqQVsPnBbxcwbIOSSorOnWgQgg7GLOruUM0eJMKiCc8xGIH6ZbELaK4Gyg0ZJUc7WO53gMwUu5siIutEpAnYJv77OhFpEpGSCmK//2JJbRkaDC6jbDihjVBFlNi9qk2owuZ71y5k32PXDNgyAW0eito9P5zaEvOGPMKNsGh5hlYFCYx0XWqmgacB11fzidbmdr9VMBh8JRAUfvvYLFqbA0zesZVQWN03bMUWVmuvj3k9Tjp5pasbYVG/Qqu+YqN2M2LVuirVQ1+9eNj7+F358OVPiHQO5C+IwZAOYfniEENHRth8O49u9oo1rNZej6au76WuboRF/QitFgEaWYBUbOGqzLJIL9z35N2ZtOMmBIIDxik1GPIiVGEzZIR3N3ilFFZzRdc09XU5h0Uz1Of5EVotCmz3J1r4buBEJCAiH4jIU16tEaoI8duXr+Dsa05BrAG64WAwZKBhSMTbvbZSCqvlqms6I5Spvq69Hl03Ao0EUSX2/3UjenqG2erzcpExIHE3PAnFEaK8APgc8PRf7+V7X+eWn9yN2iWdEGowFIASDNlYXgYwSimslouuGcKY6TxAalfGwqFOYdHk92Y7HrLKGJCENnNdpK8enIiMAQ4Fbs323r7Q1tzGb7/5Z+zoQG3ZYDA4M3FKC7e8+gV/eWmWp+uUUlgtF10zhjH76q2Wkrfbb1hYgUGuS/X79upG4CKgzstFHrzuCdQ4boYyY+SYDvY5eg3trQHvSwEyZSwWGwld65aDFb/pFUUaGmNeVPPQzEYoH2/VKVvT4ICNahsiVa5K9c3AichhwDJVfU9E9s7wvnOAcwDGjRtX0FofvPxJQccZDKXM2lVBbr96FFZA2Xa3Zi6/fR4BL//iSymsFmoDy+5t+OOhSGwLAg4Rn2gQ7aiGwLoex6ZNJHEIc3pGsZVp5InaNuJyGN3PEOVuwBEiMg+4H9hXRO5OfZOq3qKqU1V16vDhwwtaqLIm3CdFDYbSQ2lvDdDVadHRFuDDN2p5+q4B6D2kJoLUNWafHlC5Dqlel9arjRkkdQ5jdlT3OlYVtK2ulzFJu9eWDrsPbvZAaCwd+Z/rIn0zcKp6iaqOUdXxwInAy6r6DS/WGrf5htnfZDAMKHpeLDvaAjxzT5EYOLdG2Dhd1KvXZb3IS+3K7CFbS0HjxktBoxa6bgRS2dzrWBGQSodGwXnvqRVu4EqpTMOZICLuj3LyvUygP/jKNFs2GIhGiqBExkVPw/Ginmp8nC7yORoeCcS6vYgAEl/HSpOo5vR8vhmk6WTnQqknrkgVhKa4LtbvJBMAVPVV4FWv5C+bv8Ir0QZDSVBRabPfcav8ViO3FHnIbT+p0KzFdEkicVSdDSW5eEPJetsWaoMkuRFqC6ik3d8rmFIq03Ci8mDE7Q04ysSDaxheOhutBoPbhKuijN+8jaPOKoIbvWyeRuU6ZPhspKExu5eX68U75X3aPNQxqzrxXNrwZSCSfq6KLb2904ANAhqVHgXb2jTc9ZKKUirTcCSwkSdiy8LA7Xz4VL9VMBh8wQoo5121iBufnEW4qghqZdIZpWhwvYEI9M5udAo1Ol7UUz6i40W+vR5dOzK2r5a0z5ZtXy7xeu81QJtGpA+ZagBt3AxdMSHmhXrRqaTUu590vuGJ2BLxX/vGoWfvx60/6ZWgaTAMaAJBmy22b+XAE1f7rUo32jwUklPnWW+EnAxED1K9v/Z6NNQGSVmNIkkGKFOqfHJJQ+W6WA1cDkiKSxDLoIwZLdLJcPJa3SypSA3nrh3p/JmLuYyg8y1PxJaFgasdVENVXSVtTWZsjmGgo1TW2KgNYzft4Gd/m++3Qj3JVBCezcg4eH8SbnXMatRIMOYx5UBOWZXpjhWgeh04ZVEmsFOsYiZDk68RynUyQrFOe+jGmy5TZWHgWpvajHEzDHhEbCZs0c63fraEoRtEGL95kX7n03kvGZI/0u4neTl/LUdEcE4a6SbJK81kaCBvI5Rr0k7OyT2+4c1uWVkYuCf/8rzfKhgMHqIgsNm2bfzijnkMHl4iqeEpOIYvFbAFbUqzn5Sh40hWEt6S1yR9nkyGpvv31NcaGru9W22th6akyde5GviiLyOo9kRqWRi4JXNyi68bDKWJcP7VCzjk1FVY/Z025ua+Tr79LCvXgfQ2bqpkzx5M9aRyJJdklF4kG9sCDE2P9arXxc5RwsjlWh7QlxuBfqEZtVch1hBXpZZFFuU2e7g7JdZgKDaevnsone39XMjtRXuo9np0xYSeWYdpiHlDDi9o/LWRM9f/pHRMSZfQksiodKKQhu09jG1dhhvtaDAnYyMCUp1lOKoCVmT9Zx8xs/AbgX5EIwtcl1kWBq6i2vSiNAxsVJWLT9iYD9+o6bfJGb63h0rn9SSMbbwLiaPxTXNsd+eSNPRlKkO6/pcJQ+NkrLKSWh4QjV3SxUr67Fbv7E8gVrtXFPtvcSz3dSkW/9RTHvv9036rYDB4Rt3gCD//23xGT+gqLIRWKH7v66QJz2VsopxIqsjSzcTx+ALOqwhQ35hbIkdqiDbbmg7lAbGM0BwzEvMMz3qNBMa4LrMsPLg1y9b6rYLB4BlNqwOcd8BEmtZY/WfcIHPRdj8Q83hSnst2zY4bjoK8pUIRMoZtRVjv9SaHaFvrexeVazzRxCk83NCY381F0ey/ATIMkZDrYsvCwO140HaIVQSNZg0GTxD2PnINtQ39O7G+39pDZZo+UOifdUpoz0sSBkw7KtOv5WSYmkZ2G7nun3gWZbquKbne4BRdGy+XB50mKCIT7h0n/OQonvzL83S2dfqtisHgCfseu6Z/vTfonyne6erGKtci4XbnIu8sYVoZNrdbT22PeUM0NPY6JpucrL0rkwlEkEzvTedNNY1Ek8sCkuTlQ6LcAkuLr4sJgGYolO8DZWHg6ofWMniDQTTO9XCarsHgIx3tPgVjPJ7inbZuzMG4JZPOOIngXEBtW2jyuBpb0PY6qFrnnKCRw1q9SGOUurMZ8ym5yGEPsdtbVECt2DieYjRuAMHNPBFbFiHK1x5+izWNa/xWw2DwgNj06afvGkqky29dPCBDtmM6cgnVdWd7xvtQJho8dzdUbhoRqzXT9JfI5CzNbOSkUx4lF9pRnVNoVdeOBGT95yvWSd+28eAK5tUH3qCj1YQnDQOHwcO7+O7VC9n5gHUg8O5LtXz4Zi077Nl74nRJYwsEPNokC0SQuuWOYU7qlsc8074MIc0DqetddJ6plZbTZPFkEvt1xd+iK070M1QVcfnLWxYGrqLK1MEZBg6BoHLDkzMZPqqLYDzxbKf9mrGsfiwR6DeE9EPYXCDbhO4CygkK0yPNZwzECrZ7hRbT6J3s1UllU0a5xUWF68YNyiREeci3phEKu5+CajD4wTa7NNMwONpt3AACgYFo3MjpQp4LTu/Per7qGkGi/ZNlma52L12heg7yJKD5J7X4RWCcJ2LLwsBtu89kdj1qR7/VMBhcYd3qAOGq/i0JKEZ6pM/bmdtspSOjYale18NIFCLfTXp0iUlTw5fLTU7RlQgAVHgzlLosDBzA6Vcc77cKBoMrzP6kirP3msQbz9b5rYqvaONm63+WbQbRYNpWWIXgtDfnu5ecKFRvGpF/kXv8PcU36bsSqTrWE8llY+AGjWjwWwWDwSWERXMr+duVo4lG/dbFY/LplpJpX6nYQnKFYsv6cgLp6cXmRDRYZMYNCG2HVEzxRPQA+VfPziv3v+G3CgaDq6xaHsSOCAGvsgyLAMcZcelCbBlGx2jzUMdi7pJD6DXmJ9c6PNVYeUFeuDkOKR2RGajaSLaCwwLwzYMTkbEi8oqIfCYin4rIBV6t1dXZxU0X3O6VeIPBFzraArz6xCC/1fCW1G75kWDaEFvG1mHt9ejakd17dSWLqGOLrlTSJdVIuDX3tdL0u3RsmdYXdDW0P+uOrBT8DFFGgAtVdUtgZ+A8EdnSi4XeevI9ohGzKW8YeNz7+xF+q+A9uc6Iy2YM2+vBdt6nKxvyKA/I1O/S3YJxG2170AU5vfEtRKmqS4Al8d+bRORzYDTwmdtrNX613G2RBoPvhKuinPw9n9rP9UfoqhCSW4cldGxoXK9j0dV/pSc19KixpjWFN5iG/PYis5wrVwvG1Zs2PEWxByci44HtgLe9kD94xCAvxBoMPqFU1dp8+8pF7H/86v5fPl0DZHDPyPXVgFauQ+ob1/eRDEagvhFsCwL9F83p63y+5ONF4o3Z7J4DTJ3WcGo67bh3me485+qZuXLDUIVUHeWCnN74buBEpBZ4BPi+qvY6qyJyDnAOwLhxhRUDTth6LGIJapdy8N1Q7ogoFZU22+zSzJV3zCPg01+v5+2fXDCgsdZXKc9ZoLaN2tJL/2RcHRqroNFgL0OQq3xHw2VbaMRab5Q6qsFhWnjCyAHONwkZznNscGoOCtou7HIFN4Gqo/sux0m0J1JzRGIT7h4B7lHVR53eo6q3ALcATJ06tSALNX7yOEQE9bLlj8HgMQ1DI1xy83y23a3FX0U8nuTtigFNZ8Dinb/SGTHXE1CE2L5h8lMjZ+Z2aDoDY9lo4yY9n6tO73Fpo3On/kznOfd/SxdOWP0vPRl2Cv5mUQpwG/C5ql7v5VorFq3CjpokE0Nps2ZFkDEbd2D7/VX2epK31wbUcm5hlSic9jQJxY2kDKfzXMi/SabznOu/ZQZPOGfWXoR6lNrqZxblbsCpwL4i8mH85xAvFnr76fe8EGsw9DPCRV/flCXzK3xNdfd8kreHBjSd8dJ48kau9WSZHvdaM5FWP2hBLM2+L3tyac5zQf8m6c6nbYHYuX3H3Lipic5EWx/quxwHfDNwqvq6qoqqbqOq28Z/nvFirUdueMoLsQZDv9PeanHhURv7W8uVR21aITherJVYZ30366+SyLcNV3cHkUgwbZQuEQbtTqvPMKQ1l24kGdtspf6bRAVUYgY1zTlLe57F7p4f172uTe/WYG7e1DTf7I6cFHxPMukPls4zZQKG0maTya1c/OevGDm2E0vA8rvJnpeTvNvrYzYjaS+o2zB4kbGZJ8kXfqxI2rR9xzlzDqgmBpOSvdtKps/cXo+G2pDqdWApIlmSdFLPczQIEkVSOuPEElsEEDRh5WyJDYV1699Al7ojJ4WyMHBezBkyGPqL2oYurntkNlU1tv+Grb+IG1AZNhdJab+VS8KJttb3yixMeEhuXQ4kds13R1btylg2ZLY6t8p1vTMhEwYqQ5g17TlLuVFJmwCTbDABLZFLaln8uYzcaLjfKhgMBbHLQWu59T9fUl1bRsYtmUITTppGoq31PUfqtMbbdUWtomrX1R3CrF7Xq7Sh1/sS43KgdyutbANvc0nSSbOn1ssbTR7d4woBF2Wtpyw8uCPOO5Cbv3+H32oYDDmz5Q7NnH/NQsZv0VGehi1BugbKSmxvKVMheNNItGlk70Pb67Om6rtaC5cjOa0XiHR7cU5p/hlJNV4ORd6Oza3TnQs3u8LIMPdkJVEWBm7tcvc3pQ0GrzjmnGWc9uNGqmr8rgfwn9gFt7G3ZxP3egDnPSYXWolp1ALL7ndDlwkR1n/WPAxMr4SQdEXe60bEEll67MvZzt1f3BxBpCtQ7XK9Hq4s7g3fMmUChhKhsjrK6RctNcYtQXs9aO/LVMaQmVMX/JTGwInwZUbUcsW4uR0OFUuRhsacuoh0Z0CmZF+mK/KW2pW9mltr03Bvy0IAsKHrQxflxSgLA9cwtAiawBoMOTB20w6ikSJyGYoBK0djn8i4zHTxTtA0Eu2ozGx8snhIuRquHi2z8iBT6YAIsY4m2WryFHTdyN7eaz57mx6XhcTwZnJvWRi4I/7vQL9VMBhyYtmiEIFQEWVAFAO5hsISXkauF+81Y/u8bl+9s3QTudWWWOlABh26yxXiSTOOTZctnJNB8i2mz3VkUcEohLZzWWaZGLiFMxb7rYLBkBNrV4aY/b8qoqUz1cVznAqSnYm/J5+Ld4b3Zl03GszuBSZwMmLxx93ZnQ4eUuK5dIgAaqXtNwk4GnzPu9HkTQ0i7qeElIWB++DlT/xWwWDImSu/OYEvPqjGLvXp026RGiJLd07iocx0HTq0o7rXIWkv9B3VMc9HNL2H1TwU1ozNup8Xm+PWez+vu8tJ9ToItfX2kOoaM3Y/6SaRWZmOaDC2Lzls7vq2YdAPYcc8qDrCE7FlkUXZ0dLhtwoGQ1YqwlGsALQ0Wfz05I3ZaPNWrrxjHoOGmoST5IJkGT4LAg4WJeGNxTt6JBd6iwBVTWhXVdZuHtpRjVQ1OabgqwK2hbbX9himmhEl4z6iCFC9rpdu4jACx5FoMO14m4TxSpsxmTLpwBekHqn7sSeiy8LAiWU27Q3FzyU3f8WgYRE+fbeGoSO72PWgtVRWGxcuFW0a0btWKyW8JuFWx0xL6pb1NGbx8oEe3TyGzU1bX5YYOtrDAAYjmcOIVg4NmZPT//PwohKfWxoa08t18AJdnd/XV4JTcGXsjpNoT6QWGYFgWURiDSXO+EltbDihiy2ntvqtij/kWrvm0KsSiWVJdhuIdIkmliKSoX4OsteXOdTGOU3QTiXb6/kYnZgnmdQPsnalc0E8GdZ0s1C7L3S9ga75HjLk766LLosrfyBUFnbcUMJYlnLR1zdl+eIy/a7mULvWg/b6mAFU6dmxP3FMPi2nUjvuF1jAnKkcQISYQcrmqCQZnUx7e7HkkkC3Yc49EScJNwu1+4QNndPR6CLXJZeFgdtwkw38VsFgyIhtCysbQ/zuh+P8VsUXcqpdy/WYhkYIZA4b9jgmYRwbGqGuMaOxUFsyRtMy7plZmt2oJL+e1E/TkWQPLJGIk2F0T4/HvmZMOiAhiC5zXWxZGLgd9t/GbxUMhqzYUeHjN2tpby3DPeNCmiqnea3bo8vzNHZnNEKvrM3un85w4RMEspQeOLbTCmcIV6cay/Z4uYFTVmhrffFkTDqhHRDc1HWxxeKjeoopEzCUEqqCV5vuRUu6psqZPJ50x/QBEWL7YCsmxJJJUpJZnBI2ckGV9Qkt8QzPHtgW2hSbetLdRJr0Rjpd2YNjVmhi3ab89e4/FKTWdallYeD+/eCbfqtgMGRFRJk4pbUs+1A6drHPEkbT5qHZB4Q6HZdtUkCmll+FGrfW+pjHVLnOMf1fxYZQW9ryhFTSlj2At8NoPSMCkRkQmuSq1LIwcK3r2vxWwWDIghKsUC68YUHskQ/jWnwlg+fRi+Rsyzzp3kNzqqNLkPAaXcgyTGQ7SvU6CLeCOE8nEIteA1qzUVSp/i6g0a8QY+DyJxQO0tXR5bcaBkMGhJ33X8PYTWNNCcrKuCXIxfNIHfOSI90F2vEwYKq32P2+ZK/RthzHxKTefGSbFC4JY5qlXq4giiXV3w0s9wdTl0WSydjNR/utgsGQkcrqKLsfYuYWZiOXIZ9ORiTRs5H2esfWX87JF+mLvXs9tsUxecTxvW5SNKn+LhDYyHWRA+jspGfKXlvy5Tuz/FbDYHAkXBVl/Obt7H7oGr9VKX7SeCzdRi0azC0jM523GO/ZmLdnZGms+388dJpvwbfaoG31kLIHlylUXXSp/n0ijNgrIDDYVam+enAicpCIfCkis0TkYq/WWfCl+wWEBoM7KOMntXPtQ7MJujvMeGCSqft/olFxvqNgEqQWm+fjbUWD3SNlsmJLjxE3GpXYzLamkb0aIKdDtfcQ05InMMp1kb55cCISAP4M7A8sBN4VkSdU9TO311oyx/0CQoPBHYTZn1bx4eu17Lx/UedxFwW5ZFsWkpEJOYY/bUAc9uACEWTY3NgaGcoXNLm9lhMOfTHTlk8MJOMW3Bqx3C8T8NOD2wmYpapzVLUTuB840ouFxk4ye3CG4iXSZfHyI+6GZgYsuUyXLnQCdYbw53o5I3vObtOkwvJENxSHLirrPbX8vK50o38SBjXjmJxSIvIRqu4nzPi5BzcaWJD0eCHwNS8WOv0XJ/DaI295IdpgcIWyzJoslFyyLeNlB4lygh6NmNORodg8NfSo7fWxqQMp73fKrsxY8pANh8bS3WukaxZdknSh0cVI0N1WdUWfRSki54jIdBGZvnz58oJkbLTFGMI1YZc1MxjcQjnirBV+K1EapA7uTOfB5Nu8mQKmXGdJRBFhvXHsiwFK7O1Fg87NojP06ywptNN1kX4auEXA2KTHY+LP9UBVb1HVqao6dfjwwuokls5bRle7qYMzFCeBIGy5Q5mOyMmHPIxWIc2b8w5t5pKi72adWiH9OksJe43rIv0MUb4LbCYiE4gZthOBk71Y6PE/PUvZ9fYzlAQiypRdm0yIMgfSGS3Hbh6FGoM82lw5JbP0ws06tUL6dZYQEhhAhd4a21H8LvAc8DnwoKp+6sVa8z9diB01Bs5QLMS+ixWVNjV1Uc67ypSx5EQ+RqvQUoF8cCgYT8btOrW8Q6ilRmCM6yJ9Nf2q+gzwjNfrbLXrRKY/96HXyxgMOaDUD4lQWW1z8CmrOOQbKxk0NOq3UqVBHh5MoaUCOU8VT5Ds8VWuy62XZqHk06+zFInMgNAWroocGL5tFrY/YAp3Xv6g32oYypxwlc3lt89l+z2bUQWr6FO8iou8jFYhxiC1z2W+WYr90cW/JCcF5IKA1LgutSwM3LrlA6RWxFDSfO+aBWy9c0tBwzgN5G+08jQGee3xGVxGXC8RgDIxcKGw6YFk8JfK6ih7Hr6WirDZC+4TXnowAz1LsZjxoNEylEAdnBtsvM1GWIGy+KiGIqWm3nZ/VIrBXfojMcXgjDXSG7GeSC0yBo8cxL4n7e63GoYyZlVjkNbmgN9qGDIw4LMUi5muD1B7retiy8LAAZx51Yl+q2AoY1SFP126Ie2tgt17hqahGCi0h2VfybU7y4BGIep+U/yy8b1fuud1v1UwlDmvPz2Y5rVBfvqX+dQPMaUBRUl/Zymmy9wMtSHh1oFZDuBIJ2qNwu3cq7Lx4P7195f9VsFQpliWUl0XYdvdmzj9oqXGuBm6SdtSrHpdXn00BwTRGa6LLAsPbvnClSyeudRvNQxliVI3OMLd0z83GZSG3qTJ0HRqqjzgyxWiq1wXWRYe3BuPv+O3CoayRRizaQd2FGzjuBlSySdDc6CXK4R3dV1kWRg4tRUxlbWGfmW9t/bp27VcePSmvPlcPV0dPqpkKDrSDjR1YiCXK4R2w7KqXRdbFgZuk23Ho6YIydCPVISVisr16ZILZ4d558V6QmYsoSEZp8zN1vryK1eo2MkTsQP4lmA9rz7wht8qGMqMypoox5+3jFcfH0zdoAin/XgpW+zQ5rdahmLEIXNTu6oGblNlByTgTaF3WRi4WR/M81sFQ1mhbP21Fr7+nRV8/TtmUrehAAZsU+U0VB7oidiyCFFutv0EXC+wMBgcUSxLOeNik7VrMOSMvcQTsWVh4I487yAz0NvQTwhWAEaO6fRbEYOhZNA2b8aCloWB++KdWabZsqHfsCxYPDdkmisbDLnS9b4nYsviqv/fJ97FjpoGgIb+Y+ioCI0LQ6b2zWDIhc6PUHX/Gl0WBm7QiAa/VTCUCaFwlMNOW0H9YJvBwyM0Li6LPC6DoY90gLa6LrUsDNyh5+7vtwqGskA55NQVfOuy2IZ5uFIZNnKAd58wGNxA6kDcL/Qui9vLoaMGx7IozZ6IwSO22aWZA05YyeDhEayk28agGSZvMGRBoPYCRNz3t8rCwC2e3UhFOERne5ffqhgGIGf9dDGHn7GScJWNZYFtr2+WazrEGQxZCGyJVXOyJ6J9CVGKyHUi8oWIfCwij4nIIC/XG73ZBnR2GONmcJ8NJ3Rw5DdXUFVjd3tuVlkE/g0Gl7AbPRPt15/iC8BkVd0GmAFc4uVig4Y3EKooC2fV0M9M3bvJRL4Nhr6gKzzrFeyLgVPV51U1sfv+FjDG6zW33GWS10sYypC2FsGOmjikwVAw1gaeTXsphmDKN4Fn070oIueIyHQRmb58+fKCFznr16dQUWl2/A3u8sazDQQCxoczGAqm9kLPRHtm4ETkRRH5xOHnyKT3/BSIAPekk6Oqt6jqVFWdOnz48IL12eJrm/GbFy+nqq6yYBkAoXCIQCjQJxn9jRWwCFaUls6lQmtTkGvPH0c0kmGOV6lhbUXuzVv79vdk8JMiiDwEp2BVH5n9fQXimYFT1f1UdbLDzz8BROQM4DDgFO2nYW1b7TqJx1beweHnHkAgmP6CX1Xr/EcbCFqc8csT2f8be1JZ03OwVyAU4JtXn8QJFx3JmImjqBtaR1VtJYM3aECs1PnzmfUMV1fk9HlyIRQOsdGWY/jDm1ez2Q4bEwgGCIYCbL3H5mk/Z0KH2kHZ61JqB9cQCrvrGYslbDBhBOGqCiprwlTXVWIFcv9jtIIW9UPrCl4/GA4SCGb+07ACFoFgACto8cYzgzh2iyn896Wt0Kx/UiFif3YhoJ7CLzK1xIxLmuNlKIQPAiri78s2iK4CgtvA8OlYIx5DRn6KDHkIGfooVJ3h/P6a7yP1l0No+7j8fD5LQvfEMWmOrTwGBt2cg7z+6LwfBgq5UbTIfiOQi/5O3y0B2SA/XUK7Q+0lUP9rCO0MVGV4f/K/jwAVEJhM341jFVjDkcF/7KOczIgfg0BF5CDgemAvVc057jh16lSdPn26a3qoKquWruavF/6DWR/OY8LW4zj3t6cRDAW447IH+M9D/6WzvZP6oXUceOa+nPCTI6murcK2bV6+93We/MvzdLZ1ss9Ju3HE/x1EZbXzRWTdqib+95/PaVrdzOTdt6B1XStzPv6KURNGsHDmYp7+64u0t7az3bRt+MbPjmXwyEFEo1Eeuu4Jnv/Hq7Q2tROuqqCtqZ2ahio627tYtXQ1wYogk3bclCO/exCzPpjH4BENTD1wCvde9Qjvv/g/APY9eQ9O+dmx1NTHjFVbcxuBYICKygqikSiP/+lZXr73NUQsdj5sBzrbO1m3spmdDt6Orx22PU2rmpn+3Ec0frWM//3nC5rWNLPb4TtyxPkHEQwGCYWDTP/Xh7x83+s0zl9BVU2YytpKgqEAnR1dtDW1M2LcMHY5fCpvP/0eH//nc8LVFQwZNZi5H8+neW0rg4bXs9kOExg6aggjNxrOgWfuQ/2Q9Qaqo62Df93+Mv9+8L9U1oQ5+vuHsv20rXnrqfd45m8vsmjmEsJVYSprw4wYO5z9T9uLHQ/alq6OLtqa25j7yVc8csPTzHxvDiM2Gs53rj+duiG1PHHTc7SsbWXKXlux/f5b07q2lf+99gVVdVXsfNgOzP5wLnde8SCrG9eyx7E7c+z3D6WiMsSKRatYNHMpw8cOZcLkcT2/U5G5aMut0PUBUAdUQ3QeBIdBzQVYlbv3fH/7s2jznyC6FIJbIHU/RiqmoNFloC2oDIWWf0DHM2CvABkM1ScjNd8Aew3a9gB0fQnBrUEqwP4KKnZGwvsgEkAjs6DzA7CGodYgaP4DdH0O1gYQ2hoin4IEkaoToOooRJwv4Br5Cm3+A3S+BdYwpOYcpOqQpL+lDmh7Gu18GwJjoOoYJPIJ2voYaBdUbBurc7JGQOXBiFWLRlegbY9AdD5SMRUN7wdtD0Pb82DVQ825WOEdALA73oG1F4O9mJiRCQOtQBXUnInUngd2I9r6UOzzSjWEp0HFVtB6B3RMBwQCo6DyEKg6GOn4N9p6P0QWgFULMgSic0ADULE1BMeDNRQ6PwJdCeFpSPUJaMdr0PYYaAeENofwgRCajHS+jra+APZcsJdBYAjUfBcJTQJ7Lapd0Pw76PoQqImdi8AoqDkFqdgVjcyApt9BZB5YgyEwHsSG4CSoPAqxF6BtT0F0GQQ3hfCeSMXmiFShkQVoy53Q8RLYLRAYBlIPugq0Jfa9Ce+FVJ+IBMf2uAbS9S50vo1KHYR2Q3QxqgJEYuXDUgPtz8UeVx4aG07a+QbadD1E58bOtYwEa1BMN3tmXHp97N9KOiA4GWq/C5EZEJkJwa2QqsMRl6Z4i8h7qjq11/M+GbhZxL6hK+NPvaWq3852nNsGzmAwGAylTzoD50vuvKpu6se6BoPBYCgfiiGL0mAwGAwG1zEGzmAwGAwDEmPgDAaDwTAgMQbOYDAYDAMSY+AMBoPBMCAxBs5gMBgMAxJj4AwGg8EwIPGl0LtQRGQ5MN9vPbIwDFjhtxJ9oJT1N7r7Rynrb3T3D7f030hVezUrLikDVwqIyHSnivpSoZT1N7r7Rynrb3T3D6/1NyFKg8FgMAxIjIEzGAwGw4DEGDj3ucVvBfpIKetvdPePUtbf6O4fnupv9uAMBoPBMCAxHpzBYDAYBiTGwHmEiJwvIl+IyKci8hu/9SkEEblQRFREhvmtS66IyHXx8/6xiDwmIoP81ikbInKQiHwpIrNE5GK/9ckVERkrIq+IyGfx7/kFfuuULyISEJEPROQpv3XJFxEZJCIPx7/vn4vILn7rlCsi8oP4d+YTEblPRLKNPC8IY+A8QET2AY4EpqjqVsBvfVYpb0RkLHAA8JXfuuTJC8BkVd0GmAFc4rM+GZHYCO0/AwcDWwIniciW/mqVMxHgQlXdEtgZOK+EdE9wAfC530oUyO+Bf6nq5sAUSuRziMho4HvAVFWdTGxE+4lerGUMnDd8B7hGVTsAVHWZz/oUwg3ARUBJbdKq6vOqGok/fAsY46c+ObATMEtV56hqJ3A/sZujokdVl6jq+/Hfm4hdYEf7q1XuiMgY4FDgVr91yRcRaQD2BG4DUNVOVV3jq1L5EQSqRCQIVAOLvVjEGDhvmAjsISJvi8i/RWRHvxXKBxE5Elikqh/5rUsf+SbwrN9KZGE0sCDp8UJKyEgkEJHxwHbA2z6rkg83EruJs33WoxAmAMuBv8dDrLeKSI3fSuWCqi4iFtX6ClgCrFXV571YK+iF0HJARF4ENnB46afEzusQYmGbHYEHRWRjLaKU1Sz6X0osPFmUZNJdVf8Zf89PiYXQ7ulP3coREakFHgG+r6rr/NYnF0TkMGCZqr4nInv7rE4hBIHtgfNV9W0R+T1wMfBzf9XKjogMJhalmACsAR4SkW+o6t1ur2UMXIGo6n7pXhOR7wCPxg3aOyJiE+u5try/9MtGOv1FZGtiX7yPRARiIb73RWQnVV3ajyqmJdO5BxCRM4DDgGnFdFORhkXA2KTHY+LPlQQiEiJm3O5R1Uf91icPdgOOEJFDgEqgXkTuVtVv+KxXriwEFqpqwmN+mJiBKwX2A+aq6nIAEXkU2BVw3cCZEKU3PA7sAyAiE4EKSqQhqqr+T1VHqOp4VR1P7A9p+2IxbtkQkYOIhZ2OUNVWv/XJgXeBzURkgohUENtsf8JnnXJCYndAtwGfq+r1fuuTD6p6iaqOiX/HTwReLiHjRvzvcYGITIo/NQ34zEeV8uErYGcRqY5/h6bhUYKM8eC84XbgdhH5BOgETi8BT2Kg8CcgDLwQ90DfUtVv+6tSelQ1IiLfBZ4jlk12u6p+6rNaubIbcCrwPxH5MP7cpar6jH8qlRXnA/fEb4zmAGf6rE9OxEOqDwPvE9tG+ACPOpqYTiYGg8FgGJCYEKXBYDAYBiTGwBkMBoNhQGIMnMFgMBgGJMbAGQwGg2FAYgycwWAwGAYkxsAZDB4iIlER+TDpZ3wBMo7ysomxiPxLRNaUYkd9gyETpg7OYPCWNlXdto8yjgKeIo9CXhEJJjWdzsZ1xBrenpu/agZD8WI8OIOhnxGRHeJNuN8TkedEZFT8+bNF5F0R+UhEHol3etgVOAK4Lu4BbiIir4rI1Pgxw0RkXvz3M0TkCRF5GXhJRGpE5HYReSfekNdxSoGqvgQ09cuHNxj6EWPgDAZvqUoKTz4W7934R+A4Vd2BWNebq+LvfVRVd1TVxGyvs1T1TWKtu36sqtuq6uws620fl70XscbZL6vqTsRax11XKh3nDQY3MCFKg8FbeoQoRWQyMJn1rcQCxEaGAEwWkV8Bg4BaYu278uUFVV0V//0AYg2FfxR/XAmMo0QGYxoMfcUYOIOhfxHgU1XdxeG1O4CjVPWj+ESEvdPIiLA++lKZ8lpLylrHquqXBWtrMJQwJkRpMPQvXwLDRWQXiI2bEZGt4q/VAUviYcxTko5pir+WYB6wQ/z34zKs9RxwfrxjOyKyXd/VNxhKB2PgDIZ+RFU7iRmla0XkI+BDYrOwIDas8m3gDeCLpMPuB34cTxTZhNg05O+IyAfE5gym45dACPhYRD6NP+6FiLwGPARME5GFInJgoZ/PYCgmzDQBg8FgMAxIjAdnMBgMhgGJMXAGg8FgGJAYA2cwGAyGAYkxcAaDwWAYkBgDZzAYDIYBiTFwBoPBYBiQGANnMBgMhgGJMXAGg8FgGJD8Pxrm9gbLU2NLAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 504x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "import numpy as np\n",
    "\n",
    "# Create dataset for visualization\n",
    "size=40000\n",
    "feat_1 = np.random.uniform(low=-7, high=8, size=size)\n",
    "feat_2 = np.random.uniform(low=-2, high=14, size=size)\n",
    "features_vis = np.column_stack((feat_1, feat_2))\n",
    "\n",
    "labels_vis = svm.predict(features_vis)\n",
    "\n",
    "# Plot the decision boundary\n",
    "plt.figure(figsize=(7, 5))\n",
    "plt.scatter(features_vis[:, 0], features_vis[:, 1], c = labels_vis)\n",
    "# Plot original dataset\n",
    "positive_samples = [idx for idx in range(len(transformed_data_targets)) if transformed_data_targets[idx] == +1]\n",
    "negative_samples = [idx for idx in range(len(transformed_data_targets)) if transformed_data_targets[idx] == -1]\n",
    "plt.scatter(data_features[positive_samples, 0],\n",
    "            data_features[positive_samples, 1],\n",
    "            c=\"yellow\", label=\"Original positive samples\")\n",
    "plt.scatter(data_features[negative_samples, 0],\n",
    "            data_features[negative_samples, 1],\n",
    "            c=\"purple\", label=\"Original negative samples\")\n",
    "plt.title(\"SVM decision boundary\")\n",
    "plt.ylabel(\"Feature 2\")\n",
    "plt.xlabel(\"Feature 1\")\n",
    "plt.legend(loc=2)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 6. Dual approach <a class=\"anchor\" id=\"dual-approach\"></a>\n",
    "\n",
    "In the previous sections we took a detailed look at the primal SVM. To solve a primal SVM, we need to find the best values for the weights and bias. Recall that our input examples $\\mathbf{x} \\in \\mathbb{R}^D$ have $D$ features. Consequently, our weights $\\mathbf{w}$ have $D$ features, too. This can become problematic if the number of features $D$ is large.\n",
    "\n",
    "That's where the second way of formalizing SVMs (called *dual approach*) comes in handy. The optimization problem of the dual approach is independent of the number of features. Instead, the number of parameters increases with the number of examples in the training set. \n",
    "\n",
    "The dual approach uses the method of Lagrange multipliers. Lagrange multipliers allow us to find the minimum or maximum of a function if there are one or more constraints on the input values we are allowed to used.\n",
    "\n",
    "If you never heard of Lagrange multipliers I can recommend the [blog posts](https://www.khanacademy.org/math/multivariable-calculus/applications-of-multivariable-derivatives/constrained-optimization/a/lagrange-multipliers-single-constraint) and [video tutorials](https://www.youtube.com/watch?v=yuqB-d5MjZA&list=PLCg2-CTYVrQvNGLbd-FN70UxWZSeKP4wV&index=1) on the topic from Khan Academy."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 6.1 Recap Lagrange multipliers <a class=\"anchor\" id=\"recap-lagrange-multipliers\"></a>\n",
    "\n",
    "Lagrange multipliers allow us to solve constrained optimization problems. Let's say we want to maximize the function $f(x, y) = 2x + y$ under the constraint that our values of $x$ and $y$ satisfy the following equation: $g(x, y) := x^2 + y^2 = 1$. This constraint equation describes a circle of radius 1. \n",
    "\n",
    "The key insight behind the solution to this problem is that we need to find those values for $x$ and $y$ where the gradients of $f$ and $g$ are aligned. This can be expressed using a Lagrange multiplier (typically $\\lambda$):\n",
    "\n",
    "We want to find those values $x_m, y_m$ where $\\nabla f(x_m, y_m) = \\lambda \\nabla g(x_m, y_m)$.\n",
    "\n",
    "In our example the gradient vectors look as follows:\n",
    "$$ \\nabla f(x, y)=\\left[\\begin{array}{c}\n",
    "\\frac{\\partial}{\\partial x}(2 x+y) \\\\\n",
    "\\frac{\\partial}{\\partial y}(2 x+y)\n",
    "\\end{array}\\right]=\\left[\\begin{array}{c}\n",
    "2 \\\\ 1\n",
    "\\end{array}\\right]$$\n",
    " \n",
    "$$ \\nabla g(x, y)=\\left[\\begin{array}{c}\n",
    "\\frac{\\partial}{\\partial x}\\left(x^{2}+y^{2}-1\\right) \\\\\n",
    "\\frac{\\partial}{\\partial y}\\left(x^{2}+y^{2}-1\\right)\n",
    "\\end{array}\\right]=\\left[\\begin{array}{c}\n",
    "2 x \\\\ 2 y\n",
    "\\end{array}\\right]$$\n",
    "\n",
    "Therefore, the tangency condition results in \n",
    "$$ \\left[\\begin{array}{l}\n",
    "2 \\\\ 1\n",
    "\\end{array}\\right]=\\lambda \\left[\\begin{array}{l}\n",
    "2 x_{m} \\\\ 2 y_{m}\n",
    "\\end{array}\\right] $$\n",
    "\n",
    "We can rewrite the vector form into individual equations that can be solved by hand:\n",
    "- $2 = \\lambda 2 x_m $    \n",
    "- $1 = \\lambda 2 y_m $    \n",
    "- $x_m^2 + y_m^2 = 1 $\n",
    "\n",
    "Solving the equations yields \n",
    "$$ \\begin{aligned}\n",
    "\\left(x_{0}, y_{0}\\right) &=\\left(\\frac{1}{\\lambda_{0}}, \\frac{1}{2 \\lambda_{0}}\\right) \\\\\n",
    "&=\\left(\\frac{2}{\\sqrt{5}}, \\frac{1}{\\sqrt{5}}\\right) \\quad \\text { or } \\quad\\left(\\frac{-2}{\\sqrt{5}}, \\frac{-1}{\\sqrt{5}}\\right)\n",
    "\\end{aligned} $$\n",
    "\n",
    "where the first point denotes a maximum (what we wanted to find) and the second a minimum. This solves our constrained optimization problem. For more details and a full solution look at [this Khan academy post](https://www.khanacademy.org/math/multivariable-calculus/applications-of-multivariable-derivatives/constrained-optimization/a/lagrange-multipliers-single-constraint)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 6.2 Recap Lagrangian  <a class=\"anchor\" id=\"recap-lagrangian\"></a>\n",
    "\n",
    "The Lagrangian is a way to repackage the individual conditions of our constrained optimization problem into a single equation. In the example above we wanted to optimize some function $f(x, y)$ under the constraint that the inputs $x$ and $y$ satisfy the equation $g(x, y) = x^2 + y^2 = c$. In our case the constant $c$ was given by 1. We know that the solution is given by those points where the gradients of $f$ and $g$ align. The Lagrangian function puts all of this into a single equation:\n",
    "\n",
    "$$ \\mathcal{L}(x, y, \\lambda)=f(x, y)-\\lambda(g(x, y)-c) $$\n",
    "\n",
    "When computing the partial derivatives of $\\mathcal{L}$ with respect to $x, y$ and $\\lambda$ and setting them to zero, we will find that they correspond exactly to the three constraints we looked at earlier. This can be summarized by simply setting the gradient of $\\mathcal{L}$ equal to the zero vector: $\\nabla \\mathcal{L} = \\mathbf{0}$. The compact Lagrangian form is often used when solving constrained optimization problem with computers because it summarizes the elaborate problem with multiple constraints into a single equation. For more details and a full solution look at [this Khan academy post](https://www.khanacademy.org/math/multivariable-calculus/applications-of-multivariable-derivatives/constrained-optimization/a/lagrange-multipliers-single-constraint)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 6.3 Dual optimization problem   <a class=\"anchor\" id=\"dual-optimization-problem\"></a>\n",
    "For the primal soft-margin SVM we considered the following optimization problem: \n",
    "$$ \\min_{\\mathbf{w}, b, \\mathbf{\\xi}} \\frac{1}{2} \\Vert \\mathbf{w} \\Vert^2 + C \\sum_{n=1}^N \\xi_n $$\n",
    "\n",
    "$$ \\text{subject to:} $$\n",
    "\n",
    "$$ y_n(\\langle \\mathbf{w}, \\mathbf{x}_n \\rangle + b) \\ge 1 - \\xi_n $$\n",
    "\n",
    "$$ \\xi_i \\ge 0 \\text{ for all } n = 1, ..., N $$\n",
    "\n",
    "Note: the optimization problem is somewhat not displayed correctly within the GitHub version of the notebook. It should look as follows:\n",
    "\n",
    "<img src=\"figures/primal_optimization_problem.png\" width=\"220\"/>\n",
    "\n",
    "To derive the corresponding Lagrangian we will introduce two Lagrange multipliers: $\\alpha_n$ for the first constraint (that all examples are classified correctly) and $\\lambda_n$ for the second constraint (non-negativity of the slack variables). The Lagrangian is then given by:\n",
    "\n",
    "$$\n",
    "\\begin{aligned}\n",
    "\\mathfrak{L}(\\boldsymbol{w}, b, \\xi, \\alpha, \\gamma)=& \\frac{1}{2}\\|\\boldsymbol{w}\\|^{2}+C \\sum_{n=1}^{N} \\xi_{n} \\\\\n",
    "& \\underbrace{-\\sum_{n=1}^{N} \\alpha_{n}\\left(y_{n}\\left(\\left\\langle\\boldsymbol{w}, \\boldsymbol{x}_{n}\\right\\rangle+b\\right)-1+\\xi_{n}\\right)}_{\\text{first constraint}} \\underbrace{-\\sum_{n=1}^{N} \\gamma_{n} \\xi_{n}}_{\\text{second constraint}}\n",
    "\\end{aligned}\n",
    "$$\n",
    "\n",
    "Next, we have to compute the partial derivatives of the Lagrangian with respect to the variables $\\mathbf{w}, b$ and $\\xi$: $\\frac{\\partial \\mathfrak{L}}{\\partial \\mathbf{w}}, \\frac{\\partial \\mathfrak{L}}{\\partial b}, \\frac{\\partial \\mathfrak{L}}{\\partial \\xi}$. When setting the first partial derivative to zero we obtain an important interim result: \n",
    "$$ \\mathbf{w} = \\sum_{n=1}^N \\alpha_n y_n \\mathbf{x}_n $$\n",
    "\n",
    "This equation states the the optimal solution for the weight vector is given by a linear combination of our training examples. After setting the other partial derivatives to zero, using the result and simplifying the equations we end up with the following optimization problem (for details see section 12.3.1 of the [Mathematics for Machine Learning book](https://mml-book.com)):\n",
    "\n",
    "$$\\min _{\\boldsymbol{\\alpha}} \\frac{1}{2} \\sum_{i=1}^{N} \\sum_{j=1}^{N} y_{i} y_{j} \\alpha_{i} \\alpha_{j}\\left\\langle\\mathbf{x}_{i}, \\mathbf{x}_{j}\\right\\rangle-\\sum_{i=1}^{N} \\alpha_{i}$$\n",
    "$$ \\text{subject to:} $$\n",
    "$$\\sum_{i=1}^{N} y_{i} \\alpha_{i}=0$$\n",
    "$$0 \\le \\alpha_{i} \\le C \\text{ for all } i=1, \\ldots, N$$\n",
    "\n",
    "Note: the optimization problem is somewhat not displayed correctly within the GitHub version of the notebook. It should look as follows:\n",
    "\n",
    "<img src=\"figures/dual_optimization_problem.png\" width=\"270\"/>\n",
    "\n",
    "This constrained quadratic optimization problem can be solved very efficiently, for example with quadratic programming techniques. One popular library for solving dual SVMs is [libsvm](https://github.com/cjlin1/libsvm) which makes use of a decomposition method to solve the problem (see [this paper](https://www.csie.ntu.edu.tw/~cjlin/papers/libsvm.pdf) for more details). However, several other approaches exist."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 7. Primal vs. dual approach <a class=\"anchor\" id=\"primal-vs-dual\"></a>\n",
    "\n",
    "Most SVM research in the last decade has been about the dual formulation. Why this is the case is not clear. Both approaches have advantages and disadvantages. In the paper [\"Training a Support Vector Machine in the Primal\"](https://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.129.3368&rep=rep1&type=pdf) Chapelle et al. mention the following hypothesis: \n",
    "\n",
    "> We believe that it is because SVMs were first introduced in their hard margin formulation (Boser et al., 1992), for which a dual optimization (because of the constraints) seems more natural. In general, however, soft margin SVMs should be preferred, even if the training data are separable: the decision boundary is more robust because more training points are taken into account (Chapelle et al., 2000). We do not pretend that primal optimization is better in general; our main motivation wasto point out that primal and dual are two sides of the same coin and that there is no reason to look always at the same side."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 8. Kernels / non-linear SVM <a class=\"anchor\" id=\"kernel-svms\"></a>\n",
    "\n",
    "### 8.1 What is a kernel? <a class=\"anchor\" id=\"what-is-a-kernel\"></a>\n",
    "If you take another look at the optimization equation of dual SVMs you will notice that it computes the inner product $\\left\\langle\\mathbf{x}_{i}, \\mathbf{x}_{j}\\right\\rangle$ between all datapoints $\\mathbf{x}_{i}, \\mathbf{x}_{j}$. A kernel is a way to compute this inner product implicitely in some (potentially very high dimensional) feature space. To be more precise: assume we have some mapping function $\\varphi$ which maps an $n$ dimensional input vector to an $m$ dimensional output vector: $\\varphi \\, : \\, \\mathbb R^n \\to \\mathbb R^m$. Given this mapping function we can compute the dot product of two vectors $\\mathbf x$ and $\\mathbf y$ in this space as follows: $\\varphi(\\mathbf x)^T \\varphi(\\mathbf y)$.\n",
    "\n",
    "A kernel is a function $k$ that gives the same result as this dot product: $k(\\mathbf x, \\mathbf y) = \\varphi(\\mathbf x)^T \\varphi(\\mathbf y)$. In other words: the kernel function is equivalent to the dot product of the mapping function."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 8.2 What are kernels good for? <a class=\"anchor\" id=\"what-are-kernels-good-for\"></a>\n",
    "\n",
    "Until now (apart from the soft-margin SVM) our SVMs, both primal and dual, are only able to classify data that is [linearly separable](https://en.wikipedia.org/wiki/Linear_separability). However, most datasets in practice won't be of this form. We need a way to classify data that is **not** linearly separable. This is where the so called **kernel trick** comes into play.\n",
    "\n",
    "Because the objective function of the dual SVM contains inner products only between datapoints $\\mathbf{x}_i, \\mathbf{x}_j$, we can easily replace this inner product (that is, $\\left\\langle\\mathbf{x}_{i}, \\mathbf{x}_{j}\\right\\rangle$ ) with some mapping function $\\varphi(\\mathbf{x}_i)^T \\varphi(\\mathbf{x}_j)$. This mapping function can be non-linear, allowing us to compute an SVM that is non-linear with respect to the input examples. The mapping function takes our input data (which is not linearly separable) and transforms it into some higher-dimensional space where it becomes linearly separable. This is illustrated in the figure below.\n",
    "\n",
    "<img src=\"figures/feature_mapping_illustration.png\" width=\"700\"/>\n",
    "\n",
    "In theory, we could use any mapping function we like. In practice, however, computing inner products is expensive. Therefore, we use mapping functions that have a corresponding kernel function. This will allow us to map the datapoints into a higher dimensional space without ever explicitely computing the (expensive) inner products.\n",
    "\n",
    "Let's take a look at an example. "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 8.3 Example <a class=\"anchor\" id=\"kernel-example\"></a>\n",
    "Note: this example was taken from [this StackExchange post](https://stats.stackexchange.com/posts/153134).\n",
    "\n",
    "We can create a simple polynomial kernel as follows: $k(\\mathbf{x}, \\mathbf{y}) = (1 + \\mathbf x^T \\mathbf y)^2$ with $\\mathbf x, \\mathbf y \\in \\mathbb R^2$. The kernel does not seem to correspond to any mapping function $\\varphi$, it's just a function that returns a real number. Our input vectors $\\mathbf{x}, \\mathbf{y}$ are 2-dimensional: $\\mathbf x = (x_1, x_2)$ and $\\mathbf y = (y_1, y_2)$. With this knowledge we can expand the kernel computation:\n",
    "\n",
    "$\\begin{align*}\n",
    "k(\\mathbf x, \\mathbf y) & = (1 + \\mathbf x^T \\mathbf y)^2 \\\\\n",
    "&= (1 + x_1 \\, y_1  + x_2 \\, y_2)^2 \\\\\n",
    " & = 1 + x_1^2 y_1^2 + x_2^2 y_2^2 + 2 x_1 y_1 + 2 x_2 y_2 + 2 x_1 x_2 y_1 y_2\n",
    "\\end{align*}$\n",
    "\n",
    "Note that this is nothing else but a dot product between two vectors $(1, x_1^2, x_2^2, \\sqrt{2} x_1, \\sqrt{2} x_2, \\sqrt{2} x_1 x_2)$ and $(1, y_1^2, y_2^2, \\sqrt{2} y_1, \\sqrt{2} y_2, \\sqrt{2} y_1 y_2)$. This can be expressed with the following mapping function: \n",
    "$$\\varphi(\\mathbf x) = \\varphi(x_1, x_2) = (1, x_1^2, x_2^2, \\sqrt{2} x_1, \\sqrt{2} x_2, \\sqrt{2} x_1 x_2)$$\n",
    "\n",
    "So the kernel $k(\\mathbf x, \\mathbf y) = (1 + \\mathbf x^T \\mathbf y)^2 = \\varphi(\\mathbf x)^T \\varphi(\\mathbf y)$ computes a dot product in 6-dimensional space without explicitly visiting this space. The generalization from an inner product to a kernel function is known as the **kernel trick**.\n",
    "\n",
    "Several popular kernel functions exist. Popular ones are, for example, the [polyomial kernel](https://en.wikipedia.org/wiki/Polynomial_kernel) or [RBF kernel](https://en.wikipedia.org/wiki/Radial_basis_function_kernel)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 8.4 Can we also use kernels in the primal SVM? <a class=\"anchor\" id=\"kernel-in-primal-svm\"></a>\n",
    "\n",
    "Yes, the kernel trick can be applied to primal SVM's, too. It's not as straightforward as with dual SVMs but still possible. Consider [this paper by Chapelle et al.](https://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.129.3368&rep=rep1&type=pdf) as an example."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 8.5 How do I choose the right kernel for my problem? <a class=\"anchor\" id=\"choosing-the-right-kernel\"></a>\n",
    "\n",
    "The problem of choosing the right kernel has been answered in [this StackExchange post](\n",
    "https://stats.stackexchange.com/questions/18030/how-to-select-kernel-for-svm). \n",
    "\n",
    "Summary:\n",
    "- Without expert knowledge, the Radial Basis Function kernel makes a good default kernel (in case you need a non-linear model)\n",
    "- The choice of the kernel and parameters can be automated by optimising a cross-validation based model selection\n",
    "- Choosing the kernel and parameters automatically is tricky, as it is very easy to overfit the model selection criterion"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 9. Sources and further reading <a class=\"anchor\" id=\"sources\"></a>\n",
    "\n",
    "The basis for this notebook is chapter 12 of the book [Mathematics for Machine Learning](https://mml-book.github.io/). I can highly recommend to read through the entire chapter to get a deeper understanding of support vector machines.\n",
    "\n",
    "Another source I liked very much is [this MIT lecture on SVMs](https://www.youtube.com/watch?v=_PwhiWxHK8o)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}