Introduction
Using data from dogesr
(Merelo-Guervós 2022), we will, in this
vignette, analyze the change of the social network before and after the
Serrata (Puga and Trefler 2014),
revealed as a watershed event in the analysis performed in (Merelo 2023)
Set up
The main events that led to the Serrata happened during the
government of the doge Pietro Gradenigo. Let’s check the effect of that
event in the social netWe will then create two networks, one including
that doge, the next one with all the others
library(dogesr)
library(igraph)
pre.serrata <- marriage.graph.slice(to="Pietro Gradenigo")
post.serrata <- marriage.graph.slice(from="Marino Zorzi")
Social networks pre and post-Serrata
This will call the function that creates, from raw data, the
igraph
data structure with which we can work. We still need
to simplify them to be able to work with it.
weightify <- function( a.graph ) {
other.graph <- a.graph
E(other.graph)$weight <- 1
other.graph <- simplify(other.graph, edge.attr.comb=list(weight="sum"))
return(other.graph)
}
pre.serrata <- weightify(pre.serrata)
post.serrata <- weightify(post.serrata)
V(pre.serrata)$EV <- unname(unlist(eigen_centrality(pre.serrata)$vector))
V(post.serrata)$EV <- unname(unlist(eigen_centrality(post.serrata)$vector))
plot(pre.serrata,vertex.size=V(pre.serrata)$EV, layout=layout_nicely, vertex.label.cex=0.7)
![](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAASAAAAEgCAIAAACb4TnXAAAACXBIWXMAAA7DAAAOwwHHb6hkAAAgAElEQVR4nO3dd0AT5/8H8CeDPcLeQzbKFlFwIRQV11er1lntUGvF3eKv2jqwxVpxQdVqqatuVBxYtQoiAioIIqBsgbA3SYBA9v3+OHtNE0CrnIx+Xn+Ru8vUd+5y93k+DwXDMAQAIAe1t18AAAMZBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEkHAACARBAwAEtF7+wWAHsblCjdvTlZWphsYqFZVtd6/X/H06aLy8pacnKZJk6yIzdrahMOHn8nN/Uz6vp0uBO8CAjagSCTY+PGXPvnEaflyN4QQhqFp065IJNj9+xUnT76QDpi6usLjxwtk7t7pQvAu4BBxQLl8ubCqqm3ZMlf8JoWCwsP9JRLM3d1AZksMQ5mZ9W+yELwLCNiAkpJS7eqqT6VS8JtJSZVlZS2PHlUjhNrahLt3p+npbZ46dSNC6OLFfD+/KIRQSQnH2PhwVFT+xo2JWlo/jxsX9d13SVOmXDl69LnMg7e1CYcMOfF+31C/BwEbUDo6RHQ6hbhpa6s9YcIlPl+MEMIwbN06Dzo9PjlZwOFwWloENBoVIZSRUefnZz5zpv2nnzq7uRmoqNB37Bhz5szklSvjmEyO9IPDAeRbgIANKO7uBs+fN2LYq5vGxmqamooWFhoIIQ0NRQUF+p07V2xsXB48eGBtrYRvM3WqTXu7yNn5REsLn0ajUCgIIYRhiEajKCrSpB8cDiDfAgRsQPn0U2dFRdpvv2UTSyQSDCGEYdhfqaM2NAhra43270/Cbx87lh0T8/L48cB16+4jhIRCbN++9Fmzrm/fPsrERP3Mmdxp065u3Ji4bl18YmLFuHEXfHx8uFx+ePjTsLAnM2Zca20VREcX+vtf/OWXzHHjom7dKkEIyS/h88XSd3nPH0svgoC9P9HRhRTKni++uBsRkbFr15OhQ081NLRLb5CeXrt9+yOhUPLWv3aUlGipqQtzcho3bUrav/9pSMij5cvdLCw0ExIqiovZZWUtP/6YyuHwJk8ePGbMYJFI8vvvqWpqCoqK9LY2YVCQO5vNFwrFmpqKU6ZY//ZbdmkpZ+hQw7y8pp07x+7cOVJdvZ5OFzc2NoaHP21pEVhaapaUsE+dyvHyMqqoaFmxwn3LFp8TJ14ghOSX7N+fLn2XHvk8+wU4Tf/+zJplr6+vunDhYF9fc4TQxImD8F9BhMGDdf39L27cOOJdfu1oaChGRPgjhAICLjo66tJoFHX1iJCQkePGme/ZkyaRYBoaSmZmGgEBtuPGlevqNgwZoq+iQps4cRBC6Nix50eOjF+61BUhdPNmSUTE0xUr3PX1VUpLS168eMHhaNDpSjk5+XPn3ti718/amjF3riNCqKqqTVmZTqEgLS0lNpuPEKLRqDJLUlNrpO/y3wF7sPeKRvv7DIREgunoKEuvVVNTUFNTQD30a2fevMEHD34QFOROpVK2bPE5fXqyi4u+paVmYOAghJBIJKFQaH5+fvn5+WKxBL+LWPzqOBLDEIvF09NTbWpq4nA45eXlvr6+M2eO5vMlGIbMzDSOHMlECPF4IvwU5V/3wmReA7Gkq7sMeLAHe9+iogrS0mqzshoCA61MTNQPH8708DCMji48fXoysU1iYoWfX1Ramp+zs/vhw5kCgfjRo+rTpydfv/4yKqrAyUmXxxOFhfmGhqZYWTHS02snTrT63/9sZJ7o88+dZZYsXepSV9c+aVL00qV3tLWVqVTU0CCi0x1aW3MPHbrn4uKQk9N47drLjg5RVlbDkCE6o0djZ86k1taKdHWdGQxGQkIFhmF//lm6ceOIOXNifHzO+vqah4aOjo4urK9vZzI5CQkVpaWc2lpucnKlzBKZu5D+KfcZFPlvHUAeY+PDFy5M9fU1z8ioKyhgDR9ulJXVMG2azbBhp+/e/cjQUNXY+DCT+UVVVZmNzbnp09O9vTcJBBIHB+0dO1KWL3fz87P43/+uFhUtxTBs//6nQqF448YRjY0dgwZFVlQs19ZWln/GgoJmZ+eTQuFX+M0HDyr8/KIkkmA+X0xE99ixgJ9+upWU1DphgmNMzMuZM+3mzNHJycmPi0PDhllnZNQTAcZfHoWCus82IMAhYu8YOtRw6lRrc3MNLldYV8dVVKSJRBJiLZvNVlFRDg8PT02t+fjjIXPnOmZnf7pypQeNRtHXV6FQEJVKefasDj+NrqenYmmpWVLC6frZ/ubra66pqYT+edbhwoWiBQvGNjTwjI2fr1mjGBr6OC+vurDQwNLSeMkS1+3bRy1Y8AeLxSMe5MCBZ6qq9M8+c5ZfBWRAwN4riQQjfudoaCgeOZJVXt5iZqbB4fCJQ4n6+vrq6mpFRaVBgwZ189Nl6FDDpKRKhBCGIaFQ4uio84avAb/SJRNdRUW6sbH29u0hwWsXODvrWlg45eSwugrw22X7vwkC9v7cuFHc3MyLisrHT6whhOzttS9cyA8PfzpkiO7hw5mpqTWNjR1Hjz4QiSw5HP6tWyUbN454+LDKx+dsSMij4cON7t0rf/mS/fx5I0Jo1SoPExP1yMjssLAnR46Mx8+OyJNIMInkH78D8L9loisWi1taWjdt2rR5+z42m2trq9FNgN862/9B8Busb0lJSVFWVnZ3d++RR2trE+7fn75168PDh8d/8YUrlUpJSKjw84u6eXOmm5vBnDkxEgnm62u+efOwCxfuffNNSXh4QGMjT02txcpK5O09auPGR25uBiwWz8vLyN/fIjW1ZvTo81evTp8wYdD69felV/XIqx2QIGB9SFlZWV5e3vjx42k02uu37iFcLvfBgwd8vta+fWVJSfPxhfn5+aWlpb6+vqqqqu/tlQxIcIhIiocPq3R1D06aFN1VZbq8jo6OzMzM4cOH02i04mL23btM8l8m4nA49+/fd3BwMDQcJF3B5OjoaGdnd//+/dbW1vfwMgYwCBgpRo0ydXLSmzHDtqvKdHlpaWn29vY6OjoIodu3Sw8ceEb2i6yvr09ISHBzc7OxscnPb/r8c5e0tFpira2trbOzc0JCAovFIvuVDGAQMLIQRRudVqbLyMvLE4lEjo6vyog8PGTHR3aDyeTcucNE/xyvFRqasnZt/Jo18ZaWkW1tQvnixqqqqvXrb169SnF3j7G0jLS11WaxeO7uBtJbWlpaenp6JiUlNTQ0EHckCibf/BX+p2GAHOPGXfj445t796aNG3chLOwJhmE8nmj//vRdu1KnT7/a0sK/ebPYx+fs/v3pZmZblZSmhIU9mDHj6pEjmRiGJSdXjhx5dvXqezY2vx069AzDMD5ftGVL8vHjz4OCYq9ffyn9RJGRWQEBF/G/2WwehmGlpWw7u6P4kl27UgsLm4lVuJKSkqNHo21sIvHlu3alZmXVa2hEMJkcmS0xDKurq7t27VplZSV+s61NoKERweOJ5N9ydXVbT3xyAwrswUjk4qI3bZoNUZkuU1E+erRZSQl77VrPmTMlYrH3ihXDDh4MiIx8NdJEIsF+/tn//v2533zzoLGxo5tru0OHGuJ/EBWMDIZSWVnL9esvEUIrVrirqSkQq2pra7OysnJzcydO9K2oaLt27WVmZv2KFe56eipKSrT/+78H8mWQBgYGY8aMefr0aXFxMZ/PJwomZeTnN3/11X1yPsh+DAJGIgZDyc5OOzjYy9RUPSLiqXxZhoICjUJBq1Yt1dFhqKuramkpVVe34ffV01NFCJmba9jaaldXt8lf2xWJJD/+mHrhQv7Jky/wu+AVjAghbW3lLVt8Zs68vnr1PQqFYmKijq9is9mDLEynT9/T2jrom29SvvzSfdas6+PGReHbCASSvLymn35K9fOL2rr1YXj4Uy+v099+m7xo0a1bt2oPH+aPHBmqp/e5SCTCn662lrtt28Nr114uWnQLIZSQUJGb23TjRrFAIN669eGJEy9WroyLiSl+/x97nwIBI4t8ZXo3RehYFxdLxGJMLJY4OurIX9s9duy5ggJ13jxHYgAIUQaFENq82TshYe6ff5aOHn2urU2Ir2KxWBaDPhs/3nvePKf16z2XLnVJSJhLpVLwbRQUqMbG6ps2jVBSoquq0qdNs1m5cujPPz999qzO0pJRVMRZvFg/IKCDuITA5Qrd3AymTLHOzm6oq2v38DAwMVGfNs0GCqmkQTU9KRITK6Ur011c9Net8+Rw+NIV5XfvMlksXkFBc1paA4slzMpqqK3lcjj8goJmS0tNCgWdOpVTXd124sQkRUXaqlUexcXsyMhsFouH120kJFTMmmWPEFJR+fsfES+DKixk2dtrjxljlpKycOjQ02fP5i5f7oYQlp2dbW8/cf78EQghTU0le3tthJCGhgKLxT97Npd4EIlEoqhIMzfXoNEoJibqGIYwDNPXV9m9O0z6PZqbazx6VC1fSPnsWR1+1ErsbD09O6lC/o+AgJFi7FizpqZVMgvV1RUePvx7GOXkydZtbWsRQqqq3MTEkW5u+m5u+u3t6/C1MTEfSt9XUZF26FCA9BJTU/X4+PLZs+1FIon03hIhlJBQ0dTU4eNjoqurYmXF0NZWFolEIpHYzc3NwYF54sQLf3+LO3eYiYmVS5e6IETBt0EI4e1xqFRqUlIlnU5taRHgCzvdwR45ktXaKiAKKel0Kn53fGf71VfDoJAKQcD6r+Bgr3nz/iBGdjGZHCazBa9gtLJirF59b+HCIQ0N7c7OerNn2x8+fIfLFefmSjZtGjFv3g0npxMODjplZS3PntVzOHwNDUVzcw0Oh19e3rJyZRyPJxIKJcXF7Kioglmz7Kqq2nbuTMVrIF1c9PCCydhYpr299oYND9TUFPBCyi1bfJqbO/btSw8KcpfZ2fb2R9WboFSq95WXl1dXV3t7e5P0+Pn5+VVVVX5+flQq/OR+32AP1i99//3jnJxGLy8jBQVaeXnL/v3piYnz3d0N5DvL19fXFxUVBQQESKcLetC/NxCwHvPgQcXUqVemTbPx9DSqreUqK9O++85bWfn1nzCGYXg7QiaTU1DAwvvPyEhPr715s+Tbb70VFKgIITMz9a1bffBVkyZFBwV5jB5tihCSaZXD4/FSU1O9vLxUVFSkl/dsC9Gmpo7vv39sZ6dNpVLy85u//XaEkZHav3qEbgLf778Leu0S90Dk4HDsjz+K8b83bUp0czspEIhfey8mk5mSkoL9syZDhkz9hFD46mGPH39uYfFra6sAwzCJBEtIKCfuIhaL4+Pj8/Ly5B9NZst3IRSKnZ1PxMWV4TeTkyutrCK5XMG/fRyZ8pE3XNX3wUF5T6L83TMKhYaOZrP5v//++h6A2F97MKImQ55M/QSdTkUIVVe3ffXV/cjICerqCkjqQnNtLXflyus+Psu2by90dHSUbwNKbPnuLlzI5/PFH3zwakjYqFGm+vqqRD3KG+qmi1Z/bycMh4hkoVIpXl5G6em1U6daS7eOio4uPHQoc/Zs+4sXC/7v/7wmT7aOjDwaHd0cEqJBXH0WCMRddZWRXpWWVvvhh3bEISVxoZnLFfL5zCdPfrewcK2rayfagDo46Bw5kjl5srX0Jel3lJFRh19PIzg4aGdk1IWHP715s+TmzZmbNydXVLSePz9Vus2OTIesDz+0w1vxyH84+HeBRBLcI6/2/YM9GIkUFWmKijSZigf5rreXL5czmcXTp1sSNRndFEMQq5yd9dLTa7duHYkvFwjE6K9dqLm5xujR/rdvJxkY6IpEEvk2oOifO9t3oaxM53KF0ks6OkSamko+Pib4J/DBB5b4crwUU0Gh7ebNxytWHCV6Bu/b50cEXv7D6cHvgl4BASNRZmb98OFGMq2j5P+7e3jMXr16IZfLJWoyuukqg6+qq2vfvDnZzEyjqakDX376dC7660LzkSNZVVXcwEAf6V466J+NQXvq6szYsWZ5eU1EJQd+RDd2rJn8lngp5pQpZgzGMX9/ZekOWeivwJP6XdArIGA9iaioQAgdPpypqqowf/7gTltHIan/7ubmmrW16kVFRURNxmsbzqxcGaeursjjiRMTKyIiMmbPjmls7EhIqMAvNMv00pFvA0ps+e5vOTDQytfXPDz8KX7z1KmcwYN158xxUFamt7UJEEJsNg9/U3gppq2t7YEDhx0dA2UeRybwZHwX9Ar4DdZj4uLKmEzO+fN5xcXspqYOkUhy794cGo0iU/Hg6qov8989ONhr7twb27a9cHJqx2sy5CsPifoJfJWbm4GXlxHRcGbt2qH4a8CwV79VAgOtEELr1nniN/GDz+Bgr+BgL4SQkZEaseW7O3du6okTLyZMuCSRYNOn20ZHT0cIubrqa2srBwXFGRurcTj8+vp2ormvhYXo/HmjyMhsojqECHxrq0Dmw8nPb8ZXTZ5s3VMv+H2CSo6+Ijs7G8MwNze33n4hbyktrXbx4luurvqnTk1WUupy+HZjY2NWVtYHH3zwPl9bL4JDxL7Czs6utLRUKBS+ftM+ycvLKDPzk0WLhpw8+SI+vryrzXg8nsxV74ENDhH7ChUVFUNDw7KyMltb295+LV3icoWbNycrK9MNDFSrqlrv3694+nQRsVZJiTZ1aud96okilY6ODmXl7kav9PvSjX+CPRiJ7txhqqtHrF0bL73Q2/usqemRtLRa+UY0dnZ2hYWFffagXSLBxo+/5Oios3PnmPXrPXfvHmdsrIbPoPlasbFle/akoTfYgw2wmaAhYP+CTKsmmbXy7ZYmThzk52dx8uQLYuOsrIaODpGnp6GXl5H8/yQ9PT0lJaWamppOn72mhivzx3t2+XJhVVXbsmWu+E0KBYWH+79hwIgilY6Oju4D1t9LN2TAIeKbYjI5p07lFBYuQQiZmanX1LTZ2f2jgoGYn1J64YgRxlVVradO5QQFuSOEoqML58xxSE+vRX/9T8JnuyTY2dkVFRWZmJjIPHt+fvP27Y/On59K/EHGe+xeSkq1q6s+ftkKIZSUVCkQSKqr28rLW4iajI0bR0iXrYhEkrCwNGtrxsOHVfi9Wls77txhKinV4/Ucd+8yuynd2Ls3XV1d4eLFgsBAq48+sj9+/EW/mzMJ9mBvSr5Vk8wGamoKioqoqKhIeiGFglat8jh48BlCiMsVKivTifp6oiBQulYwO1uwcePGr7/+P5m+MURLmV7sLdPRIaLT/77oa2urPWHCJT5fLF2TIVO2ItM4pLGxcdWqqIqKhm7mdyZKN16+ZN+4Ubx8uZu3tzGFgqKji/pjqw8I2JuSb9XE54vDw5+GhT2ZMeNaa6tAIBA0Nzd98sknFy7kq6qG19e3c7nCs2fzxo0zZzI5GzYkTpt2RV//71bvXdUHlZeX37jBkvnPRLSU6cXeMu7uBs+fNxK/EI2N1TQ1FS0sNKRrMmTKVhISKqysGAghJSVqR0fH0aNH8/Pb6fRsqdZaXZZuGBqqNjV1PH/eqKam6Oqq30/nTIKA/QsyrZpk+hwqKipqaGhu3bp13jxHb2+TgoJmoVAyZIju1asvvbyMXr5sHjbMeP36+I4OEfGA8vVBLS3CkpISL6/pr/3P9B7+w33//eO5c2/s2ZMWEZHx9dcJK1bEikTYb7/9XSkv/wNMpmzF1FT9zp3i7OzshITE9nbe/PnzJ00a2dBggxDi8UT37pUTRf3ypRsaGoobNng1NLSvXOk+YcKgfjpnEvwGe1PyrZpSU2v27vWztmYQRbqqqqqBgYEIoXXrPH/5JTMw0GrQIMazZ3Vjxpjt2vXk+nXXmzeLGxo6ZHrUEPD/ZEpKSp6eRjJ9Y/Lzm/GWMu+zt4z8sM6dO8ds3py8aVOSgYEqh8NfvtzNwkLzxIkXRE0GXrYiFmODB+vs2pXo76+0dWtpUxPbwkJfR4eCYVqRkfMnTDg1bNjJgADr0NDRmzYN37gxqdPSjeHDjbdufejhYaCpqTRjhq18dUuPv18yQCXHm4qMzHZx0cOLxMeNiwoKck9KqlRRoYeF+fJ4ooyM+pEjTfApjJWUaBIJ5uR0wt/f8o8/iidNsqqr41pYaO7b52dre9TamlFQ0BwT82FLiwCfqqu1VbBmTXxq6sLLlwuPHMlKTp6vo6MsMwEXjyfy9j67eLHTl1+6jRx5bvFip6Ag96+/TiB1ki6RSIIPPDtx4kVIyKOcnM/wgWfdEAgEubmN3357f/lydVVVVTs7OwsLC5leIOXl5S9fvvT398cwlJhYIXOah5Cf31xQ0GxursFi8ePjy4KDvTqdhLqv65Vhnv3R3btMT89T+/alb9qUuHJlnFgsqaxsHTnyrLf3mW++eSAUilNSqun0vTduvGocf+BARl5eE4ZhfL4oKCj211+zfvop9d69sl59E2+jqqpVS+vnP/8sxW/W1LRt3Zp89WrRxx/fxDDs8uUCP7+oQ4ee+fpemDt3m6oS+vjjQ/b2R6Kinnez5Q8/XKqpqUlIKKdQdnf1vF98cSci4imbzcvPbzp2LFskkryXt9vDIGDgNaZNu/LZZ7eJmy9fsqKjCwUCsavrydpablkZx9b2N4kEi4sr8/QM9/Dw+PPP3MDAy91vOXny+djYWAzDGIyfu3reqqrWn35KXbEi9rffsvn8Tuaa6BfgNxjozpkzuenptbm5n+M3BQKxTENf6TM02tomsbEZqamvLpR3s6VAQJNIJDU1Nd2M9TIxUf/mm+Hv4T2SCgIGulRX1752bfyxY4FaWq/GFJ8+ncvlCqUb+lIoxDRor37ME6dhZFr/ymw5ZMiQFy9eDPgzABAw0CV8WGdVVeuhQ89EIiwpqdLLy8jNTV9+eNupUznLlt1RUVGoreU6OeniLX4jI7OKizlNTbxOB8LR6drZ2dzXjvXq7x3d4Cwi6BnTpl1NTKyoqlqBn2nMympYvPiWpaWmTJN9adXV1c+fP58wYQKl264AHA6fwei8LUc3q/oIuNA8wMlcLKZS9yQnV5HxRCNGGNvYaJ069apNHV51Scyj2ykTExMajVZRUdHNNv29oxscIg5wXfUA7nF41eWePelBQe541SUxrlm+C92ZM7l4fXBmZmFyQuD/Zh4bP96TWFtX137lSqGOjnJrq9DOTqt/d3Tr1XOYgHSd9gAmQ2jo4/Z2oY7Ogbi4smPHsktK2Hv2pM2YcRXDsD170nbuTMEwrKGhXU0tvLm5Iyen0cbmN4kES0xM0tWduXz5Bem1ly4VfPLJLYFAjF9IxE/lS5/lnz37Ov6k3Zzl7yPgEHGA67QHsHSZckzMHVPTQC+vY9Ktf9+uVF9Fhb5kicvBgxnl5a14jS9OvmySqA8eM2Z0YODn9vY20munTbNRVVVwdj7R0sJH/byjGwSsH2MyOXfuMOWX37xZQqXuWbbszo4dKbNnxzx5UvPll7HSPYDxMmVtbfTiRfWBA48a61Orq7nSY0beolQfbz+6YoX7H3+UDB9uhBCSSLDXdqHrdG16eu2hQwHHjweuW3cf9fOObvAbrB+LjS27eLFAfjaWKVOsjYzUPv54iK+v+eXLhTNnXpdIMOJicXs7Pz7+5cKFDAMD2h9/+FtbW2dlffn553ekdw7/dhrYe/fKz57N8/e3CAiwDApyDwy0Ki3lxMYyc3ObMjLq5Ot0790rJ+qD5dcWFbFiY8t8fEyCgtz7fUe33j5GBW8vPb22q9lYzMyO4POnREUVUCi7o6MLMQxraWnJyMgICjo5a9bpoKCbGIZ1dAjPncubPv2qs/MJDMPS0mrwB9y799XPJ4kEs7M72tZG1i+3AQ/2YP2P/Dj82lqu9EB9fKFYjGVnN6xeHaeqSg0OPv7116YMBo3FwhYsGKKkhN2+Xaqvf2jJEhcLC018Itm0tJrvvkt+8aJx9uyYc+em9MexIX1Rbye8jyotZeP1462tgsGDj2MYlpxcqaNzIDDw8rffJk6eHP3bb9nv4dk7deRIZljYEwzDkpIq8R2OTFkthmFmZkdWrYo7dOjZgwcVn376KYViaGp6QCyWiMUSmY1TUqo7rc0l7939p8BJjs4RbcaI3k+jRpk6OenNmGG7Y8eYM2cmr1wZx2SSNWqdePZOEePwickiZAbq4wtnz7YPCnIfO9bs4MGDf/xxw9xcm0qlyI/qJx62q+XgXUDAOke0GZMuFyDqEjAM0WgU/NQzqc/eKVNTdbx1LjFZhPz8EsQZPISQmpqajc3fPZhkNpauze10kgrwLmghISG9/Rr6EJFI8tNPT6qq2m7cKG5s7Fi82AkvF9i2bSRC6Pffc9hsPpPJ+eGHx1984RYYaCUQiLdvf8xkthw9mi0WYw4OOny++ODBZ4mJlbt3pzU3P1y06PvTp4Vcrig0NKW4mF1UxP7hh8ccDn/YMCP5+8o/u/w2CCFXV/3Dh7OSk6tKSzl1ddyxY80wDO3f/1QgEHO5wrKyFoFAcuzYcyqVMmqUqaqqAkLo4sWCK1eKJk2yNjRUZbF40hvPnm1/4EBGe7vI3l47PPzv5cS8XuCd9PIhah8j//MGkyoXGDfuwq5dqYWFzbt3P7GzO4oXK8jUKOzcmbJ9+6MLF/JcXE74+X2nrKisrb1fLJZUVbUaGBwSiSSVla1Dh57COqtvkH92+W164SNoeO4AABRMSURBVEMB7wAOEf9B/ucN+me5AIOhZGenHRzsZWqqHhHxVL5GAZ9mDm9LFh8fWlxaQadTHz5MVlWlKirSaDSKlpZSdXUb6qy+Qf7Z+2mvMkCAgP2D/M8bJFUuIL2ExeLp6anKVyHg08whhHg80aNH1QwGQ0lJicFg3L9/XyL5x5kD+fvKP3s/7VUGCHAd7B+Cg73mzftj6dI72trK+Fx4TGYLXi6grq6Yk9N47drLjg5RVlaDi4v+unWeiopUmetFxDRzvr7moaGj795lslg8JSUzDodbV5cXH58vFCo3N/Oo1D2ffeZcU8P19Dzt42OC31f62dvaBOvXx4eGjnnt9ah+Me7wv6u3j1EHiD//LFVTC1+z5p70whEjzpiYHH7ypAa/mNbQ0BATE5OXl4dhmLHxLwkJ5WVlnG+/TRoy5DhxF+KyW1ubQEMjgsd71ewlLa0mJOShQCDu9NnZbB5Zbwy8GzhE7BlvMpGKnp5eQEBAZWVlamoqjUZFCN2/X3HzZkl5eQtx8ElcdlNTU1BTU7hy5UpxcTFCaPBg3b170zudyqSnxh0SpcO3b5cqKOxbuPBmQ0M7QighocLSMvLcubzXPkJ5ecvt26Vv+HTyszcNSBCwHtPNkF4iAyoqKn5+fgghPp/f0NCkqkovL+esWTO0o4P3zTfx4eEPFy26sm3b3fT09KSkpI6OjgULFqxatWnbtoexsWX4xV/pmSLwoSXEJBLviLi6PWmSlZUVY+5cR7yT/rhx5gsXDl6wYPBrH+H+/YqwsCdv+HQDbB6wrkDAesybTKQiFmMHDmRWVmq3tIjmzAnbsuWImhqqqCiys4tMTy8eOlTs7q4SEVFAo6lbWFhLJLSpU8NYrEB7e50pU6yFQkl9fXs305G8I+mr2zQaRXq0/xvWIrq7G7z50/WLAf/vDgLWk+bPH1xXx713rzwqKn/+fEdiOZGBiIinRkZqs2bZa2oqmZsLZs2y1dPTOH78ExcXk5oa+tixY7/+ehKVSpFItM+fr8Qwyty5E7lcyt27zLo6LoVCkeku2CPjDkUiyY8/pl64kH/y5Ivut5Tfee7dm/7rr1kffHBx9+40/NpDW5tw9+40P7+oq1eL0D9Hdra2Cs6cyZ027erGjYnr1sUTXzq1tdxt2x5eu/Zy0aJbb/82+ioIWE/qakgv+isDSUmVBgaqCCEajRIUNH/OnDlaWip0Os3d3aCurg39dUnA0VEnNbVGRYX+0UcOy5a52dtrm5lpSCQY9o+xhj0z7lBmCq9uyOw8pefvYjJLlyz5qKGhAcOw4GCv7dtHHTv2HP01spOYgEZ6JjHiS0dmSrG3fyd9EgSsx3QzpBf9lQFra60TJ17cvl3a2Nhx+3Yrm83Dl8+f76igQIuMzA4Le6KiQldTU1BQoDY1dcTGMgcN0jx+/EV4+FM6nRIZmZ2cXCkz7pAYkvh2L7vTa+sMhhJeoIjDDxdldp7S83e1tRX9+eefZWVlGhqKFArS1VVuaupACElfdl+50kN6JjH015fOwC4yhoD1DHxIb1xcmZUVQ3pIb0ZGXUZGHZGBTZtG1NS0BQcnfPih4YkT/s+e1RcXs8vKWhITq6hUyuTJViNGGHO5wlu3Svbv9/f2Nt6xI+XRo+qCgs99fEzEYmz0aNO5cx3r6oIGDWIEB3u9fLnUyEht3DhzDAt+61G9nV5b9/ExSU6uxP9ms/n4OU8CvvOUnr8rJGRBTEyMp6cn/n2BYa++UGQuu8s8Nb7NwC4yhsajvePatWsTJ07sfjrw96O2ljtv3h+2tlra2srPntUdPTpx0CBGczNvxYpYXV0VKyuGQCD++uthysr0qKh86WmWrl2bMXXqFXz+Ljc32pQplrdvt+3dm56cPP/atZe7dqU+efIxhUKZMydGIsHwy+6Rkdnbtj2Mj5/r4qKXkFCBz95EpVI2bHiwZIlLQkKFk5Pujh1jevsj6UkQsF7AZrMfP348adKk3n4h70R6/q4TJxK3bvW1t+98pq//MjhE7AUNDQ36+vq9/Sre1f796WVlLTY2WmZm6lpaYnPzf3GO/r8DAtYLBkbAtm0b2dEh2rQpKTa2ePx4XRWVPt0jvrdAsW/PaG7mbdv2UFNTUV9ftaiIVVrKuXVrVlcbNzY2uru7v8+XRwZi/q7y8vLKSmFvv5w+CvZgPUAolAQEXBw71mzHjjHr1nnu2+enoNDlB9vS0kKj0VRVVd/nK3xHd+4w1dUj1q6Nl17o7X3W1PRIWlotm83W1tYmlqen127f/kgolHRTbfgfKUREELAecfp0Dp8v/ugjB/ymkhJt0ybvrjZ+8uSJmpra+3ppPaP7UmYWi6WlpUVsTNQld1Nt+B8pRERwiNgj0tPrnJ31pJd4ext3OqXI6dPP42JPmZlb79tnIjNXCJ8vPnw4UyAQP3pUffr05OvXX+Lzj/B4ovBw/66emsnkFBSwJk4cRIwKCwi46OioS6NRDh58FhIysrCQpaWldODAB2/4XthsfkjIQ1tbbSqVUljYvGWLj66uCkJoxAjjqqrWU6dygoLc0V+lzOnptQghmT0YPg4AIYRhKCoqv6aG++233jK7dLwQ0dd34J91hID1AD5fLF33gCPau0+bZjNoUGRFxfKhQw1DQh4t+8IiIGAcUXbk4KBz5Ejm5MnW+/enCwQSBwftM2dyT53K8fOzCAl5FBPzYffXUYju2cQ+Yd68wUuXuhQUNP/yS+aWLT4IocjI7Dd8I2Ix5u8fFRo6Gr9snZRU6e199vnzT/ECjk5nJ2pvb2ezxTt3PpVpe4oQSkys+PLLWHV1RWNjtcuXC729TWJiXs6e7bB5szcx81A3XythYb4y31D/5t+kr4BDxB7g5qafk9Mos7CrKUWOHDk8e/Ys+Zrd7quKuiLfXu7zz51ltlm61OUN38ilSwUsFp8oChkzxozBUMKrClEXpczJyckNDZxOiwnxakM1NQVHR83yctb27aPu358bGvq4qqqNKETspljxLSag6INgD9YDlixx2bcv/dSpnMWLnfAlsbFlHh6GSUmVX301jGinUVnZKn9fYgeFVxWFhfnyeKKMjHpd3e4mW5Dvnk3sE+QDKRJJZHYF1dWtS5eGjxrlYWho+dlnzgcOZJiba1y5UqSiQre01Ny27SGxO1JWpu/YkZqWVnvpUmFmZr2VFePHH1MqKlpzc5svXy4cPJjx4YdHKBQFY+P2hgbj9va2L7+MsrFRZbE6duy4OmYMQywWCQSS9evXFBZaBgffdXY2YTCUL14sWL/eEy9ETE2t2bvXz9qagZcaFxQ0418rFArl305A0TfBHqwHqKkpPHq04M4d5po18RERGYcOPXNw0F692sPERB2v35WZUgQhJF+zu3HjiIcPq3x8zoaEPBo+3Eh6Y3nyJfDdjAqT3xVER2fcvp2TkHDQzo5OjKBZv95TLMYkEkx6d6StrcTlCjZsGI7nNipqWkJChbIybf36IRxOe3Z27ejRAe7u41taUHu7KYuFrK3Nv/tugqam0t69lba2LlQqjU6nT506VU3NVE9P89NPnXV1lbZsSWKxeK8tVhwYDX9gD9YzTEzUz56dIrPw0KEA6ZtBQe74GQKE0Ny5r7IRHOwVHOyFL3z4cEGnG8tLSKiYNcsedd1eTpr8ruCLL0bHxTGzs9tfvCi6fLlm1SoPDMOGDTOaMGHQd98lSte2FxSw1NQU0tJq+XxxSQnHxkbr008ds7Nr3N3P8PmIx0OrVn2ooEDdsOEBhUJXUaFnZTUXFvJYLL6xsVpsbH1Li4BGU5ox48OwsLMvXjSFhz89cCBgzZp7ly4V4tXPMj2CIiOzu5nW6E3/MfoSCFi/hJfAz55t32l7OYSQRIJJJFhTEy8k5GFtLffkyRfFxezSUg4xw921a588elS1YcMDd3eLkydzGIwGS0t7c3N1Hk8cF8dctGhIY2PHiRPPhUIJg6EUGvrY0VFbU5P++PHjhobaUaMM9+2bOGFC9Lx5DklJlWFhvuPGmWdk1K9f75mUVDlihLFA8JWDw7GFCwd/8YWrsfFhhJCxsXp7u3D9ek98d4Svwl9qV18rioo0mW+o/ghaZ/dL8t2zMzMbfv01a/hwIzs77bY24fHjz+Pjy8+fz1u2zOWHH0bn5TU7OenFxpb98kuAg4NOfHz59esvDQ3VnJx0ly3zuHSp9PLlpvz8Ohsblqur6dGjhenpdRUVLfHx5du3j7p1q3TuXNs7d5hsdvucOfZ1dUrnzpUIhZhQKLG11UpNrTl6NJvJbFm0aIiXl9GjR9VVVW337pUtX+7m4KCTmlpz8OAzKytGQUGTk5Oe9Kre/gjfE6imH7COH3++d296Ts7f/RJTUmq8vY27uQuLxcrPz29oaLC3t29tZXz7bfLlyxNycnITEspmznS3s7Oj0d5mvovjx5+3tQl9fEy8vIze4u79GuzBBqzffnuuoaFI1JcghMzMNGprubt3p3E4gh07UmbOtIuOLly16h6XK9yw4YGBgaqrq4mxsem5c7V5eVU//BD97Fmhg0Mbh6O7enXOgQP/q6/v6Oa+dnbanU5VgRDy8DAcMcLY1FS9lz6J3gRnEQesTi9/yzTAkO9RdezYc01NtW3bpnt4UCkUSkBAwMKFPp02z5C/78C4ctWzIGC9IzQ0Ze3a+DVr4i0tI4kCv7dDNAyV0enlb5kGGPLXu4kWHaNGDff29jYxMUFdNM+Qvy9MVSEPAtYLmEzOqVM5ERH+P//sv3Kle01N27s8WlfTYS5Z4iIUSohGqPiWhw933gCD+Bs/P8nhcHg8AY2m8NdahLpunkH8PTCuXPUsOE3fCxgMpbKyluvXX06fbrtihXtrq6D77TEMY7FYOjqd/38dOtTw4sUC+eX45e8NGx6kp9fZ2GjR6ZRp02zEYsmGDQ/U1BSGDNE9fDjT1VVf5no3PgHFZ5/dYjAUqVSx9PQX9vba3d93YFy56mHvvRs+wDAM++GHx1TqnlWr4lpbBRiG1dS0bd2afPVq0ccf38Qw7PLlAj+/qEOHnvn6Xrh5s/irr76iUsaMGHGMzxdt2JAwb94NDMOEQvGOHSnnz+etWhWHz9bH54u2bEk+fvx5UFDs9esv3/q1SSSSmJgYNpvdQ+/1Pw0OEXvH5s3eCQlz//yzdPToc21twu7PH4waNcrFVYvL5eblvfD1NcUfQb5aqqfOMdTX16uoqDAYjNdvCl4HAtYLCgtZCKExY8xSUhayWPyzZ3M7PX9w61bJ+PEXHz+uzssz0dH5WFlZq6ND8umnNzkcjlgsJm86TCaTaWkJEzT3DAhYL0hIqHj8uBohhDce1NZW7vT8wZQp1rq6ykZGat995z15snVOTpO3t8f334/k84W3bt1SVxdLNwytqeH2yDkGkUhUXV1tbj7wh0K+H3CSoxdYWTFWr763cOGQhoZ2Z2e92bPt795ldnr+QCjEamq4tbVcV1d9gUD85ZexQqGYTld3cPCcNCknJKRo/vwWbW3ljg7h8uV3L1/+H3GOYckSF0PDX5YscYmI+Hs0tLf32YqK1mvXZnRTUVFVVaWvr6+s3P8GhvRNUCrVp5mb//r775P09FSCgxO8vIwmTBgkPRC4qYl96VJaQUHI4sXhmZkqoaGjORw+MSK4qIh9715ZeLj/8+cNEydaWVpqLl58y9JSMybmw26e8cGDB9bW1rAH6ylwiNjXXb1alJxctXmzz44dY2QGAjs7m2KYtr39/ICAwSYm6tOm2UiPCJZIMF1dZYkEw895nD2bS0wI2JWCgoLCwkL84jLoEXCI2NfNnm0v3Ryms4HAy1NTaxBqRAgRjQYoFEpNTduYMWYHDz4LCnI3N9fgckWGhq+6Wcl0wkhKqgwNTZkzx2HTpuN0ejWF4n3zZnFgoNXy5W63bpXgq+LiykaMMDY2VidWyXf16Y2Pp6+DPVifJj0BEq7TgcB0OlW+7NDERL29XVRXx42LK2exeNLDN2U6YYwebVZSwl671nPTJjuERn32mfPBgwF4qxxiVWTkhIMHn0mvgsrDNwEB67tu3y5taGg/dy4Pn2sLIURMgyTTX8DJSbe5uWPfvvTY2DKi0YC3t7G5ucawYUbr1t2bONHKyUmXeGT5BjsKCjQKBX399UotLQ0ajaKlpYTPWEmsYjCUFBVp0qug8vBNwCFi3zVpkpVA8JX0EnwqMPxv6YHAdDo1M/MT/O9VqzzwP3g88aFDAaWlHHv7Y2FhvkhqfyjTYMfN7e9G+d2c9JJZhV8VkO7q8y5vdqCCgA1M+ISA/v4WAQGW0hMC5uY2ZWTUyXTCuHuXyWLxCgqaS0s59fXtWVkNtbVcDodfUNBcXMzuahVUHr4JOE0PAIngNxgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJIKAAUAiCBgAJPp/xKGu1lz3p/EAAAAASUVORK5CYII=)
plot(post.serrata,vertex.size=V(post.serrata)$EV, layout=layout_nicely, vertex.label.cex=0.7)
![](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAASAAAAEgCAIAAACb4TnXAAAACXBIWXMAAA7DAAAOwwHHb6hkAAAgAElEQVR4nOydZ0AU19rHz2zvS++9N1FQEdSI2INiLARLrFGTiFhibLEkmqvRqKhEVBQbWFGxIagBFKT33jssZakL29vM+2Fy9+WCGpObNTF3fp9mz8zOzpzd/845z3kKhCAIwMDAUA+4v/oCMDD+yWACw8BQI5jAMDDUCCYwDAw1ggkMA0ONYALDwFAjmMAwMNQIJjAMDDWCCQwDQ41gAsPAUCOYwDAw1AgmMAwMNYIJDANDjWACw8BQI5jAMDDUCCYwDAw1ggkMA0ONYALDwFAjmMAwMNQIJjAMDDWCCQwDQ41gAsPAUCOYwDAw1AgmMAwMNYIJDANDjWACw8BQI5jAMDDUCCYwDAw1ggkMA0ONYALDwFAjmMAwMNQIJjAMDDWCCQwDQ41gAsPAUCOYwDAw1AgmMAwMNYIJDANDjWACw8BQI5jAMDDUCCYwDAw1ggkMA0ONYALDwFAjmMAwMNQIJjAMDDWCCQwDQ41gAsPAUCOYwDAw1AgmMAwMNYIJDANDjWACw8BQI5jAMDDUCCYwDAw1ggkMA0ONYALDwFAjmMAwMNQIJjAMDDWCCQwDQ41gAsPAUCOYwDAw1AgmMAwMNYIJDANDjWACw8BQI5jAMDDUCCYwDAw1ggkMA0ONYALDwFAjmMAwMNQIJjAMDDWCCQwDQ41gAsPAUCOYwDAw1AgmMAwMNYIJDANDjWACw8BQI5jAMDDUCCYwDAw1ggkMA0ONYALDwFAjhL/6AjB+Hz/8kFFW1j12rAGRiG9uHjh5MvfVqyUTJxr/1deF8XowgX1gmJgwvvvOC93++OPowEA3TF1/ZyAEQf7qa8D4HSgUMIGAAwBcuVK6f396WdlqBoP4V18UxhvBnmAfGKi62toEW7e+vH3bD1VXR4fw3LlCNzf9PXtum5mlrVlz6uzZIn9/uzt3qnbsGOvrayWVKs+dK5TJlOnpbdeu+TKZJNUJDx7M7OoSIQh49Kj2v5RrY2N/VVXfzJkW//Vd/nPAjBwfJF99FT9/vq3qpywUykeO1Js926q5WV5QUOvmptPSMrB+/ah9+7yuXCkFAJw8mTswIDM3Z9XX8yIjy1TnaWzsj4wsCwmZ8vPPUzZsGNXeLvhvrio+vun48Zz/5gz/PDCBfXhcv16em9tx4oQP+lImU5qaMoVCOZcrtLW1OnToMIVCplAIEAQ0NMg8nhQAkJXVvmyZ06JFDsXFqzZscFOdis0mNzUNPHpUCwBYv34Unf7bj6+enh6xWPzaXe7u+n/C7f2zwAT2gcHlijZvfnH27HQNDTLacu1aeVhYUXPzgIkJUyBQIAgik8nQXaoJtokJMyysEAAgkSjS09tUZ9PUpOzb57VgwaONGxMhCDIyYnR0CL//Pu3hw9rly+MAANHR1VOm3Dl7tnDy5Ki4uHq5XG5m9pmBwXdSqWLHjuQlS54AABQK+Mcfs27frrx6tRQ9rUym/O67tCtXSjdsSHj8uO599c3fEUxgHxgbNiQwGKTWVv6ZMwUhIfn+/o+7u8V2dpq3b1eeOpXn5KSdlCSNj6/p7BQ1NvYnJbU0NPR3dAh37RqXltbq5XVj//50Dw+DwSfcu9czKWnRs2cNEyfeFAjkqtFmcXEXlysSiRTJyc1VVb3oaJNIJE6daiOR6Bobn9HU/PXHc+lSCZGIW7zYYdEiB7Tl9OkCGo2werXLgQMTli590tcned/d9LcBE5h6SUtr1dYO/fjj6D17UmbPvn/xYslbDm5s7H/+vPHtJ7x3b25T0xerVrnU1vI6OoQTJhh3dYn27EktKVm1Zcvohw/nbdvmMmkSu65una/v/W3bxtbWrjUwoBsbM9LSlmZkfHbkyCSJROnkdAU9W3V1HwDgo49MMjM/6+uT3rhRrhptkkh4hQKeMsWMTiddvVpKIuHR0WZAwAYSieLurqdQtPf09EgkkqSkFktLNgCASv3VZlZQwCWR8AAAHR2quTmrvr7/v+3HDxZMYOplwgRjZ2edefNsDh366Pp13w0bEhob3/hre0cjAQwj06ffdXDQOnz4o6+/Hn3s2GRDQzoM/zoa1NTU5PF4DAYxI2Mp2tLeLhy8MXhXUlJLRkYbAEBbm2ppydbUpKhGm/39UnSESaMRrK01YmN/HemlpHA0Ncl0OtXGxgWHIzx//pxCkd6+nQvDsEIBK5UIAMDdXT8lhQMAQBAgl8MODlp/pO/+EWACUzt4PIRuIAjA4yH0r/21vKOR4N696tZWwbp1rgCAtLRWHZ1QPl+mekJqaGj09fUhCCgs7AQAVFb2bt36cvCGahcAwNKSvXFj4smTebt3p7i46Pj72w0ebZ47V5iayhEI5IsXO1y5UtrQ0F9fzzMzY9HpxKKirtrafqWSlJBAra7OiY6uHD36/L171Z2dwqlT77S3C+rr+8+eLZw+/e68eTY7diT/z87EsIVmtePjE2ViwnRz04uJqfP1tdq+faxq2So6uvraNV+FAj56NMfKip2W1lpZ2Rsf/2liYtKcj4+s/XIpDBvNnGk5d6714IUsExNGenqbsTHT2VlbIlG8eNFMoxF9fMwyMlpfvWpdutTy0aPoTZumHD7cB8PbwsKKzp0rPHhwYmurAN1gsUg+PlEwvO0dr//QocytW8eYmITduTO3qanfx8fs/v2a1FTOgwfzgoNz5XLlnDn0I0dOP3jgyOF81d4unDv3QU3NWgRBTp7Mk8uVu3aN6+4WW1hcaGn5UlOTotau/huCPcHeByNG6Pj5Wc+ebRUeXtzQ0D/EkDDESKBUKsPCSglkkq+vrspIMHghq6ioS0uLUlHRc/jwpBMnfPh8eXZ2+5QpZjt3jkMQJCDAGChuEQjNLBYZAODmpmdkxPDzs1ZteHubortQysvLo6Ki3n79VCphzZoRoaH5zc18dLqFgs61XFxcrl8/b2HBrq/vx+MhXV0qBAEcDsJmYgDz5FA3sbH1SUktCgUsl8MFBZ1MJun48ZwXL5p37/ZUGRKSkloWLrQDAFCpBARBUlNT+Xzm/v37Pv54LAAA/WmmpbVWVvbV1KxZtMjh/PmiH3/MNDJiQBCAIEgsluNwuMTEppcvm4lEnLGxRWVt27Ztz+Ry3oYNCZaWGgCA69fLz58vbmvjb9nyYswYA5FIsWvXK4lEcfSo99y5Z+rqcp4+pSxY4DJ3rnVcXP3Bg5kBAfYJCU3jxhkaGjIiI8t0dKjr14+ys7ukp0c7ejQ7MrLczIwpkyk7O8U3bpTX1fFmzLBE51ocDl917+hMbOvWMf/LMzHsCaZeZs+2IpHwHh6Ge/Z4LlrkUF7ebWjIWLXKZbAhwdiY8eJFMwBALJb29PQxGIzp051SU1vBICOBhQXb19cSACCRKBwctIlEfGenCP0IBoMEADJlitnkyaYKBUKhEG7f5jCZZDweOnBgwnffpQqFcnd3/cbGfnNz9okTPu7u+koljD79Tp8uGDt25JEjC48fn4E+KidONKmv523ePPrChRmhoQXm5iyhUBEcnGtpyfbwMDAyYlCphKam/oyMtl27Xnl7G3t6GllbawYEPD52zJtOJyYmNtfW8kpKugEAQUFuRkaMCxeKjx7NDgub/i6r2P88sCeYenn1iqNQwKmpradO5T192qBQIEFBbuHhRbdvV9LpRNSQsG3b2MWLn6xeHScQdJHJRG1tm40bGfX1/RcuFPf1SdCf5s6d4z7++K6X1w1vb9ODBycGBo7auzd17dpfnJ21FQqlkRGjpoZXVdWHx0OHD2cCAEgkBp/fmpXVZmHB4nKFkZFlRkb03l7xqVN52tpUGEaePq339bUqKODSaOQNGzYwGL+O4hwctIhEPAQBNptMIuGnTzevqvrcxuYiAEBPj7Z8ubOVFXvjRncAwLJlsSYmzDNnpgEALl0qNjFhAgACA0cFBo5C751EwqN7/5fBBKZeJk0yMTJiHD78kY4ONS6ufseOsUVFnTt3voLhbVyuiEzGa2lRnjyp9/Sszc8n+/qOysuTWViwo6OrKyp6nZ11nj5tGDFCBwBQV9dXVtajskxs3Tpm69Yx6PaDBzWffeb45ZcjAwNHublFODho02iElBROYqK7rS1ZoUAKClZwOPzk5JbCwpUAgKqqXk9PQ19fKwCAoyMrNpbDYDBeO4obYgBD3UGOHvWWSBT5+Z3YCPBdwAT2PnjwoMbRUXvvXq9Jk0wAAKiNISWFk5XVHh4+8+XL4i++OOnkNH3Dhq1jxlzjckVjxxqg3rr29lphYYW+vlZDLBMqXr3ilJV1P3xYKxYrioq6RozQ3bJlNImEq6vj5eYSbt5MRh+AqpHbiBE6g7fnzNEsLmYPflTGxdX39UmqqnobGvo7O0VFRV0dHcL+fmlVVe+uXeMCAh6rnqJjxujX1fEGv/d9d+uHAGamVzumpuevX/f19jZVtWhqnu7r2yiVKr/++mVCQmNgoCaLJeDzrRcuHDl//qOHD+cBALy8brS1CZYudWpt5ScmBtTV8eztL925M3fBAtt3+VA0cgSHqxkxYoSBgcGbDktISHB1ddXT0xvSLhDIPTyul5ev/kN3jPH/YEaOP5M7d6oIhODt25PFYgUAoK6O5+JyVSSSo/4NKLm5HRKJUi6HU1I4L182f/mlZkQETyAYIRAQBvtPsNnkpUsdhUIZDgcBAKytNfB43DuqC/zbKcTBwaGsrOxNx0ilUj6fr6OjM3zXYG8PjP8GTGB/JgEB9kuXOvL5MtQrz9paw9KSzefLbt6s6On5NcSju1sskSji4uobGnpdXfFstt4333gN95/o7BT5+9vFxtbX1fE6OoShoQUymTIurh49iVAo//rrl99+m3LyZN62bUmjR18bciWoU4iJiYlMJuvq6lK1D3Z3bGtrMzAwwOFe8xtAEJCR0aZyWfwD9PZKNm5M3LMn5dSpvA0bEnx9o4cckJvbceBAulwOD24UCOT/zYf+HUEw/lSys9uZzJCBASmCIKmpnHv3qoYfY2Bwtr9fGBf3NCIi5e1n8/C4fvp0PoIgO3YkS6UKiURx8mTukSOZWlqnT53KvXevyscnKjS0QFs79MmTWgRB5HLloUOZt25VBAUlTJt2B0GQH398ZmLy3RdfPNy8ObG9XTBnTvTIkRHLlsUiCHL8+BN398snT+bOnh39ww/p4eHF8+Y9CAsrRBAkKakZgo7xeBL0I86cKfD2vh0bW/eOnSCTKd3cIu7cqURfSiSKuXPvDzlGIJAxmSESiWJIO48necdP+SDAnmB/MmPHGjg6akdElAEAHj+u8/OzHhJhBQAQiUTHjx/v7KSvWpUJAJBKladO5R09mj1v3kM+X3b9ermf34Ndu15t2fIiKMgtNLRAJJKjToyoP0d7u5DPl+HxEGoLCQwcdeKEDxq5PDxyREtLzOHIxOL733xjJxDIJk0y1dGhFhd3tbcLDA2lHI5IteS1erVLaOi0CxeKAQCoTaWwsFNlblEFR78L166VSaXKTz+1R1+Syfhvv/VU7U1PT+fxeHQ6cbhdZLCf5D8DTGB/PkNUMcQxSiwW8/kDUVFRK1d+hBoGh8Tzu7vrq9ygAgLse3slgYEJ8+fbgn8HJuNw0MyZlkFB7ng8jkIhpKZyhEJ5YyP/+vXyH37ISE9v3bLlhUgkr6/nPXxYGxsLW1npLVzo/8MPqWFhTx8+LJdK5SQSnsPpfPlSIpXC8+c/xOMhpRKZPv3u3btVpaXd6EBUqYR9fKLQjzhxIvfFi6YXL5qPHctpbOz/zWDK3Fyui8t/TO08PQ3R/5Ft255MmBD+5ZebVLu4XNG5c4VRUZUXL5a8etXi4/MbflsfFpjA/nyGqGJwhJVcrszPz2ex2NHR0QAACAJgWDz/YHc+Mhm/bt2I7OyOsWMNwL9XosRiBQQBVWCyjY1mUFA8giDu7voikdzQkHHihI9AIKPRSLNnW1VX8zQ0tJ8967WwMNXRMdPVxfP5vM7Oga1bb0skCB6Pq6/n3bxZQSTiWloG1q51ZTCI6JMKj8eh+pfJlDExdfPn27LZZAgC0dE1vxlMKZUqpVLlkEb0f8Td3VJb21YmGyEQ/Jr/A12uWLDAbuJE4zetRny4YAL7I/T0iDdvfhEaWnD2bOGmTS86OoSD9w5RxeAIq9LSUrlcTqVSDQ0tnJyuoEskQ+L5RSJFUdH/myW++mrUmjUj0G00MPn588aUFM7YsQaoLUQqVZDJhO5ucV+fxNpao7y854svfklNbUMQOC+vQyxWNDYOlJQgJJJg1Cjj0lKJUEjS0VGmpvIaGnKlUsXdu3MdHbW7u8UAQGlprUKhvK1NkJTU0t8vVSjg1FROX5+kvV1w61Zlf7/U1JT5Li68I0fqlpV1q16iJpkrV0plMjg/v8vc3PTECf/791+gIvTzs6bRiC4uV7hcoZPTFQj6j1N98GaPv3oS+OEhlytdXK4kJDShL1NTOZaWF4RC2eBjmpsHjh/PQbefPq13cbly8mTu1KnXFi++mpraTCAEx8TUPnlSB8Cx2Ng6Doc/fvwNT8/rO3cmy+XKM2cKtLVPFxd3vekCJBKFo+Pl8+eLVC0sVkhlZU9lZY+n53W0JSQk7+DBDARBbG0vtrQMbN4cP2vWBalUmpnZ5uMT1dbW5uy859NPryIIIhbL09JaORy+i8sVBEFyctpR6wiCIGz2z+hGRERpYmITan4IDs6ZN+8BgiAwjNjaXhQI/uPGUQQCmbn5+YiIUgRBlErYy+vG5s0vNmxI2L49CYaRjz++l5LCOXAglkg8LhbLU1M5MIykpnK8vG7weBLVh6r4oM0e2ELz7wad6lRXr1G1jBt3Y8kShy1bRr/lXe3t7bm5uVOnTqXRaAAABAGvXrUMXn1W8ZZdKvh82d69qTQaUU+P1t8vFYnkBw5MuHKl9Pvv0168WDRihM6zZw3btyevWTMiKanF2Vl740b32bNvWVoyLS0NCgq4mzcbEQi47duLmEwW6pYRHV29adOLrKzP7t2rDgsrSk1dUlnZ6+MTFRu7wMPDcMyYa25ueiwWed48m9mzrb7++uXIkXp9fZKxYw2mTDF77RW2tQm2b0/W1qYODMgeP67Nz19OJOIDAh7DMOLqqhcS4nP3bvXKlXEhIaOYTIPGxgEvL6POTlFvr3jLlpexsQtQT6537I2/M5jAfjdbt76sru578mSBqmXFijgAgLu7fmxsfWzsgr17U1ta+LduzUGjJEUiaVRUVmCgJoI4xca2olGS8+fbolGP0dHVZ84UDk4Smpzc8rsCIt8RPp//8uXL2bNnQxAUExMzderU+Ph4X19fMvk35jyVlb1VVb2mpsy+PumLF03bto39XXGTW7e+rKnhxcTMR1+mpHBkMphIxLHZ5DVrntnZiVNTu3/4YcaqVWNUN87liu7fr9bSovD5cltbDXX0xnsD80X83VAoBKFQPrhFLFbo69O9vIxiY+tJJPzUqeZoArOTJ3NlMlgiaSku5ty+3X3mjN/x48WPH89HEASHg9DZ/Lu7Hb6d3l7J99+nsVgkXV1aTU1fQ0N/XNzCwQcwmUxNTc3m5mYCgaChocFgMPT19dvb2y0sLN5UUCI3tyM2tp7D4bPZZCIRt2yZk7W1xmuv7S2uVWKxgkD4/3mVjY2miUnYs2f+6Mvx4+W3bl06cYK2atUY1Y2rvDTr6ngODloftNkDE9jvZtIkk8uXS1Q54tGlm0OHPhp+ZFZWe3Cwj5GRm6Fh54IFCwSCX82DEASBf5sQUTv44CShql3vjlwOT5t259tvx6FLT1KpMiDg8fDD7O3t8/LyGAyGhYUFAMDQ0BAV2JsKSjg6ak+Zcqe4eOW1a+X/+lemuTn7889dVClGBvMW16pRo/SOHctBkF9vytCQzmKRzMyYEomSySStXr1KJtOIioIH37ifn/WLF80uLleuXfP9A73xtwKzIv5uZs2y9PY2PXUqD30ZGVnm6KgdEGBPoRAEAhkAgMeToM6HqHmQQqGsW7e+qWnoz2TI2HzwWP33DtvfvrCrQk9Pj0AgdHZ2GhsbAwCMjIy4XC4MwytWOKMHXLlSWl7ec+TIJPQluhZsaMjYt89LW5uyerXLmzL2vGWBeNUqFxIJHx5erGpBE2AhCIIggE6nz5gxE4Jwg288N7fjzJlply/P2rLl5R/ojb8V2BPsj3Dz5pwrV0pnzLgLw8gnn9hER38CAHB11dXUpAQGJhga0vv7pZ2doiHxHRcuFKviRFA7eFxcPZ8vG5IktLKyF92lmuj/Jq9d2B2SWufy5ZLbtytNTZGXL7kVFbl793reuFF99iw3JSUhMbHd39/+889dtm59uWiRw9mzBenpbZcuzQwJyReLFZs3v6DTWzo7Ow8ffpqTI+PxJHQ6UalEbGw0P/rIOCgo8c4dPzIZ/6aZEpmMz8r6bO/e1G+/TUFNMl9+OdLMjHX+fFFdHa+paSAxsbm5eYDD4dfW8tAb7+wUxcc3eXkZBQaOUnXUu/fG34u/0IL5oZOd3e7gcCkg4PFwh7r3zMyZ9/D44E2bEgc3jhwZoaV1Oj291cXliq3txaKizhEjrjx//suNG/dJpGAOh//sWb2+/qlvvtnW2cknk09Mn37HzS2CzT5NJp/Q1g6dP//h4cOZBgZnW1v5EHQMgGMHDz7V1g798cdMpRKuquql00/t3p3yySf39+9Pk8mULNbPjo6X/6oe+NuCDRH/OGPHGhQWrly+3Onq1VI0qcZfha+vJYWCv3q1VCD41fpSVNQFw7C5OcvUlEmhEB4/no/HQ3Q6XiqVjB7tZGND7+gQ5uVx+/okwcHHa2qKTUyYBQWdRkaM5ORF06dbyOVKEgmPDgg5HAFqpZg1yw2CAJmMx+EgOztNAwN6RETpuXPTg4NzYRjB4QAW4TIcTGD/FWQyfs4c6y+/HPmm5aD3w5o1I4hEvKYmRVWa6NixHDMzlkSi+LcHSRcAQCKRWlhYWFlZCwRSU1OSrS0bjyfs2bMHAGZjY39o6DRLS/aRI1njxhmy2WQuV4hm542IyKVS8ailgU4noo0KBdLaKggNnWpoyKBQCOCf6Kf7p4DNwT5UoqOr/f0fr1vn6uysI5EoyGScjg5l375UuRyGYdjIiCEUyuvr+4ODc7W1KZ9+GlNSsrKxkX/oUD2V2qqtTf7mm4PW1uPlcqiiwj04OJ7BIHV3i8Rixe3blVQqwdSUWV7eo69P43JFz541fvSRcWxs04ULxUKhvKaGt2rV05iYchqNoFQily+XSCSKV684/f3SN63s/dVd9ZfyV49R/4doaOA9e9YwvP3ZswY6/dSQGdS4cdd1dEJZrJBZs+7t3v3K1zc6PLx4yBt1dc8kJTWj2xs2JLS28rW0TickNF26VFxfz/vxx0wC4Tg6P2Sxfk5KqjAwCLl1qwJBkI0bfwLAnMn8icE4BcPImTNPJ0+ORM/DZv988GDGDz+kU6kn9+5NNTQMXbny2vHjORTKCYlEcfx4zuHDme3t7UQCHYc73NsrRhDEwOAs+imol1NTU7+NTTgMIwkJTf7+j/7ULvzwwIaI74831XaYOdPCx8dsyAxKLFZ4eRmNHKn3lsIR6JIUWtIBQRAjI8bg/LskEl61UgfDcHV1i1AI6+nRAAB79qzcsmX+4sW6CIJAEBg92qG9vV+hUIB/Lzrl53OXLXMqLe2aOpV15Mj/p14rKOAKhQPZ2dnXrl+ytNQY4un7lpW9/1kwgb0/3lLbYdw4Q2trDdUMKjq6OiDAHo+HfrNwxPnzxb6+0cuXx3E4gtu3Yzs7e2Ni6l69agEAwDCCriC9etUiEMj37g1xdtYKCko8ejTb3/8pk+lWUoJXKOBFix6Ehpbh8YRNm2KvXCkViRSFhZ0mJkwIAk+e1NvYUOvrYfRUMAzr6ysSE+u9vb0DAhbhcMRhad7Af778kBew/iSwOdjv4DfdkV7L4NoOAwPSAwfSP/98xNmzhdevl0+ebKpyVggKcjt+PDcwcJRQKKdQCGQyvqenp6WlOzmZER1d3dLCHzVKLzeX6+EBHT+eIxDIy8t7xGIFjyepr++HIFBe3vP4cQmTaWxoaJKbyw0Ozo2IKJXL4cOHs9DHZm83zt4e9PSIdu1KoVAAg0HPzBQQCLhp02Aul/nokbSzs9vU1FUmU0ZH14SHz7h8uYTFwnG5VD096r17VVKpctu2e59+aiISUaOimlSp2rKy2ru7xfHxjQwG6U0rewYGdLV/N39b/tIB6ofEu+SZeC1hYYVHj2YjCJKSwvHxiWIyQ8rLu6Ojq7u6RK6uVzs6hAiCHDyYIRLJB8+gjh/PMTL6FoCvxo07v3Vr4rFj2VZW4TTayby8jm3bkpYvjx0x4qqWVuh336XSaKdgGDl4MMPB4djIkcd+/PElAMcIhODq6l4Dg7Pff5+2Zs0zPP64puZJD4/LW7a8mDEjvLW19e7du0xmiKfn9Y6OjkePHvn7R2/Zcr+kpARBECeny7m5HW1tbU+fPkUfXm1tbQ8fPqyoqEAQhM+XYetd7w42RHxX3tEdaTiDC0Di8RCdTjQxYQoE8qSkZrT4A3rY8Aomrq6urq6ufD7c38/x99c1M2NSKASZDHZy0gYA1NT0KZUwgiAkEg6CQEVFj1DI0NMzdXYewOEgGEY6OkQKBRwZWbZ6tQuCgJEjtWg0xaVLRTt3Omlp6VVX/1rHWV9fXyQyffiw4cEDbm1trUQiRdP0Tp36cNWqupycjrKysry8vAkTJjg4OAAso9vvBBPYu/Imd6QhJcNdXc+OG/e1t/dtNLNFR4ewqan/0qWS5cvjVAUgw8OLX71qCQiIUWVBRN3z168f9eRJPVpDGYYRIpHMZms6ORm1t1Nqa+uam7sQBKSmcjIy2iIifBkM4sCALDOzXVVUEo/HZWVxi4rYOByAIBAdXdXTI1m50vnChWIYRgQC0NQk1dCAdu+u//rrl3Z2FIUCRl23Fi928/Y2am0VJifDu3Y9Cwubnp1dL5Eoxzk6/joAACAASURBVIwxkEhqurq6pk2bpkqfiCDg+vXy4XVxVXU03513KZn7oYPfv3//X30NHwb379fw+bLFix0GN3K5IgCgOXOs//WvjAUL7PT1aYcOJTU2Hjl0aPvt2/UBAfZcrkhbm5aby01La1UqQW+vuKNDtHq1S2Rk2cCAzMvLqKlpAADo4MHMESN0Ro/W7+kRBwW5NTYOhITk5eR0tLUJdHVpfL7i5UshDIPVqzWpVEVEREN8fBODQXJz05s50xKd9uzc6VFXxyMS8f39SgIBeHqSY2M7bGw05s61KSzk0mgSsRivUMiWLTMUCCgyGfzsWXtjo/Tu3bleXkY//5zf1ydrbh7Q04NnziRXVIATJ7KlUggAkZYWKzy8RypFtm9P1tOj2dpqvnrVsmhRzJgxBgEB9ng8ztqavWnTi9ZWwe3blRUVvY2NAxcvFiuViL39b6eqv3u3Ojy8WOVq/I8EM3K8KyNH6p4+nT+k0dSUmZ7epqr0hcfj9PW1f/75ib6+OY/XhB5AJOJu3PCdP//R7t3jjI0ZhobnZsyw8POz1tQ8jWbJBgA0NX2BboSETAEAWFqyf/nl0zt3qpYufeLqqvvDDxOoVMKZMwVBQYmenjpHjxqMGGHt6OhIJBIBAJs3u6PvPXduOgBAIJBbWYXfvDnd0/PJ4cMT2tvF+/a5hYWJamsVBw+OCQ+vHjPGasUK13PnYktLZV5eRiEheUZGDBcXHVtbjYiIUiKxTybj0WhwQ4OEQmH4+487ezaqo0M4dqzBpk0vcnK4x45l4/G4kpKuJ0/IqamtYWHTEAS5cqVEQ4OiqUkJDBxVXt6zaFHM3LnWUVF+QxadJ00ydXePXLLEAXVB3rTJ7dSpXD+/B2gQ6qlTU97LN/lewYaI78qaNSPkclhlSQcAxMc3nTs3tGQ4giCurq7Ivy3Uw2uKIwjYvj35229TJBKlKinva1O7DMkTvGGDG5VKOHRoyoIFs3g8sbV1WH19PfI6UziCIE5Otv7+ZgcPvmhs7NfQgIlEokgkcnd353LxfD4vPr7p+XOemRmzvr4/JYWDLo59+62nQKAoKJDeu3edTn+1Y8c4OzsdPB4nFCpoNMLixQ4mJsxjx7InTjRWKmGpVEmnE3E4sG7dLxAkh5XteDzU0SF4+LA2O7vD0pKdm8v96afsjg5hWVn3/Pm2ZmbMQ4cy16+Pv3PHz8VFt7Cw4+XLpk8+OVJTU1pe3n3gwAQmkzw4FVx0dPWUKXfOni2cPDkKHWwPbxmecPJvCCawoSQntzCZIUuXPgkOzt2+PXnfvlSJRAEAoNOJ6elLnz9v3LTpRUhI/pkzBfb2mvb2muHhxTjccZlMuWxZ7MKFjzo7pQMDeJWFekhO7IyMts5OEQzDM2daSCQKHx8zQ0M6DCNvshxs3Oh+82YFny8DAKSltZLJeB5PSqVSJ00al5IS0NTUlJCQ0NXVNVhmdDpxyhRTAMC+fdOysvh0ek9vb69QKGIwWGQyeeJEy5QUzqhRuhAEyWRKBwctKysNNE+bUon4+Bh2d0tpNI8FCyZDEL63VwIAkEgUJBIeQRAiEWduzrKy0kCNKG5uurNnGyoUChqNYmtnyWAQeTzZ8uVODQ08dEqZl8edM8eazSaLRHIXFx0CAeTkcOLisi9cSC4rq2IwwOjRdAaDoadHCw0toNP/IxXc8ISnw1uGJJx8D7+NPwA2RByKt7epsTHzs8+cZs+2AgDs3p3i6XkjJ2c5kYgzMmLcuDF78MFmZqzGxi+MjM5FRHzs7W1661bZzp2JTk6GTk6G27aNBQDMmmU5a5YlAABNiXPnTpWpKfPnn6ficFBb23pDQ7q9vRYMIxAEFRZ2Dk/tosoTHBTk9vhxHR6PO3Ei69mzTCsruydP6p8+XRgZmb1kybX+/kIKxWXJEuugINtvvy2Ji2sfGLihpUUcOVIjKqpl//58PJ5JJFI++ugWBEGNjfLNm3/p7JRfuDCFz5fhcFBxcdfdu1UVFc06OtKuLqWWltHWrUIrqwpzc2ZqKkehgJ8/b4BhpKGhH0EQtJYFDCvr6mpiY3twODwEQd3diLMzE0GAu3tkT49k9Gj95mZeRkaJvX0NjYafP/+WpibU14dIJMjevZVOTlojR1rHxhZOnTrl1atyHk9aUMBFV+FVqeAMDOhD3EGGO4gMGZ+/jx/H7wd7gr2GwTHqBw9O5PGkaCrsN4HH/9qNOJysu1sxuJbKEDIz21xddXE4qLKyd+vWlykpnKamgfT0NlVG21On8qZPvyuTKXfsSF6y5AkAYPly54MHMyIjy0pLu+VyOC+vLTz8WXx8Ulpaa2DgtYMHc3t7YQKBqlBAFy825OWVW1nJAQBffEGXSvmtrYKeHpFEQoZhmaMjPjOznURSmJoyXFwI+vqkjz4yFArlXl5G+/a5K5WwoyNt0SIXOh3f2Lh29WrdiorugoJOW1vNnp4gGxtNFou8bp1rePhMExOGri4ZQQCHQ3B01GUySQQCDgCopqYPgpDZs03lcmVZGVco7ORwcvF4JZ+vdHQ0qKqCRSIiHk+kUglcrjg5uQWHAyKRnMUiNzT0GxoyUCf94bX8ho+B3zL8/huCCew3wOGgsWMNcnM73j7iVyqR4uKuTZuiXVx4SiU8JLm0av4QFVXV0yMCACQltZSX9zQ09M+YcffWrcqsrHYCAcfny7y8jAAAaOYc9Mza2pT+fllCQtOqVS5EIo5EIs+eTYiO/pLFIpeWEgwMNHE4UFq68+jR6QAAW1vP+fNHAwCoVEptrYRCgb75xszGRjM7+7MFCwxgGJ4wgUmnQ/fvd3K50nXrHvn53S0rqxkYqHZy0vzsswmNjQ3GxiwPjyg7O/q4cfoczlfu7vpoJdgvvnDdudPDwYHQ3Fy/d689hYJPSuro6ZF4exvr6ZG1tQkbNxo2NHQVFjaPHq2Bx+MhSGv+/Lm6uuwrV2ZNnGguFMqDgtxGjNCFYaSkZBWBgCOR8Do6NCMj+pgxBocOTRxSzRnNqTrYHWR4y5Dh93v7SfwuMIH9Nmjo4dtH/A8e1CQlNUkksSTSi9OnC4Ykl1bNHwIC7EtKuhEEuLnpGRkxVqxwJpHw6NwGhpHBFhQVCxbYOjlpR0VVmZszlUqESiVRqT4QBMnlypYWgaenLoIAS0sLDQ0SBEFubsZovg0CgWBmRhOJ4Jcv5RCERESkcDhyCIJOn25gsXAKBU5Li7Br1wihUFxR0aGhAfX0CFNSUhgMBg4H0tOXaGpq8XjiwZchFApfvXpVWlqqpWVoY2P2yy8zZDKlmRnS3881MMAbGlI+/9xr69bxU6e6ZGR8rqFBTU5ecvq0P5NJWrbMycvLyNPTaMuW0Q8fziOTCQYG9Hv35hobM44d8z5zZhqVShisYTSybtEiBy430MKCvW3b2NratQYG9OEts2ZZlpSsQk/72qRDfwcwgf02hYWdHh4Gg1PMDx/x+/vbbdo0prk5KTExcXhyadX8ISDAAUHA4AwwCgU8d671okUODAZpwwa34ZlzcnM7HjyY98UXritXPkUb09MFjo7XJRKFk5M2iSRAEOinn3LOny8iEnFc7oC//z0AwLffNg8MyPF44vPn3bW14sxM3P37vTCM8PlwcXG/QKCoqmrevz/c1lYrLU3x449CFxfd06c5+/bV19f3T5166fbturKy/pycymvXCt3czu/c+WjcuKsZGX3l5fKIiPI5c56QSIJly6xu31549mwAicSsrRW3tgIHB23VI+XHHzOHP3NUCTaG73p/X+f7BRPYaxg8iTp3rpBGIy5Z4viWET8MI+hb2Gw2mUxGq4OD180oiETIw8OgrKw7LKyotpa3f3+6q6tOXFw9+ssLDs5VZc6pru5DM+fU1PRdvlwyZ471nj2eYvEWY2P6smWWP/zgjCBALJZSqUpdXeqSJQ7z5tkYGFD37HlsY0MHAEydSsvIEGhq0rZv96DRCJMmma5c6YLDQZMmmcyZ40Qk4sjka93diL0989KlmZ2dot5egbm5TkeH/NmzTwsKxOfPzySTcatXJ12/nlhY2PPgwS+OjqwTJ5plMiaXG8hkkuLiFPPmjdi0Kd3Skv3llyPt7LQyM9uPHMn+6adJW7aM/vrr0TduVAx/5kyebIog23x9rYbvek9f7XsHsyIOJSGhqbGx/9atiro6Xk+PWKGAExMD8HjIzk5z+/ZkOp2IjvhVY5KnTxu6ukQ3b1aMHKmrrU0FAAQFuQ2pDv7kSZ3qD7u5mX/z5hwNDbKn5w0Wi/Tw4fzFi2OeP28kkXC5ue0nTuQ6O+tQKPjt28fu2+cFAFi1ymXw5RUWrmxra6urqwsJmXL0aNr69XOPHDEXCARjx4rCw63DwwUTJhgHBRkYGBjExIgjI30NDOj379ccOvRRXh736tUyFou0fLlzXV1/UFB4d7dy9eoaGk1IIiFNTaKrV+dxOM8nTzZlscjGxsYEAtTZqXj0aNXo0ec7Oqxu3JjO5SZFRTXt2DEJAICOlr//Pu13lWz/HwQT2FCmTTOXybYObx9icFfx8ceWg49PTm6ZM+e+n581ny/r7ha/fNk8frzRokW/VsTbtm0sar4HABQWrgQAPH/eWFTUtWbNCKUSWbbMGV0bMDYO++mn7OTkxajNYzC5uR0PHlSPHcv387M6dAh+8aLfw6OstrbWycmpoYHo4dH1+HGxi4ujk5OzXJ7p4KClCnlEEITJJJaUdOPxOKlUicfj09J6yGScv79nRMRjNpvA5w+gR6JGVAiCcDhAoZCNjQ07Owd6eroJBFxbmwDdNXi0/KfkTv2ngg0R/2RUy2jffDPm2DFvpRLx9LwxpBLxYFThzKpRZ1FRl44OlUDAlZX1DD/e0VH79OlimUze0FB169aEyMis778vqKkx+OUXgYODlpeXkkzGffdd/WBzHIfDp9NPff31S6kUxuGgrKy23l7xypVVd+70yWTw1aulpqY0Ph9cuZJaWdnt4XFNJFIsXPhIIFDweNIlS540NPRDEO7evYLGxn6lErGwuNDfL01IaOzoEDY09MvlyuLirs5OUV+f9NNP/z+d8N/VbP6+wQT25/N7l9HQcOb+/l//+9FwZk1NSm5uh0ymHGzxVyqVMTHVUqnyyJGklSszU1Lqbt+eGxo6n89XGBszAwNjJBLh/PmjpFKlk5P2xInG332XJhIpli51dHbWLS7uKipamZe3vLq6b/RoA4UC0dGhfvKJxtSpyNKlI5RK6KuvZqamTjM2lstkyhkzzPF4CIKg0FDvdetceTz5+PG6+/Y5CwSye/fmPn26MD+/88WL5v5+aVhYUVeXqKtL3NUl/PxzlyHGDLX28wcBJjD18i7LaGg4c0+PaNWqVbduRff2dnZ1tSsUsq6uDh+fy5cu5WlpNTY2Nsyf/+DGjfsCQYVCAVeWhUFQRlSUMCOD/9VX8SNH6tnYyMvL++ztx4webWBkxPDzsx68WlBU1GlhwYqMLGMySSEhU4yNGb6+Gk5O2itX6hgYaH/zjVdX14ZRo/SdnZ2vXvV/9WocDtcIQVBCwjgNDfnJkz4y2dYFCzy1tbt5vA1CoRy1j+/Y4aGjQz1wYMKKFc76+jQzM/ZHH5kMMWb8JX3+twITmNoZvox240YFDnd83brnhw5l+vs/5nD4S5Y4wjDC57NSU/kLF1rR6XSRCJ4yxYZCYVCp1GnTfAIDfRgMsrPzpI8/noXDQQkvU/395xGJ+NmzrerreaWlLbm5lUwm87vvMg4ezMjN7Vi/Pj44OOfUqXwAgI4OVUODPG+ebWhoAQAATUlAIAClUqlUKnk85WDlx8Rwjh7tTk5GlEq4sLB44cJTX30VZ2R0bvPmlIiI/s8+uzN5ctSsWbeDg19RqSflckV5ea2j45Hu7p64uHIfn6jvvksLCcn39r69Y0fyxYslw+Plhjjs/uPBBKZ2hi+jTZ5samBAX7bMac8ez8WLHe7fr6FSCWw2xc1tia6unY/PyMpKCYGA/+orzxEj9AcG5HQ6HY8nIAiisviPHj3a0lJXKpVzucLeXnFRUePcuVOqqno9PAx27hzX1yf96SdvLS3KmDH6AAAEAUolsnGjO5crTExsjoqqXLLEAUEQPp/v4ODQ0sJ1cdFWLaCjJdj37vXG43FXrtTk5r7U0KiUy2XjxwMNDRweLyGTcbW1ve3t3U5ODLlc2dXVY2TUC8NKd3dNMplAoxGMjRlGRoyzZwtdXXWGrM4Pd9j9x4MJ7M/nXZbRVO6LbDaZx5P88ksjny/Lyvp1JMnjSby8jPB4aNkyZxwOXLhQfO1aqbu7Pmq0gGGkqak/JkbQ3i4QCAbEYmlxsXLDhhQSCd/Swh8zRh+CQFhY4Zkz04lEHOp/9Mkn1jo61MEpCZRKpVwud3BwcHAw4HA6B5kEf62xBADQ1fVYtcr/889X2djofPmlX28v/fPPRxGJUG1t4PHjC77/fopMBlpamI8eHdLT0zM0NIRhmETC+/lZa2pS5HK4srJ3yOr8/2BGN8xM/64IhfK9e1MpFIKeHq21lf/yZUte3vIhx+Tmdpw+nd/Y2H/tWtmaNc+//NL1Tcto4N/ui99++4pCIfJ40t27x+XlcQUC+c6dLsnJLRkZ3CNHsrS1aW5uepcvl5ibs9rbBcuXx+noUBEElJf3sFh0hQI+fjzB1VVXocDxeFItLQq6OqejQ928eXRmZtujR7UjRuhaWrJPny7YsCGhoKAzK6u9oaE/Kiq/qorDZOqNGnXN3Jyam8u1tzfp75dKpYqwsKKysu5Zs+4plUh5eX9iIoVIzG1o6M/KajcxYcbH8xEEtLZym5qUc+ZYKRQwjydlMIjo7eNwuOfPGwCAQkOnxcTUnT1byONJ+XwZ+rdSVNSZmNiMWhff5Jv7lkJ+HyiYwN4JGEamT7+7cqXzl1+OBAAgCPDzuw/DCA73H8s9jo7aDx7U8vmbyWR8f7+Uzf7/xdYhy2iRkeUPHtQ4OmqfOOEzaZJJXR2PQMDt3es1Zsw1Go0YFOT+5EndN9+MravjicVyuRz+5psxOTkdO3YkBwW5lZZ2z5tnGxV1g0LpXbFijLGxa1FRV0xMXXp6G5cr0tenoZ/o7W2KjtnweFx3t/jy5dLAwJFSqbKpacDDA2psFBEI+NLS7r17Z2dnP9+165WTk/aaNc/z87lyOcLnyxEEWbjQPiQkNyamfmBARiLhd+0aN3PmPaFQuXRpTGLi5zk5HTCM0OkENHPbiRO5EomivV14/HjO06f1QUFuxsYMHR2q6m8lKanl3LlCBoP0loxu/8CMOu8zhdWHS1RUpZnZeTSLE0pNTZ9crhx+JJpHGoYRVVLr12JiEjb4AKlUERlZ1tIyMGbMNQ6HL5Eo1q+Pt7O7mJXVtmDBw/j4RgRBCgs7J0++nZnZNm3aHQRBXF0/JZM/v3Pnjlgsi4wsO3AgjUw+yeHwkUG5rEmk4ODgnIEBqYvLFWvr8MDAhEuXinV0Tq9ceamwsJDD4ZPJJ3JzO65de+nldQVBEDe3CBubi6qrqqzsMTP7WSgUDr5yGIbj4+Obm19/d7m5HejlvRbVhb2J3+y3Dw5sDvZOqOK40JeqOK7hxneRSJSfn6+K73pToLtIJD94MFN1/iEztMElHl1ddZOTWwAAIpEcAIB6A4vF4unT/T7+eLauru7RoynNzQPffDNWqYTv3q1SnTM+vgmHw6WkcJhM0rZtYyUS5d694yZMMCYSQUWFxN7eHkEAggAHBy0DAwOJRAIAkMthNJewRKJIT28DANDpeD5fNHgtDoKglhaNOXOehIbmq25KoYB//DHr9u1KtDg1+M/AtunTrzo5OWVkZKC74uLqx4+/eepU3pw59//1r4yLF0vmz394/nwRAEDVb1Kp8tSpvKNHs+fNe8jny65fL/fze7Br16stW14MWRtU35f+p4AJ7J0YXsl7xoy7UqlyiJWstbV1YGBg//79Kme81wa6K5WAz5fl53MrK3vREw4Jbaqp6TtwIF0gkAcGjlq50vnBg9pPP32MZjhDvYFXrnwkFjNEIlhDw4RC4d++XRkeXjx1qtnFi8UBATGdnaItW1709UkkEoVcDp84kff11y8sLFi7dqUcPZopFMoZDEpkZPXOnck0GqGrS1RYyO/okFVWtrq46HR1iby8buzfn+7hYZCY2NzWJjtwIGdI9M3UqXa9vcqpU2mqm7p0qYRIxC1e/KtHGAAAdfIiEnHOzuT29vbamorOzl+LG02caFJfz9u8efSFCzNCQwtWr3YJDZ124UIxAEDVbydP5g4MyMzNWfX1vMjIMtS2efjwpBMnfIZHA72nH8EfApuDvRNvquQ9JGrd2NiYxWKtXbsWvLkSgqkp095es75+3fz5j5hMEnr+Nzk6opSWrrp3r/rTTx8bGjL6+iQPH/rFxcXNnDlTqSR4eFz/6SedtLQFLBZr+HsDAuwBAJWVvdbWbFNTZl+fNDIy7eJFW0/Pkaampl984Yoetm3bWF9fuljc6eVllJvbkZHxGdoeGDhq1CjR8eOtNja64D/j+VksWnV1lZ6eO3pTSUktCxfaAQDQ/DwoMpksPj6+o0NkZ2dXUCAjEolffXUOAIDHQ0QiHoIAm00mkfB4PKShQUa9HFX9htaPt7Jio4qtqupV1Y8fnl9g9GjKf/f1qhHsCfZOvKmS9/AYFiqVCkEQgiBvqoTwxwLd/f3t2Gwyi0XS0qJUVVWZmppmZ3cbGp51dtZKTgYbNz5XJecB/5m3Z/nyuClToqqq+qytNeh0OYUir68XUChDkxZaWFg0NzevXu2MZs5SpQQtKBA6OLAGR9+IRPITJ3IhCGdmZlZX9+sIzdiYgdb4VCjgvDzuF1883bfv+atX7b/8AsfGiktKeGiGueEM7wG0Aa0fDwYNVlW8JRrobwj2BHsnhlTyjomps7HRDAnJj46u1tGhqozvc+fa9PRIyspkZHLFmyohvCns5e2IRHIIgj77zKGxsa+urm7atGkMBsPYmDlmjMHGjSPCw8Orq6mq5DyD8/YIhXItrdCTJ3Pr63kdHW1jxmg7OmoEB+cdOTJJdfLnzxsXLnw0c6bmmDHd6elLt29PjogoS09vI5Px//qX05IlxmFhbaroG3d3/Zkz7zEYJAbDNDGxqK5O1tEh3LZt7OLFT9aufc5g4NlsnIeHaOHCkQkJif39uMrKPl1damenqKGhHy0TgcNBfX2SgoLOTz550NkpKirq6ugQ9vdLq6p629uFaL+9pX788GggtX3tfwZ/nX3lQ6WhgWdr+6up7aefsqqrexEEaWsTqA4oKipCqyj8uRgYnE1Kai4vL8/MzERbzMzCIOiYVCqlkfG6usu1tUMnT7796FEtgiDGxudcXa+eOVPg7X1bU/NnM7Pza9ZcsLE5/fDh087OzoEB6WCLKIIgc+bcZzBOhodHymQyBEFu3qxgMEL8/O5XVFQUFRUNvxLUGFhVVZWSkoI2SqXSoqKiBw8eFBQUymSyysoeAiEYQRAeTzL4s8RicWJiolQqRXf96b30dwN7gv1u2GxyU9PAo0e1n3xis379KD5fVlnZe+BA+q1bc9ADDA0NCwsLXVxc3n6ePwAMw7W1tZMm/frwodGINBqRRCJ9NPlfCgWsUFCMjelLlz7ZvXucRKKoru5dtszJ3l5rzpxoe3vipUvZJJLjwYMNY8YQwsOLW1u/On48+9Kl0hkzLJ4/b3R319PRgdatC//++8Zx4yyTknrIZJxIJOjq6g0Pbxg1qr+kpG/6dHM/PysikQjDyNdfv2SxSAUF3I6OLi43f88et6NHC+7dGzdt2jQcDicUCnt7ewEA9fUN6entLS38p0/bLSwoYrHcwCAvNLRy7dpaEsnK2Ji5e/e4f3bJWSw3/e+GSiVAELR+fXx3t3jaNAttbeq9e9VPnzaYmDAtLNihoQW5uT179yaYmPArK6GgoEShUK5K7C6TKQ8cyPhdCdxVBAfnennRjYyotra2AAC5HD53rnBgQLZ7t+ezZwOTJo1sbxeIxb1crtLYmFVd3ScWKwwN6TY2mhcvlri7G9FoioEByldfOfzww9QjR7JWrXKxsmKfP19UULBSoZBGRVXNm0cvKjKxt9c5c8anro5fUNBrYkLp6pJIpZLubgmLpdy9u1Ch6JBImq9e7eju7m1q6uLxhC4uWVlZ4Jdf+kaNogoEfImkpaio/vr18qKi9pycAQcHwtq1uZs3O1y6VA9BRG9vs+vXYWNjzatXl1Kp9HXrftm82V1Hh3r6dH5kpK+ZGevixWLUMPOPATNy/BH27vVMSlr07FnDxIk3BQI5miLKz88aNS7r6hI4LcI9e+4Pt9H/NyZmpRJpaWlBawiJxYrLl0tUu9B5f12dePZsV4FAPm0alckkGRoyNmxwQxAEAFBZObBlywIcDqLRNLhcIQQBhQImEvE4HJSfn9vX16RUgu+//wRBCJqa2mlpgtWr3YlEHA5H7ejA1dUp7OwsQkICbGw0T5/uaG01gSAIgqguLmYff2wqFLoaGwsNDAiWlqSQkC5raxcSyZ7P1wwMnAZB0FdfzWCxyHZ2Nqam7K4u2bffTra31yESLQ0NDefOtcHhoMEZgf6RDoqYwH431dV9AICPPjLJzPysr09640a5aldWVvuyZU6ffz46NW1hYuK/hv90hiScCg8vgaDja9Y8R5NUv4XHj+u6ukSJicKbN5sOHcocP/6muTmLx5MKhfK4uPqgILfeXolQqODxkJkzLR8+bOntFfT1SaKjq5OSWlAtTZmiJZPBXV0yExMmDCMymaymphaGYRKJNGKEKwRBCoVYV5eUmsppbubr6tIQBDCZRDc3vZoaIYmERxAAw4i5OSs5uZVAwGlrU7S0xNnZnQyG3ZYty6VSKCKi18yMmZ/f8vRpVkUFx8PjOoIgUqlSKlXev1/D4fDlchiCgFwONzTwpk69c/RoDgCgqqo3ICCmt1cyeXJUPw0EqgAAG55JREFUaioHfCAZ598dTGC/m6SkloyMNgCAtjbV0pKtqUkhEHBSqRIMMi67u3u0tPx/PWXk38boISbmDRtG6erSPvvMUUvrN1Zy/Pys4uKcb9yYHRTktmePZ0HBCgIB19UlMjCgP3pURyafHBiQrljhlJfHXbduZGGhpK9PwWaDoqKmxMRmGEZwOGTDhgSxGD5zpnDnzmQYRjZseHTlSq1CgfT1ady5Uy0WK3JzW4OD3UUixdWrpTt3JuvpUUUixYYNo3R0CNeulR89mn3u3HSlEhEK5TKZsr29b8UKawcH07i4pqtXyyIjfalUnEKBcLkMc3Oby5enjB7NhGFk1qxrCAKzWCSFAm5p4efkdPD5cjabQibjk5NbqFTCxInG9fV9CgW8bp3r5culDQ39jY0Df/+M8+8OJrDfjaUle+PGxJMn83bvTnFx0fH3t3N21u7tFZ84kbtjh0daWqvKE2J49r+gILchKWwHVzp/C9nZ2UQiUU9PT9UybZp5VtYyFxed8+enGxrST52a8sMPExYvdtixI7m0dDUMbystXTVhArx/v4mnp+EXX7iuWsVsa1s3MCC1s+tJSRkfH7/8/v2lDAbJ09N42zYPExMmh8Pz9rbw9jZ1cdG9evVjJyftysremhr+/v0G48cboSkMwsKmjx5NZ7HwIhFEpZqHh8+srV3j7W3a0iKGYejzz9np6RyJRHn3LnfPnqmXL0+XySCpFC4rS5HJbgOA1NXx2tsFEASNHKnn42MqFMqTkzk6OjQuN9DBQUtPj1Zbu9bdXe8t+Sc/PP5aIyYGanwf3NLeLvjuu9QHD2qWLYtFEOTevSo3t/MAjDc1/TE2tg5BELlceehQ5q1bFUFBCahnrbFx2IoVcZcvl/j53adQTigUv5rFFQpFXl6+tfXp1avP379/PyoqhkwOLi/nDP50lWtybGzswMBAaipnzpz/qD199+5dpVKJIIhUKn316lVMzHMzs7CIiFLk3765v/zSCMMIm/1zeXm5v3/EZ589QRCEz5dlZrYFB+fMm/egqKjI0NBIX/+EQCCLiChNTGxSGeg5HL6LyxUEQXJy2tF7CQnJO3gwA0EQW9uLLS0Daun09whmpv/bofJvVGUd7OuDJ00ienhoX7lS6utrpXL8MzFhHjiQDgAQCGQkEm70aP3IyDKlErl4sfjJk3q0qt38+bZ1dZK6uvSYGCOp1GTdOsvg4OyWFomrq467uwGaaWffvpRDh7LmzWM/f54dEODQ2yvetOlFXFz91q1jAgNHCYW4fftSHB3Zly/nLFxouWlTsb4+ffv25E2bXkyaZPrkSW1j4xfJyS39/dKzZ9tqaiS9vY12dpc0Ncnr1o2sqemTyZTJyUqlcs+tW3PEYsV336W5uemxWOR582zmz7d9bcb5P7AQ/7cFE9jfC7kcHl41k8GgJCTEh4Y+zMkRgNc5/kkkyitXSmtreQcOTFi/Pl5Li1pR0dPVJWpp4U+cSMfjIQcHrTlz3O/c4cfEdPj6GlRW9np6kvF4iVisiIy8b2ZGwuGgFSt0rax0Kys5Eol02zbnuLj6nTuTAwLsW1u72Gwqnd7f0QH5+48/cqROqUQ6OtYnJ7csXhzDYpHNzFgyGeztbXr69NRvvyV0d7fo6uqyWBpr146YN8/GwuLC9euzV6xwYrPJlZW9ISFTUK/IFy+aJk82HZ4x8u1umR8c2BzsL2ZwfgHU+P5aZ0UikWhpacnn88F/Ov6hb6dQ8Gw2OT+f6+amL5fDtrYadDpBKJSZm+O1tUUQBKZMGT9/vo+1tZalpVZLC2hrk1pY2I0fP5lIxE2ePMXdfRQEASaTQSIRRSIRhQJXV5ccOqSrowOFhkbu3bshMvK2i4snnU5RKGAcDrDZJAgCLBaps1MEAAIA0Nen9fSIS0q6mUzKnDmj09MbFQop+LextK6uv7CwEwBw8mRuU9OAtbWGiQnD2lrjfyH7LyawvxLU+H7hQnFoaIHK+D4kdEU1iKqpgTs6ZPn5Ndu2jS0v71m79vm9e9U4HGhs7KfRiHp6dDabsnx5bFjYdAoF39kpGDkS0GiUyZMnIwig0+loIFlfn4RGI0AQaG4eMDZmwjCCx+OZTCYAwNDQ0NHR0c7OTktLa+rUaf/X3n3HNXWvfwA/QBKGYQTDCMoUZCoCKrRyRVAUqaJ1t7XYauvA4qRX0Xq1FleL4ijF8RMU0EolaMERAaksZe8hARMIM4QEAoZA1rl/nDaXRu3Q33HU5/1XTnJCFp/XWd/v8xgbu+nq6m3YsNjObpWRkZ2Z2Shl4CUSOYPB3rr1l4iIfyGI2vCw/Pz5aien0Rs2ZKxa5SSTabS0ICdP3h8/PurmTZZUquDzB7FZXuHhnhkZLQEByStW3FiwwPbHH+uVs7xe8S+BGzUUSrC+OXp6egoKCgICAgiE/+3b377NXrDg2sSJRsHBTqdPVxUXL/v3v2/ExravXTuBwxH5+lps3pxlbKw1b57dlSv1Y8fqLVs2/siRIhsb/YkTja5eZS5bNr63V5KR0ezrSxMIZLNnWz18KHBxoR4+XFhc/LGhoSQ0NOP+/X59fYKjI41IVGcwmlEUJZE0jhyZXlTUFR9f6+Bg6OREzclplUjkjo6j16xxOXKkmMvNlUrFdnZzT5yY6ednYWBwqq8v9PDhQolEYW9POXCgYN06V19fi6Cga42Nn6GoavGFf45XeorlzZaX12ZoeCogIHnXrpzAQPq5c1V/6+lNTb137rD/7ovev3+/pqbmyfsjIh4MDkoNDE7s3/9jRMSdR496IyOLFy68hqLo+fPXExIqJBL5xIkXurpE6IjRuvr6J1EULSzkGBsf7uzsEgjEmprHsLoD2EMbNyaHhqZu3x5va/vDli0/1dTwiMSjt2+zs7NbtbSObdv2C4qigYHJWlpRiYm1paVdmZktoaGZjo6x9fX1R45ccnc/LZVKURQ1MDiJoujChdcePepTvueHD/leXol/9xt4s8Au4vObNm2MszN14ULbAwf+lZgYuHFjZnOz8K8//fZt9qlT5X/3RSdNmtTU1DQ4qHoFViKRFBTkzplDKSoiSKVkGxsD5f26unKFgvDUK0vY7MZbt250d3Nu3EijULRcXY2xbl3YQ6WlvLVrJwcE0PLzFy9aZNrY2KCurqZQKKZPH2tnR7l8uV4sltnYGCAIWlHBGzOGzOMNWlrqq6urOTg4+Pn5EQiEnJwcmUz2p7O8/qkgYC9EeY0YRRENDTVsGNRf5OZmPHKRTmeqqUWuXZt+4kTZkSNF7u7xPN5TxjFoa2vb2tpWVVWNvJPFYtXVNdFotIMH5zMYnKlTTZHfupbx+fysLGlYWPZTp3hiSzNmzDA3twgKCkJRZHBQ6uhoiD0kEAiMjTXj45nGxsZtbV0oapmdzZFKFcnJzN7eoeBg58FBiZ7eicbGXpkMPXOmctKki3y+2NRUh8f79cw7n48ODRF37rwmFA6npj7audNz5IX4u3c52Cyvv/6lvXHgNP2LystrF4mkaWmPvv56mpkZeXhYHhNTIZHI79/vSEgIzM1ti4goWLbM/uLF+x0dBV9+uT4/nx8QYI2Vf1O54mRkpD08LNfTIzGZglWrXJTFSVU4ODgwGAwej2dkZCQSiYqLi0tKhIWF0rY27VmzDEJCJgUEWLPZwoyM5ro6fl4e29XVODt7UHnWJCjIFpv4SCaTsNmNlpZGEoluenoPj9d6/Lifjg4R696QlFQZGuqyd29jVhZ7/HhCfPyHZWXjyOSu4GBrCkVr0yZ3NltIpzfSaOQ7d5b4+Vns2pV77lxVcfHHH33khPx25h1FUW3t+zExbf7+FtraxPz8/1VlCwmZFBIy6VlfbGeniEZ78xvzvep91DfbjBlXsDmX331XZGf3fyxW36FDBV9/ff/KlfoJE+K+/75MKBw2MYlWKNDw8ENqantLS8vb2gbc3eNRFM3La8OOQDicfjL5OI83qKt7Yu3aOyiK8niD2tpRAoH4Wa/L4XCOHDmSnp5+/fr1+vp6hULxrDXv3bvX3t7xx7XQamp43t6Xsds3bjxSU/vus88YEREPvL3P3rlTh6KoTCaj0+nY0VR3d/fPP//c3//rGAsHh/M3bjzCbsvlCkvLM08ei8rlckPD47dvZ7BYrD/5Qn9TX89fsSLtL678OoMt2IvS19e0s6OEhU25eZN14kRpS0v/yGotIpEUK/Cye/f22Ngzw8PikQVeqFQdBEHMzXVtbSkdHY9lMkVNTU9kZHFlJc/QUKuhoffWLZa1tX5JSdecOdZBQeO43MGUFKahoVZ+PufUKamW1t1Fi6aVlz+YMaN1x46psbE1jo7U9evT+/pCExPrkpIanJ1H19YyN282nz07RaEIQxCkuVnY0NA7Z47VyI8gFEoGBiTY7ffes8Hq5k+ZMloobNq6taC21lFDQ4NCofT09JiamhoZGU2YMCEvL2/mzJkkEmnk31G2kpk3zyYmpsLNzYROZyYkBKqrq5NIpIiIiMIC7iS3nR995JKZ2eLpSaPRyDdvPsK25ypb/nv3Wuvq+Glpj+bPH4f/b4gjOAZ7Idh13uZmIYPR3Ns7RKXqjDyOj42tPniwAEWR/Px2c/MzQqHs+PGHQUHJYrFM5Y/I5QoHB0MCQV1dXS0sbMqWLR5yOfrLLxyVyWO5uW2FhZ2zZ5u5uAxqaxP9/a0vXZofHz9XKkXnzqW7uBgtXmw3dqyusofDzp2u69ePmTXLRnlJNyOjJTKyGEGQiIiCzZuzNm3KsrQ8W1nZvXr1hOLiLmwdbNe0s7PT3NyIw+nHPqOxsbGy7pq1tTWNRrt48WJRUZHKF/JkKxnliPh169Z9tNKntVX01IJtKnXalFPscPnZXiLYgj2/nJy22tqe69ebsrI4eXntfn4WW7Z4CIXDymot4eGeNFqMuroalao9Zgy5vp5PIhmIRD1C4bCJyQ85OR+oqSHx8bUdHY/j4uaSSBo6OgQqVRsr53Lu3JwrV+pV6pPNnz/u5s2GCRMufP/9NDU1tpOTNfbo4sV26enNDAbby4smkymUPRzKykqJRAIyoiegu7vJTz81NDcL4+Nrmcw1CIKMHUueNcvSzo4y8qPJ5WhubtPVq92bNrnLZIpTpyp6e/szM5sYDIfBQVlKClNHB/nyy2QE2U8k7iwo6Lxy5WFCQmBiYl1a2iMfn7HHj5e6u5t8+eW9/n7Jjh3ZixaNRxBk2bJlixYtcXCIfWrBNpU6bYWFnS/zp8QPBOz5TZ8+ls//AkGQ0lKuQJCTkBCIIAiZ/LvjeH19zebmtZqaGlSq9v790y5cqD1wwCw4+PHnn09UV0dSU98f+QdRFAkNdffzs8AWmUxBbm7btm2TlfXJUlPLFi5EFy/2O3iwikBQy83lIghy4wYrPr7W0FCruprHYvU1NAhQFJVI5B0dj+fM/kZXz62lZQaKIgcPFtrY6OfntyO/lRVZvjwtIMC6qakvP79DJWB0eoNU2rlv3ww/Pyvs6rCjo+nFi9VxcVUmJqMYjLpVq0aFhMylUj+MjBRPnmyamtrE5Q4ymb0SieLatfejo8t++qnBx2cskai+bdvkd9+9/GTtJ/T35zOxLf+33/oMDcnKyro1NTWwKXZvOgjY85PJFN9+W6z8r0UQROVAQllXFFNZyWOxhHv2KOztW3x9p4hE0r1785UHKgkJdT094oMHC2pq+MuX26ekME1MdPr6JD/8UJGS0rhgge2aNSli8aCzs7WPDzkkZFJ4eFZf37Czc5yrq5GeHkmhQKuqesrLu4lE9ZiYitpafk/P0DRvHxbLYM+ePKFw+NGjvl27PLEB+BSKlp+fxdWrDcbGOrt3ezk7x82daxUdXYEd74nFcj8/ann5QFRUeXo6Jzm54dKlebdusfz8ZPHx2V99Zc/lileu5H3xhduFC2wuV3TxYnVb2+MPPkjDKhFkZbWMG2dw7x6HyxXJZIqIiAfGxjoczoCyYFtDg4DNFqoUbFOp0yaTKbApdlu2eLzRgzwgYM/vyWkjUVEl2FCgxMS6+PjajRvdRq7v4WFiZkY+caIURamnTiUePXrY0lK/tLQrL6996dK05cvtP/7YycmJ2tsrnjXrJ01NjZSUBe3toujoCipVs6WldcUK2gcfVLm7k2JiKhISAnfvznJ2Nly0yCk/v53HE3O5oshIHxZLSCSq9/UNc7mi6upPbGz0GQz2vHkpVlZ6w8Py99+/7udnib2Z0aO1QkLcGAx2bm6bubnu4cNFJiY6n37qMn/+uJiYyuZm7rRpVleulKWmvv/4sfT06crMzBYC4Syn5X2RyC8pabq39+W8vPbp083lcvThw14DA834+MDIyJKaGv7s2VZHj5bMnGlBJGqkp8+xsjprajpKKAzFNmKPH29GEMTe3lAi2YYgiKur0eDgFuwtjdzyEwjqFRWrXsrPiC84yfH87t1rtbbWR0ZMG8Fqcixf7lBV9YlKuhAE0dfXjIryzc5eoampk5fnTKFo9fYOWVjojR6tXVvb4+lJk0oV+/blz51r/d57NiSShr//VRpNk8Xi29gMhYXZNzSMwg5dsNMAAwNyAuHXcvBtbQM6OkQeb0hHh0Akanz99bTSUu7x4yUIgsyYYW5sPMrDw1RPj7R16+R9+/JlMpTJ7HV3N2lvH8DKiggEQ52dj7Gr5MXFXQiCXr3arK2tj1WrDg/3bGjgk0jqHh57bO3Mu7sVJ0+WEQjqISGT7O0NOZz+4GCnCROo+/c/uHSpTi5XVFf3lJdzfXzMzczIKSmNo0YRw8KmvO7lQXEDAXt+T04b+YOhQHI5ip1P8/Yea2tLkcnUNmzIEIkkWloEgUBMJKqXlnJjYvxptFFr12b09g7v3ftubGzAV19lDQ0JCgu7p06dUlTUNWoUcdMm96qqT+ztDfv6pGx2/8mTZfHxtcPDsrt3l1VWcltbB4TC4aKiTmtr/ezsNmzMhI4OYcmS8XV1/KioEoFgCEUVdDpz8mRTMzMynd5IIKivW+c6dSoNKxYSEGBtZaW7Y4c5jUZRfszlyx19fS2Skz9CUc2lS8dv2zaZTCatXOn0zjs0Ly+zd94xo9MXlJd3JyQE9vR8oatLdHc3uX+/Izp61uefT9TX1/z4Y6eX/NO8PmAX8fkp60VTKFrYtBGVAwmsM11GRrOenmZtbY9AMCQQiIVCiZubCYXSP26cwenTlW5uJhMnGuXmtrHZfceOdS9aNJ5IVDc3133woANFka4uBYVi4OXlcPZsVWNjb3f34LVrjUuWjNfRIS5ZYkwmG27a5F5Y2NnU1KenR8Iq6qBo2MCAtLm5XyD4Attu0OnM/n7Jzp1TDx0qPHDgX7t2eW7YkHn4cMqePe/ExVVraKjt2/dubGw1nz+0Zs2dy5frCYS7O3bourkpamv5uron6fQgiUR2/ny1iwvVwEDT2/vHmTMtGhoEcXHVfP5QXR3f1zfp0KHpUqlcJJLGxlbzeOKtWz3epOrWuHrVV7rfIunpzR4e8ceOlYSH52zcmCmXK27fZrm4xEVFlSxYcG3Xrhweb9DS8szChdeCg2+lpDBRFHV1vaCm9l1yMhNF0ba2gTFjTnt5Je7YkS2VypOSbk2blrhmDSMs7N7MmUlsdt/wsCwkJOPMmcrDhwtTU5uUr2tre27fvnxLyzOnT1dg9yi7+DEYbGzABI83iKLoN988WL36tqen5/r127HFNWsY2FOwwfUtLUJb23MKBZqZ2bJkyc8jHwJPBVuwl8ff39Lf/3dtnVWmxz85o37pUvvAQJuvvsoNCLAaM4a8YYPr7t1e2HMNDDTS0gIplN+dXo+OnoUgCIoiOTmtI+/38jKbOdPyvffoPj7mT21HQqVq19T0xMVVMxgBzc3f+Pv7Y4tlZcHYCs/qxoSMuMgGngTHYK+Rp86o/+YbbxOTUaGhd5WrdXWJtm+/s3r1sXnzLiFPa6Kp7BOJUShQuVzh7T1m82aPpUtTxWIZNrsZQZC+viHs6FEuRz/9lBEdPYvPbyeRTJSL+vqa2BycZ3VjevIhMBIE7DWyd++7YrEsPDw3N7d95UonbC6Mhobajz/Ou3WLdelSPbaaSCS1tCR1czNZLDFWdkqlQLeyTySCIOnpzR0dImyCyX/+8w6ZTPzwwxvGxjoUilZISCaT2SsUDnd3D0ZGFkulcqFQnJjYkJ8vwhaxY6qkpAZscP2tW6wni0ApH3pVX9rr7lXvo4I/glUIRFE0K6vFwODkJ5/cRn/rmN7Y2I11TFeWFhzZgBybQfx3MZlMZW8k8P8CtmCvNZFIit3w9bXYvn1KQUEH8luPTFtbI5UJlOgL77ax2Wxra+sXesfg9+Akx+vr7l3OpUv1fn4Ws2ZZIgiye7cXny9GEESlNOfEiUYqu20PHwqw3TZlr639+x/U1vZMmWJKJGpwOP1RUSU5OR94e48Z+XICgUAqlY6szg1eHFSVeivExlavXj0Buz13Ln3cOIPvv5+psk5paamOjo6jo+NLf3f/ZLCL+FYIDnbGbsTF1dTV8Uc2aMa0tbUVFRVZWlq+9Lf2Dwe7iG8FAkEdQZCOjsfbtv1y5cp8MpmI/H7sf3PzoYY6cVISZdkyB2UrV4lEHhFRMHJK9av+HG8e2EV8iwQFXaNStWNjA7DFkWVAXV3FRka8tDQzJvOzrCzO6dMVV68GHT1aIpXKd+707OkRW1mdbW1dR6H8SR8zoAK2YG+LxMS6kpKuurrV2KJEIleZRNze/jgjI1mlH6fKlGoPDwjY3wPHYG8FLndw8+asH37wNzD49QJ0QkLds8b+o8/ox/nUMVbgj8EW7K2wcWMmmUxqbx+Iji6XydDc3LYpU0xVxv7T6cwn+3HCoPgXBMdgAOAIdhEBwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwBEEDAAcQcAAwNF/AZcv6YKJjZ07AAAAAElFTkSuQmCC)
The pre-Serrata graph is on the top, the post-serrata on the bottom.
Before the Serrata, there were small (up to 4-5) sets of connected
components, the biggest one the “star” around the Tiepolo family (which
was an important part of the early XIV century rebellion); however, the
Serrata had the effect of connecting the network, mainly because there
were many more doges per generation (as explained in (Merelo 2023)), but also because marrying became
an tactical move to increase the reach of the social network.
It should be noted, however, that doges’ marriages were but a small
part of the complete Venetian social network; however, it is an
important sample that, besides, shows who arrived to the summit of power
(or servitude to the state, depending on how you look at it).
Social networks pre and post-Serrata
This will call the function that creates, from raw data, the
igraph
data structure with which we can work. We still need to simplify them to be able to work with it.The pre-Serrata graph is on the top, the post-serrata on the bottom. Before the Serrata, there were small (up to 4-5) sets of connected components, the biggest one the “star” around the Tiepolo family (which was an important part of the early XIV century rebellion); however, the Serrata had the effect of connecting the network, mainly because there were many more doges per generation (as explained in (Merelo 2023)), but also because marrying became an tactical move to increase the reach of the social network.
It should be noted, however, that doges’ marriages were but a small part of the complete Venetian social network; however, it is an important sample that, besides, shows who arrived to the summit of power (or servitude to the state, depending on how you look at it).