結果

問題 No.2210 equence Squence Seuence
ユーザー タコイモタコイモ
提出日時 2023-02-11 11:50:31
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,564 ms / 2,000 ms
コード長 3,343 bytes
コンパイル時間 627 ms
コンパイル使用メモリ 82,392 KB
実行使用メモリ 126,172 KB
最終ジャッジ日時 2024-07-08 02:04:44
合計ジャッジ時間 16,304 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 40 ms
56,240 KB
testcase_01 AC 40 ms
56,620 KB
testcase_02 AC 39 ms
56,328 KB
testcase_03 AC 403 ms
88,868 KB
testcase_04 AC 405 ms
90,484 KB
testcase_05 AC 694 ms
94,032 KB
testcase_06 AC 570 ms
93,636 KB
testcase_07 AC 1,236 ms
116,808 KB
testcase_08 AC 40 ms
57,532 KB
testcase_09 AC 48 ms
66,772 KB
testcase_10 AC 47 ms
65,136 KB
testcase_11 AC 47 ms
65,160 KB
testcase_12 AC 48 ms
66,256 KB
testcase_13 AC 1,527 ms
125,884 KB
testcase_14 AC 1,540 ms
126,172 KB
testcase_15 AC 1,559 ms
125,908 KB
testcase_16 AC 1,564 ms
124,624 KB
testcase_17 AC 1,537 ms
124,748 KB
testcase_18 AC 39 ms
56,412 KB
testcase_19 AC 39 ms
56,000 KB
testcase_20 AC 39 ms
56,808 KB
testcase_21 AC 40 ms
56,744 KB
testcase_22 AC 39 ms
57,232 KB
testcase_23 AC 496 ms
114,772 KB
testcase_24 AC 378 ms
103,736 KB
testcase_25 AC 442 ms
109,812 KB
testcase_26 AC 578 ms
120,340 KB
testcase_27 AC 205 ms
90,408 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

base = 37 #こっちはsr未満ならなんでもよい、modはいじれない
mod = 2305843009213693951 #2**61 - 1 制限時間に余裕があり衝突したらこっちを使う
ss = 2**31
sr = 2**30
def abm(a,b): #a*b mod m
  au = a // ss
  ad = a % ss
  bu = b // ss
  bd = b % ss
  mid = au*bd + ad*bu
  midu = mid // sr
  midd = mid % sr
  return (au*bu*2 + midu + midd*ss + ad*bd)%mod

  
class RollingHash():
    def __init__(self, s, base, mod):
        self.mod = mod
        l = len(s)
        self.pw = pw = [1]*(l+1)
        self.h = h = [0]*(l+1)
        v = 0
        vv = 1
        for i in range(l):
            ad = v % ss
            au = v // ss
            mid = au*base
            midu = mid//sr
            midd = mid%sr
            #小文字のみ
            #h[i+1] = v = (midu+midd*ss+ad*base + ord(s[i])-96) % mod
            #大文字のみ or 大文字と小文字のみ
            #h[i+1] = v = (midu+midd*ss+ad*base + ord(s[i])-64) % mod
            #通常
            #h[i+1] = v = (midu+midd*ss+ad*base + ord(s[i])+1) % mod
            #自然数
            h[i+1] = v = (midu+midd*ss+ad*base + s[i]) % mod
            ad = vv % ss
            au = vv // ss
            mid = au*base
            midu = mid//sr
            midd = mid%sr
            pw[i+1] = vv = (midu+midd*ss+ad*base)% mod
    #[l,r)
    def get(self, l, r):
        return (self.h[r] - abm(self.h[l],self.pw[r-l])) % self.mod

  
import sys
#sys.setrecursionlimit(500000)
def I(): return int(sys.stdin.readline().rstrip())
def MI(): return map(int,sys.stdin.readline().rstrip().split())
def TI(): return tuple(map(int,sys.stdin.readline().rstrip().split()))
def LI(): return list(map(int,sys.stdin.readline().rstrip().split()))
def S(): return sys.stdin.readline().rstrip()
def LS(): return list(sys.stdin.readline().rstrip())
#for i, pi in enumerate(p):
from collections import defaultdict,deque
import bisect
import itertools
dic = defaultdict(int)
d = deque()

N,K = MI()
A = LI()
S = RollingHash(A,base,mod)
#https://aotamasaki.hatenablog.com/entry/meguru_bisect
#arg桁
def is_ok(arg,x,y):
  if arg == 0:
    return True
  if arg == N:
    return False
  if arg <= x and arg <= y:
    return True
  
  if arg <= x:
    X = S.get(0,arg)
  else:
    X = (abm(S.get(0,x),S.pw[arg-x])+S.get(x+1,arg+1))%mod
  if arg <= y:
    Y = S.get(0,arg)
  else:
    Y = (abm(S.get(0,y),S.pw[arg-y])+S.get(y+1,arg+1))%mod  
  
  if X == Y:
    return True
  else:
    return False
        
def meguru_bisect(ng, ok,s,t):
    '''
    初期値のng,okを受け取り,is_okを満たす最小(最大)のokを返す
    まずis_okを定義すべし
    ng ok は  とり得る最小の値-1 とり得る最大の値+1
    最大最小が逆の場合はよしなにひっくり返す
    '''
    while (abs(ok - ng) > 1):
        mid = (ok + ng) // 2
        if is_ok(mid,s,t):
            ok = mid
        else:
            ng = mid
    return ok
def cmp(x,y):
  h = meguru_bisect(N,0,x,y)
  
  if h == N-1:
    return 0
  
  if x > h:
    f = A[h]
  else:
    f = A[h+1]
  if y > h:
    ff = A[h]
  else:
    ff = A[h+1]
  if f < ff:
    return -1
  elif f == ff:
    return 0
  else:
    return 1
    
  
D = list(range(N))
from functools import cmp_to_key
p = sorted(D, key = cmp_to_key(cmp))
x = p[K - 1]

ans = A[:x] + A[x + 1:]

print(*ans)
0