Skip to content

সমাধান — অধ্যায় ৭.৬ · Independence, 0–1 Law ও SLLN

অধ্যায় ফাইল: part-7-measure-theoretic/07-06-independence-zero-one-slln.md (§৭ অনুশীলনী)। গোটা অংশে \((\Omega,\mathcal F,\mathbb P)\) একটি probability space; random variable = measurable \(X:\Omega\to\mathbb R\) (7.3), \(\mathbb E[X]=\int_\Omega X\,d\mathbb P\) (7.4), \(X\in L^1\iff\mathbb E\lvert X\rvert<\infty\)। "i.o." = infinitely often; \(\limsup_n A_n=\bigcap_N\bigcup_{n\ge N}A_n=\{A_n\ \text{i.o.}\}\); tail σ-algebra \(\mathcal T=\bigcap_m\sigma(X_m,X_{m+1},\dots)\)

canonical তথ্য (সংখ্যাগত উত্তর reproducible, seed np.random.default_rng(20260619)): - SLLN \(\text{Exp}(1)\) (\(\mathbb E[X]=1\)): \(n=10,10^2,10^3,10^4,10^5,10^6\to\bar X_n=0.5194,\,0.9710,\,0.9603,\,0.9841,\,0.9978,\,1.0007\) (→ \(1\) a.s.); \(\text{Bernoulli}(0.3)\): \(\bar X_{10^6}=0.2999\to0.3\)। - Cauchy (\(\mathbb E\lvert X\rvert=\infty\)): \(\bar X_n\) থামে না, ঘোরে — \(1.126,\,0.851,\,-0.173\) at \(n=10^2,10^4,10^6\) (চিহ্নও বদলায়)। - BC-I \(\mathbb P(A_n)=1/n^2\): \(\sum=\pi^2/6\approx1.6449<\infty\), ≈২ ঘটনা ঘটে (saturate), \(\mathbb P(\text{i.o.})=0\)। - BC-II \(\mathbb P(A_n)=1/n\): \(\sum=\infty\), ≈৮ ঘটনা (\(N=10^5\)), \(\sim\ln N\approx11.5\) হারে বাড়ে, \(\mathbb P(\text{i.o.})=1\)। - সীমারেখা \(\mathbb P(A_n)=1/(n\log n)\): \(\sum\sim\log\log N=\infty\) (অতি-ধীর অপসারী), \(\mathbb P(\text{i.o.})=1\)


ক · ধারণাগত (conceptual)

সমাধান ১ (★)

(ক) pairwise independence (স্বাধীনতা)। \(X_1,X_2\) স্বাধীন fair coin, তাই \((X_1,X_2)\) চারটি মান \((0,0),(0,1),(1,0),(1,1)\) প্রতিটি সম্ভাবনা \(\tfrac14\)-এ নেয়। \(X_3=X_1\oplus X_2\) (XOR): \(X_3=0\) যদি \(X_1=X_2\), নয়তো \(1\)। তিনটি জোড়া যাচাই: - \(X_1,X_2\): সংজ্ঞা থেকেই স্বাধীন। - \(X_1,X_3\): \(\mathbb P(X_1{=}a,X_3{=}c)\) — স্থির \(X_1=a\)-তে \(X_3=c\) মানে \(X_2=a\oplus c\), যার সম্ভাবনা \(\tfrac12\); তাই \(\mathbb P(X_1{=}a,X_3{=}c)=\tfrac12\cdot\tfrac12=\tfrac14=\mathbb P(X_1{=}a)\mathbb P(X_3{=}c)\) (কারণ \(X_3\)-ও marginally fair: \(\mathbb P(X_3{=}c)=\tfrac12\))। স্বাধীন। - \(X_2,X_3\): symmetry-তে একইভাবে স্বাধীন।

তাই \(X_1,X_2,X_3\) জোড়ায়-জোড়ায় স্বাধীন

(খ) mutual স্বাধীনতা ব্যর্থ। \(X_3\) প্রথম দুটোর deterministic ফাংশন (\(X_3=X_1\oplus X_2\)), তাই \(X_1,X_2\) জানলে \(X_3\) সম্পূর্ণ নির্ধারিত — কোনো বাড়তি দৈবতা নেই। ফলে তিন-মুখী factorization ভাঙে; কংক্রিটভাবে: $$ \mathbb P(X_1{=}0,X_2{=}0,X_3{=}1)=0\qquad(\text{কারণ }0\oplus0=0\ne1), $$ অথচ গুণফল \(\mathbb P(X_1{=}0)\mathbb P(X_2{=}0)\mathbb P(X_3{=}1)=\tfrac12\cdot\tfrac12\cdot\tfrac12=\tfrac18\ne0\)। একটিমাত্র নির্বাচন ভাঙলেই পারস্পরিক স্বাধীনতা যায়।

