算法基础(2) | 高精度、前缀和、差分

已于 2022-07-03 16:32:11 修改
阅读量204 收藏
点赞数
分类专栏： # AcWing 算法基础课 文章标签： 算法 高精度 前缀和 差分
于 2021-07-14 10:52:14 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_46264660/article/details/118722662
版权

文章目录

一、高精度

1.1 高精度加法

在这里插入图片描述
在这里插入图片描述

#include <iostream>
#include <vector>

using namespace std;

const int N = 1e6 + 10;

vector<int> add(vector<int> &A, vector<int> &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];
        if (i < B.size()) t += B[i];
        C.push_back(t % 10);
        t /= 10;
    }
    
    if (t) C.push_back(1);
    
    return C;
}

int main()
{
    string a, b; //将长度太长的数字 以字符串读入
    vector<int> A, B;//大整数A 和 B 以数组形式存在
    
    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);
    
    for (int i = C.size() - 1; i >= 0; i -- ) printf("%d",C[i]);    
    
    return 0;
}

vector<int> add(vector<int> &A, vector<int> &B)
{
    vector<int> C;
    
    int t = 0;//进位
        
    if (A.size() < B.size()) return add(B, A);
    
    for (int i = 0; i < A.size(); i ++ )
    {
        t += A[i];
        if (i < B.size()) t += B[i];
        C.push_back(t % 10);
        t /= 10;
        
        if (t) C.push_back(1);
    }
    
    return C;
}

1.2 高精度减法

在这里插入图片描述

#include <iostream>
#include <vector>

using namespace std;

const int N = 1010;

// 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 -- )   // 长度相等情况下a 不等于 b的情况
        if (A[i] != B[i])
            return A[i] > B[i];
    return true; // a  == b的情况
}

vector<int> sub(vector<int> &A, vector<int> &B)
{
    vector<int> C;
    
    int t = 0; //表示借位
    
    for (int i = 0; i < A.size(); i ++ )
    {
        t = A[i] - t;
        if (i < B.size()) t -= B[i];
        C.push_back((t + 10) % 10);
        
        if (t < 0) t = 1;
        else t = 0;
    }
    //去除前导0
    
    while (C.size() > 1 && C.back() == 0)   C.pop_back();
    
    return C;
}

int main()
{
    string a, b;
    vector<int> A, B;
    
    cin >> a >> b;
    
    for (int i = a.size() - 1; i >= 0; i -- ) A.push_back(a[i] - '0');
    for (int i = b.size() - 1; i >= 0; i -- ) B.push_back(b[i] - '0');
    
    if (cmp(A, B))
    {
        vector<int> C = sub(A, B);
        for (int i = C.size() - 1; i >= 0; i -- )   printf("%d",C[i]);
    }else
    {
        vector<int> C = sub(B, A);
        printf("-");
        
        for (int i = C.size() - 1; i >= 0; i -- )   printf("%d",C[i]);
    }
    
    return 0;
}

1.3 高精度乘法

#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() || t; i ++ )
    {
        t += A[i] * b;
        C.push_back(t % 10);
        t /= 10;
        
    }
    return C;
}

int main()
{
    string a;
    int b;
    vector<int> A;
    
    cin >> a >> b;
    
    for (int i = a.size() - 1; i >= 0; i -- )  A.push_back(a[i] - '0');
    
    vector<int> C = mul(A, b);
    
    for (int i = C.size() - 1; i >= 0; i -- ) printf("%d",C[i]);
    
    
    return 0;
}

1.4 高精度除法

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <vector>

using namespace std;

const int N = 1e6 + 10;

//A/b 商是C 余数是r
vector<int> div(vector<int> &A, int b, int & r)
{
    vector<int> C;
    r = 0;
    for (int i = A.size() - 1; i >= 0; i--)
    {
        r = r * 10 + A[i];
        C.push_back(r / b);
        r %= b;
    }
    reverse(C.begin(),C.end());
    while(C.size() > 1 && C.back() == 0) C.pop_back();
    
    return C;
}


int main()
{
    string  a;
    int b;
    vector<int> A;
    
    cin >> a >> b;
    for ( int i = a.size() - 1; i >= 0; i--) A.push_back(a[i] - '0');
    
    int r;
    auto C = div(A, b, r);
    
    for (int i = C.size() - 1; i >= 0; i--) printf("%d",C[i]);
    cout << endl << r << endl;
    
}

高精度汇总篇

二、前缀和

2.1 一维前缀和

核心思想：
前缀和

从1开始的好处，处理边界问题，如要求[1,10]这段和的话，之间s10即可。因为s[0] = 0；

#include <iostream>
using namespace std;

const int N = 100010;

int a[N], s[N];

