結果
問題 |
No.1143 面積Nの三角形
|
ユーザー |
![]() |
提出日時 | 2022-09-22 12:50:54 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
TLE
|
実行時間 | - |
コード長 | 1,482 bytes |
コンパイル時間 | 142 ms |
コンパイル使用メモリ | 12,544 KB |
実行使用メモリ | 23,808 KB |
最終ジャッジ日時 | 2024-12-22 04:49:00 |
合計ジャッジ時間 | 19,014 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 TLE * 1 |
other | AC * 9 TLE * 9 |
ソースコード
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 def Divisors(N): divisors=[] for i in range(1,N+1): if i**2>=N: break elif N%i==0: divisors.append(i) if i**2==N: divisors+=[i]+[N//i for i in divisors[::-1]] else: divisors+=[N//i for i in divisors[::-1]] return divisors 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 N=int(readline()) divisors=Divisors(N**2) ans=0 for y in divisors: for z in divisors: if N**2%(y*z): continue def is_ok(x): return x*(x+y+z)>=N**2//(y*z) x=Bisect_Int(1<<30,-1,is_ok) if x*(x+y+z)==N**2//(y*z): if x<=y<=z: ans+=1 print(ans)