{ "cells": [ { "cell_type": "markdown", "id": "a7d04ec8", "metadata": {}, "source": [ "# measuring the relaxation time with the idling gate\n", "In order to demonstrate the simulation of decoherence noise, we build an example that simulates a Ramsey experiment as a quantum circuit run on a noisy `Processor`.\n", "The Ramsey experiment consists of a qubit that is initialized in the excited state, undergoes a $\\pi/2$ rotation around the $x$ axis, idles for a time $t$, and is finally measured after another $\\pi/2$ rotation:" ] }, { "cell_type": "code", "execution_count": 1, "id": "644d4f06", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "T2: 20.0\n", "Fitted T2: 20.05671956607771\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAc8AAAEkCAYAAABJ+7pqAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAA9hAAAPYQGoP6dpAACP/UlEQVR4nO29d5xcZ33v/36mbe+72tWqN0uyLFsuWG7YBiPZyBAglwTi8IMANwkJSSg3IQkJgbRLLsmFFHKTQAgl4EAIxAkBXImxjW25ykW2el9trzM7feb5/XHOc+bMaMopI+1Ket6vl17amTnnmefszM53vu3zFVJKNBqNRqPROCew0BvQaDQajeZ8QxtPjUaj0Whcoo2nRqPRaDQu0cZTo9FoNBqXaOOp0Wg0Go1LtPHUaDQajcYl2nhqNBqNRuMSbTw1Go1Go3FJaKE3cL4ihBDAIBBd6L1oNBqNpq60AadlFRUhbTy9MwicWuhNaDQajeassBwYqvSgNp7eiQKcPHmS9vZ2TwtkMhnuv/9+du7cSTgcruvmFgsX+jVe6NcH+hovFPQ1OmNubo4VK1ZAjaiiNp4+aW9v92U8m5ubaW9vv6DfzBfyNV7o1wf6Gi8U9DXWF10wpNFoNBqNS7Tx1Gg0Go3GJdp4ajQajUbjEm08NRqNRqNxiTaeGo1Go9G4RBvPBaZKD65Go9FoFinaeC4giWSKZ/af5pGnXiCRTC30djQajUbjEN3nuYA8++I+5pMZ/uvBn/DAI09x7bYt3Lz9Cnq7Oxd6axqNRqOpgjaeC8i1V17Kvv37mE3kGZ2Y5tGn9vDYU3vYsnEtt1x3JetWLcOQ0NVoNBrNYkKHbReQSDjMYE8bH/3Fd/LL73orm9evQgIv7z/C333tu8zF5hd6ixqNRqMpg/Y8FwFCCDatW8WmdasYnZji0d17yGRzdLS1Wsc88+I+Nq5dQVtrywLuVKPRaDSgjeeio7+3m7ff+fqi+46fGuYb/34fAJvWr2LHa69lzYqlOqSr0Wg0C4Q2nucBew8ctX7ed+g4+w4dZ7C/lxtfczlXb91EQ+TCFHnWaDSaxYrOeZ4HJFMZ6+fWlibCoSCnRyf49n/9iE999h85PjSygLvTaDSaiw9tPM8DTo+OWz+nUhk+8eH38Zadr6W3q4NAIMDSJb3W41Mzc+Tz+YXYpkaj0Vw06LDtIkdKydBIwXhmslkSiRS3Xn8VN193JRNTM0TCIevYf/j6PWSyWa6/+jK2X7mFdl1gpNFoNHVHG89FTjKVJplKA9DV0cb0bJTZaIwlvV0EhGBJT5d17OTMHPPxBPOJJD/40RPc+/BuLtu4huuv2sol61YS0AVGGo1GUxe08VzkzEZjADQ1NtDb3WEaz/L9n71dHXzyo+9nz96DPP7MSxw7NcyLrx7mxVcP09XRxtvuuIWtm9ady+1rNBrNBcl5kfMUQnxQCHFMCJEUQuwWQlxb5diHhRCyzL/v2475SpnH7z03V+OOOdNQtre2WH2fs3OxiseHQyFec8VmPvT+n+Vjv/LzvPbabTQ3NjA9G6WpscE6LplKkdO5UY1Go/HEovc8hRDvAD4LfADYDXwYuE8IsVFKOVbmlJ8GIrbbPcALwLdLjrsXeK/t9qJUZp+LxQFob2suGM8KnmcpS5f08tNvvIU3veFGXjl4lHWrllmP/fC/n+SFVw6x/cpL2X7lFro72+u/eY1Go7lAWfTGE/go8EUp5ZcBhBAfAO4E3gf8WenBUsop+20hxDuBOGcaz5SUctH3eNg9z/Y2o/hHhXKdEgmH2HbpBuu2lJJXDx1jNhrj/kee4v5HnqK1pZm33XEzWzetIxw6H94WGo1Gs3As6k9JIUQEuBr4tLpPSpkXQjwIXO9wmfcD35RSlrprtwohxoBp4EfA70spJ6vspQFosN3VBpDJZMhkMuVPqoE6r9r5M3NRAFqam2htbrTu8/qcio+8/x28vP8I3//R48xG54nNx/nn79xLU2MDN1y9ldtvqRgZd4WTazyfudCvD/Q1Xijoa3S3Ri3EYh7GLIQYBIaAG6SUT9ju/wxwi5Rye43zr8UI9W6XUj5lu195o0eBdcD/BmLA9VLKXIW1PgV8svT+u+++m+bmZpdX5pxXjo0zNjPPusEu2psbeP7QCE2RENsvXV6X9Z985RTJdBaAYECQy0tW9LWzblk3YHip6WyOhvCi/p6l0Wg0dSEej3PXXXcBdEgp5yodd6F/Ir4feMluOAGklN+03XxJCPEicBi4FXiowlqfxsi9KtqAUzt37qS93Vu+MJPJ8MADD7Bjxw7C4fISe0N3/ydjM/Ncc/VVLF/ax/OHvokIBtm1a5en57STzmT48Z4vWrfXrBjkdTdcRU9Xh5UDPXDkBP/0re9zybqVvObyTWxev5pQKOj4OZxc4/nMhX59oK/xQkFfozPm5irayyIWu/GcAHJAf8n9/UDVfKUQogV4J/AHtZ5ESnlECDEBrKeC8ZRSprAVFSlR9nA47PuNWG2NhNnj2Wartk0kUgSDQQIBf8XSoxPT2OMOM3MxLr1kbdExJ06PkZfS0tRtbmrk6q0bec0Vm1m+dIljcfp6/J4WMxf69YG+xgsFfY21z3XCojaeUsq0EOJZ4DbgHgAhRMC8/fkap/8MRo7y67WeRwixHKMqd9jPfs8GiaRhr5ubGmhqMnKeEsOotpi3vTIybtRWdXe2MzUzx9RslHw+X2SU3/i667l66yaefuEVnn5hH7PRGI8+9QKPPvUCS3q7+PX3/gytzU2+9qHRaDTnG+dDn+dngV8UQrxHCLEZ+DugBVDVt18TQny6zHnvB+4pLQISQrQKIf5cCHGdEGK1EOI24D+AQ8B9Z/VKPBBPJAFobmwkFAxaE1TU/X6YnJ4FYMPq5QSDQfL5PDNlekiX9HZx52038gcffi+//PNv5cotlxAOBRGIIgN+4MgJYvGE731pNBrNYmdRe54AUspvCSH6gD8CBoA9wB1SylHzkJVAUbe/EGIjcBOws8ySOeBy4D1AJ3AauB/4hBmaXTTkpSRp8zyN/xtJpTN1MZ6qDaajvZXuzjbGJ2eYnJ6t2PMZCATYtH4Vm9avIplKMTUzZ4Vt05kM//St75PJZtm0fhVXb93IZRvXogUBNRrNhciiN54AUsrPUyFMK6W8tcx9+6H857aUMgHcXs/9nS2SyZSVk1Qh2+amRqZno/UxnrFCD2l3ZzvjkzNMz0YdndvY0MBgf591e2YuRl9PJ6eGx3jlwFFeOXCUSDhMV0crXc1Bcvk8F3aWRaPRXEycF8bzYkUZyEg4TChoVLgqDzQe9288o5Z6UQttLUa7TWzeW9h1SU8X/+uXfo7R8SmefWkfz760n6mZOUYnphkFHnlyDztvqdpZpNFoNOcN50PO86IlXhKyNX42PND5OnuerZbxjPtas7+vm12vv4EPve9niu7v6+6wft574Cj/fu+POXZqmMXcZ6zRaDSV0J7nIiZhGki7oHtzo2E8/YZtpZSWbm5bazOtLUbFbLROBT/Hh0aLbs8nCunk3c/v5aV9h3lk9x66OtrYdukGtm3ZwIrBfsetLxqNRrOQaOO5iIknlOdZqGhtNA1pyuz/9Lx2MkUuZ4gptbc228K2/jxPxdjkdNHtkfFC0fP1V11GJBzi5f1HmJ6N8t9PPMd/P/Ec3Z3tXLnlEnbddoOeParRaBY12nguYhIpw3jaPc+mhoj5mD/jqYxkU2MDoVDI8jy95jxLGZ8wjGdneyszczFGxwt6/Zs3rGbzhtWkM1n2HTrGnr0H2XvgCFMzcxw6dqrIcI5NTtPX3ak9Uo1Gs6jQxnMRk0oZAsUNDYUJa+rnpF/P0wz7qj5NJXRQb89z66Z1PPrUC4xOTJ9xTCQc4vLN67l883pS6QyvHjxGOFyQ/osnkvyf//d12lubrePWrFjqW1lJo9Fo/KKN5yImaXqejZGC8Ww0jafvsK3Kp5rFSKpgKDqfQErp29Mbn5oBYMOa5Tz61AvE5uPkcjmCwfK6uA2RMNu2bCi67/ToBKFgkJm5GI/s3sMju/fQ2tzEZRvXcvnm9WxYs5yQHp+m0WgWAP3Js4hJpZXnWeiQbKyb51mcT1Vh21wuRyqdprGhoeK5tcjlclb4d1l/H0KAlMZg766ONsfrrF+9nD/+rV/iwJETvPjqIfbuP0IsnuDJ5/fy5PN7Gejr5hd+9k76e7s971Wj0Wi8oI3nIkYZyMaGMz3PeoVtlfEMh0IEg0FyuRzxRMqX8YyahjMgBC0tTURCQVKZHLPRmCvjCUZo97KNa7ls41pyuRyHjg/x1PN7ee7lA4yMT/Gt/3yQ33jfz7L7+b28vP8IWzet49INqy1PWqPRaM4G2nguYpSBbLCHbc3iIRXS9YpdMxeMKTHNjQ1E5+OWGL1XVP9oW2szASFoCIdIZXKWHKBXgsEgG9eu5OTQKHAAMFpiUukMz718gANHTvDy/iMIYNWKpWy5ZA2XbVxLf2+3LjjSaDR1RRvPRUwqfTY9zzMFGJqa6mM8o5bxbAEgYhYBlROd98LhE0PWz/l8nhNDI7xl5028+OphXt5/hKGRcY6dHObYyWG+/9DjLOnt4mO/8i6CutBIo9HUCW08FzHVwrapVJq8lJ77IRPJ4rAtFFpifHue0YLsH0CDaTz9ep6K8ckZwPBso7E4J4ZGue2maxjs7+OOW69jZi7K3gNH2bv/CAePnqKzva3IcH7/oZ8w0NfDpg2rfY9102g0FyfaeC5iVEVtQxnjKYF0OlNkWN1QmvMEaDaNp1/1IuV5tpt5R+V5Ruf9G89sLsfUjDHpfcuGNTz5/F4mZ2aLjulsb+PGay7nxmsuJ5VOF/WuzkZjPPjYM4CRk12zcpBLN6xh84ZVDPT16PCuRqNxhDaeixjL87TlPMOhEIFAgHw+TzKV9mE8ywgwmIbUf85Tyf4Znmc4ZBjP+TqI2U9OzyKlpCESZvWKpTz5/F7LmJajIRIpyhkLIXjDTdew98BRhscmOXx8iMPHh/jeg4/R2d7KHbdex/Yrt/jep0ajubDRxnMRo1pV7AZSCEFjQ4R4ImkWDbV6WlsZyGL1ovqEbZXnqtpfwkEjZFoP46lCtr3dndbc0WrGs5T21hbuvO1G7rztRianZ9l74CivHjrG4WOnmJmLkUim+Ksv/Suz0RjvfMsOTg6NsnnDapYu0V6pRqMpoI3nIkVKWai2bSiehNkQCRNPJC3j6oVyxUjWuLOEP+M5n0iY66k2mGDR/X6YMeeN9nS109NlTGqZmol6yv/2dHVw8/Zt3Lx9G+lMlsPHT/HKgaMcOzUMwP0/3s3h40P810M/oaOtlU3rV7F5/SouWbuy6EuHRqO5+NDGc5GSyebI5/NAcasKQCRiGNO0D+OZLCP912S2rSR8t8EY56tinHp6nvYxah3trQSEIJcz2mA627154WD0k25ev5rvP/S4dd/YxDSbN6zm0NFTzEZj7H5+L7uf30sgEGDNiqX8/Nt2WrKGGo3m4kIbz0VKJlMwjMpYKhrCxm2vnmdeyoLnaVu7qdEUnfdZMDQfL+95JhJJ8vm8L21aVbHb1tpCMBCgrbWF2WiMuZg/46n2PTQybt2Ozse56607iYTDHDk+xKuHjvHqoWOMT85w8vQobS3N1jzSp194FSnhknUr6e3q0CFejeYC57wwnkKIDwK/BQwALwC/LqV8qsKxvwB8ueTulJSy0XaMAP4Q+EWgE/gJ8CtSyoN137xH0pksYAgDlPYnWp5nxpvxtHusRZW8lgBDfXtIQyFj/xJjFJofb63geRqVvG2tzcxGY0Rj/gXtRyeMyS9KBWl6Nsro+BTrVi1j0/pVbFq/irdxCxPTs4yMTRIKhawvOT9+co8lht/V0cbGtSu5ZN1KNqxZob1TjeYCZNEbTyHEO4DPAh8AdgMfBu4TQmyUUo5VOG0O2Gi7LUse/xjwG8B7gKPAH5trXiql9B9brAMZ03hGwme+RH7DtsrrDAhB2Casrqp6/eRSs6Y2LkCzaTQCZpFTMpVmPp7waTyLe0jbWuo3DUYVI/X1dBIIBAzjOWEYTzu9XR30mvlWMPLTV162gUPHhjh2cpjp2ailvyuALRvX8v53vtn3/jQazeJh0RtP4KPAF6WUXwYQQnwAuBN4H/BnFc6RUsqRcg+YXueHgT+RUv6Hed+7gVHgrcA367l5ryivspzx9Bu2tY86s4cXVWGSn4ktKuQrMGaPqoHbzU2NpvH0991EhW3bzTaYNtMDrYfnOTFl9Iv2dXcSCoXYd+h40RzSSgghuO3Ga7jj1utJpdMcOX6a/UdOcODICYbHJq3c78TUDP/0rf8ilUpz7ZVbWL96OauW9evJMBrNecii/qsVQkSAq4FPq/uklHkhxIPA9VVObRVCHAcCwHPAx6WUe83H1mCEfx+0rTkrhNhtrlnWeAohGgB7iWUbGLnJjMfwqTqv3PnKCNlDgwoVBk0mU56eO2aKFTREwkXnq/BwMpX2fE2zpgRfU1MDuVzOWkflU6Oxec9r5/N5YmY+takxQiaTsfKqs9GY53UVY5OFsG3E/IIyNjlddd3S1zAgBOtXL2P96mXc+frrmYvNk80av4ev/dsPGR6bBODeh58EjL7dVcsHWLdykHWrlrF8cAmhCmPbFopq79MLBX2NFwb1uEan5wpV8LAYEUIMAkPADVLKJ2z3fwa4RUq5vcw51wMbgBeBDuA3gZuBLVLKU0KIGzBynINSymHbef+K4bG+o8JePgV8svT+u+++m+bm+k/wmIomePHwKC2NYV6zqThsePDUJEMTUVYu6WDtYJfrtaejCV44PEpzY5hrbWsnUhl2vzpEQAhuvmKVp33PxJLsOTRCUyTE9kuXW/fvOTTCTCzJppW9DHR7K+xJZ3I8vvckADdfsYqAEJwcm+Xw6Wn6OpvZsnqJp3UVz+w/TSyR5rI1SxACXjoyVvb373fvAI2RELl8nkw2X3Tc6oFOVg90ApDPSxB4lmDUaDTuicfj3HXXXQAdUsqKTeSL2vP0gmlk7Yb2ceBV4JeBT/hY+tMYuVdFG3Bq586dtLe3e1owk8nwwAMPsGPHDsLh4oravQeO8uLhH9LT082uXbuKHvvhw08yNPEcK1auZNfO17p+3r0HjvLC4R/SW7L2fDzB7le/TF5Kbr/jDk9C6q8cPMaeQz+w1lbXuGzpADMHj7Fp86XccPVlrtcFo6Dn8b3fpKmxgTfdeScAz+89wOH/eJC29s4zfk9ueeagUWd22+tuIRAQvHTkX0EEq65b7TW089K+w0XGs7eni1//hbczNjHN4ROnOXx8iCPHh7jjDbeyduWgdc63vvcQq1csZd3KZaxZuZTlA0sIhc6tZ+r0Gs9n9DVeGNTjGufmnImuLHbjOQHkgP6S+/uBsjnNUqSUGSHE88B68y51Xj8wbDu0H9hTZZ0UYDVAqlxhOBz2/UYst0Y+b0QEGiJnPqaUgLK5vKfnzuYMb6exoaHo/NaWgocjJZ7WzmSNHGdTU/HaSvovk8l6/n2pCuSWpkZrjc52ozJ2PpH09Trk8nnmTQ3cnq4OAkEl7JCEksKqctR6H5wcNlpgNq1bxb7DxxkenSQQCLB8sJ/lg/3cct2VSCmRFDzNk8PjpDNZDhw5yYEjhuENh4KsXDbA2pWD3Piay+lo89ee44Z6vNcXO/oaLwz8XKPT8xb1jCYpZRp4FrhN3SeECJi3n6h0nh0hRBDYSsFQHsUwoPY124HtTtc8F6iCoXIvpO9qW2taS/HaoVCoKO/phaQp7Vc6TFspGfkRYFDFRs3NBTH7FrNy17+YfdwyXC0tzTQ3NlgGc7YOo9TGzDaYyzatpbmpkVw+z3BJMZIQoihE++YdN/Gbv3wXb739ZrZuWkdLcxOZbI7Dx4d44NGnyeUKId99h47z4GNP123sm0ajqc5i9zzBCJV+VQjxDPAURqVsC2YvpxDia8CQlPJ3zdt/ADwJHMLo4fwtYBXwj2AkNYUQfwn8vhDiIIVWldPAPefommqSrtKq0mAaz5THpHi5IdvW2qZurmo3cYsyjk0lgvXWHNKk90peJb7QYmt1sSQF40mklJ7FCUoHeAN0tLcwMTXLbHSe3u5Oz/sGmJguVPIu6e3i2MlhJqZmWLG0cp42IATLBvpYNtBneabjkzMcOTHE6dEJqx8V4O7/uJ9oLM73H3qc3q4O1q4aZM2KQdauWkZfd6cWbdBo6syiN55Sym8JIfqAP8Kokt0D3CGlHDUPWQnYqy66gC+ax05jeK43SClfsR3zGQwD/AUMA/uYueai6PGEc9PnWW4ii6Wbm/K2dsI0jo2N5T3PpB/P0/Qu7TM4m5sMQ5rL50lnMmW/EDjBaoEx+0fBCAlPTM0yMxf1umXAqBKenDbyKL3dnfR0dXDs5DBT084F7cHwTJf0drGkt7hI7MiJ00WtOhPTs0xMz/LUnlcBo3r49z/0XutLgV+VJ41Gcx4YTwAp5eeBz1d47NaS2x8BPlJjPQn8gflvUZLJGsazXK5NtVF4NZ6Wrm0ZQ2MZOY+eZyFsW8Hz9NFDGldhW5vxjIRDBINBcrkc8UTSv/FsLRjPNnMeqX0eqBdm5mLkcjmCgQCd7a30mNNgSueQemX/4eNFt2+4eivNTY0cPXma46dG6DQ1gBWf/aLRjbV6+VJWLR9g9YqlWlJQo3HJeWE8L0ZU2LZcztNv2FZ5nqXTWoy1TZUhrzlP87zSqSOFnKePsK3yPG1hWyEELU2NzMXmmY8n6erwVvlsF5xXtFrqRf6M56QZsu3uaicQCFjTYCZdep6VOHbKqIHr7mxnamaOudg8P/Om1wOQzWaJ2tSXkqk0p0fGkcDQyDg/eeZFAFqbm1i1fIAtl6zleo/V0BrNxYSO3SxSVKPu2QjblhuyrbBUhjyurXKeFT1PH7NC42XCtlDwRP0UDal8amtLoWdXyQj6lf5TRTwqR1mYQ+rf88xLyYkhw3i+7vqrAIrE7UOhUNEXisaGCJ/48Pt4z9t3cct1V7Jq+QDBYJBYPMHeA0c5cmKosHY+z3d/+DC7n9/LXDxpTfnRaDTa81y0FDzP+odtC55nuZxnpOgYt6iCoDM9T/+i88o4NpcYzxaz+nbel/FUXm1hbWVIlaqRV2ajhvFsN9tKlBGdnYv5KnIy1oiSTKUJBgJccel6vvPDh5mejZJKZ6wIRSldHW10dbSxbcsGwPBOT42Mc+zkMANLeqzjRsYnefSpF6zbL/z5F1m1rJ+Vg/2sWNbP2pWDVquQRnOx4ch4CiE+W/soAynlR71vR6MoFAxVDtt6naqiioEqFQzZj3FLLc8z4cPzVOc2NRUb5np4nspA2kXrrbCtT+Op8qkdZkhYFSVlsjkSydQZXwbcoPR4e7o6aGttoaW5ifl4grHJ6aqVvHZCoRCrly9l9fKlRfdHIhHWrFzK0RNGl1cul+PIidMcOXEagDtuvY7bbzFEvuYTSY4cH2Llsv5z2nuq0SwUTj3PKx0et3i1/s4zqgnDq7CtGpjttnLS8jzLeCb18jzPqLZtLORSvXpblvEsWbu5sdCu4hVrBmmR51mfsO1sSSVvOBSiqbGBRDLFXGzel/Ecn5oBoLfbyKP293Zx5ESCsYkpx8azEr1dHcTjxV92dtx8LYlEkhOnR1m9fMC6//CxU3z5X78PQEdbCysG+1kxaHipy5b2WcVXGs2FgiPjKaV83dneiKYYJ2Fb47jMGYIEtbD6PMt4nko4IenR80zW6POUGPnUcl5vLSoZT6Ve5MerLYRtbZ5nc32qbS3P09YG097aYhjP6DwDfT2VTq2J8jxVH2pfTxdHTpy2xqv5IZFMWTNOO1oamJ1P0dTQwK7XnTmTQUrJ0iU9jIxPMRudZ3b/EV7ef8R6/L0/eyeXbzZEvmLzcTLZLJ3tbbrCV3PeonOei5RqYdtwKIhAGaKsa+NZUBiqb84zL6WVhy1dOxwKERCCvJQkUynXxjMvpbXvM4ynz0peKWWhYKhM2DaZSpPNZj2PDlM5T3s4s72tmdGJKWs+qVcmTM9TVfCqfGo9lIZOnh611uxsCjI7n+LUSPkRuldcuoErLt1AKp1haGSME0OjnDw9ysnTY4xPzRR9QXhqz6t878HHaG5qpK2lmY72FrZv28LypUvo7enUQvia8wJPnwZCiGuAn8UQKCj6FJRS/nQd9nXRU5DnO/MlEkIQiURIpdOu855SSpLpyn2eyhv1YjwzmawVt4+UrC2EoLGxgXgiSSKZptNlR0kqlbbWPsOrbVTFSN48z3QmQ9acO2ovGGpsbLC+pMSTKdpb3f+5SCnPGOANhZYY5ZV6ZXrWEHBQFbzKeKr7/XDytGEoVyxdgkwb650arjR/3qAhEmbtymWsXVmYRJNMpa1UAxgh8kAgQDyRJJ5IMjoxZWn3NkTCfOQX30l/bzdg5LEj4ZCeeapZdLh+Rwoh3gl8DbgP2AncD1yCIaz+73Xd3UVMNYUhgEgkZBhPlxW32VzOajko1bYFfwVDaZvBLWf0G5X0nwcPUYVkw6EzP0j9Sv+pkG0oGCzy9ANC0NTYQDyZIpFIFfWAOiWVzlgDwe2GWa2lvFKvWF5tu+HVdpr/+1VFgsJ804ElPUQnjd//+MQ0mWy2plC+ndIow5t33MQbbn4Nf/S5f7JSCGq9TDZnfREA+K8Hf8LuPa/Q39vF0v5eBpf0Mtjfy9L+XjraWnTYV7NgePk693HgI1LKvxVCRIEPYejD/gPFU0o0PrAUhioo/DeEw0Rx349pN1yl3iH4UxhSe4mEw2VDbyrc6kUcvlK+0++6YCsWamo848O4qamReDJFPOmtGElVAJca5rbW5qLn9kIulyNmerWdpldrhW1no77bYMZtlbyp2VEaImFS6QxT03P093V7XhfgxKkRo8XGVIfKZLN86iPvI5PNFRnm8akZ8vk8w2OTDI9N8hz7rceamxr55Efeb33BnJ6N0tzUWLFFR6OpJ16M5zrg++bPaaDFFFv/HPAjygyM1rinmjA82IQSXIZtC6Lw5Q1cwfP0bjwrfXj58RCrGU+/nmclVST7ffGEP8Pc0txUZMhUha2f3tQ5cxJMMBCgxaxm7TD7LtOZLPFk6gxBCTdM2Cp5h48Lero6OD06wcT0jG/jeXzIyKdecel6hkcnGB6b5OjJEav3VPGr7/5pZuZinB6dYHh0gtNjE5wenWB8YprGhkjR38fd99zP4WOn6O3uZGBJD1JKNq9fxSXrVtHd2a5zqZq64sV4TmMMggYYAi4DXsIQWNf16HVCKQxVCo95VRmyvMMKBq6Q8/QQts04M5719jz95jzV2qqdxk6zz0recuILxu0m83HvnmdBfKHFMgyRcIiWpkbmE0lm52KejWcylbbE5nvNYqSeznbDeE75V0ZSxUgrB/tpbmw0jefpM4ynEMISddhyyRrr/kw2e0a+ODafQGJ4q6qFR1X89nZ18Hu/8QvWscdODdPa3ER3V4c2qhpPeDGejwA7MAzmt4G/EkK83rzvoTru7aIlL6U1VLqi5xn2JqNXyzv0o21byzBbRs6X53mmgbOqbf16nmWqlq2RZx49xPmE8jwrqCL56E1Vc0ZLRQlaW5uZTyR99acqPd6WpkbrC4uq6K2P8TSLkQaXWO9F1RbjhHAoZO1H8du/+i6i83H27D3Id3/4sHV/QAi6SirUvvrtHzAzFyMcCtHf182Snk6mJ2fYe+AoKwb7z1hboynFi/H8NUB9EvwpkAFuAL4D/Emd9nVRo4qFoErO02PYtqCZW8k7NPs8vXieNQxzk4/JKokKQ7ahYJTTmQy5fN4a6O2Uan2vVj7Vp+epRqcpCp6nD+NZZowaGNNgRsenfPWnqoIju9GplyZvMpWyvObB/l4rnD067tx4VqKtpZmpmWLB/d6eTt799jdat7O5HC3NTcTmE2SyWU4Nj1lVxEf/7YesWTnIb7z3Z6zjH3vqBTraWxno66a7q8P1+0tzYeLaeEopp2w/54E/q+uONJaBg/JVq+A9bJuu0j9qrGsYkFwuRzaXIxQMOl47VaUFBuz6tu4NkbV2OQNnuy+ZSrsOVSqx+tIWGIDmRn/Sf4X+0RLP09xjKp12/XtWKONYqt6j+lOjPjxP1SeqqnehUNE767O9Rgk4tLY00djQYPWAzszFSCRTZUPzbjh41Gh7+dk33ca//eC/GZuYLvpCGgoG+c1fvsucszrLyPgUp0fGee7FlwlFmlg+UFBmSmcyfPeHD1ttUsFAwBiYHgiwce1Ktm5ax+YNq33tV3N+4rXPMwCsB5ZQMplFSvlIHfZ1UWOpC4WCFfMxKpzrNmybtkKr5V/6RpvXmEqlCTU3lT2uHKo3tHLY1ruYgSXsUGbtYDBIOBQik82S9FAko/ZTKikIBR3denueTY0NCCFMgYZkkfqQU1RYVlXuKgpzSL0bz3Ih4c52Y49+BRhUPrLPVEVqamygo62F2eg8YxPTrLLJ/rklnckyPDYJwKb1qxjo6+b06AQnT49alciKQCBAX08XfT1dbFq3kuTsKLt27SqK9qTSGa7aupGR8SnGJqbIZHPW9T/x3Muks1nLeGazWf7+6/fQ193Jkl5j3f7eLu2tXqB46fO8DrgbWAWUfrJLwP1XaE0RmSqzPBXKu3NdMFQjbBsMBgkFg2RzOdKZDC24MZ41CoZ85VMrh1bBMMyZWNZTSLiQ86y/52mNUSvxPAOBAM2NDcwnksQTCW/G0xqjVvwaqWkwUV9h2zKep2lI5+MJ172ediZMz7Ovp8u6r6erg9noPJMzs76M5/DoBPl8ntaWJjrbW1mxdAmnRyc4NTxmyQO6oa2lmXf99B2AUYvwlW9/n5dePWw9bq9JmJie5fDxIQ4fHypaIxgI0NPdwQ1Xb+WW6wyZ8Hw+z3w8QWtLs+5VPU/x8u7/e+AZ4E6Mvk4tBl9nqonCK9RjaVs4ytHaNYp61GPZRM6zV1uxGMmHelE1VSQwin2isbgnD1GFkcuGhFXBkEfPM2ZrVSmlpbmJ+UTSc95TVcNWCtv68TyV8eywGc+mxgbLw5+di1l6um4Zn1Z6vIWinO6uDo6cOM2UzwHhKne5YukShBAsX7qE3XteqamM5ASZz3P46CnA8Gr3HTrOjE3Jqb21hXf99O2MTUwzNjnN2MQ045MzZLJZxiami/6eJqZn+fTnv0ZDJEJvdwe93Z30dXfQ3NRER3sLG9euLPue0SwevBjPDcDbpZSH6r2ZSgghPgj8FjAAvAD8upTyqQrH/iLwbowWGoBngY/bjxdCfAV4T8mp90kp76jz1j1RTRReoTzHTNZtwZCZ86ziNTREwoYSkNc2mBpFTl7aYFJWUU/1kLAnz7PCDFIoTGxJeOzzrOR5gjnBZdJ7u4oyjq0tlcK2PtpgynieQgg621sZn5phNjrv2XhOmwU9diWhbjOkWlrs45axyWkA+s086tL+XgBGJ6Z9rQtGhbDqnX3jrdex79Bxjp4ctiYbNTc1cvXWTUXn5KVkZjbK2OQ0PbbrnZmNIjC+SA6NjBcNMAfYvGE1v3TXW4xj52I8unsPvd2dlqHtaG/VLTYLjBfjuRsj33lOjKcQ4h3AZ4EPmM/9YeA+IcRGKWW5r5O3Av8CPA4kgd8G7hdCbJFS2uMp9wLvtd32PpKjzlQThVeEvXqeGWeeJ3jvIa3cBuNd+i9p5TyrFyN58TwrzSAF+8QWfwVD5fKwSoTeq1BC1CoYKvZQ2qywrXfPUxnmUknCDtN4+pH/U16tPQepWkMm62Q8l5gh4SW9xv/TM3O+Qs0AJ8ze1FXLB1i+dAmNDRGSKcP4rRjsL3tOQAi6O9uLvigAXLJ2JZ/5vQ8yOTPHxNQM45MzPPjY01YU4tDRk8QTSZqbGhkem+BHjz9bdH44FKSnyzCmN77mcjatWwUo+U1ZNWqlqQ9efsN/A/xfIcQARq9n0SehlPLFemzMxkeBL0opvwwghPgARsj4fZSp9JVS/rz9thDifwL/A7gNQ5NXkZJSjtR5r3WhIM1X/7BtKl3bMDd47CGtZZj9hG0LOc8KbTCN3vOp1RSGmn0rDJ056sxaW6kMefA805ms9Tsp9Tz9hm1z+bwVpi71mNvN4qSox2kwyhMD6GwvGM96tcGMW/nUTsD4gqKM3MTUDEuX9Hpe+6QKCQ/2EwgEWLV8gP2HT3ByeKyi8axGKBSiv7eb/t5uRiem+I/7H7Vy4bF4gr0HjvKaKzbT0dbCa6+9gompGSamZpmcmSOTzTEyPsnI+GRRLvfQsVP8w9fvoa21mZ7ODrq72ukxjXdDJEzcYy+05ky8GM/vmP//k+0+CdYAiroVDAkhIsDVwKetJ5IyL4R4EDhzqGB5moEwUNpEdqsQYgxDMelHwO9LKSer7KUBsH+6toHRVpJx2WupUOeVnq+8nHAwWHFtNQA7lUq7en6V3wsFAxXPU0Y7kUi6Wzup1g6ecW2ZTIagGWZKpd3t2Vjb+KMPBsrvW30ZmI8nPKxt7DtY5ncSDhlv50w2SyKRJBQqfntXeg3B+DKhvghFwqEzjlEGPxqLu97zrOn5BQMBggFRdL7yoNOZLPPz8apRhnIooyuAUKj4tVQGfzYa8/S+j8bi5PJ5hBA0N0asNVQrz8xcjHQ67amIJpvNWWHf7o42a+2+7k5ODo8xPDphqSWVUu11VJw0JQWXLukhk8nQ39vN/sMnGBoe9/wZoHjxFSOQt2H1cpYv7eOhnzzLS/sOs+3S9fR1d/LmN9xoHZvL55mZjTIxPcvE1CwrlvZZz6+KsaKxONFYnGOniuXGA0Jw/cgYywaWcGp4jCef30t3RztdncrItp0hJXk+4eR1dLpGLYSU7up9hBCrqj0upTzuasHqzzWIIQF4g5TyCdv9nwFukVJud7DG/wNuB7ZIKZPmfe8E4hiC9uuA/w3EgOullLkK63yKMrq9d999N83N9VUlPD0Z5cDJSXram9i6tvw32onZOC8fHaOtOcLVlww6XvuVY+OMzcyzbrCLFUvKf5C8dGSUybkEl6zoYbCnrewx5dhzaISZWJLNq3rp72o94/FkOsuTr5xCCLjlitWO1wV4fO9J0pkcV1+ylLbmMz3EQ0NTnBqfY8WSdtYNutNdffTF4+Tykms3L6O5xLOVUvLjF4y39PVbltPgIhxmv96bL191xgfSidFZjgxP09/VwuZVfa72HEukeWb/aSKhIDdctuKMPT/64gnyUrJ98zKaKnjrlZhPpnl632lCwQA3bV1Z9Njx0RmODs8w0N3KppXuvbi5+RTPHRwmEg5yw5bCvnP5PI++eAKAGy9bYX1pcUM8leGpV4cIBASv3brS+n2/cnycsel51i7tYmW/N+WgvJQ8+uJxpITrLl1OYyTEyFSMfScm6Ghp4MoNSz2tq3jx8ChT0QTrl3XT1hzh+YMjhIIBbrxshStDJqUkm8uTSGdJprMkU1liyTRj04Xe3M7WRratH+D0RJQDp870FwJC0BAJsmF5D91tRhQjlc4ST2VojIRoCIcIBM5P4+qEeDzOXXfdBdAhpayYR/AiklA343i2EUL8DvBO4FZlOAGklN+0HfaSEOJF4DBGvrSSxOCnMXKvijbg1M6dO2lvdzmc0iSTyfDAAw+wY8eOoraUx55+kQMnH2P5smXs2rWz7LkHj57k5aPfo6WllV27djl+zrFv/4CxmXm2XXEF26+8tOwx0/fcz+Qrh9i0cRM3XXuF47UPf+XfmIkl2X7ttVxq9r7ZrzGby/PkK19CSti58/YzvLhqPPHKF4Ecr3/966z+QDsPPPo0p8afZnDZcna98VbH6+bzeR7e8/cA3LFzxxkhUICn9n+JRDLFDTfeZM2ZVFR6DQGGRsZ58pVv09rczJ133nnmunte4cjww3T19Lp6DQEOHx/imf3/QWdHe9lzXzz2z0zPRrnm2u2sWuau9ePIidM8ve8ea237NT6/9yBHhx+ms6vb9Z4BXtp3mOcODjOw5MxrfuaA8Xu+7vobPQnPHzp2iqdeHaKns6Po9x14+El+9PhzDAwuY9ftN5c9t9rrCEY4+JEXjhMOh3jrW95MQAiGRsbZ90/fJpOHN77xjZ69NSklT/3llwF4447XsbS/lxf/4otkc3muu+Em31KBDz32DPc98hTtrS3MxeaZiSW59robSaZSrDx4jOnZKFMzc0zNzDEXnScvJYlUluuvu451q4y5rE8+v5fv/vDHgBGRaGttoaujlc72Nlpbmmhva+aarZtpa20mL+WCFTPVeh2dMDfnLO/upc/zpyo8JDEKdA5JKY+6XbcCE0AOY1aonX6gar5SCPGbwO8Ab6iVh5VSHhFCTGAUQpU1nlLKFLaiIvWHEg6HPb9IitI1cua8zYaGSMW1VSFLNptz9fxKM7epqaHiear4JpvLu1pb5V+bmxrPOC8cDtPQWDCWeSkdry2ltIqXWluay56ncorpTNbVnu0FRm2tLWWHLjc3NZJIpoxxWRXWLvc+SJm/j9aWprLntZt9k4lEyvV7SP2uW5rP/F0b19LM9GyUZCrjeu2U9bsu3nc4HKbTLPKZTyQ9ve+VIEVHW+sZ57e3tZBIpphPuv99QKGAqquzrej8btP4zMXiNdet9Pc8ZeZpl/R0We1Sg6YSUTyRIpPLexbhn5yeJZ5IEgwEWLlsgFAoxLKBPk4MjXJ6bJIBH3lagH2HDY9+x82v4UePPsl0NMkrh45x243XsMY2tBwMoYeZuRjTs1GWL11i/S4i4TB93Z1Mz0bJ5nLMxeaZi81b03EAHtn9Ih/9xXey7/BxvvfAY3R1tNHZ0UZXexud7a10tLcavbeD/b5VpGrh53PZ6Xlecp73UMhx2rHynkKIx4C3Sil91YdLKdNCiGcxin3uAUvd6Dbg85XOE0J8DPg94HYp5TO1nkcIsRzoYZHMI601CNt4TGnbeqy2rSrA4E03t1afZzAQIBwKkskaPaRO+9jSmYzVTFyxz9MSnXdX2JOw5WnLGU772m6LhuKWulD5D1VLHD7hvmBItcBUWrvVh8qQJb5Q5vWxBBg8Fgyp9plya3e0tTA6PnXGtBSnFNprilMN6vb0rPcKYSVar6p3wfj7VN7c5PSsZ+OpCpGW9vda78GVg/2cGBrl5Okxrrpso+d9z8cTnBgy/IxN61axZ88LTEeTHD52ittuvOaM40OhkNkS01l0//Yrt7D9yi1IKYnFE0zPzDE9G+WJZ19m/5ET1nPd/8hTtLU0k0ylrRmspfzaL7zd8mj37D3Ik8+/TGebYVyVgTWMbRvNphLXYsSL8dyBIQj/e4DqnbwW+GMMYfhZjMHYfwG8vw57/CzwVSHEM+bzfRhoAVT17deAISnl75q3fxv4I+Au4JhZFQwQk1LGhBCtGLnL72B4r+uAz2C03txXh/36Ju1AYajQquJRns+BYfY+saW8gVOPZbIJV/2YSbO1RQhRcd+FcWfuqgmTVdpUFM0eJfqqjVEDaGnyLg6v2luaKnxgt/kwctXaa9ps6kVehm1X6k2FQluMV+M5PaeqeIvz7V0dxu0ZH8ZTTZIpTRl0d7UzF5tnanqOlR4qbqEgiD/Y32PdZ/WnjlesYXTE8aERJIbH3NHWQmer8ZoeOXGaXC5H0KWmshCCtpZm2lqaWT7Yz3fM6TU3b9/GI7v38MyL+/jkh9/HNVdsYno2yvRMlOnZKLPRGDNzMWajsaLX5/ToOPtNz7gcG9etIp5I8IabXkNjQ4T9R07Q0dZqebIdbS20tTS7vo564MV4/hXwS1LKx233PSSESAJfkFJuEUJ8mOJqXM9IKb8lhOjDMIgDwB7gDimlihesBPK2U34FiAD/VrLUHwKfwggDX44hktAJnAbuBz5hhmYXnMLkk/p7npZXW6UCs+EszQq1HosnXK1ttalEwhU/rJs8jjtTY8zKzfIsrO1Noq+W8bQbZdVo73ht0wturmiYvQ/brtZeo4xeLpcjmUq7Dr+p0GqppCAUpsPMxbwZz3L9o8ZtoyZhPpEknclUjbpUQrXQlOYfe7o6OHZy2Brh5gXl1fb3FozngJnzHfE5aeaEGVZdaUoetjZFrNad0YkpBvvdFarZOTU8RjQWpyES5s07buKVg0eZmJrlyIkhrrh0wxn1AeW48rJL6O3uNAzrXIyZqPH/7FyMWDzB/sNGic0/f+dert22mcefffmMNQTG+/KXfv4t9NsiA2cbL8ZzHVAuozoHrDV/Pgj4C9TbkFJ+ngphWinlrSW3V9dYK4FRfbtocaYwZDyWz+ddfYOspW0LBeOXcuHV5nI5cjkjn1opbAsFDy/potfTEoWv4h1a67qc2GKJL5QZdabwOkqtlvG0e40Jl9NgLOWiCuc0N3vX5J2vIikYCYdoiERIpdNE5+OujaclZl/F8/Q6tUV5lh0lnmdjQ4SGSJhUOsPMXMwSUHCDkg0sFTtQqkG+jKdpIPv7CvtShmd6Nkoqna4azamGCtmuMr1iIQSD/b0cOXGaU8PjvoznqwePAYZ3GAoG2XLJWn785PO8evAYV1y6ofrJJkuX9JbtvY3Nx/nDz32JbM7wi7K5HBPTs7z22m3MmgZ2Zi5GNGYUOEXn42UHO5xNvBjPZ4E/F0K8W0o5DmB6hp8BnjaP2QCcrM8WLz5qjQ0zHiu8dOlMliaHxjOddu55ugnb2j1gR2u7MEROwsGNHg2cMrblhmxba/vMp1YyMKFgkEg4TDqTIZFIujOeyephW7WWP+NZKSTcRCqdJjYfd22IYlU8TyU879XzrJTzNGQF2xidmGJ6Nup6z/l83sqX9nSVGE+fykh5KRlXkoI2T62luYnWFmPm6OjEtOeQ8OnRCQCWLy2MWbMbz2u3eVoWMEK/ABvWLDf/X8GPn3yeoyf9l47seeUg2Vye5UuX8Nbbb+bzX/k3jp8a4X3veHPRl3MlsD8bnaervZV8Pl9l1friZU7O+4E1wCkhxCEhxCHgFLAa+J/mMa3owdiecRK2DQaDVgjTaehWSulMdN5D2FaFVoOBQNXZlMoAujHMyRq6tlAwnqlUGje9yyrM68TzdJtPTaSqG08ohG7dCs+r4iV1/pnr+jGelcO2AK0+VIaq5jzbvOc8U+mM9TtUOU47KpTrJe85G50nZ4bV7SPaoGA8pzx6ntFYnEw2hxCiaPA4FPKrk1Pe1k4kU5YXP7CkYJiXDRie3unR8bLnOSEvpeXVrl6+1PzfCA2PTU5bRWde2XvAaNjYdukG1q4cpKujjVQ6w6FjxT5ZIBCgrbWF5UuXnPO8p2vjKaXcD1wKvAX4a/PfT2GIEBwwj7lHSvnP9dzoxYQTeT578YzToqFcLmcZluperXvPs5aurUIZQHdrq5xnlUIk08BJ3BVRVdO1VViep8uQcC3P0/5YwqWRq1Vt2+wj51mt2ha8C8/n8nlrP6V6vGA3njFXX4AAS2u3sSFS9ouQKlLxMot0enbOWqM0L23JCs5GPXk9Vqi5reWMmZ+WV+vRMKtK18721qLfifJwVa7VC2MTUyRTacKhkFXc1NLcZFUjHz/lXfk0k81yyJxec+klaxBCWNq9B44snoCmpwmtUsq8lPJeKeVfm//uk1KeO3/5AsdJ2Nb+uFM5qZTL0KobI+SkWMhY270GrRKSr2bgwqGQ1ZjtJnSbcpLzbCwU9rjBifEseIhuPU9nxtOf51mpDcYwfG6F51U4WAhRdt8q55nJ5lyH38vp5drp9ON5zhneW2khEphGLxgkn897NMzRimv7DQmPmJW6A0t6iu5XbSix+YTnIe9qCsyygb4io6/Cw8NjE57WBWN6TTaXo7WlySqcumSdoXR18OjiMZ6Ocp5CiN/AqKRNmj9XREr513XZ2UVMxkHBkP1xp2FbFYYNBoNVJ9t7CdumHeQljcfde55W2LaKYRZC0NgQIZ5MkUylKfNZVJaENY7MQT7VdSWvc88z7nJqS7xGta3dKLtpKclmC4LzlcK2XttgVMi2pbmxbGVxJByiqbHBDDfGXBUjlRvebUfdP+1hGsxs1Fi7vczA8kAgQHdnG+OTM0zNzJ1RUFQL5TFXNZ4ePU/VXlOa421siFj9qeOT06x0qUAFhVzqYH9xsc9ScxTcSJn+TqccM3Oma1YMWu/bNSsMCdKR8SlS6UzNCNe5wGnB0EeAb2AoCH2kynESI4yr8YGTvKTxuPI8HRpPh+t6aVWxt5NUX9v9ZBUnYVswQrfKeDrFSZ+n5XmehbCtMnJu5oXmcjnrd1LL88zn86TSmarXZ0eFVQNCVKxetMK2cXfGszBCrbIWtFIZmo3OM9DXU/G4UlRur7TSVlHwPN17h9barWcaT4Cezg7GJ2eYnJ5l/erlrtZWnmdnWePpr5JXtdeUM+h9PZ3MxeYZm5zxZDxVSLjUeCovt5w4glOUmL3KpYLh4SuDf3p03DKmC4mjsK2Uco2aOGL+XOnf2lpraWrjZJ6n8bi7nKfbcLCXnGfNsK2V83Rj4GoXDIG3ituEi7CtW88z6cR4NrovGLKH2iqtHQmHLHF1N6FbFVptbm6sqE/aahNKcEOhWKiyslSbR2WkqFmhWzp/VNHpo5J3zvQ8Kxlmdb+XFhsrbFsm3Kw8z5nZqNUG5oZJs72mnDau8kZVpa9bhk3Pc2l/8RecpabxHJ2YtmRG3aJCwisGlxTdv9y8ffJ0uTHO5x5POU87QoigEGKbEOLcdade4Fh9njUG91ph27S7sG0k4szzzOXzZB3+0daS5ius7b7aVh1baRC2wl5x63htZTwdhG3d5Icy2WxBR7ha2NbyPJ0bOBWybWyIVBVW8FI0pPKdrVUmBSnjN+/aeKo2lcprt3ms5J2LqeHd5ddWIddEMuVaWEQZxXJhWzC8IsDTgPBqnmdbawuhYJC8lEx7yafOlO9NBegzjeeYB+OZzmStMPmSEiGErs52IuEwuVyOiakZ12snkilrrNyygeIe1OXm7aGR89R4CiH+UgjxfvPnIPAI8BxwUghxa113d5HiRCQBbLlJx55nbYEE+7rgPHTruNrW6vP0EBKuEXr05Hkq77CK56nCl5ls1rEHoLxOQfV9e/E8axULWWt7KBqK1ejxtD/mdoi35XlWMcwF+T9vnmdrBePZ2BCxvoy69T5VzrO0TUXhx/OcqVIwFBCCbjN067YVJplKW1+ayhnPJb2dQGF4uBtUGLmxIXJGb3JACKstxkvoVnmdXR1tZ7y/1cDx89nzfDvwgvnzmzH6OzcBn8PQvNX4IGcqBoGTnKfxuOOcpwNdWzCa91VBkVMPseDVOjRwHhSGarfBuDeeSRfqReC817Mg+9dQdTxTk+rzdON5Js+e8aymLqRQLSzxRNJVe0Yh51llbfMxt20w0RqepxCiqBXGKVJKq9q2o4LnqYqRZl16h+lM1vqyUs54gr1oyF3Frcp3tjQ1ln1v99nCtnmXbUHKo+zt7ixbiOanaEgZ3FKvE2D5UuO+0fEp19GDs4EX49lLYRzYLuDbZn/nPwFb67WxixW7IXQyPglcVNs60LVVuG1XSTk0zF4qeZPKq63leXpog1EeYrWwbTAQsLx1pypDToqFAJpN3Vw3IeGCQIJD4+lCeL5Wm4p9XYm7kLBVbesg5+na81SyfxVynmDvI3W+djKVtvqu2yt5nm3K83RnPFWYtyESrvg+6ek0jOfE9IyrtZWx7e4qX/3b09lOIBAgncm6+jIBlUXyFX6KhsbU9JoyKlAdba20tjSRl9JXK0y98GI8R4FLzZDtHcAD5v3NGKLrGh/YezbDNYZFR0JuC4achW3Bpm/rOmzrzMB56sV06tU6XFtK6UjbFuwVt049T4fG04N3WAjb1n/t+RoCCWC0Oqlws5vQbcxBta2XnGcqnbbef5U8Tyh4jm6MnDq2ubGh4hdDFbaNzSfIZp17RDO2fGelViIVcnU7Tm2qSr4TjNew2/R23Xq1k6YhrzSkW4kweClGUjnYJWUE3oUQhYIkn4L59cCL8fwy8K/AyxhfPh80798O7KvTvi5aVJFJJByq2ZtXyHk6DdvWnhN6xtpOw7YZp6FVH+pFNapt3YZt05msFbKq1ucJBc+03p5nk4dxZ5bxbKzueXrRt405CNtCwXt0E16NOqi29TKHVBnaSDhc9cubNfLMRc5ThWIreZ1g/J6VJOWcC6NfTSBB0dVpziKd8Wo8yxs4KAwJd9sKM256nr3d5dfu6+kEjPCu25Dw2ERl4wk2daTz0XhKKT+FoWH7BeBG2xivHPBn9dvaxYnyDp1MMy/kPM+C5+myXcV5wZDx4Zb2lPOsbz5V9XgaUofV960KipwaOedhWyX9l3Zc2u86bOtCgMFJztP+uBfPs2q1bcm8UCcog9VWxesEb9q5tSptwXj/dHiQ/5uuoYoEWN6hkgh0ijKePVVEG7z2kdpznuXo6mwnGAiQyeaYdVGBrKbeQPmwLRSMp5cq4XrjZaoKUsrSWZlIKb/qfzuaQi9m7ZfGq8LQWc15Oqy2zWRz5PL5qkpHhbVrF/XYH3ea8yyIwkdqevmW5+k6bFt9z3bjmkimqoZLFYWCIWdhWzfDtp3kPKEQ1nUqAJ5KZ6z3UvWCIffzQlWlbU3j6WHkmTK0ne2VjScYIeHJ6VlXRUPVKm0VSix+LjpPNperOnTBjjKIlXKe4E3BKJvLWUa/kvEMBgL0dHcwNjHN2OSMNU+1FirM29LUWPHLW3/feex5as4uTqX5wIvCkHPD7FZlyLkwfMGYODFyuXzeCmXXWlvlLZ0aOCeVtgqvnmetXGowGLSuy2mvp/NWFfeVvE49T2U8nfZ6qnVDwWDVCIKaFwrOi4ZqVdoqOjwM27ak+Vorh23B3q7iwvOsIs2naG1uIhwKInHu1UopmTLDvNXkAlUxkhvt3KmZOaSURMKhqr9vLyIM1fKd1rrmYxPTs67yy2cDbTwXGe5Cq4unYMiptq3bNhi7ga1Vbatyoo69QwcjwxTK86x32Bbs4VWHa7sN2zqU/pNSOioYgoJxjTn0au35zlpeflurKTzvMH9oVdq2VPcO/YRtO2p4nkrByI3n6STnKYSwBBSmHRq5eCJpRWuqGk8PYdtJK99Zvk1FYYkwTMw4XtvKd1aZt9re2kJjQwQpJeMeRBjqiTaeiwwvnqdT45lx06oSdut5pp2vrcKrDnKT6phQMFgzZOUnbFuLJpderRvDXBhL5szIKU+y1toq9OrU80xnCqpINcO2SmXIob5ttTmeZ67trl2llrqQQrWUuFEZqiWQYK3tUihBSllUbVsNFfZ0WnGrPMn21paqKmUqbBuNxR1/hiiDVanSVrHELBoan3LvefZV8TyFELaRagub9zwvjKcQ4oNCiGNCiKQQYrcQ4toax/+MEGKfefxLQohdJY8LIcQfCSGGhRAJIcSDQogNZ/cqnHE2c54ph8Lw4EG9yKzkdTLtwI3KUEHXtraBc9uq4kQU3lr7LFXbgvvCHrcKQ/OJpKPiG3totVZ0osVlzjPmQCBB4Vbftpa6kMKLytCcg4IhKBhXpxJ984mkMQQb6KyxtioamnLoeU7V6PFUNDc1Wu9Pp+0qEzbPsxoFEYYZR+tCwfPsL5H8K6W/z/RqFzjv6UWe78dCiHcLIWr/FdQBIcQ7gM8CfwhchaFudJ8QYkmF428A/gX4EnAlcA9wjxDiMtthHwN+A/gARovNvLlm9U+jc4BVbVtD1xbcKwxlXIRt3Y4Oc5rzNI5x4XmmnK/r3ngWVIBqYeU869znaT/GiYcopbTCu06Np5qsUgt7vrNWaNXKeboO29b2PN2OPCvkPKsbIbcqQ3kpLeNZ2/N0V4ykvMjW1mZCNf7Wu1zOIq3V42lHHeM0dKvEGiq1qShUu8rUzJyj3GReSis/Wi1sCwU9XT/DvOuBF8/zeeAvgBEhxBeFENfVeU+lfBT4opTyy1LKVzAMXhx4X4XjPwTcK6X8cynlq1LKT2Bo7/4aGF4n8GHgT6SU/yGlfBF4NzAIvPWsXokDMq48T2/eoZuQsJMPXSml1XriLGzrfO2kw0pb4xjDCGVzOUeC9gVdWwdhW6ul5Cx4ni7Ctql0xpLEq2U8w6GQFep2YpidVtqC+z5PJ20qCrfDtpUXWU18QeFGZSg2HydvzkKtVcnbaQvbOultdFJpq1AVt1OOjadhCKu1qShU+NWpV6vaVCqpCynaWpppiBi5SSde7excjHQmSyAQsHKxlVgsYVvXrSpSyg8LIX4T+CngPcAjQohDGPJ8/yylHK3X5oQQEeBq4NO2588LIR4Erq9w2vUYnqqd+ygYxjXAAAVxB6SUs0KI3ea536ywlwbA/knYBoY357TPshR1nv18FRYMBoM11w2YzkE6nXW0B+XpBQOi5vGhoPG9KplK1Tw2ncmgPi4ConjtcteovhjMxxM1146b3lAkHHL8+wCIxeI1jYAyKE7WDtnGe9W6PigYwnCo9uuovkzE4vGaxyqPKRQMIpA1j29uamQuNs9cNFYzZGqp6TQ11rzGxkhhz+l0uqanOmcp9UQc7RkgGp2veayU0jKyTQ7WbjM95umZ2ZrXqLyx1pYm8rkc+SpfyJoaIgiMFpvZ2bmaXxKUEepoa629Z/N1m5qZc/R3rkKrHW0tNa+xq8Mw+uOT0zXXzufzVki4dO1y9PV0cGp4nOGxCbo7q39JGB41BOF7OtvJ5/NVNZPVWmMT06TS6SLt6Ep/j25weq5w2ohccQEjfPpLwO8BQeAHwF9LKX/ka2Fj7UFgCLhBSvmE7f7PALdIKbeXOScNvEdK+S+2+34V+KSUst8M6/4EGJRSDtuO+VdASinfUWEvnwI+WXr/3XffTXOVKRFuOTI8zYnRWZb1trFhefVhwPFkhqf2DREMCF57+aqaaz/5yimS6SxXbhigo6W6YTk9GeXAyUl62pvYura/6rHpTI7H954E4JYrVtX8IH356BgTs3E2LO9mWW/1b5nDk1H2n5yku62Jy9dV3wfAIy8eJ5+XbN+8jKYaikT7TkwwMhVjzdJOVvV3Vj12JpZkz6ERmhpCbN9ce+DxIy8cJy+d7eP46AxHh2cY6G5l08reqsdG4ymePTBMJBTkhstW1NzH0/uGmE9muHxdP91t1Y3nqfE5Dg1N0dfZzJbVZbMiFrlcnkdfOgHATVtXWl+2KvHC4RGmo0k2rexloLt6CHR8Zp69x8Zpb2ngqg1Lqx6byeb4ycvGe++1l6+s2Td88NQkQxNRVixpZ91g9dzaxGycl4+O0doU4ZqNtYcv/+TlE2Syea6+ZCltzdUjDoeGpjg1PsfyvnbWL6u+j0Qqw+5XhxACbr689t/X7ldPkUhluWJdP101XvOhiTkOnppy9Hfudh+vHB9nbHqetUu7WNlfPcw7ND7HwSFn+8hLyaMvHkdKHP19uSUej3PXXXcBdEgpK7rNnkQSFGbhznuBdwJjwFeAZcB/CSH+n5TyN/2sv8j4NMUebRtwaufOnbS3O2sCLiWTyfDAAw+wY8cOKwT7vQd/wonRF9iwYT27XlfJuTaYmYvx1L6vAYJdu3ZVPRbgmYNfhnSWW26++YwJ8KXs2XuQAycfoLOru+bak9OzPL73G0TCIe68886a1xj93kNMvLSfDRs2cuv1V1Zd+ydPv8j+k4+xfNkgu3bdXvManz34FaLzca6/4caa1zjxnXsZmYpxxeWXc8PVl1U9dnhsgj2H/pVgMFz0+yh3fdlsjof3/AMAu954R83Q7RPPvszR4Ufo7e1j1643Vj320LFTPHvgP+nq7HD0mh+fvIejJ0+zdesVXHHp+qrH3vfj3RwammLDurXsuv3mqtcI8OSrXyCTzXLjTa+tWX25/x+/BdEkN15/HRvXrax67NETp9l77B7Ckcaa1zg2Mc1PXv4XGhsivPlNb6p6LMDDTzzP0H8/Qd+SAXbtekPVa3zyub28fHSMFcsHHf2uD458m6GRcbZecSWXblhd9divf/c+To3PcdUVW7np2iuqHpvL5Xhq3xeQUnLzLa+rGkKWUvLYy18A4I6dbyh6Xcpd4/4jJzj4zf8iFGmqeY0Hj55k96tD9HV3nvF3Xo7wo0/zwKNP07NkgF27Xlf12Hvue5SDQ1NsvXQTu15f/XMP4MDpbzIyPsWlW7exeX3Bcaj0XnXD3JyzELZr42l6mv8fhtHcAHwP+DngPmm6sUKIrwD3An6N5wSG7F/pV5F+CpNdShmpcfyI7b7hkmP2VNqIKUNoJaXUt65wOOz5RVLY18jljHBFY0NDzXWbzRBULp8nEAgQrNHKkcnkrPNqr91onpOteazK8TREIhWPtV9jkzUfM1dz7UyuMFDaye+5sTFCdD5ONpevebzKubY6+H20mj2EyVS67LH260vacrmtLc1VB1ZDoUo0mc7U3rOZE29ubnT0+1D5w1Sm9tpqjFpba0vNa1RrT89GSTl4j6h8akdHW81jO81cXWw+UXvPZgFXe4U9n7G2mWOMVljbfo3WuLD22nsGI385NDJOLF5737Nmnranu9PR9KSOthZm5mJE5xOWJm3ZdaPzZLM5hBD09XSV/UywX6PKH07PzhEKVdfTVgO5e3u6HP0+BszRZJPTszWPnzBD5ANLehyt3d/Xw8j4FJPTc47eq25wep6XgqFTGNq2XwWWSynfLqW8VxbHf18EnvawdhFSyjTwLHCbuk8IETBvP1HhtCfsx5vssB1/FMOA2tdsx6i6rbTmOSOddSHebjsmU6OiTUppE0lwUTDkKJfqvCIW7AIMLqptHRT1gLuKW6tVpYaEHhQMfjaXq/m7LqgLRWoaTigIvDsp6lHHlA4hrry280re+YQzdSGF1a5So2hISmkZIjetKql0umYb1pxDaT6FG5UhSyChRitJYW0llFB7bTcFQ/bjpmpo3Kpioc721ppfpsEoRhKYs0VrFGhZbSo1ogwKVXHrpF1FVdr21ai0VfSbvaALWXHrxXjeJqXcbFazjpc7QEo5J6Ws7qc757PALwoh3iOE2Az8HdCCMd0FIcTXhBCfth3/V8AdQoj/JYTYZOYqrwE+b+5NAn8J/L4Q4qeEEFuBrwGnMdpaFhSVrHZSEWsUjRik0tU/ZHK5nNXr56ZVxYlIglNdW4U1d9PR2s5E4a21XQglFEQSalfENpgFIcZ51atiky4qbcHdQGyn48gUKjrhZKanm2pbsLerVDeeiWTKKgBxot3b2BCxPvhrfaBHHYrCK9yoDKkiJyWAUItCu0r1qthsLlfQzHVqPNVoshrTVVRBT60wuiIUDFp7qFUV67TSVqGMZ3Q+XlWZK5XOWK07/VUEEuwor3YhNW69TFV5VAgREkK8QQjxy0KINjCKe4QQzt5l7p7vWxjh3z/CCKtuA+6wVfWuBJbajn8cuAujiOkF4O3AW6WUL9uW/QzwNxiTYZ4GWs01nYuAniUKIgm1DZEQwjJYtSrEUrZv8PUeSZZ26Xm6URgqDMJ2aJhdKAG50bYNCGHtu1avp5s2FShUlzqR/rOk+WqMIyusbc7ddNSq4tbzVMLz1Y2n8kwbGyI1exrBeF+3OWxXsYxnDWk+hRuVIScTVeyo6Si1NGhn52JIDMPl5MsEGKFjqK0yVOjxdGaUwblAvHq8p0aPp6KxocH6UlNNSk8Z5WqC8KUooYTRiSnH03fqjZec5yqMfOZKjNaNB4Ao8Nvm7Q/Uc4MAUsrPY3qOZR67tcx93wa+XWU9CfyB+W9R4UaeTx1nTKyo/kGgjGvQQW4U3IkkOB2E7Wlth4OwrbXdGGYXEnpgfPgnU+manmfChewfFEKrqXSGXC5X9fWxpPkchm1bmpqKzqtGzKGurbW2w7BtQXvWeVV6W2uzkeOrIZRgrd3qbM9KZSiTzTIXm68agnQbtu10OJbMrmlbq2JVYc31rBG2nXQwx7MUJ0IJeSkdqwvZWdLTRTQWZ3xympWD5atolbJQNVm+Uvp6uhBCkEimmIvFHb9G9cRL2PavgGeALsD+V/PvnJlr1LjETV7SOM7ZZBU3Agn2dbO5XM05k25GnYFNYciBd2iFbZ3mPCPOcp5G7tIoRnJq5CwN2jp7nqVjyarhdByZwvIOaxhPQxRehW2dGSJVjFRLoi/mYAj2mWs7k+hT0ny11IUUTlWGstms5VHXUhdSWDM9Z2NVvaFph5q2drod6tuqnKcTdSGF5XlWEUqYi86TyWYJCGHJBTqhIBBfWdBA5S37HeY7wRAAUWIKYwuU9/RiPF+Loc5T+glyDKNNReODjIuwLdj1bWsLGYCL0KrtuFqhW7cFQ24UhtzI8wE0OpysYvcenRpPJeNXy8C59WgDgYC1h1oTUJwOwlZYurk1cp7JVNrKSzrNeTodiO1GXUjR5lAcXonCu1nbicqQWjcUDDr+XSsjm85kqr7/3BYLQcHzrKUEpPKW7oxnbc9ThVa7O9sdRa4UBYH4mYrHjDsQhC+HqhQeWaC8pxfjGcAQQyhlOUb4VuMDFX51om0LzseSpV0UIoGhcKSUO2oZOWuiikOD3+Ci2tZNXhIKHmpN45kq7NlJRSwUZPxqh23dGU/7sbXE4a2CIcc5T2eepzKADZGw4/ee04HYUQ+eZ0Hf1plXW2uiip3CUOzKnqc1x7OtxXFotSEStkLw1dZ2MsezFHVsMpWu+OUtl89bIeNaEnd2nOQ8VcWsm5AtOBOIH3Uwiqwclkzf+KSr8+qFF+N5P4Y2rEKahUJ/iKEupPGBWyNXGEtWI2zr0qO1FyPVNMxuPU8rbOum2tap5+nMeKrQa5ODNhVr7UZnxUhOB2HbsYqGanqe3sK28USyquaq25At2MeSOfM83eQ8CyHhyt5hXkpbta3znJeTdhW3+U5rbQd5T2sUmcMqXjD+ZtR7pFLodnYuRj6fJxgM0u4w1AwFDdzZuVhFEXe3rSSKQrvKdNlQdj6ftypmly6prqhWylJTBOX06ISr8+qFF+P5v4AbhRCvAI3A3RRCtr9dv61dnLgRhrcfVzvn6S4vCfbpJ/UO27op6nGZ83RYbZv0YOCU51krbFswzG6Mp1PP01vYVkpZNcccsyptnQ8WclowpAygq7Bta+3JKvYvBE76RxXK86zWrqKGWrsxQmCb61nFeE7PGo+58Tztx1cqGrIqbTvaivRea9Ha0kwkHEJS2TCrWZtLejudbxjDqw0IQSqdsULhdiamZ8nmcoRDoariD+VYvrQPgKGR8apauGcLL60qp4ArgD8FPocxZeV3gCullGP13d7FR6Ha1mnO06F36NIoA0QiZki4lvHMuC0YKnjL1d70UkrLMDvOSzrs8yyMI3PueTa59DzdeLVNjbU9z2wuZ73OTqttw6FQkRB/Jdy2qUAhbJtMpatOsSnkPN0XDFXLeapioZamRld5OJXzrDY+bM6j59nZVt3zlFJaxk8NuXZKd41eT2U8u1zkO8GIMhVCt+UN85jDcWGlhIJBa66o8l7tjIwZIdeBJd2uDL7aSyQcIp3JupobWi88DcOWUmallN+QUn5MSvmrUsp/lFI6m02kqUjOVtlad8/TxSxPhfI8axnPtOtWlcJx1Yx+Jpu1Qj1uRRKc5jydGmXjWGcFQ15ynk6UgBLmYwJnY9SstR3kPb2EbZuaGq18YLW5np5aVaxq28ofK1axkIt8JzgN25oCCR49z0rGM5lKW18IOzvcrd1Zo9fT6sN0aTyhet4zl8tZRnWJy6IesFXcljGew6bxXLqkug51OQKBAIMDhvd5avjc+21ehmG/Rwhxp+32Z4QQM0KIx80eUI1H7LJvTo2nyo3WCq1aYVs3nqdau5YAg8uwbThUKEZKVsl72r1Hx16tw5Cwl7yk8lIdh4TdFAypqtgqhlmFbJsaGxwXOUHBIFaruLU8T4ceLRjCEer4+Sq5SS+epwrbxuOJiq1SMWsItjvj6URlSBk/9wZOhW2re4etzU2uvshCQfigkvGccCliULx2e9Eadian58jn80TCIddfJqCgSFTOOxweM/KVbvOdiuUDxvSf88J4Ah/H7O8UQlyPMWT6Yxgi7p+r39YuPlRoVQjhOAzV4LDP020LDDiX6HObTy0qRqqydtJmlJ2GdNScSaeepxsPrukceJ6JKjnPeatYyLmBsx9fzauddymQoKg1FDuXy1nP6ybn2WJ6tbLK2nMx9x4tOFMZsnox293lJTstDdoa3qHL/B7U1rcd81i1aj+nXM/kmK1YyGnlcdm1q3ieA16Np5n3PDVSVin2rOLFeK4ADpk/vxX4NynlF4DfxegB1XjEnpd0+iYNO662VVW8zo2nU4k+t54n2MKrVTxE5Xk6Ddna102l0lWrS714h049T6Uw5K5gSBm4yoZZGdYmh5W2pWtXDdsm3Idtwa5vW37tmHm/EMKV0Q8EAlbxUiWhhIK6kLu8pKEyZHw5LRe6zUtpFfy4LepRRnFqeq5sdakynt0uWkkUVsFQmZynlNLy7LwYz/4+s+2jjJiBH6NsX1vlNxWZbJYJc89ePc8VgwXPs9rf+9nAi/GMAepKd2LI8wEkAXd/eZoi3IjCKwo5T2cFQw0R52s7ldFz205iHFu7sCflUtfWONYwLJIaXq0Xz1MVDFXxPHP5vPX7cNXn2VTbq3VbaasohFYr5w+9VNsax1fv9bTUhZqbXBeE1BJKiLqcqKIwVIYM77OcylA0FieXzxepETlFSe5lstmy+1a5wx4X8nkKFVqdi82fMdlnLhYnlU6bxT/uDbMSZJ+anj1j7YLn2el6XYBlZl5yamau6D04MjZJXkqamxodK0SV0t/XQ0MkTDKV5vQ59j69GM8HgH8UQvwjcAmF3s4tGC0rGo+47cUEm8JQjTFZypC48jzDzlpV3BYMGcfWNsxJD55nOBS08oHVPES3LTDgrGDI/pyuinocFAzFPYZt7b2elZifd19tC7Unq3iR5rPWtoQSyhtPr2FbqK4ypPowO9paCLrILYM5pcTMe5YrvpmcUWFb9waupbnJiqxMlCj2qErWns52R+L7pbS2NNPU2IDkTCm9EVOEQHmQbmlqbKDXzMOetOUmj58yRiuvXNbvKRwMhlb3ulWGsN3Bo6c8reEVL8bzgxhzL/uA/yGlVL741cC/1GtjFyNuReHBJpJQy8C5mBNqra1aVRwWDEXceLWWRF81z9N9RawQwlG7SsKlhJ5xbGGqSiXtUuWVRsIhV+0TTQ7CtgV1IW9hWyc5T9fGs0bO04s0n0J5lJXXVgVD7r2WaipDSgHIjfasnWptH25HhtkRQthyk8UGzq93KISwQqd20YF8Ps+weVt5kF5YYYrCnxgate47NmQYz1XLBjyvC7BhzQoAXj101Nc6bvHS5zkjpfw1KeVbpJT32u7/pJTyT+u7vYsLt6LwYDdwDkUSPLSqVPMOc7mc1eMXceV51lYZcqtrq3DSrpJ0OfnEONYwWvl8vuJAbC/FQuCwYMhDOwnUznnm8vlCUU+dw7aFNhX3nmetsK3q03QbtoXq7SqW9qzLYiFFT4UpJXkpbVNP3HueUNB/La1cVbfdKgDZWVamcnViapZ0Jks4FHI8x7Mca1cMAnDomOEdSik5euI0AKuXL614nhO2blpnrH30VNXe3Xrj3r8HhBCNwOXAEooNsJRSfq8eG7sYcatrCy6mqqiQsAtD5EQ312603eU8a4dt3U5UUTgpRkp6UAFqiISNClApSSTTZb+IJDwUIkHB80xnsmRzOUJlvFa/YdtKRT3xeAKJ0T9a77CtNRPTg3fYWmWmZzpTmHrixUOspjLkZeqJnUo9k3PReXK5HAEhPK9dqXJ13KOIgR27Yo/i1IhhSAf7e121R5Wyfs1yAI6dPE02m2VqJsrUzBzBYJA1Kwc9rwvG73vNiqUcPTnMA48+TdM5KhzyMs/zDuCfKRQN2ZGUF43XOMBLO4nbqSpuvFonIgnKwAUDgbIf+LXWrha2LeQ83XmelmGumvM0jJwbwyyEoKkhQjyZIplMlVWfsTxPF/2jxvERBMYfUCKRLFtB6nbepqLWTE9lnFqam1x/QNYaiK08O7cyd1DdwM2YodVIOOw6jG3sp7LKkFrbjfasHavitmQCijKmnR1trnOpClXYM1xSuVqQz/NvPE8Nj5PL5wkGAhw7OQwUwq5e6e/tpqOthdnoPPsOn7A85XWrBl3/fZdj5y3b+Yev38NTe16hvbmBO/J5/K9aHS+v4N8A/woslVIGSv5pw+kDbzlPl9q2rlpVagswqDmhrg2cg7yk20HYCidhW6udxK1Xa830LJ+bLHi07tYNBALW2pWEEpSBanYdWq1u4KJKuN1D+LO1hhKQV5k7sOnElslLWuLqHa2eik2qh229takoKnmefno8FUpRZ2Rs0hKPSKUz1qDqJb3einoABvp6aGyIkEqnLe/ziBlaXbvKn3cohGDblksAeGrPKzz53MtAIeTql03rVvGm224kEg7R0hj2/OXEDV6eoR/4rJRytOaRGld48Q6dattaXq2rVhXnnqebcLCxtoNqW6tVxa3xrK5Bm7e1k7gNr1qGOVl+ba85T/s5lTzEuNdeTNPApdKZsu8TPxWx9pme5Yqo7KO93NJZRWR9xmMfpkIJJZQbXO2/YKgwpcSeG1cVsn6MZ09XBw2RMNlczgrVnh4dR0pJe2uLpy8pioCtcvXQ0VPMxxNW8dBan6FVgGu3bUYAL+07zNjkNA2RCNdcvtn3uorbbrqG3/nV/481S717327wYjz/Dbi1zvsoixCiWwjxDSHEnCkB+CVz/Fm14/9GCLFfCJEQQpwQQvy1EKKj5DhZ5t87z/4VVcetKDwUFIbOikiCg5BwykObinG882pbr15tJeNpN9huPc+mGp6nH+NZSyjBi4QeGAZfhdTLtX1EPQyUVqgCo7yUZVt4VCuIF+OpDFwylT7jtZzxqACkUOLpqXS6qJAqk81av48uj2HbluYmGiIRJMXepxraPOCx5QMMSUSlA3tq2PAOT5428pLLTcEAP6jK1Zf2H+alfYeRUjLY3+tJlq+Uwf4+br3+Kuv2W29/rauCPSe0tjQRCZ+bAKiXgqFfA74thHgt8BJQ9MkqpfzremzM5BvAUmAHEAa+DHwBuKvC8YPmv98EXgFWAX9v3vf2kmPfC9xruz1Tr017xcvkk0LO01nBUEOdq229qAuBvaintrat2z+wWmFble8MBoOue+Kcep5uNHMVlmEuY4TSmQyZrFHV7NbzFELQ1trM9GyU6Hz8DM/Hy7xNRSgUoiFihPpi84miYqZUOm19AfJSMNTYEKGxIUIylWY2GqOxoWB0pj3MxLSjdFpnozEmp2YZ7DdKOFR7SWNDxPXvWSGEYKCvm+NDIwyPTTLQZ6w94lOKTrFq+QDHTg1z5MQQ11y+iRNmy8eKpf6N57Ytl/Af9z/KsZPDVr5z25YNvtdVvHnHTVy2aS0Aa1cuq9u6C4EXz/PnMJSF/gfw68BHbP8+XK+NCSE2A3cA/1NKuVtK+Zj5fO8UQpSNIUgpX5ZS/g8p5feklIellD8Cfg94sxCi9FNyRko5YvtXfZDiOcBTq4p5rH0iS9m1LZEEN20wtUOrXuaEFq3tRGHIc86zknfoLd8J1Q2c/Tn9eZ5nvhWV1xkMBj0VWLRVERzwE7aFyvlDle9siIQ9exiVppT4DdsCVuP+xPSMdZ8KrfZ2d3pu3IfCoGbVI5nOZC3RdWVMvaJCq4ePD5GXkgNHTgKwfvVyX+uC8VpuNY0bGK/dtdsu9b2uQgjB2pXLznvDCd48zz8FPgn8mZTybE4gvR7DwD1ju+9BIA9sB/7d4TodwJyUstQ1+1tTJekIhnf6ZVmp8x0QQjQA9k/ENjBk8WpJ41VCnaf+V4YkGAg4XtP+5x2PJ8p+SNkNa0DUlvJTBMzF0+nK1zifMD7Uw6Fg2WNKr1ERChZUgCqtrQxUMOj896H2os4vd54yFg0NEdevnSq4mo8nil579b9qB4mEQ67XVqL2sfn4GeeqvF9zUwPZGmpS5VBFQzNz0TPXNo1eU2ODq9dQ0dbazNjkNFMzs2QyhapMFbJsa232/DfS0drC6PgUk1MzZFYU+gHVTMy2libPayvDOzYxba2hhNF7uto9rwuFqthTI+NkMhlOnR5FSklzUwNNDWFfa68cXIIw973n5f1E5+NEwiGW9fdWXbfW66h4687XMjYxTSye4K07b6a5wvtiMeL0Gp2sUQtRxV6UP0GIKeA1UsrDHvbl5nk+DrxHSrmx5P4x4JNSyr9zsEYv8CzwdSnl79nu/wTwIyCO4UX/IfCxaiFnIcSnML40FHH33XfT3Ow+3FWOV4+PMzo9z7rBLlYscVZUIKXkxy8cB+D6LctpKONZZnN5HnvpBACvvXyl40q0dDbH4y8b32pvuWJV2W/iJ8dmOXx6miWdLVy62rkCyUwsyZ5DIzQ1hNi+ufw35idfOUUynWXb+gE6W53n+YYno+w/OUlPexNb155ZYj85G+elo2O0NkW4ZqO7Qogjw9OcGJ1lWW8bG5af6UE8f3CY2fkUm1f10t/lLqR4+PQUJ8fmyq49FU3w4uFRWhrDvGaT+2/t+05MMDIVY/VAJ6sHOosee/bAaaLxNJetWUJvh/v38ivHxxmbnmftYBcrbe/b0ekYrx6foKOlgSs3eGuE339yguHJGKv6O6xCECklj710glxecu2mZTQ3emtKOD46w9HhGfq7Wti8qq/i83lBvb8joSDXb1nO0ESUQ0NTdLc3cXmZ96RbXjg8wnS0EKHo7WjmsjX+w7YaiMfj3HXXXQAdUsryI2zw5nl+FXgH8L+9bEwI8WfAb9c4zHcJlhCiHfg+Ru7zU/bHpJR/bLv5vBCiBfgtoFq+9tPAZ22324BTO3fupL3dm1pIJpPhgQceYMeOHYTDYSa+cy+j00e4/PKtXH/VZY7XeeKVL5DOZLn55lvKVvLNxeZ57KWvIoTgTXfe6TgclclkefzlLwBw2xt2lPVqH3j0aQ6ffpo1q1exa9etNa9RcWpknD2Hvk0oHGHXrl1ln//p/f8EZHndrbe4mrrw4quH2H/yfto7usqu/dzL+3np6EMs7V9S8bkr8fCTz3Ni9An6B5aya9cbzri+fV/8JrPzKW664Xqr+MIp//3Ec5wce5KBpYPs2nVb0WN79h7kxcMP0L+kz/WeAQIPP8nI488xuGw5u26/uXjtI18D0tx682tZuezMD/ZKr6FC/uhxxp7cw7LlK9n1hhut+x/ZvYdXj0+wZtVKdu3a4XrPAE1PPMfwfz9JT18/u3a9ATAiCj9+4UsA/NSbd7mei6l4ef8Rjn7nXoKRJnbs2MEDDzxApLEFiHH99mu46rKNNdeoRCaT5aXP/iPpbI7t19/EfY/shqEpXrNtK7fddI3ndRUr9h3mn797n3X7Z3/qdlavqP4FpdbreCFQj2ucm6toL4vwYjyDwMeEELcDL3JmwdBHa5z/f4Gv1DjmCDCCoWBkYeYtu83HKiKEaMMoBooCb5NS1vLDdwOfEEI0SCnLJrPM+63HlAEKh8O+34hqDSVz19TY6LIqNkw6k0VSvpo2L9VxIVcSeqFQyFLUycvya2fNIpampup7Lv09qUb/dDpb8TxVTNTW2uzq99FiRgLSmUzZ81Jmb2pzc5Pr167VWrt43+r6VOVme1ur+7XNvGQylT7jXK+/C0WHWZU6n0gWnS+lLFLqcfMaKro6jC+P0fl40eNqHFlnu/vfhULJzU3NRq01xs2expamRuu19oJq/B+bmCYYDCKltCpiVwwO+Pq7DofDrFo2wJETpzl6cpgjx41+yXWrl9fFcF152UaGx6Z44rmXef0NV7Nh7UpXe7tQjafCzzU6Pc+L8dwKPG/+XOoe1YwBSynHgZqzY4QQTwCdQoirpZTPmne/HqPIaXeV89qB+zAM3U85LATaBkxXMpznCi8iCfbjK1XcehFIAHNodThMKp2u2K6S9CDeDoV2knTaEFkv9YYz2Sw588uE28rVWtW2KpfqRZlGzfQsVzAkpSTuUX/W2I9ZMFRm7UKbireinkoFQ6l0oYrXS6sKFFpKStV6/KgLKezzMRUql9rlUR9W0d3VQTgUIpPNMjk9RzKdJZXOEAwGrZylHzauXcmRE6f59vd/BBjGvh79kmD8bd552w3cedsNdVlP4x7XxlNK+bqzsZEyz/OqEOJe4ItCiA9gtKp8HvimlPI0gBBiGfAQ8G4p5VOm4bwfaAbeBbSb9wGMSylzQog3Ywg9PIkxg3QH8HHgL87FdVUj40HbFmrr22Y86NoqGiKm8axQceu1VUUdLzE8xNKKWvvMTLciCbX6PFU1qychgypjyVLptFWY5XYuJhjFQACJMn2eBVF49+tCQT2oVCdW3Y6Ew55l0jrajWpb1XupUNW27W3evcMe2wzLdCZDJBy2Jor4NXABIVja38OJoVFOjYwRSxjvl4G+blcTcSpx7ZVbuO/Hu60hzdu2bKjLuprFwdnXMPLHzwP7MAzkD4DHgF+yPR4GNmIYS4CrMCpxtwKHgGHbP5WAylAYq7YH+GXgoxhFQwuKlz5PqC2jl/LQAuN4bY+9mJFw2KoULrd2wjbL0+0QZaeepxpA7QZ1Tvl2EuO+cCjkKQ+nxOHLTVZRVc1eew8reZ5+21SgMCFkdi5W1C417VPIAIz2HSUKoYzm2KQRWvUjRadQ0z6OHD/NTMz4vZfL+3qhs72VO269DoHhQd9+63V1WVezOPA6VeUa4GeBlUDRp6aU8qfrsC+11hSVBRGQUh7D1q0hpXyY4u6NcufcS7E4wqKh0OfpsmfSqefp4QO9lkSf115MIQQRs7k+lUpDSRN90lLqcd8fqIxnLpcjm82eIYRQCNu69+JaHPRievUOqw3EVobZ7UQVRWtrIZ+ayWat6EbUx0BpRVtrC8FgkFwux+xcjO7OdnL5PNMzanal9/CqEIKBJT0cPj7E8Ngky5cuYXTC/wQRxbrVy3n4yec5dHyIuPn6bVq3yve6ih03X8s1V2ymuanB9d+IZnHj2vM0Zewex6iIfRuG97cFIx955uh0jWP85zzLGzgvAgkKFeqtmPNUHmKDt5AwlFcZSloerXvv0B5+LLe2Fbb14Hkqzy+dyZ4x03Peo/asQhnGTDZ35toeJ6pYazc2WC1KMVvoVuUpVejVCwEh6DZ7JtUkkdm5GHkpCQaDvnKeUBAVGBmfJC8lo+PK86yD8Vy1jHAoxOT0LIlUlkAg4LpKuhZdHW3acF6AeAnbfhz4iJTyzUAa+BCwCWPSyok67u2iw3PYtsZkFcuj9ZLzDFdXGSroz7r/cFAGt5xXaxlPD55nIBCwvOxyoVtfBUMNhTByqYfo1ztsbIjYDFzxlBKvE1UUQgirIMgeulXjtzravIdWoVC8o4ynNUGks9112L0U1aZ0enSC8Ylpoxo5FKLfh0asoqmxgeuu2mLd3r5ts6dcuObiw4vxXIfRPwmG8WwxlXk+R3E+UuMCKaWlbOG2xLowWaWS8fRmlMEu0VdBZD1l7NmL/JoyuNUMnBcJPft+ysn/KeH1Jg9GTghhGcfS4dJ+vUOlQQsQtUndSSktg9fuYWyYQq09V2Q8lcydP+9Q5T2V0VT/d/sI2SpUDvLYyWFrRNbKZf11Gzt1+63XccPVl9Hf1cLtt2yvy5qaCx8v775pTGk6YIhCu0onhcIdjUty+bxVlefV86wYtvWYSwVb2DZd3jB7nXxiP6ecYS54nt68gGpFQ6ogx4vnCZU1aP16nlC+KjaeSFqFOH5yk0qObtpWFatk/zo8CqwrVOWr6pMcM4cd9/oYv6VYNtBHU2MDyVSah594DoCVywZ8r6toaWrkrbffzOZVfb5eO83FhRfj+QhGewfAt4G/EkJ8EfgXjKpYjQfsIVf3xrN6XlIZPi+eZyFse6YRyktZKBjy4nk2VC5GsqaTeMwVKeNZ2lKSzeUsT9yL5wmF0Gmp8Ywn/BUMQfmqWOUpNjc1up4CY8cynjOFnsl6VMRCcV4SCoLoSiDdD4FAwBI9HzNnWNrFyzWahcDrSDL16fCnGK0fNwDfAf6kTvu66FCGLxAIuO4Fq5XzzNTD8yyztt3oecl5KiH0ct6hX8+zkneYsN32GhJuqRC2jfkQSFC0taoJJXbjafZLehjrZUeFVqdMUfW8lFbBkNfRXoqBJUb+cWJyhmw2y/CYYTwHl/g3ngC3XHclL+0z5LSXLulh9XJvWrkaTb3wIpIwZfs5D/yZui2E8P6pcZGT9iiQALUVhlIZ79W21UKr6r6AENYkEzdEqswLVcbTq4FrrtCPGbfmbUYIeMyZKc9S9V5aa9fB82wvk/NUYgNtPvKdYPc8DW9zdi5GLpcjEAj4Dtt2tLVaszcPnzhtGWW/sysV61Yt462338zYxDS33XSNr3FhGk098B4DsmGO6/og8DGgfsmIi4jCXEwPRT01w7beejGL1z7TMKtioYaGiKcPs8Yq80LVXEy/nmdp2Fap9/jJbSmJvPkzKmKTRY97QXme9pxnPYqF4EzPc9wMgfZ0tfsuvhFCsGbFUl49dJwf/vcTgNFK4nWOZzluue7Kuq2l0fjF8V+MEKJBCPFpIcQzQojHhRBvNe9/L3AUYxj2587ONi98lOHzUnhTSyShkJf0XtRTtp3ER7GQsZ/KFbFqYLXXD99KYdt40rs0n8Iq6olXMJ518Dxn52yeZ0x5nv7CtkonNjafIJFMMW4Ofu6rg9gAwEZTXOD4KWNuw+b1q+uyrkazGHHzdfOPgF8BjgGrgW8LIb6AYTQ/CqyWUv6fem/wYiHlUbzdOKd6ta2VP/TSi1nFO0xZAgneDFz1alulMOTNyDVVUOspeJ7ejWe5fkn7dBI/Oc9yFbFq8LN6zCtNjQ1WbnNkfJIJsyK2r7vT17qKSzesLrp92UZd1KO5cHFjPH8GQ4D97RgDpIMYYd8rpJTflFLmzsYGLxYKYVv3xrNmzlN5iB48z0g1z9OHUYZCGLls2LZunmdx2DZuyf75r4i1K/VkMllrpJwfz1MZyLlozJoqM2XmKP0aTyhUxQ6PTTI6YZQv9PV0+l7XWKeLXa+7HgHc9JrLWbfK/dBujeZ8wU2CbTnwLICU8mUhRAr4nCmQoPFJ2of+bM2wbcpHzrNqwZD3cLD9vKrVtnUO26pqW689nmAL29o8z5jpdYaCQc/DmQFaWpoJBYNkczlmo/N0d7ZbXmi3zxFcYEwM2Xf4OMOjkxwfGgUKcy3rwY6br+Wma6/QKj2aCx43nmcQQ1FIkQVi9d3OxUva42gvsHmetST0vPRiVi0Y8i7NB4XQarnxXvUqGKpUbeu1xxMKE0hi8QR5U7ygMLuyxVclaEAIy8Ocmpkjlc5Y4eB6GM/lg8Z8+Sefe5l4Ikk4FGJZHXox7WjDqbkYcON5CuArpscJRq/n3wshiibg1nOqysVEykfYtlovJkBSSej5WLtsXtKHUYbCVJNS45mzCxn4bVUprbb1qS4EBfk9KaUVFp6LGl5oR5u/oh6Azo42xqdmmJ6NWiHixoZIXYzSxnWrEEJYIeYVg0v0jEmNxgNuPM+vAmMYk1Nmga8Dp2231T+NB9J+Zm5aIglnep5SStI+jFw1eb6UD6MMlWdj2sO4XsO2qmgnYZO2g/pI6AWDQUsoQbWUFDxPf/2SUKiKHZ+csSpi6+F1gmH4164ctG5fu+3Suqyr0VxsOP6kllK+92xu5GLHCtv6qLZNlVMBymRQSWlPk09sI8nyUhZNyPAzUQUK3l8ylSafz1uiBWoQdiQc8uwVtTQ1IgAJxOMJq81DTSvxKzjQ0d7KfCJpCRgo41kPz1NNERkZn7TEJ+ohc6d4x5tv45v/+SDNTY1cc8Xmuq2r0VxM1EUkQeMfP2FbZbxyuRzZXI6QzeCokK0Qwpu2rW0/mUy26HbKxyxPKM47JlJpy5tLWipA3sOUgUCAluYmYvEE0fnEGeIDKm/plc6ONk6PTjBtjvRSRrS9jsZzeGzCei2X1kmpB4yq2F9/78/UbT2N5mKkPjN9NL7xNXPTbtBKiobSNiEDL4Us4VAIdVa6JO+Z9KFcBKoy1TDods1ZFcb1O+GitUxLifI8W31MJwHoMoXUZ2aNmjllPDvqELZVknaTU7McOzUMwGCdi3o0Go0/FrXxFEJ0CyG+IYSYE0LMCCG+JISo+ukkhHhYCCFL/v19yTErhRDfF0LEhRBjQog/F0IsqBduKQx5CNsGg0HLQylV60n6qOIF02NVRUMlYWEVtvUjwab6Le2FPfOW8fRXIGNVxZoGM5XOWL/nNp+ep6qInY0axnPKmk7i33i2tTTT3dmOxBBLEMDypUt8r6vRaOrHojaewDeALRgj0N4E3Ax8wcF5XwSW2v59TD0ghAhiDPOOYEyDeQ/wCxgKSguGH5EEsEndlXiHfttJoNBHeoZh9tmLCQUDaa+4jddhOgnYxnuZnue8rRfTz+8DipWA8nlp9WLWS3Dgiks3WD+vX7Pc1xxPjUZTfxat8RRCbAbuAP6nlHK3lPIx4NeBdwohBqufTVxKOWL7N2d7bCdwKfAuKeUeKeUPgU8AHxRC1E/F2iV+5PnAJnWXKg7b+lEXUjQ2lh8snUj6k9Czn5tIlPM8/YZtTQF302gqI9rW2ux7KkenzXgm01mklETCId9jwxTbr7zUCmnf9Jor6rKmRqOpH4u5YOh6YEZK+YztvgeBPLAd+Pcq5/68EOJdwAjwPeCPpZQq8XU98JKUctR2/H3A32F4uc+XW9CcHGO3Em1gtIeUaxFxgjovk8lYxjMYFJ7WUx+084lE0fmqNSMSDnvep/IsY/PxojWUwQuHghXXtl9jtbWjsXnrGDWOq6kh4nnPUKjmnY3GyGQyzJjeYUtzo691Abo7jPDs9GyUWKdh5Hu6Oshmy/faul+/jT/40HvJ5nI0N/nfrx9qvYYXAvoaLwzqcY1OzxWLVV1PCPFx4D1Syo0l948Bn5RS/l2F834JOI7Rg3o58H+Ap5R4gylmv0pKebvtnGZgHthleqLl1v0U8MnS+++++26am/2H1Ha/eopEKsu29QN0trr3uJ47OMzcfIotq/vo6yx4P6cnohw4NUlPexNb13qTYXvh8AjT0SSbVvYy0F3I6T3ywnHyUrJ98zKaPHq2+05MMDIVY83STlb1dwLw6vFxRqfnWbu0i5X9HZ7WBTg9GeXAyUm625u4fG1/XX4Xdh7fe5J0JkdfZzPjM3F6O5q5bI3OTWo05zPxeJy77roLoKMkalnEOfc8hRB/Bvx2jcM8N59JKe050ZeEEMPAQ0KIdVLKw17XBT4NfNZ2uw04tXPnTtrbvTWwZzIZHnjgAXbs2MHzh+8mkcpy82tv8lQccvqb32PuyEm2XLaVay7fZN3/4917OHDqcVauWM6uXTs87XPq3+9j+tXDbNy4iRtfczkA2WyOh/f8AwC73nhHxdCt/RrD5ULSP3qCkSefZ9nylezacRMAI9/6L0an57nqyit8NfHvP3KCA9/8LyINTezatYt7H36SA6cm2bhhHbtuv9nzuorTc9/jwJGTjM8YQY2rrriMN9x0je91Fxs1X8MLAH2NFwb1uMa5uYr2soiFCNv+X+ArNY45ghFyLbIiZkVst/mYU3ab/68HDpvnXltyjHJDKq4rpUwBVmJO5czC4bDvN2I4HLbk6Jqbmzytp4xXNpcvOj+bzZmPN3repxrunM5krTVUFa8AWltbisQTylHp96T6IhPJlPW4muXZ3tbq63fb123MqZyZixEOh5k1hdy7uzrq8uGxfGAJB46ctG5vXLfygv1Qgvq81xc7+hovDPxco9PzznnBkJRyXEq5r8a/NPAE0CmEuNp2+uvNPe8uu3h5tpn/D5v/PwFsFULYDfMOYA54xdNF1QE/wvD280orYv3O3ITyAu7q54aGSE3DWQ1VURudLwyWnk8YP/stGOo085LJVJpkKmXlPOsx2gtg22WXWD8LIVhZx+kkGo1mcbNoq22llK8C9wJfFEJcK4S4Efg88E0p5WkAIcQyIcQ+IcS15u11QohPCCGuFkKsFkL8FPA14BEp5Yvm0vdjGMl/FkJcIYS4HfgT4G9N7/Kck83lLP1V79W25WdjKg/Rq/4s2AZLlzGefsXKVQuGqogFu5CBv1aVhkjEKhqano1Z7SSd7fUxniuWLmHTupUAvO2OmwmFFnP9nUajqSeL1nia/DywD3gI+AHwGPBLtsfDwEZAVeykgTdgGMh9GCHi7wBvVieYQ7vfBOQwvNCvYxjYPziL11EV+ygxz32ekfJ9nn5E4RVlPc9EfYxnS4mQQTqTsVpi6tH20WmN95plds4QNKiX5wnw7v/xRm7aupLrrtxStzU1Gs3iZ1F/VZZSTgF3VXn8GFjqcUgpTwK3OFj3OLCrDlusCyrfGQgEinRp3aD6OEuNp1+FIagetvVrPNV4r1g8gZSSOTMvGQ4FfYkvKHq7Ozk9OsHL+46Qy+eJhMN1EW9XhEJBQsHF/h1Uo9HUG/1XvwiwpPl8GLjGSHkhg7rkPMuoANXNeJqeZy6XI5VOWxqxba3+hkorlg30AfDUHiOdPdjfa01v0Wg0Gq/oT5FFQNqnuhDYvcMS42nlPOsbtlVFPS0+i3oi4bAVcp6NzlsCCfWYTgIF45k3+5mXDWiBdY1G4x9tPBcBKmwbiXiPoisJPbuBA/9jw6BgIMsV9fjVn4VCVezMbNQK27b7nLepKO2ZXTagRQw0Go1/tPFcBFjjyOrgeSZTxcYz6XNgNRTGd6XSGTKm/JzSn/VbEQvFIuuW51knjdiOthbWr15u3d62ZUOVozUajcYZ2nguAtJpwyD5ynk2VAjbpuqQT22IWHlC5X3O19HzVLMxp2ejTM0Y6h4ddRjtpbjrrTvYcska3v/ON/vO0Wo0Gg0s8mrbiwU/g7AVlueZTCGlRAhBLpez1vZjNIQQtDQ3Eo3Fic0n6GxvIxZXxtNfzhMKnufMXIzxqRkA+ro7fa9bWL+d//lzP1W39TQajUZ7nouA+hQMGWHZXD5vhVbt+c96tpRAwQNtrYMofmG81xzjkzMA9PV0+V5Xo9Fozhba81wEFAqGvBvPSCSCEAIpJYlkmkg4bCkC2cOuXimdjTlfR89TDZA+dGwINeWnt9v7NBWNRqM522jPcxGg2kkawt6/ywSEsEQFkqbRjJtFPfXI8ykPMzafIJ3JWAa/HgVDywb6CAQCluHsbG/15YVrNBrN2UYbz0VAqg4SelAYLJ0wK26VhJ5fgXUoFAbNxxNWm0owGPRVxasIh0IM9hf6L9esHPS9pkaj0ZxNtPFcBChVoKYGfx5iqZhBPGl4ns1N/j1P1Xc5F523Cay31kUFCGCDrZ3k8s3r67KmRqPRnC10znMRoIynEjrwimU8E8WeZ1Ojf8/TElifjTIzV9/RXgC337qdudg8c7E4l25YU7d1NRqN5mygjeciQKkANfr0PK3QqpnrVDnP5jrkPLtsFbF2z7NeNEQivOun76jbehqNRnM20WHbRYDlefrMearK19i8IXGnqm2b6pDz7O5sB2BmNsb0TP09T41Gozmf0MZzEZBQOU/fvZhqsLThcSbqmPPsaDOmnGRzOU4NjwGFUK5Go9FcbGjjuQhI1dnzVBNPLPH2Jv/tJMFg0JqDeeL0KADdHe2+19VoNJrzEW08FwH1CtuqnktlNNWEkrY6TShZ2l88zmv5oJ5QotFoLk608Vxg8lJacnqNPsO29l5MgGidx3utW7nM+rmvu9OS7NNoNJqLjUVtPIUQ3UKIbwgh5oQQM0KILwkhKpZ4CiFWCyFkhX8/Yzuu3OPvPDdXVUwul7d+bvQhzwfF+rN5Ka3CobY6jfdat6pgPNesXFqXNTUajeZ8ZLG3qnwDWArsAMLAl4EvAHdVOP6kebydXwJ+C/hhyf3vBe613Z7xuVdPZE3jGQmHCAaDvtYqeJ5J4okkubyxdlsdJPQAVi0f4I5br2M+nuB1N1xdlzU1Go3mfGTRGk8hxGbgDuA1UspnzPt+HfiBEOI3pZSnS8+RUuaAkZJ13gb8q5QyVnL4jJRyhAVGGU+/PZ5g5DYFkM/nGR6dAAxpvlCoPi+zEILbb9lel7U0Go3mfGbRGk/gegwD94ztvgeBPLAd+PdaCwghrga2AR8s8/DfCiH+ETgC/D3wZamUycuv1QDYLVwbQCaTIWPOzHRLJpOxjGdTY4Pndey0t7UwG53nyPEhwAjl1mNdr6jnXsg9nE0u9OsDfY0XCvoa3a1RC1HFXiwoQoiPA++RUm4suX8M+KSU8u8crPH/gFullJeW3P8J4EdAHNgJ/CHwMSnlX1dZ61PAJ0vvv/vuu2n2MdNybHqeV46P09HSwJUb/OcRnz84zOx8iu72JqbmEnS2NrJt/YDvdTUajeZiIB6Pc9dddwF0SCnnKh13zj1PIcSfAb9d47DNdXieJozc6B+XPialtN/3vBCiBSMvWtF4Ap8GPmu73Qac2rlzJ+3t3vodM5kMX/mX7wCwYvkydu3yL083958P8tzLB5ieMypuN1+ynl133Ox7Xa9kMhkeeOABduzYQfgCHDN2oV8f6Gu8UNDX6Iy5uYr2soiFCNv+X+ArNY45gpG7LGokFEKEgG5K8poVeDvQDHzNwbG7gU8IIRqklKlyB5j3W4+paSLhcNjXGzGTNcK2rS3NdXlD93R1AqDiCUv7exfFH4rf39Ni50K/PtDXeKGgr7H2uU4458ZTSjkOjNc6TgjxBNAphLhaSvmseffrMdprdjt4qvcD/2k+Xy22AdOVDOfZJJPNAdDa7F9/FqCnq9gL7u/tqsu6Go1GoymwaAuGpJSvCiHuBb4ohPgARqvK54FvqkpbIcQy4CHg3VLKp9S5Qoj1wM3ArtJ1hRBvBvqBJ4EkRhvMx4G/OLtXVJ6MWTDUXAcJPYANa1YU3e7v667LuhqNRqMpsGiNp8nPYxjMhzCqbL8D/Ibt8TCwESM8a+d9wCng/jJrZjCqbz8HCOAQ8FHgi/XcuFOU59lSJ8+zu7OdttZmorE4S5f00NFWv7FhGo1GozFY1MZTSjlFZUEEpJTHMAxg6f0fx/Amy51zL8XiCAtK1sx51st4Arzrbbfz1J5XePOOm+q2pkaj0WgKLGrjeTGgwrb1mHyiuGTtSi5Zu7Ju62k0Go2mmEWtbXsxkFYFQ3WS0NNoNBrN2UcbzwUknc6QzxtNJfUaG6bRaDSas482ngtI1Jx6Eg6FiFzgfVcajUZzIaGN5wKihla3tjRZogsajUajWfxo47mAxMyh1W0tOmSr0Wg05xPaeC4gKmzboouFNBqN5rxCG88FYmxymn+/9xGgfsOqNRqNRnNu0MZzgZiLzqPGwXV3eJvKotFoNJqFQRvPBaK3u4Mbr9nKusEubrhm60JvR6PRaDQu0MZzgehsb+MtO1/LiiUdNDZEFno7Go1Go3GBNp4ajUaj0bhEG0+NRqPRaFyijadGo9FoNC7RxlOj0Wg0Gpdo46nRaDQajUu08dRoNBqNxiXaeGo0Go1G45LQQm/gfGdubs7zuZlMhng8ztzcHOELdCTZhX6NF/r1gb7GCwV9jc5w+pkulEScxh1CiGXAqYXeh0aj0WjOCsullEOVHtTG0yPCGMA5CER9LNOGYYCX+1xnMXOhX+OFfn2gr/FCQV+ju3VOyyoGUodtPWL+Uit+K3GCbQB2VErpPf67iLnQr/FCvz7Q13ihoK/RFTXP1QVDGo1Go9G4RBtPjUaj0Whcoo3nwpIC/tD8/0LlQr/GC/36QF/jhYK+xjqiC4Y0Go1Go3GJ9jw1Go1Go3GJNp4ajUaj0bhEG0+NRqPRaFyijadGo9FoNC7RxnOBEEJ8UAhxTAiRFELsFkJcu9B7qhdCiE8JIWTJv30LvS8/CCFuFkJ8Twhx2ryet5Y8LoQQfySEGBZCJIQQDwohNizQdj3h4Bq/UuZ1vXeBtusaIcTvCiGeFkJEhRBjQoh7hBAbS45pFEL8rRBiUggRE0J8RwjRv1B7dovDa3y4zOv49wu1Z7cIIX5FCPGiEGLO/PeEEOKNtsfPyWuojecCIIR4B/BZjJLqq4AXgPuEEEsWdGP1ZS+w1PbvpoXdjm9aMF6nD1Z4/GPAbwAfALYD8xivaeO52V5dqHWNAPdS/Lr+3DnYV724Bfhb4DpgBxAG7hdCtNiO+RzwZuBnzOMHge+e4336wck1AnyR4tfxY+dykz45BfwOcDVwDfAj4D+EEFvMx8/Nayil1P/O8T9gN/B52+0AhtTf7yz03up0fZ8C9iz0Ps7i9UngrbbbAhgGftN2XweQBN650PutxzWa930FuGeh91bHa+wzr/Nm22uWBt5uO2aTecx1C73felyjed/DwF8u9N7qfJ1TwPvP5WuoPc9zjBAigvGN6UF1n5Qyb96+fqH2dRbYYIb/jgghviGEWLnQGzqLrAEGKH5NZzG+JF1IrynArWY4cL8Q4u+EED0LvSEfdJj/T5n/X43hqdlfx33ACc7f17H0GhU/L4SYEEK8LIT4tBCi+VxvrB4IIYJCiHdiRE2e4By+hloY/tzTCwSB0ZL7RzG+IV0I7AZ+AdiPERL6JPCoEOIyKeWFOM1hwPy/3Gs6wIXDvRjhr6PAOuB/Az8UQlwvpcwt6M5cIoQIAH8J/ERK+bJ59wCQllLOlBx+Xr6OFa4R4G7gOHAauBz4P8BG4KfP9R69IoTYimEsG4EY8DYp5StCiG2co9dQG09N3ZFS/tB280UhxG6MP9afBb60MLvS+EVK+U3bzZeEEC8Ch4FbgYcWZFPe+VvgMs7/XHw1yl6jlPILtpsvCSGGgYeEEOuklIfP5QZ9sB/YhuFZvx34qhDilnO5AR22PfdMADmgtPqrHxg599s5+5jfAg8A6xd4K2cL9bpdNK8pgJTyCMb7+bx6XYUQnwfeBLxOSmkfaD8CRIQQnSWnnHevY5VrLMdu8//z5nWUUqallIeklM9KKX8Xo9DtQ5zD11Abz3OMlDINPAvcpu4zwyu3YYQhLjiEEK0YYb7hhd7LWeIoxh+m/TVtx6i6vSBfUwAhxHKgh/PkdTXbiT4PvA14vZTyaMkhzwIZil/HjcBKzpPX0cE1lmOb+f958TpWIAA0cA5fQx22XRg+ixFmeAZ4CvgwRsL7ywu5qXohhPgL4HsYodpBjJacHPAvC7kvP5hfAOzfzNeY+ZUpKeUJIcRfAr8vhDiIYUz/GCOndM853qpnql2j+e+TwHcwviisAz4DHALuO7c79czfAncBbwGiQgiVA5uVUiaklLNCiC8BnxVCTGEMRP4b4Akp5ZMLs2XXVL1GIcQ68/EfAJMYOc/PAY9IKV9ciA27RQjxaeCHGEVAbRjXcytw+zl9DRe6xPhi/Qf8GoZxSWGETbYv9J7qeG3fxDAcKYyerG8C6xZ6Xz6v6VaMcvfSf18xHxfAH2EYliRGtd8lC73vel0j0IRhJMcwWgGOAV8A+hd63y6ur9y1SeAXbMc0YhigKYxe3e8CAwu993pdI7AC+DGG4UwCBzG+BLUv9N5dXOOXzPdfynw/PgjsONevoR5JptFoNBqNS3TOU6PRaDQal2jjqdFoNBqNS7Tx1Gg0Go3GJdp4ajQajUbjEm08NRqNRqNxiTaeGo1Go9G4RBtPjUaj0Whcoo2nRqPRaDQu0cZTo9FoNBqXaOOp0ZwjhBAPmxq4RT+Xu32uWejnr4W5P2n+21anNb9iW/Ot9VhTc/GgjadGUwfMD+J7fCzx08An6rSdilQxkufk+cshhPicEOK7Dg79IsZw9ZdrHeiQD5nraTSu0VNVNJpFgJRy6iJ+/muB7zs4Li6lrNtMRinlLDArhKjXkpqLCO15ajRnASFEixDia0KImBBiWAjxv2ocXy6M+9dCiM8IIaaEECNCiE+VnNMmhPiGEGLefI6PVAu/CiG+AtwCfMgWrlxd+vzmz38jhPhLIcS0EGJUCPGL5jV9WQgRFUIcEkK8sWT9gBDid4UQR4UQCSHEC0KIt1e55ogQIgPcAPypuR/HY6OEEDcJITJCiEbbfavNdVaZt98uhHjJ3M+kEOJBIUSL0+fQaCqhjadGc3b4cwxD9RZgJ8a4r6tcrvEejJFK24GPAX8ghNhhe/yzwI3ATwE7gNfWeI4PYQwEVuHPpcDJKs89geEV/g3wd8C3gcfN57gf+GchRLPtnN8F3g18ANiCMSfy60KIWyo8R9bcPxgDmZcCd1TZfynbgFellEnbfVcC01LK40KIpRgzZP8J2IzxGnwXY3ycRuMLHbbVaOqMOVT6/cC7pJQPmfe9B2O2qRtelFL+ofnzQSHErwG3AQ8IIdowDNxdtud4L8Yc1bJIY1BwGmfhzxeklH9irvtp4HeACSnlF837/gj4FYxhyk8KIRqAjwNvkFI+Ya5xRAhxE/DLGDMkS/eTF0IMApNSyhdq/jbO5Arg+ZL7tgFqraUYn3HflVIeN+97ycPzaDRnoI2nRlN/1gERjCHngJFTFELsd7nOiyW3h4El5s9rgTDwlO05Zj08R83nllLmhBCTFBueUfN/tZ/1QDOGYbevE+FMA2fnSgrGzi3bgLvLrLfH/PkF4CHgJSHEfRje8r9JKac9Pp9GY6GNp0azeMmU3Jacu1RLuee27pNSStNIqv20mv/fCQyVnJuq8jzb8GA8hRBB4DLONMxXAd8x95gzw9w3YITOfx0jt7pdSnnU7XNqNHZ0zlOjqT+HMQzNdnWHEKILuKSOz3HEfI7X2J6jw8FzpIFgHfeheAXDSK6UUh4q+VcprwqwlYKn6IaNQCO2MLUQ4npgmX09afATKeUnMbzSNPA2D8+n0RShPU+Nps5IKWNCiC8Bf26GO8eAPwXydXyOqBDiq+ZzTJnP8Yfmc8gqpx4DtptVtjFgSkrpe1/mfv4C+JwQIgA8BnRgFATNSSm/WuHUALDRzH3Om+0jTthm/v/rQoi/xggb/7V5XwRACLEdI0d8P8bvZzvQB7zq5to0mnJoz1OjOTv8FvAo8D3gQQxj8mydn+OjGNWz/2U+x08wDEOyyjl/AeQwPMVxYGUd9/MJ4I8xqm5fBe7FCONWC5H+PvALGKHe33fxXNuA+zByvy9hfDn5JDAH/IZ5zBxwM/AD4ADwJ8D/klL+0MXzaDRlEVJW+5Kq0WjOF8z+xSEMA/Glhd5PPRFCPAzskVJ+2Lx9H/C0lNKNwa20tgTeJqW8x+9amosH7XlqNOcpQogrhRA/J4RYJ4S4CviG+dB/LOS+ziK/aopObMVoU/HVdiKE+HshRKw+W9NcbGjPU6M5TxFCXAn8I0bxTBojLPxRKeUF18sohFgGNJk308BxYIuU8hUfay4B2s2bw1LKeX+71FxMaOOp0Wg0Go1LdNhWo9FoNBqXaOOp0Wg0Go1LtPHUaDQajcYl2nhqNBqNRuMSbTw1Go1Go3GJNp4ajUaj0bhEG0+NRqPRaFyijadGo9FoNC7RxlOj0Wg0Gpdo46nRaDQajUu08dRoNBqNxiX/P88jQZ5efGslAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt\n", "import numpy as np\n", "import scipy\n", "from qutip import sigmaz, basis, sigmax, fidelity\n", "from qutip_qip.operations import hadamard_transform\n", "from qutip_qip.pulse import Pulse\n", "from qutip_qip.device import LinearSpinChain\n", "from qutip_qip.circuit import QubitCircuit\n", "pi = np.pi\n", "num_samples = 500\n", "amp = 0.1\n", "f = 0.5\n", "t2 = 10 / f\n", "\n", "# Define a processor.\n", "proc = LinearSpinChain(\n", " num_qubits=1, sx=amp/2, t2=t2)\n", "ham_idle = 2*pi * sigmaz()/2 * f\n", "resonant_sx = 2*pi * sigmax() - \\\n", " ham_idle / (amp/2)\n", "proc.add_drift(ham_idle, targets=0)\n", "proc.add_control(\n", " resonant_sx, targets=0, label=\"sx0\")\n", "\n", "# Define a Ramsey experiment.\n", "def ramsey(t, proc):\n", " qc = QubitCircuit(1)\n", " qc.add_gate(\"RX\", 0, arg_value=pi/2)\n", " qc.add_gate(\"IDLE\", 0, arg_value=t)\n", " qc.add_gate(\"RX\", 0, arg_value=pi/2)\n", " proc.load_circuit(qc)\n", " result = proc.run_state(\n", " init_state=basis(2, 0),\n", " e_ops = sigmaz()\n", " )\n", " return result.expect[0][-1]\n", "\n", "idle_tlist = np.linspace(0., 30., num_samples)\n", "measurements = np.asarray([ramsey(t, proc) for t in idle_tlist])\n", "\n", "rx_gate_time = 1/4/amp # pi/2\n", "total_time = 2*rx_gate_time + idle_tlist[-1]\n", "tlist = np.linspace(0., total_time, num_samples)\n", "\n", "peak_ind = scipy.signal.find_peaks(measurements)[0]\n", "decay_func = lambda t, t2, f0: f0 * np.exp(-1./t2 * t)\n", "(t2_fit, f0_fit), _ = scipy.optimize.curve_fit(decay_func, idle_tlist[peak_ind], measurements[peak_ind])\n", "print(\"T2:\", t2)\n", "print(\"Fitted T2:\", t2_fit)\n", "\n", "fig, ax = plt.subplots(figsize = (5, 3), dpi=100)\n", "ax.plot(idle_tlist[:], measurements[:], '-', label=\"Simulation\", color=\"slategray\")\n", "ax.plot(idle_tlist, decay_func(idle_tlist, t2_fit, f0_fit), '--', label=\"Theory\", color=\"slategray\")\n", "ax.set_xlabel(r\"Idling time $t$ [$\\mu$s]\")\n", "ax.set_ylabel(\"Ramsey signal\", labelpad=2)\n", "ax.set_ylim((ax.get_ylim()[0], ax.get_ylim()[1]))\n", "ax.set_position([0.18, 0.2, 0.75, 0.75])\n", "ax.grid()" ] }, { "cell_type": "markdown", "id": "0b497a70", "metadata": {}, "source": [ "In the above block, we use the linear spin chain processor just for its compiler and do not use any of its default Hamiltonians.\n", "Instead, we define an always-on drift Hamiltonian $\\sigma^z$ with frequency $f=0.5$MHz, an on-resonant $\\sigma^x$ drive with an amplitude of $0.1/2$MHz and the coherence time $T_2=10/f$.\n", "For different idling time $t$, we record the expectation value with respect to the observable $\\sigma^z$ as the solid curve.\n", "As expected, the envelope follows an exponential decay characterized by $T_2$ (dashed curve).\n", "Notice that, because $\\pi/2$-pulses are simulated as a physical process, the fitted decay does not start from 1.\n", "This demonstrates a way to include state preparation error into the simulation." ] }, { "cell_type": "code", "execution_count": 2, "id": "d7975d93", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "qutip-qip version: 0.2.0\n" ] }, { "data": { "text/html": [ "
SoftwareVersion
QuTiP4.6.3
Numpy1.22.2
SciPy1.8.0
matplotlib3.5.1
Cython0.29.27
Number of CPUs12
BLAS InfoOPENBLAS
IPython8.0.1
Python3.9.0 | packaged by conda-forge | (default, Nov 26 2020, 07:53:15) [MSC v.1916 64 bit (AMD64)]
OSnt [win32]
Sun Feb 13 12:25:17 2022 W. Europe Standard Time
" ], "text/plain": [ "" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "import qutip_qip\n", "print(\"qutip-qip version:\", qutip_qip.version.version)\n", "from qutip.ipynbtools import version_table\n", "version_table()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.9.0" } }, "nbformat": 4, "nbformat_minor": 5 }