D.1 Deutsch-Jozsa algorithm

Similar documents
D.1 Deutsch-Jozsa algorithm

92 CHAPTER III. QUANTUM COMPUTATION. Figure III.11: Diagram for swap (from NC).

Computer arithmetic. Intensive Computation. Annalisa Massini 2017/2018

D.5 Quantum error correction

Lecture 4: Elementary Quantum Algorithms

Evaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models

CSE 599d - Quantum Computing When Quantum Computers Fall Apart

QUANTUM COMPUTATION. Exercise sheet 1. Ashley Montanaro, University of Bristol H Z U = 1 2

Shor s Prime Factorization Algorithm

Impossibility of a Quantum Speed-up with a Faulty Oracle

Quantum Computing: Foundations to Frontier Fall Lecture 3

Introduction to Quantum Computing

C. QUANTUM INFORMATION 99

Compute the Fourier transform on the first register to get x {0,1} n x 0.

Fourier Sampling & Simon s Algorithm

UNM Physics 452/581: Introduction to Quantum Information, Problem Set 3, Fall 2007

Short introduction to Quantum Computing

HENSEL S LEMMA KEITH CONRAD

AMS10 HW1 Grading Rubric

Nothing Here. Fast Quantum Algorithms

Fault Tolerant Quantum Computing Robert Rogers, Thomas Sylwester, Abe Pauls

QLang: Qubit Language

Topic: Lower Bounds on Randomized Algorithms Date: September 22, 2004 Scribe: Srinath Sridhar

Lecture 1: Introduction to Quantum Computing

Ph 219b/CS 219b. Exercises Due: Wednesday 4 December 2013

LOWELL WEEKI.Y JOURINAL

Lecture 21: Quantum Communication

arxiv: v2 [quant-ph] 2 Aug 2012

Lecture 1: Overview of quantum information

Real Analysis 1 Fall Homework 3. a n.

Introduction to Quantum Algorithms Part I: Quantum Gates and Simon s Algorithm

Approximating min-max k-clustering

Deutsch Algorithm on Classical Circuits

16.2. Infinite Series. Introduction. Prerequisites. Learning Outcomes

MANY BILLS OF CONCERN TO PUBLIC

AI*IA 2003 Fusion of Multiple Pattern Classifiers PART III

Quantum Communication Complexity

Lecture 9: Connecting PH, P/poly and BPP

On the Rank of the Elliptic Curve y 2 = x(x p)(x 2)

The Euler Phi Function

MATH 2710: NOTES FOR ANALYSIS

Parallel Quantum-inspired Genetic Algorithm for Combinatorial Optimization Problem

LOWELL. MICHIGAN. WEDNESDAY, FEBRUARY NUMllEE 33, Chicago. >::»«ad 0:30am, " 16.n«l w 00 ptn Jaekten,.'''4snd4:4(>a tii, ijilwopa

6.896 Quantum Complexity Theory September 18, Lecture 5

Simon s algorithm (1994)

Pseudorandom Sequence Generation

Ph 219b/CS 219b. Exercises Due: Wednesday 22 February 2006

I Poles & zeros. I First-order systems. I Second-order systems. I E ect of additional poles. I E ect of zeros. I E ect of nonlinearities

Introduction to Quantum Information Processing

Richard Cleve David R. Cheriton School of Computer Science Institute for Quantum Computing University of Waterloo

On Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm

The Arm Prime Factors Decomposition

arxiv: v1 [quant-ph] 3 Feb 2015

Outline. EECS150 - Digital Design Lecture 26 Error Correction Codes, Linear Feedback Shift Registers (LFSRs) Simple Error Detection Coding

1-way quantum finite automata: strengths, weaknesses and generalizations

Cryptanalysis of Pseudorandom Generators

15-451/651: Design & Analysis of Algorithms October 23, 2018 Lecture #17: Prediction from Expert Advice last changed: October 25, 2018

John Weatherwax. Analysis of Parallel Depth First Search Algorithms

Simulation of quantum computers with probabilistic models

Principal Components Analysis and Unsupervised Hebbian Learning

16.2. Infinite Series. Introduction. Prerequisites. Learning Outcomes

