[CQOI2015]网络吞吐量

这题我Debug了2小时，心情复杂。
这题建模不需要说了。
讲一下如何处理即可。
最短路？
SPFA啊
先把原图存下来。
咋判断一个点是不是在最短路径上？
枚举点的邻接点
若dis[to]=dis[from]+cost即为最短路
拆点，注意1,n限制不算。
还有要开long long，inf要大一些
Luogu
BZOJ
COGS

#include <cstdio>
#include <iostream>
#include <queue>
#include <cstring> 
#define il inline
#define ll long long
using namespace std;
ll inf;
const int maxm=110000;
ll cap[maxm],capx[maxm*2],dis[maxm];
int head[maxm],to[maxm*2],net[maxm*2],deep[maxm],cnt=1;
il void add(int x,int y,ll c){cnt++,to[cnt]=y,cap[cnt]=c,net[cnt]=head[x],head[x]=cnt;}
int headx[maxm],tox[maxm*2],netx[maxm*2],cntx=1;
il void addx(int x,int y,ll c){cntx++,tox[cntx]=y,capx[cntx]=c,netx[cntx]=headx[x],headx[x]=cntx;}
queue <int> dl;
ll mcf[maxm];
bool vis[maxm];
il bool BFS(int s,int t)
{
    while(!dl.empty()) dl.pop();
    memset(deep,-1,sizeof(deep));
    dl.push(s),deep[s]=0;
    while(!dl.empty())
    {
        int x=dl.front();dl.pop();
        for(int i=head[x];i;i=net[i])
         if(cap[i]>0&&deep[to[i]]==-1)
          dl.push(to[i]),deep[to[i]]=deep[x]+1;
    }
    return deep[t]==-1?0:1;
}
ll dfs(int now,ll flow,int t)
{
    if(now==t) return flow;
    ll w,used=0;
    for(int i=head[now];i;i=net[i])
    {
        int v=to[i];
        if(deep[v]==deep[now]+1&&cap[i])
        {
            w=dfs(v,min(flow-used,cap[i]),t);
            cap[i]-=w;
            cap[i^1]+=w;
            used+=w;
            if(used==flow) return flow;
        }
    }
    if(!used) deep[now]=-1;
    return used;
}
il ll dinic(int s,int t)
{
    ll maxflow=0;
    while(BFS(s,t)) maxflow+=dfs(s,inf,t);
    return maxflow;
}
inline void adx(int x,int y,int cax)
{
    add(x,y,cax),add(y,x,0);
}
il int read()
{
    int x=0,w=1;
    char ch=0;
    while(ch<'0'||ch>'9') 
    {
        if(ch=='-') w=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
     x=(x<<3)+(x<<1)+ch-'0',ch=getchar();
    return x*w;
}
il void SPFA(int s,int t)
{
    //printf("hehe\n");
    memset(dis,127,sizeof(dis));
    vis[s]=1,dis[s]=0;
    dl.push(s);
    while(!dl.empty())
    {
        int x=dl.front();
        dl.pop();
        vis[x]=0;
        //printf("%d %d\n",x,head[x]);
        for(int i=headx[x];i;i=netx[i])
        {
         if(dis[tox[i]]>dis[x]+capx[i])
         {
            dis[tox[i]]=dis[x]+capx[i];
            if(!vis[tox[i]]) dl.push(tox[i]),vis[tox[i]]=1;
         }
        }
    }
    //printf("%d\n",dis[t]);
}
int main()
{
    //freopen("cqoi15_network.in","r",stdin);
    //freopen("cqoi15_network.out","w",stdout);
    int n=read(),m=read();
    for(int i=1;i<=m;i++)
     {
        ll cos;
        int u=read(),v=read();
        scanf("%lld",&cos);
        addx(u,v,cos),addx(v,u,cos);
     }
    SPFA(1,n);
    for(int i=1;i<=n;i++)
    {
       cin>>mcf[i];
       inf=max(inf,mcf[i]);
    }
    inf++;
    for(int i=1;i<=n;i++)
    {
        adx(i,i+n,mcf[i]);
    }
    for(int i=1;i<=n;i++)
     for(int j=headx[i];j;j=netx[j])
      {
        int tmp=tox[j];
        if(dis[tox[j]]==dis[i]+capx[j]) adx(i+n,tox[j],inf);
      }
    return printf("%lld",dinic(1+n,n))*0;
}