• [NOI2009]二叉查找树


    题目大意:

    给定一棵严格的treap,父亲节点的优先级必然小于儿子节点的。权值按照二叉树的定义,左儿子小于父亲小于右儿子。

    深度从1开始定义,每个点除优先级、数值之外,还有一个访问频度。

    访问频度所产生的代价是:访问频度*该点深度(这和事实相符)

    可以用给定的k的代价,修改任意个点的优先级为任意实数(当然,修改优先级,树的形态,各点深度就可能变化了)

    最终的总代价为:频度产生代价+修改代价。

    最小化这个总代价。

    N<=70,1<=K<=30000000

    分析:

    平衡树是一个动态的数据结构,难以抓住形态的变化,也不方便记录深度之类。所以必须抓住不变的量当做突破口。

    不管平衡树怎么转,根据二叉树的定义,它的中序遍历一定是不变的。

    所以我们可以找到这棵树的中序遍历,就把这棵树变成了一个静态的区间,只不过每个区间所代表的点的优先级可能会变。

    发现,每一个连续的子区间,都对应treap的连续一部分。可以把小的区间先建树,再把大的区间用小的区间合并。我们合并的时候枚举的划分点,就是这部分treap的树根

    区间DP顺理成章。

    除了f[l][r]之外,为了维护优先级的关系,必然要再记录一维。

    发现,只要根节点的优先级确定,子树的优先级的范围就确定了。

    所以考虑记录根节点优先级。(这里优先级只考虑相对大小,而且范围又大,所以要离散化为1~n)

    但是,朴素的f[l][r][w]中,w单单记录根节点优先级的话,由于子树所有大于w的都可以转移,还要for一遍。状态n^3,转移n^2,会爆。

    所以,我们令f[l][r][w]表示,将l~r这段区间建成treap,其中根节点优先级大于等于w的最小代价。

    根据枚举的根节点是否修改,可以设计转移方程是:

    修改:

    f[l][r][w]=min(f[l][r][w],f[l][k-1][w]+f[k+1][r][w]+K+sum[r]-sum[l-1])  ————其中,sum[i]表示,区间中,1~i的访问频度和

    当划分点的优先级ch[k]大于w时,可以不修改。

    f[l][r][w]=min(f[l][r][w],f[l][k-1][ch[k]]+f[k+1][r][ch[k]]+sum[r]-suim[l-1])

    最后答案就是:f[1][n][1];

    注意,o循环的时候,必须倒序!!因为ch[k]>=o时候,要从o更大的地方获取最小值,必须先把o较大的处理完。

    代码1(未简化):

    #include<bits/stdc++.h>
    using namespace std;
    typedef long long ll;
    const int N=77;
    const ll inf=2e18;
    ll f[N][N][N];
    int n;
    ll m;
    int a[N];
    int tot;
    int w[N],p[N],d[N];
    int prio[N];
    ll sum[N];
    int ch[N];//离散化后的优先级 
    bool cmp(int a,int b)
    {
        return w[a]<w[b];
    }
    int main()
    {
        scanf("%d%lld",&n,&m);
        for(int i=1;i<=n;i++) scanf("%d",&w[i]),a[i]=i;
        for(int i=1;i<=n;i++) scanf("%d",&p[i]),prio[i]=p[i];
        for(int i=1;i<=n;i++) scanf("%d",&d[i]);
        sort(a+1,a+n+1,cmp);
        sort(prio+1,prio+n+1);
        for(int i=1;i<=n;i++) ch[i]=lower_bound(prio+1,prio+n+1,p[a[i]])-prio;//离散化 
        
        for(int i=1;i<=n;i++)
         for(int j=1;j<=n;j++)
          for(int k=1;k<=n;k++)
            f[i][j][k]=inf;   
        for(int i=1;i<=n;i++)
         for(int k=n;k>=1;k--)
         {
             if(k<=ch[i]) f[i][i][k]=d[a[i]];
             else f[i][i][k]=d[a[i]]+m;//注意,k>ch的时候,不一定是+oo,可以通过修改改变 
         }//l=1的初值 
        
        for(int i=1;i<=n;i++)
         sum[i]=sum[i-1]+d[a[i]];//前缀和 
        for(int l=2;l<=n;l++)
         for(int i=1;i<=n;i++)
        {
            int j=l+i-1;
            if(j>n) break;
            if(l==2)//长度为二的时候,只能二并一 
            {
               for(int o=n;o>=1;o--)
                {
                    f[i][j][o]=min(f[i][j][o],f[i][i][o]+m+sum[j]-sum[i-1]);
                    f[i][j][o]=min(f[i][j][o],f[j][j][o]+m+sum[j]-sum[i-1]);
                    if(ch[i]>=o) f[i][j][o]=min(f[i][j][o],f[i][i][ch[i]]+f[j][j][ch[i]]+sum[j]-sum[i-1]-d[a[i]]);
                    if(ch[j]>=o) f[i][j][o]=min(f[i][j][o],f[i][i][ch[j]]+f[j][j][ch[j]]+sum[j]-sum[i-1]-d[a[j]]);
                }
            }
            else{
               for(int o=n;o>=1;o--)
                for(int k=i;k<=j;k++)
               {
                    if(k==i)//k在端点处,只能用端点和右边所有部分合并 
                    {
                        f[i][j][o]=min(f[i][j][o],f[i+1][j][o]+m+sum[j]-sum[i-1]);
                        if(ch[i]>=o) f[i][j][o]=min(f[i][j][o],f[i][i][ch[i]]+f[i+1][j][ch[i]]+sum[j]-sum[i-1]-d[a[i]]);
                 }
                 else if(k==j)//同理 
                 {
                     f[i][j][o]=min(f[i][j][o],f[i][j-1][o]+m+sum[j]-sum[i-1]);
                     if(ch[j]>=o) f[i][j][o]=min(f[i][j][o],f[i][j-1][ch[j]]+f[j][j][ch[j]]+sum[j]-sum[i-1]-d[a[j]]);
                 }
                 else{//正宗转移方程 
                     f[i][j][o]=min(f[i][j][o],f[i][k-1][o]+f[k+1][j][o]+m+sum[j]-sum[i-1]);
                     if(ch[k]>=o) f[i][j][o]=min(f[i][j][o],f[i][k-1][ch[k]]+f[k+1][j][ch[k]]+sum[j]-sum[i-1]);
                 }
               }
            }
        }
        printf("%lld",f[1][n][1]);
        return 0;
    }

    太恶心了。为了保证l<=r,做出了巨大的讨论。

    其实不用这么麻烦,只要让l>r的时候,赋值为0就好,相当于不存在。根本不影响答案。

    代码2(化简)

    #include<bits/stdc++.h>
    using namespace std;
    typedef long long ll;
    const int N=77;
    const ll inf=2e18;
    ll f[N][N][N];
    int n;
    ll m;
    int a[N];
    int tot;
    int w[N],p[N],d[N];
    int prio[N];
    ll sum[N];
    int ch[N];//离散化后的优先级 
    bool cmp(int a,int b)
    {
        return w[a]<w[b];
    }
    int main()
    {
        scanf("%d%lld",&n,&m);
        for(int i=1;i<=n;i++) scanf("%d",&w[i]),a[i]=i;
        for(int i=1;i<=n;i++) scanf("%d",&p[i]),prio[i]=p[i];
        for(int i=1;i<=n;i++) scanf("%d",&d[i]);
        sort(a+1,a+n+1,cmp);
        sort(prio+1,prio+n+1);
        for(int i=1;i<=n;i++) ch[i]=lower_bound(prio+1,prio+n+1,p[a[i]])-prio;
        
        for(int i=1;i<=n;i++)
         for(int j=i;j<=n;j++)
          for(int o=0;o<=n;o++)
             f[i][j][o]=inf;
           for(int i=1;i<=n;i++)
            for(int o=0;o<=n;o++)
             f[i][i-1][o]=0;//其实这步不需要,因为上面就没有给它赋值,只是在这里强调一下。 
             
        for(int i=1;i<=n;i++)
         sum[i]=sum[i-1]+d[a[i]];
          for(int o=n;o>=1;o--)
         for(int i=n;i>=1;i--)
          for(int j=i;j<=n;j++)
            for(int k=i;k<=j;k++) 
            {
              f[i][j][o]=min(f[i][j][o],f[i][k-1][o]+f[k+1][j][o]+m+sum[j]-sum[i-1]);
              if(ch[k]>=o) f[i][j][o]=min(f[i][j][o],f[i][k-1][ch[k]]+f[k+1][j][ch[k]]+sum[j]-sum[i-1]);
            }//不放心,可以考虑代入长度小于等于2的情况。0的作用就出来了。 
            //连初始化l=1都省了。 
        printf("%lld",f[1][n][1]);
        return 0;
    }

    总结:

    1.对于琢磨不透的变化,一定有不变的东西。一定要抓住其中的不变量,作为突破口。
    2.循环顺序要注意,一个是不能有后效性,一个是要保证能影响到这个状态的所有状态都处理完了。

    3.注意考虑清楚所有可能转移的方式。

  • 相关阅读:
    Jhipster 一个Spring Boot + Angular/React 全栈框架
    从企业架构到智慧油田的理论与实践
    DevExpress VCL 的 cxDBTreeList 的使用方法
    Delphi 三层框架 DataSnap 的服务器端设置
    雷达方程
    雷达方向信号产生
    目标
    流水线&并行处理
    CORDIC原理与FPGA实现(2)
    CORDIC原理与FPGA实现(1)
  • 原文地址:https://www.cnblogs.com/Miracevin/p/9062856.html
Copyright © 2020-2023  润新知