高精度压位模版

最新推荐文章于 2023-01-30 22:47:00 发布
最新推荐文章于 2023-01-30 22:47:00 发布
阅读量459 收藏 2
点赞数 2
分类专栏： 高精度
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_44584560/article/details/89405691
版权

文章目录

注：昨天的高减重载小于运算符写错了，应该是倒序从高位到低位 ——2019.4.20


bool operator <(const num &p,const num &q){//比较大小 
	if(p.a[0]<q.a[0]) return true;
	if(p.a[0]>q.a[0]) return false;
	for(int i=p.a[0];i>=1;--i){//倒序 
		if(p.a[i]!=q.a[i]) return p.a[i]<q.a[i];
	}
	return false;
}

高精度加法

/***************************************
User:Mandy.H.Y
Language:c++
Problem:luogu1601 A+B Problem（压位）
Algorhithm:High Precision Addition
****************************************/

#include<bits/stdc++.h>

using namespace std;

const int power=4; //压位压四位 
const int base=1e4;
const int maxn=505;

struct num{
	int a[maxn];
	
	num(){memset(a,0,sizeof(a));}
	
	int &operator [](int x){return a[x];}//重载中括号 
	
	num(char *s){
		
		memset(a,0,sizeof(a));
		int len=strlen(s);
		
		a[0]=(len+power-1)/power;//计算位数
		
		for(int i=0,t=0,w=0;i<len;++i,w*=10){
			if(i%power==0) {w=1,++t;}//最好写成这种形式 
			a[t]+=(s[i]-'0')*w;
			
		}
	}
	
	void print(){
		printf("%d",a[a[0]]);
		for(int i=a[0]-1;i>0;--i) printf("%04d",a[i]);//printf("%0*d",power,a[i]); 
	}
}ans,p,q;

num operator +(num &p,num &q){//重载了中括号所以不用const 
	int cur=0,len=max(p[0],q[0]);
	for(int i=1;i<=len;++i){
		cur+=i<=p[0]?p[i]:0;
		cur+=i<=q[0]?q[i]:0;
		p[i]=cur%base;cur/=base;
	}
	if(cur) p[++len]=cur;
	p[0]=len;//注意更新长度 
	return p;
}

char p1[maxn],q1[maxn];
int lenp,lenq;

int main(){
	scanf("%s%s",p1,q1);
	
	lenp=strlen(p1);lenq=strlen(q1);
	reverse(p1,p1+lenp);reverse(q1,q1+lenq);
	
	p=num(p1);q=num(q1);
	
	ans=p+q;
	ans.print();
	return 0;
}

高精度减法

/**************************************
User:Mandy.H.Y
Language:c++
Pronblem:luogu2142 A-B Problem
Algorhithm:High Precision Subtraction 
**************************************/

#include<bits/stdc++.h>

using namespace std;

const int power=4;
const int base=1e4;
const int maxn=1e4+5;

struct num{
	int a[maxn];
	
	num(){memset(a,0,sizeof(a));}
	
	num(char *s,int len){
		memset(a,0,sizeof(a));
		
		a[0]=(len+power-1)/power;
		
		for(int t=0,w=1,i=len-1;i>=0;w=(w<<1)+(w<<3),--i){//倒序省去reverse
			if((len-1-i)%power==0){++t,w=1;}
			a[t]+=(s[i]^48)*w;
		}
	}
	
	void print(){
		printf("%d",a[a[0]]);
		for(int i=a[0]-1;i>0;--i) printf("%0*d",power,a[i]);
	}
}p,q,ans;

bool operator <(const num &p,const num &q){//比较大小 
	if(p.a[0]<q.a[0]) return true;
	if(p.a[0]>q.a[0]) return false;
	for(int i=p.a[0];i>=1;--i){//倒序 
		if(p.a[i]!=q.a[i]) return p.a[i]<q.a[i];
	}
	return false;
}

num operator -(num &p,num &q){
	int len=max(p.a[0],q.a[0]);
	for(int i=1;i<=len;++i){
		p.a[i]-=i<=q.a[0]?q.a[i]:0;
		if(p.a[i]<0) p.a[i]+=base,p.a[i+1]-=1;//借位 
	}
	while(p.a[0]>1&&p.a[p.a[0]]==0) --p.a[0];
	return p;
}

