結果

問題 No.1914 Directed by Two Sequences
ユーザー chineristACchineristAC
提出日時 2021-10-31 19:02:09
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,788 ms / 3,000 ms
コード長 5,915 bytes
コンパイル時間 441 ms
コンパイル使用メモリ 87,252 KB
実行使用メモリ 193,108 KB
最終ジャッジ日時 2023-08-22 06:26:21
合計ジャッジ時間 68,119 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 93 ms
72,416 KB
testcase_01 AC 94 ms
72,292 KB
testcase_02 AC 670 ms
100,360 KB
testcase_03 AC 606 ms
101,416 KB
testcase_04 AC 608 ms
101,212 KB
testcase_05 AC 607 ms
101,740 KB
testcase_06 AC 595 ms
101,052 KB
testcase_07 AC 715 ms
99,812 KB
testcase_08 AC 701 ms
100,400 KB
testcase_09 AC 711 ms
101,176 KB
testcase_10 AC 1,578 ms
176,680 KB
testcase_11 AC 1,557 ms
176,504 KB
testcase_12 AC 1,536 ms
171,012 KB
testcase_13 AC 1,538 ms
164,564 KB
testcase_14 AC 1,621 ms
167,812 KB
testcase_15 AC 1,560 ms
167,156 KB
testcase_16 AC 1,630 ms
184,384 KB
testcase_17 AC 1,644 ms
182,912 KB
testcase_18 AC 1,592 ms
183,284 KB
testcase_19 AC 1,012 ms
137,460 KB
testcase_20 AC 1,024 ms
136,768 KB
testcase_21 AC 1,077 ms
143,608 KB
testcase_22 AC 1,669 ms
177,964 KB
testcase_23 AC 1,645 ms
179,100 KB
testcase_24 AC 1,577 ms
177,820 KB
testcase_25 AC 1,556 ms
170,916 KB
testcase_26 AC 1,467 ms
165,876 KB
testcase_27 AC 1,760 ms
170,752 KB
testcase_28 AC 1,735 ms
171,196 KB
testcase_29 AC 1,711 ms
171,868 KB
testcase_30 AC 1,629 ms
167,724 KB
testcase_31 AC 1,614 ms
179,292 KB
testcase_32 AC 1,774 ms
168,736 KB
testcase_33 AC 677 ms
116,972 KB
testcase_34 AC 1,788 ms
193,108 KB
testcase_35 AC 937 ms
147,680 KB
testcase_36 AC 1,348 ms
162,976 KB
testcase_37 AC 1,367 ms
169,724 KB
testcase_38 AC 1,508 ms
181,164 KB
testcase_39 AC 1,673 ms
188,552 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

class SegmentTree:
    def __init__(self, init_val, segfunc, ide_ele):
        n = len(init_val)
        self.segfunc = segfunc
        self.ide_ele = ide_ele
        self.num = 1 << (n - 1).bit_length()
        self.tree = [ide_ele] * 2 * self.num
        self.size = n
        for i in range(n):
            self.tree[self.num + i] = init_val[i]
        for i in range(self.num - 1, 0, -1):
            self.tree[i] = self.segfunc(self.tree[2 * i], self.tree[2 * i + 1])

    def update(self, k, x):
        k += self.num
        self.tree[k] = x
        while k > 1:
            k >>= 1
            self.tree[k] = self.segfunc(self.tree[2*k], self.tree[2*k+1])

    def query(self, l, r):
        if r==self.size:
            r = self.num

        res = self.ide_ele

        l += self.num
        r += self.num
        right = []
        while l < r:
            if l & 1:
                res = self.segfunc(res, self.tree[l])
                l += 1
            if r & 1:
                right.append(self.tree[r-1])
            l >>= 1
            r >>= 1

        for e in right[::-1]:
            res = self.segfunc(res,e)
        return res

