poj 3436 最大流

/*
题意：每一台computer都是有p部分组成，每一台机器都会对每台computer新加一些部分，删除一些部分，每一台机器每小时处理的computer的数量是固定的；给出的参数中，p个S表示对输入的p个部分的零件的要求：0表示必须没有，1表示必须有，2表示无所谓，p个D表示经过该机器后输出的p个部分零件的结果：0表示没有，1表示有；Q表示单位时间可处理的computer数量，最终求怎样安排这些机器的处理过程使得输出的效率最大；

题解：最大流
构图：将n台机器拆分成2*n个点，i->i+n表示当前机器i能处理的computer台数，然后将可以互相连接(即输入的p部分与输出的p部分需要相符合)的点连接，此处我用i+n连接j，表示机器i输出的computer可以由机器j处理，然后再加源点0，它的p部分的输出均为0，与1～n判断能否连接后再相连，权值为无限大，汇点的p部分输入均为1，n+1～2*n与汇点2*n+1判断后再相连，边权值为无限大；然后求源点到汇点的最大流即可。
*/
#include <cstdio>
#include <cstring>

#define MAXN 205
#define inf 1000000000

int max_flow(int n,int mat[][MAXN],int source,int sink,int flow[][MAXN]){
    int pre[MAXN],que[MAXN],d[MAXN],p,q,t,i,j;
    if (source==sink) return inf;
    for (i=0;i<n;i++)
        for (j=0;j<n;flow[i][j++]=0);
    for (;;){
        for (i=0;i<n;pre[i++]=0);
        pre[t=source]=source+1,d[t]=inf;
        for (p=q=0;p<=q&&!pre[sink];t=que[p++])
            for (i=0;i<n;i++)
                if (!pre[i] && (j=mat[t][i]-flow[t][i]))
                    pre[que[q++]=i]=t+1,d[i]=d[t]<j?d[t]:j;
                else if (!pre[i]&&(j=flow[i][t]))
                    pre[que[q++]=i]=-t-1,d[i]=d[t]<j?d[t]:j;
        if (!pre[sink]) break;
        for (i=sink;i!=source;)
            if (pre[i]>0)
                flow[pre[i]-1][i]+=d[sink],i=pre[i]-1;
            else
                flow[i][-pre[i]-1]-=d[sink],i=-pre[i]-1;
    }
    for (j=i=0;i<n; j+=flow[source][i++]);
    return j;
}

struct node
{
    int q;
    int s[11],d[11];
}nodes[55];

int p;

bool match(node a, node b)
{
    for(int i=0; i<p; i++)
        if (a.d[i] != b.s[i] && b.s[i] != 2)
            return false;
    return true;
}

int main(void)
{
    int n;
    int mat[MAXN][MAXN];
    int flow[MAXN][MAXN];
    while (~scanf("%d%d",&p,&n))
    {
        memset(mat,0,sizeof(mat));
        for(int i=1; i<=n; i++)
        {
            scanf("%d",&nodes[i].q);
            mat[i][i+n] = nodes[i].q; // 拆点
            for(int j=0; j<p; j++)
                scanf("%d",&nodes[i].s[j]);
            for(int j=0; j<p; j++)
                scanf("%d",&nodes[i].d[j]);
        }
        
        for(int i=0; i<p; i++)
        {
            nodes[0].d[i] = 0;
            nodes[n+1].s[i] = 1;
        }

        for(int i=1; i<=n; i++)
        {
            if (match(nodes[0],nodes[i]))
                mat[0][i] = 10005;
            if (match(nodes[i],nodes[n+1]))
                mat[i+n][2*n+1] = 10005;
        }
        for(int i=1; i<=n; i++)
            for(int j=1; j<=n; j++)
                if (i != j && match(nodes[i],nodes[j]))
                    mat[i+n][j] = nodes[i].q;

        printf("%d ",max_flow(2*n+2,mat,0,2*n+1,flow));

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