不支持 bignum 的语言, 不依赖外部模块, 如何实现 [ + - * ] 运算

substr函数

py 支持  bignum 运算, 酱子

1. 
   
2. A = 123456789987654321
   
3. B = 999999888877777619
   
4. C = 234324411000000001
   
5. D = 342561233000011001
   
6. E = A * B + C - D
   
7. 
   
8. print D # 342561233000011001
   

复制代码

问题: 不支持 bignum 的语言, 不依赖外部模块,
如何实现 [ + - * ] 运算

python 或 perl 或 awk 或 ....代码
谢谢大神

jason680

回复 1# substr函数

想解决问题...

1. bc
$ echo 123456789987654321*999999888877777619+234324411000000001-342561233000011001 | bc
123456776268861449439225475282430699

2. perl   
$ echo 123456789987654321 999999888877777619 234324411000000001 342561233000011001 | perl -Mbigint -lane '{$F[0]+=0;print ($F[0]*$F[1]+$F[2]-$F[3])}'
123456776268861449439225475282430699

3. awk 4.0 with --bignum


4. by yourself

substr函数

回复 2# jason680

J 神
不依赖外部模块, 如何实现[ ]
   

jason680

回复 3# substr函数

1.这是shell版
如果不能用bc,那用sed,awk,...的都不应该用
注：也不用编程了，反正都不能用...

2.说说你真正的问题? 这个问题才有意义...
   

substr函数

回复 4# jason680


   
谢谢 J 神
以语言描述是为了说明实现的过程

酱子, py 支持 bignum

1. 
   
2. A = 12345
   
3. B = 34567890123
   
4. 
   
5. print A + B
   
6. 
   

复制代码

py 支持 bignum, [可以假定 py 不支持],
比如

py 可以如此 实现 +

1. 
   
2. #!/usr/bin/python2
   
3. # coding: utf-8
   
4. 
   
5. A = '12345'
   
6. B = '34567890123678889999999923455555'
   
7. 
   
8. 
   
9. def ADD(a, b):
   
10.     L = len(a) - len(b)
    
11.     if L > 0:
    
12.         b = '0' * L + b
    
13.     else:
    
14.         a = '0' * -L + a
    
15. 
    
16.     N = 0
    
17.     A = [0] * len(a)
    
18.     for i in xrange(len(a) - 1, -1, -1):
    
19.         V = int(a[i]) + int(b[i]) + N
    
20.         N = V / 10
    
21.         A[i] = str(V % 10)
    
22.    
    
23.     S = ''.join(str(E) for E in A)
    
24.     if N > 0: S = str(N) + S
    
25.     return S
    
26. 
    
27. 
    
28. print ADD(B, A)
    
29. 
    
30. 
    

复制代码

减法同理


J 神, 应该如何实现 [ * ] 运算

jason680

回复 5# substr函数

以A =34,B =576为例
length(A)=2, length(B)=3
A=34 < 10^2,  B=576 < 10^3
A * B < 10^2 * 10^3 = 10^5
A * B 数组最大空间(长度)为length(A)+length(B) = 5

自行实现...

xxx34
xx576
=====
xxx24  (4*6)
xx18   (3*6)

xx28   (4*7)
x21    (3*7)

x20    (4*5)
15     (3*5)
====
19584

ly5066113

回复 1# substr函数


sed版dc计算器，满足你的需求：
http://bbs.chinaunix.net/thread-1844319-1-1.html

jason680

回复 5# substr函数

awk for multiply with bignum

$ echo 12345678901234567890 12345678901234567890 | awk '{print $1*$2}'
152415787532388344689912330787336749056

$ echo 12345678901234567890 12345678901234567890 | awk -f bigint.awk
152415787532388367501905199875019052100

$ echo 12345678901234567890*12345678901234567890 | bc
152415787532388367501905199875019052100

$ cat bigint.awk
function abs(n){
    return n>=0 ? n : -n
}
function get_num(num, pos){
    len = length(num)
    if(abs(pos) > len || pos == 0)
        return(0)
    pos = (pos>0) ? pos : len+pos+1
    return substr(num, pos, 1)
}

function mul(n1, n2, len1,len2,len,n,ret){
    len1 = length(n1)
    len2 = length(n2)
    len  = len1 + len2
    delete a
    for(c1=1; c1<=len1; c1++){
        for(c2=1; c2<=len2; c2++)  
            a[c1+c2-1] += get_num(n1,-c1) * get_num(n2,-c2);
    }
    for(n=1; n<=len; n++){
       ret = a[n]%10 ret
       a[n+1] += int(a[n]/10)
    }
    sub(/^0+/, "", ret)
    return ret
}
{
  print mul($1, $2)
}

   

substr函数

回复 7# ly5066113


谢谢 L 神
赞一个!
可惜了. ... sed 的算法不通用呀 [ ]

substr函数

回复 8# jason680


谢谢 J 神
赞一个! [   ]
这个 * 我想了 10 几天,
想不到J 神, 竟如此轻易给实现了

酱子, 翻译了下 py...赞一个!

1. 
   
2. #!/usr/bin/python2
   
3. # coding: utf-8
   
4. 
   
5. def MUL(a, b):
   
6.     A = [int(E) for E in a[::-1]]
   
7.     B = [int(E) for E in b[::-1]]
   
8.     C = [0] * (len(a) + len(b))
   
9.     for i, P in enumerate(A):
   
10.         for j, Q in enumerate(B):
    
11.             C[i + j] += P * Q
    
12.     n = 0
    
13.     for i, v in enumerate(C):
    
14.         C[i] = (n + v) % 10
    
15.         n = (n + v) / 10
    
16.    
    
17.     while C[-1] is 0: C.pop()
    
18.     return ''.join(str(E) for E in C)[::-1]
    
19. 
    
20. 
    
21. A = '99009'
    
22. B = '990199999999999999976233444444444442210506'
    
23. 
    
24. print MUL(A, B)
    
25. print int(A) * int(B)
    
26. 
    
27. # 9911112211111111111111011999877888888888888888888890
    
28. # 9911112211111111111111011999877888888888888888888890
    
29. 
    

复制代码

有了 [ + - * ] 在一个 [ / ] 就圆满了


J 神, 如何实现 [ / ] 运算 [   ]