bzoj2118 墨墨的等式 同余最短路

本文探讨了在大范围数值约束下,传统背包算法的局限性,并提出了一种新颖的解决方案。通过将问题转化为图的最短路径问题,利用Dijkstra算法在图中寻找从起点到每个节点的最短路径,从而有效解决了原问题。这种方法不仅避免了传统背包算法的复杂度瓶颈,还提供了一个更为通用的框架,适用于多种类似的问题场景。

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

https://www.lydsy.com/JudgeOnline/problem.php?id=2118

由于Bmi,Bmx太大,背包显然是不行的

我们考虑B对于a[1]取余,B=a[1]*x1+b(b=0,1,....a[1]-1)。

那么我们只要知道对于每一个b,最小体积是多少时能到达它就行了,之后就可以通过不断加a[1]来得到其他值。

那么我们对每一个b,连一条边到(b+a[1])%a[1],由于只有N<=12种数字,那么最多就N*A[1]条边,跑最短路即可。

#include<bits/stdc++.h>
using namespace std;

const int maxl=5e5+10;

int n,cnt;
int ehead[maxl];
long long bmi,bmx,ans;
long long a[maxl],dis[maxl];
bool vis[maxl];
struct ed
{
	int to,nxt,l;
}e[maxl*10*2];
typedef pair<long long,int> p;
priority_queue<p,vector<p>,greater<p> >q;

inline void add(int u,int v,int l)
{
	e[++cnt].to=v;e[cnt].l=l;
	e[cnt].nxt=ehead[u];ehead[u]=cnt;
}

inline void prework()
{
	scanf("%d%lld%lld",&n,&bmi,&bmx);
	for(int i=1;i<=n;i++)
		scanf("%lld",&a[i]);
	cnt=0;
	for(int i=0;i<a[1];i++)
		for(int j=2;j<=n;j++)
			add(i,(i+a[j])%a[1],a[j]);
}

inline void mainwork()
{
	memset(dis,0x3f,sizeof(long long)*a[1]);
	dis[0]=0;q.push(make_pair(0,0));
	p d;int u,v;
	while(!q.empty())
	{
		d=q.top();q.pop();
		u=d.second;
		if(d.first!=dis[u])
			continue;
		vis[u]=true;
		for(int i=ehead[u];i;i=e[i].nxt)
		{
			v=e[i].to;
			if(vis[v]) continue;
			if(dis[u]+e[i].l<dis[v])
			{
				dis[v]=dis[u]+e[i].l;
				q.push(make_pair(dis[v],v));
			}
		}
	}
	ans=0;
	for(int i=0;i<a[1];i++)
	if(dis[i]<bmi)
		ans+=(bmx-dis[i])/a[1]-(bmi-1-dis[i])/a[1];
	else if(dis[i]<=bmx)
		ans+=(bmx-dis[i])/a[1]+1;
}

inline void print()
{
	printf("%lld",ans);
}

int main()
{
	prework();
	mainwork();
	print();
	return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值