【BZOJ3931】【CQOI2015】网络吞吐量 最短路+网络流

链接：

#include 
int main()
{
    puts("转载请注明出处[vmurder]谢谢");
    puts("网址：blog.csdn.net/vmurder/article/details/44922137");
}

题解：

……两遍最短路然后判断哪些边可以在某条最短路上，然后加到网络流图中。
然后题意是一个点经过流量有限制，拆点就好。

然后有重边Qwq（调了好久。。。）
然后或许有自环，不过这并不影响什么。

代码：

调试过程中代码风格已经狗一样。
看到 long long i;i<=n;i++我就想哭

#include 
#include 
#include 
#include 
#include 
#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;ilong long t=INF,x=0;
        for(long long j=1;j<=n;j++)
            if(!v[j]&&f[p][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]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<return 0;
}