char p1[maxn],q1[maxn];
int lenp,lenq;

int main(){
	scanf("%s%s",p1,q1);
	lenp=strlen(p1);lenq=strlen(q1);
	p=num(p1,lenp),q=num(q1,lenq);
	if(q<p) ans=p-q;
	else putchar('-'),ans=q-p;
	ans.print();
	return 0;
}

高精度乘法

/*********************************************

User:Mandy.H.Y
Language:c++
Problem:luogu 1303 A*B Problem
Algorithm: High Precision Multiplication

**********************************************/

#include<bits/stdc++.h>

using namespace std;

const int power=4;
const int base=1e4;
const int maxn=2e3+5;

struct num{
	int a[maxn<<1];
	
	num(){memset(a,0,sizeof(a));}
	
	num(char *s,int len){
		memset(a,0,sizeof(a));
		
		a[0]=(len+power-1)/power;
		
		for(int t=0,w=1,i=len-1;i>=0;w=(w<<1)+(w<<3),--i){
			if((len-1-i)%power==0) w=1,++t;
			a[t]+=(s[i]^48)*w;//注意+= 
		}
	}
	
	void print(){
		printf("%d",a[a[0]]);
		for(int i=a[0]-1;i>0;--i) printf("%0*d",power,a[i]);
	}
}ans,p,q;

char p1[maxn],q1[maxn];
int lenp,lenq;

num operator *(const num &p,const num &q){
	num b;
	for(int i=1;i<=p.a[0];++i){
		int cur=0;
		for(int j=1;j<=q.a[0];++j)
			cur+=b.a[i+j-1]+p.a[i]*q.a[j],b.a[i+j-1]=cur%base,cur/=base;
			//cur+=b.a[i+j-1]+p.a[i]*q.a[j]累加 b.a[i+j-1]
		b.a[i+q.a[0]]=cur;
	}
	b.a[0]=q.a[0]+p.a[0];//得到大致位数 
	while(b.a[0]>1&&b.a[b.a[0]]==0) --b.a[0];//除去前导0 
	return b;
}

int main(){
	scanf("%s%s",p1,q1);
	lenp=strlen(p1);lenq=strlen(q1);
	p=num(p1,lenp);q=num(q1,lenq);
	ans=p*q;
	ans.print();
	return 0;
}

整合

/**************************
User:Mandy.H.Y
Language:c++
Problem:
Algorithm: High Precision
***************************/

#include<bits/stdc++.h>

using namespace std;

const int power=4;
const int base=1e4;
const int maxn=2e3+5;

struct num{
	int a[maxn<<1];
	
	num(){memset(a,0,sizeof(a));}
	int &operator [](int x){return a[x];}//重载中括号 
	num(char *s,int len){
		memset(a,0,sizeof(a));
		
		a[0]=(len+power-1)/power;
		
		for(int t=0,w=1,i=len-1;i>=0;w=(w<<1)+(w<<3),--i){
			if((len-1-i)%power==0) w=1,++t;
			a[t]+=(s[i]^48)*w;//注意+= 
		}
	}
	void add(int k) { if (k || a[0]) a[ ++a[0] ] = k; }     //在末尾添加一个数，除法的时候要用到  
    void re() { reverse(a+1, a+a[0]+1); } 
	void print(){
		printf("%d",a[a[0]]);
		for(int i=a[0]-1;i>0;--i) printf("%0*d",power,a[i]);
	}
}ans,p,q;

char p1[maxn],q1[maxn];
int lenp,lenq;


num operator +(num &p,num &q){//重载了中括号所以不用const 
	int cur=0,len=max(p[0],q[0]);
	for(int i=1;i<=len;++i){
		cur+=i<=p[0]?p[i]:0;
		cur+=i<=q[0]?q[i]:0;
		p[i]=cur%base;cur/=base;
	}
	if(cur) p[++len]=cur;
	p[0]=len;//注意更新长度 
	return p;
}

bool operator <(const num &p,const num &q){//比较大小 
	if(p.a[0]<q.a[0]) return true;
	if(p.a[0]>q.a[0]) return false;
	for(int i=p.a[0];i>=1;--i){//倒序 
		if(p.a[i]!=q.a[i]) return p.a[i]<q.a[i];
	}
	return false;
}

