• CF 833B


    互测题T3...

    首先有个dp是非常好想的:

    设dp[i][j]为前j个数分成i组的最大得分,则易得:dp[i][j]=max{dp[i-1][k-1]+num[k][j]},其中,num[k][j]表示从第k个数到第j个数不同值的数量

    而num数组可以预处理出来,时间复杂度O(n^2 k)

    等等,这样好像过不掉这道题啊

    我们发现,max{dp[i-1][k-1]+num[k][j]}这个东西是不是应该用什么维护一下?

    线段树!

    利用线段树,我们可以实现区间求最值!

    我们记录一个位置i上的数a[i]上一次出现的位置为las[i],那么当我们更新dp到位置i时,我们就可以进行转移,而如果一共分了k组,则最后这一组的起点一定在k以后!

    同时,利用las[i],我们可以将las[i]+1~i这一整段区间的值全部+1,因为这一段区间在更新dp时不产生重复

    最后我们在k~i这段区间上做区间查询,求最大值即为dp值

    每次更新完一组的dp值以后,都需要重新建树,类似滚动数组的原理

    还有,在建树时,考虑到转移方程中需要用到的是dp[i-1][k-1],所以在建树时赋值的下标都应当-1以便利用

    代码:

    #include <cstdio>
    #include <cmath>
    #include <cstring>
    #include <cstdlib>
    #include <iostream>
    #include <algorithm>
    #include <queue>
    #include <stack>
    #define rt1 rt<<1
    #define rt2 (rt<<1)|1
    #define ls tree[rt].lson
    #define rs tree[rt].rson
    using namespace std;
    struct Tree
    {
        int lson;
        int rson;
        int maxval;
        int lazy;
    }tree[400005];
    int dp[55][35005];
    int las[35005];
    int p[35005];
    int n,k;
    void buildtree(int l,int r,int rt,int typ)
    {
        ls=l;
        rs=r;
        tree[rt].lazy=0;
        if(l==r)
        {
            tree[rt].maxval=dp[typ][l-1];
            return;
        }
        int mid=(l+r)>>1;
        buildtree(l,mid,rt1,typ);
        buildtree(mid+1,r,rt2,typ);
        tree[rt].maxval=max(tree[rt1].maxval,tree[rt2].maxval);
    }
    void pushdown(int rt)
    {
        int t=tree[rt].lazy;
        tree[rt].lazy=0;
        tree[rt1].lazy+=t;
        tree[rt2].lazy+=t;
        tree[rt1].maxval+=t;
        tree[rt2].maxval+=t;
    }
    void ins(int l,int r,int v,int rt)
    {
        if(ls>r||rs<l)
        {
            return;
        }
        if(ls>=l&&rs<=r)
        {
            tree[rt].lazy+=v;
            tree[rt].maxval+=v;
            return;
        }
        pushdown(rt);
        int mid=(ls+rs)>>1;
        if(l<=mid)
        {
            ins(l,r,v,rt1);
        }
        if(r>mid)
        {
            ins(l,r,v,rt2);
        }
        tree[rt].maxval=max(tree[rt1].maxval,tree[rt2].maxval);
    }
    int query(int l,int r,int rt)
    {
        if(l>rs||r<ls)
        {
            return 0;
        }
        if(l<=ls&&r>=rs)
        {
            return tree[rt].maxval;
        }
        pushdown(rt);
        return max(query(l,r,rt1),query(l,r,rt2));
    }
    int main()
    {
    //    freopen("handsome.in","r",stdin);
    //    freopen("handsome.out","w",stdout);
        scanf("%d%d",&n,&k);
        for(int i=1;i<=n;i++)
        {
            int x;
            scanf("%d",&x);
            las[i]=p[x];
            p[x]=i;
        }
        buildtree(1,n,1,0);
        for(int i=1;i<=k;i++)
        {
            for(int j=1;j<=n;j++)
            {
                int lq=las[j];
                ins(lq+1,j,1,1);
                dp[i][j]=query(i,j,1);
            }
            buildtree(1,n,1,i);
        }
        printf("%d
    ",dp[k][n]);
        return 0;
    }
  • 相关阅读:
    读书是最划算的自我投资,免费送几本最好的Python学习电子书
    一个出身寒门的状元之编程
    幼儿园小朋友都在学人工智能了,已经从小学毕业的我们该如何学?
    python语音识别终极指南
    Python可视化神器——pyecharts的超详细使用指南!
    用 Python 分析过往 36 年春晚节目数据,发现一些趣事
    拉勾网数据分析职位分析
    按时分秒对数据进行分箱
    numpy
    基于协同过滤的推荐系统案列
  • 原文地址:https://www.cnblogs.com/zhangleo/p/9736344.html
Copyright © 2020-2023  润新知