글 읽기 - 스페셜 저지와 질문을 수정해 주세요

thxios457 1년 전 0

https://www.acmicpc.net/board/...

주어진 예제에서 소수점 아래 5번째 자리까지 값이 주어졌는데,

문제에서 요구하는 것은 절대오차 1e-6이라 푸는 것이 불가능합니다.

예제의 답 또한 실제 값과 차이가 있기에 질문을 작성하려다

이전에 저와 같은 질문을 쓰신 분이 있어 첨부합니다. (아마 카테고리가 잘못 설정되었던 것 같습니다)

import sys
from decimal import Decimal, getcontext
sys_input = sys.stdin.readline
sys_print = sys.stdout.write
cnt = 1

getcontext().prec = 40
ftod = lambda x: Decimal(str(x)+'0'*50)


pi = Decimal('3.14159265358979323846264338327950288419716939937510')

def d_atan2(y, x):
    if x == 0:
        return pi/2 if y > 0 else -pi/2
    return d_atan(y/x)

def d_atan(x):
    getcontext().prec += 2
    if x > 1 or x < -1:
        lasts, s, fact, num, sign = -1, 0, 1, 1/x, -1
        while s != lasts:
            lasts = s
            s += num / fact * sign
            num *= 1 / (x * x)
            sign *= -1
            fact += 2
        if x > 1:
            s += pi/2
        else:
            s -= pi/2
    else:
        lasts, s, fact, num, sign = 0, x, 1, x, 1
        while s != lasts:
            lasts = s
            fact += 2
            num *= x * x
            sign *= -1
            s += num / fact * sign
    getcontext().prec -= 2
    return +s

def d_sin(x):
    x = (x+pi) % (2*pi) - pi
    getcontext().prec += 2
    i, lasts, s, fact, num, sign = 1, 0, x, 1, x, 1
    while s != lasts:
        lasts = s
        i += 2
        fact *= i * (i-1)
        num *= x * x
        sign *= -1
        s += num / fact * sign
    getcontext().prec -= 2
    return +s

def d_cos(x):
    x = (x+pi) % (2*pi) - pi
    getcontext().prec += 2
    i, lasts, s, fact, num, sign = 0, 0, 1, 1, 1, 1
    while s != lasts:
        lasts = s
        i += 2
        fact *= i * (i-1)
        num *= x * x
        sign *= -1
        s += num / fact * sign
    getcontext().prec -= 2
    return +s

while True:
    n = int(sys_input())
    if n == 0:
        break

    inner_angle = 2 * pi / n

    x1, y1 = map(ftod, sys_input().split())
    x2, y2 = map(ftod, sys_input().split())
    x3, y3 = map(ftod, sys_input().split())

    sq1 = -(x1*x1 + y1*y1)
    sq2 = -(x2*x2 + y2*y2)
    sq3 = -(x3*x3 + y3*y3)

    a, b, s1 = x1-x3, y1-y3, sq1-sq3
    c, d, s2 = x2-x3, y2-y3, sq2-sq3

    # det_inv = 1 / (a*d - b*c)
    # A = det_inv * (d*s1 - b*s2)
    # B = det_inv * (-c*s1 + a*s2)

    det = (a*d - b*c)
    A = (d*s1 - b*s2) / det
    B = (-c*s1 + a*s2) / det

    C = sq1 - A*x1 - B*y1

    p = -A / 2
    q = -B / 2
    r_sq = p * p + q * q - C

    # p = Decimal(round(p))
    # q = Decimal(round(q))
    # r_sq = Decimal(round(r_sq))

    cur_angle = d_atan2(y1-q, x1-p)
    x_max, x_min = -1, 1
    y_max, y_min = -1, 1

    for i in range(n):
        yy, xx = d_sin(cur_angle), d_cos(cur_angle)

        if yy > y_max:
            y_max = yy
        if yy < y_min:
            y_min = yy
        if xx > x_max:
            x_max = xx
        if xx < x_min:
            x_min = xx

        cur_angle += inner_angle
        if cur_angle > pi:
            cur_angle -= 2*pi

    ans = (y_max - y_min)*(x_max - x_min)*r_sq

    print(f'Polygon {cnt}: {ans:.9f}')
    cnt += 1

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

문제

문제

출처

ICPC

스페셜 저지와 질문을 수정해 주세요

thxios457 1년 전 0

Baekjoon Online Judge

채점 현황

문제

유저 대회

고등학교 대회

출처

대학교 대회

도움말