精度问题是C语言、Python语言循环结构部分经常出现的问题:
7-2 求给定精度的简单交错序列部分和 (15 分)
本题要求编写程序,计算序列部分和 1 - 1/4 + 1/7 - 1/10 + ... 直到最后一项的绝对值不大于给定精度eps。
输入格式:
输入在一行中给出一个正实数eps。
输出格式:
在一行中按照“sum = S”的格式输出部分和的值S,精确到小数点后六位。题目保证计算结果不超过双精度范围。
输入样例1:
4E-2
输出样例1:
sum = 0.854457
输入样例2:
0.02
输出样例2:
sum = 0.826310
代码:
import math
eps = float(input())
n = 0
sum = 1.0
flag = 1
while(abs(1.0/(3*n+1))>eps):
n += 1
for i in range(1,n+1):
flag *= -1
sum += 1.0/(3*i+1)*flag
print("sum = {:.6f}".format(sum))
7-9 求误差小于输入值的e的近似值 (20 分)
代码实现(用python math库中自带的factorial函数即求阶乘函数)
import math
eps = float(input())
n = 0
while(sum(1.0/math.factorial(i) for i in range(n+2))-sum(1.0/math.factorial(i) for i in range(n+1))>=eps):
n += 1
print("{:.6f}".format(sum(1.0/math.factorial(i) for i in range(n+2))))