# 2015年美赛O奖论文A题Problem

Team Control Number

For T1 T2 T3 T4

office use only ________________ ________________ ________________ ________________

35532
Problem Chosen

A
Summary

For F1 F2 F3 F4

office use only ________________ ________________ ________________ ________________

2015 Mathematical Contest in Modeling (MCM) Summary Sheet

Team # 35532

1 / 34

Introduction
Problem background
A complex epidemic of Zaire Ebola virus (EBOV) has been affecting West Africa since approximately December 2013, with first cases likely occurring in southern Guinea [1] and facilitating several transmission chains to progress essentially unchecked in the region and to cross porous borders with neighboring Sierra Leone and Liberia and seed a limited outbreak in Nigeria via commercial airplane on 20 July 2014 [2]. Then the number of new cases appear an exponential growth. While public health interventions have been introduced in all affected countries, the numbers of infected cases and deaths from EBOV continue to increase due to the loss of effective medication. A total of 22,495 cases, with 8,981 deaths, have been reported to the World Health Organization as of 4 February 2015[3].

Our work
Since new anti-Ebola drug has been developed, a realistic, sensible, and useful mathematical model about drug allocation and delivery is necessary. Depending on the goal, we divide the spread of disease into three periods: naturally spreading period, spreading period with intervention and spreading period with effective medications. And then we develop a SEIR (susceptible-exposed-infectious-recovered) model to simulate the spread of the disease in the primary period, based on the changing trend of the numbers of infected cases and deaths. Finally, the model are improved to a SEIIR (susceptible-exposed-infectious-isolated-recovered) model to adapt to the second and third period and to predict the future trends in Guinea, Sierra Leone and Liberia. A delivery system include two main parts: numbers and locations of medical central delivering drugs decision and drug allocation of every medical center. For this reason, an evaluation system of every country’s long-term epidemic situation, two multi-objective optimization are built separately. We developed a detailed plan via our models.

Symbols and definitions
(1) Symbol for SEIIR model with effective medication
Symbol definitions susceptible individuals exposed individuals infectious individuals isolated infectious individuals non- isolated infectious individuals recover and survive individuals death case

S
E

I

I isolated

I free
R D

Team # 35532

2 / 34

N

total population of the area The transmission rate of EBOV the average durations of incubation. the average durations of infectiousness the case fatality rate the ratio of new non-isolated infective patients the non- isolated individuals fatality rate the isolated individuals fatality rate number of drug supply

?
1?

1?

f
?
f1 f2

G (t )
(2)Symbol for model II Symbol

definitions

p

q

the number of capitals of first-level administrative divisions in a country the number of medical centers for EVD set in a country;

V ? ?v1 , v2 ,
X ? ? x1 , x2 ,

, vi ,
xj

, vp?
, xq ?

the set of capitals the set of medical centers he shortest distance from x to y the metric of quantity of drug needed of the i th administrative in the coming period the population of the nth administrative division the area of the nth administrative division the population of the main city of the nth administrative division. the area of the main city of the nth administrative division. the temperature of the nth administrative division the Infected people number of the nth administrative division is represented

d ( x, y )
hi Rn Mn

rn : mn :
Bn

In

Team # 35532

3 / 34

Fn
Cn

open degree of Ebola treatment centers the temperature of the nth administrative division open degree of Ebola treatment centers the number of open ETCs in the country the number of cities in the country The economy of the n th country The education level of the n th country i age coefficient

Fn

copen
call
En

Jn
K
(3)Symbol for model III Symbol

definitions

G
Gm

The average daily production of drugs. the maximum production rate of the drugs. growth trends in drug production speed the aggregate of drug the first day of the time period of drug delivery. the time period of drug delivery
the number of non-isolated infectious individuals

rt
P tn
tn

Tp

Ii , free
Gi (t )
Ni Di

i th country’s number of the drug i th country’s the total population

i th country’s death cases

Assumptions
1. In early transmission period, EBOV spread in absence of control interventions. 2. The average duration of the incubation and infectious period were fixed to previous estimates from an outbreak of the same EBOV subtype in Congo in 1995 ( 1 ? = 5 .3 days and 1 ? =
5 .61 days) [3]. 3. Infection only occurs between the patients who are not isolated and susceptible.

Team # 35532

4 / 34

4. Even if without effective medication, the mortality rate of isolated patients is lower than that of free patients. 5. People who have recovered from Ebola will obtain the long-term immunity against Ebola. 6. Only being sent to medical center and isolated, affected people have the chance to receive therapy. 7. Ignoring patient’s different condition, the quantity of the medicine needed by all affected people whose disease is not advanced during the whole therapy is equal to one portion of drugs. 8. The population of the region we study is a constant N 9. The transmission only occurs btween only the susceptible and the infective. 10. Patients isolated will not affect the suspective. In other word,their effective contact rate equals to zero. 11 .The medication only refers to drugs used therapy of Ebola, not include vaccine. 12. If affected people want to accept therapy, they must get to a medical center set as delivering medication. 13. If a patient in whom symptoms of Ebola has appeared don’t get a medical central (usually the nearest medical center) in 24 hours, the new medication won’t affect his disease. 14. We set all the medical centers at the capitals of first-level administrative divisions. 15. We use the capital of an administrative division as a representation of the administrative

Model I
SEIR model
After transmission of the virus, susceptible individuals S enter the exposed class E before they become infectious individuals I that either recover and survive R or die D . Fig.1 shows the process.

Fig.1 the process of the model

According to the it，people can be divided into four class: ? the susceptible ( S ), who are susceptible to infection. ? the exposed( E ), who are affected but in the incubation. ? the infectious ( I ), who are infected and have the symptom. ? the recovered( R ), who recover or survive. The total population of the area is N , and N ? S ? E ? I ? R, according to the assumption 8. Numbers of all kinds of people in the initial time are shown as follows:

S (0) ? S0 ? 0 , E (0) ? E0 ? 0 , I (0) ? I 0 ? 0 , R(0) ? R0 =0
Let ? be the contact rate in absence of control interventions. S N is the proportion of the susceptible in total population. Hence, the effective contact rate is ? S N , and in a unit of time,

Team # 35532

5 / 34

the number of new patients is ? SI N .In other words, the rate of change in the number of the susceptible is

