高精度模板

　　我其实只是存个模板……
     真的只是存个模板……
     而已啦……
     支持高精度加减乘除以及快速幂……但是没有负数……（少用而且不想写……要的话自己写写符号位……）

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int bi=1e4,MN=300;
char c[1001];
struct big{
    int a[MN];
    inline big(){
        memset(a,0,sizeof(a));
        a[0]=1;
    }
    inline void read(){
        register int i,j,k;
        memset(a,0,sizeof(a));
        scanf("%s",c);k=strlen(c);
        a[0]=(k+3)/4;
        for (i=0;i<k;i++) j=(k-i+3)/4,a[j]=a[j]*10+c[i]-48;
    }
    inline void pr(){
        register int i;
        printf("%d",a[a[0]]);
        for (i=a[0]-1;i;i--) printf("%04d",a[i]);
    }
    inline big operator =(int x){
        if (x==0){
            memset(a,0,sizeof(a));
            a[0]=1;
        }
        a[0]=0;
        while (x){
            a[0]++;
            a[a[0]]=x%bi;
            x/=bi;
        }
        if (!a[0]) a[0]=1;
    }
    inline big(int x){
        *this=x;
    }
    inline void gl(){
        while(!a[a[0]]&&a[0]>1) a[0]--;
    }
    inline void operator =(big x){
        register int i;
        a[0]=x.a[0];
        for (i=1;i<=a[0];i++) a[i]=x.a[i];
    }
    inline bool operator >(big y){
        if (a[0]!=y.a[0]) return a[0]>y.a[0];
        for (register int i=a[0];i;i--){
            if (a[i]!=y.a[i]) return a[i]>y.a[i];
        }
        return 0;
    }
    inline bool operator >=(const big y){
        if (a[0]!=y.a[0]) return a[0]>y.a[0];
        for (register int i=a[0];i;i--){
            if (a[i]!=y.a[i]) return a[i]>y.a[i];
        }
        return 1;
    }
    inline bool operator <(big y){
        if (a[0]!=y.a[0]) return a[0]<y.a[0];
        for (register int i=a[0];i;i--){
            if (a[i]!=y.a[i]) return a[i]<y.a[i];
        }
        return 0;
    }
    inline bool operator <=(big y){
        if (a[0]!=y.a[0]) return a[0]<y.a[0];
        for (register int i=a[0];i;i--){
            if (a[i]!=y.a[i]) return a[i]<y.a[i];
        }
        return 1;
    }
    inline bool operator ==(big y){
        if (a[0]!=y.a[0]) return 0;
        for (register int i=a[0];i;i--){
            if (a[i]!=y.a[i]) return 0;
        }
        return 1;
    }
    inline bool operator !=(big y){
        return !(*this==y);
    }
    inline bool operator ==(int y){
        big x=y;
        return *this==x;
    }
    inline bool operator !=(int y){
        return !(*this==y);
    }
    inline void swap(big &a,big &b){
        big x=a;a=b;b=x;
    }
    inline big operator +(big x){
        big r;
        if (a[0]<x.a[0]) r.a[0]=x.a[0];else r.a[0]=a[0];
        for (register int i=1;i<=r.a[0];i++) r.a[i]=(i<=a[0]?a[i]:0)+(i<=x.a[0]?x.a[i]:0);
        for (register int i=1;i<=r.a[0];i++)
        if (r.a[i]>=bi){
            r.a[i]-=bi;r.a[i+1]++;
            if (i==r.a[0]) r.a[0]++;
        }
        return r;
    }
    inline big operator +(int x){
        big r;r=*this;
        r.a[1]+=x;
        for (register int i=1;i<=r.a[0];i++)
        if (r.a[i]>=bi){
            r.a[i]-=bi;r.a[i+1]++;
            if (i==r.a[0]) r.a[0]++;
        }else break;
        return r;
    }
    inline big operator -(big x){
        if (*this<x) swap(*this,x);
        register int i;
        big r;
        if (a[0]<x.a[0]) r.a[0]=x.a[0];else r.a[0]=a[0];
        for (i=1;i<=r.a[0];i++) r.a[i]=a[i]-(i<=x.a[0]?x.a[i]:0);
        for (i=1;i<=r.a[0];i++)
        if (r.a[i]<0){
            r.a[i+1]--;r.a[i]+=bi;
        }
        r.gl();
        return r;
    }
    inline big operator -(int x){
        register int i=1;
        big r;r=*this;
        r.a[1]-=x;
        while (r.a[i]<0) r.a[i+1]--,r.a[i]+=bi,i++;
        r.gl();
        return r;
    }
    inline big operator *(big y){
        register int i,j;
        big r;r.a[0]=a[0]+y.a[0]-1;
        for (i=1;i<=a[0];i++)
        for (j=1;j<=y.a[0];j++){
            r.a[i+j-1]+=a[i]*y.a[j];
            if (r.a[i+j-1]>=bi){
                r.a[i+j]+=r.a[i+j-1]/bi;
                r.a[i+j-1]%=bi;
                if (i+j-1==r.a[0]) r.a[0]++;
            }
        }
        for (i=1;i<=r.a[0];i++)
        if (r.a[i]>=bi){
            r.a[i+1]+=r.a[i]/bi;
            r.a[i]%=bi;
            if (i==r.a[0]) r.a[0]++;
        }
        return r;
    }
    inline big operator *(int y){
        register int i,j;
        big r;r.a[0]=a[0];
        for (i=1;i<=a[0];i++)
        r.a[i]=a[i]*y;
        for (i=1;i<=r.a[0];i++)
        if (r.a[i]>=bi){
            r.a[i+1]+=r.a[i]/bi;
            r.a[i]%=bi;
            if (i==r.a[0]) r.a[0]++;
        }
        return r;
    }
    inline big half(){
        register int i,j;big r=*this;
        for (i=r.a[0];i>1;i--) r.a[i-1]+=(r.a[i]%2)*bi,r.a[i]>>=1;
        r.a[1]>>=1;
        r.gl();
        return r;
    }
    inline big operator /(big y){
        register int i,j;
        big r,l,mid,rq=*this;
        r.a[0]=rq.a[0]+1;r.a[r.a[0]]=1;
        while(r>l){
            mid=(l+r+big(1)).half();
            if (mid*y<=rq) l=mid;else r=mid-big(1);
        }
        return l;
    }
    inline big operator %(big y){
        register int i,j;
        big rq=*this;
        return rq-(rq/y*y);
    }
    inline big operator ^(int y){
        big ans=1;
        big rq=*this;
        while(y){
            if (y&1) ans=ans*rq;
            y>>=1;
            rq=rq*rq;
        }
        return ans;
    }
    inline big operator ^(big y){
        big ans=1;
        big rq=*this;
        while(y!=0){
            if (y.a[1]&1) ans=ans*rq;
            y=y/2;
            rq=rq*rq;
        }
        return ans;
    }
};
int main(){
}

