Cod sursa(job #2779054)

Utilizator BlaugranasEnal Gemaledin Blaugranas Data 2 octombrie 2021 16:51:09
Problema Triang Scor 100
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 0.85 kb
#include<bits/stdc++.h>
using namespace std;
ifstream F("triang.in");
ofstream G("triang.out");
#define P pair<double,double>
#define x first
#define y second
P a[1501];
int i,n,j;
double r;
long long k;
int C(P a,P b)
{
    return a.x<b.x||(a.x==b.x&&a.y<b.y);
}
int B(int p,int q,double c,double d)
{
	int m=(p+q)/2;
	if(p>q)
      	return 0;
	if(fabs(c-a[m].x)<0.001&&fabs(d-a[m].y)<0.001)
      	return m;
    return a[m].x<c?B(m+1,q,c,d):B(p,m-1,c,d);
}
int main()
{
	F>>n;
	for(i=1;i<=n;++i)
      	F>>a[i].x>>a[i].y;
	sort(a+1,a+n+1,C),r=sqrt(3);
	for(i=1;i<n-1;++i)
        for(j=i+1;j<n;++j)
            k=k+(B(j+1,n,(a[i].x+a[j].x)/2+r*(a[i].y-a[j].y)/2,(a[i].y+a[j].y)/2+r*(a[j].x-a[i].x)/2)?1:0)+
            (B(j+1,n,(a[i].x+a[j].x)/2+r*(a[j].y-a[i].y)/2,(a[i].y+a[j].y)/2+r*(a[i].x-a[j].x)/2)?1:0);
	G<<k;
	return 0;
}