【BZOJ2260】商店购物 & 【BZOJ4349】最小树形图 （朱刘算法）

Description

Grant是一个个体户老板，他经营的小店因为其丰富的优惠方案深受附近居民的青睐，生意红火。小店的优惠方案十分简单有趣。Grant规定：在一次消费过程中，如果您在本店购买了精制油的话，您购买香皂时就可以享受2.00元/块的优惠价；如果您在本店购买了香皂的话，您购买可乐时就可以享受1.50元/听的优惠价……诸如此类的优惠方案就是说：如果您在本店购买了商品A的话，您就可以以P元/件的优惠价格购买商品B（购买的数量不限）。有趣的是，你需要购买同样一些商品，由于不同的购买顺序，Grant老板可能会叫你付不同数量的钱。比如你需要一块香皂（原价2.50元）、一瓶精制油（原价10.00元）、一听可乐(原价1.80元)，如果你按照可乐，精制油，香皂这样的顺序购买的话，Grant老板会问你要13.80元；而如果你按照精制油，香皂，可乐这样的顺序购买的话，您只需付13.50元。
现在该村的居民请你编写一个程序，告诉你Grant小店商品的原价，所有优惠方案及所需的商品，计算至少需要花多少钱。不允许购买任何不需要的商品，即使这样做可能使花得钱更少。

Input

，第一行为一个整数n（1<=n<=50），表示Grant小店的商品种数。接下来是n行，其中第(i+1)行由一个实数Ci (0

#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<math.h>
#define ll long long
#define inf 1e10 
using namespace std;
#define N 55
#define M 21000
int n,m,m1,cnt;
double a[N],ans,in[N];
int num[N],pos[N],pre[N],id[N],vis[N];
struct edge
{
    int x,y;
    double v;
    edge(){}
    edge(int x,int y,double v):x(x),y(y),v(v){}
}e[M];
double zhuliu(int root,int n,int m)
{
    int tn,tm;
    double ret=0;
    while(1)
    {
        for(int i=1;i<=n;i++) in[i]=inf,pre[i]=0;
        for(int i=1;i<=m;i++)
        if(e[i].v<in[e[i].y])
        {
            in[e[i].y]=e[i].v;
            pre[e[i].y]=e[i].x;
        }
        tn=tm=in[root]=0;
        for(int i=1;i<=n;i++) id[i]=vis[i]=0;
        for(int i=1,v;i<=n;i++)
        {
            ret+=in[v=i];
            while(v!=root && vis[v]!=i && !id[v]) vis[v]=i,v=pre[v];
            if(v!=root && !id[v])
            {
                id[v]=++tn;
                for(int j=pre[v];j!=v;j=pre[j]) id[j]=tn;
            }
        }
        if(!tn) break;
        for(int i=1;i<=n;i++)
        if(!id[i]) id[i]=++tn;
        for(int i=1;i<=m;i++)
        if(id[e[i].x]!=id[e[i].y]) e[++tm]=edge(id[e[i].x],id[e[i].y],e[i].v-in[e[i].y]);
        n=tn;m=tm;root=id[root];
    }
    return ret;
}
int main()
{
    double v;
    int x,y;
    scanf("%d",&n);
    cnt=1;
    for(int i=1;i<=n;i++)
    {
        scanf("%lf%d",&v,&x);
        if(x)
        {
            pos[i]=++cnt;
            e[++m]=edge(1,cnt,v);
            a[cnt]=v;num[cnt]=x-1;
        }
    }
    n=cnt;
    scanf("%d",&m1);
    for(int i=1;i<=m1;i++)
    {
        scanf("%d%d%lf",&x,&y,&v);
        if(!pos[x] || !pos[y]) continue;
        e[++m]=edge(pos[x],pos[y],v);
        a[pos[y]]=min(a[pos[y]],v);
    }
    for(int i=2;i<=cnt;i++) ans+=a[i]*num[i];
    printf("%.2lf\n",zhuliu(1,n,m)+ans);
    return 0;
}