hdu 3810 超大容量背包(n=50不是40)

先贴一个挑战上解决n=40时的超大容量背包,不过解决n=50时会TLE和MLE

#include<cstdio>
#include<cstring>
#include<vector>
#include<algorithm>
#define MIN(x,y) ((x)>(y)?(y):(x))
using namespace std;
const long long INF=5010000000000;
vector<int>Q[100];
pair<int,int>ps[1<<21];
int map[60][60];
int vis[60];
long long int w[60],v[60];
int cnt,n,tot;
long long m; 
void dfs(int x)
{
	vis[x]=1;
	Q[cnt].push_back(x);
	for(int i=1;i<=n;i++)
	{
		if(map[x][i]&&!vis[i])
		dfs(i);
	}
}
void divid()
{
	cnt=0;
	memset(vis,0,sizeof(vis));
	for(int i=1;i<=n;i++)
    {
    	if(!vis[i])
    	{
    		Q[cnt].clear();
    		dfs(i);
    		cnt++;
		}
	}
}
int find(int x)
{
	int l=0,r=tot;
	while(l<r)
	{
		int mid=(l+r)>>1;
		if(ps[mid].first+x<m)
		l=mid+1;
		else
		r=mid;
	}
	return l;
}
long long  solve(int x)
{
	int n2=Q[x].size()/2;
	if(n2>21)
	n2=21; 
	for(int i=0;i<1<<n2;i++)
	{
		long long int sw=0;
		long long int sv=0;
		for(int j=0;j<n2;j++)
		{
			if(i>>j&1)
			{
				sw+=w[Q[x][j]];
				sv+=v[Q[x][j]];
			}
		}
		ps[i]=make_pair(sw,sv);
	}
	sort(ps,ps+(1<<n2));
	tot=1;
	for(int i=1;i<1<<n2;i++)
	{
		if(ps[tot-1].second<ps[i].second)
		ps[tot++]=ps[i];
	}
//	printf("n2=%d tot=%d\n",n2,tot);
//	for(int i=0;i<tot;i++)
//	printf("first=%d second=%d\n",ps[i].first,ps[i].second); 
	long long int res=INF;
	for(int i=0;i<1<<(Q[x].size()-n2);i++)
	{
		long long int sw=0,sv=0;
		for(int j=0;j<Q[x].size()-n2;j++)
		{
			if(i>>j&1)
			{
				sw+=w[Q[x][n2+j]];
				sv+=v[Q[x][n2+j]];
			}
		}
		if(sv<res)
		{
			int lo=find(sw);
			if(lo!=tot) 
			res=MIN(res,sv+ps[lo].second);
		}
	}
	return res;
} 
long long int cal()
{
	long long int res,t;
	res=INF;
	for(int i=0;i<cnt;i++)
	{
		t=solve(i);
		if(t<res)
		res=t;
	}
	return res;
}
int main()
{
	int T,z=1;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d%lld",&n,&m);
//		printf("n=%d m=%d \n ",n,m); 
		memset(map,0,sizeof(map));
		for(int i=1;i<=n;i++)
		{
			int t;
			scanf("%lld%lld%d",&v[i],&w[i],&t);
			while(t--)
			{
				int x;
				scanf("%d",&x);
				map[i][x]=map[x][i]=1;
			}
		}
//		for(int i=1;i<=n;i++)
//		printf("v=%d w=%d\n",v[i],w[i]); 
//		for(int i=1;i<=n;i++)
//		{
//			for(int j=1;j<=n;j++)
//			printf("%d ",map[i][j]);
//			printf("\n");
//		} 
		divid(); 
//		printf("cnt=%d\n",cnt);
//		for(int i=0;i<cnt;i++)
///		{
//			printf("size=%d\n",Q[i].size()); 
//			for(int j=0;j<Q[i].size();j++)
//			printf("%d ",Q[i][j]);
//			printf("\n\n");
//		}
		long long int res=cal();
		if(res==INF)
		printf("Case %d: Poor Magina, you can't save the world all the time!\n",z++);
		else
		printf("Case %d: %lld\n",z++,res); 
	}
}


下面是用队列模拟。。。

#include<cstdio>
#include<cstring>
#include<queue>
#include<vector>
#define INF 500000010
#define MIN(x,y) ((x)>(y)?(y):(x))
using namespace std;
int n,cnt,vis[60],map[60][60];
long long m,v[60],w[60];
vector<int>Q[60];
struct node
{
	long long v,w;
	bool operator<(const node a)const
	{
		if(w!=a.w)
		return w<a.w;
		else
		return v>a.v;
	}
};
void dfs(int x)
{
	vis[x]=1;
	Q[cnt].push_back(x);
	for(int i=1;i<=n;i++)
	if(map[x][i]&&!vis[i])
	dfs(i);
}
void init()
{
	cnt=0;
	memset(vis,0,sizeof(vis));
	for(int i=1;i<=n;i++)
	if(!vis[i])
	{
		Q[cnt].clear();
		dfs(i);
		cnt++;
	}
}
long long cal()
{
	long long res=INF;
	priority_queue<node>q1,q2;
	node s; 
	for(int i=0;i<cnt;i++)
	{
		while(!q1.empty())
		q1.pop();
		while(!q2.empty())
		q2.pop();
		s.v=s.w=0; 
		q1.push(s);
		for(int j=0;j<Q[i].size();j++)
		{
			while(!q1.empty())
			{
				s=q1.top();
				q1.pop();
				q2.push(s);
				s.v+=v[Q[i][j]];
				s.w+=w[Q[i][j]];
				if(s.w>=m)
				{
					res=MIN(res,s.v);
					continue;
				}
				if(s.v>=res)
				continue;
				q2.push(s);
			}
			long long temp=INF;
			while(!q2.empty())
			{
				s=q2.top();
				q2.pop();
				if(s.v<temp)
				{
					temp=s.v;
					q1.push(s);
				}
			}
		}
	}
	 return res;
}
int main()
{
	int T,lo=1;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d%lld",&n,&m);
		memset(map,0,sizeof(map));
		for(int i=1;i<=n;i++)
		{
			int t;
			scanf("%lld%lld%d",&v[i],&w[i],&t);
			while(t--)
			{
				int x;
				scanf("%d",&x);
				map[i][x]=map[x][i]=1; 
			}
		}
		init();
		long long res=cal();
		if(res==INF)
		printf("Case %d: Poor Magina, you can't save the world all the time!\n",lo++);
		else
		printf("Case %d: %d\n",lo++,res);
	} 
}





  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值