dS ? SI . ?? dt N

(1)

We assumed that ? E is the number of people newly entering the infectious. Therefore,

1 ? is the average durations of incubation. As for the number of people newly entering the
exposed every day, it can be calculated as follows:

dE ? SI ? ? ? E (2) dt N
In the same manner, we calculate the increment of the infective people:

dI ? ? E ?? I , dt
Where ? is sum of natural cure rate and mortality of the infectious. 1 ? can be seen the average durations of infectiousness. Suppose the case fatality rate is given by f , then

(1 ? f )? I is the number of death per day. Thus, the changing rate of recovery case and death
case with time can be represented as follows respectively:

dR ? (1 ? f )? I , (3) dt dD ? f?I . dt
From the above analysis results, We described the early transmission of EBOV as a SEIR (susceptible-exposed-infectious-recovered) dynamics model.

Result
The average duration of the incubation and infectious period were fixed to previous estimates from an outbreak of the same EBOV subtype in Congo in 1995 ( 1 ? = 5 .3 days and 1 ? = 5 .61
days).The other parameter has been given by Chowell et al are listed in Table 1.
[4]

. All important parameters’ value

Table 1. Parameter estimates for the 2014 EBOV outbreak

Parameter The average durations of incubation, 1 ?

Guinea

Sierra Leone

Liberia

5.3

5.3

5.3

Team # 35532

6 / 34

the average durations of infectiousness, 1 ? Transmission rate,

5.61 0.27 0.75
2 December 2013

5.61 0.43 0.49
19 Apr 2014

5.61 0.29 0.69
16 Apr 2014

?

Case fatality rate, f Date of appearance of first infectious case, T

Based on facts, we set the initial numbers of all kinds of as follows:

S (0) ? N ? E0 ? I 0 , E (0) ? 10 , I (0) ? 1 , R(0)=0
Due to the complexity of the system of differential equations, we can’t gain its analytical solution. Therefore, we solve it numerically ode toolbox in MATLAB 2014B for three country (Guinea, Sierra Leone and Liberia) and compare the results with data in the early EBOV transmission period (shown in Fig.2), since EBOV spread in absence of control interventions in that stage.

Team # 35532

7 / 34

Fig.2.1the early EBOV transmission period of Guinea

Fig.2.2the early EBOV transmission period of Liberia

Fig.2.3the early EBOV transmission period of Sierra Leone

Fig 2. Dynamics of 2014 EBOV outbreaks in Guinea, Sierra Leone and Liberia. Outbreak data for Guinea, Sierra Leone and Liberia were based on the cumulative numbers of reported total cases (confirmed, probable and suspected) and deaths from the World Health Organization (WHO) [5].

As we can see, The model fits the reported data of cases and deaths in Guinea, Sierra Leone and Liberia well.

SEIIR model 也许大学四年，我们会一直在迷茫中度过，因为生活总是难以言说。赛氪APP与您相伴！

Team # 35532

8 / 34

Generally, when the affected people is much enough, it will draw the attention of the government, and the epidemic will be declared as a public health emergency. So does EVBO. According to WHO’s suggestion [6], the first measurement that government should take is isolate the patients, which will impede the spreading of epidemic apparently. The SEIR (susceptible-exposed-infectious-recovered) model fit the early spreading of the disease, but it don’t consider this point. Therefore, there needs to be a improvement for the it. To achieve this goal, we develop a new SEIIR(susceptible-exposed-infectious-quarantined -recovered) model based on the SEIR model, considering the isolated people and without efficient medication. In this model, as transmission of the virus, susceptible individuals S enter the exposed class E before they become insolated infectious individuals I isolated or uncontrolled infectious individuals I free that either recover and survive R or die D . Fig.3 shows the process.

Fig.3 the process of the SEIIR model

The total population:

N ? S ? E ? I free ? Iisolated ? R
Considering the assumption 10, formulas (1) changes as follows:

? SI f r e e dS ， ?? dt N
Similarly, formula (2) and (3)

dE ? SI f r e e ? ?? E ， dt N
dI free dt ? ?? E ? ? I free , (7).

And the numbers of isolated patients’ changing rate is:

The coefficient ? in Formula (7) and (8) represents the ratio of new non-isolated infective patients in new infective patients, since not all of patients are willing to be isolated out of different reason. The value of ? mainly depends on the development level of medical system. Suppose that:

dIisolated ? ? (1 ? ? ) E ? ? Iisolated , (8) dt

Team # 35532

9 / 34

(1 ? ? ) ? ? ;
If all the patients in isolation can be cured, all the infective patients will be willing to be isolated, and ? ? 0 ; The lower the mortality of isolated patients is, the smaller the value of ? is. According to it, we can describe the relationship between ? and f 2 with a function:

? ? e? f

2

Obviously, the mortality of isolated patients f1 is much lower that of non-isolated patient f 2 . So, the recovered people can be divided into two parts according to being isolated or not. Thus, the changing rate of recovery case and death case with time can be represented as follows respectively:

dR ? (1 ? f1 )? I free ? (1 ? f 2 )? Iisolated dt
dD ? f1? I free ? f 2? I isolated dt
For the initial value, S (0) ? 0 , E (0) ? E0 ? 0 , I isolated (0) ? 0 I free (0) ? 0 , R(0) ? 0 Simulation Result: According to the report of WHO, the fatality rate of Ebola range from 0.6 to 0.8. Let f1 ? 0.8 and f 2 ? 0.6 while other coefficient unchanged. Due to the naturally spreading stage before isolation measures being taken, we set the initial value with the data in the 100th day as follows:

S0 ? 12219576 , E0 ? 2209 , I 0, free ? 822 , I 0,quarantined ? 822 , R0 ? 971, D0 ? 2912
In the same manner, we get a numerical solution by using the ode toolbox in Matlab201.The results shows in Fig.4.

Team # 35532

10 / 34

Fig.4 the result of SEIIR model The SEIIR model indeed shows the blockage of isolation to the Ebola’s spreading. SEIIR model with effective medication Now since the World Medical Association has developed a new medication which could stop Ebola and cures patients whose disease is not advanced. Outside of isolation , we can eradicate Ebola with the mediation, which has more effect on Ebola’s spreading. Base on this former work, we improve our model once again and get the SEIIR model with effective medication. The model can be described by a system of differential equations as follows