(গ) এক বাক্যে। pairwise (সব জোড়া স্বাধীন) থেকে mutual (সব উপসেটে factorize) আসে না — ঠিক এই ফাঁক বন্ধ করতেই σ-algebra-র স্বাধীনতা-সংজ্ঞা "যেকোনো \(A_i\in\mathcal F_i\) নির্বাচনে \(\mathbb P(\bigcap A_i)=\prod\mathbb P(A_i)\)" দাবি করে (যা \(\Omega\) বসিয়ে সব উপসেট ঢাকে), শুধু জোড়া নয়। (← 2.2)

সমাধান ২ (★)

(ক) tail event কী। একটি ঘটনা \(A\) tail event যদি \(A\in\mathcal T=\bigcap_m\sigma(X_m,X_{m+1},\dots)\) — সমতুল্যভাবে, প্রতিটি \(m\)-এ \(A\) কেবল \(X_m,X_{m+1},\dots\)-এর উপর নির্ভর করে, অর্থাৎ প্রথম \(m-1\)টি \(X_i\) যা-খুশি বদলালেও \(A\)-র সত্য-মিথ্যা বদলায় না — সব \(m\)-এ। সংক্ষেপে: tail event "অসীম-দূরের লেজ"-এর তথ্য, কোনো সসীম উপসর্গে নির্ভরশীল নয়।

(খ) কোনগুলো tail। - (i) \(\{\sum_n X_n\ \text{converges}\}\)tail। প্রথম \(m-1\)টি পদ বদলালে কেবল একটি সসীম যোগফল সরে, যা অভিসারিতা বদলায় না। - (ii) \(\{\limsup_n\bar X_n>c\}\)tail। \(\limsup\) ও চলমান-গড়ের asymptotic আচরণ লেজ-নির্ভর; সসীম-সংখ্যক পদ বদলালে \(\bar X_n\)-এর limit-আচরণ অপরিবর্তিত। - (iii) \(\{X_1>0\}\)tail নয়। শুধু \(X_1\)-এর উপর নির্ভর; \(X_1\) বদলালেই ঘটনা বদলায়, তাই \(\sigma(X_2,X_3,\dots)\)-এ নেই। - (iv) \(\{\sum_{n=1}^{100}X_n>5\}\)tail নয়। কেবল \(X_1,\dots,X_{100}\)-এর (সসীম) ফাংশন; \(\sigma(X_{101},X_{102},\dots)\)-এ পড়ে না।

(গ) কেন \(0.5\) নয়। স্বাধীন \((X_n)\)-এ Kolmogorov 0–1 law বলে \(\mathcal T\) নিজের থেকে স্বাধীন (\(\mathcal T\perp\mathcal T\)), তাই যেকোনো \(A\in\mathcal T\)-এ \(\mathbb P(A)=\mathbb P(A\cap A)=\mathbb P(A)^2\), যার একমাত্র সমাধান \(\mathbb P(A)\in\{0,1\}\)। কাজেই \(0.5\) অসম্ভব — একটি সত্যিকারের asymptotic (tail) প্রশ্নের উত্তর নির্ধারিত (deterministic): হয় প্রায়-সব নমুনা-পথে ঘটে, নয় প্রায়-কোনোটিতেই নয়।

সমাধান ৩ (★★)

(ক) BC-I-এ স্বাধীনতা লাগে না। চাবি-ধর্ম countable subadditivity (7.2)। \(B_N=\bigcup_{n\ge N}A_n\) লিখলে $$ \mathbb P(B_N)=\mathbb P\Big(\bigcup_{n\ge N}A_n\Big)\le\sum_{n\ge N}\mathbb P(A_n), $$ আর \(\sum_n\mathbb P(A_n)<\infty\) বলে ডান পাশ (\(=\) অভিসৃত ধারার লেজ) \(N\to\infty\)-এ \(0\)-তে নামে। তাই \(\mathbb P(A_n\ \text{i.o.})=\lim_N\mathbb P(B_N)=0\) — কেবল additivity ও convergence, কোনো স্বাধীনতা নয়। স্বজ্ঞা: সম্ভাবনার যোগফল সসীম মানে গড়ে সসীম-সংখ্যক ঘটনা (\(\mathbb E[\sum_n\mathbf 1_{A_n}]=\sum_n\mathbb P(A_n)<\infty\)), তাই "অসীম-বার ঘটা" শূন্য-সম্ভাব্য।

