c算法大全分析和总结.docxVIP

一、数论算法 求两数的最大公约数 function gcd(a,b:integer):intege r; begin if b=0 then gcd:=a else gcd:=gcd (b,a mod b); end ; 求两数的最小公倍数 function lcm (a,b:intege r):integer; begin if ab then swap(a,b); lcm:=a; w hile lcm mod b0 do inc(lcm,a); end; 素数的求法 小范围内判断一个数是否为质数： function prime (n: intege r): Boolean; v ar 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; 判断 longint 范围内的数是否为素数（包含求 50000 以内的素数表）： procedure getprime; v ar i,j:longint; p:array[1..50000] of boolean; begin fillchar(p,sizeof(p),true); p[1]:=false; i:=2; w hile i50000 do begin if p[i] then begin j:=i*2; w hile j50000 do begin p[j]:=false; inc(j,i); end; end; inc(i); end; l:=0; for i:=1 to 50000 do if p[i] 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[i]=x the n break else if x mod pr[i]=0 then exit; prime:=true; end;{prime} 二、图论算法 最小生成树 Prim 算法： procedure prim (v 0:integer); var lowcost,closest:array [1..maxn] of integer; i,j,k,min:integer; begin for i:=1 to n do begin low cost[i]:=cost[v 0,i]; closest[i]:=v 0; end; for i:=1 to n-1 do begin {寻找离生成树最近的未加入顶点 k} min:=maxlongint; for j:=1 to n do if (low cost[j]min) and (low cost[j]0) then begin min:=low cost[j]; k:=j; end; low cost[k]:=0; {将顶点 k 加入生成树} {生成树中增加一条新的边 k 到 closest[k]} {修正各点的 low cost 和 clo sest 值} for j:=1 to n do if cost[k,j]lw ocost[j] then begin lowcost[j]:=cost[k,j]; closest[j]:=k; end; end; end;{prim} Kruskal 算法：(贪心) 按权值递增顺序删去图中的边，若不形成回路则将此边加入最小生成树。 function find(v :integer):integer; {返回顶点 v 所在的集合} v ar i:integer; begin i:=1; w hile (i=n) and (not v in v set[i]) do inc(i); if i=n then find:=i else find:=0; end; procedure kruskal; v ar tot,i,j:integer; begin for i:=1 to n do v set[i]:=[i];{初始化定义 n 个集合，第 I 个集合包含一个元素 I} p:=n-1; q:=1; tot:=0; {p 为尚待加入的边数， q 为边集指针} sort; {对所有边按权值递增排序，存于 e[I]中， e[I].v 1 与 e[I].v 2 为边I 所连接的两个顶点的序号， e [I].len 为第I 条边的长度} w hile p0 do begin i:=find(e[q].v 1);j:=find(