\(\newcommand{\footnotename}{footnote}\)
\(\def \LWRfootnote {1}\)
\(\newcommand {\footnote }[2][\LWRfootnote ]{{}^{\mathrm {#1}}}\)
\(\newcommand {\footnotemark }[1][\LWRfootnote ]{{}^{\mathrm {#1}}}\)
\(\newcommand \ensuremath [1]{#1}\)
\(\newcommand {\LWRframebox }[2][]{\fbox {#2}} \newcommand {\framebox }[1][]{\LWRframebox } \)
\(\newcommand {\setlength }[2]{}\)
\(\newcommand {\addtolength }[2]{}\)
\(\newcommand {\setcounter }[2]{}\)
\(\newcommand {\addtocounter }[2]{}\)
\(\newcommand {\cline }[1]{}\)
\(\newcommand {\directlua }[1]{\text {(directlua)}}\)
\(\newcommand {\luatexdirectlua }[1]{\text {(directlua)}}\)
\(\newcommand {\protect }{}\)
\(\def \LWRabsorbnumber #1 {}\)
\(\def \LWRabsorbquotenumber "#1 {}\)
\(\def \mathchar {\ifnextchar "\LWRabsorbquotenumber \LWRabsorbnumber }\)
\(\def \mathcode #1={\mathchar }\)
\(\let \delcode \mathcode \)
\(\let \delimiter \mathchar \)
\(\let \LWRref \ref \)
\(\renewcommand {\ref }{\ifstar \LWRref \LWRref }\)
\(\newcommand {\intertext }[1]{\text {#1}\notag \\}\)
\(\newcommand {\mathllap }[2][]{{#1#2}}\)
\(\newcommand {\mathrlap }[2][]{{#1#2}}\)
\(\newcommand {\mathclap }[2][]{{#1#2}}\)
\(\newcommand {\mathmbox }[1]{#1}\)
\(\newcommand {\clap }[1]{#1}\)
\(\newcommand {\LWRmathmakebox }[2][]{#2}\)
\(\newcommand {\mathmakebox }[1][]{\LWRmathmakebox }\)
\(\newcommand {\cramped }[2][]{{#1#2}}\)
\(\newcommand {\crampedllap }[2][]{{#1#2}}\)
\(\newcommand {\crampedrlap }[2][]{{#1#2}}\)
\(\newcommand {\crampedclap }[2][]{{#1#2}}\)
\(\newenvironment {crampedsubarray}[1]{}{}\)
\(\newcommand {\crampedsubstack }{}\)
\(\newcommand {\smashoperator }[2][]{#2\limits }\)
\(\newcommand {\adjustlimits }{}\)
\(\newcommand {\SwapAboveDisplaySkip }{}\)
\(\require {extpfeil}\)
\(\Newextarrow \xleftrightarrow {10,10}{0x2194}\)
\(\Newextarrow \xLeftarrow {10,10}{0x21d0}\)
\(\Newextarrow \xhookleftarrow {10,10}{0x21a9}\)
\(\Newextarrow \xmapsto {10,10}{0x21a6}\)
\(\Newextarrow \xRightarrow {10,10}{0x21d2}\)
\(\Newextarrow \xLeftrightarrow {10,10}{0x21d4}\)
\(\Newextarrow \xhookrightarrow {10,10}{0x21aa}\)
\(\Newextarrow \xrightharpoondown {10,10}{0x21c1}\)
\(\Newextarrow \xleftharpoondown {10,10}{0x21bd}\)
\(\Newextarrow \xrightleftharpoons {10,10}{0x21cc}\)
\(\Newextarrow \xrightharpoonup {10,10}{0x21c0}\)
\(\Newextarrow \xleftharpoonup {10,10}{0x21bc}\)
\(\Newextarrow \xleftrightharpoons {10,10}{0x21cb}\)
\(\newcommand {\LWRdounderbracket }[3]{\underset {#3}{\underline {#1}}}\)
\(\newcommand {\LWRunderbracket }[2][]{\LWRdounderbracket {#2}}\)
\(\newcommand {\underbracket }[1][]{\LWRunderbracket }\)
\(\newcommand {\LWRdooverbracket }[3]{\overset {#3}{\overline {#1}}}\)
\(\newcommand {\LWRoverbracket }[2][]{\LWRdooverbracket {#2}}\)
\(\newcommand {\overbracket }[1][]{\LWRoverbracket }\)
\(\newcommand {\LATEXunderbrace }[1]{\underbrace {#1}}\)
\(\newcommand {\LATEXoverbrace }[1]{\overbrace {#1}}\)
\(\newenvironment {matrix*}[1][]{\begin {matrix}}{\end {matrix}}\)
\(\newenvironment {pmatrix*}[1][]{\begin {pmatrix}}{\end {pmatrix}}\)
\(\newenvironment {bmatrix*}[1][]{\begin {bmatrix}}{\end {bmatrix}}\)
\(\newenvironment {Bmatrix*}[1][]{\begin {Bmatrix}}{\end {Bmatrix}}\)
\(\newenvironment {vmatrix*}[1][]{\begin {vmatrix}}{\end {vmatrix}}\)
\(\newenvironment {Vmatrix*}[1][]{\begin {Vmatrix}}{\end {Vmatrix}}\)
\(\newenvironment {smallmatrix*}[1][]{\begin {matrix}}{\end {matrix}}\)
\(\newenvironment {psmallmatrix*}[1][]{\begin {pmatrix}}{\end {pmatrix}}\)
\(\newenvironment {bsmallmatrix*}[1][]{\begin {bmatrix}}{\end {bmatrix}}\)
\(\newenvironment {Bsmallmatrix*}[1][]{\begin {Bmatrix}}{\end {Bmatrix}}\)
\(\newenvironment {vsmallmatrix*}[1][]{\begin {vmatrix}}{\end {vmatrix}}\)
\(\newenvironment {Vsmallmatrix*}[1][]{\begin {Vmatrix}}{\end {Vmatrix}}\)
\(\newenvironment {psmallmatrix}[1][]{\begin {pmatrix}}{\end {pmatrix}}\)
\(\newenvironment {bsmallmatrix}[1][]{\begin {bmatrix}}{\end {bmatrix}}\)
\(\newenvironment {Bsmallmatrix}[1][]{\begin {Bmatrix}}{\end {Bmatrix}}\)
\(\newenvironment {vsmallmatrix}[1][]{\begin {vmatrix}}{\end {vmatrix}}\)
\(\newenvironment {Vsmallmatrix}[1][]{\begin {Vmatrix}}{\end {Vmatrix}}\)
\(\newcommand {\LWRmultlined }[1][]{\begin {multline*}}\)
\(\newenvironment {multlined}[1][]{\LWRmultlined }{\end {multline*}}\)
\(\let \LWRorigshoveleft \shoveleft \)
\(\renewcommand {\shoveleft }[1][]{\LWRorigshoveleft }\)
\(\let \LWRorigshoveright \shoveright \)
\(\renewcommand {\shoveright }[1][]{\LWRorigshoveright }\)
\(\newenvironment {dcases}{\begin {cases}}{\end {cases}}\)
\(\newenvironment {dcases*}{\begin {cases}}{\end {cases}}\)
\(\newenvironment {rcases}{\begin {cases}}{\end {cases}}\)
\(\newenvironment {rcases*}{\begin {cases}}{\end {cases}}\)
\(\newenvironment {drcases}{\begin {cases}}{\end {cases}}\)
\(\newenvironment {drcases*}{\begin {cases}}{\end {cases}}\)
\(\newenvironment {cases*}{\begin {cases}}{\end {cases}}\)
\(\newcommand {\MoveEqLeft }[1][]{}\)
\(\def \LWRAboxed #1!|!{\fbox {\(#1\)}&\fbox {\(#2\)}} \newcommand {\Aboxed }[1]{\LWRAboxed #1&&!|!} \)
\( \newcommand {\LWRABLines }[1][\Updownarrow ]{#1 \notag \\}\newcommand {\ArrowBetweenLines }{\ifstar \LWRABLines \LWRABLines } \)
\(\newcommand {\shortintertext }[1]{\text {#1}\notag \\}\)
\(\newcommand {\vdotswithin }[1]{\hspace {.5em}\vdots }\)
\(\newcommand {\LWRshortvdotswithinstar }[1]{\vdots \hspace {.5em} & \\}\)
\(\newcommand {\LWRshortvdotswithinnostar }[1]{& \hspace {.5em}\vdots \\}\)
\(\newcommand {\shortvdotswithin }{\ifstar \LWRshortvdotswithinstar \LWRshortvdotswithinnostar }\)
\(\newcommand {\MTFlushSpaceAbove }{}\)
\(\newcommand {\MTFlushSpaceBelow }{\\}\)
\(\newcommand \lparen {(}\)
\(\newcommand \rparen {)}\)
\(\newcommand {\ordinarycolon }{:}\)
\(\newcommand {\vcentcolon }{\mathrel {\mathop \ordinarycolon }}\)
\(\newcommand \dblcolon {\vcentcolon \vcentcolon }\)
\(\newcommand \coloneqq {\vcentcolon =}\)
\(\newcommand \Coloneqq {\dblcolon =}\)
\(\newcommand \coloneq {\vcentcolon {-}}\)
\(\newcommand \Coloneq {\dblcolon {-}}\)
\(\newcommand \eqqcolon {=\vcentcolon }\)
\(\newcommand \Eqqcolon {=\dblcolon }\)
\(\newcommand \eqcolon {\mathrel {-}\vcentcolon }\)
\(\newcommand \Eqcolon {\mathrel {-}\dblcolon }\)
\(\newcommand \colonapprox {\vcentcolon \approx }\)
\(\newcommand \Colonapprox {\dblcolon \approx }\)
\(\newcommand \colonsim {\vcentcolon \sim }\)
\(\newcommand \Colonsim {\dblcolon \sim }\)
\(\newcommand {\nuparrow }{\mathrel {\cancel {\uparrow }}}\)
\(\newcommand {\ndownarrow }{\mathrel {\cancel {\downarrow }}}\)
\(\newcommand {\bigtimes }{\mathop {\Large \times }\limits }\)
\(\newcommand {\prescript }[3]{{}^{#1}_{#2}#3}\)
\(\newenvironment {lgathered}{\begin {gathered}}{\end {gathered}}\)
\(\newenvironment {rgathered}{\begin {gathered}}{\end {gathered}}\)
\(\newcommand {\splitfrac }[2]{{}^{#1}_{#2}}\)
\(\let \splitdfrac \splitfrac \)
\(\newcommand {\LWRoverlaysymbols }[2]{\mathord {\smash {\mathop {#2\strut }\limits ^{\smash {\lower 3ex{#1}}}}\strut }}\)
\(\newcommand{\alphaup}{\unicode{x03B1}}\)
\(\newcommand{\betaup}{\unicode{x03B2}}\)
\(\newcommand{\gammaup}{\unicode{x03B3}}\)
\(\newcommand{\digammaup}{\unicode{x03DD}}\)
\(\newcommand{\deltaup}{\unicode{x03B4}}\)
\(\newcommand{\epsilonup}{\unicode{x03F5}}\)
\(\newcommand{\varepsilonup}{\unicode{x03B5}}\)
\(\newcommand{\zetaup}{\unicode{x03B6}}\)
\(\newcommand{\etaup}{\unicode{x03B7}}\)
\(\newcommand{\thetaup}{\unicode{x03B8}}\)
\(\newcommand{\varthetaup}{\unicode{x03D1}}\)
\(\newcommand{\iotaup}{\unicode{x03B9}}\)
\(\newcommand{\kappaup}{\unicode{x03BA}}\)
\(\newcommand{\varkappaup}{\unicode{x03F0}}\)
\(\newcommand{\lambdaup}{\unicode{x03BB}}\)
\(\newcommand{\muup}{\unicode{x03BC}}\)
\(\newcommand{\nuup}{\unicode{x03BD}}\)
\(\newcommand{\xiup}{\unicode{x03BE}}\)
\(\newcommand{\omicronup}{\unicode{x03BF}}\)
\(\newcommand{\piup}{\unicode{x03C0}}\)
\(\newcommand{\varpiup}{\unicode{x03D6}}\)
\(\newcommand{\rhoup}{\unicode{x03C1}}\)
\(\newcommand{\varrhoup}{\unicode{x03F1}}\)
\(\newcommand{\sigmaup}{\unicode{x03C3}}\)
\(\newcommand{\varsigmaup}{\unicode{x03C2}}\)
\(\newcommand{\tauup}{\unicode{x03C4}}\)
\(\newcommand{\upsilonup}{\unicode{x03C5}}\)
\(\newcommand{\phiup}{\unicode{x03D5}}\)
\(\newcommand{\varphiup}{\unicode{x03C6}}\)
\(\newcommand{\chiup}{\unicode{x03C7}}\)
\(\newcommand{\psiup}{\unicode{x03C8}}\)
\(\newcommand{\omegaup}{\unicode{x03C9}}\)
\(\newcommand{\Alphaup}{\unicode{x0391}}\)
\(\newcommand{\Betaup}{\unicode{x0392}}\)
\(\newcommand{\Gammaup}{\unicode{x0393}}\)
\(\newcommand{\Digammaup}{\unicode{x03DC}}\)
\(\newcommand{\Deltaup}{\unicode{x0394}}\)
\(\newcommand{\Epsilonup}{\unicode{x0395}}\)
\(\newcommand{\Zetaup}{\unicode{x0396}}\)
\(\newcommand{\Etaup}{\unicode{x0397}}\)
\(\newcommand{\Thetaup}{\unicode{x0398}}\)
\(\newcommand{\Varthetaup}{\unicode{x03F4}}\)
\(\newcommand{\Iotaup}{\unicode{x0399}}\)
\(\newcommand{\Kappaup}{\unicode{x039A}}\)
\(\newcommand{\Lambdaup}{\unicode{x039B}}\)
\(\newcommand{\Muup}{\unicode{x039C}}\)
\(\newcommand{\Nuup}{\unicode{x039D}}\)
\(\newcommand{\Xiup}{\unicode{x039E}}\)
\(\newcommand{\Omicronup}{\unicode{x039F}}\)
\(\newcommand{\Piup}{\unicode{x03A0}}\)
\(\newcommand{\Varpiup}{\unicode{x03D6}}\)
\(\newcommand{\Rhoup}{\unicode{x03A1}}\)
\(\newcommand{\Sigmaup}{\unicode{x03A3}}\)
\(\newcommand{\Tauup}{\unicode{x03A4}}\)
\(\newcommand{\Upsilonup}{\unicode{x03A5}}\)
\(\newcommand{\Phiup}{\unicode{x03A6}}\)
\(\newcommand{\Chiup}{\unicode{x03A7}}\)
\(\newcommand{\Psiup}{\unicode{x03A8}}\)
\(\newcommand{\Omegaup}{\unicode{x03A9}}\)
\(\newcommand{\alphait}{\unicode{x1D6FC}}\)
\(\newcommand{\betait}{\unicode{x1D6FD}}\)
\(\newcommand{\gammait}{\unicode{x1D6FE}}\)
\(\newcommand{\digammait}{\mathit{\unicode{x03DD}}}\)
\(\newcommand{\deltait}{\unicode{x1D6FF}}\)
\(\newcommand{\epsilonit}{\unicode{x1D716}}\)
\(\newcommand{\varepsilonit}{\unicode{x1D700}}\)
\(\newcommand{\zetait}{\unicode{x1D701}}\)
\(\newcommand{\etait}{\unicode{x1D702}}\)
\(\newcommand{\thetait}{\unicode{x1D703}}\)
\(\newcommand{\varthetait}{\unicode{x1D717}}\)
\(\newcommand{\iotait}{\unicode{x1D704}}\)
\(\newcommand{\kappait}{\unicode{x1D705}}\)
\(\newcommand{\varkappait}{\unicode{x1D718}}\)
\(\newcommand{\lambdait}{\unicode{x1D706}}\)
\(\newcommand{\muit}{\unicode{x1D707}}\)
\(\newcommand{\nuit}{\unicode{x1D708}}\)
\(\newcommand{\xiit}{\unicode{x1D709}}\)
\(\newcommand{\omicronit}{\unicode{x1D70A}}\)
\(\newcommand{\piit}{\unicode{x1D70B}}\)
\(\newcommand{\varpiit}{\unicode{x1D71B}}\)
\(\newcommand{\rhoit}{\unicode{x1D70C}}\)
\(\newcommand{\varrhoit}{\unicode{x1D71A}}\)
\(\newcommand{\sigmait}{\unicode{x1D70E}}\)
\(\newcommand{\varsigmait}{\unicode{x1D70D}}\)
\(\newcommand{\tauit}{\unicode{x1D70F}}\)
\(\newcommand{\upsilonit}{\unicode{x1D710}}\)
\(\newcommand{\phiit}{\unicode{x1D719}}\)
\(\newcommand{\varphiit}{\unicode{x1D711}}\)
\(\newcommand{\chiit}{\unicode{x1D712}}\)
\(\newcommand{\psiit}{\unicode{x1D713}}\)
\(\newcommand{\omegait}{\unicode{x1D714}}\)
\(\newcommand{\Alphait}{\unicode{x1D6E2}}\)
\(\newcommand{\Betait}{\unicode{x1D6E3}}\)
\(\newcommand{\Gammait}{\unicode{x1D6E4}}\)
\(\newcommand{\Digammait}{\mathit{\unicode{x03DC}}}\)
\(\newcommand{\Deltait}{\unicode{x1D6E5}}\)
\(\newcommand{\Epsilonit}{\unicode{x1D6E6}}\)
\(\newcommand{\Zetait}{\unicode{x1D6E7}}\)
\(\newcommand{\Etait}{\unicode{x1D6E8}}\)
\(\newcommand{\Thetait}{\unicode{x1D6E9}}\)
\(\newcommand{\Varthetait}{\unicode{x1D6F3}}\)
\(\newcommand{\Iotait}{\unicode{x1D6EA}}\)
\(\newcommand{\Kappait}{\unicode{x1D6EB}}\)
\(\newcommand{\Lambdait}{\unicode{x1D6EC}}\)
\(\newcommand{\Muit}{\unicode{x1D6ED}}\)
\(\newcommand{\Nuit}{\unicode{x1D6EE}}\)
\(\newcommand{\Xiit}{\unicode{x1D6EF}}\)
\(\newcommand{\Omicronit}{\unicode{x1D6F0}}\)
\(\newcommand{\Piit}{\unicode{x1D6F1}}\)
\(\newcommand{\Rhoit}{\unicode{x1D6F2}}\)
\(\newcommand{\Sigmait}{\unicode{x1D6F4}}\)
\(\newcommand{\Tauit}{\unicode{x1D6F5}}\)
\(\newcommand{\Upsilonit}{\unicode{x1D6F6}}\)
\(\newcommand{\Phiit}{\unicode{x1D6F7}}\)
\(\newcommand{\Chiit}{\unicode{x1D6F8}}\)
\(\newcommand{\Psiit}{\unicode{x1D6F9}}\)
\(\newcommand{\Omegait}{\unicode{x1D6FA}}\)
\(\let \digammaup \Digammaup \)
\(\renewcommand {\digammait }{\mathit {\digammaup }}\)
\(\newcommand {\smallin }{\unicode {x220A}}\)
\(\newcommand {\smallowns }{\unicode {x220D}}\)
\(\newcommand {\notsmallin }{\LWRoverlaysymbols {/}{\unicode {x220A}}}\)
\(\newcommand {\notsmallowns }{\LWRoverlaysymbols {/}{\unicode {x220D}}}\)
\(\newcommand {\rightangle }{\unicode {x221F}}\)
\(\newcommand {\intclockwise }{\unicode {x2231}}\)
\(\newcommand {\ointclockwise }{\unicode {x2232}}\)
\(\newcommand {\ointctrclockwise }{\unicode {x2233}}\)
\(\newcommand {\oiint }{\unicode {x222F}}\)
\(\newcommand {\oiiint }{\unicode {x2230}}\)
\(\newcommand {\ddag }{\unicode {x2021}}\)
\(\newcommand {\P }{\unicode {x00B6}}\)
\(\newcommand {\copyright }{\unicode {x00A9}}\)
\(\newcommand {\dag }{\unicode {x2020}}\)
\(\newcommand {\pounds }{\unicode {x00A3}}\)
\(\newcommand {\iddots }{\unicode {x22F0}}\)
\(\newcommand {\utimes }{\overline {\times }}\)
\(\newcommand {\dtimes }{\underline {\times }}\)
\(\newcommand {\udtimes }{\overline {\underline {\times }}}\)
\(\newcommand {\leftwave }{\left \{}\)
\(\newcommand {\rightwave }{\right \}}\)
\(\newcommand {\toprule }[1][]{\hline }\)
\(\let \midrule \toprule \)
\(\let \bottomrule \toprule \)
\(\newcommand {\cmidrule }[2][]{}\)
\(\newcommand {\morecmidrules }{}\)
\(\newcommand {\specialrule }[3]{\hline }\)
\(\newcommand {\addlinespace }[1][]{}\)
\(\newcommand {\LWRsubmultirow }[2][]{#2}\)
\(\newcommand {\LWRmultirow }[2][]{\LWRsubmultirow }\)
\(\newcommand {\multirow }[2][]{\LWRmultirow }\)
\(\newcommand {\mrowcell }{}\)
\(\newcommand {\mcolrowcell }{}\)
\(\newcommand {\STneed }[1]{}\)
\( \newcommand {\multicolumn }[3]{#3}\)
\(\newcommand {\tothe }[1]{^{#1}}\)
\(\newcommand {\raiseto }[2]{{#2}^{#1}}\)
\(\newcommand {\ang }[2][]{(\mathrm {#2})\degree }\)
\(\newcommand {\num }[2][]{\mathrm {#2}}\)
\(\newcommand {\si }[2][]{\mathrm {#2}}\)
\(\newcommand {\LWRSI }[2][]{\mathrm {#1\LWRSInumber \,#2}}\)
\(\newcommand {\SI }[2][]{\def \LWRSInumber {#2}\LWRSI }\)
\(\newcommand {\numlist }[2][]{\mathrm {#2}}\)
\(\newcommand {\numrange }[3][]{\mathrm {#2\,\unicode {x2013}\,#3}}\)
\(\newcommand {\SIlist }[3][]{\mathrm {#2\,#3}}\)
\(\newcommand {\SIrange }[4][]{\mathrm {#2\,#4\,\unicode {x2013}\,#3\,#4}}\)
\(\newcommand {\tablenum }[2][]{\mathrm {#2}}\)
\(\newcommand {\ampere }{\mathrm {A}}\)
\(\newcommand {\candela }{\mathrm {cd}}\)
\(\newcommand {\kelvin }{\mathrm {K}}\)
\(\newcommand {\kilogram }{\mathrm {kg}}\)
\(\newcommand {\metre }{\mathrm {m}}\)
\(\newcommand {\mole }{\mathrm {mol}}\)
\(\newcommand {\second }{\mathrm {s}}\)
\(\newcommand {\becquerel }{\mathrm {Bq}}\)
\(\newcommand {\degreeCelsius }{\unicode {x2103}}\)
\(\newcommand {\coulomb }{\mathrm {C}}\)
\(\newcommand {\farad }{\mathrm {F}}\)
\(\newcommand {\gray }{\mathrm {Gy}}\)
\(\newcommand {\hertz }{\mathrm {Hz}}\)
\(\newcommand {\henry }{\mathrm {H}}\)
\(\newcommand {\joule }{\mathrm {J}}\)
\(\newcommand {\katal }{\mathrm {kat}}\)
\(\newcommand {\lumen }{\mathrm {lm}}\)
\(\newcommand {\lux }{\mathrm {lx}}\)
\(\newcommand {\newton }{\mathrm {N}}\)
\(\newcommand {\ohm }{\mathrm {\Omega }}\)
\(\newcommand {\pascal }{\mathrm {Pa}}\)
\(\newcommand {\radian }{\mathrm {rad}}\)
\(\newcommand {\siemens }{\mathrm {S}}\)
\(\newcommand {\sievert }{\mathrm {Sv}}\)
\(\newcommand {\steradian }{\mathrm {sr}}\)
\(\newcommand {\tesla }{\mathrm {T}}\)
\(\newcommand {\volt }{\mathrm {V}}\)
\(\newcommand {\watt }{\mathrm {W}}\)
\(\newcommand {\weber }{\mathrm {Wb}}\)
\(\newcommand {\day }{\mathrm {d}}\)
\(\newcommand {\degree }{\mathrm {^\circ }}\)
\(\newcommand {\hectare }{\mathrm {ha}}\)
\(\newcommand {\hour }{\mathrm {h}}\)
\(\newcommand {\litre }{\mathrm {l}}\)
\(\newcommand {\liter }{\mathrm {L}}\)
\(\newcommand {\arcminute }{^\prime }\)
\(\newcommand {\minute }{\mathrm {min}}\)
\(\newcommand {\arcsecond }{^{\prime \prime }}\)
\(\newcommand {\tonne }{\mathrm {t}}\)
\(\newcommand {\astronomicalunit }{au}\)
\(\newcommand {\atomicmassunit }{u}\)
\(\newcommand {\bohr }{\mathit {a}_0}\)
\(\newcommand {\clight }{\mathit {c}_0}\)
\(\newcommand {\dalton }{\mathrm {D}_\mathrm {a}}\)
\(\newcommand {\electronmass }{\mathit {m}_{\mathrm {e}}}\)
\(\newcommand {\electronvolt }{\mathrm {eV}}\)
\(\newcommand {\elementarycharge }{\mathit {e}}\)
\(\newcommand {\hartree }{\mathit {E}_{\mathrm {h}}}\)
\(\newcommand {\planckbar }{\mathit {\unicode {x210F}}}\)
\(\newcommand {\angstrom }{\mathrm {\unicode {x212B}}}\)
\(\let \LWRorigbar \bar \)
\(\newcommand {\bar }{\mathrm {bar}}\)
\(\newcommand {\barn }{\mathrm {b}}\)
\(\newcommand {\bel }{\mathrm {B}}\)
\(\newcommand {\decibel }{\mathrm {dB}}\)
\(\newcommand {\knot }{\mathrm {kn}}\)
\(\newcommand {\mmHg }{\mathrm {mmHg}}\)
\(\newcommand {\nauticalmile }{\mathrm {M}}\)
\(\newcommand {\neper }{\mathrm {Np}}\)
\(\newcommand {\yocto }{\mathrm {y}}\)
\(\newcommand {\zepto }{\mathrm {z}}\)
\(\newcommand {\atto }{\mathrm {a}}\)
\(\newcommand {\femto }{\mathrm {f}}\)
\(\newcommand {\pico }{\mathrm {p}}\)
\(\newcommand {\nano }{\mathrm {n}}\)
\(\newcommand {\micro }{\mathrm {\unicode {x00B5}}}\)
\(\newcommand {\milli }{\mathrm {m}}\)
\(\newcommand {\centi }{\mathrm {c}}\)
\(\newcommand {\deci }{\mathrm {d}}\)
\(\newcommand {\deca }{\mathrm {da}}\)
\(\newcommand {\hecto }{\mathrm {h}}\)
\(\newcommand {\kilo }{\mathrm {k}}\)
\(\newcommand {\mega }{\mathrm {M}}\)
\(\newcommand {\giga }{\mathrm {G}}\)
\(\newcommand {\tera }{\mathrm {T}}\)
\(\newcommand {\peta }{\mathrm {P}}\)
\(\newcommand {\exa }{\mathrm {E}}\)
\(\newcommand {\zetta }{\mathrm {Z}}\)
\(\newcommand {\yotta }{\mathrm {Y}}\)
\(\newcommand {\percent }{\mathrm {\%}}\)
\(\newcommand {\meter }{\mathrm {m}}\)
\(\newcommand {\metre }{\mathrm {m}}\)
\(\newcommand {\gram }{\mathrm {g}}\)
\(\newcommand {\kg }{\kilo \gram }\)
\(\newcommand {\of }[1]{_{\mathrm {#1}}}\)
\(\newcommand {\squared }{^2}\)
\(\newcommand {\square }[1]{\mathrm {#1}^2}\)
\(\newcommand {\cubed }{^3}\)
\(\newcommand {\cubic }[1]{\mathrm {#1}^3}\)
\(\newcommand {\per }{/}\)
\(\newcommand {\celsius }{\unicode {x2103}}\)
\(\newcommand {\fg }{\femto \gram }\)
\(\newcommand {\pg }{\pico \gram }\)
\(\newcommand {\ng }{\nano \gram }\)
\(\newcommand {\ug }{\micro \gram }\)
\(\newcommand {\mg }{\milli \gram }\)
\(\newcommand {\g }{\gram }\)
\(\newcommand {\kg }{\kilo \gram }\)
\(\newcommand {\amu }{\mathrm {u}}\)
\(\newcommand {\nm }{\nano \metre }\)
\(\newcommand {\um }{\micro \metre }\)
\(\newcommand {\mm }{\milli \metre }\)
\(\newcommand {\cm }{\centi \metre }\)
\(\newcommand {\dm }{\deci \metre }\)
\(\newcommand {\m }{\metre }\)
\(\newcommand {\km }{\kilo \metre }\)
\(\newcommand {\as }{\atto \second }\)
\(\newcommand {\fs }{\femto \second }\)
\(\newcommand {\ps }{\pico \second }\)
\(\newcommand {\ns }{\nano \second }\)
\(\newcommand {\us }{\micro \second }\)
\(\newcommand {\ms }{\milli \second }\)
\(\newcommand {\s }{\second }\)
\(\newcommand {\fmol }{\femto \mol }\)
\(\newcommand {\pmol }{\pico \mol }\)
\(\newcommand {\nmol }{\nano \mol }\)
\(\newcommand {\umol }{\micro \mol }\)
\(\newcommand {\mmol }{\milli \mol }\)
\(\newcommand {\mol }{\mol }\)
\(\newcommand {\kmol }{\kilo \mol }\)
\(\newcommand {\pA }{\pico \ampere }\)
\(\newcommand {\nA }{\nano \ampere }\)
\(\newcommand {\uA }{\micro \ampere }\)
\(\newcommand {\mA }{\milli \ampere }\)
\(\newcommand {\A }{\ampere }\)
\(\newcommand {\kA }{\kilo \ampere }\)
\(\newcommand {\ul }{\micro \litre }\)
\(\newcommand {\ml }{\milli \litre }\)
\(\newcommand {\l }{\litre }\)
\(\newcommand {\hl }{\hecto \litre }\)
\(\newcommand {\uL }{\micro \liter }\)
\(\newcommand {\mL }{\milli \liter }\)
\(\newcommand {\L }{\liter }\)
\(\newcommand {\hL }{\hecto \liter }\)
\(\newcommand {\mHz }{\milli \hertz }\)
\(\newcommand {\Hz }{\hertz }\)
\(\newcommand {\kHz }{\kilo \hertz }\)
\(\newcommand {\MHz }{\mega \hertz }\)
\(\newcommand {\GHz }{\giga \hertz }\)
\(\newcommand {\THz }{\tera \hertz }\)
\(\newcommand {\mN }{\milli \newton }\)
\(\newcommand {\N }{\newton }\)
\(\newcommand {\kN }{\kilo \newton }\)
\(\newcommand {\MN }{\mega \newton }\)
\(\newcommand {\Pa }{\pascal }\)
\(\newcommand {\kPa }{\kilo \pascal }\)
\(\newcommand {\MPa }{\mega \pascal }\)
\(\newcommand {\GPa }{\giga \pascal }\)
\(\newcommand {\mohm }{\milli \ohm }\)
\(\newcommand {\kohm }{\kilo \ohm }\)
\(\newcommand {\Mohm }{\mega \ohm }\)
\(\newcommand {\pV }{\pico \volt }\)
\(\newcommand {\nV }{\nano \volt }\)
\(\newcommand {\uV }{\micro \volt }\)
\(\newcommand {\mV }{\milli \volt }\)
\(\newcommand {\V }{\volt }\)
\(\newcommand {\kV }{\kilo \volt }\)
\(\newcommand {\W }{\watt }\)
\(\newcommand {\uW }{\micro \watt }\)
\(\newcommand {\mW }{\milli \watt }\)
\(\newcommand {\kW }{\kilo \watt }\)
\(\newcommand {\MW }{\mega \watt }\)
\(\newcommand {\GW }{\giga \watt }\)
\(\newcommand {\J }{\joule }\)
\(\newcommand {\uJ }{\micro \joule }\)
\(\newcommand {\mJ }{\milli \joule }\)
\(\newcommand {\kJ }{\kilo \joule }\)
\(\newcommand {\eV }{\electronvolt }\)
\(\newcommand {\meV }{\milli \electronvolt }\)
\(\newcommand {\keV }{\kilo \electronvolt }\)
\(\newcommand {\MeV }{\mega \electronvolt }\)
\(\newcommand {\GeV }{\giga \electronvolt }\)
\(\newcommand {\TeV }{\tera \electronvolt }\)
\(\newcommand {\kWh }{\kilo \watt \hour }\)
\(\newcommand {\F }{\farad }\)
\(\newcommand {\fF }{\femto \farad }\)
\(\newcommand {\pF }{\pico \farad }\)
\(\newcommand {\K }{\mathrm {K}}\)
\(\newcommand {\dB }{\mathrm {dB}}\)
\(\newcommand {\kibi }{\mathrm {Ki}}\)
\(\newcommand {\mebi }{\mathrm {Mi}}\)
\(\newcommand {\gibi }{\mathrm {Gi}}\)
\(\newcommand {\tebi }{\mathrm {Ti}}\)
\(\newcommand {\pebi }{\mathrm {Pi}}\)
\(\newcommand {\exbi }{\mathrm {Ei}}\)
\(\newcommand {\zebi }{\mathrm {Zi}}\)
\(\newcommand {\yobi }{\mathrm {Yi}}\)
\(\require {mhchem}\)
\(\require {cancel}\)
\(\newcommand {\fint }{âĺŊ}\)
\(\newcommand {\hdots }{\cdots }\)
\(\newcommand {\mathnormal }[1]{#1}\)
\(\newcommand {\vecs }[2]{\vec {#1}_{#2}}\)
\(\newcommand {\code }[1]{\lstinline [basicstyle=\ttfamily ]{#1}}\)
\(\renewcommand {\O }{\ensuremath {\mathcal {O}}}\)
\(\newcommand {\source }{\ensuremath {\operatorname {source}}}\)
\(\newcommand {\target }{\ensuremath {\operatorname {target}}}\)
\(\renewcommand {\num }{\ensuremath {\operatorname {num}}}\)
\(\newcommand {\indeg }{\ensuremath {\operatorname {indeg}}}\)
\(\newcommand {\mycap }{\ensuremath {\operatorname {cap}}}\)
\(\newcommand {\cost }{\ensuremath {\operatorname {cost}}}\)
\(\newcommand {\pot }{\ensuremath {\operatorname {pot}}}\)
\(\newcommand {\ccost }{\ensuremath {\operatorname {ccost}}}\)
\(\newcommand {\dcost }{\ensuremath {\operatorname {dcost}}}\)
\(\newcommand {\icost }{\ensuremath {\operatorname {icost}}}\)
\(\newcommand {\qed }{\hfill $\Box $}\)
\(\newcommand {\fracsize }[1]{{\large #1}}\)
\(\newcommand {\matrixsize }[1]{{\scriptsize #1}}\)
\(\newcommand {\name }[1]{\textsc {#1}}\)
\(\newcommand {\smallpmatrix }[1]{\left (\begin {smallmatrix}#1\end {smallmatrix}\right )}\)
\(\newcommand {\matlab }{{\fontfamily {bch}\scshape \selectfont {}Matlab}}\)
\(\newcommand {\innerproduct }[1]{\left \langle {#1}\right \rangle }\)
\(\newcommand {\norm }[1]{\left \Vert {#1}\right \Vert }\)
\(\renewcommand {\natural }{\mathbb {N}}\)
\(\newcommand {\integer }{\mathbb {Z}}\)
\(\newcommand {\rational }{\mathbb {Q}}\)
\(\newcommand {\real }{\mathbb {R}}\)
\(\newcommand {\complex }{\mathbb {C}}\)
\(\renewcommand {\d }{\mathop {}\!\mathrm {d}}\)
\(\newcommand {\dr }{\d {}r}\)
\(\newcommand {\ds }{\d {}s}\)
\(\newcommand {\dt }{\d {}t}\)
\(\newcommand {\du }{\d {}u}\)
\(\newcommand {\dv }{\d {}v}\)
\(\newcommand {\dw }{\d {}w}\)
\(\newcommand {\dx }{\d {}x}\)
\(\newcommand {\dy }{\d {}y}\)
\(\newcommand {\dz }{\d {}z}\)
\(\newcommand {\dsigma }{\d {}\sigma }\)
\(\newcommand {\dphi }{\d {}\phi }\)
\(\newcommand {\dvarphi }{\d {}\varphi }\)
\(\newcommand {\dtau }{\d {}\tau }\)
\(\newcommand {\dxi }{\d {}\xi }\)
\(\newcommand {\dtheta }{\d {}\theta }\)
\(\newcommand {\tp }{\mathrm {T}}\)
Allgemeines zu Suchbäumen
Gegeben sei eine Teilmenge \(S = \{a_1, \dotsc , a_n\}\) eines geordneten Universums \(U\).
Gesucht ist eine Organisation von \(S\) in einem binären Baum, wobei
der Baum \(n\) Knoten und eine ausgezeichnete Wurzel besitzt,
jedes \(a_i\) mit einem Knoten asoziiert ist,
jeder innere Knoten maximal einen linken und maximal einen rechten Kindknoten hat sowie
für jeden Knoten \(a_v\) des Baums sind alle Knoten im Teilbaum eines linken/rechten Kindes von \(a_v\) kleiner/größer als \(a_v\).
Die Suche nach einem Schlüssel \(x\) im Suchbaum verläuft wie folgt:
search(\(x, a_{root}\)) gibt das größte \(a_\ell \) mit \(a_\ell \le x\) im
Baum zurück, sonst \(-\infty \).
search(x, {}) = -unendlich
search(x, a)
if x <a
return search(x, a_L)
else if x = a
return a
else
return max(a, search(a, a_R))
Suchbäume können, falls der Baum ausbalanciert ist, dabei helfen, sehr „billig“ nach einem Knoten zu suchen. In diesem Fall ist die Suchzeit für ein Knoten \(\O (\log n)\). Das Problem ist, dass
das Löschen und Einfügen von Elementen den Baum wieder unbalanciert machen kann – der Extremfall ist ein „ausgearteter Baum“, bei dem in jeder Ebene nur ein Knoten ist. Daher nutzt man spezielle
Datenstrukturen, die den Baum automatisch ausbalancieren.
(2, 4)-Bäume
Gegeben sei wieder eine linear geordnete Menge \(S = \{a_1, \dotsc , a_n\}\).
\(S\) soll nur in den Blättern des Baums gespeichert werden.
Die Blätter müssen sortiert sein.
Die Blätter des Baums sollen alle die gleiche Tiefe haben.
Jeder innere Knoten hat zwischen zwei und vier Kinder.
Jede innere Knoten mit \(i\) Kindern enthält selbst \(i - 1\) Schlüssel, dabei ist der \(j\)-te Schlüssel das größte Element des \(j\)-ten Teilbaums des Knotens (\(j = 1, \dotsc
, i - 1\)).
Suche nach einem Schlüssel \(k\):
v := wurzel
while v kein Blatt do
bestimme l mit k_{l-1}(v) <k <= k_l(v) (wobei k_0(v) = -unendlich und k_grad(v) = +unendlich gesetzt wird)
v := l-tes Kind von v
od
Die Suche in einem Blatt \(v\), wobei \(Schl"ussel(linkerNachbar) < k \le Schl"ussel(v)\) ist. Ist \(k > \max S\), dann endet die Suche in dem Blatt, das am weitesten rechts liegt.
Die Laufzeit ist \(\O (h)\), wenn \(h\) die Tiefe des Baums ist.
Lemma: Sei \(T\) ein (2, 4)-Baum der Höhe \(h\) mit \(n\) Blättern.
Dann gilt \(2^h \le n \le 4^h\) und daher \(\frac {1}{2} \log n \le h \le \log n\).
Einfügen eines Elements \(k\):
Angenommen wird, dass ein Verweis auf das Blatt \(v\) mit \(Schl"ussel(linkerNachbar) < k \le Schl"ussel(v)\) vorliegt (kann durch Suchen nach \(k\) erreicht werden).
Einfügen: Füge \(k\) links von \(v\) als neues Blatt hinzu und füge \(k\) als Schlüssel in den Vaterknoten vor den Schlüssel von \(v\) ein.
Spalten: Wenn der Vaterknoten nun fünf Knoten hat, muss er aufgespaltet werden, d. h. er wird in zwei Knoten aufgeteilt, wobei der linke Knoten die ersten zwei und der rechte
Knoten die letzten drei Kinder enthält. Wenn die Wurzel gespalten werden muss, muss eine neue Wurzel erzeugt werden, sodass die Tiefe um \(1\) steigt.
v := vater(v)
while v hat fuenf Kinder do
spalte(v)
v := vater(v)
od
Die Laufzeit ist \(\O (1 + \text {Anzahl Spaltungen}) = \O (\log n)\).
Löschen eines Elements \(k\):
Löschen: Ist \(k\) kein am weitesten rechts liegendes Kind, so kann es einfach mit dem zugehörigen Schlüssel in seinem Vaterknoten gelöscht werden. Liegt
\(k\) am weitesten rechts, so muss der zugehörige Schlüssel verändert werden, der allerdings nicht im Vaterknoten liegt, sondern in einem darüber liegenden indirekten Vaterknoten.
Verschmelzen/Stehlen: Hat der Vaterknoten von \(k\) nach dem Löschen nur noch ein Kind, so muss entweder der Knoten mit einem Nachbarn verschmolzen werden oder er muss
einen Knoten von einem Nachbarknoten stehlen. Verschmolzen wird, falls der Nachbarknoten \(2\) Knoten hat, hat er \(3\) oder \(4\) Knoten, so wird gestohlen.
Dies kann bis zur Wurzel fortgesetzt werden, daher beträgt die Laufzeit \(\O (\log n)\).
Warum benutzt man nicht (2, 3)-Bäume?
Hat ein Knoten drei Kinder, so würde dieser nach einer Einfügeoperation gespalten werden. Ist jedoch der Baum „voll“, d. h. jeder Knoten hat drei Kinder, so müsste jeder gespalten werden,
sodass nun alle Knoten zwei Kinder haben. Wird nun wieder ein Knoten gelöscht, so müssten wieder alle Kinder verschmolzen werden und die Ausgangssituation wäre wiederhergestellt. Der
Zeitaufwand von (2, 3)-Bäumen ist also größer (analog zum Binärzähler).
Laufzeit: Jede beliebige Sequenz aus Einfügen und Löschen benötigt in einem (2, 4)-Baum amortisiert \(\O (1)\) Operationen.
Beweis: ?
Potential: Das Potential eines (2, 4)-Baums \(T\) mit maximal einem Knoten vom Grad 1 oder 5 ist \(\pot (T) = (2 \cdot \#1) + (1 \cdot \#2) + (0 \cdot \#3) + (2 \cdot \#4) + (4
\cdot \#5)\), wobei \(\#i\) die Anzahl der Knoten vom Grad \(i\) bedeutet.
Anwendungen von (2, 4)-Bäumen
Sortieren „leicht“ vorsortierter Folgen:
Gegeben sei eine Folge \(x_1, \dotsc , x_n\). Die Anzahl der Inversionen (Fehlstände) ist dann
\(F = \Big |\big \{(i, j) \;|\; i < j,\; x_i > x_j\big \}\Big |\), es gilt \(0 \le F \le \binom {n}{2}\).
Behauptung: Man kann mit (2, 4)-Bäumen in Zeit \(\O (n \max \{1, \log \frac {F}{n}\})\) sortieren.
Ist \(F = n\) bzw. \(F = n^2\), so kann man in \(\O (n)\) bzw. \(\O (n \log n)\) sortieren.
Beweis: Man sortiert durch Einfügen in einen (2, 4)-Baum. Angenommen, \(x_1, \dotsc , x_i\) sind schon sortiert eingefügt und \(x_{i+1}\) soll nun eingefügt werden. Der
Abstand von rechts zur eigentlich richtigen Position von \(x_{i+1}\) ist \(f_{i+1} = |\{j \;|\; j < i + 1,\; x_j > x_{i+1}\}|\), wobei \(\sum _i f_i = F\) ist.
\(x_{i+1}\) kann nun in amortisierter Zeit \(\O (1 + \max \{1, \log f_{i+1}\})\) eingefügt werden:
Man läuft vom rechten Blatt (am weitesten rechts) bis ein Knoten \(v\) erreicht wird mit \(x_{i+1} >\) größter Schlüssel \(k\) in \(v\). Dafür wird \(\O (1 + h)\) Zeit
benötigt, wenn \(h\) die Höhe von \(v\) von unten gesehen ist. Ist \(v’\) das rechte Kind von \(v\) und \(k’\) der größte Schlüssel von \(v’\), so ist \(k < x_{i+1} < k’\)
und \(x_{i+1}\) wird in einem Kind von \(v’\) eingefügt, aber nicht im rechten Knoten \(v”\) von \(v’\).
Daher gilt \(f_{i+1} \ge \text {Anzahl Bl"atter unterhalb von } v” \ge 2^{h-2}\), weil \(v”\) Höhe \(h - 2\) hat (Blätter Höhe \(0\)). Es folgt \(h \le 2 +
\log f_{i+1}\). Die Einfügung selbst (ohne Lokalisation) hat amortisierte Kosten \(\O (1)\). Daher kann \(x_{i+1}\) in \(\O (1 + \max \{1, \log f_{i+1}\})\) eingefügt werden.
Damit ist die Gesamtlaufzeit \(\O (\sum _i (1 + \max \{1, \log f_i\})) = \O (n + \sum _i \max \{1, \log f_i\}) \\ = \O (n + \sum _i (1 + \log f_i)) = \O (n + \sum _i
\log f_i) = \O (n + n \log \frac {F}{n}) = \O (n \max \{1, \log \frac {F}{n}\})\).
Die vorletzte Gleichheit erhält man mit \((\prod _i f_i)^{1/n} \le \frac {\sum _i f_i}{n} \quad |\log \)
(geometrisches Mittel ist kleiner/gleich arithmetisches Mittel).
Fingersuche (ein Finger ist ein Zeiger auf ein Blatt):
Lemma: In niveau-verbundenen (2, 4)-Bäumen kann man Fingersuche in \(\O (\log \min \{d, n - d\})\) durchführen, wobei \(d\) der Abstand des Fingers zum Ziel der Suche ist.
Niveau-verbunden heißt, dass die Kanten jeder Ebene in einer zirkulären Liste stehen, d. h. zu jedem Knoten ist der linke und rechte Nachbar bekannt (auch wenn Vaterknoten anders ist) und zu einem Knoten ganz
rechts ist der rechte Nachbar der Knoten der Ebene ganz links.
Suche von \(x\) von einem Finger aus: Laufe von dem Finger in Richtung Wurzel, bis ein Knoten \(v\) erreicht wird, sodass \(x\) unterhalb dem \(v\), dem linken Nachbar oder dem rechten Nachbar von \(v\) liegt. Dann dreht
man um und sucht ganz normal.
Die Laufzeit ist \(\O (\text {H"ohe des erreichten Knotens})\), diese ist \(\log (\min \{1 + d, n - d + 1\})\).
schnelles Mischen und Sortieren durch Mischen:
Gegeben seien sortierte Folgen \(S_1, S_2, \dotsc \) als (2, 4)-Bäume.
Ziel: Mische \(S_1\) und \(S_2\) zu \(S\) in einen (2, 4)-Baum, wobei \(|S_1| = n\) und \(|S_2| = m\) mit \(m \le n\) ist.
naiv: Füge \(S_2\) nacheinander in \(S_1\) ein. Die Laufzeit dafür ist \(\O (m \cdot \log (m + n))\), dies ist schlecht für \(m \approx n\) (gut für \(m \ll n\)).
Satz: Man kann \(S_1\) und \(S_2\) in Zeit \(\O (m \cdot \log \frac {m + n}{m}) = \O (\log \binom {m + n}{m})\) zu einem (2, 4)-Baum mischen, der \(S_1\) und \(S_2\)
enthält.
finger := "linkestes" Blatt in S_1
i := 1
while i <= m do
suche nach x_i von finger aus
fuege x_i ein
finger := Zeiger auf neues Blatt
i++
od
Alternativen für (2, 4)-Bäume
Rot-Schwarz-Bäume,
AVL-Bäume und
Skip Lists können immer dann benutzt werden, wenn Elemente mit einer Ordnung verwaltet werden sollen. Sie
ermöglichen das Suchen, Einfügen und Löschen in logarithmischer Zeit.
Manchmal kann Hashing jedoch effizienter sein, denn so ist Suchen, Einfügen und Löschen in \(\O (1)\) möglich. Dies geht aber nur, wenn die zuverwaltende Menge aus ganzen Zahlen
besteht. Außerdem können keine Anfragen der Art „größtes Element kleiner \(10\)“ beantwortet werden.
Einschub: Amortisierte Analyse
Sinn und Zweck: Man möchte zeigen, dass nicht alle Operationen auf einer bestimmten Datenstruktur teuer sind, d. h. im Durchschnitt sind die Operationen billig, auch wenn eine einzelne Operation
teuer sein kann.
Intuition: Mit jeder Operation auf der Datenstruktur wird eine konstante Zahl von Euros einbezahlt, die für den tatsächlichen Aufwand einer Operation bezahlt werden müssen, dessen Rest
aber bei billigen Operationen angespart werden kann, um teurere Operationen zu bezahlen.
Beispiel 1: Inkrementierung im Binärregister
Hier entspricht der Aufwand der Anzahl der Überträge. Im schlimmsten Fall müssen \(\log _2 n\) Überträge gemacht werden (\(n\) größte speicherbare Zahl).
Man kann zeigen: Wenn man bei Null anfängt und bei jeder Inkrementierung immer \(1\) Euro einbezahlen, so hat die Datenstruktur immer genügend Geld, um die Überträge zu bezahlen
(ein Übertrag kostet \(1\) Euro).
Hier kommt die Potentialfunktion ins Spiel: Sie ist eine untere Schranke für den Kontostand und entspricht hier der Anzahl Einsen in der aktuellen Zahl (dies müsste man zeigen).
Nach \(i\) Inkrementierungen hat man \(i\) Euro eingezahlt, der Kontostand ist nicht-negativ, d. h. man hat nicht mehr als \(i\) Euro für Überträge ausgegeben. Im Durchschnitt/amortisiert
wurden also \(\le 1\) Überträge gemacht.
Beispiel 2: Konstruktion eines Heaps in \(\O (n)\) Zeit
Beispiel 3: Spalten und Vertauschen in (2, 4)-Bäumen
(Stehlen ist uninteressant, da nicht propagierend.)
Die Potentialfunktion (untere Schranke für Kontostand) ist hier
\(\phi = 2 \cdot \#1 + 1 \cdot \#2 + 0 \cdot \#3 + 2 \cdot \#4 + 4 \cdot \#5\).
Pro Einfügen und Löschen werden \(5\) Euro auf das Konto des (2, 4)-Baums bezahlt, die gespart werden können, aber auch für Spalt-/Verschmelzoperationen ausgegeben werden
müssen. Behauptung: Der (2, 4)-Baum hat immer genügend Geld, um Spalten/Verschmelzen zu bezahlen.
Hat man dies gezeigt, so hat man, wenn man mit einem leeren Baum anfängt, nach \(i\) Operationen einen Kontostand von höchstens \(5i\) Euro. Weil der Kontostand nicht-negativ ist, ist der Aufwand
für Spalten und Verschmelzen \(\le 5i\).
Beweis:
Durch das bloße Einfügen bzw. Löschen eines Blattes erhöht sich das Potential um maximal \(2\) bzw. \(1\). Also ist \(\phi \) weiterhin eine gültige untere Schranke (\(5\) Euro wurden
eingezahlt).
Beim Spalten eines Knotens mit \(5\) Kindern erhöht sich das Potential des Vaterknotens um max. \(2\), aus dem Kind mit Potential \(4\) entstehen zwei Kinder mit Potential \(1\) und \(0\). Also nimmt das Potential
des Baums um mindestens \(1\) ab. Mit diesem Euro kann die Operation bezahlt werden und \(\phi \) bleibt untere Schranke für den Kontostand.
Beim Verschmelzen zweier Knoten mit \(1\) Kind und \(2\) Kindern erhöht sich das Potential des Vaterknotens um max. \(1\), das Potential der Kinder \(2\) und \(1\) ändert sich zu \(0\), da ein Knoten mit
drei Kindern entsteht. Also nimmt das Potential um mindestens \(2\) ab. Die Operation kann bezahlt werden und \(\phi \) bleibt untere Schranke für den Kontostand.
Beim Stehlen sind zwei Knoten mit \(1\) bzw. \(3\) oder \(4\) Kinder vorhanden. Deren Potential ändert sich von \(2\) bzw. \(0\) oder \(2\) zu \(1\) bzw. \(1\) oder \(0\). Sonst verändert sich kein
Knoten, daher auch nicht das Potential. Also nimmt bleibt das Potential gleich oder sinkt um \(3\). \(\phi \) ist weiterhin untere Schranke für Kontostand, denn es muss beim Stehlen nichts bezahlt werden.
Also bleibt \(\phi \) durchgängig untere Schranke für den Kontostand.
Beim Verschmelzen/Spalten ist immer Geld vorhanden, um die Operation zu bezahlen.