(খ) BC-II-তে স্বাধীনতা কোথায়। এখানে দেখাতে হয় "\(N\)-এর পর কিছুই ঘটে না" (\(C_N=\bigcap_{n\ge N}A_n^c\)) শূন্য-সম্ভাব্য। সসীম intersection-এ স্বাধীনতা দিয়ে গুণফলে ভাঙা হয়: $$ \mathbb P\Big(\bigcap_{n=N}^{M}A_n^c\Big)=\prod_{n=N}^{M}\big(1-\mathbb P(A_n)\big)\le\exp\Big(-\sum_{n=N}^{M}\mathbb P(A_n)\Big)\xrightarrow[M\to\infty]{}e^{-\infty}=0, $$ যেখানে \(1-x\le e^{-x}\) অসমতা ও অপসারী লেজ \(\sum_{n\ge N}\mathbb P(A_n)=\infty\) ব্যবহৃত। স্বাধীনতা ছাড়া এই factorization-টাই সম্ভব নয় — এটিই BC-I ও BC-II-এর গঠনগত পার্থক্য।

(গ) প্রতিউদাহরণ (স্বাধীনতা ছাড়া BC-II মিথ্যা)। একই ঘটনা পুনরাবৃত্তি করুন: \(A_n=A\) সব \(n\)-এ, যেখানে \(0<\mathbb P(A)<1\) (যেমন \(\mathbb P(A)=0.1\))। তখন $$ \sum_n\mathbb P(A_n)=\sum_n\mathbb P(A)=\infty\quad(\text{অপসারী}), $$ কিন্তু \(\{A_n\ \text{i.o.}\}=\bigcap_N\bigcup_{n\ge N}A=A\) (প্রতিটি union ঠিক \(A\)), তাই \(\mathbb P(A_n\ \text{i.o.})=\mathbb P(A)=0.1\ne1\)। অপসারী যোগফল সত্ত্বেও i.o.-সম্ভাবনা \(1\) নয় — কারণ ঘটনাগুলো পরস্পর-নির্ভরশীল (আসলে অভিন্ন)। এটি দেখায় BC-II-তে স্বাধীনতা সত্যিই অপরিহার্য, অলংকার নয়।


খ · গণনামূলক (computational)

সমাধান ৪ (★)

প্রতিটি ক্ষেত্রে যোগফলের আচরণ যাচাই করে Borel–Cantelli প্রয়োগ:

(ক) \(\mathbb P(A_n)=1/n^2\) (i) \(\sum_{n\ge1}1/n^2=\pi^2/6\approx1.6449<\infty\) (Basel series — অভিসারী)। (ii) যোগফল সসীম, তাই BC-I (স্বাধীনতা লাগে না)। (iii) \(\mathbb P(A_n\ \text{i.o.})=0\) — a.s. কেবল সসীম-সংখ্যক \(A_n\) ঘটে।

(খ) \(\mathbb P(A_n)=1/n\) (i) \(\sum_{n\ge1}1/n=\infty\) (harmonic series — অপসারী, বৃদ্ধি \(\sum_{n\le N}1/n\approx\ln N\))। (ii) অপসারী + independent (স্বাধীন), তাই BC-II। (iii) \(\mathbb P(A_n\ \text{i.o.})=1\) — a.s. অসীম-সংখ্যক \(A_n\) ঘটে।

(গ) \(\mathbb P(A_n)=1/(n\log n)\) (\(n\ge2\))। (i) integral test (সমাধান ৫ক): \(\int_2^\infty\frac{dx}{x\log x}=[\log\log x]_2^\infty=\infty\), তাই অপসারী (অতি-ধীর, আংশিক-যোগফল \(\sim\log\log N\))। (ii) অপসারী + স্বাধীন, তাই BC-II। (iii) \(\mathbb P(A_n\ \text{i.o.})=1\)

সারসংক্ষেপ-টেবিল:

ক্ষেত্র \(\mathbb P(A_n)\) \(\sum\mathbb P(A_n)\) lemma \(\mathbb P(A_n\ \text{i.o.})\)
(ক) \(1/n^2\) \(\pi^2/6\approx1.6449\) (অভিসারী) BC-I \(0\)
(খ) \(1/n\) \(\infty\) (\(\sim\ln N\)) BC-II \(1\)
(গ) \(1/(n\log n)\) \(\infty\) (\(\sim\log\log N\)) BC-II \(1\)

লক্ষণীয়: \(1/n^2\)\(1/n\)-এর মাঝে সীমারেখা \(1/(n\log n)\) এখনও অপসারী দিকেই — তাই i.o. \(=1\), কিন্তু একদম গা ঘেঁষে।

সমাধান ৫ (★★)

(ক) integral test দিয়ে অপসারিতা। \(f(x)=\frac1{x\log x}\) (\(x\ge2\)) ধনাত্মক, হ্রাসমান, তাই integral test প্রযোজ্য। \(u=\log x\), \(du=\frac{dx}{x}\) substitution-এ: $$ \int_2^\infty\frac{dx}{x\log x}=\int_{\log 2}^\infty\frac{du}{u}=[\log u]_{\log 2}^\infty=\infty . $$ integral অপসারী, তাই \(\sum_{n\ge2}\frac1{n\log n}\)-ও অপসারী। (তুলনায় \(\frac1{n(\log n)^2}\) অভিসারী হত, কারণ \(\int\frac{du}{u^2}<\infty\) — exponent \(\log\)-এ ঠিক \(1\)-এর সীমারেখা।)