? SI free dS , ?? dt N
dE ? SI f r e e ? ?? E , dt N
dI free dt ? ?? E ? ? I free ,

dIisolated ? ? (1 ? ? ) E ? ? Iisolated ? M (t ) , dt

dR ? (1 ? f1 )? I free ? (1 ? f 2 )? Iisolated ? M (t ) , dt dD ? f1? I free ? f 2? I isolated ? M (t ) , dt

S (0) ? 0 , E (0) ? E0 ? 0 , I isolated (0) ? 0 , I free (0) ? 0 , R(0) ? 0
Where M (t ) is the number of people who can receive the medication in a unit of time. According assumption 12, the M(t) have no effect on the changing rate of I free . And other parameters and coefficients’ definition are the same as before. In addition, as M (t ) increase, more and more patients will be being isolated, in other word,

Team # 35532

11 / 34

? becomes smaller. The results shows in Fig.5.

Fig.5 the result of SEIIR model with effective medication

The SEIIR model with effective medication indeed shows the blockage of isolation to the Ebola’s spreading.

Model II
Since the new medication still cannot be manufactured on a large scale, which indicates a shortage of the drug, a feasible delivery system is necessary based on epidemic situation of countries and districts. The setting of the delivery system can be divided into two main parts: ? the medical center number and location decision ? and the medication allocation for every affected countries and districts.

Number and location decision
Assumption:
Affected people must get to a medical center within a specified time because the new medication can only cure patients whose disease is not advanced. According to Wikipedia, we assume that if a patient in whom symptoms of Ebola has appeared don’t get a medical central (usually the nearest medical center) in 24 hours, the new medication won’t affect his disease. Accounting for the poor transport facilities in West Africa and convenience of drug managing, we set all the medical centers at the capitals of first-level administrative divisions. We simplify the distance from everywhere in an administrative division to a medical center as the distance from the capital of the administrative division to the medical center. In other words, we use the capital of an administrative division as a representation of the administrative. It’s feasible because most of people in the West Africa tend to gather in several big cities.

Model:
Generally, the aid work in the charged of some international organizations, like WHO, They, rather than every country separately, assign the medical workers and allocate the medication and other medical source. Therefore, the new medication is delivered to countries in need separately by air, then to medical centers by highway and be used for therapy of patients there.

Team # 35532

12 / 34

Our first goal is decide the number and locations of medical centers. Due to the assumption, there must be at least one medical center within a specified distance from every administrative division to ensure patients can get to a medical center. This suggests the number and location decision is a set covering problem. In set covering problems, the objective is to minimize the number or cost of facility location such that a specified level of coverage is obtained. To describe the issue more clearly, we define some parameters as follows: p , the number of capitals of first-level administrative divisions in a country; q , the number of medical centers for EVD set in a country;

V ? ?v1 , v2 ,
X ? ? x1 , x2 ,

, vi ,
xj

, v p ? , the set of capitals;
, xq ? , the set of medical centers;

d ( x, y ) , the shortest distance from x to y ;
Due to the assumption,

X ?V .
Therefore,

x j ?V , j ? 1, 2,

,q .

Then we let d (vi , X q ) be the distance from the i th administrative division to the nearest medical central:

d (vi , X q ) ? min ?d (vi , x j )? .
1? j ? q

If there is a medical central in it,

d (vi , X q )=0
Suppose that the distance patients can moved within 24h is ? , because the assumption,

max d ? vi , X q ? ? ? .
1?i ? p

?

?

Considering that healthcare works, financial resources and devices are limited, we should minimize the number of medical centers for EVD, q ,on the premise of meeting the coverage requirement. In addition, administrative division’s need in medication vary from each other. Obviously, the larger an administrative division’s quantity of the medicine needed, the closer the nearest medical center should be. So, we also should

Min

? h d (v , X
i i 1

p

q

),

Where hi is the metric of quantity of drug needed of the i th administrative in the coming period. In next part, we will clarify how to calculate hi .

Team # 35532

13 / 34

To sum up, We model the problem about medical centers number and location decision with multi-objective optimization. The formulas of this model is

Min ? hi d (vi , X q )
i ?1

p

Min q X q ? ? x1 , x2 ,
Solution:
1?i ? p

s.t. max d ? vi , X q ? ? ? , xq ? ,q x j ?V , for i ? 1, 2,

?

?

Without loss of generality, we take the most seriously effected Sierra Leone for example and locate medical centers for it. The Republic of Sierra Leone is composed of four regions: the Northern Province, Southern Province, the Eastern Province, and the Western Area. The first three provinces are further divided into 12 districts. Fig.6 shows it.

Fig.6 The 12 districts and 2 areas of Sierra Leone. Source:

According to our model, we replace every district with its capital (the western area be seen as a district and Freetown is its capital). So , the value of p equals 13. We search and calculate the shortest path and distance between any two capitals via google map, and the results be shown in Table 2. Table 2 :the shortest path and distance between any two capitals Kail Ken Koi Ma Ka Po Ka Magb Bo Bon Moy Puje Free ahu ema du ken mbi rt bal uraka the amb hun town n i a Lo a a ko Kailah 0 112 46 204 353 30 32 179 169 261 284 222 417

Team # 35532

14 / 34

un Kene ma Koidu Make ni Kamb ia Port Loko Kabal a Magb uraka Bo Mattr u jong Moya mba Pujeh un Freet own 112 46.4 204 353 300 321 179 169 261 284 222 417 0 113 192 330 227 309 167 68.5 149 173 110 305

.4 11 3 0 15 9 30 8 25 3 27 6 13 3 16 2 25 6 25 3 21 3 34 0 192 159 0 152 97. 7 118 26. 8 133 217 124 206 185 330 308 152 0 51. 1 269 176 263 298 181 337 126

0 22 7 25 3 97 .7 51 .1 0 21 6 12 2 21 0 24 4 12 8 28 4 11 9

1 30 9 27 6 11 8 26 9 21 6 0 10 5.9 14 4 33 4 24 2 32 3 30 2 167 133 26.8 176 122 105.9 0 108 191 100 216 208 68. 5 162 133 263 210 144 108 0 82. 4 106 75. 3 174 .66 149 256 217 298 244 334 191 82. 4 0 130 98. 5 237 173 253 124 181 128 242 100 106 130 0 179 114 110 213 206 337 284 323 216 75.3 98.5 179 0 312 305 340 185 126 119 302 208 174. 66 237 114 312 0

