P4392 [BOI2007]Sound 静音问题（rmq，st表）

题目描述

数字录音中，声音是用表示空气压力的数字序列描述的，序列中的每个值称为一个采样，每个采样之间间隔一定的时间。

很多声音处理任务都需要将录到的声音分成由静音隔开的几段非静音段。为了避免分成过多或者过少的非静音段，静音通常是这样定义的：m个采样的序列，该序列中采样的最大值和最小值之差不超过一个特定的阈值c。

请你写一个程序，检测n个采样中的静音。

输入格式

第一行有三个整数n，m，c（ 1<= n<=1000000，1<=m<=10000， 0<=c<=10000），分别表示总的采样数、静音的长度和静音中允许的最大噪音程度。

第2行n个整数ai (0 <= ai <= 1,000,000)，表示声音的每个采样值，每两个整数之间用空格隔开。

输出格式

列出了所有静音的起始位置i（i满足max(a[i, . . . , i+m−1]) − min(a[i, . . . , i+m−1]) <= c），每行表示一段静音的起始位置，按照出现的先后顺序输出。如果没有静音则输出NONE。

输入输出样例

输入 #1复制

7 2 0
0 1 1 2 3 2 2

输出 #1复制

2
6

思路：用两个st表记录一下区间的最大最小值，遍历一遍数组，以i为起点长度为m的区间判断最大最小值之差大小是不是<=c即可

坑点：最后一个测试点会mle，所以我们需要尽量用少的空间，

那么我们就不要再设置新的数组来存原数组了，

读入的时候直接存入st[i][0]中

然后我们看m的范围是1e4

那么我们st表的第二维只用开log2（1e4）=13就可以了

/*

 .----------------.  .----------------.  .----------------.  .----------------. 
| .--------------. || .--------------. || .--------------. || .--------------. |
| |  ________    | || |  _________   | || | ____    ____ | || |     ____     | |
| | |_   ___ `.  | || | |_   ___  |  | || ||_   \  /   _|| || |   .'    `.   | |
| |   | |   `. \ | || |   | |_  \_|  | || |  |   \/   |  | || |  /  .--.  \  | |
| |   | |    | | | || |   |  _|  _   | || |  | |\  /| |  | || |  | |    | |  | |
| |  _| |___.' / | || |  _| |___/ |  | || | _| |_\/_| |_ | || |  \  `--'  /  | |
| | |________.'  | || | |_________|  | || ||_____||_____|| || |   `.____.'   | |
| |              | || |              | || |              | || |              | |
| '--------------' || '--------------' || '--------------' || '--------------' |
 '----------------'  '----------------'  '----------------'  '----------------'

*/
#include<iostream>
#include<algorithm>
#include<cstring>
#include<vector>
#include<set>
#include<map>
#include<queue>
#include<deque>
#include<cmath>
#include<stack>
#define lowbit(x) x&(-x)
#define PI 3.1415926535
#define endl "\n"
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
int gcd(int a,int b) {
	return b? gcd(b,a%b):a;
}
/*
int dx[8]={-2,-2,-1,1,2,2,-1,1};
int dy[8]={-1,1,2,2,1,-1,-2,-2};
int dx[4]={0,-1,0,1};
int dy[4]={-1,0,1,0};
int dx[8]={-1,1,0,0,-1,-1,1,1};
int dy[8]={0,0,-1,1,-1,1,-1,1};
*/
//int e[N],ne[N],h[N],idx,w[N];
/*void add(int a,int b,int c){
	e[idx]=b;
	w[idx]=c;
	ne[idx]=h[a];
	h[a]=idx++;
}
*/
const int N=1e6+3;
int m,n,c;
int fi[N][14],fx[N][14];
void init(){
	for(int j=1;j<=20;j++)//最外层枚举的是长度
		for(int i=1;i+(1<<j)-1<=n;i++)
		{
			fi[i][j]=min(fi[i][ j-1],fi[i+(1<<( j - 1 )) ][j-1] );//右边的是i+2^(j-1)->i+(1<<j-1)
				  fx[i][j]=max(fx[i][ j-1],fx[i+(1<<( j - 1 )) ][j-1] );//右边的是i+2^(j-1)->i+(1<<j-1)
		}
}
void sove(){
	cin>>n>>m>>c;
	for(int i=1;i<=n;i++){
		int x;
		cin>>x;
		fi[i][0]=x;
		fx[i][0]=x;
	}
	init();
	int k=log2(m);
	bool f=true;
	for(int i=1;i+m-1<=n;i++){
		int l=i,r=l+m-1;
		int ansi=min(fi[l][k],fi[r-(1<<k)+1][k]);
		int ansx=max(fx[l][k],fx[r-(1<<k)+1][k]);
		if(ansx-ansi<=c){
			cout<<i<<endl;
			f=false;
		}
	}
	if(f)cout<<"NONE"<<endl;
}

int main(){
	ios::sync_with_stdio(false);
	cin.tie() ,cout.tie() ;
	int t=1;
/*	int op=log2(1e4);
	cout<<op<<endl;
	*/
//	cin>>t;
	while(t--){
		sove();
	}
	return 0;
}