最大流dinic算法,套模板1A
代码:
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#define maxn 105
using namespace std;
int map[maxn][maxn]; //邻接矩阵存储图
int dis[maxn]; //距离源点的距离
int q[maxn*maxn],f,r; //bfs队列
int n,m,ans;
int bfs()
{
memset(dis,-1,sizeof(dis)); //初始化为-1
dis[1] = 0;
f = r = 0;
q[r++] = 1;
while(f < r)
{
int j = q[f++];
for(int i = 1; i <= n; i++)
if(dis[i] < 0&&map[j][i] > 0) //如果i没有没有被别的点通到,并且j可以通到i,那么i的步数就是j的步数加一
{
dis[i] = dis[j] + 1;
q[r++] = i;
}
}
if(dis[n] > 0) return 1; //如果能够走到n
else return 0;
}
//Find代表一次增广,函数返回本次增广的流量,返回0表示无法增广
int find(int x,int low){ //
int i,a = 0;
if(x == n) return low;
for(i = 1; i <= n; i++)
if(map[x][i] > 0&&dis[i] == dis[x] + 1&&(a = find(i,min(low,map[x][i])))) //从i往后找最小的流量,就是一条增广路,并且此次流量不等于0
{
map[x][i] -= a;
map[i][x] += a; //回退边
return a;
}
return 0;
}
int main(){
int t,ca = 1;
scanf("%d",&t);
while(t--)
{
int a,b,val,tans;
scanf("%d%d",&n,&m);
memset(map,0,sizeof(map));
for(int i = 0; i < m; i++)
{
scanf("%d%d%d",&a,&b,&val);
map[a][b] += val; //可能会出现重边
}
ans = 0;
while(bfs()) //只要bfs能够进入终点
{
while(tans = find(1,0x7fffffff))
ans += tans;
}
printf("Case %d: %d\n",ca++,ans);
}
}