五月月赛 寻宝 exkmp + 主席树
: 寻宝
时间限制: Sec 内存限制: MB
提交: 解决:
[提交] [状态] [讨论版] [命题人:admin]
题目描述
采蘑菇的小西佬找到了一张上古年间的藏宝图,上面画着m座连绵不断的山,他决定去地图上记载的地点探险,可当他到达时,他发现当地其实有n座山,并且由于年代久远有些山也改变了形状,他不能准确的找到地图上藏宝的地点,于是他决定从一些连续的山中随意选择一座山作为地图上的第一个点,当他选定了起点以后,若当前的山与地图上记载山的高度相同他就会继续前往下一座山,否则就停下来,直到他到达地图上的最后一座山位置。我们假设如果他经过了大于k座山他就能找到宝藏(如果第一座山高度就不相同,视为经过0座山),现在的问题就是他有多大的可能性找到宝藏。
题目共有q次询问,每次询问有两个整数l,r代表采蘑菇的小西佬选择的区间,并且会给出第一次询问的k,之后的每个询问的k按照下面的式子算出
k=(ans*10086+666)%x; ans代表上一次询问的答案,x=min(m,r-l+1); 输入
第一行输入三个数n,m,k(<=n,m,k<=1e5); 第二行输入n个数代表n座山的高度(<=a[i]<=1e9) 第二行输入m个数代表m座山的高度(<=b[i]<=1e9) 第四行输入一个数q;(<=q<=1e5) 接下来每行输入两个整数l,r;(<=l<=r<=n) 输出
对于每个询问请输出对应的概率,以逆元的形式输出(mod=1e9+7);
题面
用exkmp预处理出每个点可以向后面跑的长度,由于k是变化的,所以需要主席树来回答。
#include <bits/stdc++.h> using namespace std;
#define pb push_back
#define fi first
#define se second
#define debug(x) cerr<<#x << " := " << x << endl;
#define bug cerr<<"-----------------------"<<endl;
typedef long long ll;
typedef long double ld;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
//const int inf = 0x3f3f3f3f;
const int mod = 1e9+;
/**********showtime************/ const int N = 1e5 + ;
int nxt[N], ex[N];
int S[N], P[N];
int a[N];
int n,m,k;
void GETNEXT() {
int i = , j, po, len=m;
nxt[] = len;
while(P[i] == P[i+] && i+ < len) i++;
nxt[] = i;
po = ;
for(i = ; i < len; i++) {
if(nxt[i-po] + i < nxt[po] + po)
nxt[i] = nxt[i-po];
else {
j=nxt[po] + po - i;
if(j < ) j = ;
while(i + j < len && P[j] == P[j+i])
j++;
nxt[i] = j;
po = i;
}
}
}
void EXKMP()
{
int i = , j, po, len = n, l2 = m;
while(S[i] == P[i] && i < l2 && i < len) i++;
ex[] = i;
po = ;
for(i = ; i < len; i++)
{
if(nxt[i-po] + i < ex[po] + po) ex[i]=nxt[i-po];
else {
j = ex[po] + po - i;
if(j < ) j = ;
while(i + j < len && j < l2 && S[j+i] == P[j]) j++;
ex[i] = j;
po = i;
}
}
} struct node{
int le, ri;
int cnt;
} tree[N * ];
int tot = ;
int root[N];
int build(int le, int ri) {
int p = ++tot;
if(le == ri) {
tree[p].cnt = ;
return p;
}
int mid = (le + ri) >> ;
tree[p].le = build(le, mid);
tree[p].ri = build(mid+, ri);
tree[p].cnt = tree[tree[p].le].cnt + tree[tree[p].ri].cnt;
return p;
}
int ins(int now, int le, int ri, int x, int val) {
int p = ++tot;
tree[p] = tree[now];
if(le == ri) {
tree[p].cnt += val;
return p;
}
int mid = (le + ri) >> ;
if(x <= mid) {
tree[p].le = ins(tree[now].le, le, mid, x, val);
}
else {
tree[p].ri = ins(tree[now].ri, mid+, ri, x, val);
}
tree[p].cnt = tree[tree[p].le].cnt + tree[tree[p].ri].cnt;
return p;
}
int query(int L, int R, int le, int ri, int rt1, int rt2) {
if(L > R) return ;
if(le >= L && ri <= R) {
return tree[rt2].cnt - tree[rt1].cnt;
} int res = ;
int mid = (le + ri) >> ;
if(mid >= L) res += query(L, R, le, mid, tree[rt1].le, tree[rt2].le);
if(mid < R) res += query(L, R, mid+, ri, tree[rt1].ri, tree[rt2].ri);
return res;
}
ll ksm(ll a, ll n) {
ll res = ;
while(n > ) {
if(n & ) res = res * a % mod;
a = a * a % mod;
n = n>>;
}
return res;
}
int main(){
scanf("%d%d%d", &n, &m, &k);
for(int i=; i<n; i++) scanf("%d", &S[i]);
for(int i=; i<m; i++) scanf("%d", &P[i]); GETNEXT();
EXKMP();
root[] = build(, n); for(int i=; i < n; i++){
root[i+] = ins(root[i], , n, ex[i]+, );
} int q; scanf("%d", &q);
int ans = ;
for(int i=; i<=q; i++) { int l, r;
scanf("%d%d", &l, &r);
if(i > ) {
int x = min(m, r - l + );
k = (1ll*ans * + ) % x;
}
int c = query(k+, n, , n, root[l-], root[r]);
// cout<<c << " " << k << endl;
ans = 1ll*c * ksm(r - l + , mod-) % mod;
printf("%d\n", ans);
} return ;
}
五月月赛 寻宝 exkmp + 主席树的更多相关文章
- BZOJ5361[Lydsy1805月赛]对称数——主席树+随机化
题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=5361 好神的一道题啊! 容易看出来是要用维护权值的数据结构,因此树链剖分首先pass掉. ...
- bzoj3207--Hash+主席树
题目大意: 给定一个n个数的序列和m个询问(n,m<=100000)和k,每个询问包含k+2个数字:l,r,b[1],b[2]...b[k],要求输出b[1]~b[k]在[l,r]中是否出现. ...
- bzoj1901--树状数组套主席树
树状数组套主席树模板题... 题目大意: 给定一个含有n个数的序列a[1],a[2],a[3]--a[n],程序必须回答这样的询问:对于给定的i,j,k,在a[i],a[i+1],a[i+2]--a[ ...
- BZOJ 3626: [LNOI2014]LCA [树链剖分 离线|主席树]
3626: [LNOI2014]LCA Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 2050 Solved: 817[Submit][Status ...
- BZOJ 1146: [CTSC2008]网络管理Network [树上带修改主席树]
1146: [CTSC2008]网络管理Network Time Limit: 50 Sec Memory Limit: 162 MBSubmit: 3522 Solved: 1041[Submi ...
- BZOJ 2588: Spoj 10628. Count on a tree [树上主席树]
2588: Spoj 10628. Count on a tree Time Limit: 12 Sec Memory Limit: 128 MBSubmit: 5217 Solved: 1233 ...
- BZOJ 1901: Zju2112 Dynamic Rankings[带修改的主席树]【学习笔记】
1901: Zju2112 Dynamic Rankings Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 7143 Solved: 2968[Su ...
- [bzoj3932][CQOI2015][任务查询系统] (主席树)
Description 最近实验室正在为其管理的超级计算机编制一套任务管理系统,而你被安排完成其中的查询部分.超级计算机中的 任务用三元组(Si,Ei,Pi)描述,(Si,Ei,Pi)表示任务从第Si ...
- [bzoj2588][count on a tree] (主席树+lca)
Description 给定一棵N个节点的树,每个点有一个权值,对于M个询问(u,v,k),你需要回答u xor lastans和v这两个节点间第K小的点权.其中lastans是上一个询问的答案,初始 ...
随机推荐
- @Validated和@Valid的区别?校验级联属性(内部类)
每篇一句 NBA里有两大笑话:一是科比没天赋,二是詹姆斯没技术 相关阅读 [小家Java]深入了解数据校验:Java Bean Validation 2.0(JSR303.JSR349.JSR380) ...
- CEPH 对象存储的系统池介绍
RGW抽象来看就是基于rados集群之上的一个rados-client实例. Object和pool简述 Rados集群网上介绍的文章很多,这里就不一一叙述,主要要说明的是object和pool.在r ...
- SpringMVC学习笔记之---简单入门
SpringMVC简单入门 (一)什么是MVC设计模式 (1)model:模型数据,业务逻辑 (3)view:呈现模型,与用户进行交互 (3)controller:负责接收并处理请求,响应客户端 (二 ...
- 浅谈单例模式及其java实现
单例模式是23种设计模式中比较简单的一种,在此聊一下单例模式. 1.什么是设计模式? 对于没有接触过设计模式的人来说,一听到设计模式这四个字就觉得这个东西很高深莫测,一下子就对这个东西产生了恐惧感,其 ...
- 对vue中nextTick()的理解及使用场景说明
异步更新队列: 首先我们要对vue的数据更新有一定理解: vue是依靠数据驱动视图更新的,该更新的过程是异步的. 即:当侦听到你的数据发生变化时, Vue将开启一个队列(该队列被Vue官方称为异步更新 ...
- react学习(一)--JSX简介
由于在中国银联实习的项目要用到react,所以不得不硬着头皮把react学习一下.这是要往全栈发展吗0.0 正文: 一个最简单的React例子如下, ReactDOM.render( <h1&g ...
- python多线程同步实例分析
进程之间通信与线程同步是一个历久弥新的话题,对编程稍有了解应该都知道,但是细说又说不清.一方面除了工作中可能用的比较少,另一方面就是这些概念牵涉到的东西比较多,而且相对较深.网络编程,服务端编程,并发 ...
- CSS3: @font-face 介绍与使用
@font-face 是CSS3中的一个模块,他主要是把自己定义的Web字体嵌入到你的网页中,随着@font-face模块的出现,我们在Web的开发中使用字体不怕只能使用Web安全字体,你们当中或许有 ...
- cmd命令行带参启动程序
cmd命令行带参启动程序 有一些程序不支持被直接启动,编写代码时,我们可以通过Process类来启动某个进程(某个软件),在不用代码调从而启动某个软件时,windows系统下,通常我们会用到cmd命令 ...
- 用 程序 解决 windows防火墙 的 弹窗 问题
今天用户反馈了一个问题,运行程序弹了个框 这个只有程序第一次运行会出来,之后就不会了. 当然改个程序名字,又会弹出来. 强烈怀疑是写到了注册表,果然被我找到了. “HKEY_LOCAL_MACHINE ...