Algorithms for bioinformatics Part 2: Data structures

Size: px
Start display at page:

Download "Algorithms for bioinformatics Part 2: Data structures"

Transcription

1 Alorithms for bioinformtics Prt 2: Dt structures Greory Kucherov LIGM/CNRS Mrne-l-Vllée

2 Pln Clssicl indexes Suffix trees DAWG nd Position heps Suffix rrys Succinct (compressed) indexes Burrows-Wheeler trnsform nd FM-index

3 Suffix rry

4 Suffix rry: definition T=ctct ctct ctct tct tct ct ct t t t t

5 Suffix rry: definition T=ctct ct ctct t tct t ct ctct t tct t

6 Suffix rry: definition T=ctct ct ctct t tct t ct ctct t tct t SA Suffix rry (SA) is permuttion of the positions of T tken in the lexicorphicl order of correspondin suffixes T[SA[i]n]< lex T[SA[i+]n] [Mnber, Myers 90] [Gonnet, Bez-Ytes, Snider 92]

7 Strin mtchin with suffix rrys T=ctct ct ctct t tct t ct ctct t tct t SA c c For ny substrin, its occurrences form n intervl of SA Intervl [L P,R P ] for P cn be found by two binry serches Serch for L P : while l<r do m=(l+r)/2 if P< lex T[SA[m]n] then l m else r m the whole serch for P tkes time O(m lo(n)+occ)

8 Strin mtchin with suffix rrys T=ctct ct ctct t tct t ct ctct t tct t SA LCP LCP[i] = lenth of the lonest common prefix between suffixes strtin t SA[i] nd SA[i-] Note: the lonest common prefix between ny two suffixes SA[i] nd SA[j] cn be computed in constnt time fter liner-time preprocessin of LCP rry. (Use rne minimum queries) usin LCP rry, intervl serch time O(m lo(n)) cn be turned into O(m+lo(n))

9 Computin LCP from SA (nd T) compute Rnk[i] = SA - [i] compute LCP[i] for i=rnk[],,rnk[n] to compute LCP[Rnk[i]], look t LCP[Rnk[i-]] let q=rnk[i], p=rnk[i-] key observtion: if h=lcp[p]>0, then LCP[q] h- h- comprisons cn be sved! k i- i h j- j h bcbdbe be h- bcbdbe bebe be SA j- i- j k i p- p q- q j=sa[p-]+

10 Computin LCP from SA (nd T) compute Rnk[i] = SA - [i] compute LCP[i] for i=rnk[],,rnk[n] to compute LCP[Rnk[i]], look t LCP[Rnk[i-]] let q=rnk[i], p=rnk[i-] key observtion: if h=lcp[p]>0, then LCP[q] h- h- comprisons cn be sved! k i- i h j- j compre successively T[i+h], T[i+h+], with respectively T[k+h], T[k+h+], h bcbdbe be h- bcbdbe bebe be SA j- i- j k i p- p q- q j=sa[p-]+

11 Computin LCP from SA (nd T) compute Rnk[i] = SA - [i] compute LCP[i] for i=rnk[],,rnk[n] to compute LCP[Rnk[i]], look t LCP[Rnk[i-]] let q=rnk[i], p=rnk[i-] key observtion: if h=lcp[p]>0, then LCP[q] h- h- comprisons cn be sved! k i- i h j- j compre successively T[i+h], T[i+h+], with respectively T[k+h], T[k+h+], compute new h=lcp[q] iterte to i+ runnin time: O(n) h bcbdbe be h- bcbdbe bebe be SA j- i- j k i j=sa[p-]+ p- p q- q

12 Strin mtchin with suffix rrys P=ccc lcp l suf l = suf r = cccc SA l r ssume we know lcp l =LCP(suf l,p), lcp r =LCP(suf r,p) nd ssume lcp l <lcp r lcp r

13 Strin mtchin with suffix rrys P=ccc lcp l suf l = suf m = ccc suf r = cccc lcp r SA l m r ssume we know lcp l =LCP(suf l,p), lcp r =LCP(suf r,p) nd ssume lcp l <lcp r consider suf m nd consider LCP(suf m,suf r )

14 Strin mtchin with suffix rrys P=ccc lcp l suf l = suf m = ccc suf r = cccc lcp r SA l m r ssume we know lcp l =LCP(suf l,p), lcp r =LCP(suf r,p) nd ssume lcp l <lcp r consider suf m nd consider LCP(suf m,suf r ) Cse : LCP(suf m,suf r )<suf r Then P > lex suf m, LCP(suf m,p)=lcp(suf m,suf r ) nd we set l m for the next itertion

15 Strin mtchin with suffix rrys P=cc lcp l suf l = suf m = ccc suf r = cccc lcp r SA l m r ssume we know lcp l =LCP(suf l,p), lcp r =LCP(suf r,p) nd ssume lcp l <lcp r consider suf m nd consider LCP(suf m,suf r ) Cse 2: LCP(suf m,suf r )>suf r Then P < lex suf m, LCP(suf m,p)=lcp r nd we set r m for the next itertion

16 Strin mtchin with suffix rrys P=ccc lcp l suf l = suf m = ccc suf r = cccc lcp r SA l m r ssume we know lcp l =LCP(suf l,p), lcp r =LCP(suf r,p) nd ssume lcp l <lcp r consider suf m nd consider LCP(suf m,suf r ) Cse 3: LCP(suf m,suf r )=suf r Then we keep comprin chrs of P with those of suf m until P[j] suf m [j] which determines if P lex suf m or P> lex suf m, nd lso the vlue LCP(suf m,p)

17 Strin mtchin with suffix rrys P=ccc lcp l suf l = suf m = ccc suf r = cccc lcp r SA l m r ssume we know lcp l =LCP(suf l,p), lcp r =LCP(suf r,p) nd ssume lcp l <lcp r consider suf m nd consider LCP(suf m,suf r ) Cse 3: LCP(suf m,suf r )=suf r Then we keep comprin chrs of P with those of suf m until P[j] suf m [j] which determines if P lex suf m or P> lex suf m, nd lso the vlue LCP(suf m,p)

18 Strin mtchin with suffix rrys P=ccc suf l = suf m = ccc lcp m suf r = cccc lcp r SA l m r ssume we know lcp l =LCP(suf l,p), lcp r =LCP(suf r,p) nd ssume lcp l <lcp r consider suf m nd consider LCP(suf m,suf r ) Cse 3: LCP(suf m,suf r )=suf r Then we keep comprin chrs of P with those of suf m until P[j] suf m [j] which determines if P lex suf m or P> lex suf m, nd lso the vlue lcp m =LCP(suf m,p)

19 Strin mtchin with suffix rrys P=ccc suf l = suf m = ccc lcp m suf r = cccc lcp r SA l m r ssume we know lcp l =LCP(suf l,p), lcp r =LCP(suf r,p) nd ssume lcp l <lcp r consider suf m nd consider LCP(suf m,suf r ) Cse 3: LCP(suf m,suf r )=suf r Then we keep comprin chrs of P with those of suf m until P[j] suf m [j] which determines if P lex suf m or P> lex suf m, nd lso the vlue lcp m =LCP(suf m,p) Conclusion: t ech step we either do binry division, or move one chr forwrd in the pttern time O( P +lo(n))

