1. 程式人生 > >高精度除以高精度

高精度除以高精度

原 題:

    高精除以高精,求它們的商和餘數。

演算法分析:

   高精除以高精是用減法模擬除法,對被除數的每一位都減去除數,一直減到當前位置的數字(包括前面的餘數)小於除數(由於每一位的數字小於10,所以對於每一位最多進行10次計算)

代 碼:

#include<iostream>

#include<cstdio>
#include<cstring>
using namespace std;
int a[101],b[101],c[101],d,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;
}