tutorial.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "1jC9CYFpLUFK",
    "outputId": "f38c5041-f815-472c-a9ad-a74ec2c056ee"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Requirement already satisfied: symusic in /usr/local/lib/python3.10/dist-packages (0.0.6)\n",
      "Requirement already satisfied: mido in /usr/local/lib/python3.10/dist-packages (1.3.0)\n",
      "Requirement already satisfied: pretty_midi in /usr/local/lib/python3.10/dist-packages (0.2.10)\n",
      "Requirement already satisfied: miditoolkit in /usr/local/lib/python3.10/dist-packages (1.0.1)\n",
      "Requirement already satisfied: music21 in /usr/local/lib/python3.10/dist-packages (9.1.0)\n",
      "Requirement already satisfied: numpy>=1.20 in /usr/local/lib/python3.10/dist-packages (from symusic) (1.23.5)\n",
      "Requirement already satisfied: packaging~=23.1 in /usr/local/lib/python3.10/dist-packages (from mido) (23.2)\n",
      "Requirement already satisfied: six in /usr/local/lib/python3.10/dist-packages (from pretty_midi) (1.16.0)\n",
      "Requirement already satisfied: matplotlib in /usr/local/lib/python3.10/dist-packages (from miditoolkit) (3.7.1)\n",
      "Requirement already satisfied: chardet in /usr/local/lib/python3.10/dist-packages (from music21) (5.2.0)\n",
      "Requirement already satisfied: joblib in /usr/local/lib/python3.10/dist-packages (from music21) (1.3.2)\n",
      "Requirement already satisfied: jsonpickle in /usr/local/lib/python3.10/dist-packages (from music21) (3.0.2)\n",
      "Requirement already satisfied: more-itertools in /usr/local/lib/python3.10/dist-packages (from music21) (10.1.0)\n",
      "Requirement already satisfied: requests in /usr/local/lib/python3.10/dist-packages (from music21) (2.31.0)\n",
      "Requirement already satisfied: webcolors>=1.5 in /usr/local/lib/python3.10/dist-packages (from music21) (1.13)\n",
      "Requirement already satisfied: contourpy>=1.0.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib->miditoolkit) (1.2.0)\n",
      "Requirement already satisfied: cycler>=0.10 in /usr/local/lib/python3.10/dist-packages (from matplotlib->miditoolkit) (0.12.1)\n",
      "Requirement already satisfied: fonttools>=4.22.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib->miditoolkit) (4.44.3)\n",
      "Requirement already satisfied: kiwisolver>=1.0.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib->miditoolkit) (1.4.5)\n",
      "Requirement already satisfied: pillow>=6.2.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib->miditoolkit) (9.4.0)\n",
      "Requirement already satisfied: pyparsing>=2.3.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib->miditoolkit) (3.1.1)\n",
      "Requirement already satisfied: python-dateutil>=2.7 in /usr/local/lib/python3.10/dist-packages (from matplotlib->miditoolkit) (2.8.2)\n",
      "Requirement already satisfied: charset-normalizer<4,>=2 in /usr/local/lib/python3.10/dist-packages (from requests->music21) (3.3.2)\n",
      "Requirement already satisfied: idna<4,>=2.5 in /usr/local/lib/python3.10/dist-packages (from requests->music21) (3.4)\n",
      "Requirement already satisfied: urllib3<3,>=1.21.1 in /usr/local/lib/python3.10/dist-packages (from requests->music21) (2.0.7)\n",
      "Requirement already satisfied: certifi>=2017.4.17 in /usr/local/lib/python3.10/dist-packages (from requests->music21) (2023.7.22)\n"
     ]
    }
   ],
   "source": [
    "!pip install symusic mido pretty_midi miditoolkit music21 matplotlib"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "oygHeqhtLUFM",
    "outputId": "09aad19a-783f-4288-8dd3-b568c9b983c7"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "--2023-11-26 16:29:28--  https://github.com/lzqlzzq/minimidi/raw/main/example/mahler.mid\n",
      "Resolving github.com (github.com)... 20.29.134.23\n",
      "Connecting to github.com (github.com)|20.29.134.23|:443... connected.\n",
      "HTTP request sent, awaiting response... 302 Found\n",
      "Location: https://raw.githubusercontent.com/lzqlzzq/minimidi/main/example/mahler.mid [following]\n",
      "--2023-11-26 16:29:28--  https://raw.githubusercontent.com/lzqlzzq/minimidi/main/example/mahler.mid\n",
      "Resolving raw.githubusercontent.com (raw.githubusercontent.com)... 185.199.108.133, 185.199.109.133, 185.199.110.133, ...\n",
      "Connecting to raw.githubusercontent.com (raw.githubusercontent.com)|185.199.108.133|:443... connected.\n",
      "HTTP request sent, awaiting response... 200 OK\n",
      "Length: 478769 (468K) [audio/midi]\n",
      "Saving to: ‘mahler.mid.4’\n",
      "\n",
      "mahler.mid.4        100%[===================>] 467.55K  --.-KB/s    in 0.03s   \n",
      "\n",
      "2023-11-26 16:29:28 (15.4 MB/s) - ‘mahler.mid.4’ saved [478769/478769]\n",
      "\n"
     ]
    }
   ],
   "source": [
    "!wget https://github.com/lzqlzzq/minimidi/raw/main/example/mahler.mid"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "id": "gPtW8uG4LUFN"
   },
   "source": [
    "# Load MIDI File\n",
    "\n",
    "* print each objects in symusic could get the corresponding summaries\n",
    "* Time unit is a quarter.\n",
    "* A *pyi* file have been generated, you can browse the members and functions of each class in it."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "lO5sIUDOLUFO",
    "outputId": "8c48c5c7-84c6-4106-97e6-e6f3fcfb36bf"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "note_num:  60411\n",
      "start_time:  780\n",
      "end_time:  1431755\n"
     ]
    }
   ],
   "source": [
    "from symusic import Score\n",
    "score = Score(\"mahler.mid\")\n",
    "print(\"note_num: \", score.note_num())\n",
    "print(\"start_time: \", score.start())\n",
    "print(\"end_time: \", score.end())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "v6b6VcS3bQLT",
    "outputId": "ff8f0c88-3bb3-475a-e008-2e1d95eda2de"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "TempoTickList[Tempo(time=0, qpm=20.00, ttype=Tick), Tempo(time=720, qpm=100.00, ttype=Tick), Tempo(time=7200, qpm=120.00, ttype=Tick), Tempo(time=12960, qpm=64.00, ttype=Tick), Tempo(time=15840, qpm=60.00, ttype=Tick), Tempo(time=18000, qpm=55.00, ttype=Tick), Tempo(time=18660, qpm=20.00, ttype=Tick), Tempo(time=18720, qpm=50.00, ttype=Tick), Tempo(time=51360, qpm=40.00, ttype=Tick), Tempo(time=53280, qpm=34.00, ttype=Tick), Tempo(time=54240, qpm=20.00, ttype=Tick), Tempo(time=54720, qpm=40.00, ttype=Tick), Tempo(time=55200, qpm=34.00, ttype=Tick), Tempo(time=56160, qpm=20.00, ttype=Tick), Tempo(time=56640, qpm=40.00, ttype=Tick), Tempo(time=57120, qpm=34.00, ttype=Tick), Tempo(time=60480, qpm=50.00, ttype=Tick), Tempo(time=155040, qpm=100.00, ttype=Tick), Tempo(time=203040, qpm=109.00, ttype=Tick), Tempo(time=220320, qpm=120.00, ttype=Tick), Tempo(time=224160, qpm=140.00, ttype=Tick), Tempo(time=228000, qpm=120.00, ttype=Tick), Tempo(time=229920, qpm=109.00, ttype=Tick), Tempo(time=231840, qpm=100.00, ttype=Tick), Tempo(time=235680, qpm=69.00, ttype=Tick), Tempo(time=241440, qpm=40.00, ttype=Tick), Tempo(time=243360, qpm=50.00, ttype=Tick), Tempo(time=250920, qpm=20.00, ttype=Tick), Tempo(time=251040, qpm=50.00, ttype=Tick), Tempo(time=266400, qpm=60.00, ttype=Tick), Tempo(time=277920, qpm=50.00, ttype=Tick), Tempo(time=335520, qpm=40.00, ttype=Tick), Tempo(time=337440, qpm=34.00, ttype=Tick), Tempo(time=339360, qpm=30.00, ttype=Tick), Tempo(time=341280, qpm=40.00, ttype=Tick), Tempo(time=343200, qpm=34.00, ttype=Tick), Tempo(time=344160, qpm=40.00, ttype=Tick), Tempo(time=345000, qpm=120.00, ttype=Tick), Tempo(time=432480, qpm=109.00, ttype=Tick), Tempo(time=444000, qpm=120.00, ttype=Tick), Tempo(time=522720, qpm=100.00, ttype=Tick), Tempo(time=561120, qpm=80.00, ttype=Tick), Tempo(time=562080, qpm=60.00, ttype=Tick), Tempo(time=562560, qpm=40.00, ttype=Tick), Tempo(time=562920, qpm=20.00, ttype=Tick), Tempo(time=563040, qpm=120.00, ttype=Tick), Tempo(time=589920, qpm=80.00, ttype=Tick), Tempo(time=591840, qpm=109.00, ttype=Tick), Tempo(time=622080, qpm=120.00, ttype=Tick), Tempo(time=684960, qpm=123.00, ttype=Tick), Tempo(time=686880, qpm=128.00, ttype=Tick), Tempo(time=688800, qpm=132.00, ttype=Tick), Tempo(time=690720, qpm=136.00, ttype=Tick), Tempo(time=692640, qpm=140.00, ttype=Tick), Tempo(time=694560, qpm=142.00, ttype=Tick), Tempo(time=696480, qpm=143.00, ttype=Tick), Tempo(time=698400, qpm=145.00, ttype=Tick), Tempo(time=700320, qpm=148.00, ttype=Tick), Tempo(time=702240, qpm=150.00, ttype=Tick), Tempo(time=704160, qpm=151.00, ttype=Tick), Tempo(time=706080, qpm=154.00, ttype=Tick), Tempo(time=708000, qpm=156.00, ttype=Tick), Tempo(time=709920, qpm=157.00, ttype=Tick), Tempo(time=711840, qpm=160.00, ttype=Tick), Tempo(time=713760, qpm=162.00, ttype=Tick), Tempo(time=719520, qpm=163.00, ttype=Tick), Tempo(time=721440, qpm=165.00, ttype=Tick), Tempo(time=723360, qpm=167.00, ttype=Tick), Tempo(time=725280, qpm=170.00, ttype=Tick), Tempo(time=727200, qpm=172.00, ttype=Tick), Tempo(time=729120, qpm=173.00, ttype=Tick), Tempo(time=731040, qpm=176.00, ttype=Tick), Tempo(time=732960, qpm=177.00, ttype=Tick), Tempo(time=734880, qpm=180.00, ttype=Tick), Tempo(time=736800, qpm=185.00, ttype=Tick), Tempo(time=738720, qpm=190.00, ttype=Tick), Tempo(time=740640, qpm=195.00, ttype=Tick), Tempo(time=742560, qpm=200.00, ttype=Tick), Tempo(time=750240, qpm=190.00, ttype=Tick), Tempo(time=752160, qpm=180.00, ttype=Tick), Tempo(time=754080, qpm=170.00, ttype=Tick), Tempo(time=756000, qpm=160.00, ttype=Tick), Tempo(time=757920, qpm=150.00, ttype=Tick), Tempo(time=759840, qpm=140.00, ttype=Tick), Tempo(time=761760, qpm=130.00, ttype=Tick), Tempo(time=763680, qpm=120.00, ttype=Tick), Tempo(time=765360, qpm=20.00, ttype=Tick), Tempo(time=765600, qpm=80.00, ttype=Tick), Tempo(time=823200, qpm=34.00, ttype=Tick), Tempo(time=834720, qpm=60.00, ttype=Tick), Tempo(time=836160, qpm=100.00, ttype=Tick), Tempo(time=844320, qpm=34.00, ttype=Tick), Tempo(time=849120, qpm=100.00, ttype=Tick), Tempo(time=850080, qpm=80.00, ttype=Tick), Tempo(time=850800, qpm=34.00, ttype=Tick), Tempo(time=854640, qpm=100.00, ttype=Tick), Tempo(time=856065, qpm=34.00, ttype=Tick), Tempo(time=860400, qpm=100.00, ttype=Tick), Tempo(time=866190, qpm=34.00, ttype=Tick), Tempo(time=866640, qpm=89.00, ttype=Tick), Tempo(time=868080, qpm=60.00, ttype=Tick), Tempo(time=868560, qpm=34.00, ttype=Tick), Tempo(time=871920, qpm=20.00, ttype=Tick), Tempo(time=872430, qpm=34.00, ttype=Tick), Tempo(time=872640, qpm=20.00, ttype=Tick), Tempo(time=873120, qpm=50.00, ttype=Tick), Tempo(time=873830, qpm=34.00, ttype=Tick), Tempo(time=874635, qpm=20.00, ttype=Tick), Tempo(time=874800, qpm=34.00, ttype=Tick), Tempo(time=876240, qpm=20.00, ttype=Tick), Tempo(time=878160, qpm=45.00, ttype=Tick), Tempo(time=880080, qpm=25.00, ttype=Tick), Tempo(time=882000, qpm=45.00, ttype=Tick), Tempo(time=883920, qpm=25.00, ttype=Tick), Tempo(time=885840, qpm=45.00, ttype=Tick), Tempo(time=894960, qpm=25.00, ttype=Tick), Tempo(time=895440, qpm=45.00, ttype=Tick), Tempo(time=897840, qpm=25.00, ttype=Tick), Tempo(time=900240, qpm=45.00, ttype=Tick), Tempo(time=903600, qpm=25.00, ttype=Tick), Tempo(time=905040, qpm=45.00, ttype=Tick), Tempo(time=917040, qpm=50.00, ttype=Tick), Tempo(time=947760, qpm=46.00, ttype=Tick), Tempo(time=949680, qpm=45.00, ttype=Tick), Tempo(time=951600, qpm=42.00, ttype=Tick), Tempo(time=953520, qpm=45.00, ttype=Tick), Tempo(time=959280, qpm=25.00, ttype=Tick), Tempo(time=960720, qpm=45.00, ttype=Tick), Tempo(time=969360, qpm=25.00, ttype=Tick), Tempo(time=970320, qpm=45.00, ttype=Tick), Tempo(time=972720, qpm=25.00, ttype=Tick), Tempo(time=975120, qpm=45.00, ttype=Tick), Tempo(time=987600, qpm=25.00, ttype=Tick), Tempo(time=988560, qpm=45.00, ttype=Tick), Tempo(time=994320, qpm=55.00, ttype=Tick), Tempo(time=1015440, qpm=50.00, ttype=Tick), Tempo(time=1021200, qpm=55.00, ttype=Tick), Tempo(time=1040400, qpm=109.00, ttype=Tick), Tempo(time=1071120, qpm=100.00, ttype=Tick), Tempo(time=1109520, qpm=102.00, ttype=Tick), Tempo(time=1111440, qpm=104.00, ttype=Tick), Tempo(time=1113360, qpm=105.00, ttype=Tick), Tempo(time=1115280, qpm=107.00, ttype=Tick), Tempo(time=1121040, qpm=109.00, ttype=Tick), Tempo(time=1145520, qpm=89.00, ttype=Tick), Tempo(time=1147440, qpm=69.00, ttype=Tick), Tempo(time=1149360, qpm=50.00, ttype=Tick), Tempo(time=1150320, qpm=30.00, ttype=Tick), Tempo(time=1151280, qpm=45.00, ttype=Tick), Tempo(time=1169490, qpm=34.00, ttype=Tick), Tempo(time=1170480, qpm=25.00, ttype=Tick), Tempo(time=1172400, qpm=45.00, ttype=Tick), Tempo(time=1174320, qpm=25.00, ttype=Tick), Tempo(time=1176240, qpm=45.00, ttype=Tick), Tempo(time=1180080, qpm=25.00, ttype=Tick), Tempo(time=1181520, qpm=60.00, ttype=Tick), Tempo(time=1184880, qpm=25.00, ttype=Tick), Tempo(time=1185330, qpm=40.00, ttype=Tick), Tempo(time=1193520, qpm=100.00, ttype=Tick), Tempo(time=1239600, qpm=94.00, ttype=Tick), Tempo(time=1247280, qpm=100.00, ttype=Tick), Tempo(time=1254960, qpm=80.00, ttype=Tick), Tempo(time=1297200, qpm=69.00, ttype=Tick), Tempo(time=1300920, qpm=20.00, ttype=Tick), Tempo(time=1301040, qpm=120.00, ttype=Tick), Tempo(time=1331760, qpm=60.00, ttype=Tick), Tempo(time=1347120, qpm=80.00, ttype=Tick), Tempo(time=1358640, qpm=69.00, ttype=Tick), Tempo(time=1360560, qpm=60.00, ttype=Tick), Tempo(time=1366320, qpm=50.00, ttype=Tick), Tempo(time=1370160, qpm=80.00, ttype=Tick), Tempo(time=1385520, qpm=120.00, ttype=Tick), Tempo(time=1408320, qpm=20.00, ttype=Tick), Tempo(time=1408560, qpm=80.00, ttype=Tick), Tempo(time=1423920, qpm=89.00, ttype=Tick), Tempo(time=1433525, qpm=89.00, ttype=Tick), Tempo(time=1435615, qpm=80.00, ttype=Tick)]\n",
      "KeySignatureTickList[KeySignature(time=0, key=0, tonality=0, ttype=Tick), KeySignature(time=7200, key=0, tonality=0, ttype=Tick), KeySignature(time=18720, key=0, tonality=0, ttype=Tick), KeySignature(time=51360, key=0, tonality=0, ttype=Tick), KeySignature(time=82080, key=0, tonality=0, ttype=Tick), KeySignature(time=103200, key=0, tonality=0, ttype=Tick), KeySignature(time=118560, key=0, tonality=0, ttype=Tick), KeySignature(time=155040, key=0, tonality=0, ttype=Tick), KeySignature(time=189600, key=0, tonality=0, ttype=Tick), KeySignature(time=220320, key=0, tonality=0, ttype=Tick), KeySignature(time=241440, key=0, tonality=0, ttype=Tick), KeySignature(time=279840, key=0, tonality=0, ttype=Tick), KeySignature(time=300960, key=0, tonality=0, ttype=Tick), KeySignature(time=320160, key=0, tonality=0, ttype=Tick), KeySignature(time=341280, key=0, tonality=0, ttype=Tick), KeySignature(time=350880, key=0, tonality=0, ttype=Tick), KeySignature(time=351840, key=0, tonality=0, ttype=Tick), KeySignature(time=371040, key=0, tonality=0, ttype=Tick), KeySignature(time=424800, key=0, tonality=0, ttype=Tick), KeySignature(time=472800, key=0, tonality=0, ttype=Tick), KeySignature(time=495840, key=0, tonality=0, ttype=Tick), KeySignature(time=538080, key=0, tonality=0, ttype=Tick), KeySignature(time=563040, key=0, tonality=0, ttype=Tick), KeySignature(time=591840, key=0, tonality=0, ttype=Tick), KeySignature(time=624480, key=0, tonality=0, ttype=Tick), KeySignature(time=625920, key=0, tonality=0, ttype=Tick), KeySignature(time=628800, key=0, tonality=0, ttype=Tick), KeySignature(time=630240, key=0, tonality=0, ttype=Tick), KeySignature(time=636000, key=0, tonality=0, ttype=Tick), KeySignature(time=637440, key=0, tonality=0, ttype=Tick), KeySignature(time=640320, key=0, tonality=0, ttype=Tick), KeySignature(time=642240, key=0, tonality=0, ttype=Tick), KeySignature(time=643680, key=0, tonality=0, ttype=Tick), KeySignature(time=647520, key=0, tonality=0, ttype=Tick), KeySignature(time=648960, key=0, tonality=0, ttype=Tick), KeySignature(time=650880, key=0, tonality=0, ttype=Tick), KeySignature(time=652320, key=0, tonality=0, ttype=Tick), KeySignature(time=658080, key=0, tonality=0, ttype=Tick), KeySignature(time=659520, key=0, tonality=0, ttype=Tick), KeySignature(time=661440, key=0, tonality=0, ttype=Tick), KeySignature(time=662880, key=0, tonality=0, ttype=Tick), KeySignature(time=664800, key=0, tonality=0, ttype=Tick), KeySignature(time=666240, key=0, tonality=0, ttype=Tick), KeySignature(time=673920, key=0, tonality=0, ttype=Tick), KeySignature(time=678240, key=0, tonality=0, ttype=Tick), KeySignature(time=682080, key=0, tonality=0, ttype=Tick), KeySignature(time=684960, key=0, tonality=0, ttype=Tick), KeySignature(time=692640, key=0, tonality=0, ttype=Tick), KeySignature(time=721440, key=0, tonality=0, ttype=Tick), KeySignature(time=765600, key=0, tonality=0, ttype=Tick), KeySignature(time=796320, key=0, tonality=0, ttype=Tick), KeySignature(time=823200, key=0, tonality=0, ttype=Tick), KeySignature(time=836640, key=0, tonality=0, ttype=Tick), KeySignature(time=850080, key=0, tonality=0, ttype=Tick), KeySignature(time=850800, key=0, tonality=0, ttype=Tick), KeySignature(time=860400, key=0, tonality=0, ttype=Tick), KeySignature(time=863280, key=0, tonality=0, ttype=Tick), KeySignature(time=866640, key=0, tonality=0, ttype=Tick), KeySignature(time=868560, key=0, tonality=0, ttype=Tick), KeySignature(time=874320, key=0, tonality=0, ttype=Tick), KeySignature(time=878160, key=0, tonality=0, ttype=Tick), KeySignature(time=891600, key=0, tonality=0, ttype=Tick), KeySignature(time=895920, key=0, tonality=0, ttype=Tick), KeySignature(time=899760, key=0, tonality=0, ttype=Tick), KeySignature(time=917040, key=0, tonality=0, ttype=Tick), KeySignature(time=934320, key=0, tonality=0, ttype=Tick), KeySignature(time=953520, key=0, tonality=0, ttype=Tick), KeySignature(time=954480, key=0, tonality=0, ttype=Tick), KeySignature(time=957360, key=0, tonality=0, ttype=Tick), KeySignature(time=959280, key=0, tonality=0, ttype=Tick), KeySignature(time=960720, key=0, tonality=0, ttype=Tick), KeySignature(time=961680, key=0, tonality=0, ttype=Tick), KeySignature(time=964560, key=0, tonality=0, ttype=Tick), KeySignature(time=965520, key=0, tonality=0, ttype=Tick), KeySignature(time=969360, key=0, tonality=0, ttype=Tick), KeySignature(time=970800, key=0, tonality=0, ttype=Tick), KeySignature(time=974640, key=0, tonality=0, ttype=Tick), KeySignature(time=980400, key=0, tonality=0, ttype=Tick), KeySignature(time=981840, key=0, tonality=0, ttype=Tick), KeySignature(time=991440, key=0, tonality=0, ttype=Tick), KeySignature(time=994320, key=0, tonality=0, ttype=Tick), KeySignature(time=1021200, key=0, tonality=0, ttype=Tick), KeySignature(time=1040400, key=0, tonality=0, ttype=Tick), KeySignature(time=1082640, key=0, tonality=0, ttype=Tick), KeySignature(time=1121040, key=0, tonality=0, ttype=Tick), KeySignature(time=1144080, key=0, tonality=0, ttype=Tick), KeySignature(time=1145520, key=0, tonality=0, ttype=Tick), KeySignature(time=1151280, key=0, tonality=0, ttype=Tick), KeySignature(time=1172400, key=0, tonality=0, ttype=Tick), KeySignature(time=1193520, key=0, tonality=0, ttype=Tick), KeySignature(time=1224240, key=0, tonality=0, ttype=Tick), KeySignature(time=1254960, key=0, tonality=0, ttype=Tick), KeySignature(time=1301040, key=0, tonality=0, ttype=Tick), KeySignature(time=1331760, key=0, tonality=0, ttype=Tick), KeySignature(time=1358640, key=0, tonality=0, ttype=Tick), KeySignature(time=1385520, key=0, tonality=0, ttype=Tick), KeySignature(time=1408560, key=0, tonality=0, ttype=Tick)]\n",
      "TimeSignatureTickList[TimeSignature(time=0, numerator=3, denominator=8, ttype=Tick), TimeSignature(time=7200, numerator=3, denominator=8, ttype=Tick), TimeSignature(time=18720, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=51360, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=82080, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=103200, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=118560, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=155040, numerator=2, denominator=2, ttype=Tick), TimeSignature(time=189600, numerator=2, denominator=2, ttype=Tick), TimeSignature(time=220320, numerator=2, denominator=2, ttype=Tick), TimeSignature(time=241440, numerator=2, denominator=2, ttype=Tick), TimeSignature(time=279840, numerator=2, denominator=2, ttype=Tick), TimeSignature(time=300960, numerator=2, denominator=2, ttype=Tick), TimeSignature(time=320160, numerator=2, denominator=2, ttype=Tick), TimeSignature(time=341280, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=350880, numerator=2, denominator=4, ttype=Tick), TimeSignature(time=351840, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=371040, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=424800, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=472800, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=495840, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=538080, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=563040, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=591840, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=624480, numerator=3, denominator=4, ttype=Tick), TimeSignature(time=625920, numerator=6, denominator=4, ttype=Tick), TimeSignature(time=628800, numerator=3, denominator=4, ttype=Tick), TimeSignature(time=630240, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=636000, numerator=3, denominator=4, ttype=Tick), TimeSignature(time=637440, numerator=6, denominator=4, ttype=Tick), TimeSignature(time=640320, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=642240, numerator=3, denominator=4, ttype=Tick), TimeSignature(time=643680, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=647520, numerator=3, denominator=4, ttype=Tick), TimeSignature(time=648960, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=650880, numerator=3, denominator=4, ttype=Tick), TimeSignature(time=652320, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=658080, numerator=3, denominator=4, ttype=Tick), TimeSignature(time=659520, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=661440, numerator=3, denominator=4, ttype=Tick), TimeSignature(time=662880, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=664800, numerator=3, denominator=4, ttype=Tick), TimeSignature(time=666240, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=673920, numerator=3, denominator=4, ttype=Tick), TimeSignature(time=678240, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=682080, numerator=3, denominator=4, ttype=Tick), TimeSignature(time=684960, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=692640, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=721440, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=765600, numerator=2, denominator=2, ttype=Tick), TimeSignature(time=796320, numerator=2, denominator=2, ttype=Tick), TimeSignature(time=823200, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=836640, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=850080, numerator=3, denominator=8, ttype=Tick), TimeSignature(time=850800, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=860400, numerator=6, denominator=4, ttype=Tick), TimeSignature(time=863280, numerator=7, denominator=4, ttype=Tick), TimeSignature(time=866640, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=868560, numerator=6, denominator=4, ttype=Tick), TimeSignature(time=874320, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=878160, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=891600, numerator=3, denominator=4, ttype=Tick), TimeSignature(time=895920, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=899760, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=917040, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=934320, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=953520, numerator=2, denominator=4, ttype=Tick), TimeSignature(time=954480, numerator=3, denominator=4, ttype=Tick), TimeSignature(time=957360, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=959280, numerator=3, denominator=4, ttype=Tick), TimeSignature(time=960720, numerator=2, denominator=4, ttype=Tick), TimeSignature(time=961680, numerator=3, denominator=4, ttype=Tick), TimeSignature(time=964560, numerator=2, denominator=4, ttype=Tick), TimeSignature(time=965520, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=969360, numerator=3, denominator=4, ttype=Tick), TimeSignature(time=970800, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=974640, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=980400, numerator=3, denominator=4, ttype=Tick), TimeSignature(time=981840, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=991440, numerator=6, denominator=4, ttype=Tick), TimeSignature(time=994320, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=1021200, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=1040400, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=1082640, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=1121040, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=1144080, numerator=3, denominator=4, ttype=Tick), TimeSignature(time=1145520, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=1151280, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=1172400, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=1193520, numerator=2, denominator=2, ttype=Tick), TimeSignature(time=1224240, numerator=2, denominator=2, ttype=Tick), TimeSignature(time=1254960, numerator=2, denominator=2, ttype=Tick), TimeSignature(time=1301040, numerator=2, denominator=2, ttype=Tick), TimeSignature(time=1331760, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=1358640, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=1385520, numerator=4, denominator=4, ttype=Tick), TimeSignature(time=1408560, numerator=2, denominator=2, ttype=Tick)]\n"
     ]
    }
   ],
   "source": [
    "print(score.tempos)\n",
    "print(score.key_signatures)\n",
    "print(score.time_signatures)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "FYkaBzW-LUFP",
    "outputId": "632e695a-c9c6-415e-d704-d71d019f4c24"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Tempo(time=0, qpm=20.00, ttype=Tick)\n",
      "KeySignature(time=0, key=0, tonality=0, ttype=Tick)\n",
      "TimeSignature(time=0, numerator=3, denominator=8, ttype=Tick)\n"
     ]
    }
   ],
   "source": [
    "print(score.tempos[0])\n",
    "print(score.key_signatures[0])\n",
    "print(score.time_signatures[0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "XdzfTIMULUFP",
    "outputId": "4eef7993-dfee-4550-84d6-ec1cd49cf5d0"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "track name:  PICCOLO\n",
      "is_drum:  False\n",
      "note_num:  1053\n",
      "start:  2160\n",
      "end:  1431755\n",
      "notes[0]:  Note(time=2160, duration=4895, pitch=106, velocity=111, ttype=Tick)\n",
      "notes[-1]:  Note(time=1423920, duration=7835, pitch=99, velocity=111, ttype=Tick)\n"
     ]
    }
   ],
   "source": [
    "print(\"track name: \", score.tracks[0].name)\n",
    "print(\"is_drum: \", score.tracks[0].is_drum)\n",
    "print(\"note_num: \", score.tracks[0].note_num())\n",
    "print(\"start: \", score.tracks[0].start())\n",
    "print(\"end: \", score.tracks[0].end())\n",
    "print(\"notes[0]: \", score.tracks[0].notes[0])\n",
    "print(\"notes[-1]: \", score.tracks[0].notes[-1])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "yrl6VKzwLUFQ",
    "outputId": "cbdb5d73-d13a-40fa-95c9-29d780adbedc"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "start:\t\t 2160\n",
      "duration:\t 4895\n",
      "pitch:\t\t 106\n",
      "velocity:\t 111\n"
     ]
    }
   ],
   "source": [
    "note = score.tracks[0].notes[0]\n",
    "print(\"start:\\t\\t\", note.start)\n",
    "print(\"duration:\\t\", note.duration)\n",
    "print(\"pitch:\\t\\t\", note.pitch)\n",
    "print(\"velocity:\\t\", note.velocity)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "id": "iaOZNxk6Yp88"
   },
   "source": [
    "# Batch Processing\n",
    "\n",
    "* sort\n",
    "* shift_time\n",
    "* shift_pitch\n",
    "* shift_velocity\n",
    "* clip(start: float, end: float, clip_end: bool)\n",
    "* filter_notes(func: Callable)\n",
    "* note_array"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "PYn1Fc9MYpGx",
    "outputId": "7763477f-c5ee-4799-acbe-a81fce3567a0"
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "IOPub data rate exceeded.\n",
      "The notebook server will temporarily stop sending output\n",
      "to the client in order to avoid crashing it.\n",
      "To change this limit, set the config variable\n",
      "`--NotebookApp.iopub_data_rate_limit`.\n",
      "\n",
      "Current values:\n",
      "NotebookApp.iopub_data_rate_limit=1000000.0 (bytes/sec)\n",
      "NotebookApp.rate_limit_window=3.0 (secs)\n",
      "\n"
     ]
    }
   ],
   "source": [
    "# inplace operation\n",
    "print(score.sort())\n",
    "print(score.tracks[0].sort())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "0u3hr6AGZFnR",
    "outputId": "db303e7b-8ea7-49d8-e6ae-f53400c4094c"
   },
   "outputs": [],
   "source": [
    "# method chaining\n",
    "score.shift_time(10) \\\n",
    "     .shift_pitch(-6) \\\n",
    "     .shift_velocity(-7) \\\n",
    "     .sort()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "GGcz6q0yZoh6",
    "outputId": "4e22d9d7-467f-41cb-df90-08118fe1897e"
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Track(name=\"PICCOLO\", program=73, is_drum=False, ttype=Tick, notes=[Note(time=12960, duration=575, pitch=94, velocity=111, ttype=Tick), Note(time=13680, duration=575, pitch=92, velocity=111, ttype=Tick), Note(time=14400, duration=575, pitch=91, velocity=111, ttype=Tick), Note(time=15120, duration=575, pitch=89, velocity=111, ttype=Tick), Note(time=15840, duration=575, pitch=91, velocity=111, ttype=Tick), Note(time=16560, duration=575, pitch=89, velocity=111, ttype=Tick), Note(time=17280, duration=575, pitch=87, velocity=111, ttype=Tick), Note(time=18000, duration=575, pitch=85, velocity=111, ttype=Tick), Note(time=18720, duration=1150, pitch=84, velocity=31, ttype=Tick), Note(time=23520, duration=575, pitch=79, velocity=31, ttype=Tick), Note(time=23520, duration=575, pitch=67, velocity=31, ttype=Tick), Note(time=24240, duration=430, pitch=60, velocity=31, ttype=Tick), Note(time=24240, duration=430, pitch=72, velocity=31, ttype=Tick), Note(time=24720, duration=190, pitch=62, velocity=31, ttype=Tick), Note(time=24720, duration=190, pitch=74, velocity=31, ttype=Tick), Note(time=24960, duration=190, pitch=76, velocity=31, ttype=Tick), Note(time=24960, duration=190, pitch=64, velocity=31, ttype=Tick), Note(time=25200, duration=190, pitch=77, velocity=31, ttype=Tick), Note(time=25200, duration=190, pitch=65, velocity=31, ttype=Tick), Note(time=25440, duration=1150, pitch=67, velocity=31, ttype=Tick), Note(time=25440, duration=1150, pitch=79, velocity=31, ttype=Tick), Note(time=28320, duration=960, pitch=79, velocity=31, ttype=Tick), Note(time=28320, duration=960, pitch=67, velocity=31, ttype=Tick), Note(time=29280, duration=770, pitch=81, velocity=31, ttype=Tick), Note(time=29280, duration=770, pitch=69, velocity=31, ttype=Tick), Note(time=30240, duration=960, pitch=67, velocity=31, ttype=Tick), Note(time=30240, duration=960, pitch=79, velocity=31, ttype=Tick), Note(time=31200, duration=770, pitch=81, velocity=31, ttype=Tick), Note(time=31200, duration=770, pitch=69, velocity=31, ttype=Tick), Note(time=32160, duration=960, pitch=72, velocity=31, ttype=Tick), Note(time=32160, duration=1535, pitch=79, velocity=31, ttype=Tick), Note(time=33120, duration=770, pitch=74, velocity=31, ttype=Tick), Note(time=34080, duration=770, pitch=79, velocity=31, ttype=Tick), Note(time=34080, duration=770, pitch=72, velocity=31, ttype=Tick), Note(time=36000, duration=770, pitch=79, velocity=31, ttype=Tick), Note(time=36000, duration=770, pitch=72, velocity=31, ttype=Tick)], controls=[], pitch_bends=[], pedals=[])"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "score.tracks[1].clip(1000, 100000, True) # start: float, end: float, clip_end: bool"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "eiU2D9XDc30t",
    "outputId": "e284d341-c07d-49dd-e64f-6aba03ae9566"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "frame_pianoroll float32 (1, 128, 5966)\n",
      "onset_pianoroll float32 (1, 128, 5966)\n",
      "both_pianoroll float32 (2, 128, 5966)\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "<matplotlib.image.AxesImage at 0x7f6e78bb0190>"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAigAAAGhCAYAAABLWk8IAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABr3klEQVR4nO3deXxcV303/s85986+al9seUnsxHH21RgCBOJikpAmJbTA45YU8pAHSHgI4WHx0yYUSjGk/ICGhoTSltCnDSm0JIFAAqmzNYnjOHac2Fm8xI537dKMZp977/n9caWRRhpJI2kk3ZE/79dLL1sz59575s5o7vee5XuEUkqBiIiIyEHkfFeAiIiIaDQGKEREROQ4DFCIiIjIcRigEBERkeMwQCEiIiLHYYBCREREjsMAhYiIiByHAQoRERE5DgMUIiIichwGKEREROQ48xqg3HXXXVi2bBm8Xi/WrFmDF154YT6rQ0RERA4xbwHKv//7v+PWW2/FV7/6VezYsQPnnnsu1q9fj87OzvmqEhERETmEmK/FAtesWYOLL74Yf//3fw8AsCwLbW1t+OxnP4uvfOUrE25rWRaOHz+OUCgEIcRcVJeIiIhmSCmFgYEBtLa2QsqJ20j0OapTkVwuh+3bt2Pjxo2Fx6SUWLduHbZs2TKmfDabRTabLfx+7NgxrF69ek7qSkRERJV15MgRLF68eMIy8xKgdHd3wzRNNDU1FT3e1NSEN954Y0z5TZs24Wtf+9qYxy/FldDhmrV6EhERUeUYyOMZ/BahUGjSsvMSoEzVxo0bceuttxZ+j8fjaGtrgw4XdMEAhYiIqCoMDiopZ3jGvAQo9fX10DQNHR0dRY93dHSgubl5THmPxwOPxzNX1SMiIqJ5Ni+zeNxuNy688EJs3ry58JhlWdi8eTPWrl07H1UiIiIiB5m3Lp5bb70V119/PS666CJccskl+P73v49kMomPf/zj81UlIiIicoh5C1A+/OEPo6urC7fffjva29tx3nnn4dFHHx0zcJaIiIhOPvOWB2Um4vE4IpEILsM1HCRLRERUJQyVx5N4CLFYDOFweMKyXIuHiIiIHIcBChERETkOAxQiIiJyHAYoRERE5DgMUIiIiMhxGKAQERGR4zBAISIiIsdhgEJERESOwwCFiIiIHIcBChERETkOAxQiIiJyHAYoRERE5DgMUIiIiMhxGKAQERGR4zBAISIiIsdhgEJERESOwwCFiIiIHIcBChERETkOAxQiIiJyHAYoRERE5DgMUIiIiMhxGKAQERGR4zBAISIiIsdhgEJERESOwwCFiIiIHIcBChERETkOAxQiIiJyHAYoRERE5DgMUIiIiMhxGKAQERGR4zBAISIiIsdhgEJERESOwwCFiIiIHIcBChERETkOAxQiIiJyHAYoRERE5DgMUIiIiMhxGKAQERGR4zBAISIiIsdhgEJERESOwwCFiIiIHIcBChERETkOAxQiIiJyHAYoRERE5DgMUIiIiMhxGKAQERGR4zBAISIiIsepeICyadMmXHzxxQiFQmhsbMS1116LPXv2FJXJZDK46aabUFdXh2AwiOuuuw4dHR2VrgoRERFVqYoHKE899RRuuukmPP/883jssceQz+fxvve9D8lkslDm85//PH7961/jF7/4BZ566ikcP34cH/zgBytdFSIiIqpSQimlZvMAXV1daGxsxFNPPYV3vetdiMViaGhowH333YcPfehDAIA33ngDZ5xxBrZs2YK3ve1tk+4zHo8jEongMlwDXbhms/pERERUIYbK40k8hFgshnA4PGHZWR+DEovFAAC1tbUAgO3btyOfz2PdunWFMqtWrcKSJUuwZcuWkvvIZrOIx+NFP0RERLRwzWqAYlkWbrnlFrzjHe/AWWedBQBob2+H2+1GNBotKtvU1IT29vaS+9m0aRMikUjhp62tbTarTURERPNsVgOUm266Cbt378b9998/o/1s3LgRsVis8HPkyJEK1ZCIiIicSJ+tHd988814+OGH8fTTT2Px4sWFx5ubm5HL5dDf31/UitLR0YHm5uaS+/J4PPB4PLNVVSIiInKYiregKKVw880344EHHsDjjz+O5cuXFz1/4YUXwuVyYfPmzYXH9uzZg8OHD2Pt2rWVrg4RERFVoYq3oNx0002477778NBDDyEUChXGlUQiEfh8PkQiEdxwww249dZbUVtbi3A4jM9+9rNYu3ZtWTN4iIiIaOGreIBy9913AwAuu+yyosd/8pOf4M///M8BAN/73vcgpcR1112HbDaL9evX44c//GGlq0JERERVatbzoMwG5kEhIiKqPo7Kg0JEREQ0VQxQiIiIyHEYoBAREZHjMEAhIiIix2GAQkRERI7DAIWIiIgchwEKEREROQ4DFCIiInIcBihERETkOAxQiIiIyHEYoBAREZHjMEAhIiIix2GAQkRERI7DAIWIiIgchwEKEREROQ4DFCIiInIcBihERETkOAxQiIiIyHEYoBAREZHjMEAhIiIix2GAQkRERI7DAIWIiIgchwEKEREROQ4DFCIiInIcBihERETkOAxQiIiIyHEYoBAREZHjMEAhchoh5rsGRETzjgEKkcNooRAgtfmuBhHRvGKAQuQwwu+DkGxFIaKTmz7fFaCTkwwEAABWMjnPNXEQIQClBv8v7VYUy5x4m3LKENFJSYtGoBa3QFgW1IHDsDKZsrfVly2BSiRhdvfMYg0nxhYUmhdq1TKo1afMdzUcQ7jckH4/AEApBeF2QQsGJt1OCwc5ZoWISkqvWYm3vurCvr/0QyxdXP6GQuDA9YvRv27l7FWuDAxQaF6kWwNILfLNdzUcQ7h0CK/H/sVSEC4d8Hgm387ns1tbiIhGiZ3iwva1/4TNl/498g3B8jcUEq3vPIreM+b3u4VdPDQv9JQJdRLe+AtdhzKMsU9YFmCahdYQZVoQWumuG+n1DjfV6jqEFFDWbNWYiKqVO67wD7HT0GcEoGUMqHI3VBbeOlGHYP8sVq4MDFBoXnj3tAMASlyqFzQZCsHs7x8eazLIyuUhTAvC7QaUBZVOQ+VyY7YXug7Z2gzrwFsAAOX3jh/0ENFJrXZbF/7t/7sCwgTqD+1H2aPVlELbfTp8J/oxn/c+DFBoXhhHj813FeaF8PuAWBxQo74qLBPKMu1xKErZAUepoEPTYEaGx6Yojw5onJJMRGOZe99E7d437f9PcVvPI9vmNTgBqjxAyV9+PjyHkzD3H4R2+gqkl0Xhf60dxpGjFT0O71DHJy4+G1p3HEhnoEwLyOcggkHAMKCyOahsFgBgpVIA7O4JWV8HwL77V4ePlRxZLlxuu4yRt8dYLJCZKiqbRcn+GCEgNM3u6gEAqUF6PRC6DiuVGv78Wapo5JhI56DMhXFuiIhGquoApeZLh3H412ehcf9BHP1AI1Z/8A0cuus0hO+rbIAiQyFYAwMMUkaTGvbc4EP9C0EEOgxoGQuu3gwGVobgSprwdGWg9SUhTKvQJSGbGtD3tkUAgGSLxOIHcrAOHSnerxCQkRAgJNTAAOBywRoYmOMXNzusgcSY7h0AEJpmd//E4tACgPR5IVoaYYX8kG8dh9nXB8AO2IQcjlBUexesXH7O6k9ENFeqevj/Py/bjMQS+8s+cYqJ+5c/jkTrLLwkt4szJUoQUuDis95Eok0g2agj3eBCvsaLRItEslFHrtYLFfDCCnoL2yi/F8kWiWSrxMByCypQeiaP8HggPG5A0yD0qo6jiwy1KI0hpP05G2wpEi4dKuiDEfXYjxd2oIoGF1sDAwumdYmIaKSq/ub/2IH3IXjI/rYOHtDwoTfXIXis8r1mQojSzfInOWUpbNt1KuqPKAQ6B1tQ+jIInnDBlTDh7s1AJNIQlir0ZYpUBoETg78pCZFMl953NmtftE1zYbdcXXI2svVeKE0g+NIxexaPrtszegDo/VlgRAuJ9HphuLXqvrOYAX3xIvSvXVyU+8WVNGF4JZQ2/Jg0FGRewfBJ6CkLvt/thMqPHXQMAFpdLeKXrYSSQGh/AolTg5B5hcDvXplSYqu5Jr1eCK8HZn8MAKCfsgz5lihMrwYowNM+APO1vRU5ltbUCLOjsyL7Imex3nk+XL0pmK/ume+qjFHVAcrAdxah5fAJmAAWP9yJnt3LUPva0crPDBECyip7gtbJwzJx+j+nx4xBiXYVj0EZeeasji7UPGsHKFG/F1ZH19j9KgUrZnfpKCNfdIFeaA5+MIjo2d2QQiHwOQ+E7oLweaFyeShNQB44XtS9JetqYQR0uEdmnT2JJM9pxYVf3gGXGG412tq1DKdHOxF1pQqPnchE0JUJ4qLoMWzrXgr5bABmX+kAxVreijO+tBsBPYvf//ISvOMPX8aRZBRyawRWu4MDlGgEqi4KDAYovW9rRsdaBa0uC9OQiLxQj8YKBSj5VYsgGaAsPEJg/wYXIq/VoenV+a7MWFUdoLg2vwRT2M3f5p79cO85+aatzje1bdfYcz74hVmKlcnAKmMGT9Hd7ugZLwuItSSNTyx/DgDwK/3tgBRQ+uCsHCHGTkl2u4paCk42mRoNf9vyHDxiuNvr/uBBvNd/FI3a8Oymvfkk9uXrcJU/g5+H3sS97jXj7tMIufE3rb9HjfTiobYLsKn193g5F8Z3PVfO6muZMY8bZsBd+DVdL1GzrBsXNR1BLO/DS0dOr9ihcmEXvCdpULzQLTu1A8f6Wue7GiVVdYAyV1RmnJkXVDWErkP6/TDj8fmuShHvK358O/t+AMDqgXbANCESKVi5HGQqB1PI4gAtb0Aayu7+UmZx0raTQPB4Du9++aPQ5PDfY3tXBJFICj73cEtbPO1FJu3GN2vjaO+O4PTsgXH36epJ4YqdH4dbNxF5Vcf6RR/HQNKLFakTs/paZkql7EHoQ5+O8GEDnTvq8VhtDUReoG5f5Y7l7UgzOFmgjr3YikhlGtoqjgFKGayBAf5xVjnh8UDU1wIOey+X/vw4rLC9Bo/Z1Q1lGDB7+qBME3osAaFpUCMGwapsDjI3fHGWUWd3Q1Sae+dBRL7eBjViDMrKTA7K5YbShgdjhwwLsPKw3CGszORgxhPj7/TNI6j/61OgNInaE0dhbI2gwTRhDc6cciqrrw9ixGKbwWcPIvhaGMrtgrAsiFiiYi3K2huHppxHg6qAUlhxbxdEMu3I3gcGKGVY0IM058N8NBVLCeVxFVoepmyW6mwMTr8eaah7SxkGIEd151gmhGENp7d3ucZsv5CZfX3A830YeVaG3pXxHpvsXbOSSeCFXRAY7CJ+a/JtnEAZRtF3k9nVBXSVGNNVAU5reaRpkNpwT8CI7zJzz/7xtxnKzzRImWZZ34MTzbwUSpU9FoMBShmk319INEYzJATk2adDHOko5PaYE/k8RDw5ra46fdkSZJfVQ3tyx7QOLTweqLxRmA6s1dfB7Okd84cu/X57WrUQ9r8BP9DTW1RGZbJwdSVgDs7ysSJB4JjGqcZENKHsFRfAFTegpXJQ28sbESvOW40Tl0UKvzdvSQDPvwIZCkF43DC7e8Zso9XVov2PT4c0AWECuYiAMAAtq2AEBBBPA//4n2Udf9YDlG9961vYuHEjPve5z+H73/8+ACCTyeALX/gC7r//fmSzWaxfvx4//OEP0dTUNNvVmRYZDjFAqRChu9Bzfg3qswYwhwGKlc1CdXROqxUkvaIBx9/lxtInp3fsoQBXZe0gQjU3QMQGigcCCwEZDtlTjHUNyuOGGfRAHSwOPKxUCuKtI4XXkWsMwL3fBSvDAIWIxiEEDl0t4Dvqg7fHi4bt5W3WdUkYt3/6Xwu//5X2p2h9HpC1UVihAFAiQEFTPS6/8XkkTQ+6swF8uHEbjudr8HqqBZdHXsPzna3Y/Y/lHX9WA5Rt27bhRz/6Ec4555yixz//+c/jN7/5DX7xi18gEong5ptvxgc/+EE8++yzs1md6eNaJxVl+ADoc3xOh9a3mQbLLZEPzqDRX5MQQhS6DZRrnNeuaXZwotv/QsqxAZVSRcnelC7scqXIweMoy1Hjboho7mmhPIygBiNV/ixAwyfwbt/wYHFzKK+mpgF66e8dpUtcHDyAAdOHE54o3u07gaOuLvhlFu/2nYAVKJ37qpRZC1ASiQQ2bNiAH//4x/jGN75ReDwWi+Gf/umfcN999+G9730vAOAnP/kJzjjjDDz//PN429veNltVmjY1YiAazYwyTdTuyUL0jj8VeTbIQACQclop831HB9D4YnT8fXu9Y8YDFFGqOIiQGNvVpJT9OdM0QGoQugaZD8CaZB0oVzwH4fEAJVr45OqVw4Mljx2f4BUS0YKmFML/7YWvV8EdL/9Grfb1HC7+/ecKv7ftsrdV8QHIvFFyMUHZO4CNv/0ohAFIQ+CfGi4F8hJaUuIbUQOqOw/gtrKOP2sByk033YSrrroK69atKwpQtm/fjnw+j3Xr1hUeW7VqFZYsWYItW7Y4MkAxYxwgVjGWCdcLe2COl/J9lshwCNC0aQUoas9B1Bz2jjuLQfh8QDo9fiBhqaLMp+Ml/jNj8cKSCkIKSMOEcLsnDFC0nkRxKvwRYmdGoeUUAkc9AAMUopNa83/stRNomlbZqxR7n9uD1XtqCr+r3n6YAMzefghZ+ibTbO/Aqu8MbaAAl27/a1qArsHIp3Gk5JZjzUqAcv/992PHjh3Ytm3bmOfa29vhdrsRjUaLHm9qakJ7e3vJ/WWzWWRHXNDiVTSiXAYCEC69kI66pJMsAZKVTELo+sSrRA91T1gVyvXh9QwnQJsilc1OHFBpcvxuFmBsa4k5znstZKGssuz/F1pHxvuMKAVEw9CGsu1qg/UwLeT9ApYLsHSJkze1G9ECI4R94zLFm7xSA1pHkqEQkM8XfddaAwOlb+osE5ClwwdlGDAmSMZpqPIzg1c8QDly5Ag+97nP4bHHHoPX6518gzJs2rQJX/va1yqyr+nQQqFpT7Pr+dA5SDcItH5367gzLWQwCCtRepXbhUqedoqdKbXU+g9CQG9tBpSC2dGJ9OVnw/ObscHuVKRPrUe21oXQ4WNT/sMuSQh78GsyabegmBakEPaU1VFUJmsvAzBI6+qHUWI2kd7aDJW0u2qE114kML+oFtrOLISuF30GhcsNoUkorwdv/mkN9GQjYAHZOgvCFHDH7ZHzSgrUvJIt+46JiJxNi4SRvXAF9M1ljnQtU8efnoXwYaO871ohIFcuh/l6BbMBllDxNce2b9+Ozs5OXHDBBdB1Hbqu46mnnsKdd94JXdfR1NSEXC6H/v7+ou06OjrQ3Nxccp8bN25ELBYr/Bw5Um4DUWWIUHDa23a/N4vTr9kL4Ro/FpQB/0m3WnJqWQTJ5ZGSzwlNg9lcA6shCqHr6LjYVdxFMg3xpW70rpIQbvfkhcsgNA0iaKdWV143hMdd+H00K5MpmrFjHG8vGYyaDVGIcBAiEoLVEIXyexE71Q/h8475DAqXDuHzQflc+OQHfo8Pf+hJXP2h5/DwB7+Lf//Q3+HDH3oSmUYL6SV5iNTcdqcR0ewRoRA6z/dUeKcC4soedFxUZl4lIZFaFq1sHUqoeAvK5Zdfjl27dhU99vGPfxyrVq3Cl7/8ZbS1tcHlcmHz5s247rrrAAB79uzB4cOHsXbt2pL79Hg88Hgq/IZMhTX9+09lSKSMSS6KM7z4ViUFiIkWYFTKHjiqFEQF8uRJw/6pGCELyYjEUOvIRK9nRBdNIcnaaPpgV5EQ9o+l7HOk66XHmUgBJYC+fAApy42cpaPX8mLA8iFueKFlBJQmIfJMNLjgSG34psdSw12DFUoqqdXUQAQDMNvtBQJVPgfp9ULWRGG0d5xUrb2OoxSUZreijrdC91TJYBBKCYgpZCuY8Pu7QioeoIRCIZx11llFjwUCAdTV1RUev+GGG3DrrbeitrYW4XAYn/3sZ7F27VpHDpAFACs+9YGVQ1p/p6O9bika8h3jF5po/MIC5d/fY3fxlHhOmSb04z2AUjDyBpqfz874C7F2dxyhI157XaUKEJqE8ttdmKovBpXLjV9HqUEG/HZfrhDQWpphHDk6ttzQ9kpBJDNAfxzRPV4Ivw+5xbWQBw8NF83lIAIBQAr8+r5LIUxAWMDjNRdAWIBrAGh9Iw/TI6E4yHvB0WoiUIsaAQAia0AkUvYAxEOVaV3u+qNV6LnQwqrvu+2/071vQp15KvZ9KIxT/zrOvFDzyIrF4RoA5PI2mHvfrMg+jfNWIPecDy3byvx+VBZ8b/bM+vIH85JJ9nvf+x6klLjuuuuKErU5ValxBeUKPfQSwkLAmizT50m2GKH55qHxn1TKvksb/L/7yZdnnHpcvfQGXFJUbtkCKaF8dsvYZBlxhaZB+H32OkBCwmyMABMFKIC9NkZPL2QqDbW8DelGN0Z2ICnDgNA1KCGw+J5dwGBm2UJrnFJ29loAZoXussg5RDiEVFsIAKCnTLi6dCivDlQoQOl+ex4/e++PcNvP/yeUEND2AollQfz5VY/jme8sKjmtneaGNTAAd0Ih1xqBVqFF/vpXetHybBrall3lfdcqBfPA4cocfAJzEqA8+eSTRb97vV7cdddduOuuu+bi8DM2k6Y0lc1O/oab1ZsFVF/aBpVIwuzpLbkkgNB1QMjC+RO6Djk0VsPlhtXXVzpoGHGxrkhQYZmlY0A5/TTxIj8cFAjd7oIZ73MiBgMHIcW4XXoiawLGiLooZXcv5g1o2cHuIZcbysgXusBkzqieAdYn2Wy1WWWYkPnBz0TegjBNwJAVW0NIJnQ8nz4VImcVPq5aVmH3QOvEXZk0J1wpyx7039Bgr8E0Q1pW2Snwp/BdK90uWNnZTQJ58vUtTJUQ0Ba3zOohrLizVtgtm9Sw/38uRvyylfaF/rRlY4poDfXQmhuLfk+94zSk3nEauj6wAlpT45ht5pIWCU9rO5XLAZ32tD3p80Fb1AytqWGcwtZwUCLkuNOMRXsXrP4YrP4YIKU9U8ftBrp6EXjD/hLSmhog3G57irZSkEe7quKzYwem0x9sTsWs/hh8ezvh29sJ91tdQGcPxLHKLRS45FET995zJVyHu+A6ZO838EYX9v7zKnbvOEB4Vw/a10h0X7WiIvur3dlnf5dMxYplFZt0MB4uFjgZIWHWhYCDs3eIav2DF1LgjHcdwL7UqQhKgXRrAJ6dxWVUOAjl0TGUmUdFQ+hbYbc2JJZZaHg+AIw/ZX52CWF3vUxjTSBlGPaCf7AXAzRrQnaLSqmem5F3nFJAKFXyTndkngIZCduDIF26fZzBWW9WNATR129vb6mK3D3NBaHrEB43MP3hXDTCuPkpKsT9+x1o/i8BY+RqyW++hfq3jnB1dwcw9x1EcJWOHlcUNZMXn3x/r+2d8o1OtjUIz0G9MmkbxsEWlMkoC8KY3fEhEy1N7WTKUni9vQnuwRx0empsV4nIGxCZ4cQ8IpuHO67gjit4eiREtnTSHqHrkIN5dKTfP9zSIoSdUKhSr2Hoy1YIOynaFMhReX5EbvwERGqo68ZS479mj8duNdF1O+GaEHb3n9QKS56L3GD3TpU1s6uhTJJzRQhoZ54OXHI2tGjp6eyljPxsVfJzVnUsc2wgMoP1rKjCLBP97SF4eit0CZ9GK6zMWxBLWoFLzobePDsL/TJAmczQjIpZJKfZzTDvlIXmf/Wi+dk+KNOEZ9/YTMBWVw/Q2T28SXcvGrb2oGFrD9r+KwHV3Vty1zIagVjcYl+cly5C/J3LB7PPuqBOW1Kh+qvCDBfhdttdUVOY8i3aWu3dGAaEZRW9ziKWCTV4t6uMPHCis2QxrbUZWl0NZE0NlEu3s0VmstAi4eHPSGc3rEzW3o8UVRPcqrwBKzF3a1oJ3YU3PhNF7K/SyJ13atnbqVXLClmMR/6fyGmW/6dC238l5u34rt40jl7RgMTXk+h97/JZOUZ1fLvNM5GZ3VkQwuebvJATKQXvr18oZCkttSDd6GZoMx4HXhue9jrR+jZmxA8hBYzaAGLLNASFhNAkMg1+VKrncyits9A0e/nwKTBrBstb9kCxiZYzKKSPnqCcFfZDCgFhmlBSApqEMgxIj9vOhQKM3VbTgGq4q7VMqGkORp4OoUn84du249vNW3BZy2dR7i1AptEL72CemnSTD77xctYQzTPX71+c1+OLZAYDq/J46ez/wNmn3lz239hUMEApg8rNboCi8uWvTXCyULkcZCoH01LQEln4unyAsqCUnZysEoSuQ/r9dtBkWXZCMyEBVd6FVEvm7ADLNCHSOchAYNwp6UMzwYTLDq1KzfYR6RyQy9vdTbm8nR5f0+xm9cEmWK2mpnhac5V19cwVZVp49M0zAACe/vIDI1fcKIwZcsfyJRd1JCL7xt17zIWvdFwIX+fs/J0wQCmDNdFCf1Ww/2pk9ccg0xnAMiEOt6M+lYVpmvaYDIGJFxoskwwGoJa0AK/EYeXy0JNpCE0r/07/mN2lZeXy0Dq67S6pPftLHysagdnVBRmxV1U2O8Z286gTnVCmCVkThYrF7dYTvx9WbAAYnJ5srFoC8Xy/nclWCLurh8ZQRh4t/+TBtshFiO46gnI/Ka69x2AOvv8j/09ExazePix9OIpn31iDple6ZyVpGwOUMszmKOW52H81GrmCsNnXVzzTRqAyaxe53DAiPnsglmVCZXN2IFCmQneLZdqtMKuWjltWeNyD/3rs5cdLGOoOk+GQPbNL0yDcruHWFiGQq3XDg8F8KkBVTDGeF0rB/bsX4QbKDk4AFAWOpYJIIrJZqRSw/VWEto/fVT9THCRbBi08u4NYpzp7pKoJAX3ZEujLl047B4qWsSB9FVgpO5+DlrEvX0LXgfrolLpMRn4uyn0PVTYLlEq3L4S91kkoBBX0Q/r9hRwDwuOB8HggPR4IY0T9XK4xM4mcaug1VA2pQV/UOt+1IJoRfdkSJK9bg/S1l0CrmdqEZOvS86CdPjbPigwEYF52AYSuQ2toKBxnNma9MUCZjBDAotmZQjVETmEaZLUTmob+S1rR8/YW5FYvntY+XPEcRG10xnWxEkloPfb6OMLnQ+qU6NS6TBYPr75d7ntoxQZg9vWPeVy43ZD1dRAtjTAawhD1tZABP4SUkOGw/RONQE/a41GUpaD8Xsj6uqqYaSKjEWg10fmuRtmk24XkeYvmuxpEM9J/SSsu2LgDy778BqxlUwi4hcD+P3PhxOVjbyJlUwMO3qAgfD7kV9nf4X1rWiFaK3+dZIAyGSFhBWb3zk+4ylzieiEQEql6iXSDRDY6vdctsnkoz8zn8SjDALLDKfizYW1KXSZWwFOYlixcLgiFSacpq3yuZJee0DQorxvK74Hltf8PXbczyrpdEG4X4HJBDGWhVRaUy17lWEyhW2q+CJcLqMB7Nmc0Dela9oBTdUvXCXyvZSu+s/gRmMEp/P0JidNWnECqZez3ofJ58I4Vb0LoOnIR+zs8XS9h+St/neRf4GSUBRlLzeqqjSozu3lWHEVZCB8xYLkEfB1Te93S64WVyUAmMsNTfGdAuNxQg8nVVC6HwImpzdaSsRTMEQHN6N9HUuk0AECrq4VKpoanHQ89b5qQyTTgcUNL5iEGUvbnYuT+hIAeSwPRCMxYHHIgaQdqlRiPM8tUJlMV9RyiDAPB41xkkapb6KiJ9+y+DhlDR33vFK5jysKB7W2oKbFYskik8My2M3B6bje8nenB4xjQYskpjfcqBwOUySgFdWRsfo9KOplm8SjDQOD5g/Zg1Gx2SoGfiISBTAaqtw/55XUz/vDKYAAqlQKUgpXOwP3ywSnVRx09MeHvI1mDCeFUWzPEiW5gdICSy8Hs7oXm80JPZ2F2dpVcRFIaBtDSCPTHYLZ3QqxcDuHSp72Y5Vypts+4yuXg2Tm1zwOR0wSePwijpxU+S0G9VSLaGI9SWHlvD0R8bNBhtnfi9H8IwEpnoL1xCCaA4HMHZ2XphaoOUISuA9bgCqlCDE4RVYUm76H/K0sNLtg2eAc3xamDo+92p1ZJYR93gmNWbfro8VannWTVWrN7nIyrI7ctsY+hrjCVycJyVeBuXNeBwZYNWGZxfpEyFK2hJMSEayoNvcemzwVXqeyvStlBxuDqxePN7LLSGcA7mEslm62aRgllGHZWYF23/x6dPn1XqcJaS0TVyuzqghhcr2uq+QbN1/aWfFxls1C737DLxOOF40BqFV+xvKoDlPZPXojWXRbEszthXnY+us71QssomF47QHHHFfJBAZlXaHwhjp5zwlAa0PCvL5UddAiXG1pjfcksqeXIXnERslENkV/ssO/YDaM40hQCWl1t0UJx1UKefTrEsU6YvX3QVq0A8gZEOgurPgKRMyASaRjHThRdjLRoBGqpPVjLiHrhevnA8HRdIWBdeh7ch7thNkahtffBOHpsOGAZpEwL3hOJKf/BAYC+qBUqlQZ0HSqTgYiEofl8EMHAmLpORqupgUqnoZSC8nmg1dTAGhgoGXAKjwcqb0AbyJYOPoSwV/vN5gCl7BlCHg8QCQL9cahE0k7als0WRo5p0QhEImPnh5lD+tI2WN294yalK0VraED/e09FskUicMJC6OdbOUWaaBYIXQeEnNtWVSGQ+NDFCB1MQm3bVbHdVnWA8tn/+Uvc+cv/gcXPAkcu9+Brf/wzvJJqwzl+e+ncX3Wfh8trX8fhXB1+jXdhxYa9iLgyOPZQaEwT+3ikz4v8skaIaQYoh64WaF7WBfFrD1BfA5nLFwUoQtOgmhuAKgxQei6oQb1hQcTi6D+nDq6UBU9PFr1n+OEZUPB15qB3dkFlR1xA62vRdXEUAJBcLLD8WC1QCFAkjr7Xh9ZnGtFzpgf1r7igt3cO/6ENXtBUPgex79C06pxf2gC9Jwnlc8N6ZQ+01iZAKeRawtBG13UyDbUQvf0QSsEKeSFFHUQuVzJAkX4/rFQKsn+g5Jo0QndB1tVADSZlE/W1sEI+pJaEEDjghhASwuO2A2VtMEJpqIPq7Z/1TMejZVY0wmtaUwpQ1KIGJD8aw8YzHsU3X38/wv+pVW/LIZGDCY8HQtNgzmmAItH/4QEkno6geVvldlvVAcp7fIfx/wXti1Y+pLDefwy1WgIXe+wLXns0gj8I7Mc+dxd+GXg3LowcRpMrhl9o55d/EClguSWmO5FTi+SwqqYTHULYCbpKrOiqPM6fJlqK4QOUy6674RUQSsLl1mD4BaQBuD0SLiFQdJ+s2c8DQD5kFSUtE1IgH1QwPRL5AGB6JPShGSqj+jIm6k6ZiOXRoFy6XW/LBHR75o7h1ab+Hrt0ey0cy27hUS571k1J2mDm13x+3Auz0jW7KxIalK7B8ugwvAKWW4emSbuuwHAyOZduD/Kd45YIw68N16VMllvD8ppOrPcfw79G+6c/mE5qkN5RswUsq/R5L/G4yhtQRh5y5PpXljWzbtz5NLJJfWjV68HPh9A0WOn0tD4fwuWGcOn29kPY4lUVhNuNqSwgJTwewCyxevVUjikF2mr6ccRf2ZQZVR2g/METn8HinfYfTcOLAheEPwuR1KEC9onWu9z4u9p1EFkNS17N4x9/uw6WS+H01GtlH0OZFly96Wl1JwBA6Bk/nt97NpbldkD2xcesu6MsBS2Vq8rBeDV7cpA9cRimiej+NGTOhIylUKsL6GkTen8a1qgPvYgNoPaNKADA161D9A0vHKhMEw07FHxHB1CvheE7OgAzP/heNtaPf/GfAvfxOEQiBSR0WF4vRDwJGAY8fvfU17WxLHtMRTYLrTsO5A2Y47VmKAVICZVKj7u+izBMKF2HCPiBbA5C+RF6KwU1ON5GDSSgRcLIeXVoQthjVea49QSwg1GVTE9ecASRt7DniVNxwYHPwr/Pg8XW9LK0ameswFsfrCtaj8k9AOT9gBoRM8k8oOWA/KjJXvW7TYS3HcPBjy0plHfHgdZ7dlRdkCJDIciAH0Z7BwBAXLgasRUB5IMCwgRSzQLLfnIAxomxq4xPJnXVeeg5Q8eyfz0MSDv4NvcfrPRLoFmQXrMCvuf2lF2+708uQPBYDvrj26d9TGWaOPrYUtTvrmyraFUHKKv+thta7AhMAHUPv4H6Z8LFd015o3CHbPX247RdYXs9k6mMNjZNyM6+aQcoLT/fA+Fywchmobq6x16cLBMiVV1fjEPcW9+w09ErBe2lvfadqGnBc6zdPuclonKjsxueLXbrh1fXYcSGAxQohZpH98BKpRB4ywMrmy2MCbEaayAHpnZRLMU6cLhwdyEjYVi9fXaAoeswp7psrTnYApRKwzxqdwGOexcy2MpiJRKl70SVPThW6DpU0A90dAMqAvH6QZjnr7QP1xeD3tYKw2e39ohcfl6WSTB8EmoK3TsAIHMGlt/XYS/ImM7AmOYg2fgZUfzgz38Elxje/tf95+NtwTfRoA9/lg7kGnEoW4/3hIpvRq7//Y1YeaIOf339vxbK/7LvIuz9f+Gyu32dQkbCMBuiwGCA0nV+ELHL0lja2ItU3oXbTnkc9/7+KmAaAcrRdQK3r/sP3Pf0lVACMIIuuBmgOJ8Q6LjQhaVPlxkoCIHEtXEMvBRB2+MzOK5SWPaTN6FS6YrebFd1gGK8dQQQ9syOMeu1lDCtaVCWNaO71JGDX8e9eBnV2H6CojEII7tcJhycZZkTvg9mf789g2L0hVfKmTcxC1FcN6kVmjbHfY8nGpVuWcPjYkyzuBtq9HZDwc94eVIsNbgQooASAsK0922lUhAKUGpoJpqw1yIC7ACpwqPmyyFMBZWf4p2SYQKd3TBmON3Y0gTa9DhcIzoOG91xtLl60CCHPzMZ3YWM5UKblijegdsCpCgq3+gawF5RhdmchQC0EYPHNcDtNhBxp6ELC62uvkkTB45HuSy06n1QUgASUFUyW4wANcWrutedR7ICuUKN9g4IjwdaOFyY3TNTVR2gzAVl2Hd8s3qMgcTkhU4GQkALhaBMc8wATNnVb8++me6uXW57ReHu7qILuvB6IJSyu15GzYYRHg9kKDj+DKv+uH2Z1HVojQ0QUsKKxWGlUtBCoaI/UpXJQpUYfzRcwLK7/3QdsrsPSpOQqRwsTYMcsD9/WjQK5A3InD1lXmVz0Bob5nxRu+je5JhzNaneGFSpNYimKLx/AFc88AUoOfweuvs0/Cj4B7Bcw49pGQmZA+4IF5/zhu0SruNd+B+//GyhvKdXw7LUzhnXba6pgQFoGF4Mse61DHqtEPaFw5AG8IlFn8LpncemNd6n6b813BT7nzjt2AkoKeByu6qyG/qkoxQat0+h61cpWI/Wo+lwZbpm8u84C72rPWi8a0tFbpwYoExCGcaszzYwZyHBTTUSmgYRDgH5PDAqQDFOdExp4NeYfXs9QE0YGJ2DxeWGtJT9Hoz6g5IeD1AbBXp6S/6xmT29EG43tJooVE0YlluHyOWAVAoiFARGBCiTjm9QCiqThfAJGL390Opq7a4/TYOMJe2cPtGwPQ4la58HlclAtLUAc73q7s49U85jMjownLZX9uH0vxnV2mGZY9cjGvqsjBpcrTIZmOk0TvvWiEHWlglzmoOu55MZiwPx4Zsb7fnX0PSSxx64rSx7te5p5nKpeehV1D7qHk4BQFXD//juMWP/JtLy/3bb6S8qcOzOCzyoX38M4p7KzNJjgOIAQju5plxKvx9WJlv6Ijdud8os3L+NDHiEBGAVH1/TYEb90Gprxk3aZeccEHYX1HjN6aO6YYSuT/p+C10b3p9hDnfvmBaEYQ3uVkC55/5PeEoLKhY2qkw3lMrn7KRQM1SJfcw7pQA1/HehstmxXaPTZA0MALxvqipaXS2EzwezY2qf7Ym6Y6TfD1k7vAqyFYtP2EUvLEApAa2lGTDshU2h6/ZMyRMdU/4eZ4Ay36QGraUZxpGj812TuSE15C9ZBc+e42NmFyjThBUfsMd2VJjKZCF6+ooulCqRtGfWmCa0mgjU6Bw1kRCO/EEIvvNOR/0/bBn7Ujwe4NQ2qFQWsCzIeBZW0r4Tt+LD+5F+v519cTAo0ZqbYBxvH/PHKkNBKKXs9XpqI1AuDerwMXtfpgkRDECl03Ad7YGhLIhQEJlGP9wu95wmZRrdfUVE8+/4/1iF2GoDZ3zXA3PfgYrsM3vpahz80PCNV9tvFsH30Avjlm/alsZxTyuOfN6CNARkHsg1GBBZiTO+hSknPGWAMs+EFLDqwsCR+a7J3BCahr7TPWjqioydXaDUrKznAAzeefcWD6Ie2e0iWpogMtmi41tBPzyX9KI/7kP9P4zdp/D5kFocgu/w4IyagWRhn0WBjs9blGfAqo9AdHRBjQxQhAB8Xns/Pg+MqA9KF5BDg4qFgDRMWJkMrMHsusrnQTaqwzOXa/EIYU+DZoBC5BxCwHxvPx4+/59wy/2fgdxXmd32rHbjufffYac1AHDZkS+i7aHxy2vPv4rF2dPx8f/3MPpNPzryEXy6dhveyAew6UcfARigVCGj8i0GTqZlMf7MpaGxBJXs0hlaD8ne8eBjcnh9JmUBpjVmEKswTQwkgrAmGuKu7HIib0Apq6zuG6WN0xVkWXZLiWlB5E0AWnHdh1p/Bv8VpgVpqPG7lmaDkFPORyO93jkZy0V00lIKyQEvtmWWQuQrdz3RMgrbso0jfp+kGvkctEQWL6WWIml40Jvz43l/Aw7kGiGmsSQHA5R5pkwTsn0B9IeXSRl51L3UD3SXGNMhNWg1ETvhWaXu0KUGLRiwBw5KAZVMAVLaa+Ok0xBee40cqz9mD84dqbMHzb+shZYdZ2qwYcAdywGdPfasFk2D1lA/NjGWUvbxhwgxdsDvYK4emCZEKg09k4UK+mG53XYXUS43ZhaT6umD70TUHgA8xbwk0yWkKLnK8oTbnLIEMpaY9npWRDS5RQ+68L2XPoTFh96afqbmURq3xfEXP/rzwu8tz5XxPdPZg9//+O0QloI0gL+IngEtB7R2vjHl4zNAmW9KVeVCgdOmFKxX3ig5aFJIAeHzAS6jYl0IQtPsrhPdTkOvcnkIXYfw++xgwOcDtFzJQbBmTy+CD20fN/MrLAtaIlvYVqurhQoHgROjy6niVgchSu6zqHsrHrcHvQ1llpVizDRdMx6H3hUvrPI8JzRt4unSJWSbQ/AAwLFZqRERAQj8ajsCQsKoYHev2vEaFr0y/P1STnoBs7sHTT9+cXADq9ACPJ21gZh+Z64IYa9voY+KCaUG6ffPT53my0QzOkZ2V8jhmSzC5R47lXS8XQzNrBlqqTBNWDVh5JfUQ7hd9lolQ39oug7our3/UlU1jPG7mwb3If1+u25Sg8jlh+swstyIP2yRzpecMi29diAlXG57fQxdhwyHhrdzF5+Dodc51YChLEPnb8zDYvAclvleeDxINbuQr/VDX9Ra6VoS0SBlGDMai6addir0luZRO1VQ+Vzhp9yu98I2g3Wabr0YoMwR6fNBa2qAjBbncNCCAYhli+epVs4mA377ogxAq6uxu2rK2S4UgvT5IH0+e2bQQAL950Rx6EofRCgI4XLBGkjYrRJ+L0TAD62uZvIdj2Ll8oAFiGWLoYWDEEE/VI89EFdGwsPl0pmixEniRGfpFqQliyCjEcjaKLTWZoiAH5kzFg0XqIsWnQNZUwPlc0PNwsBi6fEUL6g3xOWCigQhfd6y9qM1N6L9PSa6z/Oj6w+Wlh3YENEcEgJvXt+Irvctn++aFGGAMkeErkOF/HZz/UgeD/J15V14TzbC7bK7aAB7fRqPZ5ItBrfzeuxt3S77DiCbRWKRhPfMfiifncjKSmfsbKxeF5TXbe9/qgbvJvK1fgivF8rtKoydESMu7kN3EkPGy6li1gYgfD6IgB9WJADldSPVNNi8KgRUwFt0DkTAB+XWZ2WRO+F22+dv9OOahOV32+NeymCFArhw1UEkFiv0nz44hoWIHKf1kuOInTbftSjGAGWOKNOEyOSA3KiBmMqCluHsBmBwPZp8fvhibpqFnCgim7cT/5SzH8OEGjUrx9utMHAkDJHJQeXyha4fkc5BZHL2/qdBmCa0pF1nkTcgvXbLwuhVq0d2+UzYpTf0Gg0LIpODJ2ba+zIMiHSu+BzkDViu2WmRUKZZcr0dlTcgUzl7Ic4yiGwOL73VBl+ngJ4UkHW1la7qSUl4PNBqauyfcNhOqBXgjU410JoaJy80x/TmJhzuqIW301k3EBwkO0esdAbo7B47AyKXh9bZX7FR11XNMmH1x+yF8QBYiWRhMKnV01v2AnVWPD7mPDc824novhCszm47ABpaiO9EJ5RlTTwuZgIik4Ps7IEVT0ACkC1NsA4esmcFDZXRdchQyF7QUmoQi1uAvW+O3ZdhwYoPQHo9kIkUrK4eBHJ5O924aUJayv4cDb3Ovn4YZ7bCNQsLBlrpTMnWDpXNQhzvKFocciKqvQvL7q2Bu6sXPefXIL9qEeRcp+ZfgLSaKPKn2OMFZNqA1p+AcruAPfvnuWY0mdSFS+F5pGvOF/mcSOKiJah9XEP9C70VSXlfKQxQ5oo1dgE8YHBUdGJupohWg5HdFSO7Rcq9IAL2RXQ0c++bkHtR9MenDGPm4zfyRqHLRqUzUK0NY+sgJDCiu8QKlxjbAdjdUbm8/cWVzcFKJos+M6OnXlvpNEy3xKzM4bHMkksfKcOY0vos1sAA9M3bYQHwLr8EmTo3TrIh4bPD70O2zu5m01M6hGnC8s7hbC6atmSzDo+QRcsUzLdks4ba3SlYu6c+FXg2LbgARZ6zCrmGADyH+5BvjsDVm0KuMQDPoV6kV9TD0gS8//XynKYGtyumQTtlCdAXKx6DYJr2RakK6cuWwOrpGzf7q/B4SgYLs014PPYgTpcbZnc3hO4qvN9aXa097ReA2ddXkToKj51LRfj9QL7U+kIWxNBsJDl+a4ccyMBSFpDJFmbRjFynSfr9EKEgzK6ewvgXT+8cnV+pQXo9EB6PHSxaqqy/IaHrkJEwzJ5eeHrzyIfsWUrKyDvqDrLqpDNw9w1+pjMGZCIDkTe54rADCJcb2feeA9/OwyVXGtdyCjLgn7Ws2dPh7zKhxzOO+/wsuADl8FW1EBfH4H2kCT1vzyH4ah0yF6RQ87sWGNf1wu/OQ74QHHeg4myRPi+63tWMul1Be3XcQVY2O3ZcSpXoW9OK6EsuYJw/NBmNlPwDnW0yGoFqqoUR8UHb0gcZDMAc/DK3Tllkj2uxAPT3Q4bDM144TkYjULE4UBOGSI8NGJQ1IturpmG8NlR1rN1uoejrh9ZQD+nxQLjddsuJ1CAb6pA9pQGuLXFYGRNQCvqeIzDn4EIvfV7I+lqYDRHII51ALl84pxMRPh+s5a1ATy9ce49BrGiFDAdhxeLMLDsDZm8fXEPjnAzDnlGmcYaUE8hgAInPxqF9ZzG0Et9/7gELsjbqqAAltO1o0fphTrHgBsmm2gz8jxUvItUssOa0g0gtsnDFaa8i3SRw44pn8H9O/f3cJrYaJFw6Uk0C+fCofBtKzc5KvXMg3SCh/OPP5piP8zx0XDPoQS7qGtO9kg+5YYQ8MIP2+1BqpsqUj+d22wnMvK7xU/iPLD9OQDHUnVNo0ZEScNn3EEIKKK8H2RpX0YVorgJtoWlQXg+MkNt+vaPz+Yy3na4jH7Y/I1ZPL2TWsD8XYsF99cwplc3C7O6xf/pjsFIpR13wTmpuFz516tP232oJWlZBeUrnXZovxrHjjvz8LLgWlMgbOn7sfSca31J44cXTEN0r8Gvf+Wg5ZOFvt78PmsvEiuzcr8yn8gb0FJBucCE0C4Ma50PoqAEZT1V0UNXQLJgxU2eFKOqqmYjKZKD3pSAsZc/WGdGF4+5JDQaFgKUU1NBxBlPiTzXFvvJ5gP4BKMOANpAe3l9RIWs4C6ylxj1nWjQCM56AHAyalGkWWteUpSCSafg6snPS8jC6G0ZEwlBKQRjKboEqc7qwWtoCObQ2yHmrkIt64DbrgZ6+iTckqlbZLL65/QqEF2nwvf1ciOdeBmB/t8mmBgxENfjLHPB/sltwAcqiXx9D87NBaB0nUPtyCDKeQtOWIGRHL8Jv1kNJMS+RojIM+HosxJdKhMu80DpdcMtbsGLjXNClNq0gTEYj9iyeURd6obuKumomYvXHINIZaC4dpmHAGjkI+cDRwUL2RXOo/tLtAprq7e6qKdTbCvqAt45CGXlYJzpKBw8jVmlWRh7W8faxZQCgqQEinSkk81O53PDMJcuE2dUNPZGAmZv9z44MBuxZVIOfU7O5BjKWgjAse0p1mQsU9p4dQc0bCSip4fi7w/B1KqQbXIjsctJcAaLKMeMJrPw7A3tuMJBs9WP5c/bjIhJG4uwWJFskoqPW1aLSFlyAYhw8BBwEDCGAowqWEMBbyr5jHVzEbV7aLkwTekbB9KHsu0+nm3R8iTWNi5Cuo+TZkaLsboXRK+cWzQYaFZwWnpPFXUFl02WhpWGihGmF40xQTnld9oV/6HUqVTTSX2WzMOdq0LGuF31OTY8GqRSEaeePKfe9yIWEndNFCmTqFdxxwCz9DhMtDJYJtW0XtFvOR95l2akFpIBwu5ELS5gelJ3T6WS3IDuC9eYmYM3Z9myOM093RAIjZSn4j2UQ3WNNeTXYqqQsqPTU7xJUOl16O9Ms3X0ygr60DRAC8qxV0FYsLyREG/pXeDzQVp4CbeUphW20hsFpwaYJZHOFzLXlkn2Jwv61+rrhY3qLU8GPrMt4n0fZn7ATzA2dA6nZa/IAdhdUOAz9lGVTqt90qUwGsBS0mhpAanB1JSAGktB6E1DZ3KTvxZDIwbydo8NS8PQImG4BV0KNvwAj0QLh3+GDq8OFzs+sQfunL0H21EYEj2QR3W+OWfiTSluQAYrZ1ojj7wpCBgPovrgGsnbq66xUnGVC230ANf/15skxe0EpmPHElDcz+2Iw+8bm2VCGASs5cS6U9MpGCN2FnotqkDirobAeztCCezIYwMBZDYif01CYVaMWDQYouRxEMj3lmRDWkeOF/atFjfYxhYAYsRYPMLw2jwgFIWuiJfdlHu+AyucK50BoGuTg2jtC0yDqapBY3VhyEb9Ks5Ipe/xOY509QPbQUZjdPbCOHIeVTJX93vqf3w/zeDtgmQgdsWB6AG/vBAswEi0QbT8/jMh+4HM3/we+ePO/o2e1F66X30T4yf1Tyut0MltwXTwAYLk15IMKkBqMgAA0Z8RhViIBTHKRXVCmcxGaaJuJnhMCpk8DpEDeL6BnpB1sjOwykRoMn4AacX1XQ6nilZrWmBmVzRaOY3l0SE0DhBy7avVg4COkHPfzWBiXNPQ6XfqIugsol458QKK8ZfpmyDIBqUN57K6eQreUaRafJyHsGTnjvDdmLF54zpWyoKSEzLP1hBY+48hRuBOLcWXgIFwQ+LZ3bLJFmtiCDFD0zjjqdnmgMhnU7M1BOSRTq97aAqs+Auw/XDKr7IIiNUi3a1YWsisQAtLns1PW6y64BgzAUqjZm4UrkYdK2jN2lH9wfZxMBqGDaYwcAqF1x2HA7v4xFzcA3T1TrsbQcYb+L1z6uONZVDoNjDPIVXg89vpAg906I1vaxODML1dqbgaXykAAKpeH8tjjYkQgACgFsbgF5og0/fqiVhittcALu8buRAjIs1bCesXOTuk7moRrURimV0JfIDPZiCYSfCuJt23+34BUWPJGdea7mk8LMkCx3jqKSFcPzEQCni17YDqk1cJsrUP/6QHUdQQXfIAi3S57Rd9ZDFCEptkXzlQKwu2CqzcFyzTh2boXME2Yg/28VsgLSA1WIgHt5X32Y4MXR/PoCbu+Pi/iSwMIvDj17jdzxEwma2AA0u+H8pbOc2DG4uPmAJE+L1QuX+gyMrt7h2fLSHsbPWnOyYVdhoKw+mMw/Dp0TYMMBaFMC8nT6uDdd6BQh/ySenSf60fDthIBh5DoPyuK8CuDv+59C+7VZ9utQAtkJhvRROTuN3HG1+wuZaurx1Hr3FSDBRmgqHwOZv9g5lAHJZ+xdAnTJeZkDIEjzMVspRHHEHkT0usZ/z0fXCCw6KERF0mljapvuXf5g10YwlL2YNtxjj26/BhC2seU43RJKjWcU2S2DdXDso9r/y6g9MF6Ds4uUlLAcpV+n4UUsPTh56x0BtIscZ6JFigrlYJ18NB8V6NqOWNwxknCdaIPNW+k7C6BBU4ZBlS6sq0nMhQqPoalCt13KpeHyOaRe9sZYwJAraMfJVe+G7mvTBahfSMCG6lBL3NZ9KF62cdR9vIFvcUDfdWAPahUi4THncWj0mm7W2UgYf+YZmFgsMobQH8c7mNzk+BMJVMQHg9cJ/qh8gZUMgWVSiNwIF50Lt3H+1G/KzMmkJOhEGQohOjeES2FyoKlA6abAQrRVEmvd+zYtgVuVgKUY8eO4U//9E9RV1cHn8+Hs88+Gy+++GLheaUUbr/9drS0tMDn82HdunXYt2/fbFTFUcyjxyF3vAHTQa06s0UZRsVHqg91fRRYZuEYKp8D8ga6zvOMmSpsHG+ftCXEymSA3cOfQeHSYTWWN/tLBgOAEDCO2d1FKpu1F/QbYeg9F9EIRLB0gGJlMnbr38CAXd4y7YHVg6/P7OmFefhoWXWaKXNgAMLrgXXkeKFOViIB9fqBonNpHj4K15bXxmwvgwGIYADipT3DDyoF0yVgujFpwEhExUTAf9IFKBV/tX19fXjHO96B97znPXjkkUfQ0NCAffv2oaZm+Mv+jjvuwJ133omf/vSnWL58OW677TasX78er732Grzemc1R0MJhiFAQgD0oUni9UKk0VDYLEQpCCAGjo3N+BugJOZwZlMY3XtfKZF1jSkFp9t27yuUAKSE0CSudgcqahXT5kGLsCsZCTH9MRKkumXG7ccpoPRjqUin11FxNURf2LCghBBQG1+IxzaJVoYXXC2twHRhgMMfL0BRqjxtKCvt9GEGagMUGFKIpE6GgnaF7NiceOEzFA5Rvf/vbaGtrw09+8pPCY8uXLy/8XymF73//+/jLv/xLXHPNNQCAf/mXf0FTUxMefPBBfOQjH5nR8fs+sBrt77LvzsJ7dMRXGqjdqSH6Zg5H3+uGpQMrvjH3C2sJjwdyxTKIWAIqkYDZPzbXx0IivV476dg0LvpaJDz2/IwTtIwu6+lViF9+GrSMguUWSNdINGyPQb2yB9LrgWhrhRXyQux8o+hiL/1+qOzwOjcqb0B2TGEhPiGhN9bD7OsfG/yM2D/yBjBO15fweCCEPaVXBoOwkqmi9YGk3w8rnZmTHCJ6azPgdkGcuhR4bS+0lmaoWNyuixA4ev0qDJxmoO2RxQj816uw0hlkL12Ngx+0g7VlD1gQpoLnRGfRgHB/Rx75oDb+WB0iKqnj8kXw9Vnw/3Li1ddHfo9Uu4p38fzqV7/CRRddhD/+4z9GY2Mjzj//fPz4xz8uPH/w4EG0t7dj3bp1hccikQjWrFmDLVu2lNxnNptFPB4v+hlPxxrguSu/i+eu/C4867rwu6u+h9i7M+hZ7cGfXf0EvnD1ryAD/sq94DIJtxvptjCs+ghEMDjnx59rwuuZ9krBYgqZf4vKKgVfn4Wu8yS6z9HReb5E7ztzSCwP2amm/T5kF0cwcEpwTEI24fcXP2aZMDsn/iIYLmsHxKo2Yq/0W6qevsH+Y8Owg4xSZdzuwvYy4B+cpTT8WRV+vz2FeQ5Y9REojwvpxSEITYNZH7Hv4AbznoTe147nrvwuus7VC/XqOdON5678Lh6/4rsYWOxCssUN4StuEfV2pODpMzjFmGiKei+w0H3W5IkkhdsN4R1/lflqUvEA5cCBA7j77ruxcuVK/O53v8OnP/1p/O///b/x05/+FADQ3m6vh9PU1FS0XVNTU+G50TZt2oRIJFL4aWtrG/f4elpgR64eO3L1iCV82JlthZnQoWUV9iYb8Ua6Bcqch/5vy4IwBo87z1/OwuUeTqE+aweZwUer1Bo+4+1vZFnThJ6yoKUFtKz9WUBSh5YbLGNakIZVfqKwMl+Dsob3P9n6Q8q0Jhx/oYbyqZglPivKAqabIn4KM8e0hgbko167i0wNrqSs7FWMhaYBykJfwo8duXpAAiLgs+tlATty9Xg51wzTJ2D4MLbP3LLsfU2XnFqmX6I5JwT05UsLy1tUipaU0MrIkD+02OhCUPHbMcuycNFFF+Gb3/wmAOD888/H7t27cc899+D666+f1j43btyIW2+9tfB7PB4fN0hZ/HgeX+7/BACg7i0Lf/3iBrTtM+A/Eser4kxYukBzcue06jETKpuF91gcsFTZ65jMFq21Ccrrgbln/+wdRNeBaazFAwBWfGz3m5Ci5MV/ZFkrPgD/3i4sStdA5C0ol0Qu7IJ/TzdM0x5Q636rG652D8zRY4GkGHPx18LBsrri1EDCDh66e0t279iFFKBpUAMD43ZvjFyfQyWT9iyeEa9PJZLT6hoRug7h8ZSde6f7qhVINwgs/n0a7p40lLIgUllYiSRkKASzrw+R/wziyy9+AroFpFc0wHX4GMJHTHz5H+2/vXzTYKAVCQHtHcM7t+ypydOlRcIw++ZmJhPRdMhgEK/f0oLT/iUIbH+1Yvtd/LgBPTF5sreedy1Gza8qd9z5VPEApaWlBatXry567IwzzsB//ud/AgCam5sBAB0dHWhpaSmU6ejowHnnnVdynx6PB54y7/hd//USFj81+LJME2Eh7YXrTBMNr9hdDtY8JIhShgGc6ISoicKa54WirJoQjJAHcs/kZadLaBLWNFuqxr2QlrjzHll2KOeAfOuIXQcAPilgDqZnV9ksjKFZMKP2ZWdqtUY+YHcflRGgDA0SNXsmGLNiKUDKCYOEkeN1hvY5cibUdPuUha7bXS3lBijnKyifAfW4hIylYCoFkc7CTCSgt7YAfX2I/OJFRDQNHZ+8EMlWN6KWCf/xNAK/2QOhaXjjO2cBmoIVLP67FUoVLTUw5dcS8AMMUMjBhNeLz697BL94fD182yu3X/fvd5RxcIG+MwRqHlgYY7wq3sXzjne8A3v2FF/59u7di6VLlwKwB8w2Nzdj8+bNhefj8Ti2bt2KtWvXzrwClmkPdsxmIXw+qHzOngaqFFQ+N6/ZK5Vp2YMk53mAoMgZkNnZnQ2iDKPiU0lLzWDRm5vsbgQ5Yk0dyyz8KMOADAaHuxrGWXOnVLffVGfMCJd7/C4I05z8fRdiePtS+5lm94ZSyu5+KpO7X8LVq0NmDIjs4B1b3v63MIjYMKCyWbjjCq6kvW+ZzkNls7AyWbh7tMI+iuTy0NML48uTqCTLxObuVchGNOhtiwsPz3SKsL50MfSWpknLuWOlW5urUcUDlM9//vN4/vnn8c1vfhP79+/Hfffdh3/4h3/ATTfdBMC+U73lllvwjW98A7/61a+wa9cufOxjH0NrayuuvfbaitZFrVoGCIHshSscMX9cpdOw+vph5eZ5TYYTXdAOdUxebgaswURjlaJMc8yaSsLlRvf7ToGsq4U2uIpwybqcubywmvC4+x/d9aIUVGz8wdilaPW14w7AttKZMVNuRxO6q7C9FgyM+cxq4eC0ghSVn1pOmiW/T2HpIxmgvQtWr91aYcUHSp6Tuhd7EH65067/MftfWCaW/C6NpY/a+yjS0Q3XgdJjzYgWAiuRRPs/LEfvmcDRDy4pPD7Zd9BkDm5YhBNXL524kFJY9OTA/F9jKqTiV+2LL74YDzzwADZu3Iivf/3rWL58Ob7//e9jw4YNhTJf+tKXkEwmceONN6K/vx+XXnopHn300RnnQBkt3eSDV0gMLHajVtOAucohMQ5lGHOXx2ICc9GHP+5YjGnvUI3p4hCaRP9pQP0zfkBmilKwj5Rq8SF81Dfh7kt1n0y1S0WFAhD5PFBiBns5LXfCpQ/PfPJ4gGy26DMrvF6IRHLqDVOWaeeBKZN4dieEEDBHtDQVupxGnRPztb3D/x+x0KJ89mX7sVGtVWY8DnBFV1rAVDaLyH1b0X3VuYjX62gefFz4Jv4Omox+UR/6D4fRMNnxt5VYuLNKzUqzwgc+8AF84AMfGPd5IQS+/vWv4+tf//psHL7AlbC7Gbwxc/qzH8ixlFLwdgkglwfcLmjBgN1SMSoYcMcMqPzEdxTC5bbHKg0FBEJAq6stuuhORmRzMw9AR7QCDSVJG6IMA2o6n+PBqcFTyp8yOnW93z+1zMCcRkwnM6UgjvogPArmey4AAMi+DHD02LR3mTgShq99/BZUedYqGHU+6Dv2O2oNuplY0GvxeN44BiiF0EsnoIyF0eRFw1Quh9bNvbB6+2CFAsCiJmh1Y9PTe18/BmuSwa5aXU3RWj/S40HmvGVTqo/V3QsrOb2ZS4DdjVWY2myZYzLUWgOJaSVpkx4PpG8GrZNCQCxumbzcbDtZFtmkBWH5AymE90vot3XA+9UT6D1nZtN/T/15Dm2PjP899uaf1sD4y17g1PHTcFSbBR2gGIPTG41DR3hHtxApBWv3G7BSKSivDjPkBUpciI0T7ZN3Ofm8EJ4RSdY0DanmqSWas5LJmQ/CHjkFd9QFedrdZi7XuAnkymVF5j65YUn8O6YqIba8DH+XhV+f/is8fNojSDXNLMCW//0S1EvjTB8WApGzevCLVfchX1PZoRLzaf5HjhJVgBzIAJoERkzhFi435IqlUIeOTd49kckWdwMpBT09tYuh0HUIt3v6iySapt1dBXt15vme7QUAEMKeBZXMwgG1Iaoq3l4DHz90OVzSRODELM6sUQp9e2rxqcg10ONZLJQwngEKLQjquN1aZo2YKSPDQRx7XwMWP5iBdWjioMHs7Ssap6RyOfg6c+MvXFiC9PshaiKTHms8yjQLfccqnXbEejVC0yDra6E6yx+LM3uVkQtm+iSdHHyvHMGRb54GAKjbfRyzOUVixf1JdG9eDv+B/QvmZoIBCi0IpQaFCZcLqUUKqox1KUZ3nyjThJbOT+1OxKWXdazxK6GK8ow4gqZB+Tywjs//1GAxgwy0RPPBaO+A9+HBoQazfCy1bRc8wIIJToAFPgaFFjipTZjfRikF/zG7BWSytYdkiSnuMja1lhCVy0PMYJAsgFlZI2lMltypbOt2I98QgJygXjIQgL54UVFSqqnSmhqh1dVCW7F8/EJCDq97REQLHgMUqlrS65l4QS5LoenFFIRhQoYmXkF6zAJbSkEda5/SoEwrmYLZ1V12+bGV0Cat5/T2K6a9QKbwehFf6i1aVXnM7hvqkDqzBclzpj/Tx1jRCtXWjL5Lxs+UKTTpjHE5RDQn2MVD1UsIQJsgs6qyoHfG7RaUyTIJl3h+yoNdp5gQrRThmtrModkmNIm8HyXPzxDldiEX1mDpAtNt/8mH7NediUqMm2/T5eIsHqKTCFtQqoz0eiG9XgiPx04udhJThjFh+niVzUEMJKFS6UlXkFYjVl4Wum53N6w8ZUr1EboOrWZsHpayKauoHhUzg8GlKpNB5C1jwnqJgSQCxzIIHp3eNGhtsPUqF3bB229Bq68rcRCBgfWrZ5yNk4iqB1tQqoyIhAHThFAKyBsw53Hxw/mmcrkJAxQrmSpcWCfLwGr2DSdAEj4fVEsjBk6PILDvQNn1EW430Fg3/dV2lYIZT0xv24loctotD2Y8Ad/WiTNTmt290BJJOz3+dA7SWA8AyNRpCB7LQS1qBEZl8BW6C+0fzCK8iwEK0cmCAUqVEVLaM0ssBciTvD9+souuZZY/NnREhlYhBCyPBsM3xVkjQkC5ZvgnNY1MseXtd5qDSy1z0rWbVD5nT4mebt0Hz5npEpAZE6bfPTiwt/j9baqNT3tFZyKqPuziqTIqnQbSGahMxk7mdRLTwuGi2TfS6y10F0CIkjNzyqFyOWg9AwgfnOJigYYBsQDfExkITF7mnNOhNTVOa/8inkQ+ICENwPJoGFjmg754UVEZZZqIPdkM0b8w1hghoskxQKkyZiwOc2AAViIBazbGK1QRURuFGLF+jgj4IWqi9v81rei5qbCyWZjH2qHt3Del7VQuB5Gu8CrODiDDk5/HnnMjUM0lxo6UwerqRi4oIQ0Fw6chdqpEbmn9qEImlv70wMxmSRFRVWGAUmWE221n9/T7ITQNWjhsL+Y22SyVhUiIMWvXKG3wIy1k8XNToRRUPgcrmZzydgsy0+lEM6UGmR5A6dP7OlGGAQhAGgrCUlASMAJjP8/GiXbnJLAjolnHAKWKCF1H9rKzgfPPQMfHzoE4cyUO3XQWtFAI4qzTTrogRcUHoBLDQYRKpYH+uP1/0yx6bs7qNDALg1xnSkioGUzPVYnJX5Ovx4LWOfGK0ePu31JwpRSCh1LwHYmj4RUDfae5T7rPMxEVY4BSRYTbjY5L3Og9M4jwH51A/PQIPvrRxyEiYfSfET7pvtDN3r6iXCVWOg2zp3fwF3P6i/bNpE4TzHaZN1IUrTM0VWYsPmkZf3sW5olppsO3TLgSJrS9h2HtewvBp/YhvtI86T7PRFSMAUqVspTdfWGqEW+hlHa3x8lCqeKZHqV+nybhcnPGyJAyzqPM5Kff/SI1CAuAyw3h9UCl0xCGgKyvK+T8kaHQyfXZJiIGKNVE5XJoeiGH2lcTSDzYjPCeGH5+/2VQsThCh+wBs9o0B4bSMOHxQJ6yBHpr85S31SLj5kGdNzNZi6dcWvs0c78A0AfzoMTfdQpiV54J2dKElmcVDm1YCvPCVVDnnob4FWc68twS0exhG2oVUYYBz5O7ANNE02seqGwWSw74YA4MwHWgHRYAEQoC8cmb5Gl8wu1GdnEErl4PcOz4lFpiRDAI9E9vLMasEWLSRHUzZRw/Me1tVa09Nbz7HA35sIXwmyGEH3sdgc+E0dPdBj2r0H0eEN0Sct65JaJZwwClyqisPY11qDldDY55UPk8ZG0NVMBnN4VzzZIZEaaCsCxM+SyepIvZCd0FNd2sxsqeuQMAwgSMoBsupdCb9sMICFhuQOlq+rOyiKgqsYtngVDJFDr/YAn6LqyfeIVfmpTKZOE52AXR3jPlQM+KO3CQ7GwTAtqiqXeHFXT3IxfS4OsEQgck+k7zQCxphfiPOiSWWug/04LMMDghOtkwQFkgrEwGvWcrxE6VEL7pZVAlm8rnYBw+CrOjc8rbTjl3ykIgJMz66Y8PsfpjyPsFvL0WwocNJBcDucYA6n7xCtCSQf3yXsg8OEiW6CTDLp4FxB2TcMUBmAswWdgcEroO4fNBaBLmAhjzoGY7eZyyIHLG1LvDRmzvTljQMwqupAFXXIfp0aBns7C6PejOaQj2CSDPJG1EJxMGKAtI23+lIDP5ecn/sZDImhoYK1ph+nToj2+f7+rMmIrFZ28RQgBQCuLY1FubCpubJiIvd0MYJpDLoyVTj8QSH9ymiWW/MmB5JLwdcVgLIFgkovIxQFlAxHMvzyhjKNlEwIdUqxe5oETNfFemAqzM1BY9nI5CgrzpUArm3jcLv8qubuRXXQAoBdfmHYNF1PRbaIioKjFAWUgYnFRG3oB7wIQ0eT7nhWlCyyvORiM6yTFAIRrF6uuH71UALh0c9TD3lGnClVQzm7pMRFWPAQrRKFYqxXE880kpaBkLQpNQ+fmuDBHNF04zJipBhkIQHs98V+Ok5e7PcTwV0UmOAQpRCWJxM7TahTBEtjrpe49A5di9Q3QyY4BCVIIV9AJetqDMF7OnlwNkiU5yDFCISpCxFJCe/em5RERUGgfJEpWgjp6AxcylRETzhgEKUQmcxUNENL/YxUNERESOwwCFiIiIHIcBChERETkOAxQiIiJyHAYoRERE5DgMUIiIiMhxGKAQERGR4zBAISIiIsdhgEJERESOwwCFiIiIHIcBChERETkOAxQiIiJyHAYoRERE5DgVD1BM08Rtt92G5cuXw+fz4dRTT8Vf//VfQylVKKOUwu23346Wlhb4fD6sW7cO+/btq3RViIiIqEpVPED59re/jbvvvht///d/j9dffx3f/va3cccdd+AHP/hBocwdd9yBO++8E/fccw+2bt2KQCCA9evXI5PJVLo6REREVIX0Su/wueeewzXXXIOrrroKALBs2TL87Gc/wwsvvADAbj35/ve/j7/8y7/ENddcAwD4l3/5FzQ1NeHBBx/ERz7ykUpXiYiIiKpMxVtQ3v72t2Pz5s3Yu3cvAODll1/GM888gyuuuAIAcPDgQbS3t2PdunWFbSKRCNasWYMtW7aU3Gc2m0U8Hi/6ISIiooWr4i0oX/nKVxCPx7Fq1SpomgbTNPE3f/M32LBhAwCgvb0dANDU1FS0XVNTU+G50TZt2oSvfe1rla4qEREROVTFW1B+/vOf49/+7d9w3333YceOHfjpT3+K73znO/jpT3867X1u3LgRsVis8HPkyJEK1piIiIicpuItKF/84hfxla98pTCW5Oyzz8ahQ4ewadMmXH/99WhubgYAdHR0oKWlpbBdR0cHzjvvvJL79Hg88Hg8la4qEREROVTFW1BSqRSkLN6tpmmwLAsAsHz5cjQ3N2Pz5s2F5+PxOLZu3Yq1a9dWujpERERUhSregnL11Vfjb/7mb7BkyRKceeaZeOmll/Dd734Xn/jEJwAAQgjccsst+MY3voGVK1di+fLluO2229Da2oprr7220tUhIiKiKlTxAOUHP/gBbrvtNnzmM59BZ2cnWltb8b/+1//C7bffXijzpS99CclkEjfeeCP6+/tx6aWX4tFHH4XX6610dYiIiKgKCTUyxWuViMfjiEQiuAzXQBeu+a4OERERlcFQeTyJhxCLxRAOhycsy7V4iIiIyHEYoBAREZHjMEAhIiIix2GAQkRERI7DAIWIiIgchwEKEREROQ4DFCIiInIcBihERETkOAxQiIiIyHEYoBAREZHjMEAhIiIix2GAQkRERI7DAIWIiIgchwEKEREROQ4DFCIiInIcBihERETkOAxQiIiIyHEYoBAREZHjMEAhIiIix2GAQkRERI7DAIWIiIgchwEKEREROQ4DFCIiInIcBihERETkOAxQiIiIyHEYoBAREZHjMEAhIiIix2GAQkRERI7DAIWIiIgchwEKEREROQ4DFCIiInIcBihERETkOAxQiIiIyHEYoBAREZHjMEAhIiIix2GAQkRERI7DAIWIiIgchwEKEREROQ4DFCIiInIcBihERETkOAxQiIiIyHEYoBAREZHjMEAhIiIix2GAQkRERI7DAIWIiIgchwEKEREROQ4DFCIiInKcKQcoTz/9NK6++mq0trZCCIEHH3yw6HmlFG6//Xa0tLTA5/Nh3bp12LdvX1GZ3t5ebNiwAeFwGNFoFDfccAMSicSMXggREREtHFMOUJLJJM4991zcddddJZ+/4447cOedd+Kee+7B1q1bEQgEsH79emQymUKZDRs24NVXX8Vjjz2Ghx9+GE8//TRuvPHG6b8KIiIiWlCEUkpNe2Mh8MADD+Daa68FYLeetLa24gtf+AL+z//5PwCAWCyGpqYm3HvvvfjIRz6C119/HatXr8a2bdtw0UUXAQAeffRRXHnllTh69ChaW1snPW48HkckEsFluAa6cE23+kRERDSHDJXHk3gIsVgM4XB4wrIVHYNy8OBBtLe3Y926dYXHIpEI1qxZgy1btgAAtmzZgmg0WghOAGDdunWQUmLr1q0l95vNZhGPx4t+iIiIaOGqaIDS3t4OAGhqaip6vKmpqfBce3s7Ghsbi57XdR21tbWFMqNt2rQJkUik8NPW1lbJahMREZHDVMUsno0bNyIWixV+jhw5Mt9VIiIiollU0QClubkZANDR0VH0eEdHR+G55uZmdHZ2Fj1vGAZ6e3sLZUbzeDwIh8NFP0RERLRwVTRAWb58OZqbm7F58+bCY/F4HFu3bsXatWsBAGvXrkV/fz+2b99eKPP444/DsiysWbOmktUhIiKiKqVPdYNEIoH9+/cXfj948CB27tyJ2tpaLFmyBLfccgu+8Y1vYOXKlVi+fDluu+02tLa2Fmb6nHHGGXj/+9+PT37yk7jnnnuQz+dx88034yMf+UhZM3iIiIho4ZtygPLiiy/iPe95T+H3W2+9FQBw/fXX495778WXvvQlJJNJ3Hjjjejv78ell16KRx99FF6vt7DNv/3bv+Hmm2/G5ZdfDiklrrvuOtx5550VeDlERES0EMwoD8p8YR4UIiKi6jNveVCIiIiIKoEBChERETkOAxQiIiJyHAYoRERE5DgMUIiIiMhxGKAQERGR4zBAISIiIsdhgEJERESOwwCFiIiIHIcBChERETkOAxQiIiJyHAYoRERE5DgMUIiIiMhxGKAQERGR4zBAISIiIsdhgEJERESOwwCFiIiIHIcBChERETkOAxQiIiJyHAYoRERE5DgMUIiIiMhxGKAQERGR4zBAISIiIsdhgEJERESOwwCFiIiIHIcBChERETkOAxQiIiJyHAYoRERE5DgMUIiIiMhxGKAQERGR4zBAISIiIsdhgEJERESOwwCFiIiIHIcBChERETkOAxQiIiJyHAYoRERE5DgMUIiIiMhxGKAQERGR4zBAISIiIsdhgEJERESOwwCFiIiIHIcBChERETkOAxQiIiJyHAYoRERE5DgMUIiIiMhxGKAQERGR40w5QHn66adx9dVXo7W1FUIIPPjgg4Xn8vk8vvzlL+Pss89GIBBAa2srPvaxj+H48eNF++jt7cWGDRsQDocRjUZxww03IJFIzPjFEBER0cIw5QAlmUzi3HPPxV133TXmuVQqhR07duC2227Djh078Mtf/hJ79uzBH/7hHxaV27BhA1599VU89thjePjhh/H000/jxhtvnP6rICIiogVFKKXUtDcWAg888ACuvfbaccts27YNl1xyCQ4dOoQlS5bg9ddfx+rVq7Ft2zZcdNFFAIBHH30UV155JY4ePYrW1tZJjxuPxxGJRHAZroEuXNOtPhEREc0hQ+XxJB5CLBZDOByesOysj0GJxWIQQiAajQIAtmzZgmg0WghOAGDdunWQUmLr1q2zXR0iIiKqAvps7jyTyeDLX/4yPvrRjxYipfb2djQ2NhZXQtdRW1uL9vb2kvvJZrPIZrOF3+Px+OxVmoiIiObdrLWg5PN5/Mmf/AmUUrj77rtntK9NmzYhEokUftra2ipUSyIiInKiWQlQhoKTQ4cO4bHHHivqZ2pubkZnZ2dRecMw0Nvbi+bm5pL727hxI2KxWOHnyJEjs1FtIiIicoiKd/EMBSf79u3DE088gbq6uqLn165di/7+fmzfvh0XXnghAODxxx+HZVlYs2ZNyX16PB54PJ5KV5WIiIgcasoBSiKRwP79+wu/Hzx4EDt37kRtbS1aWlrwoQ99CDt27MDDDz8M0zQL40pqa2vhdrtxxhln4P3vfz8++clP4p577kE+n8fNN9+Mj3zkI2XN4CEiIqKFb8rTjJ988km85z3vGfP49ddfj7/6q7/C8uXLS273xBNP4LLLLgNgJ2q7+eab8etf/xpSSlx33XW48847EQwGy6oDpxkTERFVn6lMM55RHpT5wgCFiIio+jgqDwoRERHRVDFAISIiIsdhgEJERESOwwCFiIiIHIcBChERETkOAxQiIiJyHAYoRERE5DgMUIiIiMhxGKAQERGR4zBAISIiIsdhgEJERESOwwCFiIiIHIcBChERETkOAxQiIiJyHH2+KzAdSikAgIE8oOa5MkRERFQWA3kAw9fxiVRlgNLT0wMAeAa/neeaEBER0VQNDAwgEolMWKYqA5Ta2loAwOHDhyd9gTSxeDyOtrY2HDlyBOFweL6rU9V4LiuL57NyeC4rh+dyZpRSGBgYQGtr66RlqzJAkdIeOhOJRPgBqZBwOMxzWSE8l5XF81k5PJeVw3M5feU2LHCQLBERETkOAxQiIiJynKoMUDweD7761a/C4/HMd1WqHs9l5fBcVhbPZ+XwXFYOz+XcEaqcuT5EREREc6gqW1CIiIhoYWOAQkRERI7DAIWIiIgchwEKEREROU5VBih33XUXli1bBq/XizVr1uCFF16Y7yrNq6effhpXX301WltbIYTAgw8+WPS8Ugq33347Wlpa4PP5sG7dOuzbt6+oTG9vLzZs2IBwOIxoNIobbrgBiUSiqMwrr7yCd77znfB6vWhra8Mdd9wx2y9tzm3atAkXX3wxQqEQGhsbce2112LPnj1FZTKZDG666SbU1dUhGAziuuuuQ0dHR1GZw4cP46qrroLf70djYyO++MUvwjCMojJPPvkkLrjgAng8HqxYsQL33nvvbL+8OXX33XfjnHPOKSS0Wrt2LR555JHC8zyP0/etb30LQgjccssthcd4PsvzV3/1VxBCFP2sWrWq8DzPo4OoKnP//fcrt9ut/vmf/1m9+uqr6pOf/KSKRqOqo6Njvqs2b37729+qv/iLv1C//OUvFQD1wAMPFD3/rW99S0UiEfXggw+ql19+Wf3hH/6hWr58uUqn04Uy73//+9W5556rnn/+efXf//3fasWKFeqjH/1o4flYLKaamprUhg0b1O7du9XPfvYz5fP51I9+9KO5eplzYv369eonP/mJ2r17t9q5c6e68sor1ZIlS1QikSiU+dSnPqXa2trU5s2b1Ysvvqje9ra3qbe//e2F5w3DUGeddZZat26deumll9Rvf/tbVV9frzZu3Fgoc+DAAeX3+9Wtt96qXnvtNfWDH/xAaZqmHn300Tl9vbPpV7/6lfrNb36j9u7dq/bs2aP+7//9v8rlcqndu3crpXgep+uFF15Qy5YtU+ecc4763Oc+V3ic57M8X/3qV9WZZ56pTpw4Ufjp6uoqPM/z6BxVF6Bccskl6qabbir8bpqmam1tVZs2bZrHWjnH6ADFsizV3Nys/vZv/7bwWH9/v/J4POpnP/uZUkqp1157TQFQ27ZtK5R55JFHlBBCHTt2TCml1A9/+ENVU1OjstlsocyXv/xldfrpp8/yK5pfnZ2dCoB66qmnlFL2uXO5XOoXv/hFoczrr7+uAKgtW7YopeyAUUqp2tvbC2XuvvtuFQ6HC+fvS1/6kjrzzDOLjvXhD39YrV+/frZf0ryqqalR//iP/8jzOE0DAwNq5cqV6rHHHlPvfve7CwEKz2f5vvrVr6pzzz235HM8j85SVV08uVwO27dvx7p16wqPSSmxbt06bNmyZR5r5lwHDx5Ee3t70TmLRCJYs2ZN4Zxt2bIF0WgUF110UaHMunXrIKXE1q1bC2Xe9a53we12F8qsX78ee/bsQV9f3xy9mrkXi8UADC9QuX37duTz+aLzuWrVKixZsqTofJ599tloamoqlFm/fj3i8TheffXVQpmR+xgqs1A/x6Zp4v7770cymcTatWt5HqfppptuwlVXXTXmNfN8Ts2+ffvQ2tqKU045BRs2bMDhw4cB8Dw6TVUFKN3d3TBNs+iDAQBNTU1ob2+fp1o529B5meictbe3o7Gxseh5XddRW1tbVKbUPkYeY6GxLAu33HIL3vGOd+Css84CYL9Wt9uNaDRaVHb0+ZzsXI1XJh6PI51Oz8bLmRe7du1CMBiEx+PBpz71KTzwwANYvXo1z+M03H///dixYwc2bdo05jmez/KtWbMG9957Lx599FHcfffdOHjwIN75zndiYGCA59FhqnI1Y6K5cNNNN2H37t145pln5rsqVev000/Hzp07EYvF8B//8R+4/vrr8dRTT813tarOkSNH8LnPfQ6PPfYYvF7vfFenql1xxRWF/59zzjlYs2YNli5dip///Ofw+XzzWDMarapaUOrr66Fp2pgR1R0dHWhubp6nWjnb0HmZ6Jw1Nzejs7Oz6HnDMNDb21tUptQ+Rh5jIbn55pvx8MMP44knnsDixYsLjzc3NyOXy6G/v7+o/OjzOdm5Gq9MOBxeUF+SbrcbK1aswIUXXohNmzbh3HPPxd/93d/xPE7R9u3b0dnZiQsuuAC6rkPXdTz11FO48847oes6mpqaeD6nKRqN4rTTTsP+/fv5uXSYqgpQ3G43LrzwQmzevLnwmGVZ2Lx5M9auXTuPNXOu5cuXo7m5ueicxeNxbN26tXDO1q5di/7+fmzfvr1Q5vHHH4dlWVizZk2hzNNPP418Pl8o89hjj+H0009HTU3NHL2a2aeUws0334wHHngAjz/+OJYvX170/IUXXgiXy1V0Pvfs2YPDhw8Xnc9du3YVBX2PPfYYwuEwVq9eXSgzch9DZRb659iyLGSzWZ7HKbr88suxa9cu7Ny5s/Bz0UUXYcOGDYX/83xOTyKRwJtvvomWlhZ+Lp1mvkfpTtX999+vPB6Puvfee9Vrr72mbrzxRhWNRotGVJ9sBgYG1EsvvaReeuklBUB997vfVS+99JI6dOiQUsqeZhyNRtVDDz2kXnnlFXXNNdeUnGZ8/vnnq61bt6pnnnlGrVy5smiacX9/v2pqalJ/9md/pnbv3q3uv/9+5ff7F9w0409/+tMqEomoJ598smgaYiqVKpT51Kc+pZYsWaIef/xx9eKLL6q1a9eqtWvXFp4fmob4vve9T+3cuVM9+uijqqGhoeQ0xC9+8Yvq9ddfV3fdddeCm4b4la98RT311FPq4MGD6pVXXlFf+cpXlBBC/f73v1dK8TzO1MhZPErxfJbrC1/4gnryySfVwYMH1bPPPqvWrVun6uvrVWdnp1KK59FJqi5AUUqpH/zgB2rJkiXK7XarSy65RD3//PPzXaV59cQTTygAY36uv/56pZQ91fi2225TTU1NyuPxqMsvv1zt2bOnaB89PT3qox/9qAoGgyocDquPf/zjamBgoKjMyy+/rC699FLl8XjUokWL1Le+9a25eolzptR5BKB+8pOfFMqk02n1mc98RtXU1Ci/36/+6I/+SJ04caJoP2+99Za64oorlM/nU/X19eoLX/iCyufzRWWeeOIJdd555ym3261OOeWUomMsBJ/4xCfU0qVLldvtVg0NDeryyy8vBCdK8TzO1OgAheezPB/+8IdVS0uLcrvdatGiRerDH/6w2r9/f+F5nkfnEEopNT9tN0RERESlVdUYFCIiIjo5MEAhIiIix2GAQkRERI7DAIWIiIgchwEKEREROQ4DFCIiInIcBihERETkOAxQiIiIyHEYoBAREZHjMEAhIiIix2GAQkRERI7DAIWIiIgc5/8Hh4op5Th6M/EAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "from matplotlib import pyplot as plt\n",
    "\n",
    "frame_pianoroll = score.tracks[10].pianoroll(960 / 4, 'frame')\n",
    "print(\"frame_pianoroll\", frame_pianoroll.dtype, frame_pianoroll.shape)\n",
    "onset_pianoroll = score.tracks[10].pianoroll(960 / 4, 'onset')\n",
    "print(\"onset_pianoroll\", onset_pianoroll.dtype, onset_pianoroll.shape)\n",
    "pianoroll = score.tracks[10].pianoroll(960 / 4, 'both')\n",
    "print(\"both_pianoroll\", pianoroll.dtype, pianoroll.shape)\n",
    "plt.imshow(frame_pianoroll[0], aspect='auto')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "id": "TunMY7rGLUFQ"
   },
   "source": [
    "# Benchmark\n",
    "\n",
    "## MIDI Parsing\n",
    "\n",
    "* mido is writen in pure python, and only parsing midi file to event level\n",
    "* pretty_midi and miditoolkit is based on mido"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "ihVpllc6LUFQ",
    "outputId": "30e3d4c7-0975-422b-d984-5f8378445a3b"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Requirement already satisfied: jill in /usr/local/lib/python3.10/dist-packages (0.11.5)\n",
      "Requirement already satisfied: fire<0.6.0,>=0.5.0 in /usr/local/lib/python3.10/dist-packages (from jill) (0.5.0)\n",
      "Requirement already satisfied: jsonschema<5.0.0,>=4.17.3 in /usr/local/lib/python3.10/dist-packages (from jill) (4.19.2)\n",
      "Requirement already satisfied: python-gnupg<0.6.0,>=0.5.0 in /usr/local/lib/python3.10/dist-packages (from jill) (0.5.1)\n",
      "Requirement already satisfied: requests<3.0.0,>=2.28.2 in /usr/local/lib/python3.10/dist-packages (from jill) (2.31.0)\n",
      "Requirement already satisfied: requests-futures<2.0.0,>=1.0.0 in /usr/local/lib/python3.10/dist-packages (from jill) (1.0.1)\n",
      "Requirement already satisfied: semantic-version<3.0.0,>=2.10.0 in /usr/local/lib/python3.10/dist-packages (from jill) (2.10.0)\n",
      "Requirement already satisfied: wget<4.0,>=3.2 in /usr/local/lib/python3.10/dist-packages (from jill) (3.2)\n",
      "Requirement already satisfied: six in /usr/local/lib/python3.10/dist-packages (from fire<0.6.0,>=0.5.0->jill) (1.16.0)\n",
      "Requirement already satisfied: termcolor in /usr/local/lib/python3.10/dist-packages (from fire<0.6.0,>=0.5.0->jill) (2.3.0)\n",
      "Requirement already satisfied: attrs>=22.2.0 in /usr/local/lib/python3.10/dist-packages (from jsonschema<5.0.0,>=4.17.3->jill) (23.1.0)\n",
      "Requirement already satisfied: jsonschema-specifications>=2023.03.6 in /usr/local/lib/python3.10/dist-packages (from jsonschema<5.0.0,>=4.17.3->jill) (2023.11.1)\n",
      "Requirement already satisfied: referencing>=0.28.4 in /usr/local/lib/python3.10/dist-packages (from jsonschema<5.0.0,>=4.17.3->jill) (0.31.0)\n",
      "Requirement already satisfied: rpds-py>=0.7.1 in /usr/local/lib/python3.10/dist-packages (from jsonschema<5.0.0,>=4.17.3->jill) (0.13.0)\n",
      "Requirement already satisfied: charset-normalizer<4,>=2 in /usr/local/lib/python3.10/dist-packages (from requests<3.0.0,>=2.28.2->jill) (3.3.2)\n",
      "Requirement already satisfied: idna<4,>=2.5 in /usr/local/lib/python3.10/dist-packages (from requests<3.0.0,>=2.28.2->jill) (3.4)\n",
      "Requirement already satisfied: urllib3<3,>=1.21.1 in /usr/local/lib/python3.10/dist-packages (from requests<3.0.0,>=2.28.2->jill) (2.0.7)\n",
      "Requirement already satisfied: certifi>=2017.4.17 in /usr/local/lib/python3.10/dist-packages (from requests<3.0.0,>=2.28.2->jill) (2023.7.22)\n",
      "\u001b[1mJILL - Julia Installer 4 Linux (MacOS, Windows and FreeBSD) -- Light\u001b[0m\n",
      "\n",
      "\u001b[92mquerying release information from https://julialang-s3.julialang.org/bin/versions.json\u001b[0m\n",
      "julia 1.9.4 already installed.\n",
      "True\n"
     ]
    }
   ],
   "source": [
    "# install julia for testing MIDI.jl\n",
    "!pip install jill\n",
    "!jill install -c True"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "id": "KKKW8JljLUFR"
   },
   "outputs": [],
   "source": [
    "import mido, music21\n",
    "import pretty_midi as pm\n",
    "import miditoolkit as mtk\n",
    "p = \"mahler.mid\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {
    "id": "AwG9h0wPMsSh"
   },
   "outputs": [],
   "source": [
    "# %%timeit\n",
    "# midi_jl.load(p)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "BZ7zTxdGLUFR",
    "outputId": "e451a265-42fd-47fc-b48d-5c83a94dc14d"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "23.3 ms ± 12.3 ms per loop (mean ± std. dev. of 7 runs, 100 loops each)\n"
     ]
    }
   ],
   "source": [
    "%%timeit\n",
    "Score(p)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "PFPvSmJBLUFS",
    "outputId": "534187bb-c8c3-4a61-bf57-56fb34be2c86"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "7.15 s ± 2.11 s per loop (mean ± std. dev. of 7 runs, 1 loop each)\n"
     ]
    }
   ],
   "source": [
    "%%timeit\n",
    "# mido is writen in pure python, and only parsing midi file to event level\n",
    "mido.MidiFile(p)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "wGJ0OnE6LUFS",
    "outputId": "0479eff7-3dc4-479d-f02b-f09734e7f928"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "7.1 s ± 1.64 s per loop (mean ± std. dev. of 7 runs, 1 loop each)\n"
     ]
    }
   ],
   "source": [
    "%%timeit\n",
    "pm.PrettyMIDI(p)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "kTF9cTUBLUFT",
    "outputId": "6fbf0963-1329-47d7-ccc9-227b8914b9f1"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "6.73 s ± 1.54 s per loop (mean ± std. dev. of 7 runs, 1 loop each)\n"
     ]
    }
   ],
   "source": [
    "%%timeit\n",
    "mtk.MidiFile(p)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "4aRmyoYPLUFT",
    "outputId": "effd7541-df89-4c87-8618-325d2073664c"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "10.2 s ± 558 ms per loop (mean ± std. dev. of 7 runs, 1 loop each)\n"
     ]
    }
   ],
   "source": [
    "%%timeit\n",
    "music21.converter.parse(p)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {
    "id": "f2-NDVlELUFT"
   },
   "outputs": [],
   "source": [
    "jl_script = f\"\"\"\n",
    "import Pkg\n",
    "Pkg.add(\"MIDI\")\n",
    "Pkg.add(\"BenchmarkTools\")\n",
    "using BenchmarkTools\n",
    "using MIDI\n",
    "b = @benchmark load(\"{p}\")\n",
    "println(\"first run: \", b)\n",
    "b = @benchmark load(\"{p}\")\n",
    "println(\"scecond run: \", b)\n",
    "\"\"\"\n",
    "with open(\"bench.jl\", \"w\") as f:\n",
    "  f.write(jl_script)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "NJPsOcRDSxdg",
    "outputId": "34ba1bfc-320c-4b33-eb58-6a67a69fd31b"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[32m\u001b[1m   Resolving\u001b[22m\u001b[39m package versions...\n",
      "\u001b[32m\u001b[1m  No Changes\u001b[22m\u001b[39m to `~/.julia/environments/v1.9/Project.toml`\n",
      "\u001b[32m\u001b[1m  No Changes\u001b[22m\u001b[39m to `~/.julia/environments/v1.9/Manifest.toml`\n",
      "\u001b[?25l\u001b[?25h\u001b[2K\u001b[?25h\u001b[32m\u001b[1m   Resolving\u001b[22m\u001b[39m package versions...\n",
      "\u001b[32m\u001b[1m  No Changes\u001b[22m\u001b[39m to `~/.julia/environments/v1.9/Project.toml`\n",
      "\u001b[32m\u001b[1m  No Changes\u001b[22m\u001b[39m to `~/.julia/environments/v1.9/Manifest.toml`\n",
      "\u001b[?25l\u001b[?25h\u001b[2K\u001b[?25hfirst run: Trial(132.627 ms)\n",
      "scecond run: Trial(132.829 ms)\n"
     ]
    }
   ],
   "source": [
    "!julia bench.jl"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "id": "P4WwHyBwd2JC"
   },
   "source": [
    "## Bench Processing\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {
    "id": "nl7QrFqqeA0K"
   },
   "outputs": [],
   "source": [
    "score = Score(p)\n",
    "score2 = score.copy()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "9kkPPSLYekkt",
    "outputId": "fdba9ddb-a5f2-454a-e5fc-7423c797599e"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "276 µs ± 16.3 µs per loop (mean ± std. dev. of 7 runs, 1000 loops each)\n"
     ]
    }
   ],
   "source": [
    "%%timeit\n",
    "score.copy()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "SX7QdfpQd1eW",
    "outputId": "012d0540-fe24-4ceb-972c-8f1ac9707817"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "364 µs ± 6.46 µs per loop (mean ± std. dev. of 7 runs, 1000 loops each)\n"
     ]
    }
   ],
   "source": [
    "%%timeit\n",
    "score.shift_pitch(10)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "w3mQlDdkeNbK",
    "outputId": "36028a05-3bcd-4473-e3a4-58e785f369bc"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "10.3 ms ± 389 µs per loop (mean ± std. dev. of 7 runs, 100 loops each)\n"
     ]
    }
   ],
   "source": [
    "%%timeit\n",
    "[(t.onset_pianoroll(960 / 4), t.frame_pianoroll(960 / 4)) for t in score.tracks]"
   ]
  }
 ],
 "metadata": {
  "colab": {
   "provenance": []
  },
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.11.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}