רשת זרימה

בתורת הגרפים, רשת זרימה היא גרף מכוון, עם שני צמתים מיוחדים: צומת מקור וצומת בור, בו לכל קשת מוגדרת כמות הזרימה המקסימלית היכולה לעבור בה (פרוט בהמשך). ניתן להשתמש ברשתות זרימה כדי למדל זרימה של נוזל בצינורות, מעבר של מידע ברשתות תקשורת, מעבר של תנועה בכביש, זרם ברשתות חשמל, ועוד. לרשתות זרימה ישנם גם שימושים בפתרון בעיות נוספות כמו בעיית שידוך מקסימלי בגרפים לא מכוונים.

כל רשת זרימה מאופיינת על ידי גרף מכוון שמכיל שני צמתים מיוחדים - האחד משמש כמקור - המקום שממנו נובעת הזרימה, והשני משמש כבור - המקום שאליו מתנקזת הזרימה. כמו כן לכל קשת בגרף ישנו קיבול, שמתאר את כמות הזרימה המקסימלית שיכולה לעבור בה בפרק זמן נתון.

השאלה העיקרית שבה עוסקים כאשר חוקרים רשת זרימה היא מה הכמות הגדולה ביותר של זרימה שניתן להעביר מן המקור ועד לבור בפרק זמן נתון, ומה הדרך שבה ניתן לעשות זאת. בעיה זו מכונה "בעיית זרימת מקסימום". דרך מקובלת לפתרון שאלה זו היא הגדרה של פונקציה שמייצגת זרימה, ואז הפעלת אלגוריתם שבודק את הדרכים שבהן ניתן לשפר פונקציה זאת, עד אשר היא הופכת לאופטימלית.

==הגדרה פורמלית==
===רשת זרימה===
רשת זרימה מורכבת מגרף <math>\ G=(V,E)</math>, משני צמתים מיוחדים - המקור (Source) <math>\ s\isin V</math> והבור (Sink) <math>\ t\isin V</math>, ומפונקציית קיבול <math>\ c:V\times V\to\mathbb{R}^{+}\cup\left\{0,\infty\right\}</math> שמתאימה לכל זוג מסודר של צמתים את כמות הזרימה המקסימלית שיכולה לעבור מהראשון שבהם אל השני (כמות זו היא מספר אי שלילי, ויכולה להיות אינסוף). אם לא קיימת קשת מ-<math>\ u</math> אל <math>\ v</math>, מניחים כי <math>\ c(u,v)=0</math>.

ניתן להניח כי כל צומת ברשת נמצא על מסלול מהמקור אל הבור (אחרת ניתן להסירו בלי לשנות את קבוצת פונקציות הזרימה האפשריות). כמו כן, ההנחה כי יש רק מקור אחד ובור אחד אינה מגבילה - ניתן להפוך כל רשת זרימה עם מספר שרירותי של מקורות לרשת עם מקור יחיד על ידי הוספת צומת נוסף, "מקור-על", שיוצאת ממנו קשת עם קיבולת אינסופית אל כל אחד מהמקורות ברשת המקורית. בדרך דומה ניתן להוסיף גם "בור-על".

===זרימה===

בהינתן רשת זרימה, רוצים לתאר זרימה ברשת בצורה פורמלית. זרימה באה למדל את מעבר החומר מהמקור אל הבור, דרך שאר צמתי הגרף. התכונה המאפיינת הבסיסית של זרימה היא שכמות הזרימה שנכנסת לצומת דרך הקשתות שנכנסות אליו חייבת להיות שווה לכמות הזרימה היוצאת ממנו דרך הקשתות היוצאות ממנו, כששני היוצאים מן הכלל הם המקור (שניתן להניח שרק יוצאת ממנו זרימה) והבור (שניתן להניח שרק נכנסת אליו זרימה).

בצורה פורמלית, ניתן לתאר זרימה על ידי פונקציית זרימה <math>\ f:V\times V\to\mathbb{R}</math>. פונקציה זו מתאימה לכל זוג צמתים <math>\ (u,v)</math> את כמות הזרימה מ-<math>\ u</math> אל <math>\ v</math> (לכיוון יש חשיבות), כך שמתקיימים התנאים הבאים:
#אנטי-סימטריה: <math>\ f(u,v)=-f(v,u)</math>
#שמירה על אילוץ הקיבול: <math>f(u,v)\le c(u,v)</math>
#שימור הזרימה: לכל <math>\ u\ne s,t</math> מתקיים: <math>\ \sum_{v\in V}f(u,v)=0</math>
דוגמה לרשת זרימה. עבור כל קשת מתוארים הקיבול/גודל הזרם שמועבר.
תכונת האנטי סימטריה נועדה להקל על הטיפול המתמטי ברשת הזרימה. מכיוון שזרימה שעוברת מצומת <math>\ u</math> אל הצומת <math>\ v</math> מתוארת על ידי מספר חיובי (כמות הזרימה שעברה), ניתן לחשוב על כך כאילו הזרימה עברה מהצומת <math>\ v</math> אל <math>\ u</math> (כלומר, בכיוון השני) אך בכמות שלילית. הדבר מקל, לדוגמה, על כתיבת האילוץ השלישי, ומאפשר לדבר על הזרימה בין שני צמתים בשני הכיוונים גם יחד, ולא רק בכיוון אחד, גם אם קיימת קשת ביניהם רק בכיוון אחד.