0.6 Factoring 73. As always, the reader is encouraged to multiply out (3

1. INTRODUCTION. Fn 2 = F j F j+1 (1.1)

Principles of Computed Tomography (CT)

Quantum Communication

Extension of Minimax to Infinite Matrices

DISCRIMINANTS IN TOWERS

DRAFT - do not circulate

Quantum Computing Lecture 6. Quantum Search

Quantum Computer Algorithms: basics

Quantum Circuits and Algorithms

arxiv: v1 [physics.data-an] 26 Oct 2012

Single qubit + CNOT gates

Quantum Computing Lecture Notes, Extra Chapter. Hidden Subgroup Problem

Day 3 Selected Solutions

General Linear Model Introduction, Classes of Linear models and Estimation

Introduction to Quantum Information Processing QIC 710 / CS 768 / PH 767 / CO 681 / AM 871

Analysis of execution time for parallel algorithm to dertmine if it is worth the effort to code and debug in parallel

d A L. T O S O U LOWELL, MICHIGAN. THURSDAY, DECEMBER 5, 1929 Cadillac, Nov. 20. Indignation

Adam Paweł Zaborski. 8 Plasticity. reloading. 1. Bauschinger s effect. 2. unchanged yielding limit. 3. isotropic hardening

PHYS 301 HOMEWORK #9-- SOLUTIONS

Seminar 1. Introduction to Quantum Computing

An Introduction to Quantum Computation and Quantum Information

Proof Nets and Boolean Circuits

Unit 1 - Computer Arithmetic

b) (5 points) Give a simple quantum circuit that transforms the state

An Introduction to Quantum Computing for Non-Physicists

Convolutional Codes. Lecture 13. Figure 93: Encoder for rate 1/2 constraint length 3 convolutional code.

Introduction to Quantum Computation

Title. Author(s)Okamoto, Ryo; Hofmann, Holger F.; Takeuchi, Shigeki; CitationPhysical Review Letters, 95: Issue Date

Elementary Analysis in Q p

Positive decomposition of transfer functions with multiple poles

Introduction to Quantum Information Processing CS 467 / CS 667 Phys 667 / Phys 767 C&O 481 / C&O 681

Pythagorean triples and sums of squares

MATH342 Practice Exam

arxiv: v3 [quant-ph] 30 Sep 2009

Planck Constant and the Deutsch-Jozsa Algorithm

Tutorial on Quantum Computing. Vwani P. Roychowdhury. Lecture 1: Introduction

A CONCRETE EXAMPLE OF PRIME BEHAVIOR IN QUADRATIC FIELDS. 1. Abstract

#A47 INTEGERS 15 (2015) QUADRATIC DIOPHANTINE EQUATIONS WITH INFINITELY MANY SOLUTIONS IN POSITIVE INTEGERS

Transcription:

4 CHAPTER III. QUANTUM COMPUTATION Figure III.: Quantum circuit for Deutsch algorithm. [fig. from Nielsen & Chuang (00)] D Quantum algorithms D. Deutsch-Jozsa algorithm D..a Deutsch s algorithm In this section you will encounter your first examle of a quantum algorithm that can comute faster than a classical algorithm for the same roblem. This is a simlified version of Deutsch s original algorithm, which shows how it is ossible to extract global information about a function by using quantum arallelism and interference (Fig. III.). 8 Suose we have a function f :!, asinsec.c.5. Thegoalisto determine whether f(0) = f() with a single function evaluation. This is not a very interesting roblem (since there are only four such functions), but it is a warmu for the Deutsch-Jozsa algorithm. Simle as it is, it could be exensive to decide on a classical comuter. For examle, suose f(0) = the billionth bit of and f() = the billionth bit of e. Thentheroblemis to decide if the billionth bits of and e are the same. It is mathematically simle, but comutationally comlex. To see how we might solve this roblem, suose we have a quantum gate array U f for f; thatis,u f xi yi = xi y f(x)i. In articular, U f xi 0i = 8 This is the 998 imrovement by Cleve et al. to Deutsch s 985 algorithm (Nielsen & Chuang, 00,. 59).

