Block Crossings in Storyline Visualizations

Size: px
Start display at page:

Download "Block Crossings in Storyline Visualizations"

Transcription

1 Journal of Graph Algorithms and Appliations vol. 21, no. 5, pp (2017) DOI: /jgaa Blok Crossings in Storylin Visualizations Thomas C. van Dijk 1 Martin Fink 2 Norbrt Fishr 1 Fabian Lipp 1 Ptr Markfldr 1 Alxandr Ravsky 3 Subhash Suri 2 Alxandr Wolff 1 1 Lhrstuhl für Informatik I, Univrsität Würzburg, Grmany 2 Univrsity of California, Santa Barbara, USA 3 Pidstryhah Institut for Applid Problms of Mhanis and Mathmatis, National Aadmy of Sins of Ukrain, Lviv, Ukrain Abstrat Storylin visualizations hlp visualiz nountrs of th haratrs in a story ovr tim. Eah haratr is rprsntd by an x-monoton urv that gos from lft to right visualizing progrssion of tim. A mting is rprsntd by having th haratrs that partiipat in th mting run los togthr for som tim. In ordr to kp th visual omplxity low, rathr than just minimizing pairwis rossings of urvs, w propos to ount blok rossings, that is, pairs of intrsting bundls of lins. In a blok rossing, two bloks of paralll lins intrst ah othr, whih is lss distrating than th sam numbr of individual rossings bing sprad ovr th drawing. In this papr, w show that minimizing th numbr of blok rossings is NPhard, vn if all mtings ar of siz 2. For this spial as, w prsnt a grdy huristi, whih w valuat xprimntally. W show that th gnral as is fixd-paramtr tratabl. Our main rsults is a onstant-fator approximation algorithm for mtings of boundd siz. Th algorithm is basd on (approximatly) solving a hyprdg dltion problm on hyprgraphs that may b of indpndnt intrst. Submittd: Dmbr 2016 Rviwd: Marh 2017 Aptd: August 2017 Artil typ: Rgular papr Rvisd: May 2017 Final: August 2017 Rviwd: Jun 2017 Publishd: Otobr 2017 Communiatd by: Y. Hu and M. Nöllnburg Rvisd: July 2017 A prliminary vrsion of this papr appard in th Prodings of th 24th Intrnational Symposium on Graph Drawing and Ntwork Visualization (GD 2016) [17]. addrsss: thomas.van.dijk@uni-wurzburg.d (Thomas C. van Dijk) fink@s.usb.du (Martin Fink) norbrt.fishr@stud-mail.uni-wurzburg.d (Norbrt Fishr) fabian.lipp@uniwurzburg.d (Fabian Lipp) ptr.markfldr@stud-mail.uni-wurzburg.d (Ptr Markfldr) alxandr.ravsky@uni-wurzburg.d (Alxandr Ravsky) suri@s.usb.du (Subhash Suri) orid.org/ x (Alxandr Wolff)

2 874 T. C. van Dijk t al. Blok Crossings in Storylin Visualizations 1 Introdution A storylin visualization is a onvnint abstration for visualizing th omplx narrativ of intrations among popl, objts, or onpts. Th motivation oms from th stting of a movi, novl, or play whr th narrativ dvlops as a squn of intronntd sns, ah involving a subst of haratrs. S Fig. 1 for an xampl. Th storylin abstration of haratrs and vnts ourring ovr tim an b usd as a mtaphor for visualizing othr situations, from physial vnts involving groups of popl mting in orporat organizations, politial ladrs managing global affairs, and groups of sholars ollaborating on rsarh to abstrat o-ourrns of topis suh as a global vnt bing ovrd on th front pags of multipl lading nws outlts, or diffrnt organizations turning thir attntion to a ommon aus. A storylin visualization maps a st of haratrs of a story to a st of urvs in th plan and a squn of mtings btwn th haratrs to rgions in th plan whr th orrsponding urvs om los to ah othr. Whil Minard s visualization of Napolon s Russian ampaign [11] an b sn as an arly (and xtrmly stark) form of storylin visualization (ombining tim and loation on a map), th urrnt form of storylin visualizations sms to hav bn invntd by Munro [14] (ompar Fig. 1), who usd it to visualiz, in a ompat way, whih substs of haratrs mt ovr th ours of a movi. Eah haratr is shown as an x-monoton urv. Mtings our at rtain tims from lft to right. A mting orrsponds to a point in tim whr th haratrs that mt ar nxt to ah othr with only small gaps btwn thm. Munro highlights mtings by undrlying thm with a gray shadd rgion, whil w us a vrtial lin for that purpos. Hn, a storylin visualization an b sn as a drawing of a hyprgraph whos vrtis ar rprsntd by th urvs and whos dgs om in at spifi points in tim. A natural objtiv for th quality of a storylin visualization is to minimiz unnssary rossings among th haratr lins. Th numbr of rossings alon, howvr, is a poor masur: two bloks of loally paralll lins rossing ah othr ar far lss distrating than an qual numbr of rossings randomly sattrd throughout th drawing. Thrfor, instad of pairwis rossings, w Figur 1: Storylin visualization for Jurassi Park by xkd [14] with svral blok rossings (on of whih w highlightd by a bold grn llips).

3 JGAA, 21(5) (2017) 875 fous on minimizing th numbr of blok rossings; ah blok rossing involvs two arbitrarily larg sts of loally paralll lins forming a rossbar, with no othr lin in th rossing ara; s Fig. 1 for an xampl. Prvious Work. Th InfoVis ommunity has quikly mbrad Munro s fftiv visualization thniqu. Kim t al. [8] usd storylins to visualiz gnalogial data; mtings orrspond to marriags and spial thniqus ar usd to indiat hild parnt rlationships. Tanahashi and Ma [15] omputd storylin visualizations automatially and showd how to adjust th gomtry of individual lins to improv th asthtis of thir visualizations. Muldr t al. [13] visualizd lustrd, dynami graphs as storylins, summarizing th bhavior of loal ntworks that surround usr-sltd foi. Only rntly, a mor thortial and prinipld study was initiatd by Kostitsyna t al. [10], who onsidrd th problm of minimizing pairwis rossings (that is, not blok rossings) in storylins. Thy provd that th problm is NP-hard in gnral, and showd that it is fixd-paramtr tratabl (FPT) with rspt to th (total) numbr of haratrs. Thir FPT algorithm runs in tim O(k! 2 k log k + k! 2 n) and uss spa O(k! 2 + n), whr k is th numbr of haratrs and n is th numbr of mtings. For th spial as of 2-haratr mtings without rptitions whr th mting graph whos dgs dsrib th mtings of pairs of haratrs is a tr, thy showd that Θ(k log k) rossings always suffi and ar somtims nssary. Vry rntly, Gronmann t al. [7] formulatd an intgr linar program that minimizs th total numbr of rossings for storylins with mtings of arbitrary ardinality. Thir approah an solv instans with haratrs and up to about 50 mtings optimally in a fw sonds. Du to th natur of intgr linar programming, howvr, thir approah boms unusabl for larg instans. Our work builds on th problm formulation of Kostitsyna t al. [10] but w onsidrably xtnd thir rsults by dsigning (approximation) algorithms for gnral mtings for a diffrnt optimization goal: w minimiz th numbr of blok rossings rathr than th numbr of pairwis lin rossings. Blok rossings wr introdud by Fink t al. [5] for visualizing mtro maps; in thir stting, blok rossings happn btwn mtro lins that run on top of an mbddd graph, th mtro ntwork. Problm Dfinition. Whil a blok rossing is visually asily rognizabl by a human radr, w start by arfully dfining th problm. In th ours of doing so, w will alrady gain furthr insight in what minimizing th numbr of blok rossings amounts to. A storylin S is a pair (C, M) whr C = {1,..., k} is a st of haratrs and M = [m 1, m 2,..., m n ] with m i C and m i 2 for i = 1, 2,..., n is a squn of mtings of at last two haratrs. W all any st g C of haratrs that has at last on mting, a group (i.., g = m i for som 1 i n). W dfin th group hyprgraph H = (C, Γ) as follows: its vrtis ar th haratrs and its hyprdgs ar th groups. This hyprgraph dos

4 876 T. C. van Dijk t al. Blok Crossings in Storylin Visualizations 1 a b b+1 k Figur 2: Blok rossing (a, b, ) not inlud th tmporal aspt of th storylin it modls only th strutur of groups partiipating in th storylin. Eah group is rprsntd by a singl hyprdg, vn if it mts multipl tims. Th group hyprgraph an b built by sorting M and filtring th multiply-ourring mtings. This an b don in O(nk) tim by rprsnting th mtings as bitsts and using radix sort. Not that w do not nod th xat tims of th mtings: in a givn visualization, at any tim t, thr is a uniqu vrtial ordr π t of th haratrs. Without hanging π t by rossings, w an inras or dras vrtial gaps btwn lins. If a group g forms a ontiguous intrval in π t, thn w an bring th lins in g within a short vrtial distan δ group without any rossing, and also mak sur that all othr lins ar at a distan of δ sp > δ group ; s Fig. 12. Sin any group must b supportd at a tim just bfor its mting starts, omputing an output drawing onsists mainly of hanging th prmutation of haratrs ovr tim so that during a mting its group is supportd by th urrnt prmutation. W thrfor fous on hanging th prmutation by rossings ovr tim, and only hav to b onrnd about th ordr of mtings; th final drawing an b obtaind by a simpl post-prossing from this disrt st of prmutations. Suppos that w loally rnumbr th haratrs from top to bottom from 1 to k. Thn w an dfin, for a b <, a blok rossing (a, b, ) to b th xhang of th two onsutiv bloks a,..., b and b + 1,...,. This xhang maps th prmutation 1,..., a,..., b,...,,..., k to th prmutation 1,..., a 1, b + 1,...,, a,..., b, + 1,..., k ; s Fig. 2 for an illustration. Th sam dfinition was usd by Fink t al. [5]. A mting m fits a prmutation π (or a prmutation π supports a mting m) if th haratrs partiipating in m form an intrval in π. In othr words, thr is a prmutation of m that is part of π. If w apply a squn B of blok rossings to a prmutation π in th givn ordr, w dnot th rsulting prmutation by B(π). Problm 1 (Storylin Blok Crossing Minimization (SBCM)) Givn a storylin instan (C, M) with M = [m 1, m 2,..., m n ], find a solution onsisting of a start prmutation π 0 of C and a squn B = (B 1, B 2,..., B n ) of (possibly mpty) squns of blok rossings suh that th total numbr of blok rossings is minimizd and π i = B i (π i 1 ) supports m i for 1 i n.

