Cod sursa(job #1065513)

Utilizator mirceadinoMircea Popoveniuc mirceadino Data 23 decembrie 2013 14:02:35
Problema Patrate 3 Scor 65
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.46 kb
#include<cstdio>
#include<cmath>
#include<algorithm>

using namespace std;

inline bool equals(double A,double B)
{
    return fabs(A-B)<=1e-5;
}

struct Point
{
    double x,y;
    bool operator()(Point A,Point B)
    {
        if(equals(A.x,B.x)) return (A.y<B.y);
        return (A.x<B.x);
    }
};

int N,sol;
Point P[1005],A,B;

inline double dist(Point A,Point B)
{
    return sqrt((A.x-B.x)*(A.x-B.x) + (A.y-B.y)*(A.y-B.y));
}

void determine(Point A,Point &B,Point C,Point &D)
{
    B.x=A.x+A.y-C.y;
    B.y=A.y+C.x-A.x;
    D.x=C.x+A.y-C.y;
    D.y=C.y+C.x-A.x;
}

bool binary_s(Point A)
{
    int st=1,dr=N,md;
    for(;st<=dr;)
    {
        md=(st+dr)/2;
        if(equals(A.x,P[md].x) && equals(A.y,P[md].y)) return 1;
        if((A.x) < (P[md].x)) {dr=md-1; continue;}
        if((A.x) > (P[md].x)) {st=md+1; continue;}
        if(equals(A.x,P[md].x) && (A.y) < (P[md].y)) {dr=md-1; continue;}
        if(equals(A.x,P[md].x) && (A.y) > (P[md].y)) {st=md+1; continue;}

    }
    return 0;
}

int main()
{
    int i,j;
    freopen("patrate3.in","r",stdin);
    freopen("patrate3.out","w",stdout);
    scanf("%d",&N);
    for(i=1;i<=N;i++)
        scanf("%lf%lf",&P[i].x,&P[i].y);
    sort(P+1,P+N+1,Point());
    for(i=1;i<=N;i++)
        for(j=i+1;j<=N;j++)
        {
            determine(P[i],A,P[j],B);
            if(binary_s(A) && binary_s(B)) sol++;
        }
    printf("%d\n",sol/2);
    return 0;
}