D. QUANTUM ALGORITHMS 5 xi f(x)i and U f xi i = xi f(x)i. Usually we set y =0togettheresult f(x)i, buthereyouwillseeanalicationinwhichwewanty =. Now consider the result of alying U f to xi in the data register and to the suerosition i = ( 0i i) inthetargetregister. U f xi i = xi f(x)i xi f(x)i = xi[ f(x)i f(x)i]. Now the rightmost square bracket is 0i i if f(x) =0or i 0i if f(x) =. Therefore,wecanwrite U f xi i = xi( ) f(x) ( 0i i) =( ) f(x) xi i. (III.) [Here, ( ) x is an abbreviation for ( ) x when we want to emhasize that the sign is all that matters.] Since U f xi i =( ) f(x) xi i, theresultof alying it to an equal suerosition of x =0andx =is: U f xi i = ( ) f(x) xi i. x If f is a constant function, then f(0) = f(), and the summation is ± ( 0i+ i) i = ± +i i because both comonents have the same sign.. On the other hand, if f(0) 6= f(), then the summation is ± ( 0i i) i = ± i i because the comonents have oosite signs. That is, a constant function gives the 0i and i comonents of the data qubit the same hase, and otherwise gives them the oosite hase. Therefore, we can determine whether the function is constant or not by measuring the first qubit in the sign basis; we get +i if f(0) = f() and i otherwise. With this background, we can state Deutsch s algorithm. x algorithm Deutsch: Initial state: Begin with the qubits 0 i def = 0i. Suerosition: Transform it to a air of suerositions i def = ( 0i + i) ( 0i i) = + i. (III.)

6 CHAPTER III. QUANTUM COMPUTATION by a air of Hadamard gates. Recall that H 0i = ( 0i + i) = +i and H i = ( 0i i) = i. Function alication: Next aly U f to i = + i. As we ve seen, U f xi 0i = xi 0 f(x)i = xi f(x)i, and U f xi i = xi f(x)i = xi f(x)i. Therefore,exandEq.III.andalyU f : i def = U f i = U f ale ( 0i + i) ( 0i i) = [U f 00i U f 0i + U f 0i U f i] = [ 0,f(0)i 0, f(0)i +,f()i, f()i] There are two cases: f(0) = f() and f(0) 6= f(). Equal (constant function): If f(0) = f(), then i = [ 0,f(0)i 0, f(0)i +,f(0)i, f(0)i] = [ 0i( f(0)i f(0)i)+ i( f(0)i f(0)i)] = ( 0i + i)( f(0)i f(0)i) = ± ( 0i + i)( 0i i) = ± ( 0i + i) i = + i. The last line alies because global hase (including ±) doesn tmatter. Unequal (balanced function): If f(0) 6= f(), then i = [ 0,f(0)i 0, f(0)i +, f(0)i,f(0)i]

D. QUANTUM ALGORITHMS 7 = [ 0i( f(0)i f(0)i)+ i( f(0)i f(0)i)] = [ 0i( f(0)i f(0)i) i( f(0)i f(0)i)] = ( 0i i)( f(0)i f(0)i) = ± ( 0i i)( 0i i) = ± ( 0i i) i = i Clearly we can discriminate between the two cases by measuring the first qubit in the sign basis. In articular, note that in the unequal case, the i comonent has the oosite hase from the 0i comonent. Measurement: Therefore we can determine whether f(0) = f() or not by measuring the first bit of i in the sign basis, which we can do with the Hadamard gate (recall H +i = 0i and H i = i): 3 i def = (H I) i ± 0i i, if f(0) = f() = ± i i, if f(0) 6= f() = ± f(0) f()i i. Notice that the information we need is in the data register, not the target register. This technique is called hase kick-back (i.e., kicked back into the hase of the data register). In conclusion, we can determine whether or not f(0) = f() with a single evaluation of f, whichisquiteremarkable.ine ect,weareevaluatingf on asuerositionof 0i and i and determining how the results interfere with each other. As a result we get a definite (not robabilistic) determination of aglobalroertywithasingleevaluation. This is a clear examle where a quantum comuter can do something faster than a classical comuter. However, note that U f has to uncomute

