【香蕉OI】月之美兔（最短路）

题意

有一条线上有 $n$ 个点，有 $m$ 个人，第 $i$ 个人站在 $a_i$ 位置，可以花一个代价向左或者向右跳 $f_i$ 的距离。

现在一个东西在 $1$ 号人手里，想要把这个东西交到 $2$ 号人手里。一个人可以把东西转交给另一个人当且仅当两个人在同一个位置。

求转交的最小代价。$n,m\le 2\ast 10^4$

思路

首先有一个显而易见的结论：代价最小的转交中，一个人最多只会拿到东西一次，也就是说不会送走了这个东西再拿回来。

那么每个人对转交的贡献就是从 $a_i$ 出发跑到一个位置 $p_1$ ，然后再跑到 $p_2$ 把这个东西交给下一个人。

所以现在已经有一个显而易见的暴力了，只需要把每两个位置 $p_1,p_2$ 之间连边，边权为某个人把东西从 $p_1$ 转运到 $p_2$ 的最小代价。然后跑最短路。 $O(n^2\ast m)$

然后考虑优化连边的过程，可以把 $f_i$ 和 $a_{i}mod{f}_i$ 相同的点一起建边，假设 $\forall f_i=F,\forall a_{i}mod{f}_i=R$。方法如下：

1. 把所有标号为 $k\ast F+R$ 的点复制一遍
2. 复制的点之间连上边权为 1 的边
3. 复制的点向原点连边权为 0 的边
4. 原点向复制的点连边权为“最近的满足 $f_i=F,a_{i}mod{f}_i=R$ 的点到他的距离”的边。

走到复制的点上就表示这个东西被一个人携带。（说的不是很清楚，感性理解一下就好了）

可以证明最后点数和边数都是 $O(n\sqrt{m})$ 的。因为想要复制一组点，就需要一组新的 $F$ 和 $R$，最开始复制 $n$ 个点最少只需要 $1$ 组 $F=1,R=0$ 就行了，而后面就依次需要 $2$ 组，$3$ 组 $...$。

所以跑最短路就好了。

代码

#include<bits/stdc++.h>
using namespace std;
#define pii pair<int, int>
#define mp make_pair
#define fi first
#define se second
const int inf = 1e9 + 7;
const int N = 20000 + 10, sqrN = 200, M = (N*sqrN)<<2;
int n, m, s, t, pt;
int ncnt, dis[N*sqrN], id[N];
struct node{
	int a, p, r;
}a[N];
int h[N*sqrN], ecnt, v[M], nxt[M], w[M];
priority_queue<pii, vector<pii>, greater<pii> > que;

template<class T>inline void read(T &x){
	x = 0; bool fl = 0; char c = getchar();
	while (!isdigit(c)){if (c == '-') fl = 1; c = getchar();}
	while (isdigit(c)){x = (x<<3)+(x<<1)+c-'0'; c = getchar();}
	if (fl) x = -x;
}

void _add(int x, int y, int z){
	nxt[++ecnt] = h[x]; v[ecnt] = y; w[ecnt] = z;
	h[x] = ecnt;
}

bool cmp(node x, node y){
	if (x.p != y.p) return x.p < y.p;
	else if (x.r != y.r) return x.r < y.r;
	else return x.a < y.a;
}

int dijkstra(){
	for (int i = 0; i <= ncnt; ++ i) dis[i] = inf;
	que.push(mp(0, s)); dis[s] = 0;
	while (!que.empty()){
		int d = que.top().fi, u = que.top().se;
		que.pop();
		if (d > dis[u]) continue;
		for (int i = h[u]; i; i = nxt[i])
			if (dis[v[i]] > dis[u]+w[i]){
				dis[v[i]] = dis[u]+w[i];
				que.push(mp(dis[v[i]], v[i]));
			}
	}
	int ans = inf;
	for (int i = 0; i < n; ++ i)
		if (abs(t-i)%pt == 0)
			ans = min(ans, dis[i]+abs(t-i)/pt);
	return ans;
}

int main()
{
	read(n); read(m);
	for (int i = 1; i <= m; ++ i){
		read(a[i].a), read(a[i].p);
		a[i].r = a[i].a%a[i].p;
	}
	s = a[1].a; t = a[2].a; pt = a[2].p;
	sort(a + 1, a + m + 1, cmp);
	ncnt = n-1; ecnt = 1;
	for (int i = 1, j; i <= m; i = j){
		int p = a[i].p, r = a[i].r;
		for (j = i; a[j].p == p && a[j].r == r && j <= m; ++ j);
		if (p == 0) continue;
		for (int k = r; k < n; k += p) dis[k] = inf, id[k] = ++ncnt;
		for (int k = i; k < j; ++ k){
			for (int o = a[k].a, tmp = 0; o < n; o += p, ++tmp){
				if (k+1 < j && o == a[k+1].a) break;
				dis[o] = min(dis[o], tmp);
			}
			for (int o = a[k].a, tmp = 0; o >= 0; o -= p, ++tmp){
				if (k-1 >= i && o == a[k-1].a) break;
				dis[o] = min(dis[o], tmp);
			}
		}
		for (int k = r; k < n; k += p){
			if (k+p < n){
				_add(id[k], id[k+p], 1); _add(id[k+p], id[k], 1);
			}
			_add(k, id[k], dis[k]); _add(id[k], k, 0);
		}
	}
	printf("%d\n", dijkstra());
	return 0;
}