BZOJ.4145.[AMPPZ2014]The Prices(状压DP)

BZOJ


比较裸的状压DP。
刚开始写麻烦惹...
(f[i][s])表示考虑了前(i)家商店,所买物品状态为(s)的最小花费。
可以写求一遍一定去(i)商店的(f[i])(f[i][s]=f[i-1][s]+dis[i])),然后再和不去(i)商店的(f[i-1])取个(min)
复杂度是(O(nm2^m))吗...

可以优化,处理(f[s])表示在某家商店买(s)集合的物品的最小代价。然后令(g[s])表示考虑所有商店买(s)集合的最小代价,有(g[s]=min(f[s],g[s']+f[s ext{xor} s']))
复杂度(O(n2^m+3^m))


//27452kb	5284ms
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
#define lb(x) (x&-x)
#define gc() getchar()
typedef long long LL;
const int N=103,M=16;

int dis[N],cost[N][M+1],f[N][(1<<M)+1],ref[(1<<M)+1];

inline int read()
{
	int now=0;register char c=gc();
	for(;!isdigit(c);c=gc());
	for(;isdigit(c);now=now*10+c-48,c=gc());
	return now;
}

int main()
{
	int n=read(),m=read(),lim=(1<<m)-1;
	for(int i=0; i<m; ++i) ref[1<<i]=i;
	for(int i=1; i<=n; ++i)
		for(int j=(dis[i]=read(),0); j<m; ++j) cost[i][j]=read();
	memset(f,0x3f,sizeof f), f[0][0]=0;
	for(int i=1; i<=n; ++i)
	{
		for(int s=0; s<=lim; ++s) f[i][s]=f[i-1][s]+dis[i];
		for(int s=0; s<lim; ++s)
		{
			for(int ss=s^lim,j; ss; ss^=lb(ss))
			{
				j=ref[lb(ss)];
				f[i][s|(1<<j)]=std::min(f[i][s|(1<<j)],f[i][s]+cost[i][j]);
			}
		}
		for(int s=0; s<=lim; ++s) f[i][s]=std::min(f[i][s],f[i-1][s]);
	}
	printf("%d
",f[n][lim]);

	return 0;
}
原文地址:https://www.cnblogs.com/SovietPower/p/10752008.html