As for the value of ? , accounting for the time of preliminary judgment and poor transport facilities, a reasonable value is 92.1km. Measuring the need of every district hi

Kailahun Kenema Kono Bombali Moyamba

h 0.512665 0.553363 0.22831 0.441547 0.505654

Table.3 the revolution of city h Koinadugu 0.379487 Kambia Tonkolili 0.294881 Port Loko Bo 0.504783 Pujehun Bonthe 0.271492 Western Area

h 0.341398 0.53602 0.334563 0.436276

This optimal problem belongs to the NP-hard problems that can’t be solved in polynomial time. So, we can only get the approximate solution with a heuristic, such as genetic algorithm, simulated annealing and so on. Here, we select the genetic algorithm and write a program with Matlab2014.

Team # 35532

15 / 34

Running the program repeatedly, we obtain a relatively good results and make our plan: There are seven mental centers located in Kailahun, Makeni, Kambia, Kabala, Bo, Moyamba and Freetown respectively. The Fig. 7 shows the scope that every mental center manages.

Fig.7 locations of medical centers in Sierra Leone

locations

of

medical

centers

in

Part III
1. Specify evaluation norms As for the evaluation the Susceptible standard for Each region, there are mainly aspects that count: Population density, Temperature, Humidity, the situation of current disease infection, Open or under construction， economy， National Education and age. What follows in the section will hammer at accounting for the eight aspects. ? Population density: Population density undoubtedly accounts for key proportion in the demand level of the drug evaluation. Because Ebola is a spirited contagion disease, not only patients but also corpses can be the source of infection. Therefore, if the population density is high, the rate of spread disease will be fast in general situation. However, in some area, the pollution density of the area is low, but almost people live in a central city (in fact, the pollution density of the city is high), we deal with data preprocessing, making the aspect more reasonable evaluation model. The pollution density of the area could be calculated as follows:

a1n ?
Rn : the pollution of the number n area. M n : the acreage of the number n area.

Rn Mn

The pollution density of the main city in the area could be calculated as follows:

Team # 35532

16 / 34

a2 n ?

rn mn

rn : the pollution of the main city number n area . mn : the acreage of the main city number n area.
We define a new population density which could be calculated as follows:

An ? a1n ? b1 ? a2 n ? b2 ? b1 ? b2 ? b3 ? 1
The numerical value depend on the (a2 n ? a1n ) .

(a2 n ? a1n ) ? b3 a2 n

By the above calculation, the population density in the evaluation of the model is more reasonable and accurate. ? Temperature Form the past Information and data, the Ebola usually outbreak in the summer ,but in the winter, the entire epidemic have been more effective control, and the status of major outbreaks in Africa is the lower latitudes area, which shows that the temperature is an important aspects in the demand level of the drug evaluation. The temperature of the area is:

Bn
? Humidity Ebola outbreak in the areas which the annual precipitation is high and harmonious, especially in the tropical rain forest areas with high coverage. The World Health Organization's report shows that humid climate is conducive to the spread of the virus, therefore, we choose the humidity as an aspects in the evaluation model. The temperature of the area is

Cn

? current disease infection
Current infection situation has a great impact on the future of the epidemic. The country who has large number of infected persons face the challenges in future, and need support from World Health Organization. Every patient will require substantial and large investment Because Ebola is highly infectious diseases. Therefore, we believe that using the number of patients to measure a country's entire national epidemic of severe is reasonable. The Infected people number of the area is

In
? Open or under construction

Team # 35532

17 / 34

As we all know, controlling the flow of people is an effective measure to prevent further deterioration of infectious diseases, On the one hand, under construction can control the spread of infection to uninfected areas, on the other hand, it have the ability to reduce the risk of infection of healthy people in infected area. In a country, some cities under opening while some cities under closing (we only consider the large cities of the country). We define Open coefficient which could be calculated as follows:

Fn ?

copen call

copen ： the number of open cities in the country
call ： the number of all cities in the country.

? Economy
In general situation, the more backward economy countries is the more difficult to control the disease is. Therefore, the economic situation of a country is an important indicator to determine the national level in the medical facilities. The lower the country's GDP means the more difficult to control Ebola infection because of under-developed Medical infrastructure. We use the GDP per capita as the aspects to show the economy of the area. The economy of the area is：

En
? National Education [7] Cremation in Africa do not comply with the local traditions, which is an important reason for the large-scale outbreak of Ebola. Ebola can spread through the corpses which lead a large number of people infected at the funeral. According to the survey, the higher the education level is, the greater the proportion of accept cremation. Therefore, we choose each country received higher education as a measure of the proportion of the country's level of education. The national education of the area is：

Jn
? Age[8] According to the report of the World Health Organization, the probability of illness in young adults over the age of 14 is three times the age of 14, the probability of illness aged 65 or older is four times the age of 14. Thus, we can see that the higher proportion of the population over the age of 14 is, the greater the overall probability of illness is. We define age coefficient which could be calculated as follows:

K?

K14~65 ? K 65 K all

Model 也许大学四年，我们会一直在迷茫中度过，因为生活总是难以言说。赛氪APP与您相伴！

Team # 35532

18 / 34

Introduction:
TOPSIS (technique for order performance by similarity to ideal solution) is a useful technique in dealing with multiattribute or multi-criteria decision making (MADM/MCDM) problems in the real world[9]. It helps decision maker organize the problems to be solved, and carry out analysis, comparisons and rankings of the alternatives. Accordingly, the selection of a suitable alternative(s) will be made. However, many decision making problems within organizations will be a collaborative effort. Hence, this model will extend TOPSIS to evaluate the demand of the drug. A complete and efficient procedure for decision making will then be provided. ? Step one Since the different factors have greatly different in magnitude, we should normalized the value of this factors. the normalized value X k of the decision matrix X can be any linearij scale transformation to keep 0 ? X k ? 1 . ij We consider that the normalized value of X k is the value of the corresponding element X k ij ij divided by the operation of its column elements, i.e., vector normalization, then (take the factor as example).
k Xm 1 ?

