Jarvis模板Pascal

最新推荐文章于 2025-10-01 23:21:54 发布

原创最新推荐文章于 2025-10-01 23:21:54 发布 · 310 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#Jarvis

模板专栏收录该内容

6 篇文章

订阅专栏

本文详细介绍了一种经典的计算几何算法——Jarvis凸包算法，并提供了Pascal语言实现的源代码。文中首先给出了求解点间叉积的函数定义，进而通过逐步扫描点集的方法来确定构成凸包的顶点序列。

简介：

真的只是一个模板而已啊！！！还是Pascal的，P党水一波。

本文导读：

这只是作者是一个喜欢乱搞的人，A不了题目就直接把模板贴在了博客上，不过作者也希望可以帮到你啦，也没什么好稀奇的啦其实，就是一个Jarvis 算法。还把叉积公式另外放出来了………………………………

送上一首好听的歌《Nevada》，祝你一天好心情。

那么开始吧！！！

因为我是一个很Lazy的蒟蒻，所以把叉积公式放这里啦：

function cj(p0,p1,p2:longint):longint;
begin
  cj:=(a[p1,1]-a[p0,1])*(a[p2,2]-a[p0,2])-(a[p2,1]-a[p0,1])*(a[p1,2]-a[p0,2]);
end;

然后接下来就是Jarvis：

const
 maxn=100;
var
 a:array[0..maxn,0..2] of longint;
 n,i,rrrr:longint;
 ans:real;
function cj(p0,p1,p2:longint):longint;
begin
  cj:=(a[p1,1]-a[p0,1])*(a[p2,2]-a[p0,2])-(a[p2,1]-a[p0,1])*(a[p1,2]-a[p0,2]);
end;
procedure PROC1;
var
 x,y,i,j,k,mark,p:longint;
 longth:real;
begin
 x:=maxlongint;
 y:=maxlongint;
 for i:=1 to n do
 if (a[i,1]<x)or(a[i,1]=x)and(a[i,2]<y)then
  begin
   mark:=i;
   x:=a[i,1];
   y:=a[i,2];
  end;
 k:=mark;

 repeat
  i:=0;
   for j:=1 to n do
   if j<>k then
    begin
     if i<>0 then p:=cj(k,i,j);
     if (p<0)or (i=0)then
      begin
       longth:=sqrt(sqr(a[k,1]-a[j,1])+sqr(a[k,2]-a[j,2]));
       i:=j;
      end;
    end;
    ans:=ans+longth;
   k:=i;
 until k=mark;
end;
begin
 read(n);
 for i:=1 to n do read(a[i,1],a[i,2]);
 PROC1;
 write(ans:0:2);
end.

蒟蒻的注解：mark就是标记的意思，上过B站的肯定都知道吧，主要就是找一个在Ay最小的前提下Ax最小的点A，作为基点再Mark一下。之后通过对点的扫描实现数学模型上的旋转。什么时候结束呢？就是当枚举的点与mark的点重合的时候啦但是蒟蒻要问一下，凸包有什么现实意义吗？？？？请您在评论区告诉本蒟蒻，谢谢啦。