第25講 n進数の演算---その1 加法
第9話 プログラム解説その3
プログラム主要再掲
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;
}
のトレースの続けていきます。
i=1のとき、
c[1]=c[1]+(d[1]+e[1]);
はすでに1繰り上がっていますから、
c[1]=1+(2+1);
の計算から
c[1]=4;
次の行
c[1+1]=c[1+1]+c[1]/5;
により、
c[2]=c[2]+c[1]/5;
c[2]=0+4/5;
c[2]=0;
今回は25の位(左から3番目の位)への繰り上がりはありません。
最後の行
c[1]=c[1]%5;
は
c[1]=4%5;
ですから、
c[1]=4;
で修正されません。