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

考点在问题转化,然后就跑一个你喜欢的最大流算法(二分图可以啵?)

【poj 1087 a plug for UNIX】-LMLPHP

再来一个例子吧:

【poj 1087 a plug for UNIX】-LMLPHP

【纯手绘大米饼图片】

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

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

【poj 1087 a plug for UNIX】-LMLPHP
 1 #include<stdio.h>
2 #include<algorithm>
3 #include<queue>
4 #include<cstring>
5 #define go(i,a,b) for(int i=a;i<=b;i++)
6 #define fo(i,a,x) for(int i=a[x];i>-1;i=e[i].next)
7 #define mem(a,b) memset(a,b,sizeof(a))
8 #define inf 1000000000
9 using namespace std;const int N=90000
;
10 struct E{int v,next,flow,cap;}e[N*20
];
11 char socket[N][30],device[N][2][30],change[N][2][30
];
12 int
n,m,s,t,head[N],k,K,d[N],preE[N],preN[N],cur[N],num[N];
13 void ADD(int u,int v,int flow,int cap){e[K]=(E){v,head[u],flow,cap};head[u]=K++
;}
14 void
Big_Riced_Pancake()
15 {//m(device)--->k(change)--->n(socket)
16 go(i,1,m){ADD(s,i,0,1);ADD(i,s,0,0
);
17 go(j,1,k)if(!strcmp(device[i][1],change[j][0]))ADD(i,j+m,0,inf),ADD(j+m,i,0,0
);
18 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
);}
19 go(i,1,n)ADD(m+k+i,t,0,1),ADD(t,m+k+i,0,0);go(i,1
,k){
20 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
);
21 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
);}
22
}
23 void BFS(){queue<int>q;bool vis[N]={0};q.push(t);vis[t]=1;d[t]=0
;
24 while(!q.empty()){int u=
q.front();q.pop();fo(i,head,u)
25 {int v=e[i].v;if(!vis[v]&&!e[i].cap)vis[v]=1,d[v]=d[u]+1
,q.push(v);}}
26
}
27 int aug(){int u,a=
inf;
28 u=t;while(u!=s){int i=preE[u];a=min(a,e[i].cap-e[i].flow);u=
preN[u];}
29 u=t;while(u!=s){int i=preE[u];e[i].flow+=a;e[i^1].flow-=a;u=preN[u];}return
a;
30
}
31 int main(){scanf("%d",&n);go(i,1,n)scanf("%s"
,socket[i]);
32 scanf("%d",&m);go(i,1,m)go(j,0,1)scanf("%s"
,device[i][j]);
33 scanf("%d",&k);go(i,1,k)go(j,0,1)scanf("%s"
,change[i][j]);
34 mem(head,-1);K=0;s=0;t=n+m+k+1
;Big_Riced_Pancake();
35 BFS();go(i,0,t)num[d[i]]++,cur[i]=head[i];int u=s,flow=0
;
36 while(d[s]<t+1){u==t?flow+=aug(),u=s:1;bool retreat=1
;
37 fo(i,cur,u){int v=e[i].v;if(e[i].cap>e[i].flow&&d[u]==d[v]+1
)
38 {retreat=0;preN[v]=u;cur[u]=preE[v]=i;u=v;break;}}if(!retreat)continue
;
39 int Min=
t;
40 fo(i,head,u){int v=e[i].v;if(e[i].cap>e[i].flow)Min=
min(Min,d[v]);}
41 if(!(--num[d[u]]))break;num[d[u]=Min+1]++;cur[u]=head[u];if(u!=s)u=
preN[u];
42 }printf("%d\n",m-flow);return 0
;
43 }//Paul_Guderian

【ISAP-37msAC大米饼】

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

05-11 12:54