(খ) বৃদ্ধি-হার। integral test-এর তুলনা থেকে আংশিক-যোগফল \(\sum_{2\le n\le N}\frac1{n\log n}\approx\log\log N\) (অগ্রণী পদ)। আনুমানিক মান: - \(N=10^3\): \(\log\log 10^3=\log(3\ln 10)=\log(6.908)\approx1.93\)। - \(N=10^6\): \(\log\log 10^6=\log(6\ln 10)=\log(13.82)\approx2.63\)

তুলনায় \(\sum 1/n\approx\ln N\) দেয় \(\ln 10^3\approx6.91\), \(\ln 10^6\approx13.82\)। অর্থাৎ \(N\) হাজার-গুণ বাড়লেও \(\log\log N\) মোটে \(1.93\to2.63\) (≈\(0.7\) বৃদ্ধি), যেখানে \(\ln N\) \(6.91\to13.82\) (দ্বিগুণ) — \(\log\log\) অসম্ভব ধীরে বাড়ে, তবু সীমাহীন।

(গ) উপসংহার। অপসারী যোগফল + স্বাধীন ঘটনা ⇒ BC-II\(\mathbb P(A_n\ \text{i.o.})=1\)। স্পেকট্রামে: \(1/n^2\) (যোগফল সসীম → i.o. \(0\), "থামে") আর \(1/n\) (যোগফল অসীম → i.o. \(1\), "চলে")-এর মধ্যে \(1/(n\log n)\) পড়ে "চলে" দিকে, যদিও সীমারেখার একদম গা ঘেঁষে — অপসারিতা যত ধীরই হোক, BC-II দ্বিধাহীনভাবে \(1\) দেয়।

সমাধান ৬ (★★)

fair-coin iid \(X_i\in\{0,1\}\), \(\mathbb P(X_i{=}1)=\tfrac12\), \(\mathbb E[X_i]=\tfrac12\), \(X_i\) আবদ্ধ তাই \(\mathbb E\lvert X_i\rvert=\tfrac12<\infty\) — সর্বত্র SLLN প্রযোজ্য।

(ক) \(\bar X_n=\frac1n\sum X_i\) \(X_i\) iid, \(\mathbb E[X_i]=\tfrac12\), তাই সরাসরি SLLN: $$ \bar X_n\xrightarrow{\text{a.s.}}\tfrac12 . $$ (এটাই "আপেক্ষিক কম্পাঙ্ক → সম্ভাবনা"-র সবল রূপ: মাথা-আসার হার \(\to\tfrac12\)।)

(খ) \(\frac1n\sum_{i=1}^n X_iX_{i+1}\) \(Y_i:=X_iX_{i+1}\in\{0,1\}\), \(\mathbb E[Y_i]=\mathbb P(X_i{=}1,X_{i+1}{=}1)=\tfrac12\cdot\tfrac12=\tfrac14\) (\(X_i,X_{i+1}\) স্বাধীন)। সীমা \(\tfrac14\) a.s.। সাবধানতা: \(Y_i\)-রা স্বাধীন নয় (\(Y_i,Y_{i+1}\) overlap-করা \(X_{i+1}\) ভাগ করে), তাই মূল iid-SLLN সরাসরি লাগে না — কিন্তু \((Y_i)\) stationary ও ergodic, তাই ergodic theorem (SLLN-এর সাধারণ রূপ) দেয় \(\to\mathbb E[Y_1]=\tfrac14\)। প্রাথমিক যুক্তি: জোড়/বিজোড় সূচকে ভাঙুন — \(Y_1,Y_3,Y_5,\dots\) (অ-overlap) iid, তাই তাদের গড় \(\to\tfrac14\) a.s. (SLLN); একইভাবে \(Y_2,Y_4,\dots\); দুই উপ-গড়ের combination-ও \(\to\tfrac14\)। সীমা: $$ \frac1n\sum_{i=1}^n X_iX_{i+1}\xrightarrow{\text{a.s.}}\tfrac14 . $$

(গ) \(\frac1n\sum_{i=1}^n(2X_i-1)^2\) \(X_i\in\{0,1\}\) বলে \(2X_i-1\in\{-1,+1\}\), তাই \((2X_i-1)^2=1\) সর্বদা — একটি ধ্রুবক random variable। গড় ঠিক \(1\) প্রতিটি \(n\)-এ, তাই সীমা: $$ \frac1n\sum_{i=1}^n(2X_i-1)^2=1\xrightarrow{\text{a.s.}}1 . $$ (\(Z_i=1\) ধ্রুবক, \(\mathbb E[Z_i]=1<\infty\) — SLLN তুচ্ছভাবে মেটে।)


