Cod sursa(job #2187211)

Utilizator calinfloreaCalin Florea calinflorea Data 26 martie 2018 12:14:55
Problema Gradina Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 3.51 kb
#include <bits/stdc++.h>
#define NMax 260
using namespace std;

ifstream fin ("gradina.in");
ofstream fout ("gradina.out");

struct DB
{
    int x, y;
    short ales;
};

DB a[NMax];
int n;
double ans = LONG_MAX;
char sol[NMax];
short st[NMax], st2[NMax];
bool v[NMax];
int top, top2, N, M;

inline void Citire()
{
    int i;

    fin >> n;

    for(i = 1; i <= n; i++)
    {
        fin >> a[i].x >> a[i].y;
        a[i].ales = i;
    }
}

inline bool Compar(const DB A, const DB B)
{
    if(A.y == B.y)
        return A.x < B.x;
    return A.y < B.y;
}

inline void Initializare()
{
    int i;

    for(i = 1; i <= n; i++)
        st[i] = st2[i] = v[i] = 0;
}

///determinantul pentru determinarea semiplanului
double F(int i, int j, int p)
{
    return a[p].x * (a[i].y - a[j].y) +
           a[p].y * (a[j].x - a[i].x) +
           a[i].x * a[j].y - a[j].x * a[i].y;
}

/// calculez aria celor doua poligoane
inline void Arie()
{
    double arie = 0.0, arie2 = 0.0;
    top++;
    st[top] = st[1];

    for(int i = 2; i <= top; i++)
        arie += a[st[i - 1]].x * a[st[i]].y - a[st[i - 1]].y * a[st[i]].x;

    top2++;
    st2[top2] = st2[1];

    for(int i = 2; i <= top2; i++)
        arie2 += a[st2[i - 1]].x * a[st2[i]].y - a[st2[i - 1]].y * a[st2[i]].x;

    arie = arie / 2;
    arie2 = arie2 / 2;
    if(abs(arie - arie2) < ans)
    {
        ans = abs(arie - arie2);
        for(int i = 1; i <= top; i++)
            sol[a[st[i]].ales] = 'I';

        for(int i = 1; i <= top2; i++)
            sol[a[st2[i]].ales] = 'V';
    }


}

inline void Infasuratoare()
{
    int i;
    ///planul de jos
    st[1] = 1;
    st[2] = 2;
    top = 2;
    v[1] = v[2] = 1;
    for(i = 3; i <= N; i++)
    {
        while(top > 1 && F(st[top - 1], st[top], i) < 0)
        {
            v[st[top]] = 0;
            top--;
        }
        st[++top] = i;
        v[i] = 1;
    }

    for(i = N - 1; i >= 1; i--)
        if(v[i] == 0)
        {
            while(F(st[top - 1], st[top], i) < 0)
            {
                v[st[top]] = 0;
                top--;
            }
            st[++top] = i;
            v[i] = 1;
        }

    ///planul de sus
    st2[1] = N + 1;
    st2[2] = N + 2;
    top2 = 2;

    v[N + 1] = v[N + 2] = 1;
    for(i = N + 3; i <= n; i++)
    {
        while(top > 1 && F(st2[top2 - 1], st2[top2], i) < 0)
        {
            v[st2[top2]] = 0;
            top2--;
        }
        st2[++top2] = i;
        v[i] = 1;
    }

    for(i = n - 1; i >= N + 1; i--)
        if(v[i] == 0)
        {
            while(F(st2[top2 - 1], st2[top2], i) < 0)
            {
                v[st2[top2]] = 0;
                top2--;
            }
            st2[++top2] = i;
            v[i] = 1;
        }

    for(i = 1; i <= n; i++)
       if(v[i] == 0)
           return;

    //if(top2 < 3 || top2 > (n - 3) || top < 3 || top > (n - 3))
       // return;
    Arie();
}

///impart in doua semiplane
inline void Imparte()
{
    int i, j;

    sort(a + 1, a + n + 1, Compar);

    for(j = 1; j <= n; j++)
    {
        Initializare();
        N = j;
        M = n - j;
        i = j + 1;
        while(a[j].y == a[i].y)
            N++, M--, i++;
        j = i;

        if(N < 3 || N > (n - 3))
            continue;
        Infasuratoare();
    }
}

int main()
{
    Citire();
    Imparte();
    fout << fixed << setprecision(3) << ans << "\n";
    fout << (sol + 1);
    return 0;
}