Задача:
Найдите величину максимального потока в сети.
Входные данные.
Первая строка входного файла содержит целое число N - количество вершин сети. Вторая строка входного файла содержит целое число M - количество дуг сети. В каждой из следующих M строк содержится ровно три числа a, b, c. Эти числа описывают дугу, идущую из вершины с номером a в вершину с номером b и имеющую пропускную способность c. Последние две строки содержат целые числа s и t - номера вершин, являющихся истоком и стоком, соответственно. Вершины нумеруются последовательными натуральными числами от 1 до N. В сети между фиксированной парой вершин в фиксированном направлении может быть не более одной дуги.
Выходные данные.
Единственная строка выходного файла должна содержать одно целое число, равное величине максимального потока в сети, описанной во входном файле.
Пример.
Идея:
Повторяем, пока не найдём максимальный поток
1. Ищем путь от вершины 1 до вершины n и делаем соответствующие заполнения
2. В специальную переменную суммируем результат потока на каждом шаге
3. Заполняем все «трубы» пути на i-том шаге максимальным потоком на этом шаге
Задача о максимальном потоке заключается в нахождении потока такого, что величина потока максимальна.
Код:
type zap=record pr,dp:longint;end; //для хранения информации о потоке используем записи var c,f:array[0..900,0..900]of longint; q:array[0..900]of longint; p:array[0..900]of zap; i,s,t,v,n,m,a,b,d,pot,un,uk,x,y:longint; ud:boolean; Function min(a,b:longint):longint; //ну тут всё понятно... begin if a<b then min:=a else min:=b; end; Procedure PUD; //процедура Поиск Удлиняющей Дуги begin un:=1;uk:=2;q[1]:=s; p[s].pr:=s; p[s].dp:=1 shl 30; while (un<>uk)and(p[t].dp=0) do //повторяем, пока есть куда лить и поток не дошёл до конца begin v:=q[un];un:=un+1; //берём из очереди вершину for i:=1 to n do if ((c[v,i]>0)or(c[i,v]>0))and(p[i].dp=0) //если из i в v (или наоборот) есть путь и на этом шагу воды не добавляли,... then if (f[v,i]<c[v,i]) //то если сюда можно протечь,... then begin q[uk]:=i; uk:=uk+1; //то добавляем вершину в очередь p[i].pr:=v; // предки - это святое p[i].dp:=min(p[v].dp,c[v,i]-f[v,i]); //выбираем количество воды, которое туда потечёт (сколько текло в v или сколько можно впихнуть) end else if f[i,v]>0 then begin //иначе если течёт по обратному ребру,... q[uk]:=i; uk:=uk+1; //добавляем вершину в очередь p[i].pr:=-v; //как я уже говорил, предки - это святое, но чтобы знать, что мы получили жидкость из обратного ребра, делаем ему отрицательный знак p[i].dp:=min(p[v].dp,f[i,v]); //выбираем количество воды для течения (сколько текло в v или сколько льется) end; end; if p[t].dp=0 then ud:=false; //если ничего не дотекло, значит мы нашли максимальный поток end; Procedure VUD(v:longint); //процедура Восстановление Удлиняющей Дуги begin if abs(p[v].pr)=v then exit; //если каким-либо образом мы попали из самого себя в себя, то мы в начале if p[v].pr>0 then f[p[v].pr,v]:=f[p[v].pr,v]+p[t].dp //если есть прямой предок, то добавляем к потоку в трубе предок v и v столько, сколько долилось в последнюю вершину else if p[v].pr<0 then f[v,-p[v].pr]:=f[v,-p[v].pr]-p[t].dp; //иначе, если есть предок по обратной дуге, то уменьшаем количество воды в трубе v и -предок v на столько, сколько дотекло до последней вершины VUD(abs(p[v].pr)); //восстановление дуги от предка вершины и соответствующие заполнения end; Begin read(n,m); //ввод данных в соответствующие для них массивы for i:=1 to m do begin read(x,y); read(c[x,y]); end; read(s,t); ud:=true; //переменная, которая определяет способность пролить воду из первой в последнюю вершину while ud do //повторяем алгоритм, пока не найдём максимальный поток begin for i:=1 to n do //так как мы запускаем процедуры по несколько раз, обнуляем массивы begin p[i].pr:=0; p[i].dp:=0; end; PUD; if ud then begin pot:=pot+p[t].dp; //переменная, хранящая в себе максимальный поток VUD(t); //восстанавливаем и заполняем с конца end; end; writeln(pot); //выводим величину максимального потока End.