C++博客-Uriel's Corner-文章分类-计算几何

[LeetCode]2251. Number of Flowers in Full Bloom (Hard) Python3-2023.10.11

Uriel — Mon, 06 Nov 2023 17:53:00 GMT

给出每个花的开花时间（起点终点），输出people数组每个人到达时有多少花在开。线段树、扫描线、二分

1 #2251
2 #Runtime: 1083 ms (Beats 15.71%)
3 #Memory: 44.7 MB (Beats 32.83%)
4
5 class Solution:
6     def fullBloomFlowers(self, flowers: List[List[int]], people: List[int]) -> List[int]:
7         st = sorted(flowers, key=lambda x:x[0])
8         ed = sorted(flowers, key=lambda x:x[1])
9         ans = []
10         for x in people:
11             i = bisect_left(ed, x, key=lambda x:x[1])
12             j = bisect_right(st, x, key=lambda x:x[0])
13             ans.append(j - i)
14         return ans

Uriel 2023-11-07 01:53 发表评论

[LeetCode]149. Max Points on a Line (Hard) Python3-2023.01.08

Uriel — Sun, 08 Jan 2023 10:10:00 GMT

给出一堆二维点，问最多多少个点共线
O(n²)枚举两个点，看一样斜率的最多多少个点，2014年曾经用C++写过👉http://www.cppblog.com/Uriel/articles/205287.html
今日在Discussion看到个不错的思路（👉https://leetcode.com/problems/max-points-on-a-line/solutions/3016632/python-3-11-lines-w-explanation-and-example-t-m-95-97/），不需要折腾double型求斜率，因为点的坐标都是int型，可以求两个点dx，dy，除以GCD之后用dict统计这样的约简后的数对有多少个，因为存的是除以GCD之后的数对，所以一开始要给所有点按x值从小到大排序，保证单调增

