POJ 3436 ACM Computer Factory（拆点+前向星dinic）||（拆点+邻接矩阵dinic）||(不拆点+dinic))

题目大意：就是第一行两个数p,n分别代表电脑零件的个数，和有多少机器。接下来的n行每行2*p+1个数，第一个数为机器工作的效率，其余代表机器加工前和后电脑的变化，求机器一起工作的最大效率，并输出路径。

思路：由于源点的流入和普通链接点的流入，可能会大于当前点的容量，因此拆点.那么，怎么再判断有多少个机器用到了呢，这就要原来机器的工作效率，和剩余网络，如果原来的工作效率大于剩余网络那么一定该机器被用到了既有流。（注意建图时后来加源和汇,且机器和机器相连的时候应该注意到分开，因此就要相应的加上m值，但是应该加到i,j谁的值上，答案是i上因为如果加到j上相应的遍历的时候要减去m(Q[j].to-m)因为有逆向建图，所以在以a为下一个链接点时可能出现负值在vv[][]数组中 。）

法一：前向星

#include<iostream>
#include<cstdio>
#include<cstring>
#include<map>
#include<queue>
#include<algorithm>
#define MAX 1000
#define inf 0x3f3f3f3f
using namespace std;
struct node
{
    int w,in[52],out[52];
}q[100];

struct no
{
    int to,w,next;
}Q[100000];

int s,e,cnt,l[51100],qu[50000],head[50000];
int vv[1100][1100];
void bu(int a,int b,int c){
    Q[cnt].to=b;
    Q[cnt].w=c;
    Q[cnt].next=head[a];
    head[a]=cnt++;

    Q[cnt].to=a;
    Q[cnt].w=0;
    Q[cnt].next=head[b];
    head[b]=cnt++;
}
int bfs(){
    memset(l,-1,sizeof(l));
    l[s]=0;
    queue<int>qu;
    while(!qu.empty())
        qu.pop();
    qu.push(s);
    while(!qu.empty()){
        int v=qu.front();
        qu.pop();
        for(int i=head[v];~i;i=Q[i].next){
            if( l[Q[i].to ]==-1&&Q[i].w>0 ){
                l[Q[i].to ]=l[v]+1;
                qu.push(Q[i].to);
            }
        }
    }
    return l[e]!=-1;
}

int dfs(int x,int f){
    if(x==e||!f) return f;
    int tmp=0,a;
    for(int i=head[x];~i;i=Q[i].next){
        if(l[x]+1==l[Q[i].to ]&&Q[i].w>0 ){
            a=dfs(Q[i].to,min(Q[i].w,f-tmp));
            if(a>0){
                Q[i].w-=a;
                Q[i^1].w+=a;
                tmp+=a;
                if(tmp==inf)break;
            }
        }
    }
    if(!tmp) l[x]=-1;
    return tmp;
}

int main(){
    int n,m,j,i,k;
    while(~scanf("%d%d",&n,&m)){
        cnt=0;
        memset(vv,0,sizeof(vv));
        memset(head,-1,sizeof(head));
        s=0,e=2*m+1;
        for(i=1;i<=m;i++){
            int bj=false;
            scanf("%d",&q[i].w);
            bu(i,i+m,q[i].w);///***一般源和汇都是后来的时候才链接
            for(j=0;j<n;j++){
                scanf("%d",&q[i].in[j]);
                 if(q[i].in[j]==1){
                    bj=true;
                }
            }
            if(!bj)
                bu(0,i,q[i].w);
            bj=false;
            for(j=0;j<n;j++){
                scanf("%d",&q[i].out[j]);
                if(q[i].out[j]!=1){
                    bj=true;
                }
            }
            if(!bj)
                bu(i+m,2*m+1,q[i].w);
        }
        for(i=1;i<=m;i++){
            for(j=1;j<=m;j++){
                if(i==j)continue;
                bool vis=false;
                for(k=0;k<n;k++){
                    if( (q[i].out[k]!=q[j].in[k])&&q[j].in[k]!=2){
                        vis=true;break;
                    }
                }
                if(!vis){
                    bu( i+m,j,q[i].w);///****
                    vv[m+i][j]=1;
                }
            }
        }

        int ans=0;
        while(bfs()){
            ans+=dfs(s,inf);
        }
        int sum=0;
        for(i=1;i<=m;i++){///*********
            for(j=head[i+m];j!=-1;j=Q[j].next){
                if( vv[i+m][Q[j].to ]&&Q[j].w<q[i].w ){
                    sum++;
                }
            }
        }

        printf("%d %d\n",ans,sum);
        for(i=1;i<=m;i++){
            for(j=head[i+m];~j;j=Q[j].next){
                if(vv[i+m][Q[j].to]&&Q[j].w<q[i].w){
                    printf("%d %d %d\n",i,Q[j].to,q[i].w-Q[j].w);
                }
            }
        }
    }
}

