ACM 字符串 题目整理
AC自动机
UVa 11468 Substring
AC自动机+概率DP。
注意要补全不存在的边。
为什么要补全不存在的边呢？补全以后可以直接找到状态的转移，即从所有子节点就可以实现所有状态转移。
#include<iostream>
#include<vector>
#include<cmath>
#include<map>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<queue>
#include<string>
#define INF 1000000000LL
#define ll long long
#define maxnode 4000+5
#define sigma_size 62
using namespace std;
struct Tire
{
],end[maxnode],fail[maxnode];
int root,L;
int newnode()
{
memset(ch[L],,sizeof(ch[L]));
end[L++]=;
;
}
void init()
{
L=;
fail[]=;
end[]=;
root=newnode();
}
int idx(char c)
{
if('a'<=c&&c<='z') return c-'a';
;
;
}
void insert(char* word)
{
;
; word[i]; ++i)
{
int x=idx(word[i]);
if(!ch[now][x])
ch[now][x]=newnode();
now=ch[now][x];
}
end[now]=;
}
void build()
{
queue<int> que;
; i<sigma_size; ++i)
{
][i])
{
que.push(ch[][i]);
fail[ch[][i]]=;
}
}
while(!que.empty())
{
int q=que.front();
que.pop();
; i<sigma_size; ++i)
{
int u=ch[q][i];
if(!u)
{
ch[q][i]=ch[fail[q]][i];
continue;
}
int v=fail[q];
que.push(u);
while(v&&!ch[v][i]) v=fail[v];
fail[u]=ch[v][i];
end[u]|=end[fail[u]];
}
}
}
};
Tire ac;
];
][];
][];
];
double dp(int now,int L)
{
if(!L) return 1.0;
if(vis[now][L]) return f[now][L];
vis[now][L]=true;
f[now][L]=0.0;
; i<sigma_size; ++i)
if(cha[i])
{
int u=ac.ch[now][i];
if(!ac.end[u])
f[now][L]+=pro[i]*dp(u,L-);
}
return f[now][L];
}
int main()
{
;
//freopen("read.txt","r",stdin);
//freopen("out.txt","w",stdout);
scanf("%d",&T);
while(T--)
{
int k;
scanf("%d",&k);
ac.init();
memset(vis,,sizeof(vis));
memset(cha,,sizeof(cha));
; i<k; ++i)
{
];
scanf("%s",word);
ac.insert(word);
}
ac.build();
int n;
scanf("%d",&n);
memset(pro,,sizeof(pro));
; i<n; ++i)
{
];
double p;
scanf("%s%lf",word,&p);
]);
cha[q]=true;
pro[q]=p;
}
int L;
scanf("%d",&L);
,L);
printf("Case #%d: %.6lf\n",++kase,ans);
}
;
}
HDU 2825 Wireless Password
AC自动机+状态压缩DP
求长度为n的字符串中包含超过k个模式串的个数。
dp[i][j][k]表示长为i，当前状态为j时的字符串包括k个串时的个数，k为集合。
状态压缩使每个串唯一，而不必考虑重复。串可以重叠而考虑AC自动机。
另外注意如果当前dp为0，那就不用刷表了。
#include<iostream>
#include<vector>
#include<cmath>
#include<map>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<queue>
#include<string>
#define ll long long
#define maxnode 105
#define sigma_size 26
#define maxn 100
#define MOD 20090717
using namespace std;
int n,m,K;
int countBit(int val)
{
;
while(val)
{
cnt+=val%;
val/=;
}
return cnt;
}
struct Tire
{
int ch[maxnode][sigma_size],end[maxnode],fail[maxnode];
int root,L;
int newnode()
{
memset(ch[L],,sizeof(ch[L]));
end[L++]=;
;
}
void init()
{
L=;
fail[]=;
end[]=;
root=newnode();
}
int idx(char c)
{
return c-'a';
}
void insert(char* word,int p)
{
;
; word[i]; ++i)
{
int x=idx(word[i]);
if(!ch[now][x])
ch[now][x]=newnode();
now=ch[now][x];
}
end[now]=<<p;
}
void build()
{
queue<int> que;
; i<sigma_size; ++i)
{
][i])
{
que.push(ch[][i]);
fail[ch[][i]]=;
}
}
while(!que.empty())
{
int q=que.front();
que.pop();
; i<sigma_size; ++i)
{
int u=ch[q][i];
if(!u)
{
ch[q][i]=ch[fail[q]][i];
continue;
}
int v=fail[q];
que.push(u);
while(v&&!ch[v][i]) v=fail[v];
fail[u]=ch[v][i];
end[u]|=end[fail[u]];
}
}
}
][maxn][(<<)+];
int solve()
{
memset(dp,,sizeof(dp));
dp[][][]=;
; i<n; ++i)
{
; j<L; ++j)
{
; k<(<<m); ++k)
if(dp[i][j][k])
{
; l<; ++l)
{
int u=ch[j][l];
][u][k|end[u]];
res+=dp[i][j][k];
res%=MOD;
}
}
}
}
;
; j<(<<m); ++j)
if(countBit(j)>=K)
{
; i<L; ++i)
{
ans+=dp[n][i][j];
ans%=MOD;
}
}
return ans;
}
};
Tire ac;
int main()
{
while(scanf("%d%d%d",&n,&m,&K)!=EOF)
{
if(!n&&!m&&!K) break;
ac.init();
; i<m; ++i)
{
];
scanf("%s",word);
ac.insert(word,i);
}
ac.build();
printf("%d\n",ac.solve());
}
;
}
POJ 3691 DNA repair
AC自动机+DP
问修改最少字符使得串中不含模式子串。
dp[i][j]表示当前是当串中第j个位置是第i个节点状态时的最优解。
当c不是结束节点的时候，c表示从i转移到的状态
dp[i][j]=min{dp[c][j+1]+1}//如果此时需要修改
dp[i][j]=min{dp[c][j+1]}//如果此时不修改
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <vector>
#include <queue>
#define ll long long
using namespace std;
;
;
;
struct Tire
{
int ch[maxnode][sigma_size];
int end[maxnode],fail[maxnode];
int sz;
int newnode()
{
memset(ch[sz],,sizeof(ch[sz]));
fail[sz]=;
end[sz++]=;
;
}
void init()
{
sz=;
newnode();
}
int idx(char c)
{
;
;
;
;
}
void insert(char *word)
{
;
; word[i]; ++i)
{
int x=idx(word[i]);
if(!ch[now][x])
ch[now][x]=newnode();
now=ch[now][x];
}
end[now]=;
}
void build()
{
queue<int> que;
; i<sigma_size; ++i)
{
][i];
if(u)
{
fail[u]=;
que.push(u);
}
}
while(!que.empty())
{
int q=que.front();
que.pop();
; i<sigma_size; ++i)
{
int u=ch[q][i],v=fail[q];
if(!u)
{
ch[q][i]=ch[v][i];
continue;
}
que.push(u);
while(v&&!ch[v][i]) v=fail[v];
fail[u]=ch[v][i];
end[u]|=end[fail[u]];
}
}
}
};
][];
][];
];
int N;
Tire ac;
int dp(int now,int L)
{
;
if(vis[now][L]) return f[now][L];
vis[now][L]=true;
f[now][L]=inf;
; i<; ++i)
{
int y=ac.idx(str[L]);
int u=ac.ch[now][i];
if(!ac.end[u])
{
)!=inf)
f[now][L]=min(f[now][L],dp(u,L+)+);
)!=inf)
f[now][L]=min(f[now][L],dp(u,L+));
}
}
return f[now][L];
}
int main()
{
;
while(scanf("%d",&n)&&n)
{
];
ac.init();
memset(vis,,sizeof(vis));
; i<=n; ++i)
{
scanf("%s",word);
ac.insert(word);
}
scanf("%s",str);
ac.build();
N=strlen(str);
printf("Case %d: ",++kase);
,);
if(ans==inf) puts("-1");
else printf("%d\n",ans);
}
;
}
字典树
POJ 3630 Phone List
问有没有一个电话号是另一个电话号的前缀。
典型的trie树，时间复杂度是10n。注意不要边读入边判断。否则像12、2这样就判不出来。
#include<cstdio>
#include<vector>
#include<algorithm>
#include<queue>
#include<cstring>
#define LL unsigned int
using namespace std;
;
;
struct Tire
{
int ch[maxnode][sigma_size];
int end[maxnode];
int sz;
int newnode()
{
memset(ch[sz],,sizeof(ch[sz]));
end[sz]=;
return sz++;
}
void init()
{
sz=;
newnode();
}
bool insert(char *word,int v)
{
;
; word[i]; ++i)
{
';
if(!ch[now][x])
ch[now][x]=newnode();
now=ch[now][x];
if(end[now]&&end[now]!=v)
return false;
}
end[now]=v;
return true;
}
};
Tire tree;
][];
int main()
{
int T;
scanf("%d",&T);
while(T--)
{
int n;
scanf("%d",&n);
bool ok=true;
tree.init();
; i<=n; ++i)
{
scanf("%s",numb[i]);
if(ok&&!tree.insert(numb[i],i))
ok=false;
}
; i<=n&&ok; ++i)
{
if(!tree.insert(numb[i],i))
ok=false;
}
if(ok) puts("YES");
else puts("NO");
}
;
}
POJ 1204 Word Puzzles
在一个字符矩阵中寻找出现的字符串，并输出头字符的位置和方向。
可以用字典树来做。枚举每个位置再枚举方向，这样每次遍历一个字符串看是否出现了模式串。
#include<cstdio>
#include<vector>
#include<algorithm>
#include<queue>
#include<cstring>
#define LL unsigned int
using namespace std;
*+;
;
][];
int C,L,W;
bool judge(int x,int y)
{
<=x&&x<L&&<=y&&y<C;
}
int ax,ay,d;
][];
][]= {{-,},{-,},{,},{,},{,},{,-},{,-},{-,-}};
struct Tire
{
int ch[maxnode][sigma_size];
int end[maxnode];
int sz;
int newnode()
{
memset(ch[sz],-,sizeof(ch[sz]));
end[sz]=;
return sz++;
}
void init()
{
sz=;
newnode();
}
int idx(char c)
{
return c-'A';
}
void insert(char *word,int v)
{
;
; word[i]; ++i)
{
int x=idx(word[i]);
)
ch[now][x]=newnode();
now=ch[now][x];
}
end[now]=v;
}
void find(int x,int y)
{
;
)
{
int c=idx(mat[x][y]);
now=ch[now][c];
) break;
if(end[now])
{
ans[end[now]][]=ax;
ans[end[now]][]=ay;
ans[end[now]][]=d;
end[now]=;
}
x+=dir[d][];
y+=dir[d][];
if(!judge(x,y)) break;
}
}
void find(int x,int y,int now)
{
)return;
)
{
ans[end[now]][]=ax;
ans[end[now]][]=ay;
ans[end[now]][]=d;
end[now]=;
}
if(!judge(x,y)) return;
find(x+dir[d][],y+dir[d][],ch[now][idx(mat[x][y])]);
}
};
Tire tree;
int main()
{
while(scanf("%d%d%d",&L,&C,&W)!=EOF)
{
; i<L; ++i)
scanf("%s",mat[i]);
tree.init();
; i<=W; ++i)
{
];
scanf("%s",str);
tree.insert(str,i);
}
; i<L; ++i)
; j<C; ++j)
; k<; ++k)
{
ax=i;
ay=j;
d=k;
tree.find(i,j);
}
; i<=W; ++i)
printf(],ans[i][],ans[i][]+'A');
}
;
}
ACM 字符串 题目整理的更多相关文章
- ACM 矩阵题目整理
先从最基础的矩阵快速幂加速递推开始. HDU 1005 Number Sequence |f[n-2],f[n-1]|* |0 B| =|f[n-1], B*f[n-2]+A*f[n-1]|=|f[n ...
- 【好好补题，因为没准题目还会再出第三遍！！】ACM字符串-组合数学（官方题解是数位DP来写）
ACM字符串 .长度不能超过n .字符串中仅包含大写字母 .生成的字符串必须包含字符串“ACM”,ACM字符串要求连在一块! ok,是不是很简单?现在告诉你n的值,你来告诉我这样的字符串有多少个 输入 ...
- Noip往年题目整理
Noip往年题目整理 张炳琪 一．历年题目 按时间倒序排序 年份 T1知识点 T2知识点 T3知识点 得分 总体 2016day1 模拟 Lca,树上差分 期望dp 144 挺难的一套题目,偏思维难度 ...
- NOIp初赛题目整理
NOIp初赛题目整理 这个 blog 用来整理扶苏准备第一轮 csp 时所做的与 csp 没 有 关 系 的历年 noip-J/S 初赛题目,记录了一些我从不知道的细碎知识点,还有一些憨憨题目,不定期 ...
- ACM - 动态规划专题 题目整理
CodeForces 429B Working out 预处理出从四个顶点到某个位置的最大权值,再枚举相遇点,相遇的时候只有两种情况,取最优解即可. #include<iostream> ...
- ACM 暴力搜索题 题目整理
UVa 129 Krypton Factor 注意输出格式,比较坑爹. 每次要进行处理去掉容易的串,统计困难串的个数. #include<iostream> #include<vec ...
- ACM金牌选手整理的【LeetCode刷题顺序】
算法和数据结构知识点图 首先,了解算法和数据结构有哪些知识点,在后面的学习中有 大局观,对学习和刷题十分有帮助. 下面是我花了一天时间花的算法和数据结构的知识结构,大家可以看看. 后面是为大家 精心挑 ...
- BZOJ 题目整理
bzoj 500题纪念 总结一发题目吧,挑几道题整理一下,(方便拖板子) 1039:每条线段与前一条线段之间的长度的比例和夹角不会因平移.旋转.放缩而改变,所以将每条轨迹改为比例和夹角的序列,复制一份 ...
- C和C++字符串处理整理
在刷leetcode题目的过程中,发现自己对于c和c++字符串的处理并不是很拿手,处理起来比较费劲,而且,算法题似乎很中意字符串的处理,有很多题目都涉及到它.字符串处理比较基础,但是很重要,因此,整理 ...
随机推荐
- 将项目上传到git上，并在测试服务器上运行
上周我完成了我人生的第一个项目,并将它上传到gitHub上完成了团队合作.下面我将用倒叙的方法记录下,我的这段经历. 星期五,项目完成准备上传gitHub启动测试服务器. 1.git clone 项目 ...
- 数字图像处理作业使用OpenCV - 自定义直方图
第二次作业需要打印出来灰度直方图,当然不能使用ocv的自带calcHist函数来得到Mat对象了……结果上网搜索怎么用自己的数据创建直方图,搜到的都是直接用函数的_(:з」∠)_ 结果这个地方拖了好久 ...
- SpringMVC拦截器（资源和权限管理）
1.DispatcherServlet SpringMVC具有统一的入口DispatcherServlet,所有的请求都通过DispatcherServlet. DispatcherServle ...
- Jquery动态操作checkbox
问题的由来:html页面中使用checkbox给新建用户授权,然后提交到后台服务器,完成给用户授权.用户授权完成之后,如果用户的权限发生改变管理员需要修改多个用户的权限或查看某用户的权限,点击不同的用 ...
- Takeown--夺取文件or文件夹所有权
强制将当前目录下的所有文件及文件夹.子文件夹下的所有者更改为管理员组(administrators)命令:takeown /f * /a /r /d y 将所有d:\documents目录下的文件.子 ...
- 深入springMVC------文件上传源码解析(上篇)
最近在项目中,使用springmvc 进行上传文件时,出现了一个问题: org.springframework.web.multipart.MultipartException: The curren ...
- 使用SecureCRT远程链接Ubuntu出现 Change of username or service not allowed的问题
RT: 首先是确认ubuntu上有运行 sshd服务的 但是用SecureCRT链接时报错,默认用户名是root: 打开ssh_config检查下是否禁止直接用root登陆 这句 改成yes o ...
- 高性能 Oracle JDBC 编程
了解如何利用连接和语句池特性来提高 Oracle 驱动的 JDBC 程序的性能.作者:Yuli Vasiliev2009 年 4 月发布使用诸如连接池和语句池等池技术可以显著提高数据库密集型应用程序的 ...
- MSDN Library for vs 2010安装及使用(MSDN Library)
VS2010正式版不再有单独的MSDN Library安装选项,VS2010的ISO安装光盘里已经包含有MSDN Library,只不过要手动安装,方法如下: 1.安装完VS2010后,在开始菜单中打 ...
- Eclipse JUnit 生成报告
http://blog.sina.com.cn/s/blog_8af106960102v6qh.html 对Eclipse的工程写单元测试: 第一步: 1. 一个工程有多个测试类,将测试类放到一个测试 ...