On strong tractability of weighted multivariate integration

We prove that for every dimension and every number of points, there exists a point-set whose -weighted unanchored discrepancy is bounded from above by independently of provided that the sequence has for some (even arbitrarily large) . Here is a positive number that could be chosen arbitrarily close to zero and depends on but not on or . This result yields strong tractability of the corresponding integration problems including approximation of weighted integrals over unbounded domains such as . It also supplements the results that provide an upper bound of the form when .

     

Counting primes in residue classes

We explain how the Meissel-Lehmer-Lagarias-Miller-Odlyzko method for computing can be used for computing efficiently , the number of primes congruent to modulo up to . As an application, we computed the number of prime numbers of the form less than for several values of up to and found a new region where is less than near .

     

Rational double points on supersingular surfaces

We investigate configurations of rational double points with the total Milnor number on supersingular surfaces. The complete list of possible configurations is given. As an application, we also give the complete list of extremal (quasi-)elliptic fibrations on supersingular surfaces.

     

The Brumer-Stark conjecture in some families of extensions of specified degree

As a starting point, an important link is established between Brumer's conjecture and the Brumer-Stark conjecture which allows one to translate recent progress on the former into new results on the latter. For example, if is an abelian extension of relative degree , an odd prime, we prove the -part of the Brumer-Stark conjecture for all odd primes with belonging to a wide class of base fields. In the same setting, we study the -part and -part of Brumer-Stark with no special restriction on and are left with only two well-defined specific classes of extensions that elude proof. Extensive computations were carried out within these two classes and a complete numerical proof of the Brumer-Stark conjecture was obtained in all cases.

     

Reducing the construction cost of the component-by-component construction of good lattice rules

The construction of randomly shifted rank- lattice rules, where the number of points is a prime number, has recently been developed by Sloan, Kuo and Joe for integration of functions in weighted Sobolev spaces and was extended by Kuo and Joe and by Dick to composite numbers. To construct -dimensional rules, the shifts were generated randomly and the generating vectors were constructed component-by-component at a cost of operations. Here we consider the situation where is the product of two distinct prime numbers and . We still generate the shifts randomly but we modify the algorithm so that the cost of constructing the, now two, generating vectors component-by-component is only operations. This reduction in cost allows, in practice, construction of rules with millions of points. The rules constructed again achieve a worst-case strong tractability error bound, with a rate of convergence for 0$">.

     

Class numbers of some abelian extensions of rational function fields

Let be a monic irreducible polynomial. In this paper we generalize the determinant formula for of Bae and Kang and the formula for of Jung and Ahn to any subfields of the cyclotomic function field By using these formulas, we calculate the class numbers of all subfields of when and are small.

     

Elliptic curves with nonsplit mod representations

We calculate explicitly the -invariants of the elliptic curves corresponding to rational points on the modular curve by giving an expression defined over of the -function in terms of the function field generators and of the elliptic curve . As a result we exhibit infinitely many elliptic curves over with nonsplit mod representations.

     

An estimate for the number of integers without large prime factors

denotes the number of positive integers and free of prime factors y$">. Hildebrand and Tenenbaum provided a good approximation of . However, their method requires the solution to the equation , and therefore it needs a large amount of time for the numerical solution of the above equation for large . Hildebrand also showed approximates for , where and is the unique solution to . Let be defined by for 0$">. We show approximates , and also approximates , where . Using these approximations, we give a simple method which approximates within a factor in the range , where is any positive constant.

     

Real zeros of real odd Dirichlet -functions

Let be a real odd Dirichlet character of modulus , and let be the associated Dirichlet -function. As a consequence of the work of Low and Purdy, it is known that if and , , , then has no positive real zeros. By a simple extension of their ideas and the advantage of thirty years of advances in computational power, we are able to prove that if , then has no positive real zeros.

     

The convergence of harmonic Ritz values, harmonic Ritz vectors, and refined harmonic Ritz vectors

This paper concerns a harmonic projection method for computing an approximation to an eigenpair of a large matrix . Given a target point and a subspace that contains an approximation to , the harmonic projection method returns an approximation to . Three convergence results are established as the deviation of from approaches zero. First, the harmonic Ritz value converges to if a certain Rayleigh quotient matrix is uniformly nonsingular. Second, the harmonic Ritz vector converges to if the Rayleigh quotient matrix is uniformly nonsingular and remains well separated from the other harmonic Ritz values. Third, better error bounds for the convergence of are derived when converges. However, we show that the harmonic projection method can fail to find the desired eigenvalue --in other words, the method can miss if it is very close to . To this end, we propose to compute the Rayleigh quotient of with respect to and take it as a new approximate eigenvalue. is shown to converge to once tends to , no matter how is close to . Finally, we show that if the Rayleigh quotient matrix is uniformly nonsingular, then the refined harmonic Ritz vector, or more generally the refined eigenvector approximation introduced by the author, converges. We construct examples to illustrate our theory.

     

Some new kinds of pseudoprimes

We define some new kinds of pseudoprimes to several bases, which generalize strong pseudoprimes. We call them Sylow -pseudoprimes and elementary Abelian -pseudoprimes. It turns out that every which is a strong pseudoprime to bases 2, 3 and 5, is not a Sylow -pseudoprime to two of these bases for an appropriate prime

