微軟大樓設計方案（困難）

近日，微軟新大樓的設計方案正在廣泛征集中，其中一種方案格外引人註目。在這個方案中，大樓由 nn 棟樓組成，這些樓從左至右連成一排，編號依次為 11 到 nn，其中第 ii 棟樓有 h_ih?i?? 層。每棟樓的每一層為一個獨立的 辦公區域，可以步行 直達同層相鄰樓棟的辦公區域，以及 直達同樓棟相鄰樓層的辦公區域。

由於方案設計巧妙，上一層樓、下一層樓、向左右移動到相鄰樓棟同層的辦公區域均剛好需要 11 分鐘。在這些辦公區域中，有一些被 核心部門 占用了（一個辦公區域內最多只有一個核心部門），出於工作效率的考慮，微軟希望核心部門之間的移動時間越短越好。對於一個給定的 最大移動時間 kk，大樓的 協同值 定義為：有多少個 核心部門對 之間的移動時間不超過 kk。由於大樓門禁的限制，不可以走出整個大樓，也不可以登上天臺思考人生。你可以認為在辦公區域內的移動時間忽略不計，並且在大樓內總是按照最優方案進行移動。

對於一個給定的新大樓設計方案，你能算出方案的協同值麽？

輸入格式

第一行包含兩個正整數 n,k(1\leq k\leq 200020)n,k(1≤k≤200020)，分別表示大樓的棟數以及最大移動時間。

第二行包含 nn 個正整數 h_1,h_2,...,h_n(1\leq h_i\leq 20)h?1??,h?2??,...,h?n??(1≤h?i??≤20)，分別表示每棟樓的層數。

接下來一行包含一個正整數 mm，表示 核心部門 個數。

接下來 mm 行，每行兩個正整數 x_i,y_i(1\leq x_i\leq n,1\leq y_i\leq h_{x_i})x?i??,y?i??(1≤x?i??≤n,1≤y?i??≤h?x?i????)，表示該核心部門位於第 x_ix?i?? 棟樓的第 y_iy?i??層。

輸入數據保證 mm 個核心部門的位置不會重復。

對於簡單版本：1\leq n,m\leq 501≤n,m≤50；

對於中等版本：1\leq n\leq 200000,1\leq m\leq 20001≤n≤200000,1≤m≤2000；

對於困難版本：1\leq n,m\leq 2000001≤n,m≤200000。

輸出格式

輸出一個整數，即整個大樓的 協同值。

樣例解釋

樣例對應題目描述中的圖，核心部門 11 和核心部門 33 之間的距離為 8>78>7，因此不能計入答案。

樣例輸入

5 7
4 1 1 3 1
3
1 4
3 1
4 3

樣例輸出

2
分析：由於不同的高度僅有20個，所以考慮維護一個單調棧，將區間按到當前右端最小值分成h個部分；
　　　對於每個部分再枚舉高度；
　　　由k>=dis(A,B)=x?B??−x?A??+y?A??+y?B??-2min(h?x?A????,h?x?A??+1??,...,h?x?B????,y?A??,y?B??)算出xA範圍；
　　　最後前綴和加入答案；
　　　代碼寫的有點亂，從前往後遍歷點，同一棟樓的點多算了一次，要減去；
代碼：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <algorithm>
#include <climits>
#include <cstring>
#include <string>
#include <set>
#include <bitset>
#include <map>
#include <queue>
#include <stack>
#include <vector>
#include <cassert>
#include <ctime>
#define rep(i,m,n) for(i=m;i<=n;i++)
#define mod 1000000007
#define inf 0x3f3f3f3f
#define vi vector<int>
#define pb push_back
#define mp make_pair
#define fi first
#define se second
#define ll long long
#define pi acos(-1.0)
#define pii pair<int,int>
#define sys system("pause")
#define ls rt<<1
#define rs rt<<1|1
const int maxn=2e5+10;
const int N=5e2+10;
using namespace std;
ll gcd(ll p,ll q){return q==0?p:gcd(q,p%q);}
ll qpow(ll p,ll q){ll f=1;while(q){if(q&1)f=f*p%mod;p=p*p%mod;q>>=1;}return f;}
int n,m,k,t,a[maxn],b[maxn],id[maxn],h[maxn],qu[maxn],s[21][maxn];
ll ret;
bool cmp(int x,int y){return a[x]<a[y];}
int main()
{
    int i,j;
    scanf("%d%d",&n,&k);
    rep(i,1,n)scanf("%d",&h[i]);
    scanf("%d",&m);
    rep(i,1,m)scanf("%d%d",&a[i],&b[i]),id[i]=i,s[b[i]][a[i]]=1;
    rep(i,1,n)rep(j,1,20)rep(t,j,20)if(s[j][i]&&s[t][i]&&t-j<=k)ret--;
    rep(i,1,n)rep(j,1,20)s[j][i]+=s[j][i-1];
    sort(id+1,id+m+1,cmp);
    j=1;
    rep(i,1,m)
    {
        int pos=id[i];
        while(j<=a[pos])
        {
            while(qu[0]&&h[qu[qu[0]]]>=h[j])qu[0]--;
            qu[++qu[0]]=j;
            j++;
        }
        int x,y;
        rep(x,1,qu[0])
        {
            rep(y,1,20)
            {
                int l=(x-1>=1?qu[x-1]:0)+1,r=qu[x];
                l=max(l,a[pos]+b[pos]+y-2*min({y,b[pos],h[qu[x]]})-k);
                if(l<=r)
                {
                    ret+=s[y][r]-s[y][l-1];
                }
            }
        }
    }
    printf("%lld\n",ret);
    return 0;
}

微軟大樓設計方案