题目大意：
一些岛屿，之间有一些航线(双向的)，除了起点（最西边那个）和终点（最东边那个），航线不相交。已知每条航线的客流量的上限，问到从起点到终点最大的客流量。
裸的求最大流的问题。
数据规模比较大。
用了递归的Dinic，人工扩栈后c++提交，压着时限过了。
#include #include #include #include #include #include #include #include #include #include #include #pragma comment(linker, "/STACK:1024000000,1024000000") using namespace std; typedef long long LL; typedef unsigned long long uLL; typedef __int64 LI; typedef unsigned __int64 uLI; typedef unsigned int uI; typedef double db; #define maxn 100005 #define inf 0x3f3f3f3f struct Edge { int to,cap,next; }edge[maxn<<1]; int cnt,head[maxn],d[maxn],s,t; inline void add(int u,int v,int w) { edge[cnt].to=v; edge[cnt].cap=w; edge[cnt].next=head[u]; head[u]=cnt++; edge[cnt].to=u; edge[cnt].cap=w; edge[cnt].next=head[v]; head[v]=cnt++; } inline bool bfs() //BFS构造层次图 { int u,i; queue q; q.push(s); memset(d,0,sizeof(d)); d[s]=1; while(!q.empty()) { u=q.front();q.pop(); if(u==t) return true; //汇点的层次被算出即可停止，因为根据DFS的规则，和汇点同层或更下一层的结点是不可能走到汇点的。 for(i=head[u];~i;i=edge[i].next){ int v=edge[i].to; if(!d[v]&&edge[i].cap) { d[v]=d[u]+1; q.push(v); } } } return false; } int dfs(int u,int a) //当前结点、目前为止所有弧的最小残量 { int flow=0,f,i; if(u==t||a==0) return a; for(i=head[u];~i;i=edge[i].next) { int v=edge[i].to; if(edge[i].cap&&d[v]==d[u]+1) //满足容量不为0且另一端是下一层次 { f=dfs(v,min(a,edge[i].cap)); edge[i].cap-=f; //削减路径上各边的容量 edge[i^1].cap+=f;//添加反向边 flow+=f; //增加总流量 a-=f; //该结点之前的弧的最小残量消减f后的值 if(!a) break; //若为0，则无需再从该结点寻找增广路，因为位于该结点之前的某条弧削减f后会为0，导致无法从源点再增广到汇点 } } if(flow==0) d[u]=0; //从该点出发找不到增广路，则将该点从层次图中去掉 return flow; } int dinic() { int ans=0; while(bfs()) ans+=dfs(s,inf); return ans; } int main() { int a,b,c,T,x,y,i,n,m; scanf("%d",&T); while(T--){ scanf("%d%d",&n,&m); memset(head,-1,sizeof(head)); cnt=0; int west=inf,east=-inf; for(i=1;i<=n;++i){ scanf("%d%d",&x,&y); if(xeast) east=x,t=i; } for(i=0;i
下面是用ISAP过的，运行时间是上面的一半。
#include #include #include #include #include #include #include #include #include #include #include #pragma comment(linker, "/STACK:1024000000,1024000000") using namespace std; typedef long long LL; typedef unsigned long long uLL; typedef unsigned int uI; typedef double db; #define maxn 100005 #define inf 0x3f3f3f3f #define maxm 200005 #define maxq 100005 struct Edge{ int to,cap,next; }edge[maxm]; int Q[maxq],qhead,tail;//队列 int cnt,head[maxn],d[maxn],cur[maxn],pre[maxn],num[maxn]; int sourse,sink,nv;//源点、汇点、编号修改的上限 int n,m; inline void add(int u,int v,int w) { edge[cnt].to=v; edge[cnt].cap=w; edge[cnt].next=head[u]; head[u]=cnt++; edge[cnt].to=u; edge[cnt].cap=w; edge[cnt].next=head[v]; head[v]=cnt++; } void bfs() { memset(num,0,sizeof(num)); memset(d,-1,sizeof(d)); d[sink]=0; num[0]=1; qhead=tail=0; Q[tail++]=sink; while(qhead!=tail) { int u=Q[qhead++]; for(int i=head[u];~i;i=edge[i].next) { int v=edge[i].to; if(~d[v]) continue; d[v]=d[u]+1; Q[tail++]=v; num[d[v]]++; } } } int ISAP() { int i; for(i=0;i<=n;++i) cur[i]=head[i]; bfs(); int flow=0,u=pre[sourse]=sourse; while(d[sink]edge[cur[i]].cap) { f=edge[cur[i]].cap; neck=i; } } for(i=sourse;i!=sink;i=edge[cur[i]].to) { edge[cur[i]].cap-=f; edge[cur[i]^1].cap+=f; } flow+=f; u=neck; } for(i=cur[u];~i;i=edge[i].next) if(d[edge[i].to]+1==d[u]&&edge[i].cap) break; if(~i) { cur[u]=i; pre[edge[i].to]=u; u=edge[i].to; } else{ if(0==(--num[d[u]])) break; int mind=nv; for(i=head[u];~i;i=edge[i].next) { if(edge[i].cap&&mind>d[edge[i].to]) { cur[u]=i; mind=d[edge[i].to]; } } d[u]=mind+1; num[d[u]]++; u=pre[u]; } } return flow; } int main() { int a,b,c,T,x,y,i; scanf("%d",&T); while(T--){ scanf("%d%d",&n,&m); memset(head,-1,sizeof(head)); cnt=0; int west=inf,east=-inf; for(i=1;i<=n;++i){ scanf("%d%d",&x,&y); if(xeast) east=x,sink=i; } nv=sink+1; for(i=0;i
