二分法

引子

给n个数，m次询问，每次询问给一个数a，找到n个数中比a小的最大的数，数据保证这样的数存在。

将n个数排序后，以小于a为条件，满足条件为1，不满足为0，则对n个数一一判断，我们将得到形如 11···1100···0 的序列，题目要求即找到最后一个1所在位置，那个位置上的数就是题目所要求的数。

体现在代码上，我们可以写一个check函数

bool check (int i)
{
    if (num[i] < a)
        return 1;
    else 
        return 0;
}

由于num数组有序，则num[ i ] 满足条件——小于a，num[ i - 1 ] 也必定满足条件。
如果n很大的话，通过枚举找最后一个1的过程就很慢，但我们可以通过二分来进行查找。

void work ()
{
    int ls = 1, rs = n + 1, mid;
    while (ls + 1 < rs)
    {
        mid = (ls + rs) / 2;
        if (check(mid))
            ls = mid;
        else
            rs = mid;
    } 
    cout << num[ls];
}

我们维护一个区间使答案必定在这个区间内，然后我们继续维护，使得左端点必定合法（check为1），右端点必定不合法（check为0)然后不断通过取左右端点的中点的方式来每次将这个区间缩短一半，这样当ls+1=rs的时候，左端点就是答案。其实就是一个左闭右开的区间。
例：
每个数的位置 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
是否满足条件 | 1 | 1 | 1 | 0 | 0 | 0 | 0 | 0 |
1、ls = 1； rs = 9；mid = 5
2、ls = 1；rs = 5；mid = 3
3、ls = 3；rs = 5；mid = 4；
此时，ls + 1 == rs；退出循环，ls 是要求的位置。

例题

对订单二分，寻找最小的无法满足的订单。

• 维护一个左闭右开的区间，以能够被满足为判断条件，ls 为左区间， rs 为右区间， mid = (ls + rs) / 2，每次判断第 mid 份订单是否能够被满足，
  1、能够被满足，则说明要找的数在区间[ mid，rs），则令ls = mid
  2、不能够被满足，则说明要找的数在区间[ ls，mid），则令rs = mid
  继续判断缩小区间，直到 ls + 1 == rs，此时ls是最后一个能被满足的订单，题目要求，第一个不能被满足的订单，即位置 ls + 1。

判断一个订单是否能被满足的check函数，需要根据前 x 个订单的租用天数和教室数，计算前 x 个订单中每天被租借出去的教室总数，并对每一天判断是否超过当天可租借的教室数量，若有一天超过则订单无法被满足，若全都没有超过，则订单可被满足。

• 可以采用差分法，计算前 x 个订单中每天被租借出去的教室总数。
  创建数组s记录每天被皆的教室总数，初始化为0，因此差分数列a的值全是0，对于一个订单，l ~ r天，被借出的教室数为d，令a[ l ] = a[ l ] + d，a [ r + 1] = a [ r + 1] - d，s [ i ] = s [ i - 1 ] + a [ i ] ，在 l ~ r 天内，s 都 多加了 a[ l ] 多加的d，而 r 天后，s 多加了 a[ l ] 加的d，但也减去了 a [ r + 1 ] 减去的d，从而实现了，对 l~r天，数组 s [l] ~ s[r]，都增加了 d。

差分法：
https://blog.csdn.net/qq_41054844/article/details/111592532?utm_medium=distribute.pc_relevant.none-task-blog-baidujs_baidulandingword-0&spm=1001.2101.3001.4242

# include <iostream>
# include <cstring>

using namespace std;

const int maxx = 1e6 + 5;
int N, m;
int r[maxx], d[maxx], s[maxx], t[maxx], j[maxx], a[maxx];

bool check(int x)
{
    memset (a, 0, sizeof(a));
    for (int i = 1; i <= x; ++ i )
    {
        a[s[i]] += d[i];
        a[t[i] + 1] -= d[i];
    }
    for (int i = 1; i <= N; ++ i)
        j[i] = j[i-1] + a[i];
    for (int i = 1; i <= N; ++ i)
        if (j[i] > r[i])
            return false;
    return true;
}

int main()
{
    cin >> N >> m;
    for (int i = 1; i <= N; ++ i)
        cin >> r[i];
    for (int i = 1; i <= m; ++ i)
        cin >> d[i] >> s[i] >> t[i];
    int ls = 0, rs = m + 1, mid;
    while(ls + 1 < rs)
    {
        mid = (ls + rs) / 2;
        if (check(mid))
            ls = mid;
        else 
            rs = mid;
    }
    if (ls == m)
        cout << "0" << endl;
    else 
    {
        cout << "-1" << endl;
        cout << ls + 1 << endl;
    }
    return 0;
}