您的位置:首页 > 其它

大数除法模板

2013-08-15 12:35 148 查看
/*高精度除低精度求商模板*/
/*大数除法 ------除数为int范围*/
#include<iostream>
#define N 1000
using namespace std;
void division(char * src,int n)
{
int len = strlen(src),i,k,t=0,s=0;
char dest
;
bool flag = true;    //商是否有了第一个有效位,防止商首部一直出现0
for(i=0,k=0; i<len; i++)
{
t = s*10+(src[i]-48);    //新余数
if(t/n>0 || t==0)        //余数为0要修改商
{
dest[k++] = t/n+48,s = t%n,flag = false;
}
else                    //不够除,修改余数
{
s = t;
if(!flag)            //商已经有有效位了,补零
dest[k++] = '0';
}
}
for(i=0;i<k;i++)
cout<<dest[i];
cout<<endl;
}
int main()
{
char num
;
int n;
while(scanf("%s%d",num,&n)!=EOF)
{
division(num,n);
}
return 0;
}


大数除以大树

/*大数除法---高精度除高精度*/
/*
1.a.size<b.size 返回-1
2.a.size=b.size && a-b<0 返回-1
3.a.size=b.size && a-b=0 返回0
*/
#include<iostream>
#include<cstring>
#include<string>
#define N 2000
using namespace std;
//判断a.size 与b.size 的关系 以及做减法
int judge(char a[],int a1,char b[],int b1)
{
int i;
if(a1<b1)return -1;//a.size<b.size
bool flag=false;
if(a1==b1) //a.size==b.size && a<b
{
for(i=a1-1;i>=0;i--)
if(a[i]>b[i])
flag=true;
else if (a[i]<b[i])
{
if(!flag) return -1;
}
}
for(i=0;i<a1;i++)//前提b中b1---a1部分必须为'0'
{
a[i]=a[i]-b[i]+48;//'0'的ASCII为48
if((a[i]-'0')<0)
{
a[i]=a[i]+10;a[i+1]=a[i+1]-1;
}
}
for(i=a1-1;i>=0;i--) //返回被除数的长度
if(a[i]!='0')
return (i+1);
return 0;//a.size==b.size&&a=b的情况
}
string division(string a,string b)
{
char x1
,x2
;
int ans
;
int a_len,b_len,i,j;
a_len=a.length();
b_len=b.length();
/*初始化部分*/
/*********************************************/
memset(x1,'0',sizeof(x1));
memset(x2,'0',sizeof(x2));
memset(ans,0,sizeof(ans));
for(i=a_len-1,j=0;i>=0;i--)
x1[j++]=a[i];
for(i=b_len-1,j=0;i>=0;i--)
x2[j++]=b[i];
/*********************************************/
/*分析部分*/
/*********************************************/
if(a_len<b_len) return "0";
int temp_len=judge(x1,a_len,x2,b_len);
if(temp_len<0)return "0";
if(temp_len==0)return "1";
ans[0]++;//减掉一次,商加1
int ntimes=temp_len-b_len;
if(ntimes<0)
return "1";
else if(ntimes>0)
//扩充数位,加快减法。
{
for(i=temp_len-1;i>=0;i--)
if(i>=ntimes)
x2[i]=x2[i-ntimes];
else
x2[i]='0';
}
b_len=temp_len;
/*********************************************/
/*加快除法的部分*/
/********************************************/
for(j=0;j<=ntimes;j++)
{
int ntemp;
while((ntemp=judge(x1,temp_len,x2+j,b_len-j))>=0)
{
temp_len=ntemp;
ans[ntimes-j]++;
}
}
/*********************************************/
/*处理最后结果进位部分*/
/*********************************************/
for(i=0;i<N;i++)
if(ans[i]>=10)
{
ans[i+1]+=ans[i]/10;
ans[i]%=10;
}
/*********************************************/
/*返回string类型*/
/*********************************************/
int k=N-1;
string c="";
while(ans[k]==0&&k>0)k--;
for(i=k;i>=0;i--)
c+=(ans[i]+'0');
/*********************************************/
return c;
}
int main()
{
string a,b;
while(cin>>a>>b)
{
cout<<division(a,b)<<endl;
}
return 0;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: