【枚举优化】 优美的序列

【校内测试】优美的序列

【题目描述】
L喜欢序列（明显吃饱了没事干 ），他最近在研究一个优美的序列。假设有一个长度是N的序列，对于区间[L,R],当且仅当存在k∈[L,R]，对于任意i∈[L,R]，满足ak | ai，也就是区间里存在一个数字是该区间所有数的约数。
L想知道最长的优美区间是多少。（想知道？不告诉你！ ）
30%数据：N<=30,1<=Ai<=32
60%数据：N<=3000,1<=Ai<=1024
80%数据：N<=300 000,1<=Ai<=104876
100%数据：1<=N<=500 000,1<=Ai<2^31
【输入】
第一行是一个整数n
第二行是n个整数Ai
【输出】
第一行2个数字，m和L，表示最长的优美区间的个数，以及他们的右端点减左端点的值

接下来是一行m个数字，按升序输出这些最长优美区间的左端点
【样例输入】
5
4 6 9 3 6
【样例输出】
1 3
2

最怕考试遇见题目描述很短的题。
这道题正解似乎是一个叫单调栈的东东。这种东西蒟蒻我当然是想不到的（不知道谁吃多了能想到 ），因此我采取了最简单最粗暴最单纯的做法——枚举。逐个枚举每个点能到达的最左和最右位置，用l和r数组记录一下，再来比较就可以啦。显然，直接枚举是过不了的，但是我们可以有一个显而易见的 优化。这里以枚举第i个数最左能影响的位置为例：当a[i]是a[j]的约数时(j<i)，由于l[j]已经更新过，l[i]就可以和l[j]取一个最小值（显然l[j]更小 ），再继续向左拓展，否则停止拓展当前数，开始拓展下一个数。向右拓展也是同样的方法（向右拓展=向左拓展+ctrl c+ctrl v ）。
代码：（忽略第一句注释 ）

//枯藤老树昏鸦，小桥流水人家。 
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<string>
#include<queue>
#include<vector>
#define re register
#define lc (id<<1)
#define rc (id<<1|1)
using namespace std;
int n,m,a[500001],b,c,l[500001],r[500001];
int ans=-1;
int lm[500001];//记录最长长度的左端点
int vis[500001];//判重数组
int main()
{
	//freopen("sequence.in","r",stdin);
	//freopen("sequence.out","w",stdout);
	scanf("%d",&n);
	for(int re i=1;i<=n;i++)
		scanf("%d",&a[i]);
	for(int re i=1;i<=n;i++)
	{
		l[i]=i;
		while(a[l[i]-1]%a[i]==0&&l[i]-1>0)
		{
			l[i]=min(l[i],l[l[i]-1]);
		}
	}
	for(int re i=n;i>=1;i--)
	{
		r[i]=i;
		while(a[r[i]+1]%a[i]==0&&r[i]<n)
		{
			r[i]=max(r[i],r[r[i]+1]);
		}
	}
	for(int re i=1;i<=n;i++)
		ans=max(ans,r[i]-l[i]);
	for(int re i=1;i<=n;i++)
	{
		if(r[i]-l[i]==ans&&!vis[l[i]])
		{
			vis[l[i]]++;
			lm[++m]=l[i];
		}
	}
	sort(lm+1,lm+m+1);
	printf("%d %d\n",m,ans);
	for(int i=1;i<=m;i++)
	{
		printf("%d ",lm[i]);
	}
	return 0;
}

最后运行速度竟然比二分还快，不比单调栈慢但是思维难度降低了不少 。所以说暴力出奇迹。然而这改变不了我考试的时候忘记判重而爆零的事实。 其实这个方法本质上和单调栈是一样的。