5 JGAA, 21(5) (2017) 877 Not that in our problm dfinition no haratrs arriv aftr th start of th drawing or disappar bfor th nd; this is diffrnt in th xampl in Fig. 1. W also onsidr d-sbcm, a spial as of SBCM whr mtings involv groups of siz at most d, for an arbitrary onstant d. For xampl, 2-SBCM allows only 2-haratr mtings, a stting that was also studid by Kostitsyna t al. [10]. Our Rsults. W obsrv that a storylin has a rossing-fr visualization if and only if its group hyprgraph is an intrval hyprgraph, that is, if thr xists a prmutation of th haratr st suh that ah hyprdg orrsponds to a ontiguous blok of haratrs in this prmutation. This is quivalnt to th hyprgraph having path support [1]. A hyprgraph an b tstd for th intrval proprty in O(n 2 ) tim, whr n is th numbr of hyprdgs. W show that 2-SBCM is NP-hard (s St. 3) and that SBCM is fixd-paramtr tratabl with rspt to k (St. 4). On of our FPT algorithms an b modifid to handl pairwis rossings instad of blok rossings. Th modifid algorithm is fastr than th algorithm of Kostitsyna t al. [10] by a fator of k!. W arry out som xprimnts valuating this algorithm and othr approahs in a follow-up papr [18]. Th as of siz-2 mtings is of intrst; rall that Kim t al. [8] hav usd a storylin-lik visualization to tra gnalogial data. W invstigat strutural proprtis and prsnt a grdy algorithm for 2-SBCM that runs in O(k 3 n) tim for k haratrs. For k = 3, th grdy algorithm yilds optimal solutions. W xprimntally ompar grdy solutions to optimal solutions; s St. 5. On main rsult is a onstant-fator approximation algorithm for d-sbcm for th as that d is boundd and that mtings annot b rpatd; s St. 6. Our algorithm is basd on a solution for th following NP-omplt hyprgraph problm, whih may b of indpndnt intrst: givn a hyprgraph H, dlt th minimum numbr of hyprdgs so that th rmaindr is an intrval hyprgraph. For this problm, w dvlop a (d + 1)-approximation algorithm, whr d is th maximum ardinality of a hyprdg in H; s St. 7. Finally, w list som opn problms in St Prliminaris and Basi Obsrvations First, w onsidr th spial as whr vry mting onsists of two haratrs. For ths rstritd instans, any mting an b ralizd from any prmutation by a singl blok rossing. Obsrvation 1 Givn an instan of 2-SBCM, thr is a solution with at most on blok rossing bfor ah of th mtings. In partiular, thr is a solution with at most n blok rossings in total. Proof: Lt π b an arbitrary prmutation, and lt m = {, } M b th nxt mting. Lt i and j b th positions of th haratrs and, rsptivly, in th prmutation π, that is, π (i) = and π (j) =. Without loss of gnrality,

6 878 T. C. van Dijk t al. Blok Crossings in Storylin Visualizations π 0 π 1 π Figur 3: Optimal solution for S from th proof of Lmma 1. assum i < j. If π dos not support m, w an raliz it using th blok rossings (i, i, j 1), that is, moving th lin of dirtly abov that of. Obsrvation 1 shows that thr is a solution with at most on blok rossing bfor any mting. This raiss th qustion whthr thr is also an optimal solution that fulfills this ondition. Th answr is ngativ. Lmma 1 Thr is an instan S of 2-SBCM and a start prmutation π 0 suh that thr is no optimal solution (π 0, B) of S that starts with π 0 and uss at most on blok rossing bfor th first and btwn ah pair of onsutiv mtings. Proof: W prov th laim by ontradition. Considr th instan S = (C, M) with C = {1, 2, 3, 4, 5, 6, 7, 8} and M = [{6, 3}, {7, 2}, {1, 5}, {5, 6}, {6, 3}, {3, 4}, {4, 8}, {8, 7}]. Lt π 0 = 1, 2, 3, 4, 5, 6, 7, 8 b th start prmutation. Thr is a solution with only two blok rossings, namly (π 0, B) with B = [(2, 4, 7), (4, 5, 8)]; s Fig. 3. Lt π 1 b th prmutation aftr th first blok rossing of B on π 0, and π 2 th prmutation aftr both blok rossings. Th prmutation π 2 supports all mtings in M. Th first mting {6, 3} in M fits nithr π 0 nor π 1, that is, both blok rossings must our bfor th first mting. Now assum that thr is anothr solution (π 0, B ) with B 2 that has at most on blok rossing bfor ah mting. Starting from π 0 thr ar xatly nin fasibl blok rossings that allow th first mting. Thy yild th following prmutations (whr th bars indiat th rossing bloks): 1, 2, 4, 5, 6, 3, 7, 8 4, 5, 1, 2, 3, 6, 7, 8 1, 2, 3, 6, 4, 5, 7, 8 1, 2, 5, 6, 3, 4, 7, 8 1, 4, 5, 2, 3, 6, 7, 8 1, 2, 3, 6, 7, 4, 5, 8 1, 2, 6, 3, 4, 5, 7, 8 1, 2, 4, 5, 3, 6, 7, 8 1, 2, 3, 6, 7, 8, 4, 5 Non of ths prmutations supports th sond mting {7, 2}. So w nd th sond blok rossing bfor this mting. This sond blok rossing nds to prpar all of th rmaining mtings baus othrwis B > 2. Ths mtings an only b supportd by th prmutation σ = 1, 5, 6, 3, 4, 8, 7, 2 or

7 JGAA, 21(5) (2017) 879 m π 0 π 1 π Figur 4: Optimal solution for S from th proof of Thorm 1. its rvrs prmutation σ R. It rmains to show that non of th prmutations yildd by th fasibl first blok rossing an b transformd to σ or σ R by on additional blok rossing. All prmutations ontaining 3, 6 as a subsqun ar infasibl baus thr is only on blok rossing that swaps two nighboring haratrs and it dos not produ σ. For prmutations starting with 1, 2, thr is only on possibl blok rossing to bring haratr 2 to th nd of th prmutation whil haratr 1 stays at th first position, whih also dos not yild σ. Similarly, w an show that thr is also no blok rossing aftr any of th fasibl blok rossings for th first stp that lads to σ R. Now w modify this instan suh that w do not nd to prsrib th start prmutation. Thorm 1 Thr is an instan S of 2-SBCM suh that thr is no optimal solution (π 0, B) of S that uss at most on blok rossing btwn ah pair of onsutiv mtings. Proof: W prov th laim by onstruting an xampl. Considr th instan S = (C, M) dpitd in Fig. 4. Th instan has no solution with lss than two blok rossings, but thr is a solution with two blok rossings, namly (π 0, B) with start prmutation π 0 = 1, 2, 3, 4, 5, 6, 7, 8 and blok rossings B = [(2, 4, 7), (4, 5, 8)]; s Fig. 4. Lt π 1 b th prmutation aftr th first blok rossing of B on π 0, and lt π 2 = 1, 5, 6, 3, 4, 8, 7, 2 b th prmutation aftr th sond blok rossing in B. Th prmutation π 2 supports all mtings in th right half of M. Th first ourrn of mting {6, 3} in M (right aftr th third dashd lin in Fig. 4) fits nithr π 0 nor π 1, that is, both blok rossings our bfor that mting. Now w argu that thr is no solution with two blok rossings that ar sparatd by at last on mting. Lt m b th third ourrn of mting {1, 2} (right bfor th first dashd lin in Fig. 4). At som point in tim bfor m, ithr prmutation π 0 or π0 R must our, othrwis w would nd at last two blok rossings to raliz th mtings bfor m and at last on blok rossing aftr m to raliz th last svn mtings (baus th prmutation at m is not π 2 or π2 R ). W distinguish thr ass for th two blok rossings.

8 880 T. C. van Dijk t al. Blok Crossings in Storylin Visualizations a) Both rossings our bfor m: This is not possibl baus th rmaining mtings annot b ralizd with haratrs 1 and 2 bing adjant. b) Both rossings our aftr m: Thn th start prmutation nds to b π 0 or π R 0, whih is xatly th as handld in Lmma 1. ) Th first rossing is bfor m, th sond on aftr m. Sin w must raliz th last svn mtings with only on blok rossing, th final prmutation is π 2 or π R 2. Assum that it is π 2. Thn going bakwards from π 2, w hav to raliz m by on blok rossing. Thr ar xatly nin suh blok rossings yilding th following prmutations: 5, 6, 3, 4, 8, 7, 2, 1 6, 3, 4, 8, 7, 2, 1, 5 3, 4, 8, 7, 2, 1, 5, 6 4, 8, 7, 2, 1, 5, 6, 3 8, 7, 2, 1, 5, 6, 3, 4 7, 2, 1, 5, 6, 3, 4, 8 2, 1, 5, 6, 3, 4, 8, 7 1, 2, 5, 6, 3, 4, 8, 7 5, 6, 3, 4, 8, 7, 1, 2 It is asy to hk that non of thm ontains as a subintrval what w all a triplt; thr onsutiv numbrs in thir natural (or invrtd) ordr (for instan, 3, 4, 5 or 5, 4, 3 ). Rall that at som point in tim bfor m, prmutation π 0 or π0 R must our. Now not that any prmutation rahabl from π 0 or π0 R by on blok rossing ontains a triplt (sin w hav ight haratrs). Thus w obtain a ontradition. Th as whn th final prmutation is π2 R an b tratd similarly. Dtting Crossing-Fr Storylins. If a storylin admits a rossing-fr visualization, thn th vrtial prmutation of th haratr lins rmains th sam ovr tim, and all mtings involv groups that form ontiguous substs in that prmutation. (Th visualization an b obtaind by plaing haratrs along a vrtial lin in th orrt prmutation and for ah mting bringing its lins togthr for th duration of th mting and thn sparating thm apart again.) In othr words, a singl prmutation supports ah group of H = (C, Γ). This holds if and only if H is an intrval hyprgraph. Rall that this is th as if thr xists a prmutation π = v 1,..., v k of C suh that ah hyprdg Γ orrsponds to a ontiguous blok of haratrs in this prmutation. Not that a graph is an intrval hyprgraph if and only if all of th graph s onntd omponnts ar paths. Hn, a graph bing an intrval hyprgraph is not th sam proprty as bing an intrval graph. An intrval hyprgraph an b visualizd by plaing all of its vrtis on a lin, and drawing ah hyprdg as an intrval that inluds all vrtis of and no vrtx of V \. Chking whthr a k-vrtx hyprgraph is an intrval hyprgraph taks O(k 2 ) tim [16]. Rall that w an build H in O(nk log n) tim. Thorm 2 Givn th group hyprgraph H of an instan of SBCM with k haratrs, w an hk in O(k 2 ) tim whthr a rossing-fr solution xists. If this is th as, a prmutation ralizing a rossing-fr solution an also b found within th sam tim bound.

9 JGAA, 21(5) (2017) 881 For 2-SBCM w only nd to hk (in O(k) tim) whthr H is a olltion of vrtx-disjoint paths; this is dominatd by th tim (O(n)) for building H, as long as k = O(n). 3 NP-Compltnss of SBCM In this stion w prov that SBCM is NP-omplt. This is known for th problm Blok Crossing Minimization (BCM), introdud by Fink t al. [5]. But SBCM is not simply a gnralization of BCM baus in SBCM w an hoos an arbitrary start prmutation. Thrfor, th ida of our hardnss proof is to for a rtain start prmutation by adding som haratrs and mtings. W rdu from Sorting by Transpositions (SBT), whih has also bn usd to show th hardnss of BCM [5]. In SBT, on is givn a prmutation π and an intgr k, and th task is to did whthr thr is a squn of transpositions (whih ar quivalnt to blok rossings) of lngth at most k that transforms π to th idntity. SBT was rntly shown NP-hard by Bultau t al. [2]. W show hardnss for 2-SBCM, whih implis that SBCM is NP-hard, too. It is asy to s that th dision vrsion of SBCM is in NP: Obviously, th maximum numbr of blok rossings ndd for any numbr of haratrs and mtings is boundd by a polynomial in k and n. Thrfor also th siz of a possibl rtifiat is boundd by a polynomial. To tst th fasibility of a solution ffiintly, w simply tst whthr th prmutations btwn th blok rossings support th mtings in th right ordr from lft to right. W will us th following obvious fat. Obsrvation 2 If prmutation π nds blok rossings to b sortd, any prmutation ontaining π as a subsqun nds at last blok rossings to b sortd. Thorm 3 2-SBCM is NP-omplt. Proof: It rmains to show th NP-hardnss. W rdu from SBT. Givn an instan of SBT, that is, a prmutation π of {1,..., k}, w show how to onstrut a orrsponding instan of 2-SBCM. W xtnd th st of haratrs {1, 2,..., k} to C = {1,..., k, 1, 2,..., 2k }. Corrspondingly, w xtnd th prmutation π = π(1), π(2),..., π(k) to π = 1,..., 2k, π(1),..., π(k) and ι to ι = 1, 2,..., 2k, 1, 2,..., k. Lt M π and M ι b th squns of mtings of all nighboring pairs in π and ι, rsptivly. Lt M 1 and M 2 b th onatnations of k + 1 opis of M π and M ι, rsptivly. This yilds th instan S = (C, M) of 2-SBCM, whr M is th onatnation of M 1 and M 2 ; s Fig. 5. W show that th numbr of blok rossings ndd for th 2-SBCM instan S quals th numbr of transpositions to solv instan π of SBT, that is, to transform π to th idntity ι = 1, 2,..., k. Not that π an b sortd by at most k blok rossings. So k is an uppr bound for an optimal solution of instan π of SBT.