An
2 A ? A2 ? 2 1 2 ? An

? Step two: Construct decision matrix X The structure of the matrix can be expressed as follows:

X1 x1 x2 X? xi xm
k ? X 11 ? k ? X 21 ? ? k ? X i1 ? ? k ? ? X m1

X2
k X 12 k X 22

X X 1kj k X2 j
k X ij

Xn X 1kn ? k ? X2 n? ? k ? X in ? ? ? k X mn ? ?

X ik2
k Xm 2

k X mj

Where xi denotes the alternative i , i =1,2， , m ; X i represents the attribute or criterion j ,

j =1,2， , n ;
In the model, xi means different country, and X i means different factors in the evaluation. Observe that we can also set the outcomes of qualitative attributes from each alternative as discrete values, e.g., 1 to 10 or linguistics values, so that the quantitative values will be placed in the above decision matrix. For example, current disease infection is a factor which we can’t evaluate it directly, we can use some Index to show it[10]. ? Step three Determine the ideal and negative ideal solutions V ? and V - . k k

Team # 35532

19 / 34

Vk? = ? X 1k ? ,

k? , Xn ? ? max X ijk | j ? 1, 2, i

Vk- = ? X 1k - ,
? Step four

k, Xn

?? ? ? ?? min X
i

k ij

| j ? 1, 2,

?? , m ??
,m

Calculate the separation measure from the ideal and the negative ideal solutions, S ? and S - , i i respectively, for the group.

Si+ =

? (v
j ?1 n

n

k ij

? 2 ? vk j ) ，for alternativei, i ? 1,

,m

Si- =

? (v
j ?1

k ij

? 2 ? vk j ) ，for alternativei, i ? 1,

,m

? Step five Calculate the relative closeness C to the ideal solution for the group. Calculate the relative closeness to the ideal solution and rank the alternatives in descending order. The relative closeness of the i th alternative Ai with respect to positive idea solution can be expressed as

C?
Where 0 ? C ? 1

Si? , i ? 1, Si? ? Si?

,m

The larger the index value, the better the performance of the alternative.

2. Collect data
To evaluate h of every districts, Many data is needed, we can find every district’s data needed for our specific evaluation norms by searching from the website，like Wikipedia and WHO. We finally conclude the relative statistics of those areas and list them in a form. Table4:the data of the different district
Population density Kailahun Kenema Kono Bombali Kambia Port Loko Koinadugu Tonkolili Bo Bonthe Moyamba Pujehun 458.1595 4697.336 2618.687 267.3536 281.6377 860.1828 299.7132 1740.713 266.1477 1248.725 392.2557 1176.208 Open coefficient 100 100 0 100 50 100 50 0 100 50 100 50 Temperature (? C) 24 23 21 26 24 22 22 23 23 24 23 25 Humidity (%) 94 94 87 81 94 89 69 83 94 91 93 89 Current disease situation 4 4 3 2 3 5 4 4 4 1 4 3

Team # 35532 Western Area 25996.65 50 26 81 5

20 / 34

Note :1.the unit of the population denisty is People /per km

Then we use the data to calculate the h with TOPSIS, and the results has been shown Table 3.

The effect on ? of external factors
By analysis, we find that ? not only is associated with the nature of the epidemic but also affected by some external factors. These factors include population density, temperature, humidity, current epidemic situation ,open degree of ETC, economy, national education and age. Hence, ? vary from country to country. Some of these factors are change with time, which result ? changing with time. If we want to develop different reasonable long-term delivery system for different countries, evaluating the external factor’s impact on ? is needed. For the same reason, we can estimate the effect of external facts on ? with TOPSIS. We define a influence coefficient C .usually, C is non-negative, and as it increase, ? also appear a growth. The first four factors: population density, temperature, humidity, the current epidemic situation are defined as before. Obviously, their addition will result in increment of C. the other factors are redefined or defined as follows: ? Open degree of ETC(Ebola treatment centers) Considering there are many ETCs in a country, we redefine the open degree of ETC: h i which could be calculated as follows:

Fn ?

call copen

Where copen is the number of open ETCs in the country,

call is the number of cities in the country.
? Economy
In general situation, the more undeveloped a country is, the more difficult to control the epidemic is. Therefore, the economic situation of a country, measured by GDP per capital is an important indicator to determine the medical treatment level. The economy of the n th country is represented by：

En
? National Education

Team # 35532

21 / 34

Cremation do not comply with the local traditions in Africa, which is an important reason for the large-scale outbreak of Ebola. Ebola can spread through the corpses which lead a large number of people infected at the funeral. According to the survey, the higher the education level is, the greater the proportion of accept cremation is. Therefore, we choose the proportion of people received higher education in each country as a measure of the country's education level, whose increment results increment of ? The education level of the n th country is represented by：

Jn
? Age According to the report of the World Health Organization, the probability of being infected Ebola of young adults over the age of 14 is three times child below the age of 14, and the probability of poeple aged 65 or older is four times. Thus, we can see that the higher proportion of the population over the age of 14 is, the greater the ? will be. We define age coefficient which could be calculated as follows:

K?

K14~65 ? K 65 K all

Then we can use TOPSIS to estimate the variation of ? . Considering the epidemic mainly spread in Sierra Leone, Guinea and Liberia, we select their data to solve our models and decide their drug allocation. By searching online carefully, we final get data needed from the website of World Bank. Table 5: the data of the three country
current disease situation Guinea Liberia Sierra Leone 2917 8622 Temperature (?C) 27 25 Humidity (%) Open coefficient Economy (dollar) education (% of all) age Population density

75 91

25 50

492 297

38.13 45.16

57.69 57.11

45.5 35.5

10518

27

85

85

366

26.44

58.44

79.4

Note: the unit of the population is People /per km

After calculating ,the results are shown in Table Table 6 the initial influence coefficient of three country

C
Guinea Liberia Sierra Leone

0.28615361 0.467444301 0.728861892

Team # 35532

22 / 34

C will be used later.

Part IV
Assumptions
7 (1).Ebola drug production in the world can reach 5 ?10 doses every year[14], we can calculate

the increment rate of production speed with this data. (2).The maximal amount of drugs produced per day can finally reach 2 ?10 doses every year.
5

