高精度除法——高精度除以高精度

#include<iostream>
#include<cstring>
using namespace std;
int a[101],b[101],c[101],i;
//輸入函數 
void init(int a[]){
    string s;
    cin>>s;
    a[0]=s.length();
    for(i=1;i<=a[0];i++)
        a[i]=s[a[0]-i]-'0';//減法倒序存儲 
}
//輸出函數
void print(int a[]){
    int i;
    if(a[0]==0){
        cout<<0<<endl;
        return;
    }
    for(i=a[0];i>0;i--)
        cout<<a[i];
    cout<<endl;
    return;  //函數執行完畢回到主程序 
}
//比較函數 
int compare(int a[],int b[]){
    int i;
     if(a[0]>b[0])  
         return 1;
     if(a[0]<b[0])  
         return-1;
     for(i=a[0];i>0;i--){ //如果兩數位數相等，則按位比大小 
         if(a[i]>b[i]) 
             return 1;
          if(a[i]<b[i])
              return -1;  //按位比較若該位數相同，則判斷下一位 
     }
    return 0;//如果返回0則表示兩數相等 
}
//減法模擬除法 
void jian(int a[],int b[]){
    int flag,i;
    flag=compare(a,b);
    if(flag==0){
        a[0]=0;
        return;
    }
    if(flag==1){
        for(i=1;i<=a[0];i++){
            if(a[i]<b[i]){
                a[i+1]--;
                a[i]=a[i]+10;
            }
            a[i]-=b[i];
        }
        while(a[0]>0&&a[a[0]]==0)
            a[0]--;
        return;
    }
}
//復制數組 
void numcpy(int p[],int q[],int det){
    for(int i=1;i<=p[0];i++) 
        q[i+det-1]=p[i];
    q[0]=p[0]+det-1;
    /*
    for(int i=q[0];i>0;i--) 
        cout<<q[i];
    cout<<endl;   
    打印復制后的數字，方便理解算法，此算法主要采用低位補0做減法 
    */
}
//除法計算 
void chugao(int a[],int b[],int c[]){
    int i,tmp[101];
    c[0]=a[0]-b[0]+1;   //商的位數不超過被除數的位數－除數的位數+1 
    for(i=c[0];i>0;i--){  //每次循環確定某位商的的值，從高位開始 
        memset(tmp,0,sizeof(tmp));
        numcpy(b,tmp,i);
        while(compare(a,tmp)>=0){
            c[i]++;
            jian(a,tmp);
        }
    }
    while(c[0]>0&&c[c[0]]==0)  
        c[0]--;
    return;
} 
//主函數
int main(){
    init(a);
    init(b);
    chugao(a,b,c);
    print(c);
    print(a);
    return 0; 
}