求大数阶乘，思路应该没错，但是求15以上的阶乘时会卡住

DarrenTitor

1.   .MODEL SMALL,STDCALL
   
2.   .386
   
3.    
   
4. 
   
5.   .DATA
   
6.   digit db 0
   
7.   n     db ?
   
8.   i     db 0
   
9.   j     db 0
   
10.   num   db 0
    
11.   array db 10000 dup(0)
    
12. 
    
13. 
    
14.    
    
15.   .CODE
    
16.   main proc
    
17.     ;装载数据段地址
    
18.     mov ax,@data
    
19.     mov ds,ax
    
20. 
    
21. ;-----------------把输入存到n中
    
22.     mov  bx, 0
    
23.     Newchar:
    
24.       mov  ah, 1
    
25.       int  21h
    
26.       sub  al, 30h
    
27.       jl  endinput
    
28.       cmp  al, 9   
    
29.       jg  endinput
    
30.       cbw  
    
31.       xchg   ax, bx
    
32.       mul   cx
    
33.       xchg  ax, bx   
    
34.       add   bx, ax
    
35.       jmp   newchar
    
36.     endinput:
    
37.     mov n,bl
    
38. ;-----------------
    
39. 
    
40. 
    
41.     mov al,1
    
42.     mov array[0],al
    
43. 
    
44.     mov al,1
    
45.     mov digit,al
    
46. 
    
47.     mov ch,0
    
48.     mov cl,n
    
49.     sub cl,2
    
50.     firstloop:
    
51.       ;i = n - cx = al
    
52.       mov ah,0
    
53.       mov al,n
    
54.       sub ax,cx
    
55.       mov i,al
    
56.       ;num = 0 = dl,dh也用了
    
57.       mov dh,0
    
58.       mov dl,0
    
59.       mov num,dl
    
60. 
    
61.       ;------------ j = 0 = bx
    
62.       mov bx,0
    
63.       secondloop:
    
64.         ; --------temp=a[j]*i+num;
    
65.         mov al,i
    
66.         mov ah,array[bx]
    
67.         mul ah;->ax
    
68.         add ax,dx
    
69.         mov dh,10;暂时借用dh，一会儿要再改回dh=0
    
70.         div dh
    
71.         ; ----------a[j]=temp%10;
    
72.         mov array[bx],ah
    
73.         ; -----------num=temp/10;
    
74.         mov dl,al
    
75.         mov num,dl
    
76.         mov dh,0
    
77. 
    
78.         add bx,1
    
79.         mov j,bl
    
80.         cmp bl,digit
    
81.         jl secondloop
    
82. 
    
83.       whileloop:
    
84.         mov dl,num
    
85.         mov dh,0
    
86.         cmp dx,0
    
87.         je tonext
    
88.         mov ax,dx
    
89.         mov dh,10;暂时借用dh，一会儿要再改回dh=0
    
90.         div dh
    
91.         mov dh,0
    
92.         mov j,bl
    
93.         mov bh,0;原本bx是j，暂时改成digit
    
94.         mov bl,digit
    
95.         mov array[bx],ah
    
96.         ; num=num/10;
    
97.         mov dl,al
    
98.         ; digit++;
    
99.         add bl,1
    
100.         mov digit,bl
     
101.         mov bl,j;还原bx的值j
     
102. 
     
103.         cmp dx,0
     
104.         jne whileloop
     
105. 
     
106.       tonext:
     
107.         mov al,i
     
108.         add al,1
     
109.         mov i,al
     
110.         sub cx,1
     
111.         cmp al,n
     
112.         jle firstloop
     
113. 
     
114. 
     
115. 
     
116. 
     
117.       ; ------------逆序打印一下
     
118.       mov bh,0
     
119.       mov bl,digit
     
120.       sub bl,1
     
121.       
     
122.     printloop:
     
123.       MOV  DL,array[bx]
     
124.       add dl,30h
     
125.       mov ah,2
     
126.       int 21h
     
127.       sub bx,1
     
128.       cmp bx,0
     
129.       jnl printloop
     
130. 
     
131.       jmp exit;
     
132. 
     
133.     exit:
     
134.       mov ah,4ch
     
135.       int 21h
     
136.     main endp
     
137.     end main

复制代码

原文链接：https://blog.csdn.net/lisp1995/article/details/52403507

编译用的是MASM，我基本是完全照搬上边这篇文章的思路，只是c语言改写成汇编
(应该有些地方不用把寄存器数据写回到数据段，我这里为了跑通写的繁琐了，之后还可以简化)


问题：
求14以下结果输出正确，但是算15阶乘就会卡住
推测一，有些变量用db存储，会溢出，但是digit在算15!时应该只比10多一点，其他变量就更小了
推测二，来自dosbox的未知限制？这个不确定，我再debug试试

既然能算出14!，按照链接中这个思路应该更大的也可以算，虽然代码有点长，各位能提出建议我会很感谢的

wp231957

超范围了啊
更换算法

DarrenTitor

wp231957 发表于 2020-4-30 15:33
超范围了啊
更换算法

这个算法就是为了解决超范围的问题吧（？

wp231957

DarrenTitor 发表于 2020-4-30 15:38
这个算法就是为了解决超范围的问题吧（？

你这代码16位的？，不都是eax   ebx   ecx   edx
64位的，没玩过，不知道是啥

SВ

超过位数了，就使用多个寄存器或者多个内存单元分别存放数的高位和低位，目前没有必要做这种题目，理解循环指令的工作流程即可。