結果

問題 No.1361 [Zelkova 4th Tune *] QUADRUPLE-SEQUENCEの詩
ユーザー 👑 KazunKazun
提出日時 2020-09-14 23:07:55
言語 PyPy3
(7.3.13)
結果
AC  
実行時間 896 ms / 2,000 ms
コード長 3,587 bytes
コンパイル時間 470 ms
コンパイル使用メモリ 87,152 KB
実行使用メモリ 182,716 KB
最終ジャッジ日時 2023-08-27 16:12:48
合計ジャッジ時間 29,141 ms
ジャッジサーバーID
(参考情報)
judge13 / judge15
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 65 ms
71,172 KB
testcase_01 AC 65 ms
71,244 KB
testcase_02 AC 69 ms
71,276 KB
testcase_03 AC 65 ms
71,012 KB
testcase_04 AC 65 ms
71,188 KB
testcase_05 AC 66 ms
71,376 KB
testcase_06 AC 62 ms
71,248 KB
testcase_07 AC 65 ms
71,344 KB
testcase_08 AC 81 ms
76,224 KB
testcase_09 AC 77 ms
76,172 KB
testcase_10 AC 82 ms
76,256 KB
testcase_11 AC 81 ms
76,276 KB
testcase_12 AC 79 ms
75,856 KB
testcase_13 AC 80 ms
76,480 KB
testcase_14 AC 78 ms
75,668 KB
testcase_15 AC 86 ms
76,308 KB
testcase_16 AC 85 ms
76,116 KB
testcase_17 AC 78 ms
76,320 KB
testcase_18 AC 129 ms
78,740 KB
testcase_19 AC 102 ms
77,412 KB
testcase_20 AC 115 ms
78,556 KB
testcase_21 AC 101 ms
78,844 KB
testcase_22 AC 110 ms
78,428 KB
testcase_23 AC 123 ms
79,604 KB
testcase_24 AC 128 ms
80,312 KB
testcase_25 AC 126 ms
79,356 KB
testcase_26 AC 131 ms
79,520 KB
testcase_27 AC 109 ms
77,832 KB
testcase_28 AC 389 ms
110,568 KB
testcase_29 AC 326 ms
104,776 KB
testcase_30 AC 337 ms
99,936 KB
testcase_31 AC 431 ms
110,508 KB
testcase_32 AC 549 ms
134,276 KB
testcase_33 AC 334 ms
98,796 KB
testcase_34 AC 393 ms
120,348 KB
testcase_35 AC 423 ms
114,552 KB
testcase_36 AC 620 ms
137,168 KB
testcase_37 AC 220 ms
98,208 KB
testcase_38 AC 342 ms
109,620 KB
testcase_39 AC 282 ms
95,288 KB
testcase_40 AC 290 ms
103,212 KB
testcase_41 AC 454 ms
116,160 KB
testcase_42 AC 671 ms
135,824 KB
testcase_43 AC 403 ms
111,164 KB
testcase_44 AC 526 ms
117,924 KB
testcase_45 AC 351 ms
111,732 KB
testcase_46 AC 547 ms
124,240 KB
testcase_47 AC 330 ms
99,496 KB
testcase_48 AC 605 ms
182,440 KB
testcase_49 AC 384 ms
182,448 KB
testcase_50 AC 581 ms
182,716 KB
testcase_51 AC 476 ms
182,692 KB
testcase_52 AC 588 ms
182,624 KB
testcase_53 AC 885 ms
156,904 KB
testcase_54 AC 858 ms
156,808 KB
testcase_55 AC 896 ms
156,724 KB
testcase_56 AC 896 ms
156,740 KB
testcase_57 AC 714 ms
156,592 KB
testcase_58 AC 606 ms
153,872 KB
testcase_59 AC 587 ms
153,932 KB
testcase_60 AC 266 ms
153,976 KB
testcase_61 AC 268 ms
154,008 KB
testcase_62 AC 272 ms
153,976 KB
testcase_63 AC 120 ms
127,112 KB
testcase_64 AC 119 ms
127,180 KB
testcase_65 AC 117 ms
127,148 KB
testcase_66 AC 120 ms
127,080 KB
testcase_67 AC 121 ms
127,028 KB
testcase_68 AC 65 ms
71,372 KB
testcase_69 AC 75 ms
71,388 KB
testcase_70 AC 66 ms
71,236 KB
testcase_71 AC 64 ms
71,500 KB
testcase_72 AC 62 ms
71,432 KB
testcase_73 AC 62 ms
71,476 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

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 なる p in G,q in H を求める.
def h(x,G,H):
    if x==0:
        if 0 in G:
            return (0,H[0])
        else:
            return (G[0],0)

    H=set(H)
    for g in G:
        if g==0:
            continue

        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))

#================================================
# T=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)
0