viernes, 11 de septiembre de 2026

Complex-Exponential Coordinates for Triangles and Quadrilaterals

The present note combines two earlier strands of work by the present author. The USM paper supplies the principal-branch secant–sine identity and its boundary-value conventions, while the earlier MATINF article, “Two Identities and their Consequences” supplies the generalized opposite-angle half-angle identities for a convex quadrilateral. The new point here is to organize these ingredients through oriented complex phases, giving a multiplicative triangle closure and a two-phase quadrilateral identity whose real, imaginary, and modulus statements recover classical metric formulas.


1. USM input and the two exponential regimes

We adopt the principal-branch and boundary-value conventions of the USM paper. In particular,

$$\operatorname{arcsec}z:=\arccos(1/z),\qquad z\ne0,$$

with the principal branch of $\arccos$ and the boundary-value convention used there. The part of Theorem 3.2 needed here states that, for $0\lt y\le1$,

$$e^{i\operatorname{arcsec}y}=\tan\left(\frac12\arcsin y\right).\tag{1}$$

For $0\lt y\lt1$, the boundary value of $\operatorname{arcsec}y$ is purely imaginary; consequently, the left-hand side of (1), although written as a complex exponential, is real and positive.

Let $\theta\in(0,\pi)$. Since $\arcsin(\sin\theta)=\min\{\theta,\pi-\theta\}$, setting $y=\sin\theta$ in (1) gives the first coordinate.


Proposition 1.1 (Real-valued exponential coordinate). For $\theta\in(0,\pi)$,

$$\boxed{E(\theta):=e^{i\operatorname{arcsec}(\sin\theta)}=\tan\left(\frac{\min\{\theta,\pi-\theta\}}{2}\right)=\frac{\sin\theta}{1+|\cos\theta|}.}\tag{2}$$

Hence $E(\theta)\in(0,1]$.

To obtain a genuinely complex number from real Euclidean data, use the reciprocal argument $\csc\theta\ge1$. Define $\varepsilon_\theta=1$ for $0\lt\theta\le\pi/2$ and $\varepsilon_\theta=-1$ for $\pi/2\lt\theta\lt\pi$. Since

$$\operatorname{arcsec}(\csc\theta)=\arccos(\sin\theta)=\left|\frac{\pi}{2}-\theta\right|,$$

we obtain an oriented phase.


Proposition 1.2 (Oriented complex phase). For $\theta\in(0,\pi)$,

$$\boxed{\Phi(\theta):=e^{i\varepsilon_\theta\operatorname{arcsec}(\csc\theta)}=e^{i(\pi/2-\theta)}=\sin\theta+i\cos\theta.}\tag{3}$$

The phase has modulus $1$ and is nonreal unless $\theta=\pi/2$.

The two coordinates therefore behave in complementary ways: $E(\theta)$ uses a complex inverse secant but produces a real value, whereas $\Phi(\theta)$ uses a real inverse secant and produces a genuine complex phase.


2. Triangles

Let $ABC$ be a nondegenerate triangle with side lengths $a,b,c$ opposite $A,B,C$, respectively, semiperimeter

$$s=\frac{a+b+c}{2},$$

and area $\Delta\gt0$.

The standard metric relations

$$\sin A=\frac{2\Delta}{bc},\qquad \cos A=\frac{b^2+c^2-a^2}{2bc}\tag{4}$$

give a side-and-area form of the phase.


Proposition 2.1 (Triangle phase coordinate). Define $\Phi_A:=\Phi(A)$ and cyclically. Then

$$\boxed{\Phi_A=e^{i(\pi/2-A)}=\frac{4\Delta+i(b^2+c^2-a^2)}{2bc},}\tag{5}$$

with the analogous formulas for $\Phi_B$ and $\Phi_C$. Moreover,

$$\boxed{\Phi_A\Phi_B\Phi_C=i.}\tag{6}$$


Proof. Equation (5) follows by substituting (4) into (3). Since $A+B+C=\pi$,

$$\Phi_A\Phi_B\Phi_C=e^{i[(\pi/2-A)+(\pi/2-B)+(\pi/2-C)]}=e^{i\pi/2}=i.$$


2.1. A right-triangle specialization

Suppose $a^2+b^2=c^2$ and let $B$ be the acute angle opposite $b$. Then $\sin B=b/c$ and $\cos B=a/c$. Proposition 1.1 therefore gives

$$\boxed{e^{i\operatorname{arcsec}(b/c)}=\frac{b/c}{1+a/c}=\frac{b}{a+c}=\frac{c-a}{b}.}\tag{7}$$

The last equality follows from $b^2=c^2-a^2=(c-a)(c+a)$. Thus the real-valued USM exponential becomes a purely rational expression in the side lengths.


2.2. Heron’s formula from the phase product

Set

$$X_A=b^2+c^2-a^2,\qquad X_B=c^2+a^2-b^2,\qquad X_C=a^2+b^2-c^2.\tag{8}$$

Substitution of the metric forms into (6) yields

$$(4\Delta+iX_A)(4\Delta+iX_B)(4\Delta+iX_C)=8ia^2b^2c^2.\tag{9}$$

The right-hand side has zero real part. Hence

$$(4\Delta)^3-4\Delta(X_AX_B+X_BX_C+X_CX_A)=0.$$

Since $\Delta\gt0$,

$$16\Delta^2=X_AX_B+X_BX_C+X_CX_A.\tag{10}$$

A direct expansion gives

$$X_AX_B+X_BX_C+X_CX_A=2a^2b^2+2b^2c^2+2c^2a^2-a^4-b^4-c^4,$$

and

$$2a^2b^2+2b^2c^2+2c^2a^2-a^4-b^4-c^4=(a+b+c)(-a+b+c)(a-b+c)(a+b-c).$$

Since the four factors are $2s$, $2(s-a)$, $2(s-b)$, and $2(s-c)$, respectively, (10) gives:


Theorem 2.2 (Heron).

$$\boxed{\Delta^2=s(s-a)(s-b)(s-c).}\tag{11}$$

The role of the exponential in this proof is explicit: the additive closure $A+B+C=\pi$ becomes the multiplicative closure (6), while each phase packages the sine and cosine metric data of one vertex into a single complex number.


2.3. Comparison with Edwards’s complex proof

A different complex-number proof of Heron’s formula was published by Miles Dillon Edwards in 2007; an online reproduction of Edwards’s proof is also available. It is worth distinguishing that argument from the phase proof above, because the two share the same broad mechanism—turning an angle-sum relation into multiplication in $\mathbb C$—but use different geometric coordinates.

In Edwards’s construction, let $I$ be the incenter, let $r$ be the inradius, and write the tangent-length decomposition

$$a=y+z,\qquad b=x+z,\qquad c=x+y,$$

so that

$$s=x+y+z,\qquad x=s-a,\qquad y=s-b,\qquad z=s-c.$$

The three right triangles formed by the inradius give angles $\alpha,\beta,\gamma$ at $I$ satisfying

$$\alpha+\beta+\gamma=\pi.$$

The complex numbers $r+ix$, $r+iy$, and $r+iz$ have arguments $\alpha$, $\beta$, and $\gamma$, respectively. Their product is therefore negative real, and its imaginary part vanishes:

$$0=\operatorname{Im}\left[(r+ix)(r+iy)(r+iz)\right]=r^2(x+y+z)-xyz.\tag{12}$$

Thus

$$r^2=\frac{xyz}{s}=\frac{(s-a)(s-b)(s-c)}{s},$$

and $\Delta=rs$ gives Heron’s formula.

The contrast with the present proof is precise. Edwards uses non-unit complex numbers attached to the incenter right triangles and extracts an imaginary part because their product lies on the real axis. Here the coordinates are unit phases attached directly to the vertex angles, and after metric substitution one extracts a real part because the phase product equals $i$. The present construction does not use the incircle; more importantly for the purpose of this article, the same phase language extends to the general quadrilateral identity of Section 3 and its modulus form in Section 4.


3. General convex quadrilaterals

Let $ABCD$ be a convex quadrilateral with

$$AB=a,\qquad BC=b,\qquad CD=c,\qquad DA=d,$$

semiperimeter

$$s=\frac{a+b+c+d}{2},$$

area $K\gt0$, and opposite angles

$$\alpha=\angle BAD,\qquad \gamma=\angle BCD.$$

Define

$$P:=(s-a)(s-d),\qquad Q:=(s-b)(s-c).\tag{13}$$


3.1. The general half-angle identities

The following pair was established by the present author in “Two Identities and their Consequences”, Theorem 5, as a generalization of the cyclic half-angle identities in Theorem 1 of that article. We include a short proof for completeness.

Lemma 3.1. For every convex quadrilateral as above,

$$ad\sin^2\frac{\alpha}{2}+bc\cos^2\frac{\gamma}{2}=P,\tag{14}$$

and

$$bc\sin^2\frac{\gamma}{2}+ad\cos^2\frac{\alpha}{2}=Q.\tag{15}$$


