結果
| 問題 | 
                            No.2523 Trick Flower
                             | 
                    
| コンテスト | |
| ユーザー | 
                             | 
                    
| 提出日時 | 2023-10-27 21:40:34 | 
| 言語 | PyPy3  (7.3.15)  | 
                    
| 結果 | 
                             
                                AC
                                 
                             
                            
                         | 
                    
| 実行時間 | 1,364 ms / 4,000 ms | 
| コード長 | 1,694 bytes | 
| コンパイル時間 | 250 ms | 
| コンパイル使用メモリ | 82,572 KB | 
| 実行使用メモリ | 261,348 KB | 
| 最終ジャッジ日時 | 2024-09-25 13:45:50 | 
| 合計ジャッジ時間 | 20,106 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge2 / judge1 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 3 | 
| other | AC * 31 | 
ソースコード
import sys
from itertools import permutations
import heapq
from collections import deque
import random
input = lambda :sys.stdin.readline().rstrip()
mi = lambda :map(int,input().split())
li = lambda :list(mi())
N = int(input())
A = li()
B = li()
C = li()
C = [c-1 for c in C]
cycle = [C[i] for i in range(N)]
for _ in range(30):
    cycle = [cycle[cycle[i]] for i in range(N)]
cycle = set(cycle)
comp_cycle = [i for i in range(N)]
for i in range(N):
    if i not in cycle or comp_cycle[i]!=i:
        continue
    t = len(comp_cycle)
    tmp = [i]
    while True:
        nxt = C[tmp[-1]]
        if nxt!=i:
            tmp.append(nxt)
        else:
            break
    for v in tmp:
        comp_cycle[v] = t
    comp_cycle.append(t)
    A.append(sum([A[i] for i in tmp]))
    B.append(sum([B[i] for i in tmp]))
n = len(comp_cycle)
edge = [[] for v in range(n)]
parent = [-1] * n
for i in range(N):
    if comp_cycle[i]!=i:
        continue
    parent[comp_cycle[i]] = comp_cycle[C[i]]
    edge[comp_cycle[C[i]]].append(comp_cycle[i])
topo = []
for root in range(N,n):
    deq = deque([root])
    while deq:
        v = deq.popleft()
        topo.append(v)
        for nv in edge[v]:
            deq.append(nv)
dp = [-1] * n
def cond(k):
    for v in topo[::-1]:
        dp[v] = k * B[v]
        for nv in edge[v]:
            dp[v] += dp[nv]
        if dp[v] > A[v]:
            dp[v] -= A[v]
        else:
            dp[v] = 0
    
    for root in range(N,n):
        if dp[root]:
            return False
    return True
ok = 0
ng = 10**15 + 1
while ng-ok>1:
    mid = (ok+ng)//2
    if cond(mid):
        ok = mid
    else:
        ng = mid
print(ok)