【刷题】BZOJ 2260 商店购物

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<Ci<=1000)和一个整数Mi (0<=Mi<=100)组成，其间由一个空格分隔，分别表示第i种商品的原价和所需数量。第(n+2)行又是一个整数k，表示Grant小店的优惠方案总数。接着k行，每行有二个整数A，B（1<=A,B<=n）和一个实数P（ 0<=P<1000），表示一种优惠方案，即如果您购买了商品A，您就可以以P元/件的优惠价格购买商品B，P小于商品B的原价。所有优惠方案的（A,B）都是不同的。为了方便，Grant不收分币，所以所有价格都不会出现分。

Output

只有一个实数，表示最少需要花多少钱。输出实数须保留两位小数。

Sample Input

4
10.00 1
1.80 1
3.00 0
2.50 2
2
1 4 2.00
4 2 1.50

Sample Output

15.50

Solution

这不是双倍经验？！
与【刷题】BZOJ 4349 最小树形图一模一样

#include<bits/stdc++.h>
#define ui unsigned int
#define ll long long
#define db double
#define ld long double
#define ull unsigned long long
const int MAXN=50+5,MAXM=MAXN*MAXN;
const db inf=100000000000.00;
int n,m,times[MAXN],vis[MAXN],bel[MAXN],snt,s,pre[MAXN],nd,M[MAXN];
db in[MAXN],G[MAXN][MAXN],ans;
struct node{
	int u,v;
	db k;
};
node side[MAXM];
template<typename T> inline void read(T &x)
{
	T data=0,w=1;
	char ch=0;
	while(ch!='-'&&(ch<'0'||ch>'9'))ch=getchar();
	if(ch=='-')w=-1,ch=getchar();
	while(ch>='0'&&ch<='9')data=((T)data<<3)+((T)data<<1)+(ch^'0'),ch=getchar();
	x=data*w;
}
template<typename T> inline void write(T x,char ch='\0')
{
	if(x<0)putchar('-'),x=-x;
	if(x>9)write(x/10);
	putchar(x%10+'0');
	if(ch!='\0')putchar(ch);
}
template<typename T> inline void chkmin(T &x,T y){x=(y<x?y:x);}
template<typename T> inline void chkmax(T &x,T y){x=(y>x?y:x);}
template<typename T> inline T min(T x,T y){return x<y?x:y;}
template<typename T> inline T max(T x,T y){return x>y?x:y;}
inline db solve(int rt,int n)
{
	db res=0;
	while(true)
	{
		for(register int i=1;i<=n;++i)in[i]=inf;
		for(register int i=1;i<=snt;++i)
			if(side[i].u!=side[i].v&&in[side[i].v]>side[i].k)in[side[i].v]=side[i].k,pre[side[i].v]=side[i].u;
		for(register int i=1;i<=n;++i)
			if(i!=rt&&in[i]==inf)return -1;
		int cnt=0;
		memset(bel,0,sizeof(bel));
		memset(vis,0,sizeof(vis));
		in[rt]=0;
		for(register int i=1,j;i<=n;++i)
		{
			res+=in[i];j=i;
			while(j!=rt&&vis[j]!=i&&!bel[j])vis[j]=i,j=pre[j];
			if(j!=rt&&!bel[j])
			{
				bel[j]=++cnt;
				for(register int k=pre[j];k!=j;k=pre[k])bel[k]=cnt;
			}
		}
		if(!cnt)break;
		for(register int i=1;i<=n;++i)
			if(!bel[i])bel[i]=++cnt;
		for(register int i=1,u,v;i<=snt;++i)
		{
			u=side[i].u,v=side[i].v;
			side[i].u=bel[u];side[i].v=bel[v];
			if(bel[u]^bel[v])side[i].k-=in[v];
		}
		n=cnt;
		rt=bel[rt];
	}
	return res;
}
int main()
{
	read(n);
	s=++nd;
	for(register int i=1;i<=n;++i)
	{
		in[i]=inf;
		db cost;scanf("%lf",&cost);read(times[i]);
		if(times[i])
		{
			M[i]=++nd;
			side[++snt]=(node){s,M[i],cost};
			chkmin(in[i],cost);
		}
	}
	read(m);
	for(register int i=1;i<=m;++i)
	{
		int x,y;read(x);read(y);
		db cost;scanf("%lf",&cost);
		if(M[x]&&M[y])
		{
			side[++snt]=(node){M[x],M[y],cost};
			chkmin(in[y],cost);
		}
	}
	for(register int i=1;i<=n;++i)
		if(times[i]>1)ans+=(times[i]-1)*in[i];
	printf("%.2f\n",ans+solve(s,nd));
	return 0;
}