From dc42ec6663700862baa12e151b6f3607da4fc388 Mon Sep 17 00:00:00 2001 From: Manos Katsomallos Date: Tue, 12 Oct 2021 23:26:38 +0200 Subject: [PATCH] expt: Testing epsilon percentages --- code/expt/copenhagen-sel-eps.py | 154 +++++++++++++++++++ code/expt/hue-sel-eps.py | 149 +++++++++++++++++++ code/expt/t-drive-sel-eps.py | 179 +++++++++++++++++++++++ code/lib/lmdk_sel.py | 54 +++++++ rslt/lmdk_sel_eps/Copenhagen-sel-eps.pdf | Bin 0 -> 14281 bytes rslt/lmdk_sel_eps/HUE-sel-eps.pdf | Bin 0 -> 16859 bytes rslt/lmdk_sel_eps/T-drive-sel-eps.pdf | Bin 0 -> 16042 bytes 7 files changed, 536 insertions(+) create mode 100644 code/expt/copenhagen-sel-eps.py create mode 100644 code/expt/hue-sel-eps.py create mode 100644 code/expt/t-drive-sel-eps.py create mode 100644 rslt/lmdk_sel_eps/Copenhagen-sel-eps.pdf create mode 100644 rslt/lmdk_sel_eps/HUE-sel-eps.pdf create mode 100644 rslt/lmdk_sel_eps/T-drive-sel-eps.pdf diff --git a/code/expt/copenhagen-sel-eps.py b/code/expt/copenhagen-sel-eps.py new file mode 100644 index 0000000..818e0ec --- /dev/null +++ b/code/expt/copenhagen-sel-eps.py @@ -0,0 +1,154 @@ +#!/usr/bin/env python3 + +import sys +sys.path.insert(1, '../lib') +import argparse +import ast +from datetime import datetime +from geopy.distance import distance +import lmdk_bgt +import lmdk_lib +import lmdk_sel +import exp_mech +import math +import numpy as np +from matplotlib import pyplot as plt +import time + + +def main(args): + # Contacts for all users + cont_data = lmdk_lib.load_data(args, 'cont') + # Contacts for landmark's percentages for all users + lmdk_data = lmdk_lib.load_data(args, 'usrs_data') + # The name of the dataset + d = 'Copenhagen' + # The user's id + uid = '449' + # The landmarks percentages + lmdks_pct = [0, 20, 40, 60, 80, 100] + # The privacy budget + epsilon = 1.0 + eps_pct = [20, 40, 60, 80] + + markers = [ + '^', # 20 + 'v', # 40 + 'D', # 60 + 's' # 80 + ] + + print('\n##############################', d, '\n') + # Get user's contacts sequence + seq = cont_data[cont_data[:, 1] == float(uid)][:1000] + + # Initialize plot + lmdk_lib.plot_init() + # The x axis + x_i = np.arange(len(lmdks_pct)) + plt.xticks(x_i, np.array(lmdks_pct, int)) + plt.xlabel('Landmarks (%)') # Set x axis label. + plt.xlim(x_i.min(), x_i.max()) + # The y axis + plt.ylabel('Mean absolute error (%)') # Set y axis label. + # plt.yscale('log') + plt.ylim(0, 100) + + mae_evt = 0 + mae_usr = 0 + + for i_e, e in enumerate(eps_pct): + mae = np.zeros(len(lmdks_pct)) + + for i, pct in enumerate(lmdks_pct): + # Find landmarks + lmdks = lmdk_lib.find_lmdks_cont(lmdk_data, seq, uid, pct) + + for _ in range(args.iter): + + lmdks_sel = lmdk_sel.find_lmdks_eps(seq, lmdks, epsilon*e/100) + + # Uniform + rls_data, _ = lmdk_bgt.uniform_cont(seq, lmdks_sel, epsilon*(1 - e/100)) + mae[i] += (lmdk_bgt.mae_cont(rls_data)/args.iter)*100 + + # Calculate once + if e == eps_pct[0] and pct == lmdks_pct[0]: + # Event + rls_data_evt, _ = lmdk_bgt.uniform_cont(seq, lmdks, epsilon) + mae_evt += (lmdk_bgt.mae_cont(rls_data_evt)/args.iter)*100 + elif e == eps_pct[-1] and pct == lmdks_pct[-1]: + # User + rls_data_usr, _ = lmdk_bgt.uniform_cont(seq, lmdks, epsilon) + mae_usr += (lmdk_bgt.mae_cont(rls_data_usr)/args.iter)*100 + + # Plot line + plt.plot( + x_i, + mae, + label=str(e/100) + 'ε', + marker=markers[i_e], + markersize=lmdk_lib.marker_size, + markeredgewidth=0, + linewidth=lmdk_lib.line_width + ) + + + plt.axhline( + y = mae_evt, + color = '#212121', + linewidth=lmdk_lib.line_width + ) + plt.text(x_i[-1] + x_i[-1]*.01, mae_evt - mae_evt*.05, 'event') + + plt.axhline( + y = mae_usr, + color = '#616161', + linewidth=lmdk_lib.line_width + ) + plt.text(x_i[-1] + x_i[-1]*.01, mae_usr - mae_usr*.05, 'user') + + path = str('../../rslt/lmdk_sel_eps/' + d) + # Plot legend + lmdk_lib.plot_legend() + # # Show plot + # plt.show() + # Save plot + lmdk_lib.save_plot(path + '-sel-eps.pdf') + print('[OK]', flush=True) + + +def parse_args(): + ''' + Parse arguments. + + Optional: + res - The results archive file. + iter - The total iterations. + ''' + # Create argument parser. + parser = argparse.ArgumentParser() + + # Mandatory arguments. + + # Optional arguments. + parser.add_argument('-r', '--res', help='The results archive file.', type=str, default='/home/manos/Cloud/Data/Copenhagen/Results.zip') + parser.add_argument('-i', '--iter', help='The total iterations.', type=int, default=1) + + # Parse arguments. + args = parser.parse_args() + + return args + + +if __name__ == '__main__': + try: + start_time = time.time() + main(parse_args()) + end_time = time.time() + print('##############################') + print('Time elapsed: %s' % (time.strftime('%H:%M:%S', time.gmtime(end_time - start_time)))) + print('##############################') + except KeyboardInterrupt: + print('Interrupted by user.') + exit() diff --git a/code/expt/hue-sel-eps.py b/code/expt/hue-sel-eps.py new file mode 100644 index 0000000..101ee14 --- /dev/null +++ b/code/expt/hue-sel-eps.py @@ -0,0 +1,149 @@ +#!/usr/bin/env python3 + +import sys +sys.path.insert(1, '../lib') +import argparse +import ast +from datetime import datetime +from geopy.distance import distance +import lmdk_bgt +import lmdk_lib +import lmdk_sel +import exp_mech +import math +import numpy as np +from matplotlib import pyplot as plt +import time + + +def main(args): + # User's consumption + seq = lmdk_lib.load_data(args, 'cons') + # The name of the dataset + d = 'HUE' + # The landmarks percentages + lmdks_pct = [0, 20, 40, 60, 80, 100] + # Landmarks' thresholds + lmdks_th = [0, .54, .68, .88, 1.12, 10] + # The privacy budget + epsilon = 1.0 + eps_pct = [20, 40, 60, 80] + + markers = [ + '^', # 20 + 'v', # 40 + 'D', # 60 + 's' # 80 + ] + + print('\n##############################', d, '\n') + + # Initialize plot + lmdk_lib.plot_init() + # The x axis + x_i = np.arange(len(lmdks_pct)) + plt.xticks(x_i, np.array(lmdks_pct, int)) + plt.xlabel('Landmarks (%)') # Set x axis label. + plt.xlim(x_i.min(), x_i.max()) + # The y axis + plt.ylabel('Mean absolute error (kWh)') # Set y axis label. + plt.yscale('log') + plt.ylim(.1, 100000) + + mae_evt = 0 + mae_usr = 0 + + for i_e, e in enumerate(eps_pct): + mae = np.zeros(len(lmdks_pct)) + + for i, pct in enumerate(lmdks_pct): + # Find landmarks + lmdks = seq[seq[:, 1] < lmdks_th[i]] + + for _ in range(args.iter): + + lmdks = lmdk_sel.find_lmdks_eps(seq, lmdks, epsilon*e/100) + + # Uniform + rls_data, _ = lmdk_bgt.uniform_cons(seq, lmdks, epsilon*(1 - e/100)) + mae[i] += lmdk_bgt.mae_cons(seq, rls_data)/args.iter + + # Calculate once + if e == eps_pct[0] and pct == lmdks_pct[0]: + # Event + rls_data_evt, _ = lmdk_bgt.uniform_cons(seq, lmdks, epsilon) + mae_evt += lmdk_bgt.mae_cons(seq, rls_data_evt)/args.iter + elif e == eps_pct[-1] and pct == lmdks_pct[-1]: + # User + rls_data_usr, _ = lmdk_bgt.uniform_cons(seq, lmdks, epsilon) + mae_usr += lmdk_bgt.mae_cons(seq, rls_data_usr)/args.iter + + # Plot line + plt.plot( + x_i, + mae, + label=str(e/100) + 'ε', + marker=markers[i_e], + markersize=lmdk_lib.marker_size, + markeredgewidth=0, + linewidth=lmdk_lib.line_width + ) + + plt.axhline( + y = mae_evt, + color = '#212121', + linewidth=lmdk_lib.line_width + ) + plt.text(x_i[-1] + x_i[-1]*.01, mae_evt - mae_evt*.14, 'event') + + plt.axhline( + y = mae_usr, + color = '#616161', + linewidth=lmdk_lib.line_width + ) + plt.text(x_i[-1] + x_i[-1]*.01, mae_usr - mae_usr*.14, 'user') + + path = str('../../rslt/lmdk_sel_eps/' + d) + # Plot legend + lmdk_lib.plot_legend() + # Show plot + # plt.show() + # Save plot + lmdk_lib.save_plot(path + '-sel-eps.pdf') + print('[OK]', flush=True) + + +def parse_args(): + ''' + Parse arguments. + + Optional: + res - The results archive file. + iter - The total iterations. + ''' + # Create argument parser. + parser = argparse.ArgumentParser() + + # Mandatory arguments. + + # Optional arguments. + parser.add_argument('-r', '--res', help='The results archive file.', type=str, default='/home/manos/Cloud/Data/HUE/Results.zip') + parser.add_argument('-i', '--iter', help='The total iterations.', type=int, default=1) + + # Parse arguments. + args = parser.parse_args() + + return args + + +if __name__ == '__main__': + try: + start_time = time.time() + main(parse_args()) + end_time = time.time() + print('##############################') + print('Time elapsed: %s' % (time.strftime('%H:%M:%S', time.gmtime(end_time - start_time)))) + print('##############################') + except KeyboardInterrupt: + print('Interrupted by user.') + exit() diff --git a/code/expt/t-drive-sel-eps.py b/code/expt/t-drive-sel-eps.py new file mode 100644 index 0000000..5960074 --- /dev/null +++ b/code/expt/t-drive-sel-eps.py @@ -0,0 +1,179 @@ +#!/usr/bin/env python3 + +import sys +sys.path.insert(1, '../lib') +import argparse +from datetime import datetime +from geopy.distance import distance +import lmdk_bgt +import lmdk_lib +import lmdk_sel +import exp_mech +import numpy as np +from matplotlib import pyplot as plt +import time + + +def main(args): + # The data files + data_files = { + 'T-drive': '/home/manos/Cloud/Data/T-drive/Results.zip', + } + # Data related info + data_info = { + 'T-drive': { + 'uid': 2, + 'lmdks': { + 0: {'dist': 0, 'per': 1000}, # 0.0% + 20: {'dist': 2095, 'per': 30}, # 19.6% + 40: {'dist': 2790, 'per': 30}, # 40.2% + 60: {'dist': 3590, 'per': 30}, # 59.9% + 80: {'dist': 4825, 'per': 30}, # 79.4% + 100: {'dist': 10350, 'per': 30} # 100.0% + } + } + } + # The data sets + data_sets = {} + # Load data sets + for df in data_files: + args.res = data_files[df] + data_sets[df] = lmdk_lib.load_data(args, 'usrs_data') + # Geo-I configuration + # epsilon = level/radius + # Radius is in meters + bgt_conf = [ + {'epsilon': 1}, + ] + eps_pct = [20, 40, 60, 80] + + markers = [ + '^', # 20 + 'v', # 40 + 'D', # 60 + 's' # 80 + ] + + # The x axis + x_i = np.arange(len(list(data_info.values())[0]['lmdks'])) + + for d in data_sets: + print('\n##############################', d, '\n') + args.res = data_files[d] + data = data_sets[d] + # Truncate trajectory according to arguments + seq = data[data[:,0]==data_info[d]['uid'], :][:args.time] + + # Initialize plot + lmdk_lib.plot_init() + # The x axis + plt.xticks(x_i, np.array([key for key in data_info[d]['lmdks']]).astype(int)) + plt.xlabel('Landmarks (%)') # Set x axis label. + plt.xlim(x_i.min(), x_i.max()) + # The y axis + plt.ylabel('Mean absolute error (m)') # Set y axis label. + plt.yscale('log') + plt.ylim(1, 1000000) + + mae_evt = 0 + mae_usr = 0 + + for i_e, e in enumerate(eps_pct): + mae = np.zeros(len(data_info[d]['lmdks'])) + for i, lmdk in enumerate(data_info[d]['lmdks']): + # Find landmarks + args.dist = data_info[d]['lmdks'][lmdk]['dist'] + args.per = data_info[d]['lmdks'][lmdk]['per'] + lmdks = lmdk_lib.find_lmdks(seq, args)[:args.time] + for bgt in bgt_conf: + for _ in range(args.iter): + + lmdks = lmdk_sel.find_lmdks_eps(seq, lmdks, bgt['epsilon']*e/100) + + # Uniform + rls_data_u, _ = lmdk_bgt.uniform_r(seq, lmdks, bgt['epsilon']*(1 - e/100)) + mae[i] += lmdk_bgt.mae(seq, rls_data_u)/args.iter + + # Calculate once + if e == eps_pct[0] and lmdk == min(data_info[d]['lmdks']): + # Event + rls_data_evt, _ = lmdk_bgt.uniform_r(seq, lmdks, bgt['epsilon']) + mae_evt += lmdk_bgt.mae(seq, rls_data_evt)/args.iter + elif e == eps_pct[-1] and lmdk == max(data_info[d]['lmdks']): + # User + rls_data_usr, _ = lmdk_bgt.uniform_r(seq, lmdks, bgt['epsilon']) + mae_usr += lmdk_bgt.mae(seq, rls_data_usr)/args.iter + + # Plot line + plt.plot( + x_i, + mae, + label=str(e/100) + 'ε', + marker=markers[i_e], + markersize=lmdk_lib.marker_size, + markeredgewidth=0, + linewidth=lmdk_lib.line_width + ) + + plt.axhline( + y = mae_evt, + color = '#212121', + linewidth=lmdk_lib.line_width + ) + plt.text(x_i[-1] + x_i[-1]*.01, mae_evt - mae_evt*.05, 'event') + + plt.axhline( + y = mae_usr, + color = '#616161', + linewidth=lmdk_lib.line_width + ) + plt.text(x_i[-1] + x_i[-1]*.01, mae_usr - mae_usr*.05, 'user') + + path = str('../../rslt/lmdk_sel_eps/' + d) + # Plot legend + lmdk_lib.plot_legend() + # # Show plot + # plt.show() + # Save plot + lmdk_lib.save_plot(path + '-sel-eps.pdf') + + +def parse_args(): + ''' + Parse arguments. + + Optional: + dist - The coordinates distance threshold in meters. + per - The timestaps period threshold in mimutes. + time - The total timestamps. + iter - The total iterations. + ''' + # Create argument parser. + parser = argparse.ArgumentParser() + + # Mandatory arguments. + + # Optional arguments. + parser.add_argument('-l', '--dist', help='The coordinates distance threshold in meters.', type=int, default=200) + parser.add_argument('-p', '--per', help='The timestaps period threshold in mimutes.', type=int, default=30) + parser.add_argument('-r', '--res', help='The results archive file.', type=str, default='/home/manos/Cloud/Data/T-drive/Results.zip') + parser.add_argument('-t', '--time', help='The total timestamps.', type=int, default=1000) + parser.add_argument('-i', '--iter', help='The total iterations.', type=int, default=1) + + # Parse arguments. + args = parser.parse_args() + + return args + + +if __name__ == '__main__': + try: + start_time = time.time() + main(parse_args()) + end_time = time.time() + print('##############################') + print('Time elapsed: %s' % (time.strftime('%H:%M:%S', time.gmtime(end_time - start_time)))) + print('##############################') + except KeyboardInterrupt: + print('Interrupted by user.') + exit() diff --git a/code/lib/lmdk_sel.py b/code/lib/lmdk_sel.py index 3420005..73cbc10 100644 --- a/code/lib/lmdk_sel.py +++ b/code/lib/lmdk_sel.py @@ -391,6 +391,60 @@ def find_lmdks(seq, lmdks, epsilon): lmdks_new = seq[lmdks_seq_new - 1] return lmdks_new, epsilon - eps_sel + +def find_lmdks_eps(seq, lmdks, epsilon): + ''' + Add dummy landmarks to original landmarks. + + Parameters: + seq - All of the data points. + lmdks - The original landmarks. + epsilon - The available privacy budget. + + Returns: + lmdks_new - The new landmarks. + ''' + # The new landmarks + lmdks_new = lmdks + if len(lmdks) > 0 and len(seq) != len(lmdks): + # Get landmarks timestamps in sequence + lmdks_seq = find_lmdks_seq(seq, lmdks) + # Turn landmarks to histogram + hist, h = get_hist(get_seq(1, len(seq)), lmdks_seq) + # Find all possible options + opts = get_opts_from_top_h(get_seq(1, len(seq)), lmdks_seq) + # Get landmarks histogram with dummy landmarks + hist_new, _ = exp_mech.exponential(hist, opts, exp_mech.score, 1.0, epsilon) + # Split sequence in parts of size h + pt_idx = [] + for idx in range(1, len(seq), h): + pt_idx.append([idx, idx + h - 1]) + pt_idx[-1][1] = len(seq) + # Get new landmarks indexes + lmdks_seq_new = np.array([], dtype=int) + for i, pt in enumerate(pt_idx): + # Already landmarks + lmdks_seq_pt = lmdks_seq[(lmdks_seq >= pt[0]) & (lmdks_seq <= pt[1])] + # Sample randomly from the rest of the sequence + size = hist_new[i] - len(lmdks_seq_pt) + rglr = np.setdiff1d(np.arange(pt[0], pt[1] + 1), lmdks_seq_pt) + # Add already landmarks + lmdks_seq_new = np.concatenate([lmdks_seq_new, lmdks_seq_pt]) + # Add new landmarks + if size > 0 and len(rglr) > size: + lmdks_seq_new = np.concatenate([lmdks_seq_new, + np.random.choice( + rglr, + size = size, + replace = False + ) + ]) + # Get actual landmarks values + lmdks_new = seq[lmdks_seq_new - 1] + return lmdks_new + + + def test(): # Start and end points of the sequence # # Nonrandom diff --git a/rslt/lmdk_sel_eps/Copenhagen-sel-eps.pdf b/rslt/lmdk_sel_eps/Copenhagen-sel-eps.pdf new file mode 100644 index 0000000000000000000000000000000000000000..1e26a39b055d51494caade545712b47ea3020149 GIT binary patch literal 14281 zcmb_@2{=?=_;8WZG$fL2x%Mb#Gs7UVlQl#_8I0WwW0EaIku4EPvLyRn5|JXZmAyzR zQ7K6(iIRNx4hs4IKL79e?mXwY_rB+xcRT04?|aU9N7z7BT>>R31rsiK4&+zDPyiAj zTRXwz^7GMb= zj2sBIjwE{kH{boRzO@t4h6+g2kpE=|wNwBl5{XRp08q$z4pcyHh$IkzpU0s=As=@G zAW2n>0E8;h)6s@#r~y(3-;5~)l7}0CLL}LEFA6{39;AiPCE7X?l*wKo_YiKJ3{ny+ zB_l131*FgzNhBV809iJGx>FYwXh`%RAE(%W8iiB_YG7WYkdl|E6Jl*q|LTq|pvC}% zx(ldMRiX{q790T`BFUcW0HEfF1!Ts9N+A+lVP03N&2&B3@GRBs#T~@_S5|6HC;8m3 z8^i~x0Ec`7_g%75KHO;LqVW5-g07wVWkeq{cX)b9|331O@#E!3ly4y7ZF_Y}5Ysak zw7#k6?P`VwXY8oP8Rh%7w=blPciXT zY|zCnvoV&(eG`+>N%@mW3EG}z#a;awO!DX9r)4q)NoAq&+ZAFcM+BK~T2mA(72nA=EUL3NqZv|jJape_4N*)Qx}5;vwz-Q<%-=G zY}vL{R1)or8GQYF@eHkJMC`)k`1g%?bGUtoWbuD#v;MAdUu?vyS9)DMqeLPb%C&8u2 zW5<}myIVy@;`qJIJo2Vp4R0qSIP-U;k#=Jk58+Joj@dO|-Pc)ii9I1tBnyFYACs8O!*B7- zq5SU7XmQ+(Jt3u5w7yFHVbGTF_yTvSaek#Rmb~{$%IwF2-tJ_tI!gJWmT1}i;j~1P zj$yB(0;#VyJ&eg5_n~7CYn}D2CJ)Xlyk#!h>wj)F_Ov)(N#36P$gD&qaUwO#R!=)8 z*U0Fb;5K+=X37aA^Mr=5APpUv@q*N;!-2argyAgg#tso?vs=YEaKpxo9rb$#qsEw6 z5^tLfUH^VA!??)Y+-$34d%hwfhVlx#X9{=oQLMFggJgxRxN2ziHcKTw4xXrVZOiN^$qNVY6(7W+?IA)ox(O8 zmT30$?jw}KdiTyB2dSBrkJ7p0O+T5QZL?yTOwUcO`$BlvM*qlYLsDAT&%42_%)*tC z-<-;3y{ewgmfDY}{`&IU@tlvfy=&Ens!OoIZ69q;?v(Ff5S=mUG3dYb826<4`y1DG z1KQKJiGWOYTEf%92N8)wG<+(*(NutCbt~2W>EGH@gHR^Cj`Q$#IKzZfs0~E9f;a-O>5X zHRg2RV!)?dOk$M7a!S1Cnr%|!2I;RI+ro`lhxQ1Zn=1YCZbo0GNj$$vJG^)H%jr3g z<_o{W<0b>7Pq)aX{j)=@KYz~Hm6;dJy%=!$bnD{}-m08;Lla>h-{fD5{MHECg-Us} zwL-LYr^Z9A&tXSiv-*jsk`zw4MicRpJSja0cID`fNOmh^F zM?qnw(2Gu+<`E&FIbgX6LofB5|K)nHc*tMzntEJ_muKIsB2+|<@>Az=F98+waxY3uQ2yLiqX9Xhh87~`jCCoE0HPwEi*IC&*tp2_j1gnOJtkJ4+{+p#T_+y zbQVkYQ;bR7yV1e!h3=EeFYG3_A7LM4dy^kZlz7~*yW!TYeXkqw;rjQeXQt8*FZ84p z9Hxd>8$8ZXCf91BvUiG2=rp`^`D|U~dd4;Xv3n7f>9L&x#{6yD-FE`t8bePW-zd^8 z`|#*GYPsnPLa|)r(PTDxz-)(rR+ery?*C2F(Iw}jOjJrCV8jLe6aS1C+OS5LhBN7_RsX(m~X)!7_)n3 z?_u!g-fPca^T6cXB}-FGq;~p1t&keiIi1c^M`Y@HLweXWa~VhZ#QVu{8sEY)SY01R zzBGz@TIVO(bSuaZmg4FA?K{$rG4Z#4X~3>6|VGdjtmX6Jo$!k{C(Mgq$@e6 z&iaK{##5KThJJx6{Y^#f`~Z8wfu0wVR~z_*V$^H-_S~49-8A&sY0E0gg;o!1Q!Z+m za@9cP*6)XpzGMGcLnr&rD>_&$ymhYs{Ec8^0eM!0sqR6ogx(l=$%oo^T{13zj7{Ap zm#Jf2IURRaY?u<&A&Brhav^lrz9PZ|%ccCoO$vtG$DG4*jp|LDe+X^bFEUL%@PNeH z=6pW>tN{D7JmQ5LGJ6kS%eOFhxWEH+CyjOYXpbsAqTjG*uBpcN8timP?~}eVZ=4AE z>gL}cG<2mpiuoBFBZ>PAx+Z&j?ESVzsLJCsY#Zy*1`RusI?Xz}8t^>@PosN{jP`Q@ zXd?T!j?0~f9TdCJGLKNtrh&?s-HhDs1>~*^($5=LorZRyHe+NObFvuu&Pa{iQ~&H0CR>$(9;3^Xe42W0a?wy_~@{(S=oiY zAIlYk{Oj%>+uj&ne6rB)aQw-XuO&8kZD||6_;$8y_o{JRB z^Ye%}XIRzlx~1oHXPZrNQK~3PEY6g-L-hCowcJMRd4>@`=bQWnwI-oPA_miv9Olgr z2haJA>eavLaya_NOFXWbtm|_r;-!|zwD<7N-AT#%X|F}Ev`@Sfv%{U-W_)etM&YAg z?GI5_mrtUbq+f;=#f+ss-YTkN)8@0s<=1+xQ_gW}PvF+G(G?x+VRyBKj~u!)?HE_7 zZA!A+HuiNWKV3hfB*3_w-i|}&m;b>NVw;5{hVzSrvh74$Bl}+(zis8lpUidf#*oGj zI$nRe?YTxu(aSy@S6_o?0^ z>Z|7HD7Cok{vM|mt2rQRq^Qr#HS*N&8jn35}JZ4;abGr7RuvyEXwZHwBlhF(6?W$_L0(by+HZ;23Ex${nJ zu#M|6B;ArATILR3jznqN?{D_ujS2awuJnUPqq!kbIyeIV?4ldXBZ&;hsw%;@$gSuC zy^v@3BS!Yc5ZR8Z=pK2o^LT>uX}e&>l13T5nx`&po|o!>?KHlcX_H@{2fsP}6FBj~ z>cg(+o-X+|VJ7Vv^@~qsNg)rP<~rXT{2ldPusg?j#?en`Y_vDeU)E$wZC0^hYGT%J zcKx%r?g^_5-kP^gFqq}8AzF!<4h5$lKdF||QLn}Ql3tmoWv*C;;Z{Pc=273rXZbDR zzfS9NCVuT0exI)MGO4m)V`^`qE!I;KIni?7*J`YXzZ7`#>neG2+-$5YL6A(nzsV71 zA#4X*7q!+@XkTe9 zCI*jQE=*J{Qw;#8pG&+iD0GVDX?EVvw3&dUb%u)Hn`2#VHeT~Sd%fm~+vfOAK`*WN ztA~(UkAP~vq0djuE@mdrWfUo8%ITCI+Qer_-hCy>`x4uXYnRH@ri*dh)Kz#~^NInv zj7JO_ZHP2LVm7gcWM#pI)w7na7-}q)G43c-7^Qb$E@wdhBQx!dr^>%z)CN>?wkr zrp*Nj{^Np}tG(A6Y-=-u4xkI42N&BP&)y?o`$?$kvunEG`L1U|%-PTFFr*r>CDQtp>~00YTA72jT!%QkWv$E&%u1~G3wI|*98V&U(l<^R zo)e?Pbj0lR4fic0Bv*;`TBW{0zK};my7u_*a~-gHUXZOkd93qi+BKO;6N#kqyVL$- zJ234NLhL!a9#r%RPZ8s89Lm4foWiXza`~gunPb$U6fSfB&!_jn2^!>6C3T}Ayu)ms znz}ddPKJ!VX!e+3_Pv!_(jPtaX@~RuT~|tyzRY@Hf6EU9q>T8gtfKO@u20b5igTHP zJBSC@X7Ew(GpVn*z4XzE+mq(XW0)BBM171Wp4GH&`8mA z_+2)*_QM8;1~gZ1dC!*M8+XW$zgbyc#7!7(PuX|qn@0OfdRt-z3^$kag9KeUES&r7-c4j%TBd$uv zmWNyFyKR5@7tR#&jsgqAjQXl8Ll^x|95Id6e{BD&sqME7hetxl>4D9Prs9ZqndW7l z**v_ScCRsy^e{)iZ;?!Xd(AM0BJ&|X@h;Z@CgNb&2ZsD)sY5q+psrRYO4#O^?Qrl@ zOGtYBvp>sYQq~?(xY4lcw{h-~V#W9bJ)c$P&srmQndQ?*H|il%;C@Ecabu;;1rO~< z8h37eI(Mi>wHg_X#GKvDdW4Bz#d6Ek`s!oCB0Ec(r}u|lu!uRS0F2}(oshRe&*lU_ znWfA89vaAYj@#@lCuSAS)*Abwuqy@+bR6K!rXRm^pKrea^X8Afw~kDf%+8gCcYcP>Aa}sRSQOnUU6XkXa$)S7V=kionqHu1@bK;<;$x%OS zK>fpObRN={p{8=ERQ=YfI}axr_}x(xtBASQpcRG1F8hG7(Wo2rC`J`uy{lV3My{tH zmj8W0>-y~g?$5)fjDmxqdDb=V&2E}VWq5#&#Dg6Q8rR!D@Wile-idSHD5+nvSu%yb zxjk0*o~yLC*~QIZxfto{kx7PpepJaSGA<3^{wqe4S`h}Nbb_1NKNXwe)Pi05UJPO)0Ms2t)=yrJxW>#=9p>#zC2xw*|R5-YWS z#EFX??-TEBI?yRR@*vwic93h3iM`aT>S6uXb5qPU55F^%^u>IxurlHIC7Ts_t;Z#* zUN-01Ie7MA#}HL@)HV5SXVGq^JDZO0(pO*!A5`+X!4-*ms-~rVQX1*iU}R|`%#CT zkZSE8y_gur6%q^FWFkZ}^xnL{d>!FU+1080B_a%8knlLyPD|>>?+X1MGA1iHyls5*j}7W~a-t>bZaz7qZ2EyMNfvWP#AG*4^rA;tg$Qz6 z{1-+U-i{YRS_%Lzhnjb`AQuB_)Ze z&DF#R-XE2$?6EJu7Di;MH;)^a3ke9DTbEa~IC~Rz5ht#`9sY8}C+JLNM2jn`ktlXF zf^S4cpzPHfy-3R+dk?|$hkiG`NO$aN4bv0j=nt*W&^+4Dw1f1ZNJFvvsD$Z*TBQ_RF#E8Sk<;s6SkW;9_nE&X#?-NPX;^H^;b;CDl~_~4H28qgimR!_&{(xJ zv13DBMTq>ax-$ajUHF?qEKlz`LFZAS1i;UAM(gsae2Z#1bqT3`<+<8kmIuvdT0mjV zhAHJ@D~*q-YUPBcQ*yyd{AJnLG8>HL(><=6-wc#~{J&mjg#`lr_*cPI=?~ zZRUEjJ!?bi7lF$nvNb){0=S${&CGW%aF;6Re{bsMn{~5eyqL$td|S@>?w5zZ290iS z&D(uF?ba(ho|Z3cIwP!XNwDEg$j-X_=ytP^X8lhyKRTJmkU&5 zo65{X@TGlY2|h!upXsE3?X3!`NO;DR``{B#hHOA0sjWZm@W-&rQ6FP;(E)^tctm(qP2#PE%C(ji+oSZ7iuo?-2ZC)h*{{&Sb~xhDSSfwCAGyd!{H@%GHs& z_lZU^UnJ!|mS1UxX>-p-_I&tKd$;wa=9KOJigt@(I=%Bj)8>lrUtaZcjP}esH9CoNWO(oe{?LDpX|`s2f$QT4S&@^)#fp5n?jKgs{8|I&6`>&p z4O*tay6q}OQtNif(wnY=d96_=3bav6Lvr0{ZRR>8qw3TZYe|-76&~(X(L+zveCAv? ziMm#vY9y~4KSf76EPIEo5GU^t#!j9%uy@au6sE`#Pr<|vZWDd&Cnwk*o48+k<6q*< zm~K&Xqm@_S_;adB^$$Kvt{c+sZxb@=cQC%lY(4psy~|jBqty7(uCZaRmWmKu3i+H) zyXaQmz6pAxNL17+LajAUM6U=Gv4+j=b^GBn`LMVekJ|-}i8{OV<^}@4z2RpY+Fa8j z&yw-6y>=u5aq|1@iTkKKsUB$p6%5Bi^OU1sPv(f;j3U!@Z>iq2qai%ER`32D!WHh1 zhDMHO%xkf)+y|czzPmA`MNJe_(PR_7e_Y&_edH0^s2X*mKmApfr!C%|y$^-ylxc~2 zb3F2njC2fKS&=h-Y-dj&-X{ka!TGf%&4HxnS)Ib;)(Gp9la@Ik7a0?37x!K~!C&XO zNA-;g7lX4TjgHuMo$bDTeDW#6sTJ=wKFE9}cfYT;?Q`L0Jx0E+3GvgM@^iYBYp*gU z-;_sB!m{E9r6RpGe1*=j9}si;C0E+9SN&?k*&o{(Uv9B|9g_NPW^3L~hq+Z$yw-ph zjaqICbPe$5>XGo#S(v}odZWav{(=CZ`AzdQvXTQe41$U z^XgovyWFSHw_VcvMP-b;-LP%vTwe)fwM_+{&gTpxAQO6n!eB+i`@$;_I>Fj@_9OXz za(?7Eo2zC&IP!a+OZ5D*&IphaxTxM}mL8&!s+{b3Tfj%4{g}k@OC?3JKq=iZ4U;+J zv7PQ#iY9w0qi0k0b`RJzf|fJucDu|M#}i*>QRf5V-7@uTX|RRYc!R+w$UrY-~Qfnl@6tP>!YaX zrP>nN_H!MlEHyEEDm2jtMEk@=s_TW=gCh?YsZ{RA^4l4VJA_5u#H;oPM{fE~X&Z>0 zGaBv9xNwep|4XjNtFXD&U>6OZt1kWg_ZeuDRKZP?6yee6i zgJ-9MU(mCSz8b9|3hE7+tJe}k2GsL5dD2FgQA+37mLgQImQ@;XDRpaOlV_F2VaDH*UJcHN zy}x76JlY|CCj!$PWOSh=DZIPnq5$=new%%cx?pq5yYsuxQQJmE3!U~A+w+$>xhQ+O zS6IFkY(4-tjVmKEJu9^iV&^uyetLlXP^i@0x%Tv|`JIM#ZkxA8TK}pEwKOtRCU_8O zNofcjM{6R5Ky@UOjKEAUFa^wkK!I}7z=Stix)=~(aw$tWrKnF zTjl@;4IsemFdHI(kp>WML@=!jO!NY+EP|FY=0dasbJ=K*6i0gpum(qK?m zh=3IG-@Dyl1687j4aL!oN~QpyqoO6SLH%?IR0_y3mcic4zsvfK0i#RjFoo}*Jin!rU&(`}e>@V@mI}d5b6$6# zWDwKP$Vfn&Qcxb9y$3jUp!jQ%cr|9hj6^U2($PkVWbZ-*kYFYt!HdQsI2xd_L>N(t zt|pL9=ac&&I!hFXIQ)k-a1a&~G#0ImIUo%#PX0d>*#9dIw6qKW6)>d(z=D$yEHD`v zJb=SV11LNSz<{e}S^;5LFb#p$7KAPKTj;accjX4NHH|x zARP!T#5B}sQQ*Zp$m}1cq(QrfkzU{nWFBe{!LXn$lY;n_!Y@*v#|0WPEYxS7Bj^sp z;O1wU45*030!i^_kf4EDQfS=&nGt_0Xd?(Qy;x|o0TiD$V?g-QJ5AFN460vJ76`#0 z1{MWgd51m)+6aQ9w6t!K?&5u+?$2yk+7?pGsx%e{8NB~Z$y&O#paak{vL`ijj#I*rXn)*=K#?H2c{sX|NdS6YR)j7RJgcUOh0p~hTu_%T8UQ^M zO-H(*Rhd`UUsALliwBO2(SS4{0UBcP_kTO`2EF#O6Ys^KW3>}+j6u^7+Tb&twZ|C? zNlOm~pk*+~kki7#(PEve;zs3FBDkev^1o&c43hr0tPOQx4LY!fbrNYkdlt}_M; z>~XKk-jeY|EV74}l)>Yrz#j&O10Ivi z9cp?TgkjoF3H<6tVos0@H{Qr$%Ehs+p`1x+6|FCsk$=Qm^Mn8+vHu6Ho32Z@kVMAQ zXY9|AJl*nY=)h~bb*bAoUH+Fl$Tlp_Ru%GbFvbJFJ#9x@C^`pCI5f?{Is5l779LB! zKhOj(=?wmR0$K&ALVsJ?ttyK~qETQR0=iEr6jpi<5-o~Eii*JyDo{Hx_}vXgylyUJ zs*9sF;E9oxlEmx<9H>+`4_O4l^-q-~nPLxi0v8mv$8Es>a4hZRW@`sn6KtHp_-_G& z*b?TvC`3CL=s=M$(C;r?01k@rYfyb=b9zXFDW z%;ZWK7!&;i!-JN3Sv&9rbmJ>wPzbOBhDJi>AG-eVj6|Z5|G?0z@r(wSKg;`~&{ANi zu@VM()|D_B$O%H%AG#VhiZ{s~4OM?jO3+7~DU)l9u{+ zPD(?O-3oj%sDJy)V4Dq0-C*3Il#EFkxXe IeRbIX0n1%OvH$=8 literal 0 HcmV?d00001 diff --git a/rslt/lmdk_sel_eps/HUE-sel-eps.pdf b/rslt/lmdk_sel_eps/HUE-sel-eps.pdf new file mode 100644 index 0000000000000000000000000000000000000000..a8f0d792b92e997a5254d2285c7ba3968a551a8b GIT binary patch literal 16859 zcmb_^2|QI@)OYceAykIQkZVeWJGqyO$UH}gWXN<4;a;wJN~Oq@If+b}=O|<-GS5>P zlE_ptmJB8N&N1|)ywBJB{l44pIQ#6g_S$Q$z0Y3%z4zL2smsXm!vzG%xpLk>vdhTf z5Euk+a)n$>3<6a*w!k_-5I~Lw1ghm`hlN0;jUA0`@D>npadIrq3?vl(rGTufqr9f0 zu_G2DxK&6Fk8^||x6Fv8B!~g=-v@|xe!dXDjTg5M`nOwn|Izogtng1&e{8t4eaX7rA0|XA+l0X`e28#nS zFk3w233w+v2&k#FCIl*jb+$6aYRCh%15a86W1NGXF#(G+b=y&Xiyo*8s){wUGM2)- z0@8!L(FlY95(7gZP!K@`QUHb#LV+TygT)3>rrortW(`%gbaB zti#=(>ZwsHe7F_?T?&PIx{jT1F`tgg(_y^rw2QH$XuNDk^e}eh3 z|2C87#KcJ^x&_IXS^w;lhVf1~d!6eu^oXLM2#n9Oc>7+_$g{ioev@0eMqFM=&wS5> z<9DCjmB&4&S2xJKy<;CEid{`LzL4cHKkxss_GxyxFPgh&FjloZ;*sHmx{{+^M@nL{ zVnkc1G}&veQX|vz@W`gLWP(a+)g=2|F6B*K%)_E!*tbxv&$z<6j2WBmTkBFI3nlzb zHKkuE&wFsO%wW#G@eM&Igr(f4Mh+1MZ^fmR(D!_Ty1U!UmW@U&J0sXm&EGiJ&S=}` zB!9`xb><-tTV6}y=eOqcPIpL*&KqtN%qspF%rFT&&^oWj1g^Xh4;n0$49ZSuIEnzBR4!LC(k!&2>8 z>>xU};!yNJ?A*=k5pgT4eL2&x@*7L&sx!~?x_K=gnj+xLMi{!v`V4=DM}Bh|m0b&` zMXol?D+wL6e#zru(c3f2Pee_f+}!=3u$yhoNT$qQA_3V^+(T;7&i2U0v@Cn1&_Bjv za0GoXi(#H`u=P27@cY9v131I?8SCM=a_Epi#;8rLRtx`-htulNfmeK*0n0t=yy}Ut zW+nY-%!DV|fe6S~Vi=vROpt-f-OOjFhpl*Si{!l*O@z>oa~FDw8Z?{=2sZudoRN4* z-U6W}tX)cLP{=-p=DqCrfR}F!J?5L%s?;pJXxAzXXk6cfr*8@8D1qjPpPD|ud_G&* zuuqpZQW(nC#S+SL92LFa)>8x`Z3k;4$V>>4%w;;F+l393OtR*czFw!RkZm_Yn5NoS zJU6Pxh=vrU@yX~VzVL==Li^p-Z^eCekYrL`zpS9{f1!Onjh?<&2IF-{se8 zF(iU=k7DFMtTBLQ7KbB1#2qG*S zH&c7n(!JY)KTMUkvGyR>J-^L6e!H?b#^XynD&JH7RHknz-JCr=gEPkin|$G-)cjOb z<%ITWSQ~%Y91cbER94n%h@Xl5Qk$Jx$#?#eOa9ilE+?1lm_~lyz-e4kiIzZ$M;RUd z>8k)6hrqmYy*3%#d1FJD7LO43OLJ4&&n9|bi<#mmBhsv?O6iz6?UyITuO8=T9(%a| zw1oM&yomJW6BjKHhhG&h_jIv$USDbZGCAP=R_y&tl2`TG^cxwOBrQa2_Etxuk(A_Jp1hVm^k6XG?8UxjYLgWDFbWh zKi4lI*c7?f<~kh^0S%d1^qk*F@^jmZ3N0G2YYYQjYT*a@gsSUq2XiS)3yE`5GQaAE zC;;SD-LTs_-PfW>KOMGY&i>rR-4xM!FmT9W5@`$ehNj#FmpMr^rLR0|c4h7!N_edr4gI9U6Id%#jR3pL?_1 zjP6U1pj47y*>Z2oVP)Nz8$-8ppYv0k{Ipn#FoG#E4SSk9GGfMu)G2)wVNVO;BvUW= zB#niVj#8L7M;H;Xrc~X<_pF18s2`wMnr#U1_pg0f>S}YXE7~&c1O^TY^f>3y($E`7 zwu*bSa}8EuU$;k@eB{p!sI#O$UNY$WHr&cz^tOMX?K4!9y?43o2*UQg68wzs;idD5 zi`VVTRS8e;mipW~v8t_Yta+U;e85SDR?0oj@&bX3$w6)g-86nFnaX7OUiTz5@y7ZP z+AgNpA+^+=gOMxn+9IiA>^Onno&;{+Whfg#86sUmsh>8h&vO!YK2`k9U6*H>h%H6Xlu}rdU_56(4X(xmg>T1=S7a37ZkSDYW8EJ_!%M@*Cx(P6-9L(u9GK~BEed;(q=F6q|qmr>V*N1 zTPC6r%9xRP-lGg8Njvpdh=tJ({*3eAD%cC3YI)9bg$T6{?i9b2quaR@@m)E3 z*D-Z@GbdvDCpsqxwnvx@D8En^dfpq>$d9jOxWn)I{)o@Fb=s(*%9j=WuH){H<0Za- zF!dTJuQ6*FS$eW)QzsWNBL3}6?^g;dj~C{Xgx1$aA2aXD}Qk1az{sANb$xSPqF7^ zafbS;`R^@E^Oxf=y1^)OskfEO6OZP47yVd8XC5uzuja|Bc5&_4SMpY(m*M7*yprBa zU*$tk5;CQX)|#C|LyUbx?mrf!z4`jfPPMMM%`Pz9=w;M$@t-#dEKT~{M)LBbkJC)t z7sxYj>L5v{3Wur&c{L)5fs(ub6u07b&M)!I?#+nwi^NeqZlF zBGJcf<8toWcZEmZm0wv$Tnke{)V&S(w(<8S5Cjn;0)-UZ0p**B7tGXy80_RBenrOr z1uan+5Rd+Zl}Hpy6Q*7OBce`|t0y-tCT;q1`f|)Fr8j*+|RR0>t(oIWpXHF5Pc=qwk=L zN2eJnnPnT9nX@H_2DVYV1_hm6l`1Rnxp3MQ$J>53 z$wix4aJl63@iLLR`>7A&%p*`|_~tU(Ub`)eb{;ZRq>rso9YmE#m%NZo%!VgjT~_yD zyqug<(Dx+KB$vVC)jK&XwVU7*y7cr9_LGO`aDxu3Z!=eBHv9}!<}P3H4Cy<}Q7?(x zLw9f*^f$UAVDMe`M}?N#AB2&U@-gipn()_I&vACZ6jR%08hc1mH^d+uljE;l7tcmAtT+OizuP7KO>#cRCjkdYKAU`y|_Tk0l?hS!?m#7lY z@ilr63Yq9VECNmv|HdLnVZq%N3C9#N!6{@GF$V+^c}>O$C!Xq4eV_zq)}PyLpMUqaKP?f`cYbv2;s8}k!Ag(qqPU3O<$T#NPTy&E-dUBvsK^KO zO@b0N-9i^E8b?m9upL@kYkIw5k=T{#(^gbnx*NWQzrdl;ywI*k`FRiYKvnp;B+?gQ?Dh?xXx>hD{6ehZnB?5Zmu@?YkJ;?bWg# za}(aUJPnw+vh#gWd-hd%sm5}cWzh_`Z$9+yy%*GtOkR=3l4nI%F^`)LpF!mZtD7d< z1ce4Wu9iyBU%0L07~?heJ{aMQ;p}=-UD^I%f@q#V$IE)fHnr7 z5klx)4XJt_gCmxtm|mg#e(ib}s=T61FGkP!su^Q~@uwRVi9)#}fG-RV-km z6FY9{JY;DUyF$iF-qGSh+~Mz66N9#Rea?=mn4(UL&N=<`fvZe>7Hm~@+IQ|5=pciX zlZUH0Whw8dbX~tBT+`>w*2;-5lo+Yz?<{;K%F9z52!D`pAO({0+RuN#{R285Xd`XVjPsxu*2nT&o{sFIOYU>tKV> zL@lNRS}82cMb88 z4QSgK3@!B0Vi%)^>ZqPqjP45)6KGL>W|MM%COqM&XsU`y*=po%o=HM*CkNE?QiR`$ zGkL~~NA6``Y!ugEx@_&AsadOS{gspC?8#L}r8*pSyLD*PZFaiXS=fj|VZMtGv-R~X zBbXsQv9mqB%F~jq#Q&V$Xsq&hNPg3|@72H)H}pw-64{TB@~VQJ`D_QQg0N@QyO;WU zEj$kg%7~%m%^GSE>h;HByL7v{>oL7K{ULptnr9gy2rS)~&ih>&odk2gCk}qjjl*SO zr%0LXbMV~}LT?)Kxh(Roan1S@l9Xg<-X!L(`N}m(jY-z1MIL(e_PM{MdFdW@1y_## z%C3mrzUYiPhKg8{6tL@$vz+aiw?*e8O%^xMz6v~Yuxr7plQ$t0AwF8t-|98}++UcC z_Qc?xT~tRq9oM4n>ny4Yal6J@ocjXxCWTt9+HsYTqd9@JC#q)2Y7dqwtCM36zcP9P zPKFe2$Vek)y(KvHTyDX>#FV_H$gwp%FWaZhz>%|YU%Bz?p01@E+@I=|Vg$3Y6l9JY zX)ixaZ;ax~MH;^QF+6r+{(Qx`LkUNfs?83wX`tH>+)P6nQtD|AXjdY=Y>uIVVY+t% z4jZ)|$GcyeXPp*ROBeU(SQ@@Kc=<9$O1*dwYyG;+_g~-viP>!}R5&J_1x_KoD1HDs z!fW5<96x`bF!I# z{O1BwjIxj^Yg7l#!{?P~2`?46cYIWFNn(RuULYGkfsx20{XqdU`|^IX5Ipx0|^d3VQS`dsaJx255@D{o{IUe*0x z;9JF$t8SAFr(z$dC5~~&bu7N;F-LnH)q1#AnA_T?{4v<*z89iV=&gTV*lfbf!`v#S z?e3>-zVBDOZXGH6io#?)q_mUH|Cut^r3)3SR*_}OIym#Av!6$^lhpzXytPV*&FO`| zUp?=_L&g<2nVrX(ZqD5nH279)qKygTm1*OK#Lb_#dendPjeLCG+W|D=K)m}}9(jJj zgMGREqm@N7&LhO`asqDc z1Y7UX1BL24iF53^-4ZIl?ucQkK^UhlgN$gjkAs76cKmn&Q(^HyVq2U(C|Kq9`V=J{ld&i6vm>E=AqE zcZ`~zo~JI1a6z6~d~w}}XnlD2>7b6s`o4=3TGW5`dgR|uh2Ti!F2l)&NKz5ODTwF2 zWaB$)6`9@=OEEWZ_X!EzTDg^8|ZPb+sbuiMwHCpECYu8rA@_VM(D>sU7 ze0dTm!+15EF07F2u(vQTmW{PcP--OF6Ut8?=^sB{#bqv#6xxSiCL1=PKON^(Gh+QU zwWhP^`w=g(e#t0^bxNW|XHA)O@q}Uj55Kuj=8hwNiykz=H6My!leERRz3xd7?1T?#ZdgciKUawu#NgvO z4>{aYdY1z``oSlnPjuG9E4PDnIZFA)9ya;g*%oT|)B=e{Ys1uG@bkcCjt;D_KnsS3 zT~L6J-v}7L!2n5O7P!ZG^FVm(^~~mX%_B<0`$#*fRQIU;Z)Y|b#P04DTuoU8p&;If zdBMSX{YZa$*0;np@7R4B63b2Dwx$Oky4`+M^~#Pcs*A%_F)HZ-OtBSG$vV35O80K+ zgN>9t$y8C5q6;Lf8u(Ljv2OQh)@-|_S4i$gGCA%c}9{~FO^gm;&>cD zZYUvg-}vNPIHjLvH=E-uU$9N)2lPTDS;J%;R9i&%6t;>$1@r-&`7jB-RFqX$(FQ}RJWtpULR#1KXY{_?asL5fG%A+`; zU@&BK?wIw96LAHxpVl2vKg5Q;e2|C=Z+to{7TVRxoM)wt9?Q@5|kHHb<@on@e^WKW%b;PkG->2h*iqm7Kr_htNbJLm+*&*mC#_GC z&e0u4(m~4j!R|rolS=aVYL5WeNhOTGtig26Q|lux$0EyB%$S)3m(2!CJ{^z8AIs6# zNU1H48@+qg<&sXg+DnV?jqN{7=^dhdZw`}5=;7OG+{AO9Y$1FW%qhnuw8%%Rp60v~h9@HrR zp_O?lUm_~{ockUo`V}Vr7ZMcS-F;MpTPUxcC6$@ik|fJ%v6yOLINZN+p-QF_76L=w zK1F?rl1l;seoz-L>vDDAq@#W;zdx3#5 zxaxKV3@qsS$TuiBCcvJ)V77oCu4OLB^i(Nuu;}=#|MU~DS4>kW9v7k`9K44S`ypR}9D(!JD>CjBje&R{79}$X=o_aGaOzMzw zyz|R4sj_dfQ-DcPc&qT&9TS&-)bWuy>9Es^16`s+?Lnkj(}!*h9y^#W8}vp)!=FEg$>k zZ2)u1tW?}}Y7-mj3oJ2(Qdf&#I{S}(&i2X7Y#NFwQ~Vmq%i}b_+t;Yn#WhuzZXf=R z@f{^yk!yKN?cqBsR8=j@2MPwl7D|n@**x&Nd9M4>F*5h{m>J&PZs{C#l$o}DFwvEF zin4<24egQMENzpC&r$Zl-ZjorQ)g)7(mdsk&;*EcNG>k^Q@d`|~W5$*8z`S@&j-7BIcgnf( zA>`*NmiQB0s-FV=F*(sMGtCtR3x9k}B-Z^#^`s~+93MepJ^X=RG^w1$_0p^e$K=s@ zvakQhRb+(l*F1i8OG@V>O{@s=)=BMCXzsfX{-q~jN25NG3bS+$-B1)S@%Pq=NFxe( zOtwPTZF-BqNN_2p@8ubx;GYyxo8xVxQigYu<|xFDR&=?lK{zQy4-ql*;7Xr0*>kE_ zj+uH|pl{I$-o$Hr)5z%hzD{JDS(Jfw_QV)})iq|q3$2CDHgn|Gq4-R|((~zsN zN6Y^;TKgS~l&Aw+6N3Lx`yQJl_dPiIA6O*rI7P-DHUBMAga5WG!mn0tFAb-dMjenl z@;rIcUf}&Gb^{HWmAfXhkEo4;DMn0%+V~c$uAAG2^#Y5lZz=mB!`_dbFRcBCjXT1R z2scn46k1Q5>y@E+n5Kz<)J4a3+LFSrXx+llf3#q_zgIlToo-&k6DbiSZeig ztNWg~HX~V<8AmVBwVwJR%&LxM#&MYVVb*yym0_QC-KXy=uWsx=cnZDOxAnJZ4=J>J zj)2k-#mTTy%wTRokC1J=6JyHNu$db+uUE?~D2PdDs>BAed@ZA<3xE4LH!xMMX%xD|Vvh7JH}-TcfYR%@Zy} z4yBS1iaT8)s;tsqf}5}3gDJxv5*QBcm$*uZlMJnG&ekMV-#JKUgs1J<`>^criRb{j8CrLjmeB1EthE$m} zy;AIA3hO7WZkQ-xpRZ2ONz1xEY)no5Jb-Bdwy?-Y??cLmEBQD?*`aYag4FtE)#VEm z(PvvdLYl2$XDA{{>6I{x1tOt($D>b~P12UISE3r9==oxb24zngNPy_WTtBS_pp2B{bg7zZ` zST|qHQhBj;)A0vZmG!L*D#{zdS9@0oC_>Uy&4(DxuulS_GbM3N6ea2G`?XmZhFmi~13MJ2w| zaq*NvQ_dVQolM$l%CB5#UTWLl8^2oMMw+Z&RoKSD?)1h{8`xnsU@R1}pNLMWJw`f| z+UE6^u3Jm&py0e=_v|ENbEz*n9)CxrgZr?@z#_3`5IlGfrT!KtB9Xi2-6#zV4(O=U z>*QA<*Uqud4?q+ief4JQj-d4wW>ii&<)X$w7#mAG1M8?yKP)NXn?O=0!_sgm!&z$- z)y%>)klB1u{vGOtfL$7QT(l(eNUL zsQsmv^~cWfzgrMd6VstRpOR$5{w>sttvK_-`5a7#l!lIaQZlE+hxNqDTf!5pBJY!1 z{ZoV#&rC3)gkqHSTU#AW2g|8tSSxJZEa>=;rCi2! zdt3%RvJRKy!q~6lwn&~@(NNiC)HGf~!a-s*75RzzZh5zfPc$tTC$za}E|}V@B_p|~ z=P5tq>sq0*o1}CNJfCT2Pj%mI3Q!uTt66GH9M|}QsO1=cd9mI?+l#-YzYh_U{Y*cc zulrKr!p5@Un1|;&(P29;hd-MIn8N=ZNJjjI-%%P(_BCfI*0RYXs~n!@G{mT!IJYr; z?aMeD%_vz_v)GZ8nU3nIK&aR9y2}fAMS?>jd+7luzbvVcv84>|$H91_op9{=2&Xo)-LTp7=7*FZiXDkQ8V8WiSd z9{GVc`l`wshswfI8}-{b`OfeT6rzKVUkp^YFXS$qBW3Mg_zANCc1^ratmwq*rD3**p9HLwXa;enSI>f7#7`k!=F5F@=QP}RK-Wx z++r%*Q`8e5X_}<_l|H-g4S(-t}V0x2T8 zEU&$xHOpXcB%#emn7$?R5!ZfYvAL>4Dj)=A`;fNnuomz6n^`Dr7qJvxh;2?)Q-tXI zP~OiuOJqNA-m^;wKb5KE?_r<61^1ZU=8+8+rMdVL9VsZHgolPOaI`<7@0wzAd+Q5 z`*8S%=5$|5#2u!yZy8_iA?DvAQ3Mdn^6Ns$8Fgix46ru!p7z<<)U($@7b~*vtIj`X zkxnyscYD=mh?r)2WV#B=ac)#kt{Uq*ESE*%Ok_gD!FHenMwVgpNi}cIVC*=7v9bBJ zMy`YhIz~K8_hD^%9(r%wkKC? zq@3+b4JJ66lqhr}pI|9p7n$6kW72(ea~R*kS)^xOeREx}qP~O4bV5_{Pcy-P9p|Nt z9k82${!kSw6D+~l(F%{#1fu7GuzE{l0vP8H93lWmARy*l76(MJgJJw)V&qU6D|2%! z7-0{@#Opzj2nZC2xHrW@kU|is9To_c2ZHAT=Lg*k!?(ej1Lp%alLRXZOW+L(0=0Cr zv&7;ysjcv45J6#bs5_Q`hd}W-AXMJb1xPts60ktp9EiOK{#rQ$X$LD;2nr2>VmJ8! zF$qu=2vij!3{T)_55Nu=*6F+0gbQ| zFt-Eg^dLfjWBmVck^jH)KnMv#zzYbmgP?#O30&mD!Wak|Ed+sM;1DFR(6M;|c_<(V zYLgbo+bOqQW~c1$ck-S0Kj{G;AfL?(sBhDTU}fMX1dM3lwfPJjqY(lG1D}Dy>#tYaq)()>(= zG5{C&9=wC~Z@>Gg;3gmN{Z2U`4a&I56XXkAAml`#x9*?;z(xc4fOmlJCNGc=@ZJ9_ z$j=N=lb;V@5m2R_b_c}*FOUGZ$dTKuH(C8GwIu`+0vcwsa-bfNEhscdw4?CO9T4`X zP(r|%h!onE3lJWp2lG%U1TaJb`GT08+P8QCqZ$eWRR&=~`t z#b(+Gg4(Q_9PpPw2j0xv^xLSdD*qjJcCdsHgdZ5j zfq(zEa2>$s?~2yhfh>DP>mbz|M-l1^L?#Zm^aTtY4nUgU-OxC$uWzNl&sO4~)E?OI z*EacIiv~cf|E*}F-6-`=luqLWuw$>&_^$hw`nGsXdy4lU9e=kfJH72cth}kplb1yWj|CqIUTB=Sfuy7*T>ET?qly*?&^k!3 z!$`lczJyY853B#uBY&#(*Bt`DasLli*V&h>F93@oPC1(*aI<-2RB4Q8Uji-3{XeAx zzww>kDvfsn5F+sHDO;I=xF2u;U^fT)>=uUGI(7VSR|g_Aeklx`bp_`aGTKls5>Fu=e6Mu!5! z4u8u7vi`6YLKyUC;P0n=a0K#CIzbo)7+ZJI0h)n+=(jw$FwhNu%LDfrf6o&FCV{`@ z3Bu98%M24lfg|*8euBcNKgy!OeWl;%P-qYo>@JIj0f6rJJmeqt16ucwJ_JJx3;oFt z1I%-O<0k}1|3N2&5dK3CAvE|&?v_so10s&ydBTW4+8ibfp5^$R4h?#`-^vOL{!uRm z31W%ebQr-ub;Ae∾ZmfFa_KvT)cQ-xdrF&U|;(0SAQsy}!VL&h`f#(3F4lQ@Ah? zJ^LF!1nkdm6~OVolyxK+TiIX_ zkWfmZB>!`aKFII;@ArTGuh((Tc|Xs4Kj(SBpXa;}pSFUM5JDIY;mdgg&3*zwK;clL z)meyy1Qe!iX@_@%A_0gF6sG6pf``K7Eyk3GM+?$V!&GnOq)cs*2R;dCIBA=C4fHOlMIEaIs-iA*3R;4XH}>f z6sBu$Y2!e!gJM_FjkK)J;;qS0krm?qiG$dvfGmMPB)dTo@Kp)q0cr3A0DxPip-dv4 zbAf`I%IiX53V3%1YrKv!z#X{LBUuvMTr5d=g0V*j zg;6*-5{ZGLktks}P6PvrtPP@**A1wHcO#x7SpyaYO$JzC)uN!0H<%L?ZQcG#4o-l@ zpfDvTz@!RzYoZO{0d+jVj%*J_thxmd#*IwETRKBL<121!5ZOfz1bm)7{`0{l-+eX1 zBH{bX5f;4E{1o5!eU>fczn&0rP|f|<{5)PODreY7<+MD+;@DEFxK@u|#FUjs*Ea{n zuF=#Ct%Z5#lnKmlvpjTIp26iS4~*GHTuX)mTaD7(zG-8m7u+8Gv_AKfa-YXS*xClpIgZUS|e6bbJOkSb*?nTnJiYq8{Wy`oCHIFl2d3uh_7#8D6K_u?n@UjP? z4PN=@BuJWyQgt^Pa30@6x}T3y@Qs}>J1@R8IQm$%->E*(r{E1s45sf%Wm9P92X$3% zbHYKy?q9K9rM@}`-%238A5lW$BXiv0jY?m`JRceRSU>uqT=t4^ij8;nW61t~glJI`E02nLt4449Ff^P0O% zPaU$F*wLupOnqnvqsgui;!Q*$vyn+!AtQy7yOv?uTSkd_wCUSI9|vjtk>xW^&eOx+G9P z?T`Uem_J^kmiq^Fz^ME4tZWndY>xSYte4NTo<|h8Bxn~%!gC4$1PjLv2%}EVrTC=b zIS(1T?|RifXFaMh&rZ8%wE8?sdGGzvn*jWk{WM&Y2u0rhU8gVWPUcLHtSBwDIL$k`Z&<6W&k8PA;;a0=a)oAH+85 zPO)U7s_z*Ir{2zq=) z7~9J3g%>zs4L{j3TOlR`BPJr?Y`vbn}^bGYj>gN+lpp-?NrYRSNtfLfqy{=jA!0tmp=_=^%d0k}5;raDVa!(s$OqH}UOBf;T z$j#pYMmmGEi|1STM&FR3Z*cR~Wuky~t3L1G9Enp_7si+q8%k_rOEq8+|kwN@~R-)1W|9nrWR z=XM9{5?Ab&R_fZpDaWX1EC`A6j#*a%dTbEZ#SF!OZ-Y*eM{{Q356Fs zQpS3T8{R(NlM-)iGR{)$7}-=2Dt_yDPt3D8jjQP@w)2=2xIvuBkl72+M{yHCEDZGUvRsXXC$OXwZ_+A*;iukN+` z+y)UrUE^4QnS4}RImk`qUVH;pp)gs4R9fsx?dIHp?mQAYJ{#ijKyhTCk)tI>?X*L- z*ZDHu`^|mHov&}rMy1Ce2*f%kBZ+k<-9}l-anTjlj)%K=W4+lt;3C35bv`$^ImJ+& z^>EKDz}o5^t^NAYPxcO0JFeI|t#mr99X|0MU6zP&?&U)eM;OK$=Q3@gH7M5b4TXP7 zKpvd3xc#04Hr9qVvi6;TJHJB`#{p)wNDZ53mFR}ifm)8{P?OVceXu`mN71vI-}k6T znSOG9Q}w1K%{yj^qKis+);$j`LQoJnvABE9A$i{Xp(G^$Q{z)YNij*;ZS+9rC$DGi z2OFlIDjlqCMCC~=?*njoqrBP9!430-x)e6h?H#oLkLie&;d?-!hpV=B-fK3%5&>Wa zNTO%4pDYNbzd|;!ee=V2Kge~}&ZA24txXp$JFNMfxoNliPVNegNWC#9Z3ky{^=K$4 zu-@3Tkpn?X%JnfwNL6@#h3pw~si7sg*mV;dyvKtwG?sH5Z<$UHI3?Bq(5-HvAIp8y zQC`t->Ju?OTAafu)TevYV7Eg?!JENZ_Nno?nmO{9#kVbEld*FTpH@_qCJa=RL~6ZQ zpcH?|B&>MbO6Ks9pr*@V2^LFjKYup)!cAvf5JF~;8-{46Cg1+NS@iL^!l*BoUC7r_ z$>s<9d!;Cf`q_@!dKHpz=^?F|G}CMvKin=H_MJHwJT$$_?qrkB({P(x(?w0$G4g53 z=Vljr*qTspg^n6i!L(xSnEErRJbFz&WzBl0%>8C8BHrdgxLb~a7R6c~UU_9tynZ!ZmF&?(c{XKva10!;d=)tU+LrEz z_fbCEnyB&H&yIP*d;E8w{2VfwM>w*4np5p@tKf&}nCji+ihzR0YCa-t@98Ja|uQ8I5ak8g$ z;obWT(#7ck+8c?-=81;Ol#=av#$a;2me;MN*(>=Tm>424aNK(Cy^`31Koh9hIt;ln z0sdbZBnAhjoqw`N6b7RU*Diol97tEHh1ecC1Q7@@wYVUBje_>Q3(M_0+=b08m*{Ss z6|LT(;y3r+1@Dr!bH}}!ebRr+iUf1he@bFSHp_s~mOVhpMN_tTU1I0HeS=+^8G+xI zd#F6T+ihrQ?OSMNS@M~#=RDU{9RL*V${dbR50J=JS!z;awcCff;V82TX0%EI3_xQxdAO* z-6^7PSO4zlJlCGZ#fI0*b{Qv_7gK~XOk&1(2Zy804VtfEh(6L0Nhhf6ZQp9Vdh&(I zptKoNli@{d5-M|^2g(AS0{O;O#+IY#8`71}RT<%m@(h>XJmv+A`2 zPBT_z&R3nYU%KXz>0jDPp-e_wpN?~Xs|)l!M|HGItjT-_x!CZnWxjZbc><#Z)c6p$ z*)@8*`_~TB+Hn4kF=64ZHb2>YBVnRFsdK-Ojd{b&=e;u*eoO51zVu6i>)OJTE?Xe+mwEoEiI7jO?7_DW*%L$$3_oI-Ad zk{3#(nf0%oB**!VjE5rKal9RGsw!F&Mkzj%m;^b`TK!n$tSrBgx{V{jCG5XA5-Eb+ zMSy99=pi@6?QZ&)e&GhBaZph4ICVkG%QVfhR`Ho5Jp z;XRIeiLNu?3YfTr1$81jZnn_NzUf?bg8zVyG)BZuZ;`JlykEC zt?=DiPTmNmD$e5%=9Z|3XU^`~#<<`nz`q$6vB|iKfN@{Q&`pjr{j8)A8}|(NQ@q}? z+<)VNpB}fw0hpo2Db?7X2npdPwWm(0x2GbL_KT;fTRmBbz9uk63T@|s`J9dlJR+HA zIlJ#xwo$#54%->Wp!>Qn3><&(Qp+4&AfK!u9B6gC5p#{3>2(%9>VfD9qkGw=CiYS6 z(5^evUEOMvvdxscjxX0&df$Wi2lTw^d*p>ZO1w+^`=hc3x;>w3w?hbCQoD1mr`yhF zZ?J*{R@tWR1yZ~A(47wBj?P+KcTR73kFKr^D-?-m`qqBCL#Lf&8~DgA(7k@}NyIT4 zHrE_tXOzgBx_mynyi0u3L8QAU6&T+n=Pm{)HOP&~R%=D?DI9nfWN%xtja|Xrk$T_Fu z;my+d(l+b3xK@UgciY^cQU94UI63Xd+gJUh61}R z?unmolirnH{7C(npZm0i?2am2EPeLFGd<}Ux!zO7(pN52KRt7>?s~p&uB}mwZ{pVi zYn+;hHD^p4>tg-weyl$=QAI-;4q!CHs?LvGaaqg`FTnF2!Uur zj&}ZYCl&A4VQ%c6@Ns;|rCntZsC!gSf1h#aDDzOW`H$OhHf8_)hDWFE~U$Q8Ogi?w|fToSmW-TyRmZn9`P6^w)(` z9s;y{!DHEZycxFqEg}5_dZR6DIN$qDUMRxnQx1i_``;)h<_+{=S^E;tFXln=3letZ z_6}DRO}W3LJg+3|)e8AVb8{j8VhZ!^g>C%0aZgfib(>BIgGNJzP;fL9fr2BUI8ivb zL-)V(BNPG!f$5&JBCl+GDiKM}>-&J~dygB_C!t8Sv%6&Jqe*Fo4bC2((pVVWEsfztw<+I+E@R!|XSA&GcNG}@Ea%0G547^k1{p+g zb?@IjtGGvTOj$FV{jhT0{0Cb2I~9KM$UejKDLR%WIr|d-z+;v z(D==}L4>s(cC=}=Mt!^={Ox!Dg=WXJJQu3|!^ql)gpc1{`&|0o^h;GRh1~Gvd|$0$ z*vl97F72&`@9nZQTUBe0z%oVYL8;_U&YM+}tk|+p!a!TW*})dd;b_U0M^n!~Utzwb z7*m#ZW_M(WzE$so`xpWxYadHK4fX@p7Qt) zi()*qC9<#mZO%uAC8H@%X-} z4(4fw#Oo!yoN0Nyz{o7++e*0iHro8_G6sj-tW0PXeI=BRayjlf5AWrDy%||QlNT@E z*`Xu-y&=-sn(Cg{wZh6*F0?TnJf5mCclF__&Cm+Y;h9&)H`5Z9Q}blg#MO)RsX29s z$KvmJ-C|sH?v$UWz8TF%-bTiMUFTpnk5m*AsRP%Bqo@xAq^Cp1l+rgYwW?B~Cn%xx z1`*XURWVdM?iAAHY{U6qR~bm0=qB?>>Zpx4K@br4%`t4)a2DTpg6v&zQFqSFw%~;yTBJy4!QF)}|`;%1LDI8^3(Ja~a3E zZ9egq%{Zs<+~|DG%kV9j{q0Wc4~o8(yuuBVEK6l5GkU}-96DB}`YL7VVeck(!WkQ)ik44gZ_*)5Bh;#P^ zN#{X}H#r${b7wk!Cf^gCGZ4B{{B+^M^dVH+EH6{$k(!b|zIlA~1O4o04T)@06St>i zubv?fC$gGcnDLjSvs5NtE~uUi;22|cSJ8O*bS_}}ZG+n^gZJa4g8uN~Plp_zABiux z^JU2m^IKx@V&a6i{5B^4=k)-vRNoXCPT_z>E#0Ko)ueZ9o?6JbgLjsTBj|S@51op^ zUsJJa{JC%lV)E_gkkh^kp|-MZmZc5u2`^5ORL|3!$c7!cqxmVb_Gq{< z1fmZ;J?^^tH>w6KNl0XZC2AMLA?DKK79qU$CEBxl|T;4vE*B?KGj--O+WK#(tWugYGFKomx|^eJzr;r?`8M--B}E%WoEDH?gxi z2NNX=f8Xa|pO|DaHMnOfCTQW4lsqFYb@AkD9n3by`R9Bdv3WdVv{gyF5Ohx#nS7x7 zO2M)7{DmcvcqqqIC-p^H{P4{S9;Xc> zwO-o&s&D;m&FmH%;6F$!Z72vEPcwPsp25!HZaac%)FKs?kVda zD)>~;$KBZpX#IzW5O*u$glw{m583-D#@>1Pvp?N!PRtILOQloxTkrnqeCe22&GXxs z=wDNF(apn0BeaSBk_?T)XFXZkoF==8y2E>Wm-Qv&DHp-~dK8 z<>&I}oHE=Dv{T-XPtO%BEk6oy*HS3tsQnZcz$<9-90xz_Y20zWZ^qb9xlV7-NZX#xZhT@1Mt(LO#rn%i*u87%f)ZCaA z{;$Ljfx&EgfHC2S2b2gJd2h|Td!|BmQjSUdj#4cwy~s8*W=O*`6qsdI>Du6;0+*w^ zs4i4O}E~JmmP!NPozPNNhyTj^pNRcxoJ>$Aa?Q(=j z#l+n1Y%WB>Hah-m@`T=!17LJz!*~%iW#ddf7O65cH|Ag3KMc%U239axX@tuMteoWa zIpQA3U3mP`J!s@>w%`5x4R7O~sQ!o+6gbx>*i(P9gKwfH!!>es+$g`~J#on-a zhRUYzy9@dvW=bp!xV(wRd7eA5aSFFh*bfg~Yib`RD@-~kjCSN5qc5j@b8ugGmVwpi z*BIB(i`DLO6Osqx(|wfoFN1src!9!WXZKsOdg>aV z`!?3WvQT7Im4~U?r=D-wZYP5ID!8ef8g za5=Fr@7t=PAN>B9Olka+;Zae1Br%H4aqxqX_}wxN&(qUlJY)Mm)Bf0{RGt|wRQ>SP zRXM|tjCaIPSC1MT!}8yB3o1DZ-yidZMwFxd?G;t2;-HI$QRx)dAJWbF@E+gh#r1qE?Q4Hfucjvr4BFhvf0%?(adYWS>I%*)MZVpUwy z=;&p61V4LsbnMIN^H;7u32t;o)ZqopgE=SUxgWiIuNh+Y98cQdP|R=?@WV@|2~!%!LeVs$>gr36bfI!E6>DL3Tlj6aG8jT2L@RpEG=mi>N-( zsccG}$;I$BG_JCMeyE1Vg0;RtM^B+Bu6Z{y0oVe-(%>0jdP`H)ILg~Uf z!ZkSMzlAnlz6Do{f1`L}UrmFtDm1rp*SuW5h4NIAVzFiYWpO`75BT^AP0HYOle2x* z^^AFQ_rHDc*)7E8lSv1s64O@GOMLJ3ZLzSy?m%tQ7w+3f#VWh4xUrd^8W^5Nu@y;a zeXsB0Tyn9cxtT@JP%7^D^h?vPA>GowS;q>KAHTC@Z~VfjK5>BY4kYFK`=6cyKbE`l z=ZXd{T3-rP-)5BlYsQR%Z<+md8k~U)Ww8u`-BLIJJ=$6 zX{zc*i@85nW#pu1T^_VN0C{$uZ3aFwdxF`I<^-Ym<6HVRotsfKj{cQr^yy+{n!Upt z9pI95Q6)1GiBuMh2d(92lRx-SXu6@Z&7 zd+gZ2lHgEL+k9AEZ8`Kp_dE$hx;s(*Ax<~qi?H}qaeM`2DGM#{xDa+&$~^fu9rex6m*%|h`f(e}|7+sBB{f7Lftl%2^+EYO^o}E9l!n_-{@1J%v308;T#ca` zPNOg%Zza61L7ts0iC_Pf;`wD~YJTS@m+Oc-K3}FF7>SiL=3*u6gP4f3Cr=!YPoxi- zaOa6@XEV@Jd*#9S(!llB`wInLG%2Q)4_Y|5&%GfVRQ%vHV|^gvIvSh$;t1i1jc zq{;cY52{l5<04yU)WHBzR<+b)9$JSf5VkSsc*ZNhcUc{132+QxeT-GI5LMAh+ zhLGdI+H2@F3rQ=ltLU(3cwbD-Lv1k;{e}HzS*Mj>>_I+WSYy$r&;!0rnJHaer9!N) zUx+;Mr(voS_nHnas^{vH-|KuszTkQP(x}$Xz*nfLw+9)Q{nRw_MCa)T zGt1x2N4$NODE7Mex^3f@e@&8+Tk>~|PJ?T;4BcWjB)ZbAG^Z|3{fOrB;H7Wxxfq9O zD;p*DrB1a~O$5Vyzb|<_N0cYIC3BbTJ{OoJ7d|qV$^S5vNYS;Yg8ESH_4`$t&&w_2 z*`{=K9j=;GVcxk8y%`#RFsw?B6OdP72(*pV<5V?i@5G)F}H$~L*{;c(}B2V3^)37+!L!Wei+qzXcvb1EC6#K&(4e$APq z{Y|(yJxBHBiF*Dv_W9T39=F*%ilO2Rm*BwGbo_JW`eNr+N;$gG-^cG#AQT=?@$=uR zDiCW6Yrkx!f;wKJf;`FJCwR2t1uv6dh*6&W6B!Jbt@da8pwNdng?_&f>hGl1!N_IZ z$)41xFgBS1)|Y>a2_%^RmWd*PI)RNFC6d}|1O;GkYW(0+nKYS8H)hMTZfksg#vz|> zHgs*l?=2fOtfjc9%gd)7J{ zd40!RZ9XJe9*kXiTtya5CC7y=H@f*j&$gPw7xo0Kw?YJJDb&_*; zEioJAX*fw|82t!O|GLQP3KN@gq5mMUiMPnavC4nRq`bC`&3aT<^=~sFHhNRe(ha{- zD-Bb3u)>op$qqzM6BR3_V7lE6}Fpbi_fD^ND9NB|0f!AfWe2?$KV!PXWJ7AXT| z#U@Y`5()!~maXwnln4~&f(NRNf%;=0_`y~xnVs;qK#BFrCCS0g9%#X?^l~8DK+&R5 z*m*pO2!#;|cnFN_0bG&oNqC^**cK>T2L3y^1J`a2o=^;cgB_i2!R74DvE<*u_90e4gp00dm1YT0K)*)H7k7q*gD=C&N}YaGi1H}Pk(?0NN42$ z_^sFwWCk80z~luUEBC-h2O@xA;5~p0!2k~d2nPv@0s^e`1=?5e*82UyTcd#kUQe#T zK>7f^H6bA28I1)L0%(ECtqTeASic7Q{J92Y01mJnJcIn#+WsiGLI-SL#{;fG8CPh6 zbb$lRE|J*PGiU(NX&@c&4A5Pn1=0Z?`2Pg?0{}Joa|5D)Dy{n+6bn4S4!{9Ht&v_K z^#^NJ2ow}F%nEaW4@edi8tk;L@cJ1L_K#2^!1RX_S(6J89_$aoFc>5-8v<}NZk_uo zEnp(UfHIa0I0t*VgUHY z^NOWG7}&nySYreM1z1;j>lu6$Ry+uJX=B?u-}Uoa+n>>}u`g(tZF#Ia894t}X}B@a zt3y``4|Tw300<7uz+gl`u6_Gs1x2oX5o7}eK9X5!#RKCDyn?I+6Xfb7x*BMJ@|2ZM zE3s$g%_cz45sCumaNx=b40S74&VaL4sSu!mr&eG@V2W6~A^@Rqh0Dee2zYv>--bwt zS_>NMF%z}MAGlwMjw`Vj3|{Cp3EY4$eO9S}ArZYM!)ovaM(hfGAi%E4u^LH%aJvcz zf-8DWhI3%#U4gno0pkN#kkvR0JQN`KznOT`1_lC{@wT*YTxQE6RzGW!1^TQ9*gxN9 z!K^KzD^*o6S>$>*#Q(Y51`#%_xBrt1RBG+D^MAngr~!UlrGyaL98-mK@OI~%EJ;Fm z7dHndA_0nAGl2#isLohX3#I{RxMD^ih#{ds@LI8v1~B2Q8u33?D;VqVz1I^35vUL_ zM+1NVZ$%zJK;Km3v7V1?SL1=wt{+Bf&rn#oT{RUpbK4DV9J->zVruGOy2DwTN^YAR zZ=+BCuSEm$r~g~ghC4CZ?HI%QQKTWVTI2;=ILt(K{ zp-p%H=Nll6`(JNVfPuYPmG#l~Kdii>%A=P>1rO1MX#M)pD=7Mdq)T622n?~p*Y&8F z4Oy9Y)E3jrZe#TgBl3@0|2!aoB<=r$)eU!~m*;8|{0rPiYb_a(!FlJX1wABwV{~PKUu>B7*@Yx@@PEY{9w>+jG zhJ+&#K#GDyqL3IQ<~SV54~O#$Kw$DrC5#qEoq*bt$u4eU zFqre7CSfAU4nP8y95&~yfsYS2V!7DZLai*V9f8zvje@|2$-$F^w}k+)6$t@Sg0&wM zi@^YMu`Tou419dSAC#~LL!d=~9oH=|U<)0X5;pfki2xCOD@+tfjsAk+{_ckZO}mL6 za05j6tuQb<*aAbr|B??5N5Z$t2M0R+uP|T5f@1{C=FaN@^NZemCW6^)}z+!;ixh;ILILuamz(tUM zjX$^u224u8-yi?NMX;iOV?Y|GLnMOJ-%4oLbvTa)CXdw}Axjc@WjaBiall%Gk55rc3G)8{r<}al literal 0 HcmV?d00001