├── memo
├── timer.h
├── asura.jpeg
├── pcaps
    └── 1.pcap
├── dtreeviz
    ├── 1.png
    ├── README.md
    ├── 1.py
    └── tmp
├── reduce_local_counts_to_global_count.h
├── reduce_local_sums_to_global_sum.h
├── view.h
├── reduce_min_ind.h
├── scatter.py
├── kmeans_tbb.cpp
├── Makefile
├── reduced2-bak
├── sum_and_count.h
├── repair_empty_clusters.cpp
├── kmeans.h
├── kmeans_tbb.h
├── tmp
├── README.md
├── test_kmeans.cpp
├── asura23.cpp
└── reduced2


/memo:
--------------------------------------------------------------------------------
1 | shuf -n 1000 reduced2 > tmp
2 | 
3 | 


--------------------------------------------------------------------------------
/timer.h:
--------------------------------------------------------------------------------
https://raw.githubusercontent.com/RuoAndo/Asura/HEAD/timer.h


--------------------------------------------------------------------------------
/asura.jpeg:
--------------------------------------------------------------------------------
https://raw.githubusercontent.com/RuoAndo/Asura/HEAD/asura.jpeg


--------------------------------------------------------------------------------
/pcaps/1.pcap:
--------------------------------------------------------------------------------
https://raw.githubusercontent.com/RuoAndo/Asura/HEAD/pcaps/1.pcap


--------------------------------------------------------------------------------
/dtreeviz/1.png:
--------------------------------------------------------------------------------
https://raw.githubusercontent.com/RuoAndo/Asura/HEAD/dtreeviz/1.png


--------------------------------------------------------------------------------
/dtreeviz/README.md:
--------------------------------------------------------------------------------
1 | <pre>
2 | # pip install pip install dtreeviz
3 | # python 1.py
4 | </pre>
5 | 
6 | <img src="1.png">
7 | 
8 | 


--------------------------------------------------------------------------------
/reduce_local_counts_to_global_count.h:
--------------------------------------------------------------------------------
1 | void convert_local_counts_to_global_count( tls_type& tls, view& global ) {
2 |     global.change = 0;
3 |     for( auto i=tls.begin(); i!=tls.end(); ++i ) {
4 |         view& v = *i;
5 |         global.change += i->change;
6 |         v.change = 0;
7 |     }
8 | }
9 | 


--------------------------------------------------------------------------------
/reduce_local_sums_to_global_sum.h:
--------------------------------------------------------------------------------
 1 | void convert_local_sums_to_global_sum( size_t k, tls_type& tls, view& global ) {
 2 |     for( auto i=tls.begin(); i!=tls.end(); ++i ) {
 3 |         view& v = *i;
 4 |         for( size_t j=0; j<k; ++j ) {
 5 |             global.array[j] += v.array[j];
 6 |             v.array[j].clear();
 7 |         }
 8 |     }
 9 | }
10 | 


--------------------------------------------------------------------------------
/view.h:
--------------------------------------------------------------------------------
 1 | class view {
 2 |     view( const view& v );            // Deny copy construction
 3 |     void operator=( const view& v );  // Deny assignment
 4 | public:
 5 |     sum_and_count* array;
 6 |     size_t change;
 7 |     view( size_t k ) : array(new sum_and_count[k]), change(0) {}
 8 |     ~view() {delete[] array;}
 9 | };
10 | 
11 | typedef tbb::enumerable_thread_specific<view> tls_type;


--------------------------------------------------------------------------------
/reduce_min_ind.h:
--------------------------------------------------------------------------------
 1 | int calc_shortest_index( const point centroid[], size_t k, point value ) { 
 2 |     int min = -1;
 3 |     float mind = std::numeric_limits<float>::max();
 4 |     for( int j=0; j<k; ++j ) {
 5 |         float d = distance2(centroid[j],value);
 6 |         if( d<mind  ) {
 7 |             mind = d;
 8 |             min = j;
 9 |         }
10 |     }
11 |     return min;
12 | }
13 | 


--------------------------------------------------------------------------------
/scatter.py:
--------------------------------------------------------------------------------
 1 | import numpy as np
 2 | import matplotlib.pyplot as plt
 3 | 
 4 | import sys
 5 | 
 6 | 
 7 | args = sys.argv
 8 | 
 9 | data_set = np.loadtxt(
10 |     fname=args[1],
11 |     dtype="float",
12 |     delimiter=",",
13 | )
14 | 
15 | for data in data_set:
16 |     plt.scatter(data[0], data[1])
17 | 
18 | plt.title("packets- src=>dst")
19 | plt.xlabel("counts")
20 | plt.ylabel("bytes")
21 | plt.grid()
22 | 
23 | plt.show()
24 | 


--------------------------------------------------------------------------------
/dtreeviz/1.py:
--------------------------------------------------------------------------------
 1 | import pandas as pd
 2 | import matplotlib.pyplot as plt
 3 | from sklearn.tree import DecisionTreeRegressor
 4 | from dtreeviz.trees import *
 5 | 
 6 | df_cars = pd.read_csv("tmp")
 7 | X, y = df_cars[['COUNTS']], df_cars['BYTES']
 8 | 
 9 | dt = DecisionTreeRegressor(max_depth=3, criterion="mae")
10 | dt.fit(X, y)
11 | 
12 | fig = plt.figure()
13 | ax = fig.gca()
14 | rtreeviz_univar(dt, X, y, 'COUNTS', 'BYTES', ax=ax)
15 | plt.show()
16 | 


--------------------------------------------------------------------------------
/kmeans_tbb.cpp:
--------------------------------------------------------------------------------
 1 | #define NOMINMAX  // Required on Windows
 2 | #include <limits>
 3 | #include <cstdlib>
 4 | #include <cstdio>
 5 | #include "tbb/tbb.h"
 6 | 
 7 | #include "kmeans.h"
 8 | #include "view.h"
 9 | #include "reduce_local_counts_to_global_count.h"
10 | #include "reduce_local_sums_to_global_sum.h"
11 | #include "reduce_min_ind.h"
12 | 
13 | namespace tbb_asura {
14 | 
15 | #include "kmeans_tbb.h"
16 | 
17 | } // namespace tbb_example
18 | 


--------------------------------------------------------------------------------
/Makefile:
--------------------------------------------------------------------------------
 1 | LIBROOT 	= /usr/local/lib
 2 | CXXFLAGS 	= -Wall -Wextra -Wpedantic -std=c++11 -Os -O2 -fpermissive -lpthread -ltbb
 3 | LDFLAGS 	= -L$(LIBROOT) -I$(LIBROOT) -ltins -lpcap -L/usr/local/lib/
 4 | 
 5 | PROGS 	= asura
 6 | 
 7 | UNAME := $(shell uname -s)
 8 | ifeq ($(UNAME), $(filter $(UNAME), Darwin FreeBSD))
 9 | 	CC = clang
10 | else
11 | 	CC = gcc
12 | endif
13 | 
14 | 
15 | $(PROGS): asura23.cpp kmeans_tbb.cpp repair_empty_clusters.cpp
16 | 	$(CXX) -o $@ $^ $(CXXFLAGS) $(LDFLAGS)
17 | 	strip -s $(PROGS)
18 | clean:
19 | 	rm -f $(PROGS) *.o
20 | 


--------------------------------------------------------------------------------
/reduced2-bak:
--------------------------------------------------------------------------------
 1 | 298,388890
 2 | 64,23042
 3 | 20,9818
 4 | 546,82336
 5 | 10,646
 6 | 8,416
 7 | 18,1980
 8 | 14,820
 9 | 8,1232
10 | 14,2470
11 | 2,165
12 | 6,312
13 | 30,10468
14 | 4,762
15 | 552,38934
16 | 16,928
17 | 38,10162
18 | 4,278
19 | 12,1170
20 | 30,3534
21 | 392,429312
22 | 10,1170
23 | 24,8646
24 | 6,3300
25 | 4,232
26 | 2,202
27 | 106,99562
28 | 432,596200
29 | 458,457050
30 | 10,5610
31 | 12,1150
32 | 16,4174
33 | 52,14422
34 | 22,6986
35 | 8,4576
36 | 542,109432
37 | 352,103728
38 | 96,45634
39 | 562,39616
40 | 4,208
41 | 12,7778
42 | 30,2740
43 | 12,1848
44 | 574,39350
45 | 230,40244
46 | 10,1546
47 | 2,80
48 | 82,19624
49 | 28,3572
50 | 422,53020
51 | 


--------------------------------------------------------------------------------
/sum_and_count.h:
--------------------------------------------------------------------------------
 1 | // Structure that holds a sum of points and the number of points 
 2 | // accumulated into that sum, so that the mean can be computed.
 3 | struct sum_and_count {
 4 |     sum_and_count() : sum(), count(0) {}
 5 |     point sum;
 6 |     size_t count;
 7 |     void clear() {
 8 |         sum = point();
 9 |         count = 0;
10 |     }
11 |     void tally( const point& p ) {
12 |         sum += p;
13 |         ++count;
14 |     }
15 |     point mean() const {
16 |         return sum/count;
17 |     }
18 |     void operator+=( const sum_and_count& other ) {
19 |         sum += other.sum;
20 |         count += other.count;
21 |     };
22 | };


--------------------------------------------------------------------------------
/repair_empty_clusters.cpp:
--------------------------------------------------------------------------------
 1 | #include "kmeans.h"
 2 | 
 3 | void fix_empty_clusters( size_t n, const point points[], cluster_id id[], size_t k, point centroid[], sum_and_count sum[] ) {
 4 |     for( size_t j=0; j<k; ++j ) {
 5 |         if( sum[j].count==0 ) {
 6 | //#if HAVE_CILKPLUS
 7 | //            size_t farthest = __sec_reduce_max_ind(distance2(centroid[id[0:n]],points[0:n]));
 8 | //#else
 9 |             size_t farthest = -1;
10 |             float maxd = -1;
11 |             for( int i=0; i<n; ++i ) {
12 |                 float d = distance2(centroid[id[i]],points[i]);
13 |                 if( d>maxd  ) {
14 |                     maxd = d;
15 |                     farthest = i;
16 |                 }
17 |             }
18 | //#endif
19 |             id[farthest] = j;
20 |             sum[j].count = 1;
21 |             sum[j].sum = points[farthest];
22 |         } 
23 |     }
24 | }
25 | 


--------------------------------------------------------------------------------
/kmeans.h:
--------------------------------------------------------------------------------
 1 | #include <cstddef>
 2 | 
 3 | struct point {
 4 |     float x, y;
 5 |     void operator+=( const point& right ) {
 6 |         x += right.x;
 7 |         y += right.y;
 8 |     }
 9 |     point operator/( size_t count ) const {
10 |         point p = *this;
11 |         p.x /= count;
12 |         p.y /= count;
13 |         return p;
14 |     }
15 | };
16 | 
17 | inline float distance2( const point& a, const point& b ) {
18 |     float dx = a.x-b.x;
19 |     float dy = a.y-b.y;
20 |     return dx*dx+dy*dy;
21 | }
22 | 
23 | #include "sum_and_count.h"
24 | 
25 | typedef unsigned short cluster_id;
26 | 
27 | void fix_empty_clusters( size_t n, const point points[], cluster_id id[], size_t k, point centroid[], sum_and_count sum[] );
28 | 
29 | /*
30 | namespace cilk_example {
31 |     void compute_k_means( size_t n, const point points[], size_t k, cluster_id id[], point centroid[] ); 
32 | }
33 | */
34 | 
35 | namespace tbb_asura {
36 |     void do_k_means( size_t n, const point points[], size_t k, cluster_id id[], point centroid[] ); 
37 | }
38 | 


--------------------------------------------------------------------------------
/kmeans_tbb.h:
--------------------------------------------------------------------------------
 1 | void do_k_means( size_t n, const point points[], size_t k, cluster_id id[], point centroid[] ) {
 2 | 
 3 |     tls_type tls([&]{return k;}); 
 4 |     view global(k);
 5 |     
 6 |     tbb::parallel_for(
 7 |         tbb::blocked_range<size_t>(0,n),
 8 |         [=,&tls,&global]( tbb::blocked_range<size_t> r ) {
 9 |             view& v = tls.local();
10 |             for( size_t i=r.begin(); i!=r.end(); ++i ) {
11 |                 id[i] = i % k;  
12 |                 // Peeled "Sum step"
13 |                 v.array[id[i]].tally(points[i]);
14 |             }
15 |         }
16 |     );
17 |     
18 |     size_t change;
19 |     do {
20 |         convert_local_sums_to_global_sum( k, tls, global );
21 |         fix_empty_clusters( n, points, id, k, centroid, global.array );
22 |         
23 |         for( size_t j=0; j<k; ++j ) {
24 |             centroid[j] = global.array[j].mean();
25 |             global.array[j].clear();
26 |         }
27 | 
28 |         tbb::parallel_for(
29 |             tbb::blocked_range<size_t>(0,n),
30 |             [=,&tls,&global]( tbb::blocked_range<size_t> r ) {
31 |                 view& v = tls.local();
32 |                 for( size_t i=r.begin(); i!=r.end(); ++i ) {
33 |                     cluster_id j = calc_shortest_index(centroid, k , points[i]); 
34 |                     if( j!=id[i] ) {
35 |                         id[i] = j;
36 |                         ++v.change;
37 |                     }
38 |                     v.array[j].tally(points[i]);
39 |                 }
40 |             }
41 |         );
42 |         
43 |         convert_local_counts_to_global_count( tls, global );
44 |     } while( global.change!=0 );
45 | }
46 | 


