题目大意:
一些岛屿,之间有一些航线(双向的),除了起点(最西边那个)和终点(最东边那个),航线不相交。已知每条航线的客流量的上限,问到从起点到终点最大的客流量。
裸的求最大流的问题。
数据规模比较大。
用了递归的Dinic,人工扩栈后c++提交,压着时限过了。
#include<iostream>
#include<cstdio>
#include<cmath>
#include<map>
#include<set>
#include<algorithm>
#include<queue>
#include<stack>
#include<cstdlib>
#include<cstring>
#include<vector>
#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<int> 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(x<west) west=x,s=i;
if(x>east) east=x,t=i;
}
for(i=0;i<m;++i){
scanf("%d%d%d",&a,&b,&c);
add(a,b,c);
}
printf("%d\n",dinic());
}
return 0;
}
下面是用ISAP过的,运行时间是上面的一半。
#include<iostream>
#include<cstdio>
#include<cmath>
#include<map>
#include<set>
#include<algorithm>
#include<queue>
#include<stack>
#include<cstdlib>
#include<cstring>
#include<vector>
#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]<nv)
{
if(u==sink)
{
int f=inf,neck;
for(i=sourse;i!=sink;i=edge[cur[i]].to)
{
if(f>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(x<west) west=x,sourse=i;
if(x>east) east=x,sink=i;
}
nv=sink+1;
for(i=0;i<m;++i){
scanf("%d%d%d",&a,&b,&c);
add(a,b,c);
}
printf("%d\n",ISAP());
}
return 0;
}