python求N的阶乘
好心boy 人气:0求N的阶乘
本题要求编写程序，计算N的阶乘。
输入格式：
输入在一行中给出一个正整数 N。
输出格式：
在一行中按照“product = F”的格式输出阶乘的值F，请注意等号的左右各有一个空格。题目保证计算结果不超过双精度范围。
输入样例：
5
输出样例：
product = 120
x = int(input())
a = 1
for i in range(1, x+1):
a = a*i
print("product = %d" % float(a))实现阶乘的三种解法
问题描述：
输入一个正整数n，输出n!的值。
其中n!=123*…*n。
算法描述：
n!可能很大，而计算机能表示的整数范围有限，需要使用高精度计算的方法。使用一个数组A来表示一个大整数a，A[0]表示a的个位，A[1]表示a的十位，依次类推。
将a乘以一个整数k变为将数组A的每一个元素都乘以k，请注意处理相应的进位。
首先将a设为1，然后乘2，乘3，当乘到n时，即得到了n!的值。
输入格式：
输入包含一个正整数n，n<=1000。
输出格式：
输出n!的准确值。
样例输入：
10
样例输出：
3628800
看到这题我首先想到的是两种比较简单的解法，一是循环，二是递归。
解法一：循环
n = int(input()) ns = 1 for i in range(1,n+1): ns = ns*i print(ns)
思路比较简单，就是定义一个变量ns赋予一个初始值1，然后利用for循环直接累乘得到最终结果。
解法二：递归
def factorial(n): if n==1: return n else: return n*factorial(n-1) n = int(input()) res = factorial(n) print(res)
递归也比较好理解,当n == 2,return 2 * 1;n == 3,return 3*(2 * 1);n==4,return 4*(3*(2*1))。以此类推，再将最终的结果赋予res将其打印即可。
这两种方法都比较简单，但很显然都不符合题目要求的 “使用一个数组A来表示一个大整数a，A[0]表示a的个位，A[1]表示a的十位”，所以我们要想办法利用数组来得到n!的结果。
解法三：数组
n= int(input()) ns = [0 for i in range(10000) ] length = 1 ns[0] = length = 1 if n>=2: for i in range(2,n+1): carry = 0 for j in range(length): temp = ns[j] * i + carry carry = int(temp/10) ns[j] = temp % 10 while carry>0: ns[length] += carry%10 length+=1 carry = int(carry/10) while length>0: length -=1 print(ns[length],end='')
接下来我讲下思路：
首先定义一个ns数组用来存储n!的各个位数上的数值，利用for循环给ns加入10000个0值，以方便后面直接根据index对数组进行操作。
然后定义length作为 “数组的长度”（有真实数值的而非自动添加的0） 也即n!的结果的位数。
之后也必须用到for循环进行累乘，但跟解法一的直接累乘不同，这里是乘数(即i)跟各个位上的数分别相乘，若结果大于等于10则carry>0即向前进一位数值为carry，若j循环结束后carry>0则说明需要在当前ns的“长度”上进一位，所以length+1即位数+1，这里carry起的就是判断是否进位的作用，而length则代表着结果的位数。可能这么说有些抽象，下面我们通过分解运行过程来更直观的阐述上面的想法。
例如我们现在需要求5!,分五步,即i循环5次：
①i=1
ns[0] = length =1 , carry = 0 ∴j in range(1)
⑴ j=0
temp = ns[j] * i + carry = ns[0] * i + carry =1*1+0=1 # temp为第j位数与i相乘并加上j-1位数与i相乘后进位的值的结果 carry = int(temp/10) = 1/10 = 0 # carry=0所以不用进位 ns[j] = temp % 10 即 ns[0] = 1 % 10 =1 #只取个位数值作为第j位的值
②i=2
ns[0] = 1, length =1 , carry = 0 ∴j in range(1)
⑴ j=0
temp = ns[j] * i + carry = ns[0] * i + carry =1*2+0=2 # temp为第j位数与i相乘并加上j-1位数与i相乘后进位的值的结果 carry = int(temp/10) = 2 / 10 = 0 # carry=0所以不用进位 ns[j] = temp % 10 即 ns[0] = 2 % 10 =2 #只取个位数值作为第j位的值 #这样就已经的到2!的值了即2
③i=3
ns[0] = 2, length =1 , carry = 0 ∴j in range(1)
⑴ j=0
temp = ns[j] * i + carry = ns[0] * i + carry =2*3+0=6 # temp为第j位数与i相乘并加上j-1位数与i相乘后进位的值的结果 carry = int(temp/10) = 6 / 10 = 0 # carry=0所以不用进位 ns[j] = temp % 10 即 ns[0] = 6 % 10 =6 #只取个位数值作为第j位的值 #这样就已经的到3!的值了即6
④i=4
ns[0] = 6, length =1 , carry = 0 ∴j in range(1)
⑴ j=0
temp = ns[j] * i + carry = ns[0] * i + carry =6*4+0=24 # temp为第j位数与i相乘并加上j-1位数与i相乘后进位的值的结果 carry = int(temp/10) = 24 / 10 = 2 # carry=2>0所以需要向前进2 ns[j] = temp % 10 即 ns[0] = 24 % 10 =4 #只取个位数值作为第j位的值
j循环结束，carry>0执行while循环
while carry>0: ns[length] += carry%10 即 ns[1] += 2 % 10 = 2 #carry = 2 所以向前进2 length+=1 即 length =1+1=2 #位数加一 carry = int(carry/10) = 2 / 10 = 0 # carry = 2<10所以不需要继续进位，while循环结束 ∴length = 2 , ns[0] = 4 ,ns[1] = 2 #这样就得到4!的值ns[1]*10+ns[0] 即 24,输出时可直接倒着打印然后end=''而不需要每位数乘10*n再相加
⑤i=5
ns[0] = 4, ns[1] = 2 length =2 , carry = 0 ∴j in range(2)
⑴ j=0
temp = ns[j] * i + carry = ns[0] * i + carry =4*5+0=20 # temp为第j位数与i相乘并加上j-1位数与i相乘后进位的值的结果 carry = int(temp/10) = 20 / 10 = 2 # carry=2>0所以需要向前进2 ns[j] = temp % 10 即 ns[0] = 20 % 10 =0 #只取个位数值作为第j位的值
⑵ j=1
temp = ns[j] * i + carry = ns[1] * i + carry =2*5+2=12 # temp为第j位数与i相乘并加上j-1位数与i相乘后进位的值的结果 carry = int(temp/10) = 12 / 10 = 1 # carry=1>0所以需要向前进1 ns[j] = temp % 10 即 ns[1] = 12 % 10 =2 #只取个位数值作为第j位的值
j循环结束，carry>0执行while循环
while carry>0: ns[length] += carry%10 即 ns[2] += 1 % 10 = 1 #carry = 1 所以向前进2 length+=1 即 length =2 +1 = 3 #位数加一 carry = int(carry/10) = 1 / 10 = 0 # carry = 1<10所以不需要继续进位，while循环结束 ∴length = 3 , ns[0] = 0 , ns[1] = 2 , ns[2] = 1 # 这样就得到5!的值ns[2] ns[1] ns[0]即 120
这样看下来是否发现和小学的时候学的竖式乘法运算过程很相似，从低位数到高位数（ns[j],j in range(0,length)）依次与乘数（i）相乘，大于十则进位（carry＝temp/10＞0，若ns[length]*i＋carry ＞ 10则length＋1）。
以上为个人经验，希望能给大家一个参考，也希望大家多多支持。
加载全部内容