BUPTOJ85 Three Points On A Line

考研党200年没码过代码了，没想到一上手这么生，还是和ACM大牛们没法比啊╮(╯▽╰)╭但好歹也是技术流，码代码还算有点自信。。。为了准备复试，抱着把北邮OJ上所有题刷一遍的心情就去了。。

PS：居然在论坛听说，有考研400分的大神，不、会、写、代、码！！

BUPT83 A+B 无压力秒杀！！！First Blood！！

BUPT84 Single Number 本来以为能轻松Double Kill的。。没想到跪在超时上了，到现在还没想明白

BUPT85 Three Points On A Line 没想到连跪两题，有辱ACM业余队员称号（低调低调） =  =

BUPT86...

不行！！必须把BUPT85做出来！！决不能再这样水下去了！！！我是要成为大牛的男人！！于是——

**********************************

上面不是分割线，是我说的话。解决了看错题、拼写错误等低端失误之后，还是答案错误。最后，我只好百度了。。。看到了一大牛写的解答。。暴力归暴力，我终于知道自己哪错了。

http://blog.csdn.net/birdstorm_s/article/details/21017411

之所以称之为大牛，是因为大牛总是在字里行间体现水平。他说有O(n^2)的算法，这个可是我百度好久之后才找到的，3SUM Hard问题。。

http://cs.smith.edu/~orourke/TOPP/P11.html#Problem.11

我确定我这不是算法问题，我怀疑是float的精度问题。。。

我判断三点共线用的斜率，但斜率就要除法，两个三位小数精度的数除完之后的精度是什么来着？好像物理实验课讲过。。可惜已经忘了。。我试了10^-4、10^-5，都报错。

大牛的判断方法把斜率相等的公式推导了一下，变成了一个乘法和减法的算式。这个好，不仅效率比除法高，而且精度易计算，10^-6。

判断方法改后，直接通过了。。仅以此文记这艰难的AC，展望更加艰难的未来╮(╯▽╰)╭

/*USER_ID: test#aa3615058PROBLEM: 85SUBMISSION_TIME: 2014-03-17 23:06:13*/#include <iostream>#include <math.h>using namespace std; typedef struct Point{    float x;    float y;}Point; bool isOnALine(Point a, Point b, Point c);float slope(Point a, Point b); int main() {    int caseCount;    bool onaline;    cin>>caseCount;    for(int i = 0; i < caseCount; i++) {        int pointCount;        cin>>pointCount;         Point ps[pointCount];        for(int j = 0; j < pointCount; j++) {            cin >> ps[j].x >> ps[j].y;        }        onaline = false;        if(pointCount >= 3) {            Point a,b,c;            for(int k = 0; k < pointCount && !onaline; k++) {                a = ps[k];                for(int l = k + 1; l < pointCount && !onaline; l++) {                    b = ps[l];                    for(int g = l + 1; g < pointCount && !onaline; g++) {                        c = ps[g];                        if(isOnALine(a,b,c)) {                            onaline = true;                            break;                        }                    }                }            }        }        if(onaline) {            cout << "Yes" << endl;        } else {            cout << "No" << endl;        }    }     return 0;} bool isOnALine(Point a, Point b, Point c) {    return (fabs(a.x * c.y + b.x * a.y + b.y * c.x - b.x * c.y - a.y * c.x - b.y * a.x) <= 1e-6);} float slope(Point a, Point b) {    return (a.y - b.y) / (a.x - b.x);}