HDU 5845 Best Division

$dp$，字典树。

$dp$递推式很容易知道。dp[i]=max{dp[j]+1} a[j]^..^a[i]<=X,并且$[j,i]$长度不能超过$L$。

但是暴力来复杂度极高，所以需要用字典树维护这个东西。将前缀异或和插入到字典树中，然后不断维护$a[i]$位置之前$L$个前缀异或和就好了。

跑了$405ms$，第一次排到第一，有点激动~~~

#pragma comment(linker, "/STACK:1024000000,1024000000")
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#include<map>
#include<set>
#include<queue>
#include<stack>
#include<iostream>
using namespace std;
typedef long long LL;
const double pi=acos(-1.0),eps=1e-;
void File()
{
    freopen("D:\\in.txt","r",stdin);
    freopen("D:\\out.txt","w",stdout);
}
template <class T>
inline void read(T &x)
{
    char c = getchar(); x = ;while(!isdigit(c)) c = getchar();
    while(isdigit(c)) { x = x *  + c - ''; c = getchar();  }
}

const LL mod=(LL);
const int maxn=;
LL a[maxn],X,P,Q;
int T,n,L,dp[maxn],h;
struct Node { int cnt,mx,nx[]; }s[*maxn];
int sz; bool fail,f[maxn];

int add()
{
    s[sz].mx=; s[sz].cnt=;
    s[sz].nx[]=s[sz].nx[]=-; sz++;
    return sz-;
}

void Delete(LL num)
{
    int p=;
    for(int i=;i>=;i--)
    {
        int x; LL g=(LL)(<<i); if(num&g) x=; else x=;
        int t=p; p=s[p].nx[x]; s[p].cnt--;
        if(s[p].cnt==) { s[t].nx[x]=-; return; }
    }
}

void Insert(LL num,int val)
{
    int p=;
    for(int i=;i>=;i--)
    {
        int x; LL g=(LL)(<<i); if(num&g) x=; else x=;
        if(s[p].nx[x]==-) s[p].nx[x]=add();
        p=s[p].nx[x]; s[p].cnt++; s[p].mx=max(s[p].mx,val);
    }
}

void Find(LL num)
{
    int p=;
    for(int i=;i>=;i--)
    {
        int x; LL g=(LL)(<<i); if(num&g) x=; else x=;
        int y; if(X&g) y=; else y=;

        if(y==)
        {
            p=s[p].nx[x]; if(p==-) return;
            if(i==) h=max(h,s[p].mx),fail=;
        }
        else
        {
            if(s[p].nx[x]!=-) h=max(h,s[s[p].nx[x]].mx),fail=;
            p=s[p].nx[x^]; if(p==-) return;
            if(i==) h=max(h,s[p].mx),fail=;
        }
    }
}

int main()
{
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%lld%d",&n,&X,&L);
        scanf("%lld%lld%lld",&a[],&P,&Q);
        for(int i=;i<=n;i++) a[i]=(a[i-]*P%mod+Q)%mod;

        sz=; s[sz].cnt=; s[sz].nx[]=s[sz].nx[]=-; sz++;

        for(int i=;i<=n;i++) f[i]=;
        for(int i=;i<=n;i++)
        {
            if(i-L->=&&f[i-L-] ) Delete(a[i-L-]);
            dp[i]=; a[i]=a[i]^a[i-]; if(a[i]<=X) dp[i]=;
            h=; fail=; Find(a[i]); if(fail==) dp[i]=h+;
            if(dp[i]!=) { Insert(a[i],dp[i]); f[i]=; }
        }

        printf("%d\n",dp[n]);
    }
    return ;
}