From ed6448c1e3c1d972c1d49ebdd873ea98436390f1 Mon Sep 17 00:00:00 2001
From: J08nY
Date: Sun, 20 Apr 2025 21:40:08 +0200
Subject: Rename countermeasures pyecsca notebook to simulation.
---
analysis/countermeasures/countermeasures.ipynb | 1045 -------------
analysis/countermeasures/simulation.ipynb | 1910 ++++++++++++++++++++++++
2 files changed, 1910 insertions(+), 1045 deletions(-)
delete mode 100644 analysis/countermeasures/countermeasures.ipynb
create mode 100644 analysis/countermeasures/simulation.ipynb
diff --git a/analysis/countermeasures/countermeasures.ipynb b/analysis/countermeasures/countermeasures.ipynb
deleted file mode 100644
index 58c654d..0000000
--- a/analysis/countermeasures/countermeasures.ipynb
+++ /dev/null
@@ -1,1045 +0,0 @@
-{
- "cells": [
- {
- "cell_type": "markdown",
- "id": "bafc2f4e-05a3-4120-bcd6-5d1f5fb91cd9",
- "metadata": {},
- "source": [
- "# Distinguishing countermeasures by output"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 69,
- "id": "33ee6084-2ac3-4f95-9610-0fbc06026538",
- "metadata": {},
- "outputs": [],
- "source": [
- "import io\n",
- "import random\n",
- "import itertools\n",
- "import warnings\n",
- "\n",
- "import cypari2\n",
- "from cysignals.alarm import alarm, AlarmInterrupt\n",
- "\n",
- "from matplotlib import pyplot as plt\n",
- "from collections import Counter\n",
- "from tqdm.auto import tqdm, trange\n",
- "\n",
- "from pyecsca.misc.utils import TaskExecutor\n",
- "from pyecsca.ec.mod import mod, RandomModAction\n",
- "from pyecsca.ec.point import Point\n",
- "from pyecsca.ec.model import ShortWeierstrassModel\n",
- "from pyecsca.ec.params import load_params_ectester\n",
- "from pyecsca.ec.mult import LTRMultiplier\n",
- "from pyecsca.ec.context import local, DefaultContext\n",
- "from pyecsca.ec.countermeasures import GroupScalarRandomization, AdditiveSplitting, MultiplicativeSplitting, EuclideanSplitting, BrumleyTuveri\n",
- "\n",
- "%matplotlib ipympl"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "b1b9596c-1eba-4ace-af84-8cb279d84cc2",
- "metadata": {},
- "outputs": [],
- "source": [
- "model = ShortWeierstrassModel()\n",
- "coords = model.coordinates[\"projective\"]"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "b0afb195-8390-44c5-931e-75a70ccd4e9e",
- "metadata": {},
- "outputs": [],
- "source": [
- "add = coords.formulas[\"add-2015-rcb\"]\n",
- "dbl = coords.formulas[\"dbl-2015-rcb\"]\n",
- "mult = LTRMultiplier(add, dbl, complete=False)"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "52c877e1-5021-4ec2-9daa-dd20bec6bcb2",
- "metadata": {},
- "outputs": [],
- "source": [
- "gsr = GroupScalarRandomization(mult)\n",
- "asplit = AdditiveSplitting(mult)\n",
- "msplit = MultiplicativeSplitting(mult)\n",
- "esplit = EuclideanSplitting(mult)\n",
- "bt = BrumleyTuveri(mult)"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "27626337-dcbc-497c-a54e-02d50e2b8f34",
- "metadata": {},
- "source": [
- "## 3n test"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "c3088419-161b-4193-a1b6-6f623f217fcd",
- "metadata": {},
- "outputs": [],
- "source": [
- "key3n = 0x20959f2b437de1e522baf6d814911938157390d3ea5118660b852ab0d5387006\n",
- "params3n = load_params_ectester(io.BytesIO(b\"0xc381bb0394f34b5ed061c9107b66974f4d0a8ec89b9fe73b98b6d1368c7d974d,0x5ca6c5ee0a10097af291a8f125303fb1a3e35e8100411902245d691e0e5cb497,0x385a5a8bb8af94721f6fd10b562606d9b9df931f7fd966e96859bb9bd7c05836,0x4616af1898b92cac0f902a9daee24bbae63571cead270467c6a7886ced421f5e,0x34e896bdb1337e0ae5960fa3389fb59c2c8d6c7dbfd9aac33a844f8f98e433ef,0x412b3e5686fbc3ca4575edb0292232702ae721a7d4a230cc170a5561aa70e00f,0x01\"), \"projective\")\n",
- "bits3n = params3n.full_order.bit_length()\n",
- "point3n = Point(X=mod(0x4a48addb2e471767b7cd0f6f1d4c27fe46f4a828fc20f950bd1f72c939b36a84, params3n.curve.prime),\n",
- " Y=mod(0x13384d38c353f862832c0f067e46a3e510bb6803c20745dfb31929f4a18d890d, params3n.curve.prime),\n",
- " Z=mod(1, params3n.curve.prime), model=coords)"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "a8dde7e6-cd48-4f99-9677-23a19e4c2e5b",
- "metadata": {},
- "outputs": [],
- "source": [
- "print(f\"prime:\\t0x{params3n.curve.prime:x}\")\n",
- "print(f\"a:\\t0x{params3n.curve.parameters['a']:x}\")\n",
- "print(f\"b:\\t0x{params3n.curve.parameters['b']:x}\")\n",
- "print(f\"G:\\t[0x{params3n.generator.X:x},\\n\\t 0x{params3n.generator.Y:x}]\")\n",
- "print(f\"n:\\t0x{params3n.order:x}\")\n",
- "print(f\"3n:\\t0x{3 * params3n.order:x}\")\n",
- "print(f\"\\nP:\\t[0x{point3n.X:x},\\n\\t 0x{point3n.Y:x}]\")"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "cd6f8500-7509-45b0-8b23-471ee5014f42",
- "metadata": {},
- "outputs": [],
- "source": [
- "def generate_scalars_mod3(rem, samples):\n",
- " scalars = []\n",
- " while True:\n",
- " scalar = random.randint(0, params3n.full_order)\n",
- " if scalar % 3 == rem:\n",
- " scalars.append(scalar)\n",
- " if len(scalars) == samples:\n",
- " break\n",
- " return scalars\n",
- "\n",
- "def test_3n(countermeasure, scalars):\n",
- " ctr = Counter()\n",
- " for k in tqdm(scalars, leave=False):\n",
- " mult.init(params3n, point3n)\n",
- " kP = mult.multiply(k).to_affine()\n",
- " mult.init(params3n, point3n)\n",
- " knP = mult.multiply(k + params3n.full_order).to_affine()\n",
- " mult.init(params3n, point3n)\n",
- " k2nP = mult.multiply(k + 2 * params3n.full_order).to_affine()\n",
- "\n",
- " countermeasure.init(params3n, point3n)\n",
- " res = countermeasure.multiply(k)\n",
- " aff = res.to_affine()\n",
- " if aff.equals(kP):\n",
- " ctr[\"k\"] += 1\n",
- " elif aff.equals(knP):\n",
- " ctr[\"k + 1n\"] += 1\n",
- " elif aff.equals(k2nP):\n",
- " ctr[\"k + 2n\"] += 1\n",
- " else:\n",
- " ctr[aff] += 1\n",
- " for name, count in sorted(ctr.items()):\n",
- " print(f\"{name}:\\t{count}\")\n",
- "\n",
- "def test_3n_fixed_scalar(countermeasure, samples):\n",
- " test_3n(countermeasure, [key3n for _ in range(samples)])\n",
- "\n",
- "def test_3n_random_scalar(countermeasure, samples):\n",
- " test_3n(countermeasure, [random.randint(0, params3n.full_order) for _ in range(samples)])\n",
- "\n",
- "def test_3n_random_scalar_projected(countermeasure, samples):\n",
- " print(\"k = 0 mod 3\")\n",
- " test_3n(countermeasure, generate_scalars_mod3(0, samples))\n",
- " print()\n",
- " print(\"k = 1 mod 3\")\n",
- " test_3n(countermeasure, generate_scalars_mod3(1, samples))\n",
- " print()\n",
- " print(\"k = 2 mod 3\")\n",
- " test_3n(countermeasure, generate_scalars_mod3(2, samples))"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "46b8f74a-433d-48c9-b5b9-6bb7d2731246",
- "metadata": {},
- "source": [
- "### Fixed scalar experiments"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "fc82d4b9-91cd-423c-83aa-89721efa1ae9",
- "metadata": {},
- "source": [
- "#### Group scalar randomization"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "86532d50-2db7-4370-b449-c545b330a852",
- "metadata": {},
- "outputs": [],
- "source": [
- "test_3n_fixed_scalar(gsr, 1000)"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "aba3e713-6246-435d-93af-2e8b42ee9582",
- "metadata": {},
- "source": [
- "#### Additive splitting"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "ad421630-606f-4666-9bbf-1a446eec1b59",
- "metadata": {},
- "outputs": [],
- "source": [
- "test_3n_fixed_scalar(asplit, 1000)"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "a1284a11-4ace-437d-9826-2035cce36756",
- "metadata": {},
- "source": [
- "#### Multiplicative splitting"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "3ed5d7f3-0ba1-4b62-9635-aeb492499175",
- "metadata": {},
- "outputs": [],
- "source": [
- "test_3n_fixed_scalar(msplit, 1000)"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "5e40f8e9-d26f-41e7-adbf-a0fbdb680677",
- "metadata": {},
- "source": [
- "#### Euclidean splitting"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "314447c6-a1fb-4d3a-8988-b34c8912dd5e",
- "metadata": {},
- "outputs": [],
- "source": [
- "test_3n_fixed_scalar(esplit, 1000)"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "ff7185a2-3cd9-44d7-b1a9-982eaed561dc",
- "metadata": {},
- "source": [
- "#### Brumley and Tuveri bit-length fixing"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "f41dfc1d-1017-4aa0-bcd4-6569c53bf81e",
- "metadata": {},
- "outputs": [],
- "source": [
- "test_3n_fixed_scalar(bt, 1000)"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "28915553-b5e6-4108-bc23-e5844b6a63b8",
- "metadata": {},
- "source": [
- "### Random scalar experiments"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "566ddd10-2d0e-4b32-9b27-60770ab68155",
- "metadata": {},
- "source": [
- "#### Group scalar randomization"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "7255321a-6ad6-4938-8ec9-dd8d977686db",
- "metadata": {},
- "outputs": [],
- "source": [
- "test_3n_random_scalar(gsr, 1000)"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "376c0da8-b92d-4151-ac30-839ab5c0ceae",
- "metadata": {},
- "source": [
- "#### Additive splitting"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "b0146a9a-0803-43c4-ab29-8ba6e15934b5",
- "metadata": {},
- "outputs": [],
- "source": [
- "test_3n_random_scalar(asplit, 1000)"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "96823352-cd63-4cf8-8ab7-cf6e025837b9",
- "metadata": {},
- "source": [
- "#### Multiplicative splitting"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "5645ae6f-f5f4-419d-ba47-248532dc2114",
- "metadata": {},
- "outputs": [],
- "source": [
- "test_3n_random_scalar(msplit, 1000)"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "1784d763-932a-4f66-a1da-d7ef51cbc88a",
- "metadata": {},
- "source": [
- "#### Euclidean splitting"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "c9fc4f35-1c25-4cac-bb63-8bd70263db47",
- "metadata": {},
- "outputs": [],
- "source": [
- "test_3n_random_scalar(esplit, 1000)"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "8917dfd5-1548-414b-846a-685857bfb427",
- "metadata": {},
- "source": [
- "#### Brumley and Tuveri bit-length fixing"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "4fd6b288-08a9-4dbe-9145-e96401805315",
- "metadata": {},
- "outputs": [],
- "source": [
- "test_3n_random_scalar(bt, 1000)"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "ee9e23f6-edd9-4fc2-9e7c-b5c176991071",
- "metadata": {},
- "source": [
- "### Random scalar experiments projected to scalar divisor classes mod 3"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "e17fcb12-1e02-4bcf-a2b7-785c16c03028",
- "metadata": {},
- "source": [
- "#### Group scalar randomization"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "6c46fdbb-2ffb-4169-8e00-6d93b8407ee5",
- "metadata": {},
- "outputs": [],
- "source": [
- "test_3n_random_scalar_projected(gsr, 1000)"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "b9d635c8-f788-4700-876b-ac48e89557a7",
- "metadata": {},
- "source": [
- "#### Additive splitting"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "344a4f90-3470-40e9-a75f-b925a88c2480",
- "metadata": {},
- "outputs": [],
- "source": [
- "test_3n_random_scalar_projected(asplit, 1000)"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "78aaf83f-a689-40ad-ae88-58ab20e5e6f9",
- "metadata": {},
- "source": [
- "#### Multiplicative splitting"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "616a7726-01e6-4e9c-b7f2-fe8f14b60071",
- "metadata": {},
- "outputs": [],
- "source": [
- "test_3n_random_scalar_projected(msplit, 1000)"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "018c8f31-d7c3-497b-bdaf-580c6465753f",
- "metadata": {},
- "source": [
- "#### Euclidean splitting"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "adced4e4-37a7-43ed-97b5-01cb5d274d6b",
- "metadata": {},
- "outputs": [],
- "source": [
- "test_3n_random_scalar_projected(esplit, 1000)"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "8d7b0a42-7be7-464d-932d-8386ad912034",
- "metadata": {},
- "source": [
- "#### Brumley and Tuveri bit-length fixing"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "fe8d8295-3e69-4b60-b8c3-5710deaeb0b3",
- "metadata": {},
- "outputs": [],
- "source": [
- "test_3n_random_scalar_projected(bt, 1000)"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "43b309af-5683-4384-9623-d7633723177c",
- "metadata": {},
- "source": [
- "## Mask recovery\n",
- "Using a composite order curve we can recover the size and the actual mask values (in a known key scenario) in both GSR and multiplicative splitting. However, real-world targets do not like composite order curves and may either check the order or otherwise fail to compute on such curves. Thus, we lie to them and set the order to the next-prime of the true order, in this case $n + 92$."
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "20a26f27-620d-4d7f-92bd-b949482b5c9a",
- "metadata": {},
- "outputs": [],
- "source": [
- "pari = cypari2.Pari(256_000_000, 2_000_000_000)"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "144340bd-5372-4beb-a46e-fd60c596b254",
- "metadata": {},
- "outputs": [],
- "source": [
- "real_n = 0xa9fa3419aca88bade2cba14e317816c6828910c6ce04fcd2a2e857d25df50775\n",
- "# = 2898786277 * 2916913393 * 3067509271 * 3248233993 * 3894099889 * 4099407227 * 4101666977 * 13936975277\n",
- "real_n_facts = pari.factor(real_n)\n",
- "params92pn = load_params_ectester(io.BytesIO(b\"0xa9fa3419aca88bade2cba14e317816c79d52481d463dc9bcb12c37f45aa3b4e1,0x2ea3bfe6659f8e035735349b91fbfa2baf0cf8e640315f0fe03c1136813dec99,0x2b07c518e04b02158651e3dbbef7720015dd496bf15af02f8439f8e1503b8370,0x90fb04b1af19e8e20396ac052f260a9fb5f736b97e3cd4af08fe81a1e75dac6d,0x2302bcf700d3d5899f04d0c7441f5017c9758bfafd6ce15dbe36fb4eea76baec,0xa9fa3419aca88bade2cba14e317816c6828910c6ce04fcd2a2e857d25df507d1,0x01\"), \"projective\")\n",
- "e = pari.ellinit([int(params92pn.curve.parameters[\"a\"]), int(params92pn.curve.parameters[\"b\"])], int(params92pn.curve.prime))\n",
- "e[15][0] = real_n"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "f103129c-17d3-4217-999b-94ecb4ec523d",
- "metadata": {},
- "outputs": [],
- "source": [
- "print(f\"prime:\\t0x{params92pn.curve.prime:x}\")\n",
- "print(f\"a:\\t0x{params92pn.curve.parameters['a']:x}\")\n",
- "print(f\"b:\\t0x{params92pn.curve.parameters['b']:x}\")\n",
- "print(f\"G:\\t[0x{params92pn.generator.X:x},\\n\\t 0x{params92pn.generator.Y:x}]\")\n",
- "print(f\"n+92:\\t0x{params92pn.order:x} (fake order, given to the target, prime)\")\n",
- "print(f\"n:\\t0x{real_n:x} (real order, composite)\")"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "322d2e68-5259-4ea6-9748-2b0aa21b557f",
- "metadata": {},
- "source": [
- "### Group scalar randomization\n",
- "In GSR getting the mask out this way is quite simple. The target believes it is operating on a curve of order $n+92$ so it will use that value multiplied with the mask to randomize the scalar. Thus as a result we get:\n",
- "$$ P = [k + r(n + 92)]G $$\n",
- "\n",
- "However, the curve is truly of order $n$, thus arithmetic on its group will make this equal to: \n",
- "$$ P = [k + r 92]G $$\n",
- "\n",
- "Since this is a composite order curve, we can solve the dlog and obtain $k + r 92$ and since we assume we know $k$ we can easily compute both the mask size and mask value $r$."
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "08d99bd5-2b87-4a04-995d-7a87f9b67102",
- "metadata": {},
- "outputs": [],
- "source": [
- "key = 0x20959f2b437de1e522baf6d814911938157390d3ea5118660b852ab0d5387006 # any key works ofc\n",
- "gsr.init(params92pn, params92pn.generator)\n",
- "res = gsr.multiply(key)"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "2a869bed-8e21-46af-8f70-065f4afd6a82",
- "metadata": {},
- "outputs": [],
- "source": [
- "affine_gen = params92pn.generator.to_affine()\n",
- "affine_res = res.to_affine()"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "e440399a-bc01-488b-8822-08cc0bf1672d",
- "metadata": {},
- "outputs": [],
- "source": [
- "dlog = pari.elllog(e,\n",
- " [int(affine_res.x), int(affine_res.y)],\n",
- " [int(affine_gen.x), int(affine_gen.y)],\n",
- " real_n)"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "id": "7ea6d6ae-a6f5-4b53-8c40-787d79970cb6",
- "metadata": {},
- "outputs": [],
- "source": [
- "mask = int((dlog - key) / 92)\n",
- "mask_len = mask.bit_length()\n",
- "print(mask)\n",
- "print(mask_len)"
- ]
- },
- {
- "cell_type": "markdown",
- "id": "d40ec035-0656-4eda-8ef4-c14f9d53f49f",
- "metadata": {},
- "source": [
- "### Multiplicative splitting\n",
- "In multiplicative splitting the situation is a bit more complicated. Doing the same computation, where the target thinks the curve order is $n+92$ leads to:\n",
- "$$ P = [k r^{-1}\\pmod{n+92}][r \\mod n]G $$\n",
- "\n",
- "Since the curve is composite order we can easily compute the dlog $d$ of P to G, we get:\n",
- "$$ d = (k r^{-1})\\pmod{n+92}\\: r = k + t (n + 92) $$\n",
- "\n",
- "However, the dlog is computed $\\mod n$ so we really get: $ d = k + t 92$. We extract the $t$ out of this.\n",
- "Note that $t$ will have roughly the same size as the mask $r$, since at the left side we have $(k r^{-1}) {\\mod (n+92)}$\n",
- "of size $n$ and $r$ of size of the mask and on the right size we have $t n$ that dominates.\n",
- "Thus at this point we have recovered the mask size.\n",
- "However, $t$ is always smaller than $r$, sometimes also in bitsize.\n",
- "\n",
- "Now that we have $s$ we can go back to the original equation and get:\n",
- "$$ (k r^{-1})\\pmod{n+92}\\: r = k + t (n + 92) $$\n",
- "\n",
- "We can then factor this value, lets call it $full$, and look for divisors that are larger than $t$ but smaller than the mask length\n",
- "that we recovered before. There may be multiple candidates here and we don't know how to distinguish between\n",
- "them. It holds for all of the candidates $c$ that the rest of the value $full$ is equal to the inverse of $c \\mod (n+92)$.\n",
- "However, sometimes there is only one candidate, which is equal to the true mask value $r$."
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 80,
- "id": "b5f398fc-90d7-455e-97bd-62b682d55961",
- "metadata": {},
- "outputs": [],
- "source": [
- "def divisors(primes, powers):\n",
- " for comb in itertools.product(*[range(power+1) for power in powers]):\n",
- " value = 1\n",
- " for prime, power in zip(primes, comb):\n",
- " value *= prime**power\n",
- " yield value\n",
- "\n",
- "def pari_factor(number, flag=1+8, duration=5*60):\n",
- " # Use the 1 + 8 option to factorint to skip some costly factorization methods, we don't need the huge factors.\n",
- " # This means that some factors may remain, which in turn means that sometimes we may not find the true mask.\n",
- " pari = cypari2.Pari(256_000_000, 2_000_000_000)\n",
- " with warnings.catch_warnings(category=RuntimeWarning):\n",
- " warnings.simplefilter(\"ignore\")\n",
- " try:\n",
- " alarm(duration)\n",
- " factors = pari.factorint(number, flag)\n",
- " except AlarmInterrupt:\n",
- " return None\n",
- " primes = list(map(int, factors[0]))\n",
- " powers = list(map(int, factors[1]))\n",
- " return primes, powers\n",
- "\n",
- "def pari_dlog(params, P, G, real_n, facts_str):\n",
- " pari = cypari2.Pari(256_000_000, 2_000_000_000)\n",
- " e = pari.ellinit([int(params.curve.parameters[\"a\"]), int(params.curve.parameters[\"b\"])], int(params.curve.prime))\n",
- " e[15][0] = real_n\n",
- " facts = pari(facts_str)\n",
- " dlog = pari.elllog(e, P, G, facts)\n",
- " return int(dlog)"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 81,
- "id": "5f03e586-33df-4525-a722-f5f63d6ca28d",
- "metadata": {},
- "outputs": [
- {
- "data": {
- "application/vnd.jupyter.widget-view+json": {
- "model_id": "c3b98a7bfa2e4e568f83ee3d035266bb",
- "version_major": 2,
- "version_minor": 0
- },
- "text/plain": [
- "Collecting scalarmults: 0%| | 0/1000 [00:00= t:\n",
- " candidates.add(divisor)\n",
- " candidate_amounts.append(len(candidates))\n",
- " if len(candidates) == 1:\n",
- " candidate = candidates.pop()\n",
- " print(\"Only one candidate, we got the mask:\", candidate, candidate == r, \"incomplete factoring this may not be the r\" if not complete else \"\")\n",
- "print(f\"Total recovered masks: {len(list(filter(lambda a: a == 1, candidate_amounts)))} out of {i}\")"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 89,
- "id": "6274ff91-325f-4c6b-a4d7-d66b994d730f",
- "metadata": {},
- "outputs": [
- {
- "data": {
- "application/vnd.jupyter.widget-view+json": {
- "model_id": "3a03974d3d6c4f159824b17e4bfb0467",
- "version_major": 2,
- "version_minor": 0
- },
- "image/png": "iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjEsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvc2/+5QAAAAlwSFlzAAAPYQAAD2EBqD+naQAAQIRJREFUeJzt3X98z/X+//H72377sfkx9sOPIb+Z0TCbvkdlmY5icmr5dJB0SvkxVgrF+qFGRYTDoeKUI3IqSTU/dtCJMTY7UX4lx8R+oGyZYm3P7x9dvE9vhuH9ts3rdr1c3pfa6/18PR/P13te3nfP1y+bMcYIAAAAllGlvAcAAACA64sACAAAYDEEQAAAAIshAAIAAFgMARAAAMBiCIAAAAAWQwAEAACwGAIgAACAxRAAAQAALIYACAAAYDEEQAAAAIshAAIAAFgMARAAAMBiCIAAAAAWQwAEAACwGAIgAACAxRAAAQAALIYACAAAYDEEQAAAAIshAAIAAFgMARAAAMBiCIAAAAAWQwAEAACwGAIgAACAxRAAAQAALIYACAAAYDEEQAAAAIshAAIAAFgMARAAAMBiCIAAAAAWQwAEAACwGAIgAACAxRAAAQAALIYACAAAYDEEQAAAAIshAAIAAFgMARAAAMBiCIAAAAAWQwAEAACwGAIgAACAxRAAAQAALIYACAAAYDEEQAAAAIshAAIAAFgMARAAAMBiCIAAAAAWQwAEAACwGAIgAACAxRAAAQAALIYACAAAYDEEQAAAAIshAAIAAFgMARAAAMBiCIAAAAAWQwAEAACwGAIgAACAxbiX9wAqs5KSEh09elQ1atSQzWYr7+EAAIAyMMbop59+UnBwsKpUseZcGAHwGhw9elQNGzYs72EAAICrcPjwYTVo0KC8h1EuCIDXoEaNGpJ++wPk6+tbzqMBAABlUVBQoIYNG9q/x62IAHgNzh329fX1JQACAFDJWPn0LWse+AYAALAwAiAAAIDFEAABAAAshgAIAABgMQRAAAAAiyEAAgAAWAwBEAAAwGIIgAAAABZDAAQAALAYAiAAAIDFEAABAAAshgAIAABgMQRAAAAAiyEAAgAAWIx7eQ/AahqP+9TlNf47pbfLawAAgMqLGUAAAACLIQACAABYDAEQAADAYgiAAAAAFkMABAAAsBgCIAAAgMUQAAEAACyGAAgAAGAxBEAAAACLIQACAABYDAEQAADAYgiAAAAAFkMABAAAsBgCIAAAgMUQAAEAACyGAAgAAGAxBEAAAACLIQACAABYDAEQAADAYgiAAAAAFkMABAAAsBgCIAAAgMUQAAEAACyGAAgAAGAxBEAAAACLIQACAABYDAEQAADAYgiAAAAAFkMABAAAsBgCIAAAgMUQAAEAACyGAAgAAGAxBEAAAACLIQACAABYDAEQAADAYgiAAAAAFkMABAAAsBgCIAAAgMVUmgA4Z84cNW7cWN7e3oqIiFBaWtpF23799dfq37+/GjduLJvNphkzZlxznwAAADeKShEAly1bpoSEBCUmJiojI0NhYWGKiYlRXl5eqe1Pnz6tpk2basqUKQoMDHRKnwAAADeKShEAp0+frr/85S8aMmSI2rRpo3nz5qlq1ap6++23S23fuXNnvfrqq7r//vvl5eXllD4BAABuFBU+AJ49e1bp6emKjo62L6tSpYqio6OVmppaYfoEAACoLNzLewCXc/z4cRUXFysgIMBheUBAgPbs2XNd+zxz5ozOnDlj/7mgoOCq6gMAAJSnCj8DWJEkJSXJz8/P/mrYsGF5DwkAAOCKVfgA6O/vLzc3N+Xm5josz83NvegFHq7qc/z48crPz7e/Dh8+fFX1AQAAylOFD4Cenp4KDw9XSkqKfVlJSYlSUlIUGRl5Xfv08vKSr6+vwwsAAKCyqfDnAEpSQkKCBg8erE6dOqlLly6aMWOGCgsLNWTIEEnSoEGDVL9+fSUlJUn67SKPb775xv7/R44cUWZmpqpXr65mzZqVqU8AAIAbVaUIgHFxcTp27JgmTZqknJwcdejQQcnJyfaLOLKyslSlyv8mM48ePaqOHTvaf37ttdf02muvqXv37tqwYUOZ+gQAALhR2YwxprwHUVkVFBTIz89P+fn5ZT4c3Hjcpy4elfTfKb1dXgMAgMrqar6/bzQV/hxAAAAAOBcBEAAAwGIIgAAAABZDAAQAALAYAiAAAIDFEAABAAAshgAIAABgMQRAAAAAiyEAAgAAWEyleBQcrgxPGwEAAJfCDCAAAIDFEAABAAAshgAIAABgMQRAAAAAiyEAAgAAWAwBEAAAwGIIgAAAABZDAAQAALAYAiAAAIDFEAABAAAshgAIAABgMQRAAAAAiyEAAgAAWAwBEAAAwGIIgAAAABZDAAQAALAYAiAAAIDFEAABAAAshgAIAABgMQRAAAAAiyEAAgAAWAwBEAAAwGIIgAAAABZDAAQAALAYAiAAAIDFEAABAAAshgAIAABgMQRAAAAAiyEAAgAAWAwBEAAAwGIIgAAAABZDAAQAALAYAiAAAIDFEAABAAAshgAIAABgMQRAAAAAiyEAAgAAWAwBEAAAwGIIgAAAABZDAAQAALAYAiAAAIDFVJoAOGfOHDVu3Fje3t6KiIhQWlraJdsvX75crVq1kre3t0JDQ/XZZ585vH/q1CmNGDFCDRo0kI+Pj9q0aaN58+a5chMAAAAqhEoRAJctW6aEhAQlJiYqIyNDYWFhiomJUV5eXqntN2/erAEDBmjo0KHasWOHYmNjFRsbq127dtnbJCQkKDk5WYsXL9bu3bs1evRojRgxQitXrrxemwUAAFAubMYYU96DuJyIiAh17txZs2fPliSVlJSoYcOGGjlypMaNG3dB+7i4OBUWFmrVqlX2ZV27dlWHDh3ss3zt2rVTXFycJk6caG8THh6uO++8U5MnTy7TuAoKCuTn56f8/Hz5+vqWaZ3G4z4tU7uK7r9Tepf3EAAAuCpX8/19o6nwM4Bnz55Venq6oqOj7cuqVKmi6OhopaamlrpOamqqQ3tJiomJcWgfFRWllStX6siRIzLGaP369dq3b5969ux50bGcOXNGBQUFDi8AAIDKpsIHwOPHj6u4uFgBAQEOywMCApSTk1PqOjk5OZdtP2vWLLVp00YNGjSQp6enevXqpTlz5ugPf/jDRceSlJQkPz8/+6thw4bXsGUAAADlo8IHQFeZNWuWtmzZopUrVyo9PV3Tpk3T8OHDtW7duouuM378eOXn59tfhw8fvo4jBgAAcA738h7A5fj7+8vNzU25ubkOy3NzcxUYGFjqOoGBgZds//PPP2vChAn66KOP1Lv3b+eytW/fXpmZmXrttdcuOHx8jpeXl7y8vK51kwAAAMpVhQ+Anp6eCg8PV0pKimJjYyX9dhFISkqKRowYUeo6kZGRSklJ0ejRo+3L1q5dq8jISElSUVGRioqKVKWK4wSom5ubSkpKXLIdN5rrcTELF5oAAOAaFT4ASr/dsmXw4MHq1KmTunTpohkzZqiwsFBDhgyRJA0aNEj169dXUlKSJCk+Pl7du3fXtGnT1Lt3by1dulTbt2/X/PnzJUm+vr7q3r27xo4dKx8fH4WEhGjjxo165513NH369HLbTgAAgOuhUgTAuLg4HTt2TJMmTVJOTo46dOig5ORk+4UeWVlZDrN5UVFRWrJkiZ599llNmDBBzZs314oVK9SuXTt7m6VLl2r8+PF64IEH9MMPPygkJEQvvfSShg0bdt23DwAA4HqqFPcBrKisfB/A64FDwAAAV+A+gBa+ChgAAMCqCIAAAAAWQwAEAACwGAIgAACAxRAAAQAALIYACAAAYDEEQAAAAIshAAIAAFgMARAAAMBiCIAAAAAWQwAEAACwGAIgAACAxRAAAQAALIYACAAAYDEEQAAAAIshAAIAAFgMARAAAMBiCIAAAAAWQwAEAACwGAIgAACAxRAAAQAALIYACAAAYDEEQAAAAItxd3WB/fv3a/369crLy1NJSYnDe5MmTXJ1eQAAAJzHpQFwwYIFeuyxx+Tv76/AwEDZbDb7ezabjQAIAABQDlwaACdPnqyXXnpJTz/9tCvLAAAA4Aq49BzAH3/8Uffee68rSwAAAOAKuTQA3nvvvVqzZo0rSwAAAOAKufQQcLNmzTRx4kRt2bJFoaGh8vDwcHh/1KhRriwPAACAUtiMMcZVnTdp0uTihW02fffdd64qfV0UFBTIz89P+fn58vX1LdM6jcd96uJR3Tj+O6V3eQ8BAHADuprv7xuNS2cADx486MrucYO7HmGZkAkAsKLrdiNoY4xcONkIAACAMnJ5AHznnXcUGhoqHx8f+fj4qH379nr33XddXRYAAAAX4dJDwNOnT9fEiRM1YsQIdevWTZL05ZdfatiwYTp+/LjGjBnjyvIAAAAohUsD4KxZszR37lwNGjTIvqxPnz5q27atnnvuOQIgAABAOXDpIeDs7GxFRUVdsDwqKkrZ2dmuLA0AAICLcGkAbNasmd5///0Lli9btkzNmzd3ZWkAAABchEsPAT///POKi4vTF198YT8HcNOmTUpJSSk1GAIAAMD1XDoD2L9/f23dulX+/v5asWKFVqxYIX9/f6Wlpalfv36uLA0AAICLcOkMoCSFh4dr8eLFri4DAACAMnJ6ACwoKLA/VqWgoOCSba36+BUAAIDy5PQAWKtWLWVnZ6tevXqqWbOmbDbbBW2MMbLZbCouLnZ2eQAAAFyG0wPgv/71L9WuXVuStH79emd3DwAAgGvk9ADYvXt3+/83adJEDRs2vGAW0Bijw4cPO7s0AAAAysClVwE3adJEx44du2D5Dz/8oCZNmriyNAAAAC7CpQHw3Ll+5zt16pS8vb1dWRoAAAAX4ZLbwCQkJEiSbDabJk6cqKpVq9rfKy4u1tatW9WhQwdXlAYAAMBluCQA7tixQ9JvM4A7d+6Up6en/T1PT0+FhYXpySefdEVpAAAAXIZLAuC5q3+HDBmiN954QzVq1HBFGQAAAFwFl50DWFRUpHfffVeHDh1yVQkAAABcBZcFQA8PDzVq1IibPQMAAFQwLr0K+JlnntGECRP0ww8/XHNfc+bMUePGjeXt7a2IiAilpaVdsv3y5cvVqlUreXt7KzQ0VJ999tkFbXbv3q0+ffrIz89P1apVU+fOnZWVlXXNYwUAAKjIXBoAZ8+erS+++ELBwcFq2bKlbr75ZodXWS1btkwJCQlKTExURkaGwsLCFBMTo7y8vFLbb968WQMGDNDQoUO1Y8cOxcbGKjY2Vrt27bK3OXDggG655Ra1atVKGzZs0FdffaWJEydyexoAAHDDsxljjKs6f/755y/5fmJiYpn6iYiIUOfOnTV79mxJUklJiRo2bKiRI0dq3LhxF7SPi4tTYWGhVq1aZV/WtWtXdejQQfPmzZMk3X///fLw8NC7775b1s25QEFBgfz8/JSfny9fX98yrdN43KdXXQ/O998pvct7CACA6+xqvr9vNC65Cvicsga8Szl79qzS09M1fvx4+7IqVaooOjpaqamppa6TmppqvxfhOTExMVqxYoWk3wLkp59+qqeeekoxMTHasWOHmjRpovHjxys2NvaiYzlz5ozOnDlj/7mgoODqNwwAAKCcuDQAnpOenq7du3dLktq2bauOHTuWed3jx4+ruLhYAQEBDssDAgK0Z8+eUtfJyckptX1OTo4kKS8vT6dOndKUKVM0efJkTZ06VcnJybrnnnu0fv16h+cZ/15SUtJlZzVRubh6RpYZRgBAReTSAJiXl6f7779fGzZsUM2aNSVJJ0+e1G233aalS5eqbt26rix/USUlJZKkvn37asyYMZKkDh06aPPmzZo3b95FA+D48eMdZhYLCgrUsGFD1w8YAADAiVx6EcjIkSP1008/6euvv9YPP/ygH374Qbt27VJBQYFGjRpVpj78/f3l5uam3Nxch+W5ubkKDAwsdZ3AwMBLtvf395e7u7vatGnj0KZ169aXvArYy8tLvr6+Di8AAIDKxqUBMDk5WX/961/VunVr+7I2bdpozpw5+vzzz8vUh6enp8LDw5WSkmJfVlJSopSUFEVGRpa6TmRkpEN7SVq7dq29vaenpzp37qy9e/c6tNm3b59CQkLKNC4AAIDKyqWHgEtKSuTh4XHBcg8PD/th2LJISEjQ4MGD1alTJ3Xp0kUzZsxQYWGhhgwZIkkaNGiQ6tevr6SkJElSfHy8unfvrmnTpql3795aunSptm/frvnz59v7HDt2rOLi4vSHP/xBt912m5KTk/XJJ59ow4YN17bRAAAAFZxLZwBvv/12xcfH6+jRo/ZlR44c0ZgxY9SjR48y9xMXF6fXXntNkyZNUocOHZSZmank5GT7hR5ZWVnKzs62t4+KitKSJUs0f/58hYWF6Z///KdWrFihdu3a2dv069dP8+bN0yuvvKLQ0FC9+eab+uCDD3TLLbc4YcsBAAAqLpfeB/Dw4cPq06ePvv76a/vFEocPH1a7du20cuVKNWjQwFWlrwvuA4jL4SpgAKh4uA+giw8BN2zYUBkZGVq3bp39li2tW7dWdHS0K8sCAADgElx+H0CbzaY77rhDd9xxh6tLAQAAoAxceg6gJKWkpOiuu+7STTfdpJtuukl33XWX1q1b5+qyAAAAuAiXBsC//vWv6tWrl2rUqKH4+HjFx8fL19dXf/zjHzVnzhxXlgYAAMBFuPQQ8Msvv6zXX39dI0aMsC8bNWqUunXrppdfflnDhw93ZXkAAACUwqUzgCdPnlSvXr0uWN6zZ0/l5+e7sjQAAAAuwqUBsE+fPvroo48uWP7xxx/rrrvucmVpAAAAXIRLDwG3adNGL730kjZs2GB/DNuWLVu0adMmPfHEE3rjjTfsbcv6bGAAAABcG5feCLpJkyZlG4TNpu+++85Vw3AZbgSNy+FG0ABQ8XAjaBfPAB48eNCV3QMAAOAquPw+gOcYY+TCyUYAAACUkcsD4DvvvKPQ0FD5+PjIx8dH7du317vvvuvqsgAAALgIlx4Cnj59uiZOnKgRI0aoW7dukqQvv/xSw4YN0/HjxzVmzBhXlgcAAEApXBoAZ82apblz52rQoEH2ZX369FHbtm313HPPEQABAADKgUsPAWdnZysqKuqC5VFRUcrOznZlaQAAAFyESwNgs2bN9P7771+wfNmyZWrevLkrSwMAAOAiXHoI+Pnnn1dcXJy++OIL+zmAmzZtUkpKSqnBEAAAAK7n0hnA/v37Ky0tTf7+/lqxYoVWrFghf39/paWlqV+/fq4sDQAAgItw2QxgUVGRHn30UU2cOFGLFy92VRkAAABcIZfNAHp4eOiDDz5wVfcAAAC4Si49BBwbG6sVK1a4sgQAAACukEsvAmnevLleeOEFbdq0SeHh4apWrZrD+6NGjXJleQAAAJTCpQHwrbfeUs2aNZWenq709HSH92w2GwEQAACgHLg0AB48eND+/8YYSb8FPwAAAJQfl54DKP02C9iuXTt5e3vL29tb7dq105tvvunqsgAAALgIl84ATpo0SdOnT9fIkSMVGRkpSUpNTdWYMWOUlZWlF154wZXlAQAAUAqXBsC5c+dqwYIFGjBggH1Znz591L59e40cOZIACAAAUA5cegi4qKhInTp1umB5eHi4fv31V1eWBgAAwEW4NAAOHDhQc+fOvWD5/Pnz9cADD7iyNAAAAC7CpYeApd8uAlmzZo26du0qSdq6dauysrI0aNAgJSQk2NtNnz7d1UMBAACAXBwAd+3apZtvvlmSdODAAUmSv7+//P39tWvXLns7bg2DG1XjcZ+6vMZ/p/R2eQ0AwI3FpQFw/fr1ruweAAAAV8Hl9wEEAABAxUIABAAAsBgCIAAAgMUQAAEAACyGAAgAAGAxBEAAAACLIQACAABYDAEQAADAYgiAAAAAFkMABAAAsBgCIAAAgMUQAAEAACyGAAgAAGAxBEAAAACLIQACAABYDAEQAADAYgiAAAAAFkMABAAAsJhKFQDnzJmjxo0by9vbWxEREUpLS7tk++XLl6tVq1by9vZWaGioPvvss4u2HTZsmGw2m2bMmOHkUQMAAFQslSYALlu2TAkJCUpMTFRGRobCwsIUExOjvLy8Uttv3rxZAwYM0NChQ7Vjxw7FxsYqNjZWu3btuqDtRx99pC1btig4ONjVmwEAAFDuKk0AnD59uv7yl79oyJAhatOmjebNm6eqVavq7bffLrX9zJkz1atXL40dO1atW7fWiy++qJtvvlmzZ892aHfkyBGNHDlS//jHP+Th4XE9NgUAAKBcVYoAePbsWaWnpys6Otq+rEqVKoqOjlZqamqp66Smpjq0l6SYmBiH9iUlJRo4cKDGjh2rtm3bumbwAAAAFYx7eQ+gLI4fP67i4mIFBAQ4LA8ICNCePXtKXScnJ6fU9jk5Ofafp06dKnd3d40aNapM4zhz5ozOnDlj/7mgoKCsmwAAAFBhVIoZQFdIT0/XzJkztWjRItlstjKtk5SUJD8/P/urYcOGLh4lAACA81WKAOjv7y83Nzfl5uY6LM/NzVVgYGCp6wQGBl6y/b///W/l5eWpUaNGcnd3l7u7uw4dOqQnnnhCjRs3LrXP8ePHKz8/3/46fPjwtW8cAADAdVYpAqCnp6fCw8OVkpJiX1ZSUqKUlBRFRkaWuk5kZKRDe0lau3atvf3AgQP11VdfKTMz0/4KDg7W2LFjtXr16lL79PLykq+vr8MLAACgsqkU5wBKUkJCggYPHqxOnTqpS5cumjFjhgoLCzVkyBBJ0qBBg1S/fn0lJSVJkuLj49W9e3dNmzZNvXv31tKlS7V9+3bNnz9fklSnTh3VqVPHoYaHh4cCAwPVsmXL67txAAAA11GlCYBxcXE6duyYJk2apJycHHXo0EHJycn2Cz2ysrJUpcr/JjSjoqK0ZMkSPfvss5owYYKaN2+uFStWqF27duW1CQAAABWCzRhjynsQlVVBQYH8/PyUn59f5sPBjcd96uJRwWr+O6V3eQ8BACqVq/n+vtFUinMAAQAA4DwEQAAAAIshAAIAAFhMpbkIBEDprsd5pZxnCAA3FmYAAQAALIYACAAAYDEcAgZwWRxmBoAbCzOAAAAAFkMABAAAsBgCIAAAgMUQAAEAACyGAAgAAGAxBEAAAACLIQACAABYDAEQAADAYgiAAAAAFkMABAAAsBgCIAAAgMUQAAEAACyGAAgAAGAxBEAAAACLIQACAABYDAEQAADAYgiAAAAAFkMABAAAsBgCIAAAgMUQAAEAACyGAAgAAGAxBEAAAACLIQACAABYDAEQAADAYgiAAAAAFkMABAAAsBgCIAAAgMW4l/cAAOB6aTzuU5f2/98pvV3aPwA4CzOAAAAAFkMABAAAsBgOAQOoEFx9eBYA8D/MAAIAAFgMARAAAMBiCIAAAAAWQwAEAACwGAIgAACAxRAAAQAALIYACAAAYDEEQAAAAIvhRtAAAAfX46bcPDcZKF/MAAIAAFgMARAAAMBiCIAAAAAWQwAEAACwmEoVAOfMmaPGjRvL29tbERERSktLu2T75cuXq1WrVvL29lZoaKg+++wz+3tFRUV6+umnFRoaqmrVqik4OFiDBg3S0aNHXb0ZAAAA5arSBMBly5YpISFBiYmJysjIUFhYmGJiYpSXl1dq+82bN2vAgAEaOnSoduzYodjYWMXGxmrXrl2SpNOnTysjI0MTJ05URkaGPvzwQ+3du1d9+vS5npsFAABw3dmMMaa8B1EWERER6ty5s2bPni1JKikpUcOGDTVy5EiNGzfugvZxcXEqLCzUqlWr7Mu6du2qDh06aN68eaXW2LZtm7p06aJDhw6pUaNGlx1TQUGB/Pz8lJ+fL19f3zJtx/W4vQKA8nGj3NqE28DgRnc13983mkoxA3j27Fmlp6crOjravqxKlSqKjo5WampqqeukpqY6tJekmJiYi7aXpPz8fNlsNtWsWdMp4wYAAKiIKsWNoI8fP67i4mIFBAQ4LA8ICNCePXtKXScnJ6fU9jk5OaW2/+WXX/T0009rwIABF/3XwJkzZ3TmzBn7zwUFBVeyGQAAABVCpZgBdLWioiLdd999MsZo7ty5F22XlJQkPz8/+6thw4bXcZQAAADOUSkCoL+/v9zc3JSbm+uwPDc3V4GBgaWuExgYWKb258LfoUOHtHbt2kueCzB+/Hjl5+fbX4cPH77KLQIAACg/lSIAenp6Kjw8XCkpKfZlJSUlSklJUWRkZKnrREZGOrSXpLVr1zq0Pxf+9u/fr3Xr1qlOnTqXHIeXl5d8fX0dXgAAAJVNpTgHUJISEhI0ePBgderUSV26dNGMGTNUWFioIUOGSJIGDRqk+vXrKykpSZIUHx+v7t27a9q0aerdu7eWLl2q7du3a/78+ZJ+C39/+tOflJGRoVWrVqm4uNh+fmDt2rXl6elZPhsKAADgYpUmAMbFxenYsWOaNGmScnJy1KFDByUnJ9sv9MjKylKVKv+b0IyKitKSJUv07LPPasKECWrevLlWrFihdu3aSZKOHDmilStXSpI6dOjgUGv9+vW69dZbr8t2AQAAXG+V5j6AFRH3AQTwezfKve24DyBudNwHsBLNAAJARUdwAlBZVIqLQAAAAOA8BEAAAACLIQACAABYDAEQAADAYgiAAAAAFkMABAAAsBgCIAAAgMUQAAEAACyGG0EDQCXC04QAOAMzgAAAABZDAAQAALAYDgEDAG5Irj5cznOZUZkxAwgAAGAxBEAAAACL4RAwAAAWdj2uLOdwecXDDCAAAIDFEAABAAAshgAIAABgMQRAAAAAiyEAAgAAWAxXAQMArrsb4ZnGXD1bsVzJ76PkzGkXjqRyYAYQAADAYgiAAAAAFkMABAAAsBgCIAAAgMUQAAEAACyGAAgAAGAxBEAAAACLIQACAABYDDeCBgCggroRbpiNiokZQAAAAIshAAIAAFgMARAAAMBiCIAAAAAWQwAEAACwGAIgAACAxRAAAQAALIYACAAAYDEEQAAAAIvhSSAAAMCleKJJxcMMIAAAgMUQAAEAACyGAAgAAGAxBEAAAACLIQACAABYDAEQAADAYgiAAAAAFkMABAAAsBgCIAAAgMUQAAEAACymUgXAOXPmqHHjxvL29lZERITS0tIu2X758uVq1aqVvL29FRoaqs8++8zhfWOMJk2apKCgIPn4+Cg6Olr79+935SYAAACUu0oTAJctW6aEhAQlJiYqIyNDYWFhiomJUV5eXqntN2/erAEDBmjo0KHasWOHYmNjFRsbq127dtnbvPLKK3rjjTc0b948bd26VdWqVVNMTIx++eWX67VZAAAA153NGGPKexBlERERoc6dO2v27NmSpJKSEjVs2FAjR47UuHHjLmgfFxenwsJCrVq1yr6sa9eu6tChg+bNmydjjIKDg/XEE0/oySeflCTl5+crICBAixYt0v3333/ZMRUUFMjPz0/5+fny9fUt03bwQGwAAMpXyZnTOjzjviv6/r7RuJf3AMri7NmzSk9P1/jx4+3LqlSpoujoaKWmppa6TmpqqhISEhyWxcTEaMWKFZKkgwcPKicnR9HR0fb3/fz8FBERodTU1FID4JkzZ3TmzBn7z/n5+ZJ+C4JlVXLmdJnbAgAA5zv3XVxJ5sBcolIEwOPHj6u4uFgBAQEOywMCArRnz55S18nJySm1fU5Ojv39c8su1uZ8SUlJev755y9Y3rBhw7JtCAAAqDBOnDghPz+/8h5GuagUAbCiGD9+vMOs4smTJxUSEqKsrCyX/QEqKChQw4YNdfjwYZdNU7u6xo2wDdSoOP1To2LVuBG2gRoVp//rVSM/P1+NGjVS7dq1XdJ/ZVApAqC/v7/c3NyUm5vrsDw3N1eBgYGlrhMYGHjJ9uf+m5ubq6CgIIc2HTp0KLVPLy8veXl5XbDcz8/P5ecQ+Pr6VvoaN8I2UKPi9E+NilXjRtgGalSc/q9XjSpVKs21sE5XKbbc09NT4eHhSklJsS8rKSlRSkqKIiMjS10nMjLSob0krV271t6+SZMmCgwMdGhTUFCgrVu3XrRPAACAG0GlmAGUpISEBA0ePFidOnVSly5dNGPGDBUWFmrIkCGSpEGDBql+/fpKSkqSJMXHx6t79+6aNm2aevfuraVLl2r79u2aP3++JMlms2n06NGaPHmymjdvriZNmmjixIkKDg5WbGxseW0mAACAy1WaABgXF6djx45p0qRJysnJUYcOHZScnGy/iCMrK8thKjcqKkpLlizRs88+qwkTJqh58+ZasWKF2rVrZ2/z1FNPqbCwUI888ohOnjypW265RcnJyfL29i7TmLy8vJSYmFjqYWFnuRFq3AjbQI2K0z81KlaNG2EbqFFx+r+RalR0leY+gAAAAHCOSnEOIAAAAJyHAAgAAGAxBEAAAACLIQACAABYDAHwKs2ZM0eNGzeWt7e3IiIilJaW5tT+v/jiC919990KDg6WzWazP8PYWZKSktS5c2fVqFFD9erVU2xsrPbu3evUGnPnzlX79u3tN/OMjIzU559/7tQa55syZYr9Fj/O8txzz8lmszm8WrVq5bT+JenIkSP685//rDp16sjHx0ehoaHavn270/pv3LjxBdtgs9k0fPhwp9UoLi7WxIkT1aRJE/n4+Oimm27Siy++6PRnbf70008aPXq0QkJC5OPjo6ioKG3btu2q+7vcvmaM0aRJkxQUFCQfHx9FR0dr//79Tuv/ww8/VM+ePVWnTh3ZbDZlZmY6dRuKior09NNPKzQ0VNWqVVNwcLAGDRqko0ePOq2G9Nt+0qpVK1WrVk21atVSdHS0tm7d6tQavzds2DDZbDbNmDHDqTUefPDBC/aTXr16OXUbdu/erT59+sjPz0/VqlVT586dlZWV5bQape3rNptNr776qtNqnDp1SiNGjFCDBg3k4+OjNm3aaN68eWXuvyw1cnNz9eCDDyo4OFhVq1ZVr169rmjfK8v33C+//KLhw4erTp06ql69uvr373/BQyRuVATAq7Bs2TIlJCQoMTFRGRkZCgsLU0xMjPLy8pxWo7CwUGFhYZozZ47T+vy9jRs3avjw4dqyZYvWrl2roqIi9ezZU4WFhU6r0aBBA02ZMkXp6enavn27br/9dvXt21dff/2102r83rZt2/S3v/1N7du3d3rfbdu2VXZ2tv315ZdfOq3vH3/8Ud26dZOHh4c+//xzffPNN5o2bZpq1arltBrbtm1zGP/atWslSffee6/TakydOlVz587V7NmztXv3bk2dOlWvvPKKZs2a5bQakvTwww9r7dq1evfdd7Vz50717NlT0dHROnLkyFX1d7l97ZVXXtEbb7yhefPmaevWrapWrZpiYmL0yy+/OKX/wsJC3XLLLZo6depVjf9yNU6fPq2MjAxNnDhRGRkZ+vDDD7V371716dPHaTUkqUWLFpo9e7Z27typL7/8Uo0bN1bPnj117Ngxp9U456OPPtKWLVsUHBx8RdtQ1hq9evVy2F/ee+89p/V/4MAB3XLLLWrVqpU2bNigr776ShMnTizz7cfKUuP3Y8/Oztbbb78tm82m/v37O61GQkKCkpOTtXjxYu3evVujR4/WiBEjtHLlSqfUMMYoNjZW3333nT7++GPt2LFDISEhio6OLvP3VFm+58aMGaNPPvlEy5cv18aNG3X06FHdc889Zd6GSs3ginXp0sUMHz7c/nNxcbEJDg42SUlJLqknyXz00Ucu6fucvLw8I8ls3LjRpXVq1apl3nzzTaf3+9NPP5nmzZubtWvXmu7du5v4+Hin9Z2YmGjCwsKc1t/5nn76aXPLLbe4rP/SxMfHm5tuusmUlJQ4rc/evXubhx56yGHZPffcYx544AGn1Th9+rRxc3Mzq1atclh+8803m2eeeeaa+z9/XyspKTGBgYHm1VdftS87efKk8fLyMu+999419/97Bw8eNJLMjh07rrjfstY4Jy0tzUgyhw4dclmN/Px8I8msW7fOqTW+//57U79+fbNr1y4TEhJiXn/99avq/2I1Bg8ebPr27XvVfV6u/7i4OPPnP//ZKf1frMb5+vbta26//Xan1mjbtq154YUXHJZdy354fo29e/caSWbXrl32ZcXFxaZu3bpmwYIFV1Xj/O+5kydPGg8PD7N8+XJ7m927dxtJJjU19apqVCbMAF6hs2fPKj09XdHR0fZlVapUUXR0tFJTU8txZNcmPz9fklz2YOzi4mItXbpUhYWFLnnU3vDhw9W7d2+H34sz7d+/X8HBwWratKkeeOCBKzpcczkrV65Up06ddO+996pevXrq2LGjFixY4LT+z3f27FktXrxYDz30kGw2m9P6jYqKUkpKivbt2ydJ+s9//qMvv/xSd955p9Nq/PrrryouLr5gtsTHx8eps7LnHDx4UDk5OQ5/rvz8/BQREVHp93ebzaaaNWu6pP+zZ89q/vz58vPzU1hYmNP6LSkp0cCBAzV27Fi1bdvWaf2eb8OGDapXr55atmypxx57TCdOnHBKvyUlJfr000/VokULxcTEqF69eoqIiHD6KT6/l5ubq08//VRDhw51ar9RUVFauXKljhw5ImOM1q9fr3379qlnz55O6f/MmTOS5LCvV6lSRV5eXle9r5//PZeenq6ioiKH/btVq1Zq1KhRpd6/y4oAeIWOHz+u4uJi+xNIzgkICFBOTk45jeralJSUaPTo0erWrZvDk1KcYefOnapevbq8vLw0bNgwffTRR2rTpo1TayxdulQZGRn2xwA6W0REhBYtWqTk5GTNnTtXBw8e1P/7f/9PP/30k1P6/+677zR37lw1b95cq1ev1mOPPaZRo0bp73//u1P6P9+KFSt08uRJPfjgg07td9y4cbr//vvVqlUreXh4qGPHjho9erQeeOABp9WoUaOGIiMj9eKLL+ro0aMqLi7W4sWLlZqaquzsbKfVOefcPn0j7e+//PKLnn76aQ0YMEC+vr5O7XvVqlWqXr26vL299frrr2vt2rXy9/d3Wv9Tp06Vu7u7Ro0a5bQ+z9erVy+98847SklJ0dSpU7Vx40bdeeedKi4uvua+8/LydOrUKU2ZMkW9evXSmjVr1K9fP91zzz3auHGjE0Z/ob///e+qUaOG0w9rzpo1S23atFGDBg3k6empXr16ac6cOfrDH/7glP7PBbHx48frxx9/1NmzZzV16lR9//33V7Wvl/Y9l5OTI09Pzwv+IVSZ9+8rUWkeBQfXGT58uHbt2uWSGZSWLVsqMzNT+fn5+uc//6nBgwdr48aNTguBhw8fVnx8vNauXXtF59Bcid/PYLVv314REREKCQnR+++/75R/VZeUlKhTp056+eWXJUkdO3bUrl27NG/ePA0ePPia+z/fW2+9pTvvvPOqzp+6lPfff1//+Mc/tGTJErVt21aZmZkaPXq0goODnbod7777rh566CHVr19fbm5uuvnmmzVgwAClp6c7rcaNqqioSPfdd5+MMZo7d67T+7/tttuUmZmp48ePa8GCBbrvvvu0detW1atX75r7Tk9P18yZM5WRkeHUmevz3X///fb/Dw0NVfv27XXTTTdpw4YN6tGjxzX1XVJSIknq27evxowZI0nq0KGDNm/erHnz5ql79+7X1H9p3n77bT3wwANO//tx1qxZ2rJli1auXKmQkBB98cUXGj58uIKDg51yJMbDw0Mffvihhg4dqtq1a8vNzU3R0dG68847r+rCMld+z1VWzABeIX9/f7m5uV1wlVBubq4CAwPLaVRXb8SIEVq1apXWr1+vBg0aOL1/T09PNWvWTOHh4UpKSlJYWJhmzpzptP7T09OVl5enm2++We7u7nJ3d9fGjRv1xhtvyN3d3Sn/aj9fzZo11aJFC3377bdO6S8oKOiCQNy6dWunHmY+59ChQ1q3bp0efvhhp/c9duxY+yxgaGioBg4cqDFjxjh9Zvamm27Sxo0bderUKR0+fFhpaWkqKipS06ZNnVpHkn2fvhH293Ph79ChQ1q7dq3TZ/8kqVq1amrWrJm6du2qt956S+7u7nrrrbec0ve///1v5eXlqVGjRvZ9/dChQ3riiSfUuHFjp9QoTdOmTeXv7++U/d3f31/u7u7XbX//97//rb179zp9f//55581YcIETZ8+XXfffbfat2+vESNGKC4uTq+99prT6oSHhyszM1MnT55Udna2kpOTdeLEiSve1y/2PRcYGKizZ8/q5MmTDu0r4/59NQiAV8jT01Ph4eFKSUmxLyspKVFKSopLzm1zFWOMRowYoY8++kj/+te/1KRJk+tSt6SkxH5uhzP06NFDO3fuVGZmpv3VqVMnPfDAA8rMzJSbm5vTap1z6tQpHThwQEFBQU7pr1u3bhfcmmDfvn0KCQlxSv+/t3DhQtWrV0+9e/d2et+nT59WlSqOf6W4ubnZZz2crVq1agoKCtKPP/6o1atXq2/fvk6v0aRJEwUGBjrs7wUFBdq6dWul2t/Phb/9+/dr3bp1qlOnznWp68z9feDAgfrqq68c9vXg4GCNHTtWq1evdkqN0nz//fc6ceKEU/Z3T09Pde7c+brt72+99ZbCw8Odeh6m9Nufp6Kiouu2v/v5+alu3brav3+/tm/fXuZ9/XLfc+Hh4fLw8HDYv/fu3ausrKxKtX9fLQ4BX4WEhAQNHjxYnTp1UpcuXTRjxgwVFhZqyJAhTqtx6tQph39xHjx4UJmZmapdu7YaNWp0zf0PHz5cS5Ys0ccff6waNWrYz3fw8/OTj4/PNfcvSePHj9edd96pRo0a6aefftKSJUu0YcMGp/5lXaNGjQvOW6xWrZrq1KnjtPMZn3zySd19990KCQnR0aNHlZiYKDc3Nw0YMMAp/Y8ZM0ZRUVF6+eWXdd999yktLU3z58/X/PnzndL/OSUlJVq4cKEGDx4sd3fn7/p33323XnrpJTVq1Eht27bVjh07NH36dD300ENOrbN69WoZY9SyZUt9++23Gjt2rFq1anXV+9/l9rXRo0dr8uTJat68uZo0aaKJEycqODhYsbGxTun/hx9+UFZWlv2+fOfCQWBgYJlnIS5VIygoSH/605+UkZGhVatWqbi42L6/165dW56entdco06dOnrppZfUp08fBQUF6fjx45ozZ46OHDlyRbcautxndX5w9fDwUGBgoFq2bOmUGrVr19bzzz+v/v37KzAwUAcOHNBTTz2lZs2aKSYmxinbMHbsWMXFxekPf/iDbrvtNiUnJ+uTTz7Rhg0bnLIN574fCgoKtHz5ck2bNq3M/V5Jje7du2vs2LHy8fFRSEiINm7cqHfeeUfTp093Wo3ly5erbt26atSokXbu3Kn4+HjFxsaW+UKTy33P+fn5aejQoUpISFDt2rXl6+urkSNHKjIyUl27di3zdlRa5XkJcmU2a9Ys06hRI+Pp6Wm6dOlitmzZ4tT+169fbyRd8Bo8eLBT+i+tb0lm4cKFTunfGGMeeughExISYjw9PU3dunVNjx49zJo1a5zW/8U4+zYwcXFxJigoyHh6epr69eubuLg48+233zqtf2OM+eSTT0y7du2Ml5eXadWqlZk/f75T+zfGmNWrVxtJZu/evU7v2xhjCgoKTHx8vGnUqJHx9vY2TZs2Nc8884w5c+aMU+ssW7bMNG3a1Hh6eprAwEAzfPhwc/Lkyavu73L7WklJiZk4caIJCAgwXl5epkePHlf0GV6u/4ULF5b6fmJiolNqnLu9TGmv9evXO6XGzz//bPr162eCg4ONp6enCQoKMn369DFpaWll7r8sn9X5ruY2MJeqcfr0adOzZ09Tt25d4+HhYUJCQsxf/vIXk5OT49RteOutt0yzZs2Mt7e3CQsLMytWrHDaNpzzt7/9zfj4+Fz1vnG5GtnZ2ebBBx80wcHBxtvb27Rs2dJMmzbtim4tdbkaM2fONA0aNDAeHh6mUaNG5tlnn72iv0/K8j33888/m8cff9zUqlXLVK1a1fTr189kZ2eXuUZlZjPGybfpBwAAQIXGOYAAAAAWQwAEAACwGAIgAACAxRAAAQAALIYACAAAYDEEQAAAAIshAAIAAFgMARBAqW699VaNHj26vIdhZ4zRI488otq1a8tmsykzM7O8h2S3aNEi1axZ0/7zc889pw4dOlxynQcffLDMTxO5EZTlMwFw/RAAAVQKycnJWrRokVatWqXs7GynPerPFZ588kmH54s6w3//+98KF3wBVF48CxjAdVNcXCybzXbBQ+TL4sCBAwoKClJUVJQLRuZc1atXV/Xq1ct7GABwUcwAAhXYrbfeqlGjRumpp55S7dq1FRgYqOeee87+fmmzQidPnpTNZrM/XH7Dhg2y2WxavXq1OnbsKB8fH91+++3Ky8vT559/rtatW8vX11f/93//p9OnTzvU//XXXzVixAj5+fnJ399fEydO1O+fHnnmzBk9+eSTql+/vqpVq6aIiAiHh9qfOzS6cuVKtWnTRl5eXsrKyip1Wzdu3KguXbrIy8tLQUFBGjdunH799VdJvx0uHTlypLKysmSz2dS4ceOLfmabNm3SrbfeqqpVq6pWrVqKiYnRjz/+KOm3WcRbbrlFNWvWVJ06dXTXXXfpwIEDF3yeH374oW677TZVrVpVYWFhSk1NdaixaNEiNWrUSFWrVlW/fv104sQJh/fPP9xZXFyshIQEe92nnnpK5z+F83Jja9KkiSSpY8eOstlsuvXWW+3vvfnmm2rdurW8vb3VqlUr/fWvf73o5yP99udq5MiRGj16tGrVqqWAgAAtWLBAhYWFGjJkiGrUqKFmzZrp888/d9iGoUOHqkmTJvLx8VHLli01c+ZMh343bNigLl26qFq1aqpZs6a6deumQ4cOlTqGAwcOqGnTphoxYoSMMTp06JDuvvtu1apVS9WqVVPbtm312WefXXI7AFyD8nwQMYBL6969u/H19TXPPfec2bdvn/n73/9ubDabWbNmjTHGmIMHDxpJZseOHfZ1fvzxRyPJrF+/3hjzvweud+3a1Xz55ZcmIyPDNGvWzHTv3t307NnTZGRkmC+++MLUqVPHTJkyxaF29erVTXx8vNmzZ49ZvHixqVq1qpk/f769zcMPP2yioqLMF198Yb799lvz6quvGi8vL7Nv3z5jjDELFy40Hh4eJioqymzatMns2bPHFBYWXrCd33//valatap5/PHHze7du81HH31k/P39TWJiojHGmJMnT5oXXnjBNGjQwGRnZ5u8vLxSP68dO3YYLy8v89hjj5nMzEyza9cuM2vWLHPs2DFjjDH//Oc/zQcffGD2799vduzYYe6++24TGhpqiouLHT7PVq1amVWrVpm9e/eaP/3pTyYkJMQUFRUZY4zZsmWLqVKlipk6darZu3evmTlzpqlZs6bx8/OzjyMxMdGEhYXZf546daqpVauW+eCDD8w333xjhg4damrUqGH69u1rb3O5saWlpRlJZt26dSY7O9ucOHHCGGPM4sWLTVBQkPnggw/Md999Zz744ANTu3Zts2jRolI/o3O/2xo1apgXX3zR7Nu3z7z44ovGzc3N3HnnnWb+/Plm37595rHHHjN16tSx/77Onj1rJk2aZLZt22a+++47+5+HZcuWGWOMKSoqMn5+fubJJ5803377rfnmm2/MokWLzKFDhy74TP7zn/+YwMBA88wzz9jH1Lt3b3PHHXeYr776yhw4cMB88sknZuPGjRfdBgDXhgAIVGDdu3c3t9xyi8Oyzp07m6efftoYc2UBcN26dfY2SUlJRpI5cOCAfdmjjz5qYmJiHGq3bt3alJSU2Jc9/fTTpnXr1sYYYw4dOmTc3NzMkSNHHMbXo0cPM378eGPMbwFQksnMzLzkdk6YMMG0bNnSodacOXNM9erV7QHo9ddfNyEhIZfsZ8CAAaZbt26XbPN7x44dM5LMzp07jTH/+zzffPNNe5uvv/7aSDK7d++21/jjH//o0E9cXNwlA2BQUJB55ZVX7D8XFRWZBg0aOATAso7t979rY4y56aabzJIlSxyWvfjiiyYyMvKifZ//5+rXX3811apVMwMHDrQvy87ONpJMamrqRfsZPny46d+/vzHGmBMnThhJZsOGDaW2PfeZbNq0ydSqVcu89tprDu+Hhoaa55577qK1ADgXh4CBCq59+/YOPwcFBSkvL++a+gkICFDVqlXVtGlTh2Xn99u1a1fZbDb7z5GRkdq/f7+Ki4u1c+dOFRcXq0WLFvZz3qpXr66NGzc6HLr09PS8YBvOt3v3bkVGRjrU6tatm06dOqXvv/++zNuYmZmpHj16XPT9/fv3a8CAAWratKl8fX3th5LPPyz9+/EGBQVJkv2z2b17tyIiIhzaR0ZGXrRmfn6+srOzHdZxd3dXp06drmpsv1dYWKgDBw5o6NChDr+DyZMnO/wOSvP7bXRzc1OdOnUUGhpqXxYQEOCw3ZI0Z84chYeHq27duqpevbrmz59vH1/t2rX14IMPKiYmRnfffbdmzpyp7Oxsh5pZWVm64447NGnSJD3xxBMO740aNUqTJ09Wt27dlJiYqK+++uqS4wdwbQiAQAXn4eHh8LPNZlNJSYkk2S+mML87n6yoqOiy/dhstkv2WxanTp2Sm5ub0tPTlZmZaX/t3r3b4dwwHx8fh2DnSj4+Ppd8/+6779YPP/ygBQsWaOvWrdq6dask6ezZsw7tzv+sJF3RZ3M1yjq23zt16pQkacGCBQ6/g127dmnLli2XrFfa7/9S27106VI9+eSTGjp0qNasWaPMzEwNGTLEYXwLFy5UamqqoqKitGzZMrVo0cJhHHXr1lWXLl303nvvqaCgwKH+ww8/rO+++04DBw7Uzp071alTJ82aNeuS2wDg6hEAgUqsbt26kuQw0+LM24ScCyHnbNmyRc2bN5ebm5s6duyo4uJi5eXlqVmzZg6vwMDAK6rTunVrpaamOgTZTZs2qUaNGmrQoEGZ+2nfvv1Fb79y4sQJ7d27V88++6x69Oih1q1b2y8OudKxlva5XIyfn5+CgoIc1vn111+Vnp5+RWPz9PSU9NvFGOcEBAQoODhY33333QW/g3MXjTjLpk2bFBUVpccff1wdO3ZUs2bNSp1l7Nixo8aPH6/NmzerXbt2WrJkif09Hx8frVq1St7e3oqJidFPP/3ksG7Dhg01bNgwffjhh3riiSe0YMECp24DgP8hAAKVmI+Pj7p27aopU6Zo9+7d2rhxo5599lmn9Z+VlaWEhATt3btX7733nmbNmqX4+HhJUosWLfTAAw9o0KBB+vDDD3Xw4EGlpaUpKSlJn3766RXVefzxx3X48GGNHDlSe/bs0ccff6zExEQlJCRc0S1jxo8fr23btunxxx/XV199pT179mju3Lk6fvy4atWqpTp16mj+/Pn69ttv9a9//UsJCQlXNE7pt0OVycnJeu2117R//37Nnj1bycnJl1wnPj5eU6ZM0YoVK7Rnzx49/vjjOnnypP39soytXr168vHxUXJysnJzc5Wfny9Jev7555WUlKQ33nhD+/bt086dO7Vw4UJNnz79irftUpo3b67t27dr9erV2rdvnyZOnKht27bZ3z948KDGjx+v1NRUHTp0SGvWrNH+/fvVunVrh36qVaumTz/9VO7u7rrzzjvts5ijR4/W6tWrdfDgQWVkZGj9+vUXrAvAeQiAQCX39ttv69dff1V4eLhGjx6tyZMnO63vQYMG6eeff1aXLl00fPhwxcfH65FHHrG/v3DhQg0aNEhPPPGEWrZsqdjYWG3btk2NGjW6ojr169fXZ599prS0NIWFhWnYsGEaOnToFYfZFi1aaM2aNfrPf/6jLl26KDIyUh9//LHc3d1VpUoVLV26VOnp6WrXrp3GjBmjV1999Yr6l347L3LBggWaOXOmwsLCtGbNmsuO84knntDAgQM1ePBgRUZGqkaNGurXr5/9/bKMzd3dXW+88Yb+9re/KTg4WH379pX026HTN998UwsXLlRoaKi6d++uRYsWOX0G8NFHH9U999yjuLg4RURE6MSJE3r88cft71etWlV79uxR//791aJFCz3yyCMaPny4Hn300Qv6ql69uj7//HMZY9S7d28VFhaquLhYw4cPV+vWrdWrVy+1aNHisrezAXD1bMacdzMqAAAA3NCYAQQAALAYAiAAAIDFEAABAAAshgAIAABgMQRAAAAAiyEAAgAAWAwBEAAAwGIIgAAAABZDAAQAALAYAiAAAIDFEAABAAAshgAIAABgMf8fc4D/5Gd+91sAAAAASUVORK5CYII=",
- "text/html": [
- "\n",
- " \n",
- "
\n",
- " Figure\n",
- "
\n",
- "
\n",
- "
\n",
- "
\n",
- " Figure\n",
- "
\n",
- "
\n",
- "