글 읽기 - 삼차방정식 풀기 2 질문이요

fruity 1년 전 2

판별식 잘 사용했고 오차범위도 작게 했는데 ㅙ 틀렸을까요?

from decimal import*
getcontext().prec = 100

p = 15

for i in range(int(input())):
    a, b, c, d = map(Decimal, input().split())

    def f(x):
        return a * x ** 3 + b * x ** 2 + c * x + d
    
    def binary(s, e):
        while e - s > 10 **- p:
            m = (s + e) / 2

            if f(m) * f(s) > 0:
                s = m
            else:
                e = m
        return s
    
    if b ** 2 - 3 * a * c < 0:
        ans = [binary(Decimal(-1e+7), Decimal(1e+7))]
    else:
        x1 = (-b - (b ** 2 - 3 * a * c) ** Decimal('0.5')) / (3 * a); x2 = (-b + (b ** 2 - 3 * a * c) ** Decimal('0.5')) / (3 * a)
        
        if round(f(x1) * f(x2), p) > 0:
            ans = [binary(Decimal(-1e+7), Decimal(1e+7))]
        elif round(f(x1), p) == 0 or round(f(x2), p) == 0:
            if round(f(x1), p) == 0:
                if round(x1, p) == round(x2, p):
                    ans = [x1]
                else:
                    ans = [x1, binary(x2, Decimal(1e+7))]
            else:
                ans = [binary(Decimal(-1e+7), x1), x2]
        else:
            ans = [binary(Decimal(-1e+7), x1), binary(x1, x2), binary(x2, Decimal(1e+7))]
    ans = ['{:.15f}'.format(i) for i in ans]
    print(*ans)

jongpark1234 1년 전 2

p값을 조금 더 높여야 하고 x1이 x2보다 더 클 때 두 값을 바꿔주어야 합니다.

fruity 1년 전 2

감사합니다!

댓글을 작성하려면 로그인해야 합니다.

문제

문제

출처

ICPC

삼차방정식 풀기 2 질문이요

fruity 1년 전 2

jongpark1234 1년 전 2

fruity 1년 전 2

Baekjoon Online Judge

채점 현황

문제

유저 대회

고등학교 대회

출처

대학교 대회

도움말