tsp2.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The Cython extension is already loaded. To reload it, use:\n",
      "  %reload_ext Cython\n"
     ]
    }
   ],
   "source": [
    "%load_ext Cython"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "%%cython\n",
    "\n",
    "import random\n",
    "import operator\n",
    "import math\n",
    "\n",
    "xsize=10\n",
    "ysize=10\n",
    "\n",
    "def makeCities(n=50):\n",
    "    cities=[]\n",
    "    for i in range(0,n):\n",
    "        cities.append((i,xsize*random.random(),ysize*random.random()))\n",
    "    return cities\n",
    "    \n",
    "def connectCities(cities):\n",
    "    return list(range(0,len(cities)))\n",
    "    \n",
    "def pathLength(cities, edges):\n",
    "    length=0\n",
    "    for c in range(0,len(edges)):\n",
    "        i=edges[c]\n",
    "        j=edges[(c+1)%len(edges)]\n",
    "        x1=cities[i][1]\n",
    "        x2=cities[j][1]\n",
    "        y1=cities[i][2]\n",
    "        y2=cities[j][2]\n",
    "        length+=math.sqrt((x1-x2)**2+(y1-y2)**2)\n",
    "    return length\n",
    "    \n",
    "def maxEdgeLength(cities, edges):\n",
    "    m = 0\n",
    "    for c in range(0,len(edges)):\n",
    "        i=edges[c]\n",
    "        j=edges[(c+1)%len(edges)]\n",
    "        x1=cities[i][1]\n",
    "        x2=cities[j][1]\n",
    "        y1=cities[i][2]\n",
    "        y2=cities[j][2]\n",
    "        m = max(m, math.sqrt((x1-x2)**2+(y1-y2)**2))\n",
    "    return m\n",
    "\n",
    "def avgEdgeLength(cities, edges):\n",
    "    length = pathLength(cities, edges)\n",
    "    return length/len(cities)\n",
    "\n",
    "def optimizeStep(f, cities, edges, temperature, length):\n",
    "    newConns = list(edges)\n",
    "    for n in range(0,2):\n",
    "        i = random.randint(0,len(edges) - 1)\n",
    "        j = random.randint(0,len(edges) - 1)\n",
    "        newConns[i],newConns[j] = newConns[j],newConns[i]\n",
    "    newLength = f(cities, newConns)\n",
    "    delta = (newLength - length) / length\n",
    "    if delta > temperature:\n",
    "        return (length, edges)\n",
    "    else:\n",
    "        return (newLength, newConns)\n",
    "\n",
    "def optimize(f, cities, edges, steps=100000, decay=10000):\n",
    "    length = f(cities, edges)\n",
    "    for i in range(0, steps):\n",
    "        temperature = math.exp(-i/decay)\n",
    "        length, edges = optimizeStep(f, cities, edges, temperature, length)\n",
    "    return edges\n",
    "\n",
    "def targetValue(cities, edges):\n",
    "    return pathLength(cities, edges) + 2*maxEdgeLength(cities, edges)\n",
    "\n",
    "cities = makeCities()\n",
    "edges = connectCities(cities)\n",
    "\n",
    "edges = optimize(targetValue, cities, edges, steps=1000000)\n",
    "#edges = optimize(pathLength, cities, edges)\n",
    "#edges = optimize(maxEdgeLength, cities, edges)\n",
    "\n",
    "avgEdgeLength(cities, edges)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\n",
       "\n",
       "mpl.get_websocket_type = function() {\n",
       "    if (typeof(WebSocket) !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof(MozWebSocket) !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert('Your browser does not have WebSocket support.' +\n",
       "              'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
       "              'Firefox 4 and 5 are also supported but you ' +\n",
       "              'have to enable WebSockets in about:config.');\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n",
       "    this.id = figure_id;\n",
       "\n",
       "    this.ws = websocket;\n",
       "\n",
       "    this.supports_binary = (this.ws.binaryType != undefined);\n",
       "\n",
       "    if (!this.supports_binary) {\n",
       "        var warnings = document.getElementById(\"mpl-warnings\");\n",
       "        if (warnings) {\n",
       "            warnings.style.display = 'block';\n",
       "            warnings.textContent = (\n",
       "                \"This browser does not support binary websocket messages. \" +\n",
       "                    \"Performance may be slow.\");\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.imageObj = new Image();\n",
       "\n",
       "    this.context = undefined;\n",
       "    this.message = undefined;\n",
       "    this.canvas = undefined;\n",
       "    this.rubberband_canvas = undefined;\n",
       "    this.rubberband_context = undefined;\n",
       "    this.format_dropdown = undefined;\n",
       "\n",
       "    this.image_mode = 'full';\n",
       "\n",
       "    this.root = $('<div/>');\n",
       "    this._root_extra_style(this.root)\n",
       "    this.root.attr('style', 'display: inline-block');\n",
       "\n",
       "    $(parent_element).append(this.root);\n",
       "\n",
       "    this._init_header(this);\n",
       "    this._init_canvas(this);\n",
       "    this._init_toolbar(this);\n",
       "\n",
       "    var fig = this;\n",
       "\n",
       "    this.waiting = false;\n",
       "\n",
       "    this.ws.onopen =  function () {\n",
       "            fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n",
       "            fig.send_message(\"send_image_mode\", {});\n",
       "            fig.send_message(\"refresh\", {});\n",
       "        }\n",
       "\n",
       "    this.imageObj.onload = function() {\n",
       "            if (fig.image_mode == 'full') {\n",
       "                // Full images could contain transparency (where diff images\n",
       "                // almost always do), so we need to clear the canvas so that\n",
       "                // there is no ghosting.\n",
       "                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "            }\n",
       "            fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "        };\n",
       "\n",
       "    this.imageObj.onunload = function() {\n",
       "        this.ws.close();\n",
       "    }\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_header = function() {\n",
       "    var titlebar = $(\n",
       "        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n",
       "        'ui-helper-clearfix\"/>');\n",
       "    var titletext = $(\n",
       "        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n",
       "        'text-align: center; padding: 3px;\"/>');\n",
       "    titlebar.append(titletext)\n",
       "    this.root.append(titlebar);\n",
       "    this.header = titletext[0];\n",
       "}\n",
       "\n",
       "\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = $('<div/>');\n",
       "\n",
       "    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
       "\n",
       "    function canvas_keyboard_event(event) {\n",
       "        return fig.key_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    canvas_div.keydown('key_press', canvas_keyboard_event);\n",
       "    canvas_div.keyup('key_release', canvas_keyboard_event);\n",
       "    this.canvas_div = canvas_div\n",
       "    this._canvas_extra_style(canvas_div)\n",
       "    this.root.append(canvas_div);\n",
       "\n",
       "    var canvas = $('<canvas/>');\n",
       "    canvas.addClass('mpl-canvas');\n",
       "    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
       "\n",
       "    this.canvas = canvas[0];\n",
       "    this.context = canvas[0].getContext(\"2d\");\n",
       "\n",
       "    var rubberband = $('<canvas/>');\n",
       "    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
       "\n",
       "    var pass_mouse_events = true;\n",
       "\n",
       "    canvas_div.resizable({\n",
       "        start: function(event, ui) {\n",
       "            pass_mouse_events = false;\n",
       "        },\n",
       "        resize: function(event, ui) {\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "        stop: function(event, ui) {\n",
       "            pass_mouse_events = true;\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "    });\n",
       "\n",
       "    function mouse_event_fn(event) {\n",
       "        if (pass_mouse_events)\n",
       "            return fig.mouse_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    rubberband.mousedown('button_press', mouse_event_fn);\n",
       "    rubberband.mouseup('button_release', mouse_event_fn);\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband.mousemove('motion_notify', mouse_event_fn);\n",
       "\n",
       "    rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
       "    rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
       "\n",
       "    canvas_div.on(\"wheel\", function (event) {\n",
       "        event = event.originalEvent;\n",
       "        event['data'] = 'scroll'\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        mouse_event_fn(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.append(canvas);\n",
       "    canvas_div.append(rubberband);\n",
       "\n",
       "    this.rubberband = rubberband;\n",
       "    this.rubberband_canvas = rubberband[0];\n",
       "    this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
       "    this.rubberband_context.strokeStyle = \"#000000\";\n",
       "\n",
       "    this._resize_canvas = function(width, height) {\n",
       "        // Keep the size of the canvas, canvas container, and rubber band\n",
       "        // canvas in synch.\n",
       "        canvas_div.css('width', width)\n",
       "        canvas_div.css('height', height)\n",
       "\n",
       "        canvas.attr('width', width);\n",
       "        canvas.attr('height', height);\n",
       "\n",
       "        rubberband.attr('width', width);\n",
       "        rubberband.attr('height', height);\n",
       "    }\n",
       "\n",
       "    // Set the figure to an initial 600x600px, this will subsequently be updated\n",
       "    // upon first draw.\n",
       "    this._resize_canvas(600, 600);\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n",
       "        return false;\n",
       "    });\n",
       "\n",
       "    function set_focus () {\n",
       "        canvas.focus();\n",
       "        canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    window.setTimeout(set_focus, 100);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            // put a spacer in here.\n",
       "            continue;\n",
       "        }\n",
       "        var button = $('<button/>');\n",
       "        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
       "                        'ui-button-icon-only');\n",
       "        button.attr('role', 'button');\n",
       "        button.attr('aria-disabled', 'false');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "\n",
       "        var icon_img = $('<span/>');\n",
       "        icon_img.addClass('ui-button-icon-primary ui-icon');\n",
       "        icon_img.addClass(image);\n",
       "        icon_img.addClass('ui-corner-all');\n",
       "\n",
       "        var tooltip_span = $('<span/>');\n",
       "        tooltip_span.addClass('ui-button-text');\n",
       "        tooltip_span.html(tooltip);\n",
       "\n",
       "        button.append(icon_img);\n",
       "        button.append(tooltip_span);\n",
       "\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    var fmt_picker_span = $('<span/>');\n",
       "\n",
       "    var fmt_picker = $('<select/>');\n",
       "    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
       "    fmt_picker_span.append(fmt_picker);\n",
       "    nav_element.append(fmt_picker_span);\n",
       "    this.format_dropdown = fmt_picker[0];\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = $(\n",
       "            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n",
       "        fmt_picker.append(option)\n",
       "    }\n",
       "\n",
       "    // Add hover states to the ui-buttons\n",
       "    $( \".ui-button\" ).hover(\n",
       "        function() { $(this).addClass(\"ui-state-hover\");},\n",
       "        function() { $(this).removeClass(\"ui-state-hover\");}\n",
       "    );\n",
       "\n",
       "    var status_bar = $('<span class=\"mpl-message\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n",
       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
       "    // which will in turn request a refresh of the image.\n",
       "    this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_message = function(type, properties) {\n",
       "    properties['type'] = type;\n",
       "    properties['figure_id'] = this.id;\n",
       "    this.ws.send(JSON.stringify(properties));\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_draw_message = function() {\n",
       "    if (!this.waiting) {\n",
       "        this.waiting = true;\n",
       "        this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n",
       "    }\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    var format_dropdown = fig.format_dropdown;\n",
       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
       "    fig.ondownload(fig, format);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_resize = function(fig, msg) {\n",
       "    var size = msg['size'];\n",
       "    if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n",
       "        fig._resize_canvas(size[0], size[1]);\n",
       "        fig.send_message(\"refresh\", {});\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
       "    var x0 = msg['x0'];\n",
       "    var y0 = fig.canvas.height - msg['y0'];\n",
       "    var x1 = msg['x1'];\n",
       "    var y1 = fig.canvas.height - msg['y1'];\n",
       "    x0 = Math.floor(x0) + 0.5;\n",
       "    y0 = Math.floor(y0) + 0.5;\n",
       "    x1 = Math.floor(x1) + 0.5;\n",
       "    y1 = Math.floor(y1) + 0.5;\n",
       "    var min_x = Math.min(x0, x1);\n",
       "    var min_y = Math.min(y0, y1);\n",
       "    var width = Math.abs(x1 - x0);\n",
       "    var height = Math.abs(y1 - y0);\n",
       "\n",
       "    fig.rubberband_context.clearRect(\n",
       "        0, 0, fig.canvas.width, fig.canvas.height);\n",
       "\n",
       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_figure_label = function(fig, msg) {\n",
       "    // Updates the figure title.\n",
       "    fig.header.textContent = msg['label'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_cursor = function(fig, msg) {\n",
       "    var cursor = msg['cursor'];\n",
       "    switch(cursor)\n",
       "    {\n",
       "    case 0:\n",
       "        cursor = 'pointer';\n",
       "        break;\n",
       "    case 1:\n",
       "        cursor = 'default';\n",
       "        break;\n",
       "    case 2:\n",
       "        cursor = 'crosshair';\n",
       "        break;\n",
       "    case 3:\n",
       "        cursor = 'move';\n",
       "        break;\n",
       "    }\n",
       "    fig.rubberband_canvas.style.cursor = cursor;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_message = function(fig, msg) {\n",
       "    fig.message.textContent = msg['message'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_draw = function(fig, msg) {\n",
       "    // Request the server to send over a new figure.\n",
       "    fig.send_draw_message();\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_image_mode = function(fig, msg) {\n",
       "    fig.image_mode = msg['mode'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message(\"ack\", {});\n",
       "}\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function(fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            /* FIXME: We get \"Resource interpreted as Image but\n",
       "             * transferred with MIME type text/plain:\" errors on\n",
       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "             * to be part of the websocket stream */\n",
       "            evt.data.type = \"image/png\";\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src);\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                evt.data);\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n",
       "            fig.imageObj.src = evt.data;\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        var msg = JSON.parse(evt.data);\n",
       "        var msg_type = msg['type'];\n",
       "\n",
       "        // Call the  \"handle_{type}\" callback, which takes\n",
       "        // the figure and JSON message as its only arguments.\n",
       "        try {\n",
       "            var callback = fig[\"handle_\" + msg_type];\n",
       "        } catch (e) {\n",
       "            console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        if (callback) {\n",
       "            try {\n",
       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
       "                callback(fig, msg);\n",
       "            } catch (e) {\n",
       "                console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "}\n",
       "\n",
       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
       "mpl.findpos = function(e) {\n",
       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
       "    var targ;\n",
       "    if (!e)\n",
       "        e = window.event;\n",
       "    if (e.target)\n",
       "        targ = e.target;\n",
       "    else if (e.srcElement)\n",
       "        targ = e.srcElement;\n",
       "    if (targ.nodeType == 3) // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "\n",
       "    // jQuery normalizes the pageX and pageY\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    // offset() returns the position of the element relative to the document\n",
       "    var x = e.pageX - $(targ).offset().left;\n",
       "    var y = e.pageY - $(targ).offset().top;\n",
       "\n",
       "    return {\"x\": x, \"y\": y};\n",
       "};\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * http://stackoverflow.com/a/24161582/3208463\n",
       " */\n",
       "function simpleKeys (original) {\n",
       "  return Object.keys(original).reduce(function (obj, key) {\n",
       "    if (typeof original[key] !== 'object')\n",
       "        obj[key] = original[key]\n",
       "    return obj;\n",
       "  }, {});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.mouse_event = function(event, name) {\n",
       "    var canvas_pos = mpl.findpos(event)\n",
       "\n",
       "    if (name === 'button_press')\n",
       "    {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x;\n",
       "    var y = canvas_pos.y;\n",
       "\n",
       "    this.send_message(name, {x: x, y: y, button: event.button,\n",
       "                             step: event.step,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "\n",
       "    /* This prevents the web browser from automatically changing to\n",
       "     * the text insertion cursor when the button is pressed.  We want\n",
       "     * to control all of the cursor setting manually through the\n",
       "     * 'cursor' event from matplotlib */\n",
       "    event.preventDefault();\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    // Handle any extra behaviour associated with a key event\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.key_event = function(event, name) {\n",
       "\n",
       "    // Prevent repeat events\n",
       "    if (name == 'key_press')\n",
       "    {\n",
       "        if (event.which === this._key)\n",
       "            return;\n",
       "        else\n",
       "            this._key = event.which;\n",
       "    }\n",
       "    if (name == 'key_release')\n",
       "        this._key = null;\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.which != 17)\n",
       "        value += \"ctrl+\";\n",
       "    if (event.altKey && event.which != 18)\n",
       "        value += \"alt+\";\n",
       "    if (event.shiftKey && event.which != 16)\n",
       "        value += \"shift+\";\n",
       "\n",
       "    value += 'k';\n",
       "    value += event.which.toString();\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, {key: value,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onclick = function(name) {\n",
       "    if (name == 'download') {\n",
       "        this.handle_save(this, null);\n",
       "    } else {\n",
       "        this.send_message(\"toolbar_button\", {name: name});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n",
       "    this.message.textContent = tooltip;\n",
       "};\n",
       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to  previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
       "\n",
       "mpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n",
       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
       "    // object with the appropriate methods. Currently this is a non binary\n",
       "    // socket, so there is still some room for performance tuning.\n",
       "    var ws = {};\n",
       "\n",
       "    ws.close = function() {\n",
       "        comm.close()\n",
       "    };\n",
       "    ws.send = function(m) {\n",
       "        //console.log('sending', m);\n",
       "        comm.send(m);\n",
       "    };\n",
       "    // Register the callback with on_msg.\n",
       "    comm.on_msg(function(msg) {\n",
       "        //console.log('receiving', msg['content']['data'], msg);\n",
       "        // Pass the mpl event to the overriden (by mpl) onmessage function.\n",
       "        ws.onmessage(msg['content']['data'])\n",
       "    });\n",
       "    return ws;\n",
       "}\n",
       "\n",
       "mpl.mpl_figure_comm = function(comm, msg) {\n",
       "    // This is the function which gets called when the mpl process\n",
       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
       "\n",
       "    var id = msg.content.data.id;\n",
       "    // Get hold of the div created by the display call when the Comm\n",
       "    // socket was opened in Python.\n",
       "    var element = $(\"#\" + id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm)\n",
       "\n",
       "    function ondownload(figure, format) {\n",
       "        window.open(figure.imageObj.src);\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy,\n",
       "                           ondownload,\n",
       "                           element.get(0));\n",
       "\n",
       "    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
       "    // web socket which is closed, not our websocket->open comm proxy.\n",
       "    ws_proxy.onopen();\n",
       "\n",
       "    fig.parent_element = element.get(0);\n",
       "    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n",
       "    if (!fig.cell_info) {\n",
       "        console.error(\"Failed to find cell for figure\", id, fig);\n",
       "        return;\n",
       "    }\n",
       "\n",
       "    var output_index = fig.cell_info[2]\n",
       "    var cell = fig.cell_info[0];\n",
       "\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function(fig, msg) {\n",
       "    fig.root.unbind('remove')\n",
       "\n",
       "    // Update the output cell to use the data from the current canvas.\n",
       "    fig.push_to_output();\n",
       "    var dataURL = fig.canvas.toDataURL();\n",
       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
       "    // the notebook keyboard shortcuts fail.\n",
       "    IPython.keyboard_manager.enable()\n",
       "    $(fig.parent_element).html('<img src=\"' + dataURL + '\">');\n",
       "    fig.close_ws(fig, msg);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.close_ws = function(fig, msg){\n",
       "    fig.send_message('closing', msg);\n",
       "    // fig.ws.close()\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.push_to_output = function(remove_interactive) {\n",
       "    // Turn the data on the canvas into data in the output cell.\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\">';\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Tell IPython that the notebook contents must change.\n",
       "    IPython.notebook.set_dirty(true);\n",
       "    this.send_message(\"ack\", {});\n",
       "    var fig = this;\n",
       "    // Wait a second, then push the new image to the DOM so\n",
       "    // that it is saved nicely (might be nice to debounce this).\n",
       "    setTimeout(function () { fig.push_to_output() }, 1000);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items){\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) { continue; };\n",
       "\n",
       "        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n",
       "    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n",
       "    button.click(function (evt) { fig.handle_close(fig, {}); } );\n",
       "    button.mouseover('Stop Interaction', toolbar_mouse_event);\n",
       "    buttongrp.append(button);\n",
       "    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n",
       "    titlebar.prepend(buttongrp);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(el){\n",
       "    var fig = this\n",
       "    el.on(\"remove\", function(){\n",
       "\tfig.close_ws(fig, {});\n",
       "    });\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(el){\n",
       "    // this is important to make the div 'focusable\n",
       "    el.attr('tabindex', 0)\n",
       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
       "    // off when our div gets focus\n",
       "\n",
       "    // location in version 3\n",
       "    if (IPython.notebook.keyboard_manager) {\n",
       "        IPython.notebook.keyboard_manager.register_events(el);\n",
       "    }\n",
       "    else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    var manager = IPython.notebook.keyboard_manager;\n",
       "    if (!manager)\n",
       "        manager = IPython.keyboard_manager;\n",
       "\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which == 13) {\n",
       "        this.canvas_div.blur();\n",
       "        event.shiftKey = false;\n",
       "        // Send a \"J\" for go to next cell\n",
       "        event.which = 74;\n",
       "        event.keyCode = 74;\n",
       "        manager.command_mode();\n",
       "        manager.handle_keydown(event);\n",
       "    }\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.find_output_cell = function(html_output) {\n",
       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
       "    // IPython event is triggered only after the cells have been serialised, which for\n",
       "    // our purposes (turning an active figure into a static one), is too late.\n",
       "    var cells = IPython.notebook.get_cells();\n",
       "    var ncells = cells.length;\n",
       "    for (var i=0; i<ncells; i++) {\n",
       "        var cell = cells[i];\n",
       "        if (cell.cell_type === 'code'){\n",
       "            for (var j=0; j<cell.output_area.outputs.length; j++) {\n",
       "                var data = cell.output_area.outputs[j];\n",
       "                if (data.data) {\n",
       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
       "                    data = data.data;\n",
       "                }\n",
       "                if (data['text/html'] == html_output) {\n",
       "                    return [cell, data, j];\n",
       "                }\n",
       "            }\n",
       "        }\n",
       "    }\n",
       "}\n",
       "\n",
       "// Register the function which deals with the matplotlib target/channel.\n",
       "// The kernel may be null if the page has been refreshed.\n",
       "if (IPython.notebook.kernel != null) {\n",
       "    IPython.notebook.kernel.comm_manager.register_target('matplotlib', mpl.mpl_figure_comm);\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4nOzdeZzN1f8H8DMzmGtfE6HMhCIpVBpJ0aZEi1Iq0UKLqFFyU9JiaVGkTbZKyq2QFIa+tmixZAkpso3sjG3WO3Pv6/fHe+ZnMMy9M/d+zud+zuv5eNzH9xv3fu7bZe593bO9lSIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiilT3KKV+VkodUUr5lFLRJ/1+E6XUIqVUqlLqP6XUIEurIyIiIqKQu0FJCHxInRoAyymldimlBiulSimlGiuldiilnra4RiIiIiIKg2vUqQGwm1Jqz0m/1kcptcnCuoiIiIgoTAoKgO8qpWafdL+E3PuVs6guIiIiIgqTggLgOKXU5JPud2Hu/c6xqC4iIiIiCpNQjABGKaVqKaUq8MYbb7zxZtytlpLPASKKIAUFwAfVqWsAn1anXwNYSykF3njjjTfejL3VUkQUEaKVUrFKqRuVBMAyuf8dpWSUb6dS6nWllEspdbFSKlmdfhdwBaUUduzYgSNHjvAW4K1Xr17aa4i0G18zvl5Oes0SExOh1HtQ6ki+20okJEzDggVHcOiQ/tejsNuOHTvyAmCFsH5iEVHIdFNK+ZWEP1++/9869/cbKzknME3JkTADz3CtCkopHDlyBBS4xMRE3SVEHL5mweHrFTwrX7OHH14BpY5BqQwoBSiVhhIl+qNFi2RUrAhUqwbcfz8wcSKwZ8+Jj/V6vfB4PHC73fB4PPB6vZbVnd+RI0cYAIkMxgBYBPxwDh5fs+Dw9QqeFa+Z3w+89hpQubIfjRv3gssVh5iYHnC54pCQ0AZerxfZ2cCSJcBLLwGXXQZERwPNmgEDBgDz52ejRYvr4XLF5z4u/v8fZzUGQCKzMQAWQVJSku4SIg5fs+Dw9QpeuF8znw94+mmgZk3gzz8DH8nbtw/48kuga1egYsUMKHUUSmXnjhymwuWKg8fjCWvtBWEAJDIbAyARUSG8Xglw8fHA5s1Fv07//i8gOvr13PAnt5iYHnC73aErNkAMgERmYwAkIjqD9HSgQwfg4ouBXbuKdy2PxwOXKx5KpeYGwEyOABKRFgyARESncfgwcM01QMuWQEpK8a/n9XqRkNAGLlccoqImIirqD64BJCItGACJiAqwdy/QtCnQrh2Qmhq66+atHezT5zXExPiwaRN3AROR9YwNgHY5ioGI7Gf7dqBBA+Dee4GsrPA9z003AW+9Fb7rnwkDIJHZjAyAx6dh9B/FQET28tdfQO3awBNPADk54X2uTz+VI2J0YAAkMpuRAfDUhdhpiI1tqGUhNhHZx/LlQNWqcoaf3x/+5zt0CChVCti0KfzPdTIGQCKzGRkA3W43YmJ6nHAUg1JA6dLH0Lw5cMcdwDPPACNGANOmAStWAPv3W/OBQER6zJsHlC8PvPuutc/boQMweLC1zwkwABKZzsgAeHwEMK+NUypiYy/DkCFzMGUK8M47QJ8+wG23ySLwKlUkIJYpAzRsKIvCe/YEhgwBJk0CFi+WNUPZ2br/ZERUFN99Jz/fn31m/XN/+SXQuLH1z8sASGQ2IwNg3hrAUqWugVLZiI2tV+gawGPHgHXrgFmzgI8/BtxuoEsX4KqrZL1QdDQQEwOcdx7QujXwwAPAiy8CY8YAc+YAGzYAaWnW/RmJKDATJgBlywLTp+t5/qNHAZcLWL/e2udlACQym5EBEJAQ+OWXXyM6OhsjR/5Q7A0gXi+wZQuwcCHw+efSL/SRR4AbbpDdhC6XjCKedZb0B73zTiAxERg5UkYf/vgDOHCA08xEVnrnHaBCBWDBAr11dOoEvPyytc/JAEhkNmMDYJ4LL5RRvXDz+4E9e4Bly4BvvwWGDwd695Zp5ksvBSpXloBYtizQqJFMMz/2GDB0qEwRLV4MJCeHf1cikQn8fmDAAKBaNVnjq9u338oXRSu/ADIAEpnN+ADYsaOMwtnB0aPA2rXAzJnARx8B/fvLOWQtWwK1agFRUUCJEkDdutKdoGtX2a04diwwdy7w99/StoqITi8nR75c1akjSzPsIC1NvvytWmXdczIAEpnN+AD43HPAk0/qriIwWVnSiH7BAlms/uqrwMMPA9dfD9SvD8TGyihi9erA5ZfLtFLfvsB778n6ppUrgYMHOc1M5srKAjp3Bi64QDZu2UmXLrK22CoMgERmMz4AfvKJrNNzAp8P2L0bWLoU+OYb4O23gaeeklHOSy4BKlWSgFiuHHDRRcDNNwOPPw4MGwZ89RWwZAmwYwenmcmZUlOl80bz5sC+fbqrOdX338vovlVf0BgAicxmfABcsEB27priyBGZZv7xR+DDD4HnnwfuuQdISADOOef4NHNcnEwzP/ggMHAgMG6cTDP/8w+QkaH7T0EUnIMH5d/4tdfKz4AdZWYCFSvKFzgrMAASmc34ALhzp4Qerp0TWVnAv/8C8+dLm6pXXgEeegi47jqgXj3pWqAUcPbZwBVXAHfdBTz7LDBqlIxgrFoFpKRwmpnsY9cuOWevY0f7f3np1k1OB7ACAyCR2YwPgH6/LL5eu1Z3JZHB55MP1N9/B77+WhrZ9+ol3QyaNJERDKWko0LjxsAtt0hP1TfeACZPBn75BfjvP7nOmXi9Xng8Hrjdbng8HvZppiL5918Zze7WLTIOap89WzZ8FfbzEQoMgERmMz4AAtLtY+pU3VU4x+HDwJ9/Aj/8AHzwAdCvn0wzX3klULOmBMSSJYH4eJmS69ZNzkAbPx746Sdg/XovWrS4ES5XPGJiesDlii/0oG6ik61ZA9SoISNqVgSqUPB6pRfx4sXhfy4GQCKzMQBCwsmwYbqrMEdmpozMzJsnXRgGDQK6dwfatgXOPx+IicnJbdHn+/9WfS5XHDwej+7SKUL88otseho8OPKWI/TsKaPq4cYASGQ2BkDIWXoPPaS7CsrTv/8LiI4elhv+5BYT0wNuK8/IoIg1e7bsdP/oI92VFM28ebLGNty78RkAiczGAAhg4kTp6Uv24PF4ULLkfRwBpKBNngyUKSP/G6lyciQAzpsX3udhACQyGwMgZEND9eq6q6A8Xq8X9eoNRVTU77lrAOO4BpAK9fHHMvI3e7buSorvqaeAHj3C+xwMgERmYwCEnBGmFHDokO5KKM+rr+agVaut3AVMhfL7Za1fpUpymLkTLF4MVKkim0LChQGQyGwMgLmqVQOWLdNdBeV57DFr22JRZPL5ZJfv2WfLrl+n8PnkOJhwjmYyABKZjQEwV8uWwKRJuqugPLfcIp1KiE4nO1uOEIqLk13lTtO3r/z5woUBkMhsDIC5uneXs+jIHpo0AWbM0F0F2VVGBnDbbXLY+M6duqsJj6VL5WD1zMzwXJ8BkMhsDIC5hgwB7r1XdxWUp3JlaStHdLIjR4A2baS378GDuqsJH79fRje//z4812cAJDIbA2Cub78FmjXTXQUBQGqqbMrZv193JWQ3+/YBzZsDN94o/06czu0GunQJz7UZAInMxgCYa80aOUIi0roGONGGDYDLxb8LOlFyMnDBBUDnzkBWlu5qrLFqlfQqT0sL/bUZAInMxgCYKy1NRp1279ZdCc2dC9Svr7sKspMNG4A6dWR3eLg7ZNiJ3w80aCAzFKHGAEhkNgbAfOrUARYt0l0FjR8vfYGJAGDFCjmmacAAM0eFX34Z6NQp9NdlACQyGwNgPtddB4wdq7sKeuWV8B5/QZFjwQKgQgVg+HDdleizfr0siTh6NLTXZQAkMhsDYD5PPAH066e7CnrkEeCll3RXQbpNny7r3yZM0F2Jfo0bA19+GdprMgASmY0BMJ9335WzxUivG28EPvlEdxWk02efAWXKANOm6a7EHgYPBjp0CO01GQCJzMYAmM+PPwING+qugho1AmbN0l0F6TJiBFC+PDBvnu5K7GPTJqBUKSAlJXTXZAAkMhsDYD4bN8qbrEm7DO2ofHngzz91V0FW8/tl6r9qVfblLkizZqGdDmcAJDIbA2A+aWlexMT48Pjjb8Lj8cDr9eouyTiHD8txPIcO6a6ErOTzyRrc2rWBv/7SXY09vfUWcNNNobseAyCR2RgAc3m9XiQktEFU1GZER4+AyxWPhIQ2DIEWW7uWB3KbJitL2jA2aABs3667Gvvatg0oUSJ0HXIYAInMxgCYy+PxwOWKh1LZUApQKhUuVxw8Ho/u0owyaxbXYZokLQ1o1w5o2hTYu1d3NfZ35ZXA6NGhuRYDIJHZGABzud1uxMT0yA1/couJ6QG32627NKN88onsAibnS0kBWrYEWreWqX8q3MiRQJs2obkWAyCR2RgAcx0fAUzNDYA+xMZexhFAi730kpwDSM7j9Xrh8Xjgdrvx8cfTcfHFfnToAKSn664scuzcKdPAu3YV/1oMgERmYwDMlbcG0OWKQ3R0D0RH/4iKFZchM5NrAK3UrZt0AiFnOf7zFY/o6BcQFbUN1arNQVoaf76C1bo1MGpU8a/DAEhkNgbAfPKPUEyYMAX16vnx+uu6qzJL27bSC5ic5fgIe3ruCLsXsbHxHGEvgo8+Aq66qvjXYQAkMhsD4BmsWiXdCBYs0F2JOerXB+bM0V0Fhdqjj74LpVZzjW0I7N0r08DJycW7DgMgkdkYAAsxejRQowawZ4/uSpzP7wdKl+Y5cE6yaxfQowdQsmQ2YmI+z7fGlrvsi+OGG4Dhw4t3DQZAIrMxABbC75czyq67jh1Cwu3AARkZOnZMdyVUXEePAi+/DJQtC9x1F7B+/fE1tjExPeByxfGczWIYNw64/PLiXYMBkMhsDIABOHpUDqnl5oTwWrUKqFRJdxVUHF6vrFGrXl3Wqf32W/7fO77Glp12iufgQaBkSWDz5qJfgwGQyGwMgAFas0ZGM9igPnxmzAAuvlh3FVQUfj/w3XfABRfIbfp0dnMJt/btgaFDi/54BkAiszEABmHsWODss4Hdu3VX4kwffgjccovuKihYv/4qo33VqwMffwxkZ+uuyAwTJwKXXFL0xzMAEpmNATAIfj/wwANyEj/XA4ae2w089pjuKihQGzcCnTrJyPigQbJUgqxz5AgQGwts2FDUxzMAEpmMATBIx44BF1zgx513ruNaphDyer1o1WorWrdO4mtqc3v3Ak89BbhcQM+eoelKQUVzxx1FX5vMAEhkNgbAIHm9XlxyyUNQKhXR0SPgcsVzN2Mx5XWJiIr6HVFRE/ia2lRaGjB4MFC+PNCxI4/rsQOPB2jYsGjrLRkAiczGABik4x0NsnieWYh4PB7Exl4Mpfx8TW0oJ0e6s5xzDnDFFcCiRborojypqUCZMn689dbsoGckGACJzMYAGCS3243o6L7saBBCbrcbUVHz+ZrajN8PzJwJXHQREB8PfP01d/bajdfrRdWq8xAT82Hu+YqBj54zABKZjQEwSF9++TWio3+DUtkcrQqRV1/9H5RKy73xNbWDFStks1PVqsB77wFZWborMlNhZyd6PB6ULPkklPIF/bPDAEhkNgbAIPXtm4PSpbcgNvYidjQIgYwM2VRz7rmj2SXCBrZsAe67T1ryud3AoUO6KzJX3tpYlyv+hNG95GQvpk4Fnn0WqFVrG5TyFmn0nAGQyGwMgEH45hugYkVpa8WOBqHx4otA8+ZAejpfU50OHgT69pWdvd27A8nJuiui4+uN80bGsxEVtQ1RUX5cfLEcmfTkk7+jVKlri9RjmQGQyGwMgAFatw4oVw74/nvdlTjHmjUy0rR6te5KzJWRAbz9trTgu+km/l3Yiaw3fu+E0b3o6BF45pnj574cHyUMfvScAZDIbAyAATh8WHoBv/ii7kqcIzsbuOwyYMAA3ZWYyecDvvgCOPdc4NJLgblzdVdEJxswYEHuyF7GGUf3itpjmQGQyGwMgIXw+YDbb5fREXb/CJ3hw6VnbEaG7krM89NPQNOmEv6++EL+jZO9LFgAlC3rR716r4dtbSwDIJHZGAALMXQoULcucOCA7kqc499/gTJlgJ9/1l2JWdasAdq1k+net95i+LarJUvksO0JE4o+uhcIBkAiszEAniT/G+4LLyxA2bJ+rFypuyrn8PuBtm2BJ5/UXYk5duyQjR0ul2z04JcZ+/r9d6BCBWD06PA/FwMgkdkYAPPJf+xCdLQbSqXg/POHcUdqCI0bB9SuLY3sKbwOHwZeeEFGW7t0kSNeyL7++ENGZ0eNsub5GACJzMYAmM/xYxfyjlTw8kDiENq5Uz7gfvxRdyXOlpUlhzdXrQpcey2wfLnuiqgwa9YAVarIjmyrMAASmY0BMB85dqEHW5KFyR13yEgUhYffL+3a4uOlfdvMmWzdFgnWrwfOOgsYPNja52UAJHKm6kqpr5RSe5RSKUqpX5RSrQu4HwNgPh6PByVKvFmktkp0ZlOmyIjUvn26K3GmRYuAK64AzjlHptm5Yz0y/PMPUKMG8PLL1j83AyCRM01VSi1SSlVRSkUppfoqpY4qpSqddD8GwHymTctGdHQ6SpVqz5ZkIZSSIh9ykybprsR5/voL6NhRdo0OHgykpuquiAL1779ArVpA//56RmoZAImcabVSqne+/y6rlPIrpS4/6X4MgLnWrZPdd199lc2WZCH28MPAzTdzOjKUdu8GevaUnb29egF79+quiIKxbZucw/jMM/p+LhgAiZypi1JqvlKqhlKqpFKqv1Jqo1Iq9qT7MQBCjsWIj9czDeN0P/0ko1Pbt+uuxBmOHQMGDQLKlgU6dZIpRIosO3bI+82TT+r9UsQASORM5yqlZikZ9fMqWQuYUMD9jA+AXi/Qpo1sUGBHhNBKTQXi4oAPPtBdSeTLzgY+/hg4+2ygZUvg1191V0RFsWuXtJV89FH97zcMgETOE6WU2qyUGq+UqqiUilZKdVRKHVZKNTnpvsYHwCefBJo0kZEVCq2+fYGrrtL/QRfJ/H5g+nRpm9egAfDdd5xKj1T79gGNGgFdu9rjZ4IBkMh5qigZ+bvkpF//QynV76Rfq6CUQq9evZCYmIjExEQkJSXpfl+yzMcfA9WqAVu36q7EeZYuBUqXBjZs0F1J5PrtN6BVK6B6deCjj2S0miLTgQPyRfOee2Q0V5ekpKT/f6/v1asXAyCRA61TSo1RSpVXMiJ4q1IqQynV5qT7GTsCuGAB+9GGS1YW0Lgx8PrruiuJTJs2AXfdJev8Bg4Ejh7VXREVx6FDQLNmsszETiGeI4BEznS+Uuo7pdReJVO/a5VSjxRwPyMD4ObNcibd2LG6K3Gm116T0Y6sLN2VRJZ9+4DevWVnb48e0jmFItuRI3I+46232u/ngQGQyGzGBcCjR2V0qndv3ZU40/r1MvXL9mOBS0sDhgyR3dIdOshrSJHv2DFZA3vTTUBmpu5qTsUASGQ2owKgzwfcdhtw/fV61+E4VU4OkJAAPPus7koiQ04OMGGCHAZ82WWyLIGcIS1N+jC3bQukp+uupmAMgERmMyoADhgA1KsHHDyouxJnGjVKzjdLS9Ndib35/cCsWTISHRcHeDz22BVKoZGRAdxwA3D11fbuzMIASGQ2YwLgV18BFStK6ywKvW3bgHLlgHnzdFdib3/8IaNCVaoAI0fac2qQii4rC2jfHrjySvtv3mEAJDKbEQFw+XLZUTlzpu5KnMnvB9q1Ax55RHcl9rV1K3D//bI+sn9/2RlKzuL1ArffDjRvHhl/vwyARGZzfADctUvWWL31lu5KnGviRKBGDSAlRXcl9uD1ev+/n/S4cdOQmJgDlwvo1o0t8ZwqOxvo3Bm45JLIWWLCAEhkNkcHwIwMoEULOXmf3RPCY+9emc6cOlV3Jfbg9XqRkNAGsbEXIirqGyh1CBUrLsOyZTY6AI5CKidHRncvukiO8okUDIBEZnNsAPT7Jfi1aCFBkMLj3nuBTp10V2EfHo8HLlc8lMqGUoBSGXC54uDxeHSXRmHg8wEPPSSt+vbs0V1NcBgAiczm2AD49tsy9btrl+5KnGvGDKBSJWD3bt2V2Ifb7UZ0dGJu+JNbTEwPuN1u3aVRiPn9wOOPA+efD/z3n+5qgscASGQ2RwbAmTNl0wcPIw6fw4clYI8fr7sSe/F4PChR4nUolZMbAFM5AuhAfj/Qpw9Qt27krutkACQym+MC4F9/ARUqyLEvFD6PPw5cdx3XVp4sK8uL0qW3oGTJRMTE9IDLFYeEhDbw2qkJLBWL3w/06wfUrg1s2aK7mqJjACQym6MC4MGDctDzgAG6K3G2RYtkhHXzZt2V2M+yZUD58n58/vm3cLvd8Hg8DH8O89JLQM2awMaNuispHgZAIrM5JgBmZ0uLt9tuY1eFcEpPB+rXB955R3cl9vT440CPHrqroHB57TWgenVnHCjPAEhkNscEwD59pLWW3U/fj3RuN3D55XL0BZ0oLU2WH/z2m+5KKBzefBOoWhX480/dlYQGAyCR2RwRAMeOlTfmSF6PEwlWrpROFmvW6K7EniZNAho25LpIJxoxAqhcWX4GnIIBkMhsER8Af/4ZKFMGWLBAdyXOlp0NNG0KDByouxL7attWjh8iZ/noI+kjvmyZ7kpCiwGQyGwRHQC3bQPOOgv4+GPdlTjfG2/I6FZmpu5K7GnrVqBkycg7DJjObNw4oHx54JdfdFcSegyARGaL2AB47BjQpAnw5JO6K3G+jRtllNWJH4KhMmiQbEAi55g4EShXTna9OxEDIJHZIjIA+nzAHXcAbdoAPGEjvHw+4JprgKee0l2Jffl8wLnnAt9/r7sSChWPR446+ukn3ZWEDwMgkdkiMgC+/DIQHw8cOKC7Euf75BMJN9xdfXo//QScfTa/jDjFtGky4j1rlu5KwosBkMhsERcAv/lG1uSsW6e7Euf77z851mT2bN2V2FuXLtIZgiLfDz9I+DNhNJcBkMhsERUAV66UaZkZM3RX4nx+P9CxI/DAA7orsbeUFCA21hkHA5suKUnC37ff6q7EGgyARGaLmAC4Zw9Qpw4wbJjuSszw9ddAtWrA/v26K7G3Dz8ErrxSdxVUXPPmyZfLL7/UXYl1GACJzBYRATAzE2jZErjvPh6ya4UDB6Td1eTJuiuxv+bNgTFjdFdBwfJ6vfB4PHC73Rg06H8oV86Pzz7TXZW1GACJzGb7AOj3Aw89JO3H0tN1V2OGBx8Ebr2VYbswq1dLZxQb//g4Rv7A5vF44C3Gjhuv14uEhDZwueIRHT0USh1FfPzwYl0zEjEAEpnN9gHw3XeBmjVlQwKFX1KSbPzYsUN3Jfb39NMSlim88ge2mJgecLnikZDQ5v8Dm88HZGTITvWDB2W5yI4dwObNwN9/A2vXyvrhpUuBJUuAgQPno2TJB6FUBpQClMqEyxUHj8ej+U9qLQZAIrPZOgDmLcpeulR3JWY4dgw47zxpfUVnlpUl/acXLtRdifN5PB64XPFQKjs3sPmh1DGULJmD6Gjk/tqpt5gYwOWSLzRVqwI1asg64kqVDkCpXSfdtwfcbrfuP6qlGACJzGbbAPj330ClSsCkSborMcfTTwNXXy0jKnRmU6YA55/PaXIruN1uxMT0OCGwRUe/gp49h2PrVpkd2LtXdmQfOybh/Ez/ho8HytTc66VyBJCIjGPLAHjoENCgAdC/v+5KzPHrrzLa+s8/uiuJDLfcAgwerLuK8ArlurviODWw+RAd/TPGjZtWpOsdn1KOy51SjjthStkUDIBEZrNdAMzOBm66CWjfHsjJ0V2NGTIzgUaNgKFDdVcSGf77DyhZ0tnrJAtbd6enFglssbEXo3LlJYiP92Pt2qJf0w7hVicGQCKz2S4AJiYCDRtyZ6WVBg0CLr2UrcwCNXQo0K6d7irC6/ioW5otpklPDmyZmV4MGiRdgaZO1VJSxGMAJDKbrQLghAlA5crApk26KzHH2rVylMkff+iuJDL4/UC9etKS0MkKWndnx40S06ZJCBw4kGtXg8UASGQ22wTAX36RNWjz5umuxPnyRlOef/4F1Kt3AM8+y7n2QP38M1ClikybO9nxEcAcW4wAnsm6dbIhp0MHzhwEgwGQyGy2CIDbt0vniQ8+0FqGEfKv7YqK8iAqaiuuuOImI9dAFUX37kDv3rqrCL+8fydRUaugFFCq1NW23iiRkgLceCNw4YXcyBQoBkAis2kPgKmpsv6sZ08eqWGFU9d2pdt2ZMdujh6VUepVq3RXYg2v14v69fcjOjoHXbqstm34y5OTA/TrB1SsCMycqbsa+2MAJDKb1gDo9wN33w20bi1nd1H4ud1uREc/a/u1XXY0bhzQtKnuKqzVsqUcedOkie5KAvfll0DZssCwYfxSeSYMgERm0xoAX3sNqFsX2LdPy9MbadSoHxAVtS1fVwX7ru2ym5Ytgfff112Fta6+WjrDuFzA+vW6qwncH38A554LdO4sswx0KgZAIrNpC4BTp8ruvTVrLH9qY61bB5xzjh81akxFbGy80YfgBuvvv4HYWOk1a5JrrgEmTgQ6dZKdtpFk716ZXbjkEmDrVt3V2A8DIJHZtATA1auBcuWA776z9GmNtnSp7F4dNAjIyuIhuMHq3x+45x7dVVivTRvg00+l9V29epE3per1Ar16SS9gnjBwIgZAIrNZHgD37pWpmddft+wpjTdvnoy2jhypu5LIlJ0N1KgBzJmjuxLrXX89MH48kJ4u/4ZWrNBdUdGMGycbeEaOjLwQGy4MgERmszQAZmXJmqJ77uGbsFW++04++D77THclkeuHH4A6dcxsTXjjjcDYsfL/u3YFnn1Wbz3F8V2vAuYAACAASURBVNtvQM2aQLduQEaG7mr0YwAkMptlAdDvBx59FGjWDEhLC/vTEST0lSnDqfbiuuOOyFv/Firt2gGjR8v/nzULqF07sjtu7NwJXHklcMUV0tPZZAyARGazLACOGgWcfTaQnBz2pyLIVFf58lz3VFx79wKlSgGbN+uuRI/27WUXMCDr6apWlW4okSwzE3j4YXk/+uUX3dXowwBIZDZLAuBPP8lI1K+/hvVpCDLS+vLLsuFj6VLd1US+d96RjRCm6tDhxKNvHnsMeOIJffWEit8vnYfKlAHGjNFdjR4MgERmC3sA3LgRqFyZa9Cs4PNJm7JzzomsM9vsyu8HGjWSY1BMddttwHvvHf/vhQuBatVkY4wT5P15nnjCvMPoGQCJzBbWAHj4sPTmjOSF45HC6wUeeAA4/3yeeRYqS5fKNLrJa1bvvBMYMeL4f/t88gUjKUlfTaG2bZt0eLn6amDPHt3VWIcBkMhsYQuAOTnSQqpdOzN3T1opPV2m6po0AXbv1l2Nczz2mPSoNtlddwHDh5/4a4mJspPWSdLSgC5dZJNLpB51EywGQCKzhS0A9usHXHABcOhQyC9N+Rw5It0aWrYEUlJ0V+McaWlAhQrA77/rrkSvzp2Bt9468deWLZPXxmlHqfj98mctWxaYNEl3NeHHAEhktrAEwM8/BypVAv75J6SXpZPs2wc0bw7cdBP7nYbapEmy/s/08yq7dAHeeOPEX/P7ZanBtGl6agq3pCR5/+rb1zlrHQvCAEhktpAHwN9+k511JnZNsFJysqyvvPtu8xavW6Ft21OnPk10//3AkCGn/vqLL8q/PafatEm+ANxwg3P7PzMAEpktpAFwxw5pmcWWY+H1zz/STq9HD66vDIctW4CSJc3aEHA6Dz5YcNvGdeuA0qWBo0etr8kqR48Ct98OxMcDf/6pu5rQYwAkMlvIAmBamkxHPvIIp83CadUqoHp14Pnn+TqHy8svywc/Ad27A6++WvDvXXyx89fK+XzAK68A5coBU6boria0GACJzBaSAOj3A/feC1x1lZyyT+GxeDFQseKpa7IodHJypO/vjBm6K7GHhx+WQFyQIUOkU4gJpk+XjS8vvRTZrfDyYwAkMltIAuCQITIluXdviN6Z6BSzZskoxCef6K7E2ebOlWUMTl78H4wePST0FGTzZpkqd+oauZOtXw/UqydHLlnQPTPsGACJzFbsADh9ugSTVatC+M5EJ/B4ZGONx6O7Eue7916ZXifx2GPAgAGn//0WLcxqpZaSImebXngh8PffuqspHgZAIrMVKwD++ad0SnDa2hg7GT1aAvbs2borcb6DB4HYWGDDBt2V2McTTwD9+5/+90eMMK9Xck6OvCYVKwIzZ+qupugYAInMVuQAuH8/ULeuLJCm8Bg2TM4jW7JEdyVm+OADICFBdxX28tRTcqj76ezcCZQoAezaZV1NdjF5shwaPXRoZG7IYgAkMluRAqDXK90nOnVyzoJoO/H7ZRqyenVg9Wrd1ZijeXNg7FjdVdhLnz6F9/Ju08bco59WrpT1z3ffHXmHsTMAEpmtSAHw8ceBSy+NvDe8SJCTAzz6KHDeecDGjbqrMcfq1bLO0gmL+0PpmWfkdiaffCJrAU21bx9w7bXSi3vLFt3VBI4BkMhsQQfAjz4CzjoL2L49jO9MhsrMlJGEhg3lUG2yztNPA9266a7Cfp59VkYBz+TAAdkNHEnhJ9S8XqB3b6BqVWDePN3VBIYBkMhsAQVAr9cLj8eDLl3GIDY2GwsX8oyMUEtNBW68EbjsMllfSdbJzJQP7kWLdFdiP/36yTrAwrRvL2vhTDdhgowkjxhh/3WBDIBEZis0AHq9XiQktEGpUtdCqWMoUeJ5JCS0gdfrtfCtytlSUmTzwbXXcgpSh2+/lfPd7P6BrUP//rITuDBffCGdQUj6odesKSPKGRm6qzk9BkAisxUaAD0eD1yueCjlhVKAUqlwueLg4aF0IbF7t3xwduxo7w8LJ7v5ZjnMnE41YICcBViYo0cBl0t6BJPsir7ySuDyy+27nIMBkMhshQZAt9uNmJgeueFPbjExPeB2uy18q3KmLVuA888HunZl5wldduyQ9Wt2/ZDW7aWXZFNSIO6++/RdQ0yUmSm90c8+255HOTEAEpktiBHA1NwAmMERwBBYtw445xxZOM6jdPQZMkRGAKlgL78s/YADMXWqfKHhVPpxfj/w4YeyLtBubRwZAInMFvAaQJcrDlFRXyMqagXXABbT0qVAlSpyiDY/LPXx+yWwfPut7krs69VXge7dA7tvRgZQoQKwfHl4a4pECxfK6QmPPw5kZemuRjAAEpktqF3AiYmDULp0FhYt4nxlUc2bJ+3zTD04104WLZIgnpmpuxL7ev11WaIQqG7dgL59w1ZORNu+HWjaFGjVCtizR3c1DIBEpgv6HMB+/YA77gjju5KDffedtI76/HPdlRAgYaWwM+5MN2QIcP/9gd9/9mygVi0uazidtDSgSxegdm39I6UMgETOlaCUmqeUOqqUOqSUWlLAfYIOgDt2ALGx7FIRrM8+k/A3fbruSgiQXatlyrDVXmHeeEMCS6C8XqBaNZ6peCZ+P/D22/J+MHGivjoYAImcKUFJ6LtfKRWrlIpWSl1ewP2K1Aqua9fAzgYjMXKkrI2aP193JZRn3DigWTPdVdjfW28BnTsH95gnnpC1bnRmSUlA5cpAYqKeUwAYAImc6Wel1NsB3K9IAXDNGqB0aemBSafn98suyqpVgWXLdFdD+bVsCXzwge4q7O+dd4C77gruMYsWyb957hMr3KZNwEUXAdddJy31rMQASOQ8pZVSOUqpN5VSS5VSB5RSy5VSdxZw3yIFQEDalr3yShjelRzC55MjXmrVAtav110N5bdhgyxjOHhQdyX2N2IEcOedwT3G55N/97Nnh6cmpzl6VNZVx8XJl2urMAASOU8tpZRfKbVbKdVMyfTvHUqpLKVUi5PuW+QAOHeurPVJTw/DO1OE83qBBx6Q9mJbt+quhk72/PPAvffqriIyvPcecNttwT+ub1/gwQdDX49T+Xxy5E65ctYdS8QASOQ8FZQEwKEn/XqSUmpYAfctUgD0+4EmTYDRo8PwzhTB0tOBDh3ktdm9W3c1dDKvVzozzJ2ru5LI8MEH8u85WMuXy3FHbG8YnOnT5XUbMCD8O6kZAImcaZMKIgD26tULiYmJSExMRFJSUsBvIBMnAvXrAzk5YXyXiiBHjgDXXCPryw4d0l0NFWTGDODcc/lvNlAffQS0bx/84/x+GQGfOjX0NTnd+vXyvtq+PXD4cGivnZSU9P/v9b169WIAJHKgPkqpXUqpS5RSUUqpjkqpdKXUZSfdr8gjgICMptSqJefbmW7fPqB5c6BdOyA1VXc1dDq33y4bcygwo0fLv+mieOml4DeQkDh0SF73Cy4A/v47PM/BEUAi5+qvlEpWSh1RSq1QSt1awH2KFQABOc/qqqtC+K4UgZKTgQsvlOMy7NLmiU61Zw9QsiSwZYvuSiLH2LGy4aso1q8HXC7Z5EDBy8kB+veXI6R++CH012cAJDJbsQPg4cOyZuXXX0P4zhRB/vlHphR79uS0ot0NHw60aaO7isgyfjxw/fVFf3yTJsAXX4SuHhN5PHJo9ODBoe0dzgBIZLZiB0AAePZZoFOnEL0rRZBVq4Dq1eVbeijfmCn0/H6gUSOGkWB9+mnxQvPQocAtt4SsHGOtWgWcd55MqR87FpprMgASmS0kATA5Wc5V27QpNG9MkWDxYqBiRWmVRfa3dKlMpaWl6a4kskycKBubimrLFpl2t/qQYyfatw+49lrg4ouBv//2wuPxwO12w+PxwFuEU7cZAInMFpIACMi5d08+WezLRISZM2VKZswY3ZVQoB57TG4UnEmTgFatineNK68EPvkkNPWYzusFnnoqByVKHEbJkg8gJqYHXK54JCS0CToEMgASmS1kAXDVKmkPt39/sS9la5MnA2XKAF9/rbsSClRamoz+LV2qu5LI89VXcqxRcYwcKSNXFBoejwclSvSDUn4oBSiVCpcrDh6PJ6jrMAASmS1kARCQxeKvvRaSS9nSxx/LSf1scRVZvvhC+q1ynWbwvv5aRvCKY9cuoEQJYOfO0NRkOrfbjZiYHrnhT24xMT3gdruDug4DIJHZQhoAk5KAs85yXns4v18Ws1eqBCxZorsaClabNsA77+iuIjJ9+y1w+eXFv07bttJXmIrP4/HA5YqHUhkcASSiIgtpAPT7ZYGyk9b7+P1Av37SPmz1at3VULA2bwZKlQL27tVdSWSaOlUOOC+uMWOAK64o/nUI8Hq9SEhog1Kl7oZSB+FyxXENIBEFLaQBEAA+/xxo0CD8fSytkJMDPPqoHL+wcaPuaqgoBg4E7rhDdxWRa/p0oGnT4l/n4EHZDbx5c/GvRRICR4z4ETEx2Zg8mbuAiSh4IQ+AWVnAOecA338fsktqkZkJ3H030LAh8N9/uquhosjJAerUCU8XBVPMmCGHOYfCrbcCQ4aE5lokvceVAlJSivp4BkAik4U8AALAm28CV18d0ktaKjVV2l9ddpnzdzU72Zw5QI0aQHa27koi148/Ao0bh+ZakyaF7loky1NcLmDDhqI9ngGQyGxhCYCHDslu2d9/D+llLZGSAiQkyLEVIX5ZyGL33AM8/7zuKiLb7NnSQSUUjh2To6LWrg3N9UiWpyxcWLTHMgASmS0sARAA+vaVtkWRZPdu2cTSsSOQkaG7GiqOgwelO01RR0dIzJkDXHBB6K7XuTPw4ouhu57prrii6GeSMgASmS1sAXD7dtl9+e+/Ib90WGzZApx/PtC1K6cMneCDD4p/gDEBP/0E1K8fuutNmwbEx/NMxlDp0AEYNapoj2UAJDJb2AIgANx3H/DUU2G5dEitWycbV3r3dsbuZQKaNQPGjdNdReSbP1++GIVKRoZ0ZVm2LHTXNNmjjxZ9RJUBkMhsYQ2Af/whbdPs3Ah+6VKgShXglVc4KuEUq1bJv7ujR3VXEvkWLgTq1g3tNbt3BxITQ3tNU734ooTAomAAJDJbWAMgIB0AXn89bJcvlnnzgPLlpVcpOUefPhIyqPgWLwbOPTe010xKkhH3nJzQXtdEo0bJNHBRMAASmS3sAXD2bKB6dfttqvjuO6BsWTm4mpwjM1NGdH/+WXclzvDLL0Dt2qG9Zna2tIws6u5VOu7rr4veYYUBkMhsYQ+Afr+c/TV2bNieImiffSbhb/p03ZVQqH3zDVCvHqfzQ+W334CaNUN/3SefBB57LPTXNc3ChXIUTFEwABKZLewBEAA+/RS48EJ7bLAYOVIWoc+fr7sSCod27YChQ3VX4RzLlkkf7FD7+WegalWgCB3MKJ8NG+Qw6KJ84WEAJDKbJQEwM1NGEWbMCOvTnJHfD7z8snzocAeiMyUnS79Ztu4LnRUrZLo21Hw+mVqeNSv01zZJSoq0gyvKWzgDIJHZLAmAADBsGNC6ddifpkA+nxzxUqsWsH69nhoo/AYPBm6+WXcVzrJypaypDIfnnpNzN6no/H45b3XjxuAfywBIZDbLAmBee7ilS8P+VCfweoEHHpB1YVu3WvvcZB2/X86r+/Zb3ZU4y5o1QKVK4bn2ihXynpCeHp7rm6J2bdmtHSwGQCKzWRYAAeCZZ6QVlFXS0+WIhCZNpM0bOdeiRTK9n5mpuxJnWbtWjkoKB79fuoxMmRKe65uiefOivYYMgERmszQAbt0q0xVbtoT/uY4cAa65RtqBHToU/ucjvbp1A55+WncVzrN+veyYD5eXXwY6dQrf9U1wyy3Ahx8G/zgGQCKzWRoAAeDee2U9Xjjt2yffitu1A1JTw/tcpN+RI9L5Y/Vq3ZU4z99/A6VLh+/6f/0lu1gtfAtynIcekiAdLAZAIrNZHgBXrJAP64MHw3P95GQ5cqZzZyArKzzPQfYydqwEfgq9jRtl1D6cLrkEmDgxvM/hZG530c5UZAAkMpvlARAA2rQBhgwJ/XX/+UfaVvXsyTZTJklIKNoUGBXu33+BEiXC+xxvvMHd28UxYgRw++3BP44BkMhsWgLgzJlyuGwoF+yvWiUt5/r3ZxcIk+RNIaak6K7EmbZuBaKjw/8cJUsC+/eH93mc6quv5EtQsBgAicymJQD6fECjRsD48aG53uLFQMWKMpJAZvB6vfB4PGjRYiFattwGL1tKhMX27XLQcLi/VCUkAKNHh/c5nGrePCA+PvjHMQASmU1LAAQk/DVsWLT2cHkf/m63G/37L0LZsn6MGRP6GsmevF4vEhLaIDa2AZQ6jJIlH0BCQhuGwBDzer348MPvoRTw1VeesL6+o0bJrn0K3rp1RdupzQBIZDZtATAzE6hRA/jxx+Ael/fh73LFIypqDJRKQ/36r/DD3yAejwcuVzyUyoRSgFKpcLni4PF4dJfmGMdDdgsoBcTG1g9ryN69W9Yaso1f8Pbvl1HaYE88YAAkMpu2AAgAQ4cC114b3GOOf/hn5H74Z/DD3zButxvR0e/l/v3LLSamB9xut+7SHOP4z1mqZSH7uuuAd98N2+Udy+cDYmKAzZuDexwDIJHZtAbAgwdl6mL58sAf43a7ERPTgx/+Bhs06H9Q6li+LwEcAQw1CdlDLP05GzcOuPzysF3e0WrWBH79NbjHMAASmU1rAASAPn2Ae+4J/P6njgDyw98kf/wBVKzoR3z823C54hAT0wMuVxzXAIbQv/8CCQnboVQalMqy7OcsJUXOHPz337A9hWNdeinw3XfBPYYBkMhs2gPgli3ypr91a2D3z1ubVKpUWyjlRWxsPD/8DfHXX0C1asDw4SduBPJ4wrtBwRT798sXMpcL6N7dh2bN7rY8ZHfoAAweHNancKSbbgp+FzUDIJHZtAdAQLp2BNPH1ev1YtKkrxEVlYP335/BD38DbN0K1KoFvPSS7kqcJy1N1uNWqAC0bw+sXSu/riNkf/UVcNFFYX8ax3nwQeDVV4N7DAMgkdlsEQCXL5e1gMEe5luvHvDTT+Gpiexj1y7g/POlhzQP+Q6dnBxgwgQJ1s2bA/Pn664IOHZMeg/nhVAKTL9+wJNPBvcYBkAis9kiAAJyBtiwYcE95uab2QLM6Q4eBBo3lhGOopwZSafy+4FZs+R1jYsDJk+212t7zz3AgAG6q4gsw4cDnToF9xgGQCKz2SYA/vCDnAsYTHu4Pn2AZ54JX02k17FjQIsW0uc0O1t3Nc6wYgXQti1QpYr0kA1lO8ZQmT5dOltwtDdwX3wBtGoV3GMYAInMZpsA6PMBF14oU1KB+uAD4JZbwlcT6ZORIUHl+uvtGVIizZYtwH33yfRq//7AoUO6Kzq9zEzZ6T148Fxu8gnQ3LlA/frBPYYBkMhstgmAADB2rPQIDvSb/5w5wb/pkf15vcBttwFXXimjgFR0Bw4AiYmys7dbN+nta3derxdnnTULMTHjc3cgc6d/YdaskU08wWAAJDKbrQJgRgZw9tmyPikQW7fKCfj8XHAOnw944AGgSZPgNwXRcenpwJtvAhUrAu3aSUCIFB6PByVLdoVSPp71GaA9e+Sw7oyMwB/DAEhkNlsFQEDOAGvbNrD75uTIGYIbN4a3JrKG3w/06iW7u/fs0V1NZMrJAT7/HKhTB2jaNDJ3yUsXEje7/QQhJweIjga2bQv8MQyARGazXQA8cAAoU0Y6PgSiUSNg5szw1kTWePFFoHbt4D7E6Lg5c4BLLgHOOw+YNMleO3uDISOA3TkCGKTq1YGlSwO/PwMgkdlsFwAB4KmngC5dArvvbbfJbkaKbG+9BZx1FrBhg+5KIs+qVcANNwCVKslxIMFMA9qR1+tF3brvIzp6Llv9BeHii4EZMwK/PwMgkdlsGQA3b5ap3UBGgp57LvgDUMlePvlE1qqtXKm7ksiybRvQtavs7H3uOTkz0Sl69sxBhw5/cRdwEK6/XjbSBYoBkMhstgyAAHD33bJ7sTBjxsgbH0WmyZOBcuWAxYt1VxI5UlKk80Pp0rJhxolT5m3aAOPH664istx/f3B9lBkAicxm2wD4++8SDAo7r2zBAlnzRJHnxx9lvefs2boriQyZmcA77wCVKwPXXefsEdNzzgGWLNFdRWTp21faJQaKAZDIbLYNgABw9dXAG2+c+T47dwJRUZG/7sk0CxdK/+dvv9Vdif35fLKpo25dOR4nKcnZXTKOHJGdv/v3664ksrz5JtC5c+D3ZwAkMputA+D33wM1awJZWae/j98vQWLdOuvqouJZvlwOreUUX+H+9z85zqVOHTneJSdHd0Xht3w5ULWq7ioiz2efSU/1QDEAEpnN1gHQ5wMuuEDe2M7k0kuB776zpiYqnvXr5cOdO7fPbM0aOcC5YkUZ2UlP112RdSZNAlq21F1F5Jk9W9ppBooBkMhstg6AgGzyaNz4zFNed98tH5Jkb1u2yNquQYN0V2JfyclA9+7Sui0xUc7FNM3AgcBDD+muIvKsXCnrQwPFAEhkNtsHwIwMOeD0TBsFBgwAHn3UupooeDt3AvHxwNNPO3v9WlEdPgy43bKzt0sXCcum4he6otm5U9ZOnmnJTH4MgERms30ABIDXXpNdj6fz2WdA69bW1UPBOXAAuOgiGdmK1O4U4ZKZCYwcKdPibdrI+jfTNWkCTJ+uu4rI4/VKANyxI7D7MwASmS0iAuD+/XJcyOmOvfjlF9ksQvZz9Chw+eVAp05AdrbuauzD5wM8HiAuTpY4zJrFkVFAXpfSpdkRpqiqVgVWrAjsvgyARGaLiAAIAL16yUGnBdm/X775Hj1qbU10ZunpwLXXAjfeKCNdJBYskFBcqxYwYYIZO3sDtW0bUKKEjGZR8ILpjc4ASGS2iAmAGzZ4UaJEDp58ctgpraH8fumD6uSDcSON1wt06CC7OVNTdVdjD+vWAe3byxE4Q4cCaWm6K7IXr9eLF15YiCpV9rH9WxG1aSNfKgLBAEhktogIgF6vFwkJbRAdPQtRUXPgcsWf0hz+iitkSo308/mA++6T43kK6+Rigv/+Ax55RHb29ukD7NunuyL7yfsZL1HiVSi1usCfcSrcvfcCw4YFdl8GQCKzRUQA9Hg8cLnioVQ6lAKU8iI2thk8+RLf/fcDr7+usUgCIKOxjz8ONGgA7N2ruxq9jhwBXnxR1q927gxs2qS7IvvyeDyIjW0GpXJyf8ZT4XLFnfAzToV7+mngmWcCuy8DIJHZIiIAut1uxMT0yP1gyLul46qr/ofDh+U+r74KPPig3jpJjjKpUwfYvl13JfpkZQHvvw9Uqya703//XXdF9iZrRWdBqbQTfsZjYnrA7XbrLi+iDB0qxwgFggGQyGwREQCPjwCm/v/oQKlSd6Fhw72oUgV46y1pK3bllborNduwYcBZZwH//KO7Eut4vV54PB643W5MnuzB5MnZqFcPaNgQmDGDO3vPJCdH2tvVqQPExR1EyZL3n/AzzhHA4I0fD7RtG9h9GQCJzBYRATBvfZDLFYeYmB5wueKQkNAGWVlezJkDXHbZ8VGDQA9BpdD6+GPZiLNqle5KrHP832U8oqPfRFTUSpQsuR+jR2fzyJtCzJkDXHIJULcu8NVXQGZmwT/jXAMYnB9/lDM3A8EASGS2iAiAwIkjLQXtAv70UwmAFSvKwdA8WsM6X34JlCsn5zGaxOPxoFSpu6BUdu4XkCzExjbiqNUZrFolxwJVrgy8886JxwOd6WecArN8uSw9CAQDIJHZIiYABqJqVTkvMC5OpuCmTOEUXLjNmAGULQvMnau7Euv4/UBSElCnzuZ8G5O4bu1Mtm+XNbqlSwP9+gEpKborcqbkZCAqKrBD1xkAiczmqADYqhUwcaJMA3/8MXDOOUDz5tJHmEEw9ObPl/A3daruSqzh88mXimbNZJTlnnvWIDb2Eq5bO4NDh4Dnn5fg17WrHPRM4ZOZKV9Edu0q/L4MgERmc1QAfOghYODA4/+dng4MHy4jg1dfDfz8s77anGbpUqB8eZl6dzqvV5YVXHihdO8YOVIOtz7d2lROXUoQefddoEoV6ePNQ9qtU7lyYGtxGQCJzOaoADhsmByEerIjR4BXXpHA0q5d4L0yqWBr18oH+8iRuisJr/R04IMPgHPPBerVA8aNO7WlHdetncjnAyZPlmUYF18sU+UcfbfWhRfK614YBkAi5/tOKeVXSrUt4PccFQDzpudO58ABWX9UujRw113AX39ZV5tT/PsvULOmnLvoVEeOAG+8AVSvLjtVPR5uKgrEggWyI792bW7E0umaa+R4ncIwABI524NKqSSllE8ZEADXrJFRvsJGHHbuBJ54QlpzdesGbNliSXkR77//ZGSnb19njurs2we89JLsJG/ZEpg505l/zuI6edRz1Srv//c4HjZMRk5Jn86d5WzUwjAAEjlXbaXUttz/NWIEMC1NFkDv3h3Y/Tdvlp2JLpfsHg5k4bSp9u+XndWPPOK8UJScLC20ypQBbroJWLTIeX/GUDnx7MPnEBPjQVRUJnr3zsH+/bqrIwDo3Rt49tnC78cASORcc5RSj+T+fyMCICDTT8Fu9li/HujUSQLA88/LVDEdd+SI7Ka++25nTett3CiBNjZW/v65NrRwx7vy5LVty0apUtdw57ONDB4MPPBA4fdjACRypieVBMA8xgTAtm1lsX5RLF8uI0AVKsgat6NHQ1tbJEpPl3627do5p8vK6tXAPfdI8OveHdiwQXdFkaOgvtzR0U/z7EMbGTsWuOGGwu/HAEjkPPFKqV1KqTr5fu2MAbBXr15ITExEYmIikgLZPmZjjz0G9O9fvGssWiRnClarJsfImLqmKSsLuOUWeS3S0nRXU3xLlgDt28smoKee4pl0RXFqX+5sREWtxKefTtFdGuWaNi0b5557qMCd6UlJSf//Xt+rVy8GQCKH6aaUylRK7VNK7c+9+ZVSh5RSo0+6r+NGo5CsRwAAIABJREFUAN95B7jzzuJfx+8HZs0CmjaVA6VHj5bz4EyRkyOjZE2bAocP666m6PK6drRuLRuE3G5gzx7dVUWuk88+jI29ABUrLsPVV/sc8SUh0nm9XjRu/ASU2pd7NmX8ac+m5AggkfO4lFLnnHTzK6XuVkpVOum+jguAM2YAjRuH7no+H/Dtt3K2Vnw88MUXzloHVxC/H+jZE7jgAmDvXt3VFM3JXTsGD5auFFR8J+8CPnzYi9atZdoxI0N3dWbzeDyIjW0FpfyFdqdhACQygxHHwACynsvlkgAQStnZcrbZeecBF10ETJvmzJ2ifr+clXjeecCOHbqrCV7+rh21awPvvSddOyi8jh4FWrQAOnQwa6TcbtxuN6Kj3w2oPzUDIJHZHBcAU1O9iIry4YknhoWlM0NmpnSHqFEDuPxyYO5cZwXBIUOAs8+WHbKRJD0deP996dpRvz4wfrx1m1bYDUSkpACXXiq7xbOzdVdjpk8/nQKldkCpTI4AEtEZOSoA5q1Piorahujod864/qW40tKAN9+UvpvXXgv88kvIn8JyH3wgf541a3RXErjDh+Xw4byuHV9/be0Uff5z8Qpbc2WCffuARo3kfM1Qj8JT4Xr0yEGFCisRGxtfaH9qBkAiszkqAB7foZhd6LffUDl8GBg4EChXTnaYBtKE3Y6++EI2Sfz2m+5KArNvH/Dii9K146qr9HXtOP5vLt2yf3N2t2uX9E5+/HFnjY7b3Zw58jO8cWNgI9IMgERmc1QALPiMst6WnFG2b5+0SCtdWlox/f132J8yZKZPl0Owf/pJdyWFy9+1o1274A/9DjVZczUgoDVXJtm+XdaRJiYyBFrh8GFZ8/rJJ4E/hgGQyGyOCoCnnlGWA6UOoFu3PyxbD7Zjh5xFGBsLPPyw/c+a+9//gLJlge++013Jmf3zj7yesbHAXXcBf/yhuyLh8XgQHT03oF2Xpvn3XzlC6aWXdFfifA8/LLuwgwnbDIBEZnNUADz1jLI4XHihGxdd5EdcHPDVV9atS9q0Cbj/ftmR3Lt34P2JrfTbbzJlNHGi7kpOb9UqGVG1a9eOH37IRkzMUZQq1QlKeREbe5HRawBPtmEDcNZZsrmIwmPmTOlelJwc3OMYAInM5qgACBS8IzMnB/j0U5kiadbM2qnOtWuB22+XUbYXXpCdknawZo1s+Hj/fd2VFGzJEulCYueuHVlZQIMGwMiROfB4PKhceT+ee+5nhr+T5P1bGzFCdyXOk5Iio6zjxwf/WAZAIrM5LgCeSXo68NZbQKVKMl2ycqV1z710qTxnxYpyKPGxY9Y998k2bZKjXgYP1ldDQfJ37ahQQQKznbt2vP22HDqed+TJ008DPXrorcmuli2Tv9PRo3VX4ixdu8oXpaKss2QAJDKbUQEwT0qKHHZcujRw333Ali3WPff8+UBCgkyLjRhhfeeEHTtkcf5zz9lncX5OjnRbadpUunYMGWL/rh27dsn0+fz5x39t7lygVi37vK52s3ix7Jb//HPdlTjD9OnyZfa//4r2eAZAIrMZGQDzJCcDDz0k6/Seflp28lrB7wd+/FHOratdGxg71pruCfv2SYeMRx+1R0jxemVq/oILjnftiJR+sg8+KAce55eZKVP9kXoUkBXyNh19843uSiLbgQMyil+c9bsMgERmMzoA5lm7Frj1VpmiGjxYWodZ0d3B5wM8HulcUa9eeDepHD4s6x87d9bfyzh/144GDazt2hEKv/4qx9Bs337q791+u/2m1u1m5kx5/X74QXclkatLF6Bjx+J9kWMAJDIbA2A+CxdKP9OaNf2Ijx+O2Nj6lnR3yM6WEFSnDtCkCfD996EdoUtLA1q1Am6+WW/QOnwYGDpUunZceqmMAukOo8HKyQGaNwdee63g3x87Vqb46cymTJGRwLlzdVcSeaZMAapUKf7JAgyARGZjADyJ3w8kJi5BVNRmKOWz9Gy3zEyZBq1eXYLovHnFv2ZWlhyY3Lq1vunVffuAAQNkhPWqq4BZs+wxBV0UY8cCdevKKGZBdu4EYmKA/futrSsSTZokawIXLdJdSeTYu1fWyU6eXPxrMQASmY0BsADS3eHxkzqK9LGsu0NqqoyUVaoEtG1b9PZsOTmyTq15c0DHX3FyMtCnj326dhTXoUPy4VvYodlNm0prPSrcmDHyxWDJkuywL7mIdH4/0KmT3ELxBYoBkMhsDIAFOLWjSDaU2o1+/awdqjh0SPrdli0r633WrAn8sX4/8MgjQMOG1o9G2bVrR3H16RNYt4WBA4F777WmJid4992c3MO021uy5CJSeTzyBWTv3tBcjwGQyGwMgAUoqKPI+ee/gUqV/Lj/ftmBZ6U9e2SXsssli783bjzz/f1+6Utct27Rj4goipUrZcTRrl07imPtWnn9//qr8Pv+/ruM3uadD0hn5vF4UKLEG2yndwa7d8u6vylTQndNBkAiszEAnkZBu4B37pSRuOrVQ/tGHKjt2+UIl9hY+d/TtX56/XWgRg058NkKixfLBpPSpaXtXUG7YyOZ3y9T8X37Bnb/zEwvKlTIwP33j+Z0ZgDcbjdiYnqcsOQiJqaHZUsu7M7vB267LfSjygyARGZjAAyS3y/HtVStKtOboZqOCcY//8iHgcsFPPMM8N9/x8Nqt25/oEoVP/78M7w1+P3A7NnA1Vcf79qh47WwwpQpEvoPHy78vnmjx9HRUxEVNYvTmQE4vuQihyOABfjiCznzL9QzDwyARGZjACyiPXtkurNqVeDLL/Xsal29Gmjf3ofo6HTExHyAqKhJUOooGjd+ImyBIxK7dhRHWpqcVzhhQmD3Px5mMhhmAuT1enH55bdCqRxERz8HlyuOoTnXzp2ynGD69NBfmwGQyGwMgMWUNzrUsaO0B7Oax+NBqVJ35hs9SQ9L4MjKkhDUoIGcVxhJXTuKY9Ag4IorAj+gm9OZRTNjRjaqVz/GXcD5+P3S57dr1/BcnwGQyGwMgCFw4ABw//3yTf3TT60dDQx34EhLA0aNktDXoIGEwEjq2lEcW7fKETZLlwb+mFN3kIcnkDvNCy8A3brprsJeJkwAataU3uXhwABIZDYGwBCaMQM45xw58+50GzRC7dTAEZopx7yuHWedFbldO4rrzjulV3Qw8u8gV2oDSpToz+nMALRqJd1wSCQny9raH38M33MwABKZjQEwxA4dktBQoQLwySfhHw08+ciaYNZPFbTTee/e4107WrWK7K4dxfHTT/Ia7NkT/GPzXtfLL1+MG2/cyPBXiIwMoFSpwo83MoXfD9x4Y/BfPoLFAEhkNgbAMElKkmnTtm2BLVvC+1wFBblAHiPBMT73rMOrUKPGFJQu7cfNN0d+147i8HqBRo2Ad98t3nU+/1xCNJ3ZokWyy9XELxoFGTMGqF07sF3nxcEASGQ2BsAwOnIEePxx6Xc6alTgGwmscHzqOC136tiP6OgfMWxYku7StBsxQjqoFHfgbvVqoHx5e/2929HgwbKjnmTdafnywJw54X8uBkAiszEAWmDePCAuTs7Ms8s0V0H9jrlbVc4yrFABmDu3+NfKygJKlgT+/bf413Kym26SL0im8/lkxqBnT2uejwGQyGwMgBZJTZVesmXKAMOH699Q4fF4ULJkbyjl43l1+TzyCHDHHaG73qWXAlOnhu56TpOdLSPkq1bprkSP/Ms3Hn54Bc47z4+jR615bgZAIrMxAFps8WKgfn2gRQtg/Xp9dXi9XpQvvwYlSrwe9OYRp1q2TNrZhXLNZvfuwMCBobue06xYISOuur8Q6ZB/HW509AtQ6hgaNepr2c8gAyCR2RgANUhPB557TsLG0KEyCmK1VauAMmX8GD9+Kg/fhUy/tWgR+rA2YgTQoUNor+kkI0bIYccmOvUIJ6+lo/AMgERmYwDU6PffZbdps2bAmjXWPvfDDwOPPWbtc9rZZ5/Jru1QdzdZsEBayVHB7rwTGDZMdxV66O4awwBIZDYGQM0yM4EXX5TRwEGDrOmyceAA4HIBf/4Z/ueKBEeOyDEk33wT+munpMgH+8GDob92pPP75aDxJUt0V6LHqSOAGRwBJCLLMADaxMqVwCWXABdfLOuiwunNN4Frrw3vc0SSZ58F2rQJ3zl0554LzJ8fnmtHsg0bgNhY+RJkovyHuEdF/YyYmI8tXYfLAEhkNgZAG8nKAl57TUYDX3hBOiSEWk4OcN55wJQpob92JNqwIfyjoR07Fv9QaScaMwa45hrdVeiVtwv45pu/xYUX7rV0HS4DIJHZGABt6M8/gcsuk8OIf/sttNeePl26DOjYeGI3eS23evcO7/O8/DLQrVt4nyMSPfAA8NJLuquwh7wvIlYsAcnDAEhkNgZAm8rOBt54Q84NfPbZ0G1OuP56YMiQ0Fwr0k2fDlSrJuv0wmnqVJnepxOdd541HS8igc8HVK4sRxFZhQGQyGwMgDb311/AlVcC9eoVvz/vX3/Jmqt9+0JTWyTLyJDuLGPGhP+5Nm+WjiBWju7YXXIyEBMDyw49jgS33AKMHGnd8zEAEpmNATAC5OTIGrKyZYGnngKOHSvadXr1Ah58MLS1RarBg+X4HSsOIPb55LDj1avD/1yR4ssvZZkDHTd4MNC5s3XPxwBIZDYGwAiyaRPQujVQty7wv/8F99gjR6Tl1vLl4aktkiQnS5j+5RfrnrNVKzlrkMRjjwGJibqrsJf58+UsSqswABKZjQEwwvh8wIcfSpjr2VOCXSBGjZKpZALuuQfo2tXa53zqKeCZZ6x9Tjtr1AiYNk13FfZy7BgQHS1fUKzAAEhkNgbACLV1q2zoqF0bmD37zPf1+YAGDYBJkywpzdYWLpTwvHOntc87bpycNUhyEHlUFNeiFqRpU+Drr615LgZAIrMxAEYwvx8YO1bWl3XvfvrdrElJ0unC9E0I2dlAkyZyELbVli+XXZ7hOmw6kkyfLkcc0al69bJupJgBkMhsDIAOkJwM3HwzULMm8P33p/7+rbcC/9fevcfJXLZ/AL92drOTY1FUKLsKUdEB7eNJKZQOdHpCPJVCQgdSrUdSOatUDtGBjcIkRZRW/ZSSx/mQKFHO63xcaw+zs/P5/XHxrGV37ezOzP2d7/fzfr3mVXbHfC9j13z2nvu+rv79w1+X1YwZA1xxhZnJE+npeup127bwX9tqevfW7Qt0pk8/BRo1Cs+1GACJnI0B0Cb8fmDSJF1lat8e2L9fP/7330CpUuF/y9Nq9u/X52buXHM11KuXf0B3moYNgU8+MV2FNW3eDMTE6A8MocYASORsDIA2s2sXcO+9wIUXAtOnaxPpcLaWsKpu3XQl1KQOHXTUn5MdO8aV0MJkZXlx3nnp6NBhHDweT0hHwzEAEjkbA6AN+f2Ax6NTREQ493fVKp2vvGmT2TpGjADuv99sDaZ9/z1w6aWmq7Amr9eLhIRmcLmSERU1A253PBISmoUsBDIAEjkbA6CNDR6sAbBiRd1b5MQDCH4/0KQJ0Lev6UqA774D4uNNV2FW//66Ekpn8ng8cLvjIZIFEUAkDW53HDweT0iuxwBI5GwMgDbl9wNXX63tR778Uk8B33MPsHOn6crCa8oU4JJLij89JZj27tVA7uRvt1tuAcaPN12FNSUmJsLlGnoi/OktOroLEhMTQ3I9BkAiZ2MAtKmfftJDD8eP668PHNDmxxUqABMnOmM18NgxDX9TppiuJNfFFwMLF5quwoysLMDt1pnUlJffD3TuvBwixyGSyRVAIgo5BkCbevBB4IUXzvz4nDkailq2tP9G/MREHcFmpbDbqhUwerTpKsxYtAi44AJr/X1YQUYG8MQTQOXKftSr9yzc7jhER3eB2x3HPYBEFDIMgDa0Y4e2ftm8Of/PHz6sLzjlygHjxumkELvZuFEPfqxebbqSvPr2BTp3Nl2FGcOG6Ql1yrV9u7bFadxYv2+9Xi88Hg8SExN5CpiIQooB0Ib69QNatz77/b77Tk9kNmum/QLt5O67tfWL1Xg8wA03mK7CjLvuAkaONF2FdSxYoO2aunQx05ycAZDI2RgAbSYzU19Uvv++aPdPTQW6dwfKlAHefdceq4HffKP7Hw8cMF3JmTZsAGJjdSydk/h8uv90+XLTlZjn9wNvv61tmj74wFwdDIBEzsYAaDOTJwN16gS+z+rHH7VFSZMmwJ9/hqS0sMjM1HFvY8eariR/Pp++8K9bZ7qS8FqzBihb1nnB93THjwMPPwxUrQosXmy2FgZAImdjALSZhg115m1xpKXpIPrSpbVpsc8X3j1JwTB8OHDNNdYOGo0bW+tkcjiMHg20aGG6CrM2bwYaNABuugnYs8d0NQyARE7HAGgjS5fqwY7U1JI9zi+/ALVqAQ0b5qB+/U5wu+NPnEoM7WSCkkpJ0VWmn34yXUnhnnwy/xPadvbQQ8DAgaarMGfePG3I/vTTgFW+fRgAiZyNAdBGOnbUF5hgSE8H7rnnd4ikh20yQUl17Ai0a2e6irMbN85Zq2F+v/Y/XLDAdCXh5/cDQ4fqHtvJk01XkxcDIJGzMQDaxJ49erhgw4bgPaZOJhgUtskEJfHLL/oiu3276UrObvFiPajjlH54mzYB55yjP1Q4SWqq9uO87DJg5UrT1ZyJAZDI2RgAbWLgQG3uHEwejwexsU0g4rf0CqDPB1x7LTBokOlKiiYtDYiKAnbtMl1JeEycCPzjH6arCK8//wTq1gWaNwf27zddTf4YAImcjQHQBrxePVU4e3awH9eLuLi3ERW1JCyTCYrr/ff1BHNGhulKiq5WLWDuXNNVhEenTjqVxSlmz9aWNy++aO3DSAyARM7GAGgD06cDcXG6EhZsrVrloEOH1ZY9BXzoEFCpEvDVV6YrCcxDD+neMCe4/HLtzWh3OTnAgAF6EOmzz0xXc3YMgETOxgBoA02bAm+8EfzHTUvTfYV//BH8xw6Wnj2B22+PvP10gwcDbduariL0du3St7sPHzZdSWgdPqzTZ2rWBH77zXQ1RcMASORsDIAR7tdfdebtwYPBf+zZs/UFzarhau1awO0O7sGXcPnmG6B2bdNVhN5nnwH165uuIrTWrdPm43feGVlBlwGQyNkYACNcly5A586heeyuXYFnnw3NY5eU3w/ccgvQp4/pSoonJUVXxtLSTFcSWj17Bq81kRV9/rn23nzllcgbo8gASORsDIAR7OBBXf379dfgP7bfD1xySdFnCofb9OnARRcBkfql6/cDF1wALFliupLQuuYa/buyg1On4kyd+hn69PGhfHlg1izTlRUPAyCRszEARrA33tD9f6GwcqWubGRlhebxSyItDaheHZg0yXQlJdO8OTB+vOkqQufwYV3l3L3bdCUl5/V6kZDQDG53PFyu5+By/Yxzz92C336z1qGoQDAAEjkbA2CE8vn05G+oVldee02b2FpR//7AjTdG3ltup+vTB+jWzXQVofP117o3zg48Hg/c7vgTk3EAkWzExl5luZ6YgWAAJHI2BsAINXu29v4LVVeWhg2Bjz8OzWOXxObN+rb38uWmKym5pKRsXHHFfsu22Cmpl14CHn/cdBXB0avXAERFzY2IqThFxQBI5GwMgBGqZUud/hEKu3cD0dHA3r2hefySuPfe0B16CSev14v69TtBJA0uV1e43fGWbLJdEgkJQFKS6SpKxucDPvwQqFAhAy7X4lNWAK05FScQDIBE9jNURNaKyFERSRGRqSJSrYD7MgBGoA0btD9fqALahAn6FqvVzJunExasGEwDpWP2all+zF5xpafr/N+//jJdSfH98IO2sKlZE/j882zceGMzuN1xlp6KEwgGQCL7GSwi14pIjOg39hQRWV3AfRkAI9DTTwMdO4bu8e+7z3pzdb1eoE4d4J13TFcSHImJiYiO7mKrtxRP9eOPeorcqj0kC7Npk640ly+vB60yM/Xjp54CtsNb9gyARPZXX0RyRKRCPp9jAIwgXq8XSUkzUKpUBgYN+i4kL0AZGUCZMsCaNUF/6BJ56y2gXr3Q7XkMtzMPFdhrBfD11yNv0smRI3owx+0GnnzSHivNhWEAJLK/F0VkcwGfYwCMECfbUMTEvAKRzSHbM5acrC1WrLRys3u3rsb83/+ZriR4Tv59xsbGQWQ7YmL6Rvxbiqdq3hwYM8Z0FUWTnQ2MG6d9GW+7TSfMOAEDIJG9NReRYyLSooDPMwBGCN0zdiNEckK6YtSzJ/DUU0F9yBLr1Al44AHTVQTfybcUW7X6HNWqHUFWlj3CX3a2riKHokF5sH33HXDVVdquZvZsa/3gE2oMgET2dbeIHBaR1oXcp7yIoEePHujVqxd69eqF5ORk0/8uUT6eemoYRPaGdM+Y3w/UqKFzaq1i6VKgdGlgyxbTlYRORoauPll16kqgli0DzjvP2n0a//wTuOcePVQ0cqQ1G56HQnJy8v/+re/RowcDIJENdRANf83Pcj+uAEaAjRuBSpXSEB09CSJpIVsBXLdOe+ylpwftIUskJ0f7Eb76qulKQq9fP+Duu01XERxvvWXdP8uhQ8Bzz+k+vx49gP37TVdkDlcAieynp4gcEpEmRbgvA6DFrVunM2979/aFvA3FsGG6KmIVEycCl11mnUAaSikpQKlSegI10rVpAwwfbrqKvLKzdU9ixYraQ3PdOtMVmccASGQ/fhHJEpHUE7djJ/6bXyBkALSwVauASpV0BczvD30biiZNgPffD+pDFtuRI0DlysCMGaYrCZ/27YFnnjFdRcnk5OjX7OLFpivJ9e23QN26QO3aur3BSfv8CsMASORsDIAWtXix7qMaMSI819u/X6d/7NwZnuudTa9eeiLTSS/WS5YA5coBkfztuH69biOwwp66338HWrUCzj8fGDXKPi2EgoUBkMjZGAAtaMECDQLhbKPxySfAddeF73qFWb9e92g58W26xo2Bt982XUXxjRsHNGtmtoYDB7RZututK6oHD5qtx6oYAImcjQHQYpKTtYXGhAnhvW7btsArr4T3mqfzer2YNs2DGjU2olWrP23TEy8QU6cC8fE6gzYSPfywua8jr1cnxZx/PnDnnboCSAVjACRyNgZAC5k1S8PftGnhva7Xq+0wli0L73Xz1qCNkc85pxtEjiI2tr6tGiMXlderI9S++sp0JcVTvXr4G3b7/cDXX+sev7p19YcoOjsGQCJnYwC0iKlTNfzNmhX+a//4I1Clitm+bdro+tqQN7qOBIMGAbfearqKwG3dCsTEAGlp4bvmunV6qrdSJWDsWD3tS0XDAEjkbAyAFjBhgoY/UysXzz8PPP64mWvn1vAyRDaGtNF1pNi3T/evRdpIssmTgUaNwnOt/fuB7t31eerdW/v7UWAYAImcjQHQsNGj9cDHggXmaqhVC5g509z1fT6gYcMdiIpaEdJG15HkiSeAzp1NVxGYLl30h4lQysrSRtMVKgCtW+tEDyoeBkAiZ2MANGj4cG31YrJn2p9/agPiY8fMXN/v15Wc2rX9uOGG1iFtdB1J1q7VdiqRNKmiTp3Q7V30+/WxL79cZ/faZWyeSQyARM7GAGiA3w8MGKDzX1etMlvLyJHA7bebu/7gwXroYevW0De6jjTNmgFDhpiuomj27gWiorQFS7D9+qv2hLzwQmD8eO7zCxYGQCJnYwAMM78f6NNHx7utX2+6Gj1sMHq0mWsnJelbeb/+aub6VjdrFlC1amQ0MP7iC6BeveA+5t69QNeuuhL6wgs6HYaChwGQyNkYAMMoJ0ff7qxeHdi40XQ1+oIaEwNs3hz+a3/zDVC6NPDDD+G/dqTw+YC4OCAStkE+9xzQrVtwHiszUyfglC8P3HefPeYjWxEDIJGzMQCGic8HPPYYULOmvt1pBdOnB3/VpiiWLgXKlo2MYGPayJFAQoLpKs7u+uu1lVFJ+P3Al19qI+z69bU9EYUOAyCRszEAhoHXq5M2rrwSSEkxXU2uf/8bCHeXlY0bde9jJI87C6cjRzQsm2zSfTZHjwIuF7BjR/EfY9Uq4OabtR/lRx9F7iSUSMIASORsDIAhlpkJtGkDNGig/d2swufT5rm//BK+a+7Zo29p9ukTvmvaQc+eQIcOpqsoWHIyUKNG8X7v7t3a8ubcc/WHEf5TFD4MgETOxgAYQseP65SCxo2t16h20SKgYsXwrbSkpgLXXadBxuTEkUh0slXPrl2mK8lfv37AI48E9nsyMvSEc7lywIMPmtmH6nQMgETOxgAYIqmpQNOmektNNV3Nmfr2BTp2DM+1srI0CDdvrv9PgbvzTqB/f9NV5K9pU+DDD4t2X79f957WqKE/EPz0U2hro4IxABI5GwNgCBw6pCOxWrbUVUAruvrq8BzC8Pt1r2GDBnx7ryTmzdM+eBkZpivJKzMTiI0FNmw4+31XrAD++U9tgZSUxJVg0xgAiZyNATDI9u3TsNO6tb44WtHWrdr+5fDh0F/rpZd0tWf37tBfy878fj1ElJRkupK8Fi4EKlfW+gqSkgI8+qi2/enXz9zUGcqLAZDI2RgAgyglRV+k27a1dvPesWN1ykSojRqlB004rzU4xo3T9iiFha1wGzIEeOCB/D+Xng4MHAiUKQO0a2ed9kekGACJnI0BMEi2btUef489Zv0WFnfcAbz1VmivMX26ti9ZsiS013GStDSdHb1ggelKcrVqBbzzTt4xftOmefDJJ9moXh1o2DC8J82p6BgAiZyNATAINm3S6R7du1t/X9OxY3qiNJSrcgsW6Nt9c+aE7hpO9cILOh3DCnw+ndaxdKkXCQnN4HbHw+UagqiolShVah+SkrIt//3gZAyARM7GAFhC69cDF1+sve2s9NZcQWbNAq64InSPv3atrlIV9VQoBWbbNg3wW7aYrkSbN5crB0yd+hliY6+BiBcigEgWYmPrwsNRL5bGAEjkbAyAJbBqlU61GDAgMsIfAHTuDPTqFZrH3r4dqFoVeO210Dw+qQcfBJ5/3nTXafK4AAATDklEQVQVwLvv6naC7t2HQmT7ifCnt+joLkgM95gZCggDIJGzMQAW05IlutI1YoTpSoouJ0dXK+fPD/5jHzwI1K0LdO0aOWE4Ui1cCFSoYP407YMP6gGQ889PR3T0FIiknQiAaXC747gCaHEMgETOxgBYDAsW6FtfY8aYriQwK1bonq1gn1BOT9f+bq1bA9nZwX1sOpPfr02Ux441W8PJ1b5hw3y48cZmcLvjEB3dBW53HBISmsFr5aPwxABI5HAMgAFKTtbTrRMmmK4kcAMGAA89FNzH9Pn0UEJCgnWbXtvRpElA7dpmDh35/XrgSQSYOlU/duopYI/Hw/AXARgAiZyNATAAs2bp6daTL3qR5vrrgcmTg/d4J4NAnTrAgQPBe1w6u8xMoEoV4Ntvw3vd7GygZ08Nf253eK9NwcUASORsDIBFNG2ahr+ZM01XEjiv14v33puFqCgfPvjgy6CtzgwerHsK2eDXjAED9BBGuKSm6kziunWBm28G/vOf8F2bgo8BkMjZGACLYOJEnWaQnGy6ksB5vdqjLSYmESKb4HbHB2V/VlKS7idcsyY4dVLgdu/WObx//BH6a+3YoVNImjfXEYJxceFffaTgYgAkcjYGwLMYM0YPfPz4o+lKisfj8cDtjj+lR1sGYmPjS3RCc+5cXQ0NxWliCsy//61vw4fSqlXAJZdoCyGvF9i5E3C5AP6zEdkYAImcjQGwECNGaKuXxYtNV1J8iYmJiI7ukqdHm8gu3HHHF0hPD/zxli3TQzDTpgW/VgrcihW6On34cGgef84c/QFo+PDc9j7TpukpZIpsDIBEzsYAmA+/H3j1VW3yvGqV6WpKJncFMLdHW0xMH1SvfhgXXqj7yPbuLdpjbdoEXHghMHJkSEumADVpArz5ZvAfd9QoDZeff5734089BTz7bPCvR+HFAEjkbAyAp/H7dd7qRRfpmLdId3IP4Ok92rKyvJg3D2jZUk9zdulS+F6yPXuA+HhrTKCgvKZPB2rU0JY8weDzAc88o2E/v9Xvq64CZswIzrXIHAZAImdjADxFTg7QowdQvTqwcaPpaoLnbD3a1q4FHntMDxTcdRfwww8ahE/+vt69X0F8/EG0a5djpO8cFS47W79mv/yy5I917Bhwzz3a2mfz5jM/f/AgEBVV9FVjsi4GQCJnYwA8wecDOnUCatZ0bluTXbuAfv2AihWBBg38uPzyQYiNrQuR3+By/YLGjVuwwa9FDRumrVlKIiUFuPZa4NZbgUOH8r/P7NnagJoiHwMgkbMxAEJPNrZrB1x5pb4IOl1aGvD44ysQFbXllIMjnO9qZQcP6sns1auL9/vXrAGqVdMfgrKyCr5fnz56GpgiHwMgkbM5PgBmZgJt2miPs337TFdjHYmJiXC5uuY5PRwd3QWJiYmmS6MCPPmkvpUfqG++0ZO+Q4bknvQtSOPGOoaOIh8DIJGzOToAHj+uhyAaNSr4LS+nyu/0MFcArW39et3HGcj+vLFj9aRvUf5a09KAmBhgy5Zil0gWwgBI5GyODYCpqUDTpsBNN7GhbX4KOj3MPYDBdbYDOoFq0QJ4/fWz38/nA3r10lZHixYV7bHnz9e3ic+2SkiRgQGQyNkcGQAPHdK3slq00FVAyl+wwwnllRuy40+E7JKP6fv6a21hVNg+vrQ03fZQuzbw119Ff+wBA4D27YtdGlkMAyCRszkuAO7bBzRoALRuDWRkmK6GnGz8+Jk455yep4zp8yImZiCef34h1qwp3sp0Tg5wxRXAp5/m//ldu4DrrwduuSXwbQ+33gqMGxd4TWRNDIBEzuaoALhrF1C3LtC2rZ78JQqnfft0qkb37vp1GBWVA5Ftp43pW4nKlVNQvrz+ulIl4IYbgH/9C3jxRQ1g8+Zpn8qCVvlGjQIaNszbyzExMRHDh3+L6tX9ePTRwlcI85OVBZx7LrBuXYmfBrIIBkAiZ3NMANy2Dbj8cj0lGayJCUSF2b9fJ2b07KnTM1wu4JprdMrGzJnARx99WeBBG79fW7usWKGhccQIHcF2++1ArVpAqVLakLlaNd3H+sgjOr5w0iQ91SsCLFiQ/b+3mF2udyCSiurVJyArK/CffhYv1v6QbARuHwyARM7miAC4aRNw6aW68sIXMAqVAwd0GsczzwBXX60B7aqrgKefBr74Qj9/qpIctMnJAXbuBBYuBCZP1vD36KN6sKl69dwVxaiozRDxnfh1RrFPco8YodsmyD4YAImczfYBcP164OKLtYEtTy9SMB06BMyaBTz7rPaRdLmAevV0nOCMGUXrKxmqgza//34yACaVqJfjyfouv3w9OnZczYNANsIASORstg6Aq1drm4sBAxj+qOQOHwa++krbpzRooIHvyit1ZXn6dOvNx73++p2Ijh5b7F6OJ1coY2NrQiQNpUrdy1ZANsIASORstg2AS5YA550HDB9uuhKKVEeOAHPmAL17A9ddp4Gvdm2gWzdtnLxnj+kKCzdvXjZiYo4gNrZOsXo55jYDT2czcBtiACRyNlsGwAULdLTV6NGmK6FIcvSoHqDo00dP3rpceuCia1dg2jQ9RR5J/H6gXj0/unZdVqy3mHUcYG+OA7QpBkAiZ7NdAExO1tFWEyaYroSsLjUVmDtX26s0bAhER+tJ8c6dgSlT9JBFpPvwQz2IUpwtEFOnfgaX62eIZHMF0IYYAImczVYBcNYsDX9Tp5quhKzo2DH9AeGll3QSTHQ0ULMm8MQTwCefADt2mK4w+NLTtZfg/PmB/94BA3xwu3cgNvZqjgO0IQZAImezTQD0eDT8zZxpuhKyirQ0bZrcty+QkADExABxcUCnTto6Zds20xWGx3/+E3gLl/nzgdKlgWXLOA7QrhgAiZzNFgEwKUnDX3Ky6UrIpOPHge+/B/r1A/7xDw18l12mzb8//hjYutV0hWbs3KmNo4s693f3bqBKFeCDD0JbF5nFAEjkbBEfAMeO1QMfP/5ouhIKt/R0Xal6+WXgn/8EzjlHmyA/8ggwcSKwZYvpCq2jXTvguefOfj+fT+cEd+jA1kl2xwBI5GwRHQDfeENbvSxebLoSCoeMDOCHH4BXXtHxZ6VKAVWrAh076qGfv/9maCnI4sX6g9LZvtVffhmoU0f3S5K9MQASOVtEBkC/H3jtNd3cvnKl6WooVDIytKXPq68CN98MxMYCl1yiq1Mffqgj/hj4iq5RI+Dddwv+/MkT9L/9Fr6ayBwGQCJni7gA6Pdr246LLgLWrTNdDQVTZibw888a7ps1A9xu/Xtu3x54/31g40YGvpKYMkXb3OQ3D3vHDp2ak5QU9rLIEAZAImeLqACYkwP07Kn7vDZuNF0NlVRWFrBwITBwIHDrrRr4qlQB2rYFxo8HNmxg4AumrCydiz1nTt6Pe726h7JTJzN1kRkMgETOFjEB0OcDHn8ciI937mnOSJeVBSxaBAwaBDRvDpx7LlC5MvDQQ8B77wF//MHAF2oDB+pzf6oXX9Rm0cePm6mJzGAAJHK2iAiAXq+eYrzySiAlxXQ1VFRerx4+GDIEaNlS+8pdcAHw4IPAmDHA+vUMfOG2b5+utJ7cPjFnjh4O+eMPs3VR+DEAEjmb5QNgZibQpg1Qv76+eJF1ZWcDS5YAw4YBt9+uBwoqVQLuv1/nMv/2W/77zyi8Hn00B7fd9heeemoYypbNxKRJ2aZLIgMYAImczdIB8PhxDRKNGgGHDpmuhk6XnQ0sWwYMHw60agWULQtUrAjcd5+eNv31VwY+q/F6vbjmmicgkg6RLYiOnsLxbg7FAEjkbJYKgF5v7tippKQZaNo0BzfddPbeZRQePh+wYoX2X7zrLqB8ee3D2KYN8PbbwJo1DHxW5/F44HbHQ8QHEUAkDW53HDwej+nSKMwYAImczTIB0Ov1IiGhGdzueLhczyAqaiUqVFiOI0e4MmGKz6d9Ft96C7j7bg18FSoA99wDjBwJrFql96HIkZiYiOjoLifCn96io7sgMTHRdGkUZgyARM5mmQCYuzKRduKFKQelSrXBqFGzcfCgHiig0MrJAVav1nDXurWu7pUvr+HvzTd19Y+BL7Kd+X3GFUCnYgAkcjbLBMD8ViZEjkMk53+/dru1bUjNmsC11wJNm2o4ad8eePJJ4IUXgNdfB955R0eDff45MG8e8N//6qnHbduAw4cZYk7KydG3bd95B7j3XuD88/VE6J13AiNGAMuX6z4/so/clfY4REd3gdsdxz2ADsUASORslgmABa1MTJvmQVoasGuXNgZevhyYPx+YNQuYPBkYOxYYOhTo21ebRD/yiB5CuO02PTxSp46ODytXDnnCZenSOmXiiiuA668HbrlFV706dgSeegp46SXtVzdqlE5H+OIL4Pvv9ZTr77/r5ISjR62z5+3U/ZMejyffF/ScHGDtWv0z3XefHtgoUwa44w49yLF0KQOfExTla4XsjwGQyNksEwDDsTKRkwOkpgI7d2rfs6VLNdR9+SXw8cfaqmTwYA1/3btrGGzTRseS3XADUKuWhsYyZfKGybJlNWTWrg00bKhTLe69V8Nojx4aTocM0d53kyYBM2dqiF22TENtSgpw7Fjxe+Kdun9Sn7t4JCQ0Q1aWF+vW6Z/rgQe0B1+ZMnqyeuhQ7dHH134iZ2IAJHI2ywRAILJWJnw+4MgRYPt2bWi8eLG+3TxjBjBxorZBGThQpyx06wY8/LAenrj5ZuC663Qma5UqOg3jZJB0uXTPXbVq2vS6cWOd2nD//cBjjwFPPw3066d99t57D/j0U+Crr4D+/eejVKnWJ94yB0Qy4XJ9jfLlM1C6NNCihQbb//6XgY+IFAMgkX29JiIpInJMRBaISL187mOpAOhU2dna53DrVm2WvGgRkJwMTJ8OfPSRtlh57TWgTx+ga1fd83jXXcBNNwENGgDnnXcAIkfzrEpGRc1Ex47jkJVl+k9HRFbEAEhkTy+IyDYRqSsisSIyRER2ikjp0+7HAFgMycnJpkvIIxJOdlrtObM6Pl+B43MWGAZAInvaLCI9T/l1tIjsFZEOp92PAbAYevXqZbqEPCLhZKfVnjOr4/MVOD5ngWEAJLKf8iLiF5HGp318noi8mc99GQADZMUXGqvvn7Tic2ZlfL4Cx+csMAyARPZTTTQA1j7t4x4R+eC0j5UXEezYsQNHjx7lrYi3Hj16GK8h0m58zvh88Tmz1m3Hjh0MgEQ2E8gKYFXRfwB444033nhz5q2qEJFt5LcHcJ+cuQcwSvSbvzxvvPHGG2+Ou1UVfR0gIpvoIyJbRVu/nCsiQ0Vkh5x5CpiIiIiIbORVEdktImlScB9AIiIiIiIiIiIisrOiTAshNVRE1orIUdHnbKroiWsqmpmih5NuNV1IBEgQkfkikioih0XkF7PlWF5l0e/HPSJySEQWiUhToxVZS1sR+Vn0364cEXGd9vlrROQn0XeLdorIgLBWR0RhV9RpIaQGi8i1IhIjunl6ioisNlpR5HhERJJFX3wYAAuXIBr6Ooh+X7pEpKHRiqzvC9EAU1H0QENv0fB8nsmiLKSFaAjsJGcGwLIisktEBolIKRG5SnS/+LNhrpGIwqio00Iof/VF/zGtYLoQi6smeiDpZH9KBsDC/Swib5guIsKsEZGnT/l1GdGvNQbnvG6WMwPgo6Irp6d+7BkR2RTGuogojMpL0XsFUv5eFA3RVLh5IvLEif9nACzcuSLiE5HhIrJURA6IyHIRud9kURGgvYj8ICIXicg5IvKSiGwUXUGlXPkFwJEi8u1p90s4cb+yYaqLiMIokGkhdKbmovsmW5guxOK6iwbAkxgAC1dV9DnaLSLXib5Q3yciWXLmD2uU61IRmSv63HlFV7QSjFZkTfkFwI9EZNpp96tz4n6XhKkuIgojrgAW392ie7Ramy7E4uJF9xZVP+VjDICFO/l9OeS0jyeLHkKiM0WJyN8iMkF0O4ZL9HvziOjhBsrFFUAiEpGiTwuhXB1Ew19z04VEgEdFJFP0a2r/iZtf9Pkbb7Auq9skDICBqCj6dVX/tI+vFD3oRrnyC4CPyJl7AJ8V7gEksjVOCwlMT9EWE01MFxIh3KJvIZ1684vIv4SnMwvzjOjKaX3R1a3WIpIuIjeYLMri1oluXSkn+pzdLSIZItLMZFEW4hLdD9lSNACWPvHrKNFVvhQRGSj6PXu1iGwXngImsr1XhdNCisovuhcr9cTt2In/MhAWHdvAFM1Loi/CR0VkhWigoYLVFO0zuVf0rd/fJPfgEelqvF/0+y/nlP8/2SvxKtHT58dFf/job6BGIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIhv6f0fkYzk7AFEjAAAAAElFTkSuQmCC\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "\n",
    "%matplotlib notebook\n",
    "\n",
    "def plotCities(cities):\n",
    "    x=[city[1] for city in cities]\n",
    "    y=[city[2] for city in cities]\n",
    "    \n",
    "    plt.scatter(x,y)\n",
    "        \n",
    "def plotEdges(cities, edges):\n",
    "    x=[]\n",
    "    y=[]\n",
    "    for c in range(0,len(edges)):\n",
    "        i=edges[c]\n",
    "        j=edges[(c+1)%len(edges)]\n",
    "        x.append(cities[i][1])\n",
    "        x.append(cities[j][1])\n",
    "        y.append(cities[i][2])\n",
    "        y.append(cities[j][2])\n",
    "    plt.plot(x,y)\n",
    "\n",
    "plt.ylim((0,ysize))\n",
    "plt.xlim((0,xsize))\n",
    "plt.axes().set_aspect('equal')\n",
    "plotCities(cities)\n",
    "plotEdges(cities, edges)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "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.5.1"
  },
  "widgets": {
   "state": {},
   "version": "1.1.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}