今天做了一個(gè)毒瘤題把dinic卡了，不想因?yàn)閮?yōu)化常數(shù)改變習(xí)慣性的寫法，干脆就直接去找個(gè)效率更高的算法拿來用多好

注意這套板子的n要在全局變量中用，并且n是所有點(diǎn)的數(shù)量(包括源點(diǎn)和匯點(diǎn))，所以在有的時(shí)候需要n+=2，與我的習(xí)慣還是有一丟丟小不同，不過問題不大

%kuangbin大神

代碼：

const int N=1e5+100;//點(diǎn)數(shù)的最大值const int M=2e5+100;//邊數(shù)的最大值struct Node {int from,to,next;int cap; }edge[M];int tol,head[N],dep[N],gap[N];//gap[x]=y :說明殘留網(wǎng)絡(luò)中dep[i]==x的個(gè)數(shù)為yint n;//n是總的點(diǎn)的個(gè)數(shù)，包括源點(diǎn)和匯點(diǎn)void init() {tol=0;memset(head,-1,sizeof(head)); }void addedge(int u,int v,int w) {edge[tol].from=u;edge[tol].to=v;edge[tol].cap=w;edge[tol].next=head[u];head[u]=tol++;edge[tol].from=v;edge[tol].to=u;edge[tol].cap=0;edge[tol].next=head[v];head[v]=tol++; } void BFS(int start,int end) {memset(dep,-1,sizeof(dep));memset(gap,0,sizeof(gap));gap[0]=1;int que[N];int front,rear;front=rear=0;dep[end]=0;que[rear++]=end;while(front!=rear){int u=que[front++];if(front==N)front=0;for(int i=head[u];i!=-1;i=edge[i].next){int v=edge[i].to;if(dep[v]!=-1)continue;que[rear++]=v;if(rear==N)rear=0;dep[v]=dep[u]+1;++gap[dep[v]];}} } int SAP(int start,int end) {int res=0;BFS(start,end);int cur[N];int S[N];int top=0;memcpy(cur,head,sizeof(head));int u=start;int i;while(dep[start]<n){if(u==end){int temp=inf;int inser;for(i=0;i<top;i++)if(temp>edge[S[i]].cap){temp=edge[S[i]].cap;inser=i;}for(i=0;i<top;i++){edge[S[i]].cap-=temp;edge[S[i]^1].cap+=temp;}res+=temp;top=inser;u=edge[S[top]].from;}if(u!=end&&gap[dep[u]-1]==0)//出現(xiàn)斷層，無增廣路break;for(i=cur[u];i!=-1;i=edge[i].next)if(edge[i].cap!=0&&dep[u]==dep[edge[i].to]+1)break;if(i!=-1){cur[u]=i;S[top++]=i;u=edge[i].to;}else{int min=n;for(i=head[u];i!=-1;i=edge[i].next){if(edge[i].cap==0)continue;if(min>dep[edge[i].to]){min=dep[edge[i].to];cur[u]=i;}}--gap[dep[u]];dep[u]=min+1;++gap[dep[u]];if(u!=start)u=edge[S[--top]].from;}}return res; }

?

總結(jié)

以上是生活随笔為你收集整理的(转)网络流-最大流 SAP算法(模板)的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。