結果
問題 | No.1804 Intersection of LIS |
ユーザー |
|
提出日時 | 2024-06-26 19:54:03 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 908 ms / 2,000 ms |
コード長 | 11,343 bytes |
コンパイル時間 | 334 ms |
コンパイル使用メモリ | 82,484 KB |
実行使用メモリ | 239,396 KB |
最終ジャッジ日時 | 2024-06-26 19:54:27 |
合計ジャッジ時間 | 21,798 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 37 |
ソースコード
import os,sys,random,threading#sys.exit() 退出程序from random import randint,choice,shuffle#randint(a,b)从[a,b]范围随机选择一个数#choice(seq)seq可以是一个列表,元组或字符串,从seq中随机选取一个元素#shuffle(x)将一个可变的序列x中的元素打乱from copy import deepcopyfrom io import BytesIO,IOBasefrom types import GeneratorTypefrom functools import lru_cache,reduce#reduce(op,迭代对象)from bisect import bisect_left,bisect_right#bisect_left(x) 大于等于x的第一个下标#bisect_right(x) 大于x的第一个下标from collections import Counter,defaultdict,dequefrom itertools import accumulate,combinations,permutations#accumulate(a)用a序列生成一个累积迭代器,一般list化前面放个[0]做前缀和用#combinations(a,k)a序列选k个 组合迭代器#permutations(a,k)a序列选k个 排列迭代器from heapq import heapify,heappop,heappush#heapify将列表转为堆from typing import Generic,Iterable,Iterator,TypeVar,Union,Listfrom string import ascii_lowercase,ascii_uppercase,digits#小写字母,大写字母,十进制数字from math import ceil,floor,sqrt,pi,factorial,gcd,log,log10,log2,inf#ceil向上取整,floor向下取整 ,sqrt开方 ,factorial阶乘from decimal import Decimal,getcontext#Decimal(s) 实例化Decimal对象,一般使用字符串#getcontext().prec=100 修改精度from sys import stdin, stdout, setrecursionlimitinput = lambda: sys.stdin.readline().rstrip("\r\n")MI = lambda :map(int,input().split())li = lambda :list(MI())ii = lambda :int(input())mod = int(1e9 + 7) #998244353inf = int(1e20)py = lambda :print("YES")pn = lambda :print("NO")DIRS = [(0, 1), (1, 0), (0, -1), (-1, 0)] # 右下左上DIRS8 = [(0, 1), (1, 1), (1, 0), (1, -1), (0, -1), (-1, -1), (-1, 0),(-1, 1)] # →↘↓↙←↖↑↗class BIT2:"""区间更新,区间和查询 """__slots__ = "size", "_tree1", "_tree2"def __init__(self, n: int):self.size = nself._tree1 = [0]*(n+1)self._tree2 = [0]*(n+1)def add(self, left: int, right: int, delta: int) -> None:"""闭区间[left, right]加delta"""self._add(left, delta)self._add(right + 1, -delta)def query(self, left: int, right: int) -> int:"""闭区间[left, right]的和"""return self._query(right) - self._query(left - 1)def _add(self, index: int, delta: int) -> None:rawIndex = indexwhile index <= self.size:self._tree1[index] += deltaself._tree2[index] += (rawIndex - 1) * deltaindex += index & -indexdef _query(self, index: int) -> int:if index > self.size:index = self.sizerawIndex = indexres = 0while index > 0:res += rawIndex * self._tree1[index] - self._tree2[index]index -= index & -indexreturn resdef e(a):if a == min:return 10**15if a == max:return -10**15class LcaDoubling:def __init__(self, n, graph, root, lcaf):""""頂点は1-indexgraph[v] = [[u,w],[u,w],...]"""self.lcaf = lcafself.e = e(lcaf)# 根からの深さ、距離self.depths = [-1] * (n+1)self.distances = [-1] * (n+1)# self.ancestors[k][v]:頂点vから2^k先の親,# self.funcs[k][v]:頂点vから2^k先の親までのパス上の1辺の最大(最小)コストprev_ancestors, prev_funcs = self._init_dfs(n, graph, root)self.ancestors = [prev_ancestors]self.funcs = [prev_funcs]max_depth = max(self.depths)d = 1while d < max_depth:next_ancestors = []next_funcs = []for i in range(len(prev_ancestors)):p = prev_ancestors[i]dist = prev_funcs[i]next_ancestors.append(prev_ancestors[p])next_funcs.append(self.lcaf(dist, prev_funcs[p]))self.ancestors.append(next_ancestors)self.funcs.append(next_funcs)d <<= 1prev_ancestors = next_ancestorsprev_funcs = next_funcsdef _init_dfs(self, n, graph, root):que = [(root, -1, 0, 0)]direct_ancestors = [-1 for i in range(n+2)]direct_funcs = [self.e for i in range(n+2)]self.depths[root] = 0self.distances[root] = 0while que:crr, pre, dep, dist = que.pop()for nxt, w in graph[crr]:if nxt == pre:continuedirect_ancestors[nxt] = crrdirect_funcs[nxt] = wself.depths[nxt] = dep + 1self.distances[nxt] = dist + wque.append((nxt, crr, dep + 1, dist+w))return direct_ancestors, direct_funcsdef upstream(self, v, k):fans = self.eb = 0while k:if k & 1:fans = self.lcaf(fans, self.funcs[b][v])v = self.ancestors[b][v]k >>= 1b += 1return v, fansdef get_lca(self, u, v):du, dv = self.depths[u], self.depths[v]if du > dv:u, v = v, udu, dv = dv, dutu = utv, _ = self.upstream(v, dv-du)if tu == tv:return tufor k in range(du.bit_length())[::-1]:mu = self.ancestors[k][tu]mv = self.ancestors[k][tv]if mu != mv:tu = mutv = mvlca = self.ancestors[0][tu]return lcadef get_distance(self, u, v):lca = self.get_lca(u, v)return self.distances[u] + self.distances[v] - 2 * self.distances[lca]def get_lcaf(self, u, v):lca = self.get_lca(u, v)du = self.depths[u] - self.depths[lca]dv = self.depths[v] - self.depths[lca]_, uans = self.upstream(u, du)_, vans = self.upstream(v, dv)return self.lcaf(vans, uans)class HLD:def __init__(self, g, root):#无论是点还是dfn还是dep,都从1开始,默认0是无n=len(g)self.g=gself.fa=[0]*(n+5) #父节点,0表示无父节点self.size=[1]*(n+5) #子树大小self.dep=[0]*(n+5) #深度,根深度为1self.son=[0]*(n+5) #重儿子,0表示无儿子self.dfn=[0]*(n+5) #dfs序,子树终点的dfs序是dfn[i]+size[i]-1self.top=[0]*(n+5) #所在重链起点,起点就是自己self.rank=[0]*(n+5) #dfs序为i的节点编号fa=self.fa;size=self.size;dep=self.dep;son=self.sondfn=self.dfn;top=self.top;rank=self.rankstk=[[root,0,0]] #node,flag,fadep[root]=1while stk:u,flag,father=stk.pop()if flag:for v in g[u]:if v!=father:size[u]+=size[v]if son[u]==0 or size[v]>size[son[u]]:son[u]=velse:stk.append([u,1,father])for v in g[u]:if v!=father:stk.append([v,0,u])fa[v]=udep[v]=dep[u]+1stk=[[root,root]]tot=1while stk:u,tops=stk.pop()dfn[u]=totrank[tot]=utot+=1top[u]=topsif son[u]==0:continuefor v in g[u]:if v!=fa[u] and v!=son[u]:stk.append([v,v])stk.append([son[u],tops])def lca(self,u,v): #求u和v的最近公共祖先节点fa=self.fa;size=self.size;dep=self.dep;son=self.sondfn=self.dfn;top=self.top;rank=self.rankwhile top[u]!=top[v]:if dep[top[u]]>dep[top[v]]:u=fa[top[u]]else:v=fa[top[v]]return v if dep[u]>dep[v] else udef dis(self,u,v):dep=self.depreturn dep[u]+dep[v]-2*dep[self.lca(u,v)]def kth_fa(self,root,k): #求root节点的第k个祖先fa=self.fa;size=self.size;dep=self.dep;son=self.sondfn=self.dfn;top=self.top;rank=self.rankif k>=dep[root]: #无第k个祖先返回-1return -1while True:u=top[root]if dfn[root]-k>=dfn[u]:return rank[dfn[root]-k]k-=dfn[root]-dfn[u]+1root=fa[u]def route_query(self,u,v): #查询u到v简单路径fa=self.fa;size=self.size;dep=self.dep;son=self.sondfn=self.dfn;top=self.top;rank=self.rankroute=[]while top[u]!=top[v]:if dep[top[u]]<dep[top[v]]:u,v=v,uroute.append((dfn[top[u]],dfn[u]))u=fa[top[u]]if dep[u]>dep[v]:u,v=v,uroute.append((dfn[u],dfn[v]))return routeclass SparseTable:def __init__(self, data: list, func=gcd):# 稀疏表,O(nlogn)预处理,O(1)查询区间最值/按位与/按位或/gcd/lcm# 可以以 O(nlognlong)预处理,O(logn)查询区间lca# 下标从0开始self.func = funcself.st = st = [list(data)]i, N = 1, len(st[0])while 2 * i <= N+1:qz = st[-1]st.append([func(qz[j], qz[j + i]) for j in range(N - 2 * i + 1)])i <<= 1def query(self, begin: int, end: int): # 查询闭区间[begin, end]的最大值lg = (end - begin+1).bit_length() - 1return self.func(self.st[lg][begin], self.st[lg][end - (1 << lg) + 1])def bootstrap(f, stack=[]): #yielddef wrappedfunc(*args, **kwargs):if stack:return f(*args, **kwargs)else:to = f(*args, **kwargs)while True:if type(to) is GeneratorType:stack.append(to)to = next(to)else:stack.pop()if not stack:breakto = stack[-1].send(to)return toreturn wrappedfuncclass BIT:def __init__(self, n):self.n = nself.max_tree = [0] * (n + 1)def update(self, idx, val):while idx <= self.n:self.max_tree[idx] = max(self.max_tree[idx], val)idx += idx & (-idx)def query(self, idx):res = 0while idx > 0:res = max(res, self.max_tree[idx])idx -= idx & (-idx)return resn=ii()arr=li()dx=sorted(set(arr))d={x:i+1 for i,x in enumerate(dx)}bit=BIT(n+5)left=[0]*nright=[0]*nmx=0for i in range(n):a=arr[i]q=bit.query(d[a]-1)bit.update(d[a],q+1)mx=max(mx,q+1)left[i]=q+1cnt=[0]*(n+1)dx=dx[::-1]d={x:i+1 for i,x in enumerate(dx)}bit=BIT(n+5)for i in range(n-1,-1,-1):a=arr[i]q=bit.query(d[a]-1)bit.update(d[a],q+1)right[i]=q+1if left[i]+right[i]-1==mx:cnt[left[i]]+=1res=[]for i in range(n):if left[i]+right[i]-1==mx and cnt[left[i]]==1:res.append(arr[i])print(len(res))print(*res)