51NOD 1962 区间计数 单调栈+二分 / 线段树+扫描线

 区间计数

 

基准时间限制：1.5 秒 空间限制：262144 KB 分值: 80

 

两个数列 {An} ， {Bn} ，请求出Ans, Ans定义如下：

Ans:=Σni=1Σnj=i[max{Ai,Ai+1,...,Aj}=max{Bi,Bi+1,...,Bj}]

注：[ ]内表达式为真，则为1，否则为0.

 

 

1≤N≤3.5×1051≤Ai,Bi≤N 

 

样例解释：

7个区间分别为：（1,4），（1,5），（2,4），（2,5），（3,3），（3,5），（4,5）

Input

第一行一个整数N
第二行N个整数Ai
第三行N个整数Bi

Output

一行，一个整数Ans

Input示例

5
1 4 2 3 4
3 2 2 4 1

Output示例

 

7

 

题解：

　　第一种做法是 两个单调栈 + 二分

　　两个数都同时维护一个单调递减的 栈

　　当元素出栈是我们可以确定以其为最大值所能掌管的一段区间，那么 在对应另一个栈内通过二分找到此值也能掌管一段区间

　　求出区间交就可以了

　　第二种做法比较类似

　　枚举固定最大值为x 的区间有哪些

　　在A数组中 假设x 所在位置为  i   ，利用单调栈可以求出 其掌管范围 [L1,R1]， 对应的

　　在B数组中 假设x 所在位置为  i   ，利用单调栈可以求出 其掌管范围 [L2,R2]，

　　那么抽出 [L1,i] [i,R1],  [L2,i] [i,R2], 在二维坐标系上就是两个矩阵，求出矩阵面积交就可以了

#include <bits/stdc++.h>
inline long long read(){long long x=,f=;char ch=getchar();while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}return x*f;}
using namespace std;
#define ls i<<1
#define rs ls | 1
#define mid ((ll+rr)>>1)
#define MP make_pair
typedef long long LL;
typedef unsigned long long ULL;
const long long INF = 1e18+1LL;
const double pi = acos(-1.0);
const int  N = 4e5 + , M = 1e3, inf = 2e9;
 
int n,a[N],b[N];
int l[],r[],q[][N],pos[][N];
LL cal(int x,int p,int i) {
     int ll = l[p], rr = r[p],ok = ll;
     while(ll <= rr) {
        if(q[p][mid] > x) {
            ll = mid + ;
        } else ok = mid,rr = mid - ;
    }
  //  cout<<i<<" "<<p<<" "<<ok<<" "<<q[p][ok]<<" "<<x<<endl;
    int mmp1,mmp2;
    mmp1 = max(pos[p][ok-]+,pos[p^][r[p^]-]+);
    mmp2 = i-;
 
    if(q[p][ok] == x)
        return 1LL * max(min(pos[p^][r[p^]],pos[p][ok]) - mmp1+,)
        * max(mmp2 - max(pos[p^][r[p^]],pos[p][ok]) + ,);
    else return ;
}
int main() {
    cin >> n;
    for(int i = ; i <= n; ++i) scanf("%d",&a[i]);
    for(int i = ; i <= n; ++i) scanf("%d",&b[i]);
    l[] = l[] = ;
    r[] = r[] = ;
    pos[][] = pos[][] = ;
    LL ans = ;
    for(int i = ; i <= n; ++i) {
        LL ret = ;
        while(l[] <= r[] && a[i] >= q[][r[]]) {
            ret += cal(q[][r[]],,i);
            r[]--;
        }
        q[][++r[]] = a[i];pos[][r[]] = i;
        while(l[] <= r[] && b[i] >= q[][r[]]) {
            ret += cal(q[][r[]],,i);
            r[]--;
        }
        q[][++r[]] = b[i];pos[][r[]] = i;
        //cout<<"fuck "<<i<<" "  << ret<<endl;
        ans += ret;
    }
    while(l[] <= r[]) {
        ans += cal(q[][r[]],,n+);
        r[]--;
    }
    cout<<ans<<endl;
    return ;
}