最大流dinic模板

1. 运输问题1
   ★★☆ 输入文件：maxflowa.in 输出文件：maxflowa.out 简单对比
   时间限制：1 s 内存限制：128 MB
   【问题描述】
   一个工厂每天生产若干商品，需运输到销售部门进行销售。从产地到销地要经过某些城镇，有不同的路线可以行走，每条两城镇间的公路都有一定的流量限制。请你计算，在不考虑其它车辆使用公路的前提下，如何充分利用所有的公路，使产地运输到销地的商品最多，最多能运输多少商品。
   【输入格式】
   输入文件有若干行
   第一行，一个整数n，表示共有n个城市(2<=n<=100),产地是1号城市，销地是n号城市。
   下面有n行,每行有n个数字。第p行第q列的数字表示城镇p与城镇q之间有无公路连接。数字为0表示无，大于0表示有公路，且该数字表示该公路流量。
   【输出格式】
   输出文件有一行
   第一行，1个整数max，表示最大流量为max。
   【输入输出样例】
   输入文件名： maxflowa.in
   6
   0 4 8 0 0 0
   0 0 4 4 1 0
   0 0 0 2 2 0
   0 0 0 0 0 7
   0 0 0 6 0 9
   0 0 0 0 0 0
   输出文件名：maxflowa.out
   8

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int maxn=1010,Inf=0x7fffffff/3;
int n,s,x,tot=1,ans=0;
struct edge{
    int to,next,w;
}G[maxn*10];
int h[maxn],q[maxn*10],cur[maxn],vis[maxn];

void add(int x,int y,int z){
    ++tot; G[tot].to=y;
    G[tot].w=z;
    G[tot].next=h[x]; h[x]=tot;
}

bool bfs(){
    memset(vis,-1,sizeof(vis));
    int head=1;
    int t=1,now;
    q[1]=1; vis[1]=0;
    while (head<=t){
        now=q[head];head++;
        for (int i=h[now];i;i=G[i].next)
           if (vis[G[i].to]==-1&&G[i].w>0){
                vis[G[i].to]=vis[now]+1;
                q[++t]=G[i].to;
           }
    }
    return vis[n]!=-1;
}

int dfs(int x,int f){
    if (x==n) return f;
    int w,used=0;
    for (int i=cur[x];i;i=G[i].next)
       if (vis[G[i].to]==vis[x]+1){
            w=f-used;
            w=dfs(G[i].to,min(w,G[i].w));
            G[i].w-=w; G[i^1].w+=w;
            if (G[i].w>0) cur[x]=i;
            used+=w;
            if (used==f) return f;
       }
    if (!used) vis[x]=-1;
    return used;
}

void dinic(){
    ans=0;
    while (bfs()){
        for (int i=1;i<=n;++i) cur[i]=h[i];
        ans+=dfs(1,Inf);
    }return;
}

int main()
{
    freopen("maxflowa.in","r",stdin);
    freopen("maxflowa.out","w",stdout);
    scanf("%d",&n);
    for (int i=1;i<=n;++i)
       for (int j=1;j<=n;++j){
            int x;
            scanf("%d",&x);
            if (x>0){
               add(i,j,x);
               add(j,i,0);
            }
       }
    dinic();
    printf("%d",ans);
    return 0;
}