From b2efdd75114b0b94b06567c48aa283a83fe5b7ae Mon Sep 17 00:00:00 2001 From: Bruno Stegani Date: Tue, 21 Jul 2026 18:23:03 +0200 Subject: [PATCH 1/3] updated dihedrals: first baseline with high beta propensity --- multi-ego-basic.ff/aminoacids.rtp | 4 + multi-ego-basic.ff/ffbonded.itp | 122 ++++++++++++++++-------------- 2 files changed, 71 insertions(+), 55 deletions(-) diff --git a/multi-ego-basic.ff/aminoacids.rtp b/multi-ego-basic.ff/aminoacids.rtp index d92da36a..220febd8 100644 --- a/multi-ego-basic.ff/aminoacids.rtp +++ b/multi-ego-basic.ff/aminoacids.rtp @@ -202,6 +202,7 @@ ; ai aj ak al gromos type -CA -C N CA gd_14 -C N CA C gd_44N + -C N CA C gd_44Nb -C N CA C gd_43N N CA CB CG gd_34 N CA C +N gd_45N @@ -552,6 +553,7 @@ ; ai aj ak al gromos type -CA -C N CA gd_14 -C N CA C gd_44H + -C N CA C gd_44Hb -C N CA C gd_43H N CA CB CG gd_34 N CA C +N gd_45H @@ -710,6 +712,7 @@ ; ai aj ak al gromos type -CA -C N CA gd_14 -C N CA C gd_44K + -C N CA C gd_44Kb -C N CA C gd_43K N CA CB CG gd_34 N CA C +N gd_45K @@ -927,6 +930,7 @@ ; ai aj ak al gromos type -CA -C N CA gd_14 -C N CA C gd_44S + -C N CA C gd_44Sb -C N CA C gd_43S N CA CB OG gd_34 N CA C +N gd_45S diff --git a/multi-ego-basic.ff/ffbonded.itp b/multi-ego-basic.ff/ffbonded.itp index eac40834..b4d4c66a 100644 --- a/multi-ego-basic.ff/ffbonded.itp +++ b/multi-ego-basic.ff/ffbonded.itp @@ -391,16 +391,19 @@ ;;;;;;;;;;;;;;;;;;;; ; ASN ;;;;;;;;;;;;;;;;;;;; -#define gd_42N 130.0000 1.5000 2 +#define gd_42N 100.0000 1.5000 2 ; Backbone dihedral angle -N-CA-C-N- ASN ; -#define gd_43N 15.0000 0.600 3 +#define gd_43N -35.0000 1.700 3 ; Backbone dihedral angle -C-N-CA-C- ASN ; -#define gd_44N -110.000 3.3 1 +#define gd_44N -85.000 2.7 1 ; Backbone dihedral angle -C-N-CA-C- ASN ; -#define gd_45N 180.00 2.6 1 +#define gd_44Nb 30.000 2.2 4 +; Backbone dihedral angle -C-N-CA-C- ASN +; +#define gd_45N 0.00 1.6 1 ; Backbone dihedral angle -N-CA-C-N- ASN ; ;;;;;;;;;;;;;;;;;;;; @@ -421,31 +424,32 @@ ;;;;;;;;;;;;;;;;;;;; ; CYS ;;;;;;;;;;;;;;;;;;;; -#define gd_42C 98.0000 2.0000 2 +#define gd_42C 30.0000 2.0000 2 ; Backbone dihedral angle -N-CA-C-N- CYS ; -#define gd_43C 70.0000 4.300 3 ; lower k for higher J coupling and more population in intermediate beta region +#define gd_43C 70.0000 2.200 3 ; lower k for higher J coupling and more population in intermediate beta region ; Backbone dihedral angle -C-N-CA-C- CYS ; -#define gd_43Cb -24.0000 5.500 2 ; lower k for higher J coupling and more population in intermediate beta region +#define gd_44C 68.000 3 1 ; Backbone dihedral angle -C-N-CA-C- CYS ; -#define gd_44C 68.000 3.0 1 +#define gd_43Cb -10.0000 2.700 2 ; lower k for higher J coupling and more population in intermediate beta region ; Backbone dihedral angle -C-N-CA-C- CYS ; -#define gd_45C -125.00 1.0 1 +; +#define gd_45C 0.00 5.50 1 ; Backbone dihedral angle -N-CA-C-N- CYS ; ;;;;;;;;;;;;;;;;;;;; ;GLU ;;;;;;;;;;;;;;;;;;;; -#define gd_42E 120.0000 2.5000 2 +#define gd_42E 90.0000 2.5000 2 ; Backbone dihedral angle -N-CA-C-N- GLU ; #define gd_43E 15.0000 2.60 3 ;heigher k to shift down J coupling and have a cleaner distinction between betas ; Backbone dihedral angle -C-N-CA-C- GLU ; -#define gd_44E -120.000 1.50 1 +#define gd_44E -160.000 1.50 1 ; Backbone dihedral angle -C-N-CA-C- GLU ; #define gd_45E 180.00 2.900 1 @@ -454,28 +458,30 @@ ;;;;;;;;;;;;;;;;;;;; ; HIS ;;;;;;;;;;;;;;;;;;;; -#define gd_42H 150.0000 2.5000 2 +#define gd_42H 80.0000 2.5000 2 ; Backbone dihedral angle -N-CA-C-N- HIS ; -#define gd_43H 15.0000 0.20 3 +#define gd_43H -50.0000 2.70 3 ; Backbone dihedral angle -C-N-CA-C- HIS ; -#define gd_44H -110.000 3.10 1 +#define gd_44H -25.000 1.60 1 +; Backbone dihedral angle -C-N-CA-C- HIS +#define gd_44Hb 40.000 2.40 4 ; Backbone dihedral angle -C-N-CA-C- HIS ; -#define gd_45H 160.00 0.5 1 +#define gd_45H 0.00 3.0 1 ; Backbone dihedral angle -N-CA-C-N- HIS ; ;;;;;;;;;;;;;;;;;;;; ; LEU ;;;;;;;;;;;;;;;;;;;; -#define gd_42L 120.0000 2.5000 2 +#define gd_42L 90.0000 2.5000 2 ; Backbone dihedral angle -N-CA-C-N- LEU ; #define gd_43L 15.0000 0.600 3 ; Backbone dihedral angle -C-N-CA-C- LEU ; -#define gd_44L -140.000 2.5 1 +#define gd_44L -160.000 2.5 1 ; Backbone dihedral angle -C-N-CA-C- LEU ; #define gd_45L 180.00 2.0 1 @@ -484,7 +490,7 @@ ;;;;;;;;;;;;;;;;;;;; ; MET ;;;;;;;;;;;;;;;;;;;; -#define gd_42M 120.0000 2.5000 2 +#define gd_42M 90.0000 2.5000 2 ; Backbone dihedral angle -N-CA-C-N- MET ; #define gd_43M 15.0000 0.600 3 @@ -499,13 +505,13 @@ ;;;;;;;;;;;;;;;;;;;; ; ARG ;;;;;;;;;;;;;;;;;;;; -#define gd_42R 120.0000 2.5000 2 +#define gd_42R 95.0000 2.0000 2 ; Backbone dihedral angle -N-CA-C-N- ARG ; #define gd_43R 15.0000 1.600 3 ; Backbone dihedral angle -C-N-CA-C- ARG ; -#define gd_44R -110.000 2.50 1 +#define gd_44R -150.000 2.50 1 ; Backbone dihedral angle -C-N-CA-C- ARG ; #define gd_45R 180.000 2.50 1 @@ -514,31 +520,31 @@ ;;;;;;;;;;;;;;;;;;;; ; TRP ;;;;;;;;;;;;;;;;;;;; -#define gd_42W 98.0000 2.0000 2 +#define gd_42W 30.0000 2.0000 2 ; Backbone dihedral angle -N-CA-C-N- TRP ; -#define gd_43W 70.0000 4.300 3 ; lower k for higher J coupling and more population in intermediate beta region +#define gd_43W 80.0000 2.200 3 ; lower k for higher J coupling and more population in intermediate beta region ; Backbone dihedral angle -C-N-CA-C- TRP ; -#define gd_43Wb -24.0000 5.500 2 ; lower k for higher J coupling and more population in intermediate beta region +#define gd_43Wb -10.000 2.700 2 ; lower k for higher J coupling and more population in intermediate beta region ; Backbone dihedral angle -C-N-CA-C- TRP ; #define gd_44W 68.000 3.0 1 ; Backbone dihedral angle -C-N-CA-C- TRP ; -#define gd_45W -125.00 2.0 1 +#define gd_45W 0.00 6.0 1 ; Backbone dihedral angle -N-CA-C-N- TRP ; ;;;;;;;;;;;;;;;;;;;; ; GLN ;;;;;;;;;;;;;;;;;;;; -#define gd_42Q 120.0000 1.5000 2 +#define gd_42Q 90.0000 1.5000 2 ; Backbone dihedral angle -N-CA-C-N- GLN ; #define gd_43Q 15.0000 1.000 3 ; Backbone dihedral angle -C-N-CA-C- GLN ; -#define gd_44Q -140.000 2.60 1 +#define gd_44Q 180.000 2.60 1 ; Backbone dihedral angle -C-N-CA-C- GLN ; #define gd_45Q 180.000 2.50 1 @@ -547,122 +553,128 @@ ;;;;;;;;;;;;;;;;;;;; ; LYS ;;;;;;;;;;;;;;;;;;;; -#define gd_42K 150.0000 1.5000 2 +#define gd_42K 100.0000 2.3000 2 ; Backbone dihedral angle -N-CA-C-N- LYS ; -#define gd_43K 15.0000 2.000 3 +#define gd_43K 30.0000 2.000 3 ; Backbone dihedral angle -C-N-CA-C- LYS ; -#define gd_44K -120.000 1.00 1 +#define gd_44K 0.000 1.50 1 ; Backbone dihedral angle -C-N-CA-C- LYS ; -#define gd_45K 180.000 1.40 1 +#define gd_44Kb 65.000 1.80 4 +; Backbone dihedral angle -C-N-CA-C- LYS +; +#define gd_45K 0.000 3.0 1 ; Backbone dihedral angle -N-CA-C-N- LYS ; ;;;;;;;;;;;;;;;;;;;; ; ILE ;;;;;;;;;;;;;;;;;;;; -#define gd_42I 98.0000 2.0000 2 +#define gd_42I 30.0000 2.0000 2 ; Backbone dihedral angle -N-CA-C-N- ILE ; -#define gd_43I 70.0000 4.300 3 ; lower k for higher J coupling and more population in intermediate beta region +#define gd_43I 80.0000 2.200 3 ; lower k for higher J coupling and more population in intermediate beta region ; Backbone dihedral angle -C-N-CA-C- ILE ; -#define gd_43Ib -24.0000 5.500 2 ; lower k for higher J coupling and more population in intermediate beta region +#define gd_43Ib -10.0000 2.700 2 ; lower k for higher J coupling and more population in intermediate beta region ; Backbone dihedral angle -C-N-CA-C- ILE ; #define gd_44I 68.000 3.0 1 ; Backbone dihedral angle -C-N-CA-C- ILE ; -#define gd_45I -125.00 1.0 1 +#define gd_45I 0.00 6.0 1 ; Backbone dihedral angle -N-CA-C-N- ILE ; ;;;;;;;;;;;;;;;;;;;; ; VAL ;;;;;;;;;;;;;;;;;;;; -#define gd_42V 98.0000 2.0000 2 +#define gd_42V 30.0000 2.0000 2 ; Backbone dihedral angle -N-CA-C-N- VAL ; -#define gd_43V 70.0000 4.300 3 ; lower k for higher J coupling and more population in intermediate beta region +#define gd_43V 95.0000 2.200 3 ; lower k for higher J coupling and more population in intermediate beta region ; Backbone dihedral angle -C-N-CA-C- VAL ; -#define gd_43Vb -24.0000 5.500 2 ; lower k for higher J coupling and more population in intermediate beta region +#define gd_43Vb -10.0000 2.700 2 ; lower k for higher J coupling and more population in intermediate beta region ; Backbone dihedral angle -C-N-CA-C- VAL ; -#define gd_44V 68.000 3.0 1 +#define gd_44V 50.000 3.0 1 ; Backbone dihedral angle -C-N-CA-C- VAL ; -#define gd_45V -125.00 0.0 1 +#define gd_45V 0.00 6.0 1 ; Backbone dihedral angle -N-CA-C-N- VAL ;;;;;;;;;;;;;;;;;;;; ; PHE ;;;;;;;;;;;;;;;;;;;; -#define gd_42F 98.0000 2.0000 2 +#define gd_42F 30.0000 2.0000 2 ; Backbone dihedral angle -N-CA-C-N- PHE ; -#define gd_43F 70.0000 4.300 3 ; lower k for higher J coupling and more population in intermediate beta region +#define gd_43F 80.0000 2.200 3 ; lower k for higher J coupling and more population in intermediate beta region ; Backbone dihedral angle -C-N-CA-C- PHE ; -#define gd_43Fb -24.0000 5.500 2 ; lower k for higher J coupling and more population in intermediate beta region +#define gd_43Fb -10.0000 2.700 2 ; lower k for higher J coupling and more population in intermediate beta region ; Backbone dihedral angle -C-N-CA-C- PHE ; #define gd_44F 68.000 3.0 1 ; Backbone dihedral angle -C-N-CA-C- PHE ; -#define gd_45F -125.00 0.7 1 +#define gd_45F 0.00 5.7 1 ; Backbone dihedral angle -N-CA-C-N- PHE ;;;;;;;;;;;;;;;;;;;; ;SER ;;;;;;;;;;;;;;;;;;;; -#define gd_42S 120.0000 2.5000 2 +#define gd_42S 110.0000 2.9000 2 ; Backbone dihedral angle -N-CA-C-N- SER ; -#define gd_43S 15.0000 2.000 3 +#define gd_43S -15.0000 2.000 3 +; Backbone dihedral angle -C-N-CA-C- SER +; +#define gd_44S -10.000 1.50 1 ; Backbone dihedral angle -C-N-CA-C- SER ; -#define gd_44S -140.000 1.50 1 +#define gd_44Sb 40.000 1.70 4 ; Backbone dihedral angle -C-N-CA-C- SER ; -#define gd_45S 180.00 0.60 1 +#define gd_45S 0.00 3.00 1 ; Backbone dihedral angle -N-CA-C-N- SER ; ;;;;;;;;;;;;;;;;;;;; ;THR ;;;;;;;;;;;;;;;;;;;; -#define gd_42T 98.0000 2.0000 2 +#define gd_42T 30.0000 2.0000 2 ; Backbone dihedral angle -N-CA-C-N- THR ; -#define gd_43T 70.0000 4.300 3 ; lower k for higher J coupling and more population in intermediate beta region +#define gd_43T 80.0000 2.200 3 ; lower k for higher J coupling and more population in intermediate beta region ; Backbone dihedral angle -C-N-CA-C- THR ; -#define gd_43Tb -24.0000 5.500 2 ; lower k for higher J coupling and more population in intermediate beta region +#define gd_43Tb -10.0000 2.700 2 ; lower k for higher J coupling and more population in intermediate beta region ; Backbone dihedral angle -C-N-CA-C- THR ; #define gd_44T 68.000 3.0 1 ; Backbone dihedral angle -C-N-CA-C- THR ; -#define gd_45T -125.00 1.5 1 +#define gd_45T 0.00 6.0 1 ; Backbone dihedral angle -N-CA-C-N- THR ;;;;;;;;;;;;;;;;;;;; ; TYR ;;;;;;;;;;;;;;;;;;;; -#define gd_42Y 98.0000 2.0000 2 +#define gd_42Y 30.0000 2.0000 2 ; Backbone dihedral angle -N-CA-C-N- TYR ; -#define gd_43Y 70.0000 4.300 3 ; lower k for higher J coupling and more population in intermediate beta region +#define gd_43Y 80.0000 2.200 3 ; lower k for higher J coupling and more population in intermediate beta region ; Backbone dihedral angle -C-N-CA-C- TYR ; -#define gd_43Yb -24.0000 5.500 2 ; lower k for higher J coupling and more population in intermediate beta region +#define gd_43Yb -10.0000 2.700 2 ; lower k for higher J coupling and more population in intermediate beta region ; Backbone dihedral angle -C-N-CA-C- TYR ; #define gd_44Y 68.000 3.0 1 ; Backbone dihedral angle -C-N-CA-C- TYR ; -#define gd_45Y -125.00 1.5 1 +#define gd_45Y 0.00 5.5 1 ; Backbone dihedral angle -N-CA-C-N- TYR From f4511b70b2708ce3a111b938c295b4e1aa88650d Mon Sep 17 00:00:00 2001 From: Bruno Stegani Date: Tue, 21 Jul 2026 18:24:02 +0200 Subject: [PATCH 2/3] updated tools for better handling of datadriven interaction matrices --- src/multiego/interaction_matrix.py | 47 +++++++++++++++++++++++------- src/multiego/type_definitions.py | 9 ++++-- tools/make_mat/make_mat.py | 35 +++++++++++++++------- 3 files changed, 68 insertions(+), 23 deletions(-) diff --git a/src/multiego/interaction_matrix.py b/src/multiego/interaction_matrix.py index ce295117..f6ad4195 100644 --- a/src/multiego/interaction_matrix.py +++ b/src/multiego/interaction_matrix.py @@ -21,13 +21,13 @@ def find_class(self, module, name): class InteractionMatrix: - def __init__(self, emax=0.16, pth=None): - self.atmat = self.read_pickle_file("atdhisto.pkl") + def __init__(self, pkl_file, emax=0.16, pth=None, show=False): + self.atmat = self.read_pickle_file(pkl_file) self.emax = emax self.pth = pth self.define_p_threshold() self.nonlocal_matrix() - self.plot_energy_matrix() + self.plot_energy_matrix(show=show) def roundPartial(self, value, resolution=0.005): """ @@ -51,15 +51,17 @@ def regall(self, ps, preference, emax, pth): e = np.array([self.roundPartial(ei) for ei in e]) return e - def plot_energy_matrix(self, out_path = None): + def plot_energy_matrix(self, out_path = None, show = False): """ Plot the energy matrix as a heatmap. """ # TODO read them by gromos to avoid hardcoding them here + bkbnd_atoms = ["C", "O", "N", "CAH", "CAH2", "H"] attype_ordering_charged = [ "OM", "NL", "NZ"] - attype_ordering_polar = [ "NE", "NR", "NT", "O", "OA", "H","N", "CZ", "S", "SH", "CAH", "CH2", "CAH2", "CH1t"] - attype_ordering_apolar = ["C", "CH", "CH1", "CH3"] + attype_ordering_polar = [ "NE", "NR", "NT", "OA", "CZ", "S", "SH", "CH2r", "CH1t"] + attype_ordering_apolar = ["CH1","CH", "CH2", "CH3"] + list_atom1 = np.sort(np.array(self.atmat["attype1"].unique())) @@ -72,8 +74,8 @@ def plot_energy_matrix(self, out_path = None): list_atom2 = [at for at in list_atom2 if at not in noplot] # create ordered list of atom types based on attype_ordering - order_list_atom1 = [at for at in attype_ordering_charged if at in list_atom1] + [at for at in attype_ordering_polar if at in list_atom1] + [at for at in attype_ordering_apolar if at in list_atom1] - order_list_atom2 = [at for at in attype_ordering_charged if at in list_atom2] + [at for at in attype_ordering_polar if at in list_atom2] + [at for at in attype_ordering_apolar if at in list_atom2] + order_list_atom1 = [at for at in bkbnd_atoms if at in list_atom1] + [at for at in attype_ordering_charged if at in list_atom1] + [at for at in attype_ordering_polar if at in list_atom1] + [at for at in attype_ordering_apolar if at in list_atom1] + order_list_atom2 = [at for at in bkbnd_atoms if at in list_atom2] + [at for at in attype_ordering_charged if at in list_atom2] + [at for at in attype_ordering_polar if at in list_atom2] + [at for at in attype_ordering_apolar if at in list_atom2] energies = [] for at1 in order_list_atom1: @@ -121,7 +123,8 @@ def plot_energy_matrix(self, out_path = None): ax.grid() if out_path is not None: plt.savefig(out_path) - plt.show() + if show: + plt.show() def read_pickle_file(self, filename): @@ -149,16 +152,36 @@ def read_pickle_file(self, filename): probs = [] dists = [] cutoffs = [] + bins = [] + kde = [] + tot_reps = [] + sum_probs = [] for key in data.keys(): probs.append(data[key].p_repeats) dists.append(data[key].exp_aver) cutoffs.append(data[key].cutoff) + bins.append(data[key].xbins) + kde.append(data[key].kde) + tot_reps.append(data[key].tot_repeats) + sum_probs.append(np.sum(data[key].sum_probs)) atmat["probability"] = np.array(probs)#/np.array(cutoffs)**2 # atmat["p_water"] =np.zeros(len(probs)) atmat["exp_aver"] = np.array(dists) atmat["attype1"] = atmat["atom_pair"].str.split("_").str[0] atmat["attype2"] = atmat["atom_pair"].str.split("_").str[1] atmat["cutoff"] = np.array(cutoffs) + atmat["bins"] = bins + atmat["kde"] = kde + atmat["tot_repeats"] = tot_reps + atmat["sum_probs"] = sum_probs + # where O-H set distance to 0.195 + # atmat.loc[atmat["atom_pair"]=="O_H", "exp_aver"] = 0.195 + # atmat.loc[atmat["atom_pair"]=="O_N", "exp_aver"] = 0.29 + # atmat.loc[atmat["atom_pair"]=="C_N", "exp_aver"] = 0.41 + ps = [] + for i in range(len(atmat)): + ps.append((np.sum(kde[i]*np.diff(bins[i])[0]/bins[i]**2/ tot_reps[i] )))#* cutoffs[i]**3 )) ) + # atmat["probability"] = np.array(ps) return atmat @@ -177,7 +200,7 @@ def define_p_threshold(self): if (self.atmat.loc[self.atmat["atom_pair"]=="NL_NL", "energy"].values[0] < 0): # self.atmat.loc[self.atmat["atom_pair"]=="OM_OM", "energy"].values[0] < 0 and # self.atmat.loc[self.atmat["atom_pair"]=="O_O", "energy"].values[0] < 0 and - # self.atmat.loc[self.atmat["atom_pair"]=="O_OM", "energy"].values[0] < 0) + # self.atmat.loc[self.atmat["atom_pair"]=="OA_OA", "energy"].values[0] < 0): break pth += dp # find also the maximum pth such that the following atom-pairs have energy > 0 @@ -227,6 +250,10 @@ def define_special_nonlocal_dict(self): def nonlocal_matrix(self): self.atmat["energy"] = np.array(self.regall(self.atmat["probability"].to_numpy(), 1, self.emax, self.pth)) + self.atmat.loc[self.atmat["atom_pair"]=="O_H", "energy"] = 0.45 + self.atmat.loc[self.atmat["atom_pair"]=="O_N", "energy"] = 0.45 + self.atmat.loc[self.atmat["atom_pair"]=="O_C", "energy"] = 0.45 + self.special_nonlocal_dict = self.define_special_nonlocal_dict() diff --git a/src/multiego/type_definitions.py b/src/multiego/type_definitions.py index b66af1f5..fa3d9964 100644 --- a/src/multiego/type_definitions.py +++ b/src/multiego/type_definitions.py @@ -332,9 +332,12 @@ # from .atdhist_matrix import special_nonlocal as special_non_local from .interaction_matrix import InteractionMatrix -EMAX = 0.15 # maximum epsilon value for the colorbar in the interaction matrix plot -P_TH = 0.15 # if P_TH is None it will be chosen in the InteractionMatrix class to have NL-NL repulsive -matrix = InteractionMatrix(emax = EMAX, pth=P_TH) +EMAX = 0.15 #maximum epsilon value for the colorbar in the interaction matrix plot +P_TH = None#1.0#0.3 # 0.005 # if P_TH is None it will be chosen in the InteractionMatrix class to have NL-NL repulsive +SHOW = False +# PKL = "atdhisto_density.pkl" +PKL = "atdhisto_mindist.pkl" +matrix = InteractionMatrix(pkl_file=PKL, emax = EMAX, pth=P_TH, show=SHOW) special_non_local = matrix.special_nonlocal_dict # print(len(special_non_local_A), len(special_non_local)) # print(special_non_local_A[0]) diff --git a/tools/make_mat/make_mat.py b/tools/make_mat/make_mat.py index a2dfab9a..c71c3146 100644 --- a/tools/make_mat/make_mat.py +++ b/tools/make_mat/make_mat.py @@ -495,17 +495,33 @@ def generate_c12_values(df, types, combinations, molecule_type): # ].values[0] # print(f"c12 for {a}-{b}: {aaa}") - # TODO vectorize it as in lj.py? - c12_matrix = np.array( + print("Mapping c12 values for all atom-type combinations...") + # c12_matrix = np.array( + # [ + # type_definitions._c12_df.loc[ + # ((type_definitions._c12_df["atp1"] == a) & (type_definitions._c12_df["atp2"] == b)) + # | ((type_definitions._c12_df["atp1"] == b) & (type_definitions._c12_df["atp2"] == a)), + # "c12", + # ].values[0] + # for a, b in all_combinations + # ] + # ).reshape(len(mego_types_i), len(mego_types_j)) + # vectorized version + c12_lookup = pd.concat( [ - type_definitions._c12_df.loc[ - ((type_definitions._c12_df["atp1"] == a) & (type_definitions._c12_df["atp2"] == b)) - | ((type_definitions._c12_df["atp1"] == b) & (type_definitions._c12_df["atp2"] == a)), - "c12", - ].values[0] - for a, b in all_combinations + type_definitions._c12_df.set_index( + type_definitions._c12_df["atp1"] + "_" + type_definitions._c12_df["atp2"] + )["c12"], + type_definitions._c12_df.set_index( + type_definitions._c12_df["atp2"] + "_" + type_definitions._c12_df["atp1"] + )["c12"], ] - ).reshape(len(mego_types_i), len(mego_types_j)) + ) + # if a pair appears in both directions with different values, keep one consistently + c12_lookup = c12_lookup[~c12_lookup.index.duplicated(keep="first")] + # select only the relevant pairs which are present in the current molecule: all_combinations + all_combinations_str = [f"{a}_{b}" for a, b in all_combinations] + c12_matrix = c12_lookup.reindex(all_combinations_str).to_numpy().reshape(len(mego_types_i), len(mego_types_j)) all_c12 = np.sqrt(df["c12"].to_numpy() * df["c12"].to_numpy()[:, np.newaxis]) all_c12 = c12_matrix.copy() @@ -705,7 +721,6 @@ def main_routine(mol_i, mol_j, topology_mego, topology_ref, molecules_name, pref c12_lookup.update(zip(custom_c12_dict.name, custom_c12_dict.rc_c12)) topology_df_i["c12"] = topology_df_i["mego_type"].map(c12_lookup) - print(topology_df_i["mego_type"]) # preparing topology of molecule j topology_df_j["ref_ai"] = protein_ref_indices_j topology_df_j["ref_type"] = [a.name for a in protein_ref_j] From b01dfe365cf6c4112570eba3fdb64e11d42d251e Mon Sep 17 00:00:00 2001 From: Bruno Stegani Date: Wed, 22 Jul 2026 10:19:16 +0200 Subject: [PATCH 3/3] good Rg with hbond 0.45. Missing CH-CH interaction --- src/multiego/atdhisto.pkl | Bin 645978 -> 589072 bytes src/multiego/interaction_matrix.py | 2 ++ src/multiego/type_definitions.py | 4 ++-- 3 files changed, 4 insertions(+), 2 deletions(-) diff --git a/src/multiego/atdhisto.pkl b/src/multiego/atdhisto.pkl index 70c9b55de23e300979c2ef42307d1ecb4e65a06b..ee4c895747240d32a83c9731a31a1f464e8f5483 100644 GIT binary patch delta 62781 zcmZ5ocOaGT7q>?ud+)vXp6{V#?@@LUqHHRo8!FjEiL|sdNXXStR+3Ug%KBOvA$v=H z?={Nr`oqTgJwyF`_4n^2Q3i(|cMS;gKOCI7A9-SP)9+7l zUB(y#C;z>njO^vn&pMgsGp#a{C<9!BjszSz6oUSl&La+i8xKu~GQ=3U$hUqjNsYcr zX>WA&z)50nf&DrI?C-LOhCI9l_k*6mSOE3PV{ajE5pQNaiWgFbzC3PSn_|5s_UqiE z726xoQSB|cf5U$Co43?oG&G4!&p-~F9!JhilHsB28*(Bxb*d0H>#4nQ!HH zA_1xg>HlMxY#$Q5t!I`!^z0H<=RZBSKiiZ~zV61;vxQ>0gov|I`%Q>B@#T(JLR*aM zk!`t|7idL-ns>;(>H8~ECg={H+$}UTi5zD?Z5k0d&i>vKGL=FSdn1pveS|pke>Ail z{z<&yCPBe#>hd~G+iI>QDOco3^!_vdM&55m@oNa4nk{rwO$l&+CN6Q?B}!N#Uw)m0 z!$@$*k`zyk-upCWjkl#4ldsqWu2n@JD2P3vxy$k2(wd3`$0{!~F^cvcS!85{pqsxPP6Bqi_!ZlN0j|9Zk3QmQF|Y$^YVdI7=QMhvLMKJ zkAAYnh``T=VSz_O;vnEbX;2z04S)7)S{|j6gO2c0P3c~F7}YH4Rxefr3VpW&_s%Ip zuj|j8EE81>M9N$|ahg^gmhXjMVX4!A8Pbxgnf_XkSC;-(f42@eD+ukkHAi4&{N=(C zMgaOf-|WS|0{lFaKWf&U`gS=qvQY_DHYRW4vz2uvy>Tl+6nyDzhBPB8&DJx%}i1p0&RJoD+HN`fbQllGGp$;FlO5mR&?($I2AI! zb(D06ydMV->?v{wvafx@7E&ItFe#Pap6>x}5t^06+(*FYNsxAb;t_yLKb*R!j{sHf z?XEa8Pxxq-tcy{<=?Q`)BsZ6bJ%L~5WP+2R7jT_0HyyR{0{$04(JLps!2OTD>&X-^ zc==&PJmQKMxSSMIeVpY5MSq?Ywx@XkR+0Ch^BFIoQ_b(rKjsCVIYpU;mR_KD>PuOI zsu#G5B#$UrK%JpkV&kLOhrS_2{7KJ`U?6WA9duQpnv2c5Ud<%`a#!I2EL+p~UB;K@mSj=|(R z2bkU8m#dgfSTdc)MwHPcLS|{i`TFJth4xD>JDa5z;V4w)iBx2J;K6NfxOM2y9&FMM z3d-T4elR5*-t!=E%0$xHkfy<2) z4;ZI9@r7=pIF18@zv)U9tMHjj8zs)qR>+2rcBHEXBChvg-CbE}ttfF`8mvz?6C-%xDI@+n7qTVduY0YqjFAbj`?Y+aFyntskTg?8vm49V%Jx z2$%4*Y@t}ob7aIg1r!)Emp6gFH_H&F{Gs!e$ZJHiGPHrZxCSW`ieV*w`W~S&R&EmU zY(jbnA3U0P+KLQy28U90wj=L4qIkY#^dPclFV8FM4Iq*Q2XBT+4YrE!3o6sVUD5b!Hw|fFKUq*| zhCdy28zmJ^rO?CCyT3(t=P*EDEzer$bw+5`O)?8U&jd%~9{;AYXNDthMsl^7SRi20 zDD!9;h6Tv8u!CiatUxC!tumj@3J*h@ZN!FIK~?-rM!hf_R1JPJWHDibruSUTmELTi z8N_Wl8^Q*4jV}@+{n((w?ajlLqik@GTqt+ZjSW;e6y4Y@*}(AE>D?3(Y!G+Tk@6WG z8|0^$(^oCC!fZk@`${V-jJz{;D#E;A1?2qmB;H(BxRo?vla{~=W)=>5_oG-Lk^fx( z+heS7GGdRqVE`*w-sre85XK6>BIESh!&yOO&~xZeA}cs+m`-WkWCc-8kt>p~SRtb^ z)M03r6_zHs?%9d6fp1;IgL#~xWVb}JYS^I2?lB3CDmxq`Tn)pt*Rw<3gU7-lu^b>e zo?y|X#|b^={MFrGIbkyUsxbjE7m!;hoZnT+3D$axJyFD(lVB_Z0P% zMSr3IJ2SkyLy7IDVF2>?tePf)d%&DR){l;vET%BL==+ zGggOm7YIKaW#zqj4b! z2srrq1QUXwbk}ZD!PFpNjP}xRy%YqR%y$FAZv+7sXW($uqab+9@-a{7We}LtR321n z41(*m^MVuuL2#iz(T!#;2$E|PyhhoAp~C6J^>&S5=nZ1{$#f_fRIPJ0KgS0H@pRrg z#-}tGvWQkShWgN7V1I-T(}jR&z>>`jgg}p^D(8pL5Ma!&JD6P(0y5(RE3E{fkZ_)a z7C9UWmTTgU+U=oGmuXKoYa0d`XHHe7Jr095AvM{r{*mYV5DLj> z#7g!)LFJYzYH2YC!=cqClU&|ifxv^&y;5v^9Ebu-+9L$KVU4VhzWVoJ=og&b?_uEt z?oX*XgxjovPfhYu*F8fpiPr1?4nF8uVyQ@nAb*LWd}E$)dN! z#QENxy}0CQNa=O+nE3j3#Bd)tF3ov&NAAzmP!QFjWW|q|O%yvtLX=C%VPu(id?OP7 zC1qs&S$!F6xk=U|8X`9Si!p6;Q5XKh6S9SdCXt$1;c{giF;?IK`6o0_p1n)58QDCq zoPSmd;3E5(;{J=1m%e4uJ|So0x`@o*I57VfU*0B)ZDu99N{nVOT1>l@RGbd)XXFEw zBIBddFUWkg&tX~bHsrZ|{-55)4uno<<_|GdFESXM%61nSK=vLNofsh;M)a^hxV~Q< zL+ba2CQIf`A)d+nmahfp5NjWgZW-4lq((>hi1^(##7x?^=V~`0h>TjzJWM5piF1}U zjGrliYbf-`#XK4e&?f|VeVV5Oftv*j<~@v{P&e>&b&3W4n0b&x(4fpXJl&qG+mUze2A*3 z`=K&$2r)Scqpvp8JU`@6TyDZUc~83%oZG9q#>AxrZu9!)7nKy@GQ-Rw@y)1;jBP<)|mWF{BUlT{=B;i+if$jxIaY$Mr`TZ1gT@;*o zq(x;ZMWE#j1iD=ig1j>!pRdaa0;~P^Cfap=;9LAtMYqNW(J#)`i3{;TK+CYVa|ACq z(Yq|UAW;JxxK3djZj+K36Yj zJ&yR0(vdki_aYZ!xH4~D_>65{rFcrbS&AF9I&Lh-u0#IkZ(d2qEyoqc{oPcdSGU8) zdAh(oLhQrz1*q^qCy%KMB)q}i6r@Ym`W zrd};f>C|=LEQ>2Fr3aXt}WmOu~v4^wJ;u7iWGbNiIx{$3 zHGj)y;sXmVZL0V3d@vHfaM;_D5A*`^&U^^t1EI-^rn`B3kld-!BG`()x`;iBBjtzn zvKb3S9e%hyZK=E;!w&-`-nHFt`JwuX%8epg0g$~Lr}+$XSO6#{Vq1@v3qXi@zS4dU zLGaW~9bbzTgoJk;%H<1!pmmepx6oe*7B2`UB~1wd(ad-0mJne$&}_Je6AeQ@*Sg4{ z+af?XW8*YnFAAZl8jBz;26RM3dK8@EU}Qnt#mOxJa~GbSo0yk?2S2L23fv_D8Pihp zb?}yiyr-YPPL4>xxM=8bA65t%$-6i2>ELbHWflz$-zLAOQRNY;*7Hp&rVa ze!KF84fcHt&}{K#fZ0lm--8#afVTKDSsghM2phkDUd}y_7-!Xgm8<@N$h%i673F@| zS?5VET5e{8e|6sc?3wk!a^>wfX}n$z2`_hKt#pBph5clE@R#%^icN?kzl!dJztfU# z$n6t~qrTfQwUK4^2FeZe_-)Irz8&&4TW5#dw1(A`UYT=va<@=y9}SVvX8jm{(5d}` z?l@w`q+ME+H;XXtXD2Z2T|v}tz*++j5j-8wtIGL64!cTQCM~*X;Fpq7RcQeOtbWLP z$a{?y{KD@h-*x8(YVD6rKInQ7zWvhWy1XzrSJltHFyX-L6bwO)ylutJc(@4ILw8oxP6`_%LIpBs&J6 zBMVA-OnMMtdC#UwY%j!hxl`zi?}KjZ{5uzE^kJ|6k<5ioefVsAR=*={Kb-0-j*pNq zfQu^&%F?9`KoG;X_q-uo<7s(n){Zd*t2)b=IUXZ8ddx9w#mxxBZtBdI zT{VIOJ^B;uO-8^HF1&Byw-E%A|N2fWX$+5F^-&zNGlpx;U#_=B7{fEZ=Otpdj6sL% zzK3s>F}wS~L~=}EX7ydQafJySRlnp_JZJ*q>5YQMw5Cu`On=p0+Y~(JPem&Em;x28+LseI zOo3nQKQXFSQ!wNWXD-`i1`GNH%HIsk;J`Xpe{qr-BrBTjJ>O&o$8~DT?O4nqOuWlC z#LXNUIK@2u3e3Um@FPM@--0<@yUF*JiN^xX-D1gPx)ujgg zLu-~>7Qok{|Fhz?1-zw=Idi`seRW7D`OaVormrOhrL8O>L(hpY>ZT>gTch<;oxbFFV8M{YB1o2{rui*+g+o6KrZ6IV}lxe+BtjQsGPeTeA#rej}w4Bm7vMjpo;7 z(XWH_*X8ZUJPC}#kS(9*+MrOKL*=feG9@8wLKvu_NY^GN6s_j z6zv!v*a^0#`t^)%EBbuk_{>1^j`2ZiNI`NjHzU5HO%$!f_5OC|i?Dmt2>(ON2P-lb za-5bOjYtg+Pog{N6Ir{p(efeXGXHtG|TuhOo6Lh zQiAWbIR+Rzm8Z|sS;F$g6v-fMD`4MiO|$M}4ewLgwe_=X;E{?HePop_v@<4Z#f{m) z-Jt1LSJ)1KfJMYzIpc%ysRui#5bFR)&C@4ywT|$#r-NCR#u?UHik#Ct~~~? zUBSGZ{c%Tz8?+vbwb3&>1YOn~CN)%;!@$!`y+qe|7(~CEqVLFXhY%giOq#O?oE%hh zaTPrRUlaObITw$BU28;t(K}C2)DZd|l>u07`M{@s zdsP``U+_P0Civ~TFX*z6A;S~Lpfg@r`$f-jIKJMm5Yp-g7t~e05@0_016^epC&nWH zK9R)zQtc0bpL^2WWDf@dTlS-0HT{8b-QDbslSdHTP4x8{MmO|zy9wh}jt0ZYI=PUy zsL!Q;yFZEX4}s%9pE0)0ga9Xdm&wu4P++d@ylgNZ3eg;22hRtGL2>j`+L@^^(3>gi z8Sps)>(|m`G5P%`KuL^R+UCGXSTou4^3I2oXfbDTr9$Hr-2G*KOzXiZATXS=F{2NM zHBP0z@RQ*XaHfg+Kzlemvb1$9RgQp8#~*t|(jq{|#?8&=M+CInr3SoJh=i{X_n&`* ziG-Dt-Lq-+kzlZnjgn%Df`SMS{mXt)kSJk+2{!7E0zs~+vhri6A!bHM;UgagSf^Xo zs)jI7!hh6HQ}YZ|PZIZ1nV*4z->*{n#Lobq8}-|v(-^q^K(JizcNEm+x*o5uihxQ# zx>owmlTcz@^_0sm6zYFcQ2IX$glW^zn2*1Ffl1`t-J?Ao@WAk(O0>B%*bQ;555AYJr-ZidJ#j&o#e zT*-B9W9v42ctW?(lBdV7qMLu84-p$pg}&s*lE~~QIJ$$(3zs^d9n0`N!N!-Izx>`t z0e^R~g`#$u;}~a{7%QbqScvws$3vAi`tH5w;vc_CqQ*&6CRtSD)#KiN=Xug!j@44_ zXzC@oiq<=sKEN|_3&r}7q2mcIlr6{A>G9dWiGlIsiQj%%lClVP9=!-HUkVAk3dX44E53kX9eNh2dfg%dSMI<-m9@VJueD(V+AS|=_Ej&P3N@> zl@yE*+B$3fl!nJ~#QRPg%R+q5{rGtrIf%Tu>%hVbIVft4&3kAm4~9}2zX+SrDJ!G7 z5Qnw`^sG=1GiNKnE$^N+uMq|Koou`$Ev^V(WJn|&k0^qhPt%C%B}HHh9!t!uKpPcy zg)LsqgNpE3wEbB%nG$rGeWxZAP=X-t)9N%jN)YuRpwHD>3Hlql%AJlVfnJc{GI^j9 zs5zZBE(%wIg06Qa9`Q<`lxFLcouveXLQ6t91xnzZo=$4;RtetSn6hIYPy(t(9wk$9 zWeC{wi_la~8F=PS)U9D$l_5ww!0<_uGL-AnycjK0hF!NB^P75=fpS4)-*Yk*X!cx+ zi;z}<>BiCdKlUnM%s*GBgHeI~Cx?qh@>O7xV$byM78Ot%e%~)atP09;ylQIls&MC; z_R>=iRnReY>ghWli#0SM38`FgsndjvYfmt{ z+61&90^W+e6V!rJmuJYU-)llEz5L_uB@Gz8Gkyk81y>GPOhjBzg>Lny%Cf(dz(cpk z0n;xJ2bCBDPD#svL_WW6Mv^FWN(JmXV#5zPR7S_{maswR=k}9t3Tfe-`iJxgR~4vYb@&msiJ}pPMlIS9nIgv0jF8+TwKEF(qqWHbNxjeMJMZEmcC_9n zYEYzPTLt@M-JQJH9sH7b&*0g;bFug;Hqng&2aApmWk@gzR-@;#53<5bVtrLbv~A#O zk!{PfA_P{Z%N+BAl8_@pnIcpn4~i99gT0PYTJ?65=(hhwP zWC$!vzhn$+W-e*Dz7`;P`D{ysku~@btAz@?*#TQmS;*m=2O*1>m!yW!88ELGxvJIO zFmP96pgnN$5L~_L+NpvLj$1wodrikP z&Rh8cnG}nkP3keY9sB((@w4MFO@`Ai?f6W*kd?I}U02tj$u$=`y&`D8qo zBSA2XehwPb%z1tiORR96i&b}p?&+8gcA^PME4ynTEy@-J5@8PpM*VQjxDd# zP6C%&M-QfO@FYAAk-uD}b_%Y^W+qi8oPt3@xm`u=r{MYQJA>mw;jmZggxuWmaIn%i zHB(m_4yhs(w{zCQp**YqPPs({7#EJT&Xh$!smk1?Q!0^A+q{h9-Ui<~46zM+nq%1=fP1Q+_vU~Fm<)yq ze@fSbPjmv`Y|7+8DIvMxCoL;lq(pB_6YPvS;z#K=_BsEe)hswXZje zi*9QRA z?)M$Zy;XE4e4qBgin_$T(;aUZ_DXKTW$f^X9w}rSrQYeFyf8zbnIAvh^ zMEr^SvNEXLrByW7Lh;vf6ayD9O5nG;%S9kb5pwo?jysgTC&fW9b>3@eK_kIajEL`p0q2r1~nJtG1 zUVJ8Oq8pX=HVXt`+fEYOp@p^?ZTx0s@VC$gMq3cMjV;PDK26LfkPvAqNxH0xFLD#z z7>p?3>Ht*|Z2ev$7j&P|@Ka->d)4=$(%GuY59b=L(BUNaHJ5GI+#}??I&*agtt7ZM zcXX#Qgc_O{(kbg}knxl@B0rY>i6|wdoEE{m!VSjvif^E*CUdQ#C38o86#t+)nVq~BU+yM~ zjr>iPO=)Me-ywcO?v7^rYB>9z849K>+M zE+Rh=$@^vZ_s5MPN%S{TzJ8fO-Up@Dq}VSbLFau4eytLK{9vW?vJn|9dhz>-UfB(Q z7*FUi4$(npllwXO7A9aA{p^(3$A$sj(<~;HU%7xU=r<8FIv#qVj3G`pocb5spDMK^VovI>eyU8lEplRoaByF}wqEGoi)~e)1 zz#(3Epso6iVmCKh&v|vpoZx_*7@|@=ayDo%iu7^%!3^Q&iQU^-FibFMb;d~TF#{OP zYUa<~q=)~+2~16Y&_S?P^t-7LI{0~HR*uM=4)~hp6`I27VANPe!EKlhym#~1c*oL% zB~8%>4FmN5`mI#YUhE5I`(!E}N-Y=*dR6BA(P6689Zv)`)$b?z?_h+Mud6~%-T|a*Ak(Cb;=89S{4O)L(W;HUA@uc(X?XZ>xv#<6#rU&1wH0IHbiUX;XNjw0Y*S?MBq~{2QH3 zr3XD{Z57X=+e^X5b)j;Ay6hD`rEQ_u2qW}#hn>;xDgVvG{Qu_1LQKuMll%T=5ysZK z15CU-82&c4@Htz7eLNJgHvLy~$@w-C_;zlh8^!eBhLtR~tsPB!wnjuEVH@uFS4|`H#~&H{ z=+kbygUmWg_kzKW*+iUjUq@Sq2ELk2G&G6W8rN1NCBZ%#Crd#!Fq08)F1gK%CAbH8 z3-bhfEB?{81_vVN{n-{f9WFF?V{|bUU&AJf61~|*Ie-QadZVGq1$U!?+)<8oW9#99 z-4{v4fb+hk;|u-#L!UE}-P!oh&sX*!hw(Gg7P^@^ZbgxQBEr1veg7gLsT>(T#+vJC z_Xcs7^83J=S%Va!?OHy=dStI!c@AT0BXWL3;Kq^ZW~AV%MS2={8?wN>oL~bt0q%r=WezY0?ZJ8RUDT`z$jwBx99DzOez#e#@p` zD1BLbjr%>B;Q`*?k{!)5d?0*{HK&G?AG~Eg&4(28gVqq$r-d1QAdftElWUM4n4Sii zzl`IDg+nzjuFUhnd6#_dI88n{y!t|*&<{P)AitNPUs8q#KHo83Tg>2s+RXBEBAlFX zcagKvp_Lt&z-R1SBO9D~a^`hjH!GB^*t`4kv%-py1i>H zR4+5Zx0#;kB=j(@qeWB7ohT-FaNe!4DS-({-E==6%wd8o>znU>V42|K;d(Yq7ZXGu znOAh%%?yiHIlU;a!8v64rpus+}i6h2gX@| z$%k*2LzfjUXQnn{2+=dL1_k`qMD46VRw3;oDai)U56x6~o?wI3q~eoCZ_w0J{#bQ_ zj2#{-cFtBgvjbV_GH+*5eNw2iOR@$-YYG0OKRCQjVdXdlV7b*iy?0*Lj6S zsgH8O>ksB#9^cScF&`iPAuebdWSBLIaT(&n&&6I0qcE)U;5gWCJ;d z8WYlk%;hZlO5BMVUOG7|Z{QXaeP|a#NlSTI}HrjZSjP!38EuS+rhN zZbp{OEP4{~{DBHzw?zB{6k8~cv#fCA)JbBjp%U7z<1orTw7N;(*5Js<--we(I@<6{ zf7m07MuUWn>%3*(a+U@IzIB@@wkwPz2jyxd!Z#bWa_`@o78pDRunr&oQj%Ul{63%m zcGR@zl|Qh2CXB+DyNTi`$qJixi*Se-laVh$;r{z8k{oHBm>t@MaL(>McFe0EvAunR zkauPn5v1F{2*(a6|#7^*ugVH*Zl@TSEYCsI8QgnNdCz68_ z^-B?5kKOQk$-l^MkQNU5&?FO#F#!Gbp6?~%=t<;-g-F&L>_9H*0)F6_#wzcUoY#PAf%|?K_c9Qp-b?dg&=__SbB8SX2yy^x^<~ja=$owW?hoW zfm;$HhKc%WWTk*-GAVbAUm6U@ifKQrNJGP+^&8ezGSCzK&ia*-EQs6BSv$qZ!h5Ms zn5>e8`~GJsV^%P-us}PP^+!Su#DzV$g{|daIw+#5_JkZfZz>MyPLqTDKt|56M{=+# zoLztDr5t2B&5fpgmIIgTU!^U6$UzRn1xhPIdAJ)f;`~Mc?T%hQ(a)_f52`dNh2_WO zK~m`8vTv3=5QJyaIDe3b@dmGSHX;RJpHx&G!>B2MCfhuDZKwhW=4zEFV7(&qsY z>k7d6ZugTPCW>Gb>2+Nv2c0>v-aGblS`h-@rn89}C_&qC8QaWTN+5eFEGT*v9T`mB zAIr8uCjmwuZoer|hUej8M(3$j&@(doj+{HK0+AJCQkLr~aK`x&Pi3wuNJtqpV-juD z;Pk~*(*O9>!Heglr%I(da4gT}9V5{I47O0^#E?2rOJ!j-&DG(aFgxQ9UNvZLuc#)p zQ2~>9=X=3Zs6pI$#W~FK@H8#x{OfrsAW5NW`L9R>9!vd733!>n1TJvM zvvdmyC;r~R{O?UR5BJQDjp(_TEm!LA^x&=eEfiOuIQ*?APGV_TNtn@WUtYO+o9N%+ z0OTinp>*=K&33jA5;)6zwR%SxbncSG^5F;R_-Zy$G_>?^-E=RK5@~N~lmMa)ZTaz< zjCRGghPx?tWkeij<2ol2xks56zgF8qvAGJQ*U_CDQc}1j+b}8jwifpR*TMi-} zyV|NoWp8h*sAl@akzBMx#o6VguFQ;ZJQZ6gTF#>sgBDhU#Mrw;6osgYHa?0C6+Jr+ zMGHqyd>%{OR?%-rO=du`gT>B%iavYx33_br|L)Xop*ZdaoUh-Y!;wb9HdMv8btMv< zvwP+ z)kr_o-AF@ye|HhUBGX*_E8imvhugo0>fX{}WX_R~zh;Q*GiYq|@w_6J%#y+*A;3nfR zqT9D?7#Xt?;#MqS)6wXp6|! zNPKYfELQIXgI{3d?Tv6c@QPYxrFL7c$IW-mez$f!zt~yBU8(OsEx_{#{it;2IH*h+x}Rvi^|VIriI*eo^et;p4>Lp5ZImL~;CuK2rk{ ztbIMP8Tw3P^?kOM$lqBbi+-&Kk+RXYg4Yp~D&;##Ayz*&qdWW{zJg5@r@&!*t>z)sg|s;ZFcKTf>5WwcsQPB`8=bwR(HWSkmncE-tJ1Pm3-oHm z;F`O{yfS*@X{k@}iVEESKKs*4TNU0*CS8vkRR!sziH_tMY9Pm|ENri%4kpgNL}y0T z;f_INMFr-P2GFfNVocH31VxgXNZUnCxRuhh&U8l$4BEaw^RU;36q4WGX=FMe6*@{D z{z3=7Jm9eJQ$#?4VIt;xG6K6BZq7FjA)xt-PvC<#K%>*s$HJEYL?1fMb&Ue_Tw^@D z47!jt7{u$DtqVWQb#L=b>cZLAo(>pALp`W|{vZdrtq1$V3Lf5H*Mq?F@cS+UdZ63ZP5_j2cLCk-89noL40y<&dB;c29A4 ze8)Tr3ycL!VE$??$F3w}P$K*@p=4zQYZ0{_>H!95Ei2*&Q~TiQ=qnz=IbFa`wUc?d z>cF?l*TYBOszU6ujD7sc^3YapVe7yp4&U-AF3O|F_Tvw*+lIEF%{{n&m9G&(qoa_pV59tDvCdsF|_Yk$_`LouG!2Zy5;eAv@{nn$e(NCA&Aa0!rEY4rDVQ} zWXbxbB|5Dkv#hkodycFlE1$_C@0Sz6g^b!Yl^Q}gopX|BsgekoxU!nu3Wz};SKG$k`N9mG|9uyjBw>ZQk`&>SN^B4kCDj`5$qv;E6Hg(Z z18DU{EZN#QA>aF>h!QzBR5`S37-AH8z#_wf{gN~ zWc_^bj&r87!XIr%-4^9gCl!F&Pa(!?R|MeqSn?AsF+n(miTzIT936aI;`@(LR|xd_ z;o-;kLU75tYutIiFm&qKzKMPz42MW!&gaXD0L|yK9GuxAu($j=A<&UbDGjVlnVH?6!P0 ziT)i~tx8RgC{#VVbSJn?6dDE}&lk6f!sK(!r8)GKwa57BjR&GY6xI7STwWBUjv8g2 z3m1WdzT-{=yS;?rhx#dfp;WTw%t*;#iZOKX`ErgkE2Z#J!kTy9(5(W8yAemwtR0fVn$6D_%Q0!IMt4{5nKwh#ZB zKbEpW*5OW`+hWgO^vZwmaz`l=V{9#3n)nA_>?Vp5fH{;9HlrJq4|C`Sg~RBFy~_p# z_sAaf;PEiTQ4CX#_v7Pf9a)tt=vnj_t)I$8v4`=+ZKAk3(*f7EAtS+5kb(Afk&non z^Y8BOEUidf{mCEhBAtjZ!X!kj-H&iwcXifg9YGvxA1P8En?&>{W}Hc=<`CnIhUG=3 zWh7RF@fg zix=q8dot+OcYzTDm3--UM+TUIcd5{ks+tuPAIBOwC9p%1wa<%8O-{Hypv`OE#07_F zHLnnA@Bksz(O;`(yl_jtgc&o-3sN1@LtUwS(DICSu#1%+YVKbma=pV3r*B5sG4Tk% ztmt*x%gF*j8pkGBxheoQBZe!~K7!D@-W1c;j1dHek*SOFazfx)vpRe^MF>7DkDX6I zZ;UA$zndR669(%v`L;LciH#}kI?Z!rBEXo_y#M%<3h7c%&vLaOY)J~lzmnW@v6Kcu)jUS*Q)!6G9Mh=f zl!1OWubUSx%RmHmMF_@nk1S|P*EAA)$inN-I+LnpvJe(pe!}9HEQC+H8h8)Og2K&G znf`mSpt_HbkWmo5D7C-T|9g@&FpClt3<^lX&C=Ydb(CYYCmAzL&juKZ;ea2hr}qUuX9Tsg=`#m@?}q37}ovWbfY&CGB@b{zXrA?6G24FPphu!2g(&Efjm) zk`(vpFJaL-8*;k|Pes<;Si{Nfd-wj<_bG~PxoNu5;l;u`#*WcrDh`BElz1_0q~I3v*i(HX!}Ja^EjRwCJL^7TykE_& z>*_-!&bP@?y}kgUnJeg#<0(R>n?^j7RGuQa2Yarph7}_Na^Gifgq0wL4Z0mys!I`5 zf_^L)ZyA#2!nBXL6N_9}a-ewr;W+cgsV0zM!S zXC!<%%RVB@eaC7|$v+`aPiYyIS${^fo)`7TA80~6qsYn2Mw*bmB=WG+$<4^uAPqfs z0n8Vq_lU^Npo%ZZ3o@3$gGnuj^nfUPQ&cPR`%3JCFR@>dgw?!J!Ll|)k|QhW$AXr)5v{F=3)mTTDnJ9kgp4QP{Ab1{<<4s*ElLxV%LM5EoQjaf3g>;*>~fK zo_-%PCm0xYkPq!KI~810UF}C687Wu$V7dm71NI^6$lV`EfQM4`#N|O`NcwWpqVo`v zv-?lwg4QsiY^a%=N`=<7FO@|_2S$)%J{ZCuw|^o%?gB6OT8|>?kA^askxn4?{td+#p`Z!G z>t$4l=j;SRvb^g@qroIXCSc;o(LIUWiwKg46_`TkgN;p-{HBmob&t}s-cBKx2F{%j zk(frVjiUE5&P*eH<9b2k1JlUQ!dOW^%^BqNnfJ%3vuBXDngq$vl^NudtgrYB=UF6k zIyTVc-Yk-pKQmX?GmB*WlEXL>G0!1VB{{2Hd*={GnZEsOQFF-AipK8ud2fjM*3YV;QN=|8y}&WCf`=@s7AIy@JrmzFeSJ zT|ou;=8S8vAmQ&quKU`rAfsnK8XrEsf_VCubWxYCAZ!iK-loW|BC;8V)5ZO(2*rQz zU02$EBO|{@B9H!FL!t-|4^mG4L4wb$37$T;j`&|q*M7m9y^ef}v38~7SVt~Bw7YP6 zU=1M{oIQ(tT}5_%3XEW)T}Bl4O;;{GnnxO!eWtjlrV;JIpB&F&9JwRkTQo&Gh{QPE zOsI(OM(Rg8lRAgKAoVv~Usa^OMW4N4Oz&I?awk*vO;_a=99upm+`OQ_6f3yFY6CJ; z7c#Klaje$rBQz_d^S1|5aa#5Y_S=raBYL2pf3kEPpVc-|Y@-f(XvxxuJ)3Wnh4*|g z=QiGlll$fNp8&=4LECJ}uh7E(rM>0fA&uxNlGP{I-0Npn!S-)fpSVBYk@-(dn$FRgbqWR-qf z?yt4QL%g{=c6QMMBrhXF#_;8CqFC{n&I4) zZ58}AuX3xyb|kzrbJMT%eIoEwY@ygsOr&3N36BD84Rwb*IA2FjXP42=z05;IE}gpl zU&uY={O{==TvVhqm7UQ*X2>Wbq9US% zY(gqZcDCP3>izNc`|EP>e7$hb>zwmAkMlS-HG?8q=xOl}el6iA$oig^Q%x8a$Yz=7 ztW&p4}CwIGozyi@;0tOYF>5BA?eNx%{I0%e)4(62T&LOZWrgrEJRbKHP>h zi3QYY&AQN}y5Np|OZ`aBBv!~)bcldvqt57Y=Z>IDy4=s1Ek{w_1NtLWXf zd-Z~5L{Mb@(!-!>>pdR&FKE>}KgI*I*U@hT0v>!sd|qAB!Go(RQ_3v{JUIPGjXpGh zgTr)g$dHhUgNN4Cc2+lV;LXfi(H4XQ&z?`g`Cd54e!4(0c@77SB_u(E*KimwuG2{< z2?x4y&T+BXIN0ZKLS*D64nn(TjNE!~V5f9lQo0uhl!uJJmh|Hwtmf?TtdBUD9Qe>5 z^8p9P#SY5XzQ!(Gk<@&-3<`_8LLF&u^v}2P7i70?n_mb}k`d zh^zw#{qck|M=|F28qL^DxEvk`UwA_4obd2G^d+D5RXkXV5WNb*e)E<)E6e#Rq>vNY z6+B!>3hC{&FK~Ov;O42i6wWa+us@+z&vKL;q*wPwCd8A2kG`0vSMfP=Aj=A$XB;Ji z#6GJxL4BkkDOGv3*8~rP`zKzADU!fIpaNT<(;Bj;doc0v_X5gWyR1GIGJ_rp5-lG!=gCQ>qD-l=Yu5itWhwD$<4}a=YOJinB z#SO3BS5%YG@CX0zHv9D9h@I{~monc&7)#%|=_wOjc=(qJlkrxU+OU?4M`KS#TjLWq zJvS(=|5c?xo;s)SaK`c9xjaLLuEi0XAphW4n*)xzTg6=eTGzoZ>-BoC2}j7>;8U_6Y(-DG>6?z7>q1o~#BPRI@rVUUV(T-Dgm96sp&uEdsTApo_H4aU!p z3c_OA+4i77VMr&AVjCvj0m8R;UVU^-6e?NFH*-2dR_P&X$I4L-Bu*BB+?sG<-? z!_-!8t$rcQWv6 z`0_JyYpfdDu4l;jR#XW8wSBN#UKXB*_-gwo%7RYh8iAmSRaffgGX@D5 zfuL#jc1Sq(1#hyrac}H(dVc5=Qj-N<+up~@8nPJ3M7mO^EemdEWda&3WkGVaSoP^q zS%`kvIc#@H7WAr5QV%_pg~5ud?xhCo^X=EYadWb8Ib`r$q>vokBw=V!-75#>GwpAv zu41E!RVM@os^q|{r~f6vc0mpl_LnErDaykMK`DL8zmWi0ZCJVdz$ z+NffD0(t3|b=PAR;Oq5^C+o-*;e8(NeLOa+_|u0{-CGWJxihYVRrW9;{IFOdX8TH!GA z@*Ut~bM?*6ase<5(vx&wYphyws2gy%W?1lHf?#6q$o90+sB?ski## zhIJit)%y=~^()B;BxS{Lmw#>kO_0=(V}+w6^`V&FP}1!J!NNNhPp_Tcdgg|8!w<1; zXBw`gXo6xK1E7k0H>6kk+b^q}3bR>{%8!=FrI=Df|G9ed%2-Kk^QMF?m^sy>gve~` zB{r-jAwjsF4b89Q0B8C3+n-S!Ggcw54x<*fW4?{J?`U50@t&riKT&G!)AutSlSrdl z;x}Gv7Nwc^J73gUK&{SYA76i7L2E;0ZoH3);RuIt)Ib*==pm%wcPa(!TAECbBB6zg zdxsfX*60C;mN}hzF=7rc_wQC7j9+%FPuKJ%2LZ+dFR%|5@<3BxY0rlX{19Kq6KnNe z5G)jaGgzvKz>yPtBC~~}pz=i0#Yj^E+(}5=?=(t50#i?NmZU7)YA!ZQJt7COl<$s~ z-H`_>5B&M`LIrqZ_tbK{S`m~Act~T*l;EuLzMUW9Fexuz19ov|6;SE87~s!LPz9r^ zYH^aMsvuib6!uC>4aRpYA9pWP1NKj9;Rj^YvC$fR4dHX z2;b;i=DQQ}>MX;?3U|Uw2EG7|@7Vk&DehwslLjTL4Y;a4g?k&U zfoVg={62JF16)j>-yQvqZT;-A+9qa-a%?dBar2R(w4W`&n|dl0;yw!2Yx2f?9u z`g72hK+OG4J)z1&NPh<-TjW7!pKZyrWPLDD2ciEH3K< z`L+O9rt`qQ;?B1QX$I!yOlux*ZPM*E-Q^WDVE$mDd}q`U019njm!S~eEEk% z0cFPJ8Qs?VF*|!O#}7(-?%^7_u*EC?TGt(4Z#@9d;Ru$b?4oMy$k85PZ) zQd~`JMfEAVYo95*k#@+_*H*rtQTC{{C2zzC5>FLD`l&xqLVvFS?ZPBd^Bu!as9|z# zy~~19+)F4rHGb)Y1rc!WHsSN$i-W50em#jyGGH1?F^vq2Pe>4?eiKK&vftnIDps-s47i1%UaxX{&sO0PJ#E zlG?We17N-GbWg`&vmH@Zhx5rWz*Vu-gRtf)1i5_mH@!QApu<`>z8-`@Rr;iHPKGeZ zWFOU@q{0-mgud_hI3)tp8T=HPpG1JPPIP3S+79UWX#3>W{T-lhj!VrV5rq++7wwbI zSgYbf>&Hjf^CaH4YpRY0!&Q%@sFWNQgX&(MddKHtu%a|=u0dcDhxF;ZyPhuMFuC1= zZO>D2n0e8)Z9jzstR}@2*my_)7cpbwuT}}@PkUtfLQ4`(bj;uSP$vm~ssf2>m{_&# zG(JlaNyDjUUT>CFWuUfXtU7^T7Wle{jY)4}0a?$|y5fZ_F!Q~=%bt&A49fv=Yj(2W za--PIAyZQZ&Ma>4GbNXTD-Uc=3Nnhrz0z_+xGDlNN{8OWKNbM>W5-rx1G$0OczQBj zfEA`AYR{8OZ-ct(&23Fsj1jUD@0qa-bL8PT!`OYysD;R$^NYtHa*Y>Ua*}vBH!HRwB z=Y7Kwq=%#;PRZvY1M-&?WljZXzq_!}i95yU6z!)Mw&%*wRv2DQk_NV4oM|$`;kfbv4gC4j|8x)5R6ByvZwvsmk z9~q!k_p_atAR~~r%Q_u&WrRIB5iS>_8KH1Uk65;x5p)iZeUE5mg!?B6U8njP;gVfd z{Ky0&oX&RE&i=^=tB(rL(*I(F?)H1H>%KAq#j*aJ2cH;0uz(chVn2VdYDpp?5&L;p z53;wH7~xSt>>01x}&4v zOh75j#4YoQ2^Q|`$PE=^hQ@lM!SuV#VAM&Pw3~?qc!qt*w31k$Gl{nNyc{d|nY+00 zB(s9xV?7?K7pzbn{N?kvtE})GVzld~S>QK?_iQ1iZYnayqQnG+s{t9?%ZELztYd6R{iIN zjwR>$y+bxRci_@Rw`AU}=-FTEMt=G)^of=zF_XW>2x5&Y8%@&w|59sIckqhv-DC0p zF0-8~+o`y5ld}iDh)Mm?`)-S~|Fy0M*nR7iP8NJga|}fg=Iyh0$o~8$o5r&;d^9xM z`myei_6lG*Ajsg z``!2#Ie2)Mti`x`kP?1OqH`fK+n`A0+=)OnCeSIPmKc4&2H}}{KFUPgpwx6TSGI^B z4m@Irr1KSq3=*sKrhYMKe1~xz=cM6sgtdu_wLJ7a`kGozP$ED}Uo~CCk}CM#l#P2? zrU69;HiqE`5!@Dev9kTNHsod4$f?}X1<5h`ukhs~Y$2?~Y99<*c zJu{|YFmt_lrpp|}e)YcSwzPocvqBj}f|j7w(-4t1U4@eKtiWxEhMSGg8V>vpW~2(W2IE`WJu^+#K%RVE!j;Abj=0?_qgYZN?X!WnRFj4Xa~n`k&vdm^vw^P_vy12WZQ$~1Uw8tA4VYNr9ym`} z1C?+2%ENAJ(6nms4Q;T-G@Z|d#}lfoVL^ad;_-89NVlrV$u6}9qMa&Il;gc6~ep$oUYXLrycpKPucCx3L%LZm7Oi15jH{hW& zLaw~m2Ck-;krESZz|nk}T|LVN&UIZlrP*r(d-t+0AWmB-u9M~=6dbgL^LdR@d^xtz z%OG7TLA4h=x}PYt`0j<0%pb?XNbP_`F4jTtsT~-%*VEp=WDf@}%95rZ*$1n%XM4+e z9pF7*abvQd0|b2GMEn~RM$kas_=ech`DQ>WOa%s$CNBlqsEQ^kUC zL&EAqF7~F|?{(_|K53!EQVOs3SHRbUbnhPtyS`gTy`S*CT92yYzl0el?iu`Z<<#G? z%j1EY%FRFrRn=53&#kV!VJ$p(y~4i)*3*HW!7%T$M_tyBwIID5=yrbM-?5q& z_N7njHhbPa<}@{-t)%b34Qt7P+ql{d&uev8(|;%1jPM0FL*7QUARN&F0xL&5l0R4R zsa(AanLh|J)f(tS;?z3^)bt0@hu1%tmU6zJvz55D{PypNbKql8Xxlilp8q0cN$%!??i%w}Zd*(u*P2Ux!chPO1>BXn z&Imz+m94f-k_e2nO+5HjDGDV^{z_rr#9_Zcoro(YEY*9A;xLn-3>dy<`xJFs79Q#G zcsvS~gUr~0iqV4vdAM@6DfP670yKU*SNRH)v_4(=^+KG7B1}8d$~+xVgyGZrFBdg1 z0jRH>*LH>|0iB^&nscTSpbn`h4-A#b*rjn!uucht?oV^L*J9Hh(~&9{v$5~lP{rO5h_M<4#q|(f_c)nl|~~);3My3 zxr@O;58^|TT@ocS%UO&&fc9HrWsArIMP zA8(aik_Up*;i~Q<@(`5QzglM@50m{&!|fms18N7q(qOQZ)=uNSjN+yLqjJDOJwe&hBL`97KsC{Wf#c<$GQ6=w zeM}ybv`ujru=8x~1RW)_CuFh3Dw$qVx-ziFYCC*sOac^fgLj?QcEIB4#t%L9f-n@D zMdPuZ7j*IkoU`QFpjHe&dd+eQP7N>GE2Jq&xPVoq;=kF*z7{~b7A>+wqi|z8&)iq zW6b(LvAjgjnQR?v(mc4)Wa`c47VGg|p8!Um7!OFx-Szo*iI?e0MyxTLo!`PD>yYB) zR_EWamfTY&S&zGujAMUN_5Ym<0wnIWFOAOpJNIlP|Lca~&F7lvh&tZ<|MG(kYe{d) zABx7l+Kt^~ef!Vs-mm1ge!p+w%cW-}z6OZt__ z{!ldN31&@pT9-BnQgu}pTvszC#QB^mTW=UVUfCT^I8$nY=rwa zhEIrM;zx{yd?zAp*1w{|+l?Ntco%gH_7IR5U&*D6$9-reU5I`}Z~&3M)~D(g_>3|u zi?VpF2N7-hhn>R@hR_$$(Rl{vFNl_+No>*aD>|FXGraTOchpEZH0nz71F*gkD7=*f<{EG16#XCpeNLx7Lg+LJ$SH6eNUxH>a?bag*O;(e z{EI`IkOB^<<@;XK?8m`-?ug2eR2(=K)kv33;D9}KRV`{C9^50FS;%|wpn8eyuH*?FWjU!8rWv{I@!w(g|LG^hU%Qh0;m3cI zW2u+)KQL%)E;aCKC8jD{SIh)%SW9FLi8=9vMwyINdwz^4e2`p7f|Cp^?&KU)?xBU0 zxrDr&HYWH+YDQV|V><+u=hT_!VUj?%WHD-`Fx;AMiavHs98zK;nmCXwgu1(WUc#VZ z;z9X+T-a#AEt(wj#2-7MiRJ@Qy$ym{{^!qEI<-N$L{&;LOb@o>*&@027!tt7qtJ^} z-UQfi%%AGG%^@_U=g`R)BAvG>VnAC`wkF&XWDqz+kSXe`D5^$;sLly<&sm#m=}vIj*@G&h%L$r?&e$*!n4Dq7A^FB0Q)jT1XrK|l>I_aD z^j7!YIYWo^o8g~~E|79&m|N7r1!R}cK^?a1&dp3RJAc*%=8}TrKkB;z^XOOmj(e_< zVKqiEGVcocQT~A~_J^RGWPfR1{UL}e(a)+ebA!D-v6J{8Zop~$B0J#0VXzpa)*~>v z9sxu3Qf&yMOnypP-lkjmtI@8H^wh!e~ zCy+lf7($UkoJ7~=zoL|#TqK!aej*>M{Kt16Od-BPSrNXaIrPKturQbTGP0gK`m72O zgWQ{24qD@Qu-_H2)GSR2?mizM_Wz)T>sqeX!Kn-&M?F6ga+w*1PKpKZsKJ`U&JQM| zgVZ^pm5tf&YBvw8;tYf$?fD^vZuYK=tsqF}j=Ie{3WM;w0Womi0h7rT={rwg5HGZT z-xeSNKUCu4XPL77I)S?V9)9tv=FIIUZ0lAr*2MGMiPZi=8oUp{fNK@m=BSlw^sP=Xn|J-56MDnXhc4Z8x? z;#7`pjrxqOlgTeR#t}WD1fu%sy#i}Wkfq?IxsO>H+B9#x3*uLX(1Y>QWMax-{82{2 zPf{5^a$VZ?UONP#Mlpvt08c@M1sUTjg+_QyJ2prl3U*WzdW` zemxMg=hTU@E(j}wVZXMIowPEXa?hI|L(1^R<iY_a=mgfaK}wx{{G8uqV9`Mje&JzW_1>E}Ibu}O(cYo1YD87o}8{jCCjm=->r zdoXVyMGV1YTIEMBPaypgNBsGbmbP^e+W=>2+*9}TPeS}Z_{}7d{y9A%s z%|Ilem?@v0gLn+s5`^dqP&A!W-|FomwAyle>NERuBvt-PCQ;)BA|0}!S3CX^iB%gj z2i|{$q?~AM3QpCd4>%U{3WaybVR6ZeI`Tadq}VmE-LnO?k+TYo>$MWl-s9T%f#xk8SX@T)ux0eSGo{guaTA9wjQLD+UAmbz85*~X;<7M z-;Z$Zc^6HS29S5OhxUhVj0P(bUwzeT2o1kg|H+CG2jgw`Yt!(2MQ=rvJW0*Jqoy`> zo&2IPgrqNrmTC~jk!xF=zlr+Ybf`jv@`WR-=RDM5Yfch*l z>*Ck=&Sbp_=!*uQIjwTgB`9pQTtx`qlHs=wtBwiAKifRcxl2N5)CZ4*jMCjyA- zYHMpDf?XG0HvgO>0_S61r}xnkLs?;)yNx(8{Oq2VS}`Gp!87>WP=8{mjbuy^yGIO# zkJPm~o)Lp4eNEo)Zfxk%(K0Gzo)}`6Y#2niNkFpvyS$wV2^dwhjBz-Vz;|bkAwqvB z3EcMl9;5ew1VDwNpM)@ZP=T{sae^liY); zBRDvJDB(gC8y>`Od+_Ygz=PRb-Vs?(JV>z^MlA+nPr8E!3m%bp;E?(uUK5Rn?l5#j zE*B3bs@e%R#_+(yX+W4dtV#-;xL*Q+C8Y4mf@%1dCB|;caHYv}*t1tJj<)@XzfQsPb8rg*Fvk0gJgUm1!JX-0rV+ zy$&@yh>h^4(hyIR5Hc?rx&(z5peIa%BcSsPjku;2`o>nHur^NRM{#v%AeQk~i1Y^( znG<)`Xz&yAy)7=X7CV3@hm(w*Dn<}L_vmfG`{T$jEsf@b#tc&OC_XnVxQJqtt@XAM zh(L&so%g0K9?rPn=eUkj02lfAb=_1NsCM9tRCC3;pF4UCXGu*M;U?`TLB9nikWg;A z-(t!FuOjCa84m>0zKHkow0i);9nvP>s*lNwr4=*nf zC}GO7bS&i}1=wu!esQyg9IOJ~D$*O1gLXv5X>U$)2-)HJHex$DkQ0gfC+{GKsq(wO zPPveS{Jh?=uQ$k{Sn>l$Kp#1TihC_GYfu2qj7Q}COzb>yPcnHnN@)LGP;(=K65a|> z5#CW#0Zo(kYH27H$h}pL?jy`nVX3E30rzQY7``B8COt$Av;IHirx6X@m~oi%zDWam zm)1I8bkjf~H?iJnDO#8eqR#H{p#@@Mas8=gS}-jQnl%!ngZHwhXNf}SptZ{Cz5N$D zc;1>8U}v)pqi+;VlYZU?ZK&GUPX+_dV#oq-`qRTJrhCGKm~wh3YN55}TcHQt{Y68d z&Hy!O!hJunJ{3Pb6S%!(03(l5Vnar(&Et|PN_mJ8>YR5E&NgFho^UCz5gR5bWuA1F z!v1wT5-9lLAv3HI9CB? zvrg3fC~Cr0^$*bqe*xC8%oBf-z3xg4ddslxb@%xu z!Ko}c?f1`p&VQb|exP!b&}hJn5@Uw`-mUJlVJ(q7Px6t3kTCe-^ys~JC`it{S~=w- z%84{^Za&n7|KjYH)Bw0NUBadf)TP=zx5R*in^v3cf(RreKe6S%CkO=z`jKvo z{P0%7`Rti4ZV;|@B+7B-fJd1%+D9v}4Ybk40TK7Uv4HY9+?0MEGbrCwu_F7;1iHQz zt-V=HQ2#`>c>Wm^i21%*QrgZ8J(ap4U-FqDlj*Uw>O~f)nLYk0iNFfy=_F?|OW44_ z(|$ISY&$$U8M-DA#sP;P7~{I=xqzmq{yJwVb^{K;yyq}4yd60dn@zZa&EY#OjDZzD znAchRE71smK#u#Z8<^0T>xn!oEXKq8>(slhl?cMGl$|lFR6^h#QXEFKTL|vxD+LS( z3c-oimG`YNn4LE!MKDxY<}i68)|=hPGEoQ!a)NwT(L&f7ەkla$=#=sq z!;igBhyNsSIAb5ASWmvk4it$j&I-^Ff@E|2E@4R_sEAMq$-zcSPj`u!geeL^h=4>K z9{c`_q|EF}Swir@PWyl*i7?C^AFw-QC=723CRU^N2t!Tua*rs3FdS{+nzaZKg7VB8 zw!W7H!SS_Ji#7+w-|)J|e+l?tgkeUbc!m!{xV8Lb%&pkK`+YxJ#I%}iy=8pXK2n2Q z)-K&*jQTs8vS`^pwSf4q7ZlujI*!b#4xg&y?pt>cO8mxrg7*&cKS^Czuvx5{!TQ)D zTJ{R=t;iqv=N`%lWddw=n`E~FGNZy$xVBDn{k4{y>iIu&1pWV^r)chMwT$4I|CqR4 z#eVGdU7O&E0C%P;9%}lnZnt4A+4hcneUS4X-@o^N*h-E5?iCxm{}fv0c#ccY_;12i z0xWxk@TniRp1WbiR04M&Ec`soOj78Q6-3TyT$Dk|xn6X(R&ZqSPsL$>5HP^|nVIid zj#?|&E)W1SK_0iaDkc84u8+-|{gG!OXnoUPBTp_u&x>TPH_}ufM%x{AM;mG|DdgTq zMvva2JFX&M(#t*|3)?+|Jk_62^NZ=uJ{e3JkgrL^R%;lUO$2|JD)^47&QpXrOHLr= zv#VBd9n&c5SZc#F>3QUr^kSdt@-m`4EC+5P#4vWazt-&whKarYn#V1Ccc*VBa&)b&8h-_8k^J zPeRWMr^G&1?(JuV`q8i`s<7d6gI0KvnZjKqGcxhFyDd zPKB|7el;7-V{DhUTds7eXgKzp^61*GwzGhNp)~!FFbh;t_t8^&GQ+5I!z%tM6LfqI zer-z41OioiW`fQ$!W-Sx#AIx9s-^t4G0j^HU=~}INM_9dUem9{79<(ql#t{K!BUt3 zZm@*1$6#ComTXr_+6V@)Z$P&PuhG z37!?Vmy(w;L7;r;b(y_PpyBsCIS8X~BoE&%wNqn&V@#I6)XBC1Ctv*yK5`-WJx6GwVfPF8a7<$-hG! zflT5mBma>YGS(jnHpXX`LN?%$|5}o0%x~sQR2;{@N5adPZVjR{o4N-9w1t@$T((Lr z{k5(WcQL1M!k0KK;}tPY)$0GTp@te~|2%#+Q3a{aRw%r18;E&AOuC)_+?uD+6gK<) ze`TWnTCsz5);UEM_!9FLoIXa3(7YV=r|~iH-_#E9(RQyi|14!>Wk*;jw>7RT5 z5jj>zZoR^W^^exoA7p2Jj37vL7sa=d$<@!yRop!zd zta|H_8`eK$MvQpK^G6Vg;4JZw!6=gftY(_SroOVmqn8uuq+;9vi?e$*^#maIWmBL< zG)B3TsUb98m4KWeqLV6DF)n-d#HaR$3XtDR#7l-L!G#=K6lJVc2UTS`N883-aPM29 zz^5w!#kSSo7F%^-UwzgxiMl>aRG&UIRb~j6dxh-{9ZU$omC%}an%NxU8$UIAoVEa$ zqkL@R-z-7+WV@xYyA_bElqI;$TEp`D(&+H3w&40~h3=4q9meQflW1wThqKr1&-O$( zK=I|623dpsFibb!$VPn-#>0=j)A4kIDG{OSI9q4XvKa}yZQufzoayF|sk=h3HY_U1 z5e@!+%i%2gwM|g+PamG!*s!f?=krkaP`1l8f&qWaN~qiVi2ttT)UcIo_@g# z2ECkYi)XzcES~Ppk3*-xN+j>06yg0TkVUc81*+KI@cuh%R`om5p6~5Zs~Qds@XG4*Dsq9g9@;hke`UlrOyU z150hw9@$!7xYN|jHzwi(U!*?IJIlc|R@X*aNHQl`2@Vp|l|jABk;;(Z03FRoL}vu(htw#Is*4a{3$ zRWKO{W^TIo&%YxhZ^l}zY=TwKC8_#ejGo^5j^42TA?E*+mr5z|5g`-7oaCRV)c<(- zKotK|vFB%b$9#= zg>cq&Y8B}JHh&YWQ5V(ydb^Tg? z{HKZs2F=*OB*&M3e*~3>qRnyfn>4TdZ#ZjwPm0}oi4E(Ys=5T5?;9r}#0Qp+WIxD9 z_7%=nYc$1ZScsG~sH_a7(xxS=cT^#>?|0wm?yp5rd_Nn){Oiz_7xyH7tzfmpqNKqg zju!NCTdaTYN*lr-SAMARtqVQ6V%JfeFo4{sQZ+Np@5m-l-1yhcUua&+ z?hxIZX%rDzLrwg8o`7~xe|S#+ZWTFV`!ABmNZ^Kw!^F}YCIOL`C}QP74PSpOTfbzd zhpDy3yCy!2KpmBzqEgBXkMWAf?CM#8W-)|VEPx#>pO>7UdA%Ku9dG{OewPCVG+8QF zzHz|ljEeW2A`Xxl%Gl=5!~y;W`|7CKw}ar(hK%uEtTyO5SNNRMj}_b5tC{eNgayo) z)%woFGQp_1XP76(8zj`d>7Qm{fU>2wo_DPD;I?}HM0PYCbhtPN?KPzZU4N=BLl>-8 z@3E3MeTJ>%-S5_Y38=v(gGcxfJv9_2%lZ&uuhQoXBP_J2LH!xuT!0rfhzNJSupvC9 zhP(Q|uh!Dj!1(1e1>1vZpuF*Go8mAH2y_?W#cgPTq`H)axPcambHj}5u(i#TZ^A2@ z3+Q0fuCQE4bQ{>Ryz*y>-v+5c&ntbIup2bYQI&L{hocSKd~6fwLBs))+CS05nwRFX zln?{B7EN0w_%lE;$B%k~Su+F7m>n#l;AVuKa9-txCnIn~IbQZJXM{sMsp+}oOb~j@ zKB!|K6I@~+IyZsIBiSje?mS z>(V+{W@C9AWBcj42uwF$sdDyf+Gl1+RL=iVQ^N>hnj}H$n1W}toX{D8aVi*9%FC&; zAq5VJOSjZxRng~*NR0> zj6i|06z}3n&R(Ec#k^QSkm=7QHc1m!%yZ?7DOdjU;`bGmqk(3dQWbwX^Ose;TXnDg zTGwgHdRVT)hcCI7OzwkS;h!GO{}L2b>$GE<_5l08D>P`KRFm*c4@d=&kp1qj4oN%9N|vlk#6DcT+nmx^SoF%v0|eQw_@P+^;t8&JEQ4cuBV`ekkjZiR;G+L+g!~i}Y!taMau+vLjp) z0(UQY%~;98ouzmqW?y-DAMJG0tw<50M53g7%cNKbW9s}RS|>F)I~~&S`SnhC_jrPa zAioP%+yi8RdWodbnAy}lCnFeuQN$rOF&0C(7OB+5A8rVK!kKhedJW-N-$a!LE4JCWdP#dUSHlR- z>~*?P;AjNt+#`!K-bSGN%#gJ9iV-wD9(o^?Y6NpxocxAq@KaGbZAV;M9m? zbpjW5V_MsXn#_z~R$_W;`mPa3T1pe~cm@!Ln!~n`QIi|B=J51MbfNTx-Ecp(>HC_1 z1(3gO+wtSF1;&q3Tz&i70xX$3B^>-Lp@lzvB!Al;$Ylwg|5?2UhK*11`9@iRv%&${ zU1zM}2oD!YxS|b|c8HqZd1M0)Ro31g_H%31-!LPr3l9y8n4)4{zeV!lz z6c^9XrZMwDOCEi}DQvJ}s&%C1Wys7rA*=H*AeY)Fw#5@6Aa;q8IL4LiS(ldUC%HOYYy{2*+RcM*cYY8RZ%39gC&?f-coq zi<)D7no8ycwx@ODsOLb+xY5liBrEn#Cu#Q_I(2erv5IyH*_elX+h;@sZH|ZTs6NGk z-nBgvsV~T(zWC54QOw|OLR*@2n;k@)3v+m7xZzV%>%lxr zY|(4qy{0=Pr!lC=FnA%DLKMpEl16n|FbQmvXW>~O3;pun^9$1Dp_JzY^#vIv*fsc# z-mXRk0>$`uvDvD_$uL!qukSUWUo0fIP(lmFsb}l&J;sC;=|m=ZQvn!s-}T-M&<0up zqe!@#4ul^~8WA7TfhP)&PuKVnbm82uxt;YmJy8Dj>9|>_9?;EHkmt?oK~0iVnUiNd!29Um6^Nsq5 zKHxkep1b?#LvjIcEb9+F2=0Cs@X16EP7Z1Y)4~)WzmzTNYJuS12ey|@v>^Bng{Vq^7D#LDj^z5H1(d4#($YEze6Qp5 z2y&SSpo@zS77&*1FK6Ck1|TD$SHx&oQ2Y9udtWpHx#JD*-sEe+&iF7bQg2Og9Ku#t z_^Jb*X~pd753J97^{M=IG6j$ll9;_KDh-ot;}h{YBETmRe8v%DH&C$J|InId2Ggoq zeS8)LOsC8~-ljQ^b~2aXjVtj3|Jtln^7enUaW}PDnGf*3&m&Fy_coW1yG9)oxrv?N zXZ@0`Ds`)(#b0ZQOdQ$f8TS80G9Whd5pkrQ`FF0>1(De@>P@w=pJ}uDxuDRk?z3TC z&v;>dgf0n3@UD02*Sl^Do=P7?_I$x&6lq_O&&$<{ zvmO{(E>(;7s?|73e9HQ?TV)EJ2tQt@z&3|G2BOm$e_{$SsU9LtM~EOcOh?3e8q+Lz z_n6g0j{@u}pI+rKqlG1p9{XL<43L7RPxk)6G|P@F*p+?HVh5G)S#c}by8u`#gA(oFN4%@QbraKo-b(cYZ$P zCI??=OzZpZ$^+fbE25*B3Sf0LGU#) zdYSZ8s)2)~fYI_#H7MSBvQL^#9c+*0m{Mx0L%)ls6x$(n_&KT3efg?7lxwn3dOcSM z)_9fb%n5a%uRVQTi6FHT5}!MsQ98L3PSw`=tXA)YVdu+n?wlGx-$m-H=&u21a*IpP z_iBKyT|I-3wkAxC4cR;AX~JsIfqQoRyMVMHYKrsjF1V83ew=|?3sxzJZ>3++f<5kg z%G>A=T)e(_rNtG&HIWBlT5l0NmQXz9sEgqg+dfbej`sn)@beWYyr2zZdfaza>2zSP zoLE$;tqx=t*H_5f=|Dl%iL*o`I-qQ?yymWnO}=(N^u`*;z(}Q&ay)Yv*t>Dy!&BAa zWrpDh5>bY+-GK{tqUAu=mrX{=TO24UJy}ka2*44Kj@Q&N?C>xl?phr$E%;Pjpl zw2agTwTs=bA+z7a9P0_upOG9(FC$JR=y13B6QQT&D9G%*p;~w~>i9fp{OZUn zv}4|P=Q7N|^vbfZOYTj?J;wx)Ixp1R(kaQrjcL2M%*RxpHk&UvLQ1%5$Otr|2@(qGY1JD>cm%5P}aiMd+XX%s!6 zv|;~%|A7fFMPK*6_7f4_?PTA@`wR7$ub3Genm}bKdDlnzS+U0ID-yo(Y$5ECh7v;$~~zsnMGmld*>YX%^~UA!-5o;HdRc$ zbUaq>PrnzG z{ALvi9SnJt0*+O89KM zhH9j{ty>e<(7_WIwmltOLo(fQPhH80pyo$lCGp#`;R;I3~C8=tn6P;!l)aH z@5cHn{dcYoK4o!ulW372tkUvDU4Mo4 zFFsg-kp@&sC$Mjsv;{pNs<=yv9g!i8rfprx%`sh&j&1C`JF zLtju3!1~!SUD-uyuxoY{)G^-%R(Gq@PFOI&MGJ=BwiQNbV-a)LOJRnPX6|pPvaArF z@VkcDRc`8GA z;6iD`(syNEps*IY_@I^-WLqC;?AVU2JkoM+rp(jk18y_>QhpacxMF&b_F^C(tQL(H z;A8pVv*7*ZevEZ%rA|<|*2o9joifAifAT@r(S!BFo3x z&95S)x5o+y!Uv97*1&^;z*9I?+ngi_v(iejuZIMoX-+zXpBHOoZmUru*((GtS}%hq zZwkTw7&ELVF9bbi9~~W+h2V`%&!LQ+!tguz64`X9Fr@4zzIpw%F!b`A^z7ggf!v4t zuBC*FKu3noL&6Nz4wzFn{22CN2aq0bBuTduh2)PD>UE4_py6$J{dTMv*nA3jXJ0J_ zY=uuad#{Q?v}yP)DokU2o5-@Id=4gDvhc3a{fjUNJ=|kReoGMI4?JWY>)?Y|ZiYgSONn2J8Za}qrWPj*#&H-G{X zHKN=kU$4(>8U3Xc?At&02X*M%Mcrx?`+gnRm|MIW(3Gh2&j+GC#l`8eH-!}6fy6{= zbB8TH@Yh-rm`@Rd`PPzL(w}oT&0PKzHTqD$^FOV)HkFTbXQnn`4Zbngm11Qu>HGh) zM&qw_y};312TW#hC7K^8k}wtTe`B4Su!aW)>7tT)rT=bSp#yJG=Dewf@;%9J!Jk#U z^%5J_bxdYG2du;q5)>Hsd`1laegor z>Lz>6fiW3QnKo*-NI-ax@=x;Ja&RmL@4KR>1S3AI$`hSxpmT!mv=gQsf5SvlY3Q&v z1c&NJvR>800IDjJ^f5!wh^xrD6=VvP_2IYkw_6Y(zwLME0g^pfqmVhVaNG)HQx1Dq zeX@a1M=!e7743ymZ^a+pVfJu;wqoejE(hR$Q>c~E=Lq;y3#s>955T)huUeVLgK#n7 zq=t5c6Y!fU$CqFt9(W$+rFV~=p?I31aksb&Fx9?bB6{Egbj2aO|BtNej>oe7;`R#J zdy~EQb{!G(&5{ZmTC6%n|DKro%vMD2bWy_}YyN%v_{N6v_=hJz6 z&voC=bB%MI^F7~hJChqM5bYs9eG%sdd({SrgQnboMa{qRiHkebyPSI?-|P++!fDUs z%nm|z3hRO5&j*1$e*g1VCLX{*M|p3g!2|y690}KR^@MoauF)Gao*+46<lC6g^ghZyR3@+XWo*8v{M%Wwy zS=(EkL|nem9W_Sf(d!EulvNJ(ct1$BPUU%ReiV*({Ro_?JqqbZ^jqG$_(O3PpL55L!Q+XZus3waz&7icxwuXMgwkf~ zw9g#{#jwx1#)r@O0*Q>Ej^9CVtOCKmExYd^EQxe_2JCTxk^()227>)C^7?Vb(>+!o zEahs*IJpaEdCbxn1T{h1fTbwmj0C*yd;i1m4nzJ_OvJ`oi*Qn7SWWwHYGhQ8CcM`j zz3_L!GLN;jELU${YnEEp10q#hC+xa4pZ;S4CWgE1uo~a0+keL+?9$XKQ~UNGv^@C5 zw@|&1%}b5pEBP*ut({jV)~&1SU=>q`RV);U@#dUlj#$Z}tmV^MakE|1C*fPon#M5+ zQVLoLOiHNYCu0`~dSG*cM1g3_gKk*YuUD8DIHn>SYqvsq%`yWqeF>h(GRs2gJ{AEg zPxFxJaCSg>EFQI#3GK+({sz6MrShh$Dn{@8$T$8|P=>G&5a#k1iptq`{GX{c#h@JFVYT~O=B(NMrxSWkK zCUnw-&0-W!$e~BV&bww2S#3Yz`O}KD%I1SWXH;!V6cq zFlZcU46p8mB?|cS{5P}33?^+7|N+`4S zx_4%Q0`#-wB1pI>py}uvBOP-bIb;|a7NTyrub8U_05rzrdA3N_Zj~bVoV8lA5GpfxJ?N`47*hMhbbY0 zvsLkgI2ANIa<$M~Qh~t5HgS7UZ2iTx$Z~>PhxgG?>D3Z;xF&ctH@9R z!+a4dF3^k$9%|GuMI6V-CZ4?6)&7VIgFPt3F8;)d|1Ws)zZ|!L_r1^N4?k>!`vRP& zD5a^v&rqt5?I|^URx3CeqCf+vw?!%NISu?;?4mlmixvp4*53LJ03MkYriM~L zRii;?1|KP0m)S?6aE1_qFMjP3y|Ii~KHf4WxIKej3C@H+?H@(Mi&L?#Jw2!*xPNh0 zsu@Lu4O9g@E<+l^rD;Uoc*`54x~tc!L+eZm!EfEy&cb&4L!|e<%C4TV-L>BnAI%;7 zyY3~_iWpNoHpQ7i_dXZw{i(2ZIjvh)^X_}_n6xI@>L3VT7DsM|5jI#}<6A$<|AA{i z62Bk6C(!_^Wkkzi>JABpW zkmRqjS%ZBSP*{cRS=0D%G(kzNvHa}{Vm_bxn8!UDjW)y#tmMQZJ;hwBk;QnF(J#lG z%9Mo0y!0kxvr>)`dR3Nb}$;cl8Y=x z${o#8pW~2B;c+V^*;hzbLa@%urwHXsdS`W7y+_Qx!o2BqB^WGBw)lv48CsxC$l`AM zh_GHPP4W8A=#lTwdc%k>XnB8@O{aD(qE0vjAD+~qch!65?@!gE&kNI%1qT|DkD1{c z()cD+;T=gw{h}G|4q86(E%qxa{ayZw;UVrD3aC1nS>o1$NQxErzy02VW((fCMV@Fy zHg<_-yc4ZR| zgB^(2f1g6kN(a)VCtoZw>_l%0Dx++3JCRE7HTCd{PQ<3J9YHhIiSAV&GZw~8b|R$F zt&}x{@si%T?{S^13mM)U(q>xeLK;@$?6ki;vIptc;t;OzPY+UgYB8Cf(2Eq*MkvLI`Vhg{D9Y$7eTaijpP+@O9}&J0 z?RPrbj|RD!+oaF+qjtJ5vg77hqU?R-U=+R&-4Z^3vAwVtJ)wQLgJ*XSs^oE9nmW~m zc8#UF3Z%B9iP93W%%ZO-vvY){wX_})rEztVDORDxSN%fwsKq{@Yd1a?wel1oja?ef zFYU6?)H~U$u#}3-((*{&zYP6*T^tspSzibLTB7otdbFSFpEDs;>Wr$xSsQ7sLy|ly8pOQ>IyZb4c+Pf4q@6i{I-^Xody2p%i)B*T4<`;SNFRh37?5tvi25$@lokmt(gn%NfM4 zRnso;Omw<&HfBqZUcdTS*!5y@Ivj3x3AP{b3W|XE#2=)A_x{=XJ#1K42U+Vu zeMpSUEaHCu(7+p+g^=fw`uZarj+r`p?ikvUoqM(`qe}}kJVgE;?e+cW*Y>U$i5Hq)e;rbaYW8>XKV&FFeN7dTk588)T`zeOrTULZ zj{n2(+sO(<2wsMbytEcQFJOPAD_)28T@&9tR8tjH8Yb#^~m)`#~m@ZdSuHdIdj3U9Qe~YuYgmoM@-&FtJ2l#(Hq@4y6p<}XyyT(5Dh2xKOBunP!X(0Qilws zJ25E9BK?^eF`{~;VMwJW`>PIlUCOT|n5aXe`$R*RaRr<#m()S|P^rgB=XwWwVn@&rzNx)yl}B=eaw*C9GZ<)2p8 zb?88^8vC=HIyC%_Rn>2-4pnrjS$P}QBW`7{#j|_rk(QYKXF`#B6sth%FOpw}Jo-oLRr~x>x{FcMfMh{4Wf3B! zKO{L>KKL9d9OTR3SII(DYVD~d$?0gUGv-KX`%T1ryQ=lE#U*lT@*voYea zBHy~fr7yBf)y#$dT@OtWnbiUiJK)7sKHRn8>Oy*&oWobNdQd{g*Qo`nJ`|(t z#C7yaKjP%k^azR`KsF~XMUCLjj%h@tmlgo8Eq=g=p!mE?oSQ(bK?{PO@B zc)zqGE@J?P9El3}m|qQ`cQc;*9ViFU1*11RN{a@ObkCA$-{cUoyC)DIYBr24md}+J zFpZ*<%4aTyiH@V=i76(va+8SC`SLGqhbdI7`=ch})=!jBxy&ViL8(Rrx^C2t%%B?i z`#U0-<`D68oPm_=JgQ)tBh}GZKmv|cN*OzGi|Av2bZU{|5*j)8E`v~c8U4Dc`%Xdn z7rK|je^{3DHzMI>7wo24K}@0LWhZJ^(9!T)9gy$`rG&?rIr0-g$E;a>i$4M6S;qWK zDj)!b#a;s?LPB^vlWu=qkq{b2>`bej2w@u69e+6$`?PF5aQhn}gohB6hT!;#;Pm+W z$b};qc|uglL%JFw7^M)7e$Gt{nOCXy)A(abymPFaJT1f^@RdnGOBW-JeJyc>GnEAP z505o8l8^$MS7%03ASp;ZKaf<_Lkiv@WBGTFkU=aum|@03j;Z8X(HDLp2OS&VnRhYR z1;>r*iT+ti(8(nZ!!@3w0$cq?n~-Z5CTrq-u0FNd7d{9YqfxzaKlMF3kYRN4{x}bTGiAWV=$+T}F7W zbo+NJjtSULaCm=kX9hL#ng?S}EMOn`CA-3v6&iSc<9Nf5vVoWm!;Xn#?6B)~*weU+ z9B@o1<pkmmVkD??Ri1(IZSG|cfBQlZ z0!r<3Ka>iBHo7D1S0@OppC;G?h6F)`@XFvjGAwCIpHs3=TnI?!s5xuxgdqO*&$O;k zA;3RM{$cw@2qeTPw+H?a0)5=GV3zAxc9o^GsKGCI2fX}4bB_J_4xqWy5_@@K2XLx8 zE9;2~1IZ*(8gLN?iWFDHtM@UU-&7$<_hw9}a=iLu*ftSRm~>P9peq71B>lekF-b8> z8@0q1ESZWve(QF}gp|5=9dd|Z6$Q2J!mr<$i(-J)XqKJ`QNV>1pNccUG|rwf(WN|P z6oc3ERJl_a!Kb}}=&g!(VnBOu$jpjE9H@DdV!xik%B`dW<4hgmK+5H|D`$@ch8{Sv zqia9{j(M>r8bwROO_kbDKb)nY$>?&jz%FU{iJz3oT$G0F+#r?_JsBYLT;YDtCj&b@ zlZCZ|q@CgZz=~!AzZ6`yF}3{4ECD4)c+52Rio&}O`$Aq&?f^m_MT?fB{4hpn^^7!% z6EvTDl?7lOYU7nI)~swaz>*qz(4>_J*txj+zRb)b0YUrLz?MG5C(IV)>-2dQezLCQ zVPnJ=SF;xL{ZsM)qFyNjlA`}oM%$0fNU-_W|6}zy-mdcJ_}o_3fDJ3Q+v}}%o=6ho ztF{w*V26kH?aDQgvgX=Uh6g-0j{*35@La+jhS!>@Vwts~_ACA!e7o1K=4Car%rf1+`{>%3{3>OjXnN(B#fY$XEQ zu&#oAP5y#?&*8~?$qcY_N-A=LAOAnuM)Yet$uog%b_qU1WbFLMo-1D1%lc}|lW$np zLTtG+MAfTjiYUIpoH!NxGKzD|*P#5Z|D>3inmiO4m1GlQY=|+Jd(4*T&oXXxxo%zk zZvVpe7xGK?yd~i@iX#~jS-TnhW0Tl5bL#fKZ(TRi?AW=FDbl3G+~3+WyJ5v|iy@+3 zy=PG1QeC%2nqR3yd0DCl8b+GYkd{nCW65{q?T^zD#fnkv_CDM1r46Ea_Z8Yvf-&q$ zdAH-T#}wKpl2nz7$^9I4IXqr$<;|=I!*`ED+d-uj?qIs-!pQiDS9v~y>TQBYvyFV zsxF<>%>Wj6jr|P{GJ^A|m1Ua&M(A(t%E}C2f-ARQX0&3|n=Bqgi_$c;WOU`@XY#`LNEyrH_nc{6K8NVR7gk_UZXM#`B~AxYmmV zWU65@yTIilgd-M$UuM2H`zwUtuc8p~HQmDkqZ>|TGEH$w7YA>JxCf5cB!F4*s4av@!cU1w8Lt2-sI4n+BXpI9 zKk}!F$oI&=na6_9%`rUz!;y~=nJ5d&X%oyi7iu}k{#09_7%m52vUFSSVf1Hx_yFz1 zJMvg(Fw5ByLj^EaIUMagsQ^*e&PCz;6hSSnZoBm)VzEu2lEUi6mID<;2g#%`B75>sP2pqn+FR6x$5WVEBXBJmDyDCyYMU@0T!2YV<|sw^}BAxOORlCKHS)i`18cIHP1a%(dxI!&aZ~%z1Owt z_p~@BOwy3CbNfGVC%9T`yv5Ci6`y^K{1!&DJcfI+_APGOPJi`NYc*8=$GmFtn}JD7 zY*Q@lZyFaKLesF-@49t0U%1*euy%`Z##0?8W5M=5FJqpyBKyC~I}jObT2a|XB`9t$-T$J{kOVYx31sR>zF9d663ffLqkI^#3E_>l##}|1Wf*|=&5j33W|>kxL_%A z58X;pmkR25fU5EbT7T@%L`$RR4PIJjp=}z~91`zxP{io_WviMzv`oO^|7QF-5?J`o zk}QQsCQaVC+t6z?v_gaDd0vR}(LImX&To<1Tspqq_dQzBbZt1JT7pAgi(WUk4wWKC zQKnsX;^pZ5AK8Wpf=`HTN8$wSpGss_L^|PtHEtge=`#wNt3zDXvD?m|MkLmWU${Qr zjHtXJdbqL$jqu6L24=M*@48Oy8Gw6l%Qx<93S;RK6=R$|+ECLSl-xRouvq*ESaA4)tEPC>UH}7!cETT?J zBDgg@i+qY-jaDtqA`Y4#`0Ye<=wqU#W;@FqdPA?0%gH~72wK~y=y%Pb)L_niHYeuL ztw2ihg4YX zf=yhFbLhnw^Qmh`XVBY-+QSR?rqB!1Wjb;GF{GDqx$rLSASyeh_B&C(8&T+;XP0pQ zhL%n`(GX|XAjRJY9~1VMAU^|g>IxK*hZtR4z69S*!uuCfL}P6Me~~HQNmCyLFiozm z!w(bAvKig~X)!lC3`VkrZ0g`W+WdOdQEe-Y%!U<pr%- zlS)9;)dQ9{TBLx_DtDINNfs(9(+HL8071z7go&7Ja85h6;XD}88`p*rIQ z*I{QB7_tuD_4=+VOel>97}u&{byMwxs7ahUWY6f3&;HQ>ukPIv$+Ma;{kcz8xKRu6 ze7#u@leK|=nWa$0Lv%9sqW`t2PHk0a}!Ne32wzbOcTjk_rHC3hF33l`-+zqE98g8UQYN zA!dOX7&kqWE2eTcz*sfJTs#1vVLHRacM%{Z{js)JF+hTf!OK==U3h!pq~-nNx=>nj zb2o0^8(k2&7IclR3WFD3dTLRRNyhhQvlY&B=weL?Q-(9D0OC*RkY<+-^v&zqwp`ML zA%ZvMR}xfVNA3gdbW25u{RHWaM6%eQJI2W;MIhbOZlt4u5BR=xkJdI~{m~)EC}U1h z0+~@uH~!87$~wR(^+&%8>w6)@hFpZc~2eYE7z|Mj%52NdVA6^R{w!SRzeIc=O3KR6sur6~Nkr%+h$pa(w zKcg^uUWsOZF-HkVj8_u6Coc=j&I9;?YzHdei<{JckFnTD+50fxZYEgkV$8(%D{gmO* zvQAU&m?B8SEQkJSjB%_cHeB{D#+-5HN0jQ7HdzecL31NFQ3i%y>3#@kmIA2}$AdRK zC4nrM>%58cFREF^{d6WWQ+>rPt5(1AxwVQIgbmc zj|ppSd#iS966-o>ONuoj@sx$#2`+P#(XzmAb%jD2lOy(tqlyt8l7+H#U1?hpIf(Y} z4x`dgVgvQKj0 zl3{&2vsV`SI+a6Klrc8n8g-L^y+0%%kL&HdJ#?brBe#2R^?)E;_UjWe@Z^PVuFCYW zFszlIw)tFs3?00EvnMz;1uNrM>`LXfUPO%iPamHm!umFvwYVpZ>NfMh+)*E0i>WsC zk7Ba+kNi?!{#`QZKMSU7ST+Tgehr7|Z;)=S^ln)3@4d*Au!yjdA!|LZ{!i=%E$x2g z@tyHMPF!vCJz3kvH${YgjtANo>=w2x#SQD)QpAcm%EWkwa&iYUxA-^z+U5V&o&LQ9 znkzaK9QpI}-yGBrJRDVC&Y^!N6&)(bk$6Un0Ji!)BRUT5CY~ zSBwJM=9Sv*_a_s0!HLswO)<3%i`C3>y4m(^Mf^5rLqr8%WR zT?bq>vZ|T< z@ltaLVonY?TVMhEYSOGNuUi7~be~p;vlTR#Dxb@du!f!4^{||(qQ&~?{rnb5=9%sT%~=A6B~J><%WS zdt?WM4?>vG0T)s-51>i8Ww)Qi6Yi(J4>%}t2qNaIn5YkULB=xSgps#5#Hp@4w`DjC zBfp2FN6ZhywFv|kPP0Hh zXUIG^|48BEGdsYsByifeSc4M357R@ec@62&8C~A43)DeUDu?WqfYW6x*=Jtk9n;@wPTeTfu}@;tUDdbp>5`z)iLH}z-MV@u zGTDGqvf?~B+BOS=Pb|KtECG#?Rf^=B9If14IE3Wu{+}V)ZvbA|^!-`*iiFMW3 z8aZiot}Xv>u^(>edR|sEY@B$3f-@rw?()51S61vi4y z#-opGU>e(>4s#EVI!&R7-5I+%Jg3p*yAchR&{?E?e(yP&zIpU1&)ribX9>kD*gC!G z|Bd98)K4~M5`g86q^pNAh`{KJYcq2f2^=;e8<#UB2f9V&{>$`K;Cr{Tb((|*938fW z_VZ!FI&ThtmY$%;!G3F;-S!Da*r_NuY9fTmGaXSZ+!@CPEjQ?^g{(Q@!^c0CZpV1Q zcZ8AT#xft!?*7pwjHU2LI~;~Qxp%-6Z5V|`ya?Q6rT^`tBMz6%zTRw+lY+?3gX5Pl z;k$9CKK|jaa=>%R8PkJO0E4$>2B98GP{(Lzsi1;Wfj>MxxG@r}y{F45xgVq9kWjil zoE4x6&RG@@rvGTcj)yY8M~~{j+wBqTH{~ymn?WG<`xoJt69A9>BBWR*07e8B?6EEKCf?2*>_adhr=pAM5ICfkQZan>di{-N(#54^b6y2c@bR%v2Ubu6Z z%m~q2yQf|Ha3rfuz{_AK+_FFWY$k6f{PFb-SQ0S+#;^6EL9F(|Pp56y6nd@ZnM_%m81dQX%ul9P;r*g!6e8AQyc{M*M;$jOCdsnjNr$ zPHpbtF>Y&!j$k}^Db5<07lme;3$5YU``_(nL#^Q<_o+wp4OXBh(bN0#jwP&^zDi3T zHHTaGNV6ELjp2bvLBR*PT_C67H3mVt0Dfn)RN6Ek;I4*{%F#Cp5a2dUE&4|ccD>Da zrmyD#!ar%Fp&qo@K6N9BLTz>xz@4}ODrI`$wI;gJ{xIR$_v$q8-YtsArmL;}f#rqpjo?c&KF$I&r#TSf}^RQsl=`8Afl)=7Z{v6`^ zk(>B&C=>~^=~#}Khob{L3Ap(CR}rx@>ZrUHgW~6-y&}eLAQ$)7U-v#tKp9E2p4ZBf zQPQ(hVvPcK(eT5``j|)ek!am`!xQ62NdMkqArZk-WY3-I9e5xIz4cQ~G2+a}p(iV4 zz8NbokeNyJGnNalk%mc{K_=&0G*~hv5p$v#os>Ri63O@hrM}@z%kC*he^N)J_Wt;U zUVreZjF;~I;L$+41qLRgtr%O-|MLwYN53=pgcc-F|0Ks=pcQ>rrI-zDYef+> zE|YC>ZRi#C&$A+DFc6uYdxp5`cl7xwd)9R4ceJP9z(NSq+j+5+0&s9%-sc`N_Ckhqc zE#K@$x=FW;b27V;bwm9x#dqE4@P1dJZC|<(4`(Ik!J2MVNvu*6(9(^5oD9D2)z^*o z-BRB+-_VVI4{rDO``nFK)zf}3E_S2AN*<>*w&%*6z$eCJTU5< zj2LyT-P?;pS69m2jg|74hG9J&>-R=P%l^sgjLvpfi>z+?PgYLlg=f&@yG@72;d(hL z7kRd=GkV>McOJm#;r1GZTzatCSX-8Nx3Rul`U{??a4QFI<(1Zhn=%ZPjHmSSQQS7_FVG4&$89s zx)rZQPkdo@;Zfi+(|fHJ`dg6rP1%{;l@jZ5Hh)fcRgk_igb>j5w<%xiKKXg63WG=(bpCu&6tomRFwb4r6P3+O=3`&!{vXF z=M?8VOLk(A_&vvbevAy{Mv3bMZfyfe(^Df7c^E$`Y%@%3mSqB$dpSx(2if4OjN7z^ z4mVi!-4tm(#1GlK&mVn!P6*17qFwx$2qb;ED`uA_0TX-U2Tc#l0N<#_t$}2W7_x}n zeDJISt~<@dr^-(^Pn>HX(p3}JE?-5`1K?757fX&`P`-G7wX{c@I3uxBu)dC zJxmMQ)iH5ZwI>fA6=}kkv(dg80$Ok*)jdv>=o4bnVM+ z+Q8k>4JQ?}Azj6}EzedPenci;P7cxrj;gk6w{L0#y=OnC;xlb{aI0Q&s0jP-8D}eB zf5sX?)7W@+e9;C~dR<%+u1*_l5#=4H8f{1r=oJyfKEElP6FgUqwW^Zg6s};OAAYom zrMGB9smaor-U)3m7-X1MQ~-_8ZZbS^onFv;g-9tNISbCJlfw+PVhqpZY{aK|FR+i5*`Ky4D)=T(SK3(sy+)e z6-e#waHj%$?qeyD3`?lP>gM|ihMv_0myHo@Ukp6BmN#m5@e99m>-uWmh$il`y7ZU% zzsJY8^9%mv;!TUq`$v~#agO(vS-fGzpB1K9J(m8~fII)Ie4=b0ic^2M{ww`>7+KBG zL~gpGdKs_?zZdA);&;P}|2|CKMvTL`TTg92SA<&l!XDIOu;QVxklPOlO3^Pjpe}P- z@(emr&4b0hmowcc-!^IFi(4PsaY=lvRvzn0P(OX4YS%E*a2M=-dtnTRVucz7)K5(! zx7lCY0^NTi9{(s8`#m#Aze#se#BmN;AB%oHdkJeW{Z-d#Uc7`vUJnupZ2OIpF7B81 zKK%y;Sf4-gV}t`_M-`*DaffPF4xvC2NBm>bK z`HO7q6yT~cG$Uh*F$a4_Jz=kpqyovz7U#t3ZJ7SoLM1mhhAObVmNIgf7Sh~4Mn#rv z2Q`-phMqV&kb61w+|!62gw^d2t54EH+{g1TtKt|yf#kwPN_9raKke@ItBnz^`pfKI z3}k{Y4e2G*+nM3Gs><=w9A>a+pdib(V1cj`0+fGn6D-gu^hstRmKBr2!o^eRutDR5 z9A8);8|1XxJ?6U24i^V{S{uzcAnTAyr1~NU@cdz$J(kUh(MWuZIT^qO#xLYpLe#jS z;<$h#j-CgCvxUiKXnBEX-XqPEl@ECQ-v90|=7W%%UUkY>_(7qjW>34F0I*#@G>N-H zBM8Zld4E%!6vV(erwc=41tET1wa(W<5CnGVt7OCqfZs#UvjevL;4z#q*k;KKF>?x~ z`{X!*FjHwsI}xKVx}kf>beJAQPhRWyoT7qA4lR;cGEDNt)^>t&cmYx2TE4~a97DUZ za`VN!+fk4D4x{UarEA&Nuucj+Y}M?_`Mti;{}WUl`lFzF@k7PmRrBI=>~uo&rd+6( z+IdBLLj2a*y>7*)F;J}XC0utLTJyV!!L(oA*DR3L{I6fVs$2ebuA5ek=YIWtz7H$6 z_}#Fs=2}eGvbZ=(_tf|z&mr`nz4~zu>o}SZ&pu~z=qK`c{Lbi(`W)IW?RBl7atYD0 zHI!c@CV+$A_=|N2iJ_8iV(jrDGHB{BwQ$j=g0}G*va##5@X1cPO)Z}u=<8i}LgD*S^P(vN%8`HdWrrmj}AK`NRoV=|r+vdKz^L<3a^7D;Ao-%gHXEpvM0G>$j=W*`{!p|GeK|tPzAU zT=+#aY5*P&QqEV+>q4^)v%$cO7R(f=XnO=;>JZ^E`GV6}HaSQ2{-4?mA;?fq5R@`u z1Ln6<9v0Q)aQIhZ%QuK)0x0_P?sGj{?igVBE=CiA+>3Cd(zYX zsrhz(>nvQiuHkp^bTLHat8kX1_)!j8HTqR77^iCee!>GK9J{=K#!j}+&H_y&hH_0 z<8A6^;w>&WtZPY7tXXG{2+w0ru?h!r7s^~qfo#fr{?uxEDEZ}cLYoT(yTHNu{zG%; zwt@>btZVt7Rno18MEGO(D915>7StqL!^3b@H%iaDaS7)#vGDNV%$IXtw{7M5-LT@P zB*=`{?`35P{TxN7Vo?(Rqp?o$n`nrT{-@yWBt$rKJ$HQX9n`nGK}w$@9jWXNw4~E} zh-^)a2zR(XMy5KCaAjoAkZP4!kRL@Z5^{?w$`Hs$q^t+Cr;J`8?JAZ#^|^So@_3Qq zi+KTRR4q=?b$ElchXk&X-+zzdgkQc2Yr`t+2MA+4Bu)(^`)Cl zOc$bcd}Mv+Y!906{uGtw--p_ zdo!f86DSZEUPvcTqH$}F;G1`+5NYg!fllU6)PMYjCt2P!!gqxG&t}aa`Di(0@^luJ zC}iB;b8ZfG5HWsoRGUW}5=yd>^Ya)OLTZQwE5@3S=Fct8FQ6Fmw?2~Ei>PHZdEE2F zBJvu)b=NU&5$RC6lE&fgE~4Ou>7Aob77+=F>E!6WMf8&KRHG!Og;5~rk#{j~5!s#J z?Ry6Mdyh|<$qdCTq7b9|kuTyG(PW3FW(X`I zx{Mo6GCd3EseEShyQd4Nt}ArUtp5Tc`dOVCti6E5>jElqBD4$0!Es4$YJ48?Jz--i z!9H6=MPi@-o<~n)jf;2NEFfc>bh6;a1*GM^{58;O5p`d?BpUy05e{Znf$zf-lR>YP^ZtMmmYhF5;|rPkx$`y?yzYT?kA(7@nsUZWbsiSs~ABH z)cT4-2mAZcDO8d*?$(8RPv7!CJJF2nbHrymWUEm$yG&|yXfd*Ytm~=6Ip}nL_`99) z3G0>ky+)6&o3E{be=8JyI=d~lN;>iH38%KmI-#DpDGlO&o<}&`a4Se=!@8Qg#A8)N z$JM0~gL5l7@X)!Z7R}sF=OfH$#H9O%R`yb0Btfr!{Ny>?fvWCwKZ}d#K^dCWuO`n9 zptHX7{r>JFh@Y}!wmWVT#VEMm)g!=?WP4>#PfgFG5n@T;T=7887uFODvm#Ks?dmJO0#Ue=l|ii%E(T5Y`=@u;7r|fJN~u|><(QDMJBSq@zQ-Zs9qN0&vReVca{U* zTWzyxn5^a}Lh?PaxIB6IaM~uGIYI#xmS{!nf)%0U^0RjTlS**EOs$^H0xP>OPE`DE zR0em?hy&XRReO#OM_N)XEcrb@Z5&b?YRG zJX5leF^!)}BdcFUd%jvHp${Fxuhr%@O~j>Rg>bw4zb5YrJIE|u*mTIb1&Of)uDEPD zJ~pg)!cfXpvXX(Y|7-fV-p_O`E~xzP1mD&mUS*40`5B{I<4-k;gN(!&X9RMvW?4#J{H`)8AOHKDsTG$_DD+hQ%Z9k}Je5v)x>Z4J!-{vZpiH22 zh?g73Pzl&S82arOtS53d0oU&9EQRT5^!^<}O`a@SkYZCb>+*K&iBewUmPg;P;!`L| z#EEeV@}9z4(XFVzz~M>VjUGhr^(Z^oco@-2MtHa@BX}fZa96Y)KtVfzSM822@9cD=cl?vLw5w;5mjL87!P(#UoJD^s&v7(t$2`;U&_j2^F*D{-rn7 z;lYp0F1@p=5a}kf^!~0gxX>%kjM6K?QSNL5nHmM4Ym-l||AL7D+k4z3x=kbxs@GTq zO5Nn(jGz6dvlFtAEMA*T7by!?KYr02P?rU{vK;dK-!d@L|6Wr)Mh490$4!bfWnkL# zhi&wvH2Ap=&nCW=1_oPwakCU@h&VcL{O-0iM5TN1oO~(`ixGQ6_moM4n)JKUj$hKS zRC(vlQyhkgf72`${p+C&w9ve;dqFD;^1ku=ab~jc2vxH*-@>GF6V=LG2Qj)%H16N6 zCI`2vE-f-8$bqA9K_T6|91sW1I|Uq+hvYqGZikxXfzV#rTVk&QsG6zpcJwPiTvbKC zPK+Ytrp)LUU}d|LA9ovcnJYu@8}lTb*=ZG67WynOT%ro*g-%A*Kh=P0TuL(+tAX$l z+U6}uXh1m2futxp4cMkx^gi1hBg}Ir6c?>ig)^A&1MiE2wimaLuDP?I>y1@yjnIo2vaT7UXVrp;23> z`??ixx}UOy7`N?uUX972aJ10v;%Gb>fgbbiqx5?oh1zOF*!G&nqPx?`wZtnPahy4P z3sa9mGn$^1%6!Sl{*=QKM_3wi_@;kZ@x)!!LLpvbE{Zj+eSc(jQ1AhAAU>orlahg) z5*ym$T^^%Di8VJK+|EMAX#t7tZ*tJyn`9#tlvoyCaMuCu9iLyI>ISuf7?IcLoA%B5*Fsm3Dwv?{O9WO;)diFoqx0j(4b?W*b-jt(lVK;;1vnvom z4}*$W@Mq*%_?6!EY8Cn-P(&Pfvl_`$e6-0+tVMeCfnUGU)FZs(&o_Rl4JfZGWPgP% zt_f>8Zy0b8|B4o+a#d`(ThMKeXSwz}+YrO?Jhge;ccimRd)M6UP89KSJMY8JE;M%Z z=RS*-ZnVc|L5qZ_7YT;(50&@yqNCPFe}=#9L!FuqZ#jbcQA)v(>#H3DXgJ8r;X~d4 z5=WJciBf||cfj+=k@P{t7xnP>0o;xuI;^ z`tKcfos-S*RS0=E7z@f9xK9CwfFbRachkuo9bCI$L0=| zdOY~+H@1)9LqX)G#G!K_L-FDAJzM;4Sn;wgm{6vr5zN@HUAi_M8vXsd_sA|5{q_6& zLUO;9@TO*JXU;|HJxyd={jOW_v{_`Ue$8%Gu4Vt5OzWpImj_qRyWc5AWamy~|MiDB+&c zfj#ER%+RUPTB$aQDF`W>H>sT!24kJwckix9!H%8Plsc3OFqMDu53i64@Y!VD2&LA5 zU3nT%ctIO7kERV(W+HIs+-~~qDM00p8IO*gx}crUJ!KK73+mrza-C%`(8;v0iMa4} zfB@%99l}_G89tDI7I)AAei@?hM=V#YyxUx3emkUWuD6;Z}y9C4)~A zi4?&qWS()EQ68kI->Ee3k%hgX72GBXQjj6(A2POxwM{Y4Y1VuYg&Sj}&6Xp=aIB#qvY?${|DJ$NbRHU@0rw%MSq(^9d{R0ucMt^C{6I0k9owR_c_*(tPKO40Eal z!LpNu+vf@v>CPm4t=_u>B)<_a_{a!D+~qx@-8exJFzz^XD_jVp<66vn>8UIR+^o~f zv9{t+a{PYexhoQ|lRV{t|2Ii^rNeQMY@al&d!!ppIc;547u`G;)?|vQUF9)HH zEO&om1^2UEXL(hpVU=SUC4W2#&to{!sA?qa@l{NeiIn^POT}-Bbnj4G+N(Sa|2y-)QAKl|yt^sVwW}Sz^WnGcmYKg{ zUA=PcS?eMr!iz6bcwhoX`i$gj!E|%2+kV}{FIE>NvXy3s43rKG48JWm)GGZ{}9_PU-U*yR|XQMj$ke=>JtjPAK%*Ia>Kgj z66;3jSiKY);w=LxeKD6aZTq*nY>wA#%Wbks=r}>iE@1JvmG8~|16y6LTi4S$7)CCL z7%!Ja$%CbH{zugnbFV`B)fkF3KH!?VZsUe;Zc(xDUzdj{^x8)-qQA zuspzeW1#fW`(uB9oExm_KVC9!x(}Ni|Cn~}9rG5y8`jlVx7zq|41y_Fe;iBvjkWh~ zHuYxVF%Me(aT;DqD_4blNT@td$u4C0=aI}$n^Rl-ZCJ7Q)-zOCb!bA2uQ*O-fW3J8 z(Pi4zEZf&*9=3@>+SQ3lDwa5V@p4dt)M>1%#USy1dvs0W7KaMQm)G|EZ|23F3IrXigLEFGtv)pNIYGk!`-Py~r!e^GRt(dI6>mOCnt=U6 z(ho=2R-d%#3;x-@o9|cOmdZxEG1QOZ?=KX~dBLX2XwyN5hQ1iizuL-bxnaeAAw9#@ z=jelgm=~m=Y`WPR<;VC7h#&cW@5hj4SzxLEg zqG&cKA3df-@`nSYRRyXT&vFB2?=JDZ^9~%f*YNs&H4!!D{iD8n_*Gt8+@x09A*>mHJpSg06Z>^c)9*>LkZQ zCo#;wqNyX}9FATW6ome~ioc}`19dG*B89r}j^%pNoQp1yKYq2#Jp{mD<#!Yf#zo8| zvajx|isDQYom;WSEfONuXjdT|&s5ZVtFD9@9E_}}BPdDWP^>Oo+n~e1E z`<-x8-gQFYEC}C6DK?DOF2@`H!gnlVgkd&D0&3K2*uzatxCEUn3NI3F{GDQ{b+=gB zOg3S7P`zeghVoVksSPWZ9k*VcVb6*12VPQkV4PEOfx%mE`_B%i=a7lj{~e)CO?b+K zzWU$UqEp4jGVMgMWw~uwu@N@rFn;79$umsB%dz|&>*_LTG)a#6&{49wOo&vQMN%$Q z{vDyj#hENmXwy5WUON+XTK~?L5jL!Ad-(Z}lu<-D7AJ9WHm$p;f$9>;fN?qsb6*l5 zC4Y!q`2|UHUp+#yo?pw-sh^_7%A4{Wj@jt7@pEVQL%E2oVLNVsJs*9%R^yyZ`~q1; zSK9>reu*-l)2StAyh8VmjLLcW7a+E6F&wAaijWMqtHgNb8)RHGO{&TM4!N2<;!g~D zkKvFUN;8O!iqXKGuG=k|CCCMsv`$KvqP^!7CuH0{Al@c=+w|BneqW zt;PKb%`_eND@mwCEVUPdOaiKra*6pNPe;$2CEi{nV|6Y$kFp;z`Q(<_NDQJMdhvR4 z|6Als2V-5`Fw@nq9ik1)GL`_3nz*;}b~bQ{jKk}&6lnpn?E}|S6%ak zYx-9Q*FT%(=CA$k&~$LQ(#j{m;puax&z8miaZpN{HUGcLe}`2K@_*jQFgku)>UE~B z{vea%mKiS$(mGijix=qgE=XZ>oYTMXaH24WV`H7Nq+kH2qu1wG$|8JRj_>Uzh8+VQ zA-85i|I=l2fd$UJ=O6YtaX4;zct^WdmCbQ}Y|N{ad}hbOo!iyQBN!Zm_qyyqW%kSA zpmpYoq(`3}T7p)6I~e`OVO=ttbo7Zw4gyuag1{;9YYx3(VVe{qjyo_I_}Om|-3V$s zeP;>uI>HG#X&X^EI~>vceM$7_L|CeuW$N?gpajw=mcD^){G4a#XGlRsw@=s$9vwTN zUk(@8Me>u!t>BxC2=M}PaaQIbdH^gXbLsrho#Tu)A%REM+FA%Gc_D&e(5(t8IVsOmk{9{9UJ8rw~zEO zwcFjLB|G!!HYG=U`zCIR-4{%oPZQ_w8}1vsuV`N|moZ(7Ys6H4L9gMyoPGBD^68_! zV>d-_@{K3`wc!_elM?bK>U`R0?(|H|sh6ZM(PjLV(wm)o{#O3#-%@TT`!}Vhl486hC5fuYXO(N+_rv6kR2xm9dF9Hoq!02xJ*{_D zElE<_@v*;WZNajnlb^+HXAZypzobvUN>YM(Mpj7rizbk6t(i(qt4-9BsN|Yy&5KmR zreM1@M!KZmJNNfqE=gjxw!Z5wmQr*_msfmbp6Y7nks~Wu|DwNHRqM*{Q)(^IYS+Cr zq?V#4@oHxN$NM3&fMuGm`S)HOnYw$RbD%`Xgql}yq}sOZb=j$v#DCBltGvVXPNq=i zh+E!?e!iv>ijM(HJZmC9*)r4eF7HH-*Q+NAu08+XiIM!mQWsI*mDGRGy3BlP*rY*q6AKRr@0D{n)k4opwgAwD#RUcW}f>=#j>GtmOWK{??V{omQ@^AyOzV z?XUI__2=WRJowMFM?Hsz2_77UeKu7g^eQLM1Wsr!L9LO8qnEKVJ zpO)4q50;NkfSnsc-9c?aHHn+{MR~1DmTdVpod!gZ?$Q(nKSoaDTMA0u*}(XBbNlK2l=$1+A$LQh(> z#)d1)i+**R@3iRD!6hOEHD2ajrYy^gb_$XTUia49ZrRIa z)Jzh}s7U4~qp)XL)IY9AWK-^0@f~ZHMZJw4%a{}V8z`2X_Ml}sb;yUfm7@Njbr~wu zzfD`?;lX9@?n8IgwsX<{zMmJ}eUI%gT;@KxZbN#p#L4QXU9~Y`7%DNLO zl+SX0C=gG=-OSCf8Rd*QHwwCZuQSgfOR}ZjmA6H*V`v*Kp`|;j{T7&_M6B<=Ax%l{7)4Ge0j{N=BWv>TkT>e>~+ES z+PTHXj{s5cq~_0##^`>i6wwj49?p&YFC(hWvFk$-zZ;Vk-iaMA$RgOHo^#)?$8&bD zx2}HNRqFu1S1+P0W1K*rab}TO)`f^;buaX~-@78>*w*ctyWF9?nqu2_f(OpK?B3q; z)e}w>HKqqoc;j@=EW`10KHzhHdzycTFJcAQ9~K7tA!jy|THneaw{~wzu+<3w*Pa_P z7iBlW@=4jH9m0Vi+&gbpClrMH2Q;*|Ne1Ip;IJ8kS_oKDLN+Rz5JORa@L+hcL#9)+r5WE>P-Ov41};vpUME+VWl9&0Uspx-b7a*p3DoO%+F zC{%i*gSa^n_F6sxnLJ4l^-=z=a3u+H+M{do?2^%2sr1=xJQ+&-GnsdGreKRHx!5|X z&G2Goe35&1Gx})mudQ`S#r4VnmSbb7aLjfrRnJI6&Hdf`cCn}wKOvEi+w@yn0qEH9b@>_(=s8Wc$@j~!Ay9N zb&m2>vnsY7owF=M1<(L9w2rA@g5=a{o#l)jy4e_AUGbU~TIA7Zm` z-u>u>>HAr@@Gd}wlYc9G+>3hYoVKFduX2Z#-Bvi6REmi*ZAH<2^}dSKEZ8(}W;ytWPuvxI9lrOK0j>R!ls|*~4=laztapF>@~F*&rnMl)lz) z_5oRqVB6Q{cbw2LyDjjtygAsS56pa;(ZZxkkW;jSP@UZu3SNqmqZJf1q}Z=){%P6l z{>k1XI<-{X>f4u{&%|sFUF-Op|2v<4tp8CP>$cJo|DdG!VqF)rn%0T3#%I-T^17ZI zv=^6rsQUQqvZ%MzWhY`eCVpG^O&%t1(Y9q#x@W(hgk1MOqBv`q`wh-b#`jFpEbt&y$z{12G(IHQ+`r}Wj#g* z?!KT)#hKkmj)vY0RWV@{60=;z4*6!+iW#cr_)0Uh^lNPgOc{m0amEw=~4A*>2P5r7O6kG}t2Ww@* z;K<~bG*uafMvCw=QrzKqdhDkm(IYh+;Wpc9XIjGXu;_yQajpoQ9emc2Zxw;_{=RQt z#YI5n+)d|~+atiAPi{WEGXf?1?EC_@MBwIoB5SCB1SG@S`(_9cIC$19W`ZRG`!)yN z?0y%HcWL?k>yCya!ZP=ZY(O~1#GJT&)WT7yS)c315e}V-Y2rpAQ#hthf2cnu5suCR z>cK_haPXdrOk!CZ4(kma0)c&DP)wm#x77**tHP6>j}oB}a;$yIw?7z3J>$FOvjWk` zo>;H@)gSsw+G(+?ebMj#^?~Q)MjRBfyzr^d35^yH+~N#up>-l}we+d=cs7)A$68(= z+Pa)H$qWt3_(8Kr^}|77+@xi#@qJ18*UVd~DX9Z={}#=`Vx%mz@$|A4{$uc*i`F{T zR=bz2(0`Eq!Zyo)tnkA-tJM|7D>LW6=%3KL9MXA+idVZF<=3<0V8uf|nKQiDEPrm3 zFSRHFub*jXNRYxWVQcmFEedE$&fmUWS`Eh|uN;~x*T&w|?d#+H3}LCZU$p+FF*Gge zsP64FL-CKjtnv!hsB+sXooi>0(X#2z!AU2yxv4teo^=D$b5^4~Z!heuocndZ(3gml zf$r-=R|g>L^5&L+%0PSzsBarH4Z+R0m%POjVF=dKnouZ;glFZ`u^)S)u}E=XuhY{Q z%q zBth8h_V=o|WGpK7=+Si%Nk7P+VmmIAf}ygj>J$DcpqE{xi2M}D$)x5IkEWobl+#$F zG6hDhTa3RROhK}*%~gY}6l^&4v8Ks81-$H@^c|`xpvw`IEhJCD+Y-ikxhKico?_Y? zxHlQ1RkHQc9?9_9x$$VKR5ISEv-*CXO2WmDvQm4A_me=ju<3GUNfN$~1&C5aCt)Ia zn1jrM6o$8|eH2)e1hY9yZP(RFpbEaHIyIGu5{lV6sh&h^8VF--eVmA-q6-BVZY832 z2U%ukO(KFHzOcJ_DG|2870&6^iQvjpY`JqS5qD2YH+wut#O9MzB@(X^G0}U6NbDR+ zL|MXh?yGZ&h;#{Uk6}x~pu)AT?{Z0aby42#*SaJ)O;NqP;Y(VG<9)Wuok^Gs75`+` zm;}y%l4JGU$w)tcqJ=yw8C<8{G(Y>A46=8vi-meAn8{%^HzMtfyjtGMJ}{U(QT_EIbg>5YOT=%aDTWTd>(FBC?ly5ASC2Y{-J z?Z^IiUfBEWt-HVhCww}(zIA`E1u|ld51R22uy-BpWwP^%I68EJY2O+4C3Ct&K?<0% zcg2?N=bt=PME1^>^LE9{_CuJUrxr-~r}F&p>Uox&AS>@P|Rn8tOibW{qIb&SC@ zfi$FSSaZ^2SxQH18QTdV)umV|XVjWgZ|hr95J|%;{^6QqRN=$Pm6U(bzje%V1s+`^ zo;JSJLXlIdTIT(y`6SF@wS(SJsz^SQSFCSj?oj3V-79GSqJJ!%#41fi{K_z(F5xE0<3{aGH~dx*iR-d&C(3=Ik#c#2 zl|M8F@-;@*kzTRj=`{4*trrJ{sE|Enl<{aidO<(@YCM>aM71f~CqTF~?_qvd0^&X# zW)Zhag!fHylN%QkVI3};ZoQZYwHdXR)}(fOysZArN5(-z{@e>yc8CDY$Sk_>KAb6zu3%ka_qy1=jS7>&H|#qoAlej&Jj3oH@}x zF!FFSo?q0SP2xxe^V?F%u1%?koBvcJb3GMD6N-=IF{RdU@`xbQq1u4K*=sfz~Ph!$F=~;K!u<_3Y^_Xdb@}xi4FA zjliy1t3t$yZbj4NXd>>}+wluk5uq3?bFyZZ2$c;5*}e)HxMVKVRT-9n`Flc|Q`a(} z^@D@iZB-^%>(*!p7-r&eb#b-M+D!1I+%-b{z3U)o8FB zrI-V5`Z`rG!m$3Lq@vjRzs{stU1uc?7d0tm*0`p!Y_7fTwGD@# z9RFQw^IW$}sr|3*kh*8+=|;MLg7y1?;2POCjVmPmMe9UG$PX@c6cvl(mbw2)1cdl+ zybOy^#1Aky_+|=8DQ4Eg8MI8_tcccz2s)hlb5;w=|1xnc-_^XtSZT{ zBj*)+_c}_jbJdNNM1Rnw9)Ve1{Z8t#r4Cv$kxKnCYs61x95Bp2I%CU)WRv24W2{H(hqH4(ZhAu5d)sTgFk1{VwSJGTV3glsgf$5rOvQ z&7PRp74Ehp#RvD?Z^gZ%^F#CS;Nbdvf0QJ&j4#XwAl}?#K;J14Bdm8@#Phr`UnzPlnj9C5kc z{zDJKaf_Lkd*m}Q93{kCJfSoZpjtKej+HF}Uv?8X53xreKaTEu^Xdqwext-%T2elH z>Q_wkTsU+GCwQX2gu}CeNz?g03z^8ER zHt~->`;GKIE~7l3Dgtru83Mj@M?k#1gnpE$5`os!64GSW5qS7A*C#2Ow2)()$7&8m z;HUK(2kqtvbYJTm?_4A$+PkXn#OOuBpk5;6czPtxU0cVp?Or6J?hA%gFh(KQ{%y}u z?1i|BL=^$=Z-wvZ<7H*-~^|OO?51srIs%o|7dP-#$-YvVDo?;rY8R_$-Ur zJ3#i}M*ge+i|IZYe0eZHYlWD<=x;BxoTr=BeP*VyUOFpMe)4iz()-apHB+nNmm~>J zMA0iyuq{if32YLna_jm>5>@LM)w$IxCH+AY@~Mpfv)|psGYba9$r4^rd_QAx{iFyO z_S6R-SSJbA?p|fK+oTFSNF?Ot9%US_b$5|AR>$i8>fCvMZG0}}BP7!4<9a{i+gq%F ze83xjm2P8verIGWIj|m;WZuVOL^itB4>7O%Ir#%E^PL*!SEpb^oL zh?e~urBio1BmDUBfx#=T$f=C$x%zI|K@}HfmTIt zAack+*2`owSl(?O={vF+HhpCug}!daEm4WK!u6?eZxYn3IF<@?>Jot= zhBO47b#%B{lm_<$7ivWf(&7HgV)`EQ7EDWXotU29f`IgESL!bjaoA6|^Wg#!PTP~$ zgh})eahW>nq(5-~KHTQM8I5kHeeXP$<6;lqfhMfe3Dk(ZL}*NG)qJmIcMcyXCgk4NgZE^ zPlDB=NAgpNWSBZQeCfA*mDYpu*k;zO!#MS^+Z?hEal+_MmC{2dm8}@C% z!AaQ3=V8%yJOQOoxaBP#xX0leO*8rA*=QKslMfVbCk1k?A_Ll*P0$|v>d~g@gSSaT z6BZAhA<^z2I8E|UdUER)oGk$_#j}e%Pm~eUE5)m}i>r>ch=zPA*V*yAF-uC%R5~v% ztM5N2a>8crRVDQk%j$b5Wq$A1=Jx+PPK=|}Ea2S*7qvrLpf?oARdD6Y+?8J#AGd~0R+Y(7iq zsvI08JO7K|%Gs=uX-AG?C-d9w-zmUuwV^QDj*8UUZOwYCP7UE!!DsI)sKHxr;_CE> z8rwYeWy?-c6Y-)(tYMr$ja$PBK6k!QfhN`ZO8z-2uy$NItb2e8Z>0{~NZg>p3%-+v z{BqRDdCtMOpN0l5{Skarvow&Oy?pYb3LOIW&DwlCKo8b;@4E@~46yt4P)L3&Bgi-E zt(N+>8fDq{_x+M*CN($bM$yF_+qB)YCDI?ImC*WCO2N7(VMXj_Ug69?EI1#HkFDsiA!zW0Jd z90$&A*+p5^zyZ0umk9~toG>^S6<3|Z3F^z&sWay}v0XQi;b8z5R0L344Z^!5eBBM-1R?=CB)aKKkr`#b79i~v_o(K`}v-z)z=v&G`NBJK|znL5>|MCu{kWa#%(#7q6~@#;Bt=u5MSsOAq0g znJyLZk=-nGT%(E&TN1uRrKzISH7@K{lPXwXbEt|z4YnK4K9R6d!_i>kfgajDYKRhx zxz^aAhRn&+FK0feL4?#Y7Wt_LC8snO2O4#ho~^AZT2R9*{fUbj@70hKG}?aph8nVW zP^K2=sp0&n82f{bBtm|1O({}?w1jD_a_cu$u-%HE5_w??uhbo{y2P z#P8njefklqm>Inz#^9-nI$z&&1(vEf`(w(uPgfNmn<;A}1XRKMV8;tH8dW?n>~pZ6 zR>A(}+se9qD&Xw1Wf^)wdetp=e}7*E#5=tw(_i0ELEGYWo7syhU@gdYp*p4l-%~%& z+Zq<_<488JW)g0^B|cPRdpgXt0gFNH>x8` z$5qBdUmZE0ZzQ~Vpa$tk!h4SQs`&h}zBxWe1qM5{bbrPxL441rmYOwkNZM&EEt4*V zz?*yq^mj#}ae3p(pc8zsTV+u&^@#0O1^%# zqm{7aPszw?%=*rf|MuYSrr-*4UBZ?Ukbb3!9YuAoNoh4<{`Tvvhc|94ds#;(2U&x2 zZU6V+FCG!hFXg+h@UnkVQvc1ys;=CEvV@X2$g9C`?=eqkD>}r;mqw1Hz^RAB>C~8u z&VJP3z6$1Te)BOJEEw|TuG09o| z*Fh?{wx#Q#DYBvt5e&>MaD>WzOF6j>&eN;(u!q^f`=WHuX%7cv`>fd^Z0!W+u~O#q z8ZPii({*@7=Z36y_SUvL?r=K(olnKt10QxijWB=W3GMQnsk2<(_`*HiFu?Ca(()_? zuY*2#*q!-sn3CuV=?7wSHXD5rdW1$Rqs$kUyM8XTwE3d4g{N6~(ih(qeX{k*{cyr+ zK6rl77u5-i0YN0ZY$}54=L=sPG#R|fcE%T{>iyq)B>KXlFk#0QU0(=FOT4khSQp)HJ zw~_+GeUH2_aI`clINu8yOx>#bPF~RASUbE`#0yRI1ZC2w9Q@;I+Yfeng73>vQ}iuQ zD6Tzni2ST43eELSPagI}XV(X-VbW_MIm3y%3R20+miH~@iYMN7-&7*DbaT@IvVavZjEqUYHTx_d$D`H?rlZJ4-};a9#SGa#*nsWZ#!; zPv!8%0^L%mTz=QqSQoh*jLq&7X%Ln3L4QvQrUP$UY z`+TkbMhJeKStmr6=>oCr$xfD72W%To*I8d{jTbI7IhM2QQQK84ts_gojXHyG8-J=} zZ*qJ{ccct#Z%kAiyU7i`EuTJ~Oj)c`X{OSlTuyo9`OhyK%)33ou~i}um+DZ~wa0xm zd1jXl=9@=5?JC26Ql1ytGTOnPPpmYUKWN=TIQ0M}@#Bq(Q;f8v^Sz(3)zS+KgfUc8 zSU6EZYjDt0n}-4OpA>jIKCr?n$-HlmCO4v6xfEWX5QOxNOYA~*Vz3)=+3qYWjn3+j z#vVs`_+CouHDpnN=9(gtE>gXko<$dBcT)?jvA#8e5xO|K!<9YR*%0bAp?e6INw#8E zQlPK=+!)M!N)FLKO~7ejzkyE94CKiIl{Nkwz`OD6SB6qcJc;!ozUs7wVS)OWpln-| zj%wA)2=*9P2hx44l_eA)s$UbM-?#esxv&0#p@&#L0liuHr?d~iq zBK>}K0JlBqwg0Frn=S}RI;Pmm8FW708?xn{(ZI#dZ2idzpE+-mQxi*_5L5m=zShhMtJ}XHcbIa7Pn&YM zz;Q?XX#ZBL=;H|P1H{LavW~dB$9N%!$`MDbn+7ze9H68xW-s~O0pBdLQrKA>Q6TTm zNs%6Ix7d*p(7hV%nJGaS*`?jVpkZi{B!uN#an zZom`Sy@eK^j4{EX{zSc37ZE}%-!+7*sO5y<(BQo7yO=mA-5=|B&P!vj zsBP-oR0Zs~OZYi8r~-#Z6Z(BY8mMSmw>8*O2SPUYA1UPNgC@USXqE|3ywva}ZHqBH z&QXOnu3itRpAEJWHD>U=Ufa?ZWk)ymQvJMNeo@POhit_J;2@u`pFmAH4ik?mbvW^nr%+dB>x2zOa2+ z$l!9*7n^532~pYjA-Ruw>QlcTDZD&*8kyw}p;Go@Hn9MVub+J3-WY%%tL*l8Ms32H zUAtp=_ye&e?`K2!(?EpUj~?pX76iAXn+DZJ!7#d-#GghL0tU_(A$`w6&>c{$5mOe5 zI|6aDImGxdus<8}eB=@i+hLW&8ixod?SIYe92JRYV+WscUx~t7wyycNMbY4wywpQ< zBXvg4O^w~5h(%@T@m1|Xv3RjgUa~AZ7FnvpA8pM^uZ{Fh9c3}7y|u@wE-f0TPT3A} zzlcEL|lEy$dHtghQ7;rX9mZ7 z+o~{kJ&>hlg=93!p@+{YXQ+|wzv7ng+td4yid%ArjK8VTKVylU`j<972(DOG(=@ln z-VB+4dNVIJoi9!{VP2`Gf6#<{I?~LzH|c+TUCab|(b9;wGV!4X?Q2`nM})fZ3!i1* zKP9ZWNx}W;XDi{!L*K*d2`>q^b>ys^4LS(zC);Q4T67T{wf)v{==Tu94}bW=Y5$%O z=_!0(z^9MkyT^Hahg?75Q`*3uP8; zh?=uA@fTsog(g?S^q=gI*=qFl$r}#zJ0(#}UFF1v`c0QoFLOa%Oo?vcA~zcPPH^WC z>v`~FU8Z*B054`-cFCkO^JB4R^gaK20o2BeIVj``!Y!-rnrHcv^O!_%@5LG3 zqt`_cy1n~&u$w48G~W;${3VKCL*y6Mof5-!Ti2~eHi(0^fhF4Fqd54s4p;6v)VFS;3k)n%IkoZdV$a_Kx#p@+5XG=-L z^cc^kc@Jq&XevJ;Zj;8V?w_ahj!9!!HTSDelQf#fedmoxb>?O-j@bb^860CjZ)~+z z23@ofrOin)*iEF*ue~LMHvRGMJ}k139k}&~EKU{y_fm;PCtk_oB#VcKx4{}T-kTV0 ztzCnA4v8CS736T|ab8Z>2|1iKwfSt!EswWUpCXAl@~ANjs=PiU4`bKCLY_bcm`uz5 z2zso5^Q>KmIph?vFUIJw{Z>U#eVk9NXjMeL_|Ry$kP@0Uv+weBQG)iI5S+S2_`l}OFvyZ+~TmCzn=wCMdWC76$xD{Jy9qg_|)Xy!U)JayaH_mHRz z5fh#{={w3uSTH}N^G+GDHMjWHo+{%2ll7yu+mx}9Nu1sLvl4EUYi$);qlgaiRePPQ zqHXc~2PNLN2iN z`5kg-IY`V%bA2Ohnc5-BAN_%FI<@qgiFyr zT=9Zm2r85OXOBOeB7}awE4eCvhF~vxzJBf1S%MAar2w{qMZ(B9-80%96xgEVqkguF z8t=kyU2J(k569O{pSrA>uvpSl>^sW}N}kxd`zD;IqU5_h6vc~si5b&(^GRt)Q<6YS zwit+VMWqi4rHDv-YZB7Pjg*Ee9bCF1D3B#;h*mD` zcGbg1*tu&n{i|0dV1816msP+V>&yg-lh#|Jka>ToytobSca2>VzH5sr8S@i+kJyn? z|3cn!0ei?lVwSt1Wse^JQ#!*_b{P4Uw(qJZX+~46F)w|0rY$xLkA&%uSc5x6E$in_ zOHiIw+!@ce0rua2>2JAfiffV@23FflApc&o#Ws~RjBb8f5a)m=KYJc*TTOrz`PUvU zCw)wqEmq{rt%arFV^@=(+K88f0I= z;mK!JNGhegbvILmcIowP_Zw7T%VQ8pW~&1C7Jo6G8D%K@c{;RJDx=(M>PiBsJ(e=| zT==jCX~!;vC`3~!qbaL}gZGXS1|CFniEUAW>{aGp4n|4{UodyiCsTru$C)!KcNDSb z`}p=ZnTj|c>Lx@?kygZ({uu_%HU+%n@u1%stN_{*@*DMkk^J2bn%y=h<*|y)j3>>E z)UR8^JJt484lUZ-t`6LmL*_QdtXKO<>4fMG8}c$asH>?|3tyK5k(*}yHqz(HUwyXr z){;lKAld4JF7i;MlqH-MlE*owo!N#ta)=<@b|hoXSp&~l_jMLOWMFoQN}TSa6l7?; z%EIg=ARB9TZHu2MV(YHmr6~}Aud2TYxjz^7Y#b8ae1rw}9)7bLGN1+Ry3F!GzF&mw z)y~q+`IB`qjdYX$RjnYyMblDR4EpHYqH z53~I9eM44I|3&|J57LtzB=2GFM%h8_VG7~h2bRUWS9S{!4lrL5vv|+NJCO0e=V=<}dmuD~n(9vA-W13?{aH`f17_IY`)EILdnW;V@x3@#D|&Fb(HJKF5U+V%PWSGruSzp84)noD{=Y3TFdNrvy$J4p--1 zlfrFVxvJzM8K~0xRVxOpA2M75fnd%M&42IjEF-TB?6gPWW=?-NGvd{(g zIIAG{J(n9luPR!;?KjBnC#s@xxa0Nb%W4q3-EfF~r#gZdsLiK*G_c>^fA)}xCSE$Y zMOiA5#?C)zUTNpo#>+FCm0Ou~AjX= zuBDvt6|5oPDcx4un&$+p+go6Lq=x{(n9lF(x(Vo?K@V*U0hhl^Qo3Cuz&7-L*31C{ zL|-*B&?XS@faD^Ttq91iZ_;@yM}X6YOHCro1ms5^uiN{>5PX-;=slS*gi-3c4a9Ba z1mw+h%Jq>JGJb%1U_qOJVUu7Z_Gkjii@#edkruMvozhd80Vp_RHb)xEMUQH()~jK_ zm}aZmbc_)`k~xNYGZ|x_wA!=C+s0^Z7Cuv#xDGjsM~=;DnSf%Bi!GLXJ^IJZ!&)A! zNAlo>+w!TVn4368OwJND!y!MZz3rFGaN}gmpn$MBTHR?gSR>3)Ftm1f^*(d>#8qbL zA2P>UT~V?h+2(kTD=%1*)HdK| zD956ztr-*>sEblsO)zhKbI+F-MljvIy6X_XAxfWz7^f4p5!gh#Rh{JX&q~$1uxvgq zi}WsK>9E}*xHi4b+l7N2PF2&y>S=O3ee~`8%};}WGOwlKac{?ehQ}rH&MbR~|L0(! z)GuP$q`!695xr8NllDpYCkK04M7L#_L4SpZ_>0ysohe;9lca0>woLv1$yCB{aF*6O zssFm=45od~bpMW{J@reBRw?CNN&N?{NgtP_B3|A;a4zoQC&Fb3--34}$3ms#d(ry& zEP-xk+N%d)6zGszI~6}h2eH;H`lw(QZ2#!b{lJYAujo09-`es6mXD`aiHKrf!Eld~ zx+J<@^A5gvD~sE=RJljEltAT3Z}KHr6>r1~_;=)KV2v(?E?c)YXj5NkHrwkFG2!HX zLTShlHss388I3^B*S5UR`;Eb3=g=oa8b*72k=t!IgDIx(#T}5jUW1pQlDtF#ET%`3qEhGa9P<%@7@z@xS1=)TsmX}7r~ClKSgc9wXRD-b?BKYd;c(u<4v~43?7rr9n3exT$4)e{!+^o&yG}-Ss2n?5Vo7>6>bc70 zvCa-}B2$a@nAyS2_RYs*)^?ESwNX|h38MV^V8p=94m-QI=X(0v;d%p|plY}sb_gh5 zaZj|vCKj2Ey4iNf@4j<@{eT_(XDRp5owY+R$LEHEdOOr>weO(mwu5=RehqQsyd9#t z=RK-P$(hjjea3WMdpw+dH@Vy09+PQ36q;i%|PQHv9SXTSCpXpWE+G_Yr6= zX!XS>4iciYJ+Cbmju1u)zWLZDjuSq1+<3c6V~XJ8Y&+z0X__F>_sQ)0&!2>A!2^je z%@+ygMjviH`bvRgKWlmpf2D!qBh&EsZw#lPZPbUp=SCpuktMNo$S zDdCaFOsW`^ZD1b#u7-5UJ=af+kV0?Z>{!o`7ARaE%^AMdK@f-8u1}|Q(Vonrv3*>Gj4!ro2&9ir}d$yQXZ#s zK@YCquO4pH(#7PUxJ||_9aLvo7Jo3+#v?`hdAT7?aOh`G9I4U3c9jV>&2#G5+B94h z`cMrDXR1G>O{$`7?(5qV@~Tkj?;g8F3{gQ+b4h`4sWMibdoCrFsRYl_QWMQ*3fS;0 zHskP|987*P7wy(tgDKXc?YjzOkSA_YusACP!Om|YnkAArTYM!$!H<;HFl{gI-Y5$Wl0P6JjN33JZKWtzj?is z6Bp>WYg?(ap>1c9X>7o1m~1ksV}3+~RTFLBEbowE!F3=lUi~}aaQ<`Fgr}bgU%y@@ z%l+QFbmDXUPtHurI9FAgE$eHKZ&c3F1HTj7;rlhl8^|q}_0?Q$&nCK#e`Xw;FUPT! z2^6f<*FWfT#<@J^7)|UduD<$u?OTE)d7-pW`3J&|p#Vn9Z~Gl#+OH>jHr*4Q_zn8?x^nl6Bi%<%-vq?>BaLE0~>b^=j(%iwmRr0x{#+1R@@|p+*98SFJFXF6-5Un)lMiV9M zGn-p0e)n11UT_efs{oR|08zU$V(XH2>iB_h^>PI#dL-`9fnNXUz$ z>AJX+b2*FZhb|-vw=OCM>S0j&s`U9LJy2#~dJBaBbx$?iaX@KbG<52)0K)9N+( zG|%)gzrUJx>t}uBY1`a-Ky82sWzX&+76XJQ2CVPrH$eMm5y8)51~5EjHQA+PfIW?a zMDa652Bc}l`8zhQ25=S1aIH-+z`otP*@TXbb%oi2=;5eUOG&GMwx)y`$yaP3u_o$k#%oLTRY8VVN4hfP!PZwpdVWU|S8WN$Z$tBK@2^RTr~f`MKQaH_PVz)$F6+k(PfR zUuZrzC$s8J)k^dJgD#C=dfQl$Mx~7@i6YAO?3#Sv2!uAFn3;tc!p9MtV0oHZ!bg_7 z5%lX92!1CmDEy+x;2zXVJK;?MXlqc0EAd)Yw@Q20O$q$k{nLw`G^#YTcQW*gk1PfrA37K)B@b~b?*TY3VX0__U%t9ML7&V(WX{{*SZ4~kTuQtH*QAOeBU<;7b zrU|kCw7~VPri&5zmN=%wa^GFl3R))nqn{A3TS4M{)(urpYgGOec~eSlgFLnaX|Y#r z@Z$QK*KFFhNZ=f%+Dw|r7Bara)I?>6HmbVkJaKl&d;T@wb;J&}{e4d(1MD#gelxX? z_Gq}Bb3`uE0W>vXY!81qz`yuQ68&CBY~Vd|beFsn_G_rA8NYM_^@pTs;wo*tL2VZ>%voIpT503wPwRi|=0Y#3w=sW$??55NEOG*{qV%tPZn@rakR3GN8D1Xze;# z9vn$|lH4#x_TORHf&a!HG^v)Y;_QB5w#%l|iynlp#wlFB$% zTFD=@PC1O~z|vQwDA4@({U2}73(EGy{_bYer|-xfW0_v&PRJK+d@}IQNJXYuXCVzv%A{c}?cH>{5d~x{hj^{MS;4Txt#fGVTBSUZAuQeM3P(TOv5SM7v2n z=u73w=N~5>Wir zR^ityjmH#eYs|OE;b`e;>4^PGc$R$BuIz^jSgM~r@FGpF{j&aLo_bCbiYiKvFUabE z{aA33`bk|pe_-n<@;|n|J08pTec#??Wo7S`k-c4~y+>AMB%>6{XeSSe$jU4#l*pC_ zDH$Ot$tEOwCwt5H?o)dA>-GEVc|GULbKTGLa-Y|A-p6?y$0?`=&+o?WF*ikUCGffP z)^z})f->&suORy0u!b8K(_GGX5K*PjtYnN9bcyfG(8q z&03}^>4A*3rde^f9-OXnt_^b3hYJneEGnP$p{U2~eF?7th?0&n)HxbJfl}gIs%Qf! zy_1?ckYWI`3Eq{W_Y9!#c?XTe4Fl|lU2B!_FaX;*ZUYM)0~n5z)NZ0RfG6*^UvS{M z^EtV2uwsiUN4%}+H&yuaI4dXtSOv|@4;N@6IQ21d4>0bq3eIZZ=Qk3V)t0fD; z&n5}fK1hJ4vxi2tg#gU7@A`gZjtxGo&7VZr^UIqV3o5o8ga;zAF$SYSBR`Aaqkn*7 z9!E&sj983utNB^Q=ExzTwq)pO8U~17fxGz%_fc)0P|(r zkHq!D;CSVra4n@6_&fN04CR!>!PR|>^%WR(l<`RE)5yKj@HOCp6|0EIZGbSUzGy3GYa7H&M?S%uOe*NHIj(kQ-lLyGj`RaN}%XJnj>&f z2{!WhhbvQ*fa7a-@XZb-Ky*t-yeXBT@Z%gyo0KvbQdX2N9>6I>K(4|;)w9ar|6Q+` zK3*B#zxxtZTA&Pq#gfKDwaVZ?XdguUSs5le=2HrXl|h>V54^{ffmz{>I?Fd@IJ(Fk zX}qiqY6N3jw;}JufMix!;0Qyp^^LeYOo;Yxfm?2 z21iPcv~(R-gMi^J@{y}*pu5g<8<&}-1|wHKe@gnO2ExtR<$d#Nu$#JE$%0lLq-{Qn z&q%4md%lH}(l+W4?p(oHbx|D{UJEyG|yS+P={kr^pDz1sDtu* zI!b#Y41m6u-+rD20fl(E^Ai;WM`y27)_5SOB(^m7O-E35Ukk@hSB@azWP6%I27*&J zq6;P+pm|XTgCI6es$~XYJifq&%!dsy$MULu!2`G64fs=;Me!CGC?MG?++t z?D}RMk}-%rAG+@YdZzcSPWf>o+WW&CcTJ@kInc->F0)otd0kMX$D|z{s+w_srQC@` zwl;#fnLi_w29oDht=(u#SCJ#22h*CX$y?l%kQ=u@A8m52p2 zv?ucQn;Jve{LQL*7Qv8HJnEp#bP(p zB<9iX(U?jp(Rq}xxGR-YavsIce^|YvKacK{Jt2*-nn!!+#A5ng=h5KVpt*!=SlsGT z@NmjJvS(7BDk`2wZa&1KUGL_Rg(?Q&w_=~S2+@=7n@7gXWq7!nj(L=LyW)&)?>u^| zm}e9|Fps(#@FPwv&!c?1X1nX(u`g7x-6A!Oz0S<*?%kM2ulGn153bE4`F#U#4q&fd z0UO0BcnheYHoYc#a~=g7>C;VY&m&cbJGn6{^XO{&V*k0NdDIk;A|-epvyZtK zP*s)yjwn`U0i9V3^5s6hfNY-zT^jRSKq?;2CgvFnNbxW||mtHI)!Z6qN{^>;&>h!dDANvy0q+m?@>biu! zcix=PEn7lQVkfeCx0cYwz`Kjr%$E@jH|`_o7A_;b+>Kmr{1p^-p*)WA$O<~^t!WUI zzk;^N0tF-HRuK7j(IfY@R?$1n?S9JmRn+Nrv5|IQ6)_8YhN-ZxAsVE$=o7t$G>JjmUNM`aWCAU6u``hF@oW#bX0Se?8+Mb$SE&ME28Hl5L>yV}X?u+A`~? zO)^D9<Xr%)1k{s7tPD| z?;PR(^;}Sse0|6B!9QZH4pq-wQ|Z6NT6Gc27viJ;wcz_5EirE-!Z;bSStp_lt9lV< zOt8?KT|=loDfHBA+60PuP%-4fFo#Bk8ta(|SCGMS-EmU1jhIKft_8XhgAU&-O(BI{ zz%|Lr!8k_`xQzSoo`@Bk*|glmmbjrRNx@;kOAwY98hATv#Xv=xyS5lpm$+#?)I3Fk zfu7x$G^^7Uae!2Aihi3_fe3-S2}ho*LzGoWc@EZaTrod!jLH+B#bv!e@hU)U!|1bG zHh_3ycQOeB1h=etJtsg7l-IRrjj`=r;gU5wvzG!Klyhh$aFqpW;?8r1SEN8H{qS?| zPK+3pHkU(rS`>C)y51W2K@gl>9+w`Q#(*)SR;FWf)jUw~n9Pkx1Vhrub=x(!Il)=t zh@jD74j?W)CbwIY9VAc7u)za1Xb!V98FgZVI?cpbSz9*PEh_bXHh~R7u2p;r;bez4 z#UO3odUmM(oUH0u$N}H_{VgOPal+R)UTFnqF34+6@jJQ61)m*zo;=6h;RfeECbkS| z9*BJxE-L+u2i`atUcW5H3pd3neXMTs!jtmQeNTFMA)(?%tvMAR)D=yP;Y(tM2VX97 zKQ-b5a;q;^{s;LW??BxYi7g*|t-Z2uzyo76h-7>=@#TZ}w)5Z1qxrya!DgX4n-3Q4 zs=b@4vG>Q758^N%DtN8?ESia&AJUGq+PGldEzS_`WX;3;peC#4o1chLxISF{>M_a> zXEWbUwjU9IDm$NYUpzsWle}r2_*@VK=4uzN+6Y1XkW5t16CtRm36p9l5&~y|!c}5Z zAuyp36_d9Ugpb=7A_jV}b%*c8m)7t?zL14Vl%qU1CR~*K-TB{I1lhXr&@!Q(+F!O&0c1v%0G9-%RyiPl$UCr7#i zjBhLI{*SpZwRj^nGzRUIxfZ=s3P;%g;$GCLw*RnA;`sLf{uSN9k#{U~4&uT)jJo-j zhY(MA%iWbb<7l1TgRQ^m8{+62j@dOahf)vs?b+k7goy8!RY>Wrp-h%e>l4qm(6~20Sm+TZ;KO z;3{2nlM6o=5V(#cc>^~zh2U%G`|yC6`@YkH4|w59^Cbu0pG5LoWp8E#1R#EYPzD`N z5E2rLw0x#9k;e5wKMr3JAQ1jok-Tdktj4wOr+tMJ19>K;S4P();MCEsdch(|%%$jj zVo5q?x8eBw>SnSO9I7HcHV`Kb1UWrcM}uTwOm&_)_=GH+w?41ltR)BD3?IprF>3|( z{SWHWTI7M0w?)z2UjbrLKL-GdB9Pn4l*m0&gzzsA^Q-0cm1R?gjB>`v^7!zj_y&r3y)PmGQOH{lv)*&-S6KfJ+BI+QJG&k zu{yBphgXw?5vH^jjgCF=Mh%ogLWWY5)xo(TK1!rQ9jXsiB^L7`*v}+YV&#b-B!K7R z$6^G;yQ6N!kYGdvf3ZoNn-#!3aM(;f1h}H;V_G;5P+t4)^MIxXxF3^g?T^-g7c;jf z9(8EI_Y(3zcMeT3P-!_bQ<^Y2r`nz@rUlaE=?D3Jw6K2_ ztH0w5EeOAn*b=&}1#RbuE%bG@;boqHp+&4VM6^ocOuGBDVU${2xIjn;Xs5%I*?e>$ zdaH%xbBPWN2h}zd{uJI=4{gUO=mOz}?dl^>UAQz!$48u^3(rfK9|l+Jfr7o~x+`Ik$wyLPPXGkowVO-YGqJJ@n)n=P7;gSWsoK-O`6cKX436 z@-7A-W3}Hp^t=J6m&yvLNg06klS!i53;J+ONY-}>qc{wl%>X3_9SCenOkvs71RYsX zIerO%oA1vXGFzxXoz;2s6eI^(u9HvZq9sAW;_+yW2zK*LMJ*_QXNPq;=Bh*qS}54( z9xbba2eexcBjogk5#CeHf$M+Bpn})_5FP(x3jn^Qwy`T?>`uaaEwg=^r~iu)^|sqZ z4bQ9YpM?J%{U2ELKkg%LIGK^+)cMay<|BpXj3x*$X|9Enj1!SmC+)$t$tJY;aboHD1Ss9Y&lTWeMdtK!edT$#IPX z=Ep7?b3EgO_0gR8=li(8$;IAwFo+92hKroB!tm!&HxtNJL0} z;JD#@(y{EIV?5Bm#+-KIGY`-TEokNnMYI+v@m`z*y z96yZ!Jo#YQKJrxnO3Kx5oqQt*GedzR<*$X{6SynGa&^s@+f=x?0opBIIZ z0?&qPEMhRPZT47KRUD=_rObBY+$DfiCamy)DduW8dfL`@ToUtV4D2hxVQi|@%Ga@b zrGWDuZQ=gsQb1BZM!H}n4Iai12IbbJp~$0z0X@WYpTde9=Z$3nr(W=7Y*`l0G>f56 zsd6xSI*y6YSRSTN()Hf?E)Tw1JwOqu0F>7jAM1)LLUV23IF9bMB4}t#M(64)fi-tt z&{n4stQCLL68Fcvf`g~C!z5LJ=Ad;V(S`~hAKw+?eY&{qTRr&nJ$9#z9M z65{vjDb&FB`$5t#BC61SgjhoaYblHIT`!A2uL$P^y&qmGlmkbfCAKHt(oiQ>{3_zR z1k{*>%aqwr2*U@1^N}M-JTP*ka+`dO1*}$T7}Bz6AlyB6fKP=G#NMsZ`wPt@_sY$S z-^V+!n(?;`2!nn0+#KCG5C4+^p_$v3oO1+!9IBwYi79RTFB0=x&rp0dmz;mo<5#q# z>;`fAPlp0BT*z_Tq~53fNW_KIdF#>;YOQxNzSBI4z8B6fJi;`wXeTOcQW{o@L{wW~59&0q^L^G&5}Fc3g?8}4bJ8K${cvRZj`ff5o7KE>4x(11_z zq*3&}JwQ{XWm~?)41VU8VeE_SP`Urx(MNR{l)5J<@$Rc{SbL@RG<&*^5V#83ve6EU zf;&ill+Tv{@!TVCXKMGuI>8US(lTktEUy(XFObDNi9A15M9Kr4J)SOZs|dOt8#Sa9 z%7BN9yHFUR0*Y4r&1_o(SYPBz1P`Tzo&k$X7GJoO<1~=RY!bL6LR*8mGM2+giyhkj^(48 z@S}#nt&B?x<^+39+U&JJs#`ydH%SW;#jE_KJGFqFWce~ZlQzhf*i_SC{h~Ohp2^eY z+K^J0DL*2n10lUx5k$p0ka{wb`;{wZ8YbyRjpOCf1BIUD^7NN_(C(hH<*|*O8wtB8 zM<(<@Sa2HYVVG^MaQ=ry>`u(GR!b#)q7BPer3BmpHNlsaGG~P5NQ=;$^iw`#KJYUVA*eOPCJu(|j^+zQKoy%AxU zA=DGUAXLojg-h|)&i(#hn**)+j>&j$`0N}9nHf#Oytn>>P;1w`uH)VO*ZKT+6gv(~ zO3)^O5D~5`AOk;JVitX8K&=js7m@vTXGYH1b=0hMYoNLV4_d>tEjYc1fyP3sGG~wi z#z(xVXCKkPVLm?ZgQxa@rfCXUL?JVv(1Z6bTe3sS9l^u1YTTd|Ov3BG#R~$`j*s3p z3V>93+rC%1!f;7Z%`7isA5_kqnppTn-WxGKZzM6AKhE&PH8UFh2%|dy|oOPts4psoqS^nlAQ$?t|j{m8MLJ3q( zp@Spb|wVZ)HB&XokNhk|1r1G2!gsHu9bI(0Q_6iA;pplk{#XOz0q=$Y6g%7yMmLoJhgbgYFANrdfV^;ORMk zZSQwIIN&W@Z|JNKqfe6e-EPu{j1_0k%c=&j{o}~R=}ZH-Ihn9KiP{h(oP@e+y$vB% zja!_q*$|u@*!H~LZv>gGTpVwsjKEbV{3T9%(g?_7i0nJ{i~;2ao5nschLlHqq2qWa z@Ql2$$>^vF9Fx1TY5vv(DB~6O+lZTjepde3)Fe~5Iy~80NoEFDMCJt_>Y9Oz-)a{RxBTVosentbAtuJf=OB^4s zblt(c<{!DTEY(>+-&5V{8`&0^sxbN%4wF3!UAfHNYH|RmUTd>Bu9!g;S>8(xcM~{@ z$KXec{qpHMVIRi|wE?|Ow<*{QU>9ii-78TMDuw-$7`OMsXK{=n;1 zb+Rg%#$Gmf?7YuD(whUW(6mNX@8N>N*g-oX3@344{X%?doEub99({@G=Yc0BOA{O| zyuek-`}tf2AD9i_Sb6uHAL#K4qz>a^1b`sW{%f_HAc&VeY=12+1gl3+ynMAF1kW!n zYByyH!#EM4*)3xcXtcWTE3+s9wIgrbYobJ`4IaZmG1*$!#qMVQ!%Kf zFFE?SM+_R+^_(XTiNkGn7lMnJxvltM$(4CE3DCiOT*eb6pwDe34rj86N$j}R&kbM{ zs1rg}8Vz!iVBOC8y+L0RE_*UuQ#>FEl|en~4Q7~#=2gLn4@g3rPZL3>fF#t^Tk)t7 zOTxr1V%i635@5h7c*;;v0xn84H+|0)hi9w)4w+Koz#}1d=6F%H=F6321P8$8hW zfh;r0mJ=kJTF*vYVu2vFAg}fcddSdV(Mc2C1>0-Om}r3r<}${m1YWM8M{y^^!XJM_ zssZuuUVrVuwgc}!+X3aV)t$8ZYdaX4IWqp?^7Wmx!UeW^6ZQYS9jJZ~eN?6YZ(99| zmNb>)|8F~BGYI*!_Md4nXtCw{10LJeowb<<4;2qN{=GKsx*VyK#J_9*ik49F5G4P! zDJA#%V>t0&cN7h3$+I(l6uDI?ifP+wDbu`rdx-yC=vTCa^8x95LYzi$VrhKr3L2H8 z@HX9x)wdJU%m-TWL4XEkKi?pRLlnoobyvxO2tT&@5-~NLD7l%#bDR!vJ_l$Uv23}y z=laNZCl+u|)HDeTVF$KtRdxJ5++dF%bN@p+FECrlj165904~~)Vkhh-?mjk8yI)=m zp81&NMJHjc&}iF`FI|075L$U(KFUWHzTeF-ynj&vK4>qrKig0S63q^SqAhh;={zz2 zO-}=eUCwJJVAY@Yp3&0B4|OrKE9U#_*#@8#ClvDHy)lqUhc>?mFay`bZ!$L+Fq(8~ z)_sL^3wX7mCp*Dm1uGA`@s>lZ!A)Bx;cOeu22@LqIdx%HE+r{H#&fw3gW;w3`eCfL z7;A^rZ)43CGB)TI50=@%aV}DAEo*zIQ5zMx{nj4R^G>OcQaV6tp>gpeUI&QiTCOK% zasbk?4>CW7FcOaaW@lHtJ){oS#aN5lLxH&-h0jYn=#O#CDKW7FF6W1oKDZ@Y$VsK~ zYItD_yFz_0FWk2U$sFRap(nQBWY1hI-f0W7w;2*kzS@EXOH}8_JGO8Pl|I}{W(&pp z^`h7)4}%l^yxFy=gK(nlRGL(nHGI?LvH3u530GXsyK>XDd~=EYlG-gHTMh!e=dyT&dKiQ7?(*4UGDY3h^R%?z}#p|d0#OSJN{T+RMd zu1W+s@7&Y>(O!Z0$u50{Yk%O%!h;T-1Tdve!MR85~zTv8xf0eAiqd%{LF@IhV zmJ~4%^${Pz3Nbawos0C<-H~}{boqY-1LyAYh`NcLi?Et5p3tcM`yv5x)K8aP|NSYy zqW|+Q|6k)1WBC~p;S!Yhn^WVrBkgSBp}f6c5IZ6(X0aVYP5GwmGS(C5!LbVsUY4`y z)Sc3y)1`}8SxznPAh?dEh{AXIh~U9jW?zlM$AnOLE2z>vnG_z<+(>-;mJ%8iw$9|S z(!f-{%%$lgbYPZyA#UH3Jy1KVo}IeJ2>I2I8RErRaIjCk>QRg?D=?NHP>I)KhwjU# z6|zM*LGQME#`*JHFz-uSRc^)&3Kz)|*|;zezf3!qYZwDf@bHYO@9{$9i>Af&LwpdS zl57=E%MaoXbq!AMG56DV#-A1a1t9%|P^!0rAY3+onIyI-2wfJ{7yVudf!8LLgnJB5 z7%qHY^?Hv<&XVo#S666ZsEN;QtwxD`z%!>4c%E4dxVASe(UusrIM|=rKPV1#TQ^>W z)<{6Xy)YRqOGy}<85v0|lZ2a0wNKj`Ls%FV_aG{!iW z`vsq*;M{m=t&fT{-0(-m9=M0naNm!_P?cT=IHqZh%mQU#7N7a3H71_*dSAn0YAXvx zrV%#c?_@!^w~vrhTn^41Obpb&D+ld@lB1*4@-Rs!_u=prc?gu24dz^thaoeUg1GYv z5GnP+h-^s#$du zjby7T+=&pTvcdTIMo;+MRO!^9`O#*@valNHHLeYI$*4gr?G($RtQux)SvRXIss{G) zBe*_!el?Kd{XpZ6{Q>jR+5K~tYG6(}uj`kf2ED<%G!9U!L+05Vb$*!0?zn&V>7gGO zqKxBZnAlbahXKp|bQtB~Q`uBF>q`t4nhnxyQc#7wsI1l&d}W9ezTfncR{<&fBb0jR1=K~wbsBljMHW*lYj2XSsK}|Dl%`}`DcozVSDxW(oozLTc^nSG;$5?2s@7q^q9@AWv+YjJ-S zzBOlO{3TG_ziIj_ilwP3?4)TAA&$fTt1tdiF-qGr({S{1Ia(vMlsG6|gCsew^c>-P zk9eiDsyVp(IbsaL=DnojKnfMbjV|`Q4J+&T4p8b=KT^r<)YX;H(;6M_0JvAf|+h|p>Fsa;@pkF)!u>@E!Y^Qz{&PX#hp zG<+thFp28s_0q^7N;n>8D}BoMM4hW@LHNo-;#lImDzKOJTe&T$D!}I95Nt~DNiohFx$+qUYL$OS|C*m)3^}A!^7@jmXALjZ;IW z<%VOUBMm5pe#<)CNdr`ZVgB59w4hkBcWk4Z7I0SvD)l|-fLUfPOBH*Nc*@#wZeJWd z$PhRi-VxdZj#agLG|KjXP>Rj*c~b`XYWBADYYPK>k)OiPv}1%066GIzS{Xs|E}fmR zA`@^Zy!){?n+djTH*l=-jLZ;BsK1-{CNu2X-0OK{4-4G#bl1$tWdVuD2Oi{`?*-bI z%m%`hdtt$zIh3qpFAQ3q{%#t+7fhZH@dnhez~$RZhV7WX>M6S|E8{Ija1l$r`cY&L zP}5Cb_OYRbCB;FS;BYEP{(gNR^%4o-k404-c#a35%lvMoix?HbEVDzoUS$T|l_=eJ z7&m~}>tu)q-Zvp_ChLfke$RJ=TKE)ucJlErSA{M=2kOTc_;&Jq~)Sh}7!FFqO^}vQsqNKYZ#tI?#EvS$b^_oqQQ- zSDC(uLUNS%YkgTkr4=S?DX%vWQF3WBH5|$p@#>TG*;TjnIK3V>*j39kn<;R?o38dBB@{w3Kg^TXNvzsMmkUwqf z+hvHkE!5+msw5~jhrkFX4LU;1IF0lP--_4)$PNsepw&D8rw!?q*|1l`&gb(BiU;7_ zZ8quaoChFG_ePo94|5Qa9KxH#B!mYiudW$BGY5wN|7^2RbEqMbw9B$Lhm}vnw}vt8 zq(j%`#mAdwpf%|u(Dv30#{BP@DB_~bU~nWlU+0h+;03)l;@~oar*d75b)BY=60=BB z=4lGezLxilcAElR&<*ygFo7^1zLCKI6DSMS_||yX1YQKCfTOhu^vrgU8Tgq%?j3&d zlkZL7huZNT27ObIJ{WwEZ_*U3s#kBMU?$OD0<10|O*2S%&}!rUjLi(DpD{~0Ihw-l zczZ(k5o4^(`Mg?MVh9Wj8MNBNdf;s|)XH&P8{WU&+#2e_0HMS2Np<=vaF*&|&3%m2 z+B$8a#6~3!562IeG#2t=7Ig%Cbv8^Oz)k>Bg@>U(`2-vj>fcm!OX$tE?H_?wZ}}f36LkA#;6FA<~IR-c#H(S^bQ( zg|khGdU_GHx7zqX=~o;Y(&f(Yk{dv}JltniG6&Hu@uhHS-eEL7O-4BFl6TPAsF7A>D&)lzsihiuG< z+c+=|x*@r8PZ4&e=hc-bQ>ZT?6RPrO<<^) z9odRi*LkyAItF8^~7}0X^=zT9qVz}_sDeHkc zF%-MesmC55h8yD{vp0?rgF|ZnEnHRvF?32zJWPE^3}bTKWzDt3kZ2}?J2p%V-e=|I zHn)i(v}(L;l$iuxpK|9Vlpq0n8rLKX6%vd>x{p3mnFMa!y%-;=Ndm*DWAus=3ACLj z<)u(0fg-V%B_^PO#NLbKKGvn(O#&-tPxiVCkib(QUE8BZ0-a>hrQf_rV6(JY_}e2ApqpJh$u&m; zJOp0+33{ZMaI2_v<`F4o-o_zUyO#|1S>|uXCy>G8hO${{d~*0G9@@!zmK^quGf(Qy zl0zU4H+-|lp8~QEUzxzCr-Uovwa4lUDBGJJPk0mAcM(ivdA4((W)BKmHtoq5u<+>==f0zCh#U5)l|1fRC ziEwU=veNR$M$kaoXXNZOfkyMmUyl>aAYOazTST#QsAyoW*xO?fz1BZ*FZt#Qx{&sw z&6Z{Z>7757N8kMexf+Tmyts-Fm-O<~PL&hFA;LHr86grlb-%fdyp9asor>R(Wm9>9 z3pBSSb1(MrK$E`f41gm{b9d{1^NUNml-WtGODVXww>@|S$p?>7EY7MxM4ZNDnDA|9g$nlQF|8UMsnP0+^6yQEjH3B9$cB!)GZ+hZFm?a^9I zFcOGJkH9r(!n?z@gv%Y8&=qdeMLn(w#C{~kPYJX@G>agUT0{#T9}Q!FeOL<))J5@` z-PVGnHP(w!-?YG~ih#G>S{ojXx@doI*M=Xb=|awV=>W0ATp1phE@U%s-fG17FJoL; z%g5{VKo#fN?wo-+e$dRad<-xEr_4DL9D}hTgmwDGhp!ugt(iis_EjSoU5iyUXg30{ zOI;q3i$?Hv;_~hL{YKztV8LvSIkz6mJ1#>dW(1w*%r8-98^X3o_!B`+LyR#*>ZJa_ z05*(GqV9kJbXS{rz9BXMlC2y&;z51bSNGxgRFOW!c6ll{Xz7EroVQECKr*Iwxi7IS zb5|F%bgFrXrFDQsoBA`AjTXp-N182BX~48*PAlo48bo#mdfD12LQ?P)8EdyR^r$M9 z#bJh<>Klo(=VSO_!Aa09@zP${iz7N*i52a`9IP+LyEm~1Mebh!E+!Kb`KJ6Qr!Zq; zljz6)Lvi#dj&lvf?Bvv!hS8`PkH6D3nN|8|9^Jn=^(*>6{+!gGQ(LbjiH<<-r@@-w z{G~tFIu%K?wDU(?l!Tykk#NQ+(p~LT$c3jbZC=1?qx4Tmc zu9+sUR!2!g#LFD@-X>`vdeYZ-EEO|(XRa4i!o%o}@)KcKtcX>&^R}$j357QBB ziNhmU3U|OxPT?*;Bx$}R^vkfC_L0qS0u&K&&cg$$sQoTmHWNpCOx=SWyU9D(?Qy3 z1G(%pEhHt3Xx{aqg~^-O3jMK(2CS^!=6-E7kVZ81ayo+sf(+fUxw_cH-V?iPa`$OK z<(k-}Ij)=r7Mdum{jjIM+YARazT(ruZUbkHgHp8c!)Mq@;u0++oWIDe_?{NDy-QQ~ z%F)5T&etBo>2%=L5#d-QOb>^{=MEJ#(}RUNP13p0JurIG=3()925i^kxJ;YGhz)7m zuT+RkJ zWcu#5tJu#ATNH1$rH%PuTruyyk8p%V zu`{%r@Gpz_YHWarrR=|&$KTOk0%^dY(*=T(*kF>BpW`HjOaDN@-G5B{0Iy+*&0#R? zmSF@53wyPt%#qE?<$n|GSG43~G1f7`8oxp!bvwyBz*hE2s4n$?$*Z?Oc&_Kq?5i)s zu|qLm|9Wn%rp!0m^RV#WwSPrR?%QHoHEd%bEkS?QuJQe7Bd7Os?Y<afAk{h^wRqE@1TC4^P<6 zQJJy8CQhrJsRtWgmzy_|cYz1&^pYrwv-v@wI>BM4QwWA#+k`AIBYu|nP4bQh;*j8c zOTh08wh^4J3LZ1w4^OpDogM$Q9~yks9=uxL4|L?NZswi)AyUAfYR(B8UmYkU3MG>8 z$(SNfqCx^1gTIs)^kekh49O2$YZIa{WTExk?uRgNU#w3k#+sx(LD7QOq6Hv0?Uq%6 zAU|Y0&r(R~<^|5wg6hCB%tKM}O^8Q5Hz>NiRl7XJ1!hx^z9a}>6F*1XKA?3D&|Scj zlFMa>uyYHKlWwp<67S)K+(1^y@Dg#9zOfgm8*O4}a2N`s-lF8e*~knbyEyJY+|3L( zBV5QxT$$iX#s%4{BaCoO_ecL_FGf)I>kGNh$q2aX4mOVm89?`yw;|z61_*1DeH>rP z05WU&2ghp}U_8h`%%hb7<{xxVWv(!Q!#zIkyF84TUr_J2Nef2sH4AJu3}*z<*SIcR zS{)X!Fud16!~`eX(ytR6GC{>X_ct%%ncz#!o5pAi$9Sc%cGW_S8QLF~sik4FP5XX$ z)#>fSylXF87IC~_fh#`r^|a1=p`hHRT}hM`zHIks9pPaEMsCAtLvD6BGN&JY0VA<~ zYqt|Ha^Qe{jQK%gn2w<3@j1KDNK6>eavttqPiBLZ&<3;OM|@>(C=u;?)zAN zV*Mn(7}ftzey*OmvH6P3cPBsHb{|}(75OVa-;mxQntosZ@7lkj|07SNydmrnWQ zEthCcZ_Zk72G#cN%_8lTv`oyX>vD_&#oSyc?Nu- zGuzDnVNC=N-|0PB5hH_ybWE}oM+MXP{+Vn2G!VZKYb}b+$(gF>^sNYB0=BkOUM!ic z(AY*K#VO7O@7tK-#|C(CkXwBoy5{+ze(l&e=^e^qazzL1*qn^`lbEDx~*V1jl=IxM;w$ zx%3&qx4QgKt*-!5#g7$ohiZTW|Dxkf6-{9La4^HY4W|isI-d=54`{(uSWBE`nHF>h z1ewgzXhZA^8}X_CZP(5%ax`+HjlW0*MWt z4lpj=3QLgJ0oQ|PKGoLI_rAb-)-;0l2*_)2(Ikn3O@{1tVWZja53t z2GeIJkih8+w0cmSF}jS4bJBwj!DF2xPxQcGFkSy0jXso8y>@tXQ6Iu6dze{Q^}$y@ zDeS@(1GpVWT3p9q2(Jye3tqf1gv5u(5AvSHYVhb0p*C&IRm-5Hzh1=zo~+y^Yq2wh z1F9|tg=S{hq>{Khr3+?Y$->tb>}?KPE^ukF&>Y^pWA=#G$uP(GX6rpMm@(Odt5H$C z7tLUM__?dyRTJ#uwbs;6K+bZsWCg{7KrOC9w+DG@^g;@p3!2uOo zg;z%;KyISXfuElbGzuGuGqH3h?s8spOo&i*J?D4a${UhkEEbS>uCcXnz#Wf7I? z>QydAt)h=Ha<#88tCpi*^7NM?@!_C=ZjP2X1_ZrG-emJ9hp5aB=gFwuAe2vkVOy+4{U_u(q3e4!@MYk4ybvW;9a^Ia00p&RV*cU4V z3WHY)!fqnSt((1wmyelRvT+@!eTHdDCSo!f9wNBCD=hf9Kjyl5HI_QW7<187eA<4T z7J>ciYg7ZN>X3)S4>(PR`G`~>vf6M{1B!mVpGJvgxW?gFyXva;B-gK zRkD}LK&I`S8;-RSr+kL08V@MJZ19|DRSL%TbYZAaT~~m!y5?`qEERy~)y78CYfK1F zRdvi0qbw-N6{M4440JDJ<~y`%a+p!+*VMLAS&-r@yBU!n3)~)Fvf381;BV}@D9b4e zfp)`8Wur1sp!~`oci@E#xY6eppHIMsyLfb7WjiT@b$2w!r&VMi{zjAdJcSG#)v@Sa zZkGo0MiN!Z91Qd?Vs7dP!f2}oH^}%;N`rV%g5@qRX&@Og_F#*U2K|Quaq*b>S+=U& zHA6-jklpKGS#2N#MzQ#F$(SpYQZuQG>ZUZ6PF z7jW&^9T{=Jdp22hab6hi84nYma^lBO>0mdv_grxOuo^2dB@4K2J;W`=(L(GiAyeA| z5+F)D20a*jN<7@dc|mvrOV4cr3~6GI{9~T*-&|BQC#Y6?|FVjIM}HD1c6z4Ogt+`x@k<#A?Z__o9K+r4uc%94WN7`?C<;MN@AznDQER&Q zxh0il)Dc%dtJd>JVQcKd59jcl04);Nzbm=i9*MWyi>Gl&HE$LY3cisimGM*LPZiOK`&A60WTM9+Q*4GVx=MlDAuBK zF2l!+NiWf8h|rnI10hGD$eTxUP=KlzwwdH$d_pFr`Hn2yZ|-;ds)P;Vu8(uHU6uiY zJT9$cjWS@KY7og=DFcP}u}c@;%fRiE20RE=g2u_%@o1;iURF!H$@1c zCBAX&iXyyq(w{>SilFwYS7kU}5x$Swgp9^10%e?~+rb-(;9vhr2B#392%p6>qfeeu zgv;6u4+f7a!dNpdt;SFhm}t%SwaY0&aa!oJA3};?6S#|%5fjM?tGPGWP$@!vPFGzg zks=5+^omX3DMH8Fwy{-W?7f3;SqHfy+-xw|YGhIbQE@>$Y6(U7G=82t&|VRyVtR|x z(iGvyDHa?b6J~s2AtX0?_qGxoI@Jv4bd`Zoad2d%R~ho2Q|MDNsX(b@ZTA6;9lpQJ zz0nKX3C{^6&od}1frbuKg#kw5xE{22J&8sRx)k#*X~Z#8MH1H~$yEuc4dC>+>nRL0 zGEy}b(!3BWr|;PXteA*)IA?u>7I<(Beg~9}5<_~%-TvJ=O9(%$WH$xb7cBisf1Qp$ z(V~n$#D)JRiKF$qR;N4K&h6yj8cX!PDAT`lP=itfpIG_d6YN*?XVGa{g82h}A}r~Q zCw%*p#>AHOZs(xvf6C*b9Q$+FGMaVJstOakv2&UI!f!9k$^CT*SI(b&D%r33&t-l` zF@A)J>d$N&qiwD>%z`?g@_O*t?|&{DGC+5 zGclfBf~15KkXTU}x;rJ)ag((Q73yilfL{$d_0o^#g6TWtA9=%!W1#_^zM0nkak>c| z@)QgaztD;pL+{Q}p6o#7@q{_PCqJVDrtUZH`S;=wm&VbDv9G=&&ydc8_(=mu&-=kt z&fy_s`+4GOukZ+(%(TB6gUvK~@%|+(?b0~9M#{-^W_S`6zzW^i+BDMT@?w^i_>QoN zM0-9I&!R0hk53<~=8#FjqN-2lJZdF=y{eA4h#G@KlWZuLP-&oPuN3z(^5~en=0%Hj ztyb!?YY%O&pdN7_;zMn#=ujl*(Utf$gf8Sy(@Cx)FU{`j66NdYdd7BAuE++`;;?-8 zAZ`PFd9?V1wRQs;O&F7E_imtI8~O<2i4CNi%DKg}w1EzX&J;OpZXkusF`W76270B= z>_^wUf%c>-9O7!-K!xFI1>ufe{LGnqHh zZN-Zlr?obb1n$P*nW#-PuD7xfKC_8bPX<35blXCatwGbDseYgjiS&5R)jtrU^NR;f zQQIhiE&Q{&3?A^=K3aI3iU;pA2F9eg@xfR~XnXrMK5QuWv*Off2;i%?j0sCPW+k~1 zdN&J?5MoINxt-4vLV;e-zWRAWXr8|Odi6LF)Tf@s(;Xy&mk%^X8FYx@Y54)3>?~}m z>5H_qWFiuH_<%&N&58upxOEZ^Mq&ch!)LT&ib#MVBY@PpnFLm4{fK_FlEA{GjfF%L z2{a1K;_5!tk-%iY$CvK)Byi&NJhl(59}=KyGT@VWKmyi> zlDG9|NZ|7B4cks2g_gFvO}6QzkXV9aP2T!KijAL6)4x|i3N(=%Z#A4qA$Cn9+ad#V zDbuMw*M#i_5&LqrQ*sHxx1G;)3At_$e@$hz*Z z9Ju%2-le^yz4tDqI=A+aCZVZ~6eUC+O?x$jCK^U4v=>pLQfa5HrP3BnzY}l1U4H&} z-|M~K%Q??e=RV`UKg09tualJD#u1&dQ{~|u1Lz7<^4Bhr4n(I&Y_M8ii_-TR71Y?5 zpw3tr6FqbPH+__2eOE#Ar$^S)eFN~f+H{-vn`=$$>Hbrq?*Vbg2EmqB6IuN?_+S2; z?th}c;ZPVGSw))Hu!Zy^Iev^&PW^Q(=}9_J@La2TZByx-9(fV?EMa{XD(cN|W-K=B zU9ZO4-n5^i{&$u?(LBo;T!_PuajAWOA*355XMSb8{gvbKcV8?$HK~@byWuhZY9#s$ zZFk+jWB-Y+b;0l#uprBmW*|3MlXa6O{8gswB*o52I`8POU zSCB%4)^Sfw9x4cW$aB~JD;*rp)@g31V}hS28}~L;#Z4&=j>2qS81G^%LM{y z_iu2O@?hQOd%Ge^_@E|?{fUU4061OF4qd$<2)pPol`vHZc#aW$o}ciNfG;$V%I)TL_qb?4_6^5dzsGLq+T? z1VLD2BYG&cMF4EQ7I>JH`5{P4@n}scFTDJs=yP9x3se@_(oD(nz+Km&`drMy`O==X z^m8Qi`Q1KRtpY)}+lzH)tt72FYjOG-E^cpa1b6oIK6 z>fTzWqr)t~$98)^aTN<}GBT+RDP;l5j-c0^Z&<)JQ_wf0kp+JFD<6`dVgUuELvnjL zSYcuMj6uILD^yFqI&{;Q6*_j_e{dy<6_{^MJwMRI3Ma#f3JHUYtYA&n*3l!(1`(N0 zjAhN)z&!d1Z_8;mn15~aLF5G+IEgj3-JfHF)UN6GRbuR*cPWvd_#iuI>^aS9U(5~) ztGQms@i*H+dp02~fCHW$z4+XHoCBiOEtgIn-~^9cThFCAPAE~b_auG61uCcX0;eu; z!x3an*i?xdHzkRtBN0VBu*t>R?kuia-aOfKr$&(%T&`GT+34`XZDL+}jnXZUV^aOu z0h@um*DaTfkL89-9=z)N$2h=qPpb8i7pxdA*tRkv$^>B#ZD#U+(LfFVw3O2}GB7UD zWluW1j2I|B-&4gTaivVhI^y~(?=tn7|raTuaia`>4FF{=X|46mI4Q-(5cHMGj z;jPtM-yjO>8_vo_Z2(Wd@|)?kH9Gol%KeG{-iFq?Zds&xK5t2e*3y_H&TV~~|G5I8 z=63(0O+w*uhd*FfpjXuT?sDJH{5$5KD6XR83lm0`C9Uy)#R*!6wM@vrHXz}h=RSRZ z)`EOCcV#&5=|b-fOfJ`5?nlifLq7|D4I{RN``(<36G(|8PKB{?7Wpess2ON1Ae(>> z_eg7gp-cl#ui6$;C@$*Wm3Wv66l*QwpFX6A7iC*Kod%hp^NL1_>OMBO)nuaom6Hp^ z?)~^$!oGz7O9SLOypQ;SOZc`iOx!I(U#A6D35{&d&X) zkpfyx(%Ah&GQc9p{yMu`4%nl&ISRMQW5=4Z?|hjGkYHgtJ?5(j_XTQ_J9w3VjyF@a zH%kd5mzl;XB$c7d=op&}e&|;BC+{B~Bq#%|dI|$EzY0*ETT*|iiwPI{USwIjRbbh^ zdA1Yx*LAa9qu%XQ0UzI>!)5sYAAhS7;4Gm6#vC`tIoMR-lk<)EPBIlZTCTc1aaI}j zT(59R9ae^OcAJ#=J29bid*f$mtbg1Xt5`SxNf~_mXnOQ(m7!4Yaos}#h6&&6)#6>o z%unrI2RrJru;|^l{7Y6T%5be|^P=HRW#~vQ&UTBHZ%Tl-zDD56_tW zVhTnjvHj-09dBhtV2mt1XWMN)c;(xdcR`9Bt{%6JkbFxAB5bm|_}>yiDWj5{Nfp2owA;i2i3_aZk^b0W|xg*>5XS4 z{gXcbMb`{8_d4m~jvRjR*WMN%{rO=GF^-+TYtlD`n8F7mYV7CGR-+o*V}U=Ba<_i^ z?UEJr=(3eM)f54k6UT=*Xgfp))(+!t>#@a*poWyqBxfcZ6vq){hKJliY0W4pEz z_`mrXtrx5cx+R11CwMeKyRea!94lLcVxHV87r@XZa={WqUAUNZLCGOdAEu949bZkv z01D^1-GNVxq1%M?J6nn=SgPE6?d`DxXl3W!2Y#6YYusdyvMP44r4&!pXtjXu&P0Y* zZ;Z?wh%sHr`i4y(Ok9F}y{zDC?)R&|TCvB1xy2Zh7ACef+_RpIvW5$G(pJx~mqK(N zk&_51HiTV?v%Sh<15ctnm4E%hap10oia*BM6ff@?eY=x~bcYa*i-IiF_0Nz?YRcekffQ(>?+XS~!1N;(KqBiq41pBN# zSIpWCAu{tovieU$&ev)p#l zfM}m5+3oO=|49%V{nk&yuWYHAy!BQ8KR-^ik~Dj5SHj=hWKaFM0qH#(SN+ukkiFyD zzpMUF^!MElUnd`s5+;Y=ei@a0k4`K1)OA?bqm)&?Q_+sx){okvS_fNW+1# zbcs9TQn+`lZcE)L1$$DH+^H(0prZFztVEdNOHZQ7;a?7qRYtNl{A zrsJ=7mqZ#A*{JSt3QNPt4y(P#2*%Rj)ozwX7bFdH=11;YXJV_g!BwfY9%;xYoHtu$ zl!1znNgH`n8Q>s^8n|^;1{yQSEC#T2+ir2X7TpCI$RR7(MuF+4(}ejw3r@0Nw}WH* zkNdJ9p`LcD=c_E3@s}}g)sq9JyOTlvSj;$)#@2LnMGnRzTkG+)R~{N%=L%M%593s$V+~`UOvz4$3}XKM$TZ&XAWM+i*oIQ_S9)^wRF%%=i;sLul;C zJ3>z>LQ1%DkMz_#=}~n0#bfh3`ZK6o>DQg9_X~*H&84k6kOa)a{f_UtftR+nb)4Q# z7H}LM?jn@q*7}rJ-fDsnJdIZ1xhEtAvQlq`qIN36(~BSJBWBdV<9)b(wud$ZoH97- zP_iBRx3s)WKV$$EDR)O(XiOkF$M5)SuN?%q=}*RWO3VW0mCyOTGQ^fTMB07iBvznD zmO$8SYXy(80>}<)TY;LzYo(aCmaxBH?2Frq1?;h^4?aG-6AqxdgLe0KfVJ4+9|ToX zTpFTr_}py-lRVCwh#wfhE+gxn;0Jmjc|GZX-t2a0FxzWdbRHioOm33 zHV`eu9M^z*%X_WV;hlxj0kL0dP;7MJRBovn2z&2-pcsQK*{QjLzMW8mIJG7%6IV6Z zJ3}?eFa7xAvNODRfcWXE+66_o>Bw-M-?(8GU{OX{_a?} zk2-ukv&~lMnL1?lrHcgOw#vQzLB`MZw!+j-uQU3|TcK|yGHqdbD|GIk2#`Cd0k1LG+>AdikR7@WHf^4D$+@Qm7i;Hu2q!sp;HY?9!*gc@Zi_~($D{y` z3sxOHT@3L2SjmAY>g`aiT{cQdvK{Ig4z2FH0#GPcqXzK^ZjilGMb zNJ+Xj^l!7%4Ft5-XXwieHT(M5M}@-KA!6Az zUV!?GWI^tg$- zeO7=7q)k?~Gf9xazJg)VL(yg9(E-BphH7!Szt zs2c0>31F!b&O=@{!l34Qkjdhk7y+8w_;$yBk^=Mj`j>)RRCz@q+`9}$+sO1GRIHE7iYw$aAQ-9p;XHIgZMP7Xm%j2szmpn^};hup2o z+o9m%XObWmUD$Tj&odrlgl}^-naL9LpwI8e%_s(ac-ORr&g~Rdyytf8&wKh_AJW2W zPv*{JB!!r#H4%jY92jEo5hpi*)a#vrZwK{Z|7gV}?jn77@TE4D8~?pPvyj%0>iR&P zZmW7?L=WupgZRth^kDAAsZAZGdQd03yWDkA7mQMm_xWe)!ddkTr26){AR834Ka)Wh zex@3}G9y%BEq*%Q&9)KSp+)0t>yX}d5V@4xGC+aln{@VCJ+1-hj+C(!js=L7eL%4F z0Qj|ShhL5+0P!bvUr8*|axfmv+b<;2xb9|FUMI=U%XENVuy=N^|HSsejZWR zlzp9f*>V(yqA- z++@W+!dw8zC6NbB96$Iyup2MT|xb$?&HZ~Rv++rs4*80Q+ zi;@c~L0D5pCfSeNv4#U0&%VFouEqhU78=P2vLD#N=)~>&-uKwS)UCbuYWZ3XG%7I zxCj?C&Pmu>ka2>nD$|6j9y`STH1uvDVueGRM?IO1m_V%aJ?&CF9e4@{OT8D(!ki(! zM}BrciJ;O(@k?p@4@ATrewWK-60vl&(J9~VM{gy4X#0u zVyp9P(yJ7>)}9_3*pzJQ?BGV`A?~p6BR_=bftQ4|zzO-!-$Sx)go1t(;N{e}i1ZvM&&9s}O0N zgG7dJExJ`0_36Q>2Ba9vzyEqf6B_&)7#UN3o=h9ik)#5;)g7L?Fg;oDa9{`AG~<+n<-dJx&65ZYVK77$ z(eE#E8NpIghV3ddCQ{Fz5aRa0oPdim&jJcD+2F=@S8WDn@Hi52ID8*71Sx2&UW{gj zYBhGz)G}sR9rq?aKgtX!{zPjZ4GS=Ie4ctN#sd6F_daE6u)xz4`j3ltEKoCjTJF|i z7A%NeM9|)Lg9WaCjNd|0zyd!*T4MV8Ss-XwMd&y^D{LW&eroT%eAwW!#IrH#QZ{VL z*LM302RqQSIcMitvO{jn1$%kS33&9z$DL^M3mYVxJ92!MXM^3uvZHfMtYDN=pnvcT zGZdRW(>WQ;2vrUzd1-9u;C!DVk%Ap1hE3=UJ#;4l#`*jm9qvnr-iBG|q{b|&qV$rR z92r5YW{Yy#6kX^Rd6hTg`)X8P>_H#qQ~D| zquj-r;<@k!Okv4|)rP4|^MA+w6UEz;u90FM6BU^)Y2FP2l{N*wA*g5_K8v>}FCjyn zwy3w!+daeA<{`^gc8uxUxxPC^b7*|cW!YehRdT^D+x!1f4gQO+AwR70u%h=5(!6jN zlBBifh*-Gs`X<<082o7G2i>(zaQ@?h#)QG%l!!VrL=5&*J#kIs#FD$`^Di>T zV^v8hsw?}JWI@$b>xpnYHfhW`l0$h{0m_OG=@ws7gh+~}k&AnkKs@=(N(ikoOowUp zgrq5hVX$p&jg$)T1}xI$Zbc&^>9Gp-CdDWj*K@93(U-+WUT0tE!Be^|n zpE@*M&NFTyM5@DuFdn*y9CdiAuKVzMzdC3P1TmBGV1#z=>i`{_t?>F9ncIUn+-qnh ziEwD&3e9$6={2Gn@Z|y#Vf2^=>|{%@zEP!tQNO#KxwmM-6ZVxrJ0Dz%*xJFuRH_Nr zgrmib*tcN~ZL%=#kZn*Cx3y&}mI5ZfESgBT5RKhmPh~V-*3$;hYZs~4#9@?ooc`#7 zp`E0mZYmEdb%4UW*Y_qt2ja-af4Drwsxt>UcimCYf?Go`O-C@2^YcTMo&HklQ29f= z_91Tn==N=2iBMI5kQfrdY)L7IwAXi1#OnW17lK~soaTW&4b1Kq=b7P3*<)7z)qCWC z-cI(IugoD%;cK(N3z}N)f_Yf!e!Q_Qvv~?9!GyFR>hKYU zPP&)MFLJPN@;>__*$Po`pnu-{HBAie*rW+&V|#+AUy1IUn#JKM>-nNiEo>P4ep@cv z6$vO4de)u!P6FPJj2UH)Nx);j{(XZ#B;fEcuM+=vj4ZW3W_G__0>aD|d&wS1K*6w5 z`%YI0c-BaMw3`G=Ht9~wyc)|8hZe@uL)B*DV5F&b`r(ooFuz(*sVx!!*a4NFvoxAf2OfniuS$+jULa%=Ws~`Z?H`SQ;_(^+6rf`yQ5^EO| zmwq^W^QJB+43QmeP1l3#;sz?g3HtC!Ftki!84!%|i(msHIrIJ_Xp_C?UeR<^d(=-*3Z;83OfE5)l_9d!5 z;KRbDmU{VzcJjhmuXn5miz=Dm@pOJdr42DKwbY!XeK?E{$o~XP@f}-aVSGIYHYgU( z-G7~ID)yIzD(MrzQQ^CxMpWG4tz_-+>)#ys6J5Jw@Az}a9wQ-`r`>11+4c(GZRyj- zr7KX5tHaO4(htb>w8NgaJL}OgL+fVil18M*RwD4DvK6t08Ru>j>p~4GU$$?F??d)_ zOJ}5=2NCW3!-fzs%sz`7tR~+wh4ec(jL#R~+rW)Q%kU#VQP!hZ*9I=Fq8%ZY_(g!ULE0b2d9onpMez1$QbL`sGduC?Wc` zr#M*xk{=X1GtG#D;hTx$S)*d0#FKwk=N>^6=)&{a=p;no>%1GN<5G3|9=$Jb*#zOo ztvb!qDf~b>;`7TnmKScaofa$I#sg=w417#=xgqM#J}b{~Y;kV$y_Z6m1LBkJkLeI& z4#h{Ot8TPx&_Q7(!nKDLtn7Z+TeYxY$(QC#)6*Sf061UE{IzcK^84V%cpWoC$?iV^SUW&x^>{a1A~Szz3q z@~n{q3;3&hcQ5#}z=xsIyL>4uuyXFW=gm(naKSY#rPoZLk>nR*Pj)fC?`%50w1!%=2v>euQB%Ju5*HSKQ$X9_%`)C)MA4hCwHGwKgtGM?ixxg zVg1KwvIB3sr`TX+wTn=1Fca5ON79mo^uL*0oufVep5jPOwoh+tb) zwI}6-3vvX*%p^|GfAK}*unQOLs<39L*5QV$myZ+B*@}Ly@v2Qb+ zpqC;Z?c@fdo#*L@%D5ot#JtniMGmN0Z5$-K!v+OM3P&|7nW2!VAUaxk6PV^!bQQnH z9OyI&v*kx5_$L;)Qy1gp3v;RGU+Ybyqsn3Kk*Nd77Y4ti_%xxll=|~0Z`%j5^gkyK zzRJ;aF0X~xu5z6&;#G z*n=@C80XSUny0=>{F2-W18UI1EsOPuHqcY6BV`8(0aw>2qIl+Mck}Rus%mjiOY&2# z;eXTmPju~mXY(hmtx5B)A`0y_Q>Mk?|3UjIRkLoYxU7An+9&sFb5C)uPn5G&6VK-@8Zzh zl`GitG;?#yzA0Jk@3@WzuLyu^B^_~=pnuGLqrW+AH96Qq@n`dmpp4I z!%8%7Ap5&50kQ(l9E^ zmPU^|!@`67=SPZh1NTh+_Q+r`;%$q@ zO(4GO!QqunilAE-Y=0dsWjYwq~@S{R21ep+Fl2Z^Mi5P%iwNCPOy01+aLRU6F8cjXI#8Y0s>T;vqVp3 z@*3)>4A)ktMAY9sqRH9QXyBOh+Uoq_Gw!3-O|iZ@KTxNq9LTNsA62dJ3wx14=fA7- zzi8g|B=TF-PPl8M|5rn}vtLKJZN*`2?C;gPN%-bs*2iwuB}-Wx*ih&yJol2Bwq5<7 zmH%IKZU3=ahqjWE5kv+Xoqfdaq1U5JoI)B8k#|(V?T<#;=uOhWq(k+&=-Bwy5@qWW zR2kGZ+?rI5cBZwB-sgFThBxPlHmld5crV|5%1iaActolvk+})o-m<_udAJoZSkiti zMxRlr=S6#hN)KX@3IAr((T9|z6`5Rk2hf{5%W>p~Yh6~1*W(pGjUbLdH#R=$F?8y> zT+qPPZ|L*<9hJZblgPvEdYx(06e{%1i`Be0jb_h3x>TV(gItN8o!Ind24%dxong;8 zi-u2M^NzBaMS3hnANOIuEo)1Y5n}UM)LHqaHc)96-8PrnJ+g{z3r{~zC~BBNx7^O; zf4fbXK{P*H4$j-mAeR7l`e@P_^zlN5TI<7U)V=ort&ZO`I^9R#(ToN?Gnc=rZh)Mgn#Lp%x!M zk^pc0n-gZ#q)bV7!{4iev?=a|w{%R8qt`@*j48 zEoLW&kW~L2Y^h|xZxd}C|CAJZHmOQ~)h2=CNjr*JkM1P`-I=Ok60v2}`!;oWv&;`v zW#XDw4Kv6}_b%@{nK4ASdHbnP_EAsPi`@GRuwQyJB(XcHJ|0eFA=C-rJwkFFkbF0VE&x8wmsA{s(;?w!5b%gl9+2tA${I3tJ=uFL;i{8 zO?@FJB_%Y*bF=0?$VMi&1O^kn;gnZO{7P0-j0h#WjfsrQk;~97(x2S#P^Yt~!qM;= zRCmZ?sWPwuwbegO$`@)zhiBf^oXBoRn{>k`b~5*%?~i+V&rwo@93#p)( z`^d+gzvy6h<6=463q~+9n>A&uWd@-OaZbTl@c7)aR`Y!uc2L=)s26jG1B1+ilqGyQ zL4@v0@I)%M2z>VQ19w+n{ejcFFW0 z3rPC8X&PN;0pj*}t9`gQto1bC*tCKLK5{EQqr&KR(upl+44PR$1{@Yl%2~kplY^sN zItvhS-mmOFi>uDJZ!KE6u>kdtvppI5EEsHVNJ_XWh#e8`_&=vqVgcPLGyXgS7BF?? z*>B*=0Qn=Qn!0Il5-;(}g0(W_n_G9FeiAiT7Vvt)(or=PRrZ?b}b(3znl&G=U` z9Lj6R%m$XibF~Q=u76uM0XDy7gV;AaL~^ip&m;>^@iiT&$6YnZd5`%wCvNm_ zF2mL4lf(uWN0w0Q)=Ss+Ek#)sc8LUPc*Nj ziZ~u?=%D!ve+%dvD$W&~FRdsO{pA{Nl}jvXjM?x=De${@<)qmE8~aZ*?`b}1GbzEb zPdVg~cRC7@Bx&bJ%|wnq-z-JGWTVTIL-%&2<-c0W z(XQaUxz}$hQ6|APCP=vk4VDemn;&RE73Cwm&(1a@YhkMB49^bKw{$872qpF7B|*7RO4hgcwK`I zuD!~>#M4R-_P*bIo8&hE(w8gC6~IzdwsK}Y7Z{+(w48*V1j2WPZ`W*n4p#4m9>|b2`V(V7st~w0Y~lJ(aJGKxH6t%@uGqeDt2`W zkHs>=4U5fLURb=xKHx*=svskF&nXbOo4pCf3R~EZF>ZqAVe~V09t`k)AUBMxmL5bJ z!`meYYV_bW(>~0Qi1Vw}(dD)g9dv)8md~f51Bt^bj*i{5&{-PYTF^@io1d$R-kzid zX(5vEJQ6yH8-8|vo`DX^dzsF(a?`=P`W>&yrLl}j{DQC(D;+%F)mp8HRUO+5TZ^t0 z(SpEIP6}ltEwFs*Jy5QR|GSIWm@tM#q#xay%*T*qmZGR3=_w8S6`?*TGwI8n7DQ*JdI@1Fj#M z_%8U-z}4@^T#S2Zz*hM6boLQi_%@Y$)02!2qr>;EVV&^S72-bY|O45oB zZcUwNR>U=2QO*e8XdN2(RNKyc-G&M>V^i;ctR)AF!}3L;RU}|A+}mA1_X{-_d~djY z?I(&C${)Dxh7G&RI*p9$zoPns4CDEx1h36{f0YJcxW^{I({J=QgG3M06CRG@0+f9F=# zznj9JD1NJ38?8auK9al`rh7}_xK_oe$vV6K(R+sa1*$Sq*S5HOlvjQ{enqza(bI>b zZ!e5)5b-IHVztQ%3;TDTKhd?mJQykB50B}67cP((tdX9HlmC*Q|L3%X>SNo>?H=!5 zpQrak1pVAb4q^WCMS^A2zh#;KqH8}Q%yq!d^=nWAGD6YLY1Z!gR3zW#o!xkLV2Bt5&=qw;+fQc$MzruC>r zkw;IZIuJJ?4hiaX#=b`MntCAFf}<6kR^?a>dfS11Q5a0b4|gN`jyENBU;EIKlwiu8 zk^urznyb)VAsInj2kM4sBFE4}Kbw2YhLgy~_vhw-?F@Ux4kC?zV7VV!)^pJkhx9M>o!FOiXP@(KKsbwj_HHl&(p|(=;Nxq zPcb<(R+8N0E+7Y=CrGj)ksNlnwZy3UlEc|!#;zZ*5Q>+htwY8aGTeD66(vq4gRoqM z51ZV{;Kh4py)khzSaP_ms!U7_PtI9p#sZnAKqb}TjrUBYb024*I$h`)`T1cXC8H2!e4C4U&Tdu zlpO4pTOSa7?J1zp?5g;oR}`>g@5~^P7A5rEi+yGFjuJAx*pDb#Q9)VU7pb6bD)>n_ z!*c!xHLxugi<2AE0RKzwJhwp_`1MG7tTvVw7`m7aL{iX!&EeyVH$Ca#J6|?$a~d5G zO-JR=WY9sjG*!mbO*&AezVzAn038%bvk>e@Sm@wuzs>H}YFZ$lk`-sYO$!8-lO3D{ zT2LW3k5@ZI3r)rut^>BT;FC%3-Y`Z3s?Ae#Q-}szBt@yxD^WvZgMo14JxVD0=1qMs z8UM->a#-r*NFm$#Wa#-_L@?n)+)Yuwh+O2IAHJ!XLurRTzd4;ffg*ahIB>OMW}H_0n&1AROm@YCHd zGLFH2Z6HTlGTBSt{~hR0G|$h86mJlVH%`oP=hq-I>GxjqVNIxM@!8`j2;)^khPYpLN;-mxK764=bp_C^U4($5}x z|7C`m7OwgdyPTzD1j-2#p{u#fu#MnHj;^plj)U0KR)0W|4BA)cR7zA$-CHFTLhY5<_{L%q&xG_W?rWY;&?uDjT9qA;YTOw`Gc(^27j5xeY zcu^7(-o6R5IE>-htgn?0JZh1G$2q#D39`~q@kEOJ=NU}R2=*D)Y{sQOo4)7tJTg#1 zrE^kXuM9+6TVuoFWS!@2|YA{wxbNKYZ%D zh2%iou&npkQ8{ppRS^&RBnM~hRb0yjL%_yi@@i+&14}m{EYSusIPcDMhF{Atgx2J)j6X zDcm_NvJ}CHW;ygSnG$H8I&Rc#sRWN0juYa~C_#?*O4epCB^W8F_tK?T0kSxw$jP!j4-G=&GYH~C>iYW-@oG`?@weeckxX^`w$YoYS~0sti}2H3s)h4M}KXl z|GrA}RJJfzW(uq=^D$WWvx@BS!;zZqx`{d&@iN>7ATo);TN#dg-K zCd^Za5)08fkj|sKBJ#%HvA8=eW&L@k7;Fq&*!Y>ukQ#I^%N*eQvI)Aju(&EcVu7kl zJKd}jI6#QelDZ*tGmy@UFur`mM*yNJ9?hXg0wC2hmNIQ91pCRA6XGvo>(Q1?vM1a` zpwxahz0*e#Al&IK@u?7j?~feMXqbyY!Yz3(o*-dpS?W+Da=|}mcgg7exWo5t7jutt z7C-#_E*sk%%?nnPy;E*qH^Xkq$-OC|T<~mc*(1dj%RXK(J5_ld*9~q;w0gVovBD_T zh5W0yW7ixKF!vlQ(9)R%YLqK6fxwOfjLXgplft)gHXhnyZH`Sq zeZkRQ^eF=rTeipx1~5P#4U6EYJp(wX83txJGJtedZHVb126!@77k=&u1Nb%UQ0T?F z+wZ9epY9NLY=Ye0mP06a6YOM8=`Q48gkkq{&+YdN*V|FIAr15I8fWdjn zR*b zupy@5InEAO0)8kGu#nu9xNrK_!5lyzc{=^A1a1KaRR+9W=7edB?2@P0B)VGL`KIME zH@w^XWO9CZGu(}hbE(xp z^W*8ELxbq48A^a3O&w9ycSIIU z_QqFz5lPI_HXT0q3%OD7P{`Ghz~Yw5cZZCzO||77DlTyvU>iC2oc=fi6wdjTG^8>S zKuh=Fu6h|(kSxv}NSk1X1I*X_Yd&y6{wt4vh$XJmwm!d>QL*5{rzQ#cA1iV{ zpOJ+0(3~r6t&%YD#v*u+23Dhukk&0tmxAo@-iuz`(jY_rgg|*q8sg)Mu2kW(ZDw0- z7UN^;O*bM8pAIq}Pf7U1F2F|hbCQ2*k_v53c&_^Hl)W*oWaeb}7J)l0m7!X>OCNg1!#)7x;^~U5FY_fc_6L&Uk9Iusq zKR<$)1iVcAr0db|1L4mti3Z=DgT*YC}0H zFLkn5lIYz(S@d6YjfI$}bA`BN4SUW{cqs07`U4A=N59*4{$S`Ak|J>z70?PUO0UDZu` z#(Y(U>rY=d>u<(V3o6%xJ|H!qIo;4W?xY3`WvLv^zG|?0^UEdv2uy&iW9=nGURDE& z0QcGQST%5pYn@PwR|6uR%U>qb@BpC$$D6a&K$~99Kobi{+SCQabG4~~{hP8%V@h?9 zaI70BG*^c^0?a(mQq|$(gZ=3v^jop{wqBi+!B%)oF`?s&mF(UqQan+rQ->k_eoaOE z_nJ>hN*a`^g3nD4vnL(jmEmmM1pDZ7Md-d#XC%cg4{UpRwD(<*gzHV4rR_$9f#uRG zmlzX1h_uRQn<8O{6wRaLLiu!HM$*O4S4spG3_W$SaZ_uA#)d!O6P%agj|SHB^7EZC zrE~1mYk4{APd}8_^4Dr>isqH!(JvcLpgb<+FA0%k|K`r0XkKt3nGz`>R3P7RjxHTV zwW?H>JbQ>_gX$?6lAoab?7i%dU9u5}ugy(2)0ZfiuCH4BTpm)}MAf$#SB$P$b+q2H zeT|xoP44Jje~a3)$KQP}u0np-%in&zQj1*f{X8t-)_^X*aaUEeZ9+#9PGyBQx1vzx zBL%%y7*gEYOvN+Vjmxb=5j*6h`%s@GOO#cAKk|6pJTz`LfVL5*2DfDoqL3>yg2A-I z=))v;=`6Nz3COsveu(cY5^Z2Bl@uIB{5>pfPr}Af%FWyXTE1`S6W5+++*FgOMTh2; za>f+4?9+4JPd4)7onO4wZfaQB>{oGw zFh~tc-*0KXw55Uj)m`DA3TeP4=k6s_MqFOa;?Q(8riHwsk58N3X(9gOc+;=bw6H^# zPrfmL7DVQE~>D+Ho)U?f|xB-j{p!&}(|I6f!0dQ9q}L$h)hL?ON#J zz?NdESb7Es>5pUgxyS&s&6`~6BsYOgPTzCm%uUb{-anqywF&&%N8Vc8-2}rSW`gf} z8Q`X12dii>J?w}VFaEHP7JlxU)Ly{;Lt8E7D2Q?WDd<$ss>LN@Fx)A_b|_;B@$Vu8 zOJ>d?|D7UqqR(f?(ZRd7w+&YHqaQKfJwl{fkbhECOV)u3A6wrYkM;ZgZ_f}(_FhS5rOCOHL_?)1A}VbWqIGMb zQV0_s;4$&Uv2a*~s1G{qr6j zpZMpl9VHuwmy^ORNlAM5{X1oy3Fr84itNzvfus8`cDCEU0oFveCl9-VU^U~tJ6lpe zfc<8^!+N_S;05#J-+>&lAh10yukU_544t;GMS5v3HA!Y!-)yje=shMzZRL%1iTQ)p13(1c^lY`qtUNv+*+n)=L z=B^Lq3-iHvxT__B;RiHCYOQb9`w6|pr#Roo6+^(BTl*?T+`3+2>}%`cGT6yxluR%$ zhZUAReA(F*V6|Jh^3_x&JpJSu=&f1}b5HdT6GLhs_WiS|>zTDc^wI{f2s!v z-r7R>(Z7JF=SFG;Zy%gJ*1p5UvLDQk^{Ybd6< znF8bW=64KSr{G(xOf|7;3VcI<$nncfgSeTGOMv|}oCsrO-BLFVTN#=a@5;=;g1&A+ z;Jq0bT-g0_9Ft$pmoyQK9kM zCf3YRVGfKrbJ?GqnFEh4wdYBAP_{(bAC8jlvh*q~GJ zqtXIgE;`DnkUkHdVXaKGN9G|<$wF~g=^U8Eg{&1hGzS`{x&j@)aHrS%yMzpG%>pZ( zxbuA;jH{VX_I!YOK;Lqmkx_1(0p;^0?t;lP@Qpp?;)}r<;B(*V>b-6j=RB|X8$6x` z2|{zpiHNdUc<*vm>h;<=80B}4qOqO>s$4mq_wjR}b@a@h^o2P%aV^Q?AlE$PPpvlf zRGtU%{YKmt9Qbe6nIriO^Kj&G^D0eDN_l};Fe}M251ZwC*3NLx!ybv?qhUAa!Edin z>GAn_=#;(9SW>zGKjLkzF5|Q(;oW=vpFhN@(IM)r*eTqq_M^I1aRIJ-Xgyw{5POLg7dJN;{8X6b*R%l7WmgxnCF1-kQTmHb$s}lKSy%PEK8UZhGGzN48=Y zi^UAJVJ7xp8mL`y*x7?vL)zrB1;tlY*54&N7vt>qY#tja>PL`F9mZsKonVdf&l{ z?i<6f;=UvMi{HP2&)o8! zLDzK4#ILdNBLBChdU6JQC~f$rl)Mc;dY;9lp?zHd-LGV0jL{NA*>W_t>)Qx|s3mE0 z^6q0oXs~>*O58^wWP0b%${SfiNQ}*vl@0&;=8VLIgya4r2XAUg+6tk}WTVgq1tFxS zZ85@9CWz{N-0m(o38JUG+n>G`5Jc<^V**oY0!Y2Zr+n^^02-8@8RTi_f2TL%z?S*j*XmMSKzZt9Cu+Me97nV=ZNQk)qE0H{YT)C_LQ2 z$?y6aWWW_6Ts zBV2YTQi=!tE|+84hRYY+H6Q*=#(%rjTA4R&;z2(c2Cfl2aIY6`3C>prP8boR&hbGJ zS9^CJpO9OD>*Sv6afy!Mz8_T&6N=d|PvqUuBMl;_dC?E_?pOh->YMi^rw9&8mUirp z5Z%O&>@HXuIh@0B{o*Rw+tjdQD z%k4QLdQW=|vf)_ObgYOQt-8&ea{|9xmAFX_VQV&Y`}?g+9{4W%W}=TH-iR6*((tVk zx-`!0KCfERk+w)XDNp0|4oNGnF(H^iR@|SFN83K_$*B@6>WH;u{-B!yVKKOXXukwzmch1e42q!CB8di@n#XLYPMM1SnMEHaeb|GAlI zEkc$mJFa@JMH@fGsBOV%8MKF{Fz+@&4xQu(h+Neyhps;iFFdKU4zWFs$aM}|hZ;Jz zJAb8@M`1-cUwTR&8J!8yum3EMxV5!ZTvn_{m5(g@XH0PVrGEO%uaxz;f8WhvOdE^t z-}7zMJfnd0KI!K$hbW*b&kcS}e-x1Ol*iO{%#@Wsuf*5ms)*Vuj;iVqaOIX;{#v<4 zRwZOLGjE}TLsuI28oP$}DxoIf$jgnEN+`^02b<3uCH$MKzl6EMdyO~a@UV?Hc-D@l)-bV^Ch{ZdK zvWTIRTIq_rF;4>DRdX}pKb)xC?}xeSWkz(#=3tHCmKg}wQQ?q&*9(ILzBcxorDO+y zGL1q(`O!h`&$H|PKTMTrRsIHhIm*wNb-&Zb-*U@8<3%Y^ag5f7)R#=I_AS{&C{i?| zr7eyFW`krCQQ~(pzB^X?pL3Bgpjk%PSG!F2l9afuhVHL(L6I`0Xams}8OVI;#9{B@ zRysbhgwPUqkIk@L(V~OIRxS(IXUhvq+c$89=-bPw=>&0 z1yZjJuu5Ab!wIT&ZegZr5PT-vD)n6kY`&8o*>E8X4DBmo|0w6c^P76nR8)EJqQ(C; zvr+*8Lan%Z&iVd;E&Kh}N=p?1?Jvwleys#ts-oU_6_o-9x29aEO*wF7#zt?ItAsnE zPdAy?R>4}K^;D9cHL&9rtCG4z9dLXoo-oO-hivJeOdWE#xIPCd*M&5JW8oSN1MwEH znxondtu3Hc>wG-1q!m7G%-XJpn+?LbGJc2jl6Ek>HzQWt*#WP-97=D|PW%2((V~3Dix(!s@+@J6$HRW0k42*p=uSd;sN81Oy_4W`^{271-6Z(?(-0nNJ)DG1oG+U;tK)yrf8#!L)&$hbe`Rp) z9EUs8F-4nR|AA=~ir>}Bvj{lD z*XWUK6SBEcKt{g~%%Ajn4Pvi9`nkx_~meSGimhYFRw%gk*_NY7fQf0bdobbg}JM--`x_z=f>D!S!BuqEl=^a44cqC`ZE^ruOr4l?f} z$TM!(e#N_{cq?hfgC=>q!~OEfi?mYZQih!{MGz?d^{0|wCWn_V(xQ~Ky=oEBER6*n zky)Nzn&fIwz;`L~CB2((Bv@O#A`kQV2#-9i6y>|nb@1crTaH(WUz8G~${F|462$C? zH=R{K!iRmKU7w00;LAPp$82*kpf6|Y%3zQH2h`wab-)*RT1+(?)|3kI9h~FJ7c;;x zWXCg>^Vv|p-zm>nJ|BKc2R7&*{Q*vJ++X)+5h&U>mu4L-g?(4gy%^Fi2bH#O&$nK$ zg!7s?{nyK@;rjmM^!Gt^1Ymv>XE`_62sZcHoE2+Zz?MbEOpv`Dj@sB)XW4dwLa_wZ zH=7=49AC4hC$bOL?=E{e{d5SZTZ7t7sYl_yGwq`@x_G}RdDC}hd>oFCm+Umf;RK>< zr;IU%=w=QDyx7PVr;XG!Cui?iXM2e8&o8bp77Q4mVbm=XcORaW#(NjwPc~ zPZp#!Xl?Ob&ndw{oK~w=|_b+6>n)8Ma+ZmxDCJVVcZWTL1>%(Ib4Ax zZyQ8IXNW=ihrVpqyFLgxVfO_~OS@oiUy=I9tY(nfthD=2P6eFTE!)$vE+2HxNEdfW z5IKYCcK(fTl4*ab-5to**hn7`V#Pg7vc~UOZ!cYO*iJT+U^gj$bq0L-;%pB2$HU+dPxRqJs~Fh3dHTwz$7h%|?=5A!l?;-B z|7=_%CBw7;_o8R zpuh{>X&mfyNI7xB(v}x@J8GG-8Ck)AJg5ieCTAJYS}OnK(-n+}iH=&U1AmTOUugpP zjhIpK?zT(iJ@$%VFhCMU-zDYoX0smURywpN}gBe{{$yEPrjJc$4zEvm0=`tgY`lweIRG4ua zq0;sgE*uSUW^gWb{Qsj@MAkOp z+7`pXBbWEj&?80u!_Hfddf+2Sn;ltZfSamJkkZErj6pz@YJo!-uyMzLNZeLdp z$98@Dc+|ffoIU)nnRb`M+CM8tc}*)|i_NU=NH(DYL~kDOYT~T~*_Wnn5vG-J>t4%M zG1p4aVLp)C)>#Q(blP`aQL2K*q;1CzI#+?(9_z!l%~jwfa%rMub2ade#Uv_zsD@`W z7IS$_H4tWfD|0-i2C~9WUfgU{3)vrTIdL=MY>431LWk5kP@-lHA#K1S&VCHxOF=i7kt5l^2!J=6<4uQxx9{?ZFOI_QN9So&ZlwPdfHIWF)tN-CvE>jUM@ez_-< z`a$?{eDawG{g6A`)TEo&51h_d{e`&)faQVmVReH6pxNC?_@d%G0B@c~s27C|fOd#Y za#z#}Mx?!MoG(UB7C;0X!cD~hV1DVOBUta5*z%)snM?bh4EG9*KPCYAy z`3Wh7RV!liVRz-G&GIQ};I#9pu~tkROt9RJ9%Xq4CsQBla_oGBBj%ioSACqy>=+_h zN8_)OMA#8N5D$Uz)AOE*ZmAybv@ zArI3SA+)x^%8<0P4 z?UfD*IV`<~>KS3)p=5T_GyINt`JRe}&O7xr6Q4wEXUs)9hUKknm!u@6=szT(9r5I5 z>P__6OA$SkN?sg__uKa6vxhdR_Dkh3s9l@Z{6YHel!F$Y&W)b?@@a6)O|uW4_9Dh?GY5%e2*2mXQ&jf`^k<5 zJ*RH6@N%Lsd$;ZR4O}Qjqu%SP8#hWiBF&|v!-MKi##<)Q^P)GJuF3bJ`Oveg(_#;_ z1kgig%N^-oaauyrsP6}+Zrzr9?{MKE5oDGg)@ZOE2XQ-pocR_?5JhV19_leLilNJ< zGgKQe_iL8s+!-GSG33kVX5R1?BP8{|2pC}yzS@ViYLfU@Y5mh@X%S+`Zuj?%X)nc4 zj*zSS_G`F})}fSS(LG}5tB7}&I;NnOA&Ly!42dGTtf7Kwj6m)mu}BuaC5l9i-|-T3 zMA6Nz`|BqO9HL0hWgiG-Vn{gu_V?9iMUa$=^R9Gp5!Bja!}z{e7&&(PJCsKWBU^*- z+WeTLC>svkD!40*DlQCf;=C-37)v)izkEg*4G9%pw!)3xKlJRn*(ak^~rxR@0c5wJ5`=+0u*S*%8SnFlHwc4`@`Fd z>!0)dpa|B+mYrdl2wP4#?n{q^I+3okTIb^6SoYgo zhb@UPWzi%U8J2{pB_6%^(oBPs0qdFQhBF|-o^M29Z8lKJm{!qk%LVKH8`PO=^9fL| zRb+Fmy8uk|stcG7|A0`vhzYy9g#b2Ile1@jf<{Gl+dAnYaM6;ez`J5;Q z^M>vrM$cl{EU==-zM>eE_4iU=k}84d&*QvwJ4=Au$MIRUMG35=+45-Cxdi4GX!VP} zmO%aWPFbn$643s*M)mv*p#;wJelp!VUjqHk-7h;;N zHiCCW1858wYHVB@c+V47G1&qeYhCinFknRIffMz^eXa1KpWnmOtQE$bzD21$ZG|Yo*3!tN zR+v-XI-&We6~5>zw-;=00}~A$f{gZqHlRJ@$ySiv2H~7T3I$^A@UZ3HIh$+ka8_ZT z(bk{sFf3sEJGUKf6sowXP`5+7Y?=;lc`IzI zk|=$0qy>1t2)>_b#}C+eb*{&!1~|FKen!ewz}o=AiPLwBA>;P?`Ig-Md60G6me%w_ zI;^$cHufbr0eX}?t`jzX1TU{b=T>@e=#ylTKCm8W~Iob z#2#Gt>d0@_WyklD^l#n(hcQLzh;tulZE+zQEF>NzKf!;!B_XR)L*#RBwfBoC4B3ih zA3EByf8%nMi&9d65s@wgLNvsQ%}ln8c+-+uqfgdZxxZ{AY^IxZA^7#m@YU(p4dj`} zWCl}PIpqGA1&WDSwIlE5g=NDmNlES!E}SI|^oEvLsK7W*s-ly+o=6_(KWp=aW=|56>$i#q=1U(o`|w2^Pj z1MxSN|0~@kuNz8bJpJNk`(cRj+pG{K?@=F_y|f#HkJ@EfG~J}9LGFd9Q%=kRIQFl~ zHvyRL?Vuwo^;%Pp_3n}MnEu4EjENL-CBzdQMGwu+a4 z+70M;o>3G;)e7D9*`Xq+apdZqqrqaxXHe{&ki7)*e6eZxke(EJ@Znqqx3mnJIuUO4 zaZVOhd)z9t&Xq&X2U1;So#auLQ@L^Tb_FDzQCW8emt>jb9&A3Du7u)zCpi_Jm65ey>)w0mZLZtgbq_5fQ$;Gnd3* z=aYFt^y|zvp_|?*l~=nqq1;`g6C2K{q1pFwdRuzc5R=b{l8wz~Boed7D;;;gwz4n( znH#$WrF~ib`^cHC=(MsL{~PXYs7}5izrJxB+Q=%AHXXkmZSV*Y*COC{EHCdi6IbV| zqe9+*yoc<760d6AYc2+4^>c%Tzb`_f`8xc`dv>Dgz1+VcekbY=UnA&Ru@jT3MC&?w z??iOZqdzxyBXsu2PvU4YAWzP#C0&Z@=%9iEj45wLi4*)HE%_VK-e2E02I77Mvc(6C z1A5n@r@zC-Bwp)?BR}PP?n;%ow1>LYC%=>x#Zw*3c&0i7UOcaV4 zuZb@AbL!nRQ$#i_knKtTnK8rI8;V^fM!NGtRNArSLoG`G<^`}TLZl^X^D#7%j7i45 zko-UmZH`tp7fkA0(hmW{L|B@X)Hn9nKfzX&q@)BPVk?`P%DojIQW z(zXM06K6#IDc*0CT0?l7>HnvDQTo^UQ2F;9Cy4!3)W=9mFJ4MR9;;Q#i7lk=ZveQ> zl=}{wdat)p12vyy)?I)6sA{-Otd|#nApE;WtODGLMb{@$uW$J1(Fp}5F-+w$n-r^ zxZ&-fNQV`L(D^pE^(|MlD2)#^BY_2=xb)wIDBhl1n4IDT>wg5u1XDeM+G&crZDimt zkN4K>%T+H*iL>gAf9E*%U-KafFkNE3bhk*Sr+3UAMn3UK6NQf@vXRRbElTmL!cBbs z`s#$YCz*@_%lv_Ty>{GOa3I*wJUpS5{|>%3J0E)e_5&1bFr{vo42K{+?Iw5Q7%<&x z>bGh>9^QWT`9`z#3uuSGGMKGPg`Wpqp5-#?sGNf9;ri}-$EU$|AJ5k3jI;2N>X=Xbia9U{_?GE7NSMQg46F0w3g_We znKZ*X3@JTiaV>OVEfsQ0NZ8MMj|%y(vK}n^NrjBhubYn+#EsTB>aB-k)JOo;OCEYk zjSkei4Bt(lM(m!FNt(^ns9gB^R6g#X;LEBu?!JZw#o6h&?7@s`ocofk8#HK8dHIfM zdCaVK`a}q`CgB#Qb>KYYf8aF@+PXHqKRuHM4H`e4kej4Iv604SW2G@^!t49@g%8sr zcKYe4g?qH9#$fN^{y5CS@I7sOO)V{IJSgzNXMq+81U)ko!%2XF4a2;J+;phsMC?{I z5jymB!CHO43>{Lq&TIHZl@2{(&ij0N1szJ;6iy`UpQ1&9`=SNMt7*|Sleu+_t+eP{ z_Uyi?L0UA-NYu>aqC?#L0UBw$=@4_XpT!3=I<(NOrObi3QJFzA;@tG6k$|4{qNz`vd1JJ2&=B4ns%6u>6VX z3%}sC%SF26+E(CNh&PXKtb*@F8Y)dT-{JXM*54@aE1b^r(B{YpBscbfX#*e|DFw>s z+fVU3?E)&~t?}?GHOp@~6oI9P8{1#f30_`i6H8KT6Zook$ckwMki(bXY9pVbw z3R>#DfDpTL*Hv#{sMx6Ccj)C?$cpzjNI4x0wwKOy$9xQhncjQU+#C_GVA=okQFsg# zD7di@X5xX(Lq+ZIAQ2q99;|UoOMyY$y{4Hb1CGD>E$}oq3)T#7W2LXlg(JkxkM=bc z0L{4wpB<}y62RBVwCL?XF*K#;tkKddg5B%54Fc0CfgpL@=p}G==bt3GVX!})r1UFz7#Lh4)4BLY zpgKdvhwI@8@V|TBNk|wU0Y%O_g}GyxOOpNd%xTPa32V+r2JQb1fo#(EjH-Wwu*Lhk z!M8DIb872MF{s5q@A2VZ zV^HqN6S$6H6lgsrG6_Wkuwq|UbdgmzD3)Hmp??$COWZyHS{5}>sqyYZTtYE?(NsLP z!8H?@Gy!YRl$mxOpM~begD-SKsZjq@gRGw>v`8>KYp(6{Cz2==Ih1D>6}%s7Aml)8vdavJ3f9P|G8W zkw1^r&eg3&_C3dTtny!nCNCd7cJBOoB(a+Mg|?C+Dl|VVZSJ6iE**DP(WUn(#v#31 zkhR~s8DHmZNYTcGSO4%1MQe%mj} ztAKchHg08-l|&tGH#HyX@uGbr4R&X#nNfC{3GWVZ7e^3%Tc3qyfrYO1mRG9qfZJ*1s4{AwD6y?Q@6YF>^)KT)jWpUFWORSkBzPBqgSH;=~2M4zE+FChMj+yLya9 zq#tY)^}arYJ&{jagloHj$8D5hhwbX$%SLzn&KS>TW3Y@T^&SLof&Lh`I?f}HQaevL;TIp<*FB@qtXG?@fHVP9BCL50!Xc)lBVKj`xcRdOq>+!bPAtXm_=|Ir{wZBt4q?q@gy?Z&Vl%5Gb1|# z4i!kv{Y%6({s&EH%dZPAKp)8KoM;kto3NWA5f`>?+lV7_UbZGn(nXJtGcCkCkpUG`?BM&t$j0+^~Dc0`eL0{9~likZ7X;DfvqFH1B{oO1} zUM>ocI6@`w!f@gjvTl71^@iz2iscIZ{3GpZ#Or0cm!!lF(g$m)OS@?^!! z@7)aKy|1t`@QH^C$d|S;6fIY^C?&1PzJIkLn+Y|0d)kYeGr)dCa&?(!HoQ99{xW1B z4~F!w$D|4u!iUoi_lzw|Kp=qatUFUPnlV93fH$up~i#ED@3td=Ht*D>kw zoUaYC4mU3d^K^oJ^w{^_^F0t$TJ^r9p%2Ezj}25748et4=4W+TFthlPBh6_C#t6_W zvd>;PY8-ytEK}w&GywsBu%^E43Wlr0%*1_x9IZNJV5=&(&@f;@WFaW_5%c zXZyHi)Hkk1`-E=HXyFDP?a^=7JY9{6l!SqJknY_8}-b=KZGk%RcobDQSn_x707VUg6J}NGcr>Xs-G<^*;CbQ5I{`aa1-lzstma7kL=XLLbXgQ;K^< zNKM$$oX>{Khgp>VU23v-Jwr>7HUA)ycl8MkO7z_&=>rxlCN}17*|NnMvaS zklaX=5YP_-5%sNdzXIL^-^LSOCuTl^dr6#_3`;azdTF|9Cr>^rk=)y?UI|4n$`@RFmutNKJ}O)uSOaK-x*~i^@1bSlaWp{fuQJ{AlSw`+6E-@u8T5Hm(EKW{5|z>2yNhXSeyo^PO;IXZmqfhb}m3*1@Pt+YOqB z?8Rr_cEdyA#8;iVJ)rt*&G4cEUlD^bmSq-~w>^F4{y*+I#w3%Clfqq(oXQ+Z-24xHc4q4Aj715> z?IFZ9=lAHb<(I7`DbeH+4qs#R*@(LXSvN%;T~v@cH)l~M+3?Tn9v)+Jp)`D{2Yx5t zOfNV5MJeezVL*!B+0qh}X}+DjIPe;-b#BbCqsqteM2NPbQWnxo1{tcl*zK%ou$8VZ(_&u+#C|jidc7tKo`>cj zr}=E)`sP>lwICN(ZqD8>I+hPn*0TZ$rawTdQe#D10uHG;_-UUheN+M&GXa}oFs>pj zz#y7OuL5W{ggOe*Rl{(}``QDJ))&Lijsj zz&`x(%A{`KznDC@+MpLo?sFDwox`Yd+WsKs@*(&%RlB3tWCU8WB71%%jKXf#M;qlN z2$(j4<%pHI#UJ2bvn{7JXB@8j+{hCbn1n`ZB>(Bq6s*!*NM2=!*&Kp%s$1M$_tP+ zoH=)3(*iuoXo-@%umHnDs=wI@mlj|j`mOxKY5}@Sq-aLH7a-BasNjbC0uCVQL}iC8 zfZ*wO!QcHCpxI=5hrx>lI3e%VoP~dd>}3Ko&Mtrq-AcJsy#>&d*!Xe7F060Ffp2il z0<3rC{o*gR0E|87C;Abdh$Ppk_8ERTvWIKLMAcbG8I~vCj3w# zn+j!jz3}Yh!ZG`Ji~S42xLsdc1=ogo{O=CQR6Kf2g;IjsYvr%tfz1nohxzA#|F>AG zgupc1Y01u{eliYklgqZdL}SkRC$z6OJ{$mb1Lk*uv|X@=d2*-Oj%LuT*!Dr`LJeNS`CsIkBC(p7<1T-&NS*xcueW%n6R?Bg z#vig$WPSAT@*Dr6lz69$`Y&0oBxBJ(v-*!4uFCuLK~eRU-b%yNeB$rdn=?$;1ue6o zCF$QYI=)2Df{qkYBR+~@6egwb7;BG^)l%#PsU@49-MYNSJM}Y8HU+=g0Gu&q%hWDO zv1#JKI`-vkX^A35G{OugG5$2@W4i6^ZWgZ;thnwdj!#vK9LYuA`?~Roc&~=*VCXXBUT+>uv zvsZ^2(bnC(?G{Xj+HT~EtygD4pXn3t+{L8+SCc(sVWvb5h160~0nIYL4X|z^C z6ul0(Vs3gv7}-gqo%L~osHifk{(!Uq`ttGXcgZ|H^zxOjpw<&!#My6ToBv`B5(_l{ z-SwRZxxDZ^sHw|?q7Hld$6y9FZ93UWGs4tr6kc+<`N(D-WZi3$E|@bghLQ9dfm0WeM=BoUQ3*>Jt2tvk2=X6 ze!fBqW`Y8XFL?dXQ|pHKM?XH~&-GbYhBmT9CIB_G}20w_{|gnk|Q zGpysY%wW6SOZwPr?`dY28~m^rk> zZ*Q22NO3yZP}=`voQAeA!t6JG^46lZQ~!(VH;Sa#yK`4k+K5k!(O z-<%s+cXO4qnx*4W(lfjhe{htawjnl9T zGx6XcSvSRI5Mkqz`DVV-d*lR;bOlWmp|sG^Wr~)hq_~s;2_0~kmT=`x;yT&051^a2 zK(PGcBXFd*z3mq62wyJQ+7Ydsp#RZ>QNFcLAp3FX)>9nMKz3Y1-3>j5!`in3MT^`( zem3XN?&BWtC|~K`ey>+B#JpnLOpPzRXG!Arl6(tv9gZXB4)5S*Q;Mf>LI^-sm7ir| z7~HKoe4v9diU8+2CO^lQ#R7q|{DAR^1o$-SA)uB*gy6LEdyh${fc!0;pWLPCP_vOK z?$x_2s8q;Gs-MjThCMf2+UPJi)IBa`^;8kGyH!|OSCxWisYu=(TFa1cEXqCeQ_7=E7nQV%G?dft7*Ef z;NRA&x16Q8Xv2cxazKhkH)N>%dBSkE8*b<5#+KafhRwQLiS&26f%EOwQ}cx9-5}<= zR!csl8{RUFDWv11k@3#x&R53WaFLGwz)k*caOzGv#3GKq~H>&w|!V9{x$EL43 zLB}CvyQfnpcwL%5OxX6Z6Iw;R9&JeKg#I*(@P{&8Q18KL$nDhy2AR8RM9jKjW#|LT zbeSGN*N)pyfA4|4yAr0OKm5WgM15%c#$NdL{fvrjc`t1Hv(nV>Mjz-+SV`}T?gKH} zC%@o`b9uGFx3`UGS zgMnvPinF7F)4#0uh+YrbP@1kj6RIq4BpG3F(_VpV3Tb2`blc&j&A>!q;6W#6svM6K zm)p*w6x$AF?rC%L%8sTBih%cArID55PN*GXS5O{^skLo=x9F7lRQTL~m^XQhoS z@0qwHB`PsswhqjHA$DJk{0I#rYZQi~^hv8#Qp~0i6!_1}2&{j1dZjDHY6V&)a7_j+ zPc1A-F;GkM-(u|{zDlG&Puh5;o~n{ZqFAh|bMSJZ{}%ZvD%Hx_tw$*qD^QH9$H1U? z*?U=%5{0yJSO(u=#GC#_aVegn{K$5#{?0Gn#RqYqD3CRs?TzpFaz%?$qPiG%y6{qu z-T$|PQ-qq|&$}c%K6Z>e$s^{g=xn7Zf;w*&eznW4>sh9INlM(Gvbcj&d@U3orU!@m zjqp7;UggsBiTSdv*@(OS-gN?V^;J;g1h5Ql$;5`={tXyUAydYfroIuhq@>=?8A2FDbq< zKefVjW4D&4o|mMgkNF7cQj4pd=sRc!Nin@{ZwK-)6lbL8jQ-iy*Mi>cy|`J0FS{qN zxKG-$`7TL`?njvZ&R6PP1zAy5iJ8I9%YR9pOYhwmi97f&(XrqC!blX#`wwCK%5%fD z30sY2x|gJX17jE)%|u7M`G@fZ>0VFe`LFnbg7E3%A7QCYri0$1GmLyfUIizdg02WG zQ@JE14mi_!Gs}5({=H=5v=Yt=v>qxU&z!;%o;2?Inh+P@9rO}YS$bIuaNIs~ce&0* z>Ea$i&XZgqK0UznH(W{)F`)38-sYgUNfUG4r;af3iFsY0P%E={TdsFeO43Ul!Skgh zOhog$9xJhiGjh*NZoj$-VJSlCA9mk@19utY=&#)e0%7%v=<$aj6&-so)5a0#5FJ60 z>lp~%_HB~pbc1AxTlXzjdcX?Mx{ND=ph-~MRfF9R0t^mzts8p_Zn3iC38}$g!n0~! zxM3KGY7LB^coz-svVUHdZNMQxRpa7$FJB_?tF5|p_IV1>vfig@mCt|zRWGKE?OCAY z8?wnvHV+;hojsu|_Z_}5?__&`%S6Oo?)Db1ErIVXIxY9(N}(m?lh-4ma_Bs5`nd=f zn$j+Cz4Cfd30Y?z^_IC;;jU+;AsP|Y@MNv&DU%;HpeUox=n+e(g|*%7Yaa*H!Or4g znNFR0usCo^zk8w{PARu$3Whd7Q>IioJMI!+!QsyI!J!f8)^t|IyEa0AXT|)$d?T1! z3gy#MHGyUMfSX*)8 zU+;YhD{Jb&m!WhOt9u>zd_Ebv_hB9U6p+1YzNZfUbY-rPX{d!P-lFp1O||&)AU{!x z%SSkX_Q|QfYWQ6c9esRjHOMsN#SvujHC^29vmhfbwG~@&sQA~Ovg)sZ2-o{H=bu#o+G)R+e?Om=-TDF$B=r$a0aKVW5cz*Do4JZQ30+h8zPQf#_qY*rm>?T4h3>p(3ngdw z)Of=!?Cg8GsK&|L%`a25BwaLMyf40@B}Sd6pC=vi6q^EtPDj-*#ar1|l24vhv*qFY z-6`h#@_-BXmKOPC!z@XOjloQR%Sfby&Xea`qnFXFO4UR<FJ7pv5E7+aqOYmGp~#n#oK@d%l+ayNE+r93YUnEhFk4Nl9Wiysm>Ft@$JEiI)L)J|6NJ+HLE&=9G{#+lF#<3}v= zpAEFbl;dF|p2s>MZ^+oytFi;MPo%`VOyh#h#@hNwLoVF>e3yUmA{X8*F!f5T<-%W2 z<>ii`Rr8>$>GkiEJK?Vy4MFvXPDth_cc+;1Ky7Oi(e5e_4w5>!E!8}T*k^p|S~m}F zB&6H44)dUik-ugQfy#r->hAoCOddp&a?hKe<3VtX_98bg9=LlOeA?#00}FN{E?^}O z1UslfeR@2QU1I*IZaxo?JiQfb_?^H#akDXs-3h~^bu?a6C$#5Q^lr@UgoU3!7Q04w zf}G!Khql5_G_eXhRVTv(S|`?`ES3k?9k((GrV3q&^wGDG&Tt)m2QYHWbQ-8_(4$&>cFIf+Y%KrE;A^9q8TuNp#W)j8vWv=d zTHgb8ywcv#^$mPEsd544SXwI>_^LrbemZ25?qfz!I4dWEqz5ugd*0{u;gDG z`U71yV(Ajm))n2=2frxDk;i6XHT2>&Qe*d@=194X2o=`fwg`jcQdxrieR*T@<3w0+ zpk_0(CmB==HB8QhrUD&fPu_Ja6PEh9X5D{-HhHLwGk9*E3x%Lvelp zyT2Jony|d95rkB7LPc^FG-E%GLB83iQIAXD04v?&l%G=@{5fEluF7nOJLyY;dXI78 zl80mPCS4xL+Y=V_=5)cj!o+hFlP_S-(jjd?dLYzD-x)LGE5v*KwA_N}gYSoD+{8uo z!7ar|&+z^};2lz`x7yJUrY?*JN3Np3(IK_ZTKnO1U{JIw+PQ7Xit~&)gwO$4r*!l{ z!@vOeuGoAa@AM!jGu&}?X5ZjPaGl$j+z_M`87&Xi9R@Pdqwxd!i1r-#9;2=yfa`T} zjO0Ji6gOO#LwGM89#nZnxDVW zi6jH952GhL-DVPz9Uj;+VHAv8Ihhtsb|-qRXX}*eXiGLcj7FAm$%3i^dFho%gieSR z2ds~fSYPbXiIi|McaaEKsfEodTBG*j4^a6=dGiG?krGDJFP7S6-IDqM5sEClfTjn+ zbaSTAHop!&Zt+pqb_f#A1Lc)IUEf|Qie)o{%+FK*lt>REnn_Exs%1;b7NH1mn)HA{4dIc0?H~mmz})eqe@$ z4cc^M!^ER7KF!~l$0-$l!NWuc`{2`e#iL|YT?{tm&?Xa{3yb@obDdPAMiHUH|G`;U z5RH`~Bq^(SFFx`JPC4-DG4=)E8#46=NALt1cVE#scO3z}-qU9a*d=fPhkqP9rwo>d z5X@Dcy@ahfDG7)6Qec%8XACDTv1Fml& zD^Pydz=L;Cv0$&|(RLcFUZ>owlU+>!9NHz0uRw>B_aA&YJ3v>X1&)KY_W=orOW(;wNcCsQwFEt z3p2yK!mr5S55+Du>T{OEw;y6HS+aXNK93|bZ0(Q5>nW(K-z%GmUuqIc@!YM57u^0{ z%hOcG7pfAvYzR3j_>gAKzRXE=iT?dCjTl*`s@l^A)mQ$RW*5-_nEU!^<32Bt z_chr5UlnjDWV*W`qEq$hqrH8#p74y&4!y;+Fx&;47p_)FX6(^+lCqjjRuiZpJ zKnC%rVS0t9`#9nM&zS;|C^tm)cLl0MbNF7ParGq4i3V;qrJb9kh7qBN-Er*16)=hA zT4I7BJwK^37MhV9icVW+o|w(}dPv$`MCinoL)2TI${^dQSe9_NY$4(+Hj0E(remw- z?)Br4H@`yz03Rc-J+c~arBoB4LJzbPdf?wONE#b!k3kKsZ~eOlCt+yvto;1_uU?4_ zyJd|pTFS}HY+2t|iME0*XmN$5sjH+Ed8`X7VFzg3?+3LbNf<-4E5 zCif1Ys&7}^tCO7aLrSIy6{`NXP^PC0a%l)_gmL;GnG&jBNPSOnIlfq|xfl UN!1;?aA9C6y#$HSz(Aw_0t^efdjJ3c diff --git a/src/multiego/interaction_matrix.py b/src/multiego/interaction_matrix.py index f6ad4195..b6fa2153 100644 --- a/src/multiego/interaction_matrix.py +++ b/src/multiego/interaction_matrix.py @@ -253,6 +253,8 @@ def nonlocal_matrix(self): self.atmat.loc[self.atmat["atom_pair"]=="O_H", "energy"] = 0.45 self.atmat.loc[self.atmat["atom_pair"]=="O_N", "energy"] = 0.45 self.atmat.loc[self.atmat["atom_pair"]=="O_C", "energy"] = 0.45 + self.atmat.loc[self.atmat["atom_pair"]=="O_O", "energy"] = 0.09 + self.atmat.loc[self.atmat["atom_pair"]=="N_N", "energy"] = 0.09 self.special_nonlocal_dict = self.define_special_nonlocal_dict() diff --git a/src/multiego/type_definitions.py b/src/multiego/type_definitions.py index fa3d9964..ca5873a6 100644 --- a/src/multiego/type_definitions.py +++ b/src/multiego/type_definitions.py @@ -333,10 +333,10 @@ from .interaction_matrix import InteractionMatrix EMAX = 0.15 #maximum epsilon value for the colorbar in the interaction matrix plot -P_TH = None#1.0#0.3 # 0.005 # if P_TH is None it will be chosen in the InteractionMatrix class to have NL-NL repulsive +P_TH = 0.006#None#1.0#0.3 # 0.005 # if P_TH is None it will be chosen in the InteractionMatrix class to have NL-NL repulsive SHOW = False # PKL = "atdhisto_density.pkl" -PKL = "atdhisto_mindist.pkl" +PKL = "atdhisto.pkl" matrix = InteractionMatrix(pkl_file=PKL, emax = EMAX, pth=P_TH, show=SHOW) special_non_local = matrix.special_nonlocal_dict # print(len(special_non_local_A), len(special_non_local))