0

メストレの和を次のように定義しました。

S(E, N) = {
    my (s = 0.0);
    forprime(p = 2, N,
        my (a = ellap(E, p));
        s += (2-a) / (p+1-a)
    );

    return (s);
}

多項式を定義A(t)しB(t)、曲線の上記の合計を計算したかったのです。y^2 = x^3 + A(t)x^2 + B(t)xいくつかの場合-700 < t < 700は、それを印刷します。だから私は書いた：

for(t = -700, 700, {
    E = ellinit([0, A(t), 0, B(t), 0]);
    if(E == [], , print(t, ": ", S(E, 50000)))
})

できます。A(t, t')ここで、多項式とB(t, t')2 つの変数に対して同じことを行いたいと思います。これ-n < tは、t' < n任意の整数に対して、 (および独立している)nで曲線のみを出力します。これの正しいコードを教えてもらえますか? ありがとう。S(E, 50000) > 5tt'

4

1 に答える 1

2

私が問題を正しく理解していれば、次の PARI/GP コードが機能するはずです。

n = 100
forvec(P = [[-n, 0], [0, n]], {
    E = ellinit([0, A(P[1], P[2]), 0, B(P[1], P[2]), 0]);
    if (E == [], next());

    my (v = S(E, 50000));
    if (v > 5, print(P, ": ", v))
})

variable は、変数およびPの間隔を表すことに注意してください。0 で上限を設定し、0 で下限を設定しました。tt'tt'

于 2017-02-21T20:13:07.303 に答える