Counting points on smooth plane quartics
Edgar Costa, David Harvey, Andrew V. Sutherland
TL;DR
This work develops practical, scalable methods to count points on smooth plane quartics by computing Frobenius traces through Cartier–Manin matrices. It introduces three practical average-polynomial-time algorithms that reduce the Cartier–Manin computation to small 3×3 or compressed 16×16/28×28 matrix recurrences, enabling $O(p\log p\log\log p)$ or $O(p^{1/2}\log^2 p)$ time with modest space, and an $O(N\log^3 N)$-time, $O(N)$-space average algorithm across good primes $p\le N$. The core idea relies on a robust algebraic framework that translates the $F^m$ coefficient problem into low-dimensional linear recurrences via the differential equations $\partial_i(FG)=(m+1)(\partial_i F)G$ and carefully designed shift maps, yielding efficient computation of the Cartier–Manin matrix $A_p$ for plane quartics and, hence, $\#X(\mathbb{F}_p)=p+1-a_p$ with $a_p\equiv\operatorname{tr}(A_p)\pmod p$. The practical impact is substantial: these methods substantially extend the feasible range for point counting on genus-3 curves not cyclic covers and provide new pathways to compute $L_p(T)$ modulo $p$ and related arithmetic invariants.
Abstract
We present efficient algorithms for counting points on a smooth plane quartic curve $X$ modulo a prime $p$. We address both the case where $X$ is defined over $\mathbb F_p$ and the case where $X$ is defined over $\mathbb Q$ and $p$ is a prime of good reduction. We consider two approaches for computing $\#X(\mathbb F_p)$, one which runs in $O(p\log p\log\log p)$ time using $O(\log p)$ space and one which runs in $O(p^{1/2}\log^2\!p)$ time using $O(p^{1/2}\log p)$ space. Both approaches yield algorithms that are faster in practice than existing methods. We also present average polynomial-time algorithms for $X/\mathbb Q$ that compute $\#X(\mathbb F_p)$ for good primes $p\le N$ in $O(N\log^3\! N)$ time using $O(N)$ space. These are the first practical implementations of average polynomial-time algorithms for curves that are not cyclic covers of $\mathbb P^1$, which in combination with previous results addresses all curves of genus $g\le 3$. Our algorithms also compute Cartier-Manin/Hasse-Witt matrices that may be of independent interest.
