高精度压位模版

最新推荐文章于 2021-11-14 18:22:27 发布
最新推荐文章于 2021-11-14 18:22:27 发布
阅读量481 收藏 2
点赞数 2
分类专栏： 高精度
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_44584560/article/details/89405691
版权
本文深入探讨了高精度算术的实现，包括加法、减法、乘法和除法等核心算法，并提供了详细的C++代码示例。通过具体题目如luogu1601、luogu2142、luogu1303等，文章展示了如何应用这些算法解决实际问题。
摘要由CSDN通过智能技术生成

文章目录

注：昨天的高减重载小于运算符写错了，应该是倒序从高位到低位 ——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

高精度模板（压位版）
qq_37816449的博客
10-29 575
（这次没偷懒，加减乘除模开方都有。。） 一直不想写压位，直到。。。。。。。。。。 遇到了一道会卡不压位高精度的题，心痛到无法呼吸 而且是在考试的时候。。。。。。。。。。 于是临时推，发现压位的输入输出真的是恶心 但是压位的优点也是明显的：速度快，内存低 缺点在于，如果让你保留一定位数的话，不到万不得已别写这玩意儿 以及在高精度除以整型时，如果整型过大，是会炸掉的。 关于
高精度模板——C++实现
xiji333的博客
10-17 1555
文章目录综述朴素想法压位模板 综述 最近闲来无事写了一个比较简单的高精度模板……想来之前没有总结过这一块的知识，就写一篇博客来简单记录一下吧。关于高精度模板的资料其实挺多的，有某个方法不是很理解的话可以再搜搜资料或者自己手动模拟一下。这个模板的效率肯定不是最好的，可能还会有隐藏的bug，仅给大家提供一种实现高精度的思路。 朴素想法 我们先从最简单的开始，如果让你实现一个正数的高精度加法，你会怎么做？大部分人的做法肯定是用数组或者链表来模拟这个过程，以数组为例，从低位到高位依次存储数字的每一位，那么加法无非就
压位高精度
ankong9844的博客
02-23 408
前言 本文的所有代码用C++实现，若在代码段里发现未定义的行为，请到最后的完整代码里寻找该行为含义。若有错误敬请在评论区指正。谢谢您的阅读。 引入 普通高精度每一位存一个数字，可是对于一个int型的存储空间，可存的最大数值为2147483647，只存一个数字是极大的空间浪费。因此，为了在节约空间的同时减少时间复杂度，我们尽量在一个空间中多存几位数字。对于乘法，两数之积的位数最大可达...
压位高精度模板
08-28 247
不走程序，直接上板子。 第一个板子压四位，支持带符号的加减乘。 第二个板子压九位，支持不带符号的四则运算和取模。 都封装了。 #include <iostream> #include <cstdio> #include <cstring> using namespace std; struct intX { static const int ...
【模板】压位高精度
sdfzchy的博客
07-26 4079
【模板】高精度
压位高精度模板，支持负数
04-25
自己写的压位高精度支持负数，保证无bug，单精高精运算时不会像其他人一样偷懒直接类型转换再运算，会很慢，可以+-*/%以及快速幂，压位大小可以自己改(不过输出流重载的地方也要改一下)
C++高精度模板
03-05
支持（无限）高精度加、减、乘、除、模、大于、小于、大于等于、小于等于、位移、等于、不等于、++、--，cin>>，cout 总之，一切int类型的运算都支持
C++高精度压9位整数，带FFT，牛顿迭代除法 模板
最新发布
06-09
压9位，使用压3位FFT乘法和牛顿迭代除法，速度极快 洛谷链接：https://www.luogu.com.cn/article/p6un33cl 乘除法复杂度均为O(n*log(n)) 同时支持gcd/lcm，gcd使用stein算法实现 支持高精度开根号，采用二分 ...
C++ 高精度封装模板
08-27
此模板为已封装模板，包含 高精度的+，-，*，/，不包含小数 by LAZ
高精度类（压位）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即可~非常方便
高精度压位
unknown
07-08 430
高精度压位： http://blog.csdn.net/ssssssay/article/details/52102173
{模板}高精度压位
GH的博客
10-06 755
乘：procedure c1(x:longint); var i,j,k:longint; y:int64; begin y:=0; for i:=1 to ans[0] do begin ans[i]:=ans[i]*x+y; y:=ans[i] div mo; ans[i]:=ans[i] m
【算法专题】高精度之压位
weixin_42638946的博客
11-14 3904
高精度之压位 当数据过大时，此时long long存储不下，因此需要使用vector或者数组存储，然后计算。 一般vector或者数组中每个数据都是一个int，如果每个位置只是存储0~9一位数字的话，比较浪费空间，并且计算也会变慢。因此可以让每个位置存储连续的多位数字，这被称作压位。 这里以压4位为例，也就是说：vector或者数组中每个位置存储四个连续的数字。 1. 基本高精度运算 关于基本高精度计算可以参考：高精度运算 2. 高精度压位运算 压位和不压位的高精度计算存在三点不同点（以下提到的压
"高精度压力、差压变送器行业办公室管理制度及模板2021
最新高精度压力、差压变送器行业办公室管理制度.pdf是2021年发布的一份办公室管理制度，主要针对高精度压力、差压变送器行业，旨在规范办公室工作行为，提高工作效率和企业形象。 该文件首先介绍了制度的背景和目的...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值