1 #149
2 #Runtime: 77 ms (Beats 92.29%)
3 #Memory: 13.8 MB (Beats 94.26%)
4
5 class Solution:
6     def maxPoints(self, points: List[List[int]]) -> int:
7         points.sort()
8         ans = 0
9         for i, (x1, y1) in enumerate(points):
10             k = defaultdict(int)
11             for x2, y2 in points[i + 1 :]:
12                 dx = x2 - x1
13                 dy = y2 - y1
14                 g = gcd(dx, dy)
15                 kk = (dx // g, dy // g)
16                 k[kk] += 1
17                 ans = max(ans, k[kk])
18         return ans + 1

Uriel 2023-01-08 18:10 发表评论

[LeetCode]587. Erect the Fence (Hard) Python-2022.11.19

Uriel — Sat, 19 Nov 2022 10:52:00 GMT

给出一堆树的坐标，用最小的多边形围住这些树，求问最后哪些树在多边形的顶点上，凸包模板题
借鉴了https://leetcode.com/problems/erect-the-fence/discuss/1442266/A-Detailed-Explanation-with-Diagrams-(Graham-Scan) 的做法，先对树坐标先x后y排序，在分上下半块凸包分别构建

1 #587
2 #Runtime: 562 ms
3 #Memory Usage: 14.2 MB
4
5 class Solution(object):
6     def corssProduct(self, p1, p2, p3):
7         return (p3[1]-p2[1]) * (p2[0]-p1[0]) - (p2[1]-p1[1]) * (p3[0]-p2[0])
8
9     def buildConvexHull(self, points):
10         lower = []
11         upper = []
12         for pt in points:
13             while len(lower) >= 2 and self.corssProduct(lower[-2], lower[-1], pt) < 0:
14                 lower.pop()
15             while len(upper) >= 2 and self.corssProduct(upper[-2], upper[-1], pt) > 0:
16                 upper.pop()
17             upper.append(pt)
18             lower.append(pt)
19         return upper, lower
20
21     def outerTrees(self, trees):
22
23         """
24         :type trees: List[List[int]]
25         :rtype: List[List[int]]
26         """
27         points = sorted(trees)
28         upper, lower = self.buildConvexHull(points)
29         return set(tuple(T) for T in (upper + lower))

Uriel 2022-11-19 18:52 发表评论

POJ 1851 Map---计算几何

Uriel — Sun, 26 Sep 2010 18:23:00 GMT

摘要: ECUST ACHunter威武，今天哈尔滨现场赛一度rank6(算上旅游队)，后2h可惜了，福州继续加油~~ 继续刷计算几何。。不会... 阅读全文

Uriel 2010-09-27 02:23 发表评论

POJ 1380 Equipment Box && POJ 2172 Bricks---计算几何

Uriel — Mon, 20 Sep 2010 14:06:00 GMT

        POJ 1380是去年十月做的，今天发现POJ 2172与POJ 1380一样的。。拿原来代码直接贴了。。

        1380是问一个矩形能否塞进另一个矩形，2172其实一回事，虽然变成三维，但是只要判三次就行了。。

如图，假设小矩形正好卡进去的情况，求出边长Len需要多少，Len>a则塞不进

代码中：∠1=∠sita      ∠1+∠2=∠p

以下是2172的代码，1380判定过程跟这个一样，而且只要判一次就行

//Problem: 2172  User: Uriel
//Memory: 208K  Time: 32MS
//Language: C++  Result: Accepted
//2010.09.20

#include<math.h>
#include<stdio.h>
#include<stdlib.h>

struct point {
    double x,y,h;
}a,b;

void Swap(double & a,double & b){
    double t;
    t=a;
    a=b;
    b=t;
}

bool check(double c,double d,double a,double b){
    if(a*b<c*d)return false;
    if(a<b)Swap(a,b);
    if(c<d)Swap(c,d);
    if(b<d)return false;
    if(a>=c && b>=d)return true;
    double dis=sqrt(c*c+d*d);
    double p=asin(b/dis);
    double q=asin(d/dis);
    double sita=p-q;
    double Len=c*cos(sita)+d*sin(sita);
    if(Len<=a)return true;
    return false;
}

int main(){
    scanf("%lf %lf %lf %lf %lf",&a.x,&a.y,&a.h,&b.x,&b.y);
    if(check(a.x,a.y,b.x,b.y)||check(a.x,a.h,b.x,b.y)||check(a.y,a.h,b.x,b.y))puts("YES");
    else
        puts("NO");
    return 0;
}

Uriel 2010-09-20 22:06 发表评论

POJ 2284 That Nice Euler Circuit---计算几何

Uriel — Sat, 18 Sep 2010 17:45:00 GMT

摘要: 这几天做计算几何一直不顺，各种WA。。这题做了一天多。。WA到哭了。。这题思路很简单，用平面图的欧拉定理，V+F-E=2；其中V是顶点数，F是分割出的面数，E是边数。 &... 阅读全文

Uriel 2010-09-19 01:45 发表评论

POJ 1931 Biometrics---计算几何

Uriel — Fri, 17 Sep 2010 09:59:00 GMT

        5门专业必修+4门大类选修真纠结。。微机、计算机图像处理还行，数据库也还行，计算机组成原理算得头大，最纠结的还是编译原理，K老师的课漏听30min肯定就听不懂了，作业又一堆。。每天白天都用来上课+写作业了。。就晚上选修课结束回1教切会题。。

        看起来很水的一道计算几何卡了两天才过。。悲剧啊。。

        一开始没注意相似的点是对应的，不会说第一个多边形的第1个点对应第二个多边形的第3个点之类的，想了几个WS算法，WA。。今天注意到了这点，结果改来改去还是不对。。最后看了某解题报告，把判转向那里改了一下，总算过了。。

        代码比较挫。。很多无用的东西懒得擦掉了。。

#include<math.h>
#include<stdio.h>
#include<stdlib.h>
#include<algorithm>
using namespace std;
#define eps 1e-6
const double INF=1e20;

struct point{
    double x,y;
    point operator-(point &b){
        point c;
        c.x = x - b.x;
        c.y = y - b.y;
        return c;
    }
}p1[11],p2[11];

int n;
int dir1[11],dir2[11];
double len1[11],len2[11];

double dis(point a,point b){
    return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
}

double Cos(point p0,point p1,point p2){
    return (dis(p0,p1)+dis(p1,p2)-dis(p0,p2))/(2*sqrt(dis(p0,p1))*sqrt(dis(p1,p2)));
}

int cross_product(point p1,point p2,point p0){
    return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}

double dot_product(point p1,point p2,point p0){
    return (p1.x-p0.x)*(p2.x-p0.x)+(p1.y-p0.y)*(p2.y-p0.y);
}

double cross(point a, point b){
    return a.x * b.y - a.y * b.x;
}

int get_dir(point p1,point p2,point p3){
    double t1=cross(p2-p1,p3-p2);
    if(fabs(t1)<eps)return 1;
    if(t1<0)return 2;
    else return 3;
}

int main(){
    int f1,f2,i,j,cnt;
    bool ok;
    double minn1,len,minn2,cons;
    while(scanf("%d",&n),n){
        for(i=0;i<n;i++)scanf("%lf %lf",&p1[i].x,&p1[i].y);
        for(i=0;i<n;i++)scanf("%lf %lf",&p2[i].x,&p2[i].y);
        cons=sqrt(dis(p1[0],p1[1]))/sqrt(dis(p2[0],p2[1]));
        for(i=0;i<n;i++){
            len1[i]=sqrt(dis(p1[i],p1[(i+1)%n]))/cons;
            len2[i]=sqrt(dis(p2[i],p2[(i+1)%n]));
        }
        ok=true;
        for(i=0;i<n;i++){
            if(fabs(len1[i]-len2[i])>eps){
                ok=false;
                break;
            }
        }
        if(ok){
            for(i=0;i<n;i++){
                if(fabs(fabs(Cos(p1[i],p1[(i+1)%n],p1[(i+2)%n]))-fabs(Cos(p2[i],p2[(i+1)%n],p2[(i+2)%n])))>eps){
                    ok=false;
                    break;
                }
            }
        }
        if(ok){
            for(i=0;i<n;i++){
                if(get_dir(p1[i],p1[(i+1)%n],p1[(i+2)%n])!=get_dir(p2[i],p2[(i+1)%n],p2[(i+2)%n])){
                    ok=false;
                    break;
                }
            }
        }
        if(ok)puts("similar");
        else
            puts("dissimilar");
    }
    return 0;
}

Uriel 2010-09-17 17:59 发表评论

POJ 3129 How I Wonder What You Are!----计算几何

Uriel — Tue, 07 Sep 2010 15:45:00 GMT

       找了一晚上计算几何题，看一道一道不会的。。终于找到一道大水。。。

       题意是有天空有N颗星星，给出它们的三维坐标。。在(0,0,0)处放着M个望远镜，给出它们的仰角(tx,ty,tz)以及望远镜圆锥形的顶角phi，求总共能看见多少星星。。

      方法就是暴力N*M，一个个的判，貌似是第一还是第二次用浙大模板的三维几何~~

//Problem: 3129  User: Uriel
//Memory: 228K  Time: 47MS
//Language: C++  Result: Accepted

#include<math.h>
#include<stdio.h>
#include<stdlib.h>

struct point{
    bool flag;
    double x,y,z;
}p[550];

struct tele{
    double x,y,z,tx,ty,tz,phi;
}tel[55];

int n,m;

double vlen(point p){
    return sqrt(p.x*p.x+p.y*p.y+p.z*p.z);
}

double dmult(point u,point v){
    return u.x*v.x+u.y*v.y+u.z*v.z;
}

point subt(point u,point v){
    point ret;
    ret.x=u.x-v.x;
    ret.y=u.y-v.y;
    ret.z=u.z-v.z;
    return ret;
}

double angle_cos(point u1,point u2,point v1,point v2){
    return dmult(subt(u1,u2),subt(v1,v2))/vlen(subt(u1,u2))/vlen(subt(v1,v2));
}

int main(){
    int i,j,cnt;
    while(scanf("%d",&n),n){
        for(i=0;i<n;i++){
            p[i].flag=false;
            scanf("%lf %lf %lf",&p[i].x,&p[i].y,&p[i].z);
        }
        scanf("%d",&m);
        for(i=0;i<m;i++){
            scanf("%lf %lf %lf %lf",&tel[i].tx,&tel[i].ty,&tel[i].tz,&tel[i].phi);
            tel[i].x=0.0;
            tel[i].y=0.0;
            tel[i].y=0.0;
        }
        cnt=0;
        for(i=0;i<m;i++){
            point tmp,ori;
            ori.x=0.0;
            ori.y=0.0;
            ori.z=0.0;
            tmp.x=tel[i].tx;
            tmp.y=tel[i].ty;
            tmp.z=tel[i].tz;
            for(j=0;j<n;j++){
                if(p[j].flag)continue;
                if(acos(angle_cos(ori,tmp,ori,p[j]))<tel[i].phi){
                    cnt++;
                    p[j].flag=true;
                }
            }
        }
        printf("%d\n",cnt);
    }
    return 0;
}

Uriel 2010-09-07 23:45 发表评论

POJ 1687 Buggy Sat---计算几何

Uriel — Sat, 28 Aug 2010 15:22:00 GMT

计算几何大水题，题目求哪个多边形是最外层的，只要算每个多边形的面积，找出最大的一个就行，因为题目给的坐标标号已经是排过序的，所以求面积只要叉积即可~

//Problem: 1687  User: Uriel
//Memory: 388K  Time: 0MS
//Language: G++  Result: Accepted
//Computational Geometry -> Cross Product
//2010.08.29

#include<math.h>
#include<stdio.h>
#include<stdlib.h>

struct point{
    int x,y;
}p[55];

int cross_product(point a,point b)
{
    return a.x*b.y-a.y*b.x;
}

int main(){
    int i,j,cse,n,x,y,tmp,id,pre,sum,maxs,ans;
    scanf("%d",&cse);
    while(cse--){
        scanf("%d",&n);
        for(i=0;i<n;i++){
            scanf("%d %d",&p[i].x,&p[i].y);
        }
        scanf("%d",&x);
        maxs=-1;ans=0;
        for(i=0;i<x;i++){
            scanf("%d",&y);
            sum=0;
            for(j=0;j<y;j++){
                scanf("%d",&id);
                id--;
                if(j==0){
                    pre=id;
                    tmp=id;
                }
                else{
                    sum+=cross_product(p[pre],p[id]);
                    pre=id;
                }
            }
            sum+=cross_product(p[id],p[tmp]);
            if(sum>maxs){
                maxs=sum;
                ans=i+1;
            }
        }
        printf("%d\n",ans);
    }
    return 0;
}

Uriel 2010-08-28 23:22 发表评论

POJ 1873 The Fortified Forest---凸包+枚举

Uriel — Wed, 05 May 2010 13:04:00 GMT

摘要: 惭愧，这题从中午写到晚上，一开始一直卡着调也调不了，搞到晚上没办法换了个凸包模板，然后历经各种NC错误终于AC。。 ... 阅读全文

Uriel 2010-05-05 21:04 发表评论

POJ 1244 Slots of Fun---几何

Uriel — Tue, 10 Nov 2009 11:15:00 GMT

真是挺水的一题。。以为有什么神奇的规律，看Discuss说坐标转换就使劲往那方面想。。一天多无果。。
然后O(n^3)暴力切掉。。
贴个代码。。纪念为了找规律草稿纸上N种坐标系的图。。

/*Problem: 1244  User: Uriel
   Memory: 784K  Time: 16MS
   Language: G++  Result: Accepted*/

#include<math.h>
#include<stdio.h>
#include<stdlib.h>
#include<string.h>

struct point
{
    double x[1000];
    double y[1000];
};

point P[27];
char s;
int i,j,k,h,layer,sum[27];
bool res;

double Dis(int x,int a,int b)
{
    return sqrt((P[x].x[a]-P[x].x[b])*(P[x].x[a]-P[x].x[b])+(P[x].y[a]-P[x].y[b])*(P[x].y[a]-P[x].y[b]));
}

bool Equal(double a,double b)
{
    return fabs(a-b)<0.0001;
}

int main()
{
    while(1)
    {
        memset(sum,0,sizeof(sum));
        memset(P,0x00,sizeof(P));
        res=false;
        scanf("%d",&layer);
        if(!layer)break;
        getchar();
        for(i=1;i<=layer;++i)
        {
            for(j=1;j<=i;j++)
            {
                scanf("%c",&s);
                P[s-'a'].x[sum[s-'a']]=-i-1+2*j;
                P[s-'a'].y[sum[s-'a']]=sqrt(3.0)*(layer-i);
                sum[s-'a']++;
            }
        }
        i=0;
flag:   if(i<26)
        {
            for(j=0;j<sum[i]-2;j++)
            {
                for(k=j+1;k<sum[i]-1;k++)
                {
                    for(h=k+1;h<sum[i];h++)
                    {
                        if(Equal(Dis(i,j,k),Dis(i,k,h))&&Equal(Dis(i,j,h),Dis(i,k,h))&&Equal(Dis(i,j,h),Dis(i,k,h)))
                        {
                            printf("%c",i+'a');
                            res=true;
                            i++;
                            goto flag;
                        }
                    }
                }
            }
            i++;
            goto flag;
        }
        if(!res)
        {
            printf("LOOOOOOOOSER!\n");
        }
        else
        {
            printf("\n");
        }
    }
    system("PAUSE");
    return 0;
}

Uriel 2009-11-10 19:15 发表评论

POJ 1410 Intersection---计算几何

Uriel — Thu, 01 Oct 2009 19:13:00 GMT

判断矩形与一线段是否相交，注意线段在矩形内部的情况
其他都是模板。。没什么了。。

/*Problem: 1410  User: Uriel
   Memory: 172K  Time: 0MS
   Language: C++  Result: Accepted*/

#include<stdio.h>
#include<stdlib.h>
#include <math.h>
#define eps 1e-8
#define zero(x) (((x)>0?(x):-(x))

struct point
{
    double x,y;
};

struct line
{
    point a,b;
};

double xmult(point p1,point p2,point p0)
{
    return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}

int dot_online_in(point p,point l1,point l2)
{
    return zero(xmult(p,l1,l2))&&(l1.x-p.x)*(l2.x-p.x)<eps&&(l1.y-p.y)*(l2.y-p.y)<eps;
}

int dots_inline(point p1,point p2,point p3)
{
    return zero(xmult(p1,p2,p3));
}

int same_side(point p1,point p2,point l1,point l2)
{
    return xmult(l1,p1,l2)*xmult(l1,p2,l2)>eps;
}

int intersect_in(point u1,point u2,point v1,point v2)
{
    if (!dots_inline(u1,u2,v1)||!dots_inline(u1,u2,v2))
    {
        return !same_side(u1,u2,v1,v2)&&!same_side(v1,v2,u1,u2);
    }
    return dot_online_in(u1,v1,v2)||dot_online_in(u2,v1,v2)||dot_online_in(v1,u1,u2)||dot_online_in(v2,u1,u2);
}

point P1,P2,S1,S2,S3,S4;

int main()
{
    int t,flag;
    double temp;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%lf %lf %lf %lf %lf %lf %lf %lf",&P1.x,&P1.y,&P2.x,&P2.y,&S1.x,&S1.y,&S2.x,&S2.y);
        if(S1.x>S2.x)
        {
            temp=S1.x;
            S1.x=S2.x;
            S2.x=temp;
        }
        if(S1.y<S2.y)
        {
            temp=S1.y;
            S1.y=S2.y;
            S2.y=temp;
        }
        S3.x=S1.x;
        S3.y=S2.y;
        S4.x=S2.x;
        S4.y=S1.y;
        flag=0;
        if(intersect_in(P1,P2,S1,S3) || intersect_in(P1,P2,S1,S4) || intersect_in(P1,P2,S2,S3) || intersect_in(P1,P2,S2,S4))
        {
            flag=1;
        }
        if(P1.x>=S1.x && P1.x<=S2.x && P2.x>=S1.x && P2.x<=S2.x && P1.y<=S1.y && P1.y>=S2.y && P2.y<=S1.y && P2.y>=S2.y)
        {
            flag=1;
        }
        if(!flag)
        {
            printf("F\n");
        }
        else
        {
            printf("T\n");
        }
    }
//    system("PAUSE");
    return 0;
}

