是时候来一发高精度辣!!

第5次编辑
个人感觉啊,要是一道题目要用高精度,想必就是毒瘤题;
反正我对高精是很虚的;
c++有什么重载运算符,这种东西太高端了,辣鸡不会用啊;
所以我就搞了一个普通的高精;
希望大家指出我的程序哪里可以优化;
小辈感激不尽;


定义

#include<iostream>
#include<cstdio>
#include<cstring>
#define Ll long long
using namespace std;
const Ll NN=1e8;
const int N=8;
struct H{
    Ll a[1250],len;
    H(){memset(a,0,sizeof a);len=1;}
    void clear(){memset(a,0,sizeof a);len=1;}
};

NN就是你的进制,我这里是亿进制;
不建议用更高的进制;
N就是进制有几个0,十进制就是1;
比如一个10000位的数,压好就是1250位,这个注意;
下面是一个构造函数,说明一开始高精是0;


输入输出

void init(H &a){
    string s; cin>>s; 
    int l; 
    a.len=0;
    for(int r=s.length()-1;r>=0;r-=N){
        a.a[++a.len]=0;           
        if(r>=N-1)l=r-N+1;else l=0;
        for(int i=l;i<=r;i++)a.a[a.len]=a.a[a.len]*10+s[i]-48;
    }
}  
void outit(H a){
    printf("%lld",a.a[a.len]);
    for(int i=a.len-1;i;i--){
        for(int k=NN/10;a.a[i]<k;k/=10)printf("0");
        if(a.a[i])printf("%lld",a.a[i]); 
    }printf("\n");  
}
H s;
init(s);
outit(s);

我这里在读入时就反转了数字,输出时再反转;


赋值

H in(Ll x){
    H a;
    if(!x)return a;
    a.len=0;
    while(x){a.a[++a.len]=x%NN;x=x/NN;}
    return a;
}

加法

H jia(H a,H b){
    H c;
    int l=max(a.len,b.len);
    for(int i=1;i<=l;i++){
        c.a[i]+=a.a[i]+b.a[i];
        c.a[i+1]=c.a[i]/NN;
        c.a[i]%=NN;
    }
    if(c.a[l+1])l++;
    c.len=l;
    return c;
}

减法

H jian(H a,H b){
    if(check(a,b)==0){H c;return c;}
    for(int i=1;i<=b.len;i++)a.a[i]-=b.a[i];
    for(int i=1;i<=a.len;i++)
        if(a.a[i]<0){a.a[i]+=NN;a.a[i+1]-=1;}
    if(!a.a[a.len])a.len--;
    return a;
}

乘法

H chen(H a,H b){
    H z;
    z.len=a.len+b.len+2;
    for(int i=1;i<=a.len;i++)
    for(int j=1;j<=b.len;j++)z.a[i+j-1]+=(a.a[i]*b.a[j]);
    for(int i=1;i<=z.len;i++)z.a[i+1]+=z.a[i]/NN,z.a[i]%=NN;    
    while(z.len>1&&!z.a[z.len])z.len--;
    return z;
}

如果一边乘一边%会慢好多;这里就是为什么a数组要开Ll;


快速幂

H ksm(H a,int n){
    if(n==0)return in(1);
    if(n==1)return a;
    H c=ksm(a,n>>1);
    c=chen(c,c);
    if(n&1)c=chen(c,a);
    return c;
}

比较

int check(H a,H b){
    if(a.len>b.len)return 1;
    if(a.len<b.len)return -1;
    for(int i=a.len;i;i--)
        if(a.a[i]!=b.a[i])
            if(a.a[i]>b.a[i])return 1;else return -1;
    return 0;
}
a>b//1;
b<a//-1;
a==b//0;

注意这里要从最高位开始比较,所以要到着来;


判奇偶

bool odd(H a){if(a.a[1]&1)return 1;return 0;}

未完待续

posted @ 2017-02-26 13:12  largecube233  阅读(146)  评论(0编辑  收藏  举报