　　这是一个int128模版，用两个unsigned long long存储

#define ll unsigned long long
ll M[2315];
struct MMH{
    ll a,b;
    MMH(ll _a=0,ll _b=0):a(_a),b(_b){}
    void half(){
        b>>=1;
        if (a&1) b|=1ull<<63;
        a>>=1;
    }
    void two(){
        a<<=1;
        if (b&(1ull<<63)) a|=1;
        b<<=1;
    }
};
MMH operator +(const MMH &x,const MMH &y){
    MMH z;
    z.a=x.a+y.a;z.b=x.b+y.b;if (z.b<x.b||z.b<y.b) z.a++;
    return z;
}
MMH operator +(const MMH &x,const int &y){
    MMH z=x;
    if ((z.b+y)<z.b) z.a++;z.b+=y;
    return z;
}
MMH operator -(const MMH &x,const MMH &y){
    MMH z;
    z.a=x.a-y.a;z.b=x.b-y.b;if (y.b>x.b) z.a--;
    return z;
}
MMH operator -(const MMH &x,const int &y){
    MMH z=x;
    if (z.b<y) z.a--;z.b-=y;
    return z;
}
bool operator <(const MMH &x,const MMH &y){return x.a==y.a?x.b<y.b:x.a<y.a;}
bool operator >(const MMH &x,const MMH &y){return x.a==y.a?x.b>y.b:x.a>y.a;}
bool operator !=(const MMH &x,const MMH &y){return x.a!=y.a||x.b!=y.b;}
bool operator ==(const MMH &x,const MMH &y){return x.a==y.a&&x.b==y.b;}
MMH operator *(const MMH &x,const MMH &y){
    if (x.a){
        TTT++;
        SSS-=clock();
        MMH z=x,mmh=MMH(0,0),u=y;
        while (u.a|u.b){
            if (u.b&1) mmh=mmh+z;
            z.two();u.half();_mmh++;
        }
        SSS+=clock();
        return mmh;
    }
    ll s=(1ull<<32)-1;
    ll a=(x.b>>32)*(y.b>>32),_b1=(x.b&s)*(y.b>>32),b=_b1+(x.b>>32)*(y.b&s),c=(x.b&s)*(y.b&s),d;
    a+=(b<_b1)*1ull<<32;
    d=c+((b&s)<<32);a+=d<c;a+=b>>32;
    return MMH(a,d);
}