洛谷 5070 「Ynoi 2015」即便看不到未来

发表于 | 分类于 |

挺妙的，学到了（

考虑离线处理询问，把询问按右端点排序。考虑从左往右加入数，同时计算每个数对所有区间的贡献。

设当前正在加入第 \(i\) 个数 \(a_i\)，则 \(a_i\) 只会影响左端点在 \(a_i\) 上一次出现的位置之后的区间。
则取出值域在 \([a_i-11,a_i+11]\) 的每个在 \(i\) 之前最后出现的数（每种权值找一个），显然只有这些数是可能产生或减少贡献的。
至于为什么要 \(\pm11\) 呢，是因为要求的是极长，故要考虑删除贡献。

把这些数按照位置从大到小排好序，然后扫过去，计算以 \(a_i-1\) 为结尾和 \(a_i+1\) 为开头的最长值域连续段的长度。
这样的话把这两个接起来便可得到新的连续段，同时删去这两段单独的贡献即可。

学了学 memset0 的代码。
题目略卡常。

代码： 

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
#include <cstdio>
#include <vector>
#include <utility>
#include <algorithm>
#define lowbit(x) ((x) & -(x))
using namespace std;

const int BUFF_SIZE = 1 << 20;
char BUFF[BUFF_SIZE],*BB,*BE;
#define gc() (BB == BE ? (BE = (BB = BUFF) + fread(BUFF,1,BUFF_SIZE,stdin),BB == BE ? EOF : *BB++) : *BB++)
template<class T>
inline void read(T &x)
{
    x = 0;
    char w = 0,ch = 0;
    for(;ch < '0' || ch > '9';w |= ch == '-',ch = gc());
    for(;ch >= '0' && ch <= '9';x = (x << 3) + (x << 1) + (ch ^ '0'),ch = gc());
    w && (x = -x,1);
}

const int N = 1e6;
const int W = 10;
int n,m,a[N + 5],ls[N + 5],vis[N + 5];
int c[N + 5][W + 5],ans[N + 5][W + 5];
vector< pair<int,int> > q[N + 5],vec;
inline void update(int u,int x,int k)
{
    for(;x;x -= lowbit(x))
        c[x][u] += k;
}
inline void query(int x,int *r)
{
    for(;x <= n;x += lowbit(x))
        for(register int u = 1;u <= W;++u)
            r[u] += c[x][u];
}
int main()
{
    read(n),read(m);
    for(register int i = 1;i <= n;++i)
        read(a[i]);
    int l,r;
    for(register int i = 1;i <= m;++i)
        read(l),read(r),q[r].push_back(make_pair(l,i));
    for(register int i = 1;i <= n;++i)
    {
        vec.clear();
        for(register int j = max(1,a[i] - 11);j <= min(a[i] + 11,N);++j)
            ls[j] > ls[a[i]] && (vec.push_back(make_pair(ls[j],j)),1);
        vec.push_back(make_pair(ls[a[i]],0)),vec.push_back(make_pair(i,a[i])),sort(vec.rbegin(),vec.rend());
        for(register int j = 0,L,R,U;j < vec.size() - 1;++j)
        {
            L = R = 0,vis[vec[j].second] = 1;
            for(;vis[a[i] - L - 1];++L);
            for(;vis[a[i] + R + 1];++R);
            U = L + R + 1;
            1 <= L && L <= 10 && (update(L,vec[j].first,-1),update(L,vec[j + 1].first, 1),1),
            1 <= R && R <= 10 && (update(R,vec[j].first,-1),update(R,vec[j + 1].first, 1),1),
            1 <= U && U <= 10 && (update(U,vec[j].first, 1),update(U,vec[j + 1].first,-1),1);
        }
        for(register int j = 0;j < vec.size() - 1;vis[vec[j].second] = 0,++j);
        for(register vector< pair<int,int> >::iterator it = q[i].begin();it != q[i].end();++it)
            query(it->first,ans[it->second]);
        ls[a[i]] = i;
    }
    for(register int i = 1;i <= m;++i,puts(""))
        for(register int j = 1;j <= W;++j)
            putchar('0' + ans[i][j] % 10);
}