10 882 T. C. van Dijk t al. Blok Crossings in Storylin Visualizations M 1 M π π π 2k π 1 π k k 1 k 1 B 2k 1 k ι ι ι Figur 5: Solution for th 2-SBCM instan S orrsponding to a solution B of instan π of SBT. Th box B rprsnts th blok rossings. First, lt B b a shortst squn of blok rossings to sort π. Thn, (π, B) is a fasibl solution for S. Th start prmutation π supports all mtings in M 1 without any blok rossing. Using B, th lins ar sortd to ι, and this prmutation supports all mtings in M 2 without any furthr blok rossings; s Fig. 5. Hn, th numbr of blok rossings in any solution of π is an uppr bound for th minimum numbr of blok rossings ndd for S. For th othr dirtion, lt (π, B ) b an optimal solution for S. In th input of SBCM, w do not fix a rtain ordr of th haratrs. So any solution of 2-SBCM givs ris to a symmtri solution that is obtaind by rvrsing th ordr of th haratrs. In th following, without loss of gnrality, w assum that π (rathr than th rvrs prmutation π R ) ours in M, that is, thr is a tim t in th visualization suh that th vrtial ordr of th haratrs at t is dsribd by π. Nxt, w show that th start prmutation π ours somwhr in M 1 and that ι ours somwhr in M 2. If thr is a squn M π of mtings btwn whih thr is no blok rossing, th prmutation at this position an only b th start prmutation π or its rvrs. For a ontradition, assum that π dos not our during M 1 in th layout indud by (π, B ). Thn thr is no suh squn without any blok rossing in it. As this squn is rpatd k + 1 tims, th solution would nd at last k + 1 blok rossings. This ontradits our uppr bound, whih is k. Analogously, w an show that th prmutation ι or its rvrs ours in M 2. W now want to show that th unrvrsd vrsion of ι ours in M 2. For a ontradition, assum th opposit. W forgt about th lins 1,..., k and only onsidr th squn π = 1,..., 2k in π whih is rvrsd to ι R = 2k,..., 1 in ι R. Eriksson t al. [4] showd that w nd (l + 1)/2 blok rossings to rvrs a prmutation of l lmnts. This implis that w nd k + 1 blok rossings to transform π to ι R. As π and ι R ontain ths squns as subsquns, Obsrvation 2 implis that th transformation from π to ι R also nds at last k + 1 blok rossings. As th optimal solution uss at most k blok rossings, w know that w annot rah ι R and thus th squn of prmutations ontains π and ι.

11 JGAA, 21(5) (2017) 883 Th squn of blok rossings that transforms π to ι yilds a squn B of blok rossings of th sam lngth that transforms π to ι. Thrfor th lngth of a solution for S is an uppr bound for th lngth of an optimal solution of th orrsponding SBT instan π. Thus, th two lngths ar qual. Hardnss Without Rptitions. With arbitrarily larg mtings, w an slightly modify our hardnss proof, and show that minimizing th numbr of blok rossings is also hard without rpating th sam mting many tims. Th ida to hang our rdud SBCM instan, is to rpla th rpatd squn of 2-haratr mtings so that in ah rptition th group siz is inrasd by on for all mtings; s Fig. 6. Du to th ovrlapping strutur of th groups in a singl squn, thy an only b all supportd at th sam tim if also th 2-haratr mtings that thy rplad ar supportd. Th only thing that w hav to b arful about is that whn th groups gt largr than 3k/2, that is, half of th numbr of haratrs, thr is a growing st of haratrs in π (or ι ) that ar ontaind in xatly th sam groups, and thir rlativ ordr dos not mattr for th mtings; s Fig. 6(a). W will avoid this situation in th following way (shown in Fig. 6(b)). Sin w hav k + 1 squns of rpatd mtings at th bginning as wll as at th nd of th timlin, and w kp inrasing th group sizs, w hav groups of 2k + 3 haratrs in th nd. W rpla 1,..., 2k by a nw squn 1,..., 10k of haratrs without hanging th strutur furthr. Thn, w an inras th group siz up to 2k + 3 whil in th nd still lss than half of all haratrs ar involvd in ah group. Sin th growing mtings ompltly simulat th dsird 2-haratr mtings, th rst of th rdution and its proof stay th sam, and w gt th following rsult. Thorm 4 SBCM is NP-hard vn if mtings ar not rpatd. Not that in this rdution (diffrnt from Thorm 3) w us mtings whos siz dpnds on th input. So for this variant without rpatd mtings w do not show hardnss for d-sbcm for any fixd d. (a) Instan in whih th groups ar too larg for th numbr of haratrs. Th ara highlightd in gray shows th haratrs involvd in all mtings of a rtain siz. Thrfor, th rlativ ordr of ths haratrs dos not mattr anymor. (b) By adding haratrs (and th orrsponding mtings), w avoid haratrs bing involvd in vry mting. Figur 6: W simulat rpatd 2-haratr mtings by using groups of inrasing siz.

12 884 T. C. van Dijk t al. Blok Crossings in Storylin Visualizations 4 Exat Algorithms W prsnt two xat algorithms for SBCM. Conptually, both build up a squn of blok rossings whil kping trak of how many mtings hav alrady bn aomplishd. Th first uss polynomial spa; th sond improvs th runtim at th ost of xponntial spa. Th lattr gnralizs to improv a rsult by Kostitsyna t al. [10] about minimizing pairwis rossings (rathr than blok rossings). W start with a data strutur to kp trak of prmutations, blok rossings and mtings. It is initializd with a givn prmutation and has two oprations. Th Chk opration rturns whthr a givn mting fits th urrnt prmutation. Th BlokMov opration prforms a givn blok rossing on th prmutation and thn rturns whthr th most-rntly Chkd mting now fits. Lmma 2 A squn of arbitrarily intrlavd BlokMov and Chk oprations an b prformd in O(β + µ) tim, whr β is th numbr of alls to BlokMov and µ is sum of ardinalitis of th mtings givn to Chk. Th spa usag is O(k), whr k is th numbr of haratrs. Proof: Rprsnt th prmutation as a doubly-linkd list. 1 First onsidr a 2- mting. It taks onstant tim to hk whthr it fits: hk th prvious/nxt pointrs. Sin a blok rossing hangs at most 6 adjanis, BlokMov an updat th linkd list in onstant tim. (This rquirs giving out pointrs into th linkd-list rprsntation, sin w nd to find th haratrs at th dgs of th blok in onstant tim.) Not that all possibl blok rossings an b numratd in onstant tim ah. Now w look at a mting of ardinality m. Intrprt th linkd list as a path and onsidr th subgraph indud by th nods in th mting. If th mting fits th prmutation, this subgraph is onntd and, bing a path, has m 1 dgs; if th mting dos not fit, this subgraph has mor omponnts and thrfor fwr dgs. Th Chk opration on a mting of siz m an b prformd in O(m) tim by ounting at vry nod in th mting whthr zro, on or two of its nighbors ar also in th mting. For th amortizd runtim ovr a squn of oprations, rmmbr this ount: BlokMov an updat it in onstant tim, sin again at most six adjanis hang. In trms of spa, thr is only th doubly linkd list and th ount. Now w provid an output-snsitiv algorithm for SBCM, th runtim of whih dpnds on th numbr of blok rossings rquird by th optimum. Thorm 5 An instan S = (C, M) of SBCM an b solvd in O(k! ( k3 k 6 ) β (β + µ)) tim and O(βk) working spa if a solution with β blok rossings xists, whr µ = m M m. 1 W assum that th mtings givn to Chk ar rprsntd by rfrns to th nods in this list; if nssary, this rprsntation an b onstrutd ffiintly in prprossing.

13 JGAA, 21(5) (2017) 885 Proof: Considr a branhing algorithm that starts from a prmutation of th haratrs and kps trying all possibl blok rossings. Th instan of SBCM is solvd by finding th shortst squn of blok rossings that supports th mtings. A blok rossing an b rprsntd by indis (a, b, ) with 1 a b < k. Th numbr of possibl blok rossings is givn by adding th numbr of blok rossings with a = b and th numbr of blok rossings with a b; hn, thr ar ( ( k 2) + k ) 3 = k 3 k 6 distint blok rossings on a prmutation of lngth k. W an numrat ths in onstant tim ah by numrating all appropriat tripls (a, b, ). W us dpth-first itrativ-dpning sarh [9] from all possibl start prmutations, applying blok rossings, until w find a squn of prmutations that fulfills all mtings. Whil sarhing, th algorithm kps trak of how many mtings fit th urrnt squn of prmutations using th data strutur from Lmma 2. Corrtnss follows from th itrativ dpning of th sarh, sin w want an (unwightd) shortst squn of blok rossings. Th runtim and spa bounds follow from th standard analysis of itrativ-dpning sarh, obsrving that a nod uss O(k) spa and it taks O(β + µ) tim in total to valuat a path from root to laf. Not that µ is O(kn): thr ar n mtings and ah onsists of at most k haratrs. At th ost of xponntial spa, w an improv th runtim and gt rid of th dpndn on β, showing th problm to b fixd paramtr linar for k. Th following algorithm an asily b adaptd to optimiz pairwis rossings rathr than blok rossings. In that as w improv upon th algorithm of Kostitsyna t al. [10] by a fator of roughly k!, in trms of both running tim and spa onsumption: rall that thir FPT algorithm runs in O(k! 2 k log k + k! 2 n) tim and uss O(k! 2 + n) spa. Thorm 6 An instan of SBCM an b solvd in O(k! k 3 n) tim and O(k! k n) spa. Proof: Th algorithm is basd on th following qustion: how many blok rossings ar rquird for th first l mtings givn that w nd on prmutation π? Formally, lt f(π, l) b th optimal numbr of blok rossings in a solution to th givn instan whn rstritd to th first l mtings and to hav π as its final prmutation. Not that by dfinition th solution for th atual instan is givn by min π f(π, n), whr th minimum rangs ovr all possibl prmutations. As a bas as, f(π, 0) = 0 for all π, sin th mpty st of mtings is supportd by any prmutation. Lt π and π b prmutations that ar on blok rossing apart and lt 0 l l. If th mtings {m l+1,..., m l } all fit π, thn w hav f(π, l ) f(π, l) + 1: if w an support th first l mtings and nd on π, thn with on additional blok rossing w an support th first l mtings and nd on π. W now modl th abov obsrvation in a graph. Lt G b an unwightd dirtd graph on nods (π, l), whr π is a prmutation of haratrs and

