1.數論算法 
  求兩數的最大公約數 
  function gcd(a,b:integer):integer; 
  begin 
      if b=0 then gcd:=a 
      else gcd:=gcd (b,a mod b); 
  end ; 
    
  求兩數的最小公倍數 
  function lcm(a,b:integer):integer; 
  begin 
      if a< b then swap(a,b); 
      lcm:=a; 
      while lcm mod b >0 do inc(lcm,a); 
  end; 
    
  素數的求法 
  A.小范圍內判斷一個數是否為質數： 
  function prime (n: integer): Boolean; 
  var I: integer; 
  begin 
      for I:=2 to trunc(sqrt(n)) do 
          if n mod I=0 then 
          begin 
            prime:=false; exit; 
          end; 
      prime:=true; 
  end; 
    
  B.判斷longint范圍內的數是否為素數（包含求50000以內的素數表）： 
  procedure getprime; 
  var 
  i,j:longint; 
  p:array[1..50000] of boolean; 
  begin 
      fillchar(p,sizeof(p),true); 
      p[1]:=false; 
      i:=2; 
      while i< 50000 do 
      begin 
          if p then 
          begin 
            j:=i*2; 
            while j< 50000 do 
            begin 
                p[j]:=false; 
                inc(j,i); 
            end; 
          end; 
          inc(i); 
      end; 
      l:=0; 
      for i:=1 to 50000 do 
        if p then 
        begin 
            inc(l);
            pr[l]:=i; 
        end; 
  end;{getprime} 
  function prime(x:longint):integer; 
  var i:integer; 
  begin 
      prime:=false; 
      for i:=1 to l do 
        if pr >=x then break 
        else if x mod pr=0 then exit; 
      prime:=true; 
  end;{prime} 
    
  2． 
    
  3． 
    
    
  4.求最小生成樹 
  A.Prim算法： 
  procedure prim(v0:integer); 
  var 
  lowcost,closest:array[1..maxn] of integer; 
  i,j,k,min:integer; 
  begin 
      for i:=1 to n do 
        begin 
            lowcost:=cost[v0,i]; 
            closest:=v0; 
        end; 
      for i:=1 to n-1 do 
      begin 
          {尋找離生成樹最近的未加入頂點k} 
          min:=maxlongint; 
          for j:=1 to n do 
          if (lowcost[j]< min) and (lowcost[j]< >0) then 
          begin 
            min:=lowcost[j]; 
            k:=j; 
          end; 
          lowcost[k]:=0; {將頂點k加入生成樹} 
          {生成樹中增加一條新的邊k到closest[k]} 
          {修正各點的lowcost和closest值} 
          for j:=1 to n do 
          if cost[k,j]< lwocost[j] then 
          begin 
              lowcost[j]:=cost[k,j]; 
              closest[j]:=k; 
          end; 
      end; 
  end;{prim} 
B.Kruskal算法：(貪心) 
按權值遞增順序刪去圖中的邊，若不形成回路則將此邊加入最小生成樹。 
function find(v:integer):integer; {返回頂點v所在的集合} 
var i:integer; 
begin 
  i:=1; 
  while (i< =n) and (not v in vset) do inc(i); 
  if i< =n then find:=i 
  else find:=0; 
end; 
procedure kruskal; 
var 
tot,i,j:integer; 
begin 
  for i:=1 to n do vset:=;{初始化定義n個集合，第I個集合包含一個元素I} 
  p:=n-1; q:=1; tot:=0; {p為尚待加入的邊數，q為邊集指針} 
  sort; 
  {對所有邊按權值遞增排序，存于e[I]中，e[I].v1與e[I].v2為邊I所連接的兩個頂點的序號，e[I].len為第I條邊的長度} 
  while p >0 do 
  begin 
      i:=find(e[q].v1);j:=find(e[q].v2); 
      if i< >j then 
      begin 
          inc(tot,e[q].len); 
          vset:=vset+vset[j];vset[j]:=[]; 
          dec(p); 
      end; 
      inc(q); 
  end; 
  writeln(tot); 
end; 
    
    
  5.最短路徑 
  A.標號法求解單源點最短路徑： 
  var 
  a:array[1..maxn,1..maxn] of integer; 
  b:array[1..maxn] of integer; {b指頂點i到源點的最短路徑} 
  mark:array[1..maxn] of boolean; 
    
  procedure bhf; 
  var 
  best,best_j:integer; 
  begin 
      fillchar(mark,sizeof(mark),false); 
      mark[1]:=true; b[1]:=0;{1為源點} 
      repeat 
          best:=0; 
          for i:=1 to n do 
          If mark then {對每一個已計算出最短路徑的點} 
            for j:=1 to n do 
                if (not mark[j]) and (a[i,j] >0) then 
                    if (best=0) or (b+a[i,j]< best) then 
                    begin 
                      best:=b+a[i,j]; best_j:=j; 
                    end; 
                if best >0 then 
                begin 
                    b[best_j]:=best；mark[best_j]:=true; 
                end; 
      until best=0; 
  end;{bhf} 
    
  B.Floyed算法求解所有頂點對之間的最短路徑： 
  procedure floyed; 
  begin 
      for I:=1 to n do 
          for j:=1 to n do 
            if a[I,j] >0 then p[I,j]:=I else p[I,j]:=0; 
            {p[I,j]表示I到j的最短路徑上j的前驅結點} 
      for k:=1 to n do {枚舉中間結點} 
          for i:=1 to n do 
            for j:=1 to n do 
                if a[i,k]+a[j,k]< a[i,j] then 
                begin 
                    a[i,j]:=a[i,k]+a[k,j]; 
                    p[I,j]:=p[k,j]; 
                end; 
    end; 
C. Dijkstra 算法： 
類似標號法，本質為貪心算法。 
var 
a:array[1..maxn,1..maxn] of integer; 
b,pre:array[1..maxn] of integer; {pre指最短路徑上I的前驅結點} 
mark:array[1..maxn] of boolean; 
procedure dijkstra(v0:integer); 
begin 
  fillchar(mark,sizeof(mark),false); 
  for i:=1 to n do 
      begin 
          d:=a[v0,i]; 
          if d< >0 then pre:=v0 else pre:=0; 
      end; 
  mark[v0]:=true; 
  repeat {每循環一次加入一個離1集合最近的結點并調整其他結點的參數} 
      min:=maxint; u:=0; {u記錄離1集合最近的結點} 
      for i:=1 to n do 
          if (not mark) and (d< min) then 
          begin 
            u:=i; min:=d; 
          end; 
      if u< >0 then 
      begin 
          mark:=true; 
          for i:=1 to n do 
            if (not mark) and (a[u,i]+d< d) then 
            begin 
                d:=a[u,i]+d; 
                pre:=u; 
            end; 
      end; 
    until u=0; 
end; 
D.計算圖的傳遞閉包 
Procedure Longlink; 
Var 
T:array[1..maxn,1..maxn] of boolean; 
Begin 
  Fillchar(t,sizeof(t),false); 
  For k:=1 to n do 
      For I:=1 to n do 
          For j:=1 to n do 
            T[I,j]:=t[I,j] or (t[I,k] and t[k,j]); 
End;