在不同CAD版本中相同算式为何结果不同？

marting

ynhh 发表于 2017-3-29 15:40

你说的这个就是同余定理，网上搜索下有很多资料

1.同余定理：(a+b)%c=((a%c)+(b%c))%c

m%n举例：

123 % n =(((1%n*10%n+2%n)%n*10%n)%n+3%n)%n

利用这个性质我们可以将大整数按位分解存入数组，然后循环求出最终余数。

ynhh

marting 发表于 2017-3-29 16:21

谢谢你的指点因无法用计算器等验证
因与AUTOCAD计算的结果都不同

请你看看这个对不对？
http://bbs.mjtd.com/thread-174967-1-1.html

• (defun lrem (snum mod / r)
  
•   (while (/= "" snum)
  
•     (setq r    (rem (atoi (substr snum 1 9)) mod)
  
•     snum (substr snum 10 (- (strlen snum) 9))
  
•     snum (strcat (itoa r) snum)
  
•     )
  
•     (if (< (atoi snum) mod) (setq snum ""))
  
•   )
  
•   r
  
• )
  
  

newer

ynhh 发表于 2017-3-29 17:09

下面是根据 同余定理写的 超大整数求余数

1. (defun XD::Math:Rem (a b / str i m)
   
2.   (defun _isint (str)
   
3.     (setq c (mapcar '(lambda (x) (- x 48)) (vl-string->list str)))
   
4.     (if (vl-some '(lambda (x) (> x 9)) c)
   
5.       nil
   
6.       c
   
7.     )
   
8.   )
   
9.   (cond
   
10.     ((= (type a) 'INT)
    
11.      (setq str (itoa a))
    
12.     )
    
13.     ((= (type a) 'REAL)
    
14.      (setq str (rtos a))
    
15.     )
    
16.     ((= (type a) 'STR)
    
17.      (setq str a)
    
18.     )
    
19.   )
    
20.   (if (and str (setq str (_isint str)))
    
21.     (progn
    
22.       (setq
    
23.         i 0
    
24.         m 0
    
25.       )
    
26.       (repeat (length str)
    
27.         (setq m (rem (+ (* 10 (rem m b)) (rem (nth i str) b)) b))
    
28.         (setq i (1+ i))
    
29.       )
    
30.     )
    
31.   )
    
32.   m
    
33. )
    
34. 
    

(defun XD::Math:Rem (a b / str i m)<br />   (defun _isint (str)<br />     (setq c (mapcar '(lambda (x) (- x 48)) (vl-string->list str)))<br />     (if (vl-some '(lambda (x) (> x 9)) c)<br />       nil<br />       c<br />     )<br />   )<br />   (cond<br />     ((= (type a) 'INT)<br />      (setq str (itoa a))<br />     )<br />     ((= (type a) 'REAL)<br />      (setq str (rtos a))<br />     )<br />     ((= (type a) 'STR)<br />      (setq str a)<br />     )<br />   )<br />   (if (and str (setq str (_isint str)))<br />     (progn<br />       (setq<br />         i 0<br />         m 0<br />       )<br />       (repeat (length str)<br />         (setq m (rem (+ (* 10 (rem m b)) (rem (nth i str) b)) b))<br />         (setq i (1+ i))<br />       )<br />     )<br />   )<br />   m<br /> )<br /> <br />

命令: (xd::math:rem "20107019021211566" 56789)
51449


这个数应该是对的。

ynhh

newer 发表于 2017-3-29 21:01

感谢版主

(defun lrem (snum mod / r)
    (while (/= "" snum)
      (setq r    (rem (atoi (substr snum 1 9)) mod)
      snum (substr snum 10 (- (strlen snum) 9))
      snum (strcat (itoa r) snum)
      )
      (if (< (atoi snum) mod) (setq snum ""))
    )
    r
  )

这个据说可计算到32位
不知版主这个有没有可计算到多少位的说法
经对比结果是一样的
但这个感觉是不是更精简一些啊

newer

ynhh 发表于 2017-3-29 21:43

精简代码就几行，给你多的，是判断了你输入的字符串是否是数字组成的，和输入字符串，整数，实数都能给你算。
你贴的那个代码没这些功能。

输入的参数不符合就出错了，比如你输入了 "123a32" 是要出错的。

精简的代码，就几行。下面的代码就没纠错功能了，而且只能算字符串整数

1. (defun lrem (a b / i m c)
   
2.   (setq c (mapcar '(lambda (x) (- x 48)) (vl-string->list a)))
   
3.   (setq i 0 m 0)
   
4.   (repeat (length c)
   
5.     (setq m (rem (+ (* 10 (rem m b)) (rem (nth i c) b)) b))
   
6.     (setq i (1+ i))
   
7.   )
   
8.   m
   
9. )
   
10. 
    

(defun lrem (a b / i m c)<br />   (setq c (mapcar '(lambda (x) (- x 48)) (vl-string->list a)))<br />   (setq i 0 m 0)<br />   (repeat (length c)<br />     (setq m (rem (+ (* 10 (rem m b)) (rem (nth i c) b)) b))<br />     (setq i (1+ i))<br />   )<br />   m<br /> )<br /> <br />

多少位都行，没限制。下面是42位的。

命令: (lrem "123456789123456789012345678912345678901235" 56789)
15992

ynhh

newer 发表于 2017-3-29 21:49

感谢版主的细致讲解
你高精的技术，宽广的情怀
是晓东战无不胜的力量
衷心感谢晓东