From 63e4ec2923bfd6d378828310324524fde0350601 Mon Sep 17 00:00:00 2001 From: Bokuan Li Date: Tue, 21 Apr 2026 13:21:01 -0400 Subject: [PATCH] Recent progress. --- .vscode/project.code-snippets | 10 +++ refs.bib | 29 +++++++ spec.db | Bin 126976 -> 290816 bytes src/calculus/frobenius.tex | 77 ++++++++++++++++++ src/calculus/index.tex | 5 ++ src/calculus/integrable.tex | 147 ++++++++++++++++++++++++++++++++++ src/diffusion/martingale.tex | 73 ++++++++++++++--- src/index.tex | 6 +- src/sde/exact.tex | 103 ++++++++++++++++++++++++ src/sde/index.tex | 6 ++ src/sde/setup.tex | 90 +++++++++++++++++++++ 11 files changed, 533 insertions(+), 13 deletions(-) create mode 100644 src/calculus/frobenius.tex create mode 100644 src/calculus/index.tex create mode 100644 src/calculus/integrable.tex create mode 100644 src/sde/exact.tex create mode 100644 src/sde/index.tex create mode 100644 src/sde/setup.tex diff --git a/.vscode/project.code-snippets b/.vscode/project.code-snippets index bcded9b..ed9ea72 100644 --- a/.vscode/project.code-snippets +++ b/.vscode/project.code-snippets @@ -126,6 +126,11 @@ "prefix": "cal", "body": ["\\mathcal{$1}$0"] }, + "Mathscr": { + "scope": "latex", + "prefix": "scr", + "body": ["\\mathscr{$1}$0"] + }, "Mathfrak": { "scope": "latex", "prefix": "fk", @@ -160,5 +165,10 @@ "scope": "latex", "prefix": "filt", "body": ["$\\bracs{\\mathcal{F}_t}$"] + }, + "Path Space": { + "scope": "latex", + "prefix": "path", + "body": ["C([0, \\infty); \\real^d)"] } } diff --git a/refs.bib b/refs.bib index 4cc872a..3a77691 100644 --- a/refs.bib +++ b/refs.bib @@ -8,3 +8,32 @@ year={1997}, publisher={Springer Berlin Heidelberg} } + +@book{Baudoin, + title={Diffusion Processes and Stochastic Calculus}, + author={Baudoin, F.}, + isbn={9783037191330}, + lccn={2014395958}, + series={EMS textbooks in mathematics}, + url={https://books.google.ca/books?id=ov4kcKkzTP4C}, + year={2014}, + publisher={European Mathematical Society} +} + +@book{Lang, + title={Differential and Riemannian manifolds}, + author={Lang, Serge}, + year={2012}, + publisher={Springer Science \& Business Media} +} + +@article{Rogers, + title={Book review: Diffusions, markov processes. and martingales, volume 2, it6 calculus}, + author={Rogers, LCG and Williams, David}, + journal={Stochastics: An International Journal of Probability and Stochastic Processes}, + volume={27}, + number={1}, + pages={59--63}, + year={1989}, + publisher={Taylor \& Francis} +} diff --git a/spec.db b/spec.db index 7a09b4dedd18f234e77f9450f78f4abdb679ce4c..06e8d46d4486d842a7620ded42a2fc6c0ecbaf13 100644 GIT binary patch literal 290816 zcmeFaYiwKDmL{e{9#Pbjl5JToyIgflv#Q9lB$J14$eK|iKn(+MVx0C^xilbKEi=$-(X43Gem9~mG>f=nlv05b?M z7|djn1_t?&`IB$$eO`P?l*_)i?!CtklIQIA+H0@9_S)-l^OI|Zx+dkz)orCNMTe$` z1Yzj=k~B2r{!fR7hTg+}5gFKT*fxZ~nW7b-;wH6Z$ny0)oRrH^l1y?*h==h6r3 zpGzNpbW8g1(`(n}CI6OE+wx0wZMSZJTU4@Ik-tWs($02Ps}@R|w~(znGj;CPUPbGj z@}aWbBc)!b7d5LoG{WAzf9v`+v)No3T`$$GclI9YiB=W%#G+EG-_S~`R@GF$q$+i- zUf9;`IVwsORbSdt3Z?FDtA(nTs~7ifXnE97%4xOkJcZIAxw`VGYkLRR_1HHucXdu^ zok1ho?$yh8He&HyIG0Q$Q!yo>W#w2_4o9?fEGjDuGh8|bF!`W?x?e+j?z(T9?s zl_Ml+gH=d+GKscN+1rli5-~NK%E_r@B$)`O!&*L-O3KM}c(bZhwhFl;-L3a8!&}!& z8tq>yR|=S^kklF*ADKIa<{z*Hpl8#+rrQzTR!2pt8%TUeG_Te5w!5(;v@qeByLd_% zVzW^DstD$2L)od9+4F|p{6<9oF;g5HnOj&8TH~zi489xrdTj%knlEg2K2Np{=AF*u zcY+V>D|TLb7?M7_dh7jW8=1-&?=Y-g51+w07;Iwa_nnN9NvnM`$f~v{+Dgt=8}`-mfL3y%Vfyn+0_q z1gE+^FI6^6<#~H+_T*5j5S9A858%&OP8`ZMi)o zb8o)sZGF~0SaYCy3%9S{xb^A9YdQ@!Q>8IU|MTh1s~^560S)3o8^xNJ^oJEb@$Z-~ zHmc=E1$E^eknVRv_r;O9qeq3-S=I?=;Bl_d^T|uR2-^0VZD_icT9__P%ub`9qxv z;M30b_2*fpN}JF ztzIpcHdpmpCDzj7h4cFB3+EMNsVJpY3GdkVHLWiF=kJRrhFX)$M^BwPckVB){$XwT z*M%AU8S;)QF(nt16De7VtLbc7i-c7rnM;JD(YzJ`{}ziSWi_k-23jqq|A%g?6S1fkSMza2O@!4{BpZ(> zR5cq8YteKhs>E`UL?k7L)A5vok8CWGQsQx0jpa2(%jVQ{I6ZhL_*_ZXStr-@snZ{S zs#6pQW99^oK>;~ zFnN0tW;=((k}3uC8rE#(^^$>1Ra1)h8)`GKz;;(~HAZ&UAl{?W&8k+b6&`BEJ!xB0 zYCBar*imaZRn(-~PHs!8Z$aD^SleiXo6?FD4oFIIyIiXyv5I2gmA14JeLyl-EP_2K zZAv)3NZMCB1vHFf%+iJP!1{90Kzj`bHm!7& zm~6ZTfLIf@8E1nsY?wLj%k^gF-i7mMjfe;vu>p(>1fp9xrP#RA+^Az-E&6mQ zAVm8C0jDSig|$of(Yv4tj8Nb}XnYQ6ptuPO{0;(@>dla}&^Bbq0SLqH6U;J~wPGE} zPdRJNkXrjz13zCb0uLD6Y*PgkFXuTuE*OaEU;UWgN{p^rFYk6d8x%Q4v|lH_^HK*{ zuT%?|2C%(4P0yZ>_*ZY_SHgyt*GqnZ>E5)KH3r&CDnd*;NLK3&A?-tA%R@)d>r2YDi1zxjC+m|)<*;?0_c{?tEz>1wzT!J|N>L`2<*)Z>LvJ;_qp@&2AC4(5`s9!Ya93$me&6fs?RmXuHsAuOg=mPHZ6@U65mv9A1(mOX2WhG9JC#VDiHydA>vb zx8#rb=b8M73H_Fc6|L)Un6Pij|1p2VMSAt^qJPQXV=yRVvb9S_bI}8v(_z=%q7F&C zv^P39xO(;_)}t?r9fQ~TrS{(Fj9N!;%>9$IOzlnHYqGz0@9BSa=>LSf@a&rFIF>@l zH#-}RPDYn$AkS~%u!ON&u^K;(F2F_ub5ZE0p- z!i$Mm^h{ftSuVO5jwavepB6jaI}Pn66R-D8OC(SAPYa*yn--6~);BGdKG8EREH5U* z^6`FYkR#98)51v!%RZJ^iekr2C1&kukto449$SjX7vqWeF?(8DxoGT2zqD{FJ=2jE zN$KsyM7nNWo_q(7nBN^(JiqNCu}}M1_nQ!i!=Vnl(Uh7_ z<`gxmr1Ek)mQs>Q5{N9RhV&*vk0cT4Ap#lfRqr!E#FQ1>GgkC3uJGi$gI3LZU&n_C zu_h-^zj^N5A1iO*BY6@ZiP!M4D&b@C zI6fAR;^Q=F=%4c88k2xk9bs|FXPOEh%(+Y+dBmmZXa% z=_)KkkVNWnN>Ta&&qs$~+?1PYT%^Wf6|^qVaVqBBK*e$rPu#u*Pjt6i0b5;)n&Ed4tkdyqG zLuv(DRIHy2Rx1q-D*lReuaS-k)w38{*DH;#8E%d>pN>pT49(7-{cQG=p_!+bt|_HW zD62l8?|JE_R^8O-h)h;E^#W#^ru0TZ+g3`XCklis6`&NS)?f}J--Jd`i7cIN<_cuA zL#VU{{fQDO+FjSwg3@SypDJ&#+U@f=P#b(a#?;t)XdgIrq-)>@_Y9DDE=^7jotizH zn|*Jy6@R|HrtGLNHK7jcKj)<@i#-Ge3@Fyzn<^yeN^E_nTCR|8@O3D#0c=>_fa1+E zAXMMms@E&ErSs>hiQ3|3xx87_7IO-J{vI@&$`4}?avwZ+eCy-br3Ms=vn8yHbTWnG zX)*$R(r6g#rI^dnYicc;jzhB)?`k+du|>3g{U)<5V@<#m4Krr+qlTr$X2Td==Y({q zzVsVr`?MI90aQGJdcD$U>oiZGegm=1Zu7Ui=PQlha@_LEO%FK)JOejPHAf#7%(+RU(D0gYb#j1%75P!U(RloA%w;*G_d zi`RFG^@0jRxDqK90T0ab7Hy2A(@7Iu*R(1$sHOL{f~pl^0SNuM12E?D&gI&6dGD)B zdmR`PkHK&(8jmI{j7cZLax$Eb+VSPRD%7-#8dUm960G+ubgAV?3HGxBl%PuZt>j@D zXdsX|jCnf>4Ki9j1ICg*7=k>%g;M6IcmhW$V_f{5fl!AU8812*hS@>R%NZm%trKk< zhV`FXpckG0h0$Lb!awJa0|5sD4g?$sI1q3k;6T8Ez^jMApFR}C(|@T{U{=;8S5?k*6EKX5$*YX^wT~Ixr6|qqS1=ZsZFqxpcPf#6&=)9ns?7e z?sa74+05y<(peScg+@c!#)mypO`uWo{c|1FqxfTc`HgT#2Hn$+mFJ7jHn5&6MQzv0 zaj7#0*_2wzH@cElBr>~GwWD!*Y3u(D7zD0%ftoUvXsUEDQZ8;3+gY0~TD7u%(%F@- zcFJ`=C;GV06DGozlert}sLz~5zDTd!rMw#vn6*Ktp$@q)doJDUS1iw3x)SVNtdn-Z zaNLHt@8PNMVgkLiWxa@}zFG5%4h`6mK4>VN@W^NF$=R-CgHW*24?B~Kg>9Vx;1{iF zDN`?#O=!Np*O|AbeZ~9Cs2KwZ)(z}Hd*()cwVjH-d>~zcm|v??p+=@$#3q16js~#3 zP0@$o`?o{*=lpRX;6T8EfCB*s0uBTm2sjXMAmBj2fq(-62LcWRevTk8?mj(gNSxXI zzx&sQ+`lILft1Ne!ONRmK0~gwLdH)x-ZFd%*)rg*nPa!m(?<|_&kF>1$qCiGv`1ek z0riqkA?koqN=>@Wy>fFLGJW-Z=r>$NF@JuB<0sQ`pA0uqQS$o4oMn+`v|YG5 zt4V=RnM*ZlBa7T~OniPV@oBnW{?n{A`9+sv>=g1dl72(HXF>6iUoD0A;0 z#%52#oM;(-ktc`#yS@ML^X1>mwa&njb;#?Tfd^gu5Xp%LLSy7;IHOQ>5Y86-TFm-S&RbqtXwhTu?>T#>K%;Rkw8y2kyUbJRQQUuEs} z@dO!Ihe9562Dqk_$eA5;8wgou02p1{$wJH~tsAIjzTqXHlve`QP#zRu*qA3w17!*R zdFdF<_j%0sn{WQbnLqX||Iz2IkUrnCre$OCXg->hp+W{Pn`%~0L=uT~Qp+naXU?TG zxaNt%wPrT0YPobirpD4)O-`t>L^>Ra8^Ey5G-jZErmq-SVE%X+bAP?p+@ozNXf_l2 ziQs`b^ly&=cMM;UL?D`QD;KCRCzEm@Wa_5!*w!Z#KqF@H*0#v$B^;rlLGrTMvxjPi)eE(wjT%)i>{^-y;jCf&Yz7^jC<{Hst!j>z3!57~apVbT+_ zI@%4gWAy~B4sIJ_Fqe}+9BU%j8}I}A%n+Ag_r=%SFXPobrkIs5W*^wbC!a4STa#_L zsE;V1cekUo4KZ}2sH}E)10>Io5qO5A9<(*0meB?HN+!(L5#L(1HjaaLyC=B0n^%s_ z)_Afzbyja%bMw+8gjLvLPU4}$TgZW#K~FjMq0xH2!j$G0!+$K`TWE#L-JlS&%Y1+uWX*eaLd zAlJsd+1MbxtxY(-W$c_H1#}*YoTpYQw44V)y!?`lJv8{lD~u{x4IA%+l)yO2_v~0^ zF=v5HrLY0QlIQ$Z@G(({#{r6X1)sY`;!YVK*1?r>qlYv@J2t2^K#%C_w$X&KOzRIFz=n!&(q+ukD+KR6rq^aFk1Mhm7BY zKv4$_`rOVX?}hZ)G*fIH@1D_Z8zSj*7pl+Bj!`Dq8QQ_l;A~fxzkMeAY&`jAPnm+RTmP z8@$>sNSCok8ZTH{Mlz0i!2@%LjCyUyzxMMXe6bqUb=PJqMcaYV4c4!IrJTbkKNB(% z^hyq+ICw<4{{QKr(bGe(<)yq7l&aFabZ4XbJA(H$oJcn3&YXGl=trM?eSP`2lH^H= z4xn&xpUT6Hy^_kt(-AeDPo&j!A|BClF(sCaC^1!$lc^M3|3?u=2-}K0$yeL~GE}Lax_Gfdne%*4keuLG?EL}_^ zvX3g{9KBq9p+Ya?D3r2}Le>xve*lA%f-lso&XW^!ehu=ECa-#l12*8+9IsQa4_)sfnXMgu*&CjT z!?Fb5r{VBYB)k|+B@7=*j`N$Jg7cfU0E6UK=Xe>MOtLR01G(;gyjLf~@^VG~7K>;yl{TJa!Qe^0y#UyXG&C?%1qK6wMko>y zcnv1`f57EM;7ckvnVRC+Z8FzhV8JHAS8*B6hMTD3Bj7_S87m#=0HGfo84{^$^fivcU4^MB6YK=`?0%>RfUKIVl&7O z0uOX^#kD5 zcl?Ni^Z%QZQ2clK{`VpH{>R7A-~9CXui}p&4*i2gco6($S2!%59GRGMPo2AV=H}w? z5lNUBW{*NEHaYsy+~|pufn$@iv-9Fae0X|nbaZxj?9_<+nAXjpo1R51#!=KXJv%O*8J(Cpa%%LcGO!WxC;T|;!lDnfWVF_K> zr#2xXo_Kri#I&dN`mxB7kX!Hs=1vH!t=PHf%J{H%R2*w%t|k`(fwxWtzkjoJq;(}; z!$|xyov#Gf8^VbxTPg9_Zfo)Wqo>>6u@RSd3BLvDwV5L$GZRy1r-J?yUQzy7kWO^I znb_ea#MYeKbz)@VS$sjB0J?lkUoIhX#@kY!e(4fh-@|{^B}6^q^QuRDc1;juF>&7K z8=rpLJ?0*FyJy`K_&??LxcwN~tl(WhqM)|o=cLi8m55j z<2gDxw-kAQlm`3sIjEuju_IdEDy~>Y{5^a9wZ-yLi9J-Y1TG1tBVmm+uI!HDT(T&@)eK z?R9ZtmwkJ@ty8U(%X6=l5=1#Fh>@ru#6t)m-lJq z*LFqm>p0->15x}YH7q^@FjWddrp}Ee@3vn5Kyam8t{ImrEV{zNnB4lMv7pN}jgF25 z-5!+*YCYm_+(YSk_sEw!2y0f_oEE>ijjC?D#qC19<`JL1>2dA3#7%f9_K3eOiQ@CJ z%TwIcV6Qeec68#{v3G?rH+>7PdACsXh+iKccD*w!KK8l93qHZqI=+sN*>$(r(nQyi z$5nEykFL)>INEy08(xTxox0fict&s?5xghxQkb4hMMKfxyIx^>!uWN3%*yGS_TCMK zKVCkvB)Hg4=E9VX3F7E7$lQ`3i60!#j4Vng-wiIkJrPWwc6ra?h2R~-E37*H(ml*> zqK*P07?tW4_CV#9zw*61AwGNE?a^8f1&@jk*Xx*DzuR?oWNhwgD;*Rt0M`YKA?O;b zRG>7)F7GgdXMhJ4k4sZM-XETQx+!{}R-bOVTvY%(Ah`U(#M9Ad>a(+=XX)8lmn%2o zniE{dT|y18tpJfqKP+JaLTK!Lx3>m=Vz?Ss=O`2QNO<_0_CB9^`oXZbh0h7GHR5t< zo?&)ZHZ6Yr5hnIa45J{3Pd@?`m=;fAN>844iBE6f)hB|G6)#LpPad5a8$IqF5xrC6 zbK^%xCXdaGj7^T7^Nx-m5vC@1y6{G8`RudA@_T`|jy?SlIo zCysht{Ne3|+2fvh2~#A!O)*cd%3!--%b^^A^=3KJfdS)B4r#7Cd50Of{J+3=_^=9z3QU=0ql zFK|1qBjVTv|C{6IKmGt22siy;1I7i{6qbaDRdFIn#3goYbo*-Sn9J)S`r?^6HVO(v zFFX`o8K^Y&R2utH5s0C<#K&W<5j`ldD7qE|*SRC!*7n5@M38ltEV^XD6KnmF3(LqA zb`3K{sH5UjpvBX|h)|n$x!pqjsOz5K%D8~+o~Li5T_>KG-)oiR@IFL zRM+bhLd`=P^orYCEb`dXXT)dU7hQMVUNU5-TRXhCBDxg!NQrL_;B_9iJYa(a*K6W% zt*(^nz~^6ayH0zCt7y1B=WTuV|FU6;iX>>xYkotkH?b^BRtLv+PZ5;!qF zwtVzz>kXn_E}!5!23QY!pT613xLqrt9mj@8bk&%nu3r*d8>06id5eUi3&wrtsOz@N z_1Y-<6dap3uyZBgamU9i5q;mCys!BmgmlnkI%d@K9L&pjLg3;h{vW~ z;^{Hr*73p1021=ievO0b`P{$amfbCYcG`X0El<4Rc1?W3-r;;Pg$3aj-8jDtyTuD& z&;94SULC(=E4sKoa+=}qn0Rt*=1p;S?DSE5yzUzFQi9-_n4JV|p+w^Nku11CaJ*wK z_Vl_C8#@wsZB#iqfiuqN;<2%*bK>a9nX&iB!`QyPpy?-hsxal@>B2FbEA$*L%#`?4 za1GO?ujU#1QMz^Md*ax0FTUo)XBpA;$TQ6N6QE#&P!L#ja$Km*xm<3MG%%*cXMXHF z6<`l2uirgdg?zOR6(wR>&Ujo6myq{NJUue}G&$@UE7F$>9P7k4RnPN14;agf^zE7? zQDS+!Tt%qm+BRs`C!kp$dq#90rJm{3#Ek#s#phs}+=&UCTW9>(6~TPFT>deadw7^R zwshg7!SD7e6ofaAP?50Bs;I^zFtKE$gDiXlLMhJF8fNc8;^`k6w^zrR8Rz<=(X|DE%{ zbN+YE|0~Y<-#Py~=YQ)|@0|ag^S^Wcch3LL`QJJJJLiAr{Qn(3|KmI#VpsVt{JPA3 zx$$fGIX(jXjvv29$H@Pm==(p0eE*;C|MC3`{BZs_5O5&iK)`{30|5sD4g?$sI1q3k z;6T8EfCB*s0zW$tfVRc!V^C>1B82OEUjCz8_l%08$Aqw2sN#xybOM@Q7?I!iHl6?f zt?%DJ0`Sj#|C{fB_WfhuKl1%U-+$%%FMWU8_aFQI#P>IRzw7(4?`z+$_;!6&U(uKM z{gUqs-!0#VzW03>eapU>ZvmA!e;f!n5O5&iK)`{30|5sD4g?$sI1q3k;6T8Ez|S-U zhTn1zO+P(89Kcsg96n25-xvMvp($zj3?800L7nQgdc~*VLaR&7V+@(AuKAV902LcWR90)iNa3J77z=41R0S5vO1RMxB5axL|$x8g7vxHh@k#EqMn#O;%_dzW`D%Hz^R9gFZs}e_cj9?gMo%Nn`!v!XxgYn+9`2dFP%E>Z z)eA)x;I(4>bk>WE=wz;3+U#g>yR7aMDUBTJH}E8N7KIZOBl)9hujDMB3s%8L)oha- ztZJoRgR3j{yf^8wZ)SK*ePFapzq|TXC{>WIls1bR{MqKzGWmF=ug(D=2$50II(~Kl z0RHK}G?Llxj`VufNX4&?Gzu(W8~jar>^q>Wf{M14`c^jIxY${fS^;F_d**o98UUVj zG^kceJX?{CTOuG!8Gr{s*l6_1kI7?85)Dvh;E38?X;nqFv~W>>|6H9l(O3`(vT?a< zNvoBv3ir#b8IAKUu6K7(>*!yPEL-iByx(2PR!1dWS@K0*M-kz&GR1rmfvlvn5x^Hw ztm>I6dw~f7a@z=BuBcBdf#ONl+_n&d*Li1@erlW&<|F0t+QLTol1qb%f`hX z`!0+86&*hLs`3DFS38D~ziCgt*_FIzPwwhT{!@GMr(LCQ+mmm1CBJ7+ey=k*OLSjf z&|Mu_RiY;i}Y&Qw*dTr5}bRyVWv8sT{Sh$XUypWe zl*&H5#b}%g$?P)`>shqcyUE3P$|@S`2^`E3+i&4*waPf-?vCx~Sr{+fC(k2&(|c4< zwJn9Q+a0q|*Eem-!I{MFmP+O7Ht!CP41n+el&Pa8%7p>T_rRmD00MfpUTg>D<#Ye! z3Dx!FNsr1(yfz~T>gmEEs>fhpSo@R*ER3cU?{~2=JhP5+tOsly{jzU^D<_tRke_Wm z-Oa-RszE)(M9IVS^c2G7UDvB$o{dape*du+?x&UYcsfvswKYRhbA(Td8B&IfUNNRL56rBKIk5WEtgy>>qW z_J(D}i!Q*__O`Pumu%$1o332jV7>H8uH$D%LsgnV*1`=;Ry+Mi6c5|L=&A3(Y9{kZ zTR|7|Su?jzO77Ym^Hqa&WwX~+{Kh48lAp=hlB45QnS6QfHk+jX?{a@{$o;+X|6nNG zS7cx&{LYgR3JZ7_VFB0Ie|+|9#TFJYo5-rsL_U@d=d%b97+17JTnVRhhzXcd5R*;Q za=BzWr)pXrAE9f~E>5iOODXAuS}nNO+F zoR-bz3iliKLuQr*q+KQc0TOp|ZM1<%R7sIp2*e zM%v>-t?B8~b(Y@o_%e_?Zbf{$wWXD-+BO1L(q9CdtXCn6TBS(jEY$gGzW|gv{mJDw z>Q8vUpl}m}Xz-AfIt5#^mDl5|1UMtnP!SjTp#nx%1^OEVj~c}LzR?Sjkke{4dmQUE z3KfZXL}3K$Es_IzBv!TJo-NpMzFaJpA5lbWh{z-}YFB-UML#MRjp&wI34YNL3sVCG ziv_7x-z#b>{$in4Utqr%>U$M!Nun?|$z_W1&3_?xWz`r$`-l-SQIX6!<|ttcQGakF za$6&2S2As=9;FONf;W-ig9FQ3NM#w?cq_P$Q8q z38Bz=%f%P(pFW=ulUbr&qw^LH*cqtL>2Dor4~UL3VIm(8zhG_LfFfC2zws1 z6;e9HnMJTPeAbO%*%Y|C12enOv1efHjP_7RM-AL`Ka9jIfHOvfThNMHlU872nQ*yR zSqcqn$L3~R#hB=|jjk+h)n_tB1NM?#WstZKtfRhKsP_^%kF}%W!bT|;HcRg|^;xg1 zV#Y8FpwV4{)@jA|9#oGdTDMxgFQcutF|0`Y^&JbmW|g)tYw9!bf#GLmH3&)+$h6@? zv<(-+IxYZU(M!uO10C4G2|^nc1g7Woc0!_k+{Vl;=qIpVvDFyK=i3dr^(Qrfw=dvK(w3jbg3s%}a1!9zeug;-(3i7x{vVpA59f~q0S5vO1RMxB5O5&iK)`{30|5sD4g?$sI1q3k z@bd}*vj2C^|39y@>x|2RfCB*s0uBTm2sjXMAmBj2fq(-62LcWR90+g-IOl%{6C4OQ z5O5&iK)`{30|5sD4g?$sI1q3k;6T8Ez%K>_q?{C#?n+^4UJ6Mn{_BRW**yssK++X> zE>X6%Qe7!Zx+#zZ%fU-ZF}G9PsUdeBR{2RJ=K;>?5*=lQ5K^6k$U%JGar%)CP6-iG;5eeB7_2#I=oT9xwILY*g`E zt@1o;g93*-9%NmoirPC*py`#)8e)U%mCF0wS*KCvRn&>Wa_}jYO0l4+QoU@BY2E;q z0V_5{23!MHCOiQ_4o{v(hvzw7!$y_Cyd%#J!F@JyF?giiER=FE ziN$)XKPpRwIx3vKuiyiw5LB(;RG-mDj3xBaD7vrWJ}>nTBk`t}88HHeur zDBu@(8%Fb#Rzqt4H0m0aBl@XQaib8Cp*#jlpC(bl6A|s2jv#Pk0ALYXQ7#%SeVa57 zPa)~k{-abj6&P9{W+Z)j+umZY!BP%u(;%?X8`^~K;4_5(os(Ygxs@C|wHC@^W|{Yj zB|EY%NQG=gD!?I50rOdh|1E|AN*1fh@P4o#qDR?e#PG=xeHYU~n77#%;Z4-eqlmTf za(L{q>UIH#J?XxRDcsZnwTpxhb*+RK)ToUZ=cPR$B(Q~b|8wReBiz5!@0TBGL@MSb zpaRHBp0e(gvSo5o#(Q5X*L$tXd3c{d^k6Ml=gq^l7SR|KhhN})9`E_uWh;K8)E5N3 z*L-a!x5bB$mrB@YpX@!sH{pQ^GJ327IY5l?N!0Vn9gv@(L>mR=E^!cjzJVlk;J!ql z-0TO7M}_*9^f}MUGtvY9)#m&|QfjDrQMy`^6jcR@Vw5%5UN!{UrdK9Y9r7ko9;5bi z(&bIOya)~>h$h%$i00VfFW6G90cIU_a301Qkx@5aAQq=ZgjL2^2^5ROtLc>5Drv&P zm};!{YbY~sn`?bb$-)BzOC^#`_zb7Rz?n^2BW2b@y)vpHg7&yW3AD z|GXbNL;V(~DDUF@IAp)Qj5t_3a{FSq-8;=?Z6^QXNxH)*|;R%N{b9pr~-b^WWI0=(~NL{I2=05@~Amt41e3UqeF>PjTu z-ItD@9AY-E>O2wJCbj}(@%&9tO+0KWrQX;O=>ZIQ8GB)1=v|;Sv<8|AhYFp8H)bJ{ z?lJVmj=>*d8m?5!S*=vqsr~a`{UhVJ*4g&GrUiOu9rT!bxR`t4seh4maL=Ghkw`Z- zkiWZldv|m`Mt`XjL!I*CrJWe)n8@}kFttQ6YyRfFyXagI!NMJYS`@{q61y{f=GlxFlo$Cegc5r{F?IMdR zu^smT@)mI@P#6P5I~Z$YTT|d_OQB_x-zhQLX0MEn5D&1&k&J@#_EeJUxr;rBn1B#i zg?soWUTFn+75pY9liu^U03B^uE5u!}4D1L4ceZMI8=MVT7;F>dic>0A=!}Jd06y@Y z4|eGR?|iNdPN$G7ARmrh8+uKP(tbV$4+Tvco4tz+`ulg0n09?=7!W3Cvx>bBycmvD zHmEY-BLxob!R3}Ec62M?fMCp5m=BQyhnPKLQ2G2Ex?aW-*!pqZn%#jwrE{@|!9T*Q zF39QU)4o927qWDEwy=(4I2ca&unY55p#r+N7Ky7oIXdy-Vv^@;M2leI?|B@ zjZ#9KU+vpvqrycKuMPz)qtP97Ap-XABIBS*@V5hoS;gOTaAd{1!rrNTtk{Rjz|94H zJ$y^Rwt?}Wx{MJs3`Q^DsA#7Q=_4c^@MMa^I^1ZZ|Fm@aOB%qw=V_T3$2#z`G#8x* zO@4yEi80&k5x}| zL=|4nBkdw9PuPM-{ru2g-<3kOR%aqI<9Mg2+XWCSgX;$|1FuF$HN7!Yj3I%-Zxix? z=)q~OQbaXIDTGET@Z+BSZ~^MkOG48%2Fk?IwZ8EsMpM#HzIBX|5jbM+jdK+nKb0U* zlITRnBz?)SwM5OK2A)#i!0QEu#8w8GWziFYdVZ?M%uPZSqq;sBuqx;{jg5~90=9^c zO?b*Fe%le}H4D;b z7!QZQpeyB;d8VaNP~Ts)rnk|4uW=i14>Ruee)O)vg6!Mom`QGFPzC!eQ4O8#_mXkm zV>ocQ*3YrEatLC=g%7k43=&+o0sWxCa!Qgf;1sW6n=}Xr_(c+8OWQl}u8P%t&- z8DdaL@x@3~T2pq^GER>#xWAD^h|v^%SFwmIaFH<561%Z%w*6$j|1NH0NwJMAN#w7U zF}jQ(%wsBvRauj66XnE{A=x8g3^_ZbydbDRTw;}N;;Y&dfnNsEeBNQ$#Hj+nu5jkALoVOzIPkN~L(pWt3O)kh=47m&r{AvmC$58L z-Mm0G*ULI4hqD;Q1u0FJ{o078j->e2equ@sfkqcs!Tms+QAI z$$T;$MsR>=KAnrCw0tBHO>0Tw{+8(c&7{QEliE*?%t zlhIT>nv@es4KyaL#+0zCeJg6wgMb)!Q(!$K>7}?#cT-6@&ZVU6Q@=EE#O2=sak+VS z;#?cquS?}`aceQywUGgMxyNN~wEc?W5chSq+f!ParT`vX4HGwY+?^Ts#XXi{^jlon zS#g9q+N0aDdyVD~Pq|Nkn63Cuol=xZW%-BW#QZ!5H5-bClX5JcOUp61Vh5`ejo>?` zDfxt)P2|%r?LdArJ8w%ya2HqQ>~2XK)l2+ zoo#S3l`kQt>p{+5MNUp4qYH9elEX`p_)@!^oga>aI3&T~sCmS5 zQe`=&MZvtp5SBqr#??eL5e>OeUmFm|? zo8n_-Q~R`3g7%9leWqohDtW+T?^WA_v2oU~eJe%yi}vm85E>XiJL)%BzjXc=M0H41 znIM#3@GbD`?tU?s`rbf4ebir()wkV|?4 zc@{h9(nceTh)WUZ$QqubaV4M2<`Sv2tZ4aoQq4y5 zsYq0owQM|zFo)R$B#SKSie~9cEtK;MmWUd$SL8(_+by2EzksKI%FBq?a?UDbkQ3S4Fro63fDj~Y>wnyTk2xYkYHj;7G zZNw>Ss0L;rJ}n}a>9NdubFM*8ai4A#?cwLqR%mP1~v4-8S-;zZcvkt)fhAXB%HB!!#X zZ#9|MpNCO2BRcFTI@*LaAL-hXh9uS0@{gw$R~d*Df`cX-_bMKZ4ONX9{fjo4^Z zoB2?%N@^7p{aCFUaZ_>T?vr0K3sk89RGfKs3SWBTJwx=a!I~bb zV>a{6K<5|`d;7pCX$$AeM#r2MnA1=-U|J`LQi$VCD3uRY{;4fS>`y{Y~-7qz|*AJ1z_x`IG7uSg&s~62Ef0bOP+zz z98m`VkyyhC7~>79wUDv+@iyG#MZ!&U!_i^R-=owQR7A2?^I)a&D8&}3TEl!An^|Yw z8Qkqa57eZk)Jynxpx$Cl@Cr~~t?jeM2b*brzBj?Ze!xqFIm*ttW*uRZjjpiCHhW{k zw1fOe%pcgfpX^-jJ1cn$;0C7bGwrtJ%7KVQ3XJ)jq9(aj5r#N~45>*p`|4S7tlPD+B|m{kmMdbj61SH_c%%S>tVR{2%vV99mKH*F3Y3?{d#bLT4j~(mYJIN zRU^MqYC#2oeBI3>&dMy0SyH zk@ymqetM1A;OGX@;)`oL6|7v=3a$guR=r7{MS?paG>xnn=p#*~A2t%mwmW^&N{!F4 zpGtZz)@lhyBVfj6b8tq5i#h}St0I)7RHJIFlj3$6W;C#zCTm~#DG6X15q7dGP5|qw z<+Egf4HnBH*uVME{f3Oq6!7VpN5*jF%_iU4YFK->h+kdFC1=$bW|6Sj?gzym;v@S13b-G7u}*=5>Xkg*EYdWyrgD0e1K$gQkC_ss7bKrUsH# z)dCmJQ`#ychiq4{U)}Qi|Gyj(|I+m@b^Zh_PPGZXr)siFZ1S+H-6vj^B}<=EdM9U- z*AM+*No`yWjxjpl7_w!cBWZ>zBSZ#9!90^TPWo+dX!1Vowx`+Z9i-uN(%moaQOi9U ztB*?i-v7gdAs<&GDLJ2tYhhW9<#Ty0nT ztd>&3`EV2@!@2*@d_?#0bm>=Ba7s;faZ+HXW1rD;qgPmsg|NpFR)9Ls(Y3OwP ze#p`3u(5lG(&>1)H_ity2_gyd@}QT^&}o|g{eYKjl{V+joO$%<$8UXo-KKP)D#;X< zWJ*b82^5|Q@(m}x43t|P%rCN%3Y!KGLiw&Y*mw6*Le>Q!( zYtsaHAyhq-iU4tKT|Q_Inmm2duF-jI;L$SccKb|x2fX+6zW+}<_y5lQzjOcZ-2cPT zr*r?`)5qTn+Z))A^Sxd9wX0#+tachrI{&-ee|Kp3Pw?UV;SgZwf9Jr6`_>gdE(>0U z5rBEZcquB%uj0V?N40S>`bfj*BOd>)n_oBCf${kyQ&WcpR|@9p@mMk!({MCYQqe>% zmJO>)Hm9Vs;Vj%{N7O__j;XP1I+>Isxo8wNpJ_Rz# z804su5HreBnNO~lAwNi4@&oSYi?z{VD)e%vqwnxu=5ha92`~e=X`)*wIu}1M&(9Tl z96&D6OWU3qAoV!~w~MH3Cx2^HIOKDF;DTVDe!+ zk$>SF()K+6C*sce-#Pa<=YQw?@0|brOW!5^e|G=x9vi~H|N9r<`2Xv_Hi%ixO9*>DfJ z;qH6>2R`5+m5U{_aNv~I;)%Q#&dXXtNy+hCI0_L!BAN~3b{M#kPsMWi2sA?C(Rehf zg;R-e{Nzumx65)2ff8tdtKWbK;NB%nsvxI$Y+Tq;xW!g zP(1IgW@@|Y%#|sO!LEMQYx$!dz#TOOHQ3kxLP4t4_lg?b?$_!I?Ds+)A+?qyf^Tw} zjFI>+rVU`Zv<`3!aegp0E|&vLJ(?-_*-Xfxinewe>Y60pYe2=I+2&;hN^@NOtXNpR zP&}|d@Ri-sA99-=(Dj1?2kVu$*{DE2wU=gR?Md(n4dqD?3jT1p&41TJA<{pKURu^K zrcGtf^VIkAkii8zui~7_j-}ijf!WQpX}K{ICEZq`+I42_cd)l` z&UR|;wTR7&>_Z)8akO{bsSoTWLq_M7^_Tlb z|Lf2VxK4(ErA>w6QIhXY{r(?Ef6oo){O_FqbzE}%@;YI4dF10Bkwp6lv3hQFzj*9M zmqX-VFc37F|Hat}Qdqb?clPY$(4YDK_ye25LLSBz(MS>&7x_pG5gww^TsW7GDsX^| zcrcoril}iV7m20QY9gOYWMf$bkwuK~Tn@B3p3LF6T`9xqG5IF6&b$j`nu}1ag&gS= z#B7;9w(`FWf!SqCU}meC3C&`!RFlCrq@Lp@Y02Q`G5gyud@F^9LGJC_uIYdFJ+c}y znV$L&IxxQ*ZjBMN4y>aQ1`AG)2bnEHYGw0N(CMJ$ZXcV=dzLvU{M1;6xr2Hr`Hb{A z+8rq9rlM^kvp&rS^OQN%w-fC#R}Z0ga`z6i9VK(}JggUCs*85odc-^eT!v56i@|Cr zcsU=k*)rRl1eqY=z8Y+fjIrL!@R3QTSp~y1%fAYmQ_e49)(We7gjc(vfBz7#Y;1>% z@M^blhJi_d+3B`V1=N3-7&Ksmhf>%n>RJ^vYVKeCU8>a12|eFBc+B{-!&>B(oAcPf#z- zER~Id%m%un6{&YlYSoGK-D{17?4ZF#3yKSy){o=l(%o*qFzukpR zpY^U*ffEPPl%76c+Ln*iPukZK%b1q}Iga~zGC=PDe=pB1U2${=04^>*`r!hnrU<%A zfMk~YxbJ@p|L3?JN{-8^tHBGK)M5=VFDOOG3*8@^RlD8s zYle5VN8NN{qY7;+hTM9?-*0GBUoF*8vycDcSAhbEncfoxG9Z$5h3mMP=1QDiSSnwK zkuyY5UgZk<)O%lRY2tuim_!PLfkjh|-wr87>PV z^}P+9ixOZ+ymb=$(P6^f}XcAE(*2!Zs(aT#BH(tGB%FVn@@dHwrV~^G5Nm%6j^X z<9-GK8%)FiV{g+<)LZKGpS_>^Z0<3e`QzQ+V>aB(W7eeoXN&3okN)Zq{@sEs;7`3r zUK{#TH$GaE=a|XqpXPoyx%_)MBVHGxS0>eXKBg$?TpTeh^D1KBCSXh&jpd`MB);>p zKRa`TEU9_ziC%`B<&iHI6YP$b-JE*;U5q&}SstC4gvp~Bl%HW!IRqHx5O zPT;~vLukxoPSN7wFf3Tp(UdB4qfW$mE7TXts)i-B$ltv*OVtmoTOMwy%G!9P_d6Yf zB82~T*E#F+(iO(^t|=vA4nQujHNlkgQ&Z@=j#VAFZ9eks2<}yq;LMe>X0%F8TII&6 zx))w%0UBJO6v!^{X}WDfy*!ZO?Pcjo2yT7Jo)xC8asX?BvE{*C5f~Y44?IV>1WGX48iN&w%>d z>~2MEHdMIjq$#l4ThYBMllaV@V5v(Fc)sYvN#_f@2(KI1WB9}{3t%i#gaE5O7G|xj zbpwUrJx^i`d4RisX0_}Op1~kuW26BR(=Y!n2Tw?E4R+A>v`zlr4yq9Cm$poM&EA!P--E0zWyNtVj>Gsvb@k&?xtJjoLL~di2=-!Mhxew&e zu|vPd7W{r>j(MnCH>lnfn}2Y9J?bvPY2k5a;IvjO6xEFgbe4)yJgt<|7`@`us)Vxy zMZx>%@;a#P%cG2(00Ui#EPpG~xORB5XhF=xpUUdXp^r)*2LC_;xn?JhJon)!WVT<+ z2Nrlkm0^G0Y@~kXfl~dz`#saVsRx7Kc%v!!D}hJe?$}TKdi*`sMSk2jOqTR!OQ3` zZ)!}*O;;i(*VAj@>v4VC{pMwea(`fPRW|>j#9ATPsnz7GwKB8Zu-ZWK?YXTxhI|q1 zMFkf$!9c}eexz-sRDfJi#YPO=*Inh`GD_e%%b>0^NsLEgJH{PjKMAb&@#ASts-$3V zVnyNX8=*6CW8X$D6!?bE8XeJdvPM2PdxQ>|t8gF+c`H*W6JK--b->rWeQdUsO$YfB zd|H;b%T?Ul)V9f34Uk|%Y9D|u+Iukw(6DW$PcH*uonwG=41i)GX`=Y)%>MIt48Z>L z2L7Kp+-BFuIMZU!=@hmb_}6?HjRTdxjFBG)W(41s`g1hwAU4x-la-PF{y3sp#qvVuT? za?`9*=S*L+FiV%2Sm&i)a42vI%ecpWR;OWLZ)r)}94$CO5E88YaV# zCA1f^17;QKD5r*xAc2CI$%d% z>D>iRnk=3XoPH?9g?zP~)k=k(+P6Pd4|IpYgjjv9F<4kSh|bvcgWZRvb?>k&T??+$ zjYa#Fg`v4lu27WH9;nW4fi67sdt!Er5VY5_PMfQ$mCAK(pwn~aPC*w-k-T&rvy-QM4h?yv6^`eAlak~K)sTER0CWw*>0B(cc0XUO|peP(XpcI4rW=p z4eTSM;Fl8F`{d%hrY?!ZNd_vSlLpbXrk|a3~BgF z{?0;Ui&r^}LIzip(4*Btfv{F&r+N^E8Z0+77DT4JQ@C&2n)2&Nq5xnH4W16>YP!(4E%QtnQ_r7 zU#k@!QsiMP=AmI@%;U3J-VCLhLT(PF3IbwuFy@2nD6=ahuzp4(jv+HH8BM=FhoD(D z!(GRa8SxWs@MtrRbqtx?O`aV?W_BsqQ;8XuUmcdpKj)FOzYB)UHjvxQ4PdIo3=Z%e zI=~qOq-6+p1u@paXSUw6-bODBef4gO>lf>&slvV>^)TqyL$Pq?@yX1j%(S1n$0Eb`?M=VVm;vd(G8!#bZ|S8b^q#`vJ;7M-AJ7Ohtvjy9h;?hn-{oLrLtKr;@~P0OgLxnNHv?eex%X!;g5 zpFnJw->8v)X66+_H(<(?%GK=#<%j=iX2GK}e;$@w7Tx4y-xiR#5Uc__YoVSM%I9TR zE164iMLKNjnVxU}+V>c##E06bHTzjekrowstj3Tr@Zc^<_-F_^fFVn3uU5aa>z=bV za}+y-zRe8b3dByUY??Rq(bw$WQm{Jg>#qy- zFpMu?1j~?sIStx(C>^A#U(Z;5)Hm8fUfkN%xL3t=aS@wg?0mD^T(@Lghg~V}R0~>F zx>|yOw+1C7!_=xv6}VE{w^~c`%tYmEG4JwzSmtiM@q#=pFwCMkB=pj;3t)y~t8MXS z&vDO!02wR94wuR_w~#hN?Hd%AZqeA} zN=mj2>s_S?q~C(-Xw3o~Ef_X;Yq)K|mAow)o84uTeW~d7@TJD*_t|o*+pOXTw6lbu zTzaGxi-&R%vsNuXDixspdA&^R@{4$$Cm}i}3>^-0|MM#7vmQu}qOsxWKna9-?;@?C z7Rc}c*X=SkS;+cGXrZ5%31a=s_a(^aVSHosgaOSQ-8LWHq7;N;+!e$2L*?}vA`0n_ zR7cwdsskzaRp43#oYe;Xyf4=&>bug@!F&)lDrT)~(j^oMN~#p(uHa<54>yi@VnzZDb+A*+nfC9Ay%%GB9a28~C%Zy&%DX`?2+&2}C3y2i(4W?H~l z^>HrqY!m7lMXRq)aDzaC@hZB^d?Sl#< z+XdUORDgzvThGBQ&;-f&+t!jcK-UdgXEfIdzXhYoP9T|066||4%XpLyh^bPQ%|&}z zKBOw27(fgIr;K17g3>N(Aowj{*tgMEcOs-kZl9i5&`SG zZ`cV(zkTBK`$l}rK4e1Q1f)gjV{$}jY*|ta7999@xscEAFh*sON5F_K$Z<&yBj~|W zJiaI=5>7~}Lx!Y+9avp=$*o)Kz{|4TC9}B;@jbTP2lKtA;aC2r?%*CF=grPUd-%B5 zUD@A4P$K`5-wDL9g6-O6XOAAB1$aGB_7*g|RZ1@_X=4anHN;dD zWv@fJ>XoG@GN8g}uJ-_Y+DaElIAIt^KW;OXW}PsE5wwfiv~AHh%+LNX(lv0lSXQ_& zI1orZWm0{VM<3fvsUL^9j@p(XZK~S}?KTl!n)GgLqcrjx00xWSIRVhRW=cBfz3w=t z!w|bit*EZU4Er%urU_6MZd&8EpV z#GbGKyIAi#C}UP%Y)b3~KIl-$0K2VAb~BfL-Lv(s3)twq-5pV*eHKIa_~bm+A-P3^ zwO=RfH6{D}gznkfGQrO$xW_D1R!v7SOz*71eH_0Sv^aF^y*qd?_S^%F#RPw#?Yl|zwVfqd`;w`;s&)Kt1rI5?B z#q(lQHvppVwUHNs!=XsPiy*EsbQ`Q~E3Xq>?piJ;3~Hp)kahwQx@X~lJ}cH-?o+>B zis*nbojvD7_f?I%5_V>(m5Vz})&z}7>pW!am)PZ5g-nkVb4WSbJuMjJaq=!p*ha+_Z6UhswMpsw5^&XEF+dIvL=Y0bn~pp>Ai2_$)LLHjaV}KUFL7N U|2 z^jwu$4gQL=lcpo9U;?}AII%= z=@x`~SyCf2ZTRt;B$m54(BI)F@)Gg^Ga)_K-7k=(KBWVyKtD?`(QctnIX*`nmYCue z#B#UrjyFrl-ey&gGdZL>J=)?r+xzy2X58!^$ynCW_bW8CHs4Wr`6)}f%@a*eQAmJ; zCVXx}I{5&w3WEO316rR-wAQUEw)NphF!r-s3O; z02}gMUWUY2MtWyUumr*1I%=`!g6Wp+O?w(ZWS8@Drh1QcoAYQ!c0zpPH}A9#JCl_* zQVHtDJpBLcy=!wD$C@Un%1RUpcYxqcBt?M~C4mx2;SNHwWJ?q!%aUxJV_TLjSvt5< zuql!zKv@^_oGwuE>C>lYqhmW}dV6~&IyUAXOzelN_wu|_Mnxl!?#fEZlAfSa_C0F2^;Ccaq2c|X zOCh42gT~HPl9jjLdmDX6z2zi>B(0MaMEZG}4sChf>(FpO`!^UXa-7Mr*Jrd`>ZrCH z1ypr=0>&}W9Fq86gXR)i+Bt0OR2V+7QfC33M6tVyhPX9?s7Plz>U?@?k+1tX5t6#0 zI~U1mo#mK=2GfA_|F*#2u5kakKxJ7YzjJ1Yqg>~U?6F`?t=?0E3Vl_he&NYu?iOqp zMf3qR^rxLrQp~zQ;#LO|R%|In;#S~Wt${1OmSGc*>|oi_D-I2yCqV5KJX`YW8)($3 z$+L+oNUpskStZ=y7CKU|Vz9NH-kf8CmH=ynU~BKVvhUBw-gn%macw{Ap9hPay%pPE z8)}hV47F(?b~kppPZhp&q{YdiE1k1?&Xj^r*_EKz(Twa`%bdNgN>Lu)^%q#mI*24yEEt|H-6xw;fEhP~} z=vvVu)KdB_2kRxVp3`bhy|w0GyqCl{mHgwGAnClQRx@T@J!`KL0x%PK1r8|I8z6xq zM~_8;l)kiqklP1(rM*DWGvc-m3U_mP>&;eDYinTLHwoU0w_8ap7HAjYWPw&H*#okv zog|n?^W(8H2u3>=!DBo32iQwmIqdC$RqpGz2UIVj%R|sfaeEKK&9L7q%_Y++$HYqe|7@T;ur?~a>vDc?DD8Ihnz~ZU2cPDdui2+ct&?kNy5{w z7A7Sm@nt~}X`1_~{RP>@ND>I)5nkCzg!$N)sKi5zK!}$1HIkk>%b}D#gg_wX36Z;B z*#|3`+9G+K97fbt)ES*MusXCm<-@L&57E`n1byi3`)fv>1PRM)+f3;`X)E+ebYk{r zua4*(BeDT<#93pJUV0Ex(`#ktc9(e14m!m_hGr>6ctc*n)3(d z4ey1UY*d=$G5C|Rrtfl{d zzDwu-E4wg?B~WzBiO@vST(4hv=gj6hJ`kGpg46F9R+s(alfKSn0eb-vFalOXKV|fb z_Fj1|8M_CMQ~0*yvuO&*VWlht^=V$`<;@@r*}MG}Y@2x7C>chQ-KP4^hh!ev zBS}Z`wnoL7vsbZ;aba%R-@&1$lEvP=b?Y+(E@s=%-}&}#=K-7mnzmcTeR8JgXoGNG zK%YqGknXx8cf-!@@baK?9OmXEqV0-=l13yHdo1j8Vi6}FLihQh!@+5PaUe=`w%gKj zAHbQ!N#1S(?c(8&?Oc*+3_N|{IyCU~c>`ap^8-(xjkV%>)j|v-d8}~Ek_=_w=_7#9 z6TQ|~p%R9N5wz6dcNpO-cy*wD#4UNpX@+i_{zFwoArPw#(rpHKJ>Kco0aGw72ag4C!SU)cT+V=k<{@=PZ%<=gp1v{j@sevOrubdrQZX_C&qrgP z@rk+Nz{tn~_h(&Uug5du@=m(dVKw9(=DtZD7~xT0IP6v<(c?4Y-YM@9logutvJY(d@A&q>$vIRKKzZUE>#=jNadp(qaov*KN42- zVU(ua2pyRrOXnfH$`$6yY+~Q7%XeP;UbgHDgu)|%fZv0X!XrLDwdr6?HkhE zE_cuub$R?D?%|^oKHW2x?}!%Gyg@CixqU9h>!Na$xT{A= zN=)r4OF4@wyrMKE$;VHOjfXX3YB)JGuPWNy*tl}R$ed1}_qlWr_Zs)E3&v;1b)LGV z#Kybt_;+LuZj7mH-0k1XE~fk#F1P55Daje#D7L<0+;|uND@;ji{KdA$_iiXkis#OR zLcZ_`)vNkcbwu^!e^AxbSu|}#(HHShsTkSQF;6h>^0>8eze`}>LV4>h(m;Xzx^olrv~d~RspJ%t@t-U+?3BtGi9kus*D+R*tUYYOu!;W6`43Hqic zL*pOBm7v%B!6(+3X2jlA8OBYDF*lzmf8%PUJ)U4NJT&QJ+Qn~NY}lA--%!>x{f`ev zb*_J(AJv)C-;nsI50`=5w!O*bZuWu|>; zEGs;fR(LY4C>hRMP3Hgp`itYnIqo+m=vz^;n*QRO))CH)yNvH$AK^v`BzVY-=x@lm&b=;n7}Q5QSu;?F{izY|h4V{{cC zBdaPms+^tG*al-aP&nN8IXkMbBMSfH9n|1tKx?hsaskXj&_(vc@L9?}fOL zMDr6&*$AU0?2M{ZxE8qkj>3G(TyQ*o^~UsG}ft5T^;*6YMj&)i)ruVMdMaPVMB^Oj+;t2u$-Px$4}`>*l+zB_1Y;}Sic)j+&VXO zR$*c#W9Txk!ae7}=FTcHetq23GLf_7+?3mj-cFicn6*PdogEqvOcMHLD%K~)> zJvf7Lu%)qTMbrQKSZj^zt?kwtW7{D3oWf=mf6KF1**nIyvwO!FE4$g4!bX|W0ByH0 zkT(8$1067rm%Xp*jmGAAt*m`6Rk@;t1odU*p!V6s5%935EmFh}Sc!ySX0pjrpSPz;MLv4R}uL9^a4> z^vmbUv~lj(Uhdr0x#PpF55fNU@Ob;T0?x!oj`0)5*TZ_qlQO<=kGKO~{kZW-m<5eR zaZ?>Pjyi5n8M$#S95@oy#%A!Jdwq0??!bkDuHzn!$%_+rljg|_Oq{)&j8)cVxN!l5oni{6 zt?pO2=2w>7qYJ*%f#Hdf$aX6k@OW|~Lwk3Pr+7H3aNP%X5*Sv39`ml|@pu%!CKB_Y z=FfUs=P~45D9q(iyjs9m#2j>qFYI>M5cj?_d(?O4)^#MHu{n!vfKOpTObL!zF&-!8 zk{R|qx@HVBT_gUYMTR}#Ky*XfL@+xHC2uR{{lhv2Mv?JnUgkETKsaZM3OhZd8;>qt z=V0qB#aT+xGRAKhrV&dpmk5RzM!c0|_W|V!YuzL)HU29Dj*!YXhtHwrZ z^xfmeHuIga5SY$>SEwt@4Z6&Hr|A0g%$B0N=`%PaAF|<rCGa)$4{$Jt#~Cp=eL3~=KLXGPV$Ar}Yib)Fd-wm}Lz%w3IU z5#@5@cPcxgxwi4*=9q4L_SFZREhuc2v&-BKN0yz*>4o~fa?~f4%&Su%5F1u*n|4>zn#K( z{CbSP$YCe&D^=gzt_U2SC(L?9xSie>w}b-g;Lrtpa39}-yqX*9eN~EIwUJkAo zQh%qWR0E_9tGCsZe_Cb!Ux|A-Uj#8BRIcLu;!^oL*q+ax*>`Jn#gTPs)qO;aUoW5V zMvn3k?~y1zrkGc!2a4t&34piIL)!7 zTw;HMA^4nzAV$%IPic%XrgDnV5$1cdSb4TEdcb*SRnK<)F_y!}Rob%5XzV#tYMQ?_ z{TFVZ&)0Pv#P4r~dY^(*kKTbdhQ_Tz_ybhfp%aY zNuWs=fg?uO3%GM%xSO38z(VbN71KAd3h;H(tzfkSnBbMyfTqFJzq0u#Rw2%2vauK( z{kDru7zZrG?l`g7URRdY*WcEg4`Y2G;87--qmma&N!W*yAJAQ4oJwF&^JI4_JSDA> z1FUzc5RpdBs1gW6fMu%v3R<)v#%V_{3t%RWO&X=F;4a$bT;b=K=jeY?>VV03V~ z02*KouDx`Gj+Nww8lk1yYm|JMASIHR58`HeDyN;A+!3wf`Xs=tfa7-C0NX}RlzuSI+ICk|cGeddCtskK2!)K^ zN2#QxH~l&Ktd3enE71X9;9~PJ@8%e(orj?h)i^}da@W-VKi?Jle{J><{w?+jxZG{>{W=gjJ6bW8`a$67 z>iOCAdZW46Xg*!9E;fO0c{WBi-}!T7t1f>bXfNTibwJ_4p?qFgkoHBx>+Gz5pBp}W z$OxPno18p)^dH{+L*U%MU5nt)AzhD@(%EdHoK2>))pWU<&lWPNYO<0mEGM(cN;;7) z=JJVRzFa6&YPC$VoGYhtwOYQGPnT=u1a4zXz4=eKqoVxi1 zxl|!nO;xg)LavgjCbPLjK9MSyO4Up{Q>mm&+2V4pS}LZ?xlApcujMMmTDp`<6|;qO zzG9m7oNb2=h9ix6?_()!ZhOa^w0vXQcO3pQP6w3|jzV8A(y0@o=IXW}F~VS$*WtfT zU9$1@BJAA(%0uU?+1NwiO`o4V(>!ytR|I7sV#ea_H_rbo(UjZ7pfQp_RQjnLH3kImrE-)z2!p8L0Di?tN0<>h?2m@Ou9 z%hgIgmrXAhmJ8`@x>zk`5|v6JQ%&JJU8>cH*OW5pbgH&oO67`aoC#`$3d!(=to|&7 zq$Wld0njGxCzso4Oy|tfde#2*HYWWC_N2F8CU9^u@9kIpv=iTMu!9so&*aBoef7*2 z%;8BJ@b?aeH~4@U=rr`+DOrrQ2~WsUZT5oQnFM$?C$i4HbEb6RHNa;wrI=!$w|Nm* zk&V!Ckp0Lf_v)L?A^_7+262XVNuS(tIE#_*eg}w}kTvd`P1FD>Kyp;>O(;OkHm4-H zQQG54aOb~7lsSip!1U(NEG3o;0};VML@*E${9+=40qK8c(p%zyWDj`&_CTF}IB($l z;8mCxGvJe?jN}PI-VS zk&E5F_9o~pR}fsLxRFfBn6swQ45U^Q3u@V_#de@}FjgCLpmmJvKHL zIP{;&|E?+|7JEm8G<3O`f*`aADQ794$rUPb6yC_~?dIO4i3j*tpg@s6KV>e$t& zA1`)T1%$W|J9c1*L(uSP`!FtqWVhRwGI<+SoCIjqtp-7hlQ3;f$OIt3{CYPIkvN`^ z^lYmEA-RkqgmTSrLVz5HKrTM_wI$?~Hb!%CA%%coo#Yr&xNc{%BpP6EHgs!O#3NJU zPSd}<@!&Zk9*Ghu>e1Zows1r|yV8Cp?3>!C(OPo%!aN zk9}j$Nv=Wr#@-C_67`Wx1f99~@n#(IRl(=%m+xbdIO@U6|4=d}*HwEpJQuqWU!`@l zb0L*dY^kdjNMUUymSc5xmaKB5ix|?53q=f zLI1uZ6I6eFalVB_PM+2#;E;4Y`nT$bT|oTI|rm((9YCo48f0Jq2J zdOI_(hw4~%=u=Sz@jmpegz&b~oJ)Imhl_BqUs3cJ`5nzK6~r--G{vF>z6 zIJH7)TKSOP=^@%37IF{~RT63t5!E!b+-EJ^Ga9o~)Y2#P?;s+o99e^isO|8>GT^N9 zFxL_0WDpVcDZ-80;bJ@NI&BwaiT^iq=+Mle{rG>s6Zn4u{f`0tF9r0!nmbXcZ^ITB z&_Mundb!Ee%_^WbK2+FE)%{2yY7jodXH41F#=l>k8?E2;1=@IY&8=ut#*xs6 zkz+$gLq;uhIK;va%u9t&27zZ3RJ_Jzdcf`gI!5>I@zZBcaQ!GgnHv%ea-juuRsjGl z1?D#{f2&>5I8ZD2yBg;gfWz|QG9XJH0L5gK0bsDcwM7t1w-}4M+>e1h4Ub_>W4{jr zEy8%!TEH`IH~BL&?mqDxt~XF3pke*MX7QB#ZeD|uxam(dK9A+ z^SjRsO?bwRhx|bRuuTt@M?8l;M@EJ{p@6Q3X6KT}!cjn~={}#=>jOySN;Ig^AHcSF zw8K818~*}9|FGY$`$NOpC?KC)0hc=Oqt8i<5VzaoQiq}upF0%N#?C6*Q9ppxjsk}* z3baOJI`)>!=R-*)aimXk_-`(G{A0x<5n}|yJK_hD%E=>I{-nn<7S%m2pY8(AlBVmf zupgj?`pIK%7UVLAGQK_Ha=Al*Up$q+s*Hc*yt>5?T3q<>^xc_V?+S3@CZ-t=%`m>? z<$TGb1QURJx(VnvEl5aD3;4y!0zV406>V10jZ4}rdOS){O4F*g4rlagg#*Qo{))>=%V=M#KrOx!t_4UnfVT{i)3kdn^ z-s3+=-uF_o`%NF34UQ3^JNuN8CZvY=ZBlR{NdD+ zjOkY3n5fQj#(hO+S>rx)?`)Kbgqn8b1@OhDIGggY;`*4@2%Tga&j3g11Q2lmkn@-T za)7@#0`4-V-`u)&9VmYU^v6OR|MW`n5#VfqW~2XpR)9)5`d(EM% z0WC1^iSlV9)^S6d`0&~cX6W^leq1+_Y&2|NDyJHgBW~3nJRBVZ@R{o3hoh6@9_7?_ z&)eKM<`3>AJ!-_wx&MavW}YhidW!4y=2Jk-jR=T2d;((5$5`F13KTjpxuY1bkE#YP2#^w$`jfHQ22@z%79&4H@}kkE)SF8Ss6mvQGZ*^@h)vBV@`H6f6OP&RKC)Yc$<)GHO>`?%D=wYn11!Tn%G9&?7 zgxA!W|Gwe6(JC`G1FS6!q#HauJqpmw0)fOG8D-vg!^R;%v0jj0IeNgf)(i*(&VEwS zLaiO3y&*Hx0PY1R8oZ&5H>7}_wS^k{0HuhpH1q!PUcv3>-tSN23sfws-eYsYXD1c! ztoY)4MFvW>ra6k_9bUeR1W~i$lvQRX8RwH z;f33rR!8b$4Al8?!0j3CaZS}NM-gCp0dV*Gn&xE_M|`T{Wry*_v}_9;J<~d?6QpG8 zn1+sLKwNBK2Km0=+TIsm=v>SW{xaphM!D1v7==tPzVT z!L<2<@1@!BWOR7>JNy5rg zum^kbH3r=`VdGyKJQu#`Khd#eF3|K_uq(cq}@xyw*phVN69 z<;Q+>46J0Ft{K}bV~rhhv17miT*nUK!lxI(GtYl|QDf`CO>3P5n^S;D*Z{|l3XY9W zaBQ%#HAQ(eGrf0L1H*L(04MJ9xk7yEamO`44mav>KCFZRC&$3&=@auh1k6L1FsKRh zt?w95eWw7-F@p~^qp~N!RaA8=v>ErV2)ti@@rW|@Qcaq3dO2|VEFb+$&SSa&3||r} zq~>+C+%LW8PTHMRu<8^WJKIHAfpOYbk17*W7lu3&!=c&GUf~8-%}G#qVhG%MbU38& zW8w~v0oYkPclOkg_lAJ0={~2>D8lwT^d68kISVp?x(fEA0v>Y#@R%C}kGZ0-G-qi= z%NoDIUSDBJ#&_16b&a4hBcLPAe{i@x%j0B`)r`G#kG`q8NB)?#BhPwrUM zJ?xp#7Du)FNaYz;h9}#P6+fUzkw6*pT1l9{T3>qzETZj)o12v=`&MC(F-Yj-2q?yv zUwimR5ruykVd~EMIu?OT94MbEZ0f@N1^BX&#=r~%B=W@u4j6}3zPExUPGwvjRSU*O zYo)a?g1rbU04m#ffS*h^4uS0kjoYpFT1!~?jd4XgVjKr!#fR4zrkmQez$3=}0NWha zh7513pc_i?G-D5$|DOB_T52D+e#;n^;nuv{y8y^poZwE=ml_V~NBHz~1YA9gU%>dB z21fE_6$2cQyfj^O+o3>f1PKujhhzcMCAd`FGr-bna_FYF(i8_|yDhCYk?{l@6HeKpa48j+(BtZ7+)e9&m!sWY}h)7jCgyu20WY_1dRt_@wUjZy*ppob_g5D zvg06nzI$5_>*mYvRiA9#12FX0jpvKv)rtGfn3N3Krb5nl?>?{X%GYn><<2m0LIDnN zMG9UHcV(Ek5cX4pH2t)Z78?Iw*Q<|;8ykxsH*uV4ekDAG#lB)4mW8jadQUN?AY;2l zT?E9##%_$B#w(3^GO)~@J0HW{#IDSe*acNPXfE#~*B6u}W{3GZSgsm!iO#7sVD8wj zY*S!ulNj%-QSUaWcVz*{UupCV4JgsEUP~?!>4~CC7T4Qk5I@ex7Tc}t?rTvIG0;#I z(b?_VI7XVOzwXt|MZz6lcx~l=-oAWguheXUTO;Ka-R_vKl8BhZr0#*r7Lyn@uJSGh zR>#zCr?m%^v!ow84XYcK=h(BaZJgRQ7#TDGlf%r3j4L}m~rQdput zktPe<)#Cd7=M{h@2-j!R`uM@~(k@*77h*!a;~O;ZPAoN*C2ldOt%pUY}gc2ApDsvS~b5D)rf5THmP@YM>% zET5~*#hyx71wmkolAjB-$9rql7;I(_Oyt=6o8lY`z%T)2NTtKodS}A|B_xW%A`Ex1 zm)m{>xI$@W`M3(x6F^Ura9LJQV$@S^V0c((L(2$FIu10d-#o89*x0pjD^?EPM+Azo zzYexV66o^a`Nr;d2sbK6)^|@4&=o;sH+FY1q28>I?=_1%R2o_BV$_%=Vu2&-ib$wy zF_hle*}}|q2CUuv4xpRyoq0S&QZ)42YK49{SfJ=rv(q-7e~$GuW~LHv-N2Z|$XMGw zcqhgLX6Ic0qR6eDo?1m?o~%<;22gJsCwiC=#pYlpk5qDwfLPZyfN3eA!4I&&V#SE} zfBpRZ=M+B{ebF;b3e8Qe?aHOiv9A!TYvTSs`pB@ePq0Q_UISX^gGOBbo{!nHNpj`} z7Q$V6@5(D*Y-?R479U(09-DicSLLH&Wv>A3Hi`+JlJRfXb{{<7*(C-8E(y>YFiNu2 zAa_n6g@UysJQx#6ZI{Cw6XJpyo=vN*gbiFfuIAU>olisIY`y)%mjC_`M6U z_uH?t;dT!S*W3T~Si!{Wp`CLKVQV1~t?ySMYPD>?vqBH(xX}s2+6MG|AmD?y)gsC&vJrGLZqjy}8oOp=9?PV_GVR|PSx4&$ zJN6ZGJC2&p<6Q)+EVo<Auk-j20~7rsC#TF{j!2OPW)C#KUr!vWI91w-N{jRVsa~LA@p~T; z+7-Dt$!&GH=GYF>?Jd~z3S6=NVO|fm(BrV#+>J{g(0KwlmnXsGE>zvZ}Z5H2t6^^<&O2;+VVs$l(fK+wHvDYi$*Tcwj29AcQX z>EzhycS{gE#Ft`~qp>5UNYy+Y6%&f z`xEm_)r^(mf`)DJ1E?&#WRZA&*Ag}DJCZFrc|8f7h19pQweh@t97DO*dCX${l@7u0 zJz3QrWm^S0W{kPVnO`$!%l7;(B-g!!RzJlMz)q>>5V!{>*)*XgE#7H5<^C||p8YV6 zozAVFm@>A^oB6dSXD@twyO7@x<3NXxE}2@&lBspbGfatE&t>CPQZRER*&t~?tthYA z?Hpj|(Q^hkB%#=5FL`2yiOq!6SIu!Q1JY?#~|nX%2s zUec5V>?^G2V72#_$@v4codoi2k%alCrlZmgY(MW|$Jc;J!cHjaF1Dkz8#T29H0x~v zmdN9fO&x4+dVm`D?5z%BcJTPNW=fFL>b(BxW9f&xdOUAK@~Bk?0jt+ zim_!g&br^`{FzF%#+RX=}eKwO3+z z^&XE_LWfFq-#E0=vII2UH^1GC!G2R1CbeT0L8=6~B~&=uf|oaNrjG%P8zvBMI_yk#+|9e-_Sv6#ImV8V zcUZ&p)tn3~W~aRM`Raq_!A0tp+JXdUT?=K}S}5C#;z4)Hy;T@*q53#Wnm{{@>5qw7oSil=PWhxMDIYc2+ z&_zqHDB*$NGb9vqrr9A`k6<&Rs5VQ1i0NP(6C{*Gw8H~QY14j!*(vM*LQ>AuzXWgt zhak5zpWPMxN?sIoYoM!%oOY@mAwiQhMC0-72TyR6L}D&mdiiRx7sxN?DCkagBF-b~ zVfLjxD=`~iA)*ul$9OD_Vx-uyf5^GsYQgrTlA*xbZ9q$`@N;qqqPnQ3W6a=6$_Git zNup;{13Z?~{m>3Z&v-QG>?PCv}cYLMbS@93Tz+q6X8qY%P6Lk#_3WS#Z2jOZwEd z^Nqb4ODSx@F_De#P`JB0_?*Z8Zqc!fX4aYSAYE$@S#+{6bU5pXC5RW-OIruBc<8B# z5JNm3%b?0~pub)n(H`10WU+(mj6$tjc&nDHy!WyrIp-6E@OIYI1+ zO|d*vP$dy5ocd}zr)(l<;m}OqK%V|;=^I%0d+4NXDnT7hX_jKQwl|>~!ekRu404Vb zqPSc}U{i2a#OUZ5r!H4W4nXJ!AoK$y^q(0C{mXWgw)Dd6#?GUT*LMM-mme0laGo^L z(Lcn|arW%&O8cgWlJtJ_(K_ya+}uR!MzdbFRm9*4HWL7kpm5aHE`apsv(4h#>|h@B zz@>*x<`2fD?=lNmscutnE~n`aiiV%bdIh1UXV2axoO-h>lF1(9g3LoZD`Re`I>&~& z1mLD~Y+RG!hRE%LPHwlzKztT26n~#XG}|pbLNWgu&;P>zU;S@~)c?l*t`|Dtfuw*A zsMJCHzjqD@{<|q5t#gfBYN>*v5<* z4=|l7RtvdQxspnk)A?kjmd#YNS;Xv1CUWUSp<1mL^QCOAmMoVu>Ed#-kVpUsyO2fv zu2f~YP%>xPSrKh~u^p(6P^#1BGjq=U05t4B-CHSF@<89JI#1uxjvoqzXhW-AkK*RjsN@1KOH;w4@F~M3ti}~e5wwBFSvdbAWHpw7>qq76r^FsE^ zAb{f_fMYxKYp3-WI&Q$}eH(7P&Eod<=C_EeK-wy@_aXr3txf0(1elqqzuIkWxv6 zpAt>5Yu^zOA~ttrA$C<7gj`3CJD6@rybN9C3)S84vhk$vG zKO?OEXLG{sbYH{-tE}6mVazkpt5@LQK@Ls_sXVZ9J1V&e@a`+;VpriICeX3hoYfFj zOPN&zw5ziP*gC*5WKFtytv5C;a$>qU2zu~}-MdrYS2de<4kHmrU-9mk*iS+ayal6M z!kiHsA)MEX9kGAy3DI%e8SR&XjFPfta|5nO_tUzmM3_x>xWAWf*@4u{iaM|;d>(x4_d<}E=Ca5A7+-$hm6SMg;BQzD79z*(ob z!x(RNt<8Z&9MirZdtX}SZoHm2JuPn>Y|AocxGdgaTH&D%i?Pd870WD-%-ilk8Z58l zuT9xZbD(XmCNZ z4vV4&i$z(1AAg*%T|3G(%ZJ0^-GJ?486MW)pKVKe$ikZ6na-*(c@cgL|nMBS}LSR4B2NkzqHoxi{C;2N1mY#vhYo> zoP(Op{2}X29$6&N*(X@Mwh@!;|66Q+P-b1&F;G<-kFdB@;Z?h@&aWW}9?Pwr8JwUu zMHY({5~Ww#o`FL9xQ-3qB5eul8#|jjP}~>cd=KwIc-6xXyW#MIG2JLAsP)^Jk)QOP zkx(Q$X5=l@gv@`1prp9|OH4>2#$v5qg&#j9SDbhP=VA{38JcWjEQnc&-Oyb7;oKq| zXm`+mj^45=B~LOAgKEsoDqL;}I#f;*Dzu3RN9&twPsxeSbOfftV3&tUQLZl8K9+Ol z>_;){2i%voC-WUlX5rE!S!vuFjo2kOzH5j{14CNJwo39I@{h83K~MTr<|Qwpmb`+s z^y42r%LL{_+tR$)#lrqgX>=Siv?HO($&Wq?e01*bpJFwP?EOxM1g%m`XN#G1E|D+i z3#oi=8G^K8X}Oq9XEXS$W>dvl3gXq}T46bz$W}{7T3waQtXZ!KR$sKi#?f%B@1U$iL;5+Qa+Kg6fdbBdY4X}3z@KY zDqKjH(xq)VsLE}i{`=*1Zf#xcKn>oZhfe(n)n_|+@PGZQe~f)1ZBj5rI+dMm!l!!7 zx^(BPGxo|DLb> z;rixl&Hwzr##W!0K-FSL`{_pgn;q=;8v7~%WwWLIg%0>Hl?cyw2l}_4ME~~Vw}=Od zs4I_cy@TrRlp570&KF{LC0q?~(vQ<=|H1RCyYZ(2GRzVd71y^ON`FMe2Hb)`HQuRj zro>$l+zxhjYJ0pPn5TG2r?|~|&$%T$s4(>2@yD_E`w{`WwRDtf*%OwN^o-7cwWNRcDw#$_gaO9MeZC9j*yJ)d$oLc zts*2#RvyPT84}t#^Gzl^mz1DZXsg%;3yuLYtt<911e4M{TgWdV$DfPCrROokT{xfG z5?bu_4|E)McBv%4MN}cGGs{LqN4#WzD^1DVjPe_ILR6Pe%G8f zH8batglJFS?^~tFfl97fm(k&uMWy9030D*s_sHvJO{B1{7sa)8Mft7YG{dQBC$byh>0-r@kd?BB_NT~x#Vic|!l?Ug+_{Q7vS7Mnb$T4HNr>mz z7vQD#=zJndkmJwp6`T6=f+f=+k@=}N?5FaVeQl1r%VxE#HlB}NgWu5!6k&C{Gs_9_ z6bmK|R#VZqf0?8lHr2WpN%HByDv74i^h0pQj!Ih3s2}KeuUe2J1g#$C(zc&XOP_^> zn0bVX-M&Un3pW5pZ~H4G((Y8PeU;d4W#aSLk%sd6-?mo}3=h$rBGZ5!fU~LPXrUAa=@x zar-rPym6deNCrV$_j5>0cmQ9EgSM16cJ8bxRj-r1ER&7_tH@7#al$U9BS+^Uq)5tY z$nnL??Nc&jFYRi`lZ5k>ZOLWGUe7)nSZ5CMou;k7bwWbegwP$4U3<{mDV^1#CK0kQ^@YqW zOL~H{iG5F|qO!@I9|IB%@Bf?#XXlj%x2rfZ+wZ-NzN6l9)>o32PDIYT=0VP!gw+l` zjqG32(~$J+qsyx_xztf@IbNvh_PI5VqoSjSdJUROXjA8~ zu~U(Vsnl73v-vj3Nma6=o66m-K5?XH`g;6n|xde&LC;Z(#uasR^BfpL*tPPA2R0k8cTN4pA>j00eSD zrw|~KkEN1lGdU6hB+~;Sz`$Ai?7&*%XSUWDoP~jr@XN|yetx-=1PJ)WMG&dCBm-(!9W@qV5y`hgs(;?19>92fnTfU$$<8c@@ydOOUaRJp(FJwwlvrdNq-i5 zk3#FMeI;?^tev)fe;&FO`|6;1pwkrFX#wrNzc$n&yHI(4xeI*^p?~eK!k3P_=b=#T zOy{hgGbOgS8;(j)!L~A@ePptC!2QliefPfoH8kx_p=iE!2(fEG2ZxTGpGapn-MM>| zdnAh}5)+E9?9@0WLbHn4jEJMqSeI^I7&PSOR28<#n^5AeZ_2&Dw9cm((AfXOl!1ss znr`Y5YN^zE(+2pGI5pClk)DxS)7p(*!l}wSlRQIJNJ&aP*Mu5VUOl73g>b@IuOhTB z@-`42jZokajudHn96c`K<#bGyhg5o{^M!Ix|BADO^=2xDPPjxM`n2ww1n+f%u#(u# zvx;z@L9J94Y#m+8ZMgMh9xsa{NRZ=b=P(tD(y#M6GYFexF4FAK7Cx747{CQDpF*-&a8Mc`?IN6_?eHe((-qG679V)q<)xtRLcBXuY zl+Q%h_pTXr(xRyKt7}H-K52vUpxst~_Uef4qW1SNy;OEYciK);(`%)=y*u4qJh4Yn z=ZK*W>z~saG|x=CEE6lSvAMn@tPx?D1`cFKy@a)q#WSsZ#01SVYp3MIDgmxfn=v6q zS&F?+t{skrh9;nt-qB=z20kkJklhr5IJA4;zCRzkgwa;{Y6lQo;OCTGWjt}Rze zwe)f_lK@IUp_(f%=ZeWhHkHg&tK~u>pG+ld#d0l|&1aU=ZFsQ(SU@M#Tcuj7ZwPRN zv$p%|&j%LppN~5LBL1JO=<+r|L?Tm2mGX&lvQ#QnQ<+*SSIAcq%cWAfuv{&qi#Z?; zCgBL1Pi3m*Of6e0r84*nc=PT&)VtKh*DWsAK+{&jDnauL?ay6AL zF6UCkLLQitrR5x^eLj_6E*_v+0-D24cey?I9nb)6oCFiB;}wtI24caVD~8Ab+GPOk zVu3IWpj|9^$}Cwj4xn93C%d+J*xO)~185f$yJ!II0vTZ^%u`PomY*BirBs`}VDF&? z*ip=a(kZSBa4O+B^*_894JdpOe}d|7Ur~@>z0*#6+tXS%tLg?>qCsgW0+`_KZn#?7v}B=Z|=Ex8$d~* zp$W9JU6@}V1Q%)EIXKNBwm;lC?dr+?fBnJpfj9oZ8z1)c7UJbk?Tt_R`hgEVVGNL^ zi;%bteDKK@z5`YI4rnv*!N(Ah+UkK1K0;r_YW3=R#hNso)(#TVVc>%g1^&PXA4BA) z_QB7b&1RRfx$MANU=S;9V1EBIo8McG=zn2gi32lvoT8ltw z6cKD<9}tG9^frR8R9EzrE3pUd_#<|fju?6$Ho26^6xy4t{l=&R z!+#hXe)4pS;;p{{@cY2|dDvRyBUch#L~c zK%D<^)+sDDnorlOi%szMc8KEr(dmV1t{s%F+eQXfJp#@jPSE>9we`2dKsZYiZXc+3 zcVO88Gn;-^8TMwzVbO4AbKS%-dmDyow{7zo+rU`Z&}0xA61e(14T{cAt6=Qo%d0@mVB7O!4}LHv(D$H+ijFqPS^7T52?4oO(c+=#cY+jeRjWNkx^opYsj@(r2XOuQk3 za5?sCSDE0M~5y4XLz?4)* zdazot+gYcL&Fx3e+cj&MQ#_%jpx#Gi!Qh3&9K-eUZi*-0>zrI zoTHTJJf{Iu3_uZ;%}0njuw5?#6xccaci!go|E%xy&msgN;?aFuhl{cuMCfOm{}gbL z-o0qMcmKI@0h&qnbf@|M+TbE$GJU>(GQG9wZcU$5zv=U^h|_J!vg&IBTwGX+S|h>k z4avG{G8Fl1rq@MrUjIQeM*@^~%p7xm^|4^=JA>ygnb^LL!dC34pZ(V01~|ndf;`<= zD@M07pNPA2*?>70W} zuRbZ3n?FB)uF`AG{=qP*H`lA@XID12o+64vsk%NJdnmk3{dD-v{6#ZMwwkk!-56 zoJrLRwQM<=P9}5NQnp&DB+`{)ty)b(k6eWJ`*JnCTqvcg*+eo`gu`@kIZ@4ID}_uY zMIo2XK=WTJCvwGXF_i$wMY2+?6pG8`j@ z6ApZ|>$F9>koq?2YmbU^i9sd(ij}0D`C#_3#IP>UCE3zN5Ad zGWS&RcjlywZ4t2}YgZ4+n0+6h{SFzk!{64jWY}MLhvvVF_4U;MX9r0710?+clK!KC zDSrnAu;(%n=ic1>#onNGg|Qpx|2q*NO&GwRQ44M<-v6TY|Kj}5CZNbD`=cZ7C)dAf5&i!(KS@H0^AJ)blYi&?-s3>( z1L7itO||keBu;?bEF=r*d?A%drE;l4CQ+!A6Y1rAZMjq_RpH5y&6aY-RI!q&{lljd zROY*5V`G6sfBXDJ-XZu%mx}phEm5qNleKgroyw&Xg>s=#%;xjSLIDLOQfc4=E-&ZP z>2xVmDb}-+0vUfz? zL|~v-hy`G2Xi7}>*V_tpLij3pyV(LqVM{&F|GC7#{(s=gKXB#mIQ^!&A=ux*`Qr*K z1FIYLf&D)Wj|S)ejb^>L9hkgvZL|8evk9LkfK!isfW7T>-PA@eD0t-;12oS?veOmroVi=m*m>&UQ#<12#k7s zZhg9Xc*Gdd!~W?L7bjO0ub*x!epNhBrh_9cAJ-yU&>tSYb@I2$bfD|TkZR>+YO(x~ z>&da{cNBhRObJH9D=%-(ttdPkR@gD7J)2w6ri>$@4+`9em1t0-KOSB4Xor11H~s~J{$amg_lJhHQ9a~w z1zhU9k3J_cLfmeTOC5?veC|+48#}9LNBwSvA01M(sAo(ureklpd_I&^5=Z(phyUi1 z$3Ip)5-~0GVZTfF=qHc4S&+*d%J}w(%jFK~ zo}p9utIGH{&Z}Gepv8p`Pv4!{^{#j@Sf?2e%`m>?<$TGb1QXtXanl=TTF@I%v<3X) z1i%*6jTLQH(Tz*mEP6bOUS`v(CXW@Ic@?HBLC<@c{N9I`>9@SYdF>Nfj#+ddSV_D@!o#|GHKaN?zW-t-H;;f+Yms2VmRoIAv+^Vw8!wF35nQua; z6uva3OoziOt&!HW=kE0E+)C?Hg@-W3XIb#Pb$t}{F>dfvPebG81ivti`Lc?sjY%+k zu6MRk@l?khUW8uyvoH2%=EBQLfT zHpSVLhZWbyv_|M8(|86G?}Sfj%;?5rLsRiSHUjQ4rr+GUb)B&={(k1a}w5H zjI^%1bOWD$Zul6hyH%`0V+>60D8}pK>PD3PR$;rU)|$QUVsXymV5%!%SN9bKiV_{CawKlO)R*aANpTz9Xh06 zbS20)YD_aOp68l6a$dozIBlfe?iZTd-Fm#DwJzM?!{PRCrUl$?p1Kn7>qEit+~LE~ z(9mr3^vPggTA$J)<1D6kJyK?$gk3|Pi4)qAFA~v?d$?Pf{dV82A>As%KXLT-QR8yd zm<F?+;DWS3LM!{utvg@}S$Vx%$mcHLO4SuB9=4<~x;r?dN-78LMu-QEfEB zy7B3Ybu2++D&KpTX{}h?g@p1-!av-fiYY;#2DigcNX*LF3qD;m%vm)E^c6R z@fge@;$pL|OYi&r{t)+vM@}9I6WylrexjQf95*yG83>QES#Ng|CHB-JRye>zuHzmh z)~0Ctj<35QmoVqx;72S9#xka|6G4XQ7u5L$ALCO#RokdH9%6T;z2}Q!ar8yqDjE=F zr@=E-O*%t?qR%+{NaN((py{po)+J3h)Yc{D!Y*xdJL;V_PJmT=kI)xW>S5mcEymn? z;DTHDzo{tSY7^iCRUYMWKF4ER8^7haVp_G-01H0ERF$o-gQVXuPH{Ii|BFvryDt)4 zZ%w~QFuu2hZocW_t()kLhn$_qs6LLLs;btpiQ=t{pMPIc*#m{$$2OwGSXyBjkM`m{ zY~PH#Nv6VDMH?QSpZ0k+UPNB5D|}>L(WZPJrlnu5Gj_(Oj~|=B-s>>fI(Syp?KOTo z+bXDhY8C|Ca;bQmgIOc*dXKvlRKU1%+FgpC zj`@56PjJH3z(!+wBPDbxrf1v^N{xh}><2eeKWMyE!%! zUOzt?au0hvevJpXu5s{Qd|;F@)oAY3qa0mP@zsu(e1;SBZfN2XSjw=3TkwjTY_DrNGPq}up zHKnlY%Ca%_fpKE&=(X61*-J*@#YrBV9}2kr`lucrRxxg4#(d{Z%pGOuf|Z9Gv&Z$s zagVC|aY`DFx<|+C{=`iW?vAHB>7s|hP~09Z6xJq1M%_D;9X=lQ`$F$&ni~g?IbC^g z@->Ajabw8rt>kD#^SRvKK){2$ksf#N&%Id{n75S)zcb?EXTzKqMwFq}d(o(|Iphne z+Wlx02dl5}o4IID+Q3TV0}J{>QJ#U@XzW&N!56~j*{8WYuAocjE?1blM_u}u?iw9? zd42)QyKB*{jquSIp7sShYiO^?#*Oh8ag84qQ;7RH_DDFBaW^Ir*HwJuXsN5ZOZBRl z6K;%Nw}vT(SrbGGH>R5^E+=W?O=1&dL{)Flt@w>0@okJ7Z?SkDHckxr14j>^^o~EC z1*wKcqA{O-=H$h3|4V=CS_o(TL(-Y60xM33wTiGyj6-NY03rIxK!^smp8_7hp_lnt zY(LvWd!P}sod%W!RnEfD_y!bf6__D!VCs{_Zwacw4g7^N&|!gVn;y4qczc36kp7J< zkS$+6=mW{B_%w8P6nXp(z{G&dfWYAm^Ipkuw*h129xx1MR-lc|$L?>$QphQCSvL?8 z27qK1_JM#t03#F72Y3-GRE+g?jD8)o0kK4k8{&&gvZzF-0<#OO0IBu^hx}Mt6G1zj z>EDd1{GXUt!^ zyD(zgy(N5Hf9j6m~)|_-=Kb$PEVmh5ltZhqY53HK{;n`iwT_M#Y_`tq#{JO$S3jh=UpZ(v5 z@Ne+v$57xuWSqYLpZ-4L^!-nkYA}6HrL&b(wvwwAU>se{lwm2IgatG#{9!d+gzJAX zl}@Hp<$QWMmrVjIpjJ!*UA~fVIQE~lHi8Z)Ffx;NB(a4!J+`tvZoLh&>OTSB|NlAN zZS`)u<|k_TYOzp&XZ~`koXTfwmF0XfQL0o*8GtU-(v@UpIiIWGS0P`^ z&k%mBn4w6RwEJOI2NJj*61aU*0>b~#<*glCYX63<%Mx=7D2Sa~wOjIzbEj^Trcme?=)yo9JDT(9!>ud4$y7r*qZ`679=Y_tglya zW|%#D*LH%=^?Ga94*3lX&#`T;o2HS}_Lb^(69B%kn>KW1cSY=~45S1dXXJtX|Gv}J z8w^ByOV~B@?eX>t`=8qt`3noBXdS43VftIu_!*H8y6yiDy?Y4%3O@j(`=2j#pL+l4 zZX59~QOv=gI8iF+E16soKF5`1`2SXt>2kT8$t3dG{PJ=tTPVTPy#)BPY_?oUmEnw> zu2z<_K)!RFx(27-#{-Y|Gg}v&XXZ-1v9(@&y4aR=4e-Ldk@g38;RC$zbAP^g;XQn3rm86&8-N~woK8kBxSp+s*m0MreD9F>ep{xU-~d}nEM>&&7lu->07_ouMd5k zzhZu0WCBQ#phQZtCELjb$R#o&GUAIbzWwXUX6LdFo9!s$9X8uxvpYt>VY3}&{1GH8 zZ#Dk+|J2_2((*%=mVa%Fy@REDKPsEu6z{&|`hVBI_TbO?dDand;{UxBe zEa49ZGUakM*%YWgg0rkj)b=I*1B%F6toQ@aJ{>&6xBecY-6yX zr{5PB8Zepf{;rv4eAh0DjkbY*s*CX%4RH2(H7ecYF{1d z9P;t;zq$Fnf|%hgTn-rW1)$zLbU;ecL3vWgH-VBjr%F@3!HsHR$sFJJ0NdY6+==KbxS)#xNO! zM-1h*LN?#>ahvol3n-YffH#&d|NR>+S%6-M=acag@NhN{ zl$t}z(M&=wA!ntIYz=x2z=Oj9zY$S8H ze%Ys4KijRNF|AQVM)k&QacO$&CE)+94nm+Q7sn)`_S{GCMH3BJegWMa)cK-40Gwcah|&1e!T89{KJVE7J0cTDWa5ZSI!C|}nK&YopOwg@Wx4RO z$Nyvh)WiPC^(W=;*gr!C5dA$vvOmiIS4aLoG4c=ZM@fxekg31_=CwaAUjDmpjO|Dy z(&O<&;I^MjB(uddLafEBuqD+i@pw5`%%rNtVzQLhv&nKIS;=S1Fea@ObH!3IozCmI zY7VjO5~+-t`x{x@>j(?e_-bOkRz@~=sB4~g;_pNnyp_dr%WUGHSkcH4`SAbq#wd-& zx&9L;Zr=RA?9bzu|8~y0DtLv#$JXwW` zO^24Pl1S!?_vkK%9Q^sH$tx>lT}))i=xcQD3oP2GmDV?&v&4ctU?Eghn{3>f-0KODjPef9$T)E=$)7a@_Vd~ zCJ^~+Q?DR-yo|Yqs!LDkoxRrfGTL7cStAg<-&xpN-);|Gxk|T}mn%AHVh8(5U)-9#sI6bru%@JB zX77?;Ue=;OJir9Pl>XtMLi7%i}Htexj?|6b9katUaZ@@BS`!VPNu(JjME}8z-nkZvz22uRh zsz7P=^&-Ucra|D?%qM%!XYcE<4_wg_2+IpGKg!PFnNADB^^Mv}apOVDe6Te2o3WV> zKF$|Nhz6u3pPB>!)n=~N>zgaO_64%zA-|oryoQ7nh-tUfN`Pr%=O_Wtv=poX9j+b< z_PWjsKwD1p^+g;t$GX0`4j7l$VNye&ux*V$PV2R$mEz(CU~~f}I(QtT-gS;TXQ5J; zrES1#rf4ZW&~E9gb$#_efB(024K|qYCZa^2huljJY}snB`EuB2@O~YZMcqJ!-lt+; zVKOJq1DVw3+KrpG{>f$bEIn_p7wa7kKkBW;M3c9#9(u{gtdz?3LK5#tvoV##s-T^f zCV!9JRrh6vhxZLp$Ab8sUKyJKcE>Dk=ur!_+tS_VLMA>J;{ItiOf;z;ixPpxj zb`AyGC>j+lq9UaH(tvR4nz<=x3+vi-R4Q>0T3FN;Dw;$&+5#?5A%hp)CLc^{AC#MH zU+%W3(Y~K*AFM6u8+By;A!xuV!Bdor>pHTAHIFaj>G7^<%$gVU3x z&k%MIQ@*z^lKw@jv|2>80hOEhHbJDJ-6Nnq3%I(Ci>vsjN7pf# zi&$~g%w<|A$UnCZYDWV09B$t2sHJY!0#frQ3apdbTUAqcTZ9l`RzfGO?peG)OQu5DADTv3JOUq zqNByiCIoPw>&p*h*%g3ULY6nzR!m#LS@#uY-aBZWweMQ1jk-dWQC7f%mb$rwAbBvy zxw)=B(dHJcr|ru}0)7(Ee+6%M?Ez52%1Ut$QZ(r$yn$821?&PO@oAmX+uAFgrPV|` z1uf)`QmuzNACQX}ge{6=3jsX>uIx|`+FH5}Sf}KdqZ0g@BZX0dF~br>a+Wv-lSErF zUIBu*B9Q+I=$X}}rU(KvNcks`2Hm2RYIu|Y%FMP%42}dUWFI6Blb}z^;q=wAjwjLk zCPSI*+yLifx9$=E*VmResS?Q?Qr)!NB>JWk&=2^4amPC9!v*YjPnwo}=*%`On1Fq^ zZNqK8Y7L9Ue94e1^0MZ>`EMltA-egPZB{-oz2@sO^a^`M@0i3&)tiF($ASo0i-BV=>-i z6`81@BWH;{vtO5NDWq54v+=Y!!_LmMRhoM+aGSlfX*O*k0=i%usm`{lw_-3zva7R2 zsj_MYR=Bsd4fg-_Jz>BZ4mz_$=gQ&}8t>R!U7<0}86?z1Uq0 zuBFXvt13w!GH*R?3Se)8@dWV)wsX5fk3zH2(shCvwM1e~{nP6Bwq!~V)Railt*wix z-=p#n|C%1M)_Xj6y=kVcH=VWUSpwO#7Nf=Y!)omTc0@p5sx{UoT73Y+>dVN_vtxix9v(jLL_n0{WnG26&tL#G<>a6O$!0cYU(;lN+G&U z%OsHyy)bb|hw0xFRkTM%gO8JGn9^=(_-69OC%0<~>}{R1wpD8&I&ZXIZ(fc=^X5Ak zhYn2$Z3C79zy$HfTO>$HL3D0mD~6>cxAdq@@)TW>bs}{u!My9z+=gYW|4hFg?uwK7k!)UZY`ncTid zZZmCsg297et3%P1K1k6erC*XhlP)rc!Z~{>2-&L($?oal*SyK}wF=aOaTohDa5+KlB2p63LWItI7H~F;M>nEV-wkXTA>mZ!~x@4}) zyAmB8R)*IW;r)mJh$+1LGyvE6To(}QA}x&0JzU%(eKlYNvcmSsJ})u5Db5ot4Iazn zJyAtNIv`0oVDoOXXSFI_Y+F;2bb7zVtUj5`nTp`6P$-J9K(3f1tX{R2l-I$FR zl-)g|*6^ZsM4jQq4i%$KJ4)&dQu&YqY(ncax5F!y4tj32Tj~8;J5xW>)*RONXbkCm zkEYO6NYI1p5W>ruvQLVsydy>4zUhIAr>Pm7>(Dnyc{VAo+YiT30CNO(iA>RK>TG!p zx-V0?Ezh9!JEY|njJBGvHMeK(ExkGQ8o**#skhXQtv1o%I?x*dBV_xg9a0zdJlO%! zGqU(~(*H}l&nB5+SuWTn+p)96cu|b04l6K;3R;p3Q;=Yya>vfHWuq-nUWp|c#drGM~^#pmQQ6N+9rvQJnBDPJ4+%r9fG0d#NrSPCIxU*f=zt^khTkXrM73F zelZ3i#7jE9H+$Vv9N!y=vD0yD7!`2QgE|Ewq0gTFiu>Hh})$9w+QBmw|_Edv0$ z10m(l?>+iGeiTFojP{6L!_6-#{50n$)%Y42z^G#PeO2sS@bhuT&vOOVS*rL$P>9`m z#@55at1#hFxX-J6m0o8)4}Aq^ZjE@tE|*&k3!gXW@(c>E%Y92d?H6vh5Wa05i*+f$ z2ClV=F#eD$#o9(cKgswB#l3+E0Ryh8JFjO|)fmaTSY7pfn|XBMnou82Tw_AqP=hkFjZ}3h!rpfXaJ}xknMfvz#aWV){3X&u~8E_Uw!r^9p~B31i45 zzNpppVc*W_*qI5J%Fo0EJ|>vh9T)rtcgY}zqv|8yM?y8ue8j~XBCz!HD&xAM*3re! zbv5+e#BR{I-!B3WHMo2Yhx+2@{k`WxOy&J?=HVCoy-Y+y?!cv>kuX#) zjs+R-W4vGSu9b?V8b!wu+|Rg|g?z84bGz@{pkwn24>0=k`aQ-*=Cce z^p5rJ7Vd<(hmZ7mnL0RtR~PVoHPs-zgc(0zX0F5+a2E zsw*H|D*XyY-f+3igFz~RdqTLIm;6+cuHXm9UL5rb*RRBfE-^yaPU2oPtO}J*t_by1 zP@ufPRr9lX`c7< z1y^hON^zCyC<^fiRmun=qaLHY^Cz}ZwZgf<8muBy*O7eyS;P?bIdas55ZJ_IK+62RieqdnmQvq0enY9%B`yX{A_%GIp(`fCE3YWJIgBYO~>r?Vulxn#iNfg6dOW8pao#! zj|mfJI8XS+6=33N&If^sHNwO|0>+@?Zq({kc$triM_2LSyjwhg~~e_^Bx@zsOZ4`r|_MPchH#pm9~GalnI1hk0=YjRQ^ve}(fC zLJ&&n7dvZ&c$DA@1#cq<{~{c7(N+>|3xaR1qi1vIS)KdZ zKE%Bo1BxI4%f*f52dems>v;D*MS;7&+P8D+)9-}gHT)v5`)m5<;)_p#{1S>WnIu7;bS3xbNlegT?3!kYVhat2i(_RU(4KP+#@#@H{-uz{G#9& znR=0S5)?hBqI&g>og19v_U{;9XWXqQ*q*}TQPL%hBzGgKQw1S@hdsvfRqRGr8BYkF zU~0&t|D4oub3E2}Zc+pJ$?P zUU;z1rN4HV1Feh-L9si+Vz<}Ddlg>3FVvA8Z1~5R_!B1m+i%psx%a=`QQpNC$VC2Y z=J{&Sc!hZ$(5I-s8xR9`j|F_#FsPbqG868OMYmhJYPQQe6? zuhJJ%usZ>FxZR+b+(SSWA?zFx^%dGbgd4k=U&ZDV^3n$6^Lis*U&P~4dxhHTi@5#A z+`+JLA6Io}f1ZH0SpF-4Rvw ziT*&KFK|M5gKEF*#4$zq{k}1L`hviOa-Tdt=m{w9bE@!oRd>MSj(EZYOpU5;CbT{l zR;<@y#&6XhyBEjMVrvYT(W5-$kOIMkPART?enrJT>7mVA6)vCqGDud8-*LZ!FS_G1 z?-*0+IF$`yYy?of8^11KO@ZOYM(+YK;qn9ftN8V}Kv2*H58ZZ!Lf8`o^90cz`{34L zTV$TzlYJWVgy_@Y#<^$g0KWln^s4=5R9Ce3lsY23f!^>b;f)C2IDL8tT!ICBvJ6vC zVjJ-&D8~qBP?7~e)T1l|EeO)a%lqWHe98!cVpo-*FGvIC!yCFhLmqF);~MY>E+`=% ziUxwJ;&pq^fkp6!RkuqCV>E|^uYU-ekKdifC69}#$FZ8?$^lielUhe*?27i->bRnB ztbo#ts{CuNAo8(LPY~WX&OCptoIfW#-=@x=Q~P3sJE9XiKnmxWD16O)ccX+kd~@bk zgYE@Ds$x!~M>!VDr*f{Sq&a%LO`K|8eynG-M>^~VCUKq^X zW@=xspQ*m!FxUocg4jLNkHLdfszMNN4l-}7aFU0=fA!7M&I#jJB5Z7!M_yIIyosIOpvU;FB=%-3GO03nH8fn#Ik*jQi1#yW#7Sq*>E_x40*cISeF!9{<6 z|EaUDL;}Y%Z+;lb#`|{F3&GwYVDiD>&e~usa_UTfXe1Ko?VIc!4)@;~2)-H#9}5Sk z&wYP0@PRREoX;h47eWKU;J@@Rw%o&(ga4wu3*$Z+_w48G^YEff82`&}#wdBf+ry)y zx!fNUyMsui5Q#kcOeRuD=Tq54wF1Y>O0`rjX5)G}jRXnte7Tw}<&#Ldol8~qM6y~a zB=thBfYAN%Tqc*-3zrsRAtwL%Y0*a8od`(dJ8=!RacAq)_x+OF7CZ5X0)5&UD2Wt_}XHv)&rdKk_IHFG?e6N|7tTX=*F%?O9SH_il zE)0T`UkOIxO{K1e@D^nq;g7C8A%juVp!^^kREnS^qYc^Wfu8>w>Hihw-5&qHg~IrH z-7)k&Omd@NjG^~)S<%0s$bbLin=__}=fsJ@!GAvS*Q3C~R-#5vBvXl}ikWh`009s3 zL#EQHLbg=Ml}iO6<7%>)Nf%4`axnv3oXr)1vh&4kp^#2xQi)QslFVARNC!qSBy60v z3RxyfWRpby&ODC2IySRhEa}S>f=fntnn4^cnPy803@r5c%gwvvL~uMxAX_JbBZPQP z1V@vmirM2%+GAZh5ghI5^-3~KsuRHxLC{)J7@sDf9*rOc*=*LcFM=az`kVjl;;Y9U ze(=GB-1EntCbi!$F5fQkga3FrdIDf$C40aAN8L>JzI#_HDGQ}!I#nzpmv}OhOD0Oi zcqyJhCi8qHnXF{8gg3h?MEDOEB#pB=uflG97nq*D12z)*xL5>nVeL8YL;CR83K z;GCXaSr6N?WNZpXwO!qPyZhdI&Ue0Z?{{?13pMw=z@JJ)VIU9)wBg^GpL}|Kst50G z`0%fkP(wfwf&t|zMGLjdr{$MrS^9+3CO#!@73Rc<5L~KBmieWiP)kW5u#^bJf!swF z@8dTvT1xlg_xQDm`8~kL<8fFh&zRPzJv&~SnXVql(lbh)?CuR>$4vdc-+iO<1$*|s@N^VMN zo07>5xlDfPzXuM&5Dq5BJQb=sQ;v%R%g|l z)Xi##8Vr9m{2xm-dp;K~H*tfIN~W>{$@D-f+n!8r(sP?K=?yteli##0txIp(mduH7 z-lhu|v<;l?Fqz^nXnd_LYh;<<=6$c9#1~Js<$}b&T>Qqmc#vB?t5>pFs^E-4c|4#z zjMT02lm7K)E{j5&lQqN$>MrWPv_?18wlWx((t>2Y&D^(*CKv=Q!~T z8&3=mo~pV4otQfJ(dQvPe*^Jd)vMOmr_)M$#Z@PI7Zc4Jm$Lm)LKhNT+}m(9#Fl=# zmZKC$v~!Vf-u{}?7oYAw(Rf^c1BWljg1Bt$h9`(bxS;?yBtIV5ar$(w>7ymy$Jl{* zN5_sG_jAY1!BaVZd^o)2d-lQ(SiJU!ElZz$u^S%#1yJv~eQ@J>GaM0`>n<|Ir7!(t zn1c9ZJBdVm=n(Z~JO@Q>Ykxw+t--AD(foFAu7%aBk(-VkAKbDb9L<$!g?8=%u2$_;dN zbvN~FyWWkuJ2U%fsHS=YNUL`uP1cn=-K(^#ME{XWLc}MGHX4dnZvg2gaEP&mM2x@; zC0l@8LEuij(HH|rQ#!IA$aV-<4$d)mdov|>0_c=hW=N<~p6_;RD0LG68Aj)};pqVY zlSaf)Q>?|^ydNaqKZEjBLTyZ>vMADcOr~`;c-9w|L($&Ot_|@@r#sxQ4s=0?uRaBW z{~W&SKM7Ik%yS?}_!^bXx zBph7`d5`4D!1@MvWbn$a?z+l*&c9WR;Vl#&4;{RGFvDy;^Ixrg2HRBUUv+q6AMxe#d1K>sqARnn(!0NT7&s0SVepY20Sg@iBLJtB=gk+*&q& zlu_mmF(Q{7$DYD*l*~~U&WCC`6D~#pmknXq)N}$tw{-`Dv(tzOLI{2iXyqs}iq6rg zvR#;A^yFS5OpKMr%^BP8aUX#x)=)Xt6)WwQm4e~&RB-FMFc_wpbVAR%+X>7vIzI#G zqQ1idFY%%rm-#TS@M?(X>w6lBQ0be~DA_?cvv8OSTDf4DiGB-jQac`6S*z3lJaJ=o53=)dAJjMif*lMM8ORMj;;GL!E z86Qwx01gce#Y8{E*w&Z`ZmK=RZ;U3ycwAc(sl_-EFpW#H%!$LbVdnot!TWlJWXdwI zwh#}r4N+Aq#-Y9ag6DxoHNuxhXVB0^N+cixNa8j?+bG`?z8{_5}ON#!-=_>=#`3! zLE4TN0EZc^TuF&V=)`Ezx7k51c+^IlQbfQt-k+Q#O%@^6BwnGb9b~Q&p|=(Z>EK8U zVPhEM_B1xe31GZc(g>k2#}(~68L0&lb()N^!Gd$c*qeJB5p8pNcFJG)Q9t}*P}jUADkkdG#I@+>CJ?wpR0I~QQ-QAW z?~vexZy-uJq_Opc&W_HM7(sJjf*WVeCvR41Q0B@^k$8cx!IAn=y9daG!x74bIb%+# zN7p0L7>E7iJLGyoVuYXcO*k>d!p0e5pgu7$pBNf}tR-}+XiqUV=%*4p+=M#XXE`;SUgLYQi$D(r9Vs z7;+>K5hQ6POfZjEuL&|+h|+>RJxxQ6EeUt904&~;fT){Lwc$EvKElyn6;*IJ?&L6`%hyN?o9vd7GEKL^(XCci88QAK7VWlUvM z4vWm|&Coi&GZYCAHSYCxWaw4pG^4=eNEk?uPb#cGh}j6E9Q`s(xyeyHso(@N*nMGQ zyogsEC9Om#jg>~tsnKId&PO@UXVcK}L05ET0)Bo+tW27iihLhPAjvD9K{l&i`^jA) z!c57|jn2H!5A_o+s2(;k zaGU&|8bLVO(6^4mTodiAMZB~aOaazjn}CFS!!?StYp1nC*t2~XWqAc5lcGpZLWjF2 z3Any#(7LZ4yv8*uo81jtf#5xKTOW`d5RKA_DSM0xl}@AFOd@mpc3Q{`$k(`Np_Mg^ zI@pf%$ujN$5svS7F(Xf~(1GaO1TJl8Ym|_q6X=r7fa$Q#o#pw;kzrQZLRp08i7;+W z%#Ig*`a1a@jxXWI#?Uk2%I&Ms@_sHjjy^30(6VJZjckK-uWko8f*i6Vi0L3`j?Vh| zrXG-CsT=3WTHue`76zUjz~u?9=y0HS5PBH%68fCuxDtAZR&>IZW)8DxthZGbUG4ZN z(E46r&9X>H8sJVE&bYbYHj1>D2DiPR=d9YO*Y4Q&Qnjtrg>W6+?9#RJ8U|b-DxYGKztL@7*`nZC zOZ6JmR|4uQ%QX>rN$m_= zja1Thtc;Z^<}}03ndy{O&@|h!EGLN@ecQAO`Bc`?%yhvrb*GrgnyFMKYnhpBv8biB zPk!Lo`8&*I++p^>;-hyzxjKMePwk*=aZ~y47r!(0S!@{bN++NjTlagPyAhP-drv~^ z%I5To3-CVgz7}r8`@cnqb9#G^_u@vV^ZFiyi2wX@_e1cQ|7h6z!NV|9-O&EA2b10- zK^pR&eF5@b+aKYJ8GP}}*FnPj`TWX@-tJF9wEFn;r(h2E-GPR`e2FH#nn%dbm&v2# zH?Vx91AZmcoUe4CxUf8Q4g3bY(XEHYnc40@Qf_ zX~AFL`~^Yscl6hxYBV(T&-)%Io!!xAin)SWKu_s*Rx3D$leW_Nw3f0oO)Kg~(RTEL zleP6UA_DIqY^;=(Oy<(*oRfAe?~XfQ9b~*mT^Ks>58I^1*Sh(^fd^)>n^i0ddd_lm zJ8zmv2UR6KZ|gbJ!B1Lt#<2BbCYg1Lc1}x~x{sDZ-ZHa}X=GA)!?coV?VNrdc?)Os z9=M2rkxv)wtWor5>1ervlQQf=KA&}P#yG5;ZX`3hre|=ZhHV?ELcvIV((dg(1meRVC40S7zeJb+=n$Or z-52wAohCOePc9N%pXlCW50bC@KwF-Bh&=CmJGlJA(N%mykM|9aB$l7~82Ki6r*=Z^ zav}-e<(9v_37+CsvGZzxy^7e`Q-D_Q*zHjAZZlyQekbgmL5+UP-wvJQx_>EE5pu?L zmSYt2>Aa=o^nB7XQ%2S@@=jJyo3@cm9d@Rxz8?b4D>|`=BvwWSop+rgbm* z40L+Qd7vvmxu^r@9wV>(j({|B%f$ejPiOMz`J#i9$=goe$>|xRkV>U<22QbHYgS&* zBDo=mvX-6F^<3J_XK~Vop*y*J{<3C`9CFH6RTHN_C^Rg#Ovblq*cGz3ma*~s`Ftjw zG@PQ7)U;#{FAdYw&;vQsvXfcEbo5j*nJwfChOQY!w7j14zHQtE_YhbhV6hwV}#r%Ug>uu+sTGHaLES2wvqB`Q&R6a(@-xggzAD zzd|4S{ONn?E&ted4vK%=cpd>|is@|DMt2pmR!TD*Ep6g0$0!ui*D*=Q5&r0&YWxHl5le`qaA>HGd^=@6fz`mU?v)w?uk2 zTTG_ybk26LuVOl1EMmgWjhK-9Gp=JbWvW{Iqi(9RIq_$YvjC-tQaYmXVW1q}WxWF4a}@a?l6_jo^&W==c}&o_+y#zSV-k z2gfUK7FHgHRZsb$5!p|q2IlKb=EGyh@2-MoUhyS(!_4Q9n{t>U9S!qo&Pb*+j%{Zx zr;v1t2;y8Sha^?BipdnBDuao)X!!J9&{D-hapth??{2D#x)mW~pkz&I8&Vte)uz2{ z<&l+@(w{BD#3ikpmT}$mv)92+Ti$jJyn1%I6CvLJJ_?`lo_&^VT70@!^`5^DLVu1e5zL6u-|EvMFBP$&=%X*DJEP>sbKIeN(HfM zt>AT~0+{GQc@ifY2LO0b0+i0=VRC9pLJBrVxNC;@hkPy!UJ!fSuWD zFYJ8#HG##N`)bRVTm^W$MF0_%Lb|XAP~5lb+NGcxso`-K%p{pA9Ot!RRp8tbLEqvB zKP;3lw>W?b&~H~9xYgHoUaJk1FQYgx@{Sb;&ejGFyo<#F@0kj@t^7Nd{4Z4e%fGec z-+ET^XG*VF`_{5Q$}-mgl<3-=E_jchB(+!`x(xMST8J2Y=pqFGp;`by4IJUVOCWDp z31D5-N`SZiUeZwi=E8r-7yk3B!oRopA<|iXR|Nh`*Z0fssKDR8D)2KUxFAlR^)p(cAob9H+@B*cVB-Nm(gdL(-#{;!L(g8c?5Qq}82uL<_N zuh^I0M6z#MmF&YP*^6N%r>xOeHT(KE*6fvYHG8G1+1J+y_2JlxW*eF=x7EgIb+ 0$, $g \in C^k((-\eps, \eps) \times U \times V; F)$ ($k \ge 2$), and $(x_0, y_0) \in U \times V$, then there exists $\delta > 0$, $U_0 \in \cn_E(x_0)$, $V_0 \in \cn_F(x_0)$ and a unique $C^k$-mapping + \[ + \beta: (-\delta, \delta) \times U_0 \times V_0 \to V + \] + + such that + \begin{enumerate} + \item For all $(x, y) \in U_0 \times V_0$, $\beta(0, x, y) = y$. + \item For all $(t, x, y) \in (-\delta, \delta) \times U_0 \times V_0$, + \[ + \frac{d}{dt}\beta(t, x, y) = g(t, x, \beta(t, x, y)) + \] + \item For each fixed $y \in V_0$, let $\beta(t, x) = \beta(t, x, y)$, then for each $(t, x) \in (-\delta, \delta) \times U_0$, + \[ + (\partial_t\partial_x\beta)(t, x) = (\partial_xg)(t, x, \beta(t, x)) + (\partial_yg)(t, x, \beta(t, x)) \circ (\partial_x\beta)(t, x) + \] + \end{enumerate} +\end{lemma} +\begin{proof}[Proof, {{\cite[Proposition VI.2.1]{Lang}}}.] + Let + \[ + G: (-\eps, \eps) \times U \times V \to E \times F \quad (t, x, y) \mapsto (0, g(t, x, y)) + \] + + then by the existence and uniqueness of ODEs, there exists $\delta > 0$, $U_0 \in \cn_E(x_0)$, $V_0 \in \cn_F(y_0)$, and a unique $C^k$-mapping + \[ + B: (-\delta, \delta) \times U_0 \times V_0 \to E \times F + \] + + such that + \begin{enumerate} + \item[(a)] For each $(x, y) \in U_0 \times V_0$, $B(0, x, y) = (x, y)$. + \item[(b)] For each $(t, x, y) \in (-\delta, \delta) \times U_0 \times V_0$, + \[ + \frac{d}{dt}B(t, x, y) = (0, g(t, B(t, x, y))) + \] + \end{enumerate} + + Let $\beta(t, x, y) = \pi_2(t, x, y)$, then + \begin{enumerate} + \item For each $(x, y) \in U_0 \times V_0$, $\beta(0, x, y) = \pi_2(x, y) = y$. + \item For each $(t, x, y) \in (-\delta, \delta) \times U_0 \times V_0$, + \[ + \frac{d}{dt}\beta(t, x, y) = \pi_2(g(t, B(t, x, y))) = g(t, x, \beta(t, x, y)) + \] + \item For each fixed $y \in V_0$, let $\beta(t, x) = \beta(t, x, y)$, then for each $(t, x) \in (-\delta, \delta) \times U_0$, by the chain rule, + \[ + (\partial_t\partial_x\beta)(t, x) = (\partial_xg)(t, x, \beta(t, x)) + (\partial_yg)(t, x, \beta(t, x)) \circ (\partial_x\beta)(t, x) + \] + + \end{enumerate} +\end{proof} + +\begin{theorem} +\label{theorem:local-frobenius} + Let $E, F$ be Banach spaces, $U \subset E$ and $V \subset F$ be open subsets, and $f \in C^k(U \times V; L(E; F))$ ($k \ge 2$). + + For any $\xi, \eta \in C^k(U \times V; E)$, denote + \[ + \Xi(x, y) = (\xi(x, y), f(x, y) \cdot \xi(x, y)) \quad H(x, y) = (\eta(x, y), f(x, y) \cdot \eta(x, y)) + \] + + If for every $\xi, \eta \in C^k(U \times V; E)$ and $(x, y) \in U \times V$, + \[ + [Df(x, y) \cdot \Xi(x, y)] \cdot \eta(x, y) = [Df(x, y) \cdot H(x, y)] \cdot \xi(x, y) + \] + + then for every $(x_0, y_0) \in U \times V$, there exists $U_0 \in \cn_E(x_0)$, $V_0 \in \cn_F(y_0)$, and a unique $\alpha \in C^k(U_0 \times V_0; V)$ such that + \begin{enumerate} + \item For every $y \in V_0$, $\alpha(x_0, y) = y$. + \item For every $(x, y) \in U_0 \times V_0$, + \[ + (\partial_x \alpha)(x, y) = f(x, \alpha(x, y)) + \] + \end{enumerate} +\end{theorem} +\begin{proof} + Using translation, assume without loss of generality that $x_0 = 0$ and $y_0 = 0$. Let $B \in \cn_E(0)$ and + \[ + g: (-\eps, \eps) \times B \times V \to F \quad (t, z, y) \mapsto f(tz, y) \cdot z + \] + + then by \autoref{lemma:split-ode}, there exists $\delta > 0$, $B_0 \in \cn_E(0)$, and $\beta: (-\delta, \delta) \times B_0 \times V_0$ such that + \begin{enumerate} + \item[(a)] For each $(z, y) \in B_0 \times V_0$, $\beta(0, z, y) = y$. + \item[(b)] For each $(t, z, y) \in (-\delta, \delta) \times B_0 \times V_0$, + \[ + \frac{d}{dt}\beta(t, z, y) = f(tz, \beta(t, z, y)) \cdot z + \] + + \item[(c)] For each fixed $y \in V_0$, if $\beta(t, z) = \beta(t, z, y)$, then for each $(t, x) \in (-\delta, \delta) \times B_0$ and $h \in E$, + \begin{align*} + (\partial_t\partial_z\beta)(t, z)\cdot h &= t(\partial_xf)(tz, \beta(t, z))\cdot h \cdot z + + f(tz, \beta(t, z))\cdot h \\ + &+ (\partial_yf)(tz, \beta(t, z)) \circ (\partial_z \beta)(t, z)\cdot h \cdot z + \end{align*} + \end{enumerate} + + Following (c), let $k(t, z) = (\partial_z\beta)(t, z)\cdot h - tf(tz, \beta(t, z))\cdot h$, then + \begin{align*} + \frac{d}{dt}k(t) &= (\partial_t\partial_z\beta)(t, z)\cdot h - f(tz, \beta)\cdot h \\ + &- t [(\partial_xf)(tz, \beta)]\cdot z \cdot h - t[(\partial_y f)(tz, \beta)]f(tz, \beta) \cdot z\cdot h + \end{align*} + + where by assumption, + \begin{align*} + &[(\partial_xf)(tz, \beta)]\cdot z \cdot h - [(\partial_y f)(tz, \beta)]f(tz, \beta) \cdot z \cdot h \\ + &= Df(tz, \beta) \cdot (z, f(tz, \beta) \cdot z) \cdot h \\ + &= Df(tz, \beta) \cdot (h, f(tz, \beta) \cdot h) \cdot z \\ + &= [(\partial_xf)(tz, \beta)]\cdot h \cdot z - [(\partial_y f)(tz, \beta)]f(tz, \beta) \cdot h\cdot z + \end{align*} + + Therefore by (c), + \begin{align*} + \frac{d}{dt}k(t) &= (\partial_yf)(tz, \beta) \circ (\partial_z \beta)(t, z)\cdot h \cdot z - t[(\partial_y f)(tz, \beta)]f(tz, \beta) \cdot h\cdot z \\ + &= (\partial_y f)(tz, \beta)[\partial_z \beta(t, z) - tf(tz, \beta)] \cdot h \cdot z \\ + &= (\partial_y f)(tz, \beta) \cdot k(t) \cdot z + \end{align*} + + Since $0$ is a solution to the above equation, $k(t) = 0$ by the uniqueness of solutions to ODEs. Hence for every $t \in (-\delta, \delta)$ and $z \in B_0$, + \[ + \partial_z \beta(t, z) = tf(tz, \beta(t, z)) + \] + + By adjusting $\delta$ and $B_0$, assume without loss of generality that $\delta > 1$. In which case, if $\alpha(x) = \beta(1, x, y)$, then + \[ + \partial_x\alpha(x, y) = \partial_z\beta(1, x, y) = f(z, \beta(1, x, y)) = f(x, \alpha(x, y)) + \] +\end{proof} + + + diff --git a/src/diffusion/martingale.tex b/src/diffusion/martingale.tex index 3ddda8f..9a81676 100644 --- a/src/diffusion/martingale.tex +++ b/src/diffusion/martingale.tex @@ -30,6 +30,67 @@ is a progressively measurable process. \end{lemma} + +\begin{proposition} +\label{proposition:martingale-quadratic} + Let $(\Omega, \bracs{\cf_t|t \ge 0}, \bp)$ be a filtered probability space, + \[ + L_t(\omega)u = \frac{1}{2}\dpn{A_t(\omega), u}{\real^{d \times d}} + \dpn{B_t(\omega), u}{\real^d} + \] + + be a second-order $\bracs{\mathcal{F}_t}$-progressively measurable random differential operator, $\bracs{X_t|t \ge 0}$ be a $\bracs{\mathcal{F}_t}$-progressively measurable process with continuous sample paths, and $f \in C_b^{1, 2}([0, \infty) \times \real^d)$. If the processes + \[ + Y_t = f(t, X_t) - \int_0^t [(\partial_r + L_r)f](r, X_r)dr + \] + and + \[ + W_t = f(t, X_t)^2 - \int_0^t [(\partial_r + L_r)(f^2)](r, X_r)dr + \] + + are martingales, then + \[ + [X]_t = \frac{1}{2}\int_0^t \dpn{Df, A(r)(Df)}{\real^d}(r, X_r)dr + \] + +\end{proposition} +\begin{proof} + Firstly, since $f$ admits bounded derivatives, the integral term is of locally bounded variation. Therefore if $Z_t = f(t, X_t)$, then $[X]_t = \angles{Z}_t$. + + Let $\bracs{t_j}_0^n \subset [0, t]$ be a partition, then + \[ + Z_t^2 - Z_0^2 = \sum_{j = 1}^{n}Z_{t_j}^2 - Z_{t_{j-1}}^2 = \sum_{j = 1}^n 2Z_{t_{j-1}}(Z_{t_j} - Z_{t_{j-1}}) + \sum_{j = 1}^n (Z_{t_j} - Z_{t_{j-1}})^2 + \] + + As the above holds for every partition, + \[ + Z_t^2 - Z_0^2 = 2\int_0^t Z_rZ(dr) + \angles{Z}_t + \] + + Since $Z_t = Y_t + \int_0^t(\partial_r + L_r)fdr$, + \begin{align*} + \int_0^t Z_rZ(dr) &= \int_0^t Z_r Y(dr) + \int_0^t f(r, X_r)(\partial_r + L_r)f(r, X_r)dr \\ + Z_t^2 - Z_0^2 &= 2\int_0^t Z_r Y(dr) + 2\int_0^t f(r, X_r)(\partial_r + L_r)f(r, X_r)dr + \angles{Z}_t + \end{align*} + + On the other hand, + \[ + Z_t^2 - Z_0^2 = W_t - W_0 + \int_0^t [(\partial_r + L_r)(f^2)](r, X_r)dr + \] + + Therefore + \begin{align*} + \angles{Z}_t &= W_t - W_0 + \int_0^t [(\partial_r + L_r)(f^2)](r, X_r)dr \\ + &- 2\int_0^t Z_r Y(dr) - 2\int_0^t f(r, X_r)(\partial_r + L_r)f(r, X_r)dr \\ + &= W_t - W_0 - 2 \int_0^t f(r, X_r) Y(dr) \\ + &+ \int_0^t [L_r(f^2)(r, X_r) - 2f(r, X_r)(L_rf)(r, X_r)]dr \\ + &= W_t - W_0 - 2\int_0^t f(r, X_r)Y(dr) + \frac{1}{2}\int_0^t \dpn{Df, A(r)(Df)}{\real^d}(r, X_r)dr + \end{align*} + + where the process $t \mapsto W_t - W_0 - 2\int_0^t f(r, X_r)Y(dr)$ is of finite variation, so it must be $0$ by \cite[Lemma 5.13]{Baudoin}. + +\end{proof} + + \begin{theorem}[Integration by Parts] \label{theorem:martingale-ibp} Let $(\Omega, \bracs{\cf_t|t \ge 0}, \bp)$ be a filtered probability space, $\bracs{X_t}$ be a $\bracs{\mathcal{F}_t}$-martingale, and $\phi: [0, \infty) \times \Omega \to \complex$ be a continuous, progressively measurable function. If: @@ -81,7 +142,7 @@ \begin{align*} Y_t^{x, g} &= \exp\braks{\dpn{x, X_t - X_0 - \int_0^t b(r)dr}{\real^d} + g(t, X_t)} \\ &\cdot \exp\braks{-\frac{1}{2}\int_0^t \dpn{x + Dg, A(r)(x + Dg)}{\real^d}(r, X_r)dr} \\ - &\cdot \exp\braks{\int_0^t [(\partial_r + L_r)g](r, X_r)}dr + &\cdot \exp\braks{-\int_0^t [(\partial_r + L_r)g](r, X_r)}dr \end{align*} is a $\bracs{\mathcal{F}_t}$-martingale. @@ -179,15 +240,7 @@ (Y_{t \wedge \tau_n}^{x, g})^2 \le Y_{t \wedge \tau_n}^{2x, 2g} \cdot \exp\braks{\int_0^{t \wedge \tau_n} \dpn{x + Dg, A_s(x + Dg) }{\real^d}(s, X_s)ds} \le CY_{t \wedge \tau_n}^{2x, 2g} \] - where $\ev\braks{Y_{t \wedge \tau_n}^{2x, 2g}} = f(s, X_s)^2 \le \exp(2\norm{g}_u)$. Therefore $\bracsn{Y_{t \wedge \tau_n}^{x, g}}$ is bounded in $L^2$, uniformly integrable in $L^1$, and converges to $Y_{t}^{x, g}$ in $L^1$. - - - - - - - - + where $\ev\braks{Y_{t \wedge \tau_n}^{2x, 2g}} = f(s, X_s)^2 \le \exp(2\norm{g}_u)$. Therefore $\bracsn{Y_{t \wedge \tau_n}^{x, g}}$ is bounded in $L^2$, uniformly integrable in $L^1$, and converges to $Y_{t}^{x, g}$ in $L^1$. \end{proof} diff --git a/src/index.tex b/src/index.tex index c4c202c..f80ff8d 100644 --- a/src/index.tex +++ b/src/index.tex @@ -1,6 +1,6 @@ \part{Diffusion Processes} \label{part:diffusion} - - -\input{./diffusion/index.tex} \ No newline at end of file +\input{./diffusion/index.tex} +\input{./calculus/index.tex} +\input{./sde/index.tex} \ No newline at end of file diff --git a/src/sde/exact.tex b/src/sde/exact.tex new file mode 100644 index 0000000..d9c1b69 --- /dev/null +++ b/src/sde/exact.tex @@ -0,0 +1,103 @@ +\section{Itô Existence and Uniqueness} +\label{section:exact} + +\begin{definition}[Lipschitz Coefficient] +\label{definition:lipschitz-coefficient} + Let $\sigma: [0, \infty) \times C([0, \infty); \real^d) \to \real^n$, then $\sigma$ is \textbf{Lipschitz} if there exists $C \ge 0$ such that for any $\theta, \eta \in C([0, \infty); \real^d)$, + \[ + \norm{\sigma(t, \theta) - \sigma(t, \eta)}_{\real^n} \le C\norm{\theta - \eta}_{u, [0, t]} + \] +\end{definition} + +\begin{lemma} +\label{lemma:exact-uniform-bound} + Let $(\Omega, \bracs{\cf_t|t \ge 0}, \bp)$ be a filtered probability space, $B$ be a $\bracs{\mathcal{F}_t}$-adapted standard Brownian motion, $\sigma$ be a $\bracs{\mathcal{F}_t}$-previsible $L(\real^d; \real^n)$-valued process, $b$ be a $\bracs{\mathcal{F}_t}$-previsible $\real^n$-valued process, $\xi \in L^p(\Omega, \cf_0; \real^n)$, and + \[ + X_t = \xi + \int_0^t \sigma_s dB_s + \int_0^t b_s ds + \] + + then for any $T > 0$ and $p \ge 2$, there exists $C_{T, p, n} \ge 0$ such that for all $0 \le t \le T$, + \[ + \ev\braks{\norm{X}_{u, [0, t]}^p} \le C_{T, p, n}\braks{\ev(\norm{\xi}_{\real^n}^p) + + \ev\paren{\int_0^t \norm{\sigma_s}_{\real^n}^p + \norm{b_s}_{\real^n}^p ds }} + \] +\end{lemma} +\begin{proof}[Proof, \cite[Theorem 11.5]{Rogers}. ] + Assume without loss of generality that $\xi = 0$, then + \[ + \norm{X}_{u, [0, t]}^p \le C_p\paren{\sup_{0 \le s \le t}\int_0^s \sigma_rdB_r}^p + \paren{\int_0^t \norm{b_s}_{\real^n}ds}^p + \] + + where by Jensen's inequality, + \[ + \paren{\int_0^t \norm{b_s}ds}^p \le C_{t, p} \int_0^t \norm{b_s}_{\real^n}^p ds + \] + + and by the BDG inequality and Jensen's inequality, + \[ + \ev\braks{\paren{\sup_{0 \le s \le t}\int_0^s \sigma_rdB_r}^p} \le C_p \ev\braks{\paren{\int_0^t \norm{\sigma_s}_{\real^d}^2 ds}^{p/2}} \le C_p \int_0^t \norm{\sigma_s}_{\real^d}^p ds + \] +\end{proof} + +\begin{lemma} +\label{lemma:lipschitz-picard} + Let $\sigma: [0, \infty) \times C([0, \infty); \real^d) \to L(\real^d; \real^n)$ and $b: [0, \infty) \times C([0, \infty); \real^d) \to \real^n$ be previsible path functionals satisfying the \hyperref[Lipschitz condition]{definition:lipschitz-coefficient} with constant $K$. + + Let $(\Omega, \bracs{\cf_t|t \ge 0}, \bp)$ be a filtered probability space and $B$ be a $\bracs{\mathcal{F}_t}$-adapted standard Brownian motion. For any $\bracs{\mathcal{F}_t}$-adapted process $X: \Omega \to C([0, \infty); \real^d)$ with continuous sample paths and $\xi \in L^p(\Omega, \cf_0; \real^n)$, let + \[ + P(X, \xi)_t = \xi + \int_0^t \sigma(s, X)dB_s + \int_0^t b(s, X)ds + \] + + then for any $\bracs{\mathcal{F}_t}$-adapted process $Y: \Omega \to C([0, \infty); \real^d)$, $\eta \in L^p(\Omega, \cf_0; \real^n)$, $T > 0$, and $0 \le t \le T$, + \begin{align*} + &\ev\braks{\norm{P(X, \xi) - P(Y, \eta)}_{u, [0, t]}^p} \\ + &\le C_{K, n, T, p}\braks{\norm{\xi - \eta}_{L^p(\Omega; \real^n)}^p + \ev\paren{\int_0^t \norm{X - Y}_{u, [0, s]}^p ds}} + \end{align*} + +\end{lemma} + +\begin{theorem}[Itô] +\label{theorem:ito-existence-uniqueness} + Let $\sigma: [0, \infty) \times C([0, \infty); \real^d) \to L(\real^d; \real^n)$ and $b: [0, \infty) \times C([0, \infty); \real^d) \to \real^n$ be previsible path functionals satisfying the \hyperref[Lipschitz condition]{definition:lipschitz-coefficient}. If for each $T \ge 0$, + \[ + \sup_{0 \le s \le T}\norm{\sigma(s, 0)}_{L(\real^d; \real^n)} + \norm{b(s, 0)}_{\real^n} < \infty + \] + + then the SDE + \[ + X_t = \xi + \int_0^t \sigma(s, X)dB_s + \int_0^t b(s, X)ds + \] + + is exact and pathwise unique. +\end{theorem} +\begin{proof}[Proof of existence, \cite[Theorem 11.2]{Rogers}. ] + For each $t \ge 0$, let $\mathcal{G}_t^\circ = \sigma(\bracs{B_s|1 \le s \le t} \cup \bracs{\xi})$, $\mathcal{G}^\circ = \sigma(\bracsn{\mathcal{G}_t^\circ|t \ge 0})$, $\mathcal{N}$ be the collection of $\bp$-null sets in the completion of $\sigma(\bracsn{\mathcal{G}_t^\circ|t \ge 0})$, and $\mathcal{G}_t = \sigma(\mathcal{G}_t^\circ \cup \mathcal{N})$. + + Let $\xi \in L^\infty(\Omega, \cf_0; \real^n)$ and $X^{(0)} = \xi$. For each $m \in \natz$, define inductively + \[ + X^{(m+1)} = P(X^{(m)}, \xi)_t = \xi + \int_0^t \sigma(s, X^{(m)})dB_s + \int_0^t b(s, X^{(m)})ds + \] + + then $X^{(m)}$ has continuous sample paths, and for each $T > 0$, $\ev[{\normn{X^{(m)}}_{u, [0, s]}^2}] < \infty$. By \autoref{lemma:lipschitz-picard}, + \[ + \ev\braks{\norm{X^{(m+2)} - X^{(m+1)}}_{u, [0, T]}^2} \le C \int_0^t \ev\braks{\norm{X^{(m+1)} - X^{(m)}}_{u, [0, t]}^2} dt + \] + + Therefore + \[ + \ev\braks{\norm{X^{(m+1)} - X^{(m)}}_{u, [0, T]}^2} \le C_0 C \frac{T^n}{n!} + \] + + and $X^{(m)}$ converges uniformly in $L^2$ to a limiting process $X$, which satisfies the given equation. +\end{proof} +\begin{proof}[Proof of uniqueness. ] + Let $X$ and $Y$ be solutions of the SDE with the same setup, then by \autoref{lemma:lipschitz-picard}, + \[ + \ev\braks{\norm{X - Y}_{u, [0, T]}^2} \le C \int_0^t \ev\braks{\norm{X - Y}_{u, [0, t]}^2} dt + \] + + which implies that $X|_{[0, T]} = Y|_{[0, T]}$ almost surely. +\end{proof} + + + diff --git a/src/sde/index.tex b/src/sde/index.tex new file mode 100644 index 0000000..ca391a4 --- /dev/null +++ b/src/sde/index.tex @@ -0,0 +1,6 @@ +\chapter{Stochastic Differential Equations} +\label{chap:sde} + +\input{./setup} +\input{./exact} + diff --git a/src/sde/setup.tex b/src/sde/setup.tex new file mode 100644 index 0000000..06f7267 --- /dev/null +++ b/src/sde/setup.tex @@ -0,0 +1,90 @@ +\section{Definitions} +\label{section:sde-definitions} + +\begin{definition}[Previsible $\sigma$-Algebra] +\label{definition:previsible-sigma-algebra} + Let $(\Omega, \cf)$ be a measurable space and $\bracs{\cf_t|t \ge 0}$ be a filtration on $\Omega$, then the \textbf{previsible $\sigma$-algebra} $\mathscr{J}_\omega$ on $(0, \infty) \times \Omega$ associated with $\bracs{\mathcal{F}_t}$ is the $\sigma$-algebra generated by + \[ + \bracs{(s, t] \times A| 0 \le s < t < \infty, A \in \cf_s} + \] + + In other words, it is the smallest $\sigma$-algebra on $(0, \infty) \times \Omega$ on which every $\bracs{\mathcal{F}_t}$-adapted process with left-continuous sample paths is measurable. +\end{definition} + +\begin{definition}[Previsible Path Functional] +\label{definition:previsible-path-functional} + Let $C([0, \infty); \real^d)$ be the space of $\real^d$-valued continuous functions on $[0, \infty)$, equipped with the topology of uniform convergence. For each $t \ge 0$, let + \[ + \mathscr{X}_t = \sigma(\bracs{\pi_s|s \le t}) + \] + + and $\mathscr{J}_{C([0, \infty); \real^d)}$ be the previsible $\sigma$-algebra on $(0, \infty) \times C([0, \infty); \real^d)$, then a \textbf{previsible path functional} is a $\mathscr{J}$-measurable mapping on $(0, \infty) \times C([0, \infty); \real^d)$. +\end{definition} + +\begin{lemma} +\label{lemma:adapted-composition} + Let $(\Omega, \bracs{\cf_t})$ be a filtered space and $X: \Omega \to C([0, \infty); \real^d)$ be a $\bracs{\mathcal{F}_t}$-adapted process with continuous sample paths, then + \begin{enumerate} + \item $X$ is $(\mathscr{J}_\Omega, \mathscr{J}_{C([0, \infty); \real^d)})$-measurable. + \item For any previsible path functional $\alpha: (0, \infty) \times C([0, \infty); \real^d)$, $\alpha(t, X)$ is $\bracs{\mathcal{F}_t}$-previsible. + \end{enumerate} +\end{lemma} + +\begin{definition}[Diffusion Type SDE] +\label{definition:diffusion-sde} + Let $\sigma: \real^d \to L(\real^d; \real^n)$ and $b: \real^n \to \real^n$ be measurable functions, then a \textbf{SDE of diffusion type} is the equation + \[ + X_t = \xi + \int_0^t \sigma(X_s) dB_s + \int_0^t b(X_s)ds + \] + + under the constraint + \begin{equation} + \int_0^t \norm{\sigma(X_s)}_{\real^n}^2 + \norm{b(X_s)}_{\real^n} ds < \infty \label{equation:diffusion-constraint} + \end{equation} + + for all $t > 0$, where + \begin{itemize} + \item $B$ is a standard Brownian motion on a filtered probability space $(\Omega, \bracs{\cf_t}, \bp)$. + \item $\xi$ is a $\cf_0$-measurable random variable. + \end{itemize} +\end{definition} + +\begin{definition}[Pathwise Uniqueness] +\label{definition:pathwise-uniqueness} + Let $\sigma: \real^d \to L(\real^d; \real^n)$ and $b: \real^n \to \real^n$ be measurable functions, then the diffusion SDE + \begin{equation} + X_t = \xi + \int_0^t \sigma(X_s) dB_s + \int_0^t b(X_s)ds \label{equation:diffusion-sde} + \end{equation} + + has \textbf{pathwise uniqueness} if given + \begin{itemize} + \item A filtered probability space $(\Omega, \bracs{\cf_t}, \bp)$, + \item $\bracs{\mathcal{F}_t}$-adapted standard Brownian motion $B$, + \item Continuous $\bracs{\mathcal{F}_t}$-semimartingales $X, Y: \Omega \to C([0, \infty); \real^d)$ satisfying \autoref{equation:diffusion-sde} and \autoref{equation:diffusion-constraint}, + \end{itemize} + + then $X = Y$ almost surely. +\end{definition} + +\begin{definition}[Pathwise Exact] +\label{definition:pathwise-exact} + Let $\sigma: \real^d \to L(\real^d; \real^n)$ and $b: \real^n \to \real^n$ be measurable functions, then the diffusion SDE + \begin{equation} + X_t = \xi + \int_0^t \sigma(X_s) dB_s + \int_0^t b(X_s)ds \label{equation:diffusion-sde} + \end{equation} + + is \textbf{pathwise exact} if given + \begin{itemize} + \item A filtered probability space $(\Omega, \bracs{\cf_t}, \bp)$, + \item $\bracs{\mathcal{F}_t}$-adapted standard Brownian motion $B$, + \item $\bracs{\mathcal{F}_t}$-semimartingales $X, Y: \Omega \to C([0, \infty); \real^d)$ satisfying \autoref{equation:diffusion-sde} and \autoref{equation:diffusion-constraint}, + \end{itemize} + + then for every $t \ge 0$, $X_t = Y_t$ almost surely. +\end{definition} + + + + + +