Proof. Let $p=BD$. The law of cosines in triangles $ABD$ and $BCD$ gives

$$p^2=a^2+d^2-2ad\cos\alpha=b^2+c^2-2bc\cos\gamma,$$

so

$$ad\cos\alpha-bc\cos\gamma=\frac{a^2+d^2-b^2-c^2}{2}.\tag{16}$$

Using the half-angle formulas,

$$ad\sin^2\frac{\alpha}{2}+bc\cos^2\frac{\gamma}{2}=\frac{ad+bc+bc\cos\gamma-ad\cos\alpha}{2},$$

hence

$$ad\sin^2\frac{\alpha}{2}+bc\cos^2\frac{\gamma}{2}=\frac{ad+bc}{2}+\frac{b^2+c^2-a^2-d^2}{4}=(s-a)(s-d)=P.$$

The second identity is obtained similarly, or by subtracting the first from $ad+bc=P+Q$.

Two elementary consequences of (13) are

$$P+Q=ad+bc,\qquad P-Q=\frac{b^2+c^2-a^2-d^2}{2}.\tag{17}$$

Combining the second relation with (16) gives the useful equivalent form

$$Q-P=ad\cos\alpha-bc\cos\gamma.\tag{18}$$

This is the bridge from the generalized half-angle identities to the imaginary part of the complex relation below.


3.2. The general complex-exponential identity

Define the oriented phases

$$\Phi_\alpha:=\Phi(\alpha)=e^{i(\pi/2-\alpha)},\qquad \Phi_\gamma:=\Phi(\gamma)=e^{i(\pi/2-\gamma)}.$$

Thus

$$\Phi_\alpha=\sin\alpha+i\cos\alpha,\qquad \overline{\Phi_\gamma}=\sin\gamma-i\cos\gamma.$$

The decomposition of the quadrilateral along $BD$ supplies the real part:

$$2K=ad\sin\alpha+bc\sin\gamma.\tag{19}$$

The generalized half-angle identities, through (18), supply the imaginary part. Combining the two gives the central identity.


Theorem 3.2. For every convex quadrilateral $ABCD$,

$$\boxed{ad\,\Phi_\alpha+bc\,\overline{\Phi_\gamma}=2K+i(Q-P)=2K+\frac{i}{2}(a^2+d^2-b^2-c^2).}\tag{20}$$

Equivalently,

$$\boxed{ad\,e^{i(\pi/2-\alpha)}+bc\,e^{-i(\pi/2-\gamma)}=2K+\frac{i}{2}(a^2+d^2-b^2-c^2).}\tag{21}$$

Using the principal inverse secant, the phase factors appearing in (20) can also be written as

$$\Phi_\alpha=e^{i\varepsilon_\alpha\operatorname{arcsec}(\csc\alpha)},\qquad \overline{\Phi_\gamma}=e^{-i\varepsilon_\gamma\operatorname{arcsec}(\csc\gamma)},$$

with $\varepsilon$ defined as in Proposition 1.2.


Proof. By definition,

$$ad\Phi_\alpha+bc\overline{\Phi_\gamma}=ad(\sin\alpha+i\cos\alpha)+bc(\sin\gamma-i\cos\gamma),$$

so

$$ad\Phi_\alpha+bc\overline{\Phi_\gamma}=(ad\sin\alpha+bc\sin\gamma)+i(ad\cos\alpha-bc\cos\gamma).$$

Now apply (19) and (18).

Thus (20) packages two real metric identities into one complex relation:

$$\operatorname{Re}:\qquad 2K=ad\sin\alpha+bc\sin\gamma,$$

and

$$\operatorname{Im}:\qquad Q-P=ad\cos\alpha-bc\cos\gamma=\frac{a^2+d^2-b^2-c^2}{2}.$$


4. Bretschneider from the modulus

In the earlier MATINF article, Bretschneider’s formula was obtained by multiplying the two generalized half-angle identities and rearranging the result (Theorem 6). The calculation below packages the same underlying metric data into the modulus of the single complex relation (20).

Taking the squared modulus of (20) gives

$$4K^2+\frac14(a^2+d^2-b^2-c^2)^2=|ad\Phi_\alpha+bc\overline{\Phi_\gamma}|^2,$$

hence

$$4K^2+\frac14(a^2+d^2-b^2-c^2)^2=a^2d^2+b^2c^2+2abcd\operatorname{Re}(\Phi_\alpha\Phi_\gamma),$$

and therefore

$$4K^2+\frac14(a^2+d^2-b^2-c^2)^2=a^2d^2+b^2c^2-2abcd\cos(\alpha+\gamma).\tag{22}$$

The last equality uses

$$\Phi_\alpha\Phi_\gamma=e^{i[\pi-(\alpha+\gamma)]}.$$

From (17),

$$4PQ=(ad+bc)^2-\frac14(a^2+d^2-b^2-c^2)^2.\tag{23}$$

Substituting (23) into (22), and using

$$1+\cos(\alpha+\gamma)=2\cos^2\frac{\alpha+\gamma}{2},$$

yields the classical general area formula.


Theorem 4.1 (Bretschneider). For every convex quadrilateral,

$$\boxed{K^2=(s-a)(s-b)(s-c)(s-d)-abcd\cos^2\left(\frac{\alpha+\gamma}{2}\right).}\tag{24}$$

Thus Bretschneider’s formula is the modulus statement associated with the complex-exponential identity (20). In this sense, the modulus proof is a complex-exponential reformulation of the half-angle route used in the earlier MATINF paper.


5. Particular cases

5.1. Cyclic quadrilateral: Brahmagupta

If $ABCD$ is cyclic, then

$$\alpha+\gamma=\pi.$$

Therefore

$$\cos^2\left(\frac{\alpha+\gamma}{2}\right)=0,$$

and (24) immediately becomes:

Corollary 5.1 (Brahmagupta).

$$\boxed{K^2=(s-a)(s-b)(s-c)(s-d).}\tag{25}$$

Cyclicity makes the two phase factors appearing in (20) coincide. Since $\gamma=\pi-\alpha$,

$$\overline{\Phi_\gamma}=\Phi_\alpha,$$

so (20) becomes

$$(ad+bc)\Phi_\alpha=2K+\frac{i}{2}(a^2+d^2-b^2-c^2),$$

or

$$\boxed{e^{i(\pi/2-\alpha)}=\frac{4K+i(a^2+d^2-b^2-c^2)}{2(ad+bc)}.}\tag{26}$$

This is the cyclic quadrilateral analogue of the triangle phase coordinate.

The general half-angle identities also collapse to the familiar cyclic pair. Indeed,

$$\cos^2\frac{\gamma}{2}=\sin^2\frac{\alpha}{2},\qquad \sin^2\frac{\gamma}{2}=\cos^2\frac{\alpha}{2},$$

so (14)–(15) give

$$\boxed{\sin^2\frac{\alpha}{2}=\frac{(s-a)(s-d)}{ad+bc},\qquad \cos^2\frac{\alpha}{2}=\frac{(s-b)(s-c)}{ad+bc}.}\tag{27}$$

These are precisely the cyclic identities used as the starting point in the earlier MATINF article, Theorem 1.


5.2. Bicentric quadrilateral

If the cyclic quadrilateral is also tangential, Pitot’s relation gives

$$a+c=b+d=s.$$

Hence

$$s-a=c,\qquad s-b=d,\qquad s-c=a,\qquad s-d=b.$$

Brahmagupta therefore reduces to:

Corollary 5.2 (Bicentric area formula).

$$\boxed{K=\sqrt{abcd}.}\tag{28}$$

The same area formula was treated in the earlier MATINF paper, Theorem 4; here it appears as an immediate specialization of the Brahmagupta branch of the complex-exponential framework.


5.3. Heron as a degenerate cyclic limit

Let one side of a cyclic quadrilateral tend to zero, say $d\to0^+$, while the remaining sides tend to those of a nondegenerate triangle. Then

$$s\longrightarrow\frac{a+b+c}{2},$$

and Brahmagupta’s formula tends to

$$K^2=s(s-a)(s-b)(s-c),$$

which is Heron’s formula. This limiting route is complementary to the independent triangle phase proof in Section 2. It is a degenerate limit rather than a specialization within the class of convex quadrilaterals, since the limiting side length is zero.


5.4. What cyclicity means in phase language

For a general quadrilateral, the complex identity (20) contains the two opposite-angle phase factors $\Phi_\alpha$ and $\overline{\Phi_\gamma}$. Cyclicity is exactly the condition

$$\alpha+\gamma=\pi\quad\Longleftrightarrow\quad\overline{\Phi_\gamma}=\Phi_\alpha.$$

Thus a cyclic quadrilateral is precisely the case in which the two phase factors appearing in the general identity collapse to one. In this sense, Brahmagupta’s formula is the one-phase specialization of the two-phase Bretschneider identity.


6. Summary and outlook

