Contexto

Para Arca Continental su principal canal de distribución es el canal tradicional, es decir, las tienditas de la esquina. Esto permite que la familia de productos de la compañía Coca Cola estén siempre cerca de sus consumidores a través de estas pequeñas empresas familiares que forman parte de su propia comunidad.

Sin embargo, este tipo de formato de Retail está enfrentando una fuerte competencia por parte de los canales modernos. En un principio fueron las tiendas de conveniencia como Oxxo y 7-Eleven quienes comenzaron a crear una importante cantidad de nuevos establecimientos, cada vez más en zonas habitacionales y no solo en avenidas o gasolineras.

Instalar paquetes y llamar librerías

#install.packages("readxl") #Para importar un archivo de Excel
library(readxl)

#install.packages("DataExplorer") #Para realizar Análisis descriptivo
library(DataExplorer)

#install.packages("dplyr")
library(dplyr)

#install.packages("tidyverse") #Paquete global para manipulación y análisis de datos
library(tidyverse)

#install.packages("forecast") #Para revisar tendencias y generar pronósticos
library(forecast)

#install.packages("ggplot2") #Para generar gráficas con mejor diseño
library(ggplot2)

Importar la base de datos

#file.choose()
df <- read_excel("C:\\Users\\serva\\Downloads\\Datos Arca Continental Original.xlsx")
## Warning: Expecting numeric in N184066 / R184066C14: got a date
## Warning: Expecting numeric in O184066 / R184066C15: got a date
## Warning: Expecting numeric in P184066 / R184066C16: got a date
## Warning: Expecting numeric in Q184066 / R184066C17: got a date
## Warning: Expecting numeric in R184066 / R184066C18: got a date
## Warning: Expecting numeric in S184066 / R184066C19: got a date
## Warning: Expecting numeric in T184066 / R184066C20: got a date
## Warning: Expecting numeric in U184066 / R184066C21: got a date
## Warning: Expecting numeric in V184066 / R184066C22: got a date
## Warning: Expecting numeric in W184066 / R184066C23: got a date
## Warning: Expecting numeric in X184066 / R184066C24: got a date
## Warning: Expecting numeric in Y184066 / R184066C25: got a date
colnames(df) <- make.names(colnames(df))

Análisis descriptivo

summary(df)
##        ID              Año        Territorio        Sub.Territorio    
##  Min.   :     1   Min.   :2016   Length:466509      Length:466509     
##  1st Qu.:116628   1st Qu.:2017   Class :character   Class :character  
##  Median :233255   Median :2018   Mode  :character   Mode  :character  
##  Mean   :233255   Mean   :2018                                        
##  3rd Qu.:349882   3rd Qu.:2019                                        
##  Max.   :466509   Max.   :2019                                        
##                                                                       
##      CEDI             Cliente             Nombre          Tamaño.Cte.Industria
##  Length:466509      Length:466509      Length:466509      Length:466509       
##  Class :character   Class :character   Class :character   Class :character    
##  Mode  :character   Mode  :character   Mode  :character   Mode  :character    
##                                                                               
##                                                                               
##                                                                               
##                                                                               
##  Segmento.Det          Marca           Presentacion          Tamaño         
##  Length:466509      Length:466509      Length:466509      Length:466509     
##  Class :character   Class :character   Class :character   Class :character  
##  Mode  :character   Mode  :character   Mode  :character   Mode  :character  
##                                                                             
##                                                                             
##                                                                             
##                                                                             
##  Retornable_NR          Enero            Febrero             Marzo         
##  Length:466509      Min.   :  -19.0   Min.   :  -11.00   Min.   :  -32.00  
##  Class :character   1st Qu.:    1.0   1st Qu.:    1.00   1st Qu.:    1.00  
##  Mode  :character   Median :    2.0   Median :    2.00   Median :    3.00  
##                     Mean   :   10.1   Mean   :    9.76   Mean   :   11.36  
##                     3rd Qu.:    6.0   3rd Qu.:    6.00   3rd Qu.:    6.00  
##                     Max.   :42736.0   Max.   :42767.00   Max.   :42795.00  
##                     NA's   :233480    NA's   :231213     NA's   :227420    
##      Abril               Mayo              Junio             Julio         
##  Min.   :  -70.00   Min.   : -106.00   Min.   : -211.0   Min.   :  -60.00  
##  1st Qu.:    1.00   1st Qu.:    1.00   1st Qu.:    1.0   1st Qu.:    1.00  
##  Median :    3.00   Median :    3.00   Median :    3.0   Median :    2.00  
##  Mean   :   11.71   Mean   :   12.75   Mean   :   12.2   Mean   :   11.75  
##  3rd Qu.:    6.00   3rd Qu.:    7.00   3rd Qu.:    6.0   3rd Qu.:    6.00  
##  Max.   :42826.00   Max.   :42856.00   Max.   :42887.0   Max.   :42917.00  
##  NA's   :224057     NA's   :216910     NA's   :215753    NA's   :223411    
##      Agosto           Septiembre         Octubre          Noviembre      
##  Min.   : -211.00   Min.   : -527.0   Min.   :  -38.0   Min.   :  -25.0  
##  1st Qu.:    1.00   1st Qu.:    1.0   1st Qu.:    1.0   1st Qu.:    1.0  
##  Median :    3.00   Median :    3.0   Median :    3.0   Median :    3.0  
##  Mean   :   11.98   Mean   :   13.4   Mean   :   13.7   Mean   :   13.3  
##  3rd Qu.:    6.00   3rd Qu.:    7.0   3rd Qu.:    7.0   3rd Qu.:    6.0  
##  Max.   :42948.00   Max.   :42979.0   Max.   :43009.0   Max.   :43040.0  
##  NA's   :220242     NA's   :337314    NA's   :338386    NA's   :338460   
##    Diciembre      
##  Min.   :  -28.0  
##  1st Qu.:    1.0  
##  Median :    3.0  
##  Mean   :   14.8  
##  3rd Qu.:    7.0  
##  Max.   :43070.0  
##  NA's   :341855
str(df)
## tibble [466,509 × 25] (S3: tbl_df/tbl/data.frame)
##  $ ID                  : num [1:466509] 1 2 3 4 5 6 7 8 9 10 ...
##  $ Año                 : num [1:466509] 2016 2016 2016 2016 2016 ...
##  $ Territorio          : chr [1:466509] "Guadalajara" "Guadalajara" "Guadalajara" "Guadalajara" ...
##  $ Sub.Territorio      : chr [1:466509] "Belenes" "Belenes" "Belenes" "Belenes" ...
##  $ CEDI                : chr [1:466509] "Suc. Belenes" "Suc. Belenes" "Suc. Belenes" "Suc. Belenes" ...
##  $ Cliente             : chr [1:466509] "77737" "77737" "77737" "77737" ...
##  $ Nombre              : chr [1:466509] "ABARR" "ABARR" "ABARR" "ABARR" ...
##  $ Tamaño.Cte.Industria: chr [1:466509] "Extra Grande" "Extra Grande" "Extra Grande" "Extra Grande" ...
##  $ Segmento.Det        : chr [1:466509] "Agua Mineral" "Agua Purificada" "Agua Purificada" "Agua Saborizada" ...
##  $ Marca               : chr [1:466509] "Topo Chico A.M." "Ciel Agua Purificada" "Ciel Agua Purificada" "Ciel Exprim" ...
##  $ Presentacion        : chr [1:466509] "600 ml NR" "1 Ltro. N.R." "1.5 Lts. NR" "600 ml NR" ...
##  $ Tamaño              : chr [1:466509] "Individual" "Individual" "Individual" "Individual" ...
##  $ Retornable_NR       : chr [1:466509] "No Retornable" "No Retornable" "No Retornable" "No Retornable" ...
##  $ Enero               : num [1:466509] NA NA NA NA NA NA 1 NA 3 NA ...
##  $ Febrero             : num [1:466509] NA 2 NA NA NA NA NA 1 3 NA ...
##  $ Marzo               : num [1:466509] NA 8 3 NA NA 1 NA NA 4 NA ...
##  $ Abril               : num [1:466509] NA 4 6 NA NA NA NA 1 4 NA ...
##  $ Mayo                : num [1:466509] NA 4 3 NA NA NA 0 NA 4 NA ...
##  $ Junio               : num [1:466509] NA 2 3 NA NA NA NA 1 4 0 ...
##  $ Julio               : num [1:466509] NA 2 3 NA NA NA 0 NA 4 NA ...
##  $ Agosto              : num [1:466509] NA 2 3 NA NA NA NA 1 7 NA ...
##  $ Septiembre          : num [1:466509] NA 2 3 NA NA NA NA 1 4 NA ...
##  $ Octubre             : num [1:466509] NA 2 3 NA NA NA 0 NA 3 NA ...
##  $ Noviembre           : num [1:466509] NA 4 3 NA 0 NA NA NA 1 NA ...
##  $ Diciembre           : num [1:466509] 1 2 3 1 NA NA NA NA 3 NA ...
create_report(df)
## 
## 
## processing file: report.rmd
##   |                                             |                                     |   0%  |                                             |.                                    |   2%                                   |                                             |..                                   |   5% [global_options]                  |                                             |...                                  |   7%                                   |                                             |....                                 |  10% [introduce]                       |                                             |....                                 |  12%                                   |                                             |.....                                |  14% [plot_intro]                      |                                             |......                               |  17%                                   |                                             |.......                              |  19% [data_structure]                  |                                             |........                             |  21%                                   |                                             |.........                            |  24% [missing_profile]                 |                                             |..........                           |  26%                                   |                                             |...........                          |  29% [univariate_distribution_header]  |                                             |...........                          |  31%                                   |                                             |............                         |  33% [plot_histogram]                  |                                             |.............                        |  36%                                   |                                             |..............                       |  38% [plot_density]                    |                                             |...............                      |  40%                                   |                                             |................                     |  43% [plot_frequency_bar]              |                                             |.................                    |  45%                                   |                                             |..................                   |  48% [plot_response_bar]               |                                             |..................                   |  50%                                   |                                             |...................                  |  52% [plot_with_bar]                   |                                             |....................                 |  55%                                   |                                             |.....................                |  57% [plot_normal_qq]                  |                                             |......................               |  60%                                   |                                             |.......................              |  62% [plot_response_qq]                |                                             |........................             |  64%                                   |                                             |.........................            |  67% [plot_by_qq]                      |                                             |..........................           |  69%                                   |                                             |..........................           |  71% [correlation_analysis]            |                                             |...........................          |  74%                                   |                                             |............................         |  76% [principal_component_analysis]    |                                             |.............................        |  79%                                   |                                             |..............................       |  81% [bivariate_distribution_header]   |                                             |...............................      |  83%                                   |                                             |................................     |  86% [plot_response_boxplot]           |                                             |.................................    |  88%                                   |                                             |.................................    |  90% [plot_by_boxplot]                 |                                             |..................................   |  93%                                   |                                             |...................................  |  95% [plot_response_scatterplot]       |                                             |.................................... |  98%                                   |                                             |.....................................| 100% [plot_by_scatterplot]           
## output file: C:/Users/serva/Downloads/report.knit.md
## "C:/Program Files/RStudio/resources/app/bin/quarto/bin/tools/pandoc" +RTS -K512m -RTS "C:\Users\serva\Downloads\report.knit.md" --to html4 --from markdown+autolink_bare_uris+tex_math_single_backslash --output pandoc582861e7226c.html --lua-filter "C:\Users\serva\AppData\Local\R\win-library\4.4\rmarkdown\rmarkdown\lua\pagebreak.lua" --lua-filter "C:\Users\serva\AppData\Local\R\win-library\4.4\rmarkdown\rmarkdown\lua\latex-div.lua" --embed-resources --standalone --variable bs3=TRUE --section-divs --table-of-contents --toc-depth 6 --template "C:\Users\serva\AppData\Local\R\win-library\4.4\rmarkdown\rmd\h\default.html" --no-highlight --variable highlightjs=1 --variable theme=yeti --mathjax --variable "mathjax-url=https://mathjax.rstudio.com/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML" --include-in-header "C:\Users\serva\AppData\Local\Temp\RtmpwbXS3o\rmarkdown-str5828642136a4.html"
## 
## Output created: report.html
introduce(df)
## # A tibble: 1 × 9
##     rows columns discrete_columns continuous_columns all_missing_columns
##    <int>   <int>            <int>              <int>               <int>
## 1 466509      25               11                 14                   0
## # ℹ 4 more variables: total_missing_values <int>, complete_rows <int>,
## #   total_observations <int>, memory_usage <dbl>
plot_intro(df)