গ · প্রমাণভিত্তিক (proof-based)

সমাধান ৭ (★★)

Borel–Cantelli lemma I — পূর্ণ প্রমাণ।

(ক) monotone tail ও continuity-from-above। সংজ্ঞা থেকে \(\limsup_n A_n=\bigcap_{N\ge1}\big(\bigcup_{n\ge N}A_n\big)\)। লিখি \(B_N:=\bigcup_{n\ge N}A_n\)। বড় \(N\)-এ union-এ কম পদ, তাই $$ B_1\supseteq B_2\supseteq B_3\supseteq\cdots,\qquad \bigcap_{N}B_N=\limsup_n A_n, $$ অর্থাৎ \(B_N\downarrow\limsup_n A_n\)। যেহেতু \(\mathbb P(B_1)\le1<\infty\) (probability measure সসীম), 7.2-এর continuity from above (উপর-থেকে-ধারাবাহিকতা) বৈধ: $$ \mathbb P\big(\limsup_n A_n\big)=\mathbb P\Big(\bigcap_N B_N\Big)=\lim_{N\to\infty}\mathbb P(B_N). \tag{1} $$

(খ) subadditivity দিয়ে বাঁধা। যেকোনো স্থির \(N\)-এ, 7.2-এর countable subadditivity সরাসরি দেয়: $$ \mathbb P(B_N)=\mathbb P\Big(\bigcup_{n\ge N}A_n\Big)\ \le\ \sum_{n\ge N}\mathbb P(A_n). \tag{2} $$

(গ) লেজ \(\to0\) ও উপসংহার। ধরা শর্ত \(\sum_{n\ge1}\mathbb P(A_n)<\infty\) — অভিসৃত অঋণাত্মক ধারা। যেকোনো অভিসৃত ধারার লেজ-যোগফল শূন্যে নামে: $$ \sum_{n\ge N}\mathbb P(A_n)=\sum_{n\ge1}\mathbb P(A_n)-\sum_{n=1}^{N-1}\mathbb P(A_n)\xrightarrow[N\to\infty]{}0, $$ কারণ আংশিক-যোগফল পূর্ণ-যোগফলে অভিসৃত। (2)-এ \(0\le\mathbb P(B_N)\le\sum_{n\ge N}\mathbb P(A_n)\to0\), তাই sandwich-এ \(\mathbb P(B_N)\to0\)। (1)-এ বসিয়ে \(\mathbb P(\limsup_n A_n)=0\)\(\blacksquare\)

কোথাও independence (স্বাধীনতা) ব্যবহৃত হয়নি — lemma সম্পূর্ণ সাধারণ, কেবল additivity ও convergence-এর উপর দাঁড়ানো।

সমাধান ৮ (★★★)

\(\{\sum_n X_n\ \text{converges}\}\) একটি tail event।

(ক) অভিসারিতা প্রথম \(N-1\) পদে নিরপেক্ষ। যেকোনো \(N\ge1\)-এর জন্য দুই ধারার পার্থক্য একটি সসীম যোগফল: $$ \sum_{n=1}^\infty X_n=\underbrace{\sum_{n=1}^{N-1}X_n}{\text{সসীম, ধ্রুবক}}+\sum^\infty X_n . $$ একটি সসীম সংখ্যা যোগ করলে অভিসারিতা বদলায় না (অভিসৃত \(+\) সসীম \(=\) অভিসৃত, এবং অপসারী \(+\) সসীম \(=\) অপসারী)। তাই $$ \sum_{n=1}^\infty X_n\ \text{অভিসৃত}\quad\Longleftrightarrow\quad\sum_{n=N}^\infty X_n\ \text{অভিসৃত}. \tag{\(\ast\)} $$

(খ) প্রতিটি \(\sigma(X_N,\dots)\)-এ, তাই \(\mathcal T\)-তে। আংশিক-যোগফল \(S_k^{(N)}:=\sum_{n=N}^k X_n\) (\(k\ge N\)) প্রতিটি \(X_N,\dots,X_k\)-এর measurable ফাংশন, তাই \(\sigma(X_N,X_{N+1},\dots)\)-measurable। অভিসারিতা মানে \(\lim_k S_k^{(N)}\)-এর অস্তিত্ব, যা measurable function-অনুক্রমের limit-অস্তিত্ব — একটি measurable ঘটনা \(\big(\{(S_k^{(N)})\ \text{Cauchy}\}=\bigcap_\varepsilon\bigcup_M\bigcap_{k,l\ge M}\{\lvert S_k^{(N)}-S_l^{(N)}\rvert<\varepsilon\}\big)\) — সম্পূর্ণরূপে \(\sigma(X_N,X_{N+1},\dots)\)-এ। \((\ast)\) বলছে ঘটনা \(\{\sum_n X_n\ \text{converges}\}\) ঠিক এটির সমান, তাই $$ {\textstyle\sum_n X_n\ \text{converges}}\in\sigma(X_N,X_{N+1},\dots)\qquad\forall N . $$ সব \(N\)-এ সত্য বলে ছেদ: $$ {\textstyle\sum_n X_n\ \text{converges}}\in\bigcap_{N\ge1}\sigma(X_N,X_{N+1},\dots)=\mathcal T . $$ অর্থাৎ এটি একটি tail event

