第一章基础算法（二）高精度、前缀和、差分

最新推荐文章于 2024-09-12 19:02:20 发布

wulalalalalalalal

最新推荐文章于 2024-09-12 19:02:20 发布

阅读量570

点赞数 2

分类专栏：算法文章标签：算法 c++ 开发语言

本文链接：https://blog.csdn.net/wangtongyao_/article/details/122268633

版权

算法专栏收录该内容

24 篇文章 0 订阅

订阅专栏

这篇博客详细介绍了如何在C++中实现高精度计算，包括加法、减法、乘法和除法，并提供了相应的模板代码。同时，讲解了前缀和的概念及其在求解数组区间和中的应用，包括一维和二维数组的前缀和计算。这些算法在处理大整数运算和数组求和问题时非常实用。

摘要由CSDN通过智能技术生成

1. 高精度

1. 三类高精度问题：
（1） 大整数A + 大整数B；
（2）大整数A -  大整数B；
（3）大整数A * 小正数a；
（4）A/b
2. 大整数的位数：小于 10 ^ 6；小整数的大小：小于 10000。
3. 大整数如何存储：存储在数组中，数组的每一位存一位数字，数组的第0位存大整数的个位

期中用到的c++知识：

vector
在 c++ 中，vector 是一个十分有用的容器。它能够像容器一样存放各种类型的对象，简单地说，vector是一个能够存放任意类型的动态数组，能够增加和压缩数据。vector 是同一种类型的对象的集合，每个对象都有一个对应的整数索引值。
向量（Vector）是一个封装了动态大小数组的顺序容器（Sequence Container）。
（1）顺序序列：顺序容器中的元素按照严格的线性顺序排序。可以通过元素在序列中的位置访问对应的元素。
（2）动态数组：支持对序列中的任意元素进行快速直接访问。操供了在序列末尾相对快速地添加/删除元素的操作。

1.vector<T> v1;     保存类型为 T 对象。默认构造函数 v1 为空。
2.void push_back(const T& x):向量尾部增加一个元素X
3.void pop_back():删除向量中最后一个元素
4.back 得到数组的最后一个单元的引用

2.auto 类型推导，编译器会自己推断这个变量是什么类型的。
3.reverse函数：逆序，反转在[first,last)范围内的顺序。使用时需要包含头文件，#include <algorithm>。

一、高精度加法模板：大整数 + 大整数：

#include<iostream>
#include<vector>

using namespace std;

vector<int> add(vector<int> &A,vector<int> &B) 
// 加上引用是为了提高效率，不加引用会把整个数组拷贝一份
{
	// C=A+B
	vector<int> C;
	int t = 0; // 进位
	 
	for(int i = 0;i < A.size() || i < B.size();i++)
	{
		if(i < A.size()) t += A[i]; // A、B的位数可能不一样
		if(i < B.size()) t += B[i];
		C.push_back(t % 10); // 当前位 
		t = t / 10;
	}
	if(t) // 如果向最高位还有进位
		C.push_back(t);
	return C; 
 } 

int main()
{
	string a,b; // 将大数用字符串读进来 
	vector<int> A,B,C;
	
	cin >> a >> b; // a = "123456"
	for(int i = a.size() - 1;i >= 0;i--)
		A.push_back(a[i]-'0'); // 将字符型变为整型
	// 此时	，A = [6,5,4,3,2,1] 
	for(int i = b.size() - 1;i >= 0;i--)
		B.push_back(b[i]-'0');	 
 
	// auto C = add(A,B); // DEVC++报错
	C = add(A,B); 
	// 输出
	for(int i = C.size() - 1;i >= 0;i--)
		printf("%d",C[i]);
	return 0;
}

二、高精度减法模板：大整数 - 大整数
① A >= B：A - B = A - B；② A < B：A - B = - ( B - A )。
模板中没有考虑有负数的情况，假定A和B都是正数。（如果A，B两个数有负数，一定可以转成两个数的绝对值相加或相减）

#include<iostream>
#include<vector>
using namespace std;

