大数运算

由于编程语言提供的基本数值数据类型表示的数值范围有限，不能满足较大规模的高精度数值计算，因


此需要利用其他方法实现高精度数值的计算，于是产生了大数运算。大数运算主要有加、减、乘三种方


法。
中文名 大数运算 概 ? ?念 很大的数值的数进行一系列的运算 学 ? ?科 高等数学 应 ? ?用 空间力学
目录
1 简介
2 原理
3 应用
? HDNOJ——1002 A + B Problem II
? HDNOJ 1002 代码
简介
大数运算，顾名思义，就是很大的数值的数进行一系列的运算。
我们知道，在数学中，数值的大小是没有上限的，但是在计算机中，由于字长的限制，计算机所能表示


的范围是有限的，当我们对比较小的数进行运算时，如：1234+5678，这样的数值并没有超出计算机的表


示范围，所以可以运算。但是当我们在实际的应用中进行大量的数据处理时，会发现参与运算的数往往


超过计算机的基本数据类型的表示范围，比如说，在天文学上，如果一个星球距离我们为100万光年，那


么我们将其化简为公里，或者是米的时候，我们会发现这是一个很大的数。这样计算机将无法对其进行


直接计算。
可能我们认为实际应用中的大数也不过就是几百位而已，实际上，在某些领域里，甚至可能出现几百万


位的数据进行运算，这是我们很难想象的。如果没有计算机，那么计算效率可想而知。
既然在计算机中无法直接表示，那么大数到底如何进行运算呢，学习过数据结构的都知道线性表，将大


数拆分然后存储在线性表中，不失为一个很好的办法。
大数除法，可以类比人类手算，添位比较取商，中间结果与除数相减所得到的差参与下一轮运算，直到


结束。这里需要注意添位时，商有时需要补零，而且除法运算需要使用加、减、乘、比较运算。
原理编辑
利用数组连续性，将大数每一位上的数字单独取出放入对应的数组格中，然后再对每一位做单独的加减


乘运算。形象的说，类似于小学学习加减乘所列的式子。
应用
HDNOJ——1002 A + B Problem II
Time Limit: 2000/1000 MS (Java/Others)Memory Limit: 65536/32768 K (Java/Others)
　　Total Submission(s): 69615Accepted Submission(s): 12678
　　
　　
Problem Description
I have a very simple problem for you. Given two integers A and B,your job is to calculate?


the Sum of A + B.
Input
The first line of the input contains an integer T(1<=T<=20) which means the number of test?


cases. Then T lines follow,each line consists of two positive integers,A and B. Notice?


that the integers are very large,that means you should not process them by using 32-bit?


integer. You may assume the length of each integer will not exceed 1000.
Output
For each test case,you should output two lines. The first line is "Case #:",# means the?


number of the test case. The second line is the an equation "A + B = Sum",Sum means the?


result of A + B. Note there are some spaces int the equation. Output a blank line between?


two test cases.
Sample Input
2
1 2
112233445566778899 998877665544332211
Sample Output
Case 1:
1 + 2 = 3
Case 2:
112233445566778899 + 998877665544332211 = 2222222222222221110
Author
Ignatius.L
HDNOJ 1002 代码