20 Construction of suffix rry construction in O(n lo(n)) [Mnber, Myers 90] construction in O(n) [Kärkkäinen, Snders 03] [Ko, Aluru 03] [Kim et l 03] works on prcticl liner-time construction of suffix rry: [Non et l 09]

21 Suffix rry nd suffix tree T=ctct ct ctct t tct t ct ctct t tct t SA LCP t c t c t 7 t t c t 3 9 c t 6 t c t 2 t t 8 c t SA lef lbels left-to-riht provided tht children of ny node re ordered ccordin to the order of chrcters LCP strin lenth of lowest common ncestor (lc) 4 0

22 Enhnced suffix rrys 2 c t c t 5 t 7 t t 9 3 t c c 6 t t 2 c t t 8 t 0 4 t c T=ctct LCP LCP rry encodes the tree topoloy of ST: internl node intervl LCP[ij] s.t. min{lcp[i+j]}=q nd LCP[i]<q, LCP[j+]<q on LCP nd SA rrys, one cn simulte bottom-up nd top-down nvition in ST, s well s suffix links [Abouelhod et l 04]

23 Enhnced suffix rrys 2 c t c t 5 t 7 t t 9 3 t c c 6 t t 2 c t t 8 t 0 4 t c T=ctct LCP LCP rry encodes the tree topoloy of ST: internl node intervl LCP[ij] s.t. min{lcp[i+j]}=q nd LCP[i]<q, LCP[j+]<q on LCP nd SA rrys, one cn simulte bottom-up nd top-down nvition in ST, s well s suffix links [Abouelhod et l 04]

24 Suffix rrys: prcticl issues SA is populr dt structure tht provides spceefficient lterntive to ST Enhnced suffix rrys tke bits/chr, this results to 5-27Gb for humn enome Enhnced SAs re used in prcticl bioinformtics softwre (Vmtch,seemehl) Efficient SA construction is n ctive re of reserch (externl memory, prlleliztion, ) Compressed suffix rrys [Grossi,Vitter 00] (different structures surveyed in [Nvrro, Mkinen 07])

25 Burrows-Wheeler trnsform nd FM-index

26 Succinct nd compressed indexes succinct index tkes spce in bits proportionl to tht of the text itself previous indexes re not succinct s they tke O(n) computer words but O(n lo(n)) bits compressed index tkes spce in bits proportionl to tht of the compressed text self-index does not require storin the text

27 Burrows-Wheeler trnsform T=ctct ctct ctct ctct tctc tctc tctc ctct ctct ctct tctc tctc tctc

28 Burrows-Wheeler trnsform T=ctct ctct ctct ctct tctc tctc tctc ctct ctct ctct tctc tctc tctc SA

29 Burrows-Wheeler trnsform T=ctct ctct ctct ctct tctc tctc tctc ctct ctct ctct tctc tctc tctc

30 Burrows-Wheeler trnsform T=ctct T[SA[i]] BWT BWT[i]=T[SA[i]-] if SA[i], otherwise ctct ctc t BWT hs been defined for the purpose of ctct compression, s BWT compresses better tct c thn the input text tct c BWT is reversible! tctc c tct c tct ctc t tctc t ctc t ctc

31 Burrows-Wheeler trnsform T=ctct T[SA[i]] BWT BWT[i]=T[SA[i]-] if SA[i], otherwise ctct ctc t Obs : the first column (F) is esy to ctct reconstruct, it cn be represented by n tct c rry C[x]= y<x occ(y,t) for ech letter x tct c Ex: C=[0,,6,8,0] tctc c tct c tct ctc t tctc t ctc t ctc F L

32 Burrows-Wheeler trnsform T= T[SA[i]] BWT BWT[i]=T[SA[i]-] if SA[i], otherwise ctct ctc t Obs : the first column (F) is esy to ctct reconstruct, it cn be represented by n tct c rry C[x]= y<x occ(y,t) for ech letter x tct c Ex: C=[0,,6,8,0] tctc c tct c tct ctc t tctc t ctc t ctc F L

33 Burrows-Wheeler trnsform T= T[SA[i]] BWT BWT[i]=T[SA[i]-] if SA[i], otherwise ctct ctc t Obs : the first column (F) is esy to ctct reconstruct, it cn be represented by n tct c rry C[x]= y<x occ(y,t) for ech letter x tct c Ex: C=[0,,6,8,0] tctc c tct c tct ctc t tctc t ctc t ctc F L

34 Burrows-Wheeler trnsform T= T[SA[i]] BWT BWT[i]=T[SA[i]-] if SA[i], otherwise ctct ctc t Obs : the first column (F) is esy to ctct reconstruct, it cn be represented by n tct c rry C[x]= y<x occ(y,t) for ech letter x tct c Ex: C=[0,,6,8,0] tctc c tct c tct ctc t tctc t ctc t ctc F L

35 Burrows-Wheeler trnsform T= T[SA[i]] BWT BWT[i]=T[SA[i]-] if SA[i], otherwise ctct ctc t Obs : the first column (F) is esy to ctct reconstruct, it cn be represented by n tct c rry C[x]= y<x occ(y,t) for ech letter x tct c Ex: C=[0,,6,8,0] tctc c tct Obs 2: for identicl chrs, their reltive c tct order in F nd L is the sme j ctc t x tctc x i t ctc t ctc LF[i]=C[BWT[i]]+rnk[BWT[i],i] F L Ex: LF[]=8+

36 Burrows-Wheeler trnsform T= t T[SA[i]] BWT BWT[i]=T[SA[i]-] if SA[i], otherwise ctct ctc t Obs : the first column (F) is esy to ctct reconstruct, it cn be represented by n tct c rry C[x]= y<x occ(y,t) for ech letter x tct c Ex: C=[0,,6,8,0] tctc c tct Obs 2: for identicl chrs, their reltive c tct order in F nd L is the sme j ctc t x tctc x i t ctc t ctc LF[i]=C[BWT[i]]+rnk[BWT[i],i] F L Ex: LF[]=8+

37 Burrows-Wheeler trnsform T= t T[SA[i]] BWT BWT[i]=T[SA[i]-] if SA[i], otherwise ctct ctc t Obs : the first column (F) is esy to ctct reconstruct, it cn be represented by n tct c rry C[x]= y<x occ(y,t) for ech letter x tct c Ex: C=[0,,6,8,0] tctc c tct Obs 2: for identicl chrs, their reltive c tct order in F nd L is the sme j ctc t x tctc x i t ctc t ctc LF[i]=C[BWT[i]]+rnk[BWT[i],i] F L Ex: LF[]=8+

38 Burrows-Wheeler trnsform T= t T[SA[i]] BWT BWT[i]=T[SA[i]-] if SA[i], otherwise ctct ctc t Obs : the first column (F) is esy to ctct reconstruct, it cn be represented by n tct c rry C[x]= y<x occ(y,t) for ech letter x tct c Ex: C=[0,,6,8,0] tctc c tct Obs 2: for identicl chrs, their reltive c tct order in F nd L is the sme j ctc t x tctc x i t ctc t ctc LF[i]=C[BWT[i]]+rnk[BWT[i],i] F L Ex: LF[]=8+