// 判断A是否大于等于B 
bool cmp(vector<int> &A,vector<int> &B)
{
	if(A.size() != B.size()){
		return A.size() > B.size();
	} 
	for(int i = A.size()-1; i >= 0;i--)
	{
		if(A[i] != B[i]){
			return A[i] > B[i];
		} 
		// 找到第一位不相等的位比较大小 
	}
	return true; // A == B，也返回 true 
} 

vector<int> sub(vector<int> &A,vector<int> &B)
{
	vector<int> C;
	 // 因为前面已经判断了A，B的大小，所以这里A一定大于B 	 
	for(int i = 0,t = 0;i < A.size();i++)
	{
		t = A[i] - t;
		if(i < B.size()) t -= B[i]; // 判断B有没有第i位 
		C.push_back((t + 10) % 10);
		if(t < 0) t = 1;
		else t = 0;
	}
	while(C.size() > 1 && C.back() == 0) C.pop_back();
	//  防止出现123-120 = 003 情况，去掉前导0
	return C; 
}

int main()
{
	string a,b; // 将大数用字符串读进来 
	vector<int> A,B,C;
		
	cin >> a >> b; // a = "123456"
	
	for(int i = a.size() - 1;i >= 0;i--)
		A.push_back(a[i]-'0'); // 将字符型变为整型
	// 此时	，A = [6,5,4,3,2,1] 
	for(int i = b.size() - 1;i >= 0;i--)
		B.push_back(b[i]-'0');	 
	
	if(cmp(A,B))
	{
		C = sub(A,B);
		// 输出
		for(int i = C.size() - 1;i >= 0;i--)
			printf("%d",C[i]);
	}
	else
	{
	 	C = sub(B,A);
	 	printf("-");
	 	for(int i = C.size() - 1;i >= 0;i--)
			printf("%d",C[i]);
	 } 
	 	
	return 0;
}

三、高精度乘法模板：大整数 * 小整数
计算方法，例：123 * 12 = C₃C₂C₁C₀
C₀ = (3 * 12) % 10 = 6
t₁ = (3 * 12) / 10 = 3
C₁ = (2 * 12 + t₁) % 10 = 7
t₂ = (2 * 12 + t₁) / 10 = 2
C₂ = (1 * 12 + t₂) % 10 = 4
t₃ = (1 * 12 + t₂) / 10 = 1
C₃ = t₃ = 1
123 * 12 = C₃C₂C₁C₀ = 1476

#include<iostream>
#include<vector>
using namespace std;

vector<int> mul(vector<int> &A,int b) // 大整数 * 小整数
{
	vector<int> C;
	int t = 0;
	
	for(int i = 0;i < A.size();i++)
	{
		t += A[i] * b;
		C.push_back(t % 10);
		t /= 10;
	}
	if(t > 0) C.push_back(t);	
	while(C.size() > 1 && C.back() == 0) C.pop_back();
	// 去掉前导0
	return C;
}

int main()
{
	string a; // 将大数用字符串读进来 
	int b;
	
	vector<int> A,C;
		
	cin >> a >> b; // a = "123456"
	
	for(int i = a.size() - 1;i >= 0;i--)
		A.push_back(a[i]-'0'); // 将字符型变为整型
	// 此时	，A = [6,5,4,3,2,1] 

	C = mul(A,b);
	
	for(int i = C.size() - 1;i >= 0;i--)
		printf("%d",C[i]);	 
			
	return 0;
}

三、高精度除法模板：大整数 / 小整数
计算方法，例：1234 * 11 = C₃C₂C₁C₀
C₃ = 1 / 11
t₁ = 1 % 11
C₂ = ( t₁ * 10 + 2) / 11
t₂ = ( t₁ * 10 + 2) % 11
C₁ = ( t₂ * 10 + 3) / 11
t₃ = ( t₂ * 10 + 3) % 11
C₀ = ( t₃ * 10 + 4) / 11
t₃ = ( t₃ * 10 + 4) % 11

#include<iostream>
#include<vector>
#include<algorithm>
using namespace std;

