高精度四则运算模板

转载自：http://download.csdn.net/download/yizhenmi/3046636

[cpp]  view plain  copy

 

1.                             大数模版  
2. /*大数加法*/   
3. # include<stdio.h>  
4. # include<string.h>  
5. # include<malloc.h>  
6.   
7. void add(char* a,char* b,char* c)  
8. {  
9.     int i,j,k,max,min,n,temp;  
10.     char *s,*pmax,*pmin;  
11.     max=strlen(a);  
12.     min=strlen(b);  
13.     if (max<min)  
14.     {  
15.         temp=max;  
16.         max=min;  
17.         min=temp;  
18.         pmax=b;  
19.         pmin=a;  
20.     }  
21.     else  
22.     {  
23.         pmax=a;  
24.         pmin=b;  
25.     }  
26.     s=(char*)malloc(sizeof(char)*(max+1));  
27.     s[0]='0';  
28.     for (i=min-1,j=max-1,k=max;i>=0;i--,j--,k--)  
29.         s[k]=pmin[i]-'0'+pmax[j];  
30.        for (;j>=0;j--,k--)  
31.            s[k]=pmax[j];  
32.     for (i=max;i>=0;i--)  
33.         if (s[i]>'9')  
34.         {  
35.             s[i]-=10;  
36.             s[i-1]++;  
37.         }  
38.     if (s[0]=='0')  
39.     {  
40.         for (i=0;i<=max;i++)  
41.             c[i-1]=s[i];  
42.            c[i-1]='\0';  
43.     }  
44.     else  
45.     {  
46.         for (i=0;i<=max;i++)  
47.             c[i]=s[i];  
48.            c[i]='\0';  
49.     }  
50.     free(s);  
51. }  
52.   
53. /*大数凿法*/   
54. void subtract(char* a,char* b,char* c)  
55. {  
56.     int i,j,ca,cb;  
57.     ca=strlen(a);  
58.     cb=strlen(b);  
59.     if (ca>cb||(ca==cb&&strcmp(a,b)>=0))  
60.     {  
61.         for (i=ca-1,j=cb-1;j>=0;i--,j--)  
62.             a[i]-=(b[j]-'0');  
63.            for (i=ca-1;i>=0;i--)  
64.                if (a[i]<'0')  
65.                {  
66.                    a[i]+=10;  
67.                    a[i-1]--;  
68.             }  
69.         i=0;  
70.         while (a[i]=='0')  
71.             i++;  
72.            if (a[i]=='\0')  
73.           {  
74.               c[0]='0';  
75.               c[1]='\0';  
76.         }  
77.         else  
78.         {  
79.             for (j=0;a[i]!='\0';i++,j++)  
80.                 c[j]=a[i];  
81.                c[j]='\0';  
82.         }  
83.     }  
84.     else  
85.     {  
86.         for (i=ca-1,j=cb-1;i>=0;i--,j--)  
87.             b[j]-=(a[i]-'0');  
88.            for (j=cb-1;j>=0;j--)  
89.                if (b[j]<'0')  
90.                {  
91.                    b[j]+=10;  
92.                    b[j-1]--;  
93.             }  
94.         j=0;  
95.         while (b[j]=='0')  
96.             j++;  
97.            i=1;  
98.            c[0]='-';  
99.            for (;b[j]!='\0';i++,j++)  
100.                c[i]=b[j];  
101.           c[i]='\0';  
102.     }              
103. }  
104.   
105. /* 大数乘法*/   
106. void multiply(char* a,char* b,char* c)  
107. {  
108.     int i,j,ca,cb,* s;  
109.     ca=strlen(a);  
110.     cb=strlen(b);  
111.     s=(int*)malloc(sizeof(int)*(ca+cb));  
112.     for (i=0;i<ca+cb;i++)  
113.         s[i]=0;  
114.     for (i=0;i<ca;i++)  
115.         for (j=0;j<cb;j++)  
116.             s[i+j+1]+=(a[i]-'0')*(b[j]-'0');  
117.     for (i=ca+cb-1;i>=0;i--)  
118.         if (s[i]>=10)  
119.         {  
120.             s[i-1]+=s[i]/10;  
121.             s[i]%=10;  
122.         }  
123.     i=0;  
124.     while (s[i]==0)  
125.         i++;  
126.        for (j=0;i<ca+cb;i++,j++)  
127.            c[j]=s[i]+'0';  
128.     c[j]='\0';  
129.     free(s);  
130. }  
131.   
132. /*大数除法,返回余数*/  
133. int dividor(char* a,int b,char* c)  
134. {  
135.     int i,j,temp=0,n;  
136.     char* s;  
137.     n=strlen(a);  
138.     s=(char*)malloc(sizeof(char)*(n+1));  
139.     for (i=0;a[i]!=0;i++)  
140.        {  
141.            temp=temp*10+a[i]-'0';  
142.            s[i]=temp/b+'0';  
143.            temp%=b;  
144.     }  
145.     s[i]='\0';  
146.     for (i=0;s[i]=='0'&&s[i]!='\0';i++);  
147.     if (s[i]=='\0')  
148.     {  
149.         c[0]='0';  
150.         c[1]='\0';  
151.     }  
152.     else  
153.     {  
154.         for (j=0;s[i]!='\0';i++,j++)  
155.             c[j]=s[i];  
156.            c[j]='\0';  
157.     }  
158.     free(s);  
159.     return temp;  
160. }