Quadratic speedup for unstructured search - Grover s Al-

Similar documents
Grover s Algorithm + Quantum Zeno Effect + Vaidman

α x x 0 α x x f(x) α x x α x ( 1) f(x) x f(x) x f(x) α x = α x x 2

Chapter Newton s Method

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

IV. Performance Optimization

Solution 1 for USTC class Physics of Quantum Information

Lecture Notes on Linear Regression

Solution 1 for USTC class Physics of Quantum Information

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS

Singular Value Decomposition: Theory and Applications

14 The Postulates of Quantum mechanics

Lecture 10: May 6, 2013

Min Cut, Fast Cut, Polynomial Identities

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

1 GSW Iterative Techniques for y = Ax

5.04, Principles of Inorganic Chemistry II MIT Department of Chemistry Lecture 32: Vibrational Spectroscopy and the IR

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)

Markov Chain Monte Carlo (MCMC), Gibbs Sampling, Metropolis Algorithms, and Simulated Annealing Bioinformatics Course Supplement

C/CS/Phys C191 Grover s Quantum Search Algorithm 11/06/07 Fall 2007 Lecture 21

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017

18.1 Introduction and Recap

arxiv: v1 [quant-ph] 6 Sep 2007

Rigid body simulation

Which Separator? Spring 1

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

Calculation of time complexity (3%)

Design and Optimization of Fuzzy Controller for Inverse Pendulum System Using Genetic Algorithm

Four Bar Linkages in Two Dimensions. A link has fixed length and is joined to other links and also possibly to a fixed point.

Homework Assignment 3 Due in class, Thursday October 15

16 Reflection and transmission, TE mode

Newton s Method for One - Dimensional Optimization - Theory

Quantum Mechanics I - Session 4

Report on Image warping

Solutions to Homework 7, Mathematics 1. 1 x. (arccos x) (arccos x) 1

PHYS 215C: Quantum Mechanics (Spring 2017) Problem Set 3 Solutions

PES 1120 Spring 2014, Spendier Lecture 6/Page 1

6.842 Randomness and Computation February 18, Lecture 4

The non-negativity of probabilities and the collapse of state

This model contains two bonds per unit cell (one along the x-direction and the other along y). So we can rewrite the Hamiltonian as:

