codeforces 650D. Zip-line 线段树

题目链接

题目的意思很简单， 就是给你n个数， m个询问， 每次询问修改某一个位置的值， 然后问你修改完之后数列的lis是多少。 询问独立。

对于原数列， 我们将它离散化， 令dp1[i]为以i为结尾位置的最长上升子序列的长度， dp[2]为以i结尾的从后往前的最长下降子序列的长度。

原数列的lis显然为max(dp1[i]+dp2[i]-1)。

然后我们求出哪些位置是关键位置， 所谓关键位置， 就是说如果把这个位置的值改变， 那么lis的值也许就会减1。  求关键位置的方法看代码。

 然后对于每个询问， 令x[i]为位置， y[i]为修改后并离散化的值，我们令dp3[i]表示将x[i]位置修改为y[i]之后， 以x[i]结尾的最长上升子序列的长度， dp4[i]为最长下降的长度。

那么对于每次询问， 如果x[i]是关键节点， ans[i] = max(lis-1, dp3[i]+dp4[i]-1), 否则的话， ans[i] = max(lis, dp3[i]+dp4[i]-1)。

具体的可以看代码。

还有一点要注意的是数组不能开成4e5， 要开成8e5。

#include <iostream>
#include <vector>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <map>
#include <set>
#include <string>
#include <queue>
#include <stack>
#include <bitset>
using namespace std;
#define pb(x) push_back(x)
#define ll long long
#define mk(x, y) make_pair(x, y)
#define lson l, m, rt<<1
#define mem(a) memset(a, 0, sizeof(a))
#define rson m+1, r, rt<<1|1
#define mem1(a) memset(a, -1, sizeof(a))
#define mem2(a) memset(a, 0x3f, sizeof(a))
#define rep(i, n, a) for(int i = a; i<n; i++)
#define fi first
#define se second
typedef pair<int, int> pll;
const double PI = acos(-1.0);
const double eps = 1e-8;
const int mod = 1e9+7;
const int inf = 1061109567;
const int dir[][2] = { {-1, 0}, {1, 0}, {0, -1}, {0, 1} };
const int maxn = 8e5+5;
int sum[maxn<<2], num[maxn], dp1[maxn], dp2[maxn], dp3[maxn], dp4[maxn], a[maxn], b[maxn];
int x[maxn], y[maxn], ans[maxn];
vector <pll> v[maxn];
void update(int p, int val, int l, int r, int rt) {
    if(l == r) {
        sum[rt] = max(sum[rt], val);
        return ;
    }
    int m = l+r>>1;
    if(p<=m)
        update(p, val, lson);
    else
        update(p, val, rson);
    sum[rt] = max(sum[rt<<1],sum[rt<<1|1]);
}
int query(int L, int R, int l, int r, int rt) {
    if(L>R)
        return 0;
    if(L<=l&&R>=r) {
        return sum[rt];
    }
    int m = l+r>>1, ret = 0;
    if(L<=m)
        ret = query(L, R, lson);
    if(R>m)
        ret = max(ret, query(L, R, rson));
    return ret;
}
int main()
{
    int n, m, cnt = 0;
    cin>>n>>m;
    for(int i = 1; i<=n; i++) {
        scanf("%d", &a[i]);
        b[cnt++] = a[i];
    }
    for(int i = 0; i<m; i++) {
        scanf("%d%d", &x[i], &y[i]);
        b[cnt++] = y[i];
        v[x[i]].pb(mk(i, y[i]));
    }
    sort(b, b+cnt);
    cnt = unique(b, b+cnt)-b;
    for(int i = 1; i<=n; i++) {
        a[i] = lower_bound(b, b+cnt, a[i])-b+1;
    }
    for(int i = 1; i<=n; i++) {
        dp1[i] = query(1, a[i]-1, 1, cnt, 1)+1;
        for(int j = 0; j<v[i].size(); j++) {
            int tmp = v[i][j].fi;
            int tmpy = lower_bound(b, b+cnt, v[i][j].se)-b+1;
            dp3[tmp] = query(1, tmpy-1, 1, cnt, 1)+1;
        }
        update(a[i], dp1[i], 1, cnt, 1);
    }
    mem(sum);
    for(int i = n; i>=1; i--) {
        dp2[i] = query(a[i]+1, cnt, 1, cnt, 1)+1;
        for(int j = 0; j<v[i].size(); j++) {
            int tmp = v[i][j].fi;
            int tmpy = lower_bound(b, b+cnt, v[i][j].se)-b+1;
            dp4[tmp] = query(tmpy+1, cnt, 1, cnt, 1)+1;
        }
        update(a[i], dp2[i], 1, cnt, 1);
    }
    int maxx = 0;
    for(int i = 1; i<=n; i++) {
        maxx = max(dp1[i]+dp2[i]-1, maxx);      //求原数列lis
    }
    for(int i = 1; i<=n; i++) {
        if(dp1[i]+dp2[i]-1 == maxx) {
            num[dp1[i]]++;              //如果num[dp1[i]] == 1, 那么它就是关键节点
        }
    }
    for(int i = 0; i<m; i++) {
        int tmp = maxx;
        if(dp1[x[i]]+dp2[x[i]]-1 == maxx && num[dp1[x[i]]] == 1) {
            tmp--;
        }
        tmp = max(tmp, dp3[i]+dp4[i]-1);
        ans[i] = tmp;
    }
    for(int i = 0; i<m; i++) {
        printf("%d
", ans[i]);
    }
    return 0;
}