結果

問題 No.1324 Approximate the Matrix
コンテスト
ユーザー lam6er
提出日時 2025-03-31 17:37:10
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 2,036 bytes
コンパイル時間 197 ms
コンパイル使用メモリ 82,420 KB
実行使用メモリ 88,556 KB
最終ジャッジ日時 2025-03-31 17:38:22
合計ジャッジ時間 16,684 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 9 WA * 18 TLE * 8 -- * 7
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys

def main():
    input = sys.stdin.read().split()
    idx = 0
    N, K = int(input[idx]), int(input[idx+1])
    idx +=2
    A = list(map(int, input[idx:idx+N]))
    idx +=N
    B = list(map(int, input[idx:idx+N]))
    idx +=N
    P = []
    for i in range(N):
        row = list(map(int, input[idx:idx+N]))
        idx +=N
        P.append(row)
    
    u = [0.0]*N
    v = [0.0]*N
    
    EPS = 1e-8
    iterations = 100
    
    for _ in range(iterations):
        # Update u
        for i in range(N):
            target = 2.0 * A[i]
            D = [2.0 * P[i][j] + v[j] for j in range(N)]
            
            left = -1e20
            right = 1e20
            for __ in range(100):
                mid = (left + right)/2
                total = 0.0
                for dj in D:
                    val = mid + dj
                    if val > 0:
                        total += val
                if total < target - EPS:
                    left = mid
                else:
                    right = mid
            u[i] = (left + right)/2
        
        # Update v
        for j in range(N):
            target = 2.0 * B[j]
            C = [2.0 * P[i][j] + u[i] for i in range(N)]
            
            left = -1e20
            right = 1e20
            for __ in range(100):
                mid = (left + right)/2
                total = 0.0
                for ci in C:
                    val = mid + ci
                    if val >0:
                        total += val
                if total < target - EPS:
                    left = mid
                else:
                    right = mid
            v[j] = (left + right)/2
    
    Q = [[0]*N for _ in range(N)]
    for i in range(N):
        for j in range(N):
            Q[i][j] = max(0.0, (u[i] + v[j] + 2.0 * P[i][j]) / 2.0 )
    
    res = 0.0
    for i in range(N):
        for j in range(N):
            diff = Q[i][j] - P[i][j]
            res += diff * diff
    
    print(int(round(res)))
    
if __name__ == '__main__':
    main()
0