CN116629370B - 量子电路处理方法、装置及电子设备 - Google Patents
量子电路处理方法、装置及电子设备 Download PDFInfo
- Publication number
- CN116629370B CN116629370B CN202310761796.1A CN202310761796A CN116629370B CN 116629370 B CN116629370 B CN 116629370B CN 202310761796 A CN202310761796 A CN 202310761796A CN 116629370 B CN116629370 B CN 116629370B
- Authority
- CN
- China
- Prior art keywords
- directed
- matrix
- list
- acyclic graph
- instruction
- 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
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 18
- 239000011159 matrix material Substances 0.000 claims abstract description 321
- 238000005457 optimization Methods 0.000 claims abstract description 49
- 238000012545 processing Methods 0.000 claims abstract description 23
- 239000002096 quantum dot Substances 0.000 claims description 92
- 238000000034 method Methods 0.000 claims description 40
- 230000006870 function Effects 0.000 claims description 24
- 230000000295 complement effect Effects 0.000 claims description 9
- 238000005259 measurement Methods 0.000 description 27
- 238000004422 calculation algorithm Methods 0.000 description 18
- 238000010586 diagram Methods 0.000 description 16
- 230000008569 process Effects 0.000 description 12
- 238000004088 simulation Methods 0.000 description 11
- 238000004590 computer program Methods 0.000 description 10
- 238000004891 communication Methods 0.000 description 9
- 238000013178 mathematical model Methods 0.000 description 5
- 238000010845 search algorithm Methods 0.000 description 5
- 230000003068 static effect Effects 0.000 description 5
- 238000010276 construction Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000001351 cycling effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 125000002015 acyclic group Chemical group 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000006837 decompression Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000005283 ground state Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000001953 sensory effect Effects 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N10/00—Quantum computing, i.e. information processing based on quantum-mechanical phenomena
- G06N10/20—Models of quantum computing, e.g. quantum circuits or universal quantum computers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N10/00—Quantum computing, i.e. information processing based on quantum-mechanical phenomena
- G06N10/60—Quantum algorithms, e.g. based on quantum optimisation, quantum Fourier or Hadamard transforms
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D10/00—Energy efficient computing, e.g. low power processors, power management or thermal management
Landscapes
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Devices For Executing Special Programs (AREA)
Abstract
本公开提供了一种量子电路处理方法、装置及电子设备,涉及量子计算技术领域,具体涉及量子电路技术领域。具体实现方案为:获取第一量子电路的第一指令列表;基于第一指令列表,确定第一有向无环图;基于第一有向无环图,确定第一有向无环图等效的有向二分图的第一邻接矩阵,第一邻接矩阵用于表征有向二分图的节点集合中每两个节点之间是否存在有向通路,节点集合包括第一有向无环图中的输入节点和输出节点;基于第一邻接矩阵,确定在有向二分图中添加输出节点至输入节点的有向边数量达到极大值的优化目标下的添加有向边矩阵;基于添加有向边矩阵、第一有向无环图和第一指令列表,对第一量子电路进行等效编译,得到第二量子电路的第二指令列表。
Description
技术领域
本公开涉及量子计算技术领域,尤其涉及量子电路技术领域,具体涉及一种量子电路处理方法、装置及电子设备。
背景技术
量子计算利用量子世界中特有的运行规律,提供了一条全新的并且非常有前景的信息处理方式。在诸多特定问题上,量子算法可以带来超越经典算法的优势。例如,利用秀尔(Shor)算法,可以对大整数进行高效的分解,利用格罗弗(Grover)算法,可以更快的进行数据搜索。随着量子理论的发展,不断有新的量子算法被提出,如何高效的对这些算法进行模拟或者在真正的量子硬件上运行始终是一个重要的问题。
目前,量子算法的经典模拟或者真机运行主要受限于量子比特的数量。在经典模拟上,由于描述量子态的列向量的长度随对应比特数呈指数增长(例如一个n比特的量子态的列向量长度为2n),经典计算机很难模拟大规模的量子算法。受计算机内存和处理器能力的限制,现有的量子电路模拟方式最多能支持模拟几十个量子比特的算法。
发明内容
本公开提供了一种量子电路处理方法、装置及电子设备。
根据本公开的第一方面,提供了一种量子电路处理方法,包括:
获取第一量子电路的第一指令列表;
基于所述第一指令列表,确定第一有向无环图,所述第一有向无环图包括所述第一指令列表中指令对应的节点和至少两个第一有向边,所述第一有向边用于表征所述第一指令列表中不同指令间的时序关系;
基于所述第一有向无环图,确定所述第一有向无环图等效的有向二分图的第一邻接矩阵,所述第一邻接矩阵用于表征所述有向二分图的节点集合中每两个节点之间是否存在有向通路,所述节点集合包括所述第一有向无环图中的输入节点和输出节点;
基于所述第一邻接矩阵,确定在所述有向二分图中添加输出节点至输入节点的有向边数量达到极大值的优化目标下的添加有向边矩阵,基于所述添加有向边矩阵在所述有向二分图添加有向边的情况下,添加有向边后的所述有向二分图中不存在有向环路,且所述添加有向边矩阵满足:每个输出节点最多只能和一个输入节点相连,以及每个输入节点最多只能和一个输出节点相连;
基于所述添加有向边矩阵、所述第一有向无环图和所述第一指令列表,对所述第一量子电路进行等效编译,得到与所述第一量子电路等效的第二量子电路的第二指令列表。
根据本公开的第二方面,提供了一种量子电路处理装置,包括:
获取模块,用于获取第一量子电路的第一指令列表;
第一确定模块,用于基于所述第一指令列表,确定第一有向无环图,所述第一有向无环图包括所述第一指令列表中指令对应的节点和至少两个第一有向边,所述第一有向边用于表征所述第一指令列表中不同指令间的时序关系;
第二确定模块,用于基于所述第一有向无环图,确定所述第一有向无环图等效的有向二分图的第一邻接矩阵,所述第一邻接矩阵用于表征所述有向二分图的节点集合中每两个节点之间是否存在有向通路,所述节点集合包括所述第一有向无环图中的输入节点和输出节点;
第三确定模块,用于基于所述第一邻接矩阵,确定在所述有向二分图中添加输出节点至输入节点的有向边数量达到极大值的优化目标下的添加有向边矩阵,基于所述添加有向边矩阵在所述有向二分图添加有向边的情况下,添加有向边后的所述有向二分图中不存在有向环路,且所述添加有向边矩阵满足:每个输出节点最多只能和一个输入节点相连,以及每个输入节点最多只能和一个输出节点相连;
等效编译模块,用于基于所述添加有向边矩阵、所述第一有向无环图和所述第一指令列表,对所述第一量子电路进行等效编译,得到与所述第一量子电路等效的第二量子电路的第二指令列表。
根据本公开的第三方面,提供了一种电子设备,包括:
至少一个处理器;以及
与至少一个处理器通信连接的存储器;其中,
存储器存储有可被至少一个处理器执行的指令,该指令被至少一个处理器执行,以使至少一个处理器能够执行第一方面中的任一项方法。
根据本公开的第四方面,提供了一种存储有计算机指令的非瞬时计算机可读存储介质,该计算机指令用于使计算机执行第一方面中的任一项方法。
根据本公开的第五方面,提供了一种计算机程序产品,包括计算机程序,该计算机程序在被处理器执行时实现第一方面中的任一项方法。
根据本公开的技术解决了相关技术中量子电路的经典模拟和真机运行比较难的问题,使得能够实现对大规模量子比特的量子电路进行经典模拟和真机运行。
应当理解,本部分所描述的内容并非旨在标识本公开的实施例的关键或重要特征,也不用于限制本公开的范围。本公开的其它特征将通过以下的说明书而变得容易理解。
附图说明
附图用于更好地理解本方案,不构成对本公开的限定。其中:
图1是根据本公开第一实施例的量子电路处理方法的流程示意图;
图2是一示例的量子电路图的结构示意图;
图3是第一量子电路的结构示意图;
图4是第一有向无环图的结构示意图;
图5是第二量子电路的结构示意图;
图6是一示例的量子电路处理方法的具体流程示意图;
图7是根据本公开第二实施例的量子电路处理装置的结构示意图;
图8是用来实施本公开的实施例的示例电子设备的示意性框图。
具体实施方式
以下结合附图对本公开的示范性实施例做出说明,其中包括本公开实施例的各种细节以助于理解,应当将它们认为仅仅是示范性的。因此,本领域普通技术人员应当认识到,可以对这里描述的实施例做出各种改变和修改,而不会背离本公开的范围和精神。同样,为了清楚和简明,以下的描述中省略了对公知功能和结构的描述。
第一实施例
如图1所示,本公开提供一种量子电路处理方法,包括如下步骤:
步骤S101:获取第一量子电路的第一指令列表。
本实施例中,量子电路处理方法涉及量子计算技术领域,尤其涉及量子电路技术领域,其可以广泛应用于量子电路的经典模拟和真机运行场景下。本公开实施例的量子电路处理方法,可以由本公开实施例的量子电路处理装置执行。本公开实施例的量子电路处理装置可以配置在任意电子设备中,以执行本公开实施例的量子电路处理方法。
受计算机内存和处理器能力的限制,现有的量子电路模拟方式最多能支持模拟几十个量子比特的算法。比如,笔记本能模拟20-30个左右的量子比特,大型超级计算机和集群可以最多模拟30-40个左右的量子比特。在真机运行上,由于当前量子芯片的可扩展性问题尚未解决,导致量子计算机能提供的量子比特数非常有限。因此,量子电路优化是量子计算领域中的一个基本问题。
量子电路优化是通过一定的技术手段,可以将给定的量子电路进行简化,以降低其经典模拟和真机运行的要求,进而加速量子算法的研究和量子计算在实际场景下的落地。
而本实施例的量子电路处理可以为量子电路优化的处理,其目的在于通过对量子电路进行优化编译,可以使得编译得到的量子电路在量子比特数上对原量子电路进行大量简化。一方面,可以进一步提升量子算法经典模拟的规模,加强经典计算机对量子算法的验证能力,另一方面,也可以降低量子算法在真机运行上的比特数要求,弥补当前量子芯片可扩展性问题的不足。
更具体的,本实施例的目的是将给定的静态量子电路等效编译为动态量子电路,以减少量子电路的量子比特数量,降低其经典模拟和真机运行的要求。
其中,通过引入量子比特的中间电路测量以及在计算中将量子比特重置为其基态的能力,动态量子电路能够将经典通信和计算整合到量子电路中。这一特性将极大地增加在量子硬件上运行电路的多样性,并且对于纠错和容错量子计算的发展也至关重要。因此,动态量子电路有望成为未来许多量子应用的重要组成部分。在数学意义上和没有实验约束的情况下,静态和动态量子电路具有相同的计算能力,即任何可以用动态量子电路表示的问题也可以表示为静态问题,所以动态量子电路提供了权衡电路深度和宽度的可能性。
以下详细介绍量子电路模型。
量子电路模型是一种常用的量子计算模型。通过对初始量子态进行量子门操作完成量子态的演化,并通过量子测量提取计算结果。而量子电路图则表示了量子电路模型计算的全部过程。
图2是一示例的量子电路图的结构示意图,如图2所示,可以用一根水平线表示一个量子比特系统,从上至下依次对量子比特位进行标号,其中,量子位的标号往往从零开始。
量子电路图中时间演化的方向从左到右,最左端为初始的量子态,其中,通常每个量子比特初始化为零态,之后对初始态依次作用不同的量子门操作以完成量子态的演化。同时可以对某些量子位进行量子测量,获得测量结果。
在一些应用场景中,量子电路中的操作可能会出现对一部分量子比特进行量子测量,并根据测量结果来调控其余的量子比特的演化,此类操作称为经典控制量子操作,如图2所示的经典控制量子门201。而经过测量后的量子比特则可以被重置,该类操作可以称之为重置操作,如图2所示的重置操作202,以供后续计算继续使用。可以将包含中间测量、经典控制量子操作以及重置操作的量子电路称为动态量子电路,如图2表示的量子电路即为动态量子电路。
量子电路图中除了初始态之外其余部分,通常可以按照量子门的作用顺序用一个有序的指令列表进行表示,列表中的每一个元素代表一个量子门、经典控制量子门、量子测量或重置操作指令。具体地,可以将:
每一个单量子比特门(如H,X,Y,Z,S,T,Rx,Ry,Rz等)表示为一个包含四个元素的指令[name,which_qubit,parameters,condition]。其中,name为量子门的名称,which_qubit为量子门作用的量子位,parameters为量子门的参数(如没有参数则默认为None),condition表示该量子门操作受哪一个量子位的测量结果控制(如果没有参数则默认为None)。
例如,[Rx,2,pi,None]表示对量子位2上的量子比特作用一个Rx旋转门,旋转角度为pi。又例如,图2中经典控制量子门201为经典受控量子X门,可以表示为[X,2,None,‘a’],即作用在量子位2上的泡利Pauli X门,受控条件为测量标识ID为‘a’的测量结果,测量结果为0则不作用量子门,测量结果为1则作用量子门。
每一个双量子比特门(如控制非门CNOT,CZ门)表示为一个包含四个元素的指令[name,which_qubit,parameters,condition]。其中,name为量子门的名称,which_qubit为该双量子比特门作用的量子位构成的列表(特别地,对于受控量子门,为控制位和受控位构成的列表),parameters为量子门的参数(如没有参数则默认为None),condition表示该量子门操作受哪一个量子位的测量结果控制(如果没有参数则默认为None)。
例如,[CNOT,[1,3],None,None]表示对量子位1和量子位3作用控制非门,其中,量子位1为控制位,量子位3为受控位。[CZ,[1,2],None,None]表示在量子位1和量子位2之间作用CZ门。
每一个单比特测量表示为一个包含四个元素的指令[measure,which_qubit,basis,mid]。其中,basis由四个参数决定,其包括测量角度angle,测量平面plane,域集s,域集t,而mid为标识当前测量的标识ID。
例如,[measure,2,[0,‘YZ’,[1],[2]],‘a’]表示对量子位2进行测量,测量角度为0,测量平面为‘YZ’平面,域集s为量子位1,域集t为量子位2,当前测量指令的标识ID为‘a’。
每一个重置操作指令可以表示为一个包含四个元素的指令[reset,which_qubit,matrix,None]。其中,which_qubit为需要重置的量子位,matrix为需要重置比特的量子态矩阵,经过重置操作后的量子比特可供后续计算继续使用。
在步骤S101中,第一量子电路可以为静态量子电路,静态量子电路中仅包括量子测量操作和量子门操作,且量子测量操作位于量子门操作之后。
第一指令列表可以基于量子电路优化前的原量子电路即第一量子电路的指令列表得到,通过在第一量子电路的指令列表之前(即指令列表的首位位置)添加各个量子比特的重置操作指令,即可得到第一指令列表。
步骤S102:基于所述第一指令列表,确定第一有向无环图,所述第一有向无环图包括所述第一指令列表中指令对应的节点和至少两个第一有向边,所述第一有向边用于表征所述第一指令列表中不同指令间的时序关系。
其中,所述至少两个第一有向边构成的路径不包括有向环路。
该步骤中,第一有向无环图可以包括节点和由节点构成的第一有向边,每个节点表示第一指令列表中的一个操作指令,其中,节点标识可以为第一指令列表中的指令,也可以为第一指令列表中指令对应的标识,这里不进行具体限定。
第一有向边指的是带有指向(即方向性)的边,表示第一指令列表中不同指令间的时序关系,如节点A指向节点B的边表示节点A对应指令需要在节点B对应指令之前完成执行。
第一有向无环图中,各个第一有向边所构成的任一路径均不是有向环路,有向环路指的是带有方向性的环路,如节点A、节点B和节点C,若节点A指向节点B,节点B指向节点C、节点C再指向节点A,则节点A、节点B和节点C构成了一条有向环路。
由于量子电路中指令的时序关系可以基于量子电路的指令列表中指令的排列顺序体现,因此,可以基于第一量子电路的第一指令列表中指令的排列顺序,构建第一有向无环图。
在一可选实施方式中,可以基于第一指令列表中指令的数量,构建相应数量的节点,并基于第一指令列表中指令的排列顺序,确定各个指令之间的时序关系,基于该时序关系,在所构建的节点之间连接有向边,从而可以得到第一有向无环图。
在另一可选实施方式中,可以进行第一指令列表中指令的循环遍历,针对当前遍历的指令,在该当前遍历的指令之前查找与之存在时序关系的其他指令,相应构建存在时序关系的节点,并构建一条有向边,将其添加至图中,在循环遍历完成的情况下,即可以得到第一有向无环图。
其中,第一有向无环图中包括输入节点和输出节点,输入节点指的是第一指令列表中类型为重置操作指令对应的节点,输出节点指的是第一指令列表中类型为量子测量操作指令对应的节点。
步骤S103:基于所述第一有向无环图,确定所述第一有向无环图等效的有向二分图的第一邻接矩阵,所述第一邻接矩阵用于表征所述有向二分图的节点集合中每两个节点之间是否存在有向通路,所述节点集合包括所述第一有向无环图中的输入节点和输出节点。
步骤S104:基于所述第一邻接矩阵,确定在所述有向二分图中添加输出节点至输入节点的有向边数量达到极大值的优化目标下的添加有向边矩阵,基于所述添加有向边矩阵在所述有向二分图添加有向边的情况下,添加有向边后的所述有向二分图中不存在有向环路,且所述添加有向边矩阵满足:每个输出节点最多只能和一个输入节点相连,以及每个输入节点最多只能和一个输出节点相连。
本实施例中,可以将量子电路的优化问题转化为解决一个等价的图论优化问题,由于将一个量子电路编译成另一个量子电路,本质上是通过对已经测量的量子比特进行重置回收,以供后续指令所使用。对应第一有向无环图,相当于需要尽可能的向图中添加从输出节点指向输入节点的有向边,以表示测量比特重置被重新利用,相应的,可以减少量子电路中所需要的量子比特数。
为了保证所编译得到的量子电路与原量子电路的等效性,以及保证量子电路中指令之间的时序关系,在添加有向边时需要保证如下约束条件:
1)为了使得整个图依然符合电路指令的时序关系,需要保证在加入更多有向边之后,整个图依然是无环的,需要注意的是,图中有向边的指向表示指令的时序关系,环状结构将会破坏这种关系,也即添加有向边之后,图中依然不能存在环状结构的路径(称之为有向环路);
2)每个输出节点最多只能和一个输入节点相连,以表示输出节点测量后,由所指向的输入节点继续在该量子比特上执行操作;
3)不同输出节点不能和同一个输入节点相连。
本实施例中,通过将量子电路编译转化为对第一有向无环图等效的有向二分图的第一邻接矩阵的处理,并约束添加有向边矩阵,以求解得到在有向二分图中添加输出节点至输入节点的有向边数量达到极大值的优化目标下的添加有向边矩阵,即查找到在有向二分图中添加最多有向边的可行解。
其中,步骤S103中其目的是获取第一有向无环图等效的有向二分图的第一邻接矩阵。步骤S104中其目的是通过对第一邻接矩阵的处理,并约束添加有向边矩阵,以求解得到在有向二分图中添加输出节点至输入节点的有向边数量达到极大值的优化目标下的添加有向边矩阵。
其中,有向二分图中仅包括第一有向无环图的输入节点和输出节点,即有向二分图的两个节点子集分别为量子电路的输入节点集和输出节点集,两个节点子集构成了有向二分图的节点集合。可以舍弃第一有向无环图的中间节点,将第一有向无环图转化为只包含输入节点和输出节点的等效有向二分图。
第一有向无环图等效的有向二分图也是第一量子电路的等效有向二分图,其中,有向二分图的第一邻接矩阵为2n×2n的矩阵,n为输入节点数量,也即输出节点数量,其表征的是有向二分图的节点集合中每两个节点之间是否存在有向通路。
可以用第一目标值表征每两个节点之间存在有向通路,用第二目标值表征每两个节点之间不存在有向通路。其中,第一目标值可以为1,第二目标值可以为0,第一目标值也可以为0,第二目标值也可以为1。以下实施例中,第一目标值将以1为例进行详细说明。
在一可选实施方式中,可以基于第一有向无环图,利用深度优先搜索算法确定每个输入节点到所有输出节点是否可达,在可达的情况下,即说明输入节点到输出节点之间存在有向通路,在不可达的情况下,即说明输入节点到输出节点之间不存在有向通路。相应可以得到有向二分图的第一邻接矩阵。
在另一可选实施方式中,可以基于第一有向无环图,利用深度优先搜索算法确定每个输入节点到所有输出节点是否可达,得到可达性矩阵,将可达性矩阵作为邻接矩阵中右上角的块矩阵,从而得到有向二分图的第一邻接矩阵。该可选实施方式中,等效有向二分图的第一邻接矩阵可表示为一个2×2的反对角分块矩阵,其右上角分块可以为可连接性矩阵的互补矩阵,可连接性矩阵的互补矩阵即为可达性矩阵,其他位置的分块均为n×n的全零矩阵。其中,互补可以指的是两个矩阵中相同位置的任意元素值互补,如两个元素值分别为0和1,则表示两个元素值互补。
在又一可选实施方式中,可以基于第一有向无环图,利用深度优先搜索算法确定每个输入节点到所有输出节点是否可达,以确定是否存在可连接的有向边,并生成可连接性矩阵,基于可连接性矩阵确定有向二分图的第一邻接矩阵。
在步骤S104中,可以基于第一邻接矩阵,约束添加有向边矩阵,使得基于所述添加有向边矩阵在有向二分图添加有向边的情况下,添加有向边后的所述有向二分图中不存在有向环路,且每个输出节点最多只能和一个输入节点相连,以及每个输入节点最多只能和一个输出节点相连。之后在约束添加有向边矩阵的条件下确定满足优化目标的添加有向边矩阵。
其中,基于添加有向边矩阵在有向二分图中添加输出节点至输入节点的有向边数量达到极大值时,则说明添加有向边矩阵满足优化目标。
在对添加有向边矩阵的约束下,且满足优化目标的情况下,可以将量子电路优化问题转化为一个受约束的二元整数规划问题,利用程序语言(如python)的数学工具包(如CVXOPT)或其他整数规划问题求解器求解该二元整数规划问题,以得到二元整数规划问题的可行解,该可行解即为约束下且满足优化目标的添加有向边矩阵。
步骤S105:基于所述添加有向边矩阵、所述第一有向无环图和所述第一指令列表,对所述第一量子电路进行等效编译,得到与所述第一量子电路等效的第二量子电路的第二指令列表。
该步骤中,在一可选实施方式中,可以基于添加有向边矩阵,在第一有向无环图中添加输出节点至输入节点的第二有向边,得到第二有向无环图;之后基于第二有向无环图、第二有向边构成的第一目标列表和第一指令列表,对第一量子电路进行等效编译,得到与第一量子电路等效的第二量子电路的第二指令列表。
添加输出节点至输入节点的有向边之后,可以基于第二有向无环图,确定各个节点之间的时序关系,并基于该时序关系,对第一指令列表进行重排序,之后,可以利用第一目标列表中输入节点与输出节点之间的关系,将测量指令之后作用在其他量子比特(即与输出节点有向连接的输入节点对应指令所作用的量子比特)的指令,等效编译成作用到该测量指令所作用的量子比特上,如此,可以实现对第一量子电路的等效编译,得到与第一量子电路等效的第二量子电路的第二指令列表,第二量子电路中所需要的量子比特比第一量子电路中所需要的量子比特少。
在另一可选实施方式中,也可以基于添加有向边矩阵和第一有向无环图,直接对第一指令列表进行重新排序,以将添加有向边矩阵所表征的有向边中输入节点对应的量子比特操作指令排序到所述有向边中输出节点对应的量子比特操作指令之后,之后再进行等效编译,这样可以在等效编译过程中将已经测量的量子比特进行重置回收,以供后续指令所使用,从而使得编译得到的量子电路所需的量子比特数减少。
需要说明的是,由于添加有向边矩阵是在满足优化目标下的添加有向边矩阵,因此,可以为第一有向无环图添加最多的输出节点至输入节点的有向边,这样在测量的量子比特被重置并重新利用后,可以使得编译得到的量子电路即第二量子电路的量子比特数比原量子电路的量子比特数少,从而可以达到优化的目的。
本实施例中,通过将量子电路编译转化为量子电路的等效有向二分图的邻接矩阵的处理,以约束添加有向边矩阵,并在限定基于添加有向边矩阵在有向二分图中添加输出节点至输入节点的有向边数量达到极大值的优化目标下,求解满足该优化目标的添加有向边矩阵,找到量子电路编译问题的一个可行解,基于该可行解进行量子电路的等效编译所得到的量子电路所需的量子比特数比原量子电路的量子比特数少,因此可以减少量子电路所需的量子比特数,从而实现量子电路的优化,使得能够实现对大规模量子比特的量子电路进行经典模拟和真机运行。
可选的,所述步骤S102具体包括:
对所述第一指令列表按照从左至右的指令排列顺序进行遍历,并获取第二目标列表和当前遍历的指令所作用的第三量子位,所述第二目标列表包括所述第一指令列表中排序在所述当前遍历的指令之前的指令;
对所述第二目标列表按照从右至左的指令排列顺序进行遍历,并获取第一目标指令,所述第一目标指令所作用的量子位与所述第三量子位存在交集;
将所述第一目标指令的标识作为源节点,所述当前遍历的指令的标识作为目标节点,构建所述第一有向边;
在所述第一指令列表遍历完成的情况下,基于所构建的第一有向边,获取所述第一有向无环图。
本实施方式中,量子电路至第一有向无环图的转化过程如下:
输入:量子电路指令列表circuit_list;
输出:有向无环图(为第一有向无环图),输入节点列表,输出节点列表。
步骤1:根据量子电路指令列表circuit_list所表示的电路宽度,记录为n;循环遍历列表[0,1,…,n-1],记当前被循环的元素为idx,生成电路指令gate=[reset,idx,[[1],[0]],None],并添加至circuit_list列表最前面,获得更新后的列表circuit_list,即得到第一指令列表;
步骤2:循环遍历列表circuit_list,记当前被循环的元素为gate,且gate处于circuit_list的第i位,为gate添加一个ID属性,ID的值为i;
步骤3:初始化一个有向图graph,初始化两个空列表input_ids和output_ids;
步骤4:循环遍历列表circuit_list(即按照从左至右的指令排列顺序遍历),记当前被循环的元素为gate,且gate处于circuit_list的第i位,执行如下操作a)、b):
操作a)记录circuit_list中的前i个元素构成的列表为pre_gates(即第二目标列表);
操作b)获取gate指令所作用的量子位(即第三量子位),记当前被循环的量子位为idx;
操作c)反向循环遍历pre_gates(即按照从右至左的指令排列顺序遍历),并记当前被循环的元素为pre_gate;如果pre_gate所作用的量子位与idx有交集,则pre_gate为第一目标指令;记node1为pre_gate指令的标识ID,node2为gate指令的标识ID,向有向图graph中构建并添加节点node1和node2,添加一条从node1(源节点)指向node2(目标节点)的有向边;之后跳出本层循环;
步骤5:如果gate为重置指令,则将其对应ID加入到input_ids列表最末尾;如果gate为测量指令,则将其对应ID加入到output_ids列表最末尾;
步骤6:返回有向图graph(第一有向无环图),输入节点列表input_ids和输出节点列表output_ids。
本实施方式中,通过进行第一指令列表中指令的循环遍历,针对当前遍历的指令,在该当前遍历的指令之前查找与之存在时序关系(即两个指令作用有相同量子位)的其他指令,相应构建存在时序关系的节点,并构建一条有向边,将其添加至图中,在循环遍历完成的情况下,即可以得到第一有向无环图。如此,可以简单地实现第一有向无环图的构建。
另外,在构建第一有向无环图时,仅使用对应指令的ID即可,无需用电路指令本身作为图的节点,以减少相应运算。
在一示例中,量子电路的结构示意图如图3所示,其中,H表示H门,M表示测量measure,*-X表示CNOT门。
添加重置操作指令之后,第一量子电路的第一指令列表如下为[[reset,0,[[1],[0]],None];[reset,1,[[1],[0]],None];[reset,2,[[1],[0]],None];[h,0,None,None];[h,1,None,None];[h,2,None,None];[cnot,[0,1],None,None];[cnot,[1,2],None,None];[measure,0,[0,‘YZ’,[],[]],(0,1)];[measure,1,[0,‘YZ’,[],[]],(1,1)];[measure,2,[0,‘YZ’,[],[]],(2,1)]]。
按照量子电路至第一有向无环图的转化过程,为每个指令分配的ID依次为0,1,…,10,得到的第一有向无环图如图4所示。其中,每个节点上的数字为对应指令的ID,节点401为输入节点,节点402为输出节点。
可选的,所述第一邻接矩阵按照所述第一有向无环图中的输入节点数量进行分块,所述第一邻接矩阵为2×2的分块矩阵,所述步骤S103具体包括:
基于所述第一有向无环图,确定所述第一有向无环图的可达性矩阵,所述可达性矩阵用于表征所述第一有向无环图中输入节点到输出节点的可达性,所述可达性指示输入节点是否存在到达输出节点的有向通路;
将所述可达性矩阵作为所述第一邻接矩阵中第一行第二列的块矩阵,得到所述第一邻接矩阵。
本实施方式中,根据约束条件1):为了使得整个图依然符合电路指令的时序关系,需要保证在加入更多有向边之后,整个图依然是无环的(注意,图中的边的指向表示指令的时序关系,环状结构将会破坏这种关系)。
为了保证添加有向边后满足约束条件1),可以利用深度优先搜索算法确定每个输入节点到所有输出节点是否可达,从而避免添加从输出节点到输入节点的有向边而导致的环状结构。
可达性矩阵为二维矩阵,二维矩阵的行标表示输入节点的标号,列标为输出节点的标识,相应元素表示输入节点至输出节点是否可达,即输入节点是否存在到达输出节点的通路,若存在通路,则可达,若不存在通路,则不可达。
可连接性矩阵也可以为二维矩阵,二维矩阵的行标表示输入节点的标号,列标为输出节点的标识,相应元素可以表示输出节点与输入节点的可连接性,即是否可以在输出节点与输入节点之间添加有向边。
可连接性矩阵与可达性矩阵的指示相反,其针对相同行标和列表的元素,若可达性矩阵中元素值表示输入节点至输出节点不可达,则可连接性矩阵中元素值表示输出节点至输入节点可连接,若可达性矩阵中元素值表示输入节点至输出节点可达,则可连接性矩阵中表示输出节点至输入节点不可连接。
输出节点至输入节点的可连接性的搜索过程如下:
输入:有向无环图graph(可以为第一有向无环图),输入节点列表inputs,输出节点列表outputs;
输出:二维矩阵all_candidate(即可连接性矩阵),包含输出节点到输入节点的可连接性,元素的值为1时表示输出节点到输入节点可连接,元素的值为0时表示输出节点到输入节点不可连接。
步骤1:对输入节点列表中的每个元素,利用深度优先搜索(或者其他有效的算法)确定其到输出节点列表中每个输出节点是否可达。具体地,步骤1包括步骤c)和步骤d):
步骤c)初始化二维矩阵all_candidate,其行数等于输入节点列表的长度,列数等于输出节点列表的长度,所有元素的初始值为1(元素值为1时表示该输入节点到输出节点不可达,即输出节点到输入节点可连接;元素值为0表示该输入节点到输出节点可达,输出节点到输入节点不可连接);
步骤d)循环遍历输出节点列表outputs,记当前被遍历的元素为output:利用深度优先搜索(或者其他有效的算法)判断输入节点inputs中的每个输入节点在图graph中是否存在到达节点output的通路,若存在通路则记该输入节点为output节点的不可达节点,并将二维矩阵all_candidate对应行列下的元素标记为0;反之不改变二维矩阵all_candidate对应行列下的元素值(即保持1);
步骤2:输出矩阵all_candidate,即为可连接性矩阵。
对应图4中第一有向无环图的可连接性矩阵为如下式(1)所示,该可连接性矩阵为一个3×3的二维矩阵。
[[0,0,0]
[0,0,0](1)
[1,0,0]]
可连接性矩阵中包含了有向无环图中从任意输入节点到任意输出节点的连接性信息,由此可以舍弃中间节点,将第一有向无环图转化为只包含输入节点和输出节点的等效有向二分图。
该等效有向二分图的第一邻接矩阵包含2n×2n个元素,其中n为输入节点或输出节点的数量。第一邻接矩阵中前n行或前n列表示的节点为输入节点,第一邻接矩阵中后n行或后n列表示的节点为输出节点。
第一邻接矩阵的确定过程如下:
输入:输入节点列表inputs,输出节点列表outputs,可连接性矩阵all_candidate;
输出:等效有向二分图的邻接矩阵bipartite_adjacency(即第一邻接矩阵);
步骤1:初始化一个全零矩阵bipartite_adjacency,其行数和列数都为2n;
步骤2:将bipartite_adjacency矩阵右上角n×n个元素赋值为可连接性矩阵all_candidate中对应的n×n个元素的0/1互补值。
本实施方式中,可以基于第一有向无环图,利用深度优先搜索算法确定每个输入节点到所有输出节点是否可达,得到可达性矩阵,将可达性矩阵作为邻接矩阵中右上角的块矩阵,从而得到有向二分图的第一邻接矩阵。如此,可以简单地实现第一邻接矩阵的确定。
可选的,所述第一邻接矩阵为反对角分块矩阵,所述第一邻接矩阵中除所述可达性矩阵之外的块矩阵均为全零矩阵。
本实施方式中,当第一目标值为1时,等效有向二分图的第一邻接矩阵可表示为一个2×2的反对角分块矩阵,其右上角分块可以为可连接性矩阵的互补矩阵,其他位置的分块均为n×n的全零矩阵。其中,互补可以指的是两个矩阵中相同位置的任意元素值互补,如两个元素值分别为0和1,则表示两个元素值互补。
对应图4中的有向无环图,其等效有向二分图的第一邻接矩阵为如下式(2)所示。
[[0,0,0,1,1,1]
[0,0,0,1,1,1]
[0,0,0,0,1,1](2)
[0,0,0,0,0,0]
[0,0,0,0,0,0]
[0,0,0,0,0,0]]
如此,可以简化第一邻接矩阵中的元素值,进而简化后续基于第一邻接矩阵所构建的问题求解。
可选的,所述步骤S104具体包括:
基于所述添加有向边矩阵,构建优化目标的第一函数,所述第一函数表征所述添加有向边矩阵中元素值的和值达到极值,所述极值表征在所述有向二分图中添加输出节点至输入节点的有向边数量达到极大值;
基于所述第一邻接矩阵,构建所述添加有向边矩阵的约束条件;
基于所述约束条件和所述第一函数,求解所述优化目标下的添加有向边矩阵。
本实施方式中,可以基于添加有向边矩阵,构建优化目标的第一函数,可选的,在所述第一目标值为1的情况下,所述第一函数表征所述添加有向边矩阵中元素值的和值达到极大值,即第一函数可以为
而在第一目标值为0的情况下,第一函数可以为最小值min函数,第一函数也可以表征在有向二分图中添加输出节点至输入节点的有向边数量达到极大值。如此,可以实现优化目标的构建。
可以基于第一邻接矩阵,构建添加有向边矩阵的约束条件,使得基于所述添加有向边矩阵在所述有向二分图添加有向边的情况下,添加有向边后的所述有向二分图中不存在有向环路,且所述添加有向边矩阵满足:每个输出节点最多只能和一个输入节点相连,以及每个输入节点最多只能和一个输出节点相连。
可选的,所述约束条件包括:
所述添加有向边矩阵为所述第一有向无环图的可连接性矩阵转置的子矩阵,所述可连接性矩阵用于表征所述第一有向无环图中输出节点至输入节点的可连接性;
添加有向边后的所述有向二分图的第二邻接矩阵满足幂零;
所述添加有向边矩阵的每列元素值中最多存在一个第一目标值;
所述添加有向边矩阵的每行元素值中最多存在一个第一目标值,所述第一目标值表征在所述有向二分图中添加输出节点至输入节点的有向边。
在一可选实施方式中,如第一目标值为1时,得到等效有向二分图的第一邻接矩阵后,可以将量子电路编译问题转化为找到一个n×n元素组成的添加有向边矩阵的解Y′,使其满足下式(3)所示的优化目标,以及满足下式(4)、(5)和(6)所示的约束条件。
其中:Yij为可连接性矩阵的(j,i)元素,B为可连接性矩阵的互补矩阵,为添加输出节点至输入节点的有向边后的有向二分图的第二邻接矩阵,nilpotent为幂零(即:存在一个正整数k,使得矩阵的k次幂为全零矩阵)。
注意:添加有向边矩阵中(i,j)元素值只可为0或1,为1时表示添加从第i个输出节点到第j个输入节点的有向边,为0时表示不添加从第i个输出节点到第j个输入节点的有向边。
下面解释上述优化问题的数学模型,以及其与量子电路编译问题的等价性。
上式(3)为优化目标,表示添加有向边矩阵中需要有尽可能多的非零元素,对应量子电路编译问题需要尽可能减少量子比特数,添加尽可能多的有向边。
上式(4)表示优化问题的数学模型中的约束条件a),表示添加有向边矩阵为可连接性矩阵转置的子矩阵(即:添加有向边矩阵中非零元素只能取可连接性矩阵中非零元素的子集),对应量子电路编译问题中为了使得整个图依然符合电路指令的时序关系,需要保证在加入更多有向边之后,整个图依然是无环的。
上式(5)表示优化问题的数学模型中的约束条件b),表示基于添加有向边矩阵在有向二分图添加输出节点至输入节点的有向边后,其更新的有向二分图的第二邻接矩阵满足幂零,对应量子电路编译问题中为了使得整个图依然符合电路指令的时序关系,需要保证在加入更多有向边之后,整个图依然是无环的。
上式(6)表示优化问题的数学模型中的约束条件c),表示添加有向边矩阵中每行或每列最多只有一个非零元素,对应量子电路编译问题中每个输出节点最多只能和一个输入节点相连,以及每个输入节点最多只能和一个输出节点相连,以表示输出节点测量后,由所指向的输入节点继续在该量子比特上执行操作。
如此,可以实现约束条件的构建。
在构建优化目标和约束条件的情况下,可以将量子电路编译问题(即将原量子电路等效编译为动态量子电路,以减少所使用的量子比特数)转换为二元整数规划问题。可以利用程序语言(如python)的数学工具包(如CVXOPT)或其他整数规划问题求解器求解该二元整数规划问题,以得到二元整数规划问题的可行解,该可行解即为约束条件下且满足优化目标的添加有向边矩阵。该添加有向边矩阵可以是非全零的矩阵。如此,可以实现添加有向边矩阵的确定,从而实现输出节点至输入节点的连接方案的确定。
可选的,使添加有向边后的所述有向二分图的第二邻接矩阵满足幂零的幂为4。这样可以通过放缩法放缩幂的取值,以将量子电路编译问题转换为二次二元整数规划问题,以降低问题求解的难度和时间,从而提高量子电路编译的效率。
可选的,在所述第一目标值为1的情况下,所述约束条件包括:
BY′BY′=0;
Y′BY′B=0;
其中,Yij为所述可连接性矩阵中索引为(j,i)的元素,B为所述第一有向无环图的可达性矩阵,所述可达性矩阵中的元素值与所述可连接性矩阵中的元素值互补,Y′为添加有向边矩阵,Y′ij为添加有向边矩阵中索引为(i,j)的元素,n为所述第一有向无环图中的输入节点数量。
本实施方式中,可以通过放缩法放缩幂的取值,使得k=4,以求解添加有向边矩阵的非全零可行解,即找到一个n×n元素组成的添加有向边矩阵的非全零可行解Y′,优化问题和约束条件如下:
下面解释通过放缩法来优化问题的数学模型。
优化目标:添加有向边矩阵中有尽可能多的非零元素。
约束条件:
约束条件a),添加有向边矩阵为可连接性矩阵转置的子矩阵;
约束条件b),令k=4,即令添加有向边后的有向二分图的第二邻接矩阵的四次幂为全零矩阵。其中,令k=4是因为k<4时可能无法确保可求得问题的非全零可行解。
约束条件c),添加添加有向边矩阵每行或每列最多只有一个非零元素。
利用放缩法将量子电路编译问题转化为一个受约束的二次二元整数规划问题,可利用python CVXOPT或其他整数规划问题求解器求解,从而可以降低问题的求解难度,降低量子电路编译过程的难度和时间,从而可以提高量子电路处理的效率。
可选的,所述步骤S105具体包括:
基于所述添加有向边矩阵,在所述第一有向无环图中添加输出节点至输入节点的第二有向边,得到第二有向无环图;
基于所述第二有向无环图、所述第二有向边构成的第一目标列表和所述第一指令列表,对所述第一量子电路进行等效编译,得到与所述第一量子电路等效的第二量子电路的第二指令列表。
基于可连接性矩阵和第一邻接矩阵更新有向无环图的具体过程如下:
输入:可连接性矩阵all_candidate,有向二分图的第一邻接矩阵bipartite_adjacency,第一有向无环图;
输出:添加有向边矩阵added_edge,有向边列表added_edge_list(即第一目标列表),更新的有向无环图new_graph(即第二有向无环图)。
步骤1:利用python CVXOPT或其他整数规划问题求解器,在满足约束条件的情况下求解优化目标下的添加有向边矩阵added_edge;
步骤2:利用添加有向边矩阵added_edge,得到相应的有向边列表added_edge_list(即第二有向边构成的第一目标列表);
步骤3:基于有向边列表added_edge_list更新第一有向无环图。
之后,可以基于第二有向无环图、第二有向边构成的第一目标列表和第一指令列表,对第一量子电路进行等效编译。如此,可以基于求解得到的添加有向边矩阵实现对量子电路的等效编译。
可选的,所述基于所述第二有向无环图、所述第二有向边构成的第一目标列表和所述第一指令列表,对所述第一量子电路进行等效编译,得到与所述第一量子电路等效的第二量子电路的第二指令列表,包括:
基于所述第二有向无环图,对所述第一指令列表中的指令进行重排序,得到第三指令列表;
基于所述第一目标列表,对所述第三指令列表中的指令进行等效编译,得到所述第二指令列表。
本实施方式中,可以基于第二有向无环图,对第一量子电路的第一指令列表中指令进行重排序,得到第三指令列表。其中,第三指令列表中指令的排列顺序与第二有向无环图中节点的拓扑结构匹配。
第三指令列表中指令的排列顺序与第二有向无环图中节点的拓扑结构匹配,指的是若第二有向无环图中有一条节点A指向节点B的有向边,则第三指令列表中节点A对应指令排序在节点B对应指令之前。
在一可选实施方式中,可以通过指令之间的交换方式来基于第二有向无环图,对第一指令列表中的指令进行重排序,如进行图遍历,根据遍历到的节点和有向边,交换第一指令列表中相关的两个指令,使之与图的拓扑结构匹配。
在另一可选实施方式中,可以获取第二有向无环图对应的拓扑排序列表;基于该拓扑排序列表,对第一指令列表中的指令进行重排序。
之后,可以基于第一目标列表,对第三指令列表中的指令进行等效编译,得到第二指令列表。其中,等效编译的目的是将输入节点对应指令所作用的量子比特重新作用到输出节点测量后被重新利用的量子比特上,这样可以减少所编译得到的量子电路的量子比特数,如此实现量子电路的优化。
可选的,所述基于所述第二有向无环图,对所述第一指令列表中的指令进行重排序,得到第三指令列表,包括:
获取所述第二有向无环图对应的拓扑排序列表,所述拓扑排序列表中,针对所述第二有向无环图中的每个有向边,所述有向边的源节点排序在所述有向边的目标节点之前;
基于所述拓扑排序列表,对所述第一指令列表中的指令进行重排序,得到所述第三指令列表。
本实施方式中,拓扑排序列表指的是表征第二有向无环图拓扑结构的列表,比如,图中有一条节点A指向节点B的有向边,则在拓扑排序列表中节点A将会排在节点B的前面。可以通过遍历第二有向无环图的方式获取该拓扑排序列表。
之后,可以基于该拓扑排序列表,对第一指令列表中的指令进行重排序,比如,可以初始化一个空的指令列表,遍历该拓扑排序列表,按照该拓扑排序列表的排列顺序,将第一指令列表中该拓扑排序列表中节点对应指令排列到该指令列表,得到第三指令列表。
如此,可以实现对第一指令列表的重排序。
可选的,所述基于所述第一目标列表,对所述第三指令列表中的指令进行等效编译,得到所述第二指令列表,包括:
针对所述第一目标列表中每个所述第二有向边,对所述第三指令列表中各第二目标指令所作用的第一量子位更新为第二量子位,所述第一量子位为所述第二有向边中目标节点对应指令所作用的量子位,所述第二量子位为所述第二有向边中源节点对应指令所作用的量子位。
本实施方式中,第二有向边中,源节点指的是输出节点,目标节点指的是输入节点。可以将输入节点对应指令所作用的量子比特(即第一量子位对应量子比特)重新作用到输出节点测量后被重置的量子比特(即第二量子位对应量子比特)上,这样可以减少所编译得到的量子电路的量子比特数,如此实现量子电路的优化。
第二有向无环图至量子电路的等效编译过程如下:
输入:量子电路指令列表circuit_list(即第一指令列表),有向无环图graph(即第二有向无环图),有向边列表edges(即第一目标列表);
输出:编译后的量子电路指令列表(即第二指令列表)。
步骤1:根据有向无环图graph,获得其对应的拓扑排序列表,并记为sorted_gate_ids;
步骤2:根据sorted_gate_ids的顺序,对量子电路指令列表circuit_list进行排序;
步骤3:循环遍历有向边列表edges,记当前被循环的元素为edge,记edge的源节点对应电路指令所作用的量子位为new_idx,记edge的目标节点对应电路指令所作用的量子位为old_idx;循环遍历circuit_list,对所有作用在量子位old_idx上的指令,改写为作用在new_idx上;
步骤4:返回更新后的量子电路指令列表circuit_list作为输出。
基于如图3所示的一个静态量子电路到动态量子电路编译实例,可以获得编译后的动态量子电路为如图5所示,
对应的电路指令列表为:[[reset,0,[[1],[0]],None];[reset,1,[[1],[0]],None];[h,0,None,None];[h,1,None,None];[cnot,[0,1],None,None];[measure,0,[0,‘YZ’,[],[]],(0,1)];[reset,0,[[1],[0]],None];[h,0,None,None];[cnot,[1,0],None,None];[measure,1,[0,‘YZ’,[],[]],(1,1)];[measure,0,[0,‘YZ’,[],[]],(2,1)]]。可知,图5相对于图3,量子电路优化后,其所需要的量子比特的数量减少了。
对于一个量子电路,其完整编译过程如图6所示,具体如下:
输入:量子电路指令列表circuit_list;
输出:编译后的动态量子电路指令列表。
步骤1:以circuit_list为输入,获得第一有向无环图graph,输入节点列表inputs,输出节点列表outputs;
步骤2:以graph,inputs,outputs为输入,获得可连接性矩阵all_candidate;
步骤3:以输入节点列表inputs,输出节点列表outputs,可连接性矩阵all_candidate为输入,获得第一有向无环图的等效有向二分图的第一邻接矩阵bipartite_adjacency;
步骤4:以可连接性矩阵all_candidate,等效有向二分图的第一邻接矩阵bipartite_adjacency为输入,确定添加有向边矩阵added_edge,有向边列表edges,和更新后的有向无环图new_graph(即第二有向无环图);
步骤5:以circuit_list,new_graph,有向边列表edges为输入进行量子电路的等效编译,得到编译后的动态量子电路的指令列表作为输出。
第二实施例
如图7所示,本公开提供一种量子电路处理装置700,包括:
获取模块701,用于获取第一量子电路的第一指令列表;
第一确定模块702,用于基于所述第一指令列表,确定第一有向无环图,所述第一有向无环图包括所述第一指令列表中指令对应的节点和至少两个第一有向边,所述第一有向边用于表征所述第一指令列表中不同指令间的时序关系;
第二确定模块703,用于基于所述第一有向无环图,确定所述第一有向无环图等效的有向二分图的第一邻接矩阵,所述第一邻接矩阵用于表征所述有向二分图的节点集合中每两个节点之间是否存在有向通路,所述节点集合包括所述第一有向无环图中的输入节点和输出节点;
第三确定模块704,用于基于所述第一邻接矩阵,确定在所述有向二分图中添加输出节点至输入节点的有向边数量达到极大值的优化目标下的添加有向边矩阵,基于所述添加有向边矩阵在所述有向二分图添加有向边的情况下,添加有向边后的所述有向二分图中不存在有向环路,且所述添加有向边矩阵满足:每个输出节点最多只能和一个输入节点相连,以及每个输入节点最多只能和一个输出节点相连;
等效编译模块705,用于基于所述添加有向边矩阵、所述第一有向无环图和所述第一指令列表,对所述第一量子电路进行等效编译,得到与所述第一量子电路等效的第二量子电路的第二指令列表。
可选的,所述第一邻接矩阵按照所述第一有向无环图中的输入节点数量进行分块,所述第一邻接矩阵为2×2的分块矩阵,所述第二确定模块703,具体用于:
基于所述第一有向无环图,确定所述第一有向无环图的可达性矩阵,所述可达性矩阵用于表征所述第一有向无环图中输入节点到输出节点的可达性,所述可达性指示输入节点是否存在到达输出节点的有向通路;
将所述可达性矩阵作为所述第一邻接矩阵中第一行第二列的块矩阵,得到所述第一邻接矩阵。
可选的,所述第一邻接矩阵为反对角分块矩阵,所述第一邻接矩阵中除所述可达性矩阵之外的块矩阵均为全零矩阵。
可选的,所述第三确定模块704,具体用于:
基于所述添加有向边矩阵,构建优化目标的第一函数,所述第一函数表征所述添加有向边矩阵中元素值的和值达到极值,所述极值表征在所述有向二分图中添加输出节点至输入节点的有向边数量达到极大值;
基于所述第一邻接矩阵,构建所述添加有向边矩阵的约束条件;
基于所述约束条件和所述第一函数,求解所述优化目标下的添加有向边矩阵。
可选的,所述约束条件包括:
所述添加有向边矩阵为所述第一有向无环图的可连接性矩阵转置的子矩阵,所述可连接性矩阵用于表征所述第一有向无环图中输出节点至输入节点的可连接性;
添加有向边后的所述有向二分图的第二邻接矩阵满足幂零;
所述添加有向边矩阵的每列元素值中最多存在一个第一目标值;
所述添加有向边矩阵的每行元素值中最多存在一个第一目标值,所述第一目标值表征在所述有向二分图中添加输出节点至输入节点的有向边。
可选的,使添加有向边后的所述有向二分图的第二邻接矩阵满足幂零的幂为4。
可选的,在所述第一目标值为1的情况下,所述约束条件包括:
BY′BY′=0;
Y′BY′B=0;
其中,Yij为所述可连接性矩阵中索引为(j,i)的元素,B为所述第一有向无环图的可达性矩阵,所述可达性矩阵中的元素值与所述可连接性矩阵中的元素值互补,Y′为添加有向边矩阵,Y′ij为添加有向边矩阵中索引为(i,j)的元素,n为所述第一有向无环图中的输入节点数量。
可选的,在所述第一目标值为1的情况下,所述第一函数表征所述添加有向边矩阵中元素值的和值达到极大值。
可选的,所述等效编译模块705包括:
更新子模块,用于基于所述添加有向边矩阵,在所述第一有向无环图中添加输出节点至输入节点的第二有向边,得到第二有向无环图;
等效编译子模块,用于基于所述第二有向无环图、所述第二有向边构成的第一目标列表和所述第一指令列表,对所述第一量子电路进行等效编译,得到与所述第一量子电路等效的第二量子电路的第二指令列表。
可选的,所述等效编译子模块包括:
重排序单元,用于基于所述第二有向无环图,对所述第一指令列表中的指令进行重排序,得到第三指令列表;
等效编译单元,用于基于所述第一目标列表,对所述第三指令列表中的指令进行等效编译,得到所述第二指令列表。
可选的,所述重排序单元,具体用于:
获取所述第二有向无环图对应的拓扑排序列表,所述拓扑排序列表中,针对所述第二有向无环图中的每个有向边,所述有向边的源节点排序在所述有向边的目标节点之前;
基于所述拓扑排序列表,对所述第一指令列表中的指令进行重排序,得到所述第三指令列表。
可选的,所述等效编译单元,具体用于:
针对所述第一目标列表中每个所述第二有向边,对所述第三指令列表中各第二目标指令所作用的第一量子位更新为第二量子位,所述第一量子位为所述第二有向边中目标节点对应指令所作用的量子位,所述第二量子位为所述第二有向边中源节点对应指令所作用的量子位。
可选的,所述第一确定模块702,具体用于:
对所述第一指令列表按照从左至右的指令排列顺序进行遍历,并获取第二目标列表和当前遍历的指令所作用的第三量子位,所述第二目标列表包括所述第一指令列表中排序在所述当前遍历的指令之前的指令;
对所述第二目标列表按照从右至左的指令排列顺序进行遍历,并获取第一目标指令,所述第一目标指令所作用的量子位与所述第三量子位存在交集;
将所述第一目标指令的标识作为源节点,所述当前遍历的指令的标识作为目标节点,构建所述第一有向边;
在所述第一指令列表遍历完成的情况下,基于所构建的第一有向边,获取所述第一有向无环图。
本公开提供的量子电路处理装置700能够实现量子电路处理方法实施例实现的各个过程,且能够达到相同的有益效果,为避免重复,这里不再赘述。
本公开的技术方案中,所涉及的用户个人信息的收集、存储、使用、加工、传输、提供和公开等处理,均符合相关法律法规的规定,且不违背公序良俗。
根据本公开的实施例,本公开还提供了一种电子设备、一种可读存储介质和一种计算机程序产品。
图8示出了可以用来实施本公开的实施例的示例电子设备的示意性框图。电子设备旨在表示各种形式的数字计算机,诸如,膝上型计算机、台式计算机、工作台、个人数字助理、服务器、刀片式服务器、大型计算机、和其它适合的计算机。电子设备还可以表示各种形式的移动装置,诸如,个人数字处理、蜂窝电话、智能电话、可穿戴设备和其它类似的计算装置。本文所示的部件、它们的连接和关系、以及它们的功能仅仅作为示例,并且不意在限制本文中描述的和/或者要求的本公开的实现。
如图8所示,设备800包括计算单元801,其可以根据存储在只读存储器(ROM)802中的计算机程序或者从存储单元808加载到随机访问存储器(RAM)803中的计算机程序,来执行各种适当的动作和处理。在RAM 803中,还可存储设备800操作所需的各种程序和数据。计算单元801、ROM 802以及RAM 803通过总线804彼此相连。输入/输出(I/O)接口805也连接至总线804。
设备800中的多个部件连接至I/O接口805,包括:输入单元806,例如键盘、鼠标等;输出单元807,例如各种类型的显示器、扬声器等;存储单元808,例如磁盘、光盘等;以及通信单元809,例如网卡、调制解调器、无线通信收发机等。通信单元809允许设备800通过诸如因特网的计算机网络和/或各种电信网络与其他设备交换信息/数据。
计算单元801可以是各种具有处理和计算能力的通用和/或专用处理组件。计算单元801的一些示例包括但不限于中央处理单元(CPU)、图形处理单元(GPU)、各种专用的人工智能(AI)计算芯片、各种运行机器学习模型算法的计算单元、数字信号处理器(DSP)、以及任何适当的处理器、控制器、微控制器等。计算单元801执行上文所描述的各个方法和处理,例如量子电路处理方法。例如,在一些实施例中,量子电路处理方法可被实现为计算机软件程序,其被有形地包含于机器可读介质,例如存储单元808。在一些实施例中,计算机程序的部分或者全部可以经由ROM 802和/或通信单元809而被载入和/或安装到设备800上。当计算机程序加载到RAM 803并由计算单元801执行时,可以执行上文描述的量子电路处理方法的一个或多个步骤。备选地,在其他实施例中,计算单元801可以通过其他任何适当的方式(例如,借助于固件)而被配置为执行量子电路处理方法。
本文中以上描述的系统和技术的各种实施方式可以在数字电子电路系统、集成电路系统、场可编程门阵列(FPGA)、专用集成电路(ASIC)、专用标准产品(ASSP)、芯片上系统的系统(SOC)、负载可编程逻辑设备(CPLD)、计算机硬件、固件、软件、和/或它们的组合中实现。这些各种实施方式可以包括:实施在一个或者多个计算机程序中,该一个或者多个计算机程序可在包括至少一个可编程处理器的可编程系统上执行和/或解释,该可编程处理器可以是专用或者通用可编程处理器,可以从存储系统、至少一个输入装置、和至少一个输出装置接收数据和指令,并且将数据和指令传输至该存储系统、该至少一个输入装置、和该至少一个输出装置。
用于实施本公开的方法的程序代码可以采用一个或多个编程语言的任何组合来编写。这些程序代码可以提供给通用计算机、专用计算机或其他可编程数据处理装置的处理器或控制器,使得程序代码当由处理器或控制器执行时使流程图和/或框图中所规定的功能/操作被实施。程序代码可以完全在机器上执行、部分地在机器上执行,作为独立软件包部分地在机器上执行且部分地在远程机器上执行或完全在远程机器或服务器上执行。
在本公开的上下文中,机器可读介质可以是有形的介质,其可以包含或存储以供指令执行系统、装置或设备使用或与指令执行系统、装置或设备结合地使用的程序。机器可读介质可以是机器可读信号介质或机器可读储存介质。机器可读介质可以包括但不限于电子的、磁性的、光学的、电磁的、红外的、或半导体系统、装置或设备,或者上述内容的任何合适组合。机器可读存储介质的更具体示例会包括基于一个或多个线的电气连接、便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦除可编程只读存储器(EPROM或快闪存储器)、光纤、便捷式紧凑盘只读存储器(CD-ROM)、光学储存设备、磁储存设备、或上述内容的任何合适组合。
为了提供与用户的交互,可以在计算机上实施此处描述的系统和技术,该计算机具有:用于向用户显示信息的显示装置(例如,CRT(阴极射线管)或者LCD(液晶显示器)监视器);以及键盘和指向装置(例如,鼠标或者轨迹球),用户可以通过该键盘和该指向装置来将输入提供给计算机。其它种类的装置还可以用于提供与用户的交互;例如,提供给用户的反馈可以是任何形式的传感反馈(例如,视觉反馈、听觉反馈、或者触觉反馈);并且可以用任何形式(包括声输入、语音输入或者、触觉输入)来接收来自用户的输入。
可以将此处描述的系统和技术实施在包括后台部件的计算系统(例如,作为数据服务器)、或者包括中间件部件的计算系统(例如,应用服务器)、或者包括前端部件的计算系统(例如,具有图形用户界面或者网络浏览器的用户计算机,用户可以通过该图形用户界面或者该网络浏览器来与此处描述的系统和技术的实施方式交互)、或者包括这种后台部件、中间件部件、或者前端部件的任何组合的计算系统中。可以通过任何形式或者介质的数字数据通信(例如,通信网络)来将系统的部件相互连接。通信网络的示例包括:局域网(LAN)、广域网(WAN)和互联网。
计算机系统可以包括客户端和服务器。客户端和服务器一般远离彼此并且通常通过通信网络进行交互。通过在相应的计算机上运行并且彼此具有客户端-服务器关系的计算机程序来产生客户端和服务器的关系。服务器可以是云服务器,也可以为分布式系统的服务器,或者是结合了区块链的服务器。
应该理解,可以使用上面所示的各种形式的流程,重新排序、增加或删除步骤。例如,本公开中记载的各步骤可以并行地执行也可以顺序地执行也可以不同的次序执行,只要能够实现本公开公开的技术方案所期望的结果,本文在此不进行限制。
上述具体实施方式,并不构成对本公开保护范围的限制。本领域技术人员应该明白的是,根据设计要求和其他因素,可以进行各种修改、组合、子组合和替代。任何在本公开的精神和原则之内所作的修改、等同替换和改进等,均应包含在本公开保护范围之内。
Claims (28)
1.一种量子电路处理方法,包括:
获取第一量子电路的第一指令列表;
基于所述第一指令列表,确定第一有向无环图,所述第一有向无环图包括所述第一指令列表中指令对应的节点和至少两个第一有向边,所述第一有向边用于表征所述第一指令列表中不同指令间的时序关系;
基于所述第一有向无环图,确定所述第一有向无环图等效的有向二分图的第一邻接矩阵,所述第一邻接矩阵用于表征所述有向二分图的节点集合中每两个节点之间是否存在有向通路,所述节点集合包括所述第一有向无环图中的输入节点和输出节点;
基于所述第一邻接矩阵,确定在所述有向二分图中添加输出节点至输入节点的有向边数量达到极大值的优化目标下的添加有向边矩阵,基于所述添加有向边矩阵在所述有向二分图添加有向边的情况下,添加有向边后的所述有向二分图中不存在有向环路,且所述添加有向边矩阵满足:每个输出节点最多只能和一个输入节点相连,以及每个输入节点最多只能和一个输出节点相连;
基于所述添加有向边矩阵、所述第一有向无环图和所述第一指令列表,对所述第一量子电路进行等效编译,得到与所述第一量子电路等效的第二量子电路的第二指令列表。
2.根据权利要求1所述的方法,其中,所述第一邻接矩阵按照所述第一有向无环图中的输入节点数量进行分块,所述第一邻接矩阵为2×2的分块矩阵,所述基于所述第一有向无环图,确定所述第一有向无环图等效的有向二分图的第一邻接矩阵,包括:
基于所述第一有向无环图,确定所述第一有向无环图的可达性矩阵,所述可达性矩阵用于表征所述第一有向无环图中输入节点到输出节点的可达性,所述可达性指示输入节点是否存在到达输出节点的有向通路;
将所述可达性矩阵作为所述第一邻接矩阵中第一行第二列的块矩阵,得到所述第一邻接矩阵。
3.根据权利要求2所述的方法,其中,所述第一邻接矩阵为反对角分块矩阵,所述第一邻接矩阵中除所述可达性矩阵之外的块矩阵均为全零矩阵。
4.根据权利要求1所述的方法,其中,所述基于所述第一邻接矩阵,确定在所述有向二分图中添加输出节点至输入节点的有向边数量达到极大值的优化目标下的添加有向边矩阵,包括:
基于所述添加有向边矩阵,构建优化目标的第一函数,所述第一函数表征所述添加有向边矩阵中元素值的和值达到极值,所述极值表征在所述有向二分图中添加输出节点至输入节点的有向边数量达到极大值;
基于所述第一邻接矩阵,构建所述添加有向边矩阵的约束条件;
基于所述约束条件和所述第一函数,求解所述优化目标下的添加有向边矩阵。
5.根据权利要求4所述的方法,其中,所述约束条件包括:
所述添加有向边矩阵为所述第一有向无环图的可连接性矩阵转置的子矩阵,所述可连接性矩阵用于表征所述第一有向无环图中输出节点至输入节点的可连接性;
添加有向边后的所述有向二分图的第二邻接矩阵满足幂零;
所述添加有向边矩阵的每列元素值中最多存在一个第一目标值;
所述添加有向边矩阵的每行元素值中最多存在一个第一目标值,所述第一目标值表征在所述有向二分图中添加输出节点至输入节点的有向边。
6.根据权利要求5所述的方法,其中,使添加有向边后的所述有向二分图的第二邻接矩阵满足幂零的幂为4。
7.根据权利要求5所述的方法,其中,在所述第一目标值为1的情况下,所述约束条件包括:
BY′BY′=0;
Y′BY′B=0;
其中,Yij为所述可连接性矩阵中索引为(j,i)的元素,B为所述第一有向无环图的可达性矩阵,所述可达性矩阵中的元素值与所述可连接性矩阵中的元素值互补,Y′为添加有向边矩阵,Y′ij为添加有向边矩阵中索引为(i,j)的元素,n为所述第一有向无环图中的输入节点数量。
8.根据权利要求5所述的方法,其中,在所述第一目标值为1的情况下,所述第一函数表征所述添加有向边矩阵中元素值的和值达到极大值。
9.根据权利要求1所述的方法,其中,所述基于所述添加有向边矩阵、所述第一有向无环图和所述第一指令列表,对所述第一量子电路进行等效编译,得到与所述第一量子电路等效的第二量子电路的第二指令列表,包括:
基于所述添加有向边矩阵,在所述第一有向无环图中添加输出节点至输入节点的第二有向边,得到第二有向无环图;
基于所述第二有向无环图、所述第二有向边构成的第一目标列表和所述第一指令列表,对所述第一量子电路进行等效编译,得到与所述第一量子电路等效的第二量子电路的第二指令列表。
10.根据权利要求9所述的方法,其中,所述基于所述第二有向无环图、所述第二有向边构成的第一目标列表和所述第一指令列表,对所述第一量子电路进行等效编译,得到与所述第一量子电路等效的第二量子电路的第二指令列表,包括:
基于所述第二有向无环图,对所述第一指令列表中的指令进行重排序,得到第三指令列表;
基于所述第一目标列表,对所述第三指令列表中的指令进行等效编译,得到所述第二指令列表。
11.根据权利要求10所述的方法,其中,所述基于所述第二有向无环图,对所述第一指令列表中的指令进行重排序,得到第三指令列表,包括:
获取所述第二有向无环图对应的拓扑排序列表,所述拓扑排序列表中,针对所述第二有向无环图中的每个有向边,所述有向边的源节点排序在所述有向边的目标节点之前;
基于所述拓扑排序列表,对所述第一指令列表中的指令进行重排序,得到所述第三指令列表。
12.根据权利要求10所述的方法,其中,所述基于所述第一目标列表,对所述第三指令列表中的指令进行等效编译,得到所述第二指令列表,包括:
针对所述第一目标列表中每个所述第二有向边,对所述第三指令列表中各第二目标指令所作用的第一量子位更新为第二量子位,所述第一量子位为所述第二有向边中目标节点对应指令所作用的量子位,所述第二量子位为所述第二有向边中源节点对应指令所作用的量子位。
13.根据权利要求1所述的方法,其中,所述基于所述第一指令列表,确定第一有向无环图,包括:
对所述第一指令列表按照从左至右的指令排列顺序进行遍历,并获取第二目标列表和当前遍历的指令所作用的第三量子位,所述第二目标列表包括所述第一指令列表中排序在所述当前遍历的指令之前的指令;
对所述第二目标列表按照从右至左的指令排列顺序进行遍历,并获取第一目标指令,所述第一目标指令所作用的量子位与所述第三量子位存在交集;
将所述第一目标指令的标识作为源节点,所述当前遍历的指令的标识作为目标节点,构建所述第一有向边;
在所述第一指令列表遍历完成的情况下,基于所构建的第一有向边,获取所述第一有向无环图。
14.一种量子电路处理装置,包括:
获取模块,用于获取第一量子电路的第一指令列表;
第一确定模块,用于基于所述第一指令列表,确定第一有向无环图,所述第一有向无环图包括所述第一指令列表中指令对应的节点和至少两个第一有向边,所述第一有向边用于表征所述第一指令列表中不同指令间的时序关系;
第二确定模块,用于基于所述第一有向无环图,确定所述第一有向无环图等效的有向二分图的第一邻接矩阵,所述第一邻接矩阵用于表征所述有向二分图的节点集合中每两个节点之间是否存在有向通路,所述节点集合包括所述第一有向无环图中的输入节点和输出节点;
第三确定模块,用于基于所述第一邻接矩阵,确定在所述有向二分图中添加输出节点至输入节点的有向边数量达到极大值的优化目标下的添加有向边矩阵,基于所述添加有向边矩阵在所述有向二分图添加有向边的情况下,添加有向边后的所述有向二分图中不存在有向环路,且所述添加有向边矩阵满足:每个输出节点最多只能和一个输入节点相连,以及每个输入节点最多只能和一个输出节点相连;
等效编译模块,用于基于所述添加有向边矩阵、所述第一有向无环图和所述第一指令列表,对所述第一量子电路进行等效编译,得到与所述第一量子电路等效的第二量子电路的第二指令列表。
15.根据权利要求14所述的装置,其中,所述第一邻接矩阵按照所述第一有向无环图中的输入节点数量进行分块,所述第一邻接矩阵为2×2的分块矩阵,所述第二确定模块,具体用于:
基于所述第一有向无环图,确定所述第一有向无环图的可达性矩阵,所述可达性矩阵用于表征所述第一有向无环图中输入节点到输出节点的可达性,所述可达性指示输入节点是否存在到达输出节点的有向通路;
将所述可达性矩阵作为所述第一邻接矩阵中第一行第二列的块矩阵,得到所述第一邻接矩阵。
16.根据权利要求15所述的装置,其中,所述第一邻接矩阵为反对角分块矩阵,所述第一邻接矩阵中除所述可达性矩阵之外的块矩阵均为全零矩阵。
17.根据权利要求14所述的装置,其中,所述第三确定模块,具体用于:
基于所述添加有向边矩阵,构建优化目标的第一函数,所述第一函数表征所述添加有向边矩阵中元素值的和值达到极值,所述极值表征在所述有向二分图中添加输出节点至输入节点的有向边数量达到极大值;
基于所述第一邻接矩阵,构建所述添加有向边矩阵的约束条件;
基于所述约束条件和所述第一函数,求解所述优化目标下的添加有向边矩阵。
18.根据权利要求17所述的装置,其中,所述约束条件包括:
所述添加有向边矩阵为所述第一有向无环图的可连接性矩阵转置的子矩阵,所述可连接性矩阵用于表征所述第一有向无环图中输出节点至输入节点的可连接性;
添加有向边后的所述有向二分图的第二邻接矩阵满足幂零;
所述添加有向边矩阵的每列元素值中最多存在一个第一目标值;
所述添加有向边矩阵的每行元素值中最多存在一个第一目标值,所述第一目标值表征在所述有向二分图中添加输出节点至输入节点的有向边。
19.根据权利要求18所述的装置,其中,使添加有向边后的所述有向二分图的第二邻接矩阵满足幂零的幂为4。
20.根据权利要求18所述的装置,其中,在所述第一目标值为1的情况下,所述约束条件包括:
BY′BY′=0;
Y′BY′B=0;
其中,Yij为所述可连接性矩阵中索引为(j,i)的元素,B为所述第一有向无环图的可达性矩阵,所述可达性矩阵中的元素值与所述可连接性矩阵中的元素值互补,Y′为添加有向边矩阵,Y′ij为添加有向边矩阵中索引为(i,j)的元素,n为所述第一有向无环图中的输入节点数量。
21.根据权利要求18所述的装置,其中,在所述第一目标值为1的情况下,所述第一函数表征所述添加有向边矩阵中元素值的和值达到极大值。
22.根据权利要求14所述的装置,其中,所述等效编译模块包括:
更新子模块,用于基于所述添加有向边矩阵,在所述第一有向无环图中添加输出节点至输入节点的第二有向边,得到第二有向无环图;
等效编译子模块,用于基于所述第二有向无环图、所述第二有向边构成的第一目标列表和所述第一指令列表,对所述第一量子电路进行等效编译,得到与所述第一量子电路等效的第二量子电路的第二指令列表。
23.根据权利要求22所述的装置,其中,所述等效编译子模块包括:
重排序单元,用于基于所述第二有向无环图,对所述第一指令列表中的指令进行重排序,得到第三指令列表;
等效编译单元,用于基于所述第一目标列表,对所述第三指令列表中的指令进行等效编译,得到所述第二指令列表。
24.根据权利要求23所述的装置,其中,所述重排序单元,具体用于:
获取所述第二有向无环图对应的拓扑排序列表,所述拓扑排序列表中,针对所述第二有向无环图中的每个有向边,所述有向边的源节点排序在所述有向边的目标节点之前;
基于所述拓扑排序列表,对所述第一指令列表中的指令进行重排序,得到所述第三指令列表。
25.根据权利要求23所述的装置,其中,所述等效编译单元,具体用于:
针对所述第一目标列表中每个所述第二有向边,对所述第三指令列表中各第二目标指令所作用的第一量子位更新为第二量子位,所述第一量子位为所述第二有向边中目标节点对应指令所作用的量子位,所述第二量子位为所述第二有向边中源节点对应指令所作用的量子位。
26.根据权利要求14所述的装置,其中,所述第一确定模块,具体用于:
对所述第一指令列表按照从左至右的指令排列顺序进行遍历,并获取第二目标列表和当前遍历的指令所作用的第三量子位,所述第二目标列表包括所述第一指令列表中排序在所述当前遍历的指令之前的指令;
对所述第二目标列表按照从右至左的指令排列顺序进行遍历,并获取第一目标指令,所述第一目标指令所作用的量子位与所述第三量子位存在交集;
将所述第一目标指令的标识作为源节点,所述当前遍历的指令的标识作为目标节点,构建所述第一有向边;
在所述第一指令列表遍历完成的情况下,基于所构建的第一有向边,获取所述第一有向无环图。
27.一种电子设备,包括:
至少一个处理器;以及
与所述至少一个处理器通信连接的存储器;其中,
所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够执行权利要求1-13中任一项所述的方法。
28.一种存储有计算机指令的非瞬时计算机可读存储介质,其中,所述计算机指令用于使所述计算机执行根据权利要求1-13中任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310761796.1A CN116629370B (zh) | 2023-06-26 | 2023-06-26 | 量子电路处理方法、装置及电子设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310761796.1A CN116629370B (zh) | 2023-06-26 | 2023-06-26 | 量子电路处理方法、装置及电子设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116629370A CN116629370A (zh) | 2023-08-22 |
CN116629370B true CN116629370B (zh) | 2024-05-24 |
Family
ID=87617222
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310761796.1A Active CN116629370B (zh) | 2023-06-26 | 2023-06-26 | 量子电路处理方法、装置及电子设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116629370B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117313883B (zh) * | 2023-09-27 | 2025-01-28 | 北京百度网讯科技有限公司 | 量子电路处理方法、装置及电子设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111178532A (zh) * | 2019-12-11 | 2020-05-19 | 合肥本源量子计算科技有限责任公司 | 一种量子线路匹配的方法、装置、存储介质和电子装置 |
WO2020117552A1 (en) * | 2018-12-07 | 2020-06-11 | IonQ, Inc. | Methods and apparatuses for two-qubit gate reduction in quantum circuits |
CN116151381A (zh) * | 2023-02-20 | 2023-05-23 | 北京百度网讯科技有限公司 | 量子电路处理方法、装置及电子设备 |
CN116151384A (zh) * | 2023-02-20 | 2023-05-23 | 北京百度网讯科技有限公司 | 量子电路处理方法、装置及电子设备 |
-
2023
- 2023-06-26 CN CN202310761796.1A patent/CN116629370B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020117552A1 (en) * | 2018-12-07 | 2020-06-11 | IonQ, Inc. | Methods and apparatuses for two-qubit gate reduction in quantum circuits |
CN111178532A (zh) * | 2019-12-11 | 2020-05-19 | 合肥本源量子计算科技有限责任公司 | 一种量子线路匹配的方法、装置、存储介质和电子装置 |
CN116151381A (zh) * | 2023-02-20 | 2023-05-23 | 北京百度网讯科技有限公司 | 量子电路处理方法、装置及电子设备 |
CN116151384A (zh) * | 2023-02-20 | 2023-05-23 | 北京百度网讯科技有限公司 | 量子电路处理方法、装置及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN116629370A (zh) | 2023-08-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116151384B (zh) | 量子电路处理方法、装置及电子设备 | |
US20230267358A1 (en) | Distributed Quantum Computing Simulation Method and Apparatus | |
CN116151381B (zh) | 量子电路处理方法、装置及电子设备 | |
US11010527B2 (en) | Optimization of a quantum circuit by inserting swap gates | |
CN114881237B (zh) | 量子计算处理方法、装置及电子设备 | |
CN116611527B (zh) | 量子电路处理方法、装置及电子设备 | |
CN116187458B (zh) | 量子电路处理方法、装置及电子设备 | |
CN114970865B (zh) | 量子芯片上的量子电路处理方法、装置及电子设备 | |
CN116167446B (zh) | 量子计算处理方法、装置及电子设备 | |
CN114580645B (zh) | 随机量子测量的模拟方法、装置、设备及存储介质 | |
CN116629370B (zh) | 量子电路处理方法、装置及电子设备 | |
CN116579435B (zh) | 量子电路的分类方法、装置、电子设备、介质和产品 | |
CN115577776B (zh) | 基态能量的确定方法、装置、设备及存储介质 | |
CN109299725B (zh) | 一种张量链并行实现高阶主特征值分解的预测系统和装置 | |
CN115577782B (zh) | 量子计算方法、装置、设备及存储介质 | |
CN115458040B (zh) | 蛋白质的生成方法、装置、电子设备及存储介质 | |
CN117313881B (zh) | 量子电路的分类方法、装置及电子设备 | |
CN117313878B (zh) | 量子电路处理方法、装置及电子设备 | |
CN117313883B (zh) | 量子电路处理方法、装置及电子设备 | |
CN116167447B (zh) | 量子电路处理方法、装置及电子设备 | |
CN117313882A (zh) | 量子电路处理方法、装置及电子设备 | |
CN117313877A (zh) | 量子电路处理方法、装置及电子设备 | |
CN117313884A (zh) | 量子电路处理方法、装置及电子设备 | |
US20240135216A1 (en) | Solving optimization problems on shallow circuits using a quantum computer | |
CN117787424A (zh) | 量子电路处理方法、装置及电子设备 |
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 |