plot_missing(df)

plot_histogram(df)

plot_bar(df)
## 4 columns ignored with more than 50 categories.
## Cliente: 5249 categories
## Nombre: 1090 categories
## Marca: 56 categories
## Presentacion: 57 categories

plot_correlation(df)
## 5 features with more than 20 categories ignored!
## Cliente: 5249 categories
## Nombre: 1090 categories
## Segmento.Det: 21 categories
## Marca: 56 categories
## Presentacion: 57 categories

count(df, Territorio, sort=TRUE)
## # A tibble: 2 × 2
##   Territorio       n
##   <chr>        <int>
## 1 Guadalajara 466508
## 2 Territorio       1
count(df, Sub.Territorio, sort=TRUE)
## # A tibble: 4 × 2
##   Sub.Territorio      n
##   <chr>           <int>
## 1 Belenes        208982
## 2 Huentitán      144196
## 3 Toluquilla     113330
## 4 Sub Territorio      1
count(df, CEDI, sort=TRUE)
## # A tibble: 4 × 2
##   CEDI                 n
##   <chr>            <int>
## 1 Suc. Belenes    208982
## 2 Suc. Huentitán  144196
## 3 Suc. Toluquilla 113330
## 4 CEDI                 1
count(df, Cliente, sort=TRUE)
## # A tibble: 5,249 × 2
##    Cliente     n
##    <chr>   <int>
##  1 0286      647
##  2 2912      586
##  3 2661      537
##  4 7821      531
##  5 1859      525
##  6 5583      516
##  7 9998      508
##  8 3601      506
##  9 5879      499
## 10 0335      496
## # ℹ 5,239 more rows
count(df, Nombre, sort=TRUE)
## # A tibble: 1,090 × 2
##    Nombre     n
##    <chr>  <int>
##  1 ABARR  71186
##  2 MARIA  39816
##  3 JOSE   17479
##  4 JUAN    7580
##  5 MARTH   5759
##  6 MISCE   5700
##  7 LUIS    5585
##  8 SUPER   4565
##  9 CARLO   3991
## 10 ROSA    3890
## # ℹ 1,080 more rows
count(df, Tamaño.Cte.Industria, sort=TRUE)
## # A tibble: 5 × 2
##   Tamaño.Cte.Industria      n
##   <chr>                 <int>
## 1 Extra Grande         230190
## 2 Micro                117110
## 3 Pequeño               77875
## 4 Grande                41333
## 5 Tamaño Cte Industria      1
count(df, Segmento.Det, sort=TRUE)
## # A tibble: 21 × 2
##    Segmento.Det            n
##    <chr>               <int>
##  1 Sabores Regular    156242
##  2 Colas Regular       95720
##  3 Colas Light         43807
##  4 Jugos y Néctares    33362
##  5 Bebidas de Fruta    30641
##  6 Agua Purificada     20766
##  7 Agua Mineral        12590
##  8 Isotónicos Regular  11905
##  9 Té Regular          10062
## 10 Agua Saborizada     10056
## # ℹ 11 more rows
count(df, Marca, sort=TRUE)
## # A tibble: 56 × 2
##    Marca                    n
##    <chr>                <int>
##  1 Coca-Cola            95720
##  2 Sprite               37925
##  3 Fanta                35728
##  4 Fresca               26435
##  5 Manzana Lift         25598
##  6 Coca-Cola Light      21926
##  7 Del Valle            21325
##  8 Ciel Agua Purificada 20766
##  9 Sidral Mundet        17150
## 10 Valle Frut           15808
## # ℹ 46 more rows
count(df, Presentacion, sort=TRUE)
## # A tibble: 57 × 2
##    Presentacion         n
##    <chr>            <int>
##  1 600 ml NR        74008
##  2 1 Ltro. N.R.     36930
##  3 2 Lts. NR        36415
##  4 500 ml Ret       35165
##  5 1.5 Lts. NR      30637
##  6 Lata 235 ml      24551
##  7 400 ml NR        22877
##  8 250 ml. NR PET   21735
##  9 500 ml NR Vidrio 18758
## 10 2.5 Lts. NR      13235
## # ℹ 47 more rows
count(df, Tamaño, sort=TRUE)
## # A tibble: 3 × 2
##   Tamaño          n
##   <chr>       <int>
## 1 Individual 328513
## 2 Familiar   137995
## 3 Tamaño          1

Limpieza de base de datos

## Revisar NAs en la base de datos
# ¿Cuántos NAs tengo en la base de datos?
sum(is.na(df))
## [1] 3148501
# ¿Cuántos NAs tengo por variable?
sapply(df, function(x) sum(is.na(x)))
##                   ID                  Año           Territorio 
##                    0                    0                    0 
##       Sub.Territorio                 CEDI              Cliente 
##                    0                    0                    0 
##               Nombre Tamaño.Cte.Industria         Segmento.Det 
##                    0                    0                    0 
##                Marca         Presentacion               Tamaño 
##                    0                    0                    0 
##        Retornable_NR                Enero              Febrero 
##                    0               233480               231213 
##                Marzo                Abril                 Mayo 
##               227420               224057               216910 
##                Junio                Julio               Agosto 
##               215753               223411               220242 
##           Septiembre              Octubre            Noviembre 
##               337314               338386               338460 
##            Diciembre 
##               341855
#Reemplazar NAs con CEROS
df[is.na(df)] <- 0
sum(is.na(df))
## [1] 0
## Detectar valores atípicos
boxplot(df$Enero)

# Eliminar renglón de los totales
df <- df[df$Enero <6000, ]
boxplot(df$Enero)

Análisis de ventas

#Crear una nueva base de datos
df1 <- df
# Muestra las ventas de Enero a Junio por CEDI:
df2 <- select(df1,c(CEDI,Enero:Junio))

# Muestra los movimientos por Cedi y tamaño de tienda grande:
df3 <- df1 %>% filter(Tamaño.Cte.Industria == "Grande")

# Ordena la base de datos por Cedi, por marca y por presentación:
df4 <- df1 %>% arrange(CEDI, Marca, Presentacion)

# Agrega un campo calculado con las ventas del primer semestre y muestra las ventas del primer semestre por marca:
df5 <- df1 %>% 
  mutate(Ventas_Sem1 = Enero + Febrero + Marzo + Abril + Mayo + Junio)

ventas_sem1_por_marca <- df5 %>% 
  group_by(Marca) %>% 
  summarise(Ventas_Sem1=sum(Ventas_Sem1))

# Obtén la media de las ventas del primer semestre agrupado por marca, presentación y tamaño.
df6 <- df5 %>%
  group_by(Marca, Presentacion, Tamaño) %>%
  summarise(Ventas_Sem1=mean(Ventas_Sem1))
## `summarise()` has grouped output by 'Marca', 'Presentacion'. You can override
## using the `.groups` argument.
# Calcular medidas de tendencia central
summary(df1)
##        ID              Año        Territorio        Sub.Territorio    
##  Min.   :     1   Min.   :2016   Length:466508      Length:466508     
##  1st Qu.:116628   1st Qu.:2017   Class :character   Class :character  
##  Median :233256   Median :2018   Mode  :character   Mode  :character  
##  Mean   :233255   Mean   :2018                                        
##  3rd Qu.:349882   3rd Qu.:2019                                        
##  Max.   :466509   Max.   :2019                                        
##      CEDI             Cliente             Nombre          Tamaño.Cte.Industria
##  Length:466508      Length:466508      Length:466508      Length:466508       
##  Class :character   Class :character   Class :character   Class :character    
##  Mode  :character   Mode  :character   Mode  :character   Mode  :character    
##                                                                               
##                                                                               
##                                                                               
##  Segmento.Det          Marca           Presentacion          Tamaño         
##  Length:466508      Length:466508      Length:466508      Length:466508     
##  Class :character   Class :character   Class :character   Class :character  
##  Mode  :character   Mode  :character   Mode  :character   Mode  :character  
##                                                                             
##                                                                             
##                                                                             
##  Retornable_NR          Enero             Febrero             Marzo         
##  Length:466508      Min.   : -19.000   Min.   : -11.000   Min.   : -32.000  
##  Class :character   1st Qu.:   0.000   1st Qu.:   0.000   1st Qu.:   0.000  
##  Mode  :character   Median :   0.000   Median :   0.000   Median :   0.000  
##                     Mean   :   4.951   Mean   :   4.829   Mean   :   5.729  
##                     3rd Qu.:   2.000   3rd Qu.:   2.000   3rd Qu.:   3.000  
##                     Max.   :5333.000   Max.   :4995.000   Max.   :5636.000  
##      Abril               Mayo             Junio              Julio         
##  Min.   : -70.000   Min.   :-106.00   Min.   :-211.000   Min.   : -60.000  
##  1st Qu.:   0.000   1st Qu.:   0.00   1st Qu.:   0.000   1st Qu.:   0.000  
##  Median :   0.000   Median :   0.00   Median :   0.000   Median :   0.000  
##  Mean   :   5.992   Mean   :   6.73   Mean   :   6.464   Mean   :   6.033  
##  3rd Qu.:   3.000   3rd Qu.:   3.00   3rd Qu.:   3.000   3rd Qu.:   3.000  
##  Max.   :6164.000   Max.   :6759.00   Max.   :6033.000   Max.   :6735.000  
##      Agosto           Septiembre          Octubre           Noviembre      
##  Min.   :-211.000   Min.   :-527.000   Min.   : -38.000   Min.   : -25.00  
##  1st Qu.:   0.000   1st Qu.:   0.000   1st Qu.:   0.000   1st Qu.:   0.00  
##  Median :   0.000   Median :   0.000   Median :   0.000   Median :   0.00  
##  Mean   :   6.235   Mean   :   3.625   Mean   :   3.674   Mean   :   3.57  
##  3rd Qu.:   3.000   3rd Qu.:   1.000   3rd Qu.:   1.000   3rd Qu.:   1.00  
##  Max.   :6065.000   Max.   :6509.000   Max.   :6326.000   Max.   :5319.00  
##    Diciembre       
##  Min.   : -28.000  
##  1st Qu.:   0.000  
##  Median :   0.000  
##  Mean   :   3.858  
##  3rd Qu.:   0.000  
##  Max.   :6182.000
# Colapsar meses en una columna
df7 <- gather(df1, Mes, Ventas, Enero:Diciembre)
df7
## # A tibble: 5,598,096 × 15
##       ID   Año Territorio  Sub.Territorio CEDI         Cliente Nombre
##    <dbl> <dbl> <chr>       <chr>          <chr>        <chr>   <chr> 
##  1     1  2016 Guadalajara Belenes        Suc. Belenes 77737   ABARR 
##  2     2  2016 Guadalajara Belenes        Suc. Belenes 77737   ABARR 
##  3     3  2016 Guadalajara Belenes        Suc. Belenes 77737   ABARR 
##  4     4  2016 Guadalajara Belenes        Suc. Belenes 77737   ABARR 
##  5     5  2016 Guadalajara Belenes        Suc. Belenes 77737   ABARR 
##  6     6  2016 Guadalajara Belenes        Suc. Belenes 77737   ABARR 
##  7     7  2016 Guadalajara Belenes        Suc. Belenes 77737   ABARR 
##  8     8  2016 Guadalajara Belenes        Suc. Belenes 77737   ABARR 
##  9     9  2016 Guadalajara Belenes        Suc. Belenes 77737   ABARR 
## 10    10  2016 Guadalajara Belenes        Suc. Belenes 77737   ABARR 
## # ℹ 5,598,086 more rows
## # ℹ 8 more variables: Tamaño.Cte.Industria <chr>, Segmento.Det <chr>,
## #   Marca <chr>, Presentacion <chr>, Tamaño <chr>, Retornable_NR <chr>,
## #   Mes <chr>, Ventas <dbl>
# Agregar una columna con el Número de Mes
meses <- c("Enero"="01", "Febrero"="02","Marzo"="03", "Abril"="04", "Mayo"="05", "Junio"="06", "Julio"="07", "Agosto"="08", "Septiembre"="09", "Octubre"="10", "Noviembre"="11", "Diciembre"="12")

