結果

問題 No.1324 Approximate the Matrix
コンテスト
ユーザー roaris
提出日時 2020-12-23 22:56:34
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,076 ms / 2,000 ms
コード長 2,360 bytes
コンパイル時間 731 ms
コンパイル使用メモリ 82,176 KB
実行使用メモリ 88,960 KB
最終ジャッジ日時 2024-09-21 16:35:49
合計ジャッジ時間 18,101 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 42
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
input = sys.stdin.readline
from heapq import *

class Mincostflow:
    def __init__(self, N):
        self.N = N
        self.G = [[] for _ in range(N)]
        
    def add_edge(self, u, v, cap, cost):
        self.G[u].append([v, cap, cost, len(self.G[v])])
        self.G[v].append([u, 0, -cost, len(self.G[u])-1])
    
    def min_cost_flow(self, s, t, f):
        res = 0
        h = [0]*self.N
        
        while f>0:
            dist = [10**18]*self.N
            dist[s] = 0
            pq = [(0, s)]
            prev_v = [-1]*self.N
            prev_e = [-1]*self.N
            
            while pq:
                d, v = heappop(pq)
                
                if dist[v]<d:
                    continue
                
                for i in range(len(self.G[v])):
                    nv, cap, cost, _  = self.G[v][i]
                    
                    if cap>0 and dist[nv]>dist[v]+cost+h[v]-h[nv]:
                        dist[nv] = dist[v]+cost+h[v]-h[nv]
                        prev_v[nv] = v
                        prev_e[nv] = i
                        heappush(pq, (dist[nv], nv))
            
            if dist[t]==10**18:
                return -1
            
            for i in range(self.N):
                h[i] += dist[i]
                
            d = f
            v = t
            
            while v!=s:
                d = min(d, self.G[prev_v[v]][prev_e[v]][1])
                v = prev_v[v]
            
            f -= d
            res += d*h[t]
            v = t
            
            while v!=s:
                self.G[prev_v[v]][prev_e[v]][1] -= d
                rev = self.G[prev_v[v]][prev_e[v]][3]
                self.G[v][rev][1] += d
                v = prev_v[v]
        
        return res

N, K = map(int, input().split())
A = list(map(int, input().split()))
B = list(map(int, input().split()))
P = [list(map(int, input().split())) for _ in range(N)]
mcf = Mincostflow(2*N+2)

for i in range(1, N+1):
    mcf.add_edge(0, i, A[i-1], 0)

for i in range(N+1, 2*N+1):
    mcf.add_edge(i, 2*N+1, B[i-N-1], 0)

for i in range(N):
    for j in range(N):
        for k in range(A[i]):
            mcf.add_edge(i+1, N+1+j, 1, (k+1-P[i][j])**2-(k-P[i][j])**2)

ans = mcf.min_cost_flow(0, 2*N+1, K)

for i in range(N):
    for j in range(N):
        ans += P[i][j]**2

print(ans)
0