14 886 T. C. van Dijk t al. Blok Crossings in Storylin Visualizations 0 l n. Call a nod start nod if l = 0. Thr is an ar from (π, l) to (π, l ) if and only if π and π ar on blok rossing apart, l l, and th mtings {m l+1,..., m l } fit π. Not that w allow l = l sin w may nd to allow blok rossings that do not immdiatly ahiv an additional mting (f. Proposition 1), so G is not ayli. By onstrution w hav that f(π, l) quals th minimum graph distan from a start nod to th nod (π, l). Call a path from a start nod that ralizs this distan optimal. In G, onsidr any lngth-3 path [(π 1, l 1 ), (π 2, l 2 ), (π 3, l 3 )] with strit inquality l 3 > l 2. If mting l fits π 2, thn [(π 1, l 1 ), (π 2, l 2 + 1), (π 3, l 3 )] is also a path. Rpating this transformation shows that for all π, th nod (π, n) has an optimal path in whih vry ar maximally inrass l. Lt G b th graph whr w drop all ars from G that do not maximally inras l. Not that G still ontains a path that orrsponds to th global optimum. Th graph G has O(k! n) nods. Eah nod has outdgr O(k 3 ), sin any blok rossing ontributs at most on out-ar to a nod. Thn a bradth-first sarh from all start nods to any nod (π, n) ahivs th laimd tim and spa bounds, assuming w an numrat th outgoing ars of a nod in onstant tim ah. For a givn nod (π, l) w an numrat all possibl blok rossings in onstant tim ah, as bfor. To gnrat its outgoing ars in G, w also nd to know th maximum l suh that all mtings l + 1 up to l fit π, whr π is th prmutation rsulting from th blok rossing. Not that l only dpnds on l and π ; in partiular, it dos not dpnd on π. W an thrfor promput a tabl M(π, l) that givs this valu. Computing M(π, l) for a givn π and all l taks a total of O(kn) tim: first omput for vry m i whthr it fits π, thn omput th implid forward pointrs using a linar san. So using O(k! kn) prprossing tim and O(k! n) spa, w hav an ffiint implmntation of th bradth-first sarh. Th thorm follows. 5 A Grdy Huristi In this stion w dvlop an O(kn)-tim grdy algorithm to quikly draw good storylin visualizations for 2-SBCM. Givn an instan S = (C, M), w rsrv a list B = [ ] that th algorithm will us to stor th blok rossings. Th algorithm starts with an arbitrary prmutation π 0 of th haratrs. In vry stp th algorithm rmovs all mtings from th bginning of M that ar supportd by th urrnt prmutation π i. Subsquntly, th algorithm piks a blok rossing b suh that th rsulting prmutation π i+1 = b(π i ) supports th maximum numbr of mtings from th bginning of M, and b is appndd to th list B. This pross rpats until M is mpty. Th algorithm rturns th solution (π 0, B).

15 JGAA, 21(5) (2017) (a) Grdy solution (b) Optimal solution Figur 7: Th grdy algorithm is not optimal. Rall that, from any partiular prmutation, thr ar O(k 3 ) blok rossings. To find th appropriat blok rossings, th algorithm ould simply hk all of thm. Howvr, most of thos will rsult in prmutations that do not vn support th nxt mting, whih annot b th grdy hoi dsribd abov. Hn, our algorithm numrats only th rlvant blok rossings in th following sns: blok rossings yilding a prmutation that supports th upoming mting. Lt {, } b th nxt mting in M. If x and y ar th positions of and in th urrnt prmutation (i.., π i (x) = and π i (y) = ; without loss of gnrality, assum x < y), th rlvant blok rossings ar: {(z, x, y 1): 1 z x} {(x, z, y): x z < y} {(x + 1, y 1, z): y z k}. W s that th numbr of rlvant blok rossings in ah stp is k + 1. Lt n i b th maximum numbr of mtings at th bginning of M w an support by on of ths blok rossings. W us th data strutur in Lmma 2 and hk for ah rlvant blok rossing how many mtings an b don with this prmutation. Hn, w an idntify a blok rossing ahiving th maximum numbr in O(kn i ) tim sin w hav to hk k + 1 options, ah lasting up to n i mtings ah. Th numbrs of mtings n i in ah itration of th algorithm sum up to n and thrfor th algorithm runs in O(kn) total tim. In th dsription abov, this grdy algorithm starts with an arbitrary prmutation. Instad, w ould start with a prmutation that supports th maximum numbr of mtings bfor th first blok rossing nds to b don. In othr words, w ould to find a maximal prfix M of M suh that (C, M ) an b rprsntd without any blok rossings. W an find M in O(kn) tim: start with an mpty graph on th haratrs and sussivly add dgs for th mtings. Aftr ah addition w hk whthr th graph is still a olltion of paths, whih an b don in O(k) tim. Aftr this pross trminats, w onstrut a prmutation that supports all mtings in M, whih is asy givn th olltion of paths. S Fig. 7 for an xampl that uss th huristi start prmutation. Whil this is a snsibl huristi, w do not prov that it rdus th total numbr of blok rossings. Indd, w xprimntally obsrv that, whil this huristi for th start prmutation is gnrally good, it is not always th bst; this is disussd in th xprimntal valuation latr in this stion. Th grdy algorithm atually yilds optimal solutions for spial ass of 2-SBCM. For th following lmma, w assum that no two subsqunt mtings

16 888 T. C. van Dijk t al. Blok Crossings in Storylin Visualizations in th input ar th sam. W all an instan normal if this is th as. An instan an b normalizd by simply dropping th rpatd mtings. This dos not afft th optimal numbr of blok rossings or th bhavior of th grdy algorithm, but not that it dos lowr n. Lmma 3 A normal instan of 2-SBCM with k = 3 an b solvd using at most n/2 1 blok rossings. Proof: Not that thr ar only thr possibl mtings, namly {1, 2}, {1, 3}, and {2, 3}. Any prmutation supports prisly two of ths and not th third, and is quivalnt in this sns to its rvrs. For xampl, th prmutation 1, 2, 3 and its rvrs both support th mtings {1, 2} and {2, 3}, but not {1, 3}. Lt π and π b distint prmutations. Cas distintion shows that it is always possibl with a singl blok rossing to gt from π to ithr π or to th rvrs of π. For th analysis, w partition th squn of mtings into pohs as follows. W start from th first mting and kp going until th third distint mting ours: ths mtings form th first poh. That is, an poh altrnats btwn two diffrnt mtings. Rpating this pross partitions th ntir squn of mtings into pohs, possibly with a singl rmaining mting as final poh. A solution an hoos th start prmutation π 0 that supports th first poh. Aftr that it an always gt to a prmutation that supports th ntir nxt poh with on blok rossing. (Not that th grdy algorithm uss xatly this stratgy.) In th worst as all pohs hav lngth 2, and w nd n/2 1 blok rossings. Sin th grdy algorithm an hoos an arbitrary start prmutation it dos not nd to prform a blok rossing for th first squn of mtings baus it an build a start prmutation using th sam stratgy as dsribd abov. This lavs us with at most n/2 1 blok rossings. Thorm 7 For k = 3, th grdy algorithm produs optimal solutions. Proof: W look at th pohs from Lmma 3 again. Th grdy algorithm produs on blok rossing fwr than th numbr of pohs. Considr any poh xpt th last on and inlud th mting aftr it. By onstrution, this is th third distint mting and thrfor ths mtings togthr annot fit a singl prmutation. Thn in any solution to th problm, a blok rossing must our aftr at last on of th mtings in th poh. This holds for all pohs xpt th last on and sin thy ar disjoint, th numbr of pohs rdud by on is a lowr bound for th optimum numbr of blok rossings. Th rsult of th grdy algorithm ralizs this bound. W not hr that th grdy algorithm dos not snsibly gnraliz to SBCM with arbitrary mtings. Considr th stp that finds a blok rossing that supports th maximum numbr of subsqunt mtings. With arbitrary mtings (as opposd to only ardinality-2 mtings), it may b th as that this maximum is zro that is, thr may not xist any singl blok rossing that supports th nxt mting. In this as our algorithm has no grdy way mak progrss.

17 JGAA, 21(5) (2017) 889 Tim (in s.) n Figur 8: Runtim of th xat algorithm of Thorm 5 on random instans with k = 4( ), 5(+), 6( ), 7( ). Eah data point is th avrag of 50 random instans. Exprimntal Evaluation. In this stion, w rport on som prliminary xprimntal rsults for 2-SBCM. W hav gnratd random instans as follows. Givn n and k, w pik n pairs of haratrs as mtings, uniformly at random using rjtion sampling to nsur that onsutiv mtings ar diffrnt. (This mans that th gnratd instans ar normal in th sns of Lmma 3.) First, w onsidr th xat algorithm of Thorm 5. As xptd, its runtim dpnds havily on k (Fig. 8). Prhaps somwhat unxptdly, w obsrv xponntial runtim in n. This is atually a proprty of our random instans, in whih β tnds to inras linarly with n. Not that this xprimnt dos not invalidat th algorithm sin in pratial appliations w may b intrstd mainly in instans for whih β is indd small. W hav also gnratd instans with small solutions as follows. Pik k, n and β, thn sampl a uniformly-random start prmutation and β uniformlyrandom blok rossings. Considr th squn of prmutations rsulting from ths blok rossings. W randomly sampl n mtings by piking, for ah on indpndntly, on of th prmutations at random and thn two adjant haratrs from this prmutation. An instan onstrutd in this way might hav a solution with fwr than β blok rossings, but by onstrution th optimum is at most β. On ths instans, th runtim of th algorithm of Thorm 5 sals as xptd. Sin th onstrution of ths instans is somwhat ompliatd and arbitrarily, w now rturn to th mor natural first typ of random instans. Th xat algorithm is fasibl only for rathr small instans, so w now shift our fous to th grdy algorithm. Rall that it starts with an arbitrary

18 890 T. C. van Dijk t al. Blok Crossings in Storylin Visualizations Frquny Blok movs Figur 9: Histogram of th numbr of blok rossings usd by th grdy algorithm for all k! diffrnt start prmutations, on a singl random instan with n = 100 and k = 8. Th bst grdy solution uss 45 blok rossings; th avrag is 51.2, with standard dviation Frquny Frquny Exss blok movs Exss blok movs Figur 10: Lft: histogram of HuristiGrdy minus BstGrdy, 200 instans with with k = 7 and n = 100. Right: histogram of RandomGrdy minus HuristiGrdy, 1000 instans with k = 30 and n = 200. prmutation and prods grdily. Th histogram in Fig. 9 shows th numbr of blok rossings usd by th grdy algorithm for ah of th k! possibl start prmutations, for a singl fixd instan. This bll urv is typial for instans gnratd by our random modl, and is shown for illustration. W s, qualitativly, that thr indd xist rar start prmutations that do notiably bttr than almost all othrs (hr with 45 blok rossings). Indd, for th rportd instan, a random start prmutation dos 7.2 blok rossings wors in xptation than th bst possibl start prmutation, whil th numbr of blok rossings for diffrnt start prmutations for this instan has standard dviation of 2.4. W all th bst possibl rsult of th grdy algorithm ovr all start prmutations BstGrdy, whih w alulat by brut for. This is dtrministi

19 JGAA, 21(5) (2017) Frquny Frquny Exss blok movs Exss blok movs Figur 11: Lft: Huristi-grdy minus OPT with k = 5 and n = 12, 1000 runs. Right: th sam, xpt k = 6 and 100 runs. givn th instan. Lt RandomGrdy b th rsult of th grdy algorithm starting with a prmutation hosn uniformly at random, and lt Huristi- Grdy b th rsult of starting with th huristi prmutation that w hav dsribd abov. Th histogram in Fig. 10 (lft) shows how many mor blok rossings HuristiGrdy uss than BstGrdy on random instans. This distribution is havist nar zro (maning th huristi is oftn los to optimal), but thr ar instans whr th diffrn is larg. On avrag th huristi is 4.1 blok rossings wors than bst possibl start prmutation on ths instans, whih is an improvmnt ovr random start prmutations; s Fig. 10, right. Not that w do not know how to omput BstGrdy ffiintly, but this xprimnt dmonstrats that thr is room for improvmnt within th grdy algorithm by piking a good start prmutation. Lastly, w ompar th grdy algorithm to th optimum, whih (baus of runtim) w only do for small k and n. On 1000 random instans with k = 5 and n = 12, HuristiGrdy gav an optimal solution 56% of th tim. It ndd somtims on (38%), two (5%), or thr (1%) blok rossings mor, but nvr mor than that. With k = 6, w gt similar numbrs; s Fig. 11. This is a promising bhavior, but larly annot b xtrapolatd vrbatim to largr instans. Basd on ths xprimnts, w rommnd HuristiGrdy as an ffiint, rasonabl huristi for th 2-SBCM problm whn it ariss. W hav not run xprimnts to visualiz ral storylins (movis, books, t tra) sin th rstrition to ardinality-2 mtings is too svr. Th gnalogial visualizations of Kim t al. [8] do fit this rstrition, but rquir furthr visualization work on top of solving th SBCM instan, whih w hav not prformd. Instad, w hav fousd ths xprimnts on som ombinatorial proprtis of random instans.

