8
$\begingroup$

Consider the following extension of polynomials. The rational exponential expressions (REXes) are given by:

  1. The leaves 1 and $x$ for $x$ drawn from a class of variables; and
  2. Closed under the binary functions of addition, multiplication, exponentiation, and division.

This is a restriction to positive constants of a class of expressions studied in Brown, 1969, "Rational Exponential Expressions and a Conjecture Concerning π and e". Buchberger & Loos, 1982, "Algebraic Simplification", sect. 6, mention the existence of algorithms for finding canonical forms for these expressions, but say their correctness depends on number-theoretic conjectures that have not been settled.

I'm interested in a decision procedure for dominance, where for two REXes in one variable, $f$ and $g$, $f \leq g$ if $\exists N \in {\mathbb N}.\forall n \in {\mathbb N}. N \leq n \implies f(n) \leq g(n)$. Do we have either the existence of canonical forms, or decidability for this problem?

$\endgroup$
6
  • $\begingroup$ Could you clarify the precise sense of "monomial" that you mean here? And you only have one variable for the decision problem? And you only require dominance for integer n in that question? (And you may want to add the logic tag.) $\endgroup$ Jan 26, 2010 at 19:13
  • $\begingroup$ It looks like it would make sense either with unary exponentiation exp(x) and with binary exponentiation y^x. Which do you mean? $\endgroup$ Jan 26, 2010 at 20:04
  • $\begingroup$ @Charles: Heh. The 'poly' in polynomial does not refer to the multiplicity of variables but of terms... You'd only improve your question by replacing 'monomial' by a more usual description. Also, please do not make corrections and clarifications to the question in the comments but in the actual question: you can edit it! $\endgroup$ Jan 26, 2010 at 22:02
  • $\begingroup$ Joel: logic tag? My interest in the problem does stem from ordinal analysis, but the link is too obscure to justify it on those grounds. Why else might it count as logic? $\endgroup$ Jan 27, 2010 at 10:08
  • $\begingroup$ It seems logic related, since (i) you are considering terms in a formal language, and (ii) you are asking a decidability question. Also, the paper to which you link seems to engage with several logic issues, including the MRDP theorem. $\endgroup$ Jan 27, 2010 at 15:02

3 Answers 3

12
$\begingroup$

Yes, there are algorithms to decide asymptotic dominance - in fact for a much wider class of elementary functions. I discovered the first such algorithm circa 1980 while an undergrad member of the MIT Mathlab group researching effective algorithms for computing limits for the Macsyma symbolic computation system. Another different algorithm was discovered independently a handful of years later by John Shackell. You should be able to find references to the literature by googling the more recent buzzword "transseries".

Many computer algebra systems have (partial) implementations of these algorithms. Here's an example of my algorithm on (my generalization of) an example Rich Schroeppel proposed to attempt to stump my algorithm (he was convinced no such algorithm existed). It shows that $\rm{\: lim_{x\to\infty}\ d40 = e^a}$ alt text
Don't dare try L'Hopital's rule on that monster!

For further discussion (and the text form of the above image) see my post on sci.math, 1996/03/20, L'Hopital's rule question:
http://groups.google.com/group/sci.math/msg/05298104ac44efd2
http://groups.google.com/groups?selm=WGD.96Mar20231913%40berne.ai.mit.edu

$\endgroup$
1
  • $\begingroup$ I've finally got around to get an idea of what is going on with transseries. So, as I understand it, by tuning the construction of a class of transmonomials (typically an inverted order), you can represent certain spaces of function as infinite power series over these. I haven't seen how transseries, which seem to be uncountable, are represented in computer algebra systems, and how, given heterogenous transmonomials P < P' and Q < Q', we can decide whether PQ' < P'Q. $\endgroup$ Feb 8, 2012 at 12:43
4
$\begingroup$

The integers are a red herring. All of the rational exponential expressions are logarithmico-exponential functions in the sense of Hardy. Therefore, for any such $f$ and $g$, either $f(t) > g(t)$ for all sufficiently large real $t$, $f(t) < g(t)$ for all sufficiently large real $t$ or $f(t)=g(t)$. There is no reason to restrict to the case that $t$ is an integer.

I highly suspect that this problem is decidable, but I don't actually know this. This is the sort of thing that people who study transseries think about.

$\endgroup$
4
  • $\begingroup$ I copied these links from the related question: mathoverflow.net/questions/3057/… $\endgroup$ Jan 27, 2010 at 17:17
  • $\begingroup$ In the case of multi-variable polynomials, the difference between the reals and the integers is the difference between decidability and undecidability. This is becasue Tarski's theorem shows that the decision problem for whether p(x1,...,xn) = 0 has a solution is decidable in the reals, but MRDP shows it is undecidable in the integers. But I agree that the question is likely decidable in the one-variable case. $\endgroup$ Jan 27, 2010 at 18:17
  • $\begingroup$ This "eventual trichotomy" property is something I had been thinking about: thank you for the reference (and cf. archive.org/details/ordersofinfinity00harduoft); I must clearly study the proof. I'll take a look at the transseries intro paper; this is a new field for me. $\endgroup$ Jan 28, 2010 at 11:02
  • 2
    $\begingroup$ The crucial difficulty is that all proofs of the sign preservation at infinity for Hardy functions I know have the following logical step in some form: if the derivative preserves sign near infinity, then the function itself preserves sign near infinity. This is completely useless from the computational standpoint because it doesn't allow you to figure out which of the 2 possible signs the function preserves in finite number of steps. The transseries approach clearly has a similar problem. You have to decide if infinitely many terms are equal before you can proceed to the next order terms. $\endgroup$
    – fedja
    Jan 28, 2010 at 12:53