// A / b，商C余r 
vector<int> div(vector<int> &A,int b,int &r) // r是余数 
{
	vector<int> C;
	// 除法从最高位开始算 ，加减乘都是从最低位
	 
	for(int i = A.size() - 1;i >= 0;i--)
	{
		r = r * 10 + A[i];
		C.push_back(r / b);
		r %= b;			
	}
	
	// C第0位存的最高位 
	reverse(C.begin(),C.end());
	
	while(C.size() > 1 && C.back() == 0) C.pop_back();
	// 去掉前导0
	
	return C;
}

int main()
{
	string a; // 将大数用字符串读进来 
	int b,r = 0;
	
	vector<int> A,C;
		
	cin >> a >> b; // a = "123456"
	
	for(int i = a.size() - 1;i >= 0;i--)
		A.push_back(a[i]-'0'); // 将字符型变为整型
	// 此时	，A = [6,5,4,3,2,1] 

	C = div(A,b,r);
	
	for(int i = C.size() - 1;i >= 0;i--)
		printf("%d",C[i]);
	cout<<endl<<r<<endl; 
		
	return 0;
}

2. 前缀和

长度为n的数组中存放a₁,a₂,a₃ …… a_n，求前缀和 S_i = a₁ + a₂ + a₃ + … + a_i， S_i 表示数组前 i 项元素和。
（1）如何求 S_i；

	for(i = 1;i <= i;i++)
		S[i] = S[i-1] + a[i]

（2）求S_i 的作用：能快速求出原数组中一段数的和。
求[l,r]区间元素和：S _r - S _l-1
数组a下标一定要从1开始，S₀ = 0,这样比如求[1,10] 的元素和，就可以是S ₁₀ - S ₁

模板1

输入一个长度为 n 的整数序列。
接下来再输入 m 个询问，每个询问输入一对 l,r。
对于每个询问，输出原序列中从第 l 个数到第 r 个数的和。

#include<iostream>
using namespace std;

const int N = 100010;
int n,m;
int a[N],s[N];

int main()
{

	scanf("%d%d",&n,&m);
	for(int i = 1;i <= n;i++)
		scanf("%d",&a[i]);
	for(int i = 1;i <= n;i++)
		s[i] = s[i-1] + a[i]; // 前缀和初始化
	while(m--)
	{
		int l,r;
		scanf("%d%d",&l,&r); // 区间和计算
		printf("%d\n",s[r] - s[l-1]);
	}
	return 0;
}

模板2 子矩阵的和

求黄色阴影部分的元素和。
在这里插入图片描述
S_{（x2，y2）} = S_[x2,y2] - S_[x1,y2] - S_[x2,y1] + S_[x1,y1]
S_[i,j] = S_[i-1,j] + S_[i,j-1] - S_[i-1,j-1] + a_[i,j]

#include<iostream>
using namespace std;

const int N = 1010;
int n,m,q;
int a[N][N],s[N][N];

int main()
{
	scanf("%d%d%d",&n,&m,&q);
	for(int i = 1;i <= n;i++)
	{
		for(int j = 1;j <= m;j++)
		{
			scanf("%d",&a[i][j]);
		}
	}
	for(int i = 1;i <= n;i++)
	{
		for(int j = 1;j <= m;j++)
		{
			s[i][j] = s[i-1][j] + s[i][j-1] - s[i-1][j-1] + a[i][j];
			// 求前缀和
		}
	}
	while(q--)
	{
		int x1,x2,y1,y2;
		scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
		printf("%d\n",s[x2][y2] - s[x1-1][y2] - s[x2][y1-1] + s[x1-1][y1-1]);
		// 算子矩阵的和
	}
	return 0;
}

在这里插入图片描述

模板3 一维差分

差分：前缀和的逆运算，已知 a₁,a₂,a₃ …… a_n，构造 b₁,b₂,b₃ …… b_n，使得：
a_i = b₁ + b₂ + b₃ + …… + b_i
即 b₁ = a₁ ，b₂ = a₂ - a₁，b₃ = a₃ - a₂ … b_n = a_n - a_n-1
a数组称为b数组的前缀和，b数组称为a数组的差分（a — 前缀和数组，b — 原数组）
应用：想让a数组某个连续区间 [l,r] 中的数都加上c，我们只需要在b数组中修改两个数的值 b[l] + c 和 b[r+1] - c ，时间复杂度就由O(n)变为O(1)了。
解法：假定a数组初始的时候都是0

