F. Половина многоугольника.py

import math

RADIUS = 100
TOL = 1


def polygon(n):
    result = []
    for i in range(n):
        angle = 2 * math.pi * i / n
        x = math.cos(angle) * RADIUS
        y = math.sin(angle) * RADIUS
        result.append((
            int(round(x)),
            int(round(y))
        ))
    return result


def calc_angles(points, n):
    angles = []
    for i in range(n):
        p1 = points[i]
        p2 = points[(i + 1) % n]
        p3 = points[(i + 2) % n]
        x1 = p2[0] - p1[0]
        x2 = p3[0] - p2[0]
        y1 = p2[1] - p1[1]
        y2 = p3[1] - p2[1]
        v1 = math.sqrt(x1 ** 2 + y1 ** 2)
        v2 = math.sqrt(x2 ** 2 + y2 ** 2)
        angle = math.acos((x1 * x2 + y1 * y2) / (v1 * v2))
        angles.append(180 - angle * 180 / math.pi)
    return angles


def start1():
    n = int(input())
    for x, y in POINTS[n]:
        print(x, y)


def start2():
    n2 = int(input())
    points2 = [
        list(map(float, input().split()))
        for _ in range(n2)
    ]
    angles = calc_angles(points2, n2)
    print(angles)
    max_sovpad = -1
    answer = -1
    for n in range(n2 * 2 - 5, n2 * 2 - 1):
        if n not in ANGLES:
            continue
        sovp = 0
        for angle in angles:
            for angle2 in ANGLES[n]:
                if abs(angle - angle2) < TOL:
                    sovp += 1
                    break
        if sovp > max_sovpad:
            max_sovpad = sovp
            answer = n
    print(answer)


def main():
    global POINTS, ANGLES
    POINTS = {
        i: polygon(i)
        for i in range(10, 20 + 1)
    }
    ANGLES = {
        i: calc_angles(POINTS[i], i)
        for i in range(10, 20 + 1)
    }
    s = input()
    t = int(input())
    for _ in range(t):
        if s == 'draw':
            start1()
        else:
            start2()


main()