Skip to contents

Get a dataframe of multiple draws of only the predicted elements of the nowcast vector

Usage

get_nowcast_pred_draws(
  point_nowcast_matrix,
  reporting_triangle,
  dispersion,
  draws = 1000,
  fun_to_aggregate = sum,
  k = 1
)

Arguments

point_nowcast_matrix

Matrix of point nowcast predictions and observations, with rows representing the reference times and columns representing the delays.

reporting_triangle

Matrix of the reporting triangle, with rows representing the time points of reference and columns representing the delays. Can be a reporting matrix or incomplete reporting matrix. Can also be a ragged reporting triangle, where multiple columns are reported for the same row. (e.g. weekly reporting of daily data).

dispersion

Vector of dispersion parameters indexed by horizon from minus one to the maximum delay.

draws

Integer indicating the number of draws of the predicted nowcast vector to generate. Default is 1000.

fun_to_aggregate

Function that will operate along the nowcast vectors after summing across delays. Eventually, we can add things like mean, but for now since we are only providing a negative binomial observation model, we can only allow sum. Currently supported functions: sum.

k

Integer indicating the number of reference times to apply the fun_to_aggregate over to create target used to compute the nowcast errors.

Value

Dataframe containing the predicted point nowcast vectors indexed by reference time (pred_count), reference time (time), and the draw index (draw).

Examples

