hdu4380Farmer Greedy(多校九）

http://acm.hdu.edu.cn/showproblem.php?pid=4380

求三角形内的点的个数为奇数的三角形数 官方解题报告提供的是O（n2+m)的算法 先算两点与原地围成的三角形 最后总的有向三角形面积等于三者之和

实在不知道有向面积怎么转换成无向的 交了十几次依旧WA 在网上看到另一种做法 感觉很好

利用叉乘算出每个线段下面的点 三角形是由三条线段所围成 所以用一条减去两条（钝角）或者两条减一条（锐角）就能得到三角形区域 所以它包围的点数也就出来 了

#include <iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<string.h>
using namespace std;
#define N 1000000
struct node
{
    double x,y;
};
node q[101];
node p[1011];
int s[101][101];
int ft(double  x1,double  y1,double x2,double y2, double x3,double y3)
{
    if(x3<x1||x3>=x2)//在线段端点外的点不要
    return 0;
   if(((x1-x3)*(y2-y3)-(x2-x3)*(y1-y3))<0)//叉乘判断点在线段下方
   return 1;
   else
   return 0;
}
bool cmp(node a,node b)
{
    if(a.x==b.x)//按x排序
    return a.y<b.y;
    return a.x<b.x;
}
int main()
{
    int i,j,k,n,m,num,g,mm = 0;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        int tnum = 0;
        mm++;
        memset(s,0,sizeof(s));
        for(i = 1; i <= n ;i++)
        scanf("%lf%lf",&q[i].x,&q[i].y);
        for(j = 1; j <= m ; j++)
        scanf("%lf%lf",&p[j].x,&p[j].y);
        sort(q+1,q+n+1,cmp);
        for(i = 1; i <= n-1 ; i++)
        {
            for(j = i+1 ; j <= n ; j++)
            {
                for(g = 1 ;g <= m ; g++)
                {
                    if(ft(q[i].x,q[i].y,q[j].x,q[j].y,p[g].x,p[g].y))
                    s[i][j]++,s[j][i]++;
                }
            }
        }
        for(i = 1; i <= n-2 ; i++)//线段包围的点 互相减掉 减成一个三角形区域包围的点
        for(j = i+1; j <= n-1; j++)
        for(k = j+1; k <= n ; k++)
        if(abs(s[i][j]-s[j][k]-s[i][k])%2!=0)//锐角三角形 是两个减一个 钝角是一个减两个 直接用一个减会有负值 加上绝对值
        tnum++;
        printf("Case %d: ",mm);
        printf("%d\n",tnum);
    }
    return 0;
}