שימור הזרימה פירושו שכל הזרימה שנכנסה לצומת גם יצאה ממנו. מכיוון שזרימה שנכנסת לצומת מיוצגת על ידי מספר חיובי, וזרימה שיוצאת מהצומת מיוצגת על ידי מספר שלילי, הרי שמצב זה מתואר על ידי כך שסכום הזרימות שווה לאפס.

לכל זרימה מגדירים את ערך (או כמות או שטף) הזרימה בתור <math>\ |f|=\sum_{v\isin V}f(s,v)</math>, כלומר כמות הזרימה נטו שיוצאת מן המקור. ניתן להראות בהתבסס על תכונות הזרימה כי כמות זו שווה לכמות הזרימה נטו שנכנסת אל הבור. על כן, ניתן לחשוב על ערך הזרימה בתור כמות הזרימה שעוברת ברשת הזרימה ביחידת זמן.

יש לשים לב כי פונקציית הזרימה מגדירה את כמות הזרימה נטו שעוברת בין שני הצמתים. לדוגמה, זרימה בגודל 5 שעוברת מצומת <math>\ u</math> אל צומת <math>\ v</math> אין פירושה בהכרח כי עוברות 5 יחידות של זרימה מ-<math>\ u</math> אל <math>\ v</math> - ייתכן למשל כי בפועל עוברות 12 יחידות מ-<math>\ u</math> אל <math>\ v</math>, ועוברות 7 יחידות מ-<math>\ v</math> אל <math>\ u</math>, כך שהשינוי נטו בכמות החומר בשני הצמתים מתבטא בכך שעוברות 5 יחידות מ-<math>\ u</math> אל <math>\ v</math>.

== שיטת פורד פולקרסון ==

שיטה מקובלת למציאת הזרימה האופטימלית (כלומר, בעלת הערך הגבוה ביותר) שניתן להעביר ברשת זרימה היא שיטת פורד פולקרסון. השיטה היא סכמה כללית, וקיימים מספר אלגוריתמים, בעלי סיבוכיות זמן שונה, המממשים אותה.

הרעיון הבסיסי שעומד מאחורי השיטה הוא שיפור איטרטיבי: מתחילים עם פונקציית זרימה שנותנת ערך 0 לכל קשת, ולאחר מכן מחפשים דרכים לשפר אותה על ידי בדיקת מסלולים מהמקור אל הבור. אם מתגלה מסלול שבו ניתן להעביר עוד זרם, הפונקציה מתוקנת בהתאם לכך, והתהליך חוזר על עצמו עד שמתקבלת הפונקציה האופטימלית.

==שיטת הרם-וקדם==
שיטת הרם-וקדם (Push–relabel) מייצגת גישה שונה לפתרון בעיית זרימת המקסימום, שאינה דומה לשיטת פורד-פולקרסון. ניתן לממש אותה באופן יעיל כך שיתקבל אלגוריתם שזמן הריצה שלו הוא <math>\ O(|V|^3)</math> - שיפור יחסי לזמני הריצה של מימושי שיטת פורד-פולקרסון עבור גרפים בעלי קשתות רבות.

ההבדל המרכזי בין שיטת הרם-וקדם ושיטת פורד-פולקרסון הוא ששיטת הרם-וקדם פועלת בצורה לוקלית ובכל איטרציה עוסקת רק בצומת בודד ובשכניו, בעוד ששיטת פורד-פולקרסון היא גלובלית ובכל איטרציה עוסקת בגרף כולו ומחפשת בו מסלולי שיפור. מכאן נובע היתרון היחסי של שיטת הרם-וקדם.

