2019牛客暑期多校训练营(第二场)J Subarray

#2019牛客暑期多校训练营(第二场)J Subarray
题意: 长度为 1 e 9 1e9 1e9的区间 A A A下标为 [ 0 , 1 e 9 − 1 ] [0,1e9-1] [0,1e91],数输入 n n n个区间, [ l i , r i ] [l_i,r_i] [li,ri]区间类的值为1,其余为-1,问有多少区间和大于0.
题解: 看了来自大佬的博客,能够产生贡献的点最多只有 3 e 7 3e7 3e7个,意思是先求一个前缀和,然后画成图应该是这样。
[外链图片转存失败(img-ScePDY5G-1563932803584)(https://i.loli.net/2019/07/24/5d37b651d9d2043876.png)]
最差就是这样了,能够有影响的就只有这 3 e 7 3e7 3e7个点 (可能是分段的) ,那么问题来了,怎么求出这 3 e 7 3e7 3e7个点。
来自大佬的博客

为什么呢?在大佬眼里很简单,我画了个图才理解。

显然前面一个的 f [ i ] f[i] f[i]加后面 g [ i + 1 ] g[i+1] g[i+1]比两个区间之间的长度大就事连在一起的。(我果然太菜了)
然后处理完之后,就相当于处理

一个这样的前缀和,求所有位置有几个在他前面前缀比他小。如果范围小一点就用树状数组求一下就没了, 3 e 7 l o g ( 3 e 7 ) 3e7log(3e7) 3e7log(3e7)显然超时了。
看到这个前缀和,前后项最大只差了 1 1 1,上下界最大差值不超过 3 e 7 3e7 3e7,这再做个前缀和 sum 。用一个数组表示一个数字出现的次数,然后 s u m [ m ] = s u m [ m − 1 ] + b [ m ] sum[m]=sum[m-1]+b[m] sum[m]=sum[m1]+b[m],更新前缀和,答案就是 a n s + = s u m [ m − 1 ] ans+=sum[m-1] ans+=sum[m1]
中间有一些细节要处理

#include "bits/stdc++.h"

using namespace std;
typedef long long LL;
typedef unsigned long long uLL;
typedef pair<int, int> P;
typedef long double ld;

#define VNAME(value) (#value)
#define bug printf("*********\n");
#define debug(x) cout<<"["<<VNAME(x)<<" = "<<(x)<<"]"<<endl;
#define mid (l+r)/2
#define chl 2*k+1
#define chr 2*k+2
#define lson l,mid,chl
#define rson mid+1,r,chr
#define mem(a, b) memset(a,b,sizeof(a));

const long long mod = 1e9 + 7;
const int maxn = 1e6 + 5;
const int INF = 0x7fffffff;
const LL inf = 0x3f3f3f3f;
const double eps = 1e-8;

void fin() {
#ifndef ONLINE_JUDGE
    freopen("../data.in", "r", stdin);
#endif // ONLINE_JUDGE
}

#ifndef ONLINE_JUDGE
clock_t start = clock();
#endif // ONLINE_JUDGE
int n;
LL l[maxn], r[maxn];
LL f[maxn], g[maxn];
struct node {
    LL l, r, x;
} dat[maxn * 5];
LL sum[maxn * 30 + 20], b[maxn * 30 + 20];

int main() {
    fin();
    scanf("%d", &n);
    for (int i = 1; i <= n; i++) {
        scanf("%lld%lld", &l[i], &r[i]);
    }
    f[1] = r[1] - l[1] + 1;
    for (int i = 2; i <= n; i++) {
        f[i] = max(0LL, f[i - 1] - (l[i] - r[i - 1] - 1)) + r[i] - l[i] + 1;
    }
    g[n] = r[n] - l[n] + 1;
    for (int i = n - 1; i >= 1; --i) {
        g[i] = max(0LL, g[i + 1] - (l[i + 1] - r[i] - 1)) + r[i] - l[i] + 1;
    }
    int i = 1;
    LL ans = 0;

    while (i <= n) {
        int j = i + 1;
        LL mi = 0, mx = 0, pos = 0;
        //mi 下界,mx 上界
        while (j <= n && g[j] + f[j - 1] >= l[j] - r[j - 1] - 1) {
            j++;
        }
        j--;
        int t = i, num = 1;
        // [i,j] 区间是相互影响的
        for (; t <= j; t++) {  // 把每一段处理到 dat里面
            if (num == 1)dat[num].l = 0;
            else if (l[t] - r[t - 1] == 1) dat[num].l = pos + 1;
            else dat[num].l = pos;
            pos += r[t] - l[t] + 1;
            dat[num].r = pos;
            dat[num++].x = 1;
            mx = max(mx, pos);
            if (t != j) {
                dat[num].r = pos - 1;
                pos -= l[t + 1] - r[t] - 1;
                dat[num].l = pos + 1;
                dat[num++].x = 0;  // 0 表示下降 ,1 表示上升
                mi = min(pos, mi);
            } else {
                dat[num].r = pos - 1;
                dat[num].l = max(mi, pos - ((int) 1e9 - 1 - r[t]));
                dat[num++].x = 0;
            }
        }
        dat[0].r = min(mx, l[i]);
        dat[0].l = 1;
        dat[0].x = 0;
        for (int k = 0; k <= mx - mi + 200; k++)b[k] = sum[k] = 0;
        assert(mx - mi < maxn * 30);
        for (int k = 0; k < num; ++k) {
            dat[k].l += -mi;   //全部向上移动一个下届,保证最小值等于0
            dat[k].r += -mi;
            if (dat[k].x == 1) {
                for (int m = dat[k].l; m <= dat[k].r; ++m) {
                    b[m]++;
                    sum[m] = sum[m - 1] + b[m];
                    if (m >= 1)ans += sum[m - 1];
                }
            } else {
                LL tmp = 0;  
                if (dat[k].l > 0)tmp = sum[dat[k].l - 1];
                for (int m = dat[k].l; m <= dat[k].r; ++m) {
                    if (m >= 1)ans += tmp;
                    tmp = sum[m];  //如果是下降的 ,从小往上处理是不能把当前更新加进去因当前这个在他后面
                    b[m]++;
                    sum[m] = sum[m - 1] + b[m];
                }
            }
        }
        i = j + 1;
    }
    printf("%lld\n", ans);
#ifndef ONLINE_JUDGE
    cout << "RUNTIME:" << (1.0 * clock() - start) / 1000 << "ms" << endl;
#endif // ONLINE_JUDGE
    return 0;
}

给几组数据

2
5 9
12 17

3
3 5
6 6
7 9

3
3 5
8 8
9 10

1
999999998 999999999
200

76

49

4
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值