法二：邻接矩阵+拆点dinic

#include<iostream>
#include<cstdio>
#include<cstring>
#include<map>
#include<queue>
#include<algorithm>
#define MAX 1000
#define inf 0x3f3f3f3f
using namespace std;
struct node
{
    int w,in[52],out[52];
}q[100];

int Map[1100][1100],s,e,l[1100];
int flow[1100][1100];

int bfs(){
    memset(l,-1,sizeof(l));
    l[s]=0;
    queue<int>Q;
    while(!Q.empty())
        Q.pop();
    Q.push(s);

    while(!Q.empty()){
        int v=Q.front();
        Q.pop();
        for(int i=0;i<=e;i++){
            if(Map[v][i]&&l[i]==-1){
                l[i]=l[v]+1;
                Q.push(i);
            }
        }
    }
    return l[e]!=-1;
}
int dfs(int x,int f){
    if(x==e||!f){
        return f;
    }
    int tmp=0,a;
    for(int i=0;i<=e;i++){
        if(l[i]==l[x]+1&&Map[x][i]){
            a=dfs(i,min(Map[x][i],f-tmp));
            if(a>0){
                Map[x][i]-=a;
                Map[i][x]+=a;
                flow[x][i]+=a;
                flow[i][x]-=a;
                tmp+=a;
                if(tmp==f)break;
            }
        }
    }
    if(!tmp)l[x]=-1;
    return tmp;
}

int main()
{
    int n,m,i,j,k;
    while(~scanf("%d%d",&n,&m)){
        memset(flow,0,sizeof(flow));
        s=0,e=2*m+1;
        memset(Map,0,sizeof(Map));

        for(i=1;i<=m;i++){
            scanf("%d",&q[i].w);
            for(j=0;j<n;j++){
                scanf("%d",&q[i].in[j]);
            }
            for(j=0;j<n;j++){
                scanf("%d",&q[i].out[j]);
            }
        }

        bool bj;
        for(i=1;i<=m;i++){
            for(j=1;j<=m;j++){
                if(i!=j){
                    bj=false;
                    for(k=0;k<n;k++){
                        if( (q[i].out[k]!=q[j].in[k])&&q[j].in[k]!=2){
                            bj=true;break;
                        }
                    }
                    if(!bj){
                        Map[i+m][j]=q[i].w;///**
                    }
                }
            }
        }
        for(i=1;i<=m;i++){
            Map[i][i+m]=q[i].w;
            bj=false;
            for(j=0;j<n;j++){
                if(q[i].in[j]==1){
                    bj=true;break;
                }
            }
            if(!bj){
                Map[0][i]=q[i].w;
            }
            bj=false;

            for(j=0;j<n;j++){
                if(q[i].out[j]!=1){
                    bj=true;break;
                }
            }
            if(!bj)
                Map[i+m][2*m+1]=q[i].w;
        }
        int ans=0;
        while(bfs()){
            ans+=dfs(0,inf);
        }
        printf("%d ",ans);
        int sum=0;

        for(i=1;i<=e;i++){
            for(j=1;j<=i;j++){
                if(flow[i][j]>0)
                    sum++;
            }
        }
        printf("%d\n",sum);
        for(i=1;i<=e;i++){
            for(j=1;j<=i;j++){
                if(flow[i][j]>0)
                    printf("%d %d %d\n",i-m,j,flow[i][j]);
            }
        }
    }
    return 0;
}

法三：邻接矩阵不拆点（可能题目水了，不拆点也可以过，直接将上边代码的源点汇点和相邻的中间点去掉相应的m即可。）