arxiv: v1 [math.co] 22 Jul 2014

Transcription

1 A PROBABILISTIC APPROACH TO VALUE SETS OF POLYNOMIALS OVER FINITE FIELDS arxiv: v1 [mathco] 22 Jul 2014 ZHICHENG GAO AND QIANG WANG Abstract In this paper we study the distribution of the size of the value set for a random polynomial with degree at most q 1 over a finite field F q We obtain the exact probability distribution and show that the number of missing values tends to a normal distribution as q goes to infinity We obtain these results through a study of a random r-th order cyclotomic mappings A variation on the size of the union of some random sets is also considered 1 Introduction Let F q be the finite field of q elements with characteristic p Let γ be a fixed primitive element of F q throughout the paper The value set of a polynomial g over F q is the set V g of images when we view g as a mapping from F q to itself Clearly g is a permutation polynomial PP of F q if and only if the cardinality V g of the value set V g is q Asymptotic formulas such as V g λgq +Oq 1/2, where λg is a constant depending only on certain Galois groups associated to g, can be found in Birch and Swinnerton-Dyer [2] and Cohen [8] Later, Williams [24] proved that almost all polynomials g of degree d satisfy λg 1 1 2! + 1 3! + + 1d 1 1 d! There are also several results on explicit upper bound for V g if g is not a PP over F q ; see for example [14, 19, 20] Perhaps the most well-nown result is due to Wan [20] who proved that if a polynomial g of degree d is not a PP then 1 V g q q 1 d On the other hand, it is easy to see that V g q/d for any polynomial g over F q with degree d The polynomials achieving this lower bound are called minimal value set polynomials The classification of minimal value set polynomials over F p with 2 can be found in 2000 Mathematics Subject Classification 05A16,60E05,11T06 Key words and phrases polynomials, value sets, normal distribution, finite fields Research of authors was partially supported by NSERC of Canada 1

2 2 ZHICHENG GAO AND QIANG WANG [6, 15], and in [3] for all the minimal value set polynomials in F q [x] whose value set is a subfield of F q See [10, 21] for further results on lower bounds of V g and [13] for some classes of polynomials with small value sets More receny, algorithms and complexity in computing V g have been studied in [7] We note that all of these results mentioned above relate V g to the degreedofg It isalsowell nown thatevery polynomial g over F q such that g0 b has the form ax r fx s +b with some positive integers r,s such that s q 1 There are different ways to choose r,s in the form ax r fx s +b However, in [1], the concept of the index of a polynomial was first introduced and any non-constant polynomial g F q [x] of degree q 1 can be written uniquely as gx ax r fx q 1/l +b with index l defined below Namely, write gx ax n +a n i1 x n i 1 + +a n i x n i +b, where a, a n ij 0, j 1,, The case that 0 is trivial Thus, we shall assume that 1 Write n i r, the vanishing order of x at 0 ie, the lowest degree of x in gx b is r Then gx a x r fx q 1/l +b,wherefx x e 0 +a n i1 x e 1 + +a n i 1 x e 1 +a r, l q 1 gcdn r,n r i 1,,n r i 1,q 1 q 1 :, s and gcde 0,e 1,,e 1,l 1 The integer l q 1 s is called the index of gx From the above definition of index l, one can see that the greatest common divisor condition maes l minimal among those possible choices Clearly, the study of the value set of g over F q is equivalent to studying the value set x r fx q 1/l over F q with index l Receny Mullen, Wan and Wang [17] used an index approach to study the upper bound of the value set for any polynomial which is not a PP They proved that if g is not a PP then 2 V g q q 1 l This result improves Wan s result when the index l of a polynomial is stricy smaller than the degree d We note that the index l of a polynomial is always smaller than the degree d as long as l q 1 The above result is obtained through a study of cyclotomic mapping polynomials which were studied earlier in [11, 18, 22] The index of a polynomial is closely related to the concept of the least index of a cyclotomic mapping polynomial Recall that γ is a fixed primitive element of F q Let l q 1 and the set of all nonzero l-th powers be

