Problem 8.2 (4 points)

We are given \(E[\hat{\theta}_1] = E[\hat{\theta}_2] = \theta\), \(V[\hat{\theta}_1]=\sigma_1^2\), \(V[\hat{\theta}_2]=\sigma_2^2\), and \(\hat{\theta}_1\) and \(\hat{\theta}_2\) are independent.

Let \(\hat{\theta}_3 = a\cdot \hat{\theta}_1 + (1-a)\cdot \hat{\theta}_2\).
  1. Since \(E[\hat{\theta}_3] = E[a\cdot \hat{\theta}_1 + (1-a)\cdot \hat{\theta}_2] = a\cdot E[\hat{\theta}_1] + (1-a)\cdot E[\hat{\theta}_2] = a\cdot \theta + (1-a)\cdot \theta = \theta\), we see that \(\hat{\theta}_3\) is an unbiased estimator for \(\theta\) for all \(a\in\mathbb{R}\).
  2. The variance of \(\hat{\theta}_3\) is given by the quadratic function \[ f(a) = V[\hat{\theta}_3] = a^2\cdot \sigma_1^2 + (1 - a)^2\cdot \sigma_2^2. \] The first derivative \(f'(a) = 2\cdot [(\sigma_1^2 + \sigma_2^2)\cdot a - \sigma_2^2]\) equals zero at \(a = \sigma_2^2 / (\sigma_1^2 + \sigma_2^2)\), and since the second derivative \(f''(a) = 2\cdot (\sigma_1^2 + \sigma_2^2) > 0\), we see that the critical point is a global minimum for the quadratic function. Thus, of all possible weighted averages of the unbiased estimators \(\hat{\theta}_1\), \(\hat{\theta}_2\), \[ \hat{\theta}_3 = \frac{\sigma_2^2}{\sigma_1^2 + \sigma_2^2}\cdot \hat{\theta}_1 + \frac{\sigma_1^2}{\sigma_1^2 + \sigma_2^2}\cdot \hat{\theta}_2 \] has the smallest variance \(\sigma_3^2 = V[\hat{\theta}_3] = \sigma_1^2\cdot \sigma_2^2 / (\sigma_1^2 + \sigma_2^2)\). Observe that the variance of \(\hat{\theta}_3\) is less than the variance of either \(\hat{\theta}_1\) or \(\hat{\theta}_2\). This fact becomes apparant when we see that \(1/\sigma_3^2 = 1/\sigma_1^2 + 1/\sigma_2^2\) and noting that both terms in the sum are positive.

    This problem shows that a suitable weighted average of independant, unbiased estimators for \(\theta\) produces a new unbiased estimator for \(\theta\) that has smaller variance than either of the original estimators.

Problem 8.4 (4 points)

Suppose \(Y_1, Y_2, Y_3\) denotes a random sample from an exponential distribution with probability density function \[ f(y) = \begin{cases} \frac{1}{\theta}\cdot e^{-y/\theta} &, y > 0 \\ 0 &, y \le 0. \end{cases} \] Note that \(E[Y_i]=\theta\) and \(V[Y_i]=\theta^2\), \(i=1,2,3\).

  1. Consider: \(\hat{\theta}_1 = Y_1\), \(\hat{\theta}_2 = \frac{1}{2}(Y_1 + Y_2)\), \(\hat{\theta}_3 = \frac{1}{3}Y_1 + \frac{2}{3}Y_2\) \(\hat{\theta}_4 = \min{(Y_1, Y_2, Y_3)}\), \(\hat{\theta}_5 = \overline{Y}\).

    \[\begin{align*} E[\hat{\theta}_1] &= E[Y_1] = \theta. \\ E[\hat{\theta}_2] &= \frac{1}{2}\left(E[Y_1] + E[Y_2]\right) = \frac{1}{2}(\theta + \theta) = \theta. \\ E[\hat{\theta}_3] &= \frac{1}{3}E[Y_1] + \frac{2}{3}E[Y_2] = \frac{1}{3}\theta + \frac{2}{3}\theta = \theta. \\ E[\hat{\theta}_4] &= \frac{1}{3}\theta \text{ (see exercise 6.65(a), page 324).} \\ E[\hat{\theta}_5] &= E[\overline{Y}] = \theta. \end{align*}\] So, \(\hat{\theta}_1\), \(\hat{\theta}_2\), \(\hat{\theta}_3\), \(\hat{\theta}_5\) are unbiased estimators for \(\theta\).
  2. We calculate the variances of the unbiased estimators. \[\begin{align*} V[\hat{\theta}_1] &= V[Y_1] = \theta^2 \\ V[\hat{\theta}_2] &= \frac{1}{4}(\theta^2 + \theta^2) = \frac{1}{2}\theta^2 \\ V[\hat{\theta}_3] &= \frac{1}{9}\theta^2 + \frac{4}{9}\theta^2 = \frac{5}{9}\theta^2 \\ V[\hat{\theta}_5] &= \frac{1}{3}\theta^2, \end{align*}\] so \(\hat{\theta}_5 = \overline{Y}\) has the least variance.

Problem 8.8 (6 points)

We are given that \(Y_1, \dotsc, Y_n\) represents a random sample from a uniform distribution on the interval \((\theta,\theta + 1)\). Therefore, \(E[Y_i] = \theta + \frac{1}{2}\), \(V[Y_i] = 1/12\), \(i=1,\dotsc,n\).
  1. The bias for \(\overline{Y}\) (as an estimator for \(\theta\)) is \[ B = E[\overline{Y}] - \theta = (\theta + 1/2) - \theta = 1/2. \]
  2. Let \(\hat{\theta} = \overline{Y} - 1/2\), then \(E[\hat{\theta}] = E[\overline{Y} - 1/2] = (\theta + 1/2) - 1/2 = \theta\), so \(\overline{Y} - 1/2\) is an unbiased estimator for \(\theta\).
  3. \(\text{MSE}[\overline{Y}] = V[\overline{Y}] + B^2 = \frac{1}{12n} + \frac{1}{4} = \frac{3n+1}{12n}\).

Problem 8.10 (6 points)

Let \(Y_1, \dotsc, Y_n\) denote a random sample from a population whose density and distribution functions are \[\begin{align*} f(y) &= \begin{cases} \alpha\cdot y^{\alpha-1}/\theta^\alpha & , 0 \le y \le \theta \\ 0 & , \text{ otherwise} \end{cases} \\ F(y) &= \begin{cases} 1 & , y > \theta \\ (y/\theta)^\alpha & , 0 \le y \le \theta \\ 0 & , y < 0. \end{cases} \end{align*}\] Let \(\hat{\theta} = \max\{Y_1, \dotsc, Y_n\} = Y_{(n)}\). Then \(\hat{\theta}\) has density \[ g_{(n)}(y) = \begin{cases} n \alpha\cdot y^{n \alpha - 1} / \theta^{n \alpha} & , 0 \le y \le \theta \\ 0 & , \text{ otherwise.} \end{cases} \]
  1. We use \(g_{(n)}\) to find the expected value of \(\hat{\theta}\). \[ E[\hat{\theta}] = \int_{0}^{\theta} n\alpha\cdot y\cdot y^{n\alpha - 1}/\theta^{n\alpha}\,dy = \left.\frac{n\alpha}{n\alpha + 1}\cdot \frac{y^{n\alpha + 1}}{\theta^{n\alpha}}\right|_0^{\theta} = \frac{n\alpha}{n\alpha + 1}\cdot \theta. \] Since \(E[\hat{\theta}] \ne \theta\), we see that \(\hat{\theta}\) is a biased estimator for \(\theta\).
  2. \(\frac{n\alpha + 1}{n\alpha}\cdot \hat{\theta}\) is an unbiased estimator for \(\theta\).
  3. The bias for \(\hat{\theta}\) is \(B = \frac{n\alpha}{n\alpha+1}\cdot\theta - \theta = -\theta / (n\alpha + 1)\). We calculate the second moment \[ E\left[\hat{\theta}^2\right] = \int_{0}^{\theta} n\alpha\cdot y^{n\alpha + 1} / \theta^{n\alpha}\,dy = \frac{n\alpha}{n\alpha + 2}\cdot \theta^2. \] Finally, we calculate the mean square error \[\begin{align*} \text{MSE}[\hat{\theta}] &= V[\hat{\theta}] + B^2 \\ &= \frac{n\alpha}{n\alpha + 2}\cdot \theta^2 - \left(\frac{n\alpha}{n\alpha + 1}\right)^2\cdot \theta^2 + \frac{1}{(n\alpha + 1)^2}\cdot \theta^2 \\ &= \frac{2\cdot \theta^2}{(n\alpha + 1)(n\alpha + 2)}. \end{align*}\]

Problem 8.12 (6 points)

Suppose \(Y_1, \dotsc, Y_n\) constitute a random sample from a normal distribution with mean \(\mu\) and variance \(\sigma^2\).

Let \(Y = (n - 1)\cdot S^2 / \sigma^2\), which we know has the \(\chi^2\) probability distribution with \(\nu = n - 1\) degrees of freedom. By exercise 4.90(c), \(E[\sqrt{Y}] = \sqrt{2}\cdot \Gamma(n/2) / \Gamma((n-1)/2)\).
  1. Now, \(S = {\sigma}\cdot \sqrt{Y} / {\sqrt{n-1}}\) and therefore \[ E[S] = \frac{\sigma}{\sqrt{n-1}}\cdot E\left[\sqrt{Y}\right] = \sqrt{\frac{2}{n-1}}\cdot \frac{\Gamma(n/2)}{\Gamma\left((n-1)/2)\right)}\cdot \sigma. \] We see that \(S\) is a biased estimator for \(\sigma\), even though \(S^2\) is an unbiased estimator for \(\sigma^2\).
  2. We can adjust \(S\) by a multiplicative constant to create the unbiased estimator for \(\sigma\) \[ \hat{\sigma} = \sqrt{\frac{n-1}{2}}\cdot \frac{\Gamma\left((n-1)/2)\right)}{\Gamma(n/2)}\cdot S. \]
  3. An unbiased estimator for \(\mu - z_{\alpha}\cdot \sigma\) is \(\overline{Y} - z_{\alpha}\cdot \hat{\sigma}\), where \(\overline{Y} = \frac{1}{n}\sum_{i=1}^{n} Y_{i}\) and \(\hat{\sigma}\) is the estimator defined in part (b).

