几个简单数论的算法的实现

本文通过C++语言实现了几个经典算法,包括求最大公约数、最小公倍数及判断一个数是否为素数的方法。提供了完整的代码示例,并对算法进行了简单的说明。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

这学期上C++,无聊,技术又太烂,什么都做不了,所以就随便把一些算法用C++实现了,算法还是很经典的哈,不要应为我引用了,就不屑了哈~

素数的那两个是有点缺陷的哈,聪明的你一定可以发现哈^_^

 语言:C++

编译环境:VS2008

//求两个数的最大公约数
#include <iostream>

using namespace std;

int gcd(int a, int b)
{
 int g;
 if(b == 0)
  g = a;
 else
  g = gcd(b,a%b);
 return g;
}

int main()
{
 int a, b;
 cout << "请输入两个整数" << endl;
 cin >> a >> b;
 cout << gcd( a, b) << endl;
 system("pause");
 return 0;
}

 

//求两个数的最小公倍数
#include <iostream>

using namespace std;

int lcm(int a, int b)
{
 int temp, l;
 if(a < b)
 {
  temp = a;
  a = b;
  b = a;
 }
 l = a;
 while(l%b)//源程序是l%b > 0;
 {
  l += a;
 }
 return l;
}
int main()
{
 int a, b;
 cout << "请输入两个整数:" << endl;
 cin >> a >> b;
 cout << lcm(a, b) <<endl;
 system("pause");
 return 0;
}

 

//判断小范围内一个数是否为质数(素数)
#include <iostream>
#include <math.h>

using namespace std;

int prime(int n)
{
 int i;
 bool flag;
 for(i=2; i<sqrt((double)n); i++)
 {
  if(n%i == 0)
  {
   flag = false; break;
  }
  flag = true;
 }
 return flag;
}

int main()
{
 int n;
 cout << "请输入一个大于1的整数:" << endl;
 cin >> n;
 cout << prime(n) << endl;
 system("pause");
 return 0;
}

 

//判断LongInt范围内的数是否是素数(包含求50000以内的素数表)
#include <iostream>
#include <math.h>

using namespace std;

int main()
{
 long i, j, l, t;
 long pr[50000];
 bool a[50000], flag;
 a[1] = false;
 i = 2;
 while(i < 50000)
 {
  if(a[i])
   j = i * 2;
  while(j < 5000)
  {
   a[j] = false;
   j += i;
  }
  i++;
 }
 l = 0;
 for(i=1; i<50000; i++)
 {
  if(a[i])
   pr[l] = i;
  l++;
 }
 cin >> t;
 flag = false;
 for(i=0; i<l; i++)
 {
  if(pr[i] >= t)
   break;
  else
  {
   if(t%pr[i] == 0)
   {
    cout << flag << endl;
    system("pause");
    exit(0);
   }   
  }
 }
 flag = true;
 cout << flag << endl;
 system("pause");
 return 0;
}

算法大全(C,C++) 一、 数论算法 1.求两数的最大公约数 function gcd(a,b:integer):integer; begin if b=0 then gcd:=a else gcd:=gcd (b,a mod b); end ; 2.求两数的最小公倍数 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; 3.素数的求法 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[i] 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[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 then break else if x mod pr[i]=0 then exit; prime:=true; end;{prime} 二、图论算法 1.最小生成树 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[i]:=cost[v0,i]; closest[i]:=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[i]) 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[i]:=[i];{初始化定义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[i]:=vset[i]+vset[j];vset[j]:=[]; dec(p); end; inc(q); end; writeln(tot); end; 2.最短路径
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值