3 A PROBABILISTIC APPROACH TO VALUE SETS 3 C 0 Then C 0 is a subgroup of F q of index l The elements of the factor group F q/c 0 are the cyclotomic cosets C i : γ i C 0, i 0,1,,l 1 For any a 0,a 1,,a l 1 F q and a positive integer r, the r-th order cyclotomic mapping fa r 0,a 1,,a l 1 of index l from F q to itself see Niederreiter and Winterhof in [18] for r 1 or Wang [22] is defined by { 0, if x 0; 3 fa r 0,a 1,,a l 1 x a i x r, if x C i, 0 i l 1 It is shown that r-th order cyclotomic mappings of index l produce the polynomials of the form x r fx s where s q 1 Indeed, the polynomial presentation is given l by gx 1 l 1 l 1 a i x r ζ ji x js, l i0 where ζ γ s is a fixed primitive l-th root On the other hand, as we mentioned earlier, each polynomial fx such that f0 0 with index l can be written as x r fx q 1/l, which is an r-th order cyclotomic mapping with the least index l such that a i fζ i for i 0,,l 1 In this paper, we are interested in the probability distribution of the value set size of a random r-th order cyclotomic mapping polynomial for any given index l and any positive integer r, as defined in Equation 3 Thus this enables us to derive the probability distribution of the size of value set of a random polynomial of degree d q 1 over a finite field In Section 2 we ouine our method and a crucial result on normal distribution which is used in this paper Essentially, we are interested in the distribution of the size of the union of subsets, namely, the distribution of the random variable X l A j where A j gc j 1 for j 1,,l and t r,s In Section 3, we first consider a simplified model such that none of a i s in Equation 3 is zero Hence a i s are chosen independeny at random from F q This means that the zero is not contained in any one of the subsets A j s In particular, in Theorem 1 we obtain the distribution of the number of missing values for a random r-th order cyclotomic mapping f r a 0,a 1,,a l 1 such that none of a i s is zero Moreover, in Theorem 2, we show that this distribution is asymptotically normal In Section 4, we study any random r-th order cyclotomic mapping polynomial by choosing a i s in Equation 3 independeny at random

4 4 ZHICHENG GAO AND QIANG WANG from F q The probability distribution of value set size is given in Theorem 3 As a consequence, for l q 1, we obtain the exact probability distribution of the value set size of a random polynomial over F q with degree at most q 1 In particular, we have the following corollaries to Theorem 3 Corollary 1 Let gx be a random polynomial of degree at most q 1 over F q with g0 0 Then q 1 q 1 1+j P V g +1 1 j j q Consequeny, for oq, we have P V g +1 1! q 1 +1 q q 1 This proves that if > 1 is small compared to q then the number of polynomials over F q with degree less than or equal to q 1 such that the value set size is is always exponential in q Moreover, we have Corollary 2 Let gx be any random polynomial of degree at most q 1 over finite field F q with g0 0 Let Y q q gf q denote the number of missing nonzero values in the value set of g Let µ q q/e and σ 2 q e 1 2e 2 q Then the distribution of Y q µ q /σ q tends to the standard normal, as q Finally in Section 5 we study a variation of our model used in Section 3 We consider the case when each subset A i is chosen uniformly at random from all m i -subsets of a given n-set for i 1,,l This extends a result by Barot and Peña [4] and David [9] We also show that the size of the complement of l i1 A i is asymptotically normal 2 Methodology In [16], we obtained the following formula for the cardinality of the value set for an arbitrary polynomial Proposition 1 Proposition 23 in [16] Let gx ax r fx s +b a 0 be any polynomial over F q with index l q 1 and let gcdr,s t s Let γ be a fixed primitive element of F q Then V g c s t +1, or V g c 1 s t +1, where c {γ ir fγ si s/t i 0,,l 1}

5 A PROBABILISTIC APPROACH TO VALUE SETS 5 As discussed earlier, it is sufficient to assume that a 1 and b 0 in Proposition 1 That is, we can view gx as a r-th order cyclotomic mapping polynomial with the least index l In this case, we have gx a i x r when x C i, where a i fγ si for i 0,,l 1 Recall that C 0 is the subgroup of F q consisting of all the l-th powers of F q and welet T 0 bethesubgroupoff q consisting ofallthe-thpowers Hence T i with 0 i 1 give all the cyclotomic cosets of index We also note that x r maps C 0 onto T 0 which contains s distinct elements t So x r maps each coset C i γ i C 0 onto γ ir T 0 Therefore g maps C i onto γ ir fγ si T 0, which could be either the set {0} if a i fγ si 0 or one of the nonzero cyclotomic cosets of index We observe that c is the number of distinct cyclotomic cosets of the form γ ir fγ si T 0, possibly along with the subset {0} if one of a i s is zero Hence we have V g c s +1 or t c 1s +1, the latter happens when some of a t i s in gx a i x r equal 0 Therefore the value set problem for a random r-th order cyclotomic mapping polynomial or random polynomial g, essentially requires us to study the number c in Proposition 1, the size of union of some cyclotomic cosets and possibly the subset {0} if a i s tae zero More specifically, for 0 i l 1, each C i is mapped to A i+1 gc i which is one of T 0,,T 1 or {0} Then c is the number of distinct A j s 1 j l and the value set size is either c s +1 or t c 1s +1 More t generally, we are interested in the distribution of the random variable X l A j, whilea j arechosen independeny according toagiven distribution depending on that a i is chosen independeny at random from F q Similar problems have been studied in [4, 9] when each A j is chosen uniformly at random from all -subsets of a given n-set Let n and let D 0 {0} and D j T j 1 for 1 j 1 Let Y n be the the number of D 1,,D n 1 which are not in l A j for a random r-th order cyclotomic mapping polynomial with index l such that r,s t We will derive exact probability distributions of Y n and show that they are asymptotically normal Throughout the paper, we shall use Y n to denote the falling factorial Y n Y n 1Y n 2 Y n +1 WeuseP,E, Vtodenotetheprobability, expectation, and variance of a random variable, respectively The main tool for deriving the probability distribution of Y n in the paper is through the sieve method and the falling factorial moments Let B 1,,B n be n events in a probability space and N {1,,n} We note that PY n is the probability that exacy of the B j

6 6 ZHICHENG GAO AND QIANG WANG occur Define S h J N, J h P j J B j Then the well-nown sieve formula See eg [5, Theorm 10] gives 4 n h PY n 1 h S h, 5 h n h S PY n h Hence h 6 EY n!s We also need the following result [12, Theorem 1] in order to show that Y n is asymptotically normal Lemma 1 Let s n > µ 1 n and σ n µ n +µ 2 ns n, where µ n as n Suppose that µ n oσ 3 n, and a sequence Y n of nonnegative random variables satisfies EY n µ 2 n exp s n, 2 uniformly for all integers in the range cµ n /σ n c µ n /σ n for some constants c > c > 0 Then Y n µ n /σ n tends in distribution to the standard normal as n 3 Sizes of value sets of cyclotomic mapping polynomials with nonzero branches In this section, we study of the value set size of a random r-th order cyclotomic mapping polynomial with nonzero branches ie, none of a i s is zero This means that we choose a i in 3 independeny at random from F q and it leads to the following model Let n and D 1,D 2,D bepairwisedisjointsubsetsoff q suchthat D i s/tfor all 1 i Because a i is chosen independeny at random from F q, this means that A 1,A 2,,A l are chosen independeny and uniformly atrandomfrom{d 1,D 2,,D } Weareinterestedinthedistribution ofx n l A j Thisiscloselyrelatedtothedistributionsofllabeled balls into labeled boxes Let Y n be the number of empty boxes in a

7 A PROBABILISTIC APPROACH TO VALUE SETS 7 random distribution of l labeled balls into labeled boxes We note X n n s/ty n We first prove the following Theorem 1 Let Y n be the number of empty boxes in a random distribution of l labeled balls into n labeled boxes We have l EY n, 1 PY n 1 j j l l j Proof Let N {1,2,,} Let B j be the event that box j is empty We have S l P j J B j P B j J N, J It follows from 6 that EY n Using Equation 4, we obtain n h PY n 1 h Next we obtain h 1 h h h h 1 l S h l h h 1 h h l h 1 j j Theorem 2 Suppose t ol 1/5 as n Define µ n te 1/t l, σ 2 n te 2/t e 1/t 1 1/tl Then the distribution of Y n µ n /σ n tends to the standard normal, as n j l l

8 8 ZHICHENG GAO AND QIANG WANG Proof From Theorem 1, we have, as l, 7 l EY n 1 1 exp lln exp l exp µ n, t 1 2t 2 l +O 1 2t 2 l +O t 3 l 2 VY n EY n Y n 1+EY n EY n l l l 1exp lln 1 2 +exp lln exp 2lln exp l O exp l O exp 2l O exp l exp 2l 3 +O exp l exp l 3 2 +O exp 2l exp l 3 +O 2 8

9 A PROBABILISTIC APPROACH TO VALUE SETS 9 2 1exp l 1 2l exp l 1 l exp 2l 1 l 2 9 e 1 t e 2 t le 2 t σ 2 n Under the assumption that t ol 1/5, we can verify that µ2 n 0 as σn 6 n We note µ n /σ n is of the order t 3 l Hence for in the range specified in Lemma 1, we have 2 Ot 3 l Therefore 1 1 j 1 exp ln 1 j 1 j exp + j2 j 3 2t 2 l +O 2 t 3 l 3 exp 2 2 l exp lln exp l exp, t 2 2t 2 l 3 2t 2 l +O 2 t 3 l 3 It follows from 6 and 7 that EY n t 2 exp 2 2 µ n 2t 2 l exp 2 t+1 2t 2 l Now the result follows from Lemma 1 and the estimation s n σ2 n µ n µ 2 n as t l EY ny n 1 µ 2 n 0 when l t 2 l +O 1 t 2 l 2 > µ 1 n,

10 10 ZHICHENG GAO AND QIANG WANG The following corollary follows immediately from Theorem 1, by noting s t 1 and X n n Y n Corollary 3 Under the assumption of Theorem 1 and X n n Y n, we have PX n h 1 h n h 1 h j j n n n h j In particular, and for h on, PX n n n! n n, PX n h 1 h! nh h n n Consider l n q 1 and t 1 Because an r-th order cyclotomic mapping polynomial f r a 0,,a q 2 x always maps 0 to 0, Corollary 3 implies Corollary 4 Let gx be any random r-th order cyclotomic mapping polynomial fa r 0,,a q 2 x over finite field F q such that none of a i s is zero Then the probability of V f h+1 is given by h 1 q 1 h PX q 1 h 1 h j j q 1 q 1 q 1 h j In particular, the probability of such a random cyclotomic mapping polynomial f r a 0,,a q 2 x is a permutation polynomial is PX q 1 q 1 q 1! q 1 q 1, and for h oq, the probability of such a random cyclotomic mapping polynomial fa r 0,,a q 2 x with a value set size h+1 is PX q 1 h 1 q 1 h h! q 1h q 1 4 Sizes of value sets of random polynomials Recall q 1 ls and r is a positive integer such that r,s t In this section we consider the value set size for any random r-th order cyclotomic mapping polynomial f r a 0,,a l 1 x with index l over finite fieldf q Namely, a 0,,a l 1 areindependeny chosen atrandomfrom F q Wenotethatthevaluesetproblemforanyrandompolynomialwith

11 A PROBABILISTIC APPROACH TO VALUE SETS 11 degree at most q 1 is in fact the value set problem for a random r-th order cyclotomic mapping polynomial with index l q 1 As discussed in Section 2, we are interested in the distribution of X n l A j However, we need to include the element 0 in our analysis similar to those in Section 3 Because each a i in 3 is chosen independeny at random from F q, this leads to the following model Let us define D 0 {0} and D j T j 1 for 1 j The distribution of each random set A i is PA i D 0 1 q, and PA i D j s q tq, j 1,,l As in the case discussed earlier, we define the events B j l i1 {A i D j }, 0 j Let Y be the number of B 1,B 2,,B note B 0 is excluded which occur Then PY is the probability that exacy of cyclotomic sets T j s of index are not in the value set of g Then we have Lemma 2 Let q 1 ls and r be a positive integer such that r,s t Let gx be any random r-th order cyclotomic mapping polynomial f r a 0,,a l 1 x with index l over finite field F q Let Y be the number of cyclotomic sets of index not contained in the value set of g Then EY 1 s tq PY Proof Let N {1,2,,} Define S l, 1 j j J N, J P j J B j 1 Because A 1,,A l are mutually independent, we have S P B j Using Equation 4, we obtain 1 q + s l tq q + sj tq l 1 s l tq

12 12 ZHICHENG GAO AND QIANG WANG PY h 1 h h 1 h h h h S h 1 h q + s h tq 1 h h 1 j j l 1 q + s h tq 1 q + sj l tq Using the above lemma, we can obtain the distribution of X V g Theorem 3 Let q 1 ls and r be a positive integer such that r, s t Let fx be any random r-th order cyclotomic mapping polynomial fa r 0,,a l 1 x with index l over F q Then Proof Indeed, PX 1+s/t 1 j j 1 q + sj l tq PX 1+s/t P{Y } 1 1 j j q + sj l tq We now obtain an application of the above results on random polynomials with degree at most q 1 It is sufficient to study the distribution of the subclass of random polynomials gx such that g0 0 with degree q 1, because any polynomial fx such that f0 b can be written as gx + b with g0 0 and vice versa So we can view any random polynomial g of degree at most q 1 with g0 0 and index l as an r-th order cyclotomic mapping polynomial with the least index l, which is therefore a cyclotomic mapping polynomial with index q 1 Because there are q q 1 such polynomials with g0 0, they correspond to all the cyclotomic mapping polynomials with index q 1 So a random polynomial with degree less than or equal to q 1 and g0 0 is a random r-th order cyclotomic mapping polynomial with index l q 1 for any r 1 Therefore, Lemma 2 implies that any random polynomial with degree q 1 has expected value set size l

13 A PROBABILISTIC APPROACH TO VALUE SETS q q 1 q This verifies William s result [24] saying that almost e all the polynomials of degree q 1 is a general polynomial Moreover, applying Theorem 3 to the case l q 1 hence s t 1, we obtain exact probability distribution of the value set size for a random polynomial over finite field F q in Corollary 1 Moreover, we can drive Corollary 2 that the distribution of Y u /σ tends to the standard normal distribution as well, as l and t ol 1/5, following the same arguments as in the proof of Theorem 2 5 Size of the Union of Random Sets Recall N {1,,n} In Section 3, we considered l i1 A i where each A i a randomcyclotomic cosets of index l In Section 4, we allowed the possibility that A i could be {0}, a subset with a different size from cyclotomic cosets Earlier, Barot and Peña [4] considered the probability distribution of l i1a i where each A i is chosen independeny and uniformly at random from P m where P m is the set of all m-subsets of N In this section, we consider l i1 A i such that m i s can be distinct Let X n l i1 A i and Y n n X n In the following we establish the distributions of X n and Y n and we show that the Y n is asymptotically normal This generalizes the following result by Barot and Peña [4] because we allow m i s to be distinct Theorem 4 Barot and Peña, 2001 Let N {1,,n} Let X n l i1 A i where each A i is chosen independeny and uniformly at random from P m where P m is the set of all m-subsets of N and Y n n X n We have PX n i EX n n 1 n i m l i i i j n l 1 m j, j m l, 1 m n VX n nn 1 1 m l 1 m l EX n 2 +EX n n n 1 First of all, we extend the above results to general m i s Lemma 3 Let N {1,,n} and P mi be the set of all m i -subsets of N for 1 i l Let X n be the random variable for the size of l i1a i where A i is a random set chosen independeny and uniformly

14 14 ZHICHENG GAO AND QIANG WANG from P mi and Y n n X n We have EY n n l PY n PX n i n m j n, n h 1 h h h0 l n h i h m j i l n i 1 h i h n m j n m j h n h, Proof For each j N, let B j denote the event that j / l i1 A i We note that PY n is the probability that exacy of the B j s occur Since A 1,A 2,,A l are mutually independent, we have n S PB 1 B 2 B and hence l P1 / A j,2 / A j,, / A j n l n n m j n, EY n!s n l n m j n Nowfrom4weobtainPY n n h 1h h n l h Finally, PX n i PY n n i n l h n 1 h n+i n i h hn i h0 i h m j i l n i 1 h i h n, m j n m j h n h n m j h n h where the last equality holds after the substitution h : h n+i Next we show that Y n is asymptotically normal The condition on u i in the following theorem can be relaxed considerably, but we use the current form for the sae of simplicity

15 A PROBABILISTIC APPROACH TO VALUE SETS 15 Theorem 5 Let u j m j /n Suppose l, u j O1/l uniformly for all 1 j l, and l i1 u i > c for some positive constant c Then Y n is asymptotically normal with mean and variance, respectively, equal to µ n n σ 2 n n l 1 u i, i l l l u i i 1 u i i11 u i1 Proof For O n, We have l n m i EY n n n i1 n 1 3 l 1 exp +O 1 u 2n n 2 i ju i n j i1 µ n exp 1 3 l 1 ju i +O exp ln 1 2n n 2 1 u i1 i n j µ n exp 1 3 l 1 ju i +O 2n n 2 1 u i1 i n j µ n exp 1 l 1 ju i 3 2n 1 u i n +O because j n n 2 µ nexp 1 2n i1 1+ l i1 i1 3 u i +O n + 2 note u 2 j O1/l nl In particular we have EY n 2 µ 2 n exp 1 1+ n l i1 1 u i +O n + 1, 2 nl and for of the order n, EY n µ nexp 2 2n 1+ l u i i1

16 16 ZHICHENG GAO AND QIANG WANG It follows that s n σ2 n µ n µ 2 n EY ny n 1 exp 1 n µ 2 n 1 n 1+ i1 1 l 1 1+ u i +O n nl i1 l 1 u i +O n + 1, 2 nl and EY n µ 2 n exp s n 2 Now the theorem follows from Lemma 1 References [1] A Abary, D Ghioca, and Q Wang, On permutation polynomials of prescribed shape, Finite Fields Appl , [2] B J Birch and H P F Swinnerton-Dyer, Note on a problem of Chowla, Acta Arith , [3] H Borges and R Conceicao, On the characterization of minimal value set polynomials, J Number Theory , [4] M Barot and J Peña, Estimating the size of a union of random subsets of fixed cardinality, Elemente der Mathemati , no 4, [5] B Bollobás, Random Graphs, Academic Press, 1985 [6] L Carlitz, D J Lewis, W H Mills, and E G Straus, Polynomials over finite fields with minimal value sets, Mathematia , [7] Q Cheng, J Hill and D Wan, Counting value sets: algorithms and complexity, Tenth Algorithmic Number Theory Symposium ANTS-X, 2012, University of California at San Deigo [8] S D Cohen, The distribution of polynomials over finite fields, Acta Arith , [9] FN David, Biometria , [10] P Das and G L Mullen, Value sets of polynomials over finite fields, in Finite Fields with Applications in Coding Theory, Cryptography and Related Areas, GL Mullen, H Stichtenoth, and H Tapia-Recillas, Eds, Springer, 2002, [11] A B Evans, Orthomorphism Graphs of Groups, Lecture Notes in Mathematics, Vol 1535, Springer, Berlin, 1992

17 A PROBABILISTIC APPROACH TO VALUE SETS 17 [12] Z Gao and NC Wormald, Asymptotic normality determined by high moments, and submap counts of random maps, Probab Theory Relat Fields , [13] J Gomez-Calderon and D J Madden, Polynomials with small value set over finite fields, J Number Theory , no 2, [14] R Guralnic and D Wan, Bounds for fixed point free elements in a transitive group and applications to curves over finite fields, Israel J Math , [15] W H Mills, Polynomials with minimal value sets, Pacific J Math , [16] G L Mullen, D Wan, and Q Wang, Value sets of polynomial maps over finite fields, Quart J Math , no 4, [17] G L Mullen, D Wan, and Q Wang, An index bound on value sets of polynomial maps over finite fields, Proceedings of Worshop on the Occasion of Harald Niederreiter s 70th Birthday: Applications of Algebra and Number Theory, June 23-27, 2014 [18] H Niederreiter and A Winterhof, Cyclotomic R-orthomorphisms of finite fields, Discrete Math , [19] G Turnwald, A new criterion for permutation polynomials, Finite Fields Appl , [20] D Wan, A p-adic lifting lemma and its applications to permutation polynomials, Lecture Notes in Pure and Appl Math, Marcel Deer, New Yor, Vol 141, 1992, [21] D Wan, P J S Shiue and C S Chen, Value sets of polynomials over finite fields, Proc Amer Math Soc , [22] Q Wang, Cyclotomic mapping permutation polynomials over finite fields, Sequences, Subsequences, and ConsequencesInternational Worshop, SSC 2007, Los Angeles, CA, USA, May 31 - June 2, 2007, , Lecture Notes in Comput Sci Vol 4893, Springer, Berlin, 2007 [23] Q Wang, Cyclotomy and permutation polynomials of large indices, Finite Fields Appl , [24] K S Williams, On general polynomials, Canad Math Bull , no 4, School of Mathematics and Statistics, Carleton University, 1125 Colonel By Drive, Ottawa, ON K1S 5B6, Canada address: zgao@mathcarletonca, wang@mathcarletonca