Deutsch-Jozsa-Algorithmus

aus Wikipedia, der freien Enzyklopädie
Wechseln zu: Navigation, Suche

Der Algorithmus von Deutsch ist ein Algorithmus für Quantencomputer, mit dem man bestimmen kann, ob eine auf einem Bit operierende Funktion konstant oder balanciert ist. Diese Aufgabenstellung ist unter dem Namen Problem von Deutsch bekannt. Der Algorithmus von Deutsch ist zwar kaum von praktischem Nutzen, er war jedoch historisch der erste Quantenalgorithmus, der eine Aufgabenstellung nachweisbar schneller löst als ein klassischer Algorithmus und damit die theoretischen Möglichkeiten von Quantencomputern aufzeigt.

Eine Verallgemeinerung ist der Deutsch-Jozsa-Algorithmus. Dabei wird das Problem von Deutsch auf mehrere Bits übertragen.

Die Algorithmen wurden nach ihren Urhebern David Deutsch und Richard Jozsa benannt.

Geschichte[Bearbeiten]

In einer Arbeit von 1985 beschäftigte sich David Deutsch mit dem Spezialfall des Problems (Problem von Deutsch, Anzahl der Eingabebits: 1).[1]

1992 wurde die Idee durch David Deutsch und Richard Jozsa verallgemeinert (Problem von Deutsch-Jozsa, beliebig viele Eingabebits).[2]

Der von Deutsch ursprünglich vorgeschlagene Algorithmus war de facto nicht deterministisch. Er war mit einer Wahrscheinlichkeit von 0,5 erfolgreich. Der originale Deutsch-Jozsa-Algorithmus war zwar deterministisch, benötigte jedoch im Gegensatz zum Algorithmus von Deutsch zwei Funktionsauswertungen statt nur einer.

Durch R. Cleve, A. Ekert, C. Macchiavello und M. Mosca wurden 1998 weitere Verbesserungen vorgenommen, so dass jetzt nur noch eine Funktionsauswertung benötigt wird und der Algorithmus dennoch deterministisch bleibt.[3]

Der Deutsch-Jozsa-Algorithmus diente als Inspiration für den Shor-Algorithmus und den Grover-Algorithmus, zwei der revolutionärsten Quantenalgorithmen.[4][5]

Das Problem von Deutsch[Bearbeiten]

Problemformulierung[Bearbeiten]

Gegeben sei eine Funktion f: \lbrace0,1\rbrace \to \lbrace0,1\rbrace. Es gilt nun herauszufinden, ob die gegebene Funktion konstant ist (f(0) = f(1)) oder balanciert/gleichgewichtig (f(0) \neq f(1)). Dabei ist zu beachten, dass f nur als Black Box gegeben ist, die genaue Definition also unbekannt ist - es steht lediglich ein Bauteil zur Verfügung, welches zu einem Bit b den Wert f(b) berechnet. Um die Bedeutung des Algorithmus von Deutsch zu unterstreichen, sollte davon ausgegangen werden, dass die Benutzung dieses Bauteils teuer ist.

Zur Veranschaulichung kann man sich vorstellen, dass man entscheiden soll, ob eine gegebene Münze fair (Kopf auf der einen Seite, Zahl auf der anderen) oder unfair (Kopf oder Zahl auf beiden Seiten der Münze) ist.

Klassische Lösung[Bearbeiten]

Auf klassischem Wege bleibt nichts weiter übrig als die Funktion sowohl für f(0) als auch f(1) auszuwerten und zu vergleichen. Das entspricht dem Betrachten einer Münze auf beiden Seiten. Die Black Box (oder auch Orakel) wird also zwei mal benötigt. Im Folgenden wird gezeigt, dass der Quantenalgorithmus von Deutsch mit nur einem Aufruf auskommt. Dies bringt einen Vorteil, wenn man bedenkt, dass nach der Annahme ein Aufruf sehr teuer ist (jedoch sind die Kosten asymptotisch gleich).

Der Quantenalgorithmus[Bearbeiten]

