Computation of order conditions for symplectic - HAL

Computation of order conditions for symplectic
partitioned Runge-Kutta schemes with application to
optimal control
J. Frederic Bonnans, Julien Laurent-Varin
To cite this version:
J. Frederic Bonnans, Julien Laurent-Varin. Computation of order conditions for symplectic
partitioned Runge-Kutta schemes with application to optimal control. [Research Report] RR5398, INRIA. 2004, pp.18. <inria-00070605>
HAL Id: inria-00070605
https://hal.inria.fr/inria-00070605
Submitted on 19 May 2006
HAL is a multi-disciplinary open access
archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.
L’archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d’enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.
INSTITUT NATIONAL DE RECHERCHE EN INFORMATIQUE ET EN AUTOMATIQUE
Computation of order conditions for symplectic
partitioned Runge-Kutta schemes with application to
optimal control
J. Frédéric BONNANS — Julien LAURENT-VARIN
N° 5398
Décembre 2004
ISSN 0249-6399
ISRN INRIA/RR--5398--FR+ENG
Thème NUM
apport
de recherche
!"#$!&%'()*
+)-,./102'&+!345)!76893:(;:
<5(=;>?(
@+ACBCDFEHGEIDKJMLONQPSRRTRU V @XWYMJ[ZH\^])T'_`
aRbdcfegT`
hR
ijlk<monprqtsvuxwzy|{~}fk<monF{)z€lmoF‚„ƒ†…H€lnF{
‡C‚fˆŠ‰~n<}twH‹gŒŽ)
4’‘l‘Xˆ’‚„})“ln‚„nK” jln<‚ ” jln•X–—’˜’™Iš›uœŒrK”n<m>l‚fnž’Ÿ’Ÿ’ u-¡K𑢐У’nK{
∗
†
¤
¥+¦F§<¨#©lª’§’«S¬ nQ“|ƒ­{f”€¢{„{t}„j1no“ln<‚fƒ¯®##}fƒ¯ˆI›ˆ’°*ˆ’‚ “|nF‚”<ˆ’1“lƒ±}fƒ¯ˆI1{t°²ˆ’‚g}„jln'“lƒ­{f”‚fn}fƒ¯³K#}„ƒ­ˆ’´ˆŠ°gµ¶€l1”<ˆ’|·
{¸}„‚ Šƒ­lnK“1¹ˆ’‘l}„ƒ­mo’º;”<ˆ’H}„‚fˆ’º‘l‚fˆ’lº­n<mQ{F»+¼)j1n<½}fjln{f” jlnFmdn
°²ˆ’‚}„jln{~} #}fn
nF…H€1Š}„ƒ­ˆ’=ƒ†{ˆ’°r)€ll£In·
¾€l}¸}fS}~yH‘Xn’¿Oij1ƒ­{>‘l‚fˆ’lº­n<m2Š‘l‘XnF’‚f{}fˆ´Xn'nF{f{„n<H}„ƒ†Šº­º¯y›}„j1n'ˆ’ln
ˆŠ°4” jlnF” Àzƒ­l£OˆI‚f“|nF‚”ˆI1“|ƒ¯}„ƒ­ˆ’1{
°²ˆI‚{„yzmd‘1º¯nK”9}„ƒ†”)‘1’‚¸}fƒ±}fƒ¯ˆIlnF“Q4€ll£’n<·Á¾€l}¸}f{„” jlnFmonF{F¿ ¬ nt{„jlˆ#¼Â}„j1Š}*}fjln4}fjlnr”ˆImo‘l€|}fŠ}„ƒ­ˆ’1{*€1{„ƒ¯l£
1ƒ±·Ã”ˆIº¯ˆI€l‚„nK“Ä}„‚fn<nK{Š‚fn'1Š}„€l‚ Šº­º­yÄn<Å|‘l‚„nK{„{„nF“ƃ­=}fjlƒ­{o”F’{„n'ƒ¯Ç}„nF‚„mQ{>ˆŠ°>È#É9ʶË<Ì1ÍÃË ÎdÏ É„Ë ËÍMÉfËfË<¿›ij1ƒ­{
£Iƒ¯®InF{d¼;Ðy'}„ˆQ”ˆImd‘1€|}„n}fjln<m$HyŠŠ‘l‘l‚fˆ’‘1‚„ƒ†#}fng”<ˆ’mo‘l€|}fn<‚)‘l‚fˆ’£I‚f’m¿
t€l‚g{¸ˆ’°Ñ}~¼;’‚„nƒ†{r’lº¯n}„ˆ
”<ˆ’mo‘l€|}fn”ˆI1“|ƒ¯}„ƒ­ˆ’1{r€l‘}„ˆ'ˆI‚f“ln<‚gÒµ¶¼n>“|ƒ­{„‘lº†Ðy'}fjln<mӀl‘}fˆ'ˆ’‚ “|n<‚
Ô ¹9¿Õijln‚fnF{„€lº±} {Š‚fnrƒ­I”<”ˆI‚f“1Š1”<nr¼)ƒ¯}„j}„jlˆH{¸ngˆŠ°ÖrŠ£In<‚µ²¼)jlnF‚„nt}„j1n<y'¼;n<‚fng”<ˆ’mo‘l€|}fnF“'°²ˆ’‚ˆ’‚ “|n<‚
€1‘}fˆo H¹’{;¼;n<º­º?’{;}fjlˆI{„nˆŠ°sS€l‚f€1d¼)jln<‚fng}„jlnz€lm>¢nF‚4ˆŠ°”ˆ’¢“|ƒ±}fƒ¯ˆI1{€l‘}„ˆQˆI‚f“ln<‚r҃­{){¸}fŠ}„nF“&¿
×SØHÙÛÚfÜÝ ¨#Þ+¦’« t‘|}„ƒ­mQŠº*”<ˆ’H}„‚fˆ’º[»?ÖtŠmoƒ¯º¯}„ˆIlƒ†ŠÆ{¸y|{¸}„n<mQ{F»?‘1’‚¸}fƒ±}fƒ¯ˆIlnF“ď)€l1£’n·Ã¾€|}„}fmon}„j1ˆz“1{<»
w|yHmo‘lº­nF”}„ƒ†”{„” j1n<monF{F»|ˆ’‚ “|n<‚)”<ˆ’1“|ƒ¯}„ƒ­ˆ’¢{<»|‡Õ·Ã{„n<‚fƒ¯nK{<»lÖ4·M}f‚„nFnF{F»lt‚fƒ­n<H}„nK“'°²‚fn<n<·M}f‚„nFnF{F¿
ò åKä[àKé9ß&ê¸àÐþ áãâ¢âMäMåÐæKç)áãâ¢â[åÐèÐèŠé9ê²ä[ë~æ4ìç4íïî+ðXñò óÐærôî+ð¢õ?ö*÷áãó;ä[àKëø¯ê[ò ù)ëÃúïé9ê¶û;é ø’ä[àKëÕíïî+ðXñ*ø­ë~üãüãé¸úâMàÐáãèé ø’ä[àKë+â[ë~ý~é9óÐæ
íïàKë¸âMó#ÿ ò¸åÐþ çüã#êáã$÷ ë~#êó+¸æMò"z¸óÐê¶ö-áãý~ý,+ë&õ?%&#êôðX¶îîë~æÐîßë~.öêMláãý9î+ö+þ Ûñõ-é9²óÐ1î/óÐêMòé"zóÐ
¶ëÃâ(é9ä*'+ó1ñ<0áãóÐÑçFß|ê¶æÐá ë~é<òKê¶ý¸þ ù2ø¯éÐê*)Á÷þ+ÑóÐë~ê¶â[á á3ò09²ó+õ&4éFý<ñ<åÐçKë~â¶óÐä[ë~ý~ù5é9åÐê¶Ñä¸ó÷ä[ë+09é9êMòfù4ä[òáãé9áãóKó6ëÕ+æÐë~ëè#¢ò ê¶é9äMüãù)åÐý~ë~ë¸ó<ò ä¸åH÷l÷ôîðXõ&öK÷÷Iò óÐæ7!#¢ê¶"Héë ÿ ÷ ²î-õ Ñö•õ&éFý <åÐë~óÐý~é9åÐê²ä„9÷ +é9ù4ò áãóÐëdæK:ë ¢é9üãåÐý~ë„ò åI8÷ ;K8÷ !<"HëQíÛàÐë~â[ó#ò~=ç % ÿ åÐüãáãë~óHþ "zò åÐêMë~ó>ä MëÃäñ<çFæÐé<ý~éÐ8
Ûò ê¶áã
ó '+áãóÐê¶á òKþ ø­?ê )Áþ
∗
†
Unité de recherche INRIA Rocquencourt
Domaine de Voluceau, Rocquencourt, BP 105, 78153 Le Chesnay Cedex (France)
Téléphone : +33 1 39 63 55 11 — Télécopie : +33 1 39 63 53 30
(;(½)!7; <!$¸
!?345!7
,$./l02
&5!&%'()C+o)!" +)!$)v(;
( ; 5 +<5(
¦ « Š‚„}„ƒ†”º­nÕ}„‚ Šƒ¯}„nº¯n*”FŠº†”€lºz“|nK{?”ˆI1“|ƒ¯}„ƒ­ˆ’1{“ ˆ’‚ “|‚„n*‘¢ˆI€l‚?º†r“|ƒ†{f”‚f}„ƒ†{f#}„ƒ­ˆ’“|n‘l‚fˆ’lº­k<monF{
“ln
”ˆImdmQ’1“|nQˆ’‘l}„ƒ­mo’º¯n{„’1{”ˆIH}„‚ Šƒ­I}fnF{F»?º­ˆ’‚ {„…H€lnQº­n
{f” jl<mQnK{~}>“lnQ}~yz‘¢n
4€ll£’n<·Á¾€l}¸}f1¿SŽn
‘1‚„ˆIlº¯kFmond{„no‚„K“|€lƒ¯} ”<n<º­€lƒC“|nQº†
®I<‚fƒ ¢”F#}fƒ¯ˆI›“|nQ”ˆI1“|ƒ¯}„ƒ­ˆ’1{“ ˆI‚f“|‚fno“|nomd<}„jlˆ|“|nK{“|nQ)€ll£In·
¾€l}¸}f{„yzmd‘1º¯nK”9}„ƒ†…H€lnF{+‘1’‚¸}fƒ±}fƒ¯ˆIllFnF{F¿p4ˆI€1{moˆ’H}„‚fˆ’¢{+…H€lnº­n)”<Šº†”€1º|1’{„;{„€l‚º­nF{ՐŠ‚fl‚fnF{+lƒ†”ˆIº¯ˆI‚„nK{
Ð®InF”‚ ’”<ƒ¯ln{ nÅ|‘l‚fƒ¯mon;“|n*mQŠ1ƒ¯kF‚„nC¢#}„€1‚„nFº¯º­nCn<€|}„ƒ­º¯ƒ†{fŠH}+“|nK{’‚„1‚„nK{ﺯƒ­l‚fnF{ˆ’‚fƒ¯nFH}„F{F¿+ŽnK”ƒH‘¢nF‚„mon}
“ n<Åz‘1‚„ƒ­mon<‚4”nK{)”ˆI1“|ƒ¯}„ƒ­ˆ’1{‘1’‚)€l‘l‚„ˆI£’‚ Šmomon’¿
prˆŠ}„‚fn>”<ˆ|“|n”<’º­”<€lº­nº¯nK{t”<ˆ’1“lƒ±}fƒ¯ˆI1{Չ~€1{f…H€ Qº ˆ’‚ “|‚fndÒµ¶ƒ¯º†{t{„ˆ’H}r‘1‚„K{¸nFI}fF{C‰~€1{„…H€ 'º ˆ’‚ “|‚„n Ô ¹9¿
nF{;‚fF{„€lº±} #} {;{„ˆ’H}n<’”F”ˆ’‚ “
Ð®InF”t”n<€lœ|nÖrŠ£In<‚µ¶…H€lƒ&o”FŠº†”€lº­tº­nF{”<ˆ’1“|ƒ¯}„ƒ­ˆ’¢{+‰~€1{„…H€ >º ˆ’‚ “|‚fn
z¹n}rsS€l‚f€1µ¶…H€lƒd“l<lˆIml‚f”nF{)”<ˆ’1“lƒ±}fƒ¯ˆI1{+‰~€1{f…I€ dº ˆI‚f“|‚fnҊ¹¿
Ý §K¦ Ú ª ¦S« Žˆ’momQŠ1“|nˆ’‘|}fƒ¯mQŠº­n’»Û{„yz{¸}„kFmonF{tj1Šmoƒ­º±}fˆ’lƒ­n<1{F»Û{„” jlFmQ’{r“ln>)€l1£’n·Ã¾€|}„}fQ‘1’‚¸·
}fƒ±}fƒ¯ˆIllK{<»1{f” jlFmoI{{¸yzmo‘lº­nF”9}fƒ­…H€lnK{<»l”<ˆ’1“lƒ±}fƒ¯ˆI1{)“ ˆ’‚ “|‚„nI»|‡Õ·Ã{¸F‚„ƒ­nF{F»|Ö)·ÃŠ‚fl‚„nK{<»|’‚„l‚fnF{;ˆI‚„ƒ­n<H}„K{<¿
ɄÎIËÉ È#ÌÛÎ#ÊÑÍ[ʶÈ#Ì*ÏÈ#É
¯Ë<ÍMÊ
ŠÉ9ÍMÊÑÍ[ʶÈ#ÌÛË Î|ÌHË|Í[ÍlËQË
˜
> ?Q*
ij1n
moˆŠ}fƒ¯®##}fƒ¯ˆI=ˆŠ°)}fjlƒ­{>¼;ˆ’‚fÀĔ<ˆ’monF{°²‚„ˆIm ’½Š1’º¯y|{„ƒ­{zyÄÖtŠ£’nF‚#" Ò%$ˆŠ°4ˆI‚f“|nF‚”ˆI1“|ƒ¯}„ƒ­ˆ’1{°²ˆ’‚
ˆI‘|}„ƒ­mQŠºX”ˆ’H}f‚„ˆIºl‘l‚fˆ’lº­n<mQ{rµ¶ˆŠ°&ŠQˆI‚f“lƒ¯1’‚„yd“|ƒ'&XnF‚„nFI}fƒ­’º1nK…I€¢#}„ƒ­ˆ’X¹9¿Õijln4ƒ†“|nFƒ­{C}„ˆ>“lƒ­{f”‚fn}fƒ¯³Fn)}„jln
{¸}fŠ}„n>nF…H€1#}fƒ¯ˆIOzyS
)€l1£’n·Ã¾€|}„}f
{f” jlnFmdnI»Û¼)ƒ±}fj›
“lƒ(&Ûn<‚fn<H}t®#’º¯€lndˆŠ°*”<ˆ’H}„‚fˆ’º+’{f{¸ˆ|”<ƒ­Š}„nF“¼)ƒ¯}„j
nK’” #j )¸ƒ¯1ln<‚g{~} #}„n * ¿tÖrŠ£In<‚4ˆI1{„n<‚f®’nF{)}„j1Š}r}„j1n‚fnF{„€lº±}fƒ¯1£'ˆ’‘|}fƒ¯mQ’º¯ƒ¯}~yS{¸y|{¸}„n<m»Û#°Ñ}fn<‚t{„ˆ’mon” j¢Šl£In
ˆ’°&®#Š‚fƒ­’lº¯nI»Iƒ†{;‘1’‚¸}fƒ±}fƒ¯ˆIlnF“'4€ll£’n<·Á¾€l}¸}f>{„” j1n<mon’¿Ö4n”ˆImd‘1€|}„nK{C}fjln<ĵ²zyQj1Š¢“ï»zƒM¿ n’¿­»z¼)ƒ±}fjlˆ’€|}
”<ˆ’mo‘l€|}fn<‚4”<ˆz“lnK¹*}„jlnˆI‚f“ln<‚)”ˆI1“|ƒ¯}„ƒ­ˆ’1{°²ˆ’‚)ˆI‚f“ln<‚€l‘}„ˆo 1¿CwznFnŠº†{¸ˆd}„j1ng‚fnF{„€lº¯}f{ˆ’+° "ã,— $?Š¢-“ " Ô $&ˆI
”<ˆ’1{¸}„‚ Šƒ­lnK“´ˆI‘|}„ƒ­mQŠº*”<ˆ’H}„‚fˆ’ºÕ‘l‚fˆ’1º¯nFmo{QµM 1‚ {~}ˆI‚f“|nF‚{~} #}fn'”ˆ’¢{~}f‚f’ƒ¯lnK“´‘1‚„ˆIlº¯nFm»“|ƒ†{„”<‚„n<}„ƒ­³<nK“
z/y .*€lº¯nF‚ {4{f” jln<monI»lŠ1“o”ˆ’H}f‚„ˆIºÛ”<ˆ’1{¸}„‚ Šƒ­lnK“
‘l‚fˆ’1º¯nFm8¼)ƒ±}fjSd4€ll£’n<·Á¾€l}¸}fo“|ƒ†{f”‚fn}„ƒ­³FŠ}„ƒ­ˆ’¢¹¿
ijlnF‚„nŠ‚fnÕnF{f{¸nFI}fƒ­’º¯º­y4}~¼;ˆ)jzyz‘¢ˆ’}„jlnK{¸nK{&ƒ­}„jln;Š1’º¯y|{„ƒ­{ۈ’0° " %Ò $Á»Ðˆ’ln*ˆ’}„j1nCˆ’‚fƒ¯£Iƒ¯¢ŠºŠ‘l‚fˆ’1º¯nFm Š¢“
}fjln)ˆŠ}fjln<‚*ƒ­{C‚fnF{¸}„‚fƒ†”9}„ƒ­ˆ’dˆ’d}fjlnr{f” jln<mon’¿Ctln)j1I{+}„ˆI{„{„€lmon;}„j1Š}C}„j1n4ÖrŠmoƒ­º±}fˆ’lƒ†Šoƒ†{C{¸}„‚fˆ’l£Iº¯y
”<ˆ’z®’n<ż¿ ‚K¿ }K¿•}„j1nĔ<ˆ’H}„‚fˆ’º[»)ˆ’‚moˆ’‚fnO£’nFln<‚ Šº­º¯y }„j¢#}
}fjlnÆ{¸nK”ˆ’¢“ “|n<‚fƒ­®ÐŠ}„ƒ­®’nOˆ’°Ör’mdƒ­º¯}„ˆ’1ƒ­’
¼¿ ‚K¿ }K¿}„jlno”ˆII}f‚„ˆIº?ƒ†{gƒ¯z®’nF‚¸}fƒ¯1º¯nI¿ij1ƒ­{Šº­º­ˆ#¼4{4}fˆnFº¯ƒ­moƒ¯1Š}„n>}„jlno”ˆIH}„‚fˆ’º}„j1’lÀ|{4}fˆ
}fjln>ƒ­mo‘lº­ƒ­”<ƒ±}
}fjln<ˆI‚„nFm»I{„ˆ}„j1Š}*¼n)j¢Ð®’n)ŠQnK…I€1ƒ¯®#Šº­n<H}*{f” jln<mon°²ˆI‚Õ}„jlnr‚„nK“|€1”nK“Sµ¶{¸}fŠ}„n’»H’“Љ~ˆIƒ¯H}*{¸}fŠ}„nK¹C{¸y|{¸}„n<m¿
ij1no{„nF”<ˆ’1“›jzyz‘¢ˆ’}„jlnK{¸ƒ†{ƒ­{}„j1Š}1ˆ’lnQˆ’°*}„j1n'”ˆzn 1'”ƒ­n<H} { b {'µ²ˆ’°}fjlnQ‘1’‚¸}fƒ­”<€lº­’‚)€ll£In·Ã¾€|}¸} 
{f” jlnFmdnйƒ­{)³Fn<‚fˆ1¿
ijlnomQŠƒ­ÄnK{¸€lº¯}ˆŠ°ÖrŠ£In<2‚ " %Ò $Cƒ†{g}„j1Š}F»&ƒ¯°}fjlnoˆ’‚fƒ¯£Iƒ¯¢ŠºÕ)€ll£In·Ã¾€|}¸} {f” jlnFmdnoƒ†{ˆŠ°µ²£’º­ˆ’¢ŠºÑ¹
ˆI‚f“ln<‚ p µ²ƒ[¿ n’¿­»#¼)jln<dŠ‘l‘lº­ƒ­nF“}fˆg’€l1”<ˆ’H}„‚fˆ’º­º­nF““|'ƒ &Ûn<‚fn<H}„ƒ†ŠºHnF…H€1Š}„ƒ­ˆ’¢¹&}„jlnF}fjln‚fnF{„€lº¯}„ƒ­l£t{f” jlnFmdn
j¢’{4ˆI‚f“ln<‚ »¢¼)ƒ¯}„j´nF…H€1Šº­ƒ¯}~yƒ¯°
1€|}g{¸}„‚fƒ­”}rƒ­lnF…H€1’º¯ƒ¯}~yƒ­´{„ˆ’mon”F’{„nF{)¼)jlnFln<®In<‚ p ≥ 3 ¿
3ÃSI“l“|ƒ¯}„ƒ­ˆ’&q » ≤q =p p ƒ¯°}fjln{„” j1n<monƒ­{)p n<≤Å|‘l2º¯ƒ†”ƒ¯}4ˆŠ°+ˆ’‚ “|n<‚)Š}4mdˆH{~}) ¢¿
41ˆ’‚ˆ’‚ “|n<‚t£’‚fnFŠ}„n<‚r}„j1’´°²ˆI€l‚K»XˆIlno”<Š1lˆŠ}“|ˆ”<ˆ’mo‘l€|} #}fƒ¯ˆI1{rzySj1’1“ï¿ 3Á}ƒ­{t}fjln<Ā1{„n°²€lº+}„ˆ
‚fn<º­yĈ’Æ}„j1n'}„j1n<ˆ’‚fyĈа4ˆ’‚ “|nF‚”ˆI1“|ƒ¯}„ƒ­ˆ’1{°²ˆ’‚>‘1Š‚„}„ƒ¯}„ƒ­ˆ’lnK“=)€ll£In·Ã¾€|}¸} O{f” jln<mon’¿›ijlƒ­{}„j1n<ˆ’‚fy’»
¢’{„nF“OˆI´lƒ†”ˆIº¯ˆI€l‚g‚„ˆzˆŠ}fnF“}f‚„nFnF{g¼)ƒ¯}„jļ)jlƒ­” jЊ‚fn>I{„{„ˆ|”ƒ†#}„nK“O”nF‚¸} Šƒ­´z€lm>¢nF‚f{F»Ûƒ­{Š›nÅz}„nF1{„ƒ¯ˆI
ˆ’°ï}„jlngˆ’‚fƒ­£’ƒ­1Šº¢¼ˆI‚„Àdz2y 5;€|} ” jln<‚°²ˆ’‚µ²1ˆ’
‘¢Š‚„}„ƒ¯}„ƒ­ˆ’lnK“1¹C)€1l£’n<·Á¾€|}„}f>{„” jlnFmonF{F»H{„n<6n 5€l}f” jlnF‚ " 1»
‘?¿Cš’š $Á¿
3Á}CŠ‘1‘¢nKŠ‚ {?}fj1#}}fjln)”º†’{f{?ˆ’°¢‘1’‚¸}fƒ±}fƒ¯ˆIlnF“>)€ll£In·Ã¾€|}¸} t{f” jln<monK{+”<ˆ’moƒ­l£r°²‚fˆ’m }fjln)“|ƒ†{„”<‚„n<}„ƒ­³F#·
}fƒ¯ˆI›ˆ’°ˆI‘|}„ƒ­mQŠºC”ˆIH}„‚fˆ’º‘l‚fˆ’lº­n<mQ{Š‚fn>ƒ­›°¶’”9}‘¢Š‚„}„ƒ¯}„ƒ­ˆ’lnK“´{„yHmo‘lº­nF”}„ƒ†”o)€ll£In·Ã¾€|}¸} {f” jlnFmdnK{<»
” j¢Š‚ ’”9}fn<‚fƒ¯³FnF“ozyQ‚fn<º†#}„ƒ­ˆ’ĵ² z¹CXn<º­ˆ#¼¿Cwzˆd}„jln…H€lnF{¸}„ƒ­ˆ’¢ˆIƒ¯º†{)“|ˆ#¼)'}fˆ}fjlnˆ’lngˆŠ°nÅ|‘l‚fnF{f{¸ƒ­l£ˆI‚¸·
“ln<‚”<ˆ’1“lƒ±}fƒ¯ˆI1{C°²ˆI‚*}fjlƒ†{”<º­I{„{F¿+ij1nrmQŠƒ­'‚fnF{„€lº±};ˆŠ°ï}fjlƒ­{‘1’‘¢nF‚ƒ†{C}fj1#};¼ng”<’Qˆ’|} Šƒ­Q}„jln“|nK{¸ƒ­‚„nK“
n<Å|‘l‚„nK{„{„ƒ­ˆ’1{;€1{„ƒ¯1£' )¸”FŠº†”€lº­€17{ *gˆ’ˆI‚„ƒ­n<H}„nK“'°²‚fn<n}f‚„nFnF{F¿*iˆQ¢n{„‘¢nK”ƒ ¢”Š»¢¼)ƒ±}fjSnFI” jˆI‚„ƒ­n<H}„nK“'°²‚fn<n
}f‚„nFnt’‚„ng’{f{¸ˆ|”ƒ†#}fnF“Q{„ˆ’mong¼nFƒ¯£IjI} {<»zŠ¢“o}„jlngmQŠƒ­
ˆI‘¢nF‚fŠ}„ƒ­ˆ’'ƒ†{C}f8ˆ )„{„‘lº¯ƒ¯7} *t’zyd‚fˆzˆŠ}fnF“o}f‚„nFntƒ­H}„ˆ>
{„€lmvµ¶¼)ƒ±}fjO”ˆzn 1'”<ƒ¯nFI} { ±1¹ˆŠ°{¸€1” jˆI‚„ƒ­n<H}„nK“Q°²‚fn<n}f‚„nFnF{F¿
ijln‘1’‘¢nF‚>ƒ†{dˆ’‚f£IŠ1ƒ¯³FnF“=’{°²ˆ’º­º­ˆ#¼4{<9¿ 3Ã=}„j1nln<Åz}o{„nF”}„ƒ­ˆ’½¼;n“ln}f’ƒ¯º}„j1n“|ƒ†{f”‚fn}„ƒ­³FŠ}„ƒ­ˆ’½ˆ’°
ˆI‘|}„ƒ­mQŠº&”<ˆ’H}„‚fˆ’ºï‘l‚fˆ’1º¯nFmo{zy)€ll£In·Ã¾€|}¸} d{f” jlnFmdnK{<»l’1“{¸j1ˆ#¼ }„j1n‚„nFº­Š}„ƒ­ˆ’¼)ƒ±}fj‘¢Š‚„}„ƒ¯}„ƒ­ˆ’lnK“
{„yzmo‘lº¯nK”9}fƒ­”)€l1£’n·Ã¾€|}„}fd{„” jlnFmonF{F¿{f#}fƒ­{¸°²yzƒ¯1£Oµ² H¹¿Cijln<ƒ­{„nF”}„ƒ­ˆ’˜>¼;n‚„nF®zƒ¯nF¼ }„jln}fjln<ˆI‚„yQˆ’°
ˆI‚f“ln<‚”ˆI1“|ƒ¯}„ƒ­ˆ’1{°²ˆ’‚‘1’‚¸}fƒ±}fƒ¯ˆIlnF“ď)€l1£’n·Ã¾€|}„}fS{f” jln<monF{F¿w|nF”9}fƒ¯ˆIÆ Sƒ¯H}„‚fˆ|“|€1”<nF{ˆI‚„ƒ­n<H}„nK“O°²‚fn<n
}f‚„nFnF{F»X’1“S{„jlˆ#¼4{4jlˆ#¼ }fjlnˆI‚f“|nF‚t”<ˆ’1“|ƒ¯}„ƒ­ˆ’¢{r”<’S¢ndnÅ|‘l‚fnF{f{„nF“ƒ¯O}„n<‚fmQ{rˆ’°Õ}„jln>º­Š}¸}fn<‚K¿ 4+ƒ¯1’º¯º­y
{„nF”}„ƒ­ˆ’—“|ƒ†{f”€1{f{¸nK{C}fjlnz€lmon<‚fƒ­”FŠºÛƒ¯mo‘lº­n<monFI} #}fƒ¯ˆI&»|Š1““|ƒ­{„‘lº†Ðy|{*}fjln‚„nK{¸€1º±} {C°²ˆI‚;ˆI‚f“ln<‚;€l‘
}fˆ Ô
’1“
}„jlnz€lm>¢nF‚)ˆŠ°Õ”ˆ’¢“|ƒ±}fƒ¯ˆI1{;°²ˆ’‚)ˆI‚f“ln<‚€l‘}„ˆ
Òz¿
!
i
õõó: !;
ïÉ Î É9Ê
ʶË<Ì
!*;!; ;!&?) <5(^;>?+( +"
/
(45!
n} ’1“ ¢n °²€l1”}„ƒ­ˆ’1{
Š¢“
»X‚fnF{„‘¢nK”9}„ƒ­®’nFº¯yI¿tŽˆ’1{„ƒ†“|n<‚r}„jln
°²ˆIº¯º­ˆ#f¼)ƒ­l£>€1Φ1”ˆI1{~C}f‚f’ƒ¯1nF“ˆ’‘|}fƒ¯mQŠIRº&”<ˆ’H×}„‚fIˆ’RºÛ‘l‚f→ˆ’1Iº¯RnFm# cR → IR
m
∞
n
n
n

 Min Φ(y(T ));
ẏ(t) = f (u(t), y(t)),

y(0) = y 0 .
t ∈ [0, T ];
(P )
¬ n‚fnF{¸}„‚fƒ­”})}„jln’1Šº­y|{¸ƒ†{;}„ˆ'”<ˆ’H}„ƒ­z€lˆ’€1{r”ˆ’H}f‚„ˆIº¢°²€11”9}fƒ¯ˆI1{<¿Œtn<lˆ’}„nzy y, p) := p · f (u, y)
}fjln8
Ë1Î’È %$ >Ê ÍÃÈ#Ì¢Ê#Ì=ˆ’°+}fjln‘l‚fˆ’1º¯nFm¿4ij1n1‚ {~}rˆI‚f“ln<‚4lnK”nK{„{fŠ‚fy'ˆ’H(u,
‘l}„ƒ­mo’º¯ƒ¯}~y”ˆ’¢“|ƒ±}fƒ¯ˆI1{4ˆ’°
}fjlƒ†{)‘l‚fˆ’lº­n<m8mQÐyQXn¼)‚„ƒ¯}¸}fn<ƒ­}fjln°²ˆ’º­º­ˆ#¼)ƒ¯l£>°²ˆ’‚f&m #






ẏ(t) = f (u(t), y(t)),
ṗ(t) = Hy (u(t), y(t), p(t)),
0 = Hu (u(t), y(t), p(t)),



p(T ) = Φ0 (y(T )), y(0) = y 0 .
t ∈ [0, T ],
(OC)

¬ n{„Ðyd}fj1#} (ū, ȳ, p̄) ƒ†{)Š
n<Åz}„‚fn<mQŠºÛƒ¯°ƒ¯}){f#}fƒ­{ 1nK{ (OC) µ ū ¢nFƒ¯1£od”ˆIH}„ƒ­H€1ˆ’€1{*°²€l1”9}fƒ¯ˆI¢¹9¿ n}
¢n’n<ÅH}f‚„nFmQŠº[¿ 3Á°
(ū, ȳ, p̄)
ƒ†{)ƒ­z®’n<‚„}„ƒ­lº­ng’º¯ˆIl£d}„jln}f‚fŠ‰~nF”}„ˆ’‚fy ,
µ¸¡K¹
u 7→ H (u, y, p)
}fjln<´Hy}fjln>ƒ­mo‘lº­ƒ­”<ƒ±}g°²€l1”}„ƒ­ˆ’1{r}„jlnFˆ’‚fn<m»Ûƒ¯›
{„mo’º¯º L 1n<ƒ­£’jz¢ˆI‚„j1ˆHˆ|“ˆ’°Õ}„jlƒ†{4}f‚fŠ‰~nF”}„ˆ’‚fy’»¢¼n
j¢Ð®’nd}„j1Š}
(ƒ & u = φ(y(t), p(t)) »Õ¼)jln<‚fn φ ƒ†{ C mo’‘l‘lƒ­l£1¿Œrn 1ln
=0
}fjlnÍMÉ 1Ë $Hd(u(t),
Ê ÍÃÈ#Ì1Ê Šy(t),
Ìǐ’{ p(t))
¿dqr{„ƒ¯l£ H (φ(y(t), p(t)), y(t), p(t)) = 0 »
H(y, p) := H(φ(y, p), y, p)
ˆI|}f’ƒ¯
µMžI¹
H (y, p) = H (φ(y, p), y, p); H (y, p) = H (φ(y, p), y, p).
ŽˆI1{„nF…H€ln<H}fº¯yI»1€l1“ln<‚gjHyz‘XˆŠ}fjlnF{„ƒ­{>µ~¡Ð¹9» (OC) ƒ­{tº¯ˆ|”FŠº­º¯ynF…H€lƒ­®#Šº­n<H}4}„ˆ
}„jlnÉ„Ë Î Ë 'Î $dÊ ÍÃÈ#Ì¢Ê #Ì
,9ÍÃË ẏ(t) = H (φ(y, p), y, p), t ∈ [0, T ],
µ¶˜I¹
−ṗ(t) = H (φ(y, p), y, p), t ∈ [0, T ],
¬ nlˆ#¼ }„€l‚f
}fˆo}„jln“lƒ­{fp(T
”€1{f{„)ƒ¯ˆI=ˆ’°Φ}fjln(y(T
“|ƒ†{„”<)),‚„n<}„ƒ­³Fy(0)
#}fƒ¯ˆI=ˆ’y°?}f.jlnˆ’‘|}fƒ¯mQ’º?”ˆII}f‚„ˆIºÛ‘l‚„ˆIlº­n<m (P ) ¿
ij1n4€ll£’n<·Á¾€l}¸}fo“|ƒ†{f”‚fn}„ƒ­³FŠ}„ƒ­ˆ’”ˆI1{¸ƒ†“|nF‚„nK“'ƒ­9 " %Ò $&ƒ†{
uu
∞
∞
u
u
y
y
p
p
p
y
0
0

Min Φ(yN );


P
 y
= yk + hk Psi=1 bi f (uki , yki ),
k+1
s
yki = yk + hk j=1 aij f (ukj , ykj ),



y0 = y 0 ,
(DP1 )
°²ˆI‚Šº­º k = 0 }fˆ N − 1 Š¢“ i = 1 }fˆ s ¿'Örn<‚fn h > 0 ƒ†{g}„jln'{„ƒ¯}„nQˆŠ°;}„j1n k}„jÄ}„ƒ­mon'{~}fn<‘&»+Š1“
ƒ­{g}fjlnQ{¸n<}ˆŠ°;4€ll£’n<·Á¾€l}¸}f”ˆzn1Q”<ƒ¯nFH}f{F¿ n}€1{€l1“ln<‚fº¯ƒ­lnd}„jlnQ” jlˆIƒ­”<n>ˆ’°€1{„ƒ­l£S“|ƒ(&Ûn<‚fn<H}
(a, b)
k
²î+õ Ñö
ɄÎIËÉ È#ÌÛÎ#ÊÑÍ[ʶÈ#Ì*ÏÈ#É
¯Ë<ÍMÊ
ŠÉ9ÍMÊÑÍ[ʶÈ#ÌÛË Î|ÌHË|Í[ÍlËQË
—
®#’º¯€lnK{;ˆŠ°+”<ˆ’H}„‚fˆ’º†{ ’{f{¸ˆ|”<ƒ­Š}„nF“'¼)ƒ±}fjƒ­llnF‚){~} #}„nK{ }fjlƒ­{)”<ˆ’H}„‚ ’{¸}f{¼)ƒ¯}„jˆŠ}„j1n<‚)Š‘l‘1‚„ˆH’” jlnK{
ƒ­
¼)j1ƒ­” j
}„jln“|ƒ†{„”<‚„n<u}„ƒ­³F#}fƒ¯ˆI'ˆŠ°+”ˆ’H}f‚„ˆIº­{ƒ­{;”ˆHŠ‚ {¸nF‚C}fj1yŠ
}„jlngˆ’lngˆŠ°?}„jln{~} #}fn’»zƒM¿ n’¿­»z°²ˆ’‚”ˆIH}„‚fˆ’ºXƒ­{
} ŠÀIn<”ˆ’¢{~} ŠH}F»|ˆI‚;ˆIlº­y
 1QlnˆInK’” j}„ƒ­mdn{¸}„nF‘ĵ¶n’¿ £1¿ "­¡’»1žl»1,˜ $²¹9¿
¬ nmQÐy'‚„nF¼)‚„ƒ¯}„n (DP ) €l1“ln<‚}„j1nnF…H€lƒ­®Ð’º¯nFH};°²ˆI‚„m
kj
kj
1

Min Φ(yN );



s

X

 0=h
bi Kki + yk − yk+1 ,
k
i=1

Ps


0 = f (uki , yk + hk j=1 aij Kkj ) − Kki ,



0 = y 0 − y0 .
P
yk + hk sj=1 aij Kkj
yki
3Ã'}„jlntnÅ|‘l‚fnF{f{¸ƒ­ˆ’QXn<º­ˆ#¼¼ng”ˆ’H}f‚fI”9}
I{„{„ˆ|”ƒ†#}fnF“'¼)ƒ¯}„j
Φ(yN )+
N
−1
X
(
ƒ­H}„ˆ ¿Õijln ’£’‚ Šl£Iƒ­’°²€l¢”9}„ƒ­ˆ’
Ġ{ #
(DP )
2
pk+1 ·
hk
s
X
(DP2 )
bi Kki + yk − yk+1
!
+
s
X
)
ξki · (f (uki , yki ) − Kki ) +p0 ·(y 0 −y0 ).
rÖ n<‚fn » »’’1“ Š‚fn*}fjln ’£’‚ Šl£In*m€1º±}fƒ¯‘lº­ƒ­n<‚ {Ր’{f{¸ˆ|”ƒ†#}fnF“¼)ƒ¯}„jd}fjln)”ˆ’¢{~}f‚f’ƒ¯H}f{+ˆ’° ) ¿
C’‚„ƒ†Šp1º¯nK{ pξ ¼)ƒ­º­ºXnpƒ­H}„n<‚f‘l‚fn}fnF“S’{4}„jlnd“|ƒ†{„”<‚„n<}„ƒ­³F#}fƒ¯ˆISˆŠ°C”ˆ’·Á{¸}fŠ}„nˆ’°C”ˆ’H}fƒ¯z€lˆI€1{)°²ˆI‚„m>€l(DP
º­Š}„ƒ­ˆ’&¿
w|n}¸}fƒ¯1£}„ˆ>³<nF‚„ˆ}fjln“|n<‚fƒ¯®#Š}„ƒ­®’n4ˆ’°ï}„jlƒ†{ Š£I‚f’l£’ƒ†Šd°²€l1”}„ƒ­ˆ’¼¿ ‚K¿ }K¿?}fjlng‘l‚„ƒ­mQŠºX®#Š‚fƒ­’lº¯nK{ y » y »
°²ˆI‚ k = 1 }„ˆ N − 1 » K ’1“ u °²ˆ’‚ k = 0 . . . N − 1, i = 1 . . . s »|¼;nˆ’l}fŠƒ­
y
k=0
k+1
i=1
0
ki
i=1
2
k
0
N
k
ki
ki
= Φ0 (yN ),
0
= p
P,s
>
pk − pk+1 =
i=1 fy (uki , yki ) ξki ,
s
X
0 = hk bi pk+1 + hk
aji fy (ukj , ykj )> ξkj − ξki ,
pN
p1
j=1
tq {¸ƒ­l£
1ˆ#¼ }fjlnjzyz‘XˆŠ}„j1nF{„ƒ­{4}„j1Š} b 6= 0 »Û{¸n<} p := ξ /(h b ) °²ˆI‚tiŠ=º­º 1k .=. . s.0 }„ˆ N − 1 »X’1“
}„ˆ ¿2.Cº­ƒ¯moƒ­1#}fƒ¯1£ ξ °²‚„ˆIm }„jln'’¢ˆ#®In>‚fn<º†#}fƒ¯ˆI1{F»ï¼noˆ’l}fŠƒ­›}„jlnQnK…I€1ƒ¯®#Šº­n<H}ˆ’‘l}„ƒ­mo’º¯ƒ¯}~y
i=1
”<ˆ’1“lƒ±}fƒ¯ˆI1s{
0 = fu (uki , ykj )> ξki , k = 0 . . . N − 1,
i
ki
ki
k i
ki
Ps

 yk+1 = yk + hk Pi=1 bi f (uki , yki ),
s
p
= pk − hk i=1 b̂i Hy (uki , yki , pki ),
 k+1
0 = Hu (uki , yki , pki ),
yki
pki
y0
Ps
= yk + hk j=1 aij f (ukj , ykj ),
Ps
= pk − hk j=1 âij Hy (ukj , ykj , pkj ),
= y0,
pN = Φ0 (yN ),
(DOC)
¼)j1n<‚fn”<ˆzn1'”ƒ­n<H}f{ b̂ Š¢“ ⠐’‚„n“|n 1lnK“zyQ}fjln°²ˆ’º­º­ˆ#¼)ƒ¯l£d‚fn<º†#}„ƒ­ˆ’¢{ #
°²ˆ’‚4’º¯º i = 1, . . . , s Š¢“ j = 1, . . . , s.
µ² z¹
b
b̂ := b , â := b − a ,
b
3Á°|}fjln*Šº­£’nFl‚ Šƒ†”Õ”ˆ’¢{~}f‚f’ƒ¯H}f{
’‚„nÕº­ˆ|”<’º¯º­yrnK…I€1ƒ¯®#Šº­n<H}&}fˆ
»
}fjln< (DOC) ƒ†{gnF…H€lƒ­®Ð’º¯nFH}t}fˆ'}„Hj1no(u{fŠmo,n>y ‘1Š, p‚„}„ƒ¯}„)ƒ­ˆ’=lnK0“´)€1l£’n<·Á¾€|}„}f{„” jlnFmond’u‘l‘lº­ƒ¯=nK“φ(y}„ˆ}„jl, nop ‚fn)·
“l€1”nK“'{¸y|{¸}„nFmÓµM˜I¹¿ij1ƒ­{;Š‘l‘1‚„ˆH’” jo1’{„nF“oˆIQ°²ˆ’‚fm€1º­Š}„ƒ­ˆ’ (DP ) ƒ†{{„º­ƒ¯£IjI}fº¯yQ{„ƒ¯mo‘lº­n<‚K»Hl€|}nF…H€lƒ­®I·
’º¯nFH};}fˆd}„jlnˆIlnˆŠ°ÖrŠ£In<8‚ " %Ò $Á¿
j
i
i
ij
j
ji
i
u
ki
ki
ki
ki
2
õõó: !;
ki
ki
Ô
ïÉ Î É9Ê
ʶË<Ì
3Á}ƒ­{{fŠƒ†“}„j¢#}Րr‘1’‚¸}fƒ±}fƒ¯ˆIlnF“)€ll£In·Ã¾€|}¸} r{f” jln<montµ¶ˆ’‚+moˆ’‚fn*£’nFln<‚ Šº­º­yg’HyˆIln{¸}„nF‘>{f” jln<monK¹
ƒ†{g{¸yzmo‘lº­nF”9}fƒ­”ˆŠ°C}fjlno”ˆ’‚f‚fnF{„‘¢ˆI1“|ƒ­l£1ˆ#¼ ƒ†{t{„yzmd‘1º¯nK”9}„ƒ†”Š¿63Á}ƒ­{tÀzlˆ#¼)S}„j¢#}‘1Š‚„}„ƒ¯}„ƒ­ˆ’lnK“S)€ll£In·
¾€l}¸}fÄ{f” jln<monF{Q{f#}fƒ­{¸°²yzƒ¯1£µ¶ H¹dŠ‚fn{¸yzmo‘lº­nF”9}fƒ­”’»{„n<n9"¯¡KŸ|»ijln<ˆI‚„nFm 1¿ Ô $Á¿ ¬ nˆI|}f’ƒ¯Ç}„j1Š}Q}„jln
{f” jlnFmdnˆI|}f’ƒ¯lnK“
zy“|ƒ­{f”‚fn}fƒ¯³K#}fƒ¯ˆIˆ’°‘l‚fˆ’lº­n<m ƒ†{){¸yzmo‘lº­nF”9}fƒ­”’¿ 3Ñ1Š‚„}„ƒ†”€1º­’‚}fjln°²ˆ’º­º­ˆ#¼)ƒ¯l£
“lƒ­’£’‚ Šm8”ˆImom€|}fnF{F»|¼)jln<¼;ng€¢{¸ng}„j1nŠXˆ#®’n“|ƒ­{f(P”‚f)n}fƒ¯³K#}fƒ¯ˆI #
“lƒ­{f”‚fn}fƒ¯³K#}„ƒ­ˆ’
(P ) −−−−−−−−−−→ (DP )
ˆ’‘|}fƒ¯mQ’º¯ƒ¯}~y 
ˆ’‘|}fƒ¯mQ’º¯ƒ¯}~y 
”ˆ’¢“|ƒ±}fƒ¯ˆI1{ y
”ˆ’¢“|ƒ±}fƒ¯ˆI1{ y
(D)
“lƒ­{f”‚fn}fƒ¯³K#}„ƒ­ˆ’
(OC) −−−−−−−−−−→ (DOC)
41ˆ’‚“ln}f’ƒ¯º­nF“o‘1‚„nK{¸nFI} #}fƒ¯ˆI>ˆ’°Û‘¢Š‚„}„ƒ¯}„ƒ­ˆ’lnK“Q)€ll£In·Ã¾€|}¸} Š1“o{„yHmo‘lº­nF”}„ƒ†”)mon}fjlˆ|“l{C¼;n4‚fn°²n<‚
}fˆd}„jlnXˆzˆ’À|8{ " š1»l™ $Á¿
+ ; +! # ?9) ,./102'?Ó!?345)!
¢’{„ƒ­”}fˆHˆIº4ƒ­}„jln´{¸}„€¢“|yLj’°ˆI‚f“ln<‚'”ˆI1“|ƒ¯}„ƒ­ˆ’1{d°²ˆI‚
)€l1£’n·Ã¾€|}„}f={f” jln<monF{oƒ†{Q}fjln}„j1n<ˆ’‚fyLj’°
5·Ã{¸nF‚„ƒ­nF{gŠ¢“´’{f{„ˆz”<ƒ­Š}„nK“O”<’º­”<€lº¯€¢{tˆI´‚fˆzˆŠ}„nK“S}f‚„nFnF{F»&“|€ln>}„ˆ 5;€|}f” jlnF‚Qµ¶{„n<n " $ѹ
¿ 4lˆI‚t‘¢Š‚„}„ƒ¯}„ƒ­ˆ’lnK“
4€ll£’n<·Á¾€l}¸}fg{f” jln<monF{Ր’>n<Åz}„n<¢{¸ƒ­ˆ’dˆŠ°X}„jlƒ†{}fjln<ˆI‚„yƒ†{}„j1n)ˆ’lnˆŠ°XnÅ|‘1’1{¸ƒ­ˆ’dƒ­d‡C·Á{„n<‚fƒ¯nK{<»Š’1“>}„jln
I{„{„ˆ|”ƒ†#}fnF“”<’º­”<€lº¯€¢{4ˆ’´lƒ­”<ˆ’º­ˆ’€l‚fnF“‚fˆzˆŠ}„nK“}„‚fn<nK{<»Û{¸nF/n " š $Á¿tijlnº†#}„}„nF‚g’º¯º­ˆ#¼4{}„ˆ{~} #}fn}„jln>ˆ’‚ “|n<‚
”<ˆ’1“lƒ±}fƒ¯ˆI1{)ƒ­S}„nF‚„mQ{4ˆ’°C”ˆzn 1'”ƒ­n<H} { a » b » ⠻ېŠ1“ b̂ »¢ˆ’°}„jln°²ˆIº¯º­ˆ#¼)ƒ­l£o}~yz‘¢n #;”<n<‚„}fŠƒ­S‘Xˆ’º­yH1ˆ’moƒ­’º­{
µ¶¼)jlƒ†” jŠ‚fn*ƒ­°¶’”}{„€lm ˆŠ°¢mdˆIlˆ’moƒ†Šº†{&¼)ƒ¯}„j>€llƒ¯}”<ˆHn 1'”ƒ­n<H}f{ ¹ïƒ­}fjlnF{„n*®#Š‚fƒ†Šlº­nF{?j1Ð®’nÕ}fˆtXn*nK…H€1Šº
}fˆ4”nF‚¸} Šƒ­t°²‚ ’”}„ƒ­ˆ’1’º#z€lmXn<‚ {<¿+ijlnC{„€l1{¸}„ƒ¯}„€|}fƒ¯ˆIˆŠ° â»#Š¢“ b̂ µ¶€1{„ƒ¯l£dµ² z¹¸¹Û¼;ˆ’€lº†“g£’ƒ­®’nՔˆImo‘lº¯ƒ†”<Š}„nK“
n<Å|‘l‚„nK{„{„ƒ­ˆ’1{F»|{¸ƒ­1”nƒ­1{¸}„nK’“ˆŠ°+nF’” jmoˆ’1ˆ’moƒ­’ºï¼n¼;ˆ’€lº†“j1Ð®’nQ{„€lm"ˆ’°+‚ #}fƒ¯ˆI1ŠºÛ°²‚ ’”}„ƒ­ˆ’1{F¿ ¬ n
¼)ƒ­º­ºC{¸jlˆ#¼•}„j1Š}’moˆ’l£Šº­º+}„jlnK{¸nd}„nF‚„mQ{ƒ¯}ƒ­{{„€ 1'”ƒ­n<H}g}„ˆnÅ|‘l‚fnF{f{”<ˆ’1“|ƒ¯}„ƒ­ˆ’Ĉ’Đ )~‘l‚fƒ­1”ƒ­‘1Šº
}fn<‚fm *4{„ƒ­1”nr}„jlntˆŠ}„j1n<‚C}fn<‚fmQ{*ˆ’°ï}„jlng{¸€lm7’‚„nr’º¯‚fnFI“|y>“ln}„nF‚„moƒ­lnF“'Hyo‘l‚fn<®zƒ­ˆ’€1{”ˆI1“|ƒ¯}„ƒ­ˆ’1{F¿+ij1ƒ­{
’º¯º­ˆ#¼4{r'}„‚fn<mon<¢“|ˆ’€1{g{„ƒ¯mo‘lº­ƒ X”<#}fƒ¯ˆIS}„j¢#}‘¢nF‚„moƒ¯}f{t€1{t}„ˆ“lƒ­{„‘lº†Ðy}fjlnF{„n>”<ˆ’1“lƒ±}fƒ¯ˆI1{r€1‘O}„ˆˆ’‚ “|n<‚
Ô »|ƒ¯}fn<‚fmo{ˆ’° a Š1“ b ’º¯ˆIln’¿
3È’‚ “|nF‚}fˆ lÅlˆ’}fŠ}„ƒ­ˆ’¼n¼)ƒ­º¯º&1n<nF“{„ˆ’mon“|n 11ƒ±}fƒ¯ˆI1{<¿ ’É G ƒ†{)o‘1Šƒ­‚ (V, E) ¼)jlnF‚„n
1
j
’“› 11ƒ±}fnd”FŠ‚ “|ƒ¯¢Šº’1“
8¿ .Cº­n<mon<H} {tˆ’° µ¶‚„nK{¸‘?¿ ¹tŠ‚fn>”FŠº­º¯nK“S®’nF‚¸}fƒ­”<nF{>µ²‚fnF{„‘&¿
V
nK“|£’nK{f¹¿ÕÖ4n<‚fn P(V ) “|nFlˆŠ}fnF{;ƒ†{;E}„jl⊂nP(V
{„n})ˆŠ)°{¸€l¢{¸n<}ˆ’° V ¿CiVjln{¸n<}ˆ’° E’É ¼)ƒ­º¯ºïXn“|n<1ˆŠ}„n Ḡ ¿ 
n<º­n<mon<H}f{rˆŠ° E
ȊÉ9ʶËÌ¢ÍÁË 2Î ’É ƒ†{rQ‘¢Šƒ­‚
¼)jln<‚fn ƒ†{rI{4Xn°²ˆ’‚fn’1“SŠ1“
’‚„n”<’º¯º­nF“ˆI‚„ƒ­n<H}„GnK“'nK“|£’nK{<¿Õi(V,jlnE){„n})ˆŠ°4È#É9ʶË<VÌ1ÍÃË Î ’É S¼)ƒ­º­ºï¢n“|nFlˆŠE}fn ⊂G ¿ V × V


© Ø V = {a, b, c, d} Š1“ E = {(a, b), (c, a), (d, b), (a, d)}
”
“
‘¢#}„jƈ’°)£I‚f’‘ljă­{
 ¢lƒ±}fn
{„nF…H€ln<¢”nQˆŠ°®’nF‚¸}fƒ­”<nF{QµM#}º­nFI{~}}~¼;ˆH¹{„€1” jÄ}„j1Š}}„jlnF‚„nQn<Å|ƒ­{¸}f{
’ÆnF“|£In'¢n<}~¼nFn<Æ}~¼;ˆO{¸€¢”<”nK{„{„ƒ­®’no®’nF‚¸}fƒ­”<nF{F¿ £I‚f’‘ljƃ­{ È#Ì1ÌÛË ÍÃË Î´ƒ¯°rŠzy›‘1Šƒ­‚ˆŠ°)®In<‚„}„ƒ†”nF{Š‚fn
²î+õ Ñö
ɄÎIËÉ È#ÌÛÎ#ÊÑÍ[ʶÈ#Ì*ÏÈ#É
¯Ë<ÍMÊ
ŠÉ9ÍMÊÑÍ[ʶÈ#ÌÛË Î|ÌHË|Í[ÍlËQË
Ò
monFmXn<‚rˆ’°C#}rº­nFI{~}4ˆIln‘1Š}„j&¿ <ÊÑÉ |ÊÑÍ)ƒ­{tQ‘1Š}„jS¼)jlˆH{¸n 1‚ {¸}tŠ1“º†’{¸}rnFº¯nFmdnFH}f{r”<ˆ’ƒ­1”ƒ†“|nI¿ Í[É„Ë Ë G ƒ†{”<ˆ’l1nF”9}fnF“´£’‚ Š‘ljO¼)ƒ¯}„j1ˆ’€|}”ƒ­‚ ”€lƒ¯}f{F¿ijlno{¸n<}ˆŠ°C}f‚„nFnF{g¼)ƒ¯º­º¢nQ“|nFlˆŠ}fn T̄ ¿ ɄÈKÈ#ÍÃËfÎ
Í[É„Ë Ë ƒ­{‘1Šƒ­‚ˆŠ°
}f‚„nFnd’1“›ˆŠ°ˆIlnoˆŠ°*ƒ¯}f{®’n<‚„}„ƒ†”nK{<»&”FŠº­º¯nK“S}fjlno‚„ˆzˆŠ}K»?Š¢“›“|nFlˆŠ}fnF“ ¿'ijln
º­nFÐ®InFt{;Š‚fnr}fjln®’nF‚¸}fƒ­”<nF{;¼)ƒ±}fjˆIlº­y'ˆ’ln’“#‰¸’”nFH}nK“|£’nI»zˆŠ}„j1n<‚;}„j1’
}fjln‚fˆHˆ’}F¿ ¬ ngmQÐyQr(t)
ƒ†“|n<H}fƒ±°²y

‚fˆzˆŠ}fnF“
}„‚fn<n¼)ƒ±}fjSŠˆ’‚fƒ¯nFI}fnF“£’‚ Š‘lj?»H¼)j1n<‚fn’º¯º?nF“|£InF{;‘¢ˆIƒ¯H}4ƒ­S“|ƒ¯‚fnF”}„ƒ­ˆ’ˆŠ°º¯nKÐ®’nK{<¿
ijlnd{¸n<}rˆŠ°)ɄÈKÈ#ÍÃË Î'Í[É„Ë Ë¼)ƒ¯º­ºXn“|nFlˆŠ}fn ¿ ¬ n{¸ˆImon}„ƒ­monF{t{¸‘XnF’Àˆ’°+°²‚fn<n}„‚fn<nF{F»¢I{4ˆI‘l‘¢ˆH{¸nK“
}fˆO‚fˆHˆ’}„nK“´}f‚„nFnF{F¿ ¬ n
“|n 1lnƒ¯½S{„ƒ¯moƒ­º­’‚TmQŠl1n<‚ˆI‚„ƒ­n<H}fnF“›}„‚fn<nK{<»Šº†{¸ˆO”FŠº­º¯nK“ÆÖ)·[}„‚fn<nF{ƒ¯ "­¡K,ž $[¿
ij1n{„n}4ˆ’°ÈŠÉ9ʶËÌ¢ÍÁË ÎtÏ É„Ë ËÍ[É„Ë Ëƒ†{4“|n<1ˆŠ}„nK“ T ¿
5;ƒ†”ˆ’º­ˆ’€1‚„nK“£I‚f’‘lj1{4’‚„n£I‚f’‘lj1{}fˆ’£In}„j1n<‚4¼)ƒ¯}„j´'mo’‘l‘lƒ­l£Q}fj1#}4}fˆ'nFI” jS®’nF‚¸}fnÅ ’{f{„ˆz”<ƒ­Š}„nK{
S”<ˆ’º­ˆ’€l‚ c(v) »ˆŠ°)®#Šº­€ln B ˆ’‚ W µ¶lº­I” À›’1“ļ)jlƒ±}fnK¹¿ij1n'{¸n<}ˆŠ° <Ê fÈ ¯È |ÉfËfÎÄÈ#É9ʶË<vÌ1ÍÃË Î ’É µ¶‚„nK{¸‘?¿ <Ê È ¯È |É„Ë Î'ÉfÈFȊÍÁË ÎQÍ[É„Ë Ë K¹¼)ƒ¯º­º&Xn“|n<lˆ’}„n BG µ²‚fnF{„‘&¿ BT ¹¿
¬ nlˆ’}„n
µ¶‚„nK{¸‘&¿
¹;}„j1n{¸n<}4ˆŠ°lº†’” À›µ¶‚„nK{¸‘?¿*¼)jlƒ¯}„nй®’nF‚¸}fƒ­”<nF{
’1“ E µ²‚fnF{„‘&V¿ E= c¹*}fjl({B})
n{¸n<}4ˆŠ°nF“|£InF{;Vn<¢“|=ƒ¯1c£Qˆ’({Wlº­I})” À´µ²‚fnF{„‘&¿¼)jlƒ¯}„nй®’nF‚¸}fƒ­”<n’¿
41ˆ’‚£’ƒ­®’nFQ‚„ˆzˆŠ}fnF“o}f‚„nFn’»H¼)jlˆI{„n4®In<‚„}„ƒ†”nF{Š‚fnr’{f{„ˆz”<ƒ­Š}„nK“o¼)ƒ±}fj'º¯n<}¸}fn<‚ { i, j, . . . ƒ±};ƒ­{;”ˆIz®’n<1ƒ¯nFI}
}fˆQ“|n<lˆ’}„n
ƒ¯°+®In<‚„}„nÅ k ƒ­{¼)jlƒ¯}„n , b̂ ˆŠ}„j1n<‚f¼)ƒ­{„n’¿
µM—’¹
b̃
= b
ƒ¯°®’n<‚„}„n<Å ` ƒ­{4¼)jlƒ±}fn , â ˆ’}„jlnF‚„¼)ƒ†{„n’¿
ã
= a
Örn<‚fn i I{„{„ˆ|”ƒ†#}„nK{C¼)ƒ¯}„jnFI” j'®’nF‚¸}fnÅ k ∈ {1, . . . , #t} ®#Š‚fƒ­’lº­n i µ¶®#Š‚fyHƒ­l£°²‚„ˆIm¡r}„ˆ s ƒ¯
}„jln
{„€l1{„nF…H€lnFH}n<Å|‘l‚„nK{„{„ƒ­ˆ’1{ ¹9¿
¬ ƒ¯}„jSnK’” j
‚fˆzˆŠ}fnF“
}„‚fn<n t ƒ†{4I{„{„ˆ|”ƒ†#}„nK“
d®#Šº­€ln γ(t) ƒ¯SŠƒ­1“|€1”}„ƒ­®’n¼Ðy # γ( ) = 1 »1Š1“
µÔ ¹
γ(t) = #t × γ(t ) × · · · × γ(t ),
¼)j1n<‚fn;zy
¼;n“ln<lˆ’}„n;}„jln)1‚f’1” jlnF{ˆ’° µÑ}fjln4”ˆIllnF”}„nK“>£I‚f’‘lj1{+ˆ’l}fŠƒ­lnK“zy‚fn<moˆ#®zƒ­l£
}fjln‚fˆzˆŠ})°²t‚fˆ’,m . . t. »1, t®zƒ¯nF¼nK“S’{)‚fˆzˆŠ}fnF“}„‚fn<nK{<»1}„jln>‚„ˆztˆŠ}4Xn<ƒ­l£Q}fjln®In<‚„}„nż)jlˆI{„nˆ’1n’“Љ¸I”nFI}4nK“|£’n
j¢’{Xn<n<‚fn<moˆ#®’nK“1¹9¿
prnÅz}F»Û¼;nlnFnF“S}fˆ
‚fnF”FŠº­º&}„j1n>“|n 1lƒ¯}„ƒ­ˆ’´ˆ’°Õ}„jln>n<º­n<mon<H} Š‚fy¼;n<ƒ­£’jH}tˆ’°*Q}f‚„nFn t »X¼)jln<‚fn t ƒ†{t
1ƒ­”<ˆ’º­ˆ’€l‚fnF“
‚„ˆzˆŠ}fnF“£’‚ Š‘lj?»I°²ˆI‚)£’ƒ­®’nF)€1l£’n<·Á¾€|}„}fo”ˆzn 1Q”<ƒ¯nFH}f{ a ’1“ b #
X
µÁҊ¹
b̃ φ (t).
φ(t) =
ij1n>ŠXˆ#®’n°²€l1”9}fƒ¯ˆI1{ φ (t) Š‚fn}fjln<mQ{„n<º¯°“|n ¢lnF“Sƒ­›’Oƒ­1“|€1”}„ƒ­®’n¼Ðy’»Xzy φ ( ) = φ ( ) = 1
’1“
µ¶šI¹
φ (t) = ψ
P(t ) × · · · × ψ (t ),
¬ nj1Ð®’ng}fjln°²ˆ’º­º¯ˆ#¼)ƒ­l£d‚„nK{¸ψ€1º±(t}µ¶{„)n<n =" š|»1ij1m 3 3 39ã¿ žl¿ — φ$Ñ!¹ #
Î |Ì HË |Í[Í lË oË , ¯È ȊɄΒË<É q Ê ØzÝ ¨ Ø ŠÉ9ÍMÊÑÍ[ʶÈ#ÌÛË ÏÈ#É t ∈ BT , #t ≤ q.
µM™I¹
1
φ(t) =
,
γ(t)
wz€11{~}fƒ±}f€|}„ƒ­l£r}„jlnÕn<Åz‘1‚„nK{„{„ƒ¯ˆI1{ۈа ⠐’1“ b̂ ƒ¯'µ¶ H¹Xƒ­}„jlnC°²ˆ’‚fm€lº†;ˆŠ° φ(t) ¼;ˆ’€lº†“g£’ƒ­®’nՔˆImd‘1º¯ƒ†”<Š}„nF“
n<Å|‘l‚„nK{„{„ƒ­ˆ’1{F¿ ¬ n{„jlˆ#¼ ƒ­}fjlnln<ÅH}r{„nF”}„ƒ­ˆ’jlˆ#¼ }fˆQ{~} #}„nnK…H€lƒ¯®#’º¯nFI}K»|l€|}r{„ƒ¯mo‘lº­n<‚4”<ˆ’1“lƒ±}fƒ¯ˆI1{<¿
∗
∗
−1
B
B
W
−1
W
ik
ik
ik
ik i`
ij
ij
k
k
1
1
m
m
s
i i
i=1
i
i
i
i
i
k
1
s
jk =1
i
ijk
jk
∗
õõó: !;
m
i
š
ïÉ Î É9Ê
ʶË<Ì
(;(<!$:)>)-Ð+) ?))!
3Ã>}„j1ƒ­{Õ{„nF”}„ƒ­ˆ’¼;n” jlˆzˆI{„n;I{+”<ˆ’z®’nFH}„ƒ­ˆ’}fj1#}}„j1n{„n}Ո’°¢®’nF‚¸}fƒ­”<nF{ V ƒ­{ˆ’°l}„jln;°²ˆI‚„m # {1, . . . , #V }
¼;n)“|ˆlˆ’}Cº¯ˆH{¸n£’nFln<‚ Šº­ƒ¯}~y’¿ij1n<o¼n4{¸}f’‚¸}՝zy>{¸}fŠ}„ƒ­l£Šº¯}„nF‚„¢#}„ƒ­®’n;nÅ|‘l‚fnF{f{„ƒ¯ˆI1{ˆŠ°X}„jln)nFº¯nFmon<H}f’‚„y
¼;n<ƒ­£’jH} φ(t) »C¼)jlnF‚„n t ƒ†{>´‚„ˆzˆŠ}fnF“=1ƒ­”<ˆ’º­ˆ’€l‚fnF“Ä}f‚„nFn’»C“|n 1lnF“ǃ¯ µ¶™H¹9»°²ˆ’‚oO£In<lnF‚f’º*‘1’‚¸}fƒ±}fƒ¯ˆIlnF“
{f” jlnFmdnI¿ i+ŠÀzƒ­l£›ƒ­H}„ˆ=’”F”ˆI€lH}>}fjlnS“|n 11ƒ±}fƒ¯ˆI ˆ’°r°²€l1”}„ƒ­ˆ’1{ ψ (t ) »ˆI|}f’ƒ¯Ç}„jlnSmoˆI‚„nnÅ|‘lº­ƒ­”<ƒ±}
n<Å|‘l‚„nK{„{„ƒ­ˆ’
X
X
Y
µ¸¡FŸI¹
φ(t) =
ã .
b̃
Örn<‚fn>’£IŠƒ­ i I{„{„ˆ|”ƒ†#}„nK{r¼)ƒ¯}„j›nFI” jO®’nF‚¸}fnÅ k ∈ {1, . . . , #t} µ²®In<‚„}„nŽ¡>ƒ­{t}fjln>‚fˆzˆŠ}9¹r
®Ð’‚„ƒ†Š1º¯n
µ²®#Š‚fyzƒ¯1£>°²‚fˆ’mv¡t}fˆ s¹¿ 4lˆI‚„m>€lº­µ¸¡FŸI¹mQÐy'Xn¼)‚„ƒ¯}¸}fn<S’{
i
XX Y
Y ã
µ¸¡’¡K¹
φ(t) =
.
b̃
i
s
k
s
ik i`
i1
i1 =1
i2 ,...,i#t =1 (k,`)∈E
k
k
s
ik i`
ik
3â“|n<nK“ﻐ’º¯º
b̃i`
ƒ¯ }„j1n’¢ˆ#®Inn<Å|‘l‚„nK{„{„ƒ­ˆ’ǔFŠ1”<n<ºnÅl”nF‘|} »{„ˆ´}fj1#}Oµ¸¡FŸI¹>Š¢“ µ¸¡’¡K¹>”ˆIƒ¯¢”ƒ†“|n’¿
t ¢{¸nF‚„®In’»Xjlˆ#b̃¼nF®’nF‚F»X}„j¢#}g}„jlnQŠXˆ#®’n°²ˆ’‚fm€lº†'mQ’À’nF{{„n<1{„n>b̃ °²ˆ’‚oµ²lˆI›lnK”nK{„{fŠ‚fƒ¯º­yS”ˆIllnF”}„nK“1¹41ƒ±·
”<ˆ’º­ˆ’€l‚fnF“>ˆ’‚fƒ­n<H}„nK“£I‚f’‘lj1{F¿ zyd{¸€1” jd£’‚ Š‘1j ƒ­{* 1lƒ¯}„n4€l1ƒ¯ˆIdˆ’°Û”<ˆ’l1nF”9}fnF“d£’‚ Š‘1j1{+¼)ƒ±}fjQn<mo‘|}~y
ƒ­H}„nF‚f{„nF”}„ƒ­ˆ’Oˆ’°*®’nF‚¸}fƒ­”<nF{F»Û”FŠº­º¯nK“ È#Ì¢ÌXË<ÍÁË Î tÈ
lȊÌXË<Ì1Í >ˆŠ° t »&Š1“›“|n<lˆ’}„nK“ {t , q ∈ Q} µ²1ˆŠ}g}„ˆ
Xn”ˆI|°²€1{„nF“¼)ƒ±}fjl‚f’1” jlnK{ˆ’°d‚fˆzˆŠ}„nK“'}f‚„nFn’»l“|nFlˆŠ}fnF“ t ¹¿
Ø ½© 1Ë'Ë ­Ë QËÌ¢Í ŠÉ *Ë(Ê ,zÍtÈ¸Ï Ê È ­È |É„Ë ÎOȊÉ9ʶËÌ¢ÍÁË Î ’É t = (V, E, c) ÕÊÑÍ È#Ì ÌÛË ÍÃË 9Î È 1È#ÌÛËÌ¢Í {t , q ∈ Q} Ê Í lË ÛɄÈFÎ ÍdÈ~ÏÍ 1ËOË ­Ë QËÌ¢Í ŠÉ *Ë<(Ê ,zÍ SȸÏÊÑÍ È#Ì¢ÌXË <ÍÁË Î
È 1È#ÌÛËÌ¢Í %Ê ÑË Y
µ¸¡Kž’¹
φ(t) =
φ(t ).
ɄÈFȸ!Ï ¬ nj1Ð®Inr}fj1#}
v∈V iv =1 k∈V
(k,`)∈E
i`
i1
q
i
q
q
q∈Q
φ(t)
=
s Y
XX
v∈V iv =1 k∈V
=
s Y
XX
v∈V iv =1 q∈Q
=
Y


X
ãik i`
Y
b̃ik
b̃i`
(k,`)∈E


s
X
Y
b̃ik
k∈Vq


Y
Y
(k,`)∈Eq
b̃ik

ãik i` 
b̃i`
Y

ãik i` 
,
b̃i`
¼)j1n<‚fn)}„jlntº­I{~};nF…H€1’º¯ƒ¯}~y>€¢{¸nK{C}fjlnrƒ†“|n<H}fƒ±}~ySµ²®#’º¯ƒ†“o°²ˆ’‚Š‚flƒ¯}„‚ Š‚fy°¶Šmoƒ­º¯ƒ­nF{ˆŠ°?{„n} { I Š¢“o°²€l1”}„ƒ­ˆ’1{
¹
A




