Pagini recente » Cod sursa (job #2286027) | Rating Costea Lucian (Lucianx3) | Cod sursa (job #2030182) | Cod sursa (job #2347246) | Cod sursa (job #1065550)
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
struct Point
{
double x,y;
};
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));
}
inline bool equals(double A,double B)
{
if(fabs(A-B)<=0.00001) return 1;
return 0;
}
inline bool cmp(Point A,Point B)
{
if(equals(A.x,B.x)) return (A.y<B.y);
return (A.x<B.x);
}
void determine(Point A,Point &B,Point C,Point &D)
{
double mAC,mBD,nBD,d,a,b,c;
Point M;
M.x=(A.x+C.x)/2;
M.y=(A.y+C.y)/2;
if(A.x != C.x && A.y != C.y)
{
mAC=(A.y-C.y)/(A.x-C.x);
mBD=-1/mAC;
a=mAC;
b=-1;
c=A.y-mAC*A.x;
d=dist(A,M)*sqrt(a*a+b*b);
nBD=M.y-mBD*M.x;
B.x=(d-b*nBD-c)/(a+b*mBD);
B.y=mBD*B.x+nBD;
D.x=(-d-b*nBD-c)/(a+b*mBD);
D.y=mBD*D.x+nBD;
return;
}
if(A.x == C.x)
{
d=dist(A,M);
B.x=M.x-d;
B.y=M.y;
D.x=M.x+d;
D.y=M.y;
return;
}
if(A.y == C.y)
{
d=dist(A,M);
B.x=M.x;
B.y=M.y-d;
D.x=M.x;
D.y=M.y+d;
return;
}
}
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(cmp(A,P[md])) dr=md-1;
else st=md+1;
}
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,cmp);
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;
}