20 892 T. C. van Dijk t al. Blok Crossings in Storylin Visualizations Figur 12: Mting {, 1, 2, 3 } 6 Approximation Algorithm W now dvlop a onstant-fator approximation algorithm for d-sbcm whr d is a onstant. W initially assum that ah group mting ours xatly on, but latr show how to xtnd our rsults to th stting whr th sam group an mt a boundd numbr of tims. Ovrviw. Our approximation algorithm has th following thr main stps. 1. Rdu th input group hyprgraph H = (C, Γ) to an intrval hyprgraph H f = (C, Γ \ Γ p ) by dlting a subst Γ p Γ of th dgs of H. 2. Choos a prmutation π 0 of th haratrs that supports all groups of this intrval hyprgraph H f. Thus, π 0 is th ordr of haratrs at th bginning of th timlin. 3. Inrmntally rat support for ah dltd mting of Γ p in ordr of inrasing tim, as follows. W arbitrarily fix on haratr of th mting and mov th othr haratrs of th mting nxt to it. Aftr th mting w mov th haratrs bak to thir original positions; s Fig. 12. Stp 2 is straightforward: Stion 2 shows how to find a prmutation supporting all th groups for an intrval hyprgraph. Th main thnial parts of th algorithm ar Stp 1 and an analysis to harg at most a onstant numbr of blok rossings in Stp 3 to a blok rossing in th optimal visualization. Stp 1 rquirs solving a hyprgraph problm; this is thnially th most hallnging part, and onsums th ntir Stion 7. Bounds and Analysis. W all th dgs in Γ p paid dgs, and th dgs in Γ f = Γ \ Γ p fr dgs. Intuitivly, a fr dg an b ralizd without a blok rossing baus H f is an intrval hyprgraph, whil ah dg Γ p must b hargd to blok rossings of th optimal drawing. W initializ th drawing by plaing th haratrs in th vrtial ordr π 0, whih supports all th groups in Γ f. Now w onsidr th paid dgs in lft-to-right ordr. Suppos that th nxt mting involvs a group g Γ p. W hav g d and fix on haratr of g. To bring th rmaining haratrs of g in this haratr s viinity, w nd at most (d 1) blok rossings, on pr lin. Whn th mting is ovr, w

21 JGAA, 21(5) (2017) 893 again us up to (d 1) blok rossings to rvrt th lins bak to thir original position prsribd by π 0 ; s Fig. 12. W do this for ah paid hyprdg, giving ris to at most 2(d 1) Γ p blok rossings. W now prov that this bound is within a onstant fator of optimal. W first stablish a lowr bound on th optimal numbr of blok rossings assuming a fixd start prmutation. Lmma 4 Lt π b a prmutation of th haratrs, lt Γ f b th st of groups supportd by π, and lt Γ p = Γ \ Γ f. Any storylin visualization with start prmutation π nds at last 2 Γ p /(3d 2 ) blok rossings. Proof: Lt g Γ p. Sin g is not supportd by π, th optimal drawing dos not ontain th haratrs of g as a ontiguous blok initially. Howvr, in ordr to support this mting, ths haratrs must vntually bom ontiguous bfor th mting starts. Th ordr hangs only through (blok) rossings; w bound th numbr of groups that an bom supportd aftr ah blok rossing. Aftr a blok rossing, at most thr pairs of lins that wr not nighbors bfor an bom nighbors in th prmutation: aftr th bloks C 1, C 2 C ross, thr is on position in th prmutation whr a lin of C 1 is nxt to a lin of C 2, and two positions with a lin of C 1 (C 2, rsptivly) and a lin of C \ (C 1 C 2 ). Any group that was not supportd, but is supportd aftr th blok rossing, must ontain on of ths pairs. W an dsrib ah suh group in th nw prmutation by spifying th nw pair and th numbrs d 1 and d 2 of haratrs of th group abov and blow th nw pair in th prmutation. Sin th group siz is at most d, w hav d 1 + d 2 d. For any d 1 1, thr ar d d 1 possibl hois for d 2. Togthr with d 1, d 2 1 (sin th nw pair is inludd), w gt at most d 1 d (d d 1=1 1) = d 1 d d 1=1 1 = d(d 1)/2 d 2 /2 possibl groups for ah nw pair. Thus, th total numbr of nwly supportd groups aftr a blok rossing is at most 3d 2 /2, whih shows that th optimal numbr of blok rossings is at last 2 Γ p /(3d 2 ), omplting th proof. W now bound th loss of optimality ausd by not knowing th initial prmutation usd by th optimal solution. Th ky ida hr is to us a onstantfator approximation for th problm Intrval Hyprgraph Edg Dltion: dlt th minimum numbr of hyprdgs from H so that H boms an intrval hyprgraph. W prov th following thorm in Stion 7. Thorm 8 W an find a (d + 1)-approximation for Intrval Hyprgraph Edg Dltion on group hyprgraphs with m hyprdgs of rank d in O(m 2 ) tim. With th hlp of Thorm 8, w an show th main rsult of this papr. Thorm 9 d-sbcm admits a (3d 2 (d 2 1))-approximation algorithm that runs in O(kn) tim.

22 894 T. C. van Dijk t al. Blok Crossings in Storylin Visualizations Proof: Lt Γ p b th st of paid dgs in our algorithm, and lt Γ OPT b th st of paid dgs in th optimal solution, that is, th st of dgs that ar not supportd by th initial prmutation of th optimal solution. By Thorm 8, w hav Γ p (d + 1) Γ OPT, sin vn Γ OPT annot b smallr than an optimum solution for th intrval hyprgraph dg dltion problm. Lt ALG and OPT b th numbrs of blok rossings for our algorithm and th optimal solution, rsptivly. Our algorithm nsurs that ALG 2(d 1) Γ p 2(d 1)(d + 1) Γ OPT. On th othr hand, by Lmma 4, w hav OPT 2 Γ OPT /(3d 2 ), whih givs Γ OPT 3d 2 /2 OPT. Combining th prvious inqualitis, w gt ALG 3d 2 (d 2 1) OPT as dsird. Now w analyz th tim omplxity. W hav to onsidr th prmutation (of lngth k) of haratrs bfor and aftr ah of th n mtings, as wll as aftr ah of th O(n) blok rossings. This rsults in O(kn) tim for th last part of th algorithm, but this is dominatd by th tim (O(n 2 )) ndd for finding Γ p and for dtrmining th start prmutation. W an improv th running tim to O(kn) by a slight modifiation: using th approximation algorithm for Intrval Hyprgraph Edg Dltion is only nssary for spars instans. If H has suffiintly many dgs, any start prmutation will yild a good approximation. Sin no mting involvs mor than d haratrs, no start prmutation an support mor than dk mtings. If n 2dk, thn vn th optimal solution must thrfor rmov at last half of th dgs. Hn, taking an arbitrary start prmutation yilds an approximation fator of at most 2 < d + 1. W now hang th algorithm to us an arbitrary start prmutation if n 2dk and only us th approximation for Intrval Hyprgraph Edg Dltion othrwis. In partiular, w us th approximation only if thr ar n O(k) dgs. Hn, for spars instans w hav O(n 2 ) O(kn), and for dns instans, w skip th ostly O(n 2 )-tim stp mntiond abov. This yilds th dsird running tim of O(kn), whih is worst-as optimal sin th output omplxity is of th sam ordr. Rmark. W assumd that ah group mts only on, but w an xtnd th rsult if ah group an mt α tims, for onstant α. Our algorithm thn yilds a (α 3d 2 (d 2 1))-fator approximation; ah rptition of a mting may triggr a onstant numbr of blok rossings not prsnt in th optimal solution. Improvmnts for 2-SBCM. By using spifi struturs for 2-haratr mtings w an improv approximation fator and runtim; not that th gnral algorithm yilds a 36-approximation. For 2-haratr mtings th group hyprgraph is a graph; an intrval hyprgraph in this stting is a olltion of vrtx-disjoint paths. Our algorithm for Intrval Hyprgraph Edg Dltion for d = 2 yilds a 3-approximation. W dvlop a bttr approximation using th following obsrvation. Considr a haratr in th olltion of paths supportd in th bginning of som solution. If has two nighbors 1 and 2 in its path, but s first mting is with

Linked-List Implementation. Linked-lists for two sets. Multiple Operations. UNION Implementation. An Application of Disjoint-Set 1/9/2014

Linked-List Implementation. Linked-lists for two sets. Multiple Operations. UNION Implementation. An Application of Disjoint-Set 1/9/2014 Disjoint Sts Data Strutur (Chap. 21) A disjoint-st is a olltion ={S 1, S 2,, S k } o distint dynami sts. Eah st is idntiid by a mmbr o th st, alld rprsntativ. Disjoint st oprations: MAKE-SET(x): rat a

More information

AP Calculus BC Problem Drill 16: Indeterminate Forms, L Hopital s Rule, & Improper Intergals

AP Calculus BC Problem Drill 16: Indeterminate Forms, L Hopital s Rule, & Improper Intergals AP Calulus BC Problm Drill 6: Indtrminat Forms, L Hopital s Rul, & Impropr Intrgals Qustion No. of Instrutions: () Rad th problm and answr hois arfully () Work th problms on papr as ndd () Pik th answr

More information

Utilizing exact and Monte Carlo methods to investigate properties of the Blume Capel Model applied to a nine site lattice.

Utilizing exact and Monte Carlo methods to investigate properties of the Blume Capel Model applied to a nine site lattice. Utilizing xat and Mont Carlo mthods to invstigat proprtis of th Blum Capl Modl applid to a nin sit latti Nik Franios Writing various xat and Mont Carlo omputr algorithms in C languag, I usd th Blum Capl

More information

1 Minimum Cut Problem

1 Minimum Cut Problem CS 6 Lctur 6 Min Cut and argr s Algorithm Scribs: Png Hui How (05), Virginia Dat: May 4, 06 Minimum Cut Problm Today, w introduc th minimum cut problm. This problm has many motivations, on of which coms

More information

cycle that does not cross any edges (including its own), then it has at least

