[BZOJ4349]最小树形图

显然先选每个点都取一遍然后再取满次数最优，用最小树形图决定第一次取的顺序。

朱刘算法的流程是（总复杂度O(nm)）：

1.对除根外所有点，找到所有指向它的边中权值最小的那一条，记其权值为ind[]。

2.找到所有不包含根的、由(1)中找到的那些边构成的环，并将环缩点。若没有这样的环则结束。

3.将所有缩点后不是自环的边的权值，减去边的终点的ind。

主要思路就是贪心+调整，具体看代码实现。

#include<cstdio>
#include<algorithm>
#define rep(i,l,r) for (int i=(l); i<=(r); i++)
using namespace std;

const int N=61,M=21000;
const double inf=1e10;
int n=1,m,m1,n1,x,y,cnt,num[N],pos[N],pre[N],id[N],vis[N];
double a[N],ans,v,ind[N];
struct E{ int u,v; double w; }e[M];

double Zhuliu(int rt,int n,int m){
    int tn,tm; double res=0;
    while (1){
        rep(i,1,n) ind[i]=inf,pre[i]=id[i]=vis[i]=0;
        tn=tm=ind[rt]=0;
        rep(i,1,m) if (e[i].w<ind[e[i].v]) ind[e[i].v]=e[i].w,pre[e[i].v]=e[i].u;
        rep(i,1,n){
            int x=i; res+=ind[i];
            while (x!=rt && vis[x]!=i && !id[x]) vis[x]=i,x=pre[x];
            if (x!=rt && !id[x]){
                id[x]=++tn;
                for (int k=pre[x]; k!=x; k=pre[k]) id[k]=tn;
            }
        }
        if (!tn) break;
        rep(i,1,n) if (!id[i]) id[i]=++tn;
        rep(i,1,m) if (id[e[i].u]!=id[e[i].v]) e[++tm]=(E){id[e[i].u],id[e[i].v],e[i].w-ind[e[i].v]};
        n=tn; m=tm; rt=id[rt];
    }
    return res;
}

int main(){
    freopen("bzoj4349.in","r",stdin);
    freopen("bzoj4349.out","w",stdout);
    scanf("%d",&n1);
    rep(i,1,n1){
        scanf("%lf%d",&v,&x);
        if (x) pos[i]=++n,e[++m]=(E){1,n,v},a[n]=v,num[n]=x-1;
    }
    scanf("%d",&m1);
    rep(i,1,m1){
        scanf("%d%d%lf",&x,&y,&v);
        if (!pos[x] || !pos[y]) continue;
        e[++m]=(E){pos[x],pos[y],v}; a[pos[y]]=min(a[pos[y]],v);
    }
    rep(i,2,n) ans+=a[i]*num[i];
    printf("%.2lf
",Zhuliu(1,n,m)+ans);
    return 0;
}