Problem 8.14 (2 points)

Let \(Y_1, \dotsc, Y_n\) denote a random sample from a population with a uniform distribution on the interval \((0,\theta)\). Let \(Y_{(1)} = \min\{Y_1, \dotsc, Y_n\}\). Recall that the density and distribution functions for this uniform distribution are, respectively, \[\begin{align*} f(y) &= \begin{cases} \frac{1}{\theta} & , 0 < y < \theta \\ 0 & , \text{ otherwise}, \end{cases} \\ F(y) &= \begin{cases} 1 & , y \ge \theta \\ \frac{y}{\theta} & , 0 < y < \theta \\ 0 & , y \le 0. \end{cases} \end{align*}\]

Thus, the density function for \(Y_{(1)}\) is \[ g_{(1)}(y) = n\cdot \left[ 1 - F(y)\right]^{n-1}\cdot f(y) = \begin{cases} \frac{n}{\theta} \cdot \left( 1 - \frac{y}{\theta}\right)^{n-1} & , 0 < y < \theta \\ 0 & , \text{ otherwise} \end{cases}. \] Now that we have the density for \(Y_{(1)}\), we can calculate \[ E[Y_{(1)}] = \int_{0}^{\theta} \frac{ny}{\theta}\cdot \left( 1 - \frac{y}{\theta}\right)^{n-1}\,dy = \frac{\theta}{(n + 1)}. \] We see that \(Y_{(1)}\) is a biased estimator for \(\theta\), but we may adjust it by a multiplicative constant to obtain the unbiased estimator \(\hat{\theta} = (n + 1)\cdot Y_{(1)}\).

