結果
問題 |
No.1361 [Zelkova 4th Tune *] QUADRUPLE-SEQUENCEの詩
|
ユーザー |
![]() |
提出日時 | 2023-12-01 01:04:12 |
言語 | PyPy3 (7.3.15) |
結果 |
TLE
|
実行時間 | - |
コード長 | 2,497 bytes |
コンパイル時間 | 380 ms |
コンパイル使用メモリ | 82,560 KB |
実行使用メモリ | 118,144 KB |
最終ジャッジ日時 | 2024-09-26 15:00:32 |
合計ジャッジ時間 | 17,128 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 28 TLE * 1 -- * 45 |
ソースコード
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())) lst0=[] for a in A: for b in B: lst0.append(a*b) lst1=[] for c in C: for d in D: lst1.append(c*d) lst0.sort() def is_ok(ans): cnt=0 for cd in lst1: if cd<0: cnt+=len(lst0)-Bisect_Int(0,len(lst0)+1,lambda i:i==0 or lst0[i-1]*cd>=ans) elif cd==0: if 0<ans: cnt+=len(lst0) else: cnt+=Bisect_Int(0,len(lst0)+1,lambda i:i==0 or lst0[i-1]*cd<ans) 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])