Qiskit Aqua量子演算法實作與應用

def create_oracle(oracle_method):
    oracle_text = {"log": "~A & ~B & C", "bit": "00001000"}
    global num_iterations
    print("輸入oracle輸入字串，例如：" + oracle_text[oracle_method] + "\n或輸入'def'使用預設字串。")
    oracle_input = input('\nOracle輸入：\n')
    if oracle_input == "def":
        oracle_type = oracle_text[oracle_method]
    else:
        oracle_type = oracle_input
    num_iterations = log_length(oracle_type, oracle_method)
    return oracle_type

def create_grover(oracle_type, oracle_method):
    if oracle_method == "log":
        algorithm = Grover(LogicalExpressionOracle(oracle_type), num_iterations=num_iterations)
        oracle_circuit = Grover(LogicalExpressionOracle(oracle_type)).construct_circuit()
    else:
        algorithm = Grover(TruthTableOracle(oracle_type), num_iterations=num_iterations)
        oracle_circuit = Grover(TruthTableOracle(oracle_type)).construct_circuit()
    display(oracle_circuit.draw(output="mpl"))
    display(algorithm)
    return algorithm

def run_grover(algorithm, oracle_type, oracle_method):
    backend = Aer.get_backend('qasm_simulator')
    result = algorithm.run(backend)
    print("Oracle方法：", oracle_method)
    print("Oracle輸入：", oracle_type)
    print("Aer模擬結果：", result['top_measurement'])
    display(plot_histogram(result['measurement']))
    
    if not IBMQ.active_account():
        IBMQ.load_account()
    provider = IBMQ.get_provider()
    filtered_backend = least_busy(provider.backends(n_qubits=5, operational=True, simulator=False))
    result = algorithm.run(filtered_backend)
    print("Oracle方法：", oracle_method)
    print("Oracle輸入：", oracle_type)
    print("IBMQ " + filtered_backend.name() + " 結果：", result['top_measurement'])
    display(plot_histogram(result['measurement']))
    print(result)

def main():
    oracle_method = "log"
    while oracle_method != 0:
        print("第10章：使用Aqua執行Grover搜尋")
        print("------------------------------")
        oracle_method = input("選擇Oracle方法（log或bit）：\n")
        type = create_oracle(oracle_method)
        algorithm = create_grover(type, oracle_method)
        run_grover(algorithm, type, oracle_method)

from qiskit import Aer, IBMQ
from qiskit.aqua.algorithms import Shor
import time

def display_shor(N):
    print("建立Shor電路...")
    # Shor演算法實作

shor_circuit = Shor(N=N).construct_circuit()
print(shor_circuit)
print("Circuit data\n\nDepth: ", shor_circuit.depth(), "\nWidth: ", shor_circuit.width(), "\nSize:", shor_circuit.size())

def run_shor(N):
    if N <= 64:
        backend = Aer.get_backend('qasm_simulator')
    else:
        if not IBMQ.active_account():
            IBMQ.load_account()
        provider = IBMQ.get_provider()
        backend = provider.get_backend('ibmq_qasm_simulator')
    
    results = Shor(N=N).run(backend)
    print("\nResults:")
    if results['factors'] == []:
        print("No prime factors: ", str(N), "=", str(N))
    elif isinstance(results['factors'][0], int):
        print("Prime factors: ", str(N), "=", results['factors'][0], "^2")
    else:
        print("Prime factors: ", str(N), "=", results['factors'][0][0], "*", results['factors'][0][1])

def main():
    while True:
        number = int(input("\nEnter an odd number N > 1 (0 to exit):\n"))
        if number > 1 and number % 2 > 0:
            type = input("Enter R to run the Shor algorithm, D to display the circuit.\n")
            start_time = time.time()
            if type.upper() == "D":
                display_shor(number)
            elif type.upper() == "R":
                run_shor(number)
            end_time = time.time()
            print("Elapsed time: ", "%.2f" % (end_time - start_time), "s")
        else:
            print("The number must be odd and larger than 1.")

def load_font(font_path):
    """
    載入符合OFL授權的字型檔案，並檢查是否包含必要的著作權宣告。
    
    Args:
        font_path (str): 字型檔案路徑
    
    Returns:
        bool: 是否符合OFL規範
    """
    try:
        with open(font_path, 'r', encoding='utf-8') as font_file:
            content = font_file.read()
            # 檢查是否包含著作權宣告
            if "Copyright" in content and "SIL Open Font License" in content:
                return True
            else:
                return False
    except Exception as e:
        print(f"載入字型失敗: {e}")
        return False

# 使用範例
font_path = "path/to/font.ttf"
if load_font(font_path):
    print("字型符合OFL規範")
else:
    print("字型不符合OFL規範，請檢查授權內容")