高精度除法

/********************************************算法思想：反复做减法，看看从被除数里最多能减去多少个除数，商就是多少;所以算法核心是写一个大整数的减法函数;反复调用该函数进行减法操作;算法步骤：用数组a表示被除数，数组b表示除数，数组res表示商;先用被除数a减去除数b得到差的位数k,同时商+1;再用被除数a减去若干个除数b*(10^k);不够减了，再减去若干个除数b*(10^(k-1))...b*(10^(k-2))...;一直减到不够减为止;每成功减一次，则商的相应位+1;最后再循环处理商的进位问题;********************************************/#include<iostream>#include<cstring>#include<cstdlib>#include<cstdio>#include<climits>#include<algorithm>using namespace std;const int N=1000;int a[N],b[N],res[N*2];char x[N],y[N];int len1,len2;void init(){    len1=strlen(x);    len2=strlen(y);    memset(a,0,sizeof(a));    memset(b,0,sizeof(b));    memset(res,0,sizeof(res));    for(int i=len1-1,j=0; i>=0; i--)    {        a[j++]=x[i]-'0';    }    for(int i=len2-1,j=0; i>=0; i--)    {        b[j++]=y[i]-'0';    }}int substraction(int *p1,int *p2,int len1,int len2)//计算长度为len1的大整数减去长度为len2的大整数的结果的长度{                                                  //减的结果放在数组p1中，不够返回-1,正好返回0    if(len1<len2)        return -1;    bool flag=0;    if(len1==len2)    {        for(int i=len1-1; i>=0; i--)        {            if(p1[i]>p2[i])                flag=1;            else if(p1[i]<p2[i])            {                if(!flag)                    return -1;            }        }    }    for(int i=0; i<len1; i++)//减法    {        p1[i]-=p2[i];        if(p1[i]<0)        {            p1[i]+=10;            p1[i+1]--;        }    }    for(int i=len1-1; i>=0; i--)    {        if(p1[i])            return i+1;    }    return 0;}void output(){    for(int i=0; i<N; i++) //进位    {        if(res[i]>=10)        {            res[i+1]+=res[i]/10;            res[i]%=10;        }    }    bool flag=0;    for(int i=N-1; i>=0; i--)//输出    {        if(flag)        {            printf("%d",res[i]);        }        else if(res[i])        {            printf("%d",res[i]);            flag=1;        }    }    if(!flag)        printf("0");    printf("\n");}void solve(){    init();    if(len1<len2)    {        printf("0\n");        return;    }    len1=substraction(a,b,len1,len2);    if(len1<0)    {        printf("0\n");        return;    }    else if(len1==0)    {        printf("1\n");        return;    }    res[0]++;//减掉一次了，商+1    int k=len1-len2;    if(k<0)//减一次后不能再减了    {        output();        return;    }    else if(k>0)//将数组b乘以10的某次幂，使得其长度与数组a相同    {        for(int i=len1-1; i>=0; i--)        {            if(i>=k)                b[i]=b[i-k];            else                b[i]=0;        }    }    len2=len1;    for(int j=0; j<=k; j++)//先减去若干个b*(10^k)，不够减了再减去若干个b*(10^(k-1))...    {        int temp;        while((temp=substraction(a,b+j,len1,len2-j))>=0)//一直减到不够减为止        {            len1=temp;            res[k-j]++;//每成功减一次，则商的相应位+1        }    }    output();}int main(){    freopen("C:\\Users\\Administrator\\Desktop\\kd.txt","r",stdin);    int t;    scanf("%d",&t);    while(t--)    {        scanf("%s%s",x,y);        solve();    }    return 0;}