結果

問題 No.2523 Trick Flower
コンテスト
ユーザー MasKoaTS
提出日時 2023-10-27 15:25:15
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 950 ms / 4,000 ms
コード長 3,580 bytes
コンパイル時間 466 ms
コンパイル使用メモリ 82,240 KB
実行使用メモリ 268,484 KB
最終ジャッジ日時 2024-09-25 13:02:18
合計ジャッジ時間 15,679 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 31
権限があれば一括ダウンロードができます

ソースコード

diff # 

import itertools as iter
import collections as coll
import heapq as hq
import bisect as bis
from decimal import Decimal as dec
from functools import cmp_to_key
import math
import sys
#import pypyjit
#pypyjit.set_param('max_unroll_recursion=-1')
sys.setrecursionlimit(10 ** 6)
inp = sys.stdin.readline
input = lambda : inp()[:-1]
getN = lambda : int(inp())
getNs = lambda : map(int, inp().split())
getList = lambda :list(map(int, inp().split()))
getStrs = lambda n : [input() for _ in [0] * n]
getEdges = lambda n : [[x - 1 for x in getNs()] for _ in [0] * n]
def yexit(): print("Yes"); exit(0)
def nexit(): print("No"); exit(0)
pi = 3.141592653589793
mod = 1000000007
MOD = 998244353
INF = 4611686018427387903
dx = [1, 0, -1, 0];  dy = [0, 1, 0, -1]
#di = coll.defaultdict(int)

"""
強連結成分分解（SCC）
ref : https://manabitimes.jp/math/1250
"""

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]
        group = [[] 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)
            group[lb].append(v)
        return dag, group

 
"""
Main Code
"""

N = getN()
A, B, C = [getList() for _ in [0] * 3]

scc = SCC(N)
for i, x in enumerate(C):
    scc.add_edge(x - 1, 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):
        if(rtree[x] == -1 and dp[x] < 0):
            return False
        dp[rtree[x]] += min(dp[x], 0)
    return True

ok = 0
ng = 10 ** 14 + 1
while(abs(ok - ng) > 1):
    mid = (ok + ng) >> 1
    if check(mid):
        ok = mid
    else:
        ng = mid
print(ok)
0