(গ) এক বাক্যে। \((X_n)\) স্বাধীন বলে Kolmogorov 0–1 law প্রযোজ্য, তাই \(\mathbb P(\sum_n X_n\ \text{converges})\in\{0,1\}\) — অভিসারিতা একটি নির্ধারিত ব্যাপার (হয় a.s. অভিসৃত, নয় a.s. নয়); একইভাবে \(\{\bar X_n\ \text{অভিসৃত}\}\) tail, তাই সীমা থাকলে তা a.s. ধ্রুবক, আর SLLN শুধু সেই ধ্রুবকটি \(\mu\) বলে চিহ্নিত করে।

সমাধান ৯ (★★)

0–1 law-এর হৃৎপিণ্ড।

(ক) \(\mathbb P(A)=\mathbb P(A)^2\Rightarrow\mathbb P(A)\in\{0,1\}\) \(A\) নিজের থেকে স্বাধীন মানে \(\mathbb P(A\cap A)=\mathbb P(A)\mathbb P(A)\); কিন্তু \(A\cap A=A\), তাই \(\mathbb P(A)=\mathbb P(A)^2\)। ধরি \(p=\mathbb P(A)\): $$ p=p^2\ \Longrightarrow\ p-p^2=0\ \Longrightarrow\ p(1-p)=0\ \Longrightarrow\ p=0\ \text{বা}\ p=1 . $$ (একটি বাস্তব সংখ্যা নিজের বর্গের সমান হয় কেবল \(0\)\(1\)-এ।)

(খ) কেন \(\mathcal T\perp\mathcal T\) (যুক্তি-শৃঙ্খল)। তিন ধাপ: 1. প্রতিটি \(n\)-এ \(\mathcal T\subseteq\sigma(X_{n+1},X_{n+2},\dots)\) (ছেদের একটি সদস্য)। আর \((X_k)\) স্বাধীন বলে \(\sigma(X_1,\dots,X_n)\perp\sigma(X_{n+1},\dots)\), তাই \(\mathcal T\perp\sigma(X_1,\dots,X_n)\)প্রতিটি \(n\)-এ। 2. \(\mathcal A:=\bigcup_n\sigma(X_1,\dots,X_n)\) একটি π-system (nested, তাই intersection-বদ্ধ) এবং \(\sigma(\mathcal A)=\sigma(X_1,X_2,\dots)\)। ধাপ-১ বলছে \(\mathcal T\) এই π-system \(\mathcal A\) থেকে স্বাধীন; π-system criterion (প্রমাণ ১) দিয়ে \(\mathcal T\perp\sigma(X_1,X_2,\dots)\)। 3. কিন্তু \(\mathcal T\subseteq\sigma(X_1,X_2,\dots)\) (প্রতিটি \(\sigma(X_m,\dots)\) এর উপসেট, তাই ছেদও)। একটি σ-algebra তার ধারক থেকে স্বাধীন থাকলে বিশেষত নিজের থেকেও: \(\mathcal T\perp\mathcal T\)

(গ) tail RV a.s. ধ্রুবক। ধরা যাক \(Y\) একটি \(\mathcal T\)-measurable random variable (যেমন \(\limsup_n X_n\))। প্রতিটি \(t\)-এ \(\{Y\le t\}\in\mathcal T\), তাই (ক) দিয়ে \(F_Y(t)=\mathbb P(Y\le t)\in\{0,1\}\)। অর্থাৎ \(F_Y\) একটি অ-হ্রাসমান, ডান-অবিচ্ছিন্ন, \(\{0,1\}\)-মানের function — যা অবশ্যই একটি বিন্দু \(c=\inf\{t:F_Y(t)=1\}\)-তে \(0\) থেকে \(1\)-এ লাফ দেয়। তাই \(\mathbb P(Y=c)=F_Y(c)-F_Y(c^-)=1-0=1\)\(Y\) a.s. ধ্রুবক (\(=c\))। \(\blacksquare\)


ঘ · কোডিং (coding)

সমাধান ১০ (★★)

SLLN অভিসরণ (\(\text{Exp}(1)\to1\)) বনাম Cauchy non-convergence — পাশাপাশি।

