Cod sursa(job #170567)

Utilizator mihai_floreaFlorea Mihai Alexandru mihai_florea Data 2 aprilie 2008 21:54:54
Problema Triang Scor 100
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.64 kb
#include <cstdio>
const double cos60=0.500000000;
const double sin60=0.866025403;
int n,nr;
double x[1501],y[1501];
double dabs(double x){
       return (x>=0?x:-x);
       }
int egal (double a,double b){
    return (dabs(a-b)<0.001);
    }
void qsort(int ls,int ld){
     int i=ls,j=ld,mij=(ls+ld)/2;
     double xx=x[mij],yy=y[mij];
     do {while (x[i]<xx || (egal(x[i],xx) && y[i]<yy)) i++;
         while (x[j]>xx || (egal(x[j],xx) && y[j]>yy)) j--;
         if (i<=j) {double aux=x[i];x[i]=x[j];x[j]=aux;
                    aux=y[i];y[i]=y[j];y[j]=aux;
                    i++;j--;}
         }
       while (i<=j);
     if (j>ls) qsort(ls,j);
     if (i<ld) qsort(i,ld);
     }
int gasit(int ls,int ld,double xx,double yy){
    int mij;
    while (ls<=ld) {
          mij=(ls+ld)/2;
          if (egal(x[mij],xx) && egal(y[mij],yy)) return mij;
          if (x[mij]>xx || (egal(x[mij],xx) && y[mij]>yy)) ld=mij-1;
                        else ls=mij+1;                               
          }
    return 0;
    }
int main(){
    int i,j;
    double xx,yy;
    freopen("triang.in","r",stdin);
    freopen("triang.out","w",stdout);
    scanf("%d",&n);
    for (i=1;i<=n;i++) scanf("%lf %lf",&x[i],&y[i]);
    qsort(1,n);
    nr=0;
    for (i=1;i<=n-1;i++) 
      for (j=i+1;j<=n;j++)
       {xx=x[i]+cos60*(x[j]-x[i])+sin60*(y[i]-y[j]);
        yy=y[i]+cos60*(y[j]-y[i])-sin60*(x[i]-x[j]);
        if (gasit(j+1,n,xx,yy)) nr++;
        xx=x[j]+cos60*(x[i]-x[j])-sin60*(y[i]-y[j]);
        yy=y[j]+cos60*(y[i]-y[j])+sin60*(x[i]-x[j]);
        if (gasit(j+1,n,xx,yy)) nr++;
        }
    printf("%d",nr);
   }