df7$Número_de_Mes <- meses[df7$Mes]

# Graficar el total de ventas por mes y año

ventas.totales <- df7 %>%
  group_by(Año, Número_de_Mes) %>%
  summarise(Ventas.Totales = sum(Ventas)) %>%
  arrange(Número_de_Mes)
## `summarise()` has grouped output by 'Año'. You can override using the `.groups`
## argument.
ggplot(ventas.totales, aes(x=Número_de_Mes, y=Ventas.Totales,group=Año, color= as.factor(Año))) +
  geom_line() +
  geom_point() +
  labs(x="Mes",y="Ventas Totales(Qty)", color = "Año")

Regresión Lineal

# Elaborar Regresión Lineal

ventas_4567 <- df7 %>%
  filter(Cliente == "4567") %>%
  group_by(Año) %>%
  summarise(Ventas_Totales = sum(Ventas)) %>%
  filter(Ventas_Totales != 0)

# ventas_4567$Secuencia <- 1:nrow(ventas_4567)

ggplot(ventas_4567, aes(x=Año, y=Ventas_Totales)) +
  geom_point() +
  labs(x="Año",y="Ventas Totales(Qty)", title= "Ventas Totales del Cliente 4567")

regresion <- lm(Ventas_Totales ~ Año, data=ventas_4567)
summary(regresion)
## 
## Call:
## lm(formula = Ventas_Totales ~ Año, data = ventas_4567)
## 
## Residuals:
##     1     2     3 
## -5092 10185 -5092 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)
## (Intercept) 48918677   17799125   2.748    0.222
## Año           -24174       8820  -2.741    0.223
## 
## Residual standard error: 12470 on 1 degrees of freedom
## Multiple R-squared:  0.8825, Adjusted R-squared:  0.765 
## F-statistic: 7.512 on 1 and 1 DF,  p-value: 0.2227
# Ecuación
# y = 48918677 - 24174 * Año 

# R cuadrada ajustada
# 77%

datos <- data.frame(Año=2020:2025)
prediccion <- predict(regresion,datos)
prediccion
##          1          2          3          4          5          6 
##  87197.333  63023.333  38849.333  14675.333  -9498.667 -33672.667
ventas_4567$Tipo_de_Dato <- "Datos Reales"
datos$Ventas_Totales <- prediccion
datos$Tipo_de_Dato <- "Predicción"

datos_combinados <- rbind(ventas_4567,datos)

ggplot(datos_combinados, aes(x=Año, y=Ventas_Totales, color=Tipo_de_Dato)) +
  geom_point() +
  labs(x="Año",y="Ventas Totales(Qty)", title= "Pronóstico a 5 años de Ventas del Cliente 4567")

Series de tiempo

ventas_mensuales_4567 <- df7 %>%
  filter(Cliente == "4567") %>%
  group_by(Año, Número_de_Mes) %>%
  summarise(Ventas_Totales = sum(Ventas)) %>%
  filter(Ventas_Totales != 0)

# Confirmar que los datos que queremos modelar esten ordenados cronológicamente.

# Función de Serie de Tiempo MENSUAL, que inicia en Enero 2017
ts <- ts(data=ventas_mensuales_4567$Ventas_Totales, start = c(2017,1), frequency = 12)

# Función de Serie de Tiempo MENSUAL, que inicia en Abril 2017
# ts <- ts(data=ventas_mensuales_4567$Ventas_Totales, start = c(2017,4), frequency = 12)

# Función de Serie de Tiempo TRIMESTRAL, que inicia en Enero 2017
# ts <- ts(data=ventas_mensuales_4567$Ventas_Totales, start = c(2017,1), frequency = 4)

# Función de Serie de Tiempo TRIMESTRAL, que inicia en Octubre 2017 (Q4)
# ts <- ts(data=ventas_mensuales_4567$Ventas_Totales, start = c(2017,4), frequency = 4)

# Función de Serie de Tiempo ANUAL, que inicia en 2017
# ts <- ts(data=ventas_mensuales_4567$Ventas_Totales, start = 2017, frequency = 1)

# Crear Modelo ARIMA
# Modelo Autorregresivo Integrado de Promedio Movil.
arima <- auto.arima(ts, D=1) # D=1 por la temporalidad
arima
## Series: ts 
## ARIMA(0,0,0)(0,1,0)[12] 
## 
## sigma^2 = 10383172:  log likelihood = -189.94
## AIC=381.87   AICc=382.09   BIC=382.87
summary(arima)
## Series: ts 
## ARIMA(0,0,0)(0,1,0)[12] 
## 
## sigma^2 = 10383172:  log likelihood = -189.94
## AIC=381.87   AICc=382.09   BIC=382.87
## 
## Training set error measures:
##                     ME     RMSE     MAE       MPE     MAPE      MASE      ACF1
## Training set -36.98041 2547.446 1572.02 -1.811657 13.08459 0.6269271 0.2519497
# Generar el pronóstico de ventas
pronostico <- forecast(arima, level=95, h=36)
pronostico
##          Point Forecast      Lo 95    Hi 95
## Sep 2019          12086  5770.4219 18401.58
## Oct 2019          11427  5111.4219 17742.58
## Nov 2019          11270  4954.4219 17585.58
## Dec 2019          12227  5911.4219 18542.58
## Jan 2020           9430  3114.4219 15745.58
## Feb 2020          11312  4996.4219 17627.58
## Mar 2020          12515  6199.4219 18830.58
## Apr 2020          13334  7018.4219 19649.58
## May 2020          16286  9970.4219 22601.58
## Jun 2020          15347  9031.4219 21662.58
## Jul 2020          14005  7689.4219 20320.58
## Aug 2020          14050  7734.4219 20365.58
## Sep 2020          12086  3154.4238 21017.58
## Oct 2020          11427  2495.4238 20358.58
## Nov 2020          11270  2338.4238 20201.58
## Dec 2020          12227  3295.4238 21158.58
## Jan 2021           9430   498.4238 18361.58
## Feb 2021          11312  2380.4238 20243.58
## Mar 2021          12515  3583.4238 21446.58
## Apr 2021          13334  4402.4238 22265.58
## May 2021          16286  7354.4238 25217.58
## Jun 2021          15347  6415.4238 24278.58
## Jul 2021          14005  5073.4238 22936.58
## Aug 2021          14050  5118.4238 22981.58
## Sep 2021          12086  1147.0979 23024.90
## Oct 2021          11427   488.0979 22365.90
## Nov 2021          11270   331.0979 22208.90
## Dec 2021          12227  1288.0979 23165.90
## Jan 2022           9430 -1508.9021 20368.90
## Feb 2022          11312   373.0979 22250.90
## Mar 2022          12515  1576.0979 23453.90
## Apr 2022          13334  2395.0979 24272.90
## May 2022          16286  5347.0979 27224.90
## Jun 2022          15347  4408.0979 26285.90
## Jul 2022          14005  3066.0979 24943.90
## Aug 2022          14050  3111.0979 24988.90
plot(pronostico, main="Ventas Mensuales y Pronóstico a 3 Años del Cliente 4567", xlab="Año", ylab="Ventas (Qty)")

Preguntas de Evidencia

Pregunta 1

Preguntas Detonantes

  1. ¿Puede observarse un crecimiento en las ventas de algunos de los segmentos de productos de la familia Coca Cola en las tiendas en las que se implementó el Proyecto Siglo XXI de Arca Continental?
ventas_totales_seg <- df7 %>%
   filter(Segmento.Det == "Agua Purificada" | Segmento.Det == "Isotónicos Regular" | Segmento.Det == "Colas Regular" ) %>%
  filter(Año == 2018) %>%
  group_by(Segmento.Det, Número_de_Mes) %>%
  summarise(Ventas_Totales = sum(Ventas)) %>%
  arrange(Número_de_Mes)
## `summarise()` has grouped output by 'Segmento.Det'. You can override using the
## `.groups` argument.
ggplot(ventas_totales_seg, aes(x=Número_de_Mes, y=Ventas_Totales,group=Segmento.Det, color= as.factor(Segmento.Det))) +
  geom_line() +
  geom_point() +
  labs(x="Mes",y="Ventas Totales(Qty)", title= "Ventas Totales por Segmento", color = "Segmento")

Respuesta: En general las ventas aumentaron en 2019, pero cada segmento se comportó diferente. El incremento de ventas se puede lograr gracias a mejores campañas de mercadotecnia y lanzamiento de nuevos productos.En 2018 la venta de Coca-Cola entre Febrero y Diciembre aumentó en 100,000 unidades, mientras que las aguas isotónicas no tuvieron cambios.

Pregunta 2

  1. ¿El incremento en las ventas es similar entre los diferentes tamaños de clientes?
ventas_totales_tam_cliente <- df7 %>%
  filter(Año == 2019) %>%
  group_by(Tamaño.Cte.Industria, Número_de_Mes) %>%
  summarise(Ventas_Totales = sum(Ventas)) %>%
  arrange(Número_de_Mes) %>%
  filter(Ventas_Totales != 0)
## `summarise()` has grouped output by 'Tamaño.Cte.Industria'. You can override
## using the `.groups` argument.
ggplot(ventas_totales_tam_cliente, aes(x=Número_de_Mes, y=Ventas_Totales,group=Tamaño.Cte.Industria, color= as.factor(Tamaño.Cte.Industria))) +
  geom_line() +
  geom_point() +
  labs(x="Mes",y="Ventas Totales(Qty)", title= "Ventas Totales por Tamaño de Cliente", color = "Tamaño de Cliente")

