PTA—森森快递

        题目：森森开了一家快递公司，叫森森快递。因为公司刚刚开张，所以业务路线很简单，可以认为是一条直线上的N个城市，这些城市从左到右依次从0到(N−1)编号。由于道路限制，第i号城市（i=0,⋯,N−2）与第(i+1)号城市中间往返的运输货物重量在同一时刻不能超过C**i公斤。

        公司开张后很快接到了Q张订单，其中j张订单描述了某些指定的货物要从S**j号城市运输到T**j号城市。这里我们简单地假设所有货物都有无限货源，森森会不定时地挑选其中一部分货物进行运输。安全起见，这些货物不会在中途卸货。

        为了让公司整体效益更佳，森森想知道如何安排订单的运输，能使得运输的货物重量最大且符合道路的限制？要注意的是，发货时间有可能是任何时刻，所以我们安排订单的时候，必须保证共用同一条道路的所有货车的总重量不超载。例如我们安排1号城市到4号城市以及2号城市到4号城市两张订单的运输，则这两张订单的运输同时受2-3以及3-4两条道路的限制，因为两张订单的货物可能会同时在这些道路上运输。

        分析：这是一道要用线段树或树状数组的题目。我选择了线段树。线段树可以用来解决要进行区间修改、查询， 单点修改、查询的题目。这道题每次选出给定区间内的最小值，并将该区间内的所有值减去这个最小值，最后输出最小值之和即可。

代码：

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int M = 1e6 + 10;
const int MAX = 0x3f3f3f3f;
int n, m ,a[M];
pair<int, int> v[M];
ll tree[M * 8], tag[M * 8];

void push_down(int p)//传递函数
{
    tag[p * 2] += tag[p];
	tag[p * 2 + 1] += tag[p];
	tree[p * 2] -= tag[p];
	tree[p * 2 + 1] -= tag[p];
	tag[p] = 0;
}

void bulid(int p, int pl, int pr) //表示建线段树
{
	if (pl == pr)
	{
		tree[p] = a[pl];
		return;
	}
	int mid = (pl + pr) / 2;
	bulid(p * 2, pl, mid);
	bulid(p * 2 + 1, mid + 1, pr);
	tree[p] = min(tree[p * 2], tree[p * 2 + 1]);
}

ll query(int L, int R, int p, int pl, int pr)//查询某个区间[L, R]内的最小值
{
	if (L <= pl && pr <= R)
		return tree[p];
	if (tag[p])
		push_down(p);
	int mid = (pl + pr) / 2;
	ll res = 1e18;
	if (mid >= L)
		res = min(res, query(L, R, p * 2, pl, mid));
	if (mid < R)
		res = min(res, query(L, R, p * 2 + 1, mid + 1, pr));
	return res;
}

void update(int L, int R, int p, int pl, int pr, ll d)//区间修改函数
{
	if (L <= pl && R >= pr)
	{
		tag[p] += d;
		tree[p] -= d;
		return;
	}
	if (tag[p])
		push_down(p);
	int mid = (pl + pr) / 2;
	if (L <= mid)
		update(L, R, p * 2, pl, mid, d);
	if (mid < R)
		update(L, R, p * 2 + 1, mid + 1, pr, d);
	tree[p] = min(tree[p * 2], tree[p * 2 + 1]);
}

int main()
{
	cin >> n >> m;
	for (int i = 1; i <= n - 1; i++)
		cin >> a[i];
	bulid(1, 1, n - 1);
	for (int i = 0; i < m; i++)
	{
		int l, r;
		cin >> l >> r;
		if (l > r)
			swap(l, r);
		v[i] = {r, l};
	}
	sort(v, v + m);
	ll ans = 0;
	for (int i = 0; i < m; i++)
	{
		ll sum = query(v[i].second + 1, v[i].first, 1, 1, n - 1);
		ans += sum;
		update(v[i].second + 1, v[i].first, 1, 1, n - 1, sum);
	} 
	cout << ans;
	return 0;
}