[同余最短路]P2371

给定 $n,a_{1\dots n},l,r$，求出有多少 $b\in [l,r]$ 可以使 ${\sum }_{i=1}^n{a}_i{x}_i=b$ 存在非负整数解。
$n\le 12$，$0\le a_i\le 5\times 10^5$，$1\le l\le r\le 10^{12}$。

---

算是个思维题了（基本上没人会想到这是道最短路题）。首先建图，定义 $minn$ 为所有 $a_i$ 中最小的一个，同时将其他 $a$ 中元素前移一位。这里定义 $di{s}_i$ 是模 $minn$ 意义下的，同时 $di{s}_i$ 表示模 $minn$ 的所有符合条件的 $b$ 中最小的一个。容易得知在 $di{s}_i$ 之后的所有模 $minn$ 余 $i$ 的数都符合条件（直接加上若干个 $minn$ ）。那么我们可以得到在 $[0,r]$ 中答案的式子：
$$\sum _{i=0}^{minn-1}\frac{r-di{s}_i}{minn}+1(ifdi{s}_i\le r)$$
解释一下， $r-di{s}_i$ 表示剩下的数，除以循环节 $minn$ 得到之后还有多少符合条件的数（向下取整），再加上自己（ $+1$ ）就是答案。

再说如何建图，对于其他 $a$ 中元素而言，假设现在是第 $j$ 个数，它可以在 $a_i$ 个数之后到达 $(j+a_i)modminn$ ，所以可以建条边 $j->(j+a_i)modminn$ ，长度为 $a_i$ 。把所有边建起来就好了，然后再跑一遍最短路，起点是0，初始 $di{s}_0=0$ ，其他Inf。

处理 $[l,r]$ 区间转化成处理 $[0,l-1]$ 和 $[0,r]$ 两段，再相减即是答案。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
struct edge{
	int to,w,nxt;
}e[6000005];
ll a[25],n,l,r,minn,head[500005],ecnt;
ll dis[500005],bj[500005];
queue <int> Q;
void adde(int u,int v,int w){
	e[++ecnt].nxt=head[u];
	head[u]=ecnt;
	e[ecnt].to=v;
	e[ecnt].w=w;
}

void SPFA(){
	dis[0]=0;
	Q.push(0);
	while(!Q.empty()){
		int u=Q.front();
		Q.pop();
		bj[u]=0;
		for(int i=head[u];i;i=e[i].nxt){
			int v=e[i].to;
			if(dis[v]>dis[u]+e[i].w){
				dis[v]=dis[u]+e[i].w;
				if(bj[v]==0){
					bj[v]=1;
					Q.push(v);
				}
			}
		}
	}
}

ll ans(ll num){
	ll ret=0;
	for(int i=0;i<minn;++i){
		if(num>=dis[i])
			ret+=(num-dis[i])/minn+1;
	}
	return ret;
}

int main(){
	scanf("%lld %lld %lld",&n,&l,&r);
	scanf("%lld",&minn);
	for(int i=1;i<=n-1;++i){
		scanf("%lld",&a[i]);
		if(a[i]<minn){
			swap(a[i],minn);
		}
	}
	memset(dis,0x3f,sizeof(dis));
	for(int i=0;i<=minn-1;++i){
		for(int j=1;j<=n-1;++j){
			adde(i,(i+a[j])%minn,a[j]);
		}
	}
	SPFA();
	printf("%lld",ans(r)-ans(l-1));
	return 0;
}