Đăng ký Đăng nhập
Trang chủ Bài toán vận tải có hạn chế khả năng lưu thông và ứng dụng...

Tài liệu Bài toán vận tải có hạn chế khả năng lưu thông và ứng dụng

.PDF
51
42
79

Mô tả:

„I HÅC THI NGUY–N TR×ÍNG „I HÅC KHOA HÅC Ph¤m B¡ T¼nh B€I TON VŠN TƒI C H„N CH˜ KHƒ N‹NG L×U THÆNG V€ ÙNG DÖNG Chuy¶n Ngh nh: TON ÙNG DÖNG M¢ sè: 60.46.01.12 LUŠN V‹N TH„C Sž TON HÅC Th¡i Nguy¶n - 2013 Soá hoùa bôûi trung taâm hoïc lieäu http://www.lrc-tnu.edu.vn/ 1 Möc löc Líi nâi ¦u Ch÷ìng 1. B€I TON VŠN TƒI C KHƒ N‹NG L×U THÆNG H„N CH˜ 5 . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 1.1. B i to¡n v  t½nh ch§t . . . . . . . . . . . . . . . . . . . . 5 1.2. Ph÷ìng ¡n cüc bi¶n ban ¦u . . . . . . . . . . . . . . . . 9 1.3. Ti¶u chu©n tèi ÷u . . . . . . . . . . . . . . . . . . . . . . 10 1.4. Thuªt to¡n gi£i . . . . . . . . . . . . . . . . . . . . . . . 10 1.5. V½ dö minh håa . . . . . . . . . . . . . . . . . . . . . . . 11 Ch÷ìng 2. B€I TON VŠN TƒI VÎI L×ÑNG CUNG BÀ CHN D×ÎI 18 2.1. Nëi dung b i to¡n . . . . . . . . . . . . . . . . . . . . . . 18 2.2. B i to¡n t÷ìng ÷ìng v  t½nh ch§t . . . . . . . . . . . . 21 2.3. Thuªt to¡n gi£i . . . . . . . . . . . . . . . . . . . . . . . 25 2.4. V½ dö minh håa . . . . . . . . . . . . . . . . . . . . . . . 28 Ch÷ìng 3. B€I TON VŠN TƒI VÎI L×ÑNG CUNG BÀ CHN 32 3.1. Nëi dung b i to¡n . . . . . . . . . . . . . . . . . . . . . . 32 3.2. T½nh ch§t nghi»m cõa b i to¡n . . . . . . . . . . . . . . 34 3.3. Thuªt to¡n gi£i . . . . . . . . . . . . . . . . . . . . . . . 36 3.4. V½ dö minh håa . . . . . . . . . . . . . . . . . . . . . . . 40 3.5. Tr÷íng hñp têng qu¡t . . . . . . . . . . . . . . . . . . . 42 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49 . . . . . . . . . . . . . . . . . . . . . . 50 K¸t luªn T i li»u tham kh£o Soá hoùa bôûi trung taâm hoïc lieäu http://www.lrc-tnu.edu.vn/ 2 LÍI NÂI †U Mæ h¼nh b i to¡n vªn t£i cê iºn vîi l÷ñng cung ð c¡c nìi giao h ng (gåi l  c¡c tr¤m ph¡t) v  l÷ñng c¦u ð c¡c nìi nhªn h ng (gåi l  c¡c tr¤m thu) ành tr÷îc ¢ r§t quen thuëc trong lþ thuy¸t tèi ÷u tuy¸n t½nh. B i to¡n vªn t£i d¤ng n y câ nhi·u ùng döng rëng r¢i trong thüc ti¹n v  ¢ ÷ñc nhi·u ng÷íi quan t¥m nghi¶n cùu v  ùng döng (xem [4], [5]). Mët mð rëng tü nhi¶n cõa b i to¡n vªn t£i l  khæng qui ành tr÷îc l÷ñng cung cõa c¡c tr¤m ph¡t v  l÷ñng c¦u cõa tr¤m thu m  cho ph²p l÷ñng cung cõa méi tr¤m ph¡t hay/v  l÷ñng c¦u cõa méi tr¤m thu câ thº thay êi trong mët kho£ng cho tr÷îc. Khi â ta g°p b i to¡n vªn t£i vîi r ng buëc hai ph½a v· l÷ñng cung (cõa c¡c tr¤m ph¡t) v  l÷ñng c¦u (cõa c¡c tr¤m thu). Mæ h¼nh mð rëng n y n£y sinh tø mët sè ùng döng, trong â câ v§n · i·u h nh m¤ng xe buþt ð c¡c th nh phè. B¬ng c¡ch dòng kÿ thuªt th¶m ©n sè phö, ta câ thº ÷a b i to¡n vªn t£i vîi r ng buëc hai ph½a èi vîi l÷ñng cung (c¦u) v· b i to¡n vªn t£i cê iºn vîi mët sè bi¸n câ r ng buëc cªn tr¶n, tùc l  qui v· b i to¡n vªn t£i câ h¤n ch¸ kh£ n«ng l÷u thæng v  do â câ thº ¡p döng ph÷ìng ph¡p th¸ và ¢ bi¸t º gi£i b i to¡n. Luªn v«n n y · cªp tîi b i to¡n vªn t£i câ h¤n ch¸ kh£ n«ng l÷u thæng, tr¼nh b y thuªt to¡n th¸ và gi£i b i to¡n v  ùng döng v o xû lþ b i to¡n vªn t£i vîi r ng buëc hai ph½a trong hai tr÷íng hñp: a) l÷ñng cung (hay c¦u) bà ch°n d÷îi bði mët sè khæng ¥m; b) l÷ñng cung (hay c¦u) bà ch°n (c£ tr¶n v  d÷îi). Tr÷íng hñp khi c£ l÷ñng cung v  l÷ñng c¦u ·u thay êi công s³ ÷ñc · cªp tîi. Luªn v«n gçm líi nâi ¦u, ba ch÷ìng, k¸t luªn v  danh möc t i li»u tham kh£o. Ch÷ìng 1 vîi ti¶u · "B i to¡n vªn t£i câ h¤n ch¸ kh£ n«ng l÷u Soá hoùa bôûi trung taâm hoïc lieäu http://www.lrc-tnu.edu.vn/ 3 thæng" tr¼nh b y nëi dung v  c¡c t½nh ch§t cõa b i to¡n vªn t£i, trong â l÷ñng h ng vªn chuyºn tø méi tr¤m ph¡t ¸n méi tr¤m thu khæng ÷ñc v÷ñt qu¡ mët mùc giîi h¤n qui ành tr÷îc (do n«ng lüc ph÷ìng ti»n vªn chuyºn câ h¤n ho°c do n«ng lüc c¦u ÷íng tr¶n tuy¸n vªn chuyºn â bà h¤n ch¸ ...). Ti¸p â, · cªp tîi ph÷ìng ph¡p t¼m ph÷ìng ¡n cüc bi¶n ban ¦u cõa b i to¡n. Sau â, tr¼nh b y cì sð lþ luªn v  nëi dung thuªt to¡n th¸ và (mët bi¸n thº cõa thuªt to¡n ìn h¼nh xû lþ b i to¡n câ bi¸n bà ch°n tr¶n) gi£i b i to¡n. Cuèi ch÷ìng n¶u ra v½ dö sè minh håa thuªt to¡n gi£i ¢ tr¼nh b y. C¡c ki¸n thùc ð ch÷ìng n y s³ c¦n ¸n ð c¡c ch÷ìng sau º xû lþ v  gi£i b i to¡n vªn t£i vîi r ng buëc hai ph½a. Ch÷ìng 2 vîi ti¶u · "B i to¡n vªn t£i vîi l÷ñng cung bà ch°n d÷îi" x²t b i to¡n vªn t£i, trong â gi£ thi¸t l÷ñng cung cõa c¡c tr¤m ph¡t khæng qui ành tr÷îc m  ch¿ bà ch°n d÷îi (lîn hìn mët mùc tèi thiºu n o â, th÷íng l  sè d÷ìng), cán l÷ñng c¦u cõa c¡c tr¤m thu ¢ bi¸t tr÷îc. N¶u mæ h¼nh to¡n håc cõa b i to¡n v  n¶u c¡ch ÷a v· b i to¡n vªn t£i câ h¤n ch¸ kh£ n«ng l÷u thæng. Sau â tr¼nh b y ti¶u chu©n tèi ÷u v  thuªt to¡n gi£i b i to¡n. Cuèi ch÷ìng n¶u v½ dö sè minh håa thuªt to¡n ¢ tr¼nh b y. Ch÷ìng 3 vîi ti¶u · "B i to¡n vªn t£i vîi l÷ñng cung bà ch°n" x²t b i to¡n vªn t£i, trong â gi£ thi¸t l÷ñng cung cõa c¡c tr¤m ph¡t khæng qui ành tr÷îc m  thay êi trong mët kho£ng cho tr÷îc, tùc l  l÷ñng cung bà ch°n (c£ tr¶n v  d÷îi), cán l÷ñng c¦u cõa c¡c tr¤m thu ¢ bi¸t tr÷îc. B i to¡n x²t ð ¥y mð rëng hìn æi chót so vîi b i to¡n x²t ð ch÷ìng 2 v  do â c¡ch xû lþ công phùc t¤p hìn. N¶u mæ h¼nh to¡n håc cõa b i to¡n v  n¶u c¡ch ÷a v· b i to¡n vªn t£i câ h¤n ch¸ kh£ n«ng l÷u thæng. Sau â tr¼nh b y ti¶u chu©n tèi ÷u v  thuªt to¡n gi£i b i to¡n. Cuèi ch÷ìng n¶u v½ dö sè minh håa thuªt to¡n gi£i. Do thíi gian v  ki¸n thùc cán h¤n ch¸ n¶n ch­c ch­n luªn v«n cán câ nhúng sai sât nh§t ành, k½nh mong quþ th¦y cæ v  c¡c b¤n âng gâp þ ki¸n º t¡c gi£ ti¸p töc ho n thi»n luªn v«n n y. Nh¥n dàp n y, t¡c gi£ xin b y tä láng bi¸t ìn s¥u s­c ¸n Th¦y h÷îng Soá hoùa bôûi trung taâm hoïc lieäu http://www.lrc-tnu.edu.vn/ 4 d¨n GS - TS Tr¦n Vô Thi»u ¢ tªn t¼nh h÷îng d¨n v  gióp ï tæi trong suèt qu¡ tr¼nh l m luªn v«n. T¡c gi£ xin gûi tîi c¡c Th¦y, cæ ð Vi»n To¡n håc, c¡c Th¦y, cæ khoa To¡n, pháng  o t¤o sau ¤i håc tr÷íng ¤i håc Khoa håc- ¤i håc Th¡i Nguy¶n công nh÷ c¡c Th¦y cæ tham gia gi£ng d¤y khâa Cao håc 2011 - 2013 líi c£m ìn s¥u s­c v· cæng lao gi£ng d¤y v  t¤o måi i·u ki»n thuªn lñi cho t¡c gi£ trong suèt qu¡ tr¼nh håc tªp t¤i tr÷íng. T¡c gi£ công xin ch¥n th nh c£m ìn Sð Gi¡o döc v   o t¤o t¿nh H  Giang, Ban gi¡m hi»u, c¡c tê chùc o n thº, tê chuy¶n mæn, nhâm To¡n tr÷íng THPT Li¶n Hi»p còng b¤n b± çng nghi»p v  gia ¼nh ¢ t¤o måi i·u ki»n gióp ï, ëng vi¶n t¡c gi£ ho n th nh luªn v«n n y. Th¡i Nguy¶n, ng y 01 th¡ng 8 n«m 2013 T¡c gi£ Ph¤m B¡ T¼nh Soá hoùa bôûi trung taâm hoïc lieäu http://www.lrc-tnu.edu.vn/ 5 Ch÷ìng 1 B€I TON VŠN TƒI C KHƒ N‹NG L×U THÆNG H„N CH˜ Ch÷ìng n y · cªp tîi b i to¡n vªn t£i câ kh£ n«ng l÷u thæng h¤n ch¸. N¶u mæ h¼nh v  c¡c t½nh ch§t cõa b i to¡n. Ti¸p â, n¶u ph÷ìng ph¡p t¼m ph÷ìng ¡n cüc bi¶n ban ¦u cõa b i to¡n. Sau â, tr¼nh b y cì sð lþ luªn v  nëi dung thuªt to¡n th¸ và (mët bi¸n thº cõa thuªt to¡n ìn h¼nh) gi£i b i to¡n. Cuèi ch÷ìng n¶u v½ dö sè minh håa thuªt to¡n gi£i. Nëi dung cõa ch÷ìng ÷ñc tham kh£o chõ y¸u tø c¡c t i li»u [1], [4] v  [5]. 1.1. B i to¡n v  t½nh ch§t Gi£ sû c¦n vªn chuyºn mët lo¤i h ng (xi m«ng ch¯ng h¤n) tø m kho chùa h ng (gåi l  c¡c tr¤m ph¡t) tîi n hë ti¶u thö (gåi l  c¡c tr¤m thu). ai > 0 (i = thu j l  bj > 0 Cho bi¸t l÷ñng h ng câ (gåi l  l÷ñng cung) ð tr¤m ph¡t i l  1, 2, ... , m) v  l÷ñng h ng c¦n (gåi l  l÷ñng c¦u) ð tr¤m (j = 1, 2, ... , n). Chi ph½ vªn chuyºn mët ìn và h ng tø tr¤m ph¡t i tîi tr¤m thu j l  cij ≥ 0. Ngo i ra, do i·u ki»n v· ÷íng s¡ ho°c h¤n ch¸ v· ph÷ìng ti»n vªn chuyºn n¶n tø i tîi j ch¿ ÷ñc vªn chuyºn tèi a dij ≥ 0 ìn và h ng (dij gåi l  kh£ n«ng l÷u thæng tr¶n tuy¸n i - j). V§n · l  c¦n vªn chuyºn bao nhi¶u ìn và h ng tø méi tr¤m ph¡t tîi méi tr¤m thu sao cho khæng v÷ñt qu¡ kh£ n«ng l÷u thæng tr¶n méi tuy¸n, måi tr¤m ph¡t giao h¸t h ng, måi tr¤m thu nhªn õ h ng v  têng chi ph½ vªn chuyºn l  nhä nh§t? Mæ h¼nh to¡n håc cõa b i to¡n vªn t£i vîi kh£ n«ng l÷u thæng h¤n Soá hoùa bôûi trung taâm hoïc lieäu http://www.lrc-tnu.edu.vn/ 6 ch¸ câ d¤ng nh÷ sau (xij l  l÷ñng h ng vªn chuyºn c¦n t¼m): m X n X cij xij → min (cüc tiºu têng chi ph½ c÷îc vªn chuyºn) (1.1) i=1 j=1 vîi c¡c i·u ki»n n X j=1 m X xij = ai , i = 1, 2, ..., m (måi xij = bj, j = 1, 2, ..., n (måi tr¤m ph¡t giao h¸t h ng) tr¤m thu nhªn õ h ng) (1.2) (1.3) i=1 0 ≤ xij ≤ dij , i = 1, 2, ..., m, j = 1, 2, ..., n (l÷ñng h ng vªn chuyºn khæng ¥m v  khæng v÷ñt qu¡ kh£ n«ng l÷u thæng (1.4) dij ) Gi£ thi¸t b i to¡n (1.1) - (1.4) thäa m¢n i·u ki»n c¥n b¬ng cung c¦u: m X ai = i=1 n X bj (têng cung b¬ng têng c¦u). (1.5) j=1 xij ≥ 0 suy ra xij ≤ min {ai , bj } vîi måi i, j n¶n d¹ d ng th§y r¬ng n¸u måi kh£ n«ng l÷u thæng dij ≥ min {ai , bj } th¼ ÷ìng nhi¶n câ thº bä i·u ki»n xij ≤ dij v  khi â (1.1) - (1.5) trð th nh Tø i·u ki»n (1.2), (1.3) v  b i to¡n vªn t£i thæng th÷íng (khæng h¤n ch¸ kh£ n«ng l÷u thæng). Ch¿ g°p khâ kh«n khi câ nhúng dij ≤ min {ai , bj }. Ch÷ìng n y s³ x²t b i to¡n vªn t£i trong tr÷íng hñp n y. R§t ti¸c khæng câ i·u ki»n c¦n v  õ º b i to¡n (1.1) - (1.5) gi£i ÷ñc, nh÷ trong b i to¡n vªn t£i thæng th÷íng, m  ch¿ câ i·u ki»n c¦n v  i·u ki»n õ ri¶ng. D¹ th§y r¬ng mët i·u ki»n õ º b i to¡n (1.1) - (1.5) gi£i ÷ñc l  dij ≥ min {ai , bj } vîi måi i,j. (1.6) º ÷a ra i·u ki»n c¦n cho (1.1) - (1.5) gi£i ÷ñc, khæng gi£m têng qu¡t ta câ thº gi£ thi¸t dij ≤ min {ai , bj } Soá hoùa bôûi trung taâm hoïc lieäu vîi måi i v  j, bði v¼ vîi http://www.lrc-tnu.edu.vn/ 7 dij ≥ min {ai , bj } th¼ r ng buëc xij ≤ dij trð n¶n khæng c¦n thi¸t. B¬ng c¡ch cëng (1.4) theo måi j v  so s¡nh vîi (1.2); sau â cëng (1.4) theo måi i v  so s¡nh vîi (1.3) ta nhªn ÷ñc mët i·u ki»n c¦n º b i to¡n (1.1) - (1.5) gi£i ÷ñc l  ai ≤ n X dij , ∀i = 1, ..., m v  bj ≤ j=1 m X dij , ∀j = 1, ..., n. (1.7) i=1 N¸u vi ph¤m (1.7) (tùc dij khæng õ lîn so vîi nhu c¦u vªn chuyºn) th¼ dò câ i·u ki»n c¥n b¬ng cung c¦u (1.5), b i to¡n v¨n s³ khæng câ ph÷ìng ¡n thäa m¢n (1.2) - (1.4), do â b i to¡n s³ khæng gi£i ÷ñc. Trong thüc h nh, º b i to¡n gi£i ÷ñc th÷íng ng÷íi ta ph£i i·u ch¿nh dij d¦n c¡c h» sè mët c¡ch th½ch hñp. Tuy nhi¶n, i·u ki»n (1.7) khæng l  i·u ki»n õ º (1.1) - (1.5) gi£i ÷ñc, nh÷ ch¿ ra ð v½ dö ìn gi£n sau. V½ dö 1.1 . X²t b i to¡n vªn t£i vîi m = 3 tr¤m ph¡t v  n = 3 tr¤m thu. L÷ñng cung cõa c¡c tr¤m ph¡t v  l÷ñng c¦u cõa c¡c tr¤m thu l¦n l÷ñt l  1, 3 v  4. Kh£ n«ng l÷u thæng dij ≤ min {ai , bj } (i, j = 1, 2, 3) ÷ñc cho trong B£ng 1.1. Rã r ng · (1.2) thäa m¢n vîi i = 2, 3 ph£i câ xij = dij vîi måi i = 2, 3 v  måi j = 1, 2, 3. Nh÷ng khi â s³ vi ph¤m r ng buëc (1.3) vîi j = 1. Vªy khæng thº câ xij thäa m¢n (1.1) - (1.4) vîi måi i, j = 1, 2, 3 v  b i to¡n (1.1) - (1.5) l  khæng gi£i ÷ñc. Kþ hi»u Aij ∈ Rm+n l  v²ctì h» sè cõa bi¸n xij . D¹ th§y r¬ng v²ctì n y câ hai th nh ph¦n b¬ng 1 t¤i h ng i v  h ng m + j, cán måi th nh ph¦n kh¡c b¬ng 0. Soá hoùa bôûi trung taâm hoïc lieäu http://www.lrc-tnu.edu.vn/ 8 º cho gån, ta ghi l¤i dú li»u cõa b i to¡n d÷îi d¤ng mët b£ng chú nhªt, gåi l  b£ng vªn t£i (B£ng 1.2). B£ng gçm m h ng (i = 1, ... , m) v  n cët (j = 1, ... , n). Ché giao nhau ð h ng i, cët j gåi l  æ (i, j). Méi h ng t÷ìng ùng vîi mët tr¤m ph¡t, méi cët t÷ìng ùng vîi mët tr¤m thu. Sè ghi ð ¦u méi h ng l  l÷ñng cung, sè ghi ð ¦u méi cët l  l÷ñng c¦u. Chi ph½ vªn chuyºn dij cij ghi ð gâc tr¶n b¶n tr¡i, kh£ n«ng l÷u thæng ghi ð gâc tr¶n b¶n ph£i, l÷ñng h ng vªn chuyºn xij s³ ghi ð giúa h ng d÷îi cõa æ (i, j). Æ (i, j) biºu thà tuy¸n vªn chuyºn tø tr¤m ph¡t i tîi tr¤m thu j. °t cij = ∞ ho°c dij = 0 n¸u khæng thº vªn chuyºn h ng tø i ¸n j. Sau ¥y nh­c l¤i mët sè kh¡i ni»m quen dòng. Ma trªn x = {xij }m×n ph÷ìng ¡n cõa b i to¡n (1.1) - (1.5). Mët ph÷ìng ¡n ¤t cüc tiºu (1.1) gåi l  ph÷ìng ¡n tèi ÷u hay líi gi£i. Ph÷ìng ¡n x gåi l  ph÷ìng ¡n cüc bi¶n n¸u c¡c v²ctì Aij ùng vîi xij thäa m¢n 0 < xij < dij l  ëc lªp tuy¸n t½nh hay tªp hñp æ G = {(i, j) : 0 < xij < dij } khæng t¤o th nh chu tr¼nh. Ph÷ìng ¡n cüc bi¶n x gåi l  khæng suy bi¸n n¸u G câ óng m + n - 1 æ, tr¡i l¤i x gåi l  suy bi¸n. thäa m¢n (1.2) - (1.4) gåi l  mët Soá hoùa bôûi trung taâm hoïc lieäu http://www.lrc-tnu.edu.vn/ 9 1.2. Ph÷ìng ¡n cüc bi¶n ban ¦u V½ dö 1.2 . T¼m ph÷ìng ¡n cüc bi¶n ban ¦u cho b i to¡n vªn t£i vîi dú li»u cho trong B£ng 1.3. • H ng thù nh§t: B­t ¦u tø æ (1, 1), ta ph¥n v o æ n y l÷ñng h ng: x11 = min {a1 , b1 , d11 } = min {40, 25, 15} = 15 = d11 (15 tæ ªm, m u ä). Ti¸p â, ph¥n v o æ (1, 2) l÷ñng h ng: x12 = min {a1 − 15, b2 , d12 } = min {25, 30, 20} = 20 = d12 (20 tæ ªm, m u ä). Cuèi còng, ph¥n v o æ (1, 3) l÷ñng h ng: x13 = min {a1 − 15 − 20, b3 , d13 } = min {5, 40, 25} = 5 < d13 (æ cì sð "•"). • • H ng thù hai: L¦n l÷ñt ph¥n h ng v o c¡c æ (2. 1), (2. 2) v  (2. 3). H ng thù ba: L¦n l÷ñt ph¥n h ng v o c¡c æ (3. 3) v  (3. 4). °t xij = 0 cho t§t c£ c¡c æ cán l¤i (khæng ÷ñc ph¥n phèi h ng). K¸t qu£ l  ta ÷ñc ph÷ìng ¡n cüc bi¶n (gçm 6 æ cì sð "•") ghi ð B£ng 1.3. Trong æ (i, j): sè ghi ð gâc tr¶n b¶n tr¡i l  ph£i l  dij v  sè ghi ð h ng d÷îi l  xij cij , sè ghi ð gâc tr¶n b¶n (sè in ªm m u ä ch¿ Æ cì sð ¡nh d§u "•". Soá hoùa bôûi trung taâm hoïc lieäu http://www.lrc-tnu.edu.vn/ xij = dij ). 10 1.3. Ti¶u chu©n tèi ÷u Gi£ sû  x0 = x0ij m×n l  mët ph÷ìng ¡n cüc bi¶n cõa b i to¡n (1.1) -  G = (i, j) : 0 ≤ x0ij ≤ dij gçm m + n 1 æ khæng chùa chu tr¼nh. Gi£ sû ui , vj l  c¡c th¸ và h ng, cët thäa m¢n h» ph÷ìng tr¼nh ui +vj = cij vîi måi (i, j) ∈ G. Kþ hi»u ∆ij = ui +vj −cij vîi måi (i, j) ∈ / G (∆ij = 0, ∀ (i, j) ∈ G). ành lþ sau cho ta mët d§u 0 hi»u nhªn bi¸t khi n o x l  ph÷ìng ¡n tèi ÷u.  0 0 (D§u hi»u tèi ÷u). Ph÷ìng ¡n x = xij m×n l  tèi ÷u cõa (1.1) - (1.5) khi ∆ij ≤ 0 vîi måi (i, j) ∈ / G, x0ij = 0 v  ∆ij ≥ 0 vîi måi (i, j) ∈ / G, x0ij = dij . (1.5) t÷ìng ùng vîi tªp æ cì sð ành lþ 1.1 1.4. Thuªt to¡n gi£i B÷îc 0 (Khði t¤o). X¥y düng ph÷ìng ¡n cüc bi¶n ban ¦u x0 = x0ij m×n . Tªp æ chån cì sð t÷ìng ùng vîi x0 l  G0 =  (i, j) : 0 ≤ x0ij ≤ dij gçm (m + n - 1) ph¦n tû v  khæng chùa chu tr¼nh, ngh¾a l  h» v²ctì {Aij : (i, j) ∈ G0 } ëc lªp tuy¸n t½nh. °t ch¿ sè váng l°p k = 0. B÷îc 1 ùng vîi . T½nh c¡c th¸ và Gk , ui (i = 1, ... , m) v  vj (j = 1, ... , n) t÷ìng b¬ng c¡ch gi£i h» ph÷ìng tr¼nh (d¤ng tam gi¡c): ui + vj = cij , ∀ (i, j) ∈ Gk . B÷îc 2 câ . T½nh c¡c ÷îc l÷ñng ∆ij = 0, ∀ (i, j) ∈ Gk ) B÷îc 3 ∆ij = ui + vj − cij , ∀ (i, j) ∈ / Gk (Ta luæn . . (Kiºm tra i·u ki»n tèi ÷u) ∆ij ≤ 0 vîi måi (i, j) ∈ / Gk , xkij = 0 v  ∆ij ≥ 0 vîi måi (i, j) ∈ / Gk , xkij = dij th¼ døng thuªt to¡n: xk l  ph÷ìng ¡n tèi ÷u (ành lþ 1.1). N¸u N¸u tr¡i l¤i, chuyºn sang B÷îc 4. B÷îc 4 . (i·u ch¿nh ph÷ìng ¡n) 4a) Æ i·u ch¿nh (r,s) l  æ ¤t max trong biºu thùc  ∆ = max ∆ij vîi (i, j) ∈ / Gk , xkij = 0, −∆ij Soá hoùa bôûi trung taâm hoïc lieäu vîi (i, j) ∈ / Gk , xkij = dij > 0. http://www.lrc-tnu.edu.vn/ 11 4b) Lªp chu tr¼nh C t¤o n¶n bði æ (r, s) vîi c¡c æ thuëc h¤n b¬ng c¡ch lo¤i d¦n c¡c æ treo cõa Gk ∪ {(r, s)} Gk (ch¯ng ). C1 (tªp æ l´) v  C2 (r, s) ∈ C2 n¸u xkrs = 4c) Ph¥n ho¤ch tªp C th nh hai tªp con ríi nhau (tªp æ ch®n) vîi qui ÷îc (r, s) ∈ C1 n¸u xkrs = 0 v  drs . 4d) X¡c ành l÷ñng i·u ch¿nh h theo cæng thùc  h = min dij − xkij (i, j) ∈ C1 , xkij vîi vîi (i, j) ∈ C2 ≥ 0. Æ (p, q) ¤t min tr÷íng hñp biºu thùc tr¶n s³ bà lo¤i. B÷îc 5 . X¥y düng ph÷ìng ¡n cüc bi¶n mîi theo cæng thùc xk+1 ij  k  xij + h khi (i, j) ∈ C1 , xkij − h khi (i, j) ∈ C2 , =  k xij khi (i, j) ∈ / C. v  tªp æ chån t÷ìng ùng T«ng k ←k+1 Gk+1 = (Gk \ {(p, q)}) ∪ {(r, s)}. v  trð l¤i B÷îc 1 thüc hi»n váng l°p k mîi. 1.5. V½ dö minh håa V½ dö 1.3 . Gi£i b i to¡n vªn t£i câ kh£ n«ng l÷u thæng h¤n ch¸ vîi v²ctì cung a, v²ctì c¦u b, ma trªn c÷îc ph½ n«ng l÷u thæng D = {dij }3×4 C = {cij }3×4 v  ma trªn kh£ nh÷ sau. a = (40, 50, 60)T , b = (25, 30, 40, 55)T     2 10 12 2 15 20 25 30 C =  5 1 2 10  , D =  20 25 30 30  20 3 10 15 15 25 35 55 Gi£i. V½ dö n y câ m = 3 tr¤m ph¡t, n = 4 tr¤m thu v  thäa m¢n i·u ki»n c¥n b¬ng cung c¦u (40 + 50 + 60 = 25 + 30 + 40 + 55). Xu§t ph¡t tø ph÷ìng ¡n cüc bi¶n ban ¦u G0 x0 cho ð B£ng 1.3 vîi tªp æ cì sð = (1, 3), (2, 1), (2, 2), (2, 3). (3, 3), (3, 4), gçm 6 æ khæng chùa chu tr¼nh v  gi¡ trà h m möc ti¶u Soá hoùa bôûi trung taâm hoïc lieäu  f x0 = 1285. http://www.lrc-tnu.edu.vn/ 12 Váng l°p 0 B÷îc 1 . . C¡c th¸ và t÷ìng ùng vîi x0 : u0 = (0, −10, −2) , v 0 = (15, 11, 12, 17) B÷îc 2 ∆ij = ui + vj − cij ghi ð ma trªn ∆0 d÷îi     15 20 5 0 13 1 0 15 x0 =  10 10 30 0  , ∆0 =  0 0 0 −3  0 0 5 55 −7 6 0 0 . C¡c sè B÷îc 3 x0 ch÷a tèi ÷u v¼ cán ∆14 = 15 > 0 vîi (1, 4) ∈ / = 6 > 0 vîi (3, 2) ∈ / G0 , x32 = 0 (hai sè d÷ìng g¤ch . Ph÷ìng ¡n G0 , x14 = 0 v  ch¥n trong ∆0 ). B÷îc 4 ¥y: ∆32 . ∆14 = max {∆14 , ∆32 } = 15 > 0. b) Chu tr¼nh i·u ch¿nh gçm 4 æ C = {(1, 4), (3, 4), (3, 3), (1, 3)} c) C¡c æ l´ C1 = {(1, 4), (3, 3)} v  c¡c æ ch®n C2 = {(3, 4), (1, 3)}.   d14 − x014 = 30, x034 = 55, d) L÷ñng i·u ch¿nh h = min = 5. d33 − x033 = 30, x013 = 5 a) Æ i·u ch¿nh (r, s) = (1, 4) vîi Æ lo¤i (p, q) = (1. 3). B÷îc 5 x1 (ghi ð G1 = {(1, 4) , (2, 1) , (2, 2) , (2, 3) . (3, 3) , (3, 4)}  f x0 = 1258. . Ph÷ìng ¡n cüc bi¶n mîi Váng l°p 1 B÷îc 1 d÷îi) vîi tªp æ cì sð v  f x1  = 1210 < . . C¡c th¸ và t÷ìng ùng vîi x1 : u1 = (0, 5, 13) , v 1 = (0, −4, −3, −3, 2) . B÷îc 2 ∆ij = ui + vj − cij ghi ð ma trªn ∆1 d÷îi ¥y:     15 20 0 5 −2 −14 −15 0 x1 =  10 10 30 0  , ∆1 =  0 0 0 −3  0 0 10 50 −7 6 0 0 . C¡c sè B÷îc 3 . Ph÷ìng ¡n G1 , x132 = 0 x1 (sè d÷ìng g¤ch ch¥n trong Soá hoùa bôûi trung taâm hoïc lieäu ∆32 = 6 > 0 ∆11 = −2 < 0 ch÷a tèi ÷u v¼ cán ∆1 ), http://www.lrc-tnu.edu.vn/ vîi vîi (3, 2) ∈ / (1, 1) ∈ / 13 G1 , x111 = d11 v  g¤ch ngang tr¶n B÷îc 4 ∆12 = −14 < 0 trong ∆1 ). vîi (1, 2) ∈ / G1 , x112 = d12 (hai sè ¥m . a) Æ i·u ch¿nh (r, s) = (1, 2) vîi: −∆12 = max {−∆11 = 2, −∆12 = 14, ∆32 = 6} = 14 > 0. b) Chu tr¼nh i·u ch¿nh gçm 6 æ: C = {(1, 2), (1, 4), (3, 4), (3, 3), (2, 3), (2, 2)} . c) C¡c æ C1 = {(1, 4), (3, 3), (2, 2)}, l´ c¡c æ ch®n C2 = = 15. Æ {(1, 2), (3, 4), (2, 3)}. d) L÷ñng i·u ch¿nh  h = min x112 = 20, d14 − x114 = 25, x134 = 50, d33 − x133 = 25, x123 = 30, d22 − x122 = 15  lo¤i (p, q) = (2. 2). B÷îc 5 x2 (ghi ð d÷îi) vîi tªp æ cì sð G2 = {(1, 2), (1, 4), (2, 1), (2, 3).(3, 3), (3, 4)} v  f (x2 ) = 1000 < f (x1 ) = . Ph÷ìng ¡n cüc bi¶n mîi 1210. Váng l°p 2 B÷îc 1 . . C¡c th¸ và t÷ìng ùng vîi x2 : u2 = (0, 5, 13) , v 2 = (0, 10, −3, 2) . B÷îc 2 ∆ij = ui + vj − cij ghi ð ma trªn ∆2 d÷îi ¥y:     15 5 0 20 −2 0 −15 0 x2 =  10 25 15 0  , ∆2 =  0 9 0 −3  0 0 25 35 −7 20 0 0 . C¡c sè B÷îc 3 G1 , x232 G2 , x211 x2 ∆32 = 20 > 0 vîi (3, 2) ∈ / = 0 (sè d÷ìng g¤ch ch¥n trong ∆2 ) , v  ∆11 = −2 < 0 vîi (1, 1) ∈ / = d11 (sè ¥m g¤ch ngang tr¶n trong ∆2 ). . Ph÷ìng ¡n Soá hoùa bôûi trung taâm hoïc lieäu ch÷a tèi ÷u v¼ cán http://www.lrc-tnu.edu.vn/ 14 B÷îc 4 . a) Æ i·u ch¿nh (r, s) = (3, 2) vîi: ∆32 = max {−∆11 = 2, ∆32 = 20} = 20 > 0. b) Chu tr¼nh i·u ch¿nh gçm 4 æ: C = {(3, 2), (1, 2), (1, 4), (3, 4)} . c) C¡c æ l´ C1 = {(3, 2), (1, 4)}, c¡c æ ch®n C2 = {(1, 2), (3, 4)}. d) L÷ñng i·u ch¿nh  h = min d32 − x232 = 25, x212 = 5, d14 − x214 = 10, x234 = 35 = 5. Æ lo¤i (p, q) = (1, 2). B÷îc 5 x3 (ghi ð d÷îi) vîi tªp æ cì G3 = {(1, 4), (2, 1), (2, 3), (3, 2).(3, 3), (3, 4)} v  f (x3 ) = 900 < f (x2 ) . Ph÷ìng ¡n cüc bi¶n mîi sð = 1000. Váng l°p 3. B÷îc 1 . C¡c th¸ và t÷ìng ùng vîi x3 : u3 = (0, 5, 13) , v 3 = (0, −10, −3, 2) . B÷îc 2 ∆ij = ui + vj − cij ghi ð ma trªn ∆3 d÷îi ¥y:     15 0 0 25 −2 −20 −15 0 x3 =  10 25 15 0  , ∆3 =  0 −6 0 −3  . 0 5 25 30 −7 0 0 0 . C¡c sè B÷îc 3 x3 ch÷a tèi ÷u v¼ cán ∆11 = −2 < 0 vîi (1, 2) ∈ / 3 = −6 < 0 vîi (2, 2) ∈ / G3 , x22 = d22 (hai sè ¥m g¤ch . Ph÷ìng ¡n G3 , x311 = d11 ngang tr¶n B÷îc 4 ∆22 trong ∆3 ). v  . a) Æ i·u ch¿nh (r, s) = (2, 2) vîi: −∆22 = max {−∆11 = 2, −∆22 = 6} = 6. C = {(2, 2), (2, 3), (3, 3), (3, 2)}. C1 = {(2, 3), (3, 2)}, c¡c æ ch®n C2 = {(2, 2), (3, 3)}. b) Chu tr¼nh i·u ch¿nh gçm 6 æ c) C¡c æ l´ Soá hoùa bôûi trung taâm hoïc lieäu http://www.lrc-tnu.edu.vn/ 15 d) L÷ñng i·u ch¿nh:  h = min x322 = 25, d23 − x323 = 15, x333 = 25, d32 − x332 = 20 = 15. Æ lo¤i (p, q) = (2. 3). B÷îc 5 x4 (ghi ð d÷îi) G4 = {(1, 4), (2, 1), (2, 2), (3, 2).(3, 3), (3, 4)} v  f (x4 ) = . Ph÷ìng ¡n cüc bi¶n mîi vîi tªp æ cì sð 810 < f (x3 ) = 900. Váng l°p 4. B÷îc 1 . C¡c th¸ và t÷ìng ùng vîi x4 : u4 = (0, 11, 13) , v 4 = (−6, −10, −3, 2) B÷îc 2 ∆ij = ui + vj − cij ghi ð ma trªn ∆4 d÷îi ¥y:     15 0 0 25 −8 −20 −15 0 x4 =  10 10 30 0  , ∆4 =  0 0 6 −3  . 0 20 10 30 −13 0 0 0 . C¡c sè B÷îc 3 . Ph÷ìng ¡n G4 , x411 = d11 B÷îc 4 x4 ∆11 = −8 < 0 vîi (1, 1) ∈ / ∆4 ). ch÷a tèi ÷u v¼ cán (sè ¥m g¤ch ngang tr¶n trong . a) Æ i·u ch¿nh (r, s) = (1, 1) vîi −∆11 = 8 > 0 b) Chu tr¼nh i·u ch¿nh gçm 6 æ: C = {(1, 1), (1, 4), (3, 4), (3, 2), (2, 2), (2, 1)} . c) C¡c æ l´ C1 = {(1, 4), (3, 2), (2, 1)}, c¡c æ ch®n C2 = {(1, 1), (3, 4), (2, 2)}. d) L÷ñng i·u ch¿nh  h = min x411 = 15, d14 − x414 = 5, x434 = 30, d32 − x432 = 5, x422 = 10, d21 − x421 = 10  = 5. Æ lo¤i (p, q) = (1. 4). B÷îc 5 x5 (ghi ð d÷îi) vîi tªp æ cì G5 = {(1, 1), (2, 1), (2, 2), (3, 2), (3, 3), (3, 4)} v  f (x5 ) = 770 < f (x4 ) . Ph÷ìng ¡n cüc bi¶n mîi 810. Soá hoùa bôûi trung taâm hoïc lieäu http://www.lrc-tnu.edu.vn/ sð = 16 Váng l°p 5. B÷îc 1 . C¡c th¸ và t÷ìng ùng vîi x5 : u5 = (0, 3, 5) , v 5 = (2, −2, 5, 10) B÷îc 2 ∆ij = ui + vj − cij ghi ð ma trªn ∆5 d÷îi ¥y:     10 0 0 30 0 −12 −7 8 x5 =  15 5 30 0  , ∆5 =  0 0 6 3 . 0 25 10 25 −13 0 0 0 . C¡c sè B÷îc 3 . Ph÷ìng ¡n G5 , x524 = 0 B÷îc 4 x5 ch÷a tèi ÷u v¼ cán (sè d÷ìng g¤ch ch¥n trong ∆5 ∆24 = 3 > 0 vîi (2, 4) ∈ / ). . ∆24 = 3 > 0. b) Chu tr¼nh i·u ch¿nh gçm 4 æ C = {(2, 4), (3, 4), (3, 2), (2, 2)}. c) C¡c æ l´ C1 = {(2, 4), (3, 2)}, c¡c æ ch®n C2 = {(3, 4), (2, 2)}. a) Æ i·u ch¿nh (r, s) = (2, 4) vîi d) L÷ñng i·u ch¿nh  h = min d24 − x524 = 30, x534 = 25, d32 − x532 = 0, x522 = 5 = 0. Æ lo¤i (p, q) = (3, 2). B÷îc 5 x6 (ghi ð d÷îi) vîi tªp æ cì G6 = {(1, 1), (2, 1), (2, 2), (2, 4).(3, 3), (3, 4)} v  f (x6 ) = f (x5 ) = 770. . Ph÷ìng ¡n cüc bi¶n mîi Váng l°p 6 B÷îc 1 . C¡c th¸ và t÷ìng ùng vîi sð x6 : u6 = (0, 3, 8) , v 6 = (2, −2, 2, 7) B÷îc 2 ∆ij = ui + vj − cij ghi ð ma trªn ∆6 d÷îi ¥y:     10 0 0 30 0 −12 −10 5 x6 =  15 5 30 0  , ∆6 =  0 0 3 0 . 0 25 10 25 −10 3 0 0 . C¡c sè B÷îc 3 ∆ij ≤ 0 vîi måi (i, j) ∈ / G6 , x6ij = 0 (c¡c sè ¥m g¤ch ch¥n trong ∆6 ), ∆ij ≥ 0 vîi måi (i, j) ∈ / G6 , x6ij = dij (c¡c sè d÷ìng g¤ch ngang tr¶n trong ∆6 ).  opt Vªy x = x6 , fopt = f x6 = 770.  . Ph÷ìng ¡n x6 Soá hoùa bôûi trung taâm hoïc lieäu tèi ÷u v¼ http://www.lrc-tnu.edu.vn/ 17 Tâm l¤i, ch÷ìng n y ¢ tr¼nh b y tâm t­t c¡c k¸t qu£ v· b i to¡n vªn t£i câ kh£ n«ng l÷u thæng h¤n ch¸: mæ h¼nh v  c¡c t½nh ch§t cõa b i to¡n, ph÷ìng ph¡p t¼m ph÷ìng ¡n cüc bi¶n ban ¦u cõa b i to¡n, ti¶u chu©n tèi ÷u v  thuªt to¡n th¸ và gi£i b i to¡n. Cuèi ch÷ìng, n¶u v½ dö sè minh håa thuªt to¡n gi£i. Nâi chung, thuªt to¡n gi£i b i to¡n vªn t£i câ h¤n ch¸ kh£ n«ng l÷u thæng phùc t¤p hìn æi chót so vîi b i to¡n vªn t£i khæng h¤n ch¸ kh£ n«ng l÷u thæng. Soá hoùa bôûi trung taâm hoïc lieäu http://www.lrc-tnu.edu.vn/ 18 Ch÷ìng 2 B€I TON VŠN TƒI VÎI L×ÑNG CUNG BÀ CHN D×ÎI Ch÷ìng n y tr¼nh b y mæ h¼nh to¡n håc cõa b i to¡n vªn t£i, trong â l÷ñng cung cõa c¡c tr¤m ph¡t khæng qui ành tr÷îc m  bà ch°n tr¶n v  d÷îi, cán l÷ñng c¦u cõa c¡c tr¤m thu ¢ bi¸t tr÷îc. Sau â x²t tr÷íng hñp ìn gi£n khi l÷ñng cung cõa c¡c tr¤m ph¡t ch¿ bà ch°n d÷îi bði sè khæng ¥m. N¶u c¡ch ÷a v· b i to¡n g¦n vîi b i to¡n vªn t£i thæng th÷íng v  x²t t½nh ch§t nghi»m cõa b i to¡n. Sau â tr¼nh b y ti¶u chu©n tèi ÷u v  thuªt to¡n gi£i. Cuèi ch÷ìng n¶u v½ dö sè minh håa. Nëi dung cõa ch÷ìng ÷ñc tham kh£o chõ y¸u tø c¡c t i li»u [1], [3] v  [5]. 2.1. Nëi dung b i to¡n º thuªn ti»n cho vi»c kh£o s¡t, ð ¥y ta s³ mæ t£ b i to¡n vªn t£i vîi l÷ñng cung cõa c¡c tr¤m ph¡t bà ch°n (tr¶n v  d÷îi) nh÷ sau. Gi£ sû câ mët lo¤i h ng (ch¯ng h¤n l÷ìng thüc, xi m«ng, v.v ...) c¦n vªn chuyºn tø m iºm cung c§p (gåi l  c¡c tr¤m ph¡t), kþ hi»u i = 1, 2, ... , m, tîi n iºm ti¶u thö (gåi l  c¡c tr¤m thu), j = 1, 2, ... , n. Cho bi¸t kh£ n«ng cung c§p h ng cõa tr¤m ph¡t i thuëc kho£ng cho tr÷îc ai l  mùc (0 ≤ ai ≤ ai ) ai Trong â cung tèi thiºu v  ph¡t i , nhu c¦u ti¶u thö h ng cõa tr¤m thu j l  [ai , ai ]. l  mùc cung tèi a cõa tr¤m bj > 0 (bj cè ành) v  c÷îc ph½ vªn chuyºn mët ìn và h ng tø tr¤m ph¡t i tîi tr¤m thu j l  cij ≥ 0. B i to¡n °t ra l  h¢y t¼m mët ph÷ìng ¡n vªn chuyºn h ng tø c¡c tr¤m ph¡t tîi c¡c tr¤m thu sao cho têng chi ph½ vªn chuyºn l  nhä nh§t? Soá hoùa bôûi trung taâm hoïc lieäu http://www.lrc-tnu.edu.vn/ 19 º £m b£o cho b i to¡n gi£i ÷ñc, ta gi£ thi¸t m X ai ≤ i=1 n X bj ≤ j=1 m X ai (2.1) i=1 B i to¡n tr¶n câ thº di¹n t£ d÷îi d¤ng mæ h¼nh to¡n håc nh÷ sau: m X n X f (x) ≡ (P ) cij xij → min, (2.2) xij − xi0 = ai , i = 1, 2, ..., m, (2.3) xij = bj , j = 1, 2, ...,n, (2.4) i=1 j=1 n X j=1 m X i=1 trong â thu j, xij xi0 xij ≥ 0, i = 1, ..., m, j = 1, ..., n, (2.5) 0 ≤ xi0 ≤ ei ≡ ai − ai , i = 1, 2, ..., m, (2.6) biºu thà l÷ñng h ng c¦n vªn chuyºn tø tr¤m ph¡t i tîi tr¤m l  têng l÷ñng h ng chuyºn i tø tr¤m ph¡t i (tîi måi tr¤m ai ≤ ai + xi0 ≤ ai , i = 1, 2, ..., m. V¼ th¸, i·u ki»n (2.3) câ thº vi¸t l¤i th nh ai ≤ xi1 +...+xin ≤ ai vîi måi i. Do â, (P) cán ÷ñc gåi l  b i to¡n vªn t£i vîi r ng buëc hai ph½a v· l÷ñng cung cõa c¡c tr¤m ph¡t. thu) v÷ñt mùc cung tèi thiºu ai . Tø (2.6) suy ra Kþ hi»u s1 = a1 + a2 + ... + am (têng cung s2 = a1 + a2 + ... + am (têng cung d = b1 + b2 + ... + bn (têng c¦u). tèi thiºu), tèi a), Trong [2] c¡c t¡c gi£ ¢ sû döng mæ h¼nh tr¶n º mæ t£ b i to¡n ph¥n bê tèi ÷u c¡c trung t¥m i·u h nh (depot), i = 1, 2, ... , m, qu£n lþ c¡c tuy¸n xe buþt, j = 1, 2, ... , n, nh¬m gi£m ¸n mùc th§p nh§t têng qu¢ng ÷íng xe ch¤y khæng t£i (bi¸n t¥m i qu£n lþ v  xij = 0 xij = 1 n¸u tuy¸n xe j do trung n¸u tr¡i l¤i), nhí â cho ph²p ti¸t ki»m kho£ng 10% chi ph½ vªn h nh, trong â ai v  ai l  giîi h¤n d÷îi v  giîi h¤n tr¶n cho sè tuy¸n xe buþt m  trung t¥m i ÷ñc ph²p qu£n lþ, i·u ki»n (2.4) Soá hoùa bôûi trung taâm hoïc lieäu http://www.lrc-tnu.edu.vn/
- Xem thêm -

Tài liệu liên quan