P2662 牛场围栏（同余最短路）

题目链接：点击这里

题目大意：
给定 $n$ 种木棍的长度 $a_i$，每种木棍数量无限，每根木棍可以削短 $1$ ~ $m$ 的长度再用，然后问用这些木棍不能组成的最大长度是多少，如果所有长度都能组成或不存在这个最大长度，那么输出 $-1$

题目分析：
先求出所有数字的 $gcd$ ，若其不为 $1$ ，此时由这些数构成的数字凑出来的数字一定是这个 $gcd$ 的倍数，一定是有空缺的，即此时为无解情况应输出 $-1$ ，无解情况和有解情况也可通过如下方法得到
此题为同余最短路求最大不能构成数字的模型，我们还是选定 $a_1$ 为 $base$ ，建图然后用最短路处理出 $dis$ 数组。如果存在 $di{s}_i=inf$ ，那么此时最大值可以表示成 $i+base\times inf$ ，表示无解
对于有解情况的答案就是 $max(di{s}_i)-a_1$

具体细节见代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<set>
#include<map>
#define ll long long
#define inf 0x3f3f3f3f
#define Inf 0x3f3f3f3f3f3f3f3f
#define int ll
using namespace std;
int read()
{
	int res = 0,flag = 1;
	char ch = getchar();
	while(ch<'0' || ch>'9')
	{
		if(ch == '-') flag = -1;
		ch = getchar();
	}
	while(ch>='0' && ch<='9')
	{
		res = (res<<3)+(res<<1)+(ch^48);//res*10+ch-'0';
		ch = getchar();
	}
	return res*flag;
}
const int maxn = 5e5+5;
const int maxm = 5005;
const int mod = 998244353;
const double pi = acos(-1);
const double eps = 1e-8;
struct Edge{
	int nxt,to,val;
}edge[maxn*13];
int n,m,a[maxn],cnt,head[maxn],dis[maxn];
bool vis[maxn];
void addedge(int from,int to,int val)
{
	edge[++cnt].nxt = head[from];
	edge[cnt].to = to;
	edge[cnt].val = val;
	head[from] = cnt;
}
struct node{
	int id,dis;
	bool operator < (const node &b) const {
		return dis > b.dis;
	}
};
void dijkstra(int s)
{
	memset(dis,0x3f,sizeof(dis));
	dis[s] = 0;
	priority_queue<node>qu;
	qu.push({s,0});
	while(!qu.empty())
	{
		int h = qu.top().id;qu.pop();
		if(vis[h]) continue;
		vis[h] = true;
		for(int i = head[h];i;i = edge[i].nxt)
		{
			int to = edge[i].to,val = edge[i].val;
			if(dis[to] > dis[h]+val)
			{
				dis[to] = dis[h]+val;
				qu.push({to,dis[to]});
			}
		}
	}
}
vector<int>v; 
signed main()
{
	n = read(),m = read();
	for(int i = 1;i <= n;i++) 
	{
		int x = read();
		for(int j = 0;j <= m && x-j > 0;j++) v.push_back(x-j);
	}
	sort(v.begin(),v.end());
	v.erase(unique(v.begin(),v.end()),v.end());
	if(v[0] == 1) return 0*puts("-1");
	for(int i = 0;i < v[0];i++)
		for(int j = 1;j < v.size();j++) addedge(i,(i+v[j])%v[0],v[j]);
	dijkstra(0);int ans = 0;
	for(int i = 0;i < v[0];i++) 
		if(dis[i] == inf) return 0*puts("-1");
		else ans = max(ans,dis[i]-v[0]);
	printf("%lld\n",ans);
	return 0;
}