【BZOJ3931】【CQOI2015】网络吞吐量 最短路+网络流
2015-04-07 16:28
169 查看
链接:
#include <stdio.h> int main() { puts("转载请注明出处[vmurder]谢谢"); puts("网址:blog.csdn.net/vmurder/article/details/44922137"); }
题解:
……两遍最短路然后判断哪些边可以在某条最短路上,然后加到网络流图中。然后题意是一个点经过流量有限制,拆点就好。
然后有重边Qwq(调了好久。。。)
然后或许有自环,不过这并不影响什么。
代码:
调试过程中代码风格已经狗一样。看到 long long i;i<=n;i++我就想哭
#include <queue> #include <cstdio> #include <cstring> #include <iostream> #include <algorithm> #define N 1200 #define M 201000 #define INF 0x3f3f3f3f3f3f3f3fll using namespace std; long long map ; long long f[2] ; bool v ; long long n; void dij(long long s,bool p) { memset(v,0,sizeof v); memset(f[p],0x3f,sizeof f[p]); f[p][s]=0; for(long long i=1;i<n;i++) { long long t=INF,x=0; for(long long j=1;j<=n;j++) if(!v[j]&&f[p][j]<t) t=f[p][j],x=j; if(!x)break; v[x]=1; for(long long j=1;j<=n;j++)if(!v[j]) f[p][j]=min(f[p][j],f[p][x]+map[x][j]); } } struct Eli { long long v,next; long long len; }e[M]; long long head ,cnt; inline void add(long long u,long long v,long long len) { e[++cnt].v=v; e[cnt].len=len; e[cnt].next=head[u]; head[u]=cnt; } queue<long long>q; long long d ; bool bfs() { memset(d,0,sizeof d); while(!q.empty())q.pop(); q.push(1),d[1]=1; long long i,u,v; while(!q.empty()) { u=q.front(),q.pop(); for(i=head[u];i;i=e[i].next) { if(!d[v=e[i].v]&&e[i].len) { d[v]=d[u]+1; if(v==n)return 1; q.push(v); } } } return 0; } long long dinic(long long x,long long flow) { if(x==n)return flow; long long remain=flow,k; long long i,v; for(i=head[x];i&&remain;i=e[i].next) { if(d[v=e[i].v]==d[x]+1&&e[i].len) { k=dinic(v,min(remain,e[i].len)); if(!k)d[v]=0; e[i].len-=k,e[i^1].len+=k; remain-=k; } } return flow-remain; } long long flow ,maxflow; void build() { long long i,j,k,m; cin>>n>>m; memset(map,0x3f,sizeof map); while(m--) { cin>>i>>j>>k; map[i][j]=map[j][i]=min(map[i][j],k); } dij(1,0),dij(n,1); for(i=1;i<=n;i++)cin>>flow[i]; flow[1]=flow =INF; for(cnt=i=1;i<=n;i++) { add(i*2-1,i*2,flow[i]),add(i*2,i*2-1,0); for(j=1;j<=n;j++)if(map[i][j]<INF)if(f[0][i]+map[i][j]+f[1][j]==f[0] ) if(i!=j)add(i*2,j*2-1,INF),add(j*2-1,i*2,0); } n<<=1; } int main() { // freopen("network1.in","r",stdin); build(); while(bfs()) maxflow+=dinic(1,INF); cout<<maxflow<<endl; fclose(stdin); fclose(stdout); return 0; }
相关文章推荐
- 【BZOJ3931】【CQOI2015】网络吞吐量(最短路,网络流)
- [BZOJ3931]CQOI2015网络吞吐量|最短路|最大流
- bzoj 3931: [CQOI2015]网络吞吐量 -- 最短路+网络流
- [CQOI2015]网络吞吐量 - 最短路 - 网络流
- [bzoj3931][CQOI2015]网络吞吐量——最短路+网络流
- 【最短路】【最大流】bzoj3931 [CQOI2015]网络吞吐量
- BZOJ 3931: [CQOI2015]网络吞吐量【最短路+网络流
- [BZOJ3931] [CQOI2015]网络吞吐量 && 网络流
- [BZOJ3931][CQOI2015][网络流]网络吞吐量[水题]
- bzoj 3931: [CQOI2015]网络吞吐量(最短路+网络流)
- Bzoj3931 [CQOI2015]网络吞吐量
- bzoj3931【CQOI2015】网络吞吐量
- 3931: [CQOI2015]网络吞吐量【网络流】
- 【bzoj3931】【CQOI2015】【网络吞吐量】【spfa+最大流】
- bzoj3931 [CQOI2015]网络吞吐量(Dijkstra+最大流)
- BZOJ 3931: [CQOI2015]网络吞吐量( 最短路 + 最大流 )
- bzoj3931: [CQOI2015]网络吞吐量 网络流
- bzoj千题计划136:bzoj3931: [CQOI2015]网络吞吐量
- bzoj3931 [CQOI2015]网络吞吐量
- [bzoj3931][网络流][最短路]网络吞吐量