题目大意:
题目链接:https://jzoj.net/senior/#main/show/6294
给出序列
a
a
a,求最长的子序列
a
l
∼
a
r
a_l\sim a_r
al∼ar,使得
l
∼
r
l\sim r
l∼r中有一个位置
k
k
k能被所有
l
∼
r
l\sim r
l∼r中的数整除。
思路:
题外话:和 JZOJ3895 重题了啊。那道题还是之前B组的。
显然题目就是要我们求出最长的区间
[
l
,
r
]
[l,r]
[l,r]使得
g
c
d
(
a
l
∼
a
r
)
=
m
i
n
(
a
l
∼
a
r
)
gcd(a_l\sim a_r)=min(a_l\sim a_r)
gcd(al∼ar)=min(al∼ar)。
直接ST表求出区间
g
c
d
gcd
gcd和最小值,然后二分答案,转换成一个判定性问题。这个显然是满足单调性的。
时间复杂度
O
(
n
log
2
n
)
O(n\log^2 n)
O(nlog2n)
**出题人a范围写错,用long
long就T了,int就过了
\color{white}\texttt{**出题人a范围写错,用long long就T了,int就过了}
**出题人a范围写错,用long long就T了,int就过了
代码:
#pragma GCC optimize("Ofast")
#pragma GCC optimize("inline")
#include <queue>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int N=500010,LG=20;
int Min[N][LG+1],Gcd[N][LG+1];
int n,l,r,mid;
queue<int> q;
bool check(int len)
{
int k=log2(len);
bool ok=0;
for (int i=1;i<=n-len+1;i++)
{
int minn=min(Min[i][k],Min[i+len-(1<<k)][k]);
int GCD=__gcd(Gcd[i][k],Gcd[i+len-(1<<k)][k]);
if (minn==GCD)
{
if (!ok)
{
ok=1;
while (q.size()) q.pop();
}
q.push(i);
}
}
return ok;
}
int main()
{
freopen("point.in","r",stdin);
freopen("point.out","w",stdout);
scanf("%d",&n);
memset(Min,0x3f3f3f3f,sizeof(Min));
for (int i=1;i<=n;i++)
{
scanf("%d",&Min[i][0]);
Gcd[i][0]=Min[i][0];
}
for (int j=1;j<=LG;j++)
for (int i=1;i+(1<<(j-1))<=n;i++)
{
Min[i][j]=min(Min[i][j-1],Min[i+(1<<(j-1))][j-1]);
Gcd[i][j]=__gcd(Gcd[i][j-1],Gcd[i+(1<<(j-1))][j-1]);
}
l=1; r=n;
while (l<=r)
{
mid=(l+r)>>1;
if (check(mid)) l=mid+1;
else r=mid-1;
}
printf("%d %d\n",q.size(),l-2);
while (q.size())
{
printf("%d ",q.front());
q.pop();
}
return 0;
}