Python绘制Aliquot序列对数图时遇TypeError: int无log属性问题求助

可能的错误原因 #

• 序列中出现0或负数：Aliquot序列迭代到后期可能触达0（比如当数值为1时，真因数和为0，后续迭代都会维持0），而对数函数仅对正数有效，numpy的log ufunc无法处理0值的Python int类型，直接触发报错。
• 数据类型不统一：即便log_base转成了float，若序列里混合了Python int和numpy数值类型，也会导致ufunc处理时类型不兼容。
• 长序列触发终止条件未处理：101次迭代可能已经让序列进入终止循环或出现无效值，代码未提前过滤就直接执行对数转换。

具体解决方案 #

1. 过滤无效值再处理 #

绘图前先筛掉序列里的非正数值，只保留能计算对数的正数：

# 用列表推导式过滤
valid_seq = [num for num in sequence if num > 0]

# 或者用numpy批量处理
import numpy as np
seq_arr = np.array(sequence, dtype=np.float64)
valid_seq_arr = seq_arr[seq_arr > 0]

2. 统一序列数据类型 #

把整个序列转换成numpy浮点类型，避免Python原生int和numpy类型混合的问题：

seq_arr = np.asarray(sequence, dtype=np.float64)
# 再计算对数变换
log_seq = np.log(seq_arr) / np.log(log_base)

3. 加异常值检查逻辑 #

在计算对数前先校验序列，提前提示或处理无效值：

if any(num <= 0 for num in sequence):
    print("警告：序列包含非正数值，将跳过这些值绘图")
    sequence = [num for num in sequence if num > 0]

# 后续执行对数转换和绘图操作

4. 验证长序列的具体数值 #

打印长序列末尾的数值，确认是否出现0或异常值：

print("序列最后10个值：", sequence[-10:])

如果发现0，说明序列已终止，绘图时可以截断到最后一个有效正数的位置。

代码修改示例（graph.py片段） #

假设原对数转换代码如下：

log_sequence = np.log(sequence) / np.log(log_base)

修改为：

import numpy as np

# 转换为numpy浮点数组并过滤无效值
seq_arr = np.array(sequence, dtype=np.float64)
valid_mask = seq_arr > 0
valid_seq = seq_arr[valid_mask]

if len(valid_seq) == 0:
    raise ValueError("序列中没有可用于对数变换的正数值")

# 计算对数变换
log_seq = np.log(valid_seq) / np.log(log_base)
# 绘图时x轴索引对应过滤后的数据
x_axis = np.arange(len(sequence))[valid_mask]
plt.plot(x_axis, log_seq)

内容的提问来源于stack exchange，提问作者Dimitar Todorov