Pagini recente » Cod sursa (job #1243305) | Cod sursa (job #920520) | Cod sursa (job #1171388) | Cod sursa (job #196649) | Cod sursa (job #1884242)
#include <fstream>
#include <algorithm>
#include <iostream>
#include <vector>
#include <bitset>
#include <iomanip>
#define NMAX 255
using namespace std;
ifstream in("gradina.in");
ofstream out("gradina.out");
struct punct
{
int x,y;
int index;
} p[NMAX];
int n;
punct most_left;
inline double det(punct a, punct b,punct c)
{
return (a.x-b.x)*(b.y-c.y) - (a.y-b.y)*(b.x-c.x);
}
inline bool comparator(punct a,punct b)
{
return det(a,b,most_left) < 0;
};
inline void citire()
{
in>>n;
for(int i=0; i<n; i++)
{
in>>p[i].x>>p[i].y;
p[i].index=i;
}
}
inline vector<punct> hull(vector<punct>pct)
{
int most_left_index = 0;
//vector<punct>stiva;
static punct stiva[NMAX];
int vf = -1;
most_left = pct[0];
for(uint i=1; i<pct.size(); i++)
{
if(pct[i].x < most_left.x)
{
most_left = pct[i];
most_left_index = i;
}
}
swap(pct[most_left_index],pct[0]);
sort(pct.begin()+1,pct.end(),comparator);
//stiva.push_back(pct[0]);
// stiva.push_back(pct[1]);
stiva[++vf]=pct[0];
stiva[++vf]=pct[1];
for (uint i = 2; i < pct.size(); i++)
{
while (vf+1 >= 2 && det(stiva[vf-1], stiva[vf], pct[i]) > 0)
vf--;
//stiva.push_back(pct[i]);
stiva[++vf]=pct[i];
}
vector<punct>res;
for(int i=0;i<=vf;i++)
res.push_back(stiva[i]);
return res;
}
inline double triangle_area(punct a,punct b)
{
return (a.y+b.y)*(b.x-a.x);
}
inline double modul(double d)
{
return d>0?d:-d;
}
inline double polygon_area(vector<punct>poly)
{
double total_area = 0;
for(uint i=0; i<poly.size()-1; i++)
total_area+=triangle_area(poly[i],poly[i+1]);
total_area+=triangle_area(poly[poly.size()-1],poly[0]);
total_area/=2;
return modul(total_area);
}
double min_difference = 0x3f3f3f3f;
vector<punct>hull_ion,hull_vasile;
char who[NMAX];
int main()
{
citire();
for(int i=0; i<n; i++)
for(int j=0; j<n; j++)
{
if(i!=j)
{
vector<punct> st,dr;
for(int k=0; k<n; k++)
{
if(det(p[i],p[j],p[k])>0)
{
st.push_back(p[k]);
}
else if(det(p[i],p[j],p[k])<0)
{
dr.push_back(p[k]);
}
}
st.push_back(p[i]);
dr.push_back(p[j]);
vector<punct>hull_st = hull(st);
vector<punct>hull_dr = hull(dr);
if(hull_dr.size()==dr.size() && hull_st.size()==st.size())
{
double area_st = polygon_area(hull_st);
double area_dr = polygon_area(hull_dr);
if(modul(area_st-area_dr) < min_difference)
{
min_difference = modul(area_st-area_dr);
hull_ion = hull_dr;
hull_vasile = hull_st;
}
}
}
}
for(uint i=0; i<hull_ion.size(); i++)
who[hull_ion[i].index]='I';
for(uint i=0; i<hull_vasile.size(); i++)
who[hull_vasile[i].index]='V';
who[hull_ion.size()+hull_vasile.size()]=0;
out<<setiosflags(ios::fixed)<<setprecision(1)<<min_difference<<'\n'<<who;
return 0;
}