cycle that does not cross any edges (including its own), then it has at least W prov th following thorm: Thorm If a K n is drawn in th plan in such a way that it has a hamiltonian cycl that dos not cross any dgs (including its own, thn it has at last n ( 4 48 π + O(n crossings Th

More information

Lecture 14 (Oct. 30, 2017)

Lecture 14 (Oct. 30, 2017) Ltur 14 8.31 Quantum Thory I, Fall 017 69 Ltur 14 (Ot. 30, 017) 14.1 Magnti Monopols Last tim, w onsidrd a magnti fild with a magnti monopol onfiguration, and bgan to approah dsribing th quantum mhanis

More information

The Matrix Exponential

The Matrix Exponential Th Matrix Exponntial (with xrciss) by D. Klain Vrsion 207.0.05 Corrctions and commnts ar wlcom. Th Matrix Exponntial For ach n n complx matrix A, dfin th xponntial of A to b th matrix A A k I + A + k!

More information

Junction Tree Algorithm 1. David Barber

Junction Tree Algorithm 1. David Barber Juntion Tr Algorithm 1 David Barbr Univrsity Collg London 1 Ths slids aompany th book Baysian Rasoning and Mahin Larning. Th book and dmos an b downloadd from www.s.ul.a.uk/staff/d.barbr/brml. Fdbak and

More information

u x v x dx u x v x v x u x dx d u x v x u x v x dx u x v x dx Integration by Parts Formula

u x v x dx u x v x v x u x dx d u x v x u x v x dx u x v x dx Integration by Parts Formula 7. Intgration by Parts Each drivativ formula givs ris to a corrsponding intgral formula, as w v sn many tims. Th drivativ product rul yilds a vry usful intgration tchniqu calld intgration by parts. Starting

More information

Higher order derivatives

Higher order derivatives Robrto s Nots on Diffrntial Calculus Chaptr 4: Basic diffrntiation ruls Sction 7 Highr ordr drivativs What you nd to know alrady: Basic diffrntiation ruls. What you can larn hr: How to rpat th procss of

More information

The Matrix Exponential

The Matrix Exponential Th Matrix Exponntial (with xrciss) by Dan Klain Vrsion 28928 Corrctions and commnts ar wlcom Th Matrix Exponntial For ach n n complx matrix A, dfin th xponntial of A to b th matrix () A A k I + A + k!

More information

CPSC 665 : An Algorithmist s Toolkit Lecture 4 : 21 Jan Linear Programming

CPSC 665 : An Algorithmist s Toolkit Lecture 4 : 21 Jan Linear Programming CPSC 665 : An Algorithmist s Toolkit Lctur 4 : 21 Jan 2015 Lcturr: Sushant Sachdva Linar Programming Scrib: Rasmus Kyng 1. Introduction An optimization problm rquirs us to find th minimum or maximum) of

More information

Searching Linked Lists. Perfect Skip List. Building a Skip List. Skip List Analysis (1) Assume the list is sorted, but is stored in a linked list.

Searching Linked Lists. Perfect Skip List. Building a Skip List. Skip List Analysis (1) Assume the list is sorted, but is stored in a linked list. 3 3 4 8 6 3 3 4 8 6 3 3 4 8 6 () (d) 3 Sarching Linkd Lists Sarching Linkd Lists Sarching Linkd Lists ssum th list is sortd, but is stord in a linkd list. an w us binary sarch? omparisons? Work? What if

More information

Homework #3. 1 x. dx. It therefore follows that a sum of the

Homework #3. 1 x. dx. It therefore follows that a sum of the Danil Cannon CS 62 / Luan March 5, 2009 Homwork # 1. Th natural logarithm is dfind by ln n = n 1 dx. It thrfor follows that a sum of th 1 x sam addnd ovr th sam intrval should b both asymptotically uppr-

More information

Section 6.1. Question: 2. Let H be a subgroup of a group G. Then H operates on G by left multiplication. Describe the orbits for this operation.

Section 6.1. Question: 2. Let H be a subgroup of a group G. Then H operates on G by left multiplication. Describe the orbits for this operation. MAT 444 H Barclo Spring 004 Homwork 6 Solutions Sction 6 Lt H b a subgroup of a group G Thn H oprats on G by lft multiplication Dscrib th orbits for this opration Th orbits of G ar th right costs of H

More information

Basic Polyhedral theory

Basic Polyhedral theory Basic Polyhdral thory Th st P = { A b} is calld a polyhdron. Lmma 1. Eithr th systm A = b, b 0, 0 has a solution or thr is a vctorπ such that π A 0, πb < 0 Thr cass, if solution in top row dos not ist

More information

CS 361 Meeting 12 10/3/18

CS 361 Meeting 12 10/3/18 CS 36 Mting 2 /3/8 Announcmnts. Homwork 4 is du Friday. If Friday is Mountain Day, homwork should b turnd in at my offic or th dpartmnt offic bfor 4. 2. Homwork 5 will b availabl ovr th wknd. 3. Our midtrm

More information

Addition of angular momentum

Addition of angular momentum Addition of angular momntum April, 0 Oftn w nd to combin diffrnt sourcs of angular momntum to charactriz th total angular momntum of a systm, or to divid th total angular momntum into parts to valuat th

More information

A Propagating Wave Packet Group Velocity Dispersion

A Propagating Wave Packet Group Velocity Dispersion Lctur 8 Phys 375 A Propagating Wav Packt Group Vlocity Disprsion Ovrviw and Motivation: In th last lctur w lookd at a localizd solution t) to th 1D fr-particl Schrödingr quation (SE) that corrsponds to

More information

Introduction to Arithmetic Geometry Fall 2013 Lecture #20 11/14/2013

Introduction to Arithmetic Geometry Fall 2013 Lecture #20 11/14/2013 18.782 Introduction to Arithmtic Gomtry Fall 2013 Lctur #20 11/14/2013 20.1 Dgr thorm for morphisms of curvs Lt us rstat th thorm givn at th nd of th last lctur, which w will now prov. Thorm 20.1. Lt φ:

More information

First derivative analysis

First derivative analysis Robrto s Nots on Dirntial Calculus Chaptr 8: Graphical analysis Sction First drivativ analysis What you nd to know alrady: How to us drivativs to idntiy th critical valus o a unction and its trm points

More information

ECE602 Exam 1 April 5, You must show ALL of your work for full credit.

ECE602 Exam 1 April 5, You must show ALL of your work for full credit. ECE62 Exam April 5, 27 Nam: Solution Scor: / This xam is closd-book. You must show ALL of your work for full crdit. Plas rad th qustions carfully. Plas chck your answrs carfully. Calculators may NOT b

More information

Assignment 4 Biophys 4322/5322

Assignment 4 Biophys 4322/5322 Assignmnt 4 Biophys 4322/5322 Tylr Shndruk Fbruary 28, 202 Problm Phillips 7.3. Part a R-onsidr dimoglobin utilizing th anonial nsmbl maning rdriv Eq. 3 from Phillips Chaptr 7. For a anonial nsmbl p E

More information

(Upside-Down o Direct Rotation) β - Numbers

(Upside-Down o Direct Rotation) β - Numbers Amrican Journal of Mathmatics and Statistics 014, 4(): 58-64 DOI: 10593/jajms0140400 (Upsid-Down o Dirct Rotation) β - Numbrs Ammar Sddiq Mahmood 1, Shukriyah Sabir Ali,* 1 Dpartmnt of Mathmatics, Collg

More information

Construction of asymmetric orthogonal arrays of strength three via a replacement method

Construction of asymmetric orthogonal arrays of strength three via a replacement method isid/ms/26/2 Fbruary, 26 http://www.isid.ac.in/ statmath/indx.php?modul=prprint Construction of asymmtric orthogonal arrays of strngth thr via a rplacmnt mthod Tian-fang Zhang, Qiaoling Dng and Alok Dy

More information

Search sequence databases 3 10/25/2016

Search sequence databases 3 10/25/2016 Sarch squnc databass 3 10/25/2016 Etrm valu distribution Ø Suppos X is a random variabl with probability dnsity function p(, w sampl a larg numbr S of indpndnt valus of X from this distribution for an

More information

CS553 Lecture Register Allocation I 3

CS553 Lecture Register Allocation I 3 Low-Lvl Issus Last ltur Intrproural analysis Toay Start low-lvl issus Rgistr alloation Latr Mor rgistr alloation Instrution shuling CS553 Ltur Rgistr Alloation I 2 Rgistr Alloation Prolm Assign an unoun

More information

Addition of angular momentum

Addition of angular momentum Addition of angular momntum April, 07 Oftn w nd to combin diffrnt sourcs of angular momntum to charactriz th total angular momntum of a systm, or to divid th total angular momntum into parts to valuat

More information

Supplementary Materials

Supplementary Materials 6 Supplmntary Matrials APPENDIX A PHYSICAL INTERPRETATION OF FUEL-RATE-SPEED FUNCTION A truck running on a road with grad/slop θ positiv if moving up and ngativ if moving down facs thr rsistancs: arodynamic

More information

SECTION where P (cos θ, sin θ) and Q(cos θ, sin θ) are polynomials in cos θ and sin θ, provided Q is never equal to zero.

SECTION where P (cos θ, sin θ) and Q(cos θ, sin θ) are polynomials in cos θ and sin θ, provided Q is never equal to zero. SETION 6. 57 6. Evaluation of Dfinit Intgrals Exampl 6.6 W hav usd dfinit intgrals to valuat contour intgrals. It may com as a surpris to larn that contour intgrals and rsidus can b usd to valuat crtain

More information

On spanning trees and cycles of multicolored point sets with few intersections

On spanning trees and cycles of multicolored point sets with few intersections On spanning trs and cycls of multicolord point sts with fw intrsctions M. Kano, C. Mrino, and J. Urrutia April, 00 Abstract Lt P 1,..., P k b a collction of disjoint point sts in R in gnral position. W

More information

Derangements and Applications

Derangements and Applications 2 3 47 6 23 Journal of Intgr Squncs, Vol. 6 (2003), Articl 03..2 Drangmnts and Applications Mhdi Hassani Dpartmnt of Mathmatics Institut for Advancd Studis in Basic Scincs Zanjan, Iran mhassani@iasbs.ac.ir

More information

Abstract Interpretation: concrete and abstract semantics

Abstract Interpretation: concrete and abstract semantics Abstract Intrprtation: concrt and abstract smantics Concrt smantics W considr a vry tiny languag that manags arithmtic oprations on intgrs valus. Th (concrt) smantics of th languags cab b dfind by th funzcion

More information

EEO 401 Digital Signal Processing Prof. Mark Fowler

EEO 401 Digital Signal Processing Prof. Mark Fowler EEO 401 Digital Signal Procssing Prof. Mark Fowlr Dtails of th ot St #19 Rading Assignmnt: Sct. 7.1.2, 7.1.3, & 7.2 of Proakis & Manolakis Dfinition of th So Givn signal data points x[n] for n = 0,, -1

More information

LINEAR DELAY DIFFERENTIAL EQUATION WITH A POSITIVE AND A NEGATIVE TERM

LINEAR DELAY DIFFERENTIAL EQUATION WITH A POSITIVE AND A NEGATIVE TERM Elctronic Journal of Diffrntial Equations, Vol. 2003(2003), No. 92, pp. 1 6. ISSN: 1072-6691. URL: http://jd.math.swt.du or http://jd.math.unt.du ftp jd.math.swt.du (login: ftp) LINEAR DELAY DIFFERENTIAL

More information

Week 3: Connected Subgraphs

Week 3: Connected Subgraphs Wk 3: Connctd Subgraphs Sptmbr 19, 2016 1 Connctd Graphs Path, Distanc: A path from a vrtx x to a vrtx y in a graph G is rfrrd to an xy-path. Lt X, Y V (G). An (X, Y )-path is an xy-path with x X and y

More information

MA1506 Tutorial 2 Solutions. Question 1. (1a) 1 ) y x. e x. 1 exp (in general, Integrating factor is. ye dx. So ) (1b) e e. e c.

MA1506 Tutorial 2 Solutions. Question 1. (1a) 1 ) y x. e x. 1 exp (in general, Integrating factor is. ye dx. So ) (1b) e e. e c. MA56 utorial Solutions Qustion a Intgrating fator is ln p p in gnral, multipl b p So b ln p p sin his kin is all a Brnoulli quation -- st Sin w fin Y, Y Y, Y Y p Qustion Dfin v / hn our quation is v μ

More information

Limiting value of higher Mahler measure

Limiting value of higher Mahler measure Limiting valu of highr Mahlr masur Arunabha Biswas a, Chris Monico a, a Dpartmnt of Mathmatics & Statistics, Txas Tch Univrsity, Lubbock, TX 7949, USA Abstract W considr th k-highr Mahlr masur m k P )

More information

Fourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation.

Fourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation. Lur 7 Fourir Transforms and th Wav Euation Ovrviw and Motivation: W first discuss a fw faturs of th Fourir transform (FT), and thn w solv th initial-valu problm for th wav uation using th Fourir transform

More information

perm4 A cnt 0 for for if A i 1 A i cnt cnt 1 cnt i j. j k. k l. i k. j l. i l