הרעיון הבסיסי שמנחה את השיטה הוא זה: בתחילה, מזרימים מן המקור את כמות החומר המקסימלית האפשרית - כלומר, מרווים את כל הקשתות היוצאות ממנו. התוצאה אינה פונקציות זרימה חוקית, שכן זרימה נכנסת לשכניו של המקור אך אינה יוצאת מהן, ולכן לא מתקיים אילוץ הזרימה - לצמתים הללו ישנה "זרימה עודפת". האלגוריתם עובר על צמתים בעלי זרימה עודפת ומעביר מהן חלק מהזרימה הלאה (זהו ה-"קדם" שבשם השיטה, ובאנגלית Push), עד אשר היא מגיעה אל הבור או חוזרת אל המקור. ניתן להראות כי כאשר לא יהיה באף צומת עודף זרימה, כלומר תתקבל זרימה חוקית, יתקיימו גם תנאי משפט זרימה-מקסימלית חתך-מינימלי, ולכן הזרימה שתתקבל היא אופטימלית.

כדי להבטיח את עצירת האלגוריתם במהירות, לא מתירים העברת זרימה מכל צומת לשכן שלו, אלא מגדירים כלל נוסף, המבוסס על תכונת "גובה" של הצמתים. בגרף בעל n צמתים נותנים תחילה לכל הצמתים פרט למקור את הגובה 0, ולמקור נותנים גובה n. לאחר ההזרמה הראשונית מהמקור לשכניו, מתקיים הכלל הבא: מותר לצומת להזרים רק לצומת שקטן ממש ממנו (ניתן לדרוש אף שהפרשי הגבהים ביניהם יהיו בדיוק 1). כדי לאפשר שימוש בכלל זה מגדירים פעולה נוספת שאותה האלגוריתם מסוגל לבצע - הרמת צומת כך שיהיה גבוה מכל שכניו הרלוונטיים (זוהי פעולת ה"הרם" שבשם השיטה, ובאנגלית Lift או Relabel). הגובה נבחר בתור הגובה המינימלי הדרוש שעבורו הצומת יהיה גבוה מכל השכנים שעוד ניתן להזרים אליהם (כלומר, שהקשת אליהם אינה רוויה).

במימוש הנאיבי של שיטת הרם-וקדם, האלגוריתם הוא כדלהלן:
#הזרם מן המקור לשכניו את הכמות המקסימלית האפשרית שמתאפשרת על ידי קיבולי הקשתות המחברות ביניהם.
#קבע את גובהו של המקור ל-n ואת גובה שאר הצמתים ל-0.
#כל עוד קיימת זרימה עודפת עבור אחד הצמתים, בצע אחת משתי הפעולות הבאות עבורו:
##פעולת "קדם" לאחד משכניו, אם ניתן (הקשת בין שני הצמתים לא רוויה, והצומת המזרים גבוה משכנו).
##פעולת "הרם", אם ניתן (הצומת קטן או שווה בגובהו מכל השכנים שעוד ניתן להזרים אליהם).

מימוש זה מבטיח זמן ריצה של <math>\ O(|V|^2\cdot|E|)</math>. כדי להבטיח זמן ריצה של <math>\ O(|V|^3)</math> משכללים את דרך הבחירה של הצומת שעליו יבוצעו פעולות ההרם-וקדם, כך שהאלגוריתם מתמקד בכל פעם בצומת בודד ומרוקן אותו לחלוטין.

==יישומים==
===גרפים דו צדדיים===
דוגמה לרשת זרימה שמתקבלת לאחר הבניה
ניתן להשתמש ברשתות זרימה כדי למצוא שידוך מקסימום בגרף דו צדדי בזמן הדרוש לריצת אלגוריתם למציאת זרימת מקסימום.

בהינתן גרף דו צדדי המורכב משתי קבוצות צמתים <math>\ A,B</math>, כך שכל קשת היא בין צומת השייך ל-<math>\ A</math> וצומת השייך ל-<math>\ B</math>, בונים רשת זרימה כדלהלן: ראשית, מוסיפים צמתים חדשים שישמשו כמקור וכבור: <math>\ s,t</math>. כעת, מחברים בקשתות את <math>\ s</math> לכל צומתי <math>\ A</math>, ואת <math>\ t</math> לכל צומתי <math>\ B</math>. נותנים לכל הקשתות קיבולת 1. לא קשה להוכיח כי זרימת מקסימום ברשת זו משרה שידוך מקסימום בגרף המקורי, כאשר הקשתות שמשתתפות בשידוך הן אלו בהן עובר זרם.

==לקריאה נוספת==
* תומאס ה. קורמן, צ'ארלס א. לייזרסון, רונאלד ל. ריבסט, מבוא לאלגוריתמים, הוצאת MIT (פרק 26), תורגם לעברית על ידי האוניברסיטה הפתוחה
* שמעון אבן (1979), Graph Algorithms, הוצאת Computer Science Press

==קישורים חיצוניים==

* 

קטגוריה:רשתות זרימה
קטגוריה:תורת הגרפים
קטגוריה:אלגוריתמים בתורת הגרפים