2018.11.02【校内模拟】优美的序列（单调栈）

传送门

---

解析：

显然我们需要处理出所有数能够管辖的区间的范围。也就是处理出左右端点。

整除？这个还是可以用单调栈维护的，先讲做法，再证明。

维护栈中每个数在数列中的位置，遇到下个数先判断它是否是栈顶数的倍数，否则一直弹栈顶并判断，直到栈为空。这一轮中弹出去的数能够管辖的边界就在这个数的前一个位置。

考虑最开始栈时空的，所以这个策略得到的栈任意时刻上面的数都是下面的数的倍数。
而一个数倍数的倍数肯定还是它的倍数，所以不可能出现越出实际管辖范围的情况。
而只要是倍数就不会造成弹栈，所以不可能出现在管辖边界之前就弹出去的情况。

综上，单调栈的正确性是毋庸置疑的。

复杂度也因为单调栈而优化到了优秀的$O(n)$。

---

代码：

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define re register
#define gc getchar
#define pc putchar
#define cs const

inline int getint(){
	re int num;
	re char c;
	while(!isdigit(c=gc()));num=c^48;
	while(isdigit(c=gc()))num=(num<<1)+(num<<3)+(c^48);
	return num;
}

cs int N=500005;
int sta[N],a[N];
int lbound[N],rbound[N];
int n,top;
signed main(){
	n=getint();
	for(int re i=1;i<=n;++i){
		a[i]=getint();
	}
	
	for(int re i=1;i<=n;++i){
		while(top&&a[i]%a[sta[top]])rbound[sta[top--]]=i-1;
		sta[++top]=i;
	}
	while(top)rbound[sta[top--]]=n;
	
	for(int re i=n;i;--i){
		while(top&&a[i]%a[sta[top]])lbound[sta[top--]]=i+1;
		sta[++top]=i;
	}
	while(top)lbound[sta[top--]]=1;
	
	vector<int> ans;
	int maxn=0;
	for(int re i=1;i<=n;++i){
		if(rbound[i]-lbound[i]>maxn){
			maxn=rbound[i]-lbound[i];
			ans.clear();
		}
		if(rbound[i]-lbound[i]==maxn){
			ans.push_back(lbound[i]);
		}
	}
	sort(ans.begin(),ans.end());
	ans.erase(unique(ans.begin(),ans.end()),ans.end());
	printf("%d %d\n",ans.size(),maxn);
	for(vector<int>::iterator it=ans.begin();it!=ans.end();++it)printf("%d ",*it);
	
	return 0;
}