(3). Each patient needs 1 portion to rehabilitation, and each portion contain 10 doses drug. (4) Assumed that the drug distributed in every period of drug delivery are be used evenly every day.

Model:
Drug production model based logisitic. Because ? This model meets these characteristic that the production in initial time is small, then sharp increase in the mid-production .and production in the finally time is large and stable. ? the derivative of drug production rate increases at first and then decreases to zero. Therefore, we assume that the relationship between the drug production and time meets the logistic growth model, and the following equation can be listed:

dG G ? rt (1 ? ) dg Gm G (t0 ) ? 100

G :the daily production of drugs, namely the production speed.

Gm : the maximal production speed of drugs.

rt : growth rate of drug production speed.
amount of the drug delivered every time Our drug delivery plan think that there is only a drug delivery in a period of drug delivery

Tp . Let T ? 14days , Which is reasonable because it can save cost and don’t cause delay to
the Ebola’s controlling. Because drugs’ transporting need a certain period of time, and the drug output everyday, so we need to calculate the aggregate of drug which is the based of delivery.

Ptn ? ?

tn ?Tp

tn

Gdt

tn ：the first day of the time period of drug delivery.
Tp : the period of drug delivery
Delivery model According to former analysis, drug production and drug demand change over time.

Team # 35532

23 / 34

Comparing the drug production model and drug demand models, we can divide the whole duration into two stage: stage severely short in drugs and stage relative sufficient in drugs. now, we choose the 150th days as a time node.150th days ago, we faced with a severe drug shortage, and our purposes is to control the increment rate of the infectious .After 150th days, the production of the drugs achieve the certain number and the purpose become to reduce the number of death every day. Delivery model at the stage severely short in drugs Nonlinear programming model In the first stage of drug manufacturing, Because of the great shortage in drug, many patients can’t get effective medical treatment, and the increment of the infectious and deaths are both high. In this moment, we should give priority to the decrement of the increment of the infectious to make the epidemic becoming under the control as soon as possible which will benefit later epidemic controlling. Hence, For the first stage of drug manufacturing, we build a nonlinear programming model: ? the optimization goal which is minimizing the number of the infection sources of all countries:

Min

?I
i ?1

k

i , free

, i ? 1, 2,

,k

Where k is the number of countries needing drug assistance, and I i , free is the average number of infectious individuals non-isolated of the i th country in a period of drug delivery. ? ? the decision variables is the drug allocation distributed for each country noted by Gi (t ) , it refers the doses of drugs needed by the i th country. the resistant condition is actual number of drugs.

? G (t ) ? P , i ? 1, 2,
i ?1 i tn

k

,k

Therefore, a Nonlinear programming optimization model of drug distribution is modeled as follows:

Min
k

?I
i ?1
tn

k

i , free

,

? G (t ) ? P , i ? 1, 2,
i ?1 i

,k

i ? 1, 2

,k .

delivery model at the stage relatively sufficient in drugs As time goes on, the speed of drug manufacturing grows rapidly. After a period of time, the

Team # 35532

24 / 34

speed of drugs become high enough to meet the need of every countries to control disease. In the moment, our goal become to minimize death cases every day with the number of the infectious stable. So, a multi-objective optimization model will be developed next based on the former model for the second stage. Our goal can be described as follows:

Min
Where And

?
i ?1

k

dDi , dt

dDi is the average death cases per day in a period of drug delivery. dt

Min
At least making it stable. And there are two constraint condition, one is the drug amount restriction:

?I
i ?1

k

i , free

,

? G (t ) ? P , i ? 1, 2,
i ?1 i tn

k

,k .

The other is to consider the fairness to some extent:

(
Where

I i , free 1 k I i , free ? ? ) Ni k i ?1 Ni

1 k I i , free ? 0.2 ? k i ?1 N i

（11）.

Ni is i th country’s the total population.
I i , free Ni
is i th country’s percent of infectious individuals.

1 k I i , free is average percent ? k i ?1 N i

of uncontrolled infectious individuals. This shows the ration of difference between the proportion of patients non-isolated in all the population and the mean value of all the countries to the mean value of all the countries is less 20%. This will make the drug distribution is more evenly. The multi-objective optimization model can be described by following formulas:

min min s.t.

?D ,
i ?1 k i

k

?I
i ?1 k i ?1

i , free

? G (t ) ? P
i

tn

Team # 35532

25 / 34

(

I i , free 1 k I i , free ? ? ) Ni k i ?1 Ni

1 k I i , free ? 0.2 ? k i ?1 N i
1 k Di , free ? 0.2 ? k i ?1 N i

(

Di , free 1 k Di , free ? ? ) Ni k i ?1 N i

? Gi (t ) ? 0

Results & analysis
In the Drug production model, the formula (1 ? G Gm ) represent the resistance from the finite ability to manufacture drug. Obviously, as t goes on, rt becomes larger ， and

(1 ? G Gm ) becomes smaller ，the increment of manufacturing speed of drugs depend both
on them. According to the assumption (1) and (2), and supposed that at the 100th day, the manufacturing speed of drugs reaches the maximum, we can know that the maximal manufacturing speed of drugs

Gm = 2 ?105 .
And we solved the model with ode toolbox in Matlab2014 numerically. The result is shown in Fig.8.

Fig.8 the model based on logisitic

Since the epidemic mainly spread in Guinea, Liberia, Sierra Leone, we consider the situation where all the drugs are provided to these three country, and make a half-of-year drug delivery plan and predict the spreading trend in the half of year with the SEIIR model with efficient medication. Suppose we start to distribute the drug from the 100th day since the epidemic spread. Table.7

Team # 35532

26 / 34

shows the value of parameters that moment.
Table 7 the value of parameters that moment

S0
Guinea Liberia Sierra Leone

E0
4775 9215 7799

I 0, free
639 1228 1042

I 0,isolated
1209 2333 1079

R0
243 453 475

D0
1910 3686 3199

12181398 4129792 6223952

In our plan, the period of drug delivery is 14 days, therefore, we will make the new delivery plan every 14 days. As we know from the former analysis, within the half of year, the value of

? vary from country to country and will change with time. Due to this, to recalculate the value
of ? after a certain period time is also needed. We recalculate ? with the TOSIS mentioned above every month. Solving the optimization model with optimization toolbox in Matlab2014. We get our delivery plan for six month shown as following table and figure.

