关于同余最短路

同余最短路可以解决形如「给定 (n) 个整数,求这 (n) 个整数能拼凑出多少的其他整数(可以重复取整数)」「给定 (n) 个整数,求这 (n) 个整数不能拼凑出的最小(最大)的整数」的问题 .

1. 跳楼机

给定 (x,y,z,h),问有多少个 (kin[1,h]) 使得 (k=px+qy+rz),其中 (p,q,r) 是变量 .

不妨设 (x<y<z) .

(d_i) 表示最小的满足条件且 (kequiv ipmod x)(k) .

考虑将 (q,r) 加一(将 (p) 加一对 (d_i) 没有贡献),可以发现:

[d_{(i+y)mod x}=d_i+y ]

[d_{(i+z)mod x}=d_i+z ]

发现这与最短路中的 (dis_v=dis_u+val_{u,v}) 类似,所以对于每个 (iin[0,x)) 建一个点,将 (i o (i+y)mod x) 连一条边权为 (y) 的边,将 (i o (i+z)mod x) 连一条边权为 (z) 的边,然后跑一边最短路即可求出所有 (d_i) .

不难发现,对于 (d_ile h),每个 (d_i) 其答案的贡献为 (leftlfloordfrac{h-d_i}x ight floor+1)(+1) 是因为自己也算一个),从而,答案为:

[sum_{i=0}^{x-1}left(left(leftlfloordfrac{h-d_i}x ight floor+1 ight)[d_ile h] ight) ]

时间复杂度 (O(x))(当然这里是 (x<y<z) 前提下的,一般的,时间复杂度是 (O(min{x,y,z}))

Code:

using namespace std;
const int N=1e5+500;
typedef unsigned long long ull;
typedef long long ll;
ll x,y,z,h,dis[N];
bool vis[N];
typedef vector<pair<ll,ll> > graph[N];
graph g;
inline void addedge(ll u,ll v,ll w){g[u].push_back(make_pair(v,w));}
struct node
{
	int idx; ll dis;
	node(int _=0,ll __=0){idx=_; dis=__;}
	bool operator <(const node& u)const{return dis>u.dis;}
};
void dijkstra(int s) // spfa 爬,dijkstra 天下第一
{
	memset(dis,0x3f,sizeof dis); memset(vis,false,sizeof vis);
	priority_queue<node> q; q.push(node(s,0)); dis[s]=0;
	while (!q.empty())
	{
		node now=q.top(); q.pop(); int u=now.idx,S=g[u].size();
		if (vis[u]) continue;
		vis[u]=true;
		for (int i=0;i<S;i++)
		{
			int v=g[u][i].first,w=g[u][i].second;
			if (dis[u]+w<dis[v]){dis[v]=dis[u]+w; if (!vis[v]) q.push(node(v,dis[v]));}
		}
	}
}
int main()
{
	scanf("%lld%lld%lld%lld",&h,&x,&y,&z);
	if ((x==1)||(y==1)||(z==1)){printf("%lld",h); return 0;} // 假设 x<y<z,若 x=1,则令 q=r=0,从而 k 可以取到 [1,h] 中所有数
	if (x>y) swap(x,y); // 这六句仅仅是为了让 x<y<z
	if (y>z) swap(y,z);
	if (x>z) swap(x,z);
	if (x>y) swap(x,y);
	if (y>z) swap(y,z);
	if (x>z) swap(x,z);
	for (int i=0;i<x;i++) addedge(i,(i+y)%x,y),addedge(i,(i+z)%x,z);
	dijkstra(0); ll ans=0;
	for (int i=0;i<x;i++)
		if (h>=dis[i]) ans+=(h-dis[i]-1)/x+1; // 注意这里初始楼层是 1 所以每个 d[i] 都要 +1
	printf("%lld",ans);
	return 0;
}

仅在此题给出代码 .

2. 墨墨的等式

给定 (a_{1cdots n},l,r),问有多少个 (kin[l,r]) 使得 (k=sumlimits_{i=1}^na_it_i),其中 (t_{1cdots n}) 是变量 .

上一题的拓展,上一题是 (l=1,n=3) 的特殊情况,这个一般的情况类似 .

时间复杂度 (O(nminlimits_i{a_i}))

3. 很多序列

给定 (n) 个递增正整数 (x_{1cdots n}),求不能由这些数字线性组合表示出的最大正整数 .
(1<n<6)(x_1le 10^{6-n})(x_2ge 10^{11+n})(x_nle 10^{12+n})(gcd(x_1,x_2)=1) .

如果 (n=2),那么就是小凯的疑惑,答案就是 (x_1x_2-x_1-x_2) .

类似的,把 (i)((i+x_i)mod x_1) 连边权为 (x_i) 的边,跑一边最短路,假设最短路是 (d_{1cdots n}) 答案就是

[max_{i}d_i-x_1 ]

4. 牛场围栏

给定序列 (l_{1cdots n}),每个 (l_i) 可以减少至多 (m),问无法表示为 (l_{1cdots n}) 的线性组合的最大数 .

和上一题类似 .

5. 货币系统

https://www.luogu.com.cn/problem/P5020

问题 (3) 的拓展,两种货币系统等价,当且仅当两种货币系统的 (d) 数组完全相同,开一个新数组记录转移用到的币值,并且要求币值尽量小 .

最后统计一下币值数目即可 .

6. 巡回

给定一个 (n) 个点 (m) 条边的无向图,问有没有 (1)(n) 长度为 (t) 的路径 .
(1le n,mle 50)(1le tle 10^{18})(1le 边权le 10^4)

任选一条边 ((u,n,w)),注意到,如果有一条长度为 (d) 的路径,也必然会有长度为 (d+2w) 的路径,故考虑在 (mod 2w) 意义下做同余最短路 .

(d_{u,k}) 表示满足存在一条从 (1)(u) 长度为 (i) 的路径且 (kequiv ipmod{2w}) 的最小的 (i),转移用最短路算法实现即可 .

答案就是 ([d_{n,Tmod 2w}le T]) .

注意若不存在从 (1)(n) 的路径,那么也是不存在这样的路径的 .

原文地址:https://www.cnblogs.com/CDOI-24374/p/14436668.html