uva1659(最大费用循环流)

紫书上的一道题，做法见紫书P378,这篇博客用的第二种方法，关于正确性的证明，画图可以发现如果一个环是负环，跑最小费用流跑出的是环上的所有正边，再减去负边和即为跑一遍的负权，如果是正环，最小费用流即为负权和的相反数，再加上负权和即为0，即没有计入答案。可能这种方法在处理负权时有一定的普遍性，还是要多做题多体会啊。。。。

PS:今天费用流又写错了，写的是dinic中dfs的那种，保留了那部分代码以便日后修改，，还是记录最短路径稳啊。

#include<iostream>
#include<cstring>
#include<cmath>
#include<cstdio>
#include<algorithm>
#include<vector>
using namespace std;
const int maxn=,inf=1e9;
struct node{
    double x,y;
}a[maxn];
double calc(node a,node b){
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
struct edg{
    int nxt,to,f;
    double c;
}e[maxn];
int pe[maxn],pv[maxn],du[maxn],res,t=,n,m,last[maxn],S,T,vis[maxn],q[maxn],head,tail;
double dis[maxn],Min;
void add(int x,int y,int z,double zz){
    ++t;e[t].nxt=last[x];last[x]=t;e[t].to=y;e[t].f=z;e[t].c=zz;
    ++t;e[t].nxt=last[y];last[y]=t;e[t].to=x;e[t].f=;e[t].c=-zz;
}
int dfs(int x,int h){
    if(x==T){Min+=double(h)*dis[T];return h;}
    int tmp=,cp;
    for(int i=last[x];i;i=e[i].nxt){
        int v=e[i].to;
        if(e[i].f&&dis[v]==dis[x]+e[i].c){
            cp=dfs(v,min(h-tmp,e[i].f));
            e[i].f-=cp;e[i^].f+=cp;tmp+=cp;
            if(h==tmp){return h;}
        }
    }
    return tmp;
}
void solve(){
    while(){
        for(int i=;i<=n+;++i)dis[i]=inf;
        memset(vis,,sizeof(vis));head=tail=;
        q[++tail]=S;dis[S]=;vis[S]=;
        while(head!=tail){
            head=(head+)%maxn;
            int u=q[head];
            for(int i=last[u];i;i=e[i].nxt){
                int v=e[i].to;
                if(e[i].f&&dis[v]>dis[u]+e[i].c){
                    dis[v]=dis[u]+e[i].c;
                    pe[v]=i;pv[v]=u;
                    if(!vis[v]){
                        vis[v]=;tail=(tail+)%maxn;q[tail]=v;
                    }
                }
            }
            vis[u]=;
        }
        //cout<<dis[T]<<endl;
        //system("pause");
        if(dis[T]==inf)break;
        int tmp=inf;
        for(int u=T;u!=S;u=pv[u])tmp=min(tmp,e[pe[u]].f);
        for(int u=T;u!=S;u=pv[u]){
            e[pe[u]].f-=tmp;
            e[pe[u]^].f+=tmp;
        }
        Min+=tmp*dis[T];
        //res+=dfs(S,inf);
    }
}
double X,Y;
vector<int>bel[maxn];
int cc,cas;
int main(){
    while(scanf("%d",&n)!=EOF){
        if(!n)break;
        scanf("%lf%lf",&X,&Y);
        ++cas;
        double sum=;
        t=;memset(last,,sizeof(last));
        memset(du,,sizeof(du));res=;Min=;
        for(int i=;i<=n;++i)bel[i].clear();
        for(int i=;i<=n;++i){
            scanf("%lf%lf",&a[i].x,&a[i].y);
            while(scanf("%d",&cc)!=EOF){
                if(!cc)break;
                bel[i].push_back(cc);
            }
        }
        for(int i=;i<=n;++i){
            int len=bel[i].size();
            for(int j=;j<len;++j){
                double vv=Y-calc(a[i],a[bel[i][j]])*X;
                if(vv<){
                    sum+=vv;
                    du[bel[i][j]]++;du[i]--;
                    add(bel[i][j],i,,-vv);
                }
                else{add(i,bel[i][j],,vv);}
            }
        }
        S=n+,T=n+;
        for(int i=;i<=n;++i){
            if(!du[i])continue;
            if(du[i]>)add(S,i,du[i],);
            else{add(i,T,-du[i],);}
        }
        solve();
        printf("Case %d: %.2lf\n",cas,fabs(sum+Min));
    }
    //system("pause");
    return ;
}
/*
4 5 1
0 0 2 3 0
1 0 3 4 0
1 1 4 0 0
1 1 0
1 2 1
0 0 0
10 7 2
0 0 2 4 0
5 0 3 0
5 10 4 10 0
2 3 5 0
7 5 6 0
0 11 1 0
8 0 10 5 0
18 3 7 0
14 5 8 1 0
12 9 9 0
0
*/