def solve(N,A,B,M,E):
    from collections import deque
    A = [a-1 for a in A]
    B = [b-1 for b in B]
    res = set()

    edge = [[] for v in range(N)]
    for u,v in E:
        edge[u-1].append(v-1)
        edge[v-1].append(u-1)
    
    group = [-1 for v in range(N)]
    for v in range(N):
        mex = [False for i in range(len(edge[v])+1)]
        for nv in edge[v]:
            if nv < v and group[nv] < len(mex):
                mex[group[nv]] = True
        for i in range(len(mex)):
            if not mex[i]:
                group[v] = i
                break

    n = max(group) + 1
    clique = [[] for g in range(n)]
    for v in range(N):
        clique[group[v]].append(v)

    def direct(i,j):
        if i < j:
            return A[i] < B[j]
        else:
            return B[i] < A[j]

    def merge(a,b):
        res = []
        bi = 0
        for ai in range(len(a)):
            while bi<len(b) and direct(b[bi],a[ai]):
                res.append(b[bi])
                bi += 1
            res.append(a[ai])
        res += b[bi:]
        return res

    idx = [-1 for v in range(N)]
    for g in range(n):
        deq = deque([[v] for v in clique[g]])
        while 1 < len(deq):
            a = deq.popleft()
            b = deq.popleft()
            deq.append(merge(a,b))
        clique[g] = deq[0]

        for i in range(len(clique[g])):
            idx[clique[g][i]] = i
        for u,v in zip(clique[g],clique[g][1:]):
            res.add((u,v))
    
    memo_to = [N for v in range(N)]
    memo_from = [-1 for v in range(N)]
    ci = [i for i in range(n)]
    ci.sort(lambda g:len(clique[g]))
    idx_on_ci = [-1 for i in range(n)]
    for i in range(n):
        idx_on_ci[ci[i]] = i
    ban = [[] for v in range(N)]
    Vs = []

    

    for i in range(n):
        target = ci[i]
        Vs += clique[target]
        Vs.sort()
        val = sorted([A[v] for v in Vs]+[B[v] for v in Vs])
        comp = {e:i for i,e in enumerate(val)}
        for nv in clique[target]:
            for v in edge[nv]:
                if idx_on_ci[group[v]] <= i:
                    ban[v].append(nv)
        
        seg_to_large = SegmentTree([N]*len(comp),min,N)
        seg_to_small = SegmentTree([N]*len(comp),min,N)
        seg_from_large = SegmentTree([-1]*len(comp),max,-1)
        seg_from_small = SegmentTree([-1]*len(comp),max,-1)
        k = len(comp)

        for v in Vs[::-1]:
            for nv in ban[v]:
                if v < nv:
                    seg_to_large.update(comp[B[nv]],N)
            memo_to[v] = min(memo_to[v],seg_to_large.query(comp[A[v]],k))
            for nv in ban[v]:
                if v < nv:
                    seg_to_large.update(comp[B[nv]],idx[nv])
            if group[v]==target:
                seg_to_large.update(comp[B[v]],idx[v])
        
        for v in Vs:
            for nv in ban[v]:
                if nv < v:
                    seg_to_small.update(comp[A[nv]],N)
            memo_to[v] = min(memo_to[v],seg_to_small.query(comp[B[v]],k))
            for nv in ban[v]:
                if nv < v:
                    seg_to_small.update(comp[A[nv]],idx[nv])
            if group[v]==target:
                seg_to_small.update(comp[A[v]],idx[v])
        
        for v in Vs[::-1]:
            for nv in ban[v]:
                if v < nv:
                    seg_from_large.update(comp[B[nv]],-1)
            memo_from[v] = max(memo_from[v],seg_from_large.query(0,comp[A[v]]))
            for nv in ban[v]:
                if v < nv:
                    seg_from_large.update(comp[B[nv]],idx[nv])
            if group[v]==target:
                seg_from_large.update(comp[B[v]],idx[v])
        
        for v in Vs:
            for nv in ban[v]:
                if nv < v:
                    seg_from_small.update(comp[A[nv]],-1)
            memo_from[v] = max(memo_from[v],seg_from_small.query(0,comp[B[v]]))
            for nv in ban[v]:
                if nv < v:
                    seg_from_small.update(comp[A[nv]],idx[nv])
            if group[v]==target:
                seg_from_small.update(comp[A[v]],idx[v])
        
        for v in Vs:
            if memo_to[v]!=N:
                nv = clique[target][memo_to[v]]     
                res.add((v,nv))      
                memo_to[v] = N
            
            if memo_from[v]!=-1:
                pv = clique[target][memo_from[v]]
                res.add((pv,v))
                memo_from[v] = -1
            
            ban[v] = []
            
    
    print(len(res))
    for u,v in res:
        print(u+1,v+1)

import sys

input = lambda :sys.stdin.buffer.readline()
mi = lambda :map(int,input().split())
li = lambda :list(mi())

N,M = mi()
A = li()
B = li()
E = [tuple(mi()) for i in range(M)]
solve(N,A,B,M,E)
0