3_17_v2.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 6,
   "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": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "from qiskit import QuantumRegister, ClassicalRegister, QuantumCircuit\n",
    "from numpy import pi\n",
    "\n",
    "qreg_q = QuantumRegister(3, 'q')\n",
    "creg_c = ClassicalRegister(4, 'c')\n",
    "cal_qc = QuantumCircuit(qreg_q, creg_c)\n",
    "\n",
    "cal_qc.x(qreg_q[2])\n",
    "cal_qc.cx(qreg_q[0], qreg_q[2])\n",
    "cal_qc.cx(qreg_q[2], qreg_q[1])\n",
    "cal_qc.ccx(qreg_q[2], qreg_q[1], qreg_q[0])\n",
    "cal_qc.ccx(qreg_q[0], qreg_q[1], qreg_q[2])\n",
    "cal_qc.cx(qreg_q[1], qreg_q[2])\n",
    "# @columns [0,1,2,3,4,5]\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 = '3_17_v2.qasm'\n",
    "with open(file_path, 'w') as file:\n",
    "    file.write(compiled_cal_qc.qasm())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "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",
      "c: 1/═══\n",
      "        \n",
      "Circuit 2:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "           ░ \n",
      "q_1: ──────░─\n",
      "           ░ \n",
      "q_2: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_3: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 3:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "           ░ \n",
      "q_1: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_2: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_3: ──────░─\n",
      "           ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 4:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "           ░ \n",
      "q_1: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_2: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_3: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 5:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_1: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_2: ──────░─\n",
      "           ░ \n",
      "q_3: ──────░─\n",
      "           ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 6:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_1: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_2: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_3: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 7:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_1: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_2: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_3: ──────░─\n",
      "           ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 8:\n",
      "           ░ \n",
      "q_0: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_1: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_2: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_3: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 9:\n",
      "     ┌───┐ ░ \n",
      "q_0: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_1: ──────░─\n",
      "           ░ \n",
      "q_2: ──────░─\n",
      "           ░ \n",
      "q_3: ──────░─\n",
      "           ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 10:\n",
      "     ┌───┐ ░ \n",
      "q_0: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_1: ──────░─\n",
      "           ░ \n",
      "q_2: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_3: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 11:\n",
      "     ┌───┐ ░ \n",
      "q_0: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_1: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_2: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_3: ──────░─\n",
      "           ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 12:\n",
      "     ┌───┐ ░ \n",
      "q_0: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_1: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_2: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_3: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 13:\n",
      "     ┌───┐ ░ \n",
      "q_0: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_1: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_2: ──────░─\n",
      "           ░ \n",
      "q_3: ──────░─\n",
      "           ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 14:\n",
      "     ┌───┐ ░ \n",
      "q_0: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_1: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_2: ──────░─\n",
      "     ┌───┐ ░ \n",
      "q_3: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 15:\n",
      "     ┌───┐ ░ \n",
      "q_0: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_1: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_2: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "q_3: ──────░─\n",
      "           ░ \n",
      "c: 1/════════\n",
      "             \n",
      "Circuit 16:\n",
      "     ┌───┐ ░ \n",
      "q_0: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_1: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_2: ┤ X ├─░─\n",
      "     ├───┤ ░ \n",
      "q_3: ┤ X ├─░─\n",
      "     └───┘ ░ \n",
      "c: 1/════════\n",
      "             \n"
     ]
    }
   ],
   "source": [
    "from qiskit import QuantumRegister, ClassicalRegister, QuantumCircuit\n",
    "\n",
    "n_qbits = 4\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])\n",
    "\n",
    "    in_qc.barrier(qreg_q[0], qreg_q[1], qreg_q[2], qreg_q[3])\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": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPEAAAFvCAYAAACb0C2OAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy88F64QAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAhqklEQVR4nO3dfVQU9+Eu8Gc3oCsroIhlCSDvVECBKJAS2xi8JgUMmnsTExJj6jm2/lL1ao0/MTHmJH9woqSe5NaYtqZobNrUS4ppYgBtmkga5DYGglAPL1EMLy7uqlteZEERdvf+YdwjAZRdhx2+O8/nnB7CzszyWH34zsx3ZlZls9lsICJhqeUOQER3hiUmEhxLTCQ4lphIcCwxkeBYYiLBscREgmOJiQTHEhMJjiUmEhxLTCQ4lphIcCwxkeBYYiLBscREgmOJiQTHEhMJjiUmEhxLTCQ4lphIcCwxkeBYYiLBscREgmOJiQTHEhMJjiUmEhxLTCQ4lphIcB5yB6A709jY6ND6Fy9exPvvv4/HH38cP/jBD8a0zezZs52JRi7CkVhhLl26hLfeeguXLl2SOwpJhCUmEhxLTCQ4lphIcCyxwvj6+iI7Oxu+vr5yRyGJqGw2m03uEOQ8R89OO4Nnpyc2jsQK09/fj9bWVvT398sdhSTCEitMU1MTMjIy0NTUJHcUkghLTCQ4lphIcCwxkeBYYiLBcYpJcJxiIo7ERIJjiRWmubkZOTk5aG5uljsKSYQlVpi+vj7U1tair69P7igkEZaYSHAsMZHgWGIiwbHEChMUFIT8/HwEBQXJHYUkwnliwXGemDgSK0xHRwfee+89dHR0yB2FJMISK4zBYEBeXh4MBoPcUUgiiiixyWRCbm4uoqKioNFoEBISgo0bN6K3txerV6+GSqXCnj175I5J5BS3f3h8TU0NMjMzYTQaodVqERcXh/Pnz2P37t04e/asfbcyKSlJ3qBETnLrkdhkMiE7OxtGoxGbN2+GwWBAdXU1jEYj8vPzUVJSgsrKSqhUKiQkJMgdl8gpbl3iDRs2QK/XY/369di1axe8vb3ty3Jzc5GYmIjBwUGEhYXBx8dHxqSuo9VqsWDBAmi1WrmjkETctsQNDQ0oLCyEv78/duzYMeI68+fPBwAkJiYOeb25uRlLly6Ft7c3pk+fjmeeeQb/+c9/xj2zK4SFhaGgoABhYWFyRyGJuG2JDx48CKvVihUrVmDq1KkjrjNlyhQAQ0vc09OD9PR06PV6HDx4EG+//TbKy8vx8MMPw2q1uiT7eLJYLDCbzbBYLHJHIYm47YmtY8eOAQDS09NHXUev1wMYWuK3334b7e3t+OKLLzBr1iwAQHBwMO677z4cPnwYjzzyyPiFdoHGxkY89thjKCoqQnx8vNxxSAJuW+LW1lYAQGho6IjLBwcHUVFRAWBoiYuLi/HjH//YXmAASEtLQ0REBD7++GOnSpycnAyj0ejwdmOxatUqh9a/MT9cWlqKkydPjmmbAwcOOJiKHKXT6VBVVeXUtm5b4t7eXgDAlStXRlxeWFgIk8kEb29vhIeH21+vr6/H8uXLh60fHx+P+vp6p7IYjUa0t7c7te3tOHpf8NWrV+1fx7rteGUnabhtiXU6HTo7O1FdXY20tLQhywwGA7Zs2QIASEhIgEqlsi/r7OzEtGnThr2fn58fvvnmG6ezjBcvLy+H1tdoNPavY92WN0uMvzv5N+K2JV68eDEaGhqQn5+PBx98EDExMQCAyspKrFy5EiaTCYBrLvJwdjdpLBy9AaKurg779+9HVlbWmI+J8/LynIlGLuK2Z6dzc3MxY8YMnDt3DvHx8Zg7dy6io6ORmpqKiIgILFq0CMDw6aXp06ejq6tr2Pt1dHTAz8/PFdHHVUxMDCoqKuy/1Eh8blvi4OBglJeXY8mSJdBoNGhpaYGfnx/27t2LkpISnD59GsDwEsfGxo547FtfX4/Y2FiXZB9Pnp6e8PPzg6enp9xRSCJuW2LgeiGLi4vR09ODnp4enDhxAmvWrEFvby9aWlqgVqsxZ86cIds8/PDDOH78uH36CQBOnDiBs2fPIjs729V/BMm1tbVh7dq1aGtrkzsKScStSzyauro62Gw2REdHDzu5s2bNGgQGBmLZsmUoLi5GUVERnnzySaSmpmLZsmUyJZZOT08PysrK0NPTI3cUkogiS3zq1CkAw3elAcDHxwfHjh1DYGAgcnJy8POf/xz33XcfiouLoVYr8v8umuDc9uz0rdyqxAAQGRmJ4uJiV0Yicpoih5bblZhIJIociW9cV61EAQEB2Lp1KwICAuSOQhLh0y4Fx6ddkiJ3p5Wsu7sbR48eRXd3t9xRSCIsscLo9Xps2rRpyDw4iY0lJhIcS0wkOJaYSHAsscJoNBrExsba7ysm8XGKSXCcYiKOxESCY4kVpr6+HgkJCU4/L4wmHpZYYWw2GwYGBsCjKPfBEhMJjiUmEhxLTCQ4Rd6KqGSRkZE4fPgwQkJC5I5CEmGJFUaj0SA6OlruGCQh7k4rTHt7O7Zv386PZnEjLLHCdHV14dChQyM+IJ/ExBITCY4lJhIcS0wkOJZYYdRqNVJSUvggfDfCv0mFsVqtqKyshNVqlTsKSYQlJhIcS0wkOJaYSHAsscL4+voiOzsbvr6+ckchifAZW4LjM7aII7HC9Pf3o7W1Ff39/XJHIYmwxArT1NSEjIwMNDU1yR2FJMISEwmOJSYSHEtMJDiWmEhwnGISHKeYiCMxkeBYYoVpbm5GTk4Ompub5Y5CEmGJFaavrw+1tbXo6+uTOwpJhCUmEhxLTCQ4lphIcCyxwgQFBSE/Px9BQUFyRyGJcJ5YcJwnJo7ECtPR0YH33nsPHR0dckchibDECmMwGJCXlweDwSB3FJKIIkpsMpmQm5uLqKgoaDQahISEYOPGjejt7cXq1auhUqmwZ88euWMSOcXtP9q0pqYGmZmZMBqN0Gq1iIuLw/nz57F7926cPXvWvluZlJQkb1AiJ7n1SGwymZCdnQ2j0YjNmzfDYDCguroaRqMR+fn5KCkpQWVlJVQqFRISEuSOS+QUty7xhg0boNfrsX79euzatQve3t72Zbm5uUhMTMTg4CDCwsLg4+MjY1LX0Wq1WLBgAbRardxRSCJuW+KGhgYUFhbC398fO3bsGHGd+fPnAwASExPtr90ofWpqKiZPngyVSuWSvK4SFhaGgoIChIWFyR2FJOK2JT548CCsVitWrFiBqVOnjrjOlClTAAwtcVNTEw4dOgSdToeUlBSXZHUli8UCs9kMi8UidxSSiNuW+NixYwCA9PT0UdfR6/UAhpb4/vvvh8FgwOHDh7F48eLxDSmDxsZGpKSkuOQiEXINtz073draCgAIDQ0dcfng4CAqKioADC3xeHzkZ3JyMoxGo+TvCwCrVq1yaP0b88OlpaU4efLkmLY5cOCAg6nIUTqdDlVVVU5t67Yl7u3tBQBcuXJlxOWFhYUwmUzw9vZGeHj4uGYxGo1ob28fl/d29L7gq1ev2r+Oddvxyk7ScNsS63Q6dHZ2orq6GmlpaUOWGQwGbNmyBQCQkJAw7ievdDrduL23l5eXQ+trNBr717Fuy5slxt+d/Btx2xIvXrwYDQ0NyM/Px4MPPoiYmBgAQGVlJVauXAmTyQTANRd5OLubNBaOHtvW1dVh//79yMrKQnx8/Ji2ycvLcyYauYjbntjKzc3FjBkzcO7cOcTHx2Pu3LmIjo5GamoqIiIisGjRIgBDj4eVICYmBhUVFfZfaiQ+ty1xcHAwysvLsWTJEmg0GrS0tMDPzw979+5FSUkJTp8+DUB5Jfb09ISfnx88PT3ljkIScdvdaQCIjY1FcXHxsNfNZjNaWlqgVqsxZ84cGZLJp62tDTt37sTzzz+PWbNmyR2HJODWJR5NXV0dbDYbYmJiRjy5U1RUBACor68f8n1YWBiSk5NdF3Qc9PT0oKysDOvWrZM7CklEkSU+deoUgNF3pZcvXz7i9z/72c84Z0oTDks8Aj6xiETitie2buV2JSYSiSJH4hvXVStRQEAAtm7dioCAALmjkET4tEvB8WmXpMjdaSXr7u7G0aNH0d3dLXcUkghLrDB6vR6bNm2y34ZJ4mOJiQTHEhMJjiUmEhxLrDAajQaxsbH2+4pJfJxiEhynmIgjMZHgWGKFqa+vR0JCgv0OLRIfS6wwNpsNAwMDvMnDjbDERIJjiYkExxITCU6RtyIqWWRkJA4fPoyQkBC5o5BEWGKF0Wg0iI6OljsGSYi70wrT3t6O7du386NZ3AhLrDBdXV04dOgQurq65I5CEmGJiQTHEhMJjiUmEhxLrDBqtRopKSnj8mHqJA/+TSqM1WpFZWUlrFar3FFIIiwxkeBYYiLBscREgmOJFcbX1xfZ2dnw9fWVOwpJhM/YEhyfsUUciRWmv78fra2t6O/vlzsKSYQlVpimpiZkZGSgqalJ7igkEZaYSHAsMZHgWGIiwbHERILjFJPgOMVEHImJBMcSK0xzczNycnLQ3NwsdxSSCEusMH19faitrUVfX5/cUUgiLDGR4FhiIsGxxESCY4kVJigoCPn5+QgKCpI7CkmE88SC4zwxcSRWmI6ODrz33nvo6OiQOwpJhCVWGIPBgLy8PBgMBrmjkERYYiLBKaLEJpMJubm5iIqKgkajQUhICDZu3Ije3l6sXr0aKpUKe/bskTsmkVPc/vOJa2pqkJmZCaPRCK1Wi7i4OJw/fx67d+/G2bNn7ceGSUlJ8gYlcpJbj8QmkwnZ2dkwGo3YvHkzDAYDqqurYTQakZ+fj5KSElRWVkKlUiEhIUHuuC6h1WqxYMECaLVauaOQRNx6iumpp57CwYMHsX79erz55pvDliclJaG2thbh4eH49ttvZUh45zjFRG47Ejc0NKCwsBD+/v7YsWPHiOvMnz8fAJCYmGh/raioCI8++ihCQ0Ph5eWF2bNn48UXX4TZbHZJ7vFmsVhgNpthsVjkjkIScdsSHzx4EFarFStWrMDUqVNHXGfKlCkAhpZ4165duOuuu/Dqq6/iyJEj+OUvf4nf/e53yMjIcIsPIWtsbERKSopLRnByDbc9sXXs2DEAQHp6+qjr6PV6AENL/PHHH2PmzJn27xcuXIiZM2dixYoVOH78OO6///5xSkzkHLctcWtrKwAgNDR0xOWDg4OoqKgAMLTENxf4huTkZABAe3u7U1mSk5NhNBqd2vZ2Vq1a5dD6Ny7yKC0txcmTJ8e0zYEDBxxMRY7S6XSoqqpyalu3LXFvby8A4MqVKyMuLywshMlkgre3N8LDw2/5XmVlZQCA2NhYp7IYjUanfwHcjqM391+9etX+dazbjld2kobbllin06GzsxPV1dVIS0sbssxgMGDLli0AgISEBKhUqlHfp729HS+99BIyMjKcnkvW6XRObTcWXl5eDq2v0WjsX8e6Le94Gn938m/EbaeYNmzYgDfffBMhISH49NNPERMTAwCorKzEypUr8e2332JgYADr1q0b9Wots9mMBx54AEajEZWVlQgMDHTlH2FMHD1BNTAwgJ6eHnh7e8PT03NM23CKaWJz27PTubm5mDFjBs6dO4f4+HjMnTsX0dHRSE1NRUREBBYtWgRg6PHwza5cuYLs7Gw0Nzfjk08+mZAFdoanpyf8/PzGXGCa+Ny2xMHBwSgvL8eSJUug0WjQ0tICPz8/7N27FyUlJTh9+jSAkUs8MDCAxx57DFVVVThy5Aji4uJcHX/ctLW1Ye3atWhra5M7CknEbY+JgesnooqLi4e9bjab0dLSArVajTlz5gxZdmNu+bPPPkNpaSlSU1NdFdclenp6UFZWhnXr1skdhSTi1iUeTV1dHWw2G2JiYoad3Fm3bh3++te/4vnnn4eXlxe+/PJL+7LIyMgRp6CI5OS2u9O3curUKQAj70ofOXIEALBz506kpaUN+V9JSYlLcxKNhSJH4luVuKWlxcVpiO4MR2KFCQgIwNatWxEQECB3FJKI284TKwVvRSRFjsRK1t3djaNHj6K7u1vuKCQRllhh9Ho9Nm3aZL+Di8THEhMJjiUmEhxLTCQ4llhhNBoNYmNj7bckkvg4xSQ4TjERR2IiwbHEClNfX4+EhATU19fLHYUkwhIrjM1mw8DAAHgU5T5YYiLBscREgmOJiQSnyPuJlSwyMhKHDx9GSEiI3FFIIiyxwmg0GkRHR8sdgyTE3WmFaW9vx/bt2/mpDm6EJVaYrq4uHDp0CF1dXXJHIYmwxESCY4mJBMcSEwmOJVYYf39//OIXv4C/v7/cUUgiLLHCqFQqTJo06ZYf50piYYkV5tKlS3jrrbdw6dIluaOQRFhiIsGxxESCY4mJBMcSK4yvry+ys7Ph6+srdxSSCB+UJzg+KI84EitMf38/Wltb0d/fL3cUkghLrDBNTU3IyMhAU1OT3FFIIryfmCYkm82Gvr4+uWM4xMvLS5aLaFhimpD6+vowdepUuWM4xGw2Q6vVuvzncneaSHAsMZHguDutMPHx8WhoaJA7BkmIIzGR4FhihWlubkZOTg6am5vljkISYYkVpq+vD7W1tcJN39DoWGIiwbHERIJjiYkExxIrTFBQEPLz8xEUFCR3FJIIS6ww06ZNw9KlSzFt2jS5owhn1qxZckcYES/2UJiOjg4cOXIEmZmZ8PPzkzvOuJs9ezaSk5Mxf/58zJ49G15eXrBarbh8+TL+/e9/4+uvv8aJEydgMBhu+T7Lly/Hn//8Zzz77LN45513XJR+bFhihTEYDMjLy0NSUpLbllir1eKpp57C2rVrkZSUNOp6S5cuBQBYrVYcPXoUv/3tb3HkyBFYrdYh6y1fvhx/+ctf4OHhgYKCApw5cwbHjx8fzz+CQxSxO20ymZCbm4uoqChoNBqEhIRg48aN6O3txerVq6FSqbBnzx65Y5IEnn76abS1teHtt9++ZYFvplarkZWVheLiYpw6dQr33nuvfdnNBQaA/fv3o6KiYjyiO83tR+KamhpkZmbCaDRCq9UiLi4O58+fx+7du3H27Fl0dHQAwJj/wmliCggIwB/+8AdkZ2cPef3LL7/ERx99hK+//hq1tbXo6uqCSqWCTqfDvHnzkJKSgieffBJhYWEAgLi4OFRUVGDXrl2ora3Fu+++ay9wQUEB1qxZg4n2RCu3fsaWyWTCPffcA71ej82bN+Pll1+Gt7c3AOC1117D1q1b4eHhAYvFgq6uLvj4+Mic2HGOPmOrrq4Ojz32GIqKihAfHz+mbeR4xlZvb++Y7ycODw/Hp59+ioiICPtr77//PvLz81FdXX3b7dVqNTIyMvDKK68gJSXF/rrVaoVafX1ndSwF5v3E42DDhg3Q6/VYv349du3aZS8wAOTm5iIxMRGDg4MICwsTssDO0Gq1WLBggSz/2MZDSEgIPv/8c3uBjUYjHnnkETzxxBNjKjBwvaylpaVIS0vDCy+8gIGBAQCwF3j//v0TcgS+wW1L3NDQgMLCQvj7+2PHjh0jrjN//nwAQGJiov218vJyLF68GIGBgZg8eTKCg4PxxBNPuM3te2FhYSgoKLDvPorMw8MDH374oX3qp66uDvPmzcNHH33k1PtZLBacPXt22CN2rly5MmELDLhxiQ8ePAir1YoVK1aMuls2ZcoUAENL3NnZiblz52L37t345JNPkJ+fj7q6OqSlpUGv17sk+3iyWCwwm82wWCxyR7ljL7zwAubNmwcAOHPmDBYtWnTbqaJb+f5JrMHBQQDAunXrkJ6efueBx4nblvjYsWMAcMv/82+U8uYSL126FG+88QaWL1+OhQsXYsWKFfjggw/Q3d2NQ4cOjW9oF2hsbERKSopLnlc9nuLj4/HSSy8BuF62nJwcXLx40en3+36BCwoK8Nxzz9mX79u3z/5Lf6Jx27PTra2tAIDQ0NARlw8ODtqnCm4u8UhmzJgBAPa/YEclJyfDaDQ6te3trFq1yqH1b4xUpaWlOHny5Ji2OXDggIOp7tz352q/77//+7/h6ekJANi5c+eYj39HMlKB16xZAwB49NFHsXDhQoSHhyMnJ+eWF3pER0fbj6MdpdPpUFVV5dS2blvi3t5eANePZ0ZSWFgIk8kEb29vhIeHD1tusVhgtVrR2tqKF154ATqdDo8//rhTWYxGI9rb253a9nYcvS/46tWr9q9j3Xa8sjvLz88POTk5AK4f/rz66qtOv9doBb5xDJybm4sTJ04AANauXXvLEt/JrvydcNsS63Q6dHZ2orq6GmlpaUOWGQwGbNmyBQCQkJAw4rOCFy5caB+po6KicOzYMcycOdPpLOPFy8vLofU1Go3961i3leNmCavVOmopnn76afuf45133hn1F/Xt3K7AAPDVV1+hqqoKycnJSE5ORlJSEmpqakZ8v8DAwDsaiZ3ltiVevHgxGhoakJ+fjwcffBAxMTEAgMrKSqxcuRImkwnA6Bd57Nu3D11dXWhubsavf/1rPPTQQ6ioqHDqInhnd5PGwpl54v379yMrK2vM88R5eXnORLsjt5on/slPfmL/b2d39cdS4Jt/RnJysv1nj1biM2fOcJ5YSrm5uZgxYwbOnTuH+Ph4zJ07F9HR0UhNTUVERAQWLVoEYPTj4R/+8Ie49957kZOTg88++ww9PT147bXXXPlHGBcxMTGoqKiw/1IT0Y2pQbPZjLq6Ooe3d6TAwPWrvr7/sycSty1xcHAwysvLsWTJEmg0GrS0tMDPzw979+5FSUkJTp8+DeD2J7WA67fvRUVFucXnF3l6esLPz89+Ukg006ZNs5/DOHny5G1PgH2fowUGgFOnTuHatWsAYJ/SmkjctsQAEBsbi+LiYvT09KCnpwcnTpzAmjVr0Nvbi5aWFqjVasyZM+e273Px4kV88803iIyMdEHq8dXW1oa1a9eira1N7ihOmT59uv2/z50759C2zhQYAK5du2afvrr5508UbntMfCt1dXWw2WyIiYkZdnLn6aefRlRUFJKSkjBt2jScOXMGb7zxBjw8PLBp0yaZEkunp6cHZWVlWLdundxRnGIwGHD//fdDo9Hg0qVLY95OpVLhV7/6ldM3MzzzzDNQq9Uwm81O5R5PiizxqVOnAIy8K/2jH/0I7777Ln7zm9/g6tWrCAkJQXp6OrZt2zbqnDO5ztWrV1FeXu7wdjabDVlZWfjHP/6B2tpah6+FLisrc/hnugpL/D3r16/H+vXrXR2JXKC7uxuLFi1Cb2/vhL4W2lEsMSnKRNwdvlOKLPGN66qVKCAgAFu3bkVAQIDcUUgiiiyxkvn7+zt8vTVNbG49xUTDdXd34+jRo+ju7pY7CkmEJVYYvV6PTZs2ucW90XQdS0wkOJaYSHAsMZHgWGKF0Wg0iI2Ntd+PS+LjFJPCREZG4oMPPpA7BkmIIzGR4Fhihamvr0dCQgLq6+vljkISYYkVxmazYWBgwK1uAFA6HhPThOTl5SXpzQq/3vt/cbm3Dz5aL2z5r5xh30vB0YcWSoUlpglJpVJJ+tC5SZM1mDRgwaTJGmi12mHfi4y700SC40isMJGRkTh8+DBCQkLkjkISYYkVRqPRIDo6Wu4YJCHuTitMe3s7tm/fPuE+moWcxxIrTFdXFw4dOoSuri65o5BEWGIiwbHERIJjiYkEx7PTgps9e7ZD6/v6+uLll1/Gfffdh8DAwHFKRa7EEitMYGAgXnnlFbljkIS4O00kOJaYSHAsMZHgWGIiwbHERIJjiYkExxITCY4lJhIcS0wkOJaYSHAsMZHgWGIiwbHEE8QXX3yBZcuWITQ0FCqVCnl5eXJHotsoLS1FUlISJk+ejLCwMLz++uuy5GCJJwiz2Yy4uDi89tpr0Ol0cseh26iqqsKyZcuQmZmJmpoavPLKK9i2bRt+//vfuzwLb0WcILKyspCVlQUA2Lp1q8xp6HZef/11pKSkYMeOHQCA2NhY1NXVYefOnXj22WddmoUlJrfSfdmMC//pGvb6oMVi/3q6WT/s+5uFBQVg0iTPW/6ciooKrF69eshrGRkZ2LVrF/R6PYKDg+/gT+EYlpjciuckTxSV/hOXzb0jLu+7chX73y8d9fsfRoQgOizotj/HYDAMO+y58b3BYHBpiXlMTG7FSzMZy5csdHrbRzMXQqVSSZxqfLHE5Haiw4Jx3/w5Dm/3P3/6E/hMHdsnGwYGBsJoNA557cKFC/ZlrsQSk1vKXJiKmX7Txrz+PfHRmDs7YszrL1iwAH//+9+HvHb06FGEhoa6dFcaYIknDLPZjJqaGtTU1ODatWswGo2oqalBU1OT3NGE5OnpgSceTodafftdY19vLZY9uMCh99+0aRO++uorvPjii2hsbMQf//hHvPnmm3j++eedjew0lY0fGT8hfP7550hPTx/2+sKFC/H555+7PpCb+KyiGv84XnXLdX6R8zAiQ+92+L1LSkqwbds2NDY2QqfTYePGjXjuueecjeo0llgw5wwXMWO6L7w0k+WOIgSL1Yrf//kwzhkujrj8x8lz8fD/SHNxKmlxd1oggxYL/vy3fyD/d39BW/sFueMI4S61Gk88nA5Pz+GzqT+YMR0/XZgiQyppscQC+frUN+ju6cWkSZ4I/MEMueMIw9/PF0vSfzTktbvUajyRnQ5PD/EvlWCJv8diseBPf/oTHnroIcycOROTJ0/GrFmzkJGRgYKCAli+u9LH1QYtFhz7fycBAA/cmzTiyEKjuzcpFj+MCLF/v/jH8xEU4C9jIunwmPgmly9fxiOPPIKysjIAwN13342goCCcP38e58+fh81mQ2dnJ6ZNm+bybCdq6vG3vx+H91Qv5K7JYYmdcNnch/+z76+YOWMa/uupbKjV7jGGscQ3Wb58OYqKihAcHIx33313yNniCxcuYN++fdi4cSO0Wq1D7/vmHz9Aj/mK07lsNhvMfVdgs9kwedIkTL7Ndb00usHBQajvUkOtmlgF9p46Bf/7Z//LqW356/w7X3/9NYqKiuDh4YEjR45gzpyhV/wEBARg27ZtTr13j/nKqNfyOqr/2jX0X7smyXuRe2CJv/Phhx8CAJYsWTKswHfKe+oUp7flKKwMd/JvhCX+Tn19PQAgLU36OUNnd5MA4ERNA/7293IeC9Oo+C/iO5cvXwZw/UO4pebsMfGNURgArl0bxK/fLpQ6Gk0QPCaWgI+PDwCgu7tb8veW4piYx8I0Gpb4O/Hx8fjggw/wr3/9S/L3duZ4h8fCynInx8ScYvrOyZMnMW/ePHh6eqKmpgZxcXGy5uGxMI3VxJosk9E999yDxx9/HAMDA8jMzMQ///nPIcsvXLiAHTt2oLdXmqmiWxm0WFD2L16dRWPDkfgmly9fxrJly+y3/gUFBeHuu++GwWBAe3u7y67Y4ihMjuBIfBMfHx98+umn2LdvHx544AH09fWhtrYWarUaP/3pT7Fv3z54e3uPe45Jnh7wmarlKExjwpF4ghoYHIQKKnh43CV3FJrgWGIiwXF3mkhwLDGR4FhiIsGxxESCY4mJBMcSEwmOJSYSHEtMJDiWmEhwLDGR4FhiIsGxxESCY4mJBMcSEwmOJSYSHEtMJDiWmEhwLDGR4FhiIsGxxESCY4mJBMcSEwmOJSYSHEtMJDiWmEhwLDGR4FhiIsGxxESCY4mJBMcSEwmOJSYSHEtMJDiWmEhwLDGR4FhiIsGxxESCY4mJBMcSEwmOJSYS3P8HSmE1h8E8WuEAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 287.496x451.5 with 1 Axes>"
      ]
     },
     "execution_count": 4,
     "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])\n",
    "out_qc.measure(qreg_q[3], creg_c[0])\n",
    "\n",
    "out_qc.draw(output='mpl')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "                    ┌───┐          \n",
      "q_0: ───────■───────┤ X ├──■───────\n",
      "            │  ┌───┐└─┬─┘  │       \n",
      "q_1: ───────┼──┤ X ├──■────■────■──\n",
      "     ┌───┐┌─┴─┐└─┬─┘  │  ┌─┴─┐┌─┴─┐\n",
      "q_2: ┤ X ├┤ X ├──■────■──┤ 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": 10,
   "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.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",
    "  \n",
    "    modified_circuits_before.append(modified_qc_before)\n",
    "            \n",
    "\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",
    "    \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",
    "            # If the NOT gate is present in cal_qc for this qubit, continue without creating a new circuit\n",
    "            \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.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",
    "            new_mod_qc = QuantumCircuit(cal_qc.num_qubits)\n",
    "\n",
    "            # Add the modified circuits to the list\n",
    "            if new_depth_b == layers:\n",
    "                # Iterate through the original circuit and replace NOT gates with Identity gates\n",
    "                for instr, qargs, cargs in new_qc:\n",
    "                    if instr.name == 'id':\n",
    "                        # If it's an X gate (NOT gate), replace with Identity gate\n",
    "                        new_mod_qc.append(XGate(), qargs)\n",
    "                    else:\n",
    "                        # If it's not an X gate, keep the original instruction\n",
    "                        new_mod_qc.append(instr, qargs, cargs)\n",
    "                new_circuits1.append(new_mod_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\"h_qubit_{i}.png\")\n",
    "    compiled_new_qc = transpile(new_qc, backend)\n",
    "\n",
    "    file_path = f\"Hadamard_3_17_v2_{i}.qasm\"\n",
    "    with open(file_path, 'w') as file:\n",
    "        file.write(compiled_new_qc.qasm())\n",
    "    print()\n",
    "\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
}