num operator -(num &p,num &q){
	int len=max(p.a[0],q.a[0]);
	for(int i=1;i<=len;++i){
		p.a[i]-=i<=q.a[0]?q.a[i]:0;
		if(p.a[i]<0) p.a[i]+=base,p.a[i+1]-=1;//借位 
	}
	while(p.a[0]>1&&p.a[p.a[0]]==0) --p.a[0];
	return p;
}

num operator / (num &p,num &q)               //注意const不要冲突 
{  
    num x, y;  
    for (int i = p.a[0];i >= 1;--i)                       //从最高位开始取数  
    {  
        y.add(p.a[i]);             //把数添到末尾（最低位），这时候是高位在前，低位在后  
        y.re();                    //把数反过来，变为统一的存储方式：低位在前，高位在后  
        while ( !(y < q) )         //大于等于除数的时候，如果小于的话，其实答案上的该位就是初始的“0”  
            y = y - q, ++x.a[i];   //看能减几个除数，减几次，答案上该位就加几次。  
        y.re();                    //将数反过来，为下一次添数做准备  
    }  
    x.a[0] = p.a[0];  
    while (x.a[0] > 0 && !x.a[x.a[0]]) --x.a[0];  
    return x;  
}  

num operator *(const num &p,const num &q){
	num b;
	for(int i=1;i<=p.a[0];++i){
		int cur=0;
		for(int j=1;j<=q.a[0];++j)
			cur+=b.a[i+j-1]+p.a[i]*q.a[j],b.a[i+j-1]=cur%base,cur/=base;
		b.a[i+q.a[0]]=cur;
	}
	b.a[0]=q.a[0]+p.a[0];//得到大致位数 
	while(b.a[0]>1&&b.a[b.a[0]]==0) --b.a[0];//除去前导0 
	return b;
}

int main(){
	
	/*
	scanf("%s%s",p1,q1);
	lenp=strlen(p1);lenq=strlen(q1);
	p=num(p1,lenp);q=num(q1,lenq);
	ans=p*q;
	ans.print();
	return 0;
	*/
	
	
	scanf("%s%s",p1,q1);
	
	lenp=strlen(p1);lenq=strlen(q1);
	
	p=num(p1,lenp);q=num(q1,lenq);
	
	ans=p+q;
	ans.print();
	return 0;
	
	
	
	/*
	
	scanf("%s%s",p1,q1);
	lenp=strlen(p1);lenq=strlen(q1);
	p=num(p1,lenp),q=num(q1,lenq);
	if(q<p) ans=p-q;
	else putchar('-'),ans=q-p;
	ans.print();
	return 0;

	*/
	
}

练习

luogu1255

luogu1255

/************************
User:Mandy.H.Y
Language:c++
Problem:luogu 1255 Staircase
Algorithm:High Precision Addition
************************/

#include<bits/stdc++.h>

using namespace std;

const int power=4;
const int base=1e4;
const int maxn=5005;

struct num{
	int a[maxn];
	num(){memset(a,0,sizeof(a));}
	int &operator [](int x){return a[x];}
	void print(){
		printf("%d",a[a[0]]);
		for(int i=a[0]-1;i>0;--i) printf("%04d",a[i]);
	}
}ans[5];

int n;

num operator +(num &p,num &q){
	int cur=0,len=max(p[0],q[0]);
	for(int i=1;i<=len;++i){
		cur+=i<=p[0]?p[i]:0;
		cur+=i<=q[0]?q[i]:0;
		p[i]=cur%base,cur/=base;
	}
	p[0]=len;//注意继承长度 
	if(cur) p[++p[0]]=cur;
	return p;
}

int main(){
	scanf("%d",&n);
	if(n==1) {puts("1");return 0;}
	if(n==2) {puts("2");return 0;}
	ans[1].a[0]=1,ans[2].a[0]=1;ans[1].a[1]=1,ans[2].a[1]=2;
	for(int i=3;i<=n;++i){
		ans[i%3]=ans[(i+1)%3]+ans[(i+2)%3];
	}
	ans[n%3].print();
	return 0;
}

luogu1604

luogu1604

/***************************
User:Mandy.H.Y
Language:c++
Problem:luogu1604 Planet
Algorithm:High Precision
***************************/

#include<bits/stdc++.h>

using namespace std;

const int power=1;
const int maxn=2e4+5;

int n,base;

