HDU - 3416 Marriage Match IV(最大流+最短路)

题目链接：点击查看

题目大意：给出一个由n个点和m条边组成的无环有向图，有自环，有重边，现在给出起点st和终点ed，问从st到ed共有多少条最短路，每条路最多只能经过一次

题目分析：看完题意后感觉是要基于最短路的最大流，所以直接先跑一遍spfa再说，为什么不用迪杰斯特拉呢，主要还是因为spfa好写。。太懒了，有向图的话，存边时记得顺便存一下反向边，这样起点和终点分别跑一次spfa，然后就可以遍历一遍所有边，从而找出属于最短路的上的所有边了，把这些边都拎出来，建立最大流的图，因为每条边至多经过一次，所以我们将其流量设置为1，最后跑一边最大流就是答案了，因为是比较裸的最大流，也没什么好说的了，dinic飘过

代码：

#include<iostream> #include<cstdlib> #include<string> #include<cstring> #include<cstdio> #include<algorithm> #include<climits> #include<cmath> #include<cctype> #include<stack> #include<queue> #include<list> #include<vector> #include<set> #include<map> #include<sstream> using namespace std;typedef long long LL;const int inf=0x3f3f3f3f;const int N=1e3+100;struct Node {int to,w;Node(int TO,int W){to=TO;w=W;} };vector<Node>node[2][N];//0:正向边 1:反向边 struct Edge {int to,w,next; }edge[N*N];//边数int head[N],cnt;void addedge(int u,int v,int w) {edge[cnt].to=v;edge[cnt].w=w;edge[cnt].next=head[u];head[u]=cnt++;edge[cnt].to=u;edge[cnt].w=0;//反向边边权设置为0edge[cnt].next=head[v];head[v]=cnt++; }int d[N],now[N];//深度 当前弧优化bool bfs(int s,int t)//寻找增广路 {memset(d,0,sizeof(d));queue<int>q;q.push(s);now[s]=head[s];d[s]=1;while(!q.empty()){int u=q.front();q.pop();for(int i=head[u];i!=-1;i=edge[i].next){int v=edge[i].to;int w=edge[i].w;if(d[v])continue;if(!w)continue;d[v]=d[u]+1;now[v]=head[v];q.push(v);if(v==t)return true;}}return false; }int dinic(int x,int t,int flow)//更新答案 {if(x==t)return flow;int rest=flow,i;for(i=now[x];i!=-1&&rest;i=edge[i].next){int v=edge[i].to;int w=edge[i].w;if(w&&d[v]==d[x]+1){int k=dinic(v,t,min(rest,w));if(!k)d[v]=0;edge[i].w-=k;edge[i^1].w+=k;rest-=k;}}now[x]=i;return flow-rest; }void init() {for(int i=0;i<N;i++)for(int j=0;j<2;j++)node[j][i].clear();memset(head,-1,sizeof(head));cnt=0; }int solve(int st,int ed) {int ans=0,flow;while(bfs(st,ed))while(flow=dinic(st,ed,inf))ans+=flow;return ans; }bool vis[N];int dis[2][N];void spfa(int st,int k) {memset(vis,false,sizeof(vis));memset(dis[k],inf,sizeof(dis[k]));queue<int>q;q.push(st);vis[st]=true;dis[k][st]=0;while(q.size()){int u=q.front();q.pop();vis[u]=false;for(int i=0;i<node[k][u].size();i++){int v=node[k][u][i].to;int w=node[k][u][i].w;if(dis[k][v]>dis[k][u]+w){dis[k][v]=dis[k][u]+w;if(!vis[v]){vis[v]=true;q.push(v);}}}} }int main() { // freopen("input.txt","r",stdin); // ios::sync_with_stdio(false);int w;cin>>w;while(w--){init();int n,m;scanf("%d%d",&n,&m);while(m--){int u,v,w;scanf("%d%d%d",&u,&v,&w);if(u==v)//自环 continue;node[0][u].push_back(Node(v,w));node[1][v].push_back(Node(u,w));}int st,ed;scanf("%d%d",&st,&ed);spfa(st,0);spfa(ed,1);for(int i=1;i<=n;i++){for(int j=0;j<node[0][i].size();j++){int u=i;int v=node[0][i][j].to;int w=node[0][i][j].w;if(dis[0][u]+dis[1][v]+w==dis[0][ed])addedge(u,v,1);}}printf("%d\n",solve(st,ed));}return 0; }

 

总结

以上是生活随笔为你收集整理的HDU - 3416 Marriage Match IV(最大流+最短路)的全部内容，希望文章能够帮你解决所遇到的问题。

如果觉得生活随笔网站内容还不错，欢迎将生活随笔推荐给好友。