51Nod-1057-N的阶乘

51Nod-1057 N的阶乘

1057 N的阶乘

输入N求N的阶乘的准确值。

Input
输入N(1 <= N <= 10000)

Output
输出N的阶乘

Input示例
5
Output示例
120

就python而言，直接暴力计算就行，不必担心变量的数据范围。

代码如下：

while True:
try:
n, ans = int(input()), 1
for i in range(1, n+1):
ans *= i
print(ans)
except EOFError:
break

还有一种优化阶乘计算的解法，具体解法的方法在这里就不详细讲解了，在这里打开大佬的博客传送门——N 的阶乘的另一种解法

代码如下：

maxn = 10001
#判段素数的数组
isprime = [True]*maxn
#存入各个素数的阶数
cnt = [0]*maxn
#素数列表
prime = []

#素数打表
def Prime(n):
for i in range(2, n):
if isprime[i]:
prime.append(i)
for k in prime:
if i*k >= n:
break
isprime[i*k] = False

#计算素数的阶数
def func(n):
for p in prime:
if p > n:
break
sum = 0
num = p
while num <= n:
sum += n // num
num = num*p
cnt[p] = sum

#乘法快速幂
def power(a, b):
ans = 1
while b:
if b&1:
ans = ans*a
a = a*a
b >>= 1
return ans

Prime(maxn)
while True:
try:
n, ans= int(input()), 1
func(n)

for p in prime:
if p > n:
break
ans *= power(p, cnt[p])

print(ans)
except EOFError:
break

此题如果用C语言解题，有些棘手，因为阶乘算出的结果有些连长整型都无法存储。那么就需要分段存入一数组里，每次超过某一个数后，如以1e8为界限，把数据分成两段存入。输出的时候倒序输出数组里存储的结果。

C语言解法的大佬博客—— 51Nod－1057－N的阶乘