Y X
XX Y
µ¸¡F˜I¹


A (i) =
A (i) .
q∈Q
v∈Vq iv =1
k∈Vq
(x,y)∈Eq
q
q
q
q∈Q
i∈Iq
q
r∈Q i∈Ir
q∈Q
²î+õ Ñö
ɄÎIËÉ È#ÌÛÎ#ÊÑÍ[ʶÈ#Ì*ÏÈ#É
¯Ë<ÍMÊ
ŠÉ9ÍMÊÑÍ[ʶÈ#ÌÛË Î|ÌHË|Í[ÍlËQË
i j1n”<ˆ’1”<º¯€¢{¸ƒ­ˆ’°²ˆ’º­º¯ˆ#¼4{F¿
ƒ­®’nFOŠˆ’‚fƒ¯nFH}„nF“£’‚ Š‘1j
}fjlˆI{„nˆŠ° F ƒ†{4“|nFlˆŠ}fnF“’{
t = (V, E)
»¢’1“
F ⊂E
™
»1}„jln>{¸n<}4ˆŠ°CŠ‚ ”<{ƒ­Sˆ’‘l‘XˆI{„ƒ¯}„n“|ƒ­‚„nK”9}fƒ¯ˆI}fˆ
E > := {(x, y) ∈ V × V ; (y, x) ∈ E}.
ØzÝ ¨ Ø l Ëd˯ËQË<Ì1Í#É
1 ȯΠ#Í[Ê ;Ë
*ËÊ(,zÍ)È~Ï Ê È­È|É„Ë Î'ȊÉ9ʶËÌ¢ÍÁË Î’É φ(t) =
X
t = (V, E, c)
>
(−1)#ÊB φ(V, EW ∪ ÊB
).
¸µ ¡< H¹
lË<̵¶ H¹
µ¸¡K—’¹
ɄÈFȸ!Ï tw|€l1{¸}„ƒ¯}„€|}fƒ¯l£}„j1n>n<Åz‘1‚„nK{„{„ƒ¯ˆI1{4ˆ’° ⠐Š1“ b̂ ƒ­ µ² z¹9»Û¼;nmQÐy¼)‚„ƒ¯}„ndn<º­n<mon<H} Š‚fy¼;n<ƒ­£’jH}f{g’{
°²ˆIº¯º­ˆ#¼4{ #
XX Y
Y Y a
µ¸¡ Ô ¹
a
1−
.
φ(t) =
b
b
b
.CÅz‘¢Š1“|ƒ­l£d}„j1nº­I{~}}„nF‚„m»l£In}
Y a
XX Y
X
Y a
µ¸¡ÐҊ¹
b
(−1)
φ(t) =
.
ÊB ∈P(EB )
s
v∈V iv =1 k∈V
i` ik
ik i`
ik
(k,`)∈EW
i`
ik
(k,`)∈EB
s
#ÊB
ik
ik i`
i` ik
bi`
bik
ij1n”<ˆ’1”<º¯€¢{¸ƒ­ˆ’°²ˆ’º­º¯ˆ#¼4{F¿
t€l‚Õ‘1‚„ˆ|”nK“|€l‚fn;”<ˆ’mo‘l€|}fnF{ՈI‚f“ln<‚Õ”ˆI1“|ƒ¯}„ƒ­ˆ’1{ƒ­d‚fnF”<€l‚ {¸ƒ­®’n¼;ÐyI¿ ¬ jlnFo“lnFŠº­ƒ­l£¼)ƒ±}fjdˆI‚f“ln<‚ p »
’º¯ºl”<ˆ’1“lƒ±}fƒ¯ˆI1{ՈŠ°Ûˆ’‚ “|nF‚{„mQŠº­º¯nF‚}„j1’ j1Ð®’n)Šº­‚fnF’“ly¢nFn<Qˆ’l}fŠƒ­lnK“ï¿ ¬ n4”º†Šƒ­m•}fj1#}CnFº¯nFmon<H}f’‚„y
¼;n<ƒ­£’jH} {rˆ’°Šº­º+}„n<‚fmQ{gƒ¯›}„jlno{„€lm ˆŠ°opµ~¡Ð—’¹rj1Ð®’ndŠº­‚„nK’“|yXn<n<ƔˆImo‘l€|}„nK“ï»ïnÅl”<n<‘|}Qµ²‘Xn<‚fj1Š‘¢{f¹4°²ˆ’‚
}fjlnˆ’1ngƒ­¼)jlƒ†” j Ê = E ¿
3Ã1“ln<nF“&»#°²ˆ’‚ՐŠdŠ‚ ”Cƒ­
}fjln<>}„j1n<‚fnƒ†{+lˆg”<ˆ’H}„‚fƒ­l€|}„ƒ­ˆ’°²‚„ˆIm }„jlƒ†{Š‚ ” #&¼;n;”FŠ>ƒ¯H}„nF‚„‘1‚„n<}
}fjlƒ†{4’{Q“|nFº¯n<}„ƒ­ˆ’ˆŠ°+}„jlnEŠ‚ ”Š\»|ʼ)j1n<‚fnF’{°²ˆI‚4ŠŠ‚ ”tƒ­ E »1¼nˆ’l}fŠƒ­}„jln€1{„€1’ºÛ}„n<‚fm µ²°²ˆ’‚4dlˆI
‘¢Š‚„}„ƒ¯}„ƒ­ˆ’lnK“Q)€ll£In·Ã¾€|}¸} {„” jlnFmonK¹n<Å|”<n<‘|}}„j1Š}*¼;n4j1Ð®’n ƒ­1{¸}„nK’“Qˆ’° »H¼)jlƒ†” j'¼;ˆ’€lº†“o¢n
}fjln€1{„€1ŠºÛ}fn<‚fm8ƒ¯°+}fjln“|ƒ­‚„nK”9}„ƒ­ˆ’ˆ’°+Š‚ ”gj1’“Xn<n<S” j1’l£’nK“ïa»lŠ¢“’I{„{„ˆ|”ƒ†a#}„nK“'moƒ­H€¢{4{¸ƒ­£’?¿
prˆ#¼ ¼)jln<lnF®’nF‚tŠO’‚f”ƒ­{g“|n<º­n}fnF“ï»Ûƒ¯}t°²ˆIº¯º­ˆ#¼4{)°²‚fˆ’mº­n<momQžo}fj1#}t}„jln>}„nF‚„mӃ†{r
‘l‚„ˆ|“|€¢”9}tˆ’°
nFº¯nFmon<H}f’‚„y¼;n<ƒ­£’jH}f{t°²ˆ’‚g}f‚„nFnF{gˆŠ°{¸mQŠº­º­n<‚{¸ƒ­³<nI¿dwzƒ­1”nolº†’” Àlˆ|“|nK{j1Ð®’ndŠº­ºÕ‘¢ˆH{„{„ƒ­lº¯ndº­ˆz”F#}fƒ¯ˆI1{<»
¼)j1n<O”ˆI1“|ƒ¯}„ƒ­ˆ’1{rˆŠ°CˆI‚f“|nF‚4º­nF{f{)}„j¢Š p ¼;n<‚fn”ˆImo‘l€|}„nK“ï»X}„jln>n<º­n<mon<H} Š‚fy
¼;n<ƒ­£’jH}f{)ˆ’°*Šº­º&}f‚„nFnF{4ˆ’°
¼;n<ƒ­£’jH}4#})moˆI{¸} p − 1 j1Ð®IntXn<nFS”ˆ’mo‘l€l}„nF“&¿
ijlnƒ­“lnFSlˆ#¼ ƒ†{}fˆO°²ˆ|”€¢{ˆ’½}„jlƒ†{£I‚f’‘ljƃ¯H}fn<‚f‘l‚„n<}fŠ}„ƒ­ˆ’&»ÕŠ¢“=”ˆ’mo‘l€l}„n
ˆI‚f“ln<‚”<ˆ’1“|ƒ¯}„ƒ­ˆ’¢{<¿
¬ ƒ¯}„jÄ}fjln<ˆI‚„nFm ¡K—¼;nQ”<ŠƐ’{f{„ˆz”<ƒ­Š}„nd}„ˆSlƒ¯·Ã”ˆ’º­ˆ’€1‚„nK“S}„‚fn<nQS{„n}ˆ’°ˆ’‚fƒ­n<H}„nK“´£I‚f’‘lj´°²ˆI‚¼)jlƒ†” j
ˆIlº­y'ˆ’lnƒ†{)”ˆIllnF”}„nK“ï¿Õijlƒ­{4mdnKŠ1{;}fj1#})¼;n”<Šˆ’l}fŠƒ­ˆ’‚ “|n<‚)”ˆI1“|ƒ¯}„ƒ­ˆ’1{;°²ˆI‚)nFI” j
£’‚ Š‘1j&¿ n}
€¢{lˆ#¼ } ŠÀ’n>}„€l}„ˆ’‚fƒ†Šº&{fŠmo‘lº­nt}fˆ'{¸j1ˆ#¼ jlˆ#¼ }fjlnmon}fjlˆ|“¼;ˆ’‚fÀ|{<¿
v∈V iv =1 k∈V
ÊB ∈P (EB )
B
(k,`)∈EW
(k,`)∈ÊB
B
B
B
B
i` ik
õõó: !;
ik i`
¡KŸ
ïÉ Î É9Ê
k
5ƒ¯·Ã”ˆ’º­ˆ’€1‚„nK“Q}„‚fn<n
ʶË<Ì
l
j
t
i
φ(t) = 1/γ(t)
φ(t) = φ(t1 ) − φ(t2 )
s
X
bi âij ajk ajl =
i,j,k,l=1
s
X
bi bj ajk ajl
1
12
−
i,j,k,l=1
bj aji ajk ajl
i,j,k,l=1
k
l
j
t1 − t 2
s
X
k
l
j
−
i
i
i j1ƒ­{*nÅlŠmo‘lº­n4ƒ­º¯º­€1{¸}„‚ #}„nK{Õjlˆ#¼¼;n4”<Š'”<ˆ’mo‘l€|}fn4ˆ’‚ “|nF‚*”ˆI1“|ƒ¯}„ƒ­ˆ’1{’{f{¸ˆ|”ƒ†#}fn¼)ƒ±}fj
£’ƒ­®’nFo£’‚ Š‘lj&¿
j1’{t}~¼;ˆ”<ˆ’l1nF”9}fnF“›”ˆ’mo‘Xˆ’lnFH}f{F»?Š1“´jln<¢”n’» ) ƒ†{t}„jlno‘l‚fˆ|“|€1”9}ˆ’°*}~¼;ˆ
nFº¯nFmon<H}f’‚„y¼;n<ƒ­£’jHt}f{?}„j1Š}+j¢Ð®’n*Xn<nF>’º¯‚fnFI“|ytnF®#Šº­€1#}fnF“
µ¶¼)jln<d“|nF’º¯ƒ­φ(t
l£g¼)ƒ±}fjˆI‚f“|nF‚ g”ˆI1“|ƒ¯}„ƒ­ˆ’1{ ¹9¿
ij1n<‚fn°²ˆ’‚fn;}„jln4}„jln41n<¼”ˆ’¢“|ƒ±}fƒ¯ˆI'”<’dXn4n<Åz‘1‚„nK{„{„nF“d’{ φ(t ) = φ(t ) − 1/γ(t) »’¼)jlnF‚„n}fjln4‚fƒ¯£IjI}„·
j¢Š1“z·Ã{„ƒ­“|n‚fnF“|€¢”nF{;}fˆQ>‚ #}fƒ¯ˆI1ŠºÛz€lm>¢nF‚F¿
prˆŠ}„n'}„j1Š}>¼;nQ£’n}>ˆ’ln'ˆI‚f“ln<‚”<ˆ’1“|ƒ¯}„ƒ­ˆ’ư²ˆI‚nFI” jƝlƒ±·Ã”ˆIº¯ˆI€l‚fnF“Ä£’‚ Š‘lj?¿wzƒ­1”<nQ}„jlnF‚„nŠ‚fnQº¯nK{„{
ˆI‚„ƒ­n<H}fnF“S°²‚fn<nd}f‚„nFn>}fj1ŠƝlƒ¯·Ã”ˆ’º­ˆ’€1‚„nK“}„‚fn<nK{<»&ƒ¯}mQÐySj1Š‘1‘¢nF›}fj1#}’º¯º}fn<‚fmQ{tˆ’°}fjlnQ{¸€1mvƒ­µ¸¡K—’¹
’‚„nŠº­‚fnF’“lyQn<®#Šº­€1#}fnF“ï¿
Ø ½©|¨ ijln
’¢ˆ#®In'“|ƒ†{„”<€1{„{„ƒ­ˆ’ƃ­{>”ˆIjln<‚fn<H}¼)ƒ¯}„j=}fjln'‚fnF{„€lº±}>ˆŠ°4sS€l‚f€1 "­¡K%ž $ #>}„jlnF‚„n'n<Å|ƒ­{¸}>’{
mQ’Hy>ˆ’‚ “|n<‚”ˆ’¢“|ƒ±}fƒ¯ˆI1{Õ°²ˆI‚Cˆ’‚ “|nF‚ n °²ˆ’‚’Qƒ¯H}fn<£’‚ #}fƒ¯ˆIomdn<}„jlˆ|“QI{Õ}„j1n<‚fn4nÅ|ƒ†{~}*ˆI‚„ƒ­n<H}fnF“>°²‚„nFn)}„‚fn<n
¼)ƒ¯}„j n lˆ|“|nF{F¿
! (45)>+:-d5(=+)!(9!
3Û}„j1ƒ­{n<Å|’mo‘lº¯n
1
1
2
1
ƒ­®’nF'‚fˆHˆ’}„nK“olƒ­”<ˆ’º­ˆ’€l‚fnF“d}f‚„nFn »z“|n<lˆ’}„ntHy
}fjlntI{„{„ˆ|”ƒ†#}„nK“>°¶Šmoƒ¯º­yoˆŠ°&ˆI‚„ƒ­n<H}„nK“d£’‚ Š‘lj1{
ˆI|}f’ƒ¯1nF“zy
}fjlnˆI‘¢nF‚fŠ}„ƒ­ˆ’SˆŠ°ÕnFƒ±t}fjln<‚g”€|}„}„ƒ­l£
(σˆ’°,‚fgn<®I)n<‚ {¸ƒ­l£olº†’” ÀnF“|£InF{ “|nFlˆŠ}fnzy g }fjln>”<ˆ’|·
1nF”9}fnF“n<º­n<monFI})ˆ’°}fjlƒ­{)°¶Šmoƒ¯º­y´µ¶ˆ’|} Šƒ­lnF“zy'‚fn<®’nF‚f{„ƒ­l£>’º¯ºïlº†’” À
nF“|£InF{ ¹9¿ n<} v(g) Xn}„jln®#Šº­€ln
ˆ’°r}fjlnSn<º­n<mon<H} Š‚fy=¼nFƒ¯£IjI}K»{¸}„ˆI‚„nK“½ƒ­Æ“lŠ}f›1’{„nS”<’º¯º­nF“ © ¿ wzƒ­1”nSnFº¯nFmon<H}f’‚„y=¼;n<ƒ­£’jH}f{
i
i i
∗
²î+õ Ñö
ɄÎIËÉ È#ÌÛÎ#ÊÑÍ[ʶÈ#Ì*ÏÈ#É
¯Ë<ÍMÊ
ŠÉ9ÍMÊÑÍ[ʶÈ#ÌÛË Î|ÌHË|Í[ÍlËQË
¡I¡
’‚„n
‚fŠ}„ƒ­ˆ’1’ºCz€lmXn<‚ {<»Õ¼;n
{¸}„ˆI‚„nQ}fjln<m I{´‘1Šƒ­‚ˆ’°4ƒ¯‚f‚„nK“|€1”<ƒ¯lº­n
ƒ­H}„nF£’n<‚ {F¿Oijln‘1‚„ˆ|”nK“|€l‚fno°²ˆ’‚
ˆI|}f’ƒ¯1ƒ¯l£oˆI‚f“|nF‚)”ˆ’¢“|ƒ±}fƒ¯ˆI1{)#}4d£’ƒ­®’nFˆI‚f“ln<‚ n mQÐy'Xn¼)‚fƒ±}„}„n<’{
o¨#Þ Ø ¨ Ø µ¶¢¹
?Ý ¨
ˆ’‚ “|nF‚„nK“'zy
ƒ­1”‚fnFI{¸ƒ­l£o“|n<£I‚„nFn’»z¼)ƒ¯}„j
∈ BT
Ý t § Ø>}fjln°¶Šmoƒ­º¯yˆŠ°ˆ’‚fƒ¯nFI}fnF“
£’‚ Š‘lj1{ (σ , g )|t| ≤ n
© +©|§ ØQn<º­n<mon<H} Š‚fy'¼nFƒ¯£IjH}f{;ˆŠ°ÕŠº­ºïlˆŠ·Ã”ˆIllnK”9}„nK“ g
©
g ∈Ø ” j1nF” À'ƒ¯° P σ v(g ) + σ v(g ) = 1/γ(t)
M¦ Ø © ← © ∪ {v(g ) = σ (1/γ(t) − P σ v(g ))}
+Þ ?Ý ¨
Ø § ¨ © Ø ½©|¨ tijln ’€1{f{~· nF£’n<¢“|‚„ng)€ll£In·Ã¾€|}¸} dmdn<}„jlˆ|“Ƶ²ƒ­S¼)jlƒ†” j a Š1“ b Š‚fngnK…H€1ŠºÛ}fˆ ⠐Š1“
¹ƒ­{Àzlˆ#¼)›}„ˆSXn'{¸yzmo‘lº­nF”}„ƒ†”Š»?ˆ’°;ˆI‚f“ln<‚ ¼)j1n<‚fn ƒ­{g}fjln'z€lmXn<‚ˆ’°{¸}f’£’nK{<¿QijlnF‚„n<°²ˆ’‚fn>}fjln
b̂
” j1nF” ÀˆŠ°”<ˆ’mo‘1Š}„ƒ­lƒ¯º­ƒ¯}~yO”ˆ’¢“|ƒ±}fƒ¯ˆI1{t}„j1Š}gˆ|2s”F”€l‚ {r¼)j1n<s g Šº­‚fnF’“lyXn<º­ˆ’l£H{t}„ˆ © j¢’{r}„ˆ¢n
{f#}fƒ­{ 1nK“ï¿ ¬ n€1{„nƒ¯}4ˆ’lº­y
’{)o“|nFl€l£I£’ƒ­l£}fˆzˆ’º[¿
t€l‚rƒ¯mo‘lº­n<mon<H} #}„ƒ­ˆ’€1{„nF{)}„jln‘l‚fˆ’£I‚f’m °²ˆ’‚)}f‚„nFnF{)£In<lnF‚fŠ}„ƒ­ˆ’ˆŠ° ƒ?’1“)€1{„À’n<y "¯¡I¡ $Á¿ijlnF
¼;nl€lƒ­º­“o’º¯º|‘XˆI{f{¸ƒ­lº­n4”ˆIº¯ˆI‚fŠ}„ƒ­ˆ’1{?°²ˆ’‚Cg£’ƒ­®’nF}f‚„nFn’»’’1“oŠ‘l‘1º¯y‘1‚„ˆ|”nK“|€l‚fn Q¨#Þ Ø ¨ Ø ¿Õijln)”ˆ|“|n
ƒ†{¼)‚fƒ±}„}„n<ƒ­}„jln>Ž º­’l£’€1’£’nI¿+ijln”<ˆ’º­ˆ’‚ #}fƒ¯ˆI
‘l‚fˆ|”nF“l€l‚„ngƒ†{4’{;°²ˆIº¯º­ˆ#¼4{ #
Ý Ý ¨ Ø µ²} ¹
‚fn}„€1‚„9 " ¬ » 5+$
M¦t 2
Ø =.CτÅz‘1‚„nK{„{ t ’{º­ƒ¯lÀ|{°²‚fˆ’m"‚„ˆzˆŠ}}fˆd1‚f’1” jlnF{ t = [t ...t $
Ý ¨ ƒ >¡g}„ˆ
Žˆ’º­ˆ’‚ n<?µ ¹
Ý ´
¥ ØoŠsº­º&”ˆIº¯ˆI‚fŠ}„ƒ­ˆ’
tƒ­H}„ˆQi+Š
Ø § ¨ µMi’¢¹
∗
i
i i
i
∗
i
i
∗
i
∗
∗
∗
i
i
i
∗
n1
1
ns
s
i
“ln<lˆ’}„n'}fjln
z€lm>¢nF‚>ˆ’°r”ˆIº¯ˆI‚fŠ}„ƒ­ˆ’ƈа4l‚ Š¢” j t ¿Äijln<=}fjln
z€lm>¢nF‚>ˆ’°
Ø ½©|¨ n<}
”<ˆ’º­ˆ’‚ #}fƒ¯ˆI1{*ˆ’° t nƒ­{ (t )
c
i
i
s
Y
nc (ti )!
ni !(nc (ti ) − ni )!
i j1ƒ­{n<Åz‘1º­’ƒ¯}„jln”<ˆ’m>lƒ¯¢#}„ˆI‚„ƒ†ŠºÛn<Åz‘1º¯ˆH{¸ƒ­ˆ’Š1“}„j1nº¯ƒ­moƒ±} #}„ƒ­ˆ’ˆ’°ˆ’€1‚)‚„nK{¸€lº¯}}„ˆ n = 7 ¿
Ø ½©|¨ tsSˆH{~})ˆ’°?}„j1n”ˆ’mo‘l€l}„ƒ­l£d}„ƒ­mdnƒ†{4{„‘¢nFI})ƒ­}„jln{„nFŠ‚ ” j
ˆ’°+d£’ƒ­®’nF£I‚f’‘lj›µ¶¼)jln<1n<®’nF‚)ƒ±}
n<Å|ƒ­{¸}f{ ¹Õƒ­Q}„j1nr1’lÀÛ¿Ct€1‚*ƒ­md‘1º¯nFmdnFH}f#}fƒ¯ˆI
” jlnK” Àz{*lƒ ‰~nF”9}fƒ¯ˆI'¼)ƒ¯}„j’º¯º¢£’‚ Š‘lj¢{CˆŠ°?{fŠmon4z€lm>¢nF‚ˆ’°
®In<‚„}„ƒ†”nK{)Š1“nK“|£’nK{<¿)tH®zƒ­ˆ’€1{„º­yQ}„jlƒ†{4‘1‚„ˆ|”nK“|€l‚fn”ˆI€lº­“Xnƒ­mo‘l‚„ˆ#®InF“ï+¿ 4lˆ’‚rˆ’‚ “|n<‚4£I‚„nK#}fn<‚}„j¢Š Ô »
”<ˆ’1“lƒ±}fƒ¯ˆI1{r’‚„n}fˆzˆ
z€lmon<‚fˆ’€1{)}„ˆ
Xn>“lƒ­{„‘lº†Ðy’nF“&¿rÖ4ˆ#¼;n<®In<‚K»lƒ±}”<’S¢n>€1{„n°²€lº+}„ˆ'°²ˆ’‚fm€lº†#}fn}„jlnFm
ƒ­ˆ’‚ “|n<‚;}fˆo”<ˆ’mo‘l€|}fn}„jlnˆI‚f“ln<‚;°²ˆ’‚4{„‘¢nK”ƒ ¢”®#Šº­€lnF{ˆ’° a ’1“ b ¿
i=1
õõó: !;
¡Ðž
ïÉ Î É9Ê
i+Šlº­no¡ #Õp4€lm>¢nF‚4ˆŠ°ˆ’‚ “|nF‚”<ˆ’1“lƒ±}fƒ¯ˆI1{
t‚ “|n<‚
¡ ž ˜ — Ô
wzƒ­mo‘lº¯n
¡ ¡ ž ™ ž’Ÿ
wzyzmo‘lº­nF”9}fƒ­” ¡ ¡ ˜ š žHÒ ™1¡
‡ÕŠ‚„}„ƒ¯}„ƒ­ˆ’lnK“ ž ¡F —’ž žl¡< ™1¡ Ô
ʶË<Ì
Ò
H š
˜H—ŠŸ
1¡I¡ Ô
¬ n ¢‚f{¸}“|ƒ†{¸‘1º­Ðyt}„jlnz€lmXn<‚ˆŠ°lˆ’‚ “|nF‚?”ˆI1“|ƒ¯}„ƒ­ˆ’1{ƒ¯} Šlº­nr¡I¿t1{„n<‚f®’n}fjln*‚ Š‘1ƒ­“ƒ­1”‚fnFI{¸nՈ’°|}„jlnK{¸n
z€lm>¢nF‚f{¼)ƒ¯}„j p °²ˆ’‚4{„yzmo‘lº¯nK”9}fƒ­”{f” jln<monK{<»l’1“
nF®’nFmoˆ’‚fn¼)ƒ±}fj£’n<1n<‚ ŠºÛ‘1Š‚„}„ƒ¯}„ƒ­ˆ’1nF“{„” jlnFmonF{F¿
P
3Ão}fjln)lnÅz}C} Šlº­nF{F»’¼;n)€1{„n}„jln)€¢{¸€1’ºllˆ’}f#}fƒ¯ˆI1{
Š1“ P ¿ º¯º1ƒ¯1“lnÅ|nF{
ƒ­}„j1n>{„€lmӮА’‚„y
°²‚„ˆIm ¡}„ˆ s »X¼)jln<‚fn s ƒ­{r}„jlnH€1mXdn<‚g=ˆŠ°C{~} Šb£IanF{)ƒ­S}„j1cn>)=€l1£’n·Ãa¾€|}„}fomon}fjlˆ|“ï¿
º¯º;µ²º†#}fnÅS{¸ˆI€l‚f”<nˆŠ°f¹)}f’lº¯nK{r’‚„n£In<lnF‚fŠ}„nF“’€|}„ˆImoŠ}„ƒ†”<’º¯º­y
zy}„jln”<ˆ’mo‘l€|}fn<‚g”ˆ|“|n’¿rŽˆI1“|ƒ¯}„ƒ­ˆ’1{
°²ˆI‚)ˆ’‚ “|n<‚¡g}„ˆQ o¼;n<‚fng’º¯‚fnFI“|yQˆ’l}fŠƒ­lnK“zy
ÖtŠ£In<8‚ "ã,Ò $[¿
i+Šlº­nž #*t‚ “|‚„no¡
‚ Š‘lj ŽˆI1“|ƒ¯}„ƒ­ˆ’
j
X
i i ij
i
j
ij
bi = 1
i+Šlº­n˜ #*t‚ “|‚„nž
‚f’‘lj Žˆ’1“|ƒ¯}„ƒ­ˆ’
X
‚f’‘lj
dj =
+i Šlº­n #*t‚ “|‚„n˜
ŽˆI1“|ƒ¯}„ƒ­ˆ’
‚f’‘lj
X
cj d j =
1
6
X 1
1
d2 =
bk k
3
1
2
Žˆ’1“|ƒ¯}„ƒ­ˆ’
X
bi c2i =
1
3
²î+õ Ñö
ɄÎIËÉ È#ÌÛÎ#ÊÑÍ[ʶÈ#Ì*ÏÈ#É
¯Ë<ÍMÊ
ŠÉ9ÍMÊÑÍ[ʶÈ#ÌÛË Î|ÌHË|Í[ÍlËQË
‚ Š‘lj
+i Šlº­n—#*t‚ “|‚„n Žˆ’1“|ƒ¯}„ƒ­ˆ’
‚ Š‘1j
X 1
1
alk dk dl =
bk
8
X bi
5
aik ci dk =
bk
24
X
1
c2j dj =
12
X 1
1
ck d2k =
bk
12
‚ Š‘1j
õõó: !;
Žˆ’1“lƒ±}fƒ¯ˆI
1
24
X
ajk dj ck =
X
bi aij ci cj =
X
bi c3i =
1
8
1
4
X 1
1
d3 =
b2l l
4
i+Šlº­n Ô #*t‚ “|‚„n—
Žˆ’¢“|ƒ±}fƒ¯ˆI
‚ Š‘1j
X bi
3
aik ail dk cl =
bk
40
X
1
alk akj cj dl =
120
X bi
11
alk ail ci dk =
bk
120
X bi bj
2
ajk aik ci cj =
bk
15
X 1
1
aml alk dk dm =
bk
30
X 1
1
alm d2l dm =
bl bm
15
X 1
1
aml d2l dm =
b2l
10
X 1
7
alk dk cl dl =
bk
120
X 1
1
alk ck dk dl =
bk
40
¡K˜
Žˆ’1“lƒ±}fƒ¯ˆI
X
bi aik aij cj ck =
1
20
X
bi aik akj ci cj =
1
30
X bi
3
alk aik ci dl =
bk
40
X bi
2
aim ail dl dm =
bl bm
15
X 1
1
amk alk dl dm =
bk
20
X 1
1
akl d2k cl =
bk
60
X bi
3
ail ci d2l =
b2l
20
X
1
ajk cj dj ck =
40
X bi
7
aik ci ck dk =
bk
120
¡F ïÉ Î É9Ê
ʶË<Ì
i+Šlº­n Ô #*t‚ “|‚„n—
Žˆ’¢“|ƒ±}fƒ¯ˆI
‚ Š‘1j
‚ Š‘1j
X bi
3
aik c2i dk =
bk
20
X
1
akj c2j dk =
60
X
1
c3j dj =
20
X 1
1
c2 d 2 =
bk k k
30
X 1
1
d4m =
3
bm
5
Žˆ’1“lƒ±}fƒ¯ˆI
X
bi aij c2i cj =
1
10
X
bi aij ci c2j =
1
15
X
bi c4i =
1
5
X 1
1
cl d3l =
b2l
20
!"#
X bi
7
aim ail anl dm dn =
bl bm
144
X bi
19
aim aik alk dl cm =
bk
720
X 1
1
anm amk alk dl dn =
bk
72
X bi
13
aim amk alk ci dl =
bk
360
X 1
1
ank amn alm dk dl =
bk
144
X
alm akl ajk dj cm =
1
720
X bi
7
aik amk alm ci dl =
bk
360
X bi
1
aim ain anl dl dm =
bl bm
18
X bi
13
aim ail alk dk cm =
bk
360
X bi
19
aim aml alk ci dk =
bk
720
X bi
7
alm ail aik dk cm =
bk
360
X bi bj
61
ajm aim ail cj dl =
bl bm
720
X bi bj
7
ajl ajk aik ci cl =
bk
144
X
bi akl ajk aij ci cl =
1
144
X bi bj
bk
X
ajl alk aik ci cj =
bi ail aik akj cj cl =
1
18
1
72
²î+õ Ñö
ɄÎIËÉ È#ÌÛÎ#ÊÑÍ[ʶÈ#Ì*ÏÈ#É
¯Ë<ÍMÊ
ŠÉ9ÍMÊÑÍ[ʶÈ#ÌÛË Î|ÌHË|Í[ÍlËQË
¡Ð—
X bi
7
aik ail dk c2l =
bk
180
X
alk akj c2j dl =
1
360
X bi
13
alk ail c2i dk =
bk
180
X
bi aik aij c2j ck =
1
36
X
bi aik akj ci c2j =
1
72
X
bi ajk aij c2i ck =
1
36
X bi
19
alk aik c2i dl =
bk
360
X bi bj
7
ajk aik c2i cj =
bk
72
X bi
13
aim ail cl dl dm =
bl bm
360
X bi
1
ail aik ck dk cl =
bk
45
X 1
1
aml alk ck dk dm =
bk
180
X bi
7
ail alk ci ck dk =
bk
360
X 1
7
amk alm dk cl dl =
bk
360
X 1
1
amk alk cl dl dm =
bk
45
X
7
bi
ain aim d2m dn =
b2m bn
72
X 1
1
anm aml d2l dn =
b2l
36
X
akl ajk cj dj cl =
1
180
X bi
13
aik alk ci cl dl =
bk
360
X bi
19
aim ail d2l cm =
b2l
360
X bi
13
aim aml ci d2l =
b2l
180
X
1
1
anm aln d2l dm =
bl bm
72
X 1
1
alm akl d2k cm =
bk
360
X
1
1
anm alm d2l dn =
bl bm
36
X bi
7
aim alm ci d2l =
bl bm
180
X 1
1
alk alm dk dl cm =
bk
60
X 1
1
amk akl dk cl dm =
bk
240
X bi
7
aik ail ci dk cl =
bk
120
õõó: !;
!"#
X
ajl ajk dj ck cl =
1
120
X bi
1
aik akl ci dk cl =
bk
80
X
bi aik aij ci cj ck =
1
24
¡Ô
ïÉ Î É9Ê
ʶË<Ì
X
alj ajk cj ck dl =
1
240
X bi
11
alm ail ci dl dm =
bl bm
240
X bi bj
ajl ail ci cj dl =
b2l
11
120
X bi
1
alk aik ci ck dl =
bk
60
X
1
1
anm anl dl dm dn =
bl bm
24
!"#
X
bi aij ajk ci cj ck =
1
48
X bi
7
aml ail ci dl dm =
b2l
120
X bi
11
alk ail ci dk cl =
bk
240
X bi bj
1
ajk aik ci cj ck =
bk
24
X
1
1
anl alm dl dm dn =
bl bm
48
X bi
11
aim ail ci dl dm =
bl bm
120
X 1
1
aml alk dk cl dm =
bk
80
X 1
1
anl aml dl dm dn =
b2l
24
X 1
1
amk alk ck dl dm =
bk
120
X
1
b2m bn
amn d3m dn =
1
24
X 1
1
anm d3m dn =
b3m
12
X
1
1
alm cl d2l dm =
bl bm
40
X 1
1
aml cl d2l dm =
b2l
60
X 1
1
alk dk c2l dl =
bk
30
X 1
1
alk c2k dk dl =
bk
120
X bi
7
aik c3i dk =
bk
60
X
akj c3j dk =
1
120
X 1
1
alm d3l cm =
b2l
120
X bi
7
aim ci d3m =
b3m
60
X 1
1
akl ck d2k cl =
bk
120
X bi
1
ail ci cl d2l =
b2l
30
X
ajk c2j dj ck =
1
60
X bi
1
aik ci c2k dk =
bk
40
X
bi aij c3i cj =
1
12
X
bi aij ci c3j =
1
24
²î+õ Ñö
ɄÎIËÉ È#ÌÛÎ#ÊÑÍ[ʶÈ#Ì*ÏÈ#É
¯Ë<ÍMÊ
ŠÉ9ÍMÊÑÍ[ʶÈ#ÌÛË Î|ÌHË|Í[ÍlËQË
¡#Ò
X
!"#
1
1
alm d2l cm dm =
bl bm
90
X 1
1
akl d2k c2l =
bk
180
1
1
anm d2m d2n =
b2m bn
18
X 1
2
aml d2l cm dm =
b2l
45
X
X bi
19
ail c2i d2l =
b2l
180
X
1
90
ajk cj dj c2k =
X
bi aij c2i c2j =
X
bi c5i =
1
18
1
6
X 1 2 3
1
c l dl =
b2l
60
X 1
1
alk ck dk cl dl =
bk
72
X bi
2
aik c2i ck dk =
bk
45
X
c4j dj =
1
30
X 1 3 2
1
c d =
bk k k
60
X 1
1
cm d4m =
b3m
30
X 1 5
1
d =
b4n n
6
,)#);)!
"¯¡ $¢¿ i¿ 5;n}„}f{F¿ wz€l‚f®’nFyLj’°H€1mdnF‚„ƒ†”<’º4mon}fjlˆ|“l{'°²ˆI‚'}f‚fŠ‰~nF”9}fˆ’‚fyLjI‘|}„ƒ­moƒ¯³K#}„ƒ­ˆ’?¿ ȸÏ
|ʶΠ#Ì Ë È#Ì¢ÍMÉ„È ŠÌXÎ#Ì >Ê9»¢ž|¡#¯¡K™’˜|žŠŸzÒz»l¡K™’™Išl¿
" ž $¢¿ i¿ 5;n}„}f{F¿ ;É <ÍMÊQËÍ lÈKÎ,?ÏȊÉgÈ ÛÍMÊ È#Ì¢ÍMÉ„È 9ÊÑÌÌXÈŠÌ ÊÑÌÛË#É ÛÉ„È ŠÉ dÊÑ̒¿ïwzˆ|”ƒ­n}~y
°²ˆI‚ 3Ã1“|€1{¸}„‚fƒ†Šº?’1“ ‘l‘lº­ƒ¯nK“sO#}fjln<mQ#}fƒ­”F{µMw 3 s´¹9»¢‡Cjlƒ­º­I“|n<º­‘ljlƒ†l»1‡ »¢žŠŸIŸl¡I¿
" ˜ $¢¿ 40¿ 5;ˆ’l¢Š1{4’1“ ¿ ’€l1Ðy’¿ Š‚f£’n{f”<Šº­n“|ƒ­‚fnF”9}tˆ’‘|}fƒ¯mQŠº+”ˆII}f‚„ˆIº?Š‘1‘lº¯ƒ­nF“S}„ˆ

‚„n<·[nFI}f‚„y
‘1‚„ˆIlº¯nFm¿ ÕÈ~Ï |ʶΠŠÌ Ë È#Ì1Í[É„È #ÌÛÎ #Ì dÊ »Xž|¡ # ™’™ Ô Û¡FŸ’ŸIŸl»l¡K™’™Išl¿
" $¢¿Žt
¿ 5;€|}f” jlnF‚F¿ oË<É9Ê oË<Í lÈK,Î ÏȊÉȊɄÎ#ÊÑ0Ì #É •Î#Ê rË<ɄË<Ì1Í[Ê '
Ë ŠÍMÊ¶ÈŠÌ 9¿ ¬ ƒ¯º­n<yH·
3ÃH}fn<‚ {„”<ƒ¯nF1”n‡C€l1º¯ƒ†”<Š}„ƒ­ˆ’&¿ Iˆ’j1 ¬ ƒ­º¯nF
y wzˆ’1{ }f“ï¿­»XŽjlƒ†” jlnK{~}fn<‚K»1žŠŸ’ŸI˜l¿
" — $ ¿ ¿¢ŒrˆII} ” jln<®Š1“ ¬ ƒ¯º­º­ƒ­’m ¬ ¿1ÖtŠ£’nF‚F¿Õijln .*€lº¯nF‚4Š‘l‘1‚„ˆÐÅ|ƒ­moŠ}„ƒ­ˆ’ƒ¯S{¸}f#}fn”ˆ’¢{~}f‚f’ƒ¯lnK“
ˆI‘|}„ƒ­mQŠº?”<ˆ’H}„‚fˆ’º[¿ ŠÍ È ­»ïÒ#Ÿ1µ[žŠ˜I˜I¹ #­¡ÐҊ˜ lžŠŸ’˜1»IžŠŸIŸl¡I¿
õõó: !;
¡Kš
ïÉ Î É9Ê
ʶË<Ì
¿ŒrˆII} ” jln<®Û» ¬ ƒ¯º­º­ƒ­’m ¬ ¿ÖtŠ£’nF‚F»Š1“ rº†’“lƒ¯moƒ­‚'s=¿ Õn<º­ƒ­ˆ#®X¿ wznK”ˆI1“z·Áˆ’‚ “|n<‚o)€ll£In·
¾€l}¸}fQ’‘l‘l‚fˆÐÅ|ƒ¯mQ#}fƒ¯ˆI1{;ƒ¯O”<ˆ’H}„‚fˆ’ºï”<ˆ’1{¸}„‚ Šƒ­lnF“ˆ’‘l}„ƒ­mo’º?”ˆIH}„‚fˆ’º[¿ QËÉ Ì ­»
˜Iš1µ¸¡K¹ #㞊ŸHž lž’ž Ô µ²n<º­nF”}„‚fˆ’lƒ†”F¹»1žŠŸIŸ’Ÿl¿
"ãÒ,$ ¬ ¿ÛÖr’£’nF‚F¿)€ll£In·Ã¾€|}¸} omon}„j1ˆz“1{rƒ­Sˆ’‘|}fƒ¯mQ’º”ˆIH}„‚fˆ’º?’1“}„jln}„‚ Š1{¸°²ˆ’‚fmonF“’“Љ~ˆIƒ¯H}t{„y|{~·
}fn<m¿ oË<É ŠÍ ­»¢šzÒ|µMžI¹ #ãž# |Ò |ž’šIžl»IžŠŸIŸ’Ÿl¿
" š $8.r¿Ör’ƒ¯‚fn<‚K»4Žt¿ €1lƒ­” j?»’1“ ¿ ¬ Šl1n<‚K¿ Ë È QËÍ[É9Ê OÌ oË<É9Ê 4ÊÑÌ¢ÍÁË ’É #Í[ʶÈ#Ì¢¿ wz‘1‚„ƒ­l£’nF‚¸·
ÕnF‚„º†Š£¢» 5;n<‚fº¯ƒ­&»1žŠŸIŸIžl¿
" ™ $8.r¿CÖrŠƒ­‚fn<‚K»wۿՇ¿Cp ‚f{„n}„}F»*’1“ ¿ ¬ Šl1n<‚K¿ È FÊÑÌ ÇÈ#ÉfÎ#ÊÑÌ ŠÉ ½Î#Ê rË<ɄË<Ì1Í[Ê r
Ë ŠÍMÊ¶ÈŠÌ 9¿
w|‘l‚„ƒ­l£In<‚„· ÕnF‚„º†Š£¢» 5;n<‚fº­ƒ¯&»¢{¸nK”ˆ’¢“'nF“lƒ±}fƒ¯ˆI&»&¡F™I™’˜1¿
"­¡FŸ $6.r¿¢Ör’ƒ¯‚fn<‚tŠ1“ ¿ ¬ Š1ln<‚K¿ È KÊÑÌ SÈ#ÉfÎ#ÊÑÌ ŠÉ Î#Ê tËɄË<Ì1Í[Ê Õ
Ë ŠÍMÊ¶ÈŠÌ 9¿)wz‘l‚fƒ¯l£In<‚„· ÕnF‚„º†Š£¢»
5;nF‚„º­ƒ¯?»1{¸nK”ˆI1“'nK“|ƒ±}fƒ¯ˆI&»&¡K™’™ Ô ¿
"­¡’¡ $ ¿ ƒ;’1“ 4;¿Õ)€1{„À’n<yI¿›ijln’“|®#’I} Š£InF{ˆ’°4°²ˆ’‚f¼;’‚f“›}fjlƒ¯1ÀHƒ­l£›ƒ¯ £’n<1n<‚ #}„ƒ­l£O‚fˆzˆŠ}fnF“=Š¢“
°²‚fn<n'}„‚fn<nF{F¿ 3à #Í Ì1Ì lÈ Ê ÈŠÌ Ê É„Ë<ÍÃË HÈ#É9ÊÑÍ Ð¿
jH}„}„‘ # #¼)¼)¼¿ }fjln<ˆI‚„yI¿ ”F{„”’¿ q rƒ†”Š¿ Ž °²‚„€¢{¸ÀIn<y ˆI€jH}„}„‘ # #¼)¼)¼¿ }fjln<ˆI‚„yI¿ ”F{„”’¿ q rƒ†”Š¿ Ž ”ˆI{ z¿
"­¡Kž $ ¿XsO€l‚„€¢l¿4tSˆ’‚ “|n<‚r”<ˆ’1“lƒ±}fƒ¯ˆI1{°²ˆ’‚t‘1Š‚„}„ƒ¯}„ƒ­ˆ’1nF“{¸yzmo‘lº­nF”}„ƒ†”mon}„j1ˆz“1{<¿ oË<É 0
Ì ­»¢˜Š ¢µ Ô !¹ # žIžŠŸŠ |žIž|¡’¡I»l¡F™I™HÒ|¿
"Ô $ ¿ ²î+õ Ñö
Unité de recherche INRIA Rocquencourt
Domaine de Voluceau - Rocquencourt - BP 105 - 78153 Le Chesnay Cedex (France)
Unité de recherche INRIA Futurs : Parc Club Orsay Université - ZAC des Vignes
4, rue Jacques Monod - 91893 ORSAY Cedex (France)
Unité de recherche INRIA Lorraine : LORIA, Technopôle de Nancy-Brabois - Campus scientifique
615, rue du Jardin Botanique - BP 101 - 54602 Villers-lès-Nancy Cedex (France)
Unité de recherche INRIA Rennes : IRISA, Campus universitaire de Beaulieu - 35042 Rennes Cedex (France)
Unité de recherche INRIA Rhône-Alpes : 655, avenue de l’Europe - 38334 Montbonnot Saint-Ismier (France)
Unité de recherche INRIA Sophia Antipolis : 2004, route des Lucioles - BP 93 - 06902 Sophia Antipolis Cedex (France)
Éditeur
INRIA - Domaine de Voluceau - Rocquencourt, BP 105 - 78153 Le Chesnay Cedex (France)
ISSN 0249-6399