Respuesta: El incremento de ventas es similar independientemente de los tamaños de los clientes. Como todos los negocios están abiertos los mismos días, y están ubicados en la misma ciudad (Guadalajara), los factores externos de las ventas efectan casi igual a todos. El mes de mayo de 2019 fue el de mayor ventas en todos los tamaños de cliente, y el primer bimestre fue el más bajo también para todos los tamaños de cliente.

Pregunta 3

  1. ¿Cuál es el comportamiento observado de las unidades vendidas por mes de cada una de las marcas, independientemente de sus respectivas presentaciones?
ventas_totales_marca <- df7 %>%
  filter(Año == 2019) %>%
  filter(Marca == c("Coca-Cola", "Valle Frut", "Ciel Agua Purificada","Coca-Cola Light")) %>%
  group_by(Marca, Número_de_Mes) %>%
  summarise(Ventas_Totales = sum(Ventas)) %>%
  arrange(Número_de_Mes) %>%
  filter(Ventas_Totales != 0)
## `summarise()` has grouped output by 'Marca'. You can override using the
## `.groups` argument.
ggplot(ventas_totales_marca, aes(x=Número_de_Mes, y=Ventas_Totales,group=Marca, color= as.factor(Marca))) +
  geom_line() +
  geom_point() +
  labs(x="Mes",y="Ventas Totales(Qty)", title= "Ventas Totales por Marca", color = "Marca")

Respuesta: El comportamiento de ventas es similar entre las marcas a través de los meses, pero hay marcas que venden mucho más que otras. Los meses de más ventas de bebidas son los meses en que la temperatura es mayor, y las marcas más conocidas son las que tiene en mente el consumidor. El mes de mayo 2019 es para todas las marcas el de mayores ventas, pero la marca Coca-Cola es la que vende más que todas.

Pregunta 4

  1. ¿Se ha incrementado la venta de productos en envases retornables en los últimos dos años?
ventas_totales_envase <- df7 %>%
  filter(Retornable_NR=="Retornable") %>%
  filter(Año == c(2017,2018,2019)) %>%
  group_by(Año, Número_de_Mes) %>%
  summarise(Ventas_Totales = sum(Ventas)) %>%
  arrange(Número_de_Mes) %>%
  filter(Ventas_Totales != 0)
## `summarise()` has grouped output by 'Año'. You can override using the `.groups`
## argument.
ggplot(ventas_totales_envase, aes(x=Número_de_Mes, y=Ventas_Totales,group=Año, color= as.factor(Año))) +
  geom_line() +
  geom_point() +
  labs(x="Mes",y="Ventas Totales(Qty)", title= "Ventas Totales de Envases Retornables", color = "Año")

Respuesta: La venta de productos en envase retornable aumentó en los últimos dos años. Las campañas de conciencia ambiental han vuelto populares las opciones menos dañinas al ecosistema. De 2017 a 2018 hubo un ligero incremento en las ventas de envase retornable, y entre 2018 y 2019 el aumento fue muy notorio.

Notas: Hay otros factores a considerar como que en 2019 hay más tiendas con el proyecto Siglo XXI, o que las tiendas en general venden más productos (tanto retornables como no retornables).

Pregunta 5

  1. ¿El comportamiento de la venta de agua ha incrementado en relación al de los refrescos o las bebidas isotónicas?
ventas_totales_seg <- df7 %>%
   filter(Segmento.Det == "Agua Purificada" | Segmento.Det == "Isotónicos Regular" | Segmento.Det == "Colas Regular" ) %>%
  filter(Año == 2018) %>%
  group_by(Segmento.Det, Número_de_Mes) %>%
  summarise(Ventas_Totales = sum(Ventas)) %>%
  arrange(Número_de_Mes)
## `summarise()` has grouped output by 'Segmento.Det'. You can override using the
## `.groups` argument.
ggplot(ventas_totales_seg, aes(x=Número_de_Mes, y=Ventas_Totales,group=Segmento.Det, color= as.factor(Segmento.Det))) +
  geom_line() +
  geom_point() +
  labs(x="Mes",y="Ventas Totales(Qty)", title= "Ventas Totales por Segmento", color = "Segmento")

Respuesta: Las ventas de Agua Purificada no ha aumentado en relación a Refrescos o Bebidas Isotónicas. La persepción que el agua no tiene valor agregado en la preparación y el azúcar en la receta hacen de los refrescos la opción preferida de los consumidores mexicanos. En 2018 las ventas totales de Agua Purificada no superaron las 50,000 unidades, mientras que la Coca-Cola vendió en su mes más alto 450,000 (8 veces más que el agua).

Pregunta 6

  1. ¿Puede decirse que la venta mensual de agua está relacionada con la venta mensual de refrescos en los últimos 4 años?
ventas_totales_ayr <- df7 %>%
   filter(Segmento.Det == c("Agua Purificada", "Colas Regular")) %>%
  group_by(Segmento.Det, Año) %>%
  summarise(Ventas_Totales = sum(Ventas)) 
## `summarise()` has grouped output by 'Segmento.Det'. You can override using the
## `.groups` argument.
ggplot(ventas_totales_ayr, aes(x=Año, y=Ventas_Totales,group=Segmento.Det, color= as.factor(Segmento.Det))) +
  geom_line() +
  geom_point() +
  labs(x="Año",y="Ventas Totales(Qty)", title= "Ventas Totales por Segmento", color = "Segmento")

Respuesta: La venta de agua y refrescos se incrementa en los mismos meses, y a través de los años, pero no se puede concluir que un incremento cause el otro. Los meses calurosos y el crecimiento de la población hace que las ventas se incrementen en las bebidas, sean agua o refresco. En la gráfica obtenida se observa un incremento año con año en la venta de agua y también de refrescos.

Pregunta 7

  1. ¿A cuánto ascienden las ventas esperadas para el 2020 en la Coca Cola de 500 ml NR Vidrio?
ventas_coca_vidrio <- df7 %>%
  filter(Marca == "Coca-Cola" & Presentacion == "500 ml NR Vidrio") %>%
  group_by(Año, Número_de_Mes) %>%
  summarise(Ventas_Totales = sum(Ventas)) %>%
  filter(Ventas_Totales != 0)
## `summarise()` has grouped output by 'Año'. You can override using the `.groups`
## argument.
# Confirmar que los datos que queremos modelar esten ordenados cronológicamente.

# Función de Serie de Tiempo MENSUAL, que inicia en Enero 2017
ts <- ts(data=ventas_coca_vidrio$Ventas_Totales, start = c(2016,1), frequency = 12)

# Crear Modelo ARIMA
# Modelo Autorregresivo Integrado de Promedio Movil.
arima <- auto.arima(ts, D=1) # D=1 por la temporalidad
arima
## Series: ts 
## ARIMA(1,0,0)(0,1,0)[12] with drift 
## 
## Coefficients:
##          ar1     drift
##       0.7092  719.6837
## s.e.  0.1328  299.8533
## 
## sigma^2 = 41602728:  log likelihood = -325.42
## AIC=656.84   AICc=657.7   BIC=661.24
summary(arima)
## Series: ts 
## ARIMA(1,0,0)(0,1,0)[12] with drift 
## 
## Coefficients:
##          ar1     drift
##       0.7092  719.6837
## s.e.  0.1328  299.8533
## 
## sigma^2 = 41602728:  log likelihood = -325.42
## AIC=656.84   AICc=657.7   BIC=661.24
## 
## Training set error measures:
##                    ME     RMSE      MAE      MPE     MAPE      MASE       ACF1
## Training set 48.85447 5325.927 3245.824 -1.37305 8.085494 0.3668297 -0.1551829
# Generar el pronóstico de ventas
pronostico <- forecast(arima, level=95, h=16)
pronostico
##          Point Forecast    Lo 95    Hi 95
## Sep 2019       57009.75 44367.95 69651.56
## Oct 2019       47611.80 32113.20 63110.40
## Nov 2019       43891.36 27138.87 60643.85
## Dec 2019       49988.01 32639.01 67337.02
## Jan 2020       52649.58 35008.14 70291.03
## Feb 2020       59116.44 41329.71 76903.18
## Mar 2020       63214.35 45354.98 81073.72
## Apr 2020       45929.53 28033.74 63825.33
## May 2020       77893.53 59979.44 95807.63
## Jun 2020       66916.74 48993.45 84840.03
## Jul 2020       65103.76 47175.85 83031.68
## Aug 2020       66475.17 48544.93 84405.41
## Sep 2020       65845.94 43788.54 87903.34
## Oct 2020       56389.84 32524.67 80255.02
## Nov 2020       52628.16 27903.54 77352.78
## Dec 2020       58695.56 33549.72 83841.40
plot(pronostico, main="Ventas Mensuales y Pronóstico para 2020 de Coca-Cola 500 ml NR Vidrio", xlab="Año",  ylab="Ventas (Qty)")

Respuesta: Las ventas esperadas de Coca-Cola de 500 ml NR Vidrio para 2020 oscilan mensualmente entre 45,000 y 65,000 unidades, con una confiabilidad del 95%.

