work_16_new.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Distributed Q-Learning and SARSA\n",
    "\n",
    "The goal of this assignment is to implement both single-core and distributed versions of two temporal difference (TD) reinforcement learning algorithms, Q-Learning and SARSA. In particuar, Q-Learning and SARSA will be run in a Markov Decision Process environment in order to compute policies that optimize expected infinite horizon discounted cummulative reward.\n",
    "\n",
    "The relevant content about Q-Learning and SARSA are in the following course notes from CS533.\n",
    "\n",
    "https://oregonstate.instructure.com/courses/1719746/files/74944693/download?wrap=1\n",
    "\n",
    "\n",
    "## Recap of Q-Learning and SARSA \n",
    "\n",
    "Q-Learning and SARSA appear to be extremely similar algorithms at the code level. In particular, their update formulas are as follows:\n",
    "\n",
    "\\begin{equation}\n",
    "\\mbox{SARSA:  } Q(s,a) \\leftarrow Q(s,a) + \\alpha\\cdot \\left(r + \\beta Q(s',a') - Q(s,a)\\right) \\\\\n",
    "\\mbox{Q-Learning:  } Q(s,a) \\leftarrow Q(s,a) + \\alpha\\cdot \\left(r + \\beta \\max_a' Q(s',a') - Q(s,a)\\right) \n",
    "\\end{equation}\n",
    "where $s'$ and $r$ are the next state and reward after taking action $a$ in state $s$. Further, for SARSA, $a'$ is the action that was taken by the agent in $s'$. \n",
    "\n",
    "While these updates a are similar, the algorithms are different in a fundamental way. SARSA is an on-policy algorithm, which aims to learn the value of the policy that is being used to collect data (we call this the behavior policy). In the case of policy optimization, this behavior policy is an exploration-exploitation policy. This is why $a'$ in the SARSA algorithm is based on the action actually taken by the behavior policy. \n",
    "\n",
    "Rather, Q-Learning is an off-policy algorithm, which can learn from any set of experience tuples, regardless of how they were collected. In particular, notice that Q-Learning is not sensitive to what action (if any) was taken in $s'$, while SARSA must be given the action that was actually taken in $s'$. Even if the exploration policy were completely random, Q-Learning will converge in the limit to an optimal policy under standard assumptions on the learning rate decay. \n",
    "\n",
    "The off-policy versus on-policy difference can lead the algorithms to sometimes learn different policies in practice. The assignment below will give you the opportunity to implement both SARSA and Q-Learning and to observe and explain such differences. After first implementing a single-core version of the algorithms, you will then implement a distributed version of Q-Learning. Note that the off-policy nature of Q-Learning makes it more directly compatible with distributed implementations compared to SARSA, since the experience need not come from a particular behavior policy. \n",
    "\n",
    "To complete the exercise, you need to finish all parts of this notebook. You also need to turn in a report---the notebook exercises will indicate what should be included in the report.  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "# You will need to uncomment the following pip commands if the libraries need to be installed. \n",
    "# You may get some errors related to readchar, but they should not break the project.\n",
    "\n",
    "# !pip install --user numpy\n",
    "# !pip install --user gym\n",
    "# !pip install --user tqdm"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "import ray\n",
    "import time\n",
    "import math\n",
    "from copy import deepcopy\n",
    "import matplotlib.pyplot as plt\n",
    "from random import randint, choice, uniform\n",
    "import random\n",
    "%matplotlib inline\n",
    "import pickle\n",
    "import numpy as np\n",
    "import tqdm"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "***\n",
    "\n",
    "## FrozenLake  \n",
    "<span style=\"color:red\">(The game code and transition probability is different from assignment 2)</span>  \n",
    "We will use the FrozenLake environment as the MDP environment for this experiment. This is a type of gridworld environment, whose size (number of states) can be controlled by adjusting the grid dimensions. The environment is intended to model the process of navigating a frozen lake, while avoiding falling into holes with the objective of reaching a goal location. \n",
    "\n",
    "The environment is defined as follows:\n",
    "\n",
    "- The environment is a rectangular grid of states/cells. There are four different types of cells as indicated by the following cell labels: \n",
    "\n",
    "    - S labels the starting/initial cell, always in the top left corner\n",
    "    \n",
    "    - F labels frozen cells that are safe to step on\n",
    "\n",
    "    - H labels holes and if the agent enters a hole cell there is a pentalty of -1000 and the episode ends\n",
    "\n",
    "    - G labels the goal cell and when reached gives a reward of 1000 and the episode ends\n",
    "\n",
    "- There are four possible actions (Left, Right, Down, Up). \n",
    "\n",
    "- The transition function moves the agent in the expected direction with 0.97 probability, and there is a 0.03 probability of transitioning to one of the other randomly selected directions. \n",
    "\n",
    "- There is a reward of -1 for each action taken by the agent, which is intended to encourage the agent to reach the goal as fast as possible. \n",
    "\n",
    "- Episodes end whenever the agent falls in a hole or reaches the goal. An end-of-episode is modeled by transitioning to a zero-reward terminal state (all actions lead to that state). \n",
    "\n",
    "The version we are using for assignment 3 is a modified version of the environment at the following location.   \n",
    "  \n",
    "Source: https://github.com/openai/gym/blob/master/gym/envs/toy_text/frozen_lake.py  \n",
    "\n",
    "Execute the following cell to initialize the MDP environments. (You do not need to change the code in this part.)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "import sys\n",
    "from contextlib import closing\n",
    "\n",
    "from six import StringIO, b\n",
    "\n",
    "from gym import utils\n",
    "from gym.envs.toy_text import discrete\n",
    "\n",
    "LEFT = 0\n",
    "DOWN = 1\n",
    "RIGHT = 2\n",
    "UP = 3\n",
    "np.set_printoptions(threshold=sys.maxsize, linewidth=sys.maxsize, precision = 2)\n",
    "#TransitionProb = [0.7, 0.1, 0.1, 0.1]\n",
    "#TransitionProb = [1, 0, 0, 0]\n",
    "TransitionProb = [0.97, 0.01, 0.01, 0.01]\n",
    "def generate_row(length, h_prob):\n",
    "    row = np.random.choice(2, length, p=[1.0 - h_prob, h_prob])\n",
    "    row = ''.join(list(map(lambda z: 'F' if z == 0 else 'H', row)))\n",
    "    return row\n",
    "\n",
    "\n",
    "def generate_map(shape):\n",
    "    \"\"\"\n",
    "\n",
    "    :param shape: Width x Height\n",
    "    :return: List of text based map\n",
    "    \"\"\"\n",
    "    h_prob = 0.1\n",
    "    grid_map = []\n",
    "\n",
    "    for h in range(shape[1]):\n",
    "\n",
    "        if h == 0:\n",
    "            row = 'SF'\n",
    "            row += generate_row(shape[0] - 2, h_prob)\n",
    "        elif h == 1:\n",
    "            row = 'FF'\n",
    "            row += generate_row(shape[0] - 2, h_prob)\n",
    "\n",
    "        elif h == shape[1] - 1:\n",
    "            row = generate_row(shape[0] - 2, h_prob)\n",
    "            row += 'FG'\n",
    "        elif h == shape[1] - 2:\n",
    "            row = generate_row(shape[0] - 2, h_prob)\n",
    "            row += 'FF'\n",
    "        else:\n",
    "            row = generate_row(shape[0], h_prob)\n",
    "\n",
    "        grid_map.append(row)\n",
    "        del row\n",
    "\n",
    "    return grid_map\n",
    "\n",
    "MAPS = {\n",
    "    \n",
    "    \"4x4\": [\n",
    "        \"SFFF\",\n",
    "        \"FHFH\",\n",
    "        \"FFFF\",\n",
    "        \"HFFG\"\n",
    "    ],\n",
    "    \"8x8\": [\n",
    "        \"SFFFFFFF\",\n",
    "        \"FFFFFFFF\",\n",
    "        \"FFFHFFFF\",\n",
    "        \"FFFFFHFF\",\n",
    "        \"FFFHFFFF\",\n",
    "        \"FHHFFFHF\",\n",
    "        \"FHFFHFHF\",\n",
    "        \"FFFHFFFG\"\n",
    "    ],\n",
    "    \"Dangerous Hallway\": [\n",
    "        \"SFFFFFFF\",\n",
    "        \"FFFFFFFF\",\n",
    "        \"HFHHHFFF\",\n",
    "        \"HFHHHFFF\",\n",
    "        \"HFHHHFFF\",\n",
    "        \"HFHHHFFF\",\n",
    "        \"HFFFFFFF\",\n",
    "        \"FGFFFFFF\"\n",
    "    ],\n",
    "    \"16x16\": [\n",
    "        \"SFFFFFFFFHFFFFHF\",\n",
    "        \"FFFFFFFFFFFFFHFF\",\n",
    "        \"FFFHFFFFHFFFFFFF\",\n",
    "        \"FFFFFFFFHFFFFFFF\",\n",
    "        \"FFFFFFFFFFFFFFFF\",\n",
    "        \"FFHHFFFFFFFHFFFH\",\n",
    "        \"FFFFFFFFFFFFFFFF\",\n",
    "        \"FFFFFHFFFFFFHFFF\",\n",
    "        \"FFFFFHFFFFFFFFFH\",\n",
    "        \"FFFFFFFHFFFFFFFF\",\n",
    "        \"FFFFFFFFFFFFHFFF\",\n",
    "        \"FFFFFFHFFFFFFFFF\",\n",
    "        \"FFFFFFFFHFFFFFFF\",\n",
    "        \"FFFFFFFFFHFFFFHF\",\n",
    "        \"FFFFFFFFFFHFFFFF\",\n",
    "        \"FFFHFFFFFFFFFFFG\",\n",
    "    ],\n",
    "    \n",
    "    \"32x32\": [\n",
    "        'SFFFFFFFFFFFFFFFFFFFFFFFFFHFFFFF',\n",
    "        'FFFFFFFFHFFFFFFFFFFFFFFFFFHFFFFF',\n",
    "        'FFFHFFFFFFFFHFFHFFFFFFFFFFFFFFFF',\n",
    "        'FFFFFFFFFFFFFFHFHHFHFHFFFFFHFFFH',\n",
    "        'FFFFHFFFFFFFFFFFFFFFHFHFFFFFFFHF',\n",
    "        'FFFFFHFFFFFFFFFFHFFFFFFFFFFHFFFF',\n",
    "        'FFHHFFFFHFFFFFFFFFFFFFFFFFFFFFFF',\n",
    "        'FFFHFFFFFFFFFFHFFFHFHFFFFFFFFHFF',\n",
    "        'FFFFHFFFFFFHFFFFHFHFFFFFFFFFFFFH',\n",
    "        'FFFFHHFHFFFFHFFFFFFFFFFFFFFFFFFF',\n",
    "        'FHFFFFFFFFFFHFFFFFFFFFFFHHFFFHFH',\n",
    "        'FFFHFFFHFFFFFFFFFFFFFFFFFFFFHFFF',\n",
    "        'FFFHFHFFFFFFFFHFFFFFFFFFFFFHFFHF',\n",
    "        'FFFFFFFFFFFFFFFFHFFFFFFFHFFFFFFF',\n",
    "        'FFFFFFHFFFFFFFFHHFFFFFFFHFFFFFFF',\n",
    "        'FFHFFFFFFFFFHFFFFFFFFFFHFFFFFFFF',\n",
    "        'FFFHFFFFFFFFFHFFFFHFFFFFFHFFFFFF',\n",
    "        'FFFFFFFFFFFFFFFFFFFFFFFFFFHFFFFF',\n",
    "        'FFFFFFFFHFFFFFFFHFFFFFFFFFFFFFFH',\n",
    "        'FFHFFFFFFFFFFFFFFFHFFFFFFFFFFFFF',\n",
    "        'FFFFFFFHFFFFFFFFFFFFFFFFFFFFFFFF',\n",
    "        'FFFFFFFFFFFFFFFHFFFFHFFFFFFFHFFF',\n",
    "        'FFHFFFFHFFFFFFFFFHFFFFFFFFFFFHFH',\n",
    "        'FFFFFFFFFFHFFFFHFFFFFFFFFFFFFFFF',\n",
    "        'FFFFFFFFFFFFFFFFFHHFFHHHFFFHFFFF',\n",
    "        'FFFFFFFFFFFFFFHFFFFHFFFFFFFHFFFF',\n",
    "        'FFFFFFFHFFFFFFFFFFFFFFFFFFFFFFFF',\n",
    "        'FFFFFHFFFFFFFFFFFFFFFFHFFHFFFFFF',\n",
    "        'FFFFFFFHFFFFFFFFFHFFFFFFFFFFFFFF',\n",
    "        'FFFFFFFFFFFFFFFFFFFFFFFFHFFFFFFF',\n",
    "        'FFFFFFFFFFFFFFFFFFFFFFFFHFFFFFFF',\n",
    "        'FFFFFFFFFFFFFFFHFFFFFFFFHFFFFFFG',\n",
    "    ]\n",
    "}\n",
    "\n",
    "\n",
    "def generate_random_map(size=8, p=0.8):\n",
    "    \"\"\"Generates a random valid map (one that has a path from start to goal)\n",
    "    :param size: size of each side of the grid\n",
    "    :param p: probability that a tile is frozen\n",
    "    \"\"\"\n",
    "    valid = False\n",
    "\n",
    "    # BFS to check that it's a valid path.\n",
    "    def is_valid(arr, r=0, c=0):\n",
    "        if arr[r][c] == 'G':\n",
    "            return True\n",
    "\n",
    "        tmp = arr[r][c]\n",
    "        arr[r][c] = \"#\"\n",
    "\n",
    "        # Recursively check in all four directions.\n",
    "        directions = [(1, 0), (0, 1), (-1, 0), (0, -1)]\n",
    "        for x, y in directions:\n",
    "            r_new = r + x\n",
    "            c_new = c + y\n",
    "            if r_new < 0 or r_new >= size or c_new < 0 or c_new >= size:\n",
    "                continue\n",
    "\n",
    "            if arr[r_new][c_new] not in '#H':\n",
    "                if is_valid(arr, r_new, c_new):\n",
    "                    arr[r][c] = tmp\n",
    "                    return True\n",
    "\n",
    "        arr[r][c] = tmp\n",
    "        return False\n",
    "\n",
    "    while not valid:\n",
    "        p = min(1, p)\n",
    "        res = np.random.choice(['F', 'H'], (size, size), p=[p, 1-p])\n",
    "        res[0][0] = 'S'\n",
    "        res[-1][-1] = 'G'\n",
    "        valid = is_valid(res)\n",
    "    return [\"\".join(x) for x in res]\n",
    "\n",
    "\n",
    "class FrozenLakeEnv(discrete.DiscreteEnv):\n",
    "    \"\"\"\n",
    "    Winter is here. You and your friends were tossing around a frisbee at the park\n",
    "    when you made a wild throw that left the frisbee out in the middle of the lake.\n",
    "    The water is mostly frozen, but there are a few holes where the ice has melted.\n",
    "    If you step into one of those holes, you'll fall into the freezing water.\n",
    "    At this time, there's an international frisbee shortage, so it's absolutely imperative that\n",
    "    you navigate across the lake and retrieve the disc.\n",
    "    However, the ice is slippery, so you won't always move in the direction you intend.\n",
    "    The surface is described using a grid like the following\n",
    "\n",
    "        SFFF\n",
    "        FHFH\n",
    "        FFFH\n",
    "        HFFG\n",
    "\n",
    "    S : starting point, safe\n",
    "    F : frozen surface, safe\n",
    "    H : hole, fall to your doom\n",
    "    G : goal, where the frisbee is located\n",
    "\n",
    "    The episode ends when you reach the goal or fall in a hole.\n",
    "\n",
    "    \"\"\"\n",
    "\n",
    "    metadata = {'render.modes': ['human', 'ansi']}\n",
    "\n",
    "    def __init__(self, desc=None, map_name=\"4x4\"):\n",
    "        if desc is None and map_name is None:\n",
    "            desc = generate_random_map()\n",
    "        elif desc is None:\n",
    "            desc = MAPS[map_name]\n",
    "        self.desc = desc = np.asarray(desc,dtype='c')\n",
    "        self.nrow, self.ncol = nrow, ncol = desc.shape\n",
    "        self.reward_range = (0, 1)\n",
    "\n",
    "        nA = 4\n",
    "        nS = nrow * ncol\n",
    "        self.nS = nS\n",
    "        \n",
    "        isd = np.array(desc == b'S').astype('float64').ravel()\n",
    "        isd /= isd.sum()\n",
    "\n",
    "        rew_hole = -1000\n",
    "        rew_goal = 1000\n",
    "        rew_step = -1\n",
    "        \n",
    "        exit = nrow * ncol\n",
    "        P = {s : {a : [] for a in range(nA)} for s in range(nS + 1)}\n",
    "        \n",
    "        def to_s(row, col):\n",
    "            return row*ncol + col\n",
    "        \n",
    "        def inc(row, col, a):\n",
    "            if a == LEFT:\n",
    "                col = max(col-1,0)\n",
    "            elif a == DOWN:\n",
    "                row = min(row+1,nrow-1)\n",
    "            elif a == RIGHT:\n",
    "                col = min(col+1,ncol-1)\n",
    "            elif a == UP:\n",
    "                row = max(row-1,0)\n",
    "            return (row, col)\n",
    "\n",
    "        for row in range(nrow):\n",
    "            for col in range(ncol):\n",
    "                s = to_s(row, col)\n",
    "                for a in range(4):\n",
    "                    li = P[s][a]\n",
    "                    letter = desc[row, col]\n",
    "                    if letter in b'H':\n",
    "                        li.append((1.0, exit, -1000, True))\n",
    "                    elif letter in b'G':\n",
    "                        li.append((1.0, exit, 1000, True))\n",
    "                    else:\n",
    "                        for b, p in zip([a, (a+1)%4, (a+2)%4, (a+3)%4], TransitionProb):\n",
    "                            newrow, newcol = inc(row, col, b)\n",
    "                            newstate = to_s(newrow, newcol)\n",
    "                            newletter = desc[newrow, newcol]\n",
    "                            rew = rew_step\n",
    "                            li.append((p, newstate, rew, False))\n",
    "\n",
    "        super(FrozenLakeEnv, self).__init__(nS, nA, P, isd)\n",
    "\n",
    "    def render(self, mode='human'):\n",
    "        outfile = StringIO() if mode == 'ansi' else sys.stdout\n",
    "\n",
    "        row, col = self.s // self.ncol, self.s % self.ncol\n",
    "        desc = self.desc.tolist()\n",
    "        desc = [[c.decode('utf-8') for c in line] for line in desc]\n",
    "        if self.s < self.nS:\n",
    "            desc[row][col] = utils.colorize(desc[row][col], \"red\", highlight=True)\n",
    "        else:\n",
    "            outfile.write(\"exit\\n\")\n",
    "        if self.lastaction is not None:\n",
    "            outfile.write(\"  ({})\\n\".format([\"Left\",\"Down\",\"Right\",\"Up\"][self.lastaction]))\n",
    "        else:\n",
    "            outfile.write(\"\\n\")\n",
    "        outfile.write(\"\\n\".join(''.join(line) for line in desc)+\"\\n\")\n",
    "\n",
    "        if mode != 'human':\n",
    "            with closing(outfile):\n",
    "                return outfile.getvalue()\n",
    "    "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "***\n",
    "***\n",
    "## Initializations\n",
    "\n",
    "Run the following cell to initilize maps of different sizes. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "#map_4 = (MAPS[\"4x4\"], 4)\n",
    "#map_8 = (MAPS[\"8x8\"], 8)\n",
    "#map_32 = (MAPS[\"32x32\"], 32)\n",
    "#map_50 = (generate_map((50,50)), 50)\n",
    "#map_110 = (generate_map((110,110)), 110)\n",
    "\n",
    "map_DH = (MAPS[\"Dangerous Hallway\"], 8)\n",
    "map_16 = (MAPS[\"16x16\"], 16)\n",
    "\n",
    "MAP = map_16\n",
    "map_size = MAP[1]\n",
    "run_time = {}"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Helper Function\n",
    "The following will be used to plot the performance of learning and display the learned value function. Note that both SARSA and Q-Learning produce Q-function tables, which store a value for each state-action pair. Given a Q-function $Q(s,a)$, the corresponding value function is given by $V(s)=\\max_a Q(s,a)$. This relation is used to visualize the value function learned for the environment. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "def plot_result(total_rewards, learning_num, legend):\n",
    "    print(\"\\nLearning Performance:\\n\")\n",
    "    episodes = []\n",
    "    for i in range(len(total_rewards)):\n",
    "        episodes.append(i * learning_num + 1)\n",
    "        \n",
    "    plt.figure(num = 1)\n",
    "    fig, ax = plt.subplots()\n",
    "    plt.plot(episodes, total_rewards)\n",
    "    plt.title('performance')\n",
    "    plt.legend(legend)\n",
    "    plt.xlabel(\"Episodes\")\n",
    "    plt.ylabel(\"total rewards\")\n",
    "    plt.show()\n",
    "    \n",
    "def plot_image(q_table, MAP):\n",
    "    \n",
    "    best_value = np.max(q_table, axis = 1)[:-1].reshape((map_size,map_size))\n",
    "    best_policy = np.argmax(q_table, axis = 1)[:-1].reshape((map_size,map_size))\n",
    "    \n",
    "    print(\"\\n\\nBest Q-value and Policy:\\n\")\n",
    "    fig, ax = plt.subplots()\n",
    "    im = ax.imshow(best_value)\n",
    "\n",
    "    for i in range(best_value.shape[0]):\n",
    "        for j in range(best_value.shape[1]):\n",
    "            if MAP[i][j] in 'GH':\n",
    "                arrow = MAP[i][j]\n",
    "            elif best_policy[i, j] == 0:\n",
    "                arrow = '<'\n",
    "            elif best_policy[i, j] == 1:\n",
    "                arrow = 'v'\n",
    "            elif best_policy[i, j] == 2:\n",
    "                arrow = '>'\n",
    "            elif best_policy[i, j] == 3:\n",
    "                arrow = '^'\n",
    "            if MAP[i][j] in 'S':\n",
    "                arrow = 'S ' + arrow\n",
    "            text = ax.text(j, i, arrow,\n",
    "                           ha = \"center\", va = \"center\", color = \"black\")\n",
    "            \n",
    "    cbar = ax.figure.colorbar(im, ax = ax)\n",
    "    \n",
    "    fig.tight_layout()\n",
    "    plt.show() \n",
    "    \n",
    "    "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Simulator\n",
    "\n",
    "We now create a simulator corresponding to the environment that we want to train and test the agent in. The simulator always maintains a current state, which changes using the following functions:  \n",
    "- simulator.reset(): Reset the simulator so that the current state is the initial state.   \n",
    "- simulator.step(action): take the specified action in the current state, which returns the tuple (next_state, reward, done, info) and sets the internal current stae to be next_state. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "simulator = FrozenLakeEnv(desc = MAP[0])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "***\n",
    "# Non-Distributed RL Agents\n",
    "\n",
    "In this part, you will implement Q-learning and SARSA on two maps: Dangerous Hall and the Size 16 map. \n",
    "\n",
    "$\\epsilon$-greedy exploration will be used as the exploration strategy with a fixed value of $\\epsilon$. \n",
    "\n",
    "- Suggested learning episodes ---  Dangerous Hall map = 30000 and Size 16 map = 100,000\n",
    "\n",
    "A typical way to evaluate the learning performance of an RL agent is as follows. The agent will learning using an explore-exploit policy for \"test_interval\" number of episodes. After each \"test_interval\" learning episodes, exploration is turned off and the performance of the greedy policy is evaluated. This is done by defining the policy as $\\pi(s)=\\arg\\max_a Q(s,a)$ and then running the policy for a specified number of episodes and averaging the results. \n",
    "\n",
    "- Evaluation --- by default testing is done every \"test_interval\" episodes, but testing can be turned off by setting the do_test argument to false. \n",
    "  \n",
    "In this part of the assignment you need to do the following. \n",
    "\n",
    "- Use both SARSA and Q-Learning to learn policies for each of the two maps using learning rates of 0.001 and 0.1 (using the default value of $\\epsilon$). \n",
    "\n",
    "- Use the better of the two learning rates to learning policies using SARSA and Q-Learning on each of the two maps using $\\epsilon$ values of 0.3, 0.05.\n",
    "\n",
    "### Your final report for the asignment should include the following. \n",
    "\n",
    "1. Provide the learning curves for the above experiments. Clearly label the curves by the parameters used. \n",
    "\n",
    "2. Did you observe differences for SARSA when using the two different learning rates? If there were significant differences, what were they and how can you explain them?  \n",
    "\n",
    "3. Repeat (2) for Q-Learning.\n",
    "\n",
    "4. Did you observe differences for SARSA when using different values of $\\epsilon$? If there were significant differences, what were they and how do you explain them?\n",
    "\n",
    "5. Repeat (4) for Q-Learning. \n",
    "\n",
    "6. For the map \"Dangerous Hallway\" did you observe differences in the policies learned by SARSA and Q-Learning for the two values of epsilon (there should be differences between Q-learning and SARSA for at least one value)? If you observed a difference, give your best explanation for why Q-learning and SARSA found different solutions.\n",
    "\n",
    "\n",
    "***\n",
    "## TDagent\n",
    "The parent calss of Q-learning and SARSA"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "class agent():\n",
    "    def __init__(self, epsilon, learning_rate, learning_episodes, map_size, \n",
    "                 test_interval = 100, action_space = 4, beta = 0.999, do_test = True):\n",
    "        self.beta = beta\n",
    "        self.epsilon = epsilon\n",
    "        self.test_interval = test_interval\n",
    "        self.batch_num = learning_episodes // test_interval\n",
    "        self.action_space = action_space\n",
    "        self.state_space = map_size * map_size + 1\n",
    "        self.q_table = np.zeros((self.state_space, self.action_space))\n",
    "        self.learning_rate = learning_rate\n",
    "        self.do_test = do_test\n",
    "        \n",
    "    def explore_or_exploit_policy(self, curr_state):\n",
    "        #INSERT YOUR CODE HERE\n",
    "        q_val_for_current_state = self.q_table[curr_state]#getting the action values for all action for this state\n",
    "        \n",
    "        #select an action with epsilon-greedy policy\n",
    "        if random.uniform(0,1) > self.epsilon:\n",
    "            max_action_index = self.greedy_policy(curr_state)\n",
    "            return max_action_index\n",
    "        else:\n",
    "            return (random.choice([0,1,2,3]))\n",
    "     \n",
    "    def greedy_policy(self, curr_state):\n",
    "        #INSERT YOUR CODE HERE\n",
    "        q_val_for_current_state = self.q_table[curr_state]#getting the action values for all action for this state\n",
    "        return np.argmax(q_val_for_current_state)#returning the index of the maximum action value, this index acts as the\n",
    "                                            #action \n",
    "        \n",
    "        \n",
    "    \n",
    "    def learn_and_evaluate(self):\n",
    "        total_rewards = []\n",
    "        for i in tqdm.tqdm(range(self.batch_num), desc=\"Test Number\"):\n",
    "            \n",
    "            #INSERT YOUR CODE HERE\n",
    "            self.learn()\n",
    "            average_rewards_for_100_trails = self.evaluate(self.greedy_policy)# common evaluation method for both Q learning and SARSA\n",
    "            total_rewards.append(average_rewards_for_100_trails) \n",
    "                \n",
    "         \n",
    "        return total_rewards, self.q_table\n",
    "    \n",
    "    def evaluate(self, policy_func, trials = 100, max_steps = 1000):\n",
    "        \n",
    "        total_reward = 0\n",
    "        for _ in range(trials):\n",
    "            simulator.reset()\n",
    "            done = False\n",
    "            steps = 0\n",
    "            observation, reward, done, info = simulator.step(policy_func(0))\n",
    "            total_reward += (self.beta ** steps) * reward\n",
    "            steps += 1\n",
    "            while not done and steps < 1000:\n",
    "                observation, reward, done, info = simulator.step(policy_func(observation))\n",
    "                total_reward += (self.beta ** steps) * reward\n",
    "                steps += 1\n",
    "        \n",
    "        return total_reward / trials\n",
    "    \n",
    "    def learn(self):#allow both Q Learning agent and SARSA agent to use their own learning methods\n",
    "        pass"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Q-learning Agent"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "class QLAgent(agent):\n",
    "    def __init__(self, epsilon, learning_rate, learning_episodes, map_size,\n",
    "                 test_interval = 100, action_space = 4, beta = 0.999, do_test = True):\n",
    "        super().__init__(epsilon, learning_rate, learning_episodes, map_size, \n",
    "                         test_interval = test_interval, action_space = action_space, beta = beta, do_test = do_test)\n",
    "        self.agent_name = \"Q-learning Agent\"\n",
    "        \n",
    "    def learn(self):#learning method for Q Learning\n",
    "        for _ in range(self.test_interval): #i.e., for each episode, do this.\n",
    "            simulator.reset()\n",
    "            done = False\n",
    "            current_state = 0 #start with the initial state\n",
    "            steps = 0\n",
    "            \n",
    "            #repeat until you either reach the goal state of fall in the hole\n",
    "            while done == False and steps < 1000:\n",
    "                #INSERT YOUR CODE HERE\n",
    "                \n",
    "                #print(current_state)\n",
    "                \n",
    "                action_to_take = self.explore_or_exploit_policy(current_state)\n",
    "                \n",
    "                #take this action and obtain the reward and next state\n",
    "                next_state, reward, done, _ = simulator.step(action_to_take)#automatically sets the next state\n",
    "\n",
    "                #action to take in next state\n",
    "                \n",
    "                #take argmax acti\n",
    "                next_action = self.greedy_policy(next_state)\n",
    "                \n",
    "                #print(\"here\",self.q_table[current_state][action_to_take] )\n",
    "\n",
    "                self.q_table[current_state][action_to_take] = self.q_table[current_state][action_to_take] + self.learning_rate * (reward + self.beta * self.q_table[next_state][next_action] -self.q_table[current_state][action_to_take] )\n",
    "\n",
    "                current_state = next_state\n",
    "                steps += 1"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## SARSA Agent"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "class SARSAagent(agent):\n",
    "    def __init__(self, epsilon, learning_rate, learning_episodes, map_size, \n",
    "                 test_interval = 100, action_space = 4, beta = 0.999, do_test = True):\n",
    "        super().__init__(epsilon, learning_rate, learning_episodes, map_size,\n",
    "                         test_interval = test_interval, action_space = action_space, beta = beta, do_test = do_test)\n",
    "        self.agent_name = \"SARSA Agent\"\n",
    "        \n",
    "    def learn(self):\n",
    "        for _ in range(self.test_interval):\n",
    "            #INSERT YOUR CODE HERE\n",
    "            simulator.reset()\n",
    "            done = False\n",
    "            steps = 0\n",
    "            \n",
    "            current_state = 0 #start with the initial state\n",
    "            #repeat until you either reach the goal state of fall in the hole\n",
    "            while done == False and steps<1000:\n",
    "                #INSERT YOUR CODE HERE\n",
    "                \n",
    "                #print(current_state)\n",
    "                \n",
    "                action_to_take = self.explore_or_exploit_policy(current_state)\n",
    "                \n",
    "                #take this action and obtain the reward and next state\n",
    "                next_state, reward, done, info = simulator.step(action_to_take)#automatically sets the next state\n",
    "\n",
    "                #action to take in next state\n",
    "                \n",
    "                #take argmax acti\n",
    "                next_action = self.explore_or_exploit_policy(next_state)\n",
    "                \n",
    "                #print(\"here\",self.q_table[current_state][action_to_take] )\n",
    "\n",
    "                self.q_table[current_state][action_to_take] = self.q_table[current_state][action_to_take] + self.learning_rate * (reward + self.beta * self.q_table[next_state][next_action] -self.q_table[current_state][action_to_take] )\n",
    "\n",
    "                current_state = next_state\n",
    "                steps += 1"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "***\n",
    "# Distibuted Q-Learning\n",
    "***\n",
    "## init Ray"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "2019-05-11 15:56:16,802\tINFO node.py:423 -- Process STDOUT and STDERR is being redirected to ~/ray_tmp/logs.\n",
      "2019-05-11 15:56:16,929\tINFO services.py:363 -- Waiting for redis server at 127.0.0.1:50973 to respond...\n",
      "2019-05-11 15:56:17,052\tINFO services.py:363 -- Waiting for redis server at 127.0.0.1:33314 to respond...\n",
      "2019-05-11 15:56:17,056\tINFO services.py:760 -- Starting Redis shard with 0.5 GB max memory.\n",
      "2019-05-11 15:56:17,099\tINFO services.py:1384 -- Starting the Plasma object store with 5.0 GB memory using /dev/shm.\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "{'node_ip_address': None,\n",
       " 'redis_address': '10.9.1.3:50973',\n",
       " 'object_store_address': '/home/u25771/ray_tmp/sockets/plasma_store.13',\n",
       " 'webui_url': None,\n",
       " 'raylet_socket_name': '/home/u25771/ray_tmp/sockets/raylet.7'}"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "2019-05-11 15:56:18,215\tERROR worker.py:1780 -- The log monitor on node s001-n003 failed with the following error:\n",
      "Traceback (most recent call last):\n",
      "  File \"/home/u25771/.local/lib/python3.6/site-packages/ray/log_monitor.py\", line 268, in <module>\n",
      "    log_monitor.run()\n",
      "  File \"/home/u25771/.local/lib/python3.6/site-packages/ray/log_monitor.py\", line 217, in run\n",
      "    self.update_log_filenames()\n",
      "  File \"/home/u25771/.local/lib/python3.6/site-packages/ray/log_monitor.py\", line 92, in update_log_filenames\n",
      "    log_filenames = os.listdir(self.logs_dir)\n",
      "FileNotFoundError: [Errno 2] No such file or directory: '~/ray_tmp/logs'\n",
      "\n"
     ]
    }
   ],
   "source": [
    "ray.shutdown()\n",
    "ray.init(include_webui=False, ignore_reinit_error=True, redis_max_memory=500000000, object_store_memory=5000000000, temp_dir = '~/ray_tmp')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Distibuted Q-learning Agent \n",
    "\n",
    "Below you will implement a distributed Q-Learning algorithm, where both the learning and evaluation/testing is done in a distributed way. The goal is to demonstrate a speedup in learning (in terms of wall-clock time) by creating multiple workers to collect data and do the evaluation. We have provided a schema for how to implement this. \n",
    "\n",
    "- Main Agent: Setup all the wokrers and server and turn them on.  \n",
    "- Collector: There is a simulator inside each collector. Their job is collecting exprience from the simulator, and sending the experience the Q-table server.  \n",
    "- Evaluator: There is a simulator inside each evaluator. Their job is taking a Q-table from server, then testing the performace of the greedy policy and sending the results back to the server.  \n",
    "- Server: It maintains the Q-table. It Takes expriences from collector and performs the appropriate learning updates. It also responds to the evaluator periodically when it is time to perform tests/evaluations.   \n",
    "  \n",
    "We have provided the Evaluator code, so you don't need to implement that yourself. However, understanding that code should be helpful for your implementation of other parts. \n",
    "\n",
    "After you complete the coding for this part, you need to peform evaluations on the two maps using the following number of workers: (collector workers, evaluator workers) = (2,4), (4,4), (8,4).  You can use $\\epsilon = 0.1$ and the best of the learning rates that you considered. \n",
    "\n",
    "```python\n",
    "\n",
    "\"\"\"\n",
    "\n",
    " +---------------+\n",
    " |               |\n",
    " |     Agent     |\n",
    " |               | initializes all the wokrers and server.\n",
    " |               |\n",
    " +---------------+\n",
    "                                             \n",
    "                                             \n",
    "            +-----------+         +-----------+       +-----------+         +-----------+  \n",
    "            |           |         |           |       |           |         |           |       \n",
    "            |  Worker   |         |  Worker   |       |  Worker   |         |  Worker   |       \n",
    "            |(Collector)|--------▷|(Collector)|------▷|(Collector)|--------▷|(Collector)|               \n",
    "            |           |         |           |       |           |         |           |       \n",
    "            +-----------+         +-----------+       +-----------+         +-----------+       \n",
    "                  |                     |                   |                     |            \n",
    "                  |                     |                   |                     |             \n",
    "                  +-- batch Exprience --+---------+---------+-- batch Exprience --+             \n",
    "                                                  |                                             \n",
    "                                                  |\n",
    "                                                  ▽\n",
    "                                          +----------------+                                            \n",
    "                                          |                |                                                     \n",
    "                                          |  QL Server     |  \n",
    "                                          |   (learner)    |                   \n",
    "                                          +----------------+                   \n",
    "                                                   |\n",
    "                                                   |\n",
    "                   +-------- Q-tbale -------+---------+---------+-- Q-tbale -------+ \n",
    "                   |                     |                   |                     |              \n",
    "                   |                     |                   |                     | \n",
    "                   ▽                     ▽                   ▽                     ▽\n",
    "            +-----------+         +-----------+       +-----------+         +-----------+  \n",
    "            |           |         |           |       |           |         |           |       \n",
    "            |  Worker   |         |  Worker   |       |  Worker   |         |  Worker   |       \n",
    "            |(Evaluator)|         |(Evaluator)|       |(Evaluator)|         |(Evaluator)|               \n",
    "            |           |         |           |       |           |         |           |       \n",
    "            +-----------+         +-----------+       +-----------+         +-----------+                                                                                                                      \n",
    "\"\"\"\n",
    "\n",
    "```"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "@ray.remote    \n",
    "class QLAgent_server(agent):\n",
    "    def __init__(self, epsilon, learning_rate, learning_episodes, map_size,\n",
    "                 test_interval = 100, batch_size = 100, action_space = 4, beta = 0.999, do_test = True):\n",
    "        super().__init__(epsilon, learning_rate, learning_episodes, map_size, \n",
    "                         test_interval = test_interval, action_space = action_space, beta = beta, do_test = do_test)\n",
    "        self.collector_done = False\n",
    "        self.evaluator_done = False\n",
    "        self.learning_episodes = learning_episodes\n",
    "        self.episode = 0\n",
    "        self.reuslts = []\n",
    "        self.batch_size = batch_size\n",
    "        self.privous_q_tables = []\n",
    "        self.results = [0] * (self.batch_size + 1)\n",
    "        self.reuslt_count = 0\n",
    "        \n",
    "    def learn(self, experiences):\n",
    "        #INSERT YOUR CODE HERE\n",
    "        \n",
    "        for episode in experiences:\n",
    "            \n",
    "            #self.privous_q_tables.append(self.q_table)\n",
    "            for i, step in enumerate(episode):\n",
    "                \n",
    "                state, action, reward, next_state, next_action = step\n",
    "                \n",
    "                self.q_table[state][action] = self.q_table[state][action] + self.learning_rate*(\n",
    "                    reward + (self.beta**i)*self.q_table[next_state][next_action] - self.q_table[state][action])\n",
    "            \n",
    "            self.episode += 1\n",
    "        \n",
    "        if self.do_test:\n",
    "            if self.episode // self.test_interval + 1 > len(self.privous_q_tables):\n",
    "                self.privous_q_tables.append(self.q_table)\n",
    "        return self.get_q_table()\n",
    "        \n",
    "    def get_q_table(self):\n",
    "        if self.episode >= self.learning_episodes:\n",
    "            self.collector_done = True\n",
    "            \n",
    "        return self.q_table, self.collector_done\n",
    "    def get_final_q_values(self):\n",
    "        return self.q_table\n",
    "    # evalutor\n",
    "    def add_result(self, result, num):\n",
    "        self.results[num] = result\n",
    "    \n",
    "    def get_reuslts(self):\n",
    "        return self.results, self.q_table\n",
    "    \n",
    "    def ask_evaluation(self):\n",
    "        if len(self.privous_q_tables) > self.reuslt_count:\n",
    "            num = self.reuslt_count\n",
    "            evluation_q_table = self.privous_q_tables[num]\n",
    "            self.reuslt_count += 1\n",
    "            return evluation_q_table, False, num\n",
    "        else:\n",
    "            if self.episode >= self.learning_episodes:\n",
    "                self.evaluator_done = True\n",
    "            return [], self.evaluator_done, None\n",
    "    \n",
    "@ray.remote    \n",
    "def collecting_worker(server, simulator, epsilon, action_space = 4, batch_size = 100):\n",
    "    def explore_or_exploit_policy(curr_state):\n",
    "        #INSERT YOUR CODE HERE\n",
    "        if epsilon <= uniform(0, 1):\n",
    "            return greedy_policy(curr_state)\n",
    "        else:\n",
    "            return choice([i for i in range(0, action_space)])\n",
    "        \n",
    "    def greedy_policy(curr_state):\n",
    "        #INSERT YOUR CODE HERE\n",
    "        return q_table[curr_state].argmax()\n",
    "        \n",
    "    q_table, learn_done = ray.get(server.get_q_table.remote())\n",
    "    while True:\n",
    "        if learn_done:\n",
    "            break\n",
    "        #INSERT YOUR CODE HERE\n",
    "        experiences = []\n",
    "        for i in range(batch_size):\n",
    "            simulator.reset()\n",
    "            experiences.append([])\n",
    "            done = False\n",
    "            state = 0\n",
    "            while not(done):\n",
    "                action = explore_or_exploit_policy(state)\n",
    "                next_state, reward, done, info = simulator.step(action)\n",
    "                next_action = greedy_policy(next_state)\n",
    "                experiences[i].append((state, action, reward, next_state, next_action))\n",
    "                state = next_state\n",
    "        \n",
    "        server.learn.remote(experiences)\n",
    "        q_table, learn_done = ray.get(server.get_q_table.remote())\n",
    "        \n",
    "@ray.remote\n",
    "def evaluation_worker(server, trials = 100, action_space = 4, beta = 0.999):\n",
    "    def greedy_policy(curr_state):\n",
    "        #INSERT YOUR CODE HERE\n",
    "        return q_table[curr_state].argmax()\n",
    "    \n",
    "    while True:\n",
    "        q_table, done, num = ray.get(server.ask_evaluation.remote())\n",
    "        if done:\n",
    "            break\n",
    "        if len(q_table) == 0:\n",
    "            continue\n",
    "        total_reward = 0\n",
    "        for _ in range(trials):\n",
    "            simulator.reset()\n",
    "            done = False\n",
    "            steps = 0\n",
    "            observation, reward, done, info = simulator.step(greedy_policy(0))\n",
    "            total_reward += (beta ** steps) * reward\n",
    "            steps += 1\n",
    "            while not done:\n",
    "                observation, reward, done, info = simulator.step(greedy_policy(observation))\n",
    "                total_reward += (beta ** steps) * reward\n",
    "                steps += 1\n",
    "        server.add_result.remote(total_reward / trials, num)\n",
    "\n",
    "class distributed_QL_agent():\n",
    "    def __init__(self, epsilon, learning_rate, learning_episodes, map_size, \n",
    "                 cw_num = 4, ew_num = 4, test_interval = 100, batch_size = 100,\n",
    "                 action_space = 4, beta = 0.999, do_test = True):\n",
    "        \n",
    "        self.server = QLAgent_server.remote(epsilon, learning_rate, learning_episodes, map_size, \n",
    "                                               test_interval = test_interval, batch_size = batch_size, do_test = do_test)\n",
    "        \n",
    "        self.workers_id = []\n",
    "        self.epsilon = epsilon\n",
    "        self.batch_size = batch_size\n",
    "        self.action_space = action_space\n",
    "        self.cw_num = cw_num\n",
    "        self.ew_num = ew_num\n",
    "        self.agent_name = \"Distributed Q-learning\"\n",
    "        self.do_test = do_test\n",
    "        \n",
    "    def learn_and_evaluate(self):\n",
    "        workers_id = []\n",
    "        \n",
    "        #INSERT YOUR CODE HERE\n",
    "        for i in range(self.cw_num):\n",
    "            workers_id.append(collecting_worker.remote(self.server, simulator, \n",
    "                                                       self.epsilon, self.action_space, self.batch_size // self.cw_num))\n",
    "\n",
    "#         for i in range(self.ew_num):\n",
    "#             workers_id.append(evaluation_worker.remote(self.server, trials = 100, action_space = 4, beta = 0.999))\n",
    "            \n",
    "        while len(workers_id) > 0:\n",
    "            print(len(workers_id))\n",
    "            _, workers_id = ray.wait(workers_id, len(workers_id))\n",
    "            print(len(workers_id))\n",
    "        \n",
    "        \n",
    "        return ray.get(self.server.get_reuslts.remote())\n",
    "    def get_final_values(self):\n",
    "        return ray.get(self.server.get_final_q_values.remote())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "2019-05-11 15:56:23,764\tWARNING worker.py:330 -- WARNING: Falling back to serializing objects of type <class 'numpy.dtype'> by using pickle. This may be inefficient.\n",
      "2019-05-11 15:56:23,777\tWARNING worker.py:330 -- WARNING: Falling back to serializing objects of type <class 'mtrand.RandomState'> by using pickle. This may be inefficient.\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "4\n",
      "0\n",
      "Learning time:\n",
      "\n",
      "78.37830090522766\n",
      "\n",
      "\n",
      "Best Q-value and Policy:\n",
      "\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVoAAAEYCAYAAAAdwT4RAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsnXd4XMWht9/ZVd+V5KJqS+7dliwZ\nAzYGAwFCufdiCBeD8xAghJBCCCl8N5Byk1zIvUASOjEh9I7BgI0By9i44C432XLDBRdZsnrZXfXV\nfH/saiUbWTrnzEoryfPynMea1fntzCmMzs7OeY+QUqLRaDSa7sMW6gZoNBpNf0d3tBqNRtPN6I5W\no9Fouhnd0Wo0Gk03oztajUaj6WZ0R6vRaDTdjO5oNRpNv0YI8ZIQokQIkd/utUFCiM+FEAf8/w70\nvy6EEE8JIQ4KIXYKIaa1y9zmX/+AEOI2M23QHa1Go+nvvAJcddpr9wMrpJRjgRX+MsDVwFj/chcw\nH3wdM/BH4HzgPOCPrZ2zEXRHq9Fo+jVSyjVAxWkvzwFe9f/8KnBdu9dfkz42AgOEEKnAlcDnUsoK\nKWUl8Dnf7LzPSJjKBpjFHuuQYYkDrIWb1P4miBaluKYHCXep3a3YFCuC1JK+g71RLe+NsJZrqqrA\n6/EEdYdfealDlld4Da+/dWfDbqC+3UvPSymf7yKWLKUsApBSFgkhkvyvDwWOt1uvwP/amV43RI92\ntGGJAxj6vz+1Fj4RrVR3eHDPhT6BDOEmC4W+csiaBqW6T1wcqZTvi8QeVfvj5Bpu7WQ59txjSvV2\nRHmFl805wwyvb089UC+lnB6k6jvaEbKT1w2hhw40Gk2vQgItJv6zSLF/SAD/vyX+1wuA9HbrpQGF\nnbxuCN3RajSaXobEK1sMLxZZDLTOHLgNWNTu9Vv9sw9mANX+IYYc4NtCiIH+L8G+7X/NED3e0VZ9\nuIqC+56i4L+e5sT9z1B/8HjXIU1QkV4vJ17/F4ce/gMNxUWhbo4h1qz64ynlosKtfLV/0RnW7j+o\nHKsdrz5wSrn8q80cX/9BMJvXLfiuaKXhpSuEEG8DG4DxQogCIcQPgIeBK4QQB4Ar/GWAT4HDwEHg\nX8BPAaSUFcCDQK5/+R//a4ZQGqMVQlwFPAnYgReklA93tn5LfSO12/cz9P9+iggPw1vjQTYbH/Tu\nbmRzM7LFiy2i6zE+b10t9uiYHmhV8OstWfI+EQlJDJp9OScXvMaQW39EeLzFLyn7Od1xnM2cZ2fr\nsVIYEvgGUsp5Z/jVZR2sK4G7z/A+LwEvWWmD5StaIYQdeBbfvLNJwDwhxKROQ14vttgYRLivf7fH\nOQgbFHfKKo0nSjjx+/mBclNpJQX/9XSHb1exZAk169YFypU5OVSvWmV6WxpKiylduoivn36YxvJS\nQ5lj/3ycovdfp/bwAcw6fUuXfUzV5rZ2l61cSsW6Vd1eL0D5yhxskdEkXj2H6OEjSZpzEyffex1v\nfV2X2bLT2l3+xVIqDbYb1LZbBZV2H3/ucYres76/29PoP8+OPGXsPFM5VqqoHmsVJBKvNL70BVSu\naM8DDkopDwMIId7BNwdtz5kCIjoSb3k1x3/5ONFTRuOYmUH0pJGnrBMxNAmavTQVVxCePAjPhl04\nZkzp8P0c2VlULFpE3KxZAHjy8kj54Q8NNb6lsQHX7jyqt20CKYnPPo8Rl16JLTLKUH7Ezx/Ac2Av\nVZvXUvLJQmKnTic+61zC4uK7zMZmZFP62UcMOM/XbvfuPIbecle31wsw+NIrTylHDxtB+p33GMo6\nM7Ip+/TUdg/5nrF2g9p2t3ibyN30VKDc3FTL4MSJ3d7u4ff69/emtZQsWUjc1OnEZRvf3y2NDbjy\n86jxn2dx085j+N3GzjOVYwW+fbb3w78Hys0NtQwYNtlQVvVYq2JkSKAvodLRdjSv7PzTVxJC3IXv\nDgvsCfGkP/lr6vcdoX7315Q+9S4D532b2IunnZJxzJiCZ+MuBsy5GM+GXST9/CY6+iQROTQNr9tN\nc3U1Xo8HW0w0YQON3axx+G9/IiJ5CCnXziUiMdnoNrdtl82Gc/xknOMn0+xxU7b8Ew4//iDpP7iH\n6LThnWajUtPwetw011TjrXVji4omfICxdqvUq0pUahrNFtvdmre63TZ7OOee//NAuahwKy5XQbe3\n+/T9Xf75J3z92IOk33kPUQb299d/9Z1nyXOsnWcq2OzhTLz+14Fy+VebqS3r/n2migS8uqMNYGhe\nmX/i8PMAkaOGSmGzET1pFNGTRhExLBnXmu3f7GhnZlDyxDs4zpsMQhCemgAnOm6EIzMTz86deF0u\nHFlZhhufOvd2qrdtovCdl4nNyCYu61zCBwwynAfw1tfhyt9BzfbNCLud5Dk3EZk8xFDWOWkqrj15\neN0uYqdk91i9qsROnop7Tx7NLhexGebaDWrbrYJKu731dbjzd1CzbTP493eEwf2dctPt1GzbROHb\n/vMs2/x5FipUj7UK+oq2DdPzymRjM01FZb6OE2g4UkRYwjcH9sOTB4PNRuUHK3HMyOi0EY6sbMre\nW4DX4yH1px2OYXecGzMex5jxeGs91ORtpfDtl7DHOEi+9ibCB3b9P0LRwjeoP34U5+SppHznu0QM\nTjRcN/g+RhcvXoC31kP69423W7VeVZwZ2ZQs8rU77Q7j7W7F6narYrXdJ99v29/JN5jf3984z97y\nn2dzjJ1noUT1WFtFAk19ZOzVKCodbS4wVggxEt/15s3AdzsLSNlC6fyFtNTWg81GeMpgEu6c0+G6\nzpkZVLy5lIFPXd5pIyJSUmhpaCAsPp6wuLhO1+0Ie4yDgTNnM3DmbOoKjoLN2B0ysZOzSLluHsJu\nN10nQGRSCi0N9YTFxhEWa7zdqvWqYrXdwcpbxWq9zilZJF+vvr/bn2f1Js6zUBKqYyWR/W7oQKh8\nkyqEuAZ4At/0rpeklH/pbP3IUUOlvgW359C34J49hPIW3PoTx4N6pmVmhsvFnyYYXn9k+smtQbwF\nt1tQmkcrpfwU3wRfjUajCQq+Gxb6Fz0qldFoNJquEXg7/K6979KjHa0QYLdb+1vVGN+sVLe93vqm\n2ryhO+ih/PgfKs7Gj/4ASdutn+Ml2Yr/K6uM9QQZCbT0nuYEBX1Fq9Foeh36ilaj0Wi6Ed8NC7qj\n1Wg0mm6lpZ+NmWkfrUbTBaHSSm5Y/PtTysVHt3Ao76Meqz9UtF7RGl36Arqj1Wi6oFVVOGTeHZxc\n8BpN1VWhblK/RiLwYjO89AVC3sqyN5ZRnbM5UC5f8AWVH6/rJNE5srmZlgZrT6qTzc20NFqfLG80\n762rtVwHqCnsgq2/U9lnfWF/d4eqUHW7jaB6jqnmVWmRwvDSFwh5Rxs7KwPX+l2BsnvDbpwzjanc\n2tNUVEzlgiUU/v5vNBcbc8q20lhSTNlnizn6xCM0lZnLgnmfrapT1pmRjXvX9kDZvTsP5+Sp3Z5t\nj1m3anus+H9V6lZxyg6+9EoSr7o2UG5VFdqjzN+paLbdLd4mdnzxeGA5tneZ4bpUPbrH5z9B0YI3\nguLhNUt/HDoI+ZdhkSNT8dZ4aK6owVtTi90ZRXgHopmOaGlopHZLHp61W5BS4pw1ndQ//RJbVNfz\nMFsaG3DvyqNm6yYAYqedy7B7fm3YR6vis1V1yqoo7FSyKm5VVf+vSt2qTlkVVNpts4eT9a1fBsrF\nR7fgrjKmOVTd5uG/uB/PgX1UbVxHyccfEJd1To/tMxB4ZcivAYNKyDtaAOf5k3Fv3IO3yoXzgs5t\nXe05cd9DRKSlMujWGwhPTeo60I6vH/kzkcmpJF1vzROq4rMNhlNWRWFnNaviVlX1/6rUreqUVSFU\nPlrVbfblJ+EcP6kt//eHSP/hPUSlGX8UuBV8t+D2r462V2yNc1YG7vW7cG/ag3OG8WGDhB/fgn1A\nHKX/eI3qj5fTXF5pOJt6822ExcVT9NYrVHyxjKZKw89Z8+Xn3k5YbDyF77xM+aocmqrM5b31dVRt\n2UDhWy/SVF5q2inrzMjGtWsH7j07TX/0t5pNuel2wuLiKXz7ZcpXmttm1f2lUjf49nf1lg0Uvfki\njf79bdQpq4Jqu1VQ3WZvfR3VuRsoevMlGsvKSL5uLhHJqd3Y4nZ166GD4BOZnkRLXQNhg2IJGxhr\nOBc9eRzRk8fhdXvwbNxO6TOvYot1MPjWGwhL6Nz1GTN2PDFjfZ5Q146tFL35MvYYB0nXzzXkCVXx\n2QbDKauisLOaVXGrqvp/VepWdcqqECofreo2n3zvTV9+SibJN8zr0X0mZf8bOlDSJJolavRQmf7w\njy1lG6uN3f/e8PVx7PGxhA06dZw3orjrvyn1BcewO2O/MWZp1HVQV3CUsNg4wuM7H/N078vHMXai\nIcdpb/9Stb7gKPYOttnIrfNn2l9Gt/lMdZ+Omf3dE5yp3cF0HZje5tMOmHtvPo5xXeePzX886JrE\ncRnR8unFI7te0c9Vo/b2b01ibyRyZHrXK50B1bEno+OrzgkdP2yyL6Iyxqn6jDOjdfe2/d3d48Kg\nvs3OiaHbZ75ZB/3rirbfdbQajaav0/+GDnRHq9FoehX9cdZBz3a0Qlr20aLohFVxyo58ep/l7Nf3\nTLCcDTUqitL0HI/l7LGrHNYr7sMoO2VV6GVfBniD2B4hxHjg3XYvjQL+GxgA/BBovXvkt/6nxiCE\neAD4AeAFfi6lzFFpg76i1Wg0vYpW10HQ3k/K/UAWgBDCju9hsh8C3wcel1L+rf36QohJ+B42OxkY\nAiwXQoyTUnqttkF3tBqNplfhe9x4t3VNlwGHpJRHhTjjVfMc4B0pZQPwtRDiIHAesMFqpf1rIESj\n0fR5JAKvNL4ACUKILe2Wuzp5+5uBt9uVfyaE2CmEeEkI0TrfbihwvN06Bf7XLKM7WossL3vplPKJ\n+v3sca8NUWv6Bqs2PnhKubBkG/sPLwlRazS9mRZshhegTEo5vd3yfEfvKYSIAK4F3vO/NB8YjW9Y\noQj4e+uqHcSVbjiw3NEKIdKFECuFEHuFELuFEPeqNCRYqGgSNebpCeUfqGn7gq3866ltVqE7NIc9\ntd1SglfaDC8muBrYJqUs9tUji6WUXillC/AvfMMD4LuCbT8hPw0oVNkmlSvaZuDXUsqJwAzgbv8g\nsilKX/+cyqVtPtqyd1dSsXi96cY0FRVT+f7HFP7xr4Y1iaWnuVnLVi6lQsHNahTVenuDj1ZVdWgW\nFe2fqjKwFStqyFAdq2BtM6gpMa0haDGxmGAe7YYNhBDtxQ3XA/n+nxcDNwshIoUQI4GxwGYUsDzi\nLKUswne5jZTSJYTYi28cY4+Z94mdNYWSl5cy8CrfHxPX+t2k/f4WQ9mWhkZqt+bhWZfr0yReMJ3U\nP/wSW5Qx9V5sRjaln33EgPNmAT4369BbOhveacOLl/WVCwPlJllPYoRB85ZCveCTwpR9emp+yPeM\n5VWyqqpDb0sTm3Y8Gyg3NdeROMjY9DcV7Z9KVkVzCKE7VqqaRNXtVkFC0G9YEELEAFcAP2r38qNC\niCx/lUdafyel3C2EWICvL2sG7laZcQBBmnUghBgBZAObOvjdXcBdAGEJ3zzIUaNS8Vb7fLTNNbXY\nHNGEJxrz0Z74zYNEDE1l0Pf+k/AUc5pE8LlZvRbdrHbsXDDwhra21O+nutnYX3uVelvzofDRqqoO\n7bZwzs+6O1AuLNmGy23sE5mK9k8lq6o5DNWxUtUkhkrv2Eqwb8GVUtYCg0977XudrP8X4C/Bql+5\noxVCOIGFwC+klDWn/94/MP08QNSYIR1+homdOQnXhj00V7mJu9D4PdYJd30Pz7rNlD73Go5zp+KY\nMZ2wwcY7LADnpKm49uThdbuInWLO66qCar2h8NGmzr2d6m2bKHznZWIzsonLOpfwAd1noDodb30d\n7vwd1GzbDHa7Ke2f1WzKTbdTs20ThW/7tznb/DaH4liB2v4KxnZbRdJ3HlFjFKWOVggRjq+TfVNK\n+YHV94mdNYWTzy3GW1PLsAe/bzgXPWkc0ZP8msRN2ymd/wo2p4PBt/xnl5rEQN0Z2RQvXoC31kP6\n9+/uOhAkVOt1ZmRTssiXT7vDXN5qVlV1qIKK9k8lGwzNYSiOlaomMVR6x1a0VMaP8M32fRHYK6V8\nTKURkcOSaKlrJGxQnCkfbSt2p4O4yy4k7rILafj6GNiMHyQVr6sKqvWGwkfbij3GwcCZsxk4czZ1\nBUfB1v1XH84pWSRfP8+S6lAl20r7ba43uc2hOFbB2GZQ226rSKCln0llLPtohRAXAl8Cu4BWgUHg\nXuGOiBozRI549Edn+nWn1JXGWMoF6jbgoz0ToXQdhPITlHYdaLri2HOPBd1Hmz4lXt773gzD6/+/\nScv6r49WSrmWjif2ajQajWX64xWtdh1oNJpeR195FphRerSjFWBdk2hXm3TdEmY9r/Lxv599eWqc\nHhjL0/RPpBT6ilaj0Wi6G/2EBY1Go+lGfE9Y6F+fiHRHq9Foehn975lh/WtrNL2aVev/55RyYfE2\n9h/6OESt0fRWfLMOhOGlL9DvOloVTaIZDVx3aOg03UMoFYu96TzpS+e3F5vhpS8Q8lYWv7qcik9z\nA+WSt1dR/pH5J0Y0FRVTuWAJhb//m2FNYiuNJcWUfbqYo088QlOZseyxfz5O0fvWNHSqqsJQqfdC\npZWE0CkDVXWDKvlgKS2taA6DqVk0S6vrQF/RBpG4iyZTs253oFyzbg9xs4xpbVsaGnGvy6X4kfmU\nv7qQ8NQkUv/0SyKGdf3UiZbGBmq2bqbg+acp+XAB4YlJDLvn10QOSTNU94ifP0BsxjSqNq/l6DOP\nUL5mOc011Yayzoxs3Lu2B8ru3Xk4J081lFXNq2RjM7Jx5Z+ajTXRbm9LE5u2PxNYDh9bYTir0u7h\n9z5AbOY0qjat5ejTj1Cx2vixUsmq5lW2uaWxgeptmzn+wtMUf/QuEYnJDL/7PqJSjZ3fqtutiskn\nLPR6Qv5lWPSoVJqra2mqcOGt9mB3RhGeaMyZeeK+h4hIS2XQrTcQnmpOk/j1w38mMiWVpOutaeBO\n19CVLf+Ew48/SPoP7iG6Cw2div5ONa+aVdE72m3hnJ/9s0C5sHgbLveJbm93qBSLqnmVbVbVHKpu\ntwq+Jyz0jStVo4S8owWImzmRmvV7aK50E29Gk/jjW/CszaX0H6/hOC8LxwXnGNYkps67jZqtmyh6\n8xViM7OJzZ5u2krkra/Dlb+Dmu2bEX4NXaRBDZ2K/k41r5INlVYSQqcMVMmq5q1uczA0h6rbrUJf\nGRIwSu/oaC+aTNGzS/C6ahn+0G2Gc9GTxxE92a9J3Lid0mdexRbrYPCtN3SpSYwZO56YsT4NnGvH\nVoreeBm7w0HS9XMNdbhFC9s0dCnfMa+hU1HnqeZVsqHSSkJolIGqukHVfNCUliY1h6rtVkEiaJJq\n1rHeRq/oaKMCmsRYwgdZ1CRefiFxl19Iw9fHTWkS7TEOBlwwmwEXzKb++DE487PeTyF2chYp11nX\n0PVVTWKotJIqdYdSsaiaD6bS0ozmMFiaRSu0Tu/qT1jWJFoheswQOervP7SUdZeqqfMiSqz/TbE1\nWz/offl8UdIkfm59etCxb6spMTU9R3doEhMnJcjrX/s3w+v/69zX+q8mUaPRaLoLfQuuRqPRdCN6\n1oEiQkCYzaImUXG/C2//OnA9gcq5rj/+a1TQmkSNRqPpRvRTcDUajaYH0GO0Go1G0430x+ld/Wsg\npA8hvV5OvP4vDj38BxqKi0LdnH6P3t/mCeU+a5E2w4sRhBBHhBC7hBA7hBBb/K8NEkJ8LoQ44P93\noP91IYR4SghxUAixUwgxTXV7lDtaIYRdCLFdCLFE9b3OJkqWvE9EQhJD5t3ByQWv0VRdFeom9Wv0\n/jZPyPaZCXOXySvfS6WUWe3m3N4PrJBSjgVW+MsAVwNj/ctdwHzVTQrGFe29wF6r4aJXVlD26ZZA\n+eRbqyn9cKPlxvSUj1YlX74yB1tkNIlXzyF6+EiS5tzEyfdex1tf12U2lH5U1bpVtH8qdavs747o\nKa9rKHWawd5nZmh9lI3RRYE5wKv+n18Frmv3+mvSx0ZggBAiVaUipY5WCJEG/BvwgtX3GHDRZKq/\n3BMoV6/bS/yFE02/T1NhMZXvWvTRWvB1tqfBn//6aWP5wZdeSeJV1wbK0cNGkH7nPdijorvMhtKP\nqlq3ivZPpW6V/d2enva6hlKnGax9ZhWTV7QJQogt7Za7OnhLCSwTQmxt9/tkKWURgP/fVgXgUOB4\nu2yB/zXLqH4Z9gTwX8AZBQX+jboLIDzxm/dqR49OobnaQ1O5i+aaWuyOKCIMahJbGhqpzc3Ds3YL\nUkqcs6aTOueX2KIiu842NuDKz6Nm2yaQkrhp5zH87iuxRUYZq7uxAdfuPKr9+fjs8xhxqfG8VYbf\n+wCeA3up2rSWkiULiZs6nbjscwmLM7bPVPKqdato/1TrtorqeaLS7lDqNEOJhS/DygzcgjtLSlko\nhEgCPhdC7Otk3Y4qV3IVWO5ohRD/DpRIKbcKIS4503pSyueB5wFixg7psLHxF0ykev1emirdDJg9\n2XAbTvza76O9zYKPVtHXefhvvnzKtdbyVgmlHzUYjlKr2r9Q+VFD7XUNpU4zlAR71oGUstD/b4kQ\n4kPgPKBYCJEqpSzyDw2U+FcvANLbxdOAQpX6VYYOZgHXCiGOAO8A3xJCvGHljQbMnkzVl3uoXr+P\n+AsmGM4l/OQW7APiKH32Nao/Xk5zeaXhbMpNtxMWF0/h2y9TvjKHpqoKU21OnXs7YbHxFL7zMuWr\nzOdV8NbXUb1lA0VvvkhjeaklP6rVvGrdzoxsXLt24N6z09TH4GDUbQXV8wTU2q2yv4KRDwXBfpSN\nEMIhhIht/Rn4NpAPLAZavay3AYv8Py8GbvXPPpgBVLcOMVjF8hWtlPIB4AEA/xXtfVLKW6y8V9Sw\nRFrqGggfbE6T2KGP1ulg8G1d+2hVfZ3fyL/tz19rLG+VUPpRg+Eotar9C5UfNdRe11DqNENJkG9Y\nSAY+FD4FahjwlpRyqRAiF1gghPgBcAy40b/+p8A1wEGgFvi+agOCokls19H+e2frxYwdIsc+9gNL\ndVSVGOuAGw4fxz4glrBBA055PbIwvMtsfcFR7LFxhMefOo5lVBdYV3CUsNPywZ537d6Xj2PsRMue\nUJW8at0qhLLu0znTedIRvand3UF3aBLjxyfLGc/PM7z+skuePDs0iVLKVcCqYLyXKpGj0rte6Qyo\njvV19aywYOCcYPxRP8HOq9atQijrPh0z50lvandfoT/eGaZvwdVoNL0O3dFqNBpNN6LtXYoIJGF2\naz7asEq1po5+xdhjrTuiuYP5v0Y5do35Z6C1p5+db72e9C+s3/l0/FuKk/lVnh2keKIM3Gft/8sT\n9UrVnhHZz058fUWr0Wh6HVqTqNFoNN2IlHqMVqPRaLoZgbelfxlc+9fWmOTzr586pVzgymdP2QpD\n2ZW5fzmlXFi6nX1HPgla2zpDu1V7ltVr/nxKuahoG/u/Whyi1vQMWxb89pRy6eFcjuR+0GP1SykM\nL32Bs7qjDRWqusFQeUJDqVjsq5yt261C6zzabvDRhoyQd7QnXvqC0iVbA+XCN9dQ/MEmQ9nyT5ZQ\nvb7Nt1mxLIeq1auC3cSgc+yfj1P0vjV1nqonVMVRGkrFokq7Vb2uKqhsd1nOEqo2tW93To/5aEOK\n9I3TGl36AiHvaAdePInKdj7aqi/3MvBCY2IZZ1YWnrwdgbInLw9npnFxhlc2s67gtcBysGK98WxL\nExt3zQ8shwpWGs6O+PkDxGZMo2rzWo4+8wjla5bTXFNtKKvqCVVxlA6/9wFiM6dRtWktR59+hIrV\nxtutklVtt6rX1dvSxObcpwPL4SPLDWdVttuZkYU7v+38duf3nI+2xdtE/qePBZYTO3MM1xsMekj8\n3WOE/MuwmNEpNFXV0ljuorm6FrszmogkY47RyKFpeN1umqur8Xo82KKjCRto3LdpF2HMSrs1UC5w\n5VPTUGwsawtnRsZPAuXC0u3UeIyZ1E5X55Ut/4TDjz9I+g/u6fbbeFUcpaFULKq0W9XLareFc965\n9wTKRUXbqHEVGMqqbHfUkDSa3f52+8/vnvLR2uzhTLnmV4Fy6eFcPOXHO0kED4meR9stDLxwAlXr\n9tFU6WHgbHNPV3BkZOLZtROvy4UzK6ubWhh8vPV1uPJ3ULN9M8JuJ3nOTUR2s/KvFRVHqbe+Dnf+\nDmq2bQZ/u80oFq1mVdsdSi+rynbHTs7EvXsnzW4XsRnmz+++6aPtO2OvRukdHe3sSRx76lOaa+oY\n94g506IzK5vS9xfQ4vGQ+pO7u6mFwaVoYZs6L+U7Paf8a8WZkU3JogV4az2k3WF8n4VasWi13apZ\nFVS3+5R2/8B8u0O13ar0lbFXo/SKjjZ6eCLeukbCBzsJH+Q0lY1ISaGloQF7fDxhcX3Dtxk7OYuU\n6+aFTJ1n1VHqnJJF8vXW2q2SbUXFrRoqL6vqdkcm+87vsLj4s8pH29+GDoLiozWKY2yqnPDkHZay\nVXsGK9U9br71J1Fo18HZg3YdmCN/6RN4yoPro40eM0SOeeyHxtsw53/ODh+tRqPRBBM9RqvRaDTd\njB6jVUFAmN1rKdo8sFmp6kO3Kz2WXXOWoPzxP1SoDDsAlROsTan3rlKq9oz0tzFafUWr0Wh6FZK+\n4zAwiu5oNRpNr6OfjRzojlaj0fQypB460Gg0JpFeL4VvvUT9iWOkff+nRCanhrpJvZ9+dkmrJJUR\nQgwQQrwvhNgnhNgrhJgZrIadbcjmZloaG7o9q1WH5gjGNlvVWp6N+7sV7aM9lSeBpVLKCcBUYK/Z\nNzj+4iqKP94WKBe8vpaihZsNZSvf/xTXyg2BctWiz6nJWWO47t6g3mssLaZ06SKOPPUwjeWl3Z49\nG1WHoVJDgprW8vj8Jyha8Ib1/a2gWQy1YlFrEv0IIeKA2cCLAFLKRimlaQP1oIsnUrFmX6Bc8eU+\nBl1kTJPoOG8qtbl5gXLtlp3ETM8wXHeo1HstjQ1Ub9vM8Reepvijd4lITGb43fcRlZrWrVk4O1WH\noVJDgprWcvgv7id26jSqNq6DsGSFAAAgAElEQVTj6FOPWtjf1jWLqmpJFVrtXcG6ohVCpAshVvo/\nde8WQtzrf/1PQogTQogd/uWadpkHhBAHhRD7hRBXqm6TyhjtKKAUeFkIMRXYCtwrpfS0X0kIcRdw\nF0BE0jdvZXWMSaapyuPXJNYR5owisoP1OiJi2FC8NW6aq2pocbmxxUQTNrhnNHIq2a//+icikoeQ\nPGcuEYnJhturmoWzU3UYKjWkKr66J+EcP6mt7r8/RPoP7yEqbViXeRXNoqpaUglJsO89bwZ+LaXc\nJoSIBbYKIT73/+5xKeXf2q8shJgE3AxMBoYAy4UQ46SU1m4CQK2jDQOmAfdIKTcJIZ4E7gf+0H4l\nKeXzwPMAjnGpHV7oD7pwPBVf7qep0sOgi81pEmOmZ1C3ZRfeGheO88z/xQ2Fei/lptup2baJwrdf\nJjYjm7jscwkfMKjbs62cjarDUKkhVfHW1+HetYOa7blgs5N83VwiTHyZpqJZDKViMZhDAlLKIqDI\n/7NLCLEX6OwOpjnAO1LKBuBrIcRB4DxgQyeZTlHpaAuAAill63Nn3sfX0Zpm0MUTOfLkUppq6pj4\n6HdNZWPOnUrFawtpcdeS/P9+ZLruUKj3HGPG4xgzHm+th5q8rRS+9RL2GAfJc24ifGDnnaZKFs5e\n1WEo1JCqnHzvTV/dUzJJvmGe+v42qVkMqWLRXEebIITY0q78vP8C7xsIIUYA2cAmYBbwMyHErcAW\nfFe9lfg64Y3tYgV03jF3ieWOVkp5UghxXAgxXkq5H7gM2NNVriNiRvg0iRGDY4kYbFKTODQFWd+I\nfUAc9gE9q5FTVdDZYxwMnDmbgTNnU19wFGzGPy5ZzZ6tqsNQqCFVcU6ZSvJ3blbb3wqaxdApFgWy\nxdTQQZkRe5cQwgksBH4hpawRQswHHsTXrT8I/B24Azp8Po7SNbaSJlEIkQW8AEQAh4Hv+/8idIhj\nXKqc8vRtluo6WWDu4/HpRBaGK+U1mm5H0VeghMUx0WPPPUb9ieBqEiNHpsnU//mZ4fWP3vpAl5pE\nIUQ4sATIkVI+1sHvRwBLpJRThBAPAEgp/8//uxzgT1LKkAwdIKXcAfRqD6RGo+mDBPFvjhBC4Jsd\ntbd9JyuESPWP3wJcD+T7f14MvCWEeAzfl2FjAWNzTs+AvjNMo9H0QoJ6kTwL+B6wSwjROt/tt8A8\n/6dyCRwBfgQgpdwthFiAbyi0GbhbZcYB6I5Wo9H0RoI762AtHffcn3aS+Qvwl2C1oUc7WoEk3Gbt\nkRlRA+qV6m4psb6pojl0t/kNy/F0vdIZOHalI4gt6UH64FglQOIOpYseSrMUvnBT3Weh3Ocd0cua\no4q+otVoNL2L4N+wEHJ0R6vRaHodfcVhYBRVqYwmBKza8OAp5cLibew/tCRErdF0xsZFvzulXHI0\nl8M7PgxRa8whvV5OvP4Ch/7vDzQUF3UdCGrlJpY+gO5oNRpNh5R8vJCIhESGfPcOTr77umG9Y1CQ\nwvjSBwh5R3vkhdUULW6zBB17bR0n3s81lC19fRlVS9umt5W9u5KKxes6SZxK2dIlVG9sp4JbnkPl\nl6uMZXuBYtEKofLRqrpVlZR/ClkI3fFSrVdlu8u/yMEWFdWmd7xuLiffe8OQ3jEYCGl86QuEvKNN\nuGQCZavbNIlla/aTMHu8oWzsrAxc6/MDZfeG3cTOnGy47tjMLFw722nkduXhzOjd2j4Ab0sTm3Y8\nG1gOH/vCcDZUPlpVr6ua8s961pdXUGJ6m9ix4rHAcmzPsh6p15e3vt2Dv9WR3vFnhvSOypgZNugj\nHW3IvwxzjkmmqaqWhnI3zVW1hDkjDWsSo0al4q320FxRg7emFpsjivDEAYbrjhyShtfTppGz95BG\nTlVBZ7eFc35Wm+SjsHgbLnehoezwex/Ac2AvVZvWUrJkIXFTpxOXfS5hcfG9NguKyj+FLKgdL5s9\nnKzLfhUolxzNxV1Z0O31gvp2h46+MyRglJB3tACDLxpP+RqfJjHhEnOaROeMSbg27MFb5SZ21hTT\ndTunZOLO34nX5cKZ2XMauVAp6ELlow2G11VN+Wc968uH5nip1qu63SGjj1ypGiXkQwcAif7hg7Iv\nvyLhonGmsrEXZuBal49r4x5TwwatODOzce/c7vtYNcWcz9aZkY1r1w7ce3aats+rZFXx1tdRvWUD\nRW++SGN5qWkfbSiy0LrPtluy/atk2/I9f7xU61Xd7pChhw6CT8yIBLy1jUQkOE1rEiPTk2ipbyBs\nUCxhA2NN132KRi6u92v7VAmVjzYYXlcl5Z9CFkJ3vFTrVd3ukNFHOlCjKGkSzeIclyKznr3VUra0\nxlwHfDotX1nP96dbcN378nGMnWjJcdpjWX0LrnlCtM+OzX88+JrEYeky9Te/MLz+0Z/d16UmMdT0\niitaTc/hnGB+HDvUWc3ZR1+ZtmUU3dFqNJreRz/raHvFl2EajUbTn+lZTaIAu0VN4gCn2p1FxbEx\nlrMRVQp/jxTnAzY7++gjeBQ++8kQ/vkXLdbbrTTGCiH9vOw4YS1nawpuO1rRQwcajUbT3egbFjQa\njaYb6UPzY42iO1qNRtPrENZGGHst+suwPsiXy//7lPLJE1s4sGdRiFqj6Y/sfOGBU8rl+zZT8OUH\nPdeAfnZnmFJHK4T4pRBitxAiXwjxthAiKlgNO9uQzc20NDaEuhmdEkpNYl+lN213XzjHAuiO1ocQ\nYijwc2C6lHIKYAduNvs+h/+1hsJFbSq3I6+up+C9LYayx19cRfHH2wLlgtfXUrTQ+OPXKxd9gmvN\n+kC56tNl1KxYbSir6jhtpbG0mNKlizjy1MM0lpeazptFxXGqpEmc/wRFC96wrEn8hjt4RQ6Va1d1\nexbUjrXKPgvaOVZSTOlnizny5CM0lnX/OaaKGRdtX5mdoDp0EAZECyHCgBjAmKuvHUmXjKd09f5A\nuWz1fhJmGxPLDLp4IhVr2ly2FV/uY9BFEwzX7ZiWhWd7XqBcuz2PmOxMQ1kV12dLYwPV2zZz/IWn\nKf7oXSISkxl+931EpaYZy3ub2LL+ycBy5ODnhnK+dlt3nA6/9wFiM6dRtWktR59+hIrVy2muqTaW\n/cX9xE6dRtXGdRx96lFTWfC7g3ed5g42KAFSyYLasVbZZ0E5x/71NMWLFhCRmMTwn/2aqCHGz7F9\n7/09sJzckmMoFzT62RMWLH8ZJqU8IYT4G3AMqAOWSSm/YTUWQtwF3AUQmfRN6YtzbDKNVbU0lLlp\nqq4lLDaKqGRj8gvHmGSaqjw0lrtorq4jzBll2GULEJE+FK/LTXN1NS0un68zbFD3O06//uufiEge\nQvKcuUQkJhtubys2ezjTL7g3UD55YguuamMTIVUcp+qaxEk4x09qy/79IdJ/eA9RacO6rFvFHayS\nBbVjrbLPlM6xR/9MRHIqyddZP8cm3PjrQLl832bqSo15dINCH7lSNYrljlYIMRCYA4wEqoD3hBC3\nSCnfaL+elPJ54HmA2PEpHe6+xIvGUfblVzRWeEi8xNjTFVoZdOF4Kr70uWwHXWzOZQvgyMqgdvsu\nvDUuYs7pGcdpyk23U7NtE4Vvv0xsRjZx2ecSPmCQ2aZbRsVx6q2vw52/g5ptm8FuN61JdO/aQc32\nXLDZfZ1Acqrhup2T/e5gt3l3sEoW1LyuKvvM8jl2823UbN1E4VuvEJuZTVz29B49x1TpK0MCRlGZ\n3nU58LWUshRACPEBcAHwRqepDki8dDwHHltGU3U9mY/NNZUddPFEjjy5lKaaOiY++l2zVRNzThYV\nb7+P1+0h+d6fmMo6M7IpWbQAb62HtB/c3XXAj2PMeBxjxuOt9VCTt5XCt17CHuMgec5NhA/s/v8Z\nTmn3HcbbraRJfO9NX3ZKJsk3zLOkSXRmZlP60QK8Hg9D7zTebtUsWD/WqnrIoJ1jb77sO8eum9sj\n55gyQe5ohRBXAU/i+y7pBSnlw8GtoXNUOtpjwAwhRAy+oYPLAGPfYp2GY0QCzbVNRCQ4iTTpo40Z\nkYi3rpGIwbGmXbYAEakptNQ3YB8QT1h8zzpO7TEOBs6czcCZs6kvOAq2nhlvsuo4dU7JIvn6eZY0\nic4pU0n+zs2Wsq0ouYMVst/I99A+U6m3lVPPsWM9do4pEeQvuYQQduBZ4AqgAMgVQiyWUu4JXi1d\ntEHFRyuE+DNwE9AMbAfulFKecf5I7PgUOe0ft1iqq75Z7d6K4gMJlrOhdB0MWddoOVs4K0KpbiX6\nrOtAIaz6xUwfdB3sX/g4tSXB9dFGDU2Xw376q65X9HPg97/q1EcrhJgJ/ElKeaW//ACAlPL/VNtq\nFKXeS0r5R+CPQWqLRqPR+DD3NydBCNH+0/Tz/u+GWhkKHG9XLgDOt9448+hbcDUaTa/D5MV9WRdP\nWOjoirtHPz70rCYRSYTN2uM+IiLUHhNSlapwV1Ot9cfg2BRvxAnlx/+ha6w3/sTF1tstFW2DLQpn\ntd36SA2i2Xo21HiGWsu19A2LZwGQ3q6choU5/ypo14FGo+l9BPcW3FxgrBBipBAiAt8drIuD3+gz\no4cONBpN7yLIsw6klM1CiJ8BOfimd70kpdwdvBq6Rne0Go2m9xHkEVQp5afAp8F9V+PooQONYdas\nPHWCSVHhVr7ap/WMmm5A27s0vY2+oL9T1QV6a3uPblCVvnC8QolA27uCzoHnv+T4ojaD1qFXNnB0\nwdZuzwKUvr6MqqVtWsWyd1dSsXhdJ4k2yj9dQvWGtnUrPs+has0qQ1kVVWF7rCgWg1W3WVQ1iQXP\nPEHx229Qe9B8vmLJEmrWtW1zZU4O1atWGc5/Q7O4PIfKL43nWzF7vFQ1iSr5YCkaLaOvaINL8qXj\nKV75VaBcsvoAyReP7fYsQOysDFzr8wNl94bdxM6cbCjrnJqFO6+dwm5nHo4MYwo7FVWhqmJRqW5v\nE7kbnwosRw4Z1zOqahKH3Xc/zqxp1GxYx/HHHqVypfG8IzsLT7tj5cnLwzHVuCYxNjML187TNIsG\nj7XK8VLRJKrmVetWoh/6aEP+ZVjc2KSAJrGxuo4wZ6RhTaJKFiBqVCreag/NFTV4a2qxOaIITxxg\nKBs5NA1vq8LO7VfvDTSoWFRQFaoqFlXqttnDOXfGzwPlosKtuGqMqfNUNYnCZsMxcRKOiZPwut2U\nL/2Eow8/xNCf3ENUeuf5wLGq9usGY6IJM3isQE2zqHK8VDSJqnnVupXpIx2oUULe0QIkzx5L8ZoD\nNFbUknypOU2iShbAOWMSrg178Fa5iZ01xVw2IxP3rp14XS6cU80qFq2pCoOhWFTRJKqgqkn01tfh\nydtBzdZchM1O4g1ziUgxlndkZuLZ6TtWjizzmkTnFL9m0WVOs6h6vFT0jKp51bqV6GcdbciHDsA/\nBPDFV5SsMffRXzULEHthBq51+bg27jE8bNCKc2o27h3b8ewyPmwQyGZk49q1A/eenaY+kjnGjCd1\n7q2k33kPtqhoCt96iYJX5tNUWdHtdatw8r03Of6Px2mqLCf5hnmk33k3cdnnYgs3dmtR8TtvUvDU\n4zRVlJM0dx5Df3w3cecYzzuysvHs2I5nZx6OTPPb7MzMxr1zu+8jtImnM6geL9+x2m5qiCdYedW6\nVdBDB92Ac+RgmusaiUxwEjnY0WNZgMj0JFrqGwgbFEvYwG8+AaIzIlJSkI0N2K1o+yyqCltRUSyq\n1m0FVU2iM3MqSTdaz0ek+HWD8UHQJFrIWz1eqppElbxq3ZaRQD973LiSJtEsceOT5fnPmZdzB4Nj\nFdbHl7yHVFwHfcD/eQZC5TpQvX8+dK6DvqtJtMqx+Y9TfyK4msTolHQ5+lbjmsTdf+1ck9gb6BVX\ntBqNRnMKfe9vTqfojlaj0fQ6+uDFfafojlaj0fQ+dEdrHQHYbdZGucMs5loZmVBuObu/3vpuEsei\nLGchCON+CpyYHamQVnhEkuIXIS1R1t+gJcL6/k7ZoNY7VEzqu+P5QaUP3fFlFH1Fq9FoehWCjh+J\n0JfRHa1Go+l99LMr2l5xw4JG01/Jfe+3p5RLD+dyZMsHhvPS6+XE6y9w6P/+QENxUbCb12vpbzcs\ndNnRCiFeEkKUCCHy2702SAjxuRDigP/fHrwJWqM5eyj5eCERCYkM+e4dnHz3dZqqq0LdpJ7hLLR3\nvQJcddpr9wMrpJRjgRX+smmaXPVWYgH2zl/PkQ93BcpfvbSJw+9s7yRxKvv/uY5jH+0MlA++vJEj\n724zlK14Zyk1n28MlCsXrqD607WGssHS7rVixm+q4oVVUSwq+2gV8pUffYJrzfpAueqTZdSsWG08\nv2gJri/btrvq0xxqvlhluT1GKf8iB1tUFIlXzyF6+EiSrpvLyffewFtfZyivNYm9hy47WinlGuD0\nG7PnAK/6f34VuM5K5Rt+9A55Dy6lfNtxS47SIZeNpeiLA4Fy0cqDpF46xnA+9VtjOdlOs3hy1QGS\nLzHmS3DMyMSzqa2Tr920C8d5xqQ0Ktq99ljx0R5/7nGK3nvdkhdWRbGo6qM9/uwTnHznDWoPmc87\nzsnCs63NW1y7PY+Y7Ezj+WnZ38xnGVdL7vrsscBSsCvHcL2Dv3UliVddGyhHDxtB+p0/wx4VbSiv\nNYm9B6tfhiVLKYsApJRFQoikM60ohLgLuAsgKvlUl8BFr99K6aajHPtwJ3ufWMWQb09gyFUTiUow\ndstr/LhEGirrqC9z01hVT3hsJNHJxn0FcWOTaAzk6wiPjTKcjxwxhJYaD82VNbTUeLA5oglLMKhY\nVNDutTQ24MrPo2bbJpCSuGnnMfzuK7FFGptGNvzeB/Ac2EvVprWULFlI3NTpxGWfS1hcfJdZFcXi\n8F/cj+fAPqo2rqPk4w+IyzrHcL0Aw391P7Vf7aN6wzpKF39AbPY5xE0zlo9IH4rX5aa5qpoWt1+T\nOMj4aFdE+tCAZtGXjzGct9nDybi67XbS0sO5eCqOG65bBa1J7D10+6wDKeXzwPMA8eOTT9l9wm4j\n6YKRJF0wksaqWr56fj1r5r7Mec/eyICJKYbeP/WS0RStOkRDRS2pl5m3dyVfPIbi1QdpqKgl5Vvj\nTGVjzp1M7eZ8vFVuHDONXyGBde2eqo/W54WdjHP85DYv7GMPkn7nPUSlDe8yb1WxGBQf7YRJOCb4\nfbTLPuHIow+R9qOufbQAjuwManfswlvjImaaeeWfIyuT2h07LedDRV/VJPaVK1WjWJ11UCyESAXw\n/1titQFN7gaOf5zPtt8uobagisn/dTmxoxIM54dcNpbCFQc4ueogqRcbHzZoJfVb4zj5xVcUrz5I\nism8Y2Ymno278OTmE3OuScWiRe1eyk23ExYXT+HbL1O+MoemKuN6xFa89XVUb9lA0Zsv0lheSvKc\nm4hIHmKs3QqKRW99HdW5Gyh68yUay8os+WirN2+g8HVfPuk7xn20MedkUbt1B7Xbd5oaNgjkp/nz\nO3YSk2U+Hyr6qiaxv43RWr2iXQzcBjzs/9fSo1B3PpRD1Z4iUi4eS8Zvv40jzdhH7/bEjhyMt7aR\nyEQnUQnmNYnOkYNprm0iMsFhWrMYkZZMS10DYQPjCBuooKAzod1zjBmPY8x4vLUeavK2UvjWS9hj\nHCTPuYnwgV0LpU++/wb1x4/inDyV5Bu+S8TgRHPttqhYPPnem756p2SSfMM80/WeXPAm9cf8+Rvn\nEZFgLh+RmkJLfQP2AfGExZtX/kWk+o6XPd5aPlT0SU0i/e+KtktNohDibeASIAEoBv4IfAQsAIYB\nx4AbpZRdXlrFj0+WM5+/OVAuWXeYhPNHYAvr+sJa9Rbc5hbrU4b3F5j/iN5KWA/cgltfcBR7bBzh\n8V2Pobn35eMYO9Gy19Uwp/2f4t6bj2OcsXplB6t49uYTYzDfFKdwrrSo3IJrvVrom7fgdocmMSYx\nXU64wbgmcfs/+4EmUUo57wy/uky18qRZo1TfQgOGxlZbcU4w97ieYOGcqFavQzGv6WP0sytafWeY\nRqPpVQh6bnqXEOKvQoh9QoidQogPhRAD/K+PEELUCSF2+Jfn2mXOEULsEkIcFEI8JYTo8oped7Qa\njab30XNfhn0OTJFSZgJfAQ+0+90hKWWWf/lxu9fn45uyOta/nH5D1zfoUamMVwqqG6yNWaY6apTq\njgmzPm43dbixR2p3RB5plrMAoti6qnD0+2p3Yx29xvyXiwGk9WE7qXiZIh1ey9nYQR7L2dI681/m\ntkflMTr9DdFDj9iSUi5rV9wI/Gdn6/tnWcVJKTf4y6/hu2Hrs85y+opWo9H0Lsxczfr64wQhxJZ2\ny10Wa76DUzvMkUKI7UKI1UKIi/yvDQXaX3kV+F/rFK1J1Gg0vQ6TH2rKOpt1IIRYDnR0B9TvpJSL\n/Ov8DmgG3vT/rggYJqUsF0KcA3wkhJhMx6rcLlurO1qNRtPrUH3KRnuklJd3WpcQtwH/Dlwm/fNd\npZQNQIP/561CiEPAOHxXsO3HA9OAwq7aoIcOzjJWbnrolHJhyXb2HV4SotZoNGegh74ME0JcBfwG\nuFZKWdvu9UQhhN3/8yh8X3od9jteXEKIGf7ZBrdi4IatkHa0zYqaxMYaY2pATegJpSaxxaNWd7Ax\no7W0iur+Dik9a+96BogFPj9tGtdsYKcQIg94H/hxu5uyfgK8ABwEDtHFF2EQ4qGDXfe8hnN8KklX\nZxI3dRgGpqOdwqo7FzJwUhIj/mMiCdOGmM7vfHYTMSlOxtzg8xTsfmELYTERjP9u1/ey756/gZiU\nWEZe75tIv++lXMJiwhlzc9fyjYp3lhI2eABxV8wAfC5bW3Qk8ddcaKjdFUuWEDZwIHGzZvnyOTnY\nIiOJv+QSQ3mrlOUsIWzAQAac76u3/AtfvQNndV3v8eceJzJtGPHnzCB65BjTx6rgmSeISh9G7Lnn\nEz3KXP7kn58hYlQ6ztnnEjlxtOm6T766nIjEAQy6xjcMWPL2amzRESRcN9PU+zQWF+PK3YQnfxfJ\n37uNyKGdz0hR2t/zn/Dt7+nnW9rfKnUHhR66YUFK2aHgREq5EFh4ht9tAUzdQRPSK9qsF+8k4dKJ\nnFy8jbwfvsiJtzfQWO4ynL/i7ZtJu2Ishxbms/yWd9n/2jbqyoxPz0m/fDTHVxwOlI9/cZj0bxm7\nW23oZWM4seJgoHzii0MMuXS0oayKyxbAkZWNZ0ebH9WTl4djqjHph7eliY15/wgsh45/YbheFUfp\n8HsfIDZzGlWb1nL06UeoWL2c5ppqw3UP+/X9OKdOo2bDOo4//iiVK43nUx++D8eMLFwrNlD0u8eo\nXrKS5krj0wXjL5pM9drdgXLN2j3EzZpkKNvS2IArdzMn5j9N6cIFhCclkfbLX3fZyYLi/v7F/cRO\nnUbVxnUcfepR0/s7lD7anrxhoacI6RWtsNsYOGMMA2eMoamqlmMvr2bbLc8x5fFbcE7o2sok7DZS\nZw0nddZwGirr2P3PTeTc8Aaz51/HoEld+wkGjk+gobKOulIPDVX1RMRGEpNizIU7YFwiDVV11JV5\naKyq82V7wGULEJnW5kf1uUJjCBtozBVqt4UzY+pPA+XCku3UuE8Yyqo4SlX1jMJmwzFxEo6Jfk1i\nziccfeQhhv64a02isNmIzppIdNZEvDVuqt5fSuF9D5P8u58QOSq9y7qjR6XSXO2hqdyFt8aDzRlF\nRKIxj+7Rh/5MRGoqiTfMJSLJnDNDfX9b11KG3kfbR3pQg4R81kGzp4Hy1XspXZaPsNsY/auriRll\n3MzU5G6gYMUhjn66H1uYjez7LyF+9GDD+bRLR1Kw8mvqK2pJv9zYFWkrQy4ZTeGqQzSU1zL0MnOK\nRRWXLYAjMxPPTp/P1pHVc65QFUept74Od/4OarZtBrvdlJ6xNe/ZuYOaLbkIu51EE5rEltp6ajfn\n4V67FWG3Mej7NxCRbsx5DBB3wURq1u+lucpN/EXGlZjJt9yGK3cTxa+9gjMrG+c50w1Z1lpR3t+7\ndlCzPRdsdtNaSu2jDR4h7WgPPrIE195CBl80ntH/7xqihxo/AQG2/M8KKvKLGXLpKKb//lKc6ebv\nzEm/fDRbH/6Shup6Lnn2P0xlh142hh2PrqKxqp4Ln55jKuuYmUn5ix/hdXlI+d2dprIAjuwsyha8\nh9fjIfXun3YdCBLOjGxKFi3AW+sh7Qd3G86p6hmL3/VrEjMySZprTpNY9s93aDx0jJhzMxh851zC\nU4z7jluJv2gyhc98gtdVy4i/3Go4FzNuPDHjxuP1eHBv30rxqy9jdzhIuGEu4YO6Pt8t729FLaVK\n3cr0Ic+sUULa0Q6aPZ7R912DsFsbKh76rdFM++2lhjSLZyJ+1CCaahuJTowhOiHGVDZu5CCaa5uI\nTnSYduGquGwBIlL8rtB4cz5bVaw6Sp1Tski+fp5lPaMzYypJ/3mzpXzMeZkMvvNGJTVk1LAkWuob\nCBsUS/gg449LasXucBB/4WziL5xN/fFjCJuxL6es7++pJH/H2v5SrTsYBHMebW+gSx9tMHGOS5EZ\nz9xmKavqOlDx2TZ6rZ+seUe168AsLeFq52RTSpPlrIrroG7f2ec66A4frXNQusy8/BeG19/w3n19\n30er0Wg0PY0eo9VoNJruRKJnHajQVB9O0QHzg/IA8VPU7iIbFGn9Y3R0hPWPorPHHux6pU7wjIyw\nnG08V+3wNhy2PuyRtth63aVZao/ZER7r+RGjKy1n94Qbm/J1JqSBxxadCdUxzRarh6ubnr6jr2g1\nGo2mu9EdrUaj0XQfrXeG9Sd0R6vRaHoXUva7MVqtSdT0GBsX/e6UcsnRXA7v+NBwXnq9nHj9Xxx6\n+A80FBcFu3maXkR/cx102dEKIV4SQpQIIfLbvdbhkyPN4q1Vm+fZpKBZ1IrFvkfJkveJSEhiyLw7\nOLngNZqqqwzlVM+zvkif1iRCTz6csUcwckX7Ct98ymNnT440zMlHn6T05Tep238QKzdObPjRO+Q9\nuJTybcdN53Pu+Ij1/08nwHcAACAASURBVL2Sk1sKLdW99Zlc9r2/N1De8a/t7H4zv5NEcLIAu/6x\nkUMftNmk9ry4ha/ezusk0cae+es58mGbOWz/S5s59M52Q9nKBZ/hWrEhUK76cDk1n31psNVqlK/M\nwRYZTeLVc4gePpKkOTdx8r3X8dbXdZlVPc/2nrbPvnppE4cN7rOKj5dQs25doFy5NIfqVasMZcuW\nLqF6Y1u2fEUOlWuNZY/Pf4KiBW9Qe/iApW0u/2wJ1Rva6q5YnkPVl8bqDgZn3RWtlHINUHHaa8uk\nlM3+4kaw9qjXIf/9GxzTs3GtWUfhQ3+jOmcFzVXGVW4XvX4rqZeN59iHO1l32xscfiOX+jK3oey/\nvXsjw789mgPv7+bTee+z+5Ud1JYavytoxBWjOLr860D56PKvGX7ZiG7PAqRfPoaCFYcC5YIvDpF2\nqVG941hOfNE25axw5UFSLzUmxHGcP5Xaze30jrm7iDkvw2CrocXbxI4VjwWWY3uWdR3yM/jSK0m8\n6tpAOXrYCNLvvAd7VHSXWdXzbMhlYyn64kCgXGRmn2Vn49luTWkZm5mFa1c7VeGuPJxTekiTODUL\n92l1OzJ6RpOIBFqk8aUPEIwvw+4A3j3TL/1PpLwLwD7w1BEGYbMRkzGJmIxJeF1uqhZ/xon//l9S\nfnU3kSO6VrkJu42kC0aSdMFIGqtq+er59ayZ+zLnPXsjAyZ2bmay2W0MvXAYQy8cRn1lHXnzc1l8\n/Ttc8c//YPDkpC7rHjx+MPWV9dSW1lJfWU9EXCROg4pFlSzAgHEd6R2N3X8fPy6Rxso66ss8NFTV\nEW5C7xgxfAjeGrdP7+jyYIuJJmyw8VEjmz2crMt+FSiXHM3FXWn9Ue5GUT3P4scl0lBZR32Zm8aq\nesJjI4k2qsRsr7R0m1NaRg5Jw+tpUxXae1CTGDkkzdduC3UHhb7RfxpGqaPt4MmR30BK+TzwPEDk\nsPRv7L6Wujo8W/Nwb/Sp7wZ/90YihhpXuTW5Gzi58gAnPtuDLczG5P+6nNhRxuxMje5Gji0/zOFP\nvsIWZuO8By5iwBjjBrFhl47g6BdHqCuvY8QVIw3nVLMAQy8ZxYlVh6kvryXNpKIx9ZLRFK46SENF\nLUMvG2sqGzN9CnVbduGtduM437zeMVSonmepl4ymaNUhGipqSTW5zxxTM/Hk+ZWW2eZ0g87Jmbjz\nd+J1u3Bm9qwm0TklE/cua3Wr0leGBIxiuaPt6MmRZil79S0avj5KTPZUEm69mfAkc3eN7Xwoh6o9\nRaRcPJaM334bR5rxq6v1f1pJeX4J6d8aycz/vpjYdPN39Yz89kg2/O866qsauPK5q3ssC3694yNr\naKyuZ/Yz5vWOeY+uorG6jgueut5UNmbGVCpe/oAWl4fkB+4ylQ0VqucZ+IYPdj66kqbqOmY89R1T\nWRWlpTMzm9KPFuD1eBh6Zw9rEqdmU/LBAlo8Hobc1YOaROh307ssdbTtnhx5cfsnR5olJnsqg2+5\nybLKLeXSsUy5/wpLmsRhl41ixu8vVlIsDhg1kKbaJmISY4gxqVhUyQLEjfIrGhMcRJtUNMaOHExz\nbSNRiU7zesehycj6BuwD47EP6Fl1nlVUzzPw7TNvbSORVvZZSopPsRhnXml5iqrQRDYYmsSI5BRk\nQwP2HlZxIs9CTaIQ4m3gEiABKAb+iG+WQSRQ7l9to5Tyx11VFjksXab+5l5LDZ0w5bilXCsqroNw\nm1epbhU8zQquA6/aEHxeH3UdNA6yfrymZB61nN2zZYTlLIC9ru+5DgqeeZz6guBqEuPi0uT0839m\neP2Vyx/o+5pEKeW8Dl5+sRvaotFoND762RWtvgVXo9H0OoQeo9VoNJpupA/d8WWUHu1ohRfCq6x9\n+bT3a+NPS+2ImeMPdb3SGYiwNXe90hkIVxw8GxRu/dEqkQrtBhifVWw5u8hh/EaG00l41/yXg+0p\nvNh69pwBxyxnq6Z0ffNEZ5yssP6Fk3OV2j6rzLR2nkq14fQzvWu/m3WgpTIajabX0VO34Aoh/iSE\nOCGE2OFfrmn3uweEEAeFEPuFEFe2e/0q/2sHhRD3G6lHDx1oNJreR89e0T4upfxb+xeEEJOAm4HJ\nwBBguRBinP/XzwJXAAVArhBisZRyT2cV6CtaTZ9g04enKRaP5HJ4u3HF4tnIzhdOdT2V79tMwZcf\nhKg1JvDPozW6dBNzgHeklA1Syq+Bg8B5/uWglPKwlLIReMe/bqfojlaj0fQ+WuXfRhZIEEJsabeY\nvWXxZ37l60tCiFahw1Cg/eT9Av9rZ3q9U0La0ao6M1s81vOqPtoGhbxKtr6m0XIWYNUT+Wx7t+2L\nwbXz97D5tQOdJNpY/cROtrfLrpu/m9zXvjJcd9kby6jK2Rwoly9YSeXH6zpJBAdVH+3yx/eQ+06b\nbW3VP/ax4VVjD9088sJqiha3KRWPvbaOE+/nGsqWv5lD9bJNgXLFe19QtaT79xdA5Yef4Fq9PlCu\nWrKMmuWre6RuwKyPtkxKOb3d8nz7txJCLBdC5HewzAHmA6OBLKAI+Htr7AytOtPrnRLSjvb4P57g\n5II3qD1kzZl58s/PUPbc29TvMe8ZXXbHh2z44xcUbzlhqe6Pbv2Ulb//ksLcItP597+3lM9/t5aC\n3JOms2/f8jmf/XYjxzcXW2r3xKvS2LfsRKC87/MTTLiiyz/IAEy4Kp39y9r+mO//vIDxVxi/e8w5\nKwP3ujbvrnv9bpwzJhvKtnibyPv8scByfLdxxeLJvz5J6SvWfbSTrxrK7pzCQHnPskImfdvYLJiE\nSyZQtnpfoFy2Zj8Js8cbyjovyMC9od3+2piPw+D+At8+2/fe3wPLyS05hrOO6Vl4trbpHWu35REz\nreckQkJKw0tXSCkvl1JO6WBZJKUsllJ6pZQtwL/wDQ2A70o1vd3bpAGFnbzeKSH9Mmz4L++n9qt9\nVG9aR+mSD4ideg5x084lLM6Y4CX14fuo37kf14oNVLyxCMcF03DMOoewgV1Pk7nm3bmc3HicAwv3\nsPXv6xlx1RhGXjOO6ERj97HfuHAOBesL2b1gP+se3czYa0Yx9t9H40jseprNvA/+g2PrC8lf8BVf\nPpLLuGtGMuE/RhnK3vbh1RxZd5Id7x7ki4e3MfHfhjPp2pE4E41NLUqeMABPRQOukjrqKhuIig0n\nLtXY1KDkCQOprWjAXVJHbWUDUbERhrMAUSNTaa7x0FxRg7emFpszivBEYyIgmz2cqVe0UyweMa5Y\nHPKH31C3Zx+uNeuoeO9DnOdOwzFjOmHxBs+zifHUVjTgKqnHU9lAVFw48Qa32zkmmaaqWhrK3TRX\n1RLmjCQyydg0rsiRQ/BWu337y+XB7ogmPMGclnLCjb8OlMv3baau1Ng+i0gfitflprmqmha3X4k5\nqCc1iT3zZZgQIlVK2fpcpOuB1r9si4G3hBCP4fsybCywGd8V7VghxEjgBL4vzL7bVT0h7WiFzYZj\nwiQcEybh9bgpX/YJR/72EGl3GXNmCpuN6KyJRGdNxFvjpur9pRTe9zD/v70zD4+qOv/4551JIAkJ\nIWQlJCTsBMKmgmBkkZ2wKVi3tlq1Vq1aqm39SbGionUpKu6oUFREBFMUBGQTFEFZZN9lCySQfSPJ\nZJ2c3x8zkJjEzL0zySTA/TzPPMxdvnnvDHfeOXPOeb8ndPqDNO8QWafWZDYRHhdFeFwUxTlF7Juz\ng68mL2L4exMJ7O7Yj9ZkNtFuUATtBkVQlFPMjrd389mEpUyYN4aQHnXbNJrMJqIHRRBt1259aw8L\nxn/JTfNGERrrWNthcDgdBodjySlmy5v7mRe/glvnDyMsNtDhdQN0HRHO0fVnKcwsJmZM3e9TdbqM\naMvP65MpzCymm04tgN+A7hRsPUR5bgF+cbG69c4gJhM+sd3xibX70X5l96N9VJsfLUDMyHAOrTtH\nQWYJsWO0/QK4QOCgrmRtOkpZTiFBQ2N0aX2v7UHBtoNYcwvwvc75ucnO0OKqnlh278d6Ph+fq91o\nk6hwZwnuyyLSxx41EbgfQCl1UESWAIewWcE+pJSyAojIw8AawAz8Vyl1sLY/XJVGn9510TNz1w7E\nZCZEp2dmhaUYy/a9FGzeiZhNtL57Cs0i6zb9vkBpQSlJ609wyu5H22/aYFp11O5HW1pQysm1ify8\n4gQmDxODnhxI607avvVLCko5sfY0R746iclDGPqvAQR21qjNL+XntUkcXJ6I2UMY+VQ/gjprb+nE\njI5k9bO7KMot5Y55gzTrALqNjmTtszspyi3l1nn6KwN843qSPmc51nwLEc/crVvvLBVFRRTuct6P\ntseYcFY8sxdLbil3/TdOV+zgod04PnsNZXlF9Jx1my6t73U9SX9/GRX5FsJn3KtL6yo+V/che2EC\n1sJCQh990G1xBW1dAvWBUur3dRx7Hni+lv2rgFV64jRqok39fCHFZ6p4Zgbp88zMfO8zSk+cwadf\nTwL/eAueYdoMvwG2Pr2RzANpRA7rwLVPDdXtR7vxX5tJ359B++FRDHk6Dv922qt61v9rC2n7Muk4\noh3DnhlIKx3a1dO3krI/i84jIhk9sz8B7bQ5/VcluFNLSi3l+IV4ae5yuEBQJ39KLeX4hnjr1gI0\njwyhoqgEj9Z+eATov3ZnyPzoU0oST+PTpzdBv3fOjzakU0tKC23vmV+wly6tT3QQVkspzYJ8aRao\nfSUNgGaRdlvK1i3d9n5djB1us2g0+/vj4e9mS8zLrDLMoU1ifeLVNlJF/vnRi9uFhw/g0yVGk2dm\naduyGvssuw/h3aurJn31Etyz35+mzcBITX60fp41V9s9/V0SkXFtHeprK8E99V0yUXHhmmI3N/3y\ndZ/47izt49po1LpWgmt1Yax02bGmU4Jr2X8Q7+7dNN0ndw1xfrHJb1K1DXL9Gq6U4Po1Ugluyouv\nU3K6fm0S/VuEqwHd7tN8/tpdz176NokNSYsY1/rofPp2d1rbdlCUS7Gjhujvn7xA+yHO+7x2HKKv\nf9AAfHpqH6k3aAK4t4/WLTR6H62BgYFBdQybRAMDA4OGxki0riEVznXnNDvn6VLcHys6Oq3tEJ3u\ntDbUO99pLUCkT47TWlf7aL1MNfvFtfKn7s5XMFlmOL98D8D3vfQNVlWl9zHnbRKv73TUaS1AC3G+\n6s9rgGv/12+nD3NKt8ynyKW4tXP52SQaLVoDA4OmhcJItAYGBgYNjjEYZmBgYNCwSMXllWmvaJtE\nZbWSPns+yVOfpTQ5Vbe+wlrBvn9+weab3qHgVGYDXKFBfbBR/dK39pxK5Ija/StnGzQ6CqhQ2h+X\nAA4Trd2jMV1EDtRy7O8iokREe0lWFVy1SXRVn73gSzzDggl++E6y3vuU8uw8Xfpjs9fjExFAz5mT\nODRzJcUZ2ga+XLFodNUmschFvYFBw6PDi/YS6cvV0qL9EBhTfaeIRGJbzsHpYdqkOa+R8vkCLCed\ns0lMemc2qYuds1nMW7Yek7cXAbeNp3nnaFrfNYWs9xdRYalZBVYbiR//iLlFczr9eSj+sW3p+reR\nHH5uFeUFjpOoKxaNC+74hhVPbOPM9nSn3rO5t21i6f/t5NS2DN361a8eYetnpy9ur3/7GN9/eKoO\nxS9Z+crP/LCo8nZZ+/ZxvvswUZN27auH2f5Z5bkb3znKlo+cX3BTDwtePsfqTyp/sSx+PYXlc7XN\nRJn7YgbLF1TOHPl4diYJc7M1ad95MZulC85f3J43O4dFH2hvDLzxQi5LFlR++b/3Wh4LPjhfh6KS\nHW/u5HBCpb3jrvf3sH+hQ++U+uNKS7RKqU1AbXfGa8DjuLAwcNTUafj1uorcbZs5/eZLZH+3nvLz\n2m+kqEefwK/3VeRt3cKZ118m+1vtev9JIwi4ddzF7eadogh94gFMPtqmBkXfOZBOD1TWevr3CKfv\n67fi4dvcoTZ+8S1EjerIsf8d4us7Ejj00W6KMrStdnvvsjHExEey+7PjzJ+8lq1zD1OQrn2KzUNf\nDSc2vi07PjvFuzdu5PsPfiY/XduXS6+xbdi3OuXi9v41KfQcrc3AB6BPfBh7V1eurLt3dRq9RoVq\n0vYcG86B1ZW2nwfXpNBDoyesFStb1bqLjxPoSxhx41qxZWVlsvxhVS4D47WZ+Ayd4Md3KyuT3aZV\n+Qweq82zYMSEFnyzovK+2LCykBvitdl4Aoya4MO6ryrvjXUrLYyI11aq22FUNKfWJV7cTvwmkfbD\nXaum1MVllmidGgwTkYnAWaXUXpG658Xal5X4E4CH/y/dqcRkwrdrD3y79qC8sICsdSs59epMIv/4\nCF4Rjv9Ta9gsrllJ4n/sNouR2uzvGgNXLBpNZqHj4HA6Dg7Hkl3C92/u572xq7jjwxto09Ox85jJ\nLHQZEkaXIWEUZpew4fXDzB69jns+vp62Pet2DwuP8acwq5Tz6cUUZpfi3dKTVm20G8u0jWlJQXYp\neenFFGaX4d3Sk4Bwbfo2Mf4UZttiW3JK8dIR24yZATLy4vY5lch5tM9P7tDDh7zscrLTyjifXU4L\nfw+Cw7XN9e3Uw4vcLCtZaeXkZpfj628ipK22OeFdejQnJ8tKRlo5uVlW/PxNhLXV/pHtFtuM7Cwr\nGWlWcrKstPQ30UajPrBrIEU5xVgyLBTlFNPMrxm+YfoMcZzmQh/tZYTuRCsiPsB0YJSW8+3LSrwP\nNlOZ6setxUUUHNjD+V3bwWwmdNKtNAvV1lK5qN9nt1k0mwm56RaahWm3v2ssXLFoLMkv48iaJA4s\nT8TsYWL0jKsJ7qLdfaw4v4yDa86y58skzB7CxKd7E9pFm6FJ7KgwDqxNJT+zhF5j9b/PPUeGsn9t\nGvmZpfSJ194aBug+sg2H1qXYPWG13yP1wcAxrfhxdS65GWVcP167JSXAoLF+bPo6n5yMcoaO12cc\nM3RsC75dVUhWppXh47W3Zi8wPN6H9assZGVYGTVBn/FM9LAoTm04TVFWER1Gttcd23kUqMtr1oEz\nLdqOQHvgQms2AtglIv2VUrqG7lMTPqE46TS+PXoTOuUOmgXqs69LXbLQpo/tRejN+m0WGwtXLBpX\n/nM75/Zl0XVkBPEz+xEQpc8674tpO0nem0PMqHBufL4vgVH6Wim9xrbhixkHKMwt5b4Pr9WlBVv3\nQcKMgxTmlPHgR/10aXuODWf50/uw5JZy9/yBumO7Qty4VsyZnsT5HCszP+2kSzt0vB+z/5lKXraV\nWZ/p+6U1YkILXpqWSV62lbcW6/9iGz3Bh+eeyCY3p4L3Fzs2tK9Kh5HRbPn3jxTnlhD/3mjdsV3i\nEukS0IruRKuU2g9c/B8TkUTgGqWU7vlNvrF9CL3pdk32dbXh17M3oVNuc1rfWEQO70D/J4dosjqs\nTtdREYx99hqntADdR7Vl0sy+TutDO/lRYimnZYgXLXX6sgKEdfKlpNCKf0hzWgY77s+uSkgnP0qc\n9IR1lXZdvCkqrKB1qCcBIfrKwaO7NMdSUEFQmAeBIfo+ch26NMNSaNMG6dQCdOziSWGhIjjUTHCI\nvs9JQMcAyixl+IT44BPkmg2jLi7DrgOHfrQisggYCgQBacAMpdS8KscT0ZhovdpGqnYPPObotFpR\nZtfe+NIw5+v2L1WvAx+Ta1O5XPE68BSr01pLReN5HTx4TNvqtrXhZ3Kt7t8lrwNpJK+DO1eQeTiz\nfv1om4Wq60K1r0SxOvmNS9+PVil1u4Pj0fV2NQYGBgZgdB0YGBgYNCyXzrQtrRiJ1sDAoGmhgMvM\n6+CSSbRida0bqPlZ5/v9PP7ovEfptlcHOK0FONpZWxVRbUyJ2utS7OrrlenBlf7hUE99pdDV6XXM\n+b7SQHOB01o/F/vEW7ngH+xvcm1A+INI5/yD9zZz/v2qE6NFa2BgYNDAGInWwMDAoCG5dFy5tHLJ\n2yQqq5WzCz7gxIv/oiQtxbGgnrSuWu8pq5W09+dy5smnKE3RF9vAoKFJyyjnt39OpdO1ifQblUTc\n+CS+WNVA3QTVUaBUhebHpcAln2jTVyTQLCiE8NvvIXXJx5Tl5bpF6ypZCUvxDAkh5J67yfjoE8pz\ntcUuL3BtrmZRnvP9rhYXtIUuaG165/svC1zQ5uc5Px8Y4Hye84kgL9d5ba4LWqUUk+9OYdC13hzf\nFs2OtZF8+m4YZ1Ncm6urCzf50YrIYhHZY38kisge+/5oESmqcmxOFc3VIrJfRI6LyBviyPCFJpBo\nM9d+Re72yo74rA2rydnyrSZt1sY1mJp7Ezx2Et5R7QmZdCupny/AWuw4GbmidZXc1WsxeXnR+saJ\neHVoT+BtvyHj44VUFDmO/fNf55P4ny/J35volE3iu7dtZsnjuzm5LVO3fvYt2/nkH/s5tjVbt/bF\nm3cz72+HObI1x6nrnjH5AO88dpxDP+bp1v/9xuO8+ugZ9v1YoFt738SzPDs1jZ0/FDl13b8dn8a0\nR7LYvqVYt37iuCymPpzLD1tKdGuHxadz30PZbNqsX7thcxHNPIUH7qosDY+K9OThe/V5PLiEm9y7\nlFK3KqX6KKX6AP8DllY5fOLCMaXUA1X2v4vNKKuz/VHDRrY6DivD6pPaKsOKU5LJXPUlEfc+DMDp\nN18i/Pd/wrNV3U5S7iRqxg819q1XCfhSeSOWUUow4XSTvr8477iLsw5aVZt1oKwVnN95gux1eylO\nyqL1sFhaD++FZ2BNz4PaZh1UWBXHNqez84skMk4W0GdCBH0nRdAypGY1lY+5pIb2yPeZbF96jrST\nhVw9oQ39bgrHP6RmKW31WQcVVsXBTdn8sDSV1BMW+k8MZeBNobQKraltYarp6VthVezdlMumhAzO\nnSgiblIQgyYHExBaczZJ9Qotq1Wx67t8vknIJvlECUMnBTBsSgCtQ2uW0lafdWC1KrZ9a2Hl5/mc\nPl7G6Jt8GXuzH0GhNYc3apt1YLUqtnxbzLLFhZw6Xs64yT5M+E0LQkJrzhKoPuvAalV8u7GEzxcX\ncfx4OTdN9ubm33gTGlZTW33WgdWqWL+hmIWLLRw9Vs6tU3y4/RYf2tSitekrndDenJvLqTNlvPqs\nY++Q/qOT+Glvcf1WhpmD1EDfiZrPX3N+vsuVYfZW6RlgmFLqmIhEAyuUUrHVzmsDbFRKdbNv3w4M\nVUrdX9ffb/QWrVebCMoLCyg/n0dJ6llMXt5NKsn+Ghes9y48OtLDLXHFbMK/f2faT7+Zzi/9jpLU\nXA7e/RaFR885FmOzSew6JJQ7Zl/DHz8cSE6yhVdGbSB5v+OuC5NZ6D40mD+80ZuHPr6GrOQinhu+\nmTP7HE/HMpmFnjcEcv+bPXjsk95kJhczfdg2EvdpM6I2mYW+NwQw9e0uTP+0OxlJJTw6dDcn9jru\nNzSbhX7DWvLEO9E8/2lH0pJKuG/wYX7e63iFDrNZuG54C56fE8Zbi8M5d6ac38Sd4dAebR6+ZrMw\neLg3r7wfxNzPg0k+U078wBQO7HE8FcxsFoaP8GLOBwEsTgjkzBkrcQMy2LNbm3b0SG8+nhvIyqVB\nJJ4pp1f/VHZq0Fbn4WkZ9B1+hmvHJOnWOo2+Fm2QiPxU5fEnJyIOAtKUUseq7GsvIrtF5DsRGWTf\n1xZIrnJOsn1fnTSJWQd+PXpTcGgv5fn5+PXs61hwhWMtLCbn+8Nkr9uHeJhoN3Uc3u21OzMV55ex\nf3UKu79MwuRh4sZnehHaRZsLWFF+OXu+TmXHFymYPYRbZsbQpqs2B7Ci/HJ+WpXBj0tTMXsKv3uu\nC201agEs+eVsW5nNpqUZeHgI9z7fgchu2sxOCvOtbFmZyzcJOXh4Cg+9EEF0N22eCAXnK9iwooCv\nE/IxewiPvxRMxxjt87Lzz1ewdoWFZUsseHjCUy8H0LmbNmOa8+crWPFVMQlLivDwgJf+409MjHbt\n0uVFfLq4EE9P4Y1ZreihQdu9azOWrqz8AnvrhWAys6z0d2OiVfoKFjLratGKyHqgNk/O6UqpZfbn\ntwOLqhxLAdoppbJE5GrgSxHpAdTWenfYLdAkEq1vz76kL1uC1VJIxD0PNfblNGlOz1pG4ZGztLo+\nhqjHJtC8rTYP2wt8/sRukvbmEjuqDVP+3YfAKO0epwsfP8DpvXn0Gh3C7S/0IDhau6PT/H8c4eSe\n81w1Jpg/vNSNkGjthuEAc/5+nOO7C+g3pjX3v9yRsGjtxjGvPXaGo7stXDfWn6mzIgmP1u4aNvOv\n6RzcXczQ+Bb885UQItvrc+6aPjWbfbtKGDHOh5mvBRClQ//Xv+Sye2cZ8eO9eGW2P+3ba/+43v9I\nNjt2ljJpvDfvvt6ajh20a4dd782TLyje/SiPB+39tJYid47u128JrlJqRF3HRcQDmAxcXUVTApTY\nn+8UkRNAF2wt2Igq8gjA4c/JJpFom4eEUVFSjIdfSzz89BkjX2m0GhRDu0cnIGbnen1iR4cz+bne\nmJ2wSew9JpTb/t3dKe1VY4O584WumD2c687rPzaQ+17s6JQ+Lt6fv7wc6ZT2hnEtmDYrGA8nr3vk\neG+eeSXAKf248V7MetXfKe2NE7x5+zXn4ooIS+e34bEZmcx6O4fgQDMtfEy88GSg7r/lFO63SRwB\nHFFKXewSEJFgIFspZRWRDtgGvU4qpbJFJF9EBgDbgDuBNx0FaPTBsEuB2gbDtFLfg2F6cLUEt/pg\nmC6tC+WotQ2G6cEVu8IrtQS36mCYHhpiMKylKVAN8NBuNL6ubJFLg2Ei8iGwVSlVdQrXFOBZoByw\nYrOH/cp+7Bpsi9Z6A18DjygHibRJtGgNDAwMLqLcu5SNUuoPtez7H7bpXrWd/xMQW9uxX8NItAYG\nBk0OdZmV4BqJ1sDAoOlxiZTWasWtfbQikgGcruOUIED32mP1QGPFvVJjX4mvuTFjN2TcKKVUva6K\nKiKrsV2zVjKVUg6rsxoTtyZaR4jIT42x9k9jxb1SY1+Jr7kxYzfmazaw0eiVYQYGBgaXO0aiNTAw\nMGhgmlqiff8KecWAvQAAA8xJREFUi3ulxr4SX3Njxm7M12xAE+ujNTAwMLgcaWotWgMDA4PLDiPR\nGhgYGDQwTSLRisgYETlqXxriCTfGjRSRjSJyWEQOishUd8W2xzfb/S5XuDluKxFJEJEj9tc+0I2x\nH7W/1wdEZJGIaLfh0h/rvyKSLiIHquxrLSLrROSY/d96Nz/+lbj/sb/f+0TkCxFpkOUKaotd5djf\nRUSJiJ45qgb1QKMnWhExA28DY4HuwO0i0t1N4cuBvymlYoABwENujA0wFTjsxngXeB1YbXeJ7+2u\naxCRtsBfgGvszvVm4LYGDPkhNZcZeQL4RinVGfjGvu2OuOuAWKVUL+BnYFoDxP212IhIJDAS2yoC\nBm6m0RMt0B84rpQ6qZQqBT4DJrkjsFIqRSm1y/48H1vCceiWXh+ISAQwDpjrjnhV4rYEBgPzAJRS\npUop961KaSv79rZ7gPqgwcvTWZRSm4Dq9meTgI/szz8CbnRHXKXUWqXUBXuurfzS07RBY9t5DXgc\nDSbVBvVPU0i0bYGq1u2aloaob+xrBPXF5jHpDmZju/HdXdTdAcgA5tu7LeaKiHb3bxdQSp0FZmFr\nVaUAeUqpte6IXYVQpVSK/XpSAO1LU9Qf92Cz13MLIjIROKuUcs0308BpmkKidWppiHq9ABFfbJZo\nf1VKaVvEyrV444F0pdTOho5VCx7AVcC7Sqm+QCEN8/O5Bvb+0ElAeyAcaCEiv3NH7KaCiEzH1mW1\n0E3xfIDpwFPuiGdQO00h0SYDkVW2NS0NUV+IiCe2JLtQKbXU0fn1RBwwUUQSsXWVDBORT9wUOxlI\nVkpdaLknYEu87mAEcEoplaGUKsO2tPN1bop9gTT7SqYXVjRNd1dgEbkLGA/81pFRdD3SEdsX2177\n/RYB7BKR2tbQMmggmkKi3QF0FpH2ItIM2+DIcncEti8xPA84rJR61R0xAZRS05RSEUqpaGyvd4NS\nyi0tO6VUKpAkIl3tu4YDh9wRG1uXwQAR8bG/98Nx/2DgcuAu+/O7gGV1nFtviMgY4P+AiUopx8vv\n1hNKqf1KqRClVLT9fksGrrLfBwZuotETrX2A4GFgDbYP3RKl1EE3hY8Dfo+tRbnH/oh3U+zG5BFg\noYjsA/oA/3ZHUHsrOgHYBezHdv81WHmoiCwCfgS6ikiyiNwLvAiMFJFj2EbhX3RT3LcAP2Cd/T6b\nU+cfqd/YBo2MUYJrYGBg0MA0eovWwMDA4HLHSLQGBgYGDYyRaA0MDAwaGCPRGhgYGDQwRqI1MDAw\naGCMRGtgYGDQwBiJ1sDAwKCB+X+zVBMT13KJ/gAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# simulator.reset()\n",
    "#INSERT YOUR CODE FOR INIT PARAMS HERE\n",
    "epsilon = 0.3\n",
    "learning_rate = 0.001\n",
    "learning_episodes = 100000\n",
    "test_interval = 100\n",
    "batch_size = 100\n",
    "do_test = False\n",
    "\n",
    "start_time = time.time()\n",
    "distributed_ql_agent = distributed_QL_agent(epsilon, learning_rate, learning_episodes, map_size, \n",
    "                                            test_interval = test_interval, batch_size = batch_size,\n",
    "                                            cw_num = 4, ew_num = 4, do_test = do_test)\n",
    "total_rewards, q_table = distributed_ql_agent.learn_and_evaluate()\n",
    "run_time['Distributed Q-learning agent'] = time.time() - start_time\n",
    "print(\"Learning time:\\n\")\n",
    "print(run_time['Distributed Q-learning agent'])\n",
    "if do_test:\n",
    "    plot_result(total_rewards, test_interval, [distributed_ql_agent.agent_name])\n",
    "plot_image(q_table, MAP[0])#\n",
    "q_table = distributed_ql_agent.get_final_values()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[[ 271.79   68.57  -11.65  -11.67  -10.82  -10.14   -9.55   -9.36  -11.46 -443.06   -5.29   -3.21   -3.11   -4.86 -232.9    -3.27]\n",
      " [ 296.61  181.75    6.73  -16.24  -10.26   -9.61   -9.25   -9.22  -10.3    -7.52   -3.57   -3.38   -4.46 -358.68   -4.98   -2.6 ]\n",
      " [ 318.69  156.02  -15.32 -715.39  -12.19   -8.95   -8.75  -10.95 -501.1    -6.27   -3.61   -3.44   -3.33   -4.72   -2.76   -2.67]\n",
      " [ 339.85  187.82   12.11  -14.52   -8.31   -8.33   -8.1   -10.04 -487.95   -5.64   -3.78   -3.47   -3.22   -2.93   -2.76   -2.7 ]\n",
      " [ 360.7   190.26   -9.31  -11.84   -8.11   -7.84   -7.44   -6.88   -7.33   -4.68   -4.22   -5.12   -3.15   -2.84   -2.74   -4.3 ]\n",
      " [ 379.99   67.1  -636.69 -513.91   -9.83   -7.41   -6.87   -6.32   -5.64   -4.99   -5.54 -466.51   -4.74   -2.57   -3.13 -297.55]\n",
      " [ 401.28  194.57   -7.08   -7.93   -6.68   -9.56   -6.37   -5.83   -5.03   -3.55   -4.02   -5.37   -3.98   -2.32   -2.17   -3.44]\n",
      " [ 421.58  309.19   21.85   -5.91   -8.53 -490.     -7.99   -5.1    -2.83    2.78   -0.47   -4.11 -354.17   -3.13   -1.88   -2.35]\n",
      " [ 440.73  400.49  121.57   -0.3    -7.61 -423.21   -6.74   -7.14   -2.57   13.32   15.86   -1.57   -3.48   -1.57   -2.36 -197.57]\n",
      " [ 458.59  459.31  311.65   45.73   -4.24   -6.47   -6.93 -434.64   -2.27   32.92   70.25   19.45   -2.83   -1.31   -1.12   -1.14]\n",
      " [ 474.81  491.09  475.1   249.5    23.72   -4.05   -5.62   -5.04   -2.01   22.48  127.47  115.55 -300.35   -1.31   -0.81   -0.35]\n",
      " [ 413.55  505.91  521.42  504.88  244.55   80.66 -506.56    2.91   -2.8    28.76  213.05  344.39   49.83    0.01   -0.34    4.43]\n",
      " [ 125.07  472.83  536.48  552.75  549.27  517.16  497.06  220.86 -352.88   -2.62  125.08  486.7    82.16   17.99    0.49   30.45]\n",
      " [   8.09  214.13  529.87  569.43  588.28  607.29  626.85  624.82  307.27 -340.46  120.63  647.48  422.85  285.26 -194.35  142.6 ]\n",
      " [  -3.87   12.23  197.59  445.08  528.1   620.03  649.36  671.08  688.2   659.92 -807.53  776.04  814.21  843.36  860.87  565.65]\n",
      " [  -3.95   -3.67   -7.59 -802.26   92.84  377.74  663.12  694.99  720.15  749.77  781.3   833.03  869.96  909.8   953.34 1000.  ]]\n"
     ]
    }
   ],
   "source": [
    "table = np.zeros((16,16))\n",
    "\n",
    "for i in range(16):\n",
    "    for j in range(16):\n",
    "        table[i][j]=q_table[i*16+j].max()\n",
    "np.set_printoptions(suppress=True)\n",
    "print(table)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Comparison of different approachs\n",
    "\n",
    "Run the following cell to compare the running time of different approaches. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "from copy import deepcopy\n",
    "temp_dict = deepcopy(run_time)\n",
    "print(\"All:\")\n",
    "for _ in range(len(temp_dict)):\n",
    "    min_v = float('inf')\n",
    "    for k, v in temp_dict.items():\n",
    "        if v is None:\n",
    "            continue\n",
    "        if v < min_v:\n",
    "            min_v = v\n",
    "            name = k\n",
    "    temp_dict[name] = float('inf')\n",
    "    print(name + \": \" + str(min_v))\n",
    "    print()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Include in your report the following: \n",
    "\n",
    "7. Show the value functions learned by the distributed methods for the best policies learned with $\\epsilon$ equal to 0.3 and compare to those of the single-core method. Run the algorithm for the recommended number of episodes for each of the maps. Did the approaches produce similar results? \n",
    "\n",
    "8. Provide and compare the timing results for the single-core and distributed experiments, including the time to do the evaluations during learning. Describe the trends you observe as the number of workers increases. \n",
    "\n",
    "9. Provide and compare the timing results for the single-core and distributed experiments with the evaluation procedure turned off. That is, here you only want to provide timing results for the learning process without any interleaved evaluation. Compare the results with (8). "
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (Intel, 2019 update 2)",
   "language": "python",
   "name": "c009-intel_distribution_of_python_3_2019u2"
  },
  "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.6.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}