重温经典（xdoj1153）

省赛结束了，感觉自己还有很长的路要走。回忆过去，从我第一次接触acm，到现在不觉已经一年多了。突然找到去年第一次参加校赛时的一道未ac之题，当时可是困扰了我很久呢，现在重新拿来做了一下，算是重温经典了吧。

 

万神的线段

 

思路：寻找平行直线对数，说白了，就是利用斜率关系，寻找出相同斜率的直线，然后，根据n个对象之间可以组成n*（n-1）/2对从而求解。

至于寻找相同斜率的直线，利用斜率排序。由于斜率为浮点数，不能直接比较大小，于是进行转化，将y2/x2=y1/x1转化为y2*x1=y1*x2。最后，只要单独处理一下斜率不存在的情况就好了

 

 

/*
Author：Owen_Q
*/

#include <bits/stdc++.h>

using namespace std;

typedef struct LINE
{
    int x;
    int y;
    double k;
    bool ek;
}Line;

vector <Line> l;
vector <int> cnt;


bool cmp(const Line x,const Line y)
{
    if(x.ek&&y.ek)
    {
        return x.k<y.k;
    }
    else if(!x.ek)
    {
        return true;
    }
    else
    {
        return false;
    }
}

int main()
{
    int n;
    while(scanf("%d",&n)!=EOF)
    {
        l.clear();
        cnt.clear();
        for(int i=0;i<n;i++)
        {
            int x1,x2,y1,y2;
            scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
            Line temp;
            temp.x=x1-x2;
            temp.y=y1-y2;
            if(temp.x==0)
            {
                temp.ek=false;
                temp.k=0;
            }
            else
            {
                temp.k=(double)(temp.y)/(double)(temp.x);
                temp.ek=true;
            }
            l.push_back(temp);
     //       cout<<(temp).x<<" "<<(temp).y<<" "<<(temp).k<<endl;
        }
        sort(l.begin(),l.end(),cmp);
        vector <Line> :: iterator it;
        int tempint = 1;
        for(it=l.begin()+1;it!=l.end();it++)
        {
            if(((!(*it).ek)&&(!(*(it-1)).ek))||(((*it).ek)&&((*(it-1)).ek)&&(((*it).x)*((*(it-1)).y)==((*it).y)*((*(it-1)).x))))
            {
                tempint++;
                if(it==(l.end()-1))
                {
                    cnt.push_back(tempint);
                    tempint = 0;
                    break;
                }
            }
            else
            {
                cnt.push_back(tempint);
                tempint = 1;
                if(it==(l.end()-1))
                {
                    cnt.push_back(tempint);
                    tempint = 0;
                    break;
                }
            }
        }
        vector <int> :: iterator itcnt;
        long long sum=0;
        if(!cnt.empty())
        {
            for(itcnt=cnt.begin();itcnt!=cnt.end();itcnt++)
            {
                sum += (long long)(*itcnt)*(long long)((*itcnt)-1)/2;
            }
        }
        printf("%lld\n",sum);
    }
    return 0;
}