高精度计算(三)压位

2024-01-30 21:58
文章标签 计算 高精度 压位

本文主要是介绍高精度计算(三)压位,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

高精压位

继承咱们(一)上的例题…A+B+C 链接

有兴趣的可以用上一篇的内容做一下这道神题

不出意外的话会出现下面的字样 Time Limit Exceeded ……

我们需要用一种特殊的技巧对上算法进行处理


压位思想

还记得我们存储数字的方式?

举个栗子:1234 + 1234这是之前的存储方式

实际上我们一个数组空间(int)是可以储存 2147483647 以下的数字 这样非常浪费 而且以上存储方式需要按位相加最后的计算次数是4

现在考虑这样的存储方式

这样相对上面的存储方式来说利用度较高 

理解方式:考虑二进制 十进制 十六进制 的计算方式 满2进1 满10进1 满16进1

那么以上的存储方式请理解为100进制的计算 也就是满100进1

按照这样储存 按照(二)的规则计算 就是把1和2 3和4压到了同一个数组空间里 叫做压位

这样存储的数字计算次数就变成了2 很厉害的优化……

那我们为什么不这样呢

没错这是10000进制……

计算次数变成了1 其实仔细看一下就知道压位的思想就是把本来不用高精的计算尽量的扩大

来达到减少计算次数的目的

特殊的输出

在最后的输出答案部分 需要一些处理(下面假设写压P位的高精(也就是把P个数压到同一个数组空间里))

一句话概括:除最高位外,一切不足P位的数字输出时需要在前面补上0使得满足P位

比如压2位 2 + 99 为例 最后会用到两个数组空间(一个放1,一个放1…)但是这两个一代表的意义不同,

第二个可以原样输出因为这是最高位(即整个数的长度modP后得到的剩余部分 这一部分是独立于其他部分的 如果补上0 最后就是0101 是不对的),第一个就需要补上一个0来使得输出101而不是11

那么怎么实现呢 下面代码

代码实现