struct num{
	int a[maxn];
	num(){memset(a,0,sizeof(a));}
	int &operator [](int x){return a[x];}
	num(char *s,int len){
		a[0]=len;
		for(int i=len-1;i>=0;--i){
			if(isdigit(s[i])) a[len-i]=s[i]^48;
			else a[len-i]=s[i]-'A'+10;
		}
	}
	void print(){
		for(int i=a[0];i>0;--i){
			if(a[i]<10) printf("%d",a[i]);
			else printf("%c",a[i]-10+'A');
		}
	}
}p,q,ans;

char p1[maxn],q1[maxn];
int lenp,lenq;

num operator +(num &p,num &q){
	int cur=0,len=max(p[0],q[0]);
	for(int i=1;i<=len;++i){
		cur+=i<=p[0]?p[i]:0;
		cur+=i<=q[0]?q[i]:0;
		p[i]=cur%base,cur/=base;
	}
	if(cur) p[++len]=cur;
	p[0]=len;
	return p;
}

int main(){
	scanf("%d",&base);
	scanf("%s%s",p1,q1);
	lenp=strlen(p1);lenq=strlen(q1);
	p=num(p1,lenp);q=num(q1,lenq);
	ans=p+q;
	ans.print();
	return 0;
}

部分引自：https://www.cnblogs.com/hnqw1214/p/6351321.html

