{ "cells": [ { "cell_type": "markdown", "metadata": { "id": "7FV_m2u2Hny8" }, "source": [ "##### Copyright 2022 The Cirq Developers" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "cellView": "form", "execution": { "iopub.execute_input": "2023-05-18T09:44:34.874768Z", "iopub.status.busy": "2023-05-18T09:44:34.874322Z", "iopub.status.idle": "2023-05-18T09:44:34.878191Z", "shell.execute_reply": "2023-05-18T09:44:34.877517Z" }, "id": "2XSh1-J7HrcU" }, "outputs": [], "source": [ "# @title Licensed under the Apache License, Version 2.0 (the \"License\");\n", "# you may not use this file except in compliance with the License.\n", "# You may obtain a copy of the License at\n", "#\n", "# https://www.apache.org/licenses/LICENSE-2.0\n", "#\n", "# Unless required by applicable law or agreed to in writing, software\n", "# distributed under the License is distributed on an \"AS IS\" BASIS,\n", "# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n", "# See the License for the specific language governing permissions and\n", "# limitations under the License." ] }, { "cell_type": "markdown", "metadata": { "id": "Me7G0WVTHwHg" }, "source": [ "# Cirq basics" ] }, { "cell_type": "markdown", "metadata": { "id": "aAFfDwu7IFyx" }, "source": [ "\n", " \n", " \n", " \n", " \n", "
\n", " View on QuantumAI\n", " \n", " Run in Google Colab\n", " \n", " View source on GitHub\n", " \n", " Download notebook\n", "
" ] }, { "cell_type": "markdown", "metadata": { "id": "67gUU740MSXd" }, "source": [ "This tutorial will teach the basics of how to use Cirq. It will walk through how to use qubits, gates, and operations to create and simulate your first quantum circuit using Cirq. It will briefly introduce devices, unitary matrices, decompositions, and transformers as well.\n", "\n", "This tutorial isn’t a quantum computing 101 tutorial: it assumes familiarity of quantum computing at about the level of the textbook “Quantum Computation and Quantum Information” by Nielsen and Chuang.\n", "\n", "For more in-depth examples of quantum algorithms and experiments, see [Experiments](../experiments)." ] }, { "cell_type": "markdown", "metadata": { "id": "1dOjJlgrNUuz" }, "source": [ "To begin, please follow the instructions for [installing Cirq](./install.md)." ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "execution": { "iopub.execute_input": "2023-05-18T09:44:34.881917Z", "iopub.status.busy": "2023-05-18T09:44:34.881378Z", "iopub.status.idle": "2023-05-18T09:44:56.741422Z", "shell.execute_reply": "2023-05-18T09:44:56.740482Z" }, "id": "bd9529db1c0b" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "installing cirq...\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "installed cirq.\n" ] } ], "source": [ "try:\n", " import cirq\n", "except ImportError:\n", " print(\"installing cirq...\")\n", " !pip install --quiet cirq\n", " print(\"installed cirq.\")\n", " import cirq\n", "\n", "import cirq_google" ] }, { "cell_type": "markdown", "metadata": { "id": "xr-MMoXgNsUQ" }, "source": [ "## Qubits\n", "\n", "The first part of creating a quantum circuit is to define a set of qubits (also known as a quantum register) to act on.\n", "\n", "Cirq has three main ways of defining qubits:\n", "\n", "* `cirq.NamedQubit`: used to label qubits by an abstract name.\n", "* `cirq.LineQubit`: qubits labelled by number in a linear array.\n", "* `cirq.GridQubit`: qubits labelled by two numbers in a rectangular lattice.\n", "\n", "Here are some examples of defining each type of qubit." ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "execution": { "iopub.execute_input": "2023-05-18T09:44:56.745885Z", "iopub.status.busy": "2023-05-18T09:44:56.745389Z", "iopub.status.idle": "2023-05-18T09:44:56.750377Z", "shell.execute_reply": "2023-05-18T09:44:56.749693Z" }, "id": "PsgSo-H0Os8X" }, "outputs": [], "source": [ "# Using named qubits can be useful for abstract algorithms\n", "# as well as algorithms not yet mapped onto hardware.\n", "q0 = cirq.NamedQubit('source')\n", "q1 = cirq.NamedQubit('target')\n", "\n", "# Line qubits can be created individually\n", "q3 = cirq.LineQubit(3)\n", "\n", "# Or created in a range\n", "# This will create LineQubit(0), LineQubit(1), LineQubit(2)\n", "q0, q1, q2 = cirq.LineQubit.range(3)\n", "\n", "# Grid Qubits can also be referenced individually\n", "q4_5 = cirq.GridQubit(4, 5)\n", "\n", "# Or created in bulk in a square\n", "# This will create 16 qubits from (0,0) to (3,3)\n", "qubits = cirq.GridQubit.square(4)" ] }, { "cell_type": "markdown", "metadata": { "id": "4zE6AutyQhQ6" }, "source": [ "There are also pre-packaged sets of qubits called [Devices](../hardware/devices.ipynb). These are qubits along with a set of rules for how they can be used. A `cirq.Device` can be used to ensure that two-qubit gates are only applied to qubits that are adjacent in the hardware, and other constraints. The following example will use the `cirq_google.Sycamore` device that comes with cirq. It is a diamond-shaped grid with 54 qubits that mimics early hardware released by Google." ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "execution": { "iopub.execute_input": "2023-05-18T09:44:56.753773Z", "iopub.status.busy": "2023-05-18T09:44:56.753179Z", "iopub.status.idle": "2023-05-18T09:44:56.762696Z", "shell.execute_reply": "2023-05-18T09:44:56.762007Z" }, "id": "B0Dwgu-lQLpq" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ " (0, 5)───(0, 6)\n", " │ │\n", " │ │\n", " (1, 4)───(1, 5)───(1, 6)───(1, 7)\n", " │ │ │ │\n", " │ │ │ │\n", " (2, 3)───(2, 4)───(2, 5)───(2, 6)───(2, 7)───(2, 8)\n", " │ │ │ │ │ │\n", " │ │ │ │ │ │\n", " (3, 2)───(3, 3)───(3, 4)───(3, 5)───(3, 6)───(3, 7)───(3, 8)───(3, 9)\n", " │ │ │ │ │ │ │ │\n", " │ │ │ │ │ │ │ │\n", " (4, 1)───(4, 2)───(4, 3)───(4, 4)───(4, 5)───(4, 6)───(4, 7)───(4, 8)───(4, 9)\n", " │ │ │ │ │ │ │ │\n", " │ │ │ │ │ │ │ │\n", "(5, 0)───(5, 1)───(5, 2)───(5, 3)───(5, 4)───(5, 5)───(5, 6)───(5, 7)───(5, 8)\n", " │ │ │ │ │ │ │\n", " │ │ │ │ │ │ │\n", " (6, 1)───(6, 2)───(6, 3)───(6, 4)───(6, 5)───(6, 6)───(6, 7)\n", " │ │ │ │ │\n", " │ │ │ │ │\n", " (7, 2)───(7, 3)───(7, 4)───(7, 5)───(7, 6)\n", " │ │ │\n", " │ │ │\n", " (8, 3)───(8, 4)───(8, 5)\n", " │\n", " │\n", " (9, 4)\n" ] } ], "source": [ "print(cirq_google.Sycamore)" ] }, { "cell_type": "markdown", "metadata": { "id": "j1QTjyxLSe5c" }, "source": [ "## Gates and operations\n", "\n", "The next step is to use the qubits to create operations that can be used in the circuit. Cirq has two concepts that are important to understand here:\n", "\n", "* A `Gate` is an effect that can be applied to a set of qubits. \n", "* An `Operation` is a gate applied to a set of qubits.\n", "\n", "For instance, `cirq.H` is the quantum [Hadamard](https://en.wikipedia.org/wiki/Quantum_logic_gate#Hadamard_(H)_gate) and is a `Gate` object. `cirq.H(cirq.LineQubit(1))` is an `Operation` object and is the Hadamard gate applied to a specific qubit (line qubit number 1).\n", "\n", "Many textbook gates are included within cirq. `cirq.X`, `cirq.Y`, and `cirq.Z` refer to the single-qubit Pauli gates. `cirq.CZ`, `cirq.CNOT`, `cirq.SWAP` are a few of the common two-qubit gates. `cirq.measure` is a macro to apply a `MeasurementGate` to a set of qubits. You can find more, as well as instructions on how to create your own custom gates, on the [Gates documentation](../build/gates.ipynb) page.\n", "\n", "Here are some examples of operations that can be performed on gates and operations:" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "execution": { "iopub.execute_input": "2023-05-18T09:44:56.766138Z", "iopub.status.busy": "2023-05-18T09:44:56.765625Z", "iopub.status.idle": "2023-05-18T09:44:56.770910Z", "shell.execute_reply": "2023-05-18T09:44:56.770206Z" }, "id": "wDW-yU-fesDl" }, "outputs": [], "source": [ "# Example gates\n", "cnot_gate = cirq.CNOT\n", "pauli_z = cirq.Z\n", "\n", "# Use exponentiation to get square root gates.\n", "sqrt_x_gate = cirq.X**0.5\n", "\n", "# Some gates can also take parameters\n", "sqrt_sqrt_y = cirq.YPowGate(exponent=0.25)\n", "\n", "# Create two qubits at once, in a line.\n", "q0, q1 = cirq.LineQubit.range(2)\n", "\n", "# Example operations\n", "z_op = cirq.Z(q0)\n", "not_op = cirq.CNOT(q0, q1)\n", "sqrt_iswap_op = cirq.SQRT_ISWAP(q0, q1)\n", "\n", "# You can also use the gates you specified earlier.\n", "cnot_op = cnot_gate(q0, q1)\n", "pauli_z_op = pauli_z(q0)\n", "sqrt_x_op = sqrt_x_gate(q0)\n", "sqrt_sqrt_y_op = sqrt_sqrt_y(q0)" ] }, { "cell_type": "markdown", "metadata": { "id": "BnBGLMVvWVkz" }, "source": [ "## Circuits and moments\n", "\n", "You are now ready to construct a quantum circuit. A `Circuit` is a collection of `Moment`s. A `Moment` is a collection of `Operation`s that all act during the same abstract time slice. Each `Operation` must be applied to a disjoint set of qubits compared to each of the other `Operation`s in the `Moment`. A `Moment` can be thought of as a vertical slice of a quantum circuit diagram.\n", "\n", "Circuits can be constructed in several different ways. By default, Cirq will attempt to slide your operation into the earliest possible `Moment` when you insert it. You can use the append function in two ways: \n", "\n", "By appending each operation one-by-one: " ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "execution": { "iopub.execute_input": "2023-05-18T09:44:56.774014Z", "iopub.status.busy": "2023-05-18T09:44:56.773756Z", "iopub.status.idle": "2023-05-18T09:44:56.779756Z", "shell.execute_reply": "2023-05-18T09:44:56.779044Z" }, "id": "HEuqEZcXkz3Q" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0: ───H───\n", "\n", "1: ───H───\n", "\n", "2: ───H───\n" ] } ], "source": [ "circuit = cirq.Circuit()\n", "qubits = cirq.LineQubit.range(3)\n", "circuit.append(cirq.H(qubits[0]))\n", "circuit.append(cirq.H(qubits[1]))\n", "circuit.append(cirq.H(qubits[2]))\n", "print(circuit)" ] }, { "cell_type": "markdown", "metadata": { "id": "ccd445492b19" }, "source": [ "Or by appending some iterable of operations. A preconstructed list works:" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "execution": { "iopub.execute_input": "2023-05-18T09:44:56.783287Z", "iopub.status.busy": "2023-05-18T09:44:56.782570Z", "iopub.status.idle": "2023-05-18T09:44:56.787780Z", "shell.execute_reply": "2023-05-18T09:44:56.787110Z" }, "id": "HEuqEZcXkz3Q" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0: ───H───\n", "\n", "1: ───H───\n", "\n", "2: ───H───\n" ] } ], "source": [ "circuit = cirq.Circuit()\n", "ops = [cirq.H(q) for q in cirq.LineQubit.range(3)]\n", "circuit.append(ops)\n", "print(circuit)" ] }, { "cell_type": "markdown", "metadata": { "id": "a0a5de4d8f01" }, "source": [ "A generator that yields operations also works. This syntax will be used often in documentation, and works both with the `cirq.Circuit()` initializer and the `cirq.Circuit.append()` function." ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "execution": { "iopub.execute_input": "2023-05-18T09:44:56.791040Z", "iopub.status.busy": "2023-05-18T09:44:56.790534Z", "iopub.status.idle": "2023-05-18T09:44:56.796772Z", "shell.execute_reply": "2023-05-18T09:44:56.796116Z" }, "id": "HEuqEZcXkz3Q" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0: ───H───\n", "\n", "1: ───H───\n", "\n", "2: ───H───\n", "0: ───H───\n", "\n", "1: ───H───\n", "\n", "2: ───H───\n" ] } ], "source": [ "# Append with generator\n", "circuit = cirq.Circuit()\n", "circuit.append(cirq.H(q) for q in cirq.LineQubit.range(3))\n", "print(circuit)\n", "# Initializer with generator\n", "print(cirq.Circuit(cirq.H(q) for q in cirq.LineQubit.range(3)))" ] }, { "cell_type": "markdown", "metadata": { "id": "a6bd4d3152d5" }, "source": [ "Note that all of the Hadamard gates are pushed as far left as possible, and put into the same Moment since none overlap.\n", "\n", "If your operations are applied to the same qubits, they will be put in sequential, insertion-ordered moments. In the following example, the two-qubit gates overlap, and are placed in consecutive moments." ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "execution": { "iopub.execute_input": "2023-05-18T09:44:56.800086Z", "iopub.status.busy": "2023-05-18T09:44:56.799558Z", "iopub.status.idle": "2023-05-18T09:44:56.805086Z", "shell.execute_reply": "2023-05-18T09:44:56.804382Z" }, "id": "Lbez4guQl31P" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0: ───×───────────\n", " │\n", "1: ───×───×───────\n", " │\n", "2: ───────×───×───\n", " │\n", "3: ───────────×───\n" ] } ], "source": [ "print(cirq.Circuit(cirq.SWAP(q, q + 1) for q in cirq.LineQubit.range(3)))" ] }, { "cell_type": "markdown", "metadata": { "id": "3FC9bdlXmShh" }, "source": [ "Sometimes, you may not want Cirq to automatically shift operations all the way to the left. To construct a circuit without doing this, you can create the circuit moment-by-moment or use a different `InsertStrategy`, explained more in the [Circuit documentation](../build/circuits.ipynb)." ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "execution": { "iopub.execute_input": "2023-05-18T09:44:56.808414Z", "iopub.status.busy": "2023-05-18T09:44:56.807893Z", "iopub.status.idle": "2023-05-18T09:44:56.812877Z", "shell.execute_reply": "2023-05-18T09:44:56.812186Z" }, "id": "4AEahodTnYiI" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0: ───H───────────\n", "\n", "1: ───────H───────\n", "\n", "2: ───────────H───\n" ] } ], "source": [ "# Creates each gate in a separate moment by passing an iterable of Moments instead of Operations.\n", "print(cirq.Circuit(cirq.Moment([cirq.H(q)]) for q in cirq.LineQubit.range(3)))" ] }, { "cell_type": "markdown", "metadata": { "id": "j406AKYsobpq" }, "source": [ "### Circuits and devices\n", "\n", "One important consideration when using real quantum devices is that there are often constraints on circuits that are able to be run on the hardware. `Device` objects specify these constraints and can be used to validate your circuit to make sure that it contains no illegal operations. For more information on what constraints `Device` objects can specify and how to use them, see the [Devices](../hardware/devices.ipynb) page.\n", "\n", "The following example demonstrates this with the Sycamore Device:" ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "execution": { "iopub.execute_input": "2023-05-18T09:44:56.816256Z", "iopub.status.busy": "2023-05-18T09:44:56.815725Z", "iopub.status.idle": "2023-05-18T09:44:56.821950Z", "shell.execute_reply": "2023-05-18T09:44:56.821234Z" }, "id": "9UV-dXJOpy8B" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Qubit pair is not valid on device: (cirq.GridQubit(5, 6), cirq.GridQubit(4, 5)).\n" ] } ], "source": [ "# Create some qubits.\n", "q0 = cirq.GridQubit(5, 6)\n", "q1 = cirq.GridQubit(5, 5)\n", "q2 = cirq.GridQubit(4, 5)\n", "\n", "# Create operations using the Sycamore gate, which is supported by the Sycamore device.\n", "# However, create operations for both adjacent and non-adjacent qubit pairs.\n", "adjacent_op = cirq_google.SYC(q0, q1)\n", "nonadjacent_op = cirq_google.SYC(q0, q2)\n", "\n", "# A working circuit for the Sycamore device raises no issues.\n", "working_circuit = cirq.Circuit()\n", "working_circuit.append(adjacent_op)\n", "valid = cirq_google.Sycamore.validate_circuit(working_circuit)\n", "\n", "# A circuit using invalid operations.\n", "bad_circuit = cirq.Circuit()\n", "bad_circuit.append(nonadjacent_op)\n", "try:\n", " cirq_google.Sycamore.validate_circuit(bad_circuit)\n", "except ValueError as e:\n", " print(e)" ] }, { "cell_type": "markdown", "metadata": { "id": "xZ68bWEjoMKt" }, "source": [ "## Simulation\n", "\n", "The results of the application of a quantum circuit can be calculated by a `Simulator`. Cirq comes bundled with a simulator that can calculate the results of circuits up to about a limit of 20 qubits. It can be initialized with `cirq.Simulator()`.\n", "\n", "There are two different approaches to using a simulator:\n", "\n", "* `simulate()`: When classically simulating a circuit, a simulator can directly access and view the resulting wave function. This is useful for debugging, learning, and understanding how circuits will function. \n", "* `run()`: When using actual quantum devices, we can only access the end result of a computation and must sample the results to get a distribution of results. Running the simulator as a sampler mimics this behavior and only returns bit strings as output.\n", "\n", "Next simulate a 2-qubit \"Bell State\":" ] }, { "cell_type": "code", "execution_count": 12, "metadata": { "execution": { "iopub.execute_input": "2023-05-18T09:44:56.825158Z", "iopub.status.busy": "2023-05-18T09:44:56.824730Z", "iopub.status.idle": "2023-05-18T09:44:56.836659Z", "shell.execute_reply": "2023-05-18T09:44:56.835957Z" }, "id": "AwC4SL6CHpXm" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Simulate the circuit:\n", "measurements: (no measurements)\n", "\n", "qubits: (cirq.LineQubit(0), cirq.LineQubit(1))\n", "output vector: 0.707|00⟩ + 0.707|11⟩\n", "\n", "phase:\n", "output vector: |⟩\n" ] } ], "source": [ "# Create a circuit to generate a Bell State:\n", "# 1/sqrt(2) * ( |00⟩ + |11⟩ )\n", "bell_circuit = cirq.Circuit()\n", "q0, q1 = cirq.LineQubit.range(2)\n", "bell_circuit.append(cirq.H(q0))\n", "bell_circuit.append(cirq.CNOT(q0, q1))\n", "\n", "# Initialize Simulator\n", "s = cirq.Simulator()\n", "\n", "print('Simulate the circuit:')\n", "results = s.simulate(bell_circuit)\n", "print(results)\n", "\n", "# For sampling, we need to add a measurement at the end\n", "bell_circuit.append(cirq.measure(q0, q1, key='result'))\n", "\n", "# Sample the circuit\n", "samples = s.run(bell_circuit, repetitions=1000)" ] }, { "cell_type": "markdown", "metadata": { "id": "a9a791442b74" }, "source": [ "### Virtual Machine Simulation\n", "\n", "Cirq also provides the [Quantum Virtual Machine](../simulate/quantum_virtual_machine.ipynb), which is a simulated virtual version of quantum hardware devices. It consists of two primary components: \n", "1. A [virtual Engine interface](../simulate/virtual_engine_interface.ipynb) that enables you to verify and run circuits with the same interface that quantum hardware would have.\n", "2. A set of [noise models](../noise/representing_noise.ipynb) that try to realistically replicate the noise present in actual Google quantum hardware devices.\n", "\n", "The QVM is intended to serve as a replacement for the Google quantum hardware, in two cases: \n", "1. Running your circuit on a QVM can give an approximation of how your circuit runs under the influence of hardware-like noise. This can be useful to help you reconfigure or change your circuit to be less impacted by noise when run on actual quantum hardware.\n", "2. During the preparation process before running your circuit(s) on hardware, you can seamlessly use a QVM to test that your workflow or pipeline works before eventually switching to the nearly-identical hardware interface.\n", "\n", "Read more about the QVM and its uses in the [Quantum Virtual Machine](../simulate/quantum_virtual_machine.ipynb) page, or just start using it with the [QVM Creation Template](../simulate/qvm_builder_code.ipynb) page." ] }, { "cell_type": "markdown", "metadata": { "id": "be79fa541b8a" }, "source": [ "## Visualizing Results\n", "\n", "When you use `run()` to get a sample distribution of measurements, you can directly graph the simulated samples as a histogram with `cirq.plot_state_histogram`." ] }, { "cell_type": "code", "execution_count": 13, "metadata": { "execution": { "iopub.execute_input": "2023-05-18T09:44:56.840134Z", "iopub.status.busy": "2023-05-18T09:44:56.839594Z", "iopub.status.idle": "2023-05-18T09:44:56.993234Z", "shell.execute_reply": "2023-05-18T09:44:56.992398Z" }, "id": "69dca44091f4" }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjsAAAHHCAYAAABZbpmkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA3fUlEQVR4nO3deVyVZf7/8fdBNgUOhAqkIrgVUm7hdszUFCUzR1OnzUk02xzQScrKqXGdwrFMW0j7zrdcSn412mhm5lpSKZlimGmaOhaWLG6AooLA/fujB+fbCTSOHDx4z+v5eNyPONd93df9ue8zA2/v+7rPsRiGYQgAAMCkPNxdAAAAQG0i7AAAAFMj7AAAAFMj7AAAAFMj7AAAAFMj7AAAAFMj7AAAAFMj7AAAAFMj7AAAAFMj7ACws1gsmjZtmrvLuOpw3oC6jbADXCGLFi2SxWKxL56enmratKlGjx6tn3/+2d3lVWnr1q2aNm2a8vPzq73Nhx9+qN69eyskJEQNGjRQy5Ytddddd2nt2rX2PkePHtW0adOUmZl52bWtWbOmVgLGtGnTZLFYdPz48SrXR0ZG6o477qjxflJTUzVv3rwajwPg9xF2gCtsxowZevvtt7VgwQINHDhQ77zzjnr37q3z58+7u7RKtm7dqunTp1c77Lz44ov6wx/+IIvFosmTJ2vu3LkaPny4Dhw4oHfffdfe7+jRo5o+fXqNw8706dMve3tXOnfunJ599lmntiHsAFeOp7sLAP7bDBw4UJ07d5YkPfjgg2rUqJH+8Y9/aNWqVbrrrrvcXN3lKy0t1cyZM9W/f3+tX7++0vq8vDw3VHVl+Pr6ursEpxUVFcnPz8/dZQBXBFd2ADe75ZZbJEmHDh1yaN+3b59GjBih4OBg+fr6qnPnzlq1apVDnwsXLmj69Olq06aNfH191bBhQ/Xs2VMbNmyw9+nTp4/69OlTab+jR49WZGTkReuaNm2aJk2aJElq0aKF/fbbDz/8UGX/48ePq7CwUDfffHOV60NCQiRJmzdvVpcuXSRJY8aMsY+7aNEiSdLnn3+uP/7xj2revLl8fHwUHh6uiRMn6ty5cw61p6SkSJLDrcEK5eXlmjdvnm644Qb5+voqNDRUjzzyiE6dOnXR462J387ZOX36tB577DFFRkbKx8dHISEh6t+/v3bu3Cnpl/fko48+0o8//miv/dfvRV5ensaOHavQ0FD5+vqqQ4cOWrx4caX9njhxQvfff7+sVquCgoIUHx+vXbt2OZxP6Zfz5e/vr0OHDun2229XQECARo4cKal65/vXY2RlZemOO+6Qv7+/mjZtan8fdu/erb59+8rPz08RERFKTU110dkFao4rO4CbVYSHa665xt62Z88e3XzzzWratKmefvpp+fn56V//+peGDh2q999/X3feeaekXwJJcnKyHnzwQXXt2lWFhYXasWOHdu7cqf79+9eormHDhun777/X//t//09z585Vo0aNJEmNGzeusn9ISIjq16+vDz/8UOPHj1dwcHCV/dq2basZM2ZoypQpevjhh+1hr0ePHpKkZcuW6ezZsxo3bpwaNmyor776Sq+++qp++uknLVu2TJL0yCOP6OjRo9qwYYPefvvtSvt45JFHtGjRIo0ZM0YTJkzQ4cOH9dprr+nrr7/Wli1b5OXl9bvHf/LkySrby8vLf3fbRx99VMuXL1diYqKio6N14sQJffHFF/ruu+9000036ZlnnlFBQYF++uknzZ07V5Lk7+8v6ZdbYn369NHBgweVmJioFi1aaNmyZRo9erTy8/P1l7/8xV7H4MGD9dVXX2ncuHGKiorSBx98oPj4+CprKi0tVVxcnHr27KkXX3xRDRo0kFS9812hrKxMAwcOVK9evTR79mwtXbpUiYmJ8vPz0zPPPKORI0dq2LBhWrBggUaNGiWbzaYWLVr87vkCap0B4IpYuHChIcnYuHGjcezYMePIkSPG8uXLjcaNGxs+Pj7GkSNH7H379etntGvXzjh//ry9rby83OjRo4fRpk0be1uHDh2MQYMGXXK/vXv3Nnr37l2pPT4+3oiIiHBok2RMnTrV/vqFF14wJBmHDx+u1jFOmTLFkGT4+fkZAwcONJ577jkjIyOjUr/t27cbkoyFCxdWWnf27NlKbcnJyYbFYjF+/PFHe1tCQoJR1a+wzz//3JBkLF261KF97dq1Vbb/1tSpUw1Jl1x+e85/e94CAwONhISES+5n0KBBlc6/YRjGvHnzDEnGO++8Y28rKSkxbDab4e/vbxQWFhqGYRjvv/++IcmYN2+evV9ZWZnRt2/fSuc2Pj7ekGQ8/fTTlfZX3fNdMcbzzz9vbzt16pRRv359w2KxGO+++669fd++fZXOCeBO3MYCrrDY2Fg1btxY4eHhGjFihPz8/LRq1So1a9ZM0i9XFD755BPdddddOn36tI4fP67jx4/rxIkTiouL04EDB+xPbwUFBWnPnj06cOCAOw/Jbvr06UpNTVWnTp20bt06PfPMM4qJidFNN92k7777rlpj1K9f3/5zUVGRjh8/rh49esgwDH399de/u/2yZcsUGBio/v3728/d8ePHFRMTI39/f3366afVquP999/Xhg0bKi2hoaG/u21QUJC2bdumo0ePVmtfv7ZmzRqFhYXp3nvvtbd5eXlpwoQJOnPmjNLS0iRJa9eulZeXlx566CF7Pw8PDyUkJFx07HHjxlVqc/Z8P/jgg/afg4KCdP3118vPz89hvtn111+voKAg/ec//6nmUQO1i9tYwBWWkpKi6667TgUFBXrrrbf02WefycfHx77+4MGDMgxDf/vb3/S3v/2tyjHy8vLUtGlTzZgxQ0OGDNF1112nG2+8Ubfddpvuv/9+tW/f/kodTiX33nuv7r33XhUWFmrbtm1atGiRUlNTNXjwYH377be/O5k3KytLU6ZM0apVqyrNsSkoKPjd/R84cEAFBQX2OUK/Vd2J0r169bLfuvu16kxGnj17tuLj4xUeHq6YmBjdfvvtGjVqlFq2bPm72/74449q06aNPDwc/y3atm1b+/qK/1577bX221EVWrduXeW4np6e9kD9a86cb19f30q3MQMDA9WsWTOHOVMV7bU1RwpwFmEHuMK6du1qfxpr6NCh6tmzp+677z7t379f/v7+9jkhTzzxhOLi4qoco+IPWq9evXTo0CF98MEHWr9+vf73f/9Xc+fO1YIFC+z/ArdYLDIMo9IYZWVltXF4dlarVf3791f//v3l5eWlxYsXa9u2berdu/dFtykrK1P//v118uRJPfXUU4qKipKfn59+/vlnjR49ulrzZcrLyxUSEqKlS5dWuf5ic45c6a677tItt9yiFStWaP369XrhhRf0j3/8Q//+9781cODAWt9/VXx8fCoFKGfPd7169aoc+2LtVf3vDnAHwg7gRvXq1VNycrJuvfVWvfbaa3r66aft//r38vJSbGzs744RHBysMWPGaMyYMTpz5ox69eqladOm2cPONddcU+XthIorBJfy23+tX67OnTtr8eLFys7OvuS4u3fv1vfff6/Fixdr1KhR9vZfP132e7W1atVKGzdu1M033+xwi+ZKu/baa/XnP/9Zf/7zn5WXl6ebbrpJzz33nD3sXKz+iIgIffPNNyovL3cIJ/v27bOvr/jvp59+qrNnzzpc3Tl48GC1a3TmfANXM+bsAG7Wp08fde3aVfPmzdP58+cVEhKiPn366I033rCHg187duyY/ecTJ044rPP391fr1q1VXFxsb2vVqpX27dvnsN2uXbu0ZcuW362t4nNYqvOhgmfPnlV6enqV6z7++GNJv8zluNS4FVcIfn1FwDAMvfzyy9Wu7a677lJZWZlmzpxZaZvS0lKnPg36cpSVlVW6/RMSEqImTZo4vC9+fn5V3pa7/fbblZOTo/fee8/eVlpaqldffVX+/v72K2NxcXG6cOGC/vnPf9r7lZeX2x8Frw5nzjdwNePKDlAHTJo0SX/84x+1aNEiPfroo0pJSVHPnj3Vrl07PfTQQ2rZsqVyc3OVnp6un376Sbt27ZIkRUdHq0+fPoqJiVFwcLB27Nhhf+S5wgMPPKCXXnpJcXFxGjt2rPLy8rRgwQLdcMMNKiwsvGRdMTExkqRnnnlG99xzj7y8vDR48OAqP4zu7Nmz6tGjh7p3767bbrtN4eHhys/P18qVK/X5559r6NCh6tSpk6RfAlhQUJAWLFiggIAA+fn5qVu3boqKilKrVq30xBNP6Oeff5bVatX7779f5dyPitomTJiguLg41atXT/fcc4969+6tRx55RMnJycrMzNSAAQPk5eWlAwcOaNmyZXr55Zc1YsSIy3ujquH06dNq1qyZRowYoQ4dOsjf318bN27U9u3bNWfOHIf633vvPSUlJalLly7y9/fX4MGD9fDDD+uNN97Q6NGjlZGRocjISC1fvlxbtmzRvHnzFBAQIOmXW6Bdu3bV448/roMHDyoqKkqrVq2yPzJfnatyzpxv4KrmtufAgP8yFY+eb9++vdK6srIyo1WrVkarVq2M0tJSwzAM49ChQ8aoUaOMsLAww8vLy2jatKlxxx13GMuXL7dv9/e//93o2rWrERQUZNSvX9+IiooynnvuOaOkpMRh/Hfeecdo2bKl4e3tbXTs2NFYt25dtR49NwzDmDlzptG0aVPDw8Pjko+hX7hwwfjnP/9pDB061IiIiDB8fHyMBg0aGJ06dTJeeOEFo7i42KH/Bx98YERHRxuenp4Oj0rv3bvXiI2NNfz9/Y1GjRoZDz30kLFr165Kj1OXlpYa48ePNxo3bmxYLJZKj6H/z//8jxETE2PUr1/fCAgIMNq1a2c8+eSTxtGjR6usv0LFo+fHjh2rcn1ERMQlHz0vLi42Jk2aZHTo0MEICAgw/Pz8jA4dOhivv/66wzZnzpwx7rvvPiMoKMiQ5PBe5ObmGmPGjDEaNWpkeHt7G+3atavyMf1jx44Z9913nxEQEGAEBgYao0ePNrZs2WJIcngUPD4+3vDz86vyeKp7vi82Ru/evY0bbrihWucJcBeLYTCDDADMYuXKlbrzzjv1xRdfXPTTrIH/NoQdALhKnTt3zmESdllZmQYMGKAdO3YoJyfHrRO0gbqEOTsAcJUaP368zp07J5vNpuLiYv373//W1q1b9fzzzxN0gF/hyg4AXKVSU1M1Z84cHTx4UOfPn1fr1q01btw4hwnqAAg7AADA5PicHQAAYGqEHQAAYGpMUNYvnzp69OhRBQQEuOzj8QEAQO0yDEOnT59WkyZNKn33268RdiQdPXpU4eHh7i4DAABchiNHjqhZs2YXXU/Ykewfv37kyBFZrVY3VwMAAKqjsLBQ4eHh9r/jF0PY0f99h4zVaiXsAABwlfm9KShMUAYAAKZG2AEAAKZG2AEAAKZG2AEAAKZG2AEAAKZG2AEAAKZG2AEAAKZG2AEAAKZG2AEAAKZG2AEAAKZG2AEAAKZG2AEAAKZG2AEAAKbm1rAzbdo0WSwWhyUqKsq+/vz580pISFDDhg3l7++v4cOHKzc312GMrKwsDRo0SA0aNFBISIgmTZqk0tLSK30oAACgjvJ0dwE33HCDNm7caH/t6fl/JU2cOFEfffSRli1bpsDAQCUmJmrYsGHasmWLJKmsrEyDBg1SWFiYtm7dquzsbI0aNUpeXl56/vnnr/ixAACAusftYcfT01NhYWGV2gsKCvTmm28qNTVVffv2lSQtXLhQbdu21Zdffqnu3btr/fr12rt3rzZu3KjQ0FB17NhRM2fO1FNPPaVp06bJ29v7Sh9OJZFPf+TuEv5r/TBrkLtLAADUAW6fs3PgwAE1adJELVu21MiRI5WVlSVJysjI0IULFxQbG2vvGxUVpebNmys9PV2SlJ6ernbt2ik0NNTeJy4uToWFhdqzZ89F91lcXKzCwkKHBQAAmJNbw063bt20aNEirV27VvPnz9fhw4d1yy236PTp08rJyZG3t7eCgoIctgkNDVVOTo4kKScnxyHoVKyvWHcxycnJCgwMtC/h4eGuPTAAAFBnuPU21sCBA+0/t2/fXt26dVNERIT+9a9/qX79+rW238mTJyspKcn+urCwkMADAIBJuX3Ozq8FBQXpuuuu08GDB9W/f3+VlJQoPz/f4epObm6ufY5PWFiYvvrqK4cxKp7WqmoeUAUfHx/5+Pi4/gAAAHUWcyjdx91zKN0+Z+fXzpw5o0OHDunaa69VTEyMvLy8tGnTJvv6/fv3KysrSzabTZJks9m0e/du5eXl2fts2LBBVqtV0dHRV7x+AABQ97j1ys4TTzyhwYMHKyIiQkePHtXUqVNVr1493XvvvQoMDNTYsWOVlJSk4OBgWa1WjR8/XjabTd27d5ckDRgwQNHR0br//vs1e/Zs5eTk6Nlnn1VCQgJXbgAAgCQ3h52ffvpJ9957r06cOKHGjRurZ8+e+vLLL9W4cWNJ0ty5c+Xh4aHhw4eruLhYcXFxev311+3b16tXT6tXr9a4ceNks9nk5+en+Ph4zZgxw12HBAAA6hi3hp133333kut9fX2VkpKilJSUi/aJiIjQmjVrXF0aAAAwiTo1ZwcAAMDVCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDUCDsAAMDU6kzYmTVrliwWix577DF72/nz55WQkKCGDRvK399fw4cPV25ursN2WVlZGjRokBo0aKCQkBBNmjRJpaWlV7h6AABQV9WJsLN9+3a98cYbat++vUP7xIkT9eGHH2rZsmVKS0vT0aNHNWzYMPv6srIyDRo0SCUlJdq6dasWL16sRYsWacqUKVf6EAAAQB3l9rBz5swZjRw5Uv/85z91zTXX2NsLCgr05ptv6qWXXlLfvn0VExOjhQsXauvWrfryyy8lSevXr9fevXv1zjvvqGPHjho4cKBmzpyplJQUlZSUuOuQAABAHeL2sJOQkKBBgwYpNjbWoT0jI0MXLlxwaI+KilLz5s2Vnp4uSUpPT1e7du0UGhpq7xMXF6fCwkLt2bPnovssLi5WYWGhwwIAAMzJ0507f/fdd7Vz505t37690rqcnBx5e3srKCjIoT00NFQ5OTn2Pr8OOhXrK9ZdTHJysqZPn17D6gEAwNXAbVd2jhw5or/85S9aunSpfH19r+i+J0+erIKCAvty5MiRK7p/AABw5bgt7GRkZCgvL0833XSTPD095enpqbS0NL3yyivy9PRUaGioSkpKlJ+f77Bdbm6uwsLCJElhYWGVns6qeF3Rpyo+Pj6yWq0OCwAAMCe3hZ1+/fpp9+7dyszMtC+dO3fWyJEj7T97eXlp06ZN9m3279+vrKws2Ww2SZLNZtPu3buVl5dn77NhwwZZrVZFR0df8WMCAAB1j9vm7AQEBOjGG290aPPz81PDhg3t7WPHjlVSUpKCg4NltVo1fvx42Ww2de/eXZI0YMAARUdH6/7779fs2bOVk5OjZ599VgkJCfLx8bnixwQAAOoet05Q/j1z586Vh4eHhg8fruLiYsXFxen111+3r69Xr55Wr16tcePGyWazyc/PT/Hx8ZoxY4YbqwYAAHVJnQo7mzdvdnjt6+urlJQUpaSkXHSbiIgIrVmzppYrAwAAVyu3f84OAABAbSLsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAU3M67CxZskTFxcWV2ktKSrRkyRKXFAUAAOAqToedMWPGqKCgoFL76dOnNWbMGJcUBQAA4CpOhx3DMGSxWCq1//TTTwoMDHRJUQAAAK7iWd2OnTp1ksVikcViUb9+/eTp+X+blpWV6fDhw7rttttqpUgAAIDLVe2wM3ToUElSZmam4uLi5O/vb1/n7e2tyMhIDR8+3OUFAgAA1ES1w87UqVMlSZGRkbr77rvl6+tba0UBAAC4SrXDToX4+HhJvzx9lZeXp/Lycof1zZs3d01lAAAALuB02Dlw4IAeeOABbd261aG9YuJyWVmZy4oDAACoKafDzujRo+Xp6anVq1fr2muvrfLJLAAAgLrC6bCTmZmpjIwMRUVF1UY9AAAALuX05+xER0fr+PHjLtn5/Pnz1b59e1mtVlmtVtlsNn388cf29efPn1dCQoIaNmwof39/DR8+XLm5uQ5jZGVladCgQWrQoIFCQkI0adIklZaWuqQ+AABw9XM67PzjH//Qk08+qc2bN+vEiRMqLCx0WJzRrFkzzZo1SxkZGdqxY4f69u2rIUOGaM+ePZKkiRMn6sMPP9SyZcuUlpamo0ePatiwYfbty8rKNGjQIJWUlGjr1q1avHixFi1apClTpjh7WAAAwKQshmEYzmzg4fFLPvrtXB1XTVAODg7WCy+8oBEjRqhx48ZKTU3ViBEjJEn79u1T27ZtlZ6eru7du+vjjz/WHXfcoaNHjyo0NFSStGDBAj311FM6duyYvL29q7XPwsJCBQYGqqCgQFartUb1/1bk0x+5dDxU3w+zBrm7BAB1CL+P3ae2fh9X9++303N2Pv300xoVdjFlZWVatmyZioqKZLPZlJGRoQsXLig2NtbeJyoqSs2bN7eHnfT0dLVr184edCQpLi5O48aN0549e9SpU6daqRUAAFw9nA47vXv3dmkBu3fvls1m0/nz5+Xv768VK1YoOjpamZmZ8vb2VlBQkEP/0NBQ5eTkSJJycnIcgk7F+op1F1NcXOzwze3O3n4DAABXD6fDzmeffXbJ9b169XJqvOuvv16ZmZkqKCjQ8uXLFR8fr7S0NGfLckpycrKmT59eq/sAAAB1g9Nhp0+fPpXafj1/x9k5O97e3mrdurUkKSYmRtu3b9fLL7+su+++WyUlJcrPz3e4upObm6uwsDBJUlhYmL766iuH8Sqe1qroU5XJkycrKSnJ/rqwsFDh4eFO1Q0AAK4OTj+NderUKYclLy9Pa9euVZcuXbR+/foaF1ReXq7i4mLFxMTIy8tLmzZtsq/bv3+/srKyZLPZJEk2m027d+9WXl6evc+GDRtktVoVHR190X34+PjYH3evWAAAgDk5fWUnMDCwUlv//v3l7e2tpKQkZWRkVHusyZMna+DAgWrevLlOnz6t1NRUbd68WevWrVNgYKDGjh2rpKQkBQcHy2q1avz48bLZbOrevbskacCAAYqOjtb999+v2bNnKycnR88++6wSEhLk4+Pj7KEBAAATcjrsXExoaKj279/v1DZ5eXkaNWqUsrOzFRgYqPbt22vdunXq37+/JGnu3Lny8PDQ8OHDVVxcrLi4OL3++uv27evVq6fVq1dr3Lhxstls8vPzU3x8vGbMmOGqwwIAAFc5p8PON9984/DaMAxlZ2dr1qxZ6tixo1Njvfnmm5dc7+vrq5SUFKWkpFy0T0REhNasWePUfgEAwH8Pp8NOx44dZbFY9NvPIuzevbveeustlxUGAADgCk6HncOHDzu89vDwUOPGjeXr6+uyogAAAFzF6bATERFRG3UAAADUCqcfPZektLQ0DR48WK1bt1br1q31hz/8QZ9//rmrawMAAKgxp8POO++8o9jYWDVo0EATJkzQhAkTVL9+ffXr10+pqam1USMAAMBlc/o21nPPPafZs2dr4sSJ9rYJEybopZde0syZM3Xfffe5tEAAAICacPrKzn/+8x8NHjy4Uvsf/vCHSpOXAQAA3M3psBMeHu7wFQ4VNm7cyPdLAQCAOsfp21iPP/64JkyYoMzMTPXo0UOStGXLFi1atEgvv/yyywsEAACoCafDzrhx4xQWFqY5c+boX//6lySpbdu2eu+99zRkyBCXFwgAAFATl/XdWHfeeafuvPNOV9cCAADgck7P2dm+fbu2bdtWqX3btm3asWOHS4oCAABwFafDTkJCgo4cOVKp/eeff1ZCQoJLigIAAHAVp8PO3r17ddNNN1Vq79Spk/bu3euSogAAAFzF6bDj4+Oj3NzcSu3Z2dny9LysKUAAAAC1xumwM2DAAE2ePFkFBQX2tvz8fP31r39V//79XVocAABATTl9KebFF19Ur169FBERoU6dOkmSMjMzFRoaqrffftvlBQIAANSE02GnadOm+uabb7R06VLt2rVL9evX15gxY3TvvffKy8urNmoEAAC4bJc1ycbPz08PP/ywq2sBAABwOafn7AAAAFxNCDsAAMDUCDsAAMDUCDsAAMDUnA47LVu21IkTJyq15+fnq2XLli4pCgAAwFWcDjs//PCDysrKKrUXFxfr559/dklRAAAArlLtR89XrVpl/3ndunUKDAy0vy4rK9OmTZsUGRnp0uIAAABqqtphZ+jQoZIki8Wi+Ph4h3VeXl6KjIzUnDlzXFocAABATVU77JSXl0uSWrRooe3bt6tRo0a1VhQAAICrOP0JyocPH66NOgAAAGpFtcLOK6+8Uu0BJ0yYcNnFAAAAuFq1ws7cuXOrNZjFYiHsAACAOqVaYYdbVwAA4GrFJygDAABTc3qC8gMPPHDJ9W+99dZlFwMAAOBqToedU6dOOby+cOGCvv32W+Xn56tv374uKwwAAMAVnA47K1asqNRWXl6ucePGqVWrVi4pCgAAwFVcMmfHw8NDSUlJ1X5qCwAA4Epx2QTlQ4cOqbS01FXDAQAAuITTt7GSkpIcXhuGoezsbH300UeVvjMLAADA3ZwOO19//bXDaw8PDzVu3Fhz5sz53Se1AAAArjSnw86nn35aG3UAAADUCqfn7Jw7d05nz561v/7xxx81b948rV+/3qWFAQAAuILTYWfIkCFasmSJJCk/P19du3bVnDlzNGTIEM2fP9/lBQIAANSE02Fn586duuWWWyRJy5cvV1hYmH788UctWbLEqW9HBwAAuBKcDjtnz55VQECAJGn9+vUaNmyYPDw81L17d/34448uLxAAAKAmnA47rVu31sqVK3XkyBGtW7dOAwYMkCTl5eXJarW6vEAAAICacDrsTJkyRU888YQiIyPVtWtX2Ww2Sb9c5enUqZPLCwQAAKgJpx89HzFihHr27Kns7Gx16NDB3t6vXz/deeedLi0OAACgpi7r6yLCwsIUEBCgDRs26Ny5c5KkLl26KCoqyqXFAQAA1JTTYefEiRPq16+frrvuOt1+++3Kzs6WJI0dO1aPP/64ywsEAACoCafDzsSJE+Xl5aWsrCw1aNDA3n733Xdr7dq1Li0OAACgppyes7N+/XqtW7dOzZo1c2hv06YNj54DAIA6x+krO0VFRQ5XdCqcPHlSPj4+LikKAADAVZwOO7fccov96yIkyWKxqLy8XLNnz9att97q0uIAAABqyunbWLNnz1a/fv20Y8cOlZSU6Mknn9SePXt08uRJbdmypTZqBAAAuGxOX9m58cYb9f3336tnz54aMmSIioqKNGzYMH399ddq1apVbdQIAABw2Zy6snPhwgXddtttWrBggZ555pnaqgkAAMBlnLqy4+XlpW+++aa2agEAAHA5p29j/elPf9Kbb75ZG7UAAAC4nNNhp7S0VPPnz1fnzp31yCOPKCkpyWFxRnJysrp06aKAgACFhIRo6NCh2r9/v0Of8+fPKyEhQQ0bNpS/v7+GDx+u3Nxchz5ZWVkaNGiQGjRooJCQEE2aNEmlpaXOHhoAADAhp5/G+vbbb3XTTTdJkr7//nuHdRaLxamx0tLSlJCQoC5duqi0tFR//etfNWDAAO3du1d+fn6SfvnE5o8++kjLli1TYGCgEhMTNWzYMPuTX2VlZRo0aJDCwsK0detWZWdna9SoUfLy8tLzzz/v7OEBAACTsRiGYbi7iArHjh1TSEiI0tLS1KtXLxUUFKhx48ZKTU3ViBEjJEn79u1T27ZtlZ6eru7du+vjjz/WHXfcoaNHjyo0NFSStGDBAj311FM6duyYvL29f3e/hYWFCgwMVEFBgaxWq0uPKfLpj1w6Hqrvh1mD3F0CgDqE38fuU1u/j6v79/uyvvW8thQUFEiSgoODJUkZGRm6cOGCYmNj7X2ioqLUvHlzpaenS5LS09PVrl07e9CRpLi4OBUWFmrPnj1XsHoAAFAXOX0bq7aUl5frscce080336wbb7xRkpSTkyNvb28FBQU59A0NDVVOTo69z6+DTsX6inVVKS4uVnFxsf11YWGhqw4DAADUMXXmyk5CQoK+/fZbvfvuu7W+r+TkZAUGBtqX8PDwWt8nAABwjzoRdhITE7V69Wp9+umnDt+mHhYWppKSEuXn5zv0z83NVVhYmL3Pb5/Oqnhd0ee3Jk+erIKCAvty5MgRFx4NAACoS9wadgzDUGJiolasWKFPPvlELVq0cFgfExMjLy8vbdq0yd62f/9+ZWVlyWazSZJsNpt2796tvLw8e58NGzbIarUqOjq6yv36+PjIarU6LAAAwJzcOmcnISFBqamp+uCDDxQQEGCfYxMYGKj69esrMDBQY8eOVVJSkoKDg2W1WjV+/HjZbDZ1795dkjRgwABFR0fr/vvv1+zZs5WTk6Nnn31WCQkJ8vHxcefhAQCAOsCtYWf+/PmSpD59+ji0L1y4UKNHj5YkzZ07Vx4eHho+fLiKi4sVFxen119/3d63Xr16Wr16tcaNGyebzSY/Pz/Fx8drxowZV+owAABAHebWsFOdj/jx9fVVSkqKUlJSLtonIiJCa9ascWVpAADAJOrEBGUAAIDaQtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACm5taw89lnn2nw4MFq0qSJLBaLVq5c6bDeMAxNmTJF1157rerXr6/Y2FgdOHDAoc/Jkyc1cuRIWa1WBQUFaezYsTpz5swVPAoAAFCXuTXsFBUVqUOHDkpJSaly/ezZs/XKK69owYIF2rZtm/z8/BQXF6fz58/b+4wcOVJ79uzRhg0btHr1an322Wd6+OGHr9QhAACAOs7TnTsfOHCgBg4cWOU6wzA0b948PfvssxoyZIgkacmSJQoNDdXKlSt1zz336LvvvtPatWu1fft2de7cWZL06quv6vbbb9eLL76oJk2aXLFjAQAAdVOdnbNz+PBh5eTkKDY21t4WGBiobt26KT09XZKUnp6uoKAge9CRpNjYWHl4eGjbtm1XvGYAAFD3uPXKzqXk5ORIkkJDQx3aQ0ND7etycnIUEhLisN7T01PBwcH2PlUpLi5WcXGx/XVhYaGrygYAAHVMnb2yU5uSk5MVGBhoX8LDw91dEgAAqCV1NuyEhYVJknJzcx3ac3Nz7evCwsKUl5fnsL60tFQnT56096nK5MmTVVBQYF+OHDni4uoBAEBdUWfDTosWLRQWFqZNmzbZ2woLC7Vt2zbZbDZJks1mU35+vjIyMux9PvnkE5WXl6tbt24XHdvHx0dWq9VhAQAA5uTWOTtnzpzRwYMH7a8PHz6szMxMBQcHq3nz5nrsscf097//XW3atFGLFi30t7/9TU2aNNHQoUMlSW3bttVtt92mhx56SAsWLNCFCxeUmJioe+65hyexAACAJDeHnR07dujWW2+1v05KSpIkxcfHa9GiRXryySdVVFSkhx9+WPn5+erZs6fWrl0rX19f+zZLly5VYmKi+vXrJw8PDw0fPlyvvPLKFT8WAABQN1kMwzDcXYS7FRYWKjAwUAUFBS6/pRX59EcuHQ/V98OsQe4uAUAdwu9j96mt38fV/ftdZ+fsAAAAuAJhBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmBphBwAAmJppwk5KSooiIyPl6+urbt266auvvnJ3SQAAoA4wRdh57733lJSUpKlTp2rnzp3q0KGD4uLilJeX5+7SAACAm5ki7Lz00kt66KGHNGbMGEVHR2vBggVq0KCB3nrrLXeXBgAA3OyqDzslJSXKyMhQbGysvc3Dw0OxsbFKT093Y2UAAKAu8HR3ATV1/PhxlZWVKTQ01KE9NDRU+/btq3Kb4uJiFRcX218XFBRIkgoLC11eX3nxWZePieqpjfcTwNWL38fuU1u/jyvGNQzjkv2u+rBzOZKTkzV9+vRK7eHh4W6oBrUlcJ67KwAASLX/+/j06dMKDAy86PqrPuw0atRI9erVU25urkN7bm6uwsLCqtxm8uTJSkpKsr8uLy/XyZMn1bBhQ1ksllqt92pSWFio8PBwHTlyRFar1d3lwEV4X82L99a8eG+rZhiGTp8+rSZNmlyy31Ufdry9vRUTE6NNmzZp6NChkn4JL5s2bVJiYmKV2/j4+MjHx8ehLSgoqJYrvXpZrVb+z2VCvK/mxXtrXry3lV3qik6Fqz7sSFJSUpLi4+PVuXNnde3aVfPmzVNRUZHGjBnj7tIAAICbmSLs3H333Tp27JimTJminJwcdezYUWvXrq00aRkAAPz3MUXYkaTExMSL3rbC5fHx8dHUqVMr3fLD1Y331bx4b82L97ZmLMbvPa8FAABwFbvqP1QQAADgUgg7AADA1Ag7AADA1Ag7AADA1Ag7qFJKSooiIyPl6+urbt266auvvnJ3Saihzz77TIMHD1aTJk1ksVi0cuVKd5cEF0lOTlaXLl0UEBCgkJAQDR06VPv373d3Waih+fPnq3379vYPErTZbPr444/dXdZVibCDSt577z0lJSVp6tSp2rlzpzp06KC4uDjl5eW5uzTUQFFRkTp06KCUlBR3lwIXS0tLU0JCgr788ktt2LBBFy5c0IABA1RUVOTu0lADzZo106xZs5SRkaEdO3aob9++GjJkiPbs2ePu0q46PHqOSrp166YuXbrotddek/TL12+Eh4dr/Pjxevrpp91cHVzBYrFoxYoV9q9YgbkcO3ZMISEhSktLU69evdxdDlwoODhYL7zwgsaOHevuUq4qXNmBg5KSEmVkZCg2Ntbe5uHhodjYWKWnp7uxMgDVVVBQIOmXP4wwh7KyMr377rsqKiqSzWZzdzlXHdN8gjJc4/jx4yorK6v0VRuhoaHat2+fm6oCUF3l5eV67LHHdPPNN+vGG290dzmood27d8tms+n8+fPy9/fXihUrFB0d7e6yrjqEHQAwkYSEBH377bf64osv3F0KXOD6669XZmamCgoKtHz5csXHxystLY3A4yTCDhw0atRI9erVU25urkN7bm6uwsLC3FQVgOpITEzU6tWr9dlnn6lZs2buLgcu4O3trdatW0uSYmJitH37dr388st644033FzZ1YU5O3Dg7e2tmJgYbdq0yd5WXl6uTZs2cZ8YqKMMw1BiYqJWrFihTz75RC1atHB3Sagl5eXlKi4udncZVx2u7KCSpKQkxcfHq3PnzuratavmzZunoqIijRkzxt2loQbOnDmjgwcP2l8fPnxYmZmZCg4OVvPmzd1YGWoqISFBqamp+uCDDxQQEKCcnBxJUmBgoOrXr+/m6nC5Jk+erIEDB6p58+Y6ffq0UlNTtXnzZq1bt87dpV11ePQcVXrttdf0wgsvKCcnRx07dtQrr7yibt26ubss1MDmzZt16623VmqPj4/XokWLrnxBcBmLxVJl+8KFCzV69OgrWwxcZuzYsdq0aZOys7MVGBio9u3b66mnnlL//v3dXdpVh7ADAABMjTk7AADA1Ag7AADA1Ag7AADA1Ag7AADA1Ag7AADA1Ag7AADA1Ag7AADA1Ag7AOqcH374QRaLRZmZmRfts3nzZlksFuXn51+xugBcnQg7AK5KPXr0sH+yrCQtWrRIQUFBNRrzcscgeAF1G9+NBeCq5O3trbCwMHeXAeAqwJUdAC5VVFSkUaNGyd/fX9dee63mzJmjPn366LHHHrP3sVgsWrlypcN2QUFBlb6ja9++ferRo4d8fX114403Ki0tzb7u11dTNm/erDFjxqigoEAWi0UWi0XTpk2rsr5du3bp1ltvVUBAgKxWq2JiYrRjx45LjvH222+rc+fOCggIUFhYmO677z7l5eVJ+uWWW8V3jl1zzTWyWCz276MqLy9XcnKyWrRoofr166tDhw5avnz5ZZ9bAJeHsAPApSZNmqS0tDR98MEHWr9+vTZv3qydO3de9liPP/64vv76a9lsNg0ePFgnTpyo1K9Hjx6aN2+erFarsrOzlZ2drSeeeKLKMUeOHKlmzZpp+/btysjI0NNPPy0vL69LjnHhwgXNnDlTu3bt0sqVK/XDDz/YA014eLjef/99SdL+/fuVnZ2tl19+WZKUnJysJUuWaMGCBdqzZ48mTpyoP/3pTw6hDUDt4zYWAJc5c+aM3nzzTb3zzjvq16+fJGnx4sVq1qzZZY2XmJio4cOHS5Lmz5+vtWvX6s0339STTz7p0M/b21uBgYGyWCy/e2srKytLkyZNUlRUlCSpTZs29nUXG+OBBx6w/9yyZUu98sor6tKli86cOSN/f38FBwdLkkJCQuxzfoqLi/X8889r48aNstls9m2/+OILvfHGG+rdu/dlnBEAl4MrOwBc5tChQyopKVG3bt3sbcHBwbr++usva7yKkCBJnp6e6ty5s7777rsa1ZiUlKQHH3xQsbGxmjVrlg4dOvS722RkZGjw4MFq3ry5AgIC7EElKyvrotscPHhQZ8+eVf/+/eXv729flixZUq19AnAdwg6AK85iscgwDIe2CxcuXJF9T5s2TXv27NGgQYP0ySefKDo6WitWrLho/6KiIsXFxclqtWrp0qXavn27vX9JSclFtztz5owk6aOPPlJmZqZ92bt3L/N2gCuMsAPAZVq1aiUvLy9t27bN3nbq1Cl9//33Dv0aN26s7Oxs++sDBw7o7Nmzlcb78ssv7T+XlpYqIyNDbdu2rXLf3t7eKisrq1ad1113nSZOnKj169dr2LBhWrhw4UXH2Ldvn06cOKFZs2bplltuUVRUlH1y8q/3Lclh2+joaPn4+CgrK0utW7d2WMLDw6tVJwDXYM4OAJfx9/fX2LFjNWnSJDVs2FAhISF65pln5OHh+O+qvn376rXXXpPNZlNZWZmeeuopeXl5VRovJSVFbdq0Udu2bTV37lydOnXKYf7Mr0VGRurMmTPatGmTOnTooAYNGqhBgwYOfc6dO6dJkyZpxIgRatGihX766Sdt377dPi+oqjGaN28ub29vvfrqq3r00Uf17bffaubMmQ7jRkREyGKxaPXq1br99ttVv359BQQE6IknntDEiRNVXl6unj17qqCgQFu2bJHValV8fHxNTjUAZxgA4EKnT582/vSnPxkNGjQwQkNDjdmzZxu9e/c2/vKXv9j7/Pzzz8aAAQMMPz8/o02bNsaaNWuMwMBAY+HChYZhGMbhw4cNSUZqaqrRtWtXw9vb24iOjjY++eQT+xiffvqpIck4deqUve3RRx81GjZsaEgypk6dWqm24uJi45577jHCw8MNb29vo0mTJkZiYqJx7ty5S46RmppqREZGGj4+PobNZjNWrVplSDK+/vpr+3YzZswwwsLCDIvFYsTHxxuGYRjl5eXGvHnzjOuvv97w8vIyGjdubMTFxRlpaWk1Pc0AnGAxjN/cOAcAF+vTp486duyoefPmubsUAP+FmLMDAABMjbADAABMjdtYAADA1LiyAwAATI2wAwAATI2wAwAATI2wAwAATI2wAwAATI2wAwAATI2wAwAATI2wAwAATI2wAwAATO3/Azdiu7V0JowFAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt\n", "\n", "cirq.plot_state_histogram(samples, plt.subplot())\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": { "id": "b4f20cffe4ae" }, "source": [ "However, this histogram has some empty qubit states, which may become problematic if you work with more qubits. To graph sparse sampled data, first get the `Counts` from your results with its `histogram()` function, and pass that to `cirq.plot_state_histogram`. By collecting the results into counts, all the qubit states that were never seen are ignored." ] }, { "cell_type": "code", "execution_count": 14, "metadata": { "execution": { "iopub.execute_input": "2023-05-18T09:44:56.996669Z", "iopub.status.busy": "2023-05-18T09:44:56.996055Z", "iopub.status.idle": "2023-05-18T09:44:57.125942Z", "shell.execute_reply": "2023-05-18T09:44:57.125170Z" }, "id": "07ec16b02594" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Counter({0: 503, 3: 497})\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjsAAAHHCAYAAABZbpmkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA2R0lEQVR4nO3deVyVZf7/8fdBNgUOhAqkIrgVUm7hdszUFCQjR1OnzUk02xzQScrKqXGdwrFMWkj7zrdcSn412mRm5lpao2SKYaZp6lhYsrgBriBw//7owfl2ApWjBw/evZ6Px/2Ic93Xfd2f+6ajb+/7us+xGIZhCAAAwKQ83F0AAABAbSLsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsAAAAUyPsALCzWCyaMmWKu8u46nDegLqNsANcIfPnz5fFYrEvnp6eatq0qUaOHKmff/7Z3eVVa9OmTZoyZYoKCwtrvM1HH32k3r17KyQkRA0aNFDLli111113aeXKlfY+hw4d0pQpU5SdnX3Jta1YsaJWAsaUKVNksVh05MiRatdHRkbqjjvuuOz9ZGRkKC0t7bLHAXBxhB3gCps2bZrefvttzZ07VwMGDNA777yj3r176+zZs+4urYpNmzZp6tSpNQ47L774ov7whz/IYrFo4sSJmj17toYOHaq9e/fq3Xfftfc7dOiQpk6detlhZ+rUqZe8vSudOXNGzz77rFPbEHaAK8fT3QUAvzcDBgxQ586dJUkPPvigGjVqpH/84x9atmyZ7rrrLjdXd+nKyso0ffp0xcXFafXq1VXWFxQUuKGqK8PX19fdJTjt1KlT8vPzc3cZwBXBlR3AzW655RZJ0v79+x3ad+/erWHDhik4OFi+vr7q3Lmzli1b5tDn3Llzmjp1qtq0aSNfX181bNhQPXv21Jo1a+x9+vTpoz59+lTZ78iRIxUZGXneuqZMmaIJEyZIklq0aGG//fbDDz9U2//IkSMqLi7WzTffXO36kJAQSdL69evVpUsXSdKoUaPs486fP1+S9MUXX+iPf/yjmjdvLh8fH4WHh2v8+PE6c+aMQ+3p6emS5HBrsFJFRYXS0tJ0ww03yNfXV6GhoXrkkUd0/Pjx8x7v5fjtnJ0TJ07oscceU2RkpHx8fBQSEqK4uDht27ZN0i+/k48//lg//vijvfZf/y4KCgo0evRohYaGytfXVx06dNCCBQuq7Pfo0aO6//77ZbVaFRQUpMTERG3fvt3hfEq/nC9/f3/t379ft99+uwICAjR8+HBJNTvfvx4jJydHd9xxh/z9/dW0aVP772HHjh3q27ev/Pz8FBERoYyMDBedXeDycWUHcLPK8HDNNdfY23bu3Kmbb75ZTZs21dNPPy0/Pz/961//0uDBg/X+++/rzjvvlPRLIElNTdWDDz6orl27qri4WFu3btW2bdsUFxd3WXUNGTJE33//vf7f//t/mj17tho1aiRJaty4cbX9Q0JCVL9+fX300UcaO3asgoODq+3Xtm1bTZs2TZMmTdLDDz9sD3s9evSQJC1evFinT5/WmDFj1LBhQ3311Vd69dVX9dNPP2nx4sWSpEceeUSHDh3SmjVr9Pbbb1fZxyOPPKL58+dr1KhRGjdunA4cOKDXXntNX3/9tTZu3CgvL6+LHv+xY8eqba+oqLjoto8++qiWLFmi5ORkRUdH6+jRo/rPf/6j7777TjfddJOeeeYZFRUV6aefftLs2bMlSf7+/pJ+uSXWp08f7du3T8nJyWrRooUWL16skSNHqrCwUH/5y1/sdQwcOFBfffWVxowZo6ioKH344YdKTEystqaysjLFx8erZ8+eevHFF9WgQQNJNTvflcrLyzVgwAD16tVLM2fO1KJFi5ScnCw/Pz8988wzGj58uIYMGaK5c+dqxIgRstlsatGixUXPF1DrDABXxLx58wxJxtq1a43Dhw8bBw8eNJYsWWI0btzY8PHxMQ4ePGjv269fP6Ndu3bG2bNn7W0VFRVGjx49jDZt2tjbOnToYCQkJFxwv7179zZ69+5dpT0xMdGIiIhwaJNkTJ482f76hRdeMCQZBw4cqNExTpo0yZBk+Pn5GQMGDDCee+45Iysrq0q/LVu2GJKMefPmVVl3+vTpKm2pqamGxWIxfvzxR3tbUlKSUd0fYV988YUhyVi0aJFD+8qVK6tt/63Jkycbki64/Pac//a8BQYGGklJSRfcT0JCQpXzbxiGkZaWZkgy3nnnHXtbaWmpYbPZDH9/f6O4uNgwDMN4//33DUlGWlqavV95ebnRt2/fKuc2MTHRkGQ8/fTTVfZX0/NdOcbzzz9vbzt+/LhRv359w2KxGO+++669fffu3VXOCeBO3MYCrrDY2Fg1btxY4eHhGjZsmPz8/LRs2TI1a9ZM0i9XFD799FPdddddOnHihI4cOaIjR47o6NGjio+P1969e+1PbwUFBWnnzp3au3evOw/JburUqcrIyFCnTp20atUqPfPMM4qJidFNN92k7777rkZj1K9f3/7zqVOndOTIEfXo0UOGYejrr7++6PaLFy9WYGCg4uLi7OfuyJEjiomJkb+/vz777LMa1fH+++9rzZo1VZbQ0NCLbhsUFKTNmzfr0KFDNdrXr61YsUJhYWG699577W1eXl4aN26cTp48qQ0bNkiSVq5cKS8vLz300EP2fh4eHkpKSjrv2GPGjKnS5uz5fvDBB+0/BwUF6frrr5efn5/DfLPrr79eQUFB+u9//1vDowZqF7exgCssPT1d1113nYqKivTWW2/p888/l4+Pj339vn37ZBiG/va3v+lvf/tbtWMUFBSoadOmmjZtmgYNGqTrrrtON954o2677Tbdf//9at++/ZU6nCruvfde3XvvvSouLtbmzZs1f/58ZWRkaODAgfr2228vOpk3JydHkyZN0rJly6rMsSkqKrro/vfu3auioiL7HKHfqulE6V69etlv3f1aTSYjz5w5U4mJiQoPD1dMTIxuv/12jRgxQi1btrzotj/++KPatGkjDw/Hf4u2bdvWvr7yv9dee639dlSl1q1bVzuup6enPVD/mjPn29fXt8ptzMDAQDVr1sxhzlRle23NkQKcRdgBrrCuXbvan8YaPHiwevbsqfvuu0979uyRv7+/fU7IE088ofj4+GrHqPwLrVevXtq/f78+/PBDrV69Wv/7v/+r2bNna+7cufZ/gVssFhmGUWWM8vLy2jg8O6vVqri4OMXFxcnLy0sLFizQ5s2b1bt37/NuU15erri4OB07dkxPPfWUoqKi5Ofnp59//lkjR46s0XyZiooKhYSEaNGiRdWuP9+cI1e66667dMstt+iDDz7Q6tWr9cILL+gf//iH/v3vf2vAgAG1vv/q+Pj4VAlQzp7vevXqVTv2+dqr+/8OcAfCDuBG9erVU2pqqm699Va99tprevrpp+3/+vfy8lJsbOxFxwgODtaoUaM0atQonTx5Ur169dKUKVPsYeeaa66p9nZC5RWCC/ntv9YvVefOnbVgwQLl5uZecNwdO3bo+++/14IFCzRixAh7+6+fLrtYba1atdLatWt18803O9yiudKuvfZa/fnPf9af//xnFRQU6KabbtJzzz1nDzvnqz8iIkLffPONKioqHMLJ7t277esr//vZZ5/p9OnTDld39u3bV+ManTnfwNWMOTuAm/Xp00ddu3ZVWlqazp49q5CQEPXp00dvvPGGPRz82uHDh+0/Hz161GGdv7+/WrdurZKSEntbq1attHv3bofttm/fro0bN160tsrPYanJhwqePn1amZmZ1a775JNPJP0yl+NC41ZeIfj1FQHDMPTyyy/XuLa77rpL5eXlmj59epVtysrKnPo06EtRXl5e5fZPSEiImjRp4vB78fPzq/a23O233668vDy999579raysjK9+uqr8vf3t18Zi4+P17lz5/TPf/7T3q+iosL+KHhNOHO+gasZV3aAOmDChAn64x//qPnz5+vRRx9Venq6evbsqXbt2umhhx5Sy5YtlZ+fr8zMTP3000/avn27JCk6Olp9+vRRTEyMgoODtXXrVvsjz5UeeOABvfTSS4qPj9fo0aNVUFCguXPn6oYbblBxcfEF64qJiZEkPfPMM7rnnnvk5eWlgQMHVvthdKdPn1aPHj3UvXt33XbbbQoPD1dhYaGWLl2qL774QoMHD1anTp0k/RLAgoKCNHfuXAUEBMjPz0/dunVTVFSUWrVqpSeeeEI///yzrFar3n///WrnflTWNm7cOMXHx6tevXq655571Lt3bz3yyCNKTU1Vdna2+vfvLy8vL+3du1eLFy/Wyy+/rGHDhl3aL6oGTpw4oWbNmmnYsGHq0KGD/P39tXbtWm3ZskWzZs1yqP+9995TSkqKunTpIn9/fw0cOFAPP/yw3njjDY0cOVJZWVmKjIzUkiVLtHHjRqWlpSkgIEDSL7dAu3btqscff1z79u1TVFSUli1bZn9kviZX5Zw538BVzW3PgQG/M5WPnm/ZsqXKuvLycqNVq1ZGq1atjLKyMsMwDGP//v3GiBEjjLCwMMPLy8to2rSpcccddxhLliyxb/f3v//d6Nq1qxEUFGTUr1/fiIqKMp577jmjtLTUYfx33nnHaNmypeHt7W107NjRWLVqVY0ePTcMw5g+fbrRtGlTw8PD44KPoZ87d8745z//aQwePNiIiIgwfHx8jAYNGhidOnUyXnjhBaOkpMSh/4cffmhER0cbnp6eDo9K79q1y4iNjTX8/f2NRo0aGQ899JCxffv2Ko9Tl5WVGWPHjjUaN25sWCyWKo+h/8///I8RExNj1K9f3wgICDDatWtnPPnkk8ahQ4eqrb9S5aPnhw8frnZ9RETEBR89LykpMSZMmGB06NDBCAgIMPz8/IwOHToYr7/+usM2J0+eNO677z4jKCjIkOTwu8jPzzdGjRplNGrUyPD29jbatWtX7WP6hw8fNu677z4jICDACAwMNEaOHGls3LjRkOTwKHhiYqLh5+dX7fHU9Hyfb4zevXsbN9xwQ43OE+AuFsNgBhkAmMXSpUt155136j//+c95P80a+L0h7ADAVerMmTMOk7DLy8vVv39/bd26VXl5eW6doA3UJczZAYCr1NixY3XmzBnZbDaVlJTo3//+tzZt2qTnn3+eoAP8Cld2AOAqlZGRoVmzZmnfvn06e/asWrdurTFjxjhMUAdA2AEAACbH5+wAAABTI+wAAABTY4KyfvnU0UOHDikgIMBlH48PAABql2EYOnHihJo0aVLlu99+jbAj6dChQwoPD3d3GQAA4BIcPHhQzZo1O+96wo5k//j1gwcPymq1urkaAABQE8XFxQoPD7f/PX4+hB3933fIWK1Wwg4AAFeZi01BYYIyAAAwNcIOAAAwNcIOAAAwNcIOAAAwNcIOAAAwNcIOAAAwNcIOAAAwNcIOAAAwNcIOAAAwNcIOAAAwNcIOAAAwNcIOAAAwNcIOAAAwNbeGnSlTpshisTgsUVFR9vVnz55VUlKSGjZsKH9/fw0dOlT5+fkOY+Tk5CghIUENGjRQSEiIJkyYoLKysit9KAAAoI7ydHcBN9xwg9auXWt/7en5fyWNHz9eH3/8sRYvXqzAwEAlJydryJAh2rhxoySpvLxcCQkJCgsL06ZNm5Sbm6sRI0bIy8tLzz///BU/FgAAUPe4Pex4enoqLCysSntRUZHefPNNZWRkqG/fvpKkefPmqW3btvryyy/VvXt3rV69Wrt27dLatWsVGhqqjh07avr06Xrqqac0ZcoUeXt7X+nDqSLy6Y/dXQJQp/0wI8HdJQAwObfP2dm7d6+aNGmili1bavjw4crJyZEkZWVl6dy5c4qNjbX3jYqKUvPmzZWZmSlJyszMVLt27RQaGmrvEx8fr+LiYu3cufO8+ywpKVFxcbHDAgAAzMmtYadbt26aP3++Vq5cqTlz5ujAgQO65ZZbdOLECeXl5cnb21tBQUEO24SGhiovL0+SlJeX5xB0KtdXrjuf1NRUBQYG2pfw8HDXHhgAAKgz3Hoba8CAAfaf27dvr27duikiIkL/+te/VL9+/Vrb78SJE5WSkmJ/XVxcTOABAMCk3D5n59eCgoJ03XXXad++fYqLi1NpaakKCwsdru7k5+fb5/iEhYXpq6++chij8mmt6uYBVfLx8ZGPj4/rDwDA7xbz84Dzc/fcPLfP2fm1kydPav/+/br22msVExMjLy8vrVu3zr5+z549ysnJkc1mkyTZbDbt2LFDBQUF9j5r1qyR1WpVdHT0Fa8fAADUPW69svPEE09o4MCBioiI0KFDhzR58mTVq1dP9957rwIDAzV69GilpKQoODhYVqtVY8eOlc1mU/fu3SVJ/fv3V3R0tO6//37NnDlTeXl5evbZZ5WUlMSVGwAAIMnNYeenn37Svffeq6NHj6px48bq2bOnvvzySzVu3FiSNHv2bHl4eGjo0KEqKSlRfHy8Xn/9dfv29erV0/LlyzVmzBjZbDb5+fkpMTFR06ZNc9chAQCAOsatYefdd9+94HpfX1+lp6crPT39vH0iIiK0YsUKV5cGAABMok7N2QEAAHA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADA1wg4AADC1OhN2ZsyYIYvFoscee8zedvbsWSUlJalhw4by9/fX0KFDlZ+f77BdTk6OEhIS1KBBA4WEhGjChAkqKyu7wtUDAIC6qk6EnS1btuiNN95Q+/btHdrHjx+vjz76SIsXL9aGDRt06NAhDRkyxL6+vLxcCQkJKi0t1aZNm7RgwQLNnz9fkyZNutKHAAAA6ii3h52TJ09q+PDh+uc//6lrrrnG3l5UVKQ333xTL730kvr27auYmBjNmzdPmzZt0pdffilJWr16tXbt2qV33nlHHTt21IABAzR9+nSlp6ertLTUXYcEAADqELeHnaSkJCUkJCg2NtahPSsrS+fOnXNoj4qKUvPmzZWZmSlJyszMVLt27RQaGmrvEx8fr+LiYu3cufO8+ywpKVFxcbHDAgAAzMnTnTt/9913tW3bNm3ZsqXKury8PHl7eysoKMihPTQ0VHl5efY+vw46lesr151Pamqqpk6depnVAwCAq4HbruwcPHhQf/nLX7Ro0SL5+vpe0X1PnDhRRUVF9uXgwYNXdP8AAODKcVvYycrKUkFBgW666SZ5enrK09NTGzZs0CuvvCJPT0+FhoaqtLRUhYWFDtvl5+crLCxMkhQWFlbl6azK15V9quPj4yOr1eqwAAAAc3Jb2OnXr5927Nih7Oxs+9K5c2cNHz7c/rOXl5fWrVtn32bPnj3KycmRzWaTJNlsNu3YsUMFBQX2PmvWrJHValV0dPQVPyYAAFD3uG3OTkBAgG688UaHNj8/PzVs2NDePnr0aKWkpCg4OFhWq1Vjx46VzWZT9+7dJUn9+/dXdHS07r//fs2cOVN5eXl69tlnlZSUJB8fnyt+TAAAoO5x6wTli5k9e7Y8PDw0dOhQlZSUKD4+Xq+//rp9fb169bR8+XKNGTNGNptNfn5+SkxM1LRp09xYNQAAqEvqVNhZv369w2tfX1+lp6crPT39vNtERERoxYoVtVwZAAC4Wrn9c3YAAABqE2EHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYGmEHAACYmtNhZ+HChSopKanSXlpaqoULF7qkKAAAAFdxOuyMGjVKRUVFVdpPnDihUaNGuaQoAAAAV3E67BiGIYvFUqX9p59+UmBgoEuKAgAAcBXPmnbs1KmTLBaLLBaL+vXrJ0/P/9u0vLxcBw4c0G233VYrRQIAAFyqGoedwYMHS5Kys7MVHx8vf39/+zpvb29FRkZq6NChLi8QAADgctQ47EyePFmSFBkZqbvvvlu+vr61VhQAAICr1DjsVEpMTJT0y9NXBQUFqqiocFjfvHlz11QGAADgAk6Hnb179+qBBx7Qpk2bHNorJy6Xl5e7rDgAAIDL5XTYGTlypDw9PbV8+XJde+211T6ZBQAAUFc4HXays7OVlZWlqKio2qgHAADApZz+nJ3o6GgdOXLEJTufM2eO2rdvL6vVKqvVKpvNpk8++cS+/uzZs0pKSlLDhg3l7++voUOHKj8/32GMnJwcJSQkqEGDBgoJCdGECRNUVlbmkvoAAMDVz+mw849//ENPPvmk1q9fr6NHj6q4uNhhcUazZs00Y8YMZWVlaevWrerbt68GDRqknTt3SpLGjx+vjz76SIsXL9aGDRt06NAhDRkyxL59eXm5EhISVFpaqk2bNmnBggWaP3++Jk2a5OxhAQAAk7IYhmE4s4GHxy/56LdzdVw1QTk4OFgvvPCChg0bpsaNGysjI0PDhg2TJO3evVtt27ZVZmamunfvrk8++UR33HGHDh06pNDQUEnS3Llz9dRTT+nw4cPy9vau0T6Li4sVGBiooqIiWa3Wy6r/tyKf/til4wFm88OMBHeX4BK814Hzq633eU3//nZ6zs5nn312WYWdT3l5uRYvXqxTp07JZrMpKytL586dU2xsrL1PVFSUmjdvbg87mZmZateunT3oSFJ8fLzGjBmjnTt3qlOnTrVSKwAAuHo4HXZ69+7t0gJ27Nghm82ms2fPyt/fXx988IGio6OVnZ0tb29vBQUFOfQPDQ1VXl6eJCkvL88h6FSur1x3PiUlJQ7f3O7s7TcAAHD1cDrsfP755xdc36tXL6fGu/7665Wdna2ioiItWbJEiYmJ2rBhg7NlOSU1NVVTp06t1X0AAIC6wemw06dPnyptv56/4+ycHW9vb7Vu3VqSFBMToy1btujll1/W3XffrdLSUhUWFjpc3cnPz1dYWJgkKSwsTF999ZXDeJVPa1X2qc7EiROVkpJif11cXKzw8HCn6gYAAFcHp5/GOn78uMNSUFCglStXqkuXLlq9evVlF1RRUaGSkhLFxMTIy8tL69ats6/bs2ePcnJyZLPZJEk2m007duxQQUGBvc+aNWtktVoVHR193n34+PjYH3evXAAAgDk5fWUnMDCwSltcXJy8vb2VkpKirKysGo81ceJEDRgwQM2bN9eJEyeUkZGh9evXa9WqVQoMDNTo0aOVkpKi4OBgWa1WjR07VjabTd27d5ck9e/fX9HR0br//vs1c+ZM5eXl6dlnn1VSUpJ8fHycPTQAAGBCToed8wkNDdWePXuc2qagoEAjRoxQbm6uAgMD1b59e61atUpxcXGSpNmzZ8vDw0NDhw5VSUmJ4uPj9frrr9u3r1evnpYvX64xY8bIZrPJz89PiYmJmjZtmqsOCwAAXOWcDjvffPONw2vDMJSbm6sZM2aoY8eOTo315ptvXnC9r6+v0tPTlZ6eft4+ERERWrFihVP7BQAAvx9Oh52OHTvKYrHot59F2L17d7311lsuKwwAAMAVnA47Bw4ccHjt4eGhxo0by9fX12VFAQAAuIrTYSciIqI26gAAAKgVTj96LkkbNmzQwIED1bp1a7Vu3Vp/+MMf9MUXX7i6NgAAgMvmdNh55513FBsbqwYNGmjcuHEaN26c6tevr379+ikjI6M2agQAALhkTt/Geu655zRz5kyNHz/e3jZu3Di99NJLmj59uu677z6XFggAAHA5nL6y89///lcDBw6s0v6HP/yhyuRlAAAAd3M67ISHhzt8hUOltWvX8v1SAACgznH6Ntbjjz+ucePGKTs7Wz169JAkbdy4UfPnz9fLL7/s8gIBAAAuh9NhZ8yYMQoLC9OsWbP0r3/9S5LUtm1bvffeexo0aJDLCwQAALgcl/TdWHfeeafuvPNOV9cCAADgck7P2dmyZYs2b95cpX3z5s3aunWrS4oCAABwFafDTlJSkg4ePFil/eeff1ZSUpJLigIAAHAVp8POrl27dNNNN1Vp79Spk3bt2uWSogAAAFzF6bDj4+Oj/Pz8Ku25ubny9LykKUAAAAC1xumw079/f02cOFFFRUX2tsLCQv31r39VXFycS4sDAAC4XE5finnxxRfVq1cvRUREqFOnTpKk7OxshYaG6u2333Z5gQAAAJfD6bDTtGlTffPNN1q0aJG2b9+u+vXra9SoUbr33nvl5eVVGzUCAABcskuaZOPn56eHH37Y1bUAAAC4nNNzdgAAAK4mhB0AAGBqhB0AAGBqhB0AAGBqToedli1b6ujRo1XaCwsL1bJlS5cUBQAA4CpOh50ffvhB5eXlVdpLSkr0888/u6QoAAAAV6nxo+fLli2z/7xq1SoFBgbaX5eXl2vdunWKjIx0aXEAAACXq8ZhZ/DgwZIki8WixMREh3VeXl6KjIzUrFmzXFocAADA5apx2KmoqJAktWjRQlu2bFGjRo1qrSgAAABXcfoTlA8cOFAbdQAAANSKGoWdV155pcYDjhs37pKLAQAAcLUahZ3Zs2fXaDCLxULYAQAAdUqNwg63rgAAwNWKT1AGAACm5vQE5QceeOCC6996661LLgYAAMDVnA47x48fd3h97tw5ffvttyosLFTfvn1dVhgAAIArOB12PvjggyptFRUVGjNmjFq1auWSogAAAFzFJXN2PDw8lJKSUuOntgAAAK4Ul01Q3r9/v8rKylw1HAAAgEs4fRsrJSXF4bVhGMrNzdXHH39c5TuzAAAA3M3psPP11187vPbw8FDjxo01a9asiz6pBQAAcKU5HXY+++yz2qgDAACgVjg9Z+fMmTM6ffq0/fWPP/6otLQ0rV692qWFAQAAuILTYWfQoEFauHChJKmwsFBdu3bVrFmzNGjQIM2ZM8flBQIAAFwOp8POtm3bdMstt0iSlixZorCwMP34449auHChU9+ODgAAcCU4HXZOnz6tgIAASdLq1as1ZMgQeXh4qHv37vrxxx9dXiAAAMDlcDrstG7dWkuXLtXBgwe1atUq9e/fX5JUUFAgq9Xq8gIBAAAuh9NhZ9KkSXriiScUGRmprl27ymazSfrlKk+nTp1cXiAAAMDlcPrR82HDhqlnz57Kzc1Vhw4d7O39+vXTnXfe6dLiAAAALtclfV1EWFiYAgICtGbNGp05c0aS1KVLF0VFRbm0OAAAgMvldNg5evSo+vXrp+uuu0633367cnNzJUmjR4/W448/7vICAQAALofTYWf8+PHy8vJSTk6OGjRoYG+/++67tXLlSpcWBwAAcLmcnrOzevVqrVq1Ss2aNXNob9OmDY+eAwCAOsfpKzunTp1yuKJT6dixY/Lx8XFJUQAAAK7idNi55ZZb7F8XIUkWi0UVFRWaOXOmbr31VpcWBwAAcLmcvo01c+ZM9evXT1u3blVpaamefPJJ7dy5U8eOHdPGjRtro0YAAIBL5vSVnRtvvFHff/+9evbsqUGDBunUqVMaMmSIvv76a7Vq1ao2agQAALhkTl3ZOXfunG677TbNnTtXzzzzTG3VBAAA4DJOXdnx8vLSN998U1u1AAAAuJzTt7H+9Kc/6c0336yNWgAAAFzO6bBTVlamOXPmqHPnznrkkUeUkpLisDgjNTVVXbp0UUBAgEJCQjR48GDt2bPHoc/Zs2eVlJSkhg0byt/fX0OHDlV+fr5Dn5ycHCUkJKhBgwYKCQnRhAkTVFZW5uyhAQAAE3L6aaxvv/1WN910kyTp+++/d1hnsVicGmvDhg1KSkpSly5dVFZWpr/+9a/q37+/du3aJT8/P0m/fGLzxx9/rMWLFyswMFDJyckaMmSI/cmv8vJyJSQkKCwsTJs2bVJubq5GjBghLy8vPf/8884eHgAAMBmLYRiGu4uodPjwYYWEhGjDhg3q1auXioqK1LhxY2VkZGjYsGGSpN27d6tt27bKzMxU9+7d9cknn+iOO+7QoUOHFBoaKkmaO3eunnrqKR0+fFje3t4X3W9xcbECAwNVVFQkq9Xq0mOKfPpjl44HmM0PMxLcXYJL8F4Hzq+23uc1/fv7kr71vLYUFRVJkoKDgyVJWVlZOnfunGJjY+19oqKi1Lx5c2VmZkqSMjMz1a5dO3vQkaT4+HgVFxdr586dV7B6AABQFzl9G6u2VFRU6LHHHtPNN9+sG2+8UZKUl5cnb29vBQUFOfQNDQ1VXl6evc+vg07l+sp11SkpKVFJSYn9dXFxsasOAwAA1DF15spOUlKSvv32W7377ru1vq/U1FQFBgbal/Dw8FrfJwAAcI86EXaSk5O1fPlyffbZZw7fph4WFqbS0lIVFhY69M/Pz1dYWJi9z2+fzqp8XdnntyZOnKiioiL7cvDgQRceDQAAqEvcGnYMw1BycrI++OADffrpp2rRooXD+piYGHl5eWndunX2tj179ignJ0c2m02SZLPZtGPHDhUUFNj7rFmzRlarVdHR0dXu18fHR1ar1WEBAADm5NY5O0lJScrIyNCHH36ogIAA+xybwMBA1a9fX4GBgRo9erRSUlIUHBwsq9WqsWPHymazqXv37pKk/v37Kzo6Wvfff79mzpypvLw8Pfvss0pKSpKPj487Dw8AANQBbg07c+bMkST16dPHoX3evHkaOXKkJGn27Nny8PDQ0KFDVVJSovj4eL3++uv2vvXq1dPy5cs1ZswY2Ww2+fn5KTExUdOmTbtShwEAAOowt4admnzEj6+vr9LT05Wenn7ePhEREVqxYoUrSwMAACZRJyYoAwAA1BbCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDXCDgAAMDW3hp3PP/9cAwcOVJMmTWSxWLR06VKH9YZhaNKkSbr22mtVv359xcbGau/evQ59jh07puHDh8tqtSooKEijR4/WyZMnr+BRAACAusytYefUqVPq0KGD0tPTq10/c+ZMvfLKK5o7d642b94sPz8/xcfH6+zZs/Y+w4cP186dO7VmzRotX75cn3/+uR5++OErdQgAAKCO83TnzgcMGKABAwZUu84wDKWlpenZZ5/VoEGDJEkLFy5UaGioli5dqnvuuUffffedVq5cqS1btqhz586SpFdffVW33367XnzxRTVp0uSKHQsAAKib6uycnQMHDigvL0+xsbH2tsDAQHXr1k2ZmZmSpMzMTAUFBdmDjiTFxsbKw8NDmzdvvuI1AwCAusetV3YuJC8vT5IUGhrq0B4aGmpfl5eXp5CQEIf1np6eCg4OtvepTklJiUpKSuyvi4uLXVU2AACoY+rslZ3alJqaqsDAQPsSHh7u7pIAAEAtqbNhJywsTJKUn5/v0J6fn29fFxYWpoKCAof1ZWVlOnbsmL1PdSZOnKiioiL7cvDgQRdXDwAA6oo6G3ZatGihsLAwrVu3zt5WXFyszZs3y2azSZJsNpsKCwuVlZVl7/Ppp5+qoqJC3bp1O+/YPj4+slqtDgsAADAnt87ZOXnypPbt22d/feDAAWVnZys4OFjNmzfXY489pr///e9q06aNWrRoob/97W9q0qSJBg8eLElq27atbrvtNj300EOaO3euzp07p+TkZN1zzz08iQUAACS5Oexs3bpVt956q/11SkqKJCkxMVHz58/Xk08+qVOnTunhhx9WYWGhevbsqZUrV8rX19e+zaJFi5ScnKx+/frJw8NDQ4cO1SuvvHLFjwUAANRNFsMwDHcX4W7FxcUKDAxUUVGRy29pRT79sUvHA8zmhxkJ7i7BJXivA+dXW+/zmv79XWfn7AAAALgCYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJgaYQcAAJiaacJOenq6IiMj5evrq27duumrr75yd0kAAKAOMEXYee+995SSkqLJkydr27Zt6tChg+Lj41VQUODu0gAAgJuZIuy89NJLeuihhzRq1ChFR0dr7ty5atCggd566y13lwYAANzsqg87paWlysrKUmxsrL3Nw8NDsbGxyszMdGNlAACgLvB0dwGX68iRIyovL1doaKhDe2hoqHbv3l3tNiUlJSopKbG/LioqkiQVFxe7vL6KktMuHxMwk9p437kD73Xg/GrrfV45rmEYF+x31YedS5GamqqpU6dWaQ8PD3dDNcDvW2CauysAUNtq+31+4sQJBQYGnnf9VR92GjVqpHr16ik/P9+hPT8/X2FhYdVuM3HiRKWkpNhfV1RU6NixY2rYsKEsFkut1gv3KS4uVnh4uA4ePCir1erucgDUEt7rvx+GYejEiRNq0qTJBftd9WHH29tbMTExWrdunQYPHizpl/Cybt06JScnV7uNj4+PfHx8HNqCgoJquVLUFVarlT8Agd8B3uu/Dxe6olPpqg87kpSSkqLExER17txZXbt2VVpamk6dOqVRo0a5uzQAAOBmpgg7d999tw4fPqxJkyYpLy9PHTt21MqVK6tMWgYAAL8/pgg7kpScnHze21aA9Mvty8mTJ1e5hQnAXHiv47csxsWe1wIAALiKXfUfKggAAHAhhB0AAGBqhB0AAGBqhB0AAGBqhB38bqSnpysyMlK+vr7q1q2bvvrqK3eXBMCF5syZo/bt29s/TNBms+mTTz5xd1moAwg7+F147733lJKSosmTJ2vbtm3q0KGD4uPjVVBQ4O7SALhIs2bNNGPGDGVlZWnr1q3q27evBg0apJ07d7q7NLgZj57jd6Fbt27q0qWLXnvtNUm/fKVIeHi4xo4dq6efftrN1QGoLcHBwXrhhRc0evRod5cCN+LKDkyvtLRUWVlZio2Ntbd5eHgoNjZWmZmZbqwMQG0pLy/Xu+++q1OnTslms7m7HLiZaT5BGTifI0eOqLy8vMrXh4SGhmr37t1uqgpAbdixY4dsNpvOnj0rf39/ffDBB4qOjnZ3WXAzruwAAEzj+uuvV3Z2tjZv3qwxY8YoMTFRu3btcndZcDOu7MD0GjVqpHr16ik/P9+hPT8/X2FhYW6qCkBt8Pb2VuvWrSVJMTEx2rJli15++WW98cYbbq4M7sSVHZiet7e3YmJitG7dOntbRUWF1q1bx718wOQqKipUUlLi7jLgZlzZwe9CSkqKEhMT1blzZ3Xt2lVpaWk6deqURo0a5e7SALjIxIkTNWDAADVv3lwnTpxQRkaG1q9fr1WrVrm7NLgZYQe/C3fffbcOHz6sSZMmKS8vTx07dtTKlSurTFoGcPUqKCjQiBEjlJubq8DAQLVv316rVq1SXFycu0uDm/E5OwAAwNSYswMAAEyNsAMAAEyNsAMAAEyNsAMAAEyNsAMAAEyNsAMAAEyNsAMAAEyNsAOgzvnhhx9ksViUnZ193j7r16+XxWJRYWHhFasLwNWJsAPgqtSjRw/7J+VK0vz58xUUFHRZY17qGAQvoG7j6yIAXJW8vb351noANcKVHQAuderUKY0YMUL+/v669tprNWvWLPXp00ePPfaYvY/FYtHSpUsdtgsKCtL8+fMd2nbv3q0ePXrI19dXN954ozZs2GBf9+urKevXr9eoUaNUVFQki8Uii8WiKVOmVFvf9u3bdeuttyogIEBWq1UxMTHaunXrBcd4++231blzZwUEBCgsLEz33XefCgoKJP1yy+3WW2+VJF1zzTWyWCwaOXKkpF++cTs1NVUtWrRQ/fr11aFDBy1ZsuSSzy2AS0PYAeBSEyZM0IYNG/Thhx9q9erVWr9+vbZt23bJYz3++OP6+uuvZbPZNHDgQB09erRKvx49eigtLU1Wq1W5ubnKzc3VE088Ue2Yw4cPV7NmzbRlyxZlZWXp6aeflpeX1wXHOHfunKZPn67t27dr6dKl+uGHH+yBJjw8XO+//74kac+ePcrNzdXLL78sSUpNTdXChQs1d+5c7dy5U+PHj9ef/vQnh9AGoPZxGwuAy5w8eVJvvvmm3nnnHfXr10+StGDBAjVr1uySxktOTtbQoUMlSXPmzNHKlSv15ptv6sknn3To5+3trcDAQFkslove2srJydGECRMUFRUlSWrTpo193fnGeOCBB+w/t2zZUq+88oq6dOmikydPyt/fX8HBwZKkkJAQ+5yfkpISPf/881q7dq1sNpt92//85z9644031Lt370s4IwAuBVd2ALjM/v37VVpaqm7dutnbgoODdf3111/SeJUhQZI8PT3VuXNnfffdd5dVY0pKih588EHFxsZqxowZ2r9//0W3ycrK0sCBA9W8eXMFBATYg0pOTs55t9m3b59Onz6tuLg4+fv725eFCxfWaJ8AXIewA+CKs1gsMgzDoe3cuXNXZN9TpkzRzp07lZCQoE8//VTR0dH64IMPztv/1KlTio+Pl9Vq1aJFi7RlyxZ7/9LS0vNud/LkSUnSxx9/rOzsbPuya9cu5u0AVxhhB4DLtGrVSl5eXtq8ebO97fjx4/r+++8d+jVu3Fi5ubn213v37tXp06erjPfll1/afy4rK1NWVpbatm1b7b69vb1VXl5eozqvu+46jR8/XqtXr9aQIUM0b968846xe/duHT16VDNmzNAtt9yiqKgo++TkX+9bksO20dHR8vHxUU5Ojlq3bu2whIeH16hOAK7BnB0ALuPv76/Ro0drwoQJatiwoUJCQvTMM8/Iw8Px31V9+/bVa6+9JpvNpvLycj311FPy8vKqMl56erratGmjtm3bavbs2Tp+/LjD/Jlfi4yM1MmTJ7Vu3Tp16NBBDRo0UIMGDRz6nDlzRhMmTNCwYcPUokUL/fTTT9qyZYt9XlB1YzRv3lze3t569dVX9eijj+rbb7/V9OnTHcaNiIiQxWLR8uXLdfvtt6t+/foKCAjQE088ofHjx6uiokI9e/ZUUVGRNm7cKKvVqsTExMs51QCcYQCAC504ccL405/+ZDRo0MAIDQ01Zs6cafTu3dv4y1/+Yu/z888/G/379zf8/PyMNm3aGCtWrDACAwONefPmGYZhGAcOHDAkGRkZGUbXrl0Nb29vIzo62vj000/tY3z22WeGJOP48eP2tkcffdRo2LChIcmYPHlyldpKSkqMe+65xwgPDze8vb2NJk2aGMnJycaZM2cuOEZGRoYRGRlp+Pj4GDabzVi2bJkhyfj666/t202bNs0ICwszLBaLkZiYaBiGYVRUVBhpaWnG9ddfb3h5eRmNGzc24uPjjQ0bNlzuaQbgBIth/ObGOQC4WJ8+fdSxY0elpaW5uxQAv0PM2QEAAKZG2AEAAKbGbSwAAGBqXNkBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACmRtgBAACm9v8BSoRHB/A2nb8AAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Pull of histogram counts from the result data structure\n", "counts = samples.histogram(key='result')\n", "print(counts)\n", "\n", "# Graph the histogram counts instead of the results\n", "cirq.plot_state_histogram(counts, plt.subplot())\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": { "id": "ef90dabdd52e" }, "source": [ "A histogram over the states that were actually observed can often be more useful when analyzing results. To learn more about the available options for creating result histograms, see the [State Histograms](../simulate/state_histograms.ipynb) page. " ] }, { "cell_type": "markdown", "metadata": { "id": "06Q_7vlQSu4Z" }, "source": [ "### Using parameter sweeps\n", "\n", "Cirq circuits allow for gates to have symbols as free parameters within the circuit. This is especially useful for variational algorithms, which vary parameters within the circuit in order to optimize a cost function, but it can be useful in a variety of circumstances.\n", "\n", "For parameters, Cirq uses the library `sympy` to add `sympy.Symbol` as parameters to gates and operations. \n", "\n", "Once the circuit is complete, you can fill in the possible values of each of these parameters with a `Sweep`. There are several possibilities that can be used as a sweep:\n", "\n", "* `cirq.Points`: A list of manually specified values for one specific symbol as a sequence of floats.\n", "* `cirq.Linspace`: A linear sweep from a starting value to an ending value.\n", "* `cirq.ListSweep`: A list of manually specified values for several different symbols, specified as a list of dictionaries.\n", "* `cirq.Zip` and `cirq.Product`: Sweeps can be combined list-wise by zipping them together or through their Cartesian product.\n", "\n", "A parameterized circuit and sweep together can be run using the simulator or other sampler by changing `run()` to `run_sweep()` and adding the sweep as a parameter.\n", "\n", "Here is an example of sweeping an exponent of a X gate:" ] }, { "cell_type": "code", "execution_count": 15, "metadata": { "execution": { "iopub.execute_input": "2023-05-18T09:44:57.129507Z", "iopub.status.busy": "2023-05-18T09:44:57.128860Z", "iopub.status.idle": "2023-05-18T09:44:58.095737Z", "shell.execute_reply": "2023-05-18T09:44:58.095008Z" }, "id": "ElyizofbLGq9" }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjcAAAGwCAYAAABVdURTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABgbklEQVR4nO3de1xUdf4/8NeAXBS5iAozGoEpWog3NA0vaykk6ddLtb9Ms7xlm2mS1K66myLaJnbx0mK5lZdcN7WLlaZR3l2NchNJCVRACC+MN4RBLoIz5/cHe0ZGZuCcYYa5vZ6PB4/icM7wOQ4z8z7n8/683wpBEAQQEREROQk3Ww+AiIiIyJIY3BAREZFTYXBDREREToXBDRERETkVBjdERETkVBjcEBERkVNhcENEREROpYWtB9DcdDodLl26BF9fXygUClsPh4iIiCQQBAFlZWXo0KED3NwavjfjcsHNpUuXEBISYuthEBERkRnOnz+Pe+65p8F9XC648fX1BVD7j+Pn52fj0RAREZEUGo0GISEh+s/xhrhccCNORfn5+TG4ISIicjBSUkqYUExEREROhcENERERORUGN0RERORUGNwQERGRU2FwQ0RERE6FwQ0RERE5FQY3RERE5FQY3BAREZFTYXBDRERETsXlKhQTkfPQ6gQcyy+GurQSxeXVCGjliZKKagS29kJQay9AAVzRVKG4vHab0s8b/TsFwt2NTXOJnJlNg5vDhw/j7bffxvHjx1FUVISvvvoK48aNa/CYgwcPIiEhAb/99htCQkLw+uuvY8qUKc0yXiKyD1qdgJT9udhwNB8llTWyjg1o6YGpg8Iwe1g4gxwiJ2XT4Ka8vBy9evXCtGnT8MQTTzS6f35+PkaNGoUXX3wR//73v7Fv3z48//zzUKlUGDFiRDOMmIiszdTdGPG/F0oq8fkvF3Dz1m2zHr+ksgYr9+bgn4fP4U9/uI9BDpETUgiCINh6EEBtI6zG7tzMmzcPu3btQmZmpn7b008/jZKSEqSmpkr6PRqNBv7+/igtLWXjTCI7IQY0e7LU+DrjEorLq5vtdwe08kDyEz0QF6nSj+NKWRWCfDmFRWRP5Hx+O1TOTVpaGmJiYgy2jRgxAq+88orJY27duoVbt27pv9doNNYaHhGZITWzCEk7s1BUWmWT319SUYMXN6fj/3qq8EvBDag1d8ah8vdG4ugIxEWqbDI2IjKPQ62WUqvVCA4ONtgWHBwMjUaDyspKo8csW7YM/v7++q+QkJDmGCoRSZCaWYSZm9NtFtjU9e3JIoPABgDUpVWYuTkdqZlFNhoVEZnDoYIbcyxYsAClpaX6r/Pnz9t6SESE2qmopJ1ZsIt5cRPEsSXtzIJWZ88jJaK6HGpaSqlU4vLlywbbLl++DD8/P7Rs2dLoMV5eXvDy8mqO4RGRDMfyi+3ijk1jBABFpVU4ll+M6M5tbT0cIpLAoYKb6Oho7N6922Dbnj17EB0dbaMREZE5tDoBm9IKbD0MWdSl9ae+mYBMZJ9sGtzcvHkTubm5+u/z8/ORkZGBwMBA3HvvvViwYAEuXryITZs2AQBefPFFpKSk4C9/+QumTZuG/fv347PPPsOuXbtsdQpE1Ii7l3Y3dSm3rfzt61M4mnsNg8LbQ+nnjRvl1Vi6yzARmgnIRPbBpkvBDx48iEceeaTe9smTJ2Pjxo2YMmUKCgoKcPDgQYNj5s6di6ysLNxzzz1YuHChrCJ+XApO1HyaayWUv3cLxEYEI7pzO6MVio/mXsPuTDUqqrVWHYd4z+aDSVEMcIgsTM7nt93UuWkuDG6Imoe4EspabzCBPh54vHdHxEQoJU0HiVWN1x85h9Iq6901UgBQ+nvjyLxhnKIisiCnrXNDRI7BEiuh/hDeFkO7BtWrUGxujyh3NwXiY8Ixe1gXpOzPxcq9Z5swOtOYgExkewxuiMjiLLESaubD4VYJDsQgp5uyNeZvP4WSCnm9qaS6Umb/K8GInJXT17khoual1Qk4mnutSY+h8q+9M2NNcZEqHH89FnNjuiKgpYfFHz/I19vij0lE0vDODRE1mSV7QykAJI6OaJZ8lbpTVeKKrqW7sps8fmUzBGdEZBqDGyJqEkuuiGrTygPL/tfEsjm5uyn0U2AtPd3x4ub0Jj1ecwVnRGQcgxsiMtvuk0V46dOmBQIAENDSA1MHhWH2sHCbBwVxkSqsnRRlVj5Oay93zBhyH2IjlFYaHRFJwaXgRGSW3ScvYfaWEzC35ZLcpdzNTVw6vuFoPkoq5QU59hSsETkL1rlpAIMboqZLzSxq0tTN7Ec6Y25sN4f44BfziX74rQgbfvxd1rEBrTyQbINpNiJnJOfzm6uliEgyrU7A0ZxrmP/lqSY9zqAu7R0isAHu5OMkjonE2klRUPlLXwVVUlGDFzenY/fJS9DqBKTlXcc3GReRlnedXcaJrIh3bohIEkskDjtD9V7xTs7R3KtIOZAn6RgFAP9WHgY5POxDRSQP79wQkUWJrRQssSLK0VcSiXdywoN9JR8jAPWSk9WlVZi5OR2pmUUWHiERMbghogZZopUCUHunwpkaSja1SJ/475m0M4tTVEQWxqXgRNQgc1opKAA8Fx2KewNbmd0Lyt717xQIlb93k+5msQ8VkXUwuCGiBqk18j+810zsg5E9O1hhNPbD3U2BxNERTS74B7APFZGlcVqKiExKzSzC0m9/k7y/yt8baydFOX1gIxIL/gW0alpvKvahIrIsBjdEZJSYRFxcLr2A3cJRDzhNTo1UYgPO+OHhMGfSrTmahBK5GgY3RFSPOUnECgBLd2W7ZHKsu5sCc2O7Ys3EKNnHOvrqMSJ7xOCGiOoxJ4m4bnKsqxrZs3aaSunn1ei+CgCvDA9nHyoiK2BwQ0T1NCXB1dWTY+MiVTg6fzjmxnRtcD8BwKp9ORi8fD9r3RBZGIMbIqqnKQmuTI6tnaaKjwmX1K6BxfyILI/BDRHVI9ZwkZMJogCTY+8WF6nCkXnD8O/pAxDQ0viKKuF/X4t3/OaS+UpE1sDghojqcXdTYOGoCMkJxWIQxOTY+tzdFHBzU6CksuFVZ2rNLaTsz22mURE5NxbxIyIAdxpCXimrQsG1Cmw5Vmh0P7GmS91eSUo2gWyQ1DyklXvPAhAwe1g4g0SiJmBwQ0SyOn6/OS4SIyJV+kAoyNf5WitYmpw8pJV7c7Dl2HksHsNgkchcnJYicnFyOn6LtWwAILpzW4zt3RHRndsysGmEmMMklVrDJGOipmBwQ+TC5BbrYy0b84h9qORix3Ai8zC4IXJh5hTrA1jLxhxxkSrMjQmXvL8YSG48ms8Ah0gms4Ob6upqnDlzBrdv37bkeIioGZnT8RtgLRtzzR4WDqWfvH+7pbuyWeiPSCbZwU1FRQWmT5+OVq1aoXv37igsrF1R8fLLLyM5OdniAyQi65Db8VvEWjbmc3dTYPGYCNkNNlnoj0ge2cHNggUL8Ouvv+LgwYPw9r5zBRITE4Nt27ZZdHBEZB3mdPwWsZZN08RFqvCBxP5TInFSijk4RNLIDm6+/vprpKSkYPDgwVAo7rzBde/eHXl5eRYdHBFZnjkdvwHATQG8PzGKy5MtQGr/qbqYzE0knezg5urVqwgKCqq3vby83CDYISL7ZG4SccqEPhjZk4GNpdTtPyUWRpRCXVppxVEROQfZwU2/fv2wa9cu/fdiQPPxxx8jOjraciMjIqvYm6WWtb/K3xtrJ0VhZM8OVhqRa4uNUMK7hbvk/ZfuymbuDVEjZFcofvPNN/HYY48hKysLt2/fxurVq5GVlYUff/wRhw4dssYYicgMddspiFWE92Spse5ogaTjZz/SBYO6tGP1YSs7ll8sa9XajfJqzNycjg8mcYqQyBTZwc3gwYORkZGB5ORk9OjRAz/88AOioqKQlpaGHj16WGOMRCSTsXYKSj8vVN3WNXqsArW9oubGdmVQ0wzk1gwSUPscJe3MQmyEks8RkRFm9Zbq3LkzPvroI0uPhYgsQFwJdXfCsFpzS9LxArgiqjmZUzOoboG/KYM68bkiuotZwY1Op0Nubi6uXLkCnc7wSvAPf/iDRQZGRPKZuxKqrmmDwjjd0YzEvlPq0irZz9vSXdn4+Eg+O7IT3UV2cPPTTz9h4sSJ+P333yEIhi9FhUIBrVZrscERkTzmroSqKzZCaaHRkBRi36mZm9OhAGQHOGKBP+bgEN0he7XUiy++iH79+iEzMxPFxcW4ceOG/qu4mPUXiGypKT2fFGD1YVvRF/aT0TlcxAJ/RPXJvnOTk5ODL774Al26dLHGeIioCZra84m5NrYTF6lCbIRSv8Kt4FoFVu49K+nYugX+oju3te5AiRyA7OBmwIAByM3NZXBDZIfMzd9wUwApEzitYWvuboq7ghMBK/fmSD6e3dqJaskObl5++WW8+uqrUKvV6NGjBzw8DCtr9uzZ02KDIyJ5zM3f0AlAGx9Paw6NzBDWzkfW/uzWTlRLdnDz5JNPAgCmTZum36ZQKCAIAhOKieyAmL9xd52bxvCq3/7ICVYCfTyg1lQhLe86Cy+Sy5Md3OTn51tjHETURHdXJP7ryAfw8pYTko/nVb/9kTPNWFxeg7nbMgDUJoZzeTi5MtnBTWhoqDXGQURNYKwisdQLd7EiMVdJ2R9zpxm5PJxcneyl4ACQl5eHl19+GTExMYiJicGcOXOQl5dn6bERkQRiReK7p6DkrArmKin7ZWqZuL93C7T2Mn59yuXh5Opk37n5/vvvMWbMGPTu3RuDBg0CABw9ehTdu3fHzp07ERsba/FBEpFxTa1IHNDKA8lP9ODVvZ27e5l4kK83dIKAZz7+2eQxXB5Orkx2cDN//nzMnTsXycnJ9bbPmzePwQ1RM2pqReI1E6IwKLydBUdE1nL3MvFvMi5KOo6J4uSKZE9LZWdnY/r06fW2T5s2DVlZWRYZFBFJY+4Hl1iN+CFe0TssqQngTBQnVyQ7uGnfvj0yMjLqbc/IyEBQUJAlxkREEpnzwSVm1jDPxrGJK6lMPYNsp0GuTPa01IwZM/DCCy/g3LlzGDhwIIDanJvly5cjISHB4gMkItPMqUis5DJhp9DYSioBwMJRDzCAJZekEO5u7d0IQRCwatUqvPvuu7h06RIAoEOHDvjzn/+MOXPmQKGw7xeSRqOBv78/SktL4efnZ+vhEDWZuFqqsRfy9EFhiIlQssCbkzFWBkDEejfkTOR8fssObuoqKysDAPj6+pr7EM2OwQ05o90nL2H2lhMml3+LtWyOzBvGwMZJ1C3amH+1HKv21e9BJT7TrHdDzkDO57fsaam6HCmoIXJG4gdcVpGmwbo2XBbsXBq6W1OXgNoAJ2lnFmIjlAxsyWU0KbghItuR+gFXF5cFOz6p05AiBrbkihjcEDmg3SeL8NKn6bKP47Jgx9aUoo0MbMmVMLghcjBifo0c7B/lHJpStJGBLbkSs3pLiaqqmn4lsGbNGoSFhcHb2xsDBgzAsWPHGtx/1apV6NatG1q2bImQkBDMnTvXIuMgcgSpmUV46VPTicPGsK6N8zDn7gvr3ZArkh3c6HQ6LF26FB07dkTr1q1x7tw5AMDChQuxbt06WY+1bds2JCQkIDExEenp6ejVqxdGjBiBK1euGN3/008/xfz585GYmIjs7GysW7cO27Ztw1//+le5p0HkcMQpCbmU/t5cLeMkzLn7wno35IpkBzdvvPEGNm7ciLfeeguenp767ZGRkfj4449lPdaKFSswY8YMTJ06FREREVi7di1atWqF9evXG93/xx9/xKBBgzBx4kSEhYXh0UcfxYQJExq823Pr1i1oNBqDLyJHJHdKYvYjXbBlxkM4Mm8YAxsn0VhVYlOW7spGamaRVcZEZI9kBzebNm3Chx9+iGeeeQbu7u767b169cLp06clP051dTWOHz+OmJiYO4Nxc0NMTAzS0tKMHjNw4EAcP35cH8ycO3cOu3fvxsiRI03+nmXLlsHf31//FRISInmMRPZEzpSEyt8bc2O7IrpzW16xOxGxKjGAegFOQ8+yurQKMzenM8AhlyE7uLl48SK6dOlSb7tOp0NNTY3kx7l27Rq0Wi2Cg4MNtgcHB0OtVhs9ZuLEiViyZAkGDx4MDw8PdO7cGQ8//HCD01ILFixAaWmp/uv8+fOSx0hkT+RMSTC/xnnFRarwwaQoKP0N/x6C/bwQ0MrD6DFiilbSzixo5SRsETko2aulIiIi8J///AehoaEG27/44gv06dPHYgMz5uDBg3jzzTfx/vvvY8CAAcjNzUV8fDyWLl2KhQsXGj3Gy8sLXl5eVh0XUXOQ0kfKTQGkTGB+jbOLi1QhNkKpr1Ac5OsNnSDgmY9/NnkM692QK5Ed3CxatAiTJ0/GxYsXodPpsH37dpw5cwabNm3Ct99+K/lx2rVrB3d3d1y+fNlg++XLl6FUKo0es3DhQjz77LN4/vnnAQA9evRAeXk5XnjhBfztb3+Dm1uTFn8R2bXGGiUCQMqEPhjZk4GNK3B3UxgEKd9kXJR0HOvdkCuQHQ2MHTsWO3fuxN69e+Hj44NFixYhOzsbO3fuRGxsrOTH8fT0RN++fbFv3z79Np1Oh3379iE6OtroMRUVFfUCGDHvpwktsogchqkpCZW/N9ZOisLInh1sNDKyNanTlu18vJCWdx3fZFxEWt51TlORUzKriN+QIUOwZ8+eJv/yhIQETJ48Gf369UP//v2xatUqlJeXY+rUqQCA5557Dh07dsSyZcsAAKNHj8aKFSvQp08f/bTUwoULMXr0aIPkZiJnpdUJ8G/pib+M6Ibi8moEtvaC0s+bnb5J0rQlAMz41y+oqNbqv2fncHJGZlco/uWXX5CdnQ2gNg+nb9++sh9j/PjxuHr1KhYtWgS1Wo3evXsjNTVVn2RcWFhocKfm9ddfh0KhwOuvv46LFy+iffv2GD16NP7+97+bexpEDsNYLynxg4mBDdWdtmxI3cAGuLOSirWQyJkoBJnzORcuXMCECRNw9OhRBAQEAABKSkowcOBAbN26Fffcc481xmkxclqmE9kLU80SxZCGH0wkEttzyK1irfT3xpF5wxgok92S8/ktO+fm+eefR01NDbKzs1FcXIzi4mJkZ2dDp9PpE32JyHIaapbIJb50tzY+XrICG8BwJRWRM5A9LXXo0CH8+OOP6Natm35bt27d8I9//ANDhgyx6OCIqPHKxFziS3U1ZTUUV1KRs5B95yYkJMRosT6tVosOHbhSg8jSpH7g8IOJgKZ1/y64VmHBkRDZjuzg5u2338bLL7+MX375Rb/tl19+QXx8PN555x2LDo6IpH9YNeVDjZyHuf2nAGDV3rNs0UBOQVJCcZs2baBQ3HmplJeX4/bt22jRonZWS/x/Hx8fFBfb95wtE4rJ3ml1gkHl2etltzBnm+kEUSaD0t1MJaA3hn9LZM/kfH5LyrlZtWqVJcZFRI0wttxbCi4Hp7rEYo9y/5aYv0XOQlJwM3nyZGuPg8jlmXO1zV5SZErd/lN7stT4OuMSisurJR3L/C1ydGzGRGQHGlru3RCdALTx8bTKmMjxif2nFo3ujv/+LQYLRz0g6Tjmb5GjY3BDZAcaW+7dEF5lkxTubgpMGdSpwWRjBWqrXvfvFNicQyOyOAY3RHagKQEKr7JJKrFFA4B6AY74PfO3yBkwuCGyA+YEKLzKJnOY6iyv9PdmGw9yGmY3zgRq+0wBsPt+UkT2TqxNwlVS1BzqJhuLJQfYWZ6ciew7NzqdDkuWLIG/vz9CQ0MRGhqKgIAALF26FDqdzhpjJHJ6dacLpFDxKpuaSEw2Htu7I6I7t2VgQ05F9p2bv/3tb1i3bh2Sk5MxaNAgAMCRI0ewePFiVFVV4e9//7vFB0nkCuIiVVg7KQrzt59CSUX9FieiuTHhmD0snB9GREQmSKpQXFeHDh2wdu1ajBkzxmD7N998g5deegkXL1606AAtjRWKyd5pdQJS9udiw9F8lFTeCXJU/t5IHB3BuzVE5JIsXqG4ruLiYtx///31tt9///1233qByBG4uykQHxOO2cO6MCeCrO7udh/8OyNnIDu46dWrF1JSUvDee+8ZbE9JSUGvXr0sNjAiVyfmRBBZi7F2H4E+Hni8d0fERCgZ6JDDkj0tdejQIYwaNQr33nsvoqOjAQBpaWk4f/48du/ejSFDhlhloJbCaSkiImntPjgVSvZEzue37NVSQ4cOxdmzZ/H444+jpKQEJSUleOKJJ3DmzBm7D2yIiEh6uw91aRVmbk5HamZRs4yLyFJkT0sVFhYiJCTE6KqowsJC3HvvvRYZGJGrYe4DNRep7T4E1BaLTNqZhdgIJf8eyWHIDm46deqEoqIiBAUFGWy/fv06OnXqBK1Wa7HBEbkKY7kPnBIga5HT7kMAUFRahWP5xcwBI4che1pKEAQoFPWj95s3b8Lbmz1uiOQScx/uvpLmlABZizntPtiglRyJ5Ds3CQkJAACFQoGFCxeiVatW+p9ptVr8/PPP6N27t8UHSOTMGsp94JQAWYvY7kNdWtVo3o2onY8X0vKuc9qUHILk4ObEiRMAau/cnDp1Cp6envqfeXp6olevXnjttdcsP0IiJ9ZY7gOnBMgaxHYfMzenN7qvAoB/Kw+8+vmvUGs4bUqOQXJwc+DAAQDA1KlTsXr1ai6jJrKAvVlqSftxSoAsTewOfneuV10K1AbYte1ADFuCiNOm7HFG9kh2QvGGDRusMQ4il5OaWYR1Rwsk7WtOjgRRY+p2B9+TpcbXGZdQXF6t/3mwnxeqbuuM9jrjtCnZM9nBDRGZT1zurS6txNJd2Y3urwCg9K/NbyCyBrESdnTntvjbqAj932dxeTWu3azGB4fyTB7LaVOyVwxuiJqJseXejREAJI6O4FUxNQt3NwVKK6vx1vdnZP2dctqU7A2DG6JmIKXUvTHTBoUxn4Gajbl/p5w2JXvD4IbIyqSWujcmNkJp8fEQGWPO3ymnTcleSQpuduzYIfkBx4wZY/ZgiJyNVidg49F8Wbf4AX5oUPOT2pLhbpw2JXskKbgZN26cwfcKhQJ1m4nXrVjM9gtEtczJsQFqAxuAHxrUvOTmzQS08kDyEz04bUp2SVL7BZ1Op//64Ycf0Lt3b3z33Xf6ruC7d+9GVFQUUlNTrT1eIodgqqWCFEp/b9YOoWYnN29mzQT+jZL9kp1z88orr2Dt2rUYPHiwftuIESPQqlUrvPDCC8jObnx5K5Eza0qOzdyYcMweFs47NtTsxJYMjQXk4pTpQ1z6TXZMduPMvLw8BAQE1Nvu7++PgoICCwyJyLGZm7ugALD1v+ctPyAiCcSWDFLCak6Zkr2THdw8+OCDSEhIwOXLl/XbLl++jD//+c/o37+/RQdH5IjMrflRtyAakS2ILRlU/sanqFScMiUHIXtaav369Xj88cdx7733IiQkBABw/vx5hIeH4+uvv7b0+IgcTlNrfrAgGtlS3ZYMYqXiwNZeUPrdWb3H7uBk72QHN126dMHJkyexZ88enD59GgDwwAMPICYmxmDVFJGrEnMX1KVVZuXdsCAa2ZrYkuFuxlYABrT0wNRBYcwVI7uiEOqu6ZapqqoKXl5eDhXUaDQa+Pv7o7S0lJ3NyWrE1VIAJAc4YqLmkXnD+CFBdqex6sVcGk7WJufzW3bOjU6nw9KlS9GxY0e0bt0a+fn5AICFCxdi3bp15o2YyMmIuQvKu3IXAlp5AEC9pE3WtiF7JmUFYElFDWZuTkdqZlGzjYvIFNnBzRtvvIGNGzfirbfegqenp357ZGQkPv74Y4sOjsiRxUWqcGTeMGyZ8RBWP90bW2Y8hOOvx2KtkaCHtW3InkldASgASNqZBa3O7AkBIouQnXOzadMmfPjhhxg+fDhefPFF/fZevXrpc3CIqJax3IW6CZtMyiRHICfJXVzxZyxnh6i5yA5uLl68iC5dutTbrtPpUFNTY5FBETkirU6QHLCYStgkskdyk9y54o9sTXZwExERgf/85z8IDQ012P7FF1+gT58+FhsYkSNpbBUJAN6pIYcltXqxiCv+yNZkBzeLFi3C5MmTcfHiReh0Omzfvh1nzpzBpk2b8O2331pjjER2zdQqkpLKGqzcm4N/Hj4HzxZuKKm4c2dT5e+NxNERzLEhhyBWL37xfysATWE3e7IXshOKx44di507d2Lv3r3w8fHBokWLkJ2djZ07dyI2NtYaYySyW1JWkVRUaw0CGwBQl1ZxZQk5lLhIFdZOitKv+LsbV/yRPZF15+b27dt48803MW3aNOzZs8daYyJyGOb2kRJQ+2GQtDMLsRFKfhiQQxCT4VP252LD0XyUVN4J2pW8G0l2RHYRv9atWyMzMxNhYWFWGpJ1sYgfWdI3GRcRvzWjSY+xZcZDTC4mhyMngZ7IEuR8fsvOuRk+fDgOHTrksMENkSVZInGSK0vIETW04o+BD9ma7ODmsccew/z583Hq1Cn07dsXPj4+Bj8fM2aMxQZHZO/kriIxhitLyJkYWznIBHpqbrKnpdzcTOcgKxQKaLXaJg/KmjgtRZaWmlnU6CoSY9hLipzN7pNFeOnT+q8F8a+bVbipKazeW8rUl70HNkTW0NgqEmO4soScze6TlzB7i/EgX7yCZmsGai6ygxsiqi8uUoXjr8dibkxXBLQ0DHICWnnUC3zYS4qchVYnYPXeHLz06Qk0FLcIuNOagcjaZOfcLFmypMGfL1q0yOzBEDkydzcF4mPCMXtYl3rJlAArFJPzSc0swuIdv0GtuSX5GCbQU3OQHdx89dVXBt/X1NQgPz8fLVq0QOfOnRnckMsxtjLE2CoSLvcmZ2KqMndjmEBPzUF2cHPixIl62zQaDaZMmYLHH3/cIoMichRcGUKuSEplbmNUbM1AzcQiOTd+fn5ISkrCwoULZR+7Zs0ahIWFwdvbGwMGDMCxY8ca3L+kpASzZs2CSqWCl5cXunbtit27d5s7dCKziVeudy8DZ2sFcnbmVuZmAj01F4slFJeWlqK0tFTWMdu2bUNCQgISExORnp6OXr16YcSIEbhy5YrR/aurqxEbG4uCggJ88cUXOHPmDD766CN07NjREqdAJFlDV65cGULOTm7ejALAK8PDERuhtM6AiO4ie1rqvffeM/heEAQUFRXhX//6Fx577DFZj7VixQrMmDEDU6dOBQCsXbsWu3btwvr16zF//vx6+69fvx7FxcX48ccf4eFRu/qksUrJt27dwq1bd5LdNBqNrDESGdPYlWvdlSHMtSFnIzdvRgCwal8Otv1ynlO21CxkBzcrV640+N7NzQ3t27fH5MmTsWDBAsmPU11djePHjxsc4+bmhpiYGKSlpRk9ZseOHYiOjsasWbPwzTffoH379pg4cSLmzZsHd3d3o8csW7YMSUlJksdFJIXUK1euDCFnJFbmVpdWycq7EadsWQaBrE12cJOfn2+RX3zt2jVotVoEBwcbbA8ODsbp06eNHnPu3Dns378fzzzzDHbv3o3c3Fy89NJLqKmpQWJiotFjFixYgISEBP33Go0GISEhFjkHcj3iyqiz6jJJ+3NlCDkjdzcFEkdHYObmdCgAyQGOgNopqqSdWYiNUDL/hqxGdnBzN41Gg/3796Nbt2544IEHLDEmk3Q6HYKCgvDhhx/C3d0dffv2xcWLF/H222+bDG68vLzg5eVl1XGRazC2MsoUsbUCV4aQs4qLVOGDSVH1XhOBPh4oLq8xeRynbKk5yA5unnrqKfzhD3/A7NmzUVlZiX79+qGgoACCIGDr1q148sknJT1Ou3bt4O7ujsuXLxtsv3z5MpRK40lnKpUKHh4eBlNQDzzwANRqNaqrq+Hp6Sn3dIgkkVPTg60VyFXERaoQG6E0qPOk1lRh7raMRo/llC1Zk+zVUocPH8aQIUMA1Bb0EwQBJSUleO+99/DGG29IfhxPT0/07dsX+/bt02/T6XTYt28foqOjjR4zaNAg5ObmQqfT6bedPXsWKpWKgQ1ZjdyaHmytQK7E3U2B6M5tMbZ3R0R3bguln7SpWE7ZkjXJDm5KS0sRGFh7qz01NRVPPvkkWrVqhVGjRiEnJ0fWYyUkJOCjjz7CJ598guzsbMycORPl5eX61VPPPfecQcLxzJkzUVxcjPj4eJw9exa7du3Cm2++iVmzZsk9DSLJ5Nb0eOePvRjYkMsSk41N3bNUgMX8yPpkT0uFhIQgLS0NgYGBSE1NxdatWwEAN27cgLe3vEh8/PjxuHr1KhYtWgS1Wo3evXsjNTVVn2RcWFgIN7c78VdISAi+//57zJ07Fz179kTHjh0RHx+PefPmyT0NIsnk3j6/Vi69zw6Rs2ko2ZhTttRcZAc3r7zyCp555hm0bt0aoaGhePjhhwHUTlf16NFD9gBmz56N2bNnG/3ZwYMH622Ljo7GTz/9JPv3EJlL7u1z3m4nV2cq2VjJ1iTUTBSCIMguoXr8+HEUFhYiNjYWrVu3BgDs2rULAQEBGDRokMUHaUkajQb+/v4oLS2Fn5+frYdDDkCrEzB4+f5Gp6bEFVJH5g3jVSkR7pROUJdWori8GoGtvRDU2gtQANdu3tI3muXrhaSQ8/ltVnDjyBjckDmkrJZSAEwkJrpLYyUU2GiWpLJ6cHPhwgXs2LEDhYWFqK6uNvjZihUr5D5cs2JwQ+Zq6E2ab9BE9Um9KAB4YUCNk/P5LTvnZt++fRgzZgzuu+8+nD59GpGRkfo6N1FRUWYPmsje1a3pUfc2u9KPt9aJ7ia1hAKrFpM1yA5uFixYgNdeew1JSUnw9fXFl19+iaCgIDzzzDOIi4uzxhiJ7IZY04OIGianhAKrFpOlya5zk52djeeeew4A0KJFC1RWVqJ169ZYsmQJli9fbvEBEhGR4zGnAjGrFpOlyA5ufHx89Hk2KpUKeXl5+p9du3bNciMjsiNanYC0vOv4JuMi0vKuQ6tzqTx8ItnMKYnAMgpkKbKnpR566CEcOXIEDzzwAEaOHIlXX30Vp06dwvbt2/HQQw9ZY4xENmUskZgJxEQNEysVq0urJLUucVMAN8qrG9+RSALZd25WrFiBAQMGAACSkpIwfPhwbNu2DWFhYVi3bp3FB0hkS+Jqj7tzB9SlVZi5OR2pmUU2GhmRfRMrFQMw2YqhLp0AzPqUrymyDNa5ITJBqxMwKHk/1BrjeQAs2kfUuMbq3NTF1xQ1RM7nt+w7NwBQUlKCjz/+GAsWLEBxcTEAID09HRcvXjTn4YjsUsr+HJOBDWC4woOIjIuLVOHIvGHYMuMhzH6kc4P78jVFliI75+bkyZOIiYmBv78/CgoKMGPGDAQGBmL79u0oLCzEpk2brDFOomaVmlmElXuldbnnCg+ihoklFKS+VviaoqaSfecmISEBU6ZMQU5OjkEX8JEjR+Lw4cMWHRyRLYjFx6TiCg8iaaS+VviaoqaSHdz897//xZ/+9Kd62zt27Ai1Wm2RQRHZkpziYyr/2urERNQ4cQWVqWwaBfiaIsuQHdx4eXlBo9HU23727Fm0b9/eIoMisiU5t8QTR0cw8ZFIosZWUAkAFo56gK8pajLZwc2YMWOwZMkS1NTUAAAUCgUKCwsxb948PPnkkxYfIFFzk3pLfG5MV9a5IZIpLlKFDyZFQelv/HX2168ysXrvWRbKpCaRHdy8++67uHnzJoKCglBZWYmhQ4eiS5cu8PX1xd///ndrjJGoWTV26xwAlH5emD2sS7ONiciZxEWqsHBUhNGflVTWYOXeHPR9Yw9r3pDZzK5zc+TIEZw8eRI3b95EVFQUYmJiLD02q2CdG5JCLN4HwKC6qhjwfDApindtiMyk1QkYvHx/o7ltCvC1RnfI+fxmET8iE9h2gcg60vKuY8JHP0naV8WifvQ/cj6/Zde5AWpXTB04cABXrlyBTqcz+NmKFSvMeUgiuxMXqUJshBLH8otxpawKQb61qzj4JkvUNHKS9sWiftGd21pxRORsZAc3b775Jl5//XV069YNwcHBUCjuvNHX/X8iZyAWHyMiy5Fbx4ZF/Ugu2cHN6tWrsX79ekyZMsUKwyGyPa1O4N0aIisSk/al1pNiUT+SS3Zw4+bmhkGDBlljLEQ2xzwbIusT6928+L+k/YYE+nigb2ibZhgVORPZS8Hnzp2LNWvWWGMsRDYlrpC6+2pSXVqFmZvTuSyVyILiIlVYOykKAa08GtyvuLwGQ98+wNcfySJ7tZROp8OoUaNw9uxZREREwMPD8A9z+/btFh2gpXG1FBnT2NJUBQAlV20QWZxWJyBlfy42HM1HSWWN0X1YgoEAeZ/fsu/czJkzBwcOHEDXrl3Rtm1b+Pv7G3wROaLG+kkJuLNqg4gsx91NgfiYcBz7WwwCfTyN7iNegSftzGLlYpJEds7NJ598gi+//BKjRo2yxniImk3dxOHUTGlNX7lqg8g6jv9+A8Xl1SZ/XvcCgysYqTGyg5vAwEB07tzZGmMhajbGEoel4KoNIuuQeuGwJ0vN4IYaJXtaavHixUhMTERFRYU1xkNkdaYShxuiQO2qqf6dAq03MCIXJvXCYf3RAiYXU6Nk37l57733kJeXh+DgYISFhdVLKE5Pb3xpH5GtaHUCknZmQe6svQAgcXQEk4mJrERO7ZvFO35DbISSr0cySXZwM27cOCsMg6h5NJY4bMq0QWFcpUFkRXJq36g1t5CyPxfxMeHNMDJyRLKDm8TERGuMg6hZmJsQHBuhtPBIiOhucZEqTB8UhnVHCxrdd+Xes+imbM2LDjJKds4NkSOTmxDMXBui5hUj40KCS8PJFAY35FL6hrYxWUvjbuJsPnNtiJqPmHsjBWtPkSkMbshlpGYWYejbBxqspVGX0t+bFVGJmpmYeyMVa0+RMQxuyCWYs/x74agHGNgQ2UBcpApzJSYLs/YUGSM7uFmyZInRGjeVlZVYsmSJRQZFZElanYDFO+Qt/1YAWLorm/P5RDYye1g4lH4NBy7sGE6myA5ukpKScPPmzXrbKyoqkJSUZJFBEVlSyv4cqDXybl2zlxSRbbm7KbB4TAQUuJP/djd2DCdTZAc3giBAoaj/p/brr78iMJArSsi+pGYWYeXeHLOP53w+ke3ERarwwaQoKBtIMFaXVmHm5nQGOGRAcp2bNm3aQKFQQKFQoGvXrgYBjlarxc2bN/Hiiy9aZZBE5hCrETcF5/OJbCsuUoVh9wfjoWX7jC4GEFB7ZydpZxarFpOe5OBm1apVEAQB06ZNQ1JSEvz9/fU/8/T0RFhYGKKjo60ySCJzmFuNGKh9s1Syvg2RXWDHcJJLcnAzefJkAECnTp0wcODAej2liOyNuVNKrG9DZF+kvpY5jUwiScGNRqOBn58fAKBPnz6orKxEZWWl0X3F/YhsTeqUUmuvFrh567b+e6W/NxJHR3AZOJGdkPpa5jQyiSQFN23atEFRURGCgoIQEBBgNKFYTDTWarUWHySROcRKp+rSKpPLwJV+Xjj8l2E4/vsNXCmrQpBv7VQU79gQ2Q+pr2VOI5NIUnCzf/9+/UqoAwcOWHVARJYiVjqduTkdCsDgTVEMXRaP6Q7PFm6cpyeyYw29lkVVt3XYk6XmHVcCACgEQXCpKmUajQb+/v4oLS3lFJqLSM0sQtLOLIPkYhWnnogcTmpmEeZvP4WSipp6PxMvWNgyxXnJ+fw2K7i5ceMG1q1bh+zsbABAREQEpk6d6hB1bhjcuCatTsCx/GJOPRE5MK1OwKDk/SaLcoqrHI/MG8bXtxOS8/ktu4jf4cOHERYWhvfeew83btzAjRs38N5776FTp044fPiw2YMmsiZ3NwWiO7fF2N4dEd25Ld/4iBzQsfziBquNs7I4iSQvBRfNmjUL48ePxwcffAB3d3cAtUX8XnrpJcyaNQunTp2y+CCJzMU7NkTOg0vCSSrZwU1ubi6++OILfWADAO7u7khISMCmTZssOjiipmCuDZFz4ZJwkkr2tFRUVJQ+16au7Oxs9OrVyyKDImqq1MwizNycXq9CMfvQEDkucUl4Q/de3RTAjQaqGZNrkHTn5uTJk/r/nzNnDuLj45Gbm4uHHnoIAPDTTz9hzZo1SE5Ots4oiWQQe0oZy5RnHxoix1V3SbgpOgGY9Wk6PnDjqilXJmm1lJubGxQKBRrb1RGK+HG1lPNLy7uOCR/91Oh+W2Y8xPo2RA5o98lLmL3lBHQmPpK4aso5yfn8lnTnJj8/3yIDI2oOe7PUkvZj0iGRY2rj42UysAHYSJMkBjehoaHWHgeRRaRmFmHd0QJJ+zLpkMgxcdUUNUb2aqnGVkQ999xzZg+GqCnEXJvGiLes2YeGyDFx1RQ1RnZwEx8fb/B9TU0NKioq4OnpiVatWpkV3KxZswZvv/021Go1evXqhX/84x/o379/o8dt3boVEyZMwNixY/H111/L/r3kXH46d73e6ihjBACJoyM4F0/koBprpMkLGJK9FFysSix+3bx5E2fOnMHgwYOxZcsW2QPYtm0bEhISkJiYiPT0dPTq1QsjRozAlStXGjyuoKAAr732GoYMGSL7d5LzSc0swqx/m15BUde0QWFcRUHkwMRVUwDqLQsXv+cFjGuTHdwYEx4ejuTk5Hp3daRYsWIFZsyYgalTpyIiIgJr165Fq1atsH79epPHaLVaPPPMM0hKSsJ9993XlKGTExBr2pRU1m+mZ0xshNLKIyIia4uLVOGDSVFQ+htOPSn9vdk8k+RPS5l8oBYtcOnSJVnHVFdX4/jx41iwYIF+m5ubG2JiYpCWlmbyuCVLliAoKAjTp0/Hf/7znwZ/x61bt3Dr1i399xqNRtYYyb41VNPmbrxVTeRc4iJViI1QssUK1SM7uNmxY4fB94IgoKioCCkpKRg0aJCsx7p27Rq0Wi2Cg4MNtgcHB+P06dNGjzly5AjWrVuHjIwMSb9j2bJlSEpKkjUuchzH8osl5dmIeKuayLmITXGJ6pId3IwbN87ge4VCgfbt22PYsGF49913LTUuo8rKyvDss8/io48+Qrt27SQds2DBAiQkJOi/12g0CAkJsdYQqZlJXeoZ0MoDyU/04K1qIidVt0luOx8vQAFcu3mLd3NclOzgRqfTWeyXt2vXDu7u7rh8+bLB9suXL0OprJ8XkZeXh4KCAowePbreeFq0aIEzZ86gc+fOBsd4eXnBy8vLYmMm+1JwrVzSfmsmRGFQuLSAmIgci7EmuXWxYa7rsUhCsbk8PT3Rt29f7Nu3T79Np9Nh3759iI6Orrf//fffj1OnTiEjI0P/NWbMGDzyyCPIyMjgHRkXs/tkEVbuzWlwHwVq39ge4m1rIqdkqkluXWyY63pk37mpO8XTmBUrVkh6vMmTJ6Nfv37o378/Vq1ahfLyckydOhVAbVHAjh07YtmyZfD29kZkZKTB8QEBAQBQbzs5N7G3jBTMsyFyTlIXFIg/n//lKfh6e+Ch+9ryPcHJyQ5uTpw4gRMnTqCmpgbdunUDAJw9exbu7u6IiorS76dQSPvDGT9+PK5evYpFixZBrVajd+/eSE1N1ScZFxYWws3NpjeYyM6kZhbhpU+lBTavxHTlrWgiJyV3QUFJZQ2e+fhnTlO5AEldwetasWIFDh48iE8++QRt2rQBUFvYb+rUqRgyZAheffVVqwzUUtgV3LFpdQIGL98v+Q1t9dO9MbZ3RyuPiohs4ZuMi4jfmiH7OPHSm/VwHIucz2/Zt0TeffddLFu2TB/YAECbNm3wxhtvWH21FJHcKzX2liFyXua+vsUr+qSdWdA21F6cHJbs4Eaj0eDq1av1tl+9ehVlZWUWGRSRKWqN9MBGxYJ9RE5N7DFlTvaMAKCotArH8ostPSyyA7KDm8cffxxTp07F9u3bceHCBVy4cAFffvklpk+fjieeeMIaYyQCUJtrs/Tb3yTvz0RiIufWUI8pqaTWyiLHIju4Wbt2LR577DFMnDgRoaGhCA0NxcSJExEXF4f333/fGmMk0i/3LC5vvH+UmwJ4fyLn0olcgakeU1Jx6to5yU4oFpWXlyMvLw8A0LlzZ/j4+Fh0YNbChGLHIzeJ+P2JfTCyZwcrj4qI7IlYofho7lWkHMiTdIzK3xtH5g3jHV4HIefz2+zGmT4+PujZs6e5hxNJJjWJuK2PJ/7+eCTv2BC5ILHHlJxpJk5dOy8WkCG7J/XN6vVRDzCwIXJxUqeZ5rIGllNjcEN2T+qbldK/pZVHQkT2TsoKKqWfF2YP69JsY6Lmx+CG7F5jb1Zi/ygu+yaihlZQKf73tXhMd05HOTlJwU1UVBRu3LgBAFiyZAkqKiqsOiiiuhp7swI4d05Ed5haQaX092ZVYhchabVUy5YtkZOTg3vuuQfu7u4oKipCUFBQc4zP4rhaynGlZhYhaWeWQXIxe8QQkSniCqorZVUI8q29u8uLIMdl8dVSvXv3xtSpUzF48GAIgoB33nkHrVu3NrrvokWL5I+YqBFanQD/lp74y4huKC6vRmBrLyj9+GZFRKaJK6jI9UgKbjZu3IjExER8++23UCgU+O6779CiRf1DFQoFgxuyuIbu2DCwISKiu8ku4ufm5ga1Ws1pKWoWYmXiu/9I2dWXiMi1WLWIn06nM3tgRHJodQIW78iqF9gAtU3vFKjt6hsboeQdHCJqFHNwXIdZFYrz8vKwatUqZGdnAwAiIiIQHx+Pzp07W3Rw5NpS9uc02AW8bldfzqsTUUO4IMG1yK5z8/333yMiIgLHjh1Dz5490bNnT/z888/o3r079uzZY40xkgtKzSzCyr05kvZlV18iaog4vX13Gxd1aRVmbk5HamaRjUZG1iL7zs38+fMxd+5cJCcn19s+b948xMbGWmxw5Jq0OgFJO7Mk78+uvkRkivh+wult1yL7zk12djamT59eb/u0adOQlSX9A4nIFKmNMgFWJiaihjX2fiJOb/+Ud735BkVWJzu4ad++PTIyMuptz8jIcNgVVGRf2NWXiCxF6vvJrE85PeVMZE9LzZgxAy+88ALOnTuHgQMHAgCOHj2K5cuXIyEhweIDJNfDrr5EZClS309KKmvw4uZ0vD+xD0b27GDlUZG1yQ5uFi5cCF9fX7z77rtYsGABAKBDhw5YvHgx5syZY/EBkusRG2WqS6uMzpMD7OpLRNJIeT+pa/aWE0iBAiN78sLJkcku4ldXWVkZAMDX19diA7I2FvFzDOLqBgAGb0gs3kdEcpkqBtqQtXyPsTtyPr9l59zU5evr61CBDdk3rU5AWt51fJNxEb5eHogfHg7/lh4G+7CrLxHJJXYJD7jr/aQhSTuzoNWZfe1PNmZWET8iSzNWYKuugJYemDooDLOHhTOBmIhki4tUwdfbA898/LOk/Vkg1LE16c4NkSWYKrBVV2llDVbtzcGeLHUzjoyInMlD97WFyl96XSwWCHVcDG7IphoqsFWX+HPeKiYic7m7KZA4OkLy/iwQ6rhkBzfnzp2zxjjIRckp2Fe3lxQRkTniIlV4f2IfNDS7rQALhDo62cFNly5d8Mgjj2Dz5s2oquItO2oac2778lYxETXFyJ4dkDIhyujPxJiHBUIdm+zgJj09HT179kRCQgKUSiX+9Kc/4dixY9YYGzk5rU7AtbJbso/jrWIiaqqRPVVYOymqXg4OV2Q6B7Pr3Ny+fRs7duzAxo0bkZqaiq5du2LatGl49tln0b59e0uP02JY58Y+NLY6yhgFat94jswbxisqIrIIrU7AsfxiqEsrUVxejcDWXlD61U5J8X3Gvsj5/G5SET8AuHXrFt5//30sWLAA1dXV8PT0xFNPPYXly5dDpbK/yJfBje2ZU1CLxfuIyFqMXWyp/L2RODqC7zd2pFmK+P3yyy946aWXoFKpsGLFCrz22mvIy8vDnj17cOnSJYwdO9bchyYnJnV11N14q5iIrMFUKYqi0iq8uDkdu09estHIqClkF/FbsWIFNmzYgDNnzmDkyJHYtGkTRo4cCTe32jipU6dO2LhxI8LCwiw9VnICclZHvTK8Czq1b40gX94iJiLLk3KxxV5Tjkl2cPPBBx9g2rRpmDJlislpp6CgIKxbt67JgyPnI2el07ZfLjC/hoisRsrFlk4AXvo0HXOvhLNCugORHdzk5OQ0uo+npycmT55s1oDIuclZ6cTy50RkTXIutlbuzcGWY+exeAzzcByB7JybDRs24PPPP6+3/fPPP8cnn3xikUGR8+rfKZDlz4nILsgtK6HWVGHm5nSkZhZZaURkKbKDm2XLlqFdu3b1tgcFBeHNN9+0yKDIebH8ORHZC7kXWyK2gbF/soObwsJCdOrUqd720NBQFBYWWmRQ5NxY/pyI7IHciy2AbWAchezgJigoCCdPnqy3/ddff0XbtsyNIGlY/pyI7IGUiy1jOGVu32QHNxMmTMCcOXNw4MABaLVaaLVa7N+/H/Hx8Xj66aetMUZyUix/TkT2oKGLLVM4ZW7fZFcorq6uxrPPPovPP/8cLVrULrbS6XR47rnnsHbtWnh6elploJbCCsX2Ryx/fqWsijVtiMhmUjOLsHjHb1BrTPe8YxsY22mW9gtnz57Fr7/+ipYtW6JHjx4IDQ01a7DNjcENERGZotUJSNmfi5V7z9b7GdvA2Faz9pZyNAxubI93aojI3hnrN6X088KE/vcirJ0P37tswKrBjVarxcaNG7Fv3z5cuXIFOp3O4Of79++XP+JmxODGttigjogcRd0LsYJrFdhyrBBqDd+7bMWqjTPj4+MRHx8PrVaLyMhI9OrVy+CLyBitTsDqvTl40UiDOnUpC2MRkf1xd1MgunNbeLi5YeXeswaBDcD3Lnsm+85Nu3bt9M0yHRHv3DQ/JukRkaPaffISZm85gYZq9in9vHB0/nC+d1mZVe/ceHp6okuXLmYPjlxLamYRZm5ObzCwAVgYi4jsT2pmEV76tOHABgDUmluY/Wk60vKus3KxnZAd3Lz66qtYvXo1XCwPmcyg1QlI2pkFOX8pLIxFRPZAfP+S6rtMNSZ89BMGL9/PaSo7ILsr+JEjR3DgwAF899136N69Ozw8PAx+vn37dosNjhzbsfzievk1jWFhLCKyB+a8fwF38nC4XNy2ZAc3AQEBePzxx60xFnIycu7CiDk37CVFRPbA3LvI4p3qxTt+Q2yEknk4NiI7uNmwYYM1xkFOSO5dGPaSIiJ70dS7yGrNLaTsz0V8TLiFRkRyyM65AYDbt29j7969+Oc//4mysjIAwKVLl3Dz5k2LDo4c243yW5Ka0anYS4qI7Ez/ToFQ+XujKZdbK/eeZf6NjcgObn7//Xf06NEDY8eOxaxZs3D16lUAwPLly/Haa69ZfIDkmFIzizBLwiqDuTHhODJvGAMbIrIr7m4KJI6OAIAmBThJO7O4gsoGzCri169fP9y4cQMtW7bUb3/88cexb98+iw6OHJOUVVJuCuD9iVGIj+nKqSgisktxkSp8MCkKSn/zp6hY4sI2ZOfc/Oc//8GPP/5Yr/t3WFgYLl68aLGBkeOSsspAJwBtfOy7gzwRUVykCrERygbbMDSGJS6an+w7NzqdDlqttt72CxcuwNfX16xBrFmzBmFhYfD29saAAQNw7Ngxk/t+9NFHGDJkCNq0aYM2bdogJiamwf2p+Ul9IfMFT0SOQGzDMLZ3R8THhOPo/GHYMuMhPBYZLOl4lrhofrKDm0cffRSrVq3Sf69QKHDz5k0kJiaa1ZJh27ZtSEhIQGJiItLT09GrVy+MGDECV65cMbr/wYMHMWHCBBw4cABpaWkICQnBo48+yrtGdkTqC5kveCJyRGKwkzKxL5R+Db+PBfp4QK2pYvXiZia7t9SFCxcwYsQICIKAnJwc9OvXDzk5OWjXrh0OHz6MoKAgWQMYMGAAHnzwQaSkpACovTMUEhKCl19+GfPnz2/0eK1WizZt2iAlJQXPPfdco/uzt5T1NdaLhX2kiMhZiC1mADRajZ1dxJvGqr2l7rnnHvz666/461//irlz56JPnz5ITk7GiRMnZAc21dXVOH78OGJiYu4MyM0NMTExSEtLk/QYFRUVqKmpQWCg8eJvt27dgkajMfgi65G6Soo1bYjIGchJOmYX8eYjO6EYAFq0aIFJkyY1+Zdfu3YNWq0WwcGG85bBwcE4ffq0pMeYN28eOnToYBAg1bVs2TIkJSU1eaxkmlYn4Fh+MdSllVi6K7vRVVIpE1jThoicR92kY/F9sLi8ut5+AmrvXCftzGL1YiuTHdxs2rSpwZ9LmRqylOTkZGzduhUHDx6Et7fxqHnBggVISEjQf6/RaBASEtJcQ3R6qZlFSNqZJbkHC1dJEZEzEvNw0vKuGw1sRALuLA+P7ty2+QboYmQHN/Hx8Qbf19TUoKKiAp6enmjVqpWs4KZdu3Zwd3fH5cuXDbZfvnwZSqWywWPfeecdJCcnY+/evejZs6fJ/by8vODl5SV5TCSdONcsN0WOq6SIyFlxtah9kJ1zc+PGDYOvmzdv4syZMxg8eDC2bNki67E8PT3Rt29fg+J/Op0O+/btQ3R0tMnj3nrrLSxduhSpqano16+f3FMgC9DqBCze0XChPlO4SoqInBVXi9oHs3pL3S08PBzJycn17upIkZCQgI8++giffPIJsrOzMXPmTJSXl2Pq1KkAaqe5FixYoN9/+fLlWLhwIdavX4+wsDCo1Wqo1Wr2tWpGWp2AeV+clFXECqida1ax8zcROTEpPancFMCNBqauqOksEtwAtUnGly5dkn3c+PHj8c4772DRokXo3bs3MjIykJqaqk8yLiwsRFHRnczyDz74ANXV1fjjH/8IlUql/3rnnXcsdSrUgNTMIvR9Yw++SL8g+1gBXCVFRM6tbk8qU3QCMOtTrpqyJtl1bnbs2GHwvSAIKCoqQkpKCkJCQvDdd99ZdICWxjo35kvNLMKL/6vnYI5pg8KwaHR3C46IiMg+sd6X5cn5/JadUDxu3DiD7xUKBdq3b49hw4bh3Xfflftw5CDEZphNERvRcJI4EZGzaOPj1WC9L66asi7ZwY1Op7PGOMjOSWmGaYp4hcJcGyJydmLdr+8kTjlx1ZR1mFXEj1yPuS9A8WYrc22IyNnJrfsFcNWUtcgObuoWxGvMihUr5D482al2rc2rFaRkLxUicgFy637xjrZ1yQ5uTpw4gRMnTqCmpgbdunUDAJw9exbu7u6IiorS76dQ8CrdWaRmFmHxjt8k7z83Jhxh7XwQ5Fv7wuUdGyJyZmJOopzVOQKAkZG1LRv4Pml5soOb0aNHw9fXF5988gnatGkDoLaw39SpUzFkyBC8+uqrFh8k2Y6cq5E2rTyw7IkevEtDRC5Fbk6im6J2Ofi6owVYd7SA3cKtQPZS8I4dO+KHH35A9+6GS3ozMzPx6KOPmlXrpjlxKbh0Wp2Awcv3N/qiDWjpgamDwjB7WDivPojI5XyTcRHxWzMa3W9o1/Y4dPZqve3iu+YHk9hUuCFyPr9lF/HTaDS4erX+k3P16lWUlZXJfTiyY1KvRtZMjEJ8TFcGNkTkkqQmBZ+6WGp0u3iHIWlnFrQNrR8nyWQHN48//jimTp2K7du348KFC7hw4QK+/PJLTJ8+HU888YQ1xkg2InWF1LXyW1YeCRGR/Wqs5YICQKCPh+Ru4dR0soObtWvX4rHHHsPEiRMRGhqK0NBQTJw4EXFxcXj//fetMUayETaAIyJqXN2WC3cHOOL3fUICJD0W695YhuzgplWrVnj//fdx/fp1/cqp4uJivP/++/Dx8bHGGMlG+ncKhNLPdODCRphERLXiIlX4YFIUlP6G75lKf2+88IdO2He6fjqHMbxYtAyzi/gVFRWhqKgIf/jDH9CyZUsIgsDl305mT5YaVbe1Rn/G4nxERIbiIlWIjahd3n2lrApBvt7oG9oGQ98+0OixrHtjWbKDm+vXr+Opp57CgQMHoFAokJOTg/vuuw/Tp09HmzZt2F/KSTS2BDyAy76JiOpxd1MY9IpKy7suaWGGAF4sWpLsaam5c+fCw8MDhYWFaNWqlX77+PHjkZqaatHBkW1IKUjl1cKNjTCJiBohNYdmaNd28G/pydVSFiL7zs0PP/yA77//Hvfcc4/B9vDwcPz+++8WGxjZjpQl4GrNLXazJSJqhNQcmkNnr+HQ2Wss6Gchsu/clJeXG9yxERUXF8PLy7z+Q2Q/tDoBR3OvSdqXWf1ERA1rbJn43dSlVZi5OR2pEruKk3Gyg5shQ4Zg06ZN+u8VCgV0Oh3eeustPPLIIxYdHDWv1MwiDF6+HykHciXtz6x+IqKGNbRM3BgW9LMM2dNSb731FoYPH45ffvkF1dXV+Mtf/oLffvsNxcXFOHr0qDXGSM1ATg8pZvUTEUknLhNP2pklOblYLOjHqX/zyA5uIiMjcfbsWaSkpMDX1xc3b97EE088gVmzZkGl4hyhI5LT0ZZLwImI5Ku7TPy7zCJsSms8R5VT/+aTFdzU1NQgLi4Oa9euxd/+9jdrjYma2U/npC1VBGrv2DDZjYhIvrrLxKUEN5z6N5+snBsPDw+cPHnSWmMhG0jNLMKsf6dL3n/hqAcY2BARNYGUJGOlnxen/ptAdkLxpEmTsG7dOmuMhZqZmGdTUlkj+Zilu7KZ5EZE1ARSkow1VbeRsj8HWp0ArU5AWt51fJNxEWl51/keLIHsnJvbt29j/fr12Lt3L/r27Vuvn9SKFSssNjiyHq1OwOId0vJs6mKSGxFR04lJxvO3n0JJRf0LzIpqLVbuzcE/D5+DZws3g31YC6dxsoObzMxMREVFAQDOnj1r8DP2lnIcKftzoNaYl6zGJDcioqaLjVBi8Y4sAKbvnldUa1FRbdjjT6yF88GkKAY4JkgKbk6ePInIyEi4ubnhwIHGG4CRfdLqBBzLL8aeLDXWHy0w+3GY5EZE1HTH8ovNusgUUDudlbQzC7ERSq5cNUJScNOnTx8UFRUhKCgI9913H/773/+ibVtOSziS1MwiyTUWTGF9GyIiy2nKXXDWwmmYpITigIAA5OfnAwAKCgqg0+msOiiyLDFxuKmBDcD6NkRElmKJu+B7stQWGInzkXTn5sknn8TQoUOhUqmgUCjQr18/uLu7G9333LlzFh0gNY2cAn0NYX0bIiLLEpeEN+XCc/3RAvTvFMj35rtICm4+/PBDPPHEE8jNzcWcOXMwY8YM+Pr6WntsZAFSOnwbEz88HA/d1xZXyqoQ5Fs7FcU7NkREliMuCX9xs/RaY3dj7o1xkldLxcXFAQCOHz+O+Ph4BjcOwpw5XaWfF+YMD+cLhYjIyuIiVVjbwJLwxjD3xjjZRfw2bNjAwMaByJ3TVQBYPKY7AxsiomYSF6nC8ddj8Vik0uzHYIkOQ7KDG3Is4pyuVK/EdOXcLRFRM3N3U+C56DCzj2eJDkMMblzA0w/eK3nfsHatrDgSIiIyRUrPKWPcFMCN8mqrjMlRMbhxYqmZRRi8fD9W7j3b+M7/w+ifiMg2pPScMkYnALM+TUdqZpF1BuaAGNw4Kbm1bRSo7VfCAn1ERLYj9pxS3pVO4N+yBRrqcCQAWLzjNzbV/B/ZvaXIPomtFa6UVaGdj5esppgs0EdEZD/iIlWIjVDq39ODfL2hEwQ88/HPDR6n1tzCrH8fx4NhgQhs7QWln+uW8WBw4wSa2lqBBfqIiOyLu5vCYGn3NxkXJR2X+ttlpP52Wf+9q3YQZ3Dj4MTpJ3NuRD4XHYrHIlUuG9kTETkKc/Mhi1y0gzhzbhxYU1srPBapQnTntgxsiIjsnNyyHndL2pnlUvk4DG4cmLmtFZg8TETkWOqupJKrbhVjV8HgxoGpNebl2Ahg8jARkaOJi1Rhbky42ce7UhVjBjcOKjWzCEu//c2sY6cNCnOpuVciImcxe1g4lH7mTU+5Uh0zBjcOaPfJIry4OR3F5fKbrAFAbIT5/UuIiMh23N0UWDwmQlaRP1dMRWBw42B2n7yE2VvSzT7e1f7AiYicjVjoT2qCsQBgfL978O3JS0jLu+4SicVcCu5AUjOL8NKnJ5r0GMy1ISJyfHUL/e3JUuPrjEsobqC/1Kp9ufr/d4XaNwpBEJw/hKtDo9HA398fpaWl8PPzs/VwJNPqBAxevt/sQn1uCiBlQhRG9nTeP2YiIlclVqn/4bcibPjx9wb3FS9vHa32jZzPb05L2TmtTkBa3nWs3HPG7MAGAFIm9GFgQ0TkpNzdFOjfKdCgOrEpwv++5n95CkdzrznlNBWnpexYU9sqAK5x+5GIiOTXPiuprMEzH//slJ8TDG7sVFPaKsyNCUdYOx8E+bpu0zQiIldjbh0btRO2aGBwY4fMbavAvBoiItdlbh0b8bPmr1+dQmWNzim6iTO4sTNanYCNR/PNmopiXg0RkesS+0+pS6vMuutfXF6DudsyAAABLT0wdVAYZg8Ld8gghwnFdkKrE7B6bw76Lt2DpbuyZR8/bVAYRvbsYIWRERGRI6jbf6qp4UhJZQ1W7s1B3zf2IDWzqOmDa2YMbuxAamYR+r6xByv3nkVJJasOExGRecQCf8omdBCvq6SiBi9uTsfuk5cs8njNhdNSNiDWI1CXVuJo7jV8kX7R7MdSAFCy6jAREf1P3QJ/32UWYVNaw3VvpJj16QnEX76Jl4c7xjQVgxsrEQOYK2VVaOfjBSiAazdvoeBaBbYcKzS7o3dd4p8Xqw4TEVFd7m4KRHduCwAWCW4EAKv25WBjWgHeHBeJNj5euFJWZberchncWMjdd2P2ZF9BqZlTTFIpnbA2ARERWU5Tk4zvVlJRU68NkL93C8RGBGNQeHsEtb5zMW/LwMcu2i+sWbMGb7/9NtRqNXr16oV//OMf6N+/v8n9P//8cyxcuBAFBQUIDw/H8uXLMXLkSEm/yxrtFyxRbE+O2Y90waAu7ewyWiYiIvsi1k0DYJEARw5LFgh0qPYL27ZtQ0JCAhITE5Geno5evXphxIgRuHLlitH9f/zxR0yYMAHTp0/HiRMnMG7cOIwbNw6ZmZnNPPJa4h9NcwQ2Ytv6ubFdEd25LQMbIiJqlKWTjOUQCwQ294orm9+5GTBgAB588EGkpKQAAHQ6HUJCQvDyyy9j/vz59fYfP348ysvL8e233+q3PfTQQ+jduzfWrl3b6O+z5J2bpjazlEsBx2t0RkRE9kFMn9iTpcb6owVQoHnu5IgLX47MG9aki3KHuXNTXV2N48ePIyYmRr/Nzc0NMTExSEtLM3pMWlqawf4AMGLECJP737p1CxqNxuDLUuT28WgKlb83AxsiIjKbmGS8aHR3rDVyJ8fH063J9XGMEQAUlVbhWH6xFR7dOJsmFF+7dg1arRbBwcEG24ODg3H69Gmjx6jVaqP7q9Vqo/svW7YMSUlJlhnwXczt4yFF3QQtZyiFTURE9qPucvG6q56+z1TjpU/TrfI7rfmZeTenXy21YMECJCQk6L/XaDQICQmxyGOb28ejMXNjwh225DURETmGusvFRSN7qrDWLQrzt59CSYVlV/xa6zPTGJsGN+3atYO7uzsuX75ssP3y5ctQKo1X3FUqlbL29/LygpeXl2UGfBdxiZ2lpqacse08ERE5FvGuTsr+XKw/cg6lVbeb9Hi2KDZr05wbT09P9O3bF/v27dNv0+l02LdvH6Kjo40eEx0dbbA/AOzZs8fk/tYk9vFo6v2VqQNDsWXGQzgybxgDGyIisjl3NwXiY8KRvuhRzI3pavbj2KrYrM2npRISEjB58mT069cP/fv3x6pVq1BeXo6pU6cCAJ577jl07NgRy5YtAwDEx8dj6NChePfddzFq1Chs3boVv/zyCz788EObjF9cYmeqzk1ASw9MHhgKBRTY+GOBQe8o3qkhIiJ7JgY53ZSt633OSVltZatiszYPbsaPH4+rV69i0aJFUKvV6N27N1JTU/VJw4WFhXBzu3ODaeDAgfj000/x+uuv469//SvCw8Px9ddfIzIy0lanYJCYpS6tRHF5NQJbe9VLBH55eHi95C3m1RARkb0zloDcN7QNjv9+w+BzjxWKbcQaFYqJiIjIuhymzg0RERGRpTG4ISIiIqfC4IaIiIicCoMbIiIicioMboiIiMipMLghIiIip8LghoiIiJwKgxsiIiJyKgxuiIiIyKnYvP1CcxMLMms0GhuPhIiIiKQSP7elNFZwueCmrKwMABASEmLjkRAREZFcZWVl8Pf3b3Afl+stpdPpcOnSJfj6+kKhsGwzL41Gg5CQEJw/f94p+1Y5+/kBzn+OPD/H5+znyPNzfNY6R0EQUFZWhg4dOhg01DbG5e7cuLm54Z577rHq7/Dz83PaP1rA+c8PcP5z5Pk5Pmc/R56f47PGOTZ2x0bEhGIiIiJyKgxuiIiIyKkwuLEgLy8vJCYmwsvLy9ZDsQpnPz/A+c+R5+f4nP0ceX6Ozx7O0eUSiomIiMi58c4NERERORUGN0RERORUGNwQERGRU2FwQ0RERE6FwU0D1qxZg7CwMHh7e2PAgAE4duxYg/t//vnnuP/+++Ht7Y0ePXpg9+7dBj8XBAGLFi2CSqVCy5YtERMTg5ycHGueQqPknONHH32EIUOGoE2bNmjTpg1iYmLq7T9lyhQoFAqDr7i4OGufhklyzm/jxo31xu7t7W2wj709h3LO7+GHH653fgqFAqNGjdLvY0/P3+HDhzF69Gh06NABCoUCX3/9daPHHDx4EFFRUfDy8kKXLl2wcePGevvIfV1bk9xz3L59O2JjY9G+fXv4+fkhOjoa33//vcE+ixcvrvcc3n///VY8C9Pknt/BgweN/o2q1WqD/Rz5OTT2GlMoFOjevbt+H3t5DpctW4YHH3wQvr6+CAoKwrhx43DmzJlGj7OHz0IGNyZs27YNCQkJSExMRHp6Onr16oURI0bgypUrRvf/8ccfMWHCBEyfPh0nTpzAuHHjMG7cOGRmZur3eeutt/Dee+9h7dq1+Pnnn+Hj44MRI0agqqqquU7LgNxzPHjwICZMmIADBw4gLS0NISEhePTRR3Hx4kWD/eLi4lBUVKT/2rJlS3OcTj1yzw+orahZd+y///67wc/t6TmUe37bt283OLfMzEy4u7vj//2//2ewn708f+Xl5ejVqxfWrFkjaf/8/HyMGjUKjzzyCDIyMvDKK6/g+eefN/jwN+dvwprknuPhw4cRGxuL3bt34/jx43jkkUcwevRonDhxwmC/7t27GzyHR44cscbwGyX3/ERnzpwxGH9QUJD+Z47+HK5evdrg3M6fP4/AwMB6r0N7eA4PHTqEWbNm4aeffsKePXtQU1ODRx99FOXl5SaPsZvPQoGM6t+/vzBr1iz991qtVujQoYOwbNkyo/s/9dRTwqhRowy2DRgwQPjTn/4kCIIg6HQ6QalUCm+//bb+5yUlJYKXl5ewZcsWK5xB4+Se491u374t+Pr6Cp988ol+2+TJk4WxY8daeqhmkXt+GzZsEPz9/U0+nr09h019/lauXCn4+voKN2/e1G+zp+evLgDCV1991eA+f/nLX4Tu3bsbbBs/frwwYsQI/fdN/TezJinnaExERISQlJSk/z4xMVHo1auX5QZmIVLO78CBAwIA4caNGyb3cbbn8KuvvhIUCoVQUFCg32avz+GVK1cEAMKhQ4dM7mMvn4W8c2NEdXU1jh8/jpiYGP02Nzc3xMTEIC0tzegxaWlpBvsDwIgRI/T75+fnQ61WG+zj7++PAQMGmHxMazLnHO9WUVGBmpoaBAYGGmw/ePAggoKC0K1bN8ycORPXr1+36NilMPf8bt68idDQUISEhGDs2LH47bff9D+zp+fQEs/funXr8PTTT8PHx8dguz08f+Zo7DVoiX8ze6PT6VBWVlbvNZiTk4MOHTrgvvvuwzPPPIPCwkIbjdA8vXv3hkqlQmxsLI4eParf7ozP4bp16xATE4PQ0FCD7fb4HJaWlgJAvb+3uuzls5DBjRHXrl2DVqtFcHCwwfbg4OB6c78itVrd4P7if+U8pjWZc453mzdvHjp06GDwRxoXF4dNmzZh3759WL58OQ4dOoTHHnsMWq3WouNvjDnn161bN6xfvx7ffPMNNm/eDJ1Oh4EDB+LChQsA7Os5bOrzd+zYMWRmZuL555832G4vz585TL0GNRoNKisrLfI3b2/eeecd3Lx5E0899ZR+24ABA7Bx40akpqbigw8+QH5+PoYMGYKysjIbjlQalUqFtWvX4ssvv8SXX36JkJAQPPzww0hPTwdgmfcte3Lp0iV899139V6H9vgc6nQ6vPLKKxg0aBAiIyNN7mcvn4Uu1xWcLCM5ORlbt27FwYMHDZJun376af3/9+jRAz179kTnzp1x8OBBDB8+3BZDlSw6OhrR0dH67wcOHIgHHngA//znP7F06VIbjszy1q1bhx49eqB///4G2x35+XM1n376KZKSkvDNN98Y5KQ89thj+v/v2bMnBgwYgNDQUHz22WeYPn26LYYqWbdu3dCtWzf99wMHDkReXh5WrlyJf/3rXzYcmXV88sknCAgIwLhx4wy22+NzOGvWLGRmZtosf0su3rkxol27dnB3d8fly5cNtl++fBlKpdLoMUqlssH9xf/KeUxrMuccRe+88w6Sk5Pxww8/oGfPng3ue99996Fdu3bIzc1t8pjlaMr5iTw8PNCnTx/92O3pOWzK+ZWXl2Pr1q2S3iRt9fyZw9Rr0M/PDy1btrTI34S92Lp1K55//nl89tln9aYA7hYQEICuXbs6xHNoTP/+/fVjd6bnUBAErF+/Hs8++yw8PT0b3NfWz+Hs2bPx7bff4sCBA7jnnnsa3NdePgsZ3Bjh6emJvn37Yt++ffptOp0O+/btM7iyrys6OtpgfwDYs2ePfv9OnTpBqVQa7KPRaPDzzz+bfExrMuccgdos96VLlyI1NRX9+vVr9PdcuHAB169fh0qlssi4pTL3/OrSarU4deqUfuz29Bw25fw+//xz3Lp1C5MmTWr099jq+TNHY69BS/xN2IMtW7Zg6tSp2LJli8EyflNu3ryJvLw8h3gOjcnIyNCP3VmeQ6B2JVJubq6kiwxbPYeCIGD27Nn46quvsH//fnTq1KnRY+zms9BiqclOZuvWrYKXl5ewceNGISsrS3jhhReEgIAAQa1WC4IgCM8++6wwf/58/f5Hjx4VWrRoIbzzzjtCdna2kJiYKHh4eAinTp3S75OcnCwEBAQI33zzjXDy5Elh7NixQqdOnYTKyspmPz9BkH+OycnJgqenp/DFF18IRUVF+q+ysjJBEAShrKxMeO2114S0tDQhPz9f2Lt3rxAVFSWEh4cLVVVVdn9+SUlJwvfffy/k5eUJx48fF55++mnB29tb+O233/T72NNzKPf8RIMHDxbGjx9fb7u9PX9lZWXCiRMnhBMnTggAhBUrVggnTpwQfv/9d0EQBGH+/PnCs88+q9//3LlzQqtWrYQ///nPQnZ2trBmzRrB3d1dSE1N1e/T2L9Zc5N7jv/+97+FFi1aCGvWrDF4DZaUlOj3efXVV4WDBw8K+fn5wtGjR4WYmBihXbt2wpUrV+z+/FauXCl8/fXXQk5OjnDq1CkhPj5ecHNzE/bu3avfx9GfQ9GkSZOEAQMGGH1Me3kOZ86cKfj7+wsHDx40+HurqKjQ72Ovn4UMbhrwj3/8Q7j33nsFT09PoX///sJPP/2k/9nQoUOFyZMnG+z/2WefCV27dhU8PT2F7t27C7t27TL4uU6nExYuXCgEBwcLXl5ewvDhw4UzZ840x6mYJOccQ0NDBQD1vhITEwVBEISKigrh0UcfFdq3by94eHgIoaGhwowZM2z2piMI8s7vlVde0e8bHBwsjBw5UkhPTzd4PHt7DuX+jZ4+fVoAIPzwww/1Hsvenj9xWfDdX+I5TZ48WRg6dGi9Y3r37i14enoK9913n7Bhw4Z6j9vQv1lzk3uOQ4cObXB/Qahd/q5SqQRPT0+hY8eOwvjx44Xc3NzmPbH/kXt+y5cvFzp37ix4e3sLgYGBwsMPPyzs37+/3uM68nMoCLVLn1u2bCl8+OGHRh/TXp5DY+cFwOB1Za+fhYr/nQARERGRU2DODRERETkVBjdERETkVBjcEBERkVNhcENEREROhcENERERORUGN0RERORUGNwQERGRU2FwQ0RERE6FwQ0ROaTFixejd+/eso5RKBT4+uuvrTIeIrIfDG6IyK48/PDDeOWVVxrd77XXXqvXoI+ICABa2HoARERyCIIArVaL1q1bo3Xr1rYeDhHZId65ISK7MWXKFBw6dAirV6+GQqGAQqHAxo0boVAo8N1336Fv377w8vLCkSNH6k1L/fe//0VsbCzatWsHf39/DB06FOnp6SZ/V3V1NWbPng2VSgVvb2+EhoZi2bJlzXCWRGRtDG6IyG6sXr0a0dHRmDFjBoqKilBUVISQkBAAwPz585GcnIzs7Gz07Nmz3rFlZWWYPHkyjhw5gp9++gnh4eEYOXIkysrKjP6u9957Dzt27MBnn32GM2fO4N///jfCwsKseXpE1Ew4LUVEdsPf3x+enp5o1aoVlEolAOD06dMAgCVLliA2NtbkscOGDTP4/sMPP0RAQAAOHTqE//u//6u3f2FhIcLDwzF48GAoFAqEhoZa8EyIyJZ454aIHEK/fv0a/Pnly5cxY8YMhIeHw9/fH35+frh58yYKCwuN7j9lyhRkZGSgW7dumDNnDn744QdrDJuIbIDBDRE5BB8fnwZ/PnnyZGRkZGD16tX48ccfkZGRgbZt26K6utro/lFRUcjPz8fSpUtRWVmJp556Cn/84x+tMXQiamacliIiu+Lp6QmtViv7uKNHj+L999/HyJEjAQDnz5/HtWvXGjzGz88P48ePx/jx4/HHP/4RcXFxKC4uRmBgoFljJyL7wOCGiOxKWFgYfv75ZxQUFKB169bQ6XSSjgsPD8e//vUv9OvXDxqNBn/+85/RsmVLk/uvWLECKpUKffr0gZubGz7//HMolUoEBARY6EyIyFY4LUVEduW1116Du7s7IiIi0L59e5M5M3dbt24dbty4gaioKDz77LOYM2cOgoKCTO7v6+uLt956C/369cODDz6IgoIC7N69G25ufFskcnQKQRAEWw+CiIiIyFJ4iUJEREROhcENERERORUGN0RERORUGNwQERGRU2FwQ0RERE6FwQ0RERE5FQY3RERE5FQY3BAREZFTYXBDREREToXBDRERETkVBjdERETkVP4/yt6pXhCReNwAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import sympy\n", "\n", "# Perform an X gate with variable exponent\n", "q = cirq.GridQubit(1, 1)\n", "circuit = cirq.Circuit(cirq.X(q) ** sympy.Symbol('t'), cirq.measure(q, key='m'))\n", "\n", "# Sweep exponent from zero (off) to one (on) and back to two (off)\n", "param_sweep = cirq.Linspace('t', start=0, stop=2, length=200)\n", "\n", "# Simulate the sweep\n", "s = cirq.Simulator()\n", "trials = s.run_sweep(circuit, param_sweep, repetitions=1000)\n", "\n", "# Plot all the results\n", "x_data = [trial.params['t'] for trial in trials]\n", "y_data = [trial.histogram(key='m')[1] / 1000.0 for trial in trials]\n", "plt.scatter('t', 'p', data={'t': x_data, 'p': y_data})\n", "plt.xlabel(\"trials\")\n", "plt.ylabel(\"frequency of qubit measured to be one\")\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": { "id": "M8oLYwusz4XE" }, "source": [ "## Unitary matrices and decompositions\n", "\n", "Many quantum operations have unitary matrix representations. This matrix can be accessed by applying `cirq.unitary(operation)` to that `operation`. This can be applied to gates, operations, and circuits that support this protocol and will return the unitary matrix that represents the object. See [Protocols](../build/protocols.ipynb) for more about this and other protocols." ] }, { "cell_type": "code", "execution_count": 16, "metadata": { "execution": { "iopub.execute_input": "2023-05-18T09:44:58.099531Z", "iopub.status.busy": "2023-05-18T09:44:58.098880Z", "iopub.status.idle": "2023-05-18T09:44:58.105986Z", "shell.execute_reply": "2023-05-18T09:44:58.105276Z" }, "id": "xn9nnBA70s23" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Unitary of the X gate\n", "[[0.+0.j 1.+0.j]\n", " [1.+0.j 0.+0.j]]\n", "Unitary of SWAP operator on two qubits.\n", "[[1.+0.j 0.+0.j 0.+0.j 0.+0.j]\n", " [0.+0.j 0.+0.j 1.+0.j 0.+0.j]\n", " [0.+0.j 1.+0.j 0.+0.j 0.+0.j]\n", " [0.+0.j 0.+0.j 0.+0.j 1.+0.j]]\n", "Unitary of a sample circuit\n", "[[0.+0.j 0.+0.j 1.+0.j 0.+0.j]\n", " [1.+0.j 0.+0.j 0.+0.j 0.+0.j]\n", " [0.+0.j 0.+0.j 0.+0.j 1.+0.j]\n", " [0.+0.j 1.+0.j 0.+0.j 0.+0.j]]\n" ] } ], "source": [ "print('Unitary of the X gate')\n", "print(cirq.unitary(cirq.X))\n", "\n", "print('Unitary of SWAP operator on two qubits.')\n", "q0, q1 = cirq.LineQubit.range(2)\n", "print(cirq.unitary(cirq.SWAP(q0, q1)))\n", "\n", "print('Unitary of a sample circuit')\n", "print(cirq.unitary(cirq.Circuit(cirq.X(q0), cirq.SWAP(q0, q1))))" ] }, { "cell_type": "markdown", "metadata": { "id": "Ls6Tnx90Y94Q" }, "source": [ "### Decompositions \n", "\n", "Many gates can be decomposed into an equivalent circuit with simpler operations and gates. This is called decomposition and can be accomplished with the `cirq.decompose` protocol. \n", "\n", "For instance, a Hadamard H gate can be decomposed into X and Y gates:" ] }, { "cell_type": "code", "execution_count": 17, "metadata": { "execution": { "iopub.execute_input": "2023-05-18T09:44:58.109361Z", "iopub.status.busy": "2023-05-18T09:44:58.108803Z", "iopub.status.idle": "2023-05-18T09:44:58.113126Z", "shell.execute_reply": "2023-05-18T09:44:58.112466Z" }, "id": "u8JwZaAUfbSv" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[(cirq.Y**0.5).on(cirq.LineQubit(0)), cirq.XPowGate(global_shift=-0.25).on(cirq.LineQubit(0))]\n" ] } ], "source": [ "print(cirq.decompose(cirq.H(cirq.LineQubit(0))))" ] }, { "cell_type": "markdown", "metadata": { "id": "B8ciZZSSf2jb" }, "source": [ "Another example is the 3-qubit Toffoli gate, which is equivalent to a controlled-controlled-X gate. Many devices do not support three qubit gates, so it is important to decompose them into one and two qubit gates." ] }, { "cell_type": "code", "execution_count": 18, "metadata": { "execution": { "iopub.execute_input": "2023-05-18T09:44:58.116433Z", "iopub.status.busy": "2023-05-18T09:44:58.115873Z", "iopub.status.idle": "2023-05-18T09:44:58.125887Z", "shell.execute_reply": "2023-05-18T09:44:58.125189Z" }, "id": "bbjRWlzjgPwf" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0: ───T────────────────@─────────────────────────────────@─────────────────────────────@────────────────────────────@───────────────────────────────────────\n", " │ │ │ │\n", "1: ───T───────Y^-0.5───@───Y^0.5────@───T^-1────Y^-0.5───@────────Y^0.5───@───Y^-0.5───@──────Y^0.5────@───Y^-0.5───@──────Y^0.5────@───────────────────────\n", " │ │ │ │\n", "2: ───Y^0.5───X────────T───Y^-0.5───@───Y^0.5───T────────Y^-0.5───────────@───Y^0.5────T^-1───Y^-0.5───@───Y^0.5────T^-1───Y^-0.5───@───Y^0.5───Y^0.5───X───\n" ] } ], "source": [ "q0, q1, q2 = cirq.LineQubit.range(3)\n", "print(cirq.Circuit(cirq.decompose(cirq.TOFFOLI(q0, q1, q2))))" ] }, { "cell_type": "markdown", "metadata": { "id": "VWcik4ZwggXj" }, "source": [ "The above decomposes the Toffoli into a simpler set of one-qubit gates and two-qubit CZ gates at the cost of lengthening the circuit considerably." ] }, { "cell_type": "markdown", "metadata": { "id": "rIUbvdVQkHbX" }, "source": [ "### Transformers\n", "\n", "The last concept in this tutorial is the transformer. A transformer can take a circuit and modify it. Usually, this will entail combining or modifying operations to make it more efficient and shorter, though a transformer can, in theory, do any sort of circuit manipulation.\n", "\n", "For example, the `cirq.merge_single_qubit_gates_to_phxz` transformer will take consecutive single-qubit operations and merge them into a single `PhasedXZ` operation." ] }, { "cell_type": "code", "execution_count": 19, "metadata": { "execution": { "iopub.execute_input": "2023-05-18T09:44:58.129299Z", "iopub.status.busy": "2023-05-18T09:44:58.128745Z", "iopub.status.idle": "2023-05-18T09:44:58.136894Z", "shell.execute_reply": "2023-05-18T09:44:58.136244Z" }, "id": "5WvfOdaG5C_6" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "(1, 1): ───X^0.25───Y^0.25───T───\n", "(1, 1): ───PhXZ(a=0.304,x=0.333,z=0.142)───\n" ] } ], "source": [ "q = cirq.GridQubit(1, 1)\n", "c = cirq.Circuit(cirq.X(q) ** 0.25, cirq.Y(q) ** 0.25, cirq.Z(q) ** 0.25)\n", "print(c)\n", "c = cirq.merge_single_qubit_gates_to_phxz(c)\n", "print(c)" ] }, { "cell_type": "markdown", "metadata": { "id": "xRfQqzdx7lUI" }, "source": [ "Other transformers can assist in transforming a circuit into operations that are native operations on specific hardware devices. You can find more about transformers and how to create your own in [Transformers](../transform/transformers.ipynb)." ] }, { "cell_type": "markdown", "metadata": { "id": "8QbTGmKlYT4i" }, "source": [ "## Next steps\n", "\n", "After completing this tutorial, you should be able to use gates and operations to construct your own quantum circuits, to simulate them, and to use sweeps.\n", "\n", "There is much more to learn and try out for those who are interested:\n", "\n", "* Learn about the variety of [Gates](../build/gates.ipynb) available in cirq and more about the different ways to construct [Circuits](../build/circuits.ipynb).\n", "* Learn more about [Simulations](../simulate/simulation.ipynb) and how it works.\n", "* Learn about [Noise](../noise) and how to utilize multi-level systems using [Qudits](../build/qudits.ipynb).\n", "* Dive into some [Experiments](../experiments) and some in-depth tutorials of how to use cirq.\n", "\n", "Also, join our [cirq-announce mailing list](https://groups.google.com/forum/#!forum/cirq-announce) to hear about changes and releases or go to the [Cirq GitHub repository](https://github.com/quantumlib/Cirq/) to file issues." ] } ], "metadata": { "colab": { "collapsed_sections": [], "name": "basics.ipynb", "toc_visible": true }, "kernelspec": { "display_name": "Python 3", "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.9.16" } }, "nbformat": 4, "nbformat_minor": 0 }