昨天晚上太卡了,竟然没有全部发完,今天早上早早起来发一下,1113是一个计算几何的题目,求得是凸包问题,就是在二维平面上给定一系列的点,求一个最小的凸包将它们全部覆盖,所谓凸包就是一个多边形,在这个多边形内的任意两点间的连线都在多边形内。地址: http://acm.pku.edu.cn/JudgeOnline/problem?id=1113
#include < stdio.h >
#include < stdlib.h >
#include < math.h >
typedef struct { double x; double y; int addr; int flag; }type; type p[1005];
int stack[1005]; int top;
void inits( ) {
top = -1; }
void in( int a ) {
stack[++top] = a; }
void out() {
top --; }
void f_min( int n ) {
int i; type min;
min.x = p[0].x; min.y = p[0].y; min.addr = p[0].addr; min.flag = p[0].flag;
for ( i=1; i<n; i++ ) { if ( min.y > p[i].y ) { min.x = p[i].x; min.y = p[i].y; min.addr = p[i].addr; min.flag = p[i].flag; } else { if ( min.y == p[i].y ) { if ( min.x > p[i].x ) { min.x = p[i].x; min.y = p[i].y; min.addr = p[i].addr; min.flag = p[i].flag; } } } } p[min.addr].x = p[0].x; p[min.addr].y = p[0].y; p[min.addr].addr = p[0].addr; p[min.addr].flag = p[0].flag; p[0].x = min.x; p[0].y = min.y; p[0].addr = min.addr; p[0].flag = min.flag; }
double mul( type *a, type *b, type *c, type *d ) {
double x1 = a->x - b->x; double y1 = a->y - b->y; double x2 = c->x - d->x; double y2 = c->y - d->y;
return x1 * y2 - x2 * y1; }
double dis( type *a, type *b ) {
double x = a->x - b->x; double y = a->y - b->y; return x*x + y*y; }
int cmp( const void *a, const void *b ) {
type *c = ( type * )a; type *d = ( type * )b;
int ans; if ( mul( c, &p[0], d, &p[0] ) > 0 ) { ans = -1; } else if ( mul( c, &p[0], d, &p[0] ) < 0 ) { ans = 1; } else { if ( dis( c, &p[0] ) > dis( d, &p[0] ) ) { ans = 1; } else { if ( dis( c, &p[0] ) < dis( d, &p[0] ) ) { ans = -1; } } }
return ans; }
int aobao( int n ) {
int i; int ans = 1;
f_min( n );
qsort( p+1, n-1, sizeof( type ), cmp );
inits();
in( 0 ); p[0].flag = 1;
for ( i=1; i<n-1; i++ ) { if ( mul( &p[i], &p[0], &p[i+1], &p[0] ) != 0 ) { break; } } in( i ); p[i].flag = 1;
if ( i == n-1 ) { ans = 0; return ans; }
for ( ++i; i<n-1; i++ ) { if ( mul( &p[i], &p[0], &p[i+1], &p[0] ) != 0 ) { break; } } in( i ); p[i].flag = 1; for ( ++i; i<n; i++ ) { if ( i < n-1 ) { if ( mul( &p[i], &p[0], &p[i+1], &p[0] ) != 0 ) { while ( mul( &p[i], &p[stack[top-1]], &p[stack[top]], &p[stack[top-1]] ) >= 0 ) { p[stack[top]].flag = 0; out(); }
in( i ); p[i].flag = 1; } } else { while ( mul( &p[i], &p[stack[top-1]], &p[stack[top]], &p[stack[top-1]] ) >= 0 ) { p[stack[top]].flag = 0; out(); }
in( i ); p[i].flag = 1; } }
return ans; }
int main() {
int n; double l; double ans; int temp; int i;
while ( scanf( "%d%lf", &n, &l ) != EOF ) { for ( i=0; i<n; i++ ) { scanf( "%lf%lf", &p[i].x, &p[i].y ); p[i].addr = i; p[i].flag = 0; }
temp = aobao( n ); ans = 0; for ( i=1; i<=top; i++ ) { ans += sqrt( dis( &p[stack[i]], &p[stack[i-1]] ) ); } if ( temp ) { ans += sqrt( dis( &p[stack[top]], &p[stack[0]] ) ); } ans += 2 * l * acos(-1); printf( "%d\n", ( int )( ans+0.5 ) ); }
return 0; }
|
|
| 日 | 一 | 二 | 三 | 四 | 五 | 六 |
---|
25 | 26 | 27 | 28 | 29 | 30 | 31 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 | 17 | 18 | 19 | 20 | 21 | 22 | 23 | 24 | 25 | 26 | 27 | 28 | 29 | 30 | 1 | 2 | 3 | 4 | 5 |
|
公告
决定从线程开始!!
常用链接
留言簿(6)
随笔档案
搜索
最新评论
阅读排行榜
评论排行榜
|
|