C++博客-A code a day, keeps the girls away!-随笔分类-计算几何

hdu 3918 Beiju

AmazingCaddy — Sat, 13 Aug 2011 08:25:00 GMT

摘要: http://acm.hdu.edu.cn/showproblem.php?pid=3918一个如上图所示的杯子，一开始为空，且杯子的重量不计，沿着杯壁往里面慢慢地倒水，直到杯子倒了为止，最高能往里面倒多少水，求最后水的高度。做法：将杯身分割成梯形，每个梯形中，重心是在x轴的分量，是往一个方向偏移，也就是有单调性。求出从下往上枚举每个梯形，求出第一个使得杯子倒掉的梯形，然后在这个梯形内部二分，求出... 阅读全文

AmazingCaddy 2011-08-13 16:25 发表评论

sgu 435 UFO Circles

AmazingCaddy — Mon, 04 Oct 2010 18:08:00 GMT

摘要: http://acm.sgu.ru/problem.php?contest=0&problem=435题目大意，每个UFO会使长草的地面变成荒地，使荒地长出草来，作用范围是一个圆，求最后荒地和草地面积各为多少。圆的离散化，比赛的时候没有想仔细，没有做出来，赛后经haozi一点拨，发现可以做，就拿以前写的一个圆离散化的代码改了改，结果精度不够。重新写了一个，结果打错了一个变量，一直没有发现，... 阅读全文

AmazingCaddy 2010-10-05 02:08 发表评论

eoj 2830 Hamster 2

AmazingCaddy — Tue, 28 Sep 2010 16:27:00 GMT

摘要: http://www.acm.cs.ecnu.edu.cn/problem.php?problemid=2830题目大意：空中有很多竖直的门，一个点以V0的初速度抛出，忽略空气阻力，重力加速度为10m/s2,求此点最多能穿过几个门。解法：对于每个门，根据斜抛运动，解出初速度与x轴的夹角范围，然后求出重叠次数最大的区域，输出该次数。 eoj 2830Code highlighting produc... 阅读全文

AmazingCaddy 2010-09-29 00:27 发表评论

hdu 2966 In case of failure

AmazingCaddy — Tue, 07 Sep 2010 03:26:00 GMT

摘要: http://acm.hdu.edu.cn/showproblem.php?pid=2966题目的意思是：平面上有n个点（n<100000），求每个点的最近点到该点的平方距离。KD_Tree可以解决此题。详细资料可以参看此链接 http://en.wikipedia.org/wiki/Kd-tree，上面给出了算法。PS: 这道题时限开了恐怖的30秒。 hdu_2966Code high... 阅读全文

AmazingCaddy 2010-09-07 11:26 发表评论

hdu 3492 Segment

AmazingCaddy — Fri, 06 Aug 2010 10:26:00 GMT

摘要: 很郁闷的一题，第一次碰到使用不同的编译器，一个超时（C++），一个AC的情况（G++）。。。首先题中的要求等价于：存在一条直线l和所有的线段都相交。证明：若存在一条直线l和所有线段相交，作一条直线m和l垂直，则m就是题中要求的直线所有线段投影的一个公共点即为垂足。（l和每条线段的交点沿l投影到m上的垂足处）反过来，若存在m，所有线段在m上的投影有公共点，则过这点垂直于m作直线l， ... 阅读全文

AmazingCaddy 2010-08-06 18:26 发表评论

pku 1418 Viva Confetti

AmazingCaddy — Sun, 25 Jul 2010 11:14:00 GMT

摘要: http://124.205.79.250/JudgeOnline/problem?id=1418 最近跟着haozi大牛学习计算几何，第一道圆的离散化的题目。题目的大致意思是按照时间顺序把许多圆放在平面上，后放的圆可能将先放的圆... 阅读全文

AmazingCaddy 2010-07-25 19:14 发表评论

hdu 3471 England vs Germany

AmazingCaddy — Sun, 25 Jul 2010 04:21:00 GMT

摘要: http://acm.hdu.edu.cn/showproblem.php?pid=3471此题主要有两个问题：直线与面的交点，点是否在多边形内对于速度V的方向有三种情况，可用V与ABCD面法向量的点积判断：1 V指向ABCD面外侧，不可能进球；2 V与ABCD面平行，不可能进球；3 V指向ABCD面内侧，可能进球，分三种情况: 3.1 P在ABCD面内侧，不可能进球；... 阅读全文

AmazingCaddy 2010-07-25 12:21 发表评论

hdu 3437 Garden

AmazingCaddy — Sun, 25 Jul 2010 04:14:00 GMT

摘要: http://acm.hdu.edu.cn/showproblem.php?pid=3437AC牛的题目，精度比较恶心，最后一组数据一个误差是0.00006 我精度用了1e-8，而标程是1e-4，卡死在最后一组数据上了。。。这题看起来是三维的几何，其实完全可以压缩到二维来做。 1// Author : AmazingCaddy &n... 阅读全文

AmazingCaddy 2010-07-25 12:14 发表评论

pku 3608 Bridge Across Islands

AmazingCaddy — Mon, 12 Jul 2010 14:42:00 GMT

http://124.205.79.250/JudgeOnline/problem?id=3608

旋转卡壳的经典题目，看着网上某篇论文写的，写了一个下午还是WA，最后参考了了haozi的代码，还是WA，原来haozi的代码中将平行跟其中一种情况合并在一起写了，精度调得很高，但是我的精度太低了，调整精度之后过了。但是我觉得这题的精度是不用那么高的，原因应该是出在了，两种情况的合并上面，后来我改成三种情况，精度就1e-8过了，看来是那哥问题。代码如下：