--------------------------------------------------------------------------------
/tmp:
--------------------------------------------------------------------------------
  1 | 1,40
  2 | 1459,84866
  3 | 2,80
  4 | 2833,241780
  5 | 2988,187275
  6 | 7,448
  7 | 460,25606
  8 | 633,34543
  9 | 49328,2.21751e+06
 10 | 12,720
 11 | 12,720
 12 | 113748,9.86882e+06
 13 | 42,1766
 14 | 22,1180
 15 | 9,468
 16 | 7,448
 17 | 36654,1.48253e+06
 18 | 12,720
 19 | 2004,120240
 20 | 2301,198771
 21 | 562,24394
 22 | 87613,6.07037e+06
 23 | 4785,322049
 24 | 7737,702600
 25 | 8,336
 26 | 5,200
 27 | 5,320
 28 | 2027,368418
 29 | 144147,6.38727e+06
 30 | 18552,949675
 31 | 36,1682
 32 | 6057,313479
 33 | 3055,231404
 34 | 35,2240
 35 | 30,4740
 36 | 2005,120300
 37 | 2011,120660
 38 | 2871,208006
 39 | 203,8504
 40 | 11,704
 41 | 6434,354703
 42 | 34,1420
 43 | 2,754
 44 | 5,200
 45 | 8,2453
 46 | 1260,50400
 47 | 512,27328
 48 | 2906,159921
 49 | 3051,168545
 50 | 36734,1.61892e+06
 51 | 2,88
 52 | 8,416
 53 | 1,152
 54 | 8,6356
 55 | 336,13836
 56 | 319,16338
 57 | 3357,149826
 58 | 49002,2.05695e+06
 59 | 4090,280469
 60 | 5,320
 61 | 6,384
 62 | 13,740
 63 | 3,252
 64 | 2021,89300
 65 | 18566,2.40658e+06
 66 | 1543,118895
 67 | 2530,158477
 68 | 13,892
 69 | 7364,493109
 70 | 6577,346672
 71 | 37800,1.71083e+06
 72 | 36197,1.44788e+06
 73 | 1,40
 74 | 8,336
 75 | 1814,106880
 76 | 8,336
 77 | 3269,131080
 78 | 8,336
 79 | 142,40720
 80 | 1413,59633
 81 | 1201,54730
 82 | 5052,303120
 83 | 8,336
 84 | 5,320
 85 | 4,168
 86 | 2590,158790
 87 | 922,47954
 88 | 961,55006
 89 | 515,31285
 90 | 2663,117364
 91 | 5021,200840
 92 | 7,448
 93 | 38,2110
 94 | 1357,118210
 95 | 57283,2.3587e+06
 96 | 18704,2.40904e+06
 97 | 4389,305176
 98 | 37350,1.59014e+06
 99 | 5,1840
100 | 45555,2.00446e+06
101 | 5,320
102 | 1674,90563
103 | 3976,218128
104 | 3227,129424
105 | 12,720
106 | 3259,172496
107 | 212757,1.90418e+07
108 | 4513,184602
109 | 2551,155293
110 | 48,2486
111 | 12,720
112 | 56,3584
113 | 32821,1.44452e+06
114 | 8,336
115 | 47,2246
116 | 4,240
117 | 861,48490
118 | 103,10746
119 | 474,25486
120 | 68694,2.74806e+06
121 | 93578,6.22639e+06
122 | 4273,292429
123 | 16,960
124 | 6,240
125 | 1168,47086
126 | 1294,59918
127 | 12,768
128 | 35,2240
129 | 6104,382255
130 | 5669,345137
131 | 36205,1.59318e+06
132 | 2988,155330
133 | 18,757
134 | 584,32476
135 | 141,7598
136 | 3506,174234
137 | 8948,490344
138 | 36810,1.51835e+06
139 | 8,512
140 | 173,14420
141 | 3193,191540
142 | 1156,46479
143 | 12,504
144 | 2,128
145 | 7760,514978
146 | 32811,1.31244e+06
147 | 12,504
148 | 7,280
149 | 4,200
150 | 1,44
151 | 233,10236
152 | 58105,2.32422e+06
153 | 8,787
154 | 36979,1.67362e+06
155 | 54,3060
156 | 5256,389060
157 | 6,384
158 | 3,363
159 | 18707,2.40858e+06
160 | 1532,73227
161 | 12,504
162 | 1,84
163 | 2015,120900
164 | 1490,112296
165 | 2720,109012
166 | 3247,144870
167 | 58492,2.57383e+06
168 | 3671,223235
169 | 3377,201886
170 | 3337,200544
171 | 6955,482422
172 | 4182,167508
173 | 5,260
174 | 2028,89616
175 | 107,21192
176 | 3223,191110
177 | 36736,1.62466e+06
178 | 18535,776730
179 | 24600,1.05377e+06
180 | 68540,3.01598e+06
181 | 7479,505822
182 | 6688,400927
183 | 7,448
184 | 6,320
185 | 2,135
186 | 143,7594
187 | 4751,283937
188 | 2009,120540
189 | 9366,562608
190 | 4473,267712
191 | 2,428
192 | 3155,168995
193 | 2736,158063
194 | 651,36112
195 | 27,1080
196 | 4271,219702
197 | 2357,156725
198 | 14,768
199 | 114,7200
200 | 1186,56107
201 | 1221,55876
202 | 6334,422758
203 | 7,448
204 | 2007,120420
205 | 8034,545426
206 | 14,696
207 | 36421,1.4571e+06
208 | 6,240
209 | 2,428
210 | 1730,107629
211 | 157,8312
212 | 2339,99378
213 | 70813,4.2474e+06
214 | 4339,273937
215 | 12,504
216 | 11,484
217 | 12,1667
218 | 18548,950511
219 | 4114,229268
220 | 2,210
221 | 4348,261054
222 | 20,880
223 | 112,7896
224 | 44,2320
225 | 8,336
226 | 3146,136999
227 | 5498,252529
228 | 42,2304
229 | 8,320
230 | 92,4004
231 | 18565,946455
232 | 1158,54169
233 | 301,13452
234 | 12,504
235 | 6379,372480
236 | 48127,2.11788e+06
237 | 6,384
238 | 354,15754
239 | 3647,206914
240 | 3217,146468
241 | 28,1228
242 | 6721,390804
243 | 7020,450738
244 | 867,38144
245 | 


--------------------------------------------------------------------------------
/dtreeviz/tmp:
--------------------------------------------------------------------------------
  1 | COUNTS,BYTES
  2 | 1,40
  3 | 1459,84866
  4 | 2,80
  5 | 2833,241780
  6 | 2988,187275
  7 | 7,448
  8 | 460,25606
  9 | 633,34543
 10 | 49328,2.21751e+06
 11 | 12,720
 12 | 12,720
 13 | 113748,9.86882e+06
 14 | 42,1766
 15 | 22,1180
 16 | 9,468
 17 | 7,448
 18 | 36654,1.48253e+06
 19 | 12,720
 20 | 2004,120240
 21 | 2301,198771
 22 | 562,24394
 23 | 87613,6.07037e+06
 24 | 4785,322049
 25 | 7737,702600
 26 | 8,336
 27 | 5,200
 28 | 5,320
 29 | 2027,368418
 30 | 144147,6.38727e+06
 31 | 18552,949675
 32 | 36,1682
 33 | 6057,313479
 34 | 3055,231404
 35 | 35,2240
 36 | 30,4740
 37 | 2005,120300
 38 | 2011,120660
 39 | 2871,208006
 40 | 203,8504
 41 | 11,704
 42 | 6434,354703
 43 | 34,1420
 44 | 2,754
 45 | 5,200
 46 | 8,2453
 47 | 1260,50400
 48 | 512,27328
 49 | 2906,159921
 50 | 3051,168545
 51 | 36734,1.61892e+06
 52 | 2,88
 53 | 8,416
 54 | 1,152
 55 | 8,6356
 56 | 336,13836
 57 | 319,16338
 58 | 3357,149826
 59 | 49002,2.05695e+06
 60 | 4090,280469
 61 | 5,320
 62 | 6,384
 63 | 13,740
 64 | 3,252
 65 | 2021,89300
 66 | 18566,2.40658e+06
 67 | 1543,118895
 68 | 2530,158477
 69 | 13,892
 70 | 7364,493109
 71 | 6577,346672
 72 | 37800,1.71083e+06
 73 | 36197,1.44788e+06
 74 | 1,40
 75 | 8,336
 76 | 1814,106880
 77 | 8,336
 78 | 3269,131080
 79 | 8,336
 80 | 142,40720
 81 | 1413,59633
 82 | 1201,54730
 83 | 5052,303120
 84 | 8,336
 85 | 5,320
 86 | 4,168
 87 | 2590,158790
 88 | 922,47954
 89 | 961,55006
 90 | 515,31285
 91 | 2663,117364
 92 | 5021,200840
 93 | 7,448
 94 | 38,2110
 95 | 1357,118210
 96 | 57283,2.3587e+06
 97 | 18704,2.40904e+06
 98 | 4389,305176
 99 | 37350,1.59014e+06
100 | 5,1840
101 | 45555,2.00446e+06
102 | 5,320
103 | 1674,90563
104 | 3976,218128
105 | 3227,129424
106 | 12,720
107 | 3259,172496
108 | 212757,1.90418e+07
109 | 4513,184602
110 | 2551,155293
111 | 48,2486
112 | 12,720
113 | 56,3584
114 | 32821,1.44452e+06
115 | 8,336
116 | 47,2246
117 | 4,240
118 | 861,48490
119 | 103,10746
120 | 474,25486
121 | 68694,2.74806e+06
122 | 93578,6.22639e+06
123 | 4273,292429
124 | 16,960
125 | 6,240
126 | 1168,47086
127 | 1294,59918
128 | 12,768
129 | 35,2240
130 | 6104,382255
131 | 5669,345137
132 | 36205,1.59318e+06
133 | 2988,155330
134 | 18,757
135 | 584,32476
136 | 141,7598
137 | 3506,174234
138 | 8948,490344
139 | 36810,1.51835e+06
140 | 8,512
141 | 173,14420
142 | 3193,191540
143 | 1156,46479
144 | 12,504
145 | 2,128
146 | 7760,514978
147 | 32811,1.31244e+06
148 | 12,504
149 | 7,280
150 | 4,200
151 | 1,44
152 | 233,10236
153 | 58105,2.32422e+06
154 | 8,787
155 | 36979,1.67362e+06
156 | 54,3060
157 | 5256,389060
158 | 6,384
159 | 3,363
160 | 18707,2.40858e+06
161 | 1532,73227
162 | 12,504
163 | 1,84
164 | 2015,120900
165 | 1490,112296
166 | 2720,109012
167 | 3247,144870
168 | 58492,2.57383e+06
169 | 3671,223235
170 | 3377,201886
171 | 3337,200544
172 | 6955,482422
173 | 4182,167508
174 | 5,260
175 | 2028,89616
176 | 107,21192
177 | 3223,191110
178 | 36736,1.62466e+06
179 | 18535,776730
180 | 24600,1.05377e+06
181 | 68540,3.01598e+06
182 | 7479,505822
183 | 6688,400927
184 | 7,448
185 | 6,320
186 | 2,135
187 | 143,7594
188 | 4751,283937
189 | 2009,120540
190 | 9366,562608
191 | 4473,267712
192 | 2,428
193 | 3155,168995
194 | 2736,158063
195 | 651,36112
196 | 27,1080
197 | 4271,219702
198 | 2357,156725
199 | 14,768
200 | 114,7200
201 | 1186,56107
202 | 1221,55876
203 | 6334,422758
204 | 7,448
205 | 2007,120420
206 | 8034,545426
207 | 14,696
208 | 36421,1.4571e+06
209 | 6,240
210 | 2,428
211 | 1730,107629
212 | 157,8312
213 | 2339,99378
214 | 70813,4.2474e+06
215 | 4339,273937
216 | 12,504
217 | 11,484
218 | 12,1667
219 | 18548,950511
220 | 4114,229268
221 | 2,210
222 | 4348,261054
223 | 20,880
224 | 112,7896
225 | 44,2320
226 | 8,336
227 | 3146,136999
228 | 5498,252529
229 | 42,2304
230 | 8,320
231 | 92,4004
232 | 18565,946455
233 | 1158,54169
234 | 301,13452
235 | 12,504
236 | 6379,372480
237 | 48127,2.11788e+06
238 | 6,384
239 | 354,15754
240 | 3647,206914
241 | 3217,146468
242 | 28,1228
243 | 6721,390804
244 | 7020,450738
245 | 867,38144
246 | 