1. For $\theta\in(0,\pi)$, the USM secant–sine identity provides the real-valued exponential coordinate $E(\theta)$.

2. Passing to the reciprocal argument and orienting the principal inverse secant gives the unit phase

$$\Phi(\theta)=e^{i(\pi/2-\theta)}=\sin\theta+i\cos\theta.$$

3. For triangles, the closure $A+B+C=\pi$ becomes $\Phi_A\Phi_B\Phi_C=i$; the real part of the resulting metric product yields Heron’s formula. This phase proof is structurally related to, but geometrically distinct from, Edwards’s earlier incenter-based complex proof.

4. For general convex quadrilaterals, the generalized opposite-angle half-angle identities from the earlier MATINF paper give $Q-P=ad\cos\alpha-bc\cos\gamma$. Together with the area decomposition, this yields the two-phase relation

$$ad\Phi_\alpha+bc\overline{\Phi_\gamma}=2K+i(Q-P)=2K+\frac{i}{2}(a^2+d^2-b^2-c^2).$$

Its modulus yields Bretschneider’s formula.

5. Under cyclicity, the two phase factors $\Phi_\alpha$ and $\overline{\Phi_\gamma}$ coincide, Bretschneider reduces to Brahmagupta, and the general half-angle formulas reduce to the classical cyclic half-angle identities. The bicentric formula follows as a specialization, while Heron is recovered again as a degenerate cyclic limit.

sábado, 5 de septiembre de 2026

The Ubiquitous Half-Angle Formulas and the Riemann Hypothesis

 For some time I have been studying the recurring role of half-angle formulas in geometry, rational parametrization, and integration. Some of that background is collected in The Ubiquitous Half-Angle Formulas (PDF link). Henning Dathe’s Die Halbwinkelsubstitution und ihre Anwendungen (PDF link) gives a broad treatment of the tangent half-angle substitution and its applications. A related paper, A Unified Substitution Method for Integration (see arXiv:2505.03754) (see also the trigonometric/hyperbolic formulas for the roots of quadratic equations that gave rise to the core identities of this method, and how the tangent of half-angle emerges here as well!), develops a branch-conscious framework in which several classical substitutions and identities arise from a common pair of principal-branch exponential identities. It recovers Euler’s first and second substitutions, after the necessary scaling, sign choices, and component-dependent reciprocal reparametrizations and obtains the classical Weierstrass substitution $t=\tan(\omega/2)$ as a local unit-radius corollary of Transform 5. The same identities also contain two notable geometric specializations. For the Cayley/tangent-half-angle form, take the positive hyperbolic parametrization $y=\cosh u$. Then $\arccos(\cosh u)=iu$ and

$$\tan\!\left(\tfrac12\operatorname{arcsec}(\cosh u)\right)=\tanh(u/2),$$

so the exterior identity becomes

$$e^{-u}=\frac{1-\tanh(u/2)}{1+\tanh(u/2)}.$$

Continuing by $u=i\phi$, with $\tanh(i\phi/2)=i\tan(\phi/2)$, gives the classical Cayley half-angle representation

$$e^{-i\phi}=\frac{1-i\tan(\phi/2)}{1+i\tan(\phi/2)},\qquad e^{i\phi}=\frac{1+i\tan(\phi/2)}{1-i\tan(\phi/2)}.$$

This follows directly from Theorems 3.1–3.2 of the arXiv draft and the branch computation used in §4; the derivation is isolated explicitly in this Cayley-type half-angle note. For Lobachevsky’s angle-of-parallelism formula, the same positive hyperbolic parametrization $y=\cosh u$ in Theorem 3.1 identifies $\beta=\operatorname{arccsc}(\cosh u)=\Pi(u)$ and yields

$$\tan\!\left(\frac{\Pi(u)}2\right)=e^{-u};$$

equivalently, Theorem 3.2 gives the same result using the reciprocal parametrization $y=\operatorname{sech} u$, for which $\psi=\arcsin(\operatorname{sech} u)=\Pi(u)$. Both derivations are given explicitly in §4, “Lobachevsky’s formula as a geometric specialization,” of the arXiv draft.

The transforms also simultaneously rationalize the associated quadratic radicals and inverse-trigonometric half-angle composites, and for several natural families cancellation between the parametrized radical and the Jacobian reduces the transformed integrand all the way to a finite Laurent polynomial, allowing term-by-term integration (in particular, some quotient-radical integrals collapse directly to Laurent polynomials and require no partial-fraction decomposition; see the arXiv draft, especially Example 4 and the termwise-integration remark). 

Recent benchmarks against Mathematica’s native `Integrate` show the strongest gains in the most difficult half-angle families: on mutually solved Transform 1 cases, native `Integrate` was about $194\times$ slower in cold timing and produced expressions about $753\times$ larger on average; for Transform 4 it was about $110\times$ slower cold, with mean closed-form output more than $100\times$ larger and one native result reaching roughly $2.96$ MB. Across the five benchmark families, USM solved and verified every tested example in Transforms 2–5 and produced 95 non-timeout results out of 100 in Transform 1. 

The point of this background is not that any of it is needed for the argument below, but rather to illustrate the remarkable range and structural power of half-angle formulas. Against that background, it is especially intriguing to encounter $\tan(\theta/2)$ once again, this time in a formulation equivalent to the Riemann Hypothesis. That recurrence is what made the observation below worth recording.

I want to emphasize that I am not claiming a proof of the Riemann Hypothesis. The RH argument below was generated by ChatGPT 5.6 (OpenAI) in response to my prompts (so take it with a grain of salt).

Let

$$\xi(s)=\frac12s(s-1)\pi^{-s/2}\Gamma(s/2)\zeta(s),$$

and let $\lambda_n$ denote the Li coefficients. Li’s criterion states

$$\mathrm{RH}\iff \lambda_n\ge 0\qquad\text{for every }n\ge1.$$

The Li coefficients also have, in the standard symmetric limiting sense, the zero-sum representation

$$\lambda_n=\sum_\rho\left[1-\left(1-\frac1\rho\right)^n\right].$$

Consider the Möbius transformation

$$t(s)=1-\frac1s=\frac{s-1}{s}.$$

For finite $s$,

$$\operatorname{Re}(s)=\frac12 \iff |t(s)|=1.$$

Thus, assuming RH and writing

$$\rho=\frac12+i\gamma,\qquad \gamma>0,$$

we may write

$$1-\frac1\rho=e^{i\theta_\gamma},\qquad 0<\theta_\gamma<\pi.$$

A direct calculation gives

$$\cos\theta_\gamma=\frac{4\gamma^2-1}{4\gamma^2+1},\qquad\sin\theta_\gamma=\frac{4\gamma}{4\gamma^2+1},$$

and hence

$$\tan\frac{\theta_\gamma}{2}=\frac{\sin\theta_\gamma}{1+\cos\theta_\gamma}=\frac1{2\gamma}.$$

Equivalently,

$$\gamma=\frac1{2\tan(\theta_\gamma/2)}.$$

Under RH, a conjugate pair $\frac12\pm i\gamma$ contributes

$$2[1-\cos(n\theta_\gamma)]$$

to the Li zero sum. If $N(T)$ is the usual Riemann–von Mangoldt zero-counting function, counting upper-half-plane zeros with multiplicity, define

$$\theta(T)=2\arctan\frac1{2T},\qquad g_n(T)=2[1-\cos(n\theta(T))].$$

Then the proposed derivation is

$$\lambda_n=\int_0^\infty g_n(T)\,dN(T).$$

Stieltjes integration by parts gives

$$\lambda_n=-\int_0^\infty N(T)g_n'(T)\,dT,$$

since the boundary term vanishes using

$$N(T)=O(T\log T),\qquad g_n(T)=O_n(T^{-2}).$$

Since

$$\theta'(T)=-\frac4{4T^2+1},$$

this becomes

$$\lambda_n=8n\int_0^\infty\frac{N(T)\sin(n\theta(T))}{4T^2+1}\,dT.$$

Finally, with

$$T=\frac1{2\tan(\theta/2)}=\frac12\cot\frac\theta2,$$

we have

$$\frac{dT}{4T^2+1}=-\frac14\,d\theta,$$

and therefore, under RH,

$$\boxed{\lambda_n=2n\int_0^\pi N\left(\frac1{2\tan(\theta/2)}\right)\sin(n\theta)\,d\theta}\tag{*}$$

for every $n\ge1$.

The integral appears to converge: as $\theta\to0$,

$$\frac1{2\tan(\theta/2)}\sim\frac1\theta,$$

so

$$N\left(\frac1{2\tan(\theta/2)}\right)=O\left(\frac{\log(1/\theta)}{\theta}\right),$$

while $\sin(n\theta)\sim n\theta$. Thus the integrand is

$$O_n(\log(1/\theta)),$$

which is integrable at $0$.

The proposed converse is the following. Define

$$F(\theta)=N\left(\frac1{2\tan(\theta/2)}\right),\qquad 0<\theta<\pi.$$