#include <stdio.h>
#include <string.h>
?
int main()
{
char str1[1001],str2[1001];
int t,i,len_str1,len_str2,len_max,num=1,k;
scanf("%d",&t);
while(t--)
? ? {
? ? ? ? int a[1001]={0},b[1001]={0},c[1001]={0};
? ? ? ? scanf("%s",str1);//输入第一个大数
? ? ? ? len_str1=strlen(str1);//算出str1的长度，包括''
? ? ? ? for(i=0;i<=len_str1-1;++i)
? ? ? ? a[i]=str1[len_str1-1-i]-'0';//将str1的值逆序赋值给a
? ? ? ? scanf("%s",str2);//输入第二个大数
? ? ? ? len_str2=strlen(str2);//算出str2的长度，包括''
? ? ? ? for(i=0;i<=len_str2-1;++i)
? ? ? ? b[i]=str2[len_str2-1-i]-'0';//将str2的值逆序赋值给b
? ? ? ? if(len_str1>len_str2)//判断谁的长度长，长度长的复制到len_max
? ? ? ? ? ? len_max=len_str1;
? ? ? ? else
? ? ? ? ? ? len_max=len_str2;
? ? ? ? k=0;
? ? ? ? for(i=0;i<=len_max-1;++i)
? ? ? ? {
? ? ? ? ? ? c[i]=(a[i]+b[i]+k)%10;//取每位上两数相加的个位和
? ? ? ? ? ? k=(a[i]+b[i]+k)/10;//将两数相加的十位数赋值到k
? ? ? ? }
? ? ? ? if(k!=0)//如果k还有剩余，说明数组c还有一位，而且不超过1
? ? ? ? ? ? c[len_max]=1;
? ? ? ? printf("Case%d:n",num);//输出Case#:
? ? ? ? num++;
? ? ? ? printf("%s+%s=",str1,str2);//输出a+b=
? ? ? ? if(c[len_max]==1)//输出最高位
? ? ? ? ? ? printf("1");
? ? ? ? for(i=len_max-1;i>=0;--i)//输出最高位以下的数
? ? ? ? {
? ? ? ? ? ? printf("%d",c[i]);
? ? ? ? }
? ? ? ? printf("n");
? ? ? ? if(t>=1)
? ? ? ? ? ? printf("n");
? ? ? ? }
? ? return0;
?
}
?
//上面代码是错的，WA，格式太差懒得改了，贴个对的。
#include <stdio.h>
#include <string.h>
?
int main()
{
? ? char str1[1001],str2[1001];
? ? int t,k;
? ? scanf("%d",&t);
? ? while(t--)
? ? {
? ? ? ? int a[1001]={0},str1);
? ? ? ? len_str1=strlen(str1);
? ? ? ? for(i=0;i<=len_str1-1;++i)
? ? ? ? ? ? a[i]=str1[len_str1-1-i]-'0';
? ? ? ? scanf("%s",str2);
? ? ? ? len_str2=strlen(str2);
? ? ? ? for(i=0;i<=len_str2-1;++i)
? ? ? ? ? ? b[i]=str2[len_str2-1-i]-'0';
? ? ? ? if(len_str1>len_str2)
? ? ? ? ? ? len_max=len_str1;
? ? ? ? else
? ? ? ? ? ? len_max=len_str2;
? ? ? ? k=0;
? ? ? ? for(i=0;i<=len_max-1;++i)
? ? ? ? {
? ? ? ? ? ? c[i]=(a[i]+b[i]+k)%10;
? ? ? ? ? ? k=(a[i]+b[i]+k)/10;
? ? ? ? }
? ? ? ? if(k!=0)
? ? ? ? ? ? c[len_max]=1;
? ? ? ? printf("Case%d:n",num);
? ? ? ? num++;
? ? ? ? printf("%s+%s=",str2);
? ? ? ? if(c[len_max]==1)
? ? ? ? ? ?printf("1");
? ? ? ? for(i=len_max-1;i>=0;--i)
? ? ? ? {
? ? ? ? ? ? printf("%d",c[i]);
? ? ? ? }
? ? ? ? printf("n");
? ? ? ? //if(t>=1)
? ? ? ? //printf("n");
? ? }
? ? return 0;
}
========

大数相加

http://blog.csdn.net/hgy2011/article/details/8635251
问题描述：实现A+B=C ? 其中A、B位数超过100位
算法思想：大数使用字符串存储，每一个单元存储操作数的每一位，之后执行位相加。
基本思路：字符串反转、字符变数字、位运算、反序输出
C语言代码：
#include<stdio.h> ?
#include<string.h> ?
#define Max 101 ?
void print(char sum[]); ?
void bigNumAdd(char a[],char b[],char sum[]); ?
int main() ?
{ ?
? ? char a[Max]; ?
? ? char b[Max]; ?
? ? char sum[Max]; ?
? ? gets(a); ?
? ? gets(b); ?
? ? bigNumAdd(a,b,sum); ?
? ? print(sum); ?
? ? return 0; ?
} ?
? ?
void bigNumAdd(char a[],char sum[]) ?
{ ?
? ? int i=0; ?
? ? int c=0;//表示进位 ?
? ? ? ? ? //初始化,对以后位运算有很大帮助！ ?
? ? char m[Max]={0}; ?
? ? char n[Max]={0}; ?
? ? memset(sum,Max*sizeof(char)); //这里不能写成memset(sum,sizeof(sum));原因见注意事项1?


?
? ? //字符串反转且字符串变数字 ?
? ? int lenA=strlen(a); ?
? ? int lenB=strlen(b); ?
? ? for (i=0;i<lenA;i++) ?
? ? { ?
? ? ? ? m[i]=a[lenA-i-1]-'0'; ?
? ? } ?
? ? for (i=0;i<lenB;i++) ?
? ? { ?
? ? ? ? n[i]=b[lenB-i-1]-'0'; ?
? ? } ?
? ? //位运算 ?
? ? for (i=0;i<lenA||i<lenB;i++) ?
? ? { ?
? ? ? ? sum[i]=(m[i]+n[i]+c)%10+'0';//得到末位 ?
? ? ? ? c=(m[i]+n[i]+c)/10;//得到进位 ?
? ? } ?
} ?
? ?
void print(char sum[]) ?
{ ?
? ? int i=0; ?
? ? int j=0; ?
? ? int len = strlen(sum); ?
? ? for (i=len-1;sum[i]==0;i--); //找到第一个不为零的位置，方便输出 ?
? ? for (j=i;j>=0;j--) ?
? ? { ?
? ? ? ? printf("%c",sum[j]); ?
? ? } ?
} ?