39 Burrows-Wheeler trnsform T= t T[SA[i]] BWT BWT[i]=T[SA[i]-] if SA[i], otherwise ctct ctc t Obs : the first column (F) is esy to ctct reconstruct, it cn be represented by n tct c rry C[x]= y<x occ(y,t) for ech letter x tct c Ex: C=[0,,6,8,0] tctc c tct Obs 2: for identicl chrs, their reltive c tct order in F nd L is the sme j ctc t x tctc x i t ctc t ctc LF[i]=C[BWT[i]]+rnk[BWT[i],i] F L Ex: LF[]=8+

40 Burrows-Wheeler trnsform T= t T[SA[i]] BWT BWT[i]=T[SA[i]-] if SA[i], otherwise ctct ctc t Obs : the first column (F) is esy to ctct reconstruct, it cn be represented by n tct c rry C[x]= y<x occ(y,t) for ech letter x tct c Ex: C=[0,,6,8,0] tctc c tct Obs 2: for identicl chrs, their reltive c tct order in F nd L is the sme j ctc t x tctc x i t ctc t ctc LF[i]=C[BWT[i]]+rnk[BWT[i],i] F L Ex: LF[]=8+

41 Burrows-Wheeler trnsform T= t T[SA[i]] BWT BWT[i]=T[SA[i]-] if SA[i], otherwise ctct ctc t Obs : the first column (F) is esy to ctct reconstruct, it cn be represented by n tct c rry C[x]= y<x occ(y,t) for ech letter x tct c Ex: C=[0,,6,8,0] tctc c tct Obs 2: for identicl chrs, their reltive c tct order in F nd L is the sme j ctc t x tctc x i t ctc t ctc LF[i]=C[BWT[i]]+rnk[BWT[i],i] F L Ex: LF[]=8+

42 LF function T=ctct T[SA[i]] BWT ctct ctc t ctct tct c tct c tctc c tct c tct ctc t tctc t ctc t ctc LF[i]=C[BWT[i]]+rnk[BWT[i],i] LF[i] yields the index (in SA) of the suffix immeditely precedin (in T) the i-th suffix (in SA). Formlly, LF[SA[i]]=SA[i-] F L

43 rnk function T[SA[i]] BWT ctct ctc t ctct tct c tct c tctc c tct c tct ctc t tctc t ctc t ctc LF[i]=C[BWT[i]]+rnk[BWT[i],i] F L

44 rnk function T[SA[i]] BWT ctct ctc t ctct tct c tct c tctc c tct c tct ctc t tctc t ctc t ctc F L rnk[bwt[i],i] LF[i]=C[BWT[i]]+rnk[BWT[i],i]

45 rnk function T[SA[i]] BWT ctct ctc t ctct tct c tct c tctc c tct c tct ctc t tctc t ctc t ctc F L rnk[bwt[i],i] LF[i]=C[BWT[i]]+rnk[BWT[i],i] how bout enerl queries rnk[,i] for ny letter nd ny position i?

46 rnk/select functions iven strin T, efficiently nswer queries rnk(,i) on the number of s in T[i] rnk function (on bit vectors) turns out to be fundmentl lorithmic block for buildin succinct dt structures [Jcobson 89] rnk cn be supported in time O() usin o(n) dditionl bits of memory complementry function select(,j): output the position of the j-th occurrence of in T. select cn lso be supported in O() time

47 Strin mtchin with BWT T=ctct T[SA[i]] BWT ctct ctc t ctct tct c tct c tctc c tct c tct ctc t tctc t ctc t ctc P=tc F L

48 Strin mtchin with BWT T=ctct T[SA[i]] BWT ctct ctc t ctct tct c tct c tctc c tct c tct ctc t tctc t ctc t ctc P=tc F L

49 Strin mtchin with BWT T=ctct T[SA[i]] BWT ctct ctc t ctct tct c tct c tctc c tct c tct ctc t tctc t ctc t ctc P=tc F L

50 Strin mtchin with BWT T=ctct T[SA[i]] BWT ctct ctc t ctct tct c tct c tctc c tct c tct ctc t tctc t ctc t ctc e f P=tc [e,f] : current intervl x : letter e:= C[x]+rnk[x,e]+ f:= C[x]+rnk[x,f] F L

51 Strin mtchin with BWT T=ctct T[SA[i]] BWT ctct ctc t ctct tct c tct c tctc c tct c tct ctc t tctc t ctc t ctc e f P=tc [e,f] : current intervl x : letter e:= C[x]+rnk[x,e]+ f:= C[x]+rnk[x,f] F L

52 Strin mtchin with BWT T=ctct T[SA[i]] BWT ctct ctc t ctct tct c tct c tctc c tct c tct ctc t tctc t ctc t ctc e f P=tc [e,f] : current intervl x : letter e:= C[x]+rnk[x,e]+ f:= C[x]+rnk[x,f] F L

53 Strin mtchin with BWT T=ctct T[SA[i]] BWT ctct ctc t ctct tct c tct c tctc c tct c tct ctc t tctc t ctc t ctc e f P=tc [e,f] : current intervl x : letter e:= C[x]+rnk[x,e]+ f:= C[x]+rnk[x,f] F L

54 Strin mtchin with BWT T=ctct T[SA[i]] BWT ctct ctc t ctct tct c tct c tctc c tct c tct ctc t tctc t ctc t ctc e f P=tc [e,f] : current intervl x : letter e:= C[x]+rnk[x,e]+ f:= C[x]+rnk[x,f] F L

55 Strin mtchin with BWT T=ctct T[SA[i]] BWT ctct ctc t ctct tct c tct c tctc c tct c tct ctc t tctc t ctc t ctc e f P=tc [e,f] : current intervl x : letter e:= C[x]+rnk[x,e]+ f:= C[x]+rnk[x,f] F L

56 Strin mtchin with BWT T=ctct T[SA[i]] BWT ctct ctc t ctct tct c tct c tctc c tct c tct ctc t tctc t ctc t ctc P=tc [e,f] : current intervl x : letter e:= C[x]+rnk[x,e]+ f:= C[x]+rnk[x,f] F L

57 Strin mtchin with BWT T=ctct T[SA[i]] BWT ctct ctc t ctct tct c tct c tctc c tct c tct ctc t tctc t ctc t ctc P=tc Wht position is it?? F L

58 Strin mtchin with BWT T=ctct T[SA[i]] BWT ctct ctc t ctct tct c tct c tctc c tct c tct ctc t tctc t ctc t ctc P=tc It is position 4! F L

59 FM-index T=ctct SA T[SA[i]] BWT ctct ctc t ctct tct c tct c tctc c tct c tct ctc t tctc t ctc t ctc F L Solution: store only frction of vlues of SA. Storin one vlue over lo(n) leds to O(n lo(n)/lo(n))=o(n) bits

60 FM-index T=ctct SA T[SA[i]] BWT ctct ctc t ctct tct c tct c tctc c tct c tct ctc t tctc t ctc t ctc F L Solution: store only frction of vlues of SA. Storin one vlue over lo(n) leds to O(n lo(n)/lo(n))=o(n) bits Serch time becomes O( P lo(n)) [Ferrin, Mnzini 00] FM-index includes: BWT selection of SA vlues uxiliry structures: rry C, rnk, position mrkin

