POJ 2280&&hdu 1661

题意：给定平面上的N个点，属性分别标记为0和1，然后找一条直线，直线上的点全部溶解，一侧的1溶解，另一侧的0溶解。求出最多能溶解的点的个数。

思路：暴力枚举每个点，扫描线旋转。先做优化，如果一侧溶解0，则把属性为1的做关于当前枚举直线对称的点，这样统计一侧的点加上线上的点就是答案。O(n2).

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<set>
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define clc(a,b) memset(a,b,sizeof(a))
#define eps 1e-8
typedef long long LL;
const int mod=47000;
const int inf=0x3f3f3f3f;
using namespace std;
#define N 100010

struct node
{
    int x,y,r;
    double ang;
}v[1010],rem[1010];

int cmp(const struct node a,const struct node b)
{
    if(a.ang<b.ang)
        return 1;
    else
        return 0;
}

int sig (double a)
{
    if(fabs(a)<eps)
        return 0;
    else if(a>0)
        return 1;
    else
        return -1;
}

int cross(struct node a,struct node b,struct node c)
{
    return (b.x-a.x)*(c.y-a.y)-(c.x-a.x)*(b.y-a.y);
}

int main()
{
    int n;
    while(~scanf("%d",&n)&&n)
    {
        int ans=0;
        for(int i=0;i<n;i++)
        {
            scanf("%d%d%d",&v[i].x,&v[i].y,&v[i].r);
            rem[i]=v[i];
        }
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<n;j++)
            {
                v[j]=rem[j];
                if(v[j].r==1)
                {
                    v[j].x=rem[i].x*2-v[j].x;
                    v[j].y=rem[i].y*2-v[j].y;
                }
                v[j].ang=atan2(v[j].y-rem[i].y,v[j].x-rem[i].x);
            }
            swap(v[i],v[0]);
            sort(v+1,v+n,cmp);
            for(int j=1,t=2;j<n&&sig(v[j].ang)<=0;j++)
            {
                int on=2;
                for( ;t<n&&cross(v[0],v[j],v[t])>=0;t++)
                {
                    if(cross(v[0],v[j],v[t])==0)
                    on++;
                }
                ans=max(ans,max(t-j+1,n-(t-j+1)+on));
            }
        }
        cout<<ans<<endl;
    }
    return 0;
}