【校内测试】优美的序列
【题目描述】
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;
}
最后运行速度竟然比二分还快,不比单调栈慢但是思维难度降低了不少 。所以说暴力出奇迹。然而这改变不了我考试的时候忘记判重而爆零的事实。 其实这个方法本质上和单调栈是一样的。