注意事项:
1、要注意sizeof的用法
? ? ? char sum[100],sizeof(sum) =100 (表示整个字符串的大小)
? ? ? char* sum 或 char sum[] (为函数参数),sizeof(sum)=4 ?(表示存储一个指针的单元大小)
? ? ? 主要的区别：sum在定义的时候意义不同
2、memset(sum,Max*sizeof(char)):
? ? ? 作用：一般是用于初始化或清空一个内存块，等价于循环
? ? ?使用范围：百度说：在给char以外的数组赋值时，只能初始化为0或者-1，赋其他的值就要用循环


来实现。对字符数组可以赋任何字符
? ? ? 因为：在计算机里，0的二进制都是0，-1的二进制都是1。
3、在这里，对字符串进行初始化0操作很重要，它会影响到位操作和输出，要引起重视！
4、100位相加，和最大为101位，故申请最大为101
========

大数运算的算法

在数学运算的过程中,经常碰到这样的问题:参与运算的 数字很大或者对运算结果的精度要求很高。无


论何种计算机 语言,在描述数据类型时都有一定的精度和位数要求,比如说 16 位整型数 (int) 所能


表示的范围为 - 32768～32767,实数 (float)所能表示的精度为小数点后 8 位等,每种类型的数据都?


有其精度和位数限制。超过 20 位有效数字的数值一般就无 法表示了。所以,在上述情况下,采用一般


的程序设计无法满 足要求,必须采用高精度的数学运算才能实现。如下例: 示例 1 :
Fibonacci Numbers
A Fibonacci sequence is calculated by adding the previous
two members of the sequence,with the first two members being
both 1. f (1) = 1,f(2) = 1,f (n > 2) = f (n - 1) + f(n - 2)
Your task is to take a number as input,and print that Fi2b onacci number.
?Sample Input
100
?Sample Output
?354224848179261915075
Note: No generated Fibonacci number in excess of 1000 digits will be in the test data,i.?


e. f (20) = 6765 has 4 digits.
这个示例的计算结果可能很大,以至于所有的数据类型 都无法表达。所以我们必须用高精度数学运算相


应的算法进 行求解。
1.大数存储的实现
作为实现大数存储最常见的一类方法是利用数组。将一 个有 n 位的大数存入数组,每个数组的一个元


素表示一位十 进制数,若是 n 没有超过普通 PC 机允许的定义范围,这种算 法是简单易行的。 如果


出现超大数,则可以采用万进制的方法加以存储,在 此就不多做介绍了。
?2.大数计算的算法
?这里仅对大数加法的算法做详细描述,剩下的运算给出 算法,具体程序实现过程留给读者思考。
2. 1 大数加法
?看下面一个例子:
?122345678902345
+ 　 ?34567890012
?122380246792357
在上面的加法运算中,并没有考虑参与运算的数的位数,即使有小数也是一样。每次运算时只是利用加


法运算的规则 对参与运算的每一位进行运算,每次运算都是在 10 以内进 行,并加上了前面的进位。


通过这种运算就可以非常准确的 得到运算结果,并且可以不考虑位数和精度的问题。 考虑到计算机的


存储问题,如果我们采用数组来存储参 与运算的每个加数,则需要将上面参与运算的数看成字符,并?


将原来的数翻转,即采用下面的运算方法:
? ?543209876543221
+ 21009876543
753297642083221
我们看到运算结果与实际结果相反,因此只需要将结果 再翻转一次输出即可得到正确结果。这样做的好


处在于一旦 遇到低位向高位进位时,不会出现存储上的问题。因为最高 位存储在数组的最后一个元素?


,其后的存储单元可用于存放 进位,同时运算是从数组的第一个元素开始的。 如何用 C语言来求解上面


的加法运算呢 ? 我们须首先定 义一个存储结构来存储每一个数。由于我们采用数组来存储 每一个数,


如果数的最大位数未知,则需采用指针方式定义,如果已知则可采用数组的方式来定义,在该示例中,


其最大位 数不超过 1000 位,因此可采用数组方式来定义,但由于每个 数的实际情况各不相同,因此


还必须定义一个变量来表示每 个数实际的位数。其存储结构用 C语言表示如下:
?typedef struct fanum
{
char fnum[1001] ; / / 数 int len; / / 数的长度
}Fanum;
?基于上面的存储结果及算法,我们可采用如下的 C 语言 程序实现两个多位数的加法运算,注意该数为


