AC自动机+状态DP。

虽然很明显的AC自动机+状态DP题,但要分析问题上还是欠缺一点了。一直在犹豫枚举每一个字符选或不选的状态会不会超时,以为会达到状态有2^n,但其实根本没有。因为有很多状态是可以重复的,这是由于每一位字符选或不选,都只能转移动自动机的固定位置,所以,状态是有限的。状态压缩是为了确定必须出现的串是否全部出现。于是可以设状态:

dp[i][j][k]为当前处理到第i个字符,处于自动机的j状态,k为必须出现的字符串是否出现的情况。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm> using namespace std;
const int root=0;
const int inf=1<<30;
int trie[1650][26];
int fail[1650];
int tag[1650],tot;
int que[1650],head,tail;
int kill[1650],must[1650],score[1650];
int need;
struct Status{
int act,sc;
void init(){
act=inf,sc=-inf;
}
}dp[2][1650][260];
char str[110]; void init(int cur){
for(int i=0;i<=tot;i++){
for(int k=0;k<(1<<(need+1));k++)
dp[cur][i][k].init();
}
} void AddTrie(char *s,int sc){
int len=strlen(s);
int i=0,p=root;
while(len--){
if(trie[p][s[i]-'a']==-1){
trie[p][s[i]-'a']=++tot;
kill[tot]=must[tot]=0;
memset(trie[tot],-1,sizeof(trie[tot]));
}
p=trie[p][s[i]-'a'];
i++;
}
if(sc==999){must[p]|=(1<<need);}
else if(sc==-999) kill[p]=1;
else {
score[p]+=sc;
}
} void build_ac(){
head=tail=0;
que[tail++]=root;
while(head!=tail){
int tmp=que[head++];
int p=-1;
for(int i=0;i<26;i++){
if(trie[tmp][i]!=-1){
if(tmp==root) fail[trie[tmp][i]]=root;
else{
p=fail[tmp];
while(p!=-1){
if(trie[p][i]!=-1){
fail[trie[tmp][i]]=trie[p][i];
break;
}
p=fail[p];
}
if(p==-1) fail[trie[tmp][i]]=root;
}
if(kill[fail[trie[tmp][i]]]) kill[trie[tmp][i]]=1;
if(must[fail[trie[tmp][i]]]) must[trie[tmp][i]]|=must[fail[trie[tmp][i]]];
if(score[fail[trie[tmp][i]]]) score[trie[tmp][i]]+=score[fail[trie[tmp][i]]];
que[tail++]=trie[tmp][i];
}
else{
if(tmp==root) trie[tmp][i]=root;
else{
p=fail[tmp];
while(p!=-1){
if(trie[p][i]!=-1){
trie[tmp][i]=trie[p][i];
break;
}
p=fail[p];
}
if(p==-1) trie[tmp][i]=root;
}
}
}
}
} void Cal(){
int cur=0,next=1;
init(cur); init(next);
int len=strlen(str);
dp[cur][0][0].act=0; dp[cur][0][0].sc=0;
for(int i=0;i<len;i++){
// cout<<i<<endl;
for(int j=0;j<=tot;j++){
for(int k=0;k<(1<<(need+1));k++){
if(dp[cur][j][k].act>=inf) continue;
int pn=trie[j][str[i]-'a'];
if(kill[pn]){
if(dp[cur][j][k].act+1<dp[next][j][k].act){
dp[next][j][k].act=dp[cur][j][k].act+1;
dp[next][j][k].sc=dp[cur][j][k].sc;
}
else if(dp[cur][j][k].act+1==dp[next][j][k].act&&dp[cur][j][k].sc>dp[next][j][k].sc){
dp[next][j][k].sc=dp[cur][j][k].sc;
}
}
else{
int tmpk=k|must[pn];
if(dp[cur][j][k].act<dp[next][pn][tmpk].act){
dp[next][pn][tmpk].act=dp[cur][j][k].act;
dp[next][pn][tmpk].sc=dp[cur][j][k].sc+score[pn];
}
else if(dp[cur][j][k].act==dp[next][pn][tmpk].act&&dp[cur][j][k].sc+score[pn]>dp[next][pn][tmpk].sc){
dp[next][pn][tmpk].sc=dp[cur][j][k].sc+score[pn];
}
if(dp[cur][j][k].act+1<dp[next][j][k].act ){
dp[next][j][k].act=dp[cur][j][k].act+1;
dp[next][j][k].sc=dp[cur][j][k].sc;
}
else if(dp[cur][j][k].act+1==dp[next][j][k].act&&dp[cur][j][k].sc>dp[next][j][k].sc){
dp[next][j][k].sc=dp[cur][j][k].sc;
}
}
}
} init(cur);
cur^=1;
next^=1;
}
int minact=inf,maxsc=-inf;
for(int i=0;i<=tot;i++){
if(dp[cur][i][(1<<(need+1))-1].act<minact){
minact=dp[cur][i][(1<<(need+1))-1].act;
maxsc=dp[cur][i][(1<<(need+1))-1].sc;
}
else if(dp[cur][i][(1<<(need+1))-1].act==minact&&maxsc<dp[cur][i][(1<<(need+1))-1].sc)
maxsc=dp[cur][i][(1<<(need+1))-1].sc;
}
if(minact>=inf)
puts("Banned");
else cout<<minact<<" "<<maxsc<<endl;
} int main(){
int T,icase=0;
char s[20];int sc;
scanf("%d",&T);
while(T--){
int n;
tot=0,need=-1;
memset(trie[tot],-1,sizeof(trie[tot]));
memset(score,0,sizeof(score));
memset(kill,0,sizeof(kill));
memset(must,0,sizeof(must));
scanf("%d",&n);
for(int i=0;i<n;i++){
scanf("%s%d",s,&sc);
if(sc==999) need++;
AddTrie(s,sc);
}
// cout<<"YES"<<endl;
build_ac();
scanf("%s",str);
printf("Case %d: ",++icase);
Cal();
}
return 0;
}

  