import numpy as np

# --- (ক) Exp(1): E|X| = 1 < ∞, তাই SLLN প্রযোজ্য — X̄_n → 1 ---
rng = np.random.default_rng(20260619)
X = rng.exponential(1.0, size=10**6)          # সত্যিকারের গড় = 1
print("=== Exp(1):  X̄_n → 1 (SLLN, E|X|=1<∞) ===")
for n in [10, 100, 10**3, 10**4, 10**5, 10**6]:
    print(f"  n={n:>7}:  X̄_n = {X[:n].mean():.4f}")

# --- (খ) standard Cauchy: E|X| = ∞, তাই SLLN ভাঙে — X̄_n ঘোরে ---
rng = np.random.default_rng(20260619)
C = rng.standard_cauchy(size=10**6)           # গড় অসংজ্ঞায়িত
print("\n=== Cauchy:  X̄_n থামে না, ঘোরে (E|X|=∞) ===")
for n in [10**2, 10**4, 10**6]:
    print(f"  n={n:>7}:  X̄_n = {C[:n].mean():.3f}")

# --- (গ) পার্থক্যের গোড়া: একটিমাত্র শর্ত E|X| < ∞ ---
#     Exp-এ সত্য (গড় থামে), Cauchy-তে মিথ্যা (গড় ঘোরে)।

প্রত্যাশিত আউটপুট:

=== Exp(1):  X̄_n → 1 (SLLN, E|X|=1<∞) ===
  n=     10:  X̄_n = 0.5194
  n=    100:  X̄_n = 0.9710
  n=   1000:  X̄_n = 0.9603
  n=  10000:  X̄_n = 0.9841
  n= 100000:  X̄_n = 0.9978
  n=1000000:  X̄_n = 1.0007

=== Cauchy:  X̄_n থামে না, ঘোরে (E|X|=∞) ===
  n=    100:  X̄_n = 1.126
  n=  10000:  X̄_n = 0.851
  n=1000000:  X̄_n = -0.173

(ক) Exp(1)-এ ছোট \(n\)-এ গড় এলোমেলো (\(n=10\)-এ \(0.5194\)), কিন্তু \(n\) বাড়ার সঙ্গে দৃঢ়ভাবে \(1\)-এ গুটিয়ে আসে (\(n=10^6\)-এ \(1.0007\)) — অভিসরণ মসৃণ নয় (\(0.9710\to0.9603\) সামান্য নামে), কিন্তু সীমায় a.s. নিশ্চিত। (খ) Cauchy-তে কোনো অভিসরণ নেই: \(1.126\to0.851\to-0.173\) (\(n=10^2,10^4,10^6\)), এমনকি চিহ্ন বদলায় — heavy tail-এর বিশাল মান একটিমাত্র নমুনাতেই গড়কে টেনে সরায় (আসলে \(\bar X_n\) নিজেও হুবহু Cauchy থাকে, সরু হয় না)। (গ) পার্থক্যের গোড়ায় ঠিক একটি শর্ত: \(\mathbb E\lvert X\rvert<\infty\) — Exp-এ সত্য (\(=1\)), Cauchy-তে মিথ্যা (\(=\infty\))।

সমাধান ১১ (★★)

Borel–Cantelli I বনাম II — occurrence-count সিমুলেশন।

import numpy as np

N = 10**5
n = np.arange(1, N+1)

# --- (ক) BC-I: P(A_n)=1/n^2, সম্ভাবনা-যোগফল সসীম → i.o. = 0 (saturate) ---
rng = np.random.default_rng(20260619)
occ1 = rng.random(N) < 1.0 / n**2
print("=== BC-I: P(A_n) = 1/n^2  (∑ = π²/6 ≈ 1.6449 < ∞) ===")
print(f"  মোট ঘটনা সংখ্যা : {occ1.sum()}        (≈ 2, saturate)")
print(f"  শেষ ঘটনার সূচক  : {np.where(occ1)[0].max()+1}   (ছোট n-এ আটকে)")

# --- (খ) BC-II: P(A_n)=1/n, যোগফল অসীম, স্বাধীন → i.o. = 1 (বাড়ন্ত) ---
rng = np.random.default_rng(20260619)
occ2 = rng.random(N) < 1.0 / n
print("\n=== BC-II: P(A_n) = 1/n  (∑ = ∞, ~ln N) ===")
print(f"  মোট ঘটনা সংখ্যা : {occ2.sum()}        (≈ 8, saturate করে না)")
print(f"  তাত্ত্বিক প্রবণতা ln(N) = {np.log(N):.2f}")

প্রত্যাশিত আউটপুট:

=== BC-I: P(A_n) = 1/n^2  (∑ = π²/6 ≈ 1.6449 < ∞) ===
  মোট ঘটনা সংখ্যা : 2        (≈ 2, saturate)
  শেষ ঘটনার সূচক  : ...   (ছোট n-এ আটকে)

