BZOJ 4556 [Tjoi2016&Heoi2016]字符串 ——后缀数组 ST表 主席树 二分答案
Solution 1:
后缀数组暴力大法好
#include <map>
#include <cmath>
#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define F(i,j,k) for (int i=j;i<=k;++i)
#define D(i,j,k) for (int i=j;i>=k;--i)
#define ll long long
#define mp make_pair
#define maxn 300005 struct Suffix_Array{
int s[maxn];
int tmp[maxn],rk[maxn],sa[maxn],cnt[maxn],h[maxn];
void build(int n,int m)
{
int i,j,k; n++;
F(i,0,2*n+5) tmp[i]=rk[i]=sa[i]=h[i]=0;
F(i,0,m-1) cnt[i]=0;
F(i,0,n-1) cnt[rk[i]=s[i]]++;
F(i,1,m-1) cnt[i]+=cnt[i-1];
F(i,0,n-1) sa[--cnt[rk[i]]]=i;
for (int k=1;k<=n;k<<=1)
{
F(i,0,n-1){j=sa[i]-k;if(j<0)j+=n;tmp[cnt[rk[j]]++]=j;}
sa[tmp[cnt[0]=0]]=j=0;
F(i,1,n-1)
{
if (rk[tmp[i]]!=rk[tmp[i-1]]||rk[tmp[i]+k]!=rk[tmp[i-1]+k]) cnt[++j]=i;
sa[tmp[i]]=j;
}
memcpy(rk,sa,(n+1)*sizeof (int)); memcpy(sa,tmp,(n+1)*sizeof (int));
if (j>=n-1) break;
}
for (i=k=0;i<n;h[rk[i++]]=k)
for (k?k--:0,j=sa[rk[i]-1];s[i+k]==s[j+k];k++);
}
void work(int n,int m)
{
F(t,1,m)
{
int a,b,c,d;scanf("%d%d%d%d",&a,&b,&c,&d); a--;b--;c--;d--;
int mx=0,mn=d-c+1,mid=rk[c];
if (a<=c&&b>=c) mx=min(d-c+1,b-c+1);
for (int i=mid;i>1;--i)
{
if (mn<=mx) break;
mn=min(mn,h[i]);
if (sa[i-1]>=a&&sa[i-1]<=b)
mx=max(mx,min(mn,b-sa[i-1]+1));
}
mn=d-c+1;
for (int i=mid+1;i<=n;++i)
{
if (mn<=mx) break;
mn=min(mn,h[i]);
if (sa[i]>=a&&sa[i]<=b)
mx=max(mx,min(mn,b-sa[i]+1));
}
printf("%d\n",mx);
}
}
}SA; int n,m;
char s[maxn]; int main()
{
scanf("%d%d",&n,&m);
scanf("%s",s);
F(i,0,n-1) SA.s[i]=s[i]-'a'+1; SA.s[n]=0;
SA.build(n,30); SA.work(n,m);
}
Solution 2:
后缀数组 二分答案 主席数 ST表
每次询问二分答案,然后找出要匹配的串在SA中最左以及最右的位置,然后主席树判断即可。这样貌似是两个$\log$
可以在主席树上直接找前驱后继,然后ST表直接查询,然后就成了一个$\log$
我比较菜,写的是第一种。不过一次写对还是挺欣慰的,(废话,你慢慢写了3h)
#include <map>
#include <cmath>
#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define F(i,j,k) for (int i=j;i<=k;++i)
#define D(i,j,k) for (int i=j;i>=k;--i)
#define mp make_pair
#define maxn 200005 namespace SA{
int tmp[maxn],s[maxn],cnt[maxn],rk[maxn],sa[maxn],h[maxn];
int st[maxn][21],_log[maxn];
void build(int n,int m)
{
int i,j,k; n++;
F(i,0,2*n+1) tmp[i]=rk[i]=sa[i]=h[i]=0;
F(i,0,m-1) cnt[i]=0;
F(i,0,n-1) cnt[rk[i]=s[i]]++;
F(i,1,m-1) cnt[i]+=cnt[i-1];
F(i,0,n-1) sa[--cnt[rk[i]]]=i;
for (int k=1;k<=n;k<<=1)
{
F(i,0,n-1){j=sa[i]-k;if(j<0)j+=n;tmp[cnt[rk[j]]++]=j;}
sa[tmp[cnt[0]=0]]=j=0;
F(i,1,n-1)
{
if (rk[tmp[i]]!=rk[tmp[i-1]]||rk[tmp[i]+k]!=rk[tmp[i-1]+k]) cnt[++j]=i;
sa[tmp[i]]=j;
}
memcpy(rk,sa,(n+1)*sizeof(int)); memcpy(sa,tmp,(n+1)*sizeof(int));
if (j>=n-1) break;
}
for (i=k=0;i<n;h[rk[i++]]=k) for (k?k--:0,j=sa[rk[i]-1];s[i+k]==s[j+k];k++);
F(i,1,n-1) st[i][0]=h[i];
F(i,2,n-1) _log[i]=_log[i>>1]+1;
F(i,1,20) F(j,1,n-(1<<i)) st[j][i]=min(st[j][i-1],st[j+(1<<i-1)][i-1]);
}
int query(int a,int b,int n)
{
if (a==b) return n-sa[a]; a++;
int tmp=_log[b-a+1];
return min(st[a][tmp],st[b-(1<<tmp)+1][tmp]);
}
int lcp(int a,int b,int n)
{
a=rk[a],b=rk[b];
if (a>b) swap(a,b); if (a==b) return n-sa[a];
a++; int tmp=_log[b-a+1];
return min(st[a][tmp],st[b-(1<<tmp)+1][tmp]);
}
} namespace PT{
int ls[maxn<<4],rs[maxn<<4],sum[maxn<<4],rt[maxn],tot=0;
void modify(int o1,int & o2,int l,int r,int X,int f)
{
o2=++tot;sum[o2]=sum[o1]+f;if (l==r) return ;int mid=l+r>>1;
if (X<=mid) rs[o2]=rs[o1],modify(ls[o1],ls[o2],l,mid,X,f);
else ls[o2]=ls[o1],modify(rs[o1],rs[o2],mid+1,r,X,f);
}
int query(int o1,int o2,int l,int r,int L,int R)
{
if (L<=l&&r<=R) return sum[o2]-sum[o1];
int mid=l+r>>1,ret=0;
if (L<=mid) ret+=query(ls[o1],ls[o2],l,mid,L,R);
if (R>mid) ret+=query(rs[o1],rs[o2],mid+1,r,L,R);
return ret;
}
} int n,m;char s[maxn]; bool check(int l,int r,int x,int mid)
{
int pos=SA::rk[x]; //printf("The Postion is %d\n",pos);
int ll=1,rr=pos,posl,posr;
while (ll<rr)
{
int mmiidd=(ll+rr)/2;
if (SA::query(mmiidd,pos,n)>=mid) rr=mmiidd;
else ll=mmiidd+1;
}
posl=rr;
ll=pos,rr=n;
while (ll<rr)
{
int mmiidd=(ll+rr)/2+1;
if (SA::query(pos,mmiidd,n)>=mid) ll=mmiidd;
else rr=mmiidd-1;
}
posr=ll;
if (PT::query(PT::rt[l-1],PT::rt[r],1,n,posl,posr)) return true;
else return false;
} int main()
{
scanf("%d%d",&n,&m);
scanf("%s",s); F(i,0,n-1) SA::s[i]=s[i]-'a'+1; SA::s[n]=0;
SA::build(n,30);
F(i,0,n-1)
{
PT::modify(PT::rt[i-1],PT::rt[i],1,n,SA::rk[i],1);
}
F(i,1,m)
{
int a,b,c,d,mx;
scanf("%d%d%d%d",&a,&b,&c,&d);
a--;b--;c--;d--;
if (a<=c&&c<=b) mx=min(b-c+1,d-c+1); else mx=0;
int l=mx,r=min(d-c+1,b-a+1);
while (l<r)
{
int mid=(l+r)/2+1;
if (check(a,b-mid+1,c,mid)) l=mid;
else r=mid-1;
}
printf("%d\n",l);
}
}
BZOJ 4556 [Tjoi2016&Heoi2016]字符串 ——后缀数组 ST表 主席树 二分答案的更多相关文章
- bzoj 4556 [Tjoi2016&Heoi2016]字符串——后缀数组+主席树
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4556 本来只要查 ht[ ] 数组上的前驱和后继就行,但有长度的限制.可以二分答案解决!然后 ...
- Bzoj 4556: [Tjoi2016&Heoi2016]字符串
4556: [Tjoi2016&Heoi2016]字符串 Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 177 Solved: 92[Sub ...
- Bzoj4556: [Tjoi2016&Heoi2016]字符串 后缀数组
4556: [Tjoi2016&Heoi2016]字符串 Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 169 Solved: 87[Sub ...
- [BZOJ4556][Tjoi2016&Heoi2016]字符串 后缀数组+主席树
4556: [Tjoi2016&Heoi2016]字符串 Time Limit: 20 Sec Memory Limit: 128 MB Description 佳媛姐姐过生日的时候,她的小 ...
- 【BZOJ4556】[Tjoi2016&Heoi2016]字符串 后缀数组+二分+主席树+RMQ
[BZOJ4556][Tjoi2016&Heoi2016]字符串 Description 佳媛姐姐过生日的时候,她的小伙伴从某东上买了一个生日礼物.生日礼物放在一个神奇的箱子中.箱子外边写了一 ...
- BZOJ 4556: [Tjoi2016&Heoi2016]字符串(后缀数组 + 二分答案 + 主席树 + ST表 or 后缀数组 + 暴力)
题意 一个长为 \(n\) 的字符串 \(s\),和 \(m\) 个询问.每次询问有 \(4\) 个参数分别为 \(a,b,c,d\). 要你告诉它 \(s[a...b]\) 中的所有子串 和 \(s ...
- ●BZOJ 4556 [Tjoi2016&Heoi2016]字符串
题链: http://www.lydsy.com/JudgeOnline/problem.php?id=4556 题解: 巨恶心...但是题很好呀,可以练习好几个比较麻烦的算法~ 1).预处理 首先用 ...
- bzoj4556: [Tjoi2016&Heoi2016]字符串 (后缀数组加主席树)
题目是给出一个字符串,每次询问一个区间[a,b]中所有的子串和另一个区间[c,d]的lcp最大值,首先求出后缀数组,对于lcp的最大值肯定是rank[c]的前驱和后继,但是对于这个题会出现问题,就是题 ...
- BZOJ 2119: 股市的预测 [后缀数组 ST表]
2119: 股市的预测 Time Limit: 10 Sec Memory Limit: 259 MBSubmit: 331 Solved: 153[Submit][Status][Discuss ...
随机推荐
- UVA Planning mobile robot on Tree树上的机器人(状态压缩+bfs)
用(x,s)表示一个状态,x表示机器人的位置,s表示其他位置有没有物体.用个fa数组和act数组记录和打印路径,转移的时候判断一下是不是机器人在动. #include<bits/stdc++.h ...
- 复杂UI的组织-创建者模式-uitableview思想
复杂节目的组织-创建者模式-uitableview思想 整体说明,部件规格说明
- BFS 简单思想以及代码
BFS(广搜思想) 广度优先搜索 广度优先搜索是图论的搜索算法之一,以下便进行简单叙述 对于每一个顶点来说,都存在着三种颜色 白色,灰色,黑色 而对于每个顶点,都有三种数据类型 颜色类型,前驱或者父节 ...
- Python静态方法 类方法
通常情况下,类中函数中定义的所有函数,,都是对象的绑定方法,除此之外,还有专门的静态方法和类方法,这两个是专门给类使用的,但是对象非要调用也是不会报错的. 对象在调用的时候会把自己传递给self,也就 ...
- MFC:Unicode和多字节字符集下 CString和char的转化(MFC中)
2011-05-16 00:10 1166人阅读 评论(0) 收藏 举报 mfcdelete Unicode下 CString转TCHAR TCHAR* szMsg = new TCHAR[wcsle ...
- iframe的document操作
导语: 在我写网页代填插件的时候,有遇到拿不到input元素的时候,这时候我去看元素布局,发现有些网站登录那一块是用iframe标签写的,这时候我需要取到的那就是iframe标签下input元素 1. ...
- java基础—方法重载(overload)
一.方法的重载 方法名一样,但参数不一样,这就是重载(overload). 所谓的参数不一样,主要有两点:第一是参数的个数不一样,第二是参数的类型不一样.只要这两方面有其中的一方面不一样就可以构成方法 ...
- inner join 和 left join 的区别
1.left join.right join.inner join的区别 left join(左联接) 返回包括左表中的所有记录和右表中联结字段相等的记录 right join(右联接) 返回包括右表 ...
- OC和C++的区别
C++语言特点: 1.在C语言的基础上进行扩充和完善,使C++兼容了C语言的面向过程特点,又成为了一种面向对象的程序设计语言: 2.可以使用抽象数据类型进行基于对象的编程: 3.可以使用多继承.多态进 ...
- 函数的扩展——箭头函数this的使用
箭头函数中的this指向的是定义时的this,而不是执行时的的this . 举例: 案例中,我们的obj对象中有一个属性x和一个属性show( )方法,show( )通过this打印出x的值,结果是u ...