pku 3608

AmazingCaddy 2010-07-12 22:42 发表评论

2008-2009 world final F题 Deer-Proof Fence

AmazingCaddy — Thu, 03 Jun 2010 12:47:00 GMT

http://acmicpc-live-archive.uva.es/nuevoportal/

先预处理所有凸包，然后dp。

1#include<iostream>
2#include<cmath>
3#include<algorithm>
4#include<complex>
5#include<vector>
6#include<stdio.h>
7using namespace std;
8typedef complex<double> point;
9const int maxn = 530;
10const double eps = 1e-8;
11const double inf = 1e10;
12const double pi = acos( -1.0 );
13double dp[maxn];
14int c[10];
15bool cmp( point p1, point p2 ){ return real(p1) < real(p2) || real(p1) == real(p2) && imag(p1) < imag(p2); }
16double operator^( point p1, point p2 ){ return imag( conj( p1 ) * p2 ); }
17int dblcmp( double x ){ return x < -eps ? -1 : x > eps ; }
18double minn( double x, double y ){ return x < y ? x : y; }
19void convex_hull( vector<point> &p, vector<point> &poly )
20{
21    int sz = 0, i, k;
22    poly.clear( );
23    sort( p.begin( ), p.end( ), cmp );
24    for( i = 0; i < p.size( ); ++i )
25    {
26        while( sz >= 2 && dblcmp( poly[sz-2] - p[i] ^ poly[sz-1] - p[i] ) <= 0 )
27            poly.pop_back( ), sz--;
28        poly.push_back( p[i] );
29        sz++;
30    }
31    k = sz;
32    for( i = p.size( ) - 2; i >= 0; i-- )
33    {
34        while( sz > k && dblcmp( poly[sz-2] - p[i] ^ poly[sz-1] - p[i] ) <= 0 )
35            poly.pop_back( ), sz--;
36        poly.push_back( p[i] );
37        sz++;
38    }
39    poly.pop_back( );
40}
41
42int main( )
43{
44    double a[maxn], x, y;
45    int i, j, q, n, m, len, k = 1, sz;
46    vector<point> p, poly, poly1;
47    while( 1 )
48    {
49        p.clear( );
50        scanf("%d%d",&n,&m);
51        if( n == 0 && m == 0 )break;
52        for( i = 0; i < n; i++ )
53        {
54            scanf("%lf%lf",&x,&y);
55            p.push_back( point( x, y ) );
56        }
57        len = 1 << n;
58        for( i = 1; i < len; i++ )
59        {
60            poly.clear( );
61            poly1.clear( );
62            j = i; q = 0;
63            while( j )
64            {
65                if( j & 1 ) poly.push_back( p[q] );
66                q++;
67                j >>= 1;
68            }
69            convex_hull( poly, poly1 );
70            sz = poly1.size( );
71            for( a[i] = 2 * pi * m, j = 0; j < sz; j++ )
72                 a[i] += abs( poly1[j] - poly1[(j+1)%sz] );
73        }
74        for( dp[0] = 0.0, i = 1; i < len; i++ )
75            dp[i] = inf;
76        for( i = 0; i < len; i++ )
77            for( j = 1; j < len; j++ )
78                if( ( j & i ) == 0 )
79                    dp[i|j] = minn( dp[i] + a[j], dp[i|j] );
80        printf("Case %d: length = %.2lf\n",k++,dp[len-1]);
81    }
82    return 0;
83}
84

AmazingCaddy 2010-06-03 20:47 发表评论

fzu 1918 John’s Direction

AmazingCaddy — Mon, 31 May 2010 17:27:00 GMT

摘要: http://acm.fzu.edu.cn/problem.php?pid=19182010福州邀请赛 J 题 #include#include#include#include#includeusing namespace s... 阅读全文

AmazingCaddy 2010-06-01 01:27 发表评论

eoj 1990 Find the Point

AmazingCaddy — Fri, 14 May 2010 15:11:00 GMT

摘要: http://www.acm.cs.ecnu.edu.cn/problem.php?problemid=1990 题目的意思是求出一些点，使得这些点到给定的一些直线的距离相等，只有一个点的时候输出该点，多个点那么输出"Many"，没有符合的点那么输出"None"，大致的做法是：先用角平分线求出一定量的点（最多4个），然后枚举检查这些点，精度有点搞！！！ #include&... 阅读全文

AmazingCaddy 2010-05-14 23:11 发表评论

hdu 3408 Elune’s Arrow

AmazingCaddy — Wed, 12 May 2010 16:06:00 GMT

摘要: http://acm.hdu.edu.cn/showproblem.php?pid=3408浙江理工邀请赛D题，比赛的时候没有搞出来，赛后终于一次AC了，题目其实还是比较简单的，模拟一下就可以了。 #include#include#include#includeusing&n... 阅读全文

AmazingCaddy 2010-05-13 00:06 发表评论

poj 3525 Most Distant Point from the Sea

AmazingCaddy — Sun, 02 May 2010 16:41:00 GMT

摘要: http://acm.pku.edu.cn/JudgeOnline/problem?id=3525半平面交 #include#include#include#include#include#include#include&... 阅读全文

AmazingCaddy 2010-05-03 00:41 发表评论

hdu_3372 Build Fences

AmazingCaddy — Wed, 07 Apr 2010 17:04:00 GMT

摘要: http://acm.hdu.edu.cn/showproblem.php?pid=3372太晚了，先贴代码。。。 1#include 2#include 3#include 4#include阅读全文

AmazingCaddy 2010-04-08 01:04 发表评论