4gt13_v3.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "##4grt11 Circuit\n",
    "##All Imports Here\n",
    "\n",
    "#Non-qiskit imports\n",
    "import numpy as pi\n",
    "import matplotlib.pyplot as plt\n",
    "import time\n",
    "\n",
    "#Qiskit imports\n",
    "from qiskit import QuantumRegister, ClassicalRegister, QuantumCircuit\n",
    "from qiskit.circuit.library import XGate, YGate, ZGate, HGate, CXGate, SwapGate, CCXGate, RZGate\n",
    "from qiskit.compiler import transpile\n",
    "from qiskit.providers.fake_provider import FakeValencia\n",
    "#from qiskit.quantum_info import Operator\n",
    "from qiskit.circuit.random import random_circuit\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "from qiskit import QuantumRegister, ClassicalRegister, QuantumCircuit\n",
    "from numpy import pi\n",
    "from qiskit.circuit.library import C3XGate\n",
    "\n",
    "qreg_q = QuantumRegister(5, 'q')\n",
    "creg_c = ClassicalRegister(4, 'c')\n",
    "cal_qc = QuantumCircuit(qreg_q, creg_c)\n",
    "\n",
    "cal_qc.cx(qreg_q[2], qreg_q[1])\n",
    "cal_qc.cx(qreg_q[1], qreg_q[2])\n",
    "cal_qc.cx(qreg_q[3], qreg_q[2])\n",
    "cal_qc.cx(qreg_q[2], qreg_q[3])\n",
    "cal_qc.cx(qreg_q[4], qreg_q[3])\n",
    "cal_qc.cx(qreg_q[3], qreg_q[4])\n",
    "cal_qc.append(C3XGate(), [qreg_q[1], qreg_q[2], qreg_q[4], qreg_q[0]])\n",
    "cal_qc.ccx(qreg_q[2], qreg_q[4], qreg_q[0])\n",
    "cal_qc.ccx(qreg_q[1], qreg_q[4], qreg_q[0])\n",
    "cal_qc.cx(qreg_q[0], qreg_q[4])\n",
    "# @columns [0,1,2,3,4,5,6,7,8,9]\n",
    "backend = FakeValencia()\n",
    "\n",
    "layers = cal_qc.depth()\n",
    "compiled_cal_qc = transpile(cal_qc, backend)\n",
    "\n",
    "# Save the circuit to a QASM file\n",
    "file_path = '4gt13_v3.qasm'\n",
    "with open(file_path, 'w') as file:\n",
    "    file.write(compiled_cal_qc.qasm())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Circuit 1:\n",
      "      ░ \n",
      "q_0: ─░─\n",
      "      ░ \n",
      "q_1: ─░─\n",
      "      ░ \n",
      "q_2: ─░─\n",
      "      ░ \n",
      "q_3: ─░─\n",
      "      ░ \n",
      "q_4: ─░─\n",
      "      ░ \n",
      "c: 1/═══\n",
      "        \n",
      "Circuit 2:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "           ░ \n",
      "q_1: ──────░─\n",
      "           ░ \n",
      "q_2: ──────░─\n",
      "           ░ \n",
      "q_3: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_4: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 3:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "           ░ \n",
      "q_1: ──────░─\n",
      "           ░ \n",
      "q_2: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_3: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_4: ──────░─\n",
      "           ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 4:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "           ░ \n",
      "q_1: ──────░─\n",
      "           ░ \n",
      "q_2: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_3: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_4: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 5:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "           ░ \n",
      "q_1: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_2: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_3: ──────░─\n",
      "           ░ \n",
      "q_4: ──────░─\n",
      "           ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 6:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "           ░ \n",
      "q_1: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_2: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_3: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_4: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 7:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "           ░ \n",
      "q_1: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_2: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_3: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_4: ──────░─\n",
      "           ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 8:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "           ░ \n",
      "q_1: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_2: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_3: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_4: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 9:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_1: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_2: ──────░─\n",
      "           ░ \n",
      "q_3: ──────░─\n",
      "           ░ \n",
      "q_4: ──────░─\n",
      "           ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 10:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_1: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_2: ──────░─\n",
      "           ░ \n",
      "q_3: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_4: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 11:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_1: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_2: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_3: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_4: ──────░─\n",
      "           ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 12:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_1: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_2: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_3: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_4: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 13:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_1: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_2: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_3: ──────░─\n",
      "           ░ \n",
      "q_4: ──────░─\n",
      "           ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 14:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_1: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_2: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_3: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_4: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 15:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_1: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_2: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_3: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_4: ──────░─\n",
      "           ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 16:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_1: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_2: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_3: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_4: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "c: 1/════════\n",
      "             \n"
     ]
    }
   ],
   "source": [
    "from qiskit import QuantumRegister, ClassicalRegister, QuantumCircuit\n",
    "\n",
    "n_qbits = 5\n",
    "data = 1\n",
    "\n",
    "qreg_q = QuantumRegister(n_qbits, 'q')\n",
    "creg_c = ClassicalRegister(data, 'c')\n",
    "\n",
    "in_circuits = []\n",
    "\n",
    "for i in range(16):\n",
    "    in_qc = QuantumCircuit(qreg_q, creg_c)\n",
    "\n",
    "    # Apply X gates based on the binary representation of 'i'\n",
    "    binary_rep = format(i, '04b')\n",
    "    for j in range(len(binary_rep)):\n",
    "        if binary_rep[j] == '1':\n",
    "            in_qc.x(qreg_q[j+1])\n",
    "\n",
    "    in_qc.barrier(qreg_q[0], qreg_q[1], qreg_q[2], qreg_q[3], qreg_q[4])\n",
    "\n",
    "    in_circuits.append(in_qc)\n",
    "\n",
    "# Printing the circuit representations\n",
    "for i, in_qc in enumerate(in_circuits):\n",
    "    print(f\"Circuit {i+1}:\")\n",
    "    print(in_qc)\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPEAAAGwCAYAAAB4uIN0AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy88F64QAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAnhElEQVR4nO3dfVRU56Eu8GeI6MjwoUjCUEBGBSqgQCOYom0MXm1Bg7Y32tCjJq7l0bTRq6VWzIdps9bhVKdltetY0x4tod60iRerOakFNE0v5pZwTw2IcCxgFBcfDs5EJjDI8CXMzP0jcW4IoMy4YfPu/fzWyqLM3nt8rDy8e+937z0al8vlAhEJy0fuAET0YFhiIsGxxESCY4mJBMcSEwmOJSYSHEtMJDiWmEhwLDGR4FhiIsGxxESCY4mJBMcSEwmOJSYSHEtMJDiWmEhwLDGR4FhiIsGxxESCY4mJBMcSEwmOJSYSHEtMJDiWmEhwLDGR4FhiIsGxxESCY4mJBDdN7gD0YK5cueLR+rdu3cLJkyfxne98B4888si4tlm4cKE30WiScCRWmfb2drz22mtob2+XOwpJhCUmEhxLTCQ4lphIcCyxygQFBSErKwtBQUFyRyGJaFwul0vuEOQ9T89Oe4Nnp6c2jsQqMzAwgJaWFgwMDMgdhSTCEqtMY2MjMjIy0NjYKHcUkghLTCQ4lphIcCwxkeBYYiLBcYpJcJxiIo7ERIJjiVWmqakJ2dnZaGpqkjsKSYQlVpne3l7U1tait7dX7igkEZaYSHAsMZHgWGIiwbHEKhMeHg6j0Yjw8HC5o5BEOE8sOM4TE0dileno6MCbb76Jjo4OuaOQRFhilTGbzcjLy4PZbJY7CklEFSW2Wq3Izc1FdHQ0tFotIiMjsWfPHvT09GDbtm3QaDQ4cuSI3DGJvKL4h8fX1NQgMzMTFosFOp0O8fHxuHnzJg4fPozr16+7dyuTk5PlDUrkJUWPxFarFVlZWbBYLNi7dy/MZjOqq6thsVhgNBpRUlKCyspKaDQaJCYmyh2XyCuKLvHu3bthMpmwa9cu5OfnIyAgwL0sNzcXSUlJGBoagsFgQGBgoIxJJ49Op8Py5cuh0+nkjkISUWyJGxoaUFRUhJCQEBw8eHDUdZYsWQIASEpKGvZ6U1MT1q1bh4CAAMyePRvPPPMMPvnkkwnPPBkMBgMKCgpgMBjkjkISUWyJT5w4AafTiU2bNsHf33/UdWbOnAlgeIm7u7uRnp4Ok8mEEydO4NixYygvL8eTTz4Jp9M5KdknksPhgN1uh8PhkDsKSUSxJ7bKysoAAOnp6WOuYzKZAAwv8bFjx9DW1oa//e1vmDt3LgAgIiICy5Ytw5kzZ/Ctb31r4kJPgitXrmDDhg04deoUEhIS5I5DElBsiVtaWgAAUVFRoy4fGhpCRUUFgOElLi4uxte+9jV3gQEgLS0N8+fPx5///GevSpySkgKLxeLxduOxdetWj9a/Oz9cWlqKS5cujWub48ePe5iKPKXX61FVVeXVtootcU9PDwCgr69v1OVFRUWwWq0ICAjAvHnz3K/X19dj48aNI9ZPSEhAfX29V1ksFgva2tq82vZ+PL0vuL+/3/11vNtOVHaShmJLrNfr0dnZierqaqSlpQ1bZjabsW/fPgBAYmIiNBqNe1lnZydmzZo14v2Cg4Px0UcfeZ1lovj5+Xm0vlardX8d77a8WWLiPcjPiGJLvGrVKjQ0NMBoNGL16tWIjY0FAFRWVmLLli2wWq0AJuciD293k8bD0xsg6urqUFhYiDVr1oz7mDgvL8+baDRJFHt2Ojc3F3PmzMGNGzeQkJCAxYsXIyYmBkuXLsX8+fOxcuVKACOnl2bPng2bzTbi/To6OhAcHDwZ0SdUbGwsKioq3L/USHyKLXFERATKy8uxdu1aaLVaNDc3Izg4GEePHkVJSQmuXr0KYGSJ4+LiRj32ra+vR1xc3KRkn0i+vr4IDg6Gr6+v3FFIIootMfBpIYuLi9Hd3Y3u7m5cuHABO3bsQE9PD5qbm+Hj44NFixYN2+bJJ5/EBx984J5+AoALFy7g+vXryMrKmuy/guRaW1vx/PPPo7W1Ve4oJBFFl3gsdXV1cLlciImJGXFyZ8eOHQgLC8P69etRXFyMU6dO4bvf/S6WLl2K9evXy5RYOt3d3Th//jy6u7vljkISUWWJL1++DGDkrjQABAYGoqysDGFhYcjOzsY///M/Y9myZSguLoaPjyr/76IpTrFnp+/lXiUGgAULFqC4uHgyIxF5TZVDy/1KTCQSVY7Ed6+rVqPQ0FDs378foaGhckchifBpl4Lj0y5JlbvTatbV1YVz586hq6tL7igkEZZYZUwmE3JycobNg5PYWGIiwbHERIJjiYkExxKrjFarRVxcnPu+YhIfp5gExykm4khMJDiWWGXq6+uRmJjo9fPCaOphiVXG5XJhcHAQPIpSDpaYSHAsMZHgWGIiwanyVkQ1W7BgAc6cOYPIyEi5o5BEWGKV0Wq1iImJkTsGSYi70yrT1taGAwcO8KNZFIQlVhmbzYbTp0+P+oB8EhNLTCQ4lphIcCwxkeBYYpUJCQnB9u3bERISIncUkghLrDIajQbTp08f9pnMJDaWWGXa29vx2muvob29Xe4oJBGWmEhwLDGR4FhiIsGxxCoTFBSErKwsBAUFyR2FJMIH5QmOD8ojjsQqMzAwgJaWFgwMDMgdhSTCEqtMY2MjMjIy0NjYKHcUkghLTCQ4lphIcCwxkeBYYiLBcYpJcJxiIo7ERIJjiVWmqakJ2dnZaGpqkjsKSYQlVpne3l7U1tait7dX7igkEZaYSHAsMZHgWGIiwbHEKhMeHg6j0Yjw8HC5o5BEOE8sOM4TE0dileno6MCbb76Jjo4OuaOQRFRRYqvVitzcXERHR0Or1SIyMhJ79uxBT08Ptm3bBo1GgyNHjsgdc1KYzWbk5eXBbDbLHYUkoviPNq2pqUFmZiYsFgt0Oh3i4+Nx8+ZNHD58GNevX3ePSMnJyfIGJfKSokdiq9WKrKwsWCwW7N27F2azGdXV1bBYLDAajSgpKUFlZSU0Gg0SExPljkvkFUWXePfu3TCZTNi1axfy8/MREBDgXpabm4ukpCQMDQ3BYDAgMDBQxqRE3lNsiRsaGlBUVISQkBAcPHhw1HWWLFkCAEhKSnK/drf0S5cuxYwZMxT3cSc6nQ7Lly+HTqeTOwpJRLElPnHiBJxOJzZt2gR/f/9R15k5cyaA4SVubGzE6dOnodfrkZqaOilZJ5PBYEBBQQEMBoPcUUgiii1xWVkZACA9PX3MdUwmE4DhJX788cdhNptx5swZrFq1amJDysDhcMBut8PhcMgdhSSi2BK3tLQAAKKiokZdPjQ0hIqKCgDDS+zjo9j/SwB8enFIamrqpFwkQpNDsVNMPT09AIC+vr5RlxcVFcFqtSIgIADz5s2b0CwpKSmwWCwT8t5bt271aP2788OlpaW4dOnSuLY5fvy4h6nIU3q9HlVVVV5tq9gS6/V6dHZ2orq6GmlpacOWmc1m7Nu3DwCQmJg44SevLBYL2traJuS9Pb0vuL+/3/11vNtOVHaShmJLvGrVKjQ0NMBoNGL16tWIjY0FAFRWVmLLli2wWq0AJuciD71eP2Hv7efn59H6Wq3W/XW82/JmiYn3ID8jii1xbm4u3nrrLdy4cQMJCQlYuHAh+vv70djYiMzMTBgMBrz77rvDjocnire7SePh6bFtXV0dCgsLsWbNGiQkJIxrm7y8PG+i0SRR7FmciIgIlJeXY+3atdBqtWhubkZwcDCOHj2KkpISXL16FQAmpcRTSWxsLCoqKtx7JiQ+xY7EABAXF4fi4uIRr9vtdjQ3N8PHxweLFi2SIZl8fH19ERwcLHcMkpCiSzyWuro6uFwuxMbGjnpceOrUKQBAfX39sO8NBgNSUlImL+gEaG1txaFDh/DCCy9g7ty5cschCaiyxJcvXwYw9q70xo0bR/3+2WefFX66pbu7G+fPn8fOnTvljkISYYlHwYedkEgUe2LrXu5XYiKRqHIkvntdNZESqHIkVrPQ0FDs378foaGhckchifBpl4Lj0y6JI7HKdHV14dy5c+jq6pI7CkmEJVYZk8mEnJwc973UJD6WmEhwLDGR4FhiIsGxxCqj1WoRFxfnvq+YxMcpJsFxiok4EhMJjiVWmfr6eiQmJrpvsyTxscQq43K5MDg4yDu1FIQlJhIcS0wkOJaYSHCqvJ9YzRYsWIAzZ84gMjJS7igkEZZYZbRaLWJiYuSOQRLi7rTKtLW14cCBA/xoFgVhiVXGZrPh9OnTsNlsckchibDERIJjiYkExxITCY4lVhkfHx+kpqbCx4f/9ErBf0mVcTqdqKyshNPplDsKSYQlJhIcS0wkOJaYSHAsscoEBQUhKysLQUFBckchifAZW4LjM7aII7HKDAwMoKWlBQMDA3JHIYmwxCrT2NiIjIwMNDY2yh2FJMISEwmOJSYSHEtMJDiWmEhwnGISHKeYiCMxkeBYYpVpampCdnY2mpqa5I5CEmGJVaa3txe1tbXo7e2VOwpJhCUmEhxLTCQ4lphIcCyxyoSHh8NoNCI8PFzuKCQRzhMLjvPExJFYZTo6OvDmm2+io6ND7igkEZZYZcxmM/Ly8mA2m+WOQhJRRYmtVityc3MRHR0NrVaLyMhI7NmzBz09Pdi2bRs0Gg2OHDkid0wiryj+o01ramqQmZkJi8UCnU6H+Ph43Lx5E4cPH8b169fdu5XJycnyBiXykqJHYqvViqysLFgsFuzduxdmsxnV1dWwWCwwGo0oKSlBZWUlNBoNEhMT5Y5L5BVFl3j37t0wmUzYtWsX8vPzERAQ4F6Wm5uLpKQkDA0NwWAwIDAwUMakk0en02H58uXQ6XRyRyGJKLbEDQ0NKCoqQkhICA4ePDjqOkuWLAEAJCUluV87deoUnnrqKURFRcHPzw8LFy7Eyy+/DLvdPim5J5rBYEBBQQEMBoPcUUgiii3xiRMn4HQ6sWnTJvj7+4+6zsyZMwEML3F+fj4eeugh/PSnP8XZs2fx/e9/H7/5zW+QkZGhiM8vcjgcsNvtcDgcckchiSj2xFZZWRkAID09fcx1TCYTgOEl/vOf/4yHH37Y/f2KFSvw8MMPY9OmTfjggw/w+OOPT1DiyXHlyhVs2LABp06dQkJCgtxxSAKKLXFLSwsAICoqatTlQ0NDqKioADC8xJ8v8F0pKSkAgLa2Nq+ypKSkwGKxeLXt/WzdutWj9e/OD5eWluLSpUvj2ub48eMepiJP6fV6VFVVebWtYkvc09MDAOjr6xt1eVFREaxWKwICAjBv3rx7vtf58+cBAHFxcV5lsVgsXv8CuB9P7wvu7+93fx3vthOVnaSh2BLr9Xp0dnaiuroaaWlpw5aZzWbs27cPAJCYmAiNRjPm+7S1teGVV15BRkaG13PJer3eq+3Gw8/Pz6P1tVqt++t4t+XNEhPvQX5GFFviVatWoaGhAUajEatXr0ZsbCwAoLKyElu2bIHVagVw74s87HY71q9fj+nTp6OwsNDrLN7uJo2HpzdA1NXVobCwEGvWrBn3MXFeXp430WiSKPbsdG5uLubMmYMbN24gISEBixcvRkxMDJYuXYr58+dj5cqVAIYfD39eX18fsrKy0NTUhL/85S8ICwubzPgTJjY2FhUVFe5faiQ+xZY4IiIC5eXlWLt2LbRaLZqbmxEcHIyjR4+ipKQEV69eBTB6iQcHB7FhwwZUVVXh7NmziI+Pn+z4E8bX1xfBwcHw9fWVOwpJRJX3E9vtdgQGBkKj0aC7u3vYsaHT6UR2djbOnDmD0tJS94g9VXm6O93a2opDhw7hhRdewNy5c8e1De8nntoUOxLfS11dHVwuF2JiYkac3Nm5cyf++Mc/IicnB35+fvj73//u/q+9vV2mxNLp7u7G+fPn0d3dLXcUkogqS3z58mUAo+9Knz17FgBw6NAhpKWlDfuvpKRkUnMSjYdiz07fy71K3NzcPMlpiB4MR2IiwalyJL57XbUahYaGYv/+/QgNDZU7CklElWenlYRPuyRV7k6rWVdXF86dO4euri65o5BEWGKVMZlMyMnJcd+GSeJjiYkExxITCY4lJhIcS6wyWq0WcXFx7vuKSXycYhIcp5iIIzGR4Fhilamvr0diYiLq6+vljkISYYlVxuVyYXBwEDyKUg6WmEhwLDGR4FhiIsGp8lZENVuwYAHOnDmDyMhIuaOQRFhildFqtYiJiZE7BkmIu9Mq09bWhgMHDvCjWRSEJVYZm82G06dPw2azyR2FJMISEwmOJSYSHEtMJDiWWGV8fHyQmpoKHx/+0ysF/yVVxul0orKyEk6nU+4oJBGWmEhwLDGR4FhiIsGxxCoTFBSErKwsBAUFyR2FJMJnbAmOz9gijsQqMzAwgJaWFgwMDMgdhSTCEqtMY2MjMjIy0NjYKHcUkghLTCQ4lphIcCwxkeBYYiLBcYpJcJxiIo7ERIJjiVWmqakJ2dnZaGpqkjsKSYQlVpne3l7U1tait7dX7igkEZaYSHAsMZHgWGIiwbHEKhMeHg6j0Yjw8HC5o5BEOE8sOM4TE0dileno6MCbb76Jjo4OuaOQRFhilTGbzcjLy4PZbJY7CklEFSW2Wq3Izc1FdHQ0tFotIiMjsWfPHvT09GDbtm3QaDQ4cuSI3DGJvKL4jzatqalBZmYmLBYLdDod4uPjcfPmTRw+fBjXr19371YmJyfLG5TIS4oeia1WK7KysmCxWLB3716YzWZUV1fDYrHAaDSipKQElZWV0Gg0SExMlDsukVcUXeLdu3fDZDJh165dyM/PR0BAgHtZbm4ukpKSMDQ0BIPBgMDAQBmTTh6dTofly5dDp9PJHYUkotgSNzQ0oKioCCEhITh48OCo6yxZsgQAkJSU5H6tvLwcq1atQlhYGGbMmIGIiAg8/fTTaGhomJTcE81gMKCgoAAGg0HuKCQRxR4TnzhxAk6nE5s2bYK/v/+o68ycORPA8BJ3dnZi8eLFeO655/DII4/AZDLh4MGDSEtLwz/+8Q9ERERMSv6J4nA40NfXh5kzZ+Khhx6SOw5JQLElLisrAwCkp6ePuY7JZAIwvMTr1q3DunXrhq2XmpqKL3/5yzh9+jT27NkzAWknz5UrV7BhwwacOnUKCQkJcschCSh2d7qlpQUAEBUVNeryoaEhVFRUABhe4tHMmTMHADBtmmJ/55HAFPtT2dPTAwDo6+sbdXlRURGsVisCAgIwb968EcsdDgecTidaWlrw4osvQq/X4zvf+Y5XWVJSUmCxWLza9n62bt3q0fp3L/IoLS3FpUuXxrXN8ePHPUxFntLr9aiqqvJqW8WWWK/Xo7OzE9XV1UhLSxu2zGw2Y9++fQCAxMREaDSaEduvWLHCPVJHR0ejrKwMDz/8sFdZLBYL2travNr2fjy9ub+/v9/9dbzbTlR2koZiS7xq1So0NDTAaDRi9erViI2NBQBUVlZiy5YtsFqtAMa+yOP111+HzWZDU1MTfv7zn+Mb3/gGKioqMHfuXI+z6PV6r/8e9+Pn5+fR+lqt1v11vNvyjqeJ9yA/I4q9i8lkMiE5ORmffPIJpk2bhoULF6K/vx+NjY3IzMyE0+nEu+++i2PHjmH79u33fC+bzQaDwYDNmzdPucszPb2LaXBwEN3d3QgICICvr++4tuFdTFObYk9sRUREoLy8HGvXroVWq0VzczOCg4Nx9OhRlJSU4OrVqwDuf1ILAGbNmoXo6GhFfH6Rr68vgoODx11gmvoUuzsNAHFxcSguLh7xut1uR3NzM3x8fLBo0aL7vs+tW7fw0Ucf4bHHHpuImJOqtbUVhw4dwgsvvODVoQFNPYou8Vjq6urgcrkQGxs74rhw8+bNiI6ORnJyMmbNmoVr167hl7/8JaZNm4acnByZEkunu7sb58+fx86dO+WOQhJRZYkvX74MYPRd6a9+9at444038G//9m/o7+9HZGQk0tPT8dJLL40550wkJ5b4C3bt2oVdu3ZNdiQiryn2xNa93KvERKJR5Uh897pqNQoNDcX+/fsRGhoqdxSSiGLnidWCT7skVe5Oq1lXVxfOnTuHrq4uuaOQRFhilTGZTMjJyXHfhkniY4mJBMcSEwmOJSYSHEusMlqtFnFxce5bEkl8nGISHKeYiCMxkeBYYpWpr69HYmIi6uvr5Y5CEmGJVcblcmFwcBA8ilIOlphIcCwxkeBYYiLBqfJWRDVbsGABzpw5g8jISLmjkERYYpXRarWIiYmROwZJiLvTKtPW1oYDBw7wUx0UhCVWGZvNhtOnT8Nms8kdhSTCEhMJjiUmEhxLTCQ4llhlQkJCsH37doSEhMgdhSTCEquMRqPB9OnTR/1MZhITS6wy7e3teO2119De3i53FJIIS0wkOJaYSHAsMZHgWGKVCQoKQlZWFoKCguSOQhLhg/IExwflEUdilRkYGEBLSwsGBgbkjkISYYlVprGxERkZGWhsbJQ7CkmE9xPTlORyudDb2yt3DI/4+fnJchENS0xTUm9vL/z9/eWO4RG73Q6dTjfpfy53p4kExxITCY670yqTkJCAhoYGuWOQhDgSEwmOJVaZpqYmZGdno6mpSe4oJBGWWGV6e3tRW1sr3PQNjY0lJhIcS0wkOJaYSHAsscqEh4fDaDQiPDxc7igkEZZYZWbNmoV169Zh1qxZckcRzty5c+WOMCpe7KEyHR0dOHv2LDIzMxEcHCx3nAm3cOFCpKSkYMmSJVi4cCH8/PzgdDpx+/Zt/Nd//RcuXryICxcuwGw23/N9Nm7ciD/84Q/43ve+h9/97neTlH58WGKVMZvNyMvLQ3JysmJLrNPp8E//9E94/vnnkZycPOZ669atAwA4nU6cO3cOv/71r3H27Fk4nc5h623cuBFvvfUWpk2bhoKCAly7dg0ffPDBRP4VPKL43Wmr1Yrc3FxER0dDq9UiMjISe/bsQU9PD7Zt2waNRoMjR47IHZMksnnzZrS2tuLYsWP3LPDn+fj4YM2aNSguLsbly5fx2GOPuZd9vsAAUFhYiIqKiomI7jVFj8Q1NTXIzMyExWKBTqdDfHw8bt68icOHD+P69evo6OgAgHH/Y9PUFRoait/+9rfIysoa9vrf//53/OlPf8LFixdRW1sLm80GjUYDvV6PRx99FKmpqfjud78Lg8EAAIiPj0dFRQXy8/NRW1uLN954w13ggoIC7NixA1PtiVaKfcaW1WrFV77yFZhMJuzduxc/+clPEBAQAAD42c9+hv3792PatGlwOByw2WwIDAyUObF3PH3GVl1dHTZs2IBTp04hISFhXNvI8Yytnp6ecd9PPG/ePPz1r3/F/Pnz3a+dPHkSRqMR1dXV993ex8cHGRkZePXVV5Gamup+3el0wsfn053V8RSY9xNLbPfu3TCZTNi1axfy8/PdBQaA3NxcJCUlYWhoCAaDQdgCe0On02H58uWy/LBNhMjISLz//vvuAlssFnzrW9/C008/Pa4CA5+WtbS0FGlpaXjxxRcxODgIAO4CFxYWTskR+C5FlrihoQFFRUUICQnBwYMHR11nyZIlAICkpKQx3yczMxMajQavvvrqRMSUhcFgQEFBgXv3UWTTpk3DO++84576qaurw6OPPoo//elPXr2fw+HA9evXRzxip6+vb8oWGFBoiU+cOAGn04lNmzaNuUs2c+ZMAGOX+OTJk6ipqZmoiLJxOByw2+1wOBxyR3lgL774Ih599FEAwLVr17By5cr7ThXdyxdPYg0NDQEAdu7cifT09AcPPEEUWeKysjIAuOf/8SaTCcDoJb59+zZ+8IMfID8/f2ICyujKlStITU2dlOdVT6SEhAS88sorAD4tW3Z2Nm7duuX1+32xwAUFBfjhD3/oXv7666+7f/FPNYo8O93S0gIAiIqKGnX50NCQe5pgtBK//PLLiI2NxaZNm7B58+YHzpOSkgKLxfLA7zOarVu3erT+3ZGqtLQUly5dGtc2x48f9zDVg/viXO0X/ehHP4Kvry8A4NChQ+M+/h3NaAXesWMHAOCpp57CihUrMG/ePGRnZ9/zQo+YmBj3cbSn9Ho9qqqqvNpWkSXu6ekB8OmxzGiKiopgtVoREBCAefPmDVtWVVWF3/72t7h48aJkeSwWC9ra2iR7v8/z9L7g/v5+99fxbjtR2b0VHByM7OxsAEBnZyd++tOfev1eYxX47jFwbm4uLly4AAB4/vnn71niB9mVfxCKLLFer0dnZyeqq6uRlpY2bJnZbMa+ffsAAImJicNOYjgcDjz33HPYtWvXuKdfxptnovj5+Xm0vlardX8d77Zy3CzhdDrHLMXmzZvdf4/f/e53Y/6yvp/7FRgAPvzwQ1RVVSElJQUpKSlITk4e81xJWFjYA43E3lJkiVetWoWGhgYYjUasXr0asbGxAIDKykps2bIFVqsVwMiLPI4cOYKPP/5Y8rPR3u4mjYc388SFhYVYs2bNuH9R5eXleRPtgdxrnvjrX/+6+397u6s/ngJ//s9ISUlx/9ljlfjatWucJ5ZKbm4u5syZgxs3biAhIQGLFy9GTEwMli5divnz52PlypUAhh8PW61WvPLKK/jxj3+MoaEh2Gw22Gw2AJ/uetpstvsep4kgNjYWFRUV7l9sIro7PWi321FXV+fx9p4UGPj0qq8v/tlTiSJLHBERgfLycqxduxZarRbNzc0IDg7G0aNHUVJSgqtXrwIYXmKTyYTu7m4899xzmD17tvs/ADAajZg9ezZaW1tl+ftIydfXF8HBwe6TQqKZNWuW+zzGpUuXPP7F6mmBAeDy5cu4c+cOALintKYSRe5OA0BcXByKi4tHvG6329Hc3AwfHx8sWrTI/Xp0dDTOnz8/Yv309HQ8++yz2Lp164Qe206W1tZWHDp0CC+88MKUvT/2Xu7+YgWAGzdueLStNwUGgDt37uDWrVuIiIgY9udPFYot8Vjq6urgcrkQGxs77MSOv78/nnjiiVG3MRgMYy4TTXd3N86fP4+dO3fKHcUrZrMZjz/+OLRaLdrb28e9nUajwQ9+8AOvb2Z45pln4OPjA7vd7lXuiaS6El++fBnAvS+3pKmrv78f5eXlHm/ncrmwZs0avPfee6itrfX4WujR9tKmCpb4PqbyNbPkma6uLqxcuRI9PT2K+ndliUlVpuLu8INSXYnvXletVqGhodi/fz9CQ0PljkISUV2J1S4kJMTj661palPkPDGNraurC+fOnUNXV5fcUUgiLLHKmEwm5OTkuG/FJPGxxESCY4mJBMcSEwmOJVYZrVaLuLg49/24JD5OManMggUL8Pbbb8sdgyTEkZhIcCyxytTX1yMxMRH19fVyRyGJsMQq43K5MDg4qKgbANSOx8Q0Jfn5+Ul6s8LPj/4v3O7pRaDOD/ueyx7xvRQ8fWihVFhimpI0Go2kD52bPkOL6YMOTJ+hhU6nG/G9yLg7TSQ4jsQqs2DBApw5cwaRkZFyRyGJsMQqo9VqERMTI3cMkhB3p1Wmra0NBw4cmHIfzULeY4lVxmaz4fTp0+4H45P4WGIiwbHERIJjiYkEx7PTglu4cKFH6wcFBeEnP/kJli1bhrCwsAlKRZOJJVaZsLAwyT+6leTF3WkiwbHERIJjiYkExxITCY4lJhIcS0wkOJaYSHAsMZHgWGIiwbHERIJjiYkExxJPEX/729+wfv16REVFQaPRIC8vT+5IdB+lpaVITk7GjBkzYDAY8Itf/EKWHCzxFGG32xEfH4+f/exn0Ov1cseh+6iqqsL69euRmZmJmpoavPrqq3jppZfw7//+75OehXcxTRFr1qzBmjVrAAD79++XOQ3dzy9+8Qukpqbi4MGDAIC4uDjU1dXh0KFD+N73vjepWTgSE3mhoqICGRkZw17LyMhAS0sLTCbTpGbhSEyK0nXbjo8/sY14fcjhcH+92mQa8f3nGcJDMX267z3/HLPZPOKw5+73ZrMZERER3v4VPMYSk6L4TvfFqdL/g9v2nlGX9/b1o/Bk6Zjff3l+JGIM4ROeU0rcnSZF8dPOwMa1K7ze9qnMFdBoNPddNywsDBaLZdhrH3/8sXvZZGKJSXFiDBFYtmSRx9t9+5tfR6D/+D7ZcPny5Xj33XeHvXbu3DlERUVN6q40wN3pKcNut6OxsREAcOfOHVgsFtTU1MDf3x/R0dEypxNP5oqluNZkQnuHbVzrfyUhBosXzh/3++fk5GDZsmV4+eWXsWXLFly4cAG/+tWv8Mtf/tLLxN7TuPhp01PC+++/j/T09BGvr1ixAu+///7kB1IAk7kdv/7DO3A67/0jHhSgQ862jdDOmO7R+5eUlOCll17ClStXoNfrsWfPHvzwhz98kMheYYlJ0f53RTXe+6Dqnutsz34SC6K+NEmJpMdjYsHcMN9Cb/+A3DGE8URaMiLDHhlz+ddSFgtdYIAlFsqQw4E//Md7MP7mLbS2fSx3HCE85OODp59Mh6/vyNM/j8yZjW+uSJUhlbRY4i9wOBz4/e9/j2984xt4+OGHMWPGDMydOxcZGRkoKCiA47OLBORw8fJH6OruwfTpvgh7ZI5sOUQTEhyEtelfHfbaQz4+eDorHb7TxD+3yxJ/zu3bt7F69Wo888wzeO+99zB9+nQkJSXB6XTiL3/5C7Zv347u7m5Zsg05HCj7v5cAAE88ljzqyEJjeyw5Dl+eH+n+ftXXliA8NETGRNJhiT9n27ZtOH/+PCIiIlBWVoa2tjZ8+OGHMJlMMJvN+Nd//Vf4+t77cryJcncUDvD3w9Ikzz5/iQCNRoOnMlfATzsDUeGhWPFYktyRJMOz05+5ePEiUlJSMG3aNFy6dAmLFnl+scBYfvU/30a3vc/r7V0uF+y9fXC5XJgxfTpm3Oe6Xhrb0NAQfB7ygY9mao1fAf4z8T+e/e9ebct9ss+88847AIC1a9dKWmAA6Lb3jXktr6cG7tzBwJ07krwXKQNL/Jn6+noAQFpamuTvHeA/0+ttOQqrw4P8jLDEn7l9+zaATz+/V2re7iYBwIWaBvzHu+UI8PdD7o5sntCiEfgT8ZnAwEAAQFdXl+Tv7e0x8d1RGADu3BnCz48VSR2NpggeE0sgISEBb7/9Nv7zP/9T8veW4piYx8I0Fpb4M9/+9rfxL//yLygtLUV9fT3i4+Mle29vjnd4LKwuD3JMzCmmz3n66adx8uRJzJ07F2+88QZWrPj/N5d//PHHKCwsxO7du6HT6SY8C4+FabxY4s+5ffs21q9f7771Lzw8HF/60pdgNpvR1tYGl8uFzs5OzJo1a0JzDDkcyD9WBNttO7L+2zIsT5F2youUZWrNeMssMDAQf/3rX/H666/jiSeeQG9vL2pra+Hj44NvfvObeP311xEQEDDhOS5evgrbbTuvzqJx4Ug8BV2qu4az73+IFY8lcRSm+2KJp6jBoSFooMG0aQ/JHYWmOJaYSHA8JiYSHEtMJDiWmEhwLDGR4FhiIsGxxESCY4mJBMcSEwmOJSYSHEtMJDiWmEhwLDGR4FhiIsGxxESCY4mJBMcSEwmOJSYSHEtMJDiWmEhwLDGR4FhiIsGxxESCY4mJBMcSEwmOJSYSHEtMJDiWmEhwLDGR4FhiIsGxxESCY4mJBMcSEwmOJSYSHEtMJDiWmEhwLDGR4FhiIsGxxESCY4mJBPf/ACcaHPzdUsn0AAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 287.496x535.111 with 1 Axes>"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "##Measurement/Output circuit = out_qc\n",
    "\n",
    "out_qc = QuantumCircuit(qreg_q, creg_c)\n",
    "\n",
    "out_qc.barrier(qreg_q[0], qreg_q[1], qreg_q[2], qreg_q[3], qreg_q[4])\n",
    "out_qc.measure(qreg_q[4], creg_c[0])\n",
    "\n",
    "out_qc.draw(output='mpl')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "                                   ┌───┐┌───┐┌───┐     \n",
      "q_0: ──────────────────────────────┤ X ├┤ X ├┤ X ├──■──\n",
      "     ┌───┐                         └─┬─┘└─┬─┘└─┬─┘  │  \n",
      "q_1: ┤ X ├──■────────────────────────■────┼────■────┼──\n",
      "     └─┬─┘┌─┴─┐┌───┐                 │    │    │    │  \n",
      "q_2: ──■──┤ X ├┤ X ├──■──────────────■────■────┼────┼──\n",
      "          └───┘└─┬─┘┌─┴─┐┌───┐       │    │    │    │  \n",
      "q_3: ────────────■──┤ X ├┤ X ├──■────┼────┼────┼────┼──\n",
      "                    └───┘└─┬─┘┌─┴─┐  │    │    │  ┌─┴─┐\n",
      "q_4: ──────────────────────■──┤ X ├──■────■────■──┤ X ├\n",
      "                              └───┘               └───┘\n"
     ]
    }
   ],
   "source": [
    "from qiskit import QuantumCircuit, transpile\n",
    "from qiskit.transpiler.passes import RemoveResetInZeroState, Unroller\n",
    "from qiskit.transpiler import PassManager\n",
    "from qiskit.circuit.library.standard_gates import IGate\n",
    "\n",
    "# Create a new circuit to store the modified version\n",
    "mod_qc = QuantumCircuit(cal_qc.num_qubits)\n",
    "\n",
    "# Iterate through the original circuit and replace NOT gates with Identity gates\n",
    "for instr, qargs, cargs in cal_qc:\n",
    "    if instr.name == 'h':\n",
    "        # If it's an X gate (NOT gate), replace with Identity gate\n",
    "        mod_qc.append(IGate(), qargs)\n",
    "    else:\n",
    "        # If it's not an X gate, keep the original instruction\n",
    "        mod_qc.append(instr, qargs, cargs)\n",
    "\n",
    "print(mod_qc)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Modified Circuit (before) for qubit 0:\n",
      "\n",
      "Modified Circuit (before) for qubit 1:\n",
      "\n",
      "Modified Circuit (before) for qubit 2:\n",
      "\n",
      "Modified Circuit (before) for qubit 3:\n",
      "\n",
      "Modified Circuit (before) for qubit 4:\n",
      "\n",
      "Modified Circuit (before) for qubit 5:\n",
      "\n",
      "Modified Circuit (before) for qubit 6:\n",
      "\n",
      "Modified Circuit (before) for qubit 7:\n",
      "\n",
      "Modified Circuit (before) for qubit 8:\n",
      "\n"
     ]
    }
   ],
   "source": [
    "from qiskit import QuantumCircuit, QuantumRegister\n",
    "from qiskit.circuit.library import XGate\n",
    "\n",
    "# Initialize an empty list to store the modified circuits\n",
    "modified_circuits_before = []\n",
    "\n",
    "\n",
    "# Iterate over each qubit\n",
    "for qubit in qreg_q:\n",
    "    # Create a new circuit for the modified version\n",
    "    modified_qc_before = QuantumCircuit(qreg_q, creg_c)\n",
    "    \n",
    "   \n",
    "    # Keep track of whether a NOT gate has been inserted before and after the current qubit\n",
    "    not_before_inserted = False\n",
    "\n",
    "    # Iterate over each gate in the original circuit\n",
    "    for instruction, qargs, cargs in cal_qc:\n",
    "        \n",
    "        # Check if the gate involves the current qubit\n",
    "        if qubit not in qargs:\n",
    "            # Insert a NOT gate before the gate if not already inserted\n",
    "            if not_before_inserted:\n",
    "                modified_qc_before.append(instruction, qargs, cargs)  \n",
    "            else:\n",
    "                modified_qc_before.h(qubit)\n",
    "                modified_qc_before.append(instruction, qargs, cargs)\n",
    "                not_before_inserted = True                       \n",
    "        else:\n",
    "            not_before_inserted=False\n",
    "            modified_qc_before.append(instruction, qargs, cargs)        \n",
    "            \n",
    "    #new_depth_b =  modified_qc_before.depth()\n",
    "    # Add the modified circuits to the list\n",
    "    #if new_depth_b == layers:\n",
    "        modified_circuits_before.append(modified_qc_before)\n",
    "            \n",
    "# Initialize a list to store the modified circuits\n",
    "new_circuits1 = []\n",
    "\n",
    "for modified_qc_before in modified_circuits_before:\n",
    "    # Get the indices of the NOT gates\n",
    "    not_gate_indices = [i for i, (instruction, _, _) in enumerate(modified_qc_before) if instruction.name == 'h']\n",
    "    # Iterate over the instructions in the circuit\n",
    "    for i, (instruction, qargs, cargs) in enumerate(modified_qc_before):\n",
    "        if instruction.name == 'h' and i in not_gate_indices:\n",
    "            # Create a new circuit for each NOT gate\n",
    "            new_qc = modified_qc_before.copy()\n",
    "            # Add the corresponding NOT gate to the new circuit\n",
    "            new_qc.h(qargs[0])\n",
    "            \n",
    "            # Add the remaining gates from the original circuit to the new circuit\n",
    "            for j in reversed(range(len(new_qc))):\n",
    "                if new_qc[j][0].name == 'h' and j != i:\n",
    "                    new_qc.data.pop(j)\n",
    "            \n",
    "            new_depth_b =  new_qc.depth()\n",
    "            #print(new_depth_b)\n",
    "            # Add the modified circuits to the list\n",
    "            if new_depth_b == layers:\n",
    "                new_circuits1.append(new_qc)\n",
    "            # Add the modified circuit to the list\n",
    "           \n",
    "\n",
    "def remove_duplicates(circuit_list):\n",
    "    # Initialize an empty set to store the unique circuits\n",
    "    unique_circuits = set()\n",
    "    \n",
    "    # Initialize a new list to store the circuits without duplicates\n",
    "    new_circuit_list = []\n",
    "    \n",
    "    # Iterate over each circuit in the original list\n",
    "    for circuit in circuit_list:\n",
    "        # Convert the circuit to a custom string representation\n",
    "        circuit_str = str(circuit)\n",
    "        \n",
    "        # Check if the circuit string is already in the set\n",
    "        if circuit_str not in unique_circuits:\n",
    "            # If not, add it to the set and the new list\n",
    "            unique_circuits.add(circuit_str)\n",
    "            new_circuit_list.append(circuit)\n",
    "    \n",
    "    return new_circuit_list\n",
    "\n",
    "# Your list of circuits (assuming it is already defined)\n",
    "# circuit_list = [circuit1, circuit2, circuit3, ...]\n",
    "\n",
    "# Call the remove_duplicates function to get the new list without duplicates\n",
    "new_circuit_list = remove_duplicates(new_circuits1)        \n",
    "\n",
    "\n",
    "# Display and save the modified circuits\n",
    "#for i, modified_qc_before in enumerate(modified_circuits_before):\n",
    "    #print(f\"Modified Circuit (before) for qubit {i}:\")\n",
    "    #modified_qc_before.draw(output='mpl', filename=f\"not_before_qubit_{i}.png\")\n",
    "    #print()\n",
    "\n",
    "for i, new_qc in enumerate(new_circuit_list):\n",
    "    print(f\"Modified Circuit (before) for qubit {i}:\")\n",
    "    new_qc.draw(output='mpl', filename=f\"h_qubit_{i}.png\")\n",
    "    compiled_new_qc = transpile(new_qc, backend)\n",
    "\n",
    "    file_path = f\"Hadamard_4gt13_v3_{i}.qasm\"\n",
    "    with open(file_path, 'w') as file:\n",
    "        file.write(compiled_new_qc.qasm())\n",
    "    print()\n",
    "\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Modified Circuit (before) for qubit 0:\n",
      "\n",
      "Modified Circuit (before) for qubit 1:\n",
      "\n",
      "Modified Circuit (before) for qubit 2:\n",
      "\n",
      "Modified Circuit (before) for qubit 3:\n",
      "\n"
     ]
    }
   ],
   "source": [
    "from qiskit import QuantumCircuit, QuantumRegister\n",
    "from qiskit.circuit.library import XGate\n",
    "\n",
    "# Initialize an empty list to store the modified circuits\n",
    "modified_circuits_before = []\n",
    "\n",
    "\n",
    "# Iterate over each qubit\n",
    "for qubit in qreg_q:\n",
    "    # Create a new circuit for the modified version\n",
    "    modified_qc_before = QuantumCircuit(qreg_q, creg_c)\n",
    "    \n",
    "   \n",
    "    # Keep track of whether a NOT gate has been inserted before and after the current qubit\n",
    "    not_before_inserted = False\n",
    "\n",
    "    # Iterate over each gate in the original circuit\n",
    "    for instruction, qargs, cargs in mod_qc:\n",
    "        \n",
    "        # Check if the gate involves the current qubit\n",
    "        if qubit not in qargs:\n",
    "            # Insert a NOT gate before the gate if not already inserted\n",
    "            if not_before_inserted:\n",
    "                modified_qc_before.append(instruction, qargs, cargs)  \n",
    "            else:\n",
    "                modified_qc_before.x(qubit)\n",
    "                modified_qc_before.append(instruction, qargs, cargs)\n",
    "                not_before_inserted = True                       \n",
    "        else:\n",
    "            not_before_inserted=False\n",
    "            modified_qc_before.append(instruction, qargs, cargs)        \n",
    "            \n",
    "    #new_depth_b =  modified_qc_before.depth()\n",
    "    # Add the modified circuits to the list\n",
    "    #if new_depth_b == layers:\n",
    "        modified_circuits_before.append(modified_qc_before)\n",
    "            \n",
    "# Define a helper function to check if a qubit already has a NOT gate in a circuit\n",
    "def has_not_gate(circuit, qubit):\n",
    "    for gate, qargs, _ in circuit.data:\n",
    "        if gate.name == 'x' and qargs[0] == qubit:\n",
    "            return True\n",
    "    return False\n",
    "\n",
    "# Initialize a list to store the modified circuits\n",
    "new_circuits1 = []\n",
    "\n",
    "for modified_qc_before in modified_circuits_before:\n",
    "    # Get the indices of the NOT gates\n",
    "    not_gate_indices = [i for i, (instruction, _, _) in enumerate(modified_qc_before) if instruction.name == 'x']\n",
    "    \n",
    "    # Iterate over the instructions in the circuit\n",
    "    for i, (instruction, qargs, cargs) in enumerate(modified_qc_before):\n",
    "        if instruction.name == 'x' and i in not_gate_indices:\n",
    "            # If the NOT gate is present in cal_qc for this qubit, continue without creating a new circuit\n",
    "            if has_not_gate(cal_qc, qargs[0]):\n",
    "                continue\n",
    "\n",
    "            # Create a new circuit for each NOT gate\n",
    "            new_qc = modified_qc_before.copy()\n",
    "            # Add the corresponding NOT gate to the new circuit\n",
    "            new_qc.x(qargs[0])\n",
    "\n",
    "            # Add the remaining gates from the original circuit to the new circuit\n",
    "            for j in reversed(range(len(new_qc))):\n",
    "                if new_qc[j][0].name == 'x' and j != i and not has_not_gate(cal_qc, new_qc[j][1][0]):\n",
    "                    new_qc.data.pop(j)\n",
    "\n",
    "            new_depth_b = new_qc.depth()\n",
    "            # Add the modified circuits to the list\n",
    "            if new_depth_b == layers:\n",
    "                new_circuits1.append(new_qc)\n",
    "\n",
    "# Rest of the code remains the same\n",
    "\n",
    "\n",
    "# Rest of the code remains the same\n",
    "\n",
    "           \n",
    "\n",
    "def remove_duplicates(circuit_list):\n",
    "    # Initialize an empty set to store the unique circuits\n",
    "    unique_circuits = set()\n",
    "    \n",
    "    # Initialize a new list to store the circuits without duplicates\n",
    "    new_circuit_list = []\n",
    "    \n",
    "    # Iterate over each circuit in the original list\n",
    "    for circuit in circuit_list:\n",
    "        # Convert the circuit to a custom string representation\n",
    "        circuit_str = str(circuit)\n",
    "        \n",
    "        # Check if the circuit string is already in the set\n",
    "        if circuit_str not in unique_circuits:\n",
    "            # If not, add it to the set and the new list\n",
    "            unique_circuits.add(circuit_str)\n",
    "            new_circuit_list.append(circuit)\n",
    "    \n",
    "    return new_circuit_list\n",
    "\n",
    "# Your list of circuits (assuming it is already defined)\n",
    "# circuit_list = [circuit1, circuit2, circuit3, ...]\n",
    "\n",
    "# Call the remove_duplicates function to get the new list without duplicates\n",
    "new_circuit_list = remove_duplicates(new_circuits1)        \n",
    "\n",
    "\n",
    "# Display and save the modified circuits\n",
    "#for i, modified_qc_before in enumerate(modified_circuits_before):\n",
    "    #print(f\"Modified Circuit (before) for qubit {i}:\")\n",
    "    #modified_qc_before.draw(output='mpl', filename=f\"not_before_qubit_{i}.png\")\n",
    "    #print()\n",
    "\n",
    "for i, new_qc in enumerate(new_circuit_list):\n",
    "    print(f\"Modified Circuit (before) for qubit {i}:\")\n",
    "    new_qc.draw(output='mpl', filename=f\"not_qubit_{i}.png\")\n",
    "    #depthfinal = new_qc.depth()\n",
    "    #print(depthfinal)\n",
    "    print()\n",
    "\n",
    "\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "from itertools import combinations\n",
    "new_circuits2 = []         \n",
    "\n",
    "for modified_qc_before in modified_circuits_before:\n",
    "    # Get the indices of the NOT gates\n",
    "    not_gate_indices = [i for i, (instruction, _, _) in enumerate(modified_qc_before) if instruction.name == 'x']\n",
    "\n",
    "    # Iterate over the instructions in the circuit\n",
    "    for combination in combinations(not_gate_indices, 2):\n",
    "        new_qc = modified_qc_before.copy()                             \n",
    "            \n",
    "        # Add the remaining gates from the original circuit to the new circuit\n",
    "        for j in reversed(range(len(new_qc))):\n",
    "            if j not in combination and new_qc[j][0].name == 'x':\n",
    "                new_qc.data.pop(j)\n",
    "        \n",
    "        # Add the modified circuit to the list\n",
    "        new_circuits2.append(new_qc)\n",
    "\n",
    "\n",
    "for i, new_qc in enumerate(new_circuits2):\n",
    "    print(f\"Modified Circuit (before) for qubit {i}:\")\n",
    "    new_qc.draw(output='mpl', filename=f\"not_{i}.png\")\n",
    "    print()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Modified Circuit new for qubit 0:\n",
      "\n",
      "Modified Circuit new for qubit 1:\n",
      "\n",
      "Modified Circuit new for qubit 2:\n",
      "\n",
      "Modified Circuit new for qubit 3:\n",
      "\n",
      "Modified Circuit new for qubit 4:\n",
      "\n"
     ]
    }
   ],
   "source": [
    "from qiskit import QuantumCircuit\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "# Generate circuits with an extra NOT gate from each circuit in modified_circuits_before\n",
    "modified_cal_circuits = []\n",
    "for modified_qc_before in modified_circuits_before:\n",
    "    # Create a new circuit by combining modified_circuit with cal_qc\n",
    "    modified_cal_circuit = cal_qc.copy()\n",
    "    \n",
    "    # Add the extra NOT gate from modified_circuit to modified_cal_circuit\n",
    "    for instruction, qargs, cargs in modified_qc_before:\n",
    "        if instruction.name == 'x':\n",
    "            modified_cal_circuit.x(qargs[0])\n",
    "    \n",
    "    # Add the modified_cal_circuit to the list\n",
    "    modified_cal_circuits.append(modified_cal_circuit)\n",
    "\n",
    "\n",
    "\n",
    "for i, modified_cal_circuit in enumerate(modified_cal_circuits):\n",
    "    print(f\"Modified Circuit new for qubit {i}:\")\n",
    "    modified_cal_circuit.draw(output='mpl', filename=f\"new_{i}.png\")\n",
    "    print()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Input:0\n",
      "{'0000': 916, '0001': 84}\n"
     ]
    },
    {
     "ename": "KeyError",
     "evalue": "1",
     "output_type": "error",
     "traceback": [
      "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[1;31mKeyError\u001b[0m                                  Traceback (most recent call last)",
      "Cell \u001b[1;32mIn[14], line 22\u001b[0m\n\u001b[0;32m     19\u001b[0m count \u001b[39m=\u001b[39m \u001b[39mdict\u001b[39m(\u001b[39msorted\u001b[39m(count\u001b[39m.\u001b[39mitems()))\n\u001b[0;32m     21\u001b[0m \u001b[39mprint\u001b[39m(count)\n\u001b[1;32m---> 22\u001b[0m \u001b[39mprint\u001b[39m(count[\u001b[39m1\u001b[39;49m])\n\u001b[0;32m     24\u001b[0m \u001b[39m#draw a histogram\u001b[39;00m\n",
      "\u001b[1;31mKeyError\u001b[0m: 1"
     ]
    }
   ],
   "source": [
    "\n",
    "cal_out_qc = cal_qc.compose(out_qc, front=False, inplace=False)\n",
    "cal_out_qc.draw(output='mpl')\n",
    "## in_qc + cal_out_qc= org_qc\n",
    "for i, in_qc in enumerate(in_circuits):\n",
    "    print(f\"Input:{i}\")\n",
    "    org_qc = cal_out_qc.compose(in_qc, front=True, inplace=False)\n",
    "    org_qc.draw(output='mpl')\n",
    "    # Run the circuit and get results\n",
    "    shots = 1000    \n",
    "    backend = FakeValencia()\n",
    "    qc = transpile(org_qc, backend=backend, initial_layout= [0,1,2,3,4], optimization_level=1)\n",
    "    results = backend.run(qc,shots=shots).result()\n",
    "    \n",
    "    # Extract the counts from the results\n",
    "    count = results.get_counts(qc)\n",
    "\n",
    "    #sort the counts\n",
    "\n",
    "    count = dict(sorted(count.items()))\n",
    "\n",
    "    print(count)\n",
    "    print(count[1])\n",
    "\n",
    "    #draw a histogram\n",
    "\n",
    "    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Infected circuit 0\n",
      "Input:0\n",
      "Shots: 1000:\n",
      "{'0': 899, '1': 101}\n",
      "Shots: 1100:\n",
      "{'0': 991, '1': 109}\n",
      "Shots: 1200:\n",
      "{'0': 1074, '1': 126}\n",
      "Shots: 1300:\n",
      "{'0': 1174, '1': 126}\n",
      "Shots: 1400:\n",
      "{'0': 1236, '1': 164}\n",
      "Shots: 1500:\n",
      "{'0': 1336, '1': 164}\n",
      "Input:1\n",
      "Shots: 1000:\n",
      "{'0': 897, '1': 103}\n",
      "Shots: 1100:\n",
      "{'0': 964, '1': 136}\n",
      "Shots: 1200:\n",
      "{'0': 1062, '1': 138}\n",
      "Shots: 1300:\n",
      "{'0': 1129, '1': 171}\n",
      "Shots: 1400:\n",
      "{'0': 1255, '1': 145}\n",
      "Shots: 1500:\n",
      "{'0': 1334, '1': 166}\n",
      "Input:2\n",
      "Shots: 1000:\n",
      "{'0': 898, '1': 102}\n",
      "Shots: 1100:\n",
      "{'0': 958, '1': 142}\n",
      "Shots: 1200:\n",
      "{'0': 1066, '1': 134}\n",
      "Shots: 1300:\n",
      "{'0': 1172, '1': 128}\n",
      "Shots: 1400:\n",
      "{'0': 1230, '1': 170}\n",
      "Shots: 1500:\n",
      "{'0': 1328, '1': 172}\n",
      "Input:3\n",
      "Shots: 1000:\n",
      "{'0': 876, '1': 124}\n",
      "Shots: 1100:\n",
      "{'0': 969, '1': 131}\n",
      "Shots: 1200:\n",
      "{'0': 1046, '1': 154}\n",
      "Shots: 1300:\n",
      "{'0': 1153, '1': 147}\n",
      "Shots: 1400:\n",
      "{'0': 1225, '1': 175}\n",
      "Shots: 1500:\n",
      "{'0': 1321, '1': 179}\n",
      "Input:4\n",
      "Shots: 1000:\n",
      "{'0': 896, '1': 104}\n",
      "Shots: 1100:\n",
      "{'0': 970, '1': 130}\n",
      "Shots: 1200:\n",
      "{'0': 1082, '1': 118}\n",
      "Shots: 1300:\n",
      "{'0': 1117, '1': 183}\n",
      "Shots: 1400:\n",
      "{'0': 1248, '1': 152}\n",
      "Shots: 1500:\n",
      "{'0': 1302, '1': 198}\n",
      "Input:5\n",
      "Shots: 1000:\n",
      "{'0': 888, '1': 112}\n",
      "Shots: 1100:\n",
      "{'0': 948, '1': 152}\n",
      "Shots: 1200:\n",
      "{'0': 1051, '1': 149}\n",
      "Shots: 1300:\n",
      "{'0': 1119, '1': 181}\n",
      "Shots: 1400:\n",
      "{'0': 1230, '1': 170}\n",
      "Shots: 1500:\n",
      "{'0': 1329, '1': 171}\n",
      "Input:6\n",
      "Shots: 1000:\n",
      "{'0': 898, '1': 102}\n",
      "Shots: 1100:\n",
      "{'0': 951, '1': 149}\n",
      "Shots: 1200:\n",
      "{'0': 1059, '1': 141}\n",
      "Shots: 1300:\n",
      "{'0': 1139, '1': 161}\n",
      "Shots: 1400:\n",
      "{'0': 1235, '1': 165}\n",
      "Shots: 1500:\n",
      "{'0': 1286, '1': 214}\n",
      "Input:7\n",
      "Shots: 1000:\n",
      "{'0': 881, '1': 119}\n",
      "Shots: 1100:\n",
      "{'0': 972, '1': 128}\n",
      "Shots: 1200:\n",
      "{'0': 1046, '1': 154}\n",
      "Shots: 1300:\n",
      "{'0': 1120, '1': 180}\n",
      "Shots: 1400:\n",
      "{'0': 1219, '1': 181}\n",
      "Shots: 1500:\n",
      "{'0': 1295, '1': 205}\n",
      "Input:8\n",
      "Shots: 1000:\n",
      "{'0': 183, '1': 817}\n",
      "Shots: 1100:\n",
      "{'0': 204, '1': 896}\n",
      "Shots: 1200:\n",
      "{'0': 249, '1': 951}\n",
      "Shots: 1300:\n",
      "{'0': 269, '1': 1031}\n",
      "Shots: 1400:\n",
      "{'0': 246, '1': 1154}\n",
      "Shots: 1500:\n",
      "{'0': 292, '1': 1208}\n",
      "Input:9\n",
      "Shots: 1000:\n",
      "{'0': 200, '1': 800}\n",
      "Shots: 1100:\n",
      "{'0': 227, '1': 873}\n",
      "Shots: 1200:\n",
      "{'0': 235, '1': 965}\n",
      "Shots: 1300:\n",
      "{'0': 270, '1': 1030}\n",
      "Shots: 1400:\n",
      "{'0': 254, '1': 1146}\n",
      "Shots: 1500:\n",
      "{'0': 286, '1': 1214}\n",
      "Input:10\n",
      "Shots: 1000:\n",
      "{'0': 194, '1': 806}\n",
      "Shots: 1100:\n",
      "{'0': 234, '1': 866}\n",
      "Shots: 1200:\n",
      "{'0': 268, '1': 932}\n",
      "Shots: 1300:\n",
      "{'0': 255, '1': 1045}\n",
      "Shots: 1400:\n",
      "{'0': 291, '1': 1109}\n",
      "Shots: 1500:\n",
      "{'0': 263, '1': 1237}\n",
      "Input:11\n",
      "Shots: 1000:\n",
      "{'0': 206, '1': 794}\n",
      "Shots: 1100:\n",
      "{'0': 204, '1': 896}\n",
      "Shots: 1200:\n",
      "{'0': 252, '1': 948}\n",
      "Shots: 1300:\n",
      "{'0': 246, '1': 1054}\n",
      "Shots: 1400:\n",
      "{'0': 286, '1': 1114}\n",
      "Shots: 1500:\n",
      "{'0': 304, '1': 1196}\n",
      "Input:12\n",
      "Shots: 1000:\n",
      "{'0': 810, '1': 190}\n",
      "Shots: 1100:\n",
      "{'0': 896, '1': 204}\n",
      "Shots: 1200:\n",
      "{'0': 964, '1': 236}\n",
      "Shots: 1300:\n",
      "{'0': 1044, '1': 256}\n",
      "Shots: 1400:\n",
      "{'0': 1163, '1': 237}\n",
      "Shots: 1500:\n",
      "{'0': 1227, '1': 273}\n",
      "Input:13\n",
      "Shots: 1000:\n",
      "{'0': 807, '1': 193}\n",
      "Shots: 1100:\n",
      "{'0': 879, '1': 221}\n",
      "Shots: 1200:\n",
      "{'0': 989, '1': 211}\n",
      "Shots: 1300:\n",
      "{'0': 1068, '1': 232}\n",
      "Shots: 1400:\n",
      "{'0': 1155, '1': 245}\n",
      "Shots: 1500:\n",
      "{'0': 1229, '1': 271}\n",
      "Input:14\n",
      "Shots: 1000:\n",
      "{'0': 192, '1': 808}\n",
      "Shots: 1100:\n",
      "{'0': 171, '1': 929}\n",
      "Shots: 1200:\n",
      "{'0': 266, '1': 934}\n",
      "Shots: 1300:\n",
      "{'0': 260, '1': 1040}\n",
      "Shots: 1400:\n",
      "{'0': 289, '1': 1111}\n",
      "Shots: 1500:\n",
      "{'0': 277, '1': 1223}\n",
      "Input:15\n",
      "Shots: 1000:\n",
      "{'0': 211, '1': 789}\n",
      "Shots: 1100:\n",
      "{'0': 206, '1': 894}\n",
      "Shots: 1200:\n",
      "{'0': 232, '1': 968}\n",
      "Shots: 1300:\n",
      "{'0': 280, '1': 1020}\n",
      "Shots: 1400:\n",
      "{'0': 262, '1': 1138}\n",
      "Shots: 1500:\n",
      "{'0': 279, '1': 1221}\n",
      "Infected circuit 1\n",
      "Input:0\n",
      "Shots: 1000:\n",
      "{'0': 897, '1': 103}\n",
      "Shots: 1100:\n",
      "{'0': 990, '1': 110}\n",
      "Shots: 1200:\n",
      "{'0': 1090, '1': 110}\n",
      "Shots: 1300:\n",
      "{'0': 1176, '1': 124}\n",
      "Shots: 1400:\n",
      "{'0': 1241, '1': 159}\n",
      "Shots: 1500:\n",
      "{'0': 1323, '1': 177}\n",
      "Input:1\n",
      "Shots: 1000:\n",
      "{'0': 886, '1': 114}\n",
      "Shots: 1100:\n",
      "{'0': 961, '1': 139}\n",
      "Shots: 1200:\n",
      "{'0': 1062, '1': 138}\n",
      "Shots: 1300:\n",
      "{'0': 1145, '1': 155}\n",
      "Shots: 1400:\n",
      "{'0': 1229, '1': 171}\n",
      "Shots: 1500:\n",
      "{'0': 1330, '1': 170}\n",
      "Input:2\n",
      "Shots: 1000:\n",
      "{'0': 868, '1': 132}\n",
      "Shots: 1100:\n",
      "{'0': 958, '1': 142}\n",
      "Shots: 1200:\n",
      "{'0': 1027, '1': 173}\n",
      "Shots: 1300:\n",
      "{'0': 1138, '1': 162}\n",
      "Shots: 1400:\n",
      "{'0': 1226, '1': 174}\n",
      "Shots: 1500:\n",
      "{'0': 1323, '1': 177}\n",
      "Input:3\n",
      "Shots: 1000:\n",
      "{'0': 876, '1': 124}\n",
      "Shots: 1100:\n",
      "{'0': 977, '1': 123}\n",
      "Shots: 1200:\n",
      "{'0': 1040, '1': 160}\n",
      "Shots: 1300:\n",
      "{'0': 1150, '1': 150}\n",
      "Shots: 1400:\n",
      "{'0': 1183, '1': 217}\n",
      "Shots: 1500:\n",
      "{'0': 1307, '1': 193}\n",
      "Input:4\n",
      "Shots: 1000:\n",
      "{'0': 882, '1': 118}\n",
      "Shots: 1100:\n",
      "{'0': 971, '1': 129}\n",
      "Shots: 1200:\n",
      "{'0': 1045, '1': 155}\n",
      "Shots: 1300:\n",
      "{'0': 1138, '1': 162}\n",
      "Shots: 1400:\n",
      "{'0': 1222, '1': 178}\n",
      "Shots: 1500:\n",
      "{'0': 1312, '1': 188}\n",
      "Input:5\n",
      "Shots: 1000:\n",
      "{'0': 897, '1': 103}\n",
      "Shots: 1100:\n",
      "{'0': 978, '1': 122}\n",
      "Shots: 1200:\n",
      "{'0': 1050, '1': 150}\n",
      "Shots: 1300:\n",
      "{'0': 1153, '1': 147}\n",
      "Shots: 1400:\n",
      "{'0': 1232, '1': 168}\n",
      "Shots: 1500:\n",
      "{'0': 1310, '1': 190}\n",
      "Input:6\n",
      "Shots: 1000:\n",
      "{'0': 878, '1': 122}\n",
      "Shots: 1100:\n",
      "{'0': 959, '1': 141}\n",
      "Shots: 1200:\n",
      "{'0': 1047, '1': 153}\n",
      "Shots: 1300:\n",
      "{'0': 1135, '1': 165}\n",
      "Shots: 1400:\n",
      "{'0': 1212, '1': 188}\n",
      "Shots: 1500:\n",
      "{'0': 1318, '1': 182}\n",
      "Input:7\n",
      "Shots: 1000:\n",
      "{'0': 856, '1': 144}\n",
      "Shots: 1100:\n",
      "{'0': 926, '1': 174}\n",
      "Shots: 1200:\n",
      "{'0': 1046, '1': 154}\n",
      "Shots: 1300:\n",
      "{'0': 1109, '1': 191}\n",
      "Shots: 1400:\n",
      "{'0': 1187, '1': 213}\n",
      "Shots: 1500:\n",
      "{'0': 1285, '1': 215}\n",
      "Input:8\n",
      "Shots: 1000:\n",
      "{'0': 178, '1': 822}\n",
      "Shots: 1100:\n",
      "{'0': 206, '1': 894}\n",
      "Shots: 1200:\n",
      "{'0': 243, '1': 957}\n",
      "Shots: 1300:\n",
      "{'0': 258, '1': 1042}\n",
      "Shots: 1400:\n",
      "{'0': 253, '1': 1147}\n",
      "Shots: 1500:\n",
      "{'0': 250, '1': 1250}\n",
      "Input:9\n",
      "Shots: 1000:\n",
      "{'0': 186, '1': 814}\n",
      "Shots: 1100:\n",
      "{'0': 211, '1': 889}\n",
      "Shots: 1200:\n",
      "{'0': 249, '1': 951}\n",
      "Shots: 1300:\n",
      "{'0': 255, '1': 1045}\n",
      "Shots: 1400:\n",
      "{'0': 296, '1': 1104}\n",
      "Shots: 1500:\n",
      "{'0': 279, '1': 1221}\n",
      "Input:10\n",
      "Shots: 1000:\n",
      "{'0': 216, '1': 784}\n",
      "Shots: 1100:\n",
      "{'0': 203, '1': 897}\n",
      "Shots: 1200:\n",
      "{'0': 241, '1': 959}\n",
      "Shots: 1300:\n",
      "{'0': 210, '1': 1090}\n",
      "Shots: 1400:\n",
      "{'0': 243, '1': 1157}\n",
      "Shots: 1500:\n",
      "{'0': 299, '1': 1201}\n",
      "Input:11\n",
      "Shots: 1000:\n",
      "{'0': 209, '1': 791}\n",
      "Shots: 1100:\n",
      "{'0': 199, '1': 901}\n",
      "Shots: 1200:\n",
      "{'0': 228, '1': 972}\n",
      "Shots: 1300:\n",
      "{'0': 257, '1': 1043}\n",
      "Shots: 1400:\n",
      "{'0': 257, '1': 1143}\n",
      "Shots: 1500:\n",
      "{'0': 289, '1': 1211}\n",
      "Input:12\n",
      "Shots: 1000:\n",
      "{'0': 824, '1': 176}\n",
      "Shots: 1100:\n",
      "{'0': 899, '1': 201}\n",
      "Shots: 1200:\n",
      "{'0': 988, '1': 212}\n",
      "Shots: 1300:\n",
      "{'0': 1054, '1': 246}\n",
      "Shots: 1400:\n",
      "{'0': 1157, '1': 243}\n",
      "Shots: 1500:\n",
      "{'0': 1216, '1': 284}\n",
      "Input:13\n",
      "Shots: 1000:\n",
      "{'0': 822, '1': 178}\n",
      "Shots: 1100:\n",
      "{'0': 889, '1': 211}\n",
      "Shots: 1200:\n",
      "{'0': 973, '1': 227}\n",
      "Shots: 1300:\n",
      "{'0': 1064, '1': 236}\n",
      "Shots: 1400:\n",
      "{'0': 1151, '1': 249}\n",
      "Shots: 1500:\n",
      "{'0': 1248, '1': 252}\n",
      "Input:14\n",
      "Shots: 1000:\n",
      "{'0': 200, '1': 800}\n",
      "Shots: 1100:\n",
      "{'0': 223, '1': 877}\n",
      "Shots: 1200:\n",
      "{'0': 244, '1': 956}\n",
      "Shots: 1300:\n",
      "{'0': 241, '1': 1059}\n",
      "Shots: 1400:\n",
      "{'0': 261, '1': 1139}\n",
      "Shots: 1500:\n",
      "{'0': 295, '1': 1205}\n",
      "Input:15\n",
      "Shots: 1000:\n",
      "{'0': 199, '1': 801}\n",
      "Shots: 1100:\n",
      "{'0': 218, '1': 882}\n",
      "Shots: 1200:\n",
      "{'0': 264, '1': 936}\n",
      "Shots: 1300:\n",
      "{'0': 259, '1': 1041}\n",
      "Shots: 1400:\n",
      "{'0': 288, '1': 1112}\n",
      "Shots: 1500:\n",
      "{'0': 269, '1': 1231}\n",
      "Infected circuit 2\n",
      "Input:0\n",
      "Shots: 1000:\n",
      "{'0': 912, '1': 88}\n",
      "Shots: 1100:\n",
      "{'0': 963, '1': 137}\n",
      "Shots: 1200:\n",
      "{'0': 1062, '1': 138}\n",
      "Shots: 1300:\n",
      "{'0': 1134, '1': 166}\n",
      "Shots: 1400:\n",
      "{'0': 1256, '1': 144}\n",
      "Shots: 1500:\n",
      "{'0': 1335, '1': 165}\n",
      "Input:1\n",
      "Shots: 1000:\n",
      "{'0': 886, '1': 114}\n",
      "Shots: 1100:\n",
      "{'0': 982, '1': 118}\n",
      "Shots: 1200:\n",
      "{'0': 1039, '1': 161}\n",
      "Shots: 1300:\n",
      "{'0': 1121, '1': 179}\n",
      "Shots: 1400:\n",
      "{'0': 1234, '1': 166}\n",
      "Shots: 1500:\n",
      "{'0': 1321, '1': 179}\n",
      "Input:2\n",
      "Shots: 1000:\n",
      "{'0': 260, '1': 740}\n",
      "Shots: 1100:\n",
      "{'0': 247, '1': 853}\n",
      "Shots: 1200:\n",
      "{'0': 269, '1': 931}\n",
      "Shots: 1300:\n",
      "{'0': 279, '1': 1021}\n",
      "Shots: 1400:\n",
      "{'0': 349, '1': 1051}\n",
      "Shots: 1500:\n",
      "{'0': 349, '1': 1151}\n",
      "Input:3\n",
      "Shots: 1000:\n",
      "{'0': 222, '1': 778}\n",
      "Shots: 1100:\n",
      "{'0': 264, '1': 836}\n",
      "Shots: 1200:\n",
      "{'0': 260, '1': 940}\n",
      "Shots: 1300:\n",
      "{'0': 297, '1': 1003}\n",
      "Shots: 1400:\n",
      "{'0': 323, '1': 1077}\n",
      "Shots: 1500:\n",
      "{'0': 348, '1': 1152}\n",
      "Input:4\n",
      "Shots: 1000:\n",
      "{'0': 893, '1': 107}\n",
      "Shots: 1100:\n",
      "{'0': 997, '1': 103}\n",
      "Shots: 1200:\n",
      "{'0': 1097, '1': 103}\n",
      "Shots: 1300:\n",
      "{'0': 1187, '1': 113}\n",
      "Shots: 1400:\n",
      "{'0': 1282, '1': 118}\n",
      "Shots: 1500:\n",
      "{'0': 1365, '1': 135}\n",
      "Input:5\n",
      "Shots: 1000:\n",
      "{'0': 891, '1': 109}\n",
      "Shots: 1100:\n",
      "{'0': 974, '1': 126}\n",
      "Shots: 1200:\n",
      "{'0': 1042, '1': 158}\n",
      "Shots: 1300:\n",
      "{'0': 1168, '1': 132}\n",
      "Shots: 1400:\n",
      "{'0': 1245, '1': 155}\n",
      "Shots: 1500:\n",
      "{'0': 1316, '1': 184}\n",
      "Input:6\n",
      "Shots: 1000:\n",
      "{'0': 876, '1': 124}\n",
      "Shots: 1100:\n",
      "{'0': 981, '1': 119}\n",
      "Shots: 1200:\n",
      "{'0': 1069, '1': 131}\n",
      "Shots: 1300:\n",
      "{'0': 1127, '1': 173}\n",
      "Shots: 1400:\n",
      "{'0': 1201, '1': 199}\n",
      "Shots: 1500:\n",
      "{'0': 1306, '1': 194}\n",
      "Input:7\n",
      "Shots: 1000:\n",
      "{'0': 858, '1': 142}\n",
      "Shots: 1100:\n",
      "{'0': 972, '1': 128}\n",
      "Shots: 1200:\n",
      "{'0': 1040, '1': 160}\n",
      "Shots: 1300:\n",
      "{'0': 1132, '1': 168}\n",
      "Shots: 1400:\n",
      "{'0': 1200, '1': 200}\n",
      "Shots: 1500:\n",
      "{'0': 1294, '1': 206}\n",
      "Input:8\n",
      "Shots: 1000:\n",
      "{'0': 179, '1': 821}\n",
      "Shots: 1100:\n",
      "{'0': 197, '1': 903}\n",
      "Shots: 1200:\n",
      "{'0': 235, '1': 965}\n",
      "Shots: 1300:\n",
      "{'0': 231, '1': 1069}\n",
      "Shots: 1400:\n",
      "{'0': 255, '1': 1145}\n",
      "Shots: 1500:\n",
      "{'0': 259, '1': 1241}\n",
      "Input:9\n",
      "Shots: 1000:\n",
      "{'0': 218, '1': 782}\n",
      "Shots: 1100:\n",
      "{'0': 194, '1': 906}\n",
      "Shots: 1200:\n",
      "{'0': 215, '1': 985}\n",
      "Shots: 1300:\n",
      "{'0': 241, '1': 1059}\n",
      "Shots: 1400:\n",
      "{'0': 275, '1': 1125}\n",
      "Shots: 1500:\n",
      "{'0': 277, '1': 1223}\n",
      "Input:10\n",
      "Shots: 1000:\n",
      "{'0': 160, '1': 840}\n",
      "Shots: 1100:\n",
      "{'0': 181, '1': 919}\n",
      "Shots: 1200:\n",
      "{'0': 213, '1': 987}\n",
      "Shots: 1300:\n",
      "{'0': 240, '1': 1060}\n",
      "Shots: 1400:\n",
      "{'0': 290, '1': 1110}\n",
      "Shots: 1500:\n",
      "{'0': 260, '1': 1240}\n",
      "Input:11\n",
      "Shots: 1000:\n",
      "{'0': 183, '1': 817}\n",
      "Shots: 1100:\n",
      "{'0': 206, '1': 894}\n",
      "Shots: 1200:\n",
      "{'0': 217, '1': 983}\n",
      "Shots: 1300:\n",
      "{'0': 247, '1': 1053}\n",
      "Shots: 1400:\n",
      "{'0': 278, '1': 1122}\n",
      "Shots: 1500:\n",
      "{'0': 270, '1': 1230}\n",
      "Input:12\n",
      "Shots: 1000:\n",
      "{'0': 198, '1': 802}\n",
      "Shots: 1100:\n",
      "{'0': 229, '1': 871}\n",
      "Shots: 1200:\n",
      "{'0': 216, '1': 984}\n",
      "Shots: 1300:\n",
      "{'0': 215, '1': 1085}\n",
      "Shots: 1400:\n",
      "{'0': 280, '1': 1120}\n",
      "Shots: 1500:\n",
      "{'0': 289, '1': 1211}\n",
      "Input:13\n",
      "Shots: 1000:\n",
      "{'0': 187, '1': 813}\n",
      "Shots: 1100:\n",
      "{'0': 225, '1': 875}\n",
      "Shots: 1200:\n",
      "{'0': 230, '1': 970}\n",
      "Shots: 1300:\n",
      "{'0': 251, '1': 1049}\n",
      "Shots: 1400:\n",
      "{'0': 250, '1': 1150}\n",
      "Shots: 1500:\n",
      "{'0': 305, '1': 1195}\n",
      "Input:14\n",
      "Shots: 1000:\n",
      "{'0': 196, '1': 804}\n",
      "Shots: 1100:\n",
      "{'0': 195, '1': 905}\n",
      "Shots: 1200:\n",
      "{'0': 247, '1': 953}\n",
      "Shots: 1300:\n",
      "{'0': 247, '1': 1053}\n",
      "Shots: 1400:\n",
      "{'0': 264, '1': 1136}\n",
      "Shots: 1500:\n",
      "{'0': 264, '1': 1236}\n",
      "Input:15\n",
      "Shots: 1000:\n",
      "{'0': 204, '1': 796}\n",
      "Shots: 1100:\n",
      "{'0': 208, '1': 892}\n",
      "Shots: 1200:\n",
      "{'0': 238, '1': 962}\n",
      "Shots: 1300:\n",
      "{'0': 277, '1': 1023}\n",
      "Shots: 1400:\n",
      "{'0': 256, '1': 1144}\n",
      "Shots: 1500:\n",
      "{'0': 258, '1': 1242}\n"
     ]
    }
   ],
   "source": [
    "for i, new_qc in enumerate(new_circuits1):\n",
    "    print(f\"Infected circuit {i}\")\n",
    "    ## cal_qc + out_qc = cal_out_qc\n",
    "    new_qc.draw(output='mpl')\n",
    "    cal_out_qc = new_qc.compose(out_qc, front=False, inplace=False)\n",
    "    cal_out_qc.draw(output='mpl')\n",
    "    ## in_qc + cal_out_qc= org_qc\n",
    "    for i, in_qc in enumerate(in_circuits):\n",
    "        print(f\"Input:{i}\")\n",
    "        org_qc = cal_out_qc.compose(in_qc, front=True, inplace=False)\n",
    "        org_qc.draw(output='mpl')\n",
    "        # Run the circuit and get results\n",
    "        shots = [i for i in range(1000, 1501, 100)]\n",
    "        for i in shots:\n",
    "            print(f\"Shots: {i}:\")    \n",
    "            backend = FakeValencia()\n",
    "            qc_noisy = transpile(org_qc, backend=backend, initial_layout= [0,1,2,3,4], optimization_level=1)\n",
    "            results = backend.run(qc_noisy,shots=i).result()\n",
    "            \n",
    "            # Extract the counts from the results\n",
    "            count_noisy = results.get_counts(qc_noisy)\n",
    "\n",
    "            #sort the counts\n",
    "\n",
    "            count_noisy = dict(sorted(count_noisy.items()))\n",
    "\n",
    "            print(count_noisy)\n",
    "\n",
    "            #draw a histogram\n",
    "\n",
    "            "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Infected circuit 0\n",
      "Input:0\n",
      "{'0': 668, '1': 332}\n",
      "{'0': 644, '1': 356}\n",
      "diff=-24\n",
      "Shots: 1100:\n",
      "{'0': 653, '1': 347}\n",
      "Shots: 1200:\n",
      "{'0': 661, '1': 339}\n",
      "Shots: 1300:\n",
      "{'0': 633, '1': 367}\n",
      "Shots: 1400:\n",
      "{'0': 628, '1': 372}\n",
      "Shots: 1500:\n",
      "{'0': 654, '1': 346}\n",
      "Input:1\n",
      "{'0': 624, '1': 376}\n",
      "{'0': 426, '1': 574}\n",
      "Shots: 1100:\n",
      "{'0': 434, '1': 566}\n",
      "Shots: 1200:\n",
      "{'0': 385, '1': 615}\n",
      "Shots: 1300:\n",
      "{'0': 410, '1': 590}\n",
      "Shots: 1400:\n",
      "{'0': 420, '1': 580}\n",
      "Shots: 1500:\n",
      "{'0': 443, '1': 557}\n",
      "Input:2\n",
      "{'0': 648, '1': 352}\n",
      "{'0': 662, '1': 338}\n",
      "diff=14\n",
      "Input:3\n",
      "{'0': 613, '1': 387}\n",
      "{'0': 397, '1': 603}\n",
      "Input:4\n",
      "{'0': 643, '1': 357}\n",
      "{'0': 655, '1': 345}\n",
      "diff=12\n",
      "Input:5\n",
      "{'0': 424, '1': 576}\n",
      "{'0': 611, '1': 389}\n",
      "Input:6\n",
      "{'0': 481, '1': 519}\n",
      "{'0': 662, '1': 338}\n",
      "Input:7\n",
      "{'0': 401, '1': 599}\n",
      "{'0': 574, '1': 426}\n",
      "Input:8\n",
      "{'0': 443, '1': 557}\n",
      "{'0': 647, '1': 353}\n",
      "Input:9\n",
      "{'0': 373, '1': 627}\n",
      "{'0': 399, '1': 601}\n",
      "diff=-26\n",
      "Shots: 1100:\n",
      "{'0': 396, '1': 604}\n",
      "Shots: 1200:\n",
      "{'0': 374, '1': 626}\n",
      "Shots: 1300:\n",
      "{'0': 390, '1': 610}\n",
      "Shots: 1400:\n",
      "{'0': 400, '1': 600}\n",
      "Shots: 1500:\n",
      "{'0': 419, '1': 581}\n",
      "Input:10\n",
      "{'0': 403, '1': 597}\n",
      "{'0': 623, '1': 377}\n",
      "Shots: 1100:\n",
      "{'0': 646, '1': 354}\n",
      "Shots: 1200:\n",
      "{'0': 647, '1': 353}\n",
      "Shots: 1300:\n",
      "{'0': 624, '1': 376}\n",
      "Shots: 1400:\n",
      "{'0': 645, '1': 355}\n",
      "Shots: 1500:\n",
      "{'0': 640, '1': 360}\n",
      "Input:11\n",
      "{'0': 427, '1': 573}\n",
      "{'0': 426, '1': 574}\n",
      "diff=1\n",
      "Input:12\n",
      "{'0': 448, '1': 552}\n",
      "{'0': 659, '1': 341}\n",
      "Input:13\n",
      "{'0': 383, '1': 617}\n",
      "{'0': 424, '1': 576}\n",
      "diff=-41\n",
      "Shots: 1100:\n",
      "{'0': 405, '1': 595}\n",
      "Shots: 1200:\n",
      "{'0': 407, '1': 593}\n",
      "Shots: 1300:\n",
      "{'0': 438, '1': 562}\n",
      "Shots: 1400:\n",
      "{'0': 418, '1': 582}\n",
      "Shots: 1500:\n",
      "{'0': 395, '1': 605}\n",
      "Input:14\n",
      "{'0': 485, '1': 515}\n",
      "{'0': 623, '1': 377}\n",
      "Shots: 1100:\n",
      "{'0': 604, '1': 396}\n",
      "Shots: 1200:\n",
      "{'0': 632, '1': 368}\n",
      "Shots: 1300:\n",
      "{'0': 616, '1': 384}\n",
      "Shots: 1400:\n",
      "{'0': 635, '1': 365}\n",
      "Shots: 1500:\n",
      "{'0': 618, '1': 382}\n",
      "Input:15\n",
      "{'0': 411, '1': 589}\n",
      "{'0': 581, '1': 419}\n",
      "Shots: 1100:\n",
      "{'0': 574, '1': 426}\n",
      "Shots: 1200:\n",
      "{'0': 605, '1': 395}\n",
      "Shots: 1300:\n",
      "{'0': 592, '1': 408}\n",
      "Shots: 1400:\n",
      "{'0': 616, '1': 384}\n",
      "Shots: 1500:\n",
      "{'0': 610, '1': 390}\n",
      "Infected circuit 1\n",
      "Input:0\n",
      "{'0': 651, '1': 349}\n",
      "{'0': 415, '1': 585}\n",
      "Shots: 1100:\n",
      "{'0': 393, '1': 607}\n",
      "Shots: 1200:\n",
      "{'0': 395, '1': 605}\n",
      "Shots: 1300:\n",
      "{'0': 438, '1': 562}\n",
      "Shots: 1400:\n",
      "{'0': 435, '1': 565}\n",
      "Shots: 1500:\n",
      "{'0': 440, '1': 560}\n",
      "Input:1\n",
      "{'0': 606, '1': 394}\n",
      "{'0': 401, '1': 599}\n",
      "Shots: 1100:\n",
      "{'0': 373, '1': 627}\n",
      "Shots: 1200:\n",
      "{'0': 403, '1': 597}\n",
      "Shots: 1300:\n",
      "{'0': 403, '1': 597}\n",
      "Shots: 1400:\n",
      "{'0': 380, '1': 620}\n",
      "Shots: 1500:\n",
      "{'0': 413, '1': 587}\n",
      "Input:2\n",
      "{'0': 642, '1': 358}\n",
      "{'0': 432, '1': 568}\n",
      "Shots: 1100:\n",
      "{'0': 431, '1': 569}\n",
      "Shots: 1200:\n",
      "{'0': 423, '1': 577}\n",
      "Shots: 1300:\n",
      "{'0': 437, '1': 563}\n",
      "Shots: 1400:\n",
      "{'0': 463, '1': 537}\n",
      "Shots: 1500:\n",
      "{'0': 384, '1': 616}\n",
      "Input:3\n",
      "{'0': 595, '1': 405}\n",
      "{'0': 396, '1': 604}\n",
      "Shots: 1100:\n",
      "{'0': 402, '1': 598}\n",
      "Shots: 1200:\n",
      "{'0': 388, '1': 612}\n",
      "Shots: 1300:\n",
      "{'0': 421, '1': 579}\n",
      "Shots: 1400:\n",
      "{'0': 400, '1': 600}\n",
      "Shots: 1500:\n",
      "{'0': 394, '1': 606}\n",
      "Input:4\n",
      "{'0': 626, '1': 374}\n",
      "{'0': 456, '1': 544}\n",
      "Shots: 1100:\n",
      "{'0': 469, '1': 531}\n",
      "Shots: 1200:\n",
      "{'0': 457, '1': 543}\n",
      "Shots: 1300:\n",
      "{'0': 445, '1': 555}\n",
      "Shots: 1400:\n",
      "{'0': 449, '1': 551}\n",
      "Shots: 1500:\n",
      "{'0': 457, '1': 543}\n",
      "Input:5\n",
      "{'0': 411, '1': 589}\n",
      "{'0': 409, '1': 591}\n",
      "diff=2\n",
      "Input:6\n",
      "{'0': 432, '1': 568}\n",
      "{'0': 451, '1': 549}\n",
      "diff=-19\n",
      "Shots: 1100:\n",
      "{'0': 451, '1': 549}\n",
      "Shots: 1200:\n",
      "{'0': 462, '1': 538}\n",
      "Shots: 1300:\n",
      "{'0': 449, '1': 551}\n",
      "Shots: 1400:\n",
      "{'0': 445, '1': 555}\n",
      "Shots: 1500:\n",
      "{'0': 471, '1': 529}\n",
      "Input:7\n",
      "{'0': 406, '1': 594}\n",
      "{'0': 448, '1': 552}\n",
      "diff=-42\n",
      "Shots: 1100:\n",
      "{'0': 421, '1': 579}\n",
      "Shots: 1200:\n",
      "{'0': 439, '1': 561}\n",
      "Shots: 1300:\n",
      "{'0': 426, '1': 574}\n",
      "Shots: 1400:\n",
      "{'0': 428, '1': 572}\n",
      "Shots: 1500:\n",
      "{'0': 417, '1': 583}\n",
      "Input:8\n",
      "{'0': 409, '1': 591}\n",
      "{'0': 659, '1': 341}\n",
      "Shots: 1100:\n",
      "{'0': 645, '1': 355}\n",
      "Shots: 1200:\n",
      "{'0': 654, '1': 346}\n",
      "Shots: 1300:\n",
      "{'0': 648, '1': 352}\n",
      "Shots: 1400:\n",
      "{'0': 643, '1': 357}\n",
      "Shots: 1500:\n",
      "{'0': 640, '1': 360}\n",
      "Input:9\n",
      "{'0': 390, '1': 610}\n",
      "{'0': 642, '1': 358}\n",
      "Shots: 1100:\n",
      "{'0': 627, '1': 373}\n",
      "Shots: 1200:\n",
      "{'0': 630, '1': 370}\n",
      "Shots: 1300:\n",
      "{'0': 615, '1': 385}\n",
      "Shots: 1400:\n",
      "{'0': 632, '1': 368}\n",
      "Shots: 1500:\n",
      "{'0': 603, '1': 397}\n",
      "Input:10\n",
      "{'0': 437, '1': 563}\n",
      "{'0': 630, '1': 370}\n",
      "Shots: 1100:\n",
      "{'0': 643, '1': 357}\n",
      "Shots: 1200:\n",
      "{'0': 643, '1': 357}\n",
      "Shots: 1300:\n",
      "{'0': 625, '1': 375}\n",
      "Shots: 1400:\n",
      "{'0': 646, '1': 354}\n",
      "Shots: 1500:\n",
      "{'0': 640, '1': 360}\n",
      "Input:11\n",
      "{'0': 405, '1': 595}\n",
      "{'0': 624, '1': 376}\n",
      "Shots: 1100:\n",
      "{'0': 617, '1': 383}\n",
      "Shots: 1200:\n",
      "{'0': 617, '1': 383}\n",
      "Shots: 1300:\n",
      "{'0': 630, '1': 370}\n",
      "Shots: 1400:\n",
      "{'0': 590, '1': 410}\n",
      "Shots: 1500:\n",
      "{'0': 603, '1': 397}\n",
      "Input:12\n",
      "{'0': 450, '1': 550}\n",
      "{'0': 644, '1': 356}\n",
      "Shots: 1100:\n",
      "{'0': 653, '1': 347}\n",
      "Shots: 1200:\n",
      "{'0': 625, '1': 375}\n",
      "Shots: 1300:\n",
      "{'0': 633, '1': 367}\n",
      "Shots: 1400:\n",
      "{'0': 644, '1': 356}\n",
      "Shots: 1500:\n",
      "{'0': 607, '1': 393}\n",
      "Input:13\n",
      "{'0': 399, '1': 601}\n",
      "{'0': 430, '1': 570}\n",
      "diff=-31\n",
      "Shots: 1100:\n",
      "{'0': 406, '1': 594}\n",
      "Shots: 1200:\n",
      "{'0': 425, '1': 575}\n",
      "Shots: 1300:\n",
      "{'0': 438, '1': 562}\n",
      "Shots: 1400:\n",
      "{'0': 430, '1': 570}\n",
      "Shots: 1500:\n",
      "{'0': 421, '1': 579}\n",
      "Input:14\n",
      "{'0': 499, '1': 501}\n",
      "{'0': 463, '1': 537}\n",
      "diff=36\n",
      "Input:15\n",
      "{'0': 448, '1': 552}\n",
      "{'0': 401, '1': 599}\n",
      "diff=47\n",
      "Infected circuit 2\n",
      "Input:0\n",
      "{'0': 637, '1': 363}\n",
      "{'0': 649, '1': 351}\n",
      "diff=12\n",
      "Input:1\n",
      "{'0': 625, '1': 375}\n",
      "{'0': 403, '1': 597}\n",
      "Input:2\n",
      "{'0': 641, '1': 359}\n",
      "{'0': 473, '1': 527}\n",
      "Input:3\n",
      "{'0': 610, '1': 390}\n",
      "{'0': 402, '1': 598}\n",
      "Input:4\n",
      "{'0': 650, '1': 350}\n",
      "{'0': 675, '1': 325}\n",
      "diff=25\n",
      "Input:5\n",
      "{'0': 419, '1': 581}\n",
      "{'0': 614, '1': 386}\n",
      "Input:6\n",
      "{'0': 447, '1': 553}\n",
      "{'0': 637, '1': 363}\n",
      "Input:7\n",
      "{'0': 426, '1': 574}\n",
      "{'0': 616, '1': 384}\n",
      "Input:8\n",
      "{'0': 460, '1': 540}\n",
      "{'0': 439, '1': 561}\n",
      "diff=21\n",
      "Input:9\n",
      "{'0': 389, '1': 611}\n",
      "{'0': 411, '1': 589}\n",
      "diff=-22\n",
      "Shots: 1100:\n",
      "{'0': 388, '1': 612}\n",
      "Shots: 1200:\n",
      "{'0': 415, '1': 585}\n",
      "Shots: 1300:\n",
      "{'0': 411, '1': 589}\n",
      "Shots: 1400:\n",
      "{'0': 436, '1': 564}\n",
      "Shots: 1500:\n",
      "{'0': 375, '1': 625}\n",
      "Input:10\n",
      "{'0': 450, '1': 550}\n",
      "{'0': 497, '1': 503}\n",
      "diff=-47\n",
      "Shots: 1100:\n",
      "{'0': 469, '1': 531}\n",
      "Shots: 1200:\n",
      "{'0': 476, '1': 524}\n",
      "Shots: 1300:\n",
      "{'0': 451, '1': 549}\n",
      "Shots: 1400:\n",
      "{'0': 458, '1': 542}\n",
      "Shots: 1500:\n",
      "{'0': 450, '1': 550}\n",
      "Input:11\n",
      "{'0': 395, '1': 605}\n",
      "{'0': 414, '1': 586}\n",
      "diff=-19\n",
      "Shots: 1100:\n",
      "{'0': 392, '1': 608}\n",
      "Shots: 1200:\n",
      "{'0': 428, '1': 572}\n",
      "Shots: 1300:\n",
      "{'0': 401, '1': 599}\n",
      "Shots: 1400:\n",
      "{'0': 402, '1': 598}\n",
      "Shots: 1500:\n",
      "{'0': 420, '1': 580}\n",
      "Input:12\n",
      "{'0': 454, '1': 546}\n",
      "{'0': 431, '1': 569}\n",
      "diff=23\n",
      "Input:13\n",
      "{'0': 393, '1': 607}\n",
      "{'0': 412, '1': 588}\n",
      "diff=-19\n",
      "Shots: 1100:\n",
      "{'0': 414, '1': 586}\n",
      "Shots: 1200:\n",
      "{'0': 429, '1': 571}\n",
      "Shots: 1300:\n",
      "{'0': 410, '1': 590}\n",
      "Shots: 1400:\n",
      "{'0': 400, '1': 600}\n",
      "Shots: 1500:\n",
      "{'0': 391, '1': 609}\n",
      "Input:14\n",
      "{'0': 424, '1': 576}\n",
      "{'0': 466, '1': 534}\n",
      "diff=-42\n",
      "Shots: 1100:\n",
      "{'0': 422, '1': 578}\n",
      "Shots: 1200:\n",
      "{'0': 451, '1': 549}\n",
      "Shots: 1300:\n",
      "{'0': 419, '1': 581}\n",
      "Shots: 1400:\n",
      "{'0': 443, '1': 557}\n",
      "Shots: 1500:\n",
      "{'0': 461, '1': 539}\n",
      "Input:15\n",
      "{'0': 417, '1': 583}\n",
      "{'0': 413, '1': 587}\n",
      "diff=4\n",
      "Infected circuit 3\n",
      "Input:0\n",
      "{'0': 645, '1': 355}\n",
      "{'0': 643, '1': 357}\n",
      "diff=-2\n",
      "Shots: 1100:\n",
      "{'0': 651, '1': 349}\n",
      "Shots: 1200:\n",
      "{'0': 645, '1': 355}\n",
      "Shots: 1300:\n",
      "{'0': 641, '1': 359}\n",
      "Shots: 1400:\n",
      "{'0': 649, '1': 351}\n",
      "Shots: 1500:\n",
      "{'0': 619, '1': 381}\n",
      "Input:1\n",
      "{'0': 625, '1': 375}\n",
      "{'0': 593, '1': 407}\n",
      "diff=-32\n",
      "Shots: 1100:\n",
      "{'0': 597, '1': 403}\n",
      "Shots: 1200:\n",
      "{'0': 607, '1': 393}\n",
      "Shots: 1300:\n",
      "{'0': 602, '1': 398}\n",
      "Shots: 1400:\n",
      "{'0': 602, '1': 398}\n",
      "Shots: 1500:\n",
      "{'0': 597, '1': 403}\n",
      "Input:2\n",
      "{'0': 636, '1': 364}\n",
      "{'0': 436, '1': 564}\n",
      "Shots: 1100:\n",
      "{'0': 442, '1': 558}\n",
      "Shots: 1200:\n",
      "{'0': 424, '1': 576}\n",
      "Shots: 1300:\n",
      "{'0': 458, '1': 542}\n",
      "Shots: 1400:\n",
      "{'0': 465, '1': 535}\n",
      "Shots: 1500:\n",
      "{'0': 425, '1': 575}\n",
      "Input:3\n",
      "{'0': 590, '1': 410}\n",
      "{'0': 599, '1': 401}\n",
      "diff=9\n",
      "Input:4\n",
      "{'0': 631, '1': 369}\n",
      "{'0': 641, '1': 359}\n",
      "diff=10\n",
      "Input:5\n",
      "{'0': 405, '1': 595}\n",
      "{'0': 402, '1': 598}\n",
      "diff=3\n",
      "Input:6\n",
      "{'0': 419, '1': 581}\n",
      "{'0': 644, '1': 356}\n",
      "Input:7\n",
      "{'0': 421, '1': 579}\n",
      "{'0': 420, '1': 580}\n",
      "diff=1\n",
      "Input:8\n",
      "{'0': 409, '1': 591}\n",
      "{'0': 431, '1': 569}\n",
      "diff=-22\n",
      "Shots: 1100:\n",
      "{'0': 425, '1': 575}\n",
      "Shots: 1200:\n",
      "{'0': 421, '1': 579}\n",
      "Shots: 1300:\n",
      "{'0': 450, '1': 550}\n",
      "Shots: 1400:\n",
      "{'0': 449, '1': 551}\n",
      "Shots: 1500:\n",
      "{'0': 449, '1': 551}\n",
      "Input:9\n",
      "{'0': 424, '1': 576}\n",
      "{'0': 406, '1': 594}\n",
      "diff=18\n",
      "Input:10\n",
      "{'0': 395, '1': 605}\n",
      "{'0': 471, '1': 529}\n",
      "diff=-76\n",
      "Shots: 1100:\n",
      "{'0': 500, '1': 500}\n",
      "Shots: 1200:\n",
      "{'0': 447, '1': 553}\n",
      "Shots: 1300:\n",
      "{'0': 485, '1': 515}\n",
      "Shots: 1400:\n",
      "{'0': 470, '1': 530}\n",
      "Shots: 1500:\n",
      "{'0': 453, '1': 547}\n",
      "Input:11\n",
      "{'0': 413, '1': 587}\n",
      "{'0': 605, '1': 395}\n",
      "Shots: 1100:\n",
      "{'0': 624, '1': 376}\n",
      "Shots: 1200:\n",
      "{'0': 588, '1': 412}\n",
      "Shots: 1300:\n",
      "{'0': 579, '1': 421}\n",
      "Shots: 1400:\n",
      "{'0': 605, '1': 395}\n",
      "Shots: 1500:\n",
      "{'0': 599, '1': 401}\n",
      "Input:12\n",
      "{'0': 420, '1': 580}\n",
      "{'0': 424, '1': 576}\n",
      "diff=-4\n",
      "Shots: 1100:\n",
      "{'0': 456, '1': 544}\n",
      "Shots: 1200:\n",
      "{'0': 449, '1': 551}\n",
      "Shots: 1300:\n",
      "{'0': 435, '1': 565}\n",
      "Shots: 1400:\n",
      "{'0': 429, '1': 571}\n",
      "Shots: 1500:\n",
      "{'0': 418, '1': 582}\n",
      "Input:13\n",
      "{'0': 401, '1': 599}\n",
      "{'0': 419, '1': 581}\n",
      "diff=-18\n",
      "Shots: 1100:\n",
      "{'0': 400, '1': 600}\n",
      "Shots: 1200:\n",
      "{'0': 392, '1': 608}\n",
      "Shots: 1300:\n",
      "{'0': 416, '1': 584}\n",
      "Shots: 1400:\n",
      "{'0': 404, '1': 596}\n",
      "Shots: 1500:\n",
      "{'0': 415, '1': 585}\n",
      "Input:14\n",
      "{'0': 421, '1': 579}\n",
      "{'0': 432, '1': 568}\n",
      "diff=-11\n",
      "Shots: 1100:\n",
      "{'0': 429, '1': 571}\n",
      "Shots: 1200:\n",
      "{'0': 423, '1': 577}\n",
      "Shots: 1300:\n",
      "{'0': 430, '1': 570}\n",
      "Shots: 1400:\n",
      "{'0': 452, '1': 548}\n",
      "Shots: 1500:\n",
      "{'0': 460, '1': 540}\n",
      "Input:15\n",
      "{'0': 402, '1': 598}\n",
      "{'0': 396, '1': 604}\n",
      "diff=6\n",
      "Infected circuit 4\n",
      "Input:0\n",
      "{'0': 652, '1': 348}\n",
      "{'0': 629, '1': 371}\n",
      "diff=-23\n",
      "Shots: 1100:\n",
      "{'0': 655, '1': 345}\n",
      "Shots: 1200:\n",
      "{'0': 649, '1': 351}\n",
      "Shots: 1300:\n",
      "{'0': 649, '1': 351}\n",
      "Shots: 1400:\n",
      "{'0': 632, '1': 368}\n",
      "Shots: 1500:\n",
      "{'0': 668, '1': 332}\n",
      "Input:1\n",
      "{'0': 619, '1': 381}\n",
      "{'0': 593, '1': 407}\n",
      "diff=-26\n",
      "Shots: 1100:\n",
      "{'0': 581, '1': 419}\n",
      "Shots: 1200:\n",
      "{'0': 616, '1': 384}\n",
      "Shots: 1300:\n",
      "{'0': 605, '1': 395}\n",
      "Shots: 1400:\n",
      "{'0': 603, '1': 397}\n",
      "Shots: 1500:\n",
      "{'0': 595, '1': 405}\n",
      "Input:2\n",
      "{'0': 641, '1': 359}\n",
      "{'0': 628, '1': 372}\n",
      "diff=-13\n",
      "Shots: 1100:\n",
      "{'0': 639, '1': 361}\n",
      "Shots: 1200:\n",
      "{'0': 650, '1': 350}\n",
      "Shots: 1300:\n",
      "{'0': 627, '1': 373}\n",
      "Shots: 1400:\n",
      "{'0': 658, '1': 342}\n",
      "Shots: 1500:\n",
      "{'0': 622, '1': 378}\n",
      "Input:3\n",
      "{'0': 643, '1': 357}\n",
      "{'0': 451, '1': 549}\n",
      "Shots: 1100:\n",
      "{'0': 470, '1': 530}\n",
      "Shots: 1200:\n",
      "{'0': 486, '1': 514}\n",
      "Shots: 1300:\n",
      "{'0': 478, '1': 522}\n",
      "Shots: 1400:\n",
      "{'0': 468, '1': 532}\n",
      "Shots: 1500:\n",
      "{'0': 443, '1': 557}\n",
      "Input:4\n",
      "{'0': 626, '1': 374}\n",
      "{'0': 663, '1': 337}\n",
      "diff=37\n",
      "Input:5\n",
      "{'0': 434, '1': 566}\n",
      "{'0': 431, '1': 569}\n",
      "diff=3\n",
      "Input:6\n",
      "{'0': 453, '1': 547}\n",
      "{'0': 465, '1': 535}\n",
      "diff=-12\n",
      "Shots: 1100:\n",
      "{'0': 431, '1': 569}\n",
      "Shots: 1200:\n",
      "{'0': 448, '1': 552}\n",
      "Shots: 1300:\n",
      "{'0': 463, '1': 537}\n",
      "Shots: 1400:\n",
      "{'0': 462, '1': 538}\n",
      "Shots: 1500:\n",
      "{'0': 432, '1': 568}\n",
      "Input:7\n",
      "{'0': 441, '1': 559}\n",
      "{'0': 418, '1': 582}\n",
      "diff=23\n",
      "Input:8\n",
      "{'0': 415, '1': 585}\n",
      "{'0': 421, '1': 579}\n",
      "diff=-6\n",
      "Shots: 1100:\n",
      "{'0': 435, '1': 565}\n",
      "Shots: 1200:\n",
      "{'0': 434, '1': 566}\n",
      "Shots: 1300:\n",
      "{'0': 433, '1': 567}\n",
      "Shots: 1400:\n",
      "{'0': 444, '1': 556}\n",
      "Shots: 1500:\n",
      "{'0': 449, '1': 551}\n",
      "Input:9\n",
      "{'0': 421, '1': 579}\n",
      "{'0': 433, '1': 567}\n",
      "diff=-12\n",
      "Shots: 1100:\n",
      "{'0': 394, '1': 606}\n",
      "Shots: 1200:\n",
      "{'0': 424, '1': 576}\n",
      "Shots: 1300:\n",
      "{'0': 399, '1': 601}\n",
      "Shots: 1400:\n",
      "{'0': 409, '1': 591}\n",
      "Shots: 1500:\n",
      "{'0': 397, '1': 603}\n",
      "Input:10\n",
      "{'0': 450, '1': 550}\n",
      "{'0': 609, '1': 391}\n",
      "Shots: 1100:\n",
      "{'0': 596, '1': 404}\n",
      "Shots: 1200:\n",
      "{'0': 612, '1': 388}\n",
      "Shots: 1300:\n",
      "{'0': 597, '1': 403}\n",
      "Shots: 1400:\n",
      "{'0': 616, '1': 384}\n",
      "Shots: 1500:\n",
      "{'0': 597, '1': 403}\n",
      "Input:11\n",
      "{'0': 421, '1': 579}\n",
      "{'0': 404, '1': 596}\n",
      "diff=17\n",
      "Input:12\n",
      "{'0': 431, '1': 569}\n",
      "{'0': 454, '1': 546}\n",
      "diff=-23\n",
      "Shots: 1100:\n",
      "{'0': 456, '1': 544}\n",
      "Shots: 1200:\n",
      "{'0': 444, '1': 556}\n",
      "Shots: 1300:\n",
      "{'0': 412, '1': 588}\n",
      "Shots: 1400:\n",
      "{'0': 401, '1': 599}\n",
      "Shots: 1500:\n",
      "{'0': 439, '1': 561}\n",
      "Input:13\n",
      "{'0': 395, '1': 605}\n",
      "{'0': 376, '1': 624}\n",
      "diff=19\n",
      "Input:14\n",
      "{'0': 458, '1': 542}\n",
      "{'0': 622, '1': 378}\n",
      "Input:15\n",
      "{'0': 423, '1': 577}\n",
      "{'0': 416, '1': 584}\n",
      "diff=7\n",
      "Infected circuit 5\n",
      "Input:0\n",
      "{'0': 666, '1': 334}\n",
      "{'0': 631, '1': 369}\n",
      "diff=-35\n",
      "Shots: 1100:\n",
      "{'0': 656, '1': 344}\n",
      "Shots: 1200:\n",
      "{'0': 641, '1': 359}\n",
      "Shots: 1300:\n",
      "{'0': 655, '1': 345}\n",
      "Shots: 1400:\n",
      "{'0': 650, '1': 350}\n",
      "Shots: 1500:\n",
      "{'0': 643, '1': 357}\n",
      "Input:1\n",
      "{'0': 603, '1': 397}\n",
      "{'0': 555, '1': 445}\n",
      "diff=-48\n",
      "Shots: 1100:\n",
      "{'0': 592, '1': 408}\n",
      "Shots: 1200:\n",
      "{'0': 606, '1': 394}\n",
      "Shots: 1300:\n",
      "{'0': 588, '1': 412}\n",
      "Shots: 1400:\n",
      "{'0': 621, '1': 379}\n",
      "Shots: 1500:\n",
      "{'0': 609, '1': 391}\n",
      "Input:2\n",
      "{'0': 660, '1': 340}\n",
      "{'0': 676, '1': 324}\n",
      "diff=16\n",
      "Input:3\n",
      "{'0': 618, '1': 382}\n",
      "{'0': 592, '1': 408}\n",
      "diff=-26\n",
      "Shots: 1100:\n",
      "{'0': 585, '1': 415}\n",
      "Shots: 1200:\n",
      "{'0': 647, '1': 353}\n",
      "Shots: 1300:\n",
      "{'0': 617, '1': 383}\n",
      "Shots: 1400:\n",
      "{'0': 630, '1': 370}\n",
      "Shots: 1500:\n",
      "{'0': 615, '1': 385}\n",
      "Input:4\n",
      "{'0': 637, '1': 363}\n",
      "{'0': 479, '1': 521}\n",
      "Shots: 1100:\n",
      "{'0': 444, '1': 556}\n",
      "Shots: 1200:\n",
      "{'0': 465, '1': 535}\n",
      "Shots: 1300:\n",
      "{'0': 444, '1': 556}\n",
      "Shots: 1400:\n",
      "{'0': 480, '1': 520}\n",
      "Shots: 1500:\n",
      "{'0': 463, '1': 537}\n",
      "Input:5\n",
      "{'0': 420, '1': 580}\n",
      "{'0': 434, '1': 566}\n",
      "diff=-14\n",
      "Shots: 1100:\n",
      "{'0': 423, '1': 577}\n",
      "Shots: 1200:\n",
      "{'0': 423, '1': 577}\n",
      "Shots: 1300:\n",
      "{'0': 416, '1': 584}\n",
      "Shots: 1400:\n",
      "{'0': 425, '1': 575}\n",
      "Shots: 1500:\n",
      "{'0': 435, '1': 565}\n",
      "Input:6\n",
      "{'0': 441, '1': 559}\n",
      "{'0': 632, '1': 368}\n",
      "Shots: 1100:\n",
      "{'0': 625, '1': 375}\n",
      "Shots: 1200:\n",
      "{'0': 622, '1': 378}\n",
      "Shots: 1300:\n",
      "{'0': 637, '1': 363}\n",
      "Shots: 1400:\n",
      "{'0': 632, '1': 368}\n",
      "Shots: 1500:\n",
      "{'0': 642, '1': 358}\n",
      "Input:7\n",
      "{'0': 431, '1': 569}\n",
      "{'0': 419, '1': 581}\n",
      "diff=12\n",
      "Input:8\n",
      "{'0': 461, '1': 539}\n",
      "{'0': 452, '1': 548}\n",
      "diff=9\n",
      "Input:9\n",
      "{'0': 387, '1': 613}\n",
      "{'0': 412, '1': 588}\n",
      "diff=-25\n",
      "Shots: 1100:\n",
      "{'0': 410, '1': 590}\n",
      "Shots: 1200:\n",
      "{'0': 445, '1': 555}\n",
      "Shots: 1300:\n",
      "{'0': 394, '1': 606}\n",
      "Shots: 1400:\n",
      "{'0': 435, '1': 565}\n",
      "Shots: 1500:\n",
      "{'0': 394, '1': 606}\n",
      "Input:10\n",
      "{'0': 437, '1': 563}\n",
      "{'0': 418, '1': 582}\n",
      "diff=19\n",
      "Input:11\n",
      "{'0': 415, '1': 585}\n",
      "{'0': 388, '1': 612}\n",
      "diff=27\n",
      "Input:12\n",
      "{'0': 427, '1': 573}\n",
      "{'0': 453, '1': 547}\n",
      "diff=-26\n",
      "Shots: 1100:\n",
      "{'0': 441, '1': 559}\n",
      "Shots: 1200:\n",
      "{'0': 471, '1': 529}\n",
      "Shots: 1300:\n",
      "{'0': 456, '1': 544}\n",
      "Shots: 1400:\n",
      "{'0': 452, '1': 548}\n",
      "Shots: 1500:\n",
      "{'0': 430, '1': 570}\n",
      "Input:13\n",
      "{'0': 393, '1': 607}\n",
      "{'0': 413, '1': 587}\n",
      "diff=-20\n",
      "Shots: 1100:\n",
      "{'0': 432, '1': 568}\n",
      "Shots: 1200:\n",
      "{'0': 404, '1': 596}\n",
      "Shots: 1300:\n",
      "{'0': 429, '1': 571}\n",
      "Shots: 1400:\n",
      "{'0': 422, '1': 578}\n",
      "Shots: 1500:\n",
      "{'0': 425, '1': 575}\n",
      "Input:14\n",
      "{'0': 486, '1': 514}\n",
      "{'0': 430, '1': 570}\n",
      "diff=56\n",
      "Input:15\n",
      "{'0': 404, '1': 596}\n",
      "{'0': 406, '1': 594}\n",
      "diff=-2\n",
      "Shots: 1100:\n",
      "{'0': 403, '1': 597}\n",
      "Shots: 1200:\n",
      "{'0': 422, '1': 578}\n",
      "Shots: 1300:\n",
      "{'0': 420, '1': 580}\n",
      "Shots: 1400:\n",
      "{'0': 394, '1': 606}\n",
      "Shots: 1500:\n",
      "{'0': 390, '1': 610}\n",
      "Infected circuit 6\n",
      "Input:0\n",
      "{'0': 653, '1': 347}\n",
      "{'0': 670, '1': 330}\n",
      "diff=17\n",
      "Input:1\n",
      "{'0': 599, '1': 401}\n",
      "{'0': 600, '1': 400}\n",
      "diff=1\n",
      "Input:2\n",
      "{'0': 603, '1': 397}\n",
      "{'0': 635, '1': 365}\n",
      "diff=32\n",
      "Input:3\n",
      "{'0': 593, '1': 407}\n",
      "{'0': 599, '1': 401}\n",
      "diff=6\n",
      "Input:4\n",
      "{'0': 643, '1': 357}\n",
      "{'0': 642, '1': 358}\n",
      "diff=-1\n",
      "Shots: 1100:\n",
      "{'0': 645, '1': 355}\n",
      "Shots: 1200:\n",
      "{'0': 646, '1': 354}\n",
      "Shots: 1300:\n",
      "{'0': 654, '1': 346}\n",
      "Shots: 1400:\n",
      "{'0': 661, '1': 339}\n",
      "Shots: 1500:\n",
      "{'0': 624, '1': 376}\n",
      "Input:5\n",
      "{'0': 411, '1': 589}\n",
      "{'0': 381, '1': 619}\n",
      "diff=30\n",
      "Input:6\n",
      "{'0': 437, '1': 563}\n",
      "{'0': 438, '1': 562}\n",
      "diff=-1\n",
      "Shots: 1100:\n",
      "{'0': 461, '1': 539}\n",
      "Shots: 1200:\n",
      "{'0': 454, '1': 546}\n",
      "Shots: 1300:\n",
      "{'0': 440, '1': 560}\n",
      "Shots: 1400:\n",
      "{'0': 461, '1': 539}\n",
      "Shots: 1500:\n",
      "{'0': 460, '1': 540}\n",
      "Input:7\n",
      "{'0': 416, '1': 584}\n",
      "{'0': 409, '1': 591}\n",
      "diff=7\n",
      "Input:8\n",
      "{'0': 409, '1': 591}\n",
      "{'0': 418, '1': 582}\n",
      "diff=-9\n",
      "Shots: 1100:\n",
      "{'0': 441, '1': 559}\n",
      "Shots: 1200:\n",
      "{'0': 447, '1': 553}\n",
      "Shots: 1300:\n",
      "{'0': 458, '1': 542}\n",
      "Shots: 1400:\n",
      "{'0': 425, '1': 575}\n",
      "Shots: 1500:\n",
      "{'0': 438, '1': 562}\n",
      "Input:9\n",
      "{'0': 425, '1': 575}\n",
      "{'0': 411, '1': 589}\n",
      "diff=14\n",
      "Input:10\n",
      "{'0': 443, '1': 557}\n",
      "{'0': 426, '1': 574}\n",
      "diff=17\n",
      "Input:11\n",
      "{'0': 400, '1': 600}\n",
      "{'0': 368, '1': 632}\n",
      "diff=32\n",
      "Input:12\n",
      "{'0': 422, '1': 578}\n",
      "{'0': 594, '1': 406}\n",
      "Input:13\n",
      "{'0': 408, '1': 592}\n",
      "{'0': 416, '1': 584}\n",
      "diff=-8\n",
      "Shots: 1100:\n",
      "{'0': 381, '1': 619}\n",
      "Shots: 1200:\n",
      "{'0': 420, '1': 580}\n",
      "Shots: 1300:\n",
      "{'0': 412, '1': 588}\n",
      "Shots: 1400:\n",
      "{'0': 442, '1': 558}\n",
      "Shots: 1500:\n",
      "{'0': 437, '1': 563}\n",
      "Input:14\n",
      "{'0': 436, '1': 564}\n",
      "{'0': 599, '1': 401}\n",
      "Shots: 1100:\n",
      "{'0': 608, '1': 392}\n",
      "Shots: 1200:\n",
      "{'0': 609, '1': 391}\n",
      "Shots: 1300:\n",
      "{'0': 600, '1': 400}\n",
      "Shots: 1400:\n",
      "{'0': 618, '1': 382}\n",
      "Shots: 1500:\n",
      "{'0': 648, '1': 352}\n",
      "Input:15\n",
      "{'0': 402, '1': 598}\n",
      "{'0': 427, '1': 573}\n",
      "diff=-25\n",
      "Shots: 1100:\n",
      "{'0': 419, '1': 581}\n",
      "Shots: 1200:\n",
      "{'0': 427, '1': 573}\n",
      "Shots: 1300:\n",
      "{'0': 389, '1': 611}\n",
      "Shots: 1400:\n",
      "{'0': 395, '1': 605}\n",
      "Shots: 1500:\n",
      "{'0': 437, '1': 563}\n",
      "Infected circuit 7\n",
      "Input:0\n",
      "{'0': 642, '1': 358}\n",
      "{'0': 382, '1': 618}\n",
      "Shots: 1100:\n",
      "{'0': 368, '1': 632}\n",
      "Shots: 1200:\n",
      "{'0': 383, '1': 617}\n",
      "Shots: 1300:\n",
      "{'0': 386, '1': 614}\n",
      "Shots: 1400:\n",
      "{'0': 381, '1': 619}\n",
      "Shots: 1500:\n",
      "{'0': 407, '1': 593}\n",
      "Input:1\n",
      "{'0': 666, '1': 334}\n",
      "{'0': 434, '1': 566}\n",
      "Shots: 1100:\n",
      "{'0': 412, '1': 588}\n",
      "Shots: 1200:\n",
      "{'0': 454, '1': 546}\n",
      "Shots: 1300:\n",
      "{'0': 458, '1': 542}\n",
      "Shots: 1400:\n",
      "{'0': 416, '1': 584}\n",
      "Shots: 1500:\n",
      "{'0': 430, '1': 570}\n",
      "Input:2\n",
      "{'0': 631, '1': 369}\n",
      "{'0': 412, '1': 588}\n",
      "Shots: 1100:\n",
      "{'0': 396, '1': 604}\n",
      "Shots: 1200:\n",
      "{'0': 398, '1': 602}\n",
      "Shots: 1300:\n",
      "{'0': 402, '1': 598}\n",
      "Shots: 1400:\n",
      "{'0': 401, '1': 599}\n",
      "Shots: 1500:\n",
      "{'0': 422, '1': 578}\n",
      "Input:3\n",
      "{'0': 590, '1': 410}\n",
      "{'0': 464, '1': 536}\n",
      "Shots: 1100:\n",
      "{'0': 425, '1': 575}\n",
      "Shots: 1200:\n",
      "{'0': 447, '1': 553}\n",
      "Shots: 1300:\n",
      "{'0': 466, '1': 534}\n",
      "Shots: 1400:\n",
      "{'0': 453, '1': 547}\n",
      "Shots: 1500:\n",
      "{'0': 437, '1': 563}\n",
      "Input:4\n",
      "{'0': 644, '1': 356}\n",
      "{'0': 378, '1': 622}\n",
      "Shots: 1100:\n",
      "{'0': 374, '1': 626}\n",
      "Shots: 1200:\n",
      "{'0': 407, '1': 593}\n",
      "Shots: 1300:\n",
      "{'0': 387, '1': 613}\n",
      "Shots: 1400:\n",
      "{'0': 422, '1': 578}\n",
      "Shots: 1500:\n",
      "{'0': 409, '1': 591}\n",
      "Input:5\n",
      "{'0': 426, '1': 574}\n",
      "{'0': 629, '1': 371}\n",
      "Shots: 1100:\n",
      "{'0': 604, '1': 396}\n",
      "Shots: 1200:\n",
      "{'0': 639, '1': 361}\n",
      "Shots: 1300:\n",
      "{'0': 637, '1': 363}\n",
      "Shots: 1400:\n",
      "{'0': 636, '1': 364}\n",
      "Shots: 1500:\n",
      "{'0': 619, '1': 381}\n",
      "Input:6\n",
      "{'0': 458, '1': 542}\n",
      "{'0': 615, '1': 385}\n",
      "Shots: 1100:\n",
      "{'0': 609, '1': 391}\n",
      "Shots: 1200:\n",
      "{'0': 582, '1': 418}\n",
      "Shots: 1300:\n",
      "{'0': 625, '1': 375}\n",
      "Shots: 1400:\n",
      "{'0': 593, '1': 407}\n",
      "Shots: 1500:\n",
      "{'0': 606, '1': 394}\n",
      "Input:7\n",
      "{'0': 440, '1': 560}\n",
      "{'0': 622, '1': 378}\n",
      "Shots: 1100:\n",
      "{'0': 610, '1': 390}\n",
      "Shots: 1200:\n",
      "{'0': 627, '1': 373}\n",
      "Shots: 1300:\n",
      "{'0': 606, '1': 394}\n",
      "Shots: 1400:\n",
      "{'0': 617, '1': 383}\n",
      "Shots: 1500:\n",
      "{'0': 630, '1': 370}\n",
      "Input:8\n",
      "{'0': 453, '1': 547}\n",
      "{'0': 619, '1': 381}\n",
      "Shots: 1100:\n",
      "{'0': 620, '1': 380}\n",
      "Shots: 1200:\n",
      "{'0': 622, '1': 378}\n",
      "Shots: 1300:\n",
      "{'0': 601, '1': 399}\n",
      "Shots: 1400:\n",
      "{'0': 634, '1': 366}\n",
      "Shots: 1500:\n",
      "{'0': 604, '1': 396}\n",
      "Input:9\n",
      "{'0': 414, '1': 586}\n",
      "{'0': 679, '1': 321}\n",
      "Shots: 1100:\n",
      "{'0': 637, '1': 363}\n",
      "Shots: 1200:\n",
      "{'0': 631, '1': 369}\n",
      "Shots: 1300:\n",
      "{'0': 632, '1': 368}\n",
      "Shots: 1400:\n",
      "{'0': 647, '1': 353}\n",
      "Shots: 1500:\n",
      "{'0': 634, '1': 366}\n",
      "Input:10\n",
      "{'0': 407, '1': 593}\n",
      "{'0': 613, '1': 387}\n",
      "Shots: 1100:\n",
      "{'0': 595, '1': 405}\n",
      "Shots: 1200:\n",
      "{'0': 617, '1': 383}\n",
      "Shots: 1300:\n",
      "{'0': 615, '1': 385}\n",
      "Shots: 1400:\n",
      "{'0': 603, '1': 397}\n",
      "Shots: 1500:\n",
      "{'0': 620, '1': 380}\n",
      "Input:11\n",
      "{'0': 386, '1': 614}\n",
      "{'0': 633, '1': 367}\n",
      "Shots: 1100:\n",
      "{'0': 632, '1': 368}\n",
      "Shots: 1200:\n",
      "{'0': 641, '1': 359}\n",
      "Shots: 1300:\n",
      "{'0': 631, '1': 369}\n",
      "Shots: 1400:\n",
      "{'0': 637, '1': 363}\n",
      "Shots: 1500:\n",
      "{'0': 626, '1': 374}\n",
      "Input:12\n",
      "{'0': 448, '1': 552}\n",
      "{'0': 604, '1': 396}\n",
      "Shots: 1100:\n",
      "{'0': 618, '1': 382}\n",
      "Shots: 1200:\n",
      "{'0': 604, '1': 396}\n",
      "Shots: 1300:\n",
      "{'0': 617, '1': 383}\n",
      "Shots: 1400:\n",
      "{'0': 595, '1': 405}\n",
      "Shots: 1500:\n",
      "{'0': 594, '1': 406}\n",
      "Input:13\n",
      "{'0': 410, '1': 590}\n",
      "{'0': 630, '1': 370}\n",
      "Shots: 1100:\n",
      "{'0': 646, '1': 354}\n",
      "Shots: 1200:\n",
      "{'0': 653, '1': 347}\n",
      "Shots: 1300:\n",
      "{'0': 632, '1': 368}\n",
      "Shots: 1400:\n",
      "{'0': 667, '1': 333}\n",
      "Shots: 1500:\n",
      "{'0': 666, '1': 334}\n",
      "Input:14\n",
      "{'0': 457, '1': 543}\n",
      "{'0': 591, '1': 409}\n",
      "Shots: 1100:\n",
      "{'0': 607, '1': 393}\n",
      "Shots: 1200:\n",
      "{'0': 587, '1': 413}\n",
      "Shots: 1300:\n",
      "{'0': 571, '1': 429}\n",
      "Shots: 1400:\n",
      "{'0': 572, '1': 428}\n",
      "Shots: 1500:\n",
      "{'0': 619, '1': 381}\n",
      "Input:15\n",
      "{'0': 397, '1': 603}\n",
      "{'0': 623, '1': 377}\n",
      "Shots: 1100:\n",
      "{'0': 618, '1': 382}\n",
      "Shots: 1200:\n",
      "{'0': 610, '1': 390}\n",
      "Shots: 1300:\n",
      "{'0': 603, '1': 397}\n",
      "Shots: 1400:\n",
      "{'0': 616, '1': 384}\n",
      "Shots: 1500:\n",
      "{'0': 638, '1': 362}\n"
     ]
    }
   ],
   "source": [
    "for i, new_mod_qc in enumerate(new_circuit_list):\n",
    "    print(f\"Infected circuit {i}\")\n",
    "    ## cal_qc + out_qc = cal_out_qc\n",
    "    new_mod_qc.draw(output='mpl')\n",
    "    cal_out_qc = new_mod_qc.compose(out_qc, front=False, inplace=False)\n",
    "    cal_out_qc.draw(output='mpl')\n",
    "    ## in_qc + cal_out_qc= org_qc\n",
    "    for i, in_qc in enumerate(in_circuits):\n",
    "        print(f\"Input:{i}\")\n",
    "        org_qc = cal_out_qc.compose(in_qc, front=True, inplace=False)\n",
    "        org_qc.draw(output='mpl')\n",
    "        # Run the circuit and get results\n",
    "        #shots = [i for i in range(1000, 1501, 100)]\n",
    "        #for i in shots:\n",
    "        #print(f\"Shots: {i}:\")   \n",
    "\n",
    "        cal_clean_qc = cal_qc.compose(out_qc, front=False, inplace=False)\n",
    "        #cal_clean_qc.draw(output='mpl')\n",
    "        org_c_qc = cal_clean_qc.compose(in_qc, front=True, inplace=False)\n",
    "        #org_c_qc.draw(output='mpl')\n",
    "        # Run the circuit and get results\n",
    "        #shots = 1000    \n",
    "        backend = FakeValencia()\n",
    "        qc = transpile(org_c_qc, backend=backend, initial_layout= [0,1,2,3,4], optimization_level=1)\n",
    "        results = backend.run(qc,shots=1000).result()\n",
    "        \n",
    "        # Extract the counts from the results\n",
    "        count = results.get_counts(qc)\n",
    "\n",
    "        #sort the counts\n",
    "\n",
    "        count = dict(sorted(count.items()))\n",
    "\n",
    "        print(count) \n",
    "        \n",
    "        backend = FakeValencia()\n",
    "        qc_noisy = transpile(org_qc, backend=backend, initial_layout= [0,1,2,3,4], optimization_level=1)\n",
    "        results = backend.run(qc_noisy,shots=1000).result()\n",
    "        \n",
    "        # Extract the counts from the results\n",
    "        count_noisy = results.get_counts(qc_noisy)\n",
    "\n",
    "        #sort the counts\n",
    "\n",
    "        count_noisy = dict(sorted(count_noisy.items()))\n",
    "\n",
    "        print(count_noisy)\n",
    "        #if i == 1000:\n",
    "        \n",
    "\n",
    "        if ((count_noisy['1']/1000 > 0.5) and (count['1']/1000 > 0.5)):\n",
    "            diff = count_noisy['1']-count['1']\n",
    "            print(f\"diff={diff}\")\n",
    "            #draw a histogram\n",
    "        \n",
    "        if ((count_noisy['0']/1000 > 0.5) and (count['0']/1000 > 0.5)):\n",
    "            diff = count_noisy['0']-count['0']\n",
    "            print(f\"diff={diff}\")  \n",
    "\n",
    "        if diff < 0:\n",
    "            shots = [i for i in range(1100, 1501, 100)]\n",
    "            for i in shots:\n",
    "                print(f\"Shots: {i}:\")  \n",
    "                qc_noisy = transpile(org_qc, backend=backend, initial_layout= [0,1,2,3,4], optimization_level=1)\n",
    "                results = backend.run(qc_noisy,shots=1000).result()\n",
    "                \n",
    "                # Extract the counts from the results\n",
    "                count_noisy = results.get_counts(qc_noisy)\n",
    "\n",
    "                #sort the counts\n",
    "\n",
    "                count_noisy = dict(sorted(count_noisy.items()))\n",
    "\n",
    "                print(count_noisy)\n",
    "        "
   ]
  }
 ],
 "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.10.10"
  },
  "orig_nbformat": 4
 },
 "nbformat": 4,
 "nbformat_minor": 2
}