Pagini recente » Cod sursa (job #2850203) | Cod sursa (job #2490623) | Cod sursa (job #1495284) | Cod sursa (job #2331378) | Cod sursa (job #2187211)
#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;
}