高精度模板(转)

转的DQS大神的模板,希望他不要打我。
可以压1~9位,要改一下输出里的printf和BASE还有WIDTH的值

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int size=10010;
const int BASE=(int)1e9;
const int WIDTH=9;
typedef long long LL;
struct bign{
    int len,num[size];
    bign(){memset(num,0,sizeof(num)); len=1;}
    bign(LL x)
    {
        memset(num,0,sizeof(num)); len=0;
        do
        {
            num[++len]=x%BASE;
            x/=BASE;
        }while(x);
    }
};
void scanf(bign &ans)
{
    string in;
    cin>>in;
    int l=in.length();
    ans.len=(l-1)/WIDTH+1;
    for(int i=0;i<ans.len;i++)
    {
        int e=l-i*WIDTH;
        int s=max(0,e-WIDTH);
        sscanf(in.substr(s,e-s).c_str(),"%d",&ans.num[i+1]);
    }
}
void printf(const bign &ans)
{
    printf("%d",ans.num[ans.len]);
    for(int i=ans.len-1;i>=1;i--)
    {
        printf("%09d",ans.num[i]);
    }
}
bool operator <(const bign &a,const bign &b)
{
    if(a.len!=b.len) return a.len<b.len;
    for(int i=a.len;i>=1;i--)
        if(a.num[i]!=b.num[i]) return a.num[i]<b.num[i];
    return false;
}
bign operator +(const bign &a,const bign &b)
{
    bign ans;
    int i=1,x=0;
    while(i<=a.len||i<=b.len)
    {
        x=a.num[i]+b.num[i]+x;
        ans.num[i++]=x%BASE;
        x/=BASE;
    }
    ans.num[i]=x;
    ans.len=i;
    while(ans.len>1&&ans.num[ans.len]==0) ans.len--;
    return ans;
}
bign operator -(const bign &a,const bign &b)
{
    bign ans;
    ans.len=a.len;
    for(int i=1,x=0;i<=a.len;i++)
    {
        x=BASE+a.num[i]-b.num[i]+x;
        ans.num[i]=x%BASE;
        x=x/BASE-1;
    }
    while(ans.len>1&&ans.num[ans.len]==0) ans.len--;
    return ans;
}
bign operator *(const bign &a,const bign &b)
{
    bign ans;
    ans.len=a.len+b.len;
    for(int i=1;i<=a.len;i++)
    {
        LL x=0;
        for(int j=1;j<=b.len;j++)
        {
            x+=ans.num[i+j-1]+(LL)a.num[i]*(LL)b.num[j];
            ans.num[i+j-1]=x%BASE;
            x/=BASE;
        }
        ans.num[i+b.len]=x;
    }
    while(ans.len>1&&ans.num[ans.len]==0) ans.len--;
    return ans;
}
bool smaller(const bign &a,const bign &b,int d)
{
    if(a.len+d!=b.len) return a.len+d<b.len;
    for(int i=a.len;i>=1;i--)
        if(a.num[i]!=b.num[i+d]) return a.num[i]<b.num[i+d];
    return true;
}
void jian(bign &a,const bign &b,int d)
{
    for(int i=1,x=0;i<=a.len-d;i++)
    {
        x=BASE+a.num[i+d]-b.num[i]+x;
        a.num[i+d]=x%BASE;
        x=x/BASE-1;
    }
    while(a.len>1&&a.num[a.len]==0) a.len--;
}
bign operator /(const bign &a,const bign &b)
{
    bign ans;
    ans.len=max(1,a.len-b.len+1);
    bign mid[32],num=a;
    mid[0]=b;
    for(int i=1;i<=30;i++) mid[i]=mid[i-1]*2;
    for(int i=a.len-b.len;i>=0;i--)
    {
        int tmp=1<<30;
        for(int j=30;j>=0;j--)
        {
            if(smaller(mid[j],num,i))
            {
                jian(num,mid[j],i);
                ans.num[i+1]+=tmp;
            }
            tmp>>=1;
        }
    }
    while(ans.len>1&&ans.num[ans.len]==0) ans.len--;
    return ans;
}
bign operator %(const bign &a,const bign &b)
{
    bign ans;
    ans.len=max(1,a.len-b.len+1);
    bign mid[32],num=a;
    mid[0]=b;
    for(int i=1;i<=30;i++) mid[i]=mid[i-1]*2;
    for(int i=a.len-b.len;i>=0;i--)
    {
        int tmp=1<<30;
        for(int j=30;j>=0;j--)
        {
            if(smaller(mid[j],num,i))
            {
                jian(num,mid[j],i);
                ans.num[i+1]+=tmp;
            }
            tmp>>=1;
        }
    }
    return num;
}
int main()
{
    bign a,b;
    scanf(a);
    scanf(b);
    printf(a%b);
    return 0;
}

附上我自己打的,不要喷。。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int MAXN = 10000 + 5;
struct bign
{
	int num[MAXN],len;
};
bign operator + (bign a,bign b)
{
	bign c;
	memset(c.num,0,sizeof(c.num));
	c.len = 0;
	bool flag = false;
	int i;
	for(i = 1;i <= a.len || i <= b.len;i ++)
	{
		c.num[i] = a.num[i] + b.num[i];
		c.num[i] += flag;
		flag = false;
		if(c.num[i] >= 10)
		{
			flag = true;
			c.num[i] -= 10;
		}
	}
	if(flag)
		c.num[i] ++;
	c.len = i - 1;
	return c;
}
char tmp[MAXN];
void scanf_bign(bign &x)
{
	scanf("%s",tmp + 1);
	memset(x.num,0,sizeof(x.num));
	int n = x.len = strlen(tmp + 1);
	for(int i = n;i > 0;i --)
		x.num[i] = tmp[n - i + 1] - '0';
	return;
}
void printf_bign(bign x)
{
	for(int i = x.len;i > 0;i --)
		printf("%d",x.num[i]);
	puts("");
	return;
}
bign a,b;
int main()
{
	scanf_bign(a);
	scanf_bign(b);
	printf_bign(a + b);
	return 0;
}

 

上一篇
下一篇