61 FM-index: prcticl issues FM-index cn be implemented usin ~3 bits/chr (!!) FM-index is now widely used in prcticl bioinformtics softwre: BWA, bowtie, SOAP2 (mppin), CGA (ssembly) other succinct dt structure exist (includin compct suffix rry) nd continue to pper construction my require much more spce thn the resultin structure externl memory lorithms re importnt lorithms specilized to multi-core or GPU processor rchitectures dynmic indexes

Algorithms in Computational. Biology. More on BWT

Algorithms in Computational. Biology. More on BWT Algorithms in Computtionl Biology More on BWT tody Plese Lst clss! don't forget to submit And by next (vi emil, repo ) implementtion week or shre prgectfltw get Not I would like reding overview! Discuss

More information

Module 9: Tries and String Matching

Module 9: Tries and String Matching Module 9: Tries nd String Mtching CS 240 - Dt Structures nd Dt Mngement Sjed Hque Veronik Irvine Tylor Smith Bsed on lecture notes by mny previous cs240 instructors Dvid R. Cheriton School of Computer

More information

Module 9: Tries and String Matching

Module 9: Tries and String Matching Module 9: Tries nd String Mtching CS 240 - Dt Structures nd Dt Mngement Sjed Hque Veronik Irvine Tylor Smith Bsed on lecture notes by mny previous cs240 instructors Dvid R. Cheriton School of Computer

More information

1 APL13: Suffix Arrays: more space reduction

1 APL13: Suffix Arrays: more space reduction 1 APL13: Suffix Arrys: more spce reduction In Section??, we sw tht when lphbet size is included in the time nd spce bounds, the suffix tree for string of length m either requires Θ(m Σ ) spce or the minimum

More information

Fingerprint idea. Assume:

Fingerprint idea. Assume: Fingerprint ide Assume: We cn compute fingerprint f(p) of P in O(m) time. If f(p) f(t[s.. s+m 1]), then P T[s.. s+m 1] We cn compre fingerprints in O(1) We cn compute f = f(t[s+1.. s+m]) from f(t[s.. s+m

More information

On Suffix Tree Breadth

On Suffix Tree Breadth On Suffix Tree Bredth Golnz Bdkoeh 1,, Juh Kärkkäinen 2, Simon J. Puglisi 2,, nd Bell Zhukov 2, 1 Deprtment of Computer Science University of Wrwick Conventry, United Kingdom g.dkoeh@wrwick.c.uk 2 Helsinki

More information

Alignment of Long Sequences. BMI/CS Spring 2016 Anthony Gitter

Alignment of Long Sequences. BMI/CS Spring 2016 Anthony Gitter Alignment of Long Sequences BMI/CS 776 www.biostt.wisc.edu/bmi776/ Spring 2016 Anthony Gitter gitter@biostt.wisc.edu Gols for Lecture Key concepts how lrge-scle lignment differs from the simple cse the

More information

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies Stte spce systems nlysis (continued) Stbility A. Definitions A system is sid to be Asymptoticlly Stble (AS) when it stisfies ut () = 0, t > 0 lim xt () 0. t A system is AS if nd only if the impulse response

More information

Applied Databases. Sebastian Maneth. Lecture 16 Suffix Array, Burrows-Wheeler Transform. University of Edinburgh - March 16th, 2017

Applied Databases. Sebastian Maneth. Lecture 16 Suffix Array, Burrows-Wheeler Transform. University of Edinburgh - March 16th, 2017 Applied Dtbses Lecture 16 Suffix Arry, Burrows-Wheeler Trsform Sebsti Meth Uiversity of Ediburgh - Mrch 16th, 2017 Outlie 2 1. Suffix Arry 2. Burrows-Wheeler Trsform Outlie 3 1. Suffix Arry 2. Burrows-Wheeler

More information

Applied Databases. Sebastian Maneth. Lecture 16 Suffix Array, Burrows-Wheeler Transform. University of Edinburgh - March 10th, 2016

Applied Databases. Sebastian Maneth. Lecture 16 Suffix Array, Burrows-Wheeler Transform. University of Edinburgh - March 10th, 2016 Applied Dtbses Lecture 16 Suffix Arry, Burrows-Wheeler Trsform Sebsti Meth Uiversity of Ediburgh - Mrch 10th, 2016 2 Outlie 1. Suffix Arry 2. Burrows-Wheeler Trsform 3 Olie Strig-Mtchig how c we do Horspool

More information

Point Processing of Images. Point Processing of Images. EECE/CS 253 Image Processing. Point Processing

Point Processing of Images. Point Processing of Images. EECE/CS 253 Image Processing. Point Processing EECE/CS 253 me Processin Lecture Notes: Lecture Notes: The Point Processin of mes Richrd Aln Peters Deprtment of Electricl Enineerin nd Computer Science Fll Semester 2007 Point Processin of mes n diitl

More information

New data structures to reduce data size and search time

New data structures to reduce data size and search time New dt structures to reduce dt size nd serch time Tsuneo Kuwbr Deprtment of Informtion Sciences, Fculty of Science, Kngw University, Hirtsuk-shi, Jpn FIT2018 1D-1, No2, pp1-4 Copyright (c)2018 by The Institute

More information

Before we can begin Ch. 3 on Radicals, we need to be familiar with perfect squares, cubes, etc. Try and do as many as you can without a calculator!!!

Before we can begin Ch. 3 on Radicals, we need to be familiar with perfect squares, cubes, etc. Try and do as many as you can without a calculator!!! Nme: Algebr II Honors Pre-Chpter Homework Before we cn begin Ch on Rdicls, we need to be fmilir with perfect squres, cubes, etc Try nd do s mny s you cn without clcultor!!! n The nth root of n n Be ble

More information

Improper Integrals. Type I Improper Integrals How do we evaluate an integral such as

Improper Integrals. Type I Improper Integrals How do we evaluate an integral such as Improper Integrls Two different types of integrls cn qulify s improper. The first type of improper integrl (which we will refer to s Type I) involves evluting n integrl over n infinite region. In the grph

More information

Algorithms in Bioinformatics: A Practical Introduction. Suffix tree

Algorithms in Bioinformatics: A Practical Introduction. Suffix tree Alorithms in Bioinformtis: A Prtil Introdution Suffix tree Overview Wht is suffix tree? Simple pplition of suffix tree Liner time lorithm for onstrutin suffix tree Suffix rry FM-index 1-mismth serh Suffix

More information

2. Udi Manber, Gene Myers: Suffix arrays: A new method for online string searching, SIAM Journal on Computing 22:935-48,1993

2. Udi Manber, Gene Myers: Suffix arrays: A new method for online string searching, SIAM Journal on Computing 22:935-48,1993 7 Suffix rrys This exposition ws developed y Clemens Gröpl nd Knut Reinert. It is sed on the following sources, which re ll recommended reding: 1. Simon J. Puglisi, W. F. Smyth, nd Andrew Turpin, A txonomy

More information

20 MATHEMATICS POLYNOMIALS

20 MATHEMATICS POLYNOMIALS 0 MATHEMATICS POLYNOMIALS.1 Introduction In Clss IX, you hve studied polynomils in one vrible nd their degrees. Recll tht if p(x) is polynomil in x, the highest power of x in p(x) is clled the degree of

More information

DATA Search I 魏忠钰. 复旦大学大数据学院 School of Data Science, Fudan University. March 7 th, 2018

DATA Search I 魏忠钰. 复旦大学大数据学院 School of Data Science, Fudan University. March 7 th, 2018 DATA620006 魏忠钰 Serch I Mrch 7 th, 2018 Outline Serch Problems Uninformed Serch Depth-First Serch Bredth-First Serch Uniform-Cost Serch Rel world tsk - Pc-mn Serch problems A serch problem consists of:

More information

Where did dynamic programming come from?

Where did dynamic programming come from? Where did dynmic progrmming come from? String lgorithms Dvid Kuchk cs302 Spring 2012 Richrd ellmn On the irth of Dynmic Progrmming Sturt Dreyfus http://www.eng.tu.c.il/~mi/cd/ or50/1526-5463-2002-50-01-0048.pdf

More information

19 Optimal behavior: Game theory

19 Optimal behavior: Game theory Intro. to Artificil Intelligence: Dle Schuurmns, Relu Ptrscu 1 19 Optiml behvior: Gme theory Adversril stte dynmics hve to ccount for worst cse Compute policy π : S A tht mximizes minimum rewrd Let S (,

More information

Tries and suffixes trees

Tries and suffixes trees Trie: A dt-structure for set of words Tries nd suffixes trees Alon Efrt Comuter Science Dertment University of Arizon All words over the lhet Σ={,,..z}. In the slides, let sy tht the lhet is only {,,c,d}

More information

f(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral

f(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral Improper Integrls Every time tht we hve evluted definite integrl such s f(x) dx, we hve mde two implicit ssumptions bout the integrl:. The intervl [, b] is finite, nd. f(x) is continuous on [, b]. If one

More information

p-adic Egyptian Fractions

p-adic Egyptian Fractions p-adic Egyptin Frctions Contents 1 Introduction 1 2 Trditionl Egyptin Frctions nd Greedy Algorithm 2 3 Set-up 3 4 p-greedy Algorithm 5 5 p-egyptin Trditionl 10 6 Conclusion 1 Introduction An Egyptin frction

More information

Math& 152 Section Integration by Parts

Math& 152 Section Integration by Parts Mth& 5 Section 7. - Integrtion by Prts Integrtion by prts is rule tht trnsforms the integrl of the product of two functions into other (idelly simpler) integrls. Recll from Clculus I tht given two differentible

More information

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus Unit #9 : Definite Integrl Properties; Fundmentl Theorem of Clculus Gols: Identify properties of definite integrls Define odd nd even functions, nd reltionship to integrl vlues Introduce the Fundmentl

More information

Theoretical foundations of Gaussian quadrature

Theoretical foundations of Gaussian quadrature Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of

More information

Anatomy of a Deterministic Finite Automaton. Deterministic Finite Automata. A machine so simple that you can understand it in less than one minute

Anatomy of a Deterministic Finite Automaton. Deterministic Finite Automata. A machine so simple that you can understand it in less than one minute Victor Admchik Dnny Sletor Gret Theoreticl Ides In Computer Science CS 5-25 Spring 2 Lecture 2 Mr 3, 2 Crnegie Mellon University Deterministic Finite Automt Finite Automt A mchine so simple tht you cn

More information

Lecture 3. Introduction digital logic. Notes. Notes. Notes. Representations. February Bern University of Applied Sciences.

Lecture 3. Introduction digital logic. Notes. Notes. Notes. Representations. February Bern University of Applied Sciences. Lecture 3 Ferury 6 ern University of pplied ciences ev. f57fc 3. We hve seen tht circuit cn hve multiple (n) inputs, e.g.,, C, We hve lso seen tht circuit cn hve multiple (m) outputs, e.g. X, Y,, ; or

More information

1 Probability Density Functions

1 Probability Density Functions Lis Yn CS 9 Continuous Distributions Lecture Notes #9 July 6, 28 Bsed on chpter by Chris Piech So fr, ll rndom vribles we hve seen hve been discrete. In ll the cses we hve seen in CS 9, this ment tht our

More information

STRAND B: NUMBER THEORY

STRAND B: NUMBER THEORY Mthemtics SKE, Strnd B UNIT B Indices nd Fctors: Tet STRAND B: NUMBER THEORY B Indices nd Fctors Tet Contents Section B. Squres, Cubes, Squre Roots nd Cube Roots B. Inde Nottion B. Fctors B. Prime Fctors,

More information

The Regulated and Riemann Integrals

The Regulated and Riemann Integrals Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue

More information

Balanced binary search trees

Balanced binary search trees 02110 Inge Li Gørtz Overview Blnced binry serch trees: Red-blck trees nd 2-3-4 trees Amortized nlysis Dynmic progrmming Network flows String mtching String indexing Computtionl geometry Introduction to

More information

ODE: Existence and Uniqueness of a Solution

ODE: Existence and Uniqueness of a Solution Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =

More information

f(a+h) f(a) x a h 0. This is the rate at which

f(a+h) f(a) x a h 0. This is the rate at which M408S Concept Inventory smple nswers These questions re open-ended, nd re intended to cover the min topics tht we lerned in M408S. These re not crnk-out-n-nswer problems! (There re plenty of those in the

More information

Conservation Law. Chapter Goal. 5.2 Theory

Conservation Law. Chapter Goal. 5.2 Theory Chpter 5 Conservtion Lw 5.1 Gol Our long term gol is to understnd how mny mthemticl models re derived. We study how certin quntity chnges with time in given region (sptil domin). We first derive the very

More information

September 13 Homework Solutions

September 13 Homework Solutions College of Engineering nd Computer Science Mechnicl Engineering Deprtment Mechnicl Engineering 5A Seminr in Engineering Anlysis Fll Ticket: 5966 Instructor: Lrry Cretto Septemer Homework Solutions. Are

More information

7.2 The Definite Integral

7.2 The Definite Integral 7.2 The Definite Integrl the definite integrl In the previous section, it ws found tht if function f is continuous nd nonnegtive, then the re under the grph of f on [, b] is given by F (b) F (), where

More information

Preview 11/1/2017. Greedy Algorithms. Coin Change. Coin Change. Coin Change. Coin Change. Greedy algorithms. Greedy Algorithms

Preview 11/1/2017. Greedy Algorithms. Coin Change. Coin Change. Coin Change. Coin Change. Greedy algorithms. Greedy Algorithms Preview Greed Algorithms Greed Algorithms Coin Chnge Huffmn Code Greed lgorithms end to e simple nd strightforwrd. Are often used to solve optimiztion prolems. Alws mke the choice tht looks est t the moment,

More information

APPROXIMATE INTEGRATION

APPROXIMATE INTEGRATION APPROXIMATE INTEGRATION. Introduction We hve seen tht there re functions whose nti-derivtives cnnot be expressed in closed form. For these resons ny definite integrl involving these integrnds cnnot be

More information

Succinct Text Indexes on Large Alphabet

Succinct Text Indexes on Large Alphabet Succinct Text Indexes on Lrge Alphet Meng Zhng 1, Jijun Tng 2, Dong Guo 1, Ling Hu 1, nd Qing Li 1 1 College of Computer Science nd Technology, Jilin University, Chngchun 130012, Chin zm@mil.edu.cn, {guodg,

More information

Math 131. Numerical Integration Larson Section 4.6

Math 131. Numerical Integration Larson Section 4.6 Mth. Numericl Integrtion Lrson Section. This section looks t couple of methods for pproimting definite integrls numericlly. The gol is to get good pproimtion of the definite integrl in problems where n

More information

and that at t = 0 the object is at position 5. Find the position of the object at t = 2.

and that at t = 0 the object is at position 5. Find the position of the object at t = 2. 7.2 The Fundmentl Theorem of Clculus 49 re mny, mny problems tht pper much different on the surfce but tht turn out to be the sme s these problems, in the sense tht when we try to pproimte solutions we

More information

The problems that follow illustrate the methods covered in class. They are typical of the types of problems that will be on the tests.

The problems that follow illustrate the methods covered in class. They are typical of the types of problems that will be on the tests. ADVANCED CALCULUS PRACTICE PROBLEMS JAMES KEESLING The problems tht follow illustrte the methods covered in clss. They re typicl of the types of problems tht will be on the tests. 1. Riemnn Integrtion

More information

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite

More information

CS 188: Artificial Intelligence Spring 2007

CS 188: Artificial Intelligence Spring 2007 CS 188: Artificil Intelligence Spring 2007 Lecture 3: Queue-Bsed Serch 1/23/2007 Srini Nrynn UC Berkeley Mny slides over the course dpted from Dn Klein, Sturt Russell or Andrew Moore Announcements Assignment

More information

Problem Set 4: Solutions Math 201A: Fall 2016

Problem Set 4: Solutions Math 201A: Fall 2016 Problem Set 4: s Mth 20A: Fll 206 Problem. Let f : X Y be one-to-one, onto mp between metric spces X, Y. () If f is continuous nd X is compct, prove tht f is homeomorphism. Does this result remin true

More information

New Expansion and Infinite Series

New Expansion and Infinite Series Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University

More information

n f(x i ) x. i=1 In section 4.2, we defined the definite integral of f from x = a to x = b as n f(x i ) x; f(x) dx = lim i=1

n f(x i ) x. i=1 In section 4.2, we defined the definite integral of f from x = a to x = b as n f(x i ) x; f(x) dx = lim i=1 The Fundmentl Theorem of Clculus As we continue to study the re problem, let s think bck to wht we know bout computing res of regions enclosed by curves. If we wnt to find the re of the region below the

More information

Lecture 6. Notes. Notes. Notes. Representations Z A B and A B R. BTE Electronics Fundamentals August Bern University of Applied Sciences

Lecture 6. Notes. Notes. Notes. Representations Z A B and A B R. BTE Electronics Fundamentals August Bern University of Applied Sciences Lecture 6 epresenttions epresenttions TE52 - Electronics Fundmentls ugust 24 ern University of pplied ciences ev. c2d5c88 6. Integers () sign-nd-mgnitude representtion The set of integers contins the Nturl

More information

A recursive construction of efficiently decodable list-disjunct matrices

A recursive construction of efficiently decodable list-disjunct matrices CSE 709: Compressed Sensing nd Group Testing. Prt I Lecturers: Hung Q. Ngo nd Atri Rudr SUNY t Bufflo, Fll 2011 Lst updte: October 13, 2011 A recursive construction of efficiently decodble list-disjunct

More information

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn

More information

Math 61CM - Solutions to homework 9

Math 61CM - Solutions to homework 9 Mth 61CM - Solutions to homework 9 Cédric De Groote November 30 th, 2018 Problem 1: Recll tht the left limit of function f t point c is defined s follows: lim f(x) = l x c if for ny > 0 there exists δ

More information

1.3 Regular Expressions

1.3 Regular Expressions 56 1.3 Regulr xpressions These hve n importnt role in describing ptterns in serching for strings in mny pplictions (e.g. wk, grep, Perl,...) All regulr expressions of lphbet re 1.Ønd re regulr expressions,

More information

Energy Bands Energy Bands and Band Gap. Phys463.nb Phenomenon

Energy Bands Energy Bands and Band Gap. Phys463.nb Phenomenon Phys463.nb 49 7 Energy Bnds Ref: textbook, Chpter 7 Q: Why re there insultors nd conductors? Q: Wht will hppen when n electron moves in crystl? In the previous chpter, we discussed free electron gses,

More information

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a). The Fundmentl Theorems of Clculus Mth 4, Section 0, Spring 009 We now know enough bout definite integrls to give precise formultions of the Fundmentl Theorems of Clculus. We will lso look t some bsic emples

More information

Efficient Parallel Implementation of Classical Gram-Schmidt Orthogonalization Using Matrix Multiplication

Efficient Parallel Implementation of Classical Gram-Schmidt Orthogonalization Using Matrix Multiplication Efficient Prllel Implementtion of Clssicl Grm-Schmidt Orthogonliztion Using Mtrix Multipliction Tkuy Yokozw Tisuke Boku Disuke Tkhshi Mitsuhis Sto Grdute School of Systems nd Informtion Engineering, University

More information

Calculus of Variations

Calculus of Variations Clculus of Vritions Com S 477/577 Notes) Yn-Bin Ji Dec 4, 2017 1 Introduction A functionl ssigns rel number to ech function or curve) in some clss. One might sy tht functionl is function of nother function

More information

Riemann Sums and Riemann Integrals

Riemann Sums and Riemann Integrals Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 203 Outline Riemnn Sums Riemnn Integrls Properties Abstrct

More information

Properties of the Riemann Integral

Properties of the Riemann Integral Properties of the Riemnn Integrl Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University Februry 15, 2018 Outline 1 Some Infimum nd Supremum Properties 2

More information

12.1 Introduction to Rational Expressions

12.1 Introduction to Rational Expressions . Introduction to Rtionl Epressions A rtionl epression is rtio of polynomils; tht is, frction tht hs polynomil s numertor nd/or denomintor. Smple rtionl epressions: 0 EVALUATING RATIONAL EXPRESSIONS To

More information

Chapter 5. , r = r 1 r 2 (1) µ = m 1 m 2. r, r 2 = R µ m 2. R(m 1 + m 2 ) + m 2 r = r 1. m 2. r = r 1. R + µ m 1

Chapter 5. , r = r 1 r 2 (1) µ = m 1 m 2. r, r 2 = R µ m 2. R(m 1 + m 2 ) + m 2 r = r 1. m 2. r = r 1. R + µ m 1 Tor Kjellsson Stockholm University Chpter 5 5. Strting with the following informtion: R = m r + m r m + m, r = r r we wnt to derive: µ = m m m + m r = R + µ m r, r = R µ m r 3 = µ m R + r, = µ m R r. 4

More information

MATH STUDENT BOOK. 10th Grade Unit 5

MATH STUDENT BOOK. 10th Grade Unit 5 MATH STUDENT BOOK 10th Grde Unit 5 Unit 5 Similr Polygons MATH 1005 Similr Polygons INTRODUCTION 3 1. PRINCIPLES OF ALGEBRA 5 RATIOS AND PROPORTIONS 5 PROPERTIES OF PROPORTIONS 11 SELF TEST 1 16 2. SIMILARITY

More information

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0) 1 Tylor polynomils In Section 3.5, we discussed how to pproximte function f(x) round point in terms of its first derivtive f (x) evluted t, tht is using the liner pproximtion f() + f ()(x ). We clled this

More information

13.4 Work done by Constant Forces

13.4 Work done by Constant Forces 13.4 Work done by Constnt Forces We will begin our discussion of the concept of work by nlyzing the motion of n object in one dimension cted on by constnt forces. Let s consider the following exmple: push

More information

1 Online Learning and Regret Minimization

1 Online Learning and Regret Minimization 2.997 Decision-Mking in Lrge-Scle Systems My 10 MIT, Spring 2004 Hndout #29 Lecture Note 24 1 Online Lerning nd Regret Minimiztion In this lecture, we consider the problem of sequentil decision mking in

More information

Dynamic Fully-Compressed Suffix Trees

Dynamic Fully-Compressed Suffix Trees Motivtion Dynmic FCST s Conclusions Dynmic Fully-Compressed Suffix Trees Luís M. S. Russo Gonzlo Nvrro Arlindo L. Oliveir INESC-ID/IST {lsr,ml}@lgos.inesc-id.pt Dept. of Computer Science, University of

More information

The solutions of the single electron Hamiltonian were shown to be Bloch wave of the form: ( ) ( ) ikr

The solutions of the single electron Hamiltonian were shown to be Bloch wave of the form: ( ) ( ) ikr Lecture #1 Progrm 1. Bloch solutions. Reciprocl spce 3. Alternte derivtion of Bloch s theorem 4. Trnsforming the serch for egenfunctions nd eigenvlues from solving PDE to finding the e-vectors nd e-vlues

More information

Faster Regular Expression Matching. Philip Bille Mikkel Thorup

Faster Regular Expression Matching. Philip Bille Mikkel Thorup Fster Regulr Expression Mtching Philip Bille Mikkel Thorup Outline Definition Applictions History tour of regulr expression mtching Thompson s lgorithm Myers lgorithm New lgorithm Results nd extensions

More information

Worksheet A EXPONENTIALS AND LOGARITHMS PMT. 1 Express each of the following in the form log a b = c. a 10 3 = 1000 b 3 4 = 81 c 256 = 2 8 d 7 0 = 1

Worksheet A EXPONENTIALS AND LOGARITHMS PMT. 1 Express each of the following in the form log a b = c. a 10 3 = 1000 b 3 4 = 81 c 256 = 2 8 d 7 0 = 1 C Worksheet A Epress ech of the following in the form log = c. 0 = 000 4 = 8 c 56 = 8 d 7 0 = e = f 5 = g 7 9 = 9 h 6 = 6 Epress ech of the following using inde nottion. log 5 5 = log 6 = 4 c 5 = log 0

More information

Physics 9 Fall 2011 Homework 2 - Solutions Friday September 2, 2011

Physics 9 Fall 2011 Homework 2 - Solutions Friday September 2, 2011 Physics 9 Fll 0 Homework - s Fridy September, 0 Mke sure your nme is on your homework, nd plese box your finl nswer. Becuse we will be giving prtil credit, be sure to ttempt ll the problems, even if you

More information

THE QUADRATIC RECIPROCITY LAW OF DUKE-HOPKINS. Circa 1870, G. Zolotarev observed that the Legendre symbol ( a p

THE QUADRATIC RECIPROCITY LAW OF DUKE-HOPKINS. Circa 1870, G. Zolotarev observed that the Legendre symbol ( a p THE QUADRATIC RECIPROCITY LAW OF DUKE-HOPKINS PETE L CLARK Circ 1870, Zolotrev observed tht the Legendre symbol ( p ) cn be interpreted s the sign of multipliction by viewed s permuttion of the set Z/pZ

More information

1. Gauss-Jacobi quadrature and Legendre polynomials. p(t)w(t)dt, p {p(x 0 ),...p(x n )} p(t)w(t)dt = w k p(x k ),

1. Gauss-Jacobi quadrature and Legendre polynomials. p(t)w(t)dt, p {p(x 0 ),...p(x n )} p(t)w(t)dt = w k p(x k ), 1. Guss-Jcobi qudrture nd Legendre polynomils Simpson s rule for evluting n integrl f(t)dt gives the correct nswer with error of bout O(n 4 ) (with constnt tht depends on f, in prticulr, it depends on

More information

COMPUTER SCIENCE TRIPOS

COMPUTER SCIENCE TRIPOS CST.2011.2.1 COMPUTER SCIENCE TRIPOS Prt IA Tuesdy 7 June 2011 1.30 to 4.30 COMPUTER SCIENCE Pper 2 Answer one question from ech of Sections A, B nd C, nd two questions from Section D. Submit the nswers

More information

along the vector 5 a) Find the plane s coordinate after 1 hour. b) Find the plane s coordinate after 2 hours. c) Find the plane s coordinate

along the vector 5 a) Find the plane s coordinate after 1 hour. b) Find the plane s coordinate after 2 hours. c) Find the plane s coordinate L8 VECTOR EQUATIONS OF LINES HL Mth - Sntowski Vector eqution of line 1 A plne strts journey t the point (4,1) moves ech hour long the vector. ) Find the plne s coordinte fter 1 hour. b) Find the plne

