Pagini recente » Cod sursa (job #1523049) | Cod sursa (job #488632) | Cod sursa (job #1621881) | Cod sursa (job #892193) | Cod sursa (job #1020325)
#include <cstdio>
#include <math.h>
#include <algorithm>
#include <functional>
#include <vector>
#define pdd pair<double,double>
#define M_PI 3.14159265358979323846
#define x first
#define y second
using namespace std;
const double eps = 1e-3;
const double s60 = sin(60 * M_PI / 180.0);
const double c60 = cos(60 * M_PI / 180.0);
const int nmax = 1502;
int n;
pdd v[nmax];
inline pair<pdd,pdd> getPoints(pdd a,pdd b) {
return make_pair(
make_pair(c60 * (a.x - b.x) - s60 * (a.y - b.y) + b.x,
s60 * (a.x - b.x) + c60 * (a.y - b.y) + b.y),
make_pair(c60 * (a.x - b.x) + s60 * (a.y - b.y) + b.x,
- s60 * (a.x - b.x) + c60 * (a.y - b.y) + b.y ));
}
inline int comp(const double &a,const double &b) {
if(fabs(b - a) < eps) return 0;
return a < b ? -1 : 1;
}
inline void readData() {
scanf("%d",&n);
for(int i = 0;i < n;i++) {
scanf("%lf %lf",&v[i].x,&v[i].y);
}
}
inline bool search(pdd p) {
int pos = 0;
for(int step = 1<<11;step > 0;step >>= 1) {
if(pos + step < n) {
if(comp(p.x,v[pos + step].x) < 1 && comp(p.y,v[pos + step].y) < 1) {
pos += step;
}
}
}
return (comp(p.x,v[pos].x) == 0 && comp(p.y,v[pos].y) == 0 );
}
inline int solve() {
int ret = 0;
sort(v,v + n);
for(int i = 0;i < n;i++) {
for(int j = i + 1;j < n;j++) {
pair< pdd ,pdd > points = getPoints(v[i],v[j]);
ret += search(points.x);
ret += search(points.y);
}
}
return ret;
}
int main()
{
freopen("triang.in","r",stdin);
freopen("triang.out","w",stdout);
readData();
printf("%d\n",solve());
return 0;
}