point_nowcast_matrix <- matrix(
  c(
    80, 50, 25, 10,
    100, 50, 30, 20,
    90, 45, 25, 16.8,
    80, 40, 21.2, 19.5,
    70, 34.5, 15.4, 9.1
  ),
  nrow = 5,
  byrow = TRUE
)
reporting_triangle <- generate_triangle(point_nowcast_matrix)
disp <- c(0.8, 12.4, 9.1)
nowcast_pred_draws <- get_nowcast_pred_draws(
  point_nowcast_matrix,
  reporting_triangle,
  disp,
  draws = 5
)
nowcast_pred_draws
#>    pred_count time draw
#> 1           0    1    1
#> 2           0    2    1
#> 3          21    3    1
#> 4          49    4    1
#> 5          20    5    1
#> 6           0    1    2
#> 7           0    2    2
#> 8          15    3    2
#> 9          43    4    2
#> 10          2    5    2
#> 11          0    1    3
#> 12          0    2    3
#> 13         11    3    3
#> 14         59    4    3
#> 15         69    5    3
#> 16          0    1    4
#> 17          0    2    4
#> 18         18    3    4
#> 19         58    4    4
#> 20         17    5    4
#> 21          0    1    5
#> 22          0    2    5
#> 23         11    3    5
#> 24         41    4    5
#> 25         11    5    5
# Get nowcast pred draws over rolling sum
nowcast_pred_draws_rolling_df <- get_nowcast_pred_draws(
  point_nowcast_matrix,
  reporting_triangle,
  disp,
  500,
  fun_to_aggregate = sum,
  k = 2
)
nowcast_pred_draws_rolling_df
#>      pred_count time draw
#> 1             0    1    1
#> 2             0    2    1
#> 3            16    3    1
#> 4            53    4    1
#> 5           195    5    1
#> 6             0    1    2
#> 7             0    2    2
#> 8             9    3    2
#> 9            68    4    2
#> 10           63    5    2
#> 11            0    1    3
#> 12            0    2    3
#> 13           15    3    3
#> 14           47    4    3
#> 15            2    5    3
#> 16            0    1    4
#> 17            0    2    4
#> 18           22    3    4
#> 19           53    4    4
#> 20           51    5    4
#> 21            0    1    5
#> 22            0    2    5
#> 23            8    3    5
#> 24           30    4    5
#> 25           45    5    5
#> 26            0    1    6
#> 27            0    2    6
#> 28           14    3    6
#> 29           33    4    6
#> 30           18    5    6
#> 31            0    1    7
#> 32            0    2    7
#> 33           23    3    7
#> 34           20    4    7
#> 35           40    5    7
#> 36            0    1    8
#> 37            0    2    8
#> 38            8    3    8
#> 39           59    4    8
#> 40           13    5    8
#> 41            0    1    9
#> 42            0    2    9
#> 43           26    3    9
#> 44           37    4    9
#> 45          144    5    9
#> 46            0    1   10
#> 47            0    2   10
#> 48           24    3   10
#> 49           38    4   10
#> 50            1    5   10
#> 51            0    1   11
#> 52            0    2   11
#> 53           21    3   11
#> 54           59    4   11
#> 55          372    5   11
#> 56            0    1   12
#> 57            0    2   12
#> 58           28    3   12
#> 59           44    4   12
#> 60           17    5   12
#> 61            0    1   13
#> 62            0    2   13
#> 63           16    3   13
#> 64           29    4   13
#> 65           32    5   13
#> 66            0    1   14
#> 67            0    2   14
#> 68           10    3   14
#> 69           19    4   14
#> 70           39    5   14
#> 71            0    1   15
#> 72            0    2   15
#> 73           25    3   15
#> 74           34    4   15
#> 75            3    5   15
#> 76            0    1   16
#> 77            0    2   16
#> 78           11    3   16
#> 79           34    4   16
#> 80           17    5   16
#> 81            0    1   17
#> 82            0    2   17
#> 83           18    3   17
#> 84           36    4   17
#> 85           16    5   17
#> 86            0    1   18
#> 87            0    2   18
#> 88           31    3   18
#> 89           11    4   18
#> 90           52    5   18
#> 91            0    1   19
#> 92            0    2   19
#> 93           10    3   19
#> 94           38    4   19
#> 95           15    5   19
#> 96            0    1   20
#> 97            0    2   20
#> 98           14    3   20
#> 99           13    4   20
#> 100          18    5   20
#> 101           0    1   21
#> 102           0    2   21
#> 103          15    3   21
#> 104          30    4   21
#> 105          59    5   21
#> 106           0    1   22
#> 107           0    2   22
#> 108          21    3   22
#> 109          73    4   22
#> 110          72    5   22
#> 111           0    1   23
#> 112           0    2   23
#> 113          17    3   23
#> 114          45    4   23
#> 115         105    5   23
#> 116           0    1   24
#> 117           0    2   24
#> 118          24    3   24
#> 119          15    4   24
#> 120          11    5   24
#> 121           0    1   25
#> 122           0    2   25
#> 123          35    3   25
#> 124          32    4   25
#> 125         140    5   25
#> 126           0    1   26
#> 127           0    2   26
#> 128          16    3   26
#> 129          24    4   26
#> 130          58    5   26
#> 131           0    1   27
#> 132           0    2   27
#> 133          19    3   27
#> 134          35    4   27
#> 135           4    5   27
#> 136           0    1   28
#> 137           0    2   28
#> 138          24    3   28
#> 139          53    4   28
#> 140           4    5   28
#> 141           0    1   29
#> 142           0    2   29
#> 143          13    3   29
#> 144          54    4   29
#> 145           2    5   29
#> 146           0    1   30
#> 147           0    2   30
#> 148          23    3   30
#> 149          30    4   30
#> 150          24    5   30
#> 151           0    1   31
#> 152           0    2   31
#> 153          14    3   31
#> 154          24    4   31
#> 155          64    5   31
#> 156           0    1   32
#> 157           0    2   32
#> 158          17    3   32
#> 159          28    4   32
#> 160         190    5   32
#> 161           0    1   33
#> 162           0    2   33
#> 163          16    3   33
#> 164          36    4   33
#> 165         224    5   33
#> 166           0    1   34
#> 167           0    2   34
#> 168           4    3   34
#> 169          51    4   34
#> 170          13    5   34
#> 171           0    1   35
#> 172           0    2   35
#> 173          13    3   35
#> 174          40    4   35
#> 175         103    5   35
#> 176           0    1   36
#> 177           0    2   36
#> 178          12    3   36
#> 179          39    4   36
#> 180          86    5   36
#> 181           0    1   37
#> 182           0    2   37
#> 183          27    3   37
#> 184          24    4   37
#> 185         114    5   37
#> 186           0    1   38
#> 187           0    2   38
#> 188          10    3   38
#> 189          27    4   38
#> 190         134    5   38
#> 191           0    1   39
#> 192           0    2   39
#> 193           9    3   39
#> 194          46    4   39
#> 195           2    5   39
#> 196           0    1   40
#> 197           0    2   40
#> 198          30    3   40
#> 199          39    4   40
#> 200          70    5   40
#> 201           0    1   41
#> 202           0    2   41
#> 203          14    3   41
#> 204          36    4   41
#> 205          21    5   41
#> 206           0    1   42
#> 207           0    2   42
#> 208          27    3   42
#> 209          29    4   42
#> 210          53    5   42
#> 211           0    1   43
#> 212           0    2   43
#> 213          12    3   43
#> 214          36    4   43
#> 215           6    5   43
#> 216           0    1   44
#> 217           0    2   44
#> 218          28    3   44
#> 219          35    4   44
#> 220          39    5   44
#> 221           0    1   45
#> 222           0    2   45
#> 223          17    3   45
#> 224          32    4   45
#> 225           9    5   45
#> 226           0    1   46
#> 227           0    2   46
#> 228          16    3   46
#> 229          50    4   46
#> 230         141    5   46
#> 231           0    1   47
#> 232           0    2   47
#> 233          16    3   47
#> 234          52    4   47
#> 235          21    5   47
#> 236           0    1   48
#> 237           0    2   48
#> 238          11    3   48
#> 239          39    4   48
#> 240           0    5   48
#> 241           0    1   49
#> 242           0    2   49
#> 243          16    3   49
#> 244          36    4   49
#> 245          21    5   49
#> 246           0    1   50
#> 247           0    2   50
#> 248          12    3   50
#> 249          54    4   50
#> 250          17    5   50
#> 251           0    1   51
#> 252           0    2   51
#> 253          13    3   51
#> 254          44    4   51
#> 255          97    5   51
#> 256           0    1   52
#> 257           0    2   52
#> 258          13    3   52
#> 259          66    4   52
#> 260          52    5   52
#> 261           0    1   53
#> 262           0    2   53
#> 263          21    3   53
#> 264          35    4   53
#> 265          20    5   53
#> 266           0    1   54
#> 267           0    2   54
#> 268          11    3   54
#> 269          58    4   54
#> 270         109    5   54
#> 271           0    1   55
#> 272           0    2   55
#> 273          16    3   55
#> 274          52    4   55
#> 275          35    5   55
#> 276           0    1   56
#> 277           0    2   56
#> 278           8    3   56
#> 279          47    4   56
#> 280         152    5   56
#> 281           0    1   57
#> 282           0    2   57
#> 283          21    3   57
#> 284          49    4   57
#> 285           9    5   57
#> 286           0    1   58
#> 287           0    2   58
#> 288          22    3   58
#> 289          60    4   58
#> 290          53    5   58
#> 291           0    1   59
#> 292           0    2   59
#> 293          10    3   59
#> 294          25    4   59
#> 295          20    5   59
#> 296           0    1   60
#> 297           0    2   60
#> 298          22    3   60
#> 299          32    4   60
#> 300         165    5   60
#> 301           0    1   61
#> 302           0    2   61
#> 303          23    3   61
#> 304          56    4   61
#> 305          33    5   61
#> 306           0    1   62
#> 307           0    2   62
#> 308          19    3   62
#> 309          45    4   62
#> 310          39    5   62
#> 311           0    1   63
#> 312           0    2   63
#> 313          25    3   63
#> 314          60    4   63
#> 315           6    5   63
#> 316           0    1   64
#> 317           0    2   64
#> 318          15    3   64
#> 319          36    4   64
#> 320           3    5   64
#> 321           0    1   65
#> 322           0    2   65
#> 323          15    3   65
#> 324          49    4   65
#> 325          61    5   65
#> 326           0    1   66
#> 327           0    2   66
#> 328          23    3   66
#> 329          37    4   66
#> 330          11    5   66
#> 331           0    1   67
#> 332           0    2   67
#> 333          16    3   67
#> 334          40    4   67
#> 335         148    5   67
#> 336           0    1   68
#> 337           0    2   68
#> 338          11    3   68
#> 339          37    4   68
#> 340          83    5   68
#> 341           0    1   69
#> 342           0    2   69
#> 343           6    3   69
#> 344          36    4   69
#> 345          87    5   69
#> 346           0    1   70
#> 347           0    2   70
#> 348          19    3   70
#> 349          35    4   70
#> 350         212    5   70
#> 351           0    1   71
#> 352           0    2   71
#> 353          17    3   71
#> 354          32    4   71
#> 355          31    5   71
#> 356           0    1   72
#> 357           0    2   72
#> 358          12    3   72
#> 359          53    4   72
#> 360          42    5   72
#> 361           0    1   73
#> 362           0    2   73
#> 363          22    3   73
#> 364          40    4   73
#> 365          37    5   73
#> 366           0    1   74
#> 367           0    2   74
#> 368          42    3   74
#> 369          27    4   74
#> 370           3    5   74
#> 371           0    1   75
#> 372           0    2   75
#> 373          18    3   75
#> 374          28    4   75
#> 375          21    5   75
#> 376           0    1   76
#> 377           0    2   76
#> 378          11    3   76
#> 379          45    4   76
#> 380         340    5   76
#> 381           0    1   77
#> 382           0    2   77
#> 383          21    3   77
#> 384          59    4   77
#> 385          24    5   77
#> 386           0    1   78
#> 387           0    2   78
#> 388          15    3   78
#> 389          37    4   78
#> 390          74    5   78
#> 391           0    1   79
#> 392           0    2   79
#> 393          23    3   79
#> 394          35    4   79
#> 395          38    5   79
#> 396           0    1   80
#> 397           0    2   80
#> 398          12    3   80
#> 399          56    4   80
#> 400         201    5   80
#> 401           0    1   81
#> 402           0    2   81
#> 403          20    3   81
#> 404          44    4   81
#> 405           9    5   81
#> 406           0    1   82
#> 407           0    2   82
#> 408          14    3   82
#> 409          40    4   82
#> 410          14    5   82
#> 411           0    1   83
#> 412           0    2   83
#> 413          11    3   83
#> 414          46    4   83
#> 415           4    5   83
#> 416           0    1   84
#> 417           0    2   84
#> 418          15    3   84
#> 419          63    4   84
#> 420          34    5   84
#> 421           0    1   85
#> 422           0    2   85
#> 423          16    3   85
#> 424          40    4   85
#> 425           1    5   85
#> 426           0    1   86
#> 427           0    2   86
#> 428           7    3   86
#> 429          46    4   86
#> 430           4    5   86
#> 431           0    1   87
#> 432           0    2   87
#> 433           7    3   87
#> 434          29    4   87
#> 435          97    5   87
#> 436           0    1   88
#> 437           0    2   88
#> 438          24    3   88
#> 439          58    4   88
#> 440         163    5   88
#> 441           0    1   89
#> 442           0    2   89
#> 443          18    3   89
#> 444          28    4   89
#> 445           3    5   89
#> 446           0    1   90
#> 447           0    2   90
#> 448          18    3   90
#> 449          49    4   90
#> 450         110    5   90
#> 451           0    1   91
#> 452           0    2   91
#> 453          25    3   91
#> 454          36    4   91
#> 455         140    5   91
#> 456           0    1   92
#> 457           0    2   92
#> 458          13    3   92
#> 459          25    4   92
#> 460           3    5   92
#> 461           0    1   93
#> 462           0    2   93
#> 463          22    3   93
#> 464          24    4   93
#> 465          20    5   93
#> 466           0    1   94
#> 467           0    2   94
#> 468           8    3   94
#> 469          47    4   94
#> 470          99    5   94
#> 471           0    1   95
#> 472           0    2   95
#> 473           9    3   95
#> 474          43    4   95
#> 475          65    5   95
#> 476           0    1   96
#> 477           0    2   96
#> 478           7    3   96
#> 479          44    4   96
#> 480          11    5   96
#> 481           0    1   97
#> 482           0    2   97
#> 483           9    3   97
#> 484          32    4   97
#> 485          50    5   97
#> 486           0    1   98
#> 487           0    2   98
#> 488          14    3   98
#> 489          34    4   98
#> 490          46    5   98
#> 491           0    1   99
#> 492           0    2   99
#> 493          24    3   99
#> 494          36    4   99
#> 495           6    5   99
#> 496           0    1  100
#> 497           0    2  100
#> 498          26    3  100
#> 499          54    4  100
#> 500           7    5  100
#> 501           0    1  101
#> 502           0    2  101
#> 503          10    3  101
#> 504          62    4  101
#> 505         198    5  101
#> 506           0    1  102
#> 507           0    2  102
#> 508          32    3  102
#> 509          44    4  102
#> 510          38    5  102
#> 511           0    1  103
#> 512           0    2  103
#> 513          17    3  103
#> 514          36    4  103
#> 515           6    5  103
#> 516           0    1  104
#> 517           0    2  104
#> 518          17    3  104
#> 519          39    4  104
#> 520         150    5  104
#> 521           0    1  105
#> 522           0    2  105
#> 523          17    3  105
#> 524          32    4  105
#> 525           1    5  105
#> 526           0    1  106
#> 527           0    2  106
#> 528          16    3  106
#> 529          63    4  106
#> 530           0    5  106
#> 531           0    1  107
#> 532           0    2  107
#> 533          15    3  107
#> 534          46    4  107
#> 535          12    5  107
#> 536           0    1  108
#> 537           0    2  108
#> 538          14    3  108
#> 539          53    4  108
#> 540         231    5  108
#> 541           0    1  109
#> 542           0    2  109
#> 543          10    3  109
#> 544          64    4  109
#> 545          89    5  109
#> 546           0    1  110
#> 547           0    2  110
#> 548          16    3  110
#> 549          35    4  110
#> 550          15    5  110
#> 551           0    1  111
#> 552           0    2  111
#> 553          14    3  111
#> 554          55    4  111
#> 555          11    5  111
#> 556           0    1  112
#> 557           0    2  112
#> 558          28    3  112
#> 559          38    4  112
#> 560           1    5  112
#> 561           0    1  113
#> 562           0    2  113
#> 563          25    3  113
#> 564          44    4  113
#> 565          52    5  113
#> 566           0    1  114
#> 567           0    2  114
#> 568           9    3  114
#> 569          21    4  114
#> 570           4    5  114
#> 571           0    1  115
#> 572           0    2  115
#> 573          15    3  115
#> 574          44    4  115
#> 575          26    5  115
#> 576           0    1  116
#> 577           0    2  116
#> 578           7    3  116
#> 579          48    4  116
#> 580         146    5  116
#> 581           0    1  117
#> 582           0    2  117
#> 583          16    3  117
#> 584          35    4  117
#> 585          50    5  117
#> 586           0    1  118
#> 587           0    2  118
#> 588          20    3  118
#> 589          55    4  118
#> 590          16    5  118
#> 591           0    1  119
#> 592           0    2  119
#> 593           6    3  119
#> 594          52    4  119
#> 595          22    5  119
#> 596           0    1  120
#> 597           0    2  120
#> 598          18    3  120
#> 599          42    4  120
#> 600          66    5  120
#> 601           0    1  121
#> 602           0    2  121
#> 603          18    3  121
#> 604          43    4  121
#> 605          64    5  121
#> 606           0    1  122
#> 607           0    2  122
#> 608          18    3  122
#> 609          58    4  122
#> 610          17    5  122
#> 611           0    1  123
#> 612           0    2  123
#> 613          20    3  123
#> 614          39    4  123
#> 615           5    5  123
#> 616           0    1  124
#> 617           0    2  124
#> 618          12    3  124
#> 619          77    4  124
#> 620           8    5  124
#> 621           0    1  125
#> 622           0    2  125
#> 623          26    3  125
#> 624          81    4  125
#> 625          19    5  125
#> 626           0    1  126
#> 627           0    2  126
#> 628          14    3  126
#> 629          89    4  126
#> 630          72    5  126
#> 631           0    1  127
#> 632           0    2  127
#> 633          23    3  127
#> 634          52    4  127
#> 635          72    5  127
#> 636           0    1  128
#> 637           0    2  128
#> 638           9    3  128
#> 639          53    4  128
#> 640          29    5  128
#> 641           0    1  129
#> 642           0    2  129
#> 643          14    3  129
#> 644          58    4  129
#> 645           1    5  129
#> 646           0    1  130
#> 647           0    2  130
#> 648          10    3  130
#> 649          64    4  130
#> 650          36    5  130
#> 651           0    1  131
#> 652           0    2  131
#> 653          24    3  131
#> 654          43    4  131
#> 655          50    5  131
#> 656           0    1  132
#> 657           0    2  132
#> 658          19    3  132
#> 659          39    4  132
#> 660          24    5  132
#> 661           0    1  133
#> 662           0    2  133
#> 663          12    3  133
#> 664          46    4  133
#> 665          15    5  133
#> 666           0    1  134
#> 667           0    2  134
#> 668          16    3  134
#> 669          51    4  134
#> 670          76    5  134
#> 671           0    1  135
#> 672           0    2  135
#> 673          12    3  135
#> 674          25    4  135
#> 675         100    5  135
#> 676           0    1  136
#> 677           0    2  136
#> 678          13    3  136
#> 679          38    4  136
#> 680          21    5  136
#> 681           0    1  137
#> 682           0    2  137
#> 683          21    3  137
#> 684          25    4  137
#> 685           7    5  137
#> 686           0    1  138
#> 687           0    2  138
#> 688          18    3  138
#> 689          38    4  138
#> 690           4    5  138
#> 691           0    1  139
#> 692           0    2  139
#> 693          21    3  139
#> 694          18    4  139
#> 695         196    5  139
#> 696           0    1  140
#> 697           0    2  140
#> 698          16    3  140
#> 699          58    4  140
#> 700          66    5  140
#> 701           0    1  141
#> 702           0    2  141
#> 703          20    3  141
#> 704          35    4  141
#> 705          47    5  141
#> 706           0    1  142
#> 707           0    2  142
#> 708          23    3  142
#> 709          56    4  142
#> 710           4    5  142
#> 711           0    1  143
#> 712           0    2  143
#> 713          22    3  143
#> 714          23    4  143
#> 715          37    5  143
#> 716           0    1  144
#> 717           0    2  144
#> 718           9    3  144
#> 719          49    4  144
#> 720          62    5  144
#> 721           0    1  145
#> 722           0    2  145
#> 723           4    3  145
#> 724          33    4  145
#> 725         180    5  145
#> 726           0    1  146
#> 727           0    2  146
#> 728          21    3  146
#> 729          39    4  146
#> 730          33    5  146
#> 731           0    1  147
#> 732           0    2  147
#> 733           9    3  147
#> 734          53    4  147
#> 735          10    5  147
#> 736           0    1  148
#> 737           0    2  148
#> 738          11    3  148
#> 739          20    4  148
#> 740          28    5  148
#> 741           0    1  149
#> 742           0    2  149
#> 743          22    3  149
#> 744          43    4  149
#> 745           2    5  149
#> 746           0    1  150
#> 747           0    2  150
#> 748          23    3  150
#> 749          59    4  150
#> 750          45    5  150
#> 751           0    1  151
#> 752           0    2  151
#> 753           6    3  151
#> 754          45    4  151
#> 755          26    5  151
#> 756           0    1  152
#> 757           0    2  152
#> 758          19    3  152
#> 759          28    4  152
#> 760          81    5  152
#> 761           0    1  153
#> 762           0    2  153
#> 763          14    3  153
#> 764          51    4  153
#> 765          91    5  153
#> 766           0    1  154
#> 767           0    2  154
#> 768          20    3  154
#> 769          32    4  154
#> 770          86    5  154
#> 771           0    1  155
#> 772           0    2  155
#> 773          18    3  155
#> 774          39    4  155
#> 775          13    5  155
#> 776           0    1  156
#> 777           0    2  156
#> 778          15    3  156
#> 779          21    4  156
#> 780           4    5  156
#> 781           0    1  157
#> 782           0    2  157
#> 783          16    3  157
#> 784          31    4  157
#> 785          52    5  157
#> 786           0    1  158
#> 787           0    2  158
#> 788          17    3  158
#> 789          54    4  158
#> 790         333    5  158
#> 791           0    1  159
#> 792           0    2  159
#> 793          12    3  159
#> 794          38    4  159
#> 795          16    5  159
#> 796           0    1  160
#> 797           0    2  160
#> 798           7    3  160
#> 799          55    4  160
#> 800          14    5  160
#> 801           0    1  161
#> 802           0    2  161
#> 803           7    3  161
#> 804          38    4  161
#> 805           5    5  161
#> 806           0    1  162
#> 807           0    2  162
#> 808          20    3  162
#> 809          30    4  162
#> 810          44    5  162
#> 811           0    1  163
#> 812           0    2  163
#> 813          18    3  163
#> 814          47    4  163
#> 815          18    5  163
#> 816           0    1  164
#> 817           0    2  164
#> 818          16    3  164
#> 819          44    4  164
#> 820         117    5  164
#> 821           0    1  165
#> 822           0    2  165
#> 823          21    3  165
#> 824          37    4  165
#> 825         106    5  165
#> 826           0    1  166
#> 827           0    2  166
#> 828          11    3  166
#> 829          52    4  166
#> 830         141    5  166
#> 831           0    1  167
#> 832           0    2  167
#> 833          17    3  167
#> 834          45    4  167
#> 835         320    5  167
#> 836           0    1  168
#> 837           0    2  168
#> 838          19    3  168
#> 839          52    4  168
#> 840         204    5  168
#> 841           0    1  169
#> 842           0    2  169
#> 843          15    3  169
#> 844          48    4  169
#> 845          21    5  169
#> 846           0    1  170
#> 847           0    2  170
#> 848          21    3  170
#> 849          34    4  170
#> 850          64    5  170
#> 851           0    1  171
#> 852           0    2  171
#> 853          17    3  171
#> 854          30    4  171
#> 855          86    5  171
#> 856           0    1  172
#> 857           0    2  172
#> 858          21    3  172
#> 859          36    4  172
#> 860         116    5  172
#> 861           0    1  173
#> 862           0    2  173
#> 863           6    3  173
#> 864          18    4  173
#> 865          40    5  173
#> 866           0    1  174
#> 867           0    2  174
#> 868          11    3  174
#> 869          39    4  174
#> 870           0    5  174
#> 871           0    1  175
#> 872           0    2  175
#> 873          26    3  175
#> 874          44    4  175
#> 875         164    5  175
#> 876           0    1  176
#> 877           0    2  176
#> 878          12    3  176
#> 879          47    4  176
#> 880          22    5  176
#> 881           0    1  177
#> 882           0    2  177
#> 883          18    3  177
#> 884          48    4  177
#> 885          56    5  177
#> 886           0    1  178
#> 887           0    2  178
#> 888          11    3  178
#> 889          35    4  178
#> 890          68    5  178
#> 891           0    1  179
#> 892           0    2  179
#> 893          12    3  179
#> 894          30    4  179
#> 895          48    5  179
#> 896           0    1  180
#> 897           0    2  180
#> 898          19    3  180
#> 899          22    4  180
#> 900          10    5  180
#> 901           0    1  181
#> 902           0    2  181
#> 903          24    3  181
#> 904          49    4  181
#> 905          90    5  181
#> 906           0    1  182
#> 907           0    2  182
#> 908          19    3  182
#> 909          38    4  182
#> 910         332    5  182
#> 911           0    1  183
#> 912           0    2  183
#> 913          19    3  183
#> 914          23    4  183
#> 915           3    5  183
#> 916           0    1  184
#> 917           0    2  184
#> 918          15    3  184
#> 919          20    4  184
#> 920          30    5  184
#> 921           0    1  185
#> 922           0    2  185
#> 923          10    3  185
#> 924          47    4  185
#> 925          27    5  185
#> 926           0    1  186
#> 927           0    2  186
#> 928          17    3  186
#> 929          54    4  186
#> 930         121    5  186
#> 931           0    1  187
#> 932           0    2  187
#> 933          21    3  187
#> 934          51    4  187
#> 935          65    5  187
#> 936           0    1  188
#> 937           0    2  188
#> 938          29    3  188
#> 939          32    4  188
#> 940          45    5  188
#> 941           0    1  189
#> 942           0    2  189
#> 943          21    3  189
#> 944          54    4  189
#> 945         132    5  189
#> 946           0    1  190
#> 947           0    2  190
#> 948          15    3  190
#> 949          53    4  190
#> 950          38    5  190
#> 951           0    1  191
#> 952           0    2  191
#> 953          19    3  191
#> 954          37    4  191
#> 955          28    5  191
#> 956           0    1  192
#> 957           0    2  192
#> 958          21    3  192
#> 959          43    4  192
#> 960           3    5  192
#> 961           0    1  193
#> 962           0    2  193
#> 963          20    3  193
#> 964          40    4  193
#> 965           0    5  193
#> 966           0    1  194
#> 967           0    2  194
#> 968          12    3  194
#> 969          32    4  194
#> 970           3    5  194
#> 971           0    1  195
#> 972           0    2  195
#> 973          15    3  195
#> 974          43    4  195
#> 975         280    5  195
#> 976           0    1  196
#> 977           0    2  196
#> 978           9    3  196
#> 979          35    4  196
#> 980         218    5  196
#> 981           0    1  197
#> 982           0    2  197
#> 983           6    3  197
#> 984          42    4  197
#> 985          81    5  197
#> 986           0    1  198
#> 987           0    2  198
#> 988          12    3  198
#> 989          33    4  198
#> 990         141    5  198
#> 991           0    1  199
#> 992           0    2  199
#> 993          12    3  199
#> 994          41    4  199
#> 995         102    5  199
#> 996           0    1  200
#> 997           0    2  200
#> 998          18    3  200
#> 999          59    4  200
#> 1000         51    5  200
#> 1001          0    1  201
#> 1002          0    2  201
#> 1003         33    3  201
#> 1004         32    4  201
#> 1005          1    5  201
#> 1006          0    1  202
#> 1007          0    2  202
#> 1008          9    3  202
#> 1009         34    4  202
#> 1010         53    5  202
#> 1011          0    1  203
#> 1012          0    2  203
#> 1013         13    3  203
#> 1014         45    4  203
#> 1015        161    5  203
#> 1016          0    1  204
#> 1017          0    2  204
#> 1018         18    3  204
#> 1019         48    4  204
#> 1020        101    5  204
#> 1021          0    1  205
#> 1022          0    2  205
#> 1023         14    3  205
#> 1024         67    4  205
#> 1025          1    5  205
#> 1026          0    1  206
#> 1027          0    2  206
#> 1028         36    3  206
#> 1029         48    4  206
#> 1030          3    5  206
#> 1031          0    1  207
#> 1032          0    2  207
#> 1033         17    3  207
#> 1034         41    4  207
#> 1035         81    5  207
#> 1036          0    1  208
#> 1037          0    2  208
#> 1038         12    3  208
#> 1039         54    4  208
#> 1040         40    5  208
#> 1041          0    1  209
#> 1042          0    2  209
#> 1043         16    3  209
#> 1044         37    4  209
#> 1045         51    5  209
#> 1046          0    1  210
#> 1047          0    2  210
#> 1048         33    3  210
#> 1049         30    4  210
#> 1050        117    5  210
#> 1051          0    1  211
#> 1052          0    2  211
#> 1053          9    3  211
#> 1054         42    4  211
#> 1055         60    5  211
#> 1056          0    1  212
#> 1057          0    2  212
#> 1058         24    3  212
#> 1059         53    4  212
#> 1060         54    5  212
#> 1061          0    1  213
#> 1062          0    2  213
#> 1063         18    3  213
#> 1064         26    4  213
#> 1065        101    5  213
#> 1066          0    1  214
#> 1067          0    2  214
#> 1068         23    3  214
#> 1069         23    4  214
#> 1070         28    5  214
#> 1071          0    1  215
#> 1072          0    2  215
#> 1073         18    3  215
#> 1074         55    4  215
#> 1075         15    5  215
#> 1076          0    1  216
#> 1077          0    2  216
#> 1078         10    3  216
#> 1079         29    4  216
#> 1080         67    5  216
#> 1081          0    1  217
#> 1082          0    2  217
#> 1083          7    3  217
#> 1084         43    4  217
#> 1085         74    5  217
#> 1086          0    1  218
#> 1087          0    2  218
#> 1088         12    3  218
#> 1089         40    4  218
#> 1090          9    5  218
#> 1091          0    1  219
#> 1092          0    2  219
#> 1093         14    3  219
#> 1094         64    4  219
#> 1095         56    5  219
#> 1096          0    1  220
#> 1097          0    2  220
#> 1098         21    3  220
#> 1099         36    4  220
#> 1100          8    5  220
#> 1101          0    1  221
#> 1102          0    2  221
#> 1103         10    3  221
#> 1104         40    4  221
#> 1105         70    5  221
#> 1106          0    1  222
#> 1107          0    2  222
#> 1108         15    3  222
#> 1109         56    4  222
#> 1110         27    5  222
#> 1111          0    1  223
#> 1112          0    2  223
#> 1113         16    3  223
#> 1114         39    4  223
#> 1115         11    5  223
#> 1116          0    1  224
#> 1117          0    2  224
#> 1118          8    3  224
#> 1119         57    4  224
#> 1120          6    5  224
#> 1121          0    1  225
#> 1122          0    2  225
#> 1123         12    3  225
#> 1124         53    4  225
#> 1125         28    5  225
#> 1126          0    1  226
#> 1127          0    2  226
#> 1128         21    3  226
#> 1129         47    4  226
#> 1130         33    5  226
#> 1131          0    1  227
#> 1132          0    2  227
#> 1133         12    3  227
#> 1134         30    4  227
#> 1135          2    5  227
#> 1136          0    1  228
#> 1137          0    2  228
#> 1138         24    3  228
#> 1139         34    4  228
#> 1140          4    5  228
#> 1141          0    1  229
#> 1142          0    2  229
#> 1143         15    3  229
#> 1144         21    4  229
#> 1145         41    5  229
#> 1146          0    1  230
#> 1147          0    2  230
#> 1148          8    3  230
#> 1149         31    4  230
#> 1150         75    5  230
#> 1151          0    1  231
#> 1152          0    2  231
#> 1153         23    3  231
#> 1154         74    4  231
#> 1155         54    5  231
#> 1156          0    1  232
#> 1157          0    2  232
#> 1158          5    3  232
#> 1159         40    4  232
#> 1160          7    5  232
#> 1161          0    1  233
#> 1162          0    2  233
#> 1163         15    3  233
#> 1164         57    4  233
#> 1165        203    5  233
#> 1166          0    1  234
#> 1167          0    2  234
#> 1168         13    3  234
#> 1169         53    4  234
#> 1170         46    5  234
#> 1171          0    1  235
#> 1172          0    2  235
#> 1173         12    3  235
#> 1174         15    4  235
#> 1175         15    5  235
#> 1176          0    1  236
#> 1177          0    2  236
#> 1178         24    3  236
#> 1179         26    4  236
#> 1180         68    5  236
#> 1181          0    1  237
#> 1182          0    2  237
#> 1183         19    3  237
#> 1184         41    4  237
#> 1185          8    5  237
#> 1186          0    1  238
#> 1187          0    2  238
#> 1188         25    3  238
#> 1189         20    4  238
#> 1190         27    5  238
#> 1191          0    1  239
#> 1192          0    2  239
#> 1193         17    3  239
#> 1194         63    4  239
#> 1195         34    5  239
#> 1196          0    1  240
#> 1197          0    2  240
#> 1198         14    3  240
#> 1199         27    4  240
#> 1200          0    5  240
#> 1201          0    1  241
#> 1202          0    2  241
#> 1203         20    3  241
#> 1204         52    4  241
#> 1205          9    5  241
#> 1206          0    1  242
#> 1207          0    2  242
#> 1208         22    3  242
#> 1209         48    4  242
#> 1210         30    5  242
#> 1211          0    1  243
#> 1212          0    2  243
#> 1213          5    3  243
#> 1214         45    4  243
#> 1215          8    5  243
#> 1216          0    1  244
#> 1217          0    2  244
#> 1218         15    3  244
#> 1219         38    4  244
#> 1220        113    5  244
#> 1221          0    1  245
#> 1222          0    2  245
#> 1223         19    3  245
#> 1224         60    4  245
#> 1225         41    5  245
#> 1226          0    1  246
#> 1227          0    2  246
#> 1228         13    3  246
#> 1229         42    4  246
#> 1230        152    5  246
#> 1231          0    1  247
#> 1232          0    2  247
#> 1233         21    3  247
#> 1234         36    4  247
#> 1235         48    5  247
#> 1236          0    1  248
#> 1237          0    2  248
#> 1238         21    3  248
#> 1239         33    4  248
#> 1240        108    5  248
#> 1241          0    1  249
#> 1242          0    2  249
#> 1243         18    3  249
#> 1244         66    4  249
#> 1245         42    5  249
#> 1246          0    1  250
#> 1247          0    2  250
#> 1248         14    3  250
#> 1249         27    4  250
#> 1250          8    5  250
#> 1251          0    1  251
#> 1252          0    2  251
#> 1253         21    3  251
#> 1254         32    4  251
#> 1255        154    5  251
#> 1256          0    1  252
#> 1257          0    2  252
#> 1258         11    3  252
#> 1259         53    4  252
#> 1260          8    5  252
#> 1261          0    1  253
#> 1262          0    2  253
#> 1263         10    3  253
#> 1264         31    4  253
#> 1265          2    5  253
#> 1266          0    1  254
#> 1267          0    2  254
#> 1268         17    3  254
#> 1269         53    4  254
#> 1270         93    5  254
#> 1271          0    1  255
#> 1272          0    2  255
#> 1273         13    3  255
#> 1274         24    4  255
#> 1275         55    5  255
#> 1276          0    1  256
#> 1277          0    2  256
#> 1278         11    3  256
#> 1279         21    4  256
#> 1280        159    5  256
#> 1281          0    1  257
#> 1282          0    2  257
#> 1283         19    3  257
#> 1284         35    4  257
#> 1285        226    5  257
#> 1286          0    1  258
#> 1287          0    2  258
#> 1288         28    3  258
#> 1289         36    4  258
#> 1290          1    5  258
#> 1291          0    1  259
#> 1292          0    2  259
#> 1293         12    3  259
#> 1294         37    4  259
#> 1295          0    5  259
#> 1296          0    1  260
#> 1297          0    2  260
#> 1298         22    3  260
#> 1299         53    4  260
#> 1300         34    5  260
#> 1301          0    1  261
#> 1302          0    2  261
#> 1303         20    3  261
#> 1304         48    4  261
#> 1305         19    5  261
#> 1306          0    1  262
#> 1307          0    2  262
#> 1308         16    3  262
#> 1309         46    4  262
#> 1310         29    5  262
#> 1311          0    1  263
#> 1312          0    2  263
#> 1313         11    3  263
#> 1314         36    4  263
#> 1315         97    5  263
#> 1316          0    1  264
#> 1317          0    2  264
#> 1318         20    3  264
#> 1319         25    4  264
#> 1320         36    5  264
#> 1321          0    1  265
#> 1322          0    2  265
#> 1323         26    3  265
#> 1324         61    4  265
#> 1325         36    5  265
#> 1326          0    1  266
#> 1327          0    2  266
#> 1328         28    3  266
#> 1329         42    4  266
#> 1330          5    5  266
#> 1331          0    1  267
#> 1332          0    2  267
#> 1333         25    3  267
#> 1334         53    4  267
#> 1335         54    5  267
#> 1336          0    1  268
#> 1337          0    2  268
#> 1338         12    3  268
#> 1339         40    4  268
#> 1340         91    5  268
#> 1341          0    1  269
#> 1342          0    2  269
#> 1343          9    3  269
#> 1344         32    4  269
#> 1345          5    5  269
#> 1346          0    1  270
#> 1347          0    2  270
#> 1348          7    3  270
#> 1349         52    4  270
#> 1350         55    5  270
#> 1351          0    1  271
#> 1352          0    2  271
#> 1353         20    3  271
#> 1354         27    4  271
#> 1355         96    5  271
#> 1356          0    1  272
#> 1357          0    2  272
#> 1358         25    3  272
#> 1359         50    4  272
#> 1360         15    5  272
#> 1361          0    1  273
#> 1362          0    2  273
#> 1363         21    3  273
#> 1364         38    4  273
#> 1365         70    5  273
#> 1366          0    1  274
#> 1367          0    2  274
#> 1368         10    3  274
#> 1369         41    4  274
#> 1370        175    5  274
#> 1371          0    1  275
#> 1372          0    2  275
#> 1373         14    3  275
#> 1374         53    4  275
#> 1375          2    5  275
#> 1376          0    1  276
#> 1377          0    2  276
#> 1378         21    3  276
#> 1379         23    4  276
#> 1380         22    5  276
#> 1381          0    1  277
#> 1382          0    2  277
#> 1383         22    3  277
#> 1384         37    4  277
#> 1385         35    5  277
#> 1386          0    1  278
#> 1387          0    2  278
#> 1388         12    3  278
#> 1389         21    4  278
#> 1390        197    5  278
#> 1391          0    1  279
#> 1392          0    2  279
#> 1393         17    3  279
#> 1394         21    4  279
#> 1395         28    5  279
#> 1396          0    1  280
#> 1397          0    2  280
#> 1398         24    3  280
#> 1399         43    4  280
#> 1400         29    5  280
#> 1401          0    1  281
#> 1402          0    2  281
#> 1403          4    3  281
#> 1404         52    4  281
#> 1405         18    5  281
#> 1406          0    1  282
#> 1407          0    2  282
#> 1408         20    3  282
#> 1409         19    4  282
#> 1410         43    5  282
#> 1411          0    1  283
#> 1412          0    2  283
#> 1413         23    3  283
#> 1414         56    4  283
#> 1415        127    5  283
#> 1416          0    1  284
#> 1417          0    2  284
#> 1418         18    3  284
#> 1419         54    4  284
#> 1420          0    5  284
#> 1421          0    1  285
#> 1422          0    2  285
#> 1423         19    3  285
#> 1424         53    4  285
#> 1425         31    5  285
#> 1426          0    1  286
#> 1427          0    2  286
#> 1428         15    3  286
#> 1429         32    4  286
#> 1430        110    5  286
#> 1431          0    1  287
#> 1432          0    2  287
#> 1433         18    3  287
#> 1434         43    4  287
#> 1435        100    5  287
#> 1436          0    1  288
#> 1437          0    2  288
#> 1438         16    3  288
#> 1439         39    4  288
#> 1440         86    5  288
#> 1441          0    1  289
#> 1442          0    2  289
#> 1443         29    3  289
#> 1444         21    4  289
#> 1445         26    5  289
#> 1446          0    1  290
#> 1447          0    2  290
#> 1448         30    3  290
#> 1449         46    4  290
#> 1450         13    5  290
#> 1451          0    1  291
#> 1452          0    2  291
#> 1453         27    3  291
#> 1454         42    4  291
#> 1455          0    5  291
#> 1456          0    1  292
#> 1457          0    2  292
#> 1458         16    3  292
#> 1459         49    4  292
#> 1460          4    5  292
#> 1461          0    1  293
#> 1462          0    2  293
#> 1463         26    3  293
#> 1464         44    4  293
#> 1465         17    5  293
#> 1466          0    1  294
#> 1467          0    2  294
#> 1468          9    3  294
#> 1469         28    4  294
#> 1470          7    5  294
#> 1471          0    1  295
#> 1472          0    2  295
#> 1473         12    3  295
#> 1474         52    4  295
#> 1475          7    5  295
#> 1476          0    1  296
#> 1477          0    2  296
#> 1478         14    3  296
#> 1479         63    4  296
#> 1480         41    5  296
#> 1481          0    1  297
#> 1482          0    2  297
#> 1483         11    3  297
#> 1484         36    4  297
#> 1485         66    5  297
#> 1486          0    1  298
#> 1487          0    2  298
#> 1488          6    3  298
#> 1489         20    4  298
#> 1490         33    5  298
#> 1491          0    1  299
#> 1492          0    2  299
#> 1493         13    3  299
#> 1494         76    4  299
#> 1495         93    5  299
#> 1496          0    1  300
#> 1497          0    2  300
#> 1498          6    3  300
#> 1499         36    4  300
#> 1500        113    5  300
#> 1501          0    1  301
#> 1502          0    2  301
#> 1503         25    3  301
#> 1504         44    4  301
#> 1505        167    5  301
#> 1506          0    1  302
#> 1507          0    2  302
#> 1508         19    3  302
#> 1509         37    4  302
#> 1510          2    5  302
#> 1511          0    1  303
#> 1512          0    2  303
#> 1513         11    3  303
#> 1514         27    4  303
#> 1515        200    5  303
#> 1516          0    1  304
#> 1517          0    2  304
#> 1518         21    3  304
#> 1519         45    4  304
#> 1520        103    5  304
#> 1521          0    1  305
#> 1522          0    2  305
#> 1523         14    3  305
#> 1524         21    4  305
#> 1525         29    5  305
#> 1526          0    1  306
#> 1527          0    2  306
#> 1528          8    3  306
#> 1529         21    4  306
#> 1530         52    5  306
#> 1531          0    1  307
#> 1532          0    2  307
#> 1533         14    3  307
#> 1534         39    4  307
#> 1535          8    5  307
#> 1536          0    1  308
#> 1537          0    2  308
#> 1538         16    3  308
#> 1539         44    4  308
#> 1540         14    5  308
#> 1541          0    1  309
#> 1542          0    2  309
#> 1543         19    3  309
#> 1544         30    4  309
#> 1545        173    5  309
#> 1546          0    1  310
#> 1547          0    2  310
#> 1548         19    3  310
#> 1549         26    4  310
#> 1550         27    5  310
#> 1551          0    1  311
#> 1552          0    2  311
#> 1553         22    3  311
#> 1554         21    4  311
#> 1555         45    5  311
#> 1556          0    1  312
#> 1557          0    2  312
#> 1558         13    3  312
#> 1559         19    4  312
#> 1560        126    5  312
#> 1561          0    1  313
#> 1562          0    2  313
#> 1563         16    3  313
#> 1564         39    4  313
#> 1565         77    5  313
#> 1566          0    1  314
#> 1567          0    2  314
#> 1568         23    3  314
#> 1569         27    4  314
#> 1570          8    5  314
#> 1571          0    1  315
#> 1572          0    2  315
#> 1573          9    3  315
#> 1574         28    4  315
#> 1575         14    5  315
#> 1576          0    1  316
#> 1577          0    2  316
#> 1578         12    3  316
#> 1579         50    4  316
#> 1580        117    5  316
#> 1581          0    1  317
#> 1582          0    2  317
#> 1583         14    3  317
#> 1584         44    4  317
#> 1585         84    5  317
#> 1586          0    1  318
#> 1587          0    2  318
#> 1588          7    3  318
#> 1589         52    4  318
#> 1590          1    5  318
#> 1591          0    1  319
#> 1592          0    2  319
#> 1593         12    3  319
#> 1594         41    4  319
#> 1595          3    5  319
#> 1596          0    1  320
#> 1597          0    2  320
#> 1598         17    3  320
#> 1599         43    4  320
#> 1600         38    5  320
#> 1601          0    1  321
#> 1602          0    2  321
#> 1603          9    3  321
#> 1604         45    4  321
#> 1605         30    5  321
#> 1606          0    1  322
#> 1607          0    2  322
#> 1608         17    3  322
#> 1609         29    4  322
#> 1610        111    5  322
#> 1611          0    1  323
#> 1612          0    2  323
#> 1613          8    3  323
#> 1614         63    4  323
#> 1615          2    5  323
#> 1616          0    1  324
#> 1617          0    2  324
#> 1618          9    3  324
#> 1619         51    4  324
#> 1620         88    5  324
#> 1621          0    1  325
#> 1622          0    2  325
#> 1623         17    3  325
#> 1624         53    4  325
#> 1625        128    5  325
#> 1626          0    1  326
#> 1627          0    2  326
#> 1628         14    3  326
#> 1629         56    4  326
#> 1630         37    5  326
#> 1631          0    1  327
#> 1632          0    2  327
#> 1633          9    3  327
#> 1634         65    4  327
#> 1635         31    5  327
#> 1636          0    1  328
#> 1637          0    2  328
#> 1638         10    3  328
#> 1639         27    4  328
#> 1640        166    5  328
#> 1641          0    1  329
#> 1642          0    2  329
#> 1643         13    3  329
#> 1644         33    4  329
#> 1645         16    5  329
#> 1646          0    1  330
#> 1647          0    2  330
#> 1648         14    3  330
#> 1649         47    4  330
#> 1650         57    5  330
#> 1651          0    1  331
#> 1652          0    2  331
#> 1653         11    3  331
#> 1654         53    4  331
#> 1655         71    5  331
#> 1656          0    1  332
#> 1657          0    2  332
#> 1658         20    3  332
#> 1659         46    4  332
#> 1660         55    5  332
#> 1661          0    1  333
#> 1662          0    2  333
#> 1663         15    3  333
#> 1664         46    4  333
#> 1665          2    5  333
#> 1666          0    1  334
#> 1667          0    2  334
#> 1668          9    3  334
#> 1669         30    4  334
#> 1670         72    5  334
#> 1671          0    1  335
#> 1672          0    2  335
#> 1673          4    3  335
#> 1674         70    4  335
#> 1675          4    5  335
#> 1676          0    1  336
#> 1677          0    2  336
#> 1678         31    3  336
#> 1679         40    4  336
#> 1680          3    5  336
#> 1681          0    1  337
#> 1682          0    2  337
#> 1683          6    3  337
#> 1684         49    4  337
#> 1685         93    5  337
#> 1686          0    1  338
#> 1687          0    2  338
#> 1688         10    3  338
#> 1689         56    4  338
#> 1690         70    5  338
#> 1691          0    1  339
#> 1692          0    2  339
#> 1693         15    3  339
#> 1694         63    4  339
#> 1695         49    5  339
#> 1696          0    1  340
#> 1697          0    2  340
#> 1698          5    3  340
#> 1699         28    4  340
#> 1700          5    5  340
#> 1701          0    1  341
#> 1702          0    2  341
#> 1703         12    3  341
#> 1704         31    4  341
#> 1705         33    5  341
#> 1706          0    1  342
#> 1707          0    2  342
#> 1708         37    3  342
#> 1709         27    4  342
#> 1710         82    5  342
#> 1711          0    1  343
#> 1712          0    2  343
#> 1713         18    3  343
#> 1714         23    4  343
#> 1715         79    5  343
#> 1716          0    1  344
#> 1717          0    2  344
#> 1718         11    3  344
#> 1719         78    4  344
#> 1720        231    5  344
#> 1721          0    1  345
#> 1722          0    2  345
#> 1723         17    3  345
#> 1724         39    4  345
#> 1725          0    5  345
#> 1726          0    1  346
#> 1727          0    2  346
#> 1728         13    3  346
#> 1729         41    4  346
#> 1730         81    5  346
#> 1731          0    1  347
#> 1732          0    2  347
#> 1733         25    3  347
#> 1734         22    4  347
#> 1735          3    5  347
#> 1736          0    1  348
#> 1737          0    2  348
#> 1738         16    3  348
#> 1739         27    4  348
#> 1740         52    5  348
#> 1741          0    1  349
#> 1742          0    2  349
#> 1743         32    3  349
#> 1744         27    4  349
#> 1745        127    5  349
#> 1746          0    1  350
#> 1747          0    2  350
#> 1748         11    3  350
#> 1749         77    4  350
#> 1750         36    5  350
#> 1751          0    1  351
#> 1752          0    2  351
#> 1753         20    3  351
#> 1754         25    4  351
#> 1755         17    5  351
#> 1756          0    1  352
#> 1757          0    2  352
#> 1758         16    3  352
#> 1759         38    4  352
#> 1760        312    5  352
#> 1761          0    1  353
#> 1762          0    2  353
#> 1763         20    3  353
#> 1764         69    4  353
#> 1765          2    5  353
#> 1766          0    1  354
#> 1767          0    2  354
#> 1768         14    3  354
#> 1769         49    4  354
#> 1770         71    5  354
#> 1771          0    1  355
#> 1772          0    2  355
#> 1773         25    3  355
#> 1774         27    4  355
#> 1775         46    5  355
#> 1776          0    1  356
#> 1777          0    2  356
#> 1778         21    3  356
#> 1779         34    4  356
#> 1780          4    5  356
#> 1781          0    1  357
#> 1782          0    2  357
#> 1783         21    3  357
#> 1784         29    4  357
#> 1785        195    5  357
#> 1786          0    1  358
#> 1787          0    2  358
#> 1788          6    3  358
#> 1789         70    4  358
#> 1790          3    5  358
#> 1791          0    1  359
#> 1792          0    2  359
#> 1793         25    3  359
#> 1794         44    4  359
#> 1795          0    5  359
#> 1796          0    1  360
#> 1797          0    2  360
#> 1798         23    3  360
#> 1799         63    4  360
#> 1800         74    5  360
#> 1801          0    1  361
#> 1802          0    2  361
#> 1803         19    3  361
#> 1804         44    4  361
#> 1805         60    5  361
#> 1806          0    1  362
#> 1807          0    2  362
#> 1808         35    3  362
#> 1809         74    4  362
#> 1810        313    5  362
#> 1811          0    1  363
#> 1812          0    2  363
#> 1813         33    3  363
#> 1814         36    4  363
#> 1815        108    5  363
#> 1816          0    1  364
#> 1817          0    2  364
#> 1818         24    3  364
#> 1819         23    4  364
#> 1820        274    5  364
#> 1821          0    1  365
#> 1822          0    2  365
#> 1823         18    3  365
#> 1824         37    4  365
#> 1825         66    5  365
#> 1826          0    1  366
#> 1827          0    2  366
#> 1828          8    3  366
#> 1829         55    4  366
#> 1830          3    5  366
#> 1831          0    1  367
#> 1832          0    2  367
#> 1833         15    3  367
#> 1834         36    4  367
#> 1835         18    5  367
#> 1836          0    1  368
#> 1837          0    2  368
#> 1838         23    3  368
#> 1839         39    4  368
#> 1840         32    5  368
#> 1841          0    1  369
#> 1842          0    2  369
#> 1843         12    3  369
#> 1844         35    4  369
#> 1845          4    5  369
#> 1846          0    1  370
#> 1847          0    2  370
#> 1848         30    3  370
#> 1849         39    4  370
#> 1850         26    5  370
#> 1851          0    1  371
#> 1852          0    2  371
#> 1853         16    3  371
#> 1854         25    4  371
#> 1855          7    5  371
#> 1856          0    1  372
#> 1857          0    2  372
#> 1858          9    3  372
#> 1859         23    4  372
#> 1860         69    5  372
#> 1861          0    1  373
#> 1862          0    2  373
#> 1863         15    3  373
#> 1864         28    4  373
#> 1865          3    5  373
#> 1866          0    1  374
#> 1867          0    2  374
#> 1868         16    3  374
#> 1869         46    4  374
#> 1870         95    5  374
#> 1871          0    1  375
#> 1872          0    2  375
#> 1873         16    3  375
#> 1874         45    4  375
#> 1875          6    5  375
#> 1876          0    1  376
#> 1877          0    2  376
#> 1878         16    3  376
#> 1879         37    4  376
#> 1880         28    5  376
#> 1881          0    1  377
#> 1882          0    2  377
#> 1883         17    3  377
#> 1884         67    4  377
#> 1885        129    5  377
#> 1886          0    1  378
#> 1887          0    2  378
#> 1888         12    3  378
#> 1889         56    4  378
#> 1890         69    5  378
#> 1891          0    1  379
#> 1892          0    2  379
#> 1893         13    3  379
#> 1894         27    4  379
#> 1895          8    5  379
#> 1896          0    1  380
#> 1897          0    2  380
#> 1898         21    3  380
#> 1899         23    4  380
#> 1900          2    5  380
#> 1901          0    1  381
#> 1902          0    2  381
#> 1903         25    3  381
#> 1904         22    4  381
#> 1905         21    5  381
#> 1906          0    1  382
#> 1907          0    2  382
#> 1908         27    3  382
#> 1909         42    4  382
#> 1910        166    5  382
#> 1911          0    1  383
#> 1912          0    2  383
#> 1913         14    3  383
#> 1914         46    4  383
#> 1915         70    5  383
#> 1916          0    1  384
#> 1917          0    2  384
#> 1918         15    3  384
#> 1919         33    4  384
#> 1920         26    5  384
#> 1921          0    1  385
#> 1922          0    2  385
#> 1923         14    3  385
#> 1924         36    4  385
#> 1925          9    5  385
#> 1926          0    1  386
#> 1927          0    2  386
#> 1928         33    3  386
#> 1929         37    4  386
#> 1930          1    5  386
#> 1931          0    1  387
#> 1932          0    2  387
#> 1933          8    3  387
#> 1934         52    4  387
#> 1935         61    5  387
#> 1936          0    1  388
#> 1937          0    2  388
#> 1938         23    3  388
#> 1939         20    4  388
#> 1940         32    5  388
#> 1941          0    1  389
#> 1942          0    2  389
#> 1943         13    3  389
#> 1944         50    4  389
#> 1945         24    5  389
#> 1946          0    1  390
#> 1947          0    2  390
#> 1948         13    3  390
#> 1949         39    4  390
#> 1950        176    5  390
#> 1951          0    1  391
#> 1952          0    2  391
#> 1953         15    3  391
#> 1954         51    4  391
#> 1955        234    5  391
#> 1956          0    1  392
#> 1957          0    2  392
#> 1958         17    3  392
#> 1959         47    4  392
#> 1960         10    5  392
#> 1961          0    1  393
#> 1962          0    2  393
#> 1963         12    3  393
#> 1964         41    4  393
#> 1965         42    5  393
#> 1966          0    1  394
#> 1967          0    2  394
#> 1968         33    3  394
#> 1969         40    4  394
#> 1970         22    5  394
#> 1971          0    1  395
#> 1972          0    2  395
#> 1973          7    3  395
#> 1974         35    4  395
#> 1975         54    5  395
#> 1976          0    1  396
#> 1977          0    2  396
#> 1978         23    3  396
#> 1979         43    4  396
#> 1980         60    5  396
#> 1981          0    1  397
#> 1982          0    2  397
#> 1983         30    3  397
#> 1984         39    4  397
#> 1985         32    5  397
#> 1986          0    1  398
#> 1987          0    2  398
#> 1988         17    3  398
#> 1989         49    4  398
#> 1990         11    5  398
#> 1991          0    1  399
#> 1992          0    2  399
#> 1993         10    3  399
#> 1994         37    4  399
#> 1995          4    5  399
#> 1996          0    1  400
#> 1997          0    2  400
#> 1998         22    3  400
#> 1999         39    4  400
#> 2000         24    5  400
#> 2001          0    1  401
#> 2002          0    2  401
#> 2003         20    3  401
#> 2004         29    4  401
#> 2005         11    5  401
#> 2006          0    1  402
#> 2007          0    2  402
#> 2008         32    3  402
#> 2009         37    4  402
#> 2010         26    5  402
#> 2011          0    1  403
#> 2012          0    2  403
#> 2013          7    3  403
#> 2014         36    4  403
#> 2015         81    5  403
#> 2016          0    1  404
#> 2017          0    2  404
#> 2018         16    3  404
#> 2019         40    4  404
#> 2020         79    5  404
#> 2021          0    1  405
#> 2022          0    2  405
#> 2023         10    3  405
#> 2024         31    4  405
#> 2025         29    5  405
#> 2026          0    1  406
#> 2027          0    2  406
#> 2028         23    3  406
#> 2029         41    4  406
#> 2030         61    5  406
#> 2031          0    1  407
#> 2032          0    2  407
#> 2033         10    3  407
#> 2034         32    4  407
#> 2035          6    5  407
#> 2036          0    1  408
#> 2037          0    2  408
#> 2038         21    3  408
#> 2039         42    4  408
#> 2040         32    5  408
#> 2041          0    1  409
#> 2042          0    2  409
#> 2043         31    3  409
#> 2044         35    4  409
#> 2045         55    5  409
#> 2046          0    1  410
#> 2047          0    2  410
#> 2048         25    3  410
#> 2049         35    4  410
#> 2050         21    5  410
#> 2051          0    1  411
#> 2052          0    2  411
#> 2053          9    3  411
#> 2054         62    4  411
#> 2055         64    5  411
#> 2056          0    1  412
#> 2057          0    2  412
#> 2058         13    3  412
#> 2059         56    4  412
#> 2060         38    5  412
#> 2061          0    1  413
#> 2062          0    2  413
#> 2063         18    3  413
#> 2064         52    4  413
#> 2065          9    5  413
#> 2066          0    1  414
#> 2067          0    2  414
#> 2068         11    3  414
#> 2069         50    4  414
#> 2070          1    5  414
#> 2071          0    1  415
#> 2072          0    2  415
#> 2073         24    3  415
#> 2074         43    4  415
#> 2075         75    5  415
#> 2076          0    1  416
#> 2077          0    2  416
#> 2078         32    3  416
#> 2079         21    4  416
#> 2080         72    5  416
#> 2081          0    1  417
#> 2082          0    2  417
#> 2083         17    3  417
#> 2084         41    4  417
#> 2085         77    5  417
#> 2086          0    1  418
#> 2087          0    2  418
#> 2088         10    3  418
#> 2089         41    4  418
#> 2090          0    5  418
#> 2091          0    1  419
#> 2092          0    2  419
#> 2093         11    3  419
#> 2094         34    4  419
#> 2095          5    5  419
#> 2096          0    1  420
#> 2097          0    2  420
#> 2098         11    3  420
#> 2099         33    4  420
#> 2100        101    5  420
#> 2101          0    1  421
#> 2102          0    2  421
#> 2103         15    3  421
#> 2104         23    4  421
#> 2105          1    5  421
#> 2106          0    1  422
#> 2107          0    2  422
#> 2108         11    3  422
#> 2109         46    4  422
#> 2110        192    5  422
#> 2111          0    1  423
#> 2112          0    2  423
#> 2113         21    3  423
#> 2114         41    4  423
#> 2115         68    5  423
#> 2116          0    1  424
#> 2117          0    2  424
#> 2118         42    3  424
#> 2119         24    4  424
#> 2120         34    5  424
#> 2121          0    1  425
#> 2122          0    2  425
#> 2123         10    3  425
#> 2124         29    4  425
#> 2125        154    5  425
#> 2126          0    1  426
#> 2127          0    2  426
#> 2128         24    3  426
#> 2129         49    4  426
#> 2130         25    5  426
#> 2131          0    1  427
#> 2132          0    2  427
#> 2133         16    3  427
#> 2134         26    4  427
#> 2135         86    5  427
#> 2136          0    1  428
#> 2137          0    2  428
#> 2138         18    3  428
#> 2139         43    4  428
#> 2140         26    5  428
#> 2141          0    1  429
#> 2142          0    2  429
#> 2143          7    3  429
#> 2144         46    4  429
#> 2145         64    5  429
#> 2146          0    1  430
#> 2147          0    2  430
#> 2148         19    3  430
#> 2149         53    4  430
#> 2150          1    5  430
#> 2151          0    1  431
#> 2152          0    2  431
#> 2153         19    3  431
#> 2154         28    4  431
#> 2155        118    5  431
#> 2156          0    1  432
#> 2157          0    2  432
#> 2158          6    3  432
#> 2159         48    4  432
#> 2160         56    5  432
#> 2161          0    1  433
#> 2162          0    2  433
#> 2163         13    3  433
#> 2164         46    4  433
#> 2165         15    5  433
#> 2166          0    1  434
#> 2167          0    2  434
#> 2168         21    3  434
#> 2169         30    4  434
#> 2170        304    5  434
#> 2171          0    1  435
#> 2172          0    2  435
#> 2173         24    3  435
#> 2174         32    4  435
#> 2175        108    5  435
#> 2176          0    1  436
#> 2177          0    2  436
#> 2178         13    3  436
#> 2179         44    4  436
#> 2180         25    5  436
#> 2181          0    1  437
#> 2182          0    2  437
#> 2183         11    3  437
#> 2184         53    4  437
#> 2185        145    5  437
#> 2186          0    1  438
#> 2187          0    2  438
#> 2188         33    3  438
#> 2189         41    4  438
#> 2190         33    5  438
#> 2191          0    1  439
#> 2192          0    2  439
#> 2193         27    3  439
#> 2194         47    4  439
#> 2195        100    5  439
#> 2196          0    1  440
#> 2197          0    2  440
#> 2198         12    3  440
#> 2199         38    4  440
#> 2200         25    5  440
#> 2201          0    1  441
#> 2202          0    2  441
#> 2203         18    3  441
#> 2204         39    4  441
#> 2205         98    5  441
#> 2206          0    1  442
#> 2207          0    2  442
#> 2208         22    3  442
#> 2209         27    4  442
#> 2210        236    5  442
#> 2211          0    1  443
#> 2212          0    2  443
#> 2213         12    3  443
#> 2214         20    4  443
#> 2215         34    5  443
#> 2216          0    1  444
#> 2217          0    2  444
#> 2218         10    3  444
#> 2219         20    4  444
#> 2220          8    5  444
#> 2221          0    1  445
#> 2222          0    2  445
#> 2223         22    3  445
#> 2224         24    4  445
#> 2225         80    5  445
#> 2226          0    1  446
#> 2227          0    2  446
#> 2228         15    3  446
#> 2229         27    4  446
#> 2230         25    5  446
#> 2231          0    1  447
#> 2232          0    2  447
#> 2233         31    3  447
#> 2234         29    4  447
#> 2235          8    5  447
#> 2236          0    1  448
#> 2237          0    2  448
#> 2238         14    3  448
#> 2239         66    4  448
#> 2240         78    5  448
#> 2241          0    1  449
#> 2242          0    2  449
#> 2243         22    3  449
#> 2244         48    4  449
#> 2245         60    5  449
#> 2246          0    1  450
#> 2247          0    2  450
#> 2248          9    3  450
#> 2249         39    4  450
#> 2250         61    5  450
#> 2251          0    1  451
#> 2252          0    2  451
#> 2253         11    3  451
#> 2254         39    4  451
#> 2255        107    5  451
#> 2256          0    1  452
#> 2257          0    2  452
#> 2258         16    3  452
#> 2259         55    4  452
#> 2260         96    5  452
#> 2261          0    1  453
#> 2262          0    2  453
#> 2263         21    3  453
#> 2264         45    4  453
#> 2265        286    5  453
#> 2266          0    1  454
#> 2267          0    2  454
#> 2268          6    3  454
#> 2269         37    4  454
#> 2270         30    5  454
#> 2271          0    1  455
#> 2272          0    2  455
#> 2273         18    3  455
#> 2274         50    4  455
#> 2275         14    5  455
#> 2276          0    1  456
#> 2277          0    2  456
#> 2278         20    3  456
#> 2279         80    4  456
#> 2280          2    5  456
#> 2281          0    1  457
#> 2282          0    2  457
#> 2283         19    3  457
#> 2284         63    4  457
#> 2285         48    5  457
#> 2286          0    1  458
#> 2287          0    2  458
#> 2288          7    3  458
#> 2289         54    4  458
#> 2290          2    5  458
#> 2291          0    1  459
#> 2292          0    2  459
#> 2293         16    3  459
#> 2294         43    4  459
#> 2295         14    5  459
#> 2296          0    1  460
#> 2297          0    2  460
#> 2298         15    3  460
#> 2299         39    4  460
#> 2300          1    5  460
#> 2301          0    1  461
#> 2302          0    2  461
#> 2303         24    3  461
#> 2304         42    4  461
#> 2305        391    5  461
#> 2306          0    1  462
#> 2307          0    2  462
#> 2308         14    3  462
#> 2309         60    4  462
#> 2310         87    5  462
#> 2311          0    1  463
#> 2312          0    2  463
#> 2313         22    3  463
#> 2314         32    4  463
#> 2315          4    5  463
#> 2316          0    1  464
#> 2317          0    2  464
#> 2318         11    3  464
#> 2319         38    4  464
#> 2320         25    5  464
#> 2321          0    1  465
#> 2322          0    2  465
#> 2323         11    3  465
#> 2324         38    4  465
#> 2325          0    5  465
#> 2326          0    1  466
#> 2327          0    2  466
#> 2328         12    3  466
#> 2329         35    4  466
#> 2330         28    5  466
#> 2331          0    1  467
#> 2332          0    2  467
#> 2333         23    3  467
#> 2334         68    4  467
#> 2335         17    5  467
#> 2336          0    1  468
#> 2337          0    2  468
#> 2338         16    3  468
#> 2339         37    4  468
#> 2340         48    5  468
#> 2341          0    1  469
#> 2342          0    2  469
#> 2343         20    3  469
#> 2344         34    4  469
#> 2345          2    5  469
#> 2346          0    1  470
#> 2347          0    2  470
#> 2348         18    3  470
#> 2349         51    4  470
#> 2350          7    5  470
#> 2351          0    1  471
#> 2352          0    2  471
#> 2353         29    3  471
#> 2354         43    4  471
#> 2355         10    5  471
#> 2356          0    1  472
#> 2357          0    2  472
#> 2358         19    3  472
#> 2359         60    4  472
#> 2360          0    5  472
#> 2361          0    1  473
#> 2362          0    2  473
#> 2363         15    3  473
#> 2364         20    4  473
#> 2365         70    5  473
#> 2366          0    1  474
#> 2367          0    2  474
#> 2368         10    3  474
#> 2369         69    4  474
#> 2370          2    5  474
#> 2371          0    1  475
#> 2372          0    2  475
#> 2373         22    3  475
#> 2374         55    4  475
#> 2375         27    5  475
#> 2376          0    1  476
#> 2377          0    2  476
#> 2378         20    3  476
#> 2379         41    4  476
#> 2380        134    5  476
#> 2381          0    1  477
#> 2382          0    2  477
#> 2383          8    3  477
#> 2384         43    4  477
#> 2385          4    5  477
#> 2386          0    1  478
#> 2387          0    2  478
#> 2388         10    3  478
#> 2389          7    4  478
#> 2390          6    5  478
#> 2391          0    1  479
#> 2392          0    2  479
#> 2393         29    3  479
#> 2394         56    4  479
#> 2395          4    5  479
#> 2396          0    1  480
#> 2397          0    2  480
#> 2398         15    3  480
#> 2399         44    4  480
#> 2400         70    5  480
#> 2401          0    1  481
#> 2402          0    2  481
#> 2403         14    3  481
#> 2404         71    4  481
#> 2405         88    5  481
#> 2406          0    1  482
#> 2407          0    2  482
#> 2408         11    3  482
#> 2409         40    4  482
#> 2410         12    5  482
#> 2411          0    1  483
#> 2412          0    2  483
#> 2413         12    3  483
#> 2414         48    4  483
#> 2415        220    5  483
#> 2416          0    1  484
#> 2417          0    2  484
#> 2418         12    3  484
#> 2419         21    4  484
#> 2420         19    5  484
#> 2421          0    1  485
#> 2422          0    2  485
#> 2423          8    3  485
#> 2424         32    4  485
#> 2425          1    5  485
#> 2426          0    1  486
#> 2427          0    2  486
#> 2428         23    3  486
#> 2429         50    4  486
#> 2430         49    5  486
#> 2431          0    1  487
#> 2432          0    2  487
#> 2433         17    3  487
#> 2434         47    4  487
#> 2435        117    5  487
#> 2436          0    1  488
#> 2437          0    2  488
#> 2438         13    3  488
#> 2439         37    4  488
#> 2440        192    5  488
#> 2441          0    1  489
#> 2442          0    2  489
#> 2443         25    3  489
#> 2444         51    4  489
#> 2445         88    5  489
#> 2446          0    1  490
#> 2447          0    2  490
#> 2448         16    3  490
#> 2449         31    4  490
#> 2450        232    5  490
#> 2451          0    1  491
#> 2452          0    2  491
#> 2453         30    3  491
#> 2454         74    4  491
#> 2455          8    5  491
#> 2456          0    1  492
#> 2457          0    2  492
#> 2458         13    3  492
#> 2459         58    4  492
#> 2460         83    5  492
#> 2461          0    1  493
#> 2462          0    2  493
#> 2463         28    3  493
#> 2464         57    4  493
#> 2465         15    5  493
#> 2466          0    1  494
#> 2467          0    2  494
#> 2468          4    3  494
#> 2469         51    4  494
#> 2470        159    5  494
#> 2471          0    1  495
#> 2472          0    2  495
#> 2473          9    3  495
#> 2474         51    4  495
#> 2475         21    5  495
#> 2476          0    1  496
#> 2477          0    2  496
#> 2478         13    3  496
#> 2479         20    4  496
#> 2480        152    5  496
#> 2481          0    1  497
#> 2482          0    2  497
#> 2483         12    3  497
#> 2484         51    4  497
#> 2485         38    5  497
#> 2486          0    1  498
#> 2487          0    2  498
#> 2488         16    3  498
#> 2489         49    4  498
#> 2490        135    5  498
#> 2491          0    1  499
#> 2492          0    2  499
#> 2493         22    3  499
#> 2494         54    4  499
#> 2495          6    5  499
#> 2496          0    1  500
#> 2497          0    2  500
#> 2498         31    3  500
#> 2499         41    4  500
#> 2500         71    5  500