{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "67b87b66",
   "metadata": {},
   "source": [
    "# Quick Start"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "ceb957a2d72e49ed",
   "metadata": {},
   "source": [
    "Let us create a pulse that implements a CZ gate on two atoms in the perfect Rydberg blockade regime."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "993cd5239bf336bf",
   "metadata": {
    "ExecuteTime": {
     "end_time": "2025-11-29T17:16:03.935633Z",
     "start_time": "2025-11-29T17:16:02.997193Z"
    }
   },
   "outputs": [],
   "source": [
    "# %pip install -q --progress-bar off rydopt # Uncomment for installation on Colab\n",
    "\n",
    "import rydopt as ro\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "d9820d79460ae60",
   "metadata": {},
   "source": [
    "First, we create the target gate, by specifying the gate phases, the interaction between the atoms and the Rydberg-state decay rate."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "304aeddde0908214",
   "metadata": {
    "ExecuteTime": {
     "end_time": "2025-11-29T17:16:03.941927Z",
     "start_time": "2025-11-29T17:16:03.940236Z"
    }
   },
   "outputs": [],
   "source": [
    "gate = ro.gates.TwoQubitGate(phi=None, theta=np.pi, Vnn=float(\"inf\"), decay=0.0)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "3614e8b18eb99467",
   "metadata": {},
   "source": [
    "Then we choose a pulse ansatz that consists of a constant detuning and a sweep of the laser phase. The Rabi frequency amplitude is constant and fixed."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "51803009cfb742a8",
   "metadata": {
    "ExecuteTime": {
     "end_time": "2025-11-29T17:16:03.948572Z",
     "start_time": "2025-11-29T17:16:03.947016Z"
    }
   },
   "outputs": [],
   "source": [
    "pulse_ansatz = ro.pulses.PulseAnsatz(\n",
    "    detuning_ansatz=ro.pulses.Const(), phase_ansatz=ro.pulses.SinCrab(2)\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "3e5480112b7a3352",
   "metadata": {},
   "source": [
    "A guess of pulse parameters can serve as input for the optimization. Here, we choose to optimize four pulse parameters and package them as `ro.pulses.PulseParams(...)`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "9cb296616a8b8649",
   "metadata": {
    "ExecuteTime": {
     "end_time": "2025-11-29T17:16:03.952967Z",
     "start_time": "2025-11-29T17:16:03.951502Z"
    }
   },
   "outputs": [],
   "source": [
    "initial_params = ro.pulses.PulseParams(7.0, [0.0], [0.0, 0.0], [])"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f3e2471e0359cf73",
   "metadata": {},
   "source": [
    "Now, we perform the optimization, specifying a target gate infidelity (which also serves as tolerance for the differential equation solver)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "f68ae25336698ec5",
   "metadata": {
    "ExecuteTime": {
     "end_time": "2025-11-29T17:16:07.496526Z",
     "start_time": "2025-11-29T17:16:03.955934Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Started optimization using 1 process\n",
      "\n"
     ]
    },
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "1af3d2bed57941d0b173330f0052802e",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "proc00:   0%|          | 0/1000 [00:00<?, ?it/s]"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "=== Optimization finished using Adam ===\n",
      "\n",
      "Runtime: 8.395 seconds\n",
      "Gates with infidelity below tol=1.0e-10: 1\n",
      "\n",
      "Optimized gate:\n",
      "> infidelity <= tol\n",
      "> parameters = (7.626643597908443, -0.031259443394905904, 0.8197775118421416, 0.6638904306444476)\n",
      "> duration = 7.626643597908443\n"
     ]
    }
   ],
   "source": [
    "opt_result = ro.optimization.optimize(gate, pulse_ansatz, initial_params, tol=1e-10)\n",
    "optimized_params = opt_result.params"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "b7158c159a033f80",
   "metadata": {},
   "source": [
    "Finally, we can plot the resulting pulse."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "1dcbb0e378a481e4",
   "metadata": {
    "ExecuteTime": {
     "end_time": "2025-11-29T17:16:08.001093Z",
     "start_time": "2025-11-29T17:16:07.540615Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnAAAAHQCAYAAAAh51fQAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAYmwAAGJsBSXWDlAAAkBxJREFUeJzs3XdcleX7wPHPYcieCjhw78zcOJA0sTT9ac40K0eWlpk5srS0LDP9WmauTMtylJl7pebeAwduBWcuHICggCJwzu+Px8PhxFTO4TmHc71fL14v7mdep8fg4n7u+7o1Op1OhxBCCCGEsBp2agcghBBCCCGejCRwQgghhBBWRhI4IYQQQggrIwmcEEIIIYSVkQROCCGEEMLKSAInhBBCCGFlJIETQgghhLAyksAJIYQQQlgZSeCEEEIIIayMJHBCCCGEEFZGEjghhBBCCCsjCZwQQgghhJWRBE4IIYQQwspIAieEEEIIYWUc1A5AWBY3NzdSUlLw9/dXOxQhhBCiULp9+zaOjo4kJiY+9TUkgRNGUlJSSEtLUzsMUUC0Wi0AdnbSGW8r5JnbHnnmlscUv2clgRNG/P390Wq1XLt2Te1QRAG4ffs2gPS42hB55rZHnrnlCQwMzPc1JB0XQgghhLAyksAJIYQQQlgZSeCEEEIIIayMJHBCCCGEEFZGEjghhBBCCCsjCZwQQgghhJWRBE4IIYQQwspIApeNiIgIpk2bRu/evalZsyYODg5oNBq+/vrrfF138+bNtGnThmLFiuHi4kK1atX47LPPSEhIyPG88+fP07t3bwIDA3FyciIwMJDevXtz8eLFfMUjhBBCCOsjCVw2Zs6cyaBBg5g3bx4nT540SdXkyZMn8+KLL7JhwwZq1KhBu3btiI+P55tvvqF+/fpER0dned6ePXuoVasW8+bNw9vbm44dO+Lt7c28efN47rnn2L9/f75jE0IIIYT1kAQuG88++ywfffQRf/zxB2fOnOHNN9/M1/XCw8MZNmwY9vb2/P333+zYsYPFixdz4cIFQkNDiYiI4N133810XlJSEq+++ipJSUmMHDmSkydPsmjRIk6ePMnIkSNJTEzk1Vdf5cGDB/mKTwghhBDWQ5bSysbbb79t1M7vGnLjx49Hp9PRp08fXn755fTtrq6uzJkzhwoVKrBs2TLOnj1LtWrV0vfPnTuXGzduUKVKlUyvb7/++muWLVtGZGQk8+fPp3///vmKUQghhBDWQXrgCsCjR4/4+++/AejRo0em/WXLliU4OBiAFStWGO3Tt7t3754pibSzs6Nbt24ALF++3ORxCyGEEMIySQJXACIjI0lKSgKgfv36WR6j3x4eHm60Xd9+0vOEEEIIUXjJK9QCcOnSJQC8vb3x8PDI8pjSpUsbHQtw//59YmJiAChTpkyO5925c4fExETc3NzyHW+qNpXzd8/n+zrC8sXejwXgnuM9lSMRBUWeue2RZ255UrWpONjlLwWTBK4A3L9/HyDH5Mrd3R2Ae/cM/4Ppz8vpXP15+nPzksAFBgZmuy8qKgo7Lzs6ru6Y63WEEEII8eRiHsYQ4BqQr2vIK1QhhBBCCCsjPXAFQP/aNDExMdtj9IV8PT09M52X07kZCwBnPDcn165dy3ZfYGAgt5Ju5ek6QgghhFCHJHAFoFy5cgDExcVx//79LMfBXb161ehYUBI4X19fYmNjuXLlCrVq1cr2vGLFiplk/BuATxEfVrRfkfuBwurFxipjY3x9fVWORBQUeea2R5655Wk6smm+ryEJXAGoWrUqrq6uJCUlcejQIV544YVMxxw6dAiAunXrGm2vW7cumzdv5tChQ7Rr1y7P5+WHg50DlXwqmex6wnLdTrkNgL+Pv8qRiIIiz9z2yDO3PPmdwAAyBq5AFClShLZt2wKwcOHCTPv//fdf9u7dC0DHjsaTB/TtRYsWodVqjfZptVr++usvADp16mTyuIUQQghhmSSBM6Hp06dTrVo1evbsmWnfiBEj0Gg0/Pbbb2zYsCF9e1JSEn379iUtLY3OnTsbrcIA0Lt3b0qWLElkZCSjR4822jd69GgiIyMJDAzM8p5CCCGEKJzkFWo2jhw5woABA9LbFy5cAGDWrFmsXbs2ffuKFSsoUaIEANHR0URERFC8ePFM16tbty6TJk1i6NChtGnThmbNmuHv78+uXbuIioqiatWq/PTTT5nOc3V1ZfHixbz00kt88803rF69mmeffZaTJ09y8uRJ3NzcWLJkCS4uLqb+TyCEEEIICyUJXDbu3bvHgQMHMm2/du2a0SzO5OTkPF9zyJAh1KxZk0mTJhEWFkZiYiJlypRh5MiRjBw5Mtsiv8HBwRw7doyxY8eyefNmli1bhp+fHz179uTzzz+nYsWKT/4BhRBCCGG1NDqdTqd2EMJyBAYGotVquXHjhtqhiAJw+/bjwc3+MrjZVsgztz3yzC2PvqB+TmW9ciNj4IQQQgghrIwkcEIIIYQQVkYSOCGEEEIIKyMJnBBCCCGElZEETgghhBDCykgCJ4QQQghhZSSBE0IIIYSwMpLACSGEEEJYGUnghBBCCCGsjCRwQgghhBBWRhI4IYQQQggrIwmcEEIIIYSVkQROCCGEEMLKSAInhBBCCGFlJIETQgghhLAyksAJIYQQQlgZSeCEEEIIIayMJHBCCCGEEFZGEjghhBBCCCsjCZwQQgghhJWRBE4IIYQQwspIAieEEEIIYWUkgRNCCCGEsDKSwAkhhBBCWBlJ4IQQQgghrIwkcEIIIYQQVkYSOCGEEEIIKyMJnBBCCCGElZEETgghhBDCykgCJ4QQQghhZSSBE0IIIYSwMpLACSGEEEJYGUngcrFkyRKaN2+Oj48Pbm5u1KpVi4kTJ5KSkvJE1ylXrhwajSbXr6+++srovO3bt+d6zk8//WTKjyyEEEIIC+egdgCWbPDgwUyZMgUHBwdatGiBu7s7W7du5ZNPPmHNmjVs3LgRFxeXPF2rS5cuREdHZ7kvNjaWNWvWAPDCCy9keUxAQACtW7fOcl/VqlXzFIMQQgghCgdJ4LKxcuVKpkyZgru7Ozt27KBu3boAREdH06JFC3bv3s3o0aP57rvv8nS9nI6bOHEia9asoUqVKoSEhGR5TLVq1Zg7d+4Tfw4hhBBCFD7yCjUb33zzDQAjRoxIT94AihUrxo8//gjA9OnTiY+Pz/e9fv31VwDeeuutfF9LCCGEEIWfJHBZuH79OgcPHgSgR48emfY3bdqU0qVLk5yczLp16/J1rz179hAREYGDgwO9evXK17WEEEIIYRvkFWoWwsPDAfD19aV8+fJZHlO/fn2uXr1KeHg4r7322lPfS9/71qZNG4oXL57tcbdu3eKrr77i+vXrODs7U61aNdq2bUuZMmWe+t5CCCGEsE6SwGXh0qVLADkmR6VLlzY69mkkJiayePFiAPr27ZvjsWfPnuWLL74w2ubg4MAHH3zAxIkTcXCQRymEEELYCvmtn4X79+8D4Obmlu0x7u7uANy7d++p77N48WISEhIoXrw4bdq0yfIYLy8vBg8eTMeOHalSpQqenp5cuHCB3377jenTpzN58mQSEhKYPXt2nu8bGBiY7b6oqCj8/Py4ffv2E38eYX3i4uLUDkEUMHnmtkeeueXRarXY2eVvFJuMgVPRnDlzAOjZs2e2PWh16tRh8uTJPP/88xQvXhxXV1dq1qzJ999/z6JFiwD4+eefOXr0aEGFLYQQQgiVSQ9cFjw8PADlFWd2EhISAPD09Hyqe0RGRrJnzx7g6WefdurUidq1a3P06FHWrFlD7dq183TetWvXst0XGBiIVqvF39//qWIS1kmet+2RZ2575Jlbjvz2voH0wGWpXLlyAFy9ejXbY/T79Mc+Kf3khaZNm+arEG/16tWBnJMyIYQQQhQuksBloU6dOgDExMRkO0nh0KFDAEY14vIqLS2N+fPnA7lPXshNTEwMYOg1FEIIIUThJwlcFgIDA2nQoAEACxcuzLR/9+7dXL16FScnp2wnH+Rk3bp1REVF4eHhQdeuXZ86zuvXr7Nr1y4AgoKCnvo6QgghhLAuksBl49NPPwVgwoQJHDlyJH17TEwMAwYMAGDgwIF4eXml71uxYgXVqlUjNDQ0x2vrX5927949x5muAFOmTMlyDdXjx4/Trl07Hjx4QMWKFXnllVfy9sGEEEIIYfVkEkM2OnTowKBBg5g6dSqNGjUiNDQUNzc3tmzZQlxcHMHBwYwdO9bonPj4eCIiInj48GG21719+zZ///03kLfXp1988QXDhg2jdu3alC9fHjs7Oy5cuEB4eDharZYyZcqwZs0anJyc8veBhRBCCGE1JIHLwZQpUwgODmbGjBns3buXlJQUKlasyIgRIxgyZAhFihR54msuWLCAlJQUatSoQcOGDXM9/rPPPmPPnj2cOnWKTZs2kZiYiKenJ02aNOGVV16hf//+Mv5NCCGEsDEanU6nUzsIYTn0ZURu3LihdiiiAOgLNkt5Adshz9z2yDO3PPqC+vmpICFj4IQQQgghrIwkcEIIIYQQVkYSOCGEEEIIKyMJnBBCCCGElZEETgghhBDCykgCJ4QQQghhZQqsDlx4eDhr167lzJkzxMbGAuDr60v16tVp27btU60pKoQQQphF6iNIuAX3b0JKIuh0YOcAbn7g7g8uPqDRqB2lsGFmT+CSkpLo1asXy5Ytw93dnSpVquDj4wNAZGQka9euZcyYMXTo0IEFCxbg6upq7pCEEEIIY/HX4NwmuLIPoo5BdCTotNkf71oMiteEUnWhYiiUDgJ7x4KLV9g8sydww4cPZ9euXSxbtoz27dtjb29vtF+r1bJ69Wree+89hg8fzowZM8wdkhBCCAGJMXD8Lzi2EG6eeLJzk6Lh4jbla9ckcPKCGq9A7TeUZE5654SZmX0lBj8/P6ZOncprr72W43ELFy7kww8/5M6dO+YMR+RCVmKwLVKh3fbIMweiz8PuyXBiMaQ9Mv31/WtA08FQoxPYq79ipTxzy2OKlRjM/i/r4cOHFC1aNNfjfH19c1wEXgghhMiX+Ouw5Us4sSSb16MaKFlH6UEr/hz4VgCP4uDkqfSopSZD4h24dx1unYKoo3BpJzyMN77M7VOw/B3YNg5ajoFnOkiPnDA5sydwISEhjBkzhnr16mWbyMXExDB27FhCQkLMHY4QQghbk/IQ9k6D3d9DSlLm/WWDoXYPqNwK3P1yvpZnCSjxHFR9WWmnpcL1Q0pSeGKJcTJ39zIs6Q2lG0Hb75Qxc0KYiNlfoV64cIEXXniBmJgYWrRoQfXq1fH29gYgLi6OM2fOsG3bNnx9fdm6dSuVKlUyZzgiF/IK1bbIqxXbY3PP/PphWN4fYs4Zb7dzhNqvQZNBUKyyae6V8lAZU7dnCsRe+M/9HOD5jyFkaIFPdrC5Z24FTPEK1ewJHMC9e/f46aefWL9+PadPn+bu3bsA+Pj4UL16ddq0aUP//v3x8vIydygiF5LA2Rb5wW57bOaZp6Uqkwt2/A90acb7nusOoaPBK9A899amwfHFyuva+1HG+4rXhM6/gl8V89w7CzbzzK2I1SRwwnpIAmdb5Ae77bGJZ55wR3l1+e9u4+0lakObb5UxbgXhUSLs/E7pkcuYRBZxh1dmQI0OBRKGTTxzK2OKBE5WYhBCCFF43AiH2c2Nkzc7B3hhFLy9peCSN4AibtDyC3h7M/hVN2x/lABLesE/nym9dUI8BVUTuNjYWK5cuaJmCEIIIQqLUyvh19ZwL0OvRtFKSgLVbLh6JT1K1YX+OyCov/H2fdPhrzfgURYTK4TIhckTuOTkZGbMmMFLL71EyZIlcXFxoWTJkjRu3Jjx48cbvZobNmwYFSpUMHUIQgghbM3BX5TXpqkZylFVaQ3vbFVKg6jNwQnaTITOc8Axw4pDEetgXjtIjFYvNmGVTJrA7du3j0qVKjFo0CA2b97MzZs3SU5O5ubNmxw4cIBRo0ZRvXp1Zs6cmX6ODMETQgjx1HQ62DYe/h4GZPh98vxw6P4nOFvY5LiaXZRXuV6lDduuH4I5Lyl16oTII5P1J+/du5fQ0FCSk5OpWLEiPXr0oHbt2nh6enL37l3Cw8NZsmQJ58+fZ+DAgZw/f95UtxZCCGGLdDrYPAb2/GDYprGH9lOhzhtqRZW7gGeg7yb4oyvceryEV+wFmNsGeq0F79I5ny8EJpqF+uDBAypVqkRUVBRfffUVn332GZosqk7rdDp+/fVXBg0axMOHD3FzcyMxMZG0NBnEaSlkFqptkdlptqfQPHOdDrZ+Dbu+M2xzcIaucw1Fdi3dw3vKGLhLOwzbvMtC77XgXcZktyk0z7wQsZhZqLNnzyYqKorhw4czatSoLJM3AI1GQ9++fdmzZw8BAQEkJCSY4vZCCCFszfYJxsmboxu8sdx6kjcAZ0/o8RdUDDVsi/tXGRN3/5Z6cQmrYJIEbtWqVbi7u/PFF1/k6fjatWuzYMGCbBM9IYQQIlsHZsOOCYa2oyu8vgTKBasX09NydIHuC6HyS4Ztdy/D753hQZxaUQkrYJIE7tSpUzRu3BhXV9fcD34sNDSUdevW8euvv5oiBGFCmuT7aocghBBZO70a1n9saDu4QI/F1pm86Tk6Q7ffoVJLw7ZbJ2BRD0h5oF5cwqKZJIGLi4vDzy+XBYCz0KpVK3r16mWKEIQJaVIS4fA8tcMQQghj/+6DZW+TPttUYw/dFkD5EFXDMgkHJ3h1PgQ2MGz7dw8s7wdarXpxCYtlkgTO29ubO3fuPPF5GzduZP78+aYIQZja30Ph0i61oxBCCEX0OfizO6QlG7a1nwqVX1QvJlMr4qb0JvpVM2w7sxq2j1cvJmGxTJLA1ahRg3379pGUlPdq0tu2baNNmza89dZbpghBmJo2FRa/CbEX1Y5ECGHrHsbDn6/BwzjDthc+s+xSIU/L1VeZjOFe3LBt50Q4sVS9mIRFMkkC165dOxISEhgzZkyejj958iRvvvkmWq1WCvlasgd3YWF35YenEEKoQauF5f0h5pxhW92eSqHewsqrFLy2UCmLorfqfbh2WL2YhMUxSQLXv39/AgICmDRpEmPHjs0xKVuwYAGNGzfmxo0buLm5meL2wsR0ji6GRnSE8sNTxmAIIdSwYwJErje0yzSBNpOgsFcxKFUPOvxoaKc+hL9eh4QnH64kCieTFPIF2LVrFy+++CIpKSlUqFCB119/PX0lhvj4eMLDw1m6dCkRERHodDqGDBlCbGws8+fPl0K+FkQp5JvGja9rw5W9hh0tv4Smg9UKS5iJFPi0PVb1zM+sVZIWPc9S0G87uFtB7Kay7RvY8T9Du3wzeHMF2Nnn+RJW9cxthCkK+ZpsKa2QkBC2bNlCt27duHDhAmPHjs10jE6nw8PDg4kTJ9K/f3/69OljqtsLk9IoM7t+CoH7j1dk2PKVMjvKmqfqCyGsx93LsHKAoW3vpJTasKXkDaDZCLh5EiL+VtqXdigJ3QufqhuXUJ1JF7MPDg7m/PnzTJ06lZYtW+Lv74+DgwPFihWjQYMGjB07ljNnztC/f39ASehkDJyFcisGXX9TpukD6NJg6VuQcFvduIQQhV9aCiztC8kZxt+2mwKl6qoXk1rs7JRXqd5lDdt2TITzW9SLSVgEk7xCXbp0KW3atHmiQr4AMTExJCQkULZs2dwPFgUi01qoe6fBxlGGAyo0hzdWKD9UhNWTVyu2xyqe+cbRsHeqoV3nDXhlhnrxWIIb4TDnJUh7pLRdi8K7e8CzRK6nWsUztzEWsxbqq6++ip+fH507d+aPP/7g3r17eTqvaNGiFp+8LVmyhObNm+Pj44Obmxu1atVi4sSJpKSkPNF15s6di0ajyfFrw4YN2Z5/69YtBg4cSPny5XFyciIgIICuXbty5MiR/H7EnDUeCNX+z9C+uB32/5jt4UIIkS/nNhknb37V4OWJ6sVjKUrWgdYZlg9LioGV78kEMxtmkjFwo0ePZvny5axYsYKVK1fi6OhIaGgonTt3pn379hQrVswUtylwgwcPZsqUKTg4ONCiRQvc3d3ZunUrn3zyCWvWrGHjxo24uLjkfqEMKlasSNOmTbPcV6pUqSy3R0ZGEhISwu3bt6lQoQIdOnTg0qVLLF26lJUrV7J48WI6duz4xJ8vTzQa5S/fqOMQf0XZtuVLpSeu+LPmuacQwjYlRsOKdw1tBxfoOlcpcCug/lvK6gwnlynti9sgbDY0ejfn80ShZLJZqKAkGkuXLmX58uUcOXIEjUaDnZ0dzz//PJ07d6Zjx46UKJF7d68lWLlyJR07dsTd3Z0dO3ZQt64y9iI6OpoWLVpw4sQJhg0bxnfffZen682dO5c+ffrQq1cv5s6dm+c4dDod9erVIzw8nDfffJPffvsNe3tlXNrs2bPp378/7u7unDt3juLFi+dytdxleoWq9+9e+K0N6UvY+D8D72xT1vATVkterdgei33mOp1SPPzMGsO2dlOhniy3aOTBXZjZFO49fvXm4KzMzPWvnu0pFvvMbZjFvELVq1KlCp9++imHDh3i0qVLfPvttwQFBbF9+3YGDhxI6dKlCQ4O5vvvv+fy5cumvLXJffPNNwCMGDEiPXkDKFasGD/+qLxCnD59OvHx5i1yu379esLDw/H29ubHH39MT94A+vXrR2hoKAkJCUyZMsWscVC2CTQdYmjfPq3MTBVCCFM4scQ4eav2f0rBXmHMxQc6zgQe18FLfQjL34HU5BxPE4WP2Uaily1blqFDh7Jnzx6uX7/O9OnTadasGWFhYXz00UdUrFiR+vXr880333D27FlzhfFUrl+/zsGDBwHo0aNHpv1NmzaldOnSJCcns27dOrPGsmLFCgDat2+Pu7t7pv36+JYvX27WOABoPhJK1DK098+Ay7vNf18hROF27was+8jQdi0G//dD4S/W+7TKPw9NPjC0b54wrhUnbEKBTCUsXrw4AwYMYMuWLdy6dYtffvmF1q1bc/LkSUaNGkWNGjXy/CqyIISHhwPg6+tL+fLlszymfv36Rsfm1fnz5xk1ahT9+vVj6NCh/Prrr0RHR+cai/5+2cVx7tw5EhMTnyiWJ+ZQBDr9ooxL0Vv9ATzK+xq4QghhRKeDVQONl+z7v8ng7qdeTNagxSgIqGlo7/4Boo6pFo4oeAVeC8LX15e33nqLv//+m9u3b7NgwQI6duyIxoL+0rp06RIAZcqUyfaY0qVLGx2bV3v27GHcuHH8/PPPTJ48mb59+1K6dGn+97+s/3rKLRZ9HDqdrmBeS/tVgdDRhnbsRdj+jfnvK4QonI7MgwsZapo91w2eaa9ePNbCwQk6/gR2j+ci6tKU9VLTnqxCgrBeJluJ4WmkpqYSEhLC66+/nvvBBej+/fsAOa7Vqn+dmdeSKcWLF+ezzz6jffv2VKhQAScnJyIiIpg2bRoLFixgxIgRpKWl8emnxtW1c4sl42vVvMaiHzyZlaioKPz8/NIHvWapfCd8AhbjeOsoALp9M7hbojmpAc/l6f7CcsTFxakdgihglvTM7RJv4fvPZ+k9CWluAcTWH44up58/wsDOH7e6/XE79LhG3s0TJGz8hqR67xkdZknPXCi0Wi12+aynavIeuOTkZGbMmMFLL71EyZIlcXFxoWTJkjRu3Jjx48cbzW4cNmwYFSpUMHUIFql169Z8/fXXBAUFUaxYMTw8PKhfvz7z5s1Lf3381VdfcevWLZUjzQM7e+698A06O0cANDotnttGGgpMCiFEHrjv/hq7Rwnp7fvNx6Fz9lIxIuuTWO89Un0qpbfdDk7DPva8ihGJgmLSHrh9+/bx6quvcuPGDaMlsm7evMnNmzcJCwtjwoQJTJgwgffeU/5CsMSltDw8PAByHFOWkKD80PH09Mz3/T788EPGjx9PdHQ0Gzdu5M033zSKJTY2NttY9HE8SSw5TVvWlxHJdbq5vz80+wS2fQ2AQ2wk/mcXwAsj8xSDsCxSXsD2qP7MI9bDhQzFy2t2xbtBV/XisWadZ8GcF0GnRaNNoejuL+CtDZkWvFf9mYt0+e19AxMmcHv37iU0NJTk5GQqVqxIjx49qF27Np6enty9e5fw8HCWLFnC+fPnGThwIOfPW+5fCOXKlQPg6tWr2R6j36c/Nj/s7e2pXLky0dHRmZKrcuXKERsby5UrV3KMQ6PRFPyqFk0Hw+lVcOuE0t79PdTsAsUqF2wcQgjrknwf/s4w69TZG1qNVy0cqxdYHxoNgH3Tlfa1MGVsYf231I1LmJVJXqE+ePCArl27kpyczFdffUVkZCRffvklHTt2JDQ0lC5dujBu3DgiIiL4+eefcXZ25ocffmDZsmWmuL3J1alTB1DWas1uksKhQ4cAjGrE5UdMTAxg6P3T019ff7/s4qhcuXKWZUbMyt4RXpkOmsf/jNIeKaUALLBXVQhhQbaOMxSiBXjpa5l1ml8vfAY+5QztzV8qK1uIQsskCdzs2bOJiopi+PDhjBo1KtsZpRqNhr59+7Jnzx4CAgKMXv9ZksDAQBo0aADAwoULM+3fvXs3V69excnJiTZt2uT7fkeOHCEyMhKAoKAgo336JbJWr16d5WtUfXydOnXKdxxPpWRtCOpnaF/cDqdWqBOLEMLy3QiHsFmGdrkQZbF6kT9FXKFNhnJcD+Ng0xeqhSPMzyQJ3KpVq3B3d+eLL/L2j6V27dosWLDAokqH/Jd+NuiECROMFoyPiYlhwIABAAwcOBAvL8OA2xUrVlCtWjVCQ0ONrpWUlMSMGTPSZ5RmtHPnTjp37gwoBYL/m8C9/PLL1KlTh7i4OAYMGEBaWlr6vtmzZ7Nlyxbc3d358MMP8/mJ8+GFT8E9wND+51PlFYkQQmSk1cK64aB7vAC7vZMU7DWlyi9C9XaG9tHf4cp+9eIRZmWStVADAgKoXbs2//zzzxOd988//3Dz5k169bLMte4+/PBDpk6diqOjI6Ghobi5ubFlyxbi4uIIDg5m06ZNRovZ69c7LVu2rFFNtri4OHx8fHBycqJOnTqUKVOG1NRUIiMjOXnyJAA1a9bkn3/+yXKt2IiICEJCQrhz5w4VKlSgQYMGXLp0ibCwMBwcHEy6mH22a6Hm5sRSWNbX0G48EFqNM0lMwnxkjUTbo+ozP7oQVmYocdFshEx8MrX4azA9CFIev7Hxr8HtTkvBzkH+P7cgFrMWalxcHH5+Tz5+oVWrVhabvAFMmTKFv/76i8aNG7N3717WrVtHYGAgEyZMYOvWrUbJW05cXV0ZPXo0LVq04NatW6xfv57Vq1dz69YtWrZsyaxZszh06FCWyRtA1apVOX78OO+//z5paWmsWLGCS5cu0alTJw4cOGCy5C1fnu2svArR2z8Tbp5ULx4hhGV5GG/8Ss+rjDIRSpiWVyA0+9jQvn0KlxML1ItHmI2qPXAbN27k5s2b9OwpCxZbiqfugQO4EwEzg0H7uBJ42abQe628HrFg0gNne1R75v98ZpglCfDqAllxwVxSH8GsELijrDOudXQj5vXN+JV7RuXAhJ7F9MDVqFGDffv2kZSU9zUxt23bRps2bXjrLZnmXGj4VYUmAw3tf3fD2bXqxSOEsAx3IuDAT4Z2+WbGY7WEaTkUgbaT0pt2KYm4h/2gXjzCLEySwLVr146EhATGjBmTp+NPnjzJm2++iVartchCviIfQoYZT2jYOApSk9WLRwihLp0O1n8C2lSlbecAL0+UnnlzK9cUnnklvel8ZgncPKFiQMLUTJLA9e/fn4CAACZNmsTYsWNzTMoWLFhA48aNuXHjRo5rjQor5eQBoZ8b2ncvG//lLYSwLWf/hovbDO2g/uBfTb14bMmLX4F9EUBZ8pANI6VOZyFikjFwALt27eLFF18kJSWFChUq8Prrr6evxBAfH094eDhLly4lIiICnU7HkCFDiI2NZf78+UalMYS68jUGTk+rhdnN4OZxpe3kCR8ckUKdFkjGwNmeAn3mqY/gx4YQe1Fpu/nBB4dB1jstOJvHwO7Jhnb3hVCtrWrhCIUpxsCZbCmtkJAQtmzZQrdu3bhw4QJjx47NdIxOp8PDw4OJEyfSv39/+vTpY6rbC0tiZwetx8Pcxz8kku/BtnHQ7gdVwxJCFLDDvxmSN4AWoyV5K2hNh6I9PB+7B8pqP2wcBZVeVMbJCatmkleoesHBwZw/f56pU6fSsmVL/P39cXBwoFixYjRo0ICxY8dy5swZ+vfvDygJnYyBK6TKNYXqGWaYHZknZUWEsCUP42H7BEPb/xlZcUENzp4kNBxqaMdehLDZ6sUjTMZkr1CfRkxMDAkJCQW/CLvIlkleoerFXoIZQcoaqQAVmkPPVfm/rjAZeYVqewrsmf/31d3ry6ByS/PeU2Tp9s0ofJZ0wDFGKSuCkxd8eBRcfVWNy5ZZTBmRp1W0aFFJ3goz3/LQKEPV9Yvb4cK2bA8XQhQScVeVYt56FZpDpdBsDxdmZmdPQtPPDO3keNg1KfvjhVVQNYETNiBkGLj4GNqbxyiTHIQQhdfWryH14eOGBl4cK2VDVJZSqhFUaW3YEDZbSbSF1ZIETpiXsxeEfGRoRx2F0ytUC0cIYWZRx+D4X4Z2re5Q4jn14hEGoZ8DjxPptEewfbyq4Yj8kQROmF+Dt8Ez0NDe+jWkpagXjxDCfDZ9ATweWu3gDC1GqRqOyCCghpJQ6x37E26fUS8ekS+SwAnzc3SGFhnGX8ReVGalCiEKl0u7jIv2NnpPWVxdWI4XPk0v7otOC1u+Ujce8dQkgRMF47luShkBve3/g+QE9eIRQpiWTgdbM9T/dPaC4MGqhSOy4V0GGrxjaEesg3/3qRePeGqSwImCYWdvvMRW4m3jWWpCCOt2bhNcPWBoNxkELt6qhSNyEDJMWSFHb/MYWWLLCkkCJwpOldZQprGhvXcqJMWqF48QwjS0WuPeNzc/aPiuevGInLkVVRJsvav74dxG9eIRT0USOFFwNBpo+aWhnXwP9s1QLx4hhGmcWW1Y+xge9/C4qxePyF3jAeCWoZjztnHSC2dlVEvgKlSoQLVq1fj2229JTk5WKwxR0Mo0VNbh0zvwEyTGqBePECJ/tGnKL389z1JQT9a5tnhF3JREWy/qmDIeTlgN1RK4y5cvExkZyYgRI6hcuTK//fabWqGIgvbCSMP3jxKUV6lCCOt0/C+IjjS0m32izDwXlq9eb/AoYWhvGy+F1q2Iagncb7/9xpw5c3jvvfdwcnLi7bffVisUUdBK1YMqLxvaYbMh4bZ68Qghnk7qf4rB+laA2j3Ui0c8GUdn4164Wyfg7Br14hFPRLUErlevXvTp04fp06dz7tw5Ll68qFYoQg0Ze+FSkmDPFPViEUI8naO/Q9wVQ7v5p2DvqF484snV7WlcaF164ayGxUxikEXtbUyJWlDt/wztg7/A/ZvqxSOEeDKpj2DX94a2/zPwbGf14hFPx8EJns+w3OGdM7LcoZWwmARO2KDmGXrhUh/C7snqxSKEeDLHF0F8hsXQm30CdvIrxSrVfl0p8Ku3fYIyOUVYNPm/Tain+LPwTAdD+9BvEH9dtXCEEHmUlgI7vzO0/apB9fbqxSPyx6EIPP+xoR0dCSeXqRePyBMHU1ykRYsWprgMGo2GLVu2mORawko0HwGnVwE6SEtWeuHafpfraUIIFZ1YAnH/GtrPD5feN2tXqzvsmgR3Lynt7ROUV+J29urGJbJlkgRu+/btprgMGo3GJNcRVsS/uvJD4uRSpX1kvjIew6O4unEJIbKWlgo7vzW0i1aGGh3Vi0eYhr2j8hp85eMVNGIvwKkVULOLunGJbJkkgQNo3bo1n3zyyVOfP2HCBDZulKU8bNLzww0JXFoy7JsOL32tbkxCiKydXAaxGaoGPD9cemkKi5pdYccEuHtZae+aBDU6Se+qhTJZAle8eHGaNWv21OfPnTvXVKEIa+P/ePzMmdVK++CvEDxEWa9PCGE5tGnGvW++FWXmaWFi7wBNh8Kax+uk3j6trM5Q/f9yPk+owiRpdZUqVShRokTuB+agePHiVKlSxRThCGuUcRp7SiIcmKleLEKIrJ1aATHnDO3nP1J+6YvCo9ZrynJoeru+kzVSLZRJErizZ88ybty43A/Mwfjx4zlz5owpwhHWqEQtqPySoX1gNjyMVy8eIYQxrda4982nnPLKTRQuDkUg+END+0Y4XJDJhZZIXmwLy/H8cMP3yfEQ9rN6sQghjJ1ZBXfOGtohw2TVhcKqbk9w8zO0d0plAEskCZywHKWDoFyIob1vBjxKVC8eIYRCpzNedcGrDDzXXb14hHk5ukDjgYb2lX1weY968YgsmWTwwpUrV3I/KAdlypTJ/SBhG54fDpd3Kd8/iFWK+zYZmPM5QgjzurAVbh43tIMHKa/aROHVoK9Sl/NhnNLe+S2UC1Y1JGHMJAlcuXLlnrqGm0ajITU11RRhiMKg/PMQGATXwpT23qnQ4G1wdFY3LiFsWcZl7tz8oM4b6sUiCoaTBzR6D7aPV9oXt8G1QxBYX924RDqTvEItU6ZMll86nS79y9PTE09Pz/S2/rzSpUubIgSzWbJkCc2bN8fHxwc3Nzdq1arFxIkTSUlJeaLrhIeHM378eEJDQwkICMDR0REfHx9CQkKYMWNGttfbvn07Go0mx6+ffvrJFB/VMmg0xjNSE25B+AL14hHC1l07ZOgVB2g0QHnFJgq/oH5QxMPQlrFwFsUkPXCXL182amu1Wrp160ZiYiKfffYZvXr1wtvbG4D4+HjmzZvHuHHjqF+/Pn/99ZcpQjCLwYMHM2XKFBwcHGjRogXu7u5s3bqVTz75hDVr1rBx40ZcXHL/QZaamkrdunUBcHd3p0GDBgQEBHDt2jX27dvH7t27mT9/Pv/880/6f6f/CggIoHXr1lnuq1q16lN/RotU+SUoXhNunlDae6dCvT5SrkAINWTsfXPyVF6tCdvg6gtBbxv+DUSuh1unIKCGunEJwISFfDOaPHkya9as4dChQzz77LNG+7y8vBg0aBAtWrSgXr16TJo0ieHDh2dzJfWsXLmSKVOm4O7uzo4dO9ITsOjoaFq0aMHu3bsZPXo0332Xt79I6tWrxyeffEL79u1xcnJK337ixAlatWpFWFgYQ4cO5ddff83y/GrVqtlOsWONRpnhtqS30o67AqdXypIuQhS0OxFwdq2hXf8tcPZSLx5R8Bq9D/tnQupDpb13GnQsRG99rJhZZqHOnTuXZs2aZUreMnr22Wdp3rw58+bNM0cI+fbNN98AMGLEiPTkDaBYsWL8+OOPAEyfPp34+NxrlTk4OHDo0CG6du1qlLwB1KxZk4kTJwKwaNGiJ341W2hVbw++FQztPT9IMUkhCtruHwzf2zspr0+FbXH/z5jHE0sg7qp68Yh0ZkngLly4QLFixXI9rmjRoly8eDHX4wra9evXOXjwIAA9evTItL9p06aULl2a5ORk1q1bl+/71alTB4AHDx4QHR2d7+sVCnb20OQDQ/vmCWUmnBCiYMRdhROLDe06r4NHgHrxCPU0Hgiax+mCNlXpkROqM0sC5+bmRlhYWPpkhazodDoOHjyIm5ubOULIl/DwcAB8fX0pX758lsfUr1/f6Nj8OHdOWZqmSJEi+Pr6ZnnMrVu3+Oqrr+jfvz8ffvghM2fOzHf5FotX6zXjYpJ7pqgXixC2Zt8M5Zc1KL+8M/5BJWyLb3l4poOhfXguJMWqFY14zCwJXPPmzbl48SLDhw8nLS0t0/60tDQ+/vhjLly4QPPmzc0RQr5cunQJyLk+nX72rP7Yp6XT6dJfof7f//1fplesemfPnuWLL75g9uzZTJ06lQEDBlCxYkWGDh1aeMuwOLpAw3cN7Us7lGVdhBDmlRgDRzIMb6nR0XhIg7A9wYMM36ckwqE56sUiADNNYvjqq6/YsGEDkydPZunSpbz66qvpPVmXL19m8eLFXLlyBTc3N7788ktzhJAv9+/fB8ixd9Dd3R2Ae/fu5eteX375Jfv27cPd3Z0JEyZk2u/l5cXgwYPp2LEjVapUwdPTkwsXLvDbb78xffp0Jk+eTEJCArNnz87zPQMDA7PdFxUVhZ+fH7dv336qz2NqmnKvUNRhEnapSQA83DKRe62kJ85U4uLi1A5BFLC8PHO3sCm4pSSlt2Of6UWqhfxMEE/OJP+fO5TCO7AJRa7tBUC7bybRlbqBg9TofBparRY7u/z1oZklgatevTrr16+nR48eXLlyhUmTJhnt1+l0lCpVij/++INnnnnGHCFYhfnz5/PVV19hZ2fHr7/+SuXKlTMdU6dOnfQxcno1a9bk+++/p2nTpnTu3Jmff/6ZAQMGULt27QKKvODonL14WKMbrsd+A8Dp4gbs4/8lzausypEJUThpUhJxOWGovZhcphmpxaqrGJGwFIl13klP4OwexOB8djkPn808TlwUDLMV1mratCnnzp1j2bJlbN++nWvXrgFQqlQpmjVrRpcuXXB2tszM3cNDKVyYmJj9OpwJCQkAeHp6PtU9lixZwltvvQXAzz//TNeuXZ/4Gp06daJ27docPXqUNWvW5DmB0z+LrAQGBqLVavH393/ieMymxUdwYgFoU9HotBSNXARtJ+V+nsgzi3reokBk+8wPzIJkw+x6pxafyL+PQiLfz9GvIxycnL6smueJuXg2/0CZdCaeSH5738CMCRyAk5MTPXr0yHImpyUrV64cAFevZj9VWr9Pf+yTWL58OT169ECr1TJr1qz0RO5pVK9enaNHj+aYlFk9r0Co2RWO/am0w3+HZiOU6e1CCNPRpsH+Hw3tUvWgbBP14hGWRaOB4A9h2eNizncvwZnVyhhJUeDMMonB2ulfWcbExGQ7SeHQoUMARjXi8mLlypV0796dtLQ0Zs6cyTvvvJOvWGNiYgBDr2Gh1STDANrUhxCW9zF/Qog8OrsW7l42tJt8oPzSFkLvmQ7gnWEIy+4fpEanSiSBy0JgYCANGjQAYOHChZn27969m6tXr+Lk5ESbNm3yfN01a9bw6quvkpqaysyZM+nfv3++4rx+/Tq7dilrFAYFBeXrWhYv4Bmo3MrQDpsNyQnqxSNEYbR3muF77zJQrZ16sQjLZO9gXFIm6qjxWrmiwJj1FWpUVBSrVq0iIiKCe/fuZVkXTqPRMGeO5U1H/vTTT+nYsSMTJkzg5ZdfTu9pi4mJYcAApRr5wIED8fIyLCuzYsUKRo4cSalSpdiyZYvR9datW0eXLl1ITU3lp59+ol+/fnmKY8qUKbz++uuZCiMfP36c3r178+DBAypWrMgrr7ySn49rHYI/hHP/KN8/jFNepTZ6N8dThBB5dOUAXDtoaDcaIOsPi6zVfh22j4ck5Q0Qe6ZA+efVjckGaXQ5VdvNh2nTpjF8+HCjpaH0t9I87pLX6XRoNJosa8VZgg8//JCpU6fi6OhIaGgobm5ubNmyhbi4OIKDg9m0aZPRYvZz586lT58+lC1blsuXL6dvv337NmXKlCE5OZnAwEBCQ0Ozved3331nlKx5e3uTkJBA7dq1KV++PHZ2dly4cIHw8HC0Wi1lypRhw4YNVK9umlli+kkMN27cMMn1TEqng19awnXl9TU+5eCDIzKANh/05WJkkLrtyPaZ//UGnFmjfO/sBUNOgVMhH5phI8zy//n2CUoSpzfgAPhXM931Czl9Oa/8jF83y59XW7Zs4cMPP8TT05Nhw4axY8cO9u3bx6xZs4iMjGT58uVcvnyZwYMHU6tWLXOEYBJTpkwhODiYGTNmsHfvXlJSUqhYsSIjRoxgyJAhFClSJE/XSUpKIjk5GVAeVk7rv44ZM8Yogfvss8/Ys2cPp06dYtOmTSQmJuLp6UmTJk145ZVX6N+/f+Ef/6an0UCTgYZF7u9ehoh1UF1e8wiRLzEX4EyGRevr9ZHkTeSswduw63tIU363sf9HaD9V3ZhsjFl64Nq3b8/ff//N3r17adiwIX369GH+/PnpPW3Jycm89957LF++nCNHjlChglT4thQW3QMHkJYKU+tA/ONlxMo0hrc2qBuTFZMeONuT5TP/+yM4+LPyvZ0DDD4BniVViE6Yg9n+P189yLBih70TDD0Nbrmvgy5M0wNnlkkMYWFh1K1bl4YNG2a538nJiZkzZ+Ls7MxXX31ljhBEYWXvYDzu7co+uHZYvXiEsHZJsXD0D0O7ZldJ3kTeNBpg+D4tGQ5a3nj2wswsCdzdu3epWLFietvR0RGABw8epG9zcnIiJCQk02B/IXJV500okuH1zr7p6sUihLU79CtkWDaLxu+rF4uwLv7VoNKLhvbBnyHloXrx2BizJHC+vr5Gqxj4+PgAcOXKFaPj0tLS0uuYCZFnzp5Qr5ehfXoVxF3J/nghRNZSk41rKlZoDsVrqhaOsEIZE/7EO3BiiXqx2BizJHBlypQxWsXg2WefRafTsXatYZBsQkICu3btynFhdSGy1bA/aB7PPtWlKcv/CCGezIklkHDL0M5Y30uIvKjQHPxrGNr7Zkhh3wJilgSuWbNmnDp1ilu3lB8Mbdu2xc3NjU8//ZThw4czbdo0mjdvTmxsLK1atcrlakJkwbsMPJOh9t2R+fDwnnrxCGFtdDrYm2H4gf8zUDH7EkdCZEmjMe6Fu3MGLsjQqIJglgSua9euNG/enKNHjwLKK9Xvv/+e1NRUvv/+ewYPHsyRI0coW7YsX375pTlCELag8UDD98n3lMK+Qoi8Ob9F+WWr1/h9WTZLPJ2aXcA9wNDeN0O9WGyIWerANWjQgE2bNhlte+edd6hXrx5LliwhNjaW6tWr06dPH6OVDIR4IoH1oHQjuLpfae+fCUH9pHq8EHmxL8OyWe4ByuxTIZ6GgxMEvQNbv1baF7bCrdPKEojCbMzym+7KlStoNBpKly5ttL1u3bpPvPi7EDlq/L4hgYu/AmfXQI2O6sYkhIWzj4mAi9sNG4L6Kb+EhXha9d6CnZMg9XG1if0z4BXpiTMns7xCLVeuHN27dzfHpYUwVq0teJc1tKXrXohcuZ5YYGg4uED9t9QLRhQObkWh9muG9vHFkHBbvXhsgFkSOE9PT8qXL2+OSwthzM7euJjktYNwNUy9eISwcJqHd3GOWGnYUKsbuPqqFo8oRIwK+z6Cg7+oF4sNMEsC98wzzxiVERHCrOq8AU4ZxlLunZb9sULYOJfTi9Ho168ECOqvXjCicClWGaq0NrQP/gIpD7I/XuSLWRK4d955hz179nDw4EFzXF4IY07uUL+3oX12rRT2FSIraam4nMgwW7t8MxloLkwrY0mRpBgp7GtGZkng+vTpw4ABA3jppZf45ptviIiIIDk5OfcThXhaQf0yFPbVQtjP6sYjhCU6uwb7xJuGdqP31ItFFE7lQiAgw2oe+3+Swr5mYpYEzt7enhkzZnDv3j1Gjx7NM888g6urK/b29pm+HByk5IMwAa9AqN7O0D4yDx4lZn+8ELYo44olPuWg8kuqhSIKKY1GWSlH7/YpuLxLvXgKMbMkcDqdLs9fWq3WHCEIW5SxN+FhPBxbpF4sQliaG0fhyj5DO6i/MglICFOr2RVcixra+39SL5ZCzCwJnFarfaIvIUyidEMoUdvQPjBLuu6F0MvQ+6Z1dIM6r6sYjCjUHJ2hXh9DO2IdxF5SL55CyiwJnBCq0Gig4buGdnSEUhFcCFuXcBtOLk1vPqzWCZxlFRxhRg3eBjv9ECmdlBQxA0ngROHybCdw8ze0D0jXvRAcnqvU5XrswbNvqBeLsA2eJeCZDob2kQWQnKBaOIWRSRK4b775hr///jtf1/j777/55ptvTBGOsGUOTsZV5c9thOjz6sUjhNpSjQuqJpdpRppPBRUDEjYj4xuR5Hg49qd6sRRCJkngRo0axbJly/J1jaVLlzJ69GhThCNsXf23wM7R0A6blf2xQhR2p1dCwq305oPneqoXi7AtpRtAqXqG9oGfQMa9m4y8QhWFj0cAPNvZ0D66UJmVKoQtyjiMoGhlHpVuql4swvY0zFAdIOY8XNiiXiyFjMmKsC1dupTt27c/9fnR0dGmCkUIaPQuHH9cRuRRAoT/AY0H5HyOEIXN1YNw/bCh3bA/aOTvdlGAnnkFNo6ChMcFpPfPhMovqhtTIWGyBC4hIYGEhPwNUNRoNCaKRti8knWgdCO4ul9ph81SfnlJ3SthSzL2vjl5Qa3XID5JvXiE7XEoosxI3fa10r6wBe5Egl8VdeMqBEySwF26JPVdhAVq9K4hgbt7GSL/gWptVA1JiAJz74Yy/k2v7pvKusFIAicKWP0+sPNbSHu8pOaBn+D/vlc3pkLAJAlc2bJlTXEZIUyrWjvwDIR715T2gZmSwAnbcehX0KY+bmiUXhAh1OBWDGp2gaN/KO1jf0LoaHDxUTcuKyeDIUThZe8AQRl+aV3aCbdOqRePEAUl5SEc+s3QrtoGfMurF48QGUuKpCRB+O/qxVJISAInCre6vcDBxdCWwr7CFpxcBkkZJoZlXFxcCDWUeA7KBhvaYbNBm6ZePIWAJHCicHP1hedeNbSPL4bEGPXiEcLcdDpluICe/zNQ/nn14hFCL2MvXNwVZY1U8dQkgROFX8YfGqkP4chc1UIRwuz+3Qs3TxjaDfsr6wQLobaqbcCrjKG9X96I5IckcKLwC3gGyjcztMN+gbQU9eIRwpwyDhNw8YGar2Z/rBAF6b/jkv/dDVHH1YvHykkCJ2xDowzVwO/fgDOr1YtFCHOJuwJn1xra9XpDEVfVwhEik7o9wTHDv8kDstTh05IETtiGyq3AJ8MsPPmhIQqjg7+A7vFakxp7KR0iLI+LD9TqbmifWAKJshLT0zBJAjdr1iwuX75siktZnCVLltC8eXN8fHxwc3OjVq1aTJw4kZSUp3sFd/jwYbp27UpAQADOzs6UL1+eDz74gNu3b+d43q1btxg4cCDly5fHycmJgIAAunbtypEjR54qDptjZ2c8E+/qAbgu/+1EIfIoEQ7PM7SrtwOvQPXiESI7GcclpyUbl7wReWaSBO69997j888/N8WlLMrgwYN59dVX2bNnD0FBQbRu3ZorV67wySef0KJFCx48ePBE11u6dCmNGjVi6dKllC1blldeeQU7OzumT5/Oc889x/nz57M8LzIykueee44ZM2ZgZ2dHhw4dKFu2LEuXLqVhw4asWLHCFB+38Kv9OhRxN7SlpIgoTI4vhodxhnbGX5JCWBK/qlCxhaF9aI6MS34KJnuFmpRUuJZnWblyJVOmTMHd3Z0DBw7wzz//sGzZMs6dO0fNmjXZvXs3o0ePzvP1bty4Qa9evUhNTWXWrFmEhYXx119/ERkZyRtvvMGtW7fo0aMHOp3O6DydTkf37t25ffs2b775JpGRkfz111+EhYUxa9YsUlNT6dmzJzdv3jT1f4LCx9lTSeL0Ti6H+7fUi0cIU9HpjIcFlKgFZRqpF48Qucn4B8b9KDi9Sr1YrJTJErhbtzL/Inz33Xdp1qxZFkcbLF26lF69ejFu3DhiYiynPtc333wDwIgRI6hbt2769mLFivHjjz8CMH36dOLj4/N0vR9++IGkpCRatmxJv3790rfb29szc+ZMvLy8OHjwIBs3bjQ6b/369YSHh+Pt7c2PP/6Ivb1hMfZ+/foRGhpKQkICU6ZMeerPalMyvkbVpijLDQlh7S7tgDtnDO2G70rpEGHZKr0IvhUNbXkj8sRMlsBFRkZm2hYbG8vu3buzfdW4ZcsWXn31VRYsWMDo0aNp2LAhUVFRpgrpqV2/fp2DBw8C0KNHj0z7mzZtSunSpUlOTmbdurwVItS/5szqeu7u7rRv3x6A5cuXZ3le+/btcXd3z3Su/nr/PU9ko2hFqPySoX3oV0hNVi8eIUwhYz0tNz94trN6sQiRF/8dl3ztIFw7rF48VshkCVx0dHSmMVz6gfk3btzI8pwNGzYAMGbMGGbPnk1UVBRffvmlqUJ6auHh4QD4+vpSvnzW6wfWr1/f6Nic3L9/P/2/jf68vF5P387tvHPnzpGYmJhrLALjrvvE23BKxhAKKxZ7ESI3GNr1+oCDk3rxCJFXtV6DIh6GtvTCPRGTlhFZsGCBUfvMGaVL//jxrAv1RUVFYW9vz+eff87bb7/Nm2++mekVohouXboEQJkyZbI9pnTp0kbH5iTjDN3srpnd9XKLRX+eTqcrtDOBTa5iCyhWxdDeP1MZQySENQr7GXj879fOERr0VTUcIfLM2RPqvGFon1oB92U8d145mOpCFStWZNKkSbRt25agoCC2bdvGnTt3cHR0ZO3atXTs2DHTOXfu3MHNzS29PX78eF555RVThfTU7t+/D2AU23/pX2feu3cvz9fL6ZrZXS+3WDK+Vs1LLACBgdmXFoiKisLPzy/XsibWzuWZHnjsHKM0oo4Se/wfUkvUzemUQikuLk7tEEQ+aB4lUPTI/PS/xB9WfJl7D+zgQfb//8oztz2W/MztK3XG98BPaNCBNoXEHdNIDPpQ7bDMTqvVYmeXvz40k/XABQcHU7duXZ5//nlatWrFq6++ir29PYMHD2bRokVcvHjR6PjU1FQOHTqU3oME4OPjw8svv2yqkITI1sOqHdBm6Lp3PTFfxWiEeDrOESuwe5SQ3k56rqeK0Qjx5NK8yvKobPP0tsupP5XacCJXJumBO3v2LHFxcQQGBtK5c2c2bdoEwEcffcSgQYP48ccfefnll1m0aBF16tQBlFmed+/epVOnTqYIwaQ8PJRf7DmNKUtIUH5oenp65vl6+mt6eXnl+XoeHh7ExsZmG4v+vLzGAnDt2rVs9wUGBqLVavH398/TtaxavV6wbzoAzhc24OyUAl6lVA5KHTbxvAsbrRZO/WFoBzbAt+aLeT5dnrntsdhn/vwgWLANALsHMfjf2g21X1M5KPPKb+8bmCiBq1LFMJ5o3759XL58Ga1WS4UKFQD4/vvv6d+/P/Xr10/fdvHiRTQaDW+99ZYpQjCpcuXKAXD16tVsj9Hv0x+bk7Jly6Z/f+XKFWrWrJnn65UrV47Y2FiuXLmSYxwajcboPiIPgvrB/h+VpYd0aUoxydDCV5BaFFIXtkDsBUNbCvcKa1XhBfCrBnfOKu0DM5XltqQUTo7MshZquXLl0hM1gHfeeYelS5dSqVIlLly4wIULF7C3t2fcuHE0btzYHCHki76XMCYmJttJCocOHQIwqhGXHU9PTypVqmR0Xl6vp2/ndl7lypWzLDMicuBTFqq2MbQP/QYpT7a6hhCq2T/T8L1HCXhG/fHDQjwVjca4pEjUMWW5Q5GjAlvMvlOnTkRERHD27Fl27NjBtWvXGDFiREHd/okEBgbSoEEDABYuXJhp/+7du7l69SpOTk60adMm0/6s6CdxZHW9hIQE1qxZA5DplbL+vNWrV2f5GlV/PUt8FW0VMvZaPIiFE0vVi0WIvLoTqfTA6dXvC/aO6sUjRH491w2cMwwvyvgHisiSRvfftZsEoCyl1bFjR9zd3dmxY0d6T1hMTAwvvPACJ06cYNiwYXz33Xfp56xYsYKRI0dSqlQptmzZYnS9GzduULlyZZKSkpg9ezbvvPMOAGlpafTp04cFCxbQoEEDDhw4gCZDt7FOp6NevXqEh4fTs2dPfv311/TVGGbPnk3//v1xd3fn3LlzFC9ePN+fOzAwkEepaew+lrkwc6Gk01F2yUs4xSglb5KLVuffrhttpus+JiYWgKJFfVWORDwJ/52f4X1KWbhea+/EpTfDSHMpmqdz5ZnbHmt55sX2fY3vUaUWnE5jz6U39pHqXlLlqMzj+TrVcLDT5DgmPTeSwOXgww8/ZOrUqTg6OhIaGoqbmxtbtmwhLi6O4OBgNm3ahIuLS/rxc+fOpU+fPpQtWzbLmmxLlizhtddeIy0tjYYNG1KuXDkOHjzIxYsXCQgIYPfu3emvWjOKiIggJCSEO3fuUKFCBRo0aMClS5cICwvDwcGBxYsXZ1mm5WkEBgZyM/4hge/PM8n1rMGr9tuY6Phzertb8mgO6KqrGJEQ2fMkkX1OA3HTKDP1lqQ+z/BUGf8mrF8p7rDTaTD2GiUt+TG1PRNTu6sclXlcm9GL4l7O+UrgCuwVqjWaMmUKf/31F40bN2bv3r2sW7eOwMBAJkyYwNatW42St7zo2rUrBw4coFOnTly8eJEVK1aQlpbG+++/z7Fjx7JM3gCqVq3K8ePHef/990lLS2PFihVcunSJTp06ceDAAZMlb7ZqVVowsTrD+ME+DhtyOFoIdXW1356evAH8ltZavWCEMKHr+LFRa1h16DX7rTgjJUWyIz1wwogt9sABDHdYxPsOqwFI02lo9ugHrun8VI5KCGN2aNleZAhl7O4AcEBbjW6PZOa0KDyCNGdY7DQ2vf1Jyjv8lfaCihGZhyl64CSBE0ZsbgzcYw4JNyj/e2M0ujQAYmv1J7rJaJWjMj9rGRsjFG6XNlBqw9vp7RsvzSKhYtsnuoY8c9tjVc9cp6PMklY4x5wGINm3Kv++urnQjUs2xRg4ky2lJQoPBzsNVQI8cj+wMAmoCs+0T1/Y3jdiEb5tv4Ai2S+nVhjc1ihlU/z9bex5W6v1GXrGvUpTsmEXsH+yH+PyzG2P1T3zpu/DqvcBcIqNoMqDo1D+eXVjMjEHu/wnpDIGTgi9hu8Zvn8YD8cWqReLEP918yRc3mVoB73zxMmbEFbh2S7gmmFW9f6f1IvFgkkCJ4Re6SAoUdvQPjALZISBsBQHMtTFcnSFurLuqSikHJ2hXh9DO2Id3L2sWjiWShI4IfQ0GmiUoRcuOgIublMvHiH0EqPh+BJDu1Z3cPFRLx4hzK1BX7DT9zDrIOznHA+3RaolcBUqVKBatWp8++23JCfLNGFhIWp0BLcMCz5L172wBId/g7QMPydl3VNR2HmWNF4e7sgCSE5QLx4LpFoCd/nyZSIjIxkxYgSVK1fmt99+UysUIQwcnKD+W4b2uX8g5kL2xwthbmkpcHCOoV2xBfhVVS8eIQpKxnHJyfFw7E/1YrFAqiVwv/32G3PmzOG9997DycmJt99+O/eThCgI9d8CuwzrSobNVi8WIU6vgvtRhnbGX2pCFGaB9aFkXUP7wCzQatWLx8KolsD16tWLPn36MH36dM6dO8fFixfVCkUIYx4B8GwnQzv8D3h4T714hG3LuKh30UpQqaV6sQhRkP47LjnmHFzcql48FsZiJjGULVtW7RCEMGjY3/D9o/twdKF6sQjbde0QXD9kaAf1BzuL+bEthPk90wHcAwxtGZecTn4SCJGVUvUgMMjQDpOue6GCjL1vTp5Q+zX1YhFCDQ5FoH5fQ/v8Jog+r148FkQSOCGyk7EXLvai8oNDiIJy7wacXmlo13kTnKykkr4QplS/D9gXMbTDZqkXiwUpsDLe4eHhrF27ljNnzhAbq6zL5uvrS/Xq1Wnbti1169bN5QpCFLBnXoGNo+H+DaV94Ceo0krdmITtODgHtKmPGxpl5QUhbJG7Pzzb2TAL9ehCaDEKnL3UjUtlZu+BS0pKomvXrtSrV49vv/2WyMhI0tLSSEtLIzIykm+//ZYGDRrQuXNnkpKSzB2OEHln76gUk9S7sBXuRKgXj7AdKQ+V2m96VduAb3n14hFCbUbjkhOUyWU2zuwJ3PDhw9m1axfLli3j7t27HDp0iE2bNrFp0yYOHTpEXFwcy5YtY+/evQwfPtzc4QjxZOr1BnsnQ/uAdN2LAnBiCSTFGNqNpHCvsHEl60DpRoZ22CzQpqkXjwUwewK3ePFiJk+eTMeOHbG3t88cgJ0dHTp0YNKkSSxevNjc4QjxZNyKwXNdDe1jf8KDu+rFIwo/nU55Xa8X8CyUC1EvHiEsRcY/ZO5ehnMbVQvFEpg9gXv48CFFixbN9ThfX18ePnxo7nCEeHIZly1KSYLw39WLRRR+l3fDrZOGdsP+Sj0sIWxdtf8Dz1KGdsZZ2jbI7JMYQkJCGDNmDPXq1cs2kYuJiWHs2LGEhMhfmcICFa8JZZvCv7uVdthsaDQA7DL3KAuRbxl731x8oWbX7I8VVkWn05GWloZOpyvQ+6alKa8aU1JSCvS+ZtHgPdj1nfL99WMQdRqKVVYlFI1Gg729PRqV/sAyewI3bdo0XnjhBcqUKUOLFi2oXr063t7eAMTFxXHmzBm2bduGr68v8+bNM3c4Qjydhv0NCVzcFYhYD9X/T92YROFz9zKc/dvQrt8HHF1UC0eYRkpKCtHR0dy/fz89mSro+4PyO9fq+TaHRn6ge1yX8/JViFOvh1qj0eDq6oq7uzs+Pj4FmsyZPYGrWLEiJ0+e5KeffmL9+vXMmzePu3eVMUQ+Pj5Ur16dzz//nP79++PlZdtTgoUFq9oGvMpA/BWlfeAnSeCE6YX9DDzunbFzgAayRrS1S0lJ4cqVKzx69Ei1GBwcCqximPnZOYB/9QwldlDGjarUC6bT6UhMTEz/KlWqFHYFtFpKgTxVT09PPv74Yz7++OOCuJ0QpmfvAEFvw6bPlfblXXDzJBR/Vt24ROGRfB+OLDC0n3kFPEuqF48wiejoaB49eoRGo8Hf3x83N7cC+wWvl5qqJDuFJpFLfQjR5wxtd29w9yvwMHQ6HcnJydy/f5/4+HgSEhKIjY2lWLFiBXJ/VZ9mbGwsCQkJlClTRs0whMibuj1h+wRlIgMovXCvTFc3JlF4hP8ByfGGdsP3sj9WWAWdTsf9+/cB8Pf3x9fXV9V4HB0dVb2/yTg6gouHUg8O4NFdcCihSi9ckSJF8PDwwN7entjYWO7evUvRokUL5FWqyf8MSE5OZsaMGbz00kuULFkSFxcXSpYsSePGjRk/fjw3btxIP3bYsGFUqFDB1CEIYR4uPlCru6F9YgkkxmR/vBB5pU2D/T8a2oFBULqBevEIk9AXrQdwc3NTOZpCxi1Dj1vaI3gYn/2xBUA/BCw1NRVtAa2bbdIEbt++fVSqVIlBgwaxefNmbt68SXJyMjdv3uTAgQOMGjWK6tWrM3OmYepvQc/GESJfgjJUA099CId+VS8WUXicXQtx/xrajd9XLxZhMhl/vxX0a9NCz9nLeH3UhNvqxQJGdW4LKoEz2SvUvXv3EhoaSnJyMhUrVqRHjx7Url0bT09P7t69S3h4OEuWLOH8+fMMHDiQ8+fPm+rWQhQc/2pQsYWyrBYoJUWafACOzurGJazbvhmG773LQvV26sUihDXQaJReuHvXlXZKIjxKhCK209NpkgTuwYMHdO3aleTkZL766is+++yzTO9/u3Tpwtdff82vv/7KoEGD+OGHH6RLWVinxgMNCVzibeVVat031Y1JWK+rB+HqAUO70XtSY1CIvHAtCvdvgu5xaZaE2za1ZrBJ+nRnz55NVFQUw4cPZ9SoUdkO3tNoNPTt25c9e/YQEBBAQkKCKW4vRMGq2AL8axja+6Yr09iFeBr7MkyEcfKEOm+oF4sQ1sTOXkni9B7GQWqyauEUNJMkcKtWrcLd3Z0vvvgiT8fXrl2bBQsWqFa9WIh80WiMxyjdOQvnN6sXj7Bed/+FM6sN7Xq9wMlDvXiEsDZufkCGXCLxjmqhFDSTJHCnTp2icePGuLq65vmc0NBQ1q1bx6+/yiBwYYVqdgH3AEN77zT1YhHW68AsQ0V5jb3xurtCiNw5FAEXb0M7Kca4yG8hZpIELi4uDj+/Jy+i16pVK3r16mWKEIQoWA5OyvJaepd2QNRx9eIR1udhPByZb2jX6AhegerFI4S1cvM3fK/T2kx5J5MkcN7e3ty58+Tdlhs3bmT+/Pm5HyiEJarXBxwz9DpnnEkoRG6OLIBH9w1tKR0ihJFatWqh0WhwcnIiJiaHpKyIKxRxByD85FnsvUvywcCBT3Sv+Ph4ihYtSsOGDa2mvJlJErgaNWqwb98+kpKS8nzOtm3baNOmDW+99ZYpQhCi4Ln6Gg84P7kU7t3I/ngh9NJSlZU89Mo0gVJ11YtHCAtz8OBBjh9X3mo8evSI33//PecT3JVeuA9GTcTF2YnRH32Q6ZCEhAQcHBzw8PBIL7Cs5+XlxciRIwkLC7OajiWTJHDt2rUjISGBMWPG5On4kydP8uabb6LVaq0m0xUiS43eI30ArTZVGdMkRG7OrIL4q4Z2kyfrLRCisJszZw4ApUqVMmpny8mTpeu3s+fgUQb06oq/qy5TdYCwsDDS0tIICgoyKryrN3DgQPz8/Bg5ciTJyZY/m9UkCVz//v0JCAhg0qRJjB07NsekbMGCBTRu3JgbN25YdB24+/fv8+mnn1K1alVcXFwoVqwYbdu2ZevWrU98raSkJNauXcvAgQOpVasWHh4eFClShNKlS9O9e3f27NmT7bm9e/dGo9Hk+PXw4cP8fFSRH74VoPr/GdqHf4NkKY8jcqDTwd4MpUN8K0CV1urFI4SFSUpK4s8//wRg3rx5uLm5ceLECQ4ePJj9SRoNk39RzunbvQOkPoDk+0aH7Nu3D4DGjRtneQlnZ2d69OhBVFQUf/31V/4/iJmZJIFzdXVl8eLFODo6MmbMGKpUqcKYMWNYuXIlW7duZcWKFXz++ec888wz9O7dm8TERIYMGUKXLl1McXuTu337NvXr12f8+PHcv3+fdu3aUaNGDdavX0/Lli2ZNu3JZhwuXLiQdu3aMWPGDOLj4wkNDaV9+/Y4Ozvz119/ERISwrhx43K8RnBwML169cryK6u/JEQBapyhq/5hPITn0tUvbNvVA3DjiKHdaIAU7hUigyVLlnDv3j2qV69OaGgoXbt2BXLuhQsPD2fv/oM0qvccVSuVUzYmKstr6cuWjRo1CoBx48YZdYKsXLky/Tq9e/cGYMYMyx/TbLKltEJCQtiyZQvdunXjwoULjB07NtMxOp0ODw8PJk6cSP/+/enTp4+pbm9S/fr1IzIyktDQUFavXp1eHmXdunW0b9+ewYMH06xZM5577rk8Xc/R0ZG33nqLgQMHUqdOnfTtOp2OyZMnM2zYMEaNGkXTpk1p1qxZltd4++230/9hCQtTpiEENoBrj/863D8DGrwN9ib730sUJnumGL539obaPVQLRQhLpE/U9L/z+vTpw9y5c1m0aBGTJ0/GxcUl0zn6JKxlixaGjcn3IeUBrq6u9OrVi7/++ouHDx/SrVs3nJ0Nyx8GBQWlf1+7dm38/PwICwsjKiqKEiVKmP4DmohJf8MEBwdz/vx5fvnlF9asWcPx48eJjY3F29ub8uXL065dO/r06UPJkiUBJYGxtDFwp0+fZtWqVdjb2zNnzhyj2nZt2rShd+/ezJkzh/Hjx6d38eZG31P2XxqNhqFDh7Ju3Tq2bNnCggULsk3ghIVr8gEs7ql8H3cFzq5RykIIkdGdCIhYZ2jXf8um1m4UmT1MSeNKbN4nAD6N1BSlLpqDo+n/qCzj64qzo+l6kCMjI9m1axf29va8+aayRGFISAgVK1bkwoULLF26NH17RvrhTY1DXkAZl/w4t0i8TefOnXnxxReZP38+np6e/PnnnzkuJNC4cWNWr17Nli1beOMNy10ZxSRPc+nSpbRp0wZXV1ecnZ0ZOHAgA/MwhXfSpEl8+eWXpgjBZFasWAEoyWjZsmUz7e/Rowdz5sxhzZo1pKSk4OjomO971qlThy1btnD16tXcDxaWqdr/gU85uHtZae+dDs90UFZtEEJvz1TD9/ZOjyfBCFt2JTaJlybvVDuMp7ZxyPNUCTDd6iH64v6tW7dO7/3SaDT06dOHUaNGMWfOnCwTuPDwcACqP1sTXO2Vgr4ASXfBoySHDx9Gp9NRt27dXFeBqlGjBqtXr+bIkSMWncCZZAzcq6++ip+fH507d+aPP/7g3r17eTqvaNGiWSZJatL/I6hfv36W+/XbExMTOXfunEnuqb9OTl2127ZtY9iwYfTr14+RI0eyYsUKq5glYzPs7JWxTHrXD8GV/erFIyzPvRtwPMPA6No90ksfCCEgNTWVefPmAWQaMtSzZ0/s7OzYuXMnFy5cMNqXmJhIYmIioOQVRoV90UHiHQ4fPgxAvXr1co2jaFFlfdVbt2495ScpGCZJ4EaPHk2FChVYsWIFPXv2xN/fn7Zt2/Lrr78SHR1tilsUmEuXLgFQpkyZLPd7enri6elpdGx+nDhxgr///huAzp07Z3vc/Pnz+f777/n555+ZMGECnTp1onz58mzYsCHfMQgTqf26MqZJb88PakUiLNH+H0Gb8rihUV67CyHS/f3339y8eZOiRYvSvn17o32lS5emZcuW6HS6TEtwxsfHp3/v4eEBjs7g5GU4IDGaQ49nsGbXOZOR/nf83bt3n/ajFAiTvEL98ssv+fLLL4mMjGTp0qUsX76c9evXs2HDBvr378/zzz9P586d6dixo0UPCASlfAiQY4kTd3d37t27l+eexuwkJCTQo0cPUlNTadWqFe3atct0TK1atZgyZQqhoaGUKVOGBw8ecOzYMcaMGcPevXtp3749GzdupHnz5nm+b2Bg9sv1REVF4efnx+3bt5/mI9k8txo9cDv8o9KI3EDMmV2kFa2qblA5iIuLUzsEm6BJvkfRg7+m/8X8sGIr7qV5gAr/n8kzLzhpaWmkpKTg4OBAamrW63OW9HBk3cAm5o1DqxSttTfDbOeSHo6kpKTkfmAe/PzzzwB0794djUaT6bpvvvkmGzduZO7cuYwePTq9AkPG39exsbF4enqicSmKQ/LjxE6XxuFDSgL33HPP5RpvbGwsoBT3zetnS01NRafTkZqaSnR0dK7VIbRaLXZ2+etDM+mIxipVqvDpp5/y6aef8u+//7Js2TKWLVvG9u3b2bZtG4MGDaJhw4Z07tyZTp06Ua5cOVPeno8//pjVq1c/8Xm//PILTZs2NWksuUlJSaFr166cPHmSChUqsGDBgiyPGzJkiFHbw8ODF198kZYtW9KxY0dWrVrF4MGDOXr0aAFELXKT9FxPXI/9iiZVqc3nFv4z91p+p3JUQm0up/7ELiUxvZ1U+x0VoxGWxMnRnsoB7ma9h37VAUsuORUVFZX+Rqlnz55ZHvPKK6/g7e3NjRs3+Oeff2jTpg2glDJzc3MjMTGRmJgYPD090Tm6onV0xS4lifh797lw6TJeXl5UqlQp11j0y3YFBASY6NOZh9nqHJQtW5ahQ4cydOhQbt68yfLly1m2bBk7d+5k3759DB8+nDp16tCpUyc6depEtWrV8n3PGzduEBER8cTnJSQYCq96eCiDMfXv03M6Xt/N+qRSU1Pp3r07GzZsoGzZsmzduhU/P78nuoZGo+HLL79k1apVHDt2jKtXr1K6dOk8nXvt2rVs9wUGBqLVavH3l7E5T8cf6vaEsNkAOJ9bi/PLXykTHCyYPG8zSnkIJzIszVMuBN+aLdWL5zF55uaXkpKS3uPp4OBgkklv+aH2/XPyxx9/kJaWRq1atYzKemTk6OjIa6+9xsyZM5k3bx6vvPJK+r66deuya9cuzp07R5UqVZSNHsUh9iKnIi6i0+l4rkZ1ihQpkmssZ86cAaBBgwZP9N9Mo9Hg6OhIsWLFcj0vv71vYKIxcLkpXrw4AwYMYMuWLdy6dYtffvmF1q1bc/LkSUaNGkWNGjX47rv891L8/vvv6aVJnuSrdWtDFXR9r+CVK1eyvEfGV6dP04OYlpbG66+/zvLlyyldujTbtm176okc1atXT/8+p6RMFLAmH4Dm8V+6ujTjqvvC9hxflF5QFICmg1ULRQhLpR/Xllu9U3392LVr13Lnzp307S+88AJgWG0BACdPcHAm5fHr66SE+5mW18qK/hotMtaUs0AFksBl5Ovry1tvvcXff//N7du3+f333+nYsWOu03oLSt26yoLShw4dynK/frubm5shy8+jtLQ03njjDRYvXpyevJUvX/6pY9V384Kh51BYAO8yULOroR2+ABLuZH+8KLy0acalQwJqQsVQ9eIRwgLt2LGD8+fPA8qqCc2bN8/2a/jw4enj4zIuOt+hQwcANm3aZLiwRgPuAdR6pgo+3p4cPn6KoAb16dmzJ71792b79u2ZYgkPDycmJoagoCCLH7OvWqn4mJgY5s+fz5w5czh58qRaYWTSoUMHRo0axZ49e7hy5Uqm2agLFy4EoF27dk/UtarVaunZsyeLFi1KT94qVqyYr1gXLVoEKK9yq1a13IHyNin4Q6XnBSD1IRz4CUJHqxuTKHhn10JshpIHwR9KbUAh/iPjEllHjhzJ4cjM5w0bNgxQ6qk2adKEvXv3cubMGcMbKhcfvH2LsnbuFD7/biZHT0dy8LByj759+2a65ty5cwF4//33n/LTFJwC7YHT6XSsX7+eLl26UKpUKT766KP0d82WokaNGrzyyiukpaXRt29fHjx4kL5v/fr1zJ07Fzs7O0aOHJnp3J49e1KtWjWmTzd+ZabVaunTpw8LFy58ouTt6NGjrF69OtPsJa1Wy5w5c/j0008BGDRokEWPbbBJAc9AlZcN7YM/w8P8zVoWVkanM142y6uMrM4hRBbmz5//VMOfTp8+bXQd/aS/X375xbBRowF3f5o0qMXmv34i+sRWdA/i0el0hISEGJ3/8OFDFi5cSIkSJejWrZvZP3d+FUgP3KVLl/j111+ZO3cuN27cSF8+y1Jem/7X7NmzOX36NJs3b6ZixYqEhIRw+/ZtduzYgU6nY8qUKVmug3rlyhUiIiIy1b6bPn16eldvxYoVs1wnFqBatWqMGDEivX358mU6duyIj48PdevWJSAggLi4OE6ePJk+Ru+1117jiy++MNVHF6bUdAhErle+fxgPh+dC8CBVQxIF6NIOuH7Y0G4yUNbHFcKMunTpQnBwMLNmzeLjjz82zCJ1KQr3b4L2cWdIwi1wzjwJcdq0aURHRzN37lycnJwKMPKnY7afJsnJySxdupQ5c+awc+fO9IzZ29ub1157jT59+jB8+HB27rS8JUT8/f05dOgQ48ePZ9myZaxatQo3NzdatWrFRx99RGjok41h0deUAbJ8567XrFkzowSuVq1aDB48mEOHDnH27Fn27NmDTqcjICCALl260KdPn/Rp1MIClWkIZZrAlb1Ke98MaNgfHCz/B4MwgZ0ZJma5FoU6lrskjxCFxbRp06hfvz5jx441vA2zswM3P7gfpbQfJcCjRKN1iOPj45kwYQJBQUHZljGxNBqdiVeTP3z4MHPmzGHRokXExyvdlHZ2drRs2ZI+ffrQoUOH9Mz2hRdeYOfOnek1aoT69GVEbty4oXYohUPkRliYYUJDu6lQr5d68fyHvmCzlJQwsX/3wW+G2e2Efg4hw9SLJwN55gUnJSUlfXB+pUqVVBvqoi9Ga9NDbbSpcOsU6LRK29kLfCuY7PJP+qz1BfXzU0HCJD1wd+/e5ffff2fOnDmcOHECUMa7VaxYkT59+tCzZ88cq/8LUWhVfhH8a8DtU0p792RlyS15lVa47crQ++bsBQ2kcK8QqrJzANdihpI+D+OVGo2OzurGlQ8m+S1SokQJUlJS0Ol0uLm50bVrV/r06ZNpgKAQNkejUcbCLX9bad+9BCeXQS3LHyArntL1I3B+s6Hd8N0sx9sIIQqYuz8k3gEev3hMuGnxRdZzYpJZqI8ePQKULsHVq1fz66+/SvImhF6NjuCbYdbxzm+V+mCicNo1yfB9EXclgRNCqM/eURmPqvfgrlLmyUqZJIGrWbMmOp2Oa9eu0bJlS2rWrMmkSZO4deuWKS4vhHWzd4DnPzK0Y87BqRXqxSPM59YppfabXoO+4OqrXjxCCGPuAUCGChj3rTdPMUkCd+zYMcLCwujXrx8eHh6cOnWKjz/+mNKlS/PKK6+wYsWKTLXMhLApNbsad9Xv/A60WtXCEWaSsffNwRkaD1QvFiFEZg5FjP+oehALqcnqxZMPJivkW79+fX766SeioqL47bffCA4OJjU1lTVr1tClSxdKlizJkCFDOHbsmKluKYT1sHeEpkMN7Ttn4Owa9eIRphd9Dk4uN7Tr9VbG3AghLMt/e+ESrLMXzuQrMbi4uNCrVy927txJREREejG96Ohopk6dSt26dalXrx5Tp07l/v37pr69EJar1mvgVdrQ3jFReuEKk92TSR8cbecITaRosxAWycEJXH0M7STr7IUz61JalStXZsKECVy9epWVK1fStm1b7OzsCA8PZ8iQIYSHh5vz9kJYFociyoxUvVsnDSs1COt29zIcW2Ro13kdvEqpFo4QIhfuxTM0dJBwW7VQnlaBrIVqb29P+/btWb16NVevXmXcuHFUrFgRE9cQFsLy1XkDPEoa2jv+p6yZKazbjm9B93hmscYeggerGo4QIhcOTuCSYSxcUgykPlIvnqdQoIvZAxQvXpyRI0cSGRnJtm3beOMNWV5G2BAHJ+NeuKhjcG6jevGI/Iu5AMf+NLRrvwa+5dWLRwiRN+4BGRo6SLSusXAFnsBl1KxZM+bNm6dmCEIUvLpvGv/g2D5BeuGs2Y7/GXrf7Bzg+eHqxiOEyBtHZ3DOMBYuMQbSrKcXTtUETgib5OgCwR8a2jeOQISMhbNKdyLgxBJDu86bVl3ZXQib4/GfXjgrqgsnCZwQaqjXx3gQ7bZxMiPVGm2fYFgc276IxSxYL4TII0cXcPY2tJNirGZGqiRwQqihiKvx6gy3TsJpWZ3Bqtw6ZbyiRt1e4F06++OFEJbJo0SGhg7u31QtlCchCZwQaqnb07gu3LZvIE1WLLEa28eTXvfN3kl634QwgYiICDQaTY5fM2fONDonPDwce3t7Pvjggye6V3x8PEWLFqVh02boMo6FexALKZa/RqokcEKoxcEJmn1iaMech+OLsj9eWI6oY3Amw0oaDfqCZ4nsjxdC5ElcXBy9evWiUqVK6duKFStGkyZN6NixIz169KBVq1ZG53zwwQe4uLgwevToLK+ZkJCAg4MDHh4epKWlpW/38vJi5MiRhIWFMX/lZozXSI0y6ecyB0nghFBTrdegqOEHFdv/ZzXjL2zatvGG7x1cpO6bECYSGBjIhQsXuH79Op999hkRERHcuXOHPXv2sHz5cv744w8qVKiQfvzSpUvZs2cPAwYMwN8/66XrwsLCSEtLIygoCHt7e6N9AwcOxM/Pj5GjPifZ3t2w42EcPEoyx0c0GUnghFCTvQM0H2lox1+BI/PVi0fk7t99xitoBL3zn5lsQoin1aVLFw4cOMCWLVv4+uuvqVKlSo7HT548GYC+fftme8y+ffsAaNy4caZ9zs7O9OjRg6ioKP5avxNr6oWTBE4ItdXoBAHPGto7v7X4v/xslk4Hm78wtJ08pfdNCBO5dOkS+/fvp1mzZlkmW/8VHh7O3r17adSoEVWrVs20f8GCBWg0GkaNGgXAuHHjjMbSrVy5EoDevXsDMGPmLHDzM1wg+R48Ssz35zIXSeCEUJudHbzwmaGdcAvCZqsXj8hexHq4esDQDh4EbkXVi0eIQkT7uJTSkSNH2Lp1a67LbeoTsJYtW2a539XVlV69euHs7AxAt27d6NWrV/pXUFAQALVr18bPz4+wsDCi7mtBkyE1unfDYgutO6gdgBACqPoylKoH1w8r7V3fK7NUXX1zPk8UHG0abPnS0HYPgEYD1ItHFA4pD+HuJfPeI/Xx7HYHM/zK9ymvrGhgAhUrVmTYsGFMnjyZ0NBQ/Pz8qF69Ov7+/vj6+jJr1iyj47du3Qpk/WoUoHPnzrz44ovMnz8fT09P/vzzTzQaTZbHNm7cmNWrV7Nl+w7eaB8KCY9LiTxKgOT74Oxpks9oSpLACWEJNBpoOQbmtVPayfGw8zto/Y2qYYkMjv0Jd84a2s1HQBE39eIRhcPdS/BjI7PewtGcFx+wH/yrm+xygwcPxs7Ojh9++IE7d+5w584dAIKDgzMdGx4eDkD16tnf//Dhw+h0OurWrZtt8gZQo0YNVq9ezZEjR3ijR3dIigbt48T33nVw8lB+TlsQeYUqhKUo/zxUzjA9Pmw2xJr5L3ORNykPlDp9er4VlWWzhBAm8/PPP1OtWjWmTZvG66+/ztKlS4mMjCQhIYGdO3caHZuYmEhiojI+rWjR7IcxHD6svNWoV69ejvfWX+PWrVvKmsYZV8pJfais0GBhpAdOCEvy4pdwfpOyPJM2BbaOhS6/qh2VCPtZ+StcL/RzsDdrv4YQNmX58uX069ePihUrsn79eipXrpzj8fHx8enfe3h4ZHvcoUOHAKhfv36O1/P0VF6R3r17V9ngVhQS70Da47JO92+Ciw/Y2WdzhYInCZwQlsS/OtR5w1BK5OQyaPy+Mj5OqCMpFnZNMrRL1YNnXlEvHlG4+JRXXkOaUcrjMXCO5hoDZwIjRyrllBYtWpRr8gbg7e2d/v39+/fTE7D/ymsPnD4h9PF5vCKDxg48SxrGJ2pTlITOo3g2Vyh4ksAJYWle+AxOLIWUx6VENo6G3n9b3PgLm7F9vFLUU6/lGHkWwnQcnU06hixLKSmP72WZvcYRERFERkZSoUKFXHvK9FxdXXFzcyMxMZGYmJgsE7j4+HguXLiAl5eX0coOWYmJUV6RBgRkqOno7AWObpDyuJRIwi1wLWoxve8yBk4IS+NRHJpkWNPv3z0QuUG9eGzZ7bNwcI6hXbWNMlZRCGEyt27deqrz6tatC8Dp06ez3H/q1Cl0Oh3PPfdcjhMYAE6ePAn8p6dOo1F64fR0Wota6F4SOCEsUZNB4JZhWZiNoyD1kXrx2CKdDv75FHSP1060c4SXvlY3JiEKoapVq+Lg4MDFixf58ssvSdH3GObihRdeAAwrLfyX/jpJSbkXRtdfo0WLFsY7nNyVnji9pGhlUpMFkAROCEvk5A4vZFhiK+Y8HPhJvXhs0blNcGGLod2wPxStqF48QhRSAQEBfPXVVwCMGTOGgIAAXnrpJV5//XU++ugjNm/enOV5HTp0AGDTpk1Z7q9VqxY+Pj4cPnyYoKAgevbsSe/evdm+fbvRceHh4cTExBAUFESJEiUyX8izJEZLbMVfs4jivpLACWGp6vaCgJqG9o7/WVT3faGWlqL0vum5FoXnh6sXjxCF3MiRI9m9ezdvvfUWJUqUYP/+/SxatIhJkybx4osvsnjx4kzn1KlThyZNmhAWFsaZM2cy7ff29mbt2rWEhoZy8eJFFixYwLx58zItaD937lwA3n///ayDc3A2XmLrUQI8jM/62AIkCZwQlsrOHl7+n6H9KAE2f5n98cJ0Dv4CMecM7Rc+Axdv1cIRwhYEBwczZ84cTp06xb1790hJSWHMmDEArFmzJstzhgwZAsAvv/yS5f4mTZqwefNmoqOj0el06HQ6QkJC0vc/fPiQhQsXUqJECbp165Z9cB7FlfpweveuK6uzqEgSOCEsWblgeLazoX1sIVw9qF48tiAxRpl5quf/jNIbKoQoUHZ2dpQqVQog29IiXbp0ITg4mFmzZj3VZIhp06YRHR3N+PHjcXJyyiEYe+MJDWmPIPH2E9/PlCSBE8LSvfgVOLoa2uuHw+NFn4UZbP7c+PVIq3FgLxWXhChIsbGxvPvuuwwaNIhy5crx3nvvZXvstGnTePDgAWPHjn2ie8THxzNhwoT08XG5cvE1/ll8/xakJj/RPU1JErhs3L9/n08//ZSqVavi4uJCsWLFaNu2bfriuU+qefPmaDSabL+KF8+5OODmzZtp06YNxYoVw8XFhWrVqvHZZ5+RkJDwVPEIK+IVCCFDDe0b4XD0d/XiKcyu7IfwDP9tq/0fVGyR/fFCCLPYv38/GzZs4PXXX2f37t34+flle2ydOnVIS0tj+vTpT3QPLy8vYmJiOHDgQK5lRgClrIhXYIYNOrh344nuaUryZ2UWbt++TUhICJGRkZQoUYJ27dpx69Yt1q9fz/r165kyZQoffPBB7hfKQqtWrbJM1ry8vLI4WjF58mSGDh2KRqMhJCSEgIAAdu3axTfffMOyZcvYvXs3xYoVe6p4hJVo/AEcWQBx/yrtzWOU5MLVV9WwCpW0FFibIVF2dIXWE9SLRwgb1qZNGy5fvqx2GJkVcVMmNenXRn0YBw/vgb1LgYciCVwW+vXrR2RkJKGhoaxevRpXV6XLdN26dbRv357BgwfTrFkznnvuuSe+9ogRI2jevHmejw8PD2fYsGHY29uzZs0aXn75ZUCpa9O+fXu2bNnCu+++y9KlS584FmFFHJ2h9XhY1ENpJ8UoKzR0mKFuXIXJgVlw+5Sh3ewT8C6tXjxCCMvkUQIexBlqRMZfBe+CLzEkr1D/4/Tp06xatQp7e3vmzJmTnryB8hdB79690Wq1jB8/PoermM748ePR6XT06dMnPXkDZRmROXPmYGdnx7Jlyzh79myBxCNUVLUNVDH8G+Do73Bpl3rxFCbx140nLvhVg0YD1ItHCGG57B3BM0O9uLRHyjqpBUwSuP9YsWIFoExnLlu2bKb9PXooPSBr1qzJc7Xop/Xo0SP+/vtvo/tmVLZsWYKDgwFD3KIQ02igzbfK2nx6awerOoi20NgwQinTotd2EjgUUS8eIYRlcy1mPKEhMVpZaqsASQL3H+Hh4QDZLqir356YmMi5c+eyPCYnK1asYPDgwfTv35/PP/+cDRs2oM1mRmFkZGT6EiC5xaOPWxRy3qWhxWeGdsx52PW9evEUBqdXwZnVhvZz3aFcU/XiEUJYPo0GvMtgWKFBp/TEFeAKDTIG7j8uXboEQJkyZbLc7+npiaenJ/fu3ePSpUs888wzT3T9qVOnZtpWpUoVfv/9dxo0aJBlLN7e3nh4eGR5vdKlSxsdK2xAUH84/hdEHVPau79XasX5VVE3LmuUFAt/f2Rou/jAS09WikAIYaMcXcDdHxIe15/TaZW1Usm6Zp2pSQL3H/fv3wfAzc0t22Pc3d25d+8e9+7dy/N1Q0JCeOONN3j++ecJDAzk7t27HDx4kFGjRnHq1ClatmzJ/v37qV69+hPHAjxRLIGBgdnui4qKws/Pj9u31S1QKHLm0PRLfJZ2QqPTQtojUpb2427HP5Vik08gLi7OPAFaCY8tH+OSoRhnfPBnJCcBSYX337+tP/OClJaWRkpKCg4ODqSmpqoahzAT56I4JN2FNGUoi+7eDWIvn0Djlf3vWQCtVoudXf5eghaqBO7jjz9m9erVuR/4H7/88gtNm5r3lcl/Cwy6urpSqlQpXn75ZUJCQjh48CAjR45k5cqVZo1DFA6pfjV48FxvXI/9CoDjrXBcj80hqU4/lSOzHkX+3Y5LhGHsaHLZ5iRXbq9iREIIq6OxI82jJMQqb8E0Oi0eu8eR0GYGaMw7Sq1QJXA3btwgIiLiic/LWAxX/6oyMTEx1+M9PT2f+F7/5eTkxGeffUaHDh3YsGEDKSkpODo6mjWWa9euZbsvMDAQrVaLv79/nq8nVNJ2HFzdAbEXAHAPm4J7nc7gXz2XEzOzuef9MB52jTG0nTxx6vQj/l4BqoVU0GzumasgJSUlvcfTwcEh/We7WtS+f6Hl6EPKg3i4q9SGc4w6hP+/ayDonWxPyW/vGxSySQy///57+mK1T/LVunXr9GuUK1cOgCtXrmR5j4yvTvXH5pf+tWlycjLR0dGZYomLi0t/nfpfV69eNWkswooUcYUOMw1/5aU9ghX9lYK0ImcbPlUWo9Z76WvwKqVePEII6+ZRHMOEBmDT5xBzway3LFQJnCnUrVsXgEOHDmW5X7/dzc2NKlVMM2g8JiYm/fuMkxWqVq2aXocut3j0cQsbU6YhNMmwKkjUMZmVmpvTq4yXIivfDOrmYR1EIYTIjp092GcoPZSSBCvfA635xh9KAvcfHTp0AGDPnj1Z9sItXLgQgHbt2pmsO3rRokWA0hOnn5QAUKRIEdq2bWt034z+/fdf9u7dC0DHjh1NEouwQs0/Bb8Mr013ToTrR9SLx5LduwGrBxnaTl7wynSlJIAQQuSHnT24ZViz9eoB2Ju58oTJbme2K1upGjVq8Morr5CWlkbfvn158OBB+r7169czd+5c7OzsGDlyZKZze/bsSbVq1TItqLtt2za2b9+O7j/1YR49esSECROYNm0aAMOGDct0zREjRqDRaPjtt9/YsGFD+vakpCT69u1LWloanTt3plq1avn63MKKOTpDx5mgeTwDVZsKS99S1ucTBlotrHhXWbtQ7/++f1zLSQghTMCzJPhkWARg69dwLes3aPlVqCYxmMrs2bM5ffo0mzdvpmLFioSEhHD79m127NiBTqdjypQpWa6DeuXKFSIiIozGsQEcO3aMIUOGEBAQQO3atSlatCh37tzh+PHj3Lql1I/56KOP6Nu3b6Zr1q1bl0mTJjF06FDatGlDs2bN8Pf3Z9euXURFRVG1alV++ukn8/yHENajZB1oPgK2jVPady/B2iHQ+RfpXdLb/yNc2mFo13wVanZRLx4hROGjsYM238Hc1spaqdpUWNoH+u8CF2+T3kp64LLg7+/PoUOHGDFiBO7u7qxatYrjx4/TqlUrNm/ezKBBg3K/SAbNmjXjvffeo2zZshw/fpxly5axe/du3N3d6dmzJ3v27OHbb7/N9vwhQ4awadMmWrVqxfHjx1m1ahXu7u6MHDmSgwcPUqxYsfx+ZFEYhAyDshnK4ZxcCkf/UC8eS3I1DDZ/YWh7lYG236kXjxCi8CpZB1741NCOuwJrPjT5Kg0a3X/f6wmbpi8jcuPGDbVDEU/j3g2YGQwPYpW2oyv02w5+VbM8XF+wuVCXlEiMhlnPG2adauyg999Qtom6canEJp65hUhJSeH8+fMAVKpUSbUyHvp1u62ljEhiYiKzZ89m9erVnDp1iri4OLy9vdOHOPXr1y99gl92wsPDqV+/PgMGDEgfppQX8fHxVKhQgUqVKrF//340eXyDkelZ29vBgo7Gvf7/9wPU7wMYCurnVNYrN9IDJ0Rh4lkSOmZ4pZ6SBH+9Ybvj4bRpsPwd45IhLUbbbPImhKXbtWsXlSpVYujQoRw8eJBatWrRpUsXatWqxcGDBxkyZAiVKlViz549OV7ngw8+wMXFhdGjR2fal5CQgIODAx4eHplWqfDy8mLkyJGEhYUxf/78p/8gdvbQ6WfjSQ3rP4Hrh5/+mv+9hcmuJISwDFVaQaP3De3oSKU+nFarXkxq2TERLmw1tKu0huDBqoUjhMjenj17aNmyJTdv3mTo0KHcvHmTTZs2sXDhQjZt2sTNmzf56KOPiIqKokWLFtkmcUuXLmXPnj0MGDAgy57msLAw0tLSCAoKwt4+8/KDAwcOxM/Pj5EjR5KcnPz0H8gjwPgP6rRkWPQGJJhmqT5J4IQojF78Espk6GWKWKeUF7Elp1fDjgmGtncZ5YepCSqgCyFM68GDB3Tr1o1Hjx7xxRdfMGnSJKOyWqCs/f3tt9/y5Zdf8ujRI1599VWjShF6kydPBshyYiDAvn37AGjcuHGW+52dnenRowdRUVH89ddf+flYUKklNM9QteL+DVjcK3/XfEx+kglRGNk7wqvzwDPD6gLbx8PZderFVJCijim9jnr2ReDV+eDio15MQohszZkzh+vXr1OhQoUsX3tmNGrUKCpWrMiNGzf45ZdfjPaFh4ezd+9eGjVqRNWqxmN/FyxYgEajYdSoUQCMGzcOjUaT/pVxLfLevXsDMGPGjPx/uOc/hqptDe0re5Xl/PJJEjghCit3f+i2AOydDNuWvQ03wtWLqSDci4KF3ZXxf3rtpiozw4QQFmn16tUAvP7661m+1szIzs6OHj16ALBq1SqjffokrGXLlpnOc3V1pVevXjg7OwPQrVs3evXqlf4VFBSUfmzt2rXx8/MjLCyMqKiop/5cjwNWev+LZVi96VH2a5znldSBE6IwK1VPKVa76vGYuJRE+ONVeHsT+JRTNTSzeJQEi3ooryn0mg6B2q+pF5MQOUhOS+bqvatmvUdqaioADg6m/5Vf2rM0Thn/SHxKx48fB5TEKS/0x508edJo+9atypjXrF6Pdu7cmRdffJH58+fj6enJn3/+meMs08aNG7N69Wq2bNnCG2+8kae4suXsCd0Xws8tINk0k8okgROisKvzBkSfgz0/KO3E2/B7F+i7UdWwTC71ESx+E25kWEas2v9Bi8/Vi0mIXFy9d5WOq613KcQV7VdQyadSvq8TG6uUPvL09MzT8frj7t69a7Q9PFx5w1C9evVM5wAcPnwYnU5H3bp1cy0RUqNGDVavXs2RI0fyn8ABFKsMXX9T/og2AXmFKoQtCP0CanY1tGPOwZ/d0aTkvxvfImi1ysLR5zcbthV/DjrNlkkLQlgBb29vQKnDlhf6xM3HxzCuNTExkcRE5Wda0aJFszzv8GGljEe9evVyvYf+GvoVk0yiUktok33h/ichP9mEsAV2dvDKDCgXYth29QBe696FlMyzuKyKTgfrhysrT+j5lIPXl0ARN9XCEkLkXfny5QFDD1pu9MeVLWtYdzRj8ufh4ZHleYcOKeuS1q9fP9d7ZNfLl28N+oKTe+7H5UJeoQphKxycoNvv8FsbuH0KgCLX9+O9/j3otRwcnVUO8ClotUrydjDDTDT3AHhzJXgUVy0sIfKqtGdpVrRfYdZ7mHsMnCm0a9eOsLAwli5dytixY3N8vanValmyZEn6eXr6XjyA+/fvZ/k69kl64PQJYcZePpNxytur4pxIAieELXHxhp4rYW5bpcAvUOTaHlj0mpLcWVOPlVYLawfDkXmGbU5e8MZy8C2vWlhCPAkneyeTjCHLiTUspdWnTx8mTpxIREQE8+fPp1ev7Gul/fbbb5w/fx53d/f0ch+gzDJ1c3MjMTGRmJiYTAlcfHw8Fy5cwMvLi0qVcv9vHhMTA0BAQMDTfSgzk1eoQtgad3/ouRp8Kxi2XdgK89pDUqx6cT2J1EfKmLeMyZuzF/RcAcWfVS8uIcRTKVWqFP/73/8AGDBgQPqs1P86evQoH3zwAQDjx49PX1NUr27dugCcPn0607mnTp1Cp9Px3HPP5WmNU/0M17z01qlBEjghbJFnCei1hlTPMoZt1w/Br60hzrwlDfLtYTz80QWOLzJsc/FRktJSlvmDVgiRu/fee48hQ4aQlJRE69at00uC6G3ZsoWXX36ZBw8eMHjwYAYOHJjpGi+88AJgWG0hI31PZFJSUqZ9WdFfo0WLFk/0OQqKJHBC2CqvQOI6LSKlaDXDtugI+PkFuLJfvbhyEndFSTIv7TBscy0GvdZAydqqhSWEeHrR0dH07t2b3r17Exsbi5ubG1FRUbz44ots2LABgA0bNvDiiy9y8+ZN3NzcuHv3bvo50dHR6dfq0KEDAJs2bcp0n1q1auHj48Phw4cJCgqiZ8+e9O7dm+3bt2c6Njw8nJiYGIKCgihRooRZPnd+yRg4IWyY1tWPuA5/4LfpA2V5F4DEOzD3/6Dtd1Cvt6rxGTm3CZa/Aw8yzAjzrQhvLDV+HSyEsCoJCQnMmzcv03atVsvNmzcBuHnzJjqdDlDKhWQ8fsyYMRQrVgyAOnXq0KRJE/bu3cuZM2eM6sF5e3uzdu1aPv/8c44ePcrBgweBrNdMnTt3LgDvv/++aT6kGWh0+v8iQgCBgYFotVpu3LiR+8HC6t2+fRsAfx8PZbWGk8uMD3i2C7SdpEx+UEtaqrIo/c7/1E4q3RC6/wluWdd7EllLf+b+/ipHUvilpKRw/vx5ACpVqqTaJAJrmMRgSkuXLqVr164MHTqUSZMmPfH5Dx8+pHTp0jg6OnLp0iWcnHJfaeJJn7V+7N61a9eeOD49eYUqhABHF+g8Ryn4S4bBvSeXwsxguLRTnbhunoRfWmRO3mq9Bj1XSfImhMikS5cuBAcHM2vWrKcqwjtt2jSio6MZP358npI3tUgCJ4RQaDQQMhR6/KVMCtC7dw3mtYPl/eD+zYKJ5VEibB0Hs5tB1DHDdnsnZWH6DjOVpFMIIbIwbdo0Hjx4wNixY5/ovPj4eCZMmJA+Rs6SyRg4IYSxKq3gvX1KmY6L2wzbj/8FZ/+G4MHQsJ9StsPU0lLg2J9K8pbwn2SxWFXo/DOUqGX6+wohCpU6deqQlpb2xOd5eXml13+zdNIDJ4TIzLOEUhC39f/AMUNx30cJsO1r+KEmbP3adCVHHsTBvh9hSm1Y/YFx8qaxh5CP4N1dkrwJIcRj0gMnhMianR00eheeaQ+bPocTSwz7HsYr49J2fgeVQqF6O6jS+smWr0qKVcqBnFoJEeshLTnzMWUaw8v/k8RNCCH+QxI4IUTOPEtC51+gfl/YNg4u78qwUwfnNytfoJTzKFkX/KqCZylw9QU7B9BplYQt4RbciYBbJ+HmCeX8rBStDC2/gGr/p4zNE0IIYUQSOCFE3pRtDL3Xwr97Ye80iPwHdP8ZYxJ7Ufl6WqUbQpNBULWN0gMohBAiS5LACSGeTNkmyte9KDi2EM6shRtHnv56PuWhRkeo1V3puRNCCJErSeCEEE/HswSEDFO+7t+Ey7vhRrhS9iPuCtyPgrRHhuPtHJRlr3zLg/8zUOI5qNAcfMqp9QmEMKuMC6ZLzfzCLePz1RTQsA9J4IQQ+edRHGp2Ub70tFpIfai8ZtXpwMlDxrMJm2Jvb49Go0Gn05GcnEyRIkXUDkmYyYMHDwCws7PD3t6+QO4pg0yEEOZhZwdFXJXEzdlTkjdhczQaDa6urgDcv39f5WiEuWi1Wu7eVdZodnd3lx44IYQQwtq5u7uTmJhIfHw89vb2eHl5FVgPjV5qamqB3s9W6HQ6Hjx4wN27d9N74Ly8zFDgPBuSwAkhhBBm4uPjQ2JiIgkJCcTGxhIbG1vgMejHZxVUz5CtCggIwN3dvcDuJwmcEEIIYSYajYZSpUoRGxvL3bt3VekN09/T0dGxwO9d2NnZ2eHu7o6Xl1eBJm8gCZwQQghhVnZ2dhQrVoyiRYui1WrRarUFev/o6GgAihUrVqD3Lew0Gk36RBU1SAInhBBCFAD9L/yCHgOnv5/0wBUuMgtVCCGEEMLKSAInhBBCCGFlJIETQgghhLAyksAJIYQQQlgZSeCEEEIIIayMJHBCCCGEEFZGEjghhBBCCCuj0enX2BACKFKkCGlpaZQoUULtUEQB0BcUtbOTv+VshTxz2yPP3PJERUVhb2/Po0ePnvoa8jSFkdTU1AKvEi7Uc+vWLW7duqV2GKIAyTO3PfLMLY+9vX2+CyvLSgzCSMmSJQG4du2aypGIghAYGAjI87Yl8sxtjzzzwkl64IQQQgghrIwkcEIIIYQQVkYSOCGEEEIIKyMJnBBCCCGElZEETgghhBDCykgCJ4QQQghhZaSQrxBCCCGElZEeOCGEEEIIKyMJnBBCCCGElZEETgghhBDCykgCJ4QQQghhZSSBE0IIIYSwMpLACSGEEEJYGUnghBBCCCGsjCRwAoAlS5bQvHlzfHx8cHNzo1atWkycOJGUlBS1QxMmlJKSwpYtWxg+fDgNGjTA29sbR0dHihcvTvv27fn777/VDlEUgI8//hiNRoNGo+Hrr79WOxxhJo8ePWLq1Kk0bdoUX19fnJ2dCQwM5OWXX+avv/5SOzyRTw5qByDUN3jwYKZMmYKDgwMtWrTA3d2drVu38sknn7BmzRo2btyIi4uL2mEKE9ixYwcvvvgiAMWLF6dp06a4ublx+vRp1qxZw5o1a+jXrx8//fQTGo1G5WiFOezdu5dJkyah0WiQOu6F17Vr12jVqhWnT5+mWLFiBAcH4+bmxtWrV9m5cydubm5069ZN7TBFfuiETVuxYoUO0Lm7u+sOHz6cvv3OnTu6mjVr6gDdsGHDVIxQmNKWLVt0nTt31u3cuTPTvkWLFuns7e11gG7evHkqRCfMLTExUVe5cmVdqVKldB06dNABurFjx6odljCxpKQkXbVq1XSAbsyYMbpHjx4Z7U9MTNSFh4erE5wwGXmFauO++eYbAEaMGEHdunXTtxcrVowff/wRgOnTpxMfH69KfMK0WrRowdKlSwkJCcm0r1u3bvTu3RuA+fPnF3BkoiCMHDmSc+fOMXv2bLy8vNQOR5jJ+PHjOXv2LP369eOLL77A0dHRaL+rqyu1a9dWJzhhMpLA2bDr169z8OBBAHr06JFpf9OmTSldujTJycmsW7euoMMTKqhTpw4AV69eVTkSYWrbt29n2rRp9OzZkzZt2qgdjjCTlJQUZs6cCcDw4cNVjkaYk4yBs2Hh4eEA+Pr6Ur58+SyPqV+/PlevXiU8PJzXXnutIMMTKjh37hwAJUqUUDkSYUoJCQm89dZbBAQE8MMPP6gdjjCjI0eOEB0dTcmSJalUqRInTpxg+fLl3LhxAx8fH0JCQnj55Zexs5P+G2snCZwNu3TpEgBlypTJ9pjSpUsbHSsKr5s3bzJ37lwAOnfurG4wwqQ++ugjLl26xIoVK/Dx8VE7HGFGx48fByAwMJARI0YwceJEo8kq//vf/6hTpw4rV67M8We/sHySgtuw+/fvA+Dm5pbtMe7u7gDcu3evQGIS6khNTeWNN94gPj6emjVr0r9/f7VDEiayceNGZs2aRffu3enQoYPa4Qgzi4mJAZQ3LP/73/8YMGAAERERxMfHs2nTJqpUqUJ4eDht27aVMlFWThI4IQTvvvsuW7ZsoWjRoixdupQiRYqoHZIwgfj4ePr27Yufnx/Tpk1TOxxRAPS9bSkpKbz22mtMnz6dKlWq4OnpScuWLdm0aRPOzs6cPHmSRYsWqRytyA9J4GyYh4cHAImJidkek5CQAICnp2eBxCQK3ocffsicOXPw8fFJ/wtdFA6DBw/m2rVrTJ8+nWLFiqkdjigA+p/rQJY96WXKlKFt27YAbN68ucDiEqYnY+BsWLly5YCcZxzq9+mPFYXLsGHDmDp1Kt7e3mzcuDF9FqooHFasWIGDgwM//vhjelkgvbNnzwIwZ84cNm/eTPHixaVHphCoUKFClt9ndUxUVFSBxCTMQxI4G6b/ZR0TE8OlS5eynIl66NAhAKMacaJw+Pjjj/n+++/x8vJi48aN1K9fX+2QhBmkpqayY8eObPdfvnyZy5cvU7Zs2QKMSphL3bp101fZiI6OTp+IllF0dDRgGOMsrJO8QrVhgYGBNGjQAICFCxdm2r97926uXr2Kk5OT1I0qZEaMGMG3336Ll5cXmzZtSv93IAqXuLg4dDpdll+9evUCYOzYseh0Oi5fvqxusMIk9EvkQdavSFNSUtIT+qCgoAKNTZiWJHA27tNPPwVgwoQJHDlyJH17TEwMAwYMAGDgwIFStb0QGTVqFP/73//w9vaW5E2IQuiLL74AlBUZ9u/fn749NTWVYcOGcfHiRTw8POjTp49aIQoTkFeoNq5Dhw4MGjSIqVOn0qhRI0JDQ3Fzc2PLli3ExcURHBzM2LFj1Q5TmMjq1asZN24cAJUqVWLGjBlZHlesWDG+++67ggxNCGEioaGhjB07ltGjRxMSEkJQUBDFixfnyJEjXL58GRcXF/78808CAgLUDlXkgyRwgilTphAcHMyMGTPYu3cvKSkpVKxYkREjRjBkyBApKVGIxMbGpn9/6NCh9DGO/1W2bFlJ4ISwYqNGjSIoKIgffviBAwcOcPDgQYoXL07v3r355JNPqFatmtohinzS6DKWaBZCCCGEEBZPxsAJIYQQQlgZSeCEEEIIIayMJHBCCCGEEFZGEjghhBBCCCsjCZwQQgghhJWRBO7/27u/kKbaAI7jP01C27HOYK6/FEVJghRUSNKVo8CEhBopFKTRVVJRw1sr6SIko+iuq6igm1UYWARSF4kI6VCCioZedOWJgdvCOTDc3ouXxjvWa728ubOHvh/w5jxnj8+5+/KcnTMAAADDEHAAAACGIeAAAAAMQ8ABAAAYhoADAAAwDAEHAABgGAIOAADAMAQcAACAYQg4ADBIKpXSrVu31NTUJL/fr5UrV8rv96upqUm3b9/W/Py820sEUARl2Ww26/YiAADS3NycbNtWVVWVEomEVqxYkTc+PDystrY2OY4jj8ejxsZG1dTUKBaLaXR0VKlUSuvXr1c4HNaBAwdcugoAxcAOHACUiLdv32pxcVENDQ0F8TYyMqKDBw/KcRyFQiE5jqOhoSE9evRIQ0NDchxH3d3dmpmZUSAQ0MjIiEtXAaAYCDgAKBGjo6OSpMbGxrzj6XRa7e3tWlhY0JUrV3Tz5k1ZlpV3jmVZunHjhnp7e7WwsKC2tjal0+mirR1AcXELFQBc9vDhQ506depfx69du6aenh5t27ZN0Wi0YHfunzKZjGprazU9Pa07d+7o/Pnzy7FkAC5jBw4AXLZq1Sp1dHSosrJSktTe3q6Ojo7c35s3byRJJ0+eXDLeJKm8vFwnTpyQJD179mx5Fw7ANezAAUAJ+Pr1q2zbVnV1tRKJhMrKynJj69at05cvX/TkyRMdO3bsp3M9ffpUwWBQa9euleM4y7lsAC5hBw4ASkAkElE2m9WePXvy4k2SZmdnJUmrV6/+pbm+nxePx/OOT01NqaWlRZZlyefzqaurS6lU6jesHkCxVbi9AADA3wEnSXv37i0Ys21bsVhMyWTyl+b6Hm5erzd3LJlMKhAIaMOGDQqHw5qdnVUoFMrt7AEwCwEHACVgfHxckrRv376Csa1btyoWi2liYkLBYPCnc01MTEiStmzZkjt29+5dxWIxjY+Py+/3S5KqqqoUDAYViUR+GI4AShe3UAGgBCy1A3fkyBFJ0uPHj/Wzry1nMhmFw+G8z0nSixcvFAgEcvEmSa2trbIsS4ODg/97/QCKi4ADAJclk0lNT09rzZo12r59e8H46dOnVV1drU+fPunBgwdLznXv3j1NTU3Jsix1dnbmjn/48EF1dXV551ZUVKi2tlYfP378LdcBoHgIOABw2fv375XNZrVr166CBxgkaePGjerr65MkdXV16d27dz+cZ3JyMvfet+vXr2vTpk25sXg8Ltu2Cz7j9XpzD0kAMAcBBwAu+/btmyQt+UP0Z8+e1aVLlzQ/P6/m5ma9fv06b/zVq1c6fPiw0um0Ll68qHPnzi3rmgG4i4cYAMBlu3fvltfrVSQSUUNDg3bu3Kny8nJ1dnaqvr5e3d3duXM9Ho9mZmZ06NAhPX/+XM3NzXr58qVaWlqUzWbl8XgUj8dzt0/7+/vl8/nk9XqVSCQK/nc8HteOHTuKdKUAfhcCDgBcZtu2BgcHdfnyZU1OTmpsbEySdObMGc3Nzen+/fsFn8lkMrmX9DqOk3u4IZVK5Z1/9epV+Xw+1dXVFXzXbXFxUdFoVK2trct1aQCWCb/EAAB/gL6+PvX29urz58+qqamRJA0MDOjo0aMaGxv74etLAJQuAg4A/gCJREL19fXavHmzenp6FI/HFQqFtH//fg0MDLi9PAD/EQEHAH+IaDSqCxcuaHh4WJWVlTp+/Lj6+/tlWZbbSwPwHxFwAAAAhuE1IgAAAIYh4AAAAAxDwAEAABiGgAMAADAMAQcAAGAYAg4AAMAwBBwAAIBhCDgAAADDEHAAAACGIeAAAAAMQ8ABAAAYhoADAAAwDAEHAABgGAIOAADAMAQcAACAYQg4AAAAwxBwAAAAhiHgAAAADPMXP1PBetwMGe4AAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "ro.characterization.plot_pulse(pulse_ansatz, optimized_params);"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Multiqubit_Rydberg_Gates",
   "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.12.8"
  },
  "widgets": {
   "application/vnd.jupyter.widget-state+json": {
    "state": {
     "1abe1062e05d4382a5e7e49567965949": {
      "model_module": "@jupyter-widgets/base",
      "model_module_version": "2.0.0",
      "model_name": "LayoutModel",
      "state": {
       "flex": "2"
      }
     },
     "2e810a679bbf4d0abf10e32716aa60a2": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "2.0.0",
      "model_name": "HBoxModel",
      "state": {
       "children": [
        "IPY_MODEL_f0f69820e69a4aeba774df7edf942573",
        "IPY_MODEL_7995e7ba8a984230a22cd119cf8c5924",
        "IPY_MODEL_48dcc2e7d4f74916990aa98b77d778f1"
       ],
       "layout": "IPY_MODEL_8793459176ba45feb7ecaf6537ab9355"
      }
     },
     "3e643e4ffcc64f2aa0e53e0433bde4aa": {
      "model_module": "@jupyter-widgets/base",
      "model_module_version": "2.0.0",
      "model_name": "LayoutModel",
      "state": {}
     },
     "48dcc2e7d4f74916990aa98b77d778f1": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "2.0.0",
      "model_name": "HTMLModel",
      "state": {
       "layout": "IPY_MODEL_6d9282e02b8d4b338f5d556dcaf54561",
       "style": "IPY_MODEL_8effff66a89f4b03a832dba64ac38315",
       "value": " 1000/1000 [00:03&lt;00:00, 277.11it/s, infidelity=7.57e-11, converged=1/1]"
      }
     },
     "6d9282e02b8d4b338f5d556dcaf54561": {
      "model_module": "@jupyter-widgets/base",
      "model_module_version": "2.0.0",
      "model_name": "LayoutModel",
      "state": {}
     },
     "7995e7ba8a984230a22cd119cf8c5924": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "2.0.0",
      "model_name": "FloatProgressModel",
      "state": {
       "bar_style": "success",
       "layout": "IPY_MODEL_1abe1062e05d4382a5e7e49567965949",
       "max": 1000,
       "style": "IPY_MODEL_8ef1b0974ad748429e7290fd44338f01",
       "value": 1000
      }
     },
     "8793459176ba45feb7ecaf6537ab9355": {
      "model_module": "@jupyter-widgets/base",
      "model_module_version": "2.0.0",
      "model_name": "LayoutModel",
      "state": {
       "display": "inline-flex",
       "flex_flow": "row wrap",
       "width": "100%"
      }
     },
     "8ef1b0974ad748429e7290fd44338f01": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "2.0.0",
      "model_name": "ProgressStyleModel",
      "state": {
       "description_width": ""
      }
     },
     "8effff66a89f4b03a832dba64ac38315": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "2.0.0",
      "model_name": "HTMLStyleModel",
      "state": {
       "description_width": "",
       "font_size": null,
       "text_color": null
      }
     },
     "d8ed32c7c6f14c8da0b195252b452b2d": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "2.0.0",
      "model_name": "HTMLStyleModel",
      "state": {
       "description_width": "",
       "font_size": null,
       "text_color": null
      }
     },
     "f0f69820e69a4aeba774df7edf942573": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "2.0.0",
      "model_name": "HTMLModel",
      "state": {
       "layout": "IPY_MODEL_3e643e4ffcc64f2aa0e53e0433bde4aa",
       "style": "IPY_MODEL_d8ed32c7c6f14c8da0b195252b452b2d",
       "value": "proc00: 100%"
      }
     }
    },
    "version_major": 2,
    "version_minor": 0
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}