Quantenschaltkreis, der den Algorithmus von Deutsch implementiert

Der Quantenalgorithmus wendet die gegebene Funktion auf eine Superposition der beiden möglichen Eingaben an. Durch geschickte Auswertung erhält er somit mit nur einmaliger Anwendung des Orakels ausreichend Information über f(0) und f(1).

Da in der Quanteninformatik alle Rechenschritte umkehrbar sein müssen, wird hier eine spezielle Variante von f benötigt, unter Verwendung von \oplus (Addition mit anschließendem Modulo 2) beschrieben durch die Abbildung

U_f: \left|x\right\rangle \left|y\right\rangle \mapsto \left|x\right\rangle \left|f(x)\oplus y\right\rangle

Der Algorithmus verwendet ein Register von zwei Qubits \left|x\right\rangle \cdot \left|y\right\rangle (mit \cdot dem Tensorprodukt, welches gewöhnlich mit \otimes dargestellt wird) und besteht aus folgenden Schritten:

  1. Initialisiere wie folgt:
    \begin{align}\left|x\right\rangle\left|y\right\rangle &= \left|0\right\rangle\left|1\right\rangle\\
&= \left|a\right\rangle\end{align}
  2. Wende die Hadamard-Transformation H auf beide Qubits an:
    \begin{align}\left|x\right\rangle\left|y\right\rangle &\rightarrow H\left|x\right\rangle H\left|y\right\rangle\\
& = \frac{1}{\sqrt{2}}\bigl(\left|0\right\rangle +\left|1\right\rangle\bigr)\cdot\frac{1}{\sqrt{2}}\bigl(\left|0\right\rangle -\left|1\right\rangle\bigr)\\
& = \frac{1}{2}\bigl(\left|0\right\rangle\left|0\right\rangle - \left|0\right\rangle\left|1\right\rangle + \left|1\right\rangle\left|0\right\rangle - \left|1\right\rangle\left|1\right\rangle\bigr)\\
& = \left|b\right\rangle\end{align}
  3. Werte U_f aus:
    \begin{align}\left|x\right\rangle\left|y\right\rangle &\rightarrow U_f\left|x\right\rangle\left|y\right\rangle\\
&=\frac{1}{2}\Bigl(\left|0\right\rangle\left|0\oplus f(0)\right\rangle - \left|0\right\rangle\left|1\oplus f(0)\right\rangle + \left|1\right\rangle\left|0 \oplus f(1)\right\rangle - \left|1\right\rangle\left|1 \oplus f(1)\right\rangle\Bigr)\\
&= \frac{1}{2}\Bigl(\left|0\right\rangle \cdot \bigl(\left|f(0)\right\rangle - \left|1 \oplus f(0)\right\rangle\bigr) + \left|1\right\rangle \cdot \bigl(\left|f(1)\right\rangle - \left|1 \oplus f(1)\right\rangle\bigr)\Bigr)\\
&= \frac{1}{2}\Bigl((-1)^{f(0)}\left|0\right\rangle\cdot\bigl(\left|0\right\rangle-\left|1\right\rangle\bigr)+(-1)^{f(1)}\left|1\right\rangle\cdot\bigl(\left|0\right\rangle-\left|1\right\rangle\bigr)\Bigr)\\
&=\frac{1}{2}\bigl((-1)^{f(0)}\left|0\right\rangle + (-1)^{f(1)}\left|1\right\rangle\bigr)\cdot\bigl(\left|0\right\rangle - \left|1\right\rangle\bigr)\\
&=\left|c\right\rangle\end{align}
    \left|c\right\rangle ist im konstanten Fall \frac{1}{2}\Bigl(\pm\bigl(\left|0\right\rangle + \left|1\right\rangle\bigr)\Bigr)\cdot\bigl(\left|0\right\rangle - \left|1\right\rangle\bigr) und sonst \frac{1}{2}\Bigl(\pm\bigl(\left|0\right\rangle - \left|1\right\rangle\bigr)\Bigr)\cdot\bigl(\left|0\right\rangle - \left|1\right\rangle\bigr)
  4. Wende die Hadamard-Transformation H auf das erste Qubit an:
    \left|x\right\rangle \rightarrow H\left|x\right\rangle
  5. Miss das erste Qubit:
    Hat \left|x\right\rangle den Wert \left|0\right\rangle, so ist f konstant, ansonsten balanciert.

