{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Exercise: Molecular Dynamics with OpenMM – Solution"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"ExecuteTime": {
"end_time": "2020-07-02T08:41:57.017179Z",
"start_time": "2020-07-02T08:41:56.478203Z"
}
},
"outputs": [],
"source": [
"from IPython.core.display import display, HTML\n",
"import matplotlib as mpl\n",
"import matplotlib.pyplot as plt\n",
"\n",
"import simtk.openmm as mm # Main OpenMM functionality\n",
"import simtk.openmm.app as app # Application layer (handy interface)\n",
"import simtk.unit as unit # Unit/quantity handling"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"ExecuteTime": {
"end_time": "2020-07-02T08:41:59.476097Z",
"start_time": "2020-07-02T08:41:59.451380Z"
}
},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"# Jupyter notebook configuration\n",
"display(HTML(\"\"))"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"ExecuteTime": {
"end_time": "2020-07-02T08:41:59.700080Z",
"start_time": "2020-07-02T08:41:59.647628Z"
},
"run_control": {
"marked": true
}
},
"outputs": [],
"source": [
"# matplotlib configuration\n",
"mpl.rc_file(\"../../matplotlibrc\")"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Inspect the system"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Go to the website of the [RCSB protein data bank](https://www.rcsb.org/)\n",
"and download the *solution NMR structure of a small peptide* with the PDB-ID 6a5j as .pdb-textfile.\n",
"Open the .pdb-file with the text editor of your choice\n",
"and identify the section in which the molecular structure is described."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" - __a)__ How many residues and atoms does the molecule have?"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"*The peptide has 13 residues and 260 atoms.*"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"*Ways to find out:* \n",
" - *Inspect the .pdb text file*\n",
" - *Load the .pdb file with a molecule viewer (VMD)*\n",
" - *Load the .pdb file with OpenMM*"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"ExecuteTime": {
"end_time": "2020-06-30T07:33:28.191465Z",
"start_time": "2020-06-30T07:33:27.935342Z"
}
},
"outputs": [
{
"data": {
"text/plain": [
""
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"molecule = app.PDBFile(\"md_openmm/6a5j.pdb\")\n",
"molecule.topology"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" - Parse the .pdb file with Python"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {
"ExecuteTime": {
"end_time": "2020-07-01T07:41:44.043154Z",
"start_time": "2020-07-01T07:41:44.027736Z"
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"The structure contains:\n",
" 13 residues:\n",
" ILE LYS LYS ILE LEU SER LYS ILE LYS LYS LEU LEU LYS\n",
"\n",
" 260 atoms\n",
"\n",
" 20 frames\n",
"\n"
]
}
],
"source": [
"with open(\"md_openmm/6a5j.pdb\") as pdbfile:\n",
" print(\"The structure contains:\")\n",
" atom_count = 0 # How many atoms?\n",
" model_count = 0 # How many frames?\n",
" for line in pdbfile:\n",
" if line.startswith(\"SEQRES\"):\n",
" print(f\" {line.split()[3]} residues:\")\n",
" print(f\" {' '.join(line.split()[4:])}\\n\")\n",
" continue\n",
" if line.startswith(\"MODEL\"):\n",
" model_count += 1\n",
" atom_count = 0\n",
" continue\n",
" if line.startswith(\"ATOM\"):\n",
" atom_count += 1\n",
" continue\n",
" print(f\" {atom_count} atoms\\n\")\n",
" print(f\" {model_count} frames\\n\")"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" - __b)__ Copy the atom definitions for the second amino acid residue and include it in your report."
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {
"ExecuteTime": {
"end_time": "2020-07-01T07:49:27.799754Z",
"start_time": "2020-07-01T07:49:27.776698Z"
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"ATOM 22 N LYS A 2 -7.295 3.130 -0.024 1.00 0.08 N\n",
"ATOM 23 CA LYS A 2 -6.396 2.084 -0.500 1.00 0.07 C\n",
"ATOM 24 C LYS A 2 -6.593 0.793 0.291 1.00 0.08 C\n",
"ATOM 25 O LYS A 2 -6.656 0.809 1.521 1.00 0.09 O\n",
"ATOM 26 CB LYS A 2 -4.934 2.531 -0.384 1.00 0.08 C\n",
"ATOM 27 CG LYS A 2 -4.546 3.653 -1.334 1.00 0.11 C\n",
"ATOM 28 CD LYS A 2 -3.114 4.126 -1.099 1.00 0.12 C\n",
"ATOM 29 CE LYS A 2 -3.033 5.105 0.061 1.00 0.62 C\n",
"ATOM 30 NZ LYS A 2 -3.301 4.453 1.372 1.00 1.15 N\n",
"ATOM 31 H LYS A 2 -6.928 3.882 0.486 1.00 0.19 H\n",
"ATOM 32 HA LYS A 2 -6.625 1.895 -1.537 1.00 0.09 H\n",
"ATOM 33 HB2 LYS A 2 -4.750 2.865 0.626 1.00 0.10 H\n",
"ATOM 34 HB3 LYS A 2 -4.298 1.682 -0.589 1.00 0.10 H\n",
"ATOM 35 HG2 LYS A 2 -4.633 3.296 -2.350 1.00 0.13 H\n",
"ATOM 36 HG3 LYS A 2 -5.218 4.485 -1.184 1.00 0.12 H\n",
"ATOM 37 HD2 LYS A 2 -2.484 3.270 -0.879 1.00 0.56 H\n",
"ATOM 38 HD3 LYS A 2 -2.756 4.613 -1.994 1.00 0.49 H\n",
"ATOM 39 HE2 LYS A 2 -2.043 5.536 0.080 1.00 1.15 H\n",
"ATOM 40 HE3 LYS A 2 -3.759 5.890 -0.097 1.00 1.06 H\n",
"ATOM 41 HZ1 LYS A 2 -2.685 3.624 1.493 1.00 1.67 H\n",
"ATOM 42 HZ2 LYS A 2 -4.292 4.147 1.422 1.00 1.71 H\n",
"ATOM 43 HZ3 LYS A 2 -3.118 5.123 2.147 1.00 1.58 H\n"
]
}
],
"source": [
"print_residue = 2\n",
"with open(\"md_openmm/6a5j.pdb\") as pdbfile:\n",
" for line in pdbfile:\n",
" if not line.startswith(\"ATOM\"):\n",
" continue\n",
" if int(line[22:26]) < print_residue:\n",
" continue\n",
" if int(line[22:26]) > print_residue:\n",
" break\n",
" print(line.strip())"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Open the .pdb-file in VMD.\n",
"\n",
" - __c)__ How many structures are contained in the file?"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"*The .pdb-file contains 20 structures.*"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" - __d)__ Set the background colour to white and remove the axis icon.\n",
"Visualise the peptide structure in a representation that emphasises the secondary structure (e.g. *NewCartoon*, *NewRibbons*, ...) and choose *Secondary Structure* as colouring method.\n",
"Visualise the side chain of lysines in the structure in a representation that shows bonds and atoms (e.g. *Licorice*, *CPK*, ...) and color them by chemical element.\n",
"Render an image of the first structure frame and include it in your\n",
"report."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"![6a5j](md_openmm/6a5j_0.png)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" - __e)__ The structure file contains hydrogens. What is the total charge of the peptide?"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"__Rational approach__ *The peptide has 6 lysine residues in its sequence. All lysine side chains are protonated and the total charge of the molecule is +6*."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"__Forcefield approach__ *Get the total charge of the system from the partial charges of the atoms:*"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"ExecuteTime": {
"end_time": "2020-06-30T07:33:36.621320Z",
"start_time": "2020-06-30T07:33:36.334920Z"
}
},
"outputs": [],
"source": [
"forcefield = app.ForceField(\"amber14/protein.ff14SB.xml\")\n",
"system = forcefield.createSystem(molecule.topology)"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"ExecuteTime": {
"end_time": "2020-06-30T07:33:37.154969Z",
"start_time": "2020-06-30T07:33:37.139706Z"
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"\n",
"\n",
"\n",
"\n",
"\n"
]
}
],
"source": [
"for force in system.getForces():\n",
" print(type(force))\n",
" if isinstance(force, mm.openmm.NonbondedForce):\n",
" # For charges, we are only interested in the non-bonded forces\n",
" nonbonded = force"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"A non-bonded force entry of an atom looks like this:"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {
"ExecuteTime": {
"end_time": "2020-06-30T07:34:09.098321Z",
"start_time": "2020-06-30T07:34:09.084371Z"
}
},
"outputs": [
{
"data": {
"text/plain": [
"[Quantity(value=0.0311, unit=elementary charge),\n",
" Quantity(value=0.3249998523775958, unit=nanometer),\n",
" Quantity(value=0.7112800000000001, unit=kilojoule/mole)]"
]
},
"execution_count": 12,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"nonbonded.getParticleParameters(0)"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {
"ExecuteTime": {
"end_time": "2020-06-30T07:33:38.642730Z",
"start_time": "2020-06-30T07:33:38.638591Z"
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Partial charge 0.0311 e\n",
"Minimum energy distance 0.3249998523775958 nm\n",
"Minimum energy 0.7112800000000001 kJ/mol\n"
]
}
],
"source": [
"force_description = zip(\n",
" [\"Partial charge\", \"Minimum energy distance\", \"Minimum energy\"],\n",
" nonbonded.getParticleParameters(0)\n",
")\n",
"\n",
"for desc, comp in force_description:\n",
" print(f\"{desc:<26} {comp}\")"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {
"ExecuteTime": {
"end_time": "2020-06-30T07:33:40.437126Z",
"start_time": "2020-06-30T07:33:40.429946Z"
}
},
"outputs": [
{
"data": {
"text/plain": [
"5.999999999999995"
]
},
"execution_count": 11,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"total_charge = 0.\n",
"for i in range(nonbonded.getNumParticles()):\n",
" nb_i = clePnonbonded.getPartiarameters(i)\n",
" total_charge += nb_i[0].value_in_unit(unit.elementary_charge)\n",
"total_charge"
]
},
{
"cell_type": "code",
"execution_count": 15,
"metadata": {
"ExecuteTime": {
"end_time": "2020-06-30T07:35:42.178777Z",
"start_time": "2020-06-30T07:35:42.164249Z"
}
},
"outputs": [
{
"data": {
"text/plain": [
"6"
]
},
"execution_count": 15,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"round(total_charge)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Prepare the structure"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Save the first structure included in the .pdb-file to a new separate\n",
".pdb file, e.g. *conf.pdb*. We want to use this frame as\n",
"starting structure for our simulation.\n",
"Use OpenMM to load the structure into memory."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"*Ways to achieve this:* \n",
" - *Copy the ATOM entries for MODEL 1 manually to a new file*\n",
" - *Use File > Save Coordinates dialog in VMD*\n",
" - *Parse the file with Python*"
]
},
{
"cell_type": "code",
"execution_count": 19,
"metadata": {
"ExecuteTime": {
"end_time": "2020-07-01T08:00:13.657571Z",
"start_time": "2020-07-01T08:00:13.638562Z"
}
},
"outputs": [],
"source": [
"with open(\"md_openmm/6a5j.pdb\") as pdbfile:\n",
" with open(\"md_openmm/conf.pdb\", 'w') as newfile:\n",
" model_count = 0\n",
" for line in pdbfile:\n",
" if line.startswith(\"MODEL\"):\n",
" model_count += 1\n",
" if model_count > 1:\n",
" break\n",
" continue\n",
" if not line.startswith(\"ATOM\"):\n",
" continue\n",
" newfile.write(line)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"ExecuteTime": {
"end_time": "2020-07-02T08:42:08.207761Z",
"start_time": "2020-07-02T08:42:08.142822Z"
}
},
"outputs": [],
"source": [
"molecule = app.PDBFile(\"md_openmm/conf.pdb\")"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" - __a)__ Investigate the topology attribute of the loaded molecule. How many covalent bonds are present?"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"ExecuteTime": {
"end_time": "2020-07-02T08:42:54.695050Z",
"start_time": "2020-07-02T08:42:54.665385Z"
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"\n"
]
}
],
"source": [
"print(molecule.topology)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" - __b)__ We want to simulate the molecule using the AMBER 14SB force field in TIP3P water. Create a forcefield object suitable for this purpose."
]
},
{
"cell_type": "code",
"execution_count": 52,
"metadata": {
"ExecuteTime": {
"end_time": "2020-07-01T08:39:29.004956Z",
"start_time": "2020-07-01T08:39:28.795742Z"
}
},
"outputs": [],
"source": [
"forcefield = app.ForceField(\n",
" \"amber14/protein.ff14SB.xml\",\n",
" \"amber14/tip3p.xml\"\n",
")"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" - __c)__ Our structure file does not include solvent and has a non-zero total charge. Put the molecule in a cubic water box with a minimum distance of 1 nm to the box boarders. Make sure that also counter ions are added to neutralise the system."
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [],
"source": [
"model = app.Modeller(molecule.topology, molecule.positions)\n",
"model.addSolvent(\n",
" forcefield,\n",
" padding=1 * unit.nanometer,\n",
" neutralize=True\n",
" )"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" - __d)__ Save the modified structure (peptide plus solvent and ions) to a .pdb-file.\n",
"Load the structure into VMD and visualise the peptide as before. Additionally show ions as van-der-Waals spheres and water molecules within a distance of 2 Å in the *Licorice* representation. __Optional__ Add a transparent surface representation (e.g. *QuickSurf*) for the peptide."
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [],
"source": [
"molecule.positions = model.getPositions()\n",
"molecule.topology = model.getTopology()"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [],
"source": [
"with open(\"md_openmm/solv.pdb\", \"w\") as file_:\n",
" molecule.writeFile(\n",
" molecule.topology, molecule.positions,\n",
" file=file_\n",
" )"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"![solv](md_openmm/solv.png)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Simulation"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Setup"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" - __a)__ Create a system object from the topology of the solvated peptide with the forcefield you setup earlier. Use H-bonds constraints, and PME for nonbonded interactions (cutoff 1 nm)."
]
},
{
"cell_type": "code",
"execution_count": 50,
"metadata": {
"ExecuteTime": {
"end_time": "2020-07-01T08:39:22.211165Z",
"start_time": "2020-07-01T08:39:21.733922Z"
}
},
"outputs": [],
"source": [
"molecule = app.PDBFile(\"md_openmm/solv.pdb\")"
]
},
{
"cell_type": "code",
"execution_count": 53,
"metadata": {
"ExecuteTime": {
"end_time": "2020-07-01T08:39:34.024830Z",
"start_time": "2020-07-01T08:39:33.765468Z"
}
},
"outputs": [],
"source": [
"system = forcefield.createSystem(\n",
" molecule.topology,\n",
" nonbondedMethod=app.PME,\n",
" nonbondedCutoff=1 * unit.nanometer,\n",
" constraints=app.HBonds\n",
" )"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" - __b)__ We want to simulate in the *NVT* ensemble. Add an Andersen thermostat to the system to keep the temperature at 300 K (coupling rate 1/ps)."
]
},
{
"cell_type": "code",
"execution_count": 54,
"metadata": {
"ExecuteTime": {
"end_time": "2020-07-01T08:39:36.297499Z",
"start_time": "2020-07-01T08:39:36.279060Z"
}
},
"outputs": [
{
"data": {
"text/plain": [
"5"
]
},
"execution_count": 54,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"system.addForce(\n",
" mm.AndersenThermostat(300 * unit.kelvin, 1 / unit.picosecond)\n",
" )"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" - __c)__ We want to propagate the simulation using a Verlet integrator in timesteps of 2 fs (we need a contraint tolerance of 0.00001). Setup a simulation object with the solvated topology, the system and the integrator. Use `\"CPU\"` as the simulation platform."
]
},
{
"cell_type": "code",
"execution_count": 55,
"metadata": {
"ExecuteTime": {
"end_time": "2020-07-01T08:39:40.913769Z",
"start_time": "2020-07-01T08:39:40.880457Z"
}
},
"outputs": [],
"source": [
"integrator = mm.VerletIntegrator(2.0*unit.femtoseconds)\n",
"integrator.setConstraintTolerance(0.00001)\n",
"\n",
"simulation = app.Simulation(\n",
" molecule.topology, # Topology\n",
" system, # System\n",
" integrator, # Integrator\n",
" mm.Platform.getPlatformByName('CPU') # Platform\n",
")"
]
},
{
"cell_type": "markdown",
"metadata": {
"ExecuteTime": {
"end_time": "2020-07-01T08:10:45.239892Z",
"start_time": "2020-07-01T08:10:45.225072Z"
}
},
"source": [
"*How to figure out which platform to use?*"
]
},
{
"cell_type": "code",
"execution_count": 38,
"metadata": {
"ExecuteTime": {
"end_time": "2020-07-01T08:14:01.518063Z",
"start_time": "2020-07-01T08:14:01.503579Z"
}
},
"outputs": [
{
"data": {
"text/plain": [
"2"
]
},
"execution_count": 38,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"n_platforms = mm.Platform.getNumPlatforms() # How many platforms are supported?\n",
"n_platforms"
]
},
{
"cell_type": "code",
"execution_count": 39,
"metadata": {
"ExecuteTime": {
"end_time": "2020-07-01T08:14:23.611169Z",
"start_time": "2020-07-01T08:14:23.591267Z"
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Reference\n",
"CPU\n"
]
}
],
"source": [
"for i in range(n_platforms):\n",
" print(mm.Platform.getPlatform(i).getName())"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" - __d)__ Add the current atomic positions of the solvated peptide to the context of the simulation."
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {},
"outputs": [],
"source": [
"simulation.context.setPositions(molecule.positions)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Energy minimisation"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" - __a)__ Minimise the energy of the system. Save the new positions to a .pdb-file."
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {},
"outputs": [],
"source": [
"simulation.minimizeEnergy()"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {},
"outputs": [],
"source": [
"state = simulation.context.getState(getPositions=True)\n",
"molecule.positions = state.getPositions()"
]
},
{
"cell_type": "code",
"execution_count": 15,
"metadata": {},
"outputs": [],
"source": [
"with open(\"md_openmm/min.pdb\", \"w\") as file_:\n",
" molecule.writeFile(\n",
" molecule.topology, molecule.positions,\n",
" file=file_\n",
" )"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" - __a)__ Open the structure in VMD and render an image the same way as in the last step in which the original and the minimised atoms are overlaid. Choose different colors for the two states. Can you see the difference?"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"![min](md_openmm/min.png)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Equilibration"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" - __a)__ We now want to equilibrate our system in the *NVT* ensemble for 100 ps. Add a state reporter to the simulation object that reports the total and potential energy and the temperature at every picosecond and write it to a file (e.g. `equilibration.log`)."
]
},
{
"cell_type": "code",
"execution_count": 56,
"metadata": {
"ExecuteTime": {
"end_time": "2020-07-01T08:39:57.642065Z",
"start_time": "2020-07-01T08:39:57.114149Z"
}
},
"outputs": [],
"source": [
"molecule = app.PDBFile(\"md_openmm/min.pdb\")\n",
"simulation.context.setPositions(molecule.positions)"
]
},
{
"cell_type": "code",
"execution_count": 17,
"metadata": {},
"outputs": [],
"source": [
"run_length = 50000 # 50000 * 2 fs = 100 ps\n",
"simulation.reporters.append(\n",
" app.StateDataReporter(\n",
" \"md_openmm/equilibration.log\", 500, step=True,\n",
" potentialEnergy=True, totalEnergy=True,\n",
" temperature=True, progress=True,\n",
" remainingTime=True, speed=True,\n",
" totalSteps=run_length,\n",
" separator='\\t')\n",
" )"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" - __b)__ Run the simulation. It should take about 10 minutes to complete on a processor with 4 cores. Reduce the number of steps, if it takes much longer on your hardware."
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {
"run_control": {
"marked": false
}
},
"outputs": [],
"source": [
"simulation.step(run_length)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" - __c)__ Plot the energies against the simulation time step."
]
},
{
"cell_type": "code",
"execution_count": 45,
"metadata": {
"ExecuteTime": {
"end_time": "2020-07-01T08:38:21.270030Z",
"start_time": "2020-07-01T08:38:21.259710Z"
}
},
"outputs": [],
"source": [
"pot_e = []\n",
"tot_e = []\n",
"temperature = []\n",
"\n",
"with open(\"md_openmm/equilibration.log\") as file_:\n",
" for line in file_:\n",
" if line.startswith(\"#\"):\n",
" continue\n",
" pot_e_, tot_e_, temperature_ = line.split()[2:5]\n",
" pot_e.append(float(pot_e_))\n",
" tot_e.append(float(tot_e_))\n",
" temperature.append(float(temperature_))\n",
" \n",
" # Alternative solution:\n",
" # for l, v in zip((pot_e, tot_e, temperature), line.split()[2:5]):\n",
" # l.append(float(v))\n",
"\n",
"t = range(1, 101)"
]
},
{
"cell_type": "code",
"execution_count": 47,
"metadata": {
"ExecuteTime": {
"end_time": "2020-07-01T08:38:30.851767Z",
"start_time": "2020-07-01T08:38:29.482167Z"
}
},
"outputs": [
{
"data": {
"text/plain": [
""
]
},
"execution_count": 47,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAABNgAAAL0CAYAAADeCIoLAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAuIwAALiMBeKU/dgAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOzdd5xU9b3/8feZvruzfdkCLB1EQUAFC4KCNXajmMQSr+16NcXca25izE1MzI1pPxMTjV5LbFgSNCqWRCNSVFAsgCBFelnK9j69nN8fA4ur22d2Zgdez8djH7Nnzjnf81ndHea851sM0zRNAQAAAAAAAOgTS6oLAAAAAAAAANIZARsAAAAAAAAQBwI2AAAAAAAAIA4EbAAAAAAAAEAcCNgAAAAAAACAOBCwAQAAAAAAAHEgYAMAAAAAAADiQMAGAAAAAAAAxIGADQAAAAAAAIgDARsAAAAAAAAQBwI2AAAAAAAAIA4EbAAAAAAAAEAcCNgAAAAAAACAOBCwAQAAAAAAAHEgYAMAAAAAAADiQMAGAAAAAAAAxIGADQAAAAAAAIgDARsAAAAAAAAQBwI2AAAAAAAAIA4EbAAAAAAAAEAcCNgAAAAAAACAOBCwAQAAAAAAAHEgYAMAAAAAAADiQMAGAAAAAAAAxIGADQAAAAAAAIgDARsAAAAAAAAQBwI2AAAAAAAAIA62VBcAoPcaGxv19ttvt22Xl5fL6XSmsCIAAAAAwKEsEAiooqKibfvUU09VXl5eCisaWAjYgDT09ttv6+KLL051GQAAAACAw9T8+fN10UUXpbqMAYMhogAAAAAAAEAcCNgAAAAAAACAODBEFEhD5eXl7bbnz5+vMWPGpKgaAAAAAMChbsuWLe2mKvrifenhjoANSENfXNBgzJgxmjBhQoqqAQAAAAAcblhorz2GiAIAAAAAAABxIGADAAAAAAAA4kDABgAAAAAAAMSBgA0AAAAAAACIAwEbAAAAAAAAEAcCNgAAAAAAACAOBGwAAAAAAABAHAjYAAAAAAAAgDgQsAEAAAAAAABxIGADAAAAAAAA4kDABgAAAAAAAMSBgA0AAAAAAACIAwEbAAAAAAAAEAcCNgAAAAAAACAOtlQXAAwEXq9Xy5Yt0+7du1VVVaW8vDwNGTJE06ZNU2lpaarLAwAAAAAAAxgBG1Ju1qxZevvttxPS1qmnnqolS5b0+Pjt27frjjvu0EsvvSSPx/Ol/VarVaeddppuv/12zZ49OyE1AgAAAACAQwtDRHFIycjI6PGxTzzxhCZNmqSnn366w3BNkiKRiBYsWKDTTz9dt956qyKRSKJKBQAAAAAAhwh6sOGQcumll/bouGeffVbXXXedTNNse85ms2natGkqLy9XTU2NVqxYoebmZkmSaZq65557FAgEdP/99/dL7QAAAAAAID0RsCHl/va3v8nv9/f6vPnz5+u//uu/2razsrL0jW98o9vzVq5cqWuvvbZduHbRRRfpvvvuU3l5edtzLS0t+u1vf6u77rqr7bkHHnhAkydP1o033tjregEglfx71svzWWw4vnviWXKWjE5xRQAAAMChg4ANKdfXRQTeeOONdttf//rX5Xa7uz3vhz/8oYLBYNv2nDlzNG/ePFks7UdMZ2dn65e//KUGDRqk//zP/2x7/ic/+Ykuv/xyZWdn96luAEgm344Vql/yF3m3Lm97rv7tv6jwtJuVf8p1MizMFgEAAADEi3fVSEsVFRVasGBBu+euv/76bs9bvHixFi5c2LZdVFSkBx988Evh2ufdcsstmjVrVtt2TU2N7rnnnt4XDQBJYpqmPJvfU8VfrtPuR29oF67FDoiqbuH92jfvB4oGOp6DEgAAAEDPEbAhLT3xxBOKRqNt2+PHj9f06dO7PW/u3Lnttm+44QYVFhZ2eY5hGPrhD3/YZTsAMBCY0ahaNyxWxUNXae/cb8u/c1WXx3vWL1LFQ1crWLczSRUCAAAAhyYCNqQd0zT1+OOPt3uuJ73XIpGIXn311XbPXXvttT265tlnn62ysrK27a1bt2rNmjU9OhcA+psZjahlzRva9cA3tO/ZWxXYs77H5wZrtqniwavk2fhuP1YIAAAAHNqYgw1pZ9GiRdq+fXvbtt1u19VXX93teR999JHq6uratsvKyjRu3LgeXdNiseiUU07RvHnz2p57/fXXNWnSpF5UDgCJEQ365N+zTv6KT+Xf/an8u1Yr4qnve3v+Vu195nv752W7nnnZAGAAC9ZXyLf1A5lmVFljT5Y9f0iqSwIAiIANaeixxx5rt33BBReouLi42/PWrl3bbvukk07q1XWnT5/eLmBbt25dr84HgL4wTVOh+l2xMK3iU/kr1ihQtVmKRhJ9IdUtfED+vZ+p9NJfyOLMSmz7AJAA4ZZaeTYtlWFzyD1+lizOzFSXlDSB6q2qX/KIWte+KZmmJKnGsCjn2ItUeNpNsuV0/34YANB/CNiQVhobG/Xiiy+2e64nw0Mlaf369kOmxowZ06trjx49usv2ACCRIp4GNX74vJo/fkHh5uq428sYfYJyjrlQ9YseVKi+otPjPBsWqeKh7Sq74h45iobHfV0ASAQzGlXDsrmqX/iAzEhIklSTkavii36i7AlnpLi6/hWo2hIL1tYtaAvW2phRNa94SS1rXlfeSVcof+Y1srpY6R4AUoGADWnlmWeekd/vb9seMmSIzj777B6du2XLlnbbw4YN69W1v3j85s2be3U+APREsL5Cje89o+aVL8sM+bs/oRtZR5yi/FOvV0Z5bEh71rgZqnzudnm3vNd5DTXbtfNPF8uaXSR7QbnsBUNlzx/yue+HypqVL8Mw4q4PALoTbq5W5Qs/lW/bh+2ej/qaVPm3H8g77VIN+sr3ZXFkpKjC/hGo3KT6xQ+rdf3Cbo81Q341vPOYmj56QQWzblDu8V+TxeZIQpUAgAMI2JBWvjg89Nprr5XVau3RuY2Nje22ezKstKvjW1paFI1GZYlzrqLq6mrV1NT06pwvhoUA0p+/4lM1LJur1vWLJDPa/QldMQy5jzpDBadeL2fZEe12WTNyNPib96pu4QNqeOexThqIibTUKtJS2+FqpIYjU86yI+Q+cpbcE86QPW9wfDUDiFs0HFS4qVLhpkqZ4YAsTnfsy5UlizNbFmemDEvP3jclSrilVp7Ny+Tb+qEi3ga5hkxQznFflT2/Z68ZrRsWq2r+LxT1NnZ6TPNHL8i/8xOVfu03cpb0boRCskV8LZLM2P8Hi23/o7XdBxaBfRtVt/hheTYs6nX7UV+Tal//vRrff1aFp39b2ZPOYV7NJDFNU4G96xVq2CuLI0PWrHxZswpkzcqXxe5KdXkAkoCADWnjk08+0cqVK9u2DcPo8SqgktTa2tpuOyOjd59yfvF40zTl8XiUnR1fN/wHHnhAd955Z1xtAEhPZjQqz8Z31LBsbochVq9ZbcqeeLYKTrlOjuJRnR5mWKwqOvO7cpaNV9VLP5MZ9PX6UmbQK//OVfLvXKXaN+6Rc+hEZU84Mxa29fDGGUDvmKapYM02hep2KdxUqVDjPoUbKxVqij1GWmu7bcNwZMrqcsvizJIlM0+Zo05Q7rRLZMselJgao1EF9qyTZ9NSeTa9q8DeDe32e7e8r/p3HlPW+FnKO+lyZYw4rsPesNGgT7X/ukdNHz7fo+sGq7eq4sGrVHTOrcqddtmA62HbvOpV1S64T5GWTj5UtdhkWGOBWzTgift64cZ9qnrhJ2pYNldFZ31PWWOnx90mOhfxtajy+dvl3bysw/2GI0PWzHxZ3QWyZubL5i6QY9BIZY45SY6SsQPu9zWVogGPAlVbFKzarGDNdkVDgdjrlTNz/4cGWbEvl1sWx/7nXG5ZM/NksTtTXT4OcwRsSBuPPvpou+3Zs2dr1KjObyC/6IsBm8vVu0+SOgrkWltb4w7YABwezEhYoca9CtXuVLB2u4I1O+TbsVKhup19btPiypZr6NFylU+Sq/xoucqP7tXcO9kTz5Rj0Ejte/bWLudl64nA7rUK7F6r2n/dI+eQCcqeeCBsY3U7IBE8m5ep9l9/UrAqvikqzKBX4aC3bdu/c5Uals1Vwax/V95JV/RpWGHE1yzvlvfl2fiuvFveU8TT0E0RUXk2LJJnwyI5SsYq78RvKHvSOW1DPAOVm1T53O0K1mzr3c8WDqjm1V/Lu+UDlVx8h6yZub3+WfpD/ZJHVLfwga4PioZlRsMyuz6q14KVm7R37reVOXa6ir5yq5zFo7s/qQ8ivmb5tq9QuKVakiEZRqx3nmFIhkWGYZH2f1nsLrmGT5Ytq6Bfakm2iK9Fe568WYE9nS+AZgZ9Cgd9Cjfubb/jX3+ULadEmeNmKGvcDGWOOn5ALtwRDfrUum6BvNtXKNJaK6u7UPbcUtlyS2XLK2v7vjfDtM1oRKGG3QpUblawcrMCVbHHUMPuPtdpODL39xzMl21/70FrZn7bc67yScwvi35FwHaY+853vqP777+/36/zs5/9TD//+c/7fH4gENCzzz7b7rmeLm7Qmd5+UsQnSwCk2HCnQNVmmeFgbChnNCrTjMYmnjYjMk1TMqMyIyGF6vcoWLtdoZodCtVXtE3M3SeGRY7i0XKVT1LG/kDNXjg87qE/zpIxKr/paVU+/+NOP3nvrcCedQrsWafaf/1RztIjZB80Uva8MtnzB8uWN1j2vDLZ8soYMgP0QDTkV+2//qimD+Z1f3AfmUGv6t78k5pXvKRB5/63ssbN7P6caESeTUvV/PGL8mxe1ueVjYNVm1X98v+q9s17lXvcV2XNylPdwgdir7F95NmwSLv2rFPpZXcpY8RxfW4nXqZpqn7xQ6pf/FBC23UUj1bBrH9XxNuo+sUPdR9oSvJufk+7tn6g3GlzVHjaTbJm5sVdR7BupzyfvSPPxnfk27mqd78D+6cyKDrnVtlzS+OuJVUi/lbtnfvtLsO17oSbq9T88Qtq/vgFGVa7MkZO3R+4nSxHYerCINM05a9YreaVr6h17Zs96llpycyTPbdU1uxBsfdC4aDMSDD22O4roGjQG9ffeYc17/8AIdywR4FOjskYcayKzvpPucqPTui1AYmADWnipZdeUn19fdt2fn6+Lrnkkl614Xa72237fL0bEtXR8V9ssy++9a1v6bLLLuvVOVu2bNHFF18c97UB9JwZjajurfvVsOwpKRpO2nWt2UXKO/Fy5U69tN96Y1gzcjT4qj+pcflf1fjeMwo3VSas7UDlRgUqN3Z8XXehbHllsrkLZUbCXb4RNyU5S8aqYPaNyhw5NWH1JVo06FNg30b5965XYM96Bet2yuLMUvbRX1HOsRfxYQ16xb/3M1X9/ccK1mxPyvVCdbu096lblDlupgad+/0Ob+5DTVVqXjFfzSteUri5KmHXjvqa1LD0iR4d6xxylCItdV1eP9xcpd2P3aiCU69X/sxrk74Agmmaql/4gOrf/kvC2nSUjFXBrH+X+6jT2z5cyZlyvhqWPaWGZXO7H+4fjajpg3lqWf1PFcz+D+Ud/zUZNnuPr29GI/JXrJHns7fV+tk7CtXG8Xtpmmpdt0CeTe+q4NQblHfyN9NuUYZowKO9c78t/+5PE9amGQnJu+V9ebe8r9p//j/ZckpkdRfI4sqWNSOn3aMlI1tWV7YsmflyFA2TLbcsIf/GhFtq1PzJP9S88mWFanf06tyot1EBb6O077O46+gvvh0rVfHw1XJPOEOFZ34npSEmDj0EbEgLXxweeuWVV/Z6iGd/BGxZWVm9aqMjxcXFvV5wAUByRfytqnzuRwnr4dUTjkGjlDfj6tiwqSTcdBgWq/KnX6W8Ey9XqGGPQvW7Y18Nu9t935f52joTaa1TpLWu00+Zv8i3/SPt2f6Rck/4uorO+l7KVwyMhvwKVG5SYM/6/YHahtiQtg4WqfBt/UCBvRtUfMHtKai0f/j3bpB3y3LJMGTLKZE9t0S23BJZswel3Y1yfzKjUYUaKqSoKVtuSY9+b81oRA1L56pu0QNSJHmB/gHeTe9q59b3lT/9KhWceoMMu0veze+p6aO/y7NpafwLsfSVYSh/5nUqPO0/FA14VTX/Tnk2LO78eDOq+iWPqOnD55V30hXKPf5rSRk2apqm6hbcp4Z3H09Ie87SI1Qw+0ZljZ/1pV7LFmeWCk+7SbnT5qh+ycNq+vjFbnuSRf0tqn39bjV9+JyKvnKrso445UvBjBmNKNy4T8HaHQrW7lRg7wZ5Ni/rcrGJvjBDftW99Wc1r3xZg877QY96Tw4E0YBXe576rvwVa/r1OuHmqh4H2RZnlhyDRslRPEqO4tFyFI+Ws2SMrNmDOgzeTNOUGfIrGmhVNOBVsGqLmle9EleP1HTSuu4ttW5Yotypl6hg9o2yuQtTXdKAZEbCinjqFW6uVri5Rs1rVqS6pAHNME0z0UP9kUYWLFigFSv6/49kxowZmjFjRp/O3blzp0aNGqVo9OCbuVWrVmnKlCm9aueCCy7Qa6+91rZ9//3361vf+laPz1+7dq2OPvpgV+KcnBw1NTX1qoZEWbdunSZOnNi2vXbtWk2YMCEltQD9xb97nRqX/1URT73sBeXKOnKWMkdM7dWn7YkQqt+tvU9/r9dzAfVVxsipyp9xtTLHnDzgVn4zTVMRT71C9RXybftYLesWKFi5KSW12POHquSSnyd9+Fe4pVat695Sy9o3YzdWvbwJKb7op8qd2rse2ANNsGa7at/8kzyfvd3pMdasAtlyS2TLiYVu9oKhyjriFDkKhyWx0i8Lt9Sqdf0iebcuV8TTIKu7QM6ScXKWjpOzbJxseYPj7gES8TbJv3uN/BWfylexRoHd6xQNxOaBNax2uYYfo6yx05U5ZrocJWO+dL1Q4z5VvfBT+XbE9/7MsNrjG5a+nzW7SIbFFnfPVnvRCElGn3s92XKKVTLnrnY9WE3TVNNHz6v29d/3aKiZ4chQ7tRLlTf9KtlzS/pUR3dM01Ttv/6oxmVz427LOfhIFcz691iw1sPfy2DtTtW99We1rnurx9fJGH2Cso8+W6H63QrW7lCodmdsWoMED9/riawjTtGgc38ge8HQpF0z1LBXvl2rZcsuUsawKd2+z4gGfdr71Hfk27Gyy+MGCovLLUfRSEmxXnfRgEfRoEfRgDd1YfkAYzgylT/jauVP/+aAnAMvEUINe+XbuUrRgEdmJBT7CgelSCg2iuDAc5GgIp5GhVuqFW6uVqS1vt3vyeYan8576ODIBO5D2yNgw4D385//vN0qm8cee2yfQsFbb71V99xzT9v2D3/4Q/32t7/t8fn/+Mc/dP7557dtH3fccfr44497XUciELDhUGZGQqpb/LAa3nk0Nq/Z51hcbmWNO0VZR81W1tiT+70Hk2/HSu376/cVSfAn9l9iWOSeeKbyT/6mXEPS6285WLdTrWvfUuu6txRI9pAQw1DeiZer8Izv9OvvQthTr9Z1C9W69s1Y6BHHWyfDatfQGx6Ta+jE7g8eYMIttapb/KCaV8zvU+8Gw2pX8UU/Uc4xF/ZDdZ0LN1erdf1Cta57KzZPVBf//yzOLDlKxsYCt9JxcpSOldWVI8ncP9ei9t9omJJpxnqARIKxYcG7P5V/15peLVxizR6kzDEnxQK30SfIs/l91bz2K0X9rd2ea8sfEptY/MAE4/vnNvz8ZOPRcFBRf6uigRZF/Z5YTxV/q6KBVoXq96jx/Wfbwr9EM2wOZYyYqqwjZihz3Aw5CsplmqZ82z5Q4/t/lWfTuz3+W8o66jSVXNT5ogWBys2qfO5HPf8gxGpTzuTzlH/y1V2uuNxbpmmq9vXfq/H9Zzo/yLBo0Pk/kmvIUTKjESkSlhmNyIyGYwsdRMKSKTlKx8ieP7TPga+/4lPVvPF7+Xet7uNPkzqGzaH8Gdcof+Y1/fra7t+9TvVv/0Wez5a0PWfJyFX2pK8oZ8oFcg456kv//aNBn/Y+fYt82zu/B7DllmrodY/I4nIr4mmI9f7xNOz/PrYdrNos367Vh0VPsb4yHJlylo6VLXuQokFvLBj0HwgHPbHXyQQHhFZ3oQpm/4fcR86W1V14SEzrEA0HVb/wATW893RCft8I2LpGwIYBzTRNjRw5Ujt3Hnyz2tueZwc88sgjuvHGG9u2L7nkEr3wwgs9Pv++++7TLbfc0rZ91VVX6amnnup1HYlAwIZDVahhjyqfv13+iu7nMzFsTmWOOVHuI09T1vhTEjJh8+c1r3xFVa/8b78Mz7LllclRNFKOouFyFI9S5rgZaT3J8wHBul2xIGrdAgX2bkjade2Fw1Ty1TuVMfzLPZtN01Sobpd8O1fKt3OVAnvWKxr0ypZdJFtuWax3VW5p2/e2vFJZM/MV9TWrdcMitX76przbP0roTZAtp0TlNz8rmzvxK+hFgz7JYk3oEM1owKuGZXN7NsdTdwxDpZf9WtlHn52Y4joRatzXFqqlTcBgGD0KnGw5xSq55BfKHH1C3JcMt9apbsF9al75ctxtSbFgIWtcLFDLHHV8l+FIsL5CTR/MU/OKlzsN+Qy7S4PO/YFyjvtqtze60aBPNW/8Xs0f9fy9nSRlHTlb+TOvUUb5pF6d90WmaarmH7/tejEKw6LSOb9U9qRz4rpWb2pqXfumat/8k8KN+5JyTWt2kSzO7FjwsX8BINOM7F8IKKpoyK+or7lnbWUVyOou3P93YbY9HvwzMWVxuZVRPllZR85SxvBjYiuXdsO3Y4Xq335U3i3vd3mcY9BIZU+5QDlTzpMtpzgWrj3zn/Jt+7DTcw6Eaz3pgRfxtci7dbk8m96Vd9MyRTz13Z6TatbMPGWNP1XRgEehpkqFGysVaa2Nr1HDkL1gmJylY/d/wBH7kMOWW9ZlT/6DQ1w9iviaDgaY3s+HmbGvcHN1r1dttzizZC8cJkfRCNkLh8tRNGz/43BZnPFPE5QMgcrNqvz7/8S9+vTnEbB1jYANA9qCBQt01llntW1nZGRo7969ysvr/Y30+++/r+nTp7dtl5WVae/evV2c0d43vvENzZt38E3Tb37zG9122229riMRCNiQaKZpKtywR9GAR5aMHFkz85I+v1Xz6tdV8+qv+tabwmKNvRnLK5M9pyS2bHxOcewxt1i27EEyrD0bWmpGI7G5c5Y+2e01LXaXZFgkwxJ7E3jge8MiGYasmXmyF42Qo2iEHING7H+TNizlc4clQ6h+t7zbPozN3da4T+HGvQo17lOkpaZ/LmgYypt+lQpn36Rg3U75dq6Sf8dK+XZ9okhrXe+asjljPUv6cTGLjJFTNeTf/k+GNTHT4fp2rFTd4ofk2/6RJMlRNFKu8qPlGjpRrqFHy1E8utfXMiNhNa14KbZKYS//G3bJYtPgK/+Q8LmWoiG/Wj55TU0rX1Zg99qEtj1QuCecoeILf5LwecT8u9eq5h+/69Nk7YbNIfeEM5U77VK5hk3pdY+PaMCr5k9eU+Pyv7UbPuocfKRK59wlx6CRvWqvZd1bqp7/C0X9Lb06zzX0aOWeeLmyJ5zR66kIzGhUNf/4jZo+fL7zgyxWlc75lbKPPqvzY/pJNORX43vPqP6dRxM6j+YBzqET5T7iVGWNP0WOkrFd/g6Y0YiaPnpBdW/9udf/j7pzIPzJOnK2Mkef0G6latM05d26XPVLHpF/56reNWxYlDn6BJnhQJfDQm05xRpy/V/kKCjvde1mNBqb527TUnk2vRvXqqQJZ1iUOfZk5R57UWy+vi/8fUTDwdhccY2VCjftU6ixMjZXn9Uuw+aQxeaQceDL+rnvbQ7ZsovlKB7V7++LTNOMLRzx5p8SMrWF1V0Y+8rKly2rQNasfFkz82OPB7bdBbIXlKekF5wZjarx/WdUt+C+hEwV8HkEbF0jYMOA9sVQK55eY+FwWCUlJe1WI924caPGjRvX7bnRaFTl5eXtArnVq1dr0qT4Pu3sKwI2xMuMhBWo3CjfzlXy7Yj17PnixMWG3SVrZt7+Nw15+9845MnqLlTm6BPlHHxkQt40RAMeVb/2G7V88lr3B/eVYcjqLpQ9r0z2gmGyF5bLXjBMjsLY99aMnLZaKv/+P13OLSXFekwNvvJPcgwa0X81H6KioYDCzZUKN+5TqCEWukUDrTKsn38T7pRhs8cerXaZ4aDqlzzSs4meDcuAmFPG4nTLXjRMgT3rOz0m7+SrNegr/xXXdQL7Nqr2rT/Lu2lpl8cZdpecg4+Ua+jRcpUfLWfpuFjPkpBf0XBAZiggM+xXNBT7PupvUeMH83q9glxPGTanhvzb/QmZQy/ibVLTh8/tn7OxIQHVDTwWZ5YGnXebsqec3283a2Y0qpbVr6n2zXt7FKjai0Yqd9qlyplyfkICP9M05d/9qYLVW+UoHC7X8GP6/LOGW2rVsGyumj56QWbQ26tzre4i5U67VLnTLpUte1CXx0Z8zfJu/UAtq//Zbpjhl1hsKv3ar5U94Yxe1ZJo4ZYa1b11v5pXvRLfUHe7S5mjT4iFWeNmypZd1PtaPPWqW/BnNa+cH1ctndboyFDWmOnKOuo0WWxO1b/7WJevx/GyZg/S0OsfSdiKlOHWOvl3rVHEU6eIr0VRf7Mi/hZFfS2K+lsU8TUffOynaSzshcOVc+xFbT34DgVmNKqWNf9U3Vv3J3TF9M7YC8qVN/1K5U67LGnz6oaaKlX14s+67GkZDwK2rhGwYcCqr6/X4MGDFQgcXF9u8eLFmjVrVp/bvOaaa/Tkkwd7pfzoRz/Sr3/9627Pe+ONN3TOOQe7848ePVpbtmzpcx3xImDDAYHqrbFP4gyrLI4MGY4MWewZMuyu/dsuWewZkmEosGd9W5jm2/VJr286vihj1PEqPv9Hve5d8Hn+3etU+fztCtVXxFVLvCyZeXIUlCvia1KobleXx2aMOl5lX/9dUlaiw0ERf4tqX/99woayJZLhyJRr8Hg5Bx8l15Cj5Bx8lOwF5ZJhqOrFO7oMj0u/9ts+9WgJ1leofuH/qeXTN/rl5rQ3YkOnvD2aN+zzLM4sDbn2YbmGHNWn64Ya96rxvWfUtOKlfumVM1C4hk1R6Zxfyp4/JCnXi/hbVb/kETUuf/ZLQ+QNm0Puo86I9VaLIwBLloivWU0fPKfG5c/2Pny12r+oir8AACAASURBVJQ94UzlnXi5XOWxRa7MSEj+3Wvl3fK+vFuWy79nXfeBvtWmsq//Tu4jZ/fxp0g8/94Nqn397h5P0m/LKZF9/5QGWWOmK2PUtHa9w+KqZfdaVb/2m4HVY6uXrO4iDb3ukZR96BbxtyhYvU3B6q1tj4HqLYq09H7opiUzT+7xpyrn2IvlGjZ5wP+N91U0FFDTh/NU//ajPR6yHI+MkVNV8tU7Zc8f3ONzTNOUZ8MStXz6hsyQX/ai4XINPlLOIUfJnl/eYWDXsuZ1Vb/664T3Dj3AsNq11ePSV+5+t+057kPbI2DDgPXFOc9Gjx6tzZs3x/VCv2jRIp1++ult20VFRfrss89UWNj1ssynnXaaFi8+uAz8nXfeqTvuuKPPdcSLgA3h5mpV/+O38qxflNpCrDblz7hGBade36s322Y0qoZlc1X31v3dDsOzuGJzuUQDnnirjVvO1EtVfP5tPR5uisTzbFqqqvm/6L+hpj1kLxqh7IlnyT3h9Njwy07m/YmG/Nr98DUKVG7scL/hyFD5jU/JWTK6R9cNt9Sqfskjavr4xX4dwtoT9vyhKjzzu3JPPFOGYSga8CjcVKVQc3VsuFBTlYKVm9S6fmGnbVgz8zT0+kd7NdF8oHKTGpbOjYWLfZwbz140UlljTlSoqUrBqs39EvJbswcpo3ySXMMmyTV0kgybXd4ty+XZvCw2z2R3//8sNhWedpPyZ17To3mlEi1Ys0MNy56Ub/sK2XIGKWv8LOUcc0HC57tMhmjIr+aVr6hh2VyFG/b0+nznkAmyZQ+Sb/tHvfq3yLDaVXr53XIfcUqvr9nfTNNU6/qFanj7UQX2fSbDkfG5uaZGyFE0PDbFQeGwfl9V0YxG1bzqZdW9eW//LyqUYFZ3oYZe+3BCF8tIlIi3ScGabQpUbVGkuTo2XNOZJYvTLcOZKYvTLYszc/9zWbI4smRxuQ/ZUK0jEV+zGt55XI3Ln+33VXMNR6YGnfP9Hs0p6d36gWoX3Ndp8GxxuuU88MHe4CPlKB6l+nceU+un/+q2DovLHZvfzmaPDdm12vd/2fYP47XLcGTKljNItuxBsSlX9j9aMvO0fv167kO7QMCGAWvKlClavfrgxMR33XWXfvzjH8fd7umnn65Fiw6GEnPmzNG8efNk6aTb7r333qvvfe97bdtFRUXatm2bsrOz466lrwjYDl9mNKKmD5+PzV0yAAKnA2z5Q1R8/m3dzqkU8TWrdd1bal4xv0dz/WSMnKbSOb+UJTNPvm0fqnXDYnk2LEn+RMCGRYPO+b5yT7z8sHrjOVBFfM2qef1utax6NanXtReUyz3xLGUffVa38wx9Xqhhj3b935WK+po6brdwmMpvelpWV+f/rkT8LWpYOleN7z0tM+TvU/2JYsnIVcHsG5U7bU63iymYpqnaN36vxvc6X1XRllOsoTc83uUn+2YkJO+2D9X4/rPybn6vT3U7isfIPfEMuSecIWdx+0AzGvAoULVZgcpNClRuUrByswJVm3vcM86w2uUcPF6u8liY5ho2qcuFSyL+Vvm2fSTPlmXybn5f4cb2c8LaC4erdM5dcg3l3/ZEMiNhta57S/XvPp6QeZi6YtgcKrv87oTPNdgfzGhk//yhqf33LeJrVt3CB2Jz2SVxqL+rfJLyT71eZsCj5lWvybt1eY+ub83K15DrHvnS6wnST6ipUg3vPK7WdW/1+3vMzLHTVXLxzzocduvfs051C/4c+x3sB1njT1XxRXfEtcgS96FdI2DDgLRixQpNnTq1bdtqtWrXrl0aPLjn3Wo7s3LlSp100kkKBg9+SnHxxRfr3nvvVXn5wUlJW1pa9Lvf/U533XWXPv9n8tBDD7VbjTQVeGE7PAX2bVTVK78c0JN3u486XUXn/kD23JK256KhgDyb3lXL6tfl3fRuzyZbtVhVeNrNHfbcMKMR+SvWxFYIXL/4SzemiWZxulX69d8oa+zJ/Xod9F7rxndU/fL/9mwYjGGRs3ScMoYfI3vRcIVbahVu3KdwU+X+XleVHa4Ya8sri/VUm3hWXPMOera8r71zv93pcM6s8bNUdvnv24Z8hFtq5a9YI3/FGvl2rVZg7/p+/3S9O4bNqbyTrlD+zGtlzej5h0ymaap6/p1dDu+1F5Rr6A2PtpvzKhrwyrPlPXk2LJFn4zt9GvLiLD1C7glnxHoa9nI4uxmNxm60zGhshU/DIsmI/Q58bluGZHFk9rmXmWmaCtXulHfbBwo318hRPEruCWckdCVYtHdgwvOGpU/2yzxFhs2psiv+oKyx07s/GF8Sbq2Xd+sHsb+//a+5hrT/e6Pt0YyG5dv+sbxb3u/TBw8Zo45Xwak3KGPk1Hav7eHmajWv/qdaVr2qYM22Ds+1ZubFwrWSMb3/ATGgRbyNCtbuVKh2p4J1Bx53KVS3S2Y40H0DPWBxZcfm1Zx8rgzDULBmh+oWPqDWdQsS0v4XGY4MDTrnv3vUe6473Id2jYANA9K3v/1tPfDAA23b559/vl59NXE9FZ555hldddVV7Z6z2Ww6/vjjVV5ertraWn300Udqbm4/Jv/mm29uV1eq8MJ2eIkGfapb9KAa33+mz8OhkslwZKrwtJvkLBuvljX/VOu6t3o1N5M9f6hKL/tV25w3XTFNU8HKTfLvWadwU1XbsLRQU5XCzZVxz8tkzx+qsqv+yKfTA1jE26Saf/4/taz+R7vnDZsjtoLm8GOVMXyKXOWTZXW5O23nQJgSbqpUqCkWttkLhyVsMQ9Jqn/7UdW99edO92dPOV+KRuWvWKNQw+4+XcNVfrRyjrlIwZpt8ld8qsC+z+JfQcwwlD3lAhWefnOXvbK6YkYjqnzuR2pd91anxziKx6js8rvl27lKng2L5d36QZ9uZgxHhnKnXqLc4y9L2ITjOHQFqraq6YO/qfmT1xLSO9SeP1TFX/2ZMkdO7f5gJEQ06JN36/JYL/fP3um0t/ABmeNmquDU65UxbHKXx5mmqcDe9Wpe9apa1rzR1q69aKTKvvE7wrXDjBmNKtxcqWDtToWbqxXxNOz/qo89ehsUaa1XxNvY49eSrCNPkzUrL/YBVD+9x3eVT1LJnF/2aXXbjnAf2jUCNgw4fr9fZWVlamw8OAfDSy+9pIsvvjih13n88cd1yy23qLW1+xt/wzB0yy236Pe//72s1uTPg/JFvLAdPjyb3lX1q79WuHFfv7Rv2JxyDZ2ojBHHKmPEsXIUj1HU36qItzH2RsHbuP9NQ2PsE72qLQrs3dAvtUhS9uTzNOj8H3UZhPSUaZqK+lv3h26VseCkfreCdRUK1e1SqL6iy5t31/BjNPjy38ualR93Leh//r2fybt5mQy7S64hE+QcctSA6wFkRqPa99fvd73aYB85Bo1S4ZnfUdb4We0CwWg4GAuhd38qf8Va+Xd/2ul8Y20ruNqdsthcMpyZyhg2RblTL5Gz7Ii4a4yGg9r3zH/Ju6VvQzy7Y80qUN5JVyj3+MvaVgYGeiria1bzypfV9MFzvQ64bbmlyhxzkjJHn6isI2cNuNeew4kZCbeF9K0bFh9cKdIw5D7qdOWfcr1cg8f3vt1wSIGqzZJpylk2jrlY0aVg3U5VvXSn/DtXpa4Ii02Fs29U/sxrZVhtCWuW+9CuEbBhwPli77KSkhLt3r1bNlviXhgO2LZtm+644w7Nnz9fHs+X57OyWCw67bTT9OMf/1izZw+c1Z94YTv0hVtqVfPP/6fWtW/26HjDkSlFI9329rA43XINmxwL1IYf2+sQwjRNtXzymmr/dU/vV2TrguHIVPEFtytnyvkJa7M7ZjSqcEtNW9gWqtulYH2FDMOizNEnKufYC3kDjYSL+FtU8eA3FarbmZD2bLmlKjztZmVPOa/HQxQjvhZFfI2y2JyxQO1AsNbJXKSJFA36tOfJb8m/65OEtWkvKFf+jH9T9pTzZbE7E9YuDk9mNCLPpqVqWv63TudBMhwZyhw5TZljTlTm6JNkLxqe8vnL8GWmaSpYvU3hpn1ylh3Rbgg60N/MaESN7z+rurf+HPcUD67hxyjqbVKwdkeP5ge0F41Q6ZxfyjUk8feH3Id2jYANkOTxeLR06VLt3r1b1dXVysvL0+DBg3X88cerrKws1eV9CS9shzb/7nXa+/R3exRgGXaXCk+/WXknXiHDapMZjcgM+RUN+WUGfYoGfbHtcEDWzDw5Bo1MyGp0EW+Tat+6T80fv9jpnFI9lTHqeBVf+D9yFA6Luy4gHQSqt6rioW/GNYTZmpmn/FNviC00kGahUsTXoj2P/XunK6v2lHPIBOXPvEbuI2enZJVNHPqC1dvU9NHf5d2xQhZnljKGH6vMMScqo3yyDBsfwADoXrB6mypfvKPTFUG74hp+jIrOukUZw6ZIin1IFdi3Uf696xXYs16BvRsUrN1+8L24xarcaXNUdNb3ZHFkJPLHaMN9aNcI2IA0xAvbocu3Y4X2PHWLzKC322Mzx81U8QU/kj0v/sU/+spXsUY1r/yq1zfKttxSZU86R9mTz2UOExyWWtYuUOW8H/b6PEtmnvKO/5ryTv5mQoZSp0q4tV67H71eododvTrP4sxS5tiTlXv8ZcoYcRy9hgAAA54ZCavh3SdUt+ShDhdU+iJH6TgVnfldZY49udt/56IBrwKVmxQNtMhRMrbPc6X2FPehXUv8mDsAQJ94Ni/Tvr/+d7cTo1rdRRp03g/lnnBGym8uM8onqfymp9X4wTzVLXygy2DQ4sqWe+JZypl8jlzDjknKcDRgoMqeeKb8u69W47K5XR5nyxusjGGT5SqfJNewyXKWjE3oXCqpYnMXaMg1/6fdj1x7cI6kTljdhcoaP0vuI2cpY9TxzG8FAEgrhtWmglk3KOuImap84acKVm3u8Dh7/lAVnvEtuSee3eP3yRZnpjKGT0lkuYhD+r9DA4BDQOu6hdr3/I+6/lTLMJQ77TIVnvkdWV3ZySuuG4bVpvzpVyp74pmq+efd7ZYYN2wOZY2bqezJ5ylz3MncGAOfU3TWLZKkpuV/lRkJybA55Bx8pFzlk5VRPkmuYZMO6TmD7LmlGnLtg9rz2H8o3FzVfl9BubKOnC33UbPlGno0Q0ABAGnPWXaEym96WvWLH1bD0ielaOx9v9VdqIJZNyr3uK8y/DzNMUQUSEN0zT20NH/ymqpe+nmXy3M7Ssaq+ML/6XZJ+YEgWL1N3u0fy5ZdpIyR02TNGDhhIDAQRfwtivpaZM0uOixD6Ii3UY0fPKdQfYUcRSOUNf5UOYpHp7yHLgAA/SXibZJ3y3uy5ZTIOWRC2synyn1o1+jBBgAp1Pjh86p59VddHuOeeJZKL/1l2nyi5SgeJUfxqFSXAaQNqyt7QPVKTTZrZp4KZ9+Y6jIAAEgaa2ausiedk+oykGAEbACQIg1Ln1Ttv/7Y5TE5x16k4ot+yvAoAAAAABjACNgAIMlM01T9ogdVv+ThLo/LO/FyFZ3z3ywGAAAAAAADHAEbACRINOBRy9o35a9Yo2jAK5mmJFOmGZVMSWZUkqmov0W+HSu7bCv/lOtVeMa3mYMIAAAAANIAARsAxCka8Kjxg3lqWPaUot7GuNsrPPO7KjjlugRUBgAAAABIBgI2AOijaMCjxuV/iwVrvqaEtDnovNuUd+I3EtIWAAAAACA5CNgAoJf6I1iTYVHJxT9TzrEXJqY9AAAAAEDSELABQA9F/K1qOjAUNFHBmiRZbCq97FfKnnhm4toEAAAAACQNARsAdMMMh9S4/K+qf+exxAZrkixOt0ov+5WyjpiZ0HYBAAAAAMlDwAYAnTBNU57Plqj2jXsUqq/o8XmGI1Puo06Txe6SDItkGJJhyJCxf1uSDNnyyuQ+6nTZc0v67WcAAAAAAPQ/AjYA6ECgcrNqXr9bvm0f9vgcw5GpvBMvV/7JV8mamdeP1QEAAAAABhICNgD4nLCnXvWLHlTTRy9IZrRH5xiOTOWddIXyp19JsAYAAAAAhyECNgDQ/nnWPpyn+sUPKepv7dE5FmeWck+8nGANAAAAAA5zaR+wvfDCC/rBD34gwzC0devWVJcDIA15tryvmtd+q1Ddzh4dbzgylX/SFcqbfpWsmbn9XB0AAAAAYKBL+4CttbVVO3bskGEYqS4FQBpqeO8Z1b5+d88ONgzlHHuRCs/4jmzuwv4tDAAAAACQNtI+YAOAvmpaMb/H4VrGiGNVdM4P5Bo8vp+rAgAAAACkGwI2AIel1nULVf3y/3Z7nC1vsIrO/k+5J5xBT1kAAAAAQIdSFrDt2rUrIe3U1tYmpB0Ahw/v1uWqfP72LlcJNRwZKjjleuVNv0oWuzOJ1QEAAAAA0k3KArYRI0bQGwRA0vkrPtXeZ2+VGQl1ekz2MReo6IzvyJZTnMTKAAAAAADpKqVDRE3TTOXlARxmAlVbteep78oM+jo9puic7yt/+lVJrAoAAAAAkO5SFrAd6L1WWlqqcePG9bmdyspKbdy4MVFlAThEhRr2aM+TNyvqa+r0mIJTbyBcAwAAAAD0WsoCtjFjxmjLli0aP368Fi5c2Od2nnzySV177bUJrAzAoSbcUqs9T9ysSEtNp8fkHv81FZz+rSRWBQAAAAA4VFhSdeHjjjtOpmlq1apVqSoBwGEg4mvRnie/pVB9RafHuI/+igaddxvzQgIAAAAA+iRlAdvUqVMlSU1NTdq6dWuqygBwCIsGfdr79HcVrNrc6TGZ42ao9NJfyLCk7OUQAAAAAJDmUh6wSdLHH3+cqjIAHKJM01Tlcz+Sf9fqTo9xDT9GZV//nQyrPYmVAQAAAAAONSmbg+3YY4/V5MmTJUk1NZ3Pi9SdGTNm6PHHH09UWQAOEc0r58uz8Z1O9ztKx2nwlX+SxZGRxKoAAAAAAIeilAVsbrc7IfOvjR49WqNHj05ARQAOFeGWGtW+8YdO99sLyjXk3x6QNSM7iVUBAAAAAA5VTDoE4JBT89pvFfW3drjPllOsIdc8KJu7MMlVAQAAAAAOVQRsAA4presXqXX9wg73GTaHBl99v+z5g5NcFQAAAADgUEbABuCQEfG1qPq1X3e6v2DWjXKWjEliRQAAAACAwwEBG4BDRu2bf1KkpbbDfY7SccqfcXWSKwIAAAAAHA5StshBT23fvl3vvvuuJOmKK66QzTbgSwaQAt7tH6v54xc63mlYVHLxHTKs9uQWBQAAAAA4LAz4tOq9997TNddcI4vFoosvvlg5OTmpLgnAABMN+VX98v92uj9v+pVyDZmQxIoAAAAAAIeTtBkiappmqksAMEDVL35YobpdHe6z5Q9R4Wk3J7kiAAAAAMDhJG0CNgDoSGDfRjUsm9vp/pKLfiKLIyOJFQEAAAAADjcEbADSlhkJq2r+nVI00uH+7GMuVOboE5NcFQAAAADgcEPABiBtNb7/rAJ7N3S4z5pVoEFfuTXJFQEAAAAADkcEbADSUrC+QnWL/q/T/YPOu03WzNwkVgQAAAAAOFwRsAFIO6ZpqvrlX8oM+TvcnzX+VLknnpnkqgAAAAAAhysCNgBpp2XN6/Jt+7DDfRZnlgadf7sMw0hyVQAAAACAwxUBG4C0YoZDqlv4QKf7C8/6nuy5JUmsCAAAAABwuCNgA5BWmla8qHDDng73uYYfo9yplya5IgAAAADA4Y6ADUDaiAZ9ql/ySMc7LTaVXPgTGRZe1gAAAAAAycWdKIC00bj8r4q01nW4L+e4i+QoHpXkigAAAAAAIGADkCYivmY1vPtEh/sMm1OFs25MbkEAAAAAAOxHwAYgLTQsfVJRf0uH+3JP+LpsOcVJrggAAAAAgBgCNgADXrilVo3vP9vhPovTrYJTrk1yRQAAAAAAHETABmDAq3/7LzJD/g735Z38TVkz85JcEQAAAAAABxGwARjQQg171PTxCx3us2blK3/6lUmuCAAAAACA9mypLqA7brdbw4YNk2EYsljIA4HDTd3ih6RIuMN9+adcL4szK8kVAQAAAADQ3oBPrC666CLt2LFD27dvl9vtTnU5AJIoUL1VLZ/8o8N9ttxS5U6bk+SKAAAAAAD4sgEfsAE4fNW99YBkRjvcVzD7RlnsziRXBAAAAADAlxGwARiQ/LvXyrNhUYf77EUjlDPlgiRXBAAAAABAxwjYAAxItW/9udN9had/S4Z1wE8hCQAAAAA4TBCwARhwvNs+lG/rBx3ucw4+Uu6jTk9yRQAAAAAAdI6ADcCAYpqmahfc1+n+wjO+LYMVhQEAAAAAAwh3qQAGFM9nSxTYvbbDfRkjjlXmmOlJrggAAAAAgK4NuEmMrrvuupRc1zAMPfrooym5NoCD6t/u/O+w8IzvyjCMJFYDAAAAAED3BlzA9sQTTyT9Bto0TQI2YAAINe5TYM+6DvdljpupjOFTklwRAAAAAADdY4gogAHDs+ndTvcVnfHtJFYCAAAAAEDPDbgebFKsRxmAw49nY8cBm6N0nJxlRyS5GgAAAAAAembABWzbt29PdQkAUiAa9Mm37aMO92WNm5nkagAAAAAA6LkBF7ANHz481SUASAHf9o9lhgMd7ss6YkaSqwEAAAAAoOeYgw3AgNDZ/GuWzDy5hh6d5GoAAAAAAOg5AjYAKWeaZqfzr2WNPVmGxZrkigAAAAAA6DkCNgApF6zaonBTZYf7ssYxPBQAAAAAMLARsAFIuc6Gh8piVebY6cktBgAAAACAXhpwixwkSiQSUX19vQzDUH5+vqxWhpgBA5Vn09IOn88YNlnWjJwkVwMAAAAAQO8cMgHbzp079dRTT2np0qX66KOP1NjY2G5/Xl6epk2bphkzZuiqq67SiBEjUlMogHYi3ib5d63ucF/muJlJrgYAAAAAgN5L+4CtsrJS3/ve9/Tiiy8qGo1Kik2Y/kUNDQ1asGCBFixYoDvvvFOXXHKJ/vjHP6qsrCzZJQP4HO+W9yQz2uG+LAI2AAAAAEAaSOs52BYuXKjJkyfr73//uyKRSFuwZhhGh19SLHyLRCL6+9//rkmTJmnBggWp/BGAw15nq4fa8srkKB6V5GoAAAAAAOi9tA3YPvzwQ1144YWqqamRaZrtAjTTNFVYWKhRo0Zp1KhRKiwsbHteUtuxdXV1uvjii/XBBx+k7OcADmdmNCLP5vc63Jd1xCltf6sAAAAAAAxkaRmweb1effWrX5XP52u7ATcMQ3PmzNErr7yi2tpaVVdXa/Pmzdq8ebOqq6tVV1enV199VZdddpksFkvbOT6fT5deeqm8Xm8qfyTgsOSv+FRRX1OH+7LGzUhyNQAAAAAA9E1aBmz33HOP9u3bJ8MwZJqmxowZo+XLl+u5557T+eefr4KCgi+dk5+fr/POO0/z5s3T8uXLNWbMmLZ9+/bt0z333JPMHwGAJM+mjoeHGnaXMkZOTXI1AAAAAAD0TVoGbI888khbuDZixAi9++67mjq15zfjxx13nN555x2NGDGirZ2HH364HysG0JHO5l/LHHW8LHZXkqsBAAAAAKBv0i5g27Bhg3bt2tU279r999+vkpKSXrdTUlKiP//5z23zsu3evVvr169PdLkAOhFq3Kdg1eYO92UyPBQAAAAAkEbSLmBbvXp12/dDhgzROeec0+e2zjnnHA0dOrRte82aNXHVBqDnPJuWdrov64iZSawEAAAAAID4pF3AVlNTIym2QMHkyZPjbu/zbRxoG0D/83YSsDlKxsqeW5rkagAAAAAA6Lu0C9g8Hk/b9zk5OXG3l52d3WHbAPpPNOSXd9uHHe6j9xoAAAAAIN2kXcBWWFjY9v2+ffvibq+ysrLt+45WHwWQeL7tH8sM+TvclzWOgA0AAAAAkF7SLmArLY0NHTNNU8uXL4+r15nH49Hy5cvbtsvKyuKuD0D3Ols91JKRK1f50UmuBgAAAACA+KRdwDZjxgxZLBYZhqFAIKC77767z2394Q9/kN8f60VjsVh08sknJ6pMAJ0wTVOeTR0HbFljT5ZhsSa5IgAAAAAA4pN2AVt+fr5OPPFESbEb9V//+td65ZVXet3Oa6+9prvuukuGYcgwDJ1wwgkMEQWSIFi9TeHGjod3Zx0xI8nVAAAAAAAQv7QL2CTpxz/+sUzTlGEYCgaDmjNnjm677Ta1tLR0e25ra6tuv/12XXrppQqFQjJNU5J0++2393fZAKROe6/JsChzzPTkFgMAAAAAQALYUl1AX5x77rk655xz9Prrr8swDIXDYd199926//77df7552v69OkaN26ccnNzZRiGmpqatGnTJr333nt67bXX5PV62wI6wzB09tln67zzzkv1jwUcFjwb3+nwedewybJm5ia5GgAAAAAA4peWAZskzZs3TzNnztTq1atlGIZM05TX69Xzzz+v559/vtPzDvRYO3DOpEmT9NxzzyWrbOCwFvE1y1+xpsN9rB4KAAAAAEhXaTlEVJLcbreWLFmiOXPmtOuNJsVCtI6+JLU75tJLL9WSJUvkdrtT9nMAhxPv5vekaKTDfVlHELABAAAAANJT2gZskpSbm6vnnntO8+fP14wZM9oFaR05sH/mzJmaP3++nn/+eeXl5SWxYuDw5tm0tMPnbbmlchSPTnI1AAAAAAAkRtoOEf28Cy+8UBdeeKF27typpUuX6uOPP1Z1dbUaGhpkmqYKCgpUXFysqVOnasaMGRo+fHiqSwYOO2Y0Is/mZR3uyzpiZlvvUgAAAAAA0s0hEbAdMHz4cA0fPlxXXnllqksB8AWBPesV9TZ2uI/hoQAAAACAdJbWQ0QBpA/vjhUdPm/YXcoYOS3J1QAAAAAAkDgEbACSwrdjZYfPZwybIovdleRqAAAAAABIHAI2AP3OjEbk37Wqw30ZI45NcjUAAAAAACQWARuAfhes2qKov7XDfa7hxyS5GgAAAAAAEivtFznweDx67rnntHDhQn3yySeqqqpSc3OzwuFwr9oxDKPX5wDoGd/OjnuvGVa7XEMnJrkaAAAAAAASdGdjdAAAIABJREFUK60DtnvvvVc//elP1doa6xljmmaKKwLQkc7mX3MOmcD8awAAAACAtJeWAZtpmrr22mv11FNPtYVqhmHIMIy2/T1lGAbBHNCPTNPstAcb868BAAAAAA4FaRmw3XvvvZo7d66kgwGZaZrKyMjQ6NGjlZubK5stLX804JATqt+lSGtth/symH8NAAAAAHAISLsUKhwO6xe/+EW7YO3cc8/VbbfdphkzZrT1YgMwMHQ2PFSGRa5hk5NbDAAAAAAA/SDtArZ33nlHDQ0NbUNCb7rpJt1///2pLgtAJzobHuosHSerKzvJ1QAAAAAAkHiWVBfQWxs3bpQUm9cpOztbd999d4orAtCVznqwMf8aAAAAAOBQkXYBW0NDg6TY3GvTp09XRkZGiisC0JlQU5XCDf+fvTsPb6rM2wd+nyxt0y3d932Blr3sRUBWAQVBURhxe1FcxlFndGacd94ZGcZhXh39zTivOoqKozMq4gKigIIOi8iibGUpa3e6N13SJW2TJjm/P0IDaZM2bdMkbe/PdeVqz5JzntA25Nzn+zxPqdVtXhx/jYiIiIiIiAaJARew+fld61IWHBzswpYQUXdai2yMvwZAEc8KNiIiIiIiIhocBlzAlpaWZv6+trbWhS0hou7YGn9NHpIAmW+Qk1tDRERERERE1D8GXMA2ffp0eHt7QxRFHDt2zNXNIaIucPw1IiIiIiIiGgoGXMCmUChw//33AwBqamrw+eefu7hFRGSNoVkNXVWe1W0Kjr9GREREREREg8iAC9gA4LnnnkNkZCQA4Be/+AUqKytd3CIi6shW91CAFWxEREREREQ0uAzIgC04OBg7duxAQEAAiouLMX36dBw5csTVzSKi69gK2GTKCMgDopzcGiIiIiIiIqL+I3N1A3orIyMDP/zwA+68806cOXMG06dPx/Tp07Fw4UKkp6cjICAAEknP8sOZM2f2U2uJhh6b469x9lAiIiIiIiIaZAZswAYAqamp+Otf/4oVK1agrq4OBw8exMGDB3t1LEEQoNfrHdxCoqHJqG2Gtvyi1W2KBI6/RkRERERERIPLgA3YGhoasHr1amzbtg2AKSADAFEUXdksIgLQUnwGMBqsbmMFGxEREREREQ02AzJg02g0mD17Nk6dOgVRFBmuEbmZ1iLr3UOl3gGQhyY6uTVERERERERE/WtABmy//e1vkZWVBUEQIAgCRFGEr68vbrjhBqSmpkKpVEImG5AvjWhQsDX+mlf8eHMgTkRERERERDRYDLgUSq1W4+233zYHazKZDH/+85/xxBNPwMvLy9XNIxryjHodWkuyrW7j+GtEREREREQ0GA24gG3//v3QarXm6rV//OMfeOihh1zdLCK6Slt6HqJea3Ubx18jIiIiIiKiwUji6gb0VF5eHgDTeGtRUVEM14jcTIuN8dcED294RgxzcmuIiIiIiIiI+t+AC9iMRiMA06yhEydOdHFriKgjW+OvKeLGQZAOuKJZIiIiIiIiom4NuKvd6Oho8/fe3t4ubAn1t8bGRpw4cQI5OTlQq9Voa2uDUqlEZGQkJk6ciLi4OIedq7m5GYcOHUJJSQkqKysREBCA6OhoTJo0CREREQ47z2AnGg1ovXLa6jZFPMdfIyIiIiIiosFpwAVsKSkp5u8rKipc2BLqL7t27cIrr7yC3bt3mysWrUlNTcXDDz+Mn/3sZ1AoFL06V0FBAdauXYvPP/8cGo2m03apVIo5c+bgt7/9LWbPnt2rcwwl2oocGLVNVrcpEjj+GhEREREREQ1OA66L6OTJk5GQkABRFPHjjz+itbXV1U0iB2lqasKKFSuwaNEifP31112GawCQk5ODX//61xgzZgyOHTvW4/O99957GDNmDD744AOr4RoAGAwGfPvtt5g7dy6efvppGAyGHp9nKGm1Nf6aVA7P6JFObg0RERERERGRcwy4gA0AHnnkEQBAS0sLXnvtNRe3hhyhpaUFixYtwqeffmqxXiKRYNy4cVi+fDnuuusuzJ49G76+vhb75ObmYv78+Th50nq4Y82mTZvwwAMPoKnpWrWVTCZDZmYmVqxYgdmzZ8Pf39+8TRRFvPzyy3jyySd7+QqHBlvjr3nGjIJE7unk1hARERERERE5x4AM2J5++mlMnjwZoihi7dq1+O6771zdJOqj9evX4+DBgxbrli9fjry8PGRlZeGzzz7Dpk2bsHfvXlRVVeGll16Cp+e1wKa+vh733Xcf9Hp9t+c6efIkVq9eDVEUzeuWLl2K/Px8HD58GB9//DH27t2LkpIS/O53v7N47uuvv4633nqrj692cBJFES1FWVa3KeLZPZSIiIiIiIgGrwEZsMnlcnz11VeYMmUKWltbsWDBAqxfvx4NDQ2ubhr1QnNzM/7+979brFu1ahU+++wzJCQkdNpfoVDgV7/6Vadqt3PnzmHLli3dnu+ZZ56BTqczL99xxx3YunUrYmNjLfbz8/PD+vXrO7Xt97//PRobG7s9z1DTVl0Eg6bW6jaOv0ZERERERESD2YCb5AAAnnvuOQDAnDlzcPnyZdTV1eEPf/gD/vKXvyAzMxPp6ekIDAyERNKz/HDt2rX90Vzqxp49e9Dc3Gxe9vDwwMsvv9zt85YsWYLFixdjx44d5nXbt2/HypUrbT5n37592LNnj3k5JCQEGzZs6PJ35cknn8S2bduwf/9+AIBKpcLLL7/M35cOWmyMvwZBAq/YMc5tDBEREREREZETDciAbd26dRAEwbwsCAJEUYRGo8GePXssApSeYGDiGvn5+RbLU6ZMQVhYmF3PXbp0qUXAlpOT0+X+//73vy2W16xZg+Dg4C6fIwgCnnnmGXPA1n4c/r5Ysjn+WuRwSL18rW4jIiIiIiIiGgwGZBdRawRBsAjdeuL6sbjI+TrO4BkTE2P3czt266yrq7O5r8FgwPbt2y3WrV692q7zLFiwAJGRkeblvLw8nDlzxu52DgUcf42IiIiIiIiGqgEbsImi6LAHuVZERITFcmtrq93P7bhvUFCQzX2PHTuGmpoa83JkZCSGDRtm13kkEglmzpxpse7rr7+2u52DXVt9BfTqMqvbOP4aERERERERDXYDsovovn37XN0EcqAZM2ZYLJ88aWMsLytOnDhhsTxp0iSb+2ZnZ1ssZ2Zm2n0eAJg2bRo+/vhj8/K5c+d69PzBrKXQevUaAHjFj3NiS4iIiIiIiIicb0AGbDfeeKOrm0AOlJqaivnz5+Pbb78FABQVFWHHjh1YvHhxl8/TaDR45513LNbdd999Nvc/f/68xXJKSkqP2pmcnNzl8Yay1uLTVtd7hCZC5mO7qpCIiIiIiIhoMBiwXURpcHn99dcRGBhoXn7ggQdw6tQpm/s3NjZixYoVKCu71i1x9erVXVaw5ebmWizHxcX1qI0d9+9uQoWhRFeVZ3W9Vxyr14iIiIiIiGjwG5AVbDT4pKSkYM+ePVi+fDkKCgqgUqkwdepU3HPPPbj11luRkpICuVyOiooKHDhwABs2bEBJSYn5+YsXL8Ybb7zR5TnUarXFsr0zldrav7GxEUajERIJc2qdqsDqes/wnlUJEhEREREREQ1EDNjIbWRkZODs2bN4/fXX8d577+H8+fN45513OnUDvV5MTAx+97vf4ZFHHul2FtmmpiaLZYVC0aP2ddxfFEVoNBr4+fn16DgdVVVVQaVS9eg5HavxXMnQ0gBDU43VbfLQJCe3hoiIiIiIiMj5GLCRWzEYDAAALy+vbvcdMWIE/vSnP+HWW2/tNlwDOgds9pzjetYCuaampj4HbK+//jr++Mc/9ukYrqRT5dvc5hGa6MSWEBEREREREbkG+7YNcY8//jgEQej3x7p167ptyxdffIHk5GQ888wzds0kev78eSxfvhwpKSn48ssve/za7Qnl+rL/UGGre6jE0wcy/551wyUiIiIiIiIaiBiwkVt4//33cfvtt6O6utq8bvjw4Xj11Vdx7tw5NDQ0QKvVoqSkBNu2bcNtt91m3q+oqAhLly7F+vXruzyHr6+vxXJLS0uP2mht/47HHIpsBWzy0ESGkkRERERERDQksIsoudyFCxfw8MMPw2g0mtc99NBDeO211+Dh4WGxb3R0NKKjo7F06VJs374dK1asQGtrKwDg2WefRWpqKlauXGn1PP0RsPn4+PToGNY89thjuPPOO3v0nNzcXCxbtqzP53YEWwEbu4cSERERERHRUMGAbYhbunQpYmJi+v0806dPt7ntueeeM4dkADB79mxs2LCh29k5lyxZgtdeew1r1qwxr3vqqaewbNkyeHp6dtpfqVRaLPd0YoGqqiqLZX9/f4fMIBoWFtbjGU3dic2ALYQBGxEREREREQ0NDNiGuPnz52P+/PkuO79Op+s0ftratWvtDq5Wr16NP/3pTygqKgIAlJeXY/fu3bj11ls77Zuammqx3P4ce3Xcv+PxhiKjrgV6dZnVbaxgIyIiIiIioqGCY7CRS+Xk5KC5udm87Onp2WW1W0cSiQRz5syxWPfjjz9a3Tc9Pd1iOTc3twctBfLzLWfL7Hi8oUhXXQSIotVtHmEM2IiIiIiIiGhoYMBGLqVWqy2Wg4ODIZP1rLAyIiLCYvn6iRKuN2rUKIvlI0eO9Og8hw4d6vJ4Q5Gt7qGCVA55QLSTW0NERERERETkGgzYyKUCAgIsljUaTY+P0dTUZLFsa2bPSZMmISgoyLxcXl6Oy5cv23UOo9GI77//3mLdokWLetjSwadNlW91vTw4HoKUPdCJiIiIiIhoaGDARi4VGRlpsVxfX4+CAutVUbacOHHCYrljRVs7mUyGJUuWWKx799137TrHN998g7Kya2ONJScnY8yYMT1q52Bkc4KDsCQnt4SIiIiIiIjIdRiwkUsFBQV1Cqreeustu59/5swZ/PDDDxbrZsyYYXP/++67z2J548aNqKmp6fY8L774YpfHGapsBmyc4ICIiIiIiIiGEAZs5HIrV660WP7b3/6GvXv3dvs8tVqNe+65B0aj0bwuNjYWkydPtvmcOXPmWEyKUF1djUcffdTiGB298sor2Ldvn3k5JCQETz31VLftG+xEgx662itWtzFgIyIiIiIioqGEARu53M9//nOLrqI6nQ6LFi3Cc889h7q6uk77G41GbNu2DePHj8fZs2cttv3v//4vJJKuf61feukleHh4mJc/++wzLF++HMXFxRb7NTY24tlnn8UvfvELi/V//vOf4efnZ/frG6zaaksAg97qNgZsRERERERENJRwFHJyOR8fH2zZsgXz5s1Dc3MzAFPI9oc//AHr169HRkYG4uLi4OHhAZVKhRMnTqC2trbTcR566CHcc8893Z5v/Pjx+Oc//2mx77Zt27Bjxw5MnjwZsbGxqK6uxrFjx9DQ0GDx3J/+9Kd4+OGH+/iKBwedjQkOIEggD453bmOIiIiIiIiIXIgBG7mFzMxM7NmzB/feey9yc3PN69va2nD06FEcPXrU5nNlMhn++7//G3/84x/tPt/dd98NnU6HJ5980jwLqV6vx+HDh63uLwgCnnzySfz1r3+1+xyDna3x1+SB0ZDIPZ3cGiIiIiIiIiLXYRdRchtTp07F6dOn8eqrr9o1Q6dSqcQjjzyCU6dO4U9/+lO3XUM7Wr16NU6fPo27774bPj4+VveRSCSYN28e9uzZg7///e+QSqU9OsdgxgkOiIiIiIiIiEzcsoLtwIEDTj2fVCqFr68v/Pz8EBkZCYVC4dTz0zXe3t54/PHH8fjjj0OlUuH48eMoLi6GWq2GXq+Hv78/goODMWbMGKSnp/c4VOsoKSkJH3zwATQaDQ4ePIiSkhJUVVUhICAAUVFRmDx5ssX4cHSNrS6iDNiIiIiIiIhoqHHLgG3WrFkQBMEl55ZIJEhKSsINN9yAu+++G3PnznVJOwgIDQ3FokWLnHIuHx8fLFiwwCnnGgxEoxG66kKr2zxCk5zbGCIiIiIiIiIXc+suoqIoOv1hMBiQk5ODf/3rX7jppptwww034OLFi67+pyByK/qGSoi6FqvbWMFGREREREREQ41bB2yCILjsAZgCviNHjmDGjBm4dOmSi/81iNyHrfHXAEAemuC8hhARERERERG5AbfsIhoXF+f0LqKtra3QaDTmGSUBU8AniiJqamrwk5/8BFlZWU5tE5G7shWwSf1CIfXyc3JriIiIiIiIiFzLLQO2wsJCl51brVbj1KlT+PTTT/HPf/4TOp0OoijizJkz+OKLL7B06VKXtY3IXXCCAyIiIiIiIqJr3LqLqCsEBARg1qxZ+Mc//oH9+/fDy8vLXE23Y8cOF7eOyD3YqmBjwEZERERERERDEQO2LkyZMgWPPvooRFEEABw5csTFLSJyD7YDNs4gSkREREREREMPA7ZuLFmyBIBpwoPKykoXt4bI9fSaWhib1Va3sYKNiIiIiIiIhqIBGbBVV1c7/JhvvPGG1fXR0dEATBMeqNXWQwWioaStixlEPcJYwUZERERERERDz4AM2JYuXQqdTuew423atAlPPPGE1W0eHh4ATBVsRqPRYeckGqhsdQ+VKPwh9QlycmuIiIiIiIiIXG9ABmxHjhzB6tWrHXKsr776Cv/1X/9lHmeto5iYGBQUFKCgoAD5+dZnTiQaSrqa4KB9QhAiIiIiIiKioUTm6gb01ubNm5GSkoI//vGPvT7G999/jzvvvBN6vd5mMCCVShEfH9/rcxANNpxBlIiIiIiIiMjSgKxga7d+/Xps2rSpV889deoUlixZgtbWVge3imhw01VZr+RkwEZERERERERD1YAM2MaNGwfANC7agw8+iEOHDvXo+Tk5OViwYAEaGhogiiIEQcDzzz/fH00lGlSMWg30DdZn02XARkREREREREPVgAzYtm/fjsjISAiCAK1Wi9tuuw0FBbZnNrxeSUkJ5s+fD5VKBUEQIAgCnnnmGTzzzDP93GqigU+nKrS5zSOUM4gSERERERHR0DQgA7bo6Gh8+eWXUCgUEAQB1dXVuOWWW1BfX9/l82pqanDTTTfhypUr5jHX1qxZw+o1IjvpVNa7hwpyL8iUkU5uDREREREREZF7GJABGwCMHz8eH3zwgbkK7dKlS1i+fDkMBoPV/ZuamrBw4UJcvHjRHK7dcccdePPNN53ZbKIBzeYEByEJECQD9u2EiIiIiIiIqE8G9BXxsmXL8MILL0AURQDAvn378Mgjj3TaT6vV4tZbb8WJEyfM4dpNN92EDz/80KntJRroOIMoERERERERUWcDOmADgF/96ld48MEHIYoiRFHEu+++i7/85S/m7QaDAStXrsT+/fvN4drUqVOxdetWyGQyVzWbaEBiwEZERERERETU2YAP2ADgjTfewOzZswGYZhb93e9+h88//xwA8MADD+DLL7+EIAgQRRGjRo3CV199BYVC4comEw04Rr0ObbXFVrfJGbARERERERHREDYoAjaZTIatW7di2LBhEAQBRqMR9957L1auXIn333/fHK4lJyfjm2++gVKpdHWTiQactporgGi0uo0VbERERERERDSUDYqADQCUSiV27tyJoKAgCIKA5uZmfPbZZwBMVW1RUVH49ttvER4e7uKWEg1MtrqHQiKFR1CccxtDRERERERE5EYGTcAGAMnJydi6dSvkcrm5ag0AAgMDsWvXLiQkJLi2gUQDmE6Vb3W9PCgWgkzu5NYQERERERERuY9BFbABwIwZM7Bx40ZzuObr64uvvvoKo0aNcnHLiAY22xMcJDm5JURERERERETuxe2m0XzuuecccpykpCQUFBRgypQp2L17N3bv3t3tc9auXeuQcxMNRpxBlIiIiIiIiMg6twvY1q1bB0EQHHIsURSxd+9e7N271679GbARWScaDWirLrK6jQEbERERERERDXWDrovo9ewN6tq7kxKRdXp1OUS91uo2BmxEREREREQ01LldBRvAwIvI3WhtTHAAMGAjIiIiIiIicruAbd++fa5uAhF10GZj/DWZMgISD4WTW0NERERERETkXtwuYLvxxhtd3QQi6oAziBIRERERERHZNqjHYCMix9DZ6CLK7qFEREREREREDNiIqBuiKNquYAtjBRsRERERERERAzYi6pKhqRrG1iar21jBRkRERERERMSAjYi6Yat6DWDARkRERERERAQwYCOibtgK2KQ+gZB6Bzi5NURERERERETuhwEbEXVJV8UJDoiIiIiIiIi6woCNiLpkq4JNzoCNiIiIiIiICAADNiLqhq660Op6j1DOIEpEREREREQEMGAjoi4Y9ToYGlVWt3kExzm5NURERERERETuiQEbEdmkb6i0uU0eEOXElhARERERERG5LwZsRGSTXl1uc5ssINKJLSEiIiIiIiJyXwzYiMgmvbrC6nqJdwAkHgont4aIiIiIiIjIPTFgIyKb2mxUsMlZvUZERERERERkxoCNiGzS11sP2GTKCCe3hIiIiIiIiMh9MWAjIptYwUZERERERETUPQZsRGSTrUkOOMEBERERERER0TUM2IjIKtFohL7e+iQH7CJKREREREREdA0DNiKyyqCphWhos7qNXUSJiIiIiIiIrpG5ugE9UVlZicLCQvj5+SE5ORmenp6ubhLRoGWreyjALqJERERERERE13N5BZtOp8OGDRuwcuVKLFq0CE899RQuXbpksc+lS5cwZ84cREVFYdq0aRg9ejRCQkLwyCOPoLa21kUtJxrc2mx0DxXkXpB6Bzq5NURERERERETuy6UVbBUVFZg3bx4uXLhgXvfNN99gw4YN2LZtGxYsWIDS0lLMnj0blZWVEEXRvJ9Go8HGjRvx3Xff4cCBAwgLC3PFSyAatGxOcKCMgCAITm4NERERERERkftyaQXbHXfcgfPnz0MURYuHVqvFqlWroFKp8Oijj6KiogJeXl6YP38+7rrrLsycORNyuRyiKCInJwdr1qxx5csgGpTa1GVW18s5wQERERERERGRBZcFbNu3b8fhw4chCAKWLFmCw4cPIzs7G7/5zW8AAGq1Gi+++CJ27dqFBQsWoLCwELt378aHH36I/fv349KlS5g4cSJEUcTOnTtx7NgxV70UokHJ5gyiHH+NiIiIiIiIyILLAraPP/4YADB69Gh8/vnnmDp1KkaMGIHnn38ed911F0RRxMsvv4yYmBhs2bIFoaGhFs9PSEjAzp07ERhoGgtq8+bNTn8NRINZm40uopxBlIiIiIiIiMiSywK248ePQxAErFmzBhKJZTMefvhhAIAoinj44Yfh7e1t9RihoaG4//77IYoifvzxx35vM9FQYrOCjV1EiYiIiIiIiCy4LGArKzON7zR27NhO29LS0szfZ2ZmdnmcG2+8EQCQm5vrwNYRDW1GrQbGlgar29hFlIiIiIiIiMiSywI2rVYLAFAqlZ22XT8jaHBwcJfHiY6OBgDU19c7sHVEQ5ut7qEAu4gSERERERERdeSygK197LS6urpO2wRBMH8vlUq7PI5MJrNrPyKyn77eRsAmSCDzD7O+jYiIiIiIiGiIclnA1l6ZVlFhfZwne7UHdB0nQSCi3rNVwSbzC4UglTu5NURERERERETuTeaqEyclJeHixYsoKCiwul2lUgG4Vulmy/nz5wEA4eHhjm0g0RCmtxWwBXCCAyIiIiIiIqKOXFbBNmHCBIiiiBMnTljdHhwcjODg4E4zjHa0f/9+CIKAcePG9UcziYakNrWtGUQ5/hoRERERERFRRy4N2ABTQNZbtbW12LlzJ4Brs4kSUd/ZGoONExwQERERERERdeayLqKLFi1CcXFxn46xa9cuTJkyBQAwa9YsB7SKiABAb6uCjV1EiYiIiIiIiDpxWcAmk8kQHR3dp2OsWrUKq1atclCLiAgAREMb9I0qq9tYwUZERERERETUmcu6iBKRe9I3VAGi0eo2jsFGRERERERE1BkDNiKy0GZjBlGAFWxERERERERE1jBgIyILehsBm0ThD4mnj5NbQ0REREREROT+GLARkQV9vY0JDpSc4ICIiIiIiIjIGpdNctAXKpUK586dw/nz55GTk4P6+no0NjaiqakJAODr6wtfX18EBAQgNTUV6enpGDVqFEJDQ13cciL3Z6uLKLuHEhEREREREVk3YAK2wsJCbNq0CVu3bkVWVlavjpGRkYHbbrsNq1atQmJiooNbSDQ42Ooiygo2IiIiIiIiIuvcvotoZWUlHnvsMQwfPhzPPvssTp48CVEULR62dNzv5MmTWLt2LdLS0vCzn/0MlZWVTnwlRANDm40uovKAKCe3hIiIiIiIiGhgcOuA7fDhwxg9ejTefPNNtLW1QRRFCIJgfrTrGKR1DN6uf44oimhra8OGDRswevRoHD582BUvjcgtiaJoeww2dhElIiIiIiIissptu4h+8803WLZsGVpbWwHAHI6JogiZTIbRo0cjOTkZcXFxCAsLg0KhgEKhgCAIaGlpQUtLC6qqqlBUVIT8/HycOXMGer3e4ljV1dWYO3cuvvzyS8yfP9+VL5fILRg0dRDbWq1uk7OLKBEREREREZFVbhmw1dTU4P7770dra6s5DJPL5Vi2bBlWr16NGTNmwNvbu0fHbG5uxvfff493330X27Ztg06ngyAI0Gq1uO+++5CdnY3g4OB+ekVEA4O+3vr4awAr2IiIiIiIiIhsccsuomvXrkVlZaU5XMvMzMTZs2exefNmLFiwoMfhGgB4e3tjwYIF2Lx5M86cOYPMzExzN9KqqiqsW7fOwa+CaOCxNYOoIPOA1CfIya0hIiIiIiIiGhjcLmDT6/X4+OOPzWOsTZs2DXv27EFqaqrDzjFs2DDs3bsXmZmZAEzjTm3evNnchZRoqLI5/pp/OASJ271dEBEREREREbkFt7tiPnz4MGpra80TGmzcuBFeXl4OP4+npyc2btxoDvJqa2tx6NAhh5+HaCDR26hgY/dQIiIiIiIiItvcLmBgG77BAAAgAElEQVQrKioCYJqIICMjA2lpaf12rvT0dIwfP77TuYmGKltdROUM2IiIiIiIiIhscruArbz82gV+cnJyv58vKSnJ/H1FhfXucURDhc0uokoGbERERERERES2uF3A5uvra/6+vr6+38/X2Nho/t7Hx6ffz0fkzljBRkRERERERNRzbhewRUVFATBNPJCVldWvEw+0tbXh5MmTnc5NNBQZdS0wNqutbpMFRDi5NUREREREREQDh9sFbOPGjQNgGoNNpVLh7bff7rdzvfXWW6iqqjIvXz8eG9FQY2uCA4CTHBARERERERF1xe0CtoSEBGRmZgIwVbE99dRT2LFjh8PPs2PHDjz99NMQBAGCIGDq1KmIj493+HmIBoq2+i4CNv9wJ7aEiIiIiIiIaGBxu4ANAJ566imIoghBEKDT6bB06VKsWbMGly9f7vOxL1++jAcffBBLly5FW1sbRFEEAPzyl7/s87GJBjK92voEB1K/EEhkHk5uDREREREREdHAIXN1A6y54447sHLlSnz88ccQBAGiKOLdd9/Fu+++i8zMTMyaNQuZmZlISUlBXFwcFAqF1eO0tLTgypUryM3NxZEjR7B//34cOXIEAMwBHgCsXLkSt99+u9NeH5E74gQHRERERERERL3jlgEbAGzcuBENDQ34+uuvzUGYKIo4cuSIOSRr5+3tDYVCAS8vLwiCgJaWFrS0tKC5ubnTcduDtfbgbtGiRdi4caNTXhORO7M1BptMyYCNiIiIiIiIqCtuG7D5+Phg+/bt+M1vfoNXXnkFbW1tFkHb9TQaDTQaTbfHvD5Yk8vl+PnPf44XXngBEolb9pQlcipbY7Cxgo1ocNC1GVDXpINEEBDs7wmJRHB1k4iIiAaN8ppmfLQnD9X1rQj080RylB+So/yRHO0PX4Xc1c0jIidw24ANACQSCV566SU89thjWLduHbZs2WJRldYeuHWnPZATRREKhQJ33HEH1q1bh8TExH5pN9FAZLOCjQEbkVvTthlQU9+K6vpWVNdrUV3fipqGa8s1Da2oqdeiXqMzP8dTLkVsmA/iwnwRG276Ghfmi7hwX/h58yKAiIioJw6drcDv3zkBbZvB6vbwQAWS2gO3KH8kR/shIcIPMikLPYgGE7cO2NolJibiX//6F958803s3r0bu3btQnZ2Ni5cuIDa2tounxsYGIj09HSMGjUKCxcuxIIFC2yO2UY0VIkGPfSNKqvbZMoIJ7eGiK5X16jFmfxaVNe3QqVuD86uPRo0bT0+prbNgNzSBuSWNnTaFuDrgbgwXwyLVSI9PgAjEgIRG+rDijciIiIrdh8twfr3s2Awijb3qaxrQWVdC46cqzKv8/aUYvroCMydEI0p6aHwkEud0Vy3oFK34Gx+HYL9PTEmOcjuwhkidzcgArZ2Xl5eWLp0KZYuXWpeV1dXB7VajaamJjQ1NQEAfH194evrC6VSiaCgIFc1l2jA0DeqAKP1O27sIkrkGrWNWryyJRt7TpR1+aHd0dRNOqibanEm/9oNLF+FDOnxAUiPD8SIq6FbiNLLaW0iIvdVVq3BkXNVKK3WwCgCRqMIw9WH8erDYBRhFEUoPGUYnRiIuROi4OUxoC5DiKzacqAAf/vkLMRe/DfdrDXgm+Ol+OZ4KXy8ZJg5NgJzx0djUloo5LLBWdlW26jFW9svYsfhIrR/tEmLU+J/7hmHlGilw8+nNxiRV9qA7II6qNStiAhS4KZJMfD24vsP9Q9B7DigGRG5vXPnzmHUqFHm5ezsbIwcObLXx2spPImSdx60ui3pfw5AqvDr9bGJqGdEUcQ3x0rx98+yLbp1upvQAC8E+nmahmEQAaMoQgQgiqbXcHU1An09MH5YCG4YFY60uABWwhH1M3WTFrmlDTAYRYxKDIJPP11I1jVqsXHHRXxx6NqFsr0CfD1w56xELJ+ZCH8fD7ufd6WyCV/9UIzsglpIpRLMGR+FW6bGDohudq06PU7n1aK2QQt/bzniw30RGeID6SB5T6yqa8GuoyXQtLZheGwAZoyJGLQhEWD6f+69XTl4e8dFhx/bz1uOWeMiMWd8FCYMC3Hp77coiiivaUZWTg2ycmtwoVCNNoMRk4aHYNmMBKTG2BeKtemN+HR/Pt79+jI0rfpO22VSAQ/cPBz3zE/p0+tVN2mRXVCH7Pw6nC2oxYUiNVp1lkUEoQFe+P29GZiUFtrr8wxljr4OHWwYsBENQI5+Y2s4/RUqP/tdp/UST18k//77Xh+XaCiqb9Lhi0NFOHG5Gv7ecmSODMPsjCgoPLu/yK2sa8GLH5226EIymAT5eWLaqHDcMCocE9NC++3Cn8gZmlv1KKvRoLKuFTKpgLgwX4QHKpwaItc1anHpSj0uFqtx6YoaF6/Uo7Kuxbxd4SHFoqmxuOPGRCREOOZmWXcXyj2h8JBi6fR4/GROMsICrQ/homnVY+/JUuw8UmxRWdsuOsQbDy1Ow7wJ0W4X4JdVa3A4uwqHz1Xi5OVq6PRGi+0eMgliw3yQEOGH+AhfJET4ISHCF7FhvjAaRdP4mg1Xx9as7zC2ZoMWXh5SpMUqkTEsBBkpwQh2QWWx0Sji0+8K8Ma28xavL9jfE8umJ2DZ9HiXtKs/iaKIV7eew+a9+f1+Lj+FHDFhPogIUiAiyLvTV0ePmyqKIkpUGlOgllODU7k1Fu8pHY1JCsLtMxMwa1yk1S6uoiji4NlKvLr1HEpU3U9KmBYXgGfvy0BipH3vV/VNOhw4U45TuTXILqhDcVX352i3YnYSfnprOjw9hk7XXEdgwNY1BmxEA5Cj39hqv3sHNf95rdN6j/BUxD/+Sa+PSzSUqJu02LwnH599l49mreXdUm9PKeZOiMYtU2MxOqnzWCNGo4jPDxbijW3nOz13sJLLJMhIDcYNo8KROSIMSl/TzKZSCSCVSCCRCJAI9k9oRORooihCpW5FabUGpdXNKKtuRlm1BqU1zShVaaBu6lxh6uUhRVy4LxLCfRF/NSxJiPBDTKhPn6t5GpvbcL6oDucL63DxSj0uXVGjSt1q9/OnpIfizllJmDoirFdBlCiKOHCmAv/4/LxdF8o9IZMKWDg5FqvmJSMhwg+iKOJUbi12HrmCvVllnSpQrEmN9scjt6Yjc2SYy9439AYjTufV4nB2JQ5nV6Kossmp548P98X41GCnBW4qdQvWv5+FYxerbe4jkwqYOz4ad85OxIj4wH5tjzPoDUb8ZdNp7PyhuMv9UqL9UV7T3OcQujs+XjJEhXhjcloobpuRgKgQnx4fo01vxP5T5fj+TAVO5Vajul7b42ME+nlgybR4LJsej4ggbwBAXlkDXtmS3eXvhzVymQQPLR6Ou+amWK3w1LYZcOhsJXYfK8GRc5XQG3ofZyRE+OIP94/H8LiAXh9jsNC2Gcz/vwDA1BFhVoNOBmxdY8BGNAA5+o2t8ov1aDi+pdN672EzEH3vK70+Lg0uKnULPvg2FyUqDdLjA7D0hniEBnDSmNpGLT7ak4et3xWgxY6LwLgwH9w8NQ6LpsQgNECBK5VNeP7DUzid1/WkPfZQeEgREuCFYH8vBCs9EeLvhRDlte+DlV5o0elRXKnBlaom06OyCcVVGrsuYF1BIgBSqQRymQRxYT6YlBaKKelhGJ0U1OPAoqmlDXmlDWjTGxEV4o2IIG+3q3gh12tqacOXh4qw9UAhymqau3+CHaQSAdEh3ogL9706e68v4q7O5Bvk79kpENIbjCgob8S5gjpkF9bhXEGdw8Ka2DAf3HFjIm6eGmd3FWlOST3+b8s5nLzcswvlnhIE00VdcZWm1yHe2OQg/HRpOsYkBzu4ddaJooijF1XYfugKfrxQ1e+BSk/EXQ3cZo2LxKS0UIcGj3tOlOLFzWfQ2Gz/RDsjEwNx56xEzB4X5TbdR3NK6rFx5yWcyqmBt5cMqTHKq2OOBiA9LgBK32tdmHVtBqx77yT2nyrv8piP3pqOe29KAWCqTM8ra0ReWQPyyxqQX9aIworGPoVCtkgEYObYSNw1NxmjEgO7/Xmr1C3YdrAIXxwsQm1jz0M1W224YXQEgvw8sf1wz7uPX29kYiB+f28G4sNNFZ1n8mux62gJ9p0sQ2NLzyd4skUqEfDQ4jTcPd96oDdYtWj1OJtfi1O5tcjKqcb5IjXaOlTZjksJxrIZ8Zg19lqFIgO2rjFgIxqAHP3GVvrvn6E553Cn9crJKxC25Le9Pi4NHpeuqPHkq0csPkh7yCRYNiMB985P6Ze75FqdAUfOV6JarUVilB/GpQS71QefmvpWbNqTh8+/L+xVOCURTB9csgvqOnUbsibIzxMxoT4IUXohJMALIUpPhCoVV783PXrb5dJoFFFd34orVU0oqmzC5eJ6XChSI7+soU8fjvuTt6cU44eFYEp6GKakhyImzNdiu1ZnwOUS0+swPepwpUPXEW9PKRIj/ZEU5YekKH8kR/ohKdofQX6eznwpDtGg0eGdry5hX1Y5DEYj4sJ8kRztj6RI02tLivKDv7f9Y10NRVV1LfhkXz62HSpCsxNDEm8vmTlsC/DzRM7V39v+Dr29vWRYnBmHjJRgCAIgkQgQBAFSQYAgASSCAAHAf06WYfvhol4N4u5K00eH45El6UiO9jevE0XThAt6gxF6g+mrRBDg5y3vcfgkiiKOX6rG2zsuIrugztHNd7jxw0Lw9J2jkBTl3/3OXWhqacNfPz6L3cdKen2MEKUnFk6ORVpcAFKi/REd6vyx6LQ6A/759SVs+k9elxMJRQV7Iy3OFLgdvVjVZTWWIAC/XDEat89M7PLcbXojTlyuxp4TpThwusKhYVG7kQmBWDknCbPGRVqMadZeHbrluwJ8d7rcqZMo9YaHXIIFk2Jw/FI1yh10w8OWMUlBePb+DER3qAJsbG5DTkk9LpfU49KVeuSVNUDXZsDEtFDcd1Oqze7t7qa5VY9Tuaax9E7l1ODiFbXdP/8AXw8szozD0unxUFcWMmDrwpAO2GpqavDrX/8agKkLyjvvvOPiFhHZx9EBW9Ert0OnKui0Pnj+kwiaubrXx6XBoaC8EY+9fMjmgPuecinuuDEBd89PQYBv34MJURTx3elyvPxpNlTXdX8KUXrhponRWDglpkczTWla2pCVW4OC8kbIZaYPaoF9CFCq61vx4be5+PxgIXRt3QdjfSWVCLh/QSruW5BqdXyT/tSi1eNycT3OF6lxvqgOFwrVDqvocbSoEG9MSQ+DKIo4X6RG3tVB3nsjwNcDyVH+SI8PwPjUEIxJDnLrGccKKxrxzIaj3Vb8hAZ4mQO35Cg/jEkORkxoz7sTdcVgFFFV14LiqiYUqzQoqdKgtNrUrikjwnDL1Fi3mz0yp6QeH+3Jw7fHS93+YtOdhQcqkBqjtOzmLTG9h0kEAXqDiO/PlDu1G7wgmLrQtYdptqqGAv08MGNMJOZkRGK8HYPKn7xsCtYcUXnsTFKJgJWzk7D65uG9uiGTlVON5/6V1eWYXL3hKZciKcoPKdH+SI7yN3+9vnrMkY5fUuHFj844tKuzVCLg2fszcNPEmB49r01vxNELVdhzsgwHzlQ4PNwPD1TgzlmJmD8xBoeyK7DluwLklTU69By94e0lw81TYrHraDGaWvr/hoaHTGLXzUxvTykeuHk4dHojLhebQrWyatufe7w8pLj3plSsmpcMTyd/RrOHKIq4cEWNLw4W4dvjpX2+cSMIQHJAE97/893mdQzYLA3pgK2oqAiJiYnmO1YGg3t2jyHqyJEBmyiKyFt/A0Rd5w9LEXc+D78xC3vdTnIvoij2+A59abUGP/3bIVTXdz/Oj7enFHfOTsJdc5N7XSlTXtOMv31yFoeyK7vcLzXaHwsmx+CmSTEI6VA9p2szILugDscvVeP4JRUuFFneoVN4SPHbe8Zh3oToHrVNFEV8tCcPb+246JRgDQDS4wPwP3ePs6jAcLW6Ri0uFKlRrNLAYDBCEExjpeFqtYsgmCpfIACtWgOOX1LhZE5Np24HA4lUIiAtPgATUkOQMSwYY5KCupy0QttmgErdiqq6FlTXt8LTQ4rJaaF2TXTRUz+cr8Lafx7v9QVK+9g9k9NDMWFYqN0DZtc36ZBf3oiiikYUqzRXu/Q1obS6ucufdYCvB34yJxm3z0yAr8Kxg3P3hCiKOHZRhU178nD0gspl7RgMvD2luH/hMKyYndTtBWZDsw6fHyjEx/vyrY5hZ6+kKD8UVTT1SyAa4OuBmWOth22nc2vw9s5Lfeoq6+0pdflYmyFKLzxx+0jMmxBl1+cCXZsBb22/iI/25jm1kjEiSIHpoyMwOyMSY5L7XsXeoNHhtc/PY8eRKw5qoYmHXII/r5mEG0aF9+k42jYDfjhfhT0nynDobIVdw044U3KUH8amBCOvtKHX4bIgAIsz4/Dw4jQEK72gUrfghQ9P48h5x07uFB6owKikQIxKDMLoxECkxihx4nI1/veDrF6NMdedyGBvPHH7CNw4NtIhXbH1BiPOXf0s215pFhfui3EpQRiXEtztDe2mljZ8c6wEXxwqQk5JQ5/bc73m2is4vfkJ8zIDNksM2BJNJbyCIDBgowHDkQGboVmN/OdnW90W89B7UMSN7dVxyfUqapuxL6sce7PKUFDWAG8vGZZMi8e9N6XYVUGiUrfgp3871OOKJV+FDHfNTcaKWUnwsfMCWm8w4qM9efjnV5ehbbP/vVgiABPTQjF/QjTqmrQ4fqkap3Nr7TpG+xgp9nwQ0uoM+N8PT+Hb46V2t03p4wFBQK8uIj3lUjy8JA0rZie5VbfY3mpu1ePYJRUOZ1fi0NlKh4314ipSiYARCQHISA2BwkOKKnUrVOoWVNa1QKVutfozV/p44JcrR/c42LVFFEV8sr8Ar27Jdlg3XokApCcEmgK3tFCMTAxEi9aAgvIG5Jc3ouC6R01D336GPl4yLL8xEStnJ/WporSnmlv12HOyFJ99V+Cwiw6JAJd3pQ5ReiEtTonhsQFIi1MiOdofWTk1+HR/vnnAakfreKHcE1qdATt+uIJN/8mzu9tXkJ8nFk6OwS2ZcUiM9ENRZRPe3n4Re7PKetN8u7SHbeNTg7Hzh2Icu9i7MDYtTonMkeGYNjIc6fEB0OkNKKrUoKiiEYUVTSisaERRRROKVU1Wq+x8vGQIvjoUQIi/J0KUXgjy90RxlQZZOdWdur/ba/ywEPxyxehOA5kbjCLyShtwJq8Wp/NN3cns+ZsfmRCIp1eMxqncGmz5rsChFc+Bfh6YNTYSszOiMC41uNsqw+uJoog9J8vw8qdnUdfY+2DXGh8vGV58dAoyUh075l+rTo9zhWqUqjSoqG1BRW0zKutaUFHbgqq6ln6vthUEIDVGiYyUYIxLDca45GCLisK80gZsPVCAXUdL7A4CxyYH4Rd3jOo0oYAoith5pBj/tyW712MYRgZ7Y+aYCIxOCsKoxECb3Tbrm3R4afOZfnvfmDAsBL+4cxSSe9gV22gUkV/egGMXTTeHT+fWdBnEJ0aahk7JSA3GuJRghCi9TBX8hWp8cagI/znR92o1WxiwdY0BGwM2GoAcGbC1ll1A8RurrG5L/PVuyPzDenVc6ptWnR7fnaqAWqNDRKACiZF+do1RUlXXgn1ZZdhzsszmmDAxoT74n3vGYVyK7Q+DdY1a/Ozvh1BY0ftBtf285ZiSHoaxV++2JUb4WR1Q/nRuDV7cfAYF5c7vsrBkWhx+/ZMxXX5QV6lb8Js3j9p9kRrg64FV85Jx+4xEyGUSHM6uxM4fruDIuSq7PhBPGBaC36wa6/Cue+7CaBRxqbgeh7IrcDi7st8u/t3Vgkkx+OXK0X2q3mrTG/HXT87iy0NFDmxZZ3KZpN8rDz3lUiydHo+75iYj3MoFkVZnQHlt89VZPDWoa9RB6euBsclBV7sj2hdA55TUY9vBIuw+VtKrLlj+PnJEBfsgOsQbUSHeiArxQXSw6fuwQAW0OgOKKk1BiSkwaUJRZSNKVRqHh29Bfp4YkRCI9HhToDY8Vmkz4BJFEdkFdfhkfz72ZzluvKXxqcF4cvkoDIu1v7u+NXqDEXtPluH9b3KsdluTSgRMGxWOxZlxyBwZZvW9+kKRGhu+vNDr8Ks/eHtKMTk9DNNGhmPqyLBOlda26A1GlFY3o7K2GXKZFKFKTwQrvbqtflWpW5GVW42syzU9DtykEgE/mZOEKSPCcDa/DmfyanC2oK5HfydSiYD/WjgM9y9MNf+MDEYRh7Mr8en+fBy/5NjJMZQ+Hpg5NgKzM6IwcXjnLr3tY+0Zro4v+vdPs7utjO+NAF8PvPyzqU6fgbL9dVXUNCOnpAFbvy/o0+e1dl4eUiyYFIPpYyIwJinIropmTUsbvv6xpMs2RAQp8PhtIzE7o+vqroraZjz/4Wm7/5Z9FTLMGR+NhZNjMCYpyO5Ji0RRxO5jJfjrx2f7ZVISqUTAbTMSsOaW4fD3sezRIYoiGpvbUNOgRW2jFiUqDU5c7XHRl6re2DDTTNX5Duj+6+Uh7TKcY8DWNQZsDNhoAHJkwNZ0fh/KP3q68wapDClrf4QgcY9ZnoaSQ2cr8MKm053uFnvIJIgP90VipJ/5kRTpB7lciv2nyrHvZBnO5Ntfsn/HjYl49Nb0TmNLNbW04fH/O4zLxY4NPvx95BibbLrTNjYlCBGBCry5/SK2H3ZsV42empQWgj+vmWQ18DhXUIf/fuuoXXfuA/08cPe8FNw2I8HqxVBNfSt2HyvBjiNXrH4I9VXI8MTtI7E4M86hM725uwaNDip1q/liyGgUYRRF6A2mrwaDEUYjUFajwY8XVDhxudqpA9D3h/BABdben4GM1JAeP7e+SYf/2XgMWTk1/dAy15FJBSycHIuIIAXKqptRWq1BaXVzl93T/bzlmDAsxPQYHoL4cF+Lv50WrR57TpThi0NFOFfY80Ho/X3kWD4zEbfNSLA7IOlI12ZAsUqDogrTrL3FVU24UqnBlcomuwY2l8skGBarxKiEQIxMDMTIhEBEBCl69R5RVdeCrd8X4ouDRTbH1OxOdIg3Hr99JGaOiXDo+5Qoijhyrgpbvy/E5WI1IoK8MWtcJBZOjkGQv33/9scvqvDGlxdwoUjtsHb11NQRYfjJnCRkpIa4dJbM9sDty4NFONnP7xUxoT74w/3jMTIx0OY+eWUN2PJdAb7+saRHVer28JRLIZcJMBiuhWq9CZID/Txw89Q4FFU04kKRutv/98MDFfj7E5mID/ftcj9nMBpFHL2gwkd783oVNMeG+eD2mYm4eUqs3cMEdCSKIk7m1GDrgQIcOF0Bg1GE0scDK+ck4a45yfD0sG98MlEU8cWhIry29ZzVCi6ZVEDmyHAsnByDaaPC+zTuWUVtM9b/O6vf/kaUPh6YPjoc9RqdKVC7Gqq525AZSh8PU9fT1BCMSwlGSrQ/8ssasO37QtNNqQ4/BwZsXWPAxoCNBiBHBmx1Rzah+quXOq2XB8Yg4entvW4j9ZympQ2vbD3n1MApKtgb/333WEwcHgrAdEH61Gs/dBnUjU0Owo3jIvH+NzkO72rhKomRfvh/P52CyGBv87qvfyzGXzad7nZQ3GB/T9wzPwVLp8fb1fW2vYR/5w9XcOJyNTQteswcG4HVi4YjNMDxs7EONnqDEWfz6/DjhSocvVDVbQWct6cUw+MCkB4XcLXqJwD+3nIUVDQiv6wReWUNyC9rRH55g1N/nwUBuHteCh5anGb3hXhBeSN+veHHLgdcBoBJaaGICfFGfrnp9TljAGl3EKL0woRhIRg/LBiXi+ux62hJr6oTokK8cdecZNw8NbZfxs0DTO8D9RodrlQ24UqV5mrw1oTG5jYEK70wMsEUpqXG+Dt8chOtzoBvj5di97ES5JU1QG8wwigColE0fYUIUTRdtLeP3TksVol5E6KxfGaC0ydb6Yn2SXLe/PIiiir7XtFjr0lpIVhzSxpGJwU57Zz2EEUR/zlRhle3ZvfLuFPLpsfjidtH2v130qDRYfexEpzOrUVuWQNKqppc3r0aMHV1fvy2EeZqI1E0VYidL1Lj4hW1eTbq9pnUp4+OwC9XjrZadetqeaUN+HhfPnYfK+kyyBEEYNrIcCy/MRGT00Ltrv6yh7bNgMralj7NDlte04x/fH4eB86Uw2gUkZ4QiIWTYzB3fJRDJtRqZzSK2Lw3D29uv9ht8OXvI8fwWCViQnzwnxNl/TL7qzME+HpgwvAQZKQEIyPVdHPK1s9f06rHN8dK8Pn3hcgtNQ2rwICtawzYGLDRAOTIgE319V+hPvxBp/WKxImIeeDtXreReiYrpxrr3z/V71OQ27L0hng8vCQN69472eXdz+GxSrz682nwVcjRotVjy4ECfPhtXq+rIbqTHh+AJ24fiZIqDXYdLe7XO/HB/p548dEpGBarxIYvzuPD/+R1ub+nXIpHb03DsukJdt+ZJcera9Ti2EUVfrxQheIqDQQBSIlWYkR8ANLjAxAf4Wf3B/zaRi3yyxqQV9aIs3m1OJlT3acuG/ZIjfHHuv+a0GkMpOsZjCKOnKvEuvdOdlu9d9fcZDy2bIT5NYuiCJW69VqYWN6AnJIG8wfl/hTs74mYUB+EKL1w5Fylywd2786I+ACsmpeCG8dFDoqxD4cyURRRVNmE/LJGSCSAXCqBVCqBTCpA1uFrcZUGe7PKcOhsZY/HLBqfGowHb0lz+PhbjqZp1eOfX13CJ/vyHdJNONDPA7+9exymj47o03FadXoUlDchr7QBuWWm96W80oZ++0zRUUyoD36zaiwmDOu+mlgURdQ1aiGTSXo9kZMz1Ta0Yuv3hdh6oNDi/zE/bzmWZMbhtpkJiA5x/6EoWnX6q3+r/VsRmlfagJc/yzZPYBIW4IVhsUoMi1GavsYqER54rXpY3aTF29sv4otDRW4REttjcnoolk2Px/TRET3+92wfbmDbwUJ8+c0RHP/wcfM2BmyW3DJgS7lkAx0AACAASURBVEpKcsp59Ho9SkpKADBgo4HFkQFb+Ue/QtP5PZ3W+2UsQcTtz/W6jWQfbZsBb355AR/vy3fqzFzWdDfeUmKEH/7x1LROdw41rXp8si8fm/fkOexuno+XDI8uTcey6QkWF7rlNc345lgJdh0t6bY6wcdLhozUEExKC8HE4aHILW3An9/P6rIizVMuRVqcstvZscICvPDCI5OR5uRxV8i5RFFEQXkjsnJqcDKnGlk5Nf0SuHnIJfjZshGYOyEaJVXXKpqKqzQoqmxCqUrTbSWlTCrgmbvGYnFmnF3nrGs0TQpy7KIKRy9UoUrd/UzB1vgp5IiP8EVcmC9iwnwQG+qDmDBfxIT6wOe67ucNGh0+3V+AT/bnm6tA3MX00RFYNS8ZY5ODhlT3bLLUqtPjyLkqu8K2sclBWLM4za5gxp3klzXgb5+c7fXNKqlEwPyJ0Xj89pEI6qfJSUzVY1pk5VRj/6lyHD5X6fCZu6USAXfPS8HqRcMG/Q0ybZsBh7MrUVjRiPhwP0wbFWZXtf1Q1djcBkGA3eOk5pTU4+VPs3Eqt39uACs8pBiXGgyljwdO59X2+EZ8sL8nbsmMw63T4hDloED1h2OnkDk5w7zMgM2SWwZsEokEgiDAGU1rPw8DNhpIHBmwXdlwN7Sl5zutD5r1EILnPtbrNlL3Ll5R47l/nXTIwLRdSY32R2SINw6cruj1MaJCvPHGU9O77MLY2NyGzXvz8PG+/D6NkTVvQjSeXD6yyzGPRFHExSv12HW0GP85UYq6Rh3kMglGJwVh4vAQTBoeiuFxyk536M7k1eA3bx7r093xUYmBeP6hST2eNY8GPqNRREFFI05ersbJy9UoVmkgEQSEBSoQHuCF0EAFwgK8EBaoQFiAAqEBXigob8Qf/3USxb2c5c8eAb4eeP6hSRjbxcQlXWmv9jGFbSpk5VR3qjbz8ZJZjP1o+t4fIUrPHoVSmlY9vjhYiI/25PV5JtK+UPp44Oapsbj1hni3GEOJ3IutsG1UYiDW3DIck9JCB2wYK4oivj1Rite2nuu226hcJkF6XADGJAdhbHIQRicFdRq0vb81t+rxw/kq7Msqw+HsSrtnrLQlLS4Av717LFJj+jZBB1E7URSx92QZXvv8PCrrWvp0LKlEwKjEQEwcHopJaSEYkRBo8Vm2orYZWTk1OJVrelj7bCEIwOS0UCybnoAbRoc7vPrPkdehg5FbB2zOwoCNBhpHvrHlvzAXBk3nap2wZWuhnHBbr9tItukNRvx7dw7e/fpyt1012t8Ke/pOnRLtj9kZUZgzPsp88fjDuSq8sOlUjytVQpRe2PD0DXbf+WrR6nHy6n/+p3NrcKFIbVeXlOgQb/zqJ2MwJb1nM9eKooiaBi2UPh52jWVVUtWEX77xY68Cj5unxuKZn4xx6zGIyP20aPV4des5bDvo+Fk/k6L88NKjluMH9lWb3ogLRXWortfC21OGpCg/hAZ4OfSzmbbNgJ1HruDD/+TZdUdeEICwAAUC/TyRX9bQbUWfLRmpwVh6QzxuHBfZp8GxaejQ6gwoVjXBx0vu0L8zV9O0tOGdry5hy4FCc/W6n0KO0clBpkAtKQhp8QFu9XfSqtPjx/Mq/H/27jwuqqr/A/hnZpB9FxBBgRBXFAXFFfeyMkv7WRm5W6Zmm2bZYopmVlpZpj5WWvpkPmqmZj5pmaGJaAKKGu6KC4uyy77N3N8fPEyO3IszMMzMhc/79eLV3HvOPec7qOR8/Z5zohPTEHP6lkH/mNfCzQ5jhgTiyUGBXAJODaK0vBLf77uM7/Zd1Lvq0sm+GXw97BHa1gM92nuga5vmNQ4fq03W7VLt37ez88sQ5OuMh3q2Mlq1mhgm2Gpn8Qm2hg6PFWwkR8b6waapKMXlRX1E23wn/Qv2bXrXOUbSVVhSgeMXshB/PguxSbfuuUk5UHUAwbwJoejo54qrtwqRnJ6P5PQCJKcX4Ep6AdKzi3USb4E+Thga5oPBoT4I8Bbfz6mopAKrdp7R+4O+q6M1Vs/qJzmePkrKKpF0Nfd//9qWg6TkXJ1TxJpZKTH+gSCMH9bWZEs1bheW482vjt1zKWg1pQJ48f+CMWZwoGyrFsj8Yk7fxJKNiUZbZhrRpQUWTOqusxRTbirVGuyLT8XOmKtIzSyGu7MNfDzs4dPcHr6eDvBt7gBfT3t4u9trE+hl5WqcvpKD+PNZiD+fiXPX82rdA8fZoRke7tUaI/v51+tnGVFjlF9cjsup+XB1tKl1s3NLU1ahxrn/nfapUilgpVRApVRCpVJApfzfl6rqnoOtFVp7OfD/32QSN3OKsTX6Cq6kF8DexgruzjZo7mwDd2fb//3XBu5OVV9y/AdbJthqZ5EJNhcXFxQWFkIQBAQFBeHrrxtmo/X09HQ888wzALgHG8mLsX6wlWdexbUV4lVq/q/shLWHf51jbOoqKjX4OzkX8eczEXcuU+8qrmoj+/njxf8LrvWDc2l5Ja7eLERpuRperrYG/WtV/LlMfLDpZK2VI452VvjilX5o39q4yygqKjU4dz0PF1JuQwGgf0hLs5yeWV6hxpLvE/FbXGqt/RztrLBoSg/07mRYZR2RmJz8Uiz5/iRi/75Vr3HGPRCEaY91ZCUGqv4B48TFbMSfz0T8+SwkpxdAoQBCAt0xKiIAg0JZrUZERGQMTLDVziL/yTMsLAwHDx4EAKSkpCAiIgIqlfH/YnTtmvGXahDJScVt6T25rFxamDAS+SmrUCOvsBy5BWV3fFVdJ98sQOLF7DrtE9Lc2QZvje2Gvp3v/f23tbaq8yb7PTp44ru3B2HNrrPYdjBZZGwVPp7Ry+jJNQDavdK6BLobfWxDWDdTYcHEMLTydMA3v1wQ7ePn5YCPpvfiHk1kNO7Otlg2vSd2xlzDih+TdKo59dH5PjdMfLAt+tXz9L7GxNGuGfqHeKN/SNX3pKSsEoIAg5bZEBEREdWXRf7NIzw8XJtgKysrw+nTp9GtWzczR0XU+FTmpYneVzm4Q9mMG7hXKytX40BiOvYfT0Nyej5yC8vrtYm/lKFhPpgzJgQujqbZQNje1gqzn+qCIWE+WPafU0i+WQCgap+SBZPCENKmbhumy4lCocBzj3SAr4cDPvg+EZXqf6oMe3X0xKIpPeBkr99JUkT6UigUeLx/AMLaeWDxdyeQlJyr025vawU/L0f4eTmgtZcjWrdwgL+XI1p5Ocp6Oaip2Nnwe0RERESmZ5F/A+nRo4fOdXx8PBNsRA2gUqKCzcq1pYkjsUwXU27j59jr+PVYCgpKKhpsHif7ZpgzJgQP9PBtsDlq0y2oOTbOG4Rz129DrdGgjY9zk/uA+nCv1ghp444fDyajsKQSfTt7YUBIS9nsRUPy5N/CEV+9FoG/k3ORkVsCd2db+LdwgJuTYSdzEhEREZH5WeQnqPDwcADQ/uUyLi4Ozz33nDlDImqUKvLSRe83a8IJtqKSCuyLT8Wu2Os4dz2vQedSKIBB3Vri1Sc6w9PVrkHnuncsCnT0r9ty08bC18MBL4/ufO+OREakUCjMvlyaiIiIiOrPIhNs9913H9zd3ZGbW7VkIj4+3swRETVOlRIJNiuXprW3jyAIOH0lB7tir+OP42korcPeafrycLFFeAcP9GjvifAOnvBw4VJcIiIiIiIiubPIBBtQtUz0t99+A1B1UkV5eTmsrY27L5GVlRX8/Py4DIOaLKlDDppSBVteYRne+/cJHEnKaJDx7W2tENa2OXq090TPjp7wb+HInzlERERERESNjMUm2N566y08/PDD2uuSkhKjJ9h8fX1x9epVo45JJBeCRo3K27dE25rKHmyl5ZV4ZcURXEzNr9c4KqUCro7WcHW0hpuTDdycbBDg7Yge7T3R0d8VViqlkSImIiIiIiIiS2SxCbaBAwdi4MCB5g6DqNFSF2YDGvGTMK1cGn+CTRAELNt8yqDkmpNdMwwL90W3oOZwd65KpLk52sDJvhk3wyciIiIiImrCLDbBRkQNS2p5KAA0c238e7DtjLmGPX+l6NU3rG1zPNrXH4O6tYSNtaqBIyMiIiIiIiK5YYKNqIlSF0uckKmygtLOxbTBmFhSci6W/3C61j7NnW0wvHdrjOjjh9ZejiaKjIiIiIiIiOSICTaiJkojkWBT2bs26k34cwvK8M7aOFSqBdF2L1dbzH6qC/p2bsG904iIiIiIiEgvTLARNVFSFWwqO1cTR2I6lWoN5n+TgIy8UtH2ZlZKLJkajk4BbiaOjIiIiIiIiOSM5RlETZS6+LbofZVD402wfb37HBIuZEm2v/pEZybXiIiIiIiIyGBMsBE1UZIVbPaNM8F28GQ6vvvtkmT78F6tMSrC34QRERERERERUWPBJaJksbKzsxETE4P09HTk5OTAzc0NPj4+CAsLQ+vWrY06V3FxMQ4fPoyUlBTcunULrq6u8PX1RXh4OLy9G+eJmlIJNqV94zvg4PqtQiz+9wnJ9ratnPH60yGNeu85IiIiIiIiajhMsJHF+e9//4slS5bg6NGj0Gg0NdoVCgX69OmDV155BU899VS95kpOTsb8+fOxY8cOFBUV1WhXqVQYMmQI3nrrLQwePLhec1ka6UMOGtcSyZKySrz1dRyKSitF253smmHJ1HDYWKtMHBkRERERERE1FlwiShajoKAAI0eOxIgRIxAbGyuaXAMAQRAQGxuLMWPGYPTo0cjNza3TfOvXr0dISAg2btwomlwDALVajX379mHo0KGYPXs21Gp1neayROoiqQRb46lgEwQBH246ieT0Ask+8yeFwdfDwYRRERERERERUWPDCjayCIWFhejfvz9Onjypc9/e3h69e/dGixYtkJ+fj7i4OGRkZGjbt2/fjoyMDOzbtw+2trZ6z7dp0yZMmTIFgiBo71lZWSE8PBytW7dGZmYmEhISkJ+fD6AqUbN8+XKUlZVh1apV9Xy3lkFdInHIQSPag21L9BXsi0+VbJ/ycDv069zChBERERERERFRY8QKNrIIEydO1EmuWVlZ4b333kNmZib279+PTZs2Yffu3UhPT8fmzZvh6emp7RsTE4Pnn39e77mOHz+OyZMn6yTXRo4ciStXriA2NhZbtmzBH3/8gZSUFLzzzjs6z65evRpfffVVPd6pZRAEQfoU0UZQwSYIAtbvvYAVPyZJ9undyQuTh7c3YVRERERERETUWDHBRmb366+/Yvv27dprhUKBzZs3Y968ebC3t9fpq1QqMWbMGPz5559wc/tnr7DvvvsOR44c0Wu+N954A+Xl5drrJ554Atu3b69xcIKTkxMWL16Mzz77TOf+vHnzUFAgveRQDjRlhYBGfE8yue/BptYI+GTraXz18znJPt7udlgwMQwqJQ81ICIiIiIiovpjgo3M7oMPPtC5njJlCkaPHl3rMx06dMDy5ct17s2dO/eec0VHR2P//v3aaw8PD6xZswZKpfQfhZdffhmDBg3SXmdmZtaYW26kThAF5L1EtKxCjXfXxWP7n1cl+1hbKbFkajhcHK1NFxgRERERERE1akywkVnl5eUhJiZG597rr7+u17Njx46Fj4+P9vrQoUO4fPlyrc/8+9//1rl+7rnn0Lx581qfUSgUeOONN2odR240EstDAUAp0yWi+cXlePWLIziQmF5rvzljQtDBT75JRCIiIiIiIrI8TLCRWR0+fFjnZE5/f3+0b6/fvlhWVlYYOnSozr1t27ZJ9ler1fj555917k2ePFmvuR588EG0bNlSe3358mWcOnVKr2ctkbpY4uRVpRWUNo6mDcYIbuWWYManh3Hyck6t/V4eHYwRff1MFBURERERERE1FUywkVmlpKToXHfu3Nmg57t06aJz/d///leyb1xcHLKzs7XXLVu2RLt27fSaR6lUYsCAATr39uzZY0CklkVdJL5EVGXvAoVCXvuSXUnLx/MfH0JyuvS+eFYqBRZO7o6nh7QxYWRERERERETUVFhkgq1NmzZ4/fXXcfjwYXOHQg0sJ0e34sjV1bCle3f3P336tGTfv//+W+e6T58+Bs3Vt29fneukJOkTKi2dukTqBFF5LZ1MvJSNGZ8eRmZeqWQfe1srfDqzNx7o4WvCyIiIiIiIiKgpscgEW3JyMj799FMMGDAA3t7emDZtGvbu3YuKigpzh0ZGZm2tu9F8WVmZQc/f3T8vLw/p6eJ7cJ05c0bnOigoyKC52rTRrX66ezw5qa2CTS4OJKbh1S+OoKBE+ueCh4sN/jWrH3q09zRhZERERERERNTUWGSCrZogCMjIyMDatWvxyCOPwNPTE5GRkdi6dSsKCwvNHR4Zwd0HDEglx6SI9T9//rxo30uXLulc+/kZthfX3f0vXrxo0POWRF0inmBTyqSC7djZDLy7LgHllRrJPn4tHPHla/3RtpV8koZEREREREQkTxaZYJsxY4bO6ZCCIEAQBOTn52Pr1q2IjIyEp6cnRowYgXXr1iEzM9OM0VJ9dOzYUef6xIkTKC8v1/v5I0eO1Lh3+7b48se8PN2kkpeXl97ziPUvKCiARiOd4LFkGskKNjcTR2K41KwivPtNAtQaQbJP8H1u+HJ2BFo2tzdhZERERERERNRUWZk7ADGrVq3CqlWrEBcXhx07dmDHjh3aqiRBEKBQKFBWVoY9e/Zgz549mD59Ovr06YPHH38co0aNwn333Wfmd0D66t69OxwdHbUVicXFxfjxxx8RGRl5z2evXbuGgwcP1rhfUCC+2f3dVY92dnYGxXp3f0EQUFRUBCcnJ4PGuVtGRobBSeK7q/EMJb0Hm2VXexWXVuLNL4+hoFh6WWhElxZYNKU7bK0t8scbERERERERNUIW/Qk0PDwc4eHhWLJkCc6dO4cdO3Zg586diI+PB/BPsk2tVuPw4cM4fPgw5syZg5CQEG2yLSQkxMzvgmpjZWWFMWPGYN26ddp78+bNw6OPPgpHR8dan501a5ZoBZm+CTZbW1uDYhVLyBUWFtY7wbZ69WosXLiwXmMYSnoPNstdIioIAt7feAKX06RPC32snz/mjOkCK5VFFucSERERERFRIyWbT6EdOnTAW2+9hb/++gvXr1/HF198gSFDhkClUgGo+vBd/d9Tp05h4cKFCA0NRZs2bTBnzhyeSCrhxRdfhEKhaPCvqKgoyRjmzp0LK6t/cr1XrlzBY489JrnUU61WY/bs2dixY4dou1Kp329rhUKhV7+69rdk6mL5HXLw3W+XEH1Ceo++Jwbeh7mRIUyuERERERERkcnJ8pOor68vZs6cid9//x23bt3C+vXrMWrUKG1F0p3JtuTkZCxfvlznRNI9e/bwRFIL0rZtWyxdulTnXnR0NNq1a4dFixYhJiYGFy9eREJCAlavXo3Q0FAsX74cQFXS6879+gDA1VW8CuvuiriSkhKD4hTrf68qO0skCILkIQeWugdb7N+38OXPZyXbe7T3wMujgxtVEpSIiIiIiIjkw6KXiOrDzc0NEyZMwIQJE1BSUoJff/0VO3bswO7du5Gbmwvgn4Rb9Ymka9euhaOjI4YPH47HH38cw4cPl2WipDGZNWsW0tLS8PHHH2vvZWRkYMGCBViwYIHkc5988gl++uknpKWlae+ZMsHm4OBg0BhiXnjhBTz55JMGPXPp0iWMGjWqTvNpyooAdaVom9ICK9huZBQian0CBIkzDVo2t8d7U3qwco2IiIiIiIjMRvYJtjvZ2dlh1KhRGDVqFNRqNQ4cOICdO3fip59+QkpKCoB/km0FBQXYunUrtm7dCmtrawwdOhSjRo3CyJEj4enpac63YVIjR45Eq1atGnyeiIiIe/ZZtmwZunXrhrlz5yI1NbXWvp6enli9ejWeeOIJrF69WqfN29tb9BkXF93kkaEHC2RkZOhcOzs7670ctTZeXl4Gn2haHxqJAw4Ay9uDrai0EnO/PIbCEvGEoK21Ch8+Hw4XR2sTR0ZERERERET0j0aVYLuTSqXC0KFDMXToUHzxxRd6n0g6Y8YM9O3bFwsWLMCQIUPM/C4a3gMPPIAHHnjA3GFojR07Fk888QS2bt2KvXv3Ii4uDpmZmSguLoa3tzfatWuHJ554AmPGjIGrqytKSkpw7do17fN2dnYIDg4WHbtt27Y613c+p4+7+989nlxIHXAAWNYSUY1GwHv/Po6rNwsl+7wzrhvatrK8qjsiIiIiIiJqWhptgu1ud59IunPnTuzYsUP0RNKYmBjExMQ0iQSbJbKxscH48eMxfvz4e/aNj4/X2U8vLCxM58CEO3Xs2FHn+tKlSwbFdeXKlVrHkwup/degVEFpazlLpdfvvYA/T96UbB8/LAhDu/uaMCIiIiIiIiIicU1y06IOHTrgzTfflDyRlORj//79OteDBg2S7Nu5c2ed6yNHjhg0190n0d49nlxIVbCp7Fws5pCAQ6duYu1/z0u29+7khecflWeCk4iIiIiIiBqfJplgu9OdJ5JmZGRg/fr1GDlyJOzt7c0dGt2DRqPB+vXrtdcKhQJTpkyR7B8eHg53d3ftdXp6Oi5cuKD3XIcOHdK59/DDDxsWsIVQF0udIGoZSy2v3izAwg3HJdt9PewRNTkMKqVlJAOJiIiIiIiImnyC7U6urq6YMGECtm/fjjlz5pg7HLqH//znPzr7oj3wwAMIDAyU7G9lZYVHH31U5963336r11y//fabzkmlbdq0QUhIiIERWwaNRIJNaSEHHCzbfArFpeKHGthZq/DRtJ5wtuehBkRERERERGQ5mGAjWcrIyMDs2bO110qlEu+99949n5swYYLO9dq1a5GdnX3P55YuXVrrOHIiXcFm/gTbyUvZOHFR+tfj3YmhCPRxNmFERERERERERPfGBBtZhMpK8YolMTdv3sSQIUOQkZGhvTdr1iz07Nnzns8OGTJE5/CKrKwsTJ8+HRqNRvKZFStWIDo6Wnvt4eGBWbNm6R2vpVEX3xa9bwkJto37pA+emPxwOwzq5mPCaIiIiIiIiIj0wwQbWYR58+Zh5MiR2Lp1K27fFk8AZWdn4/PPP0enTp2QlJSkvd+9e3e9qteqLVu2DNbW/ywx3LZtG0aPHo0bN27o9CsoKMC7776LV199Vef++++/DycnJ73nszSWWsF2OS0fh/++JdrWtY07nh3e3sQREREREREREenHytwBEAFVFWy7du3Crl27oFQq0a5dOwQFBcHFxQX5+flIS0tDYmIi1Gq1znNhYWHYs2cP7Ozs9J4rLCwM33zzDcaNG6e9t3PnTuzevRs9e/ZE69atkZWVhbi4OOTn5+s8O2PGDDz//PP1e7NmZqkJtu/vUb2m5KEGREREREREZKGYYCOLo9FocO7cOZw7d67Wfs8++yw+//xzODg4GDzH2LFjUV5ejpdffhmFhYUAqpJ8sbGxov0VCgVefvllfPLJJwbPZWksMcGWnl2MffGpom3tW7sgvIOniSMiIiIiIiIi0h+XiJJFGD58OEaMGAFHR8da+9nY2OCpp55CXFwc1q5dW6fkWrXJkyfj5MmTGDt2rOQ4SqUS999/P/bv34/PPvsMKpWqzvNZAkEQajlF1MXE0fzjP/svQ60RRNvGD2sLhYLVa0RERERERGS5WMFGFqH68AG1Wo2kpCScP38eqampKCwshEqlQvPmzdG+fXv07NnToOWg9xIYGIiNGzeiqKgIMTExSElJQUZGBlxdXeHj44OePXuiZcuWRpvP3ITyEgjqCtE2c1Ww5RaU4efY66JtrTwdMLBb4/n+ExERERERUePEBBtZFJVKhZCQEISEhJh0XgcHBzz44IMmndMcpJaHAoDKwTwJth8PJqOsQi3a9sz9baDi3mtERERERERk4bhElKgJqTXBZmf6BFtxaSV+OJgs2tbc2QYP92pt4oiIiIiIiIiIDMcEG1ETIplgUyihtHUybTAAdsVeQ0Gx+JLVMUMCYdNM3nveERERERERUdPABBtRE6Iuvi16X2XnDIXStD8OKio1+M/+y6JtjnZWeDwiwKTxEBEREREREdUVE2xETYimOFf0vtLezcSRAL/FpSAzr1S07fH+AXCwa2biiIiIiIiIiIjqhgk2oiZEaomoyt7FpHFoNAI27rsk2mZtpcRTgwJNGg8RERERERFRfcgywXb27Flzh0AkS5JLRO1Ne8DBodM3ce1WoWjbI71bo7mLrUnjISIiIiIiIqoPWSbYgoODMXjwYGzZsgUVFeIbpBNRTZIVbA6mS7AJgoCNv4lXrykVQOT9QSaLhYiIiIiIiMgYZJlgA4A///wTzzzzDFq3bo23334bV69eNXdIRBZPMsFmZ7oEW+KlbCRdFd8LbkiYL1p5OpgsFiIiIiIiIiJjkG2CDaiqhMnIyMBHH32EoKAgPPLII9i9ezcEQTB3aEQWyRIq2L6TqF4DgLEPsHqNiIiIiIiI5EeWCbawsDBtEk2hUAAANBoN9u7di5EjRyIgIADvv/8+bt68ac4wiSyORmIPNqWdaQ45uHDjNo6eyRBt69XRE+1bm/awBSIiIiIiIiJjkGWCLT4+HvHx8ZgyZQrs7e0hCII20SYIAm7cuIH58+fD398fTz31FPbv32/miInMTxCEWk4RNU0Fm9TJoQAwflhbk8RAREREREREZGyyTLABVVVsa9euRVpaGlasWIHg4GCdqjZBEFBRUYEff/wRw4YNQ/v27bF8+XLk5orv/UTU2AkVpRAqy0TbTJFgu5Vbgj+Op4q2BQe4IbRt8waPgYiIiIiIiKghyDbBVs3JyQkvvvgiTp06hUOHDuGZZ56BtbU1AOhUtV28eBFz5syBr68vJk2ahKNHj5ozbCKTk6peA0yTYDtwIg0aie0Rxz0QpP3zSkRERERERCQ3sk+w3alfv37YuHEjUlNTsXTpUgQFBdWoaistLcV3332Hfv36oVu3bvjqq69QVFRk5siJGl6tCTYTHHJwIDFd9L6flwP6h3g3+PxEREREJpd1cQAAIABJREFUREREDaVRJdiqubu7Y86cOTh//jz27duH//u//4OVlRUA3aq2U6dOYcaMGfDx8cHMmTNx6tQpc4ZN1KDUEgccQKGA0tapQefOvl2KU1dyRNse7NkKSiWr14iIiIiIiEi+GmWC7U5Dhw7Ftm3bcP36dSxatAh+fn7aqjagKtFWUFCANWvWIDQ0FP369cP333+PiooKM0ZNZHyaYvH9B5V2LlAoVQ0696FTNyFILA8d1M2nQecmIiIiIiIiamiNPsFWrUWLFpg3bx6uXLmCXbt2wc/PDwqFQvslCAIEQcDRo0cxYcIEtG7dGosXL0ZhYaG5QycyCqkKNpWdS4PPfeCkxPLQFo4I8HZs8PmJiIiIiIiIGlKTSbABQG5uLj777DO8/vrruHHjhk7bnRusC4KAjIwMLFiwAIGBgdi8ebOpQyUyOqk92Bp6/7X84nIknM8SbRvUtSUPNyAiIiIiIiLZszJ3AKYQGxuLNWvWYNu2bSgrK4MgCDp7sQFA586d4ezsjNjYWAD/HIqQlZWFsWPHIikpCe+9957Z3gNRfUkm2Bq4gu3w6VtQSxwfOii0ZYPOTURERERERGQKjbaCraCgAKtXr0bXrl3Rv39/fP/99ygtLQXwT/JMpVLhySefxIEDB3Dq1CnExMQgKSkJM2fOhL29vU7fJUuW4ODBg+Z8S0T1Il3B5tag80qdHtrCzQ7tWzf88lQiIiIiIiKihtboEmwJCQmYOnUqfHx88NJLL+H06dM1DjVo2bIloqKicP36dWzZsgUDBgzQtnfs2BFffPEFrl69ikmTJulUu61cudLk74fIWCQTbPYNt0S0uLQSf53NEG0b1I3LQ4mIiIiIiKhxaBRLRIuLi7Fp0yZ8+eWXOH78OIB/ln7eeYDBwIEDMXPmTDz++ONQqWo/NbF58+b45ptvUFhYiG3btgEAjhw50rBvhKgBaSQOOVA24BLRo2cyUF6hEW0b1I3LQ4mIiIiIiKhxkHWC7fTp0/jyyy+xceNGFBQU6CTVqhNrjo6OGD9+PF544QV06tTJ4DmmTZumTbBlZmYaNX4iUzLHIQdSy0PdnWzQOdC9weYlIiIiIiIiMiVZJti+++47rFmzBkePHgUgXq3WqVMnvPDCC5gwYQIcHR3rPFdAQID2dWVlZb3iJjIn6UMOGibBVlahRuzfN0XbBnT1hkrJ5aFERERERETUOMgywTZx4kRtMu3OajWVSoVRo0Zh5syZGDhwoFHmsrKS5beISIemohRCRaloW0NVsMWfy0RxmVq0jctDiYiIiIiIqDGRffao+tCCqVOnYtq0aWjZ0rgf3G1tbY2WrCMyF7XE/msAoGqgPdikloc62TVDWDuPBpmTiIiIiIiIyBxkm2ATBAEDBgzQHlrQUJVmLVq0QHR0dIOMTWQqGonloQCgbIAKtkq1BodOiy8PjQhpAStVozvAmIiIiIiIiJowWSbYpk+fjpkzZyI4ONjcoRDJgro4V7xBoYDK1tno8yVezEZ+UYVo26BuPkafj4iIiIiIiMicZJlgW716tblDIJIVqSWiSlsnKFTG/zFw4KT48lA7axV6dvA0+nxERERERERE5sR1WkRNgOQJovbGXx6q0Qj4UyLB1ie4BWysVUafk4iIiIiIiMicmGAjagKkE2zGP+Ag6Wousm6XibYNCuXpoURERERERNT4MMFG1ARIJ9jcjD6X1OmhzayU6BPcwujzEREREREREZkbE2xETYBGag82I1ewCYKAgxIJtp4dPOFgK8ttH4mIiIiIiIhqJctPu1OmTDHaWAqFAk5OTnBxcYG3tze6d++Obt26wdra2mhzEJmbqfZgu5iSj7TsYtG2Qd24PJSIiIiIiIgaJ1km2NavXw+FQtFg49va2mLMmDF46aWXEBoa2mDzEJmKqRJsUstDVUoFIrp4G3UuIiIiIiIiIksh+yWigiDofBnS9+7+1fdKSkqwYcMG9OzZE++88w7UanVDvgWiBmeqBJvU8tDQts3h4siqUCIiIiIiImqcZJtguzNBplAotF9iSTSpvneOc2eyrbpNrVbjww8/xOTJk0387oiMyxQJtqs3C5B8s0C0bSCXhxIREREREVEjJsslosnJyQCAU6dO4dlnn0V2djYEQUBQUBCefPJJhIeHw8/PD87OzigvL0dOTg5Onz6N6Oho7Nq1C+Xl5VAoFIiMjMTixYtRVlaGvLw8nDlzBn/++Sd++OEHlJSUaBN233//PXr16oWZM2ea+Z0TGU5TUQahvES0TWXEQw6kqtcUCmBgVybYiIiIiIiIqPGSZYLN398fv/zyCyIjI1FSUoLmzZvj888/R2RkpOQzERERmDFjBm7duoWXXnoJ27Ztw+bNm5GZmYlffvkFVlZW6N27N6ZMmYJPP/0UM2bMwA8//KBNsi1evBjPPfccbGxsTPhOiepPUyJ+gigAKI1YwXbwpHiCrfN9bvBwsTXaPERERERERESWRpZLRK9evYqxY8eiuLgYHh4eOHjwYK3JtTu1aNECW7duxUsvvQRBELB//37MmTNHp4+7uzu2bNmCcePGaZeOZmRk4McffzT6eyFqaFLLQwHjLRFNzy7GueviiTxWrxEREREREVFjJ8sE29y5c3H79m0oFAp8+umn6Nixo8FjfPLJJ+jYsSMEQcAXX3yBs2fP1uizevVquLu7a/dki46OrnfsRKZWa4LNztkoc/wpUb0GAIO4/xoRERERERE1crJLsOXl5WHXrl0AqirN9K1cu5uVlRWef/557fWGDRtq9HF0dMSECRO0VWxxcXF1movInNRF4gk2pa0TFKpmRpkj7nyW6P12rV3g4+FglDmIiIiIiIiILJXsEmyxsbEoKyuDQqFAjx49oFTW/S306dNH+/qPP/4Q7TNo0CAAVaeNZmRk1HkuInNRlzTsCaJqjYCTl7JF2yK6tDDKHERERERERESWTHYJttTUVO3r5s2b12ssd3d30XHv5Ofnp32dm5tbr/mIzEGqgs1YCbZLKbdRVFop2ta9nadR5iAiIiIiIiKyZLJLsGVnZ4u+roucnBwAVdVp1a/v5uTkpH1dvVSUSE6kThFV2rsYZfzjF8X/HFpbKdEpwHinlBIRERERERFZKtkl2Dw9qypiBEFAfHw8NBpNncc6cuSI9rVUNVxhYaH2tYMD95Ii+WnoCrYTEgm24PvcYNNMZZQ5iIiIiIiIiCyZ7BJsgYGB2tc5OTnYsmVLncZRq9X4+uuvAQAKhQJt2rQR7Xft2jVtH29v7zrNRWRO0nuw1b+CTaMRcPKyeIItNKh+S7iJiIiIiIiI5EJ2Cbb+/fvD1dUVCoUCgiBg1qxZuHTpksHjvPHGGzhz5oz2+tFHHxXtd+fJoVJJOCJLpi6WSrC51Xvsy2n5KCiuEG0LbedR7/GJiIiIiIiI5EB2CTYrKytMmTIFgiBAoVAgIyMDERER+PHHH/V6PisrC+PGjcNnn30GhUIBAHB0dMT48eNF++/cuVP7ukePHvV/A0Qm1pBLRKX2X2tmpUTngPon8IiIiIiIiIjkwMrcAdRFVFQUNm/ejPT0dG2S7amnnkK7du3w5JNPokePHvD394eTkxPKy8uRm5uL06dPIzo6Gj/99BPKysq0BxYoFAosXLgQLVq0qDHPX3/9haSkJG0ibuDAgSZ9n0TG0JCHHCRezBK938nfFTbW3H+NiIiIiIiImgZZJtgcHR3x66+/YvDgwcjOztYuFz1//jzef//9Wp+trnyrfmbGjBl49dVXRfsuXLhQ+4yHhwf69+9v9PdC1JCEygpoyopE2+pbwabRCEi8JH76bre23H+NiIiIiIiImg7ZLRGtFhwcjEOHDqF79+46STOgKiEm9gVAm1hr1qwZPvjgA6xcuVJyjl9++QUajQYajQYZGRlQKmX77aImSuqAA6D+CbbkmwW4XVQu2hbWlvuvERERERERUdMh64xR+/btcfToUaxevRqdOnXSSaSJEQQBNjY2mDhxIk6cOIG5c+eaMFoi05Pafw2of4LtxAXx/ddUSgU638f914iIiIiIiKjpkOUS0TsplUpMnz4d06dPx+nTpxEbG4uTJ08iKysLeXl5sLGxgZubG/z9/dG7d29ERETAxaX+e08RyYFaYv81AFDZ1e/PwYlL4vuvdfR3hZ2N7H+0EBEREREREemtUX0K7tKlC7p06WLuMIgshrooV/S+0sYRCqtmdR5XEAQkXhKvYAvl/mtERERERETUxMguwaZWq1FU9M+m7Y6OjtwbjUiCulh8iWh9l4deu1WI3ALx/ddCuf8aERERGYlGo0F2tvg/6hERkf6aN2/O3EkDk12CbcOGDZg6dSoAQKVSITk5Gb6+vmaOisgyaYrFl4gq7eu3PPT4BfHloSqlAl0C3es1NhEREVG17OxseHl5mTsMIiLZy8jIgKenp7nDaNRkl2C7deuW9iCDsLAwJteIatFQFWwnJJaHtm/tAgdb2f1YISIiIiIiIqoX2dUHOjg4AAAUCgX8/PzMHA2RZWuIBJsgCEi8yP3XiIiIiIiIiKrJLsHm7e1t7hCIZKMhEmw3MoqQnV8m2taN+68RERERERFREyS7tVzt27fXvr5x44YZIyGyfOoS8T3YVPXYg+34RfH915QKoGsb7r9GREREDevMmTPw8OA/6hERScnKykKnTp3MHUaTI7sEW9euXREUFIRLly4hISEBubm5cHNzM3dYRBZJUyRewaasRwWb1PLQdq1d4GjXrM7jEhEREenDw8ODG3UTEZHFkd0SUQDaU0TVajU+/PBDM0dDZLnUJcZdIioIguQBB92CuP8aERERERERNU2yTLC99tpr6N27NwRBwKeffootW7aYOyQiiyOoK6ApLRRtq2uCLTWrGJl5paJtodx/jYiIiIiIiJooWSbYlEoltm/fjl69ekGtVuOZZ57B888/j8uXL5s7NCKLoS4W338NqHuC7YTE/msKBdA1iPuvERERERERUdMkuz3YAGDRokUAgMGDB+PMmTMoKCjAunXrsG7dOgQHByM0NBReXl5wcnIyaNz58+c3RLhEZlF7gq1uhxyckNh/LcjXGc721nUak4iIiIiIiEjuZJlgi4qKgkKh0F4rFAoIggAA+Pvvv5GUlFSncZlgo8ZEXZwr2VbXQw6kEmxcHkpERERERERNmSwTbGLuTLgZShCEej1PZIk0EhVsCmt7KK0MrzZLzy7GrdwS0bZQHnBARERERERETZhsE2zVFWtEJE5dbNwTRI9fEN9/DeAJokRERERERNS0yTLBFh0dbe4QiCyeZILNoY7LQy+JLw9t4+MEF0fuv0ZERERERERNlywTbAMHDjR3CEQWTzLBZle3BFsi918jIiIiIiIiEqU0dwBE1DCMuUT0Zk4x0rKLRdu4PJSIiIiIiIiaOibYiBoptcQhByp7F4PHkjo9FAC6tWWCjYiIiIiIiJo2JtiIGimNRAWbsg4VbIkS+6/d5+0Edycbg8cjIiIiIiIiakxkuQdbbVJSUnDlyhXk5OSgoKAAgiBgwoQJ5g6LyOSMeciBVAUbq9eIiIiIiIiIGkkF27Vr1/Dqq68iMDAQ/v7+GDx4MEaPHo1JkyZh8uTJos8cOnQIixYtwqJFi/DFF1+YOGKihmesQw4y80qQklkk2hbKBBsREREREclAVFQUFAqF9uvAgQPmDskiY6K6k3UFm0ajwbvvvotly5ZBrVZDEIQafRQKheizHh4e2t/MADB8+HC0adOmQeMlMhVBXQlNaYFom6EVbLXtvxbKAw6IiIiIiIiI5FvBVlFRgYceeggffvghKisra7RLJdaqdezYEYMHD9Ym5TZt2tQgcRKZg7okX7JNZWfYIQdSCTa/Fo5o7mJr0FhERERERHRvBw4c0KlsioqKMndIZhcQEKD9fgQEBJg7HKIaZJtge/bZZ/H7778DqEqmCYKA/v37Y/78+Vi8eLFoNdvdRo8erX3922+/NVisRKamLs6VbFM5uBk0ltQBB6xeIyIiIiIiIqoiywTb/v37sXHjRm1irU2bNjh27BgOHjyIqKgojB07Vq9xHnnkEQCAIAiIi4tDaWlpQ4ZNZDKa4tuSbUoDKtjyi8tx7VahaBv3XyMiIiIiIiKqIssE28KFCwFUJcb8/f0RGxuLHj16GDyOv78/XF2r9qOqqKjAuXPnjBonkblIHXCgsLaDspmN3uOcuSo+DgB0CXQ3OC4iIiIiIiJziIqKgiAI2q9BgwaZOyRqZGSXYMvJyUFsbKx27fXnn38ODw+POo/XqVMn7esLFy4YI0Qis5M8QdTesAMOzlwVX2ra3NkG3u52BsdFRERERERE1BjJLsEWExMDjUYDQRDg6emJxx57rF7j3Zmcy8jIqG94RBbBWAm2pGTxBFunALd7HiRCRERERERE1FRYmTsAQ6WnpwOoOtigLstC7+bk5KR9XVgovtcUkdyoJfZgMyTBJggCkq6JJ9iCAww7KIGIiIiIyNw0Gg2OHj2KK1euIC0tDba2tvD19cWAAQPg6elZr7EFQcDx48dx7tw5ZGRkoKysDJ6envD19UVERAQcHR2N9C7qLz09HX/99RcyMjKQnZ0NR0dHeHl5ITw8HIGBgUadSxAEHDt2DBcvXkRaWhqsra3h7e2NiIgItGrVyqhzNbTLly/j3LlzuHbtGvLz86FQKODm5gYfHx/06tWr3r+HSP5kl2DLycnRvnZzq/+H/JKSEu3rZs2a1Xs8IksgVcFmyAEHKZlFyC+qEG0Lvo8JNiIiIiKyLAcOHMDgwYO11wsWLEBUVBRKS0uxbNkyfP3117hx40aN51QqFYYNG4aPP/5YZwshfWRlZWHJkiXYuHEjMjMzRftYW1tjyJAhiIqKQq9evWodT2qVyMKFC7V7kYuJjo6udU+xiooKrFu3DqtXr8bp06cl+7Vt2xZz5szBlClTYGV173TB+vXrMXnyZO31t99+i0mTJkGj0eCzzz7DypUrkZycLPps//79sWzZslq/J5MmTcKGDRtq3L927VqtK2omTpyI9evX69yLiorS+R7e63tWVFSE3bt3Y/v27Thw4MA9V7yFhYVh9uzZePrpp6FSqWrtS42T7JaIOjs7a18XFBTUe7xbt25pX7u7c9N2ahw0UktEHfSvYEuS2H9NoQA6+Bm21JSIiIiIyBxu3ryJPn36YP78+aLJNQBQq9XYs2cPunXrhjVr1ug99s6dOxEYGIjly5dLJtcAoLy8HHv37kXv3r0xbdo0VFZWGvw+6iMhIQEdOnTAjBkzak2uAcDFixcxbdo0hIeHIzU1tU7zZWRkYMiQIXjttdckk2sAcOjQIfTv3x+bNm2q0zwNLTIyEk8//TS2bt2q13ZSx48fx7hx4zBs2LBafz9Q4yW7CrY7yy4vXrxYr7HUajVOnDihvW7ZsmW9xiOyFB4PzYZrn0ioi/KgLrn9v//mwT5A/2XVUvuvBbZ0goOt7H50EBEREVETU1ZWhkceeQSJiYkAAKVSie7duyMgIADl5eU4e/aszkF3FRUVmDFjBhQKBaZNm1br2N988w2mTp0KjUajc79Nmzbo1KkT7O3tcf36dRw7dgxqtVrb/tVXX+H69ev4+eef9aoQq6/du3djzJgxKC4u1rnfsmVLdO3aFe7u7igqKsKZM2d0Pl8nJiaiV69eOHr0qEFLOUtKSjBixAjExcUBAGxsbNCjRw/4+vpCrVbj7NmzOHPmjLZ/RUUFJk2ahM6dOyMkJKSe79a47v61dXZ2RqdOneDl5QUnJyeUlpYiJSUFp06d0lkZ98cff+Chhx5CbGwsbGxsTB02mZHsPiV36dIFQNVa7vPnzyMlJaXOa7f37Nmj/UGjUCjQu3dvo8VJZE7WHv6w9vCv1xhJV8Wr4Dpx/zUiIiKSiUq1Bpl5peYOo9HzdLWFlcryFkd9+eWXyM2t+kfjyMhILFu2DL6+vjp94uLiMGPGDCQkJGjvvfzyy4iIiEBwcLDouImJiZgxY4ZOAqZbt25YtWoV+vbtq9M3MzMT7777Lr788kvtvb1792L+/PlYsmRJjbGrK76OHj2KyMhI7f1XXnkFr776quR79fb2rnHvzJkzePrpp3WSaw899BAWLlyInj171uh/4sQJvPLKKzh06BAAIDU1FZGRkThw4IDeSx7nz5+PrKws2NnZYeHChXjhhRfg4OCg0+fYsWMYO3YsLl26BKAqyTZr1izs37+/xngff/wxoqKiAAARERHaqjpfX1/ExMRIxmGsPe+6du2KcePGYfjw4ZLLh4uLi7Fp0ya888472kq348ePY8GCBfjwww+NEgfJg+wSbB07doSvry9SU1MhCAI++eQTLF++3OBxNBqN9geaQqFA165d4erKZW9EAFBWrsbFFPGDEjpz/zUiIiKSicy8Uoye/7u5w2j0flx0P1o2tzd3GDVUJ9dee+01fPzxx6J9wsPD8eeff+Khhx7SJpbKy8sxffp07fXdnn32WZSXl2uvIyIi8Ouvv8Levub3wNPTE2vWrEFQUBBef/117f2PPvoIkZGR2gKSagEBAQCAq1ev6tx3dXXVtulDo9FgzJgxKCoq0t6LiorCggULJJ8JDQ3FH3/8gTFjxmD79u0AgJiYGHz//feYMGGCXvNmZWXBwcEB0dHRCA8PF+3Ts2dP7N+/H507d9Zu+xQdHY1Lly4hKChIp6+Hhwc8PDwAQKfiz8rKyqDvR12sXLlSrzns7e3x3HPPYdiwYYiIiNAuRf7Xv/6Ft99+W2ebK2rcLO+fGfQwduxYAFVVbCtXrsS+ffsMHuPtt9/G0aNHtddTp041WnxEcnch5TbUGkG0jRVsRERERCQX3bt3x9KlS2vtY29vj02bNulUPcXExOD48eM1+kZHR+vcd3Z2xpYtW0STa3eaM2cORowYob3WaDR1KhTR1/bt2/H3339rr5966qlak2vVrKyssGHDBnh5eWnvSSUnpSxfvlwyuVbNz89PZxmuIAiIjo42aJ6GZmgCz8/PT+fXND8/Hz///LORoyJLJssE2xtvvAFnZ2coFAqo1WqMHDkSX331lV7PZmVlYdKkSVi2bJn21BFvb29MmTKlIUMmkhWpAw7sba0Q4O1k4miIiIiIiOomKioKSuW9P/a2atUKzz//vM69u0+hBFDjRMuZM2fCx8dHr1juXi74n//8B2VlZXo9a6gVK1ZoXysUCoOWKjo6Ouokv06fPl2jok6Kr6+vzqmitbkz4QhAZ390uXr00UdhbW2tvY6NjTVjNGRqskywubu7Y8WKFRAEAQqFAqWlpZgxYwbatm2Lt956C7t27dLpf+zYMXz33XcYP3482rRpg++++w6CIEAQBKhUKnz77bc6fwiImjqpAw46+rlCpZQ+DpuIiIiIyFK4ubnhoYce0rv/M888o3MttsfX3ffGjRun9/jBwcEICwvTXpeWlurs/WYsRUVFOqu1wsPDcd999xk0xuDBg3WupZbL3u3BBx/U+/CGjh076lzrc1KnJdBoNCgoKEBaWhquXr2q85WWlgZ3d3dt37Nnz5oxUjI12e3BVm3ChAm4dOkSFi9eDIVCAUEQcPny5Rrlv4IgoE+fPjrXCoVC+8wHH3yAYcOGmTp8IosmVcEWzP3XiIiIiEgmunfvbtBJnV27doWtrS1KS6sOxjh16hTKy8u1xRi5ubm4fPmytr+rq2uNJNG99O3bV2eJaVxcXI2DEerr6NGjqKio0F4HBgbqXYFW7e4TNO9837WROghAjJub7meL27fF94A2t6KiIuzevRs//fQTEhMTcfHiRVRWVur1bPU+gNQ0yDbBBgCLFi1CmzZt8MILL6CkpES75LM6iVZNEKr2krozsWZjY4OvvvoK48ePN0vsRJYq+3YpbuaUiLYFc/81IiIiIpKJdu3aGdS/euP8c+fOAQDUajWys7PRsmVLAFUngt6pbdu2Op879dGhQwed64ao2qreZL/a5s2bsXnz5nqNmZOTo1e/u5NmtWnWrJnO9Z1JQUuxdu1avP322zV+7fVlqUlDahiyTrABwMSJEzF48GAsXboU3377LUpKqhID1Um1O1UvCR07diyioqIa/NQRIjmSql4DgOAAnrRLRERE8uHpaosfF91v7jAaPU9XW3OHIKoupze6uLjoXOfk5GgTbHdXI93dt67jG1t2drbRx6w+7fNe9NnvTi5mzZqFzz77rF5j3F0JSI2b7BNsQNVpHStXrsTSpUsRExODmJgY3LhxA9nZ2SgvL4eHhwdatGiBvn37YujQoXB1ZZKASIpUgq1lc3u4O1vmX56IiIiIxFiplGjZvPbTHanxMrS67F5j3F3EYezxjaW8vNzoY4oVsDRmP/zwQ43kWnBwMJ5++mn07NkTAQEBaNGiBWxtbWFjY6PTLyAgANeuXTNluGQhGkWCrZq9vT2GDRvGPdWI6iHpap7ofS4PJSIiIiI5qcvyvLufuXPJ452b1zfE+Mbi4eGhc71kyRK89dZbRp+nMVuwYIHO9XvvvYd33nlHr4Rofn5+Q4VFFq5RJdiIqH7UGgHnrkkccMDloUREREQkIxcuXDCof2Vlpc5hACqVSiep5unpWa/xAeD8+fM6115eXgaPcS8tWrTQua5LnE3ZxYsXdU7/HDBgAObNm6fXsyUlJcjLEy9YoMav8SyQJqJ6u5pegOIytWgbTxAlIiIiIjlJSEjQ+7RHADh58qT2BFEACAkJ0Vn+5+bmhjZt2miv8/LydBIx+oiNjdW5Dg8PF+1Xn6Wjffr00Xl+3759jWKJZ0MspxVz8eJFnesRI0bo/eyRI0caxfea6oYJNiLSktp/zUqlQNtWhm/iSkRERERkLrm5udi7d6/e/Tdt2qRzHRERUaPP3fe+//57vcc/e/YsEhIStNe2trbo3r27aN+79/UqKyvTex5PT0+EhoZqr1NTU7Fnzx69n7dUd35PDPl+GOruCjQVhc7gAAAgAElEQVRDDrNYv369kaMhOWGCjYi0pBJs7Vq5wKaZysTREBERERHVT1RUlF4nOaakpOCrr77SuTdx4sQa/SZMmKBzvXLlSty8eVOvWO7eB+3pp5+ukUirdvfBfOnp6XrNUe3FF1/UuZ4zZw4KCwsNGsPS3Pk9ycrKQkVFRYPMc/e+eOfOndPrufj4eGzevLkhQiKZYIKNiLSkEmydeMABEREREclQQkIC5s6dW2ufkpISPPPMMzoJqH79+olWlw0ZMgTdunXTXt++fRuRkZEoKSmpdY7ly5fjp59+0l4rFArMmjVLsn9gYCCsra2119HR0QYllCZMmIAOHTpor8+ePYvHH38cubnif9+XkpmZie3btxv0TEPp2LGj9nVlZSWio6MbZJ6uXbvqXG/YsOGeSdQrV67gySefbLCkH8lDo0iwRUdH491338WIESPQs2dPtG/fHoGBgQZ93bmWnqgpKiqpQHJ6gWhbZ+6/RkREREQyU12J9PHHH2Ps2LFITU2t0ScuLg79+/fHoUOHtPesra2xZs0ayXHXrVunk/w6cOAA+vfvj7/++qtG36ysLMycOROzZ8/Wuf/GG28gJCREcg5ra2v069dPe339+nU89thj+Omnn3DmzBlcvXpV5+vOveOAqgMatm3bBmdnZ+2933//HSEhIfjXv/6FggLxv/cDQE5ODrZs2YLIyEi0bt0aK1askOxrSoMHD9a5njx5MlavXo2EhARcuXJF5/uRlZVV53l8fHx0lgLn5ORg8ODBOHr0aI2+ZWVlWLt2LXr37o2rV6/C1tYWjo6OdZ6b5E3Wp4j++uuvePnll3Hp0iXtvbpuKGiqDROJLNXZ63mQ+uPDCjYiIiIikptp06Zh7969SExMxKZNm7B582aEh4fD398f5eXlOHv2bI1TPQFgxYoV6Ny5s+S4YWFhWLlyJaZPn65dfpqQkIDevXsjKCgIwcHBsLW1xY0bN3Ds2LEaBy088MADeO+99+4Z/+zZs3WqtPbu3Su5p1x0dDQGDRqkcy84OBg//vgjnnjiCdy+fRtA1VLYF154AS+99BK6dOkCPz8/ODs7o7i4GHl5ebhw4QJSUlLuGZs5PPnkk3jnnXe08aWlpWHmzJmifSdOnFiv/dA++ugjDBw4UPtrd+7cOfTp0wft27dH586dYW1tjVu3buHYsWM6lY9ffPEFFi9eLPvluFQ3sk2wLVu2DG+++SaAqqRadYKsLokynvJBBCQli5eLuzpaw9fD3sTREBERERHVj42NDX755Rc8+OCDOH36NDQaDf766y/RSjMAaNasGVasWIFp06bdc+ypU6fCzc0NkydP1kmmXLp0SacA5G5TpkzBmjVr0KxZs3vOMWLECCxevBgLFiyAWq2+Z38x999/P+Lj4xEZGYn4+HjtfbVajcTERCQmJt5zjLv3JDMXOzs77NixA6NGjRKtRjSmvn374uuvv8bzzz+vs+zz/PnzoklZlUqFTz/9FM899xwWL17coLGR5ZLlEtFff/0Vc+fO1SbGFAoFBEGAIAiws7NDy5Yt4efnp/eXv78//Pz8zPyuiMwr6Wqe6P1OAW6s8CQiIiIiWWrZsiWOHj2K+fPnw9fXV7SPSqXCQw89hMTEREyfPl3vsZ944glcvnwZr7zyCjw8PCT7NWvWDMOGDcPhw4exbt06vZJr1d555x2cOnUKb775JgYMGABvb2/Y2dnp/TwABAUF4dixY/j5559x//33Sx6scKeOHTvipZdewqFDhyxmDzYA6NGjB86cOYM1a9Zg1KhRCAoKgrOzM1Qq4x/INmnSJBw6dKhGZeCdbG1tMXr0aBw7dgwvv/yy0WMgeVEIMizfCg0NxcmTJ7WJtdatW2Pu3LkYMWIEE2XUJCQlJemUrf/9998IDg6u83iCIGDEW78it6C8RtvUEe0x+eH2dR6biIiIqK4yMzPh5eWlcy8jIwOenp5miogs2YEDB3T26VqwYAGioqK01xqNBkeOHMGVK1eQlpYGW1tb+Pr6YsCAATV+nxlKo9EgISEB586dQ2ZmJsrKyuDh4YFWrVohIiICTk5O9RrfmEpLS/HXX3/h2rVryM7ORlFRERwcHODm5oagoCB07NgRzZs3N3eYFuX69es4fPgw0tLSUFFRAS8vL/j4+KBfv34W9WtbraF+dhr7c2hjI7slopcvX9Ym1wCgV69e+O233yzyNzWRXNzMKRFNrgFAMPdfIyIiIqJGQKlUol+/fjqHBxhz7PDwcISHhxt9bGOztbXFwIEDzR2GrFSvfiOqjeyWiB45cgTAP/uu/fvf/2Zyjaie/pbYfw0AOvozwUZERERERERUG9kl2DIyMgBU7bsWGhqKtm3bmjkiIvk7c1U8wRbg7Qgne/33iCAiIiIiIiJqimSXYLtzy7igoCAzRkLUeCRJJNg6cXkoERERERER0T3JLsF258kvdT2qmIj+UVGpwYUbt0XbOjPBRkRERERERHRPskuw3XlCxY0bN8wYCVHjcDH1NsorNaJtrGAjIiIiIiIiujfZJdi6dOmCzp07QxAEJCQkIDdXenN2Irq3JIkDDmytVQj04QEiRERERERERPciuwQbALz22msAqpaIfvLJJ2aOhkjepPZf6+DnCiuVLH9EEBEREREREZmULD89T5w4EaNHj4YgCFi6dCn27Nlj7pCIZEuqgi34Pi4PJSIiIiL5GDRoEARB0H5FRUWZOyQiakJkmWADgA0bNuCxxx5DZWUlRo4cifnz5yMvL8/cYRHJSl5hGVKzikXbggNcTRwNERERERERkTxZmTuAuli0aBEAoGvXroiNjUVWVhbef/99fPrpp+jTpw86deoENzc3KJWG5Q/nz5/fEOESWaykq9JJ6eAAdxNGQkRERERERCRfskywRUVFQaFQaK8VCgUEQUBxcTH++OMP/PHHH3Ualwk2amrOSOy/5uVqC09XWxNHQ0RERERERCRPsl0iejeFQqGTdDOEIAhGjoZIHv7m/mtERERERERE9SbLCjaASTGi+tJoBJy9JpFgC2CCjYiIiIiIiEhfskywRUdHmzsEMoHi4mIcPnwYKSkpuHXrFlxdXeHr64vw8HB4e3sbda6zZ88iKSkJqampKC8vh4+PDwIDA9GrVy+D9/KTixuZRSgsqRRtY4KNiIiIiIiISH+yTLANHDjQ3CE0KYIg4Pz584iLi0NcXBzi4+Nx4sQJlJaWavsMHDgQBw4cMMp8ycnJmD9/Pnbs2IGioqIa7SqVCkOGDMFbb72FwYMH13keQRDw9ddfY9WqVTh16pRoHx8fH0yYMAHz5s2Dg4NDneeyRNdvFYreVyqA9n4uJo6GiIiIiIiISL5kmWAj09i2bRtWrVqFhIQEFBQUmGTO9evX46WXXkJhoXjyBwDUajX27duH33//Ha+++iqWLVsGlUpl0Dy3bt3CuHHj8Pvvv9faLy0tDR9++CF++OEHbN68GT169DBoHkt2I0P8e+zd3B621vzRQERERERERKQvfoomSTExMUarStPHpk2bMGXKFJ399aysrBAeHo7WrVsjMzMTCQkJyM/PB1BVgbZ8+XKUlZVh1apVes9TVFSE4cOH4/jx4zr3W7VqhZCQENja2uL8+fNISkrStl2+fBnDhg3DkSNH0L59+3q+U8uQmlmzOhAAWnk2rko9IiIiIiKi/2/vzuOiLPf/j79ZVVAQEVRcMLcjam5lLllprmmLHss0K6U8lS128vRtsb6pZR2zxWzPPGnHc7TU1DqZHjO13DWXzNQUBAIVwQWVRdb794c/5+sNMzDDDMwwvJ6PB48H1821gZ/BuT/c13UBFc07N5dChQoODlbz5s1d2ufu3bsVGxtrSq7dcccdOnr0qLZs2aIvv/xS69atU0pKil544QVT2w8//FBz5syxe6xx48aZkmt16tTRwoULlZSUpJUrV+qrr77S/v37tW3bNlMy7ezZsxo6dKhycnKc+E49R7KNBFtTEmwAAAAAADjE6xJsKSkp+umnn7RixQotWLBA//znP909pSqtZs2a6t69ux577DHNmzdPv/76q86fP68pU6a4dJxnnnlGeXl5lvKdd96pZcuWqWnTpqZ6derU0fTp0/XOO++Yrr/44ot2LWPdtGmTli5daikHBgZq3bp1Gj16dInDDLp3767NmzerZcuWlmvx8fGaPXu2Q9+bp0rhCTYAAAAAAFzCK5aIJiUladasWfrmm2+UlJRU4uv3339/iWsbN260nEYaFhamJ554osLnWdW88MILevPNN+XvX7Fhsn79ev3www+Wcv369fXxxx+XenrnxIkTtWLFCssS1vT0dM2aNUsvvfRSqWMVf/pt8uTJpe6rFh4errlz55oOU3j99df16KOPKiQkpNSxPFlufqFOnrX+JF6TyNqVPBsAAAAAAKq2Kv0EW1FRkV544QW1bt1a7733nhITE2UYhunDlvr162vq1KmaNm2a/vrXvyo+Pr4SZ141REREVHhyTVKJpwzHjx+v8PDwUtv4+PjomWeeKbWf4pKSkvTTTz9ZyrVq1dLEiRPLnF+fPn103XXXWcoZGRn65ptvymznyY6fypatlwdLRAEAAAAAcEyVTbDl5+dr8ODBmjFjhgoKCkp83cfHp9T2MTEx6tu3ryUJt3DhwgqZJ0pXWFio//znP6ZrsbGxdrUdNGiQGjVqZCnHx8dr3759NusvX77cVB42bJjCwsLsGqv4nJYtW2ZXO09la3mor4/UKDyokmcDAAAAAEDVVmUTbA8++KDWrl0r6VIyzTAM3XDDDXrppZc0ffr0Up9eu2zEiBGWz9esWVNhc4VtO3fu1OnTpy3lRo0aqU2bNna19fX11Y033mi6tmrVKpv1V69ebSr36dPH7nkWr7tmzRoVFRXZ3d7TJKdlWr3eMDxIAf5V9tcCAAAAAABuUSXvpH/44Qf961//siTWWrZsqR07dujHH3/U1KlTNWbMGLv6GTp0qCTJMAzt3LlTFy9erMhpw4r9+/ebyj179nSofa9evUzl3377rULGatu2rerVq2cpZ2VlKTEx0e72nuYYBxwAAAAAAOAyVTLBNm3aNEmXEmPR0dHasmVLqRvV2xIdHa26detKurTk9NChQy6dJ8p24MABU7lVq1YOtb/yhE9r/V12/vx5HTt2rNS2ZWnRooVdY1UFyTYSbOy/BgAAAACA46pcgu3MmTPasmWLfHx85OPjo9mzZ6t+/frl7q9du3aWzw8fPuyKKcIBcXFxpnKzZs0cal+8/pEjR+wap379+goKcmyvMXvHqgps7cHGE2wAAAAAADiu4o+IdLFNmzZZ9r6KjIzU7bff7lR/Vybn0tLSnOoLjsvIyDCVIyMjHWpfvP6FCxdUVFQkX19z7tjZcay1OXfunMN9WJOWlqb09HSH2hRPGDoiN79QJ8/mWP1ak8ja5e4XAAAAAIDqqsol2E6cOCHp0sEG5VkWWlydOnUsn2dmWt/4HRWn+M+8Vq1aDrUvXt8wDGVlZZn+XV0xjrU2Fy5ccLgPaz788EPLsufKcOJ0tmydAcITbAAAAAAAOK5KLhG9LCwszOn+cnL+70megIAAp/uDY4onvmrWrOlQe2uJMmuJUmfHsTZWVU3IJqdZXx7q6yNFhTu2bBYAAAAAAFTBBFtISIjlc1c8QXTy5EnL51eeEulOjz/+uGWPuYr8mDp1qru/1RJ8fHwqtL4z7co7lqdJTrOeGGxYL0gB/lXuVwIAAAAAF5s6darp3nHDhg3unhLg8arcEtGIiAjL585uMl9YWKg9e/ZYyo0aNXKqPziudm3znl9XPlFoD2v1i/fpinGstbE2Tnk8+uijuuuuuxxqExcXp2HDhpVrvGO2DjiIZHkoAAAAAADlUeUSbFdffbWkS3tt/f7770pJSVGTJk3K1deqVauUnZ0t6dLTST169HDZPGGfikiwBQeXTBR5coItMjKyXIculFeyjQRbU/ZfAwAAQBWWmJioq666ylIeO3as5s+f774JWTF//nzFxsZayvPmzdO4cePcNyEALlPlEmwxMTFq3Lixjh07JsMw9NZbb2nWrFkO91NUVKTXXntN0qXkWqdOnVS3bl1XT7dc7rjjjnInDR3Ru3fvCh+jLKGhoaayo6dpFj/5NSQkpMQJoq4Yx9pYnhIvjkqx9QQbCTYAAAAAAMqlyiXYJGnMmDGaOXOmDMPQ+++/ryFDhmjAgAEO9TF58mRt27bNUv7LX/7i6mmW24ABAxz+fqqq1q1bm8pJSUkOtS9ev3h/tq6np6crOztbQUH2b+pv71ieLDe/UCfPWn96r0mka57IAwAAAACguqmSO5o/88wzCgkJkY+PjwoLC3XHHXdozpw5drU9deqUxo0bpzfeeMOyaX3Dhg31wAMPVOSUYUNMTIypHBcX51D7o0ePltrfZSEhIYqKijJdi4+Pd2ishIQEu8byZCdOZ8swrH+NJ9gAAAAAACifKplgq1evnt59910ZhiEfHx9dvHhREyZMUOvWrfX888/rm2++MdXfsWOHFixYoPvuu08tW7bUggULZBiGDMOQn5+f5s2bp8DAQDd9N9Vbhw4dTOWtW7c61H7z5s2l9ueqsQ4dOqTTp09bykFBQab9HaqK5DTry0N9faSocPuf5gMAAAAAAP+nSi4RlaT7779fcXFxmj59unx8fGQYhuLj4zVz5kxTPcMw1LNnT1P58lHDhmHo73//uwYOHFjZ08f/161bN9WrV09nzpyRJJ04cUKHDx9WmzZtymxbVFSkjRs3mq7dcsstNusPHjxYa9assZQ3bNighx56yK55Fj+WetCgQVb3evN0yWmZVq83rBekAP+q9/0AAAAAAOAJqmyCTZJefvlltWzZUo8++qhycnIsSz4vJ9EuM/7/mrgrE2s1atTQnDlzdN9997ll7rjE399ft912mz7//HPLtXnz5unvf/97mW3XrFmj48ePW8otW7ZUx44dbdYfPny4Jk2aZCmvWLFCGRkZdh1WUPz0oeHDh5fZxhMds3XAQSTLQwEAAIDiCgoKtH37dsXHxys9PV2FhYWKiIhQ8+bN1atXL9WoUcPdU1R8fLwOHTqkpKQknT9/Xj4+PgoLC1NUVJS6d++uiIgId08RqBaqdIJNunT0ct++fTVz5kzNmzdPOTmXNnA3rGw0dXlJ6JgxYzR16lQ1b968kmcLa+6//35Tgm3u3Ll6+umnFR4eXmq74k8r3n///aXWb968uW644QbLU285OTmaPXu2pkyZUmq7H3/8Udu3b7eU69atq9tvv73UNp4qmRNEAQAA4GWaN29u9bC0zz//3HSfUdy8efM0btw4q19LSkrSyy+/rK+++krnzp2zWicoKEhDhw7VtGnTSt2fOTEx0eb2MrGxsYqNjbXZNiEhocR9a1ZWlr799lstW7ZMGzZsUFpams32ktS1a1dNmjRJo0aNkp+fX6l1AZRflU+wSVKzZs30/vvva+bMmdq0aZM2bdqk5ORknT59Wnl5eapfv74aNGigXr16qV+/fnY9sYTKc/PNN+vmm2/WunXrJF06iOKRRx7Rl19+aXMZ5rvvvqv169dbyvXr19dTTz1V5livvfaabrjhBlN56NChuvbaa63WP3PmjB588EHTtWeffVahoaFljuWJUmwk2JqSYAMAAAAkSR9++KEmTZqk3NzcUutlZ2dryZIlWrZsmV588UVNnTq1UuY3evRo/ec//7G7/u7du3Xvvffqs88+0xdffMETbUAF8YoE22VBQUEaOHAge6q5UGJiotXrp06dMpUvXrxos27dunXLTGq+8cYb6tmzp/Ly8iRJS5cu1YgRI/Tuu++qadOmlnoXLlzQzJkz9eqrr5rav/rqq6pTp04Z343Uu3dv3XnnnVq6dKkkKS8vT/369dMnn3yikSNHmhJ627dv19ixY02njbZs2VITJ04scxxPlJtfqJNnc6x+jSfYAACANzIKC1RwvvSne+A8/5BI+fh5x63ltGnTrCbK2rdvr9atW8vf319Hjx7Vnj17LKumCgsLNW3aNB07dkyffvpphc+xqKjIVA4JCVG7du0UGRmpOnXq6OLFi0pJSdG+ffssK7wkad26dRo8eLC2bNniEUtbAW/jHb8FUWHsPSlz+/btNutOmTKlzL/mdO3aVZ999pnuvfdey7UVK1bo22+/1XXXXaemTZvq1KlT2rlzp86fP29qO2HCBLsPK5Au7acWHx+vPXv2SJLOnz+v0aNH65lnnlGnTp0UGBiow4cPa//+/aZ2YWFhWrlypYKCquZpmydOZ8vKymlJUpPI2pU7GQAAgEpQcD5NiW8Pdfc0vF7zSSsVEBbltvE3bdqkgoICpaSkmFarjBgxQm+++abNdvXr1zeVv/vuuxL3LX379tV7772n9u3bm64nJSVp0qRJWrZsmeXa3Llz1a1btxL3Jk2aNFFCQoKkSw8S/M///I/la2+88YbuvPNOm3Ns0qSJ1eudOnXSvffeqyFDhqhdu3ZW62RnZ2vhwoV64YUXLMtId+/erSlTpmjGjBk2xwRQPiTY4DHGjBmjvLw8TZw4UZmZl067LCgo0JYtW6zW9/Hx0cSJE/XWW285NE5wcLC+++473Xvvvfrhhx8s15OTk5WcnGy1TcuWLbVo0SL96U9/cmgsT5KcZn15qK+PFBVeNZOGAAAAgK0kVO3ate3edzs3N1fjx483Xbvrrru0aNEiq/uWRUdH66uvvtITTzyh999/33L9qaee0vDhw03LMP39/S3zKJ7Uq1+/vsN7g7///vt2tQkKCtL48eM1cOBA9e7d23Kv89FHH2ny5MkKCQlxaFwApbO+wRXgJrGxsfrll180ZswYBQdbX7bo6+ur/v3764cfftA777xTro06GzZsqO+//14ff/yxrr76apv1GjVqpGeffVa//PKLunXr5vA4nsTW/msN6wUpwJ9fBQAAAKi+Fi5cqBMnTljK0dHRmjdvXpn3Gu+88446d+5sKWdnZ+ujjz6qsHlKcjgh16xZM82aNctSPn/+vEN7uAGwD0+woVTWTmOtaC1atNC//vUvZWVladOmTUpJSVFaWprq1q2rqKgoXXfddWrUqJHT4/j4+Ojhhx/Www8/rAMHDmj//v06fvy48vLyFBUVpRYtWqhHjx42D1qoalLSMq1ebxLJ/msAAACo3oqfNjp58mSbf/C/kp+fn6ZPn65bb73Vcm3+/Pl66aWXXD5HZ9x2220KDAy07Hm9ZcsWjRkzxs2zArwLCTZ4rODgYA0aNKhSxmrXrp3NvQu8RbKNJ9g44AAAAADVWX5+vnbs2GEp+/v7a9SoUXa3Hzx4sCIiIpSeni5JSkhI0PHjxxUVVbn70hUVFSkrK0sXLlywJNKuVK9ePaWmpkqSDh48WKlzA6oDEmxANWFriWhTEmwAAACoxn777TfTaZvt27d3aH8yPz8/de/eXd9++63l2s6dO3XHHXe4dJ7FZWVl6dtvv9XXX3+tvXv36siRIyooKLCr7dmzZyt0bkB1RIINqAZy8wt18myO1a/xBBsAAPBW/iGRaj5ppbun4fX8QyLdPQWnXD5h87I2bdo43Efbtm1NCbbifbra3LlzNXnyZMtTc446d+6ci2cEgAQbUA2cOJ0tW9vpNYmsXbmTAQAAqCQ+fv4KCKvcZXqoeoo/zRUaGupwH8XbnDlzxqk5leapp57SO++841QfRUVFLpoNgMtIsAHVQHKa9eWhvj5SVHhQJc8GAAAA8BzFD3bz8fFxuk9X9GHNkiVLSiTX2rdvr1GjRum6665T8+bN1aBBA9WsWVM1atQw1WvevLmSkpIqZF4ASLAB1YKt/dca1gtSgL93nJIKAAAAlEe9evVM5fIsnyzeJiwszKk52TJlyhRT+ZVXXtELL7xgV0Lv/PnzFTInAJeQYAOqgZS0TKvXm0Sy/xoAAACqt8hI8x5yhw8fdriP33//vdQ+XeHIkSOm0z9vvPFGvfjii3a1zcnJUUZGhsvnBOD/8OgKUA0k23iCjQMOAAAA4C3Kuyyzffv2qlWrlqW8f/9+h572Kiws1Pbt203XunXr5tI5SpcSbFe69dZb7W67devWEkthAbgWCTagGrC1RLQpCTYAAAB4ieJ7juXm5trVLiAgQNddd52lXFBQoMWLF9s97po1a0ynhl511VWKirJ+uEZ55yipxBNojhzGMH/+fLvrAigfEmyAl8vNL9TJszlWv8YTbAAAAPAWdevWNZVPnDhhd9v777/fVH7ttdeUnZ1dZrvCwsISyzTHjh1bIXMsvq/boUOH7Gr3888/64svvrB7HADlQ4IN8HInTmfL1tPgTSJrV+5kAAAAgApSs2ZNNW/e3FLeuXOn3fuOjRkzRg0bNrSUExISNH78eBUVFZXa7m9/+5t2795tKdeqVUsTJkywWT8mJsZUXrNmjV3zk6ROnTqZyp9//rlSU1NLbXP06FHdddddys/Pt3scAOVDgg3wcslp1peH+vpIUeFBlTwbAAAAoOL07dvX8nl2drYGDx6sL7/8Uvv371dCQoISExMtH5mZ/3cQWI0aNTRnzhxTX4sWLdKgQYNMBwtc9scff+iuu+7S7NmzTdfffvvtUg84iI6OVosWLSzlrVu3asyYMVq9erV+//130/wSExNVUFBgqRsVFaXevXtbymfOnFHfvn21bdu2EuPk5uZq7ty56tGjhxITE1WzZk3Vrs0f14GKxCmigJeztf9aw3pBCvAnxw4AAADvMXHiRC1YsMCSmNq+fbtGjRplte68efM0btw4S/m2227Tiy++qOnTp1uurV27Vu3atdPVV1+t1q1by8/PTwkJCdq1a1eJQwNiY2P1yCOPlDnHp556Sk888YSlvHDhQi1cuNBq3YSEBNNTea+//rpuuukmy/d36NAh9ezZU3/605/UoUMHBQYG6uTJk9qxY4cpgfjee+9p+vTppmsAXIu7a/hx/i4AACAASURBVMDLpaRZ/0+0MfuvAQAAwMt07txZc+bMUc2aNcvV/pVXXtG7776rwMBA0/Vff/1Vy5Yt05IlS/Tzzz+bkmu+vr564YUX9I9//MOuMR577DG7EnHW9OrVS59++qkCAgJM13///Xd99dVXWrRokdatW2dJpPn5+Wn27NkaP358ucYDYD8SbICXSzll4wTRSBJsAAAA8D6xsbE6dOiQXn75ZfXv319NmjRRcHCwfHx87Gr/xBNP6NChQxo3bpxCQkJs1qtVq5ZGjBihffv2afr06Xb37+Pjo48++kjbt2/Xk08+qZ49eyoyMtLupOC4ceO0ceNG9enTx2admjVrasSIEdqxY4cmTpxoV78AnONjFH+uFYDH++2339ShQwdLef/+/Wrfvr3Vun/+3++VeqbkKaITR7TXqJtbVtgcAQAAnJWenl5iP6u0tDRFRES4aUaobvLz87V9+3bFxcXp1KlTKiwsVEREhKKjo3X99deX+0k5V/njjz+0efNmHT9+XPn5+YqMjFRUVJSuv/561alTx61zg/tU1O9OR+5DqyP2YAO8WG5+oU6eLZlck6SmLBEFAAAAShUQEKDevXubDhfwJM2aNVOzZs3cPQ0AYoko4NVOnM6WrWdUm0RyihAAAAAAAK5Agg3wYslp1vdf8/WRosKDKnk2AAAAAAB4JxJsgBdLSbeeYGtYL0gB/rz8AQAAAABwBe6wAS+WkpZp9Xpj9l8DAAAAAMBlSLABXizllPUn2JpGkmADAAAAAMBVSLABXszWHmxNeIINAAAAAACXIcEGeKm8/EKdPJtj9WtNSbABAAAAAOAyJNgAL3X8dLYMw/rX2IMNAAAAAADXIcEGeClby0N9faSo8KBKng0AAAAAAN6LBBvgpVLSrSfYGtQLUmCAXyXPBgAAAAAA70WCDfBSKWmZVq9zwAEAAAAAAK5Fgg3wUimnrD/B1jSSBBsAAAAAAK5Egg3wUrb2YOMJNgAAAAAAXIsEG+CF8vILdfJsjtWvNSXBBgAAAACAS5FgA7zQ8dPZMgzrX2tMgg0AAAAAAJciwQZ4IVvLQ319pKjwoEqeDQAAAAAA3o0EG+CFUtKtJ9ga1AtSYIBfJc8GAAAAAADvRoIN8EIpaZlWr3PAAQAAAAAArkeCDfBCKac4QRQAAAAAgMri7+4JAHC90+dzrV5vGkmCDQAAVG2nTp1y9xQAwKPxe9I9SLABXmjB5D46m5mnY+lZSk7LVEp6tpLTM9W2WV13Tw0AAMAp7dq1c/cUAAAogQQb4IV8fHxUr04N1atTQ1e3qOfu6QAAAAAA4NXYgw0AAAAAAABwAgk2AAAAAAAAwAksEQUAAADgkcLDw5WWlubuaQBAlRceHu7uKXg9EmwAAAAAPJKvr68iIiLcPQ0AAMrEElEAAAAAAADACSTYAAAAAAAAACeQYAMAAAAAAACcQIINAAAAAAAAcAIJNgAAAAAAAMAJJNgAAAAAAAAAJ5BgAwAAAAAAAJxAgg0AAAAAAABwAgk2AAAAAAAAwAkk2AAAAAAAAAAnkGADAAAAAAAAnECCDQAAAAAAAHACCTYAAAAAAADACSTYAAAAAAAAACeQYAMAAAAAAACc4O/uCQBwXG5urqkcFxfnppkAAAAAAKqD4vedxe9LqzsSbEAVlJycbCoPGzbMTTMBAAAAAFRHycnJ6tq1q7un4TFYIgpUQRkZGe6eAgAAAACgGuO+1IwEG1AFnT9/3t1TAAAAAABUY9yXmrFEFKiCrr32WlN58eLFateunZtmA1SMuLg40/LnFStWqFWrVm6cEeB6xDmqA+Ic1QFxjurgwIEDGjlypKVc/L60uiPBBlRBISEhpnK7du3Uvn17N80GqBytWrUizuH1iHNUB8Q5qgPiHNVB8fvS6o4logAAAAAAAIATSLABAAAAAAAATiDBBgAAAAAAADiBBBsAAAAAAADgBBJsAAAAAAAAgBNIsAEAAAAAAABOIMEGAAAAAAAAOIEEGwAAAAAAAOAEEmwAAAAAAACAE0iwAQAAAAAAAE4gwQYAAAAAAAA4wd/dEwDguIiICE2ZMsVUBrwNcY7qgDhHdUCcozogzlEdEOel8zEMw3D3JAAAAAAAAICqiiWiAAAAAAAAgBNIsAEAAAAAAABOIMEGAAAAAAAAOIEEGwAAAAAAAOAEEmwAAAAAAACAE0iwAQAAAAAAAE4gwQYAAAAAAAA4gQQbAAAAAAAA4AQSbAAAAAAAAIATSLABAAAAAAAATiDBBgAAAAAAADiBBBsAAAAAAADgBH93TwCA4xISErR3714dP35cmZmZatSokaKjo9WrVy8FBAS4e3qACgsLFRcXpwMHDuj48eM6d+6catSoobCwMLVs2VLXXnutgoODXTpmfn6+Nm/erD/++EMnTpxQ7dq1FRUVpS5duqh58+YuHQtwB2IcnuTQoUP65ZdflJKSopycHNWsWVORkZFq1aqVOnXq5NTveGId7pSTk6O9e/fq4MGDOnv2rC5evKiQkBBFRkaqa9euatWqlXx8fJwehzhHVVKZ8Vql73UNAFXGkiVLjJ49exqSrH7Uq1fPmDBhgpGenu7uqaIaSkpKMmbNmmUMHTrUCAkJsRmnkgw/Pz9j8ODBxrfffuv0uGlpacaECROMevXq2RyvV69extKlS13wXQL2ufvuu0vEYXR0dLn6IsbhKTIyMoypU6ca0dHRZf6Ov+aaa4y///3vDvVPrMOdtmzZYowcOdIIDAwsNb4bN25svPTSS8bp06fLNQ5xDleIj483vvjiC+Ppp582brrpJqNOnTouec9RXGXGqzfc65JgA6qACxcuGKNGjSr1P/srPxo0aGCsXr3a3dNGNTJ69Gi747P4x6233mqkpqaWa9zvvvvOiIyMtHusMWPGGJmZmS7+7gGzr7/+2mr8lefNLjEOT7F48WIjPDzcod/vDRo0sLt/Yh3ukp+fbzz22GOGj4+Pw/G9atUqh8YizuGM9evXGwMHDiw12eXKBFtlxas33ev6GIZhCIDHKiws1O23367vvvvOdD0iIkJdunRRaGio4uPjtWfPHl35cq5Ro4bWrl2r3r17V/aUUQ1de+212rVrV4nrjRs3VuvWrdWgQQMVFBTo6NGj+uWXX1RUVGSq16ZNG/34449q2LCh3WNu2LBBgwYNUl5enuWaj4+PunbtqhYtWigjI0N79uzRqVOnTO1uu+02rVixQr6+bEMK18vIyFD79u11/PjxEl+Ljo5WYmKi3X0R4/AU06ZN09SpU0tcb9asmdq0aaOIiAhdvHhRJ06c0K+//qqsrCxJUoMGDZSamlpm/8Q63MUwDI0cOVJLly4t8bW2bdsqJiZGtWrVUnp6un7++WedPXvWVCcwMFBff/21Bg8eXOZYxDmc9c477+ipp56yq66j7zmKq6x49bp7Xbem9wCU6emnnzZl7AMCAoz33nvPyM3NNdX77bffSjxSGx4ebhw/ftxNM0d1cs0111jirkuXLsZ7771nxMXFWa2bkpJiPPTQQyX+GtW7d2+jqKjIrvGSk5ONsLAwU/vrr7/eOHDggKnexYsXjdmzZxsBAQGmus8//7zT3zNgzQMPPGCJM2eWaxDj8BRvvvlmid/Xo0ePNvbt22e1fmFhobFp0ybjqaeeMtq1a1dm/8Q63GnOnDkl4vvGG280fv311xJ18/Pzjc8++8wIDQ011Y+MjDQyMjJKHYc4hyvMmjXL6hNdNWrUMFq2bOmyJ9gqM1697V6XBBvgweLj40v8wlqxYoXN+tnZ2SV+8Tz88MOVOGNUV9dee60xdOhQY+fOnXa3+eCDD0q8QVi0aJFdba9MYkiX9n7IycmxWX/58uUl3ogkJibaPVfAHt9//70lxvz9/Uu8EXbkzS4xDk+wd+9ew9/f33Tjs2TJErvb5+fnl1mHWIc7NW/evERyLS8vr9Q2O3fuNGrUqGFqN2PGjFLbEOdwhVmzZhkBAQFG586djfHjxxuffPKJsWvXLiMvL89Yv369yxJslRWv3nivS4IN8GD333+/6RfIuHHjymzz+++/mzZn9ff3N+Lj4ythtqjOEhISytVuxIgRphgfMmRImW0OHz5s+Pn5WdoEBgYahw8fLrPd2LFjTWPFxsaWa86ANZmZmaYbtWeeeabcb3aJcXiC/Px8o2vXrqaY+uyzz1w6BrEOd9q3b1+JP/Tt3bvXrraPP/54iad7bCHO4SpnzpyxmehyVYKtMuPVG+91WcgNeKicnJwS+0E8++yzZbZr06aNhg0bZikXFBRo4cKFLp8fcKXyHs/92GOPmcrr168vs83ChQtVWFhoKf/5z39W69aty2xX/PWzePFiXbx40c6ZAqV7/vnnLXudtGjRwup+VfYixuEJlixZot27d1vK/fr1U2xsrEvHINbhTkePHjWVmzZtqk6dOtnV9o477jCVjxw5YrMucQ5XCQsLU82aNSt0jMqKV2+91yXBBnio//73v8rOzraUe/bsqbZt29rVtvgb4GXLlrl0boCrdOnSxVTOyclRRkZGqW2WL19uKtt7wxcTE6Pu3btbyllZWVqzZo2dMwVs27Jliz744ANL+ZNPPlGtWrXK3R8xDk/wySefmMqTJ092+RjEOtzp8mEclzVp0sTutk2bNjWVix9+cCXiHFVJZcWrt97rkmADPNTq1atN5T59+tjd9oYbbpC/v7+lvGfPHp08edJVUwNc5so4vezK04qKS01N1S+//GJqf/3119s9XvHX0apVq+xuC1iTm5urBx54wHIy7tixY9W/f/9y90eMwxPExcXpxx9/tJSbN2+uvn37unQMYh3uVvzkckeeDitet169elbrEeeoSiozXr31XpcEG+Ch9u/fbyr37NnT7rbBwcG6+uqrTdd+++03l8wLcKW4uDhT2d/fX/Xr17dZv/jromPHjgoODrZ7vF69epnKvC7grKlTp+r333+XdOlI+bfeesup/ohxeILiy/X79esnHx8fl45BrMPdunXrpho1aljKBw8eVE5Ojl1td+3aVaIva4hzVCWVGa/eeq9Lgg3wUAcPHjSVW7Vq5VD7li1bmsoHDhxwek6AqxXfe+Haa6+Vr6/t/5qKxzGvC7jT7t279eabb1rK77zzjsLDw53qkxiHJ9ixY4epfPnGxzAMrV27VrGxsWrXrp1CQ0MVHBys6Oho9e/fXzNmzLDsRVgWYh3uVqdOHd1///2W8sWLF/WPf/yjzHaFhYV6//33TdfGjh1rtS5xjqqkMuPVW+91SbABHujMmTM6c+aM6VqzZs0c6qN4/dI2XwXcITMzs8Qb2eHDh5fapvgTb46+LqKjo03l06dPl7pvCmBLQUGBHnjgARUUFEiSBg8erHvuucfpfolxeIKff/7ZVI6JiVFiYqL69++vAQMGaP78+Tp48KDOnz+v7Oxs/fHHH/rhhx/0/PPPq02bNnrsscdMe+tYQ6zDE8yYMcN0UNMzzzyjtWvX2qyfn5+vhx56SHv27LFcu/nmmzVixAir9YlzVCWVFa/efK9Lgg3wQMU3eQ8KCnLo8VxJioyMNJXPnTvn9LwAV3r++eeVmppqKdetW1fjx48vtU3x10bxOC9L7dq1S5y+xGsD5TFjxgzLPiXBwcH66KOPXNIvMQ5PcOLECVM5Oztb3bp107p168psm5+frw8//FC9e/cu0c+ViHV4gnr16mn9+vWWQ5dycnI0aNAg3X333VqyZIl+/fVXxcXFadu2bZo1a5auvvpqffbZZ5b21113nZYuXWpzCTVxjqqksuLVm+91S+4uDcDtMjMzTeXynEZXvM2FCxecmhPgSsuXLy+xvOLVV1+1uUnwZa56bVy5OTGvDTjqwIEDmj59uqX8yiuvmJ6AcAYxDk9Q/OYnNjZWp06dknQpofzII4/olltuUZMmTZSVlaVffvlFn332mTZt2mRps2fPHo0YMUI//vijAgICSoxBrMNTNG/eXNu3b9f8+fM1Z84c7dq1S4sXL9bixYtttgkPD9ekSZP0P//zP1bj+zLiHFVJZcWrN9/r8gQb4IGK/9Ip/pcAexT/pVO8T8BdfvnlF9OeJ5I0cOBATZgwocy2vDbgbkVFRXrwwQeVm5srSbrmmms0ceJEl/VPjMPdcnNzLfF9WUpKiiSpXbt2OnjwoN58803169dPf/rTn9S1a1fFxsZq48aNpj0JJWnr1q16/fXXrY5DrMOTFBYWqrCwUDVq1CjzQI+mTZvqzTff1KRJk0pNrknEOaqWyopXb35dkGADqoDynNzl6tO+AFf4448/NHToUNN/gtHR0frXv/5VaXHOawPOmD17trZt2ybp0qm3c+fOlZ+fX4WNR4yjshUWFlq9HhoaqtWrV6tp06Y22/7tb3/TU089Zbo2a9Ysu258iHW4y+bNmxUTE6MJEyZo8+bNKioqKrV+cnKyYmNj1axZM82dO9ehsYhzVCWe/N7cU18XJNgAD1S7dm1T2d4jw0trU7xPoLKlpaVpwIABOnbsmOVaw4YN9f333ysiIsKuPnhtwJ2OHj2qF1980VKeNGmSOnfu7NIxiHG4W1BQkNXTnCdNmlRqcu2yV155RaGhoZbymTNntGrVqhL1iHV4gh9++EH9+/c3nX7buHFjzZgxQ3v27FFGRoby8vKUmpqq1atXa+zYsfL3v7TLUnp6uv7yl7/ooYcekmEYVvsnzlGVVFa8evPrggQb4IG8+ZcOqqczZ86of//+Onz4sOVa/fr1tXbtWrVu3drufnhtwF0Mw9Bf/vIXy8mILVq00NSpU10+DjEOT2Bts+niS/tLa/vnP//ZdG3Dhg0l6hHrcLf09HSNHj3atF/UbbfdpgMHDujZZ59V586dFRoaqoCAADVo0ECDBg3S/PnztXHjRoWHh1vafPrpp5o5c6bVMYhzVCUk2JxHgg3wQFf+5Ve6dHpXVlaWQ32kpaWZynXr1nV6XkB5nDt3TgMHDtSvv/5quRYWFqbvv/9e7du3d6iv4q+N9PR0h9pnZmaW+A+Z1wbs8emnn5pOUPzkk0/KtSlvWYhxeILiMdOgQQOHDvLo0aOHqXzw4MESdYh1uNvbb79tiru2bdtq8eLFCgkJKbVdjx499OWXX5quTZs2rcR7b4k4R9VSWfHqzfe6JNgADxQeHq6wsDDTtT/++MOhPpKSkkxlR54SAlzlwoULGjx4sHbt2mW5FhISotWrV5draV3xOC4e52UpXr9evXolXmuANVOmTLF8PmTIELVq1UqJiYmlfqSmppr6KCgoKFEnLy/PVIcYhydo06aNqdyoUSOH2kdFRZnKp0+fLlGHWIe7LVmyxFR+9tln7d5svV+/frrhhhss5ZycHH3xxRcl6hHnqEoqK169+V7X390TAGBdTEyMtmzZYinHxcUpJibG7vZHjx4t0R9QmbKysjRkyBDLhvDSpce3V61apeuuu65cfRaP47i4OIfaF39dtGvXrlzzQPVz5V9kv/vuO1111VUO93Hs2LES7fbs2WNKNhPj8ATt27fXDz/8YCnXqFHDofbF61+5BO8yYh3ulJWVpfj4eNO1fv36OdRH//79tXHjRkt5+/btJeoQ56hKKjNevfVelyfYAA/VoUMHU3nr1q12t83KytK+fftK7Q+oSDk5Obr11lu1adMmy7WgoCCtXLlSvXr1Kne/xeN43759lj2x7LF58+ZS+wPcjRiHJ+jYsaOpnJGR4VD74vWv3K/qMmId7mQtphs2bOhQH8Xrnzp1qkQd4hxVSWXGq7fe65JgAzzU4MGDTWVrGwTbsnHjRhUUFFjKXbp0UYMGDVw1NaBUFy9e1O23326K2Zo1a+qbb77RjTfe6FTfjRo1Mt34FRQUmJJ4ZSn+Orrlllucmg/gasQ4PMEtt9wiHx8fS/no0aNWn0KzZf/+/aZykyZNStQh1uFO1vZrcnQPqMzMTFPZ2ibrxDmqksqMV2+91yXBBnioQYMGmTbQ3rp1qw4dOmRX2/nz55vKw4cPd+XUAJvy8vL05z//WWvXrrVcq1GjhlasWOHw0gtbisfzvHnz7Gp36NAh0/KN4OBgDRw40CVzgvfLyMiQYRgOfaxfv97UR3R0dIk61vYiJMbhblFRUerZs6elnJ+fb1oyWpbVq1ebylfuVXUlYh3uEhwcXOIwgz179jjUx5X7y0q2n4AjzlGVVFa8euu9Lgk2wEMFBQXpzjvvNF17/fXXy2x3+PBhLV++3FL29/fXPffc4/L5AcUVFBRo5MiRWrVqleVaQECAli5dqkGDBrlsnDFjxsjPz89SXrZsmY4cOVJmu+Kvn5EjR9q9mTFQmYhxeILY2FhT+e2337ar3caNG7Vjxw5L2dfXV0OGDLFal1iHO/Xp08dUnjNnjt1tU1NT9c0335iu2UokE+eoSiorXr32XtcA4LHi4+ONgIAAQ5Ll4+uvv7ZZPycnx+jVq5ep/sMPP1yJM0Z1VVBQYIwcOdIUe/7+/sayZcsqZLwHHnjANFavXr2MnJwcm/VXrFhhqh8YGGgkJiZWyNyAy9avX2+Ku+joaLvbEuNwt4KCAiMmJsYUV2+99VapbU6ePGm0bNnS1GbUqFGltiHW4S7//ve/TbHk4+NjLFiwoMx2Fy9eNPr3729qW7t2bePMmTM22xDnqGjOvOcorrLi1RvvdUmwAR7u6aefNv0SCQgIMN577z0jNzfXVO/AgQMlfuGEh4cbx48fd9PMUZ3cf//9ptiTZMycOdNISEhw+KO0/8AvS05ONsLCwkzjXX/99cbBgwdN9S5evGi8++67Jf7zfv755yvqRwFYOPNmlxiHJ1izZo3h6+triq2JEydaTSR8//33RqtWrUx1w8LCjKNHj5Y6BrEOdyksLDQ6depUIsk2ceJEm++f161bZ3Tu3LnEe55XXnml1LGIc7hKcnKy1ffPixYtMsVM48aNbb7XTk9PL3OMyopXb7vX9TEMwxAAj1VYWKjbbrvNtOxOkiIjI9W1a1fVqVNHR48e1e7du3XlyzkwMFBr1661+bg64EpXbobtrPXr15dYtmHNhg0bNGjQIOXl5Znmcc0116hFixY6d+6cdu/erfT0dFO7W2+9VStWrDA9/g5UhA0bNqhv376WcnR0tBITEx1qT4zD3d5//3098cQTpmsBAQHq0aOHGjdurJycHO3du1dJSUmmOoGBgfrmm2/s2iKAWIe7xMXF6frrr1daWprpuq+vrzp27KgWLVqoVq1aOnPmjPbs2aPU1NQSfQwZMkQrVqxQQEBAqWMR53CF5s2bl/h966ixY8eW2MesuMqKV6+713Vvfg+APS5cuGDcfffdJf5aZusjMjLSWLVqlbunjWrE3ti052P9+vV2j7ty5UojIiLC7r5Hjx5tZGZmVtwPAriCK5ZrEOPwBB9++KERFBRkdxw2aNDA2Lx5s0NjEOtwl4MHDxrXXnutw+9XfHx8jIceesjIzs62eyziHM6Kjo52+r322LFj7RqrsuLVm+51OeQAqAJq166tL774QkuWLFGPHj1s1qtXr54mTJig/fv3lzj6GPBGQ4YM0f79+/XII48oLCzMZr0ePXpo6dKlWrhwoYKDgytxhoBziHF4ggkTJmjfvn269957VadOHZv1GjZsqKlTp+r3339Xr169HBqDWIe7tG3bVlu3btXnn3+unj17lvlUfq1atTRmzBht2bJFn3zyiekkxLIQ56hKKitevelelyWiQBWUkJCg3bt36/jx48rKylLDhg0VHR2t66+/XoGBge6eHuAWeXl52rx5s5KSkpSamqrg4GA1btxYXbp00VVXXeXu6QFOI8bhCXJycrR582alpKQoNTVVgYGBioiIUKdOndSxY0eXjEGsw53OnTunn3/+WQkJCcrIyFBubq7q1KmjsLAwdejQQVdffbX8/f2dHoc4R1VSmfFale91SbABAAAAAAAATmCJKAAAAAAAAOAEEmwAAAAAAACAE0iwAQAAAAAAAE4gwQYAAAAAAAA4gQQbAAAAAAAA4AQSbAAAAAAAAIATSLABAAAAAAAATiDBBgAAAAAAADiBBBsAAAAAAADgBBJsAAAAAAAAgBNIsAEAAAAAAABOIMEGAAAAAAAAOIEEGwAAAAAAAOAEEmwAAAAAAACAE0iwAQAAAAAAAE4gwQYAAAAAAAA4gQQbAAAAAAAA4AQSbAAAAAAAAIATSLABAAAAAAAATiDBBgAAAAAAADiBBBsAAAAAAADgBBJsAAAAAAAAgBNIsAEAAAAAAABOIMEGAAAAAAAAOIEEGwAAgJfq06ePfHx8LB8AAACoGCTYAAAAgAqyZMkSU5IzKSnJ3VMCAAAVgAQbAACAB0pMTDQlZsaNG+fuKaEcVqxYYfm8S5cuio6OduNsAABARSHBBgAAAFSA/Px8rVy50lIePny4G2cDAAAqEgk2AAAAoAKsW7dO586ds5RJsAEA4L383T0BAAAAVIwNGza4ewrV2vLlyy2ft2rVSh06dHDjbAAAQEXiCTYAAADAxQzD0Ndff20p8/QaAADejQQbAAAA4GJbt25VamqqpTxs2DA3zgYAAFQ0logCAACgTAUFBdqyZYsSEhJ04sQJBQYGqm3btrrpppsUHBxcatuMjAxt2rRJR44cUU5OjsLDw9WhQwf17NlTvr7O/733zJkz2rJli1JTU3Xq1CnVrFlTERER6ty5s9q3b+90/+Vx5emhDRs2VI8ePSpt7NTUVO3YsUPHjh3TuXPn1KBBA7Vt21Y9evSQj49Pufs9e/as9u7dqyNHjujcuXPKzc1VrVq1VLduXUVHRysmJkaNGzd24XcCAEAVYgAAAMBjREdHG5Ic/pg3b16Jvm666SZTndLMmzfPan/nz583nnvuOaNBgwZWx61Tp44xZcoUIz8/v0SfiYmJxpgxY4yAgACrbaOioowFCxaU6+dUVFRkLFq0yOjevbvh6+tr8+fSuHFjY8aMGUZ2dna5CbKIeAAADhxJREFUximvVq1aWebw8MMPu7RvW/+u27ZtM4YOHWr4+fnZ/Hm/+eabVv+tSrNu3Tpj4MCBNvst/vMeP368sXfvXpd+zwAAeDqWiAIAAMCqo0eP6tprr9WMGTN08uRJq3UuXLigadOmaejQocrLy7NcX7NmjTp16qR///vfys/Pt9r2+PHjuu+++/TCCy84PK9rrrlGo0eP1vbt21VUVGSz7rFjx/Tcc8+pXbt2+u233xwap7z279+vuLg4S7ky9l/78MMP1atXL61cuVKFhYVW6xw/flxPP/20evToofT09DL7NAxDTz75pG6++WatWbPGZr9XOnbsmObOnaslS5Y4/D0AAFCVsUQUAAAAJZw9e1YDBgzQ0aNHJUkhISHq3r276tevr3Pnzmnbtm06c+aMpf6aNWv07LPPatasWdq8ebNuv/125ebmSpKioqLUpUsXhYaGKjU1VVu2bNHFixctbV977TV169bNrn3KduzYoaFDh+rUqVOm6+Hh4erSpYvq16+v3NxcHTlyRPv377d8PTExUddff702bNigzp07O/WzKcuVp4eGhoaqb9++FTrekiVL9Pjjj8swDElSZGSkunbtqrp16+rkyZPatm2bcnJyLPV37dqlvn37auPGjQoLC7PZ72uvvaZ3333XdM3f318dO3ZUdHS0goODlZOTo7Nnz+rQoUM6fvx4xXyDAABUBe5+hA4AAAD/Jzk52UhISDA2btxoWno3YsQIIyEhwebHhQsXSvTlzBLRsLAwQ5JRt25dY86cOUZeXp6pfl5enjF16lRTGz8/P2PXrl1GVFSUIclo166dsXbt2hJjnTp1yhg+fLip7VVXXWUUFRWVOscTJ04YjRo1MrXr3r27sWbNGqttjxw5YgwbNsxUv3Xr1sb58+dLHcdZXbp0sYx3zz33uLz/4v+ul/+tGjRoYHz55ZdGQUGBqX5mZqbx8ssvG4GBgaZ2Y8aMsTlGRkaGUbNmTdO/7dSpU42zZ8/abHP8+HHjs88+M/r06WO8+OKLLvt+AQCoCkiwAQAAeKCEhARTMmTs2LEO9+FMgk2SERoaavz666+ltpswYYLVZE+3bt2Mc+fO2WyXn59vdOzY0dTWWjLuSkOGDDHVf/DBB0skk6yZOHGiqd20adPKbFNeiYmJprEWL17s8jGK/7tKMho2bGgcOXKk1HbLly8vsY+arZ/5l19+aao3ZcoUh+ZY2XveAQDgbuzBBgAAAKtmz56tDh06lFrnueeeM51MefbsWQUGBuqLL75QSEiIzXb+/v7629/+Zrr2ww8/2Kz/888/67vvvrOUe/bsqTlz5sjPz6+sb0Nvv/22OnbsaCm///77luWrrnbl6aE1atTQLbfcUiHjFPfxxx+rVatWpdYZNmyYHnvsMdO12bNnW62blJRkKt95550OzadWrVoO1QcAoKojwQYAAIASmjZtqnvvvbfMes2aNVObNm1M10aOHKkWLVqU2XbgwIGm8t69e23WLb4X2GuvvSZfX/veyvr5+enJJ5+0lNPT07V161a72jrqyv3XBgwYoNq1a1fIOFe65pprdMcdd9hV96WXXlJgYKClvHLlSrsOPEhLSyv3/AAAqA5IsAEAAKCEgQMH2vV0mCS1bt3aVB48eLBd7Ro2bGhKQNk6qVSS1q5da2p300032TXGZcUPGti4caND7e1x6tQpbdq0yVKujNNDJemee+6xu254eLgpsVlUVGQ12di2bVtT+YUXXlBmZmb5JwkAgJfjFFEAAACUEBMTY3fd0NBQp9peTtycP3/eap34+HidOHHCUm7VqlWJJYxlycvLK9Gnq33zzTcqLCyUdOmpudtuu83lY1jTvXt3h+t/++23lvLOnTt1++23m+r069dPkZGRlifXtm3bptatW+vBBx/U8OHD1aVLF7ufIAQAoDogwQYAAIASwsLC7K7r729+S1m3bt1ytc3Pz7daJzk52VTetGmTrrrqKrvHsObMmTNOtbfmyuWhvXv3VkREhMvHsKb4Et2yFH/i0Nryz6CgIH388ce68847VVRUJElKTU3Vq6++qldffVV169ZVz5491bNnT91www3q0aOHatasWf5vAgCAKo4/OwEAAKAEZ55OcvWTTadPn3Zpf5J04cIFl/aXlZVlWsY6bNgwl/ZfmtIOk7Cm+BOHtpKNw4cP1+rVq60enpCRkaFVq1bppZdeUt++fRUREaExY8Zo586dDs0FAABvQYINAAAAHq348k5XMAzDpf2tWrVKFy9etJQrM8HmrCtPgS1uwIABOnjwoJYtW6ZRo0apQYMGVutlZmZq4cKFuu666/TAAw8oJyenoqYLAIBHIsEGAAAAj1a/fn1T+aGHHpJhGE59bNiwwaVzvHJ5aJcuXdS8eXOX9l8aW3vX2XLu3DlTuazlwP7+/ho+fLgWLVqk1NRUHTlyRP/85z/18MMPWz0tdt68eRozZoxDcwIAoKojwQYAAACPVvypqcOHD7tpJtbl5+dr5cqVlnJlnR56maM/jyNHjpjKkZGRDrVv1aqV7rvvPn388ceKj4/Xnj17NHr0aFOd5cuXm5bMAgDg7UiwAQAAeKDSlu1VN+3btzftG7Z161aX76HmjHXr1pmeCqvsBNu2bdscqr99+3ZTuVu3bk6N37lzZy1cuFCPPvqo6fqVT/UBAODtSLABAAB4oBo1apjKubm5bpqJ+/n5+alfv36Wcm5urhYsWODGGZldmUhq1aqVOnToUKnjL1q0yO66p0+f1po1ayxlX19f9ezZ0yXzePDBB03lhIQEl/QLAEBVQIINAADAA9WtW9dUPnHihJtm4hmeeOIJU3natGke8TMxDENff/21peyOww127dplmkNpXn75ZdOhEUOGDFFERIRL5uHv728qF08SAwDgzUiwAQAAeKCaNWuaNsrfuXOnMjIy3DchN+vTp48GDBhgKaelpWnIkCFKSUlxqJ8LFy5o4cKFLpvXtm3blJqaailX9vLQyx555BHFx8eXWmfFihX64IMPTNeefPJJq3X//e9/6+DBgw7N4Z///KepHBMT41B7AACqMhJsAAAAHqpv376Wz7OzszV48GB9+eWX2r9/vxISEpSYmGj5yMzMdONMK8fnn3+uJk2aWMp79+5Vx44dNXPmTJ06dcpmuwsXLuibb77Rgw8+qMaNG2vy5Mkum9OVy0MbNmyoHj16uKxve4WFhSk1NVW9e/fW4sWLVVRUZPp6VlaWpk+frrvvvluFhYWW6/fcc4/69+9vtc8lS5aoffv26tu3rz744AMlJibaHD89PV1PP/203n77bcs1X19f3Xvvvc59YwAAVCH+ZVcBAACAO0ycOFELFixQQUGBpEub048aNcpq3Xnz5mncuHGVOLvK16hRI3377bcaMmSIjh8/Lkk6e/asnn32WT333HOKiYlRixYtFBoaqtzcXGVkZCg+Pl6JiYkyDMPST7169Vw2pysTbHfccYd8fSv/79cfffSRRo0apdTUVN19991q0KCBrrnmGoWGhurkyZPatm2bsrOzTW3at2+v9957r9R+DcPQhg0btGHDBj3++OMKDw9X+/btFR4eruDgYGVnZ+vo0aP69ddfTYk7SZo8ebLatWvn8u8VAABPRYINAADAQ3Xu3Flz5szRo48+qosXL7p7Oh6hU6dO2rNnj+677z7TZv2GYejAgQM6cOBAmX2EhYW5ZC779+9XXFycpeyu5aF333230tLS9Ne//lVFRUU6efKkvvvuO5v1u3btqtWrVzucaDx9+rR++umnUuv4+/vrf//3f/XSSy851DcAAFUdS0QBAAA8WGxsrA4dOqSXX35Z/fv3V5MmTRQcHCwfHx93T81tIiMj9d///lc//fSThg0bpuDg4DLbXHXVVRo/frxWrVqln3/+2SXzuPLptdDQUNOS3sr2xBNP6KefftLAgQPl5+dntU5UVJTeeOMNbd++vcyDDWbPnq13331XQ4YMsSshGRISorFjx2rfvn0k1wAA1ZKPceXz8gAAAEAVU1BQoJ9//llxcXE6ffq0Lly4oKCgIIWGhqpFixaKiYlRw4YNXT5u165dtWfPHknS6NGjXXp4Qmn69OmjH3/80VIu/nb+xIkT2r59u44dO6YLFy4oMjJSbdu2VY8ePcq1hNUwDMXFxenIkSP6448/dO7cOeXn56t27dqWZaMdOnRQYGCg098bAABVFQk2AAAAwEFJSUmmU14XL16su+66q1LGLivBBgAAKh9LRAEAAAAHrVixwvJ5jRo1dMstt7hxNgAAwN045AAAAABwUKNGjTRlyhRJUpMmTVS7dm03zwgAALgTCTYAAADAQSNHjnT3FAAAgAdhiSgAAAAAAADgBBJsAAAAAAAAgBNIsAEAAAAAAABOIMEGAAAAAAAAOMHHMAzD3ZMAAAAAAAAAqiqeYAMAAAAAAACcQIINAAAAAAAAcAIJNgAAAAAAAMAJJNgAAAAAAAAAJ5BgAwAAAAAAAJxAgg0AAAAAAABwAgk2AAAAAAAAwAkk2AAAAAAAAAAnkGADAAAAAAAAnECCDQAAAAAAAHACCTYAAAAAAADACSTYAAAAAAAAACeQYAMAAAAAAACcQIINAAAAAAAAcAIJNgAAAAAAAMAJJNgAAAAAAAAAJ5BgAwAAAAAAAJxAgg0AAAAAAABwAgk2AAAAAAAAwAkk2AAAAAAAAAAnkGADAAAAAAAAnECCDQAAAAAAAHACCTYAAAAAAADACSTYAAAAAAAAACeQYAMAAAAAAACcQIINAAAAAAAAcAIJNgAAAAAAAMAJJNgAAAAAAAAAJ5BgAwAAAAAAAJzw/wDj53hN5ubGrwAAAABJRU5ErkJggg==\n",
"text/plain": [
"