前缀gcd的变化次数是log的,考虑对每一种gcd查询,问题变为查询一段区间是否存在异或前缀和=x/gcd。

  无修改的话显然可以可持久化trie,但这玩意实在没法支持修改。于是考虑分块。

  对于每一块将其中所有块内异或前缀和排序。查询时先看这块与上一块相比gcd有没有变化,如果有对其中每个位置暴力查询,否则在排序后的数组中二分。修改时暴力改每一块的前缀gcd及异或和,被修改的块暴力重构排序数组即可。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
#define N 100010
#define ll long long
ll read()
{
ll x=,f=;char c=getchar();
while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
return x*f;
}
int n,m,a[N];
int block_size,block_tot,L[N],R[N],pos[N];
int gcd_pre[N],xor_pre[N],gcd_block[N],xor_block[N];
struct data
{
int x,i;
bool operator <(const data&a) const
{
return x<a.x||x==a.x&&i<a.i;
}
}v[N];
int gcd(int n,int m){return m==?n:gcd(m,n%m);}
int main()
{
#ifndef ONLINE_JUDGE
freopen("bzoj4028.in","r",stdin);
freopen("bzoj4028.out","w",stdout);
const char LL[]="%I64d\n";
#else
const char LL[]="%lld\n";
#endif
n=read();
for (int i=;i<=n;i++) a[i]=read();
block_size=sqrt(n);block_tot=(n-)/block_size+;
for (int i=;i<=block_tot;i++)
L[i]=(i-)*block_size+,R[i]=min(n,i*block_size);
for (int i=;i<=block_tot;i++)
{
for (int j=L[i];j<=R[i];j++)
gcd_block[i]=gcd(gcd_block[i],a[j]),
xor_block[i]=xor_block[i]^a[j],
pos[j]=i;
gcd_pre[i]=gcd(gcd_pre[i-],gcd_block[i]),xor_pre[i]=xor_pre[i-]^xor_block[i];
}
for (int i=;i<=block_tot;i++)
{
v[L[i]].x=a[L[i]],v[L[i]].i=L[i];
for (int j=L[i]+;j<=R[i];j++)
v[j].x=v[j-].x^a[j],v[j].i=j;
sort(v+L[i],v+R[i]+);
}
m=read();
while (m--)
{
char c=getchar();
while (c<'A'||c>'Z') c=getchar();
if (c=='M')
{
int p=read()+,x=read();a[p]=x;p=pos[p];
gcd_block[p]=xor_block[p]=;
for (int i=L[p];i<=R[p];i++)
gcd_block[p]=gcd(gcd_block[p],a[i]),
xor_block[p]=xor_block[p]^a[i];
for (int i=p;i<=block_tot;i++)
gcd_pre[i]=gcd(gcd_pre[i-],gcd_block[i]),
xor_pre[i]=xor_pre[i-]^xor_block[i];
v[L[p]].x=a[L[p]],v[L[p]].i=L[p];
for (int i=L[p]+;i<=R[p];i++)
v[i].x=v[i-].x^a[i],v[i].i=i;
sort(v+L[p],v+R[p]+);
}
else
{
ll x=read();int ans=;
for (int i=;i<=block_tot;i++)
{
if (ans) break;
if (gcd_pre[i]==gcd_pre[i-])
{
if (x%gcd_pre[i-]==)
{
int l=L[i],r=R[i];
while (l<=r)
{
int mid=l+r>>;
if (v[mid].x>=(x/gcd_pre[i-]^xor_pre[i-])) ans=mid,r=mid-;
else l=mid+;
}
if (ans&&v[ans].x==(x/gcd_pre[i-]^xor_pre[i-])) {ans=v[ans].i;break;}
else ans=;
}
}
else
{
int p=gcd_pre[i-],q=xor_pre[i-];
for (int j=L[i];j<=R[i];j++)
{
p=gcd(p,a[j]),q^=a[j];
if (x%p==&&x/p==q) {ans=j;break;}
}
}
}
if (ans) printf("%d\n",ans-);
else printf("no\n");
}
}
return ;
}

