2-generate_pairs.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import numpy as np\n",
    "import glob,os\n",
    "import matplotlib.pylab as plt\n",
    "\n",
    "\n",
    "from scipy.spatial.distance import pdist,squareform\n",
    "\n",
    "import cooler\n",
    "import gzip\n",
    "from io import StringIO"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "def plot_matrix(mat,title,cmap,vmax,vmin=0):\n",
    "    plt.figure()\n",
    "    plt.imshow(mat,interpolation='nearest',cmap=cmap,vmax=vmax,vmin=vmin)\n",
    "    plt.colorbar()\n",
    "    plt.title(title)\n",
    "    plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "mid_positions_cell = [500000 + 1000000 * i for i in range(243)]\n",
    "mid_positions_science = [34600000 + 30000 * i for i in range(83)]\n",
    "mid_positions_cell_250kb = pd.read_csv(\"pos1.txt\",sep=\"\\t\",header=None).iloc[:,0].tolist()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "header_lines_cell= [\n",
    "        \"## pairs format v1.0\",\n",
    "        \"#sorted: chr1-chr2-pos1-pos2\",\n",
    "        \"#shape: upper triangle\",\n",
    "        \"#chromosome: chr2 243199373\",\n",
    "        \"#columns: readID chr1 pos1 chr2 pos2 strand1 strand2 phase0 phase1 phase_prob00 phase_prob01 phase_prob10 phase_prob11\"\n",
    "    ]\n",
    "header_lines_science = [\n",
    "        \"## pairs format v1.0\",\n",
    "        \"#sorted: chr1-chr2-pos1-pos2\",\n",
    "        \"#shape: upper triangle\",\n",
    "        \"#chromosome: chr21 48129895\",\n",
    "        \"#columns: readID chr1 pos1 chr2 pos2 strand1 strand2 phase0 phase1 phase_prob00 phase_prob01 phase_prob10 phase_prob11\"\n",
    "    ]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "group1\n",
      "group2\n",
      "group3\n",
      "group4\n"
     ]
    }
   ],
   "source": [
    "for group in range(1,5):\n",
    "    print(f\"group{group}\")\n",
    "    matrix_all = np.load(f\"/shareb/mliu/evaluate_impute/data/simulation_hic/Science2018/hic/matrix/Science2018_sample_matrix{group}.npy\")\n",
    "    for idx, matrix in enumerate(matrix_all):\n",
    "        copy = str(idx)\n",
    "        filename = f'/shareb/mliu/evaluate_impute/data/simulation_hic/Science2018/hic/pairs/Science2018_group{group}_copy{copy}.pairs.gz'\n",
    "\n",
    "        idx1, idx2 = np.where(np.triu(matrix) > 0)\n",
    "        values = matrix[idx1, idx2]\n",
    "        idx1_repeats = [i for i, v in zip(idx1, values) for _ in range(int(v))]\n",
    "        idx2_repeats = [i for i, v in zip(idx2, values) for _ in range(int(v))]\n",
    "\n",
    "        pos1 = [mid_positions_science[i] for i in idx1_repeats]  \n",
    "        pos2 = [mid_positions_science[j] for j in idx2_repeats]\n",
    "\n",
    "\n",
    "        pairs = pd.DataFrame({\n",
    "        'chr1': ['chr21'] * len(pos1),\n",
    "        'pos1': pos1,\n",
    "        'chr2': ['chr21'] * len(pos2),\n",
    "        'pos2': pos2,\n",
    "            })\n",
    "\n",
    "        # pairs['readID'] = \".\"\n",
    "        # pairs['strand1'] = \"+\"\n",
    "        # pairs['strand2'] = \"+\"\n",
    "        # pairs['phase0'] = 0\n",
    "        # pairs['phase1'] = 0\n",
    "        # pairs['p00'] = 1.000\n",
    "        # pairs['p01'] = 0.000\n",
    "        # pairs['p10'] = 0.000\n",
    "        # pairs['p11'] = 0.000\n",
    "        # pairs = pairs[['readID','chr1','pos1',\"chr2\",\"pos2\",\"strand1\",\"strand2\",\"phase0\",\"phase1\",\"p00\",\"p01\",\"p10\",\"p11\"]]\n",
    "        pairs = pairs[['chr1','pos1',\"chr2\",\"pos2\"]]\n",
    "        \n",
    "        # 指定文件名\n",
    "    \n",
    "        csv_buffer = StringIO()\n",
    "        for line in header_lines_science:\n",
    "            csv_buffer.write(line + '\\n')\n",
    "        pairs.to_csv(csv_buffer, index=False, header=False,sep=\"\\t\")\n",
    "\n",
    "        # 使用gzip压缩数据并写入文件\n",
    "        with gzip.open(filename, 'wt') as f:\n",
    "            f.write(csv_buffer.getvalue())\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Cell"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "idx is 0\n",
      "4610\n",
      "idx is 1\n",
      "3320\n",
      "idx is 2\n",
      "4040\n",
      "idx is 3\n",
      "4748\n",
      "idx is 4\n",
      "4805\n",
      "idx is 5\n",
      "4571\n",
      "idx is 6\n",
      "4723\n",
      "idx is 7\n",
      "3911\n",
      "idx is 8\n",
      "3554\n",
      "idx is 9\n",
      "3059\n",
      "idx is 10\n",
      "3013\n",
      "idx is 11\n",
      "4442\n",
      "idx is 12\n",
      "3224\n",
      "idx is 13\n",
      "4326\n",
      "idx is 14\n",
      "2715\n",
      "idx is 15\n",
      "2734\n",
      "idx is 16\n",
      "3299\n",
      "idx is 17\n",
      "3374\n",
      "idx is 18\n",
      "2851\n",
      "idx is 19\n",
      "4751\n",
      "idx is 20\n",
      "4024\n",
      "idx is 21\n",
      "2775\n",
      "idx is 22\n",
      "3362\n",
      "idx is 23\n",
      "4053\n",
      "idx is 24\n",
      "3102\n",
      "idx is 25\n",
      "2804\n",
      "idx is 26\n",
      "2529\n",
      "idx is 27\n",
      "3148\n",
      "idx is 28\n",
      "3815\n",
      "idx is 29\n",
      "4589\n",
      "idx is 30\n",
      "4793\n",
      "idx is 31\n",
      "4824\n",
      "idx is 32\n",
      "2800\n",
      "idx is 33\n",
      "2756\n",
      "idx is 34\n",
      "2718\n",
      "idx is 35\n",
      "3926\n",
      "idx is 36\n",
      "3635\n",
      "idx is 37\n",
      "4267\n",
      "idx is 38\n",
      "3173\n",
      "idx is 39\n",
      "3592\n",
      "idx is 40\n",
      "4441\n",
      "idx is 41\n",
      "3052\n",
      "idx is 42\n",
      "3357\n",
      "idx is 43\n",
      "2840\n",
      "idx is 44\n",
      "2792\n",
      "idx is 45\n",
      "3620\n",
      "idx is 46\n",
      "3002\n",
      "idx is 47\n",
      "4791\n",
      "idx is 48\n",
      "4709\n",
      "idx is 49\n",
      "4386\n",
      "idx is 50\n",
      "4251\n",
      "idx is 51\n",
      "4024\n",
      "idx is 52\n",
      "2661\n",
      "idx is 53\n",
      "3921\n",
      "idx is 54\n",
      "4278\n",
      "idx is 55\n",
      "4695\n",
      "idx is 56\n",
      "2684\n",
      "idx is 57\n",
      "4364\n",
      "idx is 58\n",
      "4622\n",
      "idx is 59\n",
      "2678\n",
      "idx is 60\n",
      "3770\n",
      "idx is 61\n",
      "4258\n",
      "idx is 62\n",
      "4717\n",
      "idx is 63\n",
      "2526\n",
      "idx is 64\n",
      "4612\n",
      "idx is 65\n",
      "3810\n",
      "idx is 66\n",
      "4969\n",
      "idx is 67\n",
      "3465\n",
      "idx is 68\n",
      "3198\n",
      "idx is 69\n",
      "2980\n",
      "idx is 70\n",
      "3074\n",
      "idx is 71\n",
      "3482\n",
      "idx is 72\n",
      "3541\n",
      "idx is 73\n",
      "4676\n",
      "idx is 74\n",
      "3609\n",
      "idx is 75\n",
      "2689\n",
      "idx is 76\n",
      "4652\n",
      "idx is 77\n",
      "2998\n",
      "idx is 78\n",
      "3758\n",
      "idx is 79\n",
      "4353\n",
      "idx is 80\n",
      "3760\n",
      "idx is 81\n",
      "3400\n",
      "idx is 82\n",
      "2570\n",
      "idx is 83\n",
      "3996\n",
      "idx is 84\n",
      "3442\n",
      "idx is 85\n",
      "3991\n",
      "idx is 86\n",
      "4787\n",
      "idx is 87\n",
      "3391\n",
      "idx is 88\n",
      "4227\n",
      "idx is 89\n",
      "4386\n",
      "idx is 90\n",
      "3522\n",
      "idx is 91\n",
      "3163\n",
      "idx is 92\n",
      "3159\n",
      "idx is 93\n",
      "2654\n",
      "idx is 94\n",
      "3625\n",
      "idx is 95\n",
      "3553\n",
      "idx is 96\n",
      "2659\n",
      "idx is 97\n",
      "3588\n",
      "idx is 98\n",
      "2772\n",
      "idx is 99\n",
      "3419\n",
      "idx is 100\n",
      "3235\n",
      "idx is 101\n",
      "4008\n",
      "idx is 102\n",
      "2713\n",
      "idx is 103\n",
      "4581\n",
      "idx is 104\n",
      "4194\n",
      "idx is 105\n",
      "2984\n",
      "idx is 106\n",
      "2512\n",
      "idx is 107\n",
      "3698\n",
      "idx is 108\n",
      "2709\n",
      "idx is 109\n",
      "4598\n",
      "idx is 110\n",
      "2600\n",
      "idx is 111\n",
      "4641\n",
      "idx is 112\n",
      "4843\n",
      "idx is 113\n",
      "4661\n",
      "idx is 114\n",
      "3027\n",
      "idx is 115\n",
      "4126\n",
      "idx is 116\n",
      "2840\n",
      "idx is 117\n",
      "4233\n",
      "idx is 118\n",
      "3583\n",
      "idx is 119\n",
      "3994\n",
      "idx is 120\n",
      "2986\n",
      "idx is 121\n",
      "3876\n",
      "idx is 122\n",
      "3489\n",
      "idx is 123\n",
      "3294\n",
      "idx is 124\n",
      "3199\n",
      "idx is 125\n",
      "2527\n",
      "idx is 126\n",
      "4366\n",
      "idx is 127\n",
      "4553\n",
      "idx is 128\n",
      "2965\n",
      "idx is 129\n",
      "4591\n",
      "idx is 130\n",
      "4096\n",
      "idx is 131\n",
      "4473\n",
      "idx is 132\n",
      "2835\n",
      "idx is 133\n",
      "3269\n",
      "idx is 134\n",
      "4425\n",
      "idx is 135\n",
      "3649\n",
      "idx is 136\n",
      "4737\n",
      "idx is 137\n",
      "3664\n",
      "idx is 138\n",
      "3066\n",
      "idx is 139\n",
      "2536\n",
      "idx is 140\n",
      "4980\n",
      "idx is 141\n",
      "4076\n",
      "idx is 142\n",
      "4311\n",
      "idx is 143\n",
      "3386\n",
      "idx is 144\n",
      "3625\n",
      "idx is 145\n",
      "4985\n",
      "idx is 146\n",
      "3829\n",
      "idx is 147\n",
      "2830\n",
      "idx is 148\n",
      "3087\n",
      "idx is 149\n",
      "2843\n",
      "idx is 150\n",
      "3735\n",
      "idx is 151\n",
      "4673\n",
      "idx is 152\n",
      "3801\n",
      "idx is 153\n",
      "4509\n",
      "idx is 154\n",
      "3654\n",
      "idx is 155\n",
      "4499\n",
      "idx is 156\n",
      "3411\n",
      "idx is 157\n",
      "2850\n",
      "idx is 158\n",
      "3472\n",
      "idx is 159\n",
      "3944\n",
      "idx is 160\n",
      "4159\n",
      "idx is 161\n",
      "4958\n",
      "idx is 162\n",
      "4192\n",
      "idx is 163\n",
      "4876\n",
      "idx is 164\n",
      "4922\n",
      "idx is 165\n",
      "4646\n",
      "idx is 166\n",
      "4842\n",
      "idx is 167\n",
      "4317\n",
      "idx is 168\n",
      "4051\n",
      "idx is 169\n",
      "3335\n",
      "idx is 170\n",
      "3303\n",
      "idx is 171\n",
      "3257\n",
      "idx is 172\n",
      "3291\n",
      "idx is 173\n",
      "3419\n",
      "idx is 174\n",
      "4542\n",
      "idx is 175\n",
      "4220\n",
      "idx is 176\n",
      "4168\n",
      "idx is 177\n",
      "3027\n",
      "idx is 178\n",
      "3088\n",
      "idx is 179\n",
      "4704\n",
      "idx is 180\n",
      "3610\n",
      "idx is 181\n",
      "2564\n",
      "idx is 182\n",
      "2772\n",
      "idx is 183\n",
      "2502\n",
      "idx is 184\n",
      "4623\n",
      "idx is 185\n",
      "3644\n",
      "idx is 186\n",
      "3051\n",
      "idx is 187\n",
      "2686\n",
      "idx is 188\n",
      "4920\n",
      "idx is 189\n",
      "4584\n",
      "idx is 190\n",
      "4753\n",
      "idx is 191\n",
      "2740\n",
      "idx is 192\n",
      "3022\n",
      "idx is 193\n",
      "3584\n",
      "idx is 194\n",
      "4708\n",
      "idx is 195\n",
      "4955\n",
      "idx is 196\n",
      "2883\n",
      "idx is 197\n",
      "4457\n",
      "idx is 198\n",
      "4076\n",
      "idx is 199\n",
      "3695\n",
      "idx is 200\n",
      "4316\n",
      "idx is 201\n",
      "3067\n",
      "idx is 202\n",
      "3209\n",
      "idx is 203\n",
      "4934\n",
      "idx is 204\n",
      "3989\n",
      "idx is 205\n",
      "2900\n",
      "idx is 206\n",
      "4486\n",
      "idx is 207\n",
      "2975\n",
      "idx is 208\n",
      "4658\n",
      "idx is 209\n",
      "3401\n",
      "idx is 210\n",
      "4497\n",
      "idx is 211\n",
      "2726\n",
      "idx is 212\n",
      "3597\n",
      "idx is 213\n",
      "2733\n",
      "idx is 214\n",
      "4478\n",
      "idx is 215\n",
      "4233\n",
      "idx is 216\n",
      "2784\n",
      "idx is 217\n",
      "4380\n",
      "idx is 218\n",
      "4698\n",
      "idx is 219\n",
      "3685\n",
      "idx is 220\n",
      "2738\n",
      "idx is 221\n",
      "3538\n",
      "idx is 222\n",
      "3413\n",
      "idx is 223\n",
      "4581\n",
      "idx is 224\n",
      "4663\n",
      "idx is 225\n",
      "3582\n",
      "idx is 226\n",
      "3889\n",
      "idx is 227\n",
      "4979\n",
      "idx is 228\n",
      "4004\n",
      "idx is 229\n",
      "4083\n",
      "idx is 230\n",
      "2547\n",
      "idx is 231\n",
      "3128\n",
      "idx is 232\n",
      "3485\n",
      "idx is 233\n",
      "3831\n",
      "idx is 234\n",
      "3849\n",
      "idx is 235\n",
      "3146\n",
      "idx is 236\n",
      "2832\n",
      "idx is 237\n",
      "3864\n",
      "idx is 238\n",
      "3167\n",
      "idx is 239\n",
      "3504\n",
      "idx is 240\n",
      "3096\n",
      "idx is 241\n",
      "2603\n",
      "idx is 242\n",
      "3810\n",
      "idx is 243\n",
      "4788\n",
      "idx is 244\n",
      "4398\n",
      "idx is 245\n",
      "3488\n",
      "idx is 246\n",
      "3569\n",
      "idx is 247\n",
      "4993\n",
      "idx is 248\n",
      "2963\n",
      "idx is 249\n",
      "3019\n",
      "idx is 250\n",
      "4358\n",
      "idx is 251\n",
      "2946\n",
      "idx is 252\n",
      "3519\n",
      "idx is 253\n",
      "3482\n",
      "idx is 254\n",
      "3417\n",
      "idx is 255\n",
      "2536\n",
      "idx is 256\n",
      "3652\n",
      "idx is 257\n",
      "3876\n",
      "idx is 258\n",
      "3815\n",
      "idx is 259\n",
      "4889\n",
      "idx is 260\n",
      "2729\n",
      "idx is 261\n",
      "3146\n",
      "idx is 262\n",
      "3240\n",
      "idx is 263\n",
      "3698\n",
      "idx is 264\n",
      "3545\n",
      "idx is 265\n",
      "2710\n",
      "idx is 266\n",
      "2778\n",
      "idx is 267\n",
      "4278\n",
      "idx is 268\n",
      "2799\n",
      "idx is 269\n",
      "4647\n",
      "idx is 270\n",
      "3226\n",
      "idx is 271\n",
      "3311\n",
      "idx is 272\n",
      "3748\n",
      "idx is 273\n",
      "4270\n",
      "idx is 274\n",
      "3571\n",
      "idx is 275\n",
      "4910\n",
      "idx is 276\n",
      "4043\n",
      "idx is 277\n",
      "2852\n",
      "idx is 278\n",
      "3298\n",
      "idx is 279\n",
      "4430\n",
      "idx is 280\n",
      "4799\n",
      "idx is 281\n",
      "4893\n",
      "idx is 282\n",
      "4525\n",
      "idx is 283\n",
      "2614\n",
      "idx is 284\n",
      "3762\n",
      "idx is 285\n",
      "3779\n",
      "idx is 286\n",
      "3506\n",
      "idx is 287\n",
      "2712\n",
      "idx is 288\n",
      "3248\n",
      "idx is 289\n",
      "3494\n",
      "idx is 290\n",
      "2724\n",
      "idx is 291\n",
      "4647\n",
      "idx is 292\n",
      "3182\n",
      "idx is 293\n",
      "4679\n",
      "idx is 294\n",
      "3236\n",
      "idx is 295\n",
      "3982\n",
      "idx is 296\n",
      "4141\n",
      "idx is 297\n",
      "3196\n",
      "idx is 298\n",
      "3058\n",
      "idx is 299\n",
      "3265\n",
      "idx is 300\n",
      "3523\n",
      "idx is 301\n",
      "3944\n",
      "idx is 302\n",
      "3711\n",
      "idx is 303\n",
      "3998\n",
      "idx is 304\n",
      "4161\n",
      "idx is 305\n",
      "4171\n",
      "idx is 306\n",
      "2573\n",
      "idx is 307\n",
      "4801\n",
      "idx is 308\n",
      "2780\n",
      "idx is 309\n",
      "2787\n",
      "idx is 310\n",
      "4913\n",
      "idx is 311\n",
      "3673\n",
      "idx is 312\n",
      "3941\n",
      "idx is 313\n",
      "4147\n",
      "idx is 314\n",
      "4219\n",
      "idx is 315\n",
      "4831\n",
      "idx is 316\n",
      "2589\n",
      "idx is 317\n",
      "4710\n",
      "idx is 318\n",
      "2763\n",
      "idx is 319\n",
      "2787\n",
      "idx is 320\n",
      "3378\n",
      "idx is 321\n",
      "2968\n",
      "idx is 322\n",
      "2548\n",
      "idx is 323\n",
      "4030\n",
      "idx is 324\n",
      "2587\n",
      "idx is 325\n",
      "4143\n",
      "idx is 326\n",
      "2941\n",
      "idx is 327\n",
      "4861\n",
      "idx is 328\n",
      "3381\n",
      "idx is 329\n",
      "2773\n",
      "idx is 330\n",
      "3305\n",
      "idx is 331\n",
      "3677\n",
      "idx is 332\n",
      "3610\n",
      "idx is 333\n",
      "3198\n",
      "idx is 334\n",
      "3600\n",
      "idx is 335\n",
      "4895\n",
      "idx is 336\n",
      "2526\n",
      "idx is 337\n",
      "2547\n",
      "idx is 338\n",
      "4375\n",
      "idx is 339\n",
      "4223\n",
      "idx is 340\n",
      "3354\n",
      "idx is 341\n",
      "4269\n",
      "idx is 342\n",
      "3322\n",
      "idx is 343\n",
      "3480\n",
      "idx is 344\n",
      "2670\n",
      "idx is 345\n",
      "4607\n",
      "idx is 346\n",
      "3980\n",
      "idx is 347\n",
      "3653\n",
      "idx is 348\n",
      "2661\n",
      "idx is 349\n",
      "4992\n",
      "idx is 350\n",
      "4327\n",
      "idx is 351\n",
      "3270\n",
      "idx is 352\n",
      "4136\n",
      "idx is 353\n",
      "4898\n",
      "idx is 354\n",
      "3945\n",
      "idx is 355\n",
      "4701\n",
      "idx is 356\n",
      "4327\n",
      "idx is 357\n",
      "4521\n",
      "idx is 358\n",
      "4460\n",
      "idx is 359\n",
      "3846\n",
      "idx is 360\n",
      "3703\n",
      "idx is 361\n",
      "3697\n",
      "idx is 362\n",
      "4443\n",
      "idx is 363\n",
      "4031\n",
      "idx is 364\n",
      "3438\n",
      "idx is 365\n",
      "3409\n",
      "idx is 366\n",
      "4970\n",
      "idx is 367\n",
      "2969\n",
      "idx is 368\n",
      "2861\n",
      "idx is 369\n",
      "3993\n",
      "idx is 370\n",
      "3525\n",
      "idx is 371\n",
      "4719\n",
      "idx is 372\n",
      "4204\n",
      "idx is 373\n",
      "4886\n",
      "idx is 374\n",
      "3234\n",
      "idx is 375\n",
      "3994\n",
      "idx is 376\n",
      "3756\n",
      "idx is 377\n",
      "3803\n",
      "idx is 378\n",
      "4538\n",
      "idx is 379\n",
      "4959\n",
      "idx is 380\n",
      "2600\n",
      "idx is 381\n",
      "3487\n",
      "idx is 382\n",
      "4247\n",
      "idx is 383\n",
      "3768\n",
      "idx is 384\n",
      "2897\n",
      "idx is 385\n",
      "3122\n",
      "idx is 386\n",
      "4184\n",
      "idx is 387\n",
      "2825\n",
      "idx is 388\n",
      "4613\n",
      "idx is 389\n",
      "4372\n",
      "idx is 390\n",
      "3063\n",
      "idx is 391\n",
      "3050\n",
      "idx is 392\n",
      "3720\n",
      "idx is 393\n",
      "2790\n",
      "idx is 394\n",
      "4677\n",
      "idx is 395\n",
      "4839\n",
      "idx is 396\n",
      "3516\n",
      "idx is 397\n",
      "4483\n",
      "idx is 398\n",
      "3026\n",
      "idx is 399\n",
      "3626\n",
      "idx is 400\n",
      "4561\n",
      "idx is 401\n",
      "4261\n",
      "idx is 402\n",
      "2808\n",
      "idx is 403\n",
      "3789\n",
      "idx is 404\n",
      "3027\n",
      "idx is 405\n",
      "4204\n",
      "idx is 406\n",
      "4561\n",
      "idx is 407\n",
      "4032\n",
      "idx is 408\n",
      "4961\n",
      "idx is 409\n",
      "3725\n",
      "idx is 410\n",
      "4869\n",
      "idx is 411\n",
      "4290\n",
      "idx is 412\n",
      "3434\n",
      "idx is 413\n",
      "3101\n",
      "idx is 414\n",
      "4267\n",
      "idx is 415\n",
      "3292\n",
      "idx is 416\n",
      "4410\n",
      "idx is 417\n",
      "4155\n",
      "idx is 418\n",
      "2559\n",
      "idx is 419\n",
      "4097\n",
      "idx is 420\n",
      "3004\n",
      "idx is 421\n",
      "4546\n",
      "idx is 422\n",
      "3874\n",
      "idx is 423\n",
      "4650\n",
      "idx is 424\n",
      "2934\n",
      "idx is 425\n",
      "3839\n",
      "idx is 426\n",
      "3208\n",
      "idx is 427\n",
      "4701\n",
      "idx is 428\n",
      "4177\n",
      "idx is 429\n",
      "3244\n",
      "idx is 430\n",
      "3191\n",
      "idx is 431\n",
      "3456\n",
      "idx is 432\n",
      "2832\n",
      "idx is 433\n",
      "4378\n",
      "idx is 434\n",
      "2879\n",
      "idx is 435\n",
      "4265\n",
      "idx is 436\n",
      "3426\n",
      "idx is 437\n",
      "4259\n",
      "idx is 438\n",
      "3324\n",
      "idx is 439\n",
      "3775\n",
      "idx is 440\n",
      "4713\n",
      "idx is 441\n",
      "2924\n",
      "idx is 442\n",
      "3666\n",
      "idx is 443\n",
      "4349\n",
      "idx is 444\n",
      "3920\n",
      "idx is 445\n",
      "4707\n",
      "idx is 446\n",
      "2777\n",
      "idx is 447\n",
      "3188\n",
      "idx is 448\n",
      "4054\n",
      "idx is 449\n",
      "3754\n",
      "idx is 450\n",
      "2560\n",
      "idx is 451\n",
      "4579\n",
      "idx is 452\n",
      "2920\n",
      "idx is 453\n",
      "3171\n",
      "idx is 454\n",
      "3755\n",
      "idx is 455\n",
      "4508\n",
      "idx is 456\n",
      "4582\n",
      "idx is 457\n",
      "2881\n",
      "idx is 458\n",
      "3188\n",
      "idx is 459\n",
      "4734\n",
      "idx is 460\n",
      "4144\n",
      "idx is 461\n",
      "2871\n",
      "idx is 462\n",
      "3788\n",
      "idx is 463\n",
      "2556\n",
      "idx is 464\n",
      "3706\n",
      "idx is 465\n",
      "3162\n",
      "idx is 466\n",
      "4807\n",
      "idx is 467\n",
      "3200\n",
      "idx is 468\n",
      "4655\n",
      "idx is 469\n",
      "4656\n",
      "idx is 470\n",
      "3498\n",
      "idx is 471\n",
      "4376\n",
      "idx is 472\n",
      "3244\n",
      "idx is 473\n",
      "3591\n",
      "idx is 474\n",
      "2584\n",
      "idx is 475\n",
      "3009\n",
      "idx is 476\n",
      "4623\n",
      "idx is 477\n",
      "4870\n",
      "idx is 478\n",
      "4098\n",
      "idx is 479\n",
      "3644\n",
      "idx is 480\n",
      "2513\n",
      "idx is 481\n",
      "3839\n",
      "idx is 482\n",
      "3361\n",
      "idx is 483\n",
      "4059\n",
      "idx is 484\n",
      "3803\n",
      "idx is 485\n",
      "3200\n",
      "idx is 486\n",
      "4543\n",
      "idx is 487\n",
      "4030\n",
      "idx is 488\n",
      "4114\n",
      "idx is 489\n",
      "2961\n",
      "idx is 490\n",
      "4876\n",
      "idx is 491\n",
      "3532\n",
      "idx is 492\n",
      "4192\n",
      "idx is 493\n",
      "2570\n",
      "idx is 494\n",
      "4604\n",
      "idx is 495\n",
      "3556\n",
      "idx is 496\n",
      "4044\n",
      "idx is 497\n",
      "2783\n",
      "idx is 498\n",
      "2520\n",
      "idx is 499\n",
      "3328\n",
      "idx is 500\n",
      "4259\n",
      "idx is 501\n",
      "4960\n",
      "idx is 502\n",
      "3588\n",
      "idx is 503\n",
      "4796\n",
      "idx is 504\n",
      "3279\n",
      "idx is 505\n",
      "4154\n",
      "idx is 506\n",
      "4645\n",
      "idx is 507\n",
      "4329\n",
      "idx is 508\n",
      "3378\n",
      "idx is 509\n",
      "4580\n",
      "idx is 510\n",
      "2644\n",
      "idx is 511\n",
      "3695\n",
      "idx is 512\n",
      "3920\n",
      "idx is 513\n",
      "2626\n",
      "idx is 514\n",
      "3113\n",
      "idx is 515\n",
      "3753\n",
      "idx is 516\n",
      "4546\n",
      "idx is 517\n",
      "3759\n",
      "idx is 518\n",
      "2537\n",
      "idx is 519\n",
      "3743\n",
      "idx is 520\n",
      "2920\n",
      "idx is 521\n",
      "4536\n",
      "idx is 522\n",
      "4767\n",
      "idx is 523\n",
      "3821\n",
      "idx is 524\n",
      "2944\n",
      "idx is 525\n",
      "3899\n",
      "idx is 526\n",
      "3442\n",
      "idx is 527\n",
      "3953\n",
      "idx is 528\n",
      "2900\n",
      "idx is 529\n",
      "4357\n",
      "idx is 530\n",
      "4240\n",
      "idx is 531\n",
      "3439\n",
      "idx is 532\n",
      "3430\n",
      "idx is 533\n",
      "4569\n",
      "idx is 534\n",
      "3960\n",
      "idx is 535\n",
      "4367\n",
      "idx is 536\n",
      "4287\n",
      "idx is 537\n",
      "3777\n",
      "idx is 538\n",
      "3559\n",
      "idx is 539\n",
      "2971\n",
      "idx is 540\n",
      "3364\n",
      "idx is 541\n",
      "3040\n",
      "idx is 542\n",
      "4314\n",
      "idx is 543\n",
      "4852\n",
      "idx is 544\n",
      "2785\n",
      "idx is 545\n",
      "4562\n",
      "idx is 546\n",
      "3913\n",
      "idx is 547\n",
      "2519\n",
      "idx is 548\n",
      "3180\n",
      "idx is 549\n",
      "4316\n",
      "idx is 550\n",
      "2534\n",
      "idx is 551\n",
      "4194\n",
      "idx is 552\n",
      "3941\n",
      "idx is 553\n",
      "3539\n",
      "idx is 554\n",
      "4005\n",
      "idx is 555\n",
      "2512\n",
      "idx is 556\n",
      "2769\n",
      "idx is 557\n",
      "2882\n",
      "idx is 558\n",
      "3166\n",
      "idx is 559\n",
      "2568\n",
      "idx is 560\n",
      "2941\n",
      "idx is 561\n",
      "2857\n",
      "idx is 562\n",
      "2501\n",
      "idx is 563\n",
      "2968\n",
      "idx is 564\n",
      "3147\n",
      "idx is 565\n",
      "3311\n",
      "idx is 566\n",
      "3955\n",
      "idx is 567\n",
      "3944\n",
      "idx is 568\n",
      "4207\n",
      "idx is 569\n",
      "3050\n",
      "idx is 570\n",
      "4931\n",
      "idx is 571\n",
      "4941\n",
      "idx is 572\n",
      "4585\n",
      "idx is 573\n",
      "3821\n",
      "idx is 574\n",
      "3265\n",
      "idx is 575\n",
      "4407\n",
      "idx is 576\n",
      "2594\n",
      "idx is 577\n",
      "2743\n",
      "idx is 578\n",
      "4785\n",
      "idx is 579\n",
      "3331\n",
      "idx is 580\n",
      "4411\n",
      "idx is 581\n",
      "3248\n",
      "idx is 582\n",
      "3913\n",
      "idx is 583\n",
      "3410\n",
      "idx is 584\n",
      "4282\n",
      "idx is 585\n",
      "4748\n",
      "idx is 586\n",
      "3081\n",
      "idx is 587\n",
      "4704\n",
      "idx is 588\n",
      "3655\n",
      "idx is 589\n",
      "2866\n",
      "idx is 590\n",
      "4507\n",
      "idx is 591\n",
      "4962\n",
      "idx is 592\n",
      "4737\n",
      "idx is 593\n",
      "4483\n",
      "idx is 594\n",
      "3840\n",
      "idx is 595\n",
      "3855\n",
      "idx is 596\n",
      "2647\n",
      "idx is 597\n",
      "3234\n",
      "idx is 598\n",
      "3886\n",
      "idx is 599\n",
      "4236\n",
      "idx is 600\n",
      "2791\n",
      "idx is 601\n",
      "3283\n",
      "idx is 602\n",
      "4462\n",
      "idx is 603\n",
      "2580\n",
      "idx is 604\n",
      "4989\n",
      "idx is 605\n",
      "4779\n",
      "idx is 606\n",
      "3446\n",
      "idx is 607\n",
      "3288\n",
      "idx is 608\n",
      "4103\n",
      "idx is 609\n",
      "3861\n",
      "idx is 610\n",
      "3101\n",
      "idx is 611\n",
      "4596\n",
      "idx is 612\n",
      "4785\n",
      "idx is 613\n",
      "3370\n",
      "idx is 614\n",
      "3131\n",
      "idx is 615\n",
      "3152\n",
      "idx is 616\n",
      "4579\n",
      "idx is 617\n",
      "3897\n",
      "idx is 618\n",
      "2554\n",
      "idx is 619\n",
      "3846\n",
      "idx is 620\n",
      "4516\n",
      "idx is 621\n",
      "3293\n",
      "idx is 622\n",
      "3900\n",
      "idx is 623\n",
      "4798\n",
      "idx is 624\n",
      "3955\n",
      "idx is 625\n",
      "3188\n",
      "idx is 626\n",
      "4688\n",
      "idx is 627\n",
      "4504\n",
      "idx is 628\n",
      "4475\n",
      "idx is 629\n",
      "3644\n",
      "idx is 630\n",
      "3215\n",
      "idx is 631\n",
      "3998\n",
      "idx is 632\n",
      "2563\n",
      "idx is 633\n",
      "3760\n",
      "idx is 634\n",
      "3848\n",
      "idx is 635\n",
      "4355\n",
      "idx is 636\n",
      "4883\n",
      "idx is 637\n",
      "2860\n",
      "idx is 638\n",
      "3940\n",
      "idx is 639\n",
      "4918\n",
      "idx is 640\n",
      "3372\n",
      "idx is 641\n",
      "4450\n",
      "idx is 642\n",
      "4414\n",
      "idx is 643\n",
      "3551\n",
      "idx is 644\n",
      "3704\n",
      "idx is 645\n",
      "4277\n",
      "idx is 646\n",
      "2966\n",
      "idx is 647\n",
      "3743\n",
      "idx is 648\n",
      "2859\n",
      "idx is 649\n",
      "2834\n",
      "idx is 650\n",
      "3978\n",
      "idx is 651\n",
      "3991\n",
      "idx is 652\n",
      "3914\n",
      "idx is 653\n",
      "3329\n",
      "idx is 654\n",
      "2830\n",
      "idx is 655\n",
      "4071\n",
      "idx is 656\n",
      "3635\n",
      "idx is 657\n",
      "3539\n",
      "idx is 658\n",
      "4935\n",
      "idx is 659\n",
      "3545\n",
      "idx is 660\n",
      "3433\n",
      "idx is 661\n",
      "2724\n",
      "idx is 662\n",
      "2864\n",
      "idx is 663\n",
      "4339\n",
      "idx is 664\n",
      "4694\n",
      "idx is 665\n",
      "3103\n",
      "idx is 666\n",
      "2541\n",
      "idx is 667\n",
      "3123\n",
      "idx is 668\n",
      "3834\n",
      "idx is 669\n",
      "4134\n",
      "idx is 670\n",
      "3508\n",
      "idx is 671\n",
      "4149\n",
      "idx is 672\n",
      "2764\n",
      "idx is 673\n",
      "2731\n",
      "idx is 674\n",
      "3555\n",
      "idx is 675\n",
      "3722\n",
      "idx is 676\n",
      "4224\n",
      "idx is 677\n",
      "4861\n",
      "idx is 678\n",
      "3066\n",
      "idx is 679\n",
      "4500\n",
      "idx is 680\n",
      "4459\n",
      "idx is 681\n",
      "4489\n",
      "idx is 682\n",
      "4825\n",
      "idx is 683\n",
      "3386\n",
      "idx is 684\n",
      "4911\n",
      "idx is 685\n",
      "3055\n",
      "idx is 686\n",
      "2709\n",
      "idx is 687\n",
      "4425\n",
      "idx is 688\n",
      "3844\n",
      "idx is 689\n",
      "4684\n",
      "idx is 690\n",
      "4885\n",
      "idx is 691\n",
      "2711\n",
      "idx is 692\n",
      "3138\n",
      "idx is 693\n",
      "2523\n",
      "idx is 694\n",
      "3870\n",
      "idx is 695\n",
      "3765\n",
      "idx is 696\n",
      "4208\n",
      "idx is 697\n",
      "2842\n",
      "idx is 698\n",
      "3254\n",
      "idx is 699\n",
      "2627\n",
      "idx is 700\n",
      "4338\n",
      "idx is 701\n",
      "3639\n",
      "idx is 702\n",
      "4930\n",
      "idx is 703\n",
      "3568\n",
      "idx is 704\n",
      "3869\n",
      "idx is 705\n",
      "2618\n",
      "idx is 706\n",
      "4900\n",
      "idx is 707\n",
      "2554\n",
      "idx is 708\n",
      "4487\n",
      "idx is 709\n",
      "4742\n",
      "idx is 710\n",
      "2952\n",
      "idx is 711\n",
      "2554\n",
      "idx is 712\n",
      "3488\n",
      "idx is 713\n",
      "3883\n",
      "idx is 714\n",
      "4395\n",
      "idx is 715\n",
      "4554\n",
      "idx is 716\n",
      "2527\n",
      "idx is 717\n",
      "4997\n",
      "idx is 718\n",
      "2721\n",
      "idx is 719\n",
      "3858\n",
      "idx is 720\n",
      "4810\n",
      "idx is 721\n",
      "2768\n",
      "idx is 722\n",
      "2553\n",
      "idx is 723\n",
      "4424\n",
      "idx is 724\n",
      "3600\n",
      "idx is 725\n",
      "4031\n",
      "idx is 726\n",
      "2904\n",
      "idx is 727\n",
      "4035\n",
      "idx is 728\n",
      "3937\n",
      "idx is 729\n",
      "2728\n",
      "idx is 730\n",
      "3924\n",
      "idx is 731\n",
      "3108\n",
      "idx is 732\n",
      "2756\n",
      "idx is 733\n",
      "4935\n",
      "idx is 734\n",
      "2899\n",
      "idx is 735\n",
      "2698\n",
      "idx is 736\n",
      "3770\n",
      "idx is 737\n",
      "3168\n",
      "idx is 738\n",
      "4613\n",
      "idx is 739\n",
      "4222\n",
      "idx is 740\n",
      "4640\n",
      "idx is 741\n",
      "3058\n",
      "idx is 742\n",
      "3556\n",
      "idx is 743\n",
      "3116\n",
      "idx is 744\n",
      "3890\n",
      "idx is 745\n",
      "3526\n",
      "idx is 746\n",
      "2677\n",
      "idx is 747\n",
      "2978\n",
      "idx is 748\n",
      "3783\n",
      "idx is 749\n",
      "4857\n",
      "idx is 750\n",
      "4811\n",
      "idx is 751\n",
      "2742\n",
      "idx is 752\n",
      "3179\n",
      "idx is 753\n",
      "4563\n",
      "idx is 754\n",
      "4096\n",
      "idx is 755\n",
      "2733\n",
      "idx is 756\n",
      "4395\n",
      "idx is 757\n",
      "4194\n",
      "idx is 758\n",
      "4412\n",
      "idx is 759\n",
      "2988\n",
      "idx is 760\n",
      "4796\n",
      "idx is 761\n",
      "4442\n",
      "idx is 762\n",
      "4613\n",
      "idx is 763\n",
      "4263\n",
      "idx is 764\n",
      "4547\n",
      "idx is 765\n",
      "2553\n",
      "idx is 766\n",
      "3239\n",
      "idx is 767\n",
      "4177\n",
      "idx is 768\n",
      "4823\n",
      "idx is 769\n",
      "3204\n",
      "idx is 770\n",
      "3208\n",
      "idx is 771\n",
      "4593\n",
      "idx is 772\n",
      "3130\n",
      "idx is 773\n",
      "4993\n",
      "idx is 774\n",
      "2775\n",
      "idx is 775\n",
      "4593\n",
      "idx is 776\n",
      "3217\n",
      "idx is 777\n",
      "2683\n",
      "idx is 778\n",
      "4344\n",
      "idx is 779\n",
      "4147\n",
      "idx is 780\n",
      "3258\n",
      "idx is 781\n",
      "4348\n",
      "idx is 782\n",
      "3008\n",
      "idx is 783\n",
      "3173\n",
      "idx is 784\n",
      "3789\n",
      "idx is 785\n",
      "4366\n",
      "idx is 786\n",
      "3998\n",
      "idx is 787\n",
      "4143\n",
      "idx is 788\n",
      "3049\n",
      "idx is 789\n",
      "2898\n",
      "idx is 790\n",
      "4912\n",
      "idx is 791\n",
      "2587\n",
      "idx is 792\n",
      "3948\n",
      "idx is 793\n",
      "3467\n",
      "idx is 794\n",
      "2893\n",
      "idx is 795\n",
      "3512\n",
      "idx is 796\n",
      "4783\n",
      "idx is 797\n",
      "3155\n",
      "idx is 798\n",
      "2998\n",
      "idx is 799\n",
      "3807\n",
      "idx is 800\n",
      "4171\n",
      "idx is 801\n",
      "2660\n",
      "idx is 802\n",
      "2646\n",
      "idx is 803\n",
      "4434\n",
      "idx is 804\n",
      "3013\n",
      "idx is 805\n",
      "4535\n",
      "idx is 806\n",
      "4639\n",
      "idx is 807\n",
      "3905\n",
      "idx is 808\n",
      "3102\n",
      "idx is 809\n",
      "4375\n",
      "idx is 810\n",
      "3571\n",
      "idx is 811\n",
      "4851\n",
      "idx is 812\n",
      "4593\n",
      "idx is 813\n",
      "4937\n",
      "idx is 814\n",
      "3139\n",
      "idx is 815\n",
      "3368\n",
      "idx is 816\n",
      "4586\n",
      "idx is 817\n",
      "2973\n",
      "idx is 818\n",
      "3913\n",
      "idx is 819\n",
      "3663\n",
      "idx is 820\n",
      "4994\n",
      "idx is 821\n",
      "3765\n",
      "idx is 822\n",
      "4446\n",
      "idx is 823\n",
      "4256\n",
      "idx is 824\n",
      "4701\n",
      "idx is 825\n",
      "4881\n",
      "idx is 826\n",
      "4703\n",
      "idx is 827\n",
      "2726\n",
      "idx is 828\n",
      "3955\n",
      "idx is 829\n",
      "3458\n",
      "idx is 830\n",
      "3217\n",
      "idx is 831\n",
      "3003\n",
      "idx is 832\n",
      "2787\n",
      "idx is 833\n",
      "3116\n",
      "idx is 834\n",
      "2718\n",
      "idx is 835\n",
      "4567\n",
      "idx is 836\n",
      "3230\n",
      "idx is 837\n",
      "2714\n",
      "idx is 838\n",
      "3726\n",
      "idx is 839\n",
      "4608\n",
      "idx is 840\n",
      "3565\n",
      "idx is 841\n",
      "3330\n",
      "idx is 842\n",
      "2900\n",
      "idx is 843\n",
      "4319\n",
      "idx is 844\n",
      "2509\n",
      "idx is 845\n",
      "2756\n",
      "idx is 846\n",
      "4121\n",
      "idx is 847\n",
      "3915\n",
      "idx is 848\n",
      "2667\n",
      "idx is 849\n",
      "3830\n",
      "idx is 850\n",
      "4423\n",
      "idx is 851\n",
      "2567\n",
      "idx is 852\n",
      "4586\n",
      "idx is 853\n",
      "3414\n",
      "idx is 854\n",
      "2847\n",
      "idx is 855\n",
      "2716\n",
      "idx is 856\n",
      "3477\n",
      "idx is 857\n",
      "4985\n",
      "idx is 858\n",
      "2824\n",
      "idx is 859\n",
      "4610\n",
      "idx is 860\n",
      "4240\n",
      "idx is 861\n",
      "3700\n",
      "idx is 862\n",
      "4301\n",
      "idx is 863\n",
      "4012\n",
      "idx is 864\n",
      "4311\n",
      "idx is 865\n",
      "4793\n",
      "idx is 866\n",
      "2578\n",
      "idx is 867\n",
      "2778\n",
      "idx is 868\n",
      "4185\n",
      "idx is 869\n",
      "4438\n",
      "idx is 870\n",
      "3576\n",
      "idx is 871\n",
      "4091\n",
      "idx is 872\n",
      "3888\n",
      "idx is 873\n",
      "3170\n",
      "idx is 874\n",
      "2967\n",
      "idx is 875\n",
      "3580\n",
      "idx is 876\n",
      "3245\n",
      "idx is 877\n",
      "2601\n",
      "idx is 878\n",
      "3881\n",
      "idx is 879\n",
      "3009\n",
      "idx is 880\n",
      "3916\n",
      "idx is 881\n",
      "4439\n",
      "idx is 882\n",
      "4099\n",
      "idx is 883\n",
      "4427\n",
      "idx is 884\n",
      "4178\n",
      "idx is 885\n",
      "2544\n",
      "idx is 886\n",
      "4577\n",
      "idx is 887\n",
      "3300\n",
      "idx is 888\n",
      "4571\n",
      "idx is 889\n",
      "2735\n",
      "idx is 890\n",
      "4432\n",
      "idx is 891\n",
      "4898\n",
      "idx is 892\n",
      "3892\n",
      "idx is 893\n",
      "2861\n",
      "idx is 894\n",
      "3093\n",
      "idx is 895\n",
      "4573\n",
      "idx is 896\n",
      "4867\n",
      "idx is 897\n",
      "3833\n",
      "idx is 898\n",
      "4305\n",
      "idx is 899\n",
      "4146\n",
      "idx is 900\n",
      "3228\n",
      "idx is 901\n",
      "4957\n",
      "idx is 902\n",
      "4887\n",
      "idx is 903\n",
      "3387\n",
      "idx is 904\n",
      "4940\n",
      "idx is 905\n",
      "3550\n",
      "idx is 906\n",
      "3399\n",
      "idx is 907\n",
      "2583\n",
      "idx is 908\n",
      "4901\n",
      "idx is 909\n",
      "3116\n",
      "idx is 910\n",
      "2815\n",
      "idx is 911\n",
      "3414\n",
      "idx is 912\n",
      "3731\n",
      "idx is 913\n",
      "3594\n",
      "idx is 914\n",
      "3482\n",
      "idx is 915\n",
      "4408\n",
      "idx is 916\n",
      "2542\n",
      "idx is 917\n",
      "2688\n",
      "idx is 918\n",
      "2514\n",
      "idx is 919\n",
      "3948\n",
      "idx is 920\n",
      "3988\n",
      "idx is 921\n",
      "3526\n",
      "idx is 922\n",
      "3122\n",
      "idx is 923\n",
      "4637\n",
      "idx is 924\n",
      "3245\n",
      "idx is 925\n",
      "3839\n",
      "idx is 926\n",
      "3463\n",
      "idx is 927\n",
      "4173\n",
      "idx is 928\n",
      "4429\n",
      "idx is 929\n",
      "4420\n",
      "idx is 930\n",
      "2568\n",
      "idx is 931\n",
      "4706\n",
      "idx is 932\n",
      "2569\n",
      "idx is 933\n",
      "4221\n",
      "idx is 934\n",
      "3240\n",
      "idx is 935\n",
      "2987\n",
      "idx is 936\n",
      "2524\n",
      "idx is 937\n",
      "4217\n",
      "idx is 938\n",
      "3236\n",
      "idx is 939\n",
      "4935\n",
      "idx is 940\n",
      "4048\n",
      "idx is 941\n",
      "4603\n",
      "idx is 942\n",
      "3955\n",
      "idx is 943\n",
      "3294\n",
      "idx is 944\n",
      "2854\n",
      "idx is 945\n",
      "4818\n",
      "idx is 946\n",
      "4730\n",
      "idx is 947\n",
      "4299\n",
      "idx is 948\n",
      "4933\n",
      "idx is 949\n",
      "4687\n",
      "idx is 950\n",
      "3989\n",
      "idx is 951\n",
      "4086\n",
      "idx is 952\n",
      "4724\n",
      "idx is 953\n",
      "4513\n",
      "idx is 954\n",
      "4910\n",
      "idx is 955\n",
      "2659\n",
      "idx is 956\n",
      "3323\n",
      "idx is 957\n",
      "3439\n",
      "idx is 958\n",
      "3383\n",
      "idx is 959\n",
      "3573\n",
      "idx is 960\n",
      "4457\n",
      "idx is 961\n",
      "3396\n",
      "idx is 962\n",
      "3155\n",
      "idx is 963\n",
      "4140\n",
      "idx is 964\n",
      "4765\n",
      "idx is 965\n",
      "3683\n",
      "idx is 966\n",
      "3375\n",
      "idx is 967\n",
      "4323\n",
      "idx is 968\n",
      "3186\n",
      "idx is 969\n",
      "4965\n",
      "idx is 970\n",
      "4265\n",
      "idx is 971\n",
      "3220\n",
      "idx is 972\n",
      "4459\n",
      "idx is 973\n",
      "4138\n",
      "idx is 974\n",
      "3728\n",
      "idx is 975\n",
      "4239\n",
      "idx is 976\n",
      "4872\n",
      "idx is 977\n",
      "4892\n",
      "idx is 978\n",
      "3443\n",
      "idx is 979\n",
      "3159\n",
      "idx is 980\n",
      "2980\n",
      "idx is 981\n",
      "4838\n",
      "idx is 982\n",
      "4029\n",
      "idx is 983\n",
      "4777\n",
      "idx is 984\n",
      "3276\n",
      "idx is 985\n",
      "4793\n",
      "idx is 986\n",
      "3404\n",
      "idx is 987\n",
      "4098\n",
      "idx is 988\n",
      "2633\n",
      "idx is 989\n",
      "2968\n",
      "idx is 990\n",
      "3670\n",
      "idx is 991\n",
      "4113\n",
      "idx is 992\n",
      "4958\n",
      "idx is 993\n",
      "3219\n",
      "idx is 994\n",
      "3069\n",
      "idx is 995\n",
      "4336\n",
      "idx is 996\n",
      "4666\n",
      "idx is 997\n",
      "2540\n",
      "idx is 998\n",
      "4818\n",
      "idx is 999\n",
      "3549\n",
      "idx is 1000\n",
      "2688\n",
      "idx is 1001\n",
      "3313\n",
      "idx is 1002\n",
      "3695\n",
      "idx is 1003\n",
      "4647\n",
      "idx is 1004\n",
      "3445\n",
      "idx is 1005\n",
      "4204\n",
      "idx is 1006\n",
      "2869\n",
      "idx is 1007\n",
      "3316\n",
      "idx is 1008\n",
      "2836\n",
      "idx is 1009\n",
      "4212\n",
      "idx is 1010\n",
      "4829\n",
      "idx is 1011\n",
      "4726\n",
      "idx is 1012\n",
      "4742\n",
      "idx is 1013\n",
      "4330\n",
      "idx is 1014\n",
      "4682\n",
      "idx is 1015\n",
      "4840\n",
      "idx is 1016\n",
      "3601\n",
      "idx is 1017\n",
      "3969\n",
      "idx is 1018\n",
      "3484\n",
      "idx is 1019\n",
      "3014\n",
      "idx is 1020\n",
      "3254\n",
      "idx is 1021\n",
      "3299\n",
      "idx is 1022\n",
      "4829\n",
      "idx is 1023\n",
      "3390\n",
      "idx is 1024\n",
      "4465\n",
      "idx is 1025\n",
      "2500\n",
      "idx is 1026\n",
      "3841\n",
      "idx is 1027\n",
      "3086\n",
      "idx is 1028\n",
      "4215\n",
      "idx is 1029\n",
      "2992\n",
      "idx is 1030\n",
      "4128\n",
      "idx is 1031\n",
      "4064\n",
      "idx is 1032\n",
      "2736\n",
      "idx is 1033\n",
      "3893\n",
      "idx is 1034\n",
      "3327\n",
      "idx is 1035\n",
      "3379\n",
      "idx is 1036\n",
      "4751\n",
      "idx is 1037\n",
      "3656\n",
      "idx is 1038\n",
      "4378\n",
      "idx is 1039\n",
      "4342\n",
      "idx is 1040\n",
      "4283\n",
      "idx is 1041\n",
      "2859\n",
      "idx is 1042\n",
      "3237\n",
      "idx is 1043\n",
      "4352\n",
      "idx is 1044\n",
      "2536\n",
      "idx is 1045\n",
      "2571\n",
      "idx is 1046\n",
      "4998\n",
      "idx is 1047\n",
      "2780\n",
      "idx is 1048\n",
      "2650\n",
      "idx is 1049\n",
      "4997\n",
      "idx is 1050\n",
      "4166\n",
      "idx is 1051\n",
      "4851\n",
      "idx is 1052\n",
      "4258\n",
      "idx is 1053\n",
      "4425\n",
      "idx is 1054\n",
      "3368\n",
      "idx is 1055\n",
      "3877\n",
      "idx is 1056\n",
      "2983\n",
      "idx is 1057\n",
      "2930\n",
      "idx is 1058\n",
      "3423\n",
      "idx is 1059\n",
      "2950\n",
      "idx is 1060\n",
      "4849\n",
      "idx is 1061\n",
      "2865\n",
      "idx is 1062\n",
      "3786\n",
      "idx is 1063\n",
      "3430\n",
      "idx is 1064\n",
      "3531\n",
      "idx is 1065\n",
      "3513\n",
      "idx is 1066\n",
      "4543\n",
      "idx is 1067\n",
      "4919\n",
      "idx is 1068\n",
      "4129\n",
      "idx is 1069\n",
      "2591\n",
      "idx is 1070\n",
      "4528\n",
      "idx is 1071\n",
      "2709\n",
      "idx is 1072\n",
      "2666\n",
      "idx is 1073\n",
      "3954\n",
      "idx is 1074\n",
      "4086\n",
      "idx is 1075\n",
      "2607\n",
      "idx is 1076\n",
      "3954\n",
      "idx is 1077\n",
      "3193\n",
      "idx is 1078\n",
      "4603\n",
      "idx is 1079\n",
      "4330\n",
      "idx is 1080\n",
      "4831\n",
      "idx is 1081\n",
      "4391\n",
      "idx is 1082\n",
      "2816\n",
      "idx is 1083\n",
      "4076\n",
      "idx is 1084\n",
      "4574\n",
      "idx is 1085\n",
      "2866\n",
      "idx is 1086\n",
      "4530\n",
      "idx is 1087\n",
      "4197\n",
      "idx is 1088\n",
      "3696\n",
      "idx is 1089\n",
      "4280\n",
      "idx is 1090\n",
      "4821\n",
      "idx is 1091\n",
      "4101\n",
      "idx is 1092\n",
      "4345\n",
      "idx is 1093\n",
      "3398\n",
      "idx is 1094\n",
      "4596\n",
      "idx is 1095\n",
      "2605\n",
      "idx is 1096\n",
      "2746\n",
      "idx is 1097\n",
      "4010\n",
      "idx is 1098\n",
      "3523\n",
      "idx is 1099\n",
      "3298\n",
      "idx is 1100\n",
      "4829\n",
      "idx is 1101\n",
      "2999\n",
      "idx is 1102\n",
      "4481\n",
      "idx is 1103\n",
      "3560\n",
      "idx is 1104\n",
      "3084\n",
      "idx is 1105\n",
      "4229\n",
      "idx is 1106\n",
      "3900\n",
      "idx is 1107\n",
      "4120\n",
      "idx is 1108\n",
      "4244\n",
      "idx is 1109\n",
      "4026\n",
      "idx is 1110\n",
      "4017\n",
      "idx is 1111\n",
      "3862\n",
      "idx is 1112\n",
      "3104\n",
      "idx is 1113\n",
      "4575\n",
      "idx is 1114\n",
      "3461\n",
      "idx is 1115\n",
      "4508\n",
      "idx is 1116\n",
      "3958\n",
      "idx is 1117\n",
      "2816\n",
      "idx is 1118\n",
      "4184\n",
      "idx is 1119\n",
      "2864\n",
      "idx is 1120\n",
      "3383\n",
      "idx is 1121\n",
      "4480\n",
      "idx is 1122\n",
      "4789\n",
      "idx is 1123\n",
      "3829\n",
      "idx is 1124\n",
      "2569\n",
      "idx is 1125\n",
      "4204\n",
      "idx is 1126\n",
      "3939\n",
      "idx is 1127\n",
      "3576\n",
      "idx is 1128\n",
      "2907\n",
      "idx is 1129\n",
      "4884\n",
      "idx is 1130\n",
      "3644\n",
      "idx is 1131\n",
      "2858\n",
      "idx is 1132\n",
      "3689\n",
      "idx is 1133\n",
      "2699\n",
      "idx is 1134\n",
      "2508\n",
      "idx is 1135\n",
      "2775\n",
      "idx is 1136\n",
      "2945\n",
      "idx is 1137\n",
      "4858\n",
      "idx is 1138\n",
      "4781\n",
      "idx is 1139\n",
      "4597\n",
      "idx is 1140\n",
      "2775\n",
      "idx is 1141\n",
      "3514\n",
      "idx is 1142\n",
      "4231\n",
      "idx is 1143\n",
      "4473\n",
      "idx is 1144\n",
      "2717\n",
      "idx is 1145\n",
      "3057\n",
      "idx is 1146\n",
      "3865\n",
      "idx is 1147\n",
      "3220\n",
      "idx is 1148\n",
      "4795\n",
      "idx is 1149\n",
      "2794\n",
      "idx is 1150\n",
      "3542\n",
      "idx is 1151\n",
      "4501\n",
      "idx is 1152\n",
      "4184\n",
      "idx is 1153\n",
      "3176\n",
      "idx is 1154\n",
      "3623\n",
      "idx is 1155\n",
      "3136\n",
      "idx is 1156\n",
      "2805\n",
      "idx is 1157\n",
      "3017\n",
      "idx is 1158\n",
      "4836\n",
      "idx is 1159\n",
      "4982\n",
      "idx is 1160\n",
      "3530\n",
      "idx is 1161\n",
      "3949\n",
      "idx is 1162\n",
      "3365\n",
      "idx is 1163\n",
      "2745\n",
      "idx is 1164\n",
      "3823\n",
      "idx is 1165\n",
      "4909\n",
      "idx is 1166\n",
      "2568\n",
      "idx is 1167\n",
      "2859\n",
      "idx is 1168\n",
      "2934\n",
      "idx is 1169\n",
      "4420\n",
      "idx is 1170\n",
      "2570\n",
      "idx is 1171\n",
      "2825\n",
      "idx is 1172\n",
      "2873\n",
      "idx is 1173\n",
      "4582\n",
      "idx is 1174\n",
      "2832\n",
      "idx is 1175\n",
      "4535\n",
      "idx is 1176\n",
      "4675\n",
      "idx is 1177\n",
      "2911\n",
      "idx is 1178\n",
      "2875\n",
      "idx is 1179\n",
      "3063\n",
      "idx is 1180\n",
      "2968\n",
      "idx is 1181\n",
      "2675\n",
      "idx is 1182\n",
      "2924\n",
      "idx is 1183\n",
      "3606\n",
      "idx is 1184\n",
      "3075\n",
      "idx is 1185\n",
      "4894\n",
      "idx is 1186\n",
      "3212\n",
      "idx is 1187\n",
      "3105\n",
      "idx is 1188\n",
      "2829\n",
      "idx is 1189\n",
      "2758\n",
      "idx is 1190\n",
      "3359\n",
      "idx is 1191\n",
      "4407\n",
      "idx is 1192\n",
      "4717\n",
      "idx is 1193\n",
      "3081\n",
      "idx is 1194\n",
      "2797\n",
      "idx is 1195\n",
      "3700\n",
      "idx is 1196\n",
      "4235\n",
      "idx is 1197\n",
      "2698\n",
      "idx is 1198\n",
      "2505\n",
      "idx is 1199\n",
      "3087\n",
      "idx is 1200\n",
      "3107\n",
      "idx is 1201\n",
      "2585\n",
      "idx is 1202\n",
      "3862\n",
      "idx is 1203\n",
      "4468\n",
      "idx is 1204\n",
      "4255\n",
      "idx is 1205\n",
      "3646\n",
      "idx is 1206\n",
      "3115\n",
      "idx is 1207\n",
      "3215\n",
      "idx is 1208\n",
      "4998\n",
      "idx is 1209\n",
      "2959\n",
      "idx is 1210\n",
      "3091\n",
      "idx is 1211\n",
      "3042\n",
      "idx is 1212\n",
      "4917\n",
      "idx is 1213\n",
      "3100\n",
      "idx is 1214\n",
      "3899\n",
      "idx is 1215\n",
      "3131\n",
      "idx is 1216\n",
      "4557\n",
      "idx is 1217\n",
      "4356\n",
      "idx is 1218\n",
      "4147\n",
      "idx is 1219\n",
      "2885\n",
      "idx is 1220\n",
      "4986\n",
      "idx is 1221\n",
      "3214\n",
      "idx is 1222\n",
      "3576\n",
      "idx is 1223\n",
      "4220\n",
      "idx is 1224\n",
      "2671\n",
      "idx is 1225\n",
      "2539\n",
      "idx is 1226\n",
      "3273\n",
      "idx is 1227\n",
      "3703\n",
      "idx is 1228\n",
      "4984\n",
      "idx is 1229\n",
      "3651\n",
      "idx is 1230\n",
      "3427\n",
      "idx is 1231\n",
      "3899\n",
      "idx is 1232\n",
      "4717\n",
      "idx is 1233\n",
      "2923\n",
      "idx is 1234\n",
      "3111\n",
      "idx is 1235\n",
      "4547\n",
      "idx is 1236\n",
      "4059\n",
      "idx is 1237\n",
      "2680\n",
      "idx is 1238\n",
      "3884\n",
      "idx is 1239\n",
      "4717\n",
      "idx is 1240\n",
      "3023\n",
      "idx is 1241\n",
      "2519\n",
      "idx is 1242\n",
      "3776\n",
      "idx is 1243\n",
      "4957\n",
      "idx is 1244\n",
      "2581\n",
      "idx is 1245\n",
      "3840\n",
      "idx is 1246\n",
      "3235\n",
      "idx is 1247\n",
      "4423\n",
      "idx is 1248\n",
      "4785\n",
      "idx is 1249\n",
      "3718\n",
      "idx is 1250\n",
      "3752\n",
      "idx is 1251\n",
      "2869\n",
      "idx is 1252\n",
      "4243\n",
      "idx is 1253\n",
      "4376\n",
      "idx is 1254\n",
      "4667\n",
      "idx is 1255\n",
      "4358\n",
      "idx is 1256\n",
      "4477\n",
      "idx is 1257\n",
      "4752\n",
      "idx is 1258\n",
      "4535\n",
      "idx is 1259\n",
      "3420\n",
      "idx is 1260\n",
      "4170\n",
      "idx is 1261\n",
      "2991\n",
      "idx is 1262\n",
      "4444\n",
      "idx is 1263\n",
      "4632\n",
      "idx is 1264\n",
      "4000\n",
      "idx is 1265\n",
      "2711\n",
      "idx is 1266\n",
      "3445\n",
      "idx is 1267\n",
      "3584\n",
      "idx is 1268\n",
      "4206\n",
      "idx is 1269\n",
      "2561\n",
      "idx is 1270\n",
      "4550\n",
      "idx is 1271\n",
      "2671\n",
      "idx is 1272\n",
      "3220\n",
      "idx is 1273\n",
      "3695\n",
      "idx is 1274\n",
      "4922\n",
      "idx is 1275\n",
      "4600\n",
      "idx is 1276\n",
      "4094\n",
      "idx is 1277\n",
      "4813\n",
      "idx is 1278\n",
      "4355\n",
      "idx is 1279\n",
      "4447\n",
      "idx is 1280\n",
      "3652\n",
      "idx is 1281\n",
      "3577\n",
      "idx is 1282\n",
      "3891\n",
      "idx is 1283\n",
      "3547\n",
      "idx is 1284\n",
      "4562\n",
      "idx is 1285\n",
      "4201\n",
      "idx is 1286\n",
      "2713\n",
      "idx is 1287\n",
      "3725\n",
      "idx is 1288\n",
      "4137\n",
      "idx is 1289\n",
      "2687\n",
      "idx is 1290\n",
      "2980\n",
      "idx is 1291\n",
      "2524\n",
      "idx is 1292\n",
      "4387\n",
      "idx is 1293\n",
      "4219\n",
      "idx is 1294\n",
      "3005\n",
      "idx is 1295\n",
      "3199\n",
      "idx is 1296\n",
      "4520\n",
      "idx is 1297\n",
      "3691\n",
      "idx is 1298\n",
      "4561\n",
      "idx is 1299\n",
      "3541\n",
      "idx is 1300\n",
      "3717\n",
      "idx is 1301\n",
      "2617\n",
      "idx is 1302\n",
      "4198\n",
      "idx is 1303\n",
      "4710\n",
      "idx is 1304\n",
      "2695\n",
      "idx is 1305\n",
      "2834\n",
      "idx is 1306\n",
      "4891\n",
      "idx is 1307\n",
      "2704\n",
      "idx is 1308\n",
      "2996\n",
      "idx is 1309\n",
      "2609\n",
      "idx is 1310\n",
      "4947\n",
      "idx is 1311\n",
      "4807\n",
      "idx is 1312\n",
      "3256\n",
      "idx is 1313\n",
      "3664\n",
      "idx is 1314\n",
      "4369\n",
      "idx is 1315\n",
      "4077\n",
      "idx is 1316\n",
      "3926\n",
      "idx is 1317\n",
      "3028\n",
      "idx is 1318\n",
      "2581\n",
      "idx is 1319\n",
      "4533\n",
      "idx is 1320\n",
      "2638\n",
      "idx is 1321\n",
      "4371\n",
      "idx is 1322\n",
      "4535\n",
      "idx is 1323\n",
      "3114\n",
      "idx is 1324\n",
      "4662\n",
      "idx is 1325\n",
      "3167\n",
      "idx is 1326\n",
      "2967\n",
      "idx is 1327\n",
      "3505\n",
      "idx is 1328\n",
      "3790\n",
      "idx is 1329\n",
      "3225\n",
      "idx is 1330\n",
      "3009\n",
      "idx is 1331\n",
      "4648\n",
      "idx is 1332\n",
      "4367\n",
      "idx is 1333\n",
      "4877\n",
      "idx is 1334\n",
      "4473\n",
      "idx is 1335\n",
      "4796\n",
      "idx is 1336\n",
      "4603\n",
      "idx is 1337\n",
      "3667\n",
      "idx is 1338\n",
      "3798\n",
      "idx is 1339\n",
      "4714\n",
      "idx is 1340\n",
      "3904\n",
      "idx is 1341\n",
      "3402\n",
      "idx is 1342\n",
      "3692\n",
      "idx is 1343\n",
      "3844\n",
      "idx is 1344\n",
      "2950\n",
      "idx is 1345\n",
      "2811\n",
      "idx is 1346\n",
      "3308\n",
      "idx is 1347\n",
      "4214\n",
      "idx is 1348\n",
      "3401\n",
      "idx is 1349\n",
      "2830\n",
      "idx is 1350\n",
      "4599\n",
      "idx is 1351\n",
      "3094\n",
      "idx is 1352\n",
      "4850\n",
      "idx is 1353\n",
      "4020\n",
      "idx is 1354\n",
      "4563\n",
      "idx is 1355\n",
      "3505\n",
      "idx is 1356\n",
      "4961\n",
      "idx is 1357\n",
      "3267\n",
      "idx is 1358\n",
      "4503\n",
      "idx is 1359\n",
      "2953\n",
      "idx is 1360\n",
      "4147\n",
      "idx is 1361\n",
      "2725\n",
      "idx is 1362\n",
      "3017\n",
      "idx is 1363\n",
      "2923\n",
      "idx is 1364\n",
      "2716\n",
      "idx is 1365\n",
      "4510\n",
      "idx is 1366\n",
      "3633\n",
      "idx is 1367\n",
      "2620\n",
      "idx is 1368\n",
      "4848\n",
      "idx is 1369\n",
      "4689\n",
      "idx is 1370\n",
      "4890\n",
      "idx is 1371\n",
      "4007\n",
      "idx is 1372\n",
      "3831\n",
      "idx is 1373\n",
      "2935\n",
      "idx is 1374\n",
      "3357\n",
      "idx is 1375\n",
      "4801\n",
      "idx is 1376\n",
      "3816\n",
      "idx is 1377\n",
      "4158\n",
      "idx is 1378\n",
      "4210\n",
      "idx is 1379\n",
      "3198\n",
      "idx is 1380\n",
      "2683\n",
      "idx is 1381\n",
      "2786\n",
      "idx is 1382\n",
      "4176\n",
      "idx is 1383\n",
      "4489\n",
      "idx is 1384\n",
      "3828\n",
      "idx is 1385\n",
      "4676\n",
      "idx is 1386\n",
      "4561\n",
      "idx is 1387\n",
      "2685\n",
      "idx is 1388\n",
      "2587\n",
      "idx is 1389\n",
      "4749\n",
      "idx is 1390\n",
      "4773\n",
      "idx is 1391\n",
      "4095\n",
      "idx is 1392\n",
      "4241\n",
      "idx is 1393\n",
      "4232\n",
      "idx is 1394\n",
      "3517\n",
      "idx is 1395\n",
      "2697\n",
      "idx is 1396\n",
      "4809\n",
      "idx is 1397\n",
      "3410\n",
      "idx is 1398\n",
      "4254\n",
      "idx is 1399\n",
      "3263\n",
      "idx is 1400\n",
      "4875\n",
      "idx is 1401\n",
      "3458\n",
      "idx is 1402\n",
      "3857\n",
      "idx is 1403\n",
      "3283\n",
      "idx is 1404\n",
      "3073\n",
      "idx is 1405\n",
      "3560\n",
      "idx is 1406\n",
      "4078\n",
      "idx is 1407\n",
      "3302\n",
      "idx is 1408\n",
      "3129\n",
      "idx is 1409\n",
      "2999\n",
      "idx is 1410\n",
      "2809\n",
      "idx is 1411\n",
      "3876\n",
      "idx is 1412\n",
      "4100\n",
      "idx is 1413\n",
      "4130\n",
      "idx is 1414\n",
      "3996\n",
      "idx is 1415\n",
      "3321\n",
      "idx is 1416\n",
      "3648\n",
      "idx is 1417\n",
      "2678\n",
      "idx is 1418\n",
      "3382\n",
      "idx is 1419\n",
      "4742\n",
      "idx is 1420\n",
      "2834\n",
      "idx is 1421\n",
      "4827\n",
      "idx is 1422\n",
      "3293\n",
      "idx is 1423\n",
      "3148\n",
      "idx is 1424\n",
      "2977\n",
      "idx is 1425\n",
      "2529\n",
      "idx is 1426\n",
      "3101\n",
      "idx is 1427\n",
      "4272\n",
      "idx is 1428\n",
      "4346\n",
      "idx is 1429\n",
      "4718\n",
      "idx is 1430\n",
      "2833\n",
      "idx is 1431\n",
      "3029\n",
      "idx is 1432\n",
      "3127\n",
      "idx is 1433\n",
      "4014\n",
      "idx is 1434\n",
      "3772\n",
      "idx is 1435\n",
      "3228\n",
      "idx is 1436\n",
      "4704\n",
      "idx is 1437\n",
      "3693\n",
      "idx is 1438\n",
      "3375\n",
      "idx is 1439\n",
      "2571\n",
      "idx is 1440\n",
      "3854\n",
      "idx is 1441\n",
      "4021\n",
      "idx is 1442\n",
      "4895\n",
      "idx is 1443\n",
      "3216\n",
      "idx is 1444\n",
      "3034\n",
      "idx is 1445\n",
      "3462\n",
      "idx is 1446\n",
      "3847\n",
      "idx is 1447\n",
      "2872\n",
      "idx is 1448\n",
      "4537\n",
      "idx is 1449\n",
      "3039\n",
      "idx is 1450\n",
      "3756\n",
      "idx is 1451\n",
      "4588\n",
      "idx is 1452\n",
      "3661\n",
      "idx is 1453\n",
      "2724\n",
      "idx is 1454\n",
      "3072\n",
      "idx is 1455\n",
      "2542\n",
      "idx is 1456\n",
      "2542\n",
      "idx is 1457\n",
      "2642\n",
      "idx is 1458\n",
      "2589\n",
      "idx is 1459\n",
      "2819\n",
      "idx is 1460\n",
      "4720\n",
      "idx is 1461\n",
      "3393\n",
      "idx is 1462\n",
      "4630\n",
      "idx is 1463\n",
      "3086\n",
      "idx is 1464\n",
      "4852\n",
      "idx is 1465\n",
      "3416\n",
      "idx is 1466\n",
      "3014\n",
      "idx is 1467\n",
      "3166\n",
      "idx is 1468\n",
      "3507\n",
      "idx is 1469\n",
      "2591\n",
      "idx is 1470\n",
      "3937\n",
      "idx is 1471\n",
      "3147\n",
      "idx is 1472\n",
      "4682\n",
      "idx is 1473\n",
      "2612\n",
      "idx is 1474\n",
      "3487\n",
      "idx is 1475\n",
      "2521\n",
      "idx is 1476\n",
      "4202\n",
      "idx is 1477\n",
      "3726\n",
      "idx is 1478\n",
      "2726\n",
      "idx is 1479\n",
      "3820\n",
      "idx is 1480\n",
      "3403\n",
      "idx is 1481\n",
      "4932\n",
      "idx is 1482\n",
      "4331\n",
      "idx is 1483\n",
      "3655\n",
      "idx is 1484\n",
      "4231\n",
      "idx is 1485\n",
      "4977\n",
      "idx is 1486\n",
      "3202\n",
      "idx is 1487\n",
      "4277\n",
      "idx is 1488\n",
      "4657\n",
      "idx is 1489\n",
      "3134\n",
      "idx is 1490\n",
      "3374\n",
      "idx is 1491\n",
      "2909\n",
      "idx is 1492\n",
      "3635\n",
      "idx is 1493\n",
      "2812\n",
      "idx is 1494\n",
      "4082\n",
      "idx is 1495\n",
      "4311\n",
      "idx is 1496\n",
      "4036\n",
      "idx is 1497\n",
      "4328\n",
      "idx is 1498\n",
      "2991\n",
      "idx is 1499\n",
      "3759\n",
      "idx is 1500\n",
      "3043\n",
      "idx is 1501\n",
      "2971\n",
      "idx is 1502\n",
      "4178\n",
      "idx is 1503\n",
      "2865\n",
      "idx is 1504\n",
      "4954\n",
      "idx is 1505\n",
      "4855\n",
      "idx is 1506\n",
      "3239\n",
      "idx is 1507\n",
      "4423\n",
      "idx is 1508\n",
      "4674\n",
      "idx is 1509\n",
      "3647\n",
      "idx is 1510\n",
      "2814\n",
      "idx is 1511\n",
      "4420\n",
      "idx is 1512\n",
      "2843\n",
      "idx is 1513\n",
      "3744\n",
      "idx is 1514\n",
      "2703\n",
      "idx is 1515\n",
      "4636\n",
      "idx is 1516\n",
      "3149\n",
      "idx is 1517\n",
      "2634\n",
      "idx is 1518\n",
      "4785\n",
      "idx is 1519\n",
      "2777\n",
      "idx is 1520\n",
      "3596\n",
      "idx is 1521\n",
      "2812\n",
      "idx is 1522\n",
      "3846\n",
      "idx is 1523\n",
      "3074\n",
      "idx is 1524\n",
      "4715\n",
      "idx is 1525\n",
      "4036\n",
      "idx is 1526\n",
      "3989\n",
      "idx is 1527\n",
      "3703\n",
      "idx is 1528\n",
      "3090\n",
      "idx is 1529\n",
      "3352\n",
      "idx is 1530\n",
      "4402\n",
      "idx is 1531\n",
      "4404\n",
      "idx is 1532\n",
      "3375\n",
      "idx is 1533\n",
      "4362\n",
      "idx is 1534\n",
      "4358\n",
      "idx is 1535\n",
      "2629\n",
      "idx is 1536\n",
      "3480\n",
      "idx is 1537\n",
      "3893\n",
      "idx is 1538\n",
      "4327\n",
      "idx is 1539\n",
      "4312\n",
      "idx is 1540\n",
      "2693\n",
      "idx is 1541\n",
      "4879\n",
      "idx is 1542\n",
      "3740\n",
      "idx is 1543\n",
      "4093\n",
      "idx is 1544\n",
      "4924\n",
      "idx is 1545\n",
      "3783\n",
      "idx is 1546\n",
      "4709\n",
      "idx is 1547\n",
      "3913\n",
      "idx is 1548\n",
      "3521\n",
      "idx is 1549\n",
      "4134\n",
      "idx is 1550\n",
      "3169\n",
      "idx is 1551\n",
      "3102\n",
      "idx is 1552\n",
      "3501\n",
      "idx is 1553\n",
      "3817\n",
      "idx is 1554\n",
      "4114\n",
      "idx is 1555\n",
      "3572\n",
      "idx is 1556\n",
      "3542\n",
      "idx is 1557\n",
      "3405\n",
      "idx is 1558\n",
      "3375\n",
      "idx is 1559\n",
      "4647\n",
      "idx is 1560\n",
      "4674\n",
      "idx is 1561\n",
      "4996\n",
      "idx is 1562\n",
      "2811\n",
      "idx is 1563\n",
      "4698\n",
      "idx is 1564\n",
      "2849\n",
      "idx is 1565\n",
      "3143\n",
      "idx is 1566\n",
      "3543\n",
      "idx is 1567\n",
      "4530\n",
      "idx is 1568\n",
      "3071\n",
      "idx is 1569\n",
      "4553\n",
      "idx is 1570\n",
      "4850\n",
      "idx is 1571\n",
      "3416\n",
      "idx is 1572\n",
      "3914\n",
      "idx is 1573\n",
      "3542\n",
      "idx is 1574\n",
      "4403\n",
      "idx is 1575\n",
      "2892\n",
      "idx is 1576\n",
      "4228\n",
      "idx is 1577\n",
      "3331\n",
      "idx is 1578\n",
      "2813\n",
      "idx is 1579\n",
      "3028\n",
      "idx is 1580\n",
      "2937\n",
      "idx is 1581\n",
      "3540\n",
      "idx is 1582\n",
      "3732\n",
      "idx is 1583\n",
      "4286\n",
      "idx is 1584\n",
      "3970\n",
      "idx is 1585\n",
      "3544\n",
      "idx is 1586\n",
      "3640\n",
      "idx is 1587\n",
      "3648\n",
      "idx is 1588\n",
      "3828\n",
      "idx is 1589\n",
      "2744\n",
      "idx is 1590\n",
      "3812\n",
      "idx is 1591\n",
      "4141\n",
      "idx is 1592\n",
      "3835\n",
      "idx is 1593\n",
      "2888\n",
      "idx is 1594\n",
      "4249\n",
      "idx is 1595\n",
      "2898\n",
      "idx is 1596\n",
      "4818\n",
      "idx is 1597\n",
      "3817\n",
      "idx is 1598\n",
      "3006\n",
      "idx is 1599\n",
      "4026\n",
      "idx is 1600\n",
      "3910\n",
      "idx is 1601\n",
      "4758\n",
      "idx is 1602\n",
      "3060\n",
      "idx is 1603\n",
      "2920\n",
      "idx is 1604\n",
      "3260\n",
      "idx is 1605\n",
      "4691\n",
      "idx is 1606\n",
      "3408\n",
      "idx is 1607\n",
      "4706\n",
      "idx is 1608\n",
      "3279\n",
      "idx is 1609\n",
      "4461\n",
      "idx is 1610\n",
      "3348\n",
      "idx is 1611\n",
      "4212\n",
      "idx is 1612\n",
      "3848\n",
      "idx is 1613\n",
      "3681\n",
      "idx is 1614\n",
      "2551\n",
      "idx is 1615\n",
      "4092\n",
      "idx is 1616\n",
      "3024\n",
      "idx is 1617\n",
      "2694\n",
      "idx is 1618\n",
      "3078\n",
      "idx is 1619\n",
      "2677\n",
      "idx is 1620\n",
      "4399\n",
      "idx is 1621\n",
      "4410\n",
      "idx is 1622\n",
      "4012\n",
      "idx is 1623\n",
      "4857\n",
      "idx is 1624\n",
      "2825\n",
      "idx is 1625\n",
      "3337\n",
      "idx is 1626\n",
      "3729\n",
      "idx is 1627\n",
      "3812\n",
      "idx is 1628\n",
      "2744\n",
      "idx is 1629\n",
      "3374\n",
      "idx is 1630\n",
      "3185\n",
      "idx is 1631\n",
      "4453\n",
      "idx is 1632\n",
      "3138\n",
      "idx is 1633\n",
      "3081\n",
      "idx is 1634\n",
      "4148\n",
      "idx is 1635\n",
      "3478\n",
      "idx is 1636\n",
      "4219\n",
      "idx is 1637\n",
      "3764\n",
      "idx is 1638\n",
      "2666\n",
      "idx is 1639\n",
      "2918\n",
      "idx is 1640\n",
      "3189\n",
      "idx is 1641\n",
      "3598\n",
      "idx is 1642\n",
      "4512\n",
      "idx is 1643\n",
      "3080\n",
      "idx is 1644\n",
      "4534\n",
      "idx is 1645\n",
      "4777\n",
      "idx is 1646\n",
      "4272\n",
      "idx is 1647\n",
      "4458\n",
      "idx is 1648\n",
      "3372\n",
      "idx is 1649\n",
      "4643\n",
      "idx is 1650\n",
      "3260\n",
      "idx is 1651\n",
      "4186\n",
      "idx is 1652\n",
      "4806\n",
      "idx is 1653\n",
      "3479\n",
      "idx is 1654\n",
      "3382\n",
      "idx is 1655\n",
      "3755\n",
      "idx is 1656\n",
      "4815\n",
      "idx is 1657\n",
      "4101\n",
      "idx is 1658\n",
      "4768\n",
      "idx is 1659\n",
      "2616\n",
      "idx is 1660\n",
      "3830\n",
      "idx is 1661\n",
      "4042\n",
      "idx is 1662\n",
      "4193\n",
      "idx is 1663\n",
      "4375\n",
      "idx is 1664\n",
      "3214\n",
      "idx is 1665\n",
      "4114\n",
      "idx is 1666\n",
      "4001\n",
      "idx is 1667\n",
      "3714\n",
      "idx is 1668\n",
      "4526\n",
      "idx is 1669\n",
      "3562\n",
      "idx is 1670\n",
      "3271\n",
      "idx is 1671\n",
      "2827\n",
      "idx is 1672\n",
      "4266\n",
      "idx is 1673\n",
      "4731\n",
      "idx is 1674\n",
      "3077\n",
      "idx is 1675\n",
      "3490\n",
      "idx is 1676\n",
      "4369\n",
      "idx is 1677\n",
      "4040\n",
      "idx is 1678\n",
      "4278\n",
      "idx is 1679\n",
      "3934\n",
      "idx is 1680\n",
      "4046\n",
      "idx is 1681\n",
      "2905\n",
      "idx is 1682\n",
      "4825\n",
      "idx is 1683\n",
      "4459\n",
      "idx is 1684\n",
      "4914\n",
      "idx is 1685\n",
      "2882\n",
      "idx is 1686\n",
      "3380\n",
      "idx is 1687\n",
      "2935\n",
      "idx is 1688\n",
      "3133\n",
      "idx is 1689\n",
      "2607\n",
      "idx is 1690\n",
      "3656\n",
      "idx is 1691\n",
      "2790\n",
      "idx is 1692\n",
      "2823\n",
      "idx is 1693\n",
      "3248\n",
      "idx is 1694\n",
      "2968\n",
      "idx is 1695\n",
      "2914\n",
      "idx is 1696\n",
      "3939\n",
      "idx is 1697\n",
      "4652\n",
      "idx is 1698\n",
      "3220\n",
      "idx is 1699\n",
      "2733\n",
      "idx is 1700\n",
      "3010\n",
      "idx is 1701\n",
      "3098\n",
      "idx is 1702\n",
      "3788\n",
      "idx is 1703\n",
      "4414\n",
      "idx is 1704\n",
      "4299\n",
      "idx is 1705\n",
      "3428\n",
      "idx is 1706\n",
      "4400\n",
      "idx is 1707\n",
      "2841\n",
      "idx is 1708\n",
      "3910\n",
      "idx is 1709\n",
      "4345\n",
      "idx is 1710\n",
      "2556\n",
      "idx is 1711\n",
      "2911\n",
      "idx is 1712\n",
      "4002\n",
      "idx is 1713\n",
      "2774\n",
      "idx is 1714\n",
      "2895\n",
      "idx is 1715\n",
      "4427\n",
      "idx is 1716\n",
      "3883\n",
      "idx is 1717\n",
      "3306\n",
      "idx is 1718\n",
      "3390\n",
      "idx is 1719\n",
      "4651\n",
      "idx is 1720\n",
      "4357\n",
      "idx is 1721\n",
      "4414\n",
      "idx is 1722\n",
      "4467\n",
      "idx is 1723\n",
      "4814\n",
      "idx is 1724\n",
      "3396\n",
      "idx is 1725\n",
      "4816\n",
      "idx is 1726\n",
      "3101\n",
      "idx is 1727\n",
      "2809\n",
      "idx is 1728\n",
      "2841\n",
      "idx is 1729\n",
      "4115\n",
      "idx is 1730\n",
      "4336\n",
      "idx is 1731\n",
      "2863\n",
      "idx is 1732\n",
      "4988\n",
      "idx is 1733\n",
      "3810\n",
      "idx is 1734\n",
      "4822\n",
      "idx is 1735\n",
      "2584\n",
      "idx is 1736\n",
      "4672\n",
      "idx is 1737\n",
      "4337\n",
      "idx is 1738\n",
      "3054\n",
      "idx is 1739\n",
      "4659\n",
      "idx is 1740\n",
      "3221\n",
      "idx is 1741\n",
      "3583\n",
      "idx is 1742\n",
      "3563\n",
      "idx is 1743\n",
      "3838\n",
      "idx is 1744\n",
      "4948\n",
      "idx is 1745\n",
      "2696\n",
      "idx is 1746\n",
      "3054\n",
      "idx is 1747\n",
      "2727\n",
      "idx is 1748\n",
      "3546\n",
      "idx is 1749\n",
      "4912\n",
      "idx is 1750\n",
      "3652\n",
      "idx is 1751\n",
      "3682\n",
      "idx is 1752\n",
      "4597\n",
      "idx is 1753\n",
      "3272\n",
      "idx is 1754\n",
      "2624\n",
      "idx is 1755\n",
      "3637\n",
      "idx is 1756\n",
      "3498\n",
      "idx is 1757\n",
      "3571\n",
      "idx is 1758\n",
      "4827\n",
      "idx is 1759\n",
      "4954\n",
      "idx is 1760\n",
      "4204\n",
      "idx is 1761\n",
      "3705\n",
      "idx is 1762\n",
      "3301\n",
      "idx is 1763\n",
      "3251\n",
      "idx is 1764\n",
      "3941\n",
      "idx is 1765\n",
      "3143\n",
      "idx is 1766\n",
      "4245\n",
      "idx is 1767\n",
      "3237\n",
      "idx is 1768\n",
      "4096\n",
      "idx is 1769\n",
      "4046\n",
      "idx is 1770\n",
      "4359\n",
      "idx is 1771\n",
      "2779\n",
      "idx is 1772\n",
      "4490\n",
      "idx is 1773\n",
      "3764\n",
      "idx is 1774\n",
      "4593\n",
      "idx is 1775\n",
      "2555\n",
      "idx is 1776\n",
      "2916\n",
      "idx is 1777\n",
      "3788\n",
      "idx is 1778\n",
      "3655\n",
      "idx is 1779\n",
      "4251\n",
      "idx is 1780\n",
      "4543\n",
      "idx is 1781\n",
      "3234\n",
      "idx is 1782\n",
      "3721\n",
      "idx is 1783\n",
      "3465\n",
      "idx is 1784\n",
      "4697\n",
      "idx is 1785\n",
      "2607\n",
      "idx is 1786\n",
      "4076\n",
      "idx is 1787\n",
      "3204\n",
      "idx is 1788\n",
      "3361\n",
      "idx is 1789\n",
      "4635\n",
      "idx is 1790\n",
      "2543\n",
      "idx is 1791\n",
      "4627\n",
      "idx is 1792\n",
      "2761\n",
      "idx is 1793\n",
      "2832\n",
      "idx is 1794\n",
      "2734\n",
      "idx is 1795\n",
      "3489\n",
      "idx is 1796\n",
      "4722\n",
      "idx is 1797\n",
      "4797\n",
      "idx is 1798\n",
      "4199\n",
      "idx is 1799\n",
      "3835\n",
      "idx is 1800\n",
      "2900\n",
      "idx is 1801\n",
      "4364\n",
      "idx is 1802\n",
      "4720\n",
      "idx is 1803\n",
      "3849\n",
      "idx is 1804\n",
      "3744\n",
      "idx is 1805\n",
      "3110\n",
      "idx is 1806\n",
      "3991\n",
      "idx is 1807\n",
      "2952\n",
      "idx is 1808\n",
      "4134\n",
      "idx is 1809\n",
      "3229\n",
      "idx is 1810\n",
      "2764\n",
      "idx is 1811\n",
      "2658\n",
      "idx is 1812\n",
      "4251\n",
      "idx is 1813\n",
      "4822\n",
      "idx is 1814\n",
      "4281\n",
      "idx is 1815\n",
      "4041\n",
      "idx is 1816\n",
      "4681\n",
      "idx is 1817\n",
      "3022\n",
      "idx is 1818\n",
      "3731\n",
      "idx is 1819\n",
      "4917\n",
      "idx is 1820\n",
      "2537\n",
      "idx is 1821\n",
      "4444\n",
      "idx is 1822\n",
      "3122\n",
      "idx is 1823\n",
      "4430\n",
      "idx is 1824\n",
      "3692\n",
      "idx is 1825\n",
      "4434\n",
      "idx is 1826\n",
      "4606\n",
      "idx is 1827\n",
      "4213\n",
      "idx is 1828\n",
      "4819\n",
      "idx is 1829\n",
      "3674\n",
      "idx is 1830\n",
      "2802\n",
      "idx is 1831\n",
      "4277\n",
      "idx is 1832\n",
      "4493\n",
      "idx is 1833\n",
      "4199\n",
      "idx is 1834\n",
      "2959\n",
      "idx is 1835\n",
      "3562\n",
      "idx is 1836\n",
      "2819\n",
      "idx is 1837\n",
      "3739\n",
      "idx is 1838\n",
      "4530\n",
      "idx is 1839\n",
      "4613\n",
      "idx is 1840\n",
      "4853\n",
      "idx is 1841\n",
      "4408\n",
      "idx is 1842\n",
      "4724\n",
      "idx is 1843\n",
      "3141\n",
      "idx is 1844\n",
      "4058\n",
      "idx is 1845\n",
      "4585\n",
      "idx is 1846\n",
      "3174\n",
      "idx is 1847\n",
      "4274\n",
      "idx is 1848\n",
      "2601\n",
      "idx is 1849\n",
      "2947\n",
      "idx is 1850\n",
      "2812\n",
      "idx is 1851\n",
      "3368\n",
      "idx is 1852\n",
      "4488\n",
      "idx is 1853\n",
      "3401\n",
      "idx is 1854\n",
      "4563\n",
      "idx is 1855\n",
      "4148\n",
      "idx is 1856\n",
      "2863\n",
      "idx is 1857\n",
      "4258\n",
      "idx is 1858\n",
      "3969\n",
      "idx is 1859\n",
      "3241\n",
      "idx is 1860\n",
      "4375\n",
      "idx is 1861\n",
      "3877\n",
      "idx is 1862\n",
      "3540\n",
      "idx is 1863\n",
      "4968\n",
      "idx is 1864\n",
      "3591\n",
      "idx is 1865\n",
      "4389\n",
      "idx is 1866\n",
      "3379\n",
      "idx is 1867\n",
      "3896\n",
      "idx is 1868\n",
      "2571\n",
      "idx is 1869\n",
      "3477\n",
      "idx is 1870\n",
      "4238\n",
      "idx is 1871\n",
      "3903\n",
      "idx is 1872\n",
      "4243\n",
      "idx is 1873\n",
      "2773\n",
      "idx is 1874\n",
      "3111\n",
      "idx is 1875\n",
      "3444\n",
      "idx is 1876\n",
      "2811\n",
      "idx is 1877\n",
      "4281\n",
      "idx is 1878\n",
      "3125\n",
      "idx is 1879\n",
      "4379\n",
      "idx is 1880\n",
      "3502\n",
      "idx is 1881\n",
      "2625\n",
      "idx is 1882\n",
      "2910\n",
      "idx is 1883\n",
      "4157\n",
      "idx is 1884\n",
      "3654\n",
      "idx is 1885\n",
      "3813\n",
      "idx is 1886\n",
      "2935\n",
      "idx is 1887\n",
      "4965\n",
      "idx is 1888\n",
      "3697\n",
      "idx is 1889\n",
      "2972\n",
      "idx is 1890\n",
      "3122\n",
      "idx is 1891\n",
      "4313\n",
      "idx is 1892\n",
      "4425\n",
      "idx is 1893\n",
      "4558\n",
      "idx is 1894\n",
      "4086\n",
      "idx is 1895\n",
      "4017\n",
      "idx is 1896\n",
      "4367\n",
      "idx is 1897\n",
      "3153\n",
      "idx is 1898\n",
      "2579\n",
      "idx is 1899\n",
      "4437\n",
      "idx is 1900\n",
      "4853\n",
      "idx is 1901\n",
      "2972\n",
      "idx is 1902\n",
      "3100\n",
      "idx is 1903\n",
      "4296\n",
      "idx is 1904\n",
      "3537\n",
      "idx is 1905\n",
      "3295\n",
      "idx is 1906\n",
      "3699\n",
      "idx is 1907\n",
      "2543\n",
      "idx is 1908\n",
      "3869\n",
      "idx is 1909\n",
      "4995\n",
      "idx is 1910\n",
      "3088\n",
      "idx is 1911\n",
      "3884\n",
      "idx is 1912\n",
      "3887\n",
      "idx is 1913\n",
      "3130\n",
      "idx is 1914\n",
      "2585\n",
      "idx is 1915\n",
      "4623\n",
      "idx is 1916\n",
      "2933\n",
      "idx is 1917\n",
      "3858\n",
      "idx is 1918\n",
      "3381\n",
      "idx is 1919\n",
      "2732\n",
      "idx is 1920\n",
      "4845\n",
      "idx is 1921\n",
      "3155\n",
      "idx is 1922\n",
      "2929\n",
      "idx is 1923\n",
      "4362\n",
      "idx is 1924\n",
      "3919\n",
      "idx is 1925\n",
      "3605\n",
      "idx is 1926\n",
      "4982\n",
      "idx is 1927\n",
      "3918\n",
      "idx is 1928\n",
      "3594\n",
      "idx is 1929\n",
      "2808\n",
      "idx is 1930\n",
      "4862\n",
      "idx is 1931\n",
      "4632\n",
      "idx is 1932\n",
      "4238\n",
      "idx is 1933\n",
      "4348\n",
      "idx is 1934\n",
      "4794\n",
      "idx is 1935\n",
      "4682\n",
      "idx is 1936\n",
      "3933\n",
      "idx is 1937\n",
      "2634\n",
      "idx is 1938\n",
      "3359\n",
      "idx is 1939\n",
      "3478\n",
      "idx is 1940\n",
      "2698\n",
      "idx is 1941\n",
      "4730\n",
      "idx is 1942\n",
      "2779\n",
      "idx is 1943\n",
      "3746\n",
      "idx is 1944\n",
      "4826\n",
      "idx is 1945\n",
      "3311\n",
      "idx is 1946\n",
      "3528\n",
      "idx is 1947\n",
      "4680\n",
      "idx is 1948\n",
      "3518\n",
      "idx is 1949\n",
      "4715\n",
      "idx is 1950\n",
      "4458\n",
      "idx is 1951\n",
      "2582\n",
      "idx is 1952\n",
      "3925\n",
      "idx is 1953\n",
      "4905\n",
      "idx is 1954\n",
      "2809\n",
      "idx is 1955\n",
      "4528\n",
      "idx is 1956\n",
      "3648\n",
      "idx is 1957\n",
      "4452\n",
      "idx is 1958\n",
      "4130\n",
      "idx is 1959\n",
      "3595\n",
      "idx is 1960\n",
      "3685\n",
      "idx is 1961\n",
      "2735\n",
      "idx is 1962\n",
      "4835\n",
      "idx is 1963\n",
      "3243\n",
      "idx is 1964\n",
      "2688\n",
      "idx is 1965\n",
      "3294\n",
      "idx is 1966\n",
      "2623\n",
      "idx is 1967\n",
      "2967\n",
      "idx is 1968\n",
      "3480\n",
      "idx is 1969\n",
      "2599\n",
      "idx is 1970\n",
      "2899\n",
      "idx is 1971\n",
      "2886\n",
      "idx is 1972\n",
      "3067\n",
      "idx is 1973\n",
      "3283\n",
      "idx is 1974\n",
      "4640\n",
      "idx is 1975\n",
      "3154\n",
      "idx is 1976\n",
      "4117\n",
      "idx is 1977\n",
      "3336\n",
      "idx is 1978\n",
      "4753\n",
      "idx is 1979\n",
      "3959\n",
      "idx is 1980\n",
      "4412\n",
      "idx is 1981\n",
      "3478\n",
      "idx is 1982\n",
      "2841\n",
      "idx is 1983\n",
      "3089\n",
      "idx is 1984\n",
      "3062\n",
      "idx is 1985\n",
      "4299\n",
      "idx is 1986\n",
      "4375\n",
      "idx is 1987\n",
      "4820\n",
      "idx is 1988\n",
      "2619\n",
      "idx is 1989\n",
      "4582\n",
      "idx is 1990\n",
      "2759\n",
      "idx is 1991\n",
      "2985\n",
      "idx is 1992\n",
      "3190\n",
      "idx is 1993\n",
      "4460\n",
      "idx is 1994\n",
      "4611\n",
      "idx is 1995\n",
      "3294\n",
      "idx is 1996\n",
      "4760\n",
      "idx is 1997\n",
      "4994\n",
      "idx is 1998\n",
      "3400\n",
      "idx is 1999\n",
      "3797\n",
      "idx is 2000\n",
      "3062\n",
      "idx is 2001\n",
      "4447\n",
      "idx is 2002\n",
      "4008\n",
      "idx is 2003\n",
      "4171\n",
      "idx is 2004\n",
      "2725\n",
      "idx is 2005\n",
      "2605\n",
      "idx is 2006\n",
      "3964\n",
      "idx is 2007\n",
      "3349\n",
      "idx is 2008\n",
      "4943\n",
      "idx is 2009\n",
      "3317\n",
      "idx is 2010\n",
      "4837\n",
      "idx is 2011\n",
      "4582\n",
      "idx is 2012\n",
      "3501\n",
      "idx is 2013\n",
      "2869\n",
      "idx is 2014\n",
      "3801\n",
      "idx is 2015\n",
      "3919\n",
      "idx is 2016\n",
      "3388\n",
      "idx is 2017\n",
      "4953\n",
      "idx is 2018\n",
      "4350\n",
      "idx is 2019\n",
      "4874\n",
      "idx is 2020\n",
      "2802\n",
      "idx is 2021\n",
      "2876\n",
      "idx is 2022\n",
      "3222\n",
      "idx is 2023\n",
      "3675\n",
      "idx is 2024\n",
      "4718\n",
      "idx is 2025\n",
      "3560\n",
      "idx is 2026\n",
      "3878\n",
      "idx is 2027\n",
      "4964\n",
      "idx is 2028\n",
      "3769\n",
      "idx is 2029\n",
      "3772\n",
      "idx is 2030\n",
      "4345\n",
      "idx is 2031\n",
      "3954\n",
      "idx is 2032\n",
      "4344\n",
      "idx is 2033\n",
      "4127\n",
      "idx is 2034\n",
      "3573\n",
      "idx is 2035\n",
      "3150\n",
      "idx is 2036\n",
      "3629\n",
      "idx is 2037\n",
      "3687\n",
      "idx is 2038\n",
      "4528\n",
      "idx is 2039\n",
      "2716\n",
      "idx is 2040\n",
      "3975\n",
      "idx is 2041\n",
      "3804\n",
      "idx is 2042\n",
      "4400\n",
      "idx is 2043\n",
      "3904\n",
      "idx is 2044\n",
      "3157\n",
      "idx is 2045\n",
      "3388\n",
      "idx is 2046\n",
      "3285\n",
      "idx is 2047\n",
      "3591\n",
      "idx is 2048\n",
      "4264\n",
      "idx is 2049\n",
      "4730\n",
      "idx is 2050\n",
      "4984\n",
      "idx is 2051\n",
      "4098\n",
      "idx is 2052\n",
      "3325\n",
      "idx is 2053\n",
      "3115\n",
      "idx is 2054\n",
      "4016\n",
      "idx is 2055\n",
      "4022\n",
      "idx is 2056\n",
      "2954\n",
      "idx is 2057\n",
      "3287\n",
      "idx is 2058\n",
      "3349\n",
      "idx is 2059\n",
      "2674\n",
      "idx is 2060\n",
      "3803\n",
      "idx is 2061\n",
      "2790\n",
      "idx is 2062\n",
      "3778\n",
      "idx is 2063\n",
      "2865\n",
      "idx is 2064\n",
      "4037\n",
      "idx is 2065\n",
      "4682\n",
      "idx is 2066\n",
      "4369\n",
      "idx is 2067\n",
      "3472\n",
      "idx is 2068\n",
      "2983\n",
      "idx is 2069\n",
      "3475\n",
      "idx is 2070\n",
      "4585\n",
      "idx is 2071\n",
      "3126\n",
      "idx is 2072\n",
      "4086\n",
      "idx is 2073\n",
      "4976\n",
      "idx is 2074\n",
      "3543\n",
      "idx is 2075\n",
      "3453\n",
      "idx is 2076\n",
      "4897\n",
      "idx is 2077\n",
      "4585\n",
      "idx is 2078\n",
      "2989\n",
      "idx is 2079\n",
      "3319\n",
      "idx is 2080\n",
      "4634\n",
      "idx is 2081\n",
      "3715\n",
      "idx is 2082\n",
      "4508\n",
      "idx is 2083\n",
      "3732\n",
      "idx is 2084\n",
      "3226\n",
      "idx is 2085\n",
      "2706\n",
      "idx is 2086\n",
      "4406\n",
      "idx is 2087\n",
      "4345\n",
      "idx is 2088\n",
      "3808\n",
      "idx is 2089\n",
      "4355\n",
      "idx is 2090\n",
      "3123\n",
      "idx is 2091\n",
      "4871\n",
      "idx is 2092\n",
      "3311\n",
      "idx is 2093\n",
      "2581\n",
      "idx is 2094\n",
      "2526\n",
      "idx is 2095\n",
      "4698\n",
      "idx is 2096\n",
      "2518\n",
      "idx is 2097\n",
      "2521\n",
      "idx is 2098\n",
      "3003\n",
      "idx is 2099\n",
      "4016\n",
      "idx is 2100\n",
      "4648\n",
      "idx is 2101\n",
      "3513\n",
      "idx is 2102\n",
      "3342\n",
      "idx is 2103\n",
      "3743\n",
      "idx is 2104\n",
      "4171\n",
      "idx is 2105\n",
      "3365\n",
      "idx is 2106\n",
      "4437\n",
      "idx is 2107\n",
      "4282\n",
      "idx is 2108\n",
      "4215\n",
      "idx is 2109\n",
      "4398\n",
      "idx is 2110\n",
      "3233\n",
      "idx is 2111\n",
      "4405\n",
      "idx is 2112\n",
      "4832\n",
      "idx is 2113\n",
      "3281\n",
      "idx is 2114\n",
      "3014\n",
      "idx is 2115\n",
      "2979\n",
      "idx is 2116\n",
      "2889\n",
      "idx is 2117\n",
      "2585\n",
      "idx is 2118\n",
      "3082\n",
      "idx is 2119\n",
      "3295\n",
      "idx is 2120\n",
      "4801\n",
      "idx is 2121\n",
      "3880\n",
      "idx is 2122\n",
      "3635\n",
      "idx is 2123\n",
      "3039\n",
      "idx is 2124\n",
      "3494\n",
      "idx is 2125\n",
      "4073\n",
      "idx is 2126\n",
      "4554\n",
      "idx is 2127\n",
      "4457\n",
      "idx is 2128\n",
      "2539\n",
      "idx is 2129\n",
      "4770\n",
      "idx is 2130\n",
      "4989\n",
      "idx is 2131\n",
      "4220\n",
      "idx is 2132\n",
      "4009\n",
      "idx is 2133\n",
      "4053\n",
      "idx is 2134\n",
      "3138\n",
      "idx is 2135\n",
      "4612\n",
      "idx is 2136\n",
      "3995\n",
      "idx is 2137\n",
      "3423\n",
      "idx is 2138\n",
      "4690\n",
      "idx is 2139\n",
      "3497\n",
      "idx is 2140\n",
      "3682\n",
      "idx is 2141\n",
      "3317\n",
      "idx is 2142\n",
      "3090\n",
      "idx is 2143\n",
      "3629\n",
      "idx is 2144\n",
      "4770\n",
      "idx is 2145\n",
      "4471\n",
      "idx is 2146\n",
      "2814\n",
      "idx is 2147\n",
      "3097\n",
      "idx is 2148\n",
      "3468\n",
      "idx is 2149\n",
      "4019\n",
      "idx is 2150\n",
      "3831\n",
      "idx is 2151\n",
      "2520\n",
      "idx is 2152\n",
      "4930\n",
      "idx is 2153\n",
      "4904\n",
      "idx is 2154\n",
      "3726\n",
      "idx is 2155\n",
      "4322\n",
      "idx is 2156\n",
      "3110\n",
      "idx is 2157\n",
      "2680\n",
      "idx is 2158\n",
      "2688\n",
      "idx is 2159\n",
      "4436\n",
      "idx is 2160\n",
      "4504\n",
      "idx is 2161\n",
      "4428\n",
      "idx is 2162\n",
      "4971\n",
      "idx is 2163\n",
      "3673\n",
      "idx is 2164\n",
      "4024\n",
      "idx is 2165\n",
      "3688\n",
      "idx is 2166\n",
      "3375\n",
      "idx is 2167\n",
      "3193\n",
      "idx is 2168\n",
      "4494\n",
      "idx is 2169\n",
      "4827\n",
      "idx is 2170\n",
      "2852\n",
      "idx is 2171\n",
      "2720\n",
      "idx is 2172\n",
      "4124\n",
      "idx is 2173\n",
      "4405\n",
      "idx is 2174\n",
      "3175\n",
      "idx is 2175\n",
      "3004\n",
      "idx is 2176\n",
      "3839\n",
      "idx is 2177\n",
      "4895\n",
      "idx is 2178\n",
      "3678\n",
      "idx is 2179\n",
      "4815\n",
      "idx is 2180\n",
      "3666\n",
      "idx is 2181\n",
      "2776\n",
      "idx is 2182\n",
      "3226\n",
      "idx is 2183\n",
      "3233\n",
      "idx is 2184\n",
      "3896\n",
      "idx is 2185\n",
      "3913\n",
      "idx is 2186\n",
      "4325\n",
      "idx is 2187\n",
      "4578\n",
      "idx is 2188\n",
      "3046\n",
      "idx is 2189\n",
      "4801\n",
      "idx is 2190\n",
      "3580\n",
      "idx is 2191\n",
      "3894\n",
      "idx is 2192\n",
      "3009\n",
      "idx is 2193\n",
      "3226\n",
      "idx is 2194\n",
      "4441\n",
      "idx is 2195\n",
      "2954\n",
      "idx is 2196\n",
      "3223\n",
      "idx is 2197\n",
      "3587\n",
      "idx is 2198\n",
      "2775\n",
      "idx is 2199\n",
      "3267\n",
      "idx is 2200\n",
      "2578\n",
      "idx is 2201\n",
      "3294\n",
      "idx is 2202\n",
      "3459\n",
      "idx is 2203\n",
      "3223\n",
      "idx is 2204\n",
      "2599\n",
      "idx is 2205\n",
      "3553\n",
      "idx is 2206\n",
      "4997\n",
      "idx is 2207\n",
      "3818\n",
      "idx is 2208\n",
      "4895\n",
      "idx is 2209\n",
      "4054\n",
      "idx is 2210\n",
      "3964\n",
      "idx is 2211\n",
      "3945\n",
      "idx is 2212\n",
      "2576\n",
      "idx is 2213\n",
      "4592\n",
      "idx is 2214\n",
      "2542\n",
      "idx is 2215\n",
      "3239\n",
      "idx is 2216\n",
      "4939\n",
      "idx is 2217\n",
      "3181\n",
      "idx is 2218\n",
      "3903\n",
      "idx is 2219\n",
      "3833\n",
      "idx is 2220\n",
      "3322\n",
      "idx is 2221\n",
      "3606\n",
      "idx is 2222\n",
      "3390\n",
      "idx is 2223\n",
      "3435\n",
      "idx is 2224\n",
      "4117\n",
      "idx is 2225\n",
      "3130\n",
      "idx is 2226\n",
      "4102\n",
      "idx is 2227\n",
      "2936\n",
      "idx is 2228\n",
      "4961\n",
      "idx is 2229\n",
      "4104\n",
      "idx is 2230\n",
      "3054\n",
      "idx is 2231\n",
      "2755\n",
      "idx is 2232\n",
      "3658\n",
      "idx is 2233\n",
      "4932\n",
      "idx is 2234\n",
      "2974\n",
      "idx is 2235\n",
      "3628\n",
      "idx is 2236\n",
      "4886\n",
      "idx is 2237\n",
      "3993\n",
      "idx is 2238\n",
      "3594\n",
      "idx is 2239\n",
      "4951\n",
      "idx is 2240\n",
      "3149\n",
      "idx is 2241\n",
      "3903\n",
      "idx is 2242\n",
      "2977\n",
      "idx is 2243\n",
      "4610\n",
      "idx is 2244\n",
      "3278\n",
      "idx is 2245\n",
      "4912\n",
      "idx is 2246\n",
      "4761\n",
      "idx is 2247\n",
      "4149\n",
      "idx is 2248\n",
      "3337\n",
      "idx is 2249\n",
      "4558\n",
      "idx is 2250\n",
      "2699\n",
      "idx is 2251\n",
      "3282\n",
      "idx is 2252\n",
      "2776\n",
      "idx is 2253\n",
      "3850\n",
      "idx is 2254\n",
      "3674\n",
      "idx is 2255\n",
      "3714\n",
      "idx is 2256\n",
      "4712\n",
      "idx is 2257\n",
      "4457\n",
      "idx is 2258\n",
      "3999\n",
      "idx is 2259\n",
      "3643\n",
      "idx is 2260\n",
      "4497\n",
      "idx is 2261\n",
      "4574\n",
      "idx is 2262\n",
      "2962\n",
      "idx is 2263\n",
      "2552\n",
      "idx is 2264\n",
      "3653\n",
      "idx is 2265\n",
      "3878\n",
      "idx is 2266\n",
      "3288\n",
      "idx is 2267\n",
      "2986\n",
      "idx is 2268\n",
      "4728\n",
      "idx is 2269\n",
      "2854\n",
      "idx is 2270\n",
      "4410\n",
      "idx is 2271\n",
      "2689\n",
      "idx is 2272\n",
      "4074\n",
      "idx is 2273\n",
      "2844\n",
      "idx is 2274\n",
      "4259\n",
      "idx is 2275\n",
      "4230\n",
      "idx is 2276\n",
      "3724\n",
      "idx is 2277\n",
      "4866\n",
      "idx is 2278\n",
      "2505\n",
      "idx is 2279\n",
      "4286\n",
      "idx is 2280\n",
      "4112\n",
      "idx is 2281\n",
      "3545\n",
      "idx is 2282\n",
      "2639\n",
      "idx is 2283\n",
      "3854\n",
      "idx is 2284\n",
      "4718\n",
      "idx is 2285\n",
      "2843\n",
      "idx is 2286\n",
      "2572\n",
      "idx is 2287\n",
      "3572\n",
      "idx is 2288\n",
      "3216\n",
      "idx is 2289\n",
      "4418\n",
      "idx is 2290\n",
      "2584\n",
      "idx is 2291\n",
      "3038\n",
      "idx is 2292\n",
      "4089\n",
      "idx is 2293\n",
      "4600\n",
      "idx is 2294\n",
      "2514\n",
      "idx is 2295\n",
      "3961\n",
      "idx is 2296\n",
      "3017\n",
      "idx is 2297\n",
      "3110\n",
      "idx is 2298\n",
      "3779\n",
      "idx is 2299\n",
      "2775\n",
      "idx is 2300\n",
      "2759\n",
      "idx is 2301\n",
      "3112\n",
      "idx is 2302\n",
      "2673\n",
      "idx is 2303\n",
      "2938\n",
      "idx is 2304\n",
      "4005\n",
      "idx is 2305\n",
      "4849\n",
      "idx is 2306\n",
      "2509\n",
      "idx is 2307\n",
      "3173\n",
      "idx is 2308\n",
      "4854\n",
      "idx is 2309\n",
      "3145\n",
      "idx is 2310\n",
      "3481\n",
      "idx is 2311\n",
      "3856\n",
      "idx is 2312\n",
      "3080\n",
      "idx is 2313\n",
      "4100\n",
      "idx is 2314\n",
      "4708\n",
      "idx is 2315\n",
      "3082\n",
      "idx is 2316\n",
      "3408\n",
      "idx is 2317\n",
      "4333\n",
      "idx is 2318\n",
      "4389\n",
      "idx is 2319\n",
      "4416\n",
      "idx is 2320\n",
      "4483\n",
      "idx is 2321\n",
      "4681\n",
      "idx is 2322\n",
      "4846\n",
      "idx is 2323\n",
      "4748\n",
      "idx is 2324\n",
      "4241\n",
      "idx is 2325\n",
      "3727\n",
      "idx is 2326\n",
      "4571\n",
      "idx is 2327\n",
      "3462\n",
      "idx is 2328\n",
      "2681\n",
      "idx is 2329\n",
      "3120\n",
      "idx is 2330\n",
      "2771\n",
      "idx is 2331\n",
      "2612\n",
      "idx is 2332\n",
      "4497\n",
      "idx is 2333\n",
      "3292\n",
      "idx is 2334\n",
      "4970\n",
      "idx is 2335\n",
      "2872\n",
      "idx is 2336\n",
      "4424\n",
      "idx is 2337\n",
      "3227\n",
      "idx is 2338\n",
      "4240\n",
      "idx is 2339\n",
      "3538\n",
      "idx is 2340\n",
      "3015\n",
      "idx is 2341\n",
      "2923\n",
      "idx is 2342\n",
      "3639\n",
      "idx is 2343\n",
      "4297\n",
      "idx is 2344\n",
      "3033\n",
      "idx is 2345\n",
      "3242\n",
      "idx is 2346\n",
      "2530\n",
      "idx is 2347\n",
      "2575\n",
      "idx is 2348\n",
      "3777\n",
      "idx is 2349\n",
      "2845\n",
      "idx is 2350\n",
      "2810\n",
      "idx is 2351\n",
      "4292\n",
      "idx is 2352\n",
      "2612\n",
      "idx is 2353\n",
      "4560\n",
      "idx is 2354\n",
      "3628\n",
      "idx is 2355\n",
      "4501\n",
      "idx is 2356\n",
      "3194\n",
      "idx is 2357\n",
      "3322\n",
      "idx is 2358\n",
      "3795\n",
      "idx is 2359\n",
      "2629\n",
      "idx is 2360\n",
      "3686\n",
      "idx is 2361\n",
      "3161\n",
      "idx is 2362\n",
      "3474\n",
      "idx is 2363\n",
      "4576\n",
      "idx is 2364\n",
      "2950\n",
      "idx is 2365\n",
      "3765\n",
      "idx is 2366\n",
      "4449\n",
      "idx is 2367\n",
      "4855\n",
      "idx is 2368\n",
      "3886\n",
      "idx is 2369\n",
      "2820\n",
      "idx is 2370\n",
      "3534\n",
      "idx is 2371\n",
      "3971\n",
      "idx is 2372\n",
      "3937\n",
      "idx is 2373\n",
      "2741\n",
      "idx is 2374\n",
      "2540\n",
      "idx is 2375\n",
      "2970\n",
      "idx is 2376\n",
      "3034\n",
      "idx is 2377\n",
      "3434\n",
      "idx is 2378\n",
      "3727\n",
      "idx is 2379\n",
      "4999\n",
      "idx is 2380\n",
      "2513\n",
      "idx is 2381\n",
      "3497\n",
      "idx is 2382\n",
      "4245\n",
      "idx is 2383\n",
      "3800\n",
      "idx is 2384\n",
      "2710\n",
      "idx is 2385\n",
      "2816\n",
      "idx is 2386\n",
      "2641\n",
      "idx is 2387\n",
      "4895\n",
      "idx is 2388\n",
      "4957\n",
      "idx is 2389\n",
      "3271\n",
      "idx is 2390\n",
      "2782\n",
      "idx is 2391\n",
      "3961\n",
      "idx is 2392\n",
      "3298\n",
      "idx is 2393\n",
      "4298\n",
      "idx is 2394\n",
      "2855\n",
      "idx is 2395\n",
      "2521\n",
      "idx is 2396\n",
      "4310\n",
      "idx is 2397\n",
      "3190\n",
      "idx is 2398\n",
      "4076\n",
      "idx is 2399\n",
      "4600\n",
      "idx is 2400\n",
      "3399\n",
      "idx is 2401\n",
      "4819\n",
      "idx is 2402\n",
      "3557\n",
      "idx is 2403\n",
      "2988\n",
      "idx is 2404\n",
      "3947\n",
      "idx is 2405\n",
      "4759\n",
      "idx is 2406\n",
      "2559\n",
      "idx is 2407\n",
      "3679\n",
      "idx is 2408\n",
      "3325\n",
      "idx is 2409\n",
      "4668\n",
      "idx is 2410\n",
      "3645\n",
      "idx is 2411\n",
      "3369\n",
      "idx is 2412\n",
      "2765\n",
      "idx is 2413\n",
      "4391\n",
      "idx is 2414\n",
      "4116\n",
      "idx is 2415\n",
      "4067\n",
      "idx is 2416\n",
      "4178\n",
      "idx is 2417\n",
      "3461\n",
      "idx is 2418\n",
      "3515\n",
      "idx is 2419\n",
      "2940\n",
      "idx is 2420\n",
      "4190\n",
      "idx is 2421\n",
      "2871\n",
      "idx is 2422\n",
      "3118\n",
      "idx is 2423\n",
      "4223\n",
      "idx is 2424\n",
      "2952\n",
      "idx is 2425\n",
      "3192\n",
      "idx is 2426\n",
      "2810\n",
      "idx is 2427\n",
      "4364\n",
      "idx is 2428\n",
      "2961\n",
      "idx is 2429\n",
      "3464\n",
      "idx is 2430\n",
      "4903\n",
      "idx is 2431\n",
      "3219\n",
      "idx is 2432\n",
      "4575\n",
      "idx is 2433\n",
      "3350\n",
      "idx is 2434\n",
      "4468\n",
      "idx is 2435\n",
      "4164\n",
      "idx is 2436\n",
      "4915\n",
      "idx is 2437\n",
      "4048\n",
      "idx is 2438\n",
      "2793\n",
      "idx is 2439\n",
      "3766\n",
      "idx is 2440\n",
      "4128\n",
      "idx is 2441\n",
      "2685\n",
      "idx is 2442\n",
      "3789\n",
      "idx is 2443\n",
      "4393\n",
      "idx is 2444\n",
      "3092\n",
      "idx is 2445\n",
      "3848\n",
      "idx is 2446\n",
      "3725\n",
      "idx is 2447\n",
      "3326\n",
      "idx is 2448\n",
      "2889\n",
      "idx is 2449\n",
      "3391\n",
      "idx is 2450\n",
      "4418\n",
      "idx is 2451\n",
      "3482\n",
      "idx is 2452\n",
      "4271\n",
      "idx is 2453\n",
      "3438\n",
      "idx is 2454\n",
      "4438\n",
      "idx is 2455\n",
      "2787\n",
      "idx is 2456\n",
      "3554\n",
      "idx is 2457\n",
      "4124\n",
      "idx is 2458\n",
      "4614\n",
      "idx is 2459\n",
      "3538\n",
      "idx is 2460\n",
      "4073\n",
      "idx is 2461\n",
      "2551\n",
      "idx is 2462\n",
      "3531\n",
      "idx is 2463\n",
      "3117\n",
      "idx is 2464\n",
      "3582\n",
      "idx is 2465\n",
      "3949\n",
      "idx is 2466\n",
      "4481\n",
      "idx is 2467\n",
      "4805\n",
      "idx is 2468\n",
      "2608\n",
      "idx is 2469\n",
      "3558\n",
      "idx is 2470\n",
      "4755\n",
      "idx is 2471\n",
      "3564\n",
      "idx is 2472\n",
      "2678\n",
      "idx is 2473\n",
      "2687\n",
      "idx is 2474\n",
      "3909\n",
      "idx is 2475\n",
      "4733\n",
      "idx is 2476\n",
      "2874\n",
      "idx is 2477\n",
      "2538\n",
      "idx is 2478\n",
      "4235\n",
      "idx is 2479\n",
      "4608\n",
      "idx is 2480\n",
      "3153\n",
      "idx is 2481\n",
      "3921\n",
      "idx is 2482\n",
      "2516\n",
      "idx is 2483\n",
      "3452\n",
      "idx is 2484\n",
      "4160\n",
      "idx is 2485\n",
      "2604\n",
      "idx is 2486\n",
      "4402\n",
      "idx is 2487\n",
      "3127\n",
      "idx is 2488\n",
      "2507\n",
      "idx is 2489\n",
      "4945\n",
      "idx is 2490\n",
      "2937\n",
      "idx is 2491\n",
      "2944\n",
      "idx is 2492\n",
      "2932\n",
      "idx is 2493\n",
      "2590\n",
      "idx is 2494\n",
      "4986\n",
      "idx is 2495\n",
      "4094\n",
      "idx is 2496\n",
      "2854\n",
      "idx is 2497\n",
      "4464\n",
      "idx is 2498\n",
      "3887\n",
      "idx is 2499\n",
      "2845\n",
      "idx is 2500\n",
      "4695\n",
      "idx is 2501\n",
      "3970\n",
      "idx is 2502\n",
      "4756\n",
      "idx is 2503\n",
      "2564\n",
      "idx is 2504\n",
      "2573\n",
      "idx is 2505\n",
      "4846\n",
      "idx is 2506\n",
      "3688\n",
      "idx is 2507\n",
      "4495\n",
      "idx is 2508\n",
      "3604\n",
      "idx is 2509\n",
      "4382\n",
      "idx is 2510\n",
      "4483\n",
      "idx is 2511\n",
      "3110\n",
      "idx is 2512\n",
      "2864\n",
      "idx is 2513\n",
      "3630\n",
      "idx is 2514\n",
      "4295\n",
      "idx is 2515\n",
      "2984\n",
      "idx is 2516\n",
      "4955\n",
      "idx is 2517\n",
      "2752\n",
      "idx is 2518\n",
      "4561\n",
      "idx is 2519\n",
      "4834\n",
      "idx is 2520\n",
      "3904\n",
      "idx is 2521\n",
      "3104\n",
      "idx is 2522\n",
      "3804\n",
      "idx is 2523\n",
      "2707\n",
      "idx is 2524\n",
      "4637\n",
      "idx is 2525\n",
      "3898\n",
      "idx is 2526\n",
      "4927\n",
      "idx is 2527\n",
      "2678\n",
      "idx is 2528\n",
      "4275\n",
      "idx is 2529\n",
      "3139\n",
      "idx is 2530\n",
      "4292\n",
      "idx is 2531\n",
      "4356\n",
      "idx is 2532\n",
      "3664\n",
      "idx is 2533\n",
      "2822\n",
      "idx is 2534\n",
      "4472\n",
      "idx is 2535\n",
      "3809\n",
      "idx is 2536\n",
      "2968\n",
      "idx is 2537\n",
      "4218\n",
      "idx is 2538\n",
      "2526\n",
      "idx is 2539\n",
      "4321\n",
      "idx is 2540\n",
      "3169\n",
      "idx is 2541\n",
      "3823\n",
      "idx is 2542\n",
      "3366\n",
      "idx is 2543\n",
      "2913\n",
      "idx is 2544\n",
      "3790\n",
      "idx is 2545\n",
      "2800\n",
      "idx is 2546\n",
      "4554\n",
      "idx is 2547\n",
      "3970\n",
      "idx is 2548\n",
      "4717\n",
      "idx is 2549\n",
      "4879\n",
      "idx is 2550\n",
      "3038\n",
      "idx is 2551\n",
      "4911\n",
      "idx is 2552\n",
      "4038\n",
      "idx is 2553\n",
      "4720\n",
      "idx is 2554\n",
      "2839\n",
      "idx is 2555\n",
      "3986\n",
      "idx is 2556\n",
      "4844\n",
      "idx is 2557\n",
      "3478\n",
      "idx is 2558\n",
      "2896\n",
      "idx is 2559\n",
      "4996\n",
      "idx is 2560\n",
      "4907\n",
      "idx is 2561\n",
      "3950\n",
      "idx is 2562\n",
      "3861\n",
      "idx is 2563\n",
      "2875\n",
      "idx is 2564\n",
      "4394\n",
      "idx is 2565\n",
      "4870\n",
      "idx is 2566\n",
      "4138\n",
      "idx is 2567\n",
      "3535\n",
      "idx is 2568\n",
      "3681\n",
      "idx is 2569\n",
      "4409\n",
      "idx is 2570\n",
      "3983\n",
      "idx is 2571\n",
      "4031\n",
      "idx is 2572\n",
      "2839\n",
      "idx is 2573\n",
      "3776\n",
      "idx is 2574\n",
      "3484\n",
      "idx is 2575\n",
      "3079\n",
      "idx is 2576\n",
      "4954\n",
      "idx is 2577\n",
      "4255\n",
      "idx is 2578\n",
      "3376\n",
      "idx is 2579\n",
      "3453\n",
      "idx is 2580\n",
      "4657\n",
      "idx is 2581\n",
      "3241\n",
      "idx is 2582\n",
      "2869\n",
      "idx is 2583\n",
      "3509\n",
      "idx is 2584\n",
      "3011\n",
      "idx is 2585\n",
      "4938\n",
      "idx is 2586\n",
      "4226\n",
      "idx is 2587\n",
      "3849\n",
      "idx is 2588\n",
      "2619\n",
      "idx is 2589\n",
      "3201\n",
      "idx is 2590\n",
      "3191\n",
      "idx is 2591\n",
      "4692\n",
      "idx is 2592\n",
      "3529\n",
      "idx is 2593\n",
      "3400\n",
      "idx is 2594\n",
      "2761\n",
      "idx is 2595\n",
      "4196\n",
      "idx is 2596\n",
      "3934\n",
      "idx is 2597\n",
      "4954\n",
      "idx is 2598\n",
      "3448\n",
      "idx is 2599\n",
      "2956\n",
      "idx is 2600\n",
      "3368\n",
      "idx is 2601\n",
      "3219\n",
      "idx is 2602\n",
      "3257\n",
      "idx is 2603\n",
      "3387\n",
      "idx is 2604\n",
      "4718\n",
      "idx is 2605\n",
      "4867\n",
      "idx is 2606\n",
      "3654\n",
      "idx is 2607\n",
      "3881\n",
      "idx is 2608\n",
      "4215\n",
      "idx is 2609\n",
      "4474\n",
      "idx is 2610\n",
      "4643\n",
      "idx is 2611\n",
      "3738\n",
      "idx is 2612\n",
      "3985\n",
      "idx is 2613\n",
      "2745\n",
      "idx is 2614\n",
      "3197\n",
      "idx is 2615\n",
      "3852\n",
      "idx is 2616\n",
      "2710\n",
      "idx is 2617\n",
      "4978\n",
      "idx is 2618\n",
      "2640\n",
      "idx is 2619\n",
      "4296\n",
      "idx is 2620\n",
      "4264\n",
      "idx is 2621\n",
      "4117\n",
      "idx is 2622\n",
      "3580\n",
      "idx is 2623\n",
      "3140\n",
      "idx is 2624\n",
      "3784\n",
      "idx is 2625\n",
      "3369\n",
      "idx is 2626\n",
      "3819\n",
      "idx is 2627\n",
      "4565\n",
      "idx is 2628\n",
      "4188\n",
      "idx is 2629\n",
      "3297\n",
      "idx is 2630\n",
      "4897\n",
      "idx is 2631\n",
      "3006\n",
      "idx is 2632\n",
      "4704\n",
      "idx is 2633\n",
      "4510\n",
      "idx is 2634\n",
      "3058\n",
      "idx is 2635\n",
      "3299\n",
      "idx is 2636\n",
      "2500\n",
      "idx is 2637\n",
      "4491\n",
      "idx is 2638\n",
      "2894\n",
      "idx is 2639\n",
      "3047\n",
      "idx is 2640\n",
      "2868\n",
      "idx is 2641\n",
      "4957\n",
      "idx is 2642\n",
      "4756\n",
      "idx is 2643\n",
      "3802\n",
      "idx is 2644\n",
      "3485\n",
      "idx is 2645\n",
      "4147\n",
      "idx is 2646\n",
      "3526\n",
      "idx is 2647\n",
      "2934\n",
      "idx is 2648\n",
      "3662\n",
      "idx is 2649\n",
      "4232\n",
      "idx is 2650\n",
      "3528\n",
      "idx is 2651\n",
      "4751\n",
      "idx is 2652\n",
      "3757\n",
      "idx is 2653\n",
      "4318\n",
      "idx is 2654\n",
      "3386\n",
      "idx is 2655\n",
      "3739\n",
      "idx is 2656\n",
      "3854\n",
      "idx is 2657\n",
      "3755\n",
      "idx is 2658\n",
      "3364\n",
      "idx is 2659\n",
      "3544\n",
      "idx is 2660\n",
      "4567\n",
      "idx is 2661\n",
      "4770\n",
      "idx is 2662\n",
      "3765\n",
      "idx is 2663\n",
      "3720\n",
      "idx is 2664\n",
      "3826\n",
      "idx is 2665\n",
      "4818\n",
      "idx is 2666\n",
      "4336\n",
      "idx is 2667\n",
      "2977\n",
      "idx is 2668\n",
      "2549\n",
      "idx is 2669\n",
      "4403\n",
      "idx is 2670\n",
      "4248\n",
      "idx is 2671\n",
      "2802\n",
      "idx is 2672\n",
      "4703\n",
      "idx is 2673\n",
      "4789\n",
      "idx is 2674\n",
      "4247\n",
      "idx is 2675\n",
      "2537\n",
      "idx is 2676\n",
      "4431\n",
      "idx is 2677\n",
      "4405\n",
      "idx is 2678\n",
      "4315\n",
      "idx is 2679\n",
      "3321\n",
      "idx is 2680\n",
      "3942\n",
      "idx is 2681\n",
      "4006\n",
      "idx is 2682\n",
      "3720\n",
      "idx is 2683\n",
      "3275\n",
      "idx is 2684\n",
      "2662\n",
      "idx is 2685\n",
      "4422\n",
      "idx is 2686\n",
      "2725\n",
      "idx is 2687\n",
      "3232\n",
      "idx is 2688\n",
      "4200\n",
      "idx is 2689\n",
      "3403\n",
      "idx is 2690\n",
      "4877\n",
      "idx is 2691\n",
      "3432\n",
      "idx is 2692\n",
      "4294\n",
      "idx is 2693\n",
      "2517\n",
      "idx is 2694\n",
      "4621\n",
      "idx is 2695\n",
      "4338\n",
      "idx is 2696\n",
      "4583\n",
      "idx is 2697\n",
      "2773\n",
      "idx is 2698\n",
      "4111\n",
      "idx is 2699\n",
      "4895\n",
      "idx is 2700\n",
      "3835\n",
      "idx is 2701\n",
      "2566\n",
      "idx is 2702\n",
      "2545\n",
      "idx is 2703\n",
      "4602\n",
      "idx is 2704\n",
      "3700\n",
      "idx is 2705\n",
      "3990\n",
      "idx is 2706\n",
      "3583\n",
      "idx is 2707\n",
      "4047\n",
      "idx is 2708\n",
      "4326\n",
      "idx is 2709\n",
      "2786\n",
      "idx is 2710\n",
      "3814\n",
      "idx is 2711\n",
      "3449\n",
      "idx is 2712\n",
      "3353\n",
      "idx is 2713\n",
      "2875\n",
      "idx is 2714\n",
      "3270\n",
      "idx is 2715\n",
      "4612\n",
      "idx is 2716\n",
      "2562\n",
      "idx is 2717\n",
      "4238\n",
      "idx is 2718\n",
      "2656\n",
      "idx is 2719\n",
      "4407\n",
      "idx is 2720\n",
      "4476\n",
      "idx is 2721\n",
      "4882\n",
      "idx is 2722\n",
      "3546\n",
      "idx is 2723\n",
      "3928\n",
      "idx is 2724\n",
      "3634\n",
      "idx is 2725\n",
      "3604\n",
      "idx is 2726\n",
      "2850\n",
      "idx is 2727\n",
      "3399\n",
      "idx is 2728\n",
      "2615\n",
      "idx is 2729\n",
      "3443\n",
      "idx is 2730\n",
      "3899\n",
      "idx is 2731\n",
      "4619\n",
      "idx is 2732\n",
      "4824\n",
      "idx is 2733\n",
      "3936\n",
      "idx is 2734\n",
      "4111\n",
      "idx is 2735\n",
      "4225\n",
      "idx is 2736\n",
      "4337\n",
      "idx is 2737\n",
      "2821\n",
      "idx is 2738\n",
      "3170\n",
      "idx is 2739\n",
      "3127\n",
      "idx is 2740\n",
      "4965\n",
      "idx is 2741\n",
      "3153\n",
      "idx is 2742\n",
      "3178\n",
      "idx is 2743\n",
      "4961\n",
      "idx is 2744\n",
      "4252\n",
      "idx is 2745\n",
      "3723\n",
      "idx is 2746\n",
      "4648\n",
      "idx is 2747\n",
      "4854\n",
      "idx is 2748\n",
      "2669\n",
      "idx is 2749\n",
      "3027\n",
      "idx is 2750\n",
      "4247\n",
      "idx is 2751\n",
      "4196\n",
      "idx is 2752\n",
      "4133\n",
      "idx is 2753\n",
      "3670\n",
      "idx is 2754\n",
      "3320\n",
      "idx is 2755\n",
      "4066\n",
      "idx is 2756\n",
      "2559\n",
      "idx is 2757\n",
      "3939\n",
      "idx is 2758\n",
      "4469\n",
      "idx is 2759\n",
      "3288\n",
      "idx is 2760\n",
      "3573\n",
      "idx is 2761\n",
      "3071\n",
      "idx is 2762\n",
      "4082\n",
      "idx is 2763\n",
      "4858\n",
      "idx is 2764\n",
      "4539\n",
      "idx is 2765\n",
      "2770\n",
      "idx is 2766\n",
      "4177\n",
      "idx is 2767\n",
      "2627\n",
      "idx is 2768\n",
      "4844\n",
      "idx is 2769\n",
      "3466\n",
      "idx is 2770\n",
      "2500\n",
      "idx is 2771\n",
      "3994\n",
      "idx is 2772\n",
      "4528\n",
      "idx is 2773\n",
      "3408\n",
      "idx is 2774\n",
      "3495\n",
      "idx is 2775\n",
      "3332\n",
      "idx is 2776\n",
      "3155\n",
      "idx is 2777\n",
      "3917\n",
      "idx is 2778\n",
      "3364\n",
      "idx is 2779\n",
      "4064\n",
      "idx is 2780\n",
      "4951\n",
      "idx is 2781\n",
      "4009\n",
      "idx is 2782\n",
      "4622\n",
      "idx is 2783\n",
      "3719\n",
      "idx is 2784\n",
      "4104\n",
      "idx is 2785\n",
      "2951\n",
      "idx is 2786\n",
      "3546\n",
      "idx is 2787\n",
      "4908\n",
      "idx is 2788\n",
      "3890\n",
      "idx is 2789\n",
      "3418\n",
      "idx is 2790\n",
      "4930\n",
      "idx is 2791\n",
      "3243\n",
      "idx is 2792\n",
      "3797\n",
      "idx is 2793\n",
      "4383\n",
      "idx is 2794\n",
      "4600\n",
      "idx is 2795\n",
      "2677\n",
      "idx is 2796\n",
      "3390\n",
      "idx is 2797\n",
      "2526\n",
      "idx is 2798\n",
      "2832\n",
      "idx is 2799\n",
      "3090\n",
      "idx is 2800\n",
      "4334\n",
      "idx is 2801\n",
      "2595\n",
      "idx is 2802\n",
      "2804\n",
      "idx is 2803\n",
      "3273\n",
      "idx is 2804\n",
      "3074\n",
      "idx is 2805\n",
      "3562\n",
      "idx is 2806\n",
      "3592\n",
      "idx is 2807\n",
      "2976\n",
      "idx is 2808\n",
      "3497\n",
      "idx is 2809\n",
      "4963\n",
      "idx is 2810\n",
      "4871\n",
      "idx is 2811\n",
      "3310\n",
      "idx is 2812\n",
      "3734\n",
      "idx is 2813\n",
      "2569\n",
      "idx is 2814\n",
      "3952\n",
      "idx is 2815\n",
      "4740\n",
      "idx is 2816\n",
      "3232\n",
      "idx is 2817\n",
      "3420\n",
      "idx is 2818\n",
      "2779\n",
      "idx is 2819\n",
      "3257\n",
      "idx is 2820\n",
      "4103\n",
      "idx is 2821\n",
      "3919\n",
      "idx is 2822\n",
      "2966\n",
      "idx is 2823\n",
      "2649\n",
      "idx is 2824\n",
      "3830\n",
      "idx is 2825\n",
      "4150\n",
      "idx is 2826\n",
      "4310\n",
      "idx is 2827\n",
      "4017\n",
      "idx is 2828\n",
      "3299\n",
      "idx is 2829\n",
      "4323\n",
      "idx is 2830\n",
      "2871\n",
      "idx is 2831\n",
      "3677\n",
      "idx is 2832\n",
      "3476\n",
      "idx is 2833\n",
      "3907\n",
      "idx is 2834\n",
      "3664\n",
      "idx is 2835\n",
      "3033\n",
      "idx is 2836\n",
      "3796\n",
      "idx is 2837\n",
      "4187\n",
      "idx is 2838\n",
      "4718\n",
      "idx is 2839\n",
      "4927\n",
      "idx is 2840\n",
      "3921\n",
      "idx is 2841\n",
      "4285\n",
      "idx is 2842\n",
      "3910\n",
      "idx is 2843\n",
      "3844\n",
      "idx is 2844\n",
      "2722\n",
      "idx is 2845\n",
      "2678\n",
      "idx is 2846\n",
      "3286\n",
      "idx is 2847\n",
      "3933\n",
      "idx is 2848\n",
      "4443\n",
      "idx is 2849\n",
      "2894\n",
      "idx is 2850\n",
      "4490\n",
      "idx is 2851\n",
      "3756\n",
      "idx is 2852\n",
      "2503\n",
      "idx is 2853\n",
      "3726\n",
      "idx is 2854\n",
      "3149\n",
      "idx is 2855\n",
      "3811\n",
      "idx is 2856\n",
      "3047\n",
      "idx is 2857\n",
      "4756\n",
      "idx is 2858\n",
      "4386\n",
      "idx is 2859\n",
      "4089\n",
      "idx is 2860\n",
      "3429\n",
      "idx is 2861\n",
      "4155\n",
      "idx is 2862\n",
      "2536\n",
      "idx is 2863\n",
      "3853\n",
      "idx is 2864\n",
      "3255\n",
      "idx is 2865\n",
      "4278\n",
      "idx is 2866\n",
      "3003\n",
      "idx is 2867\n",
      "3217\n",
      "idx is 2868\n",
      "3416\n",
      "idx is 2869\n",
      "2762\n",
      "idx is 2870\n",
      "3838\n",
      "idx is 2871\n",
      "2918\n",
      "idx is 2872\n",
      "3497\n",
      "idx is 2873\n",
      "4556\n",
      "idx is 2874\n",
      "4953\n",
      "idx is 2875\n",
      "2900\n",
      "idx is 2876\n",
      "3096\n",
      "idx is 2877\n",
      "4842\n",
      "idx is 2878\n",
      "4116\n",
      "idx is 2879\n",
      "4633\n",
      "idx is 2880\n",
      "4352\n",
      "idx is 2881\n",
      "3355\n",
      "idx is 2882\n",
      "4048\n",
      "idx is 2883\n",
      "4775\n",
      "idx is 2884\n",
      "4103\n",
      "idx is 2885\n",
      "3488\n",
      "idx is 2886\n",
      "2545\n",
      "idx is 2887\n",
      "4768\n",
      "idx is 2888\n",
      "2994\n",
      "idx is 2889\n",
      "3454\n",
      "idx is 2890\n",
      "2711\n",
      "idx is 2891\n",
      "3914\n",
      "idx is 2892\n",
      "2603\n",
      "idx is 2893\n",
      "3851\n",
      "idx is 2894\n",
      "2720\n",
      "idx is 2895\n",
      "3015\n",
      "idx is 2896\n",
      "4087\n",
      "idx is 2897\n",
      "3224\n",
      "idx is 2898\n",
      "4564\n",
      "idx is 2899\n",
      "3895\n",
      "idx is 2900\n",
      "4798\n",
      "idx is 2901\n",
      "4365\n",
      "idx is 2902\n",
      "2723\n",
      "idx is 2903\n",
      "3883\n",
      "idx is 2904\n",
      "3790\n",
      "idx is 2905\n",
      "3670\n",
      "idx is 2906\n",
      "4918\n",
      "idx is 2907\n",
      "4075\n",
      "idx is 2908\n",
      "3380\n",
      "idx is 2909\n",
      "3903\n",
      "idx is 2910\n",
      "4748\n",
      "idx is 2911\n",
      "4875\n",
      "idx is 2912\n",
      "2809\n",
      "idx is 2913\n",
      "3726\n",
      "idx is 2914\n",
      "4306\n",
      "idx is 2915\n",
      "4508\n",
      "idx is 2916\n",
      "2813\n",
      "idx is 2917\n",
      "3604\n",
      "idx is 2918\n",
      "3780\n",
      "idx is 2919\n",
      "2925\n",
      "idx is 2920\n",
      "3687\n",
      "idx is 2921\n",
      "2786\n",
      "idx is 2922\n",
      "3528\n",
      "idx is 2923\n",
      "3595\n",
      "idx is 2924\n",
      "3272\n",
      "idx is 2925\n",
      "3408\n",
      "idx is 2926\n",
      "2828\n",
      "idx is 2927\n",
      "2649\n",
      "idx is 2928\n",
      "3513\n",
      "idx is 2929\n",
      "3920\n",
      "idx is 2930\n",
      "2593\n",
      "idx is 2931\n",
      "4144\n",
      "idx is 2932\n",
      "4977\n",
      "idx is 2933\n",
      "4072\n",
      "idx is 2934\n",
      "2570\n",
      "idx is 2935\n",
      "2850\n",
      "idx is 2936\n",
      "3431\n",
      "idx is 2937\n",
      "2917\n",
      "idx is 2938\n",
      "3231\n",
      "idx is 2939\n",
      "4999\n",
      "idx is 2940\n",
      "2637\n",
      "idx is 2941\n",
      "3269\n",
      "idx is 2942\n",
      "2679\n",
      "idx is 2943\n",
      "3307\n",
      "idx is 2944\n",
      "4400\n",
      "idx is 2945\n",
      "2935\n",
      "idx is 2946\n",
      "4579\n",
      "idx is 2947\n",
      "3791\n",
      "idx is 2948\n",
      "4784\n",
      "idx is 2949\n",
      "3578\n",
      "idx is 2950\n",
      "4948\n",
      "idx is 2951\n",
      "4645\n",
      "idx is 2952\n",
      "4482\n",
      "idx is 2953\n",
      "4988\n",
      "idx is 2954\n",
      "3458\n",
      "idx is 2955\n",
      "4700\n",
      "idx is 2956\n",
      "3938\n",
      "idx is 2957\n",
      "3210\n",
      "idx is 2958\n",
      "3485\n",
      "idx is 2959\n",
      "4655\n",
      "idx is 2960\n",
      "4156\n",
      "idx is 2961\n",
      "3847\n",
      "idx is 2962\n",
      "4911\n",
      "idx is 2963\n",
      "4823\n",
      "idx is 2964\n",
      "4023\n",
      "idx is 2965\n",
      "3246\n",
      "idx is 2966\n",
      "4691\n",
      "idx is 2967\n",
      "3543\n",
      "idx is 2968\n",
      "2823\n",
      "idx is 2969\n",
      "4587\n",
      "idx is 2970\n",
      "3152\n",
      "idx is 2971\n",
      "3799\n",
      "idx is 2972\n",
      "4093\n",
      "idx is 2973\n",
      "3037\n",
      "idx is 2974\n",
      "4284\n",
      "idx is 2975\n",
      "3666\n",
      "idx is 2976\n",
      "2529\n",
      "idx is 2977\n",
      "3227\n",
      "idx is 2978\n",
      "2788\n",
      "idx is 2979\n",
      "3113\n",
      "idx is 2980\n",
      "3843\n",
      "idx is 2981\n",
      "4864\n",
      "idx is 2982\n",
      "3036\n",
      "idx is 2983\n",
      "2855\n",
      "idx is 2984\n",
      "3290\n",
      "idx is 2985\n",
      "3264\n",
      "idx is 2986\n",
      "3848\n",
      "idx is 2987\n",
      "4910\n",
      "idx is 2988\n",
      "2861\n",
      "idx is 2989\n",
      "4965\n",
      "idx is 2990\n",
      "4220\n",
      "idx is 2991\n",
      "4382\n",
      "idx is 2992\n",
      "3026\n",
      "idx is 2993\n",
      "2637\n",
      "idx is 2994\n",
      "3214\n",
      "idx is 2995\n",
      "3415\n",
      "idx is 2996\n",
      "2589\n",
      "idx is 2997\n",
      "4728\n",
      "idx is 2998\n",
      "4656\n",
      "idx is 2999\n",
      "3058\n",
      "idx is 3000\n",
      "2852\n",
      "idx is 3001\n",
      "2593\n",
      "idx is 3002\n",
      "3313\n",
      "idx is 3003\n",
      "2627\n",
      "idx is 3004\n",
      "2788\n",
      "idx is 3005\n",
      "4124\n",
      "idx is 3006\n",
      "4044\n",
      "idx is 3007\n",
      "4706\n",
      "idx is 3008\n",
      "2728\n",
      "idx is 3009\n",
      "4587\n",
      "idx is 3010\n",
      "4958\n",
      "idx is 3011\n",
      "3806\n",
      "idx is 3012\n",
      "3074\n",
      "idx is 3013\n",
      "4653\n",
      "idx is 3014\n",
      "3625\n",
      "idx is 3015\n",
      "4329\n",
      "idx is 3016\n",
      "4162\n",
      "idx is 3017\n",
      "4156\n",
      "idx is 3018\n",
      "4217\n",
      "idx is 3019\n",
      "4553\n",
      "idx is 3020\n",
      "4803\n",
      "idx is 3021\n",
      "4505\n",
      "idx is 3022\n",
      "4642\n",
      "idx is 3023\n",
      "4416\n",
      "idx is 3024\n",
      "3707\n",
      "idx is 3025\n",
      "2782\n",
      "idx is 3026\n",
      "4485\n",
      "idx is 3027\n",
      "4886\n",
      "idx is 3028\n",
      "4535\n"
     ]
    },
    {
     "ename": "",
     "evalue": "",
     "output_type": "error",
     "traceback": [
      "\u001b[1;31mThe Kernel crashed while executing code in the the current cell or a previous cell. Please review the code in the cell(s) to identify a possible cause of the failure. Click <a href='https://aka.ms/vscodeJupyterKernelCrash'>here</a> for more info. View Jupyter <a href='command:jupyter.viewOutput'>log</a> for further details."
     ]
    }
   ],
   "source": [
    "group = 4\n",
    "matrix_all = np.load(f\"/shareb/mliu/evaluate_impute/data/simulation_hic/Cell2020/hic/matrix/Cell2020_sample_matrix{group}_250kb.npy\")\n",
    "for idx, matrix in enumerate(matrix_all):\n",
    "    print(f\"idx is {idx}\")\n",
    "    copy = str(idx)\n",
    "    #要保存的\n",
    "    filename = f'/shareb/mliu/evaluate_impute/data/simulation_hic/Cell2020/hic/dip_pairs/Cell2020_group{group}_copy{copy}.pairs.gz'\n",
    "\n",
    "    idx1, idx2 = np.where(np.triu(matrix) > 0)\n",
    "    values = matrix[idx1, idx2]\n",
    "\n",
    "    idx1_repeats = [i for i, v in zip(idx1, values) for _ in range(int(v))]\n",
    "    idx2_repeats = [i for i, v in zip(idx2, values) for _ in range(int(v))]\n",
    "\n",
    "    pos1 = [mid_positions_cell_250kb[i] for i in idx1_repeats]  \n",
    "    pos2 = [mid_positions_cell_250kb[j] for j in idx2_repeats]\n",
    "\n",
    "\n",
    "    pairs = pd.DataFrame({\n",
    "    'chr1': ['chr2'] * len(pos1),\n",
    "    'pos1': pos1,\n",
    "    'chr2': ['chr2'] * len(pos2),\n",
    "    'pos2': pos2,\n",
    "    })\n",
    "\n",
    "    pairs['readID'] = \".\"\n",
    "    pairs['strand1'] = \"+\"\n",
    "    pairs['strand2'] = \"+\"\n",
    "    pairs['phase0'] = 0\n",
    "    pairs['phase1'] = 0\n",
    "    pairs['p00'] = 1.000\n",
    "    pairs['p01'] = 0.000\n",
    "    pairs['p10'] = 0.000\n",
    "    pairs['p11'] = 0.000\n",
    "    pairs = pairs[['readID','chr1','pos1',\"chr2\",\"pos2\",\"strand1\",\"strand2\",\"phase0\",\"phase1\",\"p00\",\"p01\",\"p10\",\"p11\"]]\n",
    "    # pairs = pairs[['chr1','pos1',\"chr2\",\"pos2\"]]\n",
    "    # 指定文件名\n",
    "    print(len(pairs))\n",
    "    csv_buffer = StringIO()\n",
    "    for line in header_lines_cell:\n",
    "        csv_buffer.write(line + '\\n')\n",
    "    pairs.to_csv(csv_buffer, index=False, header=False,sep=\"\\t\")\n",
    "\n",
    "    with gzip.open(filename, 'wt') as f:\n",
    "        f.write(csv_buffer.getvalue())\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "impute",
   "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.8.18"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}