BZOJ 1014: [JSOI2008]火星人prefix Splay+二分
1014: [JSOI2008]火星人prefix
题目连接:
http://www.lydsy.com/JudgeOnline/problem.php?id=1014
Description
火星人最近研究了一种操作:求一个字串两个后缀的公共前缀。比方说,有这样一个字符串:madamimadam,我们将这个字符串的各个字符予以标号:序号: 1 2 3 4 5 6 7 8 9 10 11 字符 m a d a m i m a d a m 现在,火星人定义了一个函数LCQ(x, y),表示:该字符串中第x个字符开始的字串,与该字符串中第y个字符开始的字串,两个字串的公共前缀的长度。比方说,LCQ(1, 7) = 5, LCQ(2, 10) = 1, LCQ(4, 7) = 0 在研究LCQ函数的过程中,火星人发现了这样的一个关联:如果把该字符串的所有后缀排好序,就可以很快地求出LCQ函数的值;同样,如果求出了LCQ函数的值,也可以很快地将该字符串的后缀排好序。 尽管火星人聪明地找到了求取LCQ函数的快速算法,但不甘心认输的地球人又给火星人出了个难题:在求取LCQ函数的同时,还可以改变字符串本身。具体地说,可以更改字符串中某一个字符的值,也可以在字符串中的某一个位置插入一个字符。地球人想考验一下,在如此复杂的问题中,火星人是否还能够做到很快地求取LCQ函数的值。
Input
第一行给出初始的字符串。第二行是一个非负整数M,表示操作的个数。接下来的M行,每行描述一个操作。操作有3种,如下所示: 1、 询问。语法:Q x y,x, y均为正整数。功能:计算LCQ(x, y) 限制:1 <= x, y <= 当前字符串长度。 2、 修改。语法:R x d,x是正整数,d是字符。功能:将字符串中第x个数修改为字符d。限制:x不超过当前字符串长度。 3、 插入:语法:I x d,x是非负整数,d是字符。功能:在字符串第x个字符之后插入字符d,如果x = 0,则在字符串开头插入。限制:x不超过当前字符串长度。
Output
对于输入文件中每一个询问操作,你都应该输出对应的答案。一个答案一行。
Sample Input
madamimadam
7
Q 1 7
Q 4 8
Q 10 11
R 3 a
Q 1 7
I 10 a
Q 2 11
Sample Output
5
1
0
2
1
Hint
数据规模:
对于100%的数据,满足:
1、 所有字符串自始至终都只有小写字母构成。
2、 M <= 150,000
3、 字符串长度L自始至终都满足L <= 100,000
4、 询问操作的个数不超过10,000个。
对于第1,2个数据,字符串长度自始至终都不超过1,000
对于第3,4,5个数据,没有插入操作。
题意
题解:
用splay去维护就好了
询问就直接暴力二分判断就行了
代码
#include<bits/stdc++.h>
using namespace std;
const int maxn = 1e6+7;
const int mod = 1e9+7;
const int P = 27;
struct Splay
{
int fa[maxn];
int ch[maxn][2];
int val[maxn];
int siz[maxn];
long long h[maxn];
long long p[maxn];
int root,tot,n;
char s[maxn];
void newnode(int &r,int father,int k)//建立新的节点
{
r=++tot;
fa[r]=father;
val[r]=k;
ch[r][0]=ch[r][1]=0;
siz[r]=1;
}
void push_up(int x)
{
int l=ch[x][0],r=ch[x][1];
siz[x]=siz[l]+siz[r]+1;
h[x]=h[l]*p[siz[r]+1]+val[x]*p[siz[r]]+h[r];
h[x]%=mod;
}
void zig(int x)
{
int y=fa[x];
ch[y][0]=ch[x][1];
fa[ch[x][1]]=y;
ch[x][1]=y;
fa[x]=fa[y];
if(y==ch[fa[y]][0])
ch[fa[y]][0]=x;
else if(y==ch[fa[y]][1])
ch[fa[y]][1]=x;
fa[y]=x;
push_up(y);
if(y==root)root=x;
}
void zag(int x)
{
int y=fa[x];
ch[y][1]=ch[x][0];
fa[ch[x][0]]=y;
ch[x][0]=y;
fa[x]=fa[y];
if(y==ch[fa[y]][0])
ch[fa[y]][0]=x;
else if(y==ch[fa[y]][1])
ch[fa[y]][1]=x;
fa[y]=x;
push_up(y);
if(y==root)root=x;
}
void splay(int x,int goal)
{
while(fa[x]!=goal)
{
int y=fa[x],z=fa[y];
if(z==goal)
{
if(x==ch[y][0])
zig(x);
else
zag(x);
break;
}
if(x==ch[y][0])
{
if(y==ch[z][0])
zig(y);
zig(x);
}
else
{
if(y==ch[z][1])
zag(y);
zag(x);
}
}
push_up(x);
}
void fi(int x,int rk,int z)
{
while(1)
{
if(rk<=siz[ch[x][0]])
x=ch[x][0];
else
{
rk-=siz[ch[x][0]];
if(rk==1)break;
rk--;
x=ch[x][1];
}
}
splay(x,z);
}
void build(int &x,int l,int r,int father)
{
if(l>r)return;
int mid = (l+r)/2;
newnode(x,father,s[mid]-'a'+1);
if(l<mid)build(ch[x][0],l,mid-1,x);
if(r>mid)build(ch[x][1],mid+1,r,x);
push_up(x);
}
void init()
{
p[0]=1;
for(int i=1;i<maxn;i++)
p[i]=p[i-1]*P%mod;
scanf("%s",s);
n=strlen(s);
newnode(root,0,1);
newnode(ch[root][1],root,1);
siz[root]=2;
build(ch[ch[root][1]][0],0,n-1,ch[root][1]);
push_up(ch[root][1]);
push_up(root);
}
long long query(int x,int len)
{
fi(root,x,0);
fi(root,x+len+1,root);
//if(len==1)return val[ch[ch[root][1]][0]];
//printf("%c %c %lld\n",val[ch[ch[root][1]][0]]+'a'-1,val[ch[root][1]]+'a'-1,h[ch[ch[root][1]][0]]);
return h[ch[ch[root][1]][0]];
}
int solve(int x,int y)
{
int l=1,r=min(siz[root]-x,siz[root]-y)-1;
int ans=0;
while(l<=r)
{
int mid = (l+r)/2;
if(query(x,mid)==query(y,mid))l=mid+1,ans=mid;
else r=mid-1;
}
return ans;
}
void change(int x,int y)
{
fi(root,x+1,0);
val[root]=y;
push_up(root);
}
void in(int x,int y)
{
fi(root,x+1,0);
fi(root,x+2,root);
newnode(ch[ch[root][1]][0],ch[root][1],y);
push_up(ch[ch[root][1]][0]);
push_up(ch[root][1]);
push_up(root);
}
}T;
int main()
{
T.init();
int q;scanf("%d",&q);
while(q--)
{
char c[10];
scanf("%s",c);
if(c[0]=='R')
{
int x;char k[10];
scanf("%d%s",&x,k);
T.change(x,k[0]-'a'+1);
}
else if(c[0]=='I')
{
int x;char k[10];
scanf("%d%s",&x,k);
T.in(x,k[0]-'a'+1);
}
else
{
int x,y;scanf("%d%d",&x,&y);
printf("%d\n",T.solve(x,y));
}
}
}
BZOJ 1014: [JSOI2008]火星人prefix Splay+二分的更多相关文章
- BZOJ 1014: [JSOI2008]火星人prefix [splay 二分+hash] 【未完】
1014: [JSOI2008]火星人prefix Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 6243 Solved: 2007[Submit] ...
- BZOJ 1014 [JSOI2008]火星人prefix (splay+二分答案+字符串hash)
题目大意:维护一个字符串,支持插入字符和替换字符的操作,以及查询该字符串两个后缀的最长公共前缀长度 乍一看以为是后缀数组,然而并没有可持久化后缀数组(雾) 看题解才知道这是一道splay题,首先要对s ...
- BZOJ 1014 [JSOI2008]火星人prefix (Splay + Hash + 二分)
1014: [JSOI2008]火星人prefix Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 8112 Solved: 2569[Submit] ...
- BZOJ 1014: [JSOI2008]火星人prefix( splay + hash )
用splay维护序列, 二分+hash来判断LCQ.. #include<bits/stdc++.h> using namespace std; typedef unsigned long ...
- bzoj 1014 [JSOI2008]火星人prefix——splay+哈希
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=1014 用splay维护字符串,每个点记录子树的哈希值,然后二分查询. 二分不是把两个点的哈希 ...
- BZOJ 1014 [JSOI2008]火星人prefix | Splay维护哈希值
题目: 题解: #include<cstdio> #include<algorithm> #include<cstring> typedef long long l ...
- bzoj 1014: [JSOI2008]火星人prefix hash && splay
1014: [JSOI2008]火星人prefix Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 3154 Solved: 948[Submit][ ...
- 求帮看!!!!BZOJ 1014 [JSOI2008]火星人prefix
1014: [JSOI2008]火星人prefix Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 4164 Solved: 1277[Submit] ...
- [BZOJ 1014] [JSOI2008] 火星人prefix 【Splay + Hash】
题目链接:BZOJ - 1014 题目分析 求两个串的 LCP ,一种常见的方法就是 二分+Hash,对于一个二分的长度 l,如果两个串的长度为 l 的前缀的Hash相等,就认为他们相等. 这里有修改 ...
随机推荐
- WebAPI初探
由于即将要接手的新项目计划用ASP.NET MVC3来开发,所以最近一段时间一直在看相关的书或文章.因为之前在大学里也曾学习过MVC2开发,也做过几个简单的MVC2的小型测试项目,不过在后来工作以后主 ...
- Delphi 7使用自定义图标关联文件类型
Delphi 7使用自定义图标关联文件类型 5.2 Delphi编程(40) 版权声明:本文为博主原创文章,未经博主允许不得转载. 在开发过程中,我们经常需要属于自己的文件类型,自定义的后缀名不仅可 ...
- 使用Drush管理Drupal站点
Drush(Drush = Drupal + Shell)就是使用命令行命令来操作Drupal站点,它的命令格式与git类似,都是双字命令(drush + 实际的命令).既然是命令行命令,也就可以使用 ...
- fopen警告处理方式
warning C4996: “fopen”被声明为否决的 问题:vs2005中编程时,遇到如下: warning C4996: “fopen”被声明为否决的. 解释:微软的警告,主要是那些都是C库的 ...
- jQuery文档加载完毕的几种写法
js中文档加载完毕.一般在body加一个onload事件或者window.onload = function () {} jQuery中有好多写法,平时也不注意,别人一问,还真觉得头大. 下面是我整理 ...
- PHP正则表达式匹配中文字符
网上有很多类似的文章,但往往都不能用 所以记录一下 preg_match_all("/([\x{4e00}-\x{9fa5}])/u", $input, $match); 注意:限 ...
- iOS tableview 优化总结
根据网络上的优化方法进行了总括.并未仔细进行语言组织.正在这些优化方法进行学习,见另一篇文章 提高app流畅度 1.cell子控件创建写在 initWithStyle:reuseIdentifier ...
- iOS OpenCV 缺少64位解决方法
- binarySearch二分查找——Javascript实现
在很早之前,我就写过了一篇也关于二分法的相关博文:JavaScript快排与原生sort的测试.当时是用二分法进行快速排序,其实和这次思路大致相当.二分查找最重要的一个条件,就是需要将数组先按照从小到 ...
- (转载)Java里快如闪电的线程间通讯
转自(http://www.infoq.com/cn/articles/High-Performance-Java-Inter-Thread-Communications) 这个故事源自一个很简单的想 ...