Because $N(T)$ is nonnegative and nondecreasing while $\theta\mapsto1/(2\tan(\theta/2))$ is decreasing, $F$ is nonnegative and nonincreasing.

The elementary lemma being used is:

If $F:(0,\pi)\to[0,\infty)$ is nonincreasing and the improper integral converges, then

$$\int_0^\pi F(\theta)\sin(n\theta)\,d\theta\ge0$$

for every $n\ge1$.

Indeed, putting $h=\pi/n$ and translating the $n$ intervals of length $h$ back to $(0,h)$ gives

$$\int_0^\pi F(\theta)\sin(n\theta)\,d\theta=\int_0^h\left[\sum_{j=0}^{n-1}(-1)^jF(x+jh)\right]\sin(nx)\,dx.$$

For fixed $x\in(0,h)$,

$$F(x)\ge F(x+h)\ge F(x+2h)\ge\cdots\ge0,$$

so the alternating sum is nonnegative, and $\sin(nx)\ge0$ on $(0,h)$.

Consequently, if $(*)$ holds for the actual Li coefficients for every $n\ge1$ without assuming RH, then

$$\lambda_n\ge0$$

for every $n$, and Li’s criterion implies RH. Thus

$$\boxed{\mathrm{RH}\iff\left[\lambda_n=2n\int_0^\pi N\left(\frac1{2\tan(\theta/2)}\right)\sin(n\theta)\,d\theta\quad\text{for every }n\ge1\right].}$$

An Elegant Collinearity Problem Involving the Incircle

A problem from a vietnamese Facebook group.

Problem
. Let $ABC$ be a triangle, and let its incircle, with center $I$, be tangent to $BC$ at $D$. Let $M$ and $N$ be the midpoints of $AD$ and $BC$, respectively. Prove that $M$, $I$, and $N$ are collinear. 


Proof 1. Let $BC=a$, $CA=b$, $AB=c$, let $s$ be the semiperimeter, $\Delta=[ABC]$, and let $r$ be the inradius. Assume $b>c$; the case $b<c$ is symmetric, while $b=c$ is immediate. Since $BD=s-b$ and $N$ is the midpoint of $BC$,
$$DN=\frac a2-(s-b)=\frac{b-c}{2}.$$
As $ID\perp BC$ and $ID=r$,
$$\tan\angle DNI=\frac{r}{DN}=\frac{2r}{b-c}.$$
Let $M'$ be the foot of the perpendicular from $M$ to $BC$. Since $M$ is the midpoint of $AD$,
$$[CDM]=\frac12[CDA]=\frac14\,b(s-c)\sin C.$$
On the other hand,
$$[CDM]=\frac12(s-c)MM',$$
hence
$$MM'=\frac{b\sin C}{2}.$$
By the cosine law in $\triangle ACD$,
$$AD^2=b^2+(s-c)^2-2b(s-c)\cos C.$$
Since $MD=AD/2$,
$$DM'^2=DM^2-MM'^2=\frac14\bigl(b\cos C-(s-c)\bigr)^2.$$
Using
$$\sin^2\frac C2=\frac{(s-a)(s-b)}{ab},\qquad\cos^2\frac C2=\frac{s(s-c)}{ab},$$
we get
$$\cos C=\cos^2\frac C2-\sin^2\frac C2=\frac{s(s-c)-(s-a)(s-b)}{ab},$$
so
$$b\cos C-(s-c)=\frac{(s-a)(b-c)}{a}.$$
Therefore
$$DM'=\frac{(s-a)(b-c)}{2a}.$$
Since $M',D,N$ occur in this order,
$$NM'=ND+DM'=\frac{b-c}{2}+\frac{(s-a)(b-c)}{2a}=\frac{s(b-c)}{2a}.$$
Thus
$$\tan\angle DNM=\frac{MM'}{NM'}=\frac{ab\sin C}{s(b-c)}=\frac{2\Delta}{s(b-c)}=\frac{2r}{b-c}.$$
Hence
$$\tan\angle DNM=\tan\angle DNI,$$
and therefore $M,I,N$ are collinear.

Proof 2This problem is essentially a degenerate case of Newton’s theorem for tangential quadrilaterals: in every tangential quadrilateral, the center of the inscribed circle lies on the line joining the midpoints of the two diagonals. As one of the vertices of the quadrilateral approaches the point of tangency $D$ on $BC$, the quadrilateral degenerates into triangle $ABC$, and the Newton line becomes precisely the line containing the midpoint $M$ of $AD$, the incenter $I$, and the midpoint $N$ of $BC$.

lunes, 17 de agosto de 2026

Problema 3 de la OMCC-2026

Sean $ABC$ un triángulo y $D$ el punto de tangencia de su incírculo con el segmento $BC$. Sean $P$ un punto del segmento $AB$ y $Q$ un punto del segmento $AC$ tales que $AP = AQ$. Sea $\Gamma_1$ la circunferencia tangente a los segmentos $AB$ y $BC$ en los puntos $P$ y $R$, respectivamente. Análogamente, sea $\Gamma_2$ la circunferencia tangente a los segmentos $AC$ y $BC$ en los puntos $Q$ y $S$, respectivamente. Los circuncírculos de $APR$ y $AQS$ se intersecan en un punto $T$ distinto de $A$. Demuestra que $A$, $D$ y $T$ están en una misma recta.


Demostración. Usando notación estándar,

$$BC=a,\qquad AC=b,\qquad AB=c,$$ 
$$\angle{BAC}=\alpha,\qquad \angle{ABC}=\beta,\qquad \angle{BCA}=\gamma$$
 y 
$$s=\frac{a+b+c}{2}.$$ 

Denotemos por $\omega_1$ y $\omega_2$ a los circuncírculos de $APR$ y $AQS$, respectivamente. Ademas, sea $AP=AQ=x$.

Asumimos que $BR<BD$ y $BS>BD$.

Lema 1. Se cumple que $DR=DS$.

Demostración. Es bien sabido, por propiedades de tangentes comunes a una circunferencia, que $CD=s-c$ y $BD=s-b$, por lo tanto, 

$$DS=CD-CS=(s-c)-(b-x)=s+x-(b+c)$$
y
$$DR=BD-BR=(s-b)-(c-x)=s+x-(b+c),$$
consecuentemente, $DS=DR$.
$\square$

Lema 2. El cuadrilátero $PRSQ$ es cíclico. 

Demostración. Al ser $BP$ y $BR$ tangentes comunes a $\Gamma_1$, esto implica que $BP=BR$, por lo que una sencilla busqueda de ángulo da $\angle{PRS}=\frac{180^{\circ}+\beta}{2}$. Por otro lado, como $APQ$ es isósceles con $AP=AQ$ y $CS=CQ$, se deduce que 

$$\begin{aligned}\angle{PQS}&=180^{\circ}-\angle{AQP}-\angle{CQS}\\&=180^{\circ}-\frac{180^{\circ}-\alpha}{2}-\frac{180^{\circ}-\gamma}{2}\\&=\frac{\alpha+\gamma}{2}\end{aligned}.$$

Así, $\angle{PRS}+\angle{PQS}=180^{\circ}.$

$\square$

De vuelta al problema original

Como $AP=AQ$ y $DR=DS$ (por el Lema 1), se sigue que $AD$ es el eje radical de $\Gamma_1$ y $\Gamma_2$. Llamemos $\rho$ a la circunferencia circunscrita de $PRSQ$ (Lema 2). La recta $PR$ es el eje radical $\Gamma_1$ y $\rho$, y de $\Gamma_1$ y $\omega_1$. Similarmente, la recta $QS$ es el eje radical de $\Gamma_2$ y $\rho$, y de $\omega_2$ y $\rho$. Supongamos que $U=PR \cap QS$. Por el teorema de ejes radicales para $\Gamma_1$, $\Gamma_2$ y $\rho$, la recta $AD$ pasa por $U$, de modo que $AU$ contiene a $D$. Considera ahora las circunferencias $\omega_1$, $\omega_2$ y $\rho$. $\omega_1$ y $\rho$ tienen como eje radical $PR$ y $\omega_2$ y $\rho$ tienen a $QS$ como eje radical, intersecándose en $U$. El eje radical entre $\omega_1$ y $\omega_2$ es $AT$ y debe contener a $U$ por el teorema de los ejes radicales. Por tanto, $AU$ contiene a $D$ y $T$ y por consiguiente, $A$, $D$ y $T$ estan alineados. 

$\square$

Nota. Para los casos $BR=AD=BS$ o $BR>BD$ y $BS<BD$, la demostracion solo requiere ajustes menores.

martes, 11 de agosto de 2026

A Problem by Ivan Pavlov

This problem was originally posted on Facebook.

