E081. Нахождение площади простого многоугольника

e-maxx algorithm original: C/C++ #algorithm #emaxx #geometry #math

Der Aufgabentext wird für die gewählte Sprache aus dem Russischen übersetzt. Code bleibt unverändert.

Static

Источник: e-maxx.ru/algo, страница PDF 236.

Пусть дан простой многоугольник (т.е. без самопересечений, но не обязательно выпуклый), заданный координатами своих вершин в порядке обхода по или против часовой стрелки. it is required find его площадь.

Способ 1

Это легко сделать, если перебрать все рёбра и сложить площади трапеций, ограниченных каждым edgeм. Площадь нужно брать с тем знаком, с каким она получится (именно благодаря знаку вся "лишняя" площадь сократится). Т.е. формула такова:

S += (X2 - X1) * (Y1 + Y2) / 2

Код:

double sq (const vector<point> & fig)

{

double res = 0;

for (unsigned i=0; i<fig.size(); i++)

{

point

p1 = i ? fig[i-1] : fig.back(),

p2 = fig[i];

res += (p1.x - p2.x) * (p1.y + p2.y);

}

return fabs (res) / 2;

}

Способ 2

Можно поступить другим образом. Выберем произвольно точку O, переберём все рёбра, прибавляя к ответу ориентированную площадь треугольника, образованного edgeм и точкой O (см. Ориентированная площадь треугольника). Опять же, благодаря знаку, вся лишняя площадь сократится, и останется только ответ. Этот способ хорош тем, что его проще обобщить на более сложные случаи (наBeispiel, когда некоторые стороны - не прямые, а дуги окружности).

C# Lösung

Auto-Entwurf, vor dem Einreichen prüfen

using System;
using System.Collections.Generic;
using System.Linq;

public static class AlgorithmDraft
{
    // Auto-generated C# draft from the original e-maxx C/C++ listing. Review before production use.
    S += (X2 - X1) * (Y1 + Y2) / 2
    double sq (const vector<point> & fig)
    {
            double res = 0;
            for (unsigned i=0; i<fig.size(); i++)
            {
                    point
                            p1 = i ? fig[i-1] : fig.back(),
                            p2 = fig[i];
                    res += (p1.x - p2.x) * (p1.y + p2.y);
            }
            return fabs (res) / 2;
    }
}

C++ Lösung

zugeordnet/original

S += (X2 - X1) * (Y1 + Y2) / 2
double sq (const vector<point> & fig)
{
        double res = 0;
        for (unsigned i=0; i<fig.size(); i++)
        {
                point
                        p1 = i ? fig[i-1] : fig.back(),
                        p2 = fig[i];
                res += (p1.x - p2.x) * (p1.y + p2.y);
        }
        return fabs (res) / 2;
}

Java Lösung

Auto-Entwurf, vor dem Einreichen prüfen

import java.util.*;
import java.math.*;

public class AlgorithmDraft {
    // Auto-generated Java draft from the original e-maxx C/C++ listing. Review before production use.
    S += (X2 - X1) * (Y1 + Y2) / 2
    double sq (const vector<point> & fig)
    {
            double res = 0;
            for (unsigned i=0; i<fig.size(); i++)
            {
                    point
                            p1 = i ? fig[i-1] : fig.back(),
                            p2 = fig[i];
                    res += (p1.x - p2.x) * (p1.y + p2.y);
            }
            return fabs (res) / 2;
    }
}

Материал разбит как Algorithmusическая Aufgabe: изучить постановку, понять асимптотику и реализовать Algorithmus на выбранном языке.

Stellen zu dieser Aufgabe

aktive Stellen with overlapping task tags are angezeigt.

Alle Stellen

Es gibt noch keine aktiven Stellen.