Uriel 2009-10-02 03:13 发表评论

PKU 1265 Area---计算几何

Uriel — Thu, 01 Oct 2009 19:01:00 GMT

求顶点为整数的任意多边形内部整点数(Pick定理)，边上整点数(GCD)，面积(叉积)

第一次知道Pick定理。。
注意。。叉积求出面积可能是负数，需要处理下

/*Problem: 1265  User: Uriel
   Memory: 172K  Time: 16MS
   Language: C++  Result: Accepted*/

#include<math.h>
#include<stdio.h>
#include<stdlib.h>
#define MAXN 110
#define ABS(x) (x<0?-x:x)
struct point
{
    int x,y;
};

point P[MAXN],T;

int t,n,i,EPoint,IPoint,Asum;

int cross_product(point a,point b)
{
    return a.x*b.y-a.y*b.x;
}

int GCD(int x,int y)
{
    if(y==0)return x;
    else
        return GCD(y,x%y);
}

int main()
{
    scanf("%d",&t);
    for(int i=1;i<=t;i++)
    {
        scanf("%d",&n);
        Asum=0;
        EPoint=0;
        P[0].x=0;
        P[0].y=0;
        for(int j=1;j<=n;j++)
        {
            scanf("%d %d",&T.x,&T.y);
            P[j].x=P[j-1].x+T.x;
            P[j].y=P[j-1].y+T.y;
            EPoint+=GCD(ABS(T.x),ABS(T.y));
            Asum+=cross_product(P[j-1],P[j]);
        }
        if(Asum<0)
        {
            Asum=-Asum;
        }
        IPoint=Asum/2+1-EPoint/2;
        printf("Scenario #%d:\n%d %d %.1f\n\n",i,IPoint,EPoint,((double)Asum)/2);
    }
//    system("PAUSE");
    return 0;
}

