COJ1249(竞争性酶抑制剂和同工酶)

题目链接

由于之前没写过网络流方面的题目，所以第一次写的时候居然想到用动态规划去做（在没有环的情况下貌似可以，有环会RE）。现在看来，这题是个求最小割的题，因为最大流量最小割相等，所以也就是求最大流量。容量网络也很明显。我用的是标号法（Edmonds-Karp算法）。

#include <stdio.h>
#include <string.h>
#define N 152
#define INF 0x7fffffff
#define MIN(a,b) ((a)<(b)?(a):(b))
int c[N][N],flow[N][N],a[N],p[N],f,n,m;
int queue[N],front,rear;
void EK(int s,int t)
{
    int u,v;
    f=0;
    memset(flow,0,sizeof(flow));
    while(1)
    {
        memset(a,0,sizeof(a));
        front=rear=0;
        queue[rear]=s,rear=(rear+1)%N;
        a[s]=INF,p[s]=0;
        while(rear!=front)
        {
            u=queue[front],front=(front+1)%N;
            for(v=1;v<=n;v++)
            {
                if(!a[v]&&c[u][v]>flow[u][v])
                {
                    a[v]=MIN(a[u],c[u][v]-flow[u][v]);
                    p[v]=u;
                    queue[rear]=v,rear=(rear+1)%N;
                }
            }
        }
        if(!a[t])   break;
        for(u=t;u!=s;u=p[u])
        {
            flow[p[u]][u]+=a[t];
            flow[u][p[u]]-=a[t];
        }
        f+=a[t];
    }
}
int main()
{
    int i,x,y,z,s,t;
    while(~scanf("%d%d",&n,&m))
    {
        memset(c,0,sizeof(c));
        for(i=0;i<m;i++)
        {
            scanf("%d%d%d",&x,&y,&z);
            c[x][y]+=z;
        }
        scanf("%d%d",&s,&t);
        EK(s,t);
        printf("%d\n",f);
    }
    return 0;
}