算法笔记 数学问题

最大公约数与最小公倍数

原题链接
原理：
最大公约数：递归 欧几里得算法
最小公倍数：得到a,b的最大公约数d 最小公倍数=a/d*b;

#include <iostream>
#include <cstdio>

int gcb(int a,int b)  //求最大公约数
{
    if (b == 0) return a;
    else return gcb(b, a % b);
}
int main()
{
   int n,m,x;
  while( scanf("%d",&n)!=EOF)
  {
     for(int i=1;i<=n;i++)
   {
       int d=1;
       scanf("%d",&m);
       for (int i=1; i<=m; i++ )
       {
           scanf("%d",&x);
           d = d/gcb(d,x)*x;     //求最大公倍
       }
       printf("%d\n",d);
   }
  }
    return 0;
}

大整数 四则运算

#include <iostream>
#include <cstdio>
#include <cstring>
struct bign
{
    int d[1000];
    int len;
    bign()
    {
        memset(d,0,sizeof(d));
        len=0;
    }
};
bign change(char str[])   //将整数转化为bign
{
    bign a;
    a.len=strlen(str);
    for(int i=0; i<a.len; i++)
    {
        a.d[i]=str[a.len-i-1]-'0';
    }
    return a;
}

int compare(bign a,bign b)
{
    if(a.len>b.len) return 1;
    else if(a.len < b.len) return -1;
    else
    {
        for(int i=a.len-1; i>=0; i--)
        {
            if(a.d[i]>b.d[i]) return 1;
            else if(a.d[i]<b.d[i]) return -1;
        }
        return 0;
    }
}

bign add(bign a,bign b)
{
    bign c;
    int carry=0;
    for(int i=0; i<a.len||i<b.len; i++) //以较长的为界限
    {
        int temp=a.d[i]+b.d[i]+carry;
        c.d[c.len++]=temp%10;
        carry=temp/10;
    }
    if(carry!=0)
    {
        c.d[c.len++] = carry;
    }
    return c;
}

bign sub(bign a,bign b)
{
    bign c;
    int len=0;
    for(int i=0; i<a.len||i<b.len; i++)
    {
        if(a.d[i]<b.d[i])
        {
            a.d[i+1]--;
            a.d[i]=10+a.d[i];
        }
        c.d[c.len++]=a.d[i]-b.d[i];
    }
    while(c.len-1>=1&&c.d[c.len-1] ==0)  //去除最高位 同时至少保留一位最低位
    {
        c.len--;
    }
    return c;
}

bign multi(bign a,int b) //高精度与高精度的乘法
{
    bign c;
    int carry=0;
    for(int i=0; i<a.len; i++)
    {
        int temp=a.d[i]*b+carry;
        c.d[c.len++]=temp%10;
        carry=temp/10;
    }
    while(carry!=0)         //乘法的进位可能不只一位
    {
        c.d[c.len++]=carry;
    }
}

bign divide(bign a,int b,int& r)
{
    bign c;
    c.len=a.len; //被除数的每一位和商的每一位都是一一对应的
    for(int i=a.len-1; i>=0; i--)
    {
        //始终是在和上一位余数加本位的和在作比较
        r=r*10+a.d[i];
        if(r<b) c.d[i]=0;
        else
        {
            c.d[i]=r/b; //商
            r=r%b; //余数
        }
    }
    //处理高位0
    while(c.len-1>=1&&c.d[c.len-1]==0)
    {
        c.len--;
    }
    return c;
}
void printf(bign a)
{
    for(int i=a.len-1; i>=0; i--)
    {
        printf("%d",a.d[i]);
    }
}
int main()
{
    char str1[1000],str2[1000];
    scanf("%s%s",str1,str2);
    bign a=change(str1);
    bign b=change(str2);
    printf(add(a,b));
    printf("\n");
    if(compare(a,b)==-1)
    {
        printf("\n-");
        printf(sub(b,a));
    }
    else
        printf(sub(a,b));

    return 0;
}