結果

問題 No.1361 [Zelkova 4th Tune *] QUADRUPLE-SEQUENCEの詩
コンテスト
ユーザー vwxyz
提出日時 2023-12-01 00:25:56
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 3,131 bytes
記録
コンパイル時間 419 ms
コンパイル使用メモリ 82,560 KB
実行使用メモリ 137,728 KB
最終ジャッジ日時 2024-09-26 14:59:19
合計ジャッジ時間 28,465 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 33 TLE * 4 -- * 37
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

import bisect
import copy
import decimal
import fractions
import heapq
import itertools
import math
import random
import sys
import time
from collections import Counter,deque,defaultdict
from functools import lru_cache,reduce
from heapq import heappush,heappop,heapify,heappushpop,_heappop_max,_heapify_max
def _heappush_max(heap,item):
    heap.append(item)
    heapq._siftdown_max(heap, 0, len(heap)-1)
def _heappushpop_max(heap, item):
    if heap and item < heap[0]:
        item, heap[0] = heap[0], item
        heapq._siftup_max(heap, 0)
    return item
from math import gcd as GCD
read=sys.stdin.read
readline=sys.stdin.readline
readlines=sys.stdin.readlines
write=sys.stdout.write
#import pypyjit
#pypyjit.set_param('max_unroll_recursion=-1')
#sys.set_int_max_str_digits(10**9)

def Bisect_Int(ok,ng,is_ok):
    while abs(ok-ng)>1:
        mid=(ok+ng)//2
        if is_ok(mid):
            ok=mid
        else:
            ng=mid
    return ok

K,L,M,N,S=map(int,readline().split())
A=list(map(int,readline().split()))
B=list(map(int,readline().split()))
C=list(map(int,readline().split()))
D=list(map(int,readline().split()))
plus0,minus0,zero0=[],[],0
for a in A:
    for b in B:
        if a*b>0:
            plus0.append(a*b)
        elif a*b<0:
            minus0.append(-a*b)
        else:
            zero0+=1
plus1,minus1,zero1=[],[],0
for c in C:
    for d in D:
        if c*d>0:
            plus1.append(c*d)
        elif c*d<0:
            minus1.append(-c*d)
        else:
            zero1+=1
minus0.sort()
minus1.sort()
plus0.sort()
plus1.sort()
def is_ok(ans):
    cnt=0
    if ans<0:
        for ab in minus0:
            cnt+=len(plus1)-bisect.bisect_left(plus1,(-ans)//ab+1)
        for ab in plus0:
            cnt+=len(minus1)-bisect.bisect_left(minus1,(-ans)//ab+1)
    elif ans==0:
        cnt+=len(plus0)*len(minus1)
        cnt+=len(plus1)*len(minus0)
    else:
        cnt+=len(plus0)*len(minus1)
        cnt+=len(plus1)*len(minus0)
        cnt+=zero0*(len(plus1)+len(minus1))+zero1*(len(plus0)+len(minus0))+zero0*zero1
        for ab in minus0:
            cnt+=bisect.bisect_right(minus1,(ans-1)//ab)
        for ab in plus0:
            cnt+=bisect.bisect_right(plus1,(ans-1)//ab)
    return cnt<S
inf=1<<60
ans=Bisect_Int(-inf,inf,is_ok)
if ans==0:
    for kk in range(K):
        if A[kk]==0:
            k,l,m,n=kk,0,0,0
            break
    else:
        for ll in range(L):
            if B[ll]==0:
                k,l,m,n=0,ll,0,0
                break
        else:
            for mm in range(M):
                if C[mm]==0:
                    k,l,m,n=0,0,mm,0
                    break
            else:
                for nn in range(N):
                    if D[nn]==0:
                        k,l,m,n=0,0,0,nn
else:
    dct={}
    for k in range(K):
        for l in range(L):
            dct[A[k]*B[l]]=(k,l)
    for m in range(M):
        for n in range(N):
            cd=C[m]*D[n]
            if cd and ans%cd==0 and ans//cd in dct:
                k,l=dct[ans//cd]
                break
        else:
            continue
        break
print(ans)
print(A[k],B[l],C[m],D[n])
0