UVA 10594-Date Flow(无向图的最小费用网络流+题目给的数据有误）

    int sa=a+skip;
    int sb=b+skip;
    save[a]=edge(from,to,cap,0,cost);
    save[b]=edge(to,from,0,0,-cost);
    save[sa]=edge(to,from,cap,0,cost);
    save[sb]=edge(from,to,0,0,-cost);
    G[from].push_back(a); G[from].push_back(sb);  // 建立临接表
    G[to].push_back(b);   G[to].push_back(sa);
}
LL add[maxn],C[maxn];   // add[]是增加的流量
int fa[maxn];           // 保存父亲边编号
bool inq[maxn];         // 入队标记
inline void init(int be)    // 初始化
{
    memset(add,0,sizeof(add));
    memset(inq,false,sizeof(inq));
    add[be]=INF;
    fa[be]=0;
    inq[be]=true;
    for(int i=0;i<=N;i++)  C[i]=INF;
    C[be]=0;
}
void MFMC(int be,int en)
{
    LL cnt=0;
    LL ret=0;           //  cnt是最大流量，ret是最小费用
    while(true)
    {
        init(be);
        queue<int> que;
        que.push(be);
        while(!que.empty())
        {
            int from=que.front();  que.pop();
            inq[from]=false;                   // spfa跑最小费，入队标记改为false
            for(int i=0;i<G[from].size();i++)
            {
                int edge_id=G[from][i];
                edge& e=save[edge_id];
                int to=e.to;
                LL cap=e.cap,flow=e.flow,cost=e.cost;
                if(cap>flow&&C[to]>C[from]+cost)    //更新
                {
                    C[to]=C[from]+cost;
                    fa[to]=edge_id;                  // 保存编号
                    add[to]=min(add[from],cap-flow);
                    if(!inq[to]){ inq[to]=true; que.push(to); }  //入队
                }
            }
        }
        if(!add[en])  break;
        cnt+=add[en];
        ret+=add[en]*C[en];
        for(int st=en;st!=be;st=save[fa[st]].from)   // 往回找
        {
            int a=fa[st],b=fa[st]^1;                 // 修改4条边
            int sa,sb;
            if(a<skip) { sa=a+skip;sb=sa^1; }
            else { sa=a-skip; sb=sa^1; }
            save[a].flow+=add[en];
            save[b].flow-=add[en];
            save[sa].flow-=add[en];
            save[sb].flow+=add[en];
        }
    }
    if(cnt==D) printf("%lld\n",ret);    
    else  printf("Impossible.\n");
}
int main()
{
    while(cin>>N>>M)
    {
        for(int i=0;i<maxn;i++)  G[i].clear();
        for(int i=1;i<=M;i++) scanf("%d%d%lld",&u[i],&v[i],&waste[i]);
        scanf("%lld%lld",&D,&K);
        save[0]=edge(0,1,D,0,0);   // 建立0到的边
        save[1]=edge(1,0,0,0,0);
        G[0].push_back(0);
        G[1].push_back(1);
        edge_num=2;
        for(int i=1;i<=M;i++)  addedge(u[i],v[i],K,waste[i]);  // 增加边
        MFMC(0,N);
    }
    return 0;
}