4069: [Apio2015]巴厘岛的雕塑

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

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

题解：

对于n<=200的：

从高到低枚举所有位，要使得对于每组该位最好都为0。

令dp[i][j]表示前i棵树分为j组是否能够满足条件。枚举k，表示令k+1~j为一组。

满足条件当且仅当dp[k][j-1]满足条件且k+1~j为一组时满足当前所有位的条件。

时间复杂度为n^3*位数。

对于A=1的：

仍然从高到低枚举每一位。由于A=1的存在，我们可以令dp[i]表示前i个数在满足条件的情况下最少分为几组。

因此复杂度降为n^2*位数。

code：

 #include<cstdio>

 #include<iostream>

 #include<cmath>

 #include<cstring>

 #include<algorithm>

 using namespace std;

 typedef long long int64;

 char ch;

 bool ok;

 void read(int &x){

     for (ok=,ch=getchar();!isdigit(ch);ch=getchar()) if (ch=='-') ok=;

     for (x=;isdigit(ch);x=x*+ch-'',ch=getchar());

     if (ok) x=-x;

 }

 const int maxn=;

 int n,a,b,v[maxn],g[maxn];

 int64 sum[maxn];

 bool f[][];

 bool check(bool flag,int64 ans){

     if (flag){

         memset(g,,sizeof(g));

         g[]=;

         for (int i=;i<=n;i++) for (int j=;j<i;j++) if (((sum[i]-sum[j])|ans)==ans) g[i]=min(g[i],g[j]+);

         return g[n]<=b;

     }

     else{

         memset(f,,sizeof(f));

         f[][]=;

         for (int i=;i<=n;i++) for (int j=;j<=b;j++) for (int k=;k<i;k++) if (((sum[i]-sum[k])|ans)==ans) f[i][j]|=f[k][j-];

         for (int i=a;i<=b;i++) if (f[n][i]) return true;

         return false;

     }

 }

 int main(){

     read(n),read(a),read(b);

     for (int i=;i<=n;i++) read(v[i]),sum[i]=sum[i-]+v[i];

     int64 ans=(~0ULL>>);

     for (int i=;i>=;i--){

         ans^=1LL<<i;

         if (!check(a==,ans)) ans^=1LL<<i;

     }

     printf("%lld\n",ans);

     return ;

 }