poj 1273 SAP 模板刷题

来源：互联网发布：山东大学网络教育官网编辑：程序博客网时间：2024/05/01 08:22

使用了小王的SAP模板

如下：

#include<stdio.h>#include<memory.h>#include<iostream>using namespace std;const int MAXNUM = 2100000000 ; //权值上限const int MAXPT = 300 ; //顶点数上限const int MAXEG = MAXPT*MAXPT*2 ; //边数上限const int MAXQUE = 200000 ; // 队列长度/* s = 1 ; // 源点 t = n ; // 汇点*/class MNF_SAP{private: //int m,n; int s,t; int d[MAXPT]; //距离标号 int pre[MAXPT]; //前置顶点 int flow[MAXPT]; //到当前点为止，所有弧的最小值 int cur[MAXPT]; //当前弧 int cnt[MAXPT]; //k标号出现次数 int que[MAXPT*2],ft,rt; bool vis[MAXPT]; void init_d() { int i,u; memset(vis,0,sizeof(vis)); ft=rt=0; d[t]=0; vis[t]=true; que[++rt]=t; while(ft!=rt) { u=que[(++ft)%MAXQUE]; for(i=gra[u]; i; i=E[i].next) if(!vis[E[i].v] && !E[i].w) { d[E[i].v]=d[u]+1; vis[E[i].v]=true; que[(++rt)%MAXQUE]=E[i].v; } } }public: struct { int v,w,next; } E[MAXEG]; int n; int gra[MAXPT]; MNF_SAP(){} ~MNF_SAP(){} int SAP() { int u,v,i,maxflow=0; //总最大流 s=1; t=n; u=s; flow[s]=MAXNUM; for(i=1; i<=n; i++) cur[i]=gra[i]; //当前弧初始化 init_d(); for(i=1; i<=n; i++) cnt[d[i]]++; cnt[0]=n; while(d[s]<n) { for(i=cur[u]; i; i=E[i].next) //找允许弧 if(E[i].w && d[E[i].v]+1==d[u]) break; if(i) //存在允许弧 { cur[u]=i; //设置当前弧 v=E[i].v; if(E[i].w<flow[u]) flow[v]=E[i].w; else flow[v]=flow[u]; //标记当前顶点为止经过的最小弧 u=v; pre[v]=i; //前置顶点边号 if(u==t) { do { E[pre[u]].w-=flow[t]; //正向弧减a[t] E[pre[u]^1].w+=flow[t]; //通过异或操作找反向弧 u=E[pre[u]^1].v; } while(u!=s); maxflow+=flow[t]; //memset(flow,0,sizeof(flow)); flow[s]=MAXNUM; } } else //不存在允许弧 { if(--cnt[d[u]]==0) break; //间隙优化 d[u]=n; cur[u]=gra[u]; for(i=gra[u]; i; i=E[i].next) if(E[i].w && d[E[i].v]+1<d[u]) d[u]=d[E[i].v]+1; //修改距离标号为最小的非允许弧加1 cnt[d[u]]++; if(u!=s) u=E[pre[u]^1].v; } } return maxflow; }} ;int main(){ MNF_SAP ob ; int k ; int i,j; int num,mum; int x,y,z; while(cin>>num>>mum) { ob.n = mum ;k = 1 ; memset(ob.gra,0,sizeof(ob.gra)); for(i=0;i<num;i++) { scanf("%d%d%d",&x,&y,&z); ob.E[++k].v = y ; ob.E[k].next = ob.gra[x] ; ob.E[k].w = z ; ob.gra[x] = k ; ob.E[++k].v = x ; ob.E[k].next = ob.gra[y] ; ob.E[k].w = 0 ; ob.gra[y] = k ; } cout<<ob.SAP()<<endl; } return 0;}