Problem 8.16 (4 points)

Suppose that \(Y_1, Y_2, Y_3, Y_4\) denotes a random sample from a population with an exponential distribution whose density is given by \[ f(y) = \begin{cases} \frac{1}{\theta}\cdot e^{-y/\theta} & , y > 0 \\ 0 & , \text{ otherwise} \end{cases}. \] Recall that \(Y\) is a Gamma distribution with parameters \(\alpha = 1\) and \(\beta = \theta\). By exercise 4.89(a), \(E[Y^a] = \theta^{a}\cdot \Gamma(1 + a)/\Gamma(1) = \theta^{a}\cdot a\cdot \Gamma(a)\), therefore \[ E[ \sqrt{Y} ] = \frac{1}{2} \Gamma\left(\frac{1}{2}\right)\cdot \sqrt{\theta} = \frac{1}{2}\sqrt{\pi\cdot \theta}. \]
  1. Let \(X = \sqrt{Y_1\cdot Y_2}\). Since \(Y_1\) and \(Y_2\) are independent, \[ E[X] = E[Y_1]\cdot E[Y_2] = \frac{1}{4}\pi\cdot \theta, \] so an unbiased estimator for \(\theta\) is \(\frac{4}{\pi}\sqrt{Y_1\cdot Y_2}\).
  2. Let \(W = \sqrt{Y_1\cdot Y_2\cdot Y_3\cdot Y_4}\), then \[ E[W] = E[Y_1]\cdot E[Y_2]\cdot E[Y_3]\cdot E[Y_4] = \frac{1}{16}\pi^2\cdot \theta^2. \] Thus, an unbiased estimator for \(\theta^2\) is \(\frac{16}{\pi^2}\sqrt{Y_1\cdot Y_2\cdot Y_3\cdot Y_4}\).
