結果

問題 No.461 三角形はいくつ?
ユーザー titiatitia
提出日時 2023-09-19 03:50:34
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 1,133 bytes
コンパイル時間 141 ms
コンパイル使用メモリ 82,168 KB
実行使用メモリ 77,792 KB
最終ジャッジ日時 2024-07-05 16:26:50
合計ジャッジ時間 20,045 ms
ジャッジサーバーID
(参考情報)
judge4 / judge1
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 45 ms
54,640 KB
testcase_01 AC 43 ms
54,384 KB
testcase_02 AC 42 ms
54,580 KB
testcase_03 AC 43 ms
54,552 KB
testcase_04 AC 57 ms
67,016 KB
testcase_05 AC 475 ms
77,468 KB
testcase_06 AC 280 ms
76,840 KB
testcase_07 AC 329 ms
76,940 KB
testcase_08 AC 251 ms
76,696 KB
testcase_09 AC 58 ms
66,892 KB
testcase_10 AC 136 ms
76,904 KB
testcase_11 AC 362 ms
77,016 KB
testcase_12 AC 210 ms
77,284 KB
testcase_13 AC 515 ms
77,256 KB
testcase_14 AC 484 ms
77,580 KB
testcase_15 AC 504 ms
77,292 KB
testcase_16 AC 79 ms
76,932 KB
testcase_17 AC 82 ms
77,260 KB
testcase_18 AC 756 ms
76,608 KB
testcase_19 AC 772 ms
77,216 KB
testcase_20 WA -
testcase_21 AC 500 ms
76,940 KB
testcase_22 WA -
testcase_23 WA -
testcase_24 WA -
testcase_25 AC 1,299 ms
76,820 KB
testcase_26 AC 1,289 ms
77,468 KB
testcase_27 AC 65 ms
70,796 KB
testcase_28 AC 66 ms
70,776 KB
testcase_29 AC 1,396 ms
76,692 KB
testcase_30 AC 1,481 ms
76,740 KB
testcase_31 AC 1,310 ms
76,880 KB
testcase_32 AC 1,251 ms
77,256 KB
testcase_33 AC 66 ms
70,808 KB
testcase_34 AC 65 ms
69,904 KB
testcase_35 AC 59 ms
68,664 KB
testcase_36 WA -
testcase_37 WA -
testcase_38 WA -
testcase_39 WA -
testcase_40 WA -
testcase_41 WA -
testcase_42 WA -
testcase_43 WA -
testcase_44 WA -
権限があれば一括ダウンロードができます

ソースコード

diff #

import sys
input = sys.stdin.readline

from bisect import bisect_left
from collections import Counter
from math import gcd

L=[[(0,1)] for i in range(3)]

N=int(input())

for i in range(N):
    p,a,b=map(int,input().split())
    G=gcd(a,b)
    L[p].append((b//G,(a+b)//G))

L[0].sort(key=lambda x:x[0]/x[1])
L[1].sort(key=lambda x:x[0]/x[1])
L[2].sort(key=lambda x:x[0]/x[1])

ANS=0
C=Counter(L[2])

for x0,x1 in L[0]:
    for y0,y1 in L[1]:
        if x0*y1+x1*y0>x1*y1:
            break

        if x0*y1>x1*y0:
            MAX=(x0,x1)
        else:
            MAX=(y0,y1)

        m0,m1=MAX

        OK=-1
        NG=len(L[2])

        while NG>OK+1:
            mid=(OK+NG)//2

            z0,z1=L[2][mid]

            if z0*m1<(m1-m0)*z1:
                OK=mid
            else:
                NG=mid

        #print(x0,x1,y0,y1,L[2],OK)

        ANS+=OK+1
        z0=x1*y1-(x0*y1+x1*y0)
        z1=x1*y1
        G=gcd(z0,z1)
        z0//=G
        z1//=G
        
        if z0*m1<(m1-m0)*z1 and (z1-z0,z1) in C:
            ANS-=C[(z0,z1)]

            #print("!",x0,x1,y0,y1,MAX,z0,z1,C[(z0,z1)])

print(ANS)


        

0