如何获取有限域上椭圆曲线的全部解点？兼释E.Np(5)计算疑问

一、为什么E.Np(5)=2是合理的 #

椭圆曲线在有限域上的点群计数默认包含无穷远点（作为群的单位元），这是你产生疑惑的核心：

1. 先把曲线方程模5简化：原曲线(E: y^2 = x^3 + 2x + 5)，在(\mathbb{F}_5)中5≡0，方程简化为(y^2 = x^3 + 2x)。
2. 遍历(\mathbb{F}_5)中所有x值验证：
   • x=0：(y^2=0) → y=0，对应仿射点(0,0)
   • x=1：(y^2=1+2=3)，(\mathbb{F}_5)的平方数只有0、1、4，3不是平方数，无解
   • x=2：(y^2=8+4=12≡2)，2不是平方数，无解
   • x=3：(y^2=27+6=33≡3)，无解
   • x=4：(y^2=64+8=72≡2)，无解
3. 仿射点仅(0,0)，加上无穷远点后总点数为2，和E.Np(5)的结果完全一致。

二、简便提取仿射点的方法 #

在Sage中，有几种快速过滤无穷远点的实用方式：

方法1：用is_infinite()筛选 #

直接通过点的方法判断并过滤：

E = EllipticCurve(Zmod(5), (2,5))
# 提取所有非无穷远点
affine_points = [p for p in E.points() if not p.is_infinite()]
# 输出直观的仿射坐标
for p in affine_points:
    print(p.xy())  # 结果：(0, 0)

方法2：手动枚举有限域元素 #

直接遍历x并计算合法的y，过程更透明：

F = GF(5)
E = EllipticCurve(F, (2,5))
affine_points = []
for x in F:
    rhs = x^3 + 2*x + 5
    if rhs.is_square():
        # 获取所有平方根（可能是0或2个）
        for y in rhs.sqrt(all=True):
            affine_points.append(E(x, y))
print(affine_points)  # 结果：[(0 : 0 : 1)]

方法3：直接调用affine_points()方法 #

部分Sage版本支持直接获取非无穷远点列表：

E = EllipticCurve(Zmod(5), (2,5))
print(E.affine_points())  # 结果：[(0, 0)]

内容的提问来源于stack exchange，提问作者Antoni Parellada