1
$\begingroup$

I don't have an answer to the one-variable decision question that you asked.

But you did introduce multi-variable expressions, and there are some fascinating results on the multi-variable analogue of your domination problem. (Perhaps you know all this already...) That is, given two rational exponential expressions $f(x_1,\ldots,x_k)$ and $g(x_1,\ldots,x_k)$, the problem is to decide whether $f(n_1,\ldots,n_k)\leq g(n_1,\ldots,n_k)$, for all positive integers $(n_1,\ldots,n_k)$ except at most finitely often.

The answer is that this decision problem is undecidable, and it is undecidable even for polynomial expressions. There is no algorithm that will tell you when one multi-variable (positive) polynomial expression eventually dominates another.

The reason is that the decision problem of Diophantine equations is encodable into this domination problem. That work famously shows that the problem to decide if a given integer polynomial $p(x_1,\ldots,x_k)$ has a zero in the integers is undecidable. This is the famous MRDP solution to Hilbert's 10th problem.

It is easy to reduce the Diophantine problem to the domination problem, as follows. First, let us restrict to non-negative integers, for which the MRDP results still apply. Suppose we are given a polynomial expression $p(x_1,\ldots,x_k)$ over the integers, and want to decide if it has a solution in the natural numbers. This expression may involve some minus signs, which your expressions do not allow, but we will take care of that by moving all the minus signs to one side. Introduce a new variable $x_0$ and consider the domination problem:

  • Does $1\leq (1+n_0)p(n_1,\ldots,n_k)^2$ for all natural numbers except finitely often?

We can expand the right hand side, and move the negative signs to the left, to arrive at an instance of your domination problem, using only positive polynomials. Now, if $p(n_1,\ldots,n_k)$ is never $0$, then the answer to the stated domination problem is Yes, since the right hand side will always be at least $1$ in this case. Conversely, if $p(n_1,\ldots,n_k)=0$ has a solution, then we arrive at infinitely many violations of domination by using any choice of $n_0$. Thus, if we could decide the domination problem, then we could decide whether $p(n_1,\ldots,n_k)=0$ has solutions in the natural numbers, which we cannot do by the MRDP theorem. QED

I think the best situation with the Diophantine equations is that it remains undecidable with nine variables, and so the domination problem I described above is undecidable with ten variables. (Perhaps Bjorn Poonen will show up here and tell us a better answer?)

Of course, this doesn't answer the one-variable question that you asked, and probably you know all this already.

My final remark is that if one can somehow represent the inverse pairing function, then one will get undecidability even in the one variable case. That is, let $f(n,m) = \frac{(n+m)(n+m+1)}{2} + m$ be one of the usual pairing functions, which is bijective between $\omega^2$ and $\omega$. Let p be the function such that $p(f(n,m)) = n$, the projection of the pairs onto the first coordinate. If the expressions are enriched to allow $p$, then one can in effect work with several variables by coding them all via pairing into one variable, and in this case, the domination problem in the one-variable case, for rational exponential expressions also allowing the function $p$, will be undecidable. It would seem speculative to suppose that $p$ is itself equivalent to a rational exponential expression, but do you know this?

$\endgroup$
5
  • $\begingroup$ I think that, if we add inverse pairing functions to Diophantines, we can only code up undecidable problems if we have two variables: we can't existentially quantify over the variable that ranges over the problem set. $\endgroup$ Jan 28, 2010 at 10:40
  • $\begingroup$ No, one variable will suffice, once we have the inverse pairing function p. First, it is easy to solve the equation to express also the other projection function q. Now, for example, one can think of triples (a,b,c) as coded by iterated pairs n=((a,b),c), and then extract the coordinates from the single variable n using a=p(p(n)), b=q(p(n)) and c=q(n). Thus, any expression in integer variables (a,b,c) can be thought of as an expression just in one variable n. $\endgroup$ Jan 28, 2010 at 13:08
  • $\begingroup$ This doesn't deal with my worry. Consider the non-primeness Diophantine set characterised by (a+2)*(b+2)-n=0. The variable n plays a different role in determining the set of non-prime numbers: the a and b are existentially quantified, while the n is the argument to the predicate. I don't see how variables of different roles in this sense can be usefully simplified with pairing/projection functions. $\endgroup$ Jan 29, 2010 at 19:08
  • $\begingroup$ What I am claiming is that for every multivariable integer domination problem f(a,b,c) <= g(a,b,c), there is an equivalent single variable problem F(n) <= G(n), using rational expressions involving the inverse pairing function, such that the original inequality holds except finitely often iff the new inequality holds except finitely often. And the proof of this is to think of n as coding the triple (a,b,c) via the polynomial coding function, and then replace all occurences of a with p(p(n)), all occurences of b with q(p(n)) and c with q(n). Your non-primeness set is not a domination problem. $\endgroup$ Jan 29, 2010 at 20:05
  • $\begingroup$ Got it! Thanks. Given the space of possible pairing/projection functions, the issue isn't easy to dismiss. $\endgroup$ Jan 29, 2010 at 21:55

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.