zoukankan      html  css  js  c++  java
  • bzoj4069 [Apio2015]巴厘岛的雕塑

    4069: [Apio2015]巴厘岛的雕塑

    Time Limit: 10 Sec  Memory Limit: 64 MB
    Submit: 468  Solved: 229
    [Submit][Status][Discuss]

    Description

    印尼巴厘岛的公路上有许多的雕塑,我们来关注它的一条主干道。
    在这条主干道上一共有 N 座雕塑,为方便起见,我们把这些雕塑从 1 到 N 连续地进行标号,其中第 i 座雕塑的年龄是 Yi 年。为了使这条路的环境更加优美,政府想把这些雕塑分成若干组,并通过在组与组之间种上一些树,来吸引更多的游客来巴厘岛。
    下面是将雕塑分组的规则:
    这些雕塑必须被分为恰好 X 组,其中 A< = X< = B,每组必须含有至少一个雕塑,每个雕塑也必须属于且只属于一个组。同一组中的所有雕塑必须位于这条路的连续一段上。
    当雕塑被分好组后,对于每个组,我们首先计算出该组所有雕塑的年龄和。
    计算所有年龄和按位取或的结果。我们这个值把称为这一分组的最终优美度。
    请问政府能得到的最小的最终优美度是多少?
    备注:将两个非负数 P 和 Q 按位取或是这样进行计算的:
    首先把 P 和 Q 转换成二进制。
    设 nP 是 P 的二进制位数,nQ 是 Q 的二进制位数,M 为 nP 和 nQ 中的最大值。P 的二进制表示为 pM−1pM−2…p1p0,Q 的二进制表示为 qM−1qM−2…q1q0,其中 pi 和 qi 分别是 P 和 Q 二进制表示下的第 i 位,第 M−1 位是数的最高位,第 0 位是数的最低位。
    P 与 Q 按位取或后的结果是: (pM−1  OR  qM−1)(pM−2 OR qM−2)…(p1 OR q1)(p0 OR q0)。其中:
    0 OR 0=0
    0 OR 1=1
    1 OR 0=1
    1 OR 1=1

    Input

    输入的第一行包含三个用空格分开的整数 N,A,B。

     
    第二行包含 N 个用空格分开的整数 Y1,Y2,…,YN。
     

    Output

    输出一行一个数,表示最小的最终优美度。

     

    Sample Input

    6 1 3
    8 1 2 1 5 4

    Sample Output

    11

    explanation

    将这些雕塑分为 2 组,(8,1,2) 和 (1,5,4),它们的和是 (11) 和 (10),最终优美度是 (11 OR 10)=11。(不难验证,这也是最终优美度的最小值。)

    HINT

     子任务 1 (9 分)


    1< = N< = 20

    1< = A< = B< = N

    0< = Yi< = 1000000000

    子任务 2 (16 分)

    1< = N< = 50

    1< = A< = B< = min{20,N}

    0< = Yi< = 10

    子任务 3 (21 分)

    1< = N< = 100

    A=1

    1< = B< = N

    0< = Yi< = 20

    子任务 4 (25 分)

    1< = N< = 100

    1< = A< = B< = N

    0< = Yi< = 1000000000

    子任务 5 (29 分)

    1< = N< = 2000

    A=1

    1< = B< = N

    0< = Yi< = 1000000000
    分析:非常好的一道题.
       分段求最值,显然是用dp.一开始我想着用常用的状态f[i][j]表示前i个数分成了j段的答案,那么f[i][j] = min{f[k][j - 1] | (sum[i] - sum[k - 1]),f[k][j]}. O(n^3)的做法. 
       这道题显然不能用上面的方法.这道题的特殊之处在于它的操作是or操作.一般遇到and,or,xor操作求最值就考虑把数变成二进制,一位一位地确定,判断可行性. 对于这道题的做法也是一样了,从高位到低位依次确定,每次确定当前位是否能为0,如果不行就为1,继续判断下一位. 那么状态就变成了f[i][j]表示前i个数分成了j段且当前位以及前位是否都符合条件.这里不单单要使得当前位满足条件,之前确定的位也必须满足条件!
       如何判断前位是否满足条件呢? 如果前位被确定是0,那么sum[i] - sum[k-1]在这一位上就必须是0,位运算|判断一下就好了. 
       最后看f[n][i]是否有一个为1就可以了(A ≤ i ≤ B),时间复杂度是O(n^2)的.
       上面这个算法过不了子问题五.子问题五的特殊点在于A = 1.之前的做法用第二维来保证最后分的段数一定在[A,B]之间.状态很多,但是状态表示的答案很少(只有0/1).现在既然A=1了,就说明不需要考虑下界了,那么满足条件的情况下使得划分的最少段数≤B就可以了.
    令g[i]表示前i个数的可行划分方案的最少段数. 如果[j,i]可以划分成一段,那么g[i] = min{g[i],g[j] + 1}.最后看g[n]是否≤B即可.
       理一下思路:由or操作想到在二进制数上考虑.  由求最值可以想到在二进制上从高位到低位逐位确定判断可行性,想到一个二维dp的方法. 子问题五由A=1可以想到让段数尽量小,调换状态和其表示的答案,其实就是可行性与最优性的转化,只不过第一个dp不能直接用最优性做,而必须要用可行性来做.
       小细节:
       1.位运算不要挤在一起写,优先级容易出问题.
       2.要得到2^i,可以用1<<i,但是当i比较大的时候,要用1LL << i. 一个比较好的做法是直接开一个数组保存.
    #include <cstdio>
    #include <cstring>
    #include <iostream>
    #include <algorithm>
    
    using namespace std;
    
    typedef long long ll;
    const ll maxn = 2010,inf = 1e17;
    ll n,A,B,a[maxn],sum[maxn],ans,f[maxn][maxn],g[maxn],len;
    
    void solve1()
    {
        while (len)
        {
            for (ll i = 1; i <= n; i++)
                for (ll j = 1; j <= i; j++)
                    f[i][j] = 0;
            f[0][0] = 1;
            for (ll i = 1; i <= n; i++)
                for (ll j = 1; j <= i; j++)
                    for (ll k = 0; k < i; k++)
                    {
                        if (f[k][j - 1])
                        {
                            ll t = sum[i] - sum[k];
                            ll temp = t & (1LL << (len - 1));
                            if (temp == 0 && (ans | (t >> len)) == ans)
                                f[i][j] = 1;
                        }
                    }
            bool flag = false;
            for (ll i = A; i <= B; i++)
                if (f[n][i])
                {
                    flag = true;
                    break;
                }
            ans *= 2;
            if (!flag)
                ans |= 1;  //这一位就是1了
            len--;
        }
    }
    
    void solve2()
    {
        while (len)
        {
            for (ll i = 1; i <= n; i++)
                g[i] = inf;
            g[0] = 0;
            for (ll i = 1; i <= n; i++)
                for (ll j = 0; j < i; j++)
                {
                    ll t = sum[i] - sum[j];
                    ll temp = t & (1LL << (len - 1));
                    if (temp == 0 && (ans | (t >> len)) == ans)
                        g[i] = min(g[i],g[j] + 1);
                }
            ans *= 2;
            if (g[n] > B)
                ans |= 1;
            len--;
        }
    }
    
    int main()
    {
        scanf("%lld%lld%lld",&n,&A,&B);
        for (ll i = 1; i <= n; i++)
        {
            scanf("%lld",&a[i]);
            sum[i] = sum[i - 1] + a[i];
        }
        ll t = sum[n];
        while (t)
        {
            len++;
            t >>= 1;
        }
        if (A != 1)
            solve1();
        else
            solve2();
        printf("%lld
    ",ans);
    
        return 0;
    }
       
  • 相关阅读:
    php随笔3-thinkphp 学习-ThinkPHP3.1快速入门(1)基础
    自己动手写PHP MVC框架
    Aptana Studio 3 官方汉化包汉化
    PHP使用手册
    关系型数据库:关系模式设计原则
    Git使用教程
    高并发与负载均衡-nginx-session一致性
    高并发与负载均衡-nginx-安装-配置-location-负载均衡
    高并发与负载均衡-nginx-反向代理概念
    Java-笔记2
  • 原文地址:https://www.cnblogs.com/zbtrs/p/8495170.html
Copyright © 2011-2022 走看看