极角排序:POJ 1696(叉积+深搜)

2023-12-02 06:18
文章标签 排序 poj 1696 叉积 极角

本文主要是介绍极角排序:POJ 1696(叉积+深搜),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

POJ1696题意:
一只很特殊的蚂蚁,只能向坐转,并且不能经过已经走过的路。一张地图上有n个食物让蚂蚁去采集,求蚂蚁经过所有食物的顺序(找出一条最长的非右拐的路径)。




样例:
Sample Input

2 (二次测试)
10 (十个食物点)
1 4 5 (点的编号,x坐标,y坐标)
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


分析:
每次找到一个基准点p0,最开始的时候是最左下角的点,然后我们再从未访问的点中找到相对于p0最小极角的那个点,比较极角用叉积来计算,每次迭代总是用此次循环中找到的最小极角来代替p0,然后再以这个新的p0,找到剩下的点中的最小极角,这样我们找到的就是一条最大的不右转的路径了。

下面是AC过的代码:
Java代码 复制代码 收藏代码
  1. import java.util.Scanner;
  2. import java.util.Arrays;
  3. class point implements Comparable {
  4. int id ;//食物原始的编号
  5. int ord ;//蚂蚁经过的顺序
  6. double x;//食物的x坐标
  7. double y;//食物的y坐标
  8. public int compareTo(Object o) {//升序排列
  9. return this.ord - ((point) o).ord;
  10. }
  11. }
  12. public class Main{
  13. point[] p;//所有的食物点
  14. boolean vis[] ;//标记是否被蚂蚁吃过了
  15. int n ;//食物的个数
  16. public Main(){
  17. }
  18. //向量(x1,y1),(x2,y2)的叉积
  19. public double det(double x1,double y1,double x2,double y2){
  20. return x1*y2 - x2*y1;
  21. }
  22. //小于0,说明向量ab的极角大于ac的极角
  23. double cross(point a,point b,point c){
  24. return det(b.x - a.x, b.y - a.y, c.x - a.x,c.y - a.y);
  25. }
  26. //求距离
  27. double dis(point a,point b) {
  28. double x1 = a.x;
  29. double y1 = a.y;
  30. double x2 = b.x;
  31. double y2 = b.y ;
  32. return (x2 - x1 )*(x2 - x1) + (y2-y1)*(y2 - y1);
  33. }
  34. //pre:蚂蚁开始的位置,num:蚂蚁经过的顺序
  35. void dfs(int pre,int num){
  36. if(num == n) return ;
  37. int mi = 1;
  38. while(vis[mi]) mi++;//找一个没有被蚂蚁吃过的点
  39. //在没有吃完的点中找出极角最小的点
  40. for(int i = mi + 1 ;i <= n;i++){
  41. if(vis[i]) continue ;//再找一个没有被蚂蚁吃过的点
  42. double w = cross(p[pre],p[mi],p[i]);//计算叉积
  43. if(w < 0) mi = i;//记录极角小的那个点,极角相同取距离近的点
  44. if((w ==0) && dis(p[pre],p[i]) < dis(p[pre],p[mi])) mi = i;
  45. }
  46. p[mi].ord = num + 1;//记录吃的顺序
  47. vis[mi] = true ;//标记为已吃过
  48. dfs(mi,num + 1) ;//递归
  49. }
  50. public static void main(String args[]){
  51. Main ma=new Main();
  52. ma.go();
  53. }
  54. public void go(){
  55. Scanner in=new Scanner(System.in);
  56. int t;
  57. t=in.nextInt();
  58. while(t-->0){
  59. n=in.nextInt();
  60. int mi = 1 ;
  61. vis=new boolean[n+1];
  62. p=new point[n+1];
  63. p[0]=new point();//这个点无用。
  64. for(int i = 1 ; i<= n;i++){//从命令行输入食物信息
  65. p[i]=new point();
  66. p[i].id=in.nextInt();//食物的编号
  67. p[i].x=in.nextDouble();//x坐标
  68. p[i].y=in.nextDouble(); //y坐标
  69. //找y最小的食物点,也即蚂蚁应该第一个吃的点
  70. if(p[mi].y > p[i].y)
  71. mi = i;
  72. }   

这篇关于极角排序:POJ 1696(叉积+深搜)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/444271

相关文章

【数据结构】——原来排序算法搞懂这些就行,轻松拿捏

前言:快速排序的实现最重要的是找基准值,下面让我们来了解如何实现找基准值 基准值的注释:在快排的过程中,每一次我们要取一个元素作为枢纽值,以这个数字来将序列划分为两部分。 在此我们采用三数取中法,也就是取左端、中间、右端三个数,然后进行排序,将中间数作为枢纽值。 快速排序实现主框架: //快速排序 void QuickSort(int* arr, int left, int rig

usaco 1.3 Mixing Milk (结构体排序 qsort) and hdu 2020(sort)

到了这题学会了结构体排序 于是回去修改了 1.2 milking cows 的算法~ 结构体排序核心: 1.结构体定义 struct Milk{int price;int milks;}milk[5000]; 2.自定义的比较函数,若返回值为正,qsort 函数判定a>b ;为负,a<b;为0,a==b; int milkcmp(const void *va,c

poj 3974 and hdu 3068 最长回文串的O(n)解法(Manacher算法)

求一段字符串中的最长回文串。 因为数据量比较大,用原来的O(n^2)会爆。 小白上的O(n^2)解法代码:TLE啦~ #include<stdio.h>#include<string.h>const int Maxn = 1000000;char s[Maxn];int main(){char e[] = {"END"};while(scanf("%s", s) != EO

hdu 2602 and poj 3624(01背包)

01背包的模板题。 hdu2602代码: #include<stdio.h>#include<string.h>const int MaxN = 1001;int max(int a, int b){return a > b ? a : b;}int w[MaxN];int v[MaxN];int dp[MaxN];int main(){int T;int N, V;s

poj 1511 Invitation Cards(spfa最短路)

题意是给你点与点之间的距离,求来回到点1的最短路中的边权和。 因为边很大,不能用原来的dijkstra什么的,所以用spfa来做。并且注意要用long long int 来存储。 稍微改了一下学长的模板。 stack stl 实现代码: #include<stdio.h>#include<stack>using namespace std;const int M

poj 3259 uva 558 Wormholes(bellman最短路负权回路判断)

poj 3259: 题意:John的农场里n块地,m条路连接两块地,w个虫洞,虫洞是一条单向路,不但会把你传送到目的地,而且时间会倒退Ts。 任务是求你会不会在从某块地出发后又回来,看到了离开之前的自己。 判断树中是否存在负权回路就ok了。 bellman代码: #include<stdio.h>const int MaxN = 501;//农场数const int

poj 1258 Agri-Net(最小生成树模板代码)

感觉用这题来当模板更适合。 题意就是给你邻接矩阵求最小生成树啦。~ prim代码:效率很高。172k...0ms。 #include<stdio.h>#include<algorithm>using namespace std;const int MaxN = 101;const int INF = 0x3f3f3f3f;int g[MaxN][MaxN];int n

poj 1287 Networking(prim or kruscal最小生成树)

题意给你点与点间距离,求最小生成树。 注意点是,两点之间可能有不同的路,输入的时候选择最小的,和之前有道最短路WA的题目类似。 prim代码: #include<stdio.h>const int MaxN = 51;const int INF = 0x3f3f3f3f;int g[MaxN][MaxN];int P;int prim(){bool vis[MaxN];

poj 2349 Arctic Network uva 10369(prim or kruscal最小生成树)

题目很麻烦,因为不熟悉最小生成树的算法调试了好久。 感觉网上的题目解释都没说得很清楚,不适合新手。自己写一个。 题意:给你点的坐标,然后两点间可以有两种方式来通信:第一种是卫星通信,第二种是无线电通信。 卫星通信:任何两个有卫星频道的点间都可以直接建立连接,与点间的距离无关; 无线电通信:两个点之间的距离不能超过D,无线电收发器的功率越大,D越大,越昂贵。 计算无线电收发器D

poj 1502 MPI Maelstrom(单源最短路dijkstra)

题目真是长得头疼,好多生词,给跪。 没啥好说的,英语大水逼。 借助字典尝试翻译了一下,水逼直译求不喷 Description: BIT他们的超级计算机最近交货了。(定语秀了一堆词汇那就省略吧再见) Valentine McKee的研究顾问Jack Swigert,要她来测试一下这个系统。 Valentine告诉Swigert:“因为阿波罗是一个分布式共享内存的机器,所以它的内存访问