Infected case Guinea

Death case

portion

Infected Case Guinea

Death case 2984 5769 5032

portion 456 679 589

Liberia Sierra Leone

8380 19334
15586

2865 5529
4798

108 235
345

Liberia Sierra Leone

8730 23986 16142

Team # 35532

27 / 34

Infected case Guinea 16661 21805 16142

Death case 3284 6169 5648

portion 685 891 981 Guinea

Infected case 6905 18954 13834

Death case 3367 6398 5781

portion 986 1278 1189

Liberia Sierra Leone

Liberia Sierra Leone

Infect ed case

Death case

portion

Infected case

Death case

portion

Guinea 4101 3475 1230 Guinea 1771 3586 1580 8747 6430 1806 2767 6529 2209 Liberia Liberia 8639 5863 1698 3910 5965 2497 Sierra Sierra Leone Leone The figures show the spreading trends of the Ebola in the three country in the next half of year. The tables show the total death cases and total infected cases at the moments delivering drugs .and the portions delivered to three country every time. According to the assumption (3), a portion of drugs equals to ten doses of drugs which can be used to cure one patient totally.

Conclusion
The isolation method and drug intervention can effectively control the EBOV spreading. And Drug manufacturing can also influence epidemic control. The difference between each country will cause the epidemic situation is different. So making epidemic control Strategy need consider a great deal.

Team # 35532

28 / 34

Considering many factors including isolated method, drug intervention, geographical conditions, the differences between different countries, the characteristics of the spread of the virus, we build those models. Model I can forecast the Infectious disease situation. Model II can determine the location of drug station. Model III can make the optimal choice of each country's number of drug. All of these help eradicating Ebola.

Model extending
In SEIIR with effective drug intervention model, we assume that medication only refers to drugs used therapy of Ebola, not include vaccine; Medications only work on patients but the vaccine can be used for individuals who live in the Susceptible area . Now we consider drugs and vaccines can simultaneously supply. The vaccines production is F dose per day. The drug production is G dose per day. Vaccine is able to let the individual to produce antibodies against the virus. So this formula (1) can be rewritten as

? SI free dS ?? ? F (t ). dt N
At the same time putting the immune into the survivor population, then this formula (3) can be rewritten as

dR ? (1 ? f1 )? I free ? (1 ? f 2 )? Iisolated ? M (t ) ? F (t ). . dt
SEIIR with effective drug intervention model was extended to SEIIR with effective drug and Vaccine intervention model. At the same time F Consistent with the Logistic block growth model，then

dF F ? r1t (1 ? ) ， dg Fm

F (t0 ) ? 100.
Fm is maximum of the vaccines production rate， r1 is vaccines production growth rate.

Sensitivity Analysis
SEIIR with effective drug intervention model
There are some coefficients related to properties of the virus in SEIIR with effective drug intervention model: 1 ? : The average durations of incubation, 1 ? : the average durations of infectiousness, ? : the transmission rate. We analysis the effect of the three factors on the death case by changing them. Table 8:the three factors on the death case
Initial value Changed value 6 Changed percent 13.21%

?

1?

5.3

Team # 35532

29 / 34

the death case

119556 5.61 119556 0.27 119556

124870 6.3 131690 0.5 9343221

4.44% 12.31% 10.15% 46% 681.4%

1?
the death case

?
the death case

Fig. 9: The

changes, and death cases

The table 8 shows that ? is the factor which has the greatest effect on the death case among the properties, and fig.9 shows the result graphicall In conclusion, ? is a key factors to control Ebola, we can take some measures, such as control the population density, isolation and so on, to reduce it. (2) Change the E0 When solving the models, we have set initial values of the initial susceptive S 0 ,exposed

E0 ,infectious I 0 ,etc. As for S 0 , I 0 , we can get relatively data from the data statistics of WHO.

Fig.10:

creases, and death cases increases

Team # 35532

30 / 34

But we only can inaccurate E0 by estimation. So it is necessary to evaluate the effectiveness by changing E0 and observing the change of result. We change E0 from 10 to 1000 by step of 10. As we can see from fig. 10, the death case increases related to E0 increases. In conclusion, the E0 can affect the result greatly. When we set model for predicting the future trend, we should determine the value of ? , E0 based on the reality situation , which will make the finally forecasting more accurate.

(2) delivery systems model
When we build the drug delivery systems model, we use the Logistic block growth model describe the change of drug production speed with time. It’s a reasonable assumption. However, in the real world ,the relationship between the output of drug and time may vary from it. Not only that, the coefficients r and Gm can also influent the results ( r is the increase speed of output growth rate; Gm is the maximum daily output).Similarly, we study their influence by changing the value of them separately, ? The sensitivity to r . The results of sensitivity analysis to r is shown on Table.9. Table 9
r 50 100(100%) 200(400%) The days use to control 250 162(35.3%) 102(59.2%) The death case 629128 124870(80.1%) 272934(56.6%)

From the Table 2, we find an obvious phenomenon that the bigger the value is ,the earlier the situation become under control.
?

The sensitivity to Gm

The Gm have been changed from 100 to 10000 by step of 100. We see the needed time to make the epidemic under the control is changed with the change of

Gm in the Fig.11

Team # 35532

31 / 34

Fig.11 Epidemic controlling is very sensitive to the maximal daily output. The more the maximum daily output is, the less the death case is less. When the maximal daily output is under 4000 dose per day, the drug production can’t meet the need of epidemic controlling. When the maximum daily output is smaller, it is more difficult to control the epidemic and the death cases become high. Therefore, a way to control the epidemic and reduce the death cases is maximize the output every day and improve the increase speed of output growth rate.

Weakness and strengths
Weakness:
In model 1, according to the sensitivity analysis, the simulation results depend on largely the initial values of parameters and coefficients. But some values of the parameters can be only estimated, which can be inaccurate. In addition, the virus has a large probability of mutation, and some parameters will changed which we ignored, which will decrease the validity of the model. In model 2, to simplify the problem, we set all the medical centers at the capitals of districts. And when calculating the shortest distance, we represent a district with its capital . In model 3, we build the optimization model without considering the vaccine.