More information

The area under the graph of f and above the x-axis between a and b is denoted by. f(x) dx. π O

The area under the graph of f and above the x-axis between a and b is denoted by. f(x) dx. π O 1 Section 5. The Definite Integrl Suppose tht function f is continuous nd positive over n intervl [, ]. y = f(x) x The re under the grph of f nd ove the x-xis etween nd is denoted y f(x) dx nd clled the

More information

Computing the Optimal Global Alignment Value. B = n. Score of = 1 Score of = a a c g a c g a. A = n. Classical Dynamic Programming: O(n )

Computing the Optimal Global Alignment Value. B = n. Score of = 1 Score of = a a c g a c g a. A = n. Classical Dynamic Programming: O(n ) Alignment Grph Alignment Mtrix Computing the Optiml Globl Alignment Vlue An Introduction to Bioinformtics Algorithms A = n c t 2 3 c c 4 g 5 g 6 7 8 9 B = n 0 c g c g 2 3 4 5 6 7 8 t 9 0 2 3 4 5 6 7 8

More information

Solution Manual. for. Fracture Mechanics. C.T. Sun and Z.-H. Jin

Solution Manual. for. Fracture Mechanics. C.T. Sun and Z.-H. Jin Solution Mnul for Frcture Mechnics by C.T. Sun nd Z.-H. Jin Chpter rob.: ) 4 No lod is crried by rt nd rt 4. There is no strin energy stored in them. Constnt Force Boundry Condition The totl strin energy