长剑凌清秋
关注
  • 2
    点赞
  • 踩
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
高精度算法-压位
LiuKairui的博客
07-31 5729
我们之前做过大整数类的运算的题目 大整数乘法 大整数加法 这个方法看似是无敌的，，， 但是那么如果是一个10000^10000位的数据呢？ 数组根本开不到这么大的。。。 有这样的题目吗？ 传送门 这时候我们就需要压位了。。。 还记得我们存储数字的方式？ 举个栗子:1234 + 1234这是之前的存储方式 实际上我们一个数组空间(int)是可以储存 2147483647 以...
高精度4位压缩法原理与实现
挽颜__
01-11 2409
高精度原理与实现
【算法专题】高精度之压位
weixin_42638946的博客
11-14 3314
高精度之压位 当数据过大时，此时long long存储不下，因此需要使用vector或者数组存储，然后计算。 一般vector或者数组中每个数据都是一个int，如果每个位置只是存储0~9一位数字的话，比较浪费空间，并且计算也会变慢。因此可以让每个位置存储连续的多位数字，这被称作压位。 这里以压4位为例，也就是说：vector或者数组中每个位置存储四个连续的数字。 1. 基本高精度运算 关于基本高精度计算可以参考：高精度运算 2. 高精度压位运算 压位和不压位的高精度计算存在三点不同点（以下提到的压
高精度计算(三)压位
sSay
08-03 1万+
压位高精
压位高精度
ankong9844的博客
02-23 373
前言 本文的所有代码用C++实现，若在代码段里发现未定义的行为，请到最后的完整代码里寻找该行为含义。若有错误敬请在评论区指正。谢谢您的阅读。 引入 普通高精度每一位存一个数字，可是对于一个int型的存储空间，可存的最大数值为2147483647，只存一个数字是极大的空间浪费。因此，为了在节约空间的同时减少时间复杂度，我们尽量在一个空间中多存几位数字。对于乘法，两数之积的位数最大可达...
高精度压位乘法
chang_yl的博客
02-24 2019
题面：地图(compress.cpp / stdin / stdout)stdin / stdoutTL: 1 s / ML: 512 MiB题目描述众而周知，8000m×8000m       的地图相当大， 却还是满足不了萌萌哒选手 compress。compress 说：“我要 xm×ym 的大地图！” 既然你无法满足 compress 的地图梦， 那你就帮他算算他要的地图到底有多大吧！输入...
压位高精度模板，支持负数
04-25
自己写的压位高精度支持负数，保证无bug，单精高精运算时不会像其他人一样偷懒直接类型转换再运算，会很慢，可以+-*/%以及快速幂，压位大小可以自己改(不过输出流重载的地方也要改一下)
C++ 高精度封装模板
08-27
此模板为已封装模板，包含 高精度的+，-，*，/，不包含小数 by LAZ
高精度模板
11-16
高精度模板，整合了多篇文档，有高精度加法，高精度减法，高精度乘法，高精度除法，高精度取模法等模板...
C++高精度压9位整数，带FFT，牛顿迭代除法 模板
最新发布
06-09
压9位，使用压3位FFT乘法和牛顿迭代除法，速度极快 洛谷链接：https://www.luogu.com.cn/article/p6un33cl 乘除法复杂度均为O(n*log(n)) 同时支持gcd/lcm，gcd使用stein算法实现 支持高精度开根号，采用二分 ...
高精度类（压位）C++实现
06-03
高精度压位 #include #include #include using namespace std; const int maxleng=500; class BigInt { private: int leng; int num[maxleng]; public:
C++ （压位）高精度四则运算
07-23
写的比较漂亮的一个高精度四则运算 压位只需改两个const即可~非常方便
2021年高精度压力、差压变送器行业财务部门表格模板汇总 .pdf
09-28
2021年各行业财务部门表格模板
C/C++高精度运算(大整数运算)详解（含压位）
热门推荐
帅比王的博客
04-22 2万+
1.高精度加法 1.1 高精度加法 高精度运算的基本运算就是加和减。和算数的加减规则一样，模拟竖式计算，考虑错位运算与进位处理。下面是我老师给的代码，目前比网上其他的代码要精简和巧妙。 #include <cstdio> #include <cstring> int main() { char a[202]={0}, b[202]={0...
C/C++ 高精度（加减乘除）算法压位优化
codeofcc的博客
01-30 957
由于上一章《C/C++ 高精度（加减乘除）算法简单实现》实现了基本的高精度计算，数组的每个元素存储一位10进制的数字。这样的存储方式并不是最优的，32位的整型其实至少可以存储9位高精度数字，数组元素存储更多的位数就是压位优化。本文将展示压位优化的原理以及压9位的实现和示例。以上就是今天要讲的内容，压位优化性能提升是比较显著的，而且实现也很容易，大部分逻辑是一致的只是底数变大了而已。从性能测试结果来看所有计算至少由4倍的提升，乘法性能提升较大有可能是测试方法不严重，这个待以后验证。
高精度计算（+-*/）（压位计算）（c++实现）
weixin_60257971的博客
11-24 939
进来的时候如果没看见除法，说明我也还不会，会了的话我会补充上去的。如果要看的除法的，得去看别人的了。 目录 初始化： 加法： 减法： 乘法： 输出：这个是重点！！！！！！！！！！！ 完整代码： 如果你是从我的“P1601 A+B Problem（高精）”这篇文章来的话，可以看一下压位计算，更加的快捷方便。 高精度计算就是因为超过了数据表达的类型，从用其他的方式来进行计算。 压位就是说将一长串的数字分割开来，来单独进行保存，从而不会发生超过范围。 初始化： 可以看到 k 其...
（压位）高精度思想教程 与 代码实现
IcEternal
07-23 5505
高精度这玩意儿真是好久没写了，时隔多年，从 Pascal 转到 C++ , 代码风格也变了 重新编了一下（经过潘神指导！），觉得比以前写的好多了~ 想了想干脆写一个小教程，希望初学者能有所收获~~ 一.高精度四则运算思想 高精度这个东西其实谁都会，上过小学数学课就知道高精度怎么做了，就是按照人的运算方式，一位一位运算。 1.加法  就是 从个位开始 两个数字相加，如
高精度-四位压缩（结构体重载运算符）+DP（洛谷P1005-方格取数）
m0_55758733的博客
03-08 879
洛谷P1005方格取数 这是一道典型的dp题 且根据题意得行与行之间没有任何的联系，只需找到每一行的最优解相加就得出最终答案。 设dp数组为f[maxn][maxn],并将每一行的数据储存在数组a[maxn]中， 分析题意dp规则为: f[q][w]=max(f[q-1][w]+pow(2,m-(w-q+1))*a[q-1],f[q][w+1]+pow(2,m-(w-q+1))a[w+1]); 但是提交后会发现只能通过60%的数据原因就下面： 其中1<=n,m<=80，也就意味着要计算2的80
高精度与高精度间的模板
07-27
引用\[2\]提供了一个适用于高精度与高精度相加的模板。这个模板使用了字符串来存储大数，并通过逆序存储的方式将字符串转化为数字。然后使用两层循环来模拟乘法运算的过程，将结果保存在数组c中。最后，通过删除前导...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值