Uriel 2009-10-02 03:01 发表评论

POJ 1228 Grandpa's Estate---计算几何

Uriel — Fri, 25 Sep 2009 15:28:00 GMT

摘要: WA无数次啊。。。看了Discuss之后又WA数次才过。。。1. n<=5,一律输出NO2. 所有点共线，NO3. 凸包上每条边上至少有三点，否则NO /**//*Problem: 1228 User: Uriel Memory: 616K Time: 0MS&... 阅读全文

Uriel 2009-09-25 23:28 发表评论

POJ 1113 Wall---计算几何

Uriel — Sat, 19 Sep 2009 16:53:00 GMT

简单凸包。。凸包周长加半径L圆的周长。。

/*Problem: 1113  User: Uriel
   Memory: 636K  Time: 47MS
   Language: G++  Result: Accepted*/

#include<math.h>
#include<stdio.h>
#include <stdlib.h>

#define eps 1e-8
#define PI 3.141592653589793
#define zero(x) (((x)>0?(x):-(x))

struct point{
    double x,y;
};

point P[3001],convex[3001];
double res;
int N,L;

//计算cross product (P1-P0)x(P2-P0)
double xmult(point p1,point p2,point p0){
    return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}
//graham算法顺时针构造包含所有共线点的凸包,O(nlogn)
point p1,p2;
int graham_cp(const void* a,const void* b){
    double ret=xmult(*((point*)a),*((point*)b),p1);
    return zero(ret)?(xmult(*((point*)a),*((point*)b),p2)>0?1:-1):(ret>0?1:-1);
}
void _graham(int n,point* p,int& s,point* ch){
    int i,k=0;
    for (p1=p2=p[0],i=1;i<n;p2.x+=p[i].x,p2.y+=p[i].y,i++)
        if (p1.y-p[i].y>eps||(zero(p1.y-p[i].y)&&p1.x>p[i].x))
            p1=p[k=i];
    p2.x/=n,p2.y/=n;
    p[k]=p[0],p[0]=p1;
    qsort(p+1,n-1,sizeof(point),graham_cp);
    for (ch[0]=p[0],ch[1]=p[1],ch[2]=p[2],s=i=3;i<n;ch[s++]=p[i++])
        for (;s>2&&xmult(ch[s-2],p[i],ch[s-1])<-eps;s--);
}

