REVIEW 3 major objections 4 minor 41 references
Probabilistic Saturations and Alt's Problem
T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read The paper argues that a probabilistic finite-field computation giving 8,652 solutions outside the base locus completes Alt's problem by ruling out any four-bar linkages beyond the 4,326 already found numerically.
desk verdict A well-executed computational study whose advertised proof of certainty for Alt's problem collapses on a false finite-field claim. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the saturation count $g_i(X,\mathbb{C}^n)=\deg V(\Theta\cdot x-1,\Lambda\cdot f)\setminus X$: the number of solutions, counted away from the base locus $X$, of $i$ general affine linear equations in $x$ and $n-i$ general linear combinations of the defining polynomials. To compute it over $\mathbb{Z}_p$, the paper uses the ideal $I_i=\langle\Theta\cdot x-1,\Lambda\cdot f,1-(\mu\cdot f)T\rangle$, whose quotient dimension equals $g_i$ by Theorem 4.1 through the Rabinowitz trick. The probabilistic engine is Corollary 4.14, which states that with entries of $\Theta,\Lambda,\mu$ chosen uniformly from $\mathbb{Z}_p$, the probability that the finite-field quotient dimension equals the true complex value is at least $((p-1)/p)^\nu(1-(g_i+2n(D_{\min}+(r+n)D_{\max})\deg V)/p)$, with $\nu$ a Gröbner-basis size bound. This inequality is what converts a successful finite-field run into a probabilistic certificate and sets the required prime size.
What would settle it
Two observations would settle the matter. For the probability guarantee, test the proof's injectivity claim by evaluating a nonconstant leading-coefficient polynomial from Example 4.6 on the random set $\{0,\ldots,p-1\}$: a polynomial like $x^2-1$ vanishes on about $2/p$ of the set, so a leading coefficient with that behavior would invalidate the bound $((p-1)/p)^\nu$. For the numerical claim, compute $g_0(X,\mathbb{C}^8)$ over $\mathbb{Z}_p$ for a prime larger than $2^{116}$, or obtain a certified rational Gröbner basis for one random trial; any value other than 8652 would refute the claimed completion of Alt's problem.
Extended reading notes
Core claim
The central claim is that Alt's problem has a complete answer: the number of distinct four-bar linkages whose coupler curve interpolates nine general points is 4,326, because the number of isolated solutions to the randomized saturation system, counted outside the base locus $X=V(f_1,\ldots,f_{15})$ of degenerate linkages, is exactly $g_0(X,\mathbb{C}^8)=8652$. The paper reports that, over the finite fields $\mathbb{Z}_p$ for the primes from $2^7+3$ through $2^{27}+29$, the dimension of the quotient $\mathbb{Z}_p[x,T]/I_0(\Theta,\Lambda,\mu)$ is always 8652, where $I_0$ is generated by eight general linear combinations of the fifteen $f_j$ together with the Rabinowitz-trick polynomial $1-(\mu\cdot f)T$. Since a zero-dimensional quotient over $\mathbb{Z}_p$ bounds the degree of the corresponding complex system when the prime is lucky, the paper concludes that no additional solutions can exist and that the upper bound $g_0/2=4326$ is sharp.
Load-bearing premise
The argument's weakest link is the claim, used in Theorem 4.7, that every nonconstant polynomial over $\mathbb{F}_p$ is injective and therefore vanishes on at most $1/p$ of the random choices; if that claim gives way, the lower bound $((p-1)/p)^\nu$ in Corollary 4.14 is not established.
Editorial extensions
If this is right
- If $g_0(X,\mathbb{C}^8)=8652$ is accepted, Alt's problem is complete: exactly 4,326 distinct four-bar linkages and 1,442 distinct coupler curves pass through nine general points, because the known constructions meet the new upper bound.
- The verified values $g_i(X,\mathbb{C}^8)$ for $i=1,\ldots,7$ certify the degrees of the varieties of linkages interpolating $9-i$ general points, upgrading the numerical table of [12] to probabilistic upper bounds for a family of synthesis problems.
- The same saturation framework computes volumes of Newton-Okounkov bodies and characteristic classes such as Euler characteristics, Chern classes, and Segre classes, so the finite-field probability analysis gives explicit confidence levels for those computations.
- Corollary 4.14 supplies an a priori recipe for prime size: for $g_0$, taking $p>2^{116}$ guarantees success probability exceeding 0.99, while the experiments show that much smaller primes already succeed in practice.
- Because each finite-field run computes a zero-dimensional quotient, a successful run over a lucky prime yields an upper bound on the complex degree, not merely a numerical suggestion.
Reading between the lines
- The numerical claim and the probability guarantee should be assessed separately: even if the injectivity assertion behind Theorem 4.7 is not valid, the empirical fact that many primes and floating-point trials return 8652 could still be true and could still make 4,326 the right answer.
- A deterministic completion of Alt's problem would follow from one certified lucky-prime computation: find one parameter set whose Gröbner basis over $\mathbb{Z}_p$ is Pauer lucky, lift the computation to $\mathbb{Q}$, and the upper bound becomes unconditional; the present paper does not carry out that lift.
- The coefficient-randomization trick is not specific to four-bar linkages: any enumerative problem whose base locus is known and whose symmetry group is understood can be bounded by the same $g_0$ computation, with the conics example in the paper serving as a small test case.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops probabilistic saturation techniques for computing the numbers gi(X,C^n), which count solutions to a general affine-linear/linear-combination system outside the base locus X. These numbers arise in Alt's four-bar linkage problem, where g0(X,C^8)/2 is proposed as an upper bound on the number of distinct linkages interpolating nine general points. The manuscript combines certified numerical Hilbert function computations with symbolic upper bounds, then analyzes finite-field Gr"obner basis computations to give probabilistic guarantees. Section 5 reports modular computations matching the known degree g0(X,C^8)=8652, and Section 6 claims that verifying this value theoretically completes Alt's problem.
Significance. If the probabilistic certificate were sound, the paper would provide a practical and rigorous method for upper bounds in Alt's problem and for related projective degree, Segre class, Chern class, and Newton-Okounkov volume computations. The Hilbert-function methodology in Section 3, which certifies lower bounds with alphaCertified and matches them with symbolic upper bounds, is a genuine strength. The examples are detailed, and the reported timings and parameter choices make the experiments reproducible in principle. However, the central finite-field probability argument in Section 4 is invalid as written, and the advertised conclusion in Section 6 substantially overstates what the computations establish.
major comments (3)
- [§4.2, Theorem 4.7, Eq. (16)] The proof of Theorem 4.7 asserts that 'a nonconstant polynomial must map Z_p to itself injectively' and uses this to conclude that each leading coefficient is nonzero with probability at least (p-1)/p. This is false over finite fields: for example, f(t)=t^2-1 is nonconstant and vanishes at two elements of F_p. The error is load-bearing because the factor ((p-1)/p)^nu in (16) is the basis for Theorem 4.9 and Corollary 4.14. Moreover, the paper's own Example 4.6 gives leading coefficients such as c1(W)=(theta_denom1 theta_num2)^3 lambda_denom1 lambda_denom2 lambda_denom3 lambda_num4, which are products of several independently chosen variables, so P(c1 != 0)=((p-1)/p)^6, strictly smaller than (p-1)/p for every p>1. The cleared-denominator leading coefficients in I_i(W) have the same multiplicative structure. Therefore the claimed lower bound in (16) is not valid, and the field-size estimate p>2^116 in Section 5.1 is unsupported.
- [§4.1, Theorem 4.5] The proof of Theorem 4.5 contains a second unsupported claim. From a Zariski dense subset D on which p divides a leading coefficient c_j(W), the proof infers that c_j(W) must be constant because 'a nonconstant polynomial must map dense sets to dense sets' and the multiples of p are not dense in Q. This is false in the classical topology even over Q: t maps to t^2 sends Q to Q_{\ge 0}, which is not dense in Q. Furthermore, the desired conclusion does not follow: a nonzero polynomial all of whose coefficients are divisible by p would satisfy p | c_j(W) for every integer W without being constant. Since Theorem 4.7 refers back to the proof of Theorem 4.5 for the nonconstancy of the leading coefficients, the modular-luckiness framework is not established as written.
- [§6 and Corollary 4.14] The concluding statement that verifying g0(X,C^8)=8652 'theoretically completes Alt's problem by showing that there could be no additional solutions' is not supported by the manuscript. Section 5 reports finitely many successful random trials over selected finite fields and numerical homotopy runs; these are empirical confirmations. Even a corrected version of Corollary 4.14 would at most provide a high-confidence probabilistic statement rather than a theoretical proof. Unless a deterministic certificate is supplied, such as an exact Gr\"obner basis computation over Q or a certified interval/arithmetic argument, the conclusion should be rephrased as probabilistic confirmation, not theoretical completion.
minor comments (4)
- [§4.1, Example 4.6] The word 'Symoblically' should be 'Symbolically'.
- [§5.1, Corollary 5.1] The bound deg(V) <= 7,620,480,000 is stated without showing the B\'ezout calculation; please include the product of the degrees or the elimination argument used, since the reader cannot verify this number from the stated D_min and D_max alone.
- [Definition 4.4] The phrase 'p is larger than any coefficient' should specify that absolute values are meant, since coefficients may be negative.
- [Table 3] It would help to state explicitly which entries are in bold because they agree with Table 2, and to explain the missing entry for g0 at p=2^27+29.
Circularity Check
No circularity: finite-field saturation values are computed independently and compared with benchmarks only after the fact.
full rationale
The claimed derivation is self-contained rather than circular. The quantities g_i(X,C^n) are defined in Eq. (2) by counting solutions outside X to a general linear combination of f, and Theorem 4.1 connects this to the dimension of R[T]/I_i(Theta,Lambda,mu); neither definition presupposes the 8652 count. The Alt upper-bound claim g_0(X,C^8)/2 follows by replacing the specialized coefficients c_j(p_i, bar p_i) with independent parameters, so it is a constructed bound, not a renamed numerical result. In Section 5 the values g_i(X,C^8)^{rand}_{Z_p} are produced by fresh modular Groebner basis computations; Table 3 even shows smaller primes giving values different from 8652 (e.g., 8492 for p=3), so the computation is not hard-wired to the benchmark. The agreement with [12, Table 1] is used only as post hoc confirmation. The authors' self-citations ([12,24,25]) are published background or benchmarks and are not load-bearing inputs forcing the conclusion. The paper's own Section 6 limitation — that Magma resources do not permit p>2^116, so the a priori probability bound cannot be directly realized — and the questionable injectivity assertion in Theorem 4.7 are correctness and confidence concerns, not circularity; they do not make the derivation equivalent to its inputs. No step reduces Eq. (10) to Eq. (2) by construction, and no fitted parameter is renamed as a prediction.
Assumptions & free parameters
assumptions (4)
- standard math Bertini's theorem and Rabinowitz trick: for general Theta, Lambda, mu, the ideal I_i(Theta,Lambda,mu) is unit or zero-dimensional, and g_i(X,C^n) equals dim_Q(R[T]/I_i).
- standard math Pauer lucky prime framework: a prime avoiding all leading coefficients of a Groebner basis over Z yields the same leading monomials over Q and F_p.
- standard math Schwartz-Zippel lemma bounds the probability that a bounded-degree polynomial vanishes on random inputs.
- ad hoc to paper Nonconstant polynomials over F_p are injective, so a leading coefficient is nonzero with probability at least (p-1)/p.
Cite this review
Pith. "Pith review of Probabilistic Saturations and Alt's Problem." pith.science (2026). https://pith.science/paper/G2JT7M7W
@misc{pith2026190806020,
author = {Pith},
title = {Pith review of: Probabilistic Saturations and Alt's Problem},
year = {2026},
howpublished = {\url{https://pith.science/paper/G2JT7M7W}},
note = {Machine review of arXiv:1908.06020}
}
read the original abstract
Alt's problem, formulated in 1923, is to count the number of four-bar linkages whose coupler curve interpolates nine general points in the plane. This problem can be phrased as counting the number of solutions to a system of polynomial equations which was first solved numerically using homotopy continuation by Wampler, Morgan, and Sommese in 1992. Since there is still not a proof that all solutions were obtained, we consider upper bounds for Alt's problem by counting the number of solutions outside of the base locus to a system arising as the general linear combination of polynomials. In particular, we derive effective symbolic and numeric methods for studying such systems using probabilistic saturations that can be employed using both finite fields and floating-point computations. We give bounds on the size of finite field required to achieve a desired level of certainty. These methods can also be applied to many other problems where similar systems arise such as computing the volumes of Newton-Okounkov bodies and computing intersection theoretic invariants including Euler characteristics, Chern classes, and Segre classes.
Figures
Reference graph
Works this paper leans on
- [1]
-
[2]
P. Aluffi. Computing characteristic classes of projective schemes. J. Symb. Comput. , 35(1), 3–19, 2003
work page 2003
-
[3]
H. Alt. ¨Uber die Erzeugung gegebener ebener Kurven mit Hilfe des Gelenkvierecks. ZAMM, 3(1), 13–19, 1923
work page 1923
-
[4]
E.A. Arnold. Modular algorithms for computing Gr¨ obner bases.J. Symb. Comput., 35(4), 403–419, 2003
work page 2003
- [5]
- [6]
-
[7]
D.J. Bates, J.D. Hauenstein, A.J. Sommese, C.W. Wampler. Bertini: Software for Numerical Algebraic Geometry. Available at bertini.nd.edu
-
[8]
D.J. Bates, J.D. Hauenstein, A.J. Sommese, C.W. Wampler. Numerically solving polynomial systems with Bertini, Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 2013
work page 2013
Show all 41 references
-
[9]
Bernstein, A.G Kuˇ snirenko, A.G
D.N. Bernstein, A.G Kuˇ snirenko, A.G. Hovanski˘ ı. Newton polyhedra, Uspehi Mat. Nauk , 31, 201–202, 1976
1976
-
[10]
L. Blum, F. Cucker, M. Shub, S. Smale. Complexity and real computation , Springer-Verlag, New York, 1998
1998
-
[11]
Bosma, J
W. Bosma, J. Cannon, C. Playoust. The Magma algebra system. I. The user language.J. Symb. Comput., 24, 235–265, 1997
1997
-
[12]
Brake, J.D
D.A. Brake, J.D. Hauenstein, A.P. Murray, D.H. Myszka, C.W. Wampler. The complete solution of Alt-Burmester synthesis problems for four-bar linkages. J. Mech. Robotics, 8(4), 041018, 2016
2016
-
[13]
Decker, G.-M
W. Decker, G.-M. Greuel, G. Pfister, M. Sch¨ onemann. Singular 4-1-2 — A computer algebra system for polynomial computations. http://www.singular.uni-kl.de (2019)
2019
-
[14]
Faug´ ere, P
J. Faug´ ere, P. Gianni, D. Lazard, T. Mora. Efficient computation of zero-dimensional Gr¨ obner bases by change of ordering. J. Symb. Comput. , 16(4), 329-344, 1993
1993
-
[15]
W. Fulton. Intersection theory, 2nd ed., Ergebnisse der Mathematik und ihrer Grenzgebiete. 3. Folge. A Series of Modern Surveys in Mathematics, vol. 2, Springer-Verlag, Berlin, 1998. PROBABILISTIC SATURATIONS AND ALT’S PROBLEM 23
1998
-
[16]
Grayson, M
D. Grayson, M. Stillman. Macaulay2, a software system for research in algebraic geometry , available at www.math.uiuc.edu/Macaulay2/
-
[17]
Griffin, J.D
Z.A. Griffin, J.D. Hauenstein, C. Peterson, A.J. Sommese. Numerical computation of the Hilbert function of a zero-scheme. Springer Proceedings in Mathematics & Statistics , 76, 235–250, 2014
2014
-
[18]
G¨ ortz, T
U. G¨ ortz, T. Wedhorn.Algebraic Geometry: Part I: Schemes . Springer Science & Business Media, 2010
2010
-
[19]
Harris, M
C. Harris, M. Helmer. Segre class computation and practical applications. To appear in Math. Comput
-
[20]
J. Harris. Algebraic geometry: a first course . Springer Science & Business Media, 2013
2013
-
[21]
Hauenstein, J.I
J.D. Hauenstein, J.I. Rodriguez. Multiprojective witness sets and a trace test. To appear in Adv. Geom
-
[22]
Hauenstein, A.J
J.D. Hauenstein, A.J. Sommese, C.W. Wampler. Regeneration homotopies for solving systems of poly- nomials. Math. Comp., 80, 345–377, 2011
2011
-
[23]
Hauenstein, F
J.D. Hauenstein, F. Sottile. Algorithm 921: alphaCertified: Certifying solutions to polynomial systems. ACM Trans. Math. Softw. , 38(4), 28, 2012
2012
-
[24]
M. Helmer. Algorithms to compute the topological Euler characteristic, Chern-Schwartz-MacPherson class and Segre class of projective varieties. J. Symb. Comput. , 73, 120–138, 2016
2016
-
[25]
M. Helmer. A direct algorithm to compute the topological Euler characteristic and Chern-Schwartz- MacPherson class of projective complete intersection varieties. Theor. Comput. Sci. , 681, 54–74, 2017
2017
-
[26]
Kaveh, A.G
K. Kaveh, A.G. Khovanskii. Newton-Okounkov bodies, semigroups of integral points, graded algebras and intersection theory. Ann. of Math. (2) , 176(2), 925–978, 2012
2012
-
[27]
Lazarsfeld, M
R. Lazarsfeld, M. Mustat ¸˘ a. Convex bodies associated to linear series.Ann. Sci. ´Ec. Norm. Sup´ er. (4), 42(5), 783–865, 2009
2009
-
[28]
Leykin, J
A. Leykin, J. Yu. Beyond polyhedral homotopies. J. Symb. Comput. , 91, 173–190, 2019
2019
-
[29]
Marco-Buzun´ ariz
M.A. Marco-Buzun´ ariz. A polynomial generalization of the Euler characteristic for algebraic sets.Journal of Singularities, 4, 114–130, 2012
2012
-
[30]
Morgan, A.J
A.P. Morgan, A.J. Sommese. Coefficient-parameter polynomial continuation. Appl. Math. Comput. , 29(2), 123–160, 1989
1989
-
[31]
Morgan, A.J
A.P. Morgan, A.J. Sommese, C.W. Wampler. A product-decomposition bound for Bezout numbers. SIAM J. Num. Anal. , 32(4), 1308–1325, 1995
1995
-
[32]
Okounkov
A. Okounkov. Brunn-Minkowski inequality for multiplicities. Invent. Math., 125(3), 405–411, 1996
1996
-
[33]
F. Pauer. On lucky ideals for Gr¨ obner basis computations. J. Symb. Comput. , 14, 471–482, 1992
1992
-
[34]
Plecnik, R.S
M.M. Plecnik, R.S. Fearing. Finding only finite roots to large kinematic synthesis systems. J. Mech. Robot., 9(2), 021005, 2017
2017
-
[35]
S. Roberts. On three-bar motion in plane space. Proc. London Mathematical Society, III, 286–319, 1875
-
[36]
B. Roth, F. Freudenstein. Synthesis of path-generating mechanisms by numerical methods. ASME Jour- nal of Engineering for Industry, Series B . 85(3):298–306, 1963
1963
-
[37]
Schwartz
J.T. Schwartz. Fast probabilistic algorithms for verification of polynomial identities. J. ACM , 27(4), 701–717, 1980
1980
-
[38]
Tari, H.-J
H. Tari, H.-J. Su, T.-Y. Li. A constrained homotopy technique for excluding unwanted solutions from polynomial equations arising in kinematics problems. Mech. Mach. Theory, 45(6), 898–910, 2010
2010
-
[39]
Traverso
C. Traverso. Gr¨ obner trace algorithms. InInternational Symposium on Symbolic and Algebraic Compu- tation. Springer, Berlin, 1988
1988
-
[40]
Wampler, A.P
C.W. Wampler, A.P. Morgan, A.J. Sommese. Complete solution of the nine-point path synthesis problem for four-bar linkages. ASME J. Mech. Des. . 114(1), 153–161, 1992
1992
-
[41]
R. Zippel. Probabilistic algorithms for sparse polynomials. LNCS, 72, 216–226, 1979. Jonathan D. Hauenstein . Department of Applied and Computational Mathematics and Statistics, University of Notre Dame, Notre Dame, IN, USA. (hauenstein@nd.edu, www.nd.edu/~jhauenst). Martin He...
1979
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.