BZOJ_P4416 [SHOI2013]阶乘字符串(状态压缩动态规划)

BZOJ传送门
Time Limit: 10 Sec Memory Limit: 256 MB
Submit: 196 Solved: 74
[Submit][Status][Discuss]
Description
给定一个由前n个小写字母组成的串S。
串S是阶乘字符串当且仅当前n个小写字母的全排列（共n!种）都作为S的子序列（可以不连续）出现。
由这个定义出发，可以得到一个简单的枚举法去验证，但是它实在太慢了。所以现在请你设计一个算法，在1秒内判断出给定的串是否是阶乘字符串。

Input
输入第1行一个整数T，表示这个文件中会有T组数据。
接下来分T个块，每块2行：
第1行一个正整数n，表示S由前n个小写字母组成。
第2行一个字符串S。

Output
对于每组数据，分别输出一行。每行是YES或者NO，表示该数据对应的串S是否是阶乘字符串。

Sample Input
2
2
bbaa
2
aba

Sample Output
NO
YES

【样例解释】
第一组数据中，ab这个串没有作为子序列出现。
N<=26
T<=5
|S|<=450

HINT
Source

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
#define N 455
char s[N];int p[N][26],f[1<<21];
int w[26];
int t,n,m,lim;char ch;
int main(){
    scanf("%d",&t);for(int i=0;i<26;i++) w[i]=1<<i;
    while(t--){
        scanf("%d",&n),m=0,lim=(1<<n);
        memset(f,0,sizeof(f));memset(p,0,sizeof(p));
        while(ch>'z'||ch<'a') ch=getchar();
        while(ch>='a'&&ch<='z') m++,s[m]=ch-'a',p[m-1][ch-'a']=m,ch=getchar();
        if(n>21){puts("NO");continue;}
        for(int i=0;i<n;++i) p[m][i]=p[m+1][i]=m+1;
        for(int j=0;j<n;j++) for(int i=m-1;~i;--i)if(!p[i][j]) p[i][j]=p[i+1][j];
        for(int S=1;S<lim;++S){
            for(int i=0;i<n;i++){
                if(w[i]&S)
                    f[S]=max(f[S],p[f[S^w[i]]][i]);
            }
        }
        if(f[lim-1]!=m+1) puts("YES");else puts("NO");
    }
    return 0; 
}