//构造凸包接口函数,传入原始点集大小n,点集p(p原有顺序被打乱!)
//返回凸包大小,凸包的点在convex中
//参数maxsize为1包含共线点,为0不包含共线点,缺省为1
//参数clockwise为1顺时针构造,为0逆时针构造,缺省为1
//在输入仅有若干共线点时算法不稳定,可能有此类情况请另行处理!
//不能去掉点集中重合的点
int graham(int n,point* p,point* convex,int maxsize=1,int dir=1){
    point* temp=new point[n];
    int s,i;
    _graham(n,p,s,temp);
    for (convex[0]=temp[0],n=1,i=(dir?1:(s-1));dir?(i<s):i;i+=(dir?1:-1))
        if (maxsize||!zero(xmult(temp[i-1],temp[i],temp[(i+1)%s])))
            convex[n++]=temp[i];
    delete []temp;
    return n;
}

double Dis(int a,int b)
{
    return sqrt((convex[a].y-convex[b].y)*(convex[a].y-convex[b].y)+(convex[a].x-convex[b].x)*(convex[a].x-convex[b].x));
}

int main()
{
    scanf("%d %d",&N,&L);
    for(int i=0;i<N;i++)
    {
        scanf("%lf %lf",&P[i].x,&P[i].y);
    }
    int M=graham(N,P,convex,1,1);
    res=0;
    for(int i=0;i<M-1;i++)
    {
//        printf("*%.2f %.2f*\n",convex[i].x,convex[i].y);
        res+=Dis(i,i+1);
    }
    res+=Dis(0,M-1);
    res+=2*PI*L;
    printf("%.0f\n",res);
//    system("PAUSE");
    return 0;
}

