Codeforces Beta Round #13 D —— Triangles 计算几何

原题：http://codeforces.com/contest/13/problem/D

题意：有n个红点，m个蓝点，在所有红点中选择三个点组成三角形，使得任意蓝点都不在该三角形内，问这样的三角形有多少个；

思路：首先，如何判断一个点是否在该三角形内 —— 我们假定三角形ABC以及所要判断的点P，如果P在ABC内，那么对边AB来说，点P和点C在边AB的同侧；对边BC来说，点P和点A在边BC的同侧，同样，边AC也是如此。我们可以用叉积来实现这个想法，叉积>0为同侧，叉积<0为异侧。

我们设一个原点O（其坐标应在输入范围以外），每次在红点中找两个点，求这两个点与O组成的三角形中有多少个蓝点，用dp[i][j]保存下来；

我们假设在红点中选三个点组成三角形，标记为ABC，那么△ABC内蓝点个数 = △ OAB + △OBC - △ OAC；

#include<bits/stdc++.h>
#define ll __int64
using namespace std;
const int maxn = 510;
int n, m;
int dp[maxn][maxn];
struct node{
    ll x, y;
    node() {}
    node(ll xi, ll yi):x(xi), y(yi){}
    node operator - (const node &N) const{
        return node(x-N.x, y-N.y);
    }
    ll operator * (const node &N) const{
        return x*N.y-y*N.x;
    }
}r[maxn], b[maxn];

int main(){
    cin>>n>>m;
    for(int i = 0;i<n;i++)  cin>>r[i].x>>r[i].y;
    for(int i = 0;i<m;i++)  cin>>b[i].x>>b[i].y;
    node o = node(-1000000001, -1000000001);
    for(int i = 0;i<n;i++){
        for(int j = 0;j<n;j++){
            if((r[i]-o)*(r[j]-o) <= 0)  continue;
            for(int k = 0;k<m;k++){
                if((r[i]-o)*(b[k]-o) > 0 && (r[j]-r[i])*(b[k]-r[i]) > 0 && (o-r[j])*(b[k]-r[j]) > 0)
                    dp[i][j]++;
            }
            dp[j][i] = -dp[i][j];
        }
    }
    int ans = 0;
    for(int i = 0;i<n;i++){
        for(int j = i+1;j<n;j++){
            for(int k = j+1;k<n;k++)
                ans += (dp[i][j]+dp[j][k]+dp[k][i] == 0);
        }
    }
    cout<<ans<<endl;
    return 0;
}