int main()
{
    int n, m;
    ios::sync_with_stdio(false);
    //原理是让cin和标准输入输出不同步 提高cin读取速度
    //副作用不能用scanf了 数据规模大于等于100w建议用scanf 否则用cin即可
    
    cin >> n >> m;
    
    for (int i = 1; i <= n; i ++ ) cin >> a[i];
    for (int i = 1; i <= n; i ++ ) s[i] = s[i - 1] + a[i];
    while (m -- )
    {
        int l, r;
        cin >> l >> r;
        cout << s[r] - s[l - 1] << endl;
    }
        
    return 0;
}

2.2 二前缀和

核心思想：
在这里插入图片描述

#include <iostream>

using namespace std;

const int N = 1010;

int a[N][N];
int s[N][N];

int main()
{
    int n, m, q;
    
    cin >> n >> m >> q;
    
    //读入数组
    for (int i = 1; i <= n; i ++ )  
        for (int j = 1; j <= m; j ++ ) 
            cin >> 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, y1, x2, y2;
        cin >> x1 >> y1 >> x2 >> y2;
        
        cout << s[x2][y2] - s[x1 - 1][y2] - s[x2][y1 - 1] + s[x1 - 1][y1 - 1] << endl;  
    }
    
    return 0;
}

都是从左上到右下。

三、差分

3.1 一维差分

在这里插入图片描述
在这里插入图片描述

#include <iostream>

using namespace std;

const int N = 100010;

int a[N], b[N];
int n, m;

void insert(int &l, int &r, int &c)
{
    b[l] += c;
    b[r + 1] -= c;
}
int main()
{
    scanf("%d%d", &n, &m);
    
    //读入数据
    for (int i = 1; i <= n; i ++ ) scanf("%d", &a[i]);
    
    //进行构造b数组
    for (int i = 1; i <= n; i ++ ) insert(i, i, a[i]);
    
    // for (int i = 1; i <= n; i ++ ) cout << b[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];
    
    //输出原数组
    for (int i = 1; i <= n; i ++ ) printf("%d ", b[i]);
    
    return 0;
}

#include <iostream>
using namespace std;

const int N = 1010;

int a[N], b[N];

int n, m;

int main()
{
    cin >> n >> m;
    
    for (int i = 1; i <= n; i ++ )  cin >> a[i];
    
    for (int i = 1; i <= n; i ++ )  b[i] = a[i] - a[i - 1];
    
    while (m -- )
    {
        int l, r, c;
        cin >> l >> r >> c;
        b[l] += c;
        b[r + 1] -= c;
    }
    
    for (int i = 1; i <= n; i ++ ) a[i] = b[i] + a[i - 1];

	for (int i = 1; i <= n; i ++ ) cout << a[i] << ' ';
    
	return 0; 
}

3.2 二维差分

在这里插入图片描述

bij是空想处理的，满足aij是b数组的前缀和，给定a数组，假想一个b数组。使得a数组是b数组的前缀和，而b数组则是a数组前缀和的逆运算。

代码:

#include <iostream>

using namespace std;

const int N = 1010;

int a[N][N], b[N][N];
int n, m, q;

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()
{
    cin >> n >> m >> q;
    
    //读入数据
    for (int i = 1; i <= n; i ++ )
        for (int j = 1; j <= m; j ++ ) 
            cin >> 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]);
    
    //做q个询问
    while ( q -- )
    {
        int x1, y1, x2, y2, c;
        cin >> x1 >> y1 >> x2 >> y2 >> c;
        insert(x1, y1, x2, y2, c);
    }
    
    //求a数组
    for (int i = 1; i <= n; i ++ )
        for (int j = 1; j <= m; j ++ )
            a[i][j] = a[i - 1][j] + a[i][j - 1] - a[i - 1][j - 1] + b[i][j];
    
    for (int i = 1; i <= n; i ++ )
    {
       for (int j = 1; j <= m; j ++ )
       {
            cout << a[i][j] << ' ';
       }
       puts("");
    }
    return 0;
}

四、习题汇总

