結果

問題 No.2523 Trick Flower
ユーザー MasKoaTSMasKoaTS
提出日時 2023-10-26 17:20:45
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 701 ms / 4,000 ms
コード長 2,894 bytes
コンパイル時間 313 ms
コンパイル使用メモリ 81,824 KB
実行使用メモリ 260,232 KB
最終ジャッジ日時 2023-10-26 17:20:57
合計ジャッジ時間 11,989 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 37 ms
53,588 KB
testcase_01 AC 36 ms
53,588 KB
testcase_02 AC 36 ms
53,588 KB
testcase_03 AC 36 ms
53,588 KB
testcase_04 AC 37 ms
53,588 KB
testcase_05 AC 36 ms
53,588 KB
testcase_06 AC 37 ms
53,588 KB
testcase_07 AC 39 ms
53,588 KB
testcase_08 AC 69 ms
72,748 KB
testcase_09 AC 68 ms
72,676 KB
testcase_10 AC 49 ms
63,928 KB
testcase_11 AC 50 ms
63,928 KB
testcase_12 AC 50 ms
63,940 KB
testcase_13 AC 64 ms
72,592 KB
testcase_14 AC 63 ms
70,536 KB
testcase_15 AC 701 ms
260,232 KB
testcase_16 AC 626 ms
257,852 KB
testcase_17 AC 273 ms
132,180 KB
testcase_18 AC 237 ms
122,524 KB
testcase_19 AC 608 ms
226,688 KB
testcase_20 AC 628 ms
202,404 KB
testcase_21 AC 568 ms
206,376 KB
testcase_22 AC 588 ms
233,552 KB
testcase_23 AC 577 ms
188,408 KB
testcase_24 AC 339 ms
154,716 KB
testcase_25 AC 566 ms
170,508 KB
testcase_26 AC 500 ms
178,596 KB
testcase_27 AC 534 ms
195,392 KB
testcase_28 AC 468 ms
186,256 KB
testcase_29 AC 170 ms
88,344 KB
testcase_30 AC 425 ms
144,324 KB
testcase_31 AC 121 ms
78,816 KB
testcase_32 AC 257 ms
118,572 KB
testcase_33 AC 565 ms
204,988 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