C++高精度运算(加减乘除)
创始人
2024-05-13 00:34:27

C++中在int 是32位bit是2^32-1。

C/C++基础数字类型的取值范围
类型    取值范围    位数
unsigned   int    0~4294967295    10位
int    2147483648~2147483647    10位
unsigned long    0~4294967295    10位
long    2147483648~2147483647    10位
long long    -9223372036854775808~9223372036854775807    19位
unsigned long long    0~1844674407370955161    19位
__int64    -9223372036854775808~9223372036854775807    19位
unsigned __int64    0~18446744073709551615    20位

用高精度可以计算的更大。

下面是加法运算:

#define _CRT_SECURE_NO_WARNINGS 1	
#include
#include
#include
using namespace std;
const int N = 100005;
int A[N], B[N], C[N];
int la, lb, lc;
void add(int A[], int B[], int C[]) {for (int i = 0; i < lc; i++) {C[i] += A[i] + B[i]; //累加C[i + 1] += C[i] / 10; // 进位C[i] %= 10; //存余数}if (C[lc]) lc++; // 注意最高位的处理 lc是位长
}
int main() {string a, b;cin >> a >> b;la = a.size(), lb = b.size(), lc = max(la, lb); // 这里是取最高位for (int i = la - 1; ~i; i--) A[la - 1 - i] = a[i] - '0';for (int i = lb - 1; ~i; i--) B[lb - 1 - i] = b[i] - '0';add(A, B, C);for (int i = lc - 1; ~i; i--) printf("%d", C[i]);return 0;
}

减法运算:

注意要用大数减小数,如果是小数减大数要加负号.在减的时候要用while确定不为0的位。

在主函数中要注意for不需要减一。

#define _CRT_SECURE_NO_WARNINGS 1	
#include
#include
#include
using namespace std;
const int N = 100005;
int A[N], B[N], C[N];
int la, lb, lc;
bool cmp(int A[], int B[]) {if (la != lb) return la > lb;for (int i = la - 1; ~i; i--) {if (A[i] != B[i]) return A[i] > B[i];}return true;
}
void sub(int A[], int B[], int C[]) {for (int i = 0; i < lc; i++) {if (A[i] < B[i]) {A[i + 1]--, A[i] += 10;}C[i] = A[i] - B[i];}while (lc && C[lc] == 0) lc--;
}
int main() {string a, b;cin >> a >> b;la = a.size(), lb = b.size(), lc = max(la, lb); // 这里是取最高位for (int i = la - 1; ~i; i--) A[la - 1 - i] = a[i] - '0';for (int i = lb - 1; ~i; i--) B[lb - 1 - i] = b[i] - '0';if (!cmp(A, B)) swap(A, B), cout << "-";sub(A, B, C);for (int i = lc ; ~i; i--) printf("%d", C[i]);return 0;
}

乘法运算:

#define _CRT_SECURE_NO_WARNINGS 1	
#include
#include
#include
using namespace std;
const int N = 100005;
int A[N], B[N], C[N];
int la, lb, lc;
void mul(int A[], int B[], int C[]) {for (int i = 0; i < la; i++) {for (int j = 0; j < lb; j++) {C[i + j] += A[i] * B[j]; // 这里要注意是加C[i + j + 1] += C[i + j] / 10;C[i + j] = C[i + j] % 10;}}while (lc && C[lc] == 0) lc--;
}
int main() {string a, b;cin >> a >> b;la = a.size(), lb = b.size(), lc = la+lb; // 这里是取最高位for (int i = la - 1; ~i; i--) A[la - 1 - i] = a[i] - '0';for (int i = lb - 1; ~i; i--) B[lb - 1 - i] = b[i] - '0';mul(A, B, C);for (int i = lc; ~i; i--)printf("%d", C[i]);return 0;
}

除法运算:

#define _CRT_SECURE_NO_WARNINGS 1	
#include
#include
#include
using namespace std;
const int N = 100005;
int A[N], B[N], C[N];
int la, lb, lc;
void mul(int A[], int B[], int C[]) {for (int i = 0; i < la; i++) {for (int j = 0; j < lb; j++) {C[i + j] += A[i] * B[j]; // 这里要注意是加C[i + j + 1] += C[i + j] / 10;C[i + j] = C[i + j] % 10;}}while (lc && C[lc] == 0) lc--;
}void div(int A[], int b, int C[]) {long long r = 0;for (int i = la - 1; ~i; i--) {r = r * 10 + A[i];C[la - 1 - i] = r / b; //从后向前r %= b;}reverse(C, C + lc);while (lc && C[lc] == 0) lc--;
}int main() {string a;int b;cin >> a >> b;la = a.size(),  lc = la; // 这里是取最高位for (int i = la - 1; ~i; i--) A[la - 1 - i] = a[i] - '0';//for (int i = lb - 1; ~i; i--) B[lb - 1 - i] = b[i] - '0';//mul(A, B, C);div(A, b, C);for (int i = lc; ~i; i--)printf("%d", C[i]);return 0;
}

相关内容

热门资讯

【硅锰】矿价强势、钢厂追涨,供... (来源:合金现货网)本周硅锰市场大体呈强势上涨运行,从消息面情绪主导到供需基本面迎来实质性改善预期,...
伊朗总统:若基础设施或经济中心... 总台记者获悉,当地时间28日,伊朗总统佩泽希齐扬表示,如果伊朗的基础设施或经济中心成为攻击目标,伊方...
@辽宁驾驶人 6月1日起,这些... (来源:东北新闻网)近日公安部正式发布《机动车驾驶人疲劳驾驶认定规则》(GA/T 2372-2026...
“增长失速、收入迷雾、家族掌控... 一家宣称在高炉本体内衬细分市场占据半壁江山的公司,一份净利润两年高增长的成绩单,一个紧扣“碳中和”热...