【POJ 1459 power network】
不可以理解的是,测评站上的0ms是怎么搞出来的。
这一题在建立超级源点和超级汇点后就变得温和可爱了。其实它本身就温和可爱。对比了能够找到的题解:
(1)艾德蒙·卡普算法(2)迪尼克算法(3)改进版艾德蒙·卡普算法(MY METHOD)
不去管那个0ms的吧,那么(3)号算法最为美妙【它的别名是:ISAP】,时间可观。
这个就算是一个ISAP的模板吧(除了输入的难看的几行外,其余均是标准的大米饼牌模板!)
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],v=e[i].v;i>-1;i=e[i].next,v=e[i].v)
7 #define mem(a,b) memset(a,b,sizeof(a))
8 using namespace std;const int N=5000;
9 struct E{int v,next,flow,cap;}e[N*N];
10 int n,n_power,n_user,m,S,T,head[N],k=0,d[N];
11 int num[N],cur[N],preE[N],preN[N];
12 void ADD(int u,int v,int flow,int cap){e[k]=(E){v,head[u],flow,cap};head[u]=k++;}
13 void BFS()
14 {
15 queue<int>q;bool vis[N]={0};q.push(T);d[T]=0;
16 while(!q.empty())
17 {
18 int u=q.front();q.pop();
19 fo(i,head,u)if(!e[i].cap&&!vis[v])
20 vis[v]=1,d[v]=d[u]+1,q.push(v);
21 }
22 }
23 int aug(){int u,a=2147483645;
24 u=T;while(u!=S){int i=preE[u];a=min(a,e[i].cap-e[i].flow);u=preN[u];}
25 u=T;while(u!=S){int i=preE[u];e[i].flow+=a;e[i^1].flow-=a;u=preN[u];}return a;}
26 int main(){while(~scanf("%d%d%d%d",&n,&n_power,&n_user,&m))
27 {
28 S=0;T=n+1;char _;mem(head,-1);
29 go(i,1,m){int u,v,cap;scanf(" %c%d%c%d%c%d",&_,&u,&_,&v,&_,&cap);
30 u++,v++;ADD(u,v,0,cap);ADD(v,u,0,0);}
31 go(i,1,n_power){int v,cap;scanf(" %c%d%c%d",&_,&v,&_,&cap);
32 v++;ADD(S,v,0,cap);ADD(v,S,0,0);}
33 go(i,1,n_user){int u,cap;scanf(" %c%d%c%d",&_,&u,&_,&cap);
34 u++;ADD(u,T,0,cap);ADD(T,u,0,0);}BFS();mem(num,0);
35
36 go(i,0,n+1)num[d[i]]++,cur[i]=head[i];
37 int u=S,flow=0;
38 while(d[S]<T+1)
39 {
40 u==T?flow+=aug(),u=S:1;bool retreat=1;
41 fo(i,cur,u) if(e[i].cap>e[i].flow&&d[u]==d[v]+1)
42 {retreat=0;cur[u]=preE[v]=i;preN[v]=u;u=v;break;}
43
44 if(!retreat)continue;
45
46 int Min=T;fo(i,head,u){if(e[i].cap>e[i].flow)Min=min(Min,d[v]);}
47 if(!(--num[d[u]]))break;num[d[u]=Min+1]++;
48 cur[u]=head[u];u==S?1:u=preN[u];
49 }
50 printf("%d\n",flow);
51 }return 0;}//Paul_Guderian
网络流两个学习要点:学会转化问题,学会几个模型?你是在质疑神奇大米饼的权威吗?那就给你个大米饼!
ZBH不给我玩以撒的结合…我写了这篇博客来表达愤怒!