Problem (Ivan Pavlov). Let \(ABC\) be a triangle with incenter \(I\). Let \(A_1,B_1,C_1\) be the reflections of \(I\) across the sides \(BC,CA,AB\), respectively, and let \(A_2,B_2,C_2\) be the reflections of \(A,B,C\) about \(I\), respectively. Prove that the lines \(A_1A_2\),\(B_1B_2\),\(C_1C_2\) are concurrent.

Proof. Let \(D,E,F\) be the intouch triangle of \(ABC\), and let \(I\) be its incenter. Since \(A_1,B_1,C_1\) are the reflections of \(I\) across \(BC,CA,AB\), respectively, we have

\[A_1B_1C_1=h_{I,2}(DEF),\]

while, since \(A_2,B_2,C_2\) are the reflections of \(A,B,C\) about \(I\), we have

\[A_2B_2C_2=h_{I,-1}(ABC).\]

Now apply the same homothety \(h_{I,\frac12}\) to both triangles. Then \(A_1B_1C_1\) is sent to \(DEF\), whereas \(A_2B_2C_2\) is sent to

\[h_{I,\frac12}\circ h_{I,-1}(ABC)=h_{I,-\frac12}(ABC).\]

By Theorem 3 of B.~Suceavă and P.~Yiu, The Feuerbach Point and Euler Lines, Forum Geometricorum 6 (2006), 191--197, every homothetic image \(h_{I,t}(ABC)\) is perspective with the intouch triangle \(DEF\); hence \(DEF\) and \(h_{I,-\frac12}(ABC)\) are perspective. Therefore their corresponding joins are concurrent. Since \(h_{I,\frac12}\) is an affine transformation, and affine transformations preserve incidence, the inverse images of those three concurrent joins, namely \(A_1A_2\),\(B_1B_2\),\(C_1C_2\) are concurrent as well.                                                   

jueves, 23 de julio de 2026

Lobachevsky’s Angle of Parallelism Inside the USM

When developing the Unified Substitution Method (USM) for integration (see arXiv:2505.03754), the primary goal was methodological: to create a single, branch-consistent calculus that unifies circular and hyperbolic substitutions. By defining everything around exponentials of principal inverse trigonometric functions, the framework maps standard radicals and half-angle integrands into rational forms while keeping the necessary domain and sign choices explicit.

However, a closer look at the core identities of the USM, specifically Theorems 1 and 2, reveals something deeper than an algebraic generalization of Euler’s substitutions. On the positive real hyperbolic components, the inverse-trigonometric angles ($\beta$ and $\psi$) appearing in these theorems specialize exactly to Lobachevsky’s angle of parallelism. The full principal-branch identities extend the same algebraic parametrization, with controlled complex branches, beyond the positive real domain where the ordinary geometric interpretation applies.

Here is a look at the geometry operating under the hood of the USM.


The Complex Engine of the USM

The foundation of the USM relies on evaluating $e^{\pm i\cos^{-1}(y)}$ and $e^{\pm i\sec^{-1}(y)}$ on their principal complex branches. This allows the framework to cover the relevant real components within one principal-branch calculus, with explicit sign and endpoint conventions.

When establishing the transforms, we defined specific angles. In Theorem 1, for the domain $\vert y\vert \ge 1$, we set:

$$\beta = \csc^{-1}(y)$$

In Theorem 2, for the domain $\vert y\vert \le 1$, we set:

$$\psi = \sin^{-1}(y)$$

In both theorems, the central algebraic identities involve the tangents of these half-angles:

$$e^{\pm i\alpha} = \tan\left(\frac{\beta}{2}\right)$$

and

$$e^{\pm i\phi} = \tan\left(\frac{\psi}{2}\right)$$

with the signs determined by the relevant real component and the principal-branch convention.

Algebraically, these identities rationalize the corresponding integrals. But geometrically, what exactly are $\beta$ and $\psi$?


The Gudermannian Connection

Restrict attention to the positive components where these variables acquire their direct hyperbolic-geometric interpretation.

Let $u \ge 0$. The Gudermannian function may be defined by:

$$\operatorname{gd}(u) = \arctan(\sinh u)$$

Because:

$$\tan(\operatorname{gd}(u)) = \sinh u$$

we obtain:

$$\cos(\operatorname{gd}(u)) = \frac{1}{\sqrt{1+\sinh^2 u}} = \frac{1}{\cosh u} = \operatorname{sech} u$$

Thus:

$$\boxed{\cos(\operatorname{gd}(u)) = \operatorname{sech} u}$$

Since:

$$0 \le \operatorname{gd}(u) < \frac{\pi}{2}$$

its complementary angle lies in the interval $(0, \pi/2]$ and satisfies:

$$\sin\left(\frac{\pi}{2}-\operatorname{gd}(u)\right) = \operatorname{sech} u$$

Therefore:

$$\boxed{\arcsin(\operatorname{sech} u) = \frac{\pi}{2}-\operatorname{gd}(u)}$$

This identity identifies the angles appearing in both USM theorems.


Theorem 2

For the positive component of Theorem 2, set:

$$y = \operatorname{sech} u, \qquad 0 < y \le 1$$

Because:

$$\psi = \sin^{-1}(y)$$

we have:

$$\psi = \arcsin(\operatorname{sech} u)$$

Consequently:

$$\boxed{\psi = \frac{\pi}{2}-\operatorname{gd}(u)}$$

Thus $\psi$ is precisely the complementary Gudermannian angle.


Theorem 1

For the positive component of Theorem 1, set:

$$y = \cosh u, \qquad y \ge 1$$

Because:

$$\beta = \csc^{-1}(y)$$

we obtain:

$$\begin{aligned} \beta &= \operatorname{arccsc}(\cosh u) \\ &= \arcsin\left(\frac{1}{\cosh u}\right) \\ &= \arcsin(\operatorname{sech} u) \end{aligned}$$

Therefore:

$$\boxed{\beta = \frac{\pi}{2}-\operatorname{gd}(u)}$$

The two USM angles are consequently identical:

$$\boxed{\beta (\operatorname{cosh} u) = \psi (\operatorname{sech} u) = \frac{\pi}{2}-\operatorname{gd}(u)}$$

This conclusion is obtained before invoking the geometry of the angle of parallelism.


Lobachevsky’s Angle of Parallelism

Lobachevsky’s angle of parallelism, denoted by $\Pi(u)$, is classically complementary to the Gudermannian:

$$\Pi(u) + \operatorname{gd}(u) = \frac{\pi}{2}$$

Equivalently:

$$\Pi(u) = \frac{\pi}{2}-\operatorname{gd}(u)$$

Combining this with the identities obtained above gives:

$$\boxed{\beta (\operatorname{cosh} u) = \psi (\operatorname{sech} u) = \Pi(u)}$$

Thus the angles appearing in the positive hyperbolic components of Theorems 1 and 2 are literally Lobachevsky’s angle of parallelism.

The complete chain is:

$$\boxed{\Pi(u) = \operatorname{arccsc}(\cosh u) = \arcsin(\operatorname{sech} u) = \frac{\pi}{2}-\operatorname{gd}(u)}$$

This also reproduces the defining trigonometric relationship:

$$\sin\Pi(u) = \operatorname{sech} u$$


The Half-Angle Parameters


The geometric interpretation becomes especially clear when we examine the half-angle parameters used by the USM.

Lobachevsky’s classical formula is:

$$\boxed{\tan\left(\frac{\Pi(u)}{2}\right) = e^{-u}}$$

In Theorem 1, after setting:

$$y = \cosh u$$

the principal inverse cosine satisfies:

$$\alpha = \cos^{-1}(\cosh u) = iu$$

Therefore, on the positive component:

$$e^{i\alpha} = e^{i(iu)} = e^{-u}$$

Theorem 1 also gives:

$$e^{i\alpha} = \tan\left(\frac{\beta}{2}\right)$$

Since $\beta = \Pi(u)$:

$$\boxed{e^{i\alpha} = \tan\left(\frac{\beta}{2}\right) = \tan\left(\frac{\Pi(u)}{2}\right) = e^{-u}}$$

Similarly, in Theorem 2, after setting:

$$y = \operatorname{sech} u$$

we have:

$$\phi = \sec^{-1}(\operatorname{sech} u) = \cos^{-1}(\cosh u) = iu$$

Hence:

$$e^{i\phi} = e^{i(iu)} = e^{-u}$$

Theorem 2 gives:

$$e^{i\phi} = \tan\left(\frac{\psi}{2}\right)$$

Since $\psi = \Pi(u)$:

$$\boxed{e^{i\phi} = \tan\left(\frac{\psi}{2}\right) = \tan\left(\frac{\Pi(u)}{2}\right) = e^{-u}}$$

Thus the two apparently different USM parameters coincide on their positive hyperbolic components:

$$\boxed{t (\operatorname{cosh} u)= r (\operatorname{sech} u) = e^{-u}}$$

More completely:

