poj3436ACM Computer Factory（最大流）

http://poj.org/problem?id=3436

刚开始做的没拆点 BFS的时候保存了下路径 恰好样例都过了 接着就是各种WA

搜了搜解题报告 是将一个电脑拆成2个点 这两点之间限制容量 电脑与电脑之间容量为INF 设置一个源点和汇点  源点全是0 汇点全是1

不用保存路径 在最后 比较一下现在的容量和原始容量 差值就为W

#include <iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
#define INF 0x3f3f3f
int n,m;
struct node
{
    int s[15],d[15];
    int q;
}ma[110];
int path[110],flow[110],map[110][110],st,en,tmap[110][110];
int bfs()
{
    int i,j,k;
    memset(path,-1,sizeof(path));
    for(i = 0 ; i <= en ; i++)
    flow[i] = INF;
    queue<int>q;
    q.push(0);
    while(!q.empty())
    {
        int tk = q.front();
        q.pop();
        if(tk==en)
        break;
        for(i = 1 ; i <= en ; i++)
        {
            if(path[i]==-1&&map[tk][i])
            {
                path[i] = tk;
                flow[i] = min(flow[tk],map[tk][i]);
                q.push(i);
            }
        }
    }
    if(path[en]==-1)
    return -1;
    return flow[en];
}
int karp()
{
    int i,j,k,now,pre,sum=0;
    while((k=bfs())!=-1)
    {
        sum+=k;
        now = en;
        while(now!=st)
        {
            pre = path[now];
            map[pre][now]-=k;
            map[now][pre]+=k;
            now = pre;
        }
    }
    return sum;
}
int main()
{
    int i,j,k,g;
    while(cin>>n>>m)
    {
        memset(map,0,sizeof(map));
        memset(tmap,0,sizeof(tmap));
        for(i = 1; i <= m ; i++)
        {
            cin>>ma[i].q;
            for(j = 1 ; j <= n ; j++)
            cin>>ma[i].s[j];
            for(j = 1; j <= n; j++)
            cin>>ma[i].d[j];
            map[i][i+m] = ma[i].q;
            tmap[i][i+m] = ma[i].q;
        }
        for(i = 1; i <= m ; i++)
        {
            int f = 1;
            for(j = 1; j <= n ; j++)
            if(ma[i].s[j]==1)
            {
                f = 0;
                break;
            }
            if(f)
            {
                map[0][i] = INF;
                tmap[0][i] = INF;
            }
        }
        for(i = 1; i <= m ; i++)
        {
            int f = 1;
            for(j = 1; j <= n ; j++)
            if(ma[i].d[j]!=1)
            {
                f = 0;
                break;
            }
            if(f)
            {
                map[i+m][2*m+1] = INF;
                tmap[i+m][2*m+1] = INF;
            }
        }
        for(i = 1; i <= m ; i++)
        for(j = 1; j <= m ; j++)
        {
            int flag = 1;
            if(i==j)
            continue;
            for(g = 1; g <= n ; g++)
            {
                if(ma[j].s[g]==2)
                continue;
                if(ma[i].d[g]!=ma[j].s[g])
                {
                    flag = 0;
                    break;
                }
            }
            if(flag)
            {
                map[i+m][j] = INF;
                tmap[i+m][j] = map[i+m][j];
            }
        }
        en = 2*m+1;st = 0;
        int num = 0;
        int ss  = karp();
        for(i = m+1; i < en ;i++)
        for(j = 1 ; j <= m ; j++)
        if(i!=j&&map[i][j]<tmap[i][j])
        num++;
        cout<<ss<<" "<<num<<endl;
        for(i = m+1; i < en ;i++)
        for(j = 1 ; j <= m ; j++)
        if(i!=j&&map[i][j]<tmap[i][j])
        {
            cout<<i-m<<" "<<j<<" "<<tmap[i][j]-map[i][j]<<endl;
        }
    }
    return 0;
}