HDU 4534的更多相关文章

  1. HDOJ 2111. Saving HDU 贪心 结构体排序

    Saving HDU Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total ...

  2. 【HDU 3037】Saving Beans Lucas定理模板

    http://acm.hdu.edu.cn/showproblem.php?pid=3037 Lucas定理模板. 现在才写,noip滚粗前兆QAQ #include<cstdio> #i ...

  3. hdu 4859 海岸线 Bestcoder Round 1

    http://acm.hdu.edu.cn/showproblem.php?pid=4859 题目大意: 在一个矩形周围都是海,这个矩形中有陆地,深海和浅海.浅海是可以填成陆地的. 求最多有多少条方格 ...

  4. HDU 4569 Special equations(取模)

    Special equations Time Limit:1000MS     Memory Limit:32768KB     64bit IO Format:%I64d & %I64u S ...

  5. HDU 4006The kth great number(K大数 +小顶堆)

    The kth great number Time Limit:1000MS     Memory Limit:65768KB     64bit IO Format:%I64d & %I64 ...

  6. HDU 1796How many integers can you find(容斥原理)

    How many integers can you find Time Limit:5000MS     Memory Limit:32768KB     64bit IO Format:%I64d ...

  7. hdu 4481 Time travel(高斯求期望)(转)

    (转)http://blog.csdn.net/u013081425/article/details/39240021 http://acm.hdu.edu.cn/showproblem.php?pi ...

  8. HDU 3791二叉搜索树解题(解题报告)

    1.题目地址: http://acm.hdu.edu.cn/showproblem.php?pid=3791 2.参考解题 http://blog.csdn.net/u013447865/articl ...

  9. hdu 4329

    problem:http://acm.hdu.edu.cn/showproblem.php?pid=4329 题意:模拟  a.     p(r)=   R'/i   rel(r)=(1||0)  R ...

随机推荐

  1. NS2学习笔记(三)

    代理(Agents) 代理是网络层数据包产生和使用的端点,用于实现多层的协议.Agent类部分用OTcl实现,部分用C++实现.C++的代码在~ns/agent.cc和~ns/agent.h中,OTc ...

  2. ASP.NET MVC5 之路由器

    这篇博客介绍的很详细 http://www.cnblogs.com/yaozhenfa/p/asp_net_mvc_route_1.html

  3. 简单的KKL诊断线~~~自己在家都可以制作obd诊断接口了 ~~

    简单的KKL诊断线~~~自己在家都可以制作~~ 适合bmw 07年以前的车型,因为新的车型使用D-can作为诊断接口,所以不能再使用kkl诊断接口不过SB开头的宝马3系还是可以使用的 更多内容欢迎查看 ...

  4. WEB文件上传之apache common upload使用(一)

    文件上传一个经常用到的功能,它有许多中实现的方案. 页面表单 + RFC1897规范 + http协议上传 页面控件(flash/html5/activeX/applet) + RFC1897规范 + ...

  5. html5——全屏显示

    基本概念 1.HTML5规范允许用户自定义网页上任一元素全屏显示 2.requestFullscreen() 开启全屏显示.cancleFullscreen() 关闭全屏显示 3.不同浏览器兼容性不一 ...

  6. 多线程-实现Runnable接口

    当一个任务或者函数多个线程同时调用时仅仅继承Thread是不行的.需要实现Runnable接口. 好处: 1.将线程的任务从线程的子类中分离出来,进行了单独的封装. 按照面向对象的思想将任务封装成对象 ...

  7. 【剑指Offer】66、机器人的运动范围

      题目描述:   地上有一个m行和n列的方格.一个机器人从坐标0,0的格子开始移动,每一次只能向左,右,上,下四个方向移动一格,但是不能进入行坐标和列坐标的数位之和大于k的格子. 例如,当k为18时 ...

  8. 27.6 Parallel的静态For,Foreach和Invoke方法

    static void Main(string[] args) { //for (int i = 0; i < 10000; i++) // DoWork(i); //Parallel.For( ...

  9. / Vijos / 题库 / 1622 / 文件查找(HOI)

    / Vijos / 题库 /1622/文件查找(HOI) 描述 WINDOWS是一个很庞大的操作系统(当然啦,看占的硬盘空间就知道了),比如说,它的文件查找系统.现在,请你用PASCAL或者C或者C+ ...

  10. springboot 使用idea打包 遇到问题

    找了很久错误后来发现添加这三句话就可以了 需要在项目的pom.xml文件中加上第47-49行的3句话