Constructing nonresidues in finite fields and the extended Riemann hypothesis

We present a new deterministic algorithm for the problem of constructing th power nonresidues in finite fields , where is prime and is a prime divisor of . We prove under the assumption of the Extended Riemann Hypothesis (ERH), that for fixed and , our algorithm runs in polynomial time. Unlike other deterministic algorithms for this problem, this polynomial-time bound holds even if is exponentially large. More generally, assuming the ERH, in time we can construct a set of elements that generates the multiplicative group . An extended abstract of this paper appeared in Proc. 23rd Ann. ACM Symp. on Theory of Computing, 1991.

     

Computing the canonical height on K3 surfaces

Let be a surface in given by the intersection of a (1,1)-form and a (2,2)-form. Then is a K3 surface with two noncommuting involutions and . In 1991 the second author constructed two height functions and which behave canonically with respect to and , and in 1993 together with the first author showed in general how to decompose such canonical heights into a sum of local heights . We discuss how the geometry of the surface is related to formulas for the local heights, and we give practical algorithms for computing the involutions , , the local heights , , and the canonical heights , .

     

Expansion and Estimation of the Range of Nonlinear Functions

Many verification algorithms use an expansion , for , where the set of matrices is usually computed as a gradient or by means of slopes. In the following, an expansion scheme is described which frequently yields sharper inclusions for . This allows also to compute sharper inclusions for the range of over a domain. Roughly speaking, has to be given by means of a computer program. The process of expanding can then be fully automatized. The function need not be differentiable. For locally convex or concave functions special improvements are described. Moreover, in contrast to other methods, may be empty without implying large overestimations for . This may be advantageous in practical applications.

     

An algorithm for matrix extension and wavelet construction

This paper gives a practical method of extending an matrix , , with Laurent polynomial entries in one complex variable , to a square matrix also with Laurent polynomial entries. If has orthonormal columns when is restricted to the torus , it can be extended to a paraunitary matrix. If has rank for each , it can be extended to a matrix with nonvanishing determinant on . The method is easily implemented in the computer. It is applied to the construction of compactly supported wavelets and prewavelets from multiresolutions generated by several univariate scaling functions with an arbitrary dilation parameter.

     

On beta expansions for Pisot numbers

Given a number , the beta-transformation is defined for by (mod 1). The number is said to be a beta-number if the orbit is finite, hence eventually periodic. In this case is the root of a monic polynomial with integer coefficients called the characteristic polynomial of . If is the minimal polynomial of , then for some polynomial . It is the factor which concerns us here in case is a Pisot number. It is known that all Pisot numbers are beta-numbers, and it has often been asked whether must be cyclotomic in this case, particularly if . We answer this question in the negative by an examination of the regular Pisot numbers associated with the smallest 8 limit points of the Pisot numbers, by an exhaustive enumeration of the irregular Pisot numbers in (an infinite set), by a search up to degree in , to degree in , and to degree in . We find the smallest counterexample, the counterexample of smallest degree, examples where is nonreciprocal, and examples where is reciprocal but noncyclotomic. We produce infinite sequences of these two types which converge to from above, and infinite sequences of with nonreciprocal which converge to from below and to the th smallest limit point of the Pisot numbers from both sides. We conjecture that these are the only limit points of such numbers in . The Pisot numbers for which is cyclotomic are related to an interesting closed set of numbers introduced by Flatto, Lagarias and Poonen in connection with the zeta function of . Our examples show that the set of Pisot numbers is not a subset of .

     

Complete solutions of a family of quartic Thue and index form equations

Continuing the recent work of the second author, we prove that the diophantine equation

for has exactly 12 solutions except when , when it has 16 solutions. If denotes one of the zeros of , then for we also find all with .

     

An asymptotic expansion for the incomplete beta function

A new asymptotic expansion is derived for the incomplete beta function , which is suitable for large , small and . This expansion is of the form

where is the incomplete Gamma function ratio and . This form has some advantages over previous asymptotic expansions in this region in which depends on as well as on and .

     

The serial test for a nonlinear pseudorandom number generator

Let and A sequence is obtained by the formula The sequence is a sequence of pseudorandom numbers of the maximal period length if and only if (mod 4), (mod 4). In this note, the uniformity is investigated by the 2-dimensional serial test for the sequence. We follow closely the method of papers by Eichenauer-Herrmann and Niederreiter.

     

Turán's Pure Power Sum Problem

Let be complex numbers, and consider the power sums , . Put , where the minimum is over all possible complex numbers satisfying the above. Turán conjectured that , for some positive absolute constant. Atkinson proved this conjecture by showing . It is now known that , for . Determining whether or approaches some other limiting value as is still an open problem. Our calculations show that an upper bound for decreases for , suggesting that decreases to a limiting value less than as .

     

On integral bases in relative quadratic extensions

