EasyChair Preprint. Computation of Some Integer Sequences in Maple

Size: px
Start display at page:

Download "EasyChair Preprint. Computation of Some Integer Sequences in Maple"

Transcription

1 EasyChair Preprit 4 Coputatio of Soe Iteger Sequeces i Maple W.L. Fa, David J. Jeffrey ad Eri Posta EasyChair preprits are iteded for rapid disseiatio of research results ad are itegrated with the rest of EasyChair. Noveber 20, 2017

2 Coputatio of Soe Iteger Sequeces i Maple W.L. Fa 1, D.J. Jeffrey 1, Eri Posta 2 1 Departet of Applied Matheatics, The Uiversity of Wester Otario, Lodo, Otario, Caada 2 Maplesoft, Waterloo wfa54@uwo.ca, deffrey@uwo.ca Abstract. We cosider soe iteger sequeces coected with cobiatorial applicatios. Specifically, we cosider Stirlig partitio ad cycle ubers, associated Stirlig partitio ad cycle ubers, ad Euleria ubers of the first ad secod ids. We cosider their evaluatio i differet cotexts. Oe cotext is the calculatio of a sigle value based o sigle iput arguets. A ore coo cotext, however, is the calculatio of a sequece of values. We copare strategies for both. Where possible, we copare with existig Maple ipleetatios. 1 Itroductio For exteded discussios of Stirlig ad Euleria ubers, we refer to [1, 2]. These ad siilar ubers arise frequetly i cobiatorial applicatios, ad have therefore bee ipleeted i several coputer algebra systes. To date, the stadard libraries of ost systes have icluded Stirlig ubers, but ot associated Stirlig ubers [3], eve though they have foud several applicatios i recet years. For exaple, they have appeared i series expasios for the Labert W fuctio [4], ad also appeared i oe for of Stirlig s series for the Gaa fuctio [2]. (Stirlig did ot defie the associated ubers.) Aother feature of ay ipleetatios is that the fuctios expect a sigle arguet, ad retur a sigle value. I practice, however, a applicatio will usually require a sequece of values, for exaple, to provide successive coefficiets i a series. The requireet of returig ultiple values has already bee recogized i soe Maple fuctios, for exaple, i the ipleetatio of Beroulli ubers: it accepts a ode paraeter. To quote fro Maple help: The ode paraeter cotrols whether or ot the beroulli routie coputes additioal Beroulli ubers i parallel with the requested oe. For exaple, if your coputer has 4 cores, the the coad beroulli(1000, sigleto=false) will copute (ad store) beroulli(1002) beroulli(1004), ad beroulli(1006). Sice i practice early all coputatios which use Beroulli ubers require ay of the, ad require the i sequece, this results i cosiderable efficiecy gais.

3 This paper addresses both the coputatio of sigle values ad of the iteger sequeces associated with the cobiatorial fuctios uder cosideratio. As a atter of teriology, we shall call a fuctio that accepts a uique arguet ad returs the correspodig uique result a sigleto fuctio, ad the correspodig operatio a sigleto coputatio. I cotrast, a fuctio acceptig a rage (explicit or iplicit) of arguets ad returig the correspodig list of values will be a sequece fuctio, ad the calculatio a sequece calculatio. 1.1 Defiitios of ubers We collect here the defiitios of all ubers cosidered. Defiitio 1. The r-associated Stirlig ubers of the first id, ore briefly Stirlig r-cycle ubers, are defied by the geeratig fuctio r 1 1 l 1 z z =! [ ] z. (1) =1 0! Rear 1. The uber [ ] gives the uber of perutatios of distict obects ito cycles, each cycle havig a iiu cardiality r [2, p 256]. Defiitio 2. The r-associated Stirlig ubers of the secod id, called ore briefly here Stirlig r-partitio ubers, are defied, usig Karaata Kuth otatio, by the geeratig fuctio r 1 e z z =! { } z. (2)! =0 0! Rear 2. The uber { } gives the uber of partitios of a set of size ito subsets, each subset havig a iiu cardiality of r [2, 5, 6]. Defiitio 3. The Euleria ubers of the first id are defied as the uber of perutatios π 1 π 2... π of {1, 2,... } that have ascets, i.e. places where π < π +1. Defiitio 4. The Euleria ubers of the secod id are defied as the uber of perutatios of the ultiset {1, 1, 2, 2,...,, } for which all ubers betwee the two occurreces of every, with 1, are greater tha, for each perutatio havig ascets, i.e. places where π < π +1. Rear 3. Note that is ot a arguet. For exaple, give the ultiset {112233}, perutatios such as or are peritted, but is ot. Aogst these peritted perutatios, we cout those with ascets. Noeclature: I [1], the ubers are called siply Euleria ubers, while the ubers are called secod-order Euleria ubers. 2

4 2 Stirlig partitio ubers The Maple 2017 ipleetatio is a sigleto fuctio, deoted stirlig2 i the cobiat pacage. It uses the forula { } = 1 ( ) ( 1). (3)! =0 For the sigleto coputatio, Table1 shows that the ties 3 are uch less usig (3). I this table, we copared the Maple fuctio stirlig2 with the ethod give below usig the recurrece relatio (7). Tiigs for a sequece calculatio, however, give i Table 2, show the ew ethod is ore efficiet. recurrece stirlig Table 1. Tiigs (sec) for geeratig a sigleto Stirlig Partitio uber. The tie usig (7) is copared with the Maple stirlig2 fuctio. 2.1 Sequece calculatio Give,, we wish to copute all Stirlig partitio ubers { i } such that i ad. We use the recurrece relatio { i } { i 1 = subect to the boudary coditios { } = 1, ad } + { i 1 1 }, (4) { } i = 1. (5) 1 Sice { i } = 0 for > i (see Fig. 1), we defie a atrix P which will ot store these zeros. { } i + 1 P i =. (6) The the recurrece relatio becoes P (i, ) = P (i 1, ) + P (i, 1). (7) The boudary coditios the becoe, respectively, P (1, ) = { } = 1, ad P (i, 1) = 1. 3 Product placeet: ties foud usig a Itel i7 i a Leovo Ultraboo. 3

5 Tiigs Table 2 shows the tiigs for fillig atrices of various sizes with iteger sequeces of Stirlig partitio ubers. The recurrece relatio (7) is copared with creatig each etry through a call to Maple s stirlig2. Fillig the square atrix P (, ) actually calculates all partitio ubers { i } with i 2 ad. This is doe for tiig coveiece, ad the atrix ca be reshaped for other applicatios. recurrece stirlig Table 2. Tiigs (sec) for geeratig sequeces of Stirlig Partitio ubers. The tie usig (7) copared with Maple stirlig2 fuctio. 3 Stirlig cycle ubers We cosider the coputatio of [ ], ipleeted i Maple 2017 as stirlig1 i the cobiat pacage. The coputatioal ethod used by stirlig1 is based o Stirlig s origial defiitio of his ubers: x = [ ] ( 1) x. (8) For give, stirlig1 costructs the product o the left, which is the collected i powers of x, so that by equatig the coefficiets of x, all ubers [ ] for 1 are deteried ad stored. Thus, a future call to [ ] with 1 will be retured by table looup, but a future call with a differet will iitiate a ew coputatio. It is iterestig that although the iterface appears to offer the user oly a sigleto coputatio, i fact a particular iteger sequece has bee coputed siletly. 3.1 Sigleto coputatio A sigleto coputatio returs the value of a fuctio for a sigle pair of iput arguets. We ipleet the ow recurrece relatio [ ] = ( 1) [ 1 4 ] + [ ] 1, (9) 1