=== BC-II: P(A_n) = 1/n  (∑ = ∞, ~ln N) ===
  মোট ঘটনা সংখ্যা : 8        (≈ 8, saturate করে না)
  তাত্ত্বিক প্রবণতা ln(N) = 11.51

(ক) \(1/n^2\)-এ যোগফল সসীম (\(\pi^2/6\approx1.6449\)), তাই \(n\le10^5\)-এ মোটে ≈২টি \(A_n\) ঘটে, আর সব ছোট \(n\)-এ (কারণ \(1/n^2\) দ্রুত শূন্যে নামে) — গণনা একটি ছোট সংখ্যায় saturate করে, \(N\) বাড়ালেও বাড়ে না। (খ) \(1/n\)-এ যোগফল অপসারী, তাই গণনা saturate করে না — ইতিমধ্যে ≈৮টি, এবং প্রত্যাশিত গণনা \(\sum_{n\le N}1/n\approx\ln N\approx11.5\)-এর হারে একটানা বাড়ে। (গ) তুলনা-টেবিল:

\(\mathbb P(A_n)\) \(\sum\mathbb P(A_n)\) স্বাধীন? count (\(N{=}10^5\)) \(\mathbb P(\text{i.o.})\)
BC-I \(1/n^2\) \(\pi^2/6\approx1.6449\) লাগে না ≈২ (saturate) \(0\)
BC-II \(1/n\) \(\infty\) লাগে ≈৮ (বাড়ন্ত) \(1\)

দুই প্রতিবেশী ক্রম \(1/n^2\)\(1/n\) — একটির যোগফল সসীম, অপরটির অসীম; এই একমাত্র পার্থক্য (সাথে স্বাধীনতা) i.o.-ভাগ্যকে \(0\) থেকে \(1\)-এ উল্টে দেয়, মাঝামাঝি নেই।

সমাধান ১২ (★★★)

একটি tail event-এর empirical সম্ভাবনা \(\{0,1\}\)-এ জমাট (vs non-tail → \(0.5\))।

import numpy as np

rng = np.random.default_rng(20260619)
RUNS, N = 200, 10**5

tail_hits, nontail_vals = [], []
for _ in range(RUNS):
    x = (rng.random(N) < 0.5).astype(float)        # iid fair coin
    # (ক) tail event C = {X̄_n → 1/2}: বড় n-এ |X̄_n - 1/2| < 0.01 ?
    tail_hits.append(abs(x.mean() - 0.5) < 0.01)
    # (খ) non-tail event {X_1 = 1}: কেবল প্রথম পদ
    nontail_vals.append(x[0])

print("=== (ক) tail event  {X̄_n → 1/2}  (0–1 law ⇒ এখানে 1) ===")
print(f"  সফল run-এর ভগ্নাংশ = {np.mean(tail_hits):.3f}   (≈ 1.0)")
print("=== (খ) non-tail event  {X_1 = 1}  (0–1 law খাটে না) ===")
print(f"  আপেক্ষিক কম্পাঙ্ক   = {np.mean(nontail_vals):.3f}   (≈ 0.5)")

প্রত্যাশিত আউটপুট:

=== (ক) tail event  {X̄_n → 1/2}  (0–1 law ⇒ এখানে 1) ===
  সফল run-এর ভগ্নাংশ = 1.000   (≈ 1.0)
=== (খ) non-tail event  {X_1 = 1}  (0–1 law খাটে না) ===
  আপেক্ষিক কম্পাঙ্ক   = 0.495   (≈ 0.5)

(ক) \(C=\{\bar X_n\to\tfrac12\}\) একটি tail event (চলমান-গড়ের সীমা লেজ-নির্ভর, প্রথম পদে নয়); SLLN বলে \(\mathbb P(C)=1\), তাই প্রতিটি বড়-\(n\) run প্রায় নিশ্চিতভাবে \(\lvert\bar X_n-\tfrac12\rvert<0.01\) মেটায় — \(200\) run-এ সফল ভগ্নাংশ ≈১.০(খ) \(\{X_1=1\}\) কেবল প্রথম পদের উপর নির্ভর — tail নয়, তাই 0–1 law খাটে না; এর সম্ভাবনা ঠিক \(\tfrac12\) থাকে, empirical কম্পাঙ্ক ≈০.৫(গ) পার্থক্যটা tail বনাম non-tail: tail event-এর সম্ভাবনা স্বাধীন অনুক্রমে \(\{0,1\}\)-এ বাঁধা (\(\mathbb P=\mathbb P^2\)), তাই এখানে \(1\)-এ জমাট; কিন্তু সসীম-উপসর্গ-নির্ভর non-tail ঘটনা সেই বাঁধনের বাইরে, ত