Uriel 2009-09-20 00:53 发表评论

POJ 2079 Triangle---计算几何

Uriel — Sat, 19 Sep 2009 16:48:00 GMT

摘要: 求凸包上的点构成的三角形最大面积。。开始抄模板构造凸包然后O(n^3)...听说有人优化就过了。。但是自己无论怎么优化都TLE。。无奈去强大的旋转卡壳。。抄了那段之后终于过了。。旋转卡壳还有些不懂，也基本不会应用。。要好好看下TLE到死的代码。。 /**//*Problem: 2079 User: Gilhirith &... 阅读全文

Uriel 2009-09-20 00:48 发表评论

POJ 2653 Pick-up sticks---计算几何

Uriel — Fri, 18 Sep 2009 12:10:00 GMT

摘要: 判断线段相交问题。。用链表做的。。（大牛们忽略下文。。）第一次完全自己写链表。。搞了一整天。。但是真的好开心啊~~ /**//*Problem: 2653 User: Uriel Memory: 8352K Time: 750MS &... 阅读全文

Uriel 2009-09-18 20:10 发表评论

POJ 1106 Transmitters---计算几何

Uriel — Thu, 17 Sep 2009 12:36:00 GMT

最近一直在搞计算几何。。不过不是钻研。。而是临时抱佛脚应付比赛。。。- -||

判断给定圆心，半径的半圆最多覆盖多少点。。
叉积应用---叉积正负判断夹角是否180度以内。。

/*Problem: 1106  User: Uriel
   Memory: 176K  Time: 0MS
   Language: C++  Result: Accepted*/