$$\boxed{\begin{aligned} t(\cosh u) &= e^{i\arccos(\cosh u)} = \tan\left(\frac{\operatorname{arccsc}(\cosh u)}{2}\right) = e^{-u}, \\[2mm] r(\operatorname{sech} u) &= e^{i\operatorname{arcsec}(\operatorname{sech} u)} = \tan\left(\frac{\arcsin(\operatorname{sech} u)}{2}\right) = e^{-u}. \end{aligned}}$$

The same imaginary principal inverse angle $iu$ appears through two reciprocal hyperbolic coordinates:

$$\cosh u \qquad \text{and} \qquad \operatorname{sech} u$$

Theorem 1 reaches it through:

$$\cos^{-1}(\cosh u) = iu$$

while Theorem 2 reaches it through:

$$\sec^{-1}(\operatorname{sech} u) = \cos^{-1}(\cosh u) = iu$$


What This Means for the USM


The parameters driving the USM transforms are not arbitrary algebraic devices. On the positive real hyperbolic components, the substitution parameter is exactly:

$$e^{-u}$$

which is equivalently the tangent of half the angle of parallelism:

$$e^{-u} = \tan\left(\frac{\Pi(u)}{2}\right)$$

The two central USM angles are therefore different inverse-function representations of the same geometric quantity:

$$\boxed{\beta = \operatorname{arccsc}(\cosh u) = \Pi(u)}$$

and

$$\boxed{\psi = \arcsin(\operatorname{sech} u) = \Pi(u)}$$

Their common Gudermannian representation is:

$$\boxed{\beta (\operatorname{cosh} u) = \psi (\operatorname{sech} u) = \Pi(u) = \frac{\pi}{2}-\operatorname{gd}(u)}$$

The full identities of Theorems 1 and 2 also include the negative real components, where principal-branch signs must be handled separately. Those cases belong to the analytic and algebraic extension of the USM rather than to the direct interpretation as the ordinary positive angle of parallelism.

On the positive hyperbolic components, however, the geometric meaning is exact: the USM parameters calculate the exponential of the negative hyperbolic distance, the tangent of half Lobachevsky’s angle of parallelism, and the half-angle parameter of the complementary Gudermannian angle—all as the same quantity.

Note. Here $u$ is the dimensionless hyperbolic distance in the standard curvature-$-1$ normalization. If the hyperbolic plane has curvature
$$K=-\frac{1}{R^2}$$
and $d$ denotes geometric distance, then $u=d/R$, and Lobachevsky’s formula becomes
$$e^{-d/R}.$$

martes, 14 de julio de 2026

Radical Integrals Made Easier: Five Families That Collapse into Laurent Polynomials

Introduction

The standard objective of a substitution involving a quadratic radical is rationalization. We replace the radical by a rational function of a new parameter and then hope that the resulting rational integral is manageable.

For several important classes of integrals, however, the Unified Substitution Method (USM) does something substantially stronger.

It does not merely produce a rational function. After the radical and the Jacobian interact, the transformed integrand becomes a finite Laurent polynomial:

\[L(p)=\sum_{k=m}^{N}c_kp^k,\]

where \(p\) is one of the USM parameters \(t\), \(r\), or \(s\).

Consequently,

\[\int L(p)\,dp = \sum_{\substack{k=m\\k\neq-1}}^{N} \frac{c_k}{k+1}p^{k+1} +c_{-1}\ln|p|+C.\]

There is no partial-fraction decomposition, no Hermite reduction, and no recursive trigonometric integration. The calculation is reduced to integrating powers of a single parameter.

This “Laurent collapse” can be viewed as the systematic principle already visible in Examples 4 and 6 of the public arXiv preprint A Unified Substitution Method for Integration (see arXiv:2505.03754v3). In those examples, a quotient radical collapses to a Laurent polynomial in \(t\), while a circular-radical integral collapses to one in \(r\). 


The parametrizations behind the collapse

Let

\[X=x+b,\qquad a>0.\]

The five families considered below arise primarily from the difference-radical and circular-radical transformations. A closely related Laurent mechanism also appears under Transform 3 for the sum radical.


The difference-radical parameter

On either exterior component \(X>a\) or \(X<-a\), let

\[R=\sqrt{X^2-a^2},\qquad \varepsilon=\operatorname{sgn}(X),\]

and set

\[t=\frac{X-\varepsilon R}{a}.\]

Then

\[X=\frac{a}{2}\left(t+t^{-1}\right),\]

\[R=-\varepsilon\frac{a}{2}\left(t-t^{-1}\right),\]

and

\[\,dx=\frac{a(t^2-1)}{2t^2}\,dt.\]

Thus \(0<t<1\) on the component \(X>a\), whereas \(-1<t<0\) on the component \(X<-a\). The formulas remain Laurent in \(t\), with the component-dependent sign carried by \(\varepsilon\).


The circular-radical parameter

For \(|X|<a\), set

\[r=\frac{X}{a+\sqrt{a^2-X^2}}.\]

Equivalently,

\[X=\frac{2ar}{1+r^2},\]

\[\sqrt{a^2-X^2} =a\frac{1-r^2}{1+r^2},\]

and

\[\,dx=2a\frac{1-r^2}{(1+r^2)^2}\,dr.\]

If \(X=a\sin\theta\) with the principal choice \(\theta=\sin^{-1}(X/a)\in(-\pi/2,\pi/2)\), then \(r=\tan(\theta/2)\), the classical tangent half-angle parameter.

The arXiv preprint explains how the corresponding USM transformations recover Euler's first and second substitutions, up to explicit sign, scaling, or reciprocal reparametrizations. Thus, the comparison below is not between unrelated methods: USM reorganizes the classical substitutions in a way that makes the cancellations easier to detect.


Transform 3: Laurent polynomials for the sum radical

The Laurent-polynomial phenomenon is not limited to the five families associated with Transforms 2 and 5. Transform 3 produces the same kind of algebraic simplification for integrals containing

\[\sqrt{X^2+a^2}.\]

Define

\[s=\frac{X+\sqrt{X^2+a^2}}{a}>0.\]

Then

\[X=\frac a2\left(s-s^{-1}\right),\]

\[\sqrt{X^2+a^2} =\frac a2\left(s+s^{-1}\right),\]

and

\[\,dx=\frac a2\left(1+s^{-2}\right)\,ds.\]

Each of the three fundamental quantities \(X\), \(\sqrt{X^2+a^2}\), and \(\,dx\) is therefore already a finite Laurent polynomial in \(s\).

It follows immediately that an expression of the form

\[Q\left(X,\sqrt{X^2+a^2}\right)\,dx,\]

where \(Q\) is a polynomial in its two arguments, becomes a finite Laurent polynomial in \(s\).

The same remains true for many integrands containing integer powers of

\[X+\sqrt{X^2+a^2}=as,\]

because multiplication or division by such a power merely introduces a Laurent monomial in \(s\).

For example, Example 3 of the arXiv preprint considers

\[\int_0^\infty \frac{\,dx}{\left(x+\sqrt{1+x^2}\right)^2}.\]

With

\[s=x+\sqrt{1+x^2},\]

we have

\[\,dx=\frac12(1+s^{-2})\,ds\]

and

\[\left(x+\sqrt{1+x^2}\right)^{-2}=s^{-2}.\]

Thus,

\[\int_0^\infty \frac{\,dx}{\left(x+\sqrt{1+x^2}\right)^2} = \frac12\int_1^\infty\left(s^{-2}+s^{-4}\right)\,ds = \frac23.\]

This is another Laurent collapse: the transformed integrand contains only two negative powers of \(s\).


Traditional comparison

The classical hyperbolic substitution

\[x=\sinh u\]

gives

\[x+\sqrt{1+x^2} =\sinh u+\cosh u =e^u.\]

The integral then becomes

\[\int_0^\infty e^{-2u}\cosh u\,du.\]

This is already manageable; a standard next step is to expand \(\cosh u\) into exponentials. Setting \(s=e^u\) then produces the same Laurent integral obtained directly by Transform 3.

Transform 3 therefore packages two classical steps,

\[x=\sinh u, \qquad s=e^u,\]

into the single algebraic parameter

\[s=x+\sqrt{1+x^2}.\]

Its advantage is not that hyperbolic substitution fails, but that the final Laurent structure is exposed from the beginning.


The five Laurent-polynomial families

The cancellation mechanism applies systematically to the following classes. Each statement is understood on a connected real component on which its integrand and substitution are defined.

For every polynomial \(P\), on either component \(X>a\) or \(X<-a\),

\[\int P(x)\sqrt{X^2-a^2}\,dx\]

becomes a finite Laurent polynomial in \(t\).

Indeed,

\[P(x) = P\left(\frac{a}{2}(t+t^{-1})-b\right)\]

is itself a Laurent polynomial, while

\[\sqrt{X^2-a^2}\,dx = -\varepsilon\frac{a^2}{4}t^{-3}(t^2-1)^2\,dt,\]

where

\[\varepsilon=\operatorname{sgn}(X)\]

is constant on the selected component.


Polynomial times a quotient radical

For every polynomial \(P\), on either component \(X>a\) or \(X<-a\),

