結果
| 問題 | No.1361 [Zelkova 4th Tune *] QUADRUPLE-SEQUENCEの詩 |
| ユーザー |
👑  Kazun
|
| 提出日時 | 2020-09-14 23:04:11 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
RE
|
| 実行時間 | - |
| コード長 | 3,447 bytes |
| コンパイル時間 | 252 ms |
| コンパイル使用メモリ | 86,876 KB |
| 実行使用メモリ | 182,868 KB |
| 最終ジャッジ日時 | 2023-08-27 16:13:32 |
| 合計ジャッジ時間 | 30,273 ms |
|
ジャッジサーバーID （参考情報） |
judge13 / judge14 |
（要ログイン）
テストケース
テストケース表示| 入力 | 結果 | 実行時間 実行使用メモリ |
|---|---|---|
| testcase_00 | AC | 64 ms
71,352 KB |
| testcase_01 | AC | 65 ms
71,316 KB |
| testcase_02 | AC | 64 ms
71,176 KB |
| testcase_03 | AC | 65 ms
71,016 KB |
| testcase_04 | AC | 64 ms
71,240 KB |
| testcase_05 | AC | 66 ms
71,280 KB |
| testcase_06 | AC | 64 ms
71,384 KB |
| testcase_07 | AC | 65 ms
71,188 KB |
| testcase_08 | AC | 82 ms
75,932 KB |
| testcase_09 | AC | 74 ms
75,924 KB |
| testcase_10 | AC | 83 ms
76,192 KB |
| testcase_11 | AC | 79 ms
76,172 KB |
| testcase_12 | AC | 80 ms
76,344 KB |
| testcase_13 | AC | 80 ms
76,408 KB |
| testcase_14 | AC | 75 ms
75,656 KB |
| testcase_15 | AC | 83 ms
76,228 KB |
| testcase_16 | AC | 84 ms
76,184 KB |
| testcase_17 | AC | 81 ms
76,268 KB |
| testcase_18 | AC | 128 ms
78,928 KB |
| testcase_19 | AC | 104 ms
77,348 KB |
| testcase_20 | AC | 122 ms
78,536 KB |
| testcase_21 | AC | 100 ms
78,488 KB |
| testcase_22 | AC | 108 ms
78,500 KB |
| testcase_23 | AC | 126 ms
79,648 KB |
| testcase_24 | AC | 127 ms
80,576 KB |
| testcase_25 | AC | 128 ms
79,644 KB |
| testcase_26 | AC | 132 ms
80,108 KB |
| testcase_27 | AC | 116 ms
77,640 KB |
| testcase_28 | RE | - |
| testcase_29 | AC | 329 ms
104,404 KB |
| testcase_30 | AC | 350 ms
99,592 KB |
| testcase_31 | AC | 442 ms
110,196 KB |
| testcase_32 | AC | 554 ms
134,272 KB |
| testcase_33 | AC | 326 ms
98,760 KB |
| testcase_34 | AC | 387 ms
120,372 KB |
| testcase_35 | AC | 418 ms
114,484 KB |
| testcase_36 | AC | 631 ms
137,436 KB |
| testcase_37 | AC | 223 ms
97,960 KB |
| testcase_38 | AC | 330 ms
109,624 KB |
| testcase_39 | AC | 280 ms
95,368 KB |
| testcase_40 | AC | 286 ms
103,292 KB |
| testcase_41 | AC | 466 ms
114,876 KB |
| testcase_42 | AC | 677 ms
135,900 KB |
| testcase_43 | AC | 405 ms
111,208 KB |
| testcase_44 | AC | 530 ms
117,860 KB |
| testcase_45 | AC | 367 ms
111,732 KB |
| testcase_46 | AC | 559 ms
124,632 KB |
| testcase_47 | AC | 323 ms
99,436 KB |
| testcase_48 | AC | 595 ms
182,868 KB |
| testcase_49 | AC | 390 ms
182,728 KB |
| testcase_50 | AC | 584 ms
182,744 KB |
| testcase_51 | AC | 478 ms
182,644 KB |
| testcase_52 | AC | 601 ms
182,700 KB |
| testcase_53 | AC | 883 ms
156,800 KB |
| testcase_54 | AC | 873 ms
156,596 KB |
| testcase_55 | AC | 874 ms
156,864 KB |
| testcase_56 | AC | 884 ms
156,488 KB |
| testcase_57 | AC | 711 ms
156,672 KB |
| testcase_58 | AC | 589 ms
153,764 KB |
| testcase_59 | AC | 595 ms
153,824 KB |
| testcase_60 | AC | 269 ms
153,648 KB |
| testcase_61 | AC | 278 ms
153,944 KB |
| testcase_62 | AC | 267 ms
153,936 KB |
| testcase_63 | RE | - |
| testcase_64 | RE | - |
| testcase_65 | RE | - |
| testcase_66 | RE | - |
| testcase_67 | RE | - |
| testcase_68 | RE | - |
| testcase_69 | RE | - |
| testcase_70 | RE | - |
| testcase_71 | AC | 63 ms
71,532 KB |
| testcase_72 | AC | 65 ms
71,388 KB |
| testcase_73 | RE | - |
ソースコード
def General_Binary_Increase_Search(L,R,cond,Integer=True,ep=1/(1<<20)):
"""条件式が単調増加であるとき,一般的な二部探索を行う.
L:解の下限
R:解の上限
cond:条件(1変数関数,広義単調減少 or 広義単調減少を満たす)
Integer:解を整数に制限するか?
ep:Integer=Falseのとき,解の許容する誤差
"""
if not(cond(R)):
return False
if Integer:
R+=1
while R-L>1:
C=L+(R-L)//2
if cond(C):
R=C
else:
L=C
return R
else:
while (R-L)>=ep:
C=L+(R-L)/2
if cond(C):
R=C
else:
L=C
return R
#================================================
def f(x): #負の時の判定
Z=0
I=0
for u in U_pos:
if I<V_negative:
while u*V_neg[I]<=x:
I+=1
if I==V_negative:
break
Z+=I
I=0
for v in V_pos:
if I<U_negative:
while v*U_neg[I]<=x:
I+=1
if I==U_negative:
break
Z+=I
return Z
def g(x): #正の時の判定
Z=0
I=V_positive
for u in U_pos:
if I>0:
while u*V_pos[I-1]>x:
I-=1
if I==0:
break
Z+=I
I=0
for v in V_neg:
if I<U_negative:
while v*U_neg[-(I+1)]<=x:
I+=1
if I==U_negative:
break
Z+=I
return Z
#================================================
# pq=x なる x in G,y in H を求める.
def h(x,G,H):
H=set(H)
for g in G:
if (x%g==0) and (x//g in H):
return (g,x//g)
return None
#================================================
#入力
K,L,M,N,S=map(int,input().split())
A=list(map(int,input().split()))
B=list(map(int,input().split()))
C=list(map(int,input().split()))
D=list(map(int,input().split()))
#================================================
# (A,B),(C,D)の2つに分けて考える.
U=[a*b for a in A for b in B]
U.sort()
V=[c*d for c in C for d in D]
V.sort()
U_pos=[u for u in U if u>0]
U_neg=[u for u in U if u<0]
V_pos=[v for v in V if v>0]
V_neg=[v for v in V if v<0]
U_positive=len(U_pos)
U_negative=len(U_neg)
U_zero=K*L-(U_positive+U_negative)
V_positive=len(V_pos)
V_negative=len(V_neg)
V_zero=M*N-(V_positive+V_negative)
#================================================
# Eの正,ゼロ,負の個数を求める
E_positive=U_positive*V_positive+U_negative*V_negative
E_negative=U_positive*V_negative+U_negative*V_positive
E_zero=K*L*M*N-(E_positive+E_negative)
#================================================
# Jを求める.
U_abs_max=abs(max(U,key=lambda u:abs(u)))
V_abs_max=abs(max(V,key=lambda v:abs(v)))
Abs_max=U_abs_max*V_abs_max+1
if S<=E_negative: #負確定
Ans=General_Binary_Increase_Search(-Abs_max,0,lambda x:f(x)>=S)
elif E_negative+1<=S<=E_negative+E_zero: #ゼロ確定
Ans=0
else: #正確定
Ans=General_Binary_Increase_Search(0,Abs_max,lambda x:g(x)>=S-(E_negative+E_zero))
#================================================
# J=abcd なる a,b,c,dを求める.
alpha,beta=h(Ans,U,V)
a,b=h(alpha,A,B)
c,d=h(beta ,C,D)
#================================================
#出力
print(Ans)
print(a,b,c,d)