第25講 n進数の演算---その1 加法
第8話 プログラム解説その2
プログラム主要再掲
int g(int* a,int* b,int* c,int n){
int i,j,k,d[10000],e[10000];
for(i=0;;i++){
if(a[i]==n){
j=i;
break;
}
}
for(i=0;;i++){
if(b[i]==n){
k=i;
break;
}
}
if(j>=k){
for(i=0;i<=j+2;i++){
c[i]=0;
d[i]=0;
e[i]=0;
}
for(i=0;i<j;i++)d[i]=a[i];
for(i=0;i<k;i++)e[i]=b[i];
for(i=0;i<j;i++){
c[i]+=(d[i]+e[i]);
c[i+1]+=c[i]/n;
c[i]=c[i]%n;
}
for(i=j+2;;i--){
if(c[i]>0){
c[i+1]=n;
return(i+1);
}
}
}
else{
for(i=0;i<=k+2;i++){
c[i]=0;
d[i]=0;
e[i]=0;
}
for(i=0;i<k;i++)d[i]=b[i];
for(i=0;i<j;i++)e[i]=a[i];
for(i=0;i<k;i++){
c[i]+=(d[i]+e[i]);
c[i+1]+=c[i]/n;
c[i]=c[i]%n;
}
for(i=k+2;;i--){
if(c[i]>0){
c[i+1]=n;
return(i+1);
}
}
}
}
プログラム解説の続き
さて、いよいよこのプログラムの核心部分
for(i=0;i<j;i++){
c[i]+=(d[i]+e[i]);
c[i+1]+=c[i]/n;
c[i]=c[i]%n;
}
に到達しました。
04023
00214
すなわち、
d[0]=3,d[1]=2,d[2]=0,d[3]=4,d[4]=0
e[0]=4,e[1]=1,e[2]=2,e[3]=0,e[4]=0
を例に取り、動きを追ってみましょう。以下の説明では、
c[i]+=(d[i]+e[i]);
c[i+1]+=c[i]/n;
は
c[i]=c[i]+(d[i]+e[i]);
c[i+1]=c[i+1]+c[i]/n;
であることに注意して、読んでください。
i=0のとき、
c[0]=c[0]+(d[0]+e[0]);
c[i]はすべて0に初期化されていますから、
c[0]=0+(d[0]+e[0]);
すなわち、
c[0]=0+(3+4);
の計算から、
c[0]=7;
現時点では、
となっていて、正確ではありませんが、心配要りません。
最後まで進行を見て頂けれければと思います。
次の行
c[i+1]=c[i+1]+c[i]/n;
は
c[0+1]=c[0+1]+c[0]/5;
c[1]=c[1]+c[0]/5;
c[1]=0+7/5;
の計算から、
c[1]=1;
これで、5の位(左から2番目の位)繰り上がり
に成功してます。
さて、最後の行の
c[i]=c[i]%n;
によって、
c[0]=c[0]%5;
から
c[0]=7%5;
c[0]=2;
と第1の位が修正され、
1の位(左から1番目の位)も正しくなり、
5の位(左から2番目の位)への繰り上がりも含めて
1番左の計算は成功裏に終わったことになります。