Let be an algebraic number field and a quadratic extension with . We describe a minimal set of elements for generating the integral elements of as an module. A consequence of this theoretical result is an algorithm for constructing such a set. The construction yields a simple procedure for computing an integral basis of as well. In the last section, we present examples of relative integral bases which were computed with the new algorithm and also give some running times.

     

Double roots of power series and related matters

For a given collection of distinct arguments , multiplicities and a real interval containing zero, we are interested in determining the smallest for which there is a power series with coefficients in , and roots of order respectively. We denote this by . We describe the usual form of the extremal series (we give a sufficient condition which is also necessary when the extremal series possesses at least non-dependent coefficients strictly inside , where is 1 or 2 as is real or complex). We focus particularly on , the size of the smallest double root of a power series lying on a given ray (of interest in connection with the complex analogue of work of Boris Solomyak on the distribution of the random series ). We computed the value of for the rationals in of denominator less than fifty. The smallest value we encountered was . For the one-sided intervals and the corresponding smallest values were and .

     

Minus class groups of the fields of the -th roots of unity

We show that for any prime number the minus class group of the field of the -th roots of unity admits a finite free resolution of length 1 as a module over the ring . Here denotes complex conjugation in . Moreover, for the primes we show that the minus class group is cyclic as a module over this ring. For these primes we also determine the structure of the minus class group.

     

Two new points in the spectrum of the absolute Mahler measure of totally positive algebraic integers

For totally positive algebraic integers of degree , we consider the set of values of , where is the Mahler measure of . C. J. Smyth has found the four smallest values of and conjectured that the fifth point is . We prove that this is so and, moreover, we give the sixth point of .

     

Spontaneous Generation of Modular Invariants

It is possible to compute and its modular equations with no perception of its related classical group structure except at . We start by taking, for prime, an unknown ``-Newtonian' polynomial equation with arbitrary coefficients (based only on Newton's polygon requirements at for and ). We then ask which choice of coefficients of leads to some consistent Laurent series solution , (where . It is conjectured that if the same Laurent series works for -Newtonian polynomials of two or more primes , then there is only a bounded number of choices for the Laurent series (to within an additive constant). These choices are essentially from the set of ``replicable functions,' which include more classical modular invariants, particularly . A demonstration for orders and is done by computation. More remarkably, if the same series works for the -Newtonian polygons of 15 special ``Fricke-Monster' values of , then is (essentially) determined uniquely. Computationally, this process stands alone, and, in a sense, modular invariants arise ``spontaneously.'

     

Optimal information for approximating periodic analytic functions

Let be a strip in the complex plane. For fixed integer let denote the class of -periodic functions , which are analytic in and satisfy in . Denote by the subset of functions from that are real-valued on the real axis. Given a function , we try to recover at a fixed point by an algorithm on the basis of the information

where , are the Fourier coefficients of . We find the intrinsic error of recovery

Furthermore the -dimensional optimal information error, optimal sampling error and -widths of in , the space of continuous functions on , are determined. The optimal sampling error turns out to be strictly greater than the optimal information error. Finally the same problems are investigated for the class , consisting of all -periodic functions, which are analytic in with -integrable boundary values. In the case sampling fails to yield optimal information as well in odd as in even dimensions.

     

Pure product polynomials and the Prouhet-Tarry-Escott problem

An -factor pure product is a polynomial which can be expressed in the form for some natural numbers . We define the norm of a polynomial to be the sum of the absolute values of the coefficients. It is known that every -factor pure product has norm at least . We describe three algorithms for determining the least norm an -factor pure product can have. We report results of our computations using one of these algorithms which include the result that every -factor pure product has norm strictly greater than if is , , , or .

     

-complete sequences of integers

An infinite sequence is -complete if every sufficiently large integer is the sum of such that no one divides the other. We investigate -completeness of sets of the form and with nonnegative.

     

The BKK root count in

The root count developed by Bernshtein, Kushnirenko and Khovanskii only counts the number of isolated zeros of a polynomial system in the algebraic torus . In this paper, we modify this bound slightly so that it counts the number of isolated zeros in . Our bound is, apparently, significantly sharper than the recent root counts found by Rojas and in many cases easier to compute. As a consequence of our result, the Huber-Sturmfels homotopy for finding all the isolated zeros of a polynomial system in can be slightly modified to obtain all the isolated zeros in .

     

On Quadrature Convergence of Extended Lagrange Interpolation

Quadrature convergence of the extended Lagrange interpolant for any continuous function is studied, where the interpolation nodes are the zeros of an orthogonal polynomial of degree and the zeros of the corresponding ``induced' orthogonal polynomial of degree . It is found that, unlike convergence in the mean, quadrature convergence does hold for all four Chebyshev weight functions. This is shown by establishing the positivity of the underlying quadrature rule, whose weights are obtained explicitly. Necessary and sufficient conditions for positivity are also obtained in cases where the nodes and interlace, and the conditions are checked numerically for the Jacobi weight function with parameters and . It is conjectured, in this case, that quadrature convergence holds for .