高精度计算采用模拟列竖式法。
通常计算用正数,对于负数,处理符号位。html
基本思想
用数组存放和表示高精度整数。一个数组元素,存放高精度整数中的一位。
为了便于计算,将数由低位到高位依次存在数组下标对应由低到高的位置上,其中,数组的第0个下标对应位置存储该数的位数。算法
#include<stdio.h> #include<string.h> #define MAX_LEN 200 int sz[MAX_LEN + 10]; char szLine[MAX_LEN + 10]; int main() { int len, i; scanf("%s", szLine); memset(sz, 0, sizeof(sz));/*数组清零*/ sz[0] = strlen(szLine); for (i = 1; i <= sz[0]; i++) { sz[i] = szLine[sz[0] - i] - '0'; } return 0; }
int cmp(int a[], int b[]) { int i; if (a[0] > b[0]) { return 1; } else if (a[0] < b[0]) { return -1; } for (i = a[0]; i >= 1; i--) { if (a[i] > b[i]) { return 1; } else if (a[i] < b[i]) { return -1; } } return 0; }
void plus(int a[], int b[]) { int i; if (a[0] < b[0]) { a[0] = b[0]; } for (i = 1; i <= a[0]; i++) { a[i] += b[i]; if (a[i] >= 10) { a[i + 1]++; a[i] -= 10; } } if (a[a[0] + 1] != 0) { a[0]++; } }
int minus(int a[], int b[]) { int i; if (cmp(a, b) == 0) { memset(a, 0, sizeof(a)); a[0] = 1; return 0; } else if (cmp(a, b) > 0) { for (i = 1; i <= a[0]; i++) { a[i] = a[i] - b[i]; if (a[i] < 0) { a[i + 1] -= 1; a[i] += 10; } } while (a[a[0]] == 0) { a[0]--; } return 1; } else { minus(b, a); for (i = 0; i <= b[0]; i++) { a[i] = b[i]; } return -1; } }/*返回值表示减法结果的符号*/
参考资料:从零开始学算法:高精度计算->强推数组