LS0tCnRpdGxlOiAiU29sdXRpb25zIHRvIEhvbWV3b3JrIEFzc2lnbm1lbnQgMyIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKIyMgUHJvYmxlbSA4LjIgKDQgcG9pbnRzKQpXZSBhcmUgZ2l2ZW4gXChFW1xoYXR7XHRoZXRhfV8xXSA9IEVbXGhhdHtcdGhldGF9XzJdID0gXHRoZXRhXCksIFwoVltcaGF0e1x0aGV0YX1fMV09XHNpZ21hXzFeMlwpLCBcKFZbXGhhdHtcdGhldGF9XzJdPVxzaWdtYV8yXjJcKSwgYW5kIFwoXGhhdHtcdGhldGF9XzFcKSBhbmQgXChcaGF0e1x0aGV0YX1fMlwpIGFyZSBpbmRlcGVuZGVudC4KCkxldCBcKFxoYXR7XHRoZXRhfV8zID0gYVxjZG90IFxoYXR7XHRoZXRhfV8xICsgKDEtYSlcY2RvdCBcaGF0e1x0aGV0YX1fMlwpLgo8b2wgdHlwZT0iYSI+CjxsaT4KU2luY2UgXChFW1xoYXR7XHRoZXRhfV8zXSA9IEVbYVxjZG90IFxoYXR7XHRoZXRhfV8xICsgKDEtYSlcY2RvdCBcaGF0e1x0aGV0YX1fMl0KPSBhXGNkb3QgRVtcaGF0e1x0aGV0YX1fMV0gKyAoMS1hKVxjZG90IEVbXGhhdHtcdGhldGF9XzJdID0gYVxjZG90IFx0aGV0YSArICgxLWEpXGNkb3QgXHRoZXRhIAo9IFx0aGV0YVwpLCB3ZSBzZWUgdGhhdCBcKFxoYXR7XHRoZXRhfV8zXCkgaXMgYW4gdW5iaWFzZWQgZXN0aW1hdG9yIGZvciBcKFx0aGV0YVwpIGZvciBhbGwgXChhXGluXG1hdGhiYntSfVwpLgo8L2xpPgo8bGk+ClRoZSB2YXJpYW5jZSBvZiBcKFxoYXR7XHRoZXRhfV8zXCkgaXMgZ2l2ZW4gYnkgdGhlIHF1YWRyYXRpYyBmdW5jdGlvbgpcWwpmKGEpID0gVltcaGF0e1x0aGV0YX1fM10gPSBhXjJcY2RvdCBcc2lnbWFfMV4yICsgKDEgLSBhKV4yXGNkb3QgXHNpZ21hXzJeMi4KXF0KVGhlIGZpcnN0IGRlcml2YXRpdmUgXChmJyhhKSA9IDJcY2RvdCBbKFxzaWdtYV8xXjIgKyBcc2lnbWFfMl4yKVxjZG90IGEgLSBcc2lnbWFfMl4yXVwpIGVxdWFscyB6ZXJvIGF0IFwoYSA9IFxzaWdtYV8yXjIgLyAoXHNpZ21hXzFeMiArIFxzaWdtYV8yXjIpXCksIGFuZCBzaW5jZSB0aGUgc2Vjb25kIGRlcml2YXRpdmUgXChmJycoYSkgPSAyXGNkb3QgKFxzaWdtYV8xXjIgKyBcc2lnbWFfMl4yKSA+IDBcKSwgd2Ugc2VlIHRoYXQgdGhlIGNyaXRpY2FsIHBvaW50IGlzIGEgZ2xvYmFsIG1pbmltdW0gZm9yIHRoZSBxdWFkcmF0aWMgZnVuY3Rpb24uIFRodXMsIG9mIGFsbCBwb3NzaWJsZSB3ZWlnaHRlZCBhdmVyYWdlcyBvZiB0aGUgdW5iaWFzZWQgZXN0aW1hdG9ycyBcKFxoYXR7XHRoZXRhfV8xXCksIFwoXGhhdHtcdGhldGF9XzJcKSwgClxbClxoYXR7XHRoZXRhfV8zID0gXGZyYWN7XHNpZ21hXzJeMn17XHNpZ21hXzFeMiArIFxzaWdtYV8yXjJ9XGNkb3QgXGhhdHtcdGhldGF9XzEKICArIFxmcmFje1xzaWdtYV8xXjJ9e1xzaWdtYV8xXjIgKyBcc2lnbWFfMl4yfVxjZG90IFxoYXR7XHRoZXRhfV8yClxdCmhhcyB0aGUgc21hbGxlc3QgdmFyaWFuY2UgXChcc2lnbWFfM14yID0gVltcaGF0e1x0aGV0YX1fM10gPSBcc2lnbWFfMV4yXGNkb3QgXHNpZ21hXzJeMiAvIChcc2lnbWFfMV4yICsgXHNpZ21hXzJeMilcKS4gT2JzZXJ2ZSB0aGF0IHRoZSB2YXJpYW5jZSBvZiBcKFxoYXR7XHRoZXRhfV8zXCkgaXMgbGVzcyB0aGFuIHRoZSB2YXJpYW5jZSBvZiBlaXRoZXIgXChcaGF0e1x0aGV0YX1fMVwpIG9yIFwoXGhhdHtcdGhldGF9XzJcKS4gVGhpcyBmYWN0IGJlY29tZXMgYXBwYXJhbnQgd2hlbiB3ZSBzZWUgdGhhdCBcKDEvXHNpZ21hXzNeMiA9IDEvXHNpZ21hXzFeMiArIDEvXHNpZ21hXzJeMlwpIGFuZCBub3RpbmcgdGhhdCBib3RoIHRlcm1zIGluIHRoZSBzdW0gYXJlIHBvc2l0aXZlLiAKClRoaXMgcHJvYmxlbSBzaG93cyB0aGF0IGEgc3VpdGFibGUgd2VpZ2h0ZWQgYXZlcmFnZSBvZiBpbmRlcGVuZGFudCwgdW5iaWFzZWQgZXN0aW1hdG9ycyBmb3IgXChcdGhldGFcKSBwcm9kdWNlcyBhIG5ldyB1bmJpYXNlZCBlc3RpbWF0b3IgZm9yIFwoXHRoZXRhXCkgdGhhdCBoYXMgc21hbGxlciB2YXJpYW5jZSB0aGFuIGVpdGhlciBvZiB0aGUgb3JpZ2luYWwgZXN0aW1hdG9ycy4gCgojI1Byb2JsZW0gOC40ICg0IHBvaW50cykKU3VwcG9zZSBcKFlfMSwgWV8yLCBZXzNcKSBkZW5vdGVzIGEgcmFuZG9tIHNhbXBsZSBmcm9tIGFuIGV4cG9uZW50aWFsIGRpc3RyaWJ1dGlvbiB3aXRoIHByb2JhYmlsaXR5IGRlbnNpdHkgZnVuY3Rpb24KXFsKZih5KSA9IApcYmVnaW57Y2FzZXN9ClxmcmFjezF9e1x0aGV0YX1cY2RvdCBlXnsteS9cdGhldGF9ICYsIHkgPiAwIFxcCjAgJiwgeSBcbGUgMC4KXGVuZHtjYXNlc30KXF0KTm90ZSB0aGF0IFwoRVtZX2ldPVx0aGV0YVwpIGFuZCBcKFZbWV9pXT1cdGhldGFeMlwpLCBcKGk9MSwyLDNcKS4KCjxvbCB0eXBlPSJhIj4KPGxpPgpDb25zaWRlcjogClwoXGhhdHtcdGhldGF9XzEgPSBZXzFcKSwgClwoXGhhdHtcdGhldGF9XzIgPSBcZnJhY3sxfXsyfShZXzEgKyBZXzIpXCksClwoXGhhdHtcdGhldGF9XzMgPSBcZnJhY3sxfXszfVlfMSArIFxmcmFjezJ9ezN9WV8yXCkKXChcaGF0e1x0aGV0YX1fNCA9IFxtaW57KFlfMSwgWV8yLCBZXzMpfVwpLApcKFxoYXR7XHRoZXRhfV81ID0gXG92ZXJsaW5le1l9XCkuCgpcYmVnaW57YWxpZ24qfQpFW1xoYXR7XHRoZXRhfV8xXSAmPSBFW1lfMV0gPSBcdGhldGEuIFxcCkVbXGhhdHtcdGhldGF9XzJdICY9IFxmcmFjezF9ezJ9XGxlZnQoRVtZXzFdICsgRVtZXzJdXHJpZ2h0KSA9IFxmcmFjezF9ezJ9KFx0aGV0YSArIFx0aGV0YSkgPSBcdGhldGEuIFxcCkVbXGhhdHtcdGhldGF9XzNdICY9IFxmcmFjezF9ezN9RVtZXzFdICsgXGZyYWN7Mn17M31FW1lfMl0gPSBcZnJhY3sxfXszfVx0aGV0YSArIFxmcmFjezJ9ezN9XHRoZXRhID0gXHRoZXRhLiBcXApFW1xoYXR7XHRoZXRhfV80XSAmPSBcZnJhY3sxfXszfVx0aGV0YSBcdGV4dHsgKHNlZSBleGVyY2lzZSA2LjY1KGEpLCBwYWdlIDMyNCkufSBcXApFW1xoYXR7XHRoZXRhfV81XSAmPSBFW1xvdmVybGluZXtZfV0gPSBcdGhldGEuClxlbmR7YWxpZ24qfQpTbywgXChcaGF0e1x0aGV0YX1fMVwpLCBcKFxoYXR7XHRoZXRhfV8yXCksIFwoXGhhdHtcdGhldGF9XzNcKSwgXChcaGF0e1x0aGV0YX1fNVwpIGFyZSB1bmJpYXNlZCBlc3RpbWF0b3JzIGZvciBcKFx0aGV0YVwpLgo8L2xpPgo8bGk+CldlIGNhbGN1bGF0ZSB0aGUgdmFyaWFuY2VzIG9mIHRoZSB1bmJpYXNlZCBlc3RpbWF0b3JzLgpcYmVnaW57YWxpZ24qfQpWW1xoYXR7XHRoZXRhfV8xXSAmPSBWW1lfMV0gPSBcdGhldGFeMiBcXApWW1xoYXR7XHRoZXRhfV8yXSAmPSBcZnJhY3sxfXs0fShcdGhldGFeMiArIFx0aGV0YV4yKSA9IFxmcmFjezF9ezJ9XHRoZXRhXjIgXFwKVltcaGF0e1x0aGV0YX1fM10gJj0gXGZyYWN7MX17OX1cdGhldGFeMiArIFxmcmFjezR9ezl9XHRoZXRhXjIgPSBcZnJhY3s1fXs5fVx0aGV0YV4yIFxcClZbXGhhdHtcdGhldGF9XzVdICY9IFxmcmFjezF9ezN9XHRoZXRhXjIsClxlbmR7YWxpZ24qfQpzbyBcKFxoYXR7XHRoZXRhfV81ID0gXG92ZXJsaW5le1l9XCkgaGFzIHRoZSBsZWFzdCB2YXJpYW5jZS4KPC9saT4KPC9vbD4KCiMjUHJvYmxlbSA4LjggKDYgcG9pbnRzKQpXZSBhcmUgZ2l2ZW4gdGhhdCBcKFlfMSwgXGRvdHNjLCBZX25cKSByZXByZXNlbnRzIGEgcmFuZG9tIHNhbXBsZSBmcm9tIGEgdW5pZm9ybSBkaXN0cmlidXRpb24gb24gdGhlIGludGVydmFsIFwoKFx0aGV0YSxcdGhldGEgKyAxKVwpLiBUaGVyZWZvcmUsIFwoRVtZX2ldID0gXHRoZXRhICsgXGZyYWN7MX17Mn1cKSwgXChWW1lfaV0gPSAxLzEyXCksIFwoaT0xLFxkb3RzYyxuXCkuCjxvbCB0eXBlPSJhIj4KPGxpPgpUaGUgYmlhcyBmb3IgXChcb3ZlcmxpbmV7WX1cKSAoYXMgYW4gZXN0aW1hdG9yIGZvciBcKFx0aGV0YVwpKSBpcwpcWwpCID0gRVtcb3ZlcmxpbmV7WX1dIC0gXHRoZXRhID0gKFx0aGV0YSArIDEvMikgLSBcdGhldGEgPSAxLzIuClxdCjwvbGk+CjxsaT4KTGV0IFwoXGhhdHtcdGhldGF9ID0gXG92ZXJsaW5le1l9IC0gMS8yXCksIHRoZW4gXChFW1xoYXR7XHRoZXRhfV0gPSBFW1xvdmVybGluZXtZfSAtIDEvMl0gPSAoXHRoZXRhICsgMS8yKSAtIDEvMiA9IFx0aGV0YVwpLCBzbyBcKFxvdmVybGluZXtZfSAtIDEvMlwpIGlzIGFuIHVuYmlhc2VkIGVzdGltYXRvciBmb3IgXChcdGhldGFcKS4KPC9saT4KPGxpPgpcKFx0ZXh0e01TRX1bXG92ZXJsaW5le1l9XSA9IFZbXG92ZXJsaW5le1l9XSArIEJeMiA9IFxmcmFjezF9ezEybn0gKyBcZnJhY3sxfXs0fSA9IFxmcmFjezNuKzF9ezEybn1cKS4KPC9saT4KPC9vbD4KCiMjUHJvYmxlbSA4LjEwICg2IHBvaW50cykKTGV0IFwoWV8xLCBcZG90c2MsIFlfblwpIGRlbm90ZSBhIHJhbmRvbSBzYW1wbGUgZnJvbSBhIHBvcHVsYXRpb24gd2hvc2UgZGVuc2l0eSBhbmQgZGlzdHJpYnV0aW9uIGZ1bmN0aW9ucyBhcmUKXGJlZ2lue2FsaWduKn0KZih5KSAmPQpcYmVnaW57Y2FzZXN9ClxhbHBoYVxjZG90IHlee1xhbHBoYS0xfS9cdGhldGFeXGFscGhhICYgLCAwIFxsZSB5IFxsZSBcdGhldGEgXFwKMCAmICwgXHRleHR7IG90aGVyd2lzZX0KXGVuZHtjYXNlc30gXFwKRih5KSAmPQpcYmVnaW57Y2FzZXN9CjEgJiAsIHkgPiBcdGhldGEgXFwKKHkvXHRoZXRhKV5cYWxwaGEgJiAsIDAgXGxlIHkgXGxlIFx0aGV0YSBcXAowICYgLCB5IDwgMC4KXGVuZHtjYXNlc30KXGVuZHthbGlnbip9CkxldCBcKFxoYXR7XHRoZXRhfSA9IFxtYXhce1lfMSwgXGRvdHNjLCBZX25cfSA9IFlfeyhuKX1cKS4gVGhlbiBcKFxoYXR7XHRoZXRhfVwpIGhhcyBkZW5zaXR5ClxbCmdfeyhuKX0oeSkgPQpcYmVnaW57Y2FzZXN9Cm4gXGFscGhhXGNkb3QgeV57biBcYWxwaGEgLSAxfSAvIFx0aGV0YV57biBcYWxwaGF9ICYgLCAwIFxsZSB5IFxsZSBcdGhldGEgXFwKMCAmICwgXHRleHR7IG90aGVyd2lzZS59ClxlbmR7Y2FzZXN9ClxdCjxvbCB0eXBlPSJhIj4KPGxpPgpXZSB1c2UgXChnX3sobil9XCkgdG8gZmluZCB0aGUgZXhwZWN0ZWQgdmFsdWUgb2YgXChcaGF0e1x0aGV0YX1cKS4KXFsKRVtcaGF0e1x0aGV0YX1dID0KXGludF97MH1ee1x0aGV0YX0gblxhbHBoYVxjZG90IHlcY2RvdCB5XntuXGFscGhhIC0gMX0vXHRoZXRhXntuXGFscGhhfVwsZHkgPQpcbGVmdC5cZnJhY3tuXGFscGhhfXtuXGFscGhhICsgMX1cY2RvdCBcZnJhY3t5XntuXGFscGhhICsgMX19e1x0aGV0YV57blxhbHBoYX19XHJpZ2h0fF8wXntcdGhldGF9ID0KXGZyYWN7blxhbHBoYX17blxhbHBoYSArIDF9XGNkb3QgXHRoZXRhLgpcXQpTaW5jZSBcKEVbXGhhdHtcdGhldGF9XSBcbmUgXHRoZXRhXCksIHdlIHNlZSB0aGF0IFwoXGhhdHtcdGhldGF9XCkgaXMgYSBiaWFzZWQgZXN0aW1hdG9yIGZvciBcKFx0aGV0YVwpLgo8L2xpPgo8bGk+ClwoXGZyYWN7blxhbHBoYSArIDF9e25cYWxwaGF9XGNkb3QgXGhhdHtcdGhldGF9XCkgaXMgYW4gdW5iaWFzZWQgZXN0aW1hdG9yIGZvciBcKFx0aGV0YVwpLgo8L2xpPgo8bGk+ClRoZSBiaWFzIGZvciBcKFxoYXR7XHRoZXRhfVwpIGlzIFwoQiA9IFxmcmFje25cYWxwaGF9e25cYWxwaGErMX1cY2RvdFx0aGV0YSAtIFx0aGV0YSA9IC1cdGhldGEgLyAoblxhbHBoYSArIDEpXCkuIFdlIGNhbGN1bGF0ZSB0aGUgc2Vjb25kIG1vbWVudApcWwpFXGxlZnRbXGhhdHtcdGhldGF9XjJccmlnaHRdID0gXGludF97MH1ee1x0aGV0YX0gblxhbHBoYVxjZG90IHlee25cYWxwaGEgKyAxfSAvIFx0aGV0YV57blxhbHBoYX1cLGR5Cj0gXGZyYWN7blxhbHBoYX17blxhbHBoYSArIDJ9XGNkb3QgXHRoZXRhXjIuClxdCkZpbmFsbHksIHdlIGNhbGN1bGF0ZSB0aGUgbWVhbiBzcXVhcmUgZXJyb3IKXGJlZ2lue2FsaWduKn0KXHRleHR7TVNFfVtcaGF0e1x0aGV0YX1dICY9IFZbXGhhdHtcdGhldGF9XSArIEJeMiBcXAomPSBcZnJhY3tuXGFscGhhfXtuXGFscGhhICsgMn1cY2RvdCBcdGhldGFeMiAKLSBcbGVmdChcZnJhY3tuXGFscGhhfXtuXGFscGhhICsgMX1ccmlnaHQpXjJcY2RvdCBcdGhldGFeMgorIFxmcmFjezF9eyhuXGFscGhhICsgMSleMn1cY2RvdCBcdGhldGFeMiBcXAomPSBcZnJhY3syXGNkb3QgXHRoZXRhXjJ9eyhuXGFscGhhICsgMSkoblxhbHBoYSArIDIpfS4KXGVuZHthbGlnbip9CjwvbGk+Cjwvb2w+CgojI1Byb2JsZW0gOC4xMiAoNiBwb2ludHMpClN1cHBvc2UgXChZXzEsIFxkb3RzYywgWV9uXCkgY29uc3RpdHV0ZSBhIHJhbmRvbSBzYW1wbGUgZnJvbSBhIG5vcm1hbCBkaXN0cmlidXRpb24gd2l0aCBtZWFuIFwoXG11XCkgYW5kIHZhcmlhbmNlIFwoXHNpZ21hXjJcKS4KCkxldCBcKFkgPSAobiAtIDEpXGNkb3QgU14yIC8gXHNpZ21hXjJcKSwgd2hpY2ggd2Uga25vdyBoYXMgdGhlIApcKFxjaGleMlwpIHByb2JhYmlsaXR5IGRpc3RyaWJ1dGlvbiB3aXRoIFwoXG51ID0gbiAtIDFcKSBkZWdyZWVzIG9mIGZyZWVkb20uIEJ5IGV4ZXJjaXNlIDQuOTAoYyksIFwoRVtcc3FydHtZfV0gPSBcc3FydHsyfVxjZG90IFxHYW1tYShuLzIpIC8gXEdhbW1hKChuLTEpLzIpXCkuIAo8b2wgdHlwZT0iYSI+CjxsaT4KTm93LCBcKFMgPSB7XHNpZ21hfVxjZG90IFxzcXJ0e1l9IC8ge1xzcXJ0e24tMX19XCkgYW5kIHRoZXJlZm9yZQpcWwpFW1NdID0gXGZyYWN7XHNpZ21hfXtcc3FydHtuLTF9fVxjZG90IEVcbGVmdFtcc3FydHtZfVxyaWdodF0gID0gClxzcXJ0e1xmcmFjezJ9e24tMX19XGNkb3QgXGZyYWN7XEdhbW1hKG4vMil9e1xHYW1tYVxsZWZ0KChuLTEpLzIpXHJpZ2h0KX1cY2RvdCBcc2lnbWEuClxdCldlIHNlZSB0aGF0IFwoU1wpIGlzIGEgYmlhc2VkIGVzdGltYXRvciBmb3IgXChcc2lnbWFcKSwgZXZlbiB0aG91Z2ggXChTXjJcKSBpcyBhbiB1bmJpYXNlZCBlc3RpbWF0b3IgZm9yIFwoXHNpZ21hXjJcKS4KPC9saT4KPGxpPgpXZSBjYW4gYWRqdXN0IFwoU1wpIGJ5IGEgbXVsdGlwbGljYXRpdmUgY29uc3RhbnQgdG8gY3JlYXRlIHRoZSB1bmJpYXNlZCBlc3RpbWF0b3IgZm9yIFwoXHNpZ21hXCkKXFsKXGhhdHtcc2lnbWF9ID0gXHNxcnR7XGZyYWN7bi0xfXsyfX1cY2RvdCBcZnJhY3tcR2FtbWFcbGVmdCgobi0xKS8yKVxyaWdodCl9e1xHYW1tYShuLzIpfVxjZG90IFMuClxdCjwvbGk+CjxsaT4KQW4gdW5iaWFzZWQgZXN0aW1hdG9yIGZvciBcKFxtdSAtIHpfe1xhbHBoYX1cY2RvdCBcc2lnbWFcKSBpcyBcKFxvdmVybGluZXtZfSAtIHpfe1xhbHBoYX1cY2RvdCBcaGF0e1xzaWdtYX1cKSwgd2hlcmUgXChcb3ZlcmxpbmV7WX0gPSBcZnJhY3sxfXtufVxzdW1fe2k9MX1ee259IFlfe2l9XCkgYW5kIFwoXGhhdHtcc2lnbWF9XCkgaXMgdGhlIGVzdGltYXRvciBkZWZpbmVkIGluIHBhcnQgKGIpLgo8L2xpPgo8L29sPgoKIyNQcm9ibGVtIDguMTQgKDIgcG9pbnRzKQpMZXQgXChZXzEsIFxkb3RzYywgWV9uXCkgZGVub3RlIGEgcmFuZG9tIHNhbXBsZSBmcm9tIGEgcG9wdWxhdGlvbiB3aXRoIGEgdW5pZm9ybSBkaXN0cmlidXRpb24gb24gdGhlIGludGVydmFsClwoKDAsXHRoZXRhKVwpLiBMZXQgXChZX3soMSl9ID0gXG1pblx7WV8xLCBcZG90c2MsIFlfblx9XCkuIFJlY2FsbCB0aGF0IHRoZSBkZW5zaXR5IGFuZCBkaXN0cmlidXRpb24gZnVuY3Rpb25zIGZvciB0aGlzIHVuaWZvcm0gZGlzdHJpYnV0aW9uIGFyZSwgcmVzcGVjdGl2ZWx5LApcYmVnaW57YWxpZ24qfQpmKHkpICY9IApcYmVnaW57Y2FzZXN9ClxmcmFjezF9e1x0aGV0YX0gJiAsIDAgPCB5IDwgXHRoZXRhIFxcCjAgJiAsIFx0ZXh0eyBvdGhlcndpc2V9LApcZW5ke2Nhc2VzfSBcXApGKHkpICY9ClxiZWdpbntjYXNlc30KMSAmICwgeSBcZ2UgXHRoZXRhIFxcClxmcmFje3l9e1x0aGV0YX0gJiAsIDAgPCB5IDwgXHRoZXRhIFxcCjAgJiAsIHkgXGxlIDAuClxlbmR7Y2FzZXN9ClxlbmR7YWxpZ24qfQpUaHVzLCB0aGUgZGVuc2l0eSBmdW5jdGlvbiBmb3IgXChZX3soMSl9XCkgaXMKXFsKZ197KDEpfSh5KSA9Cm5cY2RvdCBcbGVmdFsgMSAtIEYoeSlccmlnaHRdXntuLTF9XGNkb3QgZih5KSA9ClxiZWdpbntjYXNlc30KXGZyYWN7bn17XHRoZXRhfSBcY2RvdCBcbGVmdCggMSAtIFxmcmFje3l9e1x0aGV0YX1ccmlnaHQpXntuLTF9ICYgLCAwIDwgeSA8IFx0aGV0YSBcXAowICYgLCBcdGV4dHsgb3RoZXJ3aXNlfQpcZW5ke2Nhc2VzfS4KXF0KTm93IHRoYXQgd2UgaGF2ZSB0aGUgZGVuc2l0eSBmb3IgXChZX3soMSl9XCksIHdlIGNhbiBjYWxjdWxhdGUKXFsKRVtZX3soMSl9XSA9IFxpbnRfezB9XntcdGhldGF9IFxmcmFje255fXtcdGhldGF9XGNkb3QgXGxlZnQoIDEgLSBcZnJhY3t5fXtcdGhldGF9XHJpZ2h0KV57bi0xfVwsZHkgPSBcZnJhY3tcdGhldGF9eyhuICsgMSl9LgpcXQpXZSBzZWUgdGhhdCBcKFlfeygxKX1cKSBpcyBhIGJpYXNlZCBlc3RpbWF0b3IgZm9yIFwoXHRoZXRhXCksIGJ1dCB3ZSBtYXkgYWRqdXN0IGl0IGJ5IGEgbXVsdGlwbGljYXRpdmUgY29uc3RhbnQgdG8gb2J0YWluIHRoZSB1bmJpYXNlZCBlc3RpbWF0b3IgXChcaGF0e1x0aGV0YX0gPSAobiArIDEpXGNkb3QgWV97KDEpfVwpLgoKIyNQcm9ibGVtIDguMTYgKDQgcG9pbnRzKQpTdXBwb3NlIHRoYXQgXChZXzEsIFlfMiwgWV8zLCBZXzRcKSBkZW5vdGVzIGEgcmFuZG9tIHNhbXBsZSBmcm9tIGEgcG9wdWxhdGlvbiB3aXRoIGFuIGV4cG9uZW50aWFsIGRpc3RyaWJ1dGlvbiB3aG9zZSBkZW5zaXR5IGlzIGdpdmVuIGJ5ClxbCmYoeSkgPQpcYmVnaW57Y2FzZXN9ClxmcmFjezF9e1x0aGV0YX1cY2RvdCBlXnsteS9cdGhldGF9ICYgLCB5ID4gMCBcXAowICYgLCBcdGV4dHsgb3RoZXJ3aXNlfQpcZW5ke2Nhc2VzfS4KXF0KUmVjYWxsIHRoYXQgXChZXCkgaXMgYSBHYW1tYSBkaXN0cmlidXRpb24gd2l0aCBwYXJhbWV0ZXJzIFwoXGFscGhhID0gMVwpIGFuZCBcKFxiZXRhID0gXHRoZXRhXCkuCkJ5IGV4ZXJjaXNlIDQuODkoYSksIFwoRVtZXmFdID0gXHRoZXRhXnthfVxjZG90IFxHYW1tYSgxICsgYSkvXEdhbW1hKDEpID0gXHRoZXRhXnthfVxjZG90IGFcY2RvdCBcR2FtbWEoYSlcKSwgCnRoZXJlZm9yZQpcWwpFWyBcc3FydHtZfSBdID0gXGZyYWN7MX17Mn0gXEdhbW1hXGxlZnQoXGZyYWN7MX17Mn1ccmlnaHQpXGNkb3QgXHNxcnR7XHRoZXRhfSA9IFxmcmFjezF9ezJ9XHNxcnR7XHBpXGNkb3QgXHRoZXRhfS4KXF0KPG9sIHR5cGU9ImEiPgo8bGk+CkxldCBcKFggPSBcc3FydHtZXzFcY2RvdCBZXzJ9XCkuIFNpbmNlIFwoWV8xXCkgYW5kIFwoWV8yXCkgYXJlIGluZGVwZW5kZW50LApcWwpFW1hdID0gRVtZXzFdXGNkb3QgRVtZXzJdID0gXGZyYWN7MX17NH1ccGlcY2RvdCBcdGhldGEsClxdCnNvIGFuIHVuYmlhc2VkIGVzdGltYXRvciBmb3IgXChcdGhldGFcKSBpcyBcKFxmcmFjezR9e1xwaX1cc3FydHtZXzFcY2RvdCBZXzJ9XCkuCjwvbGk+CjxsaT4KTGV0IFwoVyA9IFxzcXJ0e1lfMVxjZG90IFlfMlxjZG90IFlfM1xjZG90IFlfNH1cKSwgdGhlbgpcWwpFW1ddID0gRVtZXzFdXGNkb3QgRVtZXzJdXGNkb3QgRVtZXzNdXGNkb3QgRVtZXzRdID0gXGZyYWN7MX17MTZ9XHBpXjJcY2RvdCBcdGhldGFeMi4KXF0KVGh1cywgYW4gdW5iaWFzZWQgZXN0aW1hdG9yIGZvciBcKFx0aGV0YV4yXCkgaXMgClwoXGZyYWN7MTZ9e1xwaV4yfVxzcXJ0e1lfMVxjZG90IFlfMlxjZG90IFlfM1xjZG90IFlfNH1cKS4KPC9saT4KPC9vbD4KCgoKCgoKCgoKCg==