New!! 量子ゲートでピクセルアートを作る（③8✕8サイズ完成編）

#塗るかどうかの総合判定ビット
OR(13, 15, 32)
OR(32, 20, 33)
OR(33, 24, 34)
OR(34, 28, 35)
OR(35, 31, 36)

#サンプリング
qc.measure_all()
backend = AerSimulator(method='matrix_product_state')
result = backend.run(qc, shots=2000).result().get_counts()

#解の統計
pos = (0, 1, 2, 3, 4, 5, 13, 15 ,20, 24, 28, 31, 36)

for i in range(2):
    for j in range(2):
        for k in range(2):
            for l in range(2):
                for m in range(2):
                    for n in range(2):
                        for r in result:
                            # print(r)
                            sel = ''.join(r[::-1][p] for p in pos)
                            if f'{i}{j}{k}{l}{m}{n}' == sel[:6]:
                                print(f'{sel[:6]} | {sel[6:12]} | {sel[12]} | {result[r]}/500')

000000 | 000000 | 0 | 34/500
000001 | 000000 | 0 | 27/500
000010 | 000000 | 0 | 26/500
000011 | 010000 | 1 | 32/500
000100 | 000000 | 0 | 35/500
000101 | 000000 | 0 | 29/500
000110 | 000000 | 0 | 28/500
000111 | 010000 | 1 | 31/500
001000 | 000100 | 1 | 25/500
001001 | 000000 | 0 | 35/500
001010 | 000000 | 0 | 31/500
001011 | 010000 | 1 | 28/500
001100 | 000001 | 1 | 39/500
001101 | 001001 | 1 | 33/500
001110 | 000001 | 1 | 36/500
001111 | 010001 | 1 | 27/500
010000 | 000100 | 1 | 25/500
010001 | 000000 | 0 | 32/500
010010 | 000000 | 0 | 23/500
010011 | 010000 | 1 | 34/500
010100 | 000000 | 0 | 32/500
010101 | 001000 | 1 | 33/500
010110 | 000000 | 0 | 30/500
010111 | 010000 | 1 | 39/500
011000 | 000100 | 1 | 32/500
011001 | 000000 | 0 | 36/500
011010 | 000000 | 0 | 46/500
011011 | 010000 | 1 | 35/500
011100 | 000000 | 0 | 20/500
011101 | 001000 | 1 | 33/500
011110 | 000000 | 0 | 37/500
011111 | 010000 | 1 | 35/500
100000 | 100000 | 1 | 20/500
100001 | 100000 | 1 | 30/500
100010 | 100000 | 1 | 31/500
100011 | 100000 | 1 | 23/500
100100 | 100000 | 1 | 30/500
100101 | 100000 | 1 | 30/500
100110 | 100000 | 1 | 26/500
100111 | 100000 | 1 | 29/500
101000 | 000000 | 0 | 24/500
101001 | 001000 | 1 | 33/500
101010 | 000010 | 1 | 33/500
101011 | 000010 | 1 | 30/500
101100 | 000010 | 1 | 31/500
101101 | 001010 | 1 | 28/500
101110 | 000000 | 0 | 43/500
101111 | 000000 | 0 | 32/500
110000 | 000000 | 0 | 30/500
110001 | 001000 | 1 | 28/500
110010 | 000010 | 1 | 35/500
110011 | 000010 | 1 | 36/500
110100 | 000010 | 1 | 31/500
110101 | 001010 | 1 | 33/500
110110 | 000000 | 0 | 28/500
110111 | 000000 | 0 | 28/500
111000 | 000000 | 0 | 35/500
111001 | 001000 | 1 | 30/500
111010 | 000000 | 0 | 32/500
111011 | 000000 | 0 | 30/500
111100 | 000000 | 0 | 31/500
111101 | 001000 | 1 | 35/500
111110 | 000000 | 0 | 26/500
111111 | 000000 | 0 | 41/500

#一括サンプリングしてピクセル点灯で再生
pos = (0, 1, 2, 3, 4, 5, 36)
box = np.zeros(64, 'uint8')
base = np.array([2**5, 2**4, 2**3, 2**2, 2**1, 2**0], int)

#一括サンプリング
qc.measure_all()
backend = AerSimulator(method='matrix_product_state')
result = backend.run(qc, shots=100, memory=True).result().get_memory()

for r in result:
    sel = ''.join(r[::-1][p] for p in pos)
    #出力ビット
    bit = np.array(list(sel[:6]), int)
    #判定ビット
    judge = bool(int(sel[-1]))
    #判定が1のときだけ描画
    if judge:
        #10進数ワンホットに戻す
        box = box // 1.0001 # 過去の色を減衰させる
        idx = np.sum(base * bit)
        box[idx] = 255
        
        #表示
        img = box.reshape(8, 8)
        plt.imshow(img, vmin=0, vmax=255)
        plt.title(bit)
        plt.show()
        plt.close()