--------------------------------------------------------------------------------
/README.md:
--------------------------------------------------------------------------------
  1 | # Massive Pcap Analyzer with Massive Multithreading (DEF CON 26)
  2 | 
  3 | “The universe is not complicated, there's just a lot of it.” - Richard Feynman
  4 | 
  5 | <img src="asura.jpeg" width=200 height=200>
  6 | 
  7 | <pre>
  8 | Required: 
  9 | libtins
 10 | http://libtins.github.io/ 
 11 | Intel TBB
 12 | 
 13 | # apt install libpcap-dev
 14 | # apt install libtins-dev
 15 | 
 16 | quick test:
 17 | # ./asura pcaps
 18 | </pre>
 19 | With the rapid increase of network traffic, we are suffering the long processing time in coping with large packet capture file in many aspects such as network diagnosis, incidents response and CTF. Wireshark and Linux grep command are basically run in single-process mode and are thus not good at analyzing tens or hundreds gigabytes of packet dump in reasonable processing time. On the other hand, high performance devices such as GPGPU are still bulky and expensive for running on our laptop.
 20 | 
 21 | Asura One is a portable anomaly detector of large packet capture file (around 100GB) leveraged by raw thread (Posix Pthreads) and highly concurrent container of Intel TBB (Threading Building Block). Asura One is designed for the lightweight anomaly detection of large PCAP file ranging from tens to hundreds million packets with the parallel clustering algorithm. Asura One takes two steps: extracting the flow vector {(srcIP,dstIP), X,Y..} and the clustering of which the outputs are {(srcIP,dstIP), Anomaly_Score}. The flow vector {(srcIP,dstIP), X,Y.. } is implemented as the concurrent hashmap of Intel TBB offering a much higher level of concurrency. Compared with the conventional C++ or Java hashmap such as STL, Intel TBB concurrent hashmap adopts lock-free techniques and is thus so faster. Another thrust of Asura One is leveraging task-decomposition by raw threads. Asura's task decomposition based multithreading can handle various types and lengths of IP packets and transforms these into well-organized input of {(srcIP,dstIP), X,Y..}. Asura One with its outputs of {(srcIP,dstIP), Anomaly_Score} should be a good partner (as pre-filter) of Wireshark. 
 22 | 
 23 | In experiment, Asura One has processed about 400,000,000 packets with 83GB PCAP file size in about 30 minutes on laptop PC. Packet dump was drastically reduced from 83GB to 7MB with about 110,000 aggregated flow vectors which are feasible for the further careful inspection by Wireshark.
 24 | 
 25 | <u>1.Lisence: Asura One is now released as open source under MIT license. </u>
 26 | 
 27 | 2.Compiler and libraries: gcc version 7.3.0 (Ubuntu 7.3.0), Posix Pthreads and Intel TBB.
 28 | 
 29 | 3.Input: PCAP files in the directory 
 30 | 
 31 | 4.Usage: ./asura DIRECTORY_NAME
 32 | 
 33 | <pre>
 34 | 	# mkdir pcap
 35 |         # cd pcap 
 36 | 	# wget -r -l 4 https://download.netresec.com/pcap/ists-12/2015-03-07/
 37 | 	# cp -r download.netresec.com/pcap/ists-12/2015-03-07/* .
 38 | 	# rm -rf download.netresec.com
 39 |         # cd ..
 40 | 	# ./build-asura.sh 
 41 |         # ./asura pcap
 42 | </pre>
 43 | 
 44 | 5.Output: {<sourceIP, destinationIP>, Anomaly_Score}
 45 | 
 46 | <pre>
 47 | 	sourceIP:sourcePort => destinationIP:destinationPort, clusterID, points(X,Y), percentage
 48 |         173.194.70.94:173.194.70.94=>10.0.2.16:10.0.2.16, clusterID:1, data(954,478932), 22%
 49 |         10.0.2.16:10.0.2.16=>195.169.125.228:195.169.125.228, clusterID:2, data(141,7144), 55%
 50 | 　　　   X.X.X.X,Z.Z.Z.Z -> 5 (73445,48),288,118644,0.242743%
 51 | </pre>
 52 | 　
 53 | If precentage is low, the pair(srcIP,srcPort->dstIP,dstPort) could be anomaly. 
 54 | 
 55 | 6.Procedure 1: extracting flow vector {<srcIP, dstIP>, X, Y}
 56 | <pre>
 57 | 	Container: typedef concurrent_hash_map<unsigned long long, int, HashCompare>
 58 | </pre>
 59 | 
 60 | <pre>
 61 | 	Main loop: 
 62 |   	pthread_create(&master, NULL, (void*)master_func, (void*)&targ[0]);
 63 |     	  for (i = 1; i < thread_num; ++i) { 
 64 |         	targ[i].id = i;
 65 |        		pthread_create(&worker[i],NULL,(void*)worker_func,(void*)&targ[i]); }
 66 |     	  for (i = 1; i < thread_num; ++i) 
 67 |  	       pthread_join(worker[i], NULL);
 68 | </pre>
 69 | 
 70 | 7.Procedure 2: calculating anomaly score {<srcIP, dstIP>, Anomaly_Score}
 71 | <pre>
 72 | 	  Main loop (K-Means):
 73 | 	  tbb::parallel_for(
 74 |           tbb::blocked_range<size_t>(0,n),
 75 |           [=,&tls,&global]( tbb::blocked_range<size_t> r ) {
 76 |               view& v = tls.local();
 77 |               for( size_t i=r.begin(); i!=r.end(); ++i ) {
 78 |                   cluster_id j = calc_shortest_index(centroid, k , points[i]); 
 79 |                   if( j!=id[i] ) {
 80 |                       id[i] = j;
 81 |                       ++v.change;
 82 |                   }
 83 |                   v.array[j].tally(points[i]);
 84 |               }
 85 |             }
 86 | </pre>
 87 | 
 88 | 8. Processing time: about 30 minutes for about 400,000,000 packets
 89 | 
 90 | <hr>
 91 | 2022-01-11
 92 | libpcap and libtins are required.
 93 | 
 94 | <pre>
 95 | # apt install libpcap-dev
 96 | # git clone https://github.com/mfontanini/libtins.git
 97 | # cd libtins/
 98 | # apt-get install libpcap-dev libssl-dev cmake
 99 | # mkdir build
100 | # cd build/
101 | # cmake ../ -DLIBTINS_ENABLE_CXX11=1
102 | # make install
103 | # ldconfig
104 | </pre>
105 | 


--------------------------------------------------------------------------------
/test_kmeans.cpp:
--------------------------------------------------------------------------------
  1 | #include "kmeans.h"
  2 | #include <cstdlib>
  3 | #include <cstdio>
  4 | #include <cstring>
  5 | #include <cmath>
  6 | 
  7 | #include <stdio.h>
  8 | #include <stdlib.h>
  9 | #include <stdbool.h>
 10 | #include <string.h>
 11 | #include <sys/types.h>
 12 | #include <sys/stat.h>
 13 | #include <unistd.h>
 14 | #include <dirent.h>
 15 | #include <pthread.h>
 16 | #include <alloca.h>
 17 | #include "timer.h"
 18 | 
 19 | #include <regex.h>
 20 | 
 21 | #include<sys/ioctl.h>
 22 | #include<arpa/inet.h>
 23 | #include<sys/socket.h>
 24 | #include<linux/if.h>
 25 | #include<net/ethernet.h>
 26 | #include<netpacket/packet.h>
 27 | #include<netinet/if_ether.h>
 28 | #include<netinet/ip.h>
 29 | #include<netinet/ip6.h>
 30 | #include<netinet/ip_icmp.h>
 31 | #include<netinet/icmp6.h>
 32 | #include<netinet/tcp.h>
 33 | #include<netinet/udp.h>
 34 | 
 35 | #include <map>
 36 | #include <iostream>
 37 | #include <fstream>
 38 | #include <string>
 39 | #include <vector>
 40 | #include <bitset>
 41 | 
 42 | #include <cstdio>
 43 | #include <cctype>
 44 | #include <iostream>
 45 | 
 46 | #include <boost/tokenizer.hpp>
 47 | 
 48 | using namespace std;
 49 | // using namespace tbb;
 50 | 
 51 | // Square root of number of intended centroids
 52 | /*
 53 | const size_t SQRT_K = 4;
 54 | const size_t K = SQRT_K*SQRT_K;
 55 | const size_t N = 533;
 56 | 
 57 | point points[N];
 58 | point centroid[K];
 59 | cluster_id id[N];
 60 | */
 61 | 
 62 | std::vector < std::vector< std::string > > parse_csv(const char* filepath)
 63 | {
 64 |     std::vector< std::vector< std::string > > cells;
 65 |     std::string line;
 66 |     std::ifstream ifs(filepath);
 67 | 
 68 |     // csvを走査
 69 |     while (std::getline(ifs, line)) {
 70 | 
 71 |         std::vector< std::string > data;
 72 | 
 73 |         // 1行を走査
 74 |         boost::tokenizer< boost::escaped_list_separator< char > > tokens(line);
 75 |         for (const std::string& token : tokens) {
 76 |             data.push_back(token);
 77 |         }
 78 | 
 79 |         // 1行読み込んだ結果を入れる
 80 |         cells.push_back(data);
 81 |     }
 82 | 
 83 |     return cells;
 84 | }
 85 | 
 86 | /*
 87 | int main(int argc, char* argv[]) {
 88 | 
 89 |   const size_t SQRT_K = 4;
 90 |   const size_t K = SQRT_K*SQRT_K;
 91 |   const size_t N = 533;
 92 | 
 93 |   point* points;
 94 |   points = (struct point *)malloc(N*sizeof(struct point));
 95 | 
 96 |   point* centroid;
 97 |   centroid = (struct point *)malloc(N*sizeof(struct point));
 98 | 
 99 |   cluster_id* id;
100 |   id = (unsigned short *)malloc(N*sizeof(unsigned short));
101 | 
102 |   std::vector<string> pair; 
103 |     
104 |     printf("Testing TBB kmeans algorithm...\n");
105 | 
106 |     const auto cells = parse_csv(argv[1]);
107 |   
108 |     int counter = 0;
109 |     for (const auto& rows : cells) {
110 | 
111 |       point& p = points[counter];
112 |       
113 |       p.x = std::stof(rows[1]);
114 |       p.y = std::stof(rows[2]);      
115 |       pair.push_back(rows[0]);
116 |       
117 |       counter = counter + 1;
118 |     }
119 |     
120 |     tbb_example::compute_k_means( N, points, K, id, centroid );
121 | 
122 | #if 1
123 | 
124 |     int* counts;
125 |     counts = (int *)malloc(K*sizeof(int));
126 | 
127 |     for( size_t i=0; i<N; ++i ) {
128 |       counts[id[i]] = 0;
129 |     }
130 |     
131 |     for( size_t i=0; i<N; ++i ) {
132 |       counts[id[i]]++;
133 |     }
134 |       
135 |     for( size_t i=0; i<N; ++i ) {
136 | 
137 | 	    unsigned long pair_long = std::stoul(pair[i]);
138 |       
139 | 	    std::bitset<64> bset_pair = std::bitset<64>(pair_long);   
140 | 	    string bset_pair_string = bset_pair.to_string();
141 | 	    
142 | 	    std::string ip1 = bset_pair_string.substr(0, 8);
143 | 	    std::string ip2 = bset_pair_string.substr(8, 8);
144 | 	    std::string ip3 = bset_pair_string.substr(16, 8);
145 | 	    std::string ip4 = bset_pair_string.substr(24, 8);
146 | 
147 | 	    std::bitset<8> ip1_bset = std::bitset<8>(ip1);
148 | 	    std::bitset<8> ip2_bset = std::bitset<8>(ip2);
149 | 	    std::bitset<8> ip3_bset = std::bitset<8>(ip3);
150 | 	    std::bitset<8> ip4_bset = std::bitset<8>(ip4);
151 | 
152 | 	    std::cout << ip1_bset.to_ulong() << "." << ip2_bset.to_ulong() << "." << ip3_bset.to_ulong() << "." << ip4_bset.to_ulong() << ",";
153 | 	    
154 | 	    std::string ip5 = bset_pair_string.substr(32, 8);
155 | 	    std::string ip6 = bset_pair_string.substr(40, 8);
156 | 	    std::string ip7 = bset_pair_string.substr(48, 8);
157 | 	    std::string ip8 = bset_pair_string.substr(56, 8);
158 | 
159 | 	    std::bitset<8> ip5_bset = std::bitset<8>(ip5);
160 | 	    std::bitset<8> ip6_bset = std::bitset<8>(ip6);
161 | 	    std::bitset<8> ip7_bset = std::bitset<8>(ip7);
162 | 	    std::bitset<8> ip8_bset = std::bitset<8>(ip8);
163 | 
164 | 	    std::cout << ip5_bset.to_ulong() << "." << ip6_bset.to_ulong() << "." << ip7_bset.to_ulong() << "." << ip8_bset.to_ulong() << " -> ";
165 | 
166 | 	    float percent = (float)counts[id[i]]/(float)N;
167 | 	    printf("%d (%g %g) counts %d / %d [%f%] \n",id[i],points[i].x,points[i].y, counts[id[i]], N, percent); 
168 |     }
169 | #endif
170 | #if 1
171 |     printf("centroids = ");
172 |     for( size_t j=0; j<K; ++j ) {
173 |         printf("(%g %g)",centroid[j].x,centroid[j].y);
174 |     }
175 |     printf("\n");
176 | #endif
177 | 
178 |     return 0;
179 | }
180 | */
181 | 


--------------------------------------------------------------------------------
/asura23.cpp:
--------------------------------------------------------------------------------
  1 | #include <pcap.h>
  2 | #include <tins/tins.h>
  3 | 
  4 | #include <stdio.h>
  5 | #include <stdlib.h>
  6 | #include <stdbool.h>
  7 | #include <string.h>
  8 | #include <sys/types.h>
  9 | #include <sys/stat.h>
 10 | #include <unistd.h>
 11 | #include <dirent.h>
 12 | #include <pthread.h>
 13 | #include <alloca.h>
 14 | 
 15 | #include <cstdio>
 16 | #include <cstdlib>
 17 | #include <cmath>
 18 | #include <ctime>
 19 | #include <iostream>
 20 | #include <fstream>
 21 | #include <sstream>
 22 | 
 23 | #include "timer.h" 
 24 | #include "kmeans.h"
 25 | 
 26 | #define MB (1048576.0)
 27 | #define UNUSE(x) (void)(x)
 28 | 
 29 | using Tins::SnifferConfiguration;
 30 | using Tins::FileSniffer;
 31 | using Tins::IP;
 32 | using Tins::DNS;
 33 | using Tins::TCP;
 34 | using Tins::UDP;
 35 | using Tins::RawPDU;
 36 | using Tins::PDU;
 37 | using Tins::Packet;
 38 | using namespace Tins;
 39 | using namespace std;
 40 | 
 41 | #include "tbb/task_scheduler_init.h"
 42 | #include "tbb/concurrent_hash_map.h"
 43 | #include "tbb/blocked_range.h"
 44 | #include "tbb/parallel_for.h"
 45 | 
 46 | using namespace std;
 47 | using namespace tbb;
 48 | 
 49 | #define N 1
 50 | #define WORKER_THREAD_NUM N
 51 | #define MAX_QUEUE_NUM N
 52 | #define END_MARK_FNAME   "///"
 53 | #define END_MARK_FLENGTH 3
 54 | 
 55 | #define DISP_FREQ 100000
 56 | 
 57 | /* srcIP, destIP */
 58 | typedef struct _addrpair {
 59 |   map<string, string> m;
 60 |   pthread_mutex_t mutex;
 61 | } addrpair_t;
 62 | addrpair_t addrpair;
 63 | 
 64 | struct HashCompare {
 65 |   static size_t hash( unsigned long long x ) {
 66 |     return (size_t)x;
 67 |     }
 68 |   static bool equal( unsigned long long x, unsigned long long y ) {
 69 |         return x==y;
 70 |     }
 71 | };
 72 | 
 73 | typedef concurrent_hash_map<std::string, int> CharTable;
 74 | static CharTable table;
 75 | 
 76 | typedef concurrent_hash_map<std::string, int> CharTable2;
 77 | static CharTable2 table2;
 78 | 
 79 | typedef struct _result {
 80 |     int num;
 81 |     char* fname;
 82 |     pthread_mutex_t mutex;    
 83 | } result_t;
 84 | result_t result;
 85 | 
 86 | typedef struct _queue {
 87 |     char* fname[MAX_QUEUE_NUM];
 88 |     int flength[MAX_QUEUE_NUM];
 89 |     int rp, wp;
 90 |     int remain;
 91 |     pthread_mutex_t mutex;
 92 |     pthread_cond_t not_full;
 93 |     pthread_cond_t not_empty;
 94 | } queue_t;
 95 | 
 96 | typedef struct _thread_arg {
 97 |     int id;
 98 |     int cpuid;
 99 |     queue_t* q;
100 |     char* srchstr;
101 |     char* dirname;
102 |     int filenum;
103 | } thread_arg_t;
104 | 
105 | 
106 | bool callback(const PDU &pdu) {
107 |   // Parse the IP header
108 |   const IP &ip = pdu.rfind_pdu<IP>(); 
109 | 
110 |   // Parse the TCP or UDP header
111 |   const TCP &tcp = pdu.rfind_pdu<TCP>();
112 | 
113 |   string source_ip = ip.src_addr().to_string();
114 |   string source_port = std::to_string(tcp.sport());
115 | 
116 |   string dest_ip = ip.dst_addr().to_string();
117 |   string dest_port = std::to_string(tcp.dport());
118 | 
119 |   string info = source_ip + ":" + source_port + "=>" + dest_ip + ":" + dest_port;
120 |   
121 |   string info2 = source_ip + "=>" + dest_ip; 
122 |   
123 |   // cout << "IP:" << "[" << ip.tot_len() << "]" << ip.src_addr() << ":" << source_ip << ':' << tcp.sport() << ":" << " -> " << ip.dst_addr() << ':' << tcp.dport() << endl;
124 | 
125 |   CharTable::accessor a;
126 | 
127 |   table.insert(a, info2);
128 |   a->second += 1;
129 | 
130 |   CharTable2::accessor a2;
131 | 
132 |   table2.insert(a2, info2);
133 |   a2->second += int(ip.tot_len());     
134 | 
135 |   return true;
136 | }
137 | 
138 | 
139 | int traverse_file(char* filename, char* srchstr, int thread_id) {
140 | 
141 |     cout << filename << endl; 
142 | 
143 |         try {
144 |                 SnifferConfiguration config;
145 |                 config.set_filter("tcp or udp");
146 |                 config.set_promisc_mode(false);
147 |                 // config.set_snap_len(65536);
148 | 
149 |                 FileSniffer sniffer(filename, config);
150 |                 sniffer.sniff_loop(callback);
151 | 		
152 |                 return 0;
153 |         }
154 |         catch (std::exception& ex) {
155 |                 std::cerr << "[X] Error: " << ex.what() << '\n';
156 |                 return -1;
157 |         }
158 |     
159 |     return 0;
160 | }
161 | 
162 | void initqueue(queue_t* q) {
163 |     int i;
164 |     q->rp = q->wp = q->remain= 0;
165 |     for (i = 0; i < MAX_QUEUE_NUM; ++i) q->fname[i] = NULL;
166 |     pthread_mutex_init(&q->mutex,    NULL);
167 |     pthread_cond_init(&q->not_full,  NULL);
168 |     pthread_cond_init(&q->not_empty, NULL);
169 |     return;
170 | }
171 | 
172 | void enqueue(queue_t* q, char* path, int size) {
173 |   
174 |     pthread_mutex_lock(&q->mutex);
175 |     while (q->remain == MAX_QUEUE_NUM) {
176 |         pthread_cond_wait(&q->not_full, &q->mutex);
177 |     }
178 |     char** fname = (char**)&q->fname[q->wp];
179 |     if (*fname != NULL) free(*fname);
180 |     *fname = (char*)malloc(size);
181 |     strcpy(*fname, path);
182 |     q->flength[q->wp] = size;
183 |     q->wp++; q->remain++;
184 | 
185 |     if (q->wp == MAX_QUEUE_NUM) q->wp = 0;
186 | 
187 |     pthread_cond_signal(&q->not_empty);
188 |     pthread_mutex_unlock(&q->mutex);
189 |     return;
190 | }
191 | 
192 | void dequeue(queue_t* q, char** fname, int* flen) {
193 | 
194 |     pthread_mutex_lock(&q->mutex);
195 |     while (q->remain == 0) 
196 |         pthread_cond_wait(&q->not_empty, &q->mutex);
197 | 
198 |     *flen  = q->flength[q->rp];
199 |     if (*fname != NULL) free(*fname);
200 |     *fname = (char*)malloc(*flen);
201 |     strcpy(*fname, q->fname[q->rp]);
202 |     q->rp++; q->remain--;
203 |     if (q->rp == MAX_QUEUE_NUM) q->rp = 0;
204 |     pthread_cond_signal(&q->not_full);
205 |     pthread_mutex_unlock(&q->mutex);
206 |     if (strcmp(*fname,"")==0) printf("rp=%d\n", q->rp-1);
207 |     return;
208 | }
209 | 
210 | int traverse_dir_thread(queue_t* q, char* dirname) {
211 |     static int cnt = 0;
212 |     struct dirent* dent;
213 |     DIR* dd = opendir(dirname);
214 | 
215 |     if (dd == NULL) {
216 |         printf("Could not open the directory %s\n", dirname); return 0;
217 |     }
218 | 
219 |     while ((dent = readdir(dd)) != NULL) {
220 |         if (strncmp(dent->d_name, ".",  2) == 0) continue;
221 |         if (strncmp(dent->d_name, "..", 3) == 0) continue;	
222 | 
223 |         int size = strlen(dirname) + strlen(dent->d_name) + 2;
224 | #if 0
225 |         char* path = (char*)malloc(size);
226 |         sprintf(path, "%s/%s", dirname, dent->d_name);
227 | 
228 |         struct stat fs;
229 |         if (stat(path, &fs) < 0)
230 |             continue;
231 |         else {
232 |             if (S_ISDIR(fs.st_mode))
233 |                 traverse_dir_thread(q, path);
234 |             else if (S_ISREG(fs.st_mode)) {
235 |                 enqueue(q, path, size);
236 |                 cnt++;
237 |             }
238 |         }
239 | #else
240 |         {
241 |             char* path = (char*)alloca(size);
242 |             sprintf(path, "%s/%s", dirname, dent->d_name);
243 | 
244 |             struct stat fs;
245 |             if (stat(path, &fs) < 0)
246 |                 continue;
247 |             else {
248 |                 if (S_ISDIR(fs.st_mode))
249 |                     traverse_dir_thread(q, path);
250 |                 else if (S_ISREG(fs.st_mode)) {
251 |                     enqueue(q, path, size);
252 |                     cnt++;
253 |                 }
254 |             }
255 |         }
256 | #endif
257 |     }
258 |     closedir(dd);
259 |     return cnt;
260 | }
261 | 
262 | void master_func(thread_arg_t* arg) {
263 |     queue_t* q = arg->q;
264 |     int i;
265 |     arg->filenum = traverse_dir_thread(q, arg->dirname);
266 | 
267 |     /* enqueue END_MARK */
268 |     for (i = 0; i < WORKER_THREAD_NUM; ++i) 
269 |         enqueue(q, END_MARK_FNAME, END_MARK_FLENGTH);
270 |     return;
271 | }
272 | 
273 | void worker_func(thread_arg_t* arg) {
274 |     int flen;
275 |     char* fname = NULL;
276 |     queue_t* q = arg->q;
277 |     char* srchstr = arg->srchstr;
278 | 
279 |     int thread_id = arg->id;
280 | 
281 |     printf("worker func %d launched \n", thread_id);
282 |     
283 | #ifdef __CPU_SET
284 |     cpu_set_t mask;    
285 |     __CPU_ZERO(&mask);
286 |     __CPU_SET(arg->cpuid, &mask);
287 |     if (sched_setaffinity(0, sizeof(mask), &mask) == -1)
288 |         printf("WARNING: faild to set CPU affinity...\n");
289 | #endif
290 | 
291 | #if 0
292 |     while (1) {
293 |         int n;
294 | 
295 |         dequeue(q, &fname, &flen));
296 | 
297 |         if (strncmp(fname, END_MARK_FNAME, END_MARK_FLENGTH + 1) == 0)
298 |             break;
299 | 
300 |         n = traverse_file(fname, srchstr, thread_id);
301 |         pthread_mutex_lock(&result.mutex);
302 | 
303 |         if (n > result.num) {
304 |             result.num = n;
305 |             if (result.fname != NULL) free(result.fname);
306 |             result.fname = (char*)malloc(flen);
307 |             strcpy(result.fname, fname);
308 |         }
309 |         pthread_mutex_unlock(&result.mutex);
310 |     }
311 | #else
312 |     char* my_result_fname;
313 |     int my_result_num = 0;
314 |     int my_result_len = 0;
315 |     while (1) {
316 |         int n;
317 | 
318 |         dequeue(q, &fname, &flen);
319 | 
320 |         if (strncmp(fname, END_MARK_FNAME, END_MARK_FLENGTH + 1) == 0)
321 |             break;
322 | 
323 |         n = traverse_file(fname, srchstr, thread_id);
324 | 
325 |         if (n > my_result_num) {
326 |             my_result_num = n;
327 |             my_result_len = flen;
328 |             my_result_fname = (char*)alloca(flen);
329 |             strcpy(my_result_fname, fname);
330 |         }
331 |     }
332 |     pthread_mutex_lock(&result.mutex);
333 |     if (my_result_num > result.num) {
334 |         result.num = my_result_num;
335 | 	if (result.fname != NULL) free(result.fname);
336 |         result.fname = (char*)malloc(my_result_len);
337 |         strcpy(result.fname, my_result_fname);	
338 |     }
339 |     pthread_mutex_unlock(&result.mutex);
340 | #endif
341 |     return;
342 | }
343 | 
344 | void print_result(thread_arg_t* arg) {
345 |     if (result.num) {
346 |         printf("Total %d files\n", arg->filenum);
347 |         printf("Max include file: %s[include %d]\n", result.fname, result.num);	
348 |     }
349 |     return;
350 | }
351 | 
352 | int main(int argc, char* argv[]) {
353 |     int i;
354 |     int thread_num = 1 + WORKER_THREAD_NUM;
355 |     unsigned int t, travdirtime;
356 |     queue_t q;
357 |     thread_arg_t targ[thread_num];
358 |     pthread_t master;
359 |     pthread_t worker[thread_num];
360 |     pthread_t worker2[thread_num];
361 |     int cpu_num;
362 | 
363 |     if (argc != 2) {
364 |         printf("Usage: ./asura [DIR] \n"); return 0;
365 |     }
366 |     cpu_num = sysconf(_SC_NPROCESSORS_CONF);
367 | 
368 |     initqueue(&q);
369 | 
370 |     for (i = 0; i < thread_num; ++i) {
371 |         targ[i].q = &q;
372 |         targ[i].dirname = argv[1];
373 |         targ[i].filenum = 0;
374 |         targ[i].cpuid = i%cpu_num;
375 |     }
376 |     result.fname = NULL;
377 | 
378 |     start_timer(&t);
379 | 
380 |     pthread_mutex_init(&result.mutex, NULL);
381 | 
382 |     pthread_create(&master, NULL, (void*)master_func, (void*)&targ[0]);
383 |     for (i = 1; i < thread_num; ++i)
384 |       { 
385 |         targ[i].id = i;
386 |         pthread_create(&worker[i], NULL, (void*)worker_func, (void*)&targ[i]);
387 |       }
388 |     for (i = 1; i < thread_num; ++i) 
389 |         pthread_join(worker[i], NULL);
390 | 
391 |     travdirtime = stop_timer(&t);
392 |     print_timer(travdirtime);
393 | 
394 |     std::vector<string> s_vec_1;
395 |     std::vector<float> s_vec_2;
396 |     std::vector<string> s_vec_3;
397 |     std::vector<float> s_vec_4;
398 |     
399 |     for( CharTable::iterator i=table.begin(); i!=table.end(); ++i )
400 |     {
401 |       s_vec_1.push_back(i->first);
402 |       s_vec_2.push_back((float)i->second);
403 |       
404 |       cout << "[counts]:" << i->first << "," << i->second << endl;
405 |     }
406 | 
407 |     for( CharTable2::iterator i=table2.begin(); i!=table2.end(); ++i )
408 |     {
409 |       s_vec_3.push_back(i->first);
410 |       s_vec_4.push_back((float)i->second);
411 |       
412 |       cout << "[size]:" << i->first << "," << i->second << endl;
413 |     }
414 | 
415 |     std::remove("reduced");
416 |     ofstream outputfile3("reduced");
417 | 
418 |     std::remove("reduced2");
419 |     ofstream outputfile3_2("reduced2");
420 | 
421 |     cout << "reduced" << endl;
422 |     
423 |     for(int i = 0; i < s_vec_1.size(); i++)
424 |       {
425 | 	/*
426 | 	double tmp_div;
427 | 	if (s_vec_4[i] == 0)
428 | 	  tmp_div = 0;
429 | 	else if (s_vec_4[i] > 0)
430 | 	  tmp_div =  (double)s_vec_2[i] / (double)s_vec_4[i];
431 | 	*/	
432 | 
433 | 	cout << s_vec_1[i] << "," << s_vec_2[i] << "," << s_vec_4[i] << endl;
434 | 	outputfile3 << s_vec_1[i] << "," << s_vec_2[i] << "," << s_vec_4[i] << endl;
435 | 	outputfile3_2 << s_vec_2[i] << "," << s_vec_4[i] << endl;
436 |       }
437 | 
438 |     outputfile3.close();
439 |     outputfile3_2.close();
440 | 
441 |     cout << "strat KMeans " << endl;
442 |     
443 |     const size_t SQRT_K = 2;
444 |     const size_t K = SQRT_K*SQRT_K;
445 |     const size_t M = s_vec_1.size();
446 | 
447 |     cout << "data size:" << M << endl;
448 |     cout << "cluster size:" << K << endl;
449 |     
450 |     point* points;
451 |     points = (struct point *)malloc(M*sizeof(struct point));
452 | 
453 |     point* centroid;
454 |     centroid = (struct point *)malloc(M*sizeof(struct point));
455 | 
456 |     cluster_id* id;
457 |     id = (unsigned short *)malloc(M*sizeof(unsigned short));
458 | 
459 |     std::vector<string> pair; 
460 |     
461 |     for (int i = 0; i < M; i++) {
462 | 
463 |       point& p = points[i];
464 | 
465 |       p.x = (float)s_vec_2[i];
466 |       p.y = (float)s_vec_4[i];
467 |       pair.push_back(s_vec_1[i]);
468 |     }
469 | 
470 |     start_timer(&t);
471 |     tbb_asura::do_k_means( M, points, K, id, centroid );
472 | 
473 |     cout << "KMeans finished" << endl;
474 |     
475 |     /*
476 |     travdirtime = stop_timer(&t);
477 |     print_timer(travdirtime);
478 |     */
479 | 
480 | 
481 |     
482 | #if 1
483 |     int* counts;
484 |     counts = (int *)malloc(K*sizeof(int));
485 | 
486 |     for( size_t i=0; i<M; ++i ) {
487 |       counts[id[i]] = 0;
488 |     }
489 |     
490 |     for( size_t i=0; i<M; ++i ) {
491 |       counts[id[i]]++;
492 |     }
493 | 
494 |     std::remove("clustered");
495 |     ofstream outputfile4("clustered");
496 | 
497 |     for( size_t i=0; i<M; ++i ) {
498 | 
499 |       float percent = (float)counts[id[i]]/(float)M;
500 |       cout << pair[i] << ", clusterID:" << id[i] << ", data(" << points[i].x << "," << points[i].y << ")" << ", " << (int)(percent * 100) << "%" << endl;
501 |       outputfile4 << pair[i] << ", clusterID:" << id[i] << ", data(" << points[i].x << "," << points[i].y << ")" << ", " << (int)(percent * 100) << "%" << endl;
502 | 
503 |       /*
504 | 	    unsigned long pair_long = std::stoul(pair[i]);
505 |       
506 | 	    std::bitset<64> bset_pair = std::bitset<64>(pair_long);   
507 | 	    string bset_pair_string = bset_pair.to_string();
508 | 	    
509 | 	    std::string ip1 = bset_pair_string.substr(0, 8);
510 | 	    std::string ip2 = bset_pair_string.substr(8, 8);
511 | 	    std::string ip3 = bset_pair_string.substr(16, 8);
512 | 	    std::string ip4 = bset_pair_string.substr(24, 8);
513 | 
514 | 	    std::bitset<8> ip1_bset = std::bitset<8>(ip1);
515 | 	    std::bitset<8> ip2_bset = std::bitset<8>(ip2);
516 | 	    std::bitset<8> ip3_bset = std::bitset<8>(ip3);
517 | 	    std::bitset<8> ip4_bset = std::bitset<8>(ip4);
518 | 
519 | 	    outputfile4 << ip1_bset.to_ulong() << "." << ip2_bset.to_ulong() << "." << ip3_bset.to_ulong() << "." << ip4_bset.to_ulong() << ",";
520 | 
521 | 	    std::string ip5 = bset_pair_string.substr(32, 8);
522 | 	    std::string ip6 = bset_pair_string.substr(40, 8);
523 | 	    std::string ip7 = bset_pair_string.substr(48, 8);
524 | 	    std::string ip8 = bset_pair_string.substr(56, 8);
525 | 
526 | 	    std::bitset<8> ip5_bset = std::bitset<8>(ip5);
527 | 	    std::bitset<8> ip6_bset = std::bitset<8>(ip6);
528 | 	    std::bitset<8> ip7_bset = std::bitset<8>(ip7);
529 | 	    std::bitset<8> ip8_bset = std::bitset<8>(ip8);
530 | 
531 | 	    outputfile4 << ip5_bset.to_ulong() << "." << ip6_bset.to_ulong() << "." << ip7_bset.to_ulong() << "." << ip8_bset.to_ulong() << " -> ";
532 | 
533 | 	    float percent = (float)counts[id[i]]/(float)M;
534 | 
535 | 	    outputfile4 << id[i] << " (" << points[i].x << "," << points[i].y << ")," << counts[id[i]] << "," << M << "," << percent << endl;
536 | 	
537 |       */    	
538 |     }
539 | 
540 |     outputfile4.close();
541 | 
542 | #endif
543 |     
544 | #if 0
545 |     printf("centroids = ");
546 |     for( size_t j=0; j<K; ++j ) {
547 |         printf("(%g %g)",centroid[j].x,centroid[j].y);
548 |     }
549 |     printf("\n");
550 | #endif
551 |     
552 |     
553 |     return 0;
554 | }
555 | 
556 | 


