量子计算机怎么编程,量子计算机编程(一)—

今天要给大家介绍的是这本书《Programming Quantum Computers —— Essential Algorithms and Code Samples》，主要讲如何在量子计算机上编程，又示例的代码可以跑 https://oreilly-qc.github.io./ ，这本书的理解方向是纯粹计算机的，连矩阵都很少涉及，一个量子计算机就像是协处理器一样，很像我们现在的GPU。

这本书一共分为三个部分，我也将分成三篇文章来写。

单量子比特

这本书中对于量子的表达方式如下circle notation：

蓝色圈的面积就是这个量子态测量后是这个状态的概率，中间的蓝色短线代表他们的相位，因为全局相位也观察不出来，所以一般$|0\rangle$就是朝上的。对于一个量子比特来说，重要的也就这两个，magnitude和relative phase(可能因为重要的是概率大小，所以他其实没有提amplitude其实是一个复数)。

这两个态实际上是等价的。

电路图呢，一般长这样

这个电路图对应的代码如下：(这个例子都可以在上面那个网站上跑起来)

qc.reset(1); // allocate one qubit

qc.write(0); // write the value zero

qc.had(); // place it into superposition of 0 and 1

var result = qc.read(); // read the result as a digital bit

第一步就是申请一个qubit，就像你要给变量分配空间一样。

第二步写0，其实写0很容易，你可以直接测量，要么0，要么1，如果结果是1的话，再做一个not操作就好；当然，你要是嫌麻烦，对于一个qubit，长时间的静置他，他也会变成0，毕竟还是基态比较稳定。

第三步就是进行一个H门的操作

第四步读，其实就是测量了

常见的作用于单量子比特操作的表达方式：

其中值得一提的是PHASE相位操作，phase操作只作用在 $|1\rangle$ 上，因为他的效果是改变相对相位，如果大家都改变就没有什么用了。可能有同学听说过绕X轴旋转和绕Y轴旋转，这些都是针对Bloch球的表达方式，与这里的circle notation的方式不要弄混了。

一个操作也可能是其他几个操作的组合，比如：

如果我们可以把not变成H+180°旋转+H，那么我们也可以把中间的180°旋转变成两个90°的旋转，中间再加上两个H，因为$HH=I$，他们可以相互抵消，这样我们还可以得到RNOT：

**COPY：**这是一个需要注意的操作，因为量子程序里没有复制，这也保障了量子传输的信息不会被窃听，因为你要窃听，你就需要去读，一读就是测量，而量子比特一测量就是坍缩。由于量子不能复制，所以，上述所有操作都是在原有的那个比特上操作的，所以操作就会被发现。

一个简单例子：

对于这里例子来说， $A_1$ 和 $A_2$ 是对qubit Hadmard门操作后测量得到的随机值，因为这个是真随机，所以就不会被窃听者提前知道或者预估，那么当我传输这个被我用红色圈出来的比特的时候，spy并不知道这个比特是否执行了H门和not操作，那他就只能猜了，25%的概率，图里面所示就是他猜有执行了，然后他再按照他的猜测如法炮制一个qubit继续传递，当B收了这个bit后，B也随机一个数据$B_2$，看是否执行H门，然后测量。这个时候消息也都收到了，测量也都测量好了，那么把 $A_2$的信息发过来也和窃听没有关系了，如果 $B_2$ $A_2$ 的结果一样，那么测量结果应该一样，如果不一样，那么一定被窃听了。

这样的成功概率有多少呢？ $B_2$ $A_2$ 一样的概率0.5，在这种情况下spy被发现的概率0.25。看起来不是很高，但是如果我们有一百个比特先检测一下这条线路，不被发现的概率将会降到百万分之一。

qc.reset(3);//申请三个qubit

qc.discard();

var a = qint.new(1, 'alice'); //给其中一个变量命名为a,但是画出来的电路图中显示alice

var fiber = qint.new(1, 'fiber');

var b = qint.new(1, 'bob');

function random_bit(q) {//对一个初始化为0的比特，进行H操作，然后测量，测量结果是随机的

q.write(0);

q.had();

return q.read();

}

// Generate two random bits

qc.label('get two random bits');

var send_had = random_bit(a);//得到是否执行H门的随机值

var send_value = random_bit(a);//得到是否not的随机值

qc.label('');

// Prepare Alice's qubit

a.write(0);//a重新赋值为0，所以前面取随机值的操作也可以在a上进行

qc.label('set value');

qc.nop();

if (send_value)

a.not();

qc.nop();

qc.label('');

qc.nop();

qc.label('apply had');

qc.nop();

if (send_had)

a.had();

qc.nop();

qc.label('');

// Send the qubit!

fiber.exchange(a);

// Activate the spy

var spy_is_present = true;

if (spy_is_present)

{

var spy_had = 1;

qc.nop();

qc.label('spy');

if (spy_had)

fiber.had();

stolen_data = fiber.read();

fiber.write(0);

if (stolen_data)

fiber.not();

if (spy_had)

fiber.had();

qc.label('');

qc.nop();

}

// Receive the qubit!

var recv_had = random_bit(b);

fiber.exchange(b);

qc.label('apply had');

qc.nop();

if (recv_had)

b.had();

qc.nop();

qc.label('');

qc.nop();

qc.label('read value');

qc.nop();

recv_val = b.read();

qc.nop();

qc.label('');

qc.nop();

// Now Alice emails Bob to tell

// him her had setting and value.

// If the had setting matches and the

// value does not, there's a spy!

if (send_had == recv_had)

if (send_value != recv_val)

qc.print('Caught a spy!\n');

多量子比特

多量子比特的表示方法和单量子比特没有什么区别，下面的数字就是把他们从01的二进制换算成了我们熟悉的十进制。0x是十六进制的表示方法，最上面的是地位，正好可以not操作表示了出来，0就是000，1就是001，2就是010……以此类推

真要说和单量子比特的区别，那主要是两方面，一个是纠缠；另一个就是受控操作。

纠缠

纠缠好说，两个纠缠的量子比特就是其中一个的测量结果会影响另一个。

受控操作

受控操作就是说有两类比特，一个是用来控制的，另一个是被控制的。

比如大家耳熟能详的CNOT，当控制比特是1的时候，就翻转被控制的比特，如果控制比特是0的时候，那就不做操作。还有CCNOT门，又叫做toffli门，有两个控制比特，只有当他们都为1的时候才翻转。

现在要介绍一个前面没有提过的CPHASE受控相位门，和CNOT相似，同样是当控制比特为1的时候才进行，但是正如前面我们提过，PHASE操作只操作在1上，也就是说，如果这个操作被执行其实只有一种情况，那就是$|11\rangle$ 。

在前面我们表示过一个旋转180°的相位操作可以表示成一个 H+NOT+H figure 2-14

那么现在受控相位操作拥有了三种表示方法：

当控制比特为0的时候，CNOT就不会发生，两个H又正好抵消了。

phase kickback

在这之所以说这么多的受控相位操作是因为这里面有个一很有趣的小技巧 QPU Trick: Phase Kickback 来看下面这个电路图

这个电路图很容易理解，一种理解方式就是把register 1当作控制比特，register 2 当作受控比特，两个H门使得register 1变成了等可能的四种情况 $\frac{1}{2}|00\rangle+\frac{1}{2}|01\rangle+\frac{1}{2}|10\rangle+\frac{1}{2}|11\rangle$ 。

现在来看看结果：