BZOJ 4145: [AMPPZ2014]The Prices

本文介绍了一个具体的状压DP问题:给定若干商店和物品种类,通过优化算法寻找购买指定物品组合所需的最小总费用。文章提供了两种不同的实现思路,并讨论了它们的时间复杂度。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

Description

你要购买m种物品各一件,一共有n家商店,你到第i家商店的路费为d[i],在第i家商店购买第j种物品的费用为c[i][j],

求最小总费用。

 

Input

第一行包含两个正整数n,m(1<=n<=100,1<=m<=16),表示商店数和物品数。

接下来n行,每行第一个正整数d[i](1<=d[i]<=1000000)表示到第i家商店的路费,接下来m个正整数,

依次表示c[i][j](1<=c[i][j]<=1000000)。

 

Output

一个正整数,即最小总费用。

 

Sample Input

3 4
5 7 3 7 9
2 1 20 3 2
8 1 20 1 1

Sample Output

16

HINT

 

在第一家店买2号物品,在第二家店买剩下的物品。

 

从m=16联想到状压dp

最暴力的写法时4^n*n

#include<cstdio> 
#include<iostream> 
#define ll long long 
using namespace std; 
inline int read() 
{ 
    char ch=getchar(); 
    int ret=0; 
    while(ch<'0'||ch>'9') ch=getchar(); 
    while(ch>='0'&&ch<='9') 
        ret=(ret<<1)+(ret<<3)+ch-'0',ch=getchar(); 
    return ret; 
} 
  
int n,m; 
const int N=105,M=(1<<16)+5; 
int a[20],d[N],lg[M]; 
ll s[N][M],f[M]; 
  
int main() 
{ 
    n=read(),m=read(); 
    lg[0]=-1; 
    for(int i=1;i<=M;i++) lg[i]=lg[i>>1]+1; 
    for(int i=1;i<=n;i++) 
    { 
        d[i]=read(); 
        for(int j=1;j<=m;j++) a[j]=read(); 
        for(int j=1;j<(1<<m);j++) 
            s[i][j]=s[i][j-(j&-j)]+a[lg[j&-j]+1]; 
    }    
    for(int i=0;i<(1<<m);i++) f[i]=1e15; 
    f[0]=0; 
    for(int i=1;i<=n;i++) 
        for(int S=0;S<(1<<m);S++) 
            for(int j=1;j<(1<<m);j++) 
            { 
                if(S&j) continue; 
                f[S|j]=min(f[S]+(ll)d[i]+s[i][j],f[S|j]); 
            } 
    printf("%lld\n",f[(1<<m)-1]); 
    return 0; 
} 

可惜超时,只有16分

再把他优化成子集枚举

恭喜,超时,变成32分

3^n*n

如果不改变状态只能干掉一个n或者把3^n优化成2^n*n

菜鸡先写3^n

对于每个状态S

求出对于每个状态所有n的最小价值{cost(i,S)+d[i]}

再子集穷举

如果存在有2堆东西从同一家商店购买,则必有一种转移会把两种东西合并,肯定比当前更优

#include<cstdio>
#include<iostream>
#define ll long long
using namespace std;
inline int read()
{
	int ret=0;
	char ch=getchar();
	while(ch<'0'||ch>'9') ch=getchar();
	while(ch>='0'&&ch<='9')
		ret=(ret<<1)+(ret<<3)+ch-'0',ch=getchar();
	return ret;
}

int n,m;
const int N=105,M=1000000;
int a[20];
ll s[M],f[M];

int main()
{
	n=read(),m=read();
	for(int S=0;S<(1<<m);S++) s[S]=1e18;
	for(int i=1;i<=n;i++)
	{
		int d=read();
		for(int j=1;j<=m;j++) a[j]=read();
		for(int S=1;S<(1<<m);S++)
		{
			ll ss=d;
			int x=S,j=0;
			while(x)
			{	
				j++;
				if(x&1) ss+=a[j];
				x>>=1;
			}
			s[S]=min(s[S],ss);
		}
	}
	for(int S1=1;S1<(1<<m);S1++)
	{
		f[S1]=s[S1];
		for(int S2=S1&(S1-1);S2;S2=(S2-1)&S1)
			f[S1]=min(f[S1],f[S2]+s[S1-S2]);
	}
	printf("%lld\n",f[(1<<m)-1]);
	return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值