結果

問題 No.5020 Averaging
ユーザー judgelawjudgelaw
提出日時 2024-02-25 16:52:09
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 944 ms / 1,000 ms
コード長 3,973 bytes
コンパイル時間 299 ms
コンパイル使用メモリ 81,700 KB
実行使用メモリ 83,312 KB
スコア 31,890,942
最終ジャッジ日時 2024-02-25 16:53:45
合計ジャッジ時間 50,205 ms
ジャッジサーバーID
（参考情報） 		judge10 / judge11
純コード判定しない問題か言語
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 940 ms
82,360 KB
testcase_01 AC 939 ms
83,072 KB
testcase_02 AC 939 ms
82,476 KB
testcase_03 AC 939 ms
82,724 KB
testcase_04 AC 940 ms
82,100 KB
testcase_05 AC 940 ms
82,336 KB
testcase_06 AC 939 ms
81,604 KB
testcase_07 AC 939 ms
82,448 KB
testcase_08 AC 939 ms
82,232 KB
testcase_09 AC 939 ms
82,160 KB
testcase_10 AC 940 ms
82,372 KB
testcase_11 AC 939 ms
82,848 KB
testcase_12 AC 939 ms
82,248 KB
testcase_13 AC 940 ms
82,632 KB
testcase_14 AC 940 ms
83,108 KB
testcase_15 AC 939 ms
81,568 KB
testcase_16 AC 943 ms
82,236 KB
testcase_17 AC 939 ms
82,696 KB
testcase_18 AC 939 ms
81,688 KB
testcase_19 AC 939 ms
81,844 KB
testcase_20 AC 939 ms
83,264 KB
testcase_21 AC 938 ms
82,356 KB
testcase_22 AC 939 ms
82,204 KB
testcase_23 AC 940 ms
82,092 KB
testcase_24 AC 939 ms
81,824 KB
testcase_25 AC 940 ms
82,372 KB
testcase_26 AC 944 ms
82,100 KB
testcase_27 AC 939 ms
82,508 KB
testcase_28 AC 939 ms
82,344 KB
testcase_29 AC 939 ms
82,044 KB
testcase_30 AC 940 ms
81,476 KB
testcase_31 AC 939 ms
82,264 KB
testcase_32 AC 943 ms
82,356 KB
testcase_33 AC 939 ms
82,848 KB
testcase_34 AC 938 ms
81,964 KB
testcase_35 AC 939 ms
82,612 KB
testcase_36 AC 939 ms
82,808 KB
testcase_37 AC 940 ms
81,648 KB
testcase_38 AC 940 ms
81,712 KB
testcase_39 AC 940 ms
82,032 KB
testcase_40 AC 939 ms
81,944 KB
testcase_41 AC 938 ms
82,144 KB
testcase_42 AC 940 ms
82,848 KB
testcase_43 AC 939 ms
81,832 KB
testcase_44 AC 939 ms
83,312 KB
testcase_45 AC 940 ms
82,340 KB
testcase_46 AC 939 ms
81,480 KB
testcase_47 AC 940 ms
82,780 KB
testcase_48 AC 939 ms
83,128 KB
testcase_49 AC 939 ms
82,084 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

from time import time
st=time()

import sys
input = lambda: sys.stdin.readline().rstrip()
from math import log10,exp
from copy import copy,deepcopy
from random import randint,random,choice



MAX_OP=50
FIRST_PHASE=0.15
TIME_LIMIT=0.9
TARGET_NUM=5*10**17
INF=10**20

# カード1のスコアだけ見ればいいのか
# 奇数回目の操作でカード1以外を操作する
# 偶数回目でカード1を操作する

# カード1のスコアを計算

def calc_front():return abs(tempAB[0][0]-TARGET_NUM)

def calc_back():return abs(tempAB[0][1]-TARGET_NUM)

def calc_Score(cnt):
    V1=calc_front()
    V2=calc_back()
    
    temp=max(V1,V2)

    if temp==0:return 2000050-cnt
    return int(2000000-100000*log10(temp+1))