More information

Lecture Note 9: Orthogonal Reduction

Lecture Note 9: Orthogonal Reduction MATH : Computtionl Methods of Liner Algebr 1 The Row Echelon Form Lecture Note 9: Orthogonl Reduction Our trget is to solve the norml eution: Xinyi Zeng Deprtment of Mthemticl Sciences, UTEP A t Ax = A

More information

MATH 101A: ALGEBRA I PART B: RINGS AND MODULES 35

MATH 101A: ALGEBRA I PART B: RINGS AND MODULES 35 MATH 101A: ALGEBRA I PART B: RINGS AND MODULES 35 9. Modules over PID This week we re proving the fundmentl theorem for finitely generted modules over PID, nmely tht they re ll direct sums of cyclic modules.

More information

5.7 Improper Integrals

5.7 Improper Integrals 458 pplictions of definite integrls 5.7 Improper Integrls In Section 5.4, we computed the work required to lift pylod of mss m from the surfce of moon of mss nd rdius R to height H bove the surfce of the

More information

5.5 The Substitution Rule

5.5 The Substitution Rule 5.5 The Substitution Rule Given the usefulness of the Fundmentl Theorem, we wnt some helpful methods for finding ntiderivtives. At the moment, if n nti-derivtive is not esily recognizble, then we re in

