De cada uno de los siguientes ejercicios, determinar cómo estará la
subpoblación en 50 generaciones. Además, deben realizar un gráfico para
visualizar el comportamiento de la población de cada organismo.
Ejercicio 1
En un parche de bosque húmedo de Sarapiquí crece una población de 300
plantas de Heliconia latispatha. Las brácteas son hojas
modificadas de color vistoso que rodean las pequeñas flores verdaderas
de la inflorescencia. En esta población, el 75 % de las plantas tiene
brácteas rojas, con genotipos RR o Rr, por lo que el alelo R es
dominante. El 25 % restante tiene brácteas amarillas, con genotipo rr, y
el alelo r es recesivo.
La reproducción de estas plantas depende de la polinización. Los
colibríes visitan las inflorescencias para alimentarse del néctar, y el
color de las brácteas es la principal señal que los atrae. Cuando un
colibrí se posa y luego vuela a otra planta, transporta polen de una
inflorescencia a otra, lo que permite la formación de frutos. Una
inflorescencia visitada con frecuencia produce más frutos y, por lo
tanto, más semillas que pueden dar origen a nuevas plantas.
Durante una temporada de floración se contó el número de frutos por
inflorescencia. Las plantas de brácteas rojas produjeron en promedio 14
frutos, mientras que las de brácteas amarillas produjeron 11.
Supuesto de tiempo: una generación equivale a 4
años.
FitnessAA <- 14/225
FitnessAA
[1] 0.06222222
```{r}
Error: attempt to use zero-length variable name
FitnessAA/Fitnessaa
[1] 0.4242424
return(list(c(dp, dq)))
Error: object 'dp' not found
```{r}
Error: attempt to use zero-length variable name
tail(sol)
time A a
[96,] 47.5 1 2.678154e-14
[97,] 48.0 1 2.008615e-14
[98,] 48.5 1 1.506461e-14
[99,] 49.0 1 1.129846e-14
[100,] 49.5 1 8.473845e-15
[101,] 50.0 1 6.355384e-15


Ejercicio 2
En una finca de Upala se estudian 200 árboles de cacao criollo. El
color de la mazorca está controlado por un gen con dos alelos: R, que
produce mazorca roja y es dominante, y r, que produce mazorca verde y es
recesivo. El 60 % de los árboles tiene mazorca roja, con genotipos RR o
Rr, y el 40 % tiene mazorca verde, con genotipo rr.
El cacao se reproduce mediante flores que aparecen directamente sobre
el tronco y las ramas gruesas. Estas flores son polinizadas
principalmente por pequeños insectos, como los mosquitos del
género Forcipomyia, que se desplazan entre flores buscando
alimento. Para que una flor se convierta en mazorca, el polen debe
llegar a su estigma y la planta debe ser compatible con ese polen, ya
que el cacao tiene un sistema de autoincompatibilidad que favorece el
cruce entre árboles distintos. Una flor bien polinizada se transforma en
una mazorca que contiene semillas, y las semillas son el producto que el
agricultor cosecha. El color de la mazorca depende del genotipo de la
planta que la produce, así que la cantidad de mazorcas sanas que produce
cada árbol se relaciona directamente con su éxito reproductivo.
Durante tres años se registró el promedio anual de mazorcas sanas por
árbol. Los árboles de mazorca roja produjeron 22 mazorcas por año, y los
de mazorca verde produjeron 18.
Supuesto de tiempo: una generación equivale a 10
años.
library(deSolve)
modelo <- function(time, state, parameters) {
with(as.list(c(state, parameters)), {
# Parámetros de fitness
w_AA <- 0.60 # Fitness relativo de AA
w_aa <- 0.40 # Fitness relativo de aa
# Frecuencias alélicas
p <- A # Frecuencia de A
q <- a # Frecuencia de a
# Fitness promedio poblacional
w_prom <- p^2 * w_AA + q^2 * w_aa
# Tasa de cambio en la frecuencia del alelo A
dp <- (p * q * (w_AA - w_aa)) / w_prom
dq <- -dp # La suma de frecuencias debe ser 1
return(list(c(dp, dq)))
})
}
tail(sol)
time A a
[96,] 47.5 1 3.779419e-08
[97,] 48.0 1 3.193754e-08
[98,] 48.5 1 2.737503e-08
[99,] 49.0 1 2.346431e-08
[100,] 49.5 1 2.011227e-08
[101,] 50.0 1 1.723909e-08