例题：
	输入一个长度为 n 的整数序列。
	接下来输入 m 个操作，每个操作包含三个整数 l,r,c，
	表示将序列中 [l,r] 之间的每个数加上 c。
	请你输出进行完所有操作后的序列。
输入格式：
	第一行包含两个整数 n 和 m。
	第二行包含 n 个整数，表示整数序列。
	接下来 m 行，每行包含三个整数 l，r，c，表示一个操作。
输出格式：
	共一行，包含 n 个整数，表示最终序列。

#include<iostream>
using namespace std;

const int N = 100010;
int n,m;
int a[N],b[N]; // a 是原数组，b 是差分数组
 
void insert(int l,int r,int c)
{
	b[l] += c;
	b[r + 1] -= c;
	/*
		比如现在读入的a数组是：1,3,6,10
		那么差分数组b应该是：1,2,3,4
		由前缀和数组得到差分数组是加上是一个插入的过程
		如何得到差分数组： 
		insert(1,1,1): b[1] = 1; b[2] = -1;  // 得到 b[i] = 1 
		insert(2,2,3): b[2] = b[2] + 3 = 2; b[3] = -3;  // 得到 b[2] = 2 
		insert(3,3,6): b[3] = -3 + 6 = 3; b[3] = -6;  // 得到 b[3] = 3
		insert(4,4,10): b[4] = -6 + 10 = 4; b[5] = -10;  // 得到 b[4] = 4 
		
	*/
}

int main()
{
	scanf("%d%d",&n,&m);
	for(int i = 1;i <= n;i++)
		scanf("%d",&a[i]);
		
	for(int i = 1;i <= n;i++) 
		insert(i,i,a[i]);
		// 由原数组求出了差分数组
	while(m--)
	{
		int l,r,c;
		scanf("%d%d%d",&l,&r,&c);
		insert(l,r,c);
	 } 
	 
	 // 所有的操作完后，再由差分数组求一边原数组
	for(int i = 1; i <= n;i++)
	  	b[i] += b[i-1];
	  	// 将b数组变为自己的前缀和 
	for(int i = 1;i <= n;i++)
		printf("%d ",b[i]);

	return 0;
}

模板4 二维差分矩阵

原矩阵 a，差分矩阵 b，给原矩阵一个子矩阵加上一个值
（x1,y1），(x2,y2）

在这里插入图片描述

#include<iostream>
using namespace std;

const int N = 1010;
int n,m,q;
int a[N][N],b[N][N]; // a 是原数组，b 是差分数组
 
void insert(int x1,int y1,int x2,int y2,int c)
{
	b[x1][y1] += c;
	b[x2+1][y1] -= c;
	b[x1][y2+1] -= c;
	b[x2+1][y2+1] += c;
}

int main()
{
	scanf("%d%d%d",&n,&m,&q);
	for(int i = 1;i <= n;i++)
	{
		for(int j = 1;j <= m;j++)
			scanf("%d",&a[i][j]);
	}
	
		for(int i = 1;i <= n;i++)
	{
		for(int j = 1;j <= m;j++)
			insert(i,j,i,j,a[i][j]);
	}

		// 由原数组求出了差分数组
	while(q--)
	{
		int x1,y1,x2,y2,c;
		scanf("%d%d%d%d%d",&x1,&y1,&x2,&y2,&c);
		insert(x1,y1,x2,y2,c);
	 } 	 
	 // 所有的操作完后，再由差分数组求一边原数组
	 	for(int i = 1;i <= n;i++)
	{
		for(int j = 1;j <= m;j++)
		{
			b[i][j] = b[i-1][j] + b[i][j-1] - b[i-1][j-1] + b[i][j];
		}
	}
	  	// 将b数组变为自己的前缀和 
		for(int i = 1;i <= n;i++)
	{
		for(int j = 1;j <= m;j++)
			printf("%d ",b[i][j]);
		printf("\n");
	}
	return 0;
}