{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Pixels-to-Control Learning\n", "\n", "> In this post, We will take a hands-on-lab of Pixels-to-Control Learning which is related on Deep Reinforcement Learning. This tutorial is the software lab, which is part of \"introduction to deep learning 2021\" offered from MIT 6.S191.\n", "\n", "- toc: true \n", "- badges: true\n", "- comments: true\n", "- author: Chanseok Kang\n", "- categories: [Python, Reinforcement_Learning, Tensorflow, MIT]\n", "- image: images/Pong-v0.gif" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Copyright Information\n", "\n", "Copyright 2021 MIT 6.S191 Introduction to Deep Learning. All Rights Reserved.\n", "\n", "Licensed under the MIT License. You may not use this file except in compliance\n", "with the License. Use and/or modification of this code outside of 6.S191 must\n", "reference:\n", "\n", "© MIT 6.S191: Introduction to Deep Learning\n", "http://introtodeeplearning.com" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Reinforcement Learning\n", "\n", "Reinforcement learning (RL) is a subset of machine learning which poses learning problems as interactions between agents and environments. It often assumes agents have no prior knowledge of a world, so they must learn to navigate environments by optimizing a reward function. Within an environment, an agent can take certain actions and receive feedback, in the form of positive or negative rewards, with respect to their decision. As such, an agent's feedback loop is somewhat akin to the idea of \"trial and error\", or the manner in which a child might learn to distinguish between \"good\" and \"bad\" actions.\n", "\n", "In practical terms, our RL agent will interact with the environment by taking an action at each timestep, receiving a corresponding reward, and updating its state according to what it has \"learned\". \n", "\n", "![alt text](image/reinforcement-learning-fig1-700.jpg)\n", "\n", "While the ultimate goal of reinforcement learning is to teach agents to act in the real, physical world, games provide a convenient proving ground for developing RL algorithms and agents. Games have some properties that make them particularly well suited for RL: \n", "\n", "1. In many cases, games have perfectly describable environments. For example, all rules of chess can be formally written and programmed into a chess game simulator;\n", "2. Games are massively parallelizable. Since they do not require running in the real world, simultaneous environments can be run on large data clusters; \n", "3. Simpler scenarios in games enable fast prototyping. This speeds up the development of algorithms that could eventually run in the real-world; and\n", "4. ... Games are fun! \n", "\n", "In previous labs, we have explored both supervised (with LSTMs, CNNs) and unsupervised / semi-supervised (with VAEs) learning tasks. Reinforcement learning is fundamentally different, in that we are training a deep learning algorithm to govern the actions of our RL agent, that is trying, within its environment, to find the optimal way to achieve a goal. The goal of training an RL agent is to determine the best next step to take to earn the greatest final payoff or return. In this lab, we focus on building a reinforcement learning algorithm to master two different environments with varying complexity. \n", "\n", "1. **Cartpole**: Balance a pole, protruding from a cart, in an upright position by only moving the base left or right. Environment with a low-dimensional observation space.\n", "2. [**Pong**](https://en.wikipedia.org/wiki/Pong): Beat your competitors (whether other AI or humans!) at the game of Pong. Environment with a high-dimensional observation space -- learning directly from raw pixels.\n", "\n", "Let's get started! First we'll import TensorFlow, the course package, and some dependencies." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "import tensorflow as tf\n", "import numpy as np\n", "import base64, io, time, gym\n", "import IPython, functools\n", "import matplotlib.pyplot as plt\n", "from tqdm import tqdm" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "To prevent memory allocation in GPU, we need to limit the memory usage." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "1 Physical GPUs, 1 Logical GPUs\n" ] } ], "source": [ "gpus = tf.config.experimental.list_physical_devices('GPU')\n", "if gpus:\n", " try:\n", " # Currently, memory growth needs to be the same across GPUs\n", " for gpu in gpus:\n", " tf.config.experimental.set_memory_growth(gpu, True)\n", " logical_gpus = tf.config.experimental.list_logical_devices('GPU')\n", " print(len(gpus), \"Physical GPUs,\", len(logical_gpus), \"Logical GPUs\")\n", " except RuntimeError as e:\n", " # Memory growth must be set before GPUs have been initialized\n", " print(e)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Before we dive in, let's take a step back and outline our approach, which is generally applicable to reinforcement learning problems in general:\n", "\n", "1. **Initialize our environment and our agent**: here we will describe the different observations and actions the agent can make in the environemnt.\n", "2. **Define our agent's memory**: this will enable the agent to remember its past actions, observations, and rewards.\n", "3. **Define a reward function**: describes the reward associated with an action or sequence of actions.\n", "4. **Define the learning algorithm**: this will be used to reinforce the agent's good behaviors and discourage bad behaviors.\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Part 1: Cartpole\n", "\n", "### 3.1 Define the Cartpole environment and agent\n", "\n", "#### Environment \n", "\n", "In order to model the environment for both the Cartpole and Pong tasks, we'll be using a toolkit developed by OpenAI called [OpenAI Gym](https://gym.openai.com/). It provides several pre-defined environments for training and testing reinforcement learning agents, including those for classic physics control tasks, Atari video games, and robotic simulations. To access the Cartpole environment, we can use `env = gym.make(\"CartPole-v0\")`, which we gained access to when we imported the `gym` package. We can instantiate different [environments](https://gym.openai.com/envs/#classic_control) by passing the enivronment name to the `make` function.\n", "\n", "One issue we might experience when developing RL algorithms is that many aspects of the learning process are inherently random: initializing game states, changes in the environment, and the agent's actions. As such, it can be helpful to set a initial \"seed\" for the environment to ensure some level of reproducibility. Much like you might use `numpy.random.seed`, we can call the comparable function in gym, `seed`, with our defined environment to ensure the environment's random variables are initialized the same each time." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[1]" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "env = gym.make('CartPole-v0')\n", "env.seed(1)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "In Cartpole, a pole is attached by an un-actuated joint to a cart, which moves along a frictionless track. The pole starts upright, and the goal is to prevent it from falling over. The system is controlled by applying a force of +1 or -1 to the cart. A reward of +1 is provided for every timestep that the pole remains upright. The episode ends when the pole is more than 15 degrees from vertical, or the cart moves more than 2.4 units from the center of the track. A visual summary of the cartpole environment is depicted below:\n", "\n", "\n", "\n", "Given this setup for the environment and the objective of the game, we can think about: 1) what observations help define the environment's state; 2) what actions the agent can take. \n", "\n", "First, let's consider the observation space. In this Cartpole environment our observations are:\n", "\n", "1. Cart position\n", "2. Cart velocity\n", "3. Pole angle\n", "4. Pole rotation rate\n", "\n", "We can confirm the size of the space by querying the environment's observation space:\n" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Environment has observation space = Box(-3.4028234663852886e+38, 3.4028234663852886e+38, (4,), float32)\n" ] } ], "source": [ "n_observations = env.observation_space\n", "print(\"Environment has observation space =\", n_observations)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Second, we consider the action space. At every time step, the agent can move either right or left. Again we can confirm the size of the action space by querying the environment:" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Number of possible actions that the agent can choose from = 2\n" ] } ], "source": [ "n_actions = env.action_space.n\n", "print(\"Number of possible actions that the agent can choose from =\", n_actions)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Cartpole agent\n", "\n", "Now that we have instantiated the environment and understood the dimensionality of the observation and action spaces, we are ready to define our agent. In deep reinforcement learning, a deep neural network defines the agent. This network will take as input an observation of the environment and output the probability of taking each of the possible actions. Since Cartpole is defined by a low-dimensional observation space, a simple feed-forward neural network should work well for our agent. We will define this using the `Sequential` API.\n", "\n", "> Note: Of course, you can define your custom model, but in this environment, simple model can handle this." ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "# Defines a feed-forward neural network\n", "def create_cartpole_model():\n", " model = tf.keras.models.Sequential([\n", " # First Dense layer\n", " tf.keras.layers.Dense(units=32, activation='relu'),\n", " tf.keras.layers.Dense(units=n_actions, activation=None)\n", " ])\n", " return model\n", "\n", "cartpole_model = create_cartpole_model()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now that we have defined the core network architecture, we will define an *action function* that executes a forward pass through the network, given a set of observations, and samples from the output. This sampling from the output probabilities will be used to select the next action for the agent. We will also add support so that the `choose_action` function can handle either a single observation or a batch of observations.\n", "\n", "**Critically, this action function is totally general -- we will use this function for both Cartpole and Pong, and it is applicable to other RL tasks, as well!**" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "### Define the agent's action function ###\n", "\n", "# Function that takes observations as input, executes a forward pass through model, \n", "# and outputs a sampled action.\n", "# Arguments:\n", "# model: the network that defines our agent\n", "# observation: observation(s) which is/are fed as input to the model\n", "# single: flag as to whether we are handling a single observation or batch of\n", "# observations, provided as an np.array\n", "# Returns:\n", "# action: choice of agent action\n", "def choose_action(model, observation, single=True):\n", " # add batch dimension to the observation if only a single example was provided\n", " observation = np.expand_dims(observation, axis=0) if single else observation\n", "\n", " # feed the observations through the model to predict the log \n", " # probabilities of each possible action.\n", " logits = model.predict(observation)\n", " \n", " # Choose an action from the categorical distribution defined by the log \n", " # probabilities of each possible action.\n", " action = tf.random.categorical(logits, num_samples=1)\n", "\n", " action = action.numpy().flatten()\n", "\n", " return action[0] if single else action" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 3.2 Define the agent's memory\n", "\n", "Now that we have instantiated the environment and defined the agent network architecture and action function, we are ready to move on to the next step in our RL workflow:\n", "1. **Initialize our environment and our agent**: here we will describe the different observations and actions the agent can make in the environemnt.\n", "2. **Define our agent's memory**: this will enable the agent to remember its past actions, observations, and rewards.\n", "3. **Define the learning algorithm**: this will be used to reinforce the agent's good behaviors and discourage bad behaviors.\n", "\n", "In reinforcement learning, training occurs alongside the agent's acting in the environment; an *episode* refers to a sequence of actions that ends in some terminal state, such as the pole falling down or the cart crashing. The agent will need to remember all of its observations and actions, such that once an episode ends, it can learn to \"reinforce\" the good actions and punish the undesirable actions via training. Our first step is to define a simple `Memory` buffer that contains the agent's observations, actions, and received rewards from a given episode. We will also add support to combine a list of `Memory` objects into a single `Memory`. This will be very useful for batching, which will help you accelerate training later on in the lab.\n", "\n", "**Once again, note the modularity of this memory buffer -- it can and will be applied to other RL tasks as well!**" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "### Agent Memory ###\n", "\n", "class Memory:\n", " def __init__(self): \n", " self.clear()\n", "\n", " # Resets/restarts the memory buffer\n", " def clear(self): \n", " self.observations = []\n", " self.actions = []\n", " self.rewards = []\n", "\n", " # Add observations, actions, rewards to memory\n", " def add_to_memory(self, new_observation, new_action, new_reward): \n", " self.observations.append(new_observation)\n", "\n", " # update the list of actions with new action\n", " self.actions.append(new_action)\n", "\n", " # update the list of rewards with new reward\n", " self.rewards.append(new_reward)\n", "\n", "# Helper function to combine a list of Memory objects into a single Memory.\n", "# This will be very useful for batching.\n", "def aggregate_memories(memories):\n", " batch_memory = Memory()\n", " \n", " for memory in memories:\n", " for step in zip(memory.observations, memory.actions, memory.rewards):\n", " batch_memory.add_to_memory(*step)\n", " \n", " return batch_memory\n", "\n", "# Instantiate a single Memory buffer\n", "memory = Memory()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 3.3 Reward function\n", "\n", "We're almost ready to begin the learning algorithm for our agent! The next step is to compute the rewards of our agent as it acts in the environment. Since we (and the agent) is uncertain about if and when the game or task will end (i.e., when the pole will fall), it is useful to emphasize getting rewards **now** rather than later in the future -- this is the idea of discounting. This is a similar concept to discounting money in the case of interest. ecall from lecture, we use reward discount to give more preference at getting rewards now rather than later in the future. The idea of discounting rewards is similar to discounting money in the case of interest.\n", "\n", "To compute the expected cumulative reward, known as the **return**, at a given timestep in a learning episode, we sum the discounted rewards expected at that time step $t$, within a learning episode, and projecting into the future. We define the return (cumulative reward) at a time step $t$, $R_{t}$ as:\n", "\n", ">$R_{t}=\\sum_{k=0}^\\infty\\gamma^kr_{t+k}$\n", "\n", "where $0 < \\gamma < 1$ is the discount factor and $r_{t}$ is the reward at time step $t$, and the index $k$ increments projection into the future within a single learning episode. Intuitively, you can think of this function as depreciating any rewards received at later time steps, which will force the agent prioritize getting rewards now. Since we can't extend episodes to infinity, in practice the computation will be limited to the number of timesteps in an episode -- after that the reward is assumed to be zero.\n", "\n", "Take note of the form of this sum -- we'll have to be clever about how we implement this function. Specifically, we'll need to initialize an array of zeros, with length of the number of time steps, and fill it with the real discounted reward values as we loop through the rewards from the episode, which will have been saved in the agents memory. What we ultimately care about is which actions are better relative to other actions taken in that episode -- so, we'll normalize our computed rewards, using the mean and standard deviation of the rewards across the learning episode.\n" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "### Reward function ###\n", "\n", "# Helper function that normalizes an np.array x\n", "def normalize(x):\n", " x -= np.mean(x)\n", " x /= np.std(x)\n", " return x.astype(np.float32)\n", "\n", "# Compute normalized, discounted, cumulative rewards (i.e., return)\n", "# Arguments:\n", "# rewards: reward at timesteps in episode\n", "# gamma: discounting factor\n", "# Returns:\n", "# normalized discounted reward\n", "def discount_rewards(rewards, gamma=0.95): \n", " discounted_rewards = np.zeros_like(rewards)\n", " R = 0\n", " for t in reversed(range(0, len(rewards))):\n", " # update the total discounted reward\n", " R = R * gamma + rewards[t]\n", " discounted_rewards[t] = R\n", " \n", " return normalize(discounted_rewards)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 3.4 Learning algorithm\n", "\n", "Now we can start to define the learing algorithm which will be used to reinforce good behaviors of the agent and discourage bad behaviours. In this lab, we will focus on *policy gradient* methods which aim to **maximize** the likelihood of actions that result in large rewards. Equivalently, this means that we want to **minimize** the negative likelihood of these same actions. We achieve this by simply **scaling** the probabilities by their associated rewards -- effectively amplifying the likelihood of actions that resujlt in large rewards.\n", "\n", "Since the log function is monotonically increasing, this means that minimizing **negative likelihood** is equivalent to minimizing **negative log-likelihood**. Recall that we can easily compute the negative log-likelihood of a discrete action by evaluting its [softmax cross entropy](https://www.tensorflow.org/api_docs/python/tf/nn/sparse_softmax_cross_entropy_with_logits). Like in supervised learning, we can use stochastic gradient descent methods to achieve the desired minimization. \n", "\n", "Let's begin by defining the loss function." ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "### Loss function ###\n", "\n", "# Arguments:\n", "# logits: network's predictions for actions to take\n", "# actions: the actions the agent took in an episode\n", "# rewards: the rewards the agent received in an episode\n", "# Returns:\n", "# loss\n", "def compute_loss(logits, actions, rewards): \n", " # compute the negative log probabilities\n", " neg_logprob = tf.nn.sparse_softmax_cross_entropy_with_logits(\n", " logits=logits, labels=actions)\n", " \n", " # scale the negative log probability by the rewards\n", " loss = tf.reduce_mean(neg_logprob * rewards)\n", " \n", " return loss" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now let's use the loss function to define a training step of our learning algorithm:" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "### Training step (forward and backpropagation) ###\n", "\n", "def train_step(model, optimizer, observations, actions, discounted_rewards):\n", " with tf.GradientTape() as tape:\n", " # Forward propagate through the agent network\n", " logits = model(observations)\n", "\n", " # compute the loss\n", " loss = compute_loss(logits, actions, discounted_rewards)\n", "\n", " # run backpropagation to minimize the loss using the tape.gradient method\n", " grads = tape.gradient(loss, model.trainable_variables)\n", " optimizer.apply_gradients(zip(grads, model.trainable_variables))\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 3.5 Run cartpole!\n", "\n", "Having had no prior knowledge of the environment, the agent will begin to learn how to balance the pole on the cart based only on the feedback received from the environment! Having defined how our agent can move, how it takes in new observations, and how it updates its state, we'll see how it gradually learns a policy of actions to optimize balancing the pole as long as possible. To do this, we'll track how the rewards evolve as a function of training -- how should the rewards change as training progresses?" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "from IPython import display as ipythondisplay\n", "\n", "class LossHistory:\n", " def __init__(self, smoothing_factor=0.0):\n", " self.alpha = smoothing_factor\n", " self.loss = []\n", " def append(self, value):\n", " self.loss.append( self.alpha*self.loss[-1] + (1-self.alpha)*value if len(self.loss)>0 else value )\n", " def get(self):\n", " return self.loss\n", "\n", "class PeriodicPlotter:\n", " def __init__(self, sec, xlabel='', ylabel='', scale=None):\n", "\n", " self.xlabel = xlabel\n", " self.ylabel = ylabel\n", " self.sec = sec\n", " self.scale = scale\n", "\n", " self.tic = time.time()\n", "\n", " def plot(self, data):\n", " if time.time() - self.tic > self.sec:\n", " plt.cla()\n", "\n", " if self.scale is None:\n", " plt.plot(data)\n", " elif self.scale == 'semilogx':\n", " plt.semilogx(data)\n", " elif self.scale == 'semilogy':\n", " plt.semilogy(data)\n", " elif self.scale == 'loglog':\n", " plt.loglog(data)\n", " else:\n", " raise ValueError(\"unrecognized parameter scale {}\".format(self.scale))\n", "\n", " plt.xlabel(self.xlabel); plt.ylabel(self.ylabel)\n", " ipythondisplay.clear_output(wait=True)\n", " ipythondisplay.display(plt.gcf())\n", "\n", " self.tic = time.time()" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEGCAYAAACKB4k+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO3deXxcVfn48c+TZDLZ96RN031faaEFyg4tuwiIG+gPQVG+KiqiqIAi7uLyRcHlK5VVBFRAAVnKTlvW0oXu+55uSbNnksx6fn/cJZNksjaTNMnzfr36ysydO3fuDeE+c55zznPEGINSSikFkNDfJ6CUUurYoUFBKaWUS4OCUkoplwYFpZRSLg0KSimlXEn9fQJHo6CgwIwdO7a/T0MppQaUlStXHjHGFMZ6bUAHhbFjx7JixYr+Pg2llBpQRGRPe69p+kgppZRLg4JSSimXBgWllFIuDQpKKaVcGhSUUkq5NCgopZRyaVBQSinl0qCglFIDzN2vbmPlnsq4HFuDglJKDSBr9lXzu1e3smzbkbgcX4OCUkoNIE+uLCU9OZEvnjE+LsfXoKCUUgPIwZpGRuenk+GNT5WiuAUFEXlARMpEZH3Utt+IyGYRWSsi/xGRnKjXbhWR7SKyRUQuiNd5KaXUQHa41s+wLG/cjh/PlsJDwIWttr0CzDTGHAdsBW4FEJHpwJXADPs9fxaRxDiem1JKDUiHa5sYlpkSt+PHLSgYY5YCla22vWyMCdlP3wNG2o8vA/5hjPEbY3YB24GT4nVuSik1EIXCEY7UD9yWQme+ALxoPy4B9kW9VmpvU0opZavwBYgYKMoagC2FjojI94EQ8KizKcZupp33Xi8iK0RkRXl5ebxOUSmljjlH6v0AFGQMopaCiFwDXAJ81hjj3PhLgVFRu40EDsR6vzFmkTFmnjFmXmFhzIWDlFJqUGoMhAFI98avy7VPg4KIXAh8D7jUGNMQ9dKzwJUi4hWRccAkYHlfnptSSh3rGoNWUEjxxC8oxG05ThF5HDgbKBCRUuAOrNFGXuAVEQF4zxjzZWPMBhH5F7ARK610gzEmHK9zU0qpgagpGAEgJWkABgVjzFUxNt/fwf4/B34er/NRSqmBzmkppCbHL8mjM5qVUmqAaLKDgjeOLQUNCkopNUD43ZaCBgWllBry+qKjWYOCUkoNEM0dzdqnoJRSQ15jMIwnUUhK1KCglFKD0rrSGraX1Xdp36ZgOK7DUSGOQ1KVUkp17qN/fAuA3Xd+pMP9jDFU+QKkxLGTGTQoKKXUgHD/W7t4+sMDpMaxkxk0faSUUgPCyxsPA80jkOJFg4JSSh0DmuuDxlaYGb/KqNE0KCilVD+JDgS1jaEO9oTCOJbLjqZBQSml+okz7wDgiM/f4b7hSMctid6iQUEppfpJQ6C5dXCkruOg4PQl/Orjs+J6ThoUlFKqnzQEmjuND9U2dbhvYyDMhMJ0Pn3i6LiekwYFpZTqJ9FBYU9FQwd7Wi2FeBbCc2hQUEqpfuKLSh/trvB1uG9DIESaJ/5TyzQoKKVUP3HWXE5MEHYf6TgoNAbCcZ/NDBoUlFKq3/j8Vkth5ogsth2uxx9qf2JaYzBMWpxnM4MGBaWU6heRiOH6R1YCcNmcEur8IV7bVNbu/g2BMGnaUlBKqcEpulzF+TOGkelN4t0dFe3vr+kjpZQavPyh5olrxdmpjM5PY19V+yOQNH2klFKDmNN/cOcVs0hMEEbnpbG3MnZQMMZYQUFbCkopNTj5naU17W//o/PSKK1sJBKjnIU/FMEY8GpLQSmlBicnfeS111selZdGIBzhcF3bmc2tA0g8aVBQSql+0GR3NHs91m24ODsFgMO1bWsgOakmJ4DEkwYFpZTqB80tBevbf4FdGrs8RmG81q2KeNLlOJVSqh+0/vbvLKITHRSMMURM1L6aPlJKqcHJ6SdwWgr5GckAHKlvDgp3vbKVCbe9QG1TyN5X00dKKTUouSkhu0/Bm5RITpqnRUvhvmW7AKioD9j7DOCgICIPiEiZiKyP2pYnIq+IyDb7Z669XUTkHhHZLiJrReSEeJ2XUkodC2J1HhdmeFsEhbC9XGelvSqb06qIp3iGnYeAC1ttuwV4zRgzCXjNfg5wETDJ/nc98H9xPC+llOp3Ta3SR2B1Nkenj5w5C0ecloJnALcUjDFLgcpWmy8DHrYfPwxcHrX9b8byHpAjIsXxOjellOpvTkshJepGn5vuobox6D53WgqDIn3UjmHGmIMA9s8ie3sJsC9qv1J7Wxsicr2IrBCRFeXl5XE9WaWUipfWQ1IBslOTqW4IuM/tmEDFIEkfdYfE2NZ2rjdgjFlkjJlnjJlXWFgY59NSSqnec7i2ib32spvO6KPkqG//uWkeqhuCGNPy9jeYWwqHnbSQ/dMpHl4KjIrabyRwoI/PTSml4uoj97zFmb95A2MM/lAYT6KQmND8nTgnzUMoYvAFWi624/QzDOg+hXY8C1xjP74GeCZq++fsUUjzgRonzaSUUoOFc3NftbcKfyjSJh2Uk2rNVajyBVoUxqvwOS2FAZw+EpHHgXeBKSJSKiLXAXcC54nINuA8+znAC8BOYDvwV+Cr8TovpZTqLyPs+kYf7quhKRhukw7KSfMAUNMYpD4QcrdXOC2FgVzmwhhzVTsvLYyxrwFuiNe5KKXUsSAx0UoV+fwhu6XQOijYLYWGABLV0+o0GgZ0UFBKKdWS07ns84doCobblMLOtVsKVQ1BKn2BFq8lJyUgEmtMTu86VkYfKaXUoOfzWymhen+IxkCY1FYrqWWmeNz9dpT7SBA4e4o1yjI5sW9u1xoUlFIqjowxPLf2AHVNQXdUkc8foiHQdnlNJ0j4/CF2ltczMjeNS2ePAKxA0hc0KCilVBw9ubKUrz22mkVLd7rb6v1hGmOkj5wg0RAIU9cUIjfNw4Uzh/fp+WqfglJKxdETK0uB5hnMAPX+II2BMMOyvC329SQmkJyUQEMgjD8UJjkpgbTkJB774skcrGm7TGc8aEtBKaV60SsbDzP19hdpsIeUOj+jO459/jANwRBpyW2/l6cnJ9IQCBEIRdzZzqdOLODjc0f2wdlrUFBKqV7125e20BSMsMcuZ+E4XGt900/xJLBufw37KhvbdDQDpCUn4fOHCYTbTm7rCxoUlFKqFyXYZSvC9uSCgJ02ctZJKMpMcfdNjbG8Zro3qqXQRyOOomlQUEqpXuTcx//w+jZ8/pAbFJyWQkNUXaPWo48AUpOT8AXCLdJHfUmDglJK9aJEe4LZSxsO89uXt7hBoarBWiehKdgcFGKlj9KTE2mwZzxrUFBKqQEuIarqaXVDkEA40uL1RVfPdR+nxUgfpWlLQSmlBjZjDI12WigpKigEw5EWQ1GTExM4ZUI+V5xgrSHWep4CWH0KjXafQl/UOmpNg4JSSh2lR9/fy7QfLmZ/dSMJUfWJwhHTIijkpnsQEbeDuXUrAppbCv6wthSUUmpAWrOvGoBFS3aQlBjdUjBunwJAXro1Wc0JCo2tFtMBq0/B6aD26ugjpZQaeAozrZv9pkN1LVoK0Z3KAHnpVsG7CUUZAORntJzRDNaIJGeEUn+0FLTMhVJKHSWnNbC3ooGpxZnu9qoGaxZzenIivkCYmkZrBNKVJ46iKNPLgqlFbY6V5m2+LevkNaWUGoCcvoFDtU0tUkJl9oS1604fB8AF063idiLCwmnDYq6PkB41TFVbCkopNQBF9xvsKPe5j51ZzMU5qWz7+UUtRia1J7oeknY0K6XUAPLhvmrKaptaBIUj9nrK0ZITE/Akdm3ltHRvVEuhHzqataWglFI9dPmf3ibDm8RZUwrJTfO4s5Zb6843fm0pKKXUAFbvDxEMRRiWlUKGN/b37O7c3KNbCjp5TSmlBiCrzHUCo/LSACjOTuG5r5/uvt6doJDq0ZaCUkoNOKGo2chOnaIxdlA4bWIBM0uy3de7MwmtRZ+CBgWllBoYostXOEFhZG4q0Hwzd/qVvZ6e9Slo+kgppQaIFkEhHMGTmECJHRR8fmsJzvz0ZACSE7s+CS26pZDeTh9FPGlQUEqpHvCHmieprS2tITkxgRE5VlBw5ifk27WOujAS1ZUSNYt5QmFGL5xp92hQUEqpHvAHW1Y4TU5KoMQOCs5chfnj8wBaFMnrTPR6DB6dp6CUUgNDdPoIrKBQnG2tv+z0C/zgkulcOLOYqcOzunXsuWNyOW1iQe+caDf1S1AQkZuALwIGWAd8HigG/gHkAauAq40xgf44P6WU6kx0+gisTuH8DC+/+/Rs5o/PB6xv+qdMyO/2sZ/6yqm9co490edtExEpAb4BzDPGzAQSgSuBXwG/M8ZMAqqA6/r63JRSqqsCrVoKwbAB4GPHj6Q4O7U/TqlXdCkoiMgEEfHaj88WkW+ISM5RfG4SkCoiSUAacBBYADxpv/4wcPlRHF8ppeLKSR9dcby1tGZjsO2COQNRV1sKTwFhEZkI3A+MAx7ryQcaY/YDvwX2YgWDGmAlUG2MCdm7lQIlsd4vIteLyAoRWVFeXt6TU1BKqS4JRwxPrSwlHDFtXnPSR84CO/4hFhQi9g37Y8DvjTE3YfUBdJuI5AKXYQWWEUA6cFGMXdv+VwCMMYuMMfOMMfMKCwt7cgpKKdUl//xgH99+Yg2PvLu7zWvO6KP8DGsuQlOw7XrLA1FXg0JQRK4CrgGes7d5eviZ5wK7jDHlxpgg8G/gVCDHTicBjAQO9PD4SinVK2qbrKqnB2qa2rzmpI+yU61b4VBLH30eOAX4uTFml4iMA/7ew8/cC8wXkTSxiosvBDYCbwCfsPe5Bnimh8dXSg1RGw/U0hAIdb5jFzllJlqvtQzN6aPs1OR29xmIuhQUjDEbjTHfMMY8bj/fZYy5sycfaIx5H6tDeRXWcNQEYBHwPeBbIrIdyMfqu1BKqS6p94e4+J5l3PzEmg73W7z+EP/6YF/MfoLWUjzW7OLYQcFqKUwssmYdXz4nZjfogNPhPAURWUc7uX0AY8xxPflQY8wdwB2tNu8ETurJ8ZRS6nCtleJZuacq5us7yuupawrx5b+vBGBsQTonjcvr8JjGvvvF6i9w+hSGZXnZ+rOL8HRj1vKxrLPJa5fYP2+wfz5i//ws0BCXM1JKqR44ZOf9o6uMRrv1qXUs313pPq+IsWxma06KqPVEteht3qTEfilxHS8dBgVjzB4AETnNGHNa1Eu3iMjbwE/ieXJKKdVVB+2gkOqJXZG0pjHY4fNYnBZCzJZCKIIIg6aF4OhqmYt0ETndGPMWgIicijWUVCmljglO+igtuW1QWL+/xi1S56juICj8/tWtVDcEyUlrf2RRvT9EenIS0p0SqANAV9s8XwD+JCK7RWQX8Gd7m1JK9bt9lQ385qUtQNtCdQCX/OEtKnwtS6n964N9zLrjJRoDLW/4oXCEh97ZzRtbytwWQl1T2xFNVb4Auek9HZl/7Oq0pSAiCcBEY8xsEckCxBhTE/9TU0qprtld4XMfO3MLHMbEHiuz84j1nkO1TYwraE58rNpbTXVDkEjEuKOOymrbzlOoagiSl5Z81Od+rOm0pWCMiQBfsx/XakBQSh1rgvZ6ybNH5VDbKi0U3XJwMj1j8tPcbc4qaY7XNh0GoLYp5M55qPAFqGloedyqhgA5QzEo2F4RkZtFZJSI5Dn/4npmSinVRYGQ1RooSE+mtinUonVQH3XTv/Wiqez65cUUZHjdba07nN/YUuY+PlTb3A+x40h9i/0qfQHy0gdfUOhqR7PTf3BD1DYDjO/d01FKqe5zWgr5GcmEI4Z6f4jMFCvfH90SyEzxICLkpDb3BbRuWZRWNVKU6aWszs+hmkYSBCIGdpb7OGF0rrtflS9A7lBtKRhjxsX4pwFBKXVMcNY2GJVrpYUOVDf3AUR3Eqd7re/B4wub+xCiWwqRiKEhEGZUnnWcgzVNTCrKBKzObIc/FMYXCJM3CDuauzzjQkRmisinRORzzr94nphSSnWV01KYVmwte7m9rDnVs7a0uRvUSSvdetE0Xv3WWUDLjmmf3YcwKtdaJKeuKURGShI5aR4qo0YvVfms9wzGPoUupY9E5A7gbGA68AJWqeu3gL/F7cyUUqqLnKAwZXgmIs1BYXtZPbf9Z527X8QOCgkJwoTCdBITpEVLwee3Rhs5LQWwl9lMT6bC19y/8Nb2I+7nDTZdbSl8Aqua6SFjzOeB2YC347copVTvC0cM9y7ZQXVD8zf3gL0UZlaqh5G5qewot4JC9Lf7k8flcdHM5mVgRISslCRqG5vTS06n9PjCdDJTrO/MKZ5E8jO8HKlvPtbzaw8wJj+NeWOa+xgGi652NDcaYyIiErLnKpShncxKqX7wwrqD/PLFzZTV+TlY00hdU4jTJhYAkJyYQHFWKmV1Vp+CL6qM9t1XHu9WPXVkpXpatRSs/bNSPJw8Lo9XN5XhTUogxZPMlkN17n6Hav1MKsocdLOZoesthRX2msx/xVo6cxWwPG5npZRSWEXuom/GAGtLqwFI8STwwrpDLNt2hKDd0exJFHLTPW7OP7o1ke5tW/4iJy2Zqqh9nKCQ7k1yWxUjclLJT/eyo9zHOzustFFtY9BdXGew6eroo68aY6qNMX8BzgOusdNISinlOu3O1/nsfe/12vEW/u+bXPD7pS22bbaDxKGa5hx/MGwVp0tMEPLSk6lsCHCgurHFPukxqqcWpCdTEZUWctJHGd4kPj53JGt/dD7fv3gaoYgVdD7z1/d5b2cF+6sb3bpIg01XO5r/BiwDlhljNsf3lJRSA9X+6kb2Vzf22vF8dl2ihkDILYnt3MTf2l7u7ucPRfAkJiAi5KYlc6Tez6l3vu6+fsakAhIS2qZ68jOSWX+geXSSk25yhq5m2XMd5o/P5/Hl+wC4cpEV9IZ0SwF4CCgG/iAiO0TkKRG5MX6npZQaKBoDYe56eUtcl6PcfaR5joAzhPRw1GzjqoYAyYnW7Sw3LZnockdZKUk8ct3JMY9bkOGloj5AxF6Frd4efdQ61XTZnBK+evaEFtuGdEvBGPO6iCwBTgTOAb4MzADujuO5KaUGgIff3c09r28nI6Wr41a6b9cRH9NHWHMQYlUsLa/zu+sa5LYqPdHRXIL8DC+hiKG2KUhOWrLbp5DhbXstk4ZltHg+WFsKXU0fvYa1fsK7WGmkE40xZR2/Syk1FNTbN+n9Vc1po3DEkBgjXdNdmd4k6vwhd96BMYa6prbrIBypD+CxWwqtZxl39I2+ICPZfr+fnLRk6pqCJEjshXpGR81dAGvk0mDU1fTRWiAAzASOA2aKSGrczkopNWBU2qN3dpQ3l6/+2fMbe+XYQbuDd4Od9/cFwkRiVMI+Uu93g0J0sTuAEdnt36qcfZ05CLuO+BiVlxZzqGnr46a1s8LbQNfV0Uc3GWPOBD4GVAAPAtXxPDGl1MCwt8LK93+4r/mW8ODbu908fU/5Q2F3kZsNB2qB5uJ1w7Ja3qDL6/zuOsmzSrK5+8o5nDTOKuR83Kjsdj/DudE7ndfr99cyw05TtdY6DZU9SPsUuhQURORrIvJP4EPgcuABrFIXSqkhrrTKCgr1rdYlaL38ZXc5fQclOansr26kpjHobhuT13I14FDEuH0KIsJlc0rciqbHleS0+xn5Uemj2qYgeysbmDEidhDJiuozueOj05k6PHbwGOi62jOUCtwFrDTGtO3lUUoNWc6IndZKqxspykrp8XGdVsHkYRnsr250b9xgLZKzfHclw7NSOGSviuakjxw3LpzE7JHZnDYxv93PyE1LJkGgot5PmT2aaWRu7HRTdErp7ClFPb6uY11X00e/ATzA1QAiUigi4+J5YkqpgcEfCpPiab6VLLp6LmCtS3A03FZBvtUqqPIF3E5mZ+W0/IxkN23UOiikJidy0aziDktROJPdyusD7pDaWJ3MrcWaHT1YdDV9dAfwPeBWe5MH+Hu8TkopNXD4gxHG5jenc2aNtNIv+48yKDjrLk+1K5FW+gJu8Tqnimm6N8kdGuoEh+7KT/dSUe93g0Lr+kixZHoHZ38CdH300ceASwEfgDHmADD4asYqpbolHDEEwpEWC98XZ6dSlOll2+G6Dt7ZuTe3lJOXnuwWu7v+kZWsseseTSi05gykJye6q6glJ/YsKBRkWjOgnU7trgSF6JbRYNPVKwsYa3UKAyAi6Z3sr5QaAvwh69v1WDsoOMHhuJHZ7g28p1btrWL++LwWQ0EffHs3iQnifk50S8HpaO6u/HQvFb5AVEuh89viYKyO6uhqR/O/ROReIEdEvoS1ZvN98TstpdRA4Hy7Lsr0ctenZrvf6meV5PDa5jLqmoLuWsnd1RgIk5XiITW55Tf33DQP6d4kUjwJZHiT3MlprfsUusopddEU6jx99MbNZ7P7iK/d1weDrnY0/xZ4EngKmAL80BhzT08/VERyRORJEdksIptE5BQRyRORV0Rkm/1z8K1eodQgE52Hv+KEkQyzRxsdNzIbY5rnF/REIByJ2U+Qn261HL6+YBKXzh7hziz29LRPISOZen+I6garEzslqf2gMK4gnXOmDt6RR9CNNZqNMa8YY75jjLkZeF1EPnsUn3s3sNgYMxVrFbdNwC3Aa8aYScBr9nOl1DGsvZSL09m8NiqFFIkYGu2qp9vL6vjhM+sJ2ctoxhIIRdx+gujUUJ5d2+iGcyZy6sQC95t9RzfzjhTa6Smnuutg7i/oig6vXkSyRORWEfmjiJwvlq8BO4FP9eQD7ZXbzgTuBzDGBIwx1cBlwMP2bg9jTZJTSh3D3M7ZVjfkggwvJTmprC1tLkt95+LNTPvhYgKhCN//z3r+9u4e1u6voT2BUHNLYe0dF1CSY80fyMtoObPYKbH08bklPboGZwLb82sPWteSPHiHm3ZFZyHxEax00Trgi8DLwCeBy4wxl/XwM8cD5cCDIrJaRO6zO66HGWMOAtg/Y7bRROR6EVkhIivKy8tj7aKU6iMd5eFnlWSzLuqm/+h7ewCrzLXTT9DeCKVIxBCKGDcopCYnMjzbSk3lt6qC+q3zpvDgtSdy6oSCHl2D05G9t9Kamd3TFsdg0VlQGG+MudYYcy9wFTAPuMQY8+FRfGYScALwf8aY47GGuXY5VWSMWWSMmWeMmVdYWHgUp6GUOlpO+sgbI+Vy3Khs9lQ0sL2snqZg2F3kpqI+QJodFNa101II2Gml6D4Fp2zGxKKWJazz0pOPKs+f36rl0dNRTINFZ0HBrVFrjAkDu4wxRzf4GEqBUmPM+/bzJ7GCxGERKQawf2ppbqWOcf5Q+2P7nZpD5961hO89tdYtpV3h87tVSffYxfTKapsYe8vzvLG5rMVxo+ceOAHIWTu5t7SufjqYh5t2RWdBYbaI1Nr/6oDjnMci0qNhBcaYQ8A+EZlib1oIbASeBa6xt10DPNOT4yul+o7f6WiOkXI5fnRzIbq3t1eQaN9s395ewfJdlQBUNwS5d8kOTv7lawD83U4xBeyg4I1qKdx/zYn87ydnU5jZ8iZ+tFI8icwe2X4l1aGmw6BgjEk0xmTZ/zKNMUlRj4+mRODXgUdFZC0wB/gFcCdwnohsA86znyuljmHNs4Db3krSvUk8cO08AMYXpLsthb8s2eHuc6C6kV++uNldPtNJF8VKH80syebjc0f2/kUAj18/Py7HHYjit35eB+w+iXkxXlrY1+eilOq5zuoFLZg6jI8cV+yO7Gmtwhdo8dxpGTgtkJ7WM+qurhTBGyqG9oBcpdRR6UoRuYL0tmskTyzK4MaFk9psP1zr55v/WE2VPZHM20cjgYZ6P0K0fmkpKKUGvsZAmB/911p2s6MJX7kxgsK04iyKstr2Dby7swKwi6zR8yJ3quc0KCilesRZ3AY6Tr846yIA7qI4WSlJ5KW1DRaO3faopL5KHwFcOnsESQnaYtAwrJTqkRp7ZbQHrp3XYfrlrMnWfKIXbzyD+eOtdZOzUj2kea3vpOML0nn7lgWcbK+pDLDXXkuhL4PCPVcdz12fntNnn3es0qCglOoRJyg4pavbc+bkQnb84mKmFWe5gSArxcOckTmcPrGAv14zj5KcVKaPaB7Q6PQp9GVQUBZNHymleqS6wRo51FlQANzhqNEF7rLTPPz9iye7+8RKJ2mfQt/T37hSqtt2H/Fx4z+sajdZXQgKDqeERChi2ryWE6ND2qsthT6nv3GlBpgfPbuBlzcc6tdz+MlzG93HXWkpOJyhq6ZtTHBbCmPy09xtmj7qe/obV2oAqfIFeOid3Vz/yMpeOd4j7+3p0VrKh6NGHnVnLsEXzxjPJ+aO5LPzR7d5LddeQS3D25zV1qDQ9/Q3rtQAsnpfVa8dKxSOcPvT63ls+d5uv9dZkKa7slM9/PaTs8mKsURnjt1SSI8OCtqn0Of0N67UALJmn1VqOvrbdE/V2vMHyur83X6vzx/qfKduclZU05ZC/9LfuFIDSJU94qfeH3JLTPSUM6S0LCoV1BWBUIRg2HDGpAIevPbEozqHaDl2+ijdm8TpE60Fc/qqzIVqpkNSlRpA6qO+oe+vbmRCYUYHe3fMGVLa3ZaCs87yWZMLe3UR+xRPIqmeRDK8ifzq43PZWe7TlkI/0KCg1ADS4G9uHZRW9TwoRCKGB9/eDVidxsaYLheFawhagSm9F1JYrf3gkmnMHJFNWnISM0t0jYP+oEFBqQHEFwgxLMvL4Vo/pVUNPT7O8+sO8uyaA4C1JkJtU6jLQ0sb7JZCWhwWuP/syWN6/Ziqe7RtptQA4vOHGF+QQVKCUFrVsxFA0Hb0UHld1/sVnPSRrkEwOGlQUGoAaQiEyUxJYkRO6lEFhV3lvhbPD9d23q8Qjhi2l9W5LYV4pI9U/9OgoNQAsLO8nm/960M2H6ojw5vEyNzUo0ofbS1rOWHtcBdGIN27dAfn3rWUlXusuRKpcUgfqf6nQUGpAeAXL2zi36v2A5DmTWR4VgplXfh2357axiBpyYl88fRxQNdGIG06aAWS9+yFcOLRp6D6nwYFpY5xxhhe31zmPk/3JlGY5aW8zo+JVUSoC/yhCBfOHM4PLplOhjepSy2FggxrctmSreXWeSRr+mgw0qCg1DGuMRgmYppLPiQnJjAsM4VAOEK1ve5AdzUFIzHd/EAAACAASURBVG5xuqIsb5daHQmthqxq+mhw0qCg1DHOZ89NKMlNBayZyM76xoe7MWoomj8UdstSl+Sksrey8/6JuqaWAUjTR4OTBgWljnENAWuyWEmOFRSqGoIUZaYA9LhfwR+KuCUkphdnseVQHYFQpN39X9pwiH+tKGVEdoq7LUVLUAxKGhSUOsY5LYWF06ySElccX0JRpt1S6GbdIrBmMwdCEVI81v/+s0ZmEwhHuO+tne2+53/sUt2FWSl854IpDMvykqCL3A9K2lOk1DHOaSlMLMpg950fASAYjuBNSmDzoe6vhRAIWy0Cp6Uwd0wuAL9evIUvnj6+w3pDjYEQN5wzkRvOmdjtz1UDg7YUlOpHXRk95HPLSjR/h/MkJnDcyGxW7+3a+gqRqOUv/UEnKFj/+xdnp/Kzy2cCsKfC1+a9ZVH9Fj3t2FYDhwYFpfrJnS9u5uRfvMbBmkYq6tvvG2jwOwXoWubwTxidy7r9NWw6WNvh59Q2BRl/2ws89PYuAJpCVpBJiSpTMXtkDgA7yuvdbX98fRtf+tsKVu5uDjxOuW01eGlQUKqP3fDYKr775Br+smQHZXV+zvz1G8z92auEwm07eneU1/OVR1cBbecFXHvaWNKSk/jLkh0dfl5FvVUi+87Fm4G2LQWA8YXpAGwvaw4Kv315K69sPMzSbUfc4bBHU6pbDQz91qcgIonACmC/MeYSERkH/APIA1YBVxtjAv11fkrFy/NrD7Z4HgxbqZ0d5T6mDM/kuoc+YNXeKj514iiGZzWP9mk9BLQ4O5UZI7LY18lw0np7hbUmOxj47ZaC19McFNK9SYzITmFHedv00dOr9zN/Qj5fPms8k4oyu3qZaoDqz5bCjcCmqOe/An5njJkEVAHX9ctZKdVP1pRWEwpHeG1zGVUNQe5dstMtPgexC9CV5KR2ul5ynb855XOk3u8Gh9ZDSicUZbjpo+hWS2MwTElOKqdOKKDQHvWkBq9+CQoiMhL4CHCf/VyABcCT9i4PA5f3x7kpFU/hSMuO5cnDmtMxGw/Utql8uvVw8+gib4xRQSW5qZTV+TucY+C0FAD++cG+mC0FsFJD28vquW/ZTtYfaNlPocFg6OivlsLvge8Czl9yPlBtjHH+ekuBklhvFJHrRWSFiKwoLy+P/5kq1YtqW3XUnjQuz328u8LHriNW+sZZvP79nZXu67FWRivJScUYWL6rss1rjjo7KBRkeHlx/UH8dgBJ8bRtKTQEwvzs+U384Ol1LV4rtOseqcGvz4OCiFwClBljVkZvjrFrzLF6xphFxph5xph5hYWFcTlHpeKlulVQmDI8y328t7LBTd/86uPHAXDInpx29pTYf+uTh1k5/pufWONuW76rkipfc3ecs67zx44fwfr9teypsPogWrc8zp1WxEdmFTOzJIv1+1u2FAoytKUwVPRHS+E04FIR2Y3VsbwAq+WQIyJO0nQkcKAfzk0NQjvK6zsc8tmXWq+BcI59s59YlEFpZSM7yuvJTvUwfURzsPjSGeN46PMnxTze7FE5XHPKGA7VNlHTGCQSMXzq3ne5+J5l7j5OzaKPHT8SgP+sLgWaJ685irNT+dNnT+CTc0e525wWi6aPho4+DwrGmFuNMSONMWOBK4HXjTGfBd4APmHvdg3wTF+fmxqcvvjwCn78341xOXZZbRP/88gKd+GZjlTU+7n6/uUAnDttGKdOyGdkbhq77/wInz9tLIFwhHd3VDCuIJ3iqBpDTp2j9pwz1Sp/selgLT579vPBmiaMscpZ/HPFPpITE5g+IosZI7L4wJ53kOKJ/b9/dEDKSrGCQm66po+GimNpnsL3gG+JyHasPob7+/l81CBRXufv0k27J/735a28tOEwD72zu9N9n/6wufF7y0VTeexL893n4wqseQK7KxoYX5jeIt/vVERtz4wR2QCsK61x+w8A9lU28tA7u9hX2eiWtvjxpTPc11u3FBxThjcPO/35FbMYX5DuFuNTg1+/1j4yxrwJvGk/3gnEbiMr1UPhiKHeH6LeH6K8zt/raRAn579hf02n+y7fVeE+zk3ztHhtyrDmG/F4O0AkJyYQCEco7CSfX5jpZVJRBq9vLuPMyc19Dyv3VnKkvuVUn3lj87jkuGKeW3uQjJTY//tnpVjndvK4PM6ZUsQ5U4o6vTY1eGhBPDWoRQ/HXLOvmnOnD+vV49fa+fqdR3xU+QJkpXpIbKd6aJUvyMySLL513mTyW93oo5/PHWONSMpISaLSF+hSIDtv+jDuXbqTgzXNQ1rvf2tXmw5jgHuuPJ5bL57m9hfEsv7HF+BJ1CqoQ9GxlD5SqtfVRi0M8+G+6t4/vr3WMcC2snrO/PUbXPrHt3h3R0WbfSsbAozKTWPB1NiBKdP+5n6yPUz1a3Yl0uHZHfcpgFXpNBwxfLDbGpqan57cIiCs/MG57uOEBOk0HZThTWo3vaQGNw0KalCLLuC2pjQOQaEpxAmjrdLTS7eWs7+6kbWlNdz54qY2+1b5Ah122L5805ks+c7Z7joFnz9tLLt+eTGZKZ523+OYVmx1Dn+wy+o7+e6FU9zXPjF3ZJuWiVLt0fSRGtScjtdxBel8uK+aSMT02uIwkYihtjHItOJMPthdycsbDwEwdXgmO8t9GGPcCWeRiKGqIUB+B0GhOLvlt/dYk9Xaf28KWSlJLLdbCudMKeKdWxbgTUqIWR5DqfZoS0ENak766MxJBdQ1hdgVY72Anrh3yQ7G3/YC/lCEnLRkphZnsfWwNfHs/OnDqPOHWnTy1jYFiRjITYvP0E4RcVsLAJkpHkbkpJKf4W0zc1mpjmhQUIOa01I4Y5I1Kmf13pYppNV7q2KWrO7M39/f4z7OSkniC6eNBax+gbljrT6BDQeaRyRV2jOM8+I43t8JCokJ0u4cBKU6o385qs9tL6vna4+toikY7nzno+TUGjp+dA756cks29ZcL2tvRQMf+/M7TPz+i6zY3X7toFjGFzQXsstK9XDp7BEsunouj39pPjNHZJGWnMhPntuIMYb/rjnAgv9dAsR3EthUe35BOGK6lXpSKpoGBdXn7nxxE8+tPcjSrfEvaOi0FLJTPSyYWsQbm8vciqIVvubSF/cubX/R+liSovol0pKTEBHOnzGcmSXZ5Gd4+erZE9hZ7qO6IchN//zQ3XdG1Gzh3naCvdby1xfo+smq5zQoqD6XlWqNpjlcF/96RLVN1pDRpMQELpk9gtqmEC+utxa5iZ79uzNqGcquqGyw0kFJCdKi/LVjql3o7o0tZYTsctnfvXBKXAvLTR6WycofnMu3z5/S+c5KtUODgupzzrKSO8q6dyPuibqmoDtD94yJBYzOS+M/q/cDzdVDL5gxjF1HfN1KZ1U3BLnkuGK2/fwixuSnt3l9bEEaAL9evIUEgfdvW8hXz47/N3gdeqqOlgYF1eec8tHxmEzWWm1jyJ0UlpAgnDwujzX7qjHGuLOdT59YQMTAn9/Y3uXjVjUEyE1Lbjd3PzLXCgqHapu4aGYxw7I6n4Cm1LFAg4Lqc06t/w/3VXe6vvDRqvMH3aAAVqnpqoYgpVWN1NkthY/OHsGls0fwlyU7OWzXMmotHDH8/tWtVPoChCOGmsZgm/pF0VI8iVxzyhi+fNYE7vr07N69KKXiSIOC6nNVDQE3D3/7M+s7XWP4aNQ2htw+DLBGIQG8t7PCbSlkeJO4+fwpBCMR/rF8HwDGGBoCzX0OH+yu5PevbuO2f6+jpjGIMZDTyZyDH182k1sumqrlItSAokGhD/zx9W1cff/7fPMfq93x6oNRlS/A6b96nf+u6Xh9pCpfgFkl1s35zS3lfO2xVXE7p7qmYIsyEdOLsyjOTuHljYep9wdJ9Vid0KPz0zh5XB7/WV2KPxTmvmW7mP7Dl9z/Xo12f8PiDYf45wdW4OispLVSA5EGhT7w25e3smzbEZ7+8AB3v7q1v0+nUzvK61tU2+yqbWX1lFY18vXHVxOJxFxNFWMMlQ0B8tI9fMMeOrnxQC3GxN7/aNU2hdyFYsCa+btwWhHvbD9CTWPL1NK1p45jd0UDd728lf+utQLbH1/fzq8Xb+ZI1EipXy3eDNBiIRylBgsNCn1s8YZD7d4wjxUL/3cJp/zy9S7vb4zhR89u4IV1B91tO9oZ4rm2tIamYITJwzL55rmTueOj0/GHIuyu6P2+BWNMm5YCwJxRufgCYf61orRF+egLZw7nrMmFvLrpsLva2QNv7+LPb+7glY2H2xx/eLYuPKMGHw0KceZMlAK4/ZLpHK71s+VwXT+eUcd6ErB2lNfz0Du7W6w+trqdkUUvbThEUoJw3vRhJCQIJ9olITYeaC7zbIzhvmU7eTXGjbg7moIRgmFDVmrLgnDHjcx2H+880rIW0vzx+ewo91HV0DLN97J9Lj+9fCYACQJFum6xGoQ0KMRZdaN1c/npZTNYaK+l2xdDMXuqLCpNcqSLi92/uaV5ZnJyYgLZqZ6Y6wkAvL+rkuNGZrudtBOLMkhMEDYfag4Kz645wM+e38RXH1vF3a9uY00Pf1/OjOXs1JYthQmFGRRkWJ8/s6TlDOOTx1tBylm+87sXTnE7p0fnpfHJuSNJ8SRQmOnFk6j/+6jBR/+q46y6wRqTn5OWzJj8NHLTPHy4t5rSqgZ+8PQ6GgPxr//THaVVzWmcdaXNBd22Ha7jU3951x1OGi06tVKY6eWCGcN4ZeNhdzLY3a9u4/MPLqfeH2Llnirm2a0DsIZujs1P4w+vb2eX/a196dYjgNXK+t2rW7nsT29z23/Wdfta1ttLZEZXDwWrYNz7t53LmjvO5/GodZIBZpVku4vmXDp7BF89eyLHj7LKR+SmeUjxJLJw6jCmDI9fuQql+tOQDAoHaxp55L09BMMR1uyr5kAch0Q6QcGZ6DR3TB7v7qzg2gc/4O/v7WVJH9T/6Y59UUHBWR8A4BcvbGL57koWbzjUYv/yOj/Ld1cywu50TU1O5PI5JdT7Qzy/1upj+N2rW3ljSzln/+YNADdl5Lh8TgkAT6zYR70/xHNrD7TpxH3s/b088u7ubo3eWr23muTEhJj1hhIThOxUT5v+Bk9iAnPtGkJO2mnhtCLSkxP50pnjAbjr07NZdPXcLp+HUgPJkFx9Y/Xeam5/ej3Ti7O4/m8rmDI8k4c+fxKlVQ2ML2xbx6annA5YgBx7otNZkwt4dVPzN+v3dlZw4czhvfaZ3T2/XUd87K9u5C9LdnDrRdM4VGOlXM6dNozn1h7klgunkZ3mYV+VFThfXH+IK08cBcAj7+0hQQRj4K5Pz2H9/hqmDM/klAn5TCzK4K/LdnLxrGL388IRwz1XHc+CqS0Xgv/6wkm8trmMD3ZX8p0n1uAPRbj21LGMyU+nwufniuNHMu2Hi7n9mQ3c/swGZozIYs6oHERgZ7mPx1p923dsPlTHpGEZ3Z4ncPX8MSzbdsRND502sYANP7nQfV3nHajBbEgGBWcN3OfWHqDCF+CdHRVc9/AHLNt2hLdvWdDp+rVdVekLsPGglStvDgpFwAZ3n7e3H+mVz3L8+c3tzBmZw/zx+Z2uMPbkylK+8+RahmelcKi2ibte2cqkogySkxK4+YLJXHT3Mn72/EbW7a9hu12naOnWcp5atZ+Tx+Xxw2ear2NCYQbzx+e7z29cOImvP76a379mDcE9d9owbrt4artB9+Txedy7xKpUesr4fK49bWyLm+81p4zh4XetNQw2HKhlQ1TH9IYDNcwYkU1rdU3BHq1fcP6M4Sy6em6LNJdSQ8WQTB/lZ3iZMiyTB9/e7W5bts26OXc28ao7aqOqcDrVMUfnpzG+wCqg9pmTR7OtrJ6yutilFbrLGMOvF2/hM/e9zzUPLu9w35qGoDve/pBd2mHV3ioqfQFy0zxMHZ7FSWPzeGJlKZsPWaOlHrz2RHLSPCzdWt7ippyWnOh23Do+MquYzJQk/vaOdSO/5aL2AwLAF04b597Av3PhlDbfxn9wyXQeuHYeRZle/vU/pxAd7y7749v8KUbdonp/qMWQ0+44f8bwuC6Io9SxakgGBYDLjy9xH586IR+nrtn9b+2irinYzru6x1ngZdHVc1ssiXjhzOEMy/K6aZiXNhzd0EuHU8sHrCDndPQaY9pUAP3D69taLBcJVv/Hyj1V7pKRHzmuuMXrM0qyOH1iAc+uOcCX/77S3d4QCLcpDJeQIJwwOpfGYJgEgVF5Hbe+hmWlsPS75/D2LQs4YXRum9c9iQksmDqM5d8/l5PG5XHjwsn88TPH8++vnkooYnjs/b1t3lPf1POgoNRQNWSDwtWnjOHMyYVcclwxf/vCSaz/0QU89ZVTKK/z8+K6Q50foAucev2tV9u66bzJvHzTWcwqyWbumFx+98pW/KGjH4VU2eomP/X2xQRCER5+ZzdTb19MRdQQ021l9aR4Evjr5+YBuJ2rO4/43FTXVSeN5sUbz+DuK+dwwYxhFGWmuKk3aL7RO0M2WzvJ3jcpIaFLefgMb1KXU3c3njuJS44bwQmjc/nJZTPYX93I3lYT4Or8ITJSNCgo1R1DNihkeJP42xdO4o+fOYGkxATSvUmcMDqXokwvS7b1zoggZ9H4rBgjXLJTPYgIN5wzgUpfoN1x/d3hLPxyz1XHu9tW7qni6Q+tlFj0fILdFT4WThvGgqlFfOu8yfzpMye4N3ynpeBJTGBacRaXzSnh3qut4HFiVFB449tns/y2hTx07Ukxz+eCGcMACPRgDeTuONNef/mF9c0zqo0x1PtDZGpLQaluGbJBIRYR4czJhTy/9iCPvr/nqMtROOmj1jNqo506oYC05EQee38vt/57LYvX97yV4rQUxuSlse5H55OUIFz11/fcyXIPvL2LKl+AdaU1lFY1Mi4/ncQE4RsLJzE8O4XTJhYA1nDN9kwustYBLsz0kpSYQFFWCtntlJCeWJTJ/PF5/OAj03p8TV0xtiCdE8fm8tj7ewnaAaghEMYYtKWgVDfp/zGtnDW5kCdXlvL9/6ynJCeVs6cUdf6mdrTXUoiW4knk6vlj3DWCH1++j3duWcCILqZRyuqa2FXuY/2BWjcFlZeeTGaKh6+cPYE/vN7cAbvlUB3H//QV9/m4gpYrhs0ZZaWBOlrjICFBWPzNM8jrpGy04x/Xn9Kl/Y7Wl8+awHUPr+DBt3dx7anj3FXVMrzt/+6VUm31eUtBREaJyBsisklENojIjfb2PBF5RUS22T/b9jb2gdPtb8sAf3/P6ryMRAw3P7GGS/6wrFud0LWNIRITxJ0h256bzpvMDy+Zzo8+Oh2w6gN11ZWL3uPTi97jp89t5NeLtwC4o2a+ff4UnviydVM+Z0oh9197ovu+m8+fzEWzWs6PmFViDes8q5NAOHV4FkXH2EpiC6YWcfrEAn7xwmY+e997bitNWwpKdU9/pI9CwLeNMdOA+cANIjIduAV4zRgzCXjNft7nctOT+fdXT+Vzp4zh1U2H+dpjq9hyuI4nV5ayfn8tt/x7HXsrGrqUWqptCpKVktTuko2OFE8iXzh9HNeeNo6pwzO7NSx2Z7mvzbboIDRvTC4/u3wmv7hiFmdNLuSJL5/CmjvO52sLJpGW3PKGmZuezOrbz+PGhZO6/PnHChHh5x+bSXF2Ch/sruLie5YBaJ+CUt3U50HBGHPQGLPKflwHbAJKgMuAh+3dHgYu7+tzc5wwOpevnD0BgOfWHuT1zWWAtWzj82sPcuZv3uDu17a1ed9Vi95j0dId7vOqhrZlmztzxQklrNpbzc+f38hXH13Js50ECGd+wPTiLN6/bSGPf2l+iyAkIvy/+WMotss8nzg2r02BuGi56ckd9ikcy8bkp/POLQu4/szxBMNW0NaWglLd068dzSIyFjgeeB8YZow5CFbgAHqezO8FxdmpLPvuOQD85iUrLXPnFbPcmjf3v7WL6oYA+yob+Ml/N3KgupF3d1bwixc20xQM0xQMs3RrObNHxR6u2R6nDtBfl+3ihXWH+OY/VrNkaznGGN7efqTNfINQxFCcncID157IsKwUTpmQH+uwQ4aIcPP5U0hOsv60CzO0vLVS3dFvX6NEJAN4CvimMaa2sxRL1PuuB64HGD16dPxOEBiVl8Z1p4/j/rd2AZDuTeL8GcN56ZtncuHdS7lv2S52lNfz4vpDPPr+Hvd9S7eWc6Q+QE1jkE/PG9WtzyzKSmFMfhp7Khp479aFXPvgcq55oHl2cm6ah7994WRmjcymKRimuiHIF88fx3BdBcyVnJTA6tvPY19VA2NbdaYrpTrWL0FBRDxYAeFRY8y/7c2HRaTYGHNQRIqBsljvNcYsAhYBzJs3L+5LmN1+yXSC4UiLlMqU4ZlcPKuYP0aVVvBHLaZz79Kd7CyvZ+6YXE6b2P1v7k/8zymU1/sZbrcAPvfAcrf2UFVDkCv+7222/PQiyu21D461Tt9jQbo3iala3lqpbuvzoCBWk+B+YJMx5q6ol54FrgHutH8+09fn1p6fXDazzbabzp3sloYGa7GWk8flU9sY5ImVpQD89pOzO+1kjqUoK8W90Y/ISeXVb51FXVOQNftq+N5Ta9lf3chdr2x1ZxL3VgE/pZSSeC2Y3u4HipwOLAPWAc7X69uw+hX+BYwG9gKfNMZUdnSsefPmmRUrVsTxbDsWiRhuf2Y9l80pcUs67Cyv59y7lvDjS2dw9Slje/0zjTF876m1/GtFKYWZXvzBMB/84Fwt56yU6jIRWWmMmRfztb4OCr2pv4NCe3z+EOlxHArZGAhzzm/f5FBtE58/bSx3fHRG3D5LKTX4dBQUdLxeHMQzIIC1utnDXziJpz/cPyDnFCiljl0aFAaoKcMz+d6FU/v7NJRSg4wWxFNKKeXSoKCUUsqlQUEppZRLg4JSSimXBgWllFIuDQpKKaVcGhSUUkq5NCgopZRyDegyFyJSDuzpdMfYCoAjvXg6A8VQvG695qFBr7nrxhhjCmO9MKCDwtEQkRXt1f4YzIbides1Dw16zb1D00dKKaVcGhSUUkq5hnJQWNTfJ9BPhuJ16zUPDXrNvWDI9ikopZRqayi3FJRSSrWiQUEppZRrSAYFEblQRLaIyHYRuaW/z6e3iMgDIlImIuujtuWJyCsiss3+mWtvFxG5x/4drBWRE/rvzHtOREaJyBsisklENojIjfb2QXvdIpIiIstFZI19zT+2t48Tkffta/6niCTb27328+3262P78/yPhogkishqEXnOfj6or1lEdovIOhH5UERW2Nvi+rc95IKCiCQCfwIuAqYDV4nI9P49q17zEHBhq223AK8ZYyYBr9nPwbr+Sfa/64H/66Nz7G0h4NvGmGnAfOAG+7/nYL5uP7DAGDMbmANcKCLzgV8Bv7OvuQq4zt7/OqDKGDMR+J2930B1I7Ap6vlQuOZzjDFzouYjxPdv2xgzpP4BpwAvRT2/Fbi1v8+rF69vLLA+6vkWoNh+XAxssR/fC1wVa7+B/A94BjhvqFw3kAasAk7GmtmaZG93/86Bl4BT7MdJ9n7S3+feg2sdad8EFwDPATIErnk3UNBqW1z/todcSwEoAfZFPS+1tw1Ww4wxBwHsn0X29kH3e7BTBMcD7zPIr9tOo3wIlAGvADuAamNMyN4l+rrca7ZfrwHy+/aMe8Xvge8CEft5PoP/mg3wsoisFJHr7W1x/dtOOoqTHagkxrahOC53UP0eRCQDeAr4pjGmViTW5Vm7xtg24K7bGBMG5ohIDvAfYFqs3eyfA/6aReQSoMwYs1JEznY2x9h10Fyz7TRjzAERKQJeEZHNHezbK9c8FFsKpcCoqOcjgQP9dC594bCIFAPYP8vs7YPm9yAiHqyA8Kgx5t/25kF/3QDGmGrgTaz+lBwRcb7oRV+Xe83269lAZd+e6VE7DbhURHYD/8BKIf2ewX3NGGMO2D/LsIL/ScT5b3soBoUPgEn2qIVk4Erg2X4+p3h6FrjGfnwNVs7d2f45e8TCfKDGaZIOJGI1Ce4HNhlj7op6adBet4gU2i0ERCQVOBer8/UN4BP2bq2v2fldfAJ43dhJ54HCGHOrMWakMWYs1v+zrxtjPssgvmYRSReRTOcxcD6wnnj/bfd3R0o/dd5cDGzFysN+v7/Ppxev63HgIBDE+tZwHVYe9TVgm/0zz95XsEZh7QDWAfP6+/x7eM2nYzWR1wIf2v8uHszXDRwHrLaveT3wQ3v7eGA5sB14AvDa21Ps59vt18f39zUc5fWfDTw32K/ZvrY19r8Nzr0q3n/bWuZCKaWUayimj5RSSrVDg4JSSimXBgWllFIuDQpKKaVcGhSUUkq5NCioIU1E6u2fY0XkM7187NtaPX+nN4+vVDxoUFDKMhboVlCwK+52pEVQMMac2s1zUqrPaVBQynIncIZdt/4mu+Dcb0TkA7s2/f8AiMjZYq3f8BjWBCFE5Gm7YNkGp2iZiNwJpNrHe9Te5rRKxD72ertW/qejjv2miDwpIptF5FF7xjYicqeIbLTP5bd9/ttRQ8ZQLIinVCy3ADcbYy4BsG/uNcaYE0XEC7wtIi/b+54EzDTG7LKff8EYU2mXnPhARJ4yxtwiIl8zxsyJ8VlXYK2DMBsosN+z1H7teGAGVs2at4HTRGQj8DFgqjHGOCUulIoHbSkoFdv5WHVkPsQqxZ2PtXgJwPKogADwDRFZA7yHVZBsEh07HXjcGBM2xhwGlgAnRh271BgTwSrZMRaoBZqA+0TkCqDhqK9OqXZoUFAqNgG+bqwVr+YYY8YZY5yWgs/dySrjfC7Wgi6zsWoSpXTh2O3xRz0OYy0gE8JqnTwFXA4s7taVKNUNGhSUstQBmVHPXwK+YpflRkQm25UqW8vGWvaxQUSmYpWwdgSd97eyFPi03W9RCJyJVbQtJnutiGxjzAvAN7FST0rFhfYpKGVZC4TsNNBDwN1YqZtVdmdvOda39NYWA18WkbVYyx++F/XaImCtiKwyVplneAYLxwAAAFlJREFUx3+wlo5cg1Xh9bvGmEN2UIklE3hGRFKwWhk39ewSleqcVklVSinl0vSRUkoplwYFpZRSLg0KSimlXBoUlFJKuTQoKKWUcmlQUEop5dKgoJRSyvX/ATdZxAe2uvDCAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEGCAYAAACKB4k+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO3deXxcVfn48c+TZDLZ96RN031faaEFyg4tuwiIG+gPQVG+KiqiqIAi7uLyRcHlK5VVBFRAAVnKTlvW0oXu+55uSbNnksx6fn/cJZNksjaTNMnzfr36ysydO3fuDeE+c55zznPEGINSSikFkNDfJ6CUUurYoUFBKaWUS4OCUkoplwYFpZRSLg0KSimlXEn9fQJHo6CgwIwdO7a/T0MppQaUlStXHjHGFMZ6bUAHhbFjx7JixYr+Pg2llBpQRGRPe69p+kgppZRLg4JSSimXBgWllFIuDQpKKaVcGhSUUkq5NCgopZRyaVBQSinl0qCglFIDzN2vbmPlnsq4HFuDglJKDSBr9lXzu1e3smzbkbgcX4OCUkoNIE+uLCU9OZEvnjE+LsfXoKCUUgPIwZpGRuenk+GNT5WiuAUFEXlARMpEZH3Utt+IyGYRWSsi/xGRnKjXbhWR7SKyRUQuiNd5KaXUQHa41s+wLG/cjh/PlsJDwIWttr0CzDTGHAdsBW4FEJHpwJXADPs9fxaRxDiem1JKDUiHa5sYlpkSt+PHLSgYY5YCla22vWyMCdlP3wNG2o8vA/5hjPEbY3YB24GT4nVuSik1EIXCEY7UD9yWQme+ALxoPy4B9kW9VmpvU0opZavwBYgYKMoagC2FjojI94EQ8KizKcZupp33Xi8iK0RkRXl5ebxOUSmljjlH6v0AFGQMopaCiFwDXAJ81hjj3PhLgVFRu40EDsR6vzFmkTFmnjFmXmFhzIWDlFJqUGoMhAFI98avy7VPg4KIXAh8D7jUGNMQ9dKzwJUi4hWRccAkYHlfnptSSh3rGoNWUEjxxC8oxG05ThF5HDgbKBCRUuAOrNFGXuAVEQF4zxjzZWPMBhH5F7ARK610gzEmHK9zU0qpgagpGAEgJWkABgVjzFUxNt/fwf4/B34er/NRSqmBzmkppCbHL8mjM5qVUmqAaLKDgjeOLQUNCkopNUD43ZaCBgWllBry+qKjWYOCUkoNEM0dzdqnoJRSQ15jMIwnUUhK1KCglFKD0rrSGraX1Xdp36ZgOK7DUSGOQ1KVUkp17qN/fAuA3Xd+pMP9jDFU+QKkxLGTGTQoKKXUgHD/W7t4+sMDpMaxkxk0faSUUgPCyxsPA80jkOJFg4JSSh0DmuuDxlaYGb/KqNE0KCilVD+JDgS1jaEO9oTCOJbLjqZBQSml+okz7wDgiM/f4b7hSMctid6iQUEppfpJQ6C5dXCkruOg4PQl/Orjs+J6ThoUlFKqnzQEmjuND9U2dbhvYyDMhMJ0Pn3i6LiekwYFpZTqJ9FBYU9FQwd7Wi2FeBbCc2hQUEqpfuKLSh/trvB1uG9DIESaJ/5TyzQoKKVUP3HWXE5MEHYf6TgoNAbCcZ/NDBoUlFKq3/j8Vkth5ogsth2uxx9qf2JaYzBMWpxnM4MGBaWU6heRiOH6R1YCcNmcEur8IV7bVNbu/g2BMGnaUlBKqcEpulzF+TOGkelN4t0dFe3vr+kjpZQavPyh5olrxdmpjM5PY19V+yOQNH2klFKDmNN/cOcVs0hMEEbnpbG3MnZQMMZYQUFbCkopNTj5naU17W//o/PSKK1sJBKjnIU/FMEY8GpLQSmlBicnfeS111selZdGIBzhcF3bmc2tA0g8aVBQSql+0GR3NHs91m24ODsFgMO1bWsgOakmJ4DEkwYFpZTqB80tBevbf4FdGrs8RmG81q2KeNLlOJVSqh+0/vbvLKITHRSMMURM1L6aPlJKqcHJ6SdwWgr5GckAHKlvDgp3vbKVCbe9QG1TyN5X00dKKTUouSkhu0/Bm5RITpqnRUvhvmW7AKioD9j7DOCgICIPiEiZiKyP2pYnIq+IyDb7Z669XUTkHhHZLiJrReSEeJ2XUkodC2J1HhdmeFsEhbC9XGelvSqb06qIp3iGnYeAC1ttuwV4zRgzCXjNfg5wETDJ/nc98H9xPC+llOp3Ta3SR2B1Nkenj5w5C0ecloJnALcUjDFLgcpWmy8DHrYfPwxcHrX9b8byHpAjIsXxOjellOpvTkshJepGn5vuobox6D53WgqDIn3UjmHGmIMA9s8ie3sJsC9qv1J7Wxsicr2IrBCRFeXl5XE9WaWUipfWQ1IBslOTqW4IuM/tmEDFIEkfdYfE2NZ2rjdgjFlkjJlnjJlXWFgY59NSSqnec7i2ib32spvO6KPkqG//uWkeqhuCGNPy9jeYWwqHnbSQ/dMpHl4KjIrabyRwoI/PTSml4uoj97zFmb95A2MM/lAYT6KQmND8nTgnzUMoYvAFWi624/QzDOg+hXY8C1xjP74GeCZq++fsUUjzgRonzaSUUoOFc3NftbcKfyjSJh2Uk2rNVajyBVoUxqvwOS2FAZw+EpHHgXeBKSJSKiLXAXcC54nINuA8+znAC8BOYDvwV+Cr8TovpZTqLyPs+kYf7quhKRhukw7KSfMAUNMYpD4QcrdXOC2FgVzmwhhzVTsvLYyxrwFuiNe5KKXUsSAx0UoV+fwhu6XQOijYLYWGABLV0+o0GgZ0UFBKKdWS07ns84doCobblMLOtVsKVQ1BKn2BFq8lJyUgEmtMTu86VkYfKaXUoOfzWymhen+IxkCY1FYrqWWmeNz9dpT7SBA4e4o1yjI5sW9u1xoUlFIqjowxPLf2AHVNQXdUkc8foiHQdnlNJ0j4/CF2ltczMjeNS2ePAKxA0hc0KCilVBw9ubKUrz22mkVLd7rb6v1hGmOkj5wg0RAIU9cUIjfNw4Uzh/fp+WqfglJKxdETK0uB5hnMAPX+II2BMMOyvC329SQmkJyUQEMgjD8UJjkpgbTkJB774skcrGm7TGc8aEtBKaV60SsbDzP19hdpsIeUOj+jO459/jANwRBpyW2/l6cnJ9IQCBEIRdzZzqdOLODjc0f2wdlrUFBKqV7125e20BSMsMcuZ+E4XGt900/xJLBufw37KhvbdDQDpCUn4fOHCYTbTm7rCxoUlFKqFyXYZSvC9uSCgJ02ctZJKMpMcfdNjbG8Zro3qqXQRyOOomlQUEqpXuTcx//w+jZ8/pAbFJyWQkNUXaPWo48AUpOT8AXCLdJHfUmDglJK9aJEe4LZSxsO89uXt7hBoarBWiehKdgcFGKlj9KTE2mwZzxrUFBKqQEuIarqaXVDkEA40uL1RVfPdR+nxUgfpWlLQSmlBjZjDI12WigpKigEw5EWQ1GTExM4ZUI+V5xgrSHWep4CWH0KjXafQl/UOmpNg4JSSh2lR9/fy7QfLmZ/dSMJUfWJwhHTIijkpnsQEbeDuXUrAppbCv6wthSUUmpAWrOvGoBFS3aQlBjdUjBunwJAXro1Wc0JCo2tFtMBq0/B6aD26ugjpZQaeAozrZv9pkN1LVoK0Z3KAHnpVsG7CUUZAORntJzRDNaIJGeEUn+0FLTMhVJKHSWnNbC3ooGpxZnu9qoGaxZzenIivkCYmkZrBNKVJ46iKNPLgqlFbY6V5m2+LevkNaWUGoCcvoFDtU0tUkJl9oS1604fB8AF063idiLCwmnDYq6PkB41TFVbCkopNQBF9xvsKPe5j51ZzMU5qWz7+UUtRia1J7oeknY0K6XUAPLhvmrKaptaBIUj9nrK0ZITE/Akdm3ltHRvVEuhHzqataWglFI9dPmf3ibDm8RZUwrJTfO4s5Zb6843fm0pKKXUAFbvDxEMRRiWlUKGN/b37O7c3KNbCjp5TSmlBiCrzHUCo/LSACjOTuG5r5/uvt6doJDq0ZaCUkoNOKGo2chOnaIxdlA4bWIBM0uy3de7MwmtRZ+CBgWllBoYostXOEFhZG4q0Hwzd/qVvZ6e9Slo+kgppQaIFkEhHMGTmECJHRR8fmsJzvz0ZACSE7s+CS26pZDeTh9FPGlQUEqpHvCHmieprS2tITkxgRE5VlBw5ifk27WOujAS1ZUSNYt5QmFGL5xp92hQUEqpHvAHW1Y4TU5KoMQOCs5chfnj8wBaFMnrTPR6DB6dp6CUUgNDdPoIrKBQnG2tv+z0C/zgkulcOLOYqcOzunXsuWNyOW1iQe+caDf1S1AQkZuALwIGWAd8HigG/gHkAauAq40xgf44P6WU6kx0+gisTuH8DC+/+/Rs5o/PB6xv+qdMyO/2sZ/6yqm9co490edtExEpAb4BzDPGzAQSgSuBXwG/M8ZMAqqA6/r63JRSqqsCrVoKwbAB4GPHj6Q4O7U/TqlXdCkoiMgEEfHaj88WkW+ISM5RfG4SkCoiSUAacBBYADxpv/4wcPlRHF8ppeLKSR9dcby1tGZjsO2COQNRV1sKTwFhEZkI3A+MAx7ryQcaY/YDvwX2YgWDGmAlUG2MCdm7lQIlsd4vIteLyAoRWVFeXt6TU1BKqS4JRwxPrSwlHDFtXnPSR84CO/4hFhQi9g37Y8DvjTE3YfUBdJuI5AKXYQWWEUA6cFGMXdv+VwCMMYuMMfOMMfMKCwt7cgpKKdUl//xgH99+Yg2PvLu7zWvO6KP8DGsuQlOw7XrLA1FXg0JQRK4CrgGes7d5eviZ5wK7jDHlxpgg8G/gVCDHTicBjAQO9PD4SinVK2qbrKqnB2qa2rzmpI+yU61b4VBLH30eOAX4uTFml4iMA/7ew8/cC8wXkTSxiosvBDYCbwCfsPe5Bnimh8dXSg1RGw/U0hAIdb5jFzllJlqvtQzN6aPs1OR29xmIuhQUjDEbjTHfMMY8bj/fZYy5sycfaIx5H6tDeRXWcNQEYBHwPeBbIrIdyMfqu1BKqS6p94e4+J5l3PzEmg73W7z+EP/6YF/MfoLWUjzW7OLYQcFqKUwssmYdXz4nZjfogNPhPAURWUc7uX0AY8xxPflQY8wdwB2tNu8ETurJ8ZRS6nCtleJZuacq5us7yuupawrx5b+vBGBsQTonjcvr8JjGvvvF6i9w+hSGZXnZ+rOL8HRj1vKxrLPJa5fYP2+wfz5i//ws0BCXM1JKqR44ZOf9o6uMRrv1qXUs313pPq+IsWxma06KqPVEteht3qTEfilxHS8dBgVjzB4AETnNGHNa1Eu3iMjbwE/ieXJKKdVVB+2gkOqJXZG0pjHY4fNYnBZCzJZCKIIIg6aF4OhqmYt0ETndGPMWgIicijWUVCmljglO+igtuW1QWL+/xi1S56juICj8/tWtVDcEyUlrf2RRvT9EenIS0p0SqANAV9s8XwD+JCK7RWQX8Gd7m1JK9bt9lQ385qUtQNtCdQCX/OEtKnwtS6n964N9zLrjJRoDLW/4oXCEh97ZzRtbytwWQl1T2xFNVb4Auek9HZl/7Oq0pSAiCcBEY8xsEckCxBhTE/9TU0qprtld4XMfO3MLHMbEHiuz84j1nkO1TYwraE58rNpbTXVDkEjEuKOOymrbzlOoagiSl5Z81Od+rOm0pWCMiQBfsx/XakBQSh1rgvZ6ybNH5VDbKi0U3XJwMj1j8tPcbc4qaY7XNh0GoLYp5M55qPAFqGloedyqhgA5QzEo2F4RkZtFZJSI5Dn/4npmSinVRYGQ1RooSE+mtinUonVQH3XTv/Wiqez65cUUZHjdba07nN/YUuY+PlTb3A+x40h9i/0qfQHy0gdfUOhqR7PTf3BD1DYDjO/d01FKqe5zWgr5GcmEI4Z6f4jMFCvfH90SyEzxICLkpDb3BbRuWZRWNVKU6aWszs+hmkYSBCIGdpb7OGF0rrtflS9A7lBtKRhjxsX4pwFBKXVMcNY2GJVrpYUOVDf3AUR3Eqd7re/B4wub+xCiWwqRiKEhEGZUnnWcgzVNTCrKBKzObIc/FMYXCJM3CDuauzzjQkRmisinRORzzr94nphSSnWV01KYVmwte7m9rDnVs7a0uRvUSSvdetE0Xv3WWUDLjmmf3YcwKtdaJKeuKURGShI5aR4qo0YvVfms9wzGPoUupY9E5A7gbGA68AJWqeu3gL/F7cyUUqqLnKAwZXgmIs1BYXtZPbf9Z527X8QOCgkJwoTCdBITpEVLwee3Rhs5LQWwl9lMT6bC19y/8Nb2I+7nDTZdbSl8Aqua6SFjzOeB2YC347copVTvC0cM9y7ZQXVD8zf3gL0UZlaqh5G5qewot4JC9Lf7k8flcdHM5mVgRISslCRqG5vTS06n9PjCdDJTrO/MKZ5E8jO8HKlvPtbzaw8wJj+NeWOa+xgGi652NDcaYyIiErLnKpShncxKqX7wwrqD/PLFzZTV+TlY00hdU4jTJhYAkJyYQHFWKmV1Vp+CL6qM9t1XHu9WPXVkpXpatRSs/bNSPJw8Lo9XN5XhTUogxZPMlkN17n6Hav1MKsocdLOZoesthRX2msx/xVo6cxWwPG5npZRSWEXuom/GAGtLqwFI8STwwrpDLNt2hKDd0exJFHLTPW7OP7o1ke5tW/4iJy2Zqqh9nKCQ7k1yWxUjclLJT/eyo9zHOzustFFtY9BdXGew6eroo68aY6qNMX8BzgOusdNISinlOu3O1/nsfe/12vEW/u+bXPD7pS22bbaDxKGa5hx/MGwVp0tMEPLSk6lsCHCgurHFPukxqqcWpCdTEZUWctJHGd4kPj53JGt/dD7fv3gaoYgVdD7z1/d5b2cF+6sb3bpIg01XO5r/BiwDlhljNsf3lJRSA9X+6kb2Vzf22vF8dl2ihkDILYnt3MTf2l7u7ucPRfAkJiAi5KYlc6Tez6l3vu6+fsakAhIS2qZ68jOSWX+geXSSk25yhq5m2XMd5o/P5/Hl+wC4cpEV9IZ0SwF4CCgG/iAiO0TkKRG5MX6npZQaKBoDYe56eUtcl6PcfaR5joAzhPRw1GzjqoYAyYnW7Sw3LZnockdZKUk8ct3JMY9bkOGloj5AxF6Frd4efdQ61XTZnBK+evaEFtuGdEvBGPO6iCwBTgTOAb4MzADujuO5KaUGgIff3c09r28nI6Wr41a6b9cRH9NHWHMQYlUsLa/zu+sa5LYqPdHRXIL8DC+hiKG2KUhOWrLbp5DhbXstk4ZltHg+WFsKXU0fvYa1fsK7WGmkE40xZR2/Syk1FNTbN+n9Vc1po3DEkBgjXdNdmd4k6vwhd96BMYa6prbrIBypD+CxWwqtZxl39I2+ICPZfr+fnLRk6pqCJEjshXpGR81dAGvk0mDU1fTRWiAAzASOA2aKSGrczkopNWBU2qN3dpQ3l6/+2fMbe+XYQbuDd4Od9/cFwkRiVMI+Uu93g0J0sTuAEdnt36qcfZ05CLuO+BiVlxZzqGnr46a1s8LbQNfV0Uc3GWPOBD4GVAAPAtXxPDGl1MCwt8LK93+4r/mW8ODbu908fU/5Q2F3kZsNB2qB5uJ1w7Ja3qDL6/zuOsmzSrK5+8o5nDTOKuR83Kjsdj/DudE7ndfr99cyw05TtdY6DZU9SPsUuhQURORrIvJP4EPgcuABrFIXSqkhrrTKCgr1rdYlaL38ZXc5fQclOansr26kpjHobhuT13I14FDEuH0KIsJlc0rciqbHleS0+xn5Uemj2qYgeysbmDEidhDJiuozueOj05k6PHbwGOi62jOUCtwFrDTGtO3lUUoNWc6IndZKqxspykrp8XGdVsHkYRnsr250b9xgLZKzfHclw7NSOGSviuakjxw3LpzE7JHZnDYxv93PyE1LJkGgot5PmT2aaWRu7HRTdErp7ClFPb6uY11X00e/ATzA1QAiUigi4+J5YkqpgcEfCpPiab6VLLp6LmCtS3A03FZBvtUqqPIF3E5mZ+W0/IxkN23UOiikJidy0aziDktROJPdyusD7pDaWJ3MrcWaHT1YdDV9dAfwPeBWe5MH+Hu8TkopNXD4gxHG5jenc2aNtNIv+48yKDjrLk+1K5FW+gJu8Tqnimm6N8kdGuoEh+7KT/dSUe93g0Lr+kixZHoHZ38CdH300ceASwEfgDHmADD4asYqpbolHDEEwpEWC98XZ6dSlOll2+G6Dt7ZuTe3lJOXnuwWu7v+kZWsseseTSi05gykJye6q6glJ/YsKBRkWjOgnU7trgSF6JbRYNPVKwsYa3UKAyAi6Z3sr5QaAvwh69v1WDsoOMHhuJHZ7g28p1btrWL++LwWQ0EffHs3iQnifk50S8HpaO6u/HQvFb5AVEuh89viYKyO6uhqR/O/ROReIEdEvoS1ZvN98TstpdRA4Hy7Lsr0ctenZrvf6meV5PDa5jLqmoLuWsnd1RgIk5XiITW55Tf33DQP6d4kUjwJZHiT3MlprfsUusopddEU6jx99MbNZ7P7iK/d1weDrnY0/xZ4EngKmAL80BhzT08/VERyRORJEdksIptE5BQRyRORV0Rkm/1z8K1eodQgE52Hv+KEkQyzRxsdNzIbY5rnF/REIByJ2U+Qn261HL6+YBKXzh7hziz29LRPISOZen+I6garEzslqf2gMK4gnXOmDt6RR9CNNZqNMa8YY75jjLkZeF1EPnsUn3s3sNgYMxVrFbdNwC3Aa8aYScBr9nOl1DGsvZSL09m8NiqFFIkYGu2qp9vL6vjhM+sJ2ctoxhIIRdx+gujUUJ5d2+iGcyZy6sQC95t9RzfzjhTa6Smnuutg7i/oig6vXkSyRORWEfmjiJwvlq8BO4FP9eQD7ZXbzgTuBzDGBIwx1cBlwMP2bg9jTZJTSh3D3M7ZVjfkggwvJTmprC1tLkt95+LNTPvhYgKhCN//z3r+9u4e1u6voT2BUHNLYe0dF1CSY80fyMtoObPYKbH08bklPboGZwLb82sPWteSPHiHm3ZFZyHxEax00Trgi8DLwCeBy4wxl/XwM8cD5cCDIrJaRO6zO66HGWMOAtg/Y7bRROR6EVkhIivKy8tj7aKU6iMd5eFnlWSzLuqm/+h7ewCrzLXTT9DeCKVIxBCKGDcopCYnMjzbSk3lt6qC+q3zpvDgtSdy6oSCHl2D05G9t9Kamd3TFsdg0VlQGG+MudYYcy9wFTAPuMQY8+FRfGYScALwf8aY47GGuXY5VWSMWWSMmWeMmVdYWHgUp6GUOlpO+sgbI+Vy3Khs9lQ0sL2snqZg2F3kpqI+QJodFNa101II2Gml6D4Fp2zGxKKWJazz0pOPKs+f36rl0dNRTINFZ0HBrVFrjAkDu4wxRzf4GEqBUmPM+/bzJ7GCxGERKQawf2ppbqWOcf5Q+2P7nZpD5961hO89tdYtpV3h87tVSffYxfTKapsYe8vzvLG5rMVxo+ceOAHIWTu5t7SufjqYh5t2RWdBYbaI1Nr/6oDjnMci0qNhBcaYQ8A+EZlib1oIbASeBa6xt10DPNOT4yul+o7f6WiOkXI5fnRzIbq3t1eQaN9s395ewfJdlQBUNwS5d8kOTv7lawD83U4xBeyg4I1qKdx/zYn87ydnU5jZ8iZ+tFI8icwe2X4l1aGmw6BgjEk0xmTZ/zKNMUlRj4+mRODXgUdFZC0wB/gFcCdwnohsA86znyuljmHNs4Db3krSvUk8cO08AMYXpLsthb8s2eHuc6C6kV++uNldPtNJF8VKH80syebjc0f2/kUAj18/Py7HHYjit35eB+w+iXkxXlrY1+eilOq5zuoFLZg6jI8cV+yO7Gmtwhdo8dxpGTgtkJ7WM+qurhTBGyqG9oBcpdRR6UoRuYL0tmskTyzK4MaFk9psP1zr55v/WE2VPZHM20cjgYZ6P0K0fmkpKKUGvsZAmB/911p2s6MJX7kxgsK04iyKstr2Dby7swKwi6zR8yJ3quc0KCilesRZ3AY6Tr846yIA7qI4WSlJ5KW1DRaO3faopL5KHwFcOnsESQnaYtAwrJTqkRp7ZbQHrp3XYfrlrMnWfKIXbzyD+eOtdZOzUj2kea3vpOML0nn7lgWcbK+pDLDXXkuhL4PCPVcdz12fntNnn3es0qCglOoRJyg4pavbc+bkQnb84mKmFWe5gSArxcOckTmcPrGAv14zj5KcVKaPaB7Q6PQp9GVQUBZNHymleqS6wRo51FlQANzhqNEF7rLTPPz9iye7+8RKJ2mfQt/T37hSqtt2H/Fx4z+sajdZXQgKDqeERChi2ryWE6ND2qsthT6nv3GlBpgfPbuBlzcc6tdz+MlzG93HXWkpOJyhq6ZtTHBbCmPy09xtmj7qe/obV2oAqfIFeOid3Vz/yMpeOd4j7+3p0VrKh6NGHnVnLsEXzxjPJ+aO5LPzR7d5LddeQS3D25zV1qDQ9/Q3rtQAsnpfVa8dKxSOcPvT63ls+d5uv9dZkKa7slM9/PaTs8mKsURnjt1SSI8OCtqn0Of0N67UALJmn1VqOvrbdE/V2vMHyur83X6vzx/qfKduclZU05ZC/9LfuFIDSJU94qfeH3JLTPSUM6S0LCoV1BWBUIRg2HDGpAIevPbEozqHaDl2+ijdm8TpE60Fc/qqzIVqpkNSlRpA6qO+oe+vbmRCYUYHe3fMGVLa3ZaCs87yWZMLe3UR+xRPIqmeRDK8ifzq43PZWe7TlkI/0KCg1ADS4G9uHZRW9TwoRCKGB9/eDVidxsaYLheFawhagSm9F1JYrf3gkmnMHJFNWnISM0t0jYP+oEFBqQHEFwgxLMvL4Vo/pVUNPT7O8+sO8uyaA4C1JkJtU6jLQ0sb7JZCWhwWuP/syWN6/Ziqe7RtptQA4vOHGF+QQVKCUFrVsxFA0Hb0UHld1/sVnPSRrkEwOGlQUGoAaQiEyUxJYkRO6lEFhV3lvhbPD9d23q8Qjhi2l9W5LYV4pI9U/9OgoNQAsLO8nm/960M2H6ojw5vEyNzUo0ofbS1rOWHtcBdGIN27dAfn3rWUlXusuRKpcUgfqf6nQUGpAeAXL2zi36v2A5DmTWR4VgplXfh2357axiBpyYl88fRxQNdGIG06aAWS9+yFcOLRp6D6nwYFpY5xxhhe31zmPk/3JlGY5aW8zo+JVUSoC/yhCBfOHM4PLplOhjepSy2FggxrctmSreXWeSRr+mgw0qCg1DGuMRgmYppLPiQnJjAsM4VAOEK1ve5AdzUFIzHd/EAAACAASURBVG5xuqIsb5daHQmthqxq+mhw0qCg1DHOZ89NKMlNBayZyM76xoe7MWoomj8UdstSl+Sksrey8/6JuqaWAUjTR4OTBgWljnENAWuyWEmOFRSqGoIUZaYA9LhfwR+KuCUkphdnseVQHYFQpN39X9pwiH+tKGVEdoq7LUVLUAxKGhSUOsY5LYWF06ySElccX0JRpt1S6GbdIrBmMwdCEVI81v/+s0ZmEwhHuO+tne2+53/sUt2FWSl854IpDMvykqCL3A9K2lOk1DHOaSlMLMpg950fASAYjuBNSmDzoe6vhRAIWy0Cp6Uwd0wuAL9evIUvnj6+w3pDjYEQN5wzkRvOmdjtz1UDg7YUlOpHXRk95HPLSjR/h/MkJnDcyGxW7+3a+gqRqOUv/UEnKFj/+xdnp/Kzy2cCsKfC1+a9ZVH9Fj3t2FYDhwYFpfrJnS9u5uRfvMbBmkYq6tvvG2jwOwXoWubwTxidy7r9NWw6WNvh59Q2BRl/2ws89PYuAJpCVpBJiSpTMXtkDgA7yuvdbX98fRtf+tsKVu5uDjxOuW01eGlQUKqP3fDYKr775Br+smQHZXV+zvz1G8z92auEwm07eneU1/OVR1cBbecFXHvaWNKSk/jLkh0dfl5FvVUi+87Fm4G2LQWA8YXpAGwvaw4Kv315K69sPMzSbUfc4bBHU6pbDQz91qcgIonACmC/MeYSERkH/APIA1YBVxtjAv11fkrFy/NrD7Z4HgxbqZ0d5T6mDM/kuoc+YNXeKj514iiGZzWP9mk9BLQ4O5UZI7LY18lw0np7hbUmOxj47ZaC19McFNK9SYzITmFHedv00dOr9zN/Qj5fPms8k4oyu3qZaoDqz5bCjcCmqOe/An5njJkEVAHX9ctZKdVP1pRWEwpHeG1zGVUNQe5dstMtPgexC9CV5KR2ul5ynb855XOk3u8Gh9ZDSicUZbjpo+hWS2MwTElOKqdOKKDQHvWkBq9+CQoiMhL4CHCf/VyABcCT9i4PA5f3x7kpFU/hSMuO5cnDmtMxGw/Utql8uvVw8+gib4xRQSW5qZTV+TucY+C0FAD++cG+mC0FsFJD28vquW/ZTtYfaNlPocFg6OivlsLvge8Czl9yPlBtjHH+ekuBklhvFJHrRWSFiKwoLy+P/5kq1YtqW3XUnjQuz328u8LHriNW+sZZvP79nZXu67FWRivJScUYWL6rss1rjjo7KBRkeHlx/UH8dgBJ8bRtKTQEwvzs+U384Ol1LV4rtOseqcGvz4OCiFwClBljVkZvjrFrzLF6xphFxph5xph5hYWFcTlHpeKlulVQmDI8y328t7LBTd/86uPHAXDInpx29pTYf+uTh1k5/pufWONuW76rkipfc3ecs67zx44fwfr9teypsPogWrc8zp1WxEdmFTOzJIv1+1u2FAoytKUwVPRHS+E04FIR2Y3VsbwAq+WQIyJO0nQkcKAfzk0NQjvK6zsc8tmXWq+BcI59s59YlEFpZSM7yuvJTvUwfURzsPjSGeN46PMnxTze7FE5XHPKGA7VNlHTGCQSMXzq3ne5+J5l7j5OzaKPHT8SgP+sLgWaJ685irNT+dNnT+CTc0e525wWi6aPho4+DwrGmFuNMSONMWOBK4HXjTGfBd4APmHvdg3wTF+fmxqcvvjwCn78341xOXZZbRP/88gKd+GZjlTU+7n6/uUAnDttGKdOyGdkbhq77/wInz9tLIFwhHd3VDCuIJ3iqBpDTp2j9pwz1Sp/selgLT579vPBmiaMscpZ/HPFPpITE5g+IosZI7L4wJ53kOKJ/b9/dEDKSrGCQm66po+GimNpnsL3gG+JyHasPob7+/l81CBRXufv0k27J/735a28tOEwD72zu9N9n/6wufF7y0VTeexL893n4wqseQK7KxoYX5jeIt/vVERtz4wR2QCsK61x+w8A9lU28tA7u9hX2eiWtvjxpTPc11u3FBxThjcPO/35FbMYX5DuFuNTg1+/1j4yxrwJvGk/3gnEbiMr1UPhiKHeH6LeH6K8zt/raRAn579hf02n+y7fVeE+zk3ztHhtyrDmG/F4O0AkJyYQCEco7CSfX5jpZVJRBq9vLuPMyc19Dyv3VnKkvuVUn3lj87jkuGKeW3uQjJTY//tnpVjndvK4PM6ZUsQ5U4o6vTY1eGhBPDWoRQ/HXLOvmnOnD+vV49fa+fqdR3xU+QJkpXpIbKd6aJUvyMySLL513mTyW93oo5/PHWONSMpISaLSF+hSIDtv+jDuXbqTgzXNQ1rvf2tXmw5jgHuuPJ5bL57m9hfEsv7HF+BJ1CqoQ9GxlD5SqtfVRi0M8+G+6t4/vr3WMcC2snrO/PUbXPrHt3h3R0WbfSsbAozKTWPB1NiBKdP+5n6yPUz1a3Yl0uHZHfcpgFXpNBwxfLDbGpqan57cIiCs/MG57uOEBOk0HZThTWo3vaQGNw0KalCLLuC2pjQOQaEpxAmjrdLTS7eWs7+6kbWlNdz54qY2+1b5Ah122L5805ks+c7Z7joFnz9tLLt+eTGZKZ523+OYVmx1Dn+wy+o7+e6FU9zXPjF3ZJuWiVLt0fSRGtScjtdxBel8uK+aSMT02uIwkYihtjHItOJMPthdycsbDwEwdXgmO8t9GGPcCWeRiKGqIUB+B0GhOLvlt/dYk9Xaf28KWSlJLLdbCudMKeKdWxbgTUqIWR5DqfZoS0ENak766MxJBdQ1hdgVY72Anrh3yQ7G3/YC/lCEnLRkphZnsfWwNfHs/OnDqPOHWnTy1jYFiRjITYvP0E4RcVsLAJkpHkbkpJKf4W0zc1mpjmhQUIOa01I4Y5I1Kmf13pYppNV7q2KWrO7M39/f4z7OSkniC6eNBax+gbljrT6BDQeaRyRV2jOM8+I43t8JCokJ0u4cBKU6o385qs9tL6vna4+toikY7nzno+TUGjp+dA756cks29ZcL2tvRQMf+/M7TPz+i6zY3X7toFjGFzQXsstK9XDp7BEsunouj39pPjNHZJGWnMhPntuIMYb/rjnAgv9dAsR3EthUe35BOGK6lXpSKpoGBdXn7nxxE8+tPcjSrfEvaOi0FLJTPSyYWsQbm8vciqIVvubSF/cubX/R+liSovol0pKTEBHOnzGcmSXZ5Gd4+erZE9hZ7qO6IchN//zQ3XdG1Gzh3naCvdby1xfo+smq5zQoqD6XlWqNpjlcF/96RLVN1pDRpMQELpk9gtqmEC+utxa5iZ79uzNqGcquqGyw0kFJCdKi/LVjql3o7o0tZYTsctnfvXBKXAvLTR6WycofnMu3z5/S+c5KtUODgupzzrKSO8q6dyPuibqmoDtD94yJBYzOS+M/q/cDzdVDL5gxjF1HfN1KZ1U3BLnkuGK2/fwixuSnt3l9bEEaAL9evIUEgfdvW8hXz47/N3gdeqqOlgYF1eec8tHxmEzWWm1jyJ0UlpAgnDwujzX7qjHGuLOdT59YQMTAn9/Y3uXjVjUEyE1Lbjd3PzLXCgqHapu4aGYxw7I6n4Cm1LFAg4Lqc06t/w/3VXe6vvDRqvMH3aAAVqnpqoYgpVWN1NkthY/OHsGls0fwlyU7OWzXMmotHDH8/tWtVPoChCOGmsZgm/pF0VI8iVxzyhi+fNYE7vr07N69KKXiSIOC6nNVDQE3D3/7M+s7XWP4aNQ2htw+DLBGIQG8t7PCbSlkeJO4+fwpBCMR/rF8HwDGGBoCzX0OH+yu5PevbuO2f6+jpjGIMZDTyZyDH182k1sumqrlItSAokGhD/zx9W1cff/7fPMfq93x6oNRlS/A6b96nf+u6Xh9pCpfgFkl1s35zS3lfO2xVXE7p7qmYIsyEdOLsyjOTuHljYep9wdJ9Vid0KPz0zh5XB7/WV2KPxTmvmW7mP7Dl9z/Xo12f8PiDYf45wdW4OispLVSA5EGhT7w25e3smzbEZ7+8AB3v7q1v0+nUzvK61tU2+yqbWX1lFY18vXHVxOJxFxNFWMMlQ0B8tI9fMMeOrnxQC3GxN7/aNU2hdyFYsCa+btwWhHvbD9CTWPL1NK1p45jd0UDd728lf+utQLbH1/fzq8Xb+ZI1EipXy3eDNBiIRylBgsNCn1s8YZD7d4wjxUL/3cJp/zy9S7vb4zhR89u4IV1B91tO9oZ4rm2tIamYITJwzL55rmTueOj0/GHIuyu6P2+BWNMm5YCwJxRufgCYf61orRF+egLZw7nrMmFvLrpsLva2QNv7+LPb+7glY2H2xx/eLYuPKMGHw0KceZMlAK4/ZLpHK71s+VwXT+eUcd6ErB2lNfz0Du7W6w+trqdkUUvbThEUoJw3vRhJCQIJ9olITYeaC7zbIzhvmU7eTXGjbg7moIRgmFDVmrLgnDHjcx2H+880rIW0vzx+ewo91HV0DLN97J9Lj+9fCYACQJFum6xGoQ0KMRZdaN1c/npZTNYaK+l2xdDMXuqLCpNcqSLi92/uaV5ZnJyYgLZqZ6Y6wkAvL+rkuNGZrudtBOLMkhMEDYfag4Kz645wM+e38RXH1vF3a9uY00Pf1/OjOXs1JYthQmFGRRkWJ8/s6TlDOOTx1tBylm+87sXTnE7p0fnpfHJuSNJ8SRQmOnFk6j/+6jBR/+q46y6wRqTn5OWzJj8NHLTPHy4t5rSqgZ+8PQ6GgPxr//THaVVzWmcdaXNBd22Ha7jU3951x1OGi06tVKY6eWCGcN4ZeNhdzLY3a9u4/MPLqfeH2Llnirm2a0DsIZujs1P4w+vb2eX/a196dYjgNXK+t2rW7nsT29z23/Wdfta1ttLZEZXDwWrYNz7t53LmjvO5/GodZIBZpVku4vmXDp7BF89eyLHj7LKR+SmeUjxJLJw6jCmDI9fuQql+tOQDAoHaxp55L09BMMR1uyr5kAch0Q6QcGZ6DR3TB7v7qzg2gc/4O/v7WVJH9T/6Y59UUHBWR8A4BcvbGL57koWbzjUYv/yOj/Ld1cywu50TU1O5PI5JdT7Qzy/1upj+N2rW3ljSzln/+YNADdl5Lh8TgkAT6zYR70/xHNrD7TpxH3s/b088u7ubo3eWr23muTEhJj1hhIThOxUT5v+Bk9iAnPtGkJO2mnhtCLSkxP50pnjAbjr07NZdPXcLp+HUgPJkFx9Y/Xeam5/ej3Ti7O4/m8rmDI8k4c+fxKlVQ2ML2xbx6annA5YgBx7otNZkwt4dVPzN+v3dlZw4czhvfaZ3T2/XUd87K9u5C9LdnDrRdM4VGOlXM6dNozn1h7klgunkZ3mYV+VFThfXH+IK08cBcAj7+0hQQRj4K5Pz2H9/hqmDM/klAn5TCzK4K/LdnLxrGL388IRwz1XHc+CqS0Xgv/6wkm8trmMD3ZX8p0n1uAPRbj21LGMyU+nwufniuNHMu2Hi7n9mQ3c/swGZozIYs6oHERgZ7mPx1p923dsPlTHpGEZ3Z4ncPX8MSzbdsRND502sYANP7nQfV3nHajBbEgGBWcN3OfWHqDCF+CdHRVc9/AHLNt2hLdvWdDp+rVdVekLsPGglStvDgpFwAZ3n7e3H+mVz3L8+c3tzBmZw/zx+Z2uMPbkylK+8+RahmelcKi2ibte2cqkogySkxK4+YLJXHT3Mn72/EbW7a9hu12naOnWcp5atZ+Tx+Xxw2ear2NCYQbzx+e7z29cOImvP76a379mDcE9d9owbrt4artB9+Txedy7xKpUesr4fK49bWyLm+81p4zh4XetNQw2HKhlQ1TH9IYDNcwYkU1rdU3BHq1fcP6M4Sy6em6LNJdSQ8WQTB/lZ3iZMiyTB9/e7W5bts26OXc28ao7aqOqcDrVMUfnpzG+wCqg9pmTR7OtrJ6yutilFbrLGMOvF2/hM/e9zzUPLu9w35qGoDve/pBd2mHV3ioqfQFy0zxMHZ7FSWPzeGJlKZsPWaOlHrz2RHLSPCzdWt7ippyWnOh23Do+MquYzJQk/vaOdSO/5aL2AwLAF04b597Av3PhlDbfxn9wyXQeuHYeRZle/vU/pxAd7y7749v8KUbdonp/qMWQ0+44f8bwuC6Io9SxakgGBYDLjy9xH586IR+nrtn9b+2irinYzru6x1ngZdHVc1ssiXjhzOEMy/K6aZiXNhzd0EuHU8sHrCDndPQaY9pUAP3D69taLBcJVv/Hyj1V7pKRHzmuuMXrM0qyOH1iAc+uOcCX/77S3d4QCLcpDJeQIJwwOpfGYJgEgVF5Hbe+hmWlsPS75/D2LQs4YXRum9c9iQksmDqM5d8/l5PG5XHjwsn88TPH8++vnkooYnjs/b1t3lPf1POgoNRQNWSDwtWnjOHMyYVcclwxf/vCSaz/0QU89ZVTKK/z8+K6Q50foAucev2tV9u66bzJvHzTWcwqyWbumFx+98pW/KGjH4VU2eomP/X2xQRCER5+ZzdTb19MRdQQ021l9aR4Evjr5+YBuJ2rO4/43FTXVSeN5sUbz+DuK+dwwYxhFGWmuKk3aL7RO0M2WzvJ3jcpIaFLefgMb1KXU3c3njuJS44bwQmjc/nJZTPYX93I3lYT4Or8ITJSNCgo1R1DNihkeJP42xdO4o+fOYGkxATSvUmcMDqXokwvS7b1zoggZ9H4rBgjXLJTPYgIN5wzgUpfoN1x/d3hLPxyz1XHu9tW7qni6Q+tlFj0fILdFT4WThvGgqlFfOu8yfzpMye4N3ynpeBJTGBacRaXzSnh3qut4HFiVFB449tns/y2hTx07Ukxz+eCGcMACPRgDeTuONNef/mF9c0zqo0x1PtDZGpLQaluGbJBIRYR4czJhTy/9iCPvr/nqMtROOmj1jNqo506oYC05EQee38vt/57LYvX97yV4rQUxuSlse5H55OUIFz11/fcyXIPvL2LKl+AdaU1lFY1Mi4/ncQE4RsLJzE8O4XTJhYA1nDN9kwustYBLsz0kpSYQFFWCtntlJCeWJTJ/PF5/OAj03p8TV0xtiCdE8fm8tj7ewnaAaghEMYYtKWgVDfp/zGtnDW5kCdXlvL9/6ynJCeVs6cUdf6mdrTXUoiW4knk6vlj3DWCH1++j3duWcCILqZRyuqa2FXuY/2BWjcFlZeeTGaKh6+cPYE/vN7cAbvlUB3H//QV9/m4gpYrhs0ZZaWBOlrjICFBWPzNM8jrpGy04x/Xn9Kl/Y7Wl8+awHUPr+DBt3dx7anj3FXVMrzt/+6VUm31eUtBREaJyBsisklENojIjfb2PBF5RUS22T/b9jb2gdPtb8sAf3/P6ryMRAw3P7GGS/6wrFud0LWNIRITxJ0h256bzpvMDy+Zzo8+Oh2w6gN11ZWL3uPTi97jp89t5NeLtwC4o2a+ff4UnviydVM+Z0oh9197ovu+m8+fzEWzWs6PmFViDes8q5NAOHV4FkXH2EpiC6YWcfrEAn7xwmY+e997bitNWwpKdU9/pI9CwLeNMdOA+cANIjIduAV4zRgzCXjNft7nctOT+fdXT+Vzp4zh1U2H+dpjq9hyuI4nV5ayfn8tt/x7HXsrGrqUWqptCpKVktTuko2OFE8iXzh9HNeeNo6pwzO7NSx2Z7mvzbboIDRvTC4/u3wmv7hiFmdNLuSJL5/CmjvO52sLJpGW3PKGmZuezOrbz+PGhZO6/PnHChHh5x+bSXF2Ch/sruLie5YBaJ+CUt3U50HBGHPQGLPKflwHbAJKgMuAh+3dHgYu7+tzc5wwOpevnD0BgOfWHuT1zWWAtWzj82sPcuZv3uDu17a1ed9Vi95j0dId7vOqhrZlmztzxQklrNpbzc+f38hXH13Js50ECGd+wPTiLN6/bSGPf2l+iyAkIvy/+WMotss8nzg2r02BuGi56ckd9ikcy8bkp/POLQu4/szxBMNW0NaWglLd068dzSIyFjgeeB8YZow5CFbgAHqezO8FxdmpLPvuOQD85iUrLXPnFbPcmjf3v7WL6oYA+yob+Ml/N3KgupF3d1bwixc20xQM0xQMs3RrObNHxR6u2R6nDtBfl+3ihXWH+OY/VrNkaznGGN7efqTNfINQxFCcncID157IsKwUTpmQH+uwQ4aIcPP5U0hOsv60CzO0vLVS3dFvX6NEJAN4CvimMaa2sxRL1PuuB64HGD16dPxOEBiVl8Z1p4/j/rd2AZDuTeL8GcN56ZtncuHdS7lv2S52lNfz4vpDPPr+Hvd9S7eWc6Q+QE1jkE/PG9WtzyzKSmFMfhp7Khp479aFXPvgcq55oHl2cm6ah7994WRmjcymKRimuiHIF88fx3BdBcyVnJTA6tvPY19VA2NbdaYrpTrWL0FBRDxYAeFRY8y/7c2HRaTYGHNQRIqBsljvNcYsAhYBzJs3L+5LmN1+yXSC4UiLlMqU4ZlcPKuYP0aVVvBHLaZz79Kd7CyvZ+6YXE6b2P1v7k/8zymU1/sZbrcAPvfAcrf2UFVDkCv+7222/PQiyu21D461Tt9jQbo3iala3lqpbuvzoCBWk+B+YJMx5q6ol54FrgHutH8+09fn1p6fXDazzbabzp3sloYGa7GWk8flU9sY5ImVpQD89pOzO+1kjqUoK8W90Y/ISeXVb51FXVOQNftq+N5Ta9lf3chdr2x1ZxL3VgE/pZSSeC2Y3u4HipwOLAPWAc7X69uw+hX+BYwG9gKfNMZUdnSsefPmmRUrVsTxbDsWiRhuf2Y9l80pcUs67Cyv59y7lvDjS2dw9Slje/0zjTF876m1/GtFKYWZXvzBMB/84Fwt56yU6jIRWWmMmRfztb4OCr2pv4NCe3z+EOlxHArZGAhzzm/f5FBtE58/bSx3fHRG3D5LKTX4dBQUdLxeHMQzIIC1utnDXziJpz/cPyDnFCiljl0aFAaoKcMz+d6FU/v7NJRSg4wWxFNKKeXSoKCUUsqlQUEppZRLg4JSSimXBgWllFIuDQpKKaVcGhSUUkq5NCgopZRyDegyFyJSDuzpdMfYCoAjvXg6A8VQvG695qFBr7nrxhhjCmO9MKCDwtEQkRXt1f4YzIbides1Dw16zb1D00dKKaVcGhSUUkq5hnJQWNTfJ9BPhuJ16zUPDXrNvWDI9ikopZRqayi3FJRSSrWiQUEppZRrSAYFEblQRLaIyHYRuaW/z6e3iMgDIlImIuujtuWJyCsiss3+mWtvFxG5x/4drBWRE/rvzHtOREaJyBsisklENojIjfb2QXvdIpIiIstFZI19zT+2t48Tkffta/6niCTb27328+3262P78/yPhogkishqEXnOfj6or1lEdovIOhH5UERW2Nvi+rc95IKCiCQCfwIuAqYDV4nI9P49q17zEHBhq223AK8ZYyYBr9nPwbr+Sfa/64H/66Nz7G0h4NvGmGnAfOAG+7/nYL5uP7DAGDMbmANcKCLzgV8Bv7OvuQq4zt7/OqDKGDMR+J2930B1I7Ap6vlQuOZzjDFzouYjxPdv2xgzpP4BpwAvRT2/Fbi1v8+rF69vLLA+6vkWoNh+XAxssR/fC1wVa7+B/A94BjhvqFw3kAasAk7GmtmaZG93/86Bl4BT7MdJ9n7S3+feg2sdad8EFwDPATIErnk3UNBqW1z/todcSwEoAfZFPS+1tw1Ww4wxBwHsn0X29kH3e7BTBMcD7zPIr9tOo3wIlAGvADuAamNMyN4l+rrca7ZfrwHy+/aMe8Xvge8CEft5PoP/mg3wsoisFJHr7W1x/dtOOoqTHagkxrahOC53UP0eRCQDeAr4pjGmViTW5Vm7xtg24K7bGBMG5ohIDvAfYFqs3eyfA/6aReQSoMwYs1JEznY2x9h10Fyz7TRjzAERKQJeEZHNHezbK9c8FFsKpcCoqOcjgQP9dC594bCIFAPYP8vs7YPm9yAiHqyA8Kgx5t/25kF/3QDGmGrgTaz+lBwRcb7oRV+Xe83269lAZd+e6VE7DbhURHYD/8BKIf2ewX3NGGMO2D/LsIL/ScT5b3soBoUPgEn2qIVk4Erg2X4+p3h6FrjGfnwNVs7d2f45e8TCfKDGaZIOJGI1Ce4HNhlj7op6adBet4gU2i0ERCQVOBer8/UN4BP2bq2v2fldfAJ43dhJ54HCGHOrMWakMWYs1v+zrxtjPssgvmYRSReRTOcxcD6wnnj/bfd3R0o/dd5cDGzFysN+v7/Ppxev63HgIBDE+tZwHVYe9TVgm/0zz95XsEZh7QDWAfP6+/x7eM2nYzWR1wIf2v8uHszXDRwHrLaveT3wQ3v7eGA5sB14AvDa21Ps59vt18f39zUc5fWfDTw32K/ZvrY19r8Nzr0q3n/bWuZCKaWUayimj5RSSrVDg4JSSimXBgWllFIuDQpKKaVcGhSUUkq5NCioIU1E6u2fY0XkM7187NtaPX+nN4+vVDxoUFDKMhboVlCwK+52pEVQMMac2s1zUqrPaVBQynIncIZdt/4mu+Dcb0TkA7s2/f8AiMjZYq3f8BjWBCFE5Gm7YNkGp2iZiNwJpNrHe9Te5rRKxD72ertW/qejjv2miDwpIptF5FF7xjYicqeIbLTP5bd9/ttRQ8ZQLIinVCy3ADcbYy4BsG/uNcaYE0XEC7wtIi/b+54EzDTG7LKff8EYU2mXnPhARJ4yxtwiIl8zxsyJ8VlXYK2DMBsosN+z1H7teGAGVs2at4HTRGQj8DFgqjHGOCUulIoHbSkoFdv5WHVkPsQqxZ2PtXgJwPKogADwDRFZA7yHVZBsEh07HXjcGBM2xhwGlgAnRh271BgTwSrZMRaoBZqA+0TkCqDhqK9OqXZoUFAqNgG+bqwVr+YYY8YZY5yWgs/dySrjfC7Wgi6zsWoSpXTh2O3xRz0OYy0gE8JqnTwFXA4s7taVKNUNGhSUstQBmVHPXwK+YpflRkQm25UqW8vGWvaxQUSmYpWwdgSd97eyFPi03W9RCJyJVbQtJnutiGxjzAvAN7FST0rFhfYpKGVZC4TsNNBDwN1YqZtVdmdvOda39NYWA18WkbVYyx++F/XaImCtiKwyVplneAYLxwAAAFlJREFUx3+wlo5cg1Xh9bvGmEN2UIklE3hGRFKwWhk39ewSleqcVklVSinl0vSRUkoplwYFpZRSLg0KSimlXBoUlFJKuTQoKKWUcmlQUEop5dKgoJRSyvX/ATdZxAe2uvDCAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "### Cartpole training! ###\n", "\n", "# Learning rate and optimizer\n", "learning_rate = 1e-3\n", "optimizer = tf.keras.optimizers.Adam(learning_rate)\n", "\n", "# instantiate cartpole agent\n", "cartpole_model = create_cartpole_model()\n", "\n", "# to track our progress\n", "smoothed_reward = LossHistory(smoothing_factor=0.9)\n", "plotter = PeriodicPlotter(sec=2, xlabel='Iterations', ylabel='Rewards')\n", "\n", "if hasattr(tqdm, '_instances'): tqdm._instances.clear() # clear if it exists\n", "for i_episode in range(500):\n", "\n", " plotter.plot(smoothed_reward.get())\n", "\n", " # Restart the environment\n", " observation = env.reset()\n", " memory.clear()\n", "\n", " while True:\n", " # using our observation, choose an action and take it in the environment\n", " action = choose_action(cartpole_model, observation)\n", " next_observation, reward, done, info = env.step(action)\n", " # add to memory\n", " memory.add_to_memory(observation, action, reward)\n", " \n", " # is the episode over? did you crash or do so well that you're done?\n", " if done:\n", " # determine total reward and keep a record of this\n", " total_reward = sum(memory.rewards)\n", " smoothed_reward.append(total_reward)\n", " \n", " # initiate training - remember we don't know anything about how the \n", " # agent is doing until it has crashed!\n", " train_step(cartpole_model, optimizer, \n", " observations=np.vstack(memory.observations),\n", " actions=np.array(memory.actions),\n", " discounted_rewards = discount_rewards(memory.rewards))\n", " \n", " # reset the memory\n", " memory.clear()\n", " break\n", " # update our observatons\n", " observation = next_observation" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "To get a sense of how our agent did, we can save a video of the trained model working on balancing the pole. Realize that this is a brand new environment that the agent has not seen before!\n", "\n", "Let's display the saved video to watch how our agent did!\n", "\n", "> Note: Actually, pyvirtualdisplay is not supported in windows. Instead, we can use videoRecoder in gym monitoring tool." ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "from gym.wrappers.monitoring import video_recorder\n", "from IPython.display import HTML\n", "from IPython import display \n", "import glob\n", "\n", "def show_video(env_name):\n", " mp4list = glob.glob('video/*.mp4')\n", " if len(mp4list) > 0:\n", " mp4 = 'video/{}.mp4'.format(env_name)\n", " video = io.open(mp4, 'r+b').read()\n", " encoded = base64.b64encode(video)\n", " display.display(HTML(data=''''''.format(encoded.decode('ascii'))))\n", " else:\n", " print(\"Could not find video\")\n", " \n", "def show_video_of_model(model, env_name):\n", " env = gym.make(env_name)\n", " vid = video_recorder.VideoRecorder(env, path=\"video/{}.mp4\".format(env_name))\n", " obs = env.reset()\n", " prev_obs = obs\n", " \n", " counter = 0\n", " done = False\n", " while not done:\n", " frame = env.render(mode='rgb_array')\n", " vid.capture_frame()\n", " if \"CartPole\" in env_name:\n", " input_obs = obs\n", " elif \"Pong\" in env_name:\n", " input_obs = pong_change(prev_obs, obs)\n", " else:\n", " raise ValueError(f\"Unknown env for saving: {env_name}\")\n", "\n", " action = model(np.expand_dims(input_obs, 0)).numpy().argmax()\n", "\n", " prev_obs = obs\n", " obs, reward, done, info = env.step(action)\n", " counter += 1\n", " \n", " env.close()" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [], "source": [ "show_video_of_model(cartpole_model, \"CartPole-v0\")" ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "show_video('CartPole-v0')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "How does the agent perform? Could you train it for shorter amounts of time and still perform well? Do you think that training longer would help even more? " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Part 2: Pong\n", "\n", "In Cartpole, we dealt with an environment that was static -- in other words, it didn't change over time. What happens if our environment is dynamic and unpredictable? Well that's exactly the case in [Pong](https://en.wikipedia.org/wiki/Pong), since part of the environment is the opposing player. We don't know how our opponent will act or react to our actions, so the complexity of our problem increases. It also becomes much more interesting, since we can compete to beat our opponent. RL provides a powerful framework for training AI systems with the ability to handle and interact with dynamic, unpredictable environments. In this part of the lab, we'll use the tools and workflow we explored in Part 1 to build an RL agent capable of playing the game of Pong.\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 3.6 Define and inspect the Pong environment\n", "\n", "As with Cartpole, we'll instantiate the Pong environment in the OpenAI gym, using a seed of 1." ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "def create_pong_env(): \n", " return gym.make(\"Pong-v0\", frameskip=5)\n", "env = create_pong_env()\n", "env.seed(1); # for reproducibility" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Let's next consider the observation space for the Pong environment. Instead of four physical descriptors of the cart-pole setup, in the case of Pong our observations are the individual video frames (i.e., images) that depict the state of the board. Thus, the observations are 210x160 RGB images (arrays of shape (210,160,3)).\n", "\n", "We can again confirm the size of the observation space by query:" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Environment has observation space = Box(0, 255, (210, 160, 3), uint8)\n" ] } ], "source": [ "print(\"Environment has observation space =\", env.observation_space)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "In Pong, at every time step, the agent (which controls the paddle) has six actions to choose from: no-op (no operation), move right, move left, fire, fire right, and fire left. Let's confirm the size of the action space by querying the environment:" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Number of possible actions that the agent can choose from = 6\n" ] } ], "source": [ "n_actions = env.action_space.n\n", "print(\"Number of possible actions that the agent can choose from =\", n_actions)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 3.7 Define the Pong agent\n", "\n", "As before, we'll use a neural network to define our agent. What network architecture do you think would be especially well suited to this game? Since our observations are now in the form of images, we'll add convolutional layers to the network to increase the learning capacity of our network. Note that you will be tasked with completing a template CNN architecture for the Pong agent -- but you should certainly experiment beyond this template to try to optimize performance!" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [], "source": [ "### Define the Pong agent ###\n", "\n", "# Functionally define layers for convenience\n", "# All convolutional layers will have ReLu activation\n", "Conv2D = functools.partial(tf.keras.layers.Conv2D, padding='same', activation='relu')\n", "Flatten = tf.keras.layers.Flatten\n", "Dense = tf.keras.layers.Dense\n", "\n", "# Defines a CNN for the Pong agent\n", "def create_pong_model():\n", " model = tf.keras.models.Sequential([\n", " # Convolutional layers\n", " Conv2D(filters=32, kernel_size=5, strides=2),\n", " Conv2D(filters=48, kernel_size=5, strides=2),\n", " Conv2D(filters=64, kernel_size=3, strides=2),\n", " Conv2D(filters=64, kernel_size=3, strides=2),\n", "\n", " Flatten(),\n", " \n", " # Fully connected layer and output\n", " Dense(units=128, activation='relu'),\n", " Dense(units=n_actions, activation=None)\n", " \n", " ])\n", " return model\n", "\n", "pong_model = create_pong_model()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Since we've already defined the action function, `choose_action(model, observation)`, we don't need to define it again. Instead, we'll be able to reuse it later on by passing in our new model we've just created, `pong_model`. This is awesome because our action function provides a modular and generalizable method for all sorts of RL agents!" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 3.8 Pong-specific functions\n", "\n", "In Part 1 (Cartpole), we implemented some key functions and classes to build and train our RL agent -- `choose_action(model, observation)` and the `Memory` class, for example. However, in getting ready to apply these to a new game like Pong, we might need to make some slight modifications. \n", "\n", "Namely, we need to think about what happens when a game ends. In Pong, we know a game has ended if the reward is +1 (we won!) or -1 (we lost unfortunately). Otherwise, we expect the reward at a timestep to be zero -- the players (or agents) are just playing eachother. So, after a game ends, we will need to reset the reward to zero when a game ends. This will result in a modified reward function." ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [], "source": [ "### Pong reward function ###\n", "\n", "# Compute normalized, discounted rewards for Pong (i.e., return)\n", "# Arguments:\n", "# rewards: reward at timesteps in episode\n", "# gamma: discounting factor. Note increase to 0.99 -- rate of depreciation will be slower.\n", "# Returns:\n", "# normalized discounted reward\n", "def discount_rewards(rewards, gamma=0.99): \n", " discounted_rewards = np.zeros_like(rewards)\n", " R = 0\n", " for t in reversed(range(0, len(rewards))):\n", " # NEW: Reset the sum if the reward is not 0 (the game has ended!)\n", " if rewards[t] != 0:\n", " R = 0\n", " # update the total discounted reward as before\n", " R = R * gamma + rewards[t]\n", " discounted_rewards[t] = R\n", " \n", " return normalize(discounted_rewards)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Additionally, we have to consider the nature of the observations in the Pong environment, and how they will be fed into our network. Our observations in this case are images. Before we input an image into our network, we'll do a bit of pre-processing to crop and scale, clean up the background colors to a single color, and set the important game elements to a single color. Let's use this function to visualize what a single observation might look like before and after pre-processing." ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [], "source": [ "import cv2\n", "\n", "def preprocess_pong(image):\n", " I = image[35:195] # Crop\n", " I = I[::2, ::2, 0] # Downsample width and height by a factor of 2\n", " I[I == 144] = 0 # Remove background type 1\n", " I[I == 109] = 0 # Remove background type 2\n", " I[I != 0] = 1 # Set remaining elements (paddles, ball, etc.) to 1\n", " I = cv2.dilate(I, np.ones((3, 3), np.uint8), iterations=1)\n", " I = I[::2, ::2, np.newaxis]\n", " return I.astype(np.float)" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAegAAADSCAYAAACBxlNzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAdTklEQVR4nO3de7RcZX3/8ffHEBIgCSFcYgiRa6SQFANECJdlkYuEiA2uqoUfYGhpj63wUwpVLrrEKm2xclF+AjYYJBQEKajQiEKMULRcAw0hEJEIgYTERK5JAAMJ398f+5mwz2Fmzpzb7H3OfF5rzTp79n5mP9/MOZlnns/s2VsRgZmZmZXLe4ouwMzMzN7NA7SZmVkJeYA2MzMrIQ/QZmZmJeQB2szMrIQ8QJuZmZWQB2gzM+sySYdJWt7DfVwj6YLeqqlZJD0u6bC+7scDtJm1PElLJb0haZ2kVZK+L2lY0XX1Z8p8QdJT6bl9TtKFkoYUXVtXVHsTERETIuLuvu7bA7SZWeZjETEM2A/4IPDljg3SoNNrr5u9vb+SuQxoAz4NDAeOAQ4HbmpmEZI2a2Z/vWmg/mGYmXVLRDwP/AyYCCDpbkn/LOl/gNeB3SRtLWmWpJWSnpd0gaRBqf0pkv5H0v+T9Kqk30g6orL/GvvbUdJtkl6StETS3+baD5J0nqTfSVor6WFJ49K2P5E0Nz3uSUmfyj1umqQn0mOel/SPaf12kuZIeiU97leVNwmpjlsk/UHSM5I+l9vfFmk2+bKkJ8jexFQlaTzwWeDEiLgvIjZExOPAXwBTJR2ea75d+jeslfTfknZO+5CkSyWtTs/jQkmV38kQSRelWfkqSd+VtEXadpik5ZLOlvR74PuSFks6NlffZpJekLRfuv+fkn6f+rlH0oS0vg04EfhiSlf+K61fKunIXC3fkrQi3b5VSQlytZyV/h0rJf1V3T/AHA/QZmY5afCbBvxvbvXJZLPB4cCzwGxgA7AHsC/wEeBvcu0PBJ4GtgPOB34kaVSd/d0ALAd2BD4B/EtuUD8TOCHVNAL4a+B1SVsBc4EfADukNldUBhdgFvCZiBhO9mbjl2n9Wamv7YHRwHlApEH6v4BHgbHAEcAZko5Ojzsf2D3djgZm1HkajwCWR8SD+ZURsQy4Hzgqt/pE4OvpuVoAXJ/WfwT4EPB+YCTwl8CLads30vpJZL+DscBXcvt8LzAK2Jnseb4hPT8VRwMvRMQj6f7PgPFkz+MjlRoiYmZa/reIGBYRH6vyb/0SMCXV8gHgANqnL+8Ftk41ngpcLmmbKvt5t4jwzTfffGvpG7AUWAe8QjZgXgFskbbdDXwt13Y0sL6yPa07AbgrLZ8CrACU2/4gcHKN/Y0DNgLDc+v+FbgmLT8JTK9S818Cv+qw7t+B89Pyc8BngBEd2nwNuBXYo8P6A4HnOqw7F/h+Wn4amJrb1kY2CFd7Pr8M3F9j243AVWn5GuDG3LZh6bkYRxaH/5Zs8HtPro2A14Ddc+sOAp5Jy4cBbwJDc9v3ANYCW6b71wNfqVHfSCCArXM1XlDl7+XItPw7YFpu29HA0lwtbwCb5bavBqY08nfpGbSZWea4iBgZETtHxGcj4o3ctmW55Z2BwcDKFBO/QjYw7pBr83ykV+PkWbLZcbX97Qi8FBFrO7Qfm5bHkQ0CHe0MHFipIdVxItmMDbI4eRrwbIqOD0rrvwksAe6U9LSkc3L727HD/s4je0NSqTNf97NVaqp4ARhTY9uYtL1i0z4jYh3wErBjRPwS+A5wObBK0kxJI8hm/lsCD+fq/HlaX/GHiPhjbr9LgMXAxyRtCfw5WfJQ+QjhwvQRwhqywReyGX0jdqT9c9Hxd/1iRGzI3X+d7I1IpzxAm5l1Lj/YLiObQW+XBvSRETEiIibk2oyVpNz995HNqqvtbwUwStLwDu2fz/W3e5WalgH/nathZGQx7N8DRMRDETGd7I3DT0gHZ0XE2og4KyJ2Az4GnJni9GVks9D8/oZHxLTU30qyNwv5Gmv5JTBO0gH5lenjgynAvNzqcbntw8ii6RWp1ssiYn9gAlmk/QWywf0NYEKuzq0jO8CvotplGisx93TgiTRoA/yftO5Isih6l0o5dfaVt4LszU1Fx991t3mANjPrgohYCdwJXCxphKT3SNpd0p/lmu0AfE7SYEmfBPYCbq+xv2XAvcC/ShoqaR+yzyorn8V+D/i6pPHpwKl9JG0LzAHeL+nk1M9gSR+UtJekzSWdKGnriHgLWEMWHSPpWEl7pDcQlfUbyWL4Nengqi3SzHKipMrBYDcB50raRtJOwP+t8xz9FvgucL2kKWlfE4BbgF9ExC9yzadJOlTS5mSfRT8QEcvSv+VASYPJIu0/Ahsj4m3gKuBSSTukf9PY3GfltdxI9rn235Nmz8lwsjdcL5LNzP+lw+NWAbvV2e8NwJclbS9pO7LPwq/rpJaGeIA2M+u6TwObA08ALwM30z7SfYDsoKMXgH8GPhERL3bcSc4JZDO3FcCPyT5Hnpu2XUI2ON5JNqDOIvv8ey3ZgHN8etzvyQ6eqnzP+GRgaYpt/w44Ka0fD/yC7DP3+4ArIuLuiNhINqOeBDyTav8e2awS4J/I4ttnUi3/0clzdHp6/HWpr5+Tff7+Fx3a/YDsALSXgP3JYnrIDoi7iuz5fZZsAL0obTubLKa/P/37fgHsWa+Y9MbqPuBg4Ie5Tdem/T9P9vu8v8NDZwF7pzj9J1V2fQEwH1gIPEZ2kFmvnHxF7T8mMTOznpB0CvA3EXFo0bVY/+YZtJmZWQl5gDYzMyuhPou4JU0Fvg0MAr4XERf2SUdmZmYDUJ/MoJWd8u5ysnOv7g2cIGnvvujLzAyySYGy010uyX2316zf6quI+wBgSUQ8HRFvkh3ePr2P+jKzFudJgQ1EfXWVj7G0P+PMcrLTyG2STkLeBrD5IPYfvdWgPirl3UYM26rTNmvWvdaESprvtdFbd9pmq1Wv9mkNL73xNuvefFudtzRr2KZJAYCkyqTgiWqNN9eQGErnrwNmHb1/n9c7bfPbhVs2vL8/8hpvxvqqr4d9NUBX66zdh92RnYR8JsD7tt4szjp4RB+V8m5HHXxQp23m3ntfEyppvvlnTuu0zeRLftqnNVx875o+3b+1pE4nBXlD2YoD37nAlFnD7rhjQadtjt5xUsP7eyDm1dzWVwP0ctqfEm4neunUZ2ZmVXQ6KcindkNpfIZjVpS+GqAfAsZL2pXs7CzHk53vtHTyM+VGZtYDSX6mPP/MjxZYiVmPdTopyKd2IzTKZ2iy0uuTg8TSlTtOB+4gu4LITZFdrNvMrC9smhSkczofD9xWcE1mPdJXM2gi4nZqnBzezKw3RcQGSZVJwSDgak8KrL/rswHazKyZPCmwgcYDtJmZWYO6coR2T/lc3GZmZiXkAdrMzKyEPECbmZmVkAdoMzOzEmr5g8Ra7eQkeT45iZlZeXkGbWZmVkIeoM3MzEqoJSPugXqlqkb09ZWqzMysd3gGbWZmVkIeoM3MzEqoFBH3qF0nctJ1tS9abQPL9484ougSzMxKrxQDtJlZT0laCqwFNgIbImJysRWZ9YwHaDMbSD4cES8UXYRZb/Bn0GZmZiXkAdrMBooA7pT0sKS2jhsltUmaL2n+W6wvoDyzrnHEbWYDxSERsULSDsBcSb+JiHsqGyNiJjATYIRGRVFFmjXKM2gzGxAiYkX6uRr4MXBAsRWZ9YwHaDPr9yRtJWl4ZRn4CLCo2KrMeqYUEfejzyziupPGA3AG8K20vrJ8Rrpfbbknbd1H9/vQmR/lROB6IC75aZf6GPfMGsx62Wjgx5Ige137QUT8vNiSzHpGEcV/FPO+rTeLsw4eUXQZ1gX5S1V29fzeF9+7hude3aDersmsUSM0Kg6UT5hjxXsg5rEmXqr6eliKiHtZbvmMKstn1FnuSVv30f0+rk/L13ezDzMzq68UA7SZmZm154jbms4RtxXNEbeVhSPuGm3dR3F9mJlZfZ5BW9N5Bm1F8wzayqL0M2gzMzNrrxQDtCPu1uvDzMzq61HEXe36q5JGAT8EdgGWAp+KiJfr7ccRd2txxG1Fc8RtZdHXEfeHI2JS7uLo5wDzImI8MC/dNzMzsy7oi4h7OjA7Lc8GjuvsAY64W68PMzOrr6cR9zPAy2TXYf33iJgp6ZWIGJlr83JEbFNvP464W4sjbiuaI24ri76MuA+JiP2AY4DTJH2o0QfmL56+7M133iT051mh+2i8rVl3SLpa0mpJi3LrRkmaK+mp9LPuhMCsv+jRAF3j+qurJI0BSD9X13jszIiYHBGT2dyTKTNryDXA1A7rfNyLDUjdjrjTNVffExFr0/Jc4GvAEcCLEXGhpHOAURHxxXr7csTdWhxxW09I2gWYExET0/0ngcMiYmWaFNwdEXvW24cjbiuLvoq4RwO/lvQo8CDw03T91QuBoyQ9BRyV7tflg8Rarw+zXjQ6IlYCpJ87FFyPWa/wqT6t6TyDtp6oMoNu6MBUSW1AG8BQttz/UE1rTsFmdfhUn2Y2kHX5uJfBDGlqgWbdUYoB2hF36/Vh1otuA2ak5RnArQXWYtZrHHFb0znitu6SdANwGLAdsAo4H/gJcBPwPuA54JMR8VK9/fggMSuLehH3Zs0uxsysuyLihBqbPNragOOI230U0oeZmdXniNuazhG3Fc0Rt5VF6Y/i9gy69fowM7P6SjFAm5mZWXuOuK3pHHFb0RxxW1k44q7R1n044jYzK6tSDNBmZmbWniNuazpH3FY0R9xWFo64a7R1H464zczKyjNoazrPoK1onkFbWZR+Bm1mZmbtlWKAdsTden2YdYekqyWtlrQot+6rkp6XtCDdfKFnGxAccVvTOeK27pL0IWAdcG1ETEzrvgqsi4iLGt2PI24rC0fcZjYgRMQ9QN1LSZoNFKUYoB1xt14fZr3sdEkLUwS+TdHFmPUGR9zWdI64rSck7QLMyUXco4EXgAC+DoyJiL+u8rg2oA1gKFvuf6g/qrYSKH3E7Rl06/Vh1lsiYlVEbIyIt4GrgANqtJsZEZMjYvJghjS3SLNuKMUAbWbWXZLG5O5+HFhUq61Zf+KI25rOEbd1l6QbgMOA7YBVwPnp/iSyiHsp8JmIWFlvPz6K28rCEXeNtu7DEbf1LxFxQkSMiYjBEbFTRMyKiJMj4k8jYp+I+PPOBmez/sIzaGs6z6CtaJ5BW1mUfgZtZmZm7ZVigHbE3Xp9mJlZfY64rekccVvRHHFbWTjiNjMz62dKMUA74m69PszMrL5OI25JVwPHAqtzp9YbBfwQ2IXse4efioiXJQn4NjANeB04JSIe6awIR9ytxRG3Fc0Rt5VFTyPua4CpHdadA8yLiPHAvHQf4BhgfLq1AVc2UqBn0K3Xh5mZ1dfpAF3j8m7TgdlpeTZwXG79tZG5HxjZ4TR8ZmZm1oCGjuKucvWYVyJiZG77yxGxjaQ5wIUR8eu0fh5wdkTMr7d/R9ytxRG3FW0gRtx3rFjQ7v7RO04qqBLrimYexV2tk6rvACS1SZovaf6yN99p0p9jW/fReFszM6uvuwP0qkp0nX6uTuuXA+Ny7XYCVlTbQf7Sb2zuyZSZmVledyPubwIvRsSFks4BRkXEFyV9FDid7CjuA4HLIqLqtVnzHHG3Fkfc1h2SxgHXAu8F3gZmRsS3a32rpN6+HHFbWfQo4k6Xd7sP2FPSckmnAhcCR0l6Cjgq3Qe4HXgaWEJ24fTPNlKgj+JuvT7MumEDcFZE7AVMAU6TtDe1v1XSZ+5YsaDTm1lP+VSf1nSeQVtvkHQr8J10OywiVqaP3O6OiD3rPbanM+hGBuBmz2A9g+6ffKpPMxtQ0sdu+wIPAKMr14BOP3eo8ZhNB6a+xfpmlWrWbaUYoB1xt14fZt0laRhwC3BGRKxp9HH5A1MHM6TvCjTrJY64G3DSdU9x3Unjiy5jwHDEbd0laTAwB7gjIi5J657EEbf1U464zazfS+f6nwUsrgzOyW3AjLQ8A7i12bWZ9YVSDNCOuFuvD7NuOAQ4GThc0oJ0m0btb5WY9WuOuK3pHHFb0RxxW1mUPuL2DLr1+jAzs/o2K7qA7tDBB3EU2Ym/jwLm3ntfwRWZmZn1rn4ZcR918EHt7nuA7l8ccVvRHHFbWQy4iHtiWp7YQNuBFg0PlD7MzKy+UgzQZmZm1p4jbms6R9xWNEfcVhaOuGvst7uPcx+OuM3M+ppn0NZ0nkFb0TyDtrIo/QzazMzM2ivFAO2Iu/X6MDOz+hxxW9M54raiOeK2sih9xO0ZdOv1YdZVksZJukvSYkmPS/p8Wv9VSc93uICGWb/XL0/1aWYtaQNwVkQ8Imk48LCkuWnbpRFxUbMK8ezYmsERtzWdI27rDZJuBb5DdhnKdV0ZoHsacZv1lgEXcU+49z7m3nvfpp/12g60aHig9GHWE5J2AfYFHkirTpe0UNLVkrYprDCzXlSKAdrMrFGShgG3AGdExBrgSmB3YBKwEri4xuPaJM2XNP8t1jetXrPu6pcRt/VvjrituyQNBuYAd0TEJVW27wLMiYiJHbflOeK2shhwEXdZY1v30Xhbs66SJGAWsDg/OEsak2v2cWBRs2sz6wueQVvTeQZt3SHpUOBXwGPA22n1ecAJZPF2AEuBz0TEynr78gzayqLeDNpfszKzfiEifg1UeyG7vdm1mDWDI273UUgfZmZWnyNuazpH3FY0R9xWFqU/SMzMzMzaK8UA7Yi79fowM7P6Oo24JV0NHAusrny3UNJXgb8F/pCanRcRt6dt5wKnAhuBz0XEHZ0V4Yi7tTjitqI54ray6GnEfQ0wtcr6SyNiUrpVBue9geOBCekxV0ga1FkHnkG3Xh9mZlZfpwN0RNwDvNTg/qYDN0bE+oh4BlgCHNCD+szMzFpSQ0dxdzx9Xoq4TwHWAPPJLgH3sqTvAPdHxHWp3SzgZxFxc739O+JuLY64rWiOuK0s+uIo7lonp6/WSdV3APkT1y97850m/Tm2dR+NtzUzs/q6NUBHxKqI2BgRbwNX8U6MvRwYl2u6E7Cixj5mRsTkiJjM5p5MmZmZ5XU34h5TOdetpH8ADoyI4yVNAH5ANmDvCMwDxkfExnr7d8TdWhxxW9EccVtZ9CjilnQDcB+wp6Tlkk4F/k3SY5IWAh8G/gEgIh4HbgKeAH4OnNbZ4Aw+irsV+zAzs/p8qk9rOs+grTskDQXuAYaQXejn5og4X9KuwI3AKOAR4OSIeLPevjyDtrLwqT7NbCBYDxweER8gO0B1qqQpwDfIzsswHniZ7ERJZv1eKQZoR9yt14dZV0VmXbo7ON0COByofJVzNnBcAeWZ9TpH3NZ0jritu9KZCR8G9gAuB75Jdu6FPdL2cWTnXphY5bFtQBvAULbc/1BNa1rdZrWUPuL2DLr1+jDrjvT1zklkX+E8ANirWrMaj9301c7BDOnLMs16RSkGaDOzroiIV4C7gSnASEmbpU01z71g1t844ramc8Rt3SFpe+CtiHhF0hbAnWQHiM0AbomIGyV9F1gYEVfU25eP4raycMRdo637cMRt/coY4K50/oWHgLkRMQc4GzhT0hJgW2BWgTWa9ZrNOm9iZla8iFgI7Ftl/dP4qnk2ADnitqZzxG1Fc8RtZeGIu0Zb9+GI28ysrDyDtqbzDNqK5hm0lUXpZ9BmZmbWXikGaEfcrdeHmZnV54jbms4RtxXNEbeVhSNuMzOzfqYUA7Qj7tbrw8zM6nPEbU3niNuK5ojbyqL0Ebdn0K3Xh5mZ1VeKAdrMzMzac8RtTeeI27pD0lDgHmAI2XUEbo6I8yVdA/wZ8GpqekpELKi3L0fcVhaOuGu0dR+OuK1fWQ8cHhEfACYBUyVNSdu+EBGT0q3u4GzWX/hqVmbWL0QW961LdwenW/ERoFkfccRtTeeI27pL0iDgYWAP4PKIODtF3AeRzbDnAedExPoqj20D2gCGsuX+h2pa0+o2q8URd4227sMRt/UvEbExIiYBOwEHSJoInAv8CfBBYBRwdo3HzoyIyRExeTBDmlazWXd5Bm1N5xm09QZJ5wOvRcRFuXWHAf8YEcfWe6wPErOyKP0M2sysM5K2lzQyLW8BHAn8RtKYtE7AccCi4qo06z2lGKAdcbdeH2bdMAa4S9JC4CFgbkTMAa6X9BjwGLAdcEGBNZr1Gkfc1nSOuK1ojritLEofcXsG3Xp9mJlZfaWYQU+aNCnmzZtXdBnWJEcccQQLFizwDNoK4xm0lUWPZtCSxkm6S9JiSY9L+nxaP0rSXElPpZ/bpPWSdJmkJZIWStqvd/85ZmZmA18jEfcG4KyI2AuYApwmaW/gHGBeRIwnnRwgtT8GGJ9ubcCVvV61mZnZANfpAB0RKyPikbS8FlgMjAWmA7NTs9lkX28grb82MvcDIytfgzAzM7PGdOkgMUm7APsCDwCjI2IlZIM4sENqNpb2x30tT+s67qtN0nxJ81988cWuV25mZjaANTxASxoG3AKcERFr6jWtsu5dR6LlT7u37bbbNlqGmZlZS2hogJY0mGxwvj4ifpRWr8qdwWcMsDqtXw6Myz18J2BF75RrZmbWGho5ilvALGBxRFyS23QbMCMtzwBuza3/dDqaewrwaiUKNzMzs8Y0cj3oQ4CTgcckLUjrzgMuBG6SdCrwHPDJtO12YBqwBHgd+KterdjMzKwFdDpAR8Svqf65MsC7vumfLqp+Wg/rMjMza2mlONWnmZmZtecB2szMrIRKcS5uSX8AXgNeKLqWKrajnHVB/61t54jYvpnFmOWl15xnKff/oQrX2HNlrq/m62EpBmgASfMjYnLRdXRU1rrAtZn1VH/4O3WNPVf2+mpxxG1mZlZCHqDNzMxKqEwD9MyiC6ihrHWBazPrqf7wd+oae67s9VVVms+gzczM7B1lmkGbmZlZUvgALWmqpCclLZF0TgnqWSrpMUkLJM1P60ZJmivpqfRzmybVcrWk1ZIW5dZVrSWd+/yy9DwulLRfAbV9VdLz6blbIGlabtu5qbYnJR3dl7WZNaJsrz3Qtf/zBdU3TtJdkhZLelzS50tY41BJD0p6NNX4T2n9rpIeSDX+UNLmRdXYqEIHaEmDgMuBY4C9gRMk7V1kTcmHI2JS7rD8c4B5ETEemJfuN8M1wNQO62rVcgwwPt3agCsLqA3g0vTcTYqI2wHS7/R4YEJ6zBXpd29WiBK/9lxD4//ni7ABOCsi9gKmAKel561MNa4HDo+IDwCTgKnpwk3fIHt9Gg+8DJxaYI0NKXoGfQCwJCKejog3gRuB6QXXVM10YHZang0c14xOI+Ie4KUGa5kOXBuZ+4GRlcuBNrG2WqYDN0bE+oh4huxCKgf0VW1mDSjla08X/883XUSsjIhH0vJaYDEwlnLVGBGxLt0dnG4BHA7cnNYXWmOjih6gxwLLcveXp3VFCuBOSQ9LakvrRlcumZl+7lBYdbVrKctzeXqK2K/OxVxlqc2soj/9TZbp9WcTSbsA+wIPULIaJQ1KV19cDcwFfge8EhEbUpMy/743KXqArnaVrKIPKz8kIvYji75Ok/ShgutpVBmeyyuB3clipZXAxWl9GWozy/PfZA9IGgbcApwREWuKrqejiNgYEZOAncjSkr2qNWtuVV1X9AC9HBiXu78TsKKgWgCIiBXp52rgx2S/3FWVuDj9XF1chTVrKfy5jIhV6T/G28BVvBNjF16bWQf96W+yTK8/SBpMNjhfHxE/SqtLVWNFRLwC3E32eflISZVLLJf5971J0QP0Q8D4dHTd5mQHEt1WVDGStpI0vLIMfARYlGqakZrNAG4tpkKoU8ttwKfT0dxTgFcrkVOzdPjM++Nkz12ltuMlDZG0K9mBbA82szazDkr12tOJ0rz+SBIwC1gcEZfkNpWpxu0ljUzLWwBHkn1WfhfwidSs6NfxxkREoTdgGvBbss8IvlRwLbsBj6bb45V6gG3Jjkx8Kv0c1aR6biCLit8ie8d/aq1ayCK7y9Pz+BgwuYDa/iP1vZDsP+yYXPsvpdqeBI4p+u/ON9/K9NqTq6nh//MF1XcoWTS8EFiQbtNKVuM+wP+mGhcBX0nrdyObGCwB/hMYUvTvu7ObzyRmZmZWQkVH3GZmZlaFB2gzM7MS8gBtZmZWQh6gzczMSsgDtJmZWQl5gDYzMyshD9BmZmYl5AHazMyshP4/A8ur9rTKjpIAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "observation = env.reset()\n", "for i in range(30):\n", " action = np.random.choice(n_actions)\n", " observation, _,_,_ = env.step(action)\n", "observation_pp = preprocess_pong(observation)\n", "\n", "f = plt.figure(figsize=(10,3))\n", "ax = f.add_subplot(121)\n", "ax2 = f.add_subplot(122)\n", "ax.imshow(observation); ax.grid(False);\n", "ax2.imshow(np.squeeze(observation_pp)); ax2.grid(False); plt.title('Preprocessed Observation');" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Let's also consider the fact that, unlike CartPole, the Pong environment has an additional element of uncertainty -- regardless of what action the agent takes, we don't know how the opponent will play. That is, the environment is changing over time, based on *both* the actions we take and the actions of the opponent, which result in motion of the ball and motion of the paddles.\n", "\n", "Therefore, to capture the dynamics, we also consider how the environment changes by looking at the difference between a previous observation (image frame) and the current observation (image frame). We've implemented a helper function, `pong_change`, that pre-processes two frames, calculates the change between the two, and then re-normalizes the values. Let's inspect this to visualize how the environment can change:" ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [], "source": [ "def pong_change(prev, curr):\n", " prev = preprocess_pong(prev)\n", " curr = preprocess_pong(curr)\n", " I = prev - curr\n", " return I" ] }, { "cell_type": "code", "execution_count": 17, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA1MAAAFhCAYAAACVlfATAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO3deZStWVke8OeFZlKGlmbQVuYhGjCiggxxShAQEGW5EhAENBEVR0YlGjXIqC4BQUwwYgRpEJwgoLIAY3SBtAMaHMCBocEWELCbphsEEdj545wLh6Lq9r3fraqz9z6/31q11qk6Qz21b39v177f851brbUAAABweq6w7QAAAAAjspkCAABYwGYKAABgAZspAACABWymAAAAFrCZAgAAWMBmamBV9WVV9TfbzgEAALvIZuoIVNVbq+qDVfX+qnpXVf1CVV39sL9Pa+1VrbV/ddivezJV9btV9aH1z3bi447HmQE4HFV1/6p67fo4fmdVvayqvnTbuU5Yz9KvOsn9X1lVH9szj156nBkB2G02U0fnXq21qyf5oiS3S/JDex9QKyP+GXx3a+3qGx/n731AVZ21jWDAqamqRyT5qSRPTHL9JDdM8t+TfN2C1/qU4/0YZ8A79syje+33IDMJgKMw4i/yQ2mtvT3Jy5LcOvn4mZ0nVNXvJ/mnJDetqmtV1c+v/2b47VX1+Kq6YlVdpaouqapbn3i9qrru+qzX9dZ/K/v3G/d93vr1L6mq11fV127c97tV9eCNz7+5ql69vl1V9dSqendVva+q/nzze56qqmpV9V1V9cYkb1x/7WlVdWFVXVpVf1JVX7bx+MdU1a9U1XlVdVlV/UVV3bKqfmCd5cKquuvG4/ddp9PNCbuuqq6V5LFJvqu19uuttQ+01v6ltfbS1tr3rR/z7Kp6/MZz9s6bt1bVo6vqz5N8oKrOOuBr51bVr1XVe6rqgqr63o3XeExV/XJV/eJ6Bry+qm67vu+5WW3wXro+4/T9p/kzPqaqfnU9Xy5N8s1V9SVVdf56Rr6zqp5RVVfeeE6rqu+sqjeu8zyuqm62fs6l66ybj/+aqnrd+vVeU1X/5jT/KAAYnM3UEauqGyS5R5L/t/HlByb5tiTXSPK2JM9J8pEkN0/yhUnumuTBrbV/TvLrSe638dz7JPm91tq793yfKyV5aZJXJLleku9J8ryqOpUa4F2TfHmSWyY5O8l9k1x0Wj/oJ9w7ye2T/Ov153+c5DZJrp3k+Ul+paquuvH4eyV5bpLPyGqNXp7Vf5efndUvez+78dh912lhTthld0xy1SQvOsPXuV+SeyY5u7X2kb1fS/KxrObSn2V1TN85ycOq6m4br/G1SV6wfvxLkjwjSVprD0zyd1mf5W+t/cSCfF+X5FfXr/28JB9N8vAk18lqDe6c5Dv3POerk3xxkjsk+f4k/zPJNya5QVZ/KXa/JKmqL0ryv5J8e5JzsppVL6mqqyzICcCgbKaOzour6pIkr07ye1lVaU54dmvt9etfPq6d5O5JHrb+2+F3J3lqkm9YP/b5+eTN1P3XX9vrDkmunuTHWmsfbq39TpLf2PPcg/xLVhu7z01SrbW/aq298ySPf/r6b2Ivqao/3XPfk1prF7fWPpgkrbXzWmsXtdY+0lp7cpKrJNnc4L2qtfby9Vr8SpLrrn+Gf8nqF6wbV9XZVXX9nHydgFN3TpJ/3NgALfX01tqFJ473fb52uyTXba09dj2X3pLk5/LJx+2rW2u/1Vr7aFZ/sfIFp5nh3I15dElV3WfjvvNbay9urX2stfbB1tqftNb+YD2P3prVBugr9rzej7fWLm2tvT7JXyZ5RWvtLa2192XVMvjC9eO+NcnPttb+sLX20dbac5L8c1azGIAdoUN+dO7dWvvtA+67cOP2jZJcKck7q+rE166w8ZjfSXK1qrp9kn/I6izPfn+bfG6SC1trH9v42tuy+tvgk2qt/U5VPSPJzyS5YVW9KMmjWmuXHvCU722tPeuA+zZ/tlTVI7M6e3Rukpbkmln9rfAJ79q4/cGsfsH76MbnyWqTeG5Ovk7AqbsoyXWq6qwz3FDtd/ztnW/nrv9i6YQrJnnVxuf/sHH7n5Jc9TRzvaO19jmnkq+qbpnkKUlum+TTsvp/4J/sec7embT3889c375Rkm+qqu/ZuP/KWc0qAHaEM1Pb0TZuX5jV32Zep7V29vrjmq21WyXJenP0y1mdYbp/kt9orV22z2u+I8kN6pPf0OKGSd6+vv2BrH55OOEzN26ntfb01toXJ7lVVnW/7zvTn219fdSjs6omfkZr7ewk70tSBzz3ZE66TsBpOT/Jh7Kq5R7kpDNjrV3O1y5McsHGMXt2a+0arbV7nGLO/V7/dOx9/v9I8tdJbtFau2aSH8yyeZSsfrYn7PnZPq219ktnkBeAwdhMbdm6TveKJE+uqmtW1RXWFzxvVk+en9V1TN+Y/St+SfKHWf3y8/1VdaWq+sqsrkd6wfr+1yX5+qr6tKq6eZJvOfHEqrpdVd1+fd3VB7L6JeujOXPXyOoap/ckOauqfiSrM1On7RTXCTgF68rajyT5maq693ouXKmq7l5VJ65Nel2Se1TVtavqM5M8bMG3+qMkl67flOJqtXpjnVtX1e1O8fnvSnLTBd/3INdIcmmS91fV5yb5jjN4rZ9L8pD17Kyq+vSqumdVXeNQkgIwBJupPjwoq3rIG5K8N6sLpj/rxJ2ttRMbpXOz6ux/itbah7O6kPvuSf4xq7c4flBr7a/XD3lqkg9n9cvJc7K6GPuEa2b1i8F7s6oGXpTkJw/h53r5Ou/frl/3QzmzWt5J1wk4da21pyR5RFb/bMN7sjo2vzvJi9cPeW5Wbxzx1qz+IuOFC77HR7P6S53bJLkgq9n0rCTXOsWXeFKSH1pfC/Wo0/3++3hUVmf4L8tq5p32z3RCa+21WV039Yys5tGbknzzmUcEYCTV2pm2KAAAAHaPM1MAAAAL2EwBAAAsYDMFAACwgM0UAADAAjZTAAAAC5x1sjur6rTe6u+nvvozziwNcCwe+rKLl/5Dpd0wn2BOM8wnYHecdDM1+y8fd7nTHQ/ldV75mvMP5XU4Oq99xD0P7bVu+5TfPLTXYjnz6dSYT/0znwDGpeYHAACwgM0UAADAAjZTAAAAC9hMAQAALGAzBQAAsMBJ381vVx307leH9e5a9OOgd746zHfXgsNkPu0O8wmgf85MAQAALGAzBQAAsIDNFAAAwAI2UwAAAAvYTAEAACxgMwUAALCAzRQAAMACNlMAAAAL2EwBAAAsYDMFAACwgM0UAADAAjZTAAAAC9hMAQAALGAzBQAAsIDNFAAAwAI2UwAALFJVz6yqH974/Duq6l1V9f6qOqeq/m1VvXH9+b23mfVMVNVdq+rFR/w9vrKq/v4UH/uYqjrvKPMsUVXXr6q/qqqrbDvLcbGZAgDgU1TVW6vqg1V1WVVdUlWvqaqHVNXHf39srT2ktfa49eOvlOQpSe7aWrt6a+2iJI9N8oz150e6GTliT0zyYyc+qaq23jSetfG1s6rq3VXVtpLwEzlOeVN2CN/rrVX1VSc+b629K8n/TfJtx/H9e3DW5T9k99zlTnfcdgSOyWsfcc9tR4DTYj7tDvOJTtyrtfbbVXWtJF+R5GlJbp/kP+3z2OsnuWqS12987UZ7Pj9lVXVWa+0jS557mKrqdkmu1Vr7gz13XZLk7kleuv78Hknem+S6xxivR89L8rNJfnrbQY6DM1MAAJxUa+19rbWXJLlvkm+qqlsnSVU9u6oeX1W3TPI364dfUlW/U1VvTnLTJC9d1/yuUlXXqqqfr6p3VtXb18+94vq1vrmqfr+qnlpVFyd5zPrr/3ldHXtvVb28qm50Itf6DNFD1lXC91bVz1RVbdz/revnXlZVb6iqL1p//dyq+rWqek9VXVBV33uSH//uSX5vn68/N8mDNj5/UJJf3HzA+vu8pKourqo3VdW3btx3tfX6vbeq3pDkdvs891QzHqiqfreqHrde28uq6hVVdZ31fTder+G3VdU71n8uj9x47rOr6vEbn3/8rFdVPTfJDfOJP9/vXz/sD5PcdPPPaWY2UwAAnJLW2h8l+fskX7bn63+b5FbrT89urf371trNkvxdVme3rt5a++ckz0nykSQ3T/KFSe6a5MEbL3X7JG9Jcr0kT1hfZ/WDSb4+qzM+r0ryS3tifU1WG5EvSHKfJHdLkqr6j1ltyB6U5JpJvjbJReua4kuT/FmSz05y5yQPq6q7HfBjf34+sVHc9OIkX15VZ1fV2es1+d97HvNL6/U6N8l/SPLEqrrz+r7/luRm64+7JfmmE09akPHy3D+rs4nXS3LlJI/ac/+/S3KLrP48/stmde8grbUH5pP/fH9i/fWPJHlTVn8e07OZAgDgdLwjybVP90lVdf2szvI8rLX2gdbau5M8Nck3bL52a+2nW2sfaa19MMm3J3lSa+2v1r+kPzHJbfac9fix1tolrbW/y+p6ndusv/7gJD/RWvvjtvKm1trbstp4Xbe19tjW2odba29J8nN7cmw6O8ll+3z9Q1lteO67fu5L1l878fPeIMmXJnl0a+1DrbXXJXlWkgeuH3KfJE9orV3cWrswydM3Xvt0M16eX2it/e16TX85n1ijE350/WfyF0l+Icn9Fn6fEy7Lat2m55opAABOx2cnuXjB826U5EpJ3rnRxLtCkgs3HnPhPs95WlU9eeNrtc7wtvXn/7Bx3z8lufr69g2SvPmAHOdW1SUbX7tiVme99vPeJNc44L5fTPKkdaZH77nv3CQXt9Y2N2JvS3Lbjfsv3HPf0oyX56A1OmFvjs9f+H1OuEZW15RNz2YKAIBTsn4zhs9O8uoFT78wyT8nuc5J3lhi7zvhXZjV2ZvnLfx+Nzvg6xe01m5xiq/z50luecB9r0ryWVnlfvWe7/eOJNeuqmtsbKhumOTt69vvzGrD9/qN+5ZmPFM3SPLXGznesb79gSSftvG4z9zzvE9558L1OxzePKuK4vR2ejP1ytecv+0IHJPbPuU3tx0BTov5tDvMJ0ZQVddM8uVZvZvfees62Glprb2zql6R5Mm1+rep3p/kJkk+p7W23xs8JMkzkzyuql7XWnv9+l0F79pa+5VT+JbPSvKUqnp1kj/NaqPzL0n+KMmlVfXorKp1H07yeUmu1lr7431e57eSvOCAn6lV1b02bm/ed2FVvSbJk6rqUVltyL4lyQPWD/nlJD9QVX+Y5NOTfM/GS59uxjP1w+s3x7hJVtdWncj4uiSPXL8JxZWTPGzP896V1ZuMbPqSJG9dVyqn55opAAAO8tKquiyrMyX/Nat/R2q/t0U/VQ/K6pfyN2RVn/vVrM7s7Ku19qIkP57kBVV1aZK/zOq6q8u13nA9Icnzs7qG58VJrt1a+2iSe2V13dAFSf4xq43XtQ54nT9N8r6quv0B97++tXbQ27/fL8mNszrT86Ik/6219sr1fT+aVaXugiSvyOrdAU+85mllPAS/l9WbRvyfJD/ZWnvF+uvPzeoM01vXGV+453lPSvJDtfp3yE68qcU3ZrUJ3gnV2sH/rtjT7n7trf6jY8DReOjLLq7Lf1TfzCeY0wzziflU1V2TfGdr7d7bznKYqurGWW3WrnQY/6ZXVV0vq43ZF7bWPnR5j5/BTtf8AADg8qzP1Lzich+449bv0Ph5285xnNT8AAAAFjhpze+iiy5So4EJnXPOOcPXaMwnmNMM8wnYHc5MAQAALOCaKQCAQ3KTpz3ZWXM4Im++7+m9SeDNXviQQ/veFzz0kfueNXdmCgAAYAGbKQAAgAVspgAAABawmQIAAFjAG1AAAADdO8w3lDgszkwBAAAsYDMFAACwgM0UAADAAjZTAAAAC9hMAQAALODd/AAAgO7d/OF/cFqPf9NT73BEST7BmSkAAIAFTnpm6rwH3OK4csDOee0j7rnv12/7lN888u/90JddfOTf46iZTzCnGeYTsDucmQIAAFjgpJuph60/zvR2DuH2YeToMVMvOXrMNFOO/TI9b/2x9/ZxZJpBT3++jsGjzdFjpply9JYJYCTOTAEAACxQrbUD73za3a998J3AGdnyNVN15N/kiJlPMKfR59NNnvZkswmOyDbfze+Chz5y39mk5rflTL3k6DHTTDn2y6Tmd2Z6+vN1DB5tjh4zzZSjt0wAI1HzAwAAWEDND7ZEze/MmE9wdMyn5dT84Oio+Z3BbVWRo83RY6aZcuyXSc3vzPT05+sYPNocPWaaKcd+mcwngFOj5gcAALCAmh/soNFrNIn5BEdJzW85NT84Omp+Z3BbVeRoc/SYaaYcvWWaQS9r6Ri0NqPn2C+Tmh/AqVHzAwAAWEDND3bQ6DWaxHyCo6Tmt5yaHxwdNb8zuK0qcrQ5esw0U47eMs2gl7V0DFqb0XPsl0nND+DUqPkBAAAsoOYHO2j0Gk1iPsFRUvNbTs0Pjo6a3xncVhU52hw9ZpopR2+ZZtDLWjoGrc3oOfbLpOYHcGrU/AAAABZQ84MdNHqNJjGf4Cip+S2n5gdHZ7iaHwAAAPtzzdSWM/WSo8dMM+XoLdMMellLx6C1GT3HfplcMwVwatT8YAeNXqNJzCc4Smp+y6n5wdHpseZ31qF9BwAAgCNymJujw6Lmt+VMveToMdNMOXrLNINe1tIxaG1Gz7FfJjU/gFOj5gc7aPQaTWI+wVFS81tOzQ/m5N38AAAADpGa35Yz9ZKjx0wz5egt0wx6WUvHoLUZPcd+mdT8AE6Nmh/soNFrNIn5BEdJzW85NT+Yk5ofAADAIVLz23KmXnL0mGmmHL1lmkEva+kYtDaj59gvk5ofwKlR84MdNHqNJjGfYFajzyc1P5iTmh8AAMAhUvPbcqZecvSYaaYcvWWaQS9r6Ri0NqPn6C0TwEjU/GAHjV6jScwnmNXo80nND+ak5gcAAHCI1Py2nKmXHD1mmilHb5lm0MtaOgatzeg5essEMBI1P9hBo9doEvMJZjX6fFLzgzmp+QEAABwiNb8tZ+olR4+ZZsrRW6YZ9LKWjkFrM3qO3jIBjMSZKQAAgAVcMwU7aPRrEhLzCWY1+nxyzRTMadE1U72c8lcVsTaj5+gt0wx6WUvHoLUZPUdvmQBGouYHAACwgJof7KDRazSJ+QSzGn0+qfnBnNT8Os3US44eM82Uo7dMM+hlLR2D1mb0HL1lAhiJmh8AAMACan6wg0av0STmE8xq9Pmk5gdzUvPrNFMvOXrMNFOO3jLNoJe1dAxam9Fz9JYJYCRqfgAAAAuo+cEOGr1Gk5hPMKvR55OaH8xJza/TTL3k6DHTTDl6yzSDXtbSMWhtRs/RWyaAkaj5AQAALKDmBzto9BpNYj7BrEafT2p+MCc1v04z9ZKjx0wz5egt0wx6WUvHoLUZPUdvmQBGouYHAACwgJof7KDRazSJ+QSzGn0+qfnBnBbV/AAAANifa6a2nKmXHD1mmilHb5lm0MtaOgatzeg5essEMBI1P9hBo9doEvMJZjX6fFLzgzmp+QEAABwiNb8tZ+olR4+ZZsrRW6YZ9LKWjkFrM3qO3jIBjETND3bQ6DWaxHyCWY0+n9T8YE5qfgAAAIdIzW/LmXrJ0WOmmXL0lmkGvaylY9DajJ6jt0wAI1Hzgx00eo0mMZ9gVqPPJzU/mJOaHwAAwCFS89typl5y9Jhpphy9ZZpBL2vpGLQ2o+foLRPASNT8YAeNXqNJzCeY1ejzSc0P5qTmBwAAcIjU/LacqZccPWaaKUdvmWbQy1o6Bq3N6Dl6ywQwEjU/2EGj12gS8wlmNfp8UvODOan5AQAAHCI1vy1n6iVHj5lmytFbphn0spaOQWszeo7eMgGMxJkpAACABVwzBTto9GsSEvMJZjX6fHLNFMxp0TVTvZzyVxWxNqPn6C3TDHpZS8egtRk9R2+ZAEai5gcAALCAmh/soNFrNIn5BLMafT6p+cGc1Pw6zdRLjh4zzZSjt0wz6GUtHYPWZvQcvWUCGImaHwAAwAJqfrCDRq/RJOYTzGr0+aTmB3NS8+s0Uy85esw0U47eMs2gl7V0DFqb0XP0lglgJGp+AAAAC6j5wQ4avUaTmE8wq9Hnk5ofzEnNr9NMveToMdNMOXrLNINe1tIxaG1Gz9FbJtjPzR/+Bwd+wDap+QEAACyg5gc7aPQaTWI+waxGn09qfkfjZGeg3vTUOxxjEnaVml+nmXrJ0WOmmXL0lmkGvaylY9DajJ6jt0wAI1HzAwAAWEDND3bQ6DWaxHyCWY0+n9T8joaaH9um5tdppl5y9Jhpphy9ZZpBL2vpGLQ2o+foLRPASNT8AAAAFlDzgx00eo0mMZ9gVqPPJzW/o6Hmx7YtqvkBAACwP9dMbTlTLzl6zDRTjt4yzaCXtXQMWpvRc/SWCWAkan6wg0av0STmE8xq9Pmk5nc01PzYtoNqfmcddxAAAJY72cZiPzNsNmb4GZiTmt+WM/WSo8dMM+XoLdMMellLx6C1GT1Hb5kARqLmBzto9BpNYj7BrEafT8dR89vFM1Owbd7NDwAA4BCp+W05Uy85esw0U47eMs2gl7V0DFqb0XP0lglgJGp+sINGr9Ek5tNRecB5b/ykz897wC22lIRdNfp8UvODOXk3PwAAuvbm+z5z36/f7IUPOeYkcGrU/LacqZccPWaaKUdvmWbQy1rOdgxed/2x7Rw9rs2sOXrLBDASNT/YQaPXaBLz6aio+bFto88nNb8z48wUvfJufgAAAIdIzW/LmXrJ0WOmmXL0lmkGvazlbMegmt/u5egtE8BI1PxgB41eo0nMJ5jV6PNJze/MqPnRK+/mBwAwgZE2RzA7Nb8tZ+olR4+ZZsrRW6YZ9LKWjkFrM3qO3jIBjMQbUAAAACzgminYQaNfk5CYTzCr0efTcVwzNTPXTNGrRW+N3sspf1URazN6jt4yzaCXtXQMWpvRc/SWCWAkan4AAAALqPnBDhq9RpOYTzCr0eeTmt+ZUfOjV2p+nWbqJUePmWbK0VumGRzHWj78TnfMXdYfm7cdg+bTTDl6y8Ruu9kLH7LvB/RKzQ8AAGABNb8j8oDz3vjx2+c94BZbTAKfavQaTXI88+kud7rjvl9/5WvOP+pvDTtr9Pmk5gdzUvPbQqbrrj+2naPHtdm1HL1lmsFxrOWtD7jtGDSfZsrRWyaAkaj5AQAALKDmd0TU/OjZ6DWaRM0PZjX6fFLzgzmp+W0hk5qfHL1mmsFxrKWaXx85esw0U47eMgGMRM0PAABgATU/2EGj12gSNT+Y1ejzSc0P5qTm12mmXnL0mGmmHL1lmsFxrKWaXx85esw0U47eMgGMRM0PAABgATU/2EGj12gSNT+Y1ejzSc0P5nRQze+s4w4CAMByb77vM0/r8Td74UOOKAmw9WumHn6nO+Yu64/N23r3rkmYKUdvmWZwHGvpmqk+cvSYaaYcvWUCGMnWa35qNHD8Rq/RJOYTzGr0+XQcNT9npuD4LXo3PwAAAPa39ZqfGk0fOXrMNFOO3jLNwHzq/78zazNGjt4yAYxEzQ920Og1msR8glmNPp/U/GBOan4AAACHSM3vDG7PlKPHTDPl6C3TDMyn/v87szZj5OgtE8BI1PxgB41eo0nMJ5jV6PNJzQ/mpOYHAABwiNT8zuD2TDl6zDRTjt4yzcB86v+/M2szRo7eMgGMRM0PdtDoNZrEfIJZjT6f1PxgTgfV/M467iAAACxncwT9UPM7g9sz5egx00w5ess0g+NYy1u95vy8cv2xedsxaD7NlKO3TAAjUfODHTR6jSY5nvkEHL/R59Nx1PyA4+fd/AAAAA6Rmt8Z3J4pR4+ZZsrRW6YZ9LKWjkFrM3qO3jIBjETND3bQ6DWaRM0PZjX6fFLzgzmp+QEAABwiNb8zuD1Tjh4zzZSjt0wz6GUtHYPWZvQcvWUCGIkzUwAAAAu4Zgp20OjXJCSumYJZjT6fXDMFc1p0zdRxnPJX8+sjR4+ZZsrRW6YZ9LKWjkFrM3qO3jIBjETNDwAAYAE1P9hBo9doEjU/mNXo80nND+bUbc3vVq85P69cf2zeVhVRo5kpR2+ZZtDLWjoGrc3oOXrLBDASNT8AAIAFtl7zA47f6DWaxHyCWY0+n9T8YE7d1vxURfrI0WOmmXL0lmkGvaylY9DajJ6jt0wAI1HzAwAAWEDND3bQ6DWaxHyCWY0+n9T8YE5qfp1m6iVHj5lmytFbphn0spaOQWszeo7eMgGMRM0PAABgATU/2EGj12gS8wlmNfp8UvODOan5dZqplxw9ZpopR2+ZZtDLWjoGrc3oOXrLBDASNT8AAIAF1PxgB41eo0nMJ5jV6PNJzQ/mtKjmBwAAwP5cM7XlTL3k6DHTTDl6yzSDXtbSMWhtRs/RWyaAkaj5wQ4avUaTmE8wq9Hnk5ofzEnNDwAA4BCp+W05Uy85esw0U47eMs2gl7V0DFqb0XP0lglgJGp+sINGr9Ek5hPMavT5pOYHc1LzAwAAOERqflvO1EuOHjPNlKO3TDPoZS0dg9Zm9By9ZQIYiZof7KDRazSJ+QSzGn0+qfnBnNT8AAAADpGa35Yz9ZKjx0wz5egt0wx6WUvHoLUZPUdvmQBGouYHO2j0Gk1iPsGsRp9Pan4wJzU/AACAQ6Tmt+VMveToMdNMOXrLNINe1tIxaG1Gz9FbJoCRqPnBDhq9RpOYTzCr0eeTmh/MSc0PAADgEKn5bTlTLzl6zDRTjt4yzaCXtXQMWpvRc/SWCWAkzkwBAAAs4Jop2EGjX5OQmE8wq9Hnk2umYE6Lrpnq5ZS/qoi1GT1Hb5lm0MtaOgatzeg5essEMBI1PwAAgAXU/GAHjV6jScwnmNXo80nND+ak5tdppl5y9Jhpphy9ZZpBL2vpGLQ2o+foLRPASNT8AAAAFlDzgx00eo0mMZ9gVqPPJzU/mJOaX6eZesnRY6aZcvSWaQa9rKVj0NqMnqO3TAAjUfMDAABYQM0PdtDoNZrEfIJZjT6f1PxgTmp+nWbqJUePmWbK0VumGfSylo5BazN6jt4yAYxEzQ8AAGABNT/YQaPXaBLzCWY1+nxS84M5qfl1mqmXHD1mmilHb5lm0MtaOgatzeg5essEMBI1PwAAgAXU/KJi5ZkAAARXSURBVGAHjV6jScwnmNXo80nND+ak5tdppl5y9Jhpphy9ZZpBL2vpGLQ2o+foLRPASNT8AAAAFlDzgx00eo0mMZ9gVqPPJzU/mNOimh8AAAD7c83UljP1kqPHTDPl6C3TDHpZS8egtRk9R2+ZAEai5gc7aPQaTWI+waxGn09qfjAnNT8AAIBDpOa35Uy95Ogx00w5ess0g17W0jFobUbP0VsmgJGo+cEOGr1Gk5hPMKvR55OaH8xJzQ8AAOAQqfltOVMvOXrMNFOO3jLNoJe1dAxam9Fz9JYJYCRqfrCDRq/RJOYTzGr0+aTmB3NS8wMAADhEan5bztRLjh4zzZSjt0wz6GUtHYPWZvQcvWUCGImaH+yg0Ws0ifkEsxp9Pqn5wZzU/AAAAA6Rmt+WM/WSo8dMM+XoLdMMellLx6C1GT1Hb5kARqLmBzto9BpNYj7BrEafT2p+MCc1PwAAgEOk5rflTL3k6DHTTDl6yzSDXtbSMWhtRs/RWyaAkTgzBQAAsIBrpmAHjX5NQmI+waxGn0+umYI5LbpmqpdT/qoi1mb0HL1lmkEva+kYtDaj5+gtE8BI1PwAAAAWUPODHTR6jSYxn2BWo88nNT+Yk5pfp5l6ydFjpply9JZpBr2spWPQ2oyeo7dMACNR8wMAAFhAzQ920Og1msR8glmNPp/U/GBOan6dZuolR4+ZZsrRW6YZ9LKWjkFrM3qO3jIBjETNDwAAYAE1P9hBo9doEvMJZjX6fFLzgzmp+XWaqZccPWaaKUdvmWbQy1o6Bq3N6Dl6ywQwEjU/AACABdT8YAeNXqNJzCeY1ejzSc0P5nRQze+skz3pxGn3nzrD2yde60xuH0aOHjP1kqPHTDPl6C3TDHpZS8egtRk9R2+ZRnfQL1zAnNT8AAAAFlDzgx00eo0mMZ9gVjPMJ2B3ODMFAACwgLdG33KmXnL0mGmmHL1lmkEva+kYtDaj5+gtE8BI1PxgB81QozGfYE4zzCdgd6j5AQAALKDmt+VMveToMdNMOXrLNINe1tIxaG1Gz9FbJoCRqPnBDpqhRmM+wZxmmE/A7lDzAwAAWEDNb8uZesnRY6aZcvSWaQa9rKVj0NqMnqO3TAAjUfODHTRDjcZ8gjnNMJ+A3aHmBwAAsICa35Yz9ZKjx0wz5egt0wx6WUvHoLUZPUdvmQBGctKa30UXXaRGAxM655xzhq/RmE8wpxnmE7A71PwAAAAWsJkCAABYwGYKAABgAZspAACABWymAAAAFrCZAgAAWMBmCgAAYAGbKQAAgAVspgAAABawmQIAAFjAZgoAAGABmykAAIAFbKYAAAAWsJkCAABYwGYKAABgAZspAACABWymAAAAFrCZAgAAWMBmCgAAYAGbKQAAgAVspgAAABawmQIAAFjAZgoAAGABmykAAIAFqrW27QwAAADDcWYKAABgAZspAACABWymAAAAFrCZAgAAWMBmCgAAYAGbKQAAgAX+P/8d4IlAXpMtAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "next_observation, _,_,_ = env.step(np.random.choice(n_actions))\n", "diff = pong_change(observation, next_observation)\n", "\n", "f, ax = plt.subplots(1, 3, figsize=(15,15))\n", "for a in ax:\n", " a.grid(False)\n", " a.axis(\"off\")\n", "ax[0].imshow(observation); ax[0].set_title('Previous Frame');\n", "ax[1].imshow(next_observation); ax[1].set_title('Current Frame');\n", "ax[2].imshow(np.squeeze(diff)); ax[2].set_title('Difference (Model Input)');" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "What do you notice? How and why might these pre-processing changes be important for training our RL algorithm? How and why might consideration of the difference between frames be important for training and performance?" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Rollout function\n", "\n", "We're now set up to define our key action algorithm for the game of Pong, which will ultimately be used to train our Pong agent. This function can be thought of as a \"rollout\", where the agent will 1) make an observation of the environment, 2) select an action based on its state in the environment, 3) execute a policy based on that action, resulting in some reward and a change to the environment, and 4) finally add memory of that action-reward to its `Memory` buffer. We will define this algorithm in the `collect_rollout` function below, and use it soon within a training block.\n", "\n", "Earlier you visually inspected the raw environment frames, the pre-processed frames, and the difference between previous and current frames. As you may have gathered, in a dynamic game like Pong, it can actually be helpful to consider the difference between two consecutive observations. This gives us information about the movement between frames -- how the game is changing. We will do this using the `pong_change` function we explored above (which also pre-processes frames for us).\n", "\n", "We will use differences between frames as the input on which actions will be selected. These observation changes will be forward propagated through our Pong agent, the CNN network model, which will then predict the next action to take based on this observation. The raw reward will be computed. The observation, action, and reward will be recorded into memory. This will loop until a particular game ends -- the rollout is completed.\n", "\n", "For now, we will define `collect_rollout` such that a batch of observations (i.e., from a batch of agent-environment worlds) can be processed serially (i.e., one at a time, in sequence). We will later utilize a parallelized version of this function that will parallelize batch processing to help speed up training! Let's get to it." ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [], "source": [ "### Rollout function ###\n", "\n", "# Key steps for agent's operation in the environment, until completion of a rollout.\n", "# An observation is drawn; the agent (controlled by model) selects an action;\n", "# the agent executes that action in the environment and collects rewards;\n", "# information is added to memory.\n", "# This is repeated until the completion of the rollout -- the Pong game ends.\n", "# Processes multiple batches serially.\n", "#\n", "# Arguments:\n", "# batch_size: number of batches, to be processed serially\n", "# env: environment\n", "# model: Pong agent model\n", "# choose_action: choose_action function\n", "# Returns:\n", "# memories: array of Memory buffers, of length batch_size, corresponding to the\n", "# episode executions from the rollout\n", "def collect_rollout(batch_size, env, model, choose_action):\n", "\n", " # Holder array for the Memory buffers\n", " memories = []\n", "\n", " # Process batches serially by iterating through them\n", " for b in range(batch_size):\n", "\n", " # Instantiate Memory buffer, restart the environment\n", " memory = Memory()\n", " next_observation = env.reset()\n", " previous_frame = next_observation\n", " done = False # tracks whether the episode (game) is done or not\n", "\n", " while not done:\n", "\n", " current_frame = next_observation\n", "\n", " # determine the observation change.\n", " frame_diff = pong_change(previous_frame, current_frame)\n", "\n", " # choose an action for the pong model, using the frame difference, and evaluate\n", " action = choose_action(model, frame_diff)\n", "\n", " # Take the chosen action\n", " next_observation, reward, done, info = env.step(action)\n", "\n", " # save the observed frame difference, the action that was taken, \n", " # and the resulting reward!\n", " memory.add_to_memory(next_observation, action, reward)\n", "\n", " previous_frame = current_frame\n", " \n", " # Add the memory from this batch to the array of all Memory buffers\n", " memories.append(memory)\n", " \n", " return memories" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "To get a sense of what is encapsulated by `collect_rollout`, we will instantiate an *untrained* Pong model, run a single rollout using this model, save the memory, and play back the observations the model sees. Note that these will be frame *differences*." ] }, { "cell_type": "code", "execution_count": 19, "metadata": {}, "outputs": [], "source": [ "def save_video_of_memory(memory, filename, size=(512,512)):\n", " import skvideo.io\n", "\n", " output_video = skvideo.io.FFmpegWriter(filename)\n", "\n", " for observation in memory.observations:\n", " output_video.writeFrame(cv2.resize(255*observation, size))\n", " \n", " output_video.close()\n", " return filename" ] }, { "cell_type": "code", "execution_count": 22, "metadata": {}, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "### Rollout with untrained Pong model ###\n", "\n", "# Model\n", "test_model = create_pong_model()\n", "\n", "# Rollout with single batch\n", "single_batch_size = 1\n", "memories = collect_rollout(single_batch_size, env, test_model, choose_action)\n", "rollout_video = save_video_of_memory(memories[0], \"video/Pong-Random-Agent.mp4\")\n", "\n", "show_video(\"Pong-Random-Agent\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 3.9 Training Pong\n", "\n", "We're now all set up to start training our RL algorithm and agent for the game of Pong! We've already defined the following:\n", "\n", "1. Loss function, `compute_loss`, and backpropagation step, `train_step`. Our loss function employs policy gradient learning. `train_step` executes a single forward pass and backpropagation gradient update.\n", "2. RL agent algorithm: `collect_rollout`. Serially processes batches of episodes, executing actions in the environment, collecting rewards, and saving these to `Memory`.\n", "\n", "We will use these functions to train the Pong agent.\n", "\n", "In the training block, episodes will be executed by agents in the environment via the RL algorithm defined in the `collect_rollout` function. Since RL agents start off with literally zero knowledge of their environment, it can often take a long time to train them and achieve stable behavior. To alleviate this, we have implemented a parallelized version of the RL algorithm, `parallelized_collect_rollout`, which you can use to accelerate training across multiple parallel batches.\n", "\n", "For training, information in the `Memory` buffer from all these batches will be aggregated (after all episodes, i.e., games, end). Discounted rewards will be computed, and this information will be used to execute a training step. Memory will be cleared, and we will do it all over again!\n", "\n", "Let's run the code block to train our Pong agent. Note that, even with parallelization, completing training and getting stable behavior will take quite a bit of time (estimated at least a couple of hours). We will again visualize the evolution of the total reward as a function of training to get a sense of how the agent is learning." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "> Note: If you want to increase the model performance, try to increase the number of iteration" ] }, { "cell_type": "code", "execution_count": 26, "metadata": {}, "outputs": [], "source": [ "### Hyperparameters and setup for training ###\n", "# Rerun this cell if you want to re-initialize the training process\n", "# (i.e., create new model, reset loss, etc)\n", "\n", "# Hyperparameters\n", "learning_rate = 1e-3\n", "MAX_ITERS = 251 # increase the maximum to train longer\n", "batch_size = 5 # number of batches to run\n", "\n", "# Model, optimizer\n", "pong_model = create_pong_model()\n", "optimizer = tf.keras.optimizers.Adam(learning_rate)\n", "iteration = 0 # counter for training steps\n", "\n", "# Plotting\n", "smoothed_reward = LossHistory(smoothing_factor=0.9)\n", "smoothed_reward.append(0) # start the reward at zero for baseline comparison\n", "plotter = PeriodicPlotter(sec=15, xlabel='Iterations', ylabel='Win Percentage (%)')\n", "\n", "# Batches and environment\n", "# To parallelize batches, we need to make multiple copies of the environment.\n", "envs = [create_pong_env() for _ in range(batch_size)] # For parallelization" ] }, { "cell_type": "code", "execution_count": 27, "metadata": {}, "outputs": [], "source": [ "def parallelized_collect_rollout(batch_size, envs, model, choose_action):\n", "\n", " assert len(envs) == batch_size, \"Number of parallel environments must be equal to the batch size.\"\n", "\n", " memories = [Memory() for _ in range(batch_size)]\n", " next_observations = [single_env.reset() for single_env in envs]\n", " previous_frames = [obs for obs in next_observations]\n", " done = [False] * batch_size\n", " rewards = [0] * batch_size\n", "\n", " while True:\n", "\n", " current_frames = [obs for obs in next_observations]\n", " diff_frames = [pong_change(prev, curr) for (prev, curr) in zip(previous_frames, current_frames)]\n", "\n", " diff_frames_not_done = [diff_frames[b] for b in range(batch_size) if not done[b]]\n", " actions_not_done = choose_action(model, np.array(diff_frames_not_done), single=False)\n", "\n", " actions = [None] * batch_size\n", " ind_not_done = 0\n", " for b in range(batch_size):\n", " if not done[b]:\n", " actions[b] = actions_not_done[ind_not_done]\n", " ind_not_done += 1\n", "\n", " for b in range(batch_size):\n", " if done[b]:\n", " continue\n", " next_observations[b], rewards[b], done[b], info = envs[b].step(actions[b])\n", " previous_frames[b] = current_frames[b]\n", " memories[b].add_to_memory(diff_frames[b], actions[b], rewards[b])\n", "\n", " if all(done):\n", " break\n", "\n", " return memories" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEGCAYAAABiq/5QAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO3dd3Rc1bXH8e9WteQmF7nKHVeMOy6YZqrpJST0QELikIQXCCmUNF4S8iAhpEEITiD0FkpMKAYDptnGtox77w3Zkotk2bLazH5/zJjIRmUkazSS5vdZa5Zm7tyyrzzec3TuufuYuyMiIvEjIdYBiIhIw1LiFxGJM0r8IiJxRolfRCTOKPGLiMSZpFgHEImOHTt67969Yx2GiEiTsmDBgl3unnnk8iaR+Hv37k12dnaswxARaVLMbHNly9XVIyISZ5T4RUTijBK/iEicUeIXEYkzSvwiInFGiV9EJM5ELfGbWQszm2dmi81suZn9b3j5Y2a20cwWhR8johWDiIh8UTRb/CXAae4+HBgBTDaz8eH3fuTuI8KPRVGMQUSkSSosLuPXr61g464D9b7vqCV+D9kffpkcfqj4v4hIBD5cs4t/fLyRvMKSet93VPv4zSzRzBYBucAMd58bfutuM1tiZn8ws9RoxiAi0hTNWLGDdunJjO7Vrt73HdXE7+4Bdx8BZAFjzWwocAcwCDgeaA/cVtm2ZjbFzLLNLDsvLy+aYYqINCplgSDvrcrltEGdSUywet9/g4zqcfd84H1gsrvnhLuBSoB/AmOr2Gaqu49x9zGZmV+oMSQi0uwcKCknEHTmrN/NvuJyzhzSKSrHiVqRNjPLBMrcPd/M0oAzgHvNrKu755iZARcDy6IVg4hIU1AeCPLAzHU88N46khMTKC4P0KZFEif1j06jN5rVObsCj5tZIqG/LF5w99fM7L3wl4IBi4AboxiDiEij9/iczfzxnbWcN6wrma1SaZuWzKWjutMyNTopOmqJ392XACMrWX5atI4pItLUBILO47M3cXzvdjx41agGOabu3BURiaGZq3LZsqeI60/o02DHVOIXEYmhx+dsokubFpx1bOcGO6YSv4hIjKzLLeSjtbu4dkIvkhMbLh0r8YuIxMjjszeTkpTAFcf3aNDjKvGLiMRAcVmAlz/dxgXDutGhVcMWMFDiFxGJgQ/X5HGgNMDFI7s1+LGV+EVEYmD68h20TUtmfN8ODX5sJX4RkQZWWh7knRU7OWNw5wa9qHuIEr+ISANau7OQy/42m33F5Zw3rEtMYlDiFxFpQPdOX83m3UX87rJhTBoYnSJsNVHiFxFpIPlFpXywJpfLj+/Bl8f0IFSrsuEp8YuINJDXl+ZQFnAuHN7wI3kqUuIXEWkAxWUBnv5kC/0yW3JstzYxjUWJX0QkyoJB56ZnPmXljn386OyBMeviOUSJX0QkyuZt2sM7K3O585zBTB7aNdbhKPGLiETb9GU7SE1K4KpxPWMdCqDELyISVcGgM33ZDk4ekBm1GbVqS4lfRCSKFm3LZ8e+Ys4ZGpubtSqjxC8iEkXvrcwlMcE4fVDDTbRSk6glfjNrYWbzzGyxmS03s/8NL+9jZnPNbK2ZPW9mKdGKQUQk1j5at4vhWW1pm54c61A+F80WfwlwmrsPB0YAk81sPHAv8Ad37w/sBW6IYgwiIjFTUFTG0m35nNg/M9ahHCZqid9D9odfJocfDpwGvBhe/jhwcbRiEBGJpdnrdxF0OKl/x1iHcpio9vGbWaKZLQJygRnAeiDf3cvDq2wDukczBhGRWHB3XluSQ6vUJEb0yIh1OIeJauJ394C7jwCygLHA4MpWq2xbM5tiZtlmlp2XlxfNMEVE6lVpeZDbX1rK60tzuGZ8w06kHokGicbd84H3gfFAhpkdGsyaBXxWxTZT3X2Mu4/JzGxc/WMiIlU5UFLO1f/4hOezt/K9047htskDYx3SF0RzVE+mmWWEn6cBZwArgZnAZeHVrgOmRSsGEZGG9s7KnczftJffXjaMW8+KfV2eykTzNrKuwONmlkjoC+YFd3/NzFYAz5nZr4GFwCNRjEFEpEEt2VZAi+QELh3ZeC9fRi3xu/sSYGQlyzcQ6u8XEWl2lmzL59hubUlqZP36FTXeyEREmphA0Fm2fR/HdW8b61CqpcQvIlJP1uft52BZgGFZSvwiInFh8dZ8AIZlNa5x+0dS4hcRqSez1u2idWoSfTu2jHUo1VLiFxGpB7n7inl9aQ5fGp1FQkLjG8JZUeOYFUBEpAnbvPsAf/9oA+VB5/oTesc6nBop8YuI1FHBwTLufGUpry/JAeC847rSu5F384ASv4hInT31yWZeX5LD907vz8R+HRjZs12sQ4qIEr+ISB2tyNlHj/Zp3HrmgFiHUiu6uCsiUkercvYxqEubWIdRa0r8IiJ1UFwWYOOuAwzu0jrWodSaEr+ISB2sy91P0GFQV7X4RUTiwsqcfQAMVItfRCQ+rN5RSGpSAr07NP7hm0dS4hcRqYPF2/IZ1KU1iY38Lt3KKPGLiNRS7r5isjfvZdKgTrEOpU6U+EVEamn68h24h+7UbYpqvIHLzDoBE4FuwEFgGZDt7sEoxyYi0ii9viSH/p1a0b9z07uwC9W0+M1skpm9BbwOnENoDt0hwE+BpWb2v2bW9MYxiYgchVcXf8bcjXu4YHi3WIdSZ9W1+M8FvunuW458w8ySgPOBM4GXKtvYzHoATwBdgCAw1d3/ZGZ3Ad8E8sKr3unub9T5DEREGsianYX88IXFjO3dnikn9411OHVWZeJ39x9V81458O8a9l0O/MDdPzWz1sACM5sRfu8P7n5fraMVEYmhN5bmUB4M8tdrRtEiOTHW4dRZxBd3zWy8mb1nZrPM7JKa1nf3HHf/NPy8EFgJdK97qCIisfXR2l0cl5VBx1apsQ7lqFTXx9/liEW3AhcCk4Ff1uYgZtYbGAnMDS+6ycyWmNmjZtY06piKSFzbV1zGoq35nHRMx1iHctSqa/H/zcx+ZmYtwq/zgauAy4F9kR7AzFoRug5wi7vvAx4C+gEjgBzg91VsN8XMss0sOy8vr7JVREQazJz1uwkEnZP6N+PE7+4XA4uA18zsWuAWQhdp04GLI9m5mSUTSvpPu/vL4f3udPdAeDjo34GxVRx/qruPcfcxmZmZtTknEZF69fKn2/j5tGW0Sk1qMpOtVKfaPn53/w9wNpABvAysdvc/u3uNTXAzM+ARYKW7319hecU7Hi4hdF+AiEijVFwW4I6Xl5LZOpVHrz+elKSmf99rdX38F5rZx8B7hJLzFcAlZvasmfWLYN8TgWuB08xsUfhxLvBbM1tqZkuAScD3j/40RESiY86G3ZSUB/nx2YMY26d9rMOpF9WN4/81MAFIA95w97HArWbWH7ib0BdBldz9Y6Cy6kUasy8ijV5JeYC1O/fzweo8WiQnNJukD9Un/gJCyT0NyD200N3XUkPSFxFp6u6fsYaHP9hAekoiE/p2aNLj9o9UXWfVJYQu5JYTGs0jIhIXygNBXlqwnZSkBIpKA5w6sGlW4axKdS3+Ynf/S3Ubm1krd99fzzGJiMTUR2t3sWt/CQ9cNZKgw9nHdo51SPWqusQ/zcwWAdOABe5+AMDM+hK6KPsVQsMxX4x6lCIiDWR/STkPvb+edunJnDWkS7MYxXOk6mr1nB4ehfMtYGL4DttyYDWhip3XufuOhglTRCT6ygJBLntoNmtz9/N/lxzXLJM+1FCPP1w1U6NwRCQuzFm/m1U7Cvn9l4fzpdFZsQ4naprn15mISB28uSyHlimJnDesac6sFSklfhERIBB03l6+k0mDOjWroZuVUeIXkbgxfdkO8otKD1v2/PwtLNqaz9yNu9l9oJRzhjbv1j5EMOcugJmdCPR393+aWSbQyt03Rjc0EZGjN3N1Ltv2HmRkjwxufGoBZx/bmetP6MN7q3Zy0Yju3PbSUob3yGBg51a0TElk0qDmXxQyksnWfwGMAQYC/wSSgacI1eIREWm0issC3PbiEnbtL+GMwaGx+G8t38nMVXmUBoI8P38rAIu35rPys31cOKIb6SkRtYebtEi6ei4hNAHLAQB3/wxomlPLi0hceenTbeQWluDA2yt2MmlgJgM7tyarfRqXjurOvuJyvjqhF6lJCZQGglw6Kj4mCYzkq63U3d3MHMDMWkY5JhGRo+buTP1wAyN6ZDCoS2uem7+Vy4/vwckDMklKSMAJTaoy+diulAWcuRt3M75Ph1iH3SAiSfwvmNnDQIaZfRP4OqE7dkVEGq2l2wvYvLuImyYdwykDM+naNo3TB3cmOfG/HR2XjAyN1f/VRcdSHnQSEiorKNz81Jj43f0+MzuT0HSLA4Gfu/uMqEcmInIU3ly2g6QE48whnclIT+HmM/pXuW5SYgJJzXsE52EiuooRTvRK9iLSJLg705ftYEK/DmSkp8Q6nEanxou7ZlZoZvuOeGw1s1fCBdtERBqVNTv3s3HXAc4+tkusQ2mUImnx3w98BjxDaEatK4AuhIq1PQqcGq3gRETqYvqyHZjBWc2snHJ9iWQ452R3f9jdC919n7tPBc519+eBpj/dvIg0O28uy2FMr3Z0at0i1qE0SpEk/qCZfcXMEsKPr1R4z6MVmIhIXWzadYBVOwrVzVONSBL/1cC1hObd3Rl+fo2ZpQE3VbWRmfUws5lmttLMlpvZzeHl7c1shpmtDf/UXw0iUi8+WpvHzc8tBGDyUCX+qkQynHMDcEEVb39czablwA/c/VMzaw0sMLMZwPXAu+5+j5ndDtwO3Fa7sEVEDldcFuC7T39K6xbJ/OrioWS1S491SI1WJLV6WgA3AMcCn3eYufvXq9vO3XOAnPDzQjNbCXQHLuK/F4QfB95HiV9EjtL0ZTvYV1zO364dzQn9OsY6nEYtkq6eJwmN4jkb+ADIAgprcxAz6w2MBOYCncNfCoe+HCqdvt7MpphZtpll5+Xl1eZwIhJHsjft4cIHPubP766lV4f0uCm7cDQiSfzHuPvPgAPu/jhwHnBcpAcws1bAS8At7r4v0u3cfaq7j3H3MZmZzb9MqojUzUPvr2fp9gI27DrAlWN7xk3ZhaMRyTj+svDPfDMbCuwAekeyczNLJpT0n3b3l8OLd5pZV3fPMbOuhC4ai4jUWk7BQWauzuXbp/Tj/GHdGNhFhYMjEUmLf2p45M1PgVeBFcC9NW1kZgY8Aqx09/srvPUqcF34+XXAtFpFLCIS9vz8rQQdrji+J0O6tSFRrf2IRNLif9fd9wIfAn0BzKxPBNtNJDT0c6mZLQovuxO4h1DFzxuALcCXax21iMS93ftLeOTjjZw2qBM9O2gET21EkvhfAkYdsexFYHR1G7n7x4RKPFTm9AiOKyJSpd/PWENRaYA7zx0U61CanCoTv5kNIjSEs62ZXVrhrTZUGNYpItLQ8gpLeH7+Vq4e15NjOqlfv7aqa/EPBM4HMjj8Bq5C4JvRDEpEpDqvLv6MQNC5dnyvWIfSJFWZ+N19GjDNzCa4+5wGjElEpFovf7qNYVlt6d9Zrf26iKSPf52Z3UloCOfn69d0566ISDS8sTSH5Z/t464LhsQ6lCYrksQ/DfgIeAcIRDccEZGqTVu0nZufW8SwrLZ8aXRWrMNpsiJJ/Onurlo6IhJzj8/exIDOrXjhWxNokRxHk+TWs0hu4HrNzM6NeiQiItXYua+YT7fkc+Hwbkr6RymSxH8zoeRfHJ5vt9DMIq65IyJSH95evgNQnf36EEk9fl02F5GYCgadlxdup19mS43brwc1tvgt5Boz+1n4dQ8zGxv90EREQu6fsYaFW/L5xkl9Yx1KsxBJV89fgQnAVeHX+4EHoxaRiEgFm3cf4IGZ6/jKmCyuOL5HrMNpFiIZ1TPO3UeZ2UIAd99rZilRjktEBICP1+0C4MZT+hEq+itHK6J6/GaWCDiAmWUCwahGJSJxbeueIn4+bRkZ6SmUBoJ0bduCPh1bxjqsZiOSxP9n4BWgk5ndDVxGqDa/iEi9O1ga4OIHZ7GnqBR3SE1K4LxhXdXar0c19vG7+9PAj4H/IzR5+sXu/q9oByYi8WnGyp3sPlDK364ZTdu0ZErKg5o8vZ5FMqpnPLDd3R909weAbWY2LvqhiUg8+vfC7XRt24IzB3fm6xP7kJRgTDxGE6jXp0hG9TxEaCTPIQfCy0RE6tXu/SV8sCaPC0d0IyHBuOm0Y5hx6yl0bZsW69CalUgSv7m7H3rh7kEiuzYgIlIrM1bsJBB0LhzeDYDEBNNF3SiIJPFvMLPvmVly+HEzsCHagYlI/HlvVS7d2rZgSNc2sQ6lWYsk8d8InABsB7YB44ApNW1kZo+aWa6ZLauw7C4z225mi8IPFX8TEQBKygN8vG4Xpw3upBE8UVZtl014/P7V7n5FHfb9GPAA8MQRy//g7vfVYX8i0gQ98vFG9hwo4UBJgAWb9/LXq0eR1S7tC8l97oY9FJUGOG1QpxhFGj+qTfzuHjCzi4A/1HbH7v6hmfWuY1wi0gxMW7SdX722AoDkRMMw7n59JQfLAuQfLOOmScfw4Zo8juvelqfmbiYtOZEJfTV0M9oiuUg7y8weAJ4nNKIHAHf/tI7HvMnMvgpkAz9w972VrWRmUwh3KfXs2bOOhxKRWNm46wA/fWUZo3pm8PjXQ3Udp364gb+8tw6AVqlJfPOJbBITjEDQSU9J5I9XjCAtRbX2o80qDNipfAWzmZUsdnc/rcadh1r8r7n70PDrzsAuQuUffgV0jWTu3jFjxnh2dnZNq4lII1FYXMYlf53N7v0lvHrTifRonw7A/pJyrn90HhcM78Zpgzrx0dpdnD+8K0u2FtC9XZpG8NQzM1vg7mO+sLymxH+UB+1NhcQf6XtHUuIXaToOlgb42mPzyN60lydvGMeEfrr5KlaqSvyR3Lnb2cweMbM3w6+HmNkNdQyia4WXlwDLqlpXRJqmH/5rMXM37uH3XxmupN9IRTKc8zHgLaBb+PUa4JaaNjKzZ4E5wEAz2xb+svitmS01syXAJOD7dYpaRBql2et38frSHL5/xgAuGtE91uFIFSK5uNvR3V8wszsA3L3czAI1beTuV1ay+JHaBigiTUP2pj3c9epyumekMeVkzZTVmEWS+A+YWQf+W49/PFAQ1ahEpEmZviyHG5/6lPSURP5y5UhaJGtkTmMWSeL/AfAq0M/MZgGZhGryi4gA8OQnm+nRPo3pN59My1SV8mrsavwXcvcFZnYKMBAwYLW7l0U9MhFptErLgzw6ayNXHN+D/SXlzF6/m1tOH6Ck30RU+a9kZv2B+4B+wFLgh+6+vaECE5HG652VO7nnzVV8smE3Xdu2AOBLo3Uxt6mo7uv5UUJ1dj4ELgT+AlzaEEGJSOM2c1UuCQbvr84D4Nrxvchqlx7jqCRS1SX+1u7+9/Dz35lZXUs0iEgzEgw676/J45zjujK0W1uy2qVxwfBuNW8ojUZ1ib+FmY0k1K8PkFbx9VHU6hGRJmxFzj7yCks4bWAnvjQ6K9bhSB1Ul/hzgPsrvN5R4bUDNdbqEZHm5dl5W/jjO2tIMDhlYGasw5E6qjLxu/ukhgxERBq3nIKD/GLaco7t3obfXjacjq1SYx2S1JHGXolIRB7+YANBd/5y5UhdyG3iIqnVIyJxbtveIp6dt4UvjcpS0m8GlPhF5HNb9xTxnacXMPbud5izfjcA7s4vpi0nMcG4+Yz+MY5Q6kNEXT1m1h3oVXF9d/8wWkGJSGzcO30VM1fl0TYtme89t5CfnjeYd1bm8u6qXH563mC6ZaTFOkSpBzUmfjO7F7gcWAEcqsrphG7sEpEmYvf+Eh6fs5mdBcX8ePJAOhxxcba4LMB7q3K5ZFR3rh3fi4sfnMXNzy2iVWoSU07uy/Un9I5N4FLvImnxXwwMdPeSaAcjItHzk1eW8daKHRihKRAfvHoU63L3M2vdLq4Z34v3V+dSVBrg3KFdGdy1DR/8aBJ7i0rp3aGl5sFtZiJJ/BuAZECJX6SJyiss4Z2VO/nmSX1pm5bM795aze6pc/h0cz6lgSAFB8tYvDWf9i1TGN+3PQBd2ragS7gOjzQvkST+ImCRmb1LheTv7t+LWlQiUq9eWbiN8qDzlTE96N0hnS27i1i1s5CLR3Yjv6iM+2esAeDm0/uTlKgxH81dJIn/1fBDRJqItTsLyT9YxrHd2vDPWZv42/vrGd2rHcd0agXAvZcN+3zd/KJS7KUlTB7ahUtGqgRDPIikHv/jDRGIiBydYNCZtX4XOwqK+fm05RwsC5CWnMjBsgCnDerEz88fUul2GekpPHztmAaOVmKpunr8L7j7V8xsKeFpFyty92GVbFZx+0eB84Fcdx8aXtYeeB7oDWwCvuLue+scvYh87tn5W/jJK8sAGNK1DV8ancXy7QVcPb4no3u1j3F00phU1+K/Ofzz/Dru+zHgAUI1/Q+5HXjX3e8xs9vDr2+r4/5F4lZRaTnZm/aya38Jew6UMrxHBn+duZ7hPTK445xBDM/K0EgcqVJ1if/y8By7C929vLY7dvcPzaz3EYsvAk4NP38ceB8lfpGIuTtPzNnMH99Zw96iL86AevclQxnft0MMIpOmpLrEnwX8CRhkZkuA2cAsYI6776nj8Tq7ew6Au+eYWaeqVjSzKcAUgJ49e9bxcCLNy//+ZwWPzd7ESf078o2T+tKzfTotUxN5bNYm9haVcsoAlUqWmpn7F7rvD1/BLAUYA5wATAg/8t298itFh2/bG3itQh9/vrtnVHh/r7u3q2k/Y8aM8ezs7JpWE2nWFm3N5+IHZ/HVCb2464JjSUiwmjeSuGZmC9z9C1fuIxmwmwa0AdqGH58Bc+sYx04z6xoOqCuQW8f9iMSdJ2ZvolVqEj+ePEhJX45KdaN6pgLHAoWEEv1s4P6jHIXzKnAdcE/457Sj2JdIs+fu/P2jDcxZv5tZ63Zz5dgetErVNBpydKr7BPUEUoG1wHZgG5Af6Y7N7FlCF3I7mtk24BeEEv4LZnYDsAX4ct3CFmn8Pss/SF5hCcN7ZNS8chXufGUpz87bSveMNFKTE7hOhdKkHlQ39eJkMzNCrf4TgB8AQ81sD6ELvL+obsfufmUVb51e12BFGruCojIe+mA95x3Xle88s4Ctew4yrk977vvycHq0r90EJu+u3Mmz87Yy5eS+3HHOIEL/HUWOXrV9/B6yDHgDeJPQqJ5+/HeMv4hUcN/bq/nbB+u54IGPyckv5ruT+rEiZx8XPPAxK3P2Rbyf4rIAd/1nOcd0asUPzxqopC/1qsrEb2bfM7PnzGwrodr75wOrgUsB3QYocoRl2wt4Zt4Wzh/WlZP6d+RXFw/lR2cP4j83nUh5wPnnrI0R7+u+t1azdc9BfnXRUFKSVDRN6ld1ffy9gReB7x8aey8ilfvTO2v583traZuWzC8vGkr7limfv9e7Y0vOGtKZt5bv5NcXBytN5PlFpdz16nIWbysgJTGBNbmFXD2uJxP66WYsqX/V9fHf2pCBiDRV+4rLeHDmOk4ZkMk9lx53WNI/5LxhXXl54XZmrd/FpIGdDtt29rrd3D9jNZt2FXHGkE4cKAnQoVUKd5w7uCFPQ+KIxoWJHKV3V+6kNBDku5OOoVObyicuObF/R1q3SOLvH25gbO/25B8s42f/XsaHa/IoDzrt0pN57OvHc0K/jg0cvcQjJX6RGry4YBtz1u+md4d0vjPpGBLDN099vHYXLy/cxtY9RXRp04KR1QzbTE1K5LbJg/j5tGWM/793KS0PkpRgfOOkvkwamMmoXu1I1gQo0kCU+EWqkbuvmDteXkJ6ShIFB8vYsqeI70w6hpc/3caDM9cRDFc8+drE3jXeTXvN+F4c06kV/164nRbJiXxtYm96dWjZAGchcjglfpFqPDFnM+VB59WbJvLigm385b11/GvBNgAuHtGN/zm9P8/N28L1E/tEtL/xfTuoeqbEnBK/SBUKDpbx9NzNnDG4M706tOTWMwdw5pDOLNlWwIgeGQzt3haAn5xXY71CkUZFiV+kEkWl5Xz9sfnsLynnu5OOAcDMGJaVwbCsupdgEGkMdDVJpBK/em0lC7fs5U9XjGTEUdTaEWmM1OIXqaC0PMh7q3J5dt4WvnlSH849rmusQxKpd0r8ImErc/Zx7SPz2LW/hF4d0rn1zIGxDkkkKpT4RYCNuw5w1d8/ITUpkT9fOZKT+3fUZOXSbCnxS1xxd+Zv2suSbfmc2L8jAzu3Jujww38tJujw3JTx9O6osfXSvCnxS1xwd95esZOH3l/Poq3/nU+od4d0WrdIZun2Av54+QglfYkLSvzS7C3ams+vXlvBgs176dk+nV9edCyTBnbiw7V5zFixk+KyALeeOYCLRnSLdagiDUKJX5otd+e3b63moffX07FVKvdcehxfHtPj81o7V4/rxdXjesU4SpGGp8QvzVJeYQkPvb+eR2dt5PIxPfjZBUM0SblIWEz+J5jZJqAQCADl7j4mFnFI87RkWz5f/tscSsqDXDm2J3dfPLTGAmoi8SSWTaBJ7r4rhseXZqg8EOT2l5aSkZ7MUzeMo3/n1rEOSaTR0d++0my4O797ezUrcvbx0NWjlPRFqhCrxO/A22bmwMPuPvXIFcxsCjAFoGfPng0cnjQlu/eX8KMXl7BrfwlLthVw5dieTB7aJdZhiTRasSrSNtHdRwHnAN81s5OPXMHdp7r7GHcfk5mZ2fARSpNQXBbgG09kM2vdLlISE/jOqf24++KhmKlPX6QqMWnxu/tn4Z+5ZvYKMBb4MBaxSMMqDwTJ219Cu/QUWiQfXUmEYNC59YVFLNqaz1+vGsU5KqgmEpEGT/xm1hJIcPfC8POzgF82dBzScIrLAjw2exNPztlMTsFBgg7dM9J44cYJdM9Iq9M+3Z1fvraCN5bu4CfnDlbSF6mFWLT4OwOvhP8UTwKecffpMYgjLpQHgiz7bB9lgSCJCcb63P2cMbgz7Vqm8O+F2/nTu2v56oTQjUwFB8tYn7e/XqcG3F9SzvWPziN7815O6t+RL43qTtv0FP44Yw3XPjKXf393Im1aJNe4n9LyIPuKy0gwY97GPby6eDtvLN3BDSf24RsnRTbtoYiEmLvHOp7xcKEAAAzaSURBVIYajRkzxrOzs2MdRpOzekch1zwyl7zCksOWt26RxLg+HXhv1U7apaew+0ApvTqkU3CwjPyiMr59aj9+fPbAz/vJ9xWXkZqUQGpS7bpmtu4p4n+eXcjS7QX84fIRXDj8vyURPtmwm6v/MZdTB2RyoLSc7hnp/ObSoYcdo6Q8wMIt+fRon84VU+ewdc/Bz99rkZzAjaf04+bT+6s/X6QKZragsvukNJyzmdq6p4jrHp2HAX+5ciStWiRRXBqgU5sW/HPWRlbtKOSsIV24//LhzN2wh3unryKzVSq9OrTkoffXs6OgmON7t+f57K0s3pqPGVw5tie/uGBIjV8AeYUl3D9jNdMWfUaiGQ9eNeoLo2zG9+3ArWcO4HdvraZVahKfbNjD+rz93HhKP04f3IlEM7737ELeWr6TxAQjLTmR2yYPwgzG9GrHsKwMUpI0gZxIXajF3wxlb9rDlCcXUB4I8vy3JjC4a5uIt3V3Hpy5jvveXgPAgM6tOH9YN3IKinl23haG98jggStHkpqUQGbr1MNa2+ty95PVLo1vP7WAWet3c8GwbtxyRn96tE+v9FjBoPPywu2c0K8D8zft4TdvrGTnvhI6tkqlfctk1uzcz3UTepFTUMzXJvZhQr/664ISiQdVtfiV+JuZWet2ccPj8+nWNo1/XDeGvpmt6ryf9JRERvTI+Dy5T1+Www9eWMyB0gAAXdq0ID08WUnH1qnM27iHzm1S2bmvhJ+eN5hvnNS3VscsDwR5f3UeryzczoHScib268g3TuqjrhyROlLijwMzV+fyrScX0LdjS576xjg6tkqt92NsyNvP9OU7aJGUyKdb9hJ0p6QsyIZdBzj72C68umg77Vul8O/vTCQpUV0xIrGkxN9IuHu9tmAf+Xgjy7YXUFhczgdrchnQuTVP3TCOdi1T6u0YtVEWCBII+lGP0ReRo6eLu3VQXBbgjaU5PDN3C9v2HuSSUd25elxPstpV3mddnZLyAHe9uoL3Vu3ktsmD+OesTYzr056fnj8EgNzCYl6Yv5UNeQdISjRSkxI5ZUAmpw7MrLLlnFNwkN+8sZJWqUm0TUvm6nG9+P4ZA2ibXvPwyGhJTkxAOV+kcVPir8Qzc7fwx3fWkF9URmkgSN+OLTm2Wxse/mA9D3+wnlvPHMBNp/Wvdh97D5Ty94828J8ln7G/uJwDpQFKy4Nktk7l1hcWk5xoLN1eQN/MVozt045LHpxNYUk53TPSCLpTWFzOk59splPrVCYe05Gi0nKKSgNcOLwbl43Owsx4Zu4Wgu7856YT6dmh9l9GIhKflPgJjS6Zv2kPndu04LHZm3hs9ibG9m7PyJEZnDIwkwl9O2BmfJZ/kLvfWMl9b68hJSmBS0Zmkdn6i/3o2/YWcfU/5rJlTxGnDexE93ZppCUnMqFfB0b2aMeTn2ziwuHd+em0Zdz5ylLat0whJSmBGd85+fOKkmWBIDNX5fL8/K3M3bCbVi2SKAs4P3pxCdMWfcaonhk89clmTh/USUlfRGpFffzAQ++v597pqwAwg2vH9+Jn5w8huZIultLyIDc8Pp+P1u4iOTE0Rv2k/plkb97DnPW7mbtxD8u2F5CalMA/v3Y8o3u1r/K4xWUB7nlzFS9kb+UfXx3DCcd0rDbOQNB5+MP1PD9/K5t3F9GjfRoPXjWKYVkZR/cLEJFmSRd3q7AyZx8XPvAxJ/fP5OQBmYzu1Y6h3dtWu00g6CzdXsAvXl3O8u0FmEFZwElKMI7Lasuonu244vgeEdeDDwT983lgI1VaHtQNTCJSLV3crcKvX19BmxbJ/O7Lw2kf4UiYxARjRI8MnvjaWO59axWtU5OY0K8Dx/duT8s6zOta26QPKOmLSJ3FdeLP3rSHWet289PzBkec9Ctqm57Mby45LgqRiYhET9w2G92d37+9hg4tU7hqnGb4EpH4EbeJ/8lPNjNnw25uOXMA6Slx/YePiMSZuEz82/NDNz6dMiCTa9TaF5E4E5eJ/8GZ6wgG4TeXHqcCYCISd+Iu8a/L3c+/srdy+fE96jztn4hIUxZXnduPfryRe95cRWpSIt8+tV+swxERiYm4afFv21vEPW+uYny/Drz1/ZPppta+iMSpmCR+M5tsZqvNbJ2Z3d4Qx7z/7TWYwT2XHqcuHhGJaw2e+M0sEXgQOAcYAlxpZkOidbz8olLueHkJLy/cztcm9lFLX0TiXiz6+McC69x9A4CZPQdcBKyo7wP9+d21PDhzHSXlQb59aj9+cOaA+j6EiEiTE4vE3x3YWuH1NmDckSuZ2RRgCkDPnnUba98tI43LRmdxzfhetZpwXESkOYtF4q9s4PwXSoS6+1RgKoSqc9blQJeNzuKy0Vl12VREpNmKxcXdbUCPCq+zgM9iEIeISFyKReKfD/Q3sz5mlgJcAbwagzhEROJSg3f1uHu5md0EvAUkAo+6+/KGjkNEJF7F5M5dd38DeCMWxxYRiXdxc+euiIiEKPGLiMQZJX4RkTijxC8iEmfMvU73RjUoM8sDNtdx847ArnoMpymK999BvJ8/6HcQr+ffy90zj1zYJBL/0TCzbHcfE+s4Yinefwfxfv6g30G8n/+R1NUjIhJnlPhFROJMPCT+qbEOoBGI999BvJ8/6HcQ7+d/mGbfxy8iIoeLhxa/iIhUoMQvIhJnmnXij8Wk7rFmZpvMbKmZLTKz7PCy9mY2w8zWhn+2i3Wc9cnMHjWzXDNbVmFZpedsIX8OfyaWmNmo2EVeP6o4/7vMbHv4c7DIzM6t8N4d4fNfbWZnxybq+mNmPcxsppmtNLPlZnZzeHncfAZqq9km/oae1L2RmeTuIyqMW74deNfd+wPvhl83J48Bk49YVtU5nwP0Dz+mAA81UIzR9BhfPH+AP4Q/ByPCFXEJ/x+4Ajg2vM1fw/9XmrJy4AfuPhgYD3w3fJ7x9BmolWab+Kkwqbu7lwKHJnWPRxcBj4efPw5cHMNY6p27fwjsOWJxVed8EfCEh3wCZJhZ14aJNDqqOP+qXAQ85+4l7r4RWEfo/0qT5e457v5p+HkhsJLQ3N5x8xmoreac+Cub1L17jGJpSA68bWYLwhPWA3R29xwI/ScBOsUsuoZT1TnH0+fipnBXxqMVuvea9fmbWW9gJDAXfQaq1JwTf0STujdDE919FKE/Z79rZifHOqBGJl4+Fw8B/YARQA7w+/DyZnv+ZtYKeAm4xd33VbdqJcuaxe8gUs058cflpO7u/ln4Zy7wCqE/43ce+lM2/DM3dhE2mKrOOS4+F+6+090D7h4E/s5/u3Oa5fmbWTKhpP+0u78cXhzXn4HqNOfEH3eTuptZSzNrfeg5cBawjNB5Xxde7TpgWmwibFBVnfOrwFfDIzvGAwWHugOakyP6rC8h9DmA0PlfYWapZtaH0AXOeQ0dX30yMwMeAVa6+/0V3orrz0C13L3ZPoBzgTXAeuAnsY6nAc63L7A4/Fh+6JyBDoRGNawN/2wf61jr+byfJdSdUUaoNXdDVedM6M/8B8OfiaXAmFjHH6XzfzJ8fksIJbquFdb/Sfj8VwPnxDr+ejj/Ewl11SwBFoUf58bTZ6C2D5VsEBGJM825q0dERCqhxC8iEmeU+EVE4owSv4hInFHiFxGJM0r8EhfMbH/4Z28zu6qe933nEa9n1+f+ReqbEr/Em95ArRJ/BNUrD0v87n5CLWMSaVBK/BJv7gFOCteo/76ZJZrZ78xsfrig2bcAzOzUcI33Zwjd5IOZ/Ttc/G75oQJ4ZnYPkBbe39PhZYf+urDwvpeF50i4vMK+3zezF81slZk9Hb77FDO7x8xWhGO5r8F/OxIXkmIdgEgDux34obufDxBO4AXufryZpQKzzOzt8LpjgaEeKl8M8HV332NmacB8M3vJ3W83s5vcfUQlx7qUUJG04UDH8DYfht8bSagm/mfALGCima0gVF5hkLu7mWXU+9mLoBa/yFmE6rYsIlTKtwOh+jUA8yokfYDvmdli4BNCRb76U70TgWc9VCxtJ/ABcHyFfW/zUBG1RYS6oPYBxcA/zOxSoOioz06kEkr8Eu8M+B//70xVfdz9UIv/wOcrmZ0KnAFMcPfhwEKgRQT7rkpJhecBIMndywn9lfESoUlDptfqTEQipMQv8aYQaF3h9VvAt8NlfTGzAeHKpkdqC+x19yIzG0Roir9Dyg5tf4QPgcvD1xEygZOpphJmuJ58Ww9Nk3gLoW4ikXqnPn6JN0uA8nCXzWPAnwh1s3wavsCaR+VTU04HbjSzJYSqWn5S4b2pwBIz+9Tdr66w/BVgAqFqqQ782N13hL84KtMamGZmLQj9tfD9up2iSPVUnVNEJM6oq0dEJM4o8YuIxBklfhGROKPELyISZ5T4RUTijBK/iEicUeIXEYkz/w/hxdEUc/flTAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "### Training Pong ###\n", "# You can run this cell and stop it anytime in the middle of training to save \n", "# a progress video (see next codeblock). To continue training, simply run this \n", "# cell again, your model will pick up right where it left off. To reset training,\n", "# you need to run the cell above. \n", "\n", "games_to_win_episode = 21 # this is set by OpenAI gym and cannot be changed.\n", "\n", "# Main training loop\n", "while iteration < MAX_ITERS:\n", "\n", " plotter.plot(smoothed_reward.get())\n", "\n", " tic = time.time()\n", " # RL agent algorithm. By default, uses serial batch processing.\n", " # memories = collect_rollout(batch_size, env, pong_model, choose_action)\n", "\n", " # Parallelized version. Uncomment line below (and comment out line above) to parallelize\n", " memories = parallelized_collect_rollout(batch_size, envs, pong_model, choose_action)\n", " print(time.time()-tic)\n", "\n", " # Aggregate memories from multiple batches\n", " batch_memory = aggregate_memories(memories)\n", "\n", " # Track performance based on win percentage (calculated from rewards)\n", " total_wins = sum(np.array(batch_memory.rewards) == 1)\n", " total_games = sum(np.abs(np.array(batch_memory.rewards)))\n", " win_rate = total_wins / total_games\n", " smoothed_reward.append(100 * win_rate)\n", " \n", " # Training!\n", " train_step(\n", " pong_model,\n", " optimizer,\n", " observations = np.stack(batch_memory.observations, 0),\n", " actions = np.array(batch_memory.actions),\n", " discounted_rewards = discount_rewards(batch_memory.rewards)\n", " )\n", "\n", " # Save a video of progress -- this can be played back later\n", " if iteration % 50 == 0:\n", " show_video_of_model(pong_model, \"Pong-v0\")\n", " \n", " iteration += 1 # Mark next episode" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "> Note: Actually, Kernel is crashed due to hardware resource limit. In order to training it more than 200 epochs, you need to use colab environment." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Finally we can put our trained agent to the test! It will play in a newly instantiated Pong environment against the \"computer\", a base AI system for Pong. Your agent plays as the green paddle. Let's watch the match instant replay!" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "show_video(\"Pong-v0\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 3.10 Conclusion and submission information\n", "\n", "That's it! Congratulations on training two RL agents and putting them to the test! We encourage you to consider the following:\n", "\n", "* How does the agent perform?\n", "* Could you train it for shorter amounts of time and still perform well?\n", "* What are some limitations of the current representation i.e., difference of current and previous frames? How is this reflected in the agent's behavior? What could be done to generate an improved representation?\n", "* How does the complexity of Pong relative to CartPole alter the rate at which the agent learns and its performance? \n", "* What are some things you could change about the agent or the learning process to potentially improve performance?\n", "\n", "Try to optimize your **Pong** model and algorithm to achieve improved performance." ] } ], "metadata": { "accelerator": "GPU", "colab": { "collapsed_sections": [], "name": "PyTorch_Tutorial.ipynb", "provenance": [] }, "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.7.10" } }, "nbformat": 4, "nbformat_minor": 4 }