examples/pennylane/0_Getting_started/0_Getting_started.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Combining PennyLane with Amazon Braket"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<div class=\"alert alert-block alert-info\">\n",
    "<b>What is PennyLane?</b> <a href=\"https://pennylane.ai/install.html\">PennyLane</a> is a Python library for differentiable programming of quantum computers, allowing you to train a quantum computer the same way as a neural network.\n",
    "</div>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<div align=\"center\"><img src=\"pl_braket.png\"/></div>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "PennyLane integrates with Amazon Braket to add additional features for quantum machine learning and optimization. This introductory tutorial walks you through how to train a quantum circuit using Amazon Braket simulators and PennyLane's automatic differentiation capabilities."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Setup"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "PennyLane is already installed on Braket notebook instances. On a local machine, PennyLane can be installed by following [these](https://pennylane.ai/install.html) instructions. It can then be imported with:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-02-19T21:42:09.404557Z",
     "start_time": "2024-02-19T21:42:03.275541Z"
    }
   },
   "outputs": [],
   "source": [
    "import pennylane as qml\n",
    "from pennylane import numpy as np"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To use Braket as a backend in PennyLane, we have to create a PennyLane device. Here we will first create a device that uses the local Braket simulator that runs on your local laptop (or on the server that hosts this notebook)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-02-19T21:42:10.351825Z",
     "start_time": "2024-02-19T21:42:09.405882Z"
    }
   },
   "outputs": [],
   "source": [
    "wires = 2  # Number of qubits\n",
    "\n",
    "dev = qml.device(\"braket.local.qubit\", wires=wires)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Below we will also show you how to scale out simulations to the AWS cloud."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Defining a circuit"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We will choose a simple two-qubit circuit with two controllable rotations and a CNOT gate."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-02-19T21:42:10.355051Z",
     "start_time": "2024-02-19T21:42:10.352639Z"
    }
   },
   "outputs": [],
   "source": [
    "@qml.qnode(dev)\n",
    "def circuit(params):\n",
    "    qml.RX(params[0], wires=0)\n",
    "    qml.RY(params[1], wires=1)\n",
    "    qml.CNOT(wires=[0, 1])\n",
    "    return qml.expval(qml.PauliZ(1))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The ``qml.qnode(dev)`` decorator binds the circuit to the local Braket device. Now, every time that ``circuit()`` is called, the quantum computation defined in the function above will be executed with Braket."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<div class=\"alert alert-block alert-info\">\n",
    "<b>Note</b> PennyLane also supports automatic differentiation with PyTorch and TensorFlow <a href=\"https://pennylane.readthedocs.io/en/stable/introduction/interfaces.html\">interfaces</a>. The choice of interface can be specified using:\n",
    "<code>@qml.qnode(dev, interface=\"&lt;interface&gt;\")</code>.\n",
    "</div>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Evaluating the circuit and accessing its gradient"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's set some values for our controllable parameters:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-02-19T21:42:10.359085Z",
     "start_time": "2024-02-19T21:42:10.356833Z"
    }
   },
   "outputs": [],
   "source": [
    "params = np.array([0.1, 0.2], requires_grad=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The circuit can be evaluated with these parameters using"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-02-19T21:42:10.393070Z",
     "start_time": "2024-02-19T21:42:10.359915Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Expectation value of circuit: 0.9751703272018161\n"
     ]
    }
   ],
   "source": [
    "print(\"Expectation value of circuit:\", circuit(params))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-02-19T21:42:10.397595Z",
     "start_time": "2024-02-19T21:42:10.394259Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Drawing of circuit:\n",
      "\n",
      "0: ──RX(0.10)─╭●─┤     \n",
      "1: ──RY(0.20)─╰X─┤  <Z>\n"
     ]
    }
   ],
   "source": [
    "print(\"Drawing of circuit:\\n\")\n",
    "print(qml.draw(circuit)(params))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "A crucial element of machine learning and optimization is accessing the gradient of a model with respect to its parameters. This functionality is built into PennyLane:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-02-19T21:42:10.400540Z",
     "start_time": "2024-02-19T21:42:10.398616Z"
    }
   },
   "outputs": [],
   "source": [
    "dcircuit = qml.grad(circuit)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Here, ``dcircuit`` is a callable function that evaluates the gradient of the circuit, i.e., its partial derivatives with respect to the controllable parameters."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-02-19T21:42:10.523442Z",
     "start_time": "2024-02-19T21:42:10.401286Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "array([-0.0978434 , -0.19767681])"
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dcircuit(params)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Training the circuit"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Suppose we now want to minimize the output of the circuit by updating its parameters. This can be done using gradient-based optimization.\n",
    "\n",
    "First, an optimizer is fixed:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-02-19T21:42:10.526039Z",
     "start_time": "2024-02-19T21:42:10.524244Z"
    }
   },
   "outputs": [],
   "source": [
    "opt = qml.GradientDescentOptimizer(stepsize=0.1)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The next step is to run the optimizer for a chosen number of iterations:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-02-19T21:42:11.062067Z",
     "start_time": "2024-02-19T21:42:10.527962Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Minimized circuit output: 0.37261070647126565\n",
      "Optimized parameters: [0.4839502  1.13630274]\n"
     ]
    },
    {
     "data": {
      "text/plain": "<Figure size 640x480 with 1 Axes>",
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjcAAAGzCAYAAADT4Tb9AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy88F64QAAAACXBIWXMAAA9hAAAPYQGoP6dpAABOkElEQVR4nO3deVxU5f4H8M+ZYRlQGEAFBkXADUFMA0XBzHLBpSyv/VIrt7ablpnSptfKtHuz1SxTqpumpbfMPcsszF1JFEFFcAlQEAcQ0AFEtpnz+4OYGhmGAWEOM3zer9e8as48Z+Z7Oo3z8TzPeR5BFEURRERERDZCJnUBRERERE2J4YaIiIhsCsMNERER2RSGGyIiIrIpDDdERERkUxhuiIiIyKYw3BAREZFNYbghIiIim8JwQ0RERDaF4YaIiIhsip2UH37gwAG8//77SEhIgFqtxtatWzFu3DiT++zfvx/R0dE4c+YMfHx88Morr2DGjBlmf6ZOp8OVK1fg4uICQRBu8wiIiIjIEkRRRHFxMXx8fCCTmb42I2m4uXHjBvr06YPHH38cDz30UL3tMzIyMGbMGDz99NNYt24dDh8+jGeffRYdOnQwa38AuHLlCnx9fW+3dCIiIpJAVlYWOnXqZLKN0FIWzhQEod4rN6+++ip++OEHpKam6rfNmDEDJ0+eRFxcnFmfo9Fo4ObmhqysLLi6ut5u2URERGQBRUVF8PX1xfXr16FUKk22lfTKTUPFxcUhKirKYNvIkSOxatUqVFZWwt7evtY+5eXlKC8v1z8vLi4GALi6ujLcEBERWRlzhpRY1YDinJwceHl5GWzz8vJCVVUV8vPzje6zZMkSKJVK/YNdUkRERLbNqsINUDux1fSq1ZXk5s+fD41Go39kZWU1e41EREQkHavqlvL29kZOTo7Btry8PNjZ2aFdu3ZG93F0dISjo6MlyiMiIqIWwKqu3ERERCA2NtZg26+//op+/foZHW9DRERErY+k4aakpARJSUlISkoCUH2rd1JSEjIzMwFUdylNnTpV337GjBm4dOkSoqOjkZqaitWrV2PVqlV46aWXpCifiIiIWiBJu6WOHz+Oe++9V/88OjoaADBt2jSsWbMGarVaH3QAICAgADt37sTcuXOxYsUK+Pj44JNPPjF7jhsiIiKyfS1mnhtLKSoqglKphEaj4a3gREREVqIhv99WNeaGiIiIqD4MN0RERGRTrOpW8JZMqxMRn1GIvOIyeLooEB7gAbmMC3MSERFZGsNNE9iVrMaiHSlQa8r021RKBRaODcaoEJWElREREbU+7Ja6TbuS1Zi57oRBsAGAHE0ZZq47gV3JaokqIyIiap0Ybm6DVidi0Y4UGLvdrGbboh0p0Ooaf0OaViciLq0A25OyEZdWcFvvRURE1BqwW+o2xGcU1rpi83ciALWmDIt3nMFd3Tugk7sTfD2c0dbRvP/s7O4iIiJqOIab25BXXHew+bu1cZewNu6S/rmbs3110HF3Rid3J3Ryd4avR/U/O7o5oY2jnb6769brNDXdXTGTQxlwiIiIjGC4uQ2eLgqz2oX7u+NmpQ6Xr5XiWmklrv/5SM4uMtre3dkeJeVVdXZ3Caju7hoR7H1bd2TxDi8iIrJFDDe3ITzAAyqlAjmaMqNBRADgrVTg239G6ENDSXkVLl8rxeXCm7h8rRRZ16r/efnaTVy+dhOam5W4Vlpp8nNrurvmbEjC4G7t0aVDG3Tp0BbuzvYQBPPCCbu8iIjIVnH5hdtU030EwCDg1ESMhnYfaW5WYv3vl/DeL+caXIvSyb466LRv++c/q0OPXztnKOzltWq+9cQ3tmYiIqLm1pDfb165uU2jQlSImRxa6yqIdyOvgiid7HFnZ3ez2o4M9kJppRbpV28g+3r1VZ/EzOtIzLxu0E4QgI5uTujSoS382zljW2J2s3d5ERERSYXhpgmMClFhRLB3k41fMbe7a+XkMP1n3KzQ4mLBDaRfvYH0qyXIyL+BtPzqfy8uq9J3ex2o57NrurziMwoR0bVdo+oHOJ6HiIikw26pFqqpurtEUUR+SQUy/gw6v6bkYs/ZvHr3c3Wyw52+7uipckGwyhU9vV3RpUMb2MvrnxqJ43mIiKipNeT3m+GmBWuOkBCXVoBH/vt7o/Z1kMvQzbOtQeDpqXJB+7aOBjVzPA8RETU1hhsTrCncAE3fvaPVibjr3T0mu7y8XBX45JE7cSGvGKnqIpxVF+NsTjFKyquMvmcHF0f09HZBT28XfH/8MjQ3jd/tVdOddujVoeyiIiKiBmG4McHawk1zaEyXlyiKuHztJlLVRUhVF+NsThHO5hTjYsENNPT/oG+fHnhb43mIiKj1YbgxgeGmWlN1ed0or8L53OorOztPq3HwQn69+wwP8sTD/Xxxp68bPF3NmwixBgcqExG1Tgw3JjDc/KWpg0JjxvN0dHNCX1833NnZDX193RDSUWkwJ8/fcaAyEVHrxXBjAsNN86lvPA9QPY/PyF5eOHVZg3O5xbW6tOxkAoJUrgaBJ6B9G/xyJocDlYmIWjGGGxMYbppXQ8bzlJRX4dTl60jKqp54MCnrOq4Wl9d6T1eFHcqqdKio0hn9TA5UJiKyfQw3JjDcNL/Gdh+Joojs6zeRlHUdSZnXkZh1HcnZGpTXEWpuxYHKRES2i+HGBIYby2iq8TwVVTp8fiANH/56vt629wZ2wJQIP/T394CLwr4xZRMRUQvFtaVIcnKZ0CRXURzsZOjn52FW273nrmLvuauQywTc0UmJyK7tENm1PcL83OscpPx3vBOLiMg2MNxQi1ffWlsA4OZsj5G9vHE0vQAXC0r1C4iu2JsGB7kMoX5uiOzaHoO6tcMdndxqLSPBO7GIiGwHu6XIKjRkoPLla6WISytAXFoBDqflI7fIcJCys4Mc4QEe+is7mQWleO5/vBOLiKgl45gbExhurFdjrq6IooiM/Bs48mfYiUsvQOGNCoM2AlDnFSHeiUVE1DIw3JjAcGPdbndcjE4n4lxu8Z9hJx+H/8jHzcr678binVhERNLigGKyWbc7UFn25ySBQSpXPHlXALaeuIy535+sd79dyWoE+7hC6cS7sIiIWjqGG2rVvJVOZrVbG3cJ649mYmCXdhge5InhwV7o5O7czNUREVFjsFuKWjVzloxo6yiHl6sCaVdvGGwPUrliRLAXRgR5IaSjKwTBePcYbzEnIrp9HHNjAsMN3crcO7Ey8m9gd0ouYlNzcfxiIXR/a+ztqsDwYE+MCPbGwC4ecLST69+bt5gTEd0+hhsTGG7ImIaGkMIbFdh7Ng+xKbk4cOEqSiu0+tfaOtphSI8O8HR1wJrDl3iLORFRE2C4MYHhhurS2O6jskot4tIKEJuai90pucgzsvjnrXiLORFRwzDcmMBwQ81JpxNxKluDNYczsC3pSr3teYs5EZF5GvL7LTP5KhE1iEwmoK+vG+7t6WlW+90puSiv0tbfkIiIzMZbwYmagaeLwqx2qw5nYNOJy7j/DhXGh3ZEaGf3Ou+6IiIi8zDcEDUDcxb7bOMoRxsHOfKKK7D+aCbWH82EXztnjOvbEf+4syP827exaM1ERLaCY26Imok5t5iPCPZGXFoBtiRexq7kHIO7rkI7u+EfoZ1wf28V3Ns4GP0MzqFDRK0FBxSbwHBDltSQW8xLK6rw65lcbEnMxqELV/Xz6NjLBdwb6InxoR1xb09PzqFDRK0Sw40JDDdkaY25upJXVIYfTl7BlhPZSFEX6bcrnexx3x0q+CgV+PDX85xDh4haDYYbExhuyNqcyynGlsTL2J54BTlFZfW25xw6RGSLeCs4kQ0J9HbB/NFBODxvKNY/NQB3d29vsr0IQK0pQ3xGoWUKJCJqYXi3FJGVkMsEDOrWHvkl5ThwIb/e9nnF9V/lISKyRbxyQ2RlzJ1DZ9uJbKRcKaq/IRGRjZE83KxcuRIBAQFQKBQICwvDwYMHTbZfsWIFgoKC4OTkhMDAQHz99dcWqpSoZaiZQ6e+0TR7z1/FmE8OYsLncfj5tBpVWp1F6iMikpqk4WbDhg2YM2cOFixYgMTERAwePBijR49GZmam0fYxMTGYP38+3nzzTZw5cwaLFi3Cc889hx07dli4ciLpyGUCFo4NBoBaAUf48/HyyB64/w4V5DIB8RmFmLn+BO5+by9i9qXh2o0KS5dMRGRRkt4tNWDAAISGhiImJka/LSgoCOPGjcOSJUtqtY+MjMSgQYPw/vvv67fNmTMHx48fx6FDh8z6TN4tRbbCnHlucjRlWH/0Ev53NBMFf4YaRzsZxvXtiGmR/gj24XeAiKxDQ36/JRtQXFFRgYSEBMybN89ge1RUFI4cOWJ0n/LycigUhuMNnJycEB8fj8rKStjb2xvdp7y8XP+8qIhjEMg2jApRYUSwt8k5dLyVCrwYFYjn7u2GH0+pseZIBpKzi7DheBY2HM9CeIAHHo/0x4hgL9jJ/7qQy5mPiciaSRZu8vPzodVq4eXlZbDdy8sLOTk5RvcZOXIkvvzyS4wbNw6hoaFISEjA6tWrUVlZifz8fKhUtSctW7JkCRYtWtQsx0AkNblMQETXdvW2U9jL8X9hnfBQaEecyLyGrw5fxK7kHMRnFCI+oxA+SgWmRPhjUn9fHM0o4MzHRGTVJL8V/NYVkEVRrHNV5Ndffx05OTkYOHAgRFGEl5cXpk+fjvfeew9yudzoPvPnz0d0dLT+eVFREXx9fZvuAIisiCAICPPzQJifh0GX1RVNGd7ddRYf/noOVbraPdU5mjLMXHeCMx8TkVWQbEBx+/btIZfLa12lycvLq3U1p4aTkxNWr16N0tJSXLx4EZmZmfD394eLiwvatzc+sZmjoyNcXV0NHkT0V5fV4XlD8cHDfdDLx8VosAH+Wvhz0Y4UaOtoQ0TUUkgWbhwcHBAWFobY2FiD7bGxsYiMjDS5r729PTp16gS5XI7vvvsO999/P2Qyye9qJ7JKNV1Wr90XbLIdZz4mImshabdUdHQ0pkyZgn79+iEiIgJffPEFMjMzMWPGDADVXUrZ2dn6uWzOnz+P+Ph4DBgwANeuXcPSpUuRnJyMtWvXSnkYRDYhr7i8/kYAcs1Y34qISEqShpuJEyeioKAAixcvhlqtRkhICHbu3Ak/Pz8AgFqtNpjzRqvV4sMPP8S5c+dgb2+Pe++9F0eOHIG/v79ER0BkO8yd+fjDX8/B0U6Gkb28IeMdVETUAnFVcCICUH37913v7kGOpgx1/aEg4K/xNz29XTBneHdEBTPkEFHz46rgRNRg5sx8/MHDfTB7WHe4ONrhbE4xZqw7gfuWH8IvZ3LQyv6eREQtGK/cEJEBc2Y+vl5agVWHMvDV4YsoKa8CAPTyccWc4T0wPMizzukciIgaqyG/3ww3RFSLuTMUXy+twJcHM/DV4QzcqNACAEI6umLOsB4YxpBDRE2I4cYEhhuipnftRgX+ezAda49c1IecOzopMWd4d9wbaBhyuLQDETUGw40JDDdEzafwbyGn9M+Q06eTEnOG98A9gR3wy5kcLu1ARI3CcGMCww1R8ysoKccXB9Px9ZFLuFlZHXL82znjYkFprbY112y4tAMRmcK7pYhIUu3aOmL+6CAcevVePHN3FyjsZEaDDcClHYio6THcEFGzadfWEfPHBGHZpL4m23FpByJqSgw3RNTsyqt0ZrXLK+bSDkR0+xhuiKjZmbu0Q/u2js1cCRG1Bgw3RNTswgM8oFIqas18fKu3fkzB0fQCi9RERLaL4YaIml19SzsAgLODHGdzijHxi98x+9tEqDU3LVojEdkOhhsisohRISrETA6Ft9Kwi8pbqcBnk0Nx6NWheGxAZwgC8MPJKxj6wX6s2PsHyv68lZyIyFyc54aILKq+GYqTszV484czOH7pGgDAr50zXr8vmMs5ELVynMTPBIYbopZPFEVsT7qCt3emIq+4HAAwpEcHvDE2GF07tJW4OiKSAsONCQw3RNajpLwKK/b+gS8PpqNSK8JeLuCJQQGYNbQbXBT2UpdHRBbEcGMCww2R9cnIv4G3fkzBnrN5AIAOLo6YP7onxvXtCNmfXVpckJPItjHcmMBwQ2S99pzNxeIdKfqlHEI7u2HRAyHIvl7KBTmJbBzDjQkMN0TWrbxKi9WHLmL5ngv6lceN4YKcRLaFC2cSkc1ytJNj5j1dsfele/Bgn7pDCxfkJGq9GG6IyCp5uSowKdzPZBsuyEnUOjHcEJHVMnehTS7ISdS6MNwQkdUyd0FOc9sRkW1guCEiq2XugpzHLhagSquzSE1EJD2GGyKyWqYW5Py7pbEX8FDMEfyRV2KZwohIUgw3RGTV6lqQU6VUIOaxUCyd0AcuCjucvKzBfZ8cxKpDGdDx7ikim8Z5bojIJpiaoVituYlXNp3CwQv5AIABAR744OE+8PVwlrJkImoATuJnAsMNUeskiiLWH83E2ztTUVqhRRsHOV67PxiT+vtytXEiK8BJ/IiIbiEIAiYP9MPPLwxGf3933KjQYv6W03h8zTHkFvFWcSJbwnBDRK2KX7s2+O6fEVgwJggOdjLsO3cVUR8dwPakbLSyC9lENovhhohaHblMwNN3d8FPz9+F3h2V0NysxAvfJeHZ9SdQUFIudXlEdJsYboio1eru5YItz0Zi7vAesJMJ+Dk5ByOXHcCvZ3KkLo2IbgPDDRG1avZyGV4Y3h3bnhuEHl5tkV9SgX9+k4Do75OguVkJoPpOrLi0AmxPykZcWgEX4iRq4Xi3FBHRn8qrtFgaex5fHEiHKFbPlfNwv07YePwy1Jq/Bh2rlAosHBuMUSF1r0pORE2Lt4KbwHBDRPVJuFSIF78/iYsFpUZfr7lxPGZyKAMOkYXwVnAiotsQ5ueBHc/fBWcHudHXa/5GuGhHCruoiFoghhsiIiOSs4tQWqGt83URgFpThviMQssVRURmYbghIjIir9i8if3MbUdElsNwQ0RkhKeLov5GDWhHRJbDcENEZER4gAdUSgVMrTolALh2o8JSJRGRmRhuiIiMkMsELBwbDAB1BhwRwLP/O4F3d51FlVZnsdqIyDSGGyKiOowKUSFmcii8lYZdTyqlAp8+eieeuisAABCzLw3Tvorn0g1ELQTnuSEiqodWJyI+oxB5xWXwdFEgPMADcln19ZwdJ6/g1c2nUFqhhY9SgZWTw9DX103agolsECfxM4Hhhoia2vncYjzzTQIy8m/AQS7Dogd74ZHwzlKXRWRTOIkfEZEF9fBywfZZgxAV7IUKrQ7zt5zGq5tOoayy7nlyiKj5SB5uVq5ciYCAACgUCoSFheHgwYMm269fvx59+vSBs7MzVCoVHn/8cRQUFFioWiIi41wV9vhschheGRUImQBsOJ6Fhz+LQ1ah8SUciKj5SBpuNmzYgDlz5mDBggVITEzE4MGDMXr0aGRmZhptf+jQIUydOhVPPvkkzpw5g40bN+LYsWN46qmnLFw5EVFtMpmAZ+/phq+fGAB3Z3ucztZg7KeHcOD8ValLI2pVJB1zM2DAAISGhiImJka/LSgoCOPGjcOSJUtqtf/ggw8QExODtLQ0/bbly5fjvffeQ1ZWllmfyTE3RGQJ2ddvYua6BJy6rIEgAC+O6IFn7+kGmczUzDlEVBerGHNTUVGBhIQEREVFGWyPiorCkSNHjO4TGRmJy5cvY+fOnRBFEbm5udi0aRPuu+++Oj+nvLwcRUVFBg8ioubW0c0J3z8TgUn9fSGKwAe/nscz6xJQVFYpdWlENk+ycJOfnw+tVgsvLy+D7V5eXsjJyTG6T2RkJNavX4+JEyfCwcEB3t7ecHNzw/Lly+v8nCVLlkCpVOofvr6+TXocRER1UdjL8c5Dd+Cd8b3hYCdDbEouHlh+COdyigFU32Iel1aA7UnZiEsr4ArjRE3ETuoCBMHwEq0oirW21UhJScHs2bPxxhtvYOTIkVCr1Xj55ZcxY8YMrFq1yug+8+fPR3R0tP55UVERAw4RWdSk8M4I9nHFzHUncLGgFONWHMYj4b74OTkHas1fC2+qlAosHBuMUSEqCaslsn6SjbmpqKiAs7MzNm7ciH/84x/67S+88AKSkpKwf//+WvtMmTIFZWVl2Lhxo37boUOHMHjwYFy5cgUqVf1/IHDMDRFJpfBGBWZ/m4hDf+Qbfb3mr3Uxk0MZcIhuYRVjbhwcHBAWFobY2FiD7bGxsYiMjDS6T2lpKWQyw5LlcjmA6is+REQtmUcbB6ye3h9tHOVGX6/5U2zRjhR2URHdBklvBY+OjsaXX36J1atXIzU1FXPnzkVmZiZmzJgBoLpLaerUqfr2Y8eOxZYtWxATE4P09HQcPnwYs2fPRnh4OHx8fKQ6DCIisyVcuoYb5XVP7icCUGvKEJ9RaLmiiGyMpGNuJk6ciIKCAixevBhqtRohISHYuXMn/Pz8AABqtdpgzpvp06ejuLgYn376KV588UW4ublh6NChePfdd6U6BCKiBskrLqu/UQPaEVFtXFuKiMiC4tIK8Mh/f6+33bdPD0RE13YWqIjIOljFmBsiotYoPMADKqUCpqbya+toh/7+7haricjWMNwQEVmQXCZg4dhgAKgz4JSUV+H17WdQpdVZrjAiG8JwQ0RkYaNCVIiZHApvpcJgu0qpwKT+vhAE4Nv4TMxYdwI3K7iyOFFDccwNEZFEtDoR8RmFyCsug6eLAuEBHpDLBOxKVmP2d0moqNIhtLMbVk3rD/c2DlKXSySphvx+M9wQEbVA8RmFeGrtMRSVVaFLhzZY+3g4fD2cpS6LSDIcUExEZOXCAzywaWYkVEoF0q/ewEMxR5ByhQv/EpmD4YaIqIXq4eWCLc9GItDLBXnF5Zj4eRyO1LF0AxH9heGGiKgFUymd8P2MCIQHeKC4vArTvorHDyevSF0WUYvGcENE1MIpnezx9RPhGNPbG5VaEbO/TcSXB9OlLouoxWK4ISKyAgp7OZY/Eorpkf4AgH//lIp//5gCHRfYJKqF4YaIyErUTAA4b3RPAMCXhzIwZ0MSyqs4Fw7R3zHcEBFZEUEQMGNIVyyd0Ad2MgE/nLyCx786huKySqlLI2oxGG6IiKzQ+NBOWD29P9o4yHEkrQATPv8deUVcSZwIYLghIrJad/fogA3PRKB9Wwekqovwj5VHkHa1BED17MdxaQXYnpSNuLQCaDk2h1oRzlBMRGTlMgtKMXX1UVwsKIWbsz2eHtwF636/BLXmrys5KqUCC8cGY1SISsJKiRqPyy+YwHBDRLaooKQcT6w5hpOXNUZfr1mBPGZyKAMOWSUuv0BE1Mq0a+uIdU8NgKOd8T/Wa/4Wu2hHCruoyOYx3BAR2Yjk7CKUV+nqfF0EoNaUIT6j0HJFEUmA4YaIyEbkFZt3t5S57YisFcMNEZGN8HRRNGk7ImvFcENEZCPCAzygUir0g4eNUSkVCA/wsFhNRFJguCEishE1yzMAqDPg9PNzh1xmKv4QWT+GGyIiGzIqRIWYyaHwVhp2Pbko7AAAO06psWLvH1KURmQxdlIXQERETWtUiAojgr0Rn1GIvOIyeLpUd0V9fiAN7+06h/d/OQdBAJ69p5vUpRI1C4YbIiIbJJcJiOjazmDbs/d0gygC7/9yDu/tOge5IOCZIV0lqpCo+bBbioioFXnu3m6IHtEDALDk57P48mC6xBURNT2GGyKiVmb2sO54YVh3AMC/f0rFqkMZEldE1LQYboiIWqE5w7tj9tDqMTdv/ZiCrw4z4JDtYLghImqFBEHA3BE98Ny91WNuFu1IwddxF6UtiqiJMNwQEbVSgiDgpahAzLynOuC8sf0Mvvn9ksRVEd0+hhsiolZMEAS8MjIQz9zdBQDw+rZk/O9opsRVEd0ehhsiolZOEATMG90TT90VAAD419bT+C6eAYesF8MNERFBEAQsuC8ITwyqDjjzt57G98eyJK6KqHEYboiICEB1wHn9/iBMj/SHKAKvbjmFTQmXpS6LqMEYboiISE8QqhffnBrhB1EEXt50EltOMOCQdWG4ISIiA4IgYNEDvTB5YGeIIvDSxpPYlpgtdVlEZuPaUkREVIsgCFj8QAi0OuDb+ExEf58EQQAe7NsRWp1Ya1FOuUyQumQiPYYbIiIySiYT8J9xIRBFEd8dy8LcDUk4na3BT6fUUGvK9O1USgUWjg3GqBCVhNUS/YXdUkREVCeZTMDb/+iNCf06QScCXx7MMAg2AJCjKcPMdSewK1ktUZVEhhhuiIjIpOorOL3hZC83+rr45z8X7UiBVicabUNkSQw3RERUr+OXruFmpbbO10UAak0Z4jMKLVcUUR0YboiIqF55xWX1N2pAO6LmxHBDRET18nRRNGk7oubEcENERPUKD/CASqmAqRu+Vcrq28KJpCZ5uFm5ciUCAgKgUCgQFhaGgwcP1tl2+vTpEASh1qNXr14WrJiIqPWRy6pnLgZQZ8B5dVRPzndDLYKk4WbDhg2YM2cOFixYgMTERAwePBijR49GZqbx1Wg//vhjqNVq/SMrKwseHh54+OGHLVw5EVHrMypEhZjJofBWGnY9CX/mme1J2ajS6iSojMiQIIqiZPftDRgwAKGhoYiJidFvCwoKwrhx47BkyZJ699+2bRvGjx+PjIwM+Pn5mfWZRUVFUCqV0Gg0cHV1bXTtRESt1a0zFNvLBUxedRRllTpM6u+LJeN7QxB4BYeaVkN+vyWbobiiogIJCQmYN2+ewfaoqCgcOXLErPdYtWoVhg8fbjLYlJeXo7y8XP+8qKiocQUTERGA6i6qiK7tDLYtfyQUz3xzHN8dy4JK6YQXhneXqDoiCbul8vPzodVq4eXlZbDdy8sLOTk59e6vVqvx888/46mnnjLZbsmSJVAqlfqHr6/vbdVNRES1jQj2wuIHQwAAH+0+j++PZ0lcEbVmkg8ovvXSpSiKZl3OXLNmDdzc3DBu3DiT7ebPnw+NRqN/ZGXxC0dE1BwmD/TDc/d2BQDM33Ia+87lSVwRtVaNCjeLFy9GaWlpre03b97E4sWLzXqP9u3bQy6X17pKk5eXV+tqzq1EUcTq1asxZcoUODg4mGzr6OgIV1dXgwcRETWPl6ICMf7O6pXDn11/AsnZGqlLolaoUeFm0aJFKCkpqbW9tLQUixYtMus9HBwcEBYWhtjYWIPtsbGxiIyMNLnv/v378ccff+DJJ580v2giImp2giDgnYfuwF3d2qO0QovpXx1DVmHtvwwTNadGhZu6uo5OnjwJDw/zJ3CKjo7Gl19+idWrVyM1NRVz585FZmYmZsyYAaC6S2nq1Km19lu1ahUGDBiAkJCQxpRPRETNyMFOhpjJoQhSuSK/pBzTvorHtRsVUpdFrUiD7pZyd3fXT5zXo0cPg4Cj1WpRUlKiDybmmDhxIgoKCrB48WKo1WqEhIRg586d+ruf1Gp1rTlvNBoNNm/ejI8//rghpRMRkQW5KOyx5vH+GL/yCNKv3sBTXx/H+qcGQFHHyuJETalB89ysXbsWoijiiSeewLJly6BUKvWvOTg4wN/fHxEREc1SaFPhPDdERJZzIbcYD8UcQVFZFUb28sLKx8I4izE1SkN+vxs1id/+/fsxaNAg2NlJNk1OozHcEBFZ1tH0AkxZFY8KrQ7TIvzw5gO9OMkfNVhDfr8bNebGxcUFqamp+ufbt2/HuHHj8K9//QsVFexXJSKivwzo0g4fTewLQQDWxl3CFwfSpS6JbFyjws0zzzyD8+fPAwDS09MxceJEODs7Y+PGjXjllVeatEAiIrJ+992hwmv3VS+8ueTns9ielC1xRWTLGhVuzp8/j759+wIANm7ciCFDhuB///sf1qxZg82bNzdlfUREZCOevCsAT94VAAB4aeNJHEnLl7gislWNvhVcp6te+XX37t0YM2YMAMDX1xf5+fyflYiIjFswJgj33aFCpVbEM18n4GwO1/ujpteocNOvXz/8+9//xjfffIP9+/fjvvvuAwBkZGTUO7swERG1XjKZgA8f7oPwAA8Ul1dh+upjUGtuSl0W2ZhGhZtly5bhxIkTmDVrFhYsWIBu3boBADZt2lTv7MJERNS6Kezl+O+Ufuju2RY5RWWYvvoYNDcrpS6LbEijbgWvS1lZGeRyOezt7ZvqLZscbwUnImoZsq/fxD9WHEZecTkGdvHA6un9cTJLg7ziMni6KBAe4ME5cUiv2ee5qZGQkIDU1FQIgoCgoCCEhoY29q0shuGGiKjlSLlShAmfx6GkvAoKexnKKnX611RKBRaODcaoEJWEFVJL0ezz3OTl5eHee+9F//79MXv2bMyaNQv9+vXDsGHDcPXq1UYVTURErU+wj6v+Dqq/BxsAyNGUYea6E9iVrJaiNLJijQo3zz//PIqLi3HmzBkUFhbi2rVrSE5ORlFREWbPnt3UNRIRkY3S6kR8fzzL6Gs13QqLdqRAq2uyERTUCjRq/YRdu3Zh9+7dCAoK0m8LDg7GihUrEBUV1WTFERGRbYvPKIRaU1bn6yIAtaYM8RmFiOjaznKFkVVr1JUbnU5ndNCwvb29fv4bIiKi+uQV1x1sGtOOCGhkuBk6dCheeOEFXLlyRb8tOzsbc+fOxbBhw5qsOCIism2eLoombUcENDLcfPrppyguLoa/vz+6du2Kbt26ISAgAMXFxVi+fHlT10hERDYqPMADKqUCpm74VimrbwsnMlejxtz4+vrixIkTiI2NxdmzZyGKIoKDgzF8+PCmro+IiGyYXCZg4dhgzFx3AgL+GkT8d5zvhhqqQVdu9uzZg+DgYBQVVa8FMmLECDz//POYPXs2+vfvj169euHgwYPNUigREdmmUSEqxEwOhbfSsOtJ6VQ9tvOHk1ew92yeFKWRlWrQJH4PPPAA7r33XsydO9fo65988gn27t2LrVu3NlmBTY2T+BERtUxanYj4jEKDGYpf356M/x3NhIvCDj/MugsB7dtIXSZJpNkm8Tt58iRGjRpV5+tRUVFISEhoyFsSEREBqO6iiujaDg/27YiIru0glwl4c2wvhPm5o7isCv/8+jhKyqukLpOsQIPCTW5ursl1o+zs7DhDMRERNRkHOxliHguFl6sjLuSV4MXvk6DjhH5UjwaFm44dO+L06dN1vn7q1CmoVFwDhIiImo6nqwIxk8PgIJfhlzO5WLH3D6lLohauQeFmzJgxeOONN1BWVnsypZs3b2LhwoW4//77m6w4IiIiAAjt7I7FD/YCACzdfR57zuZKXBG1ZA0aUJybm4vQ0FDI5XLMmjULgYGBEAQBqampWLFiBbRaLU6cOAEvL6/mrPm2cEAxEZH1em3baaz7PRMujnbYPmsQunRoK3VJZCEN+f1uULgBgEuXLmHmzJn45ZdfULOrIAgYOXIkVq5cCX9//0YXbgkMN0RE1quiSofHvvwdxy5eQzfPttj6bCRcFHWPBSXb0azhpsa1a9fwxx9/QBRFdO/eHe7u7o0q1tIYboiIrFtecRkeWH4YOUVliAr2wmeTwyDjJH82r9luBf87d3d39O/fH+Hh4VYTbIiIyPp5uigQMzkUDnIZfk3JxaccYEy3aHS4ISIiksqdnd3x73EhAICPdp/Hb6kcYEx/YbghIiKrNKG/L6YM9IMoAnO+S0La1RKpS6IWguGGiIis1uv3B6O/vzuKy6tnMC4uq5S6JGoBGG6IiMhqOdjJsPKxMHi7KpB29Qaivz/JGYyJ4YaIiKxbBxdHfDYlDA52MsSm5GL5Hg4wbu0YboiIyOr19XUzGGAcm8IBxq0Zww0REdmECf18MTXCDwAwd0MS/sjjAOPWiuGGiIhsxuv3ByPc3wMl5VX45zfHUcQBxq0Sww0REdkMe7kMKx4LhUqpQPrVG4jekMQBxq0Qww0REdmUDi6O+Gxy9QDj3al5+Pi3C9DqRMSlFWB7Ujbi0gqgZeCxaY1eW8pacW0pIqLWYePxLLy86RQAwN3ZHtdK/+qiUikVWDg2GKNCVFKVRw1kkbWliIiIWrKH+/ninh4dAMAg2ABAjqYMM9edwK5ktRSlUTNjuCEiIpuk1Yk4m1Ns9LWaLotFO1LYRWWDGG6IiMgmxWcUIqeorM7XRQBqTRniMwotVxRZBMMNERHZpLziuoNNY9qR9WC4ISIim+TpomjSdmQ9GG6IiMgmhQd4QKVUQKjjdQHVd02FB3hYsiyyAIYbIiKySXKZgIVjgwHAaMARASwcGwy5rK74Q9aK4YaIiGzWqBAVYiaHwltpvOuprFJn4YrIEiQPNytXrkRAQAAUCgXCwsJw8OBBk+3Ly8uxYMEC+Pn5wdHREV27dsXq1astVC0REVmbUSEqHHp1KL59eiA+ntQX3z49EM8P7QYAWLD1NC7m35C4QmpqdlJ++IYNGzBnzhysXLkSgwYNwueff47Ro0cjJSUFnTt3NrrPhAkTkJubi1WrVqFbt27Iy8tDVVWVhSsnIiJrIpcJiOjaTv+8v787jmYUIj6jELO+PYHNMyPhaCeXsEJqSpIuvzBgwACEhoYiJiZGvy0oKAjjxo3DkiVLarXftWsXJk2ahPT0dHh4NG4AGJdfICIiAFBrbmL0xwdxvbQST94VgNfvD5a6JDLBKpZfqKioQEJCAqKiogy2R0VF4ciRI0b3+eGHH9CvXz+899576NixI3r06IGXXnoJN2/erPNzysvLUVRUZPAgIiJSKZ3wwf/1AQCsOpSB31JzJa6Imopk4SY/Px9arRZeXl4G2728vJCTk2N0n/T0dBw6dAjJycnYunUrli1bhk2bNuG5556r83OWLFkCpVKpf/j6+jbpcRARkfUaHuyFxwf5AwBe2ngSORpO6GcLJB9QLAiGt+CJolhrWw2dTgdBELB+/XqEh4djzJgxWLp0KdasWVPn1Zv58+dDo9HoH1lZWU1+DEREZL3mje6JXj6uuFZaiRe+S+RaUzZAsnDTvn17yOXyWldp8vLyal3NqaFSqdCxY0colUr9tqCgIIiiiMuXLxvdx9HREa6urgYPIiKiGo52cix/5E44O8hxNKMQK/b+IXVJdJskCzcODg4ICwtDbGyswfbY2FhERkYa3WfQoEG4cuUKSkpK9NvOnz8PmUyGTp06NWu9RERku7p0aIt/jwsBACzbfZ6LaVo5SbuloqOj8eWXX2L16tVITU3F3LlzkZmZiRkzZgCo7lKaOnWqvv2jjz6Kdu3a4fHHH0dKSgoOHDiAl19+GU888QScnJykOgwiIrIB40M7YXxoR+hE4IXvEnHtRoXUJVEjSRpuJk6ciGXLlmHx4sXo27cvDhw4gJ07d8LPzw8AoFarkZmZqW/ftm1bxMbG4vr16+jXrx8ee+wxjB07Fp988olUh0BERDbkrQdDENC+DdSaMry86RQknC2FboOk89xIgfPcEBGRKcnZGoxfeQQVWh3eHBuM6YMCpC6JYCXz3BAREbVEIR2V+NeYngCAt3eeRXK2RuKKqKEYboiIiG4xLdIfw4O8UKHVYfa3ibhRzmV+rAnDDRER0S0EQcD7/3cHVEoF0vNv4I3tZ6QuiRqA4YaIiMgI9zYOWDaxL2QCsPnEZWw5YXw+NWp5GG6IiIjqMKBLO7wwrAcA4LVtyUi/WlLPHtQSMNwQERGZMGtoNwwI8EBphRbPf5uI8iqt1CVRPRhuiIiITJDLBHw86U64O9vjzJUivPvzOalLonow3BAREdXDW6nAhxP6AABWH87A7pRciSsiUxhuiIiIzDC0pxeevKt6Qr+XN52EWnNT4oqoLgw3REREZnplVCB6d1TiWmklXvguCVpdq5rk32ow3BAREZnJ0U6O5Y/ciTYOcsRnFOLj384jLq0A25OyEZdWwLDTQnBtKSIiogbalpiNORuSam1XKRVYODYYo0JUli/KxnFtKSIiomaksDf+85mjKcPMdSewK1lt4Yro7xhuiIiIGkCrE7FoR4rR12q6QhbtSGEXlYQYboiIiBogPqMQak1Zna+LANSaMsRnFFquKDLAcENERNQAecV1B5vGtKOmx3BDRETUAJ4uiiZtR02P4YaIiKgBwgM8oFIqIJhoo1IqEB7gYbGayBDDDRERUQPIZQIWjg0GgDoDztODu0AuMxV/qDkx3BARETXQqBAVYiaHwltp2PXkaFf9s7otKRuVWp0UpREAO6kLICIiskajQlQYEeyN+IxC5BWXwdNFAb92zhj98UGcuqxBzL40zB7WXeoyWyVeuSEiImokuUxARNd2eLBvR0R0bQcfNycsfrAXAOCT3y4gOVsjcYWtE8MNERFRE3qgjw/G9PZGlU5E9PdJKK/SSl1Sq8NwQ0RE1IQEQcBbD4agfVsHnM8twUexF6QuqdVhuCEiImpi7do64u1/9AYAfHEgDQmXOFuxJTHcEBERNYOoXt54KLQTdCLw4vcnUVpRJXVJrQbDDRERUTN5Y2wwVEoFLhaU4t2fz0pdTqvBcENERNRMlE72eO//7gAArI27hMN/5EtcUevAcENERNSMBnfvgCkD/QAAL288iaKySokrsn0MN0RERM1s3uie8GvnjCuaMry1I0Xqcmweww0REVEza+Nohw8e7gNBADYmXMbulFypS7JpDDdEREQW0N/fA/8c3AUAMG/LaRTeqJC4ItvFcENERGQhc0f0QHfPtsgvKcfr25OlLsdmMdwQERFZiMJejqUT+sJOJuCnU2rsOHlF6pJsEsMNERGRBfXupMSsod0AAK9vT0ZeUZnEFdkehhsiIiILe+7ebgjp6IrrpZWYt+U0RFGUuiSbwnBDRERkYfZyGZZO6AsHuQx7zuZh4/HLUpdkUxhuiIiIJNDDywUvRvUAACz+MQVZhaUSV2Q7GG6IiIgk8tTgLujn546S8iq8sukUdDp2TzUFhhsiIiKJyGUCPni4D5zs5YhLL8DXcRelLskmMNwQERFJyL99G/xrTE8AwDu7ziL9aonEFVk/hhsiIiKJTR7oh8Hd26OsUocXN55ElVYndUlWjeGGiIhIYoIg4N2H7oCLox0SM6/ji4PpUpdk1RhuiIiIWgAfNycsfKAXAOCj2PNIztYgLq0A25OyEZdWAC0HG5vNTuoCiIiIqNpDoR3xy5kcxKbkYtyKw6j6W6BRKRVYODYYo0JUElZoHSS/crNy5UoEBARAoVAgLCwMBw8erLPtvn37IAhCrcfZs2ctWDEREVHzEAQBw4O8AMAg2ABAjqYMM9edwK5ktRSlWRVJw82GDRswZ84cLFiwAImJiRg8eDBGjx6NzMxMk/udO3cOarVa/+jevbuFKiYiImo+Wp2IZbvPG32tJuos2pHCLqp6SBpuli5diieffBJPPfUUgoKCsGzZMvj6+iImJsbkfp6envD29tY/5HK5hSomIiJqPvEZhVBr6l5IUwSg1pQhPqPQckVZIcnCTUVFBRISEhAVFWWwPSoqCkeOHDG575133gmVSoVhw4Zh7969JtuWl5ejqKjI4EFERNQS5RWbt0K4ue1aK8nCTX5+PrRaLby8vAy2e3l5IScnx+g+KpUKX3zxBTZv3owtW7YgMDAQw4YNw4EDB+r8nCVLlkCpVOofvr6+TXocRERETcXTRdGk7Vorye+WEgTB4LkoirW21QgMDERgYKD+eUREBLKysvDBBx/g7rvvNrrP/PnzER0drX9eVFTEgENERC1SeIAHVEoFcjRlMDaqRgDgrVQgPMDD0qVZFcmu3LRv3x5yubzWVZq8vLxaV3NMGThwIC5cuFDn646OjnB1dTV4EBERtURymYCFY4MBVAeZW4kAFo4Nhlxm/CIAVZMs3Dg4OCAsLAyxsbEG22NjYxEZGWn2+yQmJkKl4j3/RERkG0aFqBAzORTeytpdT452MtzRyc3yRVkZSbuloqOjMWXKFPTr1w8RERH44osvkJmZiRkzZgCo7lLKzs7G119/DQBYtmwZ/P390atXL1RUVGDdunXYvHkzNm/eLOVhEBERNalRISqMCPZGfEYh8orL0KGtIz749RxOZF7H69uS8eW0fnUO4SCJw83EiRNRUFCAxYsXQ61WIyQkBDt37oSfnx8AQK1WG8x5U1FRgZdeegnZ2dlwcnJCr1698NNPP2HMmDFSHQIREVGzkMsERHRtp3/+rosj7vvkEH47m4cfT6kxto+PhNW1bIIoiq1qJqCioiIolUpoNBqOvyEiIqvy8e4L+Gj3ebRr44Dd0UPg3sZB6pIspiG/35Ivv0BERETmmXlPV/TwaouCGxX4z85UqctpsRhuiIiIrISDnQxLxt8BQQA2JVzGoQv5UpfUIjHcEBERWZEwP3dMi/AHAMzfego3K7TSFtQCMdwQERFZmZdGBsJHqUBW4U18VMdCm60Zww0REZGVaetoh7fGhQAAvjyYjtOXNRJX1LIw3BAREVmhYUFeGNvHBzoReHXzKVRqdVKX1GIw3BAREVmphWOD4eZsjxR1EVYdypC6nBaD4YaIiMhKtW/riNfuq16L6qPY87iYf0PiiloGhhsiIiIr9lBoR9zVrT3Kq3SYv+U0WtncvEYx3BAREVkxQRDw9j96Q2EvQ1x6ATYevyx1SZJjuCEiIrJynds548URgQCAf/+UgrziMokrkhbDDRERkQ14fJA/endUoqisCot+SJG6HEkx3BAREdkAO7kM7zzUG3KZgJ9Oq/HrmRypS5IMww0REZGN6OWjxNODuwAAXt+ejKKySokrkgbDDRERkQ2ZM7w7/No5I7eoHO/tOit1OZJguCEiIrIhCns5lozvDQBY93smjl8slLgiy2O4ISIisjGRXdtjYj9fANVLM5RXta6VwxluiIiIbNC/xgShfVtHpF29gRV706Qux6IYboiIiGyQ0tkeix/sBQCI2fcHzuUUS1yR5TDcEBER2ajRId4YEeyFSq2IVzefglbXOpZmYLghIiKyUYIg4K0HQ+DiaIekrOv4Ju6i1CVZBMMNERGRDfNWKvDK6J4AgPd+OYfs6zclrqj5MdwQERHZuMfCO6OfnztKK7RYsOUU4tLysT0pG3FpBTbZVSWIrWxt9KKiIiiVSmg0Gri6ukpdDhERkUX8kVeMUcsOouqWMKNSKrBwbDBGhagkqsw8Dfn95pUbIiKiVuCPvJJawQYAcjRlmLnuBHYlqyWoqnkw3BAREdk4rU7Eoh3GVwqviTuLdqTYTBcVww0REZGNi88ohFpTVufrIgC1pgzxGbaxVAPDDRERkY3LK6472DSmXUvHcENERGTjPF0UTdqupWO4ISIisnHhAR5QKRUQTLRRKRUID/CwWE3NieGGiIjIxsllAhaODQaAOgPOG/cHQy4zFX+sB8MNERFRKzAqRIWYyaHwVhrvepLZSLABADupCyAiIiLLGBWiwohgb8RnFCKvuAyeLgocuJCHmH3pWPTDGQzu3h7ODtYfDaz/CIiIiMhscpmAiK7t9M/7+rrhhyQ1sq/fxCe//YF5f65DZc3YLUVERNSKOTnIseiBXgCALw+m40JuscQV3T6GGyIiolZueLAXhgd5oUon4rVtybD2ZScZboiIiAgLxwZDYS/D0YxCbEvKlrqc28JwQ0RERPD1cMbzQ7sDAP7zUyo0NyslrqjxGG6IiIgIAPD04C7o2qEN8ksq8OGv56Qup9EYboiIiAgA4GAnw1sPhgAAvvn9Ek5f1khcUeMw3BAREZFeZLf2eLCvD0QReG3baWh11je4mOGGiIiIDCy4LwgujnY4eVmDb+MzpS6nwRhuiIiIyICniwIvjQwEALy36yzyS8olrqhhGG6IiIiolskD/dDLxxVFZVVYsvOs1OU0CMMNERER1SKXCfj3uBAIArD5xGUcTS+QuiSzSR5uVq5ciYCAACgUCoSFheHgwYNm7Xf48GHY2dmhb9++zVsgERFRK3VnZ3c8Et4ZAPD69mRUanUSV2QeScPNhg0bMGfOHCxYsACJiYkYPHgwRo8ejcxM04OXNBoNpk6dimHDhlmoUiIiotbplZGB8GjjgPO5JfjqcIbU5ZhF0nCzdOlSPPnkk3jqqacQFBSEZcuWwdfXFzExMSb3e+aZZ/Doo48iIiLCQpUSERG1Tm7ODpj/50rhy3ZfwJXrNyWuqH6ShZuKigokJCQgKirKYHtUVBSOHDlS535fffUV0tLSsHDhQrM+p7y8HEVFRQYPIiIiMt9DoZ3Q398dpRVaLN6RInU59ZIs3OTn50Or1cLLy8tgu5eXF3Jycozuc+HCBcybNw/r16+HnZ2dWZ+zZMkSKJVK/cPX1/e2ayciImpNZDIBb40LgVwmYNeZHOw9myd1SSZJPqBYEASD56Io1toGAFqtFo8++igWLVqEHj16mP3+8+fPh0aj0T+ysrJuu2YiIqLWpqe3K54Y5A8AWPjDGZRVaqUtyATJwk379u0hl8trXaXJy8urdTUHAIqLi3H8+HHMmjULdnZ2sLOzw+LFi3Hy5EnY2dlhz549Rj/H0dERrq6uBg8iIiJquBeG94C3qwKZhaVYuS9N6nLqJFm4cXBwQFhYGGJjYw22x8bGIjIyslZ7V1dXnD59GklJSfrHjBkzEBgYiKSkJAwYMMBSpRMREbVKbR3t8MbYYADAZ/vSkJF/Q+KKjDNv4EoziY6OxpQpU9CvXz9ERETgiy++QGZmJmbMmAGgukspOzsbX3/9NWQyGUJCQgz29/T0hEKhqLWdiIiImsfoEG/c3aMDDpy/ije2J+PrJ8KNDieRkqThZuLEiSgoKMDixYuhVqsREhKCnTt3ws/PDwCgVqvrnfOGiIiILEcQBCx+oBeilh3AwQv52Hk6B/fdoZK6LAOCKIrWt5b5bSgqKoJSqYRGo+H4GyIiokb6KPY8Pv7tArxcHfHbi/egrWPzXi9pyO+35HdLERERkfWZeU9X+LVzRm5ROZbFnpe6HAMMN0RERNRgCns5Fj3QCwDw1ZGLSFW3nElyGW6IiIioUe4J9MSY3t7Q6kS8ti0ZOl3LGOki6YBiIiIism6v3x+MfeeuIuHSNXyfkAU/jzbIKy6Dp4sC4QEekMssfycVww0RERE1mkrphLnDe+A/O1Mxf8tp/P02JZVSgYVjgzEqxLJ3U7FbioiIiG6Lyk0BALj1/uscTRlmrjuBXclqi9bDcENERESNptWJ+M9PqUZfq8k6i3akQGvB8TgMN0RERNRo8RmFUGvK6nxdBKDWlCE+o9BiNTHcEBERUaPlFdcdbBrTrikw3BAREVGjeboomrRdU2C4ISIiokYLD/CASqlAXTd8C6i+ayo8wMNiNTHcEBERUaPJZQIWjg0GgFoBp+b5wrHBFp3vhuGGiIiIbsuoEBViJofCW2nY9eStVCBmcqjF57nhJH5ERER020aFqDAi2BvxGYWcoZiIiIhsg1wmIKJrO6nLYLcUERER2RaGGyIiIrIpDDdERERkUxhuiIiIyKYw3BAREZFNYbghIiIim8JwQ0RERDaF4YaIiIhsCsMNERER2ZRWN0OxKIoAgKKiIokrISIiInPV/G7X/I6b0urCTXFxMQDA19dX4kqIiIiooYqLi6FUKk22EURzIpAN0el0uHLlClxcXCAITbuYV1FREXx9fZGVlQVXV9cmfe+WwNaPD7D9Y+TxWT9bP0Yen/VrrmMURRHFxcXw8fGBTGZ6VE2ru3Ijk8nQqVOnZv0MV1dXm/2fFrD94wNs/xh5fNbP1o+Rx2f9muMY67tiU4MDiomIiMimMNwQERGRTWG4aUKOjo5YuHAhHB0dpS6lWdj68QG2f4w8Putn68fI47N+LeEYW92AYiIiIrJtvHJDRERENoXhhoiIiGwKww0RERHZFIYbIiIisikMNw20cuVKBAQEQKFQICwsDAcPHjTZfv/+/QgLC4NCoUCXLl3w2WefWajShlmyZAn69+8PFxcXeHp6Yty4cTh37pzJffbt2wdBEGo9zp49a6GqG+bNN9+sVau3t7fJfazl/AGAv7+/0fPx3HPPGW3f0s/fgQMHMHbsWPj4+EAQBGzbts3gdVEU8eabb8LHxwdOTk645557cObMmXrfd/PmzQgODoajoyOCg4OxdevWZjqC+pk6xsrKSrz66qvo3bs32rRpAx8fH0ydOhVXrlwx+Z5r1qwxel7Lysqa+Whqq+8cTp8+vVadAwcOrPd9reUcAjB6LgRBwPvvv1/ne7aUc2jO70JL/R4y3DTAhg0bMGfOHCxYsACJiYkYPHgwRo8ejczMTKPtMzIyMGbMGAwePBiJiYn417/+hdmzZ2Pz5s0Wrrx++/fvx3PPPYfff/8dsbGxqKqqQlRUFG7cuFHvvufOnYNardY/unfvboGKG6dXr14GtZ4+fbrOttZ0/gDg2LFjBscWGxsLAHj44YdN7tdSz9+NGzfQp08ffPrpp0Zff++997B06VJ8+umnOHbsGLy9vTFixAj9+nHGxMXFYeLEiZgyZQpOnjyJKVOmYMKECTh69GhzHYZJpo6xtLQUJ06cwOuvv44TJ05gy5YtOH/+PB544IF639fV1dXgnKrVaigUiuY4BJPqO4cAMGrUKIM6d+7cafI9rekcAqh1HlavXg1BEPDQQw+ZfN+WcA7N+V1osd9DkcwWHh4uzpgxw2Bbz549xXnz5hlt/8orr4g9e/Y02PbMM8+IAwcObLYam0peXp4IQNy/f3+dbfbu3SsCEK9du2a5wm7DwoULxT59+pjd3prPnyiK4gsvvCB27dpV1Ol0Rl+3pvMHQNy6dav+uU6nE729vcV33nlHv62srExUKpXiZ599Vuf7TJgwQRw1apTBtpEjR4qTJk1q8pob6tZjNCY+Pl4EIF66dKnONl999ZWoVCqbtrgmYOz4pk2bJj744IMNeh9rP4cPPvigOHToUJNtWuo5vPV3oSV/D3nlxkwVFRVISEhAVFSUwfaoqCgcOXLE6D5xcXG12o8cORLHjx9HZWVls9XaFDQaDQDAw8Oj3rZ33nknVCoVhg0bhr179zZ3abflwoUL8PHxQUBAACZNmoT09PQ621rz+auoqMC6devwxBNP1LtArDWdvxoZGRnIyckxOD+Ojo4YMmRInd9HoO5zamqflkSj0UAQBLi5uZlsV1JSAj8/P3Tq1An3338/EhMTLVNgI+zbtw+enp7o0aMHnn76aeTl5Zlsb83nMDc3Fz/99BOefPLJetu2xHN46+9CS/4eMtyYKT8/H1qtFl5eXgbbvby8kJOTY3SfnJwco+2rqqqQn5/fbLXeLlEUER0djbvuugshISF1tlOpVPjiiy+wefNmbNmyBYGBgRg2bBgOHDhgwWrNN2DAAHz99df45Zdf8N///hc5OTmIjIxEQUGB0fbWev4AYNu2bbh+/TqmT59eZxtrO39/V/Oda8j3sWa/hu7TUpSVlWHevHl49NFHTS5G2LNnT6xZswY//PADvv32WygUCgwaNAgXLlywYLXmGT16NNavX489e/bgww8/xLFjxzB06FCUl5fXuY81n8O1a9fCxcUF48ePN9muJZ5DY78LLfl72OpWBb9dt/4tWBRFk38zNtbe2PaWZNasWTh16hQOHTpksl1gYCACAwP1zyMiIpCVlYUPPvgAd999d3OX2WCjR4/W/3vv3r0RERGBrl27Yu3atYiOjja6jzWePwBYtWoVRo8eDR8fnzrbWNv5M6ah38fG7iO1yspKTJo0CTqdDitXrjTZduDAgQaDcgcNGoTQ0FAsX74cn3zySXOX2iATJ07U/3tISAj69esHPz8//PTTTyYDgDWeQwBYvXo1HnvssXrHzrTEc2jqd6Elfg955cZM7du3h1wur5Us8/LyaiXQGt7e3kbb29nZoV27ds1W6+14/vnn8cMPP2Dv3r3o1KlTg/cfOHBgi/wbojFt2rRB796966zXGs8fAFy6dAm7d+/GU0891eB9reX81dzl1pDvY81+Dd1HapWVlZgwYQIyMjIQGxtr8qqNMTKZDP3797eK86pSqeDn52eyVms8hwBw8OBBnDt3rlHfS6nPYV2/Cy35e8hwYyYHBweEhYXp70CpERsbi8jISKP7RERE1Gr/66+/ol+/frC3t2+2WhtDFEXMmjULW7ZswZ49exAQENCo90lMTIRKpWri6ppHeXk5UlNT66zXms7f33311Vfw9PTEfffd1+B9reX8BQQEwNvb2+D8VFRUYP/+/XV+H4G6z6mpfaRUE2wuXLiA3bt3NypUi6KIpKQkqzivBQUFyMrKMlmrtZ3DGqtWrUJYWBj69OnT4H2lOof1/S606O9hkw1NbgW+++470d7eXly1apWYkpIizpkzR2zTpo148eJFURRFcd68eeKUKVP07dPT00VnZ2dx7ty5YkpKirhq1SrR3t5e3LRpk1SHUKeZM2eKSqVS3Ldvn6hWq/WP0tJSfZtbj++jjz4St27dKp4/f15MTk4W582bJwIQN2/eLMUh1OvFF18U9+3bJ6anp4u///67eP/994suLi42cf5qaLVasXPnzuKrr75a6zVrO3/FxcViYmKimJiYKAIQly5dKiYmJurvFHrnnXdEpVIpbtmyRTx9+rT4yCOPiCqVSiwqKtK/x5QpUwzuZjx8+LAol8vFd955R0xNTRXfeecd0c7OTvz9998tfnyiaPoYKysrxQceeEDs1KmTmJSUZPC9LC8v17/Hrcf45ptvirt27RLT0tLExMRE8fHHHxft7OzEo0ePtqjjKy4uFl988UXxyJEjYkZGhrh3714xIiJC7Nixo82cwxoajUZ0dnYWY2JijL5HSz2H5vwutNTvIcNNA61YsUL08/MTHRwcxNDQUINbpadNmyYOGTLEoP2+ffvEO++8U3RwcBD9/f3r/J9bagCMPr766it9m1uP79133xW7du0qKhQK0d3dXbzrrrvEn376yfLFm2nixImiSqUS7e3tRR8fH3H8+PHimTNn9K9b8/mr8csvv4gAxHPnztV6zdrOX82t6rc+pk2bJopi9W2oCxcuFL29vUVHR0fx7rvvFk+fPm3wHkOGDNG3r7Fx40YxMDBQtLe3F3v27ClpmDN1jBkZGXV+L/fu3at/j1uPcc6cOWLnzp1FBwcHsUOHDmJUVJR45MgRyx+caPr4SktLxaioKLFDhw6ivb292LlzZ3HatGliZmamwXtY8zms8fnnn4tOTk7i9evXjb5HSz2H5vwutNTvofDnARARERHZBI65ISIiIpvCcENEREQ2heGGiIiIbArDDREREdkUhhsiIiKyKQw3REREZFMYboiIiMimMNwQUavg7++PZcuWSV0GEVkAww0RNbnp06dj3LhxAIB77rkHc+bMsdhnr1mzBm5ubrW2Hzt2DP/85z8tVgcRScdO6gKIiMxRUVEBBweHRu/foUOHJqyGiFoyXrkhomYzffp07N+/Hx9//DEEQYAgCLh48SIAICUlBWPGjEHbtm3h5eWFKVOmID8/X7/vPffcg1mzZiE6Ohrt27fHiBEjAABLly5F79690aZNG/j6+uLZZ59FSUkJAGDfvn14/PHHodFo9J/35ptvAqjdLZWZmYkHH3wQbdu2haurKyZMmIDc3Fz962+++Sb69u2Lb775Bv7+/lAqlZg0aRKKi4v1bTZt2oTevXvDyckJ7dq1w/Dhw3Hjxo1m+q9JROZiuCGiZvPxxx8jIiICTz/9NNRqNdRqNXx9faFWqzFkyBD07dsXx48fx65du5Cbm4sJEyYY7L927VrY2dnh8OHD+PzzzwEAMpkMn3zyCZKTk7F27Vrs2bMHr7zyCgAgMjISy5Ytg6urq/7zXnrppVp1iaKIcePGobCwEPv370dsbCzS0tIwceJEg3ZpaWnYtm0bfvzxR/z444/Yv38/3nnnHQCAWq3GI488gieeeAKpqanYt28fxo8fDy7XRyQ9dksRUbNRKpVwcHCAs7MzvL299dtjYmIQGhqKt99+W79t9erV8PX1xfnz59GjRw8AQLdu3fDee+8ZvOffx+8EBATgrbfewsyZM7Fy5Uo4ODhAqVRCEASDz7vV7t27cerUKWRkZMDX1xcA8M0336BXr144duwY+vfvDwDQ6XRYs2YNXFxcAABTpkzBb7/9hv/85z9Qq9WoqqrC+PHj4efnBwDo3bv3bfzXIqKmwis3RGRxCQkJ2Lt3L9q2bat/9OzZE0D11ZIa/fr1q7Xv3r17MWLECHTs2BEuLi6YOnUqCgoKGtQdlJqaCl9fX32wAYDg4GC4ubkhNTVVv83f318fbABApVIhLy8PANCnTx8MGzYMvXv3xsMPP4z//ve/uHbtmvn/EYio2TDcEJHF6XQ6jB07FklJSQaPCxcu4O6779a3a9OmjcF+ly5dwpgxYxASEoLNmzcjISEBK1asAABUVlaa/fmiKEIQhHq329vbG7wuCAJ0Oh0AQC6XIzY2Fj///DOCg4OxfPlyBAYGIiMjw+w6iKh5MNwQUbNycHCAVqs12BYaGoozZ87A398f3bp1M3jcGmj+7vjx46iqqsKHH36IgQMHokePHrhy5Uq9n3er4OBgZGZmIisrS78tJSUFGo0GQUFBZh+bIAgYNGgQFi1ahMTERDg4OGDr1q1m709EzYPhhoialb+/P44ePYqLFy8iPz8fOp0Ozz33HAoLC/HII48gPj4e6enp+PXXX/HEE0+YDCZdu3ZFVVUVli9fjvT0dHzzzTf47LPPan1eSUkJfvvtN+Tn56O0tLTW+wwfPhx33HEHHnvsMZw4cQLx8fGYOnUqhgwZYrQrzJijR4/i7bffxvHjx5GZmYktW7bg6tWrDQpHRNQ8GG6IqFm99NJLkMvlCA4ORocOHZCZmQkfHx8cPnwYWq0WI0eOREhICF544QUolUrIZHX/sdS3b18sXboU7777LkJCQrB+/XosWbLEoE1kZCRmzJiBiRMnokOHDrUGJAPVV1y2bdsGd3d33H333Rg+fDi6dOmCDRs2mH1crq6uOHDgAMaMGYMePXrgtddew4cffojRo0eb/x+HiJqFIPK+RSIiIrIhvHJDRERENoXhhoiIiGwKww0RERHZFIYbIiIisikMN0RERGRTGG6IiIjIpjDcEBERkU1huCEiIiKbwnBDRERENoXhhoiIiGwKww0RERHZFIYbIiIisin/DwPsQI/yMsJ2AAAAAElFTkSuQmCC\n"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "iterations = 20\n",
    "\n",
    "costs = []\n",
    "\n",
    "for i in range(iterations):\n",
    "    params, cost = opt.step_and_cost(circuit, params)\n",
    "    costs.append(cost)\n",
    "\n",
    "# Visualize results\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "costs.append(circuit(params))\n",
    "plt.plot(costs, \"-o\")\n",
    "plt.xlabel(\"Iterations\")\n",
    "plt.ylabel(\"Cost\")\n",
    "\n",
    "print(\"Minimized circuit output:\", circuit(params))\n",
    "print(\"Optimized parameters:\", params)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<div class=\"alert alert-block alert-info\">\n",
    "<b>Note</b> The circuit considered here is very simple and can be optimized easily by hand. However, the need for PennyLane's automatic differentiation capabilities becomes apparent as we make the problem more complicated, e.g., with more gates and different types of output measurement. In later demos, we will also see how Braket can be used to parallelize evaluation of the gradient, providing a turbocharger for quantum circuit training in PennyLane.\n",
    "</div>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Running circuits on Braket's on-demand simulator, SV1"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "So far we have used the local Braket simulator. This is a great choice for quick prototyping, but it is not suitable for large circuits with many qubits and does not provide a connection to quantum hardware.\n",
    "\n",
    "Amazon Braket also provides access to on-demand, high-performance simulators and quantum processing units (QPUs) from different [providers](https://aws.amazon.com/braket/hardware-providers/). These devices can be accessed through PennyLane by changing a single line of code, unlocking the potential for machine learning and optimization on quantum hardware and high performance simulators!\n",
    "\n",
    "<div align=\"center\"><img src=\"remote-single-job.png\"/></div>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-02-19T21:42:11.066378Z",
     "start_time": "2024-02-19T21:42:11.064297Z"
    }
   },
   "outputs": [],
   "source": [
    "# Use Braket SDK Cost Tracking to estimate the cost to run this example\n",
    "from braket.tracking import Tracker\n",
    "\n",
    "braket_tasks_cost = Tracker().start()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Each remote Braket device can be selected through its [ARN](https://docs.aws.amazon.com/general/latest/gr/aws-arns-and-namespaces.html). The supported devices on Braket are listed [here](https://docs.aws.amazon.com/braket/latest/developerguide/braket-devices.html). For now, we will pick the on-demand SV1 simulator.\n",
    "\n",
    "<div class=\"alert alert-block alert-warning\">\n",
    "<b>Caution:</b> Running hybrid algorithms on a QPU can take a long time and incur high <a href=\"https://aws.amazon.com/braket/pricing/\">usage fees</a> charged to your AWS account.\n",
    "</div>"
   ]
  },
  {
   "cell_type": "markdown",
   "source": [
    "In PennyLane, all remote Braket devices are then accessed through a single PennyLane device named ``braket.aws.qubit``."
   ],
   "metadata": {
    "collapsed": false
   }
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-02-19T21:42:11.437758Z",
     "start_time": "2024-02-19T21:42:11.067501Z"
    }
   },
   "outputs": [],
   "source": [
    "from braket.devices import Devices\n",
    "dev = qml.device(\"braket.aws.qubit\", device_arn=Devices.Amazon.SV1, wires=2)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "A follow up [tutorial](./1_Parallelized_optimization_of_quantum_circuits.ipynb) shows you how to use the remote device to run multiple circuits in parallel, while the [QAOA tutorial](./2_Graph_optimization_with_QAOA.ipynb) takes a deeper dive into graph optimization, including using SV1 to optimize a 20-node graph."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's execute our circuit on SV1, as well as calculating the gradient:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-02-19T21:42:18.632136Z",
     "start_time": "2024-02-19T21:42:11.439272Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Result of circuit run on SV1: 0.37261070647126565\n",
      "Result of gradient calculation on SV1: [-0.19585986 -0.80291741]\n"
     ]
    }
   ],
   "source": [
    "@qml.qnode(dev)\n",
    "def circuit(params):\n",
    "    qml.RX(params[0], wires=0)\n",
    "    qml.RY(params[1], wires=1)\n",
    "    qml.CNOT(wires=[0, 1])\n",
    "    return qml.expval(qml.PauliZ(1))\n",
    "\n",
    "\n",
    "dcircuit = qml.grad(circuit)\n",
    "\n",
    "print(\"Result of circuit run on SV1:\", circuit(params))\n",
    "print(\"Result of gradient calculation on SV1:\", dcircuit(params))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<div class=\"alert alert-block alert-info\">\n",
    "<b>What's next?</b> Check out the other tutorials in this folder to understand how Braket and PennyLane can be combined to solve a range of problems, from graph optimization to quantum chemistry.\n",
    "</div>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-02-19T21:42:19.098494Z",
     "start_time": "2024-02-19T21:42:18.633487Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Quantum Task Summary\n",
      "{<_Amazon.SV1: 'arn:aws:braket:::device/quantum-simulator/amazon/sv1'>: {'shots': 0, 'tasks': {'COMPLETED': 2}, 'execution_duration': datetime.timedelta(microseconds=299000), 'billed_execution_duration': datetime.timedelta(seconds=6)}}\n",
      "Note: Charges shown are estimates based on your Amazon Braket simulator and quantum processing unit (QPU) task usage. Estimated charges shown may differ from your actual charges. Estimated charges do not factor in any discounts or credits, and you may experience additional charges based on your use of other services such as Amazon Elastic Compute Cloud (Amazon EC2).\n",
      "Estimated cost to run this example: 0.008 USD\n"
     ]
    }
   ],
   "source": [
    "print(\"Quantum Task Summary\")\n",
    "print(braket_tasks_cost.quantum_tasks_statistics())\n",
    "print(\n",
    "    \"Note: Charges shown are estimates based on your Amazon Braket simulator and quantum processing unit (QPU) task usage. Estimated charges shown may differ from your actual charges. Estimated charges do not factor in any discounts or credits, and you may experience additional charges based on your use of other services such as Amazon Elastic Compute Cloud (Amazon EC2).\"\n",
    ")\n",
    "print(\n",
    "    f\"Estimated cost to run this example: {braket_tasks_cost.qpu_tasks_cost() + braket_tasks_cost.simulator_tasks_cost():.3f} USD\"\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Running on a QPU with Amazon Braket Hybrid Jobs\n",
    "\n",
    "In this notebook, the classical part of the algorithm was running locally. For longer-running algorithms or those requiring more intensive compute resources, it's recommended to dispatch the algorithm to Amazon Braket Hybrid Jobs, which fully manages the classical infrastructure, allowing you to focus on the algorithm. For example, you can train a larger circuit or increase the number of iterations.\n",
    "\n",
    "The second benefit of running the algorithm as a hybrid job is that for iterative algorithms that require repeated calls to a QPU, you retain priority for that QPU. Once your quantum tasks are created in the hybrid job, they run ahead of other tasks waiting in the regular quantum task queue. This is because hybrid jobs have a separate queue from standalone tasks, ensuring that only a single hybrid job can run on a QPU at a time. This means your algorithm will not be interrupted by other quantum tasks, so it will run more efficiently and predictably. However, hybrid jobs have a separate queue from standalone tasks, so only a single hybrid job can run on a QPU at a time. For a single quantum circuit or a batch of circuits, it's recommended to create quantum tasks instead of hybrid jobs. Only iterative algorithms benefit from QPU priority queuing.\n",
    "\n",
    "Note that hybrid jobs have at least a one-minute startup time since they create a containerized environment on Amazon EC2. So for very short workloads, there is likely no need to create a hybrid job.\n",
    "\n",
    "You can run your local Python code as an Amazon Braket hybrid job by annotating your code with the `@hybrid_job`` decorator, as shown in the following code example. Only Python 3.10 is supported by default. For custom Python versions, you can choose to use a custom container from [Amazon Elastic Container Registry (ECR)](https://aws.amazon.com/ecr/) (see [BYOC](https://docs.aws.amazon.com/braket/latest/developerguide/braket-jobs-byoc.html)).\n",
    "\n",
    "\n",
    "In the following code, we create a hybrid job for 10 iterations targeting Rigetti Aspen-M-3. Since we specified Aspen-M-3 as the device, this job will run once Aspen-M-3 is available and has no jobs running ahead of it."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-02-19T21:42:19.442699Z",
     "start_time": "2024-02-19T21:42:19.099758Z"
    }
   },
   "outputs": [],
   "source": [
    "from braket.jobs import hybrid_job\n",
    "from braket.jobs.metrics import log_metric\n",
    "\n",
    "device_arn = Devices.Amazon.SV1\n",
    "# device_arn = Devices.Rigetti.AspenM3\n",
    "\n",
    "\n",
    "@hybrid_job(device=device_arn)  # set priority QPU\n",
    "def qubit_rotation(stepsize=0.1, iterations=5):\n",
    "    task_tracker = Tracker().start()  # track Braket quantum tasks costs\n",
    "    dev = qml.device(\"braket.aws.qubit\", device_arn=device_arn.value, wires=2, shots=1_000)\n",
    "\n",
    "    params = np.array([0.1, 0.2])\n",
    "\n",
    "    @qml.qnode(dev)\n",
    "    def circuit(params):\n",
    "        qml.RX(params[0], wires=0)\n",
    "        qml.RY(params[1], wires=1)\n",
    "        qml.CNOT(wires=[0, 1])\n",
    "        return qml.expval(qml.PauliZ(1))\n",
    "\n",
    "    opt = qml.GradientDescentOptimizer(stepsize)\n",
    "\n",
    "    costs = []\n",
    "    for i in range(iterations):\n",
    "        params, cost = opt.step_and_cost(circuit, params)\n",
    "        costs.append(cost)\n",
    "\n",
    "        # Record the value of the cost function with each iteration\n",
    "        log_metric(metric_name=\"cost_function\", value=cost, iteration_number=i)\n",
    "\n",
    "        # Additionally, keep track of cost in USD for Braket tasks\n",
    "        braket_task_cost = float(\n",
    "            task_tracker.qpu_tasks_cost() + task_tracker.simulator_tasks_cost()\n",
    "        )\n",
    "        log_metric(metric_name=\"braket_cost\", value=braket_task_cost, iteration_number=i)\n",
    "\n",
    "    return {\"params\": params, \"costs\": costs, \"braket_tasks_cost\": braket_task_cost}"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "How long will it take the hybrid job to run? \n",
    "Let's first check if the device is currently available with `AwsDevice(device_arn).is_available()`. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-02-19T21:42:19.614223Z",
     "start_time": "2024-02-19T21:42:19.443623Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "True"
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from braket.aws import AwsDevice\n",
    "\n",
    "AwsDevice(device_arn).is_available"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Next, we check the hybrid job queue depth with `AwsDevice(device_arn).queue_depth().jobs`. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-02-19T21:42:19.958746Z",
     "start_time": "2024-02-19T21:42:19.615073Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "'0'"
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "AwsDevice(device_arn).queue_depth().jobs"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "If the device is available and there are no hybrid jobs currently running, then it should take about 5 minutes to complete. \n",
    "\n",
    "<div class=\"alert alert-block alert-warning\">\n",
    "<b>Caution:</b> Running the following cell will only run once the QPU is available. This may take a long time and will result in <a href=\"https://aws.amazon.com/braket/pricing/\">usage fees</a> charged to your AWS account. Only uncomment the cell if you are comfortable with the potential wait-time and costs. We recommend monitoring the Billing & Cost Management Dashboard on the AWS console and being aware that hybrid jobs involving a large number of qubits can be costly.\n",
    "</div>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-02-19T21:42:29.440989Z",
     "start_time": "2024-02-19T21:42:19.960142Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "AwsQuantumJob('arn':'arn:aws:braket:us-west-2:667256736152:job/qubit-rotation-1708378939963')\n"
     ]
    }
   ],
   "source": [
    "job = qubit_rotation(stepsize=0.2, iterations=20)\n",
    "print(job)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Once the hybrid job has completed, we can retrieve the results with `job.results()`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-02-19T21:49:52.814493Z",
     "start_time": "2024-02-19T21:42:29.442078Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{'params': tensor([0.236, 2.649], requires_grad=True), 'costs': [array(0.978), array(0.954), array(0.962), array(0.93), array(0.874), array(0.848), array(0.816), array(0.744), array(0.638), array(0.528), array(0.328), array(0.24), array(0.072), array(-0.034), array(-0.202), array(-0.32), array(-0.498), array(-0.616), array(-0.672), array(-0.752)], 'braket_tasks_cost': 0.375}\n",
      "CPU times: user 160 ms, sys: 16.6 ms, total: 177 ms\n",
      "Wall time: 7min 23s\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "results = job.result()\n",
    "print(results)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-02-19T21:50:26.336487Z",
     "start_time": "2024-02-19T21:49:52.820854Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Minimized circuit output: 0.37261070647126565\n",
      "Optimized parameters: [0.4839502  1.13630274]\n"
     ]
    },
    {
     "data": {
      "text/plain": "<Figure size 640x480 with 1 Axes>",
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAksAAAGwCAYAAAC5ACFFAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy88F64QAAAACXBIWXMAAA9hAAAPYQGoP6dpAABXEUlEQVR4nO3de1wU9f4/8NewXNYLLAICi6Kg5gXxBoqAaTdFvJsdRUvMMkpP5e30rezy83LOyaxTx9Q0M81Kj3pKzSwvYd4TRAVUBE0RhXQRAVlQ5LY7vz84bK3AsuCys4Ov5+MxjxMzn5l9TxNnX8xn5vMRRFEUQUREREQ1spO6ACIiIiJbxrBEREREZALDEhEREZEJDEtEREREJjAsEREREZnAsERERERkAsMSERERkQn2UhfQFOj1ely/fh3Ozs4QBEHqcoiIiMgMoiiiqKgIPj4+sLOr/f4Rw5IFXL9+Hb6+vlKXQURERA2QlZWFtm3b1rqdYckCnJ2dAVT+y3ZxcZG4GiIiIjJHYWEhfH19Dd/jtWFYsoCqrjcXFxeGJSIiIpmp6xEaPuBNREREZALDEhEREZEJDEtEREREJjAsEREREZnAsERERERkAsMSERERkQmyCkuHDx/GqFGj4OPjA0EQ8P3339e5z6FDhxAcHAylUokOHTrgs88+q9Zm69atCAgIgJOTEwICArB9+/ZGqJ6IiIjkSFZh6c6dO+jVqxdWrFhhVvuMjAwMHz4cAwcORFJSEt566y3MnDkTW7duNbSJi4tDVFQUoqOjcfr0aURHR2PChAk4fvx4Y50GERERyYggiqIodRENIQgCtm/fjrFjx9ba5o033sAPP/yAtLQ0w7rp06fj9OnTiIuLAwBERUWhsLAQu3fvNrSJjIxEq1atsGnTJrNqKSwshEqlglar5aCUREREMmHu97es7izVV1xcHCIiIozWDR06FCdPnkR5ebnJNseOHav1uKWlpSgsLDRaLE2nFxGXnocdydcQl54HnV6WmZaIiEj2mvR0J9nZ2fDy8jJa5+XlhYqKCuTm5kKtVtfaJjs7u9bjLl68GAsXLmyUmgFgT4oGC3emQqMtMaxTq5SYPyoAkYHqRvtcIiIiqq5J31kCqs/3UtXr+Of1NbUxNU/MvHnzoNVqDUtWVpbF6t2TosGMDYlGQQkAsrUlmLEhEXtSNBb7LCIiIqpbk76z5O3tXe0OUU5ODuzt7eHu7m6yzb13m/7MyckJTk5OFq9XpxexcGcqaupwEwEIABbuTMWQAG8o7ExP+ic1nV5EQkY+copK4OmsRIi/m83XTEREVJMmHZbCwsKwc+dOo3U///wz+vbtCwcHB0Ob2NhYzJkzx6hNeHi4VWsFgISM/Gp3lP5MBKDRlmBH8jWM7d0GdjYaPtiNSERETYmswtLt27dx6dIlw88ZGRlITk6Gm5sb2rVrh3nz5uHatWv4+uuvAVS++bZixQrMnTsXMTExiIuLw9q1a43ecps1axYGDRqEJUuWYMyYMdixYwf27duHo0ePWv38copqD0p/Nve/p/Hu9yno7qNCYBsVAtu4ILCNCh1btzT77k1j3fmp6ka89+5YVTfiqslBDExERCQrsgpLJ0+exGOPPWb4ee7cuQCAZ599FuvXr4dGo0FmZqZhu7+/P3bt2oU5c+bg008/hY+PD5YtW4annnrK0CY8PBybN2/GO++8g3fffRcdO3bEli1b0L9/f+ud2P94OivNauegEHCnTIeEK/lIuJJvWK90sEOAujI4Bf4vSD3k1RIOCuNH0xrrzk9T6kYkIiKqIttxlmyJpcZZ0ulFPLxkP7K1JTUGDgGAt0qJg689iit5xUi5pkXKdS1Srmlx7nohist01fZxVNihq9oZ3X1U6NFGhcKScizZfb7a8auiy713fnR6EYV3y6G9W47Cksr/1d4tR+HdCsM/V23LzLuDs9fqHkZhU0wowjq6m/3vhYiIqDGY+/3NsGQBlhyUsqobC4BRoKktzFTR60Vk5N2pDFDXtEi5VoiU61oUlVTU6/MdFQI6tG6JopIKFN4tR1Fp/fY3x4BO7pjcvz3COrrDtbmjxY9PRERkDoYlK7L0CN6W6iYTRRFZ+Xdx9n93oI5evGnWnZ+atHBUwKWZA1TNHODSzAEuysp/VhnW2eNmUSlWHkw3+5iCAHT3ccGAjh4I7+SBfn6t0NzRvJ5hvm1HRET3i2HJihpjupPGCAM7kq9h1ubkOttNf6QjIgO94aK0N4Sje597qq3muroRXZs7YFQvH8Sl5+Fizm2j7Q4KAX3atcKAjh4Y0MkdvXxda/xcvm1HRESWwLBkRXKZGy4uPQ+T1sTX2e5+nimqTzdiTmEJjqXn4ddLufj1Ui6u3zNsQgtHBUL83TCgkwfCO3qgq7czfk7NrvFtu7q6KYmIiO7FsGRFcglL5j5AfvSNx+/rLlZD7vyIoogrecX49VIujqXnIi49D7eKy43auDV3QHG5DiXl+hqPYan6iYjowcCwZEVyCUtAwx8gr6/77UbU60WkagpxLD0Xv17KQ0JGPu6WV3/bryZ8246IiMzBsGRFcgpLgDyf+Smr0GPFgUtY9svFOtt+MrE3xvRuY4WqiIhIzsz9/pbVoJRkGZGBagwJ8JbV22SO9nYI6+BuVlgquKf7joiI6H4wLD2gFHaC7LqqQvzdoFYpa33mqsr8H87hxzPX8eKgjniiq6fNzqFHRETyUPf74EQ2QmEnYP6oAAB/PGNVpernsA7ucFAIOHHlFmK+Pokh/z6ELScyUVph3vNORERE9+IzSxYgt2eW5K6uZ65uFJbgy1+vYGP8VcMI5K2dnTA13A+T+7eHqrmDVKUTEZEN4QPeVsSwZH3mvG1XVFKOzQlZWPdrhiFYtXBUYGJIOzz/sD/auDaTonQiIrIRDEtWxLBk28oq9Nh5+jo+P3wZF24UAQDs7QSM6uWDmIEdEOBjfM04lQoR0YOBYcmKGJbkQRRFHPrtJj4/fBnH0vMM6wc+5IHpj3REeEd37D2XLbthFYiIqGEYlqyIYUl+zv6uxerD6dh1VgP9/34DfFs1Q9atu9XacioVIqKmydzvb74NRw+kHm1VWPF0EA6+9hieDWsPJ3uhxqAE/DHS+cKdqdDp+bcFEdGDhmGJHmjt3Jtj4ZhArJgUZLKdCECjLUFCRr51CiMiIpvBQSmJABSbOe/cq5sSMaCTB3q1dUXvdq4IULtA6aCo12fxAXIiInlhWCIC4OmsNKtd7u0y7Ei+jh3J1wEADgoB3dQuleHJ1xW9fF3RwaNFraOGy3FePiKiBx0f8LYAPuAtfzq9iIeX7K91KhUBgKeLE94f1xNnr2mRnFWA01kFyLtTVq2ts9LeKDz19nVFa2cn7EnRYMaGxGrH5wPkRETS4NtwVsSw1DRUhRkARoGmtjAjiiJ+v3XXEJySswpw9poWpRX6asf2USmRX1yGkvLq26o+w1ulxNE3HmeXHBGRlTAsWRHDUtNxv91k5To9LmQX4fTvBUjOLMDp3wtwMec2zP0t2xQTKrsJjomI5IphyYoYlpoWSz+AXVRSjtWH0rHiQHqdbSMDvTHtYX/0ausKR3u+rEpE1JgYlqyIYYnqEpeeh0lr4s1u39xRgRB/Nwzo6IHwTu7o5u1S60Pj9+LbdkRE5jH3+5tvwxFZQYi/G9QqZa0PkAOAqpkDBnRyR/zlfOTfKcPBCzdx8MJNAIBbC0eEdXBHeCd3DOjogfbuzSEI1QMQ37YjIrI83lmyAN5ZInOY+wC5Xi/ifHYRjqXn4tdLuUjIyMedMuNxoNq4NkN4R3cM6FR558nTWcm37YiI6ondcFbEsETmasidn3KdHqezCvDrpTz8mp6LpMxbKNcZ/9p2at0C17UlKC6reXBNvm1HRFQdw5IVMSxRfdzvM0XFZRU4ceUWjl3Kxa/puTh3vZBv2xERNQCfWSKyUQo74b4CS3NHezzSuTUe6dwaAHDrThmW77+Idb9eqXPfnKKSOtsQEZExvptMJHOtWjhiSIC3WW3NndaFiIj+wLBE1ARUvW1XV2de/OU8lJg5aTAREVViWCJqAhR2AuaPCgAAk4Hpk18uYvDHh7D3XDb4uCIRkXkYloiaiMhANVZNDoK3yrirTa1SYtUzQVg+qQ/UKiV+v3UXL31zClPWJeBSzm2JqiUikg++DWcBfBuObImpt+2Kyyrw6YFLWHM4A2U6PeztBDw3wA8zn3gIzkoHiSsnIrIuDh1gRQxLJDdXcu/g7z+m4pfzOQCA1s5OeDOyK57s08bsaVWIiOSOYcmKGJZIrg6cz8GiH1ORkXsHABDUzhWLxgQisI1K4sqIiBofw5IVMSyRnJVW6LDu6BUs338RxWU6CAIwsV87/N/QLnBr4Sh1eUREjYZhyYoYlqgpyNaWYPHuNOxIvg4AcFHa428RXfBM/3awV/zxLsj9jkBORGQrzP3+lt3bcCtXroS/vz+USiWCg4Nx5MiRWttOnToVgiBUW7p3725os379+hrblJRwpGN6sHirlPhkYh/896UwdFO7oLCkAvN/OIeRy4/i+OU8AJVz2z28ZD8mrYnHrM3JmLQmHg8v2Y89KRqJqyciajyyCktbtmzB7Nmz8fbbbyMpKQkDBw7EsGHDkJmZWWP7Tz75BBqNxrBkZWXBzc0N48ePN2rn4uJi1E6j0UCp5EjH9GAK8XfDj68+jL+P6Q5VMweczy5C1Ofx+MuqY5i+IdFoEmCg8o7UjA2JDExE1GTJqhuuf//+CAoKwqpVqwzrunXrhrFjx2Lx4sV17v/9999j3LhxyMjIQPv27QFU3lmaPXs2CgoKzK6jtLQUpaWlhp8LCwvh6+vLbjhqcm7dKcO/fr6Ajcdr/oOkioDKO1NH33icXXJEJBtNrhuurKwMp06dQkREhNH6iIgIHDt2zKxjrF27FoMHDzYEpSq3b99G+/bt0bZtW4wcORJJSUkmj7N48WKoVCrD4uvrW7+TIZKJVi0c8c8ne+C9JwNNthMBaLQlSMjIt05hRERWJJuwlJubC51OBy8vL6P1Xl5eyM7OrnN/jUaD3bt344UXXjBa37VrV6xfvx4//PADNm3aBKVSiQEDBuDixYu1HmvevHnQarWGJSsrq2EnRSQTLZzszWqXU8Rn/Yio6THv/wFtiCAY3+IXRbHaupqsX78erq6uGDt2rNH60NBQhIaGGn4eMGAAgoKCsHz5cixbtqzGYzk5OcHJyan+xRPJlKezec/wmduOiEhOZHNnycPDAwqFotpdpJycnGp3m+4liiLWrVuH6OhoODqaHjfGzs4O/fr1M3lniehBE+LvBrVKaXKSXrWqchgBIqKmRjZhydHREcHBwYiNjTVaHxsbi/DwcJP7Hjp0CJcuXcK0adPq/BxRFJGcnAy1Wn1f9RI1JQo7AfNHBQBArYGpVXNHFJdVWK8oIiIrkU1YAoC5c+fiiy++wLp165CWloY5c+YgMzMT06dPB1D5LNGUKVOq7bd27Vr0798fgYHVH1JduHAh9u7di8uXLyM5ORnTpk1DcnKy4ZhEVCkyUI1Vk4PgrTLuamvV3AGOCjukagox/rM4aLR3JaqQiKhxyOqZpaioKOTl5WHRokXQaDQIDAzErl27DG+3aTSaamMuabVabN26FZ988kmNxywoKMCLL76I7OxsqFQq9OnTB4cPH0ZISEijnw+R3EQGqjEkwLvaCN7nrmvx/PqTOJ9dhLGf/op1U/uhuw/nlyOipkFW4yzZKk53QgRk5RfjufUncCnnNlo4KrDimSA81sVT6rKIiGrV5MZZIiLb5uvWHFtnhCO8ozvulOnwwlcn8Z86BrMkIpIDhiUishhVMwesfy4ETwW1hU4v4q3tZ/H+7vPQ63kDm4jki2GJiCzK0d4O/xrfE3MGdwYAfHYoHa9uTkJJuU7iyoiIGoZhiYgsThAEzBr8ED4a3wsOCgE/ndFg8hfHkX+nTOrSiIjqjWGJiBrNU8Ft8dVzIXBW2uPk1VsYt/JXZOTekbosIqJ6YVgiokYV3skD22aEo41rM1zJK8a4lb/i5BVOuEtE8sGwRESN7iEvZ2x/ORw926pwq7gcT39xHD+euS51WUREZmFYIiKr8HRWYvOLoRjczQtlFXq88p8kfHYoHRzqjYhsHcMSEVlNc0d7rI4OxtRwPwDA+7vP453vU1Ch00tbGBGRCQxLRGRVCjsBC0Z3x7sjAyAIwMbjmXjh65O4XVoBnV5EXHoediRfQ1x6HnQcn4mIbACnO7EATndC1DB7UrIxe0sSSsr1aOvaDGU6PXKKSg3b1Sol5o8KQGSgWsIqiaip4nQnRGTzIgO9sSkmFC2d7PF7wV2joAQA2doSzNiQiD0pGokqJCJiWCIiifVs64pmDooat1Xd9l64M5VdckQkGYYlIpJUQkY+bt4urXW7CECjLUFCBsdmIiJpMCwRkaRyikos2o6IyNIYlohIUp7OSou2IyKyNIYlIpJUiL8b1ColBBNtlA526OWrslpNRER/xrBERJJS2AmYPyoAAGoNTCXlejy//gQKS8qtVxgR0f8wLBGR5CID1Vg1OQjeKuOuNrVKiVlPPIQWjgrEX85H1Op45BTy2SUisi4OSmkBHJSSyDJ0ehEJGfnIKSqBp7MSIf5uUNgJSLmmxdQvTyD3dinauDbD19NC0LF1S6nLJSKZM/f7m2HJAhiWiBpfZl4xpqw7jit5xWjV3AHrpvZDn3atpC6LiGSMI3gTUZPSzr05vpsRjp5tVbhVXI6n1xzHgfM5UpdFRA8AhiUikg2Plk7YFBOKQZ1b4265Di98fRLfnsySuiwiauIYlohIVlo42eOLKX3xZJ820OlF/N93Z7Dy4CXwiQIiaiwMS0QkO472dvhofC+8NKgDAOCDPRewcGcq9Jw/jogaAcMSEcmSnZ2AecO74Z0R3QAA649dwaubk1BaoZO4MiJqahiWiEjWXhjYAcsm9YGDQsBPZzSYuo6DVxKRZTEsEZHsje7lgy+nhqCFowJxl/M4eCURWRTDEhE1CQ8/5IEtL4XBo6UT0jSFGLfqGC7fvC11WUTUBDAsEVGTEdhGhW0zwuHn3hy/37qLv3wWh6TMW1KXRUQyx7BERE3KnwevzL9TVjl45YXKwSt1ehFx6XnYkXwNcel50PHtOSIyA6c7sQBOd0Jke+6UVmDGxkQc/u0mFHYCJvdvh59Tb0Cj/eNZJrVKifmjAhAZqJawUiKSCqc7IaIH2r2DV34Vd9UoKAFAtrYEMzYkYk+KRqIqiUgOGJaIqMlytLfDB0/1RAtHRY3bq26rL9yZyi45IqoVwxIRNWknr97CnbLaB6oUAWi0JUjIyLdeUUQkKwxLRNSk5RSZN96Sue2I6MHDsERETZqns9Ki7YjowSO7sLRy5Ur4+/tDqVQiODgYR44cqbXtwYMHIQhCteX8+fNG7bZu3YqAgAA4OTkhICAA27dvb+zTICIrCfF3g1qlhFDLdgGVb8WF+LtZsywikhFZhaUtW7Zg9uzZePvtt5GUlISBAwdi2LBhyMzMNLnfhQsXoNFoDMtDDz1k2BYXF4eoqChER0fj9OnTiI6OxoQJE3D8+PHGPh0isgKFnYD5owIAoMbAJAKYPyoACrva4hQRPehkNc5S//79ERQUhFWrVhnWdevWDWPHjsXixYurtT948CAee+wx3Lp1C66urjUeMyoqCoWFhdi9e7dhXWRkJFq1aoVNmzaZVRfHWSKyfXtSNFi4M7Xa8AG+bs1w8LXHGJaIHkBNbpylsrIynDp1ChEREUbrIyIicOzYMZP79unTB2q1Gk888QQOHDhgtC0uLq7aMYcOHWrymKWlpSgsLDRaiMi2RQaqcfSNx7EpJhSfTOyNzyYHwdlJgaz8u9gQf1Xq8ojIhskmLOXm5kKn08HLy8tovZeXF7Kzs2vcR61W4/PPP8fWrVuxbds2dOnSBU888QQOHz5saJOdnV2vYwLA4sWLoVKpDIuvr+99nBkRWYvCTkBYR3eM6d0GkYFqvD6sGwDgX3sv8G04IqqVbMJSFUEwvlUuimK1dVW6dOmCmJgYBAUFISwsDCtXrsSIESPwr3/9q8HHBIB58+ZBq9UalqysrAaeDRFJ6emQdujZVoWi0gr886c0qcshIhslm7Dk4eEBhUJR7Y5PTk5OtTtDpoSGhuLixYuGn729vet9TCcnJ7i4uBgtRCQ/CjsB/xzbA4IA7Ei+jmOXcqUuiYhskGzCkqOjI4KDgxEbG2u0PjY2FuHh4WYfJykpCWr1H5NmhoWFVTvmzz//XK9jEpF89WirQnRoewDAOztSUFahl7giIrI19lIXUB9z585FdHQ0+vbti7CwMHz++efIzMzE9OnTAVR2j127dg1ff/01AGDp0qXw8/ND9+7dUVZWhg0bNmDr1q3YunWr4ZizZs3CoEGDsGTJEowZMwY7duzAvn37cPToUUnOkYis728RXbDrbDYu37yDNUcu4+XHOkldEhHZEFmFpaioKOTl5WHRokXQaDQIDAzErl270L595V+FGo3GaMylsrIyvPbaa7h27RqaNWuG7t2746effsLw4cMNbcLDw7F582a88847ePfdd9GxY0ds2bIF/fv3t/r5EZE0VM0c8M6Ibpi9JRnLfrmI0b184OvWXOqyiMhGyGqcJVvFcZaI5E8URTy95jjiLufhia6e+OLZviZf9CAi+Wty4ywRETUmQRDw97Hd4aAQ8Mv5HMSm3pC6JCKyEQxLRET/08nTGTEDOwAAFu5MRXFZhcQVEZEtYFgiIvqTVx9/CG1cm+FawV0s++WS1OUQkQ1gWCIi+pNmjgosHN0dAPDFkcu4eKNI4oqISGoMS0RE9xgc4IXB3bxQoRfxzvcp4HswRA82hiUiohosGB0ApYMdjmfkY3vSNanLISIJMSwREdWgbavmmPnEQwCA93alQVtcLnFFRCQVhiUiolq88HAHdPJsidzbZfjw5/NSl0NEEmFYIiKqhaO9Hf4+JhAAsPF4Jk5nFUhbEBFJgmGJiMiEsI7ueLJPG4gi8M73KdDp+bA30YOGYYmIqA5vDe8GZ6U9zl7TYuPxq1KXQ0RWxrBERFSH1s5OeH1oFwDAh3svIKeoROKKiMiaGJaIiMzwdP/26NlWhaKSCrz3U5rU5RCRFTEsERGZQWEn4B9jAyEIwPfJ13EsPVfqkojIShiWiIjM1LOtKyb3bw8AePf7FJRV6CWuiIisgWGJiKgeXhvaBR4tHZF+8w7WHLksdTlEZAUMS0RE9aBq5oC3R3QDACzffxFZ+cUSV0REjY1hiYionsb2boPQDm4oKddj4c5zUpdDRI2MYYmIqJ4EofJhb3s7AfvSchCbekPqkoioETEsERE1QCdPZ8QM6gAAWPDDORSXVUhcERE1FoYlIqIGmvn4Q2jj2gzXCu5i+f5LUpdDRI2EYYmIqIGaOSqwYHR3AMCaw5dxXlOIuPQ87Ei+hrj0PM4jR9RE2EtdABGRnA0J8MLgbp7Yl5aD0SuOokz3R0BSq5SYPyoAkYFqCSskovvFO0tERPfp0S6eAGAUlAAgW1uCGRsSsSdFI0VZRGQhDEtERPdBpxfx6YGan1eqik4Ld6ayS45IxhiWiIjuQ0JGPjTaklq3iwA02hIkZORbrygisiiGJSKi+5BTVHtQakg7IrI9DEtERPfB01lp0XZEZHsYloiI7kOIvxvUKiUEE22c7O3QoXULq9VERJbFsEREdB8UdgLmjwoAgFoDU2mFHiOWHcWh325arzAishiGJSKi+xQZqMaqyUHwVhl3talVSrw7MgCdvVoi93Ypnl2XgH/8mIrSCp1ElRJRQwiiKPJ91vtUWFgIlUoFrVYLFxcXqcshIono9CISMvKRU1QCT2clQvzdoLATUFKuw3u70vB13FUAQIDaBcsm9UEnz5YSV0z0YDP3+5thyQIYlojIHPtSb+D/vjuNW8XlaOagwPxRAYjq5wtBMPXEExE1FnO/v9kNR0RkJYMDvLBn9iAM6OSOu+U6vLntLP66MREFxWVSl0ZEJjAsERFZkZeLEt883x/zhnWFvZ2A3SnZGPbJERy/nCd1aURUC4YlIiIrs7MT8NIjHbHtr+Hwc28OjbYEk9bE46OfL6Bcp5e6PCK6B8MSEZFEerZ1xU8zB2J8cFvoRWD5/kuYsDoOWfnFUpdGRH/CsEREJKEWTvb4cHwvLJvUB85KeyRlFmD4J0ewI/ma1KUR0f/ILiytXLkS/v7+UCqVCA4OxpEjR2ptu23bNgwZMgStW7eGi4sLwsLCsHfvXqM269evhyAI1ZaSEs7jRETWM7qXD3bNHIi+7VuhqLQCszYnY+6WZNwurTC00elFxKXnYUfyNcSl50Gn58vMRNZgL3UB9bFlyxbMnj0bK1euxIABA7B69WoMGzYMqampaNeuXbX2hw8fxpAhQ/Dee+/B1dUVX375JUaNGoXjx4+jT58+hnYuLi64cOGC0b5KJedxIiLr8nVrjs0vhmLFgUtY9stFbEu6hlOZt/DJxD7I1t7Fwp2p0Gj/+ENOrVJi/qgARAaqJayaqOmT1ThL/fv3R1BQEFatWmVY161bN4wdOxaLFy826xjdu3dHVFQU/t//+38AKu8szZ49GwUFBWbXUVpaitLSUsPPhYWF8PX15ThLRGQxJ67kY/bmZFwruAs7AajpJlLV6EyrJgcxMBE1QJMbZ6msrAynTp1CRESE0fqIiAgcO3bMrGPo9XoUFRXBzc3NaP3t27fRvn17tG3bFiNHjkRSUpLJ4yxevBgqlcqw+Pr61u9kiIjq0M/PDbtmDcTwHt41BiUAqFq9cGcqu+SIGpFswlJubi50Oh28vLyM1nt5eSE7O9usY3z00Ue4c+cOJkyYYFjXtWtXrF+/Hj/88AM2bdoEpVKJAQMG4OLFi7UeZ968edBqtYYlKyurYSdFRGSCqpkDokPbm2wjAtBoS5CQkW+doogeQLJ6ZglAtWkBRFE0a6qATZs2YcGCBdixYwc8PT0N60NDQxEaGmr4ecCAAQgKCsLy5cuxbNmyGo/l5OQEJyenBp4BEZH5copK624EIKeIL6UQNRbZhCUPDw8oFIpqd5FycnKq3W2615YtWzBt2jR8++23GDx4sMm2dnZ26Nevn8k7S0RE1uLpbN7LJua2I6L6k003nKOjI4KDgxEbG2u0PjY2FuHh4bXut2nTJkydOhX/+c9/MGLEiDo/RxRFJCcnQ63mw5JEJL0QfzeoVUrUdv9cQOVbcSH+brW0IKL7JZuwBABz587FF198gXXr1iEtLQ1z5sxBZmYmpk+fDqDyWaIpU6YY2m/atAlTpkzBRx99hNDQUGRnZyM7OxtardbQZuHChdi7dy8uX76M5ORkTJs2DcnJyYZjEhFJSWEnYP6oAACoMTCJAOaPCoDCru7HEYioYWQVlqKiorB06VIsWrQIvXv3xuHDh7Fr1y60b1/5AKRGo0FmZqah/erVq1FRUYGXX34ZarXasMyaNcvQpqCgAC+++CK6deuGiIgIXLt2DYcPH0ZISIjVz4+IqCaRgWqsmhwEb1X1rjZvFyUGdW4tQVVEDw5ZjbNkq8wdp4GI6H7o9CISMvKRU1SCZg4KvL39LG7eLsOkEF8sHtdT6vKIZMfc72/ZPOBNRPSgU9gJCOvobvi5pZM9nll7HJsSsvBwp9YY0ZPPWhI1Bll1wxER0R/CO3ngr492BAC8ue0Mfr9VLHFFRE0TwxIRkYzNHtwZfdq5oqikcvLdCp1e6pKImhyGJSIiGXNQ2GHZxD5wdrLHqau38MkvHCOOyNIYloiIZM7XrTneG9cDALDiwCXEpedJXBFR08KwRETUBIzq5YOovr4QRWD2liTk3ymTuiSiJqNBYWnRokUoLq7+IOHdu3exaNGi+y6KiIjqb/7oAHRo3QI3Ckvx+ndnwJFhiCyjQeMsKRQKaDQaowlpASAvLw+enp7Q6XQWK1AOOM4SEdmKc9e1ePLTYyjT6bFoTHdMCfOTuiQim2Xu93eD7iyJoghBqD60/unTp+HmxvmJiIik0t1HhXnDuwIA/vFTGlKvF0pcEZH81WtQylatWkEQBAiCgM6dOxsFJp1Oh9u3b3NONSIiiU0N98PRi7n45XwOXt2UiJ2vPozmjhyDmKih6tUN99VXX0EURTz//PNYunQpVCqVYZujoyP8/PwQFhbWKIXaMnbDEZGtyb9Thsilh5FTVMrpUIhqYe73d4OeWTp06BAGDBgAe3v+pQIwLBGRbTp2KRfPrD0OUQQ+fTqI06EQ3aNRn1lydnZGWlqa4ecdO3Zg7NixeOutt1BWxtdViYhsAadDIbKMBoWll156Cb/99hsA4PLly4iKikLz5s3x7bff4vXXX7dogURE1HCcDoXo/jUoLP3222/o3bs3AODbb7/FI488gv/85z9Yv349tm7dasn6iIjoPnA6FKL71+ChA/T6yr9O9u3bh+HDhwMAfH19kZuba7nqiIjovnE6FKL706Cw1LdvX/zjH//AN998g0OHDmHEiBEAgIyMDHh5eVm0QCIiun+cDoWo4RoUlpYuXYrExES88sorePvtt9GpUycAwHfffYfw8HCLFkhERJbB6VCIGqZBQwfUpqSkBAqFAg4ODpY6pCxw6AAikgtOh0L0h0YdOqDKqVOnsGHDBmzcuBGJiYlQKpUPXFAiIpKT7j4qvMXpUIjqpUGjSubk5CAqKgqHDh2Cq6srRFGEVqvFY489hs2bN6N169aWrpOIiCzk2XA/HOF0KERma9CdpVdffRVFRUU4d+4c8vPzcevWLaSkpKCwsBAzZ860dI1ERGRBgiDgw/G94OnshPSbd/D3H1OlLonIpjXomSWVSoV9+/ahX79+RusTEhIQERGBgoICS9UnC3xmiYjk6M/ToSyf2Acezk7IKSqBp7MSIf5uUNgJdR+ESMbM/f5u0H1XvV5f47NJDg4OhvGXiIjItlVNh/LpgXTM3JyEP//lrFYpMX9UACIDOZ8cUYO64R5//HHMmjUL169fN6y7du0a5syZgyeeeMJixRERUePqpq78a/reLoZsbQlmbEjEnhSN9YsisjENCksrVqxAUVER/Pz80LFjR3Tq1An+/v4oKirC8uXLLV0jERE1Ap1exD9/SqtxW1V4WrgzFTo9x2OiB1uDuuF8fX2RmJiI2NhYnD9/HqIoIiAgAIMHD7Z0fURE1EgSMvKh0ZbUul0EoNGWICEjH2Ed3a1XGJGNqdedpf379yMgIACFhZXjcgwZMgSvvvoqZs6ciX79+qF79+44cuRIoxRKRESWlVNUe1BqSDuipqpeYWnp0qWIiYmp8YlxlUqFl156CR9//LHFiiMiosbj6ay0aDuipqpeYen06dOIjIysdXtERAROnTp130UREVHjC/F3g1qlhKkBAhwUAh7ybGm1mohsUb3C0o0bN0xOZ2Jvb4+bN2/ed1FERNT4FHYC5o8KAIBaA1O5TsT41XG4knvHeoUR2Zh6haU2bdrg7NmztW4/c+YM1GqOyUFEJBeRgWqsmhwEb5VxV5tapcT/GxmANq7NkJF7B+NWHcOpq7ckqpJIWvUawfvVV1/FwYMHceLECSiVxr9Yd+/eRUhICB577DEsW7bM4oXaMo7gTURyp9OLSMjIrzaCd05RCV746iTO/K6Fo70dlkb1xvAe/KOYmgZzv7/rFZZu3LiBoKAgKBQKvPLKK+jSpQsEQUBaWho+/fRT6HQ6JCYmwsvLyyInIRcMS0TUlBWXVWDmpmTsS7sBAHhreFfEDOwAQeB0KCRvjRKWAODq1auYMWMG9u7di6pdBUHA0KFDsXLlSvj5+d1X4XLEsERETZ1OL+LvP6Zi/bErAIDJoe2wYFR32CsaNLYxkU1otLBU5datW7h06RJEUcRDDz2EVq1aNbhYuWNYIqIHxdqjGfjHT6kQReCxLq2x4ukgtHBq0PjGRJJr9LBEf2BYIqIHyZ6UbMzekoSScj26+7hg3dR+8HLhWEwkP+Z+f8vu/unKlSvh7+8PpVKJ4ODgOkcMP3ToEIKDg6FUKtGhQwd89tln1dps3boVAQEBcHJyQkBAALZv395Y5RMRyV5koDc2xYTCvYUjzl0vxNhPf8X57EKpyyJqNLIKS1u2bMHs2bPx9ttvIykpCQMHDsSwYcOQmZlZY/uMjAwMHz4cAwcORFJSEt566y3MnDkTW7duNbSJi4tDVFQUoqOjcfr0aURHR2PChAk4fvy4tU6LiEh2+rRrhe1/HYCOrVtAoy3BX1bF4chFjrNHTZOsuuH69++PoKAgrFq1yrCuW7duGDt2LBYvXlyt/RtvvIEffvgBaWl/zKo9ffp0nD59GnFxcQCAqKgoFBYWYvfu3YY2kZGRaNWqFTZt2mRWXeyGI6IHlba4HC9+cxLHM/JhbyfgvSd7YEI/X6nLIjJLk+uGKysrw6lTpxAREWG0PiIiAseOHatxn7i4uGrthw4dipMnT6K8vNxkm9qOCQClpaUoLCw0WoiIHkSq5g74eloInuzTBhV6Ea9vPYN/7b0AGf0dTlQn2YSl3Nxc6HS6amM4eXl5ITs7u8Z9srOza2xfUVGB3Nxck21qOyYALF68GCqVyrD4+vKvKCJ6cDnZK/DxhF6Y+XgnAMCKA5cwe0sySit0EldGZBmyCUtV7h0ETRRFkwOj1dT+3vX1Pea8efOg1WoNS1ZWltn1ExE1RYIgYG5EF3zwl56wtxOwI/k6otcmoKC4DEDlOE1x6XnYkXwNcel50Ol554nkQzaDY3h4eEChUFS745OTk1PriOHe3t41tre3t4e7u7vJNqZGIXdycoKTk1NDToOIqEmb0NcXPqpmmLHhFBIy8jFu1TE8F+6PlQcvQaMtMbRTq5SYPyoAkYGcOoVsn2zuLDk6OiI4OBixsbFG62NjYxEeHl7jPmFhYdXa//zzz+jbty8cHBxMtqntmEREZNrDD3nguxnh8FEpcfnmHby7I8UoKAFAtrYEMzYkYk+KRqIqicwnm7AEAHPnzsUXX3yBdevWIS0tDXPmzEFmZiamT58OoLJ7bMqUKYb206dPx9WrVzF37lykpaVh3bp1WLt2LV577TVDm1mzZuHnn3/GkiVLcP78eSxZsgT79u3D7NmzrX16RERNRhdvZ3w3Ixz2djU/0lDVCbdwZyq75MjmyaYbDqh8zT8vLw+LFi2CRqNBYGAgdu3ahfbt2wMANBqN0ZhL/v7+2LVrF+bMmYNPP/0UPj4+WLZsGZ566ilDm/DwcGzevBnvvPMO3n33XXTs2BFbtmxB//79rX5+RERNydW8YlSYCEIiAI22BAkZ+Qjr6G69wojqSVbjLNkqjrNERFTdjuRrmLU5uc52n0zsjTG92zR+QUT3aHLjLBERkbx4Ops3X5y57YikwrBERESNIsTfDWqVErUNxCKg8q24EH83a5ZFVG8MS0RE1CgUdgLmjwoAgBoDkwhg/qgAKGp5CJzIVjAsERFRo4kMVGPV5CB4q6p3tdkJgFrVTIKqiOqHD3hbAB/wJiIyTacXkZCRj5yiEng6O+HruCvYnXID7d2b48dXH4az0kHqEukBZO73t6yGDiAiInlS2AlGwwMEqFU48/sRXM0rxv/bcQ7/juotXXFEdWA3HBERWZ2quQM+mdgbdgKwPekatiX+LnVJRLViWCIiIkn09XPD7MGdAQDvfp+CjNw7EldEVDOGJSIikszLj3VCf3833CnTYeamJJRV6KUuiagahiUiIpKMwk7A0om94drcAWevafHh3vNSl0RUDcMSERFJSq1qhg+e6gkAWHMkAwcv5EhcEZExhiUiIpJcRHdvTAmrnBT9tW9PI6eoROKKiP7AsERERDbhreHd0NXbGbm3y/C3/56GXs9hAMk2MCwREZFNUDoosHxSHygd7HDkYi6+OHpZ6pKIADAsERGRDXnIyxnzR3UHAHyw5wJOZxVIWxARGJaIiMjGTOzni+E9vFGhFzFzcxKKSsqlLokecAxLRERkUwRBwOIne6KNazPDdChEUmJYIiIim1M1HYrCTuB0KCQ5hiUiIrJJff3cMPuJhwBwOhSSFsMSERHZrL9yOhSyAQxLRERkszgdCtkChiUiIrJpnA6FpMawRERENo/ToZCUGJaIiEgWOB0KSYVhiYiIZIHToZBUGJaIiEg2OB0KSYFhiYiIZIXToZC1MSwREZGs1DQdik4vIi49DzuSryEuPQ86Ps9EFiSIosj/ou5TYWEhVCoVtFotXFxcpC6HiOiBcOpqPiasjodOL0LVzAHau3/cYVKrlJg/KgCRgWoJKyRbZ+73N+8sERGRLAW3d8OIQG8AMApKAJCtLcGMDYnYk6KRojRqYhiWiIhIlnR6EQlXbtW4rarLZOHOVHbJ0X1jWCIiIllKyMhHdmHtg1OKADTaEiRk5FuvKGqSGJaIiEiWzB3Fm6N90/1iWCIiIlnydFZatB1RbRiWiIhIlkL83aBWKSGYaKNWKRHi72a1mqhpYlgiIiJZUtgJmD8qAABqDUytnZ1QoddbryhqkhiWiIhItiID1Vg1OQjeKuOuNrfmjrC3E3Dmdy1e/PoUSsp1ElVITQEHpbQADkpJRCQtnV5EQkY+copK4Olc2fUWfzkPL3x1EnfLdQjr4I4vnu2LFk72UpdKNqTJDUp569YtREdHQ6VSQaVSITo6GgUFBbW2Ly8vxxtvvIEePXqgRYsW8PHxwZQpU3D9+nWjdo8++igEQTBaJk6c2MhnQ0RElqSwExDW0R1jerdBWEd3KOwEDOjkga+eD0FLJ3vEXc7Ds+sSUMh55KgBZBOWnn76aSQnJ2PPnj3Ys2cPkpOTER0dXWv74uJiJCYm4t1330ViYiK2bduG3377DaNHj67WNiYmBhqNxrCsXr26MU+FiIisJMTfDd9MC4GL0h4nr97C5C+Oo6C4TOqySGZk0Q2XlpaGgIAAxMfHo3///gCA+Ph4hIWF4fz58+jSpYtZxzlx4gRCQkJw9epVtGvXDkDlnaXevXtj6dKlZtdTWlqK0tJSw8+FhYXw9fVlNxwRkY1KuaZF9NrjuFVcjm5qF2yYFgL3lk5Sl0USa1LdcHFxcVCpVIagBAChoaFQqVQ4duyY2cfRarUQBAGurq5G6zdu3AgPDw90794dr732GoqKikweZ/HixYbuQJVKBV9f33qdDxERWVdgGxW2vBQGj5ZOSNMUIurzeOSYGP2b6M9kEZays7Ph6elZbb2npyeys7PNOkZJSQnefPNNPP3000bp8ZlnnsGmTZtw8OBBvPvuu9i6dSvGjRtn8ljz5s2DVqs1LFlZWfU7ISIisrrOXs7470uh8HZR4lLObUxYHYfrBXelLotkQNKwtGDBgmoPV9+7nDx5EgAgCNVH0RBFscb19yovL8fEiROh1+uxcuVKo20xMTEYPHgwAgMDMXHiRHz33XfYt28fEhMTaz2ek5MTXFxcjBYiIrJ9HVq3xH9fCkPbVs1wJa8YE1bHISu/WOqyyMZJ+g7lK6+8UuebZ35+fjhz5gxu3LhRbdvNmzfh5eVlcv/y8nJMmDABGRkZ2L9/f53BJigoCA4ODrh48SKCgoLqPgkiIpKVdu7N8d+XwvD0mnhcySvG+M/i8J+Y/ujQuqXUpZGNkjQseXh4wMPDo852YWFh0Gq1SEhIQEhICADg+PHj0Gq1CA8Pr3W/qqB08eJFHDhwAO7u7nV+1rlz51BeXg61Wm3+iRARkaz4uDarDExfHP9fl1w8/hPTH529nKUujWyQLJ5Z6tatGyIjIxETE4P4+HjEx8cjJiYGI0eONHoTrmvXrti+fTsAoKKiAn/5y19w8uRJbNy4ETqdDtnZ2cjOzkZZWeVro+np6Vi0aBFOnjyJK1euYNeuXRg/fjz69OmDAQMGSHKuRERkHZ4uSmx+MRRdvZ2Re7sUEz+Px7nrWqnLIhski7AEVL6x1qNHD0RERCAiIgI9e/bEN998Y9TmwoUL0Gor/0P//fff8cMPP+D3339H7969oVarDUvVG3SOjo745ZdfMHToUHTp0gUzZ85EREQE9u3bB4VCYfVzJCIi6/Jo6YTNL4aiZ1sV8u+UYdLn8UjOKpC6LLIxshhnydZxuhMiInkrLCnHc1+ewKmrt9DSyR5fPtcP/fzcpC6LGlmTGmeJiIioMbkoHfD18yEI6+CO26UVmLI2Ab9eypW6LLIRDEtEREQAWvzvjtKgzq1xt1yH59afwIELOdDpRcSl52FH8jXEpedBp2eHzIOG3XAWwG44IqKmo7RCh5c3JmFf2g0o7CrvOt0q/mMCXrVKifmjAhAZyLem5Y7dcERERA3gZK/AqslBCGrnCp0eRkEJALK1JZixIRF7UjQSVUjWxrBERER0DztBwPWCmueOq+qOWbgzlV1yDwiGJSIionskZOQj28REuyIAjbYECRn51iuKJMOwREREdI+cotqDUkPakbwxLBEREd3D01lp0XYkbwxLRERE9wjxd4NapYRgoo1apUSIPweufBAwLBEREd1DYSdg/qgAAKg1MP310U5Q2JmKU9RUMCwRERHVIDJQjVWTg+CtMu5qc1RUBqQfTl/j23APCHupCyAiIrJVkYFqDAnwRkJGPnKKSuDprIRapcSIZUdw4sotrD16GS8O6ih1mdTIeGeJiIjIBIWdgLCO7hjTuw3COrrDz6MF3h1Z2UX3r72/4UJ2kcQVUmNjWCIiIqqnqH6+eLyrJ8p0eszZkoyyCr3UJVEjYlgiIiKqJ0EQ8P5TPdCquQNSNYVY9stFqUuiRsSwRERE1ACezkr888keAICVBy8hMfOWxBVRY2FYIiIiaqDhPdQY09sHehH4239P426ZTuqSqBEwLBEREd2HRaMD4eXihIzcO3h/d5rU5VAjYFgiIiK6D6rmDvjwL70AAF/FXcXRi7kSV0SWxrBERER0nwZ1bo3o0PYAgP/77jS0d8slrogsiWGJiIjIAuYN7wo/9+bQaEuw8IdzUpdDFsSwREREZAHNHe3x0YResBOAbUnXsCdFI3VJZCEMS0RERBYS3N4N0x+pnP7kre0puFlUKnFFZAkMS0RERBY0a/BD6OrtjPw7ZZi37QxEkZPtyh3DEhERkQU52Svw76jecFAI2JeWg29P/S51SXSfGJaIiIgsrJvaBXOHdAEALNqZiqz8YokrovvBsERERNQIXhzUAcHtW+F2aQX+77vT0OvZHSdXDEtERESNQGEn4KPxvdDMQYH4y/n48tgVqUuiBmJYIiIiaiR+Hi3w9ohuAIAle87j4o0iiSuihmBYIiIiakTP9G+HQZ1bo6xCj7n/PY1ynV7qkqieGJaIiIgakSAI+OCpnlA1c8DZa1p8euCS1CVRPTEsERERNTJvlRKLxnQHACzffwlnfi+QtiCqF4YlIiIiKxjdywcjeqqh04uYsyUZJeU6qUsiMzEsERERWYEgCPjHmEC0dnZC+s07+HDvBalLIjMxLBEREVlJqxaO+OCpngCAtUczcCw9V+KKyBwMS0RERFb0WFdPTArxBQD837dnUFRSLnFFVBeGJSIiIit7e0QAfN2a4VrBXSz84Rzi0vOwI/ka4tLzoONI3zZHNmHp1q1biI6OhkqlgkqlQnR0NAoKCkzuM3XqVAiCYLSEhoYatSktLcWrr74KDw8PtGjRAqNHj8bvv3PSQyIiajwtnezx0fjeAIDvEq9h0pp4zNqcjElr4vHwkv3Yk6KRtkAyIpuw9PTTTyM5ORl79uzBnj17kJycjOjo6Dr3i4yMhEajMSy7du0y2j579mxs374dmzdvxtGjR3H79m2MHDkSOh3fUiAiosaTf6e0xvXZ2hLM2JDIwGRD7KUuwBxpaWnYs2cP4uPj0b9/fwDAmjVrEBYWhgsXLqBLly617uvk5ARvb+8at2m1WqxduxbffPMNBg8eDADYsGEDfH19sW/fPgwdOtTyJ0NERA88nV7Ewp2pNW4TAQgAFu5MxZAAbyjsBKvWRtXJ4s5SXFwcVCqVISgBQGhoKFQqFY4dO2Zy34MHD8LT0xOdO3dGTEwMcnJyDNtOnTqF8vJyREREGNb5+PggMDDQ5HFLS0tRWFhotBAREZkrISMfGm1JrdtFABptCRIy8q1XFNVKFmEpOzsbnp6e1dZ7enoiOzu71v2GDRuGjRs3Yv/+/fjoo49w4sQJPP744ygtLTUc19HREa1atTLaz8vLy+RxFy9ebHh2SqVSwdfXt4FnRkRED6KcotqDUkPaUeOSNCwtWLCg2gPY9y4nT54EUDmY171EUaxxfZWoqCiMGDECgYGBGDVqFHbv3o3ffvsNP/30k8m66jruvHnzoNVqDUtWVpaZZ0xERAR4Oist2o4al6TPLL3yyiuYOHGiyTZ+fn44c+YMbty4UW3bzZs34eXlZfbnqdVqtG/fHhcvXgQAeHt7o6ysDLdu3TK6u5STk4Pw8PBaj+Pk5AQnJyezP5eIiOjPQvzdoFYpka0tQW0DBSgd7NCnnas1y6JaSHpnycPDA127djW5KJVKhIWFQavVIiEhwbDv8ePHodVqTYaae+Xl5SErKwtqtRoAEBwcDAcHB8TGxhraaDQapKSk1Ou4RERE9aGwEzB/VACAyoe5a1JSrserm5JQWsG3s6Umi2eWunXrhsjISMTExCA+Ph7x8fGIiYnByJEjjd6E69q1K7Zv3w4AuH37Nl577TXExcXhypUrOHjwIEaNGgUPDw88+eSTAACVSoVp06bhb3/7G3755RckJSVh8uTJ6NGjh+HtOCIiosYQGajGqslB8FYZd7WpVUq8/FhHONrbITb1Bl78+hQn3ZWYLIYOAICNGzdi5syZhjfXRo8ejRUrVhi1uXDhArRaLQBAoVDg7Nmz+Prrr1FQUAC1Wo3HHnsMW7ZsgbOzs2Gff//737C3t8eECRNw9+5dPPHEE1i/fj0UCoX1To6IiB5IkYFqDAnwRkJGPnKKSuDprESIvxsUdgLCOnjgha9P4NBvN/H8+hP44tm+aO4om6/tJkUQRZHjqt+nwsJCqFQqaLVauLi4SF0OERE1EQkZ+XjuywTcKdOhb/tW+PK5fnBWOkhdVpNh7ve3LLrhiIiIHkQh/m745oX+cFba4+TVW5i8NgHaYk68a20MS0RERDYsqF0rbIoJRavmDjidVYBJa+KRd7vmqVKocTAsERER2bjANipsfjEMHi0dkaopxMTP4zlgpRUxLBEREclAF29nbH4xDF4uTriYcxsTV8dDo70rdVkPBIYlIiIimejk2RL/fSkMbVyb4XLuHUxYHYes/GKpy2ryGJaIiIhkpL17C/x3ehjauzdHVv5dTFgdh4zcO1KX1aQxLBEREclMG9dm+O9LYejYugU02hJMWB2HizeKpC6ryWJYIiIikiEvFyW2vBSGrt7OuFlUiqjP45F6vVDqspokhiUiIiKZ8mjphM0vhqJHGxXy75Rh0pp4nM4qkLqsJodhiYiISMZcmztiY0x/BLVzhfZuOSZ/cRwnr+RLXVaTwrBEREQkcy5KB3w9rT/6+7uhqLQCU9Yl4Fh6rtRlNRkMS0RERE1ASyd7rH8uBAMf8kBxmQ7PfXkCBy/kQKcXEZeehx3J1xCXngednlPC1hcn0rUATqRLRES2oqRch1f+k4h9aTmwtxPgrLTHrT/NJ6dWKTF/VAAiA9USVmkbOJEuERHRA0jpoMDKZ4LRx9cVFXrRKCgBQLa2BDM2JGJPikaiCuWHYYmIiKiJUdgJ0Ghrnjuuqjtp4c5UdsmZiWGJiIioiUnIyEd2Ye0T7YoANNoSJGTwrTlzMCwRERE1MTlFtQelP7uUw1G/zcGwRERE1MR4OivNardwZyre+O4MLuXcbuSK5I1hiYiIqIkJ8XeDWqWEYKKNg0JAhV7ElpNZGPzxIbzw1QmcuJIPviRfHcMSERFRE6OwEzB/VAAAVAtMwv+W5ZP6YOuMMEQEeEEQgH1pORj/WRzGrTqGPSkaPvz9JxxnyQI4zhIREdmiPSkaLNyZavRmXE3jLKXfvI0vjlzG1sRrKKvQAwD8PVrghYH+eCqoLZQOCqvXbg3mfn8zLFkAwxIREdkqnV5EQkY+copK4OmsRIi/GxR2NXfQ3SwqxVfHruCb+KvQ3q0cn8m9hSOeDfdDdGh7tGrheF/HtzUMS1bEsERERE3JndIKbDmRhbVHM3Ct4C4AoJmDAlH9fDHtYX/4ujUHYP6dK1vFsGRFDEtERNQUlev02HVWg9WHLiNVUwgAsBOAET190MPHBYt3n8e9IaLqntKqyUE2H5gYlqyIYYmIiJoyURTx66U8rD6cjiMXc+tsLwDwVilx9I3HbbpLjnPDERERkUUIgoCHH/LAN9P646eZD2NAR3eT7ZvaCOEMS0RERGS27j4qTOjna1Zbc0cSt3UMS0RERFQv5o4Qbm47W8ewRERERPVizgjhSgc7+Hu0sFpNjYlhiYiIiOrF1AjhVUrK9Yj85DC2J/0u+ylUGJaIiIio3iID1Vg1OQjeKuOuNrVKibeGdUWA2gUFxeWYs+U0nl9/AhrtXYkqvX8cOsACOHQAERE9qGobwbtcp8fnhy/jk30XUabTw9nJHvOGd8OkEF8Igm0MJ8BxlqyIYYmIiKhmF28U4fWtZ5CUWQAACOvgjiVP9UQ79+bSFgaOs0REREQ24CEvZ3w3PRzvjOgGpYMd4i7nYejSw1h3NAM6vTzu1zAsERERUaNS2Al4YWAH7J09CKEd3HC3XIdFP6Zi/GfHcCnnttTl1YlhiYiIiKyivXsL/OeFUPzzyUC0dLJHYmYBhi87gk8PXEKFTi91ebViWCIiIiKrsbMT8Ez/9vh5ziA82qU1yir0+HDvBYxd+StSrxdKXV6NGJaIiIjI6nxcm+HLqf3w0fheUDVzQMq1QoxecRQf/3wBpRU6AJVv2sWl52FH8jXEpedJ9oyTbMLSrVu3EB0dDZVKBZVKhejoaBQUFJjcRxCEGpcPP/zQ0ObRRx+ttn3ixImNfDZEREQkCAKeCm6L2LmDMLS7Fyr0Ipbtv4RRy49i1cF0PLxkPyaticeszcmYtCYeDy/Zjz0pGuvXKZehA4YNG4bff/8dn3/+OQDgxRdfhJ+fH3bu3FnrPtnZ2UY/7969G9OmTcOlS5fQoUMHAJVhqXPnzli0aJGhXbNmzaBSqcyujUMHEBER3R9RFLHrbDb+344U5N0pq7FN1ehMqyYHITJQfd+fae73t/19f5IVpKWlYc+ePYiPj0f//v0BAGvWrEFYWBguXLiALl261Lift7e30c87duzAY489ZghKVZo3b16tLREREVmPIAgY0VONEH83DPxgP0rKqz/wLaIyMC3cmYohAd5Q2FlncEtZdMPFxcVBpVIZghIAhIaGQqVS4dixY2Yd48aNG/jpp58wbdq0ats2btwIDw8PdO/eHa+99hqKiopMHqu0tBSFhYVGCxEREd2/Szm3awxKVUQAGm0JEjLyrVaTLO4sZWdnw9PTs9p6T0/Pal1ttfnqq6/g7OyMcePGGa1/5pln4O/vD29vb6SkpGDevHk4ffo0YmNjaz3W4sWLsXDhwvqdBBEREdUpp6jEou0sQdI7SwsWLKj1Ieyq5eTJkwBQ4zwyoiiaPb/MunXr8Mwzz0CpNJ7wLyYmBoMHD0ZgYCAmTpyI7777Dvv27UNiYmKtx5o3bx60Wq1hycrKqsdZExERUW08nZV1N6pHO0uQ9M7SK6+8UuebZ35+fjhz5gxu3LhRbdvNmzfh5eVV5+ccOXIEFy5cwJYtW+psGxQUBAcHB1y8eBFBQUE1tnFycoKTk1OdxyIiIqL6CfF3g1qlRLa2BDW9gSYA8FZVTthrLZKGJQ8PD3h4eNTZLiwsDFqtFgkJCQgJCQEAHD9+HFqtFuHh4XXuv3btWgQHB6NXr151tj137hzKy8uhVt//U/ZERERUPwo7AfNHBWDGhkQIgFFgqupLmj8qwGoPdwMyecC7W7duiIyMRExMDOLj4xEfH4+YmBiMHDnS6E24rl27Yvv27Ub7FhYW4ttvv8ULL7xQ7bjp6elYtGgRTp48iStXrmDXrl0YP348+vTpgwEDBjT6eREREVF1kYFqrJocBG+VcVebt0ppsWED6kMWD3gDlW+szZw5ExEREQCA0aNHY8WKFUZtLly4AK1Wa7Ru8+bNEEURkyZNqnZMR0dH/PLLL/jkk09w+/Zt+Pr6YsSIEZg/fz4UCkXjnQwRERGZFBmoxpAAbyRk5COnqASezpVdb9a8o1RFNoNS2jIOSklERCQ/5n5/y6IbjoiIiEgqDEtEREREJjAsEREREZnAsERERERkAsMSERERkQkMS0REREQmMCwRERERmcCwRERERGQCwxIRERGRCbKZ7sSWVQ2CXlhYKHElREREZK6q7+26JjNhWLKAoqIiAICvr6/ElRAREVF9FRUVQaVS1bqdc8NZgF6vx/Xr1+Hs7AxBsNwEf4WFhfD19UVWVtYDMefcg3S+PNem60E6X55r0/WgnK8oiigqKoKPjw/s7Gp/Mol3lizAzs4Obdu2bbTju7i4NOn/WO/1IJ0vz7XpepDOl+fadD0I52vqjlIVPuBNREREZALDEhEREZEJDEs2zMnJCfPnz4eTk5PUpVjFg3S+PNem60E6X55r0/WgnW9d+IA3ERERkQm8s0RERERkAsMSERERkQkMS0REREQmMCwRERERmcCwJLGVK1fC398fSqUSwcHBOHLkiMn2hw4dQnBwMJRKJTp06IDPPvvMSpXen8WLF6Nfv35wdnaGp6cnxo4diwsXLpjc5+DBgxAEodpy/vx5K1XdMAsWLKhWs7e3t8l95HpdAcDPz6/G6/Tyyy/X2F5O1/Xw4cMYNWoUfHx8IAgCvv/+e6PtoihiwYIF8PHxQbNmzfDoo4/i3LlzdR5369atCAgIgJOTEwICArB9+/ZGOgPzmTrX8vJyvPHGG+jRowdatGgBHx8fTJkyBdevXzd5zPXr19d4rUtKShr5bEyr67pOnTq1Ws2hoaF1HtcWrytQ9/nWdI0EQcCHH35Y6zFt9do2FoYlCW3ZsgWzZ8/G22+/jaSkJAwcOBDDhg1DZmZmje0zMjIwfPhwDBw4EElJSXjrrbcwc+ZMbN261cqV19+hQ4fw8ssvIz4+HrGxsaioqEBERATu3LlT574XLlyARqMxLA899JAVKr4/3bt3N6r57NmztbaV83UFgBMnThida2xsLABg/PjxJveTw3W9c+cOevXqhRUrVtS4/YMPPsDHH3+MFStW4MSJE/D29saQIUMM80XWJC4uDlFRUYiOjsbp06cRHR2NCRMm4Pjx4411GmYxda7FxcVITEzEu+++i8TERGzbtg2//fYbRo8eXedxXVxcjK6zRqOBUqlsjFMwW13XFQAiIyONat61a5fJY9rqdQXqPt97r8+6desgCAKeeuopk8e1xWvbaESSTEhIiDh9+nSjdV27dhXffPPNGtu//vrrYteuXY3WvfTSS2JoaGij1dhYcnJyRADioUOHam1z4MABEYB469Yt6xVmAfPnzxd79epldvumdF1FURRnzZolduzYUdTr9TVul+t1BSBu377d8LNerxe9vb3F999/37CupKREVKlU4meffVbrcSZMmCBGRkYarRs6dKg4ceJEi9fcUPeea00SEhJEAOLVq1drbfPll1+KKpXKssVZWE3n+uyzz4pjxoyp13HkcF1F0bxrO2bMGPHxxx832UYO19aSeGdJImVlZTh16hQiIiKM1kdERODYsWM17hMXF1et/dChQ3Hy5EmUl5c3Wq2NQavVAgDc3NzqbNunTx+o1Wo88cQTOHDgQGOXZhEXL16Ej48P/P39MXHiRFy+fLnWtk3pupaVlWHDhg14/vnn65xUWo7X9c8yMjKQnZ1tdO2cnJzwyCOP1Po7DNR+vU3tY4u0Wi0EQYCrq6vJdrdv30b79u3Rtm1bjBw5EklJSdYp8D4dPHgQnp6e6Ny5M2JiYpCTk2OyfVO5rjdu3MBPP/2EadOm1dlWrte2IRiWJJKbmwudTgcvLy+j9V5eXsjOzq5xn+zs7BrbV1RUIDc3t9FqtTRRFDF37lw8/PDDCAwMrLWdWq3G559/jq1bt2Lbtm3o0qULnnjiCRw+fNiK1dZf//798fXXX2Pv3r1Ys2YNsrOzER4ejry8vBrbN5XrCgDff/89CgoKMHXq1FrbyPW63qvq97Q+v8NV+9V3H1tTUlKCN998E08//bTJSVa7du2K9evX44cffsCmTZugVCoxYMAAXLx40YrV1t+wYcOwceNG7N+/Hx999BFOnDiBxx9/HKWlpbXu0xSuKwB89dVXcHZ2xrhx40y2k+u1bSh7qQt40N3717coiib/Iq+pfU3rbdkrr7yCM2fO4OjRoybbdenSBV26dDH8HBYWhqysLPzrX//CoEGDGrvMBhs2bJjhn3v06IGwsDB07NgRX331FebOnVvjPk3hugLA2rVrMWzYMPj4+NTaRq7XtTb1/R1u6D62ory8HBMnToRer8fKlStNtg0NDTV6MHrAgAEICgrC8uXLsWzZssYutcGioqIM/xwYGIi+ffuiffv2+Omnn0yGCDlf1yrr1q3DM888U+ezR3K9tg3FO0sS8fDwgEKhqPZXR05OTrW/Tqp4e3vX2N7e3h7u7u6NVqslvfrqq/jhhx9w4MABtG3btt77h4aGyu4vlxYtWqBHjx611t0UrisAXL16Ffv27cMLL7xQ733leF2r3nCsz+9w1X713cdWlJeXY8KECcjIyEBsbKzJu0o1sbOzQ79+/WR3rdVqNdq3b2+ybjlf1ypHjhzBhQsXGvQ7LNdray6GJYk4OjoiODjY8OZQldjYWISHh9e4T1hYWLX2P//8M/r27QsHB4dGq9USRFHEK6+8gm3btmH//v3w9/dv0HGSkpKgVqstXF3jKi0tRVpaWq11y/m6/tmXX34JT09PjBgxot77yvG6+vv7w9vb2+jalZWV4dChQ7X+DgO1X29T+9iCqqB08eJF7Nu3r0FBXhRFJCcny+5a5+XlISsry2Tdcr2uf7Z27VoEBwejV69e9d5XrtfWbFI9WU6iuHnzZtHBwUFcu3atmJqaKs6ePVts0aKFeOXKFVEURfHNN98Uo6OjDe0vX74sNm/eXJwzZ46Ympoqrl27VnRwcBC/++47qU7BbDNmzBBVKpV48OBBUaPRGJbi4mJDm3vP99///re4fft28bfffhNTUlLEN998UwQgbt26VYpTMNvf/vY38eDBg+Lly5fF+Ph4ceTIkaKzs3OTvK5VdDqd2K5dO/GNN96otk3O17WoqEhMSkoSk5KSRADixx9/LCYlJRneAHv//fdFlUolbtu2TTx79qw4adIkUa1Wi4WFhYZjREdHG73h+uuvv4oKhUJ8//33xbS0NPH9998X7e3txfj4eKuf35+ZOtfy8nJx9OjRYtu2bcXk5GSj3+HS0lLDMe491wULFoh79uwR09PTxaSkJPG5554T7e3txePHj0txigamzrWoqEj829/+Jh47dkzMyMgQDxw4IIaFhYlt2rSR5XUVxbr/OxZFUdRqtWLz5s3FVatW1XgMuVzbxsKwJLFPP/1UbN++vejo6CgGBQUZvUr/7LPPio888ohR+4MHD4p9+vQRHR0dRT8/v1r/w7Y1AGpcvvzyS0Obe893yZIlYseOHUWlUim2atVKfPjhh8WffvrJ+sXXU1RUlKhWq0UHBwfRx8dHHDdunHju3DnD9qZ0Xavs3btXBCBeuHCh2jY5X9eqYQ7uXZ599llRFCuHD5g/f77o7e0tOjk5iYMGDRLPnj1rdIxHHnnE0L7Kt99+K3bp0kV0cHAQu3btahNB0dS5ZmRk1Po7fODAAcMx7j3X2bNni+3atRMdHR3F1q1bixEREeKxY8esf3L3MHWuxcXFYkREhNi6dWvRwcFBbNeunfjss8+KmZmZRseQy3UVxbr/OxZFUVy9erXYrFkzsaCgoMZjyOXaNhZBFP/3JCkRERERVcNnloiIiIhMYFgiIiIiMoFhiYiIiMgEhiUiIiIiExiWiIiIiExgWCIiIiIygWGJiIiIyASGJSIiIiITGJaIiBrAz88PS5culboMIrIChiUisnlTp07F2LFjAQCPPvooZs+ebbXPXr9+PVxdXautP3HiBF588UWr1UFE0rGXugAiIimUlZXB0dGxwfu3bt3agtUQkS3jnSUiko2pU6fi0KFD+OSTTyAIAgRBwJUrVwAAqampGD58OFq2bAkvLy9ER0cjNzfXsO+jjz6KV155BXPnzoWHhweGDBkCAPj444/Ro0cPtGjRAr6+vvjrX/+K27dvAwAOHjyI5557Dlqt1vB5CxYsAFC9Gy4zMxNjxoxBy5Yt4eLiggkTJuDGjRuG7QsWLEDv3r3xzTffwM/PDyqVChMnTkRRUZGhzXfffYcePXqgWbNmcHd3x+DBg3Hnzp1G+rdJROZiWCIi2fjkk08QFhaGmJgYaDQaaDQa+Pr6QqPR4JFHHkHv3r1x8uRJ7NmzBzdu3MCECROM9v/qq69gb2+PX3/9FatXrwYA2NnZYdmyZUhJScFXX32F/fv34/XXXwcAhIeHY+nSpXBxcTF83muvvVatLlEUMXbsWOTn5+PQoUOIjY1Feno6oqKijNqlp6fj+++/x48//ogff/wRhw4dwvvvvw8A0Gg0mDRpEp5//nmkpaXh4MGDGDduHDjXOZH02A1HRLKhUqng6OiI5s2bw9vb27B+1apVCAoKwnvvvWdYt27dOvj6+uK3335D586dAQCdOnXCBx98YHTMPz//5O/vj7///e+YMWMGVq5cCUdHR6hUKgiCYPR599q3bx/OnDmDjIwM+Pr6AgC++eYbdO/eHSdOnEC/fv0AAHq9HuvXr4ezszMAIDo6Gr/88gv++c9/QqPRoKKiAuPGjUP79u0BAD169LiPf1tEZCm8s0REsnfq1CkcOHAALVu2NCxdu3YFUHk3p0rfvn2r7XvgwAEMGTIEbdq0gbOzM6ZMmYK8vLx6dX+lpaXB19fXEJQAICAgAK6urkhLSzOs8/PzMwQlAFCr1cjJyQEA9OrVC0888QR69OiB8ePHY82aNbh165b5/xKIqNEwLBGR7On1eowaNQrJyclGy8WLFzFo0CBDuxYtWhjtd/XqVQwfPhyBgYHYunUrTp06hU8//RQAUF5ebvbni6IIQRDqXO/g4GC0XRAE6PV6AIBCoUBsbCx2796NgIAALF++HF26dEFGRobZdRBR42BYIiJZcXR0hE6nM1oXFBSEc+fOwc/PD506dTJa7g1If3by5ElUVFTgo48+QmhoKDp37ozr16/X+Xn3CggIQGZmJrKysgzrUlNTodVq0a1bN7PPTRAEDBgwAAsXLkRSUhIcHR2xfft2s/cnosbBsEREsuLn54fjx4/jypUryM3NhV6vx8svv4z8/HxMmjQJCQkJuHz5Mn7++Wc8//zzJoNOx44dUVFRgeXLl+Py5cv45ptv8Nlnn1X7vNu3b+OXX35Bbm4uiouLqx1n8ODB6NmzJ5555hkkJiYiISEBU6ZMwSOPPFJj119Njh8/jvfeew8nT55EZmYmtm3bhps3b9YrbBFR42BYIiJZee2116BQKBAQEIDWrVsjMzMTPj4++PXXX6HT6TB06FAEBgZi1qxZUKlUsLOr/f/mevfujY8//hhLlixBYGAgNm7ciMWLFxu1CQ8Px/Tp0xEVFYXWrVtXe0AcqLwj9P3336NVq1YYNGgQBg8ejA4dOmDLli1mn5eLiwsOHz6M4cOHo3PnznjnnXfw0UcfYdiwYeb/yyGiRiGIfC+ViIiIqFa8s0RERERkAsMSERERkQkMS0REREQmMCwRERERmcCwRERERGQCwxIRERGRCQxLRERERCYwLBERERGZwLBEREREZALDEhEREZEJDEtEREREJvx/tCdlhXSJ7kEAAAAASUVORK5CYII=\n"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(results[\"costs\"], \"o-\")\n",
    "plt.xlabel(\"Iterations\")\n",
    "plt.ylabel(\"Cost\")\n",
    "\n",
    "print(\"Minimized circuit output:\", circuit(params))\n",
    "print(\"Optimized parameters:\", params)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-02-19T21:50:26.756331Z",
     "start_time": "2024-02-19T21:50:26.338055Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Note: Charges shown are estimates based on your Amazon Braket simulator and quantum processing unit (QPU) task usage. Estimated charges shown may differ from your actual charges. Estimated charges do not factor in any discounts or credits, and you may experience additional charges based on your use of other services such as Amazon Elastic Compute Cloud (Amazon EC2).\n",
      "Estimated cost to run this example: 0.386 USD\n"
     ]
    }
   ],
   "source": [
    "job_cost = job.result()[\"braket_tasks_cost\"]\n",
    "sv1_cost = float(braket_tasks_cost.simulator_tasks_cost())\n",
    "\n",
    "print(\n",
    "    \"Note: Charges shown are estimates based on your Amazon Braket simulator and quantum processing unit (QPU) task usage. Estimated charges shown may differ from your actual charges. Estimated charges do not factor in any discounts or credits, and you may experience additional charges based on your use of other services such as Amazon Elastic Compute Cloud (Amazon EC2).\"\n",
    ")\n",
    "print(f\"Estimated cost to run this example: {job_cost + sv1_cost :.3f} USD\")"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "name": "python3",
   "language": "python",
   "display_name": "Python 3 (ipykernel)"
  },
  "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.10.13"
  },
  "vscode": {
   "interpreter": {
    "hash": "590fab68195cf107911461461f81d5c472d3d6127f579badfcfad30f03e5cab2"
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}