uvalive4513
https://vjudge.net/problem/UVALive-4513
终于做出来了。。。。。。。。。
各种sb错误,最后对拍出来了,还没改对。。。。。。。。。。。。。。。。。。。。。
快半天的时间。。。。。。。。。。
看上去第一眼就猜出了一个结论:这些相同的字符串肯定是按后缀排序的顺序存在的。
于是就想出了一个方法:求出sa和lcp,然后扫描每个长为m的区间,rmq其中最小的lcp和最大的sa作为答案(扫描是指枚举左端点,右端点是左端点+m-2)
这样复杂度nlogn,可以接受,用单调队列?(不清楚,没学好),于是我用了线段树。
(请无视)然后各种wa,各种改,搞得我怀疑自己的方法是错的,然后看了看题解,好像都是二分判定。再想了想,貌似是对的,但是不确定,就翻了翻别人的代码,看见和自己的想法差不多,就确定是对的了,然后就各种调试,看见一个又一个sb错误:update1和update搞反找了半天。。。然后又对拍拍了半天拍出来了组错误数据,发现sa[0]=0这样的串会错,于是就改掉,对了。。。
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define N 80010
struct data
{
int val,pos;
};
int n,m,k;
char s[N];
int Rank[N],sa[N],temp[N],lcp[N],tree[N<<],tree1[N<<];
bool cp(int i,int j)
{
if(Rank[i]!=Rank[j]) return Rank[i]<Rank[j];
int ri=i+k<n?Rank[i+k]:-;
int rj=j+k<n?Rank[j+k]:-;
return ri<rj;
}
void update(int l,int r,int x,int pos,int num)
{
if(l==r)
{
tree[x]=num;
return;
}
if(pos>(l+r)/) update((l+r)/+,r,x*+,pos,num);
else update(l,(l+r)/,x*,pos,num);
tree[x]=min(tree[x*],tree[x*+]);
}
int query(int l,int r,int x,int a,int b)
{
if(l>b||r<a) return N;
if(l>=a&&r<=b) return tree[x];
int i=query(l,(l+r)/,x*,a,b);
int j=query((l+r)/+,r,x*+,a,b);
return min(i,j);
}
void update1(int l,int r,int x,int pos,int Pos)
{
if(l==r)
{
tree1[x]=Pos;
return;
}
if(pos>(l+r)/) update1((l+r)/+,r,x*+,pos,Pos);
else update1(l,(l+r)/,x*,pos,Pos);
tree1[x]=max(tree1[x*],tree1[x*+]);
}
int query1(int l,int r,int x,int a,int b)
{
if(l>b||r<a) return -N;
if(l>=a&&r<=b) return tree1[x];
int i=query1(l,(l+r)/,x*,a,b);
int j=query1((l+r)/+,r,x*+,a,b);
return max(i,j);
}
void Sa()
{
n=strlen(s);
for(int i=;i<n;i++)
{
Rank[i]=s[i]; sa[i]=i;
}
for(k=;k<=n;k*=)
{
sort(sa,sa+n,cp);
temp[sa[]]=;
for(int i=;i<n;i++) temp[sa[i]]=temp[sa[i-]]+(cp(sa[i-],sa[i]));
for(int i=;i<n;i++) Rank[i]=temp[i];
}
}
void Lcp()
{
n=strlen(s);
for(int i=;i<n;i++) Rank[sa[i]]=i;
int h=;
for(int i=;i<n;i++)
{
int j=sa[Rank[i]-];
if(Rank[i]<=) continue;
if(h>) h--;
for(;((i+h<n)&&(j+h<n));h++)
{
if(s[i+h]!=s[j+h]) break;
}
lcp[Rank[i]-]=h;
}
update1(,n,,,sa[]);
for(int i=;i<n-;i++)
{
update(,n,,i+,lcp[i]);
update1(,n,,i+,sa[i+]);
}
}
void solve()
{
memset(sa,,sizeof(sa));
memset(lcp,,sizeof(lcp));
memset(tree,0x3f3f,sizeof(tree));
memset(tree1,,sizeof(tree1));
memset(Rank,,sizeof(Rank));
Sa(); Lcp();
int l=;
data ans; ans.val=; ans.pos=;
for(int i=;i+m-<n;i++)
{
data x;
x.val=query(,n,,i,i+m-);
x.pos=query1(,n,,i,i+m-);
if(x.val>ans.val||(x.val==ans.val&&x.pos>ans.pos))
{
ans.val=x.val;
ans.pos=x.pos;
}
}
if(ans.val==) printf("none\n");
else printf("%d %d\n",ans.val,ans.pos);
}
int main()
{
while()
{
scanf("%d",&m);
if(!m) break;
scanf("%s",s);
if(m==)
{
printf("%d %d\n",strlen(s),);
continue;
}
n=strlen(s);
solve();
}
return ;
}
uvalive4513的更多相关文章
- UVALive4513 Stammering Aliens(哈希法,后缀数组)
题目链接:http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=12580 [思路] 求出现次数不小于k次的最长可重叠子串和最后的出现 ...
- HASH算法小结
一.简述 HASH算法的本质是特征提取——将某种不太好表示的特征,通过某种压缩的方式映射成一个值.这样,就可以优雅解决一部分难以解决的特征统计问题. 同时考虑到hash算法的本质是个概率算法,因此并不 ...
随机推荐
- POJ 3280 Cheapest Palindrome(DP 回文变形)
题目链接:http://poj.org/problem?id=3280 题目大意:给定一个字符串,可以删除增加,每个操作都有代价,求出将字符串转换成回文串的最小代价 Sample Input 3 4 ...
- 关于php的一些开源程序
最好用的当属thinksns,目前更新到4.3.4, 社交型网站. 此网站提供大量源码,有时间可以去看看:http://down.admin5.com/
- 关于谷歌浏览器下自动填写密码的bug
有的时候我们需要用到h5中input的新类型,type=password 这个可以很好的将输入的内容变成圆点,但是这样做也有一个小bug就是在某些我们不希望自动保存密码的页面中,也会出现浏览器自动帮我 ...
- 添加标签2 jquery 和JS
TAG添加标签 做了个方法方便调用 一.JS版本 <!DOCTYPE html> <html lang="en"> <head> <met ...
- angularjs入门整理
之前发过一篇博文,从mobile angular ui的demo和其官网初识整个angularjs的大体使用,但是没很好学习,只是通过一些技术博文初步认识,陷入很多坑.所以现在在中文官网正式整理下知识 ...
- python【第十二篇】Mysql基础
内容: 1.数据库介绍及MySQL简介 2.MySQL基本操作 1 数据库介绍 1.1什么是数据库? 数据库(Database)是按照数据结构来组织.存储和管理数据的仓库,每个数据库都有一个或多个不同 ...
- Junit 源码剖析(二)
junit4 下的所有的testcase都是在Runner下执行的, 可以将Runner理解为junit运行的容器, 默认情况下junit会使用JUnit4ClassRunner作为所有testcas ...
- kernel编程规范
1. 制表符8个空格2. 每行最长80字符3. 代码块的{放在首行,但是函数的{放在次行4. 只有一行的if块,不加{}5. 不在()前后加空格6. 正常关键字后加一个空格,if, switch, c ...
- JQUERY 插件开发——LAZYLOADIMG(预加载和延迟加载图片)
开发背景 本插件开发是近期写的最后一个插件了,接下来我想把最近研究的redis最为一个系列阐述下.当然Jquery插件开发是我个人爱好,我不会停止,在将来的开发中我会继续完善,当然也会坚持写这个系列的 ...
- [BZOJ 1559] [JSOI2009] 密码 【AC自动机DP】
题目链接:BZOJ - 1559 题目分析 将给定的串建成AC自动机,然后在AC自动机上状压DP. 转移边就是Father -> Son 或 Now -> Fail. f[i][j][k] ...