perm4 A cnt 0 for for if A i 1 A i cnt cnt 1 cnt i j. j k. k l. i k. j l. i l h 4D, 4th Rank, Antisytric nsor and th 4D Equivalnt to th Cross Product or Mor Fun with nsors!!! Richard R Shiffan Digital Graphics Assoc 8 Dunkirk Av LA, Ca 95 rrs@isidu his docunt dscribs th four dinsional

More information

Codes on Planar Graphs

Codes on Planar Graphs Cods on Planar Graphs Srimathy Srinivasan, Andrw Thangaraj Dpartmnt of Eltrial Enginring Indian Institut of Thnology, Madras Email: andrw@iitm.a.in arxiv:090.0768v [s.it] 5 May 009 Astrat Cods dfind on

More information

Propositional Logic. Combinatorial Problem Solving (CPS) Albert Oliveras Enric Rodríguez-Carbonell. May 17, 2018

Propositional Logic. Combinatorial Problem Solving (CPS) Albert Oliveras Enric Rodríguez-Carbonell. May 17, 2018 Propositional Logic Combinatorial Problm Solving (CPS) Albrt Olivras Enric Rodríguz-Carbonll May 17, 2018 Ovrviw of th sssion Dfinition of Propositional Logic Gnral Concpts in Logic Rduction to SAT CNFs

More information

Mor Tutorial at www.dumblittldoctor.com Work th problms without a calculator, but us a calculator to chck rsults. And try diffrntiating your answrs in part III as a usful chck. I. Applications of Intgration

More information

EXST Regression Techniques Page 1

EXST Regression Techniques Page 1 EXST704 - Rgrssion Tchniqus Pag 1 Masurmnt rrors in X W hav assumd that all variation is in Y. Masurmnt rror in this variabl will not ffct th rsults, as long as thy ar uncorrlatd and unbiasd, sinc thy

More information

Calculus concepts derivatives

Calculus concepts derivatives All rasonabl fforts hav bn mad to mak sur th nots ar accurat. Th author cannot b hld rsponsibl for any damags arising from th us of ths nots in any fashion. Calculus concpts drivativs Concpts involving

More information

Brief Introduction to Statistical Mechanics

Brief Introduction to Statistical Mechanics Brif Introduction to Statistical Mchanics. Purpos: Ths nots ar intndd to provid a vry quick introduction to Statistical Mchanics. Th fild is of cours far mor vast than could b containd in ths fw pags.

More information

Integral Calculus What is integral calculus?

Integral Calculus What is integral calculus? Intgral Calulus What is intgral alulus? In diffrntial alulus w diffrntiat a funtion to obtain anothr funtion alld drivativ. Intgral alulus is onrnd with th opposit pross. Rvrsing th pross of diffrntiation

More information

First order differential equation Linear equation; Method of integrating factors

First order differential equation Linear equation; Method of integrating factors First orr iffrntial quation Linar quation; Mtho of intgrating factors Exampl 1: Rwrit th lft han si as th rivativ of th prouct of y an som function by prouct rul irctly. Solving th first orr iffrntial

More information

Calculus II (MAC )

Calculus II (MAC ) Calculus II (MAC232-2) Tst 2 (25/6/25) Nam (PRINT): Plas show your work. An answr with no work rcivs no crdit. You may us th back of a pag if you nd mor spac for a problm. You may not us any calculators.

More information

Answer Homework 5 PHA5127 Fall 1999 Jeff Stark

Answer Homework 5 PHA5127 Fall 1999 Jeff Stark Answr omwork 5 PA527 Fall 999 Jff Stark A patint is bing tratd with Drug X in a clinical stting. Upon admiion, an IV bolus dos of 000mg was givn which yildd an initial concntration of 5.56 µg/ml. A fw

More information

Finding low cost TSP and 2-matching solutions using certain half integer subtour vertices

Finding low cost TSP and 2-matching solutions using certain half integer subtour vertices Finding low cost TSP and 2-matching solutions using crtain half intgr subtour vrtics Sylvia Boyd and Robrt Carr Novmbr 996 Introduction Givn th complt graph K n = (V, E) on n nods with dg costs c R E,

More information

u 3 = u 3 (x 1, x 2, x 3 )

u 3 = u 3 (x 1, x 2, x 3 ) Lctur 23: Curvilinar Coordinats (RHB 8.0 It is oftn convnint to work with variabls othr than th Cartsian coordinats x i ( = x, y, z. For xampl in Lctur 5 w mt sphrical polar and cylindrical polar coordinats.

More information

That is, we start with a general matrix: And end with a simpler matrix:

That is, we start with a general matrix: And end with a simpler matrix: DIAGON ALIZATION OF THE STR ESS TEN SOR INTRO DUCTIO N By th us of Cauchy s thorm w ar abl to rduc th numbr of strss componnts in th strss tnsor to only nin valus. An additional simplification of th strss

More information

Problem Set 6 Solutions

Problem Set 6 Solutions 6.04/18.06J Mathmatics for Computr Scinc March 15, 005 Srini Dvadas and Eric Lhman Problm St 6 Solutions Du: Monday, March 8 at 9 PM in Room 3-044 Problm 1. Sammy th Shark is a financial srvic providr

More information

Steinberg s Conjecture is false

Steinberg s Conjecture is false Stinrg s Conjtur is als arxiv:1604.05108v2 [math.co] 19 Apr 2016 Vinnt Cohn-Aa Mihal Hig Danil Král Zhntao Li Estan Salgao Astrat Stinrg onjtur in 1976 that vry planar graph with no yls o lngth our or

More information

Chapter 6 Folding. Folding

Chapter 6 Folding. Folding Chaptr 6 Folding Wintr 1 Mokhtar Abolaz Folding Th folding transformation is usd to systmatically dtrmin th control circuits in DSP architctur whr multipl algorithm oprations ar tim-multiplxd to a singl

More information

ANALYSIS IN THE FREQUENCY DOMAIN

ANALYSIS IN THE FREQUENCY DOMAIN ANALYSIS IN THE FREQUENCY DOMAIN SPECTRAL DENSITY Dfinition Th spctral dnsit of a S.S.P. t also calld th spctrum of t is dfind as: + { γ }. jτ γ τ F τ τ In othr words, of th covarianc function. is dfind

More information

COMPUTER GENERATED HOLOGRAMS Optical Sciences 627 W.J. Dallas (Monday, April 04, 2005, 8:35 AM) PART I: CHAPTER TWO COMB MATH.

COMPUTER GENERATED HOLOGRAMS Optical Sciences 627 W.J. Dallas (Monday, April 04, 2005, 8:35 AM) PART I: CHAPTER TWO COMB MATH. C:\Dallas\0_Courss\03A_OpSci_67\0 Cgh_Book\0_athmaticalPrliminaris\0_0 Combath.doc of 8 COPUTER GENERATED HOLOGRAS Optical Scincs 67 W.J. Dallas (onday, April 04, 005, 8:35 A) PART I: CHAPTER TWO COB ATH

More information

u r du = ur+1 r + 1 du = ln u + C u sin u du = cos u + C cos u du = sin u + C sec u tan u du = sec u + C e u du = e u + C

u r du = ur+1 r + 1 du = ln u + C u sin u du = cos u + C cos u du = sin u + C sec u tan u du = sec u + C e u du = e u + C Tchniqus of Intgration c Donald Kridr and Dwight Lahr In this sction w ar going to introduc th first approachs to valuating an indfinit intgral whos intgrand dos not hav an immdiat antidrivativ. W bgin

More information

The Equitable Dominating Graph

The Equitable Dominating Graph Intrnational Journal of Enginring Rsarch and Tchnology. ISSN 0974-3154 Volum 8, Numbr 1 (015), pp. 35-4 Intrnational Rsarch Publication Hous http://www.irphous.com Th Equitabl Dominating Graph P.N. Vinay

More information

Function Spaces. a x 3. (Letting x = 1 =)) a(0) + b + c (1) = 0. Row reducing the matrix. b 1. e 4 3. e 9. >: (x = 1 =)) a(0) + b + c (1) = 0

Function Spaces. a x 3. (Letting x = 1 =)) a(0) + b + c (1) = 0. Row reducing the matrix. b 1. e 4 3. e 9. >: (x = 1 =)) a(0) + b + c (1) = 0 unction Spacs Prrquisit: Sction 4.7, Coordinatization n this sction, w apply th tchniqus of Chaptr 4 to vctor spacs whos lmnts ar functions. Th vctor spacs P n and P ar familiar xampls of such spacs. Othr

More information

SCHUR S THEOREM REU SUMMER 2005

SCHUR S THEOREM REU SUMMER 2005 SCHUR S THEOREM REU SUMMER 2005 1. Combinatorial aroach Prhas th first rsult in th subjct blongs to I. Schur and dats back to 1916. On of his motivation was to study th local vrsion of th famous quation

More information

Quasi-Classical States of the Simple Harmonic Oscillator

Quasi-Classical States of the Simple Harmonic Oscillator Quasi-Classical Stats of th Simpl Harmonic Oscillator (Draft Vrsion) Introduction: Why Look for Eignstats of th Annihilation Oprator? Excpt for th ground stat, th corrspondnc btwn th quantum nrgy ignstats

More information

Homotopy perturbation technique

Homotopy perturbation technique Comput. Mthods Appl. Mch. Engrg. 178 (1999) 257±262 www.lsvir.com/locat/cma Homotopy prturbation tchniqu Ji-Huan H 1 Shanghai Univrsity, Shanghai Institut of Applid Mathmatics and Mchanics, Shanghai 272,

More information

Abstract Interpretation. Lecture 5. Profs. Aiken, Barrett & Dill CS 357 Lecture 5 1

Abstract Interpretation. Lecture 5. Profs. Aiken, Barrett & Dill CS 357 Lecture 5 1 Abstract Intrprtation 1 History On brakthrough papr Cousot & Cousot 77 (?) Inspird by Dataflow analysis Dnotational smantics Enthusiastically mbracd by th community At last th functional community... At

More information

COUNTING TAMELY RAMIFIED EXTENSIONS OF LOCAL FIELDS UP TO ISOMORPHISM

COUNTING TAMELY RAMIFIED EXTENSIONS OF LOCAL FIELDS UP TO ISOMORPHISM COUNTING TAMELY RAMIFIED EXTENSIONS OF LOCAL FIELDS UP TO ISOMORPHISM Jim Brown Dpartmnt of Mathmatical Scincs, Clmson Univrsity, Clmson, SC 9634, USA jimlb@g.clmson.du Robrt Cass Dpartmnt of Mathmatics,

More information

Note If the candidate believes that e x = 0 solves to x = 0 or gives an extra solution of x = 0, then withhold the final accuracy mark.

Note If the candidate believes that e x = 0 solves to x = 0 or gives an extra solution of x = 0, then withhold the final accuracy mark. . (a) Eithr y = or ( 0, ) (b) Whn =, y = ( 0 + ) = 0 = 0 ( + ) = 0 ( )( ) = 0 Eithr = (for possibly abov) or = A 3. Not If th candidat blivs that = 0 solvs to = 0 or givs an tra solution of = 0, thn withhold

More information

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches.

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches. Subjct Chmistry Papr No and Titl Modul No and Titl Modul Tag 8/ Physical Spctroscopy / Brakdown of th Born-Oppnhimr approximation. Slction ruls for rotational-vibrational transitions. P, R branchs. CHE_P8_M

More information

Unconstrained Minimization of Quadratic Functions via Min-sum

Unconstrained Minimization of Quadratic Functions via Min-sum Unonstraind Minimization of Quadrati Funtions via Min-sum Niholas Ruozzi Computr Sin Yal Univrsity Nw Havn, CT 06520, USA Email: niholas.ruozzi@yal.du Skhar Tatikonda Eltrial Enginring Yal Univrsity Nw

More information

Strongly Connected Components

Strongly Connected Components Strongly Connctd Componnts Lt G = (V, E) b a dirctd graph Writ if thr is a path from to in G Writ if and is an quivalnc rlation: implis and implis s quivalnc classs ar calld th strongly connctd componnts

