poj2280--Amphiphilic Carbon Molecules(扫描线+极角排序+转换坐标)
2015-08-22 15:50
429 查看
题目链接:点击打开链接
题目大意:给出n个点的坐标,每个点有一个值0或者1,现在有一个隔板(无限长)去分开着n个点,一侧统计0的个数,一侧统计1的个数,如果点在板上,那么都被统计到,问最多可以统计到多少个点。
对每个点都做一次极角排序,然后按照排序后的点进行遍历板的位置,统计出在板上的点的个数,两侧的0和1的个数。然后求最大值。
但是极角排序的基准点只能是最左下角的点,所以还要进行坐标转化,如果某一个点再当前基准点的下方,那么就可以转换坐标,改变那个点的坐标为与基准点对称的点,同时点的值改变。
题目大意:给出n个点的坐标,每个点有一个值0或者1,现在有一个隔板(无限长)去分开着n个点,一侧统计0的个数,一侧统计1的个数,如果点在板上,那么都被统计到,问最多可以统计到多少个点。
对每个点都做一次极角排序,然后按照排序后的点进行遍历板的位置,统计出在板上的点的个数,两侧的0和1的个数。然后求最大值。
但是极角排序的基准点只能是最左下角的点,所以还要进行坐标转化,如果某一个点再当前基准点的下方,那么就可以转换坐标,改变那个点的坐标为与基准点对称的点,同时点的值改变。
#include <cstdio> #include <cstring> #include <cmath> #include <algorithm> using namespace std ; #define eqs 1e-9 struct node{ double x , y ; int k ; }p[1100] , q , a[1100] ; node sub(node a,node b) { a.x -= b.x ; a.y -= b.y ; return a ; } double dis(node a,node b) { double l1 = a.x*a.x+a.y*a.y ; double l2 = b.x*b.x+b.y*b.y ; if( l2-l1 >= eqs ) return 1 ; else return 0 ; } int mul(node a,node b) { double s = a.x*b.y - a.y*b.x ; if( s >= eqs ) return 1 ; else if( s <= -eqs ) return -1 ; else return 0 ; } int cmp(node a,node b) { a = sub(a,q) ; b = sub(b,q) ; int k = mul(a,b) ; if( k == 1 ) return 1 ; else if( k == -1 ) return 0 ; else { return dis(a,b) ; } } int main() { int n , m , ans , sum0 , sum1 ; int s , i , j , k ; int cnt0 , cnt1 , num0 , num1 ; while( scanf("%d", &n) && n ) { for(i = 0 ; i < n ; i++) { scanf("%lf %lf %d", &p[i].x, &p[i].y, &p[i].k) ; } q.x = q.y = 0 ; ans = 0 ; for(s = 0 ; s < n ; s++) { sum0 = sum1 = 0 ; for(i = m = 0 ; i < n ; i++) { if( i == s ) continue ; a[m].x = p[i].x - p[s].x ; a[m].y = p[i].y - p[s].y ; a[m].k = p[i].k ; if( a[m].y < 0 || (a[m].y == 0 && a[m].x < 0) ) { a[m].x = -a[m].x ; a[m].y = -a[m].y ; a[m].k = 1 - a[m].k ; } if( a[m].k ) sum1++ ; else sum0++ ; m++ ; } sort(a,a+m,cmp) ; cnt0 = cnt1 = 0 ; num0 = num1 = 0 ; for(i = j = 0 ; i < m ; i++) { if( a[i].k ) num1++ ; else num0++ ; while( mul( sub(a[j],q) , sub(a[i],q) ) ) { if( a[j].k ) cnt1++ ; else cnt0++ ; j++ ; } ans = max(ans,max(sum0-num0+cnt1,sum1-num1+cnt0)+(num0-cnt0+num1-cnt1+1) ) ; } } printf("%d\n", ans) ; } return 0 ; }
相关文章推荐
- TCP协议中的三次握手及四次握手图解
- RPC提高
- POJ3278Catch That Cow简单一维广搜
- Hibernate学习笔记
- android 判断 activity 状态
- 模式设计的秘密之模板方法模式
- C++打印地址
- map的使用
- jQuery jsonp无法捕获404、500状态错误
- struts2学习笔记
- Linux Shell 小脚本经典收藏
- 51CTO下载-千万级大型系统集成项目
- Python 初学(4)——再谈字符串
- linux内核线程睡眠与唤醒
- RPC基础
- HDU1240Asteroids!三维广搜bfs
- contiki_protothread源码
- linux基本命令:
- 电面
- Codeforces 27 E Number With The Given Amount Of Divisors