Pith. sign in

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 →

arxiv 1908.06020 v1 pith:G2JT7M7W submitted 2019-08-15 math.AG cs.SC

classification math.AGcs.SC MSC 14Q2013P1068W30
keywords Alt'sproblemfour-barlinkagesprobabilisticsaturationGröbnerbasesfinitefieldsHilbertfunctionsNewton-OkounkovbodiesChernclasses
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

Alt's 1923 problem asks how many four-bar linkages trace a coupler curve through nine prescribed points. A 1992 homotopy continuation computation found 8,652 isolated solutions, i.e. 4,326 distinct linkages, but that computation was numerical and could not rule out additional solutions. This paper tries to turn the count into a sharp upper bound by replacing the coefficients of the coupler-curve equations with independent random parameters and counting, outside the degenerate base locus $X$, the solutions to a system of general linear combinations over finite fields. The reported finite-field computations yield $g_0(X,\mathbb{C}^8)^{\mathrm{rand}}_{\mathbb{Z}_p}=8652$ for a range of primes, so the quotient dimension is an upper bound on the true complex degree; dividing by the relabeling symmetry gives at most 4,326 distinct four-bar linkages, and dividing by the full symmetry gives at most 1,442 coupler curves.

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.

Watch

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

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

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)
  1. [§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.
  2. [§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.
  3. [§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)
  1. [§4.1, Example 4.6] The word 'Symoblically' should be 'Symbolically'.
  2. [§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.
  3. [Definition 4.4] The phrase 'p is larger than any coefficient' should specify that absolute values are meant, since coefficients may be negative.
  4. [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

0 steps flagged · score 0.0 of 10

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 0 free parameters · 4 assumptions · 0 invented entities

No fitted constants are introduced: the entries of Theta, Lambda, and mu are random inputs, and p is a user-chosen prime. The derivation relies on standard algebraic geometry and Groebner basis facts, plus the false injectivity assumption in Theorem 4.7 that breaks the advertised probability guarantees.

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).
    Used in Theorem 4.1; the proof states it follows from Bertini's Theorem and the Rabinowitz trick as in [24, Thm. 4.1].
  • 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.
    Basis of the modular computations; cited from [33, Prop. 4.1] and [4].
  • standard math Schwartz-Zippel lemma bounds the probability that a bounded-degree polynomial vanishes on random inputs.
    Used in Proposition 4.13 to convert degree bounds on the discriminant into probability bounds.
  • ad hoc to paper Nonconstant polynomials over F_p are injective, so a leading coefficient is nonzero with probability at least (p-1)/p.
    Assumed in the proof of Theorem 4.7; this is false, so the theorem's probability bound does not follow.

how reviews work

0 comments
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

Figures reproduced from arXiv: 1908.06020 by the authors.

Figure 1
Figure 1. A four-bar linkage with corresponding coupler curve 2.1. Relation to Alt’s problem. A four-bar linkage is the simplest moveable planar closed￾chain linkage, one of which is shown in [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

41 extracted references · 41 canonical work pages

  1. [1]

    Abbott, M

    J. Abbott, M. Kreuzer, L. Robbiano. Computing zero-dimensional schemes. J. Symb. Comput. , 39(1), 31–49, 2005

  2. [2]

    P. Aluffi. Computing characteristic classes of projective schemes. J. Symb. Comput. , 35(1), 3–19, 2003

  3. [3]

    H. Alt. ¨Uber die Erzeugung gegebener ebener Kurven mit Hilfe des Gelenkvierecks. ZAMM, 3(1), 13–19, 1923

  4. [4]

    E.A. Arnold. Modular algorithms for computing Gr¨ obner bases.J. Symb. Comput., 35(4), 403–419, 2003

  5. [5]

    Baskar, S

    A. Baskar, S. Bandyopadhyay. An algorithm to compute the finite roots of large systems of polynomial equations arising in kinematic synthesis. Mech. Mach. Theory, 133, 493–513, 2019

  6. [6]

    Bates, W

    D.J. Bates, W. Decker, J.D. Hauenstein, C. Peterson, G. Pfister, F.-O. Schreyer, A.J. Sommese, C.W. Wampler. Comparison of probabilistic algorithms for analyzing the components of an affine algebraic variety. Appl. Math. Comput. , 231, 619–633, 2014

  7. [7]

    Bates, J.D

    D.J. Bates, J.D. Hauenstein, A.J. Sommese, C.W. Wampler. Bertini: Software for Numerical Algebraic Geometry. Available at bertini.nd.edu

  8. [8]

    Bates, J.D

    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

Show all 41 references
  1. [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

  2. [10]

    L. Blum, F. Cucker, M. Shub, S. Smale. Complexity and real computation , Springer-Verlag, New York, 1998

  3. [11]

    Bosma, J

    W. Bosma, J. Cannon, C. Playoust. The Magma algebra system. I. The user language.J. Symb. Comput., 24, 235–265, 1997

  4. [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

  5. [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)

  6. [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

  7. [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

  8. [16]

    Grayson, M

    D. Grayson, M. Stillman. Macaulay2, a software system for research in algebraic geometry , available at www.math.uiuc.edu/Macaulay2/

  9. [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

  10. [18]

    G¨ ortz, T

    U. G¨ ortz, T. Wedhorn.Algebraic Geometry: Part I: Schemes . Springer Science & Business Media, 2010

  11. [19]

    Harris, M

    C. Harris, M. Helmer. Segre class computation and practical applications. To appear in Math. Comput

  12. [20]

    J. Harris. Algebraic geometry: a first course . Springer Science & Business Media, 2013

  13. [21]

    Hauenstein, J.I

    J.D. Hauenstein, J.I. Rodriguez. Multiprojective witness sets and a trace test. To appear in Adv. Geom

  14. [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

  15. [23]

    Hauenstein, F

    J.D. Hauenstein, F. Sottile. Algorithm 921: alphaCertified: Certifying solutions to polynomial systems. ACM Trans. Math. Softw. , 38(4), 28, 2012

  16. [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

  17. [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

  18. [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

  19. [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

  20. [28]

    Leykin, J

    A. Leykin, J. Yu. Beyond polyhedral homotopies. J. Symb. Comput. , 91, 173–190, 2019

  21. [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

  22. [30]

    Morgan, A.J

    A.P. Morgan, A.J. Sommese. Coefficient-parameter polynomial continuation. Appl. Math. Comput. , 29(2), 123–160, 1989

  23. [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

  24. [32]

    Okounkov

    A. Okounkov. Brunn-Minkowski inequality for multiplicities. Invent. Math., 125(3), 405–411, 1996

  25. [33]

    F. Pauer. On lucky ideals for Gr¨ obner basis computations. J. Symb. Comput. , 14, 471–482, 1992

  26. [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

  27. [35]

    S. Roberts. On three-bar motion in plane space. Proc. London Mathematical Society, III, 286–319, 1875

  28. [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

  29. [37]

    Schwartz

    J.T. Schwartz. Fast probabilistic algorithms for verification of polynomial identities. J. ACM , 27(4), 701–717, 1980

  30. [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

  31. [39]

    Traverso

    C. Traverso. Gr¨ obner trace algorithms. InInternational Symposium on Symbolic and Algebraic Compu- tation. Springer, Berlin, 1988

  32. [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

  33. [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...

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.