#include<math.h>
#include<stdio.h>
#include<stdlib.h>
#define MAXN 1001

typedef struct point{
    int x,y;
};

point P[MAXN],a,t;
int MAX;
double r;

int cross_product(point a,point c,point b)
{
    return (a.x-b.x)*(c.y-b.y)-(c.x-b.x)*(a.y-b.y);
}

int In_Circle(point a,point b,double r)
{
    if((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)<=r*r)
    {
        return 1;
    }
    return 0;
}

int main()
{
    int i,j,k,n;
    while(scanf("%d %d %lf",&a.x,&a.y,&r)!=EOF)
    {
        if(r<0)break;
        scanf("%d",&n);
        k=0;
        for(i=1;i<=n;i++)
        {
            scanf("%d %d",&t.x,&t.y);
            if(In_Circle(a,t,r))
            {
                P[++k].x=t.x;
                P[k].y=t.y;
            }
        }
        MAX=0;
        for(i=1;i<=k;i++)
        {
            int temp=0;
            for(j=1;j<=k;j++)
            {
                if(cross_product(P[i],a,P[j])<=0)
                {
                    temp++;
                }
            }
            if(2*temp<k)                //看Discuss才想到的
            {
                temp=k-temp;
            }
            if(temp>MAX)MAX=temp;
        }
        printf("%d\n",MAX);
    }
    system("PAUSE");
    return 0;
}

Uriel 2009-09-17 20:36 发表评论