类型题目
二分AcWing 789. 数的范围
AcWing 790. 数的三次方根
高精度AcWing 791. 高精度加法
AcWing 792. 高精度减法
AcWing 793. 高精度乘法
AcWing 794. 高精度除法
前缀和与差分AcWing 795. 前缀和
AcWing 796. 子矩阵的和
AcWing 797. 差分
AcWing 798. 差分矩阵
寒夜点孤灯
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
C++算法：前缀和基础
10-15
比如：nums = {1,2,3,4},则preSum = {0,1,3,6,10}。通过preSum，我们可以求任意nums的子数组和。子数组[i,j)等于子数组[0,j)减去[0,i)，也就是子数组[i,j)的和等于preSum[j] – preSum[i]。如果i等于j，则preSum[i]-...
C++高精度算法
一个平平无奇的博客，一个平平无奇的人，一个平平无奇的约定
10-04 6156
C++高精度加、减、乘、除 大家可以直接拿走，欢迎转载
C++算法(高精度算法)
博客
05-01 3563
C++高精度算法
C++算法之高精度计算
一个程序猿~
04-30 2660
利用计算机进行计算时，有时会遇到计算精度不够的问题。但是我们可以利用程序实现高精度计算，本文介绍几种常用的高精度计算方法。
c++高精度模板
名字不能相同
07-30 333
两个大整数相乘 #include <iostream> #include <cstring> using namespace std; const int N = 50000; char a[N], b[N]; int A[N], B[N], C[2*N]; int main() { cin >> a >> b; A[0] = strlen(a); B[0] = strlen(b); for(int i=
前缀和与差分算法地深入理解
07-11
该内容包括：一维前缀和、二维前缀和、差分和差分矩阵
基于FFT谱分析算法的高精度相位差测量方法.pdf
08-29
基于FFT谱分析算法的高精度相位差测量方法
bsdiff差分算法最新版.7z
04-26
bsdiff差分算法最新版，能够实现在下位机中差分升级的开源算法，在git上也能下载到。
基于改进差分进化的高精度粒子滤波算法
03-03
基于改进差分进化的高精度粒子滤波算法
C++ 算法篇 高精度
流年的博客
11-17 4253
第一章 高精度计算 利用计算机进行数值计算，有时会遇到这样的问题：有些计算要求精度高，希望计算的数的位数可达几十位甚至几百位，虽然计算机的计算精度也算较高了，但因受到硬件的限制，往往达不到实际问题所要求的精度。我们可以利用程序设计的方法去实现这样的高精度计算。介绍常用的几种高精度计算的方法。 高精度计算中需要处理好以下几个问题： (1)数据的接收方法和存贮方法: 第一种方法是：当输入的数很长时，可采用字符串方式输入，这样可输入数...
c++中的高精度算法
2302_79027886的博客
11-19 463
虽然这些都有点复杂和蠢，但是对于刚开始接触算法的萌新来说（比如我）还是很好理解的，对于以上的代码，无一不用到了STL里的函数，所以，建议去学习一下STL里的函数，这玩意儿对于写代码还是很好用的。最后，希望本篇文章对您有帮助，虽然很简陋。如果有更好的思路，可以的话希望能提出，毕竟是学术性的交流，虽然作者很弱。。。
高精度算法（C/C++）
今を生きる
03-24 1570
高精度算法 （C语言） 做ACM题的时候，经常遇到大数的加减乘除，乘幂，阶乘的计算，这时给定的数据类型往往不够表示最后结果，这时就需要用到高精度算法。高精度算法的本质是把大数拆成若干固定长度的块，然后对每一块进行相应的运算。这里以考虑4位数字为一块为例（也可以考虑其他位，但要注意在每一块进行相应运算时不能超出数据类型的数值范围）。 1. 高精度加法 以3479957928375817+897259...
C++基础算法高精度篇
本科在读，致力于Linux/C++方向的学习，同时涉及数据结构与算法，操作系统，数据库mysql等等
07-19 1086
主要讲解了高精度算法的四种常用的计算
高精度算法（C++算法）
m0_73059696的博客
08-29 334
对高精度的加减乘除算法的介绍
C++不知算法系列之高精度数值处理算法
y6123236的专栏
11-25 1755
什么是高精度数值处理算法？高精度数值指因受限于计算机硬件的制约，超过计算机所能存储范围的数值。既然不能存储，更谈不上运算。对此类数值的加、减、乘、除运算需要提供针对性的算法方能获取到结果。此类算法的设计思路因有别于其它算法，为了研究的方便，称此类算法为高精度数值处理算法。本文将讲解如何实现对此类数值的加、减、乘、除运算。本文讲解了高精度相加、相减、相乘、相除操作。面对大数值运算时，除了要有好的算法，还需要有低层硬件的支持。在实际进行超高精度数值运算时，可能需要集群环境下并发运算能力的支持。
c++——爆肝7000字的高精度算法大合集&&实战应用
屁孩君的博客呀！！
08-27 726
学一种算法，多一种思路！！！
python前缀和差分
最新发布
04-20
Python中的前缀和和差分是两种常用的算法技巧，用于解决一些与数组或序列相关的问题。 1. 前缀和（Prefix Sum）： 前缀和是指将数组中的元素依次相加得到的新数组。通过计算前缀和，可以在O(1)的时间复杂度内得到...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值