poj1696 Space Ant(半平面交)
来源:互联网 发布:qq群优化排名多少钱 编辑:程序博客网 时间:2024/06/16 20:42
Description
The most exciting space discovery occurred at the end of the 20th century. In 1999, scientists traced down an ant-like creature in the planet Y1999 and called it M11. It has only one eye on the left side of its head and just three feet all on the right side of its body and suffers from three walking limitations:
The pictures transmitted by the Discovery space ship depicts that plants in the Y1999 grow in special points on the planet. Analysis of several thousands of the pictures have resulted in discovering a magic coordinate system governing the grow points of the plants. In this coordinate system with x and y axes, no two plants share the same x or y.
An M11 needs to eat exactly one plant in each day to stay alive. When it eats one plant, it remains there for the rest of the day with no move. Next day, it looks for another plant to go there and eat it. If it can not reach any other plant it dies by the end of the day. Notice that it can reach a plant in any distance.
The problem is to find a path for an M11 to let it live longest.
Input is a set of (x, y) coordinates of plants. Suppose A with the coordinates (xA, yA) is the plant with the least y-coordinate. M11 starts from point (0,yA) heading towards plant A. Notice that the solution path should not cross itself and all of the turns should be counter-clockwise. Also note that the solution may visit more than two plants located on a same straight line.
- It can not turn right due to its special body structure.
- It leaves a red path while walking.
- It hates to pass over a previously red colored path, and never does that.
The pictures transmitted by the Discovery space ship depicts that plants in the Y1999 grow in special points on the planet. Analysis of several thousands of the pictures have resulted in discovering a magic coordinate system governing the grow points of the plants. In this coordinate system with x and y axes, no two plants share the same x or y.
An M11 needs to eat exactly one plant in each day to stay alive. When it eats one plant, it remains there for the rest of the day with no move. Next day, it looks for another plant to go there and eat it. If it can not reach any other plant it dies by the end of the day. Notice that it can reach a plant in any distance.
The problem is to find a path for an M11 to let it live longest.
Input is a set of (x, y) coordinates of plants. Suppose A with the coordinates (xA, yA) is the plant with the least y-coordinate. M11 starts from point (0,yA) heading towards plant A. Notice that the solution path should not cross itself and all of the turns should be counter-clockwise. Also note that the solution may visit more than two plants located on a same straight line.
Input
The first line of the input is M, the number of test cases to be solved (1 <= M <= 10). For each test case, the first line is N, the number of plants in that test case (1 <= N <= 50), followed by N lines for each plant data. Each plant data consists of three integers: the first number is the unique plant index (1..N), followed by two positive integers x and y representing the coordinates of the plant. Plants are sorted by the increasing order on their indices in the input file. Suppose that the values of coordinates are at most 100.
Output
Output should have one separate line for the solution of each test case. A solution is the number of plants on the solution path, followed by the indices of visiting plants in the path in the order of their visits.
Sample Input
2
10
1 4 5
2 9 8
3 5 9
4 1 7
5 3 2
6 6 3
7 10 10
8 8 1
9 2 4
10 7 6
14
1 6 11
2 11 9
3 8 7
4 12 8
5 9 20
6 3 2
7 1 6
8 2 13
9 15 1
10 14 17
11 13 19
12 5 18
13 7 3
14 10 16
Sample Output
10 8 7 3 4 9 5 6 2 1 10
14 9 10 11 5 12 8 7 6 13 4 14 1 3 2
[Submit] [Go Back] [Status] [Discuss]
分析:
这道题的思路还算好想,但是比较难调
实际上就是一圈一圈的求凸包
但是连接点的处理比较麻烦,直接把我搞崩溃了。。。
tip
数组下标要搞明白了。。。
#include<cstdio>#include<cstring>#include<iostream>#include<algorithm>#include<cmath>using namespace std;const int N=10010;const double eps=1e-7;struct node{ int x,y; int id; node (int xx=0,int yy=0) { x=xx; y=yy; }};node p[N],sta[N],q[N];int n,m,cnt,ans[N];bool mark[N];node operator -(const node &A,const node &B) {return node(A.x-B.x,A.y-B.y);}int cmp(const node &A,const node &B){ if (A.x!=B.x) return A.x<B.x; else return A.y<B.y;}int Cross(node A,node B) {return A.x*B.y-A.y*B.x;}void TuB(){ m=0; sort(p,p+n,cmp); if (n==1) { sta[m++]=p[0]; return; } for (int i=0;i<n;i++) { while (m>1&&Cross(sta[m-1]-sta[m-2],p[i]-sta[m-2])<0) m--; //会经过一条直线上的多个点 sta[m++]=p[i]; } int k=m; for (int i=n-2;i>=0;i--) { while (m>k&&Cross(sta[m-1]-sta[m-2],p[i]-sta[m-2])<0) m--; sta[m++]=p[i]; }}int main(){ int T; scanf("%d",&T); while (T--) { scanf("%d",&n); memset(p,0,sizeof(p)); memset(q,0,sizeof(q)); memset(ans,0,sizeof(ans)); memset(mark,0,sizeof(mark)); cnt=0; int minx=100000,miny=100000,minp=0; for (int i=0;i<n;i++) { int t; scanf("%d",&t); scanf("%d%d",&p[i].x,&p[i].y); p[i].id=i; if (p[i].y==miny) //找最下方的一个点(y最小) { if (p[i].x<minx) minx=p[i].x,miny=p[i].y,minp=p[i].id; } else if (p[i].y<miny) miny=p[i].y,minx=p[i].x,minp=p[i].id; } while (1) { TuB(); int now=0; if (m==1) //凸包中就一个点 { ans[++cnt]=sta[0].id+1; break; } for (int i=0;i<m;i++) if (sta[i].id==minp&&i!=0) //起点 { if (i!=1) for (int j=i;j<m;j++) ans[++cnt]=sta[j].id+1,mark[sta[j].id+1]=1; else for (int j=i;j<m-1;j++) ans[++cnt]=sta[j].id+1,mark[sta[j].id+1]=1; now=i; break; } minp=sta[now-1].id; //minp是下一次的起点,所以我们把这一点扣去 for (int i=1;i<now-1;i++) ans[++cnt]=sta[i].id+1,mark[sta[i].id+1]=1; if (m==3) { ans[++cnt]=minp+1; break; } int k=-1; for (int i=0;i<n;i++) if (!mark[p[i].id+1]) q[++k]=p[i]; n=k; for (int i=0;i<=n;i++) p[i]=q[i]; n++; if (k==-1) break; } printf("%d ",cnt); for (int i=1;i<=cnt;i++) printf("%d ",ans[i]); printf("\n"); } return 0;}
阅读全文
0 0
- poj1696 Space Ant(半平面交)
- POJ1696-Space Ant(凸包)
- POJ1696-Space Ant
- poj1696 Space Ant
- POJ1696-Space Ant
- Poj1696 Space Ant
- poj1696 Space Ant
- poj1696 Space Ant 计算几何
- POJ1696---Space Ant(基础计算几何:极角排序)
- poj1696 Space Ant (极角排序,叉积)
- poj1696 Space Ant 点积,叉积,夹角
- POJ1696 Space Ant 【凸包-ToLeft-Test】
- 半平面交(未完成)
- poj1474(半平面交)
- poj3335(半平面交)
- BZOJ2618 (半平面交)
- 半平面交(入门)
- POJ1696 space ant(计算几何,极角排序)
- 问题和答案-内存泄露
- [LeetCode] 75. Sort Colors
- ssl1202-滑雪【记忆化搜索法】
- 让天堂的归天堂,让尘土的归尘土——谈Linux的总线、设备、驱动模型
- 似然函数个人理解
- poj1696 Space Ant(半平面交)
- ftp 简介及使用
- office365使用教程及简介
- 同步和异步
- jupyter notebook MemoryErrory
- 11. 方法 【连载 11】
- MDK 编译文件大小说明(Program Size)
- 冒泡排序法
- 从似然函数到EM算法