n=int(input())
s='0'+input()
m=int(input())
mp=[[] for i in range(n+1)]
for i in range(m):
a,b=map(int,input().split())
a+=1
b+=1
mp[a].append(b)
import sys
sys.setrecursionlimit(100000000)
dfn=[0 for i in range(n+1)]
low=[0 for i in range(n+1)]
cnt=0
stk=[]
instk=[0 for i in range(n+1)]
p=0
scc=[0 for i in range(n+1)]
def tarjan(x):
global cnt,p
cnt+=1
dfn[x]=low[x]=cnt
stk.append(x)
instk[x]=1
for i in mp[x]:
if dfn[i]==0:
tarjan(i)
low[x]=min(low[x],low[i])
elif instk[i]:
low[x]=min(low[x],dfn[i])
if dfn[x]==low[x]:
p+=1
y=stk.pop()
instk[y]=0
scc[y]=p
while x!=y:
y = stk.pop()
instk[y] = 0
scc[y] = p
for i in range(1,n+1):
if dfn[i]==0:tarjan(i)
del instk,low,dfn
ru=[0 for i in range(p+1)]
h=[0 for i in range(p+1)]
l=[0 for i in range(p+1)]
nmp=[[] for i in range(p+1)]
for i in range(1,n+1):
if s[i]=='0': h[scc[i]]+=1
else : l[scc[i]]+=1
for j in mp[i]:
if scc[i]!=scc[j]:
ru[scc[j]]+=1
nmp[scc[i]].append(scc[j])
del mp
def dfs(x,a,b):
ans=0
vis=0
for y in nmp[x]:
ans=max(dfs(y,a+h[y],b+l[y]),ans)
vis=1
if vis==0:return min(a,b)
else:return ans
ans=0
for i in range(1,p+1):
if ru[i]==0:
ans=max(ans,dfs(i,h[i],l[i]))
print(ans)
宝石收集,tarjan
WA==Wonderful Answer2024-05-27 15:06
相关推荐
纽扣66731 分钟前
【算法进阶之路】链表进阶:删除、合并、回文与排序全解析消失的旧时光-19431 小时前
统一并发模型:线程、Reactor、协程本质是一件事(从线程到协程 · 第6篇·终章)智者知已应修善业1 小时前
【51单片机不用数组动态数码管显示字符和LED流水灯】2023-10-3AI进化营-智能译站2 小时前
ROS2 C++开发系列16-智能指针管理传感器句柄|告别ROS2节点内存泄漏与野指针CS创新实验室2 小时前
从盘边到芯端——硬盘接口七十年变迁史xvhao20133 小时前
单源、多源最短路MATLAB代码顾问3 小时前
多种群协同进化算法(MPCE)求解大规模作业车间调度问题——附MATLAB代码FQNmxDG4S3 小时前
JVM内存模型详解:堆、栈、方法区与垃圾回收We་ct3 小时前
LeetCode 72. 编辑距离:动态规划经典题解AI科技星3 小时前
精细结构常数α作为SI 7大基本量纲统一耦合常数的量子几何涌现理论