st表的概念
st表可以用于处理静态(不修改)的rmq问题(区间最值问题)O(logn)预处理 ,O(1)查询。是一种操作简单的数据结构。
st表的实现
st表通过倍增,区间dp来实现。
st表预处理:
void init()
{
for(int i=1;i<=lg[n];i++)
{
for(int j=1;j+(1<<i)-1<=n;j++)
{
st[j][i]=max(st[j][i-1],st[j+(1<<(i-1))][i-1]);
st[j][i]=min(st[j][i-1],st[j+(1<<(i-1))][i-1]);
st[j][i]=gcd(g[j][i-1],g[j+(1<<(i-1))][i-1]);
}
}
}
st表查询:
int l,r;
scanf("%d %d",&l,&r);
int len=lg[r-l+1];
int x=max(h[l][len],h[r-(1<<len)+1][len]);
int y=min(s[l][len],s[r-(1<<len)+1][len]);
int gcd=gcd(st[l][len],st[r-(1<<len)+1][len]);
例题:codeforces 359D - Pair of Numbers
题意:
找到最长的区间【L,R】,并且满足区间中存在一个能将所有数整除的数,输出所有满足的答案的起始点。
题解:
我们发现区间长度存在单调性,故对区间长度进行二分。
在check函数中对该长度的每一个起点进行枚举。发现区间的最大公约数等于区间的最小值时,该区间成立。我们就只需要对该区间的最小值和最大公约数用st表维护。O(1)查询即可。
代码:
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const LL N =3e5+5;
int a[N],g[N][20],st[N][20],lg[N]={-1},n,ans;
vector<int>res;
int gcd(int a,int b)
{
return b?gcd(b,a%b):a;
}
void init()
{
for(int i=1;i<=lg[n];i++)
{
for(int j=1;j+(1<<i)-1<=n;j++)
{
st[j][i]=min(st[j][i-1],st[j+(1<<(i-1))][i-1]);
g[j][i]=gcd(g[j][i-1],g[j+(1<<(i-1))][i-1]);
}
}
}
bool check(int x)
{
vector<int>vec;
int len=lg[x+1];
for(int i=1;i+x<=n;i++)
{
int minn=min(st[i][len],st[i+x-(1<<len)+1][len]);
int tmp=gcd(g[i][len],g[i+x-(1<<len)+1][len]);
//cout<<i<<" "<<len<<" "<<minn<<" "<<tmp<<endl;
if(tmp==minn)vec.push_back(i);
}
if(vec.size())
{
if(x>=ans) res=vec;
return 1;
}
return 0;
}
void solve()
{
scanf("%lld",&n);
for(int i=1;i<=n;i++)
{
scanf("%lld",&a[i]);
lg[i]=lg[i/2]+1;
g[i][0]=st[i][0]=a[i];
res.push_back(i);
}
init();
int l=0,r=n-1;
while(l<r)
{
int mid=(l+r+1)>>1;
if(check(mid))
{
ans=max(ans,mid);
l=mid;
}
else r=mid-1;
}
printf("%d %d\n",res.size(),ans);
for(auto i:res) printf("%d ",i);
}
int main()
{
int t=1;
//scanf("%d",&t);
while(t--)
{
solve();
}
return 0;
}