2018.10.10【ZOJ2760】How Many Shortest Path（最短路）（网络流）

传送门

---

解析：

题意简述：

给出一个有向图，问从$S$到$T$的边不相交的最短路有多少条。

思路：

首先肯定要求最短路。

然后这是一个匹配问题。
我们考虑用网络流来完成。

匹配问题的一般套路就是将需要匹配的东西在流网络中用一条边串起来，设容量为1，求最大流就是匹配的数量了。

那么这道题中，什么是合法可匹配的东西？

由于我们已经求了最短路了，
那么满足$dist[v]=dist[u]+w[e]$的边就是最短路上的边，连起来就好了。

---

代码：

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define re register
#define gc getchar
#define pc putchar
#define cs const

inline int getint(){
	re int num;
	re char c;
	re bool f=0;
	while(!isdigit(c=gc()))if(c=='-')f=1;num=c^48;
	while(isdigit(c=gc()))num=(num<<1)+(num<<3)+(c^48);
	return f?-num:num;
}

cs int N=103;
cs int M=100005;
cs int INF=0x3f3f3f3f;
int n;
int last[N],nxt[M<<1],to[M<<1],ecnt=1;
int cap[M<<1];
inline void addedge(int u,int v,int val){
	nxt[++ecnt]=last[u],last[u]=ecnt,to[ecnt]=v,cap[ecnt]=val;
	nxt[++ecnt]=last[v],last[v]=ecnt,to[ecnt]=u,cap[ecnt]=0;
}

int lev[N],cur[N];
inline bool BFS(cs int &S,cs int &T){
	memset(lev,-1,sizeof lev);
	queue<int> q;
	q.push(S),lev[S]=0,cur[S]=last[S];
	while(!q.empty()){
		int u=q.front();
		q.pop();
		for(int re e=last[u],v=to[e];e;v=to[e=nxt[e]]){
			if(cap[e]&&lev[v]==-1){
				lev[v]=lev[u]+1;
				if(v==T)return true;
				cur[v]=last[v];
				q.push(v);
			}
		}
	}
	return false;
}

inline int Dinic(cs int &u,cs int &flow,cs int &T){
	if(u==T)return flow;
	int ans=0;
	for(int &e=cur[u],v=to[e];e;v=to[e=nxt[e]]){
		if(cap[e]&&lev[v]>lev[u]){
			int delta=Dinic(v,min(flow-ans,cap[e]),T);
			if(delta){
				cap[e]-=delta;
				cap[e^1]+=delta;
				ans+=delta;
				if(ans==flow)return ans;
			}
		}
	}
	return ans;
}

inline int maxflow(cs int &S,cs int &T){
	int ans=0;
	while(BFS(S,T))ans+=Dinic(S,INF,T);
	return ans;
}

int Last[N],Nxt[M],To[M],Ecnt;
int W[M];
inline void Addedge(int u,int v,int val){
	Nxt[++Ecnt]=Last[u],Last[u]=Ecnt,To[Ecnt]=v,W[Ecnt]=val;
}

int dist[N];
bitset<N> vis;
inline void Dijkstra(cs int &S){
	memset(dist,0x3f,sizeof dist);
	vis.reset();
	priority_queue<pair<int,int>,vector<pair<int,int> > ,greater<pair<int,int> > > q; 
	q.push(make_pair(dist[S]=0,S));
	while(!q.empty()){
		int u=q.top().second;
		q.pop();
		if(vis[u])continue;
		vis[u]=true;
		for(int re e=Last[u],v=To[e];e;v=To[e=Nxt[e]])
		if(dist[v]>dist[u]+W[e])q.push(make_pair(dist[v]=dist[u]+W[e],v));
	}
}

inline void init(){
	ecnt=1;
	Ecnt=0;
	memset(Last,0,sizeof Last);
	memset(last,0,sizeof last);
}

inline void build(){
	for(int re u=1;u<=n;++u)
	for(int re e=Last[u];e;e=Nxt[e])
	if(dist[To[e]]==dist[u]+W[e])addedge(u,To[e],1);
}

signed main(){
	while(~scanf("%d",&n)){
		init();
		for(int re i=1;i<=n;++i)
		for(int re j=1;j<=n;++j){
			int w=getint();
			if(w>=0)Addedge(i,j,w);
		}
		int S=getint()+1,T=getint()+1;
		if(S==T){puts("inf");continue;}
		Dijkstra(S);
		build();
		cout<<maxflow(S,T);
		pc('\n');
	}
	return 0;
}

转载于:https://www.cnblogs.com/zxyoi/p/10047239.html