[BZOJ4069][Apio2015]巴厘岛的雕塑（二分+DP）

直接DP不可做（因为按位取或不满足最优化原理）。
而求最大/小与/或/异或和，可以考虑利用二分的思想，在二进制意义下从高到低确定答案的每一位。在此题中要求最小或和，因此从高往低考虑到每一位时，这一位能为$0$就为$0$，否则就为$1$。
而如何判断能否为$0$，就可以DP了。设当前考虑到第$x$位（最低位为第$0$位），当前结果为$ans$（$ans$的第$0$位到第$x-1$位都暂时定为$0$）：
下面定义一个$or(a,b,k)$，表示$a$的第$k$位到最高位和$b$的第$k$位到最高位按位或（得出结果的第$0$位到第$k-1$位都为$0$）：

$$or(a,b,k)=(\lfloor\frac{a}{2^k}\rfloor\mbox{or}\lfloor\frac{b}{2^k}\rfloor)\times2^k$$
$f[i][j]$表示到第$i$个雕塑，分成$j$段，并且对于任意的$h>x$，如果$ans$的第$h$位为$0$，那么这$j$段的或和的第$h$位也必须为$0$，在这样的限制下，第$x$位能否为$0$。边界为$f[0][0]=true$。
转移即枚举上一段的开头$k$，如果存在一个$k$同时满足这两个条件：
1、$f[k][j-1]=true$，即如果到第$k$个雕塑，分成$j-1$段时第$x$位能为$0$，到第$i$个雕塑分成$j$段时第$x$位才有可能为$0$。
2、设$sum[]$为雕塑年龄的前缀和，条件为$sum[i]-sum[k]$的第$x$位为$0$。
3、对于任意的$h>x$，如果$ans$的第$h$位为$0$，那么$sum[i]-sum[k]$的第$h$位也必须为$0$，也就是：
$or(sum[i]-sum[k],ans,x)=ans$
满足这$3$个条件则$f[i][j]=true$。
如果最后存在一个$A\leq i\leq B$使得$f[n][i]=true$，则第$x$为定为$0$。
时间复杂度为$O(n^3\log\sum_{i=1}^nY_i)$，无法通过$2000$的数据点。
但是$2000$的数据点的特性是$A=1$，因此设$g[i]$表示到了第$i$位，满足条件并且第$x$位为$0$，最少需要分成几段，$g[0]=0$。
这时候转移条件和$f$基本相似，即枚举上一段的开头$j$，如果这个$j$满足条件，则$g[i]=\min(g[i],g[j]+1)$。
最后如果$g[n]\leq B$，那么第$x$位定为$1$。
因此$A=1$时复杂度$O(n^2\log\sum_{i=1}^nY_i)$。
代码：

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
inline int read() {
    int res = 0; bool bo = 0; char c;
    while (((c = getchar()) < '0' || c > '9') && c != '-');
    if (c == '-') bo = 1; else res = c - 48;
    while ((c = getchar()) >= '0' && c <= '9')
        res = (res << 3) + (res << 1) + (c - 48);
    return bo ? ~res + 1 : res;
}
typedef long long ll;
const int N = 2005, INF = 0x3f3f3f3f;
int n, A, B, a[N], g[N]; bool f[N][N];
ll sum[N], ans;
ll Or(ll x, ll y, int k) {
    x >>= k; y >>= k; return (x | y) << k;
}
bool _(int x) {
    int i, j, k; for (i = 0; i <= n; i++) for (j = 0; j <= n; j++)
        f[i][j] = 0; f[0][0] = 1;
    for (j = 1; j <= B; j++) for (i = 1; i <= n; i++)
    for (k = j - 1; k < i; k++) if (f[k][j - 1] &&
        Or(sum[i] - sum[k], ans, x) == ans)
            f[i][j] = 1;
    for (i = A; i <= B; i++) if (f[n][i]) return 1;
    return 0;
}
bool __(int x) {
    int i, j; memset(g, INF, sizeof(g)); g[0] = 0;
    for (i = 1; i <= n; i++) for (j = 0; j < i; j++)
        if (Or(sum[i] - sum[j], ans, x) == ans)
            g[i] = min(g[i], g[j] + 1);
    return g[n] <= B;
}
bool Try(int x) {
    if (A == 1) return __(x);
    return _(x);
}
int main() {
    int i; n = read(); A = read(); B = read();
    for (i = 1; i <= n; i++) sum[i] = sum[i - 1] + (a[i] = read());
    for (i = 43; i >= 0; i--) if (!Try(i)) ans |= 1ll << i;
    cout << ans << endl;
    return 0;
}