Calculating and Comparing Product Production Time for VS One-Health Products
\n",
"Sheldon Waugh MSc, PhD \n",
"Epidemiologist \n",
"One-Health Division \n",
"Veterinary Services and Public Health Sanitation Directorate \n",
"US Army Public Health Center \n",
"\n",
"### This Python notebook is a workable and scalable example on how to make a notebook that cleans and analyzes data in a transparent way.\n",
"We have an excel workbook that details the progress of all One-Health/Veterinary Services products to include posters, brochures and newsletters. \n",
"Taking a look at the initial structure of this spreadsheet, we have missing dates or progress for certain products. This may be due to certain products not needing certain levels of approval. It may be advantageous for us to split the data by product and then look at the time related to each product. \n",
"For my sanity, I will be getting rid of the incomplete records."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Analysis"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"First let's go ahead and load the .csv into python"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
\n",
"\n",
"
\n",
" \n",
"
\n",
"
\n",
"
Type
\n",
"
Draft Review
\n",
"
SME Content Review
\n",
"
PMD Content Review
\n",
"
VID Design Work/Edits
\n",
"
Submitted to CRC
\n",
"
CRC Review Complete
\n",
"
\n",
" \n",
" \n",
"
\n",
"
0
\n",
"
Brochure
\n",
"
NaN
\n",
"
NaN
\n",
"
5/17/2016
\n",
"
5/23/2016
\n",
"
6/1/2016
\n",
"
6/7/2016
\n",
"
\n",
"
\n",
"
1
\n",
"
Brochure
\n",
"
NaN
\n",
"
NaN
\n",
"
5/17/2016
\n",
"
5/23/2016
\n",
"
6/3/2016
\n",
"
6/13/2016
\n",
"
\n",
"
\n",
"
2
\n",
"
Brochure
\n",
"
NaN
\n",
"
NaN
\n",
"
5/17/2016
\n",
"
10/20/2016
\n",
"
11/1/2016
\n",
"
12/20/2016
\n",
"
\n",
"
\n",
"
3
\n",
"
Brochure
\n",
"
NaN
\n",
"
NaN
\n",
"
7/5/2016
\n",
"
7/5/2016
\n",
"
8/19/2016
\n",
"
9/2/2016
\n",
"
\n",
"
\n",
"
4
\n",
"
Brochure
\n",
"
5/9/2017
\n",
"
5/16/2017
\n",
"
6/6/2017
\n",
"
6/23/2017
\n",
"
8/4/2017
\n",
"
10/10/2017
\n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" Type Draft Review SME Content Review PMD Content Review \\\n",
"0 Brochure NaN NaN 5/17/2016 \n",
"1 Brochure NaN NaN 5/17/2016 \n",
"2 Brochure NaN NaN 5/17/2016 \n",
"3 Brochure NaN NaN 7/5/2016 \n",
"4 Brochure 5/9/2017 5/16/2017 6/6/2017 \n",
"\n",
" VID Design Work/Edits Submitted to CRC CRC Review Complete \n",
"0 5/23/2016 6/1/2016 6/7/2016 \n",
"1 5/23/2016 6/3/2016 6/13/2016 \n",
"2 10/20/2016 11/1/2016 12/20/2016 \n",
"3 7/5/2016 8/19/2016 9/2/2016 \n",
"4 6/23/2017 8/4/2017 10/10/2017 "
]
},
"execution_count": 1,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"import pandas as pd # import the module and alias it as pd\n",
"\n",
"timedateVOH_data = pd.read_csv('C:\\Users\\waugh\\Dropbox\\Documents\\One_Health\\VS_Product_Analysis_dates.csv', parse_dates=True)\n",
"timedateVOH_data.head()"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
\n",
"\n",
"
\n",
" \n",
"
\n",
"
\n",
"
Type
\n",
"
Draft Review
\n",
"
SME Content Review
\n",
"
PMD Content Review
\n",
"
VID Design Work/Edits
\n",
"
Submitted to CRC
\n",
"
CRC Review Complete
\n",
"
\n",
" \n",
" \n",
"
\n",
"
count
\n",
"
81
\n",
"
52
\n",
"
48
\n",
"
62
\n",
"
70
\n",
"
75
\n",
"
69
\n",
"
\n",
"
\n",
"
unique
\n",
"
3
\n",
"
36
\n",
"
29
\n",
"
31
\n",
"
33
\n",
"
33
\n",
"
21
\n",
"
\n",
"
\n",
"
top
\n",
"
Newsletter
\n",
"
6/1/2017
\n",
"
4/20/2016
\n",
"
7/11/2016
\n",
"
4/20/2016
\n",
"
8/29/2017
\n",
"
8/15/2016
\n",
"
\n",
"
\n",
"
freq
\n",
"
41
\n",
"
4
\n",
"
5
\n",
"
8
\n",
"
6
\n",
"
8
\n",
"
9
\n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" Type Draft Review SME Content Review PMD Content Review \\\n",
"count 81 52 48 62 \n",
"unique 3 36 29 31 \n",
"top Newsletter 6/1/2017 4/20/2016 7/11/2016 \n",
"freq 41 4 5 8 \n",
"\n",
" VID Design Work/Edits Submitted to CRC CRC Review Complete \n",
"count 70 75 69 \n",
"unique 33 33 21 \n",
"top 4/20/2016 8/29/2017 8/15/2016 \n",
"freq 6 8 9 "
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"timedateVOH_data.describe()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"We see that the dates we have on the csv file did not read into python as a datetime object. Without this, we're unable to determine the elapsed time for each column.\n",
"\n",
"Lets fix for the CRC columns at least so since the product wouldn't exist without CRC approval."
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"text/html": [
"
\n",
"\n",
"
\n",
" \n",
"
\n",
"
\n",
"
Type
\n",
"
Draft Review
\n",
"
SME Content Review
\n",
"
PMD Content Review
\n",
"
VID Design Work/Edits
\n",
"
Submitted to CRC
\n",
"
CRC Review Complete
\n",
"
CRCtime
\n",
"
\n",
" \n",
" \n",
"
\n",
"
0
\n",
"
Brochure
\n",
"
NaN
\n",
"
NaN
\n",
"
5/17/2016
\n",
"
5/23/2016
\n",
"
2016-06-01
\n",
"
2016-06-07
\n",
"
6 days
\n",
"
\n",
"
\n",
"
1
\n",
"
Brochure
\n",
"
NaN
\n",
"
NaN
\n",
"
5/17/2016
\n",
"
5/23/2016
\n",
"
2016-06-03
\n",
"
2016-06-13
\n",
"
10 days
\n",
"
\n",
"
\n",
"
2
\n",
"
Brochure
\n",
"
NaN
\n",
"
NaN
\n",
"
5/17/2016
\n",
"
10/20/2016
\n",
"
2016-11-01
\n",
"
2016-12-20
\n",
"
49 days
\n",
"
\n",
"
\n",
"
3
\n",
"
Brochure
\n",
"
NaN
\n",
"
NaN
\n",
"
7/5/2016
\n",
"
7/5/2016
\n",
"
2016-08-19
\n",
"
2016-09-02
\n",
"
14 days
\n",
"
\n",
"
\n",
"
4
\n",
"
Brochure
\n",
"
5/9/2017
\n",
"
5/16/2017
\n",
"
6/6/2017
\n",
"
6/23/2017
\n",
"
2017-08-04
\n",
"
2017-10-10
\n",
"
67 days
\n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" Type Draft Review SME Content Review PMD Content Review \\\n",
"0 Brochure NaN NaN 5/17/2016 \n",
"1 Brochure NaN NaN 5/17/2016 \n",
"2 Brochure NaN NaN 5/17/2016 \n",
"3 Brochure NaN NaN 7/5/2016 \n",
"4 Brochure 5/9/2017 5/16/2017 6/6/2017 \n",
"\n",
" VID Design Work/Edits Submitted to CRC CRC Review Complete CRCtime \n",
"0 5/23/2016 2016-06-01 2016-06-07 6 days \n",
"1 5/23/2016 2016-06-03 2016-06-13 10 days \n",
"2 10/20/2016 2016-11-01 2016-12-20 49 days \n",
"3 7/5/2016 2016-08-19 2016-09-02 14 days \n",
"4 6/23/2017 2017-08-04 2017-10-10 67 days "
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"df_clean = timedateVOH_data[pd.notnull(timedateVOH_data['Submitted to CRC '])]\n",
"df_clean_1 = df_clean[pd.notnull(df_clean['CRC Review Complete'])]\n",
"df_clean_1\n",
"df_clean_2 = df_clean_1.copy()\n",
"\n",
"df_clean_2['CRC Review Complete'] = pd.to_datetime(df_clean_2['CRC Review Complete'], format = '%m/%d/%Y')\n",
"df_clean_2['Submitted to CRC '] = pd.to_datetime(df_clean_2['Submitted to CRC '], format = '%m/%d/%Y')\n",
"\n",
"df_clean_2['CRCtime'] = df_clean_2['CRC Review Complete'] - df_clean_2['Submitted to CRC ']\n",
"\n",
"df_clean_2.head()\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Now that we finished the column, we now have the calculated the elapsed time!\n",
"We see, however, that certain products don't have full processing times through all columns. \n",
"This probably means that certain products may not need certain approvals. \n",
"So from now on, let's go ahead and split the created data.frame by the two main One-Health Products (Newsletters and Brochures)."
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"text/html": [
"
\n",
"\n",
"
\n",
" \n",
"
\n",
"
\n",
"
Type
\n",
"
Draft Review
\n",
"
SME Content Review
\n",
"
PMD Content Review
\n",
"
VID Design Work/Edits
\n",
"
Submitted to CRC
\n",
"
CRC Review Complete
\n",
"
CRCtime
\n",
"
VID_Edittime
\n",
"
PMD_Reviewtime
\n",
"
SME_Reviewtime
\n",
"
Draft_Reviewtime
\n",
"
\n",
" \n",
" \n",
"
\n",
"
40
\n",
"
Newsletter
\n",
"
2016-04-13
\n",
"
2016-04-20
\n",
"
2016-04-20
\n",
"
2016-04-20
\n",
"
2016-04-28
\n",
"
2016-05-09
\n",
"
11 days
\n",
"
8 days
\n",
"
0 days
\n",
"
0 days
\n",
"
7 days
\n",
"
\n",
"
\n",
"
41
\n",
"
Newsletter
\n",
"
2016-04-15
\n",
"
2016-04-20
\n",
"
2016-04-20
\n",
"
2016-04-20
\n",
"
2016-04-28
\n",
"
2016-05-09
\n",
"
11 days
\n",
"
8 days
\n",
"
0 days
\n",
"
0 days
\n",
"
5 days
\n",
"
\n",
"
\n",
"
42
\n",
"
Newsletter
\n",
"
2016-04-15
\n",
"
2016-04-20
\n",
"
2016-04-20
\n",
"
2016-04-20
\n",
"
2016-04-28
\n",
"
2016-05-09
\n",
"
11 days
\n",
"
8 days
\n",
"
0 days
\n",
"
0 days
\n",
"
5 days
\n",
"
\n",
"
\n",
"
43
\n",
"
Newsletter
\n",
"
2016-04-15
\n",
"
2016-04-20
\n",
"
2016-04-20
\n",
"
2016-04-20
\n",
"
2016-04-28
\n",
"
2016-05-09
\n",
"
11 days
\n",
"
8 days
\n",
"
0 days
\n",
"
0 days
\n",
"
5 days
\n",
"
\n",
"
\n",
"
48
\n",
"
Newsletter
\n",
"
2016-05-17
\n",
"
2016-05-17
\n",
"
2016-05-17
\n",
"
2016-05-17
\n",
"
2016-06-03
\n",
"
2016-06-07
\n",
"
4 days
\n",
"
17 days
\n",
"
0 days
\n",
"
0 days
\n",
"
0 days
\n",
"
\n",
"
\n",
"
54
\n",
"
Newsletter
\n",
"
2016-07-01
\n",
"
2016-07-05
\n",
"
2016-07-11
\n",
"
2016-07-11
\n",
"
2016-07-22
\n",
"
2016-08-15
\n",
"
24 days
\n",
"
11 days
\n",
"
0 days
\n",
"
6 days
\n",
"
4 days
\n",
"
\n",
"
\n",
"
56
\n",
"
Newsletter
\n",
"
2016-06-27
\n",
"
2016-06-27
\n",
"
2016-07-11
\n",
"
2016-07-11
\n",
"
2016-07-22
\n",
"
2016-08-15
\n",
"
24 days
\n",
"
11 days
\n",
"
0 days
\n",
"
14 days
\n",
"
0 days
\n",
"
\n",
"
\n",
"
58
\n",
"
Newsletter
\n",
"
2016-11-17
\n",
"
2016-11-23
\n",
"
2016-11-28
\n",
"
2016-11-29
\n",
"
2017-01-27
\n",
"
2017-03-30
\n",
"
62 days
\n",
"
59 days
\n",
"
1 days
\n",
"
5 days
\n",
"
6 days
\n",
"
\n",
"
\n",
"
59
\n",
"
Newsletter
\n",
"
2016-11-08
\n",
"
2016-11-23
\n",
"
2016-11-28
\n",
"
2016-11-29
\n",
"
2017-01-27
\n",
"
2017-03-30
\n",
"
62 days
\n",
"
59 days
\n",
"
1 days
\n",
"
5 days
\n",
"
15 days
\n",
"
\n",
"
\n",
"
61
\n",
"
Newsletter
\n",
"
2017-02-25
\n",
"
2017-03-15
\n",
"
2017-03-15
\n",
"
2017-03-17
\n",
"
2017-03-17
\n",
"
2017-05-08
\n",
"
52 days
\n",
"
0 days
\n",
"
2 days
\n",
"
0 days
\n",
"
18 days
\n",
"
\n",
"
\n",
"
62
\n",
"
Newsletter
\n",
"
2017-02-10
\n",
"
2017-03-13
\n",
"
2017-03-15
\n",
"
2017-03-20
\n",
"
2017-04-04
\n",
"
2017-06-01
\n",
"
58 days
\n",
"
15 days
\n",
"
5 days
\n",
"
2 days
\n",
"
31 days
\n",
"
\n",
"
\n",
"
63
\n",
"
Newsletter
\n",
"
2017-02-10
\n",
"
2017-03-10
\n",
"
2017-03-13
\n",
"
2017-03-15
\n",
"
2017-04-04
\n",
"
2017-06-01
\n",
"
58 days
\n",
"
20 days
\n",
"
2 days
\n",
"
3 days
\n",
"
28 days
\n",
"
\n",
"
\n",
"
64
\n",
"
Newsletter
\n",
"
2017-03-01
\n",
"
2017-03-16
\n",
"
2017-03-16
\n",
"
2017-03-20
\n",
"
2017-04-04
\n",
"
2017-06-01
\n",
"
58 days
\n",
"
15 days
\n",
"
4 days
\n",
"
0 days
\n",
"
15 days
\n",
"
\n",
"
\n",
"
66
\n",
"
Newsletter
\n",
"
2017-06-01
\n",
"
2017-06-23
\n",
"
2017-06-29
\n",
"
2017-07-25
\n",
"
2017-08-29
\n",
"
2017-08-29
\n",
"
0 days
\n",
"
35 days
\n",
"
26 days
\n",
"
6 days
\n",
"
22 days
\n",
"
\n",
"
\n",
"
68
\n",
"
Newsletter
\n",
"
2017-06-01
\n",
"
2017-06-23
\n",
"
2017-06-28
\n",
"
2017-07-25
\n",
"
2017-08-29
\n",
"
2017-08-29
\n",
"
0 days
\n",
"
35 days
\n",
"
27 days
\n",
"
5 days
\n",
"
22 days
\n",
"
\n",
"
\n",
"
70
\n",
"
Newsletter
\n",
"
2016-07-01
\n",
"
2016-07-05
\n",
"
2016-07-11
\n",
"
2017-08-28
\n",
"
2017-08-29
\n",
"
2017-08-29
\n",
"
0 days
\n",
"
1 days
\n",
"
413 days
\n",
"
6 days
\n",
"
4 days
\n",
"
\n",
"
\n",
"
71
\n",
"
Newsletter
\n",
"
2016-06-07
\n",
"
2016-06-30
\n",
"
2016-07-11
\n",
"
2017-08-28
\n",
"
2017-08-29
\n",
"
2017-08-29
\n",
"
0 days
\n",
"
1 days
\n",
"
413 days
\n",
"
11 days
\n",
"
23 days
\n",
"
\n",
"
\n",
"
72
\n",
"
Newsletter
\n",
"
2016-06-27
\n",
"
2016-06-27
\n",
"
2016-07-11
\n",
"
2017-08-28
\n",
"
2017-08-29
\n",
"
2017-08-29
\n",
"
0 days
\n",
"
1 days
\n",
"
413 days
\n",
"
14 days
\n",
"
0 days
\n",
"
\n",
"
\n",
"
77
\n",
"
Newsletter
\n",
"
2018-08-15
\n",
"
2018-08-24
\n",
"
2018-08-24
\n",
"
2018-08-29
\n",
"
2018-08-30
\n",
"
2018-09-04
\n",
"
5 days
\n",
"
1 days
\n",
"
5 days
\n",
"
0 days
\n",
"
9 days
\n",
"
\n",
"
\n",
"
78
\n",
"
Newsletter
\n",
"
2018-08-19
\n",
"
2018-08-30
\n",
"
2018-08-29
\n",
"
2018-08-29
\n",
"
2018-08-30
\n",
"
2018-09-04
\n",
"
5 days
\n",
"
1 days
\n",
"
0 days
\n",
"
-1 days
\n",
"
11 days
\n",
"
\n",
"
\n",
"
79
\n",
"
Newsletter
\n",
"
2018-08-15
\n",
"
2018-08-24
\n",
"
2018-08-24
\n",
"
2018-08-29
\n",
"
2018-08-30
\n",
"
2018-09-04
\n",
"
5 days
\n",
"
1 days
\n",
"
5 days
\n",
"
0 days
\n",
"
9 days
\n",
"
\n",
"
\n",
"
80
\n",
"
Newsletter
\n",
"
2018-08-16
\n",
"
2018-08-23
\n",
"
2018-08-24
\n",
"
2018-08-29
\n",
"
2018-08-30
\n",
"
2018-09-04
\n",
"
5 days
\n",
"
1 days
\n",
"
5 days
\n",
"
1 days
\n",
"
7 days
\n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" Type Draft Review SME Content Review PMD Content Review \\\n",
"40 Newsletter 2016-04-13 2016-04-20 2016-04-20 \n",
"41 Newsletter 2016-04-15 2016-04-20 2016-04-20 \n",
"42 Newsletter 2016-04-15 2016-04-20 2016-04-20 \n",
"43 Newsletter 2016-04-15 2016-04-20 2016-04-20 \n",
"48 Newsletter 2016-05-17 2016-05-17 2016-05-17 \n",
"54 Newsletter 2016-07-01 2016-07-05 2016-07-11 \n",
"56 Newsletter 2016-06-27 2016-06-27 2016-07-11 \n",
"58 Newsletter 2016-11-17 2016-11-23 2016-11-28 \n",
"59 Newsletter 2016-11-08 2016-11-23 2016-11-28 \n",
"61 Newsletter 2017-02-25 2017-03-15 2017-03-15 \n",
"62 Newsletter 2017-02-10 2017-03-13 2017-03-15 \n",
"63 Newsletter 2017-02-10 2017-03-10 2017-03-13 \n",
"64 Newsletter 2017-03-01 2017-03-16 2017-03-16 \n",
"66 Newsletter 2017-06-01 2017-06-23 2017-06-29 \n",
"68 Newsletter 2017-06-01 2017-06-23 2017-06-28 \n",
"70 Newsletter 2016-07-01 2016-07-05 2016-07-11 \n",
"71 Newsletter 2016-06-07 2016-06-30 2016-07-11 \n",
"72 Newsletter 2016-06-27 2016-06-27 2016-07-11 \n",
"77 Newsletter 2018-08-15 2018-08-24 2018-08-24 \n",
"78 Newsletter 2018-08-19 2018-08-30 2018-08-29 \n",
"79 Newsletter 2018-08-15 2018-08-24 2018-08-24 \n",
"80 Newsletter 2018-08-16 2018-08-23 2018-08-24 \n",
"\n",
" VID Design Work/Edits Submitted to CRC CRC Review Complete CRCtime \\\n",
"40 2016-04-20 2016-04-28 2016-05-09 11 days \n",
"41 2016-04-20 2016-04-28 2016-05-09 11 days \n",
"42 2016-04-20 2016-04-28 2016-05-09 11 days \n",
"43 2016-04-20 2016-04-28 2016-05-09 11 days \n",
"48 2016-05-17 2016-06-03 2016-06-07 4 days \n",
"54 2016-07-11 2016-07-22 2016-08-15 24 days \n",
"56 2016-07-11 2016-07-22 2016-08-15 24 days \n",
"58 2016-11-29 2017-01-27 2017-03-30 62 days \n",
"59 2016-11-29 2017-01-27 2017-03-30 62 days \n",
"61 2017-03-17 2017-03-17 2017-05-08 52 days \n",
"62 2017-03-20 2017-04-04 2017-06-01 58 days \n",
"63 2017-03-15 2017-04-04 2017-06-01 58 days \n",
"64 2017-03-20 2017-04-04 2017-06-01 58 days \n",
"66 2017-07-25 2017-08-29 2017-08-29 0 days \n",
"68 2017-07-25 2017-08-29 2017-08-29 0 days \n",
"70 2017-08-28 2017-08-29 2017-08-29 0 days \n",
"71 2017-08-28 2017-08-29 2017-08-29 0 days \n",
"72 2017-08-28 2017-08-29 2017-08-29 0 days \n",
"77 2018-08-29 2018-08-30 2018-09-04 5 days \n",
"78 2018-08-29 2018-08-30 2018-09-04 5 days \n",
"79 2018-08-29 2018-08-30 2018-09-04 5 days \n",
"80 2018-08-29 2018-08-30 2018-09-04 5 days \n",
"\n",
" VID_Edittime PMD_Reviewtime SME_Reviewtime Draft_Reviewtime \n",
"40 8 days 0 days 0 days 7 days \n",
"41 8 days 0 days 0 days 5 days \n",
"42 8 days 0 days 0 days 5 days \n",
"43 8 days 0 days 0 days 5 days \n",
"48 17 days 0 days 0 days 0 days \n",
"54 11 days 0 days 6 days 4 days \n",
"56 11 days 0 days 14 days 0 days \n",
"58 59 days 1 days 5 days 6 days \n",
"59 59 days 1 days 5 days 15 days \n",
"61 0 days 2 days 0 days 18 days \n",
"62 15 days 5 days 2 days 31 days \n",
"63 20 days 2 days 3 days 28 days \n",
"64 15 days 4 days 0 days 15 days \n",
"66 35 days 26 days 6 days 22 days \n",
"68 35 days 27 days 5 days 22 days \n",
"70 1 days 413 days 6 days 4 days \n",
"71 1 days 413 days 11 days 23 days \n",
"72 1 days 413 days 14 days 0 days \n",
"77 1 days 5 days 0 days 9 days \n",
"78 1 days 0 days -1 days 11 days \n",
"79 1 days 5 days 0 days 9 days \n",
"80 1 days 5 days 1 days 7 days "
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"#Create the dataframe subset\n",
"df_clean_2_News = df_clean_2.loc[df_clean_2['Type'] == 'Newsletter']\n",
"df_clean_2_News.head()\n",
"\n",
"#Erase rows with NaN values\n",
"df_clean_3_News = df_clean_2_News[pd.notnull(df_clean_2_News['VID Design Work/Edits'])]\n",
"df_clean_4_News = df_clean_3_News[pd.notnull(df_clean_3_News['PMD Content Review'])]\n",
"df_clean_5_News = df_clean_4_News[pd.notnull(df_clean_4_News['SME Content Review'])]\n",
"df_clean_6_News = df_clean_5_News[pd.notnull(df_clean_5_News['Draft Review'])]\n",
"\n",
"#Cenvert strings to datetime objects\n",
"df_clean_6_News['VID Design Work/Edits'] = pd.to_datetime(df_clean_6_News['VID Design Work/Edits'], format = '%m/%d/%Y')\n",
"df_clean_6_News['PMD Content Review'] = pd.to_datetime(df_clean_6_News['PMD Content Review'], format = '%m/%d/%Y')\n",
"df_clean_6_News['SME Content Review'] = pd.to_datetime(df_clean_6_News['SME Content Review'], format = '%m/%d/%Y')\n",
"df_clean_6_News['Draft Review'] = pd.to_datetime(df_clean_6_News['Draft Review'], format = '%m/%d/%Y')\n",
"\n",
"df_clean_6_News.head()\n",
"\n",
"#Create the elapsed time columns\n",
"df_clean_6_News['VID_Edittime'] = df_clean_6_News['Submitted to CRC '] - df_clean_6_News['VID Design Work/Edits']\n",
"df_clean_6_News['PMD_Reviewtime'] = df_clean_6_News['VID Design Work/Edits'] - df_clean_6_News['PMD Content Review']\n",
"df_clean_6_News['SME_Reviewtime'] = df_clean_6_News['PMD Content Review'] - df_clean_6_News['SME Content Review']\n",
"df_clean_6_News['Draft_Reviewtime'] = df_clean_6_News['SME Content Review'] - df_clean_6_News['Draft Review']\n",
"\n",
"df_clean_7_News = df_clean_6_News.drop(df_clean_6_News[df_clean_6_News.PMD_Reviewtime == '-10 days'].index)\n",
"df_clean_8_News = df_clean_7_News.drop(df_clean_7_News[df_clean_7_News.PMD_Reviewtime == '-6 days'].index)\n",
"#df_clean_6_News['PMD_Reviewtime'].total_seconds()\n",
"\n",
"#Good to go!\n",
"df_clean_8_News"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Now lets take a look at the summary data....."
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
\n",
"\n",
"
\n",
" \n",
"
\n",
"
\n",
"
CRCtime
\n",
"
VID_Edittime
\n",
"
PMD_Reviewtime
\n",
"
SME_Reviewtime
\n",
"
Draft_Reviewtime
\n",
"
\n",
" \n",
" \n",
"
\n",
"
count
\n",
"
22
\n",
"
22
\n",
"
22
\n",
"
22
\n",
"
22
\n",
"
\n",
"
\n",
"
mean
\n",
"
21 days 04:21:49.090909
\n",
"
14 days 08:43:38.181818
\n",
"
60 days 02:10:54.545454
\n",
"
3 days 12:00:00
\n",
"
11 days 04:21:49.090909
\n",
"
\n",
"
\n",
"
std
\n",
"
24 days 06:22:33.479019
\n",
"
17 days 14:47:37.238031
\n",
"
143 days 17:24:49.094471
\n",
"
4 days 14:02:38.929146
\n",
"
9 days 05:32:01.438379
\n",
"
\n",
"
\n",
"
min
\n",
"
0 days 00:00:00
\n",
"
0 days 00:00:00
\n",
"
0 days 00:00:00
\n",
"
-1 days +00:00:00
\n",
"
0 days 00:00:00
\n",
"
\n",
"
\n",
"
25%
\n",
"
4 days 06:00:00
\n",
"
1 days 00:00:00
\n",
"
0 days 00:00:00
\n",
"
0 days 00:00:00
\n",
"
5 days 00:00:00
\n",
"
\n",
"
\n",
"
50%
\n",
"
11 days 00:00:00
\n",
"
8 days 00:00:00
\n",
"
2 days 00:00:00
\n",
"
1 days 12:00:00
\n",
"
8 days 00:00:00
\n",
"
\n",
"
\n",
"
75%
\n",
"
45 days 00:00:00
\n",
"
16 days 12:00:00
\n",
"
5 days 00:00:00
\n",
"
5 days 18:00:00
\n",
"
17 days 06:00:00
\n",
"
\n",
"
\n",
"
max
\n",
"
62 days 00:00:00
\n",
"
59 days 00:00:00
\n",
"
413 days 00:00:00
\n",
"
14 days 00:00:00
\n",
"
31 days 00:00:00
\n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" CRCtime VID_Edittime \\\n",
"count 22 22 \n",
"mean 21 days 04:21:49.090909 14 days 08:43:38.181818 \n",
"std 24 days 06:22:33.479019 17 days 14:47:37.238031 \n",
"min 0 days 00:00:00 0 days 00:00:00 \n",
"25% 4 days 06:00:00 1 days 00:00:00 \n",
"50% 11 days 00:00:00 8 days 00:00:00 \n",
"75% 45 days 00:00:00 16 days 12:00:00 \n",
"max 62 days 00:00:00 59 days 00:00:00 \n",
"\n",
" PMD_Reviewtime SME_Reviewtime \\\n",
"count 22 22 \n",
"mean 60 days 02:10:54.545454 3 days 12:00:00 \n",
"std 143 days 17:24:49.094471 4 days 14:02:38.929146 \n",
"min 0 days 00:00:00 -1 days +00:00:00 \n",
"25% 0 days 00:00:00 0 days 00:00:00 \n",
"50% 2 days 00:00:00 1 days 12:00:00 \n",
"75% 5 days 00:00:00 5 days 18:00:00 \n",
"max 413 days 00:00:00 14 days 00:00:00 \n",
"\n",
" Draft_Reviewtime \n",
"count 22 \n",
"mean 11 days 04:21:49.090909 \n",
"std 9 days 05:32:01.438379 \n",
"min 0 days 00:00:00 \n",
"25% 5 days 00:00:00 \n",
"50% 8 days 00:00:00 \n",
"75% 17 days 06:00:00 \n",
"max 31 days 00:00:00 "
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"df_clean_8_News.describe()"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"11 days 00:00:00\n",
"8 days 00:00:00\n",
"2 days 00:00:00\n",
"1 days 12:00:00\n",
"8 days 00:00:00\n"
]
}
],
"source": [
"print(df_clean_8_News.CRCtime.median())\n",
"print(df_clean_8_News.VID_Edittime.median())\n",
"print(df_clean_8_News.PMD_Reviewtime.median())\n",
"print(df_clean_8_News.SME_Reviewtime.median())\n",
"print(df_clean_8_News.Draft_Reviewtime.median())"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Yikes, we have some weird data that led to some negative elapsed times. \n",
"We could do one of three things:\n",
"- Delete the records (rows) with the offending dates.\n",
"- Delete the PMD_Reviewtime column and analyze the other criteria.\n",
"- Leave the data unedited and ask MAJ Watkins about the dates.\n",
"Let's go ahead and do the other products.\n",
"\n",
"This type of discrepancy also demonstrates the lack of an effective way of recording data. In order to effectively produce an analysis and report of any substance, our data must be reliable. \n",
"\n",
"A recommendation for the future is that an accurate working document MUST be established in order to produce reliable production time data and reports. \n",
"\n",
"For this analysis, I've decided to delete the offending products with a negative elapsed time."
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
\n",
"\n",
"
\n",
" \n",
"
\n",
"
\n",
"
Type
\n",
"
Draft Review
\n",
"
SME Content Review
\n",
"
PMD Content Review
\n",
"
VID Design Work/Edits
\n",
"
Submitted to CRC
\n",
"
CRC Review Complete
\n",
"
CRCtime
\n",
"
VID_Edittime
\n",
"
PMD_Reviewtime
\n",
"
SME_Reviewtime
\n",
"
Draft_Reviewtime
\n",
"
\n",
" \n",
" \n",
"
\n",
"
4
\n",
"
Brochure
\n",
"
2017-05-09
\n",
"
2017-05-16
\n",
"
2017-06-06
\n",
"
2017-06-23
\n",
"
2017-08-04
\n",
"
2017-10-10
\n",
"
67 days
\n",
"
42 days
\n",
"
17 days
\n",
"
21 days
\n",
"
7 days
\n",
"
\n",
"
\n",
"
7
\n",
"
Brochure
\n",
"
2016-11-23
\n",
"
2016-12-01
\n",
"
2016-12-02
\n",
"
2016-12-06
\n",
"
2017-01-30
\n",
"
2017-03-30
\n",
"
59 days
\n",
"
55 days
\n",
"
4 days
\n",
"
1 days
\n",
"
8 days
\n",
"
\n",
"
\n",
"
9
\n",
"
Brochure
\n",
"
2016-12-07
\n",
"
2017-02-09
\n",
"
2017-02-10
\n",
"
2017-05-17
\n",
"
2017-06-06
\n",
"
2017-06-30
\n",
"
24 days
\n",
"
20 days
\n",
"
96 days
\n",
"
1 days
\n",
"
64 days
\n",
"
\n",
"
\n",
"
10
\n",
"
Brochure
\n",
"
2016-12-09
\n",
"
2017-01-26
\n",
"
2017-01-27
\n",
"
2017-01-30
\n",
"
2017-02-02
\n",
"
2017-03-30
\n",
"
56 days
\n",
"
3 days
\n",
"
3 days
\n",
"
1 days
\n",
"
48 days
\n",
"
\n",
"
\n",
"
11
\n",
"
Brochure
\n",
"
2016-09-06
\n",
"
2016-09-09
\n",
"
2016-09-09
\n",
"
2016-09-10
\n",
"
2016-11-29
\n",
"
2017-01-02
\n",
"
34 days
\n",
"
80 days
\n",
"
1 days
\n",
"
0 days
\n",
"
3 days
\n",
"
\n",
"
\n",
"
13
\n",
"
Brochure
\n",
"
2017-03-01
\n",
"
2017-03-10
\n",
"
2017-03-13
\n",
"
2017-04-10
\n",
"
2017-06-05
\n",
"
2017-06-27
\n",
"
22 days
\n",
"
56 days
\n",
"
28 days
\n",
"
3 days
\n",
"
9 days
\n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" Type Draft Review SME Content Review PMD Content Review \\\n",
"4 Brochure 2017-05-09 2017-05-16 2017-06-06 \n",
"7 Brochure 2016-11-23 2016-12-01 2016-12-02 \n",
"9 Brochure 2016-12-07 2017-02-09 2017-02-10 \n",
"10 Brochure 2016-12-09 2017-01-26 2017-01-27 \n",
"11 Brochure 2016-09-06 2016-09-09 2016-09-09 \n",
"13 Brochure 2017-03-01 2017-03-10 2017-03-13 \n",
"\n",
" VID Design Work/Edits Submitted to CRC CRC Review Complete CRCtime \\\n",
"4 2017-06-23 2017-08-04 2017-10-10 67 days \n",
"7 2016-12-06 2017-01-30 2017-03-30 59 days \n",
"9 2017-05-17 2017-06-06 2017-06-30 24 days \n",
"10 2017-01-30 2017-02-02 2017-03-30 56 days \n",
"11 2016-09-10 2016-11-29 2017-01-02 34 days \n",
"13 2017-04-10 2017-06-05 2017-06-27 22 days \n",
"\n",
" VID_Edittime PMD_Reviewtime SME_Reviewtime Draft_Reviewtime \n",
"4 42 days 17 days 21 days 7 days \n",
"7 55 days 4 days 1 days 8 days \n",
"9 20 days 96 days 1 days 64 days \n",
"10 3 days 3 days 1 days 48 days \n",
"11 80 days 1 days 0 days 3 days \n",
"13 56 days 28 days 3 days 9 days "
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"#Create the dataframe subset\n",
"df_clean_2_1 = df_clean_2.copy()\n",
"df_clean_2_bro = df_clean_2_1.loc[df_clean_2['Type'] == 'Brochure']\n",
"df_clean_2_bro.head()\n",
"pd.options.mode.chained_assignment = None\n",
"\n",
"#Erase rows with NaN values\n",
"df_clean_3_bro = df_clean_2_bro[pd.notnull(df_clean_2_bro['VID Design Work/Edits'])]\n",
"df_clean_4_bro = df_clean_3_bro[pd.notnull(df_clean_3_bro['PMD Content Review'])]\n",
"df_clean_5_bro = df_clean_4_bro[pd.notnull(df_clean_4_bro['SME Content Review'])]\n",
"df_clean_6_bro = df_clean_5_bro[pd.notnull(df_clean_5_bro['Draft Review'])]\n",
"\n",
"#Cenvert strings to datetime objects\n",
"df_clean_6_bro['VID Design Work/Edits'] = pd.to_datetime(df_clean_6_bro['VID Design Work/Edits'], format = '%m/%d/%Y')\n",
"df_clean_6_bro['PMD Content Review'] = pd.to_datetime(df_clean_6_bro['PMD Content Review'], format = '%m/%d/%Y')\n",
"df_clean_6_bro['SME Content Review'] = pd.to_datetime(df_clean_6_bro['SME Content Review'], format = '%m/%d/%Y')\n",
"df_clean_6_bro['Draft Review'] = pd.to_datetime(df_clean_6_bro['Draft Review'], format = '%m/%d/%Y')\n",
"\n",
"df_clean_6_News.head()\n",
"\n",
"#Create the elapsed time columns\n",
"df_clean_6_bro['VID_Edittime'] = df_clean_6_bro['Submitted to CRC '] - df_clean_6_bro['VID Design Work/Edits']\n",
"df_clean_6_bro['PMD_Reviewtime'] = df_clean_6_bro['VID Design Work/Edits'] - df_clean_6_bro['PMD Content Review']\n",
"df_clean_6_bro['SME_Reviewtime'] = df_clean_6_bro['PMD Content Review'] - df_clean_6_bro['SME Content Review']\n",
"df_clean_6_bro['Draft_Reviewtime'] = df_clean_6_bro['SME Content Review'] - df_clean_6_bro['Draft Review']\n",
"\n",
"#Good to go!\n",
"df_clean_6_bro"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
\n",
"\n",
"
\n",
" \n",
"
\n",
"
\n",
"
CRCtime
\n",
"
VID_Edittime
\n",
"
PMD_Reviewtime
\n",
"
SME_Reviewtime
\n",
"
Draft_Reviewtime
\n",
"
\n",
" \n",
" \n",
"
\n",
"
count
\n",
"
6
\n",
"
6
\n",
"
6
\n",
"
6
\n",
"
6
\n",
"
\n",
"
\n",
"
mean
\n",
"
43 days 16:00:00
\n",
"
42 days 16:00:00
\n",
"
24 days 20:00:00
\n",
"
4 days 12:00:00
\n",
"
23 days 04:00:00
\n",
"
\n",
"
\n",
"
std
\n",
"
19 days 09:32:32.152061
\n",
"
27 days 14:37:12.911653
\n",
"
36 days 08:56:33.395270
\n",
"
8 days 03:25:10.375190
\n",
"
26 days 00:15:41.341229
\n",
"
\n",
"
\n",
"
min
\n",
"
22 days 00:00:00
\n",
"
3 days 00:00:00
\n",
"
1 days 00:00:00
\n",
"
0 days 00:00:00
\n",
"
3 days 00:00:00
\n",
"
\n",
"
\n",
"
25%
\n",
"
26 days 12:00:00
\n",
"
25 days 12:00:00
\n",
"
3 days 06:00:00
\n",
"
1 days 00:00:00
\n",
"
7 days 06:00:00
\n",
"
\n",
"
\n",
"
50%
\n",
"
45 days 00:00:00
\n",
"
48 days 12:00:00
\n",
"
10 days 12:00:00
\n",
"
1 days 00:00:00
\n",
"
8 days 12:00:00
\n",
"
\n",
"
\n",
"
75%
\n",
"
58 days 06:00:00
\n",
"
55 days 18:00:00
\n",
"
25 days 06:00:00
\n",
"
2 days 12:00:00
\n",
"
38 days 06:00:00
\n",
"
\n",
"
\n",
"
max
\n",
"
67 days 00:00:00
\n",
"
80 days 00:00:00
\n",
"
96 days 00:00:00
\n",
"
21 days 00:00:00
\n",
"
64 days 00:00:00
\n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" CRCtime VID_Edittime \\\n",
"count 6 6 \n",
"mean 43 days 16:00:00 42 days 16:00:00 \n",
"std 19 days 09:32:32.152061 27 days 14:37:12.911653 \n",
"min 22 days 00:00:00 3 days 00:00:00 \n",
"25% 26 days 12:00:00 25 days 12:00:00 \n",
"50% 45 days 00:00:00 48 days 12:00:00 \n",
"75% 58 days 06:00:00 55 days 18:00:00 \n",
"max 67 days 00:00:00 80 days 00:00:00 \n",
"\n",
" PMD_Reviewtime SME_Reviewtime \\\n",
"count 6 6 \n",
"mean 24 days 20:00:00 4 days 12:00:00 \n",
"std 36 days 08:56:33.395270 8 days 03:25:10.375190 \n",
"min 1 days 00:00:00 0 days 00:00:00 \n",
"25% 3 days 06:00:00 1 days 00:00:00 \n",
"50% 10 days 12:00:00 1 days 00:00:00 \n",
"75% 25 days 06:00:00 2 days 12:00:00 \n",
"max 96 days 00:00:00 21 days 00:00:00 \n",
"\n",
" Draft_Reviewtime \n",
"count 6 \n",
"mean 23 days 04:00:00 \n",
"std 26 days 00:15:41.341229 \n",
"min 3 days 00:00:00 \n",
"25% 7 days 06:00:00 \n",
"50% 8 days 12:00:00 \n",
"75% 38 days 06:00:00 \n",
"max 64 days 00:00:00 "
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"df_clean_6_bro.describe()"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"45 days 00:00:00\n",
"48 days 12:00:00\n",
"10 days 12:00:00\n",
"1 days 00:00:00\n",
"8 days 12:00:00\n"
]
}
],
"source": [
"print(df_clean_6_bro.CRCtime.median())\n",
"print(df_clean_6_bro.VID_Edittime.median())\n",
"print(df_clean_6_bro.PMD_Reviewtime.median())\n",
"print(df_clean_6_bro.SME_Reviewtime.median())\n",
"print(df_clean_6_bro.Draft_Reviewtime.median())"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"
Plots and Figures
\n",
"Now lets go ahead and produce some figures to visualize our data and demonstrate which production stage takes the most and least amount of time. \n",
"First let's take a look at the brochure products..."
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjYAAAGxCAYAAABx6/zIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzs3Xt8z/X///H7285jlm2MZZicmZxKVkIO5VAJUSSH9SklNZFQOXRwqk8qnxzqm1OECpUUDSEhmZRDJDkzh2HD2Gx7/f54/fbmbQfv92ze83K7Xi6vy97v5+v1fr0f79fer/fue76er9fbZhiGIQAAAAso4u4CAAAA8gvBBgAAWAbBBgAAWAbBBgAAWAbBBgAAWAbBBgAAWAbBBgAAWAbBBgAAWAbBBgAAWAbBpgA88sgj8vPz0+nTp3Ncplu3bvLy8tLRo0clSQkJCRoyZIhq1KihokWLKjAwUNWqVVP37t31559/OvW8P//8s3x8fLRv3z6H9osXL2rSpElq1KiRAgMD5efnp+rVq2vw4MFKSEjI+wvNBxUqVFC7du2ynbdx40bZbDZNnz69wGvo2bOn/f7hw4c1YsQIbd68OcuyPXv2VLFixa66zlOnTumWW27R119/7VQNK1eulM1ms08eHh4KDQ3Vo48+qr/++svp13ItMmtYuXJlgaw/OTlZI0aMcHr9e/fuddgmRYoUUXBwsNq0aaN169YVSI051VCQ78FRo0Zl+z4p6N9Hdq58H+Y2SdL06dNls9m0d+/e61ajM06fPq2QkBDNnTvX3jZixIgs76cyZcqoTZs2+uWXX9xSZ+b7691333XL8+fk4sWLuu222/T++++7u5Q8IdgUgOjoaF24cEGff/55tvMTExO1cOFCtWvXTqGhoTp79qzuuusuTZ8+XU899ZS+/fZbzZ49W08//bT27NmT7R/YKxmGoZiYGP3nP/9R+fLl7e3Jyclq2bKl+vXrp7p162rOnDn6/vvv1b17d3388ceqW7eudu7cmW+v3QoOHz6skSNHOrXdc1KiRAn1799fL7/8slJTU51+3KhRo7Ru3Tr99NNPeuWVVxQbG6u7775bhw4dynMthUVycrJGjhzp8h/qfv36ad26dfr55581evRo/fHHH2rWrJl+//33gin0Ossp2NSrV0/r1q1TvXr1rlstmc95+VS6dGndfffdWdolqW3btlq3bp3KlClz3Wp0xsiRIxUWFqYuXbpkmbdkyRKtW7dOa9as0fjx4xUfH6+mTZtq06ZNbqi0cPLy8tKwYcP0xhtvuP2f3zwxkO/S0tKMsLAwo379+tnOnzRpkiHJWLRokWEYhjF16lRDkrFixYpsl09PT7/qc37//feGJGPHjh0O7U8//bQhyZg7d26Wx+zcudMIDAw0atasaaSlpV31OQpC+fLljbZt22Y777fffjMkGdOmTSvwGnr06OHU8/bo0cMoWrSoU+uNj483PD09jdmzZ1912Z9++smQZHz55ZcO7Z9++qkhyXjrrbdyfGxqaqpx8eJFp2pypoaffvrpmteVnePHjxuSjOHDhzu1/J49ewxJxjvvvOPQvnz5ckOS8dRTT+X42LS0NOPChQvXUq5DDQX5HixatKjD+6+wyW0fLYwSEhIMPz8/Y/LkyQ7tw4cPNyQZx48fd2jfvXu3IckYMmRIruvNr/3scjm9x/NbXvaHlJQUIygoyHj77bcLqKqCQ49NAfDw8FCPHj0UFxenLVu2ZJk/bdo0lSlTRq1bt5YkeyLO6b+eIkWu/muaNGmS7rjjDlWtWtXeFh8fr6lTp+r+++/P9j+XKlWq6JVXXtG2bdsc/mPMPDy0ZMkS1atXT35+fqpWrZqmTp2aZR3x8fF65plnVLZsWXl7eysiIkIjR45UWlraVWvOq127dqlr164qVaqUfHx8VL16dX300UcOy1y4cEEDBgxQnTp1FBgYqKCgIDVq1EjffPNNruteuXKl7rjjDklSr1697N3WI0aMcFjun3/+UZs2bVSsWDGFh4drwIABSklJcVgmNDRULVu21OTJk/P8Wu+66y5Jsh9ezDxU8Nlnn2nAgAG69dZb5ePjo3/++UeStHXrVj388MMqUaKEfH19VadOHc2YMSPLenfs2KEHHnhA/v7+CgkJUZ8+fXTmzJksy115mC5T06ZN1bRpU4e206dPa8CAAapYsaJ8fHxUqlQptWnTRjt27NDevXtVsmRJSeZ/05nbNbt1u7pNMrvzx40bp7feeksRERHy8fHRTz/9JEnav3+/nnjiCYf3y3//+19lZGQ4rPfw4cPq3LmzAgICFBgYqC5duig+Pt6p1y6ZhykrVKjg0JaSkqI33nhD1atXl6+vr4KDg9WsWTOtXbtWkmSz2XTu3DnNmDHDvk0y153Toahvv/1WjRo1kr+/vwICAtSyZcssh+YyD7ts27ZNjz/+uAIDAxUaGqrevXsrMTHRqe3sjOwORTVt2lS1atXSunXrFBUVJT8/P1WoUEHTpk2TJC1evFj16tWTv7+/IiMjtWTJkizrdWYfz62mtLS0bD/zshMYGCjJ7KXIlF/7WW77xJXee+89RUREqFixYmrUqJHWr1/vMN/Z993V9oekpCQNHDhQERER8vb21q233qqYmBidO3fOYb3e3t7q0qWLPv74Yxk32Hdle7q7AKvq3bu3xowZo6lTp2r8+PH29u3bt2vDhg0aPHiwPDw8JEmNGjWSJD355JMaOnSoGjdurODgYKefKzU1VcuWLVO/fv0c2n/66SelpaWpffv2OT62ffv2Gjp0qGJjY9WxY0d7+x9//KEBAwZo8ODBCg0N1f/93/8pOjpalSpV0r333ivJDDV33nmnihQpomHDhum2227TunXr9NZbb2nv3r32D7KrMQwj2yCUnp6epW379u2KiopSuXLl9N///lelS5fW0qVL9cILL+jEiRMaPny4JPMPysmTJzVw4EDdeuut9m3UoUMHTZs2TU8++WS2tdSrV0/Tpk1Tr1699Nprr6lt27aSpLJly9qXuXjxoh566CFFR0drwIABWr16td58800FBgZq2LBhDutr2rSphgwZotOnT+uWW25xantcLvODNDMUZBoyZIgaNWqkyZMnq0iRIipVqpR27typqKgolSpVSh9++KGCg4M1a9Ys9ezZU0ePHtWgQYMkSUePHlWTJk3k5eWliRMnKjQ0VLNnz9bzzz/vcn2Zzpw5o3vuuUd79+7VK6+8ooYNG+rs2bNavXq1jhw5oqioKC1ZskQPPPCAoqOj9dRTT2X7uq5lm3z44YeqUqWK3n33XRUvXlyVK1fW8ePHFRUVpdTUVL355puqUKGCvvvuOw0cOFC7d+/WxIkTJUnnz59XixYtdPjwYY0ePVpVqlTR4sWLnf7jmJ20tDS1bt1aP//8s2JiYnTfffcpLS1N69ev1/79+xUVFaV169bpvvvuU7NmzfT6669LkooXL57jOj///HN169ZNrVq10pw5c5SSkqJx48apadOmWr58ue655x6H5Tt27KguXbooOjpaW7Zs0ZAhQyQp239S8lN8fLx69eqlQYMGqWzZspowYYJ69+6tAwcO6KuvvtLQoUMVGBioN954Q+3bt9e///6rsLAwSc7v4zlZvHix6tatm+P+lp6errS0NGVkZGj//v167bXX5OPjo06dOmVZ9lr2s6vtE9WqVbM/z0cffaRq1arZx7S8/vrratOmjfbs2WMPXq7Kbn9ITk5WkyZNdPDgQQ0dOlS1a9fWtm3bNGzYMG3ZskXLli2zj5+SzM+vSZMmaevWrYqMjMxTHW7h7i4jK2vSpIkREhJipKam2tsGDBhgSDL+/vtvh2XfeOMNw9vb25BkSDIiIiKMPn36GH/88cdVn+fXX3/N9nDTmDFjDEnGkiVLcnzs+fPnDUlG69at7W3ly5c3fH19jX379jksFxQUZDzzzDP2tmeeecYoVqyYw3KGYRjvvvuuIcnYtm3bVWsvX768/TXnNF1+GOD+++83ypYtayQmJjqs5/nnnzd8fX2NkydPZvs8aWlpxsWLF43o6Gijbt26WWpw5VCUJOOLL75waG/Tpo1RtWrVLMvHxsYakowffvgh1+2QeRho3rx5xsWLF43k5GRj9erVRqVKlQwPDw/7+yBzuXvvvTfLOh577DHDx8fH2L9/v0N769atDX9/f+P06dOGYRjGK6+8YthsNmPz5s0Oy7Vs2TLLoagrt02mJk2aGE2aNLHff+ONNwxJRmxsbI6vMa+HosaOHWtcvHjRuHDhghEXF2fccccdhiRj8eLFDsvddtttDvuaYRjG4MGDDUnGr7/+6tD+7LPPGjabzdi5c6dhGJcOD3/zzTcOy/3nP//J8l648rVn6tGjh1G+fHn7/ZkzZxqSjE8++STX15nToagrDw2mp6cbYWFhRmRkpMPh6TNnzhilSpUyoqKi7G2Zh13GjRvnsM7nnnvO8PX1NTIyMnKt6XK5HYqaNm2aIcnYs2ePva1JkyaGJGPjxo32toSEBMPDw8Pw8/MzDh06ZG/fvHmzIcn48MMP7W153ccz+fv7G3369MnSnrlNrpyKFy9uLFiwwGHZ/NjPnNknMt+7kZGRDsMBNmzYYEgy5syZY29z9n2X2/4wevRoo0iRIsZvv/3m0P7VV18Zkozvv//eoX3Xrl2GJGPSpEk5vobCiENRBSg6OlonTpzQt99+K8n8D27WrFlq3LixKleu7LDs66+/rv3792vq1Kl65plnVKxYMU2ePFn169fXnDlzcn2ew4cPS5JKlSqV51ovT+mSVKdOHZUrV85+39fXV1WqVHE44+q7775Ts2bNFBYWprS0NPuUeYht1apVki79h5Q5XXkI4J577tFvv/2WZZo5c6bDchcuXNDy5cv1yCOPyN/f32Gdbdq00YULFxy6b7/88kvdfffdKlasmDw9PeXl5aVPP/30ms8ystlsevDBBx3aateuneVsNOnS78TZwb9dunSRl5eX/P39de+99yo9PV1fffWVateu7bDc5b1rmVasWKHmzZsrPDzcob1nz55KTk62H6746aefVLNmTd1+++0Oy3Xt2tWpGrPzww8/qEqVKmrRokWe15GTV155RV5eXvL19VX9+vW1f/9+TZkyRW3atHFY7qGHHnI4nCCZ26RGjRq68847Hdp79uwpwzC0YsUKSeY2CQgI0EMPPeSw3LVuE19fX/Xu3TvP67jczp07dfjwYXXv3t3h8HSxYsXUsWNHrV+/XsnJyQ6PufL11K5dWxcuXNCxY8fypaaclClTRvXr17ffDwoKUqlSpVSnTh17z4wkVa9eXdKlw4qu7uNXOn36tJKTk3P9LFy2bJl+++03bdiwQd99951atGihxx57TAsXLsyy7LXsZ67sE23btrX34Euy7+/ZfaY4K7v94bvvvlOtWrVUp04dh217//33Z3vY09XPr8KCQ1EFqFOnTurXr5+mTZumjh076vvvv9fRo0c1duzYbJcPDQ1Vr1691KtXL0nS6tWr1bp1a7344ot6/PHHc3ye8+fPSzLDx+Uyg8mePXtyfGzmvCt30uwOhfn4+NifSzIPaSxatCjLzpPpxIkTkqTmzZvbQ44k9ejRw+H02cDAQDVo0CDHGjMlJCQoLS1NEyZM0IQJE3J9zgULFqhz58569NFH9fLLL6t06dLy9PTUpEmTrrkb3t/fP8u29vHx0YULF7Ism7nc5dstN2PHjtV9990nDw8PhYSEZPm9ZMpuPFZCQkK27Zl/SDLHciUkJCgiIiLLcqVLl3aqxuwcP37cIQjnpxdffFFPPPGEihQpoltuuUURERFZgriU8za5ctyLlP02CQ0NzbLctW6TsLAwp8bIOSO3sXhhYWHKyMjQqVOn5O/vb2+/cj/28fGR5Pz7Ma+CgoKytHl7e2dp9/b2liT7vuPKPp6dnD4LL3f77bcrJCTEfr9169aKjIxU37599cgjjzgsey37mSv7REH8nrKr8ejRo/rnn3+u+pmdydXPr8KCYFOA/Pz89Pjjj+uTTz7RkSNHNHXqVAUEBOjRRx916vH33nuvWrVqpa+//lrHjh3L8b+QzJ305MmTDu3NmjWTp6envv76a/Xp0yfbx2YOGm7ZsqWzL8vheWvXrq2333472/mZO/qUKVMcBqZe/qHiihIlSsjDw0Pdu3dX3759s10m8w/2rFmzFBERoXnz5jn8EbxygG9By/ydOPuaK1as6FTIy+4Pe3BwsI4cOZKlPbNHL7OG4ODgbAfFZtfm6+ub7TY7ceKEw2sqWbKkDh48eNW686Js2bLXZZts2LAhy3I5bZPsBuBe+UehZMmSWrNmjTIyMvIl3GT+8cvp9RQpUkQlSpS45udxJ1f28exkbqMrPwtzU6RIEdWsWVNffvllls/Za3lP5fc+4ez7LlN2tYeEhMjPzy/Hf+6u/Jxy9fOrsOBQVAGLjo5Wenq63nnnHX3//fd67LHHHP6jkswUfeXhGck8hLNr1y75+/vnOvA0szt39+7dDu2lS5dW7969tXTpUs2bNy/L4/7++2+NHTtWNWvWzHWAcU7atWunrVu36rbbblODBg2yTJnBpmrVqg7t2f0H7Qx/f3/79Utq166d7XNmfrDZbDZ5e3s77Nzx8fFXPStKyt//av/9919JUo0aNa55XVfTvHlzrVixwv4Bm2nmzJny9/e3n03UrFkzbdu2TX/88YfDctldd6lChQpZLhD5999/Z7n2UevWrfX333/bD+1k53r1FlyuefPm2r59e5ZrlMycOVM2m03NmjWTZG6TM2fO2A8bZ8ppm/z9998OgS8hIcF+plOm1q1b68KFC1e9uN+VPaE5qVq1qm699VZ9/vnnDmepnDt3TvPnz7efKXUjc2Ufz463t7cqVqyY5bMwN+np6dqyZYt8fHxyHbidydn9zJl9whXOvu9y065dO+3evVvBwcHZbtsrP5uv5+dXfqLHpoA1aNBAtWvX1vvvvy/DMBQdHZ1lmc8++0xTpkxR165ddccddygwMFAHDx7U//3f/9lHrGd22WanbNmyqlixotavX68XXnjBYd57772nnTt36oknntDq1av14IMPysfHR+vXr9e7776rgIAAzZ8/3+H4rrPeeOMNxcbGKioqSi+88IKqVq2qCxcuaO/evfr+++81efJkh7OJ8sMHH3yge+65R40bN9azzz6rChUq6MyZM/rnn3+0aNEi+4dIu3bttGDBAj333HPq1KmTDhw4oDfffFNlypTRrl27cn2O2267TX5+fpo9e7aqV6+uYsWKKSwszGFsgLPWr1+v4ODg63JGwfDhw+3jnoYNG6agoCDNnj1bixcv1rhx4+xnV8TExGjq1Klq27at3nrrLftZUdmdgtq9e3c98cQTeu6559SxY0ft27dP48aNy3JGUkxMjObNm6eHH35YgwcP1p133qnz589r1apVateunZo1a6aAgACVL19e33zzjZo3b66goCCFhITkOeg6o3///po5c6batm2rN954Q+XLl9fixYs1ceJEPfvss6pSpYok84zE8ePH68knn9Tbb7+typUr6/vvv9fSpUuz3SZTpkzRE088of/85z9KSEjQuHHjsvxRfPzxxzVt2jT16dNHO3fuVLNmzZSRkaFff/1V1atX12OPPSZJioyM1MqVK7Vo0SKVKVNGAQEBDpdtyFSkSBGNGzdO3bp1U7t27fTMM88oJSVF77zzjk6fPq0xY8YUwBa8/pzdx3PStGlT/fDDDznOj4uLs+8LR48e1dSpU7Vjxw71798/10NYmVzZz662T7jC2fddbmJiYjR//nzde++96t+/v2rXrm0/O+zHH3/UgAED1LBhQ/vy69evl4eHh/1M2BuGmwcv3xQ++OADQ5JRo0aNbOdv377dGDBggNGgQQOjZMmShqenp1GiRAmjSZMmxmeffebUc7z++utGiRIlsr0IU2pqqvHRRx8ZDRs2NIoVK2b4+PgYVatWNQYNGmScOHEiy/I5nQWR3aj848ePGy+88IIRERFheHl5GUFBQUb9+vWNV1991Th79uxV687LBfr27Nlj9O7d27j11lsNLy8vo2TJkkZUVFSWi9iNGTPGqFChguHj42NUr17d+OSTT+xnRlxZw5VnpcyZM8eoVq2a4eXl5XAmT04X6MtuvRkZGUb58uWNfv36XXU75HSBPleX27Jli/Hggw8agYGBhre3t3H77bdne3bX9u3bjZYtWxq+vr5GUFCQER0dbXzzzTdZzorKyMgwxo0bZ1SsWNHw9fU1GjRoYKxYsSLb98KpU6eMF1980ShXrpzh5eVllCpVymjbtq3DRSOXLVtm1K1b1/Dx8TEk5XphOmcvXna15fbt22d07drVCA4ONry8vIyqVasa77zzTpYLXx48eNDo2LGjUaxYMSMgIMDo2LGjsXbt2mzfgzNmzDCqV69u+Pr6GjVq1DDmzZuX5ewUwzDPJhw2bJhRuXJlw9vb2wgODjbuu+8+Y+3atfZlNm/ebNx9992Gv7+/Icm+XXO6YOLXX39tNGzY0PD19TWKFi1qNG/e3Pjll18clsnpYnTZncV0NXk5K6pmzZpOr0eS0bdvX4c2Z/fx7GRewHHDhg0O7dmdFRUUFGQ0bNjQmDp1qsP7Ib/2s6vtE7m9dy//3MnkzPvuavvD2bNnjddee82oWrWq4e3tbQQGBhqRkZFG//79jfj4eIdlGzdubDz44IPZrqcwsxnGDXblHWTr8OHDioiI0MyZM6/p2hvIP8uXL1erVq20bds2h2tWAChYtWvX1t13361Jkya5u5Qb1u7du1W5cmUtXbo0T2Mw3YlgYyGvvPKKfvjhB23evDnfzsRA3jVr1kyVKlXSJ5984u5SgJvKkiVL9Mgjj2jXrl35fjj8ZtGrVy8dPHhQsbGx7i7FZR4jrrxWPG5Yd911lxITE1WpUqU8X60S+ePUqVOKj4/X4MGDVbRoUXeXA9xUKlWqpFtuuUU+Pj4OXwoM56SlpWn79u0aNGiQS1fBLyzosQEAAJbB8QoAAGAZBBsAAGAZBBsAAGAZlr9AX0ZGhg4fPqyAgIBsLzENAAAKH8MwdObMGZe/c83ywebw4cM5fpEgAAAo3A4cOODSafuWDzYBAQGSzA3jyqWnAQCA+yQlJSk8PNz+d9xZlg82mYefihcvTrABAOAG4+owErcOHh4xYoRsNpvDVLp0aft8wzA0YsQIhYWFyc/PT02bNtW2bdvcWDEAACjM3H5WVM2aNXXkyBH7tGXLFvu8cePG6b333tP//vc//fbbbypdurRatmypM2fOuLFiAABQWLk92Hh6eqp06dL2qWTJkpLM3pr3339fr776qjp06KBatWppxowZSk5O1ueff57j+lJSUpSUlOQwAdfi0Uclm82cHnvsUvvu3VL37lJ4uOTjI4WESI0bS199dfV1Nm16aZ2XT/fcU2AvAwBuCm4PNrt27VJYWJgiIiL02GOP6d9//5Uk7dmzR/Hx8WrVqpV9WR8fHzVp0kRr167NcX2jR49WYGCgfeKMKFyLadOyDyqGIbVsKc2aJR07JtWsKWVkSGvWmEFowwbn1l+xotSw4aWpZs38rR8AbjZuDTYNGzbUzJkztXTpUn3yySeKj49XVFSUEhISFB8fL0kKDQ11eExoaKh9XnaGDBmixMRE+3TgwIECfQ2wrt27pRdekBo1kq480/DQIWnPHvP2iBHSpk3S4sWX5u/f79xzvP66tH79pWnKlHwpHQBuWm4NNq1bt1bHjh0VGRmpFi1aaPH//8swY8YM+zJXjoY2DCPXEdI+Pj72M6A4Ewp5lZYmdesmFSkizZ4teXg4zi9TRqpUybw9YoRUr57Utq3k6Wkennr4Yeeep39/8zBWxYrS009LR4/m68sAgJuO2w9FXa5o0aKKjIzUrl277GdHXdk7c+zYsSy9OEB+GzlS+vVXaeJEKSIi63wPD+mnn6T69aXUVOn336VTp6RbbpHuvFPy8rr6cxQtao7PKVnS7P355BOzd+jcufx/PQBwsyhUwSYlJUV//fWXypQpo4iICJUuXVqxsbH2+ampqVq1apWioqLcWCWsbuNGafRo6YknzF6b7GRkSH36SHFx0osvSmfPSl98IZ04IfXrJ339de7P8f770smT0p9/SgcOSEOGmO179kgLF+bv6wGAm4lbg83AgQO1atUq7dmzR7/++qs6deqkpKQk9ejRQzabTTExMRo1apQWLlyorVu3qmfPnvL391fXrl3dWTYsbutWKT3dHDRcrJg5ZY6ZmT/fvL9o0aUxNT16mL0vjz4qZR75XLYs9+eoU0fy9jZv22zS5W9pZ8fnAACycuuVhw8ePKjHH39cJ06cUMmSJXXXXXdp/fr1Kl++vCRp0KBBOn/+vJ577jmdOnVKDRs21I8//ujy5ZWBvLhwIWtbWpo5Xbx4qW3jRqluXenvv6XMSywVLWr+PHRIat7cvD16tPTII+ZZVLNmSf/5j5T5Vp4379L6KlTI95cCADcNm2EYhruLKEhJSUkKDAxUYmIiA4mRZxUqSPv2SV26SHPnSgkJUuXK5riaIkWkGjXMw0jnzpnjazZsMHtl9u69NEZn2jSpZ89LbZ6e5gDkc+fMw1GSVL26eYaVr697XicAFBZ5/ftdqMbYADeK4GDpl1/MMThly0q7dpm9L23aSKtWmaEmJyVLSq++ai5z7Jg5LqdaNWnwYHOdhBoAyDt6bAAAQKFDjw0AALjpuXXwMJAX42P/dncJN5T+Lau4uwQAuG7osQEAAJZBsAEAAJZBsAEAAJZBsAEAAJZBsAEAAJZBsAEAAJZBsAEAAJZBsAEAAJZBsAEAAJZBsAEAAJZBsAEAAJZBsAEAAJZBsLnBPPqoZLOZ02OPXWpv2vRS++XTPfdcfZ0XL0rDh0sVK0re3lLZslJMjHTmTIG9DAAACgTf7n0DmTZN+uqr3JepWFEqWfLS/Zo1r77eHj2kOXOkIkWkypWlf/+VPvhA2rRJWrnSbAcA4EbAn6wbxO7d0gsvSI0amT0qOXn9dWn9+kvTlCm5rzcuzgw1khlmduyQ5s837//8s7RwYf7UDwDA9UCwuQGkpUndupk9J7NnSx4eOS9sj1j3AAAgAElEQVTbv7/k42P23Dz9tHT0aO7rXrLk0u2OHc2fbdtKvr7m7aVLr612AACuJ4LNDWDkSOnXX6WJE6WIiJyXK1pUCg83D0Xt2SN98onZw3PuXM6POXDg0u1SpcyfRYpIISHm7f37r71+AACuF4JNIbdxozR6tPTEE2avTU7ef186eVL6808zrAwZYrbv2ZP74STDyL3dZstb3QAAuAPBppDbulVKTzcHDRcrZk6ZvSjz55v3ExOlOnXMM5okM4x07XppHbn1upQrd+l25mGrjAwpIcG8HR6ef68FAICCRrC5QVy4YB5SOnfuUm9KWpp5PyVFeu89x9Oz5827dLtCBfPnoUNStWrmlNmL88ADl5bLPONq0SLz+a6cDwBAYUewKeR69jSDzOVT+fLmvC5dzPvJydKAAVJQkFS9utkL89Zb5jLVq0sdOpi3L16Udu40p8REs61+fenxx83b/fuboefRR837jRtL7dtft5cKAMA1I9hYQMmS0quvmoejjh2TTpwwA8rgwdIvv1w6wyknM2ZIw4aZgejff831vfCCtHgx17ABANxYbIaR0/BRa0hKSlJgYKASExNVvHhxd5eDfDA+9m93l3BD6d+yirtLAACX5fXvN/+PAwAAyyDYAAAAy+C7oq4Bh0RcwyERAEBBo8cGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYBsEGAABYRqEJNqNHj5bNZlNMTIy9LSUlRf369VNISIiKFi2qhx56SAcPHnRjlQAAoDArFMHmt99+08cff6zatWs7tMfExGjhwoWaO3eu1qxZo7Nnz6pdu3ZKT093U6UAAKAwc3uwOXv2rLp166ZPPvlEJUqUsLcnJibq008/1X//+1+1aNFCdevW1axZs7RlyxYtW7bMjRUDAIDCyu3Bpm/fvmrbtq1atGjh0B4XF6eLFy+qVatW9rawsDDVqlVLa9euzXF9KSkpSkpKcpgAAMDNwdOdTz537lzFxcVp48aNWebFx8fL29vboRdHkkJDQxUfH5/jOkePHq2RI0fme60AAKDwc1uPzYEDB/Tiiy9q9uzZ8vX1dfpxhmHIZrPlOH/IkCFKTEy0TwcOHMiPcgEAwA3AbcEmLi5Ox44dU/369eXp6SlPT0+tWrVKH374oTw9PRUaGqrU1FSdOnXK4XHHjh1TaGhojuv18fFR8eLFHSYAAHBzcFuwad68ubZs2aLNmzfbpwYNGqhbt272215eXoqNjbU/5siRI9q6dauioqLcVTYAACjE3DbGJiAgQLVq1XJoK1q0qIKDg+3t0dHRGjBggIKDgxUUFKSBAwcqMjIyy0BjAAAAyc2Dh69m/Pjx8vT0VOfOnXX+/Hk1b95c06dPl4eHh7tLAwAAhVChCjYrV650uO/r66sJEyZowoQJ7ikIAADcUNx+HRsAAID8QrABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACWQbABAACW4enqA1JSUrRhwwbt3btXycnJKlmypOrWrauIiIiCqA8AAMBpTgebtWvXasKECfr666+VmpqqW265RX5+fjp58qRSUlJUsWJFPf300+rTp48CAgIKsmYAAIBsOXUo6uGHH1anTp106623aunSpTpz5owSEhJ08OBBJScna9euXXrttde0fPlyValSRbGxsQVdNwAAQBZO9di0atVKX375pby9vbOdX7FiRVWsWFE9evTQtm3bdPjw4XwtEgAAwBlOBZu+ffs6vcKaNWuqZs2aeS4IAAAgr1wePJwpNTVVx44dU0ZGhkN7uXLlrrkoAACAvHA52OzatUu9e/fW2rVrHdoNw5DNZlN6enq+FQcAAOAKl4NNz5495enpqe+++05lypSRzWYriLoAAABc5nKw2bx5s+Li4lStWrWCqAcAACDPXL7ycI0aNXTixImCqAUAAOCauBxsxo4dq0GDBmnlypVKSEhQUlKSwwQAAOAuLh+KatGihSSpefPmDu0MHgYAAO7mcrD56aef8u3JJ02apEmTJmnv3r2SzGvgDBs2TK1bt5Zkfi/VwIEDNWfOHJ0/f17NmzfXxIkTVbZs2XyrAQAAWIfLwaZJkyb59uRly5bVmDFjVKlSJUnSjBkz9PDDD+v3339XzZo1FRMTo0WLFmnu3LkKDg7WgAED1K5dO8XFxcnDwyPf6gAAANaQ5wv0JScna//+/UpNTXVor127ttPrePDBBx3uv/3225o0aZLWr1+vsmXL6tNPP9Vnn31mP/w1a9YshYeHa9myZbr//vvzWjoAALAol4PN8ePH1atXL/3www/Zzs/rGJv09HR9+eWXOnfunBo1aqS4uDhdvHhRrVq1si8TFhamWrVqae3atTkGm5SUFKWkpNjvM6AZAICbh8tnRcXExOjUqVNav369/Pz8tGTJEs2YMUOVK1fWt99+63IBW7ZsUbFixeTj46M+ffpo4cKFqlGjhuLj4+Xt7a0SJUo4LB8aGqr4+Pgc1zd69GgFBgbap/DwcJdrAgAANyaXe2xWrFihb775RnfccYeKFCmi8uXLq2XLlipevLhGjx6ttm3burS+qlWravPmzTp9+rTmz5+vHj16aNWqVTkun3n2VU6GDBmil156yX4/KSmJcAMAwE3C5WBz7tw5lSpVSpIUFBSk48ePq0qVKoqMjNSmTZtcLsDb29s+eLhBgwb67bff9MEHH6hLly5KTU3VqVOnHHptjh07pqioqBzX5+PjIx8fH5frAAAANz6XD0VVrVpVO3fulCTVqVNHU6ZM0aFDhzR58mSVKVPmmgsyDEMpKSmqX7++vLy8FBsba5935MgRbd26NddgAwAAbl4u99jExMToyJEjkqThw4fr/vvv1+zZs+Xt7a3p06e7tK6hQ4eqdevWCg8P15kzZzR37lytXLlSS5YsUWBgoKKjozVgwAAFBwcrKChIAwcOVGRkpP0sKQAAgMu5HGy6detmv123bl3t3btXO3bsULly5RQSEuLSuo4eParu3bvryJEjCgwMVO3atbVkyRK1bNlSkjR+/Hh5enqqc+fO9gv0TZ8+nWvYAACAbOX5OjaZ/P39Va9evTw99tNPP811vq+vryZMmKAJEybkaf0AAODm4lSweemll/Tmm2+qaNGiDmccZee9997Ll8IAAABc5VSw+f3333Xx4kX77Zzkdho2AABAQXMq2Fz+xZf5+SWYAAAA+cnl070BAAAKK6d6bDp06OD0ChcsWJDnYgAAAK6FUz02l3/3UvHixbV8+XJt3LjRPj8uLk7Lly9XYGBggRUKAABwNU712EybNs1++5VXXlHnzp01efJk+/Vk0tPT9dxzz6l48eIFUyUAAIATXB5jM3XqVA0cONDhInkeHh566aWXNHXq1HwtDgAAwBUuB5u0tDT99ddfWdr/+usvZWRk5EtRAAAAeeHylYd79eql3r17659//tFdd90lSVq/fr3GjBmjXr165XuBAAAAznI52Lz77rsqXbq0xo8fb/8yzDJlymjQoEEaMGBAvhcIAADgLJeDTZEiRTRo0CANGjRISUlJksSgYQAAUCjk6QJ9aWlpWrZsmebMmWP/GoXDhw/r7Nmz+VocAACAK1zusdm3b58eeOAB7d+/XykpKWrZsqUCAgI0btw4XbhwQZMnTy6IOgEAAK7K5R6bF198UQ0aNNCpU6fk5+dnb3/kkUe0fPnyfC0OAADAFS732KxZs0a//PKLvL29HdrLly+vQ4cO5VthAAAArnK5xyYjI0Pp6elZ2g8ePKiAgIB8KQoAACAvXA42LVu21Pvvv2+/b7PZdPbsWQ0fPlxt2rTJ1+IAAABc4fKhqPHjx6tZs2aqUaOGLly4oK5du2rXrl0KCQnRnDlzCqJGAAAAp7gcbMLCwrR582bNmTNHmzZtUkZGhqKjo9WtWzeHwcQAAADXm8vBRpL8/PzUu3dv9e7dO7/rAQAAyDOng83q1audWu7ee+/NczEAAADXwulg07RpU/tVhg3DyHYZm82W7RlTAAAA14PTwaZEiRIKCAhQz5491b17d4WEhBRkXQAAAC5z+nTvI0eOaOzYsVq3bp0iIyMVHR2ttWvXqnjx4goMDLRPAAAA7uJ0sPH29laXLl20dOlS7dy5U7Vr19bzzz+v8PBwvfrqq0pLSyvIOgEAAK4qT9/uHR4ermHDhmnZsmWqUqWKxowZo6SkpPyuDQAAwCUuB5uUlBR9/vnnatGihWrVqqWQkBAtXrxYQUFBBVEfAACA05wePLxhwwZNmzZNc+fOVUREhHr27KkvvviCQAMAAAoNp4PNXXfdpXLlyumFF15Q/fr1JZnf9H2lhx56KP+qAwAAcIFLVx7ev3+/3nzzzRzncx0bAADgTk4Hm4yMjIKsAwAA4Jrl6awoAACAwsipYLNu3TqnV3ju3Dlt27YtzwUBAADklVPB5sknn1TLli31xRdf6OzZs9kus337dg0dOlSVKlXSpk2b8rVIAAAAZzg1xmb79u2aMmWKhg0bpm7duqlKlSoKCwuTr6+vTp06pR07dujcuXPq0KGDYmNjVatWrYKuGwAAIAungo2Xl5eef/55Pf/889q0aZN+/vln7d27V+fPn9ftt9+u/v37q1mzZlzTBgAAuJVLp3tLUr169VSvXr2CqAUAAOCacFYUAACwDIINAACwDIINAACwDIINAACwjHwJNqdPn86P1QAAAFwTl4PN2LFjNW/ePPv9zp07Kzg4WLfeeqv++OOPfC0OAADAFS4HmylTpig8PFySFBsbq9jYWP3www9q3bq1Xn755XwvEAAAwFkuX8fmyJEj9mDz3XffqXPnzmrVqpUqVKighg0b5nuBAAAAznK5x6ZEiRI6cOCAJGnJkiVq0aKFJMkwDKWnp+dvdQAAAC5wucemQ4cO6tq1qypXrqyEhAS1bt1akrR582ZVqlQp3wsEAABwlsvBZvz48apQoYIOHDigcePGqVixYpLMQ1TPPfdcvhcIAADgLJeDTWpqqgYOHJilPSYmJl8KAgAAyCuXx9iEhoaqd+/eWrNmTUHUAwAAkGcuB5s5c+YoMTFRzZs3V5UqVTRmzBgdPny4IGoDAABwicvB5sEHH9T8+fN1+PBhPfvss5ozZ47Kly+vdu3aacGCBUpLSyuIOgEAAK4qz1+pEBwcrP79++uPP/7Qe++9p2XLlqlTp04KCwvTsGHDlJycnJ91AgAAXJXLg4czxcfHa+bMmZo2bZr279+vTp06KTo6WocPH9aYMWO0fv16/fjjj/lZKwAAQK5cDjYLFizQtGnTtHTpUtWoUUN9+/bVE088oVtuucW+TJ06dVS3bt18LRQAAOBqXA42vXr10mOPPaZffvlFd9xxR7bLVKxYUa+++uo1FwcAAOCKPH1XlL+/f67L+Pn5afjw4XkuCgAAIC9cDjaXh5rz58/r4sWLDvOLFy9+7VUBAADkgctnRZ07d07PP/+8SpUqpWLFiqlEiRIOEwAAgLu4HGwGDRqkFStWaOLEifLx8dH//d//aeTIkQoLC9PMmTMLokYAAACnuBxsFi1apIkTJ6pTp07y9PRU48aN9dprr2nUqFGaPXu2S+saPXq07rjjDgUEBKhUqVJq3769du7c6bBMSkqK+vXrp5CQEBUtWlQPPfSQDh486GrZAADgJuBysDl58qQiIiIkmeNpTp48KUm65557tHr1apfWtWrVKvXt21fr169XbGys0tLS1KpVK507d86+TExMjBYuXKi5c+dqzZo1Onv2rNq1a6f09HRXSwcAABbn8uDhihUrau/evSpfvrxq1KihL774QnfeeacWLVrkcC0bZyxZssTh/rRp01SqVCnFxcXp3nvvVWJioj799FN99tlnatGihSRp1qxZCg8P17Jly3T//fe7Wj4AALAwl3tsevXqpT/++EOSNGTIEPtYm/79++vll1++pmISExMlSUFBQZKkuLg4Xbx4Ua1atbIvExYWplq1amnt2rXZriMlJUVJSUkOEwAAuDm43GPTv39/++1mzZppx44d2rhxo2677TbdfvvteS7EMAy99NJLuueee1SrVi1J5tc2eHt7ZznbKjQ0VPHx8dmuZ/To0Ro5cmSe6wAAADeuPH9XVKZy5cqpXLly11zI888/rz///FNr1qy56rKGYchms2U7b8iQIXrppZfs95OSkhQeHn7N9QEAgMLPpWCTkZGh6dOna8GCBdq7d69sNpsiIiLUqVMnde/ePcewcTX9+vXTt99+q9WrV6ts2bL29tKlSys1NVWnTp1y6LU5duyYoqKisl2Xj4+PfHx88lQHAAC4sTk9xsYwDD300EN66qmndOjQIUVGRqpmzZrat2+fevbsqUceecTlJzcMQ88//7wWLFigFStW2M+2ylS/fn15eXkpNjbW3nbkyBFt3bo1x2ADAABuXk732EyfPl2rV6/W8uXL1axZM4d5K1asUPv27TVz5kw9+eSTTj9537599fnnn+ubb75RQECAfdxMYGCg/Pz8FBgYqOjoaA0YMEDBwcEKCgrSwIEDFRkZaT9LCgAAIJPTPTZz5szR0KFDs4QaSbrvvvs0ePBgly/QN2nSJCUmJqpp06YqU6aMfZo3b559mfHjx6t9+/bq3Lmz7r77bvn7+2vRokXy8PBw6bkAAID1Od1j8+eff2rcuHE5zm/durU+/PBDl57cMIyrLuPr66sJEyZowoQJLq0bAADcfJzusTl58qRCQ0NznB8aGqpTp07lS1EAAAB54XSwSU9Pl6dnzh08Hh4eSktLy5eiAAAA8sLpQ1GGYahnz545nkqdkpKSb0UBAADkhdPBpkePHlddxpUzogAAAPKb08Fm2rRpBVkHAADANXP5SzABAAAKK4INAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINAACwDIINgEJj9WqpTRupZEnJZjOnyZMdl7l4URo+XKpYUfL2lsqWlWJipDNnrr7+ChUurffy6YknCuTlAHADT3cXAACZNm2SYmPN0HLiRPbL9OghzZkjFSkiVa4s/fuv9MEH5mNXrjTbr6Z6dal48Uv3K1XKl/IBFAL02AAoNLp3l5KSpKVLs58fF2eGGskMMzt2SPPnm/d//llauNC555k4UVq//tI0YsQ1lw6gkCDYACg0goMlP7+c5y9Zcul2x47mz7ZtJV9f83ZOgehKHTuaj6lSRRo0yAxTAKyBYAPghnHgwKXbpUqZP4sUkUJCzNv79199HYGB5ricwEBp1y7pnXek+++XMjLyv14A1x/BBsANwzByb7fZcn/8V19JCQnSH39Ihw6Zh74k83DU2rX5VycA9yHYALhhlCt36fbRo+bPjAwzrEhSeHjuj2/QQPLwMG97ekqdO1+a50xvD4DCj2AD4IbxwAOXbn/1lflz0SLpwgXH+Rs2SNWqmdOGDWbbtm3Sp59KKSnm/fT0S+uQzFPBAdz4CDYACo0FC8xTr5s2vdQ2bJjZ1q2bVL++9PjjZnv//mZwefRR837jxlL79ubt5GRp505zSk42244fl556yhxbU6uWdOut0owZ5rz77pMaNbouLxFAAeM6NgAKjaQkafdux7bjx82pbFnz/owZ5vVrZs40r2FTsqTUqZP01lu5X8OmenUzDC1fLu3bZ/bYREZKXbtKL7549fE5AG4MNsPIaTieNSQlJSkwMFCJiYkqfvkVufLB+Ni/83V9Vte/ZZV8WQ/b3TX5td0B4HrK699vDkUBAADLINgAAADLcOsYm9WrV+udd95RXFycjhw5ooULF6p95ug/SYZhaOTIkfr444916tQpNWzYUB999JFq1qzpxqqBmxOHAJ3H4T/AfdzaY3Pu3Dndfvvt+t///pft/HHjxum9997T//73P/32228qXbq0WrZsqTPOfI0vAAC46bg12LRu3VpvvfWWOnTokGWeYRh6//339eqrr6pDhw6qVauWZsyYoeTkZH3++eduqBYAgPx17pz5fWVVqkhFi5rfOh8ZKY0aZZ65l5MRI8wz+bKb0tKuW/mFUqE93XvPnj2Kj49Xq1at7G0+Pj5q0qSJ1q5dq2eeeSbbx6WkpCgl8wpcMkdVAwBQGPXte+l6SjVqSGfPSlu3Sq++Knl5SS+/nPvjQ0Kk225zbLvZL11QaINNfHy8JCk0NNShPTQ0VPv27cvxcaNHj9bIkSMLtDYAAPLDmjXmz1atzG+nT001w8qZM+b1lq6mbVtp+vQCLfGGU+jPirJdET0Nw8jSdrkhQ4YoMTHRPh24/OuAAQAoRBo3Nn/++KNUs6Z58ckzZ6SoKOmVV67++PnzJT8/qUwZM+T8/nvB1nsjKLQ9NqVLl5Zk9tyUKVPG3n7s2LEsvTiX8/HxkY+PT4HXBwDAtZo82fwi15kzpe3bzTZvb6lOHfOq2rnx8jIDjYeHtGOH9P335pW1162T6tYt+NoLq0LbYxMREaHSpUsrNjbW3paamqpVq1YpKirKjZUBAJA/xo+XPvtMuvtu6dgx88taAwKkiROlwYNzfly3buY33P/9t/TXX9KSJWZ7Sor00UfXp/bCyq09NmfPntU///xjv79nzx5t3rxZQUFBKleunGJiYjRq1ChVrlxZlStX1qhRo+Tv76+uXbu6sWoAAK5dcrL0+uuSYUgdO5o9NCVLmiHn22+lZctyfmzlyo73779fCg6WEhKk/fsLtu7Czq3BZuPGjWrWrJn9/ksvvSRJ6tGjh6ZPn65Bgwbp/Pnzeu655+wX6Pvxxx8VEBDgrpIBAMgXycmXTs2OizN/Xrhg9tpI5unfktS8uXTokPTII9Lo0Wbb2LHmN92XK2fej401Q40kVahwXcovtNwabJo2barcvoPTZrNpxIgRGjFixPUrCgCA6yAkRLr3Xmn1amn2bOnXX82Bw0ePmvN79DB/7t5tniF15Milx06aJA0ZYgYbf39zjI1khqGYmOv7OgqbQjvGBgAAq/v660sX6Dt82Dzdu2FDadYs6bnncn7c0KHSffeZy//7r1S+vDnuJi7OvB7OzazQnhUFAIDVlShhHlYaOzbnZfbuzdr29NPmhKzosQEAAJZBsAEAAJbBoSgAAK4wPvZvd5dww+jfsoq7S3BAjw0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0AALAMgg0A3ORWr5batJFKlpRsNnOaPNm5x37+uVSvnuTnJwUFSZ06Sbt2FWy9QG4INgBwk9u0SYqNNYOJKz7+WOrWTfr9d6lMGSk9XZo/X7r7bunw4YKpFbgagg0A3OS6d5eSkqSlS51/TEqKNHSoebtjR+nff6W//pICAqTjx6XRowumVuBqCDYAcJMLDjYPJbli40YpIcG83bGj+TMsTLrrLvO2KyEJyE8EGwCAyw4cuHS7VKlLt0NDzZ/791/feoBMBBsAgMsMI/d2m+361QJcjmADAHBZuXKXbh89eun2sWPmz/Dw61sPkIlgAwC4qubNpWrVpCFDzPt33GGOzZHMM6Ek6dAhad068/YDD1z/GgGJYAMAN70FC6RKlaSmTS+1DRtmtnXrZt7fvVvauVM6csS87+0tjRp16fEVK0o1akhnz0ohIdLgwdf1JQB2BBsAuMklJZnBZd++S23Hj5tthw7l/Linn5ZmzZLq1DGvW2OzSR06SGvXmmdIAe7g6e4CAADu1bOnOeVm797s27t1u9SrAxQG9NgAAADLINgAAADL4FAUABRi42P/dncJN4z+Lau4uwQUAvTYAAAAyyDYAAAAy7ghgs3EiRMVEREhX19f1a9fXz///LO7SwIAAIVQoQ828+bNU0xMjF599VX9/vvvaty4sVq3bq39fMMaAAC4QqEPNu+9956io7F7b0IAABZiSURBVKP11FNPqXr16nr//fcVHh6uSZMmubs0AABQyBTqs6JSU1MVFxenwVdcm7tVq1Zau3Ztto9JSUlRSkqK/X5iYqIkKSkpKd/ru3DubL6v08ry63fAdncN2/36y8/PG7a789ju7lEQf18vX6+R01fJ56BQB5sTJ04oPT1doaGhDu2hoaGKj4/P9jGjR4/WyJEjs7SH81WzbjfU3QXcpNju1x/b3D3Y7u5R0Nv9zJkzCgwMdHr5Qh1sMtlsNof7hmFkacs0ZMgQvfTSS/b7GRkZOnnypIKDg3N8jJUkJSUpPDxcBw4cUPHixd1dzk2D7e4ebHf3YLu7x8223Q3D0JkzZxTm4hePFepgExISIg8Pjyy9M8eOHcvSi5PJx8dHPj4+Dm233HJLgdVYWBUvXvymeOMXNmx392C7uwfb3T1upu3uSk9NpkI9eNjb21v169dXbGysQ3tsbKyioqLcVBUAACisCnWPjSS99NJL6t69uxo0aKBGjRrp448/1v79+9WnTx93lwYAAAoZjxEjRoxwdxG5qVWrloKDgzVq1Ci9++67On/+vD777DPdfvvt7i6t0PLw8FDTpk3l6Vnoc6ulsN3dg+3uHmx392C7X53NcPU8KgAAgEKqUI+xAQAAcAXBBgAAWAbBBgAAWAbBBgAAWAbBBgAu07NnT7Vv397dZdxUpk+fbvkLqY4YMUKhoaGy2Wz6+uuv3V1OFitXrpTNZtPp06fdXco1I9jcAOLj49WvXz9VrFhRPj4+Cg8P14MPPqjly5dLkipUqCCbzSabzSY/Pz9Vq1ZN77zzTrZfHDZ//nw1bdpUgYGBKlasmGrXrq033nhDJ0+evN4vq9B78MEH1aJFi2znrVu3TjabTZs2bZLNZtPmzZslSXv37rX/Lmw2m/5fe+ceFdV1/fEvDAyDDPIQBJ9URVGxgqPmwarAqAStVtaSBBJJokIIFcTXwi4UKtHGpE1ijCbGaBoesRgBEw1JG1w1OMjD0JYIQRhFBTWpwzIWTCWEiMz39we/ueHKjDx81LDOZ637xz13n3PP2edx9913nxlHR0f4+voiISEBZ8+evZ/Vv2OWLVsmtcPW1hZjx45FUlISvv/+ewA/tdXGxgb//ve/ZXkNBgNsbGxgZWWFCxcuyOT7q5uuedVqNfz8/JCZmXm3m40dO3bck3L7wpUrVxAXF4fRo0fDzs4Onp6eCA0NxYkTJyQZ07w/cOBAt/y+vr6wsrKStaPrOtH1+OMf/2ixHsHBwZKcUqnEuHHjsGHDBtkfDd8NIiMjUVdXd1fL7C23jnMPDw+EhIQgPT0dRqPxrtxDr9dj8+bN2LNnDwwGA+bPn99rA7pr/WxsbDB69GisWLECzc3Nd6VuJgICAmAwGPr1S78PGsKwecC5cOECpk+fjsLCQrzyyiuorq5GQUEBtFotEhISJLktW7bAYDBAr9cjKSkJGzduxN69e2VlpaSkIDIyEjNnzsRnn32GU6dOYdu2baiqqsK+ffvud9MeeGJiYlBYWIiLFy92u5aeng5/f3+4urqazXv06FEYDAZUVVXhpZdegl6vh5+fn2SM/lyYN28eDAYD6uvr8eKLL+Ltt99GUlKSTGb48OF4//33ZWlZWVkYMWKE2TLvRDcZGRlS3sjISCxfvhxHjhzpfwPN4OTk9D/3HoSHh6OqqgpZWVmoq6tDfn4+goODu72AjBo1ChkZGbK0L774Ao2NjXBwcOhWrmmd6HokJibeti6xsbEwGAw4d+4cXnnlFezatQt3++fP7O3tMXTo0LtaZl8wjfMLFy7gs88+g1arxerVq7Fw4ULcvHnTYr729vZelX/+/HkAQFhYGDw9Pbv97U9f6vfnP/8Zn3zyCeLj4/tURk8olUp4enoOjP9UpOCBZv78+RwxYgRbWlq6XWtubiZJenl5cfv27bJrGo2Gixcvls7Ly8sJgG+88YbZ+5jKEvxEe3s7PTw8+MILL8jSv//+ezo6OvLNN99kQ0MDAfDkyZMk2e3cREdHB4ODg+nl5cWbN2/etzbcCUuXLmVYWJgs7bnnnqOnpyfJn9qamprK8ePHy+R8fHz4+9//ngDY0NAgk++vbgDw0KFDsjRXV1euW7dOllZaWspZs2ZRpVJx5MiRTExMlOZPcnIyH3744W5l//KXv+SmTZvMtttoNPJPf/oTx4wZQ5VKxalTpzIvL0+6rtFo+Nprr0nnYWFhVCgU/O6770iSBoOBAHj69GmLbetKc3MzAVCn091WzsvLi8nJybSzs+OlS5ek9NjYWCYmJtLJyYkZGRky+VvXiZ4ICgri6tWrZWmLFy+mRqORpX3zzTeMiIigs7MzXV1duWjRIqnfCwoKaGdn122NSUxMZGBgIEkyIyODTk5Osuv5+fnUaDS0s7PjmDFj+MILL7C9vZ0kuW7dOi5cuFCS3b59OwHw008/ldImTJjAd955p8c2mhvnJPn5558TAN99910pDQB3797NRYsWcdCgQdy0aRNv3rzJ6Oho/uIXv6BKpeKECRNk62xaWhoByA5zaceOHet1/datW0dXV1dZ2rVr1xgbG0t3d3c6OjpSq9WysrKSJHn69GkCoF6vl+XZtm0bvby8aDQaeezYMQKQ9dPt5tLOnTs5ZcoUSfbQoUMEwLfeektKe+yxx5icnGy2XfcS4bF5gGlqakJBQQESEhLMvn2Ze6skCZ1OB71eD1tbWyk9OzsbarXaopX/v35DfRCxsbHBs88+i8zMTNlnvby8PNy4cQNRUVG9Lsva2hqrV6/GxYsXUVFRcS+qe1+wt7fv9pa6aNEiNDc3o6SkBABQUlKCpqYm/OY3v+lVmf3RTUdHB3Jzc9HU1CQb59XV1QgNDcXixYvx1VdfIScnByUlJVi5ciUAICoqCuXl5dIbNADU1NSgurraYn+mpqYiIyMDu3fvRk1NDdauXYunn34aRUVFADo/1+h0OgCd86+4uBguLi6SPo4dOwZPT0/4+Pj0qm1qtRpqtRqHDx/u8ZOPh4cHQkNDkZWVBQBobW1FTk4OoqOje3WvvlJVVYXS0lKZzltbW6HVaqFWq3H8+HGUlJRArVZj3rx5uHHjBubOnQtnZ2d8+OGHUh5T/1nS+ZEjR/D0009j1apVqK2txZ49e5CZmYmtW7cC6NR5cXGx9KmoqKgIbm5uUp80Njairq4OQUFB/W7r7Nmz4efnh48++kiWnpaWhrCwMFRXVyM6OhpGoxEjR45Ebm4uamtrsWnTJmzcuBG5ubkAgKSkJMmrZvKSJSUlISIiQvLEGAyGXv//YX19PQoKCmR9QBILFixAY2Mj/va3v6GiogIajQZz5sxBU1MTfHx8MH36dGRnZ8vK2r9/P5YsWWLWS9PTXAoODkZNTQ2uXr0KoHsf3Lx5E2VlZXfUB/3mvptSgl5j8rJ89NFHt5Xz8vKiUqmkg4MDbW1tCYAqlYqlpaWSzPz58zl16tR7XeUBh16vJwAWFhZKaYGBgXzqqadIdvdCWPJKdC0rJyfn/lT+Drn1TbG8vJxDhgxhREQESXlb16xZw+XLl5Mkly9fzrVr1/LkyZO98tiQvdONaVw7ODhQoVAQAF1dXXn27FlJ5plnnuHzzz8vy1dcXExra2v+8MMPJMmpU6dyy5Yt0vUNGzZw5syZZtvd0tJClUrFsrIyWZkxMTHSGMjPz6eTkxM7OjpYWVlJd3d3rl27luvXrydJPv/884yMjLTYLnMcPHiQLi4uVKlUDAgI4IYNG1hVVSWTMXlgDh8+zHHjxtFoNDIrK4vTpk0jSbMeG9M60fWw5CkgOz02tra2dHBwoFKpJABaW1vz4MGDksx7771HHx8fGo1GKe3HH3+kvb09jxw5QpJctWoVZ8+eLV0/cuQIlUolm5qaSHb32MyaNYsvvfSSrC779u3jsGHDSHZ6J6ytrfmvf/2LRqORQ4YM4csvvyz14/79++nh4dGjnknLHhuSjIyM5KRJk6RzAFyzZk2PZcbHxzM8PFw6N3kzenvfW+UUCgUdHByoUqkkD8/rr78uyXz++eccPHgw29raZHnHjRvHPXv2kCRff/11jh07Vrp25swZAmBNTQ1JdvPY9DSXjEYj3dzcpLHg7+/Pl19+mUOHDiVJlpWV0cbGhtevX++xjXcb4bF5gOH/ewl6881z/fr1qKysRFFREbRaLVJSUmRvACQHxrfT+8zEiRMREBCA9PR0AJ3fyouLi/v1RtyX/nxQ+PTTT6FWq6FSqfDoo48iMDAQb775Zje5mJgY5OXlobGxEXl5eX3WT291s337dlRWVuLvf/87/P39sX37dnh7e0vXKyoqkJmZKXk91Go1QkNDYTQa0dDQAKDTa2N6cyWJDz74wKLnoLa2Fm1tbQgJCZGV+f7770ten8DAQFy/fh0nT55EUVERgoKCoNVqpTdXnU7X57fW8PBwXL58Gfn5+QgNDYVOp4NGozEb1LxgwQK0tLTg+PHjSE9Pv63uTetE1+Phhx++bV2ioqJQWVmJEydOICIiAtHR0QgPD5euV1RU4Ny5c3B0dJT04+rqira2NklHUVFR0Ol0uHz5MoBOD/Kvf/1ruLi4mL1nRUUFtmzZItO5KdantbUVTk5O8Pf3h06nQ3V1NaytrREXF4eqqipcv369Xzo3h7l1c8aMGd3k3nnnHcyYMQPu7u5Qq9V49913cenSpTu+vwmtVovKykqUl5cjMTERoaGhstioiooKtLS0YMiQITKdNTQ0SH3w5JNP4uLFi/jiiy8AdPaBv78/Jk+ebPaePc0lKysrBAYGQqfT4dq1a6ipqcFvf/tbdHR0QK/XS2NWrVbfNT30FvEvWg8w48ePh5WVFfR6fY/R825ubvD29oa3tzc+/PBDeHt745FHHpF29UyYMAElJSVob2+XuTAFPRMTE4OVK1di165dyMjIgJeXF+bMmdPncvR6PQBgzJgxd7uK9wytVovdu3fD1tYWw4cPtzh2pkyZgokTJ+Kpp57CpEmTMGXKFGmnWG/orW48PT2lcZ6Xl4dp06ZhxowZ0uJsNBoRFxeHVatWdcs7evRoAMCSJUuQnJyML7/8Ej/88AO+/vprPPnkk2bvZ/rU8de//rVbMLQpALTrQ7asrAyzZ8/GrFmzUFlZibNnz6Kurg7BwcG91oUJlUqFkJAQhISEYNOmTXjuueeQlpaGZcuWyeRsbGzwzDPPIC0tDeXl5Th06JDFMk3rRF9wcnKS8vzlL3+Br68v3nvvPcTExADo1JG5zxwA4O7uDgB46KGHMG7cOBw4cAArVqzAoUOHugU9d8VoNGLz5s1YvHhxt2sqlQrAT58AlUolgoKC4OLiAl9fX5SWlkKn02HNmjV9aqc59Hp9tzF5a1hAbm4u1q5di23btuHRRx+Fo6MjXn31VZSXl9/x/bve09QHO3fuhFarxebNm/GHP/wBQKe+hg0bJn0S7YopzGDYsGHQarXYv38/HnnkEXzwwQeIi4uzeM/ezKXg4GDs3bsXxcXF8PPzg7OzMwIDA1FUVASdTtevcX83EB6bBxhXV1eEhoZi165d0hbbrlj6vQEXFxckJiYiKSlJehNesmQJWlpa8Pbbb5vNMxB+u+BeERERAYVCgf379yMrKwvLly/vs9fFaDRi586dGDNmDKZNm3aPanr3MS2oXl5ePRrE0dHR0Ol0ffbW9Fc33t7eCA8Px4YNG6Q0jUaDmpoayfjpeiiVSgDAyJEjERgYiOzsbGRnZ2Pu3Lnw8PAwe4/JkyfDzs4Oly5d6lbeqFGjJLng4GAcO3YMx48fR3BwMJydnTF58mS8+OKLGDp0KCZNmtQnnViqi7l1AOjUfVFREcLCwix6Qe4Gtra22LhxI1JTU9Ha2gqgU+dnz57F0KFDu+mo69bhJUuWIDs7G5988gmsra2xYMECi/fRaDQ4c+aM2X60tu58bJnibAoLC6UHaFBQEA4cOHDH8TUAUFhYiOrqapl3yhzFxcUICAhAfHw8pk2bBm9vb1kMlyWUSiU6Ojr6Vbe0tDS89tprkgdMo9GgsbERNjY23fTl5uYm5YuKikJOTg5OnDiB8+fPWzToTWX2NJdMcTYHDx6U9cHRo0f/d/E1gIixedCpr6+np6cnJ0+ezIMHD7Kuro61tbXcsWMHJ06cSNL8bocrV65QpVLJdm/87ne/o0Kh4Pr161lWVsYLFy7w6NGjfPzxxy3ulhJ0EhMTQxcXF1pbW/PixYtSuqUYm6NHj9JgMPD8+fP8+OOPqdVqaW9vL4vVedDpKQbg1ra3t7fz22+/lXauWIqx6a9uYGZX1FdffUUrKyv+85//JElWVVXR3t6e8fHxPHnyJOvq6vjxxx9z5cqVsnx79+7l8OHD6ebmxn379t223SkpKRwyZAgzMzN57tw5fvnll3zrrbeYmZkpyeTn51OhUNDd3V2KNVmzZg0VCgWfeOKJ27brVq5evUqtVst9+/axqqqK9fX1zM3NpYeHB6OjoyW5W+f91atX2draKp2bi7HZsmULDQaD7DDt3jKHuV1RP/74I4cNG8ZXX32VZOcuwfHjxzM4OJjHjx9nfX09dTodV61axa+//lrKV1dXRwCcOnUqY2JiZGXeGmNTUFBAGxsbpqWl8dSpU6ytreWBAweYkpIiyZjibBQKBU+dOkWSPHz4sNQPvWXp0qWcN28eDQYDv/nmG1ZUVHDr1q1Uq9VcuHChbKeeuTH4xhtvcPDgwSwoKOCZM2eYmprKwYMH08/PT5IxF2OzdetWjh49mqdPn+a3337LGzduWKyfuXk4ffp0JiQkkOzcuferX/2Kfn5+LCgoYENDA0tLS5mSkiLNDZL87rvvqFKp6Ofnxzlz5sjKuzXGpjdzyRRno1AopB1plZWVVCgUsp2B9xth2PwMuHz5MhMSEqTgvxEjRnDRokVS0J+lbZyxsbH09fVlR0eHlJaTk8PAwEA6OjrSwcFBCqQU271vT1lZGQHwsccek6VbMmxMx6BBgzhp0iTGx8fLglx/DvTVsLkVS4ZNf3Vj7qFCkiEhIZw/f750/o9//IMhISFUq9XSGN+6dassT3NzM+3s7Dho0KBuwY3mtnvv2LGDPj4+tLW1pbu7O0NDQ1lUVCTJXLt2jQqFgo8//riUZm77a29oa2tjcnIyNRoNnZycOGjQIPr4+DA1NVVmuPS0fducYdNV/6YjLi7OYhnmDBuy86Hs7u4u6c5gMPDZZ5+lm5sb7ezsOHbsWMbGxnZ7sM2cObNbMD5pfrt3QUEBAwICaG9vz8GDB/Ohhx7i3r17ZTLTp0+XGZP/+c9/aGVlJeuHnli6dKmkCxsbG7q7u3Pu3LlMT0+XrZ2k+THY1tbGZcuW0cnJic7OzlyxYgWTk5N7NGyuXLkijVP0cbs3SWZnZ1OpVEpb/f/73/8yMTGRw4cPp62tLUeNGsWoqCjZTwGQ5BNPPEEATE9Pl6Wb2+7dm7kUHh4uM2KMRiNdXV05Y8YMs+25H1iRZn6eViAQCAQCgeBniIixEQgEAoFAMGAQho1AIBAIBIIBgzBsBAKBQCAQDBiEYSMQCAQCgWDAIAwbgUAgEAgEAwZh2AgEAoFAIBgwCMNGIBAIBALBgEEYNgKBQCAQCAYMwrARCAQCgUAwYBCGjUAgEAgEggGDMGwEAoFAIBAMGP4Ph+Vjb4uXz3oAAAAASUVORK5CYII=\n",
"text/plain": [
"