A - ACM计算机工厂 POJ - 3436

就是EK建图

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<algorithm>
#include<queue>
using namespace std;
const int MAXN=110;
const int INF=0x3fffffff;
int g[MAXN][MAXN];//存边的容量，没有边的初始化为0
int path[MAXN],flow[MAXN],start,end;
int n;//点的个数，编号0-n.n包括了源点和汇点。

queue<int>q;
int bfs()
{
    int i,t;
    while(!q.empty())q.pop();//把清空队列
    memset(path,-1,sizeof(path));//每次搜索前都把路径初始化成-1
    path[start]=0;
    flow[start]=INF;//源点可以有无穷的流流进
    q.push(start);
    while(!q.empty())
    {
        t=q.front();
        q.pop();
        if(t==end)break;
        //枚举所有的点，如果点的编号起始点有变化可以改这里
        for(i=0;i<=n;i++)
        {
            if(i!=start&&path[i]==-1&&g[t][i])
            {
                flow[i]=flow[t]<g[t][i]?flow[t]:g[t][i];
                q.push(i);
                path[i]=t;
            }
        }
    }
    if(path[end]==-1)return -1;//即找不到汇点上去了。找不到增广路径了
    return flow[end];
}
int Edmonds_Karp()
{
    int max_flow=0;
    int step,now,pre;
    while((step=bfs())!=-1)
    {
        max_flow+=step;
        now=end;
        while(now!=start)
        {
            pre=path[now];
            g[pre][now]-=step;
            g[now][pre]+=step;
            now=pre;
        }
    }
    return max_flow;
}

int backup[MAXN][MAXN];//备份图
int in[MAXN][20];//输入信息
int Line[MAXN][4];

int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    int p;
    while(scanf("%d%d",&p,&n)!=EOF)
    {
        memset(g,0,sizeof(g));
        for(int i=1;i<=n;i++)
        {
            for(int j=0;j<2*p+1;j++)
              scanf("%d",&in[i][j]);
        }
        n++;
        start=0;//源点
        end=n;//汇点
        for(int i=1;i<=n-1;i++)
        {
            bool flag_s=true;
            bool flag_t=true;
            for(int j=0;j<p;j++)
            {
                if(in[i][j+1]==1)flag_s=false;//不能有1
                if(in[i][j+1+p]==0)flag_t=false;
            }
            if(flag_s)g[0][i]=in[i][0];
            if(flag_t)g[i][n]=in[i][0];
            for(int j=1;j<=n-1;j++)
               if(i!=j)
               {
                   bool flag=true;
                   for(int k=0;k<p;k++)
                     if((in[i][k+p+1]==0&&in[j][k+1]==1)||(in[i][k+p+1]==1&&in[j][k+1]==0))
                     {
                         flag=false;
                         break;
                     }
                   if(flag)g[i][j]=min(in[i][0],in[j][0]);
               }

        }
        memcpy(backup,g,sizeof(g));//先把图备份下来
        printf("%d ",Edmonds_Karp());
        int tol=0;
        for(int i=1;i<=n-1;i++)
          for(int j=1;j<=n-1;j++)
          {
              if(g[i][j]<backup[i][j])
              {
                  Line[tol][0]=i;
                  Line[tol][1]=j;
                  Line[tol][2]=backup[i][j]-g[i][j];
                  tol++;
              }
          }
        printf("%d\n",tol);
        for(int i=0;i<tol;i++)
        {
            printf("%d %d %d\n",Line[i][0],Line[i][1],Line[i][2]);
        }
    }
    return 0;
}