2.d-Multi-qubit-circuits.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Muti Qubits Circuits"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "from qiskit import *\n",
    "from math import pi\n",
    "import numpy as np\n",
    "from qiskit.visualization import *\n",
    "import matplotlib.pyplot as plt\n",
    "from qutip import *"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Controlled Phase Rotation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMoAAAB7CAYAAAA18x94AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAC1JJREFUeJzt3X9M1Heex/HnDCgMUA9npgqL2pUwgh2BVXab9lj1zEmDmyJpz3raaBq8Lf5a/6i1Zy6k7SVtkVA2TRP3uubuNpp4atTKiZWusrt1Uq/turSu7pTjxq6uiLIqRYGh1MIw9wfnrFO6x2dxfhTzeiST6Gf4Mu9JePKZHzpfSzAYDCIi/y9rvAcQGQ8UiogBhSJiQKGIGFAoIgYUiogBhSJiQKGIGFAoIgYUiogBhSJiQKGIGFAoIgYUiogBhSJiQKGIGFAoIgYUiogBhSJiQKGIGFAoIgYUiogBhSJiQKGIGFAoIgYUiogBhSJiQKGIGFAoIgYUiogBhSJiQKGIGFAoIgYUiogBhSJiQKFEWHd3N+vWrWPKlCmkpKRQXFzMyZMn4z2W3CWFEkHBYJDy8nLq6+upq6ujoaEBp9NJSUkJp0+fjvd4URH4Em71wVAg3pNEl0VnBY6cI0eOsHTpUhobG1myZAkAX375JW63G5fLRWNjY5wnjJwb7XDxFHSeH/67NREy3fDt74EtPb6zRYN2FENDQ0PU1dXhcrlITk6msLAQj8dDbm4ulZWVABw+fBiHw0FpaWnouIkTJ7JixQqampro6+uL1/gR1fEJfLQPOi/8aW1oEC6fgV/vBv/1+M0WLQrF0Jo1a3j55ZdZu3Yt77zzDsuXL2flypWcP3+eoqIiALxeL263G4vFEnbsnDlzGBwcpLW1NR6jR1R/N7T8/P/+8jWPRQZvwdkGuNcepygUA3v27GHXrl00NDSwZcsWFi1aRFVVFY888giDg4OhULq6upg8efKI4+12e+h6gKtXr/Loo4+SkpJCYWHhuHr+cvnsKBEE4fMbcKMtZiPFhEIxsG3bNkpLS1m4cGHYek5ODhMmTCA/Px8YfjL/1d0EGLG2fv168vLy+Oyzz9i4cSPLli0jEIjus2GLxRKRyzv7P2S0p7XBYJDn1r4SsduM5sWUQhlFe3s7Xq+XJ598csR1bW1tuN1ukpKSAHA4HKFd40631+x2O729vRw9epQXX3wRm81GZWUlgUCADz/8MLp3JEKs1gSDH7AgVmtCTOaJFYUyivb2dgAyMjLC1vv7+/F4PKGHXQBut5uWlpYRv3G9Xi+JiYnk5eVx7tw5HA4HTqczdH1+fj4tLS1RvBfDv+Ujcfnbx7436m1ZLFaqX/+niN1mNC+mFMoobv9A+3y+sPXa2lo6OjqYN29eaK28vJzOzk6OHTsWWhsYGGDfvn0sXryY1NRU+vr6mDRpUtj3mjRpEn6/P4r3InKmfWeUL7DAxBS4Pycm48RMYrwH+KbLzs6moKCA6upq7HY7WVlZHDx4MPSeyJ07SllZGfPnz6eiooLa2loyMzPZvn07bW1t7N27F4DU1FR6e3vDbqOnp4e0tLTY3am7kHY/zHwELnzwNVdawGIB9w/gHnvkpTccTfh8PtauXcupU6dwOBw8/fTT3HfffVRVVdHT04PNZgt97c2bN9m6dSuHDh3C7/czd+5campqWLBgAQC9vb04nU6uXLmCw+EAYObMmezevZvi4uK43L+xuPI7uPBr6L/5p7X06ZDzfUjPit9c0aJQxmj16tWcOXOGs2fP/sXHPvHEE8yYMYOamhp2795NdXU1586dIyFhfP0aDgbhlz8e/vNf/xBS7sF35G/TQ68xam5u5uGHHx7TsW+++SarVq1i8uTJuFwu3nrrrXEXCQw/zLrtXo4EFMqY+P1+fD4fGzZsGNPxU6dOpampKcJTSTQplDFIS0uL+huE8s2il4dFDCgUEQMKRcSAQhExoFBEDCgUEQMKRcSAQhExoFBEDCgUEQMKRcSAQhExoFBEDCgUEQMKRcSAQhExoFBEDCgUEQMKRcSAQhExoFBEDCgUEQP6pEgJ+UVd7G9z8ZbY3+ZYaEcRMaBQRAwoFBEDCkWM3fBfo+R5C59e/m3Y+tUbFyl53sLFq/8dp8miT6GIMd+l3zAxMZmZGXPC1lvbTpGSdB/T78+N02TRp1DEmO9SMzlZc0lICP9s9/+5dArXtCKs1nv3x+nevWcScb72ZmZNH3my09ZLp8idNvpJUMczhRJh3d3drFu3jilTppCSkkJxcTEnT56M91gR4bv80YgghoaG+PTyx6GA/nHHYv7uJSf/8YtX4jFi1CiUCAoGg5SXl1NfX09dXR0NDQ04nU5KSko4ffp0vMe7K53dl+nq6cA1rShs/fdXfkv/LT+zHxg++9jzf7+Tysfi8M5llOlEQhH09ttv4/F4aGxsZMmSJQAsWLAAt9tNVVVV6EzC49H17nYA0mzh56D71ek9zJ7xMFPSpwNwf/q0mM8WC9pRDA0NDVFXV4fL5SI5OZnCwkI8Hg+5ublUVlYCcPjwYRwOB6WlpaHjJk6cyIoVK2hqaqKvry9e49+1GVNmk5r8V+x7twZ//016Pu/iyAc/5cgH/8I//KAm3uNFnXYUQ2vWrKG+vp4XXniBoqIi3n//fVauXMn169fZvHkzAF6vF7fbjeXOs4ACc+bMYXBwkNbW1rDz0o8nqcmTeGXN2/zr0a2sfGUayRNSyP7Wd6j+4c8pyF4Q7/GiTqEY2LNnD7t27eLEiRMsXLgQgEWLFvHxxx9z6NCh0A9/V1cXDz744Ijj7XZ76HqAl156iQMHDtDa2sr+/ftZtmxZjO7J3Zkz8/u88aP/ivcYcaGHXga2bdtGaWlpKJLbcnJymDBhAvn5+cDwk/mv7ibAiDWXy8Ubb7zBQw89FL2hv2aG0S6RULd/DQc8r3G8eScv/KwsInNF82JKO8oo2tvb8Xq9PPvssyOua2trw+12k5SUBIDD4QjtGne6vXZ7Z1m1ahUAr776arTGjpsty38W7xGiQjvKKNrbh1/tycjICFvv7+/H4/GEPedwu920tLTw1f/i4/V6SUxMJC8vL/oD/xnBYHDUyzd1rmheTCmUUTidTgB8Pl/Yem1tLR0dHcybNy+0Vl5eTmdnJ8eOHQutDQwMsG/fPhYvXkxqampshpaI00OvUWRnZ1NQUEB1dTV2u52srCwOHjwYek/kzh2lrKyM+fPnU1FRQW1tLZmZmWzfvp22tjb27t0br7sgEaAdZRRWq5UDBw7gdrtZv349FRUVOJ1ONm7cSGJiIgUFBaGvtVgsNDQ0sHTpUjZv3kxZWRnXrl3j+PHj4/ZlYRmmHcXArFmzePfdd8PWVq9ezezZs7HZbGHr6enp7Nixgx07dvzZ7zcwMEAgEGBoaIiBgQG++OILkpKSIvbKUyx1dl9m/4nXWDR3JT9teBarNYFZ077L+qWvx3u0iNKOMkbNzc1j3iWeeeYZbDYb7733Hk899RQ2m42LFy9GeMLY+MjXRNGsEqamP8Bra3/F6xve46b/Ghc6fhfv0SJKO8oY+P1+fD4fGzZsGNPxO3fuZOfOnZEdKgbO/P4E/7zrcbIzC/lj1wVysuaSmpzOpse3Y0tKC31dgjURqzUhjpNGnnaUMUhLSyMQCLBp06Z4jxJT+TMXkDv9IX68/gQF2QvZ9PhPuDXweVgk56+cpbuvkwemjvwXCuOZQhFjHV3nybRnA3C9+xI3/dfJ/lZh6Pqez7vY/p8/4rkn/z1eI0aNQhFjF//4CQ9kuAkMBbBYrHzkO06RqwSAQGCQmr2reOax17BPyhjlO40/CkWM/eHqJ3x7qpuBwVvc9F/j9Ke/ZNa07wLgOXsA36Xf8G9Ht/Lcm39Dyx8+iPO0kaWPVJWQv/QjVT1n9rOwcPld3aY+UlXueXcbyXiiHUXEgHYUEQMKRcSAQhExoFBEDCgUEQMKRcSAQhExoFBEDCgUEQMKRcSAQhExoFBEDCgUEQMKRcSAQhExoFBEDCgUEQMKRcSAQhEx8L+Xr248o0/fVwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 249.228x144.48 with 1 Axes>"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "q = QuantumRegister(2)\n",
    "qc = QuantumCircuit(q)\n",
    "qc.cu1(pi/2,q[0], q[1])\n",
    "qc.draw(output='mpl')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[[1.+0.j 0.+0.j 0.+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 1.+0.j 0.+0.j]\n",
      " [0.+0.j 0.+0.j 0.+0.j 0.+1.j]]\n"
     ]
    }
   ],
   "source": [
    "backend = Aer.get_backend('unitary_simulator')\n",
    "job = execute(qc, backend)\n",
    "result = job.result()\n",
    "print(result.get_unitary(qc, decimals=3))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "-------------------"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### U operator and control phase gate"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "$ (U_3 \\otimes I \\otimes I)  \\times (CU_1 \\otimes I ) \\times (I \\otimes U_3 \\otimes I)  \\times (I \\otimes CU_1) \\times (I \\otimes I \\otimes U)$"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAggAAACoCAYAAACFfxakAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAG7NJREFUeJzt3WtUU3e6BvAnCZBwLQJCEDCCXFWgElutjlKPupaXWa5xWrVwyrFjWxltdY1Tna7RtrrGehmlndOedqaXc7pw2SoqdXpR6xQUaDvWIt6p0hS5GeQSRQLhJiQ5H2xxMhshYMgO4fl96fLNzrtfU0ge9/5nb4nZbDaDiIiI6F9IxR6AiIiIHA8DAhEREQkwIBAREZEAAwIREREJMCAQERGRAAMCERERCTAgEBERkQADAhEREQkwIBAREZEAAwIREREJMCAQERGRAAMCERERCTAgEBERkQADAhEREQkwIBAREZEAAwIREREJMCAQERGRAAMCERERCTAgEBERkQADAhEREQkwIBAREZEAAwIREREJMCAQERGRAAMCERERCTAgEBERkQADAhEREQkwIBAREZGAi9gD0NDRuX6f2CPYheuuFLFHICISHY8gEBERkQADAhEREQkwIBAREZEAAwIREREJMCCQQ6hv0cMtIxXn6yss6pV6HdwyUnHlZrU4gxERDVMMCOQQimrLoHBxxYSAMIv66dqr8HZzR4xfsEiTERENTwwI5BDO1JXhwcAxcJHKLOpFNVeRFBQOqYQ/qkRE9sR3XRvSarVYs2YNpk6dCg8PD0gkEhQXF4s91pBwprYMk5QRgvrp2qtQ91AnIqLBxYBgQ6Wlpdi/fz98fX2RnJws9jhDyrm6cqiVYy1qJrMJ5+oqMEkZgRutTZj+0SuYlfUnPLLnJeRWXBJpUnE01QFlJ4EfC4Dqi0BXh9gTEZGz45UUbWjGjBmoq6sDAGRmZuLYsWMiTzQ0VDc3oKalEUlB4Rb1C/WVMHS2Y3JwFEYovJCfshkyqRSlt2qRdvgtzB4TL9LE9nO7Bbh0GLh1zbL+wwkgKhkImyjOXETk/HgEwUomkwkZGRmIioqCQqFAYmIiCgoKEBMTgxUrVgAApFK+nANRbWgAAPjKPSzqWVdOYnJwJMJ8/CGTSiH76fVtvt2OhMDRdp/T3oydwJmDwnAAAKYu4Ifjd44mEBENBn6iWWn58uXYsmUL0tPT8cUXX2DJkiVISUlBWVkZ1Gq12OMNabF+o/CA3AO7Cj9DY3sLGtoMeO98Lt45n4OtM+7eF6G8sR6P7tuMBdnbsTBykogT20ftFaDlRu/b/PjVnbBARGRrPMVghb1792L37t3Iz8/vXlswc+ZMnD17FocOHep3QKirq0NaWhq++eYbREVFITMzExMnDt9jxT5yD3yyaB02fJWF8Hefh4eLHImBKhx+7EVMD4vr3i7cNxD5KZtRoddhzv4tWDA2ScSpB1/1RQASAOZ7b9PVDujKgKBoe01FRMMFA4IVtm/fjrlz5woWHkZGRsLV1RXx8f07F75y5UrExsbi008/xZ49e/D4449Do9FAJpP1/eQBkkgk993j9rq9NpikZ9NCY1GQuvmej3d0dULu4goA8HFzh5eb+6DNYovXyhYOvFKDEd7KPrdb9cxaHPr6v+0wERE5A7O5l391/AsGhD5otVoUFxdj7dq1gseqqqowfvx4yOVyq/s1NzfjyJEjqK6uhru7O1asWIFt27bh1KlTmDZtmi1Hdyrn6sux4assyCRSdJmMeG1mmtgjDTpDmx6+XoGQ9HENiJZ2vZ0mIqLhhAGhD1qtFgCgVFr+S66trQ0FBQWYP39+v/r9+OOP8Pf3R0BAQHctPj4ely9fHtSAYG1i7E3n+n02mGRgpoyKxoknXrHLvmzxWtnC1X8C5d/2vo1ECnyW/wHcPD6wz1BENGxwkWIffv4g12g0FvWdO3eipqYGSUn9Ow/e0tICHx8fi5qPjw8MBsP9DUpOJzQRkLnhzjqEewhJBNw87v04EdFA8QhCHyIiIpCQkIBt27bBz88PISEhyM7OxtGjRwFAsEAxOzsbAFBUVAQAyMnJQUlJCTw9PTFv3jx4enqiubnZ4jlNTU3w8vKyw9+GhhK5F5D0OHDu43+7MNJPCxdHRgHRj4o0HBE5PYnZUY6nOjCNRoP09HQUFhbC398fy5Ytg7e3NzZu3Iimpia4u99dMHevBW4qlQoVFRVobm5GQEAArl+/Dn9/fwBAeHg4PvzwQ4dfgyDmKQZ7ct2V0vdGdtTZDlwvBn7Mv/PnoFggJAEYEQY4yHpKInJCPIJghejoaOTl5VnU0tLSEBcXZxEOgL7PX3t7e2PBggXYsmULduzYgQ8//BASiQRTpkyx+dxiqm5uwOunD2Np3FSsy9sDmUQKtTICGf1cXGjLPr/9x/vQ324dcA+xuCoA1aS7ASH+l6KOQ0TDBNcgDFBRUdGAL5D0t7/9Dd9//z1GjBiBN998Ex9//PGgfsVRDMcrL2GWKh6jfQLw5ZKNyEvZhPpWPS7pqkTrkzJu2n31ICIaTngEYQAMBgM0Gg1WrVo1oOcHBQUhJyfHxlOJp6DqMhZ/+joSRqpQodchMUgFX7kn3pj1FLzcFN3buUhkkN3jK3s99fj4Vy+g4NoVm/ZRuLj12YOIiHgEYUC8vLxgNBqxevVqsUdxCNPDYvFQ8FjkPvEypofF4s1Zv0FrZ4fFh/pFXRVutjVjXECo1T3MZvOg9OmrBxERMSCQDZQ11iP8gUAAgLa5Abq2JiQEqrofb2gz4HfHM/Hu3BVW9wjx9sMFXaXN+1jTg4iIGBDIBi7f1GKcfyiMJhOkEglyKy5htmoCAKDLZMRTR9/GjuRUKD19u2t1LfpeewCweZ+eehARUc8YEOi+Xb6hxbiAUHQYO1Hf2oS8qmKolREAgOwfvkNRbRk2FOzD7KwtOHVdgwq9Dpu+OdBrjxrDLZytK7dpn556EBFRz3gdBLKatddBOFhyCotj7/21zUOaQoxQeGLm6PEO0effOdp1EH6Wm3Hnv7PXiTsHEQ0PDAhkNV4oSVwMCERkT/yaI1nNUT84iYjI9rgGgYiIiAQYEIiIiEiAAYGIiIgEGBCIiIhIgAGBiIiIBBgQiIiISIABgYiIiAQYEIiIiEiAAYGIiIgEGBCIiIhIgAGBiIiIBBgQiIiISIABgYiIiAQYEIiIiEiAAYGIiIgEJGaz2Sz2EETDTW6G/fc5e93g76Nz/b7B34nIXHeliD0CkV3wCAIREREJMCAQERGRAAMCERERCTAgEBERkQADApEDumWox5z1EpRWn7eo192qxJz1ElTWXRFpMnHUt+jhlpGK8/UVFvVKvQ5uGam4crNanMGInBgDApED0lw7DTcXBcKVEyzqJVWF8JB7I2xkjEiTiaOotgwKF1dMCAizqJ+uvQpvN3fE+AWLNBmR82JAIHJAmmtFiAyZCJnMxaL+w7VCRIWqIZUOr1/dM3VleDBwDFykMot6Uc1VJAWFQyoZXq8HkT3wt4rIAWm0RYgOe0hQL7lWiJhQYd3ZnaktwyRlhKB+uvYq1D3Uiej+MSDYkFarxZo1azB16lR4eHhAIpGguLhY7LFoCNJUnxEEAZPJhNLqs93B4Q/vzsZjmwLwUe6rYoxoV+fqyqFWjrWomcwmnKurwCRlBG60NmH6R69gVtaf8Miel5BbcUmkScVjMgIdLYDxttiTkLNw6XsTslZpaSn2798PtVqN5ORkHDt2TOyRaAi6oa9GQ1MNokLVFvWr18+jrcOAONUUAMD6pZk4+2Mubui1YoxpN9XNDahpaURSULhF/UJ9JQyd7ZgcHIURCi/kp2yGTCpF6a1apB1+C7PHxIs0sX21NQIVp4Ga7wFTFwAJEBAOqB4GRoSKPR0NZQwINjRjxgzU1dUBADIzMxkQaEB0P33ge7n7WtRPnNuLuNFTEOh7Z6HeSN/h8e5fbWgAAPjKPSzqWVdOYnJwJMJ8/C3qzbfbkRA42m7zicmgA4qygK6OfymagRvlwI0yYPx8IHicaOPREMdTDFYymUzIyMhAVFQUFAoFEhMTUVBQgJiYGKxYsQIAht3CMRocowPj4Kl4AFl5O2Boa0RTawM+//YdfP7tX/H0/B1ij2d3sX6j8IDcA7sKP0Njewsa2gx473wu3jmfg60z7t4XobyxHo/u24wF2duxMHKSiBPbh9kMXPgU6OrplMJPd9i5/AXQprfrWOREeATBSsuXL8ff//53vPzyy1Cr1Th58iRSUlKg0+nw+9//XuzxyIl4Knzw6vLDeP/Ii0h5NRQKVw9EjHoQ2545hoSIGWKPZ3c+cg98smgdNnyVhfB3n4eHixyJgSocfuxFTA+L694u3DcQ+SmbUaHXYc7+LVgwNknEqQdfQ+Wd0wu9MZuB6otA5HT7zETOhQHBCnv37sXu3buRn5+P5ORkAMDMmTNx9uxZHDp0CGq1uo8OljZt2oSDBw+ipKQEBw4cwOOPPz4YY1uQSCSDvg+yXs6u3m+iOiH8F3jj+X/adJ/2+Bm4vW7voPSdFhqLgtTN93y8o6sTchdXAICPmzu83NwHZQ7AcX6Xls/diif+44+9zmM2m3Fk/7f43YxpdpyMHJ21N3FmQLDC9u3bMXfu3O5w8LPIyEi4uroiPr5/i6GioqLwxhtv4OWXX7blmDTMZBxYjpKq79DZ1YGSqu+wZfnnYo8kmnP15djwVRZkEim6TEa8NjNN7JEGnVQqgxlmSHDvgCCRSCD7t2tHEFmLAaEPWq0WxcXFWLt2reCxqqoqjB8/HnK5vF89n3zySQDA1q1bbTKjNaxNjGQfuRn332Pdkg/6tb09fgY61+8b9H30ZMqoaJx44hW77MtRfpfqfgAuWZEJZy+cDPObjjEzDS1cVdcHrfbOinKlUmlRb2trQ0FBQb9PLxAR2cLISMDVHejlAAIAIDTRLuOQE2JA6ENAQAAAQKPRWNR37tyJmpoaJCU590IoInJMUhkwYQEgkeCeISFiKuA10q5jkRPhKYY+REREICEhAdu2bYOfnx9CQkKQnZ2No0ePAoDgCEJ2djYAoKioCACQk5ODkpISeHp6Yt68efYdnoicmv8YQL0UKP0GaLx2t+7uC4RPBkYNj2tF0SBhQOiDVCrFwYMHkZ6ejpUrV8Lf3x/Lli3Dc889h40bNyIhIcFi+8WLF1v8+eevQKpUKlRUVNhrbCIaJnxDgElLgdZG4OT/3qlNffqnIwtE94EBwQrR0dHIy8uzqKWlpSEuLg7u7pZfp7JmAVNnZyeMRiNMJhM6OzvR3t4OuVzuMF+fIsd0Q1+NA/m7MHNiCt75bC2kUhmiQydh5cK/iD1av1U3N+D104exNG4q1uXtgUwihVoZgYx+fvvAln1++4/3ob/dOuAeYvP4lwtv8q2EbIFrEAaoqKhowAsUn332Wbi7u+Prr79Gamoq3N3dUVlZaeMJydmc0eRAHT0HQb4q7Eo/gb+s+hqNhnqU1wy9GxMdr7yEWap4jPYJwJdLNiIvZRPqW/W4pKsSrU/KuGn31YPI2fAIwgAYDAZoNBqsWrVqQM/PzMxEZmambYcip3Lhaj42716EiOBE1DaUIzJkIjwVvli96C24y726t5NJXSB18O+5F1RdxuJPX0fCSBUq9DokBqngK/fEG7Oegpebons7F4kMMknP/2bpqcfHv3oBBdeu2LSPwsWtzx5EwwV/AwbAy8sLRqMRq1evFnsUclLx4TMQE/YwXluZj4SIZKxe9DY6OlstwkHZ9YvQt9yAKsix78YzPSwWDwWPRe4TL2N6WCzenPUbtHZ2WHyoX9RV4WZbM8YF9HwDqp56mM3mQenTVw+i4YIBgcgB1TSUIdgvAgCg019Do0GHiFF3v9De1NqAtz55Hi8s/j+xRrRaWWM9wh8IBABomxuga2tCQqCq+/GGNgN+dzwT785dYXWPEG8/XNBV2ryPNT2IhgsGBCIHVFn7PVTK8TCajJBIpDij+RLqqDkAAKOxCzv2PYlnf7kLfj7KPjqJ7/JNLcb5h8JoMkEqkSC34hJmqyYAALpMRjx19G3sSE6F0tO3u1bXou+1BwCb9+mpB9FwxoBA5IAq6r7HmKDx6OzqQKOhHudKjyM69M4tjAsuHoTm2mn875EX8cLfHsXlim9FnrZ3l29oMS4gFB3GTtS3NiGvqhhq5Z2jI9k/fIei2jJsKNiH2VlbcOq6BhV6HTZ9c6DXHjWGWzhbV27TPj31IBrOJGZHubA40TDS33sxFFw4gOTEJfe1z9nr7uvpVrHmXgwHS05hceyUez5+SFOIEQpPzBw93iH6/DvXXSlWb2tvP/9c2eP/NTk/fouBaAi433DgSPr6MP519MMO1YdouOIRBCIiJ8EjCGRLXINAREREAgwIREREJMCAQERERAIMCERERCTAgEBEREQCDAhEREQkwIBAREREAgwIREREJMCAQERERAIMCERERCTAgEBEREQCDAhEREQkwIBAREREAgwIREREJMCAQERERAISs9lsFnsIIiK6KzfD/vucvW7w99G5ft/g78QBuO5KEXsEm+ARBCIiIhJgQCAiIiIBBgQiIiISYEAgIiIiAQYEIqIh5pahHnPWS1Bafd6iXnerEnPWS1BZd0WkycRT36KHW0YqztdXWNQr9Tq4ZaTiys1qcQYbwhgQiIiGGM2103BzUSBcOcGiXlJVCA+5N8JGxog0mXiKasugcHHFhIAwi/rp2qvwdnNHjF+wSJMNXQwIRERDjOZaESJDJkImc7Go/3CtEFGhakilw++t/UxdGR4MHAMXqcyiXlRzFUlB4ZBKht9rcr/4itmQVqvFmjVrMHXqVHh4eEAikaC4uFjssYjIyWi0RYgOe0hQL7lWiJhQYX04OFNbhknKCEH9dO1VqHuoU98YEGyotLQU+/fvh6+vL5KTk8Ueh4iclKb6jCAImEwmlFaf7Q4Of3h3Nh7bFICPcl8VY0S7O1dXDrVyrEXNZDbhXF0FJikjcKO1CdM/egWzsv6ER/a8hNyKSyJNOnS49L0JWWvGjBmoq6sDAGRmZuLYsWMiT0REzuaGvhoNTTWIClVb1K9eP4+2DgPiVFMAAOuXZuLsj7m4odeKMaZdVTc3oKalEUlB4Rb1C/WVMHS2Y3JwFEYovJCfshkyqRSlt2qRdvgtzB4TL9LEQwOPIFjJZDIhIyMDUVFRUCgUSExMREFBAWJiYrBixQoAGJbn/YjIvnQ/feB7ufta1E+c24u40VMQ6Htnkd5I31C7zyaWakMDAMBX7mFRz7pyEpODIxHm4w+ZVArZT+/RzbfbkRA42u5zDjX8RLPS8uXLsWXLFqSnp+OLL77AkiVLkJKSgrKyMqjV6r4bEBHZwOjAOHgqHkBW3g4Y2hrR1NqAz799B59/+1c8PX+H2OOJItZvFB6Qe2BX4WdobG9BQ5sB753PxTvnc7B1xt37IpQ31uPRfZuxIHs7FkZOEnHioYGnGKywd+9e7N69G/n5+d1rC2bOnImzZ8/i0KFD/QoIHR0dSE9PR05ODlpaWpCUlIS3334bcXFxgzU+ETkRT4UPXl1+GO8feREpr4ZC4eqBiFEPYtszx5AQMUPs8UThI/fAJ4vWYcNXWQh/93l4uMiRGKjC4cdexPSwu++t4b6ByE/ZjAq9DnP2b8GCsUkiTu34GBCssH37dsydO1ew8DAyMhKurq6Ij7f+PFZXVxciIyOxdetWKJVK/PnPf8bSpUtx8eJFW49tQSKRDGp/IrKdnF2932R3Qvgv8Mbz/7TpPu3xHnF73d5B6z0tNBYFqZvv+XhHVyfkLq4AAB83d3i5uQ/aLI7+fmvtTZwZEPqg1WpRXFyMtWvXCh6rqqrC+PHjIZfLre7n6emJl156qfvPq1evxsaNG9He3g6FQmGTmYmIMg4sR0nVd+js6kBJ1XfYsvxzsUcS1bn6cmz4KgsyiRRdJiNem5km9kgOjwGhD1rtnQVBSqXSot7W1oaCggLMnz//vvqfPHkSY8aMGfRwYG1iJCLx5Wbcf491Sz7o1/b2eI/oXL9v0PdxL1NGRePEE6/YZV/O8n7LRYp9CAgIAABoNBqL+s6dO1FTU4OkpIGfw7p16xaee+45bN269b5mJCIisjUeQehDREQEEhISsG3bNvj5+SEkJATZ2dk4evQoAAgWKGZnZwMAioqKAAA5OTkoKSmBp6cn5s2b171dW1sbFi5ciKVLlyI1NdVOfxsiIiLrSMzOcixkEGk0GqSnp6OwsBD+/v5YtmwZvL29sXHjRjQ1NcHd/e5il3stTlGpVKioqABwZ6HiokWLMHLkSHzwQf8OAxKR87PFKYb+mr1u8Pch5ikGe3LdldL3RkMAjyBYITo6Gnl5eRa1tLQ0xMXFWYQDwLpzT8888wxMJhPee+89m85JROQMqpsb8Prpw1gaNxXr8vZAJpFCrYxARj8XFtqyz2//8T70t1sH3GMoYkAYoKKiIkyZMqXfz6usrMTu3buhUCjg63v3SmiXL1/G6NG8shcRDdwNfTUO5O/CzIkpeOeztZBKZYgOnYSVC/8i9mj9crzyEmap4jHaJwBfLtkIhYsb/uvIW7ikq0L8SOvfJ23ZJ2XcNPw6+uEB9xiKGBAGwGAwQKPRYNWqVf1+rkqlcpoVrkTkWM5ocqCOnoMgXxV2pZ+Am6sC2/f+J8prLiE82DHvO1BQdRmLP30dCSNVqNDrkBikgq/cE2/Megpebne/3eUikUF2j1s299Tj41+9gIJrV2zaR+Hi1mcPZ8KAMABeXl4wGo1ij0FEw9iFq/nYvHsRIoITUdtQjsiQifBU+GL1orfgLvfq3k4mdYFUKhNx0t5ND4vFQ8FjceTxP+I3R/+KV6c/gXV5eyw+1C/qqnCzrRnjAnq+v0RPPcxmM1o7O2zep68ezsT5IxARkROKD5+BmLCH8drKfCREJGP1orfR0dlqEQ7Krl+EvuUGVEHjRJy0d2WN9Qh/IBAAoG1ugK6tCQmBqu7HG9oM+N3xTLw7d4XVPUK8/XBBV2nzPtb0cCYMCEREQ1BNQxmC/SIAADr9NTQadIgYldj9eFNrA9765Hm8sPj/xBrRKpdvajHOPxRGkwlSiQS5FZcwWzUBANBlMuKpo29jR3IqlJ6+3bW6Fn2vPQDYvE9PPZwdAwIR0RBUWfs9VMrxMJqMkEikOKP5EuqoOQAAo7ELO/Y9iWd/uQt+Pso+Oonr8g0txgWEosPYifrWJuRVFUOtvBN8sn/4DkW1ZdhQsA+zs7bg1HUNKvQ6bPrmQK89agy3cLau3KZ9eurh7HgdBCIiB2PNdRD2Ht+GcapHEDt6Mlb/z2T4+4zCtqe/gFQqxYlz+/DXT9dAFTQeAPD0vO0YN+aRXvs5ynUQDpacwuLYe39D7JCmECMUnpg5erxD9OmJs1wHgQGBiMjBDORCSQUXDiA5ccmA9+koAcEZMCAQERGR0+IaBCIiIhJgQCAiIiIBBgQiIiISYEAgIiIiAQYEIiIiEmBAICIiIgEGBCIiIhJgQCAiIiIBBgQiIiISYEAgIiIiAQYEIiIiEmBAICIiIgEGBCIiIhJgQCAiIiIBBgQiIiISYEAgIiIiAQYEIiIiEmBAICIiIgEGBCIiIhL4f5MwgnwjHX9EAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 670.628x204.68 with 1 Axes>"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "q = QuantumRegister(3)\n",
    "qc = QuantumCircuit(q)\n",
    "\n",
    "qc.u3(pi/2,pi/2,pi/2,q[0])\n",
    "qc.cu1(pi/2,q[0], q[1])\n",
    "qc.u3(pi/2,pi/2,pi/2,q[1])\n",
    "qc.cu1(pi/2,q[1], q[2])\n",
    "qc.u3(pi/2,pi/2,pi/2,q[2])\n",
    "\n",
    "qc.draw(output='mpl')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[[ 0.354+0.j    -0.   -0.354j -0.   -0.354j -0.354+0.j    -0.   -0.354j\n",
      "  -0.354+0.j    -0.354+0.j     0.   +0.354j]\n",
      " [ 0.   +0.354j -0.354+0.j     0.   +0.354j -0.354+0.j     0.354-0.j\n",
      "   0.   +0.354j  0.354-0.j     0.   +0.354j]\n",
      " [ 0.   +0.354j  0.354-0.j    -0.354+0.j     0.   +0.354j  0.   +0.354j\n",
      "   0.354-0.j    -0.354+0.j     0.   +0.354j]\n",
      " [-0.354+0.j    -0.   -0.354j  0.354-0.j     0.   +0.354j -0.354+0.j\n",
      "  -0.   -0.354j  0.354-0.j     0.   +0.354j]\n",
      " [ 0.   +0.354j  0.354-0.j     0.354-0.j    -0.   -0.354j -0.354+0.j\n",
      "   0.   +0.354j  0.   +0.354j  0.354-0.j   ]\n",
      " [-0.354+0.j    -0.   -0.354j -0.354+0.j    -0.   -0.354j -0.   -0.354j\n",
      "   0.354-0.j    -0.   -0.354j  0.354-0.j   ]\n",
      " [-0.354+0.j     0.   +0.354j -0.   -0.354j -0.354+0.j     0.354-0.j\n",
      "  -0.   -0.354j  0.   +0.354j  0.354-0.j   ]\n",
      " [-0.   -0.354j  0.354-0.j     0.   +0.354j -0.354+0.j     0.   +0.354j\n",
      "  -0.354+0.j    -0.   -0.354j  0.354-0.j   ]]\n"
     ]
    }
   ],
   "source": [
    "backend = Aer.get_backend('unitary_simulator')\n",
    "job = execute(qc, backend)\n",
    "result = job.result()\n",
    "print(result.get_unitary(qc, decimals=3))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "------------"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### All Gates together in 4-qubit circut"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAApEAAADWCAYAAACJ6dk2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3XtclGX+//HXzHD0gAaomChiguIBUszcysxSM9ttt60kbc1yS7ay77e+Wrsbu+VjNc103e2obYe1X2qmbluuh93UlLbyfDY1SkBEAVFCRVFgZn5/TKAjxxuHuQd8Px8PHuJ133PNZwaYec9139d1W5xOpxMREREREQOsZhcgIiIiIo2PQqSIiIiIGKYQKSIiIiKGKUSKiIiIiGEKkSIiIiJimEKkiIiIiBimECkiIiIihilEioiIiIhhCpEiIiIiYphCpIiIiIgYphApIiIiIoYpRIqIiIiIYQqRIiIiImKYQqSIiIiIGKYQKSIiIiKGKUSKiIiIiGEKkSIiIiJimEKkiIiIiBimECkiIiIihilEioiIiIhhCpEiIiIiYphCpIiIiIgYphApIiIiIoYpRIqIiIiIYQqRIiIiImKYQqSIiIiIGKYQKSIiIiKG+ZldgDQepc98aHYJXuE/c5TZJVxxAlav8vp9lgy9w+v3WVdPLfD+ff71Ae/fp4g0bhqJFBERERHDFCJFRERExDCFSBERERExTCFSRERERAxTiBSfcOzMSQJmjWbnsUy39kMn8wmYNZr9J46YU5iIiIhUSSFSfMLW3HSC/PzpFd7RrX1L7kFaBgTTLbS9SZWJiIhIVRQixSdsy0vn2rad8bPa3Nq35hykb7torBb9qoqINGUOe8PuL56nd2YPWrt2LWPHjiUmJoZmzZrRuXNnHnzwQTIyMswuzedty02nX0SXSu1bcg+SWEW7iIg0Hafy4Ot3ofBo3fbP+xY2zoPikw1altRCIdKD5s6dy9GjR3n22WdZtWoVU6dOZdOmTSQmJpKZmWl2eT5tR14GiRHXuLU5nA525GXSL6ILx8+eYuCC57lt0Z/4yQd/YE3mHpMqNZfTaXYFIiKed2Q3nDsFO5bWHiTzvoW9y+HsD3AszTv1GXWlvFbrijUe9Oabb9KmTRu3toEDBxIdHc2cOXOYMWOGSZX5tiOnC8g5U0jfdtFu7buOHaKo9BzXt4/hqqAWrB81GZvVyvc/5DJm+esM6dzbpIq961QuZG2DY9+5Dt8Et4LIBOiQAH4BZlfnfc7iYsoefBjb/0zAOvAmV9v589h/+xyEhmL7w++xWK+8z8dH075i86cvkp+5g9JzpwkOaUu3G0Zzw31TzS5NpFbdboOy85B3wBUk+9wLra+uvF95gHQ6ofP10Kmf92utjsMOOfsgeweczgeLBUI7Q6e+ENbZ7OoaxpX3SltPDoeDWbNmERMTQ1BQEAkJCaSmptKtWzfGjx8PUClAAkRFRREeHk52dra3S240jhQVANA6sJlb+6L9X3N9+650DAnDZrVi+zEYnC45R3zbTl6v0ww538DmBZB7ABxlgBOKC+G7VNiyAErOml2h91mCg7Heew/2BR/idDpx2u3Yp04Hf39sv3v2igyQmbv+zacv30GPgWMZ98ohkt8qYMSTi2kT1cfs0kTqxGqFniOgXXewl1Q9InlpgLzmJldQ8wWOMtj5T9j/H1eAxAlOB5zIcD2Wg1+ZXWHDuPJebetp3LhxTJkyheTkZFatWsXIkSMZNWoU6enpJCYmVnu7vXv3kp+fT69evbxYbePSPfRqWgU2Y+bmZRSeO0NBcRF/27mGuTtX8+LNF65jnVF4jFs+nMydS6dzV1cf+vjZQIqOwzf/Bpw/fl3izAnY9x9vV+UbrHf9DE4U4PzyK+x/fQ3n8ePYJj+PJcDf7NJMsXvNG3SKv53YAUnY/AKw+QXQrks/YvrfY3ZpInVWU5D05QAJrpBYkPnjfy5+vf7x+4wNkH/Qy0V5gUJkHSxcuJD333+fZcuWMWnSJAYPHkxKSgo/+clPKCsrqzZElpaWkpycTHh4eMVoJUBeXh7Dhg2jWbNmJCQksGPHDm89FJ8UEtiMT+6exLbcDKLfmkDv9ybxyXdbWH7Pb7m5Y1zFftGt27J+1GS+/tVUnlo7z7yCvSR7J1WGx4sdPwhnC71Sjk+xBAdhve8e7DNn49y9G79pU7A0b1b7DZuo4JZtOLx3DTv+/QoFR/abXY5IvVUVJNM3+HaAtJf8+HpdEwtkbfVKOV6lcyLrYPr06QwfPpxBgwa5tXft2hV/f3969658bp7T6eSRRx5h8+bNrFixgrCwsIptjz32GN27d+fTTz/lgw8+4N577yUtLQ2bzVapH0+xeOAvrmTSQg9UUrUbI7uTOnpytdvPl5US6OcaZQoJCKZFQHCD1eKJ58oTPvh9BhGhnWvdL2nEYyzfMLfhC2pA/p+trN8Nz53DljQSy1VXGb6pr/ycq/K/842dlT9ozCuEtIlm3xd/54sFT9MytCPX/3IyPQc9XOc+fPn5kCuP1Wrjd/d/wOA+o0j/8VDwwrUv8vdn/mBuYVVIuOYWZv1mXc07OeGHw+DvF0CZvdQ7hdWT08CsII1E1iI7O5u9e/dy3333VdqWlZVFz549CQwMrLRtwoQJzJ8/n/nz5zNs2LCK9tOnT7NixQqef/55goODGT9+PHa7nY0bNzbo42jsdhzL4NZFf2LoR1O555M/8+fBY8wuqcH513HWjL/typtd41jzOY6PFmO5fRj2f35q6EWvKQoIbsn1d/+RB6bt5NHXc+l87Z2sfecRTp/QudjSODkcdr7c+8+K/5eWnWfj/uUmVlQ9I6/BNmvTGrtrWo+mAZRPiImIiHBrLy4uJjU1lREjRlS6zcSJE5kzZw7vvvsuSUlJbtu+++47wsLCCA8Pr2jr3bs3+/bt48Ybb2yAR+DiiTfZ0mc+9EAl9TPg6lg+v/95r9yXrwSS7UugIItaD2m/Pf8VPo56xSs1NZSA1avqvK9j8xbsr7+B7U+TscTEUPbgwzi/+C+WQTcbuk9f+TlX5akF9b9ts1Ztiel/L3vWzqHMwMwrX34+5Mpz8TmQAP5+gbzxfxuqnbVtpuKT8NXbte8X2BKKz5/1qUPxl0sjkbUoD3tpae6LUb388svk5OTQt29ft/aUlBRmz57Na6+9xsMPVz6UdObMGUJCQtzaQkJCKCoq8nDl0thFXkvNAdLiWu4n9MqYqA6A45t92Ke9hO2ZiVjje184N3L+hzgdDrPLM8WWZdM5tGc1JeeKcDocHMvYzpeLnqVT72Fc1T7W7PJEDLs0QELNs7bNFtwKwqKBWsJhx2t961xOT9BIZC26dOlCfHw806ZNIzQ0lA4dOrB06VJWrnSdw3XxpJqZM2cybdo0kpKSSExMdDtEHRISQo8ePWjevDmnT592u49Tp07RokUL7zwgaTTaXANtYiD/uyo2WlwvRnG3N70Xpeo4MzKxPz8ZW/KjWG+8oaLdetdPcSz92DUaecugGnpomsrOn+W/C57m9PEssFppGdaJuBvHkDDsSbNLEzHs0lnYmZtc7T1/POhX2zqSZom9BbYshLISqvzw37ItRDbBFbcsTh3DqFVaWhrJycls3ryZsLAwxo4dS8uWLUlJSeHUqVMEB7smedxyyy2kpqZW2cegQYNYv349p0+fJjw8nKNHj1ZMtomOjmb+/PkNejjbE8w8nO1N/jNH1b6TlzjscPC/kL0LLj4Xu0Ub6HYrXNXRvNo8ycjhbE8pGXqH1++zri7ncHZ9/fUB79+nyMWqWsZn7Z9d24ZMAocDvlnpCpK2AN8LkmdOwP41UHj4QpvFChFxrtdrv8rTJxo9jUTWQWxsLOvWuc+8GjNmDHFxcRUBEmD9+vW19tWyZUvuvPNOpkyZwksvvcT8+fOxWCwMGDDA02Wb6sjpAmZvWU5S3A1MWvcBNouVxIguzDI4IcaT/fzmP29zsuRsvfswg9UGMbdA9A2w/lVXW/9fQct2V84IpIg0fXVZB7J8+R/wzRHJ5mHQLwnOFMCG91xtA38DAU149TGdE1lPW7durXGR8ZrMmTOHb775hquuuopXX32Vf/zjHw26vI8Z1h7aw21RvekUEs5nI1NYN+oFjp09yZ78LNP6GdXjxsvqw0wXT9QOiVCAFJGm5ewPdVsH0m0dyVI4d9K7ddZF89AL3zflAAkaiayXoqIi0tLSePzxx+t1+3bt2rF69WoPV2We1Kx93PfpbOLbRJF5Mp+EdlG0DmzOK7c9RIuAoIr9/Cw2bJaqP7dU1cc/fjGR1MP7PdpP0I9prKY+RETEu6IHQKurXafo1PYhuTxIRl4LV0V6pz6pmt5F66FFixbY7XaefFInrgMM7Nid69pfw5r7/8jAjt159baHOVt63i347c7P4kTxaXqEV/0XX1UfTqezQfqprQ8REfG+0E51P8pitSpA+gKFSLls6YXHiG7VFoDs0wXkF58ivm1UxfaC4iKeWjuPt4aPr66LSn10aBnKrvxDHu+nLn2IiIhI7RQi5bLtO5FNj7BI7A4HVouFNZl7GBLVC4Ayh52HVr7BS4NGE9G8dUVb3pmTNfYBeLyfqvoQERGR+lGIlMu273g2PcIjOW8v5djZU6zL2ktiRBcAln67ia256TyX+iFDFk1h49E0Mk/m88KXi2vsI6foB7bnZXi0n6r6EBERkfrRxBq5bL8b8IuK73c8NIMlBzZi/XHSyv1xN3B/3A1u+3+ctpmkS9ou7QPg7pj+Hu2nqj5ERESkfhQixePu617zmpe/jO3vU/2IiIiIcTqcLSIiIiKGaSRS6syXLgcoTYsvX4LQDLoEoYg0BhqJFBERERHDFCJFRERExDCFSBERERExTCFSRERERAxTiBQRERERwxQiRURERMQwhUgRERERMUwhUkREREQMU4gUEREREcMUIkVERETEMIVIERERETFMIVJEREREDFOIFBERERHDFCJFRERExDCFSBERERExzM/sAkSuRGtmef+2QybV/z4bQukzH5pdQgX/maPMLqHRC1i9ypT7LRl6hyn3KyIaiRQRERGRelCIFBERERHDFCJFRERExDCFSBEf9EPRMYY+Y+H7Izvd2vN+OMTQZywcyttvUmXmOHbmJAGzRrPzWKZb+6GT+QTMGs3+E0fMKUxE5AqmECnig9IObyHAL4joiF5u7QeyNtMssCUd23QzqTJzbM1NJ8jPn17hHd3at+QepGVAMN1C25tUmYjIlUshUsQHpR3eStcOfbDZ3BdQ+PbwZmIiE7Far6w/3W156VzbtjN+Vptb+9acg/RtF43VcmU9HyIivkCvvCI+KC17K7Edr6vUfuDwZrpFVm5v6rblptMvokul9i25B0msol1ExBsKDkGhgbNpjn4DxScbrh5vU4j0oLVr1zJ27FhiYmJo1qwZnTt35sEHHyQjI8Ps0qSRSTuyrVJYdDgcfH9ke0W4fPatIdzzQjgL1kw1o0Sv2pGXQWLENW5tDqeDHXmZ9IvowvbcdAYtnMyti/7EsI+mkl6YZ1KlInKlOH0Mdv4TdvyjbkHy8A7Ytwq2LwZ7ScPX5w1abNyD5s6dS2FhIc8++yyxsbEcPnyYKVOmkJiYyPbt2+ncubPZJUojcPzkEQpO5RATmejWfvDoTorPFxEXNQCAZ5Lmsf27NRw/mW1GmV5z5HQBOWcK6dsu2q1917FDFJWe4/r2MfhZrSy/97e0DAhmVfoO/vT1P5g34nGTKva+krOQewDOnwa/QGgbC81Dza5KfIXDAcfT4dRRwAKtr4awaNBZIJeneTi06Qp5B1xBss890LpD1fse3gHfrnV93zERbAHeq7MhKUR60JtvvkmbNm3c2gYOHEh0dDRz5sxhxowZJlUmjUn+j6GwRXBrt/bPdywkrtMA2rZ2TS5p0zrS67WZ4UhRAQCtA5u5tS/a/zXXt+9Kx5Awt/YAmx9+V8i7o9MJ6V9D5iZwOi60H/zS9ebWcwT4+eCblfP8ecrG/hrryHux/fIXF9ozMiib9Dts/zsB680DTayw6SjMhj3L4XyRe3tQCPT+GbTSnLR6s1pdf2NQc5C8OEDG3gqd+nq3zoZ0ZbzSeoDD4WDWrFnExMQQFBREQkICqampdOvWjfHjxwNUCpAAUVFRhIeHk53dtEeLxHM6tY2jeVArFq17iaLiQk6dLeBfG+byrw1v8usRL5ldntd1D72aVoHNmLl5GYXnzlBQXMTfdq5h7s7VvHiz++UKz5ScI+WLRTx93U9Nqta7Mja4vi4OkOXyv3cdaqtqm9ksgYFYk+7DsfQfOEtLAXDm5lH23PNYH3pQAdJDTuXB9iVw/kzlbedOuw6rnjnh/bqakvIg2a676xD1pYe2m3KABIXIOhs3bhxTpkwhOTmZVatWMXLkSEaNGkV6ejqJiYnV3m7v3r3k5+fTq1evavcRuVjzoBCmjltO2uGtjJoaya9f7s6Xez5m2iP/JuGaQWaX53Uhgc345O5JbMvNIPqtCfR+bxKffLeF5ff8lps7xlXsV2IvY9S/XuH3A35Bz/CmP0pbchYyNta8T+FhOOGjp2Rb73QN4ThXr8VZeJKy5/6AdcRwbD+70+TKmo70r12HsnFWsdEJ9rLaf4ekdlUFyXJNOUCCDmfXycKFC3n//fdZv349gwa53sQHDx7M9u3b+fjjj6sNkaWlpSQnJxMeHl4xWgnwwgsvsGTJEg4cOMDixYu59957vfI4pPHoFX0Tr0z4yuwyfMaNkd1JHT252u12h4Mxy1/nFzH9+XnMlTF7PfdAHUYZLXBkL4RfU8t+JrAE+GO9Pwn74iVYVq7C2udabGMeMLusJqPkDBw/WMtOTsj7FroPcZ1LK/V36aHtizXVAAkaiayT6dOnM3z48IoAWa5r1674+/vTu3fvSrdxOp088sgjbN68mQULFhAWduG8rZiYGF555RX69+/f4LWXs1gs+vKhL0+YtXgcS1Jn8tnWefzxvZ81ut8BT1qatpHPMnexcN+XDFk0hafWzjN0e7Ofi/p8TX1+Jg5HLSnSCV+t2+azP0/rHbfD2WIIC8P6xGP16sPsn4OvfvWIvbZOz5/TAZ0jY02vt7bfK7PrqMuXzWbh2p/78e3hLRV1L/1iNlGJ5tdWn+e8LjQSWYvs7Gz27t3L008/XWlbVlYWPXv2JDCw8ke4CRMmMH/+fBYuXMiwYcPctv3qV78C4MUXX2yYouWKMGnke2aX4DOSut9AUvcbzC7Dq86cO1nrovMOh52i4h+8VFE92GxQch7r0NuwXGEL6De0MwYWIzxzrgktXGiyO68fT7eL1vi9o/8j/Hf3UvYd2mBiVQ1HIbIW5RNiIiIi3NqLi4tJTU1lxIgRlW4zceJE5syZw7vvvktSUpJX6qyN01nVSTFiljWzvH+fvvY7UPrMh2aXUMHXnpu6OFMAG2r5HGG12njgiSE8+1bDP76A1auM3+hwNpwtxhIbU+/7bYw/O2/Z9P/gdD5VnxMJruV+OsAPp313XdXy18rG8HO+eBJNueZBIbz+f1/XuPxPY6aPfrUIDw8HIC0tza395ZdfJicnh7593U90SElJYfbs2bz22ms8/PDDXqtTRK4szUNdy/hUywKBLVwn+/sqZ9p30Lo1lrZtzS6lSeo8gOoDJK5tnb13VlWTduks7HLVzdpuKjQSWYsuXboQHx/PtGnTCA0NpUOHDixdupSVK1cCuE2qmTlzJtOmTSMpKYnExEQ2brww7S0kJIQePXp4vX4Rabp6jnAt41N4GLDgFhgCm0Ofe31znchyzrTvsMTUlITlcrSLhXO3wHfrq9hoge63QbiuGnrZqlrGJ+1z1//rso5kY6YQWQur1cqSJUtITk7mscceIywsjLFjx/LEE0+QkpJCfHx8xb4rVqwA4KOPPuKjjz5y62fQoEGsX7/em6WLSBPnFwCJ97mW8TmyF/K/c7V3HwoRcb4dIAFsT/zG7BKavKh+EB4N2bvg8HZXW6d+EHktNGtd822ldrWtA1nXBckbK4XIOoiNjWXdunVubWPGjCEuLo7g4OCKtrqGxNLSUux2Ow6Hg9LSUs6dO0dgYKDhWVFyZTl+8giL189kcJ9RzF32NFarjdjIfjx211/MLq1BHDldwOwty0mKu4FJ6z7AZrGSGNGFWYPHmNKPr7JYXUv4hF9z4fyxyARzaxLf0jwMut16IUTG3mJqOU3Gqby6rQN5aZDc9Qnc9GjTuPShzomsp61bt9a4yHhNHn30UYKDg/nvf//L6NGjCQ4O5tChQx6uUJqabWmrSYwdSrvWUcxM/py/PP5fCouOkZGzx+zSGsTaQ3u4Lao3nULC+WxkCutGvcCxsyfZk59lSj8iIhcLaQfX3FS3dSDLg2T7XtDrzqYRIEEjkfVSVFREWloajz/+eL1uP2/ePObNm+fZoqRJ2XVwPZPfv5su7RPILciga4c+NA9qzZN3v05wYIuK/WxWP6xWm4mVekZq1j7u+3Q28W2iyDyZT0K7KFoHNueV2x6iRUBQxX5+Fhu2aq6LXVUf//jFRFIP7zfUj4hIXUUPqPu+Viv0HN5wtZhBr6L10KJFC+x2O08++aTZpUgT1Tv6Zrp17M+fH1tPfJdBPHn3G5wvPesWINOP7ubkmeNEtWv8E7YGduzOde2vYc39f2Rgx+68etvDnC097xb8dudncaL4ND2quaRhVX04nU7D/YiISN0oRIr4oJyCdNqHuqZN5p88TGFRPl2uvnCi26mzBbz+yQQm3veuWSV6VHrhMaJbuZZ5yT5dQH7xKeLbRlVsLygu4qm183hr+PjquqjUR4eWoezKP2S4HxERqRuFSBEfdCj3G6IiemJ32LFYrGxL+4zEmKEA2O1lvPThr3j0pzMJDYmopafGYd+JbHqERWJ3OLBaLKzJ3MOQqF4AlDnsPLTyDV4aNJqI5q0r2vLOnKyxD6DWfkREpP4UIkV8UGbeN3Ru15PSsvMUFh1jx/driY3sB0Dq7iWkHd7COyt+y8Q5t7Avs/FfTmvf8Wx6hEdy3l7KsbOnWJe1l8QI10js0m83sTU3nedSP2TIoilsPJpG5sl8XvhycY195BT9wPa8jBr7ERGR+tPEGhEfNPq25yq+f3viHlJ3La64TvKtfUZxa59RZpXWIH434BcV3+94aAZLDmzE+uPEl/vjbuD+OPfrYn+ctpmkS9ou7QPg7pj+NfYjIiL1pxAp0ggMShhpdgledV/3mqc8/jK2btdqq60fERGpP4VIERMMmWR2Bebzn9m0RlOvdCVD7zC7BBHxMp0TKSIiIiKGKUSKiIiIiGEKkSIiIiJimEKkiIiIiBimECkiIiIihilEioiIiIhhCpEiIiIiYphCpIiIiIgYphApIiIiIoYpRIqIiIiIYQqRIiIiImKYQqSIiIiIGKYQKSIiIiKGKUSKiIiIiGEKkSIiIiJimJ/ZBUjTV/rMh6bcr//MUabcr8iVKGD1KlPut2ToHabcr5nWzDLn9kMmXd79StOjkUgRERERMUwhUkREREQMU4gUEREREcMUIkVERETEMIVI8QnHzpwkYNZodh7LdGs/dDKfgFmj2X/iiDmFiYg0Qj8UHWPoMxa+P7LTrT3vh0MMfcbCobz9JlUmTYlCpPiErbnpBPn50yu8o1v7ltyDtAwIpltoe5MqExFpfNIObyHAL4joiF5u7QeyNtMssCUd23QzqTJpShQixSdsy0vn2rad8bPa3Nq35hykb7torBb9qoqI1FXa4a107dAHm819Jb9vD28mJjIRq9X3XlNLzzXs/uJ5vvdb1IitXbuWsWPHEhMTQ7NmzejcuTMPPvggGRkZZpfm87blptMvokul9i25B0msol1ERKqXlr2V2I7XVWo/cHgz3SIrt5vtxCH46m3IP1i3/Q/vgA3vQVF+w9YlNVOI9KC5c+dy9OhRnn32WVatWsXUqVPZtGkTiYmJZGZmml2eT9uRl0FixDVubQ6ngx15mRXhcvjiF2n/xnimbfinGSWK+DSH3ewKrhyldjhdDCVlZldSvbQj2yqFRYfDwfdHtleEy2ffGsI9L4SzYM1UM0p0U5AJZedh97Lag+ThHfDtWig5C4VHvVKeVENXrPGgN998kzZt2ri1DRw4kOjoaObMmcOMGTNMqsy3HTldQM6ZQvq2i3Zr33XsEEWl57i+fQwA79zxGz4/tJfs0wVmlCnik4oLIXML5HxzoW3nxxDVH66KNK+umjj+8xn2N+ZW3lBSAoDfP5dgCQ72clW1yymEtftgxyGwO8BqgYROcFsPiAw1u7oLjp88QsGpHGIiE93aDx7dSfH5IuKiBgDwTNI8tn+3huMns80o003Xm10fhA5vdwXJ+LugzTWV9ysPkACxt0JkgnfrFHcaiawjh8PBrFmziImJISgoiISEBFJTU+nWrRvjx48HqBQgAaKioggPDyc72/w/Ul91pMgVClsHNnNrX7T/a65v35WOIWEARLYM83ptIr6sKB82fQBHdoHjolGx4xmwbRHk7DOvtppYbx+G/7KP3b5sE58Cmw3b0//jkwHyYB7M/jdsy3AFSACHE3Yegr/8B/b70IhY/o+hsEVwa7f2z3csJK7TANq2dk1gbNPadz5lWCwQOxg69gWnveoRyUsDZKe+3q9T3ClE1tG4ceOYMmUKycnJrFq1ipEjRzJq1CjS09NJTEys9nZ79+4lPz+fXr16VbvPla576NW0CmzGzM3LKDx3hoLiIv62cw1zd67mxZt1/WuRqjidsOtTKCupaqPrn32roPikV8uqF8eatdhf/jO2/3sK6+3DzC6nklI7vPsFlNkrntoKTsDhgL//F4qr+lmYoFPbOJoHtWLRupcoKi7k1NkC/rVhLv/a8Ca/HvGS2eVVq6ogWU4B0jfpcHYdLFy4kPfff5/169czaNAgAAYPHsz27dv5+OOPqw2RpaWlJCcnEx4eXjFaef78eZKTk1m9ejVnzpyhb9++vPHGG8TFxXnt8fiakMBmfHL3JJ77YhHRb02gmV8gCW2jWH7PbxnY8cp9XkRqUnDIdSi7Jk4nHNkNXQd6p6b6cKz8N/Y35mD73TNYB95kdjlV2nEIztYQEJ24zo/cmgHOsOKiAAATGUlEQVQDfWDlnOZBIUwdt5y3V/yWUVMjCfJvRperr2XaI/8mvsvNZpdXo/IgCa5D2+UUIH2TQmQdTJ8+neHDh1cEyHJdu3bF39+f3r17V7qN0+nkkUceYfPmzaxYsYKwMNeh2LKyMrp27cqLL75IREQEM2bMICkpid27dzfoY7BYLA3af01KJi2sdZ8bI7uTOnqyR+/XzMcs0tDGDX+R+2/9fY2/506nkxUfbeCpm29s8Hr8P1tp+Db2T5bheOc9bH94DutPrq/X/Xrj73zo+PfofuMYrLbq3zId9jJeevMTbn71vgavZ/XMS8dDK+sVfROvTPjKo/fr7dfUx+76C78c+FTF/1//5Ek+feZ1r9ZQX+U/o8b4PuR01v77VU4hshbZ2dns3buXp59+utK2rKwsevbsSWBgYKVtEyZMYP78+SxcuJBhwy4cnmnevDl/+MMfKv7/5JNPkpKSwrlz5wgKCmqYB9FEPPrvt9ic8z0l9jI253zPJ798xuySRExjtdpw4sRC9W9SFosF2yVrr/oK++KlOD5YgG3y81j7+fbQksVSt+fQ4qPPdXVmLR7HgaxNlJad50DWJqaM+5fZJbnJzk+r+N5uLyW3QMvl+RqFyFqUT4iJiIhway8uLiY1NZURI0ZUus3EiROZM2cO7777LklJSTX2//XXX9O5c+cGD5BGPll4WukzH3qkn7eHJxva38zHLNLQ8r6FPXV4zx9y1/U4X234v4WA1avqvK99/kIci5dimzIZ67WXN73WG3/nqQfgn9tq3sdq8+N/xt3N8r80fD1rZnmmn0kj3zO0vzdfUy8+B7J1JBRm+/Pio8urnbXta8p/Rk39fUgTa2oRHh4OQFpamlv7yy+/TE5ODn37un+CTklJYfbs2bz22ms8/PDDNfb9ww8/8MQTT/Diiy96tmgRafLadAX/YKhhIBLwvSVQ7O/Nw7H0Y2zTp1x2gPSW67qAfy2DjFYLDOjqnXqauksn0SQm1TxrW8yjkchadOnShfj4eKZNm0ZoaCgdOnRg6dKlrFzpOv/n4kk1M2fOZNq0aSQlJZGYmMjGjRsrtoWEhNCjR4+K/xcXF3PXXXeRlJTE6NGjvfeARKRJsNqg152uNSGdUHnaMNDlBmhReeUx0zgPHsSxaDHYbNh//0cuXR/d+vOfYft1zR++zdAsAO6/Hj742pXZL36qy/9/z3UQ4nsrEzU61c3CvniyTU3rSIp3WZxNfazVA9LS0khOTmbz5s2EhYUxduxYWrZsSUpKCqdOnSL4xzXNbrnlFlJTU6vsY9CgQaxfvx5wTa65++67adOmDe+9Z+xwQmPkqcPZRvnP1PJA0vQVHoHvv4TCwxfagltD9PVwdeU5fw3GyOFsTyoZeofX7mvfEVi1Gw5fdL2Dq1vD8HiI7+i1Mjx2ONuoIZMatv/alvFxOiFtnStIWmy+HSTLf0YN/ZyZTSORdRAbG8u6devc2saMGUNcXFxFgAQqQmJtHnnkERwOB3/72988WaaIXIFad4B+SXC2EL5+x9V2w69dS6WIZ/Xo4PrKOwXTfzwf9ZkReq49oS7rQF66/M/uZZBwF4T7aJC8EihE1tPWrVsZMGCA4dsdOnSI999/n6CgIFq3vnA1gX379tGpUydPltjoHTldwOwty0mKu4FJ6z7AZrGSGNGFWYPHmF2aiM9pdtHFSRRqGla7kAvfN6bn+vjJIyxeP5PBfUYxd9nTWK02YiP78dhdfzG7NMontte2DuTFQfLILrBoZoepFCLroaioiLS0NB5//HHDt42Kimrys7U8Ze2hPdwW1ZtOIeF8NjKFIL8AHlzxOnvys+jdRoFbRMSIbWmrSYwdSrvWUcxM/pwA/yCmL3yAjJw9RLf34rkPVegQD62uhhbhte9bHiQjE6C5roZrKoXIemjRogV2+6WnhMvlSM3ax32fzia+TRSZJ/NJaBdF68DmvHLbQ7QIuLD8kZ/Fhk0fPUVEarTr4Homv383XdonkFuQQdcOfWge1Jon736d4MAWFfvZrH5YfWR9y7oEyHIWiwKkL9C7sfiEgR27c137a1hz/x8Z2LE7r972MGdLz7sFyN35WZwoPk2P8EgTKxUR8X29o2+mW8f+/Pmx9cR3GcSTd7/B+dKzbgEy/ehuTp45TlS7HjX0JFI9hUjxCemFx4hu1RaA7NMF5BefIr5tVMX2guIinlo7j7eGjzerRBGRRiOnIJ32oV0AyD95mMKifLpcfWFdzlNnC3j9kwlMvO9ds0qUJkAhUnzCvhPZ9AiLxO5wYLVYWJO5hyFRvQAoc9h5aOUbvDRoNBHNW9fSk4iIHMr9hqiIntgddiwWK9vSPiMxZigAdnsZL334Kx796UxCQyJq6UmkegqR4hP2Hc+mR3gk5+2lHDt7inVZe0mMcH2KXvrtJrbmpvNc6ocMWTSFjUfTaulNROTKlpn3DZ3b9aS07DyFRcfY8f1aYiP7AZC6ewlph7fwzorfMnHOLezL3GBytdJYaWKN+ITfDfhFxfc7HprBkgMbsf44geb+uBu4P+4Gs0oTEWl0Rt/2XMX3b0/cQ+quxVitrtfUW/uM4tY+uhiDXD6NRIpPuq+78TU4RUSkaoMSRppdgjRBGomUBqfLD4o0fd68/OCVrqlfSk8aD41EioiIiIhhCpEiIiIiYphCpIiIiIgYphApIiIiIoYpRIqIiIiIYQqRIiIiImKYQqSIiIiIGKYQKSIiIiKGKUSKiIiIiGEKkSIiIiJimEKkiIiIiBimECkiIiIihilEioiIiIhhCpEiIiIiYphCpIiIiIgY5md2ASIiIlJ3AatXmXK/JUPvuKzbP7XAQ4UY8NcHvH+fVxKNRIqIiIiIYQqRIiIiImKYDmeLiDRSZefh2PdwKgeKjl9o/2YltIyAtl0hKMS8+pqSonOwJxuyTkBu4YX2RRuhYyjEd4SWwebVJ2IGhUgRkUampBjSv4Kcb8BeWnl7zj7XV9o6CO8CXW+CFm28X2dTUHgWVuyCHZlQ5qi8feNB19fH21xB8qfXQlgLr5cpYgqFSBGRRiT/IOz/D5ScrcPOTjh+EE5kQJcboPP1YLE0eIlNxpZ0+MdWOFdFUL+U3QE7DsE3R+DnfeHGmIavT8RsCpEiIo3E0T2w7z/Gb+d0wMEvobgQ4m5XkKyLNd/A8p3Gb1dSBks2Q+EZGJGg51qaNk2s8aC1a9cyduxYYmJiaNasGZ07d+bBBx8kIyPD7NJEpJE7kVG/AHmxo3vh4Feeqacp25JevwB5sdXfwJdpnqlHfFfeAdfRgbpwOiFzMxTlN2xN3qQQ6UFz587l6NGjPPvss6xatYqpU6eyadMmEhMTyczMNLs8EWmkSs/VLUAOmeT6qknmJjiZ45m6mqLCs65D2LX56wO1r0G4bAfkn/ZMXeJ7TuXB3hWwe1ntQdLphPSv4fsvYPtSsJd4p8aGpsPZHvTmm2/Spo372esDBw4kOjqaOXPmMGPGDJMqE5HG7NBmOF/koc6ckPY5XKdFmKu0YlfdzoGsi1I7/GsHjLvZM/3Vh+M/n2F/Y27lDSWuFOP3zyVYgn1zWvk7E65mwL1T6HXLryvanE4ncx9txdDk9+l63d0mVgct20JkHzi83RUk4++CNtdU3q88QGZsACwQOxhsAV4vt0FoJLKOHA4Hs2bNIiYmhqCgIBISEkhNTaVbt26MHz8eoFKABIiKiiI8PJzs7GxvlywiTYC9DI7s8WyfJ3Ncoyjiruicaxa2J+057BrdNIv19mH4L/vY7cs28Smw2bA9/T8+GyCLCo5wpjCHNp0S3NpPHkun5Nxp2nXpZ1JlF1h+DIQd+4LTXvWI5KUBstedENHdlHIbhEJkHY0bN44pU6aQnJzMqlWrGDlyJKNGjSI9PZ3ExMRqb7d3717y8/Pp1auXF6sVkabih0NQWuz5fnP3e77Pxm5PdtXL+FwOJ65Z277CsWYt9pf/jO3/nsJ6+zCzy6lWXvoWLFYbYZE93dqPZ+2iWat2tAzraFJl7moLkk05QIJCZJ0sXLiQ999/n2XLljFp0iQGDx5MSkoKP/nJTygrK6s2RJaWlpKcnEx4eHjFaCXAAw88QLt27WjVqhX9+/dnw4YN3nooItLInMxtmH5PNVC/jVnWicbVr1GOlf/G/pdXsf3uGay3DTa7nBrlpW/hqohY/ALcR0rzs3bRNtr8UciLVRUkyzXlAAk6J7JOpk+fzvDhwxk0aJBbe9euXfH396d3796VbuN0OnnkkUfYvHkzK1asICwsrGJbSkoKf//73wkICGDFihXcc889HD16tMEfh4g0PmcaKIA0VL+N2cVXomkM/Rph/2QZjnfew/aH57D+5Hqzy6lVXvoWCvO+563fhLu1l54vot/Pfm9SVdUrD5LgOkfywoamGyBBIbJW2dnZ7N27l6effrrStqysLHr27ElgYGClbRMmTGD+/PksXLiQYcPcDxn06NEDcAVNf39/cnNzOXfuHEFBQQ3zIACLFisTaZSmjlvO9XF3urXVNgO7uu1rZl34/tTJM1gsurTKxUZN3Ubbzn3d2mqbgV3V9qcWuP9/f9r3WCyeW33c/7OVhva3L16K44MF2CY/j7Vf39pvUI3LfR/53/nOOu+bl7GV6385mbibHnRrX/D73rQzMBJpxnvfqxM2EBc1AIBFn8/g3Um/83oNl8PprPvPSYeza1E+ISYiIsKtvbi4mNTU1CoPZU+cOJE5c+bwzjvvkJSUVGW/DzzwAMHBwdx+++08/vjjDRogRaTxKik71yD9lpY2TL+NWVlJwzwnDdVvXdjnL8QxfyG2KZMvK0B6U2Hu95w/8wNR8bfTMiyy4steeo7zZwtp6wOTaqrz4LDJFQES4JcDn6r0IbAp0UhkLcLDXUPpaWlpjBgxoqL95ZdfJicnh7593f8oU1JSmD17Nq+//joPP/xwtf0uWLCAefPmsXLlSk6fbviFxIx8shAR33HwS8jY6N528YjixcpHIKvbfrFOsWF6XbjEhxth0yWzay8dVSxXPgJZ3faLDbmxF/M9+FwHrF5Vp/3s783DsWw5tulTsPbsWfsNanG5vy91ea7AdSjbL7BZpZnZOd99TYuwjjRv1a7O9+mt3/FLZ2H3HAGncuDw9kBefHR5tcv/NHYKkbXo0qUL8fHxTJs2jdDQUDp06MDSpUtZudJ1OOHikciZM2cybdo0kpKSSExMZOPGC6/8ISEhFYexy/n7+/Pzn/+chIQE+vfvT2xsrHcelIg0Gi0jat/Hl/ptzDqGVg6RnurX25wHD+JYtBhsNuy//yP2S7Zbf/4zbL+ufqDDTHnpW2gXfR1Wm3tEyfl+g6FD2d5S3TI+5edB1raOZGOmEFkLq9XKkiVLSE5O5rHHHiMsLIyxY8fyxBNPkJKSQnx8fMW+K1asAOCjjz7io48+cutn0KBBrF+/vsr7KCkpITMzUyFSRCoJiwK/QCg779l+I7p5tr+mIL4jfLwVHB4cvLIACZ0811+d7/eaawyfO+krbv7V7Crbb314jpcrqV1N60BeOtmmKQZJhcg6iI2NZd26dW5tY8aMIS4ujuCLFmqtLiRe7MSJE3z++efceeed+Pn58fbbb3P06NFKh8VFRABs/nB1L8ja5rk+W7aFkPae66+pCAl2BT5PrusYdzWEt/Rcf+I76rKQeFMPkppYU09bt26tcZHxmrz22mu0b9+edu3asWjRIlasWFFx7qWIyKU6Xw/+zTzXX+xg15ubVHZnAgR4aHjFZoW7ND7QZJ0+Vrd1IC9dR3L/f3Tt7CtaUVERaWlpPP7444ZvGxYWxhdffNEAVYlIUxXQDOKGwu5Pa96vLhNqOibCVb5xsQ+fFN4Sft4Hlmypeb+6TBIZEQ8RrTxTl/iekHYQd7vraEFt60CWB0mbP7Tp2nSuna0QWQ8tWrTAbr/0NGURkYbTNgZib4W0zy+jj1iIGVT7fle6G2Lgh7Ow5pv693FTLNzao/b9pHHrUPlaI9WyWKDrwIarxQwKkSIijUSnvhAQDAfWGJxoY4Go6+Cam8Cqk5hqZbG4Dmu3CoZlO6DUwJiBzeoagby1h04ZkKZPIVJEpBGJiHMdjv7uC8j71nWOVU1aR7pGH1tpIo0hFgsM7Abdr4Z/7YA9h6GmSdsWXJNoftYH2rf2VpUi5lKIFBFpZAJbQK8REDsIcr+FU7lQlO86Wd9qg2ah0LKd6xB4C83ZuyxtWsK4m+GHM7AzC7JOuK6FXWIHf5vrnMeOoa5Z3ZqFLVcahUgRkUYqoLnrELc0vKuaw+A4s6sQ8S06O0ZEREREDFOIFBERERHDLE5vXZ1cRERERJoMjUSKiIiIiGEKkSIiIiJimEKkiIiIiBimECkiIiIihilEioiIiIhhCpEiIiIiYphCpIiIiIgYphApIiIiIoYpRIqIiIiIYQqRIiIiImKYQqSIiIiIGKYQKSIiIiKGKUSKiIiIiGEKkSIiIiJimEKkiIiIiBimECkiIiIihilEioiIiIhhCpEiIiIiYphCpIiIiIgY9v8BIKbXRWrQ+eYAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 851.228x264.88 with 1 Axes>"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "q = QuantumRegister(4)\n",
    "qc = QuantumCircuit(q)\n",
    "\n",
    "qc.u3(pi/2,pi/2,pi/2,q[0])\n",
    "qc.cu1(pi/2,q[0], q[1])\n",
    "qc.u2(pi/2,pi/2,q[1])\n",
    "qc.u1(pi/2,q[2])\n",
    "qc.x(q[0])\n",
    "qc.ccx(q[1],q[2],q[3])\n",
    "qc.y(q[1])\n",
    "qc.z(q[2])\n",
    "qc.cx(q[2],q[3])\n",
    "qc.z(q[3])\n",
    "qc.h(q[3])\n",
    "qc.s(q[0])\n",
    "qc.cu1(pi/2,q[1], q[2])\n",
    "qc.swap(q[0],q[2])\n",
    "qc.cswap(q[0],q[1],q[3])\n",
    "\n",
    "qc.draw(output='mpl')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[[ 0.   +0.354j -0.354+0.j    -0.   -0.354j  0.354-0.j     0.   +0.j\n",
      "   0.   +0.j     0.   +0.j     0.   +0.j    -0.   -0.354j  0.354-0.j\n",
      "   0.   +0.354j -0.354+0.j     0.   +0.j     0.   +0.j     0.   +0.j\n",
      "   0.   +0.j   ]\n",
      " [ 0.   +0.j     0.   +0.j     0.   +0.j     0.   +0.j     0.354-0.j\n",
      "   0.   +0.354j -0.354+0.j    -0.   -0.354j  0.   +0.j     0.   +0.j\n",
      "   0.   +0.j     0.   +0.j    -0.354+0.j    -0.   -0.354j  0.354-0.j\n",
      "   0.   +0.354j]\n",
      " [-0.354+0.j    -0.   -0.354j -0.354+0.j    -0.   -0.354j  0.   +0.j\n",
      "   0.   +0.j     0.   +0.j     0.   +0.j     0.354-0.j     0.   +0.354j\n",
      "   0.354-0.j     0.   +0.354j  0.   +0.j     0.   +0.j     0.   +0.j\n",
      "   0.   +0.j   ]\n",
      " [ 0.   +0.j     0.   +0.j     0.   +0.j     0.   +0.j     0.354-0.j\n",
      "   0.   +0.354j -0.354+0.j    -0.   -0.354j  0.   +0.j     0.   +0.j\n",
      "   0.   +0.j     0.   +0.j     0.354-0.j     0.   +0.354j -0.354+0.j\n",
      "  -0.   -0.354j]\n",
      " [ 0.   +0.354j  0.354-0.j    -0.354+0.j     0.   +0.354j  0.   +0.j\n",
      "   0.   +0.j     0.   +0.j     0.   +0.j    -0.   -0.354j -0.354+0.j\n",
      "   0.354-0.j    -0.   -0.354j  0.   +0.j     0.   +0.j     0.   +0.j\n",
      "   0.   +0.j   ]\n",
      " [ 0.   +0.j     0.   +0.j     0.   +0.j     0.   +0.j     0.354-0.j\n",
      "  -0.   -0.354j  0.   +0.354j  0.354-0.j     0.   +0.j     0.   +0.j\n",
      "   0.   +0.j     0.   +0.j    -0.354+0.j     0.   +0.354j -0.   -0.354j\n",
      "  -0.354+0.j   ]\n",
      " [-0.354+0.j     0.   +0.354j  0.   +0.354j  0.354-0.j     0.   +0.j\n",
      "   0.   +0.j     0.   +0.j     0.   +0.j     0.354-0.j    -0.   -0.354j\n",
      "  -0.   -0.354j -0.354+0.j     0.   +0.j     0.   +0.j     0.   +0.j\n",
      "   0.   +0.j   ]\n",
      " [ 0.   +0.j     0.   +0.j     0.   +0.j     0.   +0.j     0.354-0.j\n",
      "  -0.   -0.354j  0.   +0.354j  0.354-0.j     0.   +0.j     0.   +0.j\n",
      "   0.   +0.j     0.   +0.j     0.354-0.j    -0.   -0.354j  0.   +0.354j\n",
      "   0.354-0.j   ]\n",
      " [ 0.   +0.354j -0.354+0.j    -0.   -0.354j  0.354-0.j     0.   +0.j\n",
      "   0.   +0.j     0.   +0.j     0.   +0.j     0.   +0.354j -0.354+0.j\n",
      "  -0.   -0.354j  0.354-0.j     0.   +0.j     0.   +0.j     0.   +0.j\n",
      "   0.   +0.j   ]\n",
      " [ 0.   +0.j     0.   +0.j     0.   +0.j     0.   +0.j     0.354-0.j\n",
      "   0.   +0.354j  0.354-0.j     0.   +0.354j  0.   +0.j     0.   +0.j\n",
      "   0.   +0.j     0.   +0.j    -0.354+0.j    -0.   -0.354j -0.354+0.j\n",
      "  -0.   -0.354j]\n",
      " [-0.354+0.j    -0.   -0.354j -0.354+0.j    -0.   -0.354j  0.   +0.j\n",
      "   0.   +0.j     0.   +0.j     0.   +0.j    -0.354+0.j    -0.   -0.354j\n",
      "  -0.354+0.j    -0.   -0.354j  0.   +0.j     0.   +0.j     0.   +0.j\n",
      "   0.   +0.j   ]\n",
      " [ 0.   +0.j     0.   +0.j     0.   +0.j     0.   +0.j    -0.354+0.j\n",
      "  -0.   -0.354j -0.354+0.j    -0.   -0.354j  0.   +0.j     0.   +0.j\n",
      "   0.   +0.j     0.   +0.j    -0.354+0.j    -0.   -0.354j -0.354+0.j\n",
      "  -0.   -0.354j]\n",
      " [ 0.   +0.354j  0.354-0.j    -0.354+0.j     0.   +0.354j  0.   +0.j\n",
      "   0.   +0.j     0.   +0.j     0.   +0.j     0.   +0.354j  0.354-0.j\n",
      "  -0.354+0.j     0.   +0.354j  0.   +0.j     0.   +0.j     0.   +0.j\n",
      "   0.   +0.j   ]\n",
      " [ 0.   +0.j     0.   +0.j     0.   +0.j     0.   +0.j     0.354-0.j\n",
      "  -0.   -0.354j -0.   -0.354j -0.354+0.j     0.   +0.j     0.   +0.j\n",
      "   0.   +0.j     0.   +0.j    -0.354+0.j     0.   +0.354j  0.   +0.354j\n",
      "   0.354-0.j   ]\n",
      " [-0.354+0.j     0.   +0.354j  0.   +0.354j  0.354-0.j     0.   +0.j\n",
      "   0.   +0.j     0.   +0.j     0.   +0.j    -0.354+0.j     0.   +0.354j\n",
      "   0.   +0.354j  0.354-0.j     0.   +0.j     0.   +0.j     0.   +0.j\n",
      "   0.   +0.j   ]\n",
      " [ 0.   +0.j     0.   +0.j     0.   +0.j     0.   +0.j    -0.354+0.j\n",
      "   0.   +0.354j  0.   +0.354j  0.354-0.j     0.   +0.j     0.   +0.j\n",
      "   0.   +0.j     0.   +0.j    -0.354+0.j     0.   +0.354j  0.   +0.354j\n",
      "   0.354-0.j   ]]\n"
     ]
    }
   ],
   "source": [
    "backend = Aer.get_backend('unitary_simulator')\n",
    "job = execute(qc, backend)\n",
    "result = job.result()\n",
    "print(result.get_unitary(qc, decimals=3))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}