Pith. Sign in

REVIEW

New integer sequence OEIS A392714 counts Wronskians: fast evaluation via late-growing permutations

Number Theory · Submitted 6 Oct 2026 · 24 references

Pith's verdict: Genuinely new sequence with exact values to p=18 via a clean subset recursion; the math holds up, the abstract over-claims an unpublished growth bound, and the tail c(15), c(18) need an independent run.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on October 7, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 4 major objections · 6 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2610.08636 (version 1)
Pith IDT5U3K4LY, Pith's permanent record of this paper and its review
CitationArthemy V. Kiselev, Kian C. Shah. "New integer sequence OEIS A392714 counts Wronskians: fast evaluation via late-growing permutations." arXiv:2610.08636 (2026). https://pith.science/paper/T5U3K4LY
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

The alternating composition of 2p weighted differential operators on the line equals c(p) times the Wronskian of the weights, and the paper computes the universal integers c(p) exactly up to p = 18 as the new OEIS sequence A392714.

The reading

The paper establishes an exact bridge between differential operators and combinatorics: on the real line, the fully antisymmetrised composition of $N = 2p$ weighted operators $w_j(x)\,\partial_x^p$ acts on any smooth $f$ as $c(p)\cdot\mathrm{Wronsk}(w_1,\dots,w_N)(x)\,\partial_x^p f(x)$, where $c(p)$ is a universal integer that does not depend on $f$ or on the weights. The paper's contribution is to make that constant computable: feeding in the monomial weights $w_j = x^{j-1}$ reduces the sum over all $(2p)!$ permutations to a much smaller signed sum over the late-growing permutations $\Phi_p \subseteq S_{2p-1}$, expressed through falling factorials. A memoized recursion over subsets, run modulo several primes and reconstructed by the Chinese Remainder Theorem, then yields the exact values $c(p)$ for every $p \leq 18$, the largest being $c(18) \approx 4.881\times 10^{462}$, and these integers turn out to form a sequence new to the OEIS, registered as A392714. The constants are not a curiosity: they are the coefficients of the higher-arity brackets that deform the binary Lie bracket into strongly homotopy Lie algebras, and their super-exponential growth controls how strongly such collective terms can act.

What carries the argument