8 CHAPTER III. QUANTUM COMPUTATION Figure III.3: Quantum circuit for Deutsch-Jozsa algorithm. [fig. from NC] f, whichtakesasmuchtimeascomutingit,butwewillseeothercases (Deutsch-Jozsa) where the seedu is much more than. D..b The Deutsch-Jozsa algorithm The Deutsch-Jozsa algorithm is a generalization of the Deutsch algorithm to n bits, which they ublished it in 99; here we resent the imroved version of Nielsen & Chuang (00,. 59). This is the roblem: Suose we are given an unknown function f : n! in the form of a unitary transform U f L(H n+, H) (Fig. III.3). We are told only that f is either constant or balanced, whichmeansthatitis0on half its domain and on the other half. Our task is to determine into which class the given f falls. Consider first the classical situation. We can try di erent inut bit strings x. We might (if we re lucky) discover after the second query of f that it is not constant. But we might require as many as n /+queriestoanswer the question. So we re facing O( n )functionevaluations. algorithm Deutsch-Jozsa: Initial state: As in the Deutsch algorithm, reare the initial state 0 i def = 0i n i.

D. QUANTUM ALGORITHMS 9 Suerosition: Use the Walsh-Hadamard transformation to create a suerosition of all ossible inuts: i def =(H n H) 0 i = x, i. x n n Claim: Similarly to the single qubit case (Eq. III.), we can see that U f x, i =( ) f(x) xi i, where( ) n is an abbreviation for ( ) n. From the definition of i and U f, U f x, i = xi ( f(x)i f(x)i). Since f(x), ( f(x)i f(x)i) = i if f(x) =0,andit= i if f(x) =. This establishes the claim. Function alication: Therefore, you can see that: i def = U f i = x n n ( )f(x) xi i. (III.3) In the case of a constant function, all the comonents of the data state have the same hase, otherwise they do not. To see how we can make use of this information, let s consider the state in more detail. For a single bit you can show (Exer. III.4): H xi = ( 0i +( ) x i) = z( ) xz zi = z ( ) xz zi. (This is just another way of writing H 0i = ( 0i+ i) andh i = ( 0i i).) Therefore, the general formula for the Walsh transform of n bits is: H n x,x,...,x n i = = n n z,...,z n ( ) x z + +x nz n z,z,...,z n i z n ( ) x z zi, (III.4) where x z is the bitwise inner roduct. (It doesn t matter if you do addition or since only the arity of the result is significant.) Remember this formula! Combining this and the result in Eq. III.3, 3 i def =(H n I) i = ( n z n x n )x z+f(x) zi i.

30 CHAPTER III. QUANTUM COMPUTATION Measurement: Consider the first n qubits and the amlitude of one articular basis state, z = 0i = 0i n,whichwesearatefromtherestofthe summation: 3 i = x n n ( )f(x) 0i i + z n {0} x n n ( )x z+f(x) zi i Hence, the amlitude of 0i n,theall- 0i qubit string, is P x n n ( ) f(x). Recall how in the basic Deutsch algorithm we got a sum of signs (either all the same or not) for all the function evaluations. The result is similar here, but we have n values rather than just two. We now have two cases: Constant function: If the function is constant, then all the exonents of willbethesame(eitherall0orall),andsotheamlitudewillbe±. Therefore all the other amlitudes are 0 and any measurement must yield 0 for all the qubits (since only 0i n has nonzero amlitude). Balanced function: If the function is not constant then (ex hyothesi) it is balanced, but more secifically, if it is balanced, then there must be an equal number of + and contributionstotheamlitudeof 0i n,soits amlitude is 0. Therefore, when we measure the state, at least one qubit must be nonzero (since the all-0s state has amlitude 0). The good news about the Deutsch-Jozsa algorithm is that with one quantum function evaluation we have got a result that would require between and O( n )classicalfunctionevaluations(exonentialseedu!). Thebad news is that the algorithm has no known alications! Even if it were useful, however, the roblem could be solved robabilistically on a classical comuter with only a few evaluations of f: foranerrorrobabilityof, ittakes O(log ) function evaluations. However, it illustrates rinciles of quantum comuting that can be used in more useful algorithms.