{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Zip code distribution analysis for \"geochunk\" algorithm\n",
    "\n",
    "Our goal is to group zip codes by variable length prefixes in such a way that the chunk sizes tend to be roughly equal, with no giant, leftover chunks.\n",
    "\n",
    "But first, we need to import some standard Python data libraries."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We have some zip code census data from [this blog post](https://blog.splitwise.com/2013/09/18/the-2010-us-census-population-by-zip-code-totally-free/) that we want to load and examine. We need to do this `np.object_` business to force the CSV parser to treat zip codes as strings, and not as integers (thus losing the leading zeroes)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style>\n",
       "    .dataframe thead tr:only-child th {\n",
       "        text-align: right;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: left;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Zip Code ZCTA</th>\n",
       "      <th>2010 Census Population</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>01001</td>\n",
       "      <td>16769</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>01002</td>\n",
       "      <td>29049</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>01003</td>\n",
       "      <td>10372</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>01005</td>\n",
       "      <td>5079</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>01007</td>\n",
       "      <td>14649</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>01008</td>\n",
       "      <td>1263</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>01009</td>\n",
       "      <td>741</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>01010</td>\n",
       "      <td>3609</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>01011</td>\n",
       "      <td>1370</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>01012</td>\n",
       "      <td>661</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "  Zip Code ZCTA  2010 Census Population\n",
       "0         01001                   16769\n",
       "1         01002                   29049\n",
       "2         01003                   10372\n",
       "3         01005                    5079\n",
       "4         01007                   14649\n",
       "5         01008                    1263\n",
       "6         01009                     741\n",
       "7         01010                    3609\n",
       "8         01011                    1370\n",
       "9         01012                     661"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "zippop = pd.read_csv('../src/zip2010.csv', dtype={'Zip Code ZCTA': np.object_, '2010 Census Population': np.int32})\n",
    "zippop[:10]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "308739931"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pop = zippop['2010 Census Population']\n",
    "pop.sum()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's look at the distribution of population by zip code."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.text.Text at 0x7fead3553fd0>"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAEWCAYAAACufwpNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmYHVW57/Hvz0RmJAlpOJkkUQLXwEGEFoKoF0GTMGjw\nHPQEB8JwzUXxKE4QnIIKCh6PCFcFIwTCIAEZJCKYEyaVcyHQYYYQaUIgHQNpSBickMB7/qi1odLs\nvXunU3vv7M7v8zz76aq3VlWt1QV5e1WtvUoRgZmZWRHe0OwKmJlZ/+GkYmZmhXFSMTOzwjipmJlZ\nYZxUzMysME4qZmZWGCcV2yhI2k9S13rsf46kbxRZJ8tICkk7NrseVgwnFWs4SUsl/U3SnyU9JekC\nSVs1u14lko6UdGs+FhHHRsR3mlWnvpD05vQ77vlZI+mmVOZkSRfn9glJf0nllkv6oaQBzWuFtRon\nFWuWD0bEVsAeQDvw9SbXp+VJGphfj4gnImKr/Ad4F/A34LtVDvX2VPYA4GPAp+pWaet3nFSsqSJi\nOXA9sCuApOGS5kpaJalT0qv/oKW/qq+QdJmkFyTdJentue1r3UZJPaBTyp1X0nRJj6bjPCTpwyn+\nNuAcYJ/01/qz5Y4l6VOpfqtSfYf3qMexkh6R9Kykn0hShXr01qbhkq6U1C3pMUmfK7PvxZKeB46s\n9ruW9CbgCuD0iLihWlmAiHgY+APp2pQ53i6S5qffwVOSvprim0r6kaQ/pc+PJG2a2+8rklakbUf3\nOOamkn4g6Yl0zHMkbd5bXW3D4aRiTSVpFHAQcHcKzQG6gOHAYcB3Je2f22Uy8EtgCPAL4FeS3tiH\nUz8KvAfYBvgWcLGkYRGxCDgWuC39dT+oTJ33B74HfBQYBjye6p13CPBOYLdUbmKVupRtk6Q3AL8G\n7gVGkPUcjpc0sce+VwCDgEt6afP5wCPAqb2UA0DSOLLf0d1ltm0N3AD8luxa7QjcmDZ/DRgP7A68\nHdiL1BOVNAn4MvABYCzw/h6HPg3YKe27I1m7v1lLfW0DERH++NPQD7AU+DPwLNk/yD8FNgdGAS8D\nW+fKfg+4IC2fDNye2/YGYAXwnrQewI657RcAp6Tl/YCuKnW6B5iclo8Ebu2xPX+s84Dv57ZtBbwE\njM7V49257ZcD0yuct2KbgL2BJ3qUPwk4P7fv72v8nX8p/d6HlDn/xbn1AJ4HVpMl3lOAN5Q53uHA\n3RXO9ShwUG59IrA0Lc8CTstt26l03QABfwHemtu+D/BYs/+b9af2z1r3YM0a6NDocQsm3UJaFREv\n5MKPkz1zKVlWWoiIV9KIruGsI0lHAF8ERqfQVsDQGncfDtyVq8efJT1D9lf10hR+Mlf+r+n4lVRq\nUwDDS7fgkgFkt6Ret28lkt5N1hvbLyJW9VYe2CMiOnspM4oseZQznOy6lTzOa9doOLCwx7aSNmAL\nYGHubqHI2mwtwknFNiR/AoZI2jqXWN4MLM+VGVVaSLeHRqb9IPvHe4tc2X8iu5W2Fkk7AD8nu510\nW0S8LOkesn/AIPvHvLd67pA73pbAtj3quS4qtWkN2V/pY6vsW7WukrYHLgO+HBEdfaxfOcuAKRW2\nlX4/D6b1N/PaNVpBrr1pW8nTZIMIdonsWZu1ID9TsQ1GRCwD/j/wPUmbSdoNOAa4OFdsT0n/kkY6\nHQ+8CNyett0DfEzSgHTv/n9XONWWZP8YdwNIOoq1H0Y/BYyUtEmF/S8FjpK0e3oA/V1gQUQsXbcW\n99qmO4AXJJ0oafPUrl0lvbOWg6ahwHOAmyLinD7WrZJrgWGSjk8P17eWtHfadinwdUltkoaSPRMp\nXcPLgSMljZO0BTCjdMCIeIUs2Z8habvUhhE9niHZBs5JxTY0h5PdkvoTcDUwo8dtsmuAfyO75/9J\n4F8i4qW07fPAB8me1Xwc+FW5E0TEQ8B/AreRJZB/Bv47V+Qmsr+yn5T0dJn9bwC+AVxJ9pf3W6n8\nV3styrYpIl4me+C/O/AY2V/y55INLqjFvmTPkv5Vr/+uyoO97FtV6kl+gOz3/STZAID3pc2nAB3A\nfcD9ZLcKT0n7XQ/8iOx33Jl+5p2Y4renEW03ADuvT12tsRThl3RZa5B0MtmD+E80uy5F6Y9tso2b\neypmZlYYJxUzMyuMb3+ZmVlh6tZTkTRL0kpJD/SI/7ukhyU9KOn7ufhJadqLxfnRHpImpVinpOm5\n+BhJC1L8siojdczMrEHq1lOR9F6yb01fGBGleZ3eRzaFw8ER8aKk7SJiZZoO4lKy6RyGk4342Ckd\n6o9ko0y6gDuBwyPiIUmXA1dFxBxJ5wD3RsTZvdVr6NChMXr06ELbambW3y1cuPDpiGjrrVzdvvwY\nEb+XNLpH+NNkUzS8mMqsTPHJwJwUf0xSJ1mCAeiMiCUAkuYAkyUtAvYnm0EVYDbZdBO9JpXRo0fT\n0VHkd8DMzPo/SY/3XqrxD+p3At6Tblv9LvclrhGsPd1EV4pVim8LPBsRa3rEy5I0TVKHpI7u7u6C\nmmJmZj01OqkMJJuJdTzwFeByqfyU4EWKiJkR0R4R7W1tvfbezMysjxo991cX2XOQAO6Q9ArZJH7L\nWXs+oJG8No9SufgzwCBJA1NvJV/ezMyapNE9lV+RpnKQtBOwCdnUE3OBKWkOoTFk71m4g+zB/Ng0\n0msTsqkw5qakdDPZ+zYAppJNdWFmZk1Ut56KpEvJ5h0amqbynkH2LoVZaZjxP4CpKUE8mEZzPUQ2\nM+txad4jJH0WmEc2/fWsiCjNWXQiMEfZ2/juJnvHhZmZNdFG9+XH9vb28OgvM7N1I2lhRLT3Vs7T\ntJiZWWGcVMzMrDBOKmZmVhi/TngdjJ7+m7Lxpacd3OCamJltmNxTMTOzwjipmJlZYZxUzMysME4q\nZmZWGCcVMzMrjJOKmZkVxknFzMwK46RiZmaFcVIxM7PCOKmYmVlhnFTMzKwwTipmZlYYJxUzMyuM\nk4qZmRWmbklF0ixJK9P76Htu+5KkkDQ0rUvSWZI6Jd0naY9c2amSHkmfqbn4npLuT/ucJUn1aouZ\nmdWmnj2VC4BJPYOSRgETgCdy4QOBsekzDTg7lR0CzAD2BvYCZkganPY5G/hUbr/XncvMzBqrbkkl\nIn4PrCqz6QzgBCByscnAhZG5HRgkaRgwEZgfEasiYjUwH5iUtr0pIm6PiAAuBA6tV1vMzKw2DX2m\nImkysDwi7u2xaQSwLLfelWLV4l1l4pXOO01Sh6SO7u7u9WiBmZlV07CkImkL4KvANxt1zpKImBkR\n7RHR3tbW1ujTm5ltNBrZU3krMAa4V9JSYCRwl6R/ApYDo3JlR6ZYtfjIMnEzM2uihiWViLg/IraL\niNERMZrsltUeEfEkMBc4Io0CGw88FxErgHnABEmD0wP6CcC8tO15SePTqK8jgGsa1RYzMyuvnkOK\nLwVuA3aW1CXpmCrFrwOWAJ3Az4HPAETEKuA7wJ3p8+0UI5U5N+3zKHB9PdphZma1G1ivA0fE4b1s\nH51bDuC4CuVmAbPKxDuAXdevlmZmViR/o97MzArjpGJmZoVxUjEzs8I4qZiZWWGcVMzMrDBOKmZm\nVhgnFTMzK4yTipmZFcZJxczMCuOkYmZmhXFSMTOzwjipmJlZYZxUzMysME4qZmZWGCcVMzMrjJOK\nmZkVxknFzMwK46RiZmaFqec76mdJWinpgVzsPyQ9LOk+SVdLGpTbdpKkTkmLJU3MxSelWKek6bn4\nGEkLUvwySZvUqy1mZlabevZULgAm9YjNB3aNiN2APwInAUgaB0wBdkn7/FTSAEkDgJ8ABwLjgMNT\nWYDTgTMiYkdgNXBMHdtiZmY1qFtSiYjfA6t6xP4rItak1duBkWl5MjAnIl6MiMeATmCv9OmMiCUR\n8Q9gDjBZkoD9gSvS/rOBQ+vVFjMzq00zn6kcDVyflkcAy3LbulKsUnxb4NlcgirFy5I0TVKHpI7u\n7u6Cqm9mZj01JalI+hqwBrikEeeLiJkR0R4R7W1tbY04pZnZRmlgo08o6UjgEOCAiIgUXg6MyhUb\nmWJUiD8DDJI0MPVW8uXNzKxJGtpTkTQJOAH4UET8NbdpLjBF0qaSxgBjgTuAO4GxaaTXJmQP8+em\nZHQzcFjafypwTaPaYWZm5dVzSPGlwG3AzpK6JB0D/BjYGpgv6R5J5wBExIPA5cBDwG+B4yLi5dQL\n+SwwD1gEXJ7KApwIfFFSJ9kzlvPq1RYzM6tN3W5/RcThZcIV/+GPiFOBU8vErwOuKxNfQjY6zMzM\nNhD+Rr2ZmRXGScXMzArjpGJmZoVxUjEzs8I4qZiZWWGcVMzMrDBOKmZmVhgnFTMzK4yTipmZFcZJ\nxczMCuOkYmZmhXFSMTOzwjipmJlZYXpNKpI+L+lNypwn6S5JExpROTMzay219FSOjojngQnAYOCT\nwGl1rZWZmbWkWpKK0s+DgIvSS7JUpbyZmW2kakkqCyX9F1lSmSdpa+CV+lbLzMxaUS1vfjwG2B1Y\nEhF/lbQtcFR9q2VmZq2olp5KAOOAz6X1LYHNettJ0ixJKyU9kIsNkTRf0iPp5+AUl6SzJHVKuk/S\nHrl9pqbyj0iamovvKen+tM9ZknxLzsysyWpJKj8F9gFK75x/AfhJDftdAEzqEZsO3BgRY4Eb0zrA\ngcDY9JkGnA1ZEgJmAHuTvY9+RikRpTKfyu3X81xmZtZgtSSVvSPiOODvABGxGtikt50i4vfAqh7h\nycDstDwbODQXvzAytwODJA0DJgLzI2JVOu98YFLa9qaIuD0iArgwdywzM2uSWpLKS5IGkN0GQ1Ib\nfX9Qv31ErEjLTwLbp+URwLJcua4UqxbvKhMvS9I0SR2SOrq7u/tYdTMz600tSeUs4GpgO0mnArcC\n313fE6ceRqzvcWo818yIaI+I9ra2tkac0sxso9Tr6K+IuETSQuAAsu+nHBoRi/p4vqckDYuIFekW\n1soUXw6MypUbmWLLgf16xG9J8ZFlypuZWRNV7KmkkVpD0sPylcClwC/IEsOQPp5vLlAawTUVuCYX\nPyKNAhsPPJduk80DJkganB7QTwDmpW3PSxqfRn0dkTuWmZk1SbWeykKy21MC3gysTsuDgCeAMdUO\nLOlSsl7GUEldZKO4TgMul3QM8Djw0VT8OrIvV3YCfyV9DyYiVkn6DnBnKvftiCg9/P8M2QizzYHr\n08fMzJqoYlKJiDEAkn4OXB0R16X1A6lhpFVEHF5h0wFlygZwXIXjzAJmlYl3ALv2Vg8zM2ucWh7U\njy8lFICIuB54V/2qZGZmraqWaVr+JOnrwMVp/ePAn+pXJTMza1W19FQOB9rIhhVfDWzHa9+uNzMz\ne1UtQ4pXAZ9PsxNHRPy5/tUyM7NWVMubH/9Z0t3AA8CDkhZK8gNyMzN7nVpuf/0M+GJE7BAROwBf\nAmbWt1pmZtaKakkqW0bEzaWViLiFbPp7MzOztdQy+muJpG8AF6X1TwBL6lclMzNrVbX0VI4mG/11\nFXAlMDTFzMzM1lLL6K/VvPbWRzMzs4pqGf01X9Kg3PpgSfPqWy0zM2tFtdz+GhoRz5ZWUs9lu/pV\nyczMWlUtSeUVSW8urUjagQa9XMvMzFpLLaO/vgbcKul3ZFPfvweYVtdamZlZS6rlQf1vJe0BjE+h\n4yPi6fpWy8zMWlEtPRVSErm2znUxM7MWV8szFTMzs5o4qZiZWWEqJhVJm0k6XtKPJf1fSTXdKquF\npC9IelDSA5IuTecaI2mBpE5Jl0naJJXdNK13pu2jc8c5KcUXS5pYVP3MzKxvqvVUZgPtwP3AgcB/\nFnFCSSPIvqHfHhG7AgOAKcDpwBkRsSOwGjgm7XIMsDrFz0jlkDQu7bcLMAn4qaQBRdTRzMz6plpS\nGRcRn4iInwGHkQ0lLspAYPPU+9kCWAHsD1yRts8GDk3Lk9M6afsBkpTicyLixYh4DOgE9iqwjmZm\nto6qJZWXSgsRsaaoE0bEcuAHwBNkyeQ5YCHwbO48XcCItDwCWJarx3PAtvl4mX3WImmapA5JHd3d\n3UU1xczMeqiWVN4u6fn0eQHYrbQs6fm+nlDSYLJexhhgONm7WSb19Xi1iIiZEdEeEe1tbW31PJWZ\n2Uat4sP3iKjX84n3A49FRDeApKuAfYFBkgam3shIYHkqvxwYBXSl22XbAM/k4iX5fczMrAmqjf4a\nUu2zHud8AhgvaYv0bOQA4CHgZrJnNwBTgWvS8ty0Ttp+U0REik9Jo8PGAGOBO9ajXmZmtp6qDRNe\nSDZxpMpsC+AtfTlhRCyQdAVwF7AGuJvsnfe/AeZIOiXFzku7nAdcJKkTWEU24ouIeFDS5WQJaQ1w\nXES83Jc6mZlZMard/hpTr5NGxAxgRo/wEsqM3oqIvwMfqXCcU4FTC6+gmZn1SbXbX9+TVPaptqTT\n61clMzNrVdVGf30G6JB0cJltH6hTfczMrIVVSypLgIOB76WpWjbNbSv3nMXMzDZy1ZJKRMQDZM85\nXgEWStqtMdUyM7NWVMtLuv4OfE7SgcCvJZ1V/2qZmVkrqtZTWesWV0RcD7wTeB/gHouZmb1OtZ7K\nB3sGImIlcIik8WXKm5nZRq5iTyUiuqpsu70+1TEzs1bmNz+amVlhqn35cd/0c9NKZczMzPKq9VRK\no7xua0RFzMys9VV7UP+SpJnAiHLDiCPic/WrlpmZtaJqSeUQsnefTCSbsdjMzKyqarMUP002Ff2i\niLi3gXUyM7MWVcvor2ckXS1pZfpcKWlk3WtmZmYtp5akcj7ZWxaHp8+vU8zMzGwttSSV7SLi/IhY\nkz4XAGXfs2JmZhu3WpLK05I+IWlA+nwCeKbeFTMzs9ZTS1I5Gvgo8CSwAjgMOGp9TippkKQrJD0s\naZGkfSQNkTRf0iPp5+BUVpLOktQp6T5Je+SOMzWVf0TS1PWpk5mZrb9ek0pEPB4RH4qItojYLiIO\njYgn1vO8ZwK/jYj/BbwdWARMB26MiLHAjWkd4EBgbPpMA84GkDSE7D33e5O982VGKRGZmVlzNHzu\nL0nbAO8FzgOIiH9ExLPAZGB2KjYbODQtTwYujMztwCBJw8i+PzM/IlZFxGpgPjCpgU0xM7MemjGh\n5BigGzhf0t2SzpW0JbB9RKxIZZ4Etk/LI4Bluf27UqxS/HUkTZPUIamju7u7wKaYmVleM5LKQGAP\n4OyIeAfwF1671QVk7zEGoqgTRsTMiGiPiPa2Ng9cMzOrl16TiqSv55aLmLG4C+iKiAVp/QqyJPNU\nuq1F+rkybV8OjMrtPzLFKsXNzKxJqk19f6KkfchGe5Ws94zFEfEksEzSzil0APAQ2RcsSyO4pgLX\npOW5wBFpFNh44Ll0m2weMEHS4PSAfkKKmZlZk1SbUPJh4CPAWyT9Ia1vK2nniFi8nuf9d+ASSZsA\nS8iGKL8BuFzSMcDjZMOYAa4DDgI6gb+mskTEKknfAe5M5b4dEavWs15mZrYeqiWVZ4GvAvulz9vI\negPTU2J5V19PGhH3AO1lNh1QpmwAx1U4zixgVl/rYWZmxaqWVCYC3wTeCvwQuA/4S0Ss1xcfzcys\n/6r4TCUivhoRBwBLgYuAAUCbpFsl/bpB9TMzsxZSradSMi8iOoAOSZ+OiHdLGlrvipmZWeupZZqW\nE3KrR6bY0/WqkJmZta51+vKj3wBpZmbVNOMb9WZm1k85qZiZWWGcVMzMrDBOKmZmVhgnFTMzK4yT\nipmZFcZJxczMCuOkYmZmhXFSMTOzwjipmJlZYZxUzMysME4qZmZWGCcVMzMrTNOSiqQBku6WdG1a\nHyNpgaROSZel99cjadO03pm2j84d46QUXyxpYnNaYmZmJbW8pKtePg8sAt6U1k8HzoiIOZLOAY4B\nzk4/V0fEjpKmpHL/JmkcMAXYBRgO3CBpp4h4udENGT39N2XjS087uME1MTNrrqb0VCSNBA4Gzk3r\nAvYHrkhFZgOHpuXJaZ20/YBUfjIwJyJejIjHgE5gr8a0wMzMymnW7a8fAScAr6T1bYFnI2JNWu8C\nRqTlEcAygLT9uVT+1XiZfdYiaZqkDkkd3d3dRbbDzMxyGp5UJB0CrIyIhY06Z0TMjIj2iGhva2tr\n1GnNzDY6zXimsi/wIUkHAZuRPVM5ExgkaWDqjYwElqfyy4FRQJekgcA2wDO5eEl+HzMza4KG91Qi\n4qSIGBkRo8ketN8UER8HbgYOS8WmAtek5blpnbT9poiIFJ+SRoeNAcYCdzSoGWZmVkYzR3/1dCIw\nR9IpwN3AeSl+HnCRpE5gFVkiIiIelHQ58BCwBjiuGSO/zMzsNU1NKhFxC3BLWl5CmdFbEfF34CMV\n9j8VOLV+NTQzs3Xhb9SbmVlhnFTMzKwwTipmZlYYJxUzMyuMk4qZmRXGScXMzArjpGJmZoVxUjEz\ns8I4qZiZWWGcVMzMrDBOKmZmVhgnFTMzK4yTipmZFcZJxczMCuOkYmZmhXFSMTOzwjipmJlZYZxU\nzMysMA1/nbCkUcCFwPZAADMj4kxJQ4DLgNHAUuCjEbFakoAzgYOAvwJHRsRd6VhTga+nQ58SEbMb\n2ZbejJ7+m7Lxpacd3OCamJk1RjN6KmuAL0XEOGA8cJykccB04MaIGAvcmNYBDgTGps804GyAlIRm\nAHuTvdt+hqTBjWyImZmtreFJJSJWlHoaEfECsAgYAUwGSj2N2cChaXkycGFkbgcGSRoGTATmR8Sq\niFgNzAcmNbApZmbWQ1OfqUgaDbwDWABsHxEr0qYnyW6PQZZwluV260qxSvFy55kmqUNSR3d3d2H1\nNzOztTUtqUjaCrgSOD4ins9vi4gge95SiIiYGRHtEdHe1tZW1GHNzKyHpiQVSW8kSyiXRMRVKfxU\nuq1F+rkyxZcDo3K7j0yxSnEzM2uShieVNJrrPGBRRPwwt2kuMDUtTwWuycWPUGY88Fy6TTYPmCBp\ncHpAPyHFzMysSRo+pBjYF/gkcL+ke1Lsq8BpwOWSjgEeBz6atl1HNpy4k2xI8VEAEbFK0neAO1O5\nb0fEqsY0wczMyml4UomIWwFV2HxAmfIBHFfhWLOAWcXVzszM1oe/UW9mZoVpxu2vjZ6/aW9m/ZV7\nKmZmVhgnFTMzK4yTipmZFcZJxczMCuMH9RsQP8A3s1bnnoqZmRXGPZUW4B6MmbUKJ5UWVinZgBOO\nmTWHk0o/5d6NmTWDk8pGxsnGzOrJScUAJxszK4aTilVV7blNOU5CZhs3JxUrlJOQ2cbN31MxM7PC\nuKdiTbWuPZtK3OMx2zA4qVi/4IEGZhuGlk8qkiYBZwIDgHMj4rQmV8k2IH3pCTkRmfVdSycVSQOA\nnwAfALqAOyXNjYiHmlsza2VF3ZKrxEnL+rOWTirAXkBnRCwBkDQHmAw4qdgGq95Jqy8qJTrfVrR1\n1epJZQSwLLfeBezds5CkacC0tPpnSYv7eL6hwNN93HdD1R/bBP2zXXVrk06vb/le+Fq1hh1qKdTq\nSaUmETETmLm+x5HUERHtBVRpg9Ef2wT9s139sU3QP9vVH9tUq1b/nspyYFRufWSKmZlZE7R6UrkT\nGCtpjKRNgCnA3CbXycxso9XSt78iYo2kzwLzyIYUz4qIB+t4yvW+hbYB6o9tgv7Zrv7YJuif7eqP\nbaqJIqLZdTAzs36i1W9/mZnZBsRJxczMCuOkUgNJkyQtltQpaXqz69OTpFGSbpb0kKQHJX0+xYdI\nmi/pkfRzcIpL0lmpPfdJ2iN3rKmp/COSpubie0q6P+1zliQ1sH0DJN0t6dq0PkbSglSXy9IgDSRt\nmtY70/bRuWOclOKLJU3MxRt+bSUNknSFpIclLZK0T3+4VpK+kP77e0DSpZI2a7VrJWmWpJWSHsjF\n6n5tKp2jJUWEP1U+ZAMAHgXeAmwC3AuMa3a9etRxGLBHWt4a+CMwDvg+MD3FpwOnp+WDgOsBAeOB\nBSk+BFiSfg5Oy4PTtjtSWaV9D2xg+74I/AK4Nq1fDkxJy+cAn07LnwHOSctTgMvS8rh03TYFxqTr\nOaBZ1xaYDfyftLwJMKjVrxXZF5EfAzbPXaMjW+1aAe8F9gAeyMXqfm0qnaMVP02vwIb+AfYB5uXW\nTwJOana9eqnzNWTzoS0GhqXYMGBxWv4ZcHiu/OK0/XDgZ7n4z1JsGPBwLr5WuTq3ZSRwI7A/cG36\nn/FpYGDP60M2CnCftDwwlVPPa1Yq14xrC2yT/vFVj3hLXytem91iSPrdXwtMbMVrBYxm7aRS92tT\n6Ryt+PHtr96VmwpmRJPq0qt0G+EdwAJg+4hYkTY9CWyfliu1qVq8q0y8EX4EnAC8kta3BZ6NiDVl\n6vJq/dP251L5dW1vPY0BuoHz0y29cyVtSYtfq4hYDvwAeAJYQfa7X0hrX6uSRlybSudoOU4q/Yik\nrYArgeMj4vn8tsj+BGqp8eOSDgFWRsTCZtelQAPJbq+cHRHvAP5CdrvjVS16rQaTTeY6BhgObAlM\namql6qAR16YVr3+ek0rvWmIqGElvJEsol0TEVSn8lKRhafswYGWKV2pTtfjIMvF62xf4kKSlwByy\nW2BnAoMklb64m6/Lq/VP27cBnmHd21tPXUBXRCxI61eQJZlWv1bvBx6LiO6IeAm4iuz6tfK1KmnE\ntal0jpbjpNK7DX4qmDSC5DxgUUT8MLdpLlAaeTKV7FlLKX5EGr0yHngudb3nARMkDU5/eU4gu4+9\nAnhe0vh0riNyx6qbiDgpIkZGxGiy3/tNEfFx4GbgsArtKrX3sFQ+UnxKGnE0BhhL9sC04dc2Ip4E\nlknaOYUOIHtVQ0tfK7LbXuMlbZHOW2pXy16rnEZcm0rnaD3NfqjTCh+yUR5/JBt98rVm16dM/d5N\n1l2+D7gnfQ4iu0d9I/AIcAMwJJUX2cvNHgXuB9pzxzoa6Eyfo3LxduCBtM+P6fGguQFt3I/XRn+9\nhewfmk7gl8CmKb5ZWu9M29+S2/9rqe6LyY2Gasa1BXYHOtL1+hXZCKGWv1bAt4CH07kvIhvB1VLX\nCriU7JnQS2S9ymMacW0qnaMVP56mxczMCuPbX2ZmVhgnFTMzK4yTipmZFcZJxczMCuOkYmZmhXFS\nMauBpJc+7CIZAAACSklEQVQl3ZNm4P2lpC0KPv6Rkn7cS5n9JL0rt36spCOKrIfZ+nJSMavN3yJi\n94jYFfgHcGwT6rAf8GpSiYhzIuLCJtTDrCInFbN19wdgRwBJX0y9lwckHZ9io5W9K+USZe9LuaLU\ns5G0VNLQtNwu6ZaeB5f0QWXvGLlb0g2Stk8ThR4LfCH1mN4j6WRJX0777C7p9vRej6v12js/bpF0\nuqQ7JP1R0nvq/+uxjZmTitk6SPNUHQjcL2lP4Chgb7J3ZHxK0jtS0Z2Bn0bE24Dnyd4fUqtbgfGR\nTTg5BzghIpaSvY/kjNRj+kOPfS4EToyI3ci+3T0jt21gROwFHN8jblY4JxWz2mwu6R6y6VWeIJtr\n7d3A1RHxl4j4M9kkiqWewLKI+O+0fHEqW6uRwDxJ9wNfAXapVljSNsCgiPhdCs0me9lUSWmC0YVk\n7woxq5uBvRcxM9IzlXxA1d/S23P+o9L6Gl77Y26zCvv+P+CHETFX0n7AyetU09d7Mf18Gf8/b3Xm\nnopZ3/0BODTNzLsl8OEUA3izpH3S8sfIbmkBLAX2TMv/WuG42/DalOhTc/EXyF4XvZaIeA5YnXte\n8kngdz3LmTWCk4pZH0XEXcAFZLPsLgDOjYi70+bFwHGSFpHNQnx2in8LOFNSB1nPoZyTgV9KWkj2\nmt2SXwMfLj2o77HPVOA/JN1HNgvyt9enbWZ95VmKzQqWRmpdm4Yfm21U3FMxM7PCuKdiZmaFcU/F\nzMwK46RiZmaFcVIxM7PCOKmYmVlhnFTMzKww/wMvrP33ydaj+QAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fead30cdef0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(pop, bins=50)\n",
    "plt.title(\"Population per ZIP code\")\n",
    "plt.xlabel(\"Population\")\n",
    "plt.ylabel(\"# of ZIP codes\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style>\n",
       "    .dataframe thead tr:only-child th {\n",
       "        text-align: right;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: left;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Zip Code ZCTA</th>\n",
       "      <th>2010 Census Population</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>20640</th>\n",
       "      <td>60629</td>\n",
       "      <td>113916</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27727</th>\n",
       "      <td>79936</td>\n",
       "      <td>111086</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2748</th>\n",
       "      <td>11368</td>\n",
       "      <td>109931</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30109</th>\n",
       "      <td>90650</td>\n",
       "      <td>105549</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29986</th>\n",
       "      <td>90011</td>\n",
       "      <td>103892</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30192</th>\n",
       "      <td>91331</td>\n",
       "      <td>103689</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2720</th>\n",
       "      <td>11226</td>\n",
       "      <td>101572</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30045</th>\n",
       "      <td>90201</td>\n",
       "      <td>101279</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2753</th>\n",
       "      <td>11373</td>\n",
       "      <td>100820</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2714</th>\n",
       "      <td>11220</td>\n",
       "      <td>99598</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2759</th>\n",
       "      <td>11385</td>\n",
       "      <td>98592</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2527</th>\n",
       "      <td>10467</td>\n",
       "      <td>97060</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30460</th>\n",
       "      <td>92335</td>\n",
       "      <td>95397</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26671</th>\n",
       "      <td>77084</td>\n",
       "      <td>95137</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2451</th>\n",
       "      <td>10025</td>\n",
       "      <td>94600</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2702</th>\n",
       "      <td>11208</td>\n",
       "      <td>94469</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30075</th>\n",
       "      <td>90280</td>\n",
       "      <td>94396</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26773</th>\n",
       "      <td>77449</td>\n",
       "      <td>94382</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2729</th>\n",
       "      <td>11236</td>\n",
       "      <td>93877</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27192</th>\n",
       "      <td>78521</td>\n",
       "      <td>93818</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "      Zip Code ZCTA  2010 Census Population\n",
       "20640         60629                  113916\n",
       "27727         79936                  111086\n",
       "2748          11368                  109931\n",
       "30109         90650                  105549\n",
       "29986         90011                  103892\n",
       "30192         91331                  103689\n",
       "2720          11226                  101572\n",
       "30045         90201                  101279\n",
       "2753          11373                  100820\n",
       "2714          11220                   99598\n",
       "2759          11385                   98592\n",
       "2527          10467                   97060\n",
       "30460         92335                   95397\n",
       "26671         77084                   95137\n",
       "2451          10025                   94600\n",
       "2702          11208                   94469\n",
       "30075         90280                   94396\n",
       "26773         77449                   94382\n",
       "2729          11236                   93877\n",
       "27192         78521                   93818"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "zippop.nlargest(20, '2010 Census Population')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style>\n",
       "    .dataframe thead tr:only-child th {\n",
       "        text-align: right;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: left;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Zip Code ZCTA</th>\n",
       "      <th>2010 Census Population</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>75</th>\n",
       "      <td>01199</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>383</th>\n",
       "      <td>02203</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1423</th>\n",
       "      <td>05481</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1558</th>\n",
       "      <td>06020</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1580</th>\n",
       "      <td>06059</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "     Zip Code ZCTA  2010 Census Population\n",
       "75           01199                       0\n",
       "383          02203                       0\n",
       "1423         05481                       0\n",
       "1558         06020                       0\n",
       "1580         06059                       0"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "zippop.nsmallest(5, '2010 Census Population')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now let's look at 4-digit zip code prefixes:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style>\n",
       "    .dataframe thead tr:only-child th {\n",
       "        text-align: right;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: left;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Zip Prefix</th>\n",
       "      <th>2010 Census Population</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>533</th>\n",
       "      <td>1121</td>\n",
       "      <td>710148</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>534</th>\n",
       "      <td>1122</td>\n",
       "      <td>621986</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>498</th>\n",
       "      <td>1046</td>\n",
       "      <td>605938</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3773</th>\n",
       "      <td>6062</td>\n",
       "      <td>605797</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>535</th>\n",
       "      <td>1123</td>\n",
       "      <td>588379</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>532</th>\n",
       "      <td>1120</td>\n",
       "      <td>584187</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>497</th>\n",
       "      <td>1045</td>\n",
       "      <td>560238</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3772</th>\n",
       "      <td>6061</td>\n",
       "      <td>543914</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>481</th>\n",
       "      <td>1002</td>\n",
       "      <td>505290</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3774</th>\n",
       "      <td>6063</td>\n",
       "      <td>497137</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "     Zip Prefix  2010 Census Population\n",
       "533        1121                  710148\n",
       "534        1122                  621986\n",
       "498        1046                  605938\n",
       "3773       6062                  605797\n",
       "535        1123                  588379\n",
       "532        1120                  584187\n",
       "497        1045                  560238\n",
       "3772       6061                  543914\n",
       "481        1002                  505290\n",
       "3774       6063                  497137"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def sum_by_zip_prefix(zippop, length):\n",
    "    prefixes = zippop['Zip Code ZCTA'].str[:length].rename('Zip Prefix')\n",
    "    df = pd.concat([prefixes, zippop['2010 Census Population']], axis=1)\n",
    "    return df.groupby('Zip Prefix').sum().reset_index()\n",
    "prefix_4 = sum_by_zip_prefix(zippop, 4)\n",
    "prefix_4.nlargest(10, '2010 Census Population')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.text.Text at 0x7feacf8499e8>"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY8AAAEWCAYAAACe8xtsAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcHFW99/HP10S2gKxzMWRhAkauyMMNMGJ4AJ8oKqsC\nysXEhUUkNwov5fJcIYhKXFDc8BLxggEi8CD7IlHwYkBlkXUSYgiEyASCSQhJWIPAjQR+zx91OqlM\nunu6Zqane5Lv+/XqV1efOlX1656a/tU5dbpKEYGZmVkRb2t0AGZm1v84eZiZWWFOHmZmVpiTh5mZ\nFebkYWZmhTl5mJlZYU4e1mckjZG0qAfLXyjpG70ZU38gaZKkK9L0cEl/lzSghuVqrtvfKPNLSS9K\nelDS/pLmNTquDYmTxwZK0gJJr6cvl6WSLpW0eaPjKpF0nKR78mURMSEivtOomHpK0kaS5vYkgUbE\n3yJi84h4s2hdSX+S9IUq8V2Y9of84zVJIekDqU5IeleaniTpjVTvJUn3Stqnu++toP2AjwBDI2Lv\niLg7Inbpo20bTh4buo9FxObAnkAb8PUGx9PvSRpYZfZXgeV9FUtRKTlvnn8ANwJ/BP5cYbFrUr0W\n4B7gRknqSRxdfIYlOwILIuLVnmzLus/Jw4iIxcDvgN0AJO0gaZqkFyR1SDqxVDcdbV4v6RpJr0ia\nKelfcvNXH5mm15dK+m657UqaKGl+Ws9jko5M5e8BLgT2KR3VlluXpBNTfC+keHfoFMcESU+ko+Kf\nV/pSq+E97SDpBknLJT0l6ctllr1C0grguArbGAF8Fvh+ufmd60q6M8UyHdguN681vbeBubp3pbq3\np/d5Ree6ks4G9gfOT5/p+TXE8UXgg8C4rlo6EfEGcBnwTmDbMuvq6jNeIOl0SbOBV1PMZT93SScA\nF7Nm//iWcl2iknZO+8Se6fUOaR1junrPVjsnD0PSMOAQ4OFUdDWwCNgBOAr4nqQP5RY5HLgO2Aa4\nEvi1pLd3Y9Pzyb7QtgS+BVwhaXBEzAUmAPelI+CtysT8IbIv4qOBwcDTKe68w4D3AbunegdWiaXs\ne5L0NuA3wF+AIcABwCmSDuy07PXAVsCvKqz/Z8DXgNerxFByJTCDLGl8Bzi2i7oPkn1hTwI+V65S\nRJwJ3A2cnD7Tk6sFIOl9wI+AT0XE0q4ClrQxWeJcGBHPVajW1X4zDjiU7HN8iwqfe0Rcwtr7x1md\n3ut84HSy/Wkz4JfAZRHxp67eh9XOyWPD9ut0VH8PcCdZkhgG7AucHhH/ExGzyI7yjsktNyMirk9H\nm+cCmwCji248Iq6LiGci4q2IuAZ4Ati7xsU/A0yNiJkRsRI4g+xItDVX55yIeCki/kbW9TKqyvoq\nvaf3AS0R8e2I+EdEPAlcBIzNLXtfRPw6vY91kkNqUQ2IiJu6elOShqdtfiMiVkbEXWRfotXqfjPF\ndg8wratt1BDDNmRf8t9M66zm6LQPLQT2Ao6sUrer/WZyRCxMn2Etn3tFEXER0AE8QHZwcWYty1nt\naulbtPXXERFxe74gdf28EBGv5IqfJjsnUrKwNBERb6Xugh0oSNIxwKlAayranFwXTRd2AGbm4vi7\npOfJjlIXpOJnc/VfS+uvpNJ7CmCHUtdZMoDsKH6dZTuTNAj4IVnLrtz8C8m6swC+B/wBeLFTX/7T\nwLAyi5f+Vq91iqVc3Zqkrr0rgPaIOLeGRa6NiM92XW11bEDF/Sb/Oe5I1597Vy4iS6bj0wGG9SIn\nD+vsGWAbSVvkEshwYHGuzuovp9StMzQtB9mX9Ga5uu8k6wJbi6Qdyf65DyA7cn9T0iygdF6iq8s9\nP0P2BVNa3yCyrpvFFZeortJ7WgU8FREjqyxbLdaRZMnx7nTKZSNgS0nPAqMjYgJZF0xp2zsCW0sa\nlEsgwytsYwnZ32qzXAKpljhquYT214F3sfbBQm+ptt/A2vEtpOvPvSJlIwf/E7gEmCTphoh4oTvr\nsvLcbWVriYiFwL3A9yVtIml34ASyo9GSvSR9Ip20PQVYCdyf5s0CPi1pgKSDgP9TYVODyL4slgNI\nOp50wj5ZCgyVtFGF5a8Cjpc0KvW3fw94ICIWFHvHXb6nB4FX0sncTdP72i2dE6jFHLIvzVHp8YX0\n3kZRpsUSEU8D7cC3lA3t3Q/4WLkV5+pOSnX3qVQ3WQrsVGmmpA8DpwGfjIgVNby3oqrtN5319HM/\nj6z19AXgFrIBGNaLnDysnHFkR8vPADcBZ3Xq3roZ+BTwItkJ2k+kfmyAr5B9gb1Edl7i1+U2EBGP\nAT8B7iP7UvtfrD0c9A/Ao8CzktY5AZvi+QZwA9kR+M7U2B9eQdn3lEYZHUb2Zf8U8BzZOaAta1lp\nRKyKiGdLD+AF4K30utIIpk8D7091zwIur7KJzwD7AM8D3wWuIftSLuc84ChlP6ybXGb+14BNgfu0\n7u89PtPVe61Btf1mLT353CUdDhwEfDEVnQrs2UvvwRL5ZlBWhKRJwLsK9HM3vfXpPUm6Bni88wik\nRlufPmPLuOVh1o9Jel/6XcPbUjfh4VRo7Zn1Jp8wN+vf3kn2K/BtyQYmfDEiHq6+iFnPudvKzMwK\nc7eVmZkVtt52W2233XbR2tra6DDMzPqNGTNmPBcRLbXUXW+TR2trK+3t7Y0Ow8ys35D0dK113W1l\nZmaFOXmYmVlhTh5mZlaYk4eZmRXm5GFmZoU5eZiZWWFOHmZmVpiTh5mZFebkYWZmha23vzDvidaJ\nt5QtX3DOoX0ciZlZc3LLw8zMCnPyMDOzwpw8zMysMCcPMzMrrG7JQ9JUScskzcmVXSNpVnoskDQr\nlbdKej0378LcMntJekRSh6TJklSvmM3MrDb1HG11KXA+cHmpICI+VZqW9BPg5Vz9+RExqsx6LgBO\nBB4AbgUOAn5Xh3jNzKxGdWt5RMRdwAvl5qXWw9HAVdXWIWkw8I6IuD+ym61fDhzR27GamVkxjTrn\nsT+wNCKeyJWNkPSwpDsl7Z/KhgCLcnUWpbKyJI2X1C6pffny5b0ftZmZAY1LHuNYu9WxBBgeEXsA\npwJXSnpH0ZVGxJSIaIuItpaWmm7Da2Zm3dDnvzCXNBD4BLBXqSwiVgIr0/QMSfOBdwOLgaG5xYem\nMjMza6BGtDw+DDweEau7oyS1SBqQpncCRgJPRsQSYIWk0ek8yTHAzQ2I2czMcuo5VPcq4D5gF0mL\nJJ2QZo1l3RPlHwBmp6G71wMTIqJ0sv1LwMVABzAfj7QyM2u4unVbRcS4CuXHlSm7AbihQv12YLde\nDc7MzHrEvzA3M7PCnDzMzKwwJw8zMyvMycPMzApz8jAzs8KcPMzMrDAnDzMzK8zJw8zMCnPyMDOz\nwpw8zMysMCcPMzMrzMnDzMwKc/IwM7PCnDzMzKwwJw8zMyvMycPMzApz8jAzs8KcPMzMrDAnDzMz\nK6xuyUPSVEnLJM3JlU2StFjSrPQ4JDfvDEkdkuZJOjBXflAq65A0sV7xmplZ7erZ8rgUOKhM+U8j\nYlR63AogaVdgLPDetMx/SRogaQDwc+BgYFdgXKprZmYNNLBeK46IuyS11lj9cODqiFgJPCWpA9g7\nzeuIiCcBJF2d6j7Wy+GamVkBjTjncbKk2alba+tUNgRYmKuzKJVVKi9L0nhJ7ZLaly9f3ttxm5lZ\n0tfJ4wJgZ2AUsAT4SW+uPCKmRERbRLS1tLT05qrNzCynbt1W5UTE0tK0pIuA36aXi4FhuapDUxlV\nys3MrEH6tOUhaXDu5ZFAaSTWNGCspI0ljQBGAg8CDwEjJY2QtBHZSfVpfRmzmZmtq24tD0lXAWOA\n7SQtAs4CxkgaBQSwAPg3gIh4VNK1ZCfCVwEnRcSbaT0nA7cBA4CpEfFovWI2M7Pa1HO01bgyxZdU\nqX82cHaZ8luBW3sxNDMz6yH/wtzMzApz8jAzs8KcPMzMrDAnDzMzK8zJw8zMCnPyMDOzwpw8zMys\nMCcPMzMrzMnDzMwKc/IwM7PCnDzMzKwwJw8zMyvMycPMzApz8jAzs8KcPMzMrDAnDzMzK8zJw8zM\nCnPyMDOzwpw8zMysMCcPMzMrrG7JQ9JUScskzcmV/UjS45JmS7pJ0lapvFXS65JmpceFuWX2kvSI\npA5JkyWpXjGbmVlt6tnyuBQ4qFPZdGC3iNgd+CtwRm7e/IgYlR4TcuUXACcCI9Oj8zrNzKyP1S15\nRMRdwAudyn4fEavSy/uBodXWIWkw8I6IuD8iArgcOKIe8ZqZWe0aec7j88Dvcq9HSHpY0p2S9k9l\nQ4BFuTqLUllZksZLapfUvnz58t6P2MzMgAYlD0lnAquAX6WiJcDwiNgDOBW4UtI7iq43IqZERFtE\ntLW0tPRewGZmtpaBfb1BSccBhwEHpK4oImIlsDJNz5A0H3g3sJi1u7aGpjIzM2ugLlsekv5V0hZp\n+uuSbpS0Z3c2Jukg4DTg4xHxWq68RdKANL0T2YnxJyNiCbBC0ug0yuoY4ObubNvMzHpPLd1W34iI\nVyTtB3wYuIRsBFRVkq4C7gN2kbRI0gnA+cAWwPROQ3I/AMyWNAu4HpgQEaWT7V8CLgY6gPmsfZ7E\nzMwaoJZuqzfT86HAlIi4RdJ3u1ooIsaVKb6kQt0bgBsqzGsHdqshTjMz6yO1JI/Fkn4BfAT4gaSN\n2UB/md468Zay5QvOObSPIzEza6xaksDRwG3AgRHxErAN8NW6RmVmZk2ty+SRTmwvA/ZLRauAJ+oZ\nlJmZNbdaRludBZzOmkuJvB24op5BmZlZc6ul2+pI4OPAqwAR8QzZiCkzM9tA1ZI8/pF+zBcAkgbV\nNyQzM2t2tSSPa9Noq60knQjcTva7CzMz20B1OVQ3In4s6SPACmAX4JsRMb3ukZmZWdPqMnlI2jUl\ni+m5sjER8ad6BmZmZs2r1m6r05TZVNLPgO/XOzAzM2tetSSP9wPDgXuBh4BngH3rGZSZmTW3WpLH\nG8DrwKbAJsBTEfFWXaMyM7OmVkvyeIgsebwP2B8YJ+m6ukZlZmZNrZYLI56QrmwL2R3/Dpf0uTrG\nZGZmTa5i8pD0johYATwpaZtOs8tfXtbMzDYI1VoeV5LdLnYG2a/LlZsXwE51jMvMzJpYteRxTnp+\nT0T8T18EY2Zm/UO1E+bnped7+yIQMzPrP6q1PN6QNAUYKmly55kR8eX6hWVmZs2sWsvjMOAPZMN0\nZ5R5dEnSVEnLJM3JlW0jabqkJ9Lz1qlckiZL6pA0W9KeuWWOTfWfkHRs8bdpZma9qWLLIyKeA66W\nNDci/tLN9V8KnA9cniubCNwREedImphenw4cDIxMj/cDFwDvTyO9zgLayE7Uz5A0LSJe7GZMZmbW\nQ7X8SPB1SXeUWg+Sdpf09VpWHhF3AS90Kj4cuCxNXwYckSu/PDL3k10CfjBwIDA9Il5ICWM6cFAt\n2zczs/qoJXlcRHYL2jcAImI2MLYH29w+Ipak6WeB7dP0EGBhrt6iVFapfB2Sxktql9S+fPnyHoRo\nZmbV1JI8NouIBzuVreqNjefvUNhL65sSEW0R0dbS0tJbqzUzs05qSR7PSdqZNbehPYrsMiXdtTR1\nR5Gel6XyxcCwXL2hqaxSuZmZNUgtyeMk4BfAP0taDJwCTOjBNqcBpRFTxwI358qPSaOuRgMvp+6t\n24CPSto6jcz6aCozM7MGqXphRElvA9oi4sOSBgFvi4hXal25pKuAMcB2khaRjZo6h+wGUycATwNH\np+q3AocAHcBrwPEAEfGCpO+QXd0X4NsR0fkkvJmZ9aGqySMi3pJ0GnBtRLxadOURMa7CrAPK1A2y\nVk659UwFphbdvpmZ1Uct3Va3S/oPScPSD/y2KXOVXTMz24DUcj+PT6XnfKvAV9U1M9uAdZk8ImJE\nXwRiZmb9R8VuK0kjJd0saY6kqySV/WGemZlteKqd85gK/Bb4JDAT+FmfRGRmZk2vWrfVFhFxUZr+\nkaSZfRGQmZk1v2rJYxNJe7Dm9rOb5l9HhJOJmdkGqlryWAKcm3v9bO51AB+qV1BmZtbcqt3P44N9\nGYiZmfUftfxI0MzMbC1OHmZmVli133nsm5437rtwzMysP6jW8picnu/ri0DMzKz/qDba6g1JU4Ah\nkiZ3nhkRX65fWGZm1syqJY/DgA8DBwIz+iYcMzPrD6oN1X0OuFrS3Ij4Sx/GZGZmTa6W0VbPS7pJ\n0rL0uEHS0LpHZmZmTauW5PFLsvuL75Aev0llZma2gaolefxTRPwyIlalx6VAS53jMjOzJlZL8nhO\n0mclDUiPzwLPd3eDknaRNCv3WCHpFEmTJC3OlR+SW+YMSR2S5kk6sLvbNjOz3lHLbWg/T3Yvj5+S\nXRDxXuD47m4wIuYBowAkDQAWAzeldf40In6cry9pV2As8F6ybrPbJb07It7sbgxmZtYztdyG9mng\n43Xa/gHA/Ih4WlKlOocDV0fESuApSR3A3vjHi2ZmDdPoa1uNBa7KvT5Z0mxJUyVtncqGAAtzdRal\nMjMza5CGJQ9JG5G1aK5LRRcAO5N1aS0BftKNdY6X1C6pffny5b0Wq5mZra2RLY+DgZkRsRQgIpZG\nxJsR8RZwEVnXFGTnRIbllhuaytYREVMioi0i2lpaPCDMzKxeukwekr6em+7NK+yOI9dlJWlwbt6R\nwJw0PQ0YK2ljSSOAkcCDvRiHmZkVVO2S7KdL2gc4KlfcKyepJQ0CPgLcmCv+oaRHJM0GPgj8O0BE\nPApcCzwG/DdwkkdamZk1VrXRVo8D/wrsJOnu9HpbSbuk4bbdFhGvAtt2KvtclfpnA2f3ZJv11Drx\nlrLlC845tI8jMTPrG9W6rV4CvgZ0AGOA81L5REn31jkuMzNrYtVaHgcC3yQbAXUuMBt4NSK6/QNB\nMzNbP1RseUTE1yLiAGAB8P+AAUCLpHsk/aaP4jMzsyZUy+VJbouIdqBd0hcjYj9J29U7MDMza15d\nDtWNiNNyL49LZc/VKyAzM2t+hX4k6DsKmpkZNP7aVmZm1g85eZiZWWFOHmZmVpiTh5mZFebkYWZm\nhTl5mJlZYU4eZmZWmJOHmZkV5uRhZmaFOXmYmVlhTh5mZlaYk4eZmRXm5GFmZoU5eZiZWWENSx6S\nFkh6RNIsSe2pbBtJ0yU9kZ63TuWSNFlSh6TZkvZsVNxmZtb4lscHI2JURLSl1xOBOyJiJHBHeg1w\nMDAyPcYDF/R5pGZmtlqjk0dnhwOXpenLgCNy5ZdH5n5gK0mDGxGgmZk1NnkE8HtJMySNT2XbR8SS\nNP0ssH2aHgIszC27KJWtRdJ4Se2S2pcvX16vuM3MNngDG7jt/SJisaR/AqZLejw/MyJCUhRZYURM\nAaYAtLW1FVrWzMxq17DkERGL0/MySTcBewNLJQ2OiCWpW2pZqr4YGJZbfGgqa2qtE28pW77gnEP7\nOBIzs97VkG4rSYMkbVGaBj4KzAGmAcemascCN6fpacAxadTVaODlXPeWmZn1sUa1PLYHbpJUiuHK\niPhvSQ8B10o6AXgaODrVvxU4BOgAXgOO7/uQzcyspCHJIyKeBP6lTPnzwAFlygM4qQ9CMzOzGjTb\nUF0zM+sHnDzMzKwwJw8zMyvMycPMzApz8jAzs8KcPMzMrDAnDzMzK8zJw8zMCnPyMDOzwpw8zMys\nMCcPMzMrrJH389hg+VLtZtbfueVhZmaFOXmYmVlhTh5mZlaYk4eZmRXm5GFmZoU5eZiZWWFOHmZm\nVlifJw9JwyT9UdJjkh6V9JVUPknSYkmz0uOQ3DJnSOqQNE/SgX0ds5mZra0RPxJcBfzfiJgpaQtg\nhqTpad5PI+LH+cqSdgXGAu8FdgBul/TuiHizT6M2M7PV+rzlERFLImJmmn4FmAsMqbLI4cDVEbEy\nIp4COoC96x+pmZlV0tDLk0hqBfYAHgD2BU6WdAzQTtY6eZEssdyfW2wRFZKNpPHAeIDhw4fXLe56\n8WVLzKy/aNgJc0mbAzcAp0TECuACYGdgFLAE+EnRdUbElIhoi4i2lpaWXo3XzMzWaEjykPR2ssTx\nq4i4ESAilkbEmxHxFnARa7qmFgPDcosPTWVmZtYgjRhtJeASYG5EnJsrH5yrdiQwJ01PA8ZK2ljS\nCGAk8GBfxWtmZutqxDmPfYHPAY9ImpXKvgaMkzQKCGAB8G8AEfGopGuBx8hGap3kkVZmZo3V58kj\nIu4BVGbWrVWWORs4u25BmZlZIf6FuZmZFeY7CfYDHsJrZs3GLQ8zMyvMycPMzApzt1U/5u4sM2sU\ntzzMzKwwJw8zMyvM3VbroUrdWeAuLTPrHW55mJlZYU4eZmZWmJOHmZkV5nMeGxgP7zWz3uCWh5mZ\nFeaWhwFukZhZMU4eVpWTipmV424rMzMrzC0P6xa3SMw2bE4e1quq/bq9HCcbs/7JycMayi0Ys/7J\n5zzMzKywftPykHQQcB4wALg4Is5pcEhWR73V/eWLRJrVR79IHpIGAD8HPgIsAh6SNC0iHmtsZNYs\niiabass4qZh1rV8kD2BvoCMingSQdDVwOODkYb2uO4momTj5WV/oL8ljCLAw93oR8P7OlSSNB8an\nl3+XNK+b29sOeK6by/Y1x9r7+kucUCZW/aBBkXStX3+uTaw3Y92x1or9JXnUJCKmAFN6uh5J7RHR\n1gsh1Z1j7X39JU5wrPXiWLvWX0ZbLQaG5V4PTWVmZtYA/SV5PASMlDRC0kbAWGBag2MyM9tg9Ytu\nq4hYJelk4DayobpTI+LROm6yx11ffcix9r7+Eic41npxrF1QRDRiu2Zm1o/1l24rMzNrIk4eZmZW\nXET4kR7AQcA8oAOYWOdtTQWWAXNyZdsA04En0vPWqVzA5BTXbGDP3DLHpvpPAMfmyvcCHknLTGZN\nF2XZbVSJcxjwR7IfZD4KfKWJY90EeBD4S4r1W6l8BPBAWv81wEapfOP0uiPNb82t64xUPg84sKt9\npNI2atgPBgAPA79t5liBBelvNAtob9Z9IC2zFXA98DgwF9inGWMFdkmfZ+mxAjilGWMtG39PvgDX\npwfZP/F8YCdgI7IvoF3ruL0PAHuydvL4YemfHJgI/CBNHwL8Lu08o4EHcjvAk+l56zRd2tEeTHWV\nlj242jaqxDm4tJMCWwB/BXZt0lgFbJ6m3072BTkauBYYm8ovBL6Ypr8EXJimxwLXpOld099/Y7Iv\n2vlp/6i4j1TaRg37wanAlaxJHk0ZK1ny2K5TWdPtA6neZcAX0vRGZMmkKWPt9P3zLNmP9Jo61tUx\n9/RLcH15kB2d3JZ7fQZwRp232crayWMeMDhNDwbmpelfAOM61wPGAb/Ilf8ilQ0GHs+Vr65XaRsF\nYr6Z7BpjTR0rsBkwk+xKBM8BAzv/nclG7+2Tpgemeur8ty/Vq7SPpGXKbqOLGIcCdwAfAn5bbT1N\nEOsC1k0eTbcPAFsCT5GOsJs51k7xfRT4c3+ItfTwOY81yl0CZUgfx7B9RCxJ088C26fpSrFVK19U\nprzaNrokqRXYg+yIviljlTRA0iyyLsHpZEffL0XEqjLrXx1Tmv8ysG033sO2VbZRzX8CpwFvpdfV\n1tPoWAP4vaQZ6TJA0Jz7wAhgOfBLSQ9LuljSoCaNNW8scFUX62mWWAGfMG9akR0SRLNsQ9LmwA3A\nKRGxorvr6a5atxERb0bEKLKj+r2Bf65nXN0l6TBgWUTMaHQsNdovIvYEDgZOkvSB/Mwm2gcGknUH\nXxARewCvknXLFF1PjxT839oI+DhwXU/W013d3YaTxxrNcAmUpZIGA6TnZV3EVq18aJnyatuoSNLb\nyRLHryLixmaOtSQiXiI70b8PsJWk0g9i8+tfHVOavyXwfDfew/NVtlHJvsDHJS0AribrujqvSWMl\nIhan52XATWSJuRn3gUXAooh4IL2+niyZNGOsJQcDMyNiaRfraYZYV3PyWKMZLoEyjWzUBOn55lz5\nMcqMBl5OTc7bgI9K2lrS1mT9preleSskjZYk4JhO6yq3jbLS8pcAcyPi3CaPtUXSVml6U7JzM3PJ\nkshRFWItrf8o4A/pKGwaMFbSxpJGACPJTjyW3UfSMpW2UVZEnBERQyOiNa3nDxHxmWaMVdIgSVuU\npsn+dnNown0gIp4FFkraJRUdQDZSsOlizRnHmi6rautphljXKHqSZH1+kI1m+CtZP/mZdd7WVcAS\n4A2yo6UTyPqj7yAbPnc7sE2qK7KbYc0nG3bXllvP58mG4XUAx+fK28j+wecD57NmiF7ZbVSJcz+y\nJu1s1gwpPKRJY92dbNjr7LS+b6bynci+UDvIugY2TuWbpNcdaf5OuXWdmeKZRxqhUm0fqbSNGveF\nMawZbdV0sab6f2HNEOgzq/19GrkPpGVGAe1pP/g12QikZo11EFlrcMtcWVPG2vnhy5OYmVlh7rYy\nM7PCnDzMzKwwJw8zMyvMycPMzApz8jAzs8KcPMxyJL0paZakOZKuk7RZL6//OEnnd1FnjKT/nXs9\nQdIxvRmHWU85eZit7fWIGBURuwH/ACY0IIYxwOrkEREXRsTlDYjDrCInD7PK7gbeBSDp1NQamSPp\nlFTWKulxSb+SNFfS9aWWiqQFkrZL022S/tR55ZI+JukBZRfwu13S9souPjkB+PfUAtpf0iRJ/5GW\nGSXpfkmzJd2UflGMpD9J+oGkByX9VdL+9f94bEPm5GFWhrJrPx0MPCJpL+B4ssu7jwZOlLRHqroL\n8F8R8R6ym/l8qcBm7gFGR3YBv6uB0yJiAdk9Nn6aWkB3d1rmcuD0iNid7FfGZ+XmDYyIvcluKHQW\nZnXk5GG2tk2VXdK9Hfgb2XW99gNuiohXI+LvwI1A6ch+YUT8OU1fkerWaihwm6RHgK8C761WWdKW\nwFYRcWcquozspmIlpYtWziC7V4xZ3QzsuorZBuX1yC7pvlp2TbmKOl/fp/R6FWsOzjapsOzPgHMj\nYpqkMcCkQpGua2V6fhP/b1udueVh1rW7gSMkbZauKntkKgMYLmmfNP1psq4oyO68t1ea/mSF9W7J\nmktkH5srf4Xslr9riYiXgRdz5zM+B9zZuZ5ZX3DyMOtCRMwELiW7Eu0DwMUR8XCaPY/s5khzya7e\nekEq/xa7Vv4aAAAAaElEQVRwnqR2spZAOZOA6yTNILs1bMlvgCNLJ8w7LXMs8CNJs8muHvvtnrw3\ns+7yVXXNuimNjPptGtZrtkFxy8PMzApzy8PMzApzy8PMzApz8jAzs8KcPMzMrDAnDzMzK8zJw8zM\nCvv/3Fm3sH/MEGgAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fead34c2438>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(prefix_4['2010 Census Population'], bins=50)\n",
    "plt.title(\"Population per 4-digit ZIP prefix\")\n",
    "plt.xlabel(\"Population\")\n",
    "plt.ylabel(\"# of Prefixes\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "And various shorter prefix lengths. We tried to use 3-digit prefixes for the first Pachynorm run, which as you can, gave us some huge chunks with just under 3 million people."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style>\n",
       "    .dataframe thead tr:only-child th {\n",
       "        text-align: right;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: left;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Zip Prefix</th>\n",
       "      <th>2010 Census Population</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>699</th>\n",
       "      <td>770</td>\n",
       "      <td>2906700</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>555</th>\n",
       "      <td>606</td>\n",
       "      <td>2680484</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>91</th>\n",
       "      <td>112</td>\n",
       "      <td>2504700</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>803</th>\n",
       "      <td>900</td>\n",
       "      <td>2404395</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>843</th>\n",
       "      <td>945</td>\n",
       "      <td>2274102</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>276</th>\n",
       "      <td>300</td>\n",
       "      <td>2177710</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>679</th>\n",
       "      <td>750</td>\n",
       "      <td>2079278</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>818</th>\n",
       "      <td>917</td>\n",
       "      <td>1924074</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>306</th>\n",
       "      <td>331</td>\n",
       "      <td>1885845</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59</th>\n",
       "      <td>070</td>\n",
       "      <td>1683417</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    Zip Prefix  2010 Census Population\n",
       "699        770                 2906700\n",
       "555        606                 2680484\n",
       "91         112                 2504700\n",
       "803        900                 2404395\n",
       "843        945                 2274102\n",
       "276        300                 2177710\n",
       "679        750                 2079278\n",
       "818        917                 1924074\n",
       "306        331                 1885845\n",
       "59         070                 1683417"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "prefix_3 = sum_by_zip_prefix(zippop, 3)\n",
    "prefix_3.nlargest(10, '2010 Census Population')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.text.Text at 0x7feacf6cd3c8>"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAEWCAYAAACufwpNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHu9JREFUeJzt3XucZOO97/HP14w7McY0GXNJu8whtiNMOjIOyXHCdk+Q\niBBh3DLHDkdEshlyIfskws4Om9zHJUjEnSCX4xYSEpf0MBiGGAxmzJgeDIKN4Xf+WE9T06qqq6vW\nqurq/r5fr3rVqmc9a63f06u6fvU8a9VaigjMzMzysEKrAzAzs6HDScXMzHLjpGJmZrlxUjEzs9w4\nqZiZWW6cVMzMLDdOKpY7SdtLmt/A8j+T9M08Y2oHkg6WdHvJ639I2rDGZWuu224kfUfSEkmLJE1M\nbR3R6risPCeVIU7SPEmvpX/EZyWdL2mNVsfVq+8HKUBEHBER/7dVMdVL0lckPS7pJUnPSDpD0sh6\n1xcRa0TE4wOtm/bxd6rEeWJ6P5Q+XpEUkg5KdeZJ2jFNHyzprVTvJUmzJO1Rb7sGQtJE4KvAZhHx\n/oh4KrX1rWZs3wbOSWV4+GRErAFMBrqAb7Q4nrZXIVlcC0yOiPcBmwMfAo5uamA1iIhT0gfzOw/g\nDOAh4MoKi92R6o0CzgUuk7R2I3Eo099n0ETguYhY3Mi2rHmcVIaRiFgA/IHsAw9J60u6VtLzkuZK\n+mJvXUknS7pC0qWSXpZ0j6QPlcwPSRuXvK747VjSdEmPpfU8JGnvVP5B4GfANulb8NJy65L0xRTf\n8yne9fvEcYSkRyUtlfRjSaoQR39tWl/SlZJ6JD0h6egyy/5K0kvAwWX+vo9FxNLeRYC3gY371itZ\n5zqpPS9JuhvYqM/8d/7Gqe51qe7f0pDQ7X3rSpoGHAAcl/6m11Xafsmyu5Elv30i4pVqdSPibeA8\nYNW+8aZ1HSzpL5J+JOlFSQ9L2qFk/q2SvivpL8CrwIaS1pJ0rqSFkhakto1IPaUbgfVTW86X1Jna\nOlLSaEnzJX0yrXuN9D45qL82W3GcVIYRSROA3YB7U9ElwHxgfWAf4BRJnyhZZE/gcmA08GvgN5JW\nrGPTjwEfA9YCvg38StLYiJgDHEH6FhwRo8rE/Ange8C+wFjgyRR3qT2AjwBbpHo7V4mlbJvSN+br\ngPuAccAOwDGSdu6z7BVk39YvKrdySZ9PSWcJWU/l51Vi+THwX6ldh6ZHtbqvAO8HpqbHe0TEjBTb\nv6e/6SerrBNJncAvgS+m/VFV6qEdDvwDeLRCtY+S7fMxwEnAVZJGl8w/EJgGrEm2P88HlpEl4K2A\nnYDDI+ImYFfgmdSWg/u09Xmyv9nZktYl623NiogL+2uHFcdJZXj4TeoF3A78iSx5TAC2BY6PiP+K\niFnAOUDpt7yZEXFFRLwJnA6sAkwZ6MYj4vKIeCYi3o6IS8k+jLaucfEDgPMi4p6IeB04gaxn01lS\n59SIWBoRTwG3AFtWWV+lNn0E6IiIf4uIN9LxibOB/UqWvSMifpPa8VqFtv46DX/9N7Je2LPl6ik7\n0PwZ4FsR8UpEzAYu6KfuSRHxakQ8VKnuQEhamSxJXpT2SzVT0ntoEbA/sHdEvFih7mLgPyPizbTe\nR4DdS+afHxEPRsQysuS+G3BM+jssJksO+71nrWVExA1kXxJuTuv537UsZ8Wp+yCitZW90re+d6Qh\npOcj4uWS4ifJjrn0erp3IiLeVnZG1/oMUBqOOBboTEVrkH2LrcX6wD0lcfxD0nNkvYl5qXhRSf1X\n0/orqdSmIBtmWVpSdwRwW7ll+xMRj0p6EPgJ8GlJJwInptm/Ak4m+/8rXeeTFVbXUaZuzbFUcSbw\nJtmB8P7cGRHb1bjeBbH8lWqfZPn3TWnsHwBWBBaWjFquwMDaNwM4CjglIp4bwHJWACeV4esZYLSk\nNUsSy0RgQUmdCb0TaXhofFoOsg/v1Urqvp9sKG05kj5A9o1/B7Jv+m9JmkV2zAGyD/P+4vxAyfpW\nB9bpE+dAVGrTMuCJiJhUZdmBXtJ7JOm4Q0ScApxSsu0RaZsTgIdT8cQK6+lJdccDf+/bjnrilHQg\nWe9ncuq15WmcJJUklolkJzGUi+9p4HVgTOq5DEj6O84ALgS+JOkXETG3zrgtBx7+GqYi4mngr8D3\nJK0iaQvgMLJv0b0+LOnTaRz9GLJ//jvTvFnA59MB1V2A/1lhU6uTfYj0AEg6hHSiQPIsMF7SShWW\nvxg4RNKWabjmFOCuiJg3sBb326a7gZclHS9p1dSuzSV9pNYVSzo8je0jaTOyobqby9VNp8ReBZws\nabVUv9Jxkr51N2X5Ycq+ngUq/mZF0uZkPagD0vsgb+sCR6djVZ8FPgj8vlzFiFgI3AD8QNL7JK0g\naSNJld5PfZ1I9v46FPg+cKH8G5aWclIZ3vYnG5J6BriabMy+dJjsGuBzwAtkB1c/XfKt9svAJ4Gl\nZMc9flNuA2n8/wfAHWQfdv8d+EtJlT8CDwKLJC0ps/xNwDfJTnVdSPbNv6bx9grKtil9cO9Bdjzm\nCbID7eeQnVxQq22BByS9QvYh+nveHfIq5yiyobpFZAerf9FP3bVS3V+SJdvXK9Q9F9hM2dlw5fbL\nsWTJ/iq99/cq1eKt1V3AJLK/4XfJziqrNix1ELAS2SnNL5Ad5xnb30YkfZisLQel/XcaWYKZ3lD0\n1hD5Jl1WjqSTgY0j4gutjiUvQ6lNkk4D3h8RZXs3rSLpYLIzt2o9/mJDjHsqZm1A0qaStlBma7Kh\nyqtbHZdZXz5Qb9Ye1iQb8lqfbBjxB2RDeWaDioe/zMwsNx7+MjOz3LT18NeYMWOis7Oz1WGYmbWV\nmTNnLomIjiLW3dZJpbOzk+7u7laHYWbWViRVunpDwzz8ZWZmuXFSMTOz3DipmJlZbpxUzMwsN04q\nZmaWGycVMzPLjZOKmZnlprCkIuk8SYslzS4z76uSQtKY9FqSzpI0V9L9kiYXFZeZmRWnyJ7K+cAu\nfQvTvdF3Ap4qKd6V7P4Lk4BpwE8LjMvMzApS2C/qI+LPkjrLzDoDOI7lr7C6J3Bhuv3onZJGSRqb\n7grXdJ3Tf1e2fN6puzc5EjOz9tLUYyqS9gQWRMR9fWaNI7tXda/5qazcOqZJ6pbU3dPTU1CkZmZW\nj6YlFUmrkd1a9VuNrCciZkREV0R0dXQUcj00MzOrUzMvKLkRsAFwnySA8cA96S52C4AJJXXHpzIz\nM2sjTeupRMQDEbFuRHRGRCfZENfkiFgEXAsclM4CmwK82KrjKWZmVr8iTym+GLgD2ETSfEmHVan+\ne+BxYC5wNvClouIyM7PiFHn21/79zO8smQ7gyKJiMTOz5vAv6s3MLDdOKmZmlpu2vp1ws/lHkWZm\n1bmnYmZmuXFSMTOz3DipmJlZbpxUzMwsN04qZmaWGycVMzPLjZOKmZnlxknFzMxy46RiZma5cVIx\nM7PcOKmYmVlunFTMzCw3TipmZpYbJxUzM8uNk4qZmeXG91PJge+zYmaWcU/FzMxyU1hSkXSepMWS\nZpeUfV/Sw5Lul3S1pFEl806QNFfSI5J2LiouMzMrTpE9lfOBXfqU3QhsHhFbAH8HTgCQtBmwH/BP\naZmfSBpRYGxmZlaAwpJKRPwZeL5P2Q0RsSy9vBMYn6b3BC6JiNcj4glgLrB1UbGZmVkxWnlM5VDg\nD2l6HPB0ybz5qczMzNpIS5KKpK8Dy4CL6lh2mqRuSd09PT35B2dmZnVrelKRdDCwB3BAREQqXgBM\nKKk2PpW9R0TMiIiuiOjq6OgoNFYzMxuYpiYVSbsAxwGfiohXS2ZdC+wnaWVJGwCTgLubGZuZmTWu\nsB8/SroY2B4YI2k+cBLZ2V4rAzdKArgzIo6IiAclXQY8RDYsdmREvFVUbGZmVozCkkpE7F+m+Nwq\n9b8LfLeoeMzMrHj+Rb2ZmeXGScXMzHLjpGJmZrlxUjEzs9w4qZiZWW6cVMzMLDdOKmZmlhsnFTMz\ny42TipmZ5cZJxczMcuOkYmZmuXFSMTOz3DipmJlZbpxUzMwsN04qZmaWGycVMzPLjZOKmZnlxknF\nzMxy46RiZma5cVIxM7PcOKmYmVluCksqks6TtFjS7JKy0ZJulPRoel47lUvSWZLmSrpf0uSi4jIz\ns+IU2VM5H9ilT9l04OaImATcnF4D7ApMSo9pwE8LjMvMzApSWFKJiD8Dz/cp3hO4IE1fAOxVUn5h\nZO4ERkkaW1RsZmZWjGYfU1kvIham6UXAeml6HPB0Sb35qew9JE2T1C2pu6enp7hIzcxswFp2oD4i\nAog6lpsREV0R0dXR0VFAZGZmVq9mJ5Vne4e10vPiVL4AmFBSb3wqMzOzNtLspHItMDVNTwWuKSk/\nKJ0FNgV4sWSYzMzM2sTIolYs6WJge2CMpPnAScCpwGWSDgOeBPZN1X8P7AbMBV4FDikqLjMzK05h\nSSUi9q8wa4cydQM4sqhYzMysOfyLejMzy42TipmZ5abfpCLps5LWTNPfkHSVL6NiZmbl1NJT+WZE\nvCxpO2BH4Fx8GRUzMyujlqTyVnreHZgREb8DViouJDMza1e1JJUFkn4OfA74vaSVa1zOzMyGmVqS\nw77A9cDOEbEUGA38a6FRmZlZW+o3qUTEq2SXU9kuFS0DHi0yKDMza0+1nP11EnA8cEIqWhH4VZFB\nmZlZe6pl+Gtv4FPAKwAR8QywZpFBmZlZe6olqbxRepl6SasXG5KZmbWrWpLKZensr1GSvgjcBJxT\nbFhmZtaO+r2gZET8h6R/Bl4CNgG+FRE3Fh6ZmZm1nX6TiqTNUhK5saRs+4i4tcjAzMys/dQ6/HVc\nuoHWqpJ+CHyv6MDMzKz91JJUPgpMBP4K/A14Bti2yKDMzKw91ZJU3gReA1YFVgGeiIi3C43KzMza\nUi1J5W9kSeUjwMeA/SVdXmhUZmbWlmq5nfBhEdGdphcCe0o6sMCYzMysTVVMKpLeFxEvAY9LGt1n\n9u+KDcvMzNpRtZ7Kr4E9gJlkv6ZXybwANqx3o5K+Ahye1vMAcAgwFrgEWCdt88CIeKPebfSnc7rz\noplZ3qodUzk1PX8wIjaMiA1KHo0klHHA0UBXRGwOjAD2A04DzoiIjYEXgMPq3YaZmbVGtaRyZnr+\nawHbHQmsKmkksBrZsZpPAFek+RcAexWwXTMzK1C14a83Jc0Axks6q+/MiDi6ng1GxAJJ/wE8RXZW\n2Q1kw11LI2JZqjYfGFdueUnTgGkAEydOrCcEMzMrSLWeyh7AH8k++GeWedRF0trAnsAGwPrA6sAu\ntS4fETMioisiujo6OuoNw8zMClCxpxIRS4BLJM2JiPty3OaOZD+g7AGQdBXZL/RHSRqZeivjgQU5\nbtPMzJqglh8/vibpZkmzASRtIekbDWzzKWCKpNUkCdgBeAi4Bdgn1ZkKXNPANszMrAVqSSpnk91K\n+E2AiLif7GytukTEXWQH5O8hO514BWAG2S2Lj5U0l+y04nPr3YaZmbVGLb+oXy0i7s46Fe9YVqly\nLSLiJOCkPsWPA1s3sl4zM2utWnoqSyRtxLu3E96H7BRgMzOz5dTSUzmSbHhqU0kLgCeAAwqNyszM\n2lLVpCJpBbJfvu8oaXVghYh4uTmhmZlZu6k6/JXum3Jcmn7FCcXMzKqp5ZjKTZK+JmmCpNG9j8Ij\nMzOztlPLMZXPpecjS8oaukqxmZkNTf0mlYjYoBmBmJlZ+6s4/CVpkqRrJM2WdHG6ZL2ZmVlF1Y6p\nnAf8FvgM2a/ff9iUiMzMrG1VG/5aMyLOTtPfl3RPMwIyM7P2VS2prCJpK969jfCqpa8jwknGzMyW\nUy2pLAROL3m9qOR1kN2p0czM7B3V7qfyv5oZiJmZtb9afvxoZmZWEycVMzPLTbXfqWybnlduXjhm\nZtbOqvVUzkrPdzQjEDMza3/Vzv56U9IMYJyks/rOjIijiwvLzMzaUbWksgewI7AzMLM54QwtndN/\nV7Z83qm7NzkSM7PmqHZK8RLgEklzIuK+JsZkZmZtqpazv56TdLWkxelxpaTxhUdmZmZtp5ak8gvg\nWmD99LguldVN0ihJV0h6WNIcSdukm3/dKOnR9Lx2I9swM7PmqyWprBsRv4iIZelxPtDR4HbPBP5f\nRGwKfAiYA0wHbo6IScDN6bWZmbWRWpLKEklfkDQiPb4APFfvBiWtBXwcOBcgIt6IiKXAnsAFqdoF\nwF71bsPMzFqjlqRyKLAv2QUlFwL7AIc0sM0NgB7gF5LulXSOpNWB9SJiYaqzCFiv3MKSpknqltTd\n09PTQBhmZpa3fpNKRDwZEZ+KiI6IWDci9oqIpxrY5khgMvDTiNgKeIU+Q10REWRXQi4Xz4yI6IqI\nro6ORkfhzMwsT6249td8YH5E3JVeX0GWZJ6VNBYgPS9uQWxmZtaApieViFgEPC1pk1S0A/AQ2Rlm\nU1PZVOCaZsdmZmaNqfaL+iL9H+AiSSsBj5Mdo1kBuEzSYcCTZMdxzMysjfSbVCR9IyK+k6ZXjojX\nG91oRMwCusrM2qHRdZuZWetUu/T98ZK2ITvbq5evWGxmZhVV66k8DHwW2FDSben1OpI2iYhHmhLd\nEOULTZrZUFXtQP1S4ERgLrA92a/gAaZL+mvBcZmZWRuq1lPZGfgWsBFwOnA/8EpENPLDRzMzG8Iq\n9lQi4sSI2AGYB/wSGAF0SLpd0nVNis/MzNpILacUXx8R3UC3pH+JiO0kjSk6MDMzaz+1XKbluJKX\nB6eyJUUFZGZm7WtAv6j3HSDNzKyaVlz7y8zMhignFTMzy42TipmZ5aZVF5S0MvxLezNrd+6pmJlZ\nbpxUzMwsN04qZmaWGycVMzPLjZOKmZnlxknFzMxy46RiZma5cVIxM7PctCypSBoh6V5Jv02vN5B0\nl6S5ki6VtFKrYjMzs/q0sqfyZWBOyevTgDMiYmPgBeCwlkRlZmZ1a0lSkTQe2B04J70W8AngilTl\nAmCvVsRmZmb1a1VP5T+B44C30+t1gKURsSy9ng+MK7egpGmSuiV19/T0FB+pmZnVrOlJRdIewOKI\nmFnP8hExIyK6IqKro6Mj5+jMzKwRrbhK8bbApyTtBqwCvA84ExglaWTqrYwHFrQgNjMza0DTeyoR\ncUJEjI+ITmA/4I8RcQBwC7BPqjYVuKbZsZmZWWMG0+9UjgeOlTSX7BjLuS2Ox8zMBqilN+mKiFuB\nW9P048DWrYzHzMwaM5h6KmZm1uacVMzMLDdOKmZmlhsnFTMzy42TipmZ5cZJxczMcuOkYmZmuWnp\n71SsMZ3Tf1e2fN6puzc5EjOzjHsqZmaWG/dU2kClHomZ2WDjnoqZmeXGScXMzHLjpGJmZrlxUjEz\ns9w4qZiZWW6cVMzMLDc+pXgIGugpyP6xpJnlxT0VMzPLjZOKmZnlxknFzMxy0/RjKpImABcC6wEB\nzIiIMyWNBi4FOoF5wL4R8UKz47N3VTs24+MwZlZOK3oqy4CvRsRmwBTgSEmbAdOBmyNiEnBzem1m\nZm2k6UklIhZGxD1p+mVgDjAO2BO4IFW7ANir2bGZmVljWnpMRVInsBVwF7BeRCxMsxaRDY+VW2aa\npG5J3T09PU2J08zMatOypCJpDeBK4JiIeKl0XkQE2fGW94iIGRHRFRFdHR0dTYjUzMxq1ZKkImlF\nsoRyUURclYqflTQ2zR8LLG5FbGZmVr+mJxVJAs4F5kTE6SWzrgWmpumpwDXNjs3MzBrTisu0bAsc\nCDwgaVYqOxE4FbhM0mHAk8C+LYjNzMwa0PSkEhG3A6owe4dmxmJmZvnyL+rNzCw3TipmZpYbJxUz\nM8uN76diA77/iplZJe6pmJlZbpxUzMwsN04qZmaWGycVMzPLjQ/UW10qHdyvdPOugdY3s/bknoqZ\nmeXGScXMzHLj4S/LlX/zYja8OalYS7XqWIuP8ZgVw8NfZmaWGycVMzPLjYe/bFDyKctm7clJxYa0\nZpw44IRm9i4Pf5mZWW7cUzGrgU+VNquNk4pZCQ+XmTXGw19mZpabQddTkbQLcCYwAjgnIk5tcUg2\niAzlYaiBts09GxuMBlVSkTQC+DHwz8B84G+Sro2Ih1obmdnAtUsCzCuZVVtPXqeCD7b6ldTztxjo\nugbrl4rBNvy1NTA3Ih6PiDeAS4A9WxyTmZnVSBHR6hjeIWkfYJeIODy9PhD4aEQcVVJnGjAtvdwE\neKTOzY0BljQQ7mA01Nrk9gx+Q61NQ609UL5NH4iIjiI2NqiGv2oRETOAGY2uR1J3RHTlENKgMdTa\n5PYMfkOtTUOtPdD8Ng224a8FwISS1+NTmZmZtYHBllT+BkyStIGklYD9gGtbHJOZmdVoUA1/RcQy\nSUcB15OdUnxeRDxY0OYaHkIbhIZam9yewW+otWmotQea3KZBdaDezMza22Ab/jIzszbmpGJmZrkZ\nlklF0i6SHpE0V9L0VscDIGmepAckzZLUncpGS7pR0qPpee1ULklnpfjvlzS5ZD1TU/1HJU0tKf9w\nWv/ctKyqbaOO+M+TtFjS7JKylsVfbRsNtulkSQvSfpolabeSeSek7T0iaeeS8rLvt3RCyl2p/NJ0\ncgqSVk6v56b5nf1to8b2TJB0i6SHJD0o6cv1/g0Hw36q0p623EeSVpF0t6T7Unu+nXcMebazoogY\nVg+yEwAeAzYEVgLuAzYbBHHNA8b0Kft3YHqang6clqZ3A/4ACJgC3JXKRwOPp+e10/Taad7dqa7S\nsrtW20Yd8X8cmAzMHgzxV9pGDm06GfhambqbpffSysAG6T02otr7DbgM2C9N/wz4lzT9JeBnaXo/\n4NJq2xhAe8YCk9P0msDf0zrbcj9VaU9b7qP0N1gjTa8I3JX+JrnEkGc7q7ajng+Qdn4A2wDXl7w+\nAThhEMQ1j/cmlUeAsWl6LPBImv45sH/fesD+wM9Lyn+eysYCD5eUv1Ov0jbqbEMny38Atyz+StvI\noU0nU/4Da7n3EdkZjNtUer+RfYAsAUb2fV/2LpumR6Z6qrSNBvbXNWTX2Wv7/dSnPW2/j4DVgHuA\nj+YVQ57trBb7cBz+Ggc8XfJ6fiprtQBukDRT2aVoANaLiIVpehGwXpqu1IZq5fPLlFfbRh5aGX+R\n+/moNFRznt4dLhxom9YBlkbEsjLxvbNMmv9iqp9bm9IwxlZk34bbfj/1aQ+06T6SNELSLGAxcCNZ\nzyKvGPJsZ0XDMakMVttFxGRgV+BISR8vnRnZV4VCz/8uchvtHn+JnwIbAVsCC4EfFLy93ElaA7gS\nOCYiXiqd1477qUx72nYfRcRbEbEl2dVEtgY2bXFIAzYck8qgvBRMRCxIz4uBq8neUM9KGguQnhen\n6pXaUK18fJlyqmwjD62Mv5D9HBHPpn/8t4GzyfZTPW16DhglaWSf8uXWleavleo33CZJK5J9AF8U\nEVel4rbdT+Xa0+77KLVhKXAL2VBUXjHk2c6KhmNSGXSXgpG0uqQ1e6eBnYDZKa7eM2umko0Zk8oP\nSmfOTAFeTEML1wM7SVo7dfl3IhsbXQi8JGmKJAEH9VlXuW3koZXxV9pGQ3o/GJO9yfZT7/b2S2fL\nbABMIjtoXfb9lr6t3wLsUyH23jbtA/wx1a+0jVpjF3AuMCciTi+Z1Zb7qVJ72nUfSeqQNCpNr0p2\nfGhOjjHk2c7K6jkg1u4PsjNO/k42Xvn1QRDPhmRnYtwHPNgbE9nY5c3Ao8BNwOhULrKbmT0GPAB0\nlazrUGBuehxSUt5F9s/1GPAj3r2aQtlt1NGGi8mGGt4kG5M9rJXxV9tGg236ZVrf/ekfbmxJ/a+n\n7T1COuup2vst7fe7U1svB1ZO5auk13PT/A3720aN7dmObNjpfmBWeuzWrvupSnvach8BWwD3prhn\nA9/KO4Y821np4cu0mJlZbobj8JeZmRXEScXMzHLjpGJmZrlxUjEzs9w4qZiZWW6cVGzYkfSWsqvX\nzpZ0uaTVcl7/wZJ+1E+d7SX9j5LXR0g6KM84zFrBScWGo9ciYsuI2Bx4AziiBTFsD7yTVCLiZxFx\nYQviMMuVk4oNd7cBGwNIOjb1XmZLOiaVdUp6WNJFkuZIuqK3Z6PsHjhj0nSXpFv7rlzSJ5Xdh+Je\nSTdJWk/ZxQ+PAL6SekwfU3YPkK+lZbaUdKeyCyJerXfvPXKrpNOU3XPj75I+Vvyfx2xgnFRs2ErX\nMtoVeEDSh4FDyC41PgX4oqStUtVNgJ9ExAeBl8juMVGr24EpEbEVcAlwXETMI7tnxRmpx3Rbn2Uu\nBI6PiC3Ifhl+Usm8kRGxNXBMn3KzQcFJxYajVZVdXrwbeIrs+lHbAVdHxCsR8Q/gKqC3J/B0RPwl\nTf8q1a3VeOB6SQ8A/wr8U7XKktYCRkXEn1LRBWQ3C+vVexHImWT3ejEbVEb2X8VsyHktssuLvyO7\nNmFFfa9l1Pt6Ge9+MVulwrI/BE6PiGslbU92A6lGvJ6e38L/vzYIuadilrkN2EvSaulK0XunMoCJ\nkrZJ058nG9KC7G6dH07Tn6mw3rV49zLiU0vKXya7Be5yIuJF4IWS4yUHAn/qW89ssHJSMQMi4h7g\nfLIrsd4FnBMR96bZj5DdOG0O2T3Zf5rKvw2cKambrOdQzsnA5ZJmkt2Ktdd1wN69B+r7LDMV+L6k\n+8luNPVvjbTNrJl8lWKzKtKZWr9Npx+bWT/cUzEzs9y4p2JmZrlxT8XMzHLjpGJmZrlxUjEzs9w4\nqZiZWW6cVMzMLDf/H7nT5MpRysfwAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7feacf6d0400>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(prefix_3['2010 Census Population'], bins=50)\n",
    "plt.title(\"Population per 3-digit ZIP prefix\")\n",
    "plt.xlabel(\"Population\")\n",
    "plt.ylabel(\"# of Prefixes\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style>\n",
       "    .dataframe thead tr:only-child th {\n",
       "        text-align: right;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: left;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Zip Prefix</th>\n",
       "      <th>2010 Census Population</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>31</th>\n",
       "      <td>33</td>\n",
       "      <td>9285786</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>90</th>\n",
       "      <td>92</td>\n",
       "      <td>9033701</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>58</th>\n",
       "      <td>60</td>\n",
       "      <td>8783167</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>11</td>\n",
       "      <td>7579159</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>30</td>\n",
       "      <td>7254703</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>75</th>\n",
       "      <td>77</td>\n",
       "      <td>7012758</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>93</th>\n",
       "      <td>95</td>\n",
       "      <td>6786586</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>46</th>\n",
       "      <td>48</td>\n",
       "      <td>6634334</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>76</th>\n",
       "      <td>78</td>\n",
       "      <td>6358782</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30</th>\n",
       "      <td>32</td>\n",
       "      <td>6263654</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   Zip Prefix  2010 Census Population\n",
       "31         33                 9285786\n",
       "90         92                 9033701\n",
       "58         60                 8783167\n",
       "9          11                 7579159\n",
       "28         30                 7254703\n",
       "75         77                 7012758\n",
       "93         95                 6786586\n",
       "46         48                 6634334\n",
       "76         78                 6358782\n",
       "30         32                 6263654"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "prefix_2 = sum_by_zip_prefix(zippop, 2)\n",
    "prefix_2.nlargest(10, '2010 Census Population')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style>\n",
       "    .dataframe thead tr:only-child th {\n",
       "        text-align: right;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: left;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Zip Prefix</th>\n",
       "      <th>2010 Census Population</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>67</th>\n",
       "      <td>69</td>\n",
       "      <td>184566</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>49</th>\n",
       "      <td>51</td>\n",
       "      <td>486400</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>80</th>\n",
       "      <td>82</td>\n",
       "      <td>520260</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>86</th>\n",
       "      <td>88</td>\n",
       "      <td>620867</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>05</td>\n",
       "      <td>625741</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>84</th>\n",
       "      <td>86</td>\n",
       "      <td>644295</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>56</th>\n",
       "      <td>58</td>\n",
       "      <td>671781</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>26</td>\n",
       "      <td>795935</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55</th>\n",
       "      <td>57</td>\n",
       "      <td>813464</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>79</th>\n",
       "      <td>81</td>\n",
       "      <td>832537</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   Zip Prefix  2010 Census Population\n",
       "67         69                  184566\n",
       "49         51                  486400\n",
       "80         82                  520260\n",
       "86         88                  620867\n",
       "4          05                  625741\n",
       "84         86                  644295\n",
       "56         58                  671781\n",
       "24         26                  795935\n",
       "55         57                  813464\n",
       "79         81                  832537"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "prefix_2.nsmallest(10, '2010 Census Population')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.text.Text at 0x7feacf5d1f60>"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEWCAYAAABrDZDcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGmtJREFUeJzt3Xm4HFWZx/HvDyKyRdarEsL1AjIIgyhwFRxAGYIDsrqg\nwrBFGDI4AjIqiIwKMzoIw4jjzkR2QVAQEEZn2ATZdwJCIoIQIBAMAcISUAi880edkE57u2/dpaty\n+/w+z9PP7a4+VfVW3e5+65yqOkcRgZmZ5WupugMwM7N6ORGYmWXOicDMLHNOBGZmmXMiMDPLnBOB\nmVnmnAhs2CRtI2nWCOY/SdJXRjOmsUDSZEnXNbx+QdI6JectXXaskfR1SXMlPSGpN23r0nXHlQMn\ngi4haaakl9KX54+STpe0Yt1xLdT84wcQEQdFxNfqimm4JB0u6R5Jz0t6SNLhI1leRKwYEQ8OtWz6\nH3+9TZxHpc9D42O+pJC0byozU9J26flkSa+mcs9JmiZp55FsW1mSeoHPAxtGxFsj4pG0ra9Wsf7c\nORF0l10iYkVgU6Af+HLN8Yx5ksYNNBnYF1gF2AE4WNIelQZWQkQcm35MX38A3wKmAz9vMduNqdzK\nwCnAzyStMpI4VBjst6YXeCoi5oxkXTY8TgRdKCIeA/4X2AhA0gRJF0t6WtIDkg5cWFbSMZLOl/TT\ndIR7h6R3Nbwfkt7e8LrlUaikIyX9IS1nuqSPpOkbACcB70tHm/MGWpakA1N8T6d4JzTFcZCk+yXN\nk/R9SWoRx2DbNEHSzyU9mY7oDx1g3rMkPQdMHmD//kdE3BERCyLiPuAXwJat/h+SVkvb85ykW4B1\nm95/fR+nspeksrem5pLrmstKmgLsBRyR9uklrdbfMO+OwKHA7hExv13ZiHgNOBVYrjnetKzJkq6X\n9D1Jz0r6naRJDe9fLenfJV0PvAisI2klSadImi3psbRtS6cayeXAhLQtp0vqS9s6TtKqkmZJ2iUt\ne8X0Odl3sG22cpwIupCktYAdgTvTpHOBWcAEYHfgWEnbNsyyG3AesCrwE+AiSW8Yxqr/AGwNrAT8\nK3CWpDUiYgZwEOloMyJWHiDmbYFvAJ8A1gAeTnE32hl4D7BxKrd9m1gG3KZ0ZHoJcBewJjAJOEzS\n9k3znk9xVHx2uw1OyWhr4N42xb4P/Clt1/7p0a7sfOCtwH7p8RciYmqK7T/SPt1lkDj7gB8DB6b/\nR1upJvQPwAvA/S2KbU7xP18dOBq4QNKqDe/vA0wBxlP8P08HFgBvBzYB/g74h4i4AvgQ8HjalslN\n2/o0xT77kaQ3U9RqpkXEmYNth5XjRNBdLkpH29cBv6H4wV+L4mj1ixHxp4iYBpxM0bSx0O0RcX5E\nvAKcCCwLbDHUlUfEeRHxeES8FhE/pfgBeW/J2fcCTk1H2n8GvkRRg+hrKHNcRMyLiEeAq4B3t1le\nq216D9ATEf8WES+n9vYfAY1NOzdGxEVpO14aJO5jKL5Hpw30poqTnR8DvhoR8yPiHuCMQcoeHREv\nRsT0VmWHQtIbKRLb2en/0s4W6TP0BLAn8JGIeLZF2TnAf0XEK2m59wE7Nbx/ekTcGxELKBLyjsBh\naT/MofhBL9WkFhGXUST2K9Ny/rHMfFbOQO2fNnZ9OB1dvS41rzwdEc83TH6Y4hzCQo8ufBIRr6m4\nEmgCQ5Sq6p8D+tKkFSmOFsuYANzREMcLkp6iOGqfmSY/0VD+xbT8VlptU1A0QcxrKLs0cO1A87Yj\n6WCKhLp1Sl5IOgo4KhU5iyJRjGta5sMtFtkzQNlSsQzi28ArFCdjB3NTRGxVcrmPxeK9Vj7M4p+b\nxtjfBrwBmN3QorcUQ9u+qcDBwLER8dQQ5rNBOBF0v8eBVSWNb0gGvcBjDWXWWvgkNZ1MTPNB8YO7\nfEPZt1I0My1G0tsojqwnURxRvyppGsWJVSh+gAeL820Ny1sBWK0pzqFotU0LgIciYr028w7aJa+k\n/YEjgfdHxOv7IyKOBY5tKLd0WudawO/S5N4Wi30ylZ0I/L55O4YZ5z4UtYxNU+1oNK0pSQ3JoBe4\nuEV8jwJ/BlZPNYQhSftxKnAm8E+STouIB4YZtzVx01CXi4hHgRuAb0haVtLGwAEUR6sLbSbpo6ld\n+DCKL+xN6b1pwN+nk3o7AB9osaoVKL74TwJI+hTpZHXyR2CipGVazH8O8ClJ705NGccCN0fEzKFt\n8aDbdAvwvKQvSloubddGkt5TdsGS9krxfXCwyz7T5Y8XAMdIWl7ShrRu928u+w4Wb8Jr9keg5T0F\nkjYCfgDslT4Ho+3NwKHp3MvHgQ2AXw1UMCJmA5cB35T0JklLSVpXUqvPU7OjKD5f+wMnAGfK9xiM\nGieCPOxJ0VzzOHAhRRt0YxPSL4BPAs9QnOD7aMPR42eBXYB5FO34Fw20gtSe/U3gRoofqHcC1zcU\n+TXFCdUnJM0dYP4rgK9QXNY4m+JKlZFckjngNqUf250pzi88BMylOGey0hCW/XWK2sqtWnR9/klt\nyh9M0Yz1BMUJ0wHPJzSUXSmV/TFFgvxzi7KnABuquIpqoP/L5ygS9AX6y/sJjhqg/FDdDKxHsQ//\nneJqpHZNNvsCy1BcvvoMxXmLNQZbiaTNKLZl3/T/O54iKRw5oujtdfLANHmTdAzw9ojYu+5YRks3\nbZOk44G3RsSAtYi6SJpMccVP2fMJtgRzjcBsCSLpHZI2VuG9FM14F9Ydl3U3nyw2W7KMp2gOmkDR\nxPZNimYus45x05CZWebcNGRmlrkx0TS0+uqrR19fX91hmJmNKbfffvvciOgZrNyYSAR9fX3cdttt\ndYdhZjamSGp1F/ti3DRkZpY5JwIzs8w5EZiZZc6JwMwsc04EZmaZcyIwM8tcxxKBpFMlzZF0T8O0\nE9LYpndLulDSXwxZaGZm1epkjeB0YIemaZcDG0XExhQDb3ypg+s3M7MSOpYIIuIa4OmmaZc1jE50\nE8VITGZmVqM67yzeH2g5kLakKcAUgN7eViP7da++I3857HlnHrfT4IXMzJJaThZL+heKsVnPblUm\nIqZGRH9E9Pf0DNpVhpmZDVPlNYI0stHOwKRwH9hmZrWrNBGkwc+PAD4QES9WuW4zMxtYJy8fPYdi\nIPP1Jc2SdADwPYoRmC6XNG2QAb/NzKwCHasRRMSeA0w+pVPrMzOz4fGdxWZmmXMiMDPLnBOBmVnm\nnAjMzDLnRGBmljknAjOzzDkRmJllzonAzCxzTgRmZplzIjAzy5wTgZlZ5uocmMY6xIPamNlQuEZg\nZpY5JwIzs8w5EZiZZc6JwMwsc04EZmaZcyIwM8ucE4GZWeacCMzMMudEYGaWOScCM7PMORGYmWXO\nicDMLHNOBGZmmetYIpB0qqQ5ku5pmLaqpMsl3Z/+rtKp9ZuZWTmdrBGcDuzQNO1I4MqIWA+4Mr02\nM7MadSwRRMQ1wNNNk3cDzkjPzwA+3Kn1m5lZOVUPTPOWiJidnj8BvKVVQUlTgCkAvb29FYQ2+kYy\nQIyZWVVqO1kcEQFEm/enRkR/RPT39PRUGJmZWV6qTgR/lLQGQPo7p+L1m5lZk6oTwcXAfun5fsAv\nKl6/mZk16eTlo+cANwLrS5ol6QDgOOCDku4HtkuvzcysRh07WRwRe7Z4a1Kn1mlmZkPnO4vNzDLn\nRGBmljknAjOzzDkRmJllzonAzCxzTgRmZplzIjAzy5wTgZlZ5pwIzMwy50RgZpY5JwIzs8xVPTCN\nLeFGMpjOzON2GsVIzKwqrhGYmWXOicDMLHNOBGZmmXMiMDPLnBOBmVnmnAjMzDLnRGBmljknAjOz\nzDkRmJllzonAzCxzTgRmZplzIjAzy5wTgZlZ5pwIzMwyV0sikPTPku6VdI+kcyQtW0ccZmZWQyKQ\ntCZwKNAfERsBSwN7VB2HmZkVBk0Ekj4uaXx6/mVJF0jadITrHQcsJ2kcsDzw+AiXZ2Zmw1RmhLKv\nRMR5krYCtgNOAH4IbD6cFUbEY5L+E3gEeAm4LCIuay4naQowBaC3t3c4qxqxkYzWZWY2VpRpGno1\n/d0JmBoRvwSWGe4KJa0C7AasDUwAVpC0d3O5iJgaEf0R0d/T0zPc1ZmZ2SDKJILHJP038EngV5Le\nWHK+VrYDHoqIJyPiFeAC4G9GsDwzMxuBMj/onwAuBbaPiHnAqsDhI1jnI8AWkpaXJGASMGMEyzMz\nsxEYNBFExIvAHGCrNGkBcP9wVxgRNwPnA3cAv00xTB3u8szMbGQGPVks6WigH1gfOA14A3AWsOVw\nVxoRRwNHD3d+MzMbPWWahj4C7ArMB4iIx4HxnQzKzMyqUyYRvBwRAQSApBU6G5KZmVWpTCL4Wbpq\naGVJBwJXACd3NiwzM6vKoOcIIuI/JX0QeI7iPMFXI+LyjkdmZmaVKHOyeMP0w395w7RtIuLqTgZm\nZmbVKNs0dIQKy0n6LvCNTgdmZmbVKJMINgd6gRuAWyk6iBv2paNmZrZkKZMIXqHoHG45YFmK7iFe\n62hUZmZWmTKJ4FaKRPAeYGtgT0nndTQqMzOrTJluqA+IiNvS89nAbpL26WBMZmZWoZaJQNKbIuI5\n4EFJqza97Y76zcy6RLsawU+AnYHbKe4qVsN7AazTwbjMKjOSAYhmHrfTKEZiVo92ieC49HeDiPhT\nFcGYmVn12p0s/nb6e0MVgZiZWT3a1QhekTQVmCjpO81vRsShnQvLzMyq0i4R7EwxrOT2FOcJzMys\nC7VMBBExFzhX0oyIuKvCmMzMrEJlbih7SdKVku4BkLSxpC93OC4zM6tImUTwI+BLFF1NEBF3A3t0\nMigzM6tOmUSwfETc0jRtQSeCMTOz6pVJBHMlrcuioSp3p+hqwszMukCZvoY+A0wF3iHpMeAhYK+O\nRmVmZpVpmwgkLQX0R8R2adD6pSLi+WpCMzOzKrRtGkrjDhyRns93EjAz6z5lzhFcIekLktaStOrC\nR8cjMzOzSpQ5R/DJ9PczDdNG1PuopJWBk4GN0rL2j4gbh7s8MzMbvkETQUSs3YH1fhv4v4jYXdIy\nwPIdWIeZmZXQsmlI0nqSfiHpHknnSFpzNFYoaSXg/cApABHxckTMG41lm5nZ0LWrEZwKnAlcA+wK\nfBf46Cisc23gSeA0Se+i6NDusxExv7GQpCnAFIDe3t5RWK112kgGeAEP8mJWl3Yni8dHxI8i4r6I\nOAHoG6V1jgM2BX4YEZsA84EjmwtFxNSI6I+I/p6enlFatZmZNWtXI1hW0iYsGqJyucbXEXHHMNc5\nC5gVETen1+czQCIwM7NqtEsEs4ETG14/0fA6gG2Hs8KIeELSo5LWj4j7gEnA9OEsy8zMRq7deAR/\n28H1HgKcna4YehD4VAfXZWZmbZS5j2DURcQ0oL+OdZuZ2eLK3FlsZmZdrN19BFumv2+sLhwzM6ta\nuxrBd9Jfd/1gZtbF2p0jeEXSVGBNSd9pfjMiDu1cWGZmVpV2iWBnYDtge4q7f83MrAu1u3x0LnCu\npBkRcVeFMZmZWYXKXDX0lKQLJc1Jj59LmtjxyMzMrBJlEsFpwMXAhPS4JE0zM7MuUCYRvDkiTouI\nBelxOuBe4MzMukSZRDBX0t6Slk6PvYGnOh2YmZlVo0wi2B/4BEWnc7OB3XHfQGZmXaPMUJUPUwxM\nY2ZmXaiWTueqNNJRs6w6/l+Z1cOdzpmZZc6JwMwsc4MmAklfbnjunkjNzLpMu26ovyjpfRRXCS3k\nnkjNzLpMu5PFvwM+Dqwj6dr0erWGsYbNzKwLtGsamgccBTwAbAN8O00/UtINHY7LzMwq0q5GsD3w\nVWBd4ETgbmB+RPhmMjOzLtKyRhARR0XEJGAm8GNgaaBH0nWSLqkoPjMz67AyN5RdGhG3AbdJ+nRE\nbCVp9U4HZmZm1Rj08tGIOKLh5eQ0bW6nAjIzs2oN6YYyj1RmZtZ9fGexmVnmnAjMzDJXWyJIg9zc\nKel/6orBzMzqrRF8FphR4/rNzIyaEoGkicBOwMl1rN/MzBapa2Ca/wKOAMa3KiBpCjAFoLe3t6Kw\nzIZmpIPpzDxup1GKxGz4Kq8RSNoZmBMRt7crFxFTI6I/Ivp7enoqis7MLD91NA1tCewqaSZwLrCt\npLNqiMPMzKghEUTElyJiYkT0AXsAv46IvauOw8zMCr6PwMwsc3WdLAYgIq4Grq4zBjOz3LlGYGaW\nOScCM7PMORGYmWXOicDMLHNOBGZmmXMiMDPLnBOBmVnmnAjMzDLnRGBmljknAjOzzDkRmJllrta+\nhsxyN9KBbYbLA+JYI9cIzMwy50RgZpY5JwIzs8w5EZiZZc6JwMwsc04EZmaZcyIwM8ucE4GZWeac\nCMzMMudEYGaWOScCM7PMORGYmWXOicDMLHOVJwJJa0m6StJ0SfdK+mzVMZiZ2SJ1dEO9APh8RNwh\naTxwu6TLI2J6DbGYmWWv8hpBRMyOiDvS8+eBGcCaVcdhZmaFWgemkdQHbALcPMB7U4ApAL29vZXG\nZdbt6hoQZ6TG4oA6I93XVWxzbSeLJa0I/Bw4LCKea34/IqZGRH9E9Pf09FQfoJlZJmpJBJLeQJEE\nzo6IC+qIwczMCnVcNSTgFGBGRJxY9frNzGxxddQItgT2AbaVNC09dqwhDjMzo4aTxRFxHaCq12tm\nZgPzncVmZplzIjAzy5wTgZlZ5pwIzMwy50RgZpY5JwIzs8w5EZiZZc6JwMwsc04EZmaZcyIwM8uc\nE4GZWeacCMzMMlfrCGVmZkMxktG+RjLS11gd0a0s1wjMzDLnRGBmljknAjOzzDkRmJllzonAzCxz\nTgRmZplzIjAzy5wTgZlZ5pwIzMwy50RgZpY5JwIzs8w5EZiZZc6JwMwsc7UkAkk7SLpP0gOSjqwj\nBjMzK1SeCCQtDXwf+BCwIbCnpA2rjsPMzAp11AjeCzwQEQ9GxMvAucBuNcRhZmbUMzDNmsCjDa9n\nAZs3F5I0BZiSXr4g6b5Blrs6MHdUIhz7vC8W8b5YJOt9oeMXezlm9kVT3EP1tjKFltgRyiJiKjC1\nbHlJt0VEfwdDGjO8LxbxvljE+2IR74vF1dE09BiwVsPriWmamZnVoI5EcCuwnqS1JS0D7AFcXEMc\nZmZGDU1DEbFA0sHApcDSwKkRce8oLLp0M1IGvC8W8b5YxPtiEe+LBoqIumMwM7Ma+c5iM7PMORGY\nmWVuzCWCwbqnkPRGST9N798sqa/6KKtRYl98TtJ0SXdLulJSqWuKx6Ky3ZZI+pikkNSVlw6W2Q+S\nPpE+F/dK+knVMValxPejV9JVku5M35Ed64hziRARY+ZBcXL5D8A6wDLAXcCGTWX+CTgpPd8D+Gnd\ncde4L/4WWD49/3TO+yKVGw9cA9wE9Ncdd02fifWAO4FV0us31x13jftiKvDp9HxDYGbdcdf1GGs1\ngjLdU+wGnJGenw9MkqQKY6zKoPsiIq6KiBfTy5so7tnoRmW7LfkacDzwpyqDq1CZ/XAg8P2IeAYg\nIuZUHGNVyuyLAN6Unq8EPF5hfEuUsZYIBuqeYs1WZSJiAfAssFol0VWrzL5odADwvx2NqD6D7gtJ\nmwJrRcQvqwysYmU+E38F/JWk6yXdJGmHyqKrVpl9cQywt6RZwK+AQ6oJbcmzxHYxYaNH0t5AP/CB\numOpg6SlgBOByTWHsiQYR9E8tA1FDfEaSe+MiHm1RlWPPYHTI+Kbkt4H/FjSRhHxWt2BVW2s1QjK\ndE/xehlJ4yiqfE9VEl21SnXVIWk74F+AXSPizxXFVrXB9sV4YCPgakkzgS2Ai7vwhHGZz8Qs4OKI\neCUiHgJ+T5EYuk2ZfXEA8DOAiLgRWJaiM7rsjLVEUKZ7iouB/dLz3YFfRzob1GUG3ReSNgH+myIJ\ndGtbMAyyLyLi2YhYPSL6IqKP4nzJrhFxWz3hdkyZ78dFFLUBJK1O0VT0YJVBVqTMvngEmAQgaQOK\nRPBkpVEuIcZUIkht/gu7p5gB/Cwi7pX0b5J2TcVOAVaT9ADwOaArR0AruS9OAFYEzpM0TVJX9ulU\ncl90vZL74VLgKUnTgauAwyOi62rMJffF54EDJd0FnANM7tKDxkG5iwkzs8yNqRqBmZmNPicCM7PM\nORGYmWXOicDMLHNOBGZmSxhJp0qaI+meEmW/la4KnCbp95KGfHOgrxqyriXpVeC3FHfTzgD2a+h7\naTSWP5mi87qD25TZBng5Im5Irw8CXoyIM0crDus+kt4PvACcGREbDWG+Q4BNImL/oazPNQLrZi9F\nxLvTF+ll4KAaYtgG+JuFLyLiJCcBG0xEXAM83ThN0rqS/k/S7ZKulfSOAWbdk+KeiCFxIrBcXAu8\nHV4fp+Ge9DgsTeuT9DtJZ0uaIel8Scun92amu3CR1C/p6uaFS9oljX9xp6QrJL0ljYVxEPDPqdq+\ntaRjJH0hzfPu1PHb3ZIulLRKmn61pOMl3ZKq+lt3fvfYGDAVOCQiNgO+APyg8c003sjawK+HumAn\nAut6qc+pDwG/lbQZ8Clgc4o+hw5MXXEArA/8ICI2AJ6jGNuirOuALSJiE4ouj4+IiJnAScC3Us3k\n2qZ5zgS+GBEbUzRhHd3w3riIeC9wWNN0y5CkFSlqludJmkbRdcwaTcX2AM6PiFeHunwnAutmy6Uv\nzW0U/cqcAmwFXBgR8yPiBeACYOER96MRcX16flYqW9ZE4FJJvwUOB/66XWFJKwErR8Rv0qQzgPc3\nFLkg/b0d6BtCHNadlgLmpQOKhY8NmsrswTCahRYu3KxbvdTwpTkkDVDSTvOVEwtfL2DRd2XZFvN+\nF/heRLwT+Mc25cpa2FPsq7i7+OxFxHPAQ5I+DqDCuxa+n84XrALcOJzlOxFYbq4FPixpeUkrAB9J\n0wB6U7/0AH9P0dwDMBPYLD3/WIvlrsSibo73a5j+PEU32IuJiGeBZxra//cBftNczvIk6RyKH/X1\nJc2SdACwF3BA6iTvXhYfcW0P4NzhdprnIw3LSkTcIel04JY06eSIuDOd2L0P+IykU4HpwA9TmX8F\nTpH0NeDqFos+hqL99hmKk3Vrp+mXAOdL2o2/HAFrP+CkdFL6QYpzF2ZExJ4t3hpwRLmIOGYk6/N9\nBGYUVw0B/zOUa7bNuoWbhszMMucagZlZ5lwjMDPLnBOBmVnmnAjMzDLnRGBmljknAjOzzP0/zqUB\nf18naXoAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7feacf5e1e10>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(prefix_2['2010 Census Population'], bins=20)\n",
    "plt.title(\"Population per 2-digit ZIP prefix\")\n",
    "plt.xlabel(\"Population\")\n",
    "plt.ylabel(\"# of Prefixes\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style>\n",
       "    .dataframe thead tr:only-child th {\n",
       "        text-align: right;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: left;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Zip Prefix</th>\n",
       "      <th>2010 Census Population</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>23236989</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>32978141</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2</td>\n",
       "      <td>30390138</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>3</td>\n",
       "      <td>42575790</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>4</td>\n",
       "      <td>32249466</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>5</td>\n",
       "      <td>16513346</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>6</td>\n",
       "      <td>23499661</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>7</td>\n",
       "      <td>36345281</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>8</td>\n",
       "      <td>21075552</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>9</td>\n",
       "      <td>49875567</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "  Zip Prefix  2010 Census Population\n",
       "0          0                23236989\n",
       "1          1                32978141\n",
       "2          2                30390138\n",
       "3          3                42575790\n",
       "4          4                32249466\n",
       "5          5                16513346\n",
       "6          6                23499661\n",
       "7          7                36345281\n",
       "8          8                21075552\n",
       "9          9                49875567"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sum_by_zip_prefix(zippop, 1)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Variable length chunks\n",
    "\n",
    "Now that we've explored the data, we're ready to experiment with variable length chunks.  First, we want to build a table that will allow us to quickly look up the population of any prefix from length 0 (the entire country) to length 5 (a single zip)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[2906700, 10538]"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "prefix_populations = [ sum_by_zip_prefix(zippop, prefix).set_index('Zip Prefix') for prefix in range(0,6) ]\n",
    "def prefix_population(prefix):\n",
    "    return prefix_populations[len(prefix)]['2010 Census Population'].get(prefix, 0)\n",
    "[prefix_population('770'), prefix_population('03755')]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['9020',\n",
       " '9021',\n",
       " '9022',\n",
       " '9023',\n",
       " '9024',\n",
       " '9025',\n",
       " '9026',\n",
       " '9027',\n",
       " '9028',\n",
       " '9029']"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def sub_prefixes(prefix):\n",
    "    return [ prefix + str(n) for n in range(0,10) ]\n",
    "sub_prefixes('902')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style>\n",
       "    .dataframe thead tr:only-child th {\n",
       "        text-align: right;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: left;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>2010 Census Population</th>\n",
       "      <th>Prefixes</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Chunk</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>010_0</th>\n",
       "      <td>77922</td>\n",
       "      <td>[0100]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>010_1</th>\n",
       "      <td>96062</td>\n",
       "      <td>[0101, 0102]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>010_2</th>\n",
       "      <td>76753</td>\n",
       "      <td>[0103, 0104]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>010_3</th>\n",
       "      <td>73002</td>\n",
       "      <td>[0105, 0106]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>010_4</th>\n",
       "      <td>37535</td>\n",
       "      <td>[0107]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>010_5</th>\n",
       "      <td>88283</td>\n",
       "      <td>[0108]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>010_6</th>\n",
       "      <td>19818</td>\n",
       "      <td>[0109]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0110_0</th>\n",
       "      <td>92014</td>\n",
       "      <td>[01100, 01101, 01102, 01103, 01104, 01105, 011...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0110_1</th>\n",
       "      <td>30250</td>\n",
       "      <td>[01109]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>011_0</th>\n",
       "      <td>46571</td>\n",
       "      <td>[0111, 0112, 0113, 0114, 0115, 0116, 0117, 011...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>012_0</th>\n",
       "      <td>89837</td>\n",
       "      <td>[0120, 0121, 0122, 0123]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>012_1</th>\n",
       "      <td>40852</td>\n",
       "      <td>[0124, 0125, 0126, 0127, 0128, 0129]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>014_0</th>\n",
       "      <td>82626</td>\n",
       "      <td>[0140, 0141, 0142, 0143, 0144]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>014_1</th>\n",
       "      <td>60658</td>\n",
       "      <td>[0145]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>014_2</th>\n",
       "      <td>68528</td>\n",
       "      <td>[0146, 0147, 0148, 0149]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>015_0</th>\n",
       "      <td>82551</td>\n",
       "      <td>[0150, 0151]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>015_1</th>\n",
       "      <td>85393</td>\n",
       "      <td>[0152, 0153]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>015_2</th>\n",
       "      <td>76897</td>\n",
       "      <td>[0154, 0155]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>015_3</th>\n",
       "      <td>79853</td>\n",
       "      <td>[0156, 0157]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>015_4</th>\n",
       "      <td>49025</td>\n",
       "      <td>[0158, 0159]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0160_0</th>\n",
       "      <td>77412</td>\n",
       "      <td>[01600, 01601, 01602, 01603, 01604]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0160_1</th>\n",
       "      <td>79958</td>\n",
       "      <td>[01605, 01606, 01607, 01608, 01609]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>016_0</th>\n",
       "      <td>31032</td>\n",
       "      <td>[0161, 0162, 0163, 0164, 0165, 0166, 0167, 016...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>017_0</th>\n",
       "      <td>73888</td>\n",
       "      <td>[0170, 0171]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>017_1</th>\n",
       "      <td>52571</td>\n",
       "      <td>[0172, 0173]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>017_2</th>\n",
       "      <td>81318</td>\n",
       "      <td>[0174]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>017_3</th>\n",
       "      <td>82505</td>\n",
       "      <td>[0175]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>017_4</th>\n",
       "      <td>88745</td>\n",
       "      <td>[0176, 0177, 0178, 0179]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0184_0</th>\n",
       "      <td>76383</td>\n",
       "      <td>[01840, 01841, 01842, 01843]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0184_1</th>\n",
       "      <td>75601</td>\n",
       "      <td>[01844, 01845, 01846, 01847, 01848, 01849]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>988_1</th>\n",
       "      <td>68451</td>\n",
       "      <td>[9882, 9883]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>988_2</th>\n",
       "      <td>58751</td>\n",
       "      <td>[9884, 9885, 9886, 9887, 9888, 9889]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9890_0</th>\n",
       "      <td>91008</td>\n",
       "      <td>[98900, 98901, 98902, 98903, 98904, 98905, 989...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9890_1</th>\n",
       "      <td>35240</td>\n",
       "      <td>[98908, 98909]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>989_0</th>\n",
       "      <td>77532</td>\n",
       "      <td>[9891, 9892, 9893]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>989_1</th>\n",
       "      <td>80875</td>\n",
       "      <td>[9894, 9895, 9896, 9897, 9898, 9899]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>990_0</th>\n",
       "      <td>82533</td>\n",
       "      <td>[9900, 9901]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>990_1</th>\n",
       "      <td>59475</td>\n",
       "      <td>[9902, 9903, 9904, 9905, 9906, 9907, 9908, 9909]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>991_0</th>\n",
       "      <td>67655</td>\n",
       "      <td>[9910, 9911, 9912, 9913, 9914, 9915]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>991_1</th>\n",
       "      <td>47583</td>\n",
       "      <td>[9916, 9917, 9918, 9919]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9920_0</th>\n",
       "      <td>61795</td>\n",
       "      <td>[99200, 99201, 99202, 99203, 99204]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9920_1</th>\n",
       "      <td>76939</td>\n",
       "      <td>[99205, 99206]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9920_2</th>\n",
       "      <td>80047</td>\n",
       "      <td>[99207, 99208, 99209]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>992_0</th>\n",
       "      <td>76513</td>\n",
       "      <td>[9921]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>992_1</th>\n",
       "      <td>47562</td>\n",
       "      <td>[9922, 9923, 9924, 9925, 9926, 9927, 9928, 9929]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>993_0</th>\n",
       "      <td>99875</td>\n",
       "      <td>[9930, 9931, 9932]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>993_1</th>\n",
       "      <td>91450</td>\n",
       "      <td>[9933]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>993_2</th>\n",
       "      <td>33446</td>\n",
       "      <td>[9934]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>993_3</th>\n",
       "      <td>79378</td>\n",
       "      <td>[9935]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>993_4</th>\n",
       "      <td>44835</td>\n",
       "      <td>[9936, 9937, 9938, 9939]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9950_0</th>\n",
       "      <td>97248</td>\n",
       "      <td>[99500, 99501, 99502, 99503, 99504]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9950_1</th>\n",
       "      <td>87630</td>\n",
       "      <td>[99505, 99506, 99507, 99508, 99509]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>995_0</th>\n",
       "      <td>97276</td>\n",
       "      <td>[9951, 9952, 9953, 9954, 9955, 9956]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>995_1</th>\n",
       "      <td>38506</td>\n",
       "      <td>[9957, 9958, 9959]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>996_0</th>\n",
       "      <td>78544</td>\n",
       "      <td>[9960, 9961, 9962, 9963, 9964]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>996_1</th>\n",
       "      <td>89277</td>\n",
       "      <td>[9965, 9966, 9967]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>996_2</th>\n",
       "      <td>16262</td>\n",
       "      <td>[9968, 9969]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>997_0</th>\n",
       "      <td>96393</td>\n",
       "      <td>[9970, 9971]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>997_1</th>\n",
       "      <td>37796</td>\n",
       "      <td>[9972, 9973, 9974, 9975, 9976, 9977, 9978, 9979]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>99_0</th>\n",
       "      <td>92694</td>\n",
       "      <td>[994, 998, 999]</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>4444 rows × 2 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "        2010 Census Population  \\\n",
       "Chunk                            \n",
       "010_0                    77922   \n",
       "010_1                    96062   \n",
       "010_2                    76753   \n",
       "010_3                    73002   \n",
       "010_4                    37535   \n",
       "010_5                    88283   \n",
       "010_6                    19818   \n",
       "0110_0                   92014   \n",
       "0110_1                   30250   \n",
       "011_0                    46571   \n",
       "012_0                    89837   \n",
       "012_1                    40852   \n",
       "014_0                    82626   \n",
       "014_1                    60658   \n",
       "014_2                    68528   \n",
       "015_0                    82551   \n",
       "015_1                    85393   \n",
       "015_2                    76897   \n",
       "015_3                    79853   \n",
       "015_4                    49025   \n",
       "0160_0                   77412   \n",
       "0160_1                   79958   \n",
       "016_0                    31032   \n",
       "017_0                    73888   \n",
       "017_1                    52571   \n",
       "017_2                    81318   \n",
       "017_3                    82505   \n",
       "017_4                    88745   \n",
       "0184_0                   76383   \n",
       "0184_1                   75601   \n",
       "...                        ...   \n",
       "988_1                    68451   \n",
       "988_2                    58751   \n",
       "9890_0                   91008   \n",
       "9890_1                   35240   \n",
       "989_0                    77532   \n",
       "989_1                    80875   \n",
       "990_0                    82533   \n",
       "990_1                    59475   \n",
       "991_0                    67655   \n",
       "991_1                    47583   \n",
       "9920_0                   61795   \n",
       "9920_1                   76939   \n",
       "9920_2                   80047   \n",
       "992_0                    76513   \n",
       "992_1                    47562   \n",
       "993_0                    99875   \n",
       "993_1                    91450   \n",
       "993_2                    33446   \n",
       "993_3                    79378   \n",
       "993_4                    44835   \n",
       "9950_0                   97248   \n",
       "9950_1                   87630   \n",
       "995_0                    97276   \n",
       "995_1                    38506   \n",
       "996_0                    78544   \n",
       "996_1                    89277   \n",
       "996_2                    16262   \n",
       "997_0                    96393   \n",
       "997_1                    37796   \n",
       "99_0                     92694   \n",
       "\n",
       "                                                 Prefixes  \n",
       "Chunk                                                      \n",
       "010_0                                              [0100]  \n",
       "010_1                                        [0101, 0102]  \n",
       "010_2                                        [0103, 0104]  \n",
       "010_3                                        [0105, 0106]  \n",
       "010_4                                              [0107]  \n",
       "010_5                                              [0108]  \n",
       "010_6                                              [0109]  \n",
       "0110_0  [01100, 01101, 01102, 01103, 01104, 01105, 011...  \n",
       "0110_1                                            [01109]  \n",
       "011_0   [0111, 0112, 0113, 0114, 0115, 0116, 0117, 011...  \n",
       "012_0                            [0120, 0121, 0122, 0123]  \n",
       "012_1                [0124, 0125, 0126, 0127, 0128, 0129]  \n",
       "014_0                      [0140, 0141, 0142, 0143, 0144]  \n",
       "014_1                                              [0145]  \n",
       "014_2                            [0146, 0147, 0148, 0149]  \n",
       "015_0                                        [0150, 0151]  \n",
       "015_1                                        [0152, 0153]  \n",
       "015_2                                        [0154, 0155]  \n",
       "015_3                                        [0156, 0157]  \n",
       "015_4                                        [0158, 0159]  \n",
       "0160_0                [01600, 01601, 01602, 01603, 01604]  \n",
       "0160_1                [01605, 01606, 01607, 01608, 01609]  \n",
       "016_0   [0161, 0162, 0163, 0164, 0165, 0166, 0167, 016...  \n",
       "017_0                                        [0170, 0171]  \n",
       "017_1                                        [0172, 0173]  \n",
       "017_2                                              [0174]  \n",
       "017_3                                              [0175]  \n",
       "017_4                            [0176, 0177, 0178, 0179]  \n",
       "0184_0                       [01840, 01841, 01842, 01843]  \n",
       "0184_1         [01844, 01845, 01846, 01847, 01848, 01849]  \n",
       "...                                                   ...  \n",
       "988_1                                        [9882, 9883]  \n",
       "988_2                [9884, 9885, 9886, 9887, 9888, 9889]  \n",
       "9890_0  [98900, 98901, 98902, 98903, 98904, 98905, 989...  \n",
       "9890_1                                     [98908, 98909]  \n",
       "989_0                                  [9891, 9892, 9893]  \n",
       "989_1                [9894, 9895, 9896, 9897, 9898, 9899]  \n",
       "990_0                                        [9900, 9901]  \n",
       "990_1    [9902, 9903, 9904, 9905, 9906, 9907, 9908, 9909]  \n",
       "991_0                [9910, 9911, 9912, 9913, 9914, 9915]  \n",
       "991_1                            [9916, 9917, 9918, 9919]  \n",
       "9920_0                [99200, 99201, 99202, 99203, 99204]  \n",
       "9920_1                                     [99205, 99206]  \n",
       "9920_2                              [99207, 99208, 99209]  \n",
       "992_0                                              [9921]  \n",
       "992_1    [9922, 9923, 9924, 9925, 9926, 9927, 9928, 9929]  \n",
       "993_0                                  [9930, 9931, 9932]  \n",
       "993_1                                              [9933]  \n",
       "993_2                                              [9934]  \n",
       "993_3                                              [9935]  \n",
       "993_4                            [9936, 9937, 9938, 9939]  \n",
       "9950_0                [99500, 99501, 99502, 99503, 99504]  \n",
       "9950_1                [99505, 99506, 99507, 99508, 99509]  \n",
       "995_0                [9951, 9952, 9953, 9954, 9955, 9956]  \n",
       "995_1                                  [9957, 9958, 9959]  \n",
       "996_0                      [9960, 9961, 9962, 9963, 9964]  \n",
       "996_1                                  [9965, 9966, 9967]  \n",
       "996_2                                        [9968, 9969]  \n",
       "997_0                                        [9970, 9971]  \n",
       "997_1    [9972, 9973, 9974, 9975, 9976, 9977, 9978, 9979]  \n",
       "99_0                                      [994, 998, 999]  \n",
       "\n",
       "[4444 rows x 2 columns]"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "class ZipChunk:\n",
    "    \"A group of zip codes roughly equal in size to all other such chunks.\"\n",
    " \n",
    "    def __init__(self, id):\n",
    "        self.id = id\n",
    "        self.prefixes = []\n",
    "\n",
    "    def add_prefix(self, prefix):\n",
    "        self.prefixes.append(prefix)\n",
    "    \n",
    "    def population(self):\n",
    "        return sum([ prefix_population(prefix) for prefix in self.prefixes ])\n",
    "    \n",
    "def build_chunks_recursive(population, chunks, prefix):\n",
    "    lower_limit = population / 2\n",
    "    upper_limit = population\n",
    "    if lower_limit < prefix_population(prefix) <= upper_limit or len(prefix) >= 5:\n",
    "        chunk = ZipChunk(prefix)\n",
    "        chunk.add_prefix(prefix)\n",
    "        chunks.append(chunk)\n",
    "    else:\n",
    "        leftovers = []\n",
    "        for child in sub_prefixes(prefix): \n",
    "            if prefix_population(child) >= upper_limit:\n",
    "                build_chunks_recursive(population, chunks, child)\n",
    "            else:\n",
    "                leftovers.append(child)\n",
    "        i = 0\n",
    "        chunk = ZipChunk(\"%s_%i\" % (prefix, i))\n",
    "        while len(leftovers) > 0:\n",
    "            child = leftovers.pop(0)\n",
    "            assert(prefix_population(child) <= upper_limit)\n",
    "            if chunk.population() + prefix_population(child) > upper_limit:\n",
    "                i += 1\n",
    "                chunks.append(chunk)\n",
    "                chunk = ZipChunk(\"%s_%i\" % (prefix, i))\n",
    "            chunk.add_prefix(child)\n",
    "        if chunk.population() > 0:\n",
    "            chunks.append(chunk)\n",
    "    \n",
    "def build_chunks(population):\n",
    "    chunks = []\n",
    "    build_chunks_recursive(population, chunks, '')\n",
    "    chunks = sorted(chunks, key=lambda chunk: chunk.id)\n",
    "    result = pd.DataFrame([ (chunk.id, chunk.population(), chunk.prefixes) for chunk in chunks ])\n",
    "    result.columns = ['Chunk', '2010 Census Population', 'Prefixes']\n",
    "    result.set_index('Chunk', inplace=True)\n",
    "    return result\n",
    "\n",
    "chunks = build_chunks(100000)\n",
    "chunks"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.text.Text at 0x7feacf516b00>"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEWCAYAAAB8LwAVAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGo1JREFUeJzt3X28ZmVd7/HPV1BUMAGZaHhqUMlX5FHUiSDtHI6WKGbg\nQxwo5UGU7OgrKU8KZomWhafCoycPOIoJqeADgkSWCZpJKTaDJAiSowwy4wCDiuJDCPg7f6xry81m\n7b3vPcy973vv/Xm/XvfrXutaT7+1F8zvvq5rrWulqpAkaboHjDsASdJkMkFIknqZICRJvUwQkqRe\nJghJUi8ThCSplwlCS0qSQ5JsvB/bn5nkD7dlTJPg/v5d5tj3cUkuG8W+NV4mCI1Mkg1JfpDku0lu\nTvLuJDuNO64pff+wVdVLq+qPxxXT/ZHkwCQfTXJbkm8m+VyS48cdlxYvE4RG7dlVtRPwRGA18Nox\nx7PoJdm+p+xg4BPAp4BHA48Afht45sJGp6XEBKEFUVWbgL8HHguQZI8kF7VfuuuTvGRq3SSnJvlQ\nkvcnuT3JFUkeP7C8kjx6YP7dSf6k77hJTk7ylbafa5I8p5X/LHAmcHCr4dzWt68kL2nxfbPFu8e0\nOF6a5MvtV/vbkmSGOOY6pz2SnJ9kS5Lrk/xOz7bvSfId4LieQ/w5cHZVvamqbq3Ouqo6clocr0xy\nS5LNg7WLJP+U5MUD8/eqXc3zXP88yWVJHt63XIuHCUILIsnewGHA51vRecBGYA/g+cCfJnnqwCaH\nAx8EdgXeB1yY5IFbceivAL8EPBx4PfCeJCur6lrgpcBnqmqnqtq5J+anAn8GHAmsBG5ocQ/6VeDn\ngce19Q6dJZbec0ryAOBvgX8H9gSeBpyU5NBp234I2Bl477Q4Hwoc3JbP5qfa32FP4ATgbUl2mWOb\nQbOea5IHJHlHW/70qvr2PPatCWSC0Khd2H6dX0bX/PGnLVk8GXh1Vf1nVV0JvBM4ZmC7dVX1oaq6\nEzgdeDBw0HwPXlUfrKqvV9WPqur9wJeBA4fc/DeBd1XVFVV1B3AKXY1j1cA6p1XVbVX1NeCTwAGz\n7G+mc/p5YEVVvaGqflhVXwXeARw1sO1nqurCdh4/mLbfXej+X948x/ncCbyhqu6sqo8C3wUeM8c2\ng2Y71wcC59Ilv2dX1ffnsV9NqPu0ZUrb2BFVdclgQWum+WZV3T5QfANdH8WUG6cmqupH7Q6cPZin\nJMcAvwesakU7AbsNufkewBUDcXw3yTfofoFvaMU3Daz//bb/mcx0TgXsMdXM1WwHfLpv2x7fAn5E\nV8v50izrfaOq7ppHvNPNdq6PBh4PHFhVP5zHPjXBrEFoHL4O7JrkYQNl+wCbBub3nppoTTB7te2g\n+8fpoQPr/lTfQZL8NN0v8ZcDj2jNSFcDU23ncw1l/HXgpwf2tyNd5++mGbeY3UzndCNwfVXtPPB5\nWFUdNrDtjLG2X+ufAZ63lXEBfI8h/qazuBY4Hvj7JPOplWiCmSC04KrqRuBfgT9L8uAkj6NrE3/P\nwGpPSvLcdsfOScAdwGfbsiuB30iyXZJnAP9thkPtSPcP6xaA1in72IHlNwN7JXnQDNufCxyf5IAk\nOwB/ClxeVRvmd8ZzntPngNuTvDrJQ9p5PTbJz89j368Cjkvy+0keAZDk8Umm95nM5ErguUke2m4A\nOGEexwagqs4FXgNckuRR891ek8cEoXE5mq7Z5+vABcDrpjVFfQT4H3TNJy8Entva7gFeATwbuI2u\nn+DCvgNU1TXAX9L9ur4Z+C/Avwys8gngi8BNSW7t2f4S4A+B8+na9x/FvfsF5qv3nKrqbroO4AOA\n64Fb6fpkhr4LqKr+FXhq+3w1yTeBNcBHh9zFm4Ef0v2dzmZaR/g84jgbeAPwiWl9NVqE4guDNGmS\nnAo8uqpeMO5YtpWleE5a+qxBSJJ6mSAkSb1sYpIk9bIGIUnqtagflNttt91q1apV4w5DkhaVdevW\n3VpVK+Zab1EniFWrVrF27dpxhyFJi0qSG4ZZb2RNTEn2TvLJNoLmF5O8opWfmmRTkivb57CBbU5p\nI2deN22gMknSAhtlDeIu4JVVdUUbUmFdko+3ZW+uqr8YXDnJ/nQPIf0c3fg0lyT5mfYQkSRpgY2s\nBlFVm6vqijZ9O91YLXvOssnhwHlVdUdVXQ+sZ/hRNyVJ29iC3MXUHrl/AnB5K3p5ki8kedfAePR7\ncu8RKzfSk1CSnJhkbZK1W7ZsGWHUkrS8jTxBpHsH8fnASVX1HeAMujFtDqAb3+Yv57O/qlpTVaur\navWKFXN2wkuSttJIE0R7A9j5wHur6sMAVXVzVd1dVT+iG4p5qhlpEwPDIdMNhby1wypLku6nUd7F\nFOAs4NqqOn2gfOXAas+hG58f4CLgqCQ7JNkX2I9uGGRJ0hiM8i6mJ9MNaXxVkitb2WuAo5McQDdO\n/wbgtwCq6otJPgBcQ3cH1Mu8g0mSxmdkCaKqLuOeN3cNmnF8+qp6I/DGUcUkSRreon6SWtLysurk\nv+st33DasxY4kuXBwfokSb1MEJKkXiYISVIvE4QkqZcJQpLUywQhSerlba6Slh1vlx2ONQhJUi8T\nhCSplwlCktTLPghJi559CqNhDUKS1MsEIUnqZYKQJPUyQUiSepkgJEm9TBCSpF4mCElSLxOEJKmX\nCUKS1MsnqSUtWTM9Ya3hWIOQJPWyBiFpbBxDabJZg5Ak9TJBSJJ6mSAkSb1MEJKkXiYISVIvE4Qk\nqZcJQpLUywQhSeplgpAk9RrZk9RJ9gbOAXYHClhTVW9JsivwfmAVsAE4sqq+lSTAW4DDgO8Dx1XV\nFaOKT5Km88nuextlDeIu4JVVtT9wEPCyJPsDJwOXVtV+wKVtHuCZwH7tcyJwxghjkyTNYWQJoqo2\nT9UAqup24FpgT+Bw4Oy22tnAEW36cOCc6nwW2DnJylHFJ0ma3YIM1pdkFfAE4HJg96ra3BbdRNcE\nBV3yuHFgs42tbPNAGUlOpKthsM8++4wsZknj4zDdk2HkndRJdgLOB06qqu8MLquqouufGFpVramq\n1VW1esWKFdswUknSoJEmiCQPpEsO762qD7fim6eajtr3La18E7D3wOZ7tTJJ0hiMLEG0u5LOAq6t\nqtMHFl0EHNumjwU+MlB+TDoHAd8eaIqSJC2wUfZBPBl4IXBVkitb2WuA04APJDkBuAE4si37KN0t\nruvpbnM9foSxSZLmMLIEUVWXAZlh8dN61i/gZaOKR9JozdaxvFyfI1jsfJJaktTLBCFJ6mWCkCT1\nMkFIknqZICRJvUwQkqReJghJUi8ThCSplwlCktTLBCFJ6mWCkCT1MkFIknotyBvlJC1vviFucbIG\nIUnqZYKQJPUyQUiSepkgJEm9TBCSpF4mCElSLxOEJKmXCUKS1MsEIUnqZYKQJPUyQUiSepkgJEm9\nTBCSpF4mCElSLxOEJKnXnO+DSPLrwD9U1e1JXgs8EfiTqrpi5NFJGpuZ3uGw4bRnLXAkGpdhahB/\n2JLDU4BfBs4CzhhtWJKkcRsmQdzdvp8FrKmqvwMeNLqQJEmTYJhXjm5K8nbgV4A3JdkB+y6kZcvX\nhy4fw/xDfyTwMeDQqroN2BX4/ZFGJUkau2FqEEdX1VlTM1W1OckrgH8cXViStpady9pWhqlBPC/J\nb07NJHkbsGKujZK8K8ktSa4eKDs1yaYkV7bPYQPLTkmyPsl1SQ6d74lIkratYWoQzwMuSvIj4BnA\nbVV1whDbvRv4K+CcaeVvrqq/GCxIsj9wFPBzwB7AJUl+pqruRpI0FjPWIJLsmmRX4CHAi4FXAbcD\nr2/ls6qqfwa+OWQchwPnVdUdVXU9sB44cMhtJUkjMFsNYh1QQAa+n9U+BTxyK4/58iTHAGuBV1bV\nt4A9gc8OrLOxld1HkhOBEwH22WefrQxBkjSXGRNEVe07guOdAfwxXYL5Y+AvgRfNZwdVtQZYA7B6\n9era1gFKS5Wd15qvYfogSPKLwKrB9atqet/CnKrq5oF9vgO4uM1uAvYeWHWvViZJGpNhxmL6G+BR\nwJXc81R1cd/O5zklWVlVm9vsc4CpO5wuAt6X5HS6Tur9gM/Nd/+SpG1nmBrEamD/qppXc06Sc4FD\ngN2SbAReBxyS5AC6BLMB+C2Aqvpikg8A1wB3AS/zDiZJGq9hEsTVwE8Bm+dacVBVHd1TfFZP2dT6\nbwTeOJ9jSJJGZ5gEsRtwTZLPAXdMFVbVr40sKknS2A2TIE4ddRCSpMkzZ4Koqk8tRCCSpMkyzF1M\nt9N1KkP3HogHAt+rqp8YZWCSpPEapgbxsKnpJKEbFuOgUQYlSRq/eb34pzoXAo62KklL3DBNTM8d\nmH0A3XMR/zmyiCRJE2GYu5iePTB9F90DboePJBpJ0sQYpg/i+IUIRJI0WYZpYloBvIT7DtY3r1FY\nJUmLyzBNTB8BPg1cwj2D9UmSlrhhEsRDq+rVI49EkjRRhrnN9eIkh408EknSRJmxBjHwBHWA1yS5\nA7izzZdPUkvS0jbbK0cfNtMySUvHTK8ilWZsYkpyaJLn95Q/L8mvjDYsSdK4zdYH8UdA30iunwLe\nMJpwJEmTYra7mHaoqi3TC6vq1iQ7jjAmSUOwaUijNluC+Ikk21fVXYOFSR4IPGS0YUmaYiLQuMzW\nxPRh4B2DtYUkOwFntmWSpCVstgTxWuBm4IYk65KsA64HtrRlkqQlbLbbXO8CTk7yeuDRrXh9Vf1g\nQSKTJI3VMKO5/gC4agFikSRNkHm9UU6StHzM9qDck9v3DgsXjiRpUsxWg3hr+/7MQgQiSZoss/VB\n3JlkDbBnkrdOX1hVvzO6sKTlx+cdNGlmSxC/CvwycCiwbmHCkSRNitluc70VOC/JtVX17wsYkyRp\nAgxzF9M3klyQ5Jb2OT/JXiOPTJI0VsO8cvSvgfcBv97mX9DKHPJb2gr2NWixGKYG8ZNV9ddVdVf7\nvBtYMeK4JEljNkyCuDXJC5Js1z4vAL4x6sAkSeM1TIJ4EXAkcBOwGXg+cPxcGyV5V+uzuHqgbNck\nH0/y5fa9SytPkrcmWZ/kC0meuHWnI0naVuZMEFV1Q1X9WlWtqKqfrKojquprQ+z73cAzppWdDFxa\nVfsBl7Z5gGcC+7XPicAZw56AJGk0RjYWU1X9M/DNacWHA2e36bOBIwbKz6nOZ4Gdk6wcVWySpLkt\n9GB9u1fV5jZ9E7B7m94TuHFgvY2tTJI0JmMbzbWqCqj5bpfkxCRrk6zdsuU+r8yWJG0jcyaIJK8d\nmL6/I7vePNV01L5vaeWbgL0H1turld1HVa2pqtVVtXrFCu+2laRRmW2471cnOZjurqUp93dk14uA\nY9v0scBHBsqPaXczHQR8e6ApSpI0BrM9Sf0luqenH5nk023+EUkeU1XXzbXjJOcChwC7JdkIvA44\nDfhAkhOAG+hunwX4KHAYsB74PkPcRitJGq3ZEsRtwGvo/pE/BPhZ4Ol076l+TFX94mw7rqqjZ1j0\ntJ51C3jZEPFKkhbIbAniUOCPgEcBpwNfAL5XVf66l6RlYMY+iKp6TVU9DdgA/A2wHbAiyWVJ/naB\n4pMkjckwo7l+rKrWAmuT/HZVPSXJbqMOTJI0XnMmiKp61cDsca3s1lEFJC2EmYbc3nDasxY4Emly\nDVOD+DHfLCfdl8lGS9XYnqSWJE22edUgJA3PN8dpsbMGIUnqZQ1CGoK1AS1H1iAkSb1MEJKkXiYI\nSVIvE4QkqZcJQpLUywQhSerlba5aVBzWQlo41iAkSb1MEJKkXiYISVIv+yA0kbbV0BYOkSFtPWsQ\nkqReJghJUi+bmKQBNklJ97AGIUnqZYKQJPWyiUljZZOONLmsQUiSepkgJEm9bGLSjGZr/plpcLxx\nDaZnU5W07VmDkCT1MkFIknqZICRJvUwQkqRedlJrq9gpLC19Y0kQSTYAtwN3A3dV1eokuwLvB1YB\nG4Ajq+pb44hPkjTeJqb/XlUHVNXqNn8ycGlV7Qdc2uYlSWMySX0QhwNnt+mzgSPGGIskLXvjShAF\n/GOSdUlObGW7V9XmNn0TsHvfhklOTLI2ydotW7YsRKyStCyNq5P6KVW1KclPAh9P8qXBhVVVSapv\nw6paA6wBWL16de86kqT7byw1iKra1L5vAS4ADgRuTrISoH3fMo7YJEmdBa9BJNkReEBV3d6mnw68\nAbgIOBY4rX1/ZKFjW668ZVVSn3E0Me0OXJBk6vjvq6p/SPJvwAeSnADcABw5htg0IiYhafFZ8ARR\nVV8FHt9T/g3gaQsdjySp3yTd5ipJmiAmCElSL8diWsTG9XIeScuDNQhJUi8ThCSpl01Mi4C3iEoa\nBxPEMmKikTQfNjFJknqZICRJvUwQkqReJghJUi87qZcgO6MlbQsmiDHwCWhJi4FNTJKkXiYISVIv\nm5gmiH0HkiaJNQhJUi8ThCSplwlCktTLPogRsk9B0mJmDUKS1MsaxDz4gJuk5cQEsQ3YlCRpKbKJ\nSZLUywQhSeplgpAk9TJBSJJ6LdtO6tk6lr0rSZKsQUiSZmCCkCT1MkFIknqZICRJvZZtJ/VsfDJa\nkqxBSJJmMHE1iCTPAN4CbAe8s6pOG3NIkpa55TpQ50QliCTbAW8DfgXYCPxbkouq6prxRiZJ97Ut\nE8ckJqGJShDAgcD6qvoqQJLzgMMBE4SkRWNbPog7zsQxaQliT+DGgfmNwC8MrpDkRODENvvdJNdt\n5bF2A27dym0n2VI8r6V4TrA0z2spnhNsw/PKm7bFXu73fn56mJUmLUHMqarWAGvu736SrK2q1dsg\npImyFM9rKZ4TLM3zWornBEv3vOYyaXcxbQL2Hpjfq5VJkhbYpCWIfwP2S7JvkgcBRwEXjTkmSVqW\nJqqJqaruSvJy4GN0t7m+q6q+OKLD3e9mqgm1FM9rKZ4TLM3zWornBEv3vGaVqhp3DJKkCTRpTUyS\npAlhgpAk9VqWCSLJM5Jcl2R9kpPHHc90SfZO8skk1yT5YpJXtPJdk3w8yZfb9y6tPEne2s7nC0me\nOLCvY9v6X05y7ED5k5Jc1bZ5a5Is0Lltl+TzSS5u8/smubzF8f52cwJJdmjz69vyVQP7OKWVX5fk\n0IHysVzXJDsn+VCSLyW5NsnBi/1aJfnd9t/e1UnOTfLgxXitkrwryS1Jrh4oG/m1mekYi05VLasP\nXef3V4BHAg8C/h3Yf9xxTYtxJfDENv0w4D+A/YH/DZzcyk8G3tSmDwP+HghwEHB5K98V+Gr73qVN\n79KWfa6tm7btMxfo3H4PeB9wcZv/AHBUmz4T+O02/T+BM9v0UcD72/T+7ZrtAOzbruV247yuwNnA\ni9v0g4CdF/O1ontg9XrgIQPX6LjFeK2A/wo8Ebh6oGzk12amYyy2z9gDWPAThoOBjw3MnwKcMu64\n5oj5I3TjU10HrGxlK4Hr2vTbgaMH1r+uLT8aePtA+dtb2UrgSwPl91pvhOexF3Ap8FTg4vY/1a3A\n9tOvDd2dbAe36e3bepl+vabWG9d1BR7e/jHNtPJFe624Z0SDXdvf/mLg0MV6rYBV3DtBjPzazHSM\nxfZZjk1MfcN57DmmWObUqutPAC4Hdq+qzW3RTcDubXqmc5qtfGNP+aj9H+BVwI/a/COA26rqrp44\nfhx7W/7ttv58z3XU9gW2AH/dms7emWRHFvG1qqpNwF8AXwM20/3t17H4r9WUhbg2Mx1jUVmOCWLR\nSLITcD5wUlV9Z3BZdT9NFs09ykl+FbilqtaNO5ZtbHu6JowzquoJwPfomhR+bBFeq13oBsncF9gD\n2BF4xliDGpGFuDaL7foPWo4JYlEM55HkgXTJ4b1V9eFWfHOSlW35SuCWVj7TOc1WvldP+Sg9Gfi1\nJBuA8+iamd4C7Jxk6oHNwTh+HHtb/nDgG8z/XEdtI7Cxqi5v8x+iSxiL+Vr9MnB9VW2pqjuBD9Nd\nv8V+raYsxLWZ6RiLynJMEBM/nEe7E+Is4NqqOn1g0UXA1B0Ux9L1TUyVH9PuwjgI+Har3n4MeHqS\nXdqvwqfTtf1uBr6T5KB2rGMG9jUSVXVKVe1VVavo/uafqKrfBD4JPH+Gc5o61+e39auVH9XunNkX\n2I+uo3As17WqbgJuTPKYVvQ0uuHpF+21omtaOijJQ9sxp85pUV+rAQtxbWY6xuIy7k6QcXzo7lb4\nD7o7Kf5g3PH0xPcUuirpF4Ar2+cwunbdS4EvA5cAu7b1Q/eipa8AVwGrB/b1ImB9+xw/UL4auLpt\n81dM62Qd8fkdwj13MT2S7h+N9cAHgR1a+YPb/Pq2/JED2/9Bi/s6Bu7oGdd1BQ4A1rbrdSHdnS6L\n+loBrwe+1I77N3R3Ii26awWcS9ePciddbe+Ehbg2Mx1jsX0cakOS1Gs5NjFJkoZggpAk9TJBSJJ6\nmSAkSb1MEJKkXiYILTtJ7k5yZRup9INJHrqN939ckr+aY51DkvziwPxLkxyzLeOQ7i8ThJajH1TV\nAVX1WOCHwEvHEMMhwI8TRFWdWVXnjCEOaUYmCC13nwYeDZDk91qt4uokJ7WyVene8/DedO96+NBU\njSPJhiS7tenVSf5p+s6TPDvdOxI+n+SSJLu3ARhfCvxuq8n8UpJTk/yvts0BST7b3klwQe55X8E/\nJXlTks8l+Y8kvzT6P4+WMxOElq02btAzgauSPAk4HvgFuvH9X5LkCW3VxwD/r6p+FvgO3fsPhnUZ\ncFB1A/mdB7yqqjbQvU/hza0m8+lp25wDvLqqHkf3RO/rBpZtX1UHAidNK5e2OROElqOHJLmSbniM\nr9GNe/UU4IKq+l5VfZdugLqpX+g3VtW/tOn3tHWHtRfwsSRXAb8P/NxsKyd5OLBzVX2qFZ1N99Kb\nKVMDN66je8+BNDLbz72KtOT8oKoOGCzI7G/xnD4ezdT8XdzzI+vBM2z7f4HTq+qiJIcAp84r0vu6\no33fjf//asSsQUidTwNHtBFMdwSe08oA9klycJv+DbpmI4ANwJPa9PNm2O/DuWcI6GMHym+ne53s\nvVTVt4FvDfQvvBD41PT1pIVggpCAqroCeDfdaKSXA++sqs+3xdcBL0tyLd1IrWe08tcDb0mylu4X\nfZ9TgQ8mWUf3Ks4pfws8Z6qTeto2xwJ/nuQLdCPFvuH+nJu0tRzNVZpFu+Po4nZLrLSsWIOQJPWy\nBiFJ6mUNQpLUywQhSeplgpAk9TJBSJJ6mSAkSb3+P6HC8/9lpHlRAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7feacf51d2e8>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(chunks['2010 Census Population'], bins=50)\n",
    "plt.title(\"Population per Chunk\")\n",
    "plt.xlabel(\"Population\")\n",
    "plt.ylabel(\"# of Chunks\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "308739931"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "chunks['2010 Census Population'].sum()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(4444, 2)"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "chunks.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style>\n",
       "    .dataframe thead tr:only-child th {\n",
       "        text-align: right;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: left;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>2010 Census Population</th>\n",
       "      <th>Prefixes</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Chunk</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>010_0</th>\n",
       "      <td>210857</td>\n",
       "      <td>[0100, 0101, 0102, 0103]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>010_1</th>\n",
       "      <td>238700</td>\n",
       "      <td>[0104, 0105, 0106, 0107, 0108]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>010_2</th>\n",
       "      <td>19818</td>\n",
       "      <td>[0109]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>015_0</th>\n",
       "      <td>244841</td>\n",
       "      <td>[0150, 0151, 0152, 0153, 0154, 0155]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>015_1</th>\n",
       "      <td>128878</td>\n",
       "      <td>[0156, 0157, 0158, 0159]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>017_0</th>\n",
       "      <td>207777</td>\n",
       "      <td>[0170, 0171, 0172, 0173, 0174]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>017_1</th>\n",
       "      <td>171250</td>\n",
       "      <td>[0175, 0176, 0177, 0178, 0179]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>018_0</th>\n",
       "      <td>184413</td>\n",
       "      <td>[0180, 0181, 0182]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>018_1</th>\n",
       "      <td>227505</td>\n",
       "      <td>[0183, 0184]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>018_2</th>\n",
       "      <td>210817</td>\n",
       "      <td>[0185, 0186, 0187]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>018_3</th>\n",
       "      <td>90391</td>\n",
       "      <td>[0188, 0189]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>019_0</th>\n",
       "      <td>231079</td>\n",
       "      <td>[0190, 0191, 0192]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>019_1</th>\n",
       "      <td>245235</td>\n",
       "      <td>[0193, 0194, 0195, 0196, 0197, 0198, 0199]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>01_0</th>\n",
       "      <td>168835</td>\n",
       "      <td>[011]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>01_1</th>\n",
       "      <td>213800</td>\n",
       "      <td>[012, 013]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>01_2</th>\n",
       "      <td>211812</td>\n",
       "      <td>[014]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>01_3</th>\n",
       "      <td>188402</td>\n",
       "      <td>[016]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>020_0</th>\n",
       "      <td>184275</td>\n",
       "      <td>[0200, 0201, 0202, 0203, 0204]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>020_1</th>\n",
       "      <td>216096</td>\n",
       "      <td>[0205, 0206, 0207, 0208, 0209]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0212_0</th>\n",
       "      <td>243585</td>\n",
       "      <td>[02120, 02121, 02122, 02123, 02124, 02125, 021...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0212_1</th>\n",
       "      <td>16439</td>\n",
       "      <td>[02129]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0213_0</th>\n",
       "      <td>220173</td>\n",
       "      <td>[02130, 02131, 02132, 02133, 02134, 02135, 021...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0213_1</th>\n",
       "      <td>36349</td>\n",
       "      <td>[02139]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>021_0</th>\n",
       "      <td>136539</td>\n",
       "      <td>[0210, 0211]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>021_1</th>\n",
       "      <td>207555</td>\n",
       "      <td>[0214]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>021_2</th>\n",
       "      <td>220030</td>\n",
       "      <td>[0215, 0216]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>021_3</th>\n",
       "      <td>203272</td>\n",
       "      <td>[0217, 0218, 0219]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>023_0</th>\n",
       "      <td>241125</td>\n",
       "      <td>[0230, 0231, 0232, 0233, 0234]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>023_1</th>\n",
       "      <td>199938</td>\n",
       "      <td>[0235, 0236, 0237, 0238, 0239]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>024_0</th>\n",
       "      <td>242552</td>\n",
       "      <td>[0240, 0241, 0242, 0243, 0244, 0245, 0246]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>981_3</th>\n",
       "      <td>207457</td>\n",
       "      <td>[9815, 9816, 9817, 9818, 9819]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>982_0</th>\n",
       "      <td>168159</td>\n",
       "      <td>[9820, 9821]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>982_1</th>\n",
       "      <td>230928</td>\n",
       "      <td>[9822, 9823]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>982_2</th>\n",
       "      <td>138281</td>\n",
       "      <td>[9824, 9825, 9826]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>982_3</th>\n",
       "      <td>249350</td>\n",
       "      <td>[9827, 9828]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>982_4</th>\n",
       "      <td>90767</td>\n",
       "      <td>[9829]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>983_0</th>\n",
       "      <td>240218</td>\n",
       "      <td>[9830, 9831, 9832, 9833, 9834]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>983_1</th>\n",
       "      <td>148730</td>\n",
       "      <td>[9835, 9836]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>983_2</th>\n",
       "      <td>162860</td>\n",
       "      <td>[9837]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>983_3</th>\n",
       "      <td>161800</td>\n",
       "      <td>[9838, 9839]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>984_0</th>\n",
       "      <td>214804</td>\n",
       "      <td>[9840, 9841, 9842, 9843]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>984_1</th>\n",
       "      <td>185311</td>\n",
       "      <td>[9844, 9845, 9846, 9847, 9848, 9849]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>985_0</th>\n",
       "      <td>241456</td>\n",
       "      <td>[9850, 9851, 9852]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>985_1</th>\n",
       "      <td>234098</td>\n",
       "      <td>[9853, 9854, 9855, 9856, 9857, 9858, 9859]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>986_0</th>\n",
       "      <td>227225</td>\n",
       "      <td>[9860, 9861, 9862, 9863, 9864, 9865]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>986_1</th>\n",
       "      <td>191810</td>\n",
       "      <td>[9866, 9867]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>986_2</th>\n",
       "      <td>154295</td>\n",
       "      <td>[9868, 9869]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>989_0</th>\n",
       "      <td>203780</td>\n",
       "      <td>[9890, 9891, 9892, 9893]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>989_1</th>\n",
       "      <td>80875</td>\n",
       "      <td>[9894, 9895, 9896, 9897, 9898, 9899]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>98_0</th>\n",
       "      <td>219761</td>\n",
       "      <td>[987, 988]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>992_0</th>\n",
       "      <td>218781</td>\n",
       "      <td>[9920]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>992_1</th>\n",
       "      <td>124075</td>\n",
       "      <td>[9921, 9922, 9923, 9924, 9925, 9926, 9927, 992...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>993_0</th>\n",
       "      <td>224771</td>\n",
       "      <td>[9930, 9931, 9932, 9933, 9934]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>993_1</th>\n",
       "      <td>124213</td>\n",
       "      <td>[9935, 9936, 9937, 9938, 9939]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>995_0</th>\n",
       "      <td>184878</td>\n",
       "      <td>[9950]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>995_1</th>\n",
       "      <td>135782</td>\n",
       "      <td>[9951, 9952, 9953, 9954, 9955, 9956, 9957, 995...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>99_0</th>\n",
       "      <td>142008</td>\n",
       "      <td>[990]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>99_1</th>\n",
       "      <td>136934</td>\n",
       "      <td>[991, 994]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>99_2</th>\n",
       "      <td>184083</td>\n",
       "      <td>[996]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>99_3</th>\n",
       "      <td>205187</td>\n",
       "      <td>[997, 998, 999]</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>1780 rows × 2 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "        2010 Census Population  \\\n",
       "Chunk                            \n",
       "010_0                   210857   \n",
       "010_1                   238700   \n",
       "010_2                    19818   \n",
       "015_0                   244841   \n",
       "015_1                   128878   \n",
       "017_0                   207777   \n",
       "017_1                   171250   \n",
       "018_0                   184413   \n",
       "018_1                   227505   \n",
       "018_2                   210817   \n",
       "018_3                    90391   \n",
       "019_0                   231079   \n",
       "019_1                   245235   \n",
       "01_0                    168835   \n",
       "01_1                    213800   \n",
       "01_2                    211812   \n",
       "01_3                    188402   \n",
       "020_0                   184275   \n",
       "020_1                   216096   \n",
       "0212_0                  243585   \n",
       "0212_1                   16439   \n",
       "0213_0                  220173   \n",
       "0213_1                   36349   \n",
       "021_0                   136539   \n",
       "021_1                   207555   \n",
       "021_2                   220030   \n",
       "021_3                   203272   \n",
       "023_0                   241125   \n",
       "023_1                   199938   \n",
       "024_0                   242552   \n",
       "...                        ...   \n",
       "981_3                   207457   \n",
       "982_0                   168159   \n",
       "982_1                   230928   \n",
       "982_2                   138281   \n",
       "982_3                   249350   \n",
       "982_4                    90767   \n",
       "983_0                   240218   \n",
       "983_1                   148730   \n",
       "983_2                   162860   \n",
       "983_3                   161800   \n",
       "984_0                   214804   \n",
       "984_1                   185311   \n",
       "985_0                   241456   \n",
       "985_1                   234098   \n",
       "986_0                   227225   \n",
       "986_1                   191810   \n",
       "986_2                   154295   \n",
       "989_0                   203780   \n",
       "989_1                    80875   \n",
       "98_0                    219761   \n",
       "992_0                   218781   \n",
       "992_1                   124075   \n",
       "993_0                   224771   \n",
       "993_1                   124213   \n",
       "995_0                   184878   \n",
       "995_1                   135782   \n",
       "99_0                    142008   \n",
       "99_1                    136934   \n",
       "99_2                    184083   \n",
       "99_3                    205187   \n",
       "\n",
       "                                                 Prefixes  \n",
       "Chunk                                                      \n",
       "010_0                            [0100, 0101, 0102, 0103]  \n",
       "010_1                      [0104, 0105, 0106, 0107, 0108]  \n",
       "010_2                                              [0109]  \n",
       "015_0                [0150, 0151, 0152, 0153, 0154, 0155]  \n",
       "015_1                            [0156, 0157, 0158, 0159]  \n",
       "017_0                      [0170, 0171, 0172, 0173, 0174]  \n",
       "017_1                      [0175, 0176, 0177, 0178, 0179]  \n",
       "018_0                                  [0180, 0181, 0182]  \n",
       "018_1                                        [0183, 0184]  \n",
       "018_2                                  [0185, 0186, 0187]  \n",
       "018_3                                        [0188, 0189]  \n",
       "019_0                                  [0190, 0191, 0192]  \n",
       "019_1          [0193, 0194, 0195, 0196, 0197, 0198, 0199]  \n",
       "01_0                                                [011]  \n",
       "01_1                                           [012, 013]  \n",
       "01_2                                                [014]  \n",
       "01_3                                                [016]  \n",
       "020_0                      [0200, 0201, 0202, 0203, 0204]  \n",
       "020_1                      [0205, 0206, 0207, 0208, 0209]  \n",
       "0212_0  [02120, 02121, 02122, 02123, 02124, 02125, 021...  \n",
       "0212_1                                            [02129]  \n",
       "0213_0  [02130, 02131, 02132, 02133, 02134, 02135, 021...  \n",
       "0213_1                                            [02139]  \n",
       "021_0                                        [0210, 0211]  \n",
       "021_1                                              [0214]  \n",
       "021_2                                        [0215, 0216]  \n",
       "021_3                                  [0217, 0218, 0219]  \n",
       "023_0                      [0230, 0231, 0232, 0233, 0234]  \n",
       "023_1                      [0235, 0236, 0237, 0238, 0239]  \n",
       "024_0          [0240, 0241, 0242, 0243, 0244, 0245, 0246]  \n",
       "...                                                   ...  \n",
       "981_3                      [9815, 9816, 9817, 9818, 9819]  \n",
       "982_0                                        [9820, 9821]  \n",
       "982_1                                        [9822, 9823]  \n",
       "982_2                                  [9824, 9825, 9826]  \n",
       "982_3                                        [9827, 9828]  \n",
       "982_4                                              [9829]  \n",
       "983_0                      [9830, 9831, 9832, 9833, 9834]  \n",
       "983_1                                        [9835, 9836]  \n",
       "983_2                                              [9837]  \n",
       "983_3                                        [9838, 9839]  \n",
       "984_0                            [9840, 9841, 9842, 9843]  \n",
       "984_1                [9844, 9845, 9846, 9847, 9848, 9849]  \n",
       "985_0                                  [9850, 9851, 9852]  \n",
       "985_1          [9853, 9854, 9855, 9856, 9857, 9858, 9859]  \n",
       "986_0                [9860, 9861, 9862, 9863, 9864, 9865]  \n",
       "986_1                                        [9866, 9867]  \n",
       "986_2                                        [9868, 9869]  \n",
       "989_0                            [9890, 9891, 9892, 9893]  \n",
       "989_1                [9894, 9895, 9896, 9897, 9898, 9899]  \n",
       "98_0                                           [987, 988]  \n",
       "992_0                                              [9920]  \n",
       "992_1   [9921, 9922, 9923, 9924, 9925, 9926, 9927, 992...  \n",
       "993_0                      [9930, 9931, 9932, 9933, 9934]  \n",
       "993_1                      [9935, 9936, 9937, 9938, 9939]  \n",
       "995_0                                              [9950]  \n",
       "995_1   [9951, 9952, 9953, 9954, 9955, 9956, 9957, 995...  \n",
       "99_0                                                [990]  \n",
       "99_1                                           [991, 994]  \n",
       "99_2                                                [996]  \n",
       "99_3                                      [997, 998, 999]  \n",
       "\n",
       "[1780 rows x 2 columns]"
      ]
     },
     "execution_count": 21,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "chunks_250000 = build_chunks(250000)\n",
    "chunks_250000"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.text.Text at 0x7feacf55c320>"
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEWCAYAAABrDZDcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGw1JREFUeJzt3Xu4HFWd7vHvC4EgEElCMjEkxHA7KHrkFhEGVAYQBERA\nGA6oECLKwTMcRWeEgDhc9HFgdHB0hgNGQaIgKOHqZWTk+oA6MAmEa0AuJpCQhAQIBEQg8Dt/1GpS\nbLp7997Z1b13r/fzPP3s6lVrVa3V3bt/tVZVr1JEYGZm+Vqr0xUwM7POciAwM8ucA4GZWeYcCMzM\nMudAYGaWOQcCM7PMORBYv0jaXdLCNSh/vqSvDWSdBoM1fV162fbRkm7rQ/5tJM2WpCrq026Shkt6\nUNLYTtel2zgQdAFJ8yW9JOkFSUslXSRpw07Xq6beF1hEHBcRX+9UndaEpJ0k/VrSCknPSLpD0rRO\n16uOrwPfjohIX6IXSFogaaWkuZL2rWWUNFlSpM9Q7fG10vrhki6U9LykJZK+3GzHkr6U8j2fyg3v\nsa+bJP05fbHv1UrZiHgZuBCYPkCvjyUOBN3jgIjYENgBmAKc2uH6DHmShtVJ2wW4EbgF2BLYGPg8\nsG/PvJ0kaTzwN8DVKWkY8ATwYWAjis/HzyVN7lF0ZERsmB7lQH06sBXwzrTdEyV9tMG+96H4st4z\n5d8cOKOU5VLgLorX7qvArNpRfgtlfwpMLQcWGwAR4ccQfwDzgb1Kz78F/DItbwJcCzwDPAJ8rpTv\ndGAW8DNgJXAnsG1pfQBblp5fBHwjLe8OLCytmw48mrbzAHBwSn838BfgNeAFYEXPbaXnn0v1eybV\nd5Me9TgOeBhYAZwLqMFr0VubNgGuAJYBfwK+UKfsxcDzwGfrbP824Nwm78XuwELg74GngMXAtNL6\nm8vbBY4GbmulrXXyfivVZ6M69TgKuL6Xz809wCFpeXLa97AGeZ8E9i49/zpwWYO8PwW+WXq+J7Ak\nLf8P4GVgRGn9rcBxvZUtpT0MfLjT/3fd9HCPoMtI2hTYj+KIC+Ayii+mTYBDgW9K2qNU5EDgcmA0\nxT/h1ZLW6ceuHwU+SHG0eQZwsaTxETGP4ovtD1EcZY6sU+c9gH8CDgPGAwtSvcs+BrwfeF/Kt0+T\nutRtk6S1gF8AdwMTKL5kTkhHoeWys4CRwCU96rk+sEta38w70uswATgGOFfSqF7KlDVtq6S1JP0g\nrd87Ip6rs43/CTzUaAeSxlF8Kd/fY9UCSQsl/UjSmJR3FMX7cncp393Aexps/j118o6TtHFa91hE\nrGywrWZla+YB2zZqm/WdA0H3uFrSCoojxFsovvA3BXYFToqIv0TEXOCHFEeLNXMiYlZEvAqcA6wH\n7NzXnUfE5RHxZES8HhE/ozhq26nF4p8CLoyIO6MYBz4Z2KXHsMVZEbEiIh4HbgK2a7K9Rm16PzA2\nIs6MiFci4jHgB8DhpbJ/iIirUzte6rHdURT/M4t7ac+rwJkR8WpE/JqiJ7R1L2XKmrV1HYqhldEU\nw4F/brCNkRQ9ordIgf4SYGZEPJiSl1O8Pu8EdgRGsDoQ1s43lQPOcylPPRvWyUvK33Ndz201K1uz\nkqJ9NkDeMgZqQ9ZBEXF9OUHSJsAzPY6+FlCcQ6h5orYQEa+nK1426evOJR0FfJliiAGKf+gxLRbf\nhGIIp1aPFyQ9TXFEPT8lLynl/zOrv5zqadSmADZJAbNmbYqhibeUreNZ4HWKo+MHm+R7OiJW9aG+\nPTVr65YUR8M7RcQrvdT1LV/UqVf0E+AV4PhaekS8AMxOT5dKOh5YLGkERSADeDvFMF9tuW6gSfnf\nXnpeW15ZZ13PbTUrWzOCYtjMBoh7BN3tSWB0+meumQQsKj3ftLaQviQmpnJQfAmtX8r7jno7kfRO\niiPr44GN0/DPfUDtssXeprh9kuJItLa9DShOJC5qWKK5Rm16AvhTRIwsPUZExH6lsg3rmo6+/wAc\n0s96AbxIC69pE/OAacB/SGrWy7iHYujnDeky0guAcRTnBl5tUr72OqwVEc9S9ILKwzHb8tZhpZr7\n6+RdGhFPp3Wb9/hMlrfVrGzNu3nz8JGtIQeCLhYRTwC/B/5J0nqS3kcxZn1xKduOkj6RrpA5geJE\n3n+ldXOBT0paO10h8uEGu9qA4otjGUC6lPK9pfVLgYmS1m1Q/lJgmqTt0tUg3wRuj4j5fWtxr226\nA1gp6SRJb0vteq+k9/dh2ycCR0v6Sm3cWtK2knqe02hkLvAJSetL2pLi/eiTiLgUOAW4XtIWDbL9\nFthB0nqltPMovkQP6DnsJekDkrZO5x82Br4H3Fw6//Bj4FRJoyS9i+Lk/kWl8iFp91LeY1T8jmEk\nxRVKF6W6/zG9Bqelz+TBFOc6ruitbNrPBIphsdpn1AaAA0H3O4JiuOZJ4CrgtB5DSNcA/4tiKOFI\n4BOlI8UvAgdQdMM/xepLEd8kIh4A/oXiaHkpxYnK35Wy3EhxpLdE0vI65a8HvkbxZbAY2II3j9v3\nVd02RcRrFCdit6O4Ymg5xTmTjVrdcET8HtgjPR6T9AwwA/h1i5v4DsWwzFJgJj1OSPehHjOBM4Eb\n61wCSkQspXjdD4Q3em3/m6LtS0q/FfhUKrI58BuKIZj7KILnEaVNnkZxQcACinNQ34qI36Rtb5rK\n3Zv2/RvgnynObzyeypxW2tbhFMOTzwJnAYdGxLIWy36S4tzGyy2/WNar2mVpliFJp1NcHvrpTtdl\noHRjm/pL0jYUwWanqPAfXdKngfdExMlV7SPtZzjFkNCHIuKpKveVG58sNutSqafWl2Gv/u7n4t5z\nDch+Xgbe1Y595cZDQ2ZmmfPQkJlZ5twjMDPL3JA4RzBmzJiYPHlyp6thZjakzJkzZ3lE9Dpt95AI\nBJMnT2b27Nm9ZzQzszdIWtBKPg8NmZllzoHAzCxzDgRmZplzIDAzy5wDgZlZ5hwIzMwy50BgZpY5\nBwIzs8w5EJiZZW5I/LLYzGwomTz9V3XT55+1f5tr0hr3CMzMMucegZlZMtSO5AeKewRmZplzIDAz\ny5wDgZlZ5hwIzMwy50BgZpY5BwIzs8z58lEzsw7r9GWr7hGYmWXOgcDMLHMOBGZmmXMgMDPLnAOB\nmVnmHAjMzDLnQGBmljkHAjOzzDkQmJllzoHAzCxzlU4xIelLwGeBAO4FpgHjgcuAjYE5wJER8UqV\n9TAzWxOdngKiapX1CCRNAL4ATImI9wJrA4cDZwPfiYgtgWeBY6qqg5mZ9a7qoaFhwNskDQPWBxYD\newCz0vqZwEEV18HMzJqobGgoIhZJ+jbwOPAS8J8UQ0ErImJVyrYQmFCvvKRjgWMBJk2aVFU1zcza\nptEQU6dVOTQ0CjgQ2AzYBNgA+Gir5SNiRkRMiYgpY8eOraiWZmZW5dDQXsCfImJZRLwKXAnsCoxM\nQ0UAE4FFFdbBzMx6UWUgeBzYWdL6kgTsCTwA3AQcmvJMBa6psA5mZtaLKs8R3C5pFnAnsAq4C5gB\n/Aq4TNI3UtoFVdXBzPJW9WWfg3XMv68q/R1BRJwGnNYj+TFgpyr3a2ZmrfMvi83MMueb15tZdrpl\nSGeguEdgZpY5BwIzs8w5EJiZZc6BwMwscw4EZmaZcyAwM8ucA4GZWeYcCMzMMudAYGaWOQcCM7PM\nORCYmWXOgcDMLHMOBGZmmXMgMDPLnAOBmVnmHAjMzDLnQGBmljkHAjOzzDkQmJllzoHAzCxzDgRm\nZplzIDAzy5wDgZlZ5hwIzMwy50BgZpY5BwIzs8w5EJiZZc6BwMwscw4EZmaZG9bpCpiZranJ03/V\n6SoMae4RmJllzoHAzCxzDgRmZplzIDAzy5wDgZlZ5ioNBJJGSpol6UFJ8yTtImm0pN9Kejj9HVVl\nHczMrLmqewTfBX4TEe8CtgXmAdOBGyJiK+CG9NzMzDqkskAgaSPgQ8AFABHxSkSsAA4EZqZsM4GD\nqqqDmZn1rsoewWbAMuBHku6S9ENJGwDjImJxyrMEGFevsKRjJc2WNHvZsmUVVtPMLG+9BgJJfytp\nRFo+VdKVknZoYdvDgB2A8yJie+BFegwDRUQAUa9wRMyIiCkRMWXs2LEt7M7MzPqjlR7B1yJipaTd\ngL0ohnrOa6HcQmBhRNyens+iCAxLJY0HSH+f6nu1zcxsoLQSCF5Lf/cHZkTEr4B1eysUEUuAJyRt\nnZL2BB4ArgWmprSpwDV9qrGZmQ2oViadWyTp+8BHgLMlDaf1cwv/F7hE0rrAY8C0VPbnko4BFgCH\n9b3aZpYjTy5XjVYCwWHAR4FvR8SKNJzzlVY2HhFzgSl1Vu3ZehXNzKxKrRzZHxERV0bEwwDpip89\nqq2WmZm1Sys9gkMk/SUiLgGQdC6wXrXVMjOzdmkpEADXSnqdYohoRUQcU221zMysXRoGAkmjS08/\nC1wN/A44Q9LoiHim6sqZmVn1mvUI5lD82Eulv/unRwCbV147MzOrXMNAEBGbtbMiZjY0NLqEc/5Z\n+1e6fatOSzevl/TXwORy/oj4cUV1MjOzNuo1EEj6CbAFMJfVvzIOwIHAzKwLtNIjmAJskyaIMzOz\nLtPKD8ruA95RdUXMzKwzWukRjAEekHQH8HItMSI+XlmtzMysbVoJBKdXXQkzM+ucXgNBRNzSjoqY\nWbWqvuzThq5Wrhpayeq7iK0LrAO8GBFvr7JiZmbWHq30CEbUliWJ4ubzO1dZKTMbetzjGLr6dPP6\nKFwN7FNRfczMrM1aGRr6ROnpWhS/K/hLZTUyM7O2auWqoQNKy6uA+RTDQ2Zm1gVaOUcwrR0VMTOz\nzmhlaGgs8DneOuncZ6qrlpl1C59EHvxaGRq6BrgVuJ7Vk86ZmVmXaCUQrB8RJ1VeEzMz64hWLh/9\npaT9Kq+JmZl1RLN7Ftd+USzgFEkvA6+m5+FfFpuZdYdmt6oc0WidmZl1j4ZDQ5L2kXRonfRDJH2k\n2mqZmVm7NDtH8I9AvZlHbwHOrKY6ZmbWbs0CwfCIWNYzMSKWAxtUVyUzM2unZpePvl3SsIhYVU6U\ntA7wtmqrZWbdrtEPzaz9mvUIrgR+IOmNo39JGwLnp3VmZtYFmgWCU4GlwAJJcyTNAf4ELEvrzMys\nCzS7fHQVMF3SGcCWKfmRiHipLTUzM7O2aGX20ZeAe9tQFzMz64A+3aHMzMy6T7MflO2a/g5vX3XM\nzKzdmg0NfQ/YEfgDsEN7qmNm7eb7BVizQPCqpBnABEnf67kyIr5QXbXMzKxdmgWCjwF7AfsAc/q7\nA0lrA7OBRRHxMUmbAZcBG6ftHhkRr/R3+2ZmtmaaXT66HLhM0ryIuHsN9vFFYB5Qm7b6bOA7EXGZ\npPOBY4Dz1mD7Zma2Blq5auhpSVdJeio9rpA0sZWNp3z7Az9MzwXsAcxKWWYCB/Wj3mZmNkBauVXl\nj4CfAn+bnn86pbUyFfW/AicCtXsbbAysKM1ftBCYUK+gpGOBYwEmTZrUwq7MulOnTuZ6LqB8tNIj\n+KuI+FFErEqPi4CxvRWS9DHgqYjo1/mFiJgREVMiYsrYsb3uzszM+qmVHsFySZ8GLk3PjwCebqHc\nrsDH0/2O16M4R/BdYGRpVtOJwKK+V9vMzAZKKz2CzwCHAUuAxcChwLTeCkXEyRExMSImA4cDN0bE\np4Cb0jYApgLX9KPeZmY2QFqZa2gB8PEB3OdJFFcjfQO4C7hgALdtZmZ91MrQ0BqLiJuBm9PyY8BO\n7divmZn1zpPOmZllzoHAzCxzvQYCSaeWlj0TqZlZl2k2DfVJknZh9RU+UMxEamZmXaTZyeIHKX5N\nvLmkW9PzjSVtHREPtaV2ZmZWuWZDQyuAU4BHgN0pfgwGxX2Mf19xvczMrE2a9Qj2Af4R2AI4B7gH\neDEiev0xmZmZDR0NewQRcUpE7AnMB34CrA2MlXSbpF+0qX5mZlaxVn5Qdl1EzAZmS/p8ROwmaUzV\nFTMzs/bo9fLRiDix9PTolLa8qgqZmVl79WmKiTW8U5lZNnxDeBtK/MtiM7PMtWXSOTNrrj93A3Ov\nwwaKewRmZplzIDAzy5wDgZlZ5hwIzMwy55PFZi3oz8lcs6HCPQIzs8w5EJiZZc6BwMwscw4EZmaZ\ncyAwM8ucA4GZWeZ8+ah1tcE2H48vQ7XByD0CM7PMuUdgVtINR+zd0AZrL/cIzMwy50BgZpY5BwIz\ns8w5EJiZZc6BwMwscw4EZmaZ8+Wj1tBg+zGWmVXDPQIzs8y5RzCE+YjdzAaCewRmZpmrLBBI2lTS\nTZIekHS/pC+m9NGSfivp4fR3VFV1MDOz3lXZI1gF/H1EbAPsDPydpG2A6cANEbEVcEN6bmZmHVJZ\nIIiIxRFxZ1peCcwDJgAHAjNTtpnAQVXVwczMeteWk8WSJgPbA7cD4yJicVq1BBjXoMyxwLEAkyZN\nqr6SNqR5xk2z/qv8ZLGkDYErgBMi4vnyuogIIOqVi4gZETElIqaMHTu26mqamWWr0h6BpHUogsAl\nEXFlSl4qaXxELJY0HniqyjrY4OZLYM06r8qrhgRcAMyLiHNKq64FpqblqcA1VdXBzMx6V2WPYFfg\nSOBeSXNT2inAWcDPJR0DLAAOq7AO1oLBOL5edU9hMLbZrFMqCwQRcRugBqv3rGq/ZmbWN/5lsZlZ\n5hwIzMwy50BgZpY5zz6akapPkLbjUlCf5DUbeO4RmJllzj2CDujrkbOPgs2sSu4RmJllzoHAzCxz\nHhqqSH+GczwEZGad4B6BmVnm3COwPnPPxay7uEdgZpY59wiscu5BmA1u7hGYmWXOgcDMLHMOBGZm\nmXMgMDPLnAOBmVnmHAjMzDLny0db1I659geKL9c0s75wj8DMLHPuEawhH32b2VDnHoGZWeYcCMzM\nMpft0NBQOvlrZlYl9wjMzDLnQGBmljkHAjOzzDkQmJllzoHAzCxzDgRmZplzIDAzy5wDgZlZ5rL9\nQVkjnjvIzHLjHoGZWea6vkfgI3wzs+bcIzAzy1xHAoGkj0p6SNIjkqZ3og5mZlZoeyCQtDZwLrAv\nsA1whKRt2l0PMzMrdKJHsBPwSEQ8FhGvAJcBB3agHmZmRmdOFk8Anig9Xwh8oGcmSccCx6anL0h6\nqI/7GQMs71cNhy63OR85tju7NuvsNW7zO1vJNGivGoqIGcCM/paXNDsipgxglQY9tzkfObbbba5O\nJ4aGFgGblp5PTGlmZtYBnQgE/w1sJWkzSesChwPXdqAeZmZGB4aGImKVpOOB64C1gQsj4v4KdtXv\nYaUhzG3OR47tdpsroohox37MzGyQ8i+Lzcwy50BgZpa5rgsE3TB9haT5ku6VNFfS7JQ2WtJvJT2c\n/o5K6ZL0vdTeeyTtUNrO1JT/YUlTS+k7pu0/ksqq/a0ESRdKekrSfaW0ytvZaB8dbPPpkhal93uu\npP1K605O9X9I0j6l9Lqf83QRxu0p/WfpggwkDU/PH0nrJ7enxSBpU0k3SXpA0v2SvpjSu/a9btLm\nwfleR0TXPChOPj8KbA6sC9wNbNPpevWjHfOBMT3S/hmYnpanA2en5f2A/wAE7AzcntJHA4+lv6PS\n8qi07o6UV6nsvh1q54eAHYD72tnORvvoYJtPB/6hTt5t0md4OLBZ+myv3exzDvwcODwtnw98Pi3/\nH+D8tHw48LM2tnk8sENaHgH8MbWta9/rJm0elO912//5K37xdwGuKz0/GTi50/XqRzvm89ZA8BAw\nvvQheygtfx84omc+4Ajg+6X076e08cCDpfQ35etAWyfz5i/FytvZaB8dbHOjL4c3fX4prrTbpdHn\nPH0JLgeGpfQ38tXKpuVhKZ869J5fA3wkh/e6TpsH5XvdbUND9aavmNChuqyJAP5T0hwVU20AjIuI\nxWl5CTAuLTdqc7P0hXXSB4t2tLPRPjrp+DQMcmFp+KKvbd4YWBERq3qkv2lbaf1zKX9bpWGK7YHb\nyeS97tFmGITvdbcFgm6xW0TsQDFD699J+lB5ZRShvuuv+21HOwfJa3kesAWwHbAY+JfOVqcakjYE\nrgBOiIjny+u69b2u0+ZB+V53WyDoiukrImJR+vsUcBXFjK1LJY0HSH+fStkbtblZ+sQ66YNFO9rZ\naB8dERFLI+K1iHgd+AHF+w19b/PTwEhJw3qkv2lbaf1GKX9bSFqH4gvxkoi4MiV39Xtdr82D9b3u\ntkAw5KevkLSBpBG1ZWBv4D6KdtSukphKMeZISj8qXWmxM/Bc6gpfB+wtaVTqfu5NMYa4GHhe0s7p\nyoqjStsaDNrRzkb76IjaF1VyMMX7DUU9D09XgWwGbEVxUrTu5zwd8d4EHJrK93z9am0+FLgx5a9c\nev0vAOZFxDmlVV37Xjdq86B9rztx4qTikzL7UZyhfxT4aqfr04/6b05xZcDdwP21NlCM8d0APAxc\nD4xO6aK40c+jwL3AlNK2PgM8kh7TSulT0gfwUeDf6dxJw0spusevUoxxHtOOdjbaRwfb/JPUpnvS\nP/H4Uv6vpvo/ROnqrkaf8/T5uSO9FpcDw1P6eun5I2n95m1s824UQzL3AHPTY79ufq+btHlQvtee\nYsLMLHPdNjRkZmZ95EBgZpY5BwIzs8w5EJiZZc6BwMwscw4E1rUkvZZmeLxP0uWS1h/g7R8t6d97\nybO7pL8uPT9O0lEDWQ+zNeVAYN3spYjYLiLeC7wCHNeBOuwOvBEIIuL8iPhxB+ph1pADgeXiVmBL\nAElfTr2E+ySdkNImS3pQ0iWS5kmaVetBqLg/xJi0PEXSzT03LumANPf7XZKulzQuTTZ2HPCl1DP5\noIr56P8hldlO0n+lCciu0ur5+G+WdLakOyT9UdIHq395LGcOBNb10nwr+wL3StoRmAZ8gGL++s9J\n2j5l3Rr4fxHxbuB5inndW3UbsHNEbA9cBpwYEfMp5on/TuqZ3NqjzI+BkyLifRS/Nj2ttG5YROwE\nnNAj3WzAORBYN3ubpLnAbOBxirlfdgOuiogXI+IF4EqgdsT9RET8Li1fnPK2aiJwnaR7ga8A72mW\nWdJGwMiIuCUlzaS4aU1NbWK2ORT3LzCrzLDes5gNWS9FxHblBDW/K2fP+VZqz1ex+qBpvQZl/w04\nJyKulbQ7xQ1I1sTL6e9r+P/UKuYegeXmVuAgSeun2V0PTmkAkyTtkpY/STHcA8Ud43ZMy4c02O5G\nrJ4GeGopfSXFrQrfJCKeA54tjf8fCdzSM59ZOzgQWFYi4k7gIopZGW8HfhgRd6XVD1HcCGgexT1x\nz0vpZwDflTSb4gi9ntOByyXNobg1YM0vgINrJ4t7lJkKfEvSPRQ3KjlzTdpm1l+efdSMN24n+Mt0\nqalZVtwjMDPLnHsEZmaZc4/AzCxzDgRmZplzIDAzy5wDgZlZ5hwIzMwy9/8BacD3vVNflwMAAAAA\nSUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fead3277f28>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(chunks_250000['2010 Census Population'], bins=50)\n",
    "plt.title(\"Population per Chunk (250,000)\")\n",
    "plt.xlabel(\"Population\")\n",
    "plt.ylabel(\"# of Chunks\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "chunks_250000.to_csv('geochunk_zip_250000.csv')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style>\n",
       "    .dataframe thead tr:only-child th {\n",
       "        text-align: right;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: left;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>2010 Census Population</th>\n",
       "      <th>Prefixes</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Chunk</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>018_0</th>\n",
       "      <td>411918</td>\n",
       "      <td>[0180, 0181, 0182, 0183, 0184]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>018_1</th>\n",
       "      <td>301208</td>\n",
       "      <td>[0185, 0186, 0187, 0188, 0189]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>01_0</th>\n",
       "      <td>469375</td>\n",
       "      <td>[010]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>01_1</th>\n",
       "      <td>382635</td>\n",
       "      <td>[011, 012, 013]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>01_2</th>\n",
       "      <td>211812</td>\n",
       "      <td>[014]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>01_3</th>\n",
       "      <td>373719</td>\n",
       "      <td>[015]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>01_4</th>\n",
       "      <td>188402</td>\n",
       "      <td>[016]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>01_5</th>\n",
       "      <td>379027</td>\n",
       "      <td>[017]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>01_6</th>\n",
       "      <td>476314</td>\n",
       "      <td>[019]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>021_0</th>\n",
       "      <td>396563</td>\n",
       "      <td>[0210, 0211, 0212]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>021_1</th>\n",
       "      <td>464077</td>\n",
       "      <td>[0213, 0214]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>021_2</th>\n",
       "      <td>423302</td>\n",
       "      <td>[0215, 0216, 0217, 0218, 0219]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>027_0</th>\n",
       "      <td>461337</td>\n",
       "      <td>[0270, 0271, 0272, 0273, 0274, 0275, 0276, 0277]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>027_1</th>\n",
       "      <td>64999</td>\n",
       "      <td>[0278, 0279]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>028_0</th>\n",
       "      <td>451970</td>\n",
       "      <td>[0280, 0281, 0282, 0283, 0284, 0285, 0286, 0287]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>028_1</th>\n",
       "      <td>221007</td>\n",
       "      <td>[0288, 0289]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>02_0</th>\n",
       "      <td>428586</td>\n",
       "      <td>[020, 022]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>02_1</th>\n",
       "      <td>441063</td>\n",
       "      <td>[023]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>02_2</th>\n",
       "      <td>408771</td>\n",
       "      <td>[024]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>02_3</th>\n",
       "      <td>264504</td>\n",
       "      <td>[025, 026]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>02_4</th>\n",
       "      <td>379604</td>\n",
       "      <td>[029]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>03_0</th>\n",
       "      <td>377547</td>\n",
       "      <td>[030]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>03_1</th>\n",
       "      <td>481457</td>\n",
       "      <td>[031, 032, 033, 034]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>03_2</th>\n",
       "      <td>457569</td>\n",
       "      <td>[035, 036, 037, 038]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>03_3</th>\n",
       "      <td>48241</td>\n",
       "      <td>[039]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>04_0</th>\n",
       "      <td>449276</td>\n",
       "      <td>[040, 041]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>04_1</th>\n",
       "      <td>484241</td>\n",
       "      <td>[042, 043, 044, 045]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>04_2</th>\n",
       "      <td>346497</td>\n",
       "      <td>[046, 047, 048, 049]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>05_0</th>\n",
       "      <td>480908</td>\n",
       "      <td>[050, 051, 052, 053, 054, 055, 056]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>05_1</th>\n",
       "      <td>144833</td>\n",
       "      <td>[057, 058, 059]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>972_0</th>\n",
       "      <td>392534</td>\n",
       "      <td>[9720, 9721]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>972_1</th>\n",
       "      <td>449177</td>\n",
       "      <td>[9722, 9723, 9724, 9725, 9726, 9727, 9728, 9729]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>973_0</th>\n",
       "      <td>462185</td>\n",
       "      <td>[9730, 9731, 9732, 9733, 9734, 9735]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>973_1</th>\n",
       "      <td>135381</td>\n",
       "      <td>[9736, 9737, 9738, 9739]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>974_0</th>\n",
       "      <td>394715</td>\n",
       "      <td>[9740, 9741, 9742, 9743, 9744, 9745, 9746]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>974_1</th>\n",
       "      <td>161028</td>\n",
       "      <td>[9747, 9748, 9749]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>97_0</th>\n",
       "      <td>319301</td>\n",
       "      <td>[971]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>97_1</th>\n",
       "      <td>355650</td>\n",
       "      <td>[975, 976]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>97_2</th>\n",
       "      <td>389914</td>\n",
       "      <td>[977, 978, 979]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>980_0</th>\n",
       "      <td>341896</td>\n",
       "      <td>[9800, 9801]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>980_1</th>\n",
       "      <td>226002</td>\n",
       "      <td>[9802]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>980_2</th>\n",
       "      <td>382521</td>\n",
       "      <td>[9803, 9804]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>980_3</th>\n",
       "      <td>398103</td>\n",
       "      <td>[9805, 9806, 9807, 9808, 9809]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>981_0</th>\n",
       "      <td>428717</td>\n",
       "      <td>[9810, 9811]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>981_1</th>\n",
       "      <td>432100</td>\n",
       "      <td>[9812, 9813, 9814, 9815, 9816, 9817, 9818, 9819]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>982_0</th>\n",
       "      <td>447430</td>\n",
       "      <td>[9820, 9821, 9822, 9823, 9824]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>982_1</th>\n",
       "      <td>430055</td>\n",
       "      <td>[9825, 9826, 9827, 9828, 9829]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>983_0</th>\n",
       "      <td>388948</td>\n",
       "      <td>[9830, 9831, 9832, 9833, 9834, 9835, 9836]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>983_1</th>\n",
       "      <td>324660</td>\n",
       "      <td>[9837, 9838, 9839]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>986_0</th>\n",
       "      <td>419035</td>\n",
       "      <td>[9860, 9861, 9862, 9863, 9864, 9865, 9866, 9867]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>986_1</th>\n",
       "      <td>154295</td>\n",
       "      <td>[9868, 9869]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>98_0</th>\n",
       "      <td>400115</td>\n",
       "      <td>[984]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>98_1</th>\n",
       "      <td>475554</td>\n",
       "      <td>[985, 987]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>98_2</th>\n",
       "      <td>219761</td>\n",
       "      <td>[988]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>98_3</th>\n",
       "      <td>284655</td>\n",
       "      <td>[989]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>99_0</th>\n",
       "      <td>257246</td>\n",
       "      <td>[990, 991]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>99_1</th>\n",
       "      <td>342856</td>\n",
       "      <td>[992]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>99_2</th>\n",
       "      <td>370680</td>\n",
       "      <td>[993, 994]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>99_3</th>\n",
       "      <td>320660</td>\n",
       "      <td>[995]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>99_4</th>\n",
       "      <td>389270</td>\n",
       "      <td>[996, 997, 998, 999]</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>875 rows × 2 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "       2010 Census Population  \\\n",
       "Chunk                           \n",
       "018_0                  411918   \n",
       "018_1                  301208   \n",
       "01_0                   469375   \n",
       "01_1                   382635   \n",
       "01_2                   211812   \n",
       "01_3                   373719   \n",
       "01_4                   188402   \n",
       "01_5                   379027   \n",
       "01_6                   476314   \n",
       "021_0                  396563   \n",
       "021_1                  464077   \n",
       "021_2                  423302   \n",
       "027_0                  461337   \n",
       "027_1                   64999   \n",
       "028_0                  451970   \n",
       "028_1                  221007   \n",
       "02_0                   428586   \n",
       "02_1                   441063   \n",
       "02_2                   408771   \n",
       "02_3                   264504   \n",
       "02_4                   379604   \n",
       "03_0                   377547   \n",
       "03_1                   481457   \n",
       "03_2                   457569   \n",
       "03_3                    48241   \n",
       "04_0                   449276   \n",
       "04_1                   484241   \n",
       "04_2                   346497   \n",
       "05_0                   480908   \n",
       "05_1                   144833   \n",
       "...                       ...   \n",
       "972_0                  392534   \n",
       "972_1                  449177   \n",
       "973_0                  462185   \n",
       "973_1                  135381   \n",
       "974_0                  394715   \n",
       "974_1                  161028   \n",
       "97_0                   319301   \n",
       "97_1                   355650   \n",
       "97_2                   389914   \n",
       "980_0                  341896   \n",
       "980_1                  226002   \n",
       "980_2                  382521   \n",
       "980_3                  398103   \n",
       "981_0                  428717   \n",
       "981_1                  432100   \n",
       "982_0                  447430   \n",
       "982_1                  430055   \n",
       "983_0                  388948   \n",
       "983_1                  324660   \n",
       "986_0                  419035   \n",
       "986_1                  154295   \n",
       "98_0                   400115   \n",
       "98_1                   475554   \n",
       "98_2                   219761   \n",
       "98_3                   284655   \n",
       "99_0                   257246   \n",
       "99_1                   342856   \n",
       "99_2                   370680   \n",
       "99_3                   320660   \n",
       "99_4                   389270   \n",
       "\n",
       "                                               Prefixes  \n",
       "Chunk                                                    \n",
       "018_0                    [0180, 0181, 0182, 0183, 0184]  \n",
       "018_1                    [0185, 0186, 0187, 0188, 0189]  \n",
       "01_0                                              [010]  \n",
       "01_1                                    [011, 012, 013]  \n",
       "01_2                                              [014]  \n",
       "01_3                                              [015]  \n",
       "01_4                                              [016]  \n",
       "01_5                                              [017]  \n",
       "01_6                                              [019]  \n",
       "021_0                                [0210, 0211, 0212]  \n",
       "021_1                                      [0213, 0214]  \n",
       "021_2                    [0215, 0216, 0217, 0218, 0219]  \n",
       "027_0  [0270, 0271, 0272, 0273, 0274, 0275, 0276, 0277]  \n",
       "027_1                                      [0278, 0279]  \n",
       "028_0  [0280, 0281, 0282, 0283, 0284, 0285, 0286, 0287]  \n",
       "028_1                                      [0288, 0289]  \n",
       "02_0                                         [020, 022]  \n",
       "02_1                                              [023]  \n",
       "02_2                                              [024]  \n",
       "02_3                                         [025, 026]  \n",
       "02_4                                              [029]  \n",
       "03_0                                              [030]  \n",
       "03_1                               [031, 032, 033, 034]  \n",
       "03_2                               [035, 036, 037, 038]  \n",
       "03_3                                              [039]  \n",
       "04_0                                         [040, 041]  \n",
       "04_1                               [042, 043, 044, 045]  \n",
       "04_2                               [046, 047, 048, 049]  \n",
       "05_0                [050, 051, 052, 053, 054, 055, 056]  \n",
       "05_1                                    [057, 058, 059]  \n",
       "...                                                 ...  \n",
       "972_0                                      [9720, 9721]  \n",
       "972_1  [9722, 9723, 9724, 9725, 9726, 9727, 9728, 9729]  \n",
       "973_0              [9730, 9731, 9732, 9733, 9734, 9735]  \n",
       "973_1                          [9736, 9737, 9738, 9739]  \n",
       "974_0        [9740, 9741, 9742, 9743, 9744, 9745, 9746]  \n",
       "974_1                                [9747, 9748, 9749]  \n",
       "97_0                                              [971]  \n",
       "97_1                                         [975, 976]  \n",
       "97_2                                    [977, 978, 979]  \n",
       "980_0                                      [9800, 9801]  \n",
       "980_1                                            [9802]  \n",
       "980_2                                      [9803, 9804]  \n",
       "980_3                    [9805, 9806, 9807, 9808, 9809]  \n",
       "981_0                                      [9810, 9811]  \n",
       "981_1  [9812, 9813, 9814, 9815, 9816, 9817, 9818, 9819]  \n",
       "982_0                    [9820, 9821, 9822, 9823, 9824]  \n",
       "982_1                    [9825, 9826, 9827, 9828, 9829]  \n",
       "983_0        [9830, 9831, 9832, 9833, 9834, 9835, 9836]  \n",
       "983_1                                [9837, 9838, 9839]  \n",
       "986_0  [9860, 9861, 9862, 9863, 9864, 9865, 9866, 9867]  \n",
       "986_1                                      [9868, 9869]  \n",
       "98_0                                              [984]  \n",
       "98_1                                         [985, 987]  \n",
       "98_2                                              [988]  \n",
       "98_3                                              [989]  \n",
       "99_0                                         [990, 991]  \n",
       "99_1                                              [992]  \n",
       "99_2                                         [993, 994]  \n",
       "99_3                                              [995]  \n",
       "99_4                               [996, 997, 998, 999]  \n",
       "\n",
       "[875 rows x 2 columns]"
      ]
     },
     "execution_count": 24,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "chunks_500000 = build_chunks(500000)\n",
    "chunks_500000"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.text.Text at 0x7feacf3a7eb8>"
      ]
     },
     "execution_count": 25,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEWCAYAAABrDZDcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAG8dJREFUeJzt3Xm4XFWZ7/HvjwBhSCCExBgIGBDabvRCwJgGwb5pUEAQ\nQaVpBiUgbS7aXuHiFQKNLdA8LVxvY2s3LcShiYKMAgloNzJfcIBOGANhCJgIMWQAAgGZEt77x16H\nbIqqOnWSs6vOqfX7PE89tffa07vq1Km31lq79lZEYGZm+Vqv0wGYmVlnORGYmWXOicDMLHNOBGZm\nmXMiMDPLnBOBmVnmnAhsrUmaLOnpddj+Aklf78+YBoJ1fV162fcxku7sw/o7SZotSVXE026Sxkia\nJ2lop2PpJk4EXULSAkmvSHpJ0hJJF0ka1um4etT7AIuI4yPiHzoV07qQNEnSLyStkPScpLslHdvp\nuOr4B+D/RvrBkKTbJL2a3icvSXq0vLKkIyUtlPSypGsljSwtGynpmrRsoaQjGx1UhXMlPZse55aT\nkaQJkuZI+mN6ntDKthGxBLgVmNpvr5A5EXSZgyJiGLAbMBE4vcPxDHqS1q9TtgdwC3A7sAOwJfBF\n4OPtja45SWOBvwSurVn05YgYlh7vK63/fuBC4HPAGOCPwL+VtjsfeD0tOwr4XtqmnqnAIcAuwM7A\nQcD/SMfZEJgJXAxsAcwAZqbyptsml9TM27qKCD+64AEsAD5amv8WcH2a3gqYBTwHzAe+UFrvDOAq\n4HJgJXAPsEtpeQA7lOYvAs5O05OBp0vLpgFPpP08DHwqlf8Z8CqwGngJWFG7rzT/hRTfcynerWri\nOB54HFhB8aGkBq9Fb3XaCvgZsAz4HfCVOtteDLwI/E2d/d8JnN/kbzEZeBr4KrAUWAwcW1p+W3m/\nwDHAna3Utc6630rxbF4njqOBm2rK3nbsmmX/CPy0NP9eig/+4cCmafpPSst/ApzTYF+/BqaW5o8D\nfpum9wUWlf9+wO+B/XvbNs2vT5Gk3tPp/7tuebhF0IUkbQMcANybii6j+GDaCjgU+EdJe5c2ORi4\nEhgJ/BS4VtIGa3HoJ4CPAJsDZwIXSxobEfMoPth+E8W30BF1Yt4b+CZwGDAWWJjiLvsE8CGKb4mH\nAfs1iaVunSStB1wH3A9sDewDnChpv5ptrwJGUHz7LMe5CbBHWt7Mu9PrsDXFB9n5krboZZuypnWV\ntJ6k76fl+0bEC3X28d+AR+uUf1PSckm/kjS5VP5+itcFgIh4gvThnx6rIuKx0vr3p23qedu+atZ9\nP/BApE/15IGa5Y22JSJWUXxh2KXBsa2PnAi6y7WSVlB8Q7yd4gN/G2BP4JSIeDUi7gN+QPFtscec\niLgqIt4AzgM2Anbv68Ej4sqI+ENEvBkRl1N8o53U4uZHAT+KiHsi4jXgVGAPSeNL65wTESsi4vcU\n/cQT3rmbtzSq04eA0RFxVkS8HhFPAt8HDi9t+5uIuDbV45Wa/W5B8X+zuJf6vAGcFRFvRMQvKFpC\n7+tlm7Jmdd0AuJQiyR0UEX9ssI8RFC2islOA7SkS1HTgOknvTcuGAbUJ5QWKFsEwihZSvWX11O7r\nBWBY6utvdpzetu2xMtXP+sE7+j9tUDskIm4qF0jaCnguIsofCAspxhB6PNUzERFvpjNeturrwSUd\nDZwEjE9Fw4BRLW6+FUUXTk8cL0l6luIDa0Eqfqa0/h/T/htpVKcAtkoJs8cQ4I5629bxPPAmRavl\nkSbrPZu+ubYab61mdd2B4tvwpIh4vZdY3/ZBHRF3lWZnSDqCovX4LxTJarOafWxG8aH7ZpNl9dTu\nazPgpYgISc2O03TbUtlwim4z6wduEXS/PwAjJZU/ELal6KPtsU3PROo6GZe2g+JDaJPSuu+udxBJ\n76H4Zv1lYMvU/TMX6PkW19tlbv8AvKe0v00pBmEXNdyiuUZ1egr4XUSMKD2GR8QBpW0bxpq+ff8G\n+MxaxgXwMi28pk3MA44F/kNSs1bGAxRdOs0Ea/5GD1HqbpG0PTAUeCw91pe0Y2nbXdI29bxtXzXr\nPgTsXPMNf+ea5Y227RnA34G3dx/ZOnAi6HIR8RTF4Ns3JW0kaWeKPuuLS6t9UNKn0z/YicBrwG/T\nsvuAIyUNkbQ/8N8bHGpTig+VZQDpVMoPlJYvAcaVzgypdSlwbDqtcCjFwOVdEbGgbzXutU53Aysl\nnSJp41SvD0j6UB/2fTJwjKSvSdoSQNIukmrHNBq5D/i0pE0k7UDx9+iTiLgUOA24qdS1U+tGYDdJ\nG6UYR0jaL70P1pd0FPAXwH+m9S8BDpL0kZSIzwKujoiVEfEycDVwlqRNJe1JMZbyk7Tv8ZKi1JX3\nY+AkSVunVulXKU4OgGLAejXwFUlDJX05ld/SwrZQdDcuiIiFfXrRrCEngjwcQdFd8wfgGuAbNV1I\nM4G/puhK+Bzw6dS3DnACxel7Kyj68WtPRQQgIh4G/oni2/ISioHKX5VWuYXiW90zkpbX2f4m4OsU\nZ/Mspjhj5fDa9fqgbp0iYjXFQOwEijOGllOMmWze6o4j4tfA3unxpKTnKPrbf9HiLr5NMQi7hOLU\nyUuar94wjhkUH9a31Iyl9CxfQvG6H5yKNgDOpkjWy4H/SdGd+Fha/yGKQf1LKM52Gg58qbTLLwEb\np2WXAl9M20DRAlvImhbchRSD8g9StAx/nspI3VmHUIxTrQA+n+J4vbdtk6OAC1p5jaw1enu3m+VG\n0hkUp4d+ttOx9JdurNPakrQTRbKZFBX+s0s6HVgWERf2uvK6HeddFCdC7BoRr1Z5rJx4sNisi6WW\nWl+6vdb2OGdXfYx0nKUUv0uxfuSuITOzzLlryMwsc24RmJllblCMEYwaNSrGjx/f6TDMzAaVOXPm\nLI+I0b2tNygSwfjx45k9e3anwzAzG1QktfRbC3cNmZllzonAzCxzTgRmZplzIjAzy5wTgZlZ5pwI\nzMwy50RgZpY5JwIzs8w5EZiZZW5Q/LLYzKwbjJ/287rlC845sM2RvJ1bBGZmmXMiMDPLnBOBmVnm\nnAjMzDJX6WCxpAXASmA1sCoiJkoaCVwOjAcWAIdFxPNVxmFmZo21o0XwlxExISImpvlpwM0RsSNw\nc5o3M7MO6UTX0MHAjDQ9AzikAzGYmVlSdSII4JeS5kiamsrGRMTiNP0MMKbehpKmSpotafayZcsq\nDtPMLF9V/6Bsr4hYJOldwI2SHikvjIiQFPU2jIjpwHSAiRMn1l3HzMzWXaUtgohYlJ6XAtcAk4Al\nksYCpOelVcZgZmbNVZYIJG0qaXjPNLAvMBeYBUxJq00BZlYVg5mZ9a7KrqExwDWSeo7z04j4T0n/\nBVwh6ThgIXBYhTGYmVkvKksEEfEksEud8meBfao6rpmZ9Y1/WWxmljknAjOzzPl+BGZmHdbp+xS4\nRWBmljknAjOzzDkRmJllzonAzCxzTgRmZplzIjAzy5wTgZlZ5pwIzMwy50RgZpY5JwIzs8w5EZiZ\nZc6JwMwsc04EZmaZcyIwM8ucE4GZWeacCMzMMudEYGaWOScCM7PMORGYmWXO9yw2M+tnje5BPFC5\nRWBmljknAjOzzDkRmJllzonAzCxzTgRmZplzIjAzy5wTgZlZ5ipPBJKGSLpX0vVpfjtJd0maL+ly\nSRtWHYOZmTXWjhbBCcC80vy5wLcjYgfgeeC4NsRgZmYNVJoIJI0DDgR+kOYF7A1clVaZARxSZQxm\nZtZc1S2CfwZOBt5M81sCKyJiVZp/Gti64hjMzKyJyhKBpE8ASyNizlpuP1XSbEmzly1b1s/RmZlZ\njypbBHsCn5S0ALiMokvoO8AIST0XuxsHLKq3cURMj4iJETFx9OjRFYZpZpa3yhJBRJwaEeMiYjxw\nOHBLRBwF3AocmlabAsysKgYzM+tdJ35HcApwkqT5FGMGP+xADGZmlrTlfgQRcRtwW5p+EpjUjuOa\nmVnv/MtiM7PMORGYmWXOt6o0M+vFYLv1ZF+5RWBmljknAjOzzDkRmJllzmMEZta1GvXtLzjnwD6t\n3+3cIjAzy5wTgZlZ5pwIzMwy5zECMxtQmvXTN+rbt3XjFoGZWeacCMzMMudEYGaWOScCM7PMORGY\nmWXOicDMLHO9JgJJfyVpeJo+XdLVknarPjQzM2uHVloEX4+IlZL2Aj5KcY/h71UblpmZtUsriWB1\nej4QmB4RPwc2rC4kMzNrp1YSwSJJFwJ/DfxC0tAWtzMzs0GglQ/0w4AbgP0iYgUwEvhapVGZmVnb\ntHKtoSMi4oc9MxGxWNIJwC+rC8vMrDq53negkVYSwWckvRoRlwBIOh/YqNqwzMysXVpKBMAsSW8C\n+wMrIuK4asMyM7N2aZgIJI0szf4NcC3wK+BMSSMj4rmqgzMzs+o1axHMAQJQ6fnA9Ahg+8qjMzMr\n6es9iK01DRNBRGzXzkDMzKwzWrpDmaQPA+PL60fEjyuKyczM2qjXRCDpJ8B7gftY8yvjAJwIzMy6\nQCstgonAThERVQdjZmbt18ovi+cC7+7rjiVtJOluSfdLekjSmal8O0l3SZov6XJJvm6RmVkHtZII\nRgEPS7pB0qyeRwvbvQbsHRG7ABOA/SXtDpwLfDsidgCeB/ybBDOzDmqla+iMtdlx6kp6Kc1ukB4B\n7A0cmcpnpP37stZmZh3SayKIiNvXdueShlD8HmEH4HzgCYpfJq9KqzwNbN1g26nAVIBtt912bUMw\nsw7zuf8DXyt3KFsp6cX0eFXSakkvtrLziFgdEROAccAk4E9bDSwipkfExIiYOHr06FY3MzOzPmql\nRTC8Z1qSgIOB3ftykIhYIelWYA9ghKT1U6tgHLCobyGbmVl/6tMNZqJwLbBfb+tKGi1pRJreGPgY\nMA+4FTg0rTYFmNmniM3MrF+18oOyT5dm16P4XcGrLex7LDAjjROsB1wREddLehi4TNLZwL0U90A2\nM7MOaeWsoYNK06uABRTdQ01FxAPArnXKn6QYLzAzswGglTGCY9sRiJmZdUYrXUOjgS/wzovOfb66\nsMzMrF1a6RqaCdwB3MSai86ZmVmXaCURbBIRp1QeiZmZdUQrp49eL+mAyiMxM7OOaHbP4pWsuUXl\naZJeA95I8xERm7UnRDMzq1KzW1UOb7TMzKxWo2sKdfuxu0HDriFJ+0k6tE75ZyR9rNqwzMysXZqN\nEfw9UO/Ko7cDZ1UTjpmZtVuzRDA0IpbVFkbEcmDT6kIyM7N2apYINpP0jjEESRsAG1cXkpmZtVOz\nRHA18H1Jb337lzQMuCAtMzOzLtAsEZwOLAEWSpojaQ7wO2BZWmZmZl2g2emjq4Bpks6kuNUkwPyI\neKUtkZmZWVu0cvXRV4AH2xCLmdXhe/5a1fp0hzIzM+s+zX5Qtmd6Htq+cMzMrN2atQi+m55/045A\nzMysM5qNEbwhaTqwtaTv1i6MiK9UF5aZmbVLs0TwCeCjwH7AnPaEY2Zm7dbs9NHlwGWS5kXE/W2M\nyczM2qiVs4aelXSNpKXp8TNJ4yqPzMzM2qKVRPDvwCxgq/S4LpWZmVkXaCURvCsi/j0iVqXHRcDo\niuMyM7M2aSURLJf0WUlD0uOzwLNVB2ZmZu3RSiL4PHAY8AywGDgUOLbKoMzMrH1audbQQuCTbYjF\nzMw6wNcaMjPLnBOBmVnmnAjMzDLX6xiBpNMj4uw0PTQiXqs+LDNbW4Pl/gWN4rT2a3YZ6lMk7UFx\nllCPlq9EKmkbSbdKeljSQ5JOSOUjJd0o6fH0vMXah29mZuuqWdfQI8BfAdtLukPS94EtJb2vxX2v\nAr4aETsBuwN/K2knYBpwc0TsCNyc5s3MrEOaJYIVwGnAfGAy8J1UPk3Sr3vbcUQsjoh70vRKYB6w\nNXAwMCOtNgM4ZK0iNzOzftEsEewH/Bx4L3Ae8OfAyxFxbER8uC8HkTQe2BW4CxgTEYvTomeAMQ22\nmSpptqTZy5Yt68vhzMysDxomgog4LSL2ARYAPwGGAKMl3SnpulYPIGkY8DPgxIh4seYYAUSD40+P\niIkRMXH0aF/ayMysKr2eNQTcEBGzgdmSvhgRe0ka1crOJW1AkQQuiYirU/ESSWMjYrGkscDStQvd\nzMz6Q6+/I4iIk0uzx6Sy5b1tJ0nAD4F5EXFeadEsYEqangLMbDVYMzPrf620CN7SxzuV7Ql8DnhQ\n0n2p7DTgHOAKSccBCykuaGdmA8xg+T2Crbs+JYK+iIg7ATVYvE9VxzUzs77xJSbMzDLnRGBmlrnK\nuobMrHW+7o51klsEZmaZcyIwM8ucE4GZWeY8RmBZ8jnya8/jGd3HLQIzs8w5EZiZZc6JwMwscx4j\nMMuEx0WsEbcIzMwy50RgZpY5JwIzs8x5jMBsHficeusGbhGYmWXOicDMLHNOBGZmmfMYgVnmPM5h\nbhGYmWXOicDMLHNOBGZmmfMYgVkL3I9u3cwtAjOzzDkRmJllzonAzCxzHiMwG6Q8bmH9xS0CM7PM\nORGYmWXOicDMLHMeI7CO6q/76Lq/3GztVdYikPQjSUslzS2VjZR0o6TH0/MWVR3fzMxaU2XX0EXA\n/jVl04CbI2JH4OY0b2ZmHVRZIoiI/wc8V1N8MDAjTc8ADqnq+GZm1pp2jxGMiYjFafoZYEyjFSVN\nBaYCbLvttm0IbfDpr/71/jLQ4jGz1nTsrKGICCCaLJ8eERMjYuLo0aPbGJmZWV7anQiWSBoLkJ6X\ntvn4ZmZWo92JYBYwJU1PAWa2+fhmZlajsjECSZcCk4FRkp4GvgGcA1wh6ThgIXBYVccfyKruS+/r\n/gdi336nYvLvESxHlSWCiDiiwaJ9qjqmmZn1nS8xYWaWOScCM7PM+VpDFXFfs5kNFm4RmJllzonA\nzCxzTgRmZplzIjAzy5wTgZlZ5pwIzMwy50RgZpY5/44gI339bcNA/C3EQIzJbLBzi8DMLHNOBGZm\nmXMiMDPLnMcIWuS+6bXn185sYHOLwMwsc04EZmaZcyIwM8ucxwhqDMT+7IEYk5l1D7cIzMwy50Rg\nZpY5JwIzs8xlO0bgfnczs4JbBGZmmXMiMDPLnBOBmVnmun6MwGMBZmbNuUVgZpY5JwIzs8w5EZiZ\nZa7rxwgGE49nmFkndKRFIGl/SY9Kmi9pWidiMDOzQtsTgaQhwPnAx4GdgCMk7dTuOMzMrNCJFsEk\nYH5EPBkRrwOXAQd3IA4zM6MzYwRbA0+V5p8G/rx2JUlTgalp9iVJj5YWjwKWVxbhwOa658l1z5DO\nXee6v6eVlQbsYHFETAem11smaXZETGxzSAOC6+6658Z1r77unegaWgRsU5ofl8rMzKwDOpEI/gvY\nUdJ2kjYEDgdmdSAOMzOjA11DEbFK0peBG4AhwI8i4qE+7qZul1EmXPc8ue55akvdFRHtOI6ZmQ1Q\nvsSEmVnmnAjMzDI3qBLBYL40haQfSVoqaW6pbKSkGyU9np63SOWS9N1Uzwck7VbaZkpa/3FJU0rl\nH5T0YNrmu5LU7BjtJGkbSbdKeljSQ5JOyKX+kjaSdLek+1Pdz0zl20m6K8V7eTpxAklD0/z8tHx8\naV+npvJHJe1XKq/7f9HoGO0maYikeyVd3yyubqu7pAXpPXmfpNmpbGC+5yNiUDwoBpafALYHNgTu\nB3bqdFx9iP8vgN2AuaWy/wNMS9PTgHPT9AHAfwACdgfuSuUjgSfT8xZpeou07O60rtK2H292jDbX\nfSywW5oeDjxGcXmRrq9/imdYmt4AuCvFeQVweCq/APhimv4ScEGaPhy4PE3vlN7zQ4Ht0v/CkGb/\nF42O0YG//0nAT4Hrm8XVbXUHFgCjasoG5Hu+7W+KdXhR9wBuKM2fCpza6bj6WIfxvD0RPAqMTdNj\ngUfT9IXAEbXrAUcAF5bKL0xlY4FHSuVvrdfoGB1+HWYCH8ut/sAmwD0Uv6RfDqyfyt96b1OcTbdH\nml4/rafa93vPeo3+L9I2dY/R5jqPA24G9gaubxZXF9Z9Ae9MBAPyPT+YuobqXZpi6w7F0l/GRMTi\nNP0MMCZNN6prs/Kn65Q3O0ZHpOb+rhTfjLOof+oauQ9YCtxI8S12RUSsSquU432rjmn5C8CW9P01\n2bLJMdrpn4GTgTfTfLO4uq3uAfxS0hwVl8yBAfqeH7CXmMhNRISkSs/lbccxmpE0DPgZcGJEvJi6\nNNsWW6fqHxGrgQmSRgDXAH/a7hg6QdIngKURMUfS5E7H0wF7RcQiSe8CbpT0SHnhQHrPD6YWQTde\nmmKJpLEA6XlpKm9U12bl4+qUNztGW0nagCIJXBIRV/cSW9fVHyAiVgC3UnRVjJDU80WsHO9bdUzL\nNweepe+vybNNjtEuewKflLSA4irDewPfaRJXN9WdiFiUnpdSfAGYxAB9zw+mRNCNl6aYBfScBTCF\nou+8p/zodCbB7sALqal3A7CvpC3SmQD7UvR9LgZelLR7OnPg6Jp91TtG26SYfgjMi4jzSou6vv6S\nRqeWAJI2phgbmUeREA6tE1c53kOBW6Lo7J0FHJ7OrNkO2JFisLDu/0XaptEx2iIiTo2IcRExPsV1\nS0Qc1SSurqm7pE0lDe+ZpnivzmWgvufbPYCyjoMvB1CccfIE8HedjqePsV8KLAbeoOjPO46iL/Nm\n4HHgJmBkWlcUN+95AngQmFjaz+eB+elxbKl8YnqjPQH8K2t+NV73GG2u+14U/aUPAPelxwE51B/Y\nGbg31X0u8PepfHuKD7P5wJXA0FS+UZqfn5ZvX9rX36X6PUo6Q6TZ/0WjY3To/T+ZNWcNdX3d0/Hv\nT4+HemIbqO95X2LCzCxzg6lryMzMKuBEYGaWOScCM7PMORGYmWXOicDMLHNOBNa1JK1OV36cK+lK\nSZv08/6PkfSvvawzWdKHS/PHSzq6P+MwW1dOBNbNXomICRHxAeB14PgOxDAZeCsRRMQFEfHjDsRh\n1pATgeXiDmAHAEknpVbCXEknprLxkh6RdImkeZKu6mlBqLiu/Kg0PVHSbbU7l3SQiuvf3yvpJklj\n0gX2jgf+V2qZfETSGZL+d9pmgqTfqrj+/DVac2362ySdq+I+Bo9J+kj1L4/lzInAul665szHgQcl\nfRA4luJS0LsDX5C0a1r1fcC/RcSfAS9SXB+/VXcCu0fErhTX1Tk5IhZQXAv/26llckfNNj8GTomI\nnSl+TfqN0rL1I2IScGJNuVm/cyKwbraxiss/zwZ+T3G9o72AayLi5Yh4Cbga6PnG/VRE/CpNX5zW\nbdU44AZJDwJfA97fbGVJmwMjIuL2VDSD4uZFPXouzDeH4j4WZpXxZaitm70SERPKBSpd+rqO2uut\n9MyvYs2Xpo0abPsvwHkRMStdcvmMPkX6Tq+l59X4/9Qq5haB5eYO4BBJm6SrQn4qlQFsK2mPNH0k\nRXcPFHea+mCa/kyD/W7OmssATymVr6S4PefbRMQLwPOl/v/PAbfXrmfWDk4ElpWIuAe4iOLKlHcB\nP4iIe9PiR4G/lTSP4v6w30vlZwLfUXED8tUNdn0GcKWkORS3SexxHfCpnsHimm2mAN+S9AAwAThr\nXepmtrZ89VEz3rqF5vXpVFOzrLhFYGaWObcIzMwy5xaBmVnmnAjMzDLnRGBmljknAjOzzDkRmJll\n7v8DUCHvnUJ7nGUAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7feacf1ea828>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(chunks_500000['2010 Census Population'], bins=50)\n",
    "plt.title(\"Population per Chunk (500,000)\")\n",
    "plt.xlabel(\"Population\")\n",
    "plt.ylabel(\"# of Chunks\")"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.5.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}