bzoj1132 [POI2008]Tro（叉积）

bzoj1132 [POI2008]Tro

原题地址：http://www.lydsy.com/JudgeOnline/problem.php?id=1132

题意：
平面上有N个点. 求出所有以这N个点为顶点的三角形的面积和 N<=3000

数据范围
N在[3,3000] ，N个点的坐标值在[0,10000]
保留一位小数,误差不超过0.1

题解：
好题。
我想了两种做法，
第一种：考虑一般求多边形面积的方法。

S=S蓝1+S蓝2-S红
其实就是计算每条边的贡献。
对于x以右的点极角排序，
ab 这样的边贡献+1
（保存一个两两边贡献的后缀和）
x与a这样的点连边的贡献就（前面的点的个数-后面的点的个数）

正解：
考虑叉积到底是个什么。
x1y2-x2y1
如果a对于x的坐标（x1,y1），b对于x的坐标（x2,y2）
提个公因数，
x1贡献 ∑y2
y1贡献-∑x2

因为要用叉积算，要保证>0，于是让它有序就行。
先按x排序，枚举每个点，他之右的极角排序，然后存个前缀和算。

---

一个困扰我很久的点：
第一遍排序:bool cmp1(const Poi &A,const Poi &B){return A.x<B.x;}
结果：WA
第一遍排序:bool cmp1(const Poi &A,const Poi &B){return A.x<B.x||(A.x==B.x&&A.y<B.y);}
结果：AC

一组数据：
（10，1）（10，2）（10，3） ，（15，2）
正确答案是10.0
但对于第一种排序，显然它对于前三个点的顺序不作要求，
当顺序为：（10，2）（10，3） （10，1），（15，2）时，
答案为 5.0，错了！！！
而对称的（10，2）（10，1） （10，3），（15，2）
答案则是正确的 10.0
为什么？感谢Doggu非常耐心的答疑。
打印出过程发现，同样第一个数取（10，2），那么剩下三个点对于他的坐标分别为：
（5，0）（0，1）（0，-1）
第一种的极角排序：（0，-1）（0，1）（5，0）
正常的极角排序：（0，-1）（5，0）（0，1）
极角排序不对！！！
考虑到我们为了追求常数小，比较极角用的不是atan2，而是叉积：

bool cmp2(const Poi &A,const Poi &B){return A.x*B.y-B.x*A.y>0;}

那么对于（0，-1）和（0，1）他们的叉积为0，其实是互相小于的，一旦比较顺序不对，排序就不对。

因此必须按x为第一关键字，y为第二关键字排，才不会出现这种上下都有的情况。

下次用叉积也要小心了。

---

代码：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#define LL long long
using namespace std;
const int N=3005;
struct Poi
{
    int x,y;
    Poi(){}
    Poi(int x,int y):x(x),y(y){}
}P[N],v[N];
Poi operator-(const Poi &A,const Poi &B){return Poi(A.x-B.x,A.y-B.y);}
bool cmp1(const Poi &A,const Poi &B){return A.x<B.x||(A.x==B.x&&A.y<B.y);}
bool cmp2(const Poi &A,const Poi &B){return A.x*B.y-B.x*A.y>0;}
int n;
LL ret=0;
void solve()
{
    sort(P+1,P+n+1,cmp1);
    for(int i=1;i<n-1;i++)
    {
        int cnt=0; LL sx=0,sy=0;
        for(int j=i+1;j<=n;j++)
        {
            v[++cnt]=P[j]-P[i];
        }
        sort(v+1,v+cnt+1,cmp2);
        for(int j=1;j<=cnt;j++)
        {       
            ret-=1LL*v[j].x*sy;
            ret+=1LL*v[j].y*sx;
            sx+=v[j].x; sy+=v[j].y;
        }
    }
    printf("%lld.%lld\n",ret>>1,5LL*(ret&1));
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++) scanf("%d%d",&P[i].x,&P[i].y);
    solve();
    return 0;
}