結果

問題 No.2523 Trick Flower
ユーザー MasKoaTSMasKoaTS
提出日時 2023-10-26 17:20:45
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 599 ms / 4,000 ms
コード長 2,894 bytes
コンパイル時間 178 ms
コンパイル使用メモリ 82,152 KB
実行使用メモリ 260,332 KB
最終ジャッジ日時 2024-09-25 12:26:09
合計ジャッジ時間 10,231 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 40 ms
54,104 KB
testcase_01 AC 37 ms
53,696 KB
testcase_02 AC 36 ms
52,612 KB
testcase_03 AC 36 ms
53,216 KB
testcase_04 AC 38 ms
52,704 KB
testcase_05 AC 37 ms
53,400 KB
testcase_06 AC 36 ms
52,552 KB
testcase_07 AC 38 ms
53,492 KB
testcase_08 AC 66 ms
72,964 KB
testcase_09 AC 65 ms
72,340 KB
testcase_10 AC 50 ms
62,672 KB
testcase_11 AC 51 ms
63,564 KB
testcase_12 AC 50 ms
65,216 KB
testcase_13 AC 65 ms
71,840 KB
testcase_14 AC 65 ms
71,728 KB
testcase_15 AC 599 ms
260,332 KB
testcase_16 AC 562 ms
258,032 KB
testcase_17 AC 270 ms
132,276 KB
testcase_18 AC 231 ms
122,596 KB
testcase_19 AC 565 ms
220,788 KB
testcase_20 AC 565 ms
202,444 KB
testcase_21 AC 526 ms
206,424 KB
testcase_22 AC 539 ms
238,892 KB
testcase_23 AC 579 ms
188,496 KB
testcase_24 AC 322 ms
158,188 KB
testcase_25 AC 481 ms
170,576 KB
testcase_26 AC 457 ms
178,656 KB
testcase_27 AC 481 ms
195,512 KB
testcase_28 AC 423 ms
186,608 KB
testcase_29 AC 161 ms
88,704 KB
testcase_30 AC 386 ms
144,252 KB
testcase_31 AC 122 ms
79,232 KB
testcase_32 AC 256 ms
118,716 KB
testcase_33 AC 520 ms
205,340 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
input = sys.stdin.readline


class SCC:
    def __init__(self, n):
        self.n = n
        self.route = [[] for _ in [0] * n]
        self.rev_route = [[] for _ in [0] * n]
        self.label = [-1] * n
        self.count = 0
    
    def add_edge(self, v_first, v_second):
        self.route[v_first].append(v_second)
        self.rev_route[v_second].append(v_first)

    def _build(self):
        self.v_list = []
        self.visited = set([])
        for v in range(self.n):
            if(v in self.visited):
                continue
            self._dfs(v)
        for v in reversed(self.v_list):
            if(self.label[v] != -1):
                continue
            self._rev_dfs(v)
            self.count += 1

    def _dfs(self, v_start):
        stack = [v_start, 0]
        while(stack):
            v, id = stack[-2:]
            if(id == 0 and v in self.visited):
                stack.pop()
                stack.pop()
                continue
            self.visited.add(v)
            if(id < len(self.route[v])):
                stack[-1] += 1
                stack.extend([self.route[v][id], 0])
                continue
            stack.pop()
            self.v_list.append(stack.pop())

    def _rev_dfs(self, v_start):
        stack = [v_start]
        self.label[v_start] = self.count
        while(stack):
            v_now = stack.pop()
            for v_next in self.rev_route[v_now]:
                if(self.label[v_next] != -1):
                    continue
                self.label[v_next] = self.count
                stack.append(v_next)

    def construct(self):
        self._build()
        dag = [[] for _ in [0] * self.count]
        groups = [[] for _ in [0] * self.count]
        for v, lb in enumerate(self.label):
            for v_next in self.route[v]:
                lb_next = self.label[v_next]
                if(lb == lb_next):
                    continue
                dag[lb].append(lb_next)
            groups[lb].append(v)
        return dag, groups


"""
Main Code
"""

N = int(input())
A = list(map(int, input().split()))
B = list(map(int, input().split()))
C = [x - 1 for x in map(int, input().split())]

scc = SCC(N)
for i, x in enumerate(C):
    if(i == x):
        continue
    scc.add_edge(x, i)
dag, groups = scc.construct()

n = len(groups)
a = [sum(A[x] for x in group) for group in groups]
b = [sum(B[x] for x in group) for group in groups]
rtree = [-1] * n
for x, ys in enumerate(dag):
    for y in ys:
        rtree[y] = x

def check(mid):
    dp = [x - mid * y for x, y in zip(a, b)]
    for x in range(n - 1, -1, -1):
        y = rtree[x]
        if(y == -1 and dp[x] < 0):
            return False
        dp[y] += min(dp[x], 0)
    return True

ok = 0
ng = sum(a) // sum(b) + 1
while(abs(ok - ng) > 1):
    mid = (ok + ng) >> 1
    if check(mid):
        ok = mid
    else:
        ng = mid
print(ok)
0