6 subect to boudary coditios [ ] =1, for 1, (10) [ ] =( 1)!. (11) 1 We defie the vector u (i) = [ ] i + 1. I Fig. 1, we see that for fixed i, u (i) describes ubers alog the ith diagoal lie, coutig fro the left. The recurrece relatio (9) ca be writte i ters of u as u (i) = (i + 2)u (i 1) + u (i) 1, with u (i) 1 = (i 1)!. We ote that oce u (i 1) is used, it does ot eed to be stored further, so we ca overwrite storage. Our iteratio schee is thus (Maple otatio for the ith eleet of a vector is u[i]) u[] = (i + 2)u[] + u[ 1]. Therefore, we iitialize u[1] = u (i) 1 = [ i 1] = (i 1)! ad fill i diagoal lies successively. Coplexity The ai of this subsectio is to gai isight ito the best ways to test the ipleetatios, by idetifyig the worse cases for the ethods. A full bit coplexity is beyod the scope of this paper, ad will require ore wor o estiates for the sizes of Stirlig ubers. As poited out by Wilf [9], the available estiates are for [ ] whe is fixed ad, whereas the preset algoriths require owledge of the opposite case. I order to calculate the uber [ ], a vector of legth ust be recoputed (overwritte) ties. Each iteratio requires oe ultiplicatio ad 3 additios. Therefore the coplexity is ( ). We ca therefore expect that the worst case for the ethod will be = /2. Sice Maple s approach ad the preset oe calculate differet sets of ubers, a direct copariso is ot very eaigful, ad so we siply ae a brief copariso betwee oe-tie calculatios. Notice that i Table 3, the ties tae by stirlig1 are approxiately idepedet of as expected. 3.2 A fiite su For copleteess, we etio that a sigleto cycle uber ca be foud fro a fiite su, as was doe for a sigleto partitio uber. We have [ ] = =0 ( 1) + ( )( 2 ){ } +. (12) 5

7 Fig. 1. Schee for calculatig sigleto Stirlig cycle [ { ] or partitio ubers }. The coputatio proceeds fro left to right ad botto to top. At each stage oly which is progressively overwritte. The ope circles show the base of each successive loop. The blac filled circles show the recurrece relatio used. The larger circle is calculated fro the two saller oes. The triagles lie show the poits coputed by oe call to stirlig1. the ubers o oe diagoal eed to be stored i the vector u (i) stirlig1 Preset schee Table 3. Ties for a sigle call to Maple s stirlig1 ad the preset sigleto coputatio. Tiigs (sec) based o 10 trials, with eory beig cleared before each call. 6

8 Cobiig this with (3), we ca express a cycle uber as a double su. This, however, is too slow to warrat further cosideratio. 3.3 Sequece calculatio The ethod used above for partitio ubers ca be readily adapted for cycle ubers. Give,, we copute all Stirlig cycle ubers [ i ] such that i ad. We use the recurrece relatio [ ] [ ] [ ] i i 1 i 1 = (i 1) +, (13) 1 subect to the boudary coditios [ ] = 1, ad [ ] i = (i 1)!. (14) 1 Sice [ i ] = 0 for > i (see Fig. 1), we defie a atrix C which will ot store these zeros. [ ] i + 1 C i =. (15) The the boudary coditios are C(1, ) = [ ] = 1, ad C(i, 1) = (i 1)!. The recurrece relatio becoes C(i, ) = (i + 2) C(i 1, ) + C(i, 1). (16) Tiigs Table 4 shows the tiig for fillig atrices of various sizes with iteger sequeces of Stirlig cycle ubers. The recurrece relatio (16) is copared with creatig each etry through a call to Maple s stirlig1. Fillig the square atrix C(, ) actually calculated all cycle ubers [ i ] with i 2. This is doe for tiig purposes, ad the atrix ca be reshaped for other applicatios. The copariso is to copute the sae ubers usig the sequece calculatio fuctio stirlig1. Larger values of (, ) are ot tabulated because a bug i Maple 2016 (ad earlier) caused larger arguets to fail. This will be corrected i Maple Associated Stirlig ubers There are o ow aalogues of (3) or (8) for the associated Stirlig ubers for r 2; hece we ust use either the geeratig fuctios (2) ad (1), or the followig recurrece relatios. { } { } ( ){ } + 1 r + 1 = +, (17) r 1 1 [ ] + 1 [ ] = + r 1 [ r ]. (18) 7

9 recurrece stirlig Table 4. Tiigs (sec) for geeratig sequeces of Stirlig cycle ubers. The tie usig (16) copared with Maple s stirlig1 fuctio. Note that 0 = 1. The boudary cases are { } = 1, r, (19) 1 [ ] = ( 1)!, r, (20) 1 { } r = (r)! (r!)!, 1, (21) [ ] r = (r)! r!, 1. (22) 4.1 Sigleto Stirlig 2-partitio ad 2-cycle The two coputatios have the sae structure, ad ca be described i parallel. We choose to ipleet { } { } { } 1 2 = + ( 1), (23) [ ] [ ] 1 = ( 1) We also have boudary coditios We defie the vector [ ] 2 = [ ] { } 2 1 [ ] 2 + ( 1) 1 = (2)! = (2 1)!!!2 { } (2 + 1)! = 2 3( 1)!2 = 2 u (i) = [ ] i ,. (24)

10 ad siilarly for 2-partitio ubers. I Fig. 2, we see that if we fix i, the u (i) describes ubers alog the ith diagoal lie. Now u (i) 1 = i! ad u (i) = (i + 2 2)u (i 1) + (i + 2 2)u (i) 1. We ote that oce u (i 1) is used, it does ot eed to be stored further, so we ca overwrite storage. Our iteratio schee is thus u[] = (i + 2 2)(u[] + u[ 1]). For iitializatio, we ca use a special case of (24): [ ] [ ] = u (1) = (2 + 1) = (2 + 1)u (1). Therefore, we iitialize u to i = 1 usig u (1) = u[] = 1 ad fill i oe lie at a tie by fixig i ad loopig over. Each loop starts settig u (i) 1 = i! = iu (i 1) 1. We the loop over i. Fig. 2. Calculatig 2-partitio ad 2-cycle ubers. As with the r = 1 case, oly ubers o oe slopig lie eed to be ept at ay stage of the coputatio. The sae covetio for illustratig the recurrece relatio is used. 4.2 Sequece calculatio of 2-partitio ad 2-cycle ubers Give,, we copute all Stirlig 2-partitio ubers { } i or 2-cycle ubers [ i ] such that i ad. We use the recurrece relatios (23) or (24) as 9

11 appropriate. Sice { i } = [ ] i = 0 for 2 > i (see Fig. 2), we defie a atrix C which will ot store these zeros. [ ] i C i =. (25) The the recurrece relatio for 2-partitio becoes C(i, ) = C(i 1, ) + (i + 2 2)C(i 1, 1). (26) The recurrece relatio for 2-cycle becoes C(i, ) = (i + 2 2) C(i 1, ) + (i + 2 2)C(i 1, 1). (27) 4.3 Sigleto Stirlig r-partitio ad r-cycle ubers Fro the above discussio of 1-associated ad 2-associated ubers, the geeralizatio is clear. We have to ipleet [ ] [ ] [ ] + 1 r + 1 = + ( 1)( 2)... ( r + 2). (28) 1 We defie the vector u (i) = [ ] i + r 1 The geeralizatio of Fig. 2 to oe cotaiig lies of slope 1/r is ot show. For fixed i, u (i) describes ubers alog oe of the lies, with u (i) 1 = (i + r 2)! ad u (i) = (i + r 2)u (i 1). + (i + r 2)(i + r 3)... (i + r r)u (i) 1. We ote that oce u (i 1) is used, it does ot eed to be stored further, so we ca overwrite storage. Our iteratio schee is thus u[] = (i + r 2)u[] + (i + r 2)... (i + r r)u[ 1]. For iitializatio, we ca use a special case of (28): [ ] [ ] r + r = u (1) r =(r + r 1)(r + r 2)... (r + 1) =(r + r 1)(r + r 2)... (r + 1)u (1). We iitialize u usig u (1) = u[] = 1 ad fill i each lie by fixig i ad loopig over. We start each loop by settig u (i) 1 = (i + r 2)! = (i + r 2)u (i 1) 1. We the loop over i. 10