The load-bearing object is the late-growing permutation set $\Phi_p\subseteq S_{2p-1}$ together with the closed formula (8) that expresses $c(p)$ through it. Late-growing permutations are those whose partial sums never exceed the partial sums of the average permutation; after the monomial substitution they are exactly the permutations whose iterated action keeps the exponent of $x$ nonnegative, so they replace the full group $S_{2p}$ in the alternating sum. Formula (8) reads $c(p)\prod_{k=0}^{2p-1}k! = \sum_{\sigma\in\Phi_p}(-1)^\sigma\prod_{k=1}^{2p-1}(E_k)_p$, where $E_k$ is the running exponent built by the right-to-left operator action and $(E_k)_p$ is the falling factorial. The decisive acceleration is Lemma 4.2: the exponent $E(A) = p + (p+1)|A| - \sum_{j\in A} j$ depends only on the set $A$ of labels still to be placed, never on their order, so the sum obeys the subset recursion $\Psi(A) = \sum_{e\in A,\ E(A\setminus\{e\})\geq p}(-1)^{\#\{f\in A:\ f>e\}}(E(A\setminus\{e\}))_p\,\Psi(A\setminus\{e\})$ with $\Psi(\{2,\dots,2p\}) = c(p)\prod_{k=0}^{2p-1}k!$. Memoising this recursion over bit masks and sweeping it modulo primes, with a Chinese Remainder Theorem reconstruction, is Algorithm 2; the underlying identity (2) is supported by the appendix's argument that antisymmetry plus the derivative budget of $2p^2$ forces the unique surviving Leibniz profile $\{0,1,\dots,N-1\}$, making the alternating sum proportional to the Wronskian.

Extended reading notes

Core claim

The central claim is identity (2): for any smooth $f$ and any weights, $A_p[w_1,\dots,w_N] = c(p)\cdot\mathrm{Wronsk}(w_1,\dots,w_N)\cdot\partial_x^p$, where $A_p$ is the signed sum $\sum_{\sigma\in S_{2p}}(-1)^\sigma w_{\sigma(1)}\partial_x^p\circ\cdots\circ w_{\sigma(N)}\partial_x^p$ and the Wronskian is the determinant of the matrix of iterated derivatives of the weights. The paper computes the universal constants by specialising to the monomials $1,x,\dots,x^{2p-1}$: after that substitution a permutation contributes precisely when it is late-growing, and the constant is the signed sum $c(p)\prod_{k=0}^{2p-1}k! = \sum_{\sigma\in\Phi_p}(-1)^\sigma\prod_{k=1}^{2p-1}(E_k)_p$ of falling factorials $(E)_p = E!/(E-p)! = E(E-1)\cdots(E-p+1)$. Evaluating this formula gives $c(1)=1$, $c(2)=2$, $c(3)=90$, $c(4)=586656$, and so on up to $c(18)\approx 4.881\times10^{462}$, a set of values that matches no existing entry in the OEIS and is therefore registered there as the new sequence A392714. The growth is $\log c(p) \simeq 2p^2\log p - bp^2 + o(p^2)$ with subleading constant $b\geq 2.6744$, and the algorithm that reaches $p=18$ exploits the fact that the recursion's state is just the set of not-yet-placed labels, stored as bit masks, with all arithmetic done modulo primes.

Load-bearing premise

Everything computed inherits from the appendix's assertion, in proving Theorem 1.1, that the coefficient of the unique surviving Leibniz profile {0,1,...,N-1} is a universal integer depending only on p and not on the test function or the weight functions; the appendix shows the profile is unique but does not work out the coefficient's universality in detail.

Editorial extensions

If this is right

  • The exact integers c(p) for 1 ≤ p ≤ 18 are now the registered sequence A392714, and no existing OEIS entry matched the computed values, so the sequence is new to the registry.
  • The fitted growth law predicts that c(19) through c(22) have about 532, 606, 685, and 770 decimal digits, so running the same algorithm on larger machines should produce the next exact values.
  • Since log c(p) ~ 2p^2 log p, the constants outrun the factorial (2p)!; in the strongly homotopy Lie models that motivate identity (2), that super-exponential growth means the collective N-ary brackets at arity N > 4 are numerically suppressed relative to the binary Lie bracket.
  • The prime factorisations follow a visible pattern — small primes in high powers, very few very large primes — and v_p(c(p)) = p − 1 for every prime p ≤ 17, supporting Conjecture 5.1; the factorisation of c(14) yields a newly confirmed 206-decimal-digit prime.
  • The counts of even and odd permutations in every Phi_p differ by exactly (−1)^(p+1), a parity law the paper states and announces a proof of in a sequel.

Reading between the lines

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

  • Extension: because the exponent E(A) in Lemma 4.2 depends on A only through |A| and the sum of its labels, the memoized subset table may collapse into a two-parameter recurrence that produces c(p) without enumerating any subsets; the paper does not pursue this.
  • Extension: the same template — monomial weights, falling-factorial sums, subset-state recursion, modular reconstruction — should transfer to the multidimensional Wronskian brackets the authors list as an open problem, yielding constants c(p_1,...,p_d) indexed by per-variable differential orders.
  • Extension: the next term is a cheap joint test — running Algorithm 2 for p = 19 should produce a value of about 532 digits whose 19-adic valuation would be 18 if Conjecture 5.1 holds; a mismatch would expose either the algorithm or the valuation pattern.
Original abstract

The alternating composition of $N = 2p$ weighted differential operators $w_j(x)\cdot\partial_x^{\,p}$ of strict order $p$ on the line $\mathbb{R} \ni x$ is again an operator of order $p$; its coefficient is the universal constant $c(p)$ times the Wronskian of the weights $w_1,\ldots,w_N$. Lie brackets of vector fields fix $c(p=1)=1$; we want to find $c(p \geqslant 2)$: e.g., $c(2) = 2$ or $c(3) = 90$. Direct symbolic expansion (over $|S_{2p}| =(2p)!$ permutations) fails for $p \geqslant 4$. Taking the monomials $w_j = x^{j-1}$ reduces the summation to the much smaller set $\Phi_p \subseteq S_{2p-1} \subsetneq S_{2p}$ of late-growing permutations. Expressing $c(p)$ as a signed sum of products of falling factorials, we implement and speed up the algorithm that gains all the integer values up to $c(18) = 4.881\ldots \cdot 10^{462}$. The resulting sequence is new, now registered as OEIS A392714; its (sub)leading-order growth rate is $\log c(p) \simeq 2p^2\log p -b p^2 + \overline{o}(p^2)$ for $p\gg 1$, with $b\geqslant 2.6744$.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-10-07: 1 request, 5,499 input tokens.

What the abstract itself states, as probabilities: claims first 0.35 · states limitation 0.80 · human participants 0.01 · controlled comparison 0.03 · randomized assignment 0.01 · states quantitative result 0.98 · compares with prior results 0.74 · mentions code or data availability 0.17.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

The alternating composition of $N = 2p$ weighted differential operators $w_j(x)\cdot\partial_x^{\,p}$ of strict order $p$ on the line $\mathbb{R} \ni x$ is again an operator of order $p$; its coefficient is the universal constant $c(p)$ times the Wronskian of the weights $w_1,\ldots,w_N$. Lie brackets of vector fields fix $c(p=1)=1$; we want to find $c(p \geqslant 2)$: e.g., $c(2) = 2$ or $c(3) = 90$. Direct symbolic expansion (over $|S_{2p}| =(2p)!$ permutations) fails for $p \geqslant 4$. Taking the monomials $w_j = x^{j-1}$ reduces the summation to the much smaller set $\Phi_p \subseteq S_{2p-1} \subsetneq S_{2p}$ of late-growing permutations. Expressing $c(p)$ as a signed sum of products of falling factorials, we implement and speed up the algorithm that gains all the integer values up to $c(18) = 4.881\ldots \cdot 10^{462}$. The resulting sequence is new, now registered as OEIS A392714; its (sub)leading-order growth rate is $\log c(p) \simeq 2p^2\log p -b p^2 + \overline{o}(p^2)$ for $p\gg 1$, with $b\geqslant 2.6744$.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 The alternating composition of $N = 2p$ weighted differential operators $w_j(x)\cdot\partial_x^{\,p}$ of strict order $p$ on the line $\mathbb{R} \ni x$ is again an operator of order $p$; its coefficient is the universal constant $c(p)$ times the Wronskian of the weights $w_1,\ldots,w_N$. background 0.71 0.39 0.29 asserted 1.00
s1 Lie brackets of vector fields fix $c(p=1)=1$; we want to find $c(p \geqslant 2)$: e.g., $c(2) = 2$ or $c(3) = 90$. gap 0.70 0.44 0.07 asserted 0.96
s2 Direct symbolic expansion (over $|S_{2p}| =(2p)!$ permutations) fails for $p \geqslant 4$. limitation 0.52 0.60 0.03 asserted 1.00
s3 Taking the monomials $w_j = x^{j-1}$ reduces the summation to the much smaller set $\Phi_p \subseteq S_{2p-1} \subsetneq S_{2p}$ of late-growing permutations. method 0.97 0.81 0.07 asserted 1.00
s4 Expressing $c(p)$ as a signed sum of products of falling factorials, we implement and speed up the algorithm that gains all the integer values up to $c(18) = 4.881\ldots \cdot 10^{462}$. method 0.50 0.95 0.17 asserted 1.00
s5 The resulting sequence is new, now registered as OEIS A392714; its (sub)leading-order growth rate is $\log c(p) \simeq 2p^2\log p -b p^2 + \overline{o}(p^2)$ for $p\gg 1$, with $b\geqslant 2.6744$. result 0.91 0.92 0.63 asserted 0.98

What kind of paper the classifier reports this is: method or tool 0.85 · evidence from simulation or computation 0.79.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.