博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
POJ 1226 Substrings
阅读量:7193 次
发布时间:2019-06-29

本文共 2687 字,大约阅读时间需要 8 分钟。

POJ_1226

    这个题目可以枚举长度最小的字符串的所有子串,然后暴搜。但为了练一练后缀数组,所以就用后缀数组写了。

#include
#include
#define MAXD 20100 char b[110]; int N, M, h[110], f[MAXD], r[MAXD], rank[MAXD], sa[MAXD], height[MAXD], ws[MAXD], wv[MAXD], wa[MAXD], wb[MAXD]; void init() {
int i, j, k, x = 128; k = 0; for(i = 0; i < N; i ++) {
scanf("%s", b); for(j = 0; b[j]; j ++) {
r[k] = b[j]; f[k] = i; ++ k; } r[k ++] = x ++; for(-- j; j >= 0; j --) {
r[k] = b[j]; f[k] = i; ++ k; } r[k ++] = x ++; } M = k - 1; r[M] = 0; } int cmp(int *p, int x, int y, int l) {
return p[x] == p[y] && p[x + l] == p[y + l]; } void da(int n, int m) {
int i, j, p, *x = wa, *y = wb, *t; for(i = 0; i < m; i ++) ws[i] = 0; for(i = 0; i < n; i ++) ++ ws[x[i] = r[i]]; for(i = 1; i < m; i ++) ws[i] += ws[i - 1]; for(i = n - 1; i >= 0; i --) sa[-- ws[x[i]]] = i; for(j = p = 1; p < n; j *= 2, m = p) {
for(p = 0, i = n - j; i < n; i ++) y[p ++] = i; for(i = 0; i < n; i ++) if(sa[i] >= j) y[p ++] = sa[i] - j; for(i = 0; i < n; i ++) wv[i] = x[y[i]]; for(i = 0; i < m; i ++) ws[i] = 0; for(i = 0; i < n; i ++) ++ ws[wv[i]]; for(i = 1; i < m; i ++) ws[i] += ws[i - 1]; for(i = n - 1; i >= 0; i --) sa[-- ws[wv[i]]] = y[i]; for(t = x, x = y, y = t, x[sa[0]] = 0, i = p = 1; i < n; i ++) x[sa[i]] = cmp(y, sa[i - 1], sa[i], j) ? p - 1 : p ++; } } void calheight(int n) {
int i, j, k = 0; for(i = 1; i <= n; i ++) rank[sa[i]] = i; for(i = 0; i < n; height[rank[i ++]] = k) for(k ? -- k : 0, j = sa[rank[i] - 1]; r[i + k] == r[j + k]; k ++); } int check(int mid) {
int i, j, k = 0; for(i = 1; i <= M; i ++) {
if(height[i] < mid) {
if(k && k == N) return 1; k = 0; } else {
if(k == 0) {
memset(h, 0, sizeof(h)); h[f[sa[i - 1]]] = 1; ++ k; } if(!h[f[sa[i]]]) {
h[f[sa[i]]] = 1; ++ k; } } } return 0; } void solve() {
int i, j, k, min, max, mid; da(M + 1, 400); calheight(M); min = 0, max = 100; for(;;) {
mid = (max + min) / 2; if(mid == min) break; if(check(mid)) min = mid; else max = mid; } printf("%d\n", mid); } int main() {
int t; scanf("%d", &t); while(t --) {
scanf("%d", &N); init(); solve(); } return 0; }

转载地址:http://vvtkm.baihongyu.com/

你可能感兴趣的文章
poj 2154 Color——带优化的置换
查看>>
Yii 框架里数据库操作详解-[增加、查询、更新、删除的方法 'AR模式'](转)
查看>>
asp.net 服务器 上传文件到 FTP服务器
查看>>
ZooKeeper学习第一期---Zookeeper简单介绍
查看>>
Mybatis(spring)(多个参数)(插入数据返回id)
查看>>
操作系统学习笔记三 进程
查看>>
Map的嵌套使用
查看>>
实习公司的开发环境搭建教程
查看>>
Linux busybox mount -a fstab
查看>>
Dijkstra
查看>>
property测试代码:
查看>>
[C# | WinCE | Solution] 在 WinCE 上访问 SSL 加密后的 WCF SOAP 服务接口出现“未能与远程服务器建立信任关系”...
查看>>
css2D转换和3D转换
查看>>
synchronized关键字的用法总结
查看>>
用C实现一个简单的对拍器——致每个曾经为求AC披星戴月的程序员们
查看>>
C++的const类成员函数
查看>>
js 事件驱动机制
查看>>
进程相关内容
查看>>
框架,样式
查看>>
SqlServer2005通过出生日期算年龄函数
查看>>