LS0tDQp0aXRsZTogIkFyY2EgQ29udGluZW50YWwgUHJveWVjdG8gU2lnbG8gWFhJIg0KYXV0aG9yOiAiU2VydmFuZG9CRF9BMDA4MzQ2NDciDQpkYXRlOiAiMjAyNC0wOS0xMiINCm91dHB1dDogDQogIGh0bWxfZG9jdW1lbnQ6DQogICAgdG9jOiBUUlVFDQogICAgdG9jX2Zsb2F0OiBUUlVFDQogICAgY29kZV9kb3dubG9hZDogVFJVRQ0KICAgIHRoZW1lOiBjZXJ1bGVhbg0KLS0tDQoNCmBgYHtyIHNldHVwLCBpbmNsdWRlPUZBTFNFfQ0Ka25pdHI6Om9wdHNfY2h1bmskc2V0KGVjaG8gPSBUUlVFKQ0KYGBgDQoNCiFbXShDOlxcVXNlcnNcXHNlcnZhXFxEb3dubG9hZHNcXEFyY2EtQ29udGluZW50YWwtYWJyaXJhLXVuYS1udWV2YS1wbGFudGEtZW4tRXN0YWRvcy1Vbmlkb3MtR1IuanBnKQ0KDQojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+Q29udGV4dG88L3NwYW4+DQpQYXJhICoqQXJjYSBDb250aW5lbnRhbCoqIHN1IHByaW5jaXBhbCBjYW5hbCBkZSBkaXN0cmlidWNpw7NuIGVzIGVsIGNhbmFsIHRyYWRpY2lvbmFsLCBlcyBkZWNpciwgbGFzIHRpZW5kaXRhcyBkZSBsYSBlc3F1aW5hLiBFc3RvIHBlcm1pdGUgcXVlIGxhIGZhbWlsaWEgZGUgcHJvZHVjdG9zIGRlIGxhIGNvbXBhw7HDrWEgQ29jYSBDb2xhIGVzdMOpbiBzaWVtcHJlIGNlcmNhIGRlIHN1cyBjb25zdW1pZG9yZXMgYSB0cmF2w6lzIGRlIGVzdGFzIHBlcXVlw7FhcyBlbXByZXNhcyBmYW1pbGlhcmVzIHF1ZSBmb3JtYW4gcGFydGUgZGUgc3UgcHJvcGlhIGNvbXVuaWRhZC4NCg0KU2luIGVtYmFyZ28sIGVzdGUgdGlwbyBkZSBmb3JtYXRvIGRlIFJldGFpbCBlc3TDoSBlbmZyZW50YW5kbyB1bmEgZnVlcnRlIGNvbXBldGVuY2lhIHBvciBwYXJ0ZSBkZSBsb3MgY2FuYWxlcyBtb2Rlcm5vcy4gRW4gdW4gcHJpbmNpcGlvIGZ1ZXJvbiBsYXMgdGllbmRhcyBkZSBjb252ZW5pZW5jaWEgY29tbyBPeHhvIHkgNy1FbGV2ZW4gcXVpZW5lcyBjb21lbnphcm9uIGEgY3JlYXIgdW5hIGltcG9ydGFudGUgY2FudGlkYWQgZGUgbnVldm9zIGVzdGFibGVjaW1pZW50b3MsIGNhZGEgdmV6IG3DoXMgZW4gem9uYXMgaGFiaXRhY2lvbmFsZXMgeSBubyBzb2xvIGVuIGF2ZW5pZGFzIG8gZ2Fzb2xpbmVyYXMuDQoNCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWU7Ij5JbnN0YWxhciBwYXF1ZXRlcyB5IGxsYW1hciBsaWJyZXLDrWFzPC9zcGFuPg0KYGBge3IgbWVzc2FnZT1GQUxTRSwgd2FybmluZz1GQUxTRX0NCiNpbnN0YWxsLnBhY2thZ2VzKCJyZWFkeGwiKSAjUGFyYSBpbXBvcnRhciB1biBhcmNoaXZvIGRlIEV4Y2VsDQpsaWJyYXJ5KHJlYWR4bCkNCg0KI2luc3RhbGwucGFja2FnZXMoIkRhdGFFeHBsb3JlciIpICNQYXJhIHJlYWxpemFyIEFuw6FsaXNpcyBkZXNjcmlwdGl2bw0KbGlicmFyeShEYXRhRXhwbG9yZXIpDQoNCiNpbnN0YWxsLnBhY2thZ2VzKCJkcGx5ciIpDQpsaWJyYXJ5KGRwbHlyKQ0KDQojaW5zdGFsbC5wYWNrYWdlcygidGlkeXZlcnNlIikgI1BhcXVldGUgZ2xvYmFsIHBhcmEgbWFuaXB1bGFjacOzbiB5IGFuw6FsaXNpcyBkZSBkYXRvcw0KbGlicmFyeSh0aWR5dmVyc2UpDQoNCiNpbnN0YWxsLnBhY2thZ2VzKCJmb3JlY2FzdCIpICNQYXJhIHJldmlzYXIgdGVuZGVuY2lhcyB5IGdlbmVyYXIgcHJvbsOzc3RpY29zDQpsaWJyYXJ5KGZvcmVjYXN0KQ0KDQojaW5zdGFsbC5wYWNrYWdlcygiZ2dwbG90MiIpICNQYXJhIGdlbmVyYXIgZ3LDoWZpY2FzIGNvbiBtZWpvciBkaXNlw7FvDQpsaWJyYXJ5KGdncGxvdDIpDQpgYGANCg0KIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZTsiPkltcG9ydGFyIGxhIGJhc2UgZGUgZGF0b3M8L3NwYW4+DQpgYGB7cn0NCiNmaWxlLmNob29zZSgpDQpkZiA8LSByZWFkX2V4Y2VsKCJDOlxcVXNlcnNcXHNlcnZhXFxEb3dubG9hZHNcXERhdG9zIEFyY2EgQ29udGluZW50YWwgT3JpZ2luYWwueGxzeCIpDQpjb2xuYW1lcyhkZikgPC0gbWFrZS5uYW1lcyhjb2xuYW1lcyhkZikpDQpgYGANCg0KIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZTsiPkFuw6FsaXNpcyBkZXNjcmlwdGl2bzwvc3Bhbj4NCmBgYHtyfQ0Kc3VtbWFyeShkZikNCnN0cihkZikNCmNyZWF0ZV9yZXBvcnQoZGYpDQppbnRyb2R1Y2UoZGYpDQpwbG90X2ludHJvKGRmKQ0KcGxvdF9taXNzaW5nKGRmKQ0KcGxvdF9oaXN0b2dyYW0oZGYpDQpwbG90X2JhcihkZikNCnBsb3RfY29ycmVsYXRpb24oZGYpDQoNCmNvdW50KGRmLCBUZXJyaXRvcmlvLCBzb3J0PVRSVUUpDQpjb3VudChkZiwgU3ViLlRlcnJpdG9yaW8sIHNvcnQ9VFJVRSkNCmNvdW50KGRmLCBDRURJLCBzb3J0PVRSVUUpDQpjb3VudChkZiwgQ2xpZW50ZSwgc29ydD1UUlVFKQ0KY291bnQoZGYsIE5vbWJyZSwgc29ydD1UUlVFKQ0KY291bnQoZGYsIFRhbWHDsW8uQ3RlLkluZHVzdHJpYSwgc29ydD1UUlVFKQ0KY291bnQoZGYsIFNlZ21lbnRvLkRldCwgc29ydD1UUlVFKQ0KY291bnQoZGYsIE1hcmNhLCBzb3J0PVRSVUUpDQpjb3VudChkZiwgUHJlc2VudGFjaW9uLCBzb3J0PVRSVUUpDQpjb3VudChkZiwgVGFtYcOxbywgc29ydD1UUlVFKQ0KYGBgDQoNCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWU7Ij5MaW1waWV6YSBkZSBiYXNlIGRlIGRhdG9zPC9zcGFuPg0KYGBge3J9DQojIyBSZXZpc2FyIE5BcyBlbiBsYSBiYXNlIGRlIGRhdG9zDQojIMK/Q3XDoW50b3MgTkFzIHRlbmdvIGVuIGxhIGJhc2UgZGUgZGF0b3M/DQpzdW0oaXMubmEoZGYpKQ0KDQojIMK/Q3XDoW50b3MgTkFzIHRlbmdvIHBvciB2YXJpYWJsZT8NCnNhcHBseShkZiwgZnVuY3Rpb24oeCkgc3VtKGlzLm5hKHgpKSkNCg0KI1JlZW1wbGF6YXIgTkFzIGNvbiBDRVJPUw0KZGZbaXMubmEoZGYpXSA8LSAwDQpzdW0oaXMubmEoZGYpKQ0KDQojIyBEZXRlY3RhciB2YWxvcmVzIGF0w61waWNvcw0KYm94cGxvdChkZiRFbmVybykNCg0KIyBFbGltaW5hciByZW5nbMOzbiBkZSBsb3MgdG90YWxlcw0KZGYgPC0gZGZbZGYkRW5lcm8gPDYwMDAsIF0NCmJveHBsb3QoZGYkRW5lcm8pDQpgYGANCg0KIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZTsiPkFuw6FsaXNpcyBkZSB2ZW50YXM8L3NwYW4+DQpgYGB7cn0NCiNDcmVhciB1bmEgbnVldmEgYmFzZSBkZSBkYXRvcw0KZGYxIDwtIGRmDQojIE11ZXN0cmEgbGFzIHZlbnRhcyBkZSBFbmVybyBhIEp1bmlvIHBvciBDRURJOg0KZGYyIDwtIHNlbGVjdChkZjEsYyhDRURJLEVuZXJvOkp1bmlvKSkNCg0KIyBNdWVzdHJhIGxvcyBtb3ZpbWllbnRvcyBwb3IgQ2VkaSB5IHRhbWHDsW8gZGUgdGllbmRhIGdyYW5kZToNCmRmMyA8LSBkZjEgJT4lIGZpbHRlcihUYW1hw7FvLkN0ZS5JbmR1c3RyaWEgPT0gIkdyYW5kZSIpDQoNCiMgT3JkZW5hIGxhIGJhc2UgZGUgZGF0b3MgcG9yIENlZGksIHBvciBtYXJjYSB5IHBvciBwcmVzZW50YWNpw7NuOg0KZGY0IDwtIGRmMSAlPiUgYXJyYW5nZShDRURJLCBNYXJjYSwgUHJlc2VudGFjaW9uKQ0KDQojIEFncmVnYSB1biBjYW1wbyBjYWxjdWxhZG8gY29uIGxhcyB2ZW50YXMgZGVsIHByaW1lciBzZW1lc3RyZSB5IG11ZXN0cmEgbGFzIHZlbnRhcyBkZWwgcHJpbWVyIHNlbWVzdHJlIHBvciBtYXJjYToNCmRmNSA8LSBkZjEgJT4lIA0KICBtdXRhdGUoVmVudGFzX1NlbTEgPSBFbmVybyArIEZlYnJlcm8gKyBNYXJ6byArIEFicmlsICsgTWF5byArIEp1bmlvKQ0KDQp2ZW50YXNfc2VtMV9wb3JfbWFyY2EgPC0gZGY1ICU+JSANCiAgZ3JvdXBfYnkoTWFyY2EpICU+JSANCiAgc3VtbWFyaXNlKFZlbnRhc19TZW0xPXN1bShWZW50YXNfU2VtMSkpDQoNCiMgT2J0w6luIGxhIG1lZGlhIGRlIGxhcyB2ZW50YXMgZGVsIHByaW1lciBzZW1lc3RyZSBhZ3J1cGFkbyBwb3IgbWFyY2EsIHByZXNlbnRhY2nDs24geSB0YW1hw7FvLg0KZGY2IDwtIGRmNSAlPiUNCiAgZ3JvdXBfYnkoTWFyY2EsIFByZXNlbnRhY2lvbiwgVGFtYcOxbykgJT4lDQogIHN1bW1hcmlzZShWZW50YXNfU2VtMT1tZWFuKFZlbnRhc19TZW0xKSkNCg0KIyBDYWxjdWxhciBtZWRpZGFzIGRlIHRlbmRlbmNpYSBjZW50cmFsDQpzdW1tYXJ5KGRmMSkNCg0KIyBDb2xhcHNhciBtZXNlcyBlbiB1bmEgY29sdW1uYQ0KZGY3IDwtIGdhdGhlcihkZjEsIE1lcywgVmVudGFzLCBFbmVybzpEaWNpZW1icmUpDQpkZjcNCg0KIyBBZ3JlZ2FyIHVuYSBjb2x1bW5hIGNvbiBlbCBOw7ptZXJvIGRlIE1lcw0KbWVzZXMgPC0gYygiRW5lcm8iPSIwMSIsICJGZWJyZXJvIj0iMDIiLCJNYXJ6byI9IjAzIiwgIkFicmlsIj0iMDQiLCAiTWF5byI9IjA1IiwgIkp1bmlvIj0iMDYiLCAiSnVsaW8iPSIwNyIsICJBZ29zdG8iPSIwOCIsICJTZXB0aWVtYnJlIj0iMDkiLCAiT2N0dWJyZSI9IjEwIiwgIk5vdmllbWJyZSI9IjExIiwgIkRpY2llbWJyZSI9IjEyIikNCg0KZGY3JE7Dum1lcm9fZGVfTWVzIDwtIG1lc2VzW2RmNyRNZXNdDQoNCiMgR3JhZmljYXIgZWwgdG90YWwgZGUgdmVudGFzIHBvciBtZXMgeSBhw7FvDQoNCnZlbnRhcy50b3RhbGVzIDwtIGRmNyAlPiUNCiAgZ3JvdXBfYnkoQcOxbywgTsO6bWVyb19kZV9NZXMpICU+JQ0KICBzdW1tYXJpc2UoVmVudGFzLlRvdGFsZXMgPSBzdW0oVmVudGFzKSkgJT4lDQogIGFycmFuZ2UoTsO6bWVyb19kZV9NZXMpDQoNCmdncGxvdCh2ZW50YXMudG90YWxlcywgYWVzKHg9TsO6bWVyb19kZV9NZXMsIHk9VmVudGFzLlRvdGFsZXMsZ3JvdXA9QcOxbywgY29sb3I9IGFzLmZhY3RvcihBw7FvKSkpICsNCiAgZ2VvbV9saW5lKCkgKw0KICBnZW9tX3BvaW50KCkgKw0KICBsYWJzKHg9Ik1lcyIseT0iVmVudGFzIFRvdGFsZXMoUXR5KSIsIGNvbG9yID0gIkHDsW8iKQ0KYGBgDQoNCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWU7Ij5SZWdyZXNpw7NuIExpbmVhbDwvc3Bhbj4NCmBgYHtyfQ0KIyBFbGFib3JhciBSZWdyZXNpw7NuIExpbmVhbA0KDQp2ZW50YXNfNDU2NyA8LSBkZjcgJT4lDQogIGZpbHRlcihDbGllbnRlID09ICI0NTY3IikgJT4lDQogIGdyb3VwX2J5KEHDsW8pICU+JQ0KICBzdW1tYXJpc2UoVmVudGFzX1RvdGFsZXMgPSBzdW0oVmVudGFzKSkgJT4lDQogIGZpbHRlcihWZW50YXNfVG90YWxlcyAhPSAwKQ0KDQojIHZlbnRhc180NTY3JFNlY3VlbmNpYSA8LSAxOm5yb3codmVudGFzXzQ1NjcpDQoNCmdncGxvdCh2ZW50YXNfNDU2NywgYWVzKHg9QcOxbywgeT1WZW50YXNfVG90YWxlcykpICsNCiAgZ2VvbV9wb2ludCgpICsNCiAgbGFicyh4PSJBw7FvIix5PSJWZW50YXMgVG90YWxlcyhRdHkpIiwgdGl0bGU9ICJWZW50YXMgVG90YWxlcyBkZWwgQ2xpZW50ZSA0NTY3IikNCg0KcmVncmVzaW9uIDwtIGxtKFZlbnRhc19Ub3RhbGVzIH4gQcOxbywgZGF0YT12ZW50YXNfNDU2NykNCnN1bW1hcnkocmVncmVzaW9uKQ0KDQojIEVjdWFjacOzbg0KIyB5ID0gNDg5MTg2NzcgLSAyNDE3NCAqIEHDsW8gDQoNCiMgUiBjdWFkcmFkYSBhanVzdGFkYQ0KIyA3NyUNCg0KZGF0b3MgPC0gZGF0YS5mcmFtZShBw7FvPTIwMjA6MjAyNSkNCnByZWRpY2Npb24gPC0gcHJlZGljdChyZWdyZXNpb24sZGF0b3MpDQpwcmVkaWNjaW9uDQoNCnZlbnRhc180NTY3JFRpcG9fZGVfRGF0byA8LSAiRGF0b3MgUmVhbGVzIg0KZGF0b3MkVmVudGFzX1RvdGFsZXMgPC0gcHJlZGljY2lvbg0KZGF0b3MkVGlwb19kZV9EYXRvIDwtICJQcmVkaWNjacOzbiINCg0KZGF0b3NfY29tYmluYWRvcyA8LSByYmluZCh2ZW50YXNfNDU2NyxkYXRvcykNCg0KZ2dwbG90KGRhdG9zX2NvbWJpbmFkb3MsIGFlcyh4PUHDsW8sIHk9VmVudGFzX1RvdGFsZXMsIGNvbG9yPVRpcG9fZGVfRGF0bykpICsNCiAgZ2VvbV9wb2ludCgpICsNCiAgbGFicyh4PSJBw7FvIix5PSJWZW50YXMgVG90YWxlcyhRdHkpIiwgdGl0bGU9ICJQcm9uw7NzdGljbyBhIDUgYcOxb3MgZGUgVmVudGFzIGRlbCBDbGllbnRlIDQ1NjciKQ0KYGBgDQoNCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWU7Ij5TZXJpZXMgZGUgdGllbXBvPC9zcGFuPg0KYGBge3IgbWVzc2FnZT1GQUxTRSwgd2FybmluZz1GQUxTRX0NCg0KdmVudGFzX21lbnN1YWxlc180NTY3IDwtIGRmNyAlPiUNCiAgZmlsdGVyKENsaWVudGUgPT0gIjQ1NjciKSAlPiUNCiAgZ3JvdXBfYnkoQcOxbywgTsO6bWVyb19kZV9NZXMpICU+JQ0KICBzdW1tYXJpc2UoVmVudGFzX1RvdGFsZXMgPSBzdW0oVmVudGFzKSkgJT4lDQogIGZpbHRlcihWZW50YXNfVG90YWxlcyAhPSAwKQ0KDQojIENvbmZpcm1hciBxdWUgbG9zIGRhdG9zIHF1ZSBxdWVyZW1vcyBtb2RlbGFyIGVzdGVuIG9yZGVuYWRvcyBjcm9ub2zDs2dpY2FtZW50ZS4NCg0KIyBGdW5jacOzbiBkZSBTZXJpZSBkZSBUaWVtcG8gTUVOU1VBTCwgcXVlIGluaWNpYSBlbiBFbmVybyAyMDE3DQp0cyA8LSB0cyhkYXRhPXZlbnRhc19tZW5zdWFsZXNfNDU2NyRWZW50YXNfVG90YWxlcywgc3RhcnQgPSBjKDIwMTcsMSksIGZyZXF1ZW5jeSA9IDEyKQ0KDQojIEZ1bmNpw7NuIGRlIFNlcmllIGRlIFRpZW1wbyBNRU5TVUFMLCBxdWUgaW5pY2lhIGVuIEFicmlsIDIwMTcNCiMgdHMgPC0gdHMoZGF0YT12ZW50YXNfbWVuc3VhbGVzXzQ1NjckVmVudGFzX1RvdGFsZXMsIHN0YXJ0ID0gYygyMDE3LDQpLCBmcmVxdWVuY3kgPSAxMikNCg0KIyBGdW5jacOzbiBkZSBTZXJpZSBkZSBUaWVtcG8gVFJJTUVTVFJBTCwgcXVlIGluaWNpYSBlbiBFbmVybyAyMDE3DQojIHRzIDwtIHRzKGRhdGE9dmVudGFzX21lbnN1YWxlc180NTY3JFZlbnRhc19Ub3RhbGVzLCBzdGFydCA9IGMoMjAxNywxKSwgZnJlcXVlbmN5ID0gNCkNCg0KIyBGdW5jacOzbiBkZSBTZXJpZSBkZSBUaWVtcG8gVFJJTUVTVFJBTCwgcXVlIGluaWNpYSBlbiBPY3R1YnJlIDIwMTcgKFE0KQ0KIyB0cyA8LSB0cyhkYXRhPXZlbnRhc19tZW5zdWFsZXNfNDU2NyRWZW50YXNfVG90YWxlcywgc3RhcnQgPSBjKDIwMTcsNCksIGZyZXF1ZW5jeSA9IDQpDQoNCiMgRnVuY2nDs24gZGUgU2VyaWUgZGUgVGllbXBvIEFOVUFMLCBxdWUgaW5pY2lhIGVuIDIwMTcNCiMgdHMgPC0gdHMoZGF0YT12ZW50YXNfbWVuc3VhbGVzXzQ1NjckVmVudGFzX1RvdGFsZXMsIHN0YXJ0ID0gMjAxNywgZnJlcXVlbmN5ID0gMSkNCg0KIyBDcmVhciBNb2RlbG8gQVJJTUENCiMgTW9kZWxvIEF1dG9ycmVncmVzaXZvIEludGVncmFkbyBkZSBQcm9tZWRpbyBNb3ZpbC4NCmFyaW1hIDwtIGF1dG8uYXJpbWEodHMsIEQ9MSkgIyBEPTEgcG9yIGxhIHRlbXBvcmFsaWRhZA0KYXJpbWENCnN1bW1hcnkoYXJpbWEpDQoNCiMgR2VuZXJhciBlbCBwcm9uw7NzdGljbyBkZSB2ZW50YXMNCnByb25vc3RpY28gPC0gZm9yZWNhc3QoYXJpbWEsIGxldmVsPTk1LCBoPTM2KQ0KcHJvbm9zdGljbw0KcGxvdChwcm9ub3N0aWNvLCBtYWluPSJWZW50YXMgTWVuc3VhbGVzIHkgUHJvbsOzc3RpY28gYSAzIEHDsW9zIGRlbCBDbGllbnRlIDQ1NjciLCB4bGFiPSJBw7FvIiwgeWxhYj0iVmVudGFzIChRdHkpIikNCmBgYA0KDQojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+UHJlZ3VudGFzIGRlIEV2aWRlbmNpYTwvc3Bhbj4NCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWU7Ij5QcmVndW50YSAxPC9zcGFuPg0KIyMgUHJlZ3VudGFzIERldG9uYW50ZXMNCg0KMS4gwr9QdWVkZSBvYnNlcnZhcnNlIHVuIGNyZWNpbWllbnRvIGVuIGxhcyB2ZW50YXMgZGUgYWxndW5vcyBkZSBsb3Mgc2VnbWVudG9zIGRlIHByb2R1Y3RvcyBkZSBsYSBmYW1pbGlhIENvY2EgQ29sYSBlbiBsYXMgdGllbmRhcyBlbiBsYXMgcXVlIHNlIGltcGxlbWVudMOzIGVsIFByb3llY3RvIFNpZ2xvIFhYSSBkZSBBcmNhIENvbnRpbmVudGFsPw0KYGBge3J9DQp2ZW50YXNfdG90YWxlc19zZWcgPC0gZGY3ICU+JQ0KICAgZmlsdGVyKFNlZ21lbnRvLkRldCA9PSAiQWd1YSBQdXJpZmljYWRhIiB8IFNlZ21lbnRvLkRldCA9PSAiSXNvdMOzbmljb3MgUmVndWxhciIgfCBTZWdtZW50by5EZXQgPT0gIkNvbGFzIFJlZ3VsYXIiICkgJT4lDQogIGZpbHRlcihBw7FvID09IDIwMTgpICU+JQ0KICBncm91cF9ieShTZWdtZW50by5EZXQsIE7Dum1lcm9fZGVfTWVzKSAlPiUNCiAgc3VtbWFyaXNlKFZlbnRhc19Ub3RhbGVzID0gc3VtKFZlbnRhcykpICU+JQ0KICBhcnJhbmdlKE7Dum1lcm9fZGVfTWVzKQ0KDQpnZ3Bsb3QodmVudGFzX3RvdGFsZXNfc2VnLCBhZXMoeD1Ow7ptZXJvX2RlX01lcywgeT1WZW50YXNfVG90YWxlcyxncm91cD1TZWdtZW50by5EZXQsIGNvbG9yPSBhcy5mYWN0b3IoU2VnbWVudG8uRGV0KSkpICsNCiAgZ2VvbV9saW5lKCkgKw0KICBnZW9tX3BvaW50KCkgKw0KICBsYWJzKHg9Ik1lcyIseT0iVmVudGFzIFRvdGFsZXMoUXR5KSIsIHRpdGxlPSAiVmVudGFzIFRvdGFsZXMgcG9yIFNlZ21lbnRvIiwgY29sb3IgPSAiU2VnbWVudG8iKQ0KYGBgDQoqKlJlc3B1ZXN0YTogRW4gZ2VuZXJhbCBsYXMgdmVudGFzIGF1bWVudGFyb24gZW4gMjAxOSwgcGVybyBjYWRhIHNlZ21lbnRvIHNlIGNvbXBvcnTDsyBkaWZlcmVudGUuIEVsIGluY3JlbWVudG8gZGUgdmVudGFzIHNlIHB1ZWRlIGxvZ3JhciBncmFjaWFzIGEgbWVqb3JlcyBjYW1wYcOxYXMgZGUgbWVyY2Fkb3RlY25pYSB5IGxhbnphbWllbnRvIGRlIG51ZXZvcyBwcm9kdWN0b3MuRW4gMjAxOCBsYSB2ZW50YSBkZSBDb2NhLUNvbGEgZW50cmUgRmVicmVybyB5IERpY2llbWJyZSBhdW1lbnTDsyBlbiAxMDAsMDAwIHVuaWRhZGVzLCBtaWVudHJhcyBxdWUgbGFzIGFndWFzIGlzb3TDs25pY2FzIG5vIHR1dmllcm9uIGNhbWJpb3MuKioNCg0KIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZTsiPlByZWd1bnRhIDI8L3NwYW4+DQoyLiDCv0VsIGluY3JlbWVudG8gZW4gbGFzIHZlbnRhcyBlcyBzaW1pbGFyIGVudHJlIGxvcyBkaWZlcmVudGVzIHRhbWHDsW9zIGRlIGNsaWVudGVzPw0KYGBge3J9DQp2ZW50YXNfdG90YWxlc190YW1fY2xpZW50ZSA8LSBkZjcgJT4lDQogIGZpbHRlcihBw7FvID09IDIwMTkpICU+JQ0KICBncm91cF9ieShUYW1hw7FvLkN0ZS5JbmR1c3RyaWEsIE7Dum1lcm9fZGVfTWVzKSAlPiUNCiAgc3VtbWFyaXNlKFZlbnRhc19Ub3RhbGVzID0gc3VtKFZlbnRhcykpICU+JQ0KICBhcnJhbmdlKE7Dum1lcm9fZGVfTWVzKSAlPiUNCiAgZmlsdGVyKFZlbnRhc19Ub3RhbGVzICE9IDApDQoNCmdncGxvdCh2ZW50YXNfdG90YWxlc190YW1fY2xpZW50ZSwgYWVzKHg9TsO6bWVyb19kZV9NZXMsIHk9VmVudGFzX1RvdGFsZXMsZ3JvdXA9VGFtYcOxby5DdGUuSW5kdXN0cmlhLCBjb2xvcj0gYXMuZmFjdG9yKFRhbWHDsW8uQ3RlLkluZHVzdHJpYSkpKSArDQogIGdlb21fbGluZSgpICsNCiAgZ2VvbV9wb2ludCgpICsNCiAgbGFicyh4PSJNZXMiLHk9IlZlbnRhcyBUb3RhbGVzKFF0eSkiLCB0aXRsZT0gIlZlbnRhcyBUb3RhbGVzIHBvciBUYW1hw7FvIGRlIENsaWVudGUiLCBjb2xvciA9ICJUYW1hw7FvIGRlIENsaWVudGUiKQ0KYGBgDQoNCioqUmVzcHVlc3RhOiBFbCBpbmNyZW1lbnRvIGRlIHZlbnRhcyBlcyBzaW1pbGFyIGluZGVwZW5kaWVudGVtZW50ZSBkZSBsb3MgdGFtYcOxb3MgZGUgbG9zIGNsaWVudGVzLiBDb21vIHRvZG9zIGxvcyBuZWdvY2lvcyBlc3TDoW4gYWJpZXJ0b3MgbG9zIG1pc21vcyBkw61hcywgeSBlc3TDoW4gdWJpY2Fkb3MgZW4gbGEgbWlzbWEgY2l1ZGFkIChHdWFkYWxhamFyYSksIGxvcyBmYWN0b3JlcyBleHRlcm5vcyBkZSBsYXMgdmVudGFzIGVmZWN0YW4gY2FzaSBpZ3VhbCBhIHRvZG9zLiBFbCBtZXMgZGUgbWF5byBkZSAyMDE5IGZ1ZSBlbCBkZSBtYXlvciB2ZW50YXMgZW4gdG9kb3MgbG9zIHRhbWHDsW9zIGRlIGNsaWVudGUsIHkgZWwgcHJpbWVyIGJpbWVzdHJlIGZ1ZSBlbCBtw6FzIGJham8gdGFtYmnDqW4gcGFyYSB0b2RvcyBsb3MgdGFtYcOxb3MgZGUgY2xpZW50ZS4qKg0KDQojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+UHJlZ3VudGEgMzwvc3Bhbj4NCjMuIMK/Q3XDoWwgZXMgZWwgY29tcG9ydGFtaWVudG8gb2JzZXJ2YWRvIGRlIGxhcyB1bmlkYWRlcyB2ZW5kaWRhcyBwb3IgbWVzIGRlIGNhZGEgdW5hIGRlIGxhcyBtYXJjYXMsIGluZGVwZW5kaWVudGVtZW50ZSBkZSBzdXMgcmVzcGVjdGl2YXMgcHJlc2VudGFjaW9uZXM/DQpgYGB7cn0NCnZlbnRhc190b3RhbGVzX21hcmNhIDwtIGRmNyAlPiUNCiAgZmlsdGVyKEHDsW8gPT0gMjAxOSkgJT4lDQogIGZpbHRlcihNYXJjYSA9PSBjKCJDb2NhLUNvbGEiLCAiVmFsbGUgRnJ1dCIsICJDaWVsIEFndWEgUHVyaWZpY2FkYSIsIkNvY2EtQ29sYSBMaWdodCIpKSAlPiUNCiAgZ3JvdXBfYnkoTWFyY2EsIE7Dum1lcm9fZGVfTWVzKSAlPiUNCiAgc3VtbWFyaXNlKFZlbnRhc19Ub3RhbGVzID0gc3VtKFZlbnRhcykpICU+JQ0KICBhcnJhbmdlKE7Dum1lcm9fZGVfTWVzKSAlPiUNCiAgZmlsdGVyKFZlbnRhc19Ub3RhbGVzICE9IDApDQoNCmdncGxvdCh2ZW50YXNfdG90YWxlc19tYXJjYSwgYWVzKHg9TsO6bWVyb19kZV9NZXMsIHk9VmVudGFzX1RvdGFsZXMsZ3JvdXA9TWFyY2EsIGNvbG9yPSBhcy5mYWN0b3IoTWFyY2EpKSkgKw0KICBnZW9tX2xpbmUoKSArDQogIGdlb21fcG9pbnQoKSArDQogIGxhYnMoeD0iTWVzIix5PSJWZW50YXMgVG90YWxlcyhRdHkpIiwgdGl0bGU9ICJWZW50YXMgVG90YWxlcyBwb3IgTWFyY2EiLCBjb2xvciA9ICJNYXJjYSIpDQpgYGANCg0KKipSZXNwdWVzdGE6IEVsIGNvbXBvcnRhbWllbnRvIGRlIHZlbnRhcyBlcyBzaW1pbGFyIGVudHJlIGxhcyBtYXJjYXMgYSB0cmF2w6lzIGRlIGxvcyBtZXNlcywgcGVybyBoYXkgbWFyY2FzIHF1ZSB2ZW5kZW4gbXVjaG8gbcOhcyBxdWUgb3RyYXMuIExvcyBtZXNlcyBkZSBtw6FzIHZlbnRhcyBkZSBiZWJpZGFzIHNvbiBsb3MgbWVzZXMgZW4gcXVlIGxhIHRlbXBlcmF0dXJhIGVzIG1heW9yLCB5IGxhcyBtYXJjYXMgbcOhcyBjb25vY2lkYXMgc29uIGxhcyBxdWUgdGllbmUgZW4gbWVudGUgZWwgY29uc3VtaWRvci4gRWwgbWVzIGRlIG1heW8gMjAxOSBlcyBwYXJhIHRvZGFzIGxhcyBtYXJjYXMgZWwgZGUgbWF5b3JlcyB2ZW50YXMsIHBlcm8gbGEgbWFyY2EgQ29jYS1Db2xhIGVzIGxhIHF1ZSB2ZW5kZSBtw6FzIHF1ZSB0b2Rhcy4qKg0KDQojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+UHJlZ3VudGEgNDwvc3Bhbj4NCjQuIMK/U2UgaGEgaW5jcmVtZW50YWRvIGxhIHZlbnRhIGRlIHByb2R1Y3RvcyBlbiBlbnZhc2VzIHJldG9ybmFibGVzIGVuIGxvcyDDumx0aW1vcyBkb3MgYcOxb3M/DQpgYGB7cn0NCnZlbnRhc190b3RhbGVzX2VudmFzZSA8LSBkZjcgJT4lDQogIGZpbHRlcihSZXRvcm5hYmxlX05SPT0iUmV0b3JuYWJsZSIpICU+JQ0KICBmaWx0ZXIoQcOxbyA9PSBjKDIwMTcsMjAxOCwyMDE5KSkgJT4lDQogIGdyb3VwX2J5KEHDsW8sIE7Dum1lcm9fZGVfTWVzKSAlPiUNCiAgc3VtbWFyaXNlKFZlbnRhc19Ub3RhbGVzID0gc3VtKFZlbnRhcykpICU+JQ0KICBhcnJhbmdlKE7Dum1lcm9fZGVfTWVzKSAlPiUNCiAgZmlsdGVyKFZlbnRhc19Ub3RhbGVzICE9IDApDQoNCmdncGxvdCh2ZW50YXNfdG90YWxlc19lbnZhc2UsIGFlcyh4PU7Dum1lcm9fZGVfTWVzLCB5PVZlbnRhc19Ub3RhbGVzLGdyb3VwPUHDsW8sIGNvbG9yPSBhcy5mYWN0b3IoQcOxbykpKSArDQogIGdlb21fbGluZSgpICsNCiAgZ2VvbV9wb2ludCgpICsNCiAgbGFicyh4PSJNZXMiLHk9IlZlbnRhcyBUb3RhbGVzKFF0eSkiLCB0aXRsZT0gIlZlbnRhcyBUb3RhbGVzIGRlIEVudmFzZXMgUmV0b3JuYWJsZXMiLCBjb2xvciA9ICJBw7FvIikNCmBgYA0KDQoqKlJlc3B1ZXN0YTogTGEgdmVudGEgZGUgcHJvZHVjdG9zIGVuIGVudmFzZSByZXRvcm5hYmxlIGF1bWVudMOzIGVuIGxvcyDDumx0aW1vcyBkb3MgYcOxb3MuIExhcyBjYW1wYcOxYXMgZGUgY29uY2llbmNpYSBhbWJpZW50YWwgaGFuIHZ1ZWx0byBwb3B1bGFyZXMgbGFzIG9wY2lvbmVzIG1lbm9zIGRhw7FpbmFzIGFsIGVjb3Npc3RlbWEuIERlIDIwMTcgYSAyMDE4IGh1Ym8gdW4gbGlnZXJvIGluY3JlbWVudG8gZW4gbGFzIHZlbnRhcyBkZSBlbnZhc2UgcmV0b3JuYWJsZSwgeSBlbnRyZSAyMDE4IHkgMjAxOSBlbCBhdW1lbnRvIGZ1ZSBtdXkgbm90b3Jpby4qKiAgDQoNCioqTm90YXM6IEhheSBvdHJvcyBmYWN0b3JlcyBhIGNvbnNpZGVyYXIgY29tbyBxdWUgZW4gMjAxOSBoYXkgbcOhcyB0aWVuZGFzIGNvbiBlbCBwcm95ZWN0byBTaWdsbyBYWEksIG8gcXVlIGxhcyB0aWVuZGFzIGVuIGdlbmVyYWwgdmVuZGVuIG3DoXMgcHJvZHVjdG9zICh0YW50byByZXRvcm5hYmxlcyBjb21vIG5vIHJldG9ybmFibGVzKS4qKg0KDQojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+UHJlZ3VudGEgNTwvc3Bhbj4NCjUuIMK/RWwgY29tcG9ydGFtaWVudG8gZGUgbGEgdmVudGEgZGUgYWd1YSBoYSBpbmNyZW1lbnRhZG8gZW4gcmVsYWNpw7NuIGFsIGRlIGxvcyByZWZyZXNjb3MgbyBsYXMgYmViaWRhcyBpc290w7NuaWNhcz8NCmBgYHtyfQ0KdmVudGFzX3RvdGFsZXNfc2VnIDwtIGRmNyAlPiUNCiAgIGZpbHRlcihTZWdtZW50by5EZXQgPT0gIkFndWEgUHVyaWZpY2FkYSIgfCBTZWdtZW50by5EZXQgPT0gIklzb3TDs25pY29zIFJlZ3VsYXIiIHwgU2VnbWVudG8uRGV0ID09ICJDb2xhcyBSZWd1bGFyIiApICU+JQ0KICBmaWx0ZXIoQcOxbyA9PSAyMDE4KSAlPiUNCiAgZ3JvdXBfYnkoU2VnbWVudG8uRGV0LCBOw7ptZXJvX2RlX01lcykgJT4lDQogIHN1bW1hcmlzZShWZW50YXNfVG90YWxlcyA9IHN1bShWZW50YXMpKSAlPiUNCiAgYXJyYW5nZShOw7ptZXJvX2RlX01lcykNCg0KZ2dwbG90KHZlbnRhc190b3RhbGVzX3NlZywgYWVzKHg9TsO6bWVyb19kZV9NZXMsIHk9VmVudGFzX1RvdGFsZXMsZ3JvdXA9U2VnbWVudG8uRGV0LCBjb2xvcj0gYXMuZmFjdG9yKFNlZ21lbnRvLkRldCkpKSArDQogIGdlb21fbGluZSgpICsNCiAgZ2VvbV9wb2ludCgpICsNCiAgbGFicyh4PSJNZXMiLHk9IlZlbnRhcyBUb3RhbGVzKFF0eSkiLCB0aXRsZT0gIlZlbnRhcyBUb3RhbGVzIHBvciBTZWdtZW50byIsIGNvbG9yID0gIlNlZ21lbnRvIikNCmBgYA0KDQoqKlJlc3B1ZXN0YTogTGFzIHZlbnRhcyBkZSBBZ3VhIFB1cmlmaWNhZGEgbm8gaGEgYXVtZW50YWRvIGVuIHJlbGFjacOzbiBhIFJlZnJlc2NvcyBvIEJlYmlkYXMgSXNvdMOzbmljYXMuIExhIHBlcnNlcGNpw7NuIHF1ZSBlbCBhZ3VhIG5vIHRpZW5lIHZhbG9yIGFncmVnYWRvIGVuIGxhIHByZXBhcmFjacOzbiB5IGVsIGF6w7pjYXIgZW4gbGEgcmVjZXRhIGhhY2VuIGRlIGxvcyByZWZyZXNjb3MgbGEgb3BjacOzbiBwcmVmZXJpZGEgZGUgbG9zIGNvbnN1bWlkb3JlcyBtZXhpY2Fub3MuIEVuIDIwMTggbGFzIHZlbnRhcyB0b3RhbGVzIGRlIEFndWEgUHVyaWZpY2FkYSBubyBzdXBlcmFyb24gbGFzIDUwLDAwMCB1bmlkYWRlcywgbWllbnRyYXMgcXVlIGxhIENvY2EtQ29sYSB2ZW5kacOzIGVuIHN1IG1lcyBtw6FzIGFsdG8gNDUwLDAwMCAoOCB2ZWNlcyBtw6FzIHF1ZSBlbCBhZ3VhKS4qKg0KDQojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+UHJlZ3VudGEgNjwvc3Bhbj4NCjYuIMK/UHVlZGUgZGVjaXJzZSBxdWUgbGEgdmVudGEgbWVuc3VhbCBkZSBhZ3VhIGVzdMOhIHJlbGFjaW9uYWRhIGNvbiBsYSB2ZW50YSBtZW5zdWFsIGRlIHJlZnJlc2NvcyBlbiBsb3Mgw7psdGltb3MgNCBhw7Fvcz8NCmBgYHtyfQ0KdmVudGFzX3RvdGFsZXNfYXlyIDwtIGRmNyAlPiUNCiAgIGZpbHRlcihTZWdtZW50by5EZXQgPT0gYygiQWd1YSBQdXJpZmljYWRhIiwgIkNvbGFzIFJlZ3VsYXIiKSkgJT4lDQogIGdyb3VwX2J5KFNlZ21lbnRvLkRldCwgQcOxbykgJT4lDQogIHN1bW1hcmlzZShWZW50YXNfVG90YWxlcyA9IHN1bShWZW50YXMpKSANCiAgDQpnZ3Bsb3QodmVudGFzX3RvdGFsZXNfYXlyLCBhZXMoeD1Bw7FvLCB5PVZlbnRhc19Ub3RhbGVzLGdyb3VwPVNlZ21lbnRvLkRldCwgY29sb3I9IGFzLmZhY3RvcihTZWdtZW50by5EZXQpKSkgKw0KICBnZW9tX2xpbmUoKSArDQogIGdlb21fcG9pbnQoKSArDQogIGxhYnMoeD0iQcOxbyIseT0iVmVudGFzIFRvdGFsZXMoUXR5KSIsIHRpdGxlPSAiVmVudGFzIFRvdGFsZXMgcG9yIFNlZ21lbnRvIiwgY29sb3IgPSAiU2VnbWVudG8iKQ0KYGBgDQoNCioqUmVzcHVlc3RhOiBMYSB2ZW50YSBkZSBhZ3VhIHkgcmVmcmVzY29zIHNlIGluY3JlbWVudGEgZW4gbG9zIG1pc21vcyBtZXNlcywgeSBhIHRyYXbDqXMgZGUgbG9zIGHDsW9zLCBwZXJvIG5vIHNlIHB1ZWRlIGNvbmNsdWlyIHF1ZSB1biBpbmNyZW1lbnRvIGNhdXNlIGVsIG90cm8uIExvcyBtZXNlcyBjYWx1cm9zb3MgeSBlbCBjcmVjaW1pZW50byBkZSBsYSBwb2JsYWNpw7NuIGhhY2UgcXVlIGxhcyB2ZW50YXMgc2UgaW5jcmVtZW50ZW4gZW4gbGFzIGJlYmlkYXMsIHNlYW4gYWd1YSBvIHJlZnJlc2NvLiBFbiBsYSBncsOhZmljYSBvYnRlbmlkYSBzZSBvYnNlcnZhIHVuIGluY3JlbWVudG8gYcOxbyBjb24gYcOxbyBlbiBsYSB2ZW50YSBkZSBhZ3VhIHkgdGFtYmnDqW4gZGUgcmVmcmVzY29zLiAqKg0KDQojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+UHJlZ3VudGEgNzwvc3Bhbj4NCjcuIMK/QSBjdcOhbnRvIGFzY2llbmRlbiBsYXMgdmVudGFzIGVzcGVyYWRhcyBwYXJhIGVsIDIwMjAgZW4gbGEgQ29jYSBDb2xhIGRlIDUwMCBtbCBOUiBWaWRyaW8/DQpgYGB7cn0NCnZlbnRhc19jb2NhX3ZpZHJpbyA8LSBkZjcgJT4lDQogIGZpbHRlcihNYXJjYSA9PSAiQ29jYS1Db2xhIiAmIFByZXNlbnRhY2lvbiA9PSAiNTAwIG1sIE5SIFZpZHJpbyIpICU+JQ0KICBncm91cF9ieShBw7FvLCBOw7ptZXJvX2RlX01lcykgJT4lDQogIHN1bW1hcmlzZShWZW50YXNfVG90YWxlcyA9IHN1bShWZW50YXMpKSAlPiUNCiAgZmlsdGVyKFZlbnRhc19Ub3RhbGVzICE9IDApDQoNCiMgQ29uZmlybWFyIHF1ZSBsb3MgZGF0b3MgcXVlIHF1ZXJlbW9zIG1vZGVsYXIgZXN0ZW4gb3JkZW5hZG9zIGNyb25vbMOzZ2ljYW1lbnRlLg0KDQojIEZ1bmNpw7NuIGRlIFNlcmllIGRlIFRpZW1wbyBNRU5TVUFMLCBxdWUgaW5pY2lhIGVuIEVuZXJvIDIwMTcNCnRzIDwtIHRzKGRhdGE9dmVudGFzX2NvY2FfdmlkcmlvJFZlbnRhc19Ub3RhbGVzLCBzdGFydCA9IGMoMjAxNiwxKSwgZnJlcXVlbmN5ID0gMTIpDQoNCiMgQ3JlYXIgTW9kZWxvIEFSSU1BDQojIE1vZGVsbyBBdXRvcnJlZ3Jlc2l2byBJbnRlZ3JhZG8gZGUgUHJvbWVkaW8gTW92aWwuDQphcmltYSA8LSBhdXRvLmFyaW1hKHRzLCBEPTEpICMgRD0xIHBvciBsYSB0ZW1wb3JhbGlkYWQNCmFyaW1hDQpzdW1tYXJ5KGFyaW1hKQ0KDQojIEdlbmVyYXIgZWwgcHJvbsOzc3RpY28gZGUgdmVudGFzDQpwcm9ub3N0aWNvIDwtIGZvcmVjYXN0KGFyaW1hLCBsZXZlbD05NSwgaD0xNikNCnByb25vc3RpY28NCnBsb3QocHJvbm9zdGljbywgbWFpbj0iVmVudGFzIE1lbnN1YWxlcyB5IFByb27Ds3N0aWNvIHBhcmEgMjAyMCBkZSBDb2NhLUNvbGEgNTAwIG1sIE5SIFZpZHJpbyIsIHhsYWI9IkHDsW8iLCAgeWxhYj0iVmVudGFzIChRdHkpIikNCmBgYA0KDQoqKlJlc3B1ZXN0YTogTGFzIHZlbnRhcyBlc3BlcmFkYXMgZGUgQ29jYS1Db2xhIGRlIDUwMCBtbCBOUiBWaWRyaW8gcGFyYSAyMDIwIG9zY2lsYW4gbWVuc3VhbG1lbnRlIGVudHJlIDQ1LDAwMCB5IDY1LDAwMCB1bmlkYWRlcywgY29uIHVuYSBjb25maWFiaWxpZGFkIGRlbCA5NSUuKioNCg==