--------------------------------------------------------------------------------
/reduced2:
--------------------------------------------------------------------------------
   1 | 4,168
   2 | 118,140493
   3 | 18,760
   4 | 4,168
   5 | 18,760
   6 | 4,168
   7 | 453,630296
   8 | 4,168
   9 | 4,168
  10 | 9407,1.37391e+06
  11 | 4,168
  12 | 4,168
  13 | 12,576
  14 | 4,168
  15 | 6,288
  16 | 29,2822
  17 | 4,168
  18 | 1369,78888
  19 | 4140,182128
  20 | 4,168
  21 | 4068,178960
  22 | 5157,324939
  23 | 4088,179840
  24 | 4163,183140
  25 | 4073,179180
  26 | 871,39123
  27 | 4216,185472
  28 | 4199,184724
  29 | 4028,177200
  30 | 4,168
  31 | 4,168
  32 | 4,168
  33 | 4,168
  34 | 85996,1.46221e+07
  35 | 4,168
  36 | 4,168
  37 | 4,168
  38 | 4,168
  39 | 4,168
  40 | 4,168
  41 | 4,168
  42 | 4,168
  43 | 4,168
  44 | 4,168
  45 | 4,168
  46 | 4,168
  47 | 3830,177052
  48 | 348,19298
  49 | 4,168
  50 | 4,168
  51 | 4,168
  52 | 3188,253089
  53 | 4,168
  54 | 84,14228
  55 | 4,168
  56 | 4,168
  57 | 447,619324
  58 | 4,168
  59 | 4,168
  60 | 4,168
  61 | 4,168
  62 | 4,168
  63 | 4,168
  64 | 4,168
  65 | 4,168
  66 | 4,168
  67 | 4,168
  68 | 4,168
  69 | 194,17206
  70 | 4,168
  71 | 4,168
  72 | 4,168
  73 | 4095,180148
  74 | 4112,180896
  75 | 4,168
  76 | 77544,4.4225e+06
  77 | 4,168
  78 | 4009,176364
  79 | 4437,195196
  80 | 36,5967
  81 | 4296,188992
  82 | 4421,194492
  83 | 4410,194008
  84 | 4,168
  85 | 643,140468
  86 | 4,168
  87 | 4,168
  88 | 486,101318
  89 | 4336,190752
  90 | 4,168
  91 | 1960,498049
  92 | 4585,201708
  93 | 4,168
  94 | 7112,1.4379e+06
  95 | 4,168
  96 | 1430,300119
  97 | 4,168
  98 | 4,168
  99 | 3135,676129
 100 | 4,168
 101 | 4149,182524
 102 | 4,168
 103 | 4069,179004
 104 | 4066,178872
 105 | 4232,186176
 106 | 63,6811
 107 | 4048,178080
 108 | 4116,181072
 109 | 4,168
 110 | 4102,180456
 111 | 4,168
 112 | 4109,180764
 113 | 500,729036
 114 | 4074,179224
 115 | 4069,179004
 116 | 4,168
 117 | 4,168
 118 | 4,168
 119 | 315,17343
 120 | 4,168
 121 | 12,504
 122 | 4,168
 123 | 12,504
 124 | 4,168
 125 | 543,25719
 126 | 4,168
 127 | 4,168
 128 | 944,57860
 129 | 12,504
 130 | 4,168
 131 | 12,504
 132 | 4,168
 133 | 198,183857
 134 | 12,504
 135 | 12,504
 136 | 12,504
 137 | 12,504
 138 | 4810,683934
 139 | 12,504
 140 | 12,504
 141 | 352,497500
 142 | 350,500054
 143 | 216,22476
 144 | 164,10745
 145 | 54,2810
 146 | 18,760
 147 | 6,240
 148 | 18,760
 149 | 18,760
 150 | 18,760
 151 | 65,89328
 152 | 18,760
 153 | 18,760
 154 | 18,760
 155 | 18,760
 156 | 4,168
 157 | 24,1152
 158 | 18,760
 159 | 18,760
 160 | 54312,1.11198e+07
 161 | 11,460
 162 | 4,168
 163 | 11,460
 164 | 4,168
 165 | 5,630
 166 | 4,168
 167 | 12,504
 168 | 4,168
 169 | 4,168
 170 | 12,504
 171 | 4,168
 172 | 12,504
 173 | 4,168
 174 | 12,504
 175 | 4,168
 176 | 12,504
 177 | 4,168
 178 | 11,464
 179 | 4,168
 180 | 1466,73453
 181 | 4,168
 182 | 12,504
 183 | 4,168
 184 | 12,504
 185 | 242,338160
 186 | 4,168
 187 | 1737,2.53597e+06
 188 | 4,168
 189 | 4,168
 190 | 240,338080
 191 | 4,168
 192 | 4,168
 193 | 4,168
 194 | 4,168
 195 | 20796,2.90624e+06
 196 | 492,86934
 197 | 87,5710
 198 | 868,1.28023e+06
 199 | 899,1.28734e+06
 200 | 190,217122
 201 | 872,1.28039e+06
 202 | 4,168
 203 | 4,168
 204 | 4,168
 205 | 53,56981
 206 | 4,168
 207 | 4,168
 208 | 18,760
 209 | 4,168
 210 | 18,760
 211 | 4,168
 212 | 4,168
 213 | 4,168
 214 | 4,168
 215 | 4,168
 216 | 4,168
 217 | 18,760
 218 | 4,168
 219 | 18,760
 220 | 4,168
 221 | 18,760
 222 | 18,760
 223 | 4,168
 224 | 18,760
 225 | 4,168
 226 | 18,760
 227 | 4,168
 228 | 18,760
 229 | 26,2719
 230 | 4,168
 231 | 18,760
 232 | 4,168
 233 | 18,760
 234 | 4,168
 235 | 18,760
 236 | 4,168
 237 | 18,760
 238 | 4,168
 239 | 18,760
 240 | 4,168
 241 | 18,760
 242 | 4,168
 243 | 18,760
 244 | 4,168
 245 | 44,52161
 246 | 18,760
 247 | 4,168
 248 | 18,760
 249 | 4,168
 250 | 4,168
 251 | 144,16278
 252 | 18,760
 253 | 18,760
 254 | 4,168
 255 | 4,168
 256 | 4,168
 257 | 4,168
 258 | 4,168
 259 | 55,6391
 260 | 4,168
 261 | 4,168
 262 | 26604,1.38782e+06
 263 | 4,168
 264 | 4,168
 265 | 85,9002
 266 | 12,504
 267 | 12,504
 268 | 105,109989
 269 | 108,114552
 270 | 1369,506165
 271 | 416,37204
 272 | 12,504
 273 | 11,464
 274 | 204,8168
 275 | 12,504
 276 | 40,4641
 277 | 12,504
 278 | 12,504
 279 | 128,109526
 280 | 12,504
 281 | 12,504
 282 | 12,504
 283 | 4,168
 284 | 4,168
 285 | 4,168
 286 | 4,168
 287 | 6,288
 288 | 38,19594
 289 | 4,168
 290 | 4,168
 291 | 4,168
 292 | 4,168
 293 | 4,168
 294 | 38,19594
 295 | 4,168
 296 | 4,168
 297 | 4,168
 298 | 4,168
 299 | 423,27181
 300 | 4,168
 301 | 4,168
 302 | 49,5021
 303 | 4,168
 304 | 4,168
 305 | 4,168
 306 | 4,168
 307 | 4,168
 308 | 30,1440
 309 | 30,1440
 310 | 4,168
 311 | 4,168
 312 | 158,25384
 313 | 486,53273
 314 | 30,1440
 315 | 4,168
 316 | 4886,654593
 317 | 30,1440
 318 | 4,168
 319 | 30,1440
 320 | 4,168
 321 | 12,504
 322 | 30,1440
 323 | 4,168
 324 | 12,504
 325 | 30,1440
 326 | 4,168
 327 | 12,504
 328 | 30,1440
 329 | 4,168
 330 | 12,504
 331 | 30,1440
 332 | 4,168
 333 | 12,504
 334 | 30,1440
 335 | 4,168
 336 | 12,504
 337 | 4,168
 338 | 4,168
 339 | 4,168
 340 | 4,168
 341 | 4,168
 342 | 4,168
 343 | 4,168
 344 | 4,168
 345 | 4,168
 346 | 4,168
 347 | 4,168
 348 | 4,168
 349 | 4,168
 350 | 4,168
 351 | 4,168
 352 | 4,168
 353 | 4,168
 354 | 4,168
 355 | 304,16280
 356 | 4,168
 357 | 4,168
 358 | 4,168
 359 | 4,168
 360 | 4,168
 361 | 4,168
 362 | 4,168
 363 | 4,168
 364 | 4,168
 365 | 117,14980
 366 | 4,168
 367 | 4,168
 368 | 4,168
 369 | 100,14352
 370 | 18,760
 371 | 4,168
 372 | 18,760
 373 | 18,760
 374 | 18,760
 375 | 18,760
 376 | 18,760
 377 | 18,760
 378 | 18,760
 379 | 4,168
 380 | 384,23830
 381 | 4,168
 382 | 4,168
 383 | 59,70246
 384 | 18,760
 385 | 4,168
 386 | 8,1680
 387 | 18,760
 388 | 4,168
 389 | 4,168
 390 | 4,168
 391 | 4,168
 392 | 59,70247
 393 | 54,12592
 394 | 4,168
 395 | 2039,82424
 396 | 4,168
 397 | 4,168
 398 | 4,168
 399 | 30,1440
 400 | 4,168
 401 | 30,1440
 402 | 4,168
 403 | 909,45973
 404 | 1.23905e+06,7.18492e+07
 405 | 4,168
 406 | 4,168
 407 | 4,168
 408 | 4,168
 409 | 14868,2.22137e+07
 410 | 30,1440
 411 | 4,168
 412 | 30,1440
 413 | 4,168
 414 | 30,1440
 415 | 4,168
 416 | 4,168
 417 | 129,121069
 418 | 30,1440
 419 | 4,168
 420 | 4,168
 421 | 30,1440
 422 | 4,168
 423 | 30,1440
 424 | 4,168
 425 | 5.90949e+06,3.39681e+08
 426 | 30,1440
 427 | 4,168
 428 | 30,1440
 429 | 4,168
 430 | 4,168
 431 | 4127,181556
 432 | 4,168
 433 | 4087,179796
 434 | 4,168
 435 | 4078,179400
 436 | 4,168
 437 | 4135,181908
 438 | 2219,3.21221e+06
 439 | 4,168
 440 | 4144,182304
 441 | 4,168
 442 | 4176,183712
 443 | 4,168
 444 | 4091,179972
 445 | 4,168
 446 | 4113,180940
 447 | 4073,179180
 448 | 4061,178652
 449 | 4,168
 450 | 4,168
 451 | 18,760
 452 | 26616,1.40172e+06
 453 | 18,760
 454 | 33,19996
 455 | 718,72734
 456 | 27,2742
 457 | 18,760
 458 | 18,760
 459 | 18,760
 460 | 1820,98832
 461 | 18,760
 462 | 73,6342
 463 | 18,760
 464 | 18,760
 465 | 18,760
 466 | 6,607
 467 | 18,760
 468 | 4,168
 469 | 4,168
 470 | 4,168
 471 | 4,168
 472 | 4,168
 473 | 4,168
 474 | 4,168
 475 | 4,168
 476 | 26601,1.31683e+06
 477 | 26601,1.31683e+06
 478 | 26601,1.31683e+06
 479 | 26607,1.31712e+06
 480 | 26595,1.31653e+06
 481 | 4,168
 482 | 4,168
 483 | 26610,1.31725e+06
 484 | 3926,172712
 485 | 4040,177728
 486 | 12,504
 487 | 687,48431
 488 | 551,49942
 489 | 12,504
 490 | 99,11348
 491 | 2088,84384
 492 | 48,4441
 493 | 12,504
 494 | 5,566
 495 | 54,55257
 496 | 4090,179928
 497 | 12,504
 498 | 4086,179752
 499 | 12,504
 500 | 4007,176276
 501 | 165,9876
 502 | 12,504
 503 | 4014,176584
 504 | 12,504
 505 | 52,56941
 506 | 3925,172668
 507 | 12,504
 508 | 3960,174208
 509 | 129,121912
 510 | 12,504
 511 | 4044,177904
 512 | 12,504
 513 | 3994,175704
 514 | 30,1440
 515 | 30,1440
 516 | 30,1440
 517 | 30,1440
 518 | 30,1440
 519 | 30,1440
 520 | 30,1440
 521 | 30,1440
 522 | 304,16280
 523 | 30,1440
 524 | 30,1440
 525 | 4,840
 526 | 12,504
 527 | 12,504
 528 | 12,504
 529 | 12,504
 530 | 4,168
 531 | 12,504
 532 | 4,168
 533 | 12,504
 534 | 12,504
 535 | 11,464
 536 | 4,168
 537 | 12,504
 538 | 4,168
 539 | 12,504
 540 | 4,168
 541 | 4,168
 542 | 4,168
 543 | 4,168
 544 | 4,168
 545 | 4,168
 546 | 4,168
 547 | 4,168
 548 | 4,168
 549 | 4,168
 550 | 4,168
 551 | 4,168
 552 | 4,168
 553 | 12,504
 554 | 4,168
 555 | 12,504
 556 | 12,504
 557 | 12,504
 558 | 12,504
 559 | 12,504
 560 | 12,504
 561 | 12,504
 562 | 42,4282
 563 | 135,13409
 564 | 4,168
 565 | 12,504
 566 | 4,168
 567 | 12,504
 568 | 4,168
 569 | 4,168
 570 | 4,168
 571 | 2044,82624
 572 | 4,168
 573 | 4,168
 574 | 4,168
 575 | 4,168
 576 | 4,168
 577 | 4,168
 578 | 56,6189
 579 | 4,168
 580 | 12,504
 581 | 606,33165
 582 | 4283,188420
 583 | 12,504
 584 | 4094,180104
 585 | 30,1440
 586 | 12,504
 587 | 4042,177816
 588 | 30,1440
 589 | 12,504
 590 | 4154,182744
 591 | 12,504
 592 | 3990,175528
 593 | 12,504
 594 | 4066,178872
 595 | 12,504
 596 | 4067,178916
 597 | 12,504
 598 | 4101,180412
 599 | 60,2880
 600 | 179,235040
 601 | 481,680393
 602 | 60,2880
 603 | 232,264143
 604 | 30,1440
 605 | 4071,179092
 606 | 30,1440
 607 | 958,38320
 608 | 260,294707
 609 | 476,681466
 610 | 4105,180588
 611 | 30,1440
 612 | 12,504
 613 | 474,681386
 614 | 30,1440
 615 | 12,504
 616 | 30,1440
 617 | 481,683084
 618 | 30,1440
 619 | 34,12907
 620 | 159,132093
 621 | 1,44
 622 | 4,168
 623 | 4,168
 624 | 4,168
 625 | 4,168
 626 | 4,168
 627 | 47,4401
 628 | 4,168
 629 | 4,168
 630 | 4,168
 631 | 4,168
 632 | 4,168
 633 | 297,15528
 634 | 453,21744
 635 | 453,21744
 636 | 42,7114
 637 | 1245,70621
 638 | 43,51620
 639 | 465,22320
 640 | 465,22320
 641 | 56,6189
 642 | 465,22320
 643 | 1581,2.21561e+06
 644 | 465,22320
 645 | 1549,2.21953e+06
 646 | 465,22320
 647 | 1414,2.06754e+06
 648 | 26590,1.41368e+06
 649 | 465,22320
 650 | 465,22320
 651 | 453,21744
 652 | 3029,377012
 653 | 68,11494
 654 | 4,168
 655 | 3983,175220
 656 | 4,168
 657 | 4033,177420
 658 | 4,168
 659 | 227,315434
 660 | 4193,184460
 661 | 4,168
 662 | 226,314862
 663 | 4092,180016
 664 | 4,168
 665 | 3977,174956
 666 | 4,168
 667 | 4004,176144
 668 | 4,168
 669 | 4039,177684
 670 | 4,168
 671 | 4082,179576
 672 | 4,168
 673 | 52,56941
 674 | 4,168
 675 | 29496,6.85539e+06
 676 | 170,233210
 677 | 1906,97870
 678 | 4111,180852
 679 | 4120,181248
 680 | 702,997554
 681 | 4,168
 682 | 4,168
 683 | 453,21744
 684 | 1682,84203
 685 | 289,18915
 686 | 453,21744
 687 | 453,21744
 688 | 4,168
 689 | 453,21744
 690 | 4,168
 691 | 453,21744
 692 | 4,168
 693 | 453,21744
 694 | 4,168
 695 | 533,25319
 696 | 453,21744
 697 | 4,168
 698 | 453,21744
 699 | 4,168
 700 | 453,21744
 701 | 4,168
 702 | 304,16392
 703 | 453,21744
 704 | 4,168
 705 | 4176,956546
 706 | 542,566525
 707 | 56,13205
 708 | 23,14410
 709 | 12,504
 710 | 60,22144
 711 | 12,504
 712 | 24,960
 713 | 77,43467
 714 | 1,40
 715 | 12,504
 716 | 12,504
 717 | 12,504
 718 | 2639,3.84796e+06
 719 | 12,504
 720 | 12,504
 721 | 12,504
 722 | 317,17423
 723 | 12,504
 724 | 12,504
 725 | 465,22320
 726 | 465,22320
 727 | 465,22320
 728 | 15369,2.22537e+07
 729 | 465,22320
 730 | 465,22320
 731 | 465,22320
 732 | 465,22320
 733 | 465,22320
 734 | 206,219426
 735 | 465,22320
 736 | 465,22320
 737 | 568,30510
 738 | 30,1440
 739 | 30,1440
 740 | 213,13697
 741 | 215,21450
 742 | 30,1440
 743 | 30,1440
 744 | 203,22603
 745 | 30,1440
 746 | 30,1440
 747 | 30,1440
 748 | 30,1440
 749 | 30,1440
 750 | 30,1440
 751 | 218,21570
 752 | 4,168
 753 | 4,168
 754 | 3,144
 755 | 40,4366
 756 | 4,168
 757 | 4,168
 758 | 4,168
 759 | 4,168
 760 | 4,168
 761 | 4,168
 762 | 4,168
 763 | 4,168
 764 | 1853,2.57577e+06
 765 | 3190,4.6198e+06
 766 | 1,40
 767 | 3040,4.38847e+06
 768 | 2848,4.19397e+06
 769 | 5529,325628
 770 | 1525,2.2242e+06
 771 | 2900,4.23202e+06
 772 | 85792,1.45366e+07
 773 | 29736,6.88492e+06
 774 | 30,1440
 775 | 30,1440
 776 | 30,1440
 777 | 30,1440
 778 | 30,1440
 779 | 30,1440
 780 | 30,1440
 781 | 30,1440
 782 | 26617,1.4151e+06
 783 | 30,1440
 784 | 30,1440
 785 | 4,168
 786 | 2214,1.73712e+06
 787 | 4,168
 788 | 482,53658
 789 | 30,1440
 790 | 30,1440
 791 | 30,1440
 792 | 30,1440
 793 | 30,1440
 794 | 465,22320
 795 | 30,1440
 796 | 465,22320
 797 | 30,1440
 798 | 30,1440
 799 | 465,22320
 800 | 465,22320
 801 | 54,55257
 802 | 30,1440
 803 | 465,22320
 804 | 30,1440
 805 | 465,22320
 806 | 465,22320
 807 | 465,22320
 808 | 465,22320
 809 | 352,350158
 810 | 465,22320
 811 | 510,740343
 812 | 5.90957e+06,3.37835e+08
 813 | 511,740383
 814 | 510,740343
 815 | 4,160
 816 | 193,26738
 817 | 59,5897
 818 | 30,1440
 819 | 315,17366
 820 | 30,1440
 821 | 30,1440
 822 | 30,1440
 823 | 36,3096
 824 | 30,1440
 825 | 30,1440
 826 | 30,1440
 827 | 30,1440
 828 | 294,16969
 829 | 30,1440
 830 | 30,1440
 831 | 4648,3.26558e+06
 832 | 1041,77699
 833 | 185,26410
 834 | 4,168
 835 | 1049,1.38881e+06
 836 | 4,168
 837 | 1044,1.39191e+06
 838 | 4,168
 839 | 4,168
 840 | 2163,3.08215e+06
 841 | 4,168
 842 | 274,383637
 843 | 4,168
 844 | 4,168
 845 | 285,400311
 846 | 4,168
 847 | 4,168
 848 | 4,168
 849 | 4,168
 850 | 4,168
 851 | 184,255289
 852 | 4,168
 853 | 4,168
 854 | 195,260068
 855 | 4,168
 856 | 4,168
 857 | 195,260068
 858 | 4,168
 859 | 4,168
 860 | 4,168
 861 | 311,367401
 862 | 85546,1.45409e+07
 863 | 119,9071
 864 | 209,20621
 865 | 22,3299
 866 | 48,4981
 867 | 12,504
 868 | 4,168
 869 | 4,168
 870 | 12,504
 871 | 4,168
 872 | 4,168
 873 | 12,504
 874 | 4,168
 875 | 4,168
 876 | 12,504
 877 | 4,168
 878 | 12,504
 879 | 4,168
 880 | 4,168
 881 | 8,336
 882 | 4,168
 883 | 4,168
 884 | 12,504
 885 | 4,168
 886 | 8,336
 887 | 4,168
 888 | 4,168
 889 | 1,40
 890 | 12,504
 891 | 4,168
 892 | 289,18915
 893 | 12,504
 894 | 4,168
 895 | 2041,82504
 896 | 6,607
 897 | 82,10670
 898 | 79,6582
 899 | 73,9344
 900 | 4,168
 901 | 4,168
 902 | 4,168
 903 | 4,168
 904 | 546,565980
 905 | 4,168
 906 | 4,168
 907 | 97,10002
 908 | 4,168
 909 | 4,168
 910 | 4,168
 911 | 4,168
 912 | 71960,1.94603e+07
 913 | 4,168
 914 | 4,168
 915 | 46788,7.60035e+06
 916 | 4,168
 917 | 4,168
 918 | 72184,1.94766e+07
 919 | 4,168
 920 | 4,168
 921 | 72350,1.95208e+07
 922 | 4,168
 923 | 4,168
 924 | 1305,55131
 925 | 42,7114
 926 | 4,168
 927 | 4,168
 928 | 4,168
 929 | 910,112385
 930 | 4,168
 931 | 3,144
 932 | 4,168
 933 | 4,168
 934 | 4,168
 935 | 12,504
 936 | 4,168
 937 | 12,504
 938 | 4,168
 939 | 12,504
 940 | 4,168
 941 | 12,504
 942 | 4,168
 943 | 73028,4.11658e+06
 944 | 4,168
 945 | 12,504
 946 | 4,168
 947 | 12,504
 948 | 4,168
 949 | 57,3735
 950 | 1812,87860
 951 | 12,504
 952 | 4,168
 953 | 4,168
 954 | 192,19375
 955 | 4,168
 956 | 12,504
 957 | 190,22569
 958 | 12,504
 959 | 1940,2.73366e+06
 960 | 165,191172
 961 | 4,168
 962 | 4,168
 963 | 12,504
 964 | 4,168
 965 | 135,13409
 966 | 12,504
 967 | 4,168
 968 | 558,804560
 969 | 12,504
 970 | 4,168
 971 | 4,168
 972 | 1164,51257
 973 | 500,725816
 974 | 12,504
 975 | 4,168
 976 | 4,168
 977 | 590,814766
 978 | 12,504
 979 | 4,168
 980 | 4,168
 981 | 12,504
 982 | 4,168
 983 | 4,168
 984 | 12,504
 985 | 4,168
 986 | 4,168
 987 | 12,504
 988 | 4,168
 989 | 4,168
 990 | 209,14591
 991 | 500,725816
 992 | 12,504
 993 | 4,168
 994 | 4,168
 995 | 500,725816
 996 | 12,504
 997 | 4,168
 998 | 4,168
 999 | 4,168
