2018.09.27 bzoj2118: 墨墨的等式(最短路+背包)

传送门
好题啊。
首先找到最小的一个非零系数记做a1a_1a1,然后如果WWW modmodmod a1=W′a_1=W'a1=W modmodmod a1a_1a1,且WWW是方程的一个可行解,那么显然W′W'W也是一个可行解。
自然会想到我们用完全背包的思想对每一个余数求出要达到这个余数的最小可行解,这样整个范围中到达这个余数的可行解个数就可以统计出来了。
因此我们把所有BBBa1a_1a1的余数看做点,那么aaa数列就能够被当做边权,在建出的图上面跑最短路就行了。
代码:

#include<bits/stdc++.h>
#define ll long long
#define N 7500005
using namespace std;
int n,first[N],cnt=-1;
struct edge{int v,next;ll w;}e[N];
ll bmx,bmn,dis[N],a[20],mn=2e18;
bool in[N];
inline void add(int u,int v,ll w){e[++cnt].v=v,e[cnt].w=w,e[cnt].next=first[u],first[u]=cnt;}
inline ll query(ll x){
	ll ret=0;
	for(register int i=0;i<mn;++i)if(dis[i]<=x)ret+=(x-dis[i])/mn+1;
	return ret;
}
int main(){
	memset(first,-1,sizeof(first));
	scanf("%d%lld%lld",&n,&bmn,&bmx);
	for(register int i=1;i<=n;++i){
		scanf("%lld",&a[i]);
		if(!a[i]){--i,--n;continue;}
		mn=min(mn,a[i]);
	}
	for(register int i=0;i<mn;++i)for(register int j=1;j<=n;++j)add(i,(i+a[j])%mn,a[j]);
	queue<int>q;
	memset(dis,127,sizeof(dis));
	dis[0]=0,in[0]=1,q.push(0);
	while(!q.empty()){
		int x=q.front();
		q.pop(),in[x]=false;
		for(register int i=first[x];~i;i=e[i].next){
			int v=e[i].v;
			if(dis[v]>dis[x]+e[i].w){
				dis[v]=dis[x]+e[i].w;
				if(!in[v])in[v]=true,q.push(v);
			}
		}
	}
	printf("%lld",query(bmx)-query(bmn-1));
	return 0;
}

转载于:https://www.cnblogs.com/ldxcaicai/p/9738191.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值