Der Trick besteht also darin, dass wir die Funktionswerte in die Amplitude verlagern.

Realisierung[Bearbeiten]

Über die erste experimentelle Realisierung wurde von S. Gulde in Nature 421, 48 (2003) berichtet. Die dafür benötigten Qubits wurden durch den elektronischen und vibratorischen Freiheitsgrad eines Ca^{+}-Ions in einer Falle implementiert.

Das Problem von Deutsch-Jozsa[Bearbeiten]

Problemformulierung[Bearbeiten]

Nun wird die Funktion f auf n Eingabebits verallgemeinert: f: \lbrace0,1\rbrace^n \to \lbrace0,1\rbrace Es wird zugesichert, dass die Funktion entweder konstant (alle Eingaben werden auf ein und denselben Wert abgebildet) oder balanciert (die Hälfte der Eingaben wird auf 0 und die andere Hälfte auf 1 abgebildet) ist. Herauszufinden ist nun, welche der beiden Alternativen zutrifft. Erneut ist zu beachten, dass f nur als Black Box bzw. Orakel gegeben ist.

Klassische Lösung[Bearbeiten]

Ein klassischer Computer müsste die Funktion im schlimmsten Fall für mehr als die Hälfte aller möglichen Eingaben auswerten. So kann es beispielsweise passieren, dass selbst bei noch so geschickter Wahl die ersten 2^{n-1} (halb so viel wie insgesamt möglich) getesteten Eingaben zum Wert 0 führen und trotzdem noch nicht bekannt ist, welche Alternative gilt. Ergibt auch die nächste Auswertung 0, so ist die Funktion konstant, ergibt sie hingegen 1, ist sie balanciert.

Ein klassischer Algorithmus benötigt also im worst case 2^{n-1}+1 Auswertungen von f. Wie im Folgenden gezeigt wird, benötigt der Algorithmus von Deutsch-Jozsa nur eine einzige Auswertung. Dies stellt einen exponentiellen Laufzeitgewinn dar, was im Gegensatz zum Problem / Algorithmus von Deutsch eine echte Verbesserung ist.

Der Quantenalgorithmus[Bearbeiten]

Quantenschaltkreis, der den Algorithmus von Deutsch-Jozsa implementiert

Der Quantenalgorithmus wendet die gegebene Funktion auf eine Superposition aller möglichen Eingaben an. Durch geschickte Auswertung erhält er somit mit nur einmaliger Anwendung des Orakels ausreichend Information über alle Funktionswerte.

Auf Grund der zu erhaltenden Umkehrbarkeit wird wieder eine spezielle Variante von f benötigt, beschrieben durch die Abbildung

U_f: \left|x\right\rangle \left|y\right\rangle \mapsto \left|x\right\rangle \left|f(x)\oplus y\right\rangle

wobei \left|x\right\rangle die Eingabe der Länge n ist.

Der Algorithmus benutzt n+1 Qubits und durchläuft folgende Schritte:

  1. Initialisierung: die ersten n Bits auf \left|0\right\rangle und das letzte Bit auf \left|1\right\rangle setzen:
    \begin{align}\left|x\right\rangle\left|y\right\rangle &= \left|0\right\rangle^{\otimes n} \left|1\right\rangle\\& = \left|a\right\rangle\end{align}
  2. Wende die Hadamard-Transformation H_{n+1} (also H_2 auf alle Qubits) an:
    \begin{align}\left|x\right\rangle\left|y\right\rangle &\rightarrow H_{n+1}\left|x\right\rangle\left|y\right\rangle\\