整数。
?void strAdd(Fanum3sf1,Fanum3sf2)
{
int i,len,flag = 0 ;
Fanum sf3 ;
len = sf1 - > len > sf2 - > len ? sf1 - > len:sf2 - > len;
sf3. len = len; for(i = 0 ;i < len;i + + ){
sf3.fnum[i] = sf1 - > fnum[i] - ’0’ + sf2 - > fnum[i] - ’0 ’
+flag;
flag = sf3.fnum[i]/ 10 ;
sf3.fnum[i] = sf3.fnum[i] %10 +’0 ’;
} if (flag = = 1){
sf3.fnum[i] =’1 ’;
sf3. len + + ;
} for(i = 0 ;i < sf3. len;i + + ){
sf1 - >fnum[i] = sf2 - >fnum[i] ;
sf1 - > len = sf2 - > len;
sf2 - >fnum[i] = sf3.fnum[i] ;
sf2 - > len = sf3.len;
}
return;
}
注意此处用 sf1 存储加法运算的结果,sf2 用于存储运算 前 sf1 中存储的变量,主要是为了能够与上


面的示例更相吻 合。 有了上面的函数,要完成该题所规定的任务,只需要如下 的主程序便可实现。
?int main()
{
FILE*in,*out ;
Fanum f1,f2 ;
long num,i ;
?in =fopen(“ in. txt”,“r” ) ;
out =fopen(“out. txt”,“w” ) ;
?while( fscanf (in,“%d”,&num) ! = EOF)
{ f1. len = 1 ; f2. len = 1 ; for(i = 0 ;i < 1001 ;i + + )
{ f1.fnum[i] =’0 ’; f2.fnum[i] =’0 ’; }
f1.fnum[0] =’1 ’;
f2.fnum[0] =’1 ’;
?while(num ! = 2)
{ strAdd( &f1,&f2) ; num - - ; }
?for(i =f2. len - 1 ;i > = 0 ;i - - )
fprintf (out,f2.fnum[i] - ’0 ’) ;
fprintf (out,“ n” ) ; }
?fclose(in) ; fclose(out) ;
return 1 ;
}
该题输入数据存放在 in. txt 文件中,而输出结果存放在 out. txt 文件中。 大数加法运算实现算法


如下:
?(1)将 A、 B 按位对齐;
?(2)低位开始逐位相加;
(3)对结果做进位调整。
2. 2 大数减法
?大数减法运算实现算法如下:
(1)将 A、 B 按位对齐;
(2)低位开始逐位相减;
(3)对结果做借位调整。
2. 3 大数乘法
?大数乘法运算实现算法如下:
?(1)引入 sum2 、 sum1 中间过渡量;
(2)在 n 的每一位上处理 m;
(3)通过每一层循环,实现乘法的加法化;
(4)对结果做进位调整。
?2. 4 大数除法
大数除法运算实现算法如下:
?(1)引入 al 来标识 a 的长度,bl 来标识 b 的长度;
?(2)测算 a 和 b 的长度;
(3)高位开始,对位做减法,并完成借位;
?(4)高位开始逐位计算商; (5)整理商,产生余数。
?2. 5 大数取模
在取模运算中用到了上面的除法运算,只需返回余数即 可。
?3.结论
大数运算实际上都是建立在位运算的基础上的,实际上 就是对存储大数的数组里每一个数组项进行操作


。将大数翻 转过来存储是为了防止在计算过程中出现进位不好表示和溢 出的问题。这种处理也可以看


作是对字符串的处理,即把这 个数组里的每一位数字转化成字符进行处理。同时如果涉及 到小数,如


果是加法运算,可以把小数点前后的数据分开处 理;如果是乘法运算,可将其转换成整数运算之后,再


转换成 小数。
参考文献: [1]王永祥.超高精度超大数算法与程序设计[M].陕西:西安交通 大学出版社. 1990. 6. [2]


王昌锐.大数论[M]. 台北:徐氏基金会. 1970. [3]王金荣 等.大数模乘算法的分析与研究[J ]. 计算机


工程与应 用. 2004. 24. [4]Robert L. Kruse,Alexander J. Ryba.“Data Structures And Pro2g?


ram Design In C+ + ” . Higher Education Press Pearson Education. 5,2001.


}
========

大数运算（加减乘除，取模，乘方）

