題目大意

給定一個序列\(a_1,a_2,...,a_n\)定義它的一個子序列\(a_{i_j},a_{i_2},...,a_{i_k}(i_1<i_2<...<i_k)\)的美麗值為

求子序列中美麗值最大的

「NOIP2021模擬賽8.18 B」序列（seq）

問題解決

我們很容易寫出一個\(DP\)轉移方程是\(dp[i]=max(dp[j]+\left\lfloor\dfrac{a_{i}}{a_j}\right\rfloor)\)

然後發現\(\left\lfloor\dfrac{a_{i}}{a_j}\right\rfloor\)的取值實際只有\(\sqrt{a_i}\)個，也就是要想辦法取出\(\left\lfloor\dfrac{a_{i}}{a_j}\right\rfloor\)值相同的\(a_j\)的最大值

如何求解\(\left\lfloor\dfrac{a_{i}}{a_j}\right\rfloor\)值相同的\(a_j\)使用數論分塊就可以了，複雜度是\(n\sqrt{n}\)的

然後思考查詢和修改，這裡要思考一個問題，由於\(\left\lfloor\dfrac{a_{i}}{a_j}\right\rfloor\)

顯然，可以用樹狀陣列來維護，但是這樣讓時間複雜度到達\(n\sqrt{n}logn\)

每次\(\sqrt{n}\)來修改，自然而然想到分塊，先建立一個\(a[i]\)到\(dp[i]\)的對映，將分成\(\sqrt{n}\)塊，記錄下前\(i\)塊的最大值\(S[i]\)和每個數在塊中的前面的最大值\(s[i]\)所以每次求出\(dp[i]\)

code

#include<bits/stdc++.h>
using namespace std;
const int maxn=100005,maxm=330;
int N,F[maxn],a[maxn],b[maxn];
int tot,c[maxn],M[maxm],m[maxn],max_x;
int st[maxm],ed[maxm],bel[maxn],Ans,sq;
inline int read(){
	int ret=0,f=1;char ch=getchar();
	while(ch>'9'||ch<'0'){if(ch=='-')f=-f;ch=getchar();}
	while(ch<='9'&&ch>='0')ret=ret*10+ch-'0',ch=getchar();
	return ret*f;
}
void updata(int x,int val){
	for(int j=x;j<=ed[bel[x]];j++)m[j]=max(m[j],val);
	for(int j=bel[x];j<=sq;j++)M[j]=max(M[j],val);
	
}
int main(){
	freopen("seq.in","r",stdin);
	freopen("seq.out","w",stdout);
	N=read();
	for(int i=1;i<=N;i++){a[i]=read();max_x=max(max_x,a[i]);}
	sq=sqrt(max_x);
	for(int i=1;i<=sq;i++){
		st[i]=max_x/sq*(i-1)+1;
		ed[i]=max_x/sq*i;
	}ed[sq]=max_x;
	for(int i=1;i<=sq;i++)
		for(int j=st[i];j<=ed[i];j++)
			bel[j]=i;
	memset(M,128,sizeof M);
	memset(m,128,sizeof m);
	updata(a[1],0);
	for(int i=2;i<=N;i++){
		F[i]=F[i-1];
		for(int l=1,r=0;l<=a[i];l++){r=min(a[i]/(a[i]/l),a[i]);int ret=max(M[bel[r]-1],m[r]);F[i]=max(F[i],ret+(a[i]/r));l=r;}
		if(F[i]>m[a[i]]){
			updata(a[i],F[i]);
		}
	}
	printf("%d\n",F[N]);
	return 0;
}