传送门
解析:
显然我们需要处理出所有数能够管辖的区间的范围。也就是处理出左右端点。
整除?这个还是可以用单调栈维护的,先讲做法,再证明。
维护栈中每个数在数列中的位置,遇到下个数先判断它是否是栈顶数的倍数,否则一直弹栈顶并判断,直到栈为空。这一轮中弹出去的数能够管辖的边界就在这个数的前一个位置。
考虑最开始栈时空的,所以这个策略得到的栈任意时刻上面的数都是下面的数的倍数。
而一个数倍数的倍数肯定还是它的倍数,所以不可能出现越出实际管辖范围的情况。
而只要是倍数就不会造成弹栈,所以不可能出现在管辖边界之前就弹出去的情况。
综上,单调栈的正确性是毋庸置疑的。
复杂度也因为单调栈而优化到了优秀的。
代码:
#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;
}