Cod sursa(job #1828379)

Utilizator Alexa2001Alexa Tudose Alexa2001 Data 13 decembrie 2016 10:45:24
Problema Laser Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.46 kb
#include <bits/stdc++.h>

using namespace std;

const double eps = 1e-8, PI = 3.14159265358979323846264338327950288;
const int Nmax = 515;

int X1, Y1, X2, Y2, i, n, j, p[Nmax], xr;
bitset<Nmax> a[Nmax];
double angle1, angle2;
pair<double, double> A[Nmax];
vector<double> ans;

double get_angle(int x, int y)
{
    double angle = atan2(y,x);
    if(angle < 0) angle += 2*PI;
    return angle;
}

bool intersect(double x, double y, double z)
{
    if(y-x > PI) return (z<x+eps || z>y-eps);
    return (x<z+eps && z-eps<y);
}

int main()
{
    freopen("laser.in", "r", stdin);
    freopen("laser.out", "w", stdout);

    scanf("%d", &n);

    for(i=1; i<=n; ++i)
    {
        scanf("%d%d%d%d", &X1, &Y1, &X2, &Y2);

        angle1 = get_angle(X1, Y1);
        angle2 = get_angle(X2, Y2);
        if(angle1 > angle2) swap(angle1, angle2);
        A[i] = {angle1, angle2};
    }

    for(i=1; i<=n; ++i)
    {
        for(j=1; j<=n; ++j)
            a[i][j] = intersect(A[i].first, A[i].second, A[j].first);
        scanf("%d", &xr);
        a[i][n+1] = xr;
    }

    for(i=1; i<=n; ++i)
    {
        for(j=1; j<=n+1; ++j)
            if(a[i][j]) break;

        if(j==n+2) continue;
        p[i] = j;

        for(j=1; j<=n; ++j)
        if(i!=j && a[j][p[i]])
            a[j] ^= a[i];
    }

    for(i=1; i<=n; ++i)
    if(p[i] && a[i][n+1])
        ans.push_back(A[p[i]].first);

    printf("%d\n", ans.size());
    for(auto x : ans)
        printf("%.6lf\n", x/PI*180 + eps);

    return 0;
}