Strengths
In model 1, we considered the influence of isolated people and improved our basic model through adding isolated people and drug intervention, which make the simulation more real and effective. The reasonable setting of the parameter ( ? ) , will bring the model the ability of selfregulation. In model 2, when choosing locations of medical centers, considering the drug distribution point will not change once selected, we get our result from a long-term perspective. Which will make the results more reasonable. In model 3, we building a logistic block growth model about drug production. This model is

Team # 35532

32 / 34

reasonable. We consider two situations. The one is the condition of severe lacking drug, other one is the condition of relatively sufficient drug. At last we assume the disease is under control when the growth rate of infectious individuals develop into zero.

Letter for the world medical association
To the world medical association, Ebola has received wide attention of people all over the world since 2014, and until now, the situation of the Ebola in Africa is still not optimistic. Recently, the number of the infected people appeared an increment again. With the coming warm rainy season, we must pay attention to the development trend of Ebola disease. To prevent the spread of epidemic, we propose the following suggestion Based on our study: National measures to control the epidemic infection Build strong and effective isolation to prevent deterioration of the condition, and focus on the patient in the incubation to avoid them a new source of infection before being isolated. Provide the public with the right to preventive methods by publicity, and reduce large gatherings in the public. The government should popularize cremation, even take lethal force cremation for death for Ebola patients with known. The government should control the regional migration. For description of drug distribution Because the drug is in shortage, the need of each country medicines cannot be met. Therefore: Before drug production reach a certain amount, according to the different circumstances of each country and a reasonable prediction, we set a drug delivery system whose purpose is the overall controlled objectives, and the overall distribution of the situation is adjusted every two weeks. After drug production speed reach the necessary requirement .we use minimazing the number of deaths as the target for distribution. If you want to know more detailed analysis and the data, welcome to read our paper. We wish that Ebola will be controlled soon. Yours sincerely

References
[1] Baize S, Pannetier D, Oestereich L, et al. Emergence of Zaire Ebola virus disease in Guinea[J]. New England Journal of Medicine, 2014, 371(15): 1418-1425. World Health Organization: Ebola Virus Disease, West Africa –Update on 27.July 2014. 2014. [2] Gire S K, Goba A, Andersen K G, et al. Genomic surveillance elucidates Ebola virus origin and transmission during the 2014 outbreak[J]. Science, 2014, 345(6202): 1369-1372.

Team # 35532

33 / 34

[3] Ebola Situation Report - 4 February 2015

http://apps.who.int/ebola/en/ebola-situation-report/situation-reports/ebola-situation-report-4february-2015 Chowell G, Hengartner N W, Castillo-Chavez C, et al. The basic reproductive number of Ebola and the effects of public health measures: the cases of Congo and Uganda[J]. Journal of Theoretical Biology, 2004, 229(1): 119-126. [4] Althaus C L. Estimating the reproduction number of Zaire ebolavirus (EBOV) during the 2014 outbreak in West Africa[J]. arXiv preprint arXiv:1408.3505, 2014. [5] Lekone, P. E. & Finkenstadt, B. F. 2006 Statistical inference in a stochastic epidemic SEIR model with control intervention: Ebola as a case study. Biometrics 62, 1170– 7.doi:10.1111/j.1541-0420.2006.00609.x. [6] Hand hygiene for health workers caring for Ebola patients http://www.who.int/entity/csr/disease/ebola/hand-hygiene/en/index.html [7] Owen S H, Daskin M S. Strategic facility location: A review[J]. European Journal of Operational Research, 1998, 111(3): 423-447. [8] Ebola Situation Report - 4 February 2015 http://apps.who.int/ebola/en/ebola-situation-report/situation-reports/ebola-situation-report-4february-2015 [9] S.M. Belenson, K.C. Kapur, An algorithm for solving multicriterion linear programming problems with examples, Operational Research Quarterly 24 (1) (1973) 65–77. [10] C.L. Hwang, K. Yoon, Multiple Attribute Decision Making, Springer-Verlag, Berlin, 1981. [11] introduction about ebola http://www.who.int/csr/disease/ebola/en/ [12] National Education Data.worldbank.com [13] Ebola Situation Report - 4 February 2015 http://apps.who.int/ebola/en/ebola-situation-report/situation-reports/ebola-situationreport-4-february-2015 [14] WHO approves Ebola trial drug http://english.cntv.cn/2014/08/13/VIDE1407889322218675.shtml

2015 美赛题目及翻译
2015 美赛题目及翻译_理学_高等教育_教育专区。2015 Contest Problems MCM PROBLEMS PROBLEM A: Eradicating Ebola The world medical association has a ...
2017美赛MCM A题M奖获奖论文_图文
2017美赛MCM A题M奖获奖论文_理学_高等教育_教育...Its construction was controversial, and a 2015 ...Therefore, the research problems are: a) ...
2017美赛B题O奖论文56731中英对照版
2017美赛B题O奖论文56731中英对照版 - For office use only T1 T2 T3 T4 Team Control Number 56731 Problem Chos...
2014年美赛数学建模A题翻译版论文
2014年美赛数学建模A题翻译版论文_理学_高等教育_教育专区。2014年美国赛数学建模A题中文翻译 数学建模竞赛(MCM / ICM)汇总表 基于细胞的高速公路交通模型 自动机...
2015美赛题目
Click the title below to download a PDF of the 2015 ICM Problem C. 问题...2015美赛题目 1页 免费 2015美赛建模B题优秀论文... 24页 1下载券 2015...
2017年建模美赛A题带翻译
2017年建模美赛A题带翻译 - Problem A: Managing The Zambezi River The Kariba Dam on the Zambezi River is one...
2016数学建模美赛F题论文_图文
2016数学建模美赛F题论文_理学_高等教育_教育专区。...在2016年美国数学建模竞赛当中获得H奖的论文,希望对...The essence of this problem is kind of a case...
2015美赛H奖论文
2015美赛H奖论文_工学_高等教育_教育专区。2015美国...Problem Chosen D 2015 Mathematical Contest in ...o Economic growth reached a limit in the VM ?...

2014美赛数模A题H奖
2014美赛数模A题H奖_IT认证_资格考试/认证_教育专区...2014美赛A题 北师O奖 暂无评价 29页 免费 2014数学...©2015 Baidu 使用百度前必读 | 文库协议 | 网...