L3-028 森森旅游（单源最短路 + multiset 维护多重集最小值

L3-028 森森旅游 (30 分)
思路：
首先不考虑查询和修改操作，只考虑原始图中的答案
把路径分为前半段现金支付和后半段旅游金支付，可以发现去枚举一个决策点，维护 $min(cos{t}_1+cos{t}_2)$ 即可
前半段花费显然是单源最短路，后半段建反图也同样是单源最短路
由于后半段是旅游金，需要转化为现金，答案即为 $min(cos{t}_1+\frac{(cos{t}_2+a_i-1)}{a_i})$
这样我们就可以得到原始图中的最小花费
每次修改是将 $a[x]$ 修改为 $d$，每次需要插入删除一个元素，我们可以用 $multiset$ 维护，为什么是 $multi$，因为可能存在多个点的答案相同，而只修改其中一个点时使得这个点的花费更劣时，但答案不变。
注意，如果删除元素 $x$，$se.erase(x)$ 是将里边所有值为 $x$ 的元素删除，而只删除其中一个 $x$ 方法为 $se.erase(se.find(x))$
code：

#include<bits/stdc++.h>
#define endl '\n'
#define ll long long
#define ull unsigned long long
#define ld long double
#define all(x) x.begin(), x.end()
#define mem(x, d) memset(x, d, sizeof(x))
#define eps 1e-6
using namespace std;
const int maxn = 1e6 + 9;
const int mod = 1e9 + 7;
const int inf = 0x3f3f3f3f;
const ll INF = 0x3f3f3f3f3f3f3f3f;
ll n, m, q;
struct node{
	int next, to, w;
}e[maxn];
int head[2][maxn], cnt;
ll dis[2][maxn];
bool vis[2][maxn];
int a[maxn];

void add(int x, int y, int z, int op){
	e[++cnt].to = y;
	e[cnt].next = head[op][x];
	e[cnt].w = z;
	head[op][x] = cnt;
}
void dij(int op, int s){
	priority_queue<pair<ll,int>,vector<pair<ll,int> >,greater<pair<ll,int> > > q;
	q.push({0, s});
	dis[op][s] = 0;
	while(!q.empty()){
		int now = q.top().second;q.pop();
		if(vis[op][now]) continue;
		vis[op][now] = 1;
		for(int i = head[op][now]; i; i = e[i].next){
			int to = e[i].to;
			if(dis[op][to] > dis[op][now] + e[i].w){
				dis[op][to] = dis[op][now] + e[i].w;
				q.push({dis[op][to], to});
			}
		}
	}
}

void work()
{
	mem(dis, 0x3f);
	cin >> n >> m >> q;
	for(int i = 1; i <= m; ++i){
		int a, b, c, d;cin >> a >> b >> c >> d;
		add(a, b, c, 0);
		add(b, a, d, 1);
	}
	for(int i = 1; i <= n; ++i) cin >> a[i];
	dij(0, 1);
	dij(1, n);
	multiset <ll> se;
	for(int i = 1; i <= n; ++i){
		if(dis[0][i] == INF || dis[1][i] == INF) continue;
			se.insert(dis[0][i] + (dis[1][i] + a[i] - 1) / a[i]);
	}
	while(q--){
		ll x, d;cin >> x >> d;
		if(dis[0][x] == INF || dis[1][x] == INF) {
			cout << *se.begin() << endl;
			continue;
		}
		se.erase(se.find(dis[0][x] + (dis[1][x] + a[x] - 1ll) / a[x]));
		a[x] = d;
		se.insert(dis[0][x] + (dis[1][x] + a[x] - 1ll) / a[x]);
		cout << *se.begin() << endl;
	}
}

int main()
{
	ios::sync_with_stdio(0);
//	int TT;cin>>TT;while(TT--)
	work();
	return 0;
}