We also give examples of strong pseudoprimes to many bases which are not Sylow -pseudoprimes to two bases only, where or

     

-estimate for the discrete Plateau Problem

In this paper we prove the convergence rates for a fully discrete finite element procedure for approximating minimal, possibly unstable, surfaces.

Originally this problem was studied by G. Dziuk and J. Hutchinson. First they provided convergence rates in the and norms for the boundary integral method. Subsequently they obtained the convergence estimates using a fully discrete finite element method. We use the latter framework for our investigation.

     

Bounds for computing the tame kernel

The tame kernel of the of a number field  is the kernel of some explicit map , where the product runs over all finite primes  of  and is the residue class field at . When is a set of primes of , containing the infinite ones, we can consider the -unit group  of . Then has a natural image in . The tame kernel is contained in this image if  contains all finite primes of  up to some bound. This is a theorem due to Bass and Tate. An explicit bound for imaginary quadratic fields was given by Browkin. In this article we give a bound, valid for any number field, that is smaller than Browkin's bound in the imaginary quadratic case and has better asymptotics. A simplified version of this bound says that we only have to include in  all primes with norm up to  , where  is the discriminant of . Using this bound, one can find explicit generators for the tame kernel, and a ``long enough' search would also yield all relations. Unfortunately, we have no explicit formula to describe what ``long enough' means. However, using theorems from Keune, we can show that the tame kernel is computable.

     

A sensitive algorithm for detecting the inequivalence of Hadamard matrices

A Hadamard matrix of side is an matrix with every entry either or , which satisfies . Two Hadamard matrices are called equivalent if one can be obtained from the other by some sequence of row and column permutations and negations. To identify the equivalence of two Hadamard matrices by a complete search is known to be an NP hard problem when increases. In this paper, a new algorithm for detecting inequivalence of two Hadamard matrices is proposed, which is more sensitive than those known in the literature and which has a close relation with several measures of uniformity. As an application, we apply the new algorithm to verify the inequivalence of the known inequivalent Hadamard matrices of order ; furthermore, we show that there are at least pairwise inequivalent Hadamard matrices of order . The latter is a new discovery.

     

Solutions of the congruence

To supplement existing data, solutions of are tabulated for primes with and . For , five new solutions 2^{32}$"> are presented. One of these, for , also satisfies the ``reverse' congruence . An effective procedure for searching for such ``double solutions' is described and applied to the range , . Previous to this, congruences are generally considered for any and fixed prime to see where the smallest prime solution occurs.

     

Linear algebra algorithms for divisors on an algebraic curve

We use an embedding of the symmetric th power of any algebraic curve of genus into a Grassmannian space to give algorithms for working with divisors on , using only linear algebra in vector spaces of dimension , and matrices of size . When the base field is finite, or if has a rational point over , these give algorithms for working on the Jacobian of that require field operations, arising from the Gaussian elimination. Our point of view is strongly geometric, and our representation of points on the Jacobian is fairly simple to deal with; in particular, none of our algorithms involves arithmetic with polynomials. We note that our algorithms have the same asymptotic complexity for general curves as the more algebraic algorithms in Florian Hess' 1999 Ph.D. thesis, which works with function fields as extensions of . However, for special classes of curves, Hess' algorithms are asymptotically more efficient than ours, generalizing other known efficient algorithms for special classes of curves, such as hyperelliptic curves (Cantor 1987), superelliptic curves (Galbraith, Paulus, and Smart 2002), and curves (Harasawa and Suzuki 2000); in all those cases, one can attain a complexity of .

     

Deciding the nilpotency of the Galois group by computing elements in the centre

We present a new algorithm for computing the centre of the Galois group of a given polynomial along with its action on the set of roots of , without previously computing the group. We show that every element in the centre is representable by a family of polynomials in . For computing such polynomials, we use quadratic Newton-lifting and truncated expressions of the roots of over a -adic number field. As an application we give a method for deciding the nilpotency of the Galois group. If is irreducible with nilpotent Galois group, an algorithm for computing it is proposed.

     

The Steiner triple systems of order 19

Using an orderly algorithm, the Steiner triple systems of order are classified; there are pairwise nonisomorphic such designs. For each design, the order of its automorphism group and the number of Pasch configurations it contains are recorded; of the designs are anti-Pasch. There are three main parts of the classification: constructing an initial set of blocks, the seeds; completing the seeds to triple systems with an algorithm for exact cover; and carrying out isomorph rejection of the final triple systems. Isomorph rejection is based on the graph canonical labeling software nauty supplemented with a vertex invariant based on Pasch configurations. The possibility of using the (strongly regular) block graphs of these designs in the isomorphism tests is utilized. The aforementioned value is in fact a lower bound on the number of pairwise nonisomorphic strongly regular graphs with parameters .

     

Biquadratic reciprocity and a Lucasian primality test

Let be the sequence defined from a given initial value, the seed, , by the recurrence . Then, for a suitable seed , the number (where is odd) is prime iff . In general depends both on and on . We describe a slight modification of this test which determines primality of numbers with a seed which depends only on , provided . In particular, when , odd, we have a test with a single seed depending only on , in contrast with the unmodified test, which, as proved by W. Bosma in Explicit primality criteria for , Math. Comp. 61 (1993), 97-109, needs infinitely many seeds. The proof of validity uses biquadratic reciprocity.