Ejercicio 3
En un embalse de la vertiente del Pacífico vive una población de 500
individuos de guapote. Esta especie de cíclido presenta dos morfos de
coloración que dependen de un gen con dos alelos. El morfo dominante
tiene coloración oscura con franjas laterales, con genotipos AA o Aa. El
morfo recesivo tiene coloración clara y uniforme, con genotipo aa. En la
población, el 84 % pertenece al morfo dominante y el 16 % al morfo
recesivo.
El guapote tiene cuidado parental biparental. Después de la
reproducción, la pareja construye un nido, defiende los huevos y luego
protege a las crías durante sus primeras semanas. Esta defensa reduce la
depredación, pero las crías siguen enfrentando riesgos como la
competencia por alimento, los depredadores de mayor tamaño y los cambios
en la calidad del agua. La supervivencia de las crías hasta los 6 meses
de edad determina cuántos individuos llegan a reproducirse. En este
caso, sobrevivió el 92 % de las crías del morfo dominante y el 78 % de
las del morfo recesivo.
Supuesto de tiempo: una generación equivale a 3
años.
modelo <- function(time, state, parameters) {
with(as.list(c(state, parameters)), {
# Parámetros de fitness
w_AA <- 0.84 # Fitness relativo de AA
w_aa <- 0.16 # Fitness relativo de aa
# Frecuencias alélicas
p <- A # Frecuencia de A
q <- a # Frecuencia de a
# Fitness promedio poblacional
w_prom <- p^2 * w_AA + q^2 * w_aa
# Tasa de cambio en la frecuencia del alelo A
dp <- (p * q * (w_AA - w_aa)) / w_prom
dq <- -dp # La suma de frecuencias debe ser 1
return(list(c(dp, dq)))
})
}
tail(sol)
time A a
[96,] 47.5 1 8.837337e-17
[97,] 48.0 1 6.290985e-17
[98,] 48.5 1 4.478329e-17
[99,] 49.0 1 3.187963e-17
[100,] 49.5 1 2.269397e-17
[101,] 50.0 1 1.615503e-17