1000 | 4,168
1001 | 4,168
1002 | 4,168
1003 | 653,32718
1004 | 795,138203
1005 | 101,108631
1006 | 4,168
1007 | 47,5422
1008 | 4,168
1009 | 4,168
1010 | 4,168
1011 | 4,168
1012 | 4,168
1013 | 4,168
1014 | 4,168
1015 | 12,504
1016 | 4,168
1017 | 4,168
1018 | 11,460
1019 | 4,168
1020 | 4,168
1021 | 12,504
1022 | 4,168
1023 | 4,168
1024 | 102,126956
1025 | 12,504
1026 | 4,168
1027 | 4,168
1028 | 12,504
1029 | 4,168
1030 | 4,168
1031 | 102,126956
1032 | 12,504
1033 | 4,168
1034 | 4,168
1035 | 52,59242
1036 | 12,504
1037 | 4,168
1038 | 4,168
1039 | 43,51620
1040 | 12,504
1041 | 4,168
1042 | 4,168
1043 | 4,168
1044 | 1013,133262
1045 | 4,168
1046 | 12,504
1047 | 4,168
1048 | 12,504
1049 | 4,168
1050 | 4,168
1051 | 4,168
1052 | 5649,436116
1053 | 4,168
1054 | 26604,1.38782e+06
1055 | 4,168
1056 | 2785,176230
1057 | 4,168
1058 | 4,168
1059 | 4,168
1060 | 1687,98894
1061 | 4,168
1062 | 4,168
1063 | 4,168
1064 | 12,504
1065 | 12,504
1066 | 4,168
1067 | 4,168
1068 | 4,168
1069 | 8211,331896
1070 | 4,168
1071 | 4,168
1072 | 206,9048
1073 | 4,168
1074 | 18,760
1075 | 4,168
1076 | 18,760
1077 | 12,504
1078 | 18,760
1079 | 4109,312070
1080 | 12,504
1081 | 18,760
1082 | 12,504
1083 | 18,760
1084 | 102,95125
1085 | 10,420
1086 | 18,760
1087 | 72,88375
1088 | 12,504
1089 | 18,760
1090 | 3890,277099
1091 | 12,504
1092 | 18,760
1093 | 260,25663
1094 | 12,504
1095 | 4,168
1096 | 18,760
1097 | 12,504
1098 | 4,168
1099 | 18,760
1100 | 750,37158
1101 | 26596,1.414e+06
1102 | 649,241055
1103 | 26611,1.38374e+06
1104 | 264,13257
1105 | 4,168
1106 | 4,168
1107 | 59,57552
1108 | 4,168
1109 | 4,168
1110 | 4,168
1111 | 4,168
1112 | 133,8771
1113 | 264,13257
1114 | 4,168
1115 | 324,21049
1116 | 4,168
1117 | 475,477211
1118 | 4,168
1119 | 12,1214
1120 | 4,168
1121 | 12,504
1122 | 4177,183756
1123 | 12,504
1124 | 4114,180984
1125 | 4072,179136
1126 | 12,504
1127 | 4135,181908
1128 | 12,504
1129 | 12,504
1130 | 4104,180544
1131 | 26607,1.31714e+06
1132 | 4043,177860
1133 | 12,504
1134 | 4031,177332
1135 | 12,504
1136 | 12,504
1137 | 4176,183712
1138 | 12,504
1139 | 124,9271
1140 | 12,504
1141 | 825,73045
1142 | 132,7612
1143 | 3993,175660
1144 | 4034,177464
1145 | 3550,264115
1146 | 4,168
1147 | 19,1700
1148 | 4,168
1149 | 4,168
1150 | 4,168
1151 | 4,168
1152 | 4,168
1153 | 4,168
1154 | 4,168
1155 | 4,168
1156 | 4,168
1157 | 4,168
1158 | 4,168
1159 | 4,168
1160 | 4,168
1161 | 4,168
1162 | 4,168
1163 | 56,54482
1164 | 932,47869
1165 | 729,74125
1166 | 4,168
1167 | 4,168
1168 | 4,168
1169 | 4,168
1170 | 127,176937
1171 | 4,168
1172 | 18,760
1173 | 4,168
1174 | 18,760
1175 | 4,168
1176 | 4,168
1177 | 18,760
1178 | 290,367372
1179 | 18,760
1180 | 18,760
1181 | 18,760
1182 | 4,168
1183 | 18,760
1184 | 4,168
1185 | 18,760
1186 | 18,760
1187 | 290,367369
1188 | 18,760
1189 | 255,354326
1190 | 1975,2.93181e+06
1191 | 12,504
1192 | 1414,2.09475e+06
1193 | 12,504
1194 | 14766,3.15398e+06
1195 | 1415,2.09479e+06
1196 | 12,504
1197 | 1,40
1198 | 12,504
1199 | 12,504
1200 | 12,504
1201 | 12,504
1202 | 1462,2.10532e+06
1203 | 12,504
1204 | 1220,77946
1205 | 26,1704
1206 | 4,168
1207 | 261,14920
1208 | 4,168
1209 | 1471,84593
1210 | 4,168
1211 | 12,504
1212 | 29656,6.89764e+06
1213 | 4,168
1214 | 12,504
1215 | 3189,188670
1216 | 4,168
1217 | 1549,214527
1218 | 38,3478
1219 | 506,32845
1220 | 18,760
1221 | 18,760
1222 | 18,760
1223 | 86,103240
1224 | 18,760
1225 | 18,760
1226 | 1086,43440
1227 | 18,760
1228 | 18,760
1229 | 2072,1.20336e+06
1230 | 18,760
1231 | 18,760
1232 | 18,760
1233 | 26599,1.38754e+06
1234 | 4,168
1235 | 4,168
1236 | 301,17198
1237 | 4,168
1238 | 4,168
1239 | 4,168
1240 | 4,168
1241 | 4,168
1242 | 28,2782
1243 | 4,168
1244 | 4,168
1245 | 4,168
1246 | 36,3127
1247 | 4,168
1248 | 2,80
1249 | 4,168
1250 | 4,168
1251 | 4,168
1252 | 4,168
1253 | 4,168
1254 | 4,168
1255 | 4,168
1256 | 4,168
1257 | 4,168
1258 | 4,168
1259 | 4,168
1260 | 4,168
1261 | 1,40
1262 | 4,168
1263 | 4,168
1264 | 4,168
1265 | 4,168
1266 | 4,168
1267 | 4,168
1268 | 4,168
1269 | 4,168
1270 | 18,760
1271 | 18,760
1272 | 18,760
1273 | 4817,6.64982e+06
1274 | 18,760
1275 | 210,13484
1276 | 4,168
1277 | 18,760
1278 | 4,168
1279 | 18,760
1280 | 4,168
1281 | 4,168
1282 | 4,168
1283 | 4,168
1284 | 4,168
1285 | 4,168
1286 | 195,183728
1287 | 4,168
1288 | 4,168
1289 | 4,168
1290 | 1288,61891
1291 | 4,168
1292 | 4,168
1293 | 4,168
1294 | 4,168
1295 | 4,168
1296 | 10,440
1297 | 4,168
1298 | 12,504
1299 | 2005,88212
1300 | 4,168
1301 | 11,460
1302 | 4,168
1303 | 12,504
1304 | 4,168
1305 | 12,504
1306 | 11,464
1307 | 12,504
1308 | 12,504
1309 | 12,504
1310 | 12,504
1311 | 12,504
1312 | 2493,1.80318e+06
1313 | 18,760
1314 | 18,760
1315 | 18,760
1316 | 18,760
1317 | 18,760
1318 | 18,760
1319 | 18,760
1320 | 18,760
1321 | 18,760
1322 | 18,760
1323 | 18,760
1324 | 12,504
1325 | 317,17423
1326 | 12,504
1327 | 18,760
1328 | 18,760
1329 | 12,504
1330 | 1,60
1331 | 18,760
1332 | 11,464
1333 | 4,240
1334 | 870,46752
1335 | 12,504
1336 | 18,760
1337 | 18,760
1338 | 12,504
1339 | 18,760
1340 | 12,504
1341 | 18,760
1342 | 12,504
1343 | 12,504
1344 | 12,504
1345 | 18,760
1346 | 18,760
1347 | 4,168
1348 | 4,168
1349 | 4,168
1350 | 4,168
1351 | 4,168
1352 | 4,168
1353 | 36658,4.79932e+06
1354 | 4,168
1355 | 4,168
1356 | 4,168
1357 | 4,168
1358 | 30,1440
1359 | 4,168
1360 | 30,1440
1361 | 4,168
1362 | 4,168
1363 | 4,168
1364 | 30,1440
1365 | 4,168
1366 | 4,168
1367 | 30,1440
1368 | 4,168
1369 | 4,168
1370 | 30,1440
1371 | 4,168
1372 | 30,1440
1373 | 4,168
1374 | 4,168
1375 | 4,168
1376 | 4,168
1377 | 29492,6.82413e+06
1378 | 4,168
1379 | 418,607569
1380 | 4,168
1381 | 4,168
1382 | 4,168
1383 | 4,168
1384 | 581,62186
1385 | 4,168
1386 | 4,168
1387 | 4,168
1388 | 4,168
1389 | 4,168
1390 | 4,168
1391 | 4,168
1392 | 4,168
1393 | 289,23327
1394 | 4,168
1395 | 4,168
1396 | 4,168
1397 | 4,168
1398 | 4,168
1399 | 4,168
1400 | 4543,6.50118e+06
1401 | 3,144
1402 | 62,5829
1403 | 4,168
1404 | 4,168
1405 | 4,168
1406 | 4,168
1407 | 4,168
1408 | 4,168
1409 | 4,168
1410 | 4,168
1411 | 4,168
1412 | 4,168
1413 | 4,168
1414 | 4,168
1415 | 4089,179884
1416 | 4,168
1417 | 4,168
1418 | 4172,183536
1419 | 4,168
1420 | 4,168
1421 | 4,168
1422 | 226,242201
1423 | 4,168
1424 | 4,168
1425 | 4164,183184
1426 | 4,168
1427 | 4202,184856
1428 | 4,168
1429 | 4304,189344
1430 | 4,168
1431 | 4351,191412
1432 | 4,168
1433 | 4,168
1434 | 4242,186616
1435 | 4,168
1436 | 4,168
1437 | 4197,184636
1438 | 4,168
1439 | 4182,183976
1440 | 4,168
1441 | 4151,182612
1442 | 18,760
1443 | 18,760
1444 | 18,760
1445 | 18,760
1446 | 18,760
1447 | 160,10304
1448 | 18,760
1449 | 76,39188
1450 | 18,760
1451 | 18,760
1452 | 14,592
1453 | 18,760
1454 | 4,168
1455 | 4095,180148
1456 | 4,168
1457 | 4066,178872
1458 | 4,168
1459 | 4056,178432
1460 | 4,168
1461 | 4089,179884
1462 | 4,168
1463 | 1724,79006
1464 | 4032,177376
1465 | 4,168
1466 | 4094,180104
1467 | 4,168
1468 | 4,168
1469 | 4125,181468
1470 | 4,168
1471 | 4,168
1472 | 3986,175352
1473 | 4,168
1474 | 4,168
1475 | 4005,176188
1476 | 4,168
1477 | 4,168
1478 | 4132,181776
1479 | 4,168
1480 | 4,168
1481 | 4,168
1482 | 4,168
1483 | 4,168
1484 | 4,168
1485 | 4,168
1486 | 3946,173592
1487 | 12,504
1488 | 3972,174736
1489 | 12,504
1490 | 4,168
1491 | 4078,179400
1492 | 12,504
1493 | 4,168
1494 | 4039,177684
1495 | 12,504
1496 | 4,168
1497 | 4080,179488
1498 | 124,16544
1499 | 12,504
1500 | 4,168
1501 | 4134,181864
1502 | 12,504
1503 | 4,168
1504 | 193,10718
1505 | 4084,179664
1506 | 12,504
1507 | 4,168
1508 | 4058,178520
1509 | 12,504
1510 | 4,168
1511 | 4027,177156
1512 | 12,504
1513 | 4,168
1514 | 3999,175924
1515 | 12,504
1516 | 4,168
1517 | 4,168
1518 | 4,168
1519 | 4,168
1520 | 30,1440
1521 | 4,168
1522 | 4,168
1523 | 30,1440
1524 | 4,168
1525 | 4,168
1526 | 4,168
1527 | 4,168
1528 | 4,168
1529 | 4,168
1530 | 4,168
1531 | 4,168
1532 | 4,168
1533 | 4,168
1534 | 4,168
1535 | 177,12807
1536 | 30,1440
1537 | 30,1440
1538 | 30,1440
1539 | 4,168
1540 | 4,168
1541 | 785,41754
1542 | 30,1440
1543 | 4,168
1544 | 4,168
1545 | 4,168
1546 | 30,1440
1547 | 4,168
1548 | 71,89568
1549 | 30,1440
1550 | 5210,2.05989e+06
1551 | 311,17649
1552 | 30,1440
1553 | 30,1440
1554 | 4081,179532
1555 | 12,504
1556 | 4025,177068
1557 | 12,504
1558 | 4139,182084
1559 | 12,504
1560 | 1347,82715
1561 | 12,504
1562 | 4103,180500
1563 | 4059,178564
1564 | 12,504
1565 | 4171,183492
1566 | 12,504
1567 | 12,504
1568 | 3993,175660
1569 | 491,36764
1570 | 4084,179664
1571 | 12,504
1572 | 4,840
1573 | 5115,6.73883e+06
1574 | 3999,175924
1575 | 12,504
1576 | 4005,176188
1577 | 12,504
1578 | 30,1440
1579 | 4,168
1580 | 30,1440
1581 | 4,168
1582 | 30,1440
1583 | 4,168
1584 | 30,1440
1585 | 4,168
1586 | 1531,2.22107e+06
1587 | 30,1440
1588 | 30,1440
1589 | 30,1440
1590 | 4,168
1591 | 30,1440
1592 | 4,168
1593 | 30,1440
1594 | 30,1440
1595 | 4,168
1596 | 4,168
1597 | 4,168
1598 | 4,168
1599 | 4,168
1600 | 4,168
1601 | 4,168
1602 | 9993,6.76356e+06
1603 | 4,168
1604 | 736,946189
1605 | 115,6168
1606 | 4,168
1607 | 4,168
1608 | 4,168
1609 | 174,183500
1610 | 4,168
1611 | 4,168
1612 | 4,168
1613 | 4,168
1614 | 30,1440
1615 | 4,168
1616 | 30,1440
1617 | 4,168
1618 | 30,1440
1619 | 3,144
1620 | 30,1440
1621 | 3,144
1622 | 30,1440
1623 | 4,168
1624 | 30,1440
1625 | 4,168
1626 | 500,729036
1627 | 30,1440
1628 | 256,313154
1629 | 30,1440
1630 | 30,1440
1631 | 30,1440
1632 | 28,7331
1633 | 26613,1.38827e+06
1634 | 4,168
1635 | 82,8359
1636 | 4,168
1637 | 165,191172
1638 | 4,168
1639 | 4,168
1640 | 4,168
1641 | 4,168
1642 | 4,168
1643 | 4,168
1644 | 4,168
1645 | 4,168
1646 | 4,168
1647 | 4,168
1648 | 4,168
1649 | 4,168
1650 | 4,168
1651 | 4,168
1652 | 4,168
1653 | 4,168
1654 | 26604,1.38782e+06
1655 | 4,168
1656 | 134,14435
1657 | 4,168
1658 | 453,21744
1659 | 453,21744
1660 | 2,120
1661 | 453,21744
1662 | 1,60
1663 | 453,21744
1664 | 1,60
1665 | 23,2443
1666 | 30,1440
1667 | 1,60
1668 | 30,1440
1669 | 3,180
1670 | 453,21744
1671 | 2,120
1672 | 453,21744
1673 | 225,246519
1674 | 30,1440
1675 | 30,1440
1676 | 170,233210
1677 | 171,233250
1678 | 2648,3.79677e+06
1679 | 4058,178520
1680 | 12,504
1681 | 12,504
1682 | 4015,176628
1683 | 12,504
1684 | 4007,176276
1685 | 4129,181644
1686 | 12,504
1687 | 4203,184900
1688 | 12,504
1689 | 12,504
1690 | 4009,176364
1691 | 12,504
1692 | 4056,178432
1693 | 12,504
1694 | 4188,184240
1695 | 4104,180544
1696 | 12,504
1697 | 4049,178124
1698 | 12,504
1699 | 127,176937
1700 | 163,20564
1701 | 60,5809
1702 | 269,24698
1703 | 4,168
1704 | 4,168
1705 | 4,168
1706 | 4,168
1707 | 4,168
1708 | 4,168
1709 | 4,168
1710 | 4,168
1711 | 4,168
1712 | 4,168
1713 | 465,22320
1714 | 465,22320
1715 | 465,22320
1716 | 465,22320
1717 | 465,22320
1718 | 465,22320
1719 | 465,22320
1720 | 465,22320
1721 | 4142,601439
1722 | 465,22320
1723 | 465,22320
1724 | 146,12397
1725 | 30,1440
1726 | 30,1440
1727 | 30,1440
1728 | 30,1440
1729 | 30,1440
1730 | 30,1440
1731 | 30,1440
1732 | 240,94186
1733 | 30,1440
1734 | 251,276143
1735 | 30,1440
1736 | 30,1440
1737 | 43,51620
1738 | 497,25521
1739 | 6262,9.22066e+06
1740 | 453,21744
1741 | 453,21744
1742 | 453,21744
1743 | 453,21744
1744 | 453,21744
1745 | 453,21744
1746 | 232,27186
1747 | 453,21744
1748 | 453,21744
1749 | 453,21744
1750 | 102,13440
1751 | 453,21744
1752 | 100,13753
1753 | 465,22320
1754 | 465,22320
1755 | 465,22320
1756 | 465,22320
1757 | 200,282832
1758 | 465,22320
1759 | 465,22320
1760 | 400,565664
1761 | 465,22320
1762 | 29522,6.86393e+06
1763 | 465,22320
1764 | 200,282832
1765 | 465,22320
1766 | 465,22320
1767 | 30,1440
1768 | 4,168
1769 | 12,504
1770 | 30,1440
1771 | 4,168
1772 | 12,504
1773 | 30,1440
1774 | 861,38723
1775 | 4,168
1776 | 12,504
1777 | 30,1440
1778 | 958,41840
1779 | 12,504
1780 | 30,1440
1781 | 4,168
1782 | 12,504
1783 | 30,1440
1784 | 4,168
1785 | 12,504
1786 | 30,1440
1787 | 4,168
1788 | 12,504
1789 | 30,1440
1790 | 192,13911
1791 | 4,168
1792 | 12,504
1793 | 12,504
1794 | 12,504
1795 | 30,1440
1796 | 30,1440
1797 | 75,63484
1798 | 48,2522
1799 | 3,120
1800 | 386,303933
1801 | 943,47876
1802 | 835,1.08907e+06
1803 | 1086,47460
1804 | 30,1440
1805 | 30,1440
1806 | 43,51620
1807 | 30,1440
1808 | 1136,53701
1809 | 30,1440
1810 | 30,1440
1811 | 1553,2.25495e+06
1812 | 26595,1.40952e+06
1813 | 9385,487247
1814 | 30,1440
1815 | 30,1440
1816 | 30,1440
1817 | 30,1440
1818 | 30,1440
1819 | 2295,308226
1820 | 447,619324
1821 | 26592,1.31643e+06
1822 | 26592,1.31644e+06
1823 | 293,15368
1824 | 26586,1.31614e+06
1825 | 26613,1.31748e+06
1826 | 87,5710
1827 | 7002,537259
1828 | 30,1440
1829 | 30,1440
1830 | 30,1440
1831 | 465,22320
1832 | 77188,1.21898e+07
1833 | 30,1440
1834 | 465,22320
1835 | 30,1440
1836 | 72575,4.20078e+06
1837 | 30,1440
1838 | 30,1440
1839 | 2003,80128
1840 | 30,1440
1841 | 30,1440
1842 | 465,22320
1843 | 30,1440
1844 | 465,22320
1845 | 465,22320
1846 | 465,22320
1847 | 465,22320
1848 | 465,22320
1849 | 465,22320
1850 | 465,22320
1851 | 721,55107
1852 | 4,168
1853 | 253,14784
1854 | 440,55325
1855 | 121,7172
1856 | 319,40861
1857 | 576,66725
1858 | 610,55827
1859 | 104,11551
1860 | 185,19698
1861 | 4,168
1862 | 390,27407
1863 | 4,168
1864 | 1030,1.35056e+06
1865 | 996,1.33378e+06
1866 | 4,168
1867 | 1118,1.37797e+06
1868 | 4,168
1869 | 1001,1.349e+06
1870 | 916,1.31493e+06
1871 | 12,504
1872 | 12,504
1873 | 12,504
1874 | 119,115207
1875 | 12,504
1876 | 12,504
1877 | 12,504
1878 | 465,22320
1879 | 12,504
1880 | 79,101588
1881 | 465,22320
1882 | 12,504
1883 | 465,22320
1884 | 12,504
1885 | 465,22320
1886 | 574,775425
1887 | 12,504
1888 | 465,22320
1889 | 1269,1.85024e+06
1890 | 465,22320
1891 | 510,24480
1892 | 284,404708
1893 | 510,24480
1894 | 510,24480
1895 | 5560,274416
1896 | 3717,5.55415e+06
1897 | 22838,3.68029e+06
1898 | 3717,5.55415e+06
1899 | 23070,3.71692e+06
1900 | 906,51816
1901 | 23270,3.71748e+06
1902 | 3722,5.55435e+06
1903 | 23224,3.7227e+06
1904 | 23106,3.6983e+06
1905 | 3713,5.55253e+06
1906 | 23122,3.70954e+06
1907 | 4,168
1908 | 12,504
1909 | 4,168
1910 | 12,504
1911 | 18,760
1912 | 4,168
1913 | 12,504
1914 | 18,760
1915 | 4,168
1916 | 12,504
1917 | 4,168
1918 | 12,504
1919 | 165,9876
1920 | 4,168
1921 | 12,504
1922 | 4,168
1923 | 12,504
1924 | 4,168
1925 | 12,504
1926 | 18,760
1927 | 136,8891
1928 | 18,760
1929 | 18,760
1930 | 4,168
1931 | 18,13364
1932 | 18,760
1933 | 4,168
1934 | 12,504
1935 | 18,760
1936 | 12,504
1937 | 18,760
1938 | 2,80
1939 | 111,108548
1940 | 18,760
1941 | 18,760
1942 | 418,607569
1943 | 1.23904e+06,7.16688e+07
1944 | 671,793890
1945 | 118,9031
1946 | 12,504
1947 | 12,504
1948 | 5,630
1949 | 4,168
1950 | 12,504
1951 | 5,566
1952 | 4,168
1953 | 12,504
1954 | 4,168
1955 | 12,504
1956 | 4,168
1957 | 12,504
1958 | 129,153694
1959 | 4,168
1960 | 12,504
1961 | 4,168
1962 | 12,504
1963 | 151,166172
1964 | 4,168
1965 | 12,504
1966 | 4,168
1967 | 12,504
1968 | 79,8602
1969 | 44,53014
1970 | 971,71618
1971 | 681,45528
1972 | 36662,1.81746e+07
1973 | 29182,6.79439e+06
1974 | 11,460
1975 | 12,504
1976 | 12,504
1977 | 3,144
1978 | 12,504
1979 | 11,460
1980 | 239,269301
1981 | 12,504
1982 | 1247,257205
1983 | 12,504
1984 | 12,504
1985 | 12,504
1986 | 12,504
1987 | 246,23809
1988 | 4,168
1989 | 700,772162
1990 | 


--------------------------------------------------------------------------------