洛谷P2371 [国家集训队]墨墨的等式

不愧是国家集训队的题目，神仙题啊
题目所让求的实际上是$({\sum }_{i=1}^n{a}_i{x}_i)$在所给定的$[l,r]$区间内的取值个数
我们令$({\sum }_{i=1}^n{a}_i{x}_i)$为$S$
我们考虑对于$S$膜$a_i$的值，假设$S=k\ast a_i+rest$，如果存在方式可以表示出S的话，那么一定存在$S^{{}^{\prime }}=(k+1)\ast a_i+rest$，所以在这种情况下，S提供的贡献就是下式
$$ANS+=\{\begin{array}{ll}{B}_{max}-S)/a_i+1;& \text{S ∈[Bmin,Bmax] }\\ (B_{max}-S)/a_i-(B_{min}-S)/a_i& \text{S ∈[0,Bmin)}\end{array}$$

那么接下来我们需要考虑两件事

如何获取S

如何确定$a_i$的值

通过上式不难看出$a_i$是越小越好的，所以取$a$数组中最小值就好
那么S呢，因为对$a_i$取了模，所以对我们更重要的实际上是能达到$rest$的最小值是多少，因为取模后总共只有$a_i-1$种取值，比如膜$a_i$后是$rest$，那么$k_1\ast a_i+rest$和$k_2\ast a_i+rest$在$k_1<k_2$的情况下前者一定优于后者，前者对于答案的贡献一定包含后者，所以后者是无用的
所以我们要求的总共只有$a_i$种情况
那么这又如何求出呢，跑一次膜意义下最短路就行了
令$dis[i]$表示余数为$i$的最小数，那么连边很显然
对于$a_i(i>=2)$，我们连接$dis[i]$和$(dis[i]+a_i)\mathrm{m}\mathrm{o}\mathrm{d}{a}_1$，边权为$a_i$
完事，开跑

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<queue>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn = 5e5+7;
const ll INF = 1e13;
ll mod,dis[maxn],n,bmax,bmin,val[200];
bool cmp1(ll a,ll b){return a<b;}
int vis[maxn];

void spfa(){
	for(int i=0;i<mod;i++)dis[i]=INF;
	dis[0]=0;vis[0]=1;
	queue<ll>q;q.push(0);
	while(!q.empty()){
		ll f1=q.front();q.pop();vis[f1]=0;
		for(int i=2;i<=n;i++){
			ll to=(f1+val[i])%mod;
			if(dis[to]>dis[f1]+val[i]){
				dis[to]=dis[f1]+val[i];
				if(!vis[to])vis[to]=1,q.push(to);
			}
		}
	}
}

int main(){
	scanf("%lld%lld%lld",&n,&bmin,&bmax);bmin--;
	for(ll i=1;i<=n;i++){
		scanf("%lld",&val[i]);
	}
	sort(val+1,val+1+n,cmp1);
	mod=val[1];
	spfa();
	ll ans=0;
	for(int i=0;i<mod;i++){
		if(dis[i]>bmax)continue;
		else if(dis[i]>=bmin&&dis[i]<=bmax){
			ans+=(bmax-dis[i])/mod+1;
		}
		else {
			ans+=(bmax-dis[i])/mod-(bmin-dis[i])/mod;
		}
	}
	cout<<ans<<endl;
	return 0;
}