BZOJ4028 HEOI2015公约数数列(分块)的更多相关文章

  1. [BZOJ4028][HEOI2015]公约数数列(分块)

    先发掘性质: 1.xor和gcd均满足交换律与结合率. 2.前缀gcd最多只有O(log)个. 但并没有什么数据结构能同时利用这两个性质,结合Q=10000,考虑分块. 对每块记录这几个信息: 1.块 ...

  2. 【BZOJ4028】[HEOI2015]公约数数列 分块

    [BZOJ4028][HEOI2015]公约数数列 Description 设计一个数据结构. 给定一个正整数数列 a_0, a_1, ..., a_{n - 1},你需要支持以下两种操作: 1. M ...

  3. BZOJ 4028: [HEOI2015]公约数数列 分块

    4028: [HEOI2015]公约数数列 题目连接: http://www.lydsy.com/JudgeOnline/problem.php?id=4028 Description 设计一个数据结 ...

  4. bzoj4028: [HEOI2015]公约数数列

    Description 设计一个数据结构. 给定一个正整数数列 a_0, a_1, ..., a_{n - 1},你需要支持以下两种操作: 1. MODIFY id x: 将 a_{id} 修改为 x ...

  5. [BZOJ4028][HAOI2015]公约数数列[分块+分析暴力]

    题意 题目链接 分析 首先明确 \(xor\) 运算和 \(\rm gcd\) 没有联系! 注意到一个数字取 \(\rm gcd\) 且保证每次取 \(\rm gcd\) 值都会变小的话,最多取 \( ...

  6. 【BZOJ4028】[HEOI2015]公约数数列(分块)

    [BZOJ4028][HEOI2015]公约数数列(分块) 题面 BZOJ 洛谷 题解 看一道题目就不会做系列 首先\(gcd\)最多只会有\(log\)种取值,所以我们可以暴力枚举出所有可能的\(g ...

  7. BZOJ 4028: [HEOI2015]公约数数列 【分块 + 前缀GCD】

    任意门:https://www.lydsy.com/JudgeOnline/problem.php?id=4028 4028: [HEOI2015]公约数数列 Time Limit: 10 Sec   ...

  8. 洛谷 P4108 / loj 2119 [HEOI2015] 公约数数列 题解【分块】

    看样子分块题应该做的还不够. 题目描述 设计一个数据结构. 给定一个正整数数列 \(a_0, a_1, \ldots , a_{n-1}\),你需要支持以下两种操作: MODIFY id x: 将 \ ...

  9. luogu P4108 [HEOI2015]公约数数列——solution

    -by luogu 不会啊.... 然后%了一发题解, 关键是 考虑序列{$a_n$}的前缀gcd序列, 它是单调不升的,且最多只会改变$log_2N$次,因为每变一次至少除2 于是,当我们询问x时: ...

随机推荐

  1. 《C++ Primer》第II部分:C++标准库

    <C++ Primer>第II部分:C++标准库 前言 把<C++ Primer>读薄系列笔记.本篇为第II部分C++标准库,包含全书第8-12章重难点: IO库 顺序容器 范 ...

  2. ubuntu网卡

    查看网卡类型  http://blog.csdn.net/eddy_liu/article/details/6578819 qii@ubuntu:~$ lspci | grep -i net 03:0 ...

  3. underscore.js 分析 第一天

    Underscore 是一个非常实用的Javascript类库. 通过研究他能提高自身的JS水平. 我们看到整个代码被 (function() { /*  代码 */ }.call(this)); 包 ...

  4. Spark性能优化--开发调优与资源调优

    参考: https://tech.meituan.com/spark-tuning-basic.html https://zhuanlan.zhihu.com/p/22024169 一.开发调优 1. ...

  5. Elastic stack ——X-Pack安装

    X-Pack是一个Elastic Stack的扩展,将安全,警报,监视,报告和图形功能包含在一个易于安装的软件包中.在Elasticsearch 5.0.0之前,您必须安装单独的Shield,Watc ...

  6. 微信小程序—day03

    昨日问题 接着上一篇,昨天遇到的scroll-view组件不能滚动的问题. 今天经过调试,发现是由于:图片的实际宽高,大于给image设定的宽高导致的. 解决办法: 减小图片的实际宽高,使之小于ima ...

  7. python基础之全局局部变量及函数参数

    1.局部变量和全局变量 1.1局部变量 局部变量是在函数内部定义的变量,只能在定义函数的内部使用 函数执行结束后,函数内部的局部变量会被系统收回 不同函数可以定义相同名字的局部变量,但是各用个的互不影 ...

  8. 第5章 Linux网络编程基础

    第5章 Linux网络编程基础 5.1 socket地址与API 一.理解字节序 主机字节序一般为小端字节序.网络字节序一般为大端字节序.当格式化的数据在两台使用了不同字节序的主机之间直接传递时,接收 ...

  9. 【WXS全局对象】Math

    Math对象用于执行数学任务. 属性: 名称 说明 Math.E 代表算术常量 e,即自然对数的底数,其值近似于 2.71828. Math.LN10 就是 loge10,即 10 的自然对数,其值近 ...

  10. js设计模式:工厂模式、构造函数模式、原型模式、混合模式

    一.js面向对象程序 var o1 = new Object();     o1.name = "宾宾";     o1.sex = "男";     o1.a ...