高精度全家桶 High Accuracy Arithmetics

未验证

#include<bits/stdc++.h>
using namespace std;

struct Ha //没有说负负得正。。。 
{
    vector<int> dat; bool neg;
    inline int sz()
    {
        return dat.size();
    }
    inline int&operator[](int i)
    {
        return dat[i];
    }
    inline void resz(int n)
    {
        dat.resize(n);
    }
    friend Ha operator + (Ha a, Ha b)
    {
        Ha c; int len = max(a.sz(), b.sz());
        a.resz(len), b.resz(len), c.resz(len+1);
        for(int i=0; i<len; i++)
        {
            c[i]+=a[i]+b[i];
            if(c[i]>=10) c[i]-=10, c[i+1]++;
        }
        if(!c[len]) c.resz(len);
        return c;
    }
    friend bool operator == (Ha a, Ha b)
    {
        if(a.sz()^b.sz()) return 0;
        for(int i=a.sz()-1; ~i; --i)
            if(a[i]^b[i]) return 0;
        return 1;
    }
    friend bool operator < (Ha a, Ha b)
    {
        if(a.sz()^b.sz()) return a.sz()<b.sz();
        for(int i=a.sz()-1; ~i; --i)
            if(a[i]^b[i]) return a[i]<b[i];
        return 0;
    }
    friend Ha operator - (Ha a, Ha b)
    {
        Ha c;
        if(a==b)
        {
            c.dat.push_back(0);
            return c;
        }
        if(a<b)
        {
            swap(a,b);
            c.neg = 1;
        }
        int len = a.sz(); 
        b.resz(len); c.resz(len);
        for(int i=0; i<len; i++)
        {
            c[i]+=a[i]-b[i];
            if(c[i]<0) c[i]+=10, c[i+1]--;
        }
        while(!c[c.sz()-1]) c.dat.pop_back(); 
        return c;
    }
    friend Ha operator * (Ha a, Ha b)
    {
        Ha c; int len = a.sz() + b.sz();
        c.resz(len);
        for(int i=0; i<a.sz(); i++)
            for(int j=0; j<b.sz(); j++)
                c[i+j] += a[i]*b[j];
        for(int i=0; i<c.sz(); i++)
        {
            if(c[i]>10)
            {
                int k = c[i]/10;
                c[i] -= k;
                c[i+1] += k;
            }
        }
        if(!c[len-1]) c.dat.pop_back();
        return c;
    }
    friend Ha operator / (Ha a, int b)
    {
        Ha c; int len=a.sz(), t=0;
        for(int i=len-1; ~i; --i)
        {
            t = t*10+a[i];
            if(t>=10)
            {
                int n = t/b;
                c.dat.push_back(n);
                t -= n*b;
            }else if(c.sz()) c.dat.push_back(0);
        }
        reverse(c.dat.begin(), c.dat.end());
        return c;
    }
    void read()
    {
        char c=getchar();
        while(c<'0'||c>'9') c=getchar();
        while(c>='0'&&c<='9') dat.push_back(c&15), c=getchar();
        reverse(dat.begin(), dat.end());
    }
    void print()
    {
        if(neg) printf("-");
        for(int i=sz()-1; ~i; --i) printf("%d", dat[i]);
    }
}a;

int main()
{
    a.read();
    a.print();
}

FSYの

