链接:
#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[N][N];
long long f[2][N];
bool v[N];
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[N],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[N];
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[N],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[n]=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][n])
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;
}
时间: 2024-10-05 04:09:24