FitnessAA
[1] 0.9190476
Fitnessaa
[1] 0.775
LS0tDQp0aXRsZTogIkVqZXJjaWNpb3MgQ2FtYmlvcyBHZW7DqXRpY29zIFBvYmxhY2lvbmFsZXMiDQpvdXRwdXQ6IA0KICBodG1sX25vdGVib29rOiANCiAgICB0b2M6IHRydWUNCiAgICB0aGVtZTogZmxhdGx5DQotLS0NCg0KRGUgY2FkYSB1bm8gZGUgbG9zIHNpZ3VpZW50ZXMgZWplcmNpY2lvcywgZGV0ZXJtaW5hciBjw7NtbyBlc3RhcsOhIGxhIHN1YnBvYmxhY2nDs24gZW4gNTAgZ2VuZXJhY2lvbmVzLiBBZGVtw6FzLCBkZWJlbiByZWFsaXphciB1biBncsOhZmljbyBwYXJhIHZpc3VhbGl6YXIgZWwgY29tcG9ydGFtaWVudG8gZGUgbGEgcG9ibGFjacOzbiBkZSBjYWRhIG9yZ2FuaXNtby4NCg0KIyBFamVyY2ljaW8gMQ0KDQpFbiB1biBwYXJjaGUgZGUgYm9zcXVlIGjDum1lZG8gZGUgU2FyYXBpcXXDrSBjcmVjZSB1bmEgcG9ibGFjacOzbiBkZSAzMDAgcGxhbnRhcyBkZcKgKkhlbGljb25pYSBsYXRpc3BhdGhhKi4gTGFzIGJyw6FjdGVhcyBzb24gaG9qYXMgbW9kaWZpY2FkYXMgZGUgY29sb3IgdmlzdG9zbyBxdWUgcm9kZWFuIGxhcyBwZXF1ZcOxYXMgZmxvcmVzIHZlcmRhZGVyYXMgZGUgbGEgaW5mbG9yZXNjZW5jaWEuIEVuIGVzdGEgcG9ibGFjacOzbiwgZWwgNzUgJSBkZSBsYXMgcGxhbnRhcyB0aWVuZSBicsOhY3RlYXMgcm9qYXMsIGNvbiBnZW5vdGlwb3MgUlIgbyBSciwgcG9yIGxvIHF1ZSBlbCBhbGVsbyBSIGVzIGRvbWluYW50ZS4gRWwgMjUgJSByZXN0YW50ZSB0aWVuZSBicsOhY3RlYXMgYW1hcmlsbGFzLCBjb24gZ2Vub3RpcG8gcnIsIHkgZWwgYWxlbG8gciBlcyByZWNlc2l2by4NCg0KTGEgcmVwcm9kdWNjacOzbiBkZSBlc3RhcyBwbGFudGFzIGRlcGVuZGUgZGUgbGEgcG9saW5pemFjacOzbi4gTG9zIGNvbGlicsOtZXMgdmlzaXRhbiBsYXMgaW5mbG9yZXNjZW5jaWFzIHBhcmEgYWxpbWVudGFyc2UgZGVsIG7DqWN0YXIsIHkgZWwgY29sb3IgZGUgbGFzIGJyw6FjdGVhcyBlcyBsYSBwcmluY2lwYWwgc2XDsWFsIHF1ZSBsb3MgYXRyYWUuIEN1YW5kbyB1biBjb2xpYnLDrSBzZSBwb3NhIHkgbHVlZ28gdnVlbGEgYSBvdHJhIHBsYW50YSwgdHJhbnNwb3J0YSBwb2xlbiBkZSB1bmEgaW5mbG9yZXNjZW5jaWEgYSBvdHJhLCBsbyBxdWUgcGVybWl0ZSBsYSBmb3JtYWNpw7NuIGRlIGZydXRvcy4gVW5hIGluZmxvcmVzY2VuY2lhIHZpc2l0YWRhIGNvbiBmcmVjdWVuY2lhIHByb2R1Y2UgbcOhcyBmcnV0b3MgeSwgcG9yIGxvIHRhbnRvLCBtw6FzIHNlbWlsbGFzIHF1ZSBwdWVkZW4gZGFyIG9yaWdlbiBhIG51ZXZhcyBwbGFudGFzLg0KDQpEdXJhbnRlIHVuYSB0ZW1wb3JhZGEgZGUgZmxvcmFjacOzbiBzZSBjb250w7MgZWwgbsO6bWVybyBkZSBmcnV0b3MgcG9yIGluZmxvcmVzY2VuY2lhLiBMYXMgcGxhbnRhcyBkZSBicsOhY3RlYXMgcm9qYXMgcHJvZHVqZXJvbiBlbiBwcm9tZWRpbyAxNCBmcnV0b3MsIG1pZW50cmFzIHF1ZSBsYXMgZGUgYnLDoWN0ZWFzIGFtYXJpbGxhcyBwcm9kdWplcm9uIDExLg0KDQoqKlN1cHVlc3RvIGRlIHRpZW1wbzoqKsKgdW5hIGdlbmVyYWNpw7NuIGVxdWl2YWxlIGEgNCBhw7Fvcy4NCg0KYGBge3J9DQpGaXRuZXNzQUEgPC0gMTQvMjI1IA0KYGBgDQpgYGB7cn0NCkZpdG5lc3NhYSA8LSAxMS83NQ0KYGBgDQpgYGB7cn0NCkZpdG5lc3NBQS9GaXRuZXNzYWENCmBgYA0KYGBge3J9DQpsaWJyYXJ5KGRlU29sdmUpDQoNCm1vZGVsbyA8LSBmdW5jdGlvbih0aW1lLCBzdGF0ZSwgcGFyYW1ldGVycykgew0KICB3aXRoKGFzLmxpc3QoYyhzdGF0ZSwgcGFyYW1ldGVycykpLCB7DQogICAgIyBQYXLDoW1ldHJvcyBkZSBmaXRuZXNzDQogICAgd19BQSA8LSAwLjc1ICAjIEZpdG5lc3MgcmVsYXRpdm8gZGUgQUENCiAgICB3X2FhIDwtIDAuMjUgIyBGaXRuZXNzIHJlbGF0aXZvIGRlIGFhDQogICAgDQogICAgIyBGcmVjdWVuY2lhcyBhbMOpbGljYXMNCiAgICBwIDwtIEEgICMgRnJlY3VlbmNpYSBkZSBBDQogICAgcSA8LSBhICAjIEZyZWN1ZW5jaWEgZGUgYQ0KICAgIA0KICAgICMgRml0bmVzcyBwcm9tZWRpbyBwb2JsYWNpb25hbA0KICAgIHdfcHJvbSA8LSBwXjIgKiB3X0FBICsgcV4yICogd19hYSAgDQogICAgDQogICAgIyBUYXNhIGRlIGNhbWJpbyBlbiBsYSBmcmVjdWVuY2lhIGRlbCBhbGVsbyBBDQogICAgZHAgPC0gKHAgKiBxICogKHdfQUEgLSB3X2FhKSkgLyB3X3Byb20NCiAgICBkcSA8LSAtZHAgICMgTGEgc3VtYSBkZSBmcmVjdWVuY2lhcyBkZWJlIHNlciAxDQogICAgDQogICAgcmV0dXJuKGxpc3QoYyhkcCwgZHEpKSkNCiAgfSkNCn0NCmBgYA0KYGBge3J9DQojIENvbmRpY2lvbmVzIGluaWNpYWxlcw0KZXN0YWRvX2luaWNpYWwgPC0gYyhBID0gMC43NSwgYSA9IDAuMjUpDQojIFBhcsOhbWV0cm9zDQpwYXJhbWV0cm9zIDwtIE5VTEwNCiMgVGllbXBvDQp0aWVtcG8gPC0gc2VxKDAsIDUwLCBieSA9IDAuNSkNCg0KIyMgUmVzb2x1Y2nDs24gZGVsIG1vZGVsbw0Kc29sIDwtIG9kZSh5ID0gZXN0YWRvX2luaWNpYWwsIHRpbWVzID0gdGllbXBvLCBmdW5jID0gbW9kZWxvLCBwYXJtcyA9IHBhcmFtZXRyb3MpICMjIEVzdGEgZsOzcm11bGEgbWUgcmVzdWVsdmUgbGEgZWN1YWNpw7NuIGRpZmVyZW5jaWFsDQoNCmBgYA0KYGBge3J9DQpoZWFkKHNvbCkNCg0KdGFpbChzb2wpDQpgYGANCmBgYHtyfQ0KcGxvdChzb2wsDQogICAgIHhsYWIgPSAiVGllbXBvIiwgeWxhYiA9ICJGcmVjdWVuY2lhIGRlIGFsZWxvcyIsDQogICAgIGNvbCA9IGMoImJsdWUiLCAicmVkIiksIGx3ZCA9IDIpDQoNCg0KbGlicmFyeShnZ3Bsb3QyKQ0KYGBgDQpgYGB7cn0NCmRmIDwtIGFzLmRhdGEuZnJhbWUoc29sKQ0KZGYgPC0gY2JpbmQodGllbXBvID0gdGllbXBvLCBkZikNCg0KZ2dwbG90KGRmLCBhZXMoeCA9IHRpZW1wbykpICsNCiAgZ2VvbV9saW5lKGFlcyh5ID0gbG9nKEEpLCBjb2xvciA9ICJBbGVsbyBBIikpICsgDQogIGdlb21fbGluZShhZXMoeSA9IGxvZyhhKSwgY29sb3IgPSAiQWxlbG8gYSIpKSArDQogIGxhYnMoeCA9ICJUaWVtcG8iLCB5ID0gIkZyZWN1ZW5jaWEgZGUgYWxlbG9zIikgKw0KICBzY2FsZV9jb2xvcl9tYW51YWwobmFtZSA9ICJBbGVsb3MiLCB2YWx1ZXMgPSBjKCJibHVlIiwgInJlZCIpKSArDQogIHRoZW1lX21pbmltYWwoKQ0KDQpgYGANCg0KDQojIEVqZXJjaWNpbyAyDQoNCkVuIHVuYSBmaW5jYSBkZSBVcGFsYSBzZSBlc3R1ZGlhbiAyMDAgw6FyYm9sZXMgZGUgY2FjYW8gY3Jpb2xsby4gRWwgY29sb3IgZGUgbGEgbWF6b3JjYSBlc3TDoSBjb250cm9sYWRvIHBvciB1biBnZW4gY29uIGRvcyBhbGVsb3M6IFIsIHF1ZSBwcm9kdWNlIG1hem9yY2Egcm9qYSB5IGVzIGRvbWluYW50ZSwgeSByLCBxdWUgcHJvZHVjZSBtYXpvcmNhIHZlcmRlIHkgZXMgcmVjZXNpdm8uIEVsIDYwICUgZGUgbG9zIMOhcmJvbGVzIHRpZW5lIG1hem9yY2Egcm9qYSwgY29uIGdlbm90aXBvcyBSUiBvIFJyLCB5IGVsIDQwICUgdGllbmUgbWF6b3JjYSB2ZXJkZSwgY29uIGdlbm90aXBvIHJyLg0KDQpFbCBjYWNhbyBzZSByZXByb2R1Y2UgbWVkaWFudGUgZmxvcmVzIHF1ZSBhcGFyZWNlbiBkaXJlY3RhbWVudGUgc29icmUgZWwgdHJvbmNvIHkgbGFzIHJhbWFzIGdydWVzYXMuIEVzdGFzIGZsb3JlcyBzb24gcG9saW5pemFkYXMgcHJpbmNpcGFsbWVudGUgcG9yIHBlcXVlw7FvcyBpbnNlY3RvcywgY29tbyBsb3MgbW9zcXVpdG9zIGRlbCBnw6luZXJvwqAqRm9yY2lwb215aWEqLCBxdWUgc2UgZGVzcGxhemFuIGVudHJlIGZsb3JlcyBidXNjYW5kbyBhbGltZW50by4gUGFyYSBxdWUgdW5hIGZsb3Igc2UgY29udmllcnRhIGVuIG1hem9yY2EsIGVsIHBvbGVuIGRlYmUgbGxlZ2FyIGEgc3UgZXN0aWdtYSB5IGxhIHBsYW50YSBkZWJlIHNlciBjb21wYXRpYmxlIGNvbiBlc2UgcG9sZW4sIHlhIHF1ZSBlbCBjYWNhbyB0aWVuZSB1biBzaXN0ZW1hIGRlIGF1dG9pbmNvbXBhdGliaWxpZGFkIHF1ZSBmYXZvcmVjZSBlbCBjcnVjZSBlbnRyZSDDoXJib2xlcyBkaXN0aW50b3MuIFVuYSBmbG9yIGJpZW4gcG9saW5pemFkYSBzZSB0cmFuc2Zvcm1hIGVuIHVuYSBtYXpvcmNhIHF1ZSBjb250aWVuZSBzZW1pbGxhcywgeSBsYXMgc2VtaWxsYXMgc29uIGVsIHByb2R1Y3RvIHF1ZSBlbCBhZ3JpY3VsdG9yIGNvc2VjaGEuIEVsIGNvbG9yIGRlIGxhIG1hem9yY2EgZGVwZW5kZSBkZWwgZ2Vub3RpcG8gZGUgbGEgcGxhbnRhIHF1ZSBsYSBwcm9kdWNlLCBhc8OtIHF1ZSBsYSBjYW50aWRhZCBkZSBtYXpvcmNhcyBzYW5hcyBxdWUgcHJvZHVjZSBjYWRhIMOhcmJvbCBzZSByZWxhY2lvbmEgZGlyZWN0YW1lbnRlIGNvbiBzdSDDqXhpdG8gcmVwcm9kdWN0aXZvLg0KDQpEdXJhbnRlIHRyZXMgYcOxb3Mgc2UgcmVnaXN0csOzIGVsIHByb21lZGlvIGFudWFsIGRlIG1hem9yY2FzIHNhbmFzIHBvciDDoXJib2wuIExvcyDDoXJib2xlcyBkZSBtYXpvcmNhIHJvamEgcHJvZHVqZXJvbiAyMiBtYXpvcmNhcyBwb3IgYcOxbywgeSBsb3MgZGUgbWF6b3JjYSB2ZXJkZSBwcm9kdWplcm9uIDE4Lg0KDQoqKlN1cHVlc3RvIGRlIHRpZW1wbzoqKsKgdW5hIGdlbmVyYWNpw7NuIGVxdWl2YWxlIGEgMTAgYcOxb3MuDQoNCmBgYHtyfQ0KbGlicmFyeShkZVNvbHZlKQ0KDQptb2RlbG8gPC0gZnVuY3Rpb24odGltZSwgc3RhdGUsIHBhcmFtZXRlcnMpIHsNCiAgd2l0aChhcy5saXN0KGMoc3RhdGUsIHBhcmFtZXRlcnMpKSwgew0KICAgICMgUGFyw6FtZXRyb3MgZGUgZml0bmVzcw0KICAgIHdfQUEgPC0gMC42MCAgIyBGaXRuZXNzIHJlbGF0aXZvIGRlIEFBDQogICAgd19hYSA8LSAwLjQwICMgRml0bmVzcyByZWxhdGl2byBkZSBhYQ0KICAgIA0KICAgICMgRnJlY3VlbmNpYXMgYWzDqWxpY2FzDQogICAgcCA8LSBBICAjIEZyZWN1ZW5jaWEgZGUgQQ0KICAgIHEgPC0gYSAgIyBGcmVjdWVuY2lhIGRlIGENCiAgICANCiAgICAjIEZpdG5lc3MgcHJvbWVkaW8gcG9ibGFjaW9uYWwNCiAgICB3X3Byb20gPC0gcF4yICogd19BQSArIHFeMiAqIHdfYWEgIA0KICAgIA0KICAgICMgVGFzYSBkZSBjYW1iaW8gZW4gbGEgZnJlY3VlbmNpYSBkZWwgYWxlbG8gQQ0KICAgIGRwIDwtIChwICogcSAqICh3X0FBIC0gd19hYSkpIC8gd19wcm9tDQogICAgZHEgPC0gLWRwICAjIExhIHN1bWEgZGUgZnJlY3VlbmNpYXMgZGViZSBzZXIgMQ0KICAgIA0KICAgIHJldHVybihsaXN0KGMoZHAsIGRxKSkpDQogIH0pDQp9DQpgYGANCmBgYHtyfQ0KIyBDb25kaWNpb25lcyBpbmljaWFsZXMNCmVzdGFkb19pbmljaWFsIDwtIGMoQSA9IDAuNjAsIGEgPSAwLjQwKQ0KIyBQYXLDoW1ldHJvcw0KcGFyYW1ldHJvcyA8LSBOVUxMDQojIFRpZW1wbw0KdGllbXBvIDwtIHNlcSgwLCA1MCwgYnkgPSAwLjUpDQoNCiMjIFJlc29sdWNpw7NuIGRlbCBtb2RlbG8NCnNvbCA8LSBvZGUoeSA9IGVzdGFkb19pbmljaWFsLCB0aW1lcyA9IHRpZW1wbywgZnVuYyA9IG1vZGVsbywgcGFybXMgPSBwYXJhbWV0cm9zKSAjIyBFc3RhIGbDs3JtdWxhIG1lIHJlc3VlbHZlIGxhIGVjdWFjacOzbiBkaWZlcmVuY2lhbA0KaGVhZChzb2wpDQoNCnRhaWwoc29sKQ0KDQpgYGANCmBgYHtyfQ0KcGxvdChzb2wsDQogICAgIHhsYWIgPSAiVGllbXBvIiwgeWxhYiA9ICJGcmVjdWVuY2lhIGRlIGFsZWxvcyIsDQogICAgIGNvbCA9IGMoImJsdWUiLCAicmVkIiksIGx3ZCA9IDIpDQpgYGANCmBgYHtyfQ0KZGYgPC0gYXMuZGF0YS5mcmFtZShzb2wpDQpkZiA8LSBjYmluZCh0aWVtcG8gPSB0aWVtcG8sIGRmKQ0KDQpnZ3Bsb3QoZGYsIGFlcyh4ID0gdGllbXBvKSkgKw0KICBnZW9tX2xpbmUoYWVzKHkgPSBsb2coQSksIGNvbG9yID0gIkFsZWxvIEEiKSkgKyANCiAgZ2VvbV9saW5lKGFlcyh5ID0gbG9nKGEpLCBjb2xvciA9ICJBbGVsbyBhIikpICsNCiAgbGFicyh4ID0gIlRpZW1wbyIsIHkgPSAiRnJlY3VlbmNpYSBkZSBhbGVsb3MiKSArDQogIHNjYWxlX2NvbG9yX21hbnVhbChuYW1lID0gIkFsZWxvcyIsIHZhbHVlcyA9IGMoImJsdWUiLCAicmVkIikpICsNCiAgdGhlbWVfbWluaW1hbCgpDQoNCmBgYA0KDQoNCiMgRWplcmNpY2lvIDMNCg0KRW4gdW4gZW1iYWxzZSBkZSBsYSB2ZXJ0aWVudGUgZGVsIFBhY8OtZmljbyB2aXZlIHVuYSBwb2JsYWNpw7NuIGRlIDUwMCBpbmRpdmlkdW9zIGRlIGd1YXBvdGUuIEVzdGEgZXNwZWNpZSBkZSBjw61jbGlkbyBwcmVzZW50YSBkb3MgbW9yZm9zIGRlIGNvbG9yYWNpw7NuIHF1ZSBkZXBlbmRlbiBkZSB1biBnZW4gY29uIGRvcyBhbGVsb3MuIEVsIG1vcmZvIGRvbWluYW50ZSB0aWVuZSBjb2xvcmFjacOzbiBvc2N1cmEgY29uIGZyYW5qYXMgbGF0ZXJhbGVzLCBjb24gZ2Vub3RpcG9zIEFBIG8gQWEuIEVsIG1vcmZvIHJlY2VzaXZvIHRpZW5lIGNvbG9yYWNpw7NuIGNsYXJhIHkgdW5pZm9ybWUsIGNvbiBnZW5vdGlwbyBhYS4gRW4gbGEgcG9ibGFjacOzbiwgZWwgODQgJSBwZXJ0ZW5lY2UgYWwgbW9yZm8gZG9taW5hbnRlIHkgZWwgMTYgJSBhbCBtb3JmbyByZWNlc2l2by4NCg0KRWwgZ3VhcG90ZSB0aWVuZSBjdWlkYWRvIHBhcmVudGFsIGJpcGFyZW50YWwuIERlc3B1w6lzIGRlIGxhIHJlcHJvZHVjY2nDs24sIGxhIHBhcmVqYSBjb25zdHJ1eWUgdW4gbmlkbywgZGVmaWVuZGUgbG9zIGh1ZXZvcyB5IGx1ZWdvIHByb3RlZ2UgYSBsYXMgY3LDrWFzIGR1cmFudGUgc3VzIHByaW1lcmFzIHNlbWFuYXMuIEVzdGEgZGVmZW5zYSByZWR1Y2UgbGEgZGVwcmVkYWNpw7NuLCBwZXJvIGxhcyBjcsOtYXMgc2lndWVuIGVuZnJlbnRhbmRvIHJpZXNnb3MgY29tbyBsYSBjb21wZXRlbmNpYSBwb3IgYWxpbWVudG8sIGxvcyBkZXByZWRhZG9yZXMgZGUgbWF5b3IgdGFtYcOxbyB5IGxvcyBjYW1iaW9zIGVuIGxhIGNhbGlkYWQgZGVsIGFndWEuIExhIHN1cGVydml2ZW5jaWEgZGUgbGFzIGNyw61hcyBoYXN0YSBsb3MgNiBtZXNlcyBkZSBlZGFkIGRldGVybWluYSBjdcOhbnRvcyBpbmRpdmlkdW9zIGxsZWdhbiBhIHJlcHJvZHVjaXJzZS4gRW4gZXN0ZSBjYXNvLCBzb2JyZXZpdmnDsyBlbCA5MiAlIGRlIGxhcyBjcsOtYXMgZGVsIG1vcmZvIGRvbWluYW50ZSB5IGVsIDc4ICUgZGUgbGFzIGRlbCBtb3JmbyByZWNlc2l2by4NCg0KKipTdXB1ZXN0byBkZSB0aWVtcG86KirCoHVuYSBnZW5lcmFjacOzbiBlcXVpdmFsZSBhIDMgYcOxb3MuDQoNCmBgYHtyfQ0KbW9kZWxvIDwtIGZ1bmN0aW9uKHRpbWUsIHN0YXRlLCBwYXJhbWV0ZXJzKSB7DQogIHdpdGgoYXMubGlzdChjKHN0YXRlLCBwYXJhbWV0ZXJzKSksIHsNCiAgICAjIFBhcsOhbWV0cm9zIGRlIGZpdG5lc3MNCiAgICB3X0FBIDwtIDAuODQgICMgRml0bmVzcyByZWxhdGl2byBkZSBBQQ0KICAgIHdfYWEgPC0gMC4xNiAjIEZpdG5lc3MgcmVsYXRpdm8gZGUgYWENCiAgICANCiAgICAjIEZyZWN1ZW5jaWFzIGFsw6lsaWNhcw0KICAgIHAgPC0gQSAgIyBGcmVjdWVuY2lhIGRlIEENCiAgICBxIDwtIGEgICMgRnJlY3VlbmNpYSBkZSBhDQogICAgDQogICAgIyBGaXRuZXNzIHByb21lZGlvIHBvYmxhY2lvbmFsDQogICAgd19wcm9tIDwtIHBeMiAqIHdfQUEgKyBxXjIgKiB3X2FhICANCiAgICANCiAgICAjIFRhc2EgZGUgY2FtYmlvIGVuIGxhIGZyZWN1ZW5jaWEgZGVsIGFsZWxvIEENCiAgICBkcCA8LSAocCAqIHEgKiAod19BQSAtIHdfYWEpKSAvIHdfcHJvbQ0KICAgIGRxIDwtIC1kcCAgIyBMYSBzdW1hIGRlIGZyZWN1ZW5jaWFzIGRlYmUgc2VyIDENCiAgICANCiAgICByZXR1cm4obGlzdChjKGRwLCBkcSkpKQ0KICB9KQ0KfQ0KYGBgDQpgYGB7cn0NCiMgQ29uZGljaW9uZXMgaW5pY2lhbGVzDQplc3RhZG9faW5pY2lhbCA8LSBjKEEgPSAwLjg0LCBhID0gMC4xNikNCiMgUGFyw6FtZXRyb3MNCnBhcmFtZXRyb3MgPC0gTlVMTA0KIyBUaWVtcG8NCnRpZW1wbyA8LSBzZXEoMCwgNTAsIGJ5ID0gMC41KQ0KDQojIyBSZXNvbHVjacOzbiBkZWwgbW9kZWxvDQpzb2wgPC0gb2RlKHkgPSBlc3RhZG9faW5pY2lhbCwgdGltZXMgPSB0aWVtcG8sIGZ1bmMgPSBtb2RlbG8sIHBhcm1zID0gcGFyYW1ldHJvcykgIyMgRXN0YSBmw7NybXVsYSBtZSByZXN1ZWx2ZSBsYSBlY3VhY2nDs24gZGlmZXJlbmNpYWwNCmhlYWQoc29sKQ0KDQp0YWlsKHNvbCkNCmBgYA0KYGBge3J9DQpwbG90KHNvbCwNCiAgICAgeGxhYiA9ICJUaWVtcG8iLCB5bGFiID0gIkZyZWN1ZW5jaWEgZGUgYWxlbG9zIiwNCiAgICAgY29sID0gYygiYmx1ZSIsICJyZWQiKSwgbHdkID0gMikNCmBgYA0KYGBge3J9DQpkZiA8LSBhcy5kYXRhLmZyYW1lKHNvbCkNCmRmIDwtIGNiaW5kKHRpZW1wbyA9IHRpZW1wbywgZGYpDQoNCmdncGxvdChkZiwgYWVzKHggPSB0aWVtcG8pKSArDQogIGdlb21fbGluZShhZXMoeSA9IGxvZyhBKSwgY29sb3IgPSAiQWxlbG8gQSIpKSArIA0KICBnZW9tX2xpbmUoYWVzKHkgPSBsb2coYSksIGNvbG9yID0gIkFsZWxvIGEiKSkgKw0KICBsYWJzKHggPSAiVGllbXBvIiwgeSA9ICJGcmVjdWVuY2lhIGRlIGFsZWxvcyIpICsNCiAgc2NhbGVfY29sb3JfbWFudWFsKG5hbWUgPSAiQWxlbG9zIiwgdmFsdWVzID0gYygiYmx1ZSIsICJyZWQiKSkgKw0KICB0aGVtZV9taW5pbWFsKCkNCmBgYA0KYGBge3J9DQojIGVzcGVyYW56YSBkZSB2aWRhIGd1YXBvdGUgQUEgZG9uZGUgbGEgcG9ibGFjacOzbiBkZSBHdWFwb3RlIGFsZWxvcyBBQSBlcyBkZSA0MjAgeSAzODYgc29icmV2aXZlbg0KRml0bmVzc0FBIDwtIDM4Ni80MjANCkZpdG5lc3NBQQ0KYGBgDQpgYGB7cn0NCiMgZXNwZXJhbnphIGRlIHZpZGEgZ3VhcG90ZSBhYSBkb25kZSBsYSBwb2JsYWNpw7NuIGRlIEd1YXBvdGUgYWxlbG9zIGFhIGVzIGRlIDgwIHkgNjIgc29icmV2aXZlbg0KRml0bmVzc2FhIDwtIDYyLzgwDQpGaXRuZXNzYWENCmBgYA0KDQo=