\(O(10 n \log n)\) 能过,居然不卡常,青结了。
感觉是比较套路的一道 Ynoi 了 qwq。
首先看题目,需要找的就是一段长度为 \(1 \dots 10\) 的极长连续的(公差为 1)的等差数列,考虑暴力把一个个数丢进去,会造成怎样的结果,无非这两种情况:
-
无法拓展,单独成段;
-
连接上了在自己左右两边的极长连续段(记为 \(S_1\) 和 \(S_2\)),那么新的极长连续段为 \(|S_1|+|S_2|+1\),同时原来的两个对答案的贡献消失。
把询问离线下来,都丢到右端点上,对于第二种情况,相当于去掉在 \(S_1\) \(S_2\) 上的覆盖,同时加上 \(S_1+mid+S_2\) 的覆盖,对于询问,实际上就是 \(l\) 这个点被某一个长度为 \(x_i\) 的线段覆盖的次数,注意到段长只有 \(10\),直接暴力分成 \(10\) 个线段做事没问题的,区间修改,单点询问肯定用常数小的树状数组啦。
至于怎么在暴力拓展时保留单点询问的正确性,对于一个新数 \(x\),设 \(last_x\) 为上一次出现的次数,则 \(x\) 只对 \((last_x, i]\) 产生贡献(前面的去重后选用 \(last_x\) 的就行了),于是我们从 \(i\) 往左拓展,每次遇到上面的情况就进行区间覆盖操作,注意若 10 能拓展到 11 不算数。
复杂度 \(O(10nk \log n)\),\(k\) 是一个小常数。
Code:
#include <stdio.h>
#include <algorithm>
#include <string.h>
#include <cctype>
#include <vector>
using namespace std;
const int N=1e6+1;
inline int read(){
char ch=getchar();int x=0, f=1;
while(!isdigit(ch)) f=(ch=='-'?-1:f), ch=getchar();
while(isdigit(ch)) x=(x<<3)+(x<<1)+ch-'0', ch=getchar();
return x*f;
}
typedef pair <int, int> Pii;
int n, m, a[N], las[N+11], Mx;
char ans[N][11];
Pii tmp[25];
vector < Pii > Q[N];
struct BITtree{
int sum[N];
#define lowbit(x) ((x)&(-(x)))
void add(int x, int v){for(; x; x-=lowbit(x)) sum[x]+=v;}
int query(int x){int res=0;for(; x<=n; x+=lowbit(x)) res+=sum[x];return res;}
#undef lowbit
}bit[10];
void chafen(int flg, int l, int r, int v){
if(flg>10||flg<1) return ;
bit[flg-1].add(l, v), bit[flg-1].add(r, -v);
}
signed main(){
n=read(), m=read();
for(int i=1; i<=n; i++) Mx=max(Mx, a[i]=read());
for(int i=1, l, r; i<=m; i++)
l=read(), r=read(), Q[r].push_back(make_pair(l, i));
int pos=1;
for(int i=1; i<=n; i++){
int x=a[i], L=max(1, x-11), R=min(x+11, Mx), top=0;
for(int j=L; j<=R; j++) if(las[j]>las[x]) tmp[++top]=make_pair(-las[j], j);
tmp[++top]=make_pair(-i, x);sort(tmp+1, tmp+top+1);
for(int j=1; j<=top; j++) tmp[j].first*=-1;tmp[top+1].first=las[x], L=R=0;//找到所有的极长区间
for(int j=1; j<=top; j++){
int nxt=tmp[j].second, least=tmp[j].first;
if(nxt<x) while(L<=10&&(x-L-1==nxt||las[x-L-1]>least)) L++;
if(nxt>x) while(R<=10&&(x+R+1==nxt||las[x+R+1]>least)) R++;
//找到匹配的最长段
chafen(L, least, tmp[j+1].first, -1);
chafen(R, least, tmp[j+1].first, -1);
//消掉之前的贡献
chafen(L+R+1, least, tmp[j+1].first, 1);
//加上新的贡献
}
for(int j=0; j<Q[i].size(); j++)
for(int k=1; k<=10; k++)
ans[Q[i][j].second][k]=bit[k-1].query(Q[i][j].first)%10+'0';
las[x]=i;
}
for(int i=1; i<=m; i++) printf("%s\n", ans[i]+1);
return 0;
}
标签:tmp,10,ch,即便,int,Ynoi2015,++,看不到,include
From: https://www.cnblogs.com/sizeof127/p/16609235.html