結果
問題 | No.550 夏休みの思い出(1) |
ユーザー |
![]() |
提出日時 | 2021-11-12 16:04:56 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 135 ms / 2,000 ms |
コード長 | 1,436 bytes |
コンパイル時間 | 85 ms |
コンパイル使用メモリ | 12,800 KB |
実行使用メモリ | 10,880 KB |
最終ジャッジ日時 | 2024-11-25 03:31:47 |
合計ジャッジ時間 | 8,947 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 55 |
ソースコード
import sys # sys.setrecursionlimit(200005) int1 = lambda x: int(x)-1 p2D = lambda x: print(*x, sep="\n") def II(): return int(sys.stdin.readline()) def LI(): return list(map(int, sys.stdin.readline().split())) def LLI(rows_number): return [LI() for _ in range(rows_number)] def LI1(): return list(map(int1, sys.stdin.readline().split())) def LLI1(rows_number): return [LI1() for _ in range(rows_number)] def SI(): return sys.stdin.readline().rstrip() # dij = [(0, 1), (-1, 0), (0, -1), (1, 0)] # dij = [(0, 1), (-1, 0), (0, -1), (1, 0), (1, 1), (1, -1), (-1, 1), (-1, -1)] # inf = 18446744073709551615 inf = 4294967295 md = 10**9+7 # md = 998244353 def factors(a): if a <= 0: return [0] dd0, dd1 = [], [] for d in range(1, min(round(a**0.5), 10**6)+1): if a%d: continue dd0.append(d) return dd0 def PolDiv(aa, v): res = [aa[-1]]*(len(aa)-1) for i in range(len(aa)-3, -1, -1): res[i] = res[i+1]*v+aa[i+1] rem = res[0]*v+aa[0] return res, rem a, b, c = LI() ans = [] if c == 0: a, b, c = 1, a, b ans.append(0) else: ff = factors(abs(c)) ff += [-f for f in ff] for f in ff: if f**3+a*f**2+b*f+c == 0: ans.append(f) break abc, _ = PolDiv([c, b, a, 1], ans[0]) c, b, a = abc # print(ans, a, b, c) # print(b**2-4*a*c) d = round((b**2-4*a*c)**0.5) ans.append((-b+d)//a//2) ans.append((-b-d)//a//2) ans.sort() print(*ans)