首页 > 编程语言 >论高精度算法

论高精度算法

时间:2023-07-21 15:48:01浏览次数:27  
标签:string 高精度 int void -- 算法 s1

一、概念:高精度也可以称之为大整数,我们对于超出整型(int)甚至是(long long)数据范围的数称为高精度数。(注int范围:-2147483648~2147483647      long long范围:-9223372036854775808~9223372036854775808)

二、用途:对于这些高精度数计算机无法快速,准确无误地算出来,甚至有时无法正常存储,这时候就要用到高精度算法。

三、基本方法:

(1)、输入与存储

数据的输入,需要考虑位数问题,很可能会有上千位,同时也要考虑输入的数字会带字母,如20进制,18进制等等。所以我们采用字符串输入或字符数组输入。下面提供两种方式:

string a,b;
cin>>a>>b;
//方式一
char a[n];
char b[m];
for(int i=1;i<=n;i++)
{
    cin>>a[i];
}
for(int i=1;i<=m;i++)
{
    cin>>b[i];
}
//方式二,使用时n,m根据不同情况取不同值,一定要保证有足够的分配空间

(2) 、数据的转换与存储

输入进来的并不是数,而是字符串,不能直接进行运算。所以我们需要将字符串打散,微分成整数数组。

比如说,在计算123+4567时。首先输入a,b字符串,内容分别为123,4567。接着逆序存入数组储存:

 代码就大概长这个笋样:

string a,b;//字符串类型接收 
int n[10005],m[10005];
void tr(){
    int len1=a.length();
    n[0]=len1;//将输入数据的位数存储在整型数组的第一个元素,便于开启下面的循环 
    for(int i=0;i<len1;i++){
        n[i+1]=n[len1-1-i]-'0';//逆序存储 
    }
    int len2=b.length();
    m[0]=len1;
    for(int i=0;i<len2;i++){
        m[i+1]=m[len1-1-i]-'0';//逆序存储 
    }
}
char a[1010],b[1010];
int n[1010],m[1010];
void tr(){
    int len1=st(a);
    n[0]=len1;//将输入数据的位数存储在整型数组的第一个元素,便于开启下面的循环
    for(int i=0;i<len1;i++)
    {
        n[i+1]=n[len1-1-i]-'0';//逆序存储
    }
    int len2=st(b);
    m[0]=len2;
    for(int i=0;i<len2;i++)
    {
        m[i+1]=m[len1-1-i]-'0';//逆序存储 
    }
} 

通过这种方式,我们就能将输入的数据很好地放入数组中了。

(3)数组的计算,输出

高精度算法主要有加法、减法、乘法和除法。

通过类似竖式进行计算,以达到快速运算的目的。

(注:输出时去掉下标为0的数组)

四、计算代码

(1)加法

 

 

 将两个正整数相加,c[i]=a[i]+b[i](注:数组下标为0表示数组个数,下标为一的数才是个位),对应的位置相加,得到的新数字>9时,留下个位,向高位进1,模拟加法的竖式运算。

 

