hdu5715XOR 游戏

链接：http://acm.hdu.edu.cn/showproblem.php?pid=5715

题意：中文题。

分析：二分+dp。要求最小值最大，很容易想到二分答案，但是对于答案x怎么判断是否可行呢？我们可以用dp[i][j]表示前i个元素分成j组是否可行来判断。显然dp[i][j]=max(a[i]^a[k])并且要求i-k<=l，k<i和dp[k][j-1]是可行状态。我们只要维护一个a[i-l]~a[i-1]的trie树然后每次用a[i]去找最大的a[i]^a[k]即可。详见代码。O(30*t*n*m)

代码：

#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<bitset>
#include<math.h>
#include<vector>
#include<string>
#include<stdio.h>
#include<cstring>
#include<iostream>
#include<algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;
const int N=10010;
const int mod=100000000;
const int MOD1=1000000007;
const int MOD2=1000000009;
const double EPS=0.00000001;
typedef long long ll;
const ll MOD=1000000007;
const int MAX=2000000010;
const ll INF=1ll<<55;
const double pi=acos(-1.0);
typedef double db;
typedef unsigned long long ull;
int n,m,l,tot,a[N],e[30],f[30*N][2],sum[30*N];
bool q[2][N];
void add(int x) {
    int i,d,w=0;sum[0]++;
    for (i=29;i>=0;i--) {
        d=(x&e[i])>>i;
        if (f[w][d]) { w=f[w][d];sum[w]++; }
        else { f[w][d]=++tot;w=tot;sum[w]=1; }
    }
}
void del(int x) {
    int i,d,w=0;sum[0]--;
    for (i=29;i>=0;i--) {
        d=(x&e[i])>>i;
        w=f[w][d];sum[w]--;
    }
}
int get(int x) {
    int i,d,w=0,ret=0;
    for (i=29;i>=0;i--) {
        d=(x&e[i])>>i;
        if (f[w][d^1]&&sum[f[w][d^1]]) { ret+=e[i];w=f[w][d^1]; }
        else w=f[w][d];
    }
    return ret;
}
int pd(int x) {
    memset(q,0,sizeof(q));
    int i,j,g,now,pre;
    now=pre=1;q[now][0]=1;
    for (i=1;i<=m;i++) {
        pre=now;now^=1;g=0;
        memset(q[now],0,sizeof(q[now]));
        for (j=0;j<=n;j++) {
            if (j>l&&q[pre][j-l-1]) del(a[j-l-1]);
            if (get(a[j])>=x) { g++;q[now][j]=1; }
            if (q[pre][j]) add(a[j]);
        }
        for (j=n-l;j<=n;j++)
        if (q[pre][j]) del(a[j]);
        if (!g) return 0;
    }
    return q[now][n];
}
int main()
{
    int i,t,ca,L,R,mid;
    scanf("%d", &t);tot=0;
    for (i=0;i<30;i++) e[i]=1<<i;
    for (ca=1;ca<=t;ca++) {
        scanf("%d%d%d", &n, &m, &l);
        for (i=1;i<=n;i++) scanf("%d", &a[i]),a[i]^=a[i-1];
        for (i=0;i<=tot;i++) f[i][0]=f[i][1]=sum[i]=0;
        tot=0;L=0;R=1<<30;mid=(L+R)>>1;
        while (L+1<R)
        if (pd(mid)) { L=mid;mid=(L+R)>>1; }
        else { R=mid;mid=(L+R)>>1; }
        printf("Case #%d:\n%d\n", ca, L);
    }
    return 0;
}