Grover_4bits_trace.nb

$ket[x_String] := Apply[CircleTimes, Table[If[StringTake[x, {i, i}] == "0", ket[0], ket[1], "NotDef"], {i, StringLength[x]}]]$

$Clear[MatrixDisplay] MatrixDisplay[tensor_CircleTimes] := Map[MatrixForm, tensor] MatrixDispla ... matrix_] := Map[MatrixForm, matrix]/;MatrixQ[matrix] MatrixDisplay[ket__] := MatrixDisplay[{ket}]$

$bitQ[x_] := If[VectorQ[x] && Dimensions[x] == {2}, True, False, False]$

$Clear[Tensor] Tensor[tensorList_CircleTimes] := Apply[Tensor, tensorList] Tensor[{ket_}] := ke ... Table[fm[[i, j]] tensor2, {j, Length[fm[[i]]]}]], {i, Length[fm]}]]] /;VectorQ[ketList, MatrixQ]$

Clear[WalshHadamard] WalshHadamard[ket_] := WalshHadamardMatrix . ket /;bitQ[ket] WalshHadamar ... List] /; VectorQ[ketList, bitQ] WalshHadamard[ketList_CircleTimes] := Map[WalshHadamard, ketList]

$Clear[GroverSearchMain] GroverSearchMain[f_, BasicState_List, iteration_Integer] := Module[{Rf ... {n}]]] ; amplitude[k_Integer] := amplitude[k] = Uf . amplitude[k - 1] ; amplitude[iteration] ]$

$Clear[GroverSearchMainTrace] GroverSearchMainTrace[f_, BasicState_List, iteration_Integer] := ... e[k_Integer] := amplitude[k] = Uf . amplitude[k - 1] ; Table[amplitude[i], {i, 0, iteration}] ]$

$<<Graphics`Graphics` Clear[GroverShowGraphics] GroverShowGraphics[amplitude_List] := Mo ... esOrigin -> {0, 0}] ; Do[ShowGraphics[i], {i, Dimensions[amplitude][[1]]}] ]/;MatrixQ[amplitude]$

$GroverSearchSolve[amplitude_List] := Module[{r, i, s}, r = Random[] ; s = amplitude[[1]]^2 ; ... i - 1]]^2, "=", N[amplitude[[i - 1]]^2, 3], "で観測されました"] ]/;VectorQ[amplitude]$

$GroverSearchSolveRetNum[amplitude_List] := Module[{r, i, s}, r = Random[] ; s = amplitude[[1]]^2 ; For[i = 2, s< r, i ++, s = s + amplitude[[i]]^2] ; i - 2 ]/;VectorQ[amplitude]$

$Clear[GroverSearch] GroverSearch[f_, Number_Integer, t_Integer] := Module[{n, stateX, i}, n = ... ket[i, n], {i, 0, Num - 1}] ; i = Floor[N[π]/4 Num/t^(1/2)] ; GroverSearchMain[f, stateX, i] ]$

$Clear[GroverSearchTrace] GroverSearchTrace[f_, Number_Integer, t_Integer] := Module[{n, stateX ... eiling[i] - i) >1/2, i = Floor[i], i = Ceiling[i]] ; *) GroverSearchMainTrace[f, stateX, i] ]$

[Graphics:HTMLFiles/index_19.gif]

[Graphics:HTMLFiles/index_20.gif]

[Graphics:HTMLFiles/index_21.gif]

[Graphics:HTMLFiles/index_22.gif]

Created by Mathematica (February 3, 2005)