Case A. P k = Ni ( 2L i k 1 ) + (# big cells) 10d 2 P k.

CME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

Problem Set 9 Solutions

CHAPTER 4. Vector Spaces

Dynamic Programming. Preview. Dynamic Programming. Dynamic Programming. Dynamic Programming (Example: Fibonacci Sequence)

Lecture 20: November 7

Online Classification: Perceptron and Winnow

Section 8.3 Polar Form of Complex Numbers

Finding Dense Subgraphs in G(n, 1/2)

Topic 5: Non-Linear Regression

The Fourier Transform

ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00

Poisson brackets and canonical transformations

Finding Primitive Roots Pseudo-Deterministically

STATISTICAL MECHANICS

Solution Set #1

5 The Rational Canonical Form

Dynamic Programming 4/5/12. Dynamic programming. Fibonacci numbers. Fibonacci: a first attempt. David Kauchak cs302 Spring 2012

Mathematical Preparations

Markov Chain Monte Carlo Lecture 6

MMA and GCMMA two methods for nonlinear optimization

APPENDIX A Some Linear Algebra

Nice plotting of proteins II

A 2D Bounded Linear Program (H,c) 2D Linear Programming

find (x): given element x, return the canonical element of the set containing x;

Answers Problem Set 2 Chem 314A Williamsen Spring 2000

Feature Selection: Part 1

Hashing. Alexandra Stefan

Note 10. Modeling and Simulation of Dynamic Systems

EPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski

Structured quantum search in NP-complete problems using the cumulative density of states

Homework Notes Week 7

Stanford University CS254: Computational Complexity Notes 7 Luca Trevisan January 29, Notes for Lecture 7

Structure and Drive Paul A. Jensen Copyright July 20, 2003

A linear imaging system with white additive Gaussian noise on the observed data is modeled as follows:

Alternative Coins for Quantum Random Walk Search Optimized for a Hypercube

Review of Taylor Series. Read Section 1.2

CIS526: Machine Learning Lecture 3 (Sept 16, 2003) Linear Regression. Preparation help: Xiaoying Huang. x 1 θ 1 output... θ M x M

THE ARIMOTO-BLAHUT ALGORITHM FOR COMPUTATION OF CHANNEL CAPACITY. William A. Pearlman. References: S. Arimoto - IEEE Trans. Inform. Thy., Jan.

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Introduction to Algorithms

CSE 599d - Quantum Computing Introduction to Quantum Error Correction

Title: Radiative transitions and spectral broadening

2.3 Nilpotent endomorphisms

The Minimum Universal Cost Flow in an Infeasible Flow Network

THEOREMS OF QUANTUM MECHANICS

COS 521: Advanced Algorithms Game Theory and Linear Programming

Tracking with Kalman Filter

From Biot-Savart Law to Divergence of B (1)

The Geometry of Logit and Probit

Expected Value and Variance

Department of Electrical & Electronic Engineeing Imperial College London. E4.20 Digital IC Design. Median Filter Project Specification

Outline. Communication. Bellman Ford Algorithm. Bellman Ford Example. Bellman Ford Shortest Path [1]

Chapter Introduction

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011

10.34 Fall 2015 Metropolis Monte Carlo Algorithm

Appendix B: Resampling Algorithms

VQ widely used in coding speech, image, and video

IMGS-261 Solutions to Homework #9

Transcription:

Quadratc speedup for unstructured search - Grover s Al- CS 94- gorthm /8/07 Sprng 007 Lecture 11 001 Unstructured Search Here s the problem: You are gven a boolean functon f : {1,,} {0,1}, and are promsed that for eactly one a {1,, }, f(a) = 1 Thn of ths as a table of sze, where eactly one element has value 1, and all the others are 0 Snce we assume f can be computed classcally n polynomal tme, we can also compute t n superposton: α 0 α f() As we saw before, we can use crcut for f to put nformaton about f() n the phase by effectng the transformaton: α α ( 1) f() Here s another way of creatng ths phase state: ( α 0 1 ) ( f() f() ) α ( ) = α f() f() = α ( 1) f() ( 0 1 ) ow, we mght as well assume f s a blac bo or oracle All we need to do s desgn an algorthm that fnds a : f(a) = 1 00 Grover s algorthm Grover s algorthm fnds a n O( ) steps Consder the dmensonal Hlbert space spanned by 1,, We wsh to fnd a There s a state that we can create: = 1 Consder the two dmensonal subspace spanned by a and Let e be the state orthogonal to a n ths subspace Let θ be the angle between and e Then sn θ = 1/ and therefore θ 1/ See Fgure?? for an llustraton of these vectors a s the target, so we want to ncrease θ But how? One way to rotate a vector s to mae two reflectons In partcular, we can rotate a vector v by θ by reflectng about e and then reflectng about Ths transformaton s also llustrated n Fgure?? Each step of our algorthm s a rotaton by θ (we dscuss the mplementaton below) Ths means that we need π/ θ teratons for the algorthm to complete ow, what s θ? CS 94-, Sprng 007, Lecture 11 0-1

(apts) (v pts) (v pts) θ φ θ θ + φ (ψ 0 pts) = 1 =1 ( pts) (e pts) Fgure 01: To rotate v by θ to v, we reflect around e and then reflect around ψ 0 ψ 0 a = cos(π/ θ) = sn(θ) = 1 Snce sn θ θ, we now that θ 1 Thus, we need O( ) teratons for the algorthm to complete In the end, we get very close to a, and then wth hgh probablty, a measurement of the state yelds a How do you mplement the two reflectons? 1 Reflecton about e s easy We can reflect about the hyperplane orthogonal to a by flppng the phase of the component n the drecton of a ; e carry out the transformaton α α ( 1) f() For the reflecton about u, we wll actually reflect about n the dmensonal space as follows: apply the Hadamard transform H n to transform to 0 n ow apply a phase flp f the regster contents are anythng other than 0 n And apply the Hadamard transform to swtch bac from the Hadamard bass 003 Another approach Let s loo at the search algorthm dfferently, wth all superpostons The rotaton about u, D, s an nverson about the mean : (a) For = n, D can be decomposed and rewrtten as: CS 94-, Sprng 007, Lecture 11 0-

1 0 0 0 1 0 D = H H 0 0 1 0 0 = H + I H 0 0 = H H + I / / / / / / = + I / / / / + 1 / / / / + 1 / = / / / + 1 Observe that D s epressed as the product of three untary matrces (two Hadamard matrces separated by a condtonal phase shft matr) Therefore, D s also untary Regardng the mplementaton, both the Hadamard and the condtonal phase shft transforms can be effcently realzed wthn O(n) gates (b) Consder D operatng on a vector α to generate another vector β : D α 1 α = If we let be the mean ampltude, then the epresson descrbes a reflecton of about the mean Thus, the ampltude of β = j α j + = + can be consdered an nverson about the mean wth respect to The quantum search algorthm teratvely mproves the probablty of measurng a soluton Here s how: (a) Start state s ψ 0 = 1 (b) Invert the phase of a usng f (c) Then nvert about the mean usng D (d) Repeat steps and 3 O( ) tmes, so n each teraton α a ncreases by β 1 β β CS 94-, Sprng 007, Lecture 11 0-3

(a) (b) (c) Fgure 0: The frst three steps of Grover s algorthm We start wth a unform superposton of all bass vectors n (a) In (b), we have used the functon f to nvert the phase of α After runnng the dffuson operator D, we amplfy α whle decreasng all other ampltudes CS 94-, Sprng 007, Lecture 11 0-4

Ths process s llustrated n Fgure?? Suppose we just want to fnd a wth probablty 1 Untl ths pont, the rest of the bass vectors wll have 1 ampltude at least In each teraton of the algorthm, α a ncreases by at least = Eventually, α a = 1 The number of teratons to get to ths α a s 004 More applcatons Grover s algorthm s often called a database search algorthm Ths msnomer has been the cause of a lot of confuson, snce essental that the algorthm be able to query n superposton But there are a number of applcatons of unstructured search: (a) Fnd the mnmum n O( ) steps Eercse (b) Appromately count elements, or generate random ones (c) O( 1/3 ) algorthm for the collson problem (d) Speed up the test for matr multplcaton In ths problem we are gven three matrces, A, B, and C, and are told that the product of the frst two equals the thrd We wsh to verfy that ths s ndeed true An effcent (randomzed) way of dong ths s pcng a random array r, and checng to see whether Cr = ABr = A(Br) Classcally, we can do the chec n O(n ) tme, but usng a smlar approach to Grover s algorthm we can speed t up to O(n 175 ) tme CS 94-, Sprng 007, Lecture 11 0-5