結果

問題 No.5020 Averaging
コンテスト
ユーザー Klavis
提出日時 2024-02-25 14:53:46
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 60 ms / 1,000 ms
コード長 1,665 bytes
コンパイル時間 171 ms
コンパイル使用メモリ 81,828 KB
実行使用メモリ 66,496 KB
スコア 14,429,657
最終ジャッジ日時 2024-02-25 14:53:53
合計ジャッジ時間 4,524 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge11
純コード判定しない問題か言語
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 50
権限があれば一括ダウンロードができます

ソースコード

diff # 

import random
import math

N = int(input())  # カードの枚数
A = [0] * N  # 各カードの表面の数値
B = [0] * N  # 各カードの裏面の数値

for i in range(N):
    A[i], B[i] = map(int, input().split())

# 目標値
target = 5 * 10**17

# 初期温度、終了温度、冷却率
T0 = 10000000
Tf = 1000
alpha = 0.99

def evaluate(A, B, target):
    """評価関数：カード1の表裏の数値と目標値との差の絶対値の和"""
    return abs(A[0] - target) + abs(B[0] - target)

def simulated_annealing(A, B, T0, Tf, alpha, max_operations=50):
    operations = []
    T = T0
    current_score = evaluate(A, B, target)
    while T > Tf and len(operations) < max_operations:
        # ランダムに2枚のカードを選択
        u, v = random.sample(range(N), 2)
        # 操作を試みる
        new_A_u, new_A_v = (A[u] + A[v]) // 2, (A[u] + A[v]) // 2
        new_B_u, new_B_v = (B[u] + B[v]) // 2, (B[u] + B[v]) // 2
        # 一時的に更新
        A[u], A[v], B[u], B[v] = new_A_u, new_A_v, new_B_u, new_B_v
        new_score = evaluate(A, B, target)
        delta = new_score - current_score
        if delta < 0 or random.random() < math.exp(-delta / T):
            # 更新を確定
            current_score = new_score
            operations.append((u+1, v+1))
        else:
            # 更新を取り消し
            A[u], A[v], B[u], B[v] = A[u], A[v], B[u], B[v]
        # 温度の更新
        T *= alpha
    return operations

operations = simulated_annealing(A, B, T0, Tf, alpha)

# 出力
print(len(operations))
for u, v in operations:
    print(u, v)


print(abs(A[0] - target) + abs(B[0] - target))
0