【poj 1087 a plug for UNIX】

在大米饼的帮助下，终于找到了大米饼程序中如同大米饼一般的错误！

考点在问题转化，然后就跑一个你喜欢的最大流算法（二分图可以啵？）

再来一个例子吧：

【纯手绘大米饼图片】

其中有的边权是1，否则就是inf，所以就将问题转化为求超级源点(0)到超级汇点(13)的最大流。我依旧使用ISAP，很开心用ISAP做完了所有老师要求用迪尼克或者艾德蒙·卡普算法做的几个题目，开心点在哪里呢…在乎大米饼之中也。最后一个困扰了me20min的错误是：

把device[i][1]错写成device[i][2],使得CMP函数失效。

#include<stdio.h>
#include<algorithm>
#include<queue>
#include<cstring>
#define go(i,a,b) for(int i=a;i<=b;i++)
#define fo(i,a,x) for(int i=a[x];i>-1;i=e[i].next)
#define mem(a,b) memset(a,b,sizeof(a))
#define inf 1000000000
using namespace std;const int N=90000;
struct E{int v,next,flow,cap;}e[N*20];
char socket[N][30],device[N][2][30],change[N][2][30];
int n,m,s,t,head[N],k,K,d[N],preE[N],preN[N],cur[N],num[N];
void ADD(int u,int v,int flow,int cap){e[K]=(E){v,head[u],flow,cap};head[u]=K++;}
void Big_Riced_Pancake()
{//m(device)--->k(change)--->n(socket) 
    go(i,1,m){ADD(s,i,0,1);ADD(i,s,0,0);
    go(j,1,k)if(!strcmp(device[i][1],change[j][0]))ADD(i,j+m,0,inf),ADD(j+m,i,0,0);
    go(j,1,n)if(!strcmp(device[i][1],socket[j]))ADD(i,m+k+j,0,1),ADD(m+k+j,i,0,0);}           
    go(i,1,n)ADD(m+k+i,t,0,1),ADD(t,m+k+i,0,0);go(i,1,k){
    go(j,1,k)if((i!=j)&&(!strcmp(change[i][1],change[j][0])))ADD(i+m,j+m,0,inf),ADD(j+m,i+m,0,0);
    go(j,1,n)if(!strcmp(change[i][1],socket[j]))ADD(i+m,m+k+j,0,inf),ADD(m+k+j,i+m,0,0);}   
}
void BFS(){queue<int>q;bool vis[N]={0};q.push(t);vis[t]=1;d[t]=0;
    while(!q.empty()){int u=q.front();q.pop();fo(i,head,u)
    {int v=e[i].v;if(!vis[v]&&!e[i].cap)vis[v]=1,d[v]=d[u]+1,q.push(v);}}
}
int aug(){int u,a=inf;
    u=t;while(u!=s){int i=preE[u];a=min(a,e[i].cap-e[i].flow);u=preN[u];}
    u=t;while(u!=s){int i=preE[u];e[i].flow+=a;e[i^1].flow-=a;u=preN[u];}return a;
}
int main(){scanf("%d",&n);go(i,1,n)scanf("%s",socket[i]);
    scanf("%d",&m);go(i,1,m)go(j,0,1)scanf("%s",device[i][j]);
    scanf("%d",&k);go(i,1,k)go(j,0,1)scanf("%s",change[i][j]);
    mem(head,-1);K=0;s=0;t=n+m+k+1;Big_Riced_Pancake();
    BFS();go(i,0,t)num[d[i]]++,cur[i]=head[i];int u=s,flow=0;
    while(d[s]<t+1){u==t?flow+=aug(),u=s:1;bool retreat=1;
        fo(i,cur,u){int v=e[i].v;if(e[i].cap>e[i].flow&&d[u]==d[v]+1)
        {retreat=0;preN[v]=u;cur[u]=preE[v]=i;u=v;break;}}if(!retreat)continue;
        int Min=t;
        fo(i,head,u){int v=e[i].v;if(e[i].cap>e[i].flow)Min=min(Min,d[v]);}
        if(!(--num[d[u]]))break;num[d[u]=Min+1]++;cur[u]=head[u];if(u!=s)u=preN[u];       
    }printf("%d
",m-flow);return 0;
}//Paul_Guderian

感谢帮我给手稿拍照的同学！