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