高明的最长公共子序列

国家集训队论文。？

https://wenku.baidu.com/view/ed99e4f77c1cfad6195fa776.html

XJOI题 重复串。 虽然不满的1000^3复杂度能水过，但还是学了一下。。

以后压位别压太死了(用unsigned压32位这种)，可能会引出很多麻烦，比如自然溢出之类的。。

放代码。

#include <bits/stdc++.h>
#define U unsigned
#define ll long long
const U h=2147483648;
using namespace std;
U g[26][34],u[34],X[34],H[32];
int n,o,mx,f[65536]; char S[1005];
int main(){
    scanf("%d%s",&n,S+1);
    H[0]=1; for (int i=1;i<32;++i) H[i]=H[i-1]<<1;
    for (int i=1;i<65536;++i) f[i]=f[i>>1]+(i&1);
    for (int I=1;I<=n;++I){
        memset(g,0,sizeof g); memset(u,0,sizeof u);
        for (int i=I,j=0;i<=n;++i,++j) g[S[i]-'a'][j>>5]|=H[j&31];
        for (int k=1;k<I;++k){
            o=S[k]-'a';
            for (int i=0;i<34;++i) X[i]=u[i]|g[o][i];
            for (int i=33;~i;--i) u[i+1]|=(bool)(u[i]&h),u[i]<<=1;
            u[0]|=1; o=n-I+2; X[o>>5]|=H[o&31];
            (u[o>>5]|=H[o&31])^=H[o&31];
            for (int i=0,j=0;i<34;++i)
                if ((ll)X[i]<(ll)u[i]+j) u[i]=X[i]-u[i]-j,j=1;
                else u[i]=X[i]-u[i]-j,j=0;
            (u[o>>5]|=H[o&31])^=H[o&31];
            (X[o>>5]|=H[o&31])^=H[o&31];
            for (int i=0;i<34;++i) (u[i]^=X[i])&=X[i];
        }
        o=0;
        for (int i=0;i<34;++i) o+=f[u[i]&65535]+f[u[i]>>16];
        mx=max(mx,o);
    }
    printf("%d
",n-mx*2);
    return 0;
}