\[\int P(x)\sqrt{\frac{X-a}{X+a}}\,dx\]

also becomes a finite Laurent polynomial in \(t\).

The decisive cancellation is

\[\sqrt{\frac{X-a}{X+a}}\,dx = -\frac a2(1-t^{-1})^2\,dt.\]

The linear factor \(1+t\) contributed by the denominator of the radical cancels a matching factor in the Jacobian.


Reciprocal powers divided by a circular radical

For every integer \(n\geq1\), on a connected component of \(0<|X|<a\),

\[\int\frac{\,dx}{X^n\sqrt{a^2-X^2}}\]

becomes

\[2^{1-n}a^{-n} \int r^{-n}(1+r^2)^{n-1}\,dr.\]

Because \(n-1\geq0\), the factor \((1+r^2)^{n-1}\) expands finitely.


A circular radical divided by reciprocal powers

For every integer \(n\geq3\), on a connected component of \(0<|X|<a\),

\[\int\frac{\sqrt{a^2-X^2}}{X^n}\,dx\]

becomes

\[2^{1-n}a^{2-n} \int r^{-n}(1-r^2)^2(1+r^2)^{n-3}\,dr.\]

The restriction \(n\geq3\) is exactly what prevents a residual denominator involving \(1+r^2\).


Reciprocal powers multiplied by a root ratio

For every integer \(n\geq2\), on a connected component of \(0<|X|<a\) and with principal square roots,

\[\int \frac1{X^n} \sqrt{\frac{a\pm X}{a\mp X}}\,dx\]

becomes

\[2^{1-n}a^{1-n} \int r^{-n}(1\pm r)^2(1+r^2)^{n-2}\,dr.\]

Once again, the exponent \(n-2\) is nonnegative, so the transformed expression is a finite Laurent polynomial.

These formulas show that the restrictions on \(n\) are structural rather than arbitrary: they are precisely the thresholds at which all remaining powers of \(1+r^2\) move into the numerator.


Five examples: USM versus the traditional routes

Example 1: A polynomial times a difference radical

Consider

\[I_1=\int x^2\sqrt{x^2-1}\,dx, \qquad x>1.\]

Set

\[t=x-\sqrt{x^2-1}.\]

Then

\[x=\frac12(t+t^{-1}),\]

and direct substitution gives

\[I_1 = \int\left( -\frac{t^3}{16} +\frac1{8t} -\frac1{16t^5} \right)\,dt.\]

The transformed integrand contains only three Laurent monomials. Therefore,

\[I_1 = -\frac{t^4}{64} +\frac18\ln t +\frac{t^{-4}}{64} +C.\]

Using

\[t^{-1}=x+\sqrt{x^2-1},\]

we obtain

\[\boxed{ I_1= \frac18\left[ x\sqrt{x^2-1}(2x^2-1) -\ln\left(x+\sqrt{x^2-1}\right) \right]+C }.\]


Traditional comparison

The hyperbolic substitution \(x=\cosh u\) produces

\[\int \cosh^2u\,\sinh^2u\,du,\]

which is commonly handled with multiple-angle identities before back-substitution.

The trigonometric substitution \(x=\sec\theta\) produces an integral involving \(\sec^3\theta\tan^2\theta\), normally handled through identities and reduction formulas.

Both methods are valid. The USM advantage is that the complete algebraic structure is exposed immediately: three powers of \(t\), integrated term by term.


Example 2: The quotient-radical cancellation

Consider the integral appearing as Example 4 in the arXiv preprint:

\[I_2=\int\sqrt{\frac{x+1}{x+3}}\,dx, \qquad x>-1.\]

Here

\[X=x+2, \qquad a=1,\]

so the USM parameter is

\[t=x+2-\sqrt{x^2+4x+3}.\]

The radical and Jacobian combine to give

\[I_2 = -\frac12\int\left(1-2t^{-1}+t^{-2}\right)\,dt.\]

Thus,

\[I_2 = -\frac t2+\ln t+\frac1{2t}+C.\]

Since

\[\frac1{2t}-\frac t2 =\sqrt{x^2+4x+3},\]

the result is

\[\boxed{ I_2= \sqrt{x^2+4x+3} +\ln\left(x+2-\sqrt{x^2+4x+3}\right)+C }.\]

The important feature is not merely that the integral is solvable. It is that the transformed integrand is already a Laurent polynomial.


Traditional comparison

The natural direct rationalization is

\[u=\sqrt{\frac{x+1}{x+3}}.\]

Solving for \(x\) and differentiating gives

\[I_2=\int\frac{4u^2}{(u^2-1)^2}\,du.\]

The radical has disappeared, but the price is a repeated quadratic denominator. A standard continuation uses partial fractions or introduces a second hyperbolic substitution.

The USM parameter packages both steps into one transformation:

\[-\frac12\left(1-2t^{-1}+t^{-2}\right).\]

Moreover, multiplying the original integrand by any polynomial \(P(x)\) preserves the Laurent form, because

\[x=\frac12(t+t^{-1})-2.\]

The direct \(u\)-substitution, by contrast, generally increases the degree of the rational denominator as the polynomial factor becomes more complicated.

This is the first seed of the general Laurent-polynomial phenomenon: Example 4 demonstrates the cancellation in one particular integral, while the corresponding family shows that it persists for every polynomial multiplier.


Example 3: Reciprocal powers and a circular radical

Now consider Example 6 from the arXiv preprint:

\[I_3=\int\frac{\,dx}{x^3\sqrt{4-x^2}}, \qquad 0<x<2.\]

Set

\[r=\frac{2-\sqrt{4-x^2}}{x} =\frac{x}{2+\sqrt{4-x^2}}.\]

Transform 5 gives

\[I_3 = \frac1{32}\int\left(r^{-3}+2r^{-1}+r\right)\,dr.\]

Therefore,

\[I_3 = -\frac1{64r^2} +\frac1{16}\ln r +\frac{r^2}{64} +C.\]

After combining the reciprocal powers,

\[\boxed{ I_3= \frac1{16}\left[ \ln\left(\frac{2-\sqrt{4-x^2}}{x}\right) -\frac{2\sqrt{4-x^2}}{x^2} \right]+C }.\]


Traditional comparison

The usual circular substitution

\[x=2\sin\theta\]

produces

\[I_3=\frac18\int\csc^3\theta\,d\theta.\]

That integral is standard and is commonly evaluated by a reduction formula or an integration-by-parts derivation, followed by the reconstruction of both \(\csc\theta\) and \(\cot\theta\) in terms of \(x\).

The USM route bypasses the reduction formula entirely. The radical and Jacobian cancel before integration, leaving exactly three Laurent monomials.

This is the second worked example in the arXiv preprint from which the general Laurent classification naturally emerges.


Example 4: A circular radical over a high reciprocal power

Consider

\[I_4=\int\frac{\sqrt{1-x^2}}{x^5}\,dx, \qquad 0<x<1.\]

Use

\[r=\frac{x}{1+\sqrt{1-x^2}}.\]

This is the fourth family with \(a=1\) and \(n=5\). The general formula gives

\[I_4 = \frac1{16}\int r^{-5}(1-r^2)^2(1+r^2)^2\,dr.\]

But

\[(1-r^2)^2(1+r^2)^2=(1-r^4)^2,\]

so

\[I_4 = \frac1{16}\int\left(r^{-5}-2r^{-1}+r^3\right)\,dr.\]

Hence

\[I_4 = \frac{r^4-r^{-4}}{64} -\frac18\ln r+C.\]

Returning to \(x\),

\[\boxed{ I_4= \frac18\left[ \ln\left(\frac{1+\sqrt{1-x^2}}{x}\right) -\frac{(2-x^2)\sqrt{1-x^2}}{x^4} \right]+C }.\]


Traditional comparison

With \(x=\sin\theta\), one obtains

\[I_4 = \int\frac{\cos^2\theta}{\sin^5\theta}\,d\theta = \int\left(\csc^5\theta-\csc^3\theta\right)\,d\theta.\]

A standard evaluation uses reduction formulas for two odd powers of \(\csc\theta\). It is entirely feasible, but the method treats the integral as a trigonometric-reduction problem.

USM instead detects the algebraic factorization

\[r^{-5}-2r^{-1}+r^3\]

automatically. The same pattern persists for every \(n\geq3\); increasing \(n\) merely enlarges a finite binomial expansion.


Example 5: Reciprocal powers times a root ratio

Consider

\[I_5= \int\frac1{x^3}\sqrt{\frac{1+x}{1-x}}\,dx, \qquad 0<x<1.\]

Again set

\[r=\frac{x}{1+\sqrt{1-x^2}}.\]

Because

\[\sqrt{\frac{1+x}{1-x}} =\frac{1+r}{1-r},\]

the fifth Laurent family gives

\[I_5 = \frac14\int r^{-3}(1+r)^2(1+r^2)\,dr.\]

Expanding,

