将一个长度为n的序列分为k段
使得总价值最大一段区间的价值表示为区间内不同数字的个数
n<=35000,k<=50
这题的dp是十分显然的，用dp[i][j]表示前i个数字分成j段的最大值
状态转移方程就是 dp[i][j] = max(dp[i - 1][k] + dis[k + 1][j] | (0 <= k <= j - 1))
但是乍一想可能是个O(knn)的算法，仔细一想果然是。
考虑线段树的优化。显然根据状态转移方程，线段树维护的是dp[i - 1][k] + dis[k + 1][j]的最大值
也就是对每一层都build一次线段树，然后以此查询更新，时间复杂度是O(knlnn)，看着十分科学
其中一个难点在于对颜色种类的维护，dis这个35k * 35k的数组肯定是存不下的，仔细一想其实不需要什么花里胡哨的操作，用一个pre数组存储上一个这个颜色出现的地方，当遇到这个颜色的时候，在线段树pre[x] + 1,x的位置全部加上1，意在pre[x] + 1到x出现了这个颜色，种类++；
然后整个思路就出来了，对于每一层状态，用一个build更新上一层的dp[i - 1][k]，dis在查询的过程中动态更新，这种听起来就很像正解的做法一般就是正解。
#include <map>
#include <set>
#include <cmath>
#include <queue>
#include <stack>
#include <vector>
#include <string>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <sstream>
#include <iostream>
#include <algorithm>
#include <functional>
using namespace std;
const int MAXBUF=;char buf[MAXBUF],*ps=buf,*pe=buf+;
inline bool isdigit(const char& n) {return (n>=''&&n<='');}
inline void rnext(){if(++ps==pe)pe=(ps=buf)+fread(buf,sizeof(char),sizeof(buf)/sizeof(char),stdin);}
template <class T> inline bool in(T &ans){
#ifdef VSCode
ans=;T f=;register char c;
do{c=getchar();if ('-'==c)f=-;}while(!isdigit(c)&&c!=EOF);
if(c==EOF)return false;do{ans=(ans<<)+(ans<<)+c-;
c=getchar();}while(isdigit(c)&&c!=EOF);ans*=f;return true;
#endif
#ifndef VSCode
ans =;T f=;if(ps==pe)return false;do{rnext();if('-'==*ps)f=-;}
while(!isdigit(*ps)&&ps!=pe);if(ps==pe)return false;do{ans=(ans<<)+(ans<<)+*ps-;
rnext();}while(isdigit(*ps)&&ps!=pe);ans*=f;return true;
#endif
}const int MAXOUT=;
char bufout[MAXOUT], outtmp[],*pout = bufout, *pend = bufout+MAXOUT;
inline void write(){fwrite(bufout,sizeof(char),pout-bufout,stdout);pout = bufout;}
inline void out_char(char c){*(pout++)=c;if(pout==pend)write();}
inline void out_str(char *s){while(*s){*(pout++)=*(s++);if(pout==pend)write();}}
template <class T>inline void out_int(T x) {if(!x){out_char('');return;}
if(x<)x=-x,out_char('-');int len=;while(x){outtmp[len++]=x%+;x/=;}outtmp[len]=;
for(int i=,j=len-;i<j;i++,j--) swap(outtmp[i],outtmp[j]);out_str(outtmp);}
template<typename T, typename... T2>
inline int in(T& value, T2&... value2) { in(value); return in(value2...); }
#define For(i, x, y) for(int i=x;i<=y;i++)
#define _For(i, x, y) for(int i=x;i>=y;i--)
#define Mem(f, x) memset(f,x,sizeof(f))
#define Sca(x) scanf("%d", &x)
#define Scl(x) scanf("%lld",&x);
#define Pri(x) printf("%d\n", x)
#define Prl(x) printf("%lld\n",x);
#define CLR(u) for(int i=0;i<=N;i++)u[i].clear();
#define LL long long
#define ULL unsigned long long
#define mp make_pair
#define PII pair<int,int>
#define PIL pair<int,long long>
#define PLL pair<long long,long long>
#define pb push_back
#define fi first
#define se second
#define Vec Point
typedef vector<int> VI;
const double eps = 1e-;
const int maxn = ;
const int INF = 0x3f3f3f3f;
const int mod = 1e9 + ;
int N,M,tmp,K;
int a[maxn];
int dp[maxn][]; //以这个结尾
int pre[maxn];
int display[maxn];
struct Tree{
int l,r;
int MAX,lazy;
}tree[maxn * ];
void Pushdown(int t){
if(tree[t].lazy){
tree[t << ].lazy += tree[t].lazy; tree[t << ].MAX += tree[t].lazy;
tree[t << | ].lazy += tree[t].lazy; tree[t << | ].MAX += tree[t].lazy;
tree[t].lazy = ;
}
}
void Pushup(int t){
tree[t].MAX = max(tree[t << ].MAX,tree[t << | ].MAX);
}
void Build(int t,int l,int r,int flag){
tree[t].l = l; tree[t].r = r;
tree[t].lazy = ;
if(l == r){
tree[t].MAX = dp[l - ][flag];
return;
}
int m = (l + r) >> ;
Build(t << ,l,m,flag);
Build(t << | ,m + ,r,flag);
Pushup(t);
}
void update(int t,int l,int r){
if(l <= tree[t].l && tree[t].r <= r){
tree[t].MAX++;
tree[t].lazy++;
return;
}
Pushdown(t);
int m = (tree[t].l + tree[t].r) >> ;
if(r <= m) update(t << ,l,r);
else if(l > m) update(t << | ,l,r);
else{
update(t << ,l,m);
update(t << | ,m + ,r);
}
Pushup(t);
}
int query(int t,int l,int r){
if(l <= tree[t].l && tree[t].r <= r){
return tree[t].MAX;
}
Pushdown(t);
int m = (tree[t].l + tree[t].r) >> ;
if(r <= m){
return query(t << ,l,r);
}else if(l > m){
return query(t << | ,l,r);
}else{
return max(query(t << | ,m + ,r),query(t << ,l,m));
}
}
int main()
{
in(N,K);
For(i,,N) display[i] = ;
For(i,,N){
in(a[i]);
pre[i] = display[a[i]];
display[a[i]] = i;
}
For(i,,K){ //dp[i][j] = max(dp[i - 1][k] + dis[k + 1][j]) (0 <= k <= j - 1));
Build(,,N,i - );
For(j,,N){
update(,pre[j] + ,j);
dp[j][i] = query(,,j);
// cout << dp[j][i] << " ";
}
// cout << endl;
}
Pri(dp[N][K]);
#ifdef VSCode
write();
system("pause");
#endif
return ;
}