12 4.4 Sequece calculatio of r-partitio ad r-cycle ubers Give,, we copute all Stirlig r-partitio ubers { } i or r-cycle ubers [ i ] such that i ad. The recurrece relatio applied here ca refer to (17) ad (18), which is subect to the boudary coditios Sice { i } = [ ] i these zeros. { } 1 = 1 [ ] 1 = 1. (29) 1 = 0 for r > i, we defie a atrix C which will ot store [ ] i + r 1 C i =. (30) The the recurrece relatio for r-partitio becoes ( ) i + r 2 C(i, ) = C(i 1, ) + C(i r + 1, 1). (31) r 1 The recurrece relatio for r-cycle becoes C(i, ) = (i+r 2) C(i 1, )+(i+r 2)(i+r 3)... (i+r r) C(i r+1, 1) (32) 4.5 Ipleetatio i Maple I our ipleetatio of Stirlig ubers, we provide procedures for users to copute either a sigleto Stirlig uber or a sequece of Stirlig ubers. The procedures are 1. StirligRCycle: to calculate a sigleto Stirlig r-cycle uber. 2. StirligRCycleMatrix: to calculate a sequece of Stirlig r-cycle ubers. 3. StirligRPartitio: to calculate a sigleto Stirlig r-partitio uber. 4. StirligRPartitioMatrix: to calculate a sequece of Stirlig r-partitio ubers. Neither Maple or Matheatica has a ipleetatio with which to copare our progras. Therefore we have prograed the recurrece relatios, as well as the geeratig fuctios i Maple. I Table 5 below, we copared our ew schee for coputig a sigleto r-associated Stirlig cycle uber with usig the geeratig fuctio. The geeratig fuctio for Stirlig r-cycle ubers is: StirRCycleGe := proc(,,r) local t, z, p; t:=series((l(1/(1-z)) - add(z^p/p, p=1..r-1))^, z=0,+1);!*coeff(t, z, )/!; ed proc; 11

13 Sigleto schee Geeratig fuctio Table 5. Tiigs i secods of coputatios of sigle Stirlig r-cycle uber. Colu headigs give the fuctios used. The ubers tested were [ ] Table 5 shows that the sigleto schee is uch faster tha the geeratig fuctio for the coputatio of sigle r-associated Stirlig cycle uber. For the coputatio of a sequece of r-associated Stirlig cycle ubers, we copared three ethods: (1) a loop callig the sigleto fuctio; (2) a loop callig the geeratig fuctio; (3) the sequece procedure. The results are collected i Tables 6 ad 7, ad show that the sequece procedure is fastest. Sigleto schee Geeratig fuctio Table 6. Tiigs i secods of coputatios of a sequece of r-associated Stirlig cycle ubers. Colu headigs give the fuctios used. The iput arguet is, ad the retur is a atrix. Siilar tests were perfored for r-associated Stirlig partitio ubers. The geeratig fuctio for Stirlig r-partitio ubers is: StirRPartGe := proc(,,r) local t, z, p; t:=series((exp(z) - add(z^p/p!, p=0..r-1))^, z=0, +1);!*coeff(t, z, )/!; ed proc; The test data are collected i Tables 8, 9, 10. Sice the patter is siilar to that for cycle ubers, the discussio ad tables are abbreviated. 12

14 Sigleto schee Sequece schee Table 7. Tiigs i secods of coputatios of a sequece of r-associated Stirlig cycle ubers. Colu headigs give the fuctios used. The iput arguet is, ad the retur is a atrix. Sigleto schee Geeratig fuctio Table 8. Tiigs i secods of coputatios of sigle r-associated Stirlig partitio uber. Colu headigs give the fuctios beig used. The ubers tested were { 1000 } 18. Sigleto schee Geeratig fuctio Table 9. Tiigs i secods of coputatios of a sequece of r-associated Stirlig partitio uber. Colu headigs give the fuctios beig used. The iput arguet is, ad the retur is a atrix. Sigleto schee Sequece schee Table 10. Tiigs i secods of coputatios of a sequece of r-associated Stirlig partitio ubers. Colu headigs give the fuctios used. The iput arguet is, ad the retur is a atrix. 13

15 5 A ultiple threads approach to sequece calculatios The Maple help for Beroulli ubers, quoted i the itroductio, states that additioal values of Beroulli ubers are calculated i parallel. This sectio explored ways i which parallel coputatio could be applied to Stirlig ubers. For this, we use the Threads pacage i Maple. Whe we geerate the ubers iside a atrix, istead of fillig the atrix row by row ad colu by colu, we fill each diagoal fro left to right. Here is the ai part i the sequetial code to fill Stirlig r-cycle ubers i the atrix by diagoal with give iput arguets (, r) where is the size of atrix. for N fro 3 to do for fro 2 to N-1 do pd := ul(n-+r-1-l, l = 1.. r-1); A(N-+r, ) := pd*a(n-, -1)+(N-+r-2)*A(N-+r-1, ); ed do; ed do; Accordig to the recurrece relatio, we ow that we ca divide such diagoal ito a left half ad right half. So we defie two subrouties accordigly. fileft := proc (N, r) local, Nsplit, pd, l; global A; Nsplit := floor((1/2)*n+1/2); for fro 2 to Nsplit do pd := ul(n-+r-1-l, l = 1.. r-1); A(N-+r, ) := pd*a(n-, -1)+(N-+r-2)*A(N-+r-1, ); ed do; ed proc; ad filrght := proc (N, r) local, Nsplit, pd, l; global A; Nsplit := floor((1/2)*n+1/2); for fro Nsplit+1 to N-1 do pd := ul(n-+r-1-l, l = 1.. r-1); A(N-+r, ) := pd*a(n-, -1)+(N-+r-2)*A(N-+r-1, ); ed do; ed proc; Ad for each half of the diagoal, we ca establish a idepedet thread to fulfill the tas. We ipleeted this approach i Maple. Threaded := proc (, r) local N,, Nsplit; global A; A := Matrix(,, fill = 0); A(1, 1) := 1; for N fro 3 to do Threads:-Tas:-Start(ull, Tas = [fileft, N, r], Tas = [filrght, N, r]) ed do; ed proc; 14

16 Table 11 copares the threaded schee with the sequetial schee i the coputatio of a atrix of Stirlig cycle ubers. The table reflects the liitatio that there is a overhead cost to settig up ew threads, ad the beefit of the threaded approach is felt oly whe the aout of wor achieved withi a thread outweighs the overhead. I this ipleetatio, ew threads are created for each loop. We are explorig ew ethods of calculatio which will allow the threads to wor ore efficietly, with less overhead. Threaded schee Sequetial schee Table 11. Tiigs i secods of copariso of threaded code with sequetial code i geeratig sequeces of Stirlig Cycle ubers. The tests were ade o a AMD 8-core processor. 6 Ipleetatio of Euleria ubers The Euleria ubers share ay siilarities with the Stirlig ubers, ad all the ethods described above ca be applied to their case. The ubers obey the followig recurrece relatios [1]. 1 1 = ( + 1) + ( ), (33) 1 = ( + 1) 1 + (2 1) 1 1. (34) The preset Maple fuctios euleria1 ad euleria2 are recursively prograed ipleetatios of these equatios. As a cosequece, they are very slow for large arguets. The ew ipleetatio of these ubers cosists of 4 fuctios, which follow the patters of the Stirlig uber ipleetatios. 1. Euleria1: calculates a sigleto Euleria uber of the first id. As with Stirlig partitio ubers, a fiite su is ow which is distictly the fastest ethod for a sigleto coputatio [8]: +1 ( ) + 1 = ( 1) ( + 1). (35) =0 15

