hdu 3685 Rotational Painting

题目：有一副玻璃画，要立在桌子上，可以在垂直于桌面的平面旋转，问有几种稳定的摆放状态。

分析：计算几何、凸包、重心。摆放条件：首先、能够稳定摆放一定是用凸包的边与桌面接触；其次、如果稳定，那么重心必须落在对应的凸包的边内（边上不行）。

            1.重心计算：将多边形分割成三角形（分割形式任意），然后每个三角形用自己的重心代替，看成质点，质量为三角形面积（利用叉乘计算），最后求出所有质点的平均值即可；此算法和面积计算一样与凸凹无关。

            2.稳定状态：对于给已求得的重心枚举凸包的边即可；对于稳定的判断，可以得到一些判断条件：

               (1).两底角为锐角，会出现精度问题。

               (2).利用叉乘计算，稍稍有点麻烦。

               (3).利用解析几何的可行域求解，即对于给定直线方程带入两点结果符号相反。

               这里利用方案3，计算方便，也不会出现精度问题。已知重心c，凸包底边ab，则(b.x-a.x)*(x-c.x)+(b.y-a.y)*(y-c.y)=0即为目标方程。

注意：精度问题、过多的浮点运算会出现精度问题。

#include <algorithm> #include <iostream>#include <cstdlib>#include <cstdio>#include <cmath>using namespace std;typedef struct pnode{double x,y,d;}point; point P[ 50005 ];point S[ 50005 ];point MP;double crossproduct( point a, point b, point c )//ab到ac {return (b.x-a.x)*(c.y-a.y)-(c.x-a.x)*(b.y-a.y);}double dist( point a, point b ){return sqrt((b.x-a.x)*(b.x-a.x)+(b.y-a.y)*(b.y-a.y));}point center( int N ){double sum = 0.0,s;point  p; p.x = p.y = 0.0;if ( N == 0 ) return p;if ( N == 1 ) return P[0];if ( N == 2 ) {p.x = (P[0].x+P[1].x)/2;p.y = (P[0].y+P[1].y)/2;return p;}for ( int i = 2 ; i < N ; ++ i ) {s = crossproduct( P[0], P[i-1], P[i] );p.x += (P[0].x+P[i-1].x+P[i].x)*s;p.y += (P[0].y+P[i-1].y+P[i].y)*s;sum += s;}p.x /= 3*sum; p.y /= 3*sum;return p;}bool cmp1( point a, point b ){if ( a.x == b.x ) return a.y < b.y;else return a.x < b.x;}bool cmp2( point a, point b ){return crossproduct( P[0], a, b ) > 0;}bool cmp3( point a, point b ){double cp = crossproduct( P[0], a, b );if ( cp == 0.0 ) {if ( !crossproduct( P[0], a, MP ) )return a.d < b.d;else return a.d > b.d;}else return cp > 0;}int Graham( int N ){point c = center( N );sort( P+0, P+N, cmp1 );sort( P+1, P+N, cmp2 );//处理共线，去掉共线的点 for ( int i = 1 ; i < N ; ++ i )P[i].d = dist( P[0], P[i] );MP = P[N-1];sort( P+1, P+N, cmp3 );int top = -1;if ( N > 0 ) S[++ top] = P[0];if ( N > 1 ) S[++ top] = P[1];if ( N > 2 ) {S[++ top] = P[2];for ( int i = 3 ; i < N ; ++ i ) {while ( crossproduct( S[top-1], S[top], P[i] ) < 0 ) -- top;S[++ top] = P[i];}}S[++ top] = P[0];int    count = 0;double A,B;for ( int i = 0 ; i < top ; ++ i ) {A = S[i+1].x-S[i].x;B = S[i+1].y-S[i].y;if ( (A*(S[i].x-c.x)+B*(S[i].y-c.y))*(A*(S[i+1].x-c.x)+B*(S[i+1].y-c.y)) < 0 )count ++;}return count;}int main(){int  N,T; while ( scanf("%d",&T) != EOF ) while ( T -- ) {scanf("%d",&N);for ( int i = 0 ; i < N ; ++ i ) scanf("%lf%lf",&P[i].x,&P[i].y);printf("%d\n",Graham( N ));}return 0;}