More information

INCOMPLETE KLOOSTERMAN SUMS AND MULTIPLICATIVE INVERSES IN SHORT INTERVALS. xy 1 (mod p), (x, y) I (j)

INCOMPLETE KLOOSTERMAN SUMS AND MULTIPLICATIVE INVERSES IN SHORT INTERVALS. xy 1 (mod p), (x, y) I (j) INCOMPLETE KLOOSTERMAN SUMS AND MULTIPLICATIVE INVERSES IN SHORT INTERVALS T D BROWNING AND A HAYNES Abstract W invstigat th solubility of th congrunc xy (mod ), whr is a rim and x, y ar rstrictd to li

More information

Final Exam Solutions

Final Exam Solutions CS 2 Advancd Data Structurs and Algorithms Final Exam Solutions Jonathan Turnr /8/20. (0 points) Suppos that r is a root of som tr in a Fionacci hap. Assum that just for a dltmin opration, r has no childrn

More information

Network Congestion Games

Network Congestion Games Ntwork Congstion Gams Assistant Profssor Tas A&M Univrsity Collg Station, TX TX Dallas Collg Station Austin Houston Bst rout dpnds on othrs Ntwork Congstion Gams Travl tim incrass with congstion Highway

More information

Computing and Communications -- Network Coding

Computing and Communications -- Network Coding 89 90 98 00 Computing and Communications -- Ntwork Coding Dr. Zhiyong Chn Institut of Wirlss Communications Tchnology Shanghai Jiao Tong Univrsity China Lctur 5- Nov. 05 0 Classical Information Thory Sourc

More information

Thus, because if either [G : H] or [H : K] is infinite, then [G : K] is infinite, then [G : K] = [G : H][H : K] for all infinite cases.

Thus, because if either [G : H] or [H : K] is infinite, then [G : K] is infinite, then [G : K] = [G : H][H : K] for all infinite cases. Homwork 5 M 373K Solutions Mark Lindbrg and Travis Schdlr 1. Prov that th ring Z/mZ (for m 0) is a fild if and only if m is prim. ( ) Proof by Contrapositiv: Hr, thr ar thr cass for m not prim. m 0: Whn

More information

Physics 506 Winter 2006 Homework Assignment #12 Solutions. Textbook problems: Ch. 14: 14.2, 14.4, 14.6, 14.12

Physics 506 Winter 2006 Homework Assignment #12 Solutions. Textbook problems: Ch. 14: 14.2, 14.4, 14.6, 14.12 Physis 56 Wintr 6 Homwork Assignmnt # Solutions Ttbook problms: Ch. 4: 4., 4.4, 4.6, 4. 4. A partil of harg is moving in narly uniform nonrlativisti motion. For tims nar t = t, its vtorial position an

More information

Data Assimilation 1. Alan O Neill National Centre for Earth Observation UK

Data Assimilation 1. Alan O Neill National Centre for Earth Observation UK Data Assimilation 1 Alan O Nill National Cntr for Earth Obsrvation UK Plan Motivation & basic idas Univariat (scalar) data assimilation Multivariat (vctor) data assimilation 3d-Variational Mthod (& optimal

More information

Differential Equations

Differential Equations Prfac Hr ar m onlin nots for m diffrntial quations cours that I tach hr at Lamar Univrsit. Dspit th fact that ths ar m class nots, th should b accssibl to anon wanting to larn how to solv diffrntial quations

More information

Square of Hamilton cycle in a random graph

Square of Hamilton cycle in a random graph Squar of Hamilton cycl in a random graph Andrzj Dudk Alan Friz Jun 28, 2016 Abstract W show that p = n is a sharp thrshold for th random graph G n,p to contain th squar of a Hamilton cycl. This improvs

More information

Cramér-Rao Inequality: Let f(x; θ) be a probability density function with continuous parameter

Cramér-Rao Inequality: Let f(x; θ) be a probability density function with continuous parameter WHEN THE CRAMÉR-RAO INEQUALITY PROVIDES NO INFORMATION STEVEN J. MILLER Abstract. W invstigat a on-paramtr family of probability dnsitis (rlatd to th Parto distribution, which dscribs many natural phnomna)

More information

1 Isoparametric Concept

1 Isoparametric Concept UNIVERSITY OF CALIFORNIA BERKELEY Dpartmnt of Civil Enginring Spring 06 Structural Enginring, Mchanics and Matrials Profssor: S. Govindj Nots on D isoparamtric lmnts Isoparamtric Concpt Th isoparamtric

More information

Title: Vibrational structure of electronic transition

Title: Vibrational structure of electronic transition Titl: Vibrational structur of lctronic transition Pag- Th band spctrum sn in th Ultra-Violt (UV) and visibl (VIS) rgions of th lctromagntic spctrum can not intrprtd as vibrational and rotational spctrum

More information

Equidistribution and Weyl s criterion

Equidistribution and Weyl s criterion Euidistribution and Wyl s critrion by Brad Hannigan-Daly W introduc th ida of a sunc of numbrs bing uidistributd (mod ), and w stat and prov a thorm of Hrmann Wyl which charactrizs such suncs. W also discuss

More information

Solutions to Homework 5

Solutions to Homework 5 Solutions to Homwork 5 Pro. Silvia Frnánz Disrt Mathmatis Math 53A, Fall 2008. [3.4 #] (a) Thr ar x olor hois or vrtx an x or ah o th othr thr vrtis. So th hromati polynomial is P (G, x) =x (x ) 3. ()

More information

On the irreducibility of some polynomials in two variables

On the irreducibility of some polynomials in two variables ACTA ARITHMETICA LXXXII.3 (1997) On th irrducibility of som polynomials in two variabls by B. Brindza and Á. Pintér (Dbrcn) To th mmory of Paul Erdős Lt f(x) and g(y ) b polynomials with intgral cofficints

More information

Lecture 16: Bipolar Junction Transistors. Large Signal Models.

Lecture 16: Bipolar Junction Transistors. Large Signal Models. Whits, EE 322 Ltur 16 Pag 1 of 8 Ltur 16: Bipolar Juntion Transistors. Larg Signal Modls. Transistors prform ky funtions in most ltroni iruits. This is rtainly tru in RF iruits, inluding th NorCal 40A.

More information

Differentiation of Exponential Functions

Differentiation of Exponential Functions Calculus Modul C Diffrntiation of Eponntial Functions Copyright This publication Th Northrn Albrta Institut of Tchnology 007. All Rights Rsrvd. LAST REVISED March, 009 Introduction to Diffrntiation of

More information

Alpha and beta decay equation practice

Alpha and beta decay equation practice Alpha and bta dcay quation practic Introduction Alpha and bta particls may b rprsntd in quations in svral diffrnt ways. Diffrnt xam boards hav thir own prfrnc. For xampl: Alpha Bta α β alpha bta Dspit

More information

ph People Grade Level: basic Duration: minutes Setting: classroom or field site

ph People Grade Level: basic Duration: minutes Setting: classroom or field site ph Popl Adaptd from: Whr Ar th Frogs? in Projct WET: Curriculum & Activity Guid. Bozman: Th Watrcours and th Council for Environmntal Education, 1995. ph Grad Lvl: basic Duration: 10 15 minuts Stting:

More information

The second condition says that a node α of the tree has exactly n children if the arity of its label is n.

The second condition says that a node α of the tree has exactly n children if the arity of its label is n. CS 6110 S14 Hanout 2 Proof of Conflunc 27 January 2014 In this supplmntary lctur w prov that th λ-calculus is conflunt. This is rsult is u to lonzo Church (1903 1995) an J. arkly Rossr (1907 1989) an is

More information

Recall that by Theorems 10.3 and 10.4 together provide us the estimate o(n2 ), S(q) q 9, q=1

Recall that by Theorems 10.3 and 10.4 together provide us the estimate o(n2 ), S(q) q 9, q=1 Chaptr 11 Th singular sris Rcall that by Thorms 10 and 104 togthr provid us th stimat 9 4 n 2 111 Rn = SnΓ 2 + on2, whr th singular sris Sn was dfind in Chaptr 10 as Sn = q=1 Sq q 9, with Sq = 1 a q gcda,q=1

More information

Deift/Zhou Steepest descent, Part I

Deift/Zhou Steepest descent, Part I Lctur 9 Dift/Zhou Stpst dscnt, Part I W now focus on th cas of orthogonal polynomials for th wight w(x) = NV (x), V (x) = t x2 2 + x4 4. Sinc th wight dpnds on th paramtr N N w will writ π n,n, a n,n,

More information

4037 ADDITIONAL MATHEMATICS

4037 ADDITIONAL MATHEMATICS CAMBRIDGE INTERNATIONAL EXAMINATIONS GCE Ordinary Lvl MARK SCHEME for th Octobr/Novmbr 0 sris 40 ADDITIONAL MATHEMATICS 40/ Papr, maimum raw mark 80 This mark schm is publishd as an aid to tachrs and candidats,

More information

Linear Non-Gaussian Structural Equation Models

Linear Non-Gaussian Structural Equation Models IMPS 8, Durham, NH Linar Non-Gaussian Structural Equation Modls Shohi Shimizu, Patrik Hoyr and Aapo Hyvarinn Osaka Univrsity, Japan Univrsity of Hlsinki, Finland Abstract Linar Structural Equation Modling

More information

COMP108 Algorithmic Foundations

COMP108 Algorithmic Foundations Grdy mthods Prudn Wong http://www.s.liv..uk/~pwong/thing/omp108/01617 Coin Chng Prolm Suppos w hv 3 typs of oins 10p 0p 50p Minimum numr of oins to mk 0.8, 1.0, 1.? Grdy mthod Lrning outoms Undrstnd wht

More information

Bifurcation Theory. , a stationary point, depends on the value of α. At certain values

Bifurcation Theory. , a stationary point, depends on the value of α. At certain values Dnamic Macroconomic Thor Prof. Thomas Lux Bifurcation Thor Bifurcation: qualitativ chang in th natur of th solution occurs if a paramtr passs through a critical point bifurcation or branch valu. Local

More information

Lecture 37 (Schrödinger Equation) Physics Spring 2018 Douglas Fields

Lecture 37 (Schrödinger Equation) Physics Spring 2018 Douglas Fields Lctur 37 (Schrödingr Equation) Physics 6-01 Spring 018 Douglas Filds Rducd Mass OK, so th Bohr modl of th atom givs nrgy lvls: E n 1 k m n 4 But, this has on problm it was dvlopd assuming th acclration

More information

N1.1 Homework Answers

N1.1 Homework Answers Camrig Essntials Mathmatis Cor 8 N. Homwork Answrs N. Homwork Answrs a i 6 ii i 0 ii 3 2 Any pairs of numrs whih satisfy th alulation. For xampl a 4 = 3 3 6 3 = 3 4 6 2 2 8 2 3 3 2 8 5 5 20 30 4 a 5 a

More information

What are those βs anyway? Understanding Design Matrix & Odds ratios

What are those βs anyway? Understanding Design Matrix & Odds ratios Ral paramtr stimat WILD 750 - Wildlif Population Analysis of 6 What ar thos βs anyway? Undrsting Dsign Matrix & Odds ratios Rfrncs Hosmr D.W.. Lmshow. 000. Applid logistic rgrssion. John Wily & ons Inc.

More information

COHORT MBA. Exponential function. MATH review (part2) by Lucian Mitroiu. The LOG and EXP functions. Properties: e e. lim.

COHORT MBA. Exponential function. MATH review (part2) by Lucian Mitroiu. The LOG and EXP functions. Properties: e e. lim. MTH rviw part b Lucian Mitroiu Th LOG and EXP functions Th ponntial function p : R, dfind as Proprtis: lim > lim p Eponntial function Y 8 6 - -8-6 - - X Th natural logarithm function ln in US- log: function

More information