\[I_5 = \frac14\int\left(r^{-3}+2r^{-2}+2r^{-1}+2+r\right)\,dr.\]

Thus,

\[I_5 = -\frac1{8r^2} -\frac1{2r} +\frac12\ln r +\frac r2 +\frac{r^2}{8} +C.\]

Combining reciprocal pairs gives

\[\boxed{ I_5= -\frac{\sqrt{1-x^2}}{2x^2} -\frac{\sqrt{1-x^2}}{x} +\frac12\ln\left(\frac{x}{1+\sqrt{1-x^2}}\right)+C }.\]

The conjugate family

\[\frac1{x^n}\sqrt{\frac{1-x}{1+x}}\]

is handled identically, with \((1+r)^2\) replaced by \((1-r)^2\).


Traditional comparison

The substitution \(x=\sin\theta\) gives

\[\sqrt{\frac{1+\sin\theta}{1-\sin\theta}} =\frac{1+\sin\theta}{\cos\theta}\]

on the relevant component. Consequently,

\[I_5 = \int\left(\csc^3\theta+\csc^2\theta\right)\,d\theta.\]

A standard evaluation combines the odd-cosecant reduction formula with a separate elementary integration of \(\csc^2\theta\).

USM reaches the same trigonometric half-angle parameter without first converting the problem into trigonometric notation. More importantly, it reveals immediately that the entire class for \(n\geq2\) is Laurent-polynomial after substitution.

Compressing \(t^n\pm t^{-n}\): the binomial difference and sum formulas

Let \(n\in\mathbb N\) and \(y\in(-\infty,-1]\cup[1,\infty)\). Termwise integration often produces reciprocal pairs such as

\[t^n-t^{-n} \qquad\text{and}\qquad t^n+t^{-n}.\]

Let

\[w=\sqrt{y^2-1}.\]

For \(y\geq1\), one has \(t=y-w\) and \(t^{-1}=y+w\); for \(y\leq-1\), these two reciprocal roots exchange roles. Treating the two roots as an unordered pair avoids unnecessary expression swell.

Indeed,

\[(y-w)(y+w)=1,\]

so

\[\{y-w,y+w\}=\{t,t^{-1}\}.\]

Define

\[D_n(y)=(y-w)^n-(y+w)^n.\]

By the binomial theorem, all even powers of \(w\) cancel, leaving only the odd terms:

\[D_n(y) = -2\sum_{j=0}^{\lfloor(n-1)/2\rfloor} \binom{n}{2j+1} y^{n-2j-1}w^{2j+1}.\]

Since

\[w^{2j+1}=(y^2-1)^jw,\]

this may also be written as

\[D_n(y) = -2\sqrt{y^2-1} \sum_{j=0}^{\lfloor(n-1)/2\rfloor} \binom{n}{2j+1} y^{n-2j-1}(y^2-1)^j.\]

For the USM parameter chosen as the small root on both exterior components,

\[t= \begin{cases} y-\sqrt{y^2-1},&y\geq1,\\[4pt] y+\sqrt{y^2-1},&y\leq-1, \end{cases}\]

the result is

\[t^n-t^{-n}=\sigma D_n(y),\]

where

\[\sigma= \begin{cases} 1,&y\geq1,\\ -1,&y\leq-1. \end{cases}\]

The factor \(\sigma\) is essential: the roles of \(y-w\) and \(y+w\) are interchanged on the negative exterior component.

The first few identities are

\[D_1(y)=-2\sqrt{y^2-1},\]

\[D_2(y)=-4y\sqrt{y^2-1},\]

\[D_3(y)=-2(4y^2-1)\sqrt{y^2-1},\]

and

\[D_4(y)=-8y(2y^2-1)\sqrt{y^2-1}.\]

There is an equally useful companion formula for the sum. Define

\[S_n(y)=(y-w)^n+(y+w)^n.\]

The odd powers now cancel, giving

\[S_n(y) = 2\sum_{j=0}^{\lfloor n/2\rfloor} \binom{n}{2j} y^{n-2j}w^{2j},\]

or

\[S_n(y) = 2\sum_{j=0}^{\lfloor n/2\rfloor} \binom{n}{2j} y^{n-2j}(y^2-1)^j.\]

Unlike the difference, the sum is unchanged when the two reciprocal roots are interchanged. Therefore,

\[t^n+t^{-n}=S_n(y)\]

on both exterior components, with no additional sign factor. Equivalently, in terms of the Chebyshev polynomials \(T_n\) and \(U_{n-1}\),

\[S_n(y)=2T_n(y),\qquad D_n(y)=-2\sqrt{y^2-1}\,U_{n-1}(y).\]

For example,

\[t^2+t^{-2}=4y^2-2,\]

and

\[t^4+t^{-4}=16y^4-16y^2+2.\]

These formulas shorten the final back-substitution substantially. In Example 1, the Laurent primitive contains

\[\frac{t^{-4}-t^4}{64}.\]

For \(x>1\), the difference formula gives

\[t^4-t^{-4} =-8x(2x^2-1)\sqrt{x^2-1}.\]

Hence

\[\frac{t^{-4}-t^4}{64} = \frac18x(2x^2-1)\sqrt{x^2-1},\]

which recovers the algebraic part of the antiderivative in one step.

The arXiv preprint explicitly presents the binomial-difference identity as a device for simplifying terms of the form \(t^n-t^{-n}\) during back-substitution (see arXiv:2505.03754v3).


Why Laurent polynomials are better than general rational functions

The difference is computationally meaningful.

A general rational function

\[\frac{A(p)}{B(p)}\]

may require factorization of \(B\), polynomial division, repeated-factor decomposition, irreducible quadratic terms, or Hermite reduction.

A Laurent polynomial

\[\sum c_kp^k\]

requires none of these operations. Its integration complexity is visible at inspection:

• every exponent except \(-1\) produces another power;

• the exponent \(-1\) produces a logarithm;

• no new algebraic denominator is introduced;

• the number of terms is controlled by finite polynomial or binomial expansions.

There is also an advantage during back-substitution. Positive and negative powers often occur in pairs such as

\[t^m-t^{-m} \qquad\text{or}\qquad t^m+t^{-m}.\]

The binomial difference and sum formulas convert these pairs directly into expressions involving \(y\) and \(\sqrt{y^2-1}\). This avoids separately expanding two large reciprocal powers and then simplifying the result.

That is why the final answers in the examples above can be written using only the original radical and one logarithm, despite passing through several Laurent powers.


What USM is—and is not—claiming

The Laurent-polynomial property is a precise structural statement. For the stated families and domain restrictions, the transformed integrand has a finite Laurent expansion.

It does \emph{not} imply that every individual integral will have a shorter final antiderivative than every classical derivation. For small exponents, a familiar trigonometric substitution may be equally short. Nor does it constitute a universal theorem about computer-algebra performance.

The real advantage is systematic predictability:

1. the appropriate parameter is selected from the radical geometry;

2. the branch and real component are stated explicitly;

3. cancellation occurs before integration;

4. the resulting integral is termwise;

5. reciprocal powers can be recombined before back-substitution;

6. the procedure scales uniformly with \(P(x)\) or \(n\).

Traditional substitutions usually teach these examples as separate categories: Euler substitutions for difference radicals, circular substitutions for \(\sqrt{a^2-x^2}\), hyperbolic substitutions for \(\sqrt{x^2+a^2}\), half-angle identities for root ratios, and reduction formulas for powers of secant or cosecant.

USM shows that they are manifestations of a common algebraic mechanism.


Final perspective

Examples 4 and 6 of the arXiv preprint display two striking simplifications:

\[\sqrt{\frac{x+1}{x+3}}\,dx \quad\longrightarrow\quad -\frac12\left(1-2t^{-1}+t^{-2}\right)\,dt,\]

and

\[\frac{\,dx}{x^3\sqrt{4-x^2}} \quad\longrightarrow\quad \frac1{32}\left(r^{-3}+2r^{-1}+r\right)\,dr.\]

Transform 3 supplies a parallel sum-radical example:

\[\frac{1}{\left(x+\sqrt{1+x^2}\right)^2}\,dx \quad\longrightarrow\quad \frac12\left(s^{-2}+s^{-4}\right)\,ds.\]

The five systematic families explain the deeper pattern behind the first two examples, while Transform 3 shows that the Laurent principle extends naturally beyond those families.

The best substitution does not merely remove the radical. It makes the radical, the powers of the centered variable, and the Jacobian simplify one another before the integration begins.

Once the integration has been performed, the binomial difference and sum formulas complete the process by compressing reciprocal pairs such as

\[t^n\pm t^{-n}\]

directly into the original variable and radical.

For these classes, USM moves the problem beyond rationalization. What initially appears to be a radical integral becomes finite arithmetic with powers—and, at worst, one logarithm.


Reference

Emmanuel Antonio José García, A Unified Substitution Method for Integration, arXiv:2505.03754v3 [math.GM], revised May 23, 2026.