第11話 プログラム解説その５
プログラム主要再掲
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=3のとき、
c[3]=c[3]+(d[3]+e[3]);
今回も繰り上がりがありませんでしたから、
c[3]=0+(4+0);
の計算から
c[3]=4;
次の行
c[3+1]=c[3+1]+c[3]/5;
により、
c[4]=c[4]+c[4]/5;
c[4]=0+4/5;
c[4]=0;
今回も625の位(左から５番目の位）への繰り上がりはありません。
最後の行
c[4]=c[4]%5;
は
c[4]=4%5;
ですから、
c[4]=4;
で修正されません。
i=4のとき、
c[4]=c[4]+(d[4]+e[4]);
今回も繰り上がりがありませんでしたから、
c[4]=0+(0+0);
の計算から
c[0]=0;
次の行
c[4+1]=c[4+1]+c[4]/5;
により、
c[5]=c[5]+c[4]/5;
c[5]=0+0/5;
c[5]=0;
今回も左から６番目の位への繰り上がりはありません。
最後の行
c[5]=c[5]%5;
は
c[5]=0%5;
ですから、
c[5]=5;
で修正されません。
以上の進行により
が見事に計算されています。