More information

Duality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below.

Duality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below. Dulity #. Second itertion for HW problem Recll our LP emple problem we hve been working on, in equlity form, is given below.,,,, 8 m F which, when written in slightly different form, is 8 F Recll tht we

More information

Math 113 Exam 2 Practice

Math 113 Exam 2 Practice Mth Em Prctice Februry, 8 Em will cover sections 6.5, 7.-7.5 nd 7.8. This sheet hs three sections. The first section will remind you bout techniques nd formuls tht you should know. The second gives number

More information

University of Alabama Department of Physics and Astronomy. PH126: Exam 1

University of Alabama Department of Physics and Astronomy. PH126: Exam 1 University of Albm Deprtment of Physics nd Astronomy PH 16 LeClir Fll 011 Instructions: PH16: Exm 1 1. Answer four of the five questions below. All problems hve equl weight.. You must show your work for

More information

Riemann Sums and Riemann Integrals

Riemann Sums and Riemann Integrals Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 2013 Outline 1 Riemnn Sums 2 Riemnn Integrls 3 Properties

More information

Nondeterminism and Nodeterministic Automata

Nondeterminism and Nodeterministic Automata Nondeterminism nd Nodeterministic Automt 61 Nondeterminism nd Nondeterministic Automt The computtionl mchine models tht we lerned in the clss re deterministic in the sense tht the next move is uniquely

More information

Chapter 0. What is the Lebesgue integral about?

Chapter 0. What is the Lebesgue integral about? Chpter 0. Wht is the Lebesgue integrl bout? The pln is to hve tutoril sheet ech week, most often on Fridy, (to be done during the clss) where you will try to get used to the ides introduced in the previous

More information

Read section 3.3, 3.4 Announcements:

Read section 3.3, 3.4 Announcements: Dte: 3/1/13 Objective: SWBAT pply properties of exponentil functions nd will pply properties of rithms. Bell Ringer: 1. f x = 3x 6, find the inverse, f 1 x., Using your grphing clcultor, Grph 1. f x,f

More information

Mapping the delta function and other Radon measures

Mapping the delta function and other Radon measures Mpping the delt function nd other Rdon mesures Notes for Mth583A, Fll 2008 November 25, 2008 Rdon mesures Consider continuous function f on the rel line with sclr vlues. It is sid to hve bounded support

More information

Non-Linear & Logistic Regression

Non-Linear & Logistic Regression Non-Liner & Logistic Regression If the sttistics re boring, then you've got the wrong numbers. Edwrd R. Tufte (Sttistics Professor, Yle University) Regression Anlyses When do we use these? PART 1: find

More information

Solution for Assignment 1 : Intro to Probability and Statistics, PAC learning

Solution for Assignment 1 : Intro to Probability and Statistics, PAC learning Solution for Assignment 1 : Intro to Probbility nd Sttistics, PAC lerning 10-701/15-781: Mchine Lerning (Fll 004) Due: Sept. 30th 004, Thursdy, Strt of clss Question 1. Bsic Probbility ( 18 pts) 1.1 (

More information

HW3, Math 307. CSUF. Spring 2007.

HW3, Math 307. CSUF. Spring 2007. HW, Mth 7. CSUF. Spring 7. Nsser M. Abbsi Spring 7 Compiled on November 5, 8 t 8:8m public Contents Section.6, problem Section.6, problem Section.6, problem 5 Section.6, problem 7 6 5 Section.6, problem

More information

The Properties of Stars

The Properties of Stars 10/11/010 The Properties of Strs sses Using Newton s Lw of Grvity to Determine the ss of Celestil ody ny two prticles in the universe ttrct ech other with force tht is directly proportionl to the product

More information

Is there an easy way to find examples of such triples? Why yes! Just look at an ordinary multiplication table to find them!

Is there an easy way to find examples of such triples? Why yes! Just look at an ordinary multiplication table to find them! PUSHING PYTHAGORAS 009 Jmes Tnton A triple of integers ( bc,, ) is clled Pythgoren triple if exmple, some clssic triples re ( 3,4,5 ), ( 5,1,13 ), ( ) fond of ( 0,1,9 ) nd ( 119,10,169 ). + b = c. For

More information

I. INTEGRAL THEOREMS. A. Introduction

I. INTEGRAL THEOREMS. A. Introduction 1 U Deprtment of Physics 301A Mechnics - I. INTEGRAL THEOREM A. Introduction The integrl theorems of mthemticl physics ll hve their origin in the ordinry fundmentl theorem of clculus, i.e. xb x df dx dx

More information

Math 360: A primitive integral and elementary functions

Math 360: A primitive integral and elementary functions Mth 360: A primitive integrl nd elementry functions D. DeTurck University of Pennsylvni October 16, 2017 D. DeTurck Mth 360 001 2017C: Integrl/functions 1 / 32 Setup for the integrl prtitions Definition:

More information

Uninformed Search Lecture 4

Uninformed Search Lecture 4 Lecture 4 Wht re common serch strtegies tht operte given only serch problem? How do they compre? 1 Agend A quick refresher DFS, BFS, ID-DFS, UCS Unifiction! 2 Serch Problem Formlism Defined vi the following

More information

Torsion in Groups of Integral Triangles

Torsion in Groups of Integral Triangles Advnces in Pure Mthemtics, 01,, 116-10 http://dxdoiorg/1046/pm011015 Pulished Online Jnury 01 (http://wwwscirporg/journl/pm) Torsion in Groups of Integrl Tringles Will Murry Deprtment of Mthemtics nd Sttistics,

More information

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying Vitli covers 1 Definition. A Vitli cover of set E R is set V of closed intervls with positive length so tht, for every δ > 0 nd every x E, there is some I V with λ(i ) < δ nd x I. 2 Lemm (Vitli covering)

More information

MAC-solutions of the nonexistent solutions of mathematical physics

MAC-solutions of the nonexistent solutions of mathematical physics Proceedings of the 4th WSEAS Interntionl Conference on Finite Differences - Finite Elements - Finite Volumes - Boundry Elements MAC-solutions of the nonexistent solutions of mthemticl physics IGO NEYGEBAUE

More information