#include<bits/stdc++.h>
#define cs const
using namespace std;
#define poly vector<int> 
cs int N = 4e5 + 5;
cs double PI = acos(-1.0);
struct node{
    double x, y; 
    node(double _x = 0, double _y = 0){ x = _x; y = _y; }
    node operator + (cs node &a){ return node(x + a.x, y + a.y); }
    node operator - (cs node &a){ return node(x - a.x, y - a.y); }
    node operator * (cs node &a){ return node(x * a.x - y * a.y, x * a.y + y * a.x); }  
}f[N], g[N]; 
int bit, up, rev[N];
void init(int len){ 
    bit = 0; up = 1; while(up <= len) up <<= 1, bit++;
    for(int i = 0; i < up; i++) rev[i] = (rev[i>>1]>>1) | ((i&1) << (bit-1)); 
}
void FFT(node *a, int typ){
    for(int i = 0; i < up; i++) if(i < rev[i]) swap(a[i], a[rev[i]]);
    for(int i = 1; i < up; i <<= 1){
        node wn(cos(PI / i), typ * sin(PI / i));
        for(int j = 0; j < up; j += (i<<1)){
            node w(1, 0);
            for(int k = 0; k < i; k++, w = w * wn){
                node x = a[k+j], y = a[k+j+i] * w;
                a[k+j] = x + y; a[k+j+i] = x - y;
            }
        }
    }
}
struct Bignum{
    poly a; int flg; // + -  
    #define pb push_back
    friend Bignum Mul(Bignum A, Bignum B){
        Bignum C;
        int n = A.a.size(), m = B.a.size();
        int len = n + m - 1; init(len);
        for(int i = 0; i < up; i++) f[i] = g[i] = node(0, 0);
        for(int i = 0; i < n; i++) f[i].x = A.a[i];
        for(int i = 0; i < m; i++) g[i].x = B.a[i];
        FFT(f, 1); FFT(g, 1); 
        for(int i = 0; i < up; i++) f[i] = f[i] * g[i];
        FFT(f, -1);
        for(int i = 0; i < up; i++){
            int v = (int)(f[i].x / up + 0.5);
            C.a.pb(v);
        } C.a.resize(len); return C;
    }
    friend Bignum operator + (Bignum A, Bignum B){
        Bignum C; int len = max(A.a.size(), B.a.size());
        A.a.resize(len); B.a.resize(len); C.a.resize(len + 1);
        for(int i = 0; i < len; i++){
            C.a[i] += A.a[i] + B.a[i];
            if(C.a[i] >= 10) C.a[i] -= 10, C.a[i + 1]++;
        } if(!C.a[len]) C.a.resize(len); return C;
    }
    friend bool operator == (Bignum A, Bignum B){
        if(A.a.size() != B.a.size()) return false;
        for(int i = A.a.size() - 1; i >= 0; i--){
            if(A.a[i] != B.a[i]) return false;
        } return true;
    }
    friend bool operator < (Bignum A, Bignum B){
        if(A.a.size() != B.a.size()) return A.a.size() < B.a.size();
        for(int i = A.a.size() - 1; i >= 0; i--){
            if(A.a[i] != B.a[i]) return A.a[i] < B.a[i]; 
        } 
    } 
    friend Bignum operator - (Bignum A, Bignum B){
        Bignum C;
        if(A == B){ C.a.pb(0); return C;}
        if(A < B) swap(A, B), C.flg = -1;
        int len = A.a.size();
        C.a.resize(len); B.a.resize(len);
        for(int i = 0; i < len; i++){
            C.a[i] += A.a[i] - B.a[i];
            if(C.a[i] < 0) C.a[i] += 10, C.a[i + 1]--; 
        } while(!C.a[C.a.size() - 1]) C.a.pop_back(); return C;
    }
    friend Bignum operator * (Bignum A, Bignum B){
        Bignum C; int len = A.a.size() + B.a.size();
        if(len <= 512){
            C.a.resize(len); 
            for(int i = 0; i < A.a.size(); i++){
                for(int j = 0; j < B.a.size(); j++){
                    C.a[i + j] += A.a[i] * B.a[j];
                }
            } 
        } 
        else{ C = Mul(A, B); C.a.resize(len); }
        for(int i = 0; i < C.a.size(); i++){
            if(C.a[i] >= 10){
                int k = C.a[i] / 10;
                C.a[i] -= k * 10;
                C.a[i + 1] += k; 
            } 
        } if(!C.a[len - 1]) C.a.pop_back(); return C;
    }
    friend int operator % (Bignum A, int b){
        int nx = 0;
        for(int i = A.a.size() - 1; i >= 0; i--){
            nx = (nx * 10 + A.a[i]) % b;
        } return nx;
    }
    friend Bignum operator / (Bignum A, int b){
        Bignum B; int len = A.a.size(), nx = 0;
        for(int i = len - 1; i >= 0; i--){
            nx = nx * 10 + A.a[i];
            if(nx >= b){
                int k = nx / b; B.a.pb(k); 
                nx -= k * b;
            } else if(B.a.size()) B.a.pb(0);
        } reverse(B.a.begin(), B.a.end()); return B;
    }
    void read(){
        char ch = 0; 
        while(!isdigit(ch)) ch = getchar();
        while(isdigit(ch)) a.pb(ch - '0'), ch = getchar();
        reverse(a.begin(), a.end()); // 小的位在前面  
    }
    void print(){ if(flg == -1) putchar('-'); for(int i = a.size()-1; i >= 0; i--) cout << a[i]; }
}A, B; int x;
int main(){
    A.read(); B.read(); 
    (A + B).print(); puts("");
    (A - B).print(); puts("");
    (A * B).print(); puts("");
    A.read(); scanf("%d", &x);
    int y = A % x, ans = y; 
    while(y) B.a.push_back(y % 10), y /= 10; A = A - B;
    (A / x).print(); puts(""); cout << ans << '\n'; 
    return 0; 
}
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 203,271评论 5 476
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,275评论 2 380
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 150,151评论 0 336
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,550评论 1 273
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,553评论 5 365
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,559评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,924评论 3 395
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,580评论 0 257
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,826评论 1 297
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,578评论 2 320
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,661评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,363评论 4 318
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,940评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,926评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,156评论 1 259
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,872评论 2 349
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,391评论 2 342

推荐阅读更多精彩内容