ICC - Ejemplo didactico - Laurel o Yanny.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<img src=\"../images/qiskit-heading.gif\" alt=\"Note: In order for images to show up in this jupyter notebook you need to select File => Trusted Notebook\" width=\"500 px\" align=\"left\">"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# _*Laurel or Yanny?*_ \n",
    "\n",
    "The latest version of this notebook is available on https://github.com/QISKit/qiskit-tutorial.\n",
    "\n",
    "For more information about how to use the IBM Q Experience (QX), consult the [tutorials](https://quantumexperience.ng.bluemix.net/qstage/#/tutorial?sectionId=c59b3710b928891a1420190148a72cce&pageIndex=0), or check out the [community](https://quantumexperience.ng.bluemix.net/qstage/#/community).\n",
    "\n",
    "***\n",
    "### Contributors\n",
    "Adam D. Perruzzi, IBM Q Consulting\n",
    "#### Copyright Attribution\n",
    "The audio files used in this notebook are equalized versions of the recording found at [https://www.vocabulary.com/dictionary/laurel](https://www.vocabulary.com/dictionary/laurel), used under Fair Use provisions and cited as follows: \"Text from Vocabulary.com, Copyright ©1998-2018 Thinkmap, Inc. All rights reserved.\"\n",
    "***\n",
    "#### Prerequisite\n",
    "To run this tutorial, you must have audio-related libraries installed. \n",
    "In MacOS, they can be installed by:\n",
    "```\n",
    "pip install pydub\n",
    "brew install portaudio\n",
    "pip install pyaudio\n",
    "```\n",
    "***"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This program is a 'Hello World' type introduction to the world of quantum computing. It introduces some of the basic concepts of quantum computing in a way that is informative, fun, and engaging.\n",
    "\n",
    "One of the functions of quantum computers - more specifically of the qubits that comprise the quantum computer - is the ability to generate truly random numbers. Classical random number generators are often [pseudo-random](https://www.random.org/randomness), meaning that they are not truly random. Instead, these numbers are often generated using a seed, such as the curent timestamp.\n",
    "\n",
    "One can obtain a random number from a quantum computer by polling a single qubit. In doing so, the qubit 'chooses' to be either a 0 or a 1 in a purely random way. Repeating this polling many times and averaging the results will tend towards an average value of 0.5. That is the method we'll use for this program.\n",
    "\n",
    "And what can we use our random number for? Certainly a great number of important things, but also to try and settle one of the great debates of our time: Laurel or Yanny?\n",
    "\n",
    "In May 2018, an audio clip was discovered on Vocabulary.com accompanying the word [laurel](https://www.vocabulary.com/dictionary/laurel). The clip went viral as part of a tweet, and listeners from around the world weighed in on what they thought they heard. Some heard the intended word 'laurel' while others heard 'yanny'. There are some interesting notes on why this is the case in [this Wired article](https://www.wired.com/story/yanny-and-laurel-true-history).\n",
    "\n",
    "As the jury is still out on exactly what word people are hearing, let's ask a quantum computer what it thinks! The first step is to import the required libraries and initialize the quantum and classical registers that will be used for this program. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<qiskit._measure.Measure at 0x1628b179f60>"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from qiskit import ClassicalRegister, QuantumRegister\n",
    "from qiskit import QuantumCircuit, execute\n",
    "from qiskit.tools.visualization import plot_histogram\n",
    "\n",
    "#Set up the quantum and classical registers, and combine them into a circuit\n",
    "qr = QuantumRegister(1)\n",
    "cr = ClassicalRegister(1)\n",
    "qc = QuantumCircuit(qr, cr)\n",
    "\n",
    "qc.h(qr[0]) #Create a superposition on the single quantum bit\n",
    "qc.measure(qr[0], cr[0]) #Measure the single bit, and store its value in the clasical bit"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We also need to set up the credentials to access the IBM quantum computers. Make sure you've got an IBM Q Experience account set up, and you retrieve your API to paste below (or into the Qconfig.py file, if desired)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "'''from qiskit import register, available_backends, get_backend\n",
    "\n",
    "#Import the config file (Qconfig.py) to retrieve the API token and API url\n",
    "try:\n",
    "    import sys\n",
    "    sys.path.append('../') #Parent directory\n",
    "    import Qconfig\n",
    "    qx_config = {\n",
    "        'APItoken': Qconfig.APItoken,\n",
    "        'url': Qconfig.config['url']}\n",
    "except Exception as e:\n",
    "    print(e)\n",
    "    qx_config = {\n",
    "        'APItoken':'YOUR_TOKEN_HERE',\n",
    "        'url':'https://quantumexperience.ng.bluemix.net/api'}\n",
    "\n",
    "#Setup API\n",
    "register(qx_config['APItoken'], qx_config['url'])'''"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can now run the program on the quantum simulator. You can replace the selected simulation backend with the real quantum computer, if you don't mind a bit of a wait in the queue.\n",
    "\n",
    "You can also change the number of 'shots' taken, or the number of times the program is run. The smaller this number is, the more likely the average will favor one outcome over the other. The larger this number is, the more likely the results will be split 50-50. You'll be able to see the effects of this difference in the eventual audio output."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEZCAYAAACNebLAAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAHlZJREFUeJzt3X2cVmW97/HPN5BdpKE8aCMDAc3oFpIsB5TTfqmJ4wDm0C4zPC8DM2Xnlo0nK6Ozk5NWW1KrndunyKfBMnxqy5xE4JyKjm5FGLAoMIOUZHBCBHzKBxz4nT/WzXgzrGHmxpl138x8368XL9e67ute949xZr5c61rrWooIzMzMWntXsQswM7PS5IAwM7NUDggzM0vlgDAzs1QOCDMzS+WAMDOzVA4IMzNL5YAwM7NUDggzM0vVu9gFvBMDBw6MYcOGFbsMM7MDysqVK1+IiEHt9TugA2LYsGE0NDQUuwwzswOKpL90pJ9PMZmZWSoHhJmZpXJAmJlZKgeEmZmlckCYmVkqB4SZmaVyQJiZWSoHhJmZpXJAmJlZKgeEmZmlckCYmVmqzAJC0gRJT0laL2lWG33OlrRW0hpJd2VVm5mZ7S2Txfok9QJuAKqBRmCFpPqIWJvXpxL4OvCxiNgu6fAsajMzs3RZjSDGAusj4umI2AHMBya36nMhcENEbAeIiOczqs3MzFJkFRCDgY15+425tnxHAUdJ+i9JyyRNSDuQpOmSGiQ1bNmypYvKNTOzrAJCKW3Rar83UAmcApwD3CLp0L3eFDE3IqoiomrQoHafd2FmZvspq4BoBIbk7ZcDz6X0WRARb0XEM8BTJIFhZtamRYsWcfTRR1NRUcGcOXP2ev2OO+5g0KBBHHfccRx33HHccsste7z+8ssvM3jwYGbMmNHStnLlSo499lgqKiqYOXMmEcm/Z7dt20Z1dTWVlZVUV1ezfft2ACKCmTNnUlFRwejRo1m1alXLserq6qisrKSyspK6urqu+BJ0nYjo8j8ko4OngeFAH+B3wKhWfSYAdbntgSSnpAbs67jHH398mFnP1dzcHCNGjIg///nP8eabb8bo0aNjzZo1e/S5/fbb4+KLL27zGDNnzoxzzjlnjz5jxoyJRx99NHbt2hUTJkyIhQsXRkTEV7/61bjqqqsiIuKqq66Kyy67LCIiHnzwwZgwYULs2rUrHnvssRg7dmxERGzdujWGDx8eW7dujW3btsXw4cNj27Ztnfo12B9AQ3Tgd3cmI4iIaAZmAIuBJ4F7ImKNpCsl1ea6LQa2SloL/Br4akRszaI+MzswLV++nIqKCkaMGEGfPn2YMmUKCxYs6PD7V65cyebNmzn99NNb2pqamnj55ZcZN24ckpg6dSoPPPAAAAsWLGDatGkATJs2bY/2qVOnIokTTzyRF198kaamJhYvXkx1dTX9+/fnsMMOo7q6mkWLFnXiV6BrZXYfREQsjIijIuKDEfGdXNvsiKjPbUdEXBoRIyPi2IiYn1VtZnZg2rRpE0OGvH32ury8nE2bNu3V7/7772f06NGcddZZbNyYXC+za9cuvvzlL3PNNdfsdczy8vLUY27evJmysjIAysrKeP755/dZR0frK1W+k9rMDlgRra91AWnPa2LOPPNMNmzYwOrVqznttNNaRgA33ngjkyZN2uMXeEeP2dE69udYpSSTG+XMzLpCeXl5y4gAoLGxkSOPPHKPPgMGDGjZvvDCC/na174GwGOPPcbDDz/MjTfeyKuvvsqOHTs4+OCDueSSS2hsbEw95hFHHEFTUxNlZWU0NTVx+OGH77OO8vJyli5dukf7Kaec0ml//67mEYSZHbDGjBnDunXreOaZZ9ixYwfz58+ntrZ2jz5NTU0t2/X19RxzzDEA/PSnP+XZZ59lw4YNXHvttUydOpU5c+ZQVlbGIYccwrJly4gI5s2bx+TJyX29tbW1LVci1dXV7dE+b948IoJly5bRr18/ysrKqKmpYcmSJWzfvp3t27ezZMkSampqsvjSdAqPIMzsgNW7d2+uv/56ampq2LlzJ+effz6jRo1i9uzZVFVVUVtby3XXXUd9fT29e/emf//+3HHHHe0e96abbuK8887j9ddfZ+LEiUycOBGAWbNmcfbZZ3PrrbcydOhQ7r33XgAmTZrEwoULqaiooG/fvtx+++0A9O/fn8svv5wxY8YAMHv2bPr37981X4wuoLRzZAeKqqqqaGhoKHYZZmYHFEkrI6KqvX4+xWRmZqkcEGZmlsoBYWZmqRwQZmaWygFhZmapHBBmZpbKAWFmZqkcEGZmlsoBYWZmqbzUhpkVzYevnVLsEg5Yv/tK1z8RwSMIMzNL5YAwM7NUDggzM0vlgDAzs1QOCDMzS+WAMDOzVA4IMzNL5YAwM7NUDggzM0uVWUBImiDpKUnrJc1Kef08SVsk/Tb354KsansnFi1axNFHH01FRQVz5sxps999992HJPKfob169WrGjRvHqFGjOPbYY3njjTcAWLlyJcceeywVFRXMnDmT3c8N37ZtG9XV1VRWVlJdXc327dsBiAhmzpxJRUUFo0ePZtWqVS2fUVdXR2VlJZWVldTV1XXFl8DMuqlMAkJSL+AGYCIwEjhH0siUrndHxHG5P7dkUds7sXPnTi6++GIeeugh1q5dy89+9jPWrl27V79XXnmF6667jhNOOKGlrbm5mXPPPZebb76ZNWvWsHTpUg466CAALrroIubOncu6detYt24dixYtAmDOnDmMHz+edevWMX78+JZAeuihh1r6zp07l4suughIAuWKK67g8ccfZ/ny5VxxxRUtoWJm1p6sRhBjgfUR8XRE7ADmA5Mz+uwus3z5cioqKhgxYgR9+vRhypQpLFiwYK9+l19+OZdddhnvfve7W9qWLFnC6NGj+fCHPwzAgAED6NWrF01NTbz88suMGzcOSUydOpUHHngAgAULFjBt2jQApk2btkf71KlTkcSJJ57Iiy++SFNTE4sXL6a6upr+/ftz2GGHUV1d3RI2ZmbtySogBgMb8/Ybc22tfVrSakn3SRqSTWn7b9OmTQwZ8naZ5eXlbNq0aY8+TzzxBBs3buQTn/jEHu1/+tOfkERNTQ0f/ehHufrqq1uOWV5ennrMzZs3U1ZWBkBZWRnPP//8PuvoSH1mZm3JajVXpbRFq/3/DfwsIt6U9EWgDjh1rwNJ04HpAEOHDu3sOguye24gn/T2X3XXrl186Utf4o477tirX3NzM4888ggrVqygb9++jB8/nuOPP573ve99+zxmIXW0V5+Z2b5kNYJoBPJHBOXAc/kdImJrRLyZ2/0xcHzagSJibkRURUTVoEGDuqTYjiovL2fjxrcHRo2NjRx55JEt+6+88gp/+MMfOOWUUxg2bBjLli2jtraWhoYGysvLOfnkkxk4cCB9+/Zl0qRJrFq1ivLychobG1OPecQRR9DU1ARAU1MThx9++D7raK8+M7N9ySogVgCVkoZL6gNMAerzO0gqy9utBZ7MqLb9NmbMGNatW8czzzzDjh07mD9/PrW1tS2v9+vXjxdeeIENGzawYcMGTjzxROrr66mqqqKmpobVq1fz2muv0dzczG9+8xtGjhxJWVkZhxxyCMuWLSMimDdvHpMnJ9M1tbW1LVci1dXV7dE+b948IoJly5bRr18/ysrKqKmpYcmSJWzfvp3t27ezZMkSampqsv9CmdkBKZNTTBHRLGkGsBjoBdwWEWskXQk0REQ9MFNSLdAMbAPOy6K2d6J3795cf/311NTUsHPnTs4//3xGjRrF7Nmzqaqq2iMsWjvssMO49NJLGTNmDJKYNGkSZ5xxBgA33XQT5513Hq+//joTJ05k4sSJAMyaNYuzzz6bW2+9laFDh3LvvfcCMGnSJBYuXEhFRQV9+/bl9ttvB6B///5cfvnljBkzBoDZs2fTv3//rvySmFk3orTz1AeKqqqqyL+vwMwOLH6i3P57J0+Uk7QyIqra6+c7qc3MLJUDwszMUjkgzMwslQPCzMxSOSDMzCyVA8LMzFI5IMzMLJUDwszMUmW1WF/J8Q06+++d3KBjZgcOjyDMzCyVA8LMzFI5IMzMLJUDwszMUjkgzMwslQPCzMxSOSDMzCyVA8LMzFI5IMzMLJUDwszMUjkgzMwsVYcDQtIgSQfntntJ+rykqZIcMmZm3VAhv9x/AVTmtr8DfAW4FPheZxdlZmbFV8hqrkcBv81tnwv8N+BVYA3wpU6uy8zMiqyQgNgJ9JF0FPBSRDybO710cNeUZmZmxVRIQDwE3AMMAHY/EGAksKmzizIzs+IrZA7iAuBB4FbgqlzbQOCbHXmzpAmSnpK0XtKsffQ7S1JIqiqgNjMz62QdHkFExJvA3NxppSOApohY2pH3SuoF3ABUA43ACkn1EbG2Vb9DgJnA4x2ty8zMukYhl7keKuku4A1gfa6tVtK3O/D2scD6iHg6InaQnKKanNLvW8DVuc8wM7MiKuQU083AS8AHgB25tseAz3bgvYOBjXn7jbm2FpI+AgyJiF/s60CSpktqkNSwZcuWjtZuZmYFKmSSejxwZES8JSkAImKLpMM78F6ltEXLi8lpqx8A57V3oIiYC8wFqKqqina6m5nZfipkBPESyaR0C0lDgaYOvLcRGJK3Xw48l7d/CPAhYKmkDcCJQL0nqs3MiqeQgLgFuF/Sx4F3SRoH1JGcemrPCqBS0nBJfYApQP3uFyPipYgYGBHDImIYsAyojYiGAuozM7NOVMgppu+STB7fABwE3Ab8CPhhe2+MiGZJM4DFQC/gtohYI+lKoCEi6vd9BDMzy1ohl7kG8O+5PwWLiIXAwlZts9voe8r+fIaZmXWefQaEpJMi4v/ltk9tq19E/KqzCzMzs+JqbwRxI8nkMSR3UKcJYESnVWRmZiVhnwERER/K2x7e9eWYmVmpKORO6gVttP+888oxM7NSUchlrh9vo/2UTqjDzMxKTLtXMeUuRYXkWRBXtnp5BPCXTq/KzMyKriOXue6+A/pd7Hk3dJCsr/TNTq7JzMxKQLsBERGfB5D0aET8uOtLMjOzUtDefRDDImJDbveXklIvZ42Ipzu7MDMzK672RhC/J1lID5JnQAR7r8waJMtnmJlZN9LefRCH5G0XcsWTmZkd4PxL38zMUrU3B/EweQ/2aUtEnNRpFZmZWUlobw7ilkyqMDOzktPeHERdVoWYmVlpae8U0+ci4s7c9vlt9YuI2zq7MDMzK672TjGdA9yZ2/5cG32C5OlyZmbWjbR3imlS3nZbi/WZmVk3VMgzqZF0KHAGcCTwHPBgRLzYFYWZmVlxFfI8iFOBDcBMYAzwL8AGSeO7pjQzMyumQkYQ1wPTI+Ke3Q2SPgPcAPx9ZxdmZmbFVcid1EcC97dq+0/g/Z1XjpmZlYpCAmIecHGrtoty7WZm1s0UstTGu4CLJF0GbAIGA0cAy7q0QjMzK4pCl9rY7wcGSZoA/JBkafBbImJOq9e/SDJC2Qm8SjLfsXZ/P8/MzN6ZTJbakNSLZDK7GmgEVkiqbxUAd0XEzbn+tcD3gQmd8flmZla4Qu+DOAIYCwwk78FBHVhqYyywfveT5yTNByYDLQERES/n9X8vHVhF1szMuk6HA0LSJ4GfAOuAUcAa4EPAI7S/1MZgYGPefiNwQspnXAxcCvQBTu1obWZm1vkKuYrp28DnI+IjwN9y/50OrOzAe1s/phRSRggRcUNEfBD4GvCN1ANJ0yU1SGrYsmVLx6s3M7OCFBIQQyPi3lZtdcDUDry3ERiSt19OslRHW+YDn0x7ISLmRkRVRFQNGjSoAx9tZmb7o5CAeD43BwHJEhvjgA+SXJXUnhVApaThkvoAU4D6/A6SKvN2zyA5lWVmZkVSyCT1j4F/ILmb+gfAr4FdwPfae2NENEuaASwmCZTbImKNpCuBhoioB2ZIOg14C9gOTCvob2JmZp2qwwEREd/N254naSnw3oh4soPvXwgsbNU2O2/7ko7WYmZmXa/Qy1x7ASfy9nLfvovazKybKuQy19HAA8C7SSady4E3JP1jRPyui+ozM7MiKWSS+jaSu6EHR8RYknsbrsePGzUz65YKCYijgH+PiADI/feHQOU+32VmZgekQgJiIVDbqu1M4MHOK8fMzEpFe8t938nbdzz3AuZLWkmybMYQ4HhgQZdWaGZmRdHeJPX6Vvt/yNteS3Jfg5mZdUPtLfd9RVaFmJlZaSn0PoiPA58juYJpE/CTiPhVVxRmZmbF1eFJakkXAHcDfwV+DjQBd0m6sItqMzOzIipkBHEZUJ1/U5yku0nWZtrvR5GamVlpKuQy1wHkPQEu5ymgf+eVY2ZmpaKQgHgE+L6kvgCS3gtcAzzaFYWZmVlxFRIQXwSOBV6StBl4Efgw8E9dUZiZmRVXh+YgJAl4D3Aa8H5yq7lGRGMX1mZmZkXUoYCIiJD0e+CQXCg4GMzMurlCTjE9QbJgn5mZ9QCFXOa6FFgk6Q6StZh2r9FERHjJbzOzbqaQgPgY8Axwcqv2wM+EMDPrdtoNiNxlrd8AXgVWAf8WEW92dWFmZlZcHZmDuJ7kuQ9PAp8Gru3SiszMrCR0JCAmAqdHxGW57U90bUlmZlYKOhIQ742IJoCI2Aj069qSzMysFHRkkrp3bplvtbGPl/w2M+t+OhIQz7PnVUpbW+0HMKIzizIzs+JrNyAiYlhnfJCkCcAPSZ5tfUtEzGn1+qXABUAzsAU4PyL+0hmfbWZmhSvkTur9JqkXcAPJJPdI4BxJI1t1ewKoiojRwH3A1VnUZmZm6TIJCGAssD4ino6IHcB8YHJ+h4j4dUS8lttdBpRnVJuZmaXIKiAGkyzPsVtjrq0tXwAeSntB0nRJDZIatmzZ0oklmplZvqwCQiltkdKGpHOBKpKHEe39poi5EVEVEVWDBg3qxBLNzCxfIWsxvRONwJC8/XLgudadJJ0G/CtwspfzMDMrrqxGECuASknDJfUBpgD1+R0kfQT4EVAbEc9nVJeZmbUhk4CIiGZgBrCYZE2neyJijaQrJdXmul0DHAzcK+m3kurbOJyZmWUgq1NMRMRCYGGrttl526dlVYuZmbUvq1NMZmZ2gHFAmJlZKgeEmZmlckCYmVkqB4SZmaVyQJiZWSoHhJmZpXJAmJlZKgeEmZmlckCYmVkqB4SZmaVyQJiZWSoHhJmZpXJAmJlZKgeEmZmlckCYmVkqB4SZmaVyQJiZWSoHhJmZpXJAmJlZKgeEmZmlckCYmVkqB4SZmaVyQJiZWSoHhJmZpcosICRNkPSUpPWSZqW8fpKkVZKaJZ2VVV1mZpYuk4CQ1Au4AZgIjATOkTSyVbdngfOAu7KoyczM9q13Rp8zFlgfEU8DSJoPTAbW7u4QERtyr+3KqCYzM9uHrE4xDQY25u035toKJmm6pAZJDVu2bOmU4szMbG9ZBYRS2mJ/DhQRcyOiKiKqBg0a9A7LMjOztmQVEI3AkLz9cuC5jD7bzMz2Q1YBsQKolDRcUh9gClCf0Webmdl+yCQgIqIZmAEsBp4E7omINZKulFQLIGmMpEbgM8CPJK3JojYzM0uX1VVMRMRCYGGrttl52ytITj2ZmVkJ8J3UZmaWygFhZmapHBBmZpbKAWFmZqkcEGZmlsoBYWZmqRwQZmaWygFhZmapHBBmZpbKAWFmZqkcEGZmlsoBYWZmqRwQZmaWygFhZmapHBBmZpbKAWFmZqkcEGZmlsoBYWZmqRwQZmaWygFhZmapHBBmZpbKAWFmZqkcEGZmlsoBYWZmqTILCEkTJD0lab2kWSmv/52ku3OvPy5pWFa1mZnZ3jIJCEm9gBuAicBI4BxJI1t1+wKwPSIqgB8A382iNjMzS5fVCGIssD4ino6IHcB8YHKrPpOButz2fcB4ScqoPjMzayWrgBgMbMzbb8y1pfaJiGbgJWBAJtWZmdleemf0OWkjgdiPPkiaDkzP7b4q6al3WFupGgi8UOwi0uirdxe7BLMslOzPILzjn8MPdKRTVgHRCAzJ2y8HnmujT6Ok3kA/YFvrA0XEXGBuF9VZMiQ1RERVsesw66n8M5jdKaYVQKWk4ZL6AFOA+lZ96oFpue2zgF9FxF4jCDMzy0YmI4iIaJY0A1gM9AJui4g1kq4EGiKiHrgVuFPSepKRw5QsajMzs3TyP9JLk6TpudNpZlYE/hl0QJiZWRu81IaZmaVyQJiZWSoHhJmZpXJAmJm1QdJBkobntnvc78se9xcudZIOKnYNZtbiNOCbknpHxK5iF5M1B0TpeULSTySN64n/YjErMX8FjgFWSKoGUE5xy8qGfwGVnm8AhwEPA3+WdKWkyiLXZNYjRcQTETEWuB+4RNJJkVPs2rLg+yBKVO4ZGl8AZgAfAlaR3G3+84jYXMzazHoaSX2BC4EvAdcCN0bELknqzmHhgCgxuWBQbsnz3W2HAxcDU0lWYbw8Ir5TpBLNegRJ5wIHA38iWdn1LeAfSJYBugW4MiJ2dueQcECUkN3faLnzm0OAzwJrgVUR0ZTrMxZ4PSJ+X8RSzbo1SRXA90hWlT4EWEnyNMxdwPNAJdAAzIiI14tVZ1dzQJSQ3JUSuxc2PAPYAZwJ3BER50vqGxGvFbdKs+4t7x9qh0bEi7m2gRHxQu7CkXcDg4CbgV9GxLXFrLcreZK6hOSdVpoF3BQRk4GFJPMPANdKOrUoxZn1ELlwGAScLmlcru2F3H93AW9ExF+A3wG1xau063kEUWIkHQPcFREfkXQo8DQwNCJelfRHYEpE/La4VZp1X5LOBK4AXgcOzTU/ANweEevz+p1Ccrr38cyLzEhWT5SzjmsEtucmyAYDy3PhcBLQ1+Fg1uW+AtwH3Am8F/g4MInkHqU7I+KfASJiadEqzIgDooTkJqd3Af8BnENyxcT3JX2W5MqJHxexPLNuT9J7gK3A3N2nlXIPMbubZJL65Vxbr4jYWbRCM+KAKAF532yfAsYBXwcOB94HVAN/A+4hGeaaWdc5nuQKwlslzQL+mJsb3AY8srtTTwgH8BxESci7auJHwO8j4vq814YCW7rzpXRmpULS50gmngeTzP89BTxB8nP5l2LWVgwOiBKRG9reBKxrfRNcTxnOmpUCSX2AKuBUklUM+gDNwFci4tli1pY1n2IqHaOBMcBnJX0AuDMiHoaeM5w1K6bcP9KOBBoj4lHg0dyVhCcBH+1p4QAeQZQMSQeTrBpZRfINOQR4heQOzpsiYlMRyzPr9iT9B/Bpkn84/wr4BfDrnvyz54Aoory5h75AOfBMRLwlaSDJ0HYMyeV1X4+IZcWs1aw7k1QD/Az4KPB+knXPTgPeAH4L/HNEvFq8CovDAVFEeQFxAck35GKS9V1WR0RTbuG+sohoLGqhZt2cpE8BH4yIa1q1fxo4JSL+pTiVFZcDogRIGkayMN/HSO7c3Az8F8nVEw0R8beiFWfWA0j6JcmppdqIeKmNPt121da2OCBKjKQTgPOAz5MExKd2r+RqZp1P0hHA1SRzf31Jboq7zasWOCCKZvelq5LelfasW0lXAETE/8q+OrOeI7eCQV+SZz6MBz4BHEeyqsE3I+InRSyvqLyaa5HkXbq6TdL/kTShVZcPkDwP18y6UO4Jon8DnouI24CzSG6W+wnJsx/oqc+H9wiiyCSdDHyO5JtSJJfWbQKmA8dFxIbiVWfWvUn6IHAByf0PzwBrgCfyV23tyRwQJSJ3xVItySJ9O4FbI+L/Frcqs+5N0sNAE/Bn4GySNZf+CvwRuDcilhexvKJzQBSRpMOAE0hGDr/x0+LMsiNpFLA0Igbl9reRjNy/RvIzeW5E/LEnXr20m5fayFje5PQnSb4Rd5Ks8/INSSuA+3cvsWFmXerjwEMAkj4PPB0R90n6GzApIv4IyRxFEWssqh458VJku69Y+jdgHjCRJCC2AmOB/5T0zeKUZtajNJHcnApwFLD7lG4tud+NPXVyejePIDKWu3P6KGBARNwEIOk4kmU1PkZye/9tRSzRrNuTVAmcDvxrrulB4IHc8vrjSG5c7fF6dDpmSdLBuYlogCNIHgCEpFrg2Yh4BvgNcExPXDXSLGPvIbnX4UFJp0bEIySPGg1gxu7J6bR7lHoST1JnQFI/4KfAApIlNLYBZcBqkseKXgLcSDJB1hwR/71IpZr1KJK+Dfw98K2I+F2x6yk1HkFk41WS00ZnAz8HvgUMJ1nB9fck/x9+mPvvd9o4hpl1vm8DfwAelvQ/csvukzfa79E8gshY7tLWzwAzgbdIzn3+iuT50w9HxNYilmfW7Uk6jeQ57y+RPHflSZKfyUuAayPiuiKWV1IcEEWUW8V1OvCPJKOMmRHxWDFrMuvOcs9e+SnJRPSfSJbS+AeSZ09DciXhEuCfIqLHL3XjgCgBuUvpxgBPRcSLxa7HrLvKLcw3ICJekDQEeC0itkoaDGwnef70dcDyiLi+mLWWAgeEmRl73MR6MfA/I2JwsWsqNk9Sm5mxxwrLa4DvFrOWUuERhJmZpfIIwszMUjkgzMwslQPCzMxSOSDMzCyVA8LMzFI5IMzMLNX/B98MBqYADJYkAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1628b1e7ba8>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "backend = 'local_qasm_simulator' #Replace 'ibmq_qasm_simulator' with 'ibmqx5' to run on the quantum computer\n",
    "shots_sim = 100 #Adjust this number as desired, with effects as described above\n",
    "\n",
    "job_sim = execute(qc, backend, shots=shots_sim) #Run job on chosen backend for chosen number of shots\n",
    "stats_sim = job_sim.result().get_counts() #Retrieve results\n",
    "\n",
    "#Select '0' to represent 'laurel'\n",
    "if '0' not in stats_sim.keys():\n",
    "    stats_sim['laurel'] = 0\n",
    "else:\n",
    "    stats_sim['laurel'] = stats_sim.pop('0')\n",
    "\n",
    "#Which leaves '1' to represent 'yanny'\n",
    "if '1' not in stats_sim.keys():\n",
    "    stats_sim['yanny'] = 0\n",
    "else:\n",
    "    stats_sim['yanny'] = stats_sim.pop('1')\n",
    "\n",
    "plot_histogram(stats_sim)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This graph shows how often the qubit was measured to be a 0, and how often it was measured to be a 1. We can use these results to help us find the answer to the laurel/yanny debate.\n",
    "\n",
    "We can do this by starting with two separate audio tracks, one for laurel and one for yanny. These tracks have been equalized to highlight one name or the other. First we have to install the required libraries. If you're having trouble, follow the guide for pydub installation available [here](https://github.com/jiaaro/pydub#installation)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "from pydub import AudioSegment\n",
    "from pydub.playback import play\n",
    "\n",
    "AudioSegment.converter = r\"D:\\\\Program Files\\\\ffmpeg\\\\bin\\\\ffmpeg.exe\"\n",
    "#Import two tracks\n",
    "laurel = AudioSegment.from_wav('laurel_or_yanny_audio_files/laurel.wav')\n",
    "yanny = AudioSegment.from_wav('laurel_or_yanny_audio_files/yanny.wav')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "If you've never been able to hear one name or the other, take a listen to either or both of the two tracks below. I can't guarantee you'll hear both distinctly, but these tracks will get you as close as possible!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "play(laurel) #Listen to the laurel-specific track"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "play(yanny) #Listen to the yanny-specific track"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now we can apply our results from the quantum computer to these audio tracks. We can do this by proportionally reducing the volume for the less common track, and increasing the volume for the more common track."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "ename": "NameError",
     "evalue": "name 'stats_sim' is not defined",
     "output_type": "error",
     "traceback": [
      "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[1;31mNameError\u001b[0m                                 Traceback (most recent call last)",
      "\u001b[1;32m<ipython-input-5-92f78a7d4d8f>\u001b[0m in \u001b[0;36m<module>\u001b[1;34m()\u001b[0m\n\u001b[0;32m      1\u001b[0m \u001b[1;31m#Modify the volumes based on the results of the experiment\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m----> 2\u001b[1;33m \u001b[0mlaurel\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mlaurel\u001b[0m \u001b[1;33m+\u001b[0m \u001b[1;33m(\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;36m100\u001b[0m\u001b[1;33m*\u001b[0m\u001b[0mstats_sim\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;34m'laurel'\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m/\u001b[0m\u001b[0mshots_sim\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m-\u001b[0m\u001b[1;36m50\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;31m#Laurel\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m      3\u001b[0m \u001b[0myanny\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0myanny\u001b[0m \u001b[1;33m+\u001b[0m \u001b[1;33m(\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;36m100\u001b[0m\u001b[1;33m*\u001b[0m\u001b[0mstats_sim\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;34m'yanny'\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m/\u001b[0m\u001b[0mshots_sim\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m-\u001b[0m\u001b[1;36m50\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;31m#Yanny\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m      4\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m      5\u001b[0m \u001b[1;31m#Mix the two together and play the result\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;31mNameError\u001b[0m: name 'stats_sim' is not defined"
     ]
    }
   ],
   "source": [
    "#Modify the volumes based on the results of the experiment\n",
    "laurel = laurel + ((100*stats_sim['laurel']/shots_sim)-50) #Laurel\n",
    "yanny = yanny + ((100*stats_sim['yanny']/shots_sim)-50) #Yanny\n",
    "\n",
    "#Mix the two together and play the result\n",
    "mixed = laurel.overlay(yanny)\n",
    "#mixed = yanny.overlay(laurel)\n",
    "play(mixed)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "If you're having trouble getting the above clip to play in the notebook, you can run the code below to download the clip using your audio player of choice. The audio file will be placed in the same folder as the other audio files."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "mixed.export('laurel_or_yanny_audio_files/quantumLaurelYanny.wav', format='wav')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "If the results above are roughly 50-50, the combination will sound similar to the original track. If the experiment was conducted with fewer shots, the clip may lean more towards one sound or another. Go back earlier in the program to change the number of shots, re-run, and see if you can get one or the other to sound stronger or weaker.\n",
    "\n",
    "You can also change the backend to run this on an actual quantum computer, instead of a simulation. It may take a while to run, as there may be a queue for the machine.\n",
    "\n",
    "***\n",
    "\n",
    "And there you have it! Clearly the quantum computer can't make up its mind what it hears either, so the question will just have to go unanswered for the time being :)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "print(\"Installed packages are as the following\")\n",
    "!python --version\n",
    "print()\n",
    "!conda list 'qiskit|IBMQuantumExperience|numpy|scipy'"
   ]
  }
 ],
 "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.6.4"
  },
  "toc": {
   "colors": {
    "hover_highlight": "#DAA520",
    "running_highlight": "#FF0000",
    "selected_highlight": "#FFD700"
   },
   "moveMenuLeft": true,
   "nav_menu": {
    "height": "48px",
    "width": "252px"
   },
   "navigate_menu": true,
   "number_sections": true,
   "sideBar": true,
   "threshold": 4,
   "toc_cell": false,
   "toc_section_display": "block",
   "toc_window_display": false
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}