http://blog.csdn.net/agoago_2009/article/details/6862290 #include "Big__CALC.h" ? ?? int main() ? { ? ? ? char str1[200]="9876543210987"; ? ?? ? ? char str2[200]="1234567890"; ? ? ?? ? ? char result[200]={0}; ? ?? ? ? printf("str1 ? ? ? ?= %st%dn",strlen(str1)-1); ? ? ? printf("str2 ? ? ? ?= %st%dn",str2,strlen(str2)-1); ? ?? ? ? int ret = Compare(str1,str2); ? ? ? if (ret==1) ? ? ? { ? ? ? ? ? printf("str1 > str2n"); ? ? ? } ? ? ? else if(ret==-1) printf("str1 < str2n"); ? ? ? else printf("str1 = str2n"); ? ? ? ?? ? ? Add(str1,result); ? ? ? printf("str1 + str2 = %st%dn",result,strlen(result)-1); ? ? ? ?? ?? ? ? Sub(str1,result); ? ? ? printf("str1 - str2 = %st%dn",strlen(result)-1); ? ? ? ?? ?? ? ? Chen(str1,result); ? ? ? printf("str1 * str2 = %st%dn",strlen(result)-1); ? ?? ? ? char div_result[200]={0}; ? ? ? Mod(str1,div_result); ? ? ? printf("str1 M str2 = %st%dn",strlen(result)-1); ? ? ? printf("str1 / str2 = %st%dn",div_result,strlen(div_result)-1); ? ? ? ? ? ? ? int num=10; ? ? ? char str[200]="999"; ? ? ? ChenFang(str,num,result); ? ? ? ThrowAway_0(str); ? ? ? printf("%s ^ %d = %st%dn",str,strlen(result)-1); ? ?? ? ? getchar(); ? ?? ? ? return 1; ? } ? ?? ?? void Trans(char *str_num1,char *str_num2,char *tempbuf1,char *tempbuf2) ? { ? ? ? int len_num1=0; ? ? ? int len_num2=0; ? ? ? int i=0; ? ? ? while(str_num1[i]!='') ? ? ? { ? ? ? ? ? len_num1++; ? ? ? ? ? i++; ? ? ? } ? ?? // ?printf("字符串1的长度: length1=%dn",len_num1); ? ? ? i=0; ? ? ? while(str_num2[i]!='') ? ? ? { ? ? ? ? ? len_num2++; ? ? ? ? ? i++; ? ? ? } ? ?? // ?printf("字符串2的长度: length2=%dnn",len_num2); ? ?? ? ? tempbuf1[0]='0'; ? ? ? tempbuf2[0]='0'; ? ?? //======================================================================= ? ? ? if(len_num2>=len_num1) ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? //补成相同长度 ? ? ? { ? ? ? ? ? for(i=1;i<=(len_num2-len_num1);i++) ? ? ? ? ? { ? ? ? ? ? ? ? tempbuf1[i]='0'; ? ? ? ? ? } ? ? ? ? ? for(i=len_num2-len_num1+1;i<=len_num2;i++) ? ? ? ? ? { ? ? ? ? ? ? ? tempbuf1[i]=str_num1[i-(len_num2-len_num1+1)]; ? ? ? ? ? } ? ? ? ? ? for(i=1;i<=len_num2;i++) ? ? ? ? ? { ? ? ? ? ? ? ? tempbuf2[i]=str_num2[i-1]; ? ? ? ? ? } ? ? ? } ? //------------------------------------------ ? ? ? else if(len_num2<len_num1) ? ? ? { ? ? ? ? ? for(i=1;i<=(len_num1-len_num2);i++) ? ? ? ? ? { ? ? ? ? ? ? ? tempbuf2[i]='0'; ? ? ? ? ? } ? ? ? ? ? for(i=len_num1-len_num2+1;i<=len_num1;i++) ? ? ? ? ? { ? ? ? ? ? ? ? tempbuf2[i]=str_num2[i-(len_num1-len_num2+1)]; ? ? ? ? ? } ? ? ? ? ? for(i=1;i<=len_num1;i++) ? ? ? ? ? { ? ? ? ? ? ? ? tempbuf1[i]=str_num1[i-1]; ? ? ? ? ? } ? ? ? } ? /* ?printf("-------------------------------------------------------n转换之后:n");? ? ? printf("str1=: %sn",tempbuf1);? ? ? printf("str2=: %sn",tempbuf2);? //==========================================================================? ? ? int len=0;? ? ? i=0;? ? ? while(tempbuf1[i]!='')? ? ? {? ? ? ? ? len++;? ? ? ? ? i++;? ? ? }? */ ? } ? ?? //========================================================================================= =========== ? //extern "C" __declspec(dllexport) ? void ?Chen(char *tempbuf1,? char *tempbuf2,char *result) ? { ? ? ? char buf1[200]={0}; ? ? ? char buf2[200]={0}; ? ? ? Trans(tempbuf1,tempbuf2,buf1,buf2); ? ? ? strcpy(tempbuf1,buf1); ? ? ? strcpy(tempbuf2,buf2); ? ?? ? ? int len=0; ? ? ? int i=0; ? ? ? int j=0; ? ? ? int n=0; ? ? ? int jinwei=0; ? ?? ? ? while(tempbuf1[i]!='') ?// 字符串长度 ? ? ? { ? ? ? ? ? len++; ? ? ? ? ? i++; ? ? ? } ? ?? ? ? int temp[200]={0}; ? ? ? ?// 该数组用来存储各次方的系数 10为底 ? ?? ? ? int max=2*len; ? ? ? for(i=0;i<=max;i++) ? ? ? { ? ? ? ? ? if(i<len-1) ? ? ? ? ? { ? ? ? ? ? ? ? for(j=0;(j<=len-1)&&(j>=i+1-len)&&(j<=i);j++) ? ? ? ? ? ? ? { ? ? ? ? ? ?? ? ? ? ? ? ? ? ? temp[i]+=((int)tempbuf1[len-1-j]-48 )*((int)tempbuf2[len-1-i+j]-48); ? ? ? ? ? ? ?? ? ? ? ? ? ? } ? ? ? ? ? ? ? temp[i]+=jinwei; ? ? ? ? ? ? ? ?? ? ? ? ? ? ? if (temp[i]>=10) ?//&&temp[i]<100 ? ? ? ? ? ? ? { ? ? ? ? ? ? ? ? ? jinwei ?= temp[i]/10; ? ? ? ? ? ? ? ? ? temp[i] = temp[i]%10; ? ? ? ? ? ? ? } ? ? ? ? ? ? ? else jinwei=0; ? ? ? ? ? } ? ? ? ? ? else if (i>=len-1) ? ? ? ? ? { ? ? ? ? ? ? ? for(j=i-len+2;(j<=len-1)&&(j>=i+1-len)&&(j<=i);j++) ? ? ? ? ? ? ? { ? ? ? ? ? ?? ? ? ? ? ? ? ? ? temp[i]+=((int)tempbuf1[len-1-j]-48 )*((int)tempbuf2[len-1-i+j]-48); ? ? ? ? ? ? ?? ? ? ? ? ? ? } ? ? ? ? ? ? ? temp[i]+=jinwei; ? ? ? ? ? ? ? ?? ? ? ? ? ? ? if (temp[i]>=10) ?//&&temp[i]<100 ? ? ? ? ? ? ? { ? ? ? ? ? ? ? ? ? jinwei ?= temp[i]/10; ? ? ? ? ? ? ? ? ? temp[i] = temp[i]%10; ? ? ? ? ? ? ? } ? ? ? ? ? ? ? else jinwei=0; ? ? ? ? ? } ? ? ? ? ? ?? ? ? ? ? ?? ? ? } ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?? ? ? i=max; ? ? ? while(i>=0) ? ? ? { ? ? ? ? ? if(temp[i]!=0) ? ? ? ? ? ? ? break; ? ? ? ? ? ?? ? ? ? ? else i--; ? ? ? ? ? ? ? ?? ? ? } ? // ?printf("str1 * str2= "); ? ? ? int num=i; // 科学计数法次数 ? /* ?for (j=num;j>=0;j--)? ? ? {? ? ? ? ? ? ? ? ? ? printf("%d",temp[j]);? ? ? ? ? ? ? ? } ? ? ? ? printf("t%dn",num);? */ ? //=========================================== ? ? ? memset(result,200); ? ? ? for(i=num;i>=0;i--) ? ? ? { ? ? ? ? ? result[num-i]=(char)(temp[i]+48); ? ?? ? ? } ? // ?printf("result = %st%dn",strlen(result)-1); ? ? ? ?? } ? ?? ?? //===================================================================== ? ?? //extern "C" __declspec(dllexport) ?? void ?ChenFang(char *tempbuf1,? int num,char *result) ? ?// 大数乘方 ? { ? ? ? static int count=0; ? ?? ? ? if (count==0) ? ? ? { ? ? ? ? ? char c[200]={'1',0}; ? ? ? ? ? memset(result,200); ? ? ? ? ? Chen( tempbuf1,c,result); ? ?? ? ? } ? ? ? count++; ? ?? ? ? if(num==1) ? ? ? { ? ? ? ? ? return ; ? ? ? } ? ?? ? ? else if(num>1) ? ? ? { ? ? ? ? ? Chen( tempbuf1,result); ? ? ? ? ? num--; ? ? ? ? ? ChenFang(tempbuf1,result); ? ? ? } ? ?? ? ? ThrowAway_0 (result); ? ?? } ? ?? ?? //================================================= ? ?? //extern "C" __declspec(dllexport) ?? void ?Mod(char *tempbuf1,char *result,char *div_result ) ? ?// 大数求 余 ? { ? ? ? memset(result,200); ? ? ? ThrowAway_0 (tempbuf1); ? ? ? ThrowAway_0 (tempbuf2); ? ? ? int max1 = strlen(tempbuf1); ? ? ? int max2 = strlen(tempbuf2); ? ? ? int temp=0; ? ? ? int ret=1; ? ? ? char div[200]={0}; ? ? ? int ?count=0; ? ? ? char AfterSub[200]={0}; ? ?? ? ? if(max1-max2>=0) ? ? ? ? ? ? ? { ? ? ? ? ? while(ret==1) ? ? ? ? ? ? { ? ? ? ? ? ? ? ? ? ? ? ? ? ?? ? ? ? ? ? ? Sub (tempbuf1,AfterSub); ? //减去之后在判断是否得到结果----即： tempbuf1 < tempbuf2 ? ? ? ? ? ? ? memset(tempbuf1,200); ? ? ? ? ? ? ? strcpy (tempbuf1,AfterSub); ?? ? ? ? ? ? ? ret=Compare(tempbuf1,tempbuf2); ? ? ? ? ? ? ? ThrowAway_0 (tempbuf1); ? ? ? ? ? ? ? ThrowAway_0 (tempbuf2); ?? ? ? ? ? ? ? count++; ? ? ? ? ? ? ? if(count>=1000) ? ? ? ? ? ? ? { ? ? ? ? ? ? ? ? ? itoa(count,div,10); ? ? ? ? ? ? ? ? ? Add(div_result,div_result); ? ? ? ? ? ? ? ? ? count=0; ? ? ? ? ? ? ? } ? ? ? ? ? ? ? ?? ? ? ? ? } ? ? ? ? ? ? ? ?? ? ? ? ? strcpy (result,tempbuf1); ? ? ? ? ? ThrowAway_0 (result); ? ? ? ? ? itoa(count,10); ? ? ? ? ? Add(div_result,div_result); ? ? ? } ? ?? ? ? else ?// if (Compare(tempbuf1,tempbuf2)==-1) // 后者大的话，mod 就是 tempbuf1 ? ? ? { ? ? ? ? ? strcpy (result,tempbuf1); ? ? ? ? ? ThrowAway_0 (result); ? ? ? } ? } ? ?? //======================================================== ? ?? //extern "C" __declspec(dllexport) ?? void Add(char *tempbuf1,char *result) // 大数相加 result = tempbuf1 +? tempbuf2 ? { ? ? ? char buf1[200]={0}; ? ? ? char buf2[200]={0}; ? ? ? Trans(tempbuf1,buf2); ? ?? ? ? int i=0; ? ? ? int temp=0; ? ? ? int jinwei=0; ? ? ? int len=0; ? ? ? while(tempbuf1[i]!='') ? ? ? { ? ? ? ? ? len++; ? ? ? ? ? i++; ? ? ? } ? ? ? for(i=len-1;i>=0;i--) ? ? ? { ? ?? ? ? ? ? temp=(int)(tempbuf1[i]+tempbuf2[i]+jinwei-96); ? ? ? ? ? if(temp>=10) ? ? ? ? ? { ? ? ? ? ? ? ? temp=temp-10; ? ? ? ? ? ? ? jinwei=1; ? ? ? ? ? } ? ? ? ? ? else jinwei=0; ? ? ? ? ? result[i]=(char)(temp+48); ? ? ? ? } ? ? ? ThrowAway_0 (result); ? ?? } ? ?? //================================================ ? ?? void ?ThrowAway_0 (char *tempbuf ) ?// 去除结果前面的 连续的无意义的 "0" ? { ? ? ? char buf[200]={0}; ? ?? ? ? int n = strlen(tempbuf)-1; ? ? ? int i=0; ? ? ? while(i<n) ? ? ? { ? ? ? ? ? if (tempbuf[i]!='0') ? ? ? ? ? { ? ? ? ? ? ? ? break; ? ? ? ? ? } ? ? ? ? ? else ?? ? ? ? ? { ? ? ? ? ? ? ? i++; ? ? ? ? ? } ? ? ? } ? ? ? int Throw = i; ? ? ? for (i=0;i<=n-Throw;i++) ? ? ? { ? ? ? ? ? buf[i]=tempbuf[i+Throw]; ? ? ? } ? ?? ? ? strcpy(tempbuf,buf); ? ?? } ? ?? //======================================================= ? ?? //extern "C" __declspec(dllexport) ?? void Sub(char *tempbuf1,char *tempbuf2,char *result) ? ? // 大数相减 result = tempbuf1 -? tempbuf2 ? { ? ? ? ThrowAway_0 (tempbuf1); ? ? ? ThrowAway_0 (tempbuf2); ? ? ? memset(result,200); ? ? ? char buf1[200]={0}; ? ? ? char buf2[200]={0}; ? ? ? Trans(tempbuf1,buf2); ? ? ? memset(tempbuf1,200); ? ? ? memset(tempbuf2,200); ? ?? ? ? strcpy(tempbuf1,buf2); ? ?? ? ? int i=0; ? ? ? int temp=0; ? ? ? int jiewei=0; ? ? ? int len=0; ? ? ? int ret=1; ? ?? ? ? while(tempbuf1[i]!='') ? ?// tempbuf1 和 tempbuf2 的长度相等 ? ? ? { ? ? ? ? ? len++; ? ? ? ? ? i++; ? ? ? } ? ? ? ?? ? ? ret = Compare(tempbuf1,tempbuf2); ? ? ? if(ret==1) ? ? ? { ? ? ? ? ? for(i=len-1;i>=0;i--) ? ? ? ? ? { ? ? ? ? ? ? ? ? temp = (int)tempbuf1[i] - (int)tempbuf2[i] - jiewei; ? ? ? ? ? ? ? if (temp>=0) ? ? ? ? ? ? ? { ? ? ? ? ? ? ? ? ? result[i]=(char)(temp+48); ? ? ? ? ? ? ? ? ? jiewei=0; ? ? ? ? ? ? ? } ? ? ? ? ? ? ? else if (temp<0) ? ? ? ? ? ? ? { ? ? ? ? ? ? ? ? ? result[i]=(char)(temp+10+48); ? ? ? ? ? ? ? ? ? jiewei=1; ? ? ? ? ? ? ? } ? ? ? ? ? } ? ? ? ? ? ThrowAway_0 (result); ? ? ? } ? ? ? else if(ret==0) ?? ? ? { ? ? ? ? ? memset(result,200); ? ? ? ? ? result[0]='0'; ? ? ? } ? ? ? else if(ret==-1) ?? ? ? { ? ? ? ? ? for(i=len-1;i>=0;i--) ? ? ? ? ? { ? ? ? ? ? ? ? ? temp = (int)tempbuf2[i] - (int)tempbuf1[i] - jiewei; ? ? ? ? ? ? ? if (temp>=0) ? ? ? ? ? ? ? { ? ? ? ? ? ? ? ? ? result[i]=(char)(temp+48); ? ? ? ? ? ? ? ? ? jiewei=0; ? ? ? ? ? ? ? } ? ? ? ? ? ? ? else if (temp<0) ? ? ? ? ? ? ? { ? ? ? ? ? ? ? ? ? result[i]=(char)(temp+10+48); ? ? ? ? ? ? ? ? ? jiewei=1; ? ? ? ? ? ? ? } ? ? ? ? ? } ? ? ? ? ? ThrowAway_0 (result); ? ? ? ? ? memset(buf1,200); ? ? ? ? ? sprintf(buf1,"-%s",result); ?? ? ? ? ? strcpy(result,buf1); ? ? ? } ? ? ? ?? } ? ?? //====================================================================== ? ?? //=================================================================================== ? ?? int Compare(char *tempbuf1,char *tempbuf2) ? { ? ? ? ThrowAway_0 (tempbuf1); ? ? ? ThrowAway_0 (tempbuf2); ? ? ? char buf1[200]={0}; ? ? ? char buf2[200]={0}; ? ? ? Trans(tempbuf1,200); ? ? ? strcpy(tempbuf1,buf2); ? ? ? ?? ?? ? ? int ret=1; ? ? ? int count=0; ? ? ? ? ? while(count<200) ? ? ? { ? ? ? ? ? ?? ? ? ? ? int m=(int)tempbuf1[count]-48; ? ? ? ? ? int n=(int)tempbuf2[count]-48; ? ? ? ? ? if(m==n) ? ? ? ? ? { ? ? ? ? ? ? ? count++; ? ? ? ? ? ? ? ret=0; ? ? ? ? ? } ? ? ? ? ? else if(m>n) ? ? ? ? ? { ? ? ? ? ? // ?printf("tempbuf1>tempbuf2n"); ? ? ? ? ? ? ? ret=1; ? ? ? ? ? ? ? break; ? ? ? ? ? ? ? ?? ? ? ? ? } ? ? ? ? ? else if(m<n) ? ? ? ? ? { ? ? ? ? ? // ?printf("tempbuf1<tempbuf2n"); ? ? ? ? ? ? ? ret=-1; ? ? ? ? ? ? ? break; ? ? ? ? ? } ? ? ? } ? ? ? return ret; ? } ? 2.头文件 [cpp] view plain copy #include <stdio.h> ?? #include <string.h> ? #include <stdlib.h> ? ?? //extern "C" __declspec(dllexport) ?? void Add (char *tempbuf1,char *result); ?// 大数相加 result = tempbuf1 +? tempbuf2 ? ?? //extern "C" __declspec(dllexport) ?? void Sub(char *tempbuf1,char *result); ? ? // 大数相减 result = tempbuf1 -? tempbuf2 (默认前者大) ? ?? int Compare(char *tempbuf1,char *tempbuf2); ? ?? //extern "C" __declspec(dllexport) ? ? void Chen (char *tempbuf1,char *result); // 大数相乘 result = tempbuf1 *? tempbuf2 ? ?? void Trans (char *str_num1,? char *str_num2,char *tempbuf2); // 大数转化 为等长，且最前面添加 “0” ? ?? ?? //extern "C" __declspec(dllexport) ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?// 大数除法? div_result = tempbuf1 / tempbuf2 ? void ?Mod(char *tempbuf1,char *div_result ); ? ?// 大数求 余 result = tempbuf1 % tempbuf2 ? ? ?? ?? //extern "C" __declspec(dllexport) ?? void ?ChenFang (char *tempbuf1,char *result); ? ?// 大数乘方 ? result =? tempbuf1 ^ num ? ?? void ?ThrowAway_0 (char *tempbuf ); ?// 去除结果前面的 连续的无意义的 "0" ? ========