17 2. Euleria1Matrix: calculates a sequece of Euleria ubers of the first id. This follows the sequece calculatio of Stirlig ubers, usig (33). 3. Euleria2: calculates sigleto Euleria ubers of the secod id. This follows the sipleto ethod used earlier for Stirlig cycle ubers. 4. Euleria2Matrix: calculates a sequece of Euleria ubers of the secod id. This follows the sequece calculatio of Stirlig ubers. 6.1 Tiigs for Euleria uber calculatios I view of the siilarities with Stirlig ubers, we shall ot labour the coparisos betwee ethods, sice they for the sae processio of speeds see before. Table 12 copares the ew ipleetatios, followig the patters set above. Euleria1 Euleria1Matrix Euleria2 Euleria2Matrix Table 12. Tiigs i secods of coputatios of Euleria ubers. Colu headigs give the fuctios used. Refereces 1. Graha, R.L., Kuth, D.E., Patashi, O., Cocrete Matheatics, Addiso- Wesley Publishig Co., Readig, Massachusetts, Cotet, L., Advaced Cobiatorics, D. Reidel Publishig Co., Dordrecht, Hollad, Howard, F. T., Associated Stirlig Nubers, The Fiboacci Quarterly, Vol. 18(4), , Corless, R.M., Jeffrey, D.J., Kuth, D.E., A Sequece of Series for the Labert W Fuctio, Proceedigs of ISSAC 1997, ed. W.W. Kuechli, ACM Press, Karaata, J., Theoree sur la soabilite expoetielle et d autres soabilites rattachat, Matheatica, Clu, Roaia, vol. 9, , Kuth, D.E., Two Notes o Notatio, The Aerica Matheatical Mothly, vol. 99, , Stirlig, J., Methodus Differetialis, Lodo, Leher, D. H., Geeralized Euleria Nubers, Joural of Cobiatorial Theory, Series A 32, , Wilf, H. S., The Asyptotic Behavior of the Stirlig Nubers of the First Kid, Joural of Cobiatorial Theory, Series A 64, ,

Binomial transform of products

