Hoj 11619 Balance 切割多边形

void sovle()

{

this->turn_back();

polygon a, b;

int i, j, k;

point tmp;

p[n] = p[0];

point p1( 0 , 0 ), p2( 0 , 0 );

for ( i = 0 ; i < n ; i ++ )

{

if ( p[i].y >= 0 )

a.add( p[i] );

if ( p[i].y <= 0 )

b.add( p[i] );

if( sig( p[i].y * p[i + 1 ].y ) < 0 )

{

tmp.x = p[i].x + ( p[i].y /( p[i].y - p[i+1].y ))*( p[i+1].x - p[i].x );

// 如此可包含 + - 号的问题

tmp.y = 0;

a.add( tmp );

b.add( tmp );

}

p1 = a.poly_centroid();

p2 = b.poly_centroid();

if ( sig( p1.x - p2.x ) == 0 )

printf("Balanced.\n");

else if ( sig ( p1.x - p2.x ) < 0 )

printf("CE is aft of CLR by %.2lf units.\n",fabs(p1.x-p2.x));

else

printf("CE is forward of CLR by %.2lf units.\n",fabs(p1.x-p2.x));

}

int main()

{

polygon poly;

while( poly.read() != EOF )

{

poly.sovle();

}

return 0;

}

posted on 2009-08-28 17:59 Huicpc217 阅读(103) 评论(0) 编辑收藏引用所属分类: 计算几何

只有注册用户登录后才能发表评论。
【推荐】100%开源！大型工业跨平台软件C++源码提供，建模，组态！

相关文章: 三维——面旋转 Hoj 11619 Balance 切割多边形计算几何模板

网站导航: 博客园 IT新闻 BlogJava 知识库博问管理

常用链接

留言簿(2)

随笔档案

文章分类

文章档案

Hnu Acmer

搜索

最新评论

阅读排行榜

评论排行榜