【模板】高精度

闲来无事发个高精度QwQ

可以压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;
}