一、高精度
1、计算形式
(1)、A+B(两个大整数相加)位数为1+e6
模拟加法计算过程如下:
(2)、A-B(两个大整数相减)位数为1+e6
模拟加法计算过程如下:
(3)、A*a(一个大整数乘一个小整数) 大整数位数小于等于1+e6,小整数小于等于1000
(4)、计算A/a的商
大整数的储存(加减乘除大整数的存储都是一致的)
(1)大整数逆序存入数组,数组第0位存个位数,数组每一位存一个数字
(2)模拟计算过程如下
二、高精度代码实列
1、计算A + B
代码示例如下:
#include <iostream>
#include <vector>
using namespace std;
const int N = 1e6 + 10;//防止越界
//模板 C = A + B
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];//超过范围时该位不存在记为0则不用加入
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;//vector数组里面自带获取长度的函数size()
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);//编译器推断C的类型,等价于vector<int>
for(int i = C.size() - 1; i >= 0; i --) printf("%d", c[i]);
return 0;
}
注:计算大整数相加的代码模板(理解记忆模板)
大整数相加模板一:
//模板 C = A + B
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];//超过范围时该位不存在记为0则不用加入
if(i < B.size()) t += B[i];
C.push_back(t % 10);
t /= 10;//判断是否向前进位
}
if(t) C.push_back(1);//最高位是否有进位
return C;
}
大整数相加模板二:改变一下判断条件
//模板 C = A + B
vector<int> add(vector<int> &A, vector<int> &B)//提高效率
{
vector<int> C;
if(A.size() < B.size()) return add(B, A);//进入循环判断时可以少判断一个数组的长度是否超过
int t = 0;// 上一个的进位
for(int i = 0; i < A.size() || i < B.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;
}
2、计算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两个数相等
}
//模板 C = A - B
vector<int> add(vector<int> &A, vector<int> &B)//提高效率
{
vector<int> C;
for(int i = 0,t = 0;i < A.size(); i ++)//定义t为进位
{
t = A[i] - t;//计算当前这一位的值
if(i < B.size()) t -= B[i];//有B这一位才进行相减
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;//vector数组里面自带获取长度的函数size()
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))
{
auto C = sub(A, B);
for(int i = C.size() - 1; i >= 0; i --) printf("%d", C[i]);
}
else
{
auto C = sub(B, A);
printf("-")
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;
for(int i = 0,t = 0;i < A.size(); i ++)//定义t为进位
{
t = A[i] - t;//计算当前这一位的值
if(i < B.size()) t -= B[i];//有B这一位才进行相减
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;
}
3、计算A * b
代码示例如下:
#include <iostream>
#include <vector>
using namespace std;
//C = A * b
vector<int>mul(vector<int> &A, int b)//一个高精度的整数乘以一个低精度的整数
{
vector<int> C;//定义一个答案
int t = 0;//进位
for(int i = 0; i < A.size() | t ; i ++)//A没有循环完或者t不为零就继续做
{
if(i < A.size()) t += A[i] * b;
C.push_back(t % 10);
t /= 10;
}
return C;
}
int main()
{
string a;
int b;
cin >> a >> b;
vector<int> A;
for(int i = a.size() - 1; i >= 0; i --) A.push_back(a[i] - '0'); //从高位进行转化 ,-'0'将字符转化为整数
auto C = mul(A, b);
for(int i = C.size() - 1; i >= 0; i --) printf("%d", C[i]);
return 0;
}
大整数乘以小整数模板
//C = A * b
vector<int>mul(vector<int> &A, int b)//一个高精度的整数乘以一个低精度的整数
{
vector<int> C;//定义一个答案
int t = 0;//进位
for(int i = 0; i < A.size() | t ; i ++)//A没有循环完或者t不为零就继续做
{
if(i < A.size()) t += A[i] * b;
C.push_back(t % 10);
t /= 10;
}
return C;
}
4、计算A / b
代码示例如下:
#include <iostream>
#include <vector>
using namespace std;
//C = 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.beegin(), C.end());//从高位开始储存
while(C.size() > 1 && C.back() == 0) C.pop_back();//长度大于一,并且高位有零,去除前导零
return C;
}
int main()
{
string a;
int b;
cin >> a >> b;
vector<int> A;
for(int i = a.size() - 1; i >= 0; i --) A.push_back(a[i] - '0'); //从高位进行转化 ,-'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;
return 0;
}
大整数除以小整数模板
/C = A / b
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.beegin(), C.end());//从高位开始储存
while(C.size() > 1 && C.back() == 0) C.pop_back();//长度大于一,并且高位有零,去除前导零
return C;
}
二、前缀和
1、一维前缀和(利用公式)
时间复杂度处理前为O(n),相减时为O(1),一维前缀和的计算公式 s[r] - s[l -1]
模板分析如下:
代码实例
一维前缀和(利用公式)
时间复杂度处理前为O(n),相减时为O(1)
一维前缀和的计算公式 s[r] - s[l -1]
代码实例
#include <isotream>
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;
}
注: ios::sync_with_stdio(false);cin和标准输入输出不同步,可以提高cin的读取速度但是速度还是没有scanf快
2、二维前缀和
模板分析如下:
代码示例
#include <iostream>
using namespace std;
const int N = 1010;
int n, m, q;
int 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", &s[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];//计算每一个前缀和
while (q -- )
{
int x1, y1, x2, 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;
}
三、差分(前缀和的逆运算)
1、一维差分
模板原理分析如下:
代码示例如下:
#include <iostream>
using namespace std;
const int N = 100010;
int n, m;
int a[N], b[N];
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]);
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];//用前缀和求原数组A
for (int i = 1; i <= n; i ++ ) printf("%d ", b[i]);
return 0;
}
一维差分核心公式:
将原数组看作为一个全零的数组在进行数的插入
void insert(int l, int r, int c)
{
b[l] += c;
b[r + 1] -= c;
}
2、二维差分
模板原理分析如下:
代码示例如下:
#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 ++ )//先假定为0在进行插入
for (int j = 1; j <= m; j ++ )
insert(i, j, i, j, a[i][j]);
while (q -- )// 进行q次操作
{
int x1, y1, x2, y2, c;
cin >> 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];
for (int i = 1; i <= n; i ++ )
{
for (int j = 1; j <= m; j ++ ) printf("%d ", b[i][j]);
puts("");
}
return 0;
}
二维差分核心代码
与一维差分处理类似先将二维数组看成一个全零数组然后再进行插入
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;
}
代码转载于
作者:yxc
链接:https://www.acwing.com/activity/content/code/content/39800/
来源:AcWing