Cod sursa(job #1598509)

Utilizator AeroHHorea Stefan AeroH Data 12 februarie 2016 23:05:30
Problema Rubarba Scor 70
Compilator cpp Status done
Runda Arhiva de probleme Marime 0.97 kb
#include <bits/stdc++.h>
#define x first
#define y second
#define eps 0.000000000000001//10^-14
#define ld long double
using namespace std;
string z="rubarba.";
ifstream f(z+"in");
ofstream g(z+"out");

int i,n;
ld x,y,s,c,mx,Mx,my,My,m,l=0,r=0.78539816339744830962 * 2;
pair<ld,ld> v[100002];

ld e(ld u)
{
    mx=my=1e9;
    Mx=My=-1e9;
    for(int i=1;i<=n;++i)
    {
        s=sin(u),c=cos(u);
        x = v[i].x * c - v[i].y * s;
        y = v[i].x * s + v[i].y * c;
        mx=min(mx,x);
        Mx=max(Mx,x);
        my=min(my,y);
        My=max(My,y);
    }
    return (Mx-mx)*(My-my);
}

double ts(double l, double r)
{
    if (r - l < eps)
        return l;

    double m1 = l + (r - l) / 3.0;
    double m2 = r - (r - l) / 3.0;

    if (e(m1) < e(m2))
        return ts(l, m2);
    else
        return ts(m1, r);
}

int main()
{
    f>>n;
    for(i=1;i<=n;++i)
        f>>v[i].x>>v[i].y;
    g<<fixed<<setprecision(3)<<e(ts(l,r));
}