int add(int a[],int b[]) //计算a=a+b
{
int i,k; k=a[0]>b[0]?a[0]:b[0]; //k是a和b中位数最大的一个的位数 for(i=1;i<=k;i++) {a[i+1]+=(a[i]+b[i])/10; //若有进位,则先进位 a[i]=(a[i]+b[i])%10;} //计算当前位数字,注意:这条语句与上一条不能交换。 if(a[k+1]>0) a[0]=k+1; //修正新的a的位数(a+b最多只能的一个进位) else a[0]=k; return 0; }

(2)减法

同理我们得到减法运算程序,将加法的逢十进一换为不够借一当十。

int add(int a[],int b[]);//计算a=a-b,返加符号位0:正数 1:负数
{
 int flag,i
  flag=compare(a,b); //调用比较函数判断大小
if (falg==0)//相等
  {memset(a,0,sizeof(a));return 0;}  //若a=b,则a=0,也可在return前加一句a[0]=1,表示是 1位数0
if(flag==1) //大于  
  {  for(i=1;i<=a[0];i++)
      {  if(a[i]<b[i]){ a[i+1]--;a[i]+=10;} //若不够减则向上借一位
        a[i]=a[i]-b[i];}
     while(a[a[0]]==0) a[0]--; //修正a的位数
    return 0;}
if (flag==-1)//小于  则用a=b-a,返回-1
    { for(i=1;i<=b[0];i++)       {  if(b[i]<a[i]){ b[i+1]--;b[i]+=10;} //若不够减则向上借一位
        a[i]=b[i]-a[i];}
      a[0]=b[0];
     while(a[a[0]]==0) a[0]--; //修正a的位数
    return -1;}
}

(3)乘法

乘法可以分为2种,一个是高精度乘单精度,一个是高精度乘高精度

高精度乘单精度

例如我们要算123*5,给定a=123,b=5,我们需要将高精度数每一位乘上单精度数,然后处理进位。

 代码如下:

#include <bits/stdc++.h>
using namespace std;
int main()
{
    int a[251]={};
    string s1;
    int b;
    cin>>s1>>b;
    for(int i=0;i<s1.size();i++){
        a[i]=s1[s1.size()-i-1]-'0';
    }
    for(int i=0;i<s1.size();i++){
        a[i]=a[i]*b;
    }
    for(int i=0;i<s1.size()+4;i++){
        if(a[i]>=10){
            a[i+1]+=a[i]/10;
            a[i]%=10;
        }
    }
    int point=0;
    for(int i=s1.size()+4;i>=0;i--){
        if(a[i]!=0){
            point=i;
            break;
        }
    }
    for(int i=point;i>=0;i--){
        cout<<a[i];
    }
    return 0;
}

高精度乘高精度

例如计算123*456

 我们通过竖式计算法,用一个数去乘另一个数的每一位(单精度),所以要用到双重循环。

#include <bits/stdc++.h>
using namespace std;
int main(){
 string s1,s2;
    int a[4002]={},b[4002]={},result[8004]={};
    cin>>s1>>s2;
    for(int i=0;i<s1.size();i++) a[i]=s1[s1.size()-i-1]-'0';
    for(int i=0;i<s2.size();i++) b[i]=s2[s2.size()-i-1]-'0';
    for(int i=0;i<s1.size();i++)
        {
        for(int j=0;j<s2.size();j++)
                {
            result[i+j]+=a[i]*b[j];
               result[i+j+1]+=result[i+j]/10,
            result[i+j]%=10;
        }
    }
     int p=0;
    for(int i=s1.size()+s2.size()-1;i>=0;i--)
        {
        if(result[i]>0){
            p=i;break;
        }
    }
    for(int i=p;i>=0;i--)
        {
        cout<<result[i];
    }
    return 0;
}      

(4)除法

最后也是我认为最烦的———除法。

根据上面的经验,我们仍采用竖式计算的方法

 我们让计算机进行通过减法进行试商,模拟我们日常除法,主要有以下几步:

1、将除数移动和被除数对齐,位数不够时,补0,

2、利用被除数减去除数,一直减到被除数小于除数,减的次数,就是“试商”的结果,每移动一次。

3、重复上述步骤,一直到被除数和除数的位数相等为止。

#include <iostream>
#include <cstring>
using namespace std;
const int N =1001;
int aa[N],bb[N],cc[N]; 
void inputNum(string ss,int a[]);
void printArr(int a[]);  
void jian(int a[],int b[]);
void jisuan(int a[],int b[],int c[]);
void movei(int a[],int b[],int i);
 
int compare (int a[],int b[]); 
 
int main()
{
    string s1 ="525353422";
    string s2 ="253532";
    inputNum(s1,aa);
    printArr(aa);
    
    inputNum(s2,bb);
    printArr(bb);
    
    if (compare(aa,bb) ==0) cout<<1<<endl;
    else if (compare(aa,bb))
    {
        jisuan(aa,bb,cc);
    }
    else if (compare(aa,bb)<0)
    {
        jisuan(bb,aa,cc);
    }
    
    printArr(cc);
    printArr(aa);
    return 0;
}
int compare (int a[],int b[])
{
    int i;
    if (a[0]>b[0]) return 1;
    if (a[0]<b[0]) return -1;
    for (i=a[0];i>0;i--)
    {
        if (a[i]>b[i]) return 1;
        if (a[i]<b[i]) return -1;
    }
    return 0;
}
 
void movei(int a[],int b[],int j)
{
    for (int i=1;i<=a[0];i++)
    {
        b[i+j-1] =a[i];
    }
    b[0] = a[0]+j-1;
}
 
void jisuan(int a[],int b[],int c[])
{
  c[0] = a[0]-b[0]+1;  
  int temp[1001];
  for (int i=c[0];i>0;i--) 
  {
      memset(temp,0,sizeof(temp)); 
      movei(b,temp,i); 
        while (compare(a,temp)>=0)
      {
          c[i]++; 
        jian(a,temp);
      }
  }
  int m=c[0];
  while (c[0]>0 && c[m]==0) 
  {
     c[0]--;
  }
    
}
 
 
void jian(int a[],int b[])
{
    for (int i=1;i<=a[0];i++)
    {
        if (a[i]<b[i])
        {
            a[i+1]--;
            a[i]+=10;
        }
        a[i]-=b[i];
    }
    int i=a[0];
    while (a[i]==0)
    {
        i--;
    }
    a[0]=i;
}
 
void inputNum(string ss,int a[])
{
    int len = ss.length();
    a[0] = len;
    for (int i=0;i<len;i++)
    {
        a[len-i] = ss[i] -48;
    }
}
void printArr(int a[])
{
 
    for (int i=a[0];i>0;i--)
    {
        cout<<a[i];
    }
    cout<<endl;
}

标签:string,高精度,int,void,--,算法,s1
From: https://www.cnblogs.com/ltrh/p/17571527.html

相关文章

  • Python数学建模系列(六):蒙特卡洛算法
    文章目录前言往期文章1、蒙特卡洛算法样例1样例2样例32、三门问题3、M*M豆问题结语前言Hello!小伙伴!非常感谢您阅读海轰的文章,倘若文中有错误的地方,欢迎您指出~ 自我介绍ଘ(੭ˊᵕˋ)੭昵称:海轰标签:程序猿|C++选手|学生简介:因C语言结识编程,随后转入计算机专业,有幸拿过一些国奖......
  • codility算法题:找出不在数组中的最小正整数
    1.题目读题   考查点 2.解法思路 代码逻辑 具体实现解法一:publicclassSolution{publicstaticvoidmain(String[]args){System.out.println(solution(newint[]{1,3,6,4,1,2}));System.out.println(solution(newint[]{1,......
  • codility算法题:猫过桥问题
    1.题目读题  考查点 2.解法思路 代码逻辑 具体实现 publicclassSolutions{publicstaticvoidmain(String[]args){System.out.println(solution(10,newint[]{2,3,4,8},newint[]{2,5}));System.out.println(solution(10,......
  • 请享用美味的快速幂算法-通俗易懂版
    一、算法整体思路第1步按照最直接、最好理解的方式看,2的n次幂是n个2相乘,即有如下公式例如:第2步然而为了节省大量时间,通过简单的思考和严格数学推理,我们不难理解以下结论: 1.偶数幂的情况:通过幂函数运算法则,有2n=(2n/2)2,即有如下等式:例如24......
  • python算法的时间复杂度怎么算
    项目方案:计算列表中元素的平方和1.项目背景在很多应用中,我们需要对一个列表中的元素进行一些计算操作。例如,计算一个列表中所有元素的平方和。这个项目方案就是要实现这样的功能。2.问题定义给定一个列表nums,计算列表中所有元素的平方和。即,对于列表中的每个元素num,计算nu......
  • bellman-ford算法理解
    bellman-ford算法理解从本题谈起再回归到最短路。本题为限制边数的最短路,是这个算法优势领域的题目。为什么它能解决?最外层每循坏一次,就是各点向外走一条边,内层对边的遍历是对所有边进行松弛操作,每次进行该操作时,需要用到备份数组,目的是防止连锁反应,保证每次每个点到起点的距离......
  • LntonCEC算法算力云平台服务通过EasyNTS内网穿透到公网上的具体操作流程
    算法算力云平台的主要特点包括高性能、高可靠性、高可扩展性和低成本。LntonCEC算法算力云平台是一种为用户提供高效、强大的算法计算服务的云计算平台。它可以帮助用户快速、灵活地运行各种复杂的计算模型和算法,包括机器学习、人工智能、大数据分析、图像识别等领域。算法算力云平......
  • 揭秘绿幕抠图算法技术
     绿幕抠图为什么是“绿幕”呢?人眼的感光系统和摄像机的感光芯片采集的色彩中,最常见的就是红、蓝、绿三原色。红色在演员服饰和物体中较为常见,不利于背景分离;绿色则是人体肤色最少的颜色,人眼对绿色最为敏感,而相机也是模仿人眼的设计,CMOS采集信息是按照RGGB,其信号最强、噪波最少......
  • 遗传算法 深度学习
    遗传算法与深度学习1.遗传算法和深度学习的概述遗传算法和深度学习是两种不同的优化算法,它们在解决问题时有着不同的应用场景和方法。遗传算法是一种通过模拟生物进化过程中的自然选择和遗传机制来寻找最优解的优化算法。而深度学习则是一种通过神经网络模型来学习和识别复杂的......
  • 1307:高精度乘法
    1307:【例1.3】高精度乘法时间限制:1000ms      内存限制:65536KB【题目描述】输入两个高精度正整数M和N(M和N均小于100位)。求这两个高精度数的积。【输入】输入两个高精度正整数M和N。【输出】求这两个高精度数的积。【输入样例】363【输出样例】......