#笛卡尔树,dp#洛谷 7244 章节划分 分析 代码
考虑段数受到答案限制,而答案为最大值的约数,那么枚举答案,
设(dp[i])表示前(i)个位置分完最多可以分多少段只要(dp[n]geq k)即合法。
那么(dp[i]=max{dp[j]}+1,ans|max{a[jsim i]}),这可以用数据结构实现,
当然可以用笛卡尔树来做,维护一个大根堆,如果最大值为答案的约数那就将两边分开,
否则将一边合并到左右两边,然后判断能否分出(k)段
代码
#include <cstdio>
#include <cctype>
#include <cmath>
#define rr register
using namespace std;
const int N=100011; int a[N],n,m,x,bl;
inline signed iut(){
rr int ans=0; rr char c=getchar();
while (!isdigit(c)) c=getchar();
while (isdigit(c)) ans=(ans<<3)+(ans<<1)+(c^48),c=getchar();
return ans;
}
inline signed max(int a,int b){return a>b?a:b;}
struct Descartes{
int st[N],ls[N],rs[N],top,last,root;
inline void build(){
top=0,last=0;
for (rr int i=1;i<=n;++i){
while (a[st[top]]<a[i]) --top;
if (top<last) ls[i]=st[top+1];
if (top) rs[st[top]]=i;
st[last=++top]=i;
}
root=st[1];
}
inline signed query(int now,int l,int r){
if (!now) return 0;
if (a[now]%x==0) return query(ls[now],l,now-1)+query(rs[now],now+1,r)+1;
rr int ans=0;
if (l>1) ans=max(ans,query(rs[now],now+1,r));
if (r<n) ans=max(ans,query(ls[now],l,now-1));
return ans;
}
}Tre;
signed main(){
n=iut(),m=iut();
for (rr int i=1;i<=n;++i) a[0]=max(a[0],a[i]=iut());
Tre.build(),bl=sqrt(a[0]);
for (rr int i=1;i<=bl;++i)
if (a[0]%i==0){
x=a[0]/i;
if (Tre.query(Tre.root,1,n)>=m)
return !printf("%d",x);
}
for (rr int i=bl;i;--i)
if (a[0]%i==0){
x=i;
if (Tre.query(Tre.root,1,n)>=m)
return !printf("%d",x);
}
return 0;
}