# 二枚のカードの表裏の平均値を返す
def calc_ave(u,v):

    retf=(tempAB[u][0]+tempAB[v][0])//2
    retb=(tempAB[u][1]+tempAB[v][1])//2

    return retf,retb

def rewrite(u,v):
    tempf,tempb=calc_ave(u,v)
    tempAB[u][0]=tempf
    tempAB[v][0]=tempf

    tempAB[u][1]=tempb
    tempAB[v][1]=tempb


def Judge(N,AB,X,UV):
    def WrongAnswer(i,s):
        print(str(i) + "th Wrong Answer (" + s + ")")
        exit(0)

    # Step 1. テストケースの入力を受け取る

    A = [0] * (N + 1)
    B = [0] * (N + 1)
    for i in range(1, N+1):
        A[i], B[i] = AB[i-1]

    # Step 2. あなたの出力を受け取る
    if X < 0 or X > 50:
        WrongAnswer("X is out of range")

    U = [0] * (X + 1)
    V = [0] * (X + 1)
    for i in range(1, X+1):
        U[i], V[i] = UV[i-1]
        if U[i] == V[i]:
            WrongAnswer(i,"U[i] and V[i] are same")
        elif U[i] < 1 or U[i] > N:
            WrongAnswer(i,"U[i] is out of range")
        elif V[i] < 1 or V[i] > N:
            WrongAnswer(i,"V[i] is out of range")

    # Step 3. シミュレーションを行う
    for i in range(1, X+1):
        avgA = (A[U[i]] + A[V[i]]) // 2
        avgB = (B[U[i]] + B[V[i]]) // 2
        A[U[i]] = avgA
        A[V[i]] = avgA
        B[U[i]] = avgB
        B[V[i]] = avgB

    # Step 4. 出力
    ErrorA = abs(500000000000000000 - A[1])
    ErrorB = abs(500000000000000000 - B[1])
    Score = (int)(2000000.0 - 100000.0 * log10(1.0 * max(ErrorA, ErrorB) + 1.0))
    if ErrorA == 0 and ErrorB == 0:
        Score = 2000050 - X
    print("Accepted!")
    print("Error of A[1] = " + str(ErrorA))
    print("Error of B[1] = " + str(ErrorB))
    print("Score = " + str(Score))

N=int(input())
AB=[list(map(int, input().split()))  for _ in range(N)]


NowScore=0
UV=[]
ScoreMove=[]


while time()-st<=TIME_LIMIT:
    tempAB=deepcopy(AB)    
    tempUV=[]
    preop=-1

    opcnt=MAX_OP

    # 奇数回目は1にマッチさせるものを選ぶ
    # 偶数回目で1とマッチさせる
    for cnt in range(1,opcnt+1):

        if cnt&1:
            tarf=2*TARGET_NUM-tempAB[0][0]
            tarb=2*TARGET_NUM-tempAB[0][1]

            pretemp=10**20

            u,v=-1,-1
            # 適当にu,vの平均を計算してaとマッチさせられるか考える
            for _ in range(50):
                tu=randint(1,N-1)
                tv=randint(1,N-1)
                while tu==tv:tv=randint(1,N-1)

                tf,tb=calc_ave(tu,tv)
                temp=abs(tf-tarf)+abs(tb-tarb)
                if temp<pretemp:
                    pretemp=temp
                    u=tu
                    v=tv


            rewrite(u,v)
            tempUV.append((u+1,v+1))
            preop=u

        else:
            # 1とさっき計算したものをマッチさせる
            rewrite(0,preop)
            tempUV.append((1,preop+1))


    tempScore=calc_Score(opcnt)

    if tempScore>NowScore:
        NowScore=tempScore
        UV=deepcopy(tempUV)

    # else:
    #     T=(time()-st)/TIME_LIMIT
    #     if exp((tempScore-NowScore)/T)>random():
    #         NowScore=tempScore
    #         UV=deepcopy(tempUV)


    ScoreMove.append(NowScore)


L=len(UV)

print(L)
for u,v in UV:print(u,v)
0