BigNumCalculator
好久没写日志了，
其实草稿写了不少，
无论是感悟还是技术，
到最后却总是删掉，
总是觉得少了一股劲，
总越来越觉得自己还不够透彻。。。
上周时候，偶然同事讨论起了大数运算，
想起来自己以前做过一个BigNumCalculator 的阶乘计算器，
代码半天没找着，
打回家让老妈开了qq远程控制家里的电脑，
良久，才在一个很深的目录找到了代码，
忽然自己惊讶了，目录里还有很多自己过往写的代码，
从Web，各种DCC工具的练习，C++， DX。。。
好久未有这般感慨了，往日的一幕幕，浮现在眼前，
大学伊始，不准带电脑，
试过图书馆阅览室站一个早上看技术书，
然后破解一台角落处查书用的电脑，
然后站一个下午用记事本写代码，
像JS之类的就直接用IE测试，
抑或厚着脸皮买了鞋套，偷偷混进机房，
搞得几位机房的老师都认得我了；
有半年断网，
咬咬牙，买了一个NB的8G优盘，
然后在图书馆机房充值，
每周末去一次性准备好自己想看的资料，
然后回寝室慢慢啃。
现在想想挺酸楚的，
可仔细回想当时充满勇气和努力的自己，
我100%确定我当时是开心的，毫无保留地付出的，
反观今日的自己，物质上有了很大的充实，
同时也背负起了一些沉重遥远的现实目标，更多的时候，自己开始计较了。。。
这一切是我想要的么？
庆幸，
我现在还能找到当时的代码，
以及尚存的勇气和热血，
毋须多言，
翻开青涩的代码，
留下一个有趣的注释标头，
卯足劲，
一切从头开始，
明天还有好远，
继续，努力！
1
/*
2 091111
3 大数阶乘运算
4 节日快乐
5 -----Zephyr
6 */
7 /*
8 基本思路，化而治之
9 123456789
10 * 211
11 ---------
12 123456789
13 +1234567890
14 +24691357800
15 ------------
16 =26049382479
17 */
18 #include <iostream>
19 #include <stdio.h>
20 #include <windows.h>
21 #include <assert.h>
22 #include <conio.h>
23
24 #pragma comment(lib,"winmm.lib")
25 using namespace std;
26 const int MaxLength= 128;
27 const int Max=MaxLength- 1;
28
29 void Copy( int *des, int *src)
30 {
31 for( int j=Max;j>= 0;j--)
32 {
33 des[j]=src[j];
34 }
35 }
36
37 void Zero( int *src)
38 {
39 for( int j= 0;j<MaxLength;j++)
40 {
41 src[j]= 0;
42 }
43 }
44
45 int GetLength( int *a)
46 {
47 int i=Max;
48 while( a[i] <= 0 && i> 0 )
49 --i;
50 return ++i;
51 }
52
53 void Increase( int *a1)
54 {
55 a1[Max]++;
56 for( int j=Max;j> 0;j--)
57 {
58 if(a1[j]> 10)
59 {
60 a1[j]-= 10;
61 a1[j- 1]++;
62 }
63 else
64 break;
65 }
66 }
67
68 // 相加的时候向前提了k个位置
69 void Sum( int *sum, int * a, int k= 0)
70 {
71 int j2;
72 // 有漏洞，但可以把Max设大
73 for( int j= 0 ;j < GetLength(a) ;++j)
74 {
75 j2= j+k;
76 sum[ j2 ]+=a[j];
77 // 加法最大值应该只有 8+8 == 19
78 if(sum[ j2 ] >= 10)
79 {
80 if( j2 >= MaxLength )
81 assert( " Over flow! ");
82 if( sum[ j2 ]< 0 || sum[ j2 ] > 19 )
83 assert( " Oh no ");
84 sum[j2]-= 10;
85 if(sum[j2]>= 10)
86 sum[j2]= 0;
87 sum[j2+ 1]++;
88 }
89 }
90 }
91
92 void Multi( int *sum, int *aO, int k) // 十进制乘法
93 {
94 if( k< 0 || k> 9 )
95 assert( " Wrong Number ");
96 // a0*k送入结果sum中
97 int t= 0;
98 // 不能改变a0
99 Copy(sum,aO);
100 int j;
101 for( j= 0 ; j < GetLength(aO) ; ++j )
102 {
103 sum[j]*=k;
104 if(t)
105 {
106 sum[j]+=t;
107 t= 0;
108 }
109 // 进位
110 if(sum[j] >= 10)
111 {
112 t=sum[j]/ 10;
113 sum[j]=sum[j]% 10;
114 }
115 }
116 // 最后一位还可能会进位
117 if(t)
118 {
119 sum[j]+=t;
120 }
121 }
122
123 void MultiX( int *total, int *t) // 大数相乘
124 {
125 int temp[MaxLength];
126 int tempSum[MaxLength];
127 Zero( tempSum );
128 for( int j= 0; j<GetLength( t ) ;++j)
129 {
130 Multi(temp,total,t[j]);
131 Sum(tempSum,temp,j);
132 }
133 Copy(total,tempSum);
134 }
135
136 void AnotherWay( int n)
137 {
138 // 精确到12!看来你也不行哈。。。-_-!
139 // 寻址问题
140 SIZE_T x= 1;
141 int i= 1;
142 while(i<=n)
143 {
144 x*=i;
145 i++;
146 }
147 cout<< " K.O: "<<endl<<x<<endl;
148 /* x=479001600;
149 x*=13;
150 cout<<x; */
151 }
152
153 void main0()
154 {
155 int total[MaxLength],t[MaxLength];
156 int n;
157 float t1= 0;
158 bool debug= 0;
159 cout<< " 输入阶乘级数： "<<endl;
160 cin>>n;
161 while(n> 0)
162 {
163 t1=( float)timeGetTime();
164 Zero(total);Zero(t);
165 total[Max]=t[Max]= 1;
166 AnotherWay(n);
167 n--;
168 int i= 2;
169 while(n> 0) // 阶乘，从1*2*3*...*n
170 {
171 Increase(t);
172 MultiX(total,t);
173 if(debug)
174 {
175 cout<< " 第 "<<i++<< " 层: "<<endl; // DeBug使用，i既是当前乘数
176 for( int j=MaxLength-GetLength(total);j<MaxLength;j++)
177 cout<<total[j];
178 cout<<endl;
179 }
180 n--;
181 }
182 cout<< " 结果为： "<<endl;
183 for( int j=MaxLength-GetLength(total);j<MaxLength;j++) // 从不为零的数位开始输出结果
184 cout<<total[j];
185 cout<<endl<< " 用时: "<<( float)timeGetTime()-t1<< " ms "<<endl;
186 cout<< " -----------Over-------------- "<<endl;
187 cout<< " 输入阶乘级数(0结束)： "<<endl;
188 cin>>n;
189 }
190 }
191
192 void GetABigNumber( int InVec[] )
193 {
194 int n;
195 int i= 0;
196 char c;
197 // cin>>n;
198 // cin.get(c);
199 c=getch();
200 n=atoi(&c);
201 while ( n>= 0 && n < 10 )
202 {
203 if( i > 0)
204 {
205 int i1=i;
206 while( i1 > 0)
207 {
208 InVec[ i1 ]=InVec[i1- 1];
209 --i1;
210 }
211 }
212 InVec[ 0]=n;
213 ++i;
214 c=getch();
215 n=atoi(&c);
216 }
217 }
218 void PraseString2BigNumber( int inVec[], char Srt[])
219 {
220 int length = strlen( Srt );
221 Zero(inVec);
222 char ch;
223 for ( int i=length- 1, j= 0;j<length;++j)
224 {
225 ch=Srt[i--];
226 inVec[ j ]=atoi( &ch );
227 }
228 }
229 void main()
230 {
231 unsigned int p1= 10;
232 unsigned int p2= 0xa0000000;
233 float pw=p2/p1;
234 int total[MaxLength],t[MaxLength];
235 int n;
236 float t1= 0;
237 bool debug= 0;
238 /*
239 Zero(total);
240 Zero(t);
241 cout<<"输入第一个大数："<<endl;
242 GetABigNumber(total);
243 cout<<"输入第二个大数："<<endl;
244 GetABigNumber(t);
245 */
246 // char str1[]="12888";
247 // char str2[]="12888";
248 // 16384
249 char str1[]= " 145682565648569854785693568458967890 ";
250 char str2[]= " 64318641564154513745146843514 ";
251 PraseString2BigNumber(total,str1);
252 PraseString2BigNumber(t,str2);
253 cout<< " 第一个大数： "<<str1<<endl;
254 cout<< " 乘以第二个大数： "<<str2<<endl;
255 t1=( float)timeGetTime();
256 MultiX(total,t);
257 cout<< " 结果为： "<<endl;
258 for( int j=GetLength(total)- 1 ;j>= 0;--j) // 从不为零的数位开始输出结果
259 cout<<total[j];
260 cout<<endl<< " 用时: "<<( float)timeGetTime()-t1<< " ms "<<endl;
261 cout<< " -----------Over-------------- "<<endl;
262 cin>>n;
263 }
2 091111
3 大数阶乘运算
4 节日快乐
5 -----Zephyr
6 */
7 /*
8 基本思路，化而治之
9 123456789
10 * 211
11 ---------
12 123456789
13 +1234567890
14 +24691357800
15 ------------
16 =26049382479
17 */
18 #include <iostream>
19 #include <stdio.h>
20 #include <windows.h>
21 #include <assert.h>
22 #include <conio.h>
23
24 #pragma comment(lib,"winmm.lib")
25 using namespace std;
26 const int MaxLength= 128;
27 const int Max=MaxLength- 1;
28
29 void Copy( int *des, int *src)
30 {
31 for( int j=Max;j>= 0;j--)
32 {
33 des[j]=src[j];
34 }
35 }
36
37 void Zero( int *src)
38 {
39 for( int j= 0;j<MaxLength;j++)
40 {
41 src[j]= 0;
42 }
43 }
44
45 int GetLength( int *a)
46 {
47 int i=Max;
48 while( a[i] <= 0 && i> 0 )
49 --i;
50 return ++i;
51 }
52
53 void Increase( int *a1)
54 {
55 a1[Max]++;
56 for( int j=Max;j> 0;j--)
57 {
58 if(a1[j]> 10)
59 {
60 a1[j]-= 10;
61 a1[j- 1]++;
62 }
63 else
64 break;
65 }
66 }
67
68 // 相加的时候向前提了k个位置
69 void Sum( int *sum, int * a, int k= 0)
70 {
71 int j2;
72 // 有漏洞，但可以把Max设大
73 for( int j= 0 ;j < GetLength(a) ;++j)
74 {
75 j2= j+k;
76 sum[ j2 ]+=a[j];
77 // 加法最大值应该只有 8+8 == 19
78 if(sum[ j2 ] >= 10)
79 {
80 if( j2 >= MaxLength )
81 assert( " Over flow! ");
82 if( sum[ j2 ]< 0 || sum[ j2 ] > 19 )
83 assert( " Oh no ");
84 sum[j2]-= 10;
85 if(sum[j2]>= 10)
86 sum[j2]= 0;
87 sum[j2+ 1]++;
88 }
89 }
90 }
91
92 void Multi( int *sum, int *aO, int k) // 十进制乘法
93 {
94 if( k< 0 || k> 9 )
95 assert( " Wrong Number ");
96 // a0*k送入结果sum中
97 int t= 0;
98 // 不能改变a0
99 Copy(sum,aO);
100 int j;
101 for( j= 0 ; j < GetLength(aO) ; ++j )
102 {
103 sum[j]*=k;
104 if(t)
105 {
106 sum[j]+=t;
107 t= 0;
108 }
109 // 进位
110 if(sum[j] >= 10)
111 {
112 t=sum[j]/ 10;
113 sum[j]=sum[j]% 10;
114 }
115 }
116 // 最后一位还可能会进位
117 if(t)
118 {
119 sum[j]+=t;
120 }
121 }
122
123 void MultiX( int *total, int *t) // 大数相乘
124 {
125 int temp[MaxLength];
126 int tempSum[MaxLength];
127 Zero( tempSum );
128 for( int j= 0; j<GetLength( t ) ;++j)
129 {
130 Multi(temp,total,t[j]);
131 Sum(tempSum,temp,j);
132 }
133 Copy(total,tempSum);
134 }
135
136 void AnotherWay( int n)
137 {
138 // 精确到12!看来你也不行哈。。。-_-!
139 // 寻址问题
140 SIZE_T x= 1;
141 int i= 1;
142 while(i<=n)
143 {
144 x*=i;
145 i++;
146 }
147 cout<< " K.O: "<<endl<<x<<endl;
148 /* x=479001600;
149 x*=13;
150 cout<<x; */
151 }
152
153 void main0()
154 {
155 int total[MaxLength],t[MaxLength];
156 int n;
157 float t1= 0;
158 bool debug= 0;
159 cout<< " 输入阶乘级数： "<<endl;
160 cin>>n;
161 while(n> 0)
162 {
163 t1=( float)timeGetTime();
164 Zero(total);Zero(t);
165 total[Max]=t[Max]= 1;
166 AnotherWay(n);
167 n--;
168 int i= 2;
169 while(n> 0) // 阶乘，从1*2*3*...*n
170 {
171 Increase(t);
172 MultiX(total,t);
173 if(debug)
174 {
175 cout<< " 第 "<<i++<< " 层: "<<endl; // DeBug使用，i既是当前乘数
176 for( int j=MaxLength-GetLength(total);j<MaxLength;j++)
177 cout<<total[j];
178 cout<<endl;
179 }
180 n--;
181 }
182 cout<< " 结果为： "<<endl;
183 for( int j=MaxLength-GetLength(total);j<MaxLength;j++) // 从不为零的数位开始输出结果
184 cout<<total[j];
185 cout<<endl<< " 用时: "<<( float)timeGetTime()-t1<< " ms "<<endl;
186 cout<< " -----------Over-------------- "<<endl;
187 cout<< " 输入阶乘级数(0结束)： "<<endl;
188 cin>>n;
189 }
190 }
191
192 void GetABigNumber( int InVec[] )
193 {
194 int n;
195 int i= 0;
196 char c;
197 // cin>>n;
198 // cin.get(c);
199 c=getch();
200 n=atoi(&c);
201 while ( n>= 0 && n < 10 )
202 {
203 if( i > 0)
204 {
205 int i1=i;
206 while( i1 > 0)
207 {
208 InVec[ i1 ]=InVec[i1- 1];
209 --i1;
210 }
211 }
212 InVec[ 0]=n;
213 ++i;
214 c=getch();
215 n=atoi(&c);
216 }
217 }
218 void PraseString2BigNumber( int inVec[], char Srt[])
219 {
220 int length = strlen( Srt );
221 Zero(inVec);
222 char ch;
223 for ( int i=length- 1, j= 0;j<length;++j)
224 {
225 ch=Srt[i--];
226 inVec[ j ]=atoi( &ch );
227 }
228 }
229 void main()
230 {
231 unsigned int p1= 10;
232 unsigned int p2= 0xa0000000;
233 float pw=p2/p1;
234 int total[MaxLength],t[MaxLength];
235 int n;
236 float t1= 0;
237 bool debug= 0;
238 /*
239 Zero(total);
240 Zero(t);
241 cout<<"输入第一个大数："<<endl;
242 GetABigNumber(total);
243 cout<<"输入第二个大数："<<endl;
244 GetABigNumber(t);
245 */
246 // char str1[]="12888";
247 // char str2[]="12888";
248 // 16384
249 char str1[]= " 145682565648569854785693568458967890 ";
250 char str2[]= " 64318641564154513745146843514 ";
251 PraseString2BigNumber(total,str1);
252 PraseString2BigNumber(t,str2);
253 cout<< " 第一个大数： "<<str1<<endl;
254 cout<< " 乘以第二个大数： "<<str2<<endl;
255 t1=( float)timeGetTime();
256 MultiX(total,t);
257 cout<< " 结果为： "<<endl;
258 for( int j=GetLength(total)- 1 ;j>= 0;--j) // 从不为零的数位开始输出结果
259 cout<<total[j];
260 cout<<endl<< " 用时: "<<( float)timeGetTime()-t1<< " ms "<<endl;
261 cout<< " -----------Over-------------- "<<endl;
262 cin>>n;
263 }