习惯写很多注释 在这里不多做解释

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>#define p 8//要压的位数 
#define carry 100000000//相应的10的P次方 用于进位 
//这样如果要改变压位的位数只要改这里就可以 
using namespace std;const int Maxn=50001;char s1[Maxn],s2[Maxn];
int a[Maxn],b[Maxn],ans[Maxn];int change(char s[],int n[])//压位的核心部分 
{char temp[Maxn];//中间变量 记录每p位的数 int len=strlen(s+1),cur=0;while(len/p){//如果len大于等于p strncpy(temp,s+len-p+1,p);//从后面截出来p位数 n[++cur]=atoi(temp);//把temp搞成数字 len-=p;//继续下p位 }if(len){//如果最后len不是正好p的倍数 也就是还剩下点不够p位的 memset(temp,0,sizeof(temp));strncpy(temp,s+1,len);//全截下来 n[++cur]=atoi(temp);//赋上 }return cur;//返回一个位数 
}
//这里就是(二)中的内容了 计算正常计算就行 
int add(int a[],int b[],int c[],int l1,int l2)
{int x=0,l3=max(l1,l2);for(int i=1;i<=l3;i++){c[i]=a[i]+b[i]+x;x=c[i]/carry;//进位 c[i]%=carry;}while(x>0){c[++l3]=x%10;x/=10;}return l3;//返回答案的位数 
}void print(int a[],int len)
{ printf("%d",a[len]);//处理高位 for(int i=len-1;i>=1;i--)printf("%0*d",p,a[i]);//输出p位 不足补0 printf("\n");
}int main()
{
//	freopen("t.in","r",stdin);
//	freopen("t.out","w",stdout);scanf("%s%s",s1+1,s2+1);//读入两个字符串 int la=change(s1,a);//将s1这个字符串转化为a这个整型数组 int lb=change(s2,b);//同上 int len=add(a,b,ans,la,lb); //计算长度为la的a数组和长度为lb的b数组最后把答案赋给ans数组 并顺便计算出ans的长度(便于输出) print(ans,len);//输出函数 
}
欢迎指出Bug


End。

这篇关于高精度计算(三)压位的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/661669

相关文章

poj 1113 凸包+简单几何计算

题意: 给N个平面上的点,现在要在离点外L米处建城墙,使得城墙把所有点都包含进去且城墙的长度最短。 解析: 韬哥出的某次训练赛上A出的第一道计算几何,算是大水题吧。 用convexhull算法把凸包求出来,然后加加减减就A了。 计算见下图: 好久没玩画图了啊好开心。 代码: #include <iostream>#include <cstdio>#inclu

uva 1342 欧拉定理(计算几何模板)

题意: 给几个点,把这几个点用直线连起来,求这些直线把平面分成了几个。 解析: 欧拉定理: 顶点数 + 面数 - 边数= 2。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#inc

uva 11178 计算集合模板题

题意: 求三角形行三个角三等分点射线交出的内三角形坐标。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <stack>#include <vector>#include <

XTU 1237 计算几何

题面: Magic Triangle Problem Description: Huangriq is a respectful acmer in ACM team of XTU because he brought the best place in regional contest in history of XTU. Huangriq works in a big compa

音视频入门基础:WAV专题(10)——FFmpeg源码中计算WAV音频文件每个packet的pts、dts的实现

一、引言 从文章《音视频入门基础:WAV专题(6)——通过FFprobe显示WAV音频文件每个数据包的信息》中我们可以知道,通过FFprobe命令可以打印WAV音频文件每个packet(也称为数据包或多媒体包)的信息,这些信息包含该packet的pts、dts: 打印出来的“pts”实际是AVPacket结构体中的成员变量pts,是以AVStream->time_base为单位的显

计算数组的斜率,偏移,R2

模拟Excel中的R2的计算。         public bool fnCheckRear_R2(List<double[]> lRear, int iMinRear, int iMaxRear, ref double dR2)         {             bool bResult = true;             int n = 0;             dou

GPU 计算 CMPS224 2021 学习笔记 02

并行类型 (1)任务并行 (2)数据并行 CPU & GPU CPU和GPU拥有相互独立的内存空间,需要在两者之间相互传输数据。 (1)分配GPU内存 (2)将CPU上的数据复制到GPU上 (3)在GPU上对数据进行计算操作 (4)将计算结果从GPU复制到CPU上 (5)释放GPU内存 CUDA内存管理API (1)分配内存 cudaErro

Java - BigDecimal 计算分位(百分位)

日常开发中,如果使用数据库来直接查询一组数据的分位数,就比较简单,直接使用对应的函数就可以了,例如:         PERCENT_RANK() OVER(PARTITION BY 分组列名 ORDER BY 目标列名) AS 目标列名_分位数         如果是需要在代码逻辑部分进行分位数的计算,就需要我们自己写一个工具类来支持计算了 import static ja

OpenStack离线Train版安装系列—2计算节点-环境准备

本系列文章包含从OpenStack离线源制作到完成OpenStack安装的全部过程。 在本系列教程中使用的OpenStack的安装版本为第20个版本Train(简称T版本),2020年5月13日,OpenStack社区发布了第21个版本Ussuri(简称U版本)。 OpenStack部署系列文章 OpenStack Victoria版 安装部署系列教程 OpenStack Ussuri版

新一代车载(E/E)架构下的中央计算载体---HPC软件架构简介

老规矩,分享一段喜欢的文字,避免自己成为高知识低文化的工程师: 屏蔽力是信息过载时代一个人的特殊竞争力,任何消耗你的人和事,多看一眼都是你的不对。非必要不费力证明自己,无利益不试图说服别人,是精神上的节能减排。 无人问津也好,技不如人也罢,你都要试着安静下来,去做自己该做的事.而不是让内心的烦躁、焦虑、毁掉你本就不多的热情和定力。 时间不知不觉中,快要来到夏末秋初。一年又过去了一大半,成