资源描述：
《常见算法的MATLAB实现八.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、常见算法的MATLAB实现八1.数论算法求两数的最大公约数functiongcd(a,b:integer):integer;beginifb=0thengcd:=aelsegcd:=gcd(b,amodb);end;求两数的最小公倍数functionlcm(a,b:integer):integer;beginifa<bthenswa(a,b);lcm:=a;whilelcmmodb>0doinc(lcm,a);end;素数的求法A.小范围内判断一个数是否为质数：functionrime(n:integer):Boolean;varI:in
2、teger;beginforI:=2totrunc(sqrt(n))doifnmodI=0thenbeginrime:=false;exit;end;rime:=true;end;B.判断longint范围内的数是否为素数（包含求50000以内的素数表）：roceduregetrime;vari,j:longint;:array[1..50000]ofboolean;beginfillchar(,sizeof(),true);[1]:=false;i:=2;whilei<50000dobeginifthenbeginj:=i*2;whilej&
3、lt;50000dobegin[j]:=false;inc(j,i);end;end;inc(i);end;l:=0;fori:=1to50000doifthenbegininc(l);r[l]:=i;end;end;{getrime}functionrime(x:longint):integer;vari:integer;beginrime:=false;fori:=1toldoifr>=xthenbreakelseifxmodr=0thenexit;rime:=true;end;{rime}2．3．4.求最小生成树A.rim算法：roced
4、urerim(v0:integer);varlowcost,closest:array[1..maxn]ofinteger;i,j,k,min:integer;beginfori:=1tondobeginlowcost:=cost[v0,i];closest:=v0;end;fori:=1ton-1dobegin{寻找离生成树最近的未加入顶点k}min:=maxlongint;forj:=1tondoif(lowcost[j]<min)and(lowcost[j]<>0)thenbeginmin:=lowcost[j];k:=j;e
5、nd;lowcost[k]:=0;{将顶点k加入生成树}{生成树中增加一条新的边k到closest[k]}{修正各点的lowcost和closest值}forj:=1tondoifcost[k,j]<lwocost[j]thenbeginlowcost[j]:=cost[k,j];closest[j]:=k;end;end;end;{rim}B.Kruskal算法：(贪心)按权值递增顺序删去图中的边若不形成回则将此边加入最小生成树。functionfind(v:integer):integer;{返回顶点v所在的集合}vari:integer;
6、begini:=1;while(i<=n)and(notvinvset)doinc(i);ifi<=nthenfind:=ielsefind:=0;end;rocedurekruskal;vartot,i,j:integer;beginfori:=1tondovset:=;{初始化定义n个集合第I个集合包含一个素I}:=n-1;q:=1;tot:=0;{为尚待加入的边数q为边集指针}sort;{对所有边按权值递增排序存于e[I]中e[I].v1与e[I].v2为边I所连接的两个顶点的序号e[I].len为第I条边的长度}while>
7、0dobegini:=find(e[q].v1);j:=find(e[q].v2);ifi<>jthenbegininc(tot,e[q].len);vset:=vset+vset[j];vset[j]:=[];dec();end;inc(q);end;writeln(tot);end;5.最短径A.标号法求解单源点最短径：vara:array[1..maxn,1..maxn]ofinteger;b:array[1..maxn]ofinteger;{b指顶点i到源点的最短径}mark:array[1..maxn]ofboolean;roc
8、edurebhf;varbest,best_j:integer;beginfillchar(mark,siz
