├── f1max-master ├── f1max ├── test.sh ├── README.md └── f1max.c ├── data_example └── graph.txt ├── logs.txt ├── graph_building.py ├── Propagation_Phase.py ├── README.md ├── LICENSE ├── seed_set_expansion.py ├── Graclus_centers.py ├── init.py ├── compare.py ├── com1.txt └── com2.txt /f1max-master/f1max: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/greyfoxamine/Overlapping-Community-Detection-Using-Seed-Set-Expansion-/HEAD/f1max-master/f1max -------------------------------------------------------------------------------- /data_example/graph.txt: -------------------------------------------------------------------------------- 1 | 1 2 2 | 2 3 3 | 1 6 4 | 1 5 5 | 3 9 6 | 2 1 7 | 3 2 8 | 9 3 9 | 3 5 10 | 6 7 11 | 3 7 12 | 2 7 13 | 5 8 14 | 2 9 15 | 3 5 16 | 9 1 17 | 10 3 18 | 1 5 -------------------------------------------------------------------------------- /f1max-master/test.sh: -------------------------------------------------------------------------------- 1 | s1="$(./f1max $1 $2)" 2 | s2="$(./f1max $2 $1)" 3 | s3=$(bc <<<"scale=6; ($s1+$s2)/2") 4 | s4=$(bc <<<"scale=6; 2*$s1*$s2/($s1+$s2)") 5 | echo $s1, $s2, $s3, $s4 6 | -------------------------------------------------------------------------------- /logs.txt: -------------------------------------------------------------------------------- 1 | Filtering phase in :0.7024049758911133 2 | seeding phase in :44.99510598182678 3 | seed set expansion phase in :0.03710007667541504 4 | building graph with community colors in :0.03859090805053711 5 | -------------------------------------------------------------------------------- /graph_building.py: -------------------------------------------------------------------------------- 1 | 2 | import matplotlib.pyplot as plt 3 | import networkx as nx 4 | import sys 5 | 6 | def file_graph_show( g,i ): 7 | plt.figure(i) 8 | 9 | sp=nx.spring_layout(g) 10 | 11 | plt.axis('off') 12 | 13 | nx.draw_networkx(g,pos=sp,with_labels=False,node_size=35) 14 | return plt 15 | 16 | def file_graph_building( path ): 17 | 18 | g=nx.read_edgelist(path,create_using=nx.Graph(),nodetype=int) 19 | 20 | print nx.info(g) 21 | return g 22 | 23 | 24 | 25 | 26 | 27 | 28 | 29 | 30 | -------------------------------------------------------------------------------- /Propagation_Phase.py: -------------------------------------------------------------------------------- 1 | import matplotlib.pyplot as plt 2 | import networkx as nx 3 | import sys 4 | 5 | def difference(S, R): 6 | #R = G.copy() 7 | #R.remove_nodes_from(n for n in G if n in H) 8 | #return R 9 | print " prolongation phase processing ..." 10 | DIF=nx.Graph() 11 | for edge in S.edges_iter(): 12 | if not R.has_edge(edge[0],edge[1]): 13 | DIF.add_edge(edge[0],edge[1]) 14 | 15 | for edge in DIF.edges_iter(): 16 | R.add_edge(edge[0],edge[1]) 17 | print "after prolongation phase done" 18 | return R -------------------------------------------------------------------------------- /README.md: -------------------------------------------------------------------------------- 1 | # Overlapping-Community-Detection-Using-Seed-Set-Expansion- 2 | # How to execute 3 | Dans votre terminal ,installer les librairies suivantes avec pip en tapant les commandes suivantes: 4 | pip install networkx 5 | pip install community 6 | 7 | python init.py data_example/graph.txt :pour un petit graph 8 | 9 | python init.py data_example/facebook_combined.txt :pour un grand graph 10 | # Test et comparaison 11 | 12 | python compare.py data_example/graph.txt :pour un petit graph 13 | 14 | python compare.py data_example/facebook_combined.txt :pour un grand graph 15 | 16 | cd f1max-master/ 17 | 18 | ./f1max ../com1.txt ../com2.txt 19 | 20 | 21 | -------------------------------------------------------------------------------- /LICENSE: -------------------------------------------------------------------------------- 1 | MIT License 2 | 3 | Copyright (c) 2016 Amine Dourid 4 | 5 | Permission is hereby granted, free of charge, to any person obtaining a copy 6 | of this software and associated documentation files (the "Software"), to deal 7 | in the Software without restriction, including without limitation the rights 8 | to use, copy, modify, merge, publish, distribute, sublicense, and/or sell 9 | copies of the Software, and to permit persons to whom the Software is 10 | furnished to do so, subject to the following conditions: 11 | 12 | The above copyright notice and this permission notice shall be included in all 13 | copies or substantial portions of the Software. 14 | 15 | THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR 16 | IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, 17 | FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE 18 | AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER 19 | LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, 20 | OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE 21 | SOFTWARE. 22 | -------------------------------------------------------------------------------- /seed_set_expansion.py: -------------------------------------------------------------------------------- 1 | from __future__ import division 2 | 3 | import random as rdm 4 | 5 | import networkx as nx 6 | def seed_set_expansion(G,seeds): 7 | #star random walk 8 | expansion ={} 9 | compteur= 0 #Start execution counter 10 | execution = 0 11 | for seed in seeds: 12 | 13 | vertexid = seed 14 | VisitedVertex = {} 15 | 16 | #Execute the random walk with size 10000 (10000 steps) 17 | while compteur < G.number_of_nodes()/len(seeds): 18 | #Accumulate the amount of times each vertex is visited 19 | if vertexid in VisitedVertex: 20 | VisitedVertex[vertexid] += 1 21 | else: 22 | VisitedVertex[vertexid] = 1 23 | #Visualize the vertex neighborhood 24 | Vertex_Neighbors = G.neighbors(vertexid) 25 | #Choose a vertex from the vertex neighborhood to start the next random walk 26 | vertexid = rdm.choice(Vertex_Neighbors) 27 | compteur = compteur + 1 28 | mostvisited = sorted(VisitedVertex, key = VisitedVertex.get,reverse = True) 29 | expansion.update({seed: mostvisited}) 30 | compteur = 0 31 | 32 | return expansion 33 | def color_building_list(G,expansion): 34 | partition ={} 35 | 36 | 37 | valeur=0 38 | for liste_value in expansion.values(): 39 | for node in liste_value: 40 | partition.update({node: valeur}) 41 | valeur = valeur+1 42 | 43 | values = [partition.get(node) for node in G.nodes()] 44 | print values 45 | for i, val in enumerate(values): 46 | if val == None: 47 | values[i]=values[i-1] 48 | return values 49 | -------------------------------------------------------------------------------- /Graclus_centers.py: -------------------------------------------------------------------------------- 1 | from __future__ import division 2 | 3 | import networkx as nx 4 | import community 5 | 6 | sigma =1 7 | 8 | def minimum_of_float_list(liste): 9 | minimum=0 10 | for i in liste: 11 | print i 12 | if float(i) 2 | #include 3 | #include 4 | 5 | #define NNODES 10000000 //max number of nodes: will increase if needed 6 | #define NCOMS 10000000 //max number of communities: will increase if needed 7 | #define SCOM 10000000 //max size of community: will increase if needed 8 | #define NCPN 10 //max number of community per nodes: will increase if needed 9 | 10 | typedef struct { 11 | unsigned s;//size of the community 12 | unsigned smax;//maximum size of the community 13 | unsigned *nodes;//nodes in the community 14 | } community; 15 | 16 | typedef struct { 17 | unsigned n;//number of nodes 18 | unsigned m;//max number of nodes 19 | 20 | unsigned **c;//c[i]=communities containing node i. 21 | unsigned *s;//s[i]=number of communities containing node i. 22 | unsigned *smax;//smax[i]=max number of communities 23 | 24 | unsigned nc;//number of communities 25 | unsigned mc;//max number of communities 26 | unsigned *size;//size[i]=size of the community i 27 | 28 | unsigned *tmp;//tmp[i] number of nodes shared by the current community and community i 29 | unsigned *list;//list of the community sharing at least one node with the current community 30 | unsigned nlist;//length of list 31 | } compare; 32 | 33 | 34 | compare *alloccompare(){ 35 | compare *comp=malloc(sizeof(compare)); 36 | comp->n=0; 37 | comp->m=NNODES; 38 | comp->c=calloc(comp->m,sizeof(unsigned*)); 39 | comp->s=calloc(comp->m,sizeof(unsigned)); 40 | comp->smax=calloc(comp->m,sizeof(unsigned)); 41 | 42 | comp->nc=0; 43 | comp->mc=NCOMS; 44 | comp->size=malloc(comp->mc*sizeof(unsigned)); 45 | 46 | comp->tmp=calloc(comp->mc,sizeof(unsigned)); 47 | comp->list=calloc(comp->mc,sizeof(unsigned)); 48 | comp->nlist=0; 49 | return comp; 50 | } 51 | 52 | community *alloccom(){ 53 | community *com=malloc(sizeof(community)); 54 | com->smax=SCOM; 55 | com->nodes=malloc(SCOM*sizeof(unsigned)); 56 | return com; 57 | } 58 | 59 | bool readlinecom(FILE* file,community* com){ 60 | char c; 61 | 62 | com->s=0; 63 | while(fscanf(file,"%u%c",com->nodes+com->s,&c)==2){ 64 | if ( ++(com->s) == com->smax) { 65 | com->smax+=SCOM; 66 | com->nodes=realloc(com->nodes,com->smax*sizeof(unsigned)); 67 | } 68 | if (c=='\n') { 69 | return 1; 70 | } 71 | } 72 | return 0; 73 | } 74 | 75 | void com2comp(community* com,compare* comp){ 76 | unsigned i,j,tmp; 77 | 78 | if (comp->nc==comp->mc){ 79 | comp->mc+=NCOMS; 80 | comp->size=realloc(comp->size,comp->mc*sizeof(unsigned)); 81 | comp->tmp=realloc(comp->tmp,comp->mc*sizeof(unsigned)); 82 | comp->list=realloc(comp->size,comp->mc*sizeof(unsigned)); 83 | for (i=comp->mc-NCOMS;imc;i++){ 84 | comp->tmp[i]=0; 85 | } 86 | } 87 | comp->size[comp->nc]=com->s; 88 | for (i=0;is;i++){ 89 | if (com->nodes[i]>comp->m){ 90 | tmp=com->nodes[i]-comp->m+NCOMS; 91 | comp->m+=tmp; 92 | comp->c=realloc(comp->c,comp->m*sizeof(unsigned*)); 93 | comp->s=realloc(comp->s,comp->m*sizeof(unsigned)); 94 | comp->smax=realloc(comp->smax,comp->m*sizeof(unsigned)); 95 | for (j=comp->m-tmp;jm;j++){ 96 | comp->c[j]=NULL; 97 | comp->s[j]=0; 98 | comp->smax[j]=0; 99 | } 100 | } 101 | if (comp->s[com->nodes[i]]==comp->smax[com->nodes[i]]){ 102 | if (comp->s[com->nodes[i]]==0){ 103 | comp->c[com->nodes[i]]=malloc(NCPN*sizeof(unsigned)); 104 | comp->smax[com->nodes[i]]=NCPN; 105 | } 106 | else { 107 | comp->smax[com->nodes[i]]+=NCPN; 108 | comp->c[com->nodes[i]]=realloc(comp->c[com->nodes[i]],comp->smax[com->nodes[i]]*sizeof(unsigned)); 109 | } 110 | } 111 | comp->c[com->nodes[i]][comp->s[com->nodes[i]]++]=comp->nc; 112 | } 113 | comp->nc++; 114 | } 115 | 116 | double comsim(community* com,compare *comp){ 117 | double sim=0,sim2; 118 | unsigned i,j; 119 | for (i=0;is;i++){ 120 | for (j=0;js[com->nodes[i]];j++){ 121 | if (comp->tmp[comp->c[com->nodes[i]][j]]==0){ 122 | comp->list[comp->nlist++]=comp->c[com->nodes[i]][j]; 123 | } 124 | comp->tmp[comp->c[com->nodes[i]][j]]++; 125 | } 126 | } 127 | 128 | for (i=0;inlist;i++){ 129 | sim2=((double)(comp->tmp[comp->list[i]]))/(com->s+comp->size[comp->list[i]]); 130 | if (sim2>sim){ 131 | sim=sim2; 132 | } 133 | comp->tmp[comp->list[i]]=0; 134 | } 135 | comp->nlist=0; 136 | return 2*sim; 137 | } 138 | 139 | int main(int argc,char** argv){ 140 | 141 | unsigned ncom=0; 142 | double score=0; 143 | community *com=alloccom(); 144 | compare *comp=alloccompare(); 145 | FILE* file; 146 | 147 | file=fopen(argv[2],"r"); 148 | while(readlinecom(file,com)){ 149 | com2comp(com,comp); 150 | } 151 | fclose(file); 152 | 153 | file=fopen(argv[1],"r"); 154 | while(readlinecom(file,com)){ 155 | score+=comsim(com,comp); 156 | ncom++; 157 | } 158 | fclose(file); 159 | 160 | score/=(double)ncom; 161 | 162 | printf("%lf\n",score); 163 | 164 | return 0; 165 | } 166 | -------------------------------------------------------------------------------- /com1.txt: -------------------------------------------------------------------------------- 1 | 2118 2624 2233 2325 2326 2408 1917 1923 2022 2559 2563 2078 2083 2607 2098 2142 2625 2116 2646 2139 2140 2654 1979 2171 2176 2190 2218 2237 2600 1918 2303 2308 2323 2376 2386 2404 2418 1912 1921 1927 2467 1966 2482 1980 2498 2500 2521 2020 2538 2542 2543 2551 2560 2061 2575 2136 2068 2478 2582 2073 2074 2593 2594 2084 2086 2088 2601 2090 2092 2093 2097 2611 2612 2615 2616 2617 2619 2108 2109 2110 2117 2122 2123 2638 2128 2131 2648 2138 2150 2154 2155 2164 2174 2175 2177 2189 2199 2201 2203 2213 2216 2226 2234 2240 2244 2253 2254 2257 2261 2266 2268 2270 2271 2272 2275 2283 1986 2307 2309 2324 2328 2329 2337 2339 2340 2347 2351 2356 2359 2103 2381 2104 2392 2393 2394 2395 2409 2410 2423 2428 2430 1925 1929 2445 2446 2447 1938 2030 1943 1945 1947 2460 1951 2544 2469 1958 1963 2026 1971 2484 2491 2495 1984 1985 1993 2027 2000 2005 2518 2011 2526 2640 2018 2019 2539 2032 2546 2549 2552 2553 2554 2556 2557 2046 2 | 0 9 91 239 8 26 55 56 67 170 339 21 31 80 109 113 122 142 168 203 245 259 261 277 304 320 322 334 40 75 77 94 98 101 104 115 117 119 123 129 134 135 172 181 187 197 213 223 232 271 282 284 291 308 312 325 332 347 3 4 7 10 13 14 22 24 25 27 28 36 41 48 59 63 79 85 88 110 128 139 140 141 150 151 152 156 164 166 169 176 184 186 188 190 195 200 201 202 211 212 224 226 229 231 234 235 244 247 252 260 264 265 272 275 276 280 281 285 286 298 303 318 323 328 330 331 338 340 344 3 | 578 107 659 632 1684 2719 414 3038 3100 3101 3106 577 583 3154 595 597 615 2685 640 647 2703 658 661 670 675 3252 3280 3288 2783 3309 3198 3318 3332 3344 1810 3382 2889 3422 2972 2991 3169 2881 3188 1110 1542 3081 3086 1555 1046 1048 3247 3102 1567 2831 553 2712 555 3116 557 3118 3120 1586 1588 1086 3253 1603 582 1609 1098 1615 592 596 1622 599 603 2833 3168 609 3172 2661 617 1642 620 3182 627 1140 630 1144 1146 1659 2686 3199 619 1156 2693 650 2701 654 662 1688 665 666 2715 2724 677 680 681 2730 1199 1717 2742 3256 628 2748 3263 3264 2849 3274 3275 2765 2770 3283 3291 2780 3296 635 3281 3311 1268 3320 2687 1793 3330 1579 2670 3299 3337 3343 3345 1813 3351 3355 1825 2853 3366 3369 2858 1839 1847 3210 3391 3393 1864 1867 3404 3406 1360 3411 3419 1376 2707 2938 2951 2959 400 664 2713 2819 418 422 2992 2718 438 3145 3002 3005 3008 3403 3013 3152 3017 3019 1187 2894 1494 1505 3240 1524 3062 1193 3156 4 | 107 1235 1554 1067 1636 1680 1204 1717 1258 1100 1014 1036 1070 1611 1617 1128 1643 1647 1659 1729 1644 1279 1835 1877 1888 1390 1598 376 1435 1454 946 1459 456 1004 512 1537 514 1539 517 525 526 1551 1556 1559 1563 1056 545 546 1573 1583 1584 1076 1078 1110 1086 1603 1095 1610 1612 1613 1622 1624 1469 1115 1125 1132 1652 1653 1146 1150 641 1669 1670 1159 1672 649 1169 1682 1685 1688 1306 1051 1704 1706 1715 969 1723 1725 1214 373 1730 1229 1742 1745 1746 1238 1242 1756 1245 1764 1259 1260 1261 1099 1273 1280 1282 1284 1799 1800 1804 1298 1667 1818 1819 1820 1309 560 1826 1158 1830 1833 1330 1843 1343 1857 1858 1347 1862 1352 1420 1866 1871 1873 1366 1879 1508 1883 1884 350 1377 1891 1380 1893 1897 1389 1396 374 1911 896 898 1412 389 1345 1772 906 395 396 397 916 917 1432 924 417 419 1444 1350 934 431 1462 1465 1467 957 1475 1476 966 455 971 461 462 1487 977 978 981 1192 1503 994 1019 484 489 1515 1005 1009 507 1021 510 5 | 3342 2951 3086 3154 2661 3207 2728 2730 2780 2974 3022 3082 3103 3108 3114 3115 2815 3152 3195 2687 3223 2717 2739 3259 2750 2754 3276 3293 2782 2784 3299 3327 2834 2836 2839 2853 2869 2874 3387 3395 3410 2787 3423 3426 2960 2973 2981 2986 2992 3033 3036 3073 3085 3089 3090 2820 3098 3104 3107 3109 3110 2908 3119 3120 3137 3339 3149 3157 3162 3169 3174 3176 3177 3178 3179 2669 3432 2675 2676 3189 3190 2679 2689 3203 3208 2698 2699 2702 3219 1684 2713 3226 3228 3233 3234 2723 2724 2727 2729 3246 2741 2677 3265 3274 2763 3279 3283 2773 3386 3295 2785 2786 2854 3302 3309 2801 3319 3320 3323 3324 3325 3332 3335 3337 2827 3340 2829 2835 3356 3360 2849 3363 3366 2856 3369 2696 2866 2867 3382 2871 2873 2879 2880 3394 3399 3214 2891 2892 2893 2901 2904 3417 2906 2907 3420 3421 3422 2912 2915 2917 2920 2931 2943 3201 2971 2972 2985 2994 3000 3001 3002 3013 3016 3024 3027 3029 3034 3035 3047 2813 3240 3411 3067 6 | 1684 2661 3051 2828 3320 3035 3082 3089 3178 3201 3247 2742 2780 3321 2813 3351 2863 3381 3034 3038 3101 3154 2665 3179 3185 3209 3219 3222 2719 3263 2758 3288 2778 2781 2784 3302 2795 2799 2818 3340 3344 3350 3374 2865 3397 3403 3411 3434 2958 2992 3005 2901 3079 3330 3416 3106 3107 3115 3120 3122 3138 3144 3152 3160 2668 3261 3347 2677 3198 3200 3207 3211 2700 3214 2869 3225 2720 2727 2728 3245 3249 2739 3255 3256 3258 3260 2762 3267 3272 3273 3274 2764 2765 3278 2768 3285 3287 3291 3298 3299 2793 3308 2798 3313 2803 3316 2787 3326 3332 2823 3337 2826 3342 2833 2947 3353 3355 2847 3360 2849 2853 2861 3378 2697 3386 2878 3391 3393 3395 2892 3406 3409 3410 3412 2702 3415 2904 3419 3422 2936 2708 2939 3307 2951 2796 1419 2968 2970 2973 2991 3001 3009 3013 3018 3020 3021 3024 3041 3046 3048 3049 3053 3070 3062 3067 2881 7 | 414 653 579 946 1006 1367 1069 1086 584 588 589 1621 599 1626 618 630 632 1673 1677 658 678 680 1238 107 1255 434 1827 1327 394 1361 1882 348 1767 573 906 916 422 423 439 983 473 1522 1530 1025 1538 515 517 1542 520 1549 527 1554 1557 1047 1048 538 1052 1114 542 544 1569 1574 1063 1799 1583 1078 1079 1080 1597 1179 1602 1603 583 590 592 1618 1620 601 606 1631 611 612 1125 1126 619 620 1645 1135 627 1653 637 1663 640 643 1674 1679 1683 1685 1113 1175 1688 367 670 672 681 1707 1710 1199 1714 1209 629 1142 1735 1736 1742 1243 523 1763 1768 551 1772 1773 1789 1791 1793 1284 1750 1287 1800 1804 1805 1811 1823 1826 1322 1329 1331 1844 1341 1001 566 567 1358 1874 1879 1881 655 1885 1888 353 354 1891 1900 1682 366 369 370 403 374 1911 376 381 384 388 395 664 1155 1431 1435 1565 413 1349 419 896 428 596 435 948 952 444 962 452 966 1015 1487 465 978 1491 475 483 484 999 488 489 1520 503 1529 1021 8 | 1684 3136 2700 3327 3397 3406 1042 3120 3171 107 641 3224 2719 2730 348 3291 2780 2811 1298 1824 3385 3386 3412 374 376 1423 417 428 1453 465 3060 1021 515 3077 3078 3079 1549 3086 3087 531 3093 3097 538 1051 1052 3206 3248 1575 553 555 3118 3360 2995 570 3168 3140 1093 3425 3035 1104 1617 3342 3162 1628 1631 1120 3170 3174 614 3176 1129 1132 3182 1649 3189 1655 1656 3198 3179 1158 3207 2696 1163 2754 1680 1172 2926 3223 2713 2716 2721 3237 1190 3241 3244 2733 3279 2735 2736 1713 3255 2819 1215 1729 1730 3267 3331 2763 1231 3280 2936 1746 3287 2778 1757 1246 1247 2784 3299 2793 2801 3315 3113 3320 2812 1277 3330 1067 3335 3337 3067 2832 3347 3348 2847 1315 3365 2854 2859 2864 2868 2782 1334 3391 1347 3396 1349 3403 2894 1872 3410 2787 1366 1369 1884 1082 3422 353 354 3427 1893 366 2928 378 1408 2947 1413 2950 2951 2955 2958 913 403 2964 406 3225 418 1443 1444 431 432 1459 436 3001 444 957 1475 3233 456 458 975 470 3034 475 3038 1531 1508 999 3051 492 497 2966 1015 1529 507 3070 511 9 | 3980 4030 3982 4023 4014 4020 4038 3986 3985 3989 3994 3998 4004 4018 3981 4002 4009 4017 4037 3988 3995 4031 3991 3997 4000 4019 4021 4026 4027 3999 4033 3992 4003 4005 4011 4036 3993 3996 4007 4013 4034 594 10 | 2507 2600 2150 2299 2058 2059 2030 2611 2619 2118 2184 2188 2266 2290 2324 2326 2331 2354 2404 2464 1979 1985 2040 2553 2563 2573 2063 2064 2575 2078 2591 2593 2086 2103 2108 2121 2123 2124 2131 2140 2142 2172 2074 2218 2590 2257 2275 2430 2090 2308 2311 2369 2370 2381 2392 2469 1962 1966 2482 2485 1984 1993 2005 2542 2033 2552 2560 2561 2055 2056 2060 2574 2579 2073 2586 2083 2084 2602 2093 2607 2098 2104 2229 2112 2631 2136 2139 2576 2154 2165 2495 2190 2200 2210 2216 2220 2233 2253 2258 2262 2276 2088 2300 2303 2549 2323 2339 2347 2348 2353 2363 2374 2615 2383 2386 2387 2391 2020 2407 2408 2414 2623 2428 1918 1929 2446 1938 1941 2454 2462 1953 2467 2477 2489 2492 1983 1986 2499 2504 2520 2526 2532 2536 2044 2539 2546 2037 2550 2551 2043 2556 2045 2559 11 | 2561 2069 2088 2095 2218 2233 2324 2363 2374 1912 2059 2090 2611 2617 2655 2112 2172 2176 2212 2235 2507 2290 2306 2309 2184 2356 2103 2430 2624 2446 1946 2461 2467 2460 1966 1979 2549 1993 1995 2018 2026 2552 2506 2056 2060 2574 2575 2576 2579 2072 2074 2075 2589 2078 2591 2592 2083 2084 1577 2604 1970 2609 2142 2615 2616 2620 2109 2623 2625 2115 2116 2630 2631 2121 2126 2639 2640 2131 2133 2646 2648 2139 2653 2654 2143 2145 2149 2154 2182 2183 136 2185 2190 2586 2210 2221 2222 2237 2242 2244 2246 2463 2257 2258 2275 2276 2278 2282 2283 2287 2291 2292 1918 2294 2601 2300 2307 2308 2311 2319 2331 2332 2333 2353 2344 2345 2347 2352 2366 2369 2370 2376 2377 2378 2381 2104 2386 2040 2398 2619 2404 2409 2410 2638 2420 2427 1916 1917 2436 1925 1926 1928 2093 1938 2451 2542 1943 1944 2543 1948 1951 1952 2472 1962 1964 2477 2482 2484 2485 2492 1984 1985 2499 2504 1994 2001 2003 2005 2006 2007 2127 2524 2525 2526 2024 2028 2029 2030 2031 2032 2033 2037 2550 2045 2046 2047 12 | 2083 2184 2206 1983 2554 2563 2052 2060 2072 2074 2084 2095 2619 2142 2655 2148 2150 2183 2210 2216 2226 2229 2240 2247 2564 2261 2266 2271 2273 2328 1925 2352 2369 2374 2381 2393 2395 1912 1920 1938 1941 1946 2462 2117 2464 1971 1985 1986 1993 2506 2007 2526 2549 2560 2561 2063 2064 2578 2073 2077 2590 2593 2598 2087 2600 1577 2602 2608 2609 2098 2654 2615 2108 2622 2623 2112 2625 2629 2118 2121 2123 2124 2287 2136 2138 2653 2154 2164 2165 2169 2495 2172 2182 2187 2188 2194 2201 2586 2213 2218 2220 2223 2233 2235 2237 2630 2242 2243 2246 2252 2253 2257 2267 2274 2276 2278 1917 2088 2290 2294 2037 2303 2306 2309 2310 2323 2329 2607 2333 2339 2340 2348 2610 2356 2359 2363 2376 2103 2386 2387 2407 2410 2414 2417 2418 1916 2429 2624 2434 1924 2454 1926 1929 2449 1942 1943 1945 1947 2461 1953 2468 2469 1962 1963 1966 2482 2484 2485 2489 2463 1984 2379 2501 2507 2508 1997 2512 2040 2002 2518 2521 2021 2030 2543 2550 2551 2552 2043 2045 2046 13 | 107 1604 677 348 988 479 1573 553 561 573 1611 590 1620 1628 1117 1146 1158 1160 1199 1315 1345 1347 1377 1404 896 1423 424 453 975 469 995 997 486 1019 515 518 520 1036 525 1550 532 536 538 1563 1571 1060 1064 1067 1580 605 435 1589 1078 1591 1594 1596 1088 579 1605 524 1100 1613 592 598 1636 606 1125 609 443 1124 1637 1126 1127 1641 1130 1643 1645 625 617 1659 1148 1663 641 1671 1673 1165 1680 465 1172 1177 666 1179 1703 1192 1715 1716 1205 1722 1214 1728 1730 1219 1732 1222 1735 1227 1741 1743 1238 1751 1235 1255 1771 1261 1263 1267 1767 1273 1786 1279 1796 1799 1288 1289 1291 1298 559 1820 1824 1827 1317 1830 1833 1322 1323 402 1849 566 1359 1877 1881 1882 1372 1082 357 1898 1899 1388 365 1390 369 370 375 381 897 902 1857 1419 397 408 1426 1432 580 1436 929 932 497 498 596 947 436 437 1462 439 440 1467 444 957 1470 1471 960 1473 1099 417 456 680 460 395 1488 612 467 505 1871 989 459 484 487 1512 423 494 509 1520 1009 1522 1017 1018 1021 14 | 686 805 705 711 719 747 814 842 713 724 741 762 780 781 810 818 823 687 726 728 736 738 745 766 797 800 813 817 828 829 833 834 847 853 688 691 694 703 706 708 710 718 727 734 760 763 771 772 777 782 784 787 795 796 799 806 815 816 824 838 840 848 850 856 695 697 698 701 707 715 717 731 732 733 737 746 752 753 754 756 761 764 770 773 774 783 785 786 788 794 802 803 804 819 821 825 826 827 830 832 835 836 837 844 846 849 854 860 864 893 15 | 1680 348 350 369 373 475 513 524 547 1631 651 107 626 1655 1155 654 1179 680 684 1205 1767 455 1323 1864 366 1432 1459 483 1512 1001 504 596 506 1024 515 598 1543 520 1547 525 1039 1554 1049 1563 541 544 1059 548 551 553 1067 557 558 1078 1591 569 1082 1596 1085 1086 523 1604 438 1100 1102 1618 1619 1108 1621 1622 601 1115 604 1630 609 1125 1647 1649 630 1146 1149 1665 644 1673 1163 1173 1175 1689 1180 1184 994 1700 1702 1703 1707 1198 1204 1718 1211 1056 1730 1221 1222 1742 433 1315 1238 1768 1774 1777 1794 1273 1277 1790 384 1282 1153 1289 1807 367 1813 623 1820 1826 1827 1316 561 1330 1842 1331 1903 566 482 1358 1873 1877 1370 352 353 489 1382 359 1384 361 362 1388 1391 368 371 374 1399 376 896 388 389 391 1516 395 397 400 408 403 916 407 412 419 420 934 936 428 430 432 1457 1460 1462 440 442 1467 445 1525 449 450 967 456 457 458 461 463 465 980 469 471 527 1506 484 1190 490 491 1004 1518 495 496 1192 579 1013 503 507 16 | 107 1211 397 418 481 517 524 3842 531 541 555 1078 3683 1147 646 1388 3722 651 1172 1804 1813 1158 1335 1352 1886 1376 1898 376 1173 1604 423 1017 1530 513 514 515 519 3592 522 3599 528 1554 1047 1796 1563 547 3505 556 1582 559 1589 1080 1085 1086 1603 1092 3437 1608 586 587 1102 921 1106 1620 1622 1112 1892 604 1117 1124 1126 1137 1299 1140 1146 1148 1149 3710 1665 1669 1160 1164 1165 1167 1902 1049 1185 1734 1191 683 1216 1198 1199 3553 3501 1728 1730 1222 1314 1827 1751 1239 1754 1243 1402 1758 1317 1766 1767 3817 1261 3824 3830 1271 1272 981 1280 3456 1794 1284 1798 1799 1800 3850 1291 1752 1298 1811 3802 1157 1312 1826 1315 1318 3948 1837 1839 1845 1846 1850 1857 993 1866 1847 1358 568 995 1879 348 1338 1380 1893 361 363 1900 1389 366 367 1395 1909 967 916 378 896 902 391 392 394 398 400 404 405 1431 917 1435 416 420 3526 1351 428 429 430 432 433 946 947 436 952 1470 3521 966 455 1480 463 1870 982 1495 1496 987 1502 483 1509 489 3563 1516 1006 1520 504 17 | -------------------------------------------------------------------------------- /com2.txt: -------------------------------------------------------------------------------- 1 | 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 2704 2740 2814 2838 2885 3003 3290 2 | 34 173 198 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 579 580 581 584 585 586 587 588 589 590 591 592 593 597 598 601 602 603 604 605 606 607 608 609 610 611 612 613 614 616 617 618 619 620 621 622 623 624 625 626 629 630 631 633 634 636 637 638 639 641 642 644 645 646 648 649 651 652 653 654 655 656 657 660 663 664 666 667 668 669 671 672 673 674 676 677 678 679 680 682 683 684 685 905 910 912 915 930 931 936 939 948 962 976 1001 1012 1013 1025 1064 1080 1094 1095 1113 1114 1122 1142 1154 1155 1162 1179 1210 1212 1223 1234 1237 1263 1294 1295 1300 1313 1318 1320 1332 1342 1349 1356 1357 1358 1372 1374 1397 1400 1408 1422 1425 1427 1430 1443 1446 1453 1455 1464 1478 1487 1490 1497 1506 1512 1514 1526 1543 1545 1549 1565 1566 1574 1599 1602 1606 1616 1625 1631 1645 1654 1657 1660 1671 1673 1677 1679 1692 1694 1695 1696 1731 1740 1743 1744 1762 1763 1766 1767 1773 1776 1777 1778 1790 1801 1808 1814 1829 1831 1844 1853 1870 1878 1881 1882 1901 1910 1967 3 | 107 897 899 904 906 907 911 916 918 920 921 922 925 926 927 932 934 946 947 952 953 959 960 961 966 967 972 973 978 980 982 983 991 993 995 996 997 998 999 1002 1003 1004 1006 1008 1017 1022 1024 1026 1027 1028 1029 1034 1038 1039 1040 1046 1047 1048 1049 1050 1054 1055 1056 1058 1059 1063 1068 1069 1071 1074 1075 1076 1077 1078 1079 1083 1084 1086 1087 1091 1092 1096 1101 1105 1107 1110 1112 1116 1117 1119 1123 1124 1125 1126 1128 1130 1132 1133 1135 1145 1146 1149 1150 1151 1153 1156 1157 1160 1161 1163 1164 1171 1172 1173 1175 1180 1181 1182 1184 1185 1191 1194 1195 1196 1198 1199 1200 1201 1202 1205 1206 1207 1208 1209 1211 1214 1219 1220 1222 1224 1230 1238 1239 1242 1243 1250 1253 1255 1256 1262 1265 1266 1267 1269 1270 1271 1272 1276 1278 1280 1283 1285 1287 1288 1289 1290 1291 1293 1297 1301 1302 1305 1312 1321 1323 1326 1327 1329 1330 1331 1334 1335 1336 1337 1339 1340 1341 1344 1346 1351 1352 1359 1361 1365 1367 1369 1370 1375 1376 1377 1380 1383 1386 1387 1388 1389 1390 1391 1393 1395 1398 1399 1401 1402 1405 1407 1409 1411 1416 1419 1420 1431 1434 1437 1439 1440 1441 1442 1447 1449 1456 1457 1458 1460 1461 1463 1466 1467 1470 1471 1476 1477 1480 1483 1484 1485 1486 1488 1491 1498 1501 1509 1513 1516 1517 1518 1519 1520 1521 1522 1523 1524 1527 1528 1530 1532 1535 1537 1538 1539 1540 1542 1544 1547 1551 1552 1554 1556 1557 1558 1559 1560 1561 1563 1564 1570 1571 1572 1580 1581 1584 1585 1586 1589 1590 1594 1597 1598 1600 1603 1604 1605 1608 1609 1610 1612 1613 1614 1617 1618 1619 1620 1621 1622 1623 1626 1627 1632 1637 1639 1643 1644 1651 1652 1653 1659 1662 1663 1665 1668 1669 1670 1674 1675 1678 1683 1685 1688 1689 1697 1700 1702 1707 1710 1711 1712 1714 1717 1721 1722 1723 1724 1730 1734 1735 1736 1737 1741 1746 1750 1752 1753 1754 1757 1761 1765 1768 1769 1771 1772 1774 1775 1779 1780 1782 1789 1791 1792 1793 1795 1796 1797 1799 1800 1804 1805 1809 1810 1811 1813 1816 1817 1819 1821 1822 1823 1826 1827 1832 1833 1834 1835 1836 1839 1842 1843 1845 1846 1849 1851 1852 1854 1856 1858 1860 1861 1863 1864 1865 1867 1868 1874 1877 1879 1883 1886 1888 1890 1891 1898 1900 1902 1905 1909 1911 4 | 136 1465 1577 1718 1912 1913 1915 1916 1920 1922 1923 1926 1930 1932 1933 1934 1936 1937 1939 1940 1941 1944 1945 1947 1948 1949 1952 1954 1955 1956 1959 1960 1961 1964 1965 1969 1970 1974 1975 1977 1978 1981 1982 1987 1990 1994 1996 2002 2003 2007 2008 2010 2015 2017 2023 2026 2028 2029 2031 2032 2034 2035 2036 2038 2039 2041 2042 2047 2048 2049 2051 2052 2053 2054 2057 2062 2065 2066 2068 2070 2071 2072 2075 2079 2080 2081 2082 2085 2087 2089 2091 2094 2096 2099 2100 2101 2102 2105 2106 2107 2110 2111 2113 2114 2117 2119 2120 2125 2126 2127 2128 2129 2130 2132 2133 2134 2135 2137 2138 2141 2143 2144 2145 2146 2148 2149 2151 2152 2153 2155 2156 2158 2160 2161 2162 2163 2166 2167 2168 2169 2173 2174 2175 2176 2177 2178 2180 2181 2182 2183 2185 2186 2187 2189 2191 2193 2194 2195 2196 2197 2198 2199 2203 2204 2205 2208 2209 2211 2214 2215 2217 2219 2221 2222 2223 2224 2226 2227 2228 2230 2231 2232 2235 2236 2238 2239 2241 2242 2243 2245 2246 2247 2248 2249 2250 2251 2252 2254 2256 2259 2260 2263 2264 2265 2267 2268 2269 2270 2272 2273 2274 2277 2279 2280 2281 2282 2283 2285 2286 2289 2291 2292 2293 2294 2295 2296 2298 2301 2302 2304 2305 2310 2312 2313 2314 2315 2316 2317 2318 2319 2321 2322 2325 2327 2328 2330 2332 2333 2335 2336 2338 2341 2342 2343 2344 2345 2351 2355 2357 2358 2360 2361 2362 2366 2367 2368 2371 2372 2373 2375 2377 2379 2380 2382 2384 2385 2388 2389 2390 2394 2396 2397 2398 2399 2400 2402 2403 2405 2406 2411 2412 2413 2415 2416 2417 2419 2420 2421 2422 2424 2425 2431 2432 2434 2436 2437 2438 2439 2440 2441 2445 2448 2450 2451 2452 2453 2455 2456 2457 2458 2461 2463 2465 2470 2471 2472 2473 2474 2475 2479 2480 2481 2483 2487 2488 2491 2493 2496 2497 2498 2501 2502 2503 2505 2508 2509 2510 2511 2512 2514 2516 2517 2518 2519 2522 2523 2525 2527 2528 2529 2530 2531 2533 2534 2535 2537 2540 2541 2544 2547 2548 2555 2557 2558 2562 2565 2566 2567 2568 2569 2570 2571 2572 2577 2580 2581 2582 2587 2588 2589 2592 2594 2595 2596 2597 2598 2599 2603 2605 2608 2609 2610 2612 2614 2616 2617 2618 2621 2628 2629 2632 2634 2635 2637 2641 2642 2643 2644 2649 2651 2652 2653 2656 5 | 389 596 896 898 902 913 914 917 919 924 933 935 937 940 941 942 943 944 945 949 954 957 964 965 968 969 971 974 975 977 981 984 985 986 988 989 994 1000 1005 1007 1009 1010 1011 1014 1015 1018 1019 1020 1021 1032 1035 1036 1037 1041 1042 1044 1051 1052 1060 1062 1065 1066 1067 1070 1072 1081 1082 1089 1090 1093 1099 1100 1102 1103 1104 1106 1108 1109 1115 1120 1121 1127 1129 1131 1134 1139 1141 1143 1147 1148 1152 1158 1159 1166 1167 1168 1169 1170 1176 1183 1188 1189 1190 1192 1204 1213 1215 1217 1221 1225 1226 1227 1228 1229 1231 1233 1235 1236 1241 1244 1245 1246 1247 1248 1249 1257 1258 1259 1260 1261 1264 1268 1273 1275 1277 1279 1281 1282 1284 1286 1292 1296 1298 1299 1303 1304 1306 1308 1309 1310 1311 1315 1316 1322 1324 1338 1343 1345 1347 1348 1350 1354 1355 1362 1364 1366 1373 1379 1381 1382 1385 1392 1396 1404 1406 1412 1413 1414 1415 1417 1418 1423 1426 1428 1429 1432 1435 1436 1438 1444 1448 1451 1454 1459 1462 1469 1472 1473 1474 1475 1479 1481 1482 1489 1492 1495 1496 1500 1502 1503 1504 1507 1508 1510 1515 1525 1529 1531 1536 1541 1550 1569 1573 1575 1576 1578 1582 1583 1587 1591 1592 1595 1596 1611 1624 1628 1630 1633 1634 1635 1636 1638 1640 1646 1647 1649 1650 1655 1658 1661 1664 1667 1672 1680 1681 1682 1686 1690 1691 1699 1703 1704 1706 1709 1715 1716 1720 1725 1727 1728 1729 1732 1738 1739 1742 1745 1748 1749 1751 1755 1756 1759 1764 1770 1781 1783 1785 1786 1794 1807 1812 1815 1818 1820 1824 1828 1830 1840 1841 1847 1848 1850 1855 1857 1859 1862 1869 1871 1872 1873 1875 1876 1884 1885 1887 1893 1894 1897 1899 1903 1904 1906 1907 6 | 576 577 578 582 583 595 599 600 615 627 628 632 635 640 643 647 650 658 659 661 662 665 670 675 681 7 | 594 3980 3981 3982 3983 3984 3985 3986 3987 3988 3989 3990 3991 3992 3993 3994 3995 3996 3997 3998 3999 4000 4001 4002 4003 4004 4005 4006 4007 4008 4009 4010 4011 4012 4013 4014 4015 4016 4017 4018 4019 4020 4021 4022 4023 4024 4025 4026 4027 4028 4029 4030 4031 4032 4033 4034 4035 4036 4037 4038 8 | 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 843 844 845 846 847 848 849 850 851 852 853 854 855 856 858 859 860 861 863 864 866 867 869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 892 893 894 895 9 | 857 862 865 868 1085 3437 3438 3439 3440 3441 3442 3443 3444 3445 3446 3447 3448 3449 3450 3451 3452 3453 3454 3455 3456 3457 3458 3459 3460 3461 3462 3463 3464 3465 3466 3467 3468 3469 3470 3471 3472 3473 3474 3475 3476 3477 3478 3479 3480 3481 3482 3483 3484 3485 3486 3487 3488 3489 3490 3491 3492 3493 3494 3495 3496 3497 3498 3499 3500 3501 3502 3503 3504 3505 3506 3507 3508 3509 3510 3511 3512 3513 3514 3515 3516 3517 3518 3519 3520 3521 3522 3523 3524 3525 3526 3527 3528 3529 3530 3531 3532 3533 3534 3535 3536 3537 3538 3539 3540 3541 3542 3543 3544 3545 3546 3547 3548 3549 3550 3551 3552 3553 3554 3555 3556 3557 3558 3559 3560 3561 3562 3563 3564 3565 3566 3567 3568 3569 3570 3571 3572 3573 3574 3575 3576 3577 3578 3579 3580 3581 3582 3583 3584 3585 3586 3587 3588 3589 3590 3591 3592 3593 3594 3595 3596 3597 3598 3599 3600 3601 3602 3603 3604 3605 3606 3607 3608 3609 3610 3611 3612 3613 3614 3615 3616 3617 3618 3619 3620 3621 3622 3623 3624 3625 3626 3627 3628 3629 3630 3631 3632 3633 3634 3635 3636 3637 3638 3639 3640 3641 3642 3643 3644 3645 3646 3647 3648 3649 3650 3651 3652 3653 3654 3655 3656 3657 3658 3659 3660 3661 3662 3663 3664 3665 3666 3667 3668 3669 3670 3671 3672 3673 3674 3675 3676 3677 3678 3679 3680 3681 3682 3683 3684 3685 3686 3687 3688 3689 3690 3691 3692 3693 3694 3695 3696 3697 3698 3699 3700 3701 3702 3703 3704 3705 3706 3707 3708 3709 3710 3711 3712 3713 3714 3715 3716 3717 3718 3719 3720 3721 3722 3723 3724 3725 3726 3727 3728 3729 3730 3731 3732 3733 3734 3735 3736 3737 3738 3739 3740 3741 3742 3743 3744 3745 3746 3747 3748 3749 3750 3751 3752 3753 3754 3755 3756 3757 3758 3759 3760 3761 3762 3763 3764 3765 3766 3767 3768 3769 3770 3771 3772 3773 3774 3775 3776 3777 3778 3779 3780 3781 3782 3783 3784 3785 3786 3787 3788 3789 3790 3791 3792 3793 3794 3795 3796 3797 3798 3799 3800 3801 3802 3803 3804 3805 3806 3807 3808 3809 3810 3811 3812 3813 3814 3815 3816 3817 3818 3819 3820 3821 3822 3823 3824 3825 3826 3827 3828 3829 3830 3831 3832 3833 3834 3835 3836 3837 3838 3839 3840 3841 3842 3843 3844 3845 3846 3847 3848 3849 3850 3851 3852 3853 3854 3855 3856 3857 3858 3859 3860 3861 3862 3863 3864 3865 3866 3867 3868 3869 3870 3871 3872 3873 3874 3875 3876 3877 3878 3879 3880 3881 3882 3883 3884 3885 3886 3887 3888 3889 3890 3891 3892 3893 3894 3895 3896 3897 3898 3899 3900 3901 3902 3903 3904 3905 3906 3907 3908 3909 3910 3911 3912 3913 3914 3915 3916 3917 3918 3919 3920 3921 3922 3923 3924 3925 3926 3927 3928 3929 3930 3931 3932 3933 3934 3935 3936 3937 3938 3939 3940 3941 3942 3943 3944 3945 3946 3947 3948 3949 3950 3951 3952 3953 3954 3955 3956 3957 3958 3959 3960 3961 3962 3963 3964 3965 3966 3967 3968 3969 3970 3971 3972 3973 3974 3975 3976 3977 3978 3979 10 | 900 901 903 908 909 923 928 929 938 950 951 955 956 958 963 970 979 987 992 1016 1023 1030 1031 1033 1043 1045 1053 1057 1061 1073 1088 1097 1098 1111 1118 1136 1137 1138 1144 1165 1174 1177 1178 1186 1187 1193 1197 1203 1216 1218 1232 1240 1251 1252 1254 1274 1307 1314 1317 1319 1325 1328 1333 1353 1360 1363 1368 1371 1378 1384 1394 1403 1410 1421 1424 1433 1445 1452 1468 1493 1494 1499 1511 1533 1546 1548 1553 1555 1562 1567 1568 1579 1588 1593 1601 1607 1615 1629 1641 1648 1676 1687 1693 1698 1701 1705 1708 1713 1719 1733 1747 1760 1784 1787 1788 1798 1802 1803 1806 1825 1837 1838 1866 1880 1889 1892 1895 1896 1908 11 | 990 1140 1450 1505 1534 1642 1656 1666 1684 1726 1758 2661 2662 2665 2667 2668 2671 2672 2674 2675 2677 2678 2682 2684 2685 2686 2687 2688 2690 2691 2692 2693 2696 2697 2700 2701 2702 2707 2708 2709 2711 2712 2713 2714 2715 2716 2718 2719 2720 2721 2722 2723 2724 2725 2727 2728 2730 2732 2733 2734 2735 2736 2737 2739 2742 2744 2747 2748 2751 2752 2753 2758 2760 2762 2764 2765 2768 2770 2771 2772 2774 2775 2776 2778 2779 2780 2781 2783 2784 2785 2787 2788 2790 2791 2792 2793 2795 2796 2797 2799 2801 2802 2803 2804 2805 2808 2811 2812 2813 2816 2817 2818 2819 2820 2822 2823 2824 2825 2826 2828 2829 2830 2831 2832 2833 2836 2840 2841 2842 2843 2844 2845 2846 2847 2848 2849 2853 2855 2856 2857 2858 2859 2860 2861 2863 2865 2868 2869 2870 2871 2875 2876 2878 2881 2882 2883 2886 2892 2893 2894 2895 2898 2900 2901 2902 2903 2904 2914 2918 2921 2922 2923 2926 2930 2932 2933 2934 2935 2936 2937 2938 2939 2941 2942 2947 2949 2950 2951 2954 2955 2957 2958 2962 2964 2965 2968 2970 2971 2973 2975 2976 2979 2980 2983 2984 2989 2990 2991 2992 2995 2996 2997 2998 2999 3001 3002 3005 3006 3007 3009 3011 3012 3013 3015 3016 3017 3018 3019 3020 3021 3024 3027 3028 3031 3032 3034 3035 3037 3038 3039 3041 3043 3046 3048 3049 3050 3051 3053 3055 3057 3058 3060 3061 3062 3063 3066 3067 3068 3069 3070 3071 3072 3074 3075 3076 3077 3079 3081 3082 3083 3086 3087 3088 3089 3091 3093 3094 3097 3099 3100 3101 3102 3103 3105 3106 3107 3109 3110 3111 3112 3113 3115 3116 3118 3120 3121 3122 3124 3125 3126 3127 3128 3131 3132 3133 3134 3136 3138 3140 3141 3144 3145 3146 3147 3148 3150 3151 3152 3154 3155 3156 3157 3159 3160 3161 3162 3164 3165 3167 3168 3169 3170 3171 3172 3173 3175 3176 3177 3178 3179 3180 3181 3182 3183 3185 3186 3187 3188 3189 3190 3191 3192 3194 3196 3197 3198 3199 3200 3201 3203 3204 3205 3206 3207 3208 3209 3210 3211 3212 3214 3215 3216 3217 3218 3219 3221 3222 3223 3224 3225 3230 3233 3235 3236 3239 3240 3241 3242 3244 3245 3246 3247 3248 3249 3251 3252 3253 3254 3255 3256 3258 3259 3260 3261 3262 3263 3264 3265 3266 3267 3268 3270 3271 3272 3273 3274 3275 3276 3277 3278 3279 3280 3281 3282 3285 3287 3288 3289 3291 3292 3295 3296 3297 3298 3299 3302 3305 3306 3307 3308 3310 3312 3313 3316 3319 3320 3321 3322 3323 3326 3327 3328 3329 3330 3331 3332 3334 3335 3336 3337 3338 3340 3341 3342 3343 3344 3345 3346 3347 3349 3350 3351 3352 3353 3354 3355 3356 3357 3358 3359 3360 3361 3365 3366 3367 3368 3369 3370 3371 3372 3373 3374 3375 3376 3377 3378 3379 3380 3381 3383 3385 3386 3389 3391 3393 3395 3396 3397 3398 3401 3402 3403 3404 3405 3406 3407 3408 3409 3410 3411 3412 3415 3416 3418 3419 3420 3421 3422 3424 3425 3427 3429 3431 3432 3434 12 | 1914 1919 1921 1927 1928 1931 1935 1950 1951 1957 1958 1968 1972 1973 1976 1980 1991 1995 1998 1999 2000 2001 2004 2006 2009 2011 2012 2013 2014 2016 2018 2022 2024 2025 2027 2050 2061 2097 2116 2157 2159 2171 2192 2202 2225 2284 2297 2337 2346 2364 2365 2378 2435 2447 2459 2494 2538 2583 2585 2620 2626 2627 2633 2636 2639 2640 2645 2647 2648 2657 2658 2659 2660 13 | 1917 1918 1925 1929 1938 1942 1943 1946 1953 1962 1963 1966 1971 1979 1983 1984 1985 1986 1989 1993 1997 2005 2019 2020 2021 2030 2033 2037 2040 2043 2044 2045 2046 2055 2056 2058 2059 2060 2063 2064 2067 2069 2073 2074 2077 2078 2083 2084 2086 2088 2090 2092 2093 2095 2098 2103 2104 2108 2109 2112 2115 2118 2121 2122 2123 2124 2131 2136 2139 2140 2142 2147 2150 2154 2164 2165 2172 2179 2184 2188 2190 2200 2201 2206 2210 2212 2213 2216 2218 2220 2229 2233 2234 2237 2240 2244 2253 2257 2258 2261 2262 2266 2271 2275 2276 2278 2287 2288 2290 2299 2300 2303 2306 2307 2308 2309 2311 2323 2324 2326 2329 2331 2334 2339 2340 2347 2348 2350 2352 2353 2354 2356 2359 2363 2369 2370 2374 2376 2381 2383 2386 2387 2391 2392 2393 2395 2404 2407 2408 2409 2410 2414 2418 2423 2427 2428 2429 2430 2433 2442 2443 2446 2449 2454 2460 2462 2464 2467 2468 2469 2476 2477 2478 2482 2484 2485 2489 2492 2495 2499 2500 2504 2506 2507 2520 2521 2524 2526 2532 2536 2539 2542 2543 2545 2546 2549 2550 2551 2552 2553 2554 2556 2559 2560 2561 2563 2564 2573 2574 2575 2576 2578 2579 2584 2586 2590 2591 2593 2600 2601 2602 2604 2606 2607 2611 2613 2615 2619 2623 2624 2625 2630 2631 2638 2646 2654 2655 14 | 1924 1988 1992 2076 2170 2207 2255 2320 2349 2401 2426 2444 2466 2486 2490 2513 2515 2622 2650 15 | 2663 2664 2666 2669 2673 2676 2679 2680 2681 2683 2689 2694 2695 2698 2705 2706 2710 2717 2726 2729 2731 2738 2741 2743 2745 2746 2749 2750 2754 2755 2756 2757 2759 2761 2763 2766 2769 2773 2777 2782 2786 2789 2794 2798 2800 2806 2807 2809 2810 2815 2821 2827 2835 2837 2839 2850 2851 2852 2854 2862 2864 2866 2867 2872 2873 2874 2877 2880 2884 2887 2888 2890 2891 2896 2897 2899 2905 2906 2907 2908 2909 2910 2911 2912 2913 2915 2916 2917 2919 2920 2924 2925 2927 2928 2929 2931 2940 2943 2944 2945 2946 2948 2952 2953 2956 2960 2961 2963 2966 2967 2969 2974 2977 2978 2981 2985 2986 2987 2988 2993 2994 3000 3004 3010 3014 3022 3023 3025 3026 3029 3030 3033 3036 3040 3042 3044 3045 3047 3052 3054 3056 3059 3064 3065 3073 3078 3080 3084 3085 3090 3092 3095 3096 3098 3104 3108 3114 3117 3119 3123 3129 3130 3135 3137 3139 3142 3143 3149 3153 3158 3163 3166 3174 3184 3193 3195 3202 3213 3220 3226 3227 3228 3229 3231 3232 3234 3237 3238 3243 3250 3257 3269 3284 3286 3293 3294 3300 3301 3303 3304 3315 3317 3324 3333 3339 3348 3362 3363 3364 3384 3387 3388 3390 3392 3394 3399 3400 3413 3414 3417 3426 3428 3430 3433 3435 3436 16 | 2670 2699 2703 2767 2834 2879 2889 2959 2972 2982 3008 3283 3309 3311 3314 3318 3325 3382 3423 17 | --------------------------------------------------------------------------------