Binomial transform of products Jauary 02 207 Bioial trasfor of products Khristo N Boyadzhiev Departet of Matheatics ad Statistics Ohio Norther Uiversity Ada OH 4580 USA -boyadzhiev@ouedu Abstract Give the bioial trasfors { b } ad {

More information

Automated Proofs for Some Stirling Number Identities

Automated Proofs for Some Stirling Number Identities Autoated Proofs for Soe Stirlig Nuber Idetities Mauel Kauers ad Carste Scheider Research Istitute for Sybolic Coputatio Johaes Kepler Uiversity Altebergerstraße 69 A4040 Liz, Austria Subitted: Sep 1, 2007;

More information

Double Derangement Permutations

Double Derangement Permutations Ope Joural of iscrete Matheatics, 206, 6, 99-04 Published Olie April 206 i SciRes http://wwwscirporg/joural/ojd http://dxdoiorg/04236/ojd2066200 ouble erageet Perutatios Pooya aeshad, Kayar Mirzavaziri

More information

Chapter 2. Asymptotic Notation

Chapter 2. Asymptotic Notation Asyptotic Notatio 3 Chapter Asyptotic Notatio Goal : To siplify the aalysis of ruig tie by gettig rid of details which ay be affected by specific ipleetatio ad hardware. [1] The Big Oh (O-Notatio) : It

More information

Discrete Mathematics: Lectures 8 and 9 Principle of Inclusion and Exclusion Instructor: Arijit Bishnu Date: August 11 and 13, 2009

Discrete Mathematics: Lectures 8 and 9 Principle of Inclusion and Exclusion Instructor: Arijit Bishnu Date: August 11 and 13, 2009 Discrete Matheatics: Lectures 8 ad 9 Priciple of Iclusio ad Exclusio Istructor: Arijit Bishu Date: August ad 3, 009 As you ca observe by ow, we ca cout i various ways. Oe such ethod is the age-old priciple

More information

Bertrand s postulate Chapter 2

Bertrand s postulate Chapter 2 Bertrad s postulate Chapter We have see that the sequece of prie ubers, 3, 5, 7,... is ifiite. To see that the size of its gaps is ot bouded, let N := 3 5 p deote the product of all prie ubers that are

More information

A Pair of Operator Summation Formulas and Their Applications

A Pair of Operator Summation Formulas and Their Applications A Pair of Operator Suatio Forulas ad Their Applicatios Tia-Xiao He 1, Leetsch C. Hsu, ad Dogsheg Yi 3 1 Departet of Matheatics ad Coputer Sciece Illiois Wesleya Uiversity Blooigto, IL 6170-900, USA Departet

More information

Math 4707 Spring 2018 (Darij Grinberg): midterm 2 page 1. Math 4707 Spring 2018 (Darij Grinberg): midterm 2 with solutions [preliminary version]

Math 4707 Spring 2018 (Darij Grinberg): midterm 2 page 1. Math 4707 Spring 2018 (Darij Grinberg): midterm 2 with solutions [preliminary version] Math 4707 Sprig 08 Darij Griberg: idter page Math 4707 Sprig 08 Darij Griberg: idter with solutios [preliiary versio] Cotets 0.. Coutig first-eve tuples......................... 3 0.. Coutig legal paths

More information

Some results on the Apostol-Bernoulli and Apostol-Euler polynomials

Some results on the Apostol-Bernoulli and Apostol-Euler polynomials Soe results o the Apostol-Beroulli ad Apostol-Euler polyoials Weipig Wag a, Cagzhi Jia a Tiaig Wag a, b a Departet of Applied Matheatics, Dalia Uiversity of Techology Dalia 116024, P. R. Chia b Departet

More information

6.4 Binomial Coefficients

6.4 Binomial Coefficients 64 Bioial Coefficiets Pascal s Forula Pascal s forula, aed after the seveteeth-cetury Frech atheaticia ad philosopher Blaise Pascal, is oe of the ost faous ad useful i cobiatorics (which is the foral ter

More information

19.1 The dictionary problem

19.1 The dictionary problem CS125 Lecture 19 Fall 2016 19.1 The dictioary proble Cosider the followig data structural proble, usually called the dictioary proble. We have a set of ites. Each ite is a (key, value pair. Keys are i

More information

distinct distinct n k n k n! n n k k n 1 if k n, identical identical p j (k) p 0 if k > n n (k)

distinct distinct n k n k n! n n k k n 1 if k n, identical identical p j (k) p 0 if k > n n (k) THE TWELVEFOLD WAY FOLLOWING GIAN-CARLO ROTA How ay ways ca we distribute objects to recipiets? Equivaletly, we wat to euerate equivalece classes of fuctios f : X Y where X = ad Y = The fuctios are subject

More information

The Hypergeometric Coupon Collection Problem and its Dual

The Hypergeometric Coupon Collection Problem and its Dual Joural of Idustrial ad Systes Egieerig Vol., o., pp -7 Sprig 7 The Hypergeoetric Coupo Collectio Proble ad its Dual Sheldo M. Ross Epstei Departet of Idustrial ad Systes Egieerig, Uiversity of Souther

More information

Bernoulli Polynomials Talks given at LSBU, October and November 2015 Tony Forbes

Bernoulli Polynomials Talks given at LSBU, October and November 2015 Tony Forbes Beroulli Polyoials Tals give at LSBU, October ad Noveber 5 Toy Forbes Beroulli Polyoials The Beroulli polyoials B (x) are defied by B (x), Thus B (x) B (x) ad B (x) x, B (x) x x + 6, B (x) dx,. () B 3

More information

CS 70 Second Midterm 7 April NAME (1 pt): SID (1 pt): TA (1 pt): Name of Neighbor to your left (1 pt): Name of Neighbor to your right (1 pt):

CS 70 Second Midterm 7 April NAME (1 pt): SID (1 pt): TA (1 pt): Name of Neighbor to your left (1 pt): Name of Neighbor to your right (1 pt): CS 70 Secod Midter 7 April 2011 NAME (1 pt): SID (1 pt): TA (1 pt): Nae of Neighbor to your left (1 pt): Nae of Neighbor to your right (1 pt): Istructios: This is a closed book, closed calculator, closed

More information

(s)h(s) = K( s + 8 ) = 5 and one finite zero is located at z 1

(s)h(s) = K( s + 8 ) = 5 and one finite zero is located at z 1 ROOT LOCUS TECHNIQUE 93 should be desiged differetly to eet differet specificatios depedig o its area of applicatio. We have observed i Sectio 6.4 of Chapter 6, how the variatio of a sigle paraeter like

More information

Chapter 9 Computation of the Discrete. Fourier Transform

Chapter 9 Computation of the Discrete. Fourier Transform Chapter 9 Coputatio of the Discrete Fourier Trasfor Itroductio Efficiet Coputatio of the Discrete Fourier Trasfor Goertzel Algorith Deciatio-I-Tie FFT Algoriths Deciatio-I-Frequecy FFT Algoriths Ipleetatio

More information

q-fibonacci polynomials and q-catalan numbers Johann Cigler [ ] (4) I don t know who has observed this well-known fact for the first time.

q-fibonacci polynomials and q-catalan numbers Johann Cigler [ ] (4) I don t know who has observed this well-known fact for the first time. -Fiboacci polyoials ad -Catala ubers Joha Cigler The Fiboacci polyoials satisfy the recurrece F ( x s) = s x = () F ( x s) = xf ( x s) + sf ( x s) () with iitial values F ( x s ) = ad F( x s ) = These

More information

Lecture 11. Solution of Nonlinear Equations - III

Lecture 11. Solution of Nonlinear Equations - III Eiciecy o a ethod Lecture Solutio o Noliear Equatios - III The eiciecy ide o a iterative ethod is deied by / E r r: rate o covergece o the ethod : total uber o uctios ad derivative evaluatios at each step

More information

Generalized Fibonacci-Like Sequence and. Fibonacci Sequence

Generalized Fibonacci-Like Sequence and. Fibonacci Sequence It. J. Cotep. Math. Scieces, Vol., 04, o., - 4 HIKARI Ltd, www.-hiari.co http://dx.doi.org/0.88/ijcs.04.48 Geeralized Fiboacci-Lie Sequece ad Fiboacci Sequece Sajay Hare epartet of Matheatics Govt. Holar

More information

) is a square matrix with the property that for any m n matrix A, the product AI equals A. The identity matrix has a ii

) is a square matrix with the property that for any m n matrix A, the product AI equals A. The identity matrix has a ii square atrix is oe that has the sae uber of rows as colus; that is, a atrix. he idetity atrix (deoted by I, I, or [] I ) is a square atrix with the property that for ay atrix, the product I equals. he

More information

AVERAGE MARKS SCALING

AVERAGE MARKS SCALING TERTIARY INSTITUTIONS SERVICE CENTRE Level 1, 100 Royal Street East Perth, Wester Australia 6004 Telephoe (08) 9318 8000 Facsiile (08) 95 7050 http://wwwtisceduau/ 1 Itroductio AVERAGE MARKS SCALING I

More information

Integrals of Functions of Several Variables

Integrals of Functions of Several Variables Itegrals of Fuctios of Several Variables We ofte resort to itegratios i order to deterie the exact value I of soe quatity which we are uable to evaluate by perforig a fiite uber of additio or ultiplicatio

More information

#A18 INTEGERS 11 (2011) THE (EXPONENTIAL) BIPARTITIONAL POLYNOMIALS AND POLYNOMIAL SEQUENCES OF TRINOMIAL TYPE: PART I

#A18 INTEGERS 11 (2011) THE (EXPONENTIAL) BIPARTITIONAL POLYNOMIALS AND POLYNOMIAL SEQUENCES OF TRINOMIAL TYPE: PART I #A18 INTEGERS 11 (2011) THE (EXPONENTIAL) BIPARTITIONAL POLYNOMIALS AND POLYNOMIAL SEQUENCES OF TRINOMIAL TYPE: PART I Miloud Mihoubi 1 Uiversité des Scieces et de la Techologie Houari Bouediee Faculty

More information

Probabilistic Analysis of Rectilinear Steiner Trees

Probabilistic Analysis of Rectilinear Steiner Trees Probabilistic Aalysis of Rectiliear Steier Trees Chuhog Che Departet of Electrical ad Coputer Egieerig Uiversity of Widsor, Otario, Caada, N9B 3P4 E-ail: cche@uwidsor.ca Abstract Steier tree is a fudaetal

More information

CALCULATION OF FIBONACCI VECTORS

CALCULATION OF FIBONACCI VECTORS CALCULATION OF FIBONACCI VECTORS Stuart D. Aderso Departmet of Physics, Ithaca College 953 Daby Road, Ithaca NY 14850, USA email: saderso@ithaca.edu ad Dai Novak Departmet of Mathematics, Ithaca College

More information

On an Algorithm for Isomorphism-Free Generations of Combinatorial Objects

On an Algorithm for Isomorphism-Free Generations of Combinatorial Objects O a Algorith for Isoorphis-Free Geeratios of Cobiatorial Objects Krasiir Yordzhev Faculty of Matheatics ad Natural Scieces South-West Uiversity, Blagoevgrad, Bulgaria yordzhev@swubg Abstract: I the wor

More information

2. F ; =(,1)F,1; +F,1;,1 is satised by thestirlig ubers of the rst kid ([1], p. 824). 3. F ; = F,1; + F,1;,1 is satised by the Stirlig ubers of the se

2. F ; =(,1)F,1; +F,1;,1 is satised by thestirlig ubers of the rst kid ([1], p. 824). 3. F ; = F,1; + F,1;,1 is satised by the Stirlig ubers of the se O First-Order Two-Diesioal Liear Hoogeeous Partial Dierece Equatios G. Neil Have y Ditri A. Gusev z Abstract Aalysis of algoriths occasioally requires solvig of rst-order two-diesioal liear hoogeeous partial

More information

COMP 2804 Solutions Assignment 1

COMP 2804 Solutions Assignment 1 COMP 2804 Solutios Assiget 1 Questio 1: O the first page of your assiget, write your ae ad studet uber Solutio: Nae: Jaes Bod Studet uber: 007 Questio 2: I Tic-Tac-Toe, we are give a 3 3 grid, cosistig

More information

Axioms of Measure Theory

Axioms of Measure Theory MATH 532 Axioms of Measure Theory Dr. Neal, WKU I. The Space Throughout the course, we shall let X deote a geeric o-empty set. I geeral, we shall ot assume that ay algebraic structure exists o X so that

More information

subcaptionfont+=small,labelformat=parens,labelsep=space,skip=6pt,list=0,hypcap=0 subcaption ALGEBRAIC COMBINATORICS LECTURE 8 TUESDAY, 2/16/2016

subcaptionfont+=small,labelformat=parens,labelsep=space,skip=6pt,list=0,hypcap=0 subcaption ALGEBRAIC COMBINATORICS LECTURE 8 TUESDAY, 2/16/2016 subcaptiofot+=small,labelformat=pares,labelsep=space,skip=6pt,list=0,hypcap=0 subcaptio ALGEBRAIC COMBINATORICS LECTURE 8 TUESDAY, /6/06. Self-cojugate Partitios Recall that, give a partitio λ, we may

More information

x !1! + 1!2!

x !1! + 1!2! 4 Euler-Maclauri Suatio Forula 4. Beroulli Nuber & Beroulli Polyoial 4.. Defiitio of Beroulli Nuber Beroulli ubers B (,,3,) are defied as coefficiets of the followig equatio. x e x - B x! 4.. Expreesio

More information

Evaluation of Bessel Functions Using a Computer Program

Evaluation of Bessel Functions Using a Computer Program Evaluatio of Bessel Fuctios Usig a Coputer Progra P. S. Yeh, Ph.D. Abstract I cylidrical coordiate, there are two types of Bessel fuctios. These fuctios are the Bessel fuctio ad the odified Bessel fuctio.

More information

On the Fibonacci-like Sequences of Higher Order

On the Fibonacci-like Sequences of Higher Order Article Iteratioal Joural of oder atheatical Scieces, 05, 3(): 5-59 Iteratioal Joural of oder atheatical Scieces Joural hoepage: wwwoderscietificpressco/jourals/ijsaspx O the Fiboacci-like Sequeces of

More information

Math 155 (Lecture 3)

Math 155 (Lecture 3) Math 55 (Lecture 3) September 8, I this lecture, we ll cosider the aswer to oe of the most basic coutig problems i combiatorics Questio How may ways are there to choose a -elemet subset of the set {,,,

More information

Commutativity in Permutation Groups

Commutativity in Permutation Groups Commutativity i Permutatio Groups Richard Wito, PhD Abstract I the group Sym(S) of permutatios o a oempty set S, fixed poits ad trasiet poits are defied Prelimiary results o fixed ad trasiet poits are

More information

A Stirling Encounter with Harmonic Numbers

A Stirling Encounter with Harmonic Numbers VOL 75, NO, APRIL 00 95 A Stirlig Ecouter with Haroic Nubers ARTHUR T BENJAMIN GREGORY O PRESTON Harvey Mudd College Clareot, CA 91711 bejai@hcedu gpresto@hcedu JENNIFER J QUINN Occidetal College 1600

More information

Recursive Algorithms. Recurrences. Recursive Algorithms Analysis

Recursive Algorithms. Recurrences. Recursive Algorithms Analysis Recursive Algorithms Recurreces Computer Sciece & Egieerig 35: Discrete Mathematics Christopher M Bourke cbourke@cseuledu A recursive algorithm is oe i which objects are defied i terms of other objects

More information

A string of not-so-obvious statements about correlation in the data. (This refers to the mechanical calculation of correlation in the data.

A string of not-so-obvious statements about correlation in the data. (This refers to the mechanical calculation of correlation in the data. STAT-UB.003 NOTES for Wedesday 0.MAY.0 We will use the file JulieApartet.tw. We ll give the regressio of Price o SqFt, show residual versus fitted plot, save residuals ad fitted. Give plot of (Resid, Price,

More information

Recurrence Relations

Recurrence Relations Recurrece Relatios Aalysis of recursive algorithms, such as: it factorial (it ) { if (==0) retur ; else retur ( * factorial(-)); } Let t be the umber of multiplicatios eeded to calculate factorial(). The

More information

NICK DUFRESNE. 1 1 p(x). To determine some formulas for the generating function of the Schröder numbers, r(x) = a(x) =

NICK DUFRESNE. 1 1 p(x). To determine some formulas for the generating function of the Schröder numbers, r(x) = a(x) = AN INTRODUCTION TO SCHRÖDER AND UNKNOWN NUMBERS NICK DUFRESNE Abstract. I this article we will itroduce two types of lattice paths, Schröder paths ad Ukow paths. We will examie differet properties of each,

More information

JORGE LUIS AROCHA AND BERNARDO LLANO. Average atchig polyoial Cosider a siple graph G =(V E): Let M E a atchig of the graph G: If M is a atchig, the a

JORGE LUIS AROCHA AND BERNARDO LLANO. Average atchig polyoial Cosider a siple graph G =(V E): Let M E a atchig of the graph G: If M is a atchig, the a MEAN VALUE FOR THE MATCHING AND DOMINATING POLYNOMIAL JORGE LUIS AROCHA AND BERNARDO LLANO Abstract. The ea value of the atchig polyoial is coputed i the faily of all labeled graphs with vertices. We dee

More information

Exercise 8 CRITICAL SPEEDS OF THE ROTATING SHAFT

Exercise 8 CRITICAL SPEEDS OF THE ROTATING SHAFT Exercise 8 CRITICA SEEDS OF TE ROTATING SAFT. Ai of the exercise Observatio ad easureet of three cosecutive critical speeds ad correspodig odes of the actual rotatig shaft. Copariso of aalytically coputed

More information

Sequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence

Sequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece 1, 1, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet

More information

The Rand and block distances of pairs of set partitions

The Rand and block distances of pairs of set partitions The Rad ad block distaces of pairs of set partitios Frak Ruskey 1 ad Jeifer Woodcock 1 Dept. of Computer Sciece, Uiversity of Victoria, CANADA Abstract. The Rad distaces of two set partitios is the umber

More information

Jacobi symbols. p 1. Note: The Jacobi symbol does not necessarily distinguish between quadratic residues and nonresidues. That is, we could have ( a

Jacobi symbols. p 1. Note: The Jacobi symbol does not necessarily distinguish between quadratic residues and nonresidues. That is, we could have ( a Jacobi sybols efiitio Let be a odd positive iteger If 1, the Jacobi sybol : Z C is the costat fuctio 1 1 If > 1, it has a decopositio ( as ) a product of (ot ecessarily distict) pries p 1 p r The Jacobi

More information

Ma/CS 6a Class 22: Power Series

Ma/CS 6a Class 22: Power Series Ma/CS 6a Class 22: Power Series By Ada Sheffer Power Series Mooial: ax i. Polyoial: a 0 + a 1 x + a 2 x 2 + + a x. Power series: A x = a 0 + a 1 x + a 2 x 2 + Also called foral power series, because we

More information

The Binomial Multi-Section Transformer

The Binomial Multi-Section Transformer 4/15/2010 The Bioial Multisectio Matchig Trasforer preset.doc 1/24 The Bioial Multi-Sectio Trasforer Recall that a ulti-sectio atchig etwork ca be described usig the theory of sall reflectios as: where:

More information

Counting Well-Formed Parenthesizations Easily

Counting Well-Formed Parenthesizations Easily Coutig Well-Formed Parethesizatios Easily Pekka Kilpeläie Uiversity of Easter Filad School of Computig, Kuopio August 20, 2014 Abstract It is well kow that there is a oe-to-oe correspodece betwee ordered

More information

Define a Markov chain on {1,..., 6} with transition probability matrix P =

Define a Markov chain on {1,..., 6} with transition probability matrix P = Pla Group Work 0. The title says it all Next Tie: MCMC ad Geeral-state Markov Chais Midter Exa: Tuesday 8 March i class Hoework 4 due Thursday Uless otherwise oted, let X be a irreducible, aperiodic Markov

More information

X. Perturbation Theory

X. Perturbation Theory X. Perturbatio Theory I perturbatio theory, oe deals with a ailtoia that is coposed Ĥ that is typically exactly solvable of two pieces: a referece part ad a perturbatio ( Ĥ ) that is assued to be sall.

More information

42 Dependence and Bases

42 Dependence and Bases 42 Depedece ad Bases The spa s(a) of a subset A i vector space V is a subspace of V. This spa ay be the whole vector space V (we say the A spas V). I this paragraph we study subsets A of V which spa V

More information

Wavelet Transform Theory. Prof. Mark Fowler Department of Electrical Engineering State University of New York at Binghamton

Wavelet Transform Theory. Prof. Mark Fowler Department of Electrical Engineering State University of New York at Binghamton Wavelet Trasfor Theory Prof. Mark Fowler Departet of Electrical Egieerig State Uiversity of New York at Bighato What is a Wavelet Trasfor? Decopositio of a sigal ito costituet parts Note that there are

More information

Recursive Algorithm for Generating Partitions of an Integer. 1 Preliminary

Recursive Algorithm for Generating Partitions of an Integer. 1 Preliminary Recursive Algorithm for Geeratig Partitios of a Iteger Sug-Hyuk Cha Computer Sciece Departmet, Pace Uiversity 1 Pace Plaza, New York, NY 10038 USA scha@pace.edu Abstract. This article first reviews the

More information

Sums Involving Moments of Reciprocals of Binomial Coefficients

Sums Involving Moments of Reciprocals of Binomial Coefficients 1 2 3 47 6 23 11 Joural of Iteger Sequeces, Vol. 14 (2011, Article 11.6.6 Sums Ivolvig Momets of Reciprocals of Biomial Coefficiets Hacèe Belbachir ad Mourad Rahmai Uiversity of Scieces ad Techology Houari

More information

Classification of problem & problem solving strategies. classification of time complexities (linear, logarithmic etc)

Classification of problem & problem solving strategies. classification of time complexities (linear, logarithmic etc) Classificatio of problem & problem solvig strategies classificatio of time complexities (liear, arithmic etc) Problem subdivisio Divide ad Coquer strategy. Asymptotic otatios, lower boud ad upper boud:

More information

MT5821 Advanced Combinatorics

MT5821 Advanced Combinatorics MT5821 Advaced Combiatorics 9 Set partitios ad permutatios It could be said that the mai objects of iterest i combiatorics are subsets, partitios ad permutatios of a fiite set. We have spet some time coutig

More information

FACTORS OF SUMS OF POWERS OF BINOMIAL COEFFICIENTS. Dedicated to the memory of Paul Erdős. 1. Introduction. n k. f n,a =

FACTORS OF SUMS OF POWERS OF BINOMIAL COEFFICIENTS. Dedicated to the memory of Paul Erdős. 1. Introduction. n k. f n,a = FACTORS OF SUMS OF POWERS OF BINOMIAL COEFFICIENTS NEIL J. CALKIN Abstract. We prove divisibility properties for sus of powers of bioial coefficiets ad of -bioial coefficiets. Dedicated to the eory of

More information

Lecture 19. Curve fitting I. 1 Introduction. 2 Fitting a constant to measured data

Lecture 19. Curve fitting I. 1 Introduction. 2 Fitting a constant to measured data Lecture 9 Curve fittig I Itroductio Suppose we are preseted with eight poits of easured data (x i, y j ). As show i Fig. o the left, we could represet the uderlyig fuctio of which these data are saples

More information

5.6 Binomial Multi-section Matching Transformer

5.6 Binomial Multi-section Matching Transformer 4/14/21 5_6 Bioial Multisectio Matchig Trasforers 1/1 5.6 Bioial Multi-sectio Matchig Trasforer Readig Assiget: pp. 246-25 Oe way to axiize badwidth is to costruct a ultisectio Γ f that is axially flat.

More information

Harmonic Number Identities Via Euler s Transform

Harmonic Number Identities Via Euler s Transform 1 2 3 47 6 23 11 Joural of Iteger Sequeces, Vol. 12 2009), Article 09.6.1 Harmoic Number Idetities Via Euler s Trasform Khristo N. Boyadzhiev Departmet of Mathematics Ohio Norther Uiversity Ada, Ohio 45810

More information

The Binomial Multi- Section Transformer

The Binomial Multi- Section Transformer 4/4/26 The Bioial Multisectio Matchig Trasforer /2 The Bioial Multi- Sectio Trasforer Recall that a ulti-sectio atchig etwork ca be described usig the theory of sall reflectios as: where: ( ω ) = + e +

More information

Hoggatt and King [lo] defined a complete sequence of natural numbers

Hoggatt and King [lo] defined a complete sequence of natural numbers REPRESENTATIONS OF N AS A SUM OF DISTINCT ELEMENTS FROM SPECIAL SEQUENCES DAVID A. KLARNER, Uiversity of Alberta, Edmoto, Caada 1. INTRODUCTION Let a, I deote a sequece of atural umbers which satisfies

More information

TR/46 OCTOBER THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION A. TALBOT

TR/46 OCTOBER THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION A. TALBOT TR/46 OCTOBER 974 THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION by A. TALBOT .. Itroductio. A problem i approximatio theory o which I have recetly worked [] required for its solutio a proof that the

More information

Sequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018

Sequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018 CSE 353 Discrete Computatioal Structures Sprig 08 Sequeces, Mathematical Iductio, ad Recursio (Chapter 5, Epp) Note: some course slides adopted from publisher-provided material Overview May mathematical

More information

Generating Functions and Their Applications

Generating Functions and Their Applications Geeratig Fuctios ad Their Applicatios Agustius Peter Sahaggau MIT Matheatics Departet Class of 2007 18.104 Ter Paper Fall 2006 Abstract. Geeratig fuctios have useful applicatios i ay fields of study. I

More information

A new sequence convergent to Euler Mascheroni constant

A new sequence convergent to Euler Mascheroni constant You ad Che Joural of Iequalities ad Applicatios 08) 08:7 https://doi.org/0.86/s3660-08-670-6 R E S E A R C H Ope Access A ew sequece coverget to Euler Mascheroi costat Xu You * ad Di-Rog Che * Correspodece:

More information

A symmetrical Eulerian identity

A symmetrical Eulerian identity Joural of Combiatorics Volume 17, Number 1, 29 38, 2010 A symmetrical Euleria idetity Fa Chug, Ro Graham ad Do Kuth We give three proofs for the followig symmetrical idetity ivolvig biomial coefficiets

More information

Math 2784 (or 2794W) University of Connecticut

Math 2784 (or 2794W) University of Connecticut ORDERS OF GROWTH PAT SMITH Math 2784 (or 2794W) Uiversity of Coecticut Date: Mar. 2, 22. ORDERS OF GROWTH. Itroductio Gaiig a ituitive feel for the relative growth of fuctios is importat if you really

More information

Orthogonal Functions

Orthogonal Functions Royal Holloway Uiversity of odo Departet of Physics Orthogoal Fuctios Motivatio Aalogy with vectors You are probably failiar with the cocept of orthogoality fro vectors; two vectors are orthogoal whe they

More information

SOME TRIBONACCI IDENTITIES

SOME TRIBONACCI IDENTITIES Mathematics Today Vol.7(Dec-011) 1-9 ISSN 0976-38 Abstract: SOME TRIBONACCI IDENTITIES Shah Devbhadra V. Sir P.T.Sarvajaik College of Sciece, Athwalies, Surat 395001. e-mail : drdvshah@yahoo.com The sequece

More information

A Combinatoric Proof and Generalization of Ferguson s Formula for k-generalized Fibonacci Numbers

A Combinatoric Proof and Generalization of Ferguson s Formula for k-generalized Fibonacci Numbers Jue 5 00 A Combiatoric Proof ad Geeralizatio of Ferguso s Formula for k-geeralized Fiboacci Numbers David Kessler 1 ad Jeremy Schiff 1 Departmet of Physics Departmet of Mathematics Bar-Ila Uiversity, Ramat

More information

A New Type of q-szász-mirakjan Operators

A New Type of q-szász-mirakjan Operators Filoat 3:8 07, 567 568 https://doi.org/0.98/fil7867c Published by Faculty of Scieces ad Matheatics, Uiversity of Niš, Serbia Available at: http://www.pf.i.ac.rs/filoat A New Type of -Szász-Miraka Operators

More information

The heat equation. But how to evaluate: F F? T T. We first transform this analytical equation into the corresponding form in finite-space:

The heat equation. But how to evaluate: F F? T T. We first transform this analytical equation into the corresponding form in finite-space: L9_0/Cop. Astro.-, HS 0, Uiversität Basel/AAH The heat equatio T t T x. We first trasfor this aalytical equatio ito the correspodig for i fiite-space:.. t T δt T T F F F F F x x x x T T T F x But how to

More information

Problem. Consider the sequence a j for j N defined by the recurrence a j+1 = 2a j + j for j > 0

Problem. Consider the sequence a j for j N defined by the recurrence a j+1 = 2a j + j for j > 0 GENERATING FUNCTIONS Give a ifiite sequece a 0,a,a,, its ordiary geeratig fuctio is A : a Geeratig fuctios are ofte useful for fidig a closed forula for the eleets of a sequece, fidig a recurrece forula,

More information

arxiv: v1 [math.co] 23 Mar 2016

arxiv: v1 [math.co] 23 Mar 2016 The umber of direct-sum decompositios of a fiite vector space arxiv:603.0769v [math.co] 23 Mar 206 David Ellerma Uiversity of Califoria at Riverside August 3, 208 Abstract The theory of q-aalogs develops

More information

Some Explicit Formulae of NAF and its Left-to-Right. Analogue Based on Booth Encoding

Some Explicit Formulae of NAF and its Left-to-Right. Analogue Based on Booth Encoding Vol.7, No.6 (01, pp.69-74 http://dx.doi.org/10.1457/ijsia.01.7.6.7 Some Explicit Formulae of NAF ad its Left-to-Right Aalogue Based o Booth Ecodig Dog-Guk Ha, Okyeo Yi, ad Tsuyoshi Takagi Kookmi Uiversity,

More information

Analysis of Algorithms. Introduction. Contents

Analysis of Algorithms. Introduction. Contents Itroductio The focus of this module is mathematical aspects of algorithms. Our mai focus is aalysis of algorithms, which meas evaluatig efficiecy of algorithms by aalytical ad mathematical methods. We

More information

CSI 2101 Discrete Structures Winter Homework Assignment #4 (100 points, weight 5%) Due: Thursday, April 5, at 1:00pm (in lecture)

CSI 2101 Discrete Structures Winter Homework Assignment #4 (100 points, weight 5%) Due: Thursday, April 5, at 1:00pm (in lecture) CSI 101 Discrete Structures Witer 01 Prof. Lucia Moura Uiversity of Ottawa Homework Assigmet #4 (100 poits, weight %) Due: Thursday, April, at 1:00pm (i lecture) Program verificatio, Recurrece Relatios

More information

A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence

A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece,, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet as

More information

Linear chord diagrams with long chords

Linear chord diagrams with long chords Liear chord diagrams with log chords Everett Sulliva Departmet of Mathematics Dartmouth College Haover New Hampshire, U.S.A. everett..sulliva@dartmouth.edu Submitted: Feb 7, 2017; Accepted: Oct 7, 2017;

More information

SphericalHarmonicY. Notations. Primary definition. Specific values. Traditional name. Traditional notation. Mathematica StandardForm notation

SphericalHarmonicY. Notations. Primary definition. Specific values. Traditional name. Traditional notation. Mathematica StandardForm notation SphericalHaroicY Notatios Traditioal ae Spherical haroic Traditioal otatio Y, φ Matheatica StadardFor otatio SphericalHaroicY,,, φ Priary defiitio 05.0.0.000.0 Y, φ φ P cos ; 05.0.0.000.0 Y 0 k, φ k k

More information

The multiplicative structure of finite field and a construction of LRC

The multiplicative structure of finite field and a construction of LRC IERG6120 Codig for Distributed Storage Systems Lecture 8-06/10/2016 The multiplicative structure of fiite field ad a costructio of LRC Lecturer: Keeth Shum Scribe: Zhouyi Hu Notatios: We use the otatio

More information

Combinatorial and Automated Proofs of Certain Identities

Combinatorial and Automated Proofs of Certain Identities Cobiatorial ad Autoated Proofs of Certai Idetities The MIT Faculty has ade this article opely available Please share how this access beefits you Your story atters Citatio As Published Publisher Brereto,

More information

SZEGO S THEOREM STARTING FROM JENSEN S THEOREM

SZEGO S THEOREM STARTING FROM JENSEN S THEOREM UPB Sci Bull, Series A, Vol 7, No 3, 8 ISSN 3-77 SZEGO S THEOREM STARTING FROM JENSEN S THEOREM Cǎli Alexe MUREŞAN Mai îtâi vo itroduce Teorea lui Jese şi uele coseciţe ale sale petru deteriarea uǎrului

More information

arxiv: v1 [math.fa] 3 Apr 2016

arxiv: v1 [math.fa] 3 Apr 2016 Aticommutator Norm Formula for Proectio Operators arxiv:164.699v1 math.fa] 3 Apr 16 Sam Walters Uiversity of Norther British Columbia ABSTRACT. We prove that for ay two proectio operators f, g o Hilbert

More information

International Journal of Multidisciplinary Research and Modern Education (IJMRME) ISSN (Online): (www.rdmodernresearch.

International Journal of Multidisciplinary Research and Modern Education (IJMRME) ISSN (Online): (www.rdmodernresearch. (wwwrdoderresearchco) Volue II, Issue II, 2016 PRODUC OPERAION ON FUZZY RANSIION MARICES V Chiadurai*, S Barkavi**, S Vayabalaji*** & J Parthiba**** * Departet of Matheatics, Aaalai Uiversity, Aaalai Nagar,

More information

6 Integers Modulo n. integer k can be written as k = qn + r, with q,r, 0 r b. So any integer.

6 Integers Modulo n. integer k can be written as k = qn + r, with q,r, 0 r b. So any integer. 6 Itegers Modulo I Example 2.3(e), we have defied the cogruece of two itegers a,b with respect to a modulus. Let us recall that a b (mod ) meas a b. We have proved that cogruece is a equivalece relatio

More information

AN EFFICIENT ESTIMATION METHOD FOR THE PARETO DISTRIBUTION

AN EFFICIENT ESTIMATION METHOD FOR THE PARETO DISTRIBUTION Joural of Statistics: Advaces i Theory ad Applicatios Volue 3, Nuber, 00, Pages 6-78 AN EFFICIENT ESTIMATION METHOD FOR THE PARETO DISTRIBUTION Departet of Matheatics Brock Uiversity St. Catharies, Otario

More information

Math 475, Problem Set #12: Answers

Math 475, Problem Set #12: Answers Math 475, Problem Set #12: Aswers A. Chapter 8, problem 12, parts (b) ad (d). (b) S # (, 2) = 2 2, sice, from amog the 2 ways of puttig elemets ito 2 distiguishable boxes, exactly 2 of them result i oe

More information

CALCULATING FIBONACCI VECTORS

CALCULATING FIBONACCI VECTORS THE GENERALIZED BINET FORMULA FOR CALCULATING FIBONACCI VECTORS Stuart D Aderso Departmet of Physics, Ithaca College 953 Daby Road, Ithaca NY 14850, USA email: saderso@ithacaedu ad Dai Novak Departmet

More information

P. Z. Chinn Department of Mathematics, Humboldt State University, Arcata, CA

P. Z. Chinn Department of Mathematics, Humboldt State University, Arcata, CA RISES, LEVELS, DROPS AND + SIGNS IN COMPOSITIONS: EXTENSIONS OF A PAPER BY ALLADI AND HOGGATT S. Heubach Departmet of Mathematics, Califoria State Uiversity Los Ageles 55 State Uiversity Drive, Los Ageles,

More information

Chapter 6. Advanced Counting Techniques

Chapter 6. Advanced Counting Techniques Chapter 6 Advaced Coutig Techiques 6.: Recurrece Relatios Defiitio: A recurrece relatio for the sequece {a } is a equatio expressig a i terms of oe or more of the previous terms of the sequece: a,a2,a3,,a

More information

A Tabu Search Method for Finding Minimal Multi-Homogeneous Bézout Number

A Tabu Search Method for Finding Minimal Multi-Homogeneous Bézout Number Joural of Matheatics ad Statistics 6 (): 105-109, 010 ISSN 1549-3644 010 Sciece Publicatios A Tabu Search Method for Fidig Miial Multi-Hoogeeous Bézout Nuber Hassa M.S. Bawazir ad Ali Abd Raha Departet

More information

A Note on the Symmetric Powers of the Standard Representation of S n

A Note on the Symmetric Powers of the Standard Representation of S n A Note o the Symmetric Powers of the Stadard Represetatio of S David Savitt 1 Departmet of Mathematics, Harvard Uiversity Cambridge, MA 0138, USA dsavitt@mathharvardedu Richard P Staley Departmet of Mathematics,

More information

5.6 Binomial Multi-section Matching Transformer

5.6 Binomial Multi-section Matching Transformer 4/14/2010 5_6 Bioial Multisectio Matchig Trasforers 1/1 5.6 Bioial Multi-sectio Matchig Trasforer Readig Assiget: pp. 246-250 Oe way to axiize badwidth is to costruct a ultisectio Γ f that is axially flat.

More information

The Random Walk For Dummies

The Random Walk For Dummies The Radom Walk For Dummies Richard A Mote Abstract We look at the priciples goverig the oe-dimesioal discrete radom walk First we review five basic cocepts of probability theory The we cosider the Beroulli

More information

Expression for Restricted Partition Function through Bernoulli Polynomials

Expression for Restricted Partition Function through Bernoulli Polynomials Expressio for Restricted Partitio Fuctio through Beroulli Polyomials Boris Y. Rubistei Departmet of Mathematics, Uiversity of Califoria, Davis, Oe Shields Dr., Davis, CA 9566, U.S.A. February 28, 2005

More information

A New Recursion for Space-Filling Geometric Fractals

A New Recursion for Space-Filling Geometric Fractals A New Recursio for Space-Fillig Geometric Fractals Joh Shier Abstract. A recursive two-dimesioal geometric fractal costructio based upo area ad perimeter is described. For circles the radius of the ext

More information

Some remarks on the paper Some elementary inequalities of G. Bennett

Some remarks on the paper Some elementary inequalities of G. Bennett Soe rears o the paper Soe eleetary iequalities of G. Beett Dag Ah Tua ad Luu Quag Bay Vieta Natioal Uiversity - Haoi Uiversity of Sciece Abstract We give soe couterexaples ad soe rears of soe of the corollaries

More information