CN110162536B - 一种量子搜索方法、系统、电子装置及存储介质 - Google Patents

一种量子搜索方法、系统、电子装置及存储介质 Download PDF

Info

Publication number
CN110162536B
CN110162536B CN201910286148.9A CN201910286148A CN110162536B CN 110162536 B CN110162536 B CN 110162536B CN 201910286148 A CN201910286148 A CN 201910286148A CN 110162536 B CN110162536 B CN 110162536B
Authority
CN
China
Prior art keywords
state
computer system
quantum
measurement result
module
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910286148.9A
Other languages
English (en)
Other versions
CN110162536A (zh
Inventor
王平
马利芬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen University
Original Assignee
Shenzhen University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenzhen University filed Critical Shenzhen University
Priority to CN201910286148.9A priority Critical patent/CN110162536B/zh
Publication of CN110162536A publication Critical patent/CN110162536A/zh
Application granted granted Critical
Publication of CN110162536B publication Critical patent/CN110162536B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/24Querying
    • G06F16/242Query formulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N10/00Quantum computing, i.e. information processing based on quantum-mechanical phenomena

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Artificial Intelligence (AREA)
  • Computational Linguistics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computing Systems (AREA)
  • Software Systems (AREA)
  • Superconductor Devices And Manufacturing Methods Thereof (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开了一种量子搜索方法、系统、电子装置及存储介质,用于量子计算,解决了现有技术中的量子搜索方法有待提高的问题,其包括:在接收需要搜索的问题后,初始化用于量子计算的计算机系统,得到计算机系统的初始状态;对计算机系统的初始状态进行转换,得到计算机系统的均衡叠加态;使用G算子迭代并更新计算机系统的均衡叠加态,得到计算机系统的系统状态;对计算机系统的系统状态直接进行测量,得到测量结果;使用经典Oracle函数对所述测量结果进行验证,在验证所述测量结果为问题正确解的目标状态后,输出所述测量结果;从而降低了量子搜索方法线路的复杂性,提高了量子搜索方法的速率。

Description

一种量子搜索方法、系统、电子装置及存储介质
技术领域
本发明涉及信息安全技术领域,尤其涉及一种量子搜索方法、系统、电子装置及存储介质。
背景技术
由于量子特性在信息处理过程中有着特殊的性质,因而可能在提高算法速率、提供安全性、增大信息存储量和提高精准度等这几个方面突破经典计算的极限。因此,将量子力学首次应用到计算机科学领域中,诞生一门新的交叉学科,即量子计算,对国家安全、军事和产业等领域有着颠覆性影响的研究方向。
在计算机科学中,从无序的数据库中找到目标数据是数据库搜索问题。该问题为:假设有N个无序条目,将其标号,确定里面是否存在一个‘ABC’条目,如果存在这个‘ABC’条目,找出这个条目的下标。经典算法需经过O(N)次运算才能找到这个条目。1996年,Grover提出了运算速度得到平方增长的量子搜索方法,即Grover算法,Grover算法通过反复使用G算子,放大目标状态的概率幅,减小非目标状态的概率幅,以接近1的概率找到目标状态。该算法将oracle调用次数从经典算法的O(N)降低到
Figure BDA0002023327770000011
但是,虽然Grover算法能以高概率测量到搜索问题解的状态,但是应用了更多次的G量子线路,使整个算法的量子线路相较于经典算法更为复杂,因此现有的量子搜索方法,运算效率仍有待提高。
发明内容
本发明的主要目的在于提供一种量子搜索方法,旨在解决现有技术中的量子搜索方法有待提高的技术问题。
为实现上述目的,本发明第一方面提供一种量子搜索方法,包括:在接收需要搜索的问题后,初始化用于量子计算的计算机系统,得到计算机系统的初始状态;对计算机系统的初始状态进行转换,得到计算机系统的均衡叠加态;使用G算子迭代并更新计算机系统的均衡叠加态,得到计算机系统的系统状态;对计算机系统的系统状态直接进行测量,得到测量结果;使用经典Oracle函数对所述测量结果进行验证,在验证所述测量结果为问题正确解的目标状态后,输出所述测量结果。
进一步地,所述对计算机系统的初始状态进行转换,得到计算机系统的均衡叠加态包括:将计算机系统的初始状态经过Hadamard变换,得到计算机系统的第一均衡叠加态;将计算机系统归一化的目标状态及归一化的非目标状态作为二维空间中的状态向量,并在计算机系统的第一均衡叠加态的基础上,重新构造计算机系统的系统状态,得到适用于G算子的计算机系统的系统状态。
进一步地,所述使用G算子迭代并更新计算机系统的均衡叠加态,得到计算机系统的系统状态包括:使用Oracle算子查找计算机系统的第一均衡叠加态中是否存在问题解的状态,若存在所述问题解的状态,则标记所述问题解的状态;对经过Oracle算子查找之后的计算机系统的系统状态进行条件相移操作,以更新计算机系统的均衡叠加态;使用G算子对更新后的计算机系统的均衡叠加态进行
Figure BDA0002023327770000021
次迭代并在每次迭代后更新计算机系统的系统状态。
进一步地,所述使用经典Oracle函数对所述测量结果进行验证,在验证所述测量结果为问题正确解的目标状态后,输出所述测量结果包括:使用经典Oracle函数对所述测量结果进行验证;若所述测量结果经过经典Oracle函数验证后与问题正确解的目标状态一致,则输出所述测量结果作为问题正确解;若所述测量结果经过经典Oracle函数验证后与问题正确解的目标状态不一致,则重新初始化用于量子计算的计算机系统,得到计算机系统的初始状态,并重复后续步骤,直至所述测量结果与问题正确解的目标状态一致,输出所述测量结果。
本发明第二方面提供一种量子搜索系统,包括:量子态制备模块,用于初始化用于量子计算的计算机系统,得到计算机系统的初始状态;量子态变换模块,用于对所述量子态制备模块得到的计算机系统的初始状态进行转换,得到计算机系统的均衡叠加态;量子态更新模块,用于使用G算子迭代并更新所述量子态变换模块得到的计算机系统的均衡叠加态,得到计算机系统的系统状态;测量模块,用于对经过所述量子态更新模块更新后的计算机系统的系统状态直接进行测量,得到测量结果;结果验证模块,用于使用经典Oracle函数对所述测量模块得到的测量结果进行验证,在验证所述测量结果为问题正确解的目标状态后,输出所述测量结果。
进一步地,所述量子态变换模块包括:第一变换单元,用于将所述量子态制备模块得到的计算机系统的初始状态经过Hadamard变换,得到计算机系统的第一均衡叠加态;第二变换单元,用于将目标状态划分单元划分的计算机系统归一化的目标状态及归一化的非目标状态作为二维空间中的状态向量,并在所述第一变换单元得到的计算机系统的第一均衡叠加态的基础上,重新构造计算机系统的系统状态,得到适用于G算子的计算机系统的系统状态。
进一步地,所述量子态更新模块包括:状态标记单元,用于在Oracle中查找量子的均衡叠加态中是否存在搜索问题解的状态,若存在所述问题解的状态,则标记所述问题解的状态;相位转移单元,用于对经过所述状态标记单元的Oracle算子查找之后的计算机系统的系统状态进行条件相移操作,以更新计算机系统的均衡叠加态;迭代更新单元,用于使用G算子对所述相位转移单元更新后的计算机系统的均衡叠加态进行
Figure BDA0002023327770000031
次迭代,并在每次迭代后更新计算机系统的系统状态。
进一步地,所述结果验证模块包括:Oracle验证单元,用于使用经典Oracle函数对所述测量模块得到的测量结果进行验证;测量结果输出单元,用于在所述测量模块得到的测量结果经过经典Oracle函数验证后与所述量子态变换模块构造的问题正确解的目标状态一致,则输出所述测量结果作为问题正确解;循环单元,用于若所述测量模块得到的测量结果经过所述Oracle验证单元的验证后与所述量子态变换模块构造的问题正确解的目标状态不一致,则重新调用量子态制备模块、量子态变换模块、量子态更新模块、测量模块及结果验证模块,直至所述测量结果与问题正确解的目标状态一致,输出所述测量结果。
本发明第三方面提供一种电子装置,包括:存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时,实现上述中的任意一项所述方法。
本发明第四方面提供一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时,实现上述中的任意一项所述方法。
本发明提供一种量子搜索方法、系统、电子装置及存储介质,有益效果在于:在计算计算机系统的均衡叠加态并更新计算机系统的均衡叠加态后,通过多次测量及验证的方法来确定问题的正确解,减少了G量子线路的调用,从而显著降低了量子搜索方法线路的复杂性,提高了量子搜索方法的运算效率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例量子搜索方法的流程示意框图;
图2为本发明实施例量子搜索方法的量子线路简图;
图3为本发明实施例量子搜索方法测量的量子线路和经典Oracle函数线路;
图4为本发明实施例量子搜索方法的G算子迭代的量子线路简图;
图5为本发明实施例量子搜索方法的几何过程描述图;
图6为本发明实施例量子搜索方法经过x次G算子迭代的几何描述图;
图7为本发明实施例量子搜索系统的结构示意框图;
图8为本发明电子装置的结构示意框图。
具体实施方式
为使得本发明的发明目的、特征、优点能够更加的明显和易懂,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而非全部实施例。基于本发明中的实施例,本领域技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参阅图1,为一种量子搜索方法,包括:S1、在接收需要搜索的问题后,初始化用于量子计算的计算机系统,得到计算机系统的初始状态;S2、对计算机系统的初始状态进行转换,得到计算机系统的均衡叠加态;S3、使用G算子迭代并更新计算机系统的均衡叠加态,得到计算机系统的系统状态;S4、对计算机系统的系统状态直接进行测量,得到测量结果;S5、使用经典Oracle函数对测量结果进行验证,在验证测量结果为问题正确解的目标状态后,输出测量结果。
对计算机系统的初始状态进行转换,得到计算机系统的均衡叠加态包括:将计算机系统的初始状态经过Hadamard变换,得到计算机系统的第一均衡叠加态;将计算机系统归一化的目标状态及归一化的非目标状态作为二维空间中的状态向量,并在计算机系统的第一均衡叠加态的基础上,系统状态,得到适用于G算子的计算机系统的系统状态。
使用G算子迭代并更新计算机系统的均衡叠加态,得到计算机系统的系统状态包括:使用Oracle算子查找计算机系统的第一均衡叠加态中是否存在问题解的状态,若存在问题解的状态,则标记问题解的状态;对经过Oracle算子查找之后的计算机系统的系统状态进行条件相移操作,以更新计算机系统的均衡叠加态;使用G算子对更新后的计算机系统的均衡叠加态进行
Figure BDA0002023327770000051
次迭代并在每次迭代后更新计算机系统的系统状态。
使用经典Oracle函数对测量结果进行验证,在验证测量结果为问题正确解的目标状态后,输出测量结果包括:使用经典Oracle函数对测量结果进行验证;若测量结果经过经典Oracle函数验证后与问题正确解的目标状态一致,则输出测量结果作为问题正确解;若测量结果经过经典Oracle函数验证后与问题正确解的目标状态不一致,则重新初始化用于量子计算的计算机系统,得到计算机系统的初始状态,并重复后续步骤,直至测量结果与问题正确解的目标状态一致,输出测量结果。
请参阅图2,本申请提供的量子搜索方法,使用了两个寄存器,第一个寄存器使用了n个量子比特;第二个寄存器使用少量量子比特作为Oracle的工作空间。
请参阅图3,在图3中(a)为测量的量子线路,(b)为经典Oracle函数线路,因此本申请提供的量子搜索方法,先制备计算机系统的初始状态,
Figure BDA0002023327770000061
其次应用Hadamard变换,把初始状态
Figure BDA0002023327770000066
变成均衡叠加态,即如公式1所示,公式1表示如下:
Figure BDA0002023327770000062
随后使用G算子迭代,由图2可知,本申请提供的量子搜索方法应用
Figure BDA0002023327770000063
次G算子进行迭代,采用一次测量和经典Oracle函数验证能以1/2的概率得到搜索问题的解,因此,经过n次测量和验证,成功找到搜索问题解的概率为1一(1/2)n。完成G算子迭代的量子线路如图4所示,可分为如下步骤:
(1)应用oracle操作O,查找叠加状态|ψ>中是否存在搜索问题解的状态。如果存在这个状态,就标记这个状态。
(2)对(1)中的结果使用Hadamard变换
Figure BDA0002023327770000064
(3)实施条件相移操作,使状态|0>之外的其他状态获得-1的相位,即如公式2所示,公式2表示如下:
Figure BDA0002023327770000067
(4)对(3)中的状态,实施Hadamard变换
Figure BDA0002023327770000065
由此可见,实施一次G算子的迭代操作要经过(1)-(4)的四个步骤,其中步骤(2)至步骤(4)总的作用如公式3所示,公式3表示如下:
Figure BDA0002023327770000071
公式3中的状态|ψ>为计算机系统的均衡叠加状态。所以G算子迭代操作可以使用公式4表示,公式4表示如下:
G=(2|ψ><ψ|-I)O
在公式4中,O是酉矩阵,表述为O=I-2|τ><τ|;实施公式4中的迭代操作,输出状态如公式5所示,公式5表示如下:
Figure BDA0002023327770000072
对公式5中的状态再次实施G算子的迭代操作,将公式5的状态进行变换,变换后的状态如公式6所示,公式6表示如下:
Figure BDA0002023327770000073
对公式6表示的状态再次进行G算子的迭代操作,将公式6的状态进行变换,变换后的状态如公式7所示,公式7表示如下:
Figure BDA0002023327770000074
使用G算子进行三次迭代后可发现,均衡叠加状态|ψ>概率幅随着迭代次数的增加而减少,搜索问题的解状态|τ>的概率幅随着迭代次数的增加而增加。按照这样的规律,初始状态经过x次G算子迭代之后,计算机系统的状态更新为|ψx>=Gx|ψ)。
使用G算子迭代并更新计算机系统的均衡叠加态包括:使用Oracle算子查找计算机系统的第一均衡叠加态中是否存在问题解的状态,若存在问题解的状态,则标记问题解的状态;对标记的问题解的状态进行Hadamard变换,得到问题解的均衡叠加态;对问题解的均衡叠加态进行条件相移,使得相位为0之外的其他状态的问题解的均衡叠加态相位转变为-1;对相位为O及相位为-1的问题解的均衡叠加态进行Hadamard变换,以更新计算机系统的均衡叠加态;使用G算子对更新后的计算机系统的均衡叠加态进行
Figure BDA0002023327770000081
次迭代并在每次迭代后更新计算机系统的均衡叠加态。
由公式4可知量子搜索方法的迭代算子G=(2|ψ><ψ|-I)O;其中,算子2|ψ><ψ|-I和O可以看作是在二维空间中关于对应量子态的两次反射变换。因此该概率型量子搜索算法的迭代算子G可以看成是由搜索问题解的状态和非搜索问题解成的二维空间中的一个旋转算子。该概率型量子搜索算法过程的几何思想是:应用
Figure BDA0002023327770000082
次的旋转算子G,使均衡叠加状态|ψ>以1/2的概率接近目标状态|β>。这里描述的目标状态即为搜索问题解的状态。因此该概率型量子搜索算法的几何过程如下:
首先将公式1表示的计算机系统的均衡叠加态分为归一化的目标状态和归一化的非目标状态,归一化的目标状态及归一化的非目标状态分别如公式8及公式9所示,公式8表示如下:
Figure BDA0002023327770000083
公式9表示如下:
Figure BDA0002023327770000084
在公式8及公式9中,N是所有状态的总数,M是搜索问题的解的个数,因此,公式1中的计算机系统的初始状态|ψ>可以重新定义,重新定义的|ψ>如公式10所示,公式10表示如下:
Figure BDA0002023327770000085
在公式10中,计算机系统的状态|ψ>看作是由状态向量|α>和向量|β>所形成的二维空间里的向量。公式10中|ψ>分量的系数如公式11所示,公式11表示如下:
Figure BDA0002023327770000091
Figure BDA0002023327770000092
则对公式10中的|ψ>再次构造为公式12,公式12表示如下:
Figure BDA0002023327770000093
如图5所示,公式12中的|ψ>为二维平面中的初始状态。迭代算子G是图5中的一个旋转算子,即G=(2|ψ><ψ|-I)O。其中,O算子是通过改变目标状态的相位,标记了这个目标状态的酉算子。且这个算子为|α>和|β>所张成的二维平面上相对于向量|α>实施了一次反射变换,即O(a|α>+b|β>)=a|α>-b|β>。类似地,在|α>和|β>所形成的平面上,2|ψ><ψ|-I算子是相对于向量|ψ>进行了一次反射变换。上述两次反射变换的乘积形成了一次旋转操作,即为旋转算子G。
本申请提供的量子搜索方法过程是:反复应用
Figure BDA0002023327770000094
次G算子,采用多次测量和经典Oracle函数的程序验证技术,能以大于1/2的概率找到搜索问题的解。通过一次测量和经典oracle函数验证能以1/2的概率找到搜索问题的解。如若经过n次测量和验证,成功找到搜索问题解的概率为1-(1/2)n
如图5所示,计算机系统的均衡叠加状态|ψ>经过一次旋转算子G作用后,公式12中的|ψ>发生变换,得到|ψ1>,|ψ1>如公式13所示,公式13表示如下:
Figure BDA0002023327770000095
第二次实施旋转算子G后,公式13中的|ψ1>发生变化,|ψ1>转变为|ψ2>,|ψ2>如公式14表示,公式14表示如下:
Figure BDA0002023327770000096
第三次实施旋转算子G后,公式14中的|ψ2>转变为|ψ3>,|ψ3>如公式15表示,公式15表示如下:
Figure BDA0002023327770000097
在对|ψ>进行了三次迭代时,|ψ>每经过一次旋转算子G后,状态相位增加一个旋转角度θ;如图6所示,按照迭代的规律,在x次旋转算子G迭代操作之后,|ψ>状态变换为|ψx>,|ψx>如公式16所示,公式16表示如下:
Figure BDA0002023327770000101
若公式16中,
Figure BDA0002023327770000102
对|ψx>进行测量能以1/2的概率测量到搜索问题的解状态|β>。在图中,此时量子状态|ψx>与向量|α>之间的夹角为π/4。
就量子搜索方法的量子线路而言,在原始Grover算法运行过程中,G算子线路是由单量子比特门的复合而成的量子线路。G算子的量子线路比测量的量子线路和经典Oracle函数验证线路更复杂。在原始Grover算法中,均衡叠加状态经过
Figure BDA0002023327770000103
次G迭代之后,对输出状态进行一次测量能以1的概率测量到搜索问题解的状态。该算法虽然能以高概率测量到搜索问题解的状态,但是应用了更多次的G量子线路,使整个算法的量子线路更为复杂。相比于G量子线路的复杂度,测量和验证技术的复杂度要小很多。因此,在概率型量子搜索算法中,实施少于原始Grover算法的G算子迭代,应用多次测量和经典Oracle函数程序验证技术降低G量子线路的调用次数,以此降低整个算法过程的量子线路的复杂度和降低算法调用Oracle的次数。该改进算法经过
Figure BDA0002023327770000104
次G算子迭代之后,对输出状态进行一次测量和经典Oracle函数验证能以1/2的概率找到搜索问题的解。所以该算法是具有一定成功率的算法。因此经过n次测量和验证,成功找到搜索问题解的概率为1-(1/2)n。如果测量结果验证成功,则测量结果就是搜索问题的解,算法成功。否则算法失败,重新进行该算法运算,直至验证的测量结果是搜索问题的解为止。相比于原始Grover算法,Grover改进算法以多次测量和验证技术为代价,减少了
Figure BDA0002023327770000105
次G量子线路的调用,显著降低了计算复杂性。
根据算法几何过程中旋转角度的变换规律,应用一次G算子,系统状态的相位会增加θ角度。随着迭代次数的增加,叠加状态|ψ>的相位会相应的增加。若经过x次G迭代之后,对|ψx>进行测量能以1/2的概率测量到目标状态|β>;则公式16中目标状态的系数满足公式17所示的条件,公式17表示如下:
Figure BDA0002023327770000111
通过公式17可得,本申请提供的量子搜索方法找到搜索问题的解的最少迭代次数为公式18所示,公式18表示如下:
Figure BDA0002023327770000112
在公式18中,CI(x)表示取实数x最接近的整数;由公式11可知
Figure BDA0002023327770000113
当M<<N,则
Figure BDA0002023327770000114
因此,本申请提供的量子搜索方法的最少迭代次数如公式19所示,公式19表示如下:
Figure BDA0002023327770000115
当M=1时,
Figure BDA0002023327770000116
本申请提供的量子搜索方法的迭代次数如公式20所示,公式20表示如下:
Figure BDA0002023327770000117
而在经过
Figure BDA0002023327770000118
次G算子迭代之后,通过测量和采用经典Oracle函数验证,就能以1/2的概率找到搜索问题的解。如果测量结果验证成功,则测量结果就是搜索问题的解,算法成功。否则算法失败,重新进行该算法运算,直至验证的测量结果是搜索问题的解为止。由此,Grover改进算法是最多调
Figure BDA0002023327770000119
次Oracle算子,然后经过多次测量和经典Oracle函数程序验证技术,能以大于1/2的概率找到搜索问题的解的改进算法。因为该算法是具有一定成功率的算法,所以经过n次测量和验证之后,成功找到搜索问题解的概率为1-(1/2)n
若假设经过y次G算子迭代,本申请提供的量子搜索方法通过多次测量和经典Oracle验证能以概率P找到搜索问题的解。由公式16可知,经过y次G迭代之后的计算机系统的状态|ψy>可以如公式21所示,公式21表示如下:
Figure BDA0002023327770000121
在满足公式22的条件下,本申请提供的量子搜索方法的解的最少迭代次数可以如公式23所示;
公式22表示如下:
Figure BDA0002023327770000122
公式23表示如下:
Figure BDA0002023327770000123
CI(y)表示取实数y最接近的整数,由公式11可知
Figure BDA0002023327770000124
当M<<N,则
Figure BDA0002023327770000125
令M=1,则
Figure BDA0002023327770000126
则本申请提供的量子搜索方法的G算子迭代次数如公式24所示,公式24表示如下:
Figure BDA0002023327770000127
从而可知,本申请提供的量子搜索方法经过
Figure BDA0002023327770000128
次G算子迭代之后,能以P的概率测量到搜索问题的解状态,然后采用经典Oracle函数验证,就能以P的概率找到搜索问题的解。如果测量结果验证成功,则测量结果就是搜索问题的解,算法成功。否则算法失败,重新进行该方法进行运算,直至验证的测量结果是搜索问题的解为止。因此,量子搜索方法最多调
Figure BDA0002023327770000129
次Oracle算子,然后经过多次测量和经典Oracle函数验证的方法,以大于P的概率找到搜索问题的解的改进算法。因为该算法是具有一定成功率的方法,所以经过n次测量和验证之后,成功找到搜索问题解的概率为1-(1-P)n。由上述分析可知,若以P的概率找到搜索问题的解,则申请提供的量子搜索方法的计算复杂度为
Figure BDA0002023327770000131
相比于原始Grover算法,申请提供的量子搜索方法的G算子迭代次数从
Figure BDA0002023327770000132
降低到
Figure BDA0002023327770000133
同时整个算法的量子线路减少了
Figure BDA0002023327770000134
Figure BDA0002023327770000135
个G算子的量子线路。而通过分析可得,当户=1时,本申请提供的量子搜索方法的计算复杂度为
Figure BDA0002023327770000136
此种情况下,是本申请提供的量子搜索方法的一个特殊情况,即为原始Grover算法。当P=1/2,本申请提供的量子搜索方法的计算复杂度
Figure BDA0002023327770000137
本申请实施例提供一种量子搜索系统,请参阅图7,包括:量子态制备模块1、量子态变换模块2、量子态更新模块3、测量模块4及结果验证模块5;量子态制备模块1用于初始化用于量子计算的计算机系统,得到计算机系统的初始状态;量子态变换模块2用于对量子态制备模块1得到的计算机系统的初始状态进行转换,得到计算机系统的均衡叠加态;量子态更新模块3用于使用G算子迭代并更新量子态变换模块2得到的计算机系统的均衡叠加态,得到计算机系统的系统状态;测量模块4用于对经过量子态更新模块3更新后的计算机系统的系统状态直接进行测量,得到测量结果;结果验证模块5用于使用经典Oracle函数对测量模块4得到的测量结果进行验证,在验证测量结果为问题正确解的目标状态后,输出测量结果。
量子态变换模块2包括:第一变换单元及第二变换单元;第一变换单元用于将量子态制备模块1得到的计算机系统的初始状态经过Hadamard变换,得到计算机系统的第一均衡叠加态;第二变换单元用于将目标状态划分单元划分的计算机系统归一化的目标状态及归一化的非目标状态作为二维空间中的状态向量,并在第一变换单元得到的计算机系统的第一均衡叠加态的基础上,重新构造计算机系统的系统状态,得到适用于G算子的计算机系统的系统状态。
量子态更新模块3包括:状态标记单元、相位转移单元及迭代更新单元;状态标记单元用于在Oracle中查找量子的均衡叠加态中是否存在搜索问题解的状态,若存在问题解的状态,则标记问题解的状态;相位转移单元,用于对经过状态标记单元的Oracle算子查找之后的计算机系统的系统状态进行条件相移操作,以更新计算机系统的均衡叠加态;迭代更新单元用于使用G算子对相位转移单元更新后的计算机系统的均衡叠加态进行
Figure BDA0002023327770000141
次迭代,并在每次迭代后更新计算机系统的系统状态。
结果验证模块5包括:Oracle验证单元、测量结果输出单元及循环单元;Oracle验证单元用于使用经典Oracle函数对测量模块4得到的测量结果进行验证;测量结果输出单元用于在测量模块4得到的测量结果经过经典Oracle函数验证后与量子态变换模块2构造的问题正确解的目标状态一致,则输出测量结果作为问题正确解;循环单元用于若测量模块4得到的测量结果经过Oracle验证单元的验证后与量子态变换模块2构造的问题正确解的目标状态不一致,则重新调用量子态制备模块1、量子态变换模块2、量子态更新模块、测量模块4及结果验证模块5,直至测量结果与问题正确解的目标状态一致,输出测量结果。
本申请实施例提供一种电子装置,请参阅8,该电子装置包括:存储器601、处理器602及存储在存储器601上并可在处理器602上运行的计算机程序,处理器602执行该计算机程序时,实现前述中描述的量子搜索方法。
进一步的,该电子装置还包括:至少一个输入设备603以及至少一个输出设备604。
上述存储器601、处理器602、输入设备603以及输出设备604,通过总线605连接。
其中,输入设备603具体可为摄像头、触控面板、物理按键或者鼠标等等。输出设备604具体可为显示屏。
存储器601可以是高速随机存取记忆体(RAM,Random Access Memory)存储器,也可为非不稳定的存储器(non-volatile memory),例如磁盘存储器。存储器601用于存储一组可执行程序代码,处理器602与存储器601耦合。
进一步的,本申请实施例还提供了一种计算机可读存储介质,该计算机可读存储介质可以是设置于上述各实施例中的电子装置中,该计算机可读存储介质可以是前述中的存储器601。该计算机可读存储介质上存储有计算机程序,该程序被处理器602执行时实现前述实施例中描述的量子搜索方法。
进一步的,该计算机可存储介质还可以是U盘、移动硬盘、只读存储器601(ROM,Read-Only Memory)、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
在本申请所提供的几个实施例中,应该理解到,所揭露的装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,所述模块的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个模块或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或模块的间接耦合或通信连接,可以是电性,机械或其它的形式。
所述作为分离部件说明的模块可以是或者也可以不是物理上分开的,作为模块显示的部件可以是或者也可以不是物理模块,即可以位于一个地方,或者也可以分布到多个网络模块上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能模块可以集成在一个处理模块中,也可以是各个模块单独物理存在,也可以两个或两个以上模块集成在一个模块中。上述集成的模块既可以采用硬件的形式实现,也可以采用软件功能模块的形式实现。
所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。需要说明的是,对于前述的各方法实施例,为了简便描述,故将其都表述为一系列的动作组合,但是本领域技术人员应该知悉,本发明并不受所描述的动作顺序的限制,因为依据本发明,某些步骤可以采用其它顺序或者同时进行。其次,本领域技术人员也应该知悉,说明书中所描述的实施例均属于优选实施例,所涉及的动作和模块并不一定都是本发明所必须的。
在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述的部分,可以参见其它实施例的相关描述。
以上为对本发明所提供的一种量子搜索方法、系统、电子装置及存储介质的描述,对于本领域的技术人员,依据本发明实施例的思想,在具体实施方式及应用范围上均会有改变之处,综上,本说明书内容不应理解为对本发明的限制。

Claims (8)

1.一种量子搜索方法,其特征在于,包括:
在接收需要搜索的问题后,初始化用于量子计算的计算机系统,得到计算机系统的初始状态;
对计算机系统的初始状态进行转换,得到计算机系统的均衡叠加态;
使用G算子迭代并更新计算机系统的均衡叠加态,得到计算机系统的系统状态,其中包括:
将计算机系统的初始状态经过Hadamard变换,得到计算机系统的第一均衡叠加态,
将计算机系统归一化的目标状态及归一化的非目标状态作为二维空间中的状态向量,并在计算机系统的第一均衡叠加态的基础上,重新构造计算机系统的系统状态,得到适用于G算子的计算机系统的系统状态;
对计算机系统的系统状态直接进行测量,得到测量结果;
使用经典Oracle函数对所述测量结果进行验证,在验证所述测量结果为问题正确解的目标状态后,输出所述测量结果。
2.根据权利要求1所述的量子搜索方法,其特征在于,
所述使用G算子迭代并更新计算机系统的均衡叠加态,得到计算机系统的系统状态包括:
使用Oracle算子查找计算机系统的第一均衡叠加态中是否存在问题解的状态,若存在所述问题解的状态,则标记所述问题解的状态;
对经过Oracle算子查找之后的计算机系统的系统状态进行条件相移操作,以更新计算机系统的均衡叠加态;
使用G算子对更新后的计算机系统的均衡叠加态进行
Figure FDA0002986869260000011
次迭代并在每次迭代后更新计算机系统的系统状态,其中,O是算子通过改变目标状态的相位目标状态的酉算子,N是所有状态的总数。
3.根据权利要求1所述的量子搜索方法,其特征在于,
所述使用经典Oracle函数对所述测量结果进行验证,在验证所述测量结果为问题正确解的目标状态后,输出所述测量结果包括:
使用经典Oracle函数对所述测量结果进行验证;
若所述测量结果经过经典Oracle函数验证后与问题正确解的目标状态一致,则输出所述测量结果作为问题正确解;
若所述测量结果经过经典Oracle函数验证后与问题正确解的目标状态不一致,则重新初始化用于量子计算的计算机系统,得到计算机系统的初始状态,并重复后续步骤,直至所述测量结果与问题正确解的目标状态一致,输出所述测量结果。
4.一种量子搜索系统,其特征在于,包括:
量子态制备模块,用于初始化用于量子计算的计算机系统,得到计算机系统的初始状态;
量子态变换模块,用于对所述量子态制备模块得到的计算机系统的初始状态进行转换,得到计算机系统的均衡叠加态,其中包括:
将计算机系统的初始状态经过Hadamard变换,得到计算机系统的第一均衡叠加态,
将计算机系统归一化的目标状态及归一化的非目标状态作为二维空间中的状态向量,并在计算机系统的第一均衡叠加态的基础上,重新构造计算机系统的系统状态,得到适用于G算子的计算机系统的系统状态;
量子态更新模块,用于使用G算子迭代并更新所述量子态变换模块得到的计算机系统的均衡叠加态,得到计算机系统的系统状态;
测量模块,用于对经过所述量子态更新模块更新后的计算机系统的系统状态直接进行测量,得到测量结果;
结果验证模块,用于使用经典Oracle函数对所述测量模块得到的测量结果进行验证,在验证所述测量结果为问题正确解的目标状态后,输出所述测量结果。
5.根据权利要求4所述的量子搜索系统,其特征在于,
所述量子态更新模块包括:
状态标记单元,用于在Oracle中查找量子的均衡叠加态中是否存在搜索问题解的状态,若存在所述问题解的状态,则标记所述问题解的状态;
相位转移单元,用于对经过所述状态标记单元的Oracle算子查找之后的计算机系统的系统状态进行条件相移操作,以更新计算机系统的均衡叠加态;
迭代更新单元,用于使用G算子对所述相位转移单元更新后的计算机系统的均衡叠加态进行
Figure FDA0002986869260000031
次迭代,并在每次迭代后更新计算机系统的系统状态,其中,O是算子通过改变目标状态的相位目标状态的酉算子,N是所有状态的总数。
6.根据权利要求4所述的量子搜索系统,其特征在于,
所述结果验证模块包括:
Oracle验证单元,用于使用经典Oracle函数对所述测量模块得到的测量结果进行验证;
测量结果输出单元,用于在所述测量模块得到的测量结果经过经典Oracle函数验证后与所述量子态变换模块构造的问题正确解的目标状态一致,则输出所述测量结果作为问题正确解;
循环单元,用于若所述测量模块得到的测量结果经过所述Oracle验证单元的验证后与所述量子态变换模块构造的问题正确解的目标状态不一致,则重新调用量子态制备模块、量子态变换模块、量子态更新模块、测量模块及结果验证模块,直至所述测量结果与问题正确解的目标状态一致,输出所述测量结果。
7.一种电子装置,包括:存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时,实现权利要求1至4中的任意一项所述方法。
8.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时,实现权利要求1至3中的任意一项所述方法。
CN201910286148.9A 2019-04-10 2019-04-10 一种量子搜索方法、系统、电子装置及存储介质 Active CN110162536B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910286148.9A CN110162536B (zh) 2019-04-10 2019-04-10 一种量子搜索方法、系统、电子装置及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910286148.9A CN110162536B (zh) 2019-04-10 2019-04-10 一种量子搜索方法、系统、电子装置及存储介质

Publications (2)

Publication Number Publication Date
CN110162536A CN110162536A (zh) 2019-08-23
CN110162536B true CN110162536B (zh) 2021-07-16

Family

ID=67639067

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910286148.9A Active CN110162536B (zh) 2019-04-10 2019-04-10 一种量子搜索方法、系统、电子装置及存储介质

Country Status (1)

Country Link
CN (1) CN110162536B (zh)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110738321B (zh) * 2019-10-15 2022-04-29 北京百度网讯科技有限公司 一种量子信号处理方法及装置
CN111062482B (zh) * 2019-11-28 2023-07-04 南方科技大学 一种量子态重构方法、装置、系统和存储介质
CN111598245B (zh) * 2020-03-12 2023-04-07 西北大学 一种基于Grover算法的大数据集搜索方法、量子计算机
EP3886003A1 (en) 2020-03-25 2021-09-29 Terra Quantum AG Method, apparatus and computer program product for determining the component of a magnetic field in a predetermined direction
CN111241356B (zh) 2020-04-26 2020-08-11 腾讯科技(深圳)有限公司 基于模拟量子算法的数据搜索方法、装置及设备
CN111626422B (zh) * 2020-05-11 2023-08-22 深圳大学 一种用于求解集合交集的量子系统和算法
CN112132286B (zh) * 2020-08-14 2024-04-19 西北大学 基于Grover算法的大数据集搜索多个目标项的方法
CN112182494B (zh) * 2020-09-27 2023-04-28 中国人民解放军战略支援部队信息工程大学 基于Grover量子计算搜索算法的整数分解优化仿真方法及系统
CN112270411B (zh) * 2020-10-22 2022-05-24 苏州浪潮智能科技有限公司 量子退相干程度检测方法、装置及电子设备和存储介质
CN112966456B (zh) * 2021-03-05 2023-08-25 深圳大学 一种提高量子计算效率的电路设计方法和系统
CN113434723B (zh) * 2021-05-26 2023-10-10 中国电子技术标准化研究院 基于量子Grover算法的图像检索方法、装置和介质
CN116136968B (zh) * 2021-11-16 2024-06-21 本源量子计算科技(合肥)股份有限公司 Grover算法的量子线路优化方法及相关装置
CN116702914B (zh) * 2023-07-07 2024-06-25 北京师范大学 一种基于量子电路的模式匹配方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105531725A (zh) * 2013-06-28 2016-04-27 D-波系统公司 用于对数据进行量子处理的系统和方法
CN107742156A (zh) * 2017-10-18 2018-02-27 浙江工商大学 基于优化量子随机行走搜索算法的量子相干性方法
WO2018220368A1 (en) * 2017-05-30 2018-12-06 Gtn Ltd Tensor network machine learning system
CN109214517A (zh) * 2018-10-16 2019-01-15 深圳职业技术学院 量子状态搜索方法及装置

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10318881B2 (en) * 2013-06-28 2019-06-11 D-Wave Systems Inc. Systems and methods for quantum processing of data
CN103679751B (zh) * 2013-12-03 2016-08-17 大连大学 一种基于量子计算的碰撞检测方法
CN106023242A (zh) * 2015-04-09 2016-10-12 广东易富网络科技有限公司 一种基于量子均值漂移的抗遮挡多运动车辆追踪方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105531725A (zh) * 2013-06-28 2016-04-27 D-波系统公司 用于对数据进行量子处理的系统和方法
WO2018220368A1 (en) * 2017-05-30 2018-12-06 Gtn Ltd Tensor network machine learning system
CN107742156A (zh) * 2017-10-18 2018-02-27 浙江工商大学 基于优化量子随机行走搜索算法的量子相干性方法
CN109214517A (zh) * 2018-10-16 2019-01-15 深圳职业技术学院 量子状态搜索方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Quantum Rough Counting and Its Application to Grover"s Search Algorithm;Naphan Benchasattabuse等;《2018 3rd International Conference on Computer and Communication Systems》;20180430;第21-24页 *
基于Grover算法的量子处理架构的设计与仿真研究;代永涛;《中国优秀硕士学位论文全文数据库 信息科技辑(月刊)2017年》;20171115(第11期);I137-72 *

Also Published As

Publication number Publication date
CN110162536A (zh) 2019-08-23

Similar Documents

Publication Publication Date Title
CN110162536B (zh) 一种量子搜索方法、系统、电子装置及存储介质
Zhang et al. A novel quantum representation for log-polar images
Beilina et al. Numerical linear algebra: Theory and applications
WO2018075655A1 (en) Exact quantum circuits and circuit syntheses for qudit and multiple qubit circuits
Huber et al. A randomized tensor train singular value decomposition
Li et al. Design of threshold segmentation method for quantum image
Oh et al. High-performance tucker factorization on heterogeneous platforms
Beggs et al. Poisson–Riemannian geometry
Neumayer et al. Approximation of Lipschitz functions using deep spline neural networks
Mohammadisiahroudi et al. Efficient use of quantum linear system algorithms in interior point methods for linear optimization
Chiang A Note on the⊤‐Stein Matrix Equation
CN111291892A (zh) 一种量子并行搜索方法
Wang et al. An efficient numerical method for the symmetric positive definite second-order cone linear complementarity problem
Acus et al. The characteristic polynomial in calculation of exponential and elementary functions in Clifford algebras
Zubair et al. Multigrid for high-dimensional elliptic partial differential equations on non-equidistant grids
Liu et al. Preintegration via active subspace
Song et al. A unified scheme of resnet and softmax
Zhang et al. Accelerating exact nearest neighbor search in high dimensional Euclidean space via block vectors
Zhang et al. Accelerated PALM for nonconvex low-rank matrix recovery with theoretical analysis
Ren et al. Randomized quaternion QLP decomposition for low-rank approximation
Cui et al. Quantum radial basis function method for the Poisson equation
Zhang et al. A generalized floating-point quantum representation of 2-D data and their applications
CN111832734B (zh) 一种量子图像乘法运算的设计方法及其仿真实现方法
Joswig et al. Plücker coordinates and lines in space
Bujanović et al. A contribution to the theory and practice of the block Kogbetliantz method for computing the SVD

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant