Code
3.1 Basic Distributional Quantities
Teori Resiko
**
Dalam Materi ini akan mempelajari: - moments, - percentiles, and - generating functions.
3.1.1 Moments
Dengan memisalkan X merupakan suatu variabel acak kontinu dengan probability density function (pdf) \(f_X(x)\) dan fungsi ditribusi \(F_X(x)\) .
Momen baku ke-k dari X yang dinotasikan dengan \(μ'_k\) adalah nilai ekspektasi dari pangkat ke-k dari X yang diharapkan, asalkan mempunyai nilai. Momen pertama \(μ'_1\) adalah nilai tengah (mean) dari X yang biasanya dilambangkan dengan μ . Rumus untuk \(μ'_k\) adalah sebagai berikut :
\[
\begin{align}
μ'_k = E(X^k) = \int_{0}^{∞}x^kf_X(x)dx
\end{align}
\]
Untuk mendukung dari variabel acak X dapat diasumsikan sebagai nonnegatif dikarenakan pada kenyataannya jarang sekali bernilai negatif. Salah satu contohnya yang menunjukkan bahwa raw moments untuk variabel nonnegatif dapat dihitung menggunakan rumus :
\[
\begin{align}
μ'_k = \int_{0}^{∞}kx^{k-1}[1-F_X(x)]dx,
\end{align}
\]
yang didasarkan pada fungsi survival \(S_X(x)=1−F_X(x)\) . Pada rumus ini berguna disaat k=1.
Central moment ke -k dari X yang dinotasikan dengan \(μ_k\) yang merupakan nilai yang diharapkan dari pangkat ke-k dari deviasi x dan dari mean μ. Maka untuk rumus \(μ_k\) adalah
\[
\begin{align}
μ_k = E[(X-μ)^k] = \int_{0}^{∞}(x-μ)^kf_X(x)dx
\end{align}
\]
Momen pusat kedua \(μ_2\) mendefinisikan varians dari X yang dinotasikan dengan \(σ^2\) . Akar kuadrat dari varians adalah simpangan baku \(σ\) .
Rasio momen sentral ketiga terhadap pangkat tiga dari deviasi standar \((μ_3/σ^3)\) yang mendefinisikan koefisien kemiringan yang merupakan ukuran simetri. Koefisien kemencengan yang positif menunjukkan bahwa distribusi condong ke kanan (condong ke kanan).
Rasio momen sentral keempat dengan pangkat empat dari deviasi standar \((μ_4/σ^4)\) mendefinisikan koefisien kurtosis. Distribusi normal memiliki koefisien kurtosis 3. Distribusi dengan koefisien kurtosis lebih besar dari 3 memiliki ekor yang lebih berat daripada distribusi normal, sedangkan distribusi dengan koefisien kurtosis kurang dari 3 memiliki ekor yang lebih ringan dan lebih datar.
Example 3.1.1
Mengasumsikan bahwa variabel acak X memiliki distribusi gamma dengan rata-rata 9 dan skewness 1. Maka dicari perhitungan varians dari X.
pdf dari X adalah :
\[
\begin{align}
f_{X}\left( x \right) = \frac{\left( x / \theta \right)^{\alpha}}{x ~\Gamma\left( \alpha \right)} e^{- x / \theta}
\end{align}
\] Untuk \(x>0\) , dan \(α>0\) , maka momen baku ke-k adalah
\[
\begin{align}
\mu_{k}^{\prime} = \mathrm{E}\left( X^{k} \right) = \int_{0}^{\infty}{\frac{1}{\Gamma\left( \alpha \right)\theta^{\alpha}}x^{k + \alpha - 1}e^{- x / \theta} dx} = \frac{\Gamma\left( k + \alpha \right)}{\Gamma\left( \alpha \right)}\theta^{k}
\end{align}
\]
Dengan memberikan \(\Gamma\left( r + 1 \right) = r\Gamma\left( r \right)\) dan \(\Gamma\left( 1 \right) = 1\) . Maka \(\mu_{1}^{\prime} = \mathrm{E}\left( X \right) = \alpha\theta\)
$$
\[\begin{align}
\mu_{2}^{\prime} &= \mathrm{E}\left( X^{2} \right) = \left( \alpha + 1 \right)\alpha\theta^{2} =
\mu_{3}^{\prime} = \mathrm{E}\left( X^{3} \right) = \left( \alpha + 2 \right)\left( \alpha + 1 \right)\alpha\theta^{3}\\
\mathrm{Var}\left( X \right) &= (\alpha + 1)\alpha\theta^2 - (\alpha\theta)^2 = \alpha\theta^{2}\\
\end{align}\] \[
\]
\[\begin{array}{ll}
\text{Skewness} &= \frac{\mathrm{E}\left\lbrack {(X - \mu_{1}^{\prime})}^{3} \right\rbrack}{{\left( \mathrm{Var}X \right)}^{3/2}} = \frac{\mu_{3}^{\prime} - 3\mu_{2}^{\prime}\mu_{1}^{\prime} + 2{\mu_{1}^{\prime}}^{3}}{{\left(\mathrm{Var} X \right)}^{3/2}} \\
&= \frac{\left( \alpha + 2 \right)\left( \alpha + 1 \right)\alpha\theta^{3} - 3\left( \alpha + 1 \right)\alpha^{2}\theta^{3} + 2\alpha^{3}\theta^{3}}{\left( \alpha\theta^{2} \right)^{3/2}} \\
&= \frac{2}{\alpha^{1/2}} = 1.
\end{array}\]
$$
Maka didapatkan hasil \[
\begin{align}
α&=4\\
E(X)&=αθ=8\\
θ&=2\\
Var(X)&=αθ^2=16
\end{align}
\]
3.1.2 Quantiles
Kuantil dapat digunakan dalam menggambarkan karakteristik distribusi X. Ketika distribusi X kontinu, untuk suatu pecahan tertentu 0≤p≤1 kuantil yang sesuai adalah solusi dari persamaan
\[
\begin{align}
F_X(π_p)=p
\end{align}
\] Sebagai contoh, titik tengah distribusi, \(π_{0.5}\) adalah median. Persentil adalah jenis kuantil; persentil \(100p\) yang merupakan angka sedemikian rupa sehingga \(100×p\) persen dari data berada di bawahnya.
3.1.3 Moment Generating Function
Fungsi pembangkit momen (mgf)dilambangkan dengan MX(t) secara unik mencirikan distribusi dari X . Meskipun ada kemungkinan dua distribusi yang berbeda memiliki momen yang sama namun tetap berbeda, tidak demikian halnya dengan fungsi pembangkit momen. Artinya, jika dua variabel acak memiliki fungsi pembangkit momen yang sama, maka keduanya memiliki distribusi yang sama. Fungsi pembangkit momen diberikan oleh untuk semua nilai t yang memiliki nilai ekspektasi.
\[
\begin{align}
M_X(t) = E(e^{tX}) = \int_{0}^{∞}e^{tX}f_X(x)dx
\end{align}
\] MGF adalah fungsi real yang turunan ke-k pada nol sama dengan momen mentah ke-k dari X . Dalam simbol, ini adalah
\[
\begin{align}
\frac{d^k}{dt^k}MX(t)\Bigr|_{t=0} = E(X^{k})
\end{align}
\]
Example 3.1.3
Variabel acak X memiliki distribusi eksponensial dengan mean \(1/b\) . Maka dapat menacari b jika \(M_{X}\left( - b^{2} \right) = 0.2\) .
\[
\begin{align}
M_{X}(t) = \mathrm{E}\left( e^{tX} \right) = \int_{0}^{\infty}{e^{\text{tx}}be^{- bx} dx} = \int_{0}^{\infty}{be^{- x\left( b - t \right)} dx} = \frac{b}{\left( b - t \right)}.
\end{align}
\] Maka
\[
\begin{align}
M_{X}\left( - b^{2} \right) = \frac{b}{\left( b + b^{2} \right)} = \frac{1}{\left( 1 + b \right)} = 0.2
\end{align}
\]
Maka akan didapatkan
\[
\begin{align}
\frac{1}{\left( 1 + b \right)} &= 0.2\\
1&=0.2(1+b)\\
1&=0.2+0.2b\\
0.8&=0.2b\\
4&=b\\
\end{align}
\] Kita juga dapat menggunakan fungsi pembangkit momen untuk menghitung fungsi pembangkit probabilitas dengan
\[
\begin{align}
P_X(z)= E(z^X)=M_X(logz)
\end{align}
\]
LS0tDQp0aXRsZTogIjMuMSBCYXNpYyBEaXN0cmlidXRpb25hbCBRdWFudGl0aWVzICINCnN1YnRpdGxlOiAiVGVvcmkgUmVzaWtvIg0KYXV0aG9yOiAiQ2xhcmEgRGVsbGEgRXZhbmlhICgyMDIwNDkyMDAxOCkiDQpkYXRlOiAgImByIGZvcm1hdChTeXMuRGF0ZSgpLCAnJUIgJWQsICVZJylgIg0Kb3V0cHV0Og0KICBybWRmb3JtYXRzOjpyb2JvYm9vazogICAjIGh0dHBzOi8vZ2l0aHViLmNvbS9qdWJhL3JtZGZvcm1hdHMNCiAgICBzZWxmX2NvbnRhaW5lZDogdHJ1ZQ0KICAgIHRodW1ibmFpbHM6IHRydWUNCiAgICBsaWdodGJveDogdHJ1ZQ0KICAgIGdhbGxlcnk6IHRydWUNCiAgICBsaWJfZGlyOiBsaWJzDQogICAgZGZfcHJpbnQ6ICJwYWdlZCINCiAgICBjb2RlX2ZvbGRpbmc6ICJzaG93Ig0KICAgIGNvZGVfZG93bmxvYWQ6IHllcw0KICAgIGNzczogInN0eWxlLmNzcyINCg0KLS0tDQoNCg0KYGBge3IgaW5jbHVkZT1GQUxTRX0NCmtuaXRyOjpvcHRzX2NodW5rJHNldChjbGFzcy5zb3VyY2UgPSAibm9jb3B5IiwNCiAgICAgICAgICAgICAgICAgICAgICBjbGFzcy5vdXRwdXQgPSAibm9jb3B5IiwNCiAgICAgICAgICAgICAgICAgICAgICBtZXNzYWdlID0gRiwNCiAgICAgICAgICAgICAgICAgICAgICB3YXJuaW5nID0gRikNCmBgYA0KDQoNCjxicj4NCg0KDQo8aW1nIHN0eWxlPSJmbG9hdDogcmlnaHQ7IG1hcmdpbjogLTUwcHggNTBweCAwcHggNTBweDsgd2lkdGg6MjUlIiBzcmM9Im1lLmpwZWciLz4gDQoNCnwNCjotLS0tIHw6LS0tLQ0KS29udGFrfCAqOiogJFxkb3duYXJyb3ckDQpFbWFpbHwgY2xhcmEuZXZhbmlhQHN0dWRlbnQubWF0YW5hdW5pdmVyc2l0eS5hYy5pZA0KSW5zdGFncmFtIHwgaHR0cHM6Ly93d3cuaW5zdGFncmFtLmNvbS9jbGFyYWV2YW5pYS8gDQpSUHVicyAgfCBodHRwczovL3JwdWJzLmNvbS9jbGFyYWRlbGxhZXZhbmlhLyANCg0KKioNCg0KRGFsYW0gTWF0ZXJpIGluaSBha2FuIG1lbXBlbGFqYXJpOg0KICAtIG1vbWVudHMsDQogIC0gcGVyY2VudGlsZXMsIGFuZA0KICAtIGdlbmVyYXRpbmcgZnVuY3Rpb25zLg0KDQojIDMuMS4xIE1vbWVudHMNCg0KRGVuZ2FuIG1lbWlzYWxrYW4gYFhgIG1lcnVwYWthbiBzdWF0dSB2YXJpYWJlbCBhY2FrIGtvbnRpbnUgZGVuZ2FuIGBwcm9iYWJpbGl0eSBkZW5zaXR5IGZ1bmN0aW9uIChwZGYpYCAkZl9YKHgpJCBkYW4gZnVuZ3NpIGRpdHJpYnVzaSAkRl9YKHgpJC4gDQoNCk1vbWVuIGJha3Uga2UtayBkYXJpIGBYYCB5YW5nIGRpbm90YXNpa2FuIGRlbmdhbiAkzrwnX2skIGFkYWxhaCBuaWxhaSBla3NwZWt0YXNpIGRhcmkgcGFuZ2thdCBrZS1rIGRhcmkgYFhgIHlhbmcgZGloYXJhcGthbiwgYXNhbGthbiBtZW1wdW55YWkgbmlsYWkuIE1vbWVuIHBlcnRhbWEgJM68J18xJCBhZGFsYWggbmlsYWkgdGVuZ2FoIChtZWFuKSBkYXJpIGBYYCB5YW5nIGJpYXNhbnlhIGRpbGFtYmFuZ2thbiBkZW5nYW4gYM68YCAuIFJ1bXVzIHVudHVrICTOvCdfayQgYWRhbGFoIHNlYmFnYWkgYmVyaWt1dCA6DQoNCiQkDQpcYmVnaW57YWxpZ259DQrOvCdfayA9IEUoWF5rKSA9IFxpbnRfezB9XnviiJ59eF5rZl9YKHgpZHgNClxlbmR7YWxpZ259DQokJA0KDQpVbnR1ayBtZW5kdWt1bmcgZGFyaSB2YXJpYWJlbCBhY2FrIGBYYCBkYXBhdCBkaWFzdW1zaWthbiBzZWJhZ2FpIG5vbm5lZ2F0aWYgZGlrYXJlbmFrYW4gcGFkYSBrZW55YXRhYW5ueWEgamFyYW5nIHNla2FsaSBiZXJuaWxhaSBuZWdhdGlmLiBTYWxhaCBzYXR1IGNvbnRvaG55YSB5YW5nIG1lbnVuanVra2FuIGJhaHdhIHJhdyBtb21lbnRzIHVudHVrIHZhcmlhYmVsIG5vbm5lZ2F0aWYgZGFwYXQgZGloaXR1bmcgbWVuZ2d1bmFrYW4gcnVtdXMgOg0KDQokJA0KXGJlZ2lue2FsaWdufQ0KzrwnX2sgPSAgXGludF97MH1ee+KInn1reF57ay0xfVsxLUZfWCh4KV1keCwNClxlbmR7YWxpZ259DQokJA0KDQp5YW5nIGRpZGFzYXJrYW4gcGFkYSBmdW5nc2kgc3Vydml2YWwgJFNfWCh4KT0x4oiSRl9YKHgpJC4gUGFkYSBydW11cyBpbmkgYmVyZ3VuYSBkaXNhYXQgaz0xLg0KDQpDZW50cmFsIG1vbWVudCBrZSAtayBkYXJpIGBYYCB5YW5nIGRpbm90YXNpa2FuIGRlbmdhbiAkzrxfayQgeWFuZyBtZXJ1cGFrYW4gbmlsYWkgeWFuZyBkaWhhcmFwa2FuIGRhcmkgcGFuZ2thdCBrZS1rIGRhcmkgZGV2aWFzaSB4IGRhbiBkYXJpIG1lYW4gYM68YC4gTWFrYSB1bnR1ayBydW11cyAkzrxfayQgYWRhbGFoIA0KDQokJA0KXGJlZ2lue2FsaWdufQ0KzrxfayA9IEVbKFgtzrwpXmtdID0gXGludF97MH1ee+KInn0oeC3OvClea2ZfWCh4KWR4DQpcZW5ke2FsaWdufQ0KJCQNCg0KTW9tZW4gcHVzYXQga2VkdWEgJM68XzIkIG1lbmRlZmluaXNpa2FuIHZhcmlhbnMgZGFyaSBgWGAgeWFuZyBkaW5vdGFzaWthbiBkZW5nYW4gJM+DXjIkLiBBa2FyIGt1YWRyYXQgZGFyaSB2YXJpYW5zIGFkYWxhaCBzaW1wYW5nYW4gYmFrdSAkz4MkIC4gDQoNClJhc2lvIG1vbWVuIHNlbnRyYWwga2V0aWdhIHRlcmhhZGFwIHBhbmdrYXQgdGlnYSBkYXJpIGRldmlhc2kgc3RhbmRhciAkKM68XzMvz4NeMykkIHlhbmcgbWVuZGVmaW5pc2lrYW4ga29lZmlzaWVuIGtlbWlyaW5nYW4geWFuZyBtZXJ1cGFrYW4gdWt1cmFuIHNpbWV0cmkuIEtvZWZpc2llbiBrZW1lbmNlbmdhbiB5YW5nIHBvc2l0aWYgbWVudW5qdWtrYW4gYmFod2EgZGlzdHJpYnVzaSBjb25kb25nIGtlIGthbmFuIChjb25kb25nIGtlIGthbmFuKS4gDQoNClJhc2lvIG1vbWVuIHNlbnRyYWwga2VlbXBhdCBkZW5nYW4gcGFuZ2thdCBlbXBhdCBkYXJpIGRldmlhc2kgc3RhbmRhciAkKM68XzQvz4NeNCkkIG1lbmRlZmluaXNpa2FuIGtvZWZpc2llbiBrdXJ0b3Npcy4gRGlzdHJpYnVzaSBub3JtYWwgbWVtaWxpa2kga29lZmlzaWVuIGt1cnRvc2lzIDMuIERpc3RyaWJ1c2kgZGVuZ2FuIGtvZWZpc2llbiBrdXJ0b3NpcyBsZWJpaCBiZXNhciBkYXJpIDMgbWVtaWxpa2kgZWtvciB5YW5nIGxlYmloIGJlcmF0IGRhcmlwYWRhIGRpc3RyaWJ1c2kgbm9ybWFsLCBzZWRhbmdrYW4gZGlzdHJpYnVzaSBkZW5nYW4ga29lZmlzaWVuIGt1cnRvc2lzIGt1cmFuZyBkYXJpIDMgbWVtaWxpa2kgZWtvciB5YW5nIGxlYmloIHJpbmdhbiBkYW4gbGViaWggZGF0YXIuIA0KDQojIyBFeGFtcGxlIDMuMS4xIA0KTWVuZ2FzdW1zaWthbiBiYWh3YSB2YXJpYWJlbCBhY2FrIFggbWVtaWxpa2kgZGlzdHJpYnVzaSBnYW1tYSBkZW5nYW4gcmF0YS1yYXRhIDkgZGFuIHNrZXduZXNzIDEuIE1ha2EgZGljYXJpIHBlcmhpdHVuZ2FuIHZhcmlhbnMgZGFyaSBYLg0KDQpwZGYgZGFyaSBYIGFkYWxhaCA6DQoNCiQkDQpcYmVnaW57YWxpZ259DQpmX3tYfVxsZWZ0KCB4IFxyaWdodCkgPSBcZnJhY3tcbGVmdCggeCAvIFx0aGV0YSBccmlnaHQpXntcYWxwaGF9fXt4IH5cR2FtbWFcbGVmdCggXGFscGhhIFxyaWdodCl9IGVeey0geCAvIFx0aGV0YX0NClxlbmR7YWxpZ259DQokJA0KVW50dWsgJHg+MCQsIGRhbiAkzrE+MCQsIG1ha2EgbW9tZW4gYmFrdSBrZS1rIGFkYWxhaA0KDQokJA0KXGJlZ2lue2FsaWdufQ0KXG11X3trfV57XHByaW1lfSA9IFxtYXRocm17RX1cbGVmdCggWF57a30gXHJpZ2h0KSA9IFxpbnRfezB9XntcaW5mdHl9e1xmcmFjezF9e1xHYW1tYVxsZWZ0KCBcYWxwaGEgXHJpZ2h0KVx0aGV0YV57XGFscGhhfX14XntrICsgXGFscGhhIC0gMX1lXnstIHggLyBcdGhldGF9IGR4fSA9IFxmcmFje1xHYW1tYVxsZWZ0KCBrICsgXGFscGhhIFxyaWdodCl9e1xHYW1tYVxsZWZ0KCBcYWxwaGEgXHJpZ2h0KX1cdGhldGFee2t9DQpcZW5ke2FsaWdufQ0KJCQNCg0KRGVuZ2FuIG1lbWJlcmlrYW4gJFxHYW1tYVxsZWZ0KCByICsgMSBccmlnaHQpID0gclxHYW1tYVxsZWZ0KCByIFxyaWdodCkkIGRhbiAkXEdhbW1hXGxlZnQoIDEgXHJpZ2h0KSA9IDEkLiBNYWthICRcbXVfezF9XntccHJpbWV9ID0gXG1hdGhybXtFfVxsZWZ0KCBYIFxyaWdodCkgPSBcYWxwaGFcdGhldGEkDQoNCg0KJCQNClxiZWdpbnthbGlnbn0NClxtdV97Mn1ee1xwcmltZX0gJj0gXG1hdGhybXtFfVxsZWZ0KCBYXnsyfSBccmlnaHQpID0gXGxlZnQoIFxhbHBoYSArIDEgXHJpZ2h0KVxhbHBoYVx0aGV0YV57Mn0gPQ0KXG11X3szfV57XHByaW1lfSA9IFxtYXRocm17RX1cbGVmdCggWF57M30gXHJpZ2h0KSA9IFxsZWZ0KCBcYWxwaGEgKyAyIFxyaWdodClcbGVmdCggXGFscGhhICsgMSBccmlnaHQpXGFscGhhXHRoZXRhXnszfVxcDQoNClxtYXRocm17VmFyfVxsZWZ0KCBYIFxyaWdodCkgJj0gKFxhbHBoYSArIDEpXGFscGhhXHRoZXRhXjIgLSAoXGFscGhhXHRoZXRhKV4yID0gXGFscGhhXHRoZXRhXnsyfVxcDQoNClxlbmR7YWxpZ259DQokJA0KJCRcYmVnaW57YXJyYXl9e2xsfQ0KXHRleHR7U2tld25lc3N9ICAmPSBcZnJhY3tcbWF0aHJte0V9XGxlZnRcbGJyYWNrIHsoWCAtIFxtdV97MX1ee1xwcmltZX0pfV57M30gXHJpZ2h0XHJicmFja317e1xsZWZ0KCBcbWF0aHJte1Zhcn1YIFxyaWdodCl9XnszLzJ9fSA9IFxmcmFje1xtdV97M31ee1xwcmltZX0gLSAzXG11X3syfV57XHByaW1lfVxtdV97MX1ee1xwcmltZX0gKyAye1xtdV97MX1ee1xwcmltZX19XnszfX17e1xsZWZ0KFxtYXRocm17VmFyfSBYIFxyaWdodCl9XnszLzJ9fSBcXA0KICY9IFxmcmFje1xsZWZ0KCBcYWxwaGEgKyAyIFxyaWdodClcbGVmdCggXGFscGhhICsgMSBccmlnaHQpXGFscGhhXHRoZXRhXnszfSAtIDNcbGVmdCggXGFscGhhICsgMSBccmlnaHQpXGFscGhhXnsyfVx0aGV0YV57M30gKyAyXGFscGhhXnszfVx0aGV0YV57M319e1xsZWZ0KCBcYWxwaGFcdGhldGFeezJ9IFxyaWdodCleezMvMn19IFxcDQogJj0gXGZyYWN7Mn17XGFscGhhXnsxLzJ9fSA9IDEuDQogXGVuZHthcnJheX0NCiAkJA0KDQpNYWthIGRpZGFwYXRrYW4gaGFzaWwgDQokJA0KXGJlZ2lue2FsaWdufQ0KzrEmPTRcXA0KRShYKSY9zrHOuD04XFwNCs64Jj0yXFwNClZhcihYKSY9zrHOuF4yPTE2DQpcZW5ke2FsaWdufQ0KJCQNCg0KDQojIDMuMS4yIFF1YW50aWxlcw0KDQpLdWFudGlsIGRhcGF0IGRpZ3VuYWthbiBkYWxhbSBtZW5nZ2FtYmFya2FuIGthcmFrdGVyaXN0aWsgZGlzdHJpYnVzaSBgWGAuIEtldGlrYSBkaXN0cmlidXNpIGBYYCBrb250aW51LCB1bnR1ayBzdWF0dSBwZWNhaGFuIHRlcnRlbnR1IGAw4omkcOKJpDFgIGt1YW50aWwgeWFuZyBzZXN1YWkgYWRhbGFoIHNvbHVzaSBkYXJpIHBlcnNhbWFhbiANCg0KJCQNClxiZWdpbnthbGlnbn0NCkZfWCjPgF9wKT1wDQpcZW5ke2FsaWdufQ0KJCQNClNlYmFnYWkgY29udG9oLCB0aXRpayB0ZW5nYWggZGlzdHJpYnVzaSwgJM+AX3swLjV9JCBhZGFsYWggbWVkaWFuLiBQZXJzZW50aWwgYWRhbGFoIGplbmlzIGt1YW50aWw7IHBlcnNlbnRpbCAkMTAwcCQgeWFuZyBtZXJ1cGFrYW4gYW5na2Egc2VkZW1pa2lhbiBydXBhIHNlaGluZ2dhICQxMDDDl3AkIHBlcnNlbiBkYXJpIGRhdGEgYmVyYWRhIGRpIGJhd2FobnlhLg0KDQojIDMuMS4zIE1vbWVudCBHZW5lcmF0aW5nIEZ1bmN0aW9uDQoNCkZ1bmdzaSBwZW1iYW5na2l0IG1vbWVuIChtZ2YpZGlsYW1iYW5na2FuIGRlbmdhbiBgTVgodClgIHNlY2FyYSB1bmlrIG1lbmNpcmlrYW4gZGlzdHJpYnVzaSBkYXJpIGBYYCAuIE1lc2tpcHVuIGFkYSBrZW11bmdraW5hbiBkdWEgZGlzdHJpYnVzaSB5YW5nIGJlcmJlZGEgbWVtaWxpa2kgbW9tZW4geWFuZyBzYW1hIG5hbXVuIHRldGFwIGJlcmJlZGEsIHRpZGFrIGRlbWlraWFuIGhhbG55YSBkZW5nYW4gZnVuZ3NpIHBlbWJhbmdraXQgbW9tZW4uIEFydGlueWEsIGppa2EgZHVhIHZhcmlhYmVsIGFjYWsgbWVtaWxpa2kgZnVuZ3NpIHBlbWJhbmdraXQgbW9tZW4geWFuZyBzYW1hLCBtYWthIGtlZHVhbnlhIG1lbWlsaWtpIGRpc3RyaWJ1c2kgeWFuZyBzYW1hLiBGdW5nc2kgcGVtYmFuZ2tpdCBtb21lbiBkaWJlcmlrYW4gb2xlaCB1bnR1ayBzZW11YSBuaWxhaSB0IHlhbmcgbWVtaWxpa2kgbmlsYWkgZWtzcGVrdGFzaS4NCg0KJCQNClxiZWdpbnthbGlnbn0NCk1fWCh0KSA9IEUoZV57dFh9KSA9IFxpbnRfezB9XnviiJ59ZV57dFh9Zl9YKHgpZHgNClxlbmR7YWxpZ259DQokJA0KTUdGIGFkYWxhaCBmdW5nc2kgcmVhbCB5YW5nIHR1cnVuYW4ga2UtayBwYWRhIG5vbCBzYW1hIGRlbmdhbiBtb21lbiBtZW50YWgga2UtayBkYXJpIFggLiBEYWxhbSBzaW1ib2wsIGluaSBhZGFsYWggDQoNCiQkDQpcYmVnaW57YWxpZ259DQpcZnJhY3tkXmt9e2R0Xmt9TVgodClcQmlncnxfe3Q9MH0gPSBFKFhee2t9KQ0KXGVuZHthbGlnbn0NCiQkDQoNCiMjIEV4YW1wbGUgMy4xLjMNClZhcmlhYmVsIGFjYWsgWCBtZW1pbGlraSBkaXN0cmlidXNpIGVrc3BvbmVuc2lhbCBkZW5nYW4gbWVhbiAkMS9iJC4gTWFrYSBkYXBhdCBtZW5hY2FyaSBiIGppa2EgJE1fe1h9XGxlZnQoIC0gYl57Mn0gXHJpZ2h0KSA9IDAuMiQuDQoNCiQkDQpcYmVnaW57YWxpZ259DQpNX3tYfSh0KSA9IFxtYXRocm17RX1cbGVmdCggZV57dFh9IFxyaWdodCkgPSBcaW50X3swfV57XGluZnR5fXtlXntcdGV4dHt0eH19YmVeey0gYnh9IGR4fSA9IFxpbnRfezB9XntcaW5mdHl9e2JlXnstIHhcbGVmdCggYiAtIHQgXHJpZ2h0KX0gZHh9ID0gXGZyYWN7Yn17XGxlZnQoIGIgLSB0IFxyaWdodCl9Lg0KXGVuZHthbGlnbn0NCiQkDQpNYWthIA0KDQokJA0KXGJlZ2lue2FsaWdufQ0KTV97WH1cbGVmdCggLSBiXnsyfSBccmlnaHQpID0gXGZyYWN7Yn17XGxlZnQoIGIgKyBiXnsyfSBccmlnaHQpfSA9IFxmcmFjezF9e1xsZWZ0KCAxICsgYiBccmlnaHQpfSA9IDAuMg0KXGVuZHthbGlnbn0NCiQkDQoNCk1ha2EgYWthbiBkaWRhcGF0a2FuIA0KDQokJA0KXGJlZ2lue2FsaWdufQ0KIFxmcmFjezF9e1xsZWZ0KCAxICsgYiBccmlnaHQpfSAmPSAwLjJcXA0KIDEmPTAuMigxK2IpXFwNCiAxJj0wLjIrMC4yYlxcDQogMC44Jj0wLjJiXFwNCiA0Jj1iXFwNClxlbmR7YWxpZ259DQokJA0KS2l0YSBqdWdhIGRhcGF0IG1lbmdndW5ha2FuIGZ1bmdzaSBwZW1iYW5na2l0IG1vbWVuIHVudHVrIG1lbmdoaXR1bmcgZnVuZ3NpIHBlbWJhbmdraXQgcHJvYmFiaWxpdGFzIGRlbmdhbg0KDQokJA0KXGJlZ2lue2FsaWdufQ0KUF9YKHopPSBFKHpeWCk9TV9YKGxvZ3opDQpcZW5ke2FsaWdufQ0KJCQNCg0KDQojIFJlZmVyZW5zaSANCi0gaHR0cHM6Ly9vcGVuYWN0dGV4dHMuZ2l0aHViLmlvL0xvc3MtRGF0YS1BbmFseXRpY3MvQ2hhcEZyZXF1ZW5jeS1Nb2RlbGluZy5odG1sI1M6Z29vZG5lc3Mtb2YtZml0