&=\left(\frac{1}{\sqrt{2^{n}}}\sum_{x=0}^{2^n-1} \left|x\right\rangle\right) \cdot \frac{1}{\sqrt{2}} \bigl(\left|0\right\rangle - \left|1\right\rangle \bigr)\\
&=\frac{1}{\sqrt{2^{n+1}}}\sum_{x=0}^{2^n-1} \left|x\right\rangle \bigl(\left|0\right\rangle - \left|1\right\rangle \bigr)\\
&=\left|b\right\rangle\end{align}
  3. Werte U_f aus:
    \begin{align}\left|x\right\rangle\left|y\right\rangle &\rightarrow U_f\left|x\right\rangle\left|y\right\rangle\\
&= \frac{1}{\sqrt{2^{n+1}}}\sum_{x=0}^{2^n-1} \left|x\right\rangle \bigl(\left|f(x)\right\rangle - \left|1\oplus f(x)\right\rangle \bigr) \\
&= \frac{1}{\sqrt{2^{n+1}}}\sum_{x=0}^{2^n-1} (-1)^{f(x)} \left|x\right\rangle \bigl(\left|0\right\rangle - \left|1\right\rangle \bigr) \\
&= \left(\frac{1}{\sqrt{2^{n}}}\sum_{x=0}^{2^n-1} (-1)^{f(x)} \left|x\right\rangle\right) \cdot \frac{1}{\sqrt{2}}\bigl(\left|0\right\rangle - \left|1\right\rangle \bigr) \\ 
&= \left|c\right\rangle\end{align}
  4. Wende die Hadamard-Transformation auf \left|x\right\rangle an:
    \begin{align} \left|x\right\rangle\left|y\right\rangle &\rightarrow \bigl(H_n\left|x\right\rangle\bigr)\left|y\right\rangle\\ 
&= \left(\frac{1}{2^n}\sum_{z=0}^{2^n-1}\sum_{x=0}^{2^n-1} (-1)^{\vec x\cdot \vec z +f(x)}\left|z\right\rangle\right)\cdot
\frac{1}{\sqrt{2}}\bigl(\left|0\right\rangle-\left|1\right\rangle\bigr) \\
&= \left|d\right\rangle \end{align}
  5. Miss das Register \left|x\right\rangle:
    Ist es \left|0\ldots 0\right\rangle, so ist f konstant, ansonsten balanciert.

Die Interpretation des Messergebnisses ist folgendermaßen einzusehen: Ist f balanciert, so gleichen sich die Vorzeichen für \left|z\right\rangle=\left|0\right\rangle aus, so dass mit einer Wahrscheinlichkeit von 0 der Wert \left|0\right\rangle gemessen wird. Im anderen Fall, also wenn f konstant ist, sind genau die Amplituden für \left|z\right\rangle\neq\left|0\right\rangle gleich 0, da sich für solche z die x immer so einteilen lassen, dass die Hälfte skalarmultipliziert mit z einen geraden Wert ergibt, die andere Hälfte aber einen ungeraden.

Einzelnachweise[Bearbeiten]

  1. David Deutsch: The Church-Turing principle and the universal quantum computer. (PDF) In: Proceedings of the Royal Society of London A. 400, 1985, S. 97.
  2. David Deutsch und Richard Jozsa: Rapid solutions of problems by quantum computation. In: Proceedings of the Royal Society of London A. 439, 1992, S. 553.
  3. R. Cleve, A. Ekert, C. Macchiavello und M. Mosca: Quantum algorithms revisited. (PDF) In: Proceedings of the Royal Society of London A. 454, 1998, S. 339-354.
  4. Lov K. Grover: A fast quantum mechanical algorithm for database search. (PDF) In: Proceedings of the Twenty-Eighth Annual ACM Symposium on Theory of Computing. 1996, S. 212-219.
  5. Peter W. Shor: Algorithms for Quantum Computation: Discrete Logarithms and Factoring. (PDF) In: IEEE Symposium on Foundations of Computer Science. 1994, S. 124-134.

Literatur[Bearbeiten]

  • Matthias Homeister: Quantum computing verstehen. Vieweg, Wiesbaden 2005, ISBN 3-528-05921-4, S. 33–37, 62–65