search-4e.ipynb 258 ko
Newer Older
3001 3002 3003 3004 3005 3006 3007 3008 3009 3010 3011 3012 3013 3014 3015 3016 3017 3018 3019 3020 3021 3022 3023 3024 3025 3026 3027 3028 3029 3030 3031 3032 3033 3034 3035 3036 3037 3038 3039 3040 3041 3042 3043 3044 3045 3046 3047 3048 3049 3050 3051 3052 3053 3054 3055 3056 3057 3058 3059 3060 3061 3062 3063 3064 3065 3066 3067 3068 3069 3070 3071 3072 3073 3074 3075 3076 3077 3078 3079 3080 3081 3082 3083 3084 3085 3086 3087 3088 3089 3090 3091 3092 3093 3094 3095 3096 3097 3098 3099 3100 3101 3102 3103 3104 3105 3106 3107 3108 3109 3110 3111 3112 3113 3114 3115 3116 3117 3118 3119 3120 3121 3122 3123 3124 3125 3126 3127 3128 3129 3130 3131 3132 3133 3134 3135 3136 3137 3138 3139 3140 3141 3142 3143 3144 3145 3146 3147 3148 3149 3150 3151 3152 3153 3154 3155 3156 3157 3158 3159 3160 3161 3162 3163 3164 3165 3166 3167 3168 3169 3170 3171 3172 3173 3174 3175 3176 3177 3178 3179 3180 3181 3182 3183 3184 3185 3186 3187 3188 3189 3190 3191 3192 3193 3194 3195 3196 3197 3198 3199 3200 3201 3202 3203 3204 3205 3206 3207 3208 3209 3210 3211 3212 3213 3214 3215 3216 3217 3218 3219 3220 3221 3222 3223 3224 3225 3226 3227 3228 3229 3230 3231 3232 3233 3234 3235 3236 3237 3238 3239 3240 3241 3242 3243 3244 3245 3246 3247 3248 3249 3250 3251 3252 3253 3254 3255 3256 3257 3258 3259 3260 3261 3262 3263 3264 3265 3266 3267 3268 3269 3270 3271 3272 3273 3274 3275 3276 3277 3278 3279 3280 3281 3282 3283 3284 3285 3286 3287 3288 3289 3290 3291 3292 3293 3294 3295 3296 3297 3298 3299 3300 3301 3302 3303 3304 3305 3306 3307 3308 3309 3310 3311 3312 3313 3314 3315 3316 3317 3318 3319 3320 3321 3322 3323 3324 3325 3326 3327 3328 3329 3330 3331 3332 3333 3334 3335 3336 3337 3338 3339 3340 3341 3342 3343 3344 3345 3346 3347 3348 3349 3350 3351 3352 3353 3354 3355 3356 3357 3358 3359 3360 3361 3362 3363 3364 3365 3366 3367 3368 3369 3370 3371 3372 3373 3374 3375 3376 3377 3378 3379 3380 3381 3382 3383 3384 3385 3386 3387 3388 3389 3390 3391 3392 3393 3394 3395 3396 3397 3398 3399 3400 3401 3402 3403 3404 3405 3406 3407 3408 3409 3410 3411 3412 3413 3414 3415 3416 3417 3418 3419 3420 3421 3422 3423 3424 3425 3426 3427 3428 3429 3430 3431 3432 3433 3434 3435 3436 3437 3438 3439 3440 3441 3442 3443 3444 3445 3446 3447 3448 3449 3450 3451 3452 3453 3454 3455 3456 3457 3458 3459 3460 3461 3462 3463 3464 3465 3466 3467 3468 3469 3470 3471 3472 3473 3474 3475 3476 3477 3478 3479 3480 3481 3482 3483 3484 3485 3486 3487 3488 3489 3490 3491 3492 3493 3494 3495 3496 3497 3498 3499 3500 3501 3502 3503 3504 3505 3506 3507 3508 3509 3510 3511 3512 3513 3514 3515 3516 3517 3518 3519 3520 3521 3522 3523 3524 3525 3526 3527 3528 3529 3530 3531 3532 3533 3534 3535 3536 3537 3538 3539 3540 3541 3542 3543 3544 3545 3546 3547 3548 3549 3550 3551 3552 3553 3554 3555 3556 3557 3558 3559 3560 3561 3562 3563 3564 3565 3566 3567 3568 3569 3570 3571 3572 3573 3574 3575 3576 3577 3578 3579 3580 3581 3582 3583 3584 3585 3586 3587 3588 3589 3590 3591 3592 3593 3594 3595 3596 3597 3598 3599 3600 3601 3602 3603 3604 3605 3606 3607 3608 3609 3610 3611 3612 3613 3614 3615 3616 3617 3618 3619 3620 3621 3622 3623 3624 3625 3626 3627 3628 3629 3630 3631 3632 3633 3634 3635 3636 3637 3638 3639 3640 3641 3642 3643 3644 3645 3646 3647 3648 3649 3650 3651 3652 3653 3654 3655 3656 3657 3658 3659 3660 3661 3662 3663 3664 3665 3666 3667 3668 3669 3670 3671 3672 3673 3674 3675 3676 3677 3678 3679 3680 3681 3682 3683 3684 3685 3686 3687 3688 3689 3690 3691 3692 3693 3694 3695 3696 3697 3698 3699 3700 3701 3702 3703 3704 3705 3706 3707 3708 3709 3710 3711 3712 3713 3714 3715 3716 3717 3718 3719 3720 3721 3722 3723 3724 3725 3726 3727 3728 3729 3730 3731 3732 3733 3734 3735 3736 3737 3738 3739 3740 3741 3742 3743 3744 3745 3746 3747 3748 3749 3750 3751 3752 3753 3754 3755 3756 3757 3758 3759 3760 3761 3762 3763 3764 3765 3766 3767 3768 3769 3770 3771 3772 3773 3774 3775 3776 3777 3778 3779 3780 3781 3782 3783 3784 3785 3786 3787 3788 3789 3790 3791 3792 3793 3794 3795 3796 3797 3798
       "              '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",
       "            if (mpl.ratio != 1) {\n",
       "                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n",
       "            }\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",
       "        fig.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 backingStore = this.context.backingStorePixelRatio ||\n",
       "\tthis.context.webkitBackingStorePixelRatio ||\n",
       "\tthis.context.mozBackingStorePixelRatio ||\n",
       "\tthis.context.msBackingStorePixelRatio ||\n",
       "\tthis.context.oBackingStorePixelRatio ||\n",
       "\tthis.context.backingStorePixelRatio || 1;\n",
       "\n",
       "    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\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 * mpl.ratio);\n",
       "        canvas.attr('height', height * mpl.ratio);\n",
       "        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\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'] / mpl.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n",
       "    var x1 = msg['x1'] / mpl.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\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 * mpl.ratio;\n",
       "    var y = canvas_pos.y * mpl.ratio;\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",
       "    var width = fig.canvas.width/mpl.ratio\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 + '\" width=\"' + width + '\">');\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 width = this.canvas.width/mpl.ratio\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\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,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4nOzdeVhU9f4H8A+bZKCIO2rivo0wyKJmppWllVct8VZa3UzS8tpyK5Os27XbTbtmVrds0UrLJU3LMi1/loVLpeAOuIEm7qKibCrrvH9/fGEYEBA4cM4s79fzzKPMOXPOZ5iZM2/O+S4iRERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERM7MW0TmiEiyiOwVkcVF93cWkT9EJElE4kSkh81jKltGRERERHbuHRF5T0Tcin4OKPr3VxEZW/T/USKyxeYxlS0jIiIiIjvmIyLpIuJb5v7mRfd7Fv3sJiJnRKTdNZYRERERkZ0LFpHDIvJfEdkuIptFZJCIhInIvjLrxonIgGssIyIiIiI7FyYiEJG/Ff1sFpFzItJHVHtAW9ukJABWtKys50TkRPHNw8MDrVu35o033njjjTfeHOgmIrm1kjrIbjQVkUIR8bC5L05E7hORDKn4EnBFyyrVunVrEBERkWMRdSKHnMxPInJ30f8DRZ0BDBCRDVK6o8dWm8dUtqxCDIBERESORxgAnVIHUYEuQUR2i8i9Rfd3FdW7N0lU+0CTzWMqW1YhBkAiIiLHIwyApAUDIBERkeMRBkDSggGQiIjI8QgDIGnBAEhEROR4hAGQtGAAJCIicjzCAEhaMAASERE5HmEAJC0YAImIiByPMACSFgyAREREjkcYAEkLBkAiIiLHIwyApAUDIBERkeMRBkDSggGQiIjI8QgDIGnBAEhEROR4hAGQtGAAJCIicjzCAEhaOHsAzM3NxZQpU9CxY0d069YNJpMJ8+fPr5VtBwYGIiEhoVa2RUREVB3CAEhaOHsAHD16NEaOHIns7GwAwJEjR9CtWzfMnTv3qnXz8/OrtW0GQCIiMoowAJIWzhwAk5OTUb9+fZw/f77U/T/88APatGmDmJgYmM1mPPXUU+jbty+WL1+OJUuWoHfv3ggJCYHZbMYPP/xgfdymTZvQs2dPREREYNKkSWjbtq01ACYlJeHuu+9GeHg4goOD8cEHH1gf9+CDDyIsLAxBQUEYOnQoUlNT9fkFEBGR0xIGQNLCmQPgV199heDg4Kvuv3DhAkQEy5cvh5ubGzZv3mxddv78eVgsFgDqbGFAQADy8vKQk5ODVq1aISYmxrptEUFCQgIKCgoQHh6O/fv3AwAuXbqEoKAg7NixAwBw7tw56/bfeOMNTJo0qa6eMhERuQhhACQtnD0Ams3mq+63DYBdu3YttWzbtm0YMmQITCYTzGYzvLy8kJycjD179qBTp06l1vXz80NCQgL27t2L+vXrw2w2W2/t2rXDokWLAADvvvsuwsLC0LNnT3Ts2BE33XRT3T1pIiJyCcIASFo4cwBMSkq65iXgsLCwUss6duyIb7/91vqzv78/EhISsHv37goDYGJiIm644YZya9i8eTM6deqEs2fPAgBWrVp11T6JiIiqSxgASQtnDoAAcN9992HkyJG4dOkSgJJOIB999FG5AbBx48bYtWsXAGDRokXWy7w5OTlo3bo1Nm7cCABYsWKFdVl+fj66du2KL774wrqd5ORkpKWl4fvvv0doaCgKCgqQm5uLu+++mwGQiIg0EwZA0sLZA2BOTg4mT56MDh06WIeB+eSTTwCg3AC4aNEitGvXDjfddBMmT55cqqOHbSeQF1544apOIEOHDkVQUBB69OiBgQMH4sSJE8jPz8d9992HTp064ZZbbsFLL73EAEhERJoJAyBp4ewBsEIWCxAbC0RHA+PGqX9jY42uioiIqEqEAZC0cMkAmJAAmEyAlxfg6QmIqH+9vNT9iYlGV0hERFQpYQAkLVwuACYkAD4+gLu7Cn5lb+7uajlDIBFRlWQeP250CS5JGABJC5cKgBaLOsNXUfizDYEmk9HVEpELs1gsSFy4EPHz52P/0qXYt2QJUtavr5Vt71+6VPM2CvPzcWD5cuxdvBibpk69avm5vXsRN2sW9i9dar2dLWozfensWeyeOxf7ly7F3sWLsWfePFgKCwEABbm5iJ8/H/uWLLFuv7CCWZr2L1uGvYsXY//Spdj98ce4eOhQuesd/OabWnnO9kYYAEkLlwqAsbHqMm9l4a/45uUFxMUZXTERuaiEBQuQXWbWoDM7d+LkH39o3nZtBcD8K1cq3F7Kzz/jks0g+Lbi58+3DrgPAJfOncPBr78GACR+8QXyLl+2Lsu7fBlJNkNzFTu2YcNVgW/3xx9ftd75ffuQsn49AyBRWS4VAKOjS9r8Xevm6anWJyLSmcViQaLNsFK2DttMT3nwm28Q/9ln2Pfll9j14YdIO3jQuixl/XrsW7IE+5ctQ/z8+Ti3dy9y0tOxf+lS/PL009i/dCkKcnNrpd7ywtX+ZctKhTxb+5Ysueq+vYsXq2VffnnVsoPffFOtbRQryM3Fro8+qrBGRycMgKSFSwXAceOqFv6Kb1FRRldMRC4o+8wZHPnpp0rXORoTg1Nbt5a6b/s776AgLw/A1WHINjBVFIYK8/NLXbItvl08fLjSWsrb3t5Fi3D4hx+w/6uvsO/LL3G+aKpMQIW87NOnrT+fio3Fid9+u6rOYvGffXbVfReSknBo9Wrrz7lZWdgzb16pdfZ88glys7MrrNHRCQMgaeFSAZBnAInIAaQfOYKTW7ZUuk55oSjz+HFrO8Ed//sfcjMzy31sbYeh8rZ3Nj4eeUUD8ANA0sqVuJCcDAAoLCjAr889h7i33sLmf/4Tuz780Lpeyvr11raCAJC6axfWP/lkufs9sGIFYp5/Hjveew8/TZxoDXsAcOL333HKphkPAyBRGS4VANkGkIgcQPbp0zgaE1PpOuWdKbNYLNZLqIUFBTiybh2SVq7Ewa+/LnW5V48AWJ7i2vbMm1cqrF38889SZ/OOxsRYO76c3rYNCQsWXLWto7/8grN79lh/LsjNRfynn6KwoAA56elXPYYBkKgMlwqA7AVMRA7AUlh41SXcYofWrAGgLrGWlZ6SgqO//HLV/XmXLmFP0QxIgD6XgM/t3Xv1el99hctpafhz7dqrlpUXaAEVaovb8V1r/fP79uHMzp2q9/CiRaWewy9PP439X31V6fNwNMIASFq4VAAEgMREFNb3QYGUHwILxB0WjgNIRAZLXLjwql60Z+PjceL33wEAxzZuRGrRvOXFtr/zDgrz85Fx9CjOxseXWpbw+efW/+txBjD+s8+sQ7sAQNapU0j5+WdYLBbEz59fal2LxYL4Tz8td9u7P/643HEGExcuvKoTy4EVK3DlwoUq1+johAGQtHC5AAjg5XsSkSAmFHqWngmkwMMLCWLCW48y/BGRsSwWC/YtWaLGuVu2DPuWLLmqY8ih779H4sKF2LtoEXbOmYMLSUkly1avxv6vvsL+pUsR/9lnOLNjh3XZ3sWLVUeMM2dqpdbywlVOeroaw3DZMjVen80Zy9Tdu5Hw+efWZXs++QSXz5+3Lj+9fbv1/rQDB6z3Xzx82BqA8y5fRvynn6ozfMuWIeHzz63LqlqjoxMGQNLC1QJgUhLg4QHccgtUG7/oaNXbNzoaBVviEBqqmv/ZHHOIiIjsjjAAkhauFgAfflid8Nu4sfzlv/+ult95p2oySEREZI+EAZC0cKUAeOCA6t9x++2Vr1ccEr//Xp+6iIiIqksYAEkLVwqAY8aoYFdJMxEAwKlTgK8v0KEDUDTTERERkV0RBkDSwlUC4N69gJuburRbFW++qcLi9Ol1WxcREVFNCAMgaeEqAfC++1Sgi42t2vq5uUCXLsD11wPHjtVtbURERNUlDICkhSsEwPh4Ff6GDq3e49auVY+7//66qYuIiKimhAGQtHCFADhypApy27dX/7HDh6vHbthQ+3URERHVlDAAkhbOHgB37VIBbsSImj3+8GHA2xsICgLy82u3NiIiopoSBkDSwtkD4IgRKgCWmTGpWl5+WW1jzpzaq4uIiEgLYQAkLZw5AG7froJbZKS27WRnA23aAP7+QJmpOYmIiAwhDICkhTMHwKFD1dAvCQnat7VsmQqTjz+ufVtERERaCQOgU0oRkQMisrvodv817hcR6Swif4hIkojEiUiPquzIWQPg1q2124PXYgEGDlSB0mZOdSIiIkMIA6BTShGRntW4X0TkVxEZW/T/USKypSo7ctYAOGSICmt799beNvfsUVPJ9evHeYKJiMhYwgDolFKkegGwuYiki4hn0c9uInJGRNpda0fOGAB//12d/Rszpva3/eSTatuLFtX+tomIiKpKGACdUoqI7BGRBBH5VESaXeP+MBHZV2YbcSIy4Fo7csYAePvt6kzdgQO1v+20NKBJEyAgAMjMrP3tExERVYUwADqltkX/eonITBH58Rr3h4nI3jLb2CblB8DnRL1pTojICT8/P6Pfw7Vq40Z1hu5vf6u7fcydq/YxZUrd7YOIqi8wMBBdu3ZFvs2gnWFhYYiJianzfR85cgRz586t8/0QFRMGQKcXICJZ17i/uYhkCC8B45ZbAA8PIDm57vZRUACEhgJeXnVzlpGIaiYwMBDt2rXDxx9/bL1PrwAYExODsLCwOt8PUTFhAHQ6PiLSyObn50RkUyX3F9sgpTuBbK3KzpwpAP76qzozN25c3e+ruJ3hnXeyQwiRvQgMDMQXX3yB1q1b49KlSwBKAmBmZiYee+wxREREICgoCI8//jjy8vIAALNnz0Z4eDhCQkIQERGBrVu3WrcpIpgxYwYiIiLQvn17/Pzzz3jxxRcREhKCHj16IDExEQDQtWtX1K9fH2azGcOGDdP/yZPLEQZAp9NBRHaJSLyotn6rRJ3Jq+j+Yl1F9fxNEpHtImKqys6cJQBaLED//oCnJ/Dnn/rs8+GHVQj8/nt99kdElQsMDERCQgJGjx6N119/HUBJABw/fjwWLlwIALBYLIiKisLbb78NADh79qx1G1u2bIHJZLL+LCKYUzQN0PLly3H99ddjzZo1AICZM2di9OjRAHgGkPQnDICkhbMEwJ9+UmFs/Hj99nnqFODrC3ToAFy5ot9+iah8xQHw8OHDaNq0Kc6fP28NgM2aNUNQUBDMZjPMZjO6dOmCiRMnAgDWrVuHAQMGwGQywWw2w83NDbm5uQBUADxXNAXQoUOH4Ovra93f+vXr0adPHwAMgKQ/YQAkLZwhAFosQN++qk1eSoq++37zTRU8p0/Xd79EdLXiAAgATz31FP7xj39YA2DTpk1x+PDhqx6Tm5sLX19fbN++HQCQkZEBEUFWVhYAlPr/kSNH0KRJE+tjbUMfAyDpTRgASQt7CICrVq2y/lVefGvVqhW8vb2r9Pi1a1UImzgRuOuuu3Do0KFq1/DII4/g/fffr/bjcnOBLl2A668Hjh2r9sOJqBbZBsCzZ8+iadOmCAgIQExMDMaNG4fx48dbewhfuHABycnJyMjIgLe3N1JTUwEAr7/+eo0C4I4dO9CpUyfdniuRMACSFvYQAMu6ePEiOnXqhE8++eSa61osQEQE4OVViKNHC2u8z5oGQKAkgD7wQI13T0S1wDYAAsBrr70GEVGdQDIyMDEyEqamTRHk74/QFi3wc9FnfubMmQgMDMTNN9+MWbNm1SgA5ufnY+jQoTCZTOwEQroQBkDSwt4CYGFhIe6++248/vjj1vtmzZqFiIgI9OrVC3fddReOFZ1qmzZtGm655SGI3IsmTYJw4sSJUl8AAwcOxJQpU9C/f3906NCh1DZPnDiB2267DUFBQRg+fDiGDh1a4wAIAMOHqxC4YUONN0FEdSUhATCZVDsRT0/1YfX0VD+bTEBRT14iRyIMgKSFvQXAf/7zn+jbt6+1AfaSJUswfvx4FBQUAAAWLlyI4cOHAwD+9a9p8PJqDW/vVJw8qR5fNgBGRkaioKAAly9fRrt27fDHH38AAEaOHIlXX30VAHD48GH4+vpqCoCHDgHe3kBwMGAzBi0RGS0hAfDxUdMDiVx9c3dXyxkCycEIAyBpYU8B8LvvvkNAQABOFqc5AH/961/Rvn17a9vAnj17wmw2AwDuv38aRB7Hs8+WbKNsAFyxYoV12T333INFRZP4+vv748SJE9ZlI0aM0BQAAeDll9X3SdGIEURkNItFneGrKPzZhkCboV+IHIEwAJIW9hIADxw4gCZNmmDz5s2l7h81ahQ+++yzq9YvLASaN58GT8/nceZMyf1lA+Dq1autyyIjI7FgwQIAdRMAs7OBNm0Af3+gaNQIIjJSbKy6zFtZ+Cu+eXkBcXFGV0xUZcIASFrYQwDMzMxE9+7drYOt2lq4cCHCwsKQlpYGAMjLy8POnTvx9deAyDSEhz9fav2qBsDIyEi89tprAFTDbq2XgIstW6a+S2yaGxKRUaKjS9r8Xevm6anWJ3IQwgBIWthDAHzjjTfg7u5+1VAwZrMZJ0+exDtvv42gTp0Q3KwZejRqhLduvRX3t4+Fl9c0TJxYswBo2wkkMjIS99xzT60EQIsFGDgQcHMDduzQvDki0mLcuKqFv+JbVJTRFRNVmTAAkhb2EAArVU7vvUJ3T+SKF840tc/ee3v2qCZF/fpxnmAiQ1XjDKCFZwDJwQgDIGlh1wHwGr33LHbce+/JJ1WZRX1OiMgIsbGwVLENYL6bF3Z8zDaA5DiEAZC0sNsA6OC999LSgCZNgIAAIDPT6GqIXNOVyxak+JpQIJUfRwrFHYligggwZAhQNCsckV0TBkDSwm4DoBP03ps7V5U3ZYrRlRC5nitXgDvvBHpIInI8fdQVg0rGATz6YyIefFC13xUBRo0C9u0z+lkQVUwYAEkLuw2ATtB7r6AACA1V+fTAAaOrIXIdly8Dd9xR0q+jMD6xyjOBxMcDI0aUZMOxY4GUFAOfDFEFhAGQtLDbAOgkvfd++02Vd+ed7BBCpIdLl4BBg9TnbsIENWaoVVyc+mMxKkr9W8mVgy1bgFtvLbnI8NRTKDXmKJHRhAGQtLDbAOgEZwCLPfSQKvP7742uhMi5ZWeXhLaJE8uEvxpavx7o3Vtt8/rrgZdeAi5c0L5dIq2EAZC0sNsA6ARtAIudOgX4+gIdOqh2SURU+7KygAED1CHhySdr94y7xQJ8+626WiwCNGoEzJihAieRUYQBkLSw2wDo4L2Ay3rzTVXu9OlGV0LkfDIzgf791WfsmWfqrrlFQQGwcCHQvr3aV4sWwHvvATk5dbM/osoIAyBpYbcBEFANs318UOhWee89exwHsKzcXKBLF3UJ6fhxo6shch4ZGWrQdRHguef0aWubmwt8+KEa5kkECAwEFiwA8vPrft9ExYQBkLSw6wAIAImJ+PN6E3LFS43UX0nvPXu3dq0q/4EHjK6EyDmkpwN9+6rP1Qsv6N/R6tIldXa/cWNVQ7duwIoV7PBF+hAGQNLC3gNgUpI6sL48uOq99+zZsGHq+WzYYHQlRI7t4sWSzhkvvmhs6EpPB155RV2QEAHCwoD/+z8GQapbwgBIWth7AJw6VR1Qf/7Z6Epqx6FDQL16QHAwLxcR1dSFC0B4uDo2/POf9hO0UlOBZ58FvL1VbQMGqKGgiOqCMACSFvYcAPPzgVatVPua2hjOwV689JL6cpgzx+hKiBxPWpoaYF0EmDbNfsKfrWPHgMceAzw8VJ1DhwK7dhldFTkbYQAkLew5AK5Zow6er75qdCW1KzsbaNMG8PcHzp0zuhoix3H+PBASoo4L//630dVc28GDwP33l/Rbu/9+dR9RbRAGQNLCngPgyJFqXk5nnIZp2TL1hfD440ZXQuQYzp5VTSdEgNdfN7qa6tm1S50FFFFnBR97TJ0lJNJCGABJC3sNgKmpqrPvHXcYXUndsFiAgQNVwN250+hqiOxbairQs6cKUG+8YXQ1NffbbyWDVderB/zjH+q5EdWEMACSFvYaAGfPVgfJpUuNrqTu7NmjhjK86Sb7bMdEZA/OnAF69FDHg1mzjK5GO4tF9RAubsfo66t6EKenG10ZORphACQt7DEAWizqgO/v7/xTp02apL4EFi0yuhIi+3PqlBpbTwR4+22jq6ldFosaM7D4+fn7AzNnqrEFqyMwMBAJCQnV3v+2bdswZsyYaj+O7IcwAJIW9hgAt24tmc/T2aWlAU2aqBkFMjONrobIfpw8qWbPEQHefdfoaupOfr6aRaRtW/VcAwLULCO5uVV7fE0DIDk+YQAkLewxAE6YoA6ErjJswscfq+c7ZYrRlRDZh+PHgc6dXWu4pJwcNa9w8+bqebdvr+YdLiio/HHFAXD27NkIDw9HSEgIIiIisHXrVgDA5cuXcd9996F79+4IDg7GHUUNq2NiYhAWFlbXT4vqkDAAkhb2FgCzs4EGDYBevYyuRD8FBao9kJcXh4ggOnYM6NhRhaCPPjK6Gv1lZwMzZgCNGqnfgckEfPttxe2EiwPg2bNnrfdt2bIFJpMJALBy5Upr6AOAtLQ0AAyAzkAYAEkLewuAX3zhWn/1F/vtN/W877qLHULIdaWkqDNfIsDcuUZXY6wLF9Sg8ddfr34fERHlz4hUHADXrVuHAQMGwGQywWw2w83NDbm5uTh8+DBuuOEGTJw4EcuWLUNmUVsTBkDHJwyApIW9BcABA9Q0ShcuGF2J/h56SB3oV682uhIi/R05ArRrp4ZG+vRTo6uxH6dPA089pa4QiAC33QZs2VKyPDAwEHv27IGvry+2b98OAMjIyICIICsrCwCQlZWFVatW4amnnkLbtm1x4cIFBkAnIAyApIU9BcCkJHWAGz3a6EqMcfKkGhKiY0fn7/1MZOvwYdUJws1NdYigqx05Aowdq4aOEgGGDwfi40sCoLe3N1KLBhV8/fXXrQHw+PHjyM7OBgDk5ubihhtuwJ49exgAnYAwAJIW9hQAp05VB7byLnO4ipkz1e9g+nSjKyHSR3KymhrR3V11eqDK7dsHjBqljhNuboCPTyB+/DEBM2fORGBgIG6++WbMmjVLBcDMTPz4zjsIbtYMQf7+MDVtipfGjgXAS8DOQBgASQt7CYD5+UCrVkBgIFBYaHQ1xsnNVUNfXH+96glJ5MySkoDWrVX4W7LE6Gocy/btwJAhKgh6eqppJU+csFkhIUH1IPHyUisUr+jlpe5PTDSsdqodwgBIWthLAFyzRh2fXn3V6EqM9+OP6nfxwANGV0JUdw4cUGPeeXg494w/dW3DBqBfP3XMuO46YPJk4MKmBMDHp+R6cdmbu7tazhDo0IQBkLSwlwA4cqS6nJGSYnQl9mHYMHWc3rjR6EqIat++fUDLlir8LV9udDWOz2JRf0SbzYCIBXvdTSiUCsKfbQgsGiqGHJMwAJIW9hAAU1PVlQmboapc3qFDarL44GB1eZzIWSQmAi1aqM/8118bXY1zKSwEfno9FnniVXn4K755eQFxcUaXTTUkDICkhT0EwNmz1bGIl4FKe+kl9Xv54AOjKyGqHfHxQLNmKnd8+63R1Tip6GhYitv8Xevm6QlERxtdMdWQMACSFkYHQIsF6NFDTYTOoU9Ky85WDeT9/YFz54yuhkib3buBpk1V+Pv+e6OrcWLjxlUt/BXfoqKMrphqSBgASYvaDICZmZnw8fFBVDUOKFu3qmPQk0+Wv1xsBjN1RUuXqt/P448bXQlRze3aBTRurJo1rFljdDVOLjq6pNcvzwA6NWEAJC1qMwDOmzcPN998Mxo1alRhaCsoM7P5hAnqOLRrV4VvcJcOgBaLmh3FzQ3YudPoaoiqb8cOdRbb21v1cKc6FhtbMm0I2wA6NWEAdEopInJARHYX3e4vur+ziPwhIkkiEiciPWweU9myCtVmAOzTpw/WrFmDYcOG4bPPPgMALFiwAEOGDMHDDz+MsLAw/PHHH5g9ezbCw8MRHBwCd/cIdOmy1bqNb775Bl27dkXfvn3x2muvuXwABNSlM3d34KabOE8wOZa4OKBRIzU8yf/9n9HVuAiLRfXurWgIGPYCdhrCAOiUUkSkZzn3/yoiY4v+P0pEtlRxWYVqKwAmJiaiVatWKCgowHfffYd+/foBUAHQx8cHSUlJ1nXPnj0LAPjiC0BkCwIC1EEoNTUVjRs3xoEDBwAAM2fOZAAsMmmSOmYvXmx0JURVs3Ur4OcH1K/v2rP7GCIxsdJxAAvdOA6gMxAGQKeUIlcHwOYiki4inkU/u4nIGRFpd41llaqtAPjss8/ixRdfBADk5+ejZcuW2Ldvn/UMoK1169ZhwIABuP56E9zczHBzc0Nubi5WrVqF22+/3brexYsXGQCLpKUBTZqogXMzM42uhqhyf/wBNGyowt8vvxhdjYtKTLxqJhCLpyfyxAsH65mQv5vhz9EJA6BTShGRPSKSICKfikgzEQkTkX1l1osTkQHXWFap2giAeXl5aNasGVq1aoXAwEAEBgaiYcOGmDx5MhYsWIDIyEjrurm5ufD19cXKldshAkRGZlhD3nfffccAWImPP1Z/wLPNNtmz338HGjRQ0xlu2GB0NYS4OHXQiIoCoqMx97E4iAALFhhdGGklDIBOqW3Rv14iMlNEfhQV8vaWWW+blATAipaV9ZyoN80JETnh5+en+U24YsUK9OnTp9R9CQkJaN68OebNm1cqAGZkZMDb2xtPP50KEWDs2NetIS81NRVNmjTBwYMHAaBkQnMGQABAQQHQq5f6g77oV0RkVzZtAnx91dXFTZuMrobKk/nDS3oAACAASURBVJ6uOuW0bw/k5RldDWkhDIBOL0BEskRd5s2Qii8BV7SsUrVxBvDOO+/Ee++9d9X9ISEhmD17dqkACABvvDETHh6B8Pa+GW++WTrkffPNN+jSpQtuvPFGzJ49mwGwjN9+U2cB77qLHULIvmzYoIJfgwbqfUr2a/p0dRyZO9foSkgLYQB0Oj4i0sjm5+dEZFPR/zdI6Y4eW23Wq2xZhXQZCNpiUUMTREcD48YheVQ0IiQWr75a97t2Rg89pA7eq1cbXQmR8uuv6pJvgwaq/R/Zt6wsNSh3mzYcgN+RCQOg0+kgIrtEJF5UG8BVUnImr6uo3r1JIrJdREw2j6tsWYXqPAAmJFzVEDnfzRO54oXczib2QquBkyfVZbaOHXnwJuOtX686ezRsqHr+kmOYNUv9IVnOxRtyEMIASFrUaQBMSKh0KAK4cyiCmpo5U/0KZ8wwuhJyZevWqTH+GjXieMKO5tIloGVLdbt0yehqqCaEAZC0qLMAyMFI61RuLtCli7rsdvy40dWQK1q7Vs3u4e8PbN9udDVUE//7nzoMv/WW0ZVQTQgDIGlRZwGQ0xHVuR9/VL++0aONroRczZo1al7fxo05RaEju3JFtQNs2lS1CyTHIgyApEWdBUBOSK6LYcPUr3DjRqMrIVfx/ffqb7YmTdQ0heTYiscXZXMSxyMMgKRFnQXAceOqFv6Kb1FRdVOHk0tOVmdigoOB/HyjqyFn9913Kvw1awbExxtdDdWG3FygXTt1KT893ehqqDqEAZC04BlAx/fSS+rX+MEHRldCzuybb9RHtXlz9ttyNvPnq2MIh+ZyLMIASFqwDaDjy8oCWrdWf8GfO2d0NeSMli8HPDxUj9F9+4yuhmpbfj7QubMayictzehqqKqEAZC0YC9g57B0qfp1PvGE0ZWQs1m2TIW/gADgwAGjq6G6snixOoZMnWp0JVRVwgBIWtTpOICJiRwHUCcWCzBgAODmxl6ZVHuWLFEf01atOP+0sysoAHr0UIfks2eNroaqQhgASYs6nwkkMdE6E0ihh2oTWODuqS77mjgTSG3avVt9Wd90E+cJJu0WLlTvpzZtVGcjcn7Ll6u/zZ9/3uhKqCqEAZC00GUuYACIi8P58dGYJ1HYeGM02/zVkUmT1AF88WKjKyFHtmCBOpvcti1w+LDR1ZBeCgsBs1nN7nLqlNHV0LUIAyBpoVsABJCSwr8u61pamhqfLSAAyMw0uhpyRJ99psJfYCBw5IjR1ZDeVq1Sx+mnnjK6EroWYQAkLfQMgGfOqAPLpEm67dIlFQ/sypF1qLrmzVPvnfbt1R9s5HosFiAiQo0veuyY0dVQZYQBkLTQMwCmp6svl3HjdNulSyooAHr1Us0s2XCfquqjj9Tns0MH4OhRo6shI61dq94LEyYYXQlVRhgASQs9A2BurjqojBmj2y5d1m+/qd/13XcbXQk5gjlz1PulUyfg+HGjqyGjWSxAv35q4G+2AbVfwgBIWugZAC0W1bZo5EjddunSHnpIfamvXm10JWTP/vc/9T7p0gU4edLoashe/PKLel+MHWt0JVQRYQAkLfQMgABQvz7PSunl5EnA1xfo2BHIyTG6GrJHb7+tvuS7dWOvT7raLbeooYDYlMQ+CQMgaaF3APT3B269VdddurSZM9UX/IwZRldC9mbWLPXe6NEDOH3a6GrIHm3ezGY79kwYAEkLvQNgq1bAjTfqukuXlpOj5vi8/nq27aIS//2v+mI3mYDUVKOrIXs2ZIhqusMx++2PMACSFnoHwA4dVA9V0s+PP6ov+9Gjja6E7MHrr6v3Q1AQp/yia4uNVe+XUaOMroTKEgZA0kLvANijB9C9u667JAB/+Ys6iG/caHQlZKR//1u9D0JCgHPnjK6GHMWwYep9s2uX0ZWQLWEAJC30DoBhYUC7drrukqDmcq1XDwgOBvLzja6G9GaxANOmqS/xXr3UjDFEVbVrl3rvDBtmdCVkSxgASQu9A+BNNwEtW+q6Syoydao6iH/wgdGVkJ4sFuCf/1SvfVgYcOGC0RWRI4qMVO+h2FijK6FiwgBIWugdAAcNAho10nWXVCQrC2jdWvXEPn/e6GpIDxZLSfDv3Ru4eNHoishRJSaqziBDhhhdCRUTBkDSQu8AOHQo4O2t6y7JxtKlKgw88YTRlVBds1iAF15Qr3ffvmoqRiItxoxR76fffjO6EgIYAEkjvQPgqFHqAGKx6LpbKmKxAAMGqL/kd+40uhqqKxYL8Nxz6rPWrx+QkWF0ReQMDh5UA0NzLFf7IAyApIXeAbB4erIrV3TdLdnYvVsdxG+6iUHcGVkswDPPqM9Z//5AZqbRFZEzGTtWvbd+/dXoSkgYAEkLvQPgY4+pgwfbIhlryJBvIBKKtm3N6NatG2677TYUFhbirrvuwqFDhwAAAwcOxOoKJhKOiorCpk2b9CyZqsBiAZ58Un3GBgxQ7T6JatOffwKenvwD0h4IAyBpoXcALP5y4ryjxjl9+jSaNm0GP78UBASoM0Q7duyApczRvLIASPansBCYOFF9vm69FcjONroiclYTJqj32dq1Rlfi2oQBkLTQOwAWN0r/809dd0s2du7ciYCAALzzThZEgOjokmWBgYFISEgAoALgM888g4EDB6JTp06YPHmyNSTahsNHHnkE77//vnUbzz//PKZNm6bb8yEV/oq/lAcNAi5dMroicmbHjqlxRcPDeRbQSMIASFroHQBfeUV9Se3bp+tuyUZhYSFGjhwJf39/+PndA3f3N7Fp0wkAVwfAO+64A3l5ebh06RLCwsLw1VdfWZcxANqHwsKSphV33AFcvmx0ReQKiq/mrFpldCWuSxgASQu9A+D06eqgwR6oxtu/fz8mT/4YIiPg6emH5OTkqwLg4sWLreu/8847GD9+vHUZA6DxCgqARx9Vn6khQxj+SD+nTgHXXQeYzeqPENKfMACSFnoHwNmz1ZfV77/ruluqxIMPAiJDEBU1+5oBcMKECdZlxQEwKioK7777rnW9SZMmMQDqoKAA+Nvf1Ofp7rvZs5709/zz6v23YoXRlbgmYQAkLfQOgB9+qA4Yv/yi627JxokTJ/CbzUiue/degJtbF7Rs+R3ati0dAIcMGYL8/HxcvnwZERERWL58uXVZcQCcMWMGRo8eDQA4f/482rdvzwBYx/Lzi4O7mp81J8foisgVnT0L+PgAPXqoP0hIX8IASFroHQDnz1dfWj/8oOtuyUZKSgoGDx6Mzp07w2w2w2QyYciQ6RABGjUKRMKXXwLR0RjYsiWm9u2LgaGhlXYCSUtLw6233oqePXtixIgReOCBBxgA61B+PvDAA+pzNGIEkJtrdEXkyoqnGlyyxOhKXI8wAJIWegfAL79UB4tvvtF1t3QNOTnA3W0TsNfNBIunlxroS0T96+UFmExqMlCqssDAQHTt2hVmsxlmsxmPP/645m3m5QH33ademnvvZfgj46WlAQ0bAp07qz9OSD/CAEha6B0Av/1WfXnZNC0je5CQgHxvHxSIu3qByt7c3dW1HobAKrNtT1kb8vKAyEj1cowapX4uK5/fwGSAadPU+3L+fKMrcS3CAEha6B0A165VB4pPP9V1t1QZi0Wd4XOvIPzZhkCTyehqHUZ5AXDhwoWIiIhASEgIQkJCsNZmJN0NGzbAZDIhKCgIzzzzDFq1aoX9+/cDAJ566h9o1CgcImY0b34LDhxQs7UkJyejRYsW+Ne//oWbbroJ7733HtatW4e+ffuiV69eMJlMWLhwoXUfr776aqmzksePH9fhN0HOLj0d8PcH2rXjWWk9CQMgaaF3AIyJUVlizhxdd0uViY1Vl3krC3/FNy8vIC7O6IodQtlLwCtXrsS5c+es7SgPHz6MgIAAFBQU4MqVKwgICMDvRd3jly9fDhHB/v37kZMDDBlyDiLAmDHA558vwogRIwCoACgi1vEZAeDChQsoKGqRf+7cObRt2xZnzpzB2bNn4e/vjytF3YUvXbqEHPYeoVpSPMTXxx8bXYnrEAZA0kLvALhlizpIvPWWrrulykRHl7T5u9bN07P01CFUofLOAMbGxmLw4MEwmUwwm83w8vLCkSNHsGPHDnTr1q3Uur6+vtizZz+GDgVEFqFp0z4wmUzo3r07AgMDAagA6OPjU+pxBw8exD333GPdh4+PD2JiYpCfn49evXph1KhRmDt3Lk6cOFGnz59cS1YW0LQp0KYNhyTSizAAkhZ6B8Ddu1WO+M9/dN0tVWbcuKqFv+JbVJTRFTuE8gJgYGBgqfmVGzRogP3792P79u3o3r17qXV9fX1x8837IXIY113XDIcOHQGg5m1u0aIFgJJLwLYGDhyIOXPmWM80BgUFWS815+fnY+PGjfj3v/+N1q1bW884EtXUnXfeaR0I/q231CHivfeA4OBgrFy5strbq2rb2VdeeQXLli2r9vadiTAAkhZ6B8ADB9QB4uWXdd0tVYZnAOtEeV9kfn5+1vsWLFhgvcx75coVtGzZElu2bAEAfPnl1xARiOzHiBE70aZNG1y5cgUWiwXjxo2rNAAGBwdbQ+avv/4KNzc3rF27FhkZGTh16pR1vcjISHzwwQd19vzJNaxYsQKhoaEA1BzULVsCTZpsQ7NmzZFXXk+la6jtzlPOTBgASQu9A2BKisoRzz+v626pMmwDWCfK+yJbsGABAgMD0b9/f0yZMgWtW7e2dvT45ZdfYDKZEB7eG+3aTYVIEzz44CkUFqrZVdq3b48BAwbgtddeqzQArl27Fp06dULfvn0RFRWF3r17Y+3atUhJSUHv3r3Rs2dP9OzZE3/961+RmZmpzy+DnFZubi6aNWuG3bt3A1Bn/0QmYuDAyYiPj0f//v3Rq1cvdO/eHTNmzLA+7pNPPkH37t1hNpvRs2dPbN26FcDV85HbnjGPjIzEggULAFw9BaUrEgZA0kLvAHjmjMoRf/+7rrulyrAXsP4sFhW8o6PVJfjoaGT++iuys4HbbgNEfkKDBu04xyo5hGeffRbPPPMMAODixStwc/NHo0b7cPJkprWj0eXLlxESEoJt27YBABo2bIiTJ08CAPLy8pCVlQWAAbA6hAHQqU0TEYhIz6KfU0TkgIjsLrrdb7NuZxH5Q0SSRCRORHpUZQd6B8CMDJUlxo3Tdbd0LYmJapy/CkKgheMA1p6EBBWkvUoPuP2Jhwe6uXnDWzqjWbNw/P77H0ZXSlQliYmJaNq0KXJzc7FkyRJ06NAPIsDUqal46KGH0LNnTwQHB8Pf3x+ffPIJAOCee+7BoEGD8O677+LgwYPWbTEAVp0wADqtUBFZKyJHpXQA7FnB+r+KyNii/48SkS1V2YneATA3F9bhLMjOJCZeFUwK3D2RK15Ia8WZQGpFQkKlQbtA3JHj6QNLAn/X5Fj69OmDr7/+GoMGDcLcuZ+hXTugXr0oPPnk89YByu+9915raLNYLNi2bRtmzZqFjh07YunSpQBKB8BBgwbhu+++s+5j6NChDIA2hAHQKXmLCnDtpXTos/2/reYiki4inkU/u4nIGRFpd60d6R0ALRbAzQ0YOVLX3VJ1xMWpS5NRUch9Nhq3+8WhXbvyZ56gaqjipXYLL7WTA5o3bx5CQ0PRsGFDZGVlYcECQGQkhgx5BwBw4MAB+Pr64v3330d+fj6Sk5Otj42Ojsazzz4LoHQAnDBhAqZOnQoA+PPPP+Hn58cAaEMYAJ3STBGZVPT/FCkdAPeISIKIfCoizYruDxORfWW2ESciA661I70DIABcfz1w112675Zq6PXXVTYpOu5STbGzDTmxzMxM+Pj4YFxR+578fKBt251wdzchJCQco0aNwvDhw/H+++8jJycH/fv3t45Vefvtt1tnpbENgIcPH0Z4eDhCQ0MxevRoDB48mAHQhjAAOp0bRV3OdSv6OUVKAmDbon+9RIXEH4t+DhORvWW2s03KD4DPiXrTnBCRE35+frq/aRs3Bm69VffdUg2lpwONGqnJ3osmmKCa4HA75GKWLAFELPjgkdIdnhAba3RpTkEYAJ3OiyJySlTwSxGRAhE5KSJ3lVkvQESyiv7fXEQyxAEuAQNAq1bAjTfqvlvS4F//UrlkyRKjK3FgHHCbXEzB7gQkeZuQK16weJR0eIKXl2rmwHbFmggDoNNLEXUG0EdEGtnc/5yIbLL5eYOU7gSytSobNyIAdugA9Oql+25Jg7Q0wNcX6NEDHJqkhixTolHozjOA5CKKOjxZ3Cpo88qRBTQTBkCnlyIqAHYQkV0iEi+qDeAqKX2Gr6uojiNJIrJdRExV2bgRAdBkAspMe0oOIDpaHbdXrDC6EscTHw9EBcUiV9gGkFwAxxbVhTAAkhZGBMCwMKBdO913SxqlpgL16wNmszq+07VlZalZbzw8AHc3C041NqlevvxSJGfGDk+6EAZA0sKIAHjTTWq+SHI8zz6rjtmrVhldiX2zWICVK4E2bdTvKyIC2L4d1xxwm5fFyCmww5MuhAGQtDAiAA4apHqVkuM5dQrw9gbCw3kWsCJ//gkMHaq+2/z8gA8+KNN7upwBt9kwnpwKOzzpQhgASQsjAuBf/qJCBDmmSZPUMXvtWqMrsS+5ucD06eoyuQjw4IPA6dOVPMBmwG1ER/MyGDkPngHUhTAAkhZGBMBRo9TnnmeQHNOxY+pkVb9+fA2LxcSojk0iQJcuwPr1RldEZCC2AdSFMACSFkYEwIceUp/7K1d03zXVkvHj1Wv4yy9GV2Ks1FTg4YfV7+K664D//AfIyTG6KiKDsRewLoQBkLQwIgA+9pj67F+4oPuuqZYcPqx6tt5yi9GVGKOwEPjoI9WWVURNbXj4sNFVEdkRdniqc8IASFoYEQCfekp9/k+d0n3XVIseeUS9jps3G12JvnbuBHr3Vs+9dWvg6695KZyoXOV0eMoTTxS4s8NTbRAGQNLCiAD4wgvqy/PPP3XfNdWigwfVH/FDhhhdiT4yMoCnn1bP2cNDDYmTmWl0VUQOoKjDkyUqCu9eF42/dWebv9ogDICkhREB8JVXVADct0/3XVMtGz1avZbOPLe7xQJ89RUQEKCea9++wO7dRldF5JhGjFB/QGVlGV2J4xMGQNLCiAA4Y4b6It2xQ/ddUy1LTFSv5bBhRldSN5KTgcGD1XP09wfmzeNcyERazJqlPk/sKa+dMACSFkYEwLffVgeA33/XfddUByIj1eu5a5fRldSeK1eAV19V41WKqPaOZ88aXRWR49uyRX2mXn3V6EocnzAAkhZGBMAPP+QQIs5k1y71ekZGGl1J7fjpJ6BzZ/WcevQANm40uiIi55Gbq4ZMuv12oytxfMIASFoYEQDnz1dfrj/8oPuuqY4MG6ZeU0fu1HfqFPDAA+p51K8P/Pe/6suKiGrXwIFqBJj8fKMrcWzCAEhaGBEAly5VX7Jff637rqmOxMaq13T0aKMrqb6CAuC994CGDdVzGD4cSEkxuioi5/Xyy+qztn270ZU4NmEAJC2MCIDffqs+/IsX675rqkODB6shUg4eNLqSqouLA0JD1fuxbVvgu++MrojI+a1dqz5z775rdCWOTRgASQsjAmDxh//TT3XfNdWhzZtLOkzYu4sXgb//HXBzU+PTTpkCZGcbXRWRa0hPV5+9UaOMrsSxCQMgaWFEAIyJUUFhzhzdd0117JZb1Bhf9jotmsWizjy3aKHegzffDCQkGF0Vkesxm4GWLTmLjhbCAEhaGBEAt25VX76zZum+a6pjv/yiXtvx442u5GoHDgC33abqa9oUWLCAXz5ERpk0SX0WDx0yuhLHJQyApIURAXD3bvXB/89/dN811TGLBejXT039efSo0dUoly8D//wnUK+eet899hhw/rzRVRG5tuLOgJ9/bnQljksYAEkLIwLggQPqg//yy7rvmnRQ3MZz0iSjKwF+/BHo0EHVExQE/Pab0RUREQAcP17yBxnVjDAAkhZGBMCjR9UH//nndd816cBiAcLD1Swap04ZU8Px4yUzlPj4AG+9BeTlGVMLEZWvXTugWzejq3BcwgBIWhgRAFNT1Rfz3/+u+65JJ6tWqdf42Wf13W9+vppq0NdX7X/kSODYMX1rIKKqeegh9TnlNIs1IwyApIURATAjQ33ox43TfdekE4tF9fKrX1+/g/sff6h9igDt23OmGSJ79/HH6vP67bdGV+KYhAGQtKjtAGg2m2E2m9G9e3d4eHhYf77vvvuwbds2jBkzBrm56kM/ZkzVtrlq1SpMnjy5VuukuhMYGIiuXbsiMNAMkS4YOPANALC+/rVp27ZtGDVqDMaPV+8pLy/gpZeAS5dqdTdEVAcSE9kcSAthACQt6uoM4JEjR9CkSZNyl1ksasaIe++tk13XWD4npqwVgYGBSEhIQGEh0LnzSYj44aefYmt9PxaL6kHYtKn6ErnlFmDfvlrfDRHVkcJCwN8f6NPH6EockzAAkhZ6BsCYmBiEhYUBAK677gi8vJrg5ZdfRkhICLp27Ypt27Zh/PjxCAoKQkREBE6ePAkAWLBgASIjIwEASUlJ6NevH4KDg9GzZ0+8XNSVOCsrC48++ihMJhNMJhNeffVV635nz56N8PBwhISEICIiAlu3brX9AOGtt97CwIEDMXnyZMTHx6N///7o1asXunfvjhkzZtTJ78eZFQdAAFiyBBDpjVGjlpd6/YvfH1V9/QcPHowHHngAPXv2RFhYGNauPYwBAwCRGHh6hmHRIhUI16xZg/DwcAQHB8NsNltfaxFBVlaWtcYmTZrgyJEj+v5iiOgqf/mLmo2HZ+2rTxgASQujAqCf3xGICNasWQMAePPNN+Hn54ddu3YBACZOnIipU6cCKB0An376aUyfPt26zbS0NADAlClT8OCDD6KwsBDZ2dkICQnB8uXLAQBnbRqhbdmyBSaTyfqziJTaXmZmJnJycgAAly9fRkhICLZt21YLvxHXYRsAExP3w8urAxo0OIvVq0sHwOq8/n5+fkhJSUF2NtCnTzTc3CbAzQ0YPjwGISFqmwcPHkSLFi1wsGgy4ry8PKSnpwNgACSyV//9rzqDHxNjdCWORxgASQujAmDz5kfg7u5rXbZ+/fpSwezTTz/F/fffD6B0APz666/Rvn17vPTSS1i3bh0KCwsBAKGhodi8ebP18W+//TYmTJgAAFi3bh0GDBgAk8kEs9kMNzc35ObmAlDB4PTp09bHpaam4qGHHkLPnj0RHBwMf39/fPLJJ7X2e3EFxW0Au3XrBjc3N4we/W5Rp5/SAdDXt+qv/9ChQ7FqFRAYCIh8B1/fQYiNLf2emjNnDh599NFyaxIGQCK79NtvnBigpoQBkLQwKgC2bXsEHh5Nyl0GlA59tv8HgFOnTuGLL77AsGHDcNdddwEAevXqhd9sRvl9++238fjjjyM3Nxe+vr7Yvn07ACAjI6NUGJAywSAqKgrPP/+8tT3gvffei/fff1/z78OV2J4B/Pnnn3HdddehVat4NGwYg169Sl8CLlbZ6//WWwsQEBAJEaBBA2D8+NUYMGDgVY+rLAB6eHhYzwYCgI+PDwMgkR3IyVFjhg4ZYnQljkcYAEkLowJg5841C4BJSUkoKCgAAOzfvx/+/v4A1CXgRx55BBaLBdnZ2QgNDcWKFSuQkZEBb29vpKamAgBef/31SgPgyJEj8c477wAADhw4AF9fXwbAarINgADw7LPPIjh4OERi0Lp11QPgvfdGYuZMoF69BRCJxP33AydPAqtXr8bAgQOvelxycjJatmxZ7iXgLl26YN26dQCAb775BiLCAEhkJ/r3V3/cFR3aqYqEAZC0MCQAxsZiVtMn0EC8gehoIDa2ygFw+vTp6NGjB0JCQmA2m/HVV18BUJ1Axo4dW24nkJkzZyIwMBA333wzZs2aVWkA3LlzJ0wmE8LDwzFq1CgMHz6cAbCaygbACxcuoHHjxvDzewuenmG4sjEWR554Ak28K379p05dgIYN1Vm/5s0XoH//kjPAFQVAAPjhhx8QGhqKoKAghISEIDZW9T7+8ccf0alTJ/Tr1w8vv/wyLwET2ZEXX1SXgYuaAFMVCQMgaaHrQNAJCYDJBHh5IV881Sfe01MN3mYyqUGhyGktfjEBCWJCgYeXet3Lef3PngXGjlWL6tUDpk0DrlwxunIiqktr1qjPPP/Wrh5hACQtdAuACQlqUlZ3d/VJL3tzd1fLGQKdU0ICLNf7oEDKf/0t7u7Iq+eDGxsmQgS44w4gKcnooolIDxcuAG5uQFG/L6oiYQAkLXQJgBaLOsNTUfizDYE2PUHJSVTx9S8Qd+z3NGHZMvUQInIdQUFAq1b87FeHMACSFroEwNhYdZmvsvBXfPPyAuLi6r4m0k81Xn+LJ19/Ilc0caI6DLBpbtUJAyBpoUsAjI4uafN1rZunp1qfnAdffyK6BjVrELBokdGVOA5hACQtdAmA48ZV7cu/+BYVVfc1kX74+hPRNRw9qj7+jz9udCWOQxgASQt7OwNo4Rkg58MzgERUBTfcwGbg1SEMgKSFvbUBzBMvLH8hjkN/OBO2ASWiKhg9Wh0Gzp83uhLHIAyApIU99QIudHPHfg8TRFRvsPfeAy5frvvyqI6xFzgRVcEHH6hDwfffG12JYxAGQNJCt3EAExOrNA7gpbhEzJ4NtGih7g4IAN55h0HQ4VXx9ec4kESua88edTiYMsXoShyDMACSFrrOBJKYaJ0JpKKZIIpduqSCX8uWarUWLYDZs4HsbP3KpVpWjdefiFxPQQHg5wf062d0JY5BGABJC10DYLG4ONXQPypK/VtJm6/Ll9Wl4FatVF5o1gx4803AZvpecjTVeP2JyLXcdZf6m5BXfa5NGACd2jQRgYj0LPq5s4j8ISJJIhInIj1s1q1sWYUMCYA1cOUKMGcO0KaNCoJNmwJvvAFkZhpdGRER1Zbp09UxftMmoyuxf8IA6LRCRWStiByVkgD4q4iMLfr/KBHZ/MwTwgAAIABJREFUYrN+Zcsq5CgBsFhODvDRR0Dbtuog0bixOmBkZBhdGRERabVxozq2z5hhdCX2TxgAnZK3qADXXkRSRAXA5iKSLiKeReu4icgZEWl3jWWVcrQAWCw3F5g3D2jXTh0s/P2B114D0tONroyIiGrq8mV1Cfjuu42uxP4JA6BTmikik4r+nyIqAIaJyL4y68WJyIBrLKuUowbAYnl5wKefAu3bqyDo5wdMmwZcvGh0ZUREVBM33qiO5YWFRldi34QB0OncKOpyrlvRzylSEgD3lll3m5QEwIqWlfWcqDfNCRE54efnZ/R7uFbk5QELFgAdO6og2LAh8MorQFqa0ZUREVF1TJmijuPx8UZXYt+EAdDpvCgip0QFvxQRKRCRkyLyiIhkSMWXgCtaVilHPwNYVn4+sHAh0KWLOoA0aAC89BJHlicichTff6+O3x9+aHQl9k0YAJ1eipR0AtkgpTt6bLVZr7JlFXK2AFisoABYvBjo2lUdSHx9gRdfBM6dM7oyIiKqzPnz6rg9ZozRldg3YQB0eilSEgC7iuockiQi20XEZLNeZcsq5KwBsFhBAfDll0D37uqA4uMDvPACkJpqdGVERFSRHj3UaA9UMWEAJC2cPQAWKywEvvoK6NlTBcH69YHnngNOnza6MiIiKmvCBHWsPnrU6ErslzAAkhauEgCLFRYCX38NBAerg8t11wH/+Adw6pTRlRERUbGFC9Ux+ssvja7EfgkDIGnhagGwWGEhsHIlEBJSEgSffho4ccLoyoiI6M8/1bF54kSjK7FfwgBIWrhqACxmsQCrVgFhYepg4+0NTJoEHDtmdGVERK7LYlFzwAcFGV2J/RIGQNLC1QNgMYsFWLMGiIhQQbBePeCJJ2re/iQwMBAJCQm1W2Q5zGYzLnPWdCJyQvfdB7i5ARcuGF2JfRIGQNKCAbA0iwX48UegTx8VBL28VGPkI0eqtx29AiARkbN67z11HP7hB6MrsU/CAEhaMACWz2IB1q0D+vVTByBPTyAqCjh8uGqPLy8Alr0vLCwMMTExAIC9e/eid+/eMJlMGD16NPr06YPVq1cDAJKTkzFo0CAEBQXBbDbj22+/tW5DRJCVlaXtyRIR2aGdO9Xxd+pUoyuxT8IASFowAFbOYgHWrwduvlkdiDw8gEcfBZKTK39cdQNgaGgoFi1aBADYvn073N3drQGwd+/emDt3LgAgKSkJjRs3xrGiRorCAEhETio/X83mdPPNRldin4QBkLRgAKwaiwX49Vdg4MCSIPi3vwEHD5a/fnUCYEZGBurVq4dCm5nPQ0JCsHr1amRmZqJevXooKCiwLhs+fDi+LBobQRgAiciJDR6sOufl5Bhdif0RBkDSggGw+jZsAG67TQVBd3fgwQeB/ftLr1NeAOzYsSN2795t/dlkMiEmJgbp6enw9vaGxWKxLisOgOWFw+HDh2Pp0qUAGACJyLm99po61v7+u9GV2B9hACQtGABrbvNm4Pbb1cHJzQ0YPRrYu1ctKy8ADh482HopNzY2Fp6entZLwL169cKSJUsAADt37oSHh0epS8Dz588HABw6dAhNmjThJWAicgkxMeoYO3Om0ZXYH2EAJC0YALX7/XdgyJCSIHjffUCrVoFo0aIFWrdubb1t3rwZPXr0QO/evfHEE0/AbDZbA2BCQgLCw8MRGhqKqKgohIaGYtOmTQBUJ5DbbrsNwcHB5XYCyc7ONuJpExHVuUuXVCe8YcOMrsT+CAMgacEAWHu2bgXuvlsFQREgMhLYs6ecFS0WIDYWiI4Gxo0DoqORvWGD9RLw3r170axZM1y4xuBXqampuO6660pdOiYicjZ9+gCNG6sZnKiEMACSFgyAtS8uDvjLX0qC4L33Art2FS1MSABMJjXAoKendYyZdR4eCPb2RlDnzggKCsJ33313jX3EoX379pgxY0bdPyEiIgM9/7w6VBY3sSFFGABJCwbAurNjBzBiREkQfPKWBBRc56N6jhTfaXtzdwd8fIDERKNLJyKyG99+qw6RRU2oqYgwAJIWDIB1b9cuYOS9FiSICQVSQfizDYEmk9ElExHZjbNn1eHx4YeNrsS+CAMgacEAqJPYWFg8vSoPf8U3Ly91HZlcUmBgILp27Qqz2Yzu3btjzpw52LZtG8aMGQMAuHjxImZq7BJZ3jaioqKsHY+I7E3XrkD79kZXYV+EAZC0YADUSXR0SZu/a908PdX65JJshxA6duwY/Pz8sMemN9GRI0fQpEkTTfuojW0Q6SkqSh0eT5wwuhL7IQyApAUDoE7Gjata+Cu+RUUZXTEZpOwYkhEREVixYgXCwsIAAEOGDIGHhwfMZrP1vtOnT+Ovf/0rIiIiEBQUhFdeeQUAUFhYiEmTJqFr164IDg5GaGgorly5Uu42Bg4caB17ksjeLFigDo3Llhldif0QBkDSggFQJzwDSFVkGwDj4+PRoEEDzJs3zxrUyjt7N3jwYGzcuBEAkJ+fjyFDhmDlypXYuXMnunXrZp1JJj09HYWFheVugwGQ7FlyclFnuieNrsR+CAMgacEAqJPYWNW2j20A6Rps2wDeeOONWLFiBWJiYioMgNnZ2fD09ITZbLbeOnbsiBkzZiA9PR0dO3bEI488gs8//xznzp0rdxsAAyDZN4sFaNECCAkxuhL7IQyApAUDoE4sFtW7t6IhYNgLmIqUN41gZQEwMzMTXl5eyMvLK3d7OTk5WLduHV588UUEBAQgOTmZAZAcUmSkOkSmpxtdiX0QBkDSggFQR4mJapw/jgNIlbhWAExLS0P9+vWRn59vXX7bbbfhtddes/588uRJHD9+HGfPnkVaWhoAwGKx4MYbb8SqVavK3QYDINm7d95Rh8r/b+/O46Is9zaA/1Bww53IpRQNFXAYRllc0tS0LJdU1PJYFAUu7Z3EojL3rHNOm2Ud67VXzSVTK7e3tLRITQ3NFSyDDFLS3HEhZZvr/eOeGQYFBB5mnpnh+n4+85FZn3vueXzmmvu5lw0b9C6JaxAGQNKCAdDJUlNtK4GYLX0C87281Wlfg4Hhj8oOgJZlBMeEhSG4USNENG8OJCfj+PHjGD16NEJDQxEaGopu3bph37592L17N8LDw2E0GmEwGDB27FhbS+GYMWMQHBzMQSDkNnbtUgHw5Zf1LolrEAZA0oIBUCc7dwKJiVjoHY9lAYns80fXV8oygvzxQNVFfr46SdKnj94lcQ3CAEhaMADqq3VrwNIAQ1S6lBR2HyACcMcdQN26QG6u3iXRnzAAkhYMgPoymYBbbtG7FOTSOICIyGbaNLW7//ij3iXRnzAAkhYMgPrq0wdo0kTvUpBL4xRCRDabNqld/Y039C6J/oQBkLRgANRXdDTg5QVY5ukluhYnESeyuXgRqFkTGDZM75LoTxgASQsGQH1ZV4g7d07vkpDL4jKCRMVERgI33KB6R1RnwgBIWjAA6ishQX1n//673iUhl8UWQKJi/vlPtbsfOqR3SfQlDICkBQOgvl55RR3Idu/WuyTkstgHkKiYzz5Tu/tHH+ldEn0JAyBpwQCor/ffVweyTZv0Lgm5LLMZ5o4GFAhHARMBwPHjapePjdW7JPoSBkDSggFQX598og5kK1fqXRJyZQufS8VF8UVhKSHQzHkAqZpp1w4IDNS7FPoSBkDSggFQX+vXq+/wefP0Lgm5quRk1bVvUJtUFIQUXwmkoIY3csUHZ1pyJRB3JiK4ePFisdtKWhLwetasWYOJEydWZdFs/vzzT/TRYQmOkuqhd+/e6NdvHUQmY+7cT51eJlchDICkBQOgvn78UQXA//xH75KQK8rOBtq2BWrXBvbts9xoWUYQ8fG48mwi+jXciYAAwLLEL7khqUAAzM/Pd1axXEJpAfCpp9aVevaksLAQhdVgbi1hACQtGAD1deiQCoAvvqh3ScjVmM3Avfeq/eO//y39cbNmsUO8u5PrBMCAgAC88sor6NOnD+6//35Mnz4dJpMJJpMJoaGhEBFkZmZiwYIFGDFihO015s+fD5PJhLCwMERERCAjIwMAsGHDBvTo0QPh4eHo0qULNm/eDABISkqCyWTCY489hrCwMHTs2BG7du0CAGRkZMDPz88JtVFcaQFw7tx1EIlFnz5zAABTp05FTEwMoqOjYTQakZWVhSeeeAJBQUEICwtDeHg4Ll++DABYtGgRQkNDYTQaMXDgQGRlZQEA2rdvj59++sm2nfnz5yM6OhoA8OabbyIyMhKdOnVCVFQUfnSBpUiEAZC0YADU14kT6sv70Uf1Lgm5mg8+UPvGiBFlz3d2/rxaTaZtW7YCuispRwAcN24czCXsCGPHjkVcXBwAFAuASUlJCAwMxLFjxwAAOTk5yMnJweHDh9G9e3ecP38eAJCeno6WLVsiLy8PSUlJ8Pb2toW+uXPnon///gD0DYBBQUG2wGsymeDr64u1a9ehdu1YtGpVFABvuukmnDhxAgCwZ88eBAcH21oCs7OzUVhYiJSUFDRr1swW+l555RUMHDgQADBr1iw88cQTtm336tULa9euBQCcPHnSdvuOHTtgcIEBV8IASFowAOorL099yY8apXdJyJUcOADUqQO0aVO+ScJnzFD70fz5ji8bVT0pJQCmWvp1BgQEYMeOHdc8b+bMmbjrrrtsp4XtA+DEiRMxffr0a57z/vvvw9/fv1igatmyJQ4fPoykpKRiwWbfvn0ItIy0cLUWwHXr1qF161h4ec3BhQsqAI4fP972mOzsbAQGBiI2NhYLFy7EqVOnAADvvvsu4u0mSz979izq1KkDs9mMrKws+Pv7Izc3F7/99huaN29uq9uvv/4avXr1gsFggMlkgpeXF3Jzc51QA6UTBkDSggFQf76+gOVHNhEuXQKCg9U4j/KeZcrOBho3VqMiq1kXMY/g7+9vOz1r5evra2t1KikELVq0COHh4cWCY3kC4HvvvYcHH3ywxHIkJSUhIiLCdj0lJQUBAQEAXDMARkbGQmQONm5UATAhIaHY465cuYKvv/4aL7zwAlq0aIH09HS88847GDNmjO0x586dswVAALjrrruwcuVKvPzyy3juuecAALm5uahfv77t9PD58+dLDO3OJgyApAUDoP5uvhmIitK7FOQqHnmkcgODpk1Tz/v4Y8eUixwnJiYGzz//vO36xx9/jM6dO9uuXx2CNm3ahHbt2tlO71rZB8DNmzcjMDAQx48fB1B0CjgtLQ3+/v7FXi85ORmA+wXAQYNUAJw69doAePLkSZw5cwYAYDab0b17d6xZswapqalo2bKlrV5ee+01DBo0yPa85cuXY8CAAWjdujUOHjwIQAW+2rVr204vv/LKKwyA5P4YAPVnNKo5rYgWL1YhbsAAoKKDGM+dAxo2BNq3Zyuguzl9+jRiYmJgNBphMplw9913Iy0tzXb/1SGoT58+aNmyZbHTuH/++ec1g0AWLFgAo9GIsLAwREZGIjMzE4A6ndm1a1eEhYUhODgY999/PwD3C4APPhgLH5856Nfv2gC4e/duhIeHw2g0wmAwYOzYscizdJL9+OOPYTQarxkEAqhWw6ZNm6Jr167Ftvnvf/8bAQEBuO222/D6668zAJL7YwDUX69egA7HVXIxv/6qugO0aAHY9TevkMmTVYBcvLhqy0Y6MpvVZJCJiUBcnPrX0mJX7ZnNeCwiGW/4JKLg4epXN8IASFowAOpv6FCgZs2yR3qSZ7t8GejUSa3mlpRU+dc5cwZo0AAICgIKCqqseKSXlBS1vJ/d5N/w9lbXDdV88m9L3RTU8EGeVM+6EQZAj/SNiBwQkX0islVEOlluzxSRQ5bb94nIKLvntBeR7SKSJiI7RaRjeTbEAKi/hx9Wxy7LrAxUDT35pNoHpk7V/lqTJqnX+uQT7a9FOkpJUU3CNUpZA7o6L//HugHAAOipGtv9PUxE9lj+zhSR0FKe852IPGz5e6SI7CjPhhgA9ffss+qYZemeQ9XMF1+oz79376pptTt9GqhfHwgJYSug2zKbVStWaQHHPui4wHx0TsW6sREGQI8XKyI/Wf7OlJID4I0iki0i3pbrXiLyl4i0ud6LMwA6XkBAAPz9/W0dkAHg22+/hYggISHBNofb3r3XPrekqQ0qIzY2FnPmqAlTJ0+ejE8/rb7rZ7qSzEw1fcsNNwB2/dA1e+EFtU/xY3ZTycnqVGZZAcd68fFRywNWF6wbG2EA9FiLROSo5WKw3JYpIvtFJEVEPhIRf8vtESLy81XP3ykiva63EQZAxwsICEBERAQ+++wz220PPPAAIiMjkZCQgDlz1LHqu++ufa4jAiC5hrw8oHt39dl/+WXVvvapU+oMmMFQ8dHE5AISE4v6/F3v4u2tHl9dsG5shAHQ48WKyFeWv1tb/vURkX/b3R4hIgevet4uKTkAThC102SJSFajRo303oc9XkBAAN5//33bckPWGepffvllJCQkYMkSQOR1tGsXhc6dO2PAgAE4cuQIABUAR48ejcGDByMkJAS33367bW6rAwcOoGfPnujcuTNCQkLw6quv2raZlZWFvn37wmg0YsiQIRg0aJAtANqHwU2bNqFbt27o1KkTDAYD5nMpCadJTFTfUVWQ70v0/PPq9VeudMzrkwPFxZUv4FgvditbeDzWjY0wAFYLl0XE76rbWojIRcvfN4rIeeEpYJdknccqODgYWVlZmDt3Ll544QVb695zzy2FyFjMm6c6bC1atAhDhgwBoALgLbfcYgt9o0aNsgW9Cxcu4MqVKwCAv//+G506dbKt4Tl8+HBMmzYNAHD48GHUr1+/xAB49uxZFFg6ip05cwYBAQHXTC5LVW/DBvXd1KUL4KjVpE6eBOrVU/NMshXQffz6K/BVWGLRyNZq3sp1DbYA2ggDoMdpKCIt7a5Hi/qQfaX44JAJIrLF7vr3UnwQyI/l2RgDoONZA+CsWbMwa9YsREVFIS0tzRYA+/a9FyJt0bKlmtA1NDQUJpMJgAqA9ouTz54927aO5YkTJxATE4PQ0FCEhYWhSZMmmDdvHgCgSZMmxSY3HTp0aIkBMD09HdHR0bb1LX19fbFx40an1Et1dewY4O8PNGoE/P67Y7c1caL6Hvz8c8duh7RLTgaGDwe8vIAoSUaeV/n6uRXU8MGFbz23n9s12AfQRhgAPU4rUf33UkT199skahqYW0Rkr6jpYVJEZI0Ub+ELEjXyN03UoBGDlAMDoONZA+Cff/6JZs2a4bbbbgNQ1L+vf/+REPlfTJp07XOv7gM4Z84cxMbGAgDi4+ORkJBgW6w8OjraFuzKGwD79euH2bNn29bB7Ny5M9atW1e1FUA2BQVA377OOzX7119A3bqAycR5Jl2R2axag2+/vSizDB4MbN1SvpGuBVIDKWLAjTeqyb+rxWfMUcA2wgBIWjAAOp79Ukbz58/H9u3bARSFu3ffXQSRCMTFqdO8eXl52LNnT7HHWNkHwOHDh+Ptt98GABw6dKjYad4RI0ZgxowZANQSTqWdAg4PD8eqVasAqLVDa9SowQDoQDNnqu+mRx913jat0wytXu28bVLZ8vPVPI2dOhWdqXzoITW9nU1q6nXnujP7+mL1rFTccIO6qXdvwLJ8rWe7Tt2YvTgPINF1MQA6XklrWQKWcDdhAvI+nI9npDVCatRCWL166NiyJd58882ix5QSAPfs2QODwYDIyEiMHDkSQ4YMsQU7+0EgI0aMwLBhw0oMgN988w3atWuHrl274pFHHkH37t0ZAB1kyxb1fRUWBvz9t/O2e/w4UKcO0LlzNWkhcmE5OcB77wFt26qsUq8e8MwzwB9/lPKE1NRyrQRy5oz6UeHlVdTt7dIl570vXZRQN+aa3sgVH2T4ciUQoutiANTRF18AtWurA9fVv2Jr12aTjQc5dQq46Sb1hf/LL87f/jPPqN1q7Vrnb5tUQJs5U/X9FFFrf0+bpibtLpedO1Wqi49X/5bSry05GQgPV9to3RpYtaoahP6r6iax706IAD//rHfBHE8YAEkLBkCdWJd/uN6FIdDtmc2qX5cIsGCBPmX480/1myIyshoEAhdy5Ig6Be/rqz7/gADg3Xcd2zpXUKBaGRs2VNscNAg4fNhx23M1W7eq9/3YY3qXxPGEAZC0YADUQWGhreXvupfatfUuLWn01lvqo4yJ0Td8WdcbrupJp+laBw8CsbFFZ22NRmDJEjX5t7McP672ORHVBWDmTMAya5RHM5tVK2i9esDZs3qXxrGEAZC0YADUwfz5FZvIdOFCvUtMlbRrl+qi1KEDcOGCvmU5ehSoVUvNPchWQMfYtg0YMqTov26vXsBXX+lb3999BwQHq/J06ABUh1mePv5Yvd/XX9e7JI4lDICkBQOgDrp2rVgA7NpV7xJTJWRnA7fcohpxS1rnWQ+PP652qfXr9S6JNgEBAQgKCoLJZEJwcDBGjx6NS5U8r5qRkQE/P79Kl6WwEFi3DujZs+i/7LBhwI4d5Xv+1KlTkeuo2cAtcnOB115TUwKJAKNGqW4BnurKFaBZM3XK3TJLlkcSBkDSggFQB+3bVywAduigd4mpgsxm9SUrovpjuYo//lAtkt26uXcroP3IerPZjEGDBuG9SlZ0ZQNgXp5qaTIYrINz8xEXV/FBPiKCixcvVnj7lZGZCQwdqsrboAEwe7bnBqSpU9X79ORJ0IUBkLRgANQBWwA93v/8j/rooqNdL2iNH6/K9s03epek8uwD4OXLl9GnTx8sXboUCxYswIgRI2yPW7duHXr37m27Pn/+fJhMJoSFhSEiIgIZGRm2ADh58mSEh4cjMDAQX9p1lHzggQcQEREBo9GIQYMG4fffT2D2bODGG5MgYoKPz1No0aIbPvhgBS5cuIAxY8YgKioKRqMR48ePR56l49/MmTMRHBwMk0mt+JOZmYnx48dDRGA0GmEymXDixAmn1N/atUCbNmo/MJnUqWtPc/y4+rHTq5feJXEcYQAkLRgAdcA+gB7twAHV6T4gwDU7oWdmqi/GHj1cL5yWl/0p4IYNG+L2229Hfn5+mQEwKSkJgYGBtrWuc3JykJOTg4yMDIgIVltG3K9fvx4d7FrdT506ZfkX6Nv3NdSu/QREgMaNk+Dl5YWvvtpqe+zYsWOxaNEiAKplMj4+Hm+99RbOnj2LRo0a4W/LBJA5OTm4fPkyAOe2ANrLyQEmTSpaVS0+Xr1HT/Lgg+q9uUoXjKomDICkBQOgDjgK2GNdugSEhAA1awKWBV9c0tixavfatEnvklSOfQtgfn4+4uLiMGHChDID4MSJEzF9+vRrXisjIwO+vr6269nZ2ahZs6bt+pQps3HjjRHw8gqFSCBq1+6B//4X2LAhCUFBQcVey9/f39aaZzKZ0KFDBzz22GMoKChAZGQkhg8fjg8++ABHjx61PUd0CoBWhw4B/fqp/aFpU2DePHWI8gQ//aTe18MP610SxxAGQNKCAVAnq1dzHkAPFBenPrZ//UvvkpTt99/VFCXuenrs6tV1vvzyS3Ts2BGLFy/GsGHDbLevXLmyXAHQvg/gxYsXISLYvx+4886tEGkHkZPo3BmYOHENwsMjAKgWxYiIiGKvdcMNN+BwKZPuFRQUYMuWLZgxYwZuuukmbNmyBYD+ARBQLcHLlgHNm6v9t1s3z2k169FDjX530tl1pxIGQNKCAVBHq1cDtWvDLFwJxBMsXao+urvuco8WFGtYTUrSuyQVd3UATEhIQHR0NLZv3462bdvi8uXLyM/PR3R0tC0Abt68GYGBgTh+/DiA4qeArQHQbAbWr1cBUP1XXIsGDcKxfn0BrlzJxcCBA22hr6QAGBcXh7FjxyLfMrLi7NmzSE9Px4ULF/Cn3bDbUaNG4Z133gEANGjQAFlZWY6pqArKzlarxtSooS7PPAOcP693qbRZsULt55al0T2KMACSFgyA+rv8wUJsl6441qCDGvDBPn9uJy0NqF8faNHCfVoaDh9Wp6r79NG7JBVn3wewY8eOGDhwII4cOQIAePTRRxEYGIh+/fph4sSJxQaBLFiwAEajEWFhYYiMjERmZqYtAK5apVq+RFQAHDkS2LEjH/fddx/atWuHPn364KWXXiozAF64cAGPPfYYDAYDjEYjwsPDsXHjRhw9ehRdu3ZFaGgojEYjhg8fjuzsbADAtGnT0KFDB6cOArmevXutdaH26WXL3Le/aH4+0KqVat108Gw7TicMgKQFA6D+jhxRB9pnntG7JFQZV64AnTsDXl7At9/qXZqKefhhte9t3qx3SaqQ2awWxU1MVM2ciYnqeglyc4H//d+iiZJr1QLGjVOBvrorLFT9AZs2VXXTr5/qL+iO/vUv9R6WLNG7JFVLGABJCwZA/e3frw5O06bpXRKqjKeeUp/f5Ml6l6Ti0tLUqb5+/fQuSRVJSVET8/n4FK3D5u2trhsMQGoqALUqyxtvAC1bqoc0bKhyomWAMNk5dUqNEBZR1ThpkhpB7E7OnFGTYEdFuW9LZkmEAZC0YADU3/ffq4OrpUsQuZFVq9Rnd9tt7juhrnWqjB9+0LskGqWkAL6+KtGWNKCqRg0U1vPF7LGpaNxY3dS8OfDvf6u+b1S2bduAsDBVb23aqNVP3Mm4carsrjw6v6KEAZC0YADUnzVEWKYPIzfxxx9AkyaAn59aZ9dd/fqrykx33ql3STQwm1ULX2nhz3IpkBpIEQPat1enN69c0bvg7iU/H3j7bdXf1brkXWam3qUqn9TUomXwPIUwAJIWDID6s84L7W6/qKuzvDzg1ls953N74AE3bx1JTi6a0fg6l8KaPijYsVPvEru1rKyipQ7r1VN97NxhgMUdd6iBT+78g82eMACSFgyA+nvzTQ85BVeNvPii+swmTNC7JFXj55/VIJa779a7JJWUmFjU5+96F29v9XjS7JtvipY2Dwlx/SmF1q1TZX3pJb1LUjWEAZC0YADU38svq4PSwYN6l4TK45tvVFiKinJCfxYfAAAetElEQVSPVo/y+sc/1H5YyoBZ12ad1LC8l/h4vUvsMa5cAWbOVMsfigAxMWodXldUWAi0a6e6bVhW5XNrwgBIWjAA6u+JJ9SBkyMQXd/x48CNN6pRo6Us+OC2Dh5UwXbgQL1LUglsAdTd4cNq3xEBGjUC3nsPKCjQu1TXeucdVcaPPtK7JNoJAyBpwQCov/vvVwcky9rw5KIKCorWTF2+XO/SOMZ996n3t9PdushVoA8gfHzc8A26B7NZDWpr1UpVdXi467Uonz8PNGgAGI3uPyWMMACSFgyA+hswQJ0+Idf2yivqS23cOL1L4jgHDqj3eM89epekgso5Chg1aqjHkUNdulTUKOvlBTz6KHD2rN6lKvL002p3+O47vUuijTAAkhYMgPrr1k0tt0Sua+tWlR1CQz2j71BZRoxQX467d+tdkgpKTQV8fVHoVfo8gPD1tU0GTY538CDQu7eqfn9/tcqlK7S6paerYDp0qN4l0UYYAEkLBkD9BQcDHTvqXQoqzenTwM03q+kuqsNAHevKNG755Ziait/qGJArPjCXsRIIOY/ZDCxerPrOWidNT0nRu1TA4MEqBLpzX15hACQtGAD116wZ0KOH3qWgkpjNwJAh6otr/ny9S+M80dHqPe/dq3dJKmbfPlXumUN2qnOQ8fHqX/b50925c2rAm5eXmotv4kTg4kX9yrNxo9pXnn1WvzJoJQyApAUDoL7MZrUA/eDBepeESjJ7tvqSeOAB1zh15Sx79qj3PXy43iWpGOu6zG47oXU1sGsXEBmpPqebbwY++0yf/1vWbqMNG6q1od2RMACSFgyA+srJKZo7i1zLrl3qzGH79u77BaGFteVz/369S1I+ly+rpflCQqpXWHdHBQXA3Lmwrcl8992qX56zffih2v6cOc7fdlUQBkDSggFQX1lZ6gD01FN6l4TsnT8PBAaq1tk9e/QujT5++kntmyNH6l2S8lm2TJX3jTf0LgmV14kTQGys+txq1wamTXPudFg5OUDTpupHXmGh87ZbVYQBkLRgANRXSoo6+E2ZondJyMpsLloV49139S6NvgYPVvXgCp32r+eOO9R4jxMn9C4JVdTmzep0rIj64bVhg/O2nZiotvvll87bZlURBkDSggFQX1u2qIPP22/rXRKymjdPfSbDhvFU4s6dqi5GjdK7JGXLyFDlHDFC75JQZeXlAf/5jxptb215PnrU8ds9ckQNSunf3/HbqmrCAEhaMADqa80adbBbuFDvkhCgZgmpWxdo3Ro4c0bv0riGAQPUyE1XngJnyhT1/+irr/QuCWl15IgafCQC1K8PvPmmCoeOdO+9ans//+zY7VQ1YQAkLRgA9bVwoTrwrFmjd0koJ0fNx1izJrBtm96lcR0//qj20dGj9S5JyQoK1NJjN9/smmvPUuV8+SVwyy1q3zMa1WTsjvLDD2o7jz7quG04gjAAkhYMgPp6+2114NmyRe+S0Jgx6rN47TW9S+J67rpLtQL+8oveJbnWhg3qc3v5Zb1LQlXt779V626tWuozfuQR4OTJqt+O2QxERKjTz660ZN31CAMgacEA6BgBAQEICgpCWFgY2rVrhyFDhmBbCc1K1lNXpXWyz8jIgJ+fX4W3n5GRgQ8//LDCz6uuPvlEfQ79+7vnaEBH27bNdacrsp6+c+cVHahsaWnq/6aImurnww+r/v/pokXq9V9/vWpf15GEAZC0YAB0jICAAKTYpbrVq1ejUaNG+PHHH4s97sknCyFSiKyskl+nsgEwKSkJERERFX5edZServoaNW8O/PWX3qVxXXfcoZbT/fVXvUtS5NQpNVdj3756l4QczWwGVqwAWrZUQa1Ll6pdr/rKFbUqU0AAkJ9fda/rSMIASFowADrG1QEQAF588UWMHDkSU6dORUxMDKKjo9G4sREiWdi6dRe6desGo9GIqKgo/PDDDwCKAmBCQgK6dOmCjh074ttvv7W95oYNG9CjRw+Eh4ejS5cu2Lx5MwAgKCgIdevWhclkwj333AMAmDhxIiIjI2EymdCrVy+kpaU5qTZc15UrQHi4Or25aZPepXFtW7eqL96HHtK7JEWsXSiWLtW7JOQsFy4AEyaovro1aqg5VLOzq+a1p01T+9Pnn1fN6zmaMACSFgyAjlFSAPziiy8QEhKCqVOn4qabbsKJEycweDDg45OLVq1aYYNl8qutW7eiefPmuHTpEjIyMiAiWGgZJrxjxw40a9YMly5dwuHDh9G9e3ecP38eAJCeno6WLVsiLy+vxBbAU6dO2f5etmwZBg0a5MgqcAvPPKMO+JMm6V0S99C3r/ri1WPVhqtZl/Jq3Fj1FaPqZf9+4NZb1f/fZs3UjwCt0zb99Zfqb9irV9WU0dGEAZC0YAB0jJIC4Oeff24LgOPHjwcA9OgBNG16AIGBgcUeGxYWhu3btyMjIwO1atVCoV2HF5PJhG3btuH999+Hv78/TCaT7dKyZUscPny4xAC4dOlSdOvWDQaDASEhIajun711Cp6ePd3nlI/evv++qDO+3qyjk594Qu+SkF4KC4H58wE/P7Uv3H679qlcHnpIvdbevVVTRkcSBkDSorqHAEcpKQC+8MILuPfeezF16lQkJCQAUNOOtGmzH+3atSv22LCwMOzYsaPUALh9+3a89957ePDBB0vc/tUB8I8//oCfnx8OW3rK79+/v1J9Cz3FH3+ozuRNmzpnsllP0ru3agXUe9DF2LHqi7q6LtVHRU6fBsaNU/uDjw/w4otqWqfK2L1bvc7DD1dtGR1BGABJCwZAxyhtEMiOHTuKBcAWLYCuXdUpYGvfvm3btl1zCnjRokUAgOTkZNsp4LS0NPj7+xfbTnJyMgBg9+7dxULlgQMH0KJFC+Tk5MBsNmPs2LHVNgDm56uWVxFg7Vq9S+N+vvtO1V18vH5luHQJaNAA6NxZvzKQ69mxA+jUSe2fAQGVn1+1Z091KtjVlxUUBkDSggHQMeyngQkMDMQ999xjmwbGPgDWqaNWWti5c2exQSBbLbOeWgeBTJkypcRBIF9//TW6du2KsLAwBAcH4/777wcA5OfnY9CgQTAYDLZBIE8//TTatGmDXr16YebMmdU2AL70kvqC+Oc/9S6JezKbgdtuU+vuZmToU4YFC9Rn+N57+myfXFd+vlrDu2FDtY/cc0/59lPrMdtkMiEgwAQRE5588iBWr16N4OBgmEwmHDhwoEJlmTx5Mj799NPKvZFyEAZA0oIBUAdmM5CcjPyERHwkcVjTMRGwtNyRY23cqEb8RkSoEcBUORs3qi/XceP02X7PnkDt2u41aS8517FjwP33q/20bl1g1qyy/8/bn7XJz1eryzRvDvTvfzdWrFhR4e3nO6FjsTAAeqRvROSAiOwTka0i0slye3sR2S4iaSKyU0Q62j2nrPtKxQDoZCkpauiijw/MNb0BERR4eauOKwaDWoyWHOKvv9RowQYNgN9+07s07s1sViMwfXxUf0pnOnRIfalbGruJyvTtt0BQkNpngoJKn+7p6m47//oXIPIUatf2Rdu2bdG9e3cAwPr169G5c2cYjUb06tULBy2LZCclJcFkMuGpp55Ct27dsGLFCsTGxmLOnDkAgLy8PCQmJiIqKgomkwmjRo3CuXPnNL03YQD0SI3t/h4mInssf38nIg9b/h4pIjvsHlfWfaViAHSilBTA11dNXiVy7aVGDXU/Q2CVKywE7rxTVfOyZXqXxjN8/bU+66c+/7zarl1PCKIyXbkCvPqqagm0rmt97Fjxx9ifAjaZTAgNNaFOnVw0aNAba9euAwCcOHECfn5+tlPBS5YsgcFgAKACoJeXl637DoBiAXDWrFmYOXOm7b4ZM2bg6aef1vS+hAHQ48WKyE8icqOIZIuIt+V2LxH5S0TaXOe+MjEAOol10rLSwp99CLQcUKjqvPqqqt6xY/Uuiecwm4Fu3VQr4JEjztlmXp5qxb3lFi7ZRxWXkaH6BIqoPoLvvFM0BVRJMzeokcW98frrKgCuXbsW/fr1K/aYRo0a4dixY0hKSkJQUFCx++wDYFRUVLGAGRISgsGDB2t6P8IA6LEWichRy8UgIhEi8vNVj9kpIr2uc1+ZGACdJDlZfVOWFf6sFx8fYOdOvUvsEvLz8zF9+nQEBQWhY8eOCAoKwtixYyt06uSHHwCRAWjf/rdKTw1BJVu/3rlz8a1erbb3yivO2R55pjVr1ChhETVqeMeOkgNgaqoKgLfdts7yvDW44447ij2mUaNGOH78eIlzr9oHwMjIyGID+KqCMAB6vFgR+UpUyDt41X27pCgAlnbf1SaI2mmyRCSrUaNGVbpDUikSE9WwyXIEQHNNbxRMTNS7xC7hoYcewuDBg3HW0tu/sLAQK1assM1naFVah+szZ1Rn7rp1AUtXHapCZjMQFaWmzChtPeuqdM89qpHcGdsiz3bpkpov0Pq7vH79AGzdmnLN4xo37o0aNdbh6FHg5MmTuOGGG/CzZbbpZcuWFTsFXFYAnDFjBgYNGoQcy6/QnJwcpGrs7iMMgNXCZRFpJiLnpfRTwKXdVya2ADpJXFz5Wv8sl3kSj7Ztgf79gccfB956S81Z98sv1Wf0anp6OurWrVtsCTurkjpcL126FF26dEGnTp1gMpnwf//3JYYMUVXq51f06z49PR39+vWD0WiEyWTCqlWrAAAzZ87Ek08+advGxYsX0aRJE5w6dQoHDhxAz5490blzZ4SEhODVV191TiW4gf/7P1XHTz3l2O0cO6YmoB440LHboerl55/VCiIiAahZMwitWhWtrLRlyxaEhvaGyFq8H5sMJCZiff/+6NysGcLatbtmEEhZATAvLw+TJk1CaGgojEYjjEYjlixZoqnswgDocRqKSEu769GiPmQvEfleig/0+NHucWXdVyoGQCepQAtggZc3Pm+fCKOxqNOy/cXLS52+uOMO1QH/jTfUqbGDB4HLl/V+o1Vn+fLlCAsLK/G+kjpcnz59GmbLYqAZGRlo2LAFRPIwenTx0ztdunTBhx9+CABIS0tD06ZNceTIERw9ehT+/v7Izc0FAMyfPx/Dhw8HAFy4cAFXLMn777//RqdOnbBr1y7HvHE3YzYDkZFqWpY//3Tcdl57Te3/n3/uuG1Q9WQ2q7WEmzVT+1iPHmqtYQAo3J+CX30MyBUfmK3HcG/XmLlBGAA9TitR/fdSRGS/iGySomlggkSN7k0TNTDEYPe8su4rFQOgk1SyD6DZrL5Uv/8e+Ogj4IUXgJEjAZNJDRguKRy2aqV+0Y4bB/znP8AXXwAHDlR+aSS9XC8AXt3heteuXbjrrrtgMBjQoYMJIj5o3Tod588XBcALFy6gVq1aKCgosD1vyJAh+OSTTwAA/fv3x8qVKwEAt912G9atKxr9FxMTg9DQUISFhaFJkyaYN2+eI962W1q71rGTa5vNQPv2gL8/YMnnRFUuO1u1ZNeooVqb//NQCsz1fFHo5ZozNwgDIGnBAOgkDhgFbDYDx48DW7aoBdFfegm47z61PFaDBiW//E03qbVcx4xR81x99hmwb5/qD+NqrKeAT58+fc19JZ1uCQwMxKpVq3DhAtCuHSDSBJ9+qlr9rAHw/Pnz16ytPGTIECyzzA2zbNkyDBo0COnp6WjRooUtKMbHxyMhIcHW1zA6Otp2aofUvti5s1rZ5vjxqn/9zZvV/mtZQIfIoXbvBrpEmZEiBhSI687cIAyApAUDoBOlpjptHkCzWa1j+cMPwMKFwMsvA//4h1oBw7pE0tWXFi3UEl+PPKKmTVmxAtizB7hwoQreeyXFxMRg6NChtlG/ZrMZH3/8MT766KNrAmDTpk2xZ89ejB4NiCyGiNhO+159Cnj+/PkAgN9++w1+fn44YpnH5MqVK/D390d8fDyef/5522sPHz4cb7/9NgDg0KFDqF+/PgPgVVatUvvRhAlV/9oPPaRe29L3nsjhCncko6Cma8/cIAyApAUDoJOlptpWAoFO/UnMZuDUKTX1waJFwJQpalWFqCigceOSj2/Nmql+MbGxagqOTz8FfvpJnTJxpLy8PEyZMgUdOnRAx44dERISgvHjx2PVF18gIjhY9a2MiwMSE7F42jTccEMbiPRAu3YT0bp16xIDYHp6Ovr27YuwsLBig0CsHn/8cYgIfvnlF9tte/bsgcFgQGRkJEaOHIkhQ4YwAF7FbFZdE+rWVauuVJXsbPWaloUYiJyjAv224e2tHu9kwgBIWjAA6mTnTnXAiI9X/7rQvH9nzqgui0uWANOmATExQNeugJ9fycc+f3/15fzgg8D06cAnn6i347B1Wu2W07MP0WZvHxz0MqBf81ScOeOgbVOZPv9cfRzPPVd1rzl3rnrNjz6qutckuq4KztyA+HinF1EYAEkLBkCqiLNngV27VMibMUOdmrv1VhUCSzomNm2qwuMDDwBTp6pQ+eOPwOnTqsWowq6znF6B1EBBHS6np5fCQiA0FKhXDzh5smpeMzJSfeR6dkWgaogtgOTpGACpqmRnq87Ty5er08QPPwz07Fk0tcLVl8aN1Wnn0aOByZOBjz8Gtm9XwaHEcMjl9NzCihXqY6iK78P9+9VrxcVpfy2iCnGD1ZuEAZC0YAAkZ7hwAdi7F1i5Ug0wiYsDevVSA09KOp42bKgGrIwaBUyaBCxYAOz7n2SYvV37gEyqFdBgUK12JczhXSFPP60+zm3bqqZsROXmBj84hQGQtGAAJL1duqRaej7/XE1NM2YM0KcPcPPNxY+zr0ki8sS1T8mQ8umn6mN46aXKv8aVK6oLQXBwJbsLEGnlxJkbKkMYAEkLBkByZTk5qtvfqlXA/kjX75RNSkGBCm4NGqDSA3KsIfL116u2bEQV4gIzN5RGGABJCwZAchtu0Cmbiixdqj6Kl1+u3PPvvFN9jCdOVG25iCrFBWduEAZA0oIBkNyGG3TKpiIFBUBQkOrPWdEpgTIz1bKGlqWYiagEwgBIWjAAkttwg07ZVNzixeojmTq1Ys+bOlU978svHVEqIs8gDICkBQMguRUX75RNxeXnq3WZGzUCLKv5XVdBAdC6tVq32rIUMxGVQBgASQsGQHI7Ltwpm661cKH6iKZPL9/jv/5aPX7SJMeWi8jdCQMgacEASG7LBTtlU3FqLedp8PEJQo0aHWE0dsLQoUOxd+9erFmzBhMnTrzmOffdpwLg4cOV2+bcuXPx1ltvaSw5kesTBkDSggGQiBzlgQcewNChQzFnzlmIqBVi1q5diyVLlpT4+NOngVq1gNtvL/018/PzHVRaIvciDICkBQMgETlCWloa6tWrhzNnziAvD2jbVk3sbF3Td8GCBRgxYgQAICkpCREREZg92zLp92spCAgIAABkZGTAz88P06dPR8+ePfHWW28hOzsbw4cPR1BQEG6//XbExMQgISEBADB16lTb3/bbAIB169ahd+/eTqsDIkcSBkDSggGQiBxh+fLlCAsLs12fN0+Fu1dfVddLCoChoWrAyK5dxQOgiGDp0qW215owYQLiLRN9nz17Fm3atGEApGpHGABJCwZAInKE5cuXw2Qy2a7//PNv8PExoUaNDoiNHXNNAAwOjoAI8PjjQEpK8QBYp04dmO3Wg+vcuTN++OEH2/V//vOfDIBU7QgDIGnBAEhEjmA9BXzWbhboDz8ERBYgNHREsXC2detW+Pl1ggiwezewa9eua04B2+vUqRO2bdtmu15aAFy8eDGGDRtme9zKlSsZAMljCAMgacEASESOMnr0aERHR+OcZRLA3FygceP3UKvWCHzwQVEATEvLgkgjhIaeBgA8/fTTZQbAZ599FuPGjQMAnDt3Dm3bti0xAG7fvh1t27bF5cuXkZ+fj+joaAZA8hjCAEhaMAASkaPk5uZiypQp6NChA0JCQnDrrbfCZBoGkR0Ye8dkjOjQAYiLw4GBibhJxqNp0zbo06cPpkyZUmYAPHfuHIYNG4aQkBAMGDAAMTExmDFjBoDiARAAHn30UQQGBqJfv36YOHEiAyB5DGEAJC0YAInImXJ3p+CQtwF54gOzZSLvfPFGrvigILh8E3nn5eXh8uXLAIDz58/DaDRi48aNji46kUsRBkDSggGQiJwmJQXw9UWhl7al/E6cOIHw8HCYTCZ06NAB08u7zAiRBxEGQNKCAZCInMJsVkv1lbaOs30INBj0Li2RyxMGQNKCAZCInCI5Wa3XXFb4s158fLi0H9F1CAMgacEASEROkZgIWPr8Xffi7a0eT0SlEgZA0oIBkIicIi6ufOHPerGs9EFEJRMGQNKCAZCInIItgERVShgASQsGQCJyCvYBJKpSwgBIWjAAEpFTcBQwUZUSBkDSggGQiJwmNVXN81daCCznPIBExABIGjEAEpFTpaaqFj4fn6I+gd7e6rqhfCuBEBEDIGnEAEhEuti5Uw30iI9X/7LPH1GFCAMgacEASERE5H6EAZC0YAAkIiJyP8IASFowABIREbkfYQAkLRgAiYiI3I8wAJIWDIBERETuRxgASQsGQCIiIvcjDICkBQMgERGR+xEGQNKCAZCIiMj9CAMgacEASERE5H6EAZC0YAAkIiJyP8IA6HHqiMhqEUkTkX0iskFE2lju+15Efrfcvk9EnrV73o2Wx6aLSKqI9CzPxhgAiYiI3I8wAHqcOiIyUES8LNefFJFvLH9/LyKDS3nefBGZZvk7SkT+EBHv622MAZCIiMj9CAOgx4sUkd8sf38vpQfASyLib3d9p4j0ud6LMwASERG5H2EA9HiLRORty9/fi8gvIpIiIstF5BbL7X4i8vdVz1shIg9d78UZAImIiNyPMAB6tJdEZIeI1LNcb2X510vUqeGfLdf9RCTnqueulJID4ARRO02WiGQ1atRI732YiIiIKkgYAD3WRBH5SUQal/GYK6LCn4gKgDwFTEREVA0IA6BHmiAiu0Wkid1t3iLSzO76CFEDPawWSvFBIEeEg0CIiIg8kjAAepybRQQicliKpntJFhFfUS2CKSKyX0S+FRGT3fOaiRotnC4iB0Wkd3k2xgBIRETkfoQBkLRgACQiInI/wgBIGuWK3aAQF7lku0AZqtOF9c269tQL65v17cmXAiHyMFl6F6CaYX07D+vauVjfzsX6di7WN3kc7tTOxfp2Hta1c7G+nYv17Vysb/I43Kmdi/XtPKxr52J9Oxfr27lY3+RxJuhdgGqG9e08rGvnYn07F+vbuVjfRERERERERERERERERERE5GzvikimqJVOQu1uby8i20UkTdQaxh3LeR+VrY6IrBZVd/tEZIOItLHcd6PlerqIpIpIT7vnlXUfle0bETkgqr63ikgny+3cxx1nqhQ/prCuHSNTRA5J0QpVoyy3s74do7aIvCdFq3stsdzO+ia31EvUUneZUjwAficiD1v+HikiO8p5H5WtjogMFBEvy/UnRQUUEZH5Unzt6D+kaO3osu6jsjW2+3uYiOyx/M193DHCRWS9qH3UekxhXTtGphQ/bluxvh3jbVGNJtbjdwvLv6xvcmuZUnQguVHU7PHWgOElIn+Jaqkq6z6quEgR+c3y9yUR8be7b6eI9CnHfVR+saLW7+Y+7hi1RX3BtZWiYwrr2nEy5doAyPp2DF9RdVf/qttZ3+T2MqXoQBIhIj9fdf9OUa2FZd1HFbdI1K9KPxH5+6r7VojIQ9e5j8pnkYgctVwMwn3cUf4tIk9Y/s4UdUxhXTtOpojsF5EUEflI1I9E1rdjhInIYRH5l6gfkVtFpJ+wvskDZErxAHjwqvt3SdEOXdp9VDEviWotqScq5OVcdf9KKQqApd1HFRMrIl8J93FH6C7qdJf19FimFAVA1rVjtLb86yMqfHPfdpwIUf1arcddk4icEpGuwvomN5cpxU8Bn5fSm7RLu4/Kb6KoX5H2/dNypPTTvGXdRxVzWUSaCffxqvaCiBwTdSzJFJECEflTVOhmXTteCxG5KDx+O8oNIlIoIjXtbtspIvcJ65vcXKYU70vyvRTvuPpjOe+j65sgIrtFpMlVty+U4gM9jkjRgaOs+6h0DUWkpd31aFFLNXkJ93FHy5SiY8r3wrquar5S/AfkBBHZYvn7e2F9O8I3ogbxiYgEiGoBbCGsb3JT74v6QiwQ9cvEOiAhSNTpyTRRLVUGu+eUdR+V7WZRpxEOS9HUDcmW+5qJOsBYpxjobfe8su6j0rUS9Ss9RVRfqU1SNA0M93HHypSiAMi6rnq3iMheUVMcpYjIGilqWWJ9O8YtogJdiqhjd7TldtY3ERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERVa3/BweTteOYikFZAAAAAElFTkSuQmCC\" width=\"799.9999880790713\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "next_state = cities\n",
    "states = []\n",
    "\n",
    "# creating plotting area\n",
    "fig = plt.figure(figsize = (8,6))\n",
    "ax = plt.axes(xlim=(60, 600), ylim=(245, 600))\n",
    "line, = ax.plot([], [], c=\"b\",linewidth = 1.5, marker = 'o', markerfacecolor = 'r', markeredgecolor = 'r',markersize = 10)\n",
    "text = ax.text(450, 565, \"\", fontdict = font)\n",
    "\n",
    "# to plot only the final states of every simulated annealing iteration\n",
    "for iterations in range(100):\n",
    "    tsp_problem = TSP_problem(next_state)  \n",
    "    states.append(simulated_annealing(tsp_problem))\n",
    "    next_state = states[-1]\n",
    "    \n",
    "anim = animation.FuncAnimation(fig, animate, init_func=init,\n",
    "                           frames=len(states),interval=len(states), blit=True, repeat = False)\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",
  },
  "widgets": {
   "state": {},
   "version": "1.1.1"