高精度計算1

高精度計算在最近幾天上去，現貼上的只有高精度階乘、加法、減法，加減法支持負數及小數，階乘最多支持10000位，要增加請更改a數組的大小，加減法最多支持100位，若要修改，修改MAXN的大小

// 趣味C.cpp : 定義控制檯應用程序的入口點。 // 高精度的總結~ #include <stdio.h> #include <string.h> #include <stdlib.h> #include <memory.h> #define MAXN 100 void addition(char *num1, char *num2); void Subtract(char *num1, char *num2); // 階乘 void fact(int n) { int a[10000]; int nBit = 1; // 位數 int temp; // 階數的任一元素與臨時結果的某位的乘積結果 int i, j, carry; // 進位 // 從1*2開始 a[0] = 1; for (i = 2; i <= n; i++) { carry = 0; // 每次算完都把進位清零 for (j = 0; j < nBit; j++) // 按被乘數的位數來循環 { temp = a[j] * i + carry; a[j] = temp % 10; carry = temp / 10; // 存進位 } // 如果有進位digit位數就相加 while (carry) { ++nBit; a[nBit - 1] = carry % 10; // 存進位 carry /= 10; // 看還能不能進位 } } printf("%d! = ", n); for (j = nBit; j >= 1; j--) printf("%d", a[j-1]); printf("/n/n"); } // 初始化數據 void InitData(char *num1, char *num2) { // 計算小數點的位置 size_t PosDotNum1 = 0, PosDotNum2 = 0; int i, diff; for (i = 0; num1[i] != '.' && num1[i] != '/0'; i++) PosDotNum1++; if (num1[i] == '/0') { num1[i] = '.'; num1[i + 1] = '/0'; } for (i = 0; num2[i] != '.' && num2[i] != '/0'; i++) PosDotNum2++; if (num2[i] == '/0') { num2[i] = '.'; num2[i + 1] = '/0'; } // 使兩個數的小數點對齊 if (PosDotNum1 < PosDotNum2) { diff = PosDotNum2 - PosDotNum1; for (i = strlen(num1); i >= 0; i--) num1[i + diff] = num1[i]; // 前導補0 for (i = 0; i < diff; i++) num1[i] = '0'; } else if (PosDotNum2 < PosDotNum1) { diff = PosDotNum1 - PosDotNum2; for (i = strlen(num2); i >= 0; i--) num2[i + diff] = num2[i]; // 前導補0 for (i = 0; i < diff; i++) num2[i] = '0'; } // 後繼補0 for (i = strlen(num1); i < MAXN - 1; i++) num1[i] = '0'; num1[MAXN - 1] = '/0'; for (i = strlen(num2); i < MAXN - 1; i++) num2[i] = '0'; num2[MAXN - 1] = '/0'; } // 打印結果 void PrintResult(char *num) { int i, n; for (i = 0; num[i] == '0'; i++) ; // 去除後繼零 for (n = MAXN - 1; num[n] == '0'; n--) ; if (-1 == n) { printf("0/n"); return; } if ('.' == num[n]) num[n] = '/0'; else num[n + 1] = '/0'; if ('.' == num[i]) i--; for (; num[i] != '/0'; i++) printf("%c", num[i]); printf("/n"); } // 高精度加法 void addition(char *num1, char *num2) { char result[MAXN + 1]; int i, temp; if ('-' == num1[0] && num2[0] != '-') { Subtract(num2, num1 + 1); return; } else if (num1[0] != '-' && num2[0] == '-') { Subtract(num1, num2 + 1); return; } else if (num1[0] == '-' && num2[0] == '-') { printf("-"); addition(num1 + 1, num2 + 1); return; } // 初始化兩個加數與和 InitData(num1, num2); memset(result, '0', sizeof(result) / sizeof(char)); result[MAXN] = '/0'; // 計算和 for (i = MAXN - 2; i >= 0; i--) { if ('.' == num1[i]) { result[i + 1] = '.'; continue; } temp = (result[i + 1] - '0') + (num1[i] - '0') + (num2[i] - '0'); result[i + 1] = temp % 10 + '0'; if ('.' == num1[i - 1] && (temp / 10)) result[i - 1] += temp / 10; result[i] += temp / 10; } // 打印結果 PrintResult(result); } // 高精度減法 void Subtract(char *num1, char *num2) { char result[MAXN + 1]; char temp1[MAXN], temp2[MAXN]; int i, temp; if (num1[0] == '-' && num2[0] != '-') { printf("-"); addition(num1 + 1, num2); return; } else if (num1[0] != '-' && num2[0] == '-') { strcpy(temp2, num2 + 1); addition(num1, temp2); return; } else if (num1[0] == '-' && num2[0] == '-') { strcpy(temp1, num1 + 1); strcpy(temp2, num2 + 1); Subtract(temp1, temp2); return; } // 初始化兩個減數與差 InitData(num1, num2); memset(result, '0', sizeof(result) / sizeof(char)); result[MAXN] = '/0'; if (strcmp(num1, num2) < 0) { printf("-"); Subtract(num2, num1); return; } // 計算差 for (i = MAXN - 2; i >= 0; i--) { if ('.' == num1[i]) { result[i + 1] = '.'; continue; } temp = (num1[i] - '0') - (num2[i] - '0'); if (temp < 0 && i > 0) { temp += 10; if ('.' == num1[i - 1]) num1[i - 2]--; else num1[i - 1]--; } result[i + 1] = temp + '0'; } // 打印計算結果 PrintResult(result); } // 高精度乘法計算 void multiply(char *num1, char *num2) { } int main(int argc, char* argv[]) { int sel; int n; char num1[MAXN], num2[MAXN]; printf("===========================================/n"); printf(" 歡迎使用高精度計算器 Dos版 VS.JQCTBU/n"); printf("===========================================/n"); while (1) { printf("/n/n1.階乘/n"); printf("2.加法/n"); printf("3.減法/n"); printf("4.乘法/n"); printf("5.除法/n"); printf("6.多項式計算（簡單）/n"); printf("7.多項式計算（高精度）/n"); printf("0.退出/n/n"); printf("請輸入你的選擇："); scanf("%d", &sel); num1[0] = num2[0] = '/0'; switch (sel) { case 1: printf("/n請輸入一個數：/n"); scanf("%d", &n); fact(n); break; case 2: printf("/n請輸入兩個數/n第一個數："); scanf("%s", num1); printf("第二個數："); scanf("%s", num2); printf("%s + %s = /n", num1, num2); addition(num1, num2); break; case 3: printf("/n請輸入兩個數/n被減數："); scanf("%s", num1); printf("減數："); scanf("%s", num2); printf("%s - %s = /n", num1, num2); Subtract(num1, num2); break; case 4: printf("/n請輸入兩個需要相乘的數/n第一個數："); scanf("%s", num1); printf("第二個數："); scanf("%s", num2); multiply(num1, num2); break; case 0: exit(0); break; } } return 0; }