4028: [HEOI2015]公约数数列
4028: [HEOI2015]公约数数列
Description
设计一个数据结构. 给定一个正整数数列 a_0, a_1, ..., a_{n - 1},你需要支持以下两种操作:
Input
输入数据的第一行包含一个正整数 n.
Output
对于每个 QUERY 询问,在单独的一行中输出结果。如果不存在这样的 p,输出 no.
Sample Input
1353600 5821200 10752000 1670400 3729600 6844320 12544000 117600 59400 640
10
MODIFY 7 20321280
QUERY 162343680
QUERY 1832232960000
MODIFY 0 92160
QUERY 1234567
QUERY 3989856000
QUERY 833018560
MODIFY 3 8600
MODIFY 5 5306112
QUERY 148900352
Sample Output
0
no
2
8
8
HINT
对于 100% 的数据,n <= 100000,q <= 10000,a_i <= 10^9 (0 <= i < n),QUERY x 中的 x <= 10^18,MODIFY id x 中的 0 <= id < n,1 <= x <= 10^9.
代码
#include <cstdio> #include <cmath> #include <algorithm> #define N 100050 #define tp t[n] #define L(i) (i-1)*k+1 #define R(i) min(i*k,n) #define ff() for(int i=L(x);i<=R(x);i++) #define f(a,b,c) for(int a=b;a<=c;a++) using namespace std; typedef long long LL; int n,m,t[N],a[N],k; LL g[N],xo[N],q; struct Mo{LL x;int r;}b[N]; bool cmp(Mo p1,Mo p2) {return p1.x==p2.x?p1.r<p2.r:p1.x<p2.x;} LL gcd(LL x,LL y){return y==0?x:gcd(y,x%y);} void rebuild(int x) {g[x] = xo[x] = 0;ff()g[x]=gcd(g[x],a[i]);ff()xo[x]^=a[i];b[L(x)]=(Mo){a[L(x)],L(x)};f(i,L(x)+1,R(x))b[i]=(Mo){b[i-1].x^a[i],i};sort(b+L(x),b+R(x)+1,cmp); } void modify(){int id,x;scanf("%d%d",&id,&x);a[id+1]=x;rebuild(t[id+1]);} int lb(int l,int r,int p) {while(l<r){int mid=(l+r)>>1;b[mid].x>=p?r=mid:l=mid+1;}return l;} void solve() {scanf("%lld",&q);LL sum=0,sg=0;f(i,1,tp){if(gcd(g[i],sg)==sg){if(q%sg!=0){sum^=xo[i];sg=gcd(sg,g[i]);continue;}LL p =(q/sg)^sum;int pos=lb(L(i),R(i),p);if(b[pos].x==p){printf("%d\n",b[pos].r-1);return ;}} else{LL p=sg,cur=sum;f(j,L(i),R(i)){p=gcd(p,a[j]),cur^=a[j];if(1LL*p*cur==q){printf("%d\n",j-1);return;}}}sum^=xo[i],sg=gcd(sg,g[i]);}puts("no");} int main() {scanf("%d",&n);k=sqrt(n);f(i,1,n)t[i]=(i-1)/k+1,scanf("%d",&a[i]);f(i,1,tp)rebuild(i);scanf("%d",&m);f(i,1,m){char s[10];scanf("%s",s+1);s[1]=='M'?modify():solve();}return 0; }
转载于:https://www.cnblogs.com/muzu/p/7898396.html
4028: [HEOI2015]公约数数列相关推荐
- 【BZOJ4028】[HEOI2015]公约数数列(分块/数量级很小法)
[BZOJ4028][HEOI2015]公约数数列 https://www.luogu.com.cn/problem/P4108 求解最靠前的一个位置前缀gcd和前缀异或的乘积恰好等于x,并且要求支持 ...
- luogu P4108 [HEOI2015]公约数数列(分块、gcd性质)
题解链接 #include <cstdio> #include <algorithm> #include <cstring> #include <vector ...
- [暑假的bzoj刷水记录]
(这篇我就不信有网站来扣) 这个暑假打算刷刷题啥的 但是写博客好累啊 堆一起算了 隔一段更新一下. 7月27号之前刷的的就不写了 , 写的累 代码不贴了,可以找我要啊.. 2017.8.27upd ...
- Noip前的大抱佛脚----赛前任务
赛前任务 tags:任务清单 前言 现在xzy太弱了,而且他最近越来越弱了,天天被爆踩,天天被爆踩 题单不会在作业部落发布,所以可(yi)能(ding)会不及时更新 省选前的练习莫名其妙地成为了Noi ...
- 退役前的做题记录1.0
退役前的做题记录1.0 租酥雨最近很懒qwq,具体表现在写题的时候不想发题解了. 但是想想这样也不太好,就决定发个一句话(半句话到几句话不等)题解上来. 2018-09.18-2018-09.28 [ ...
- 费氏(Fibonacci)数列、最大公约数,最小公倍数
费式数列 说明 Fibonacci为1200年代的欧洲数学家,在他的着作中曾经提到:「若有一只免子每个月生一只小免子,一个月后小免子也开始生产.起初只有一只免子,一个月后就有两只免子,二个月后有三只免 ...
- 【勇敢牛牛,不怕困难】有手就行栏目:头歌教学平台 - 湖南工业大学刘强老师的C语言函数实战课堂作业答案 - > - > {求和+回文数计算+编写函数求表达式的值+阶乘数列+亲密数+公约数公倍数求解}
第一关:求和 任务描述 题目描述:给你一个n,要求你编写一个函数求1+2+-+n. 输入 输入一个n 输出 输出1+2+-+n的和 测试说明 样例输入: 100 样例输出: 5050 分析: 这个是真 ...
- 求水仙花数 以及 最大公约数最小公倍数 以及 冒泡法 以及字母大小写的转换 以及简单选择排序法 以及斐波那契数列
什么是水仙花数呢? 假设一个数ABC是一个三位数. 水仙花数就是 A*A*A +B*B*B +C*C*C = ABC; 水仙花的具体代码. 接着是最大公约数以及最小公倍数 如果两个数是ab,最大公约数 ...
- C语言求数列最小公倍数,C语言复习---获取最小公倍数(公式法:两个数相乘等于最小公倍数乘以最大公约数)...
Android 图片浏览器 从原来位置放大至全屏显示 android 图片浏览器 特点: 1.从网络加载图片,只需要传图片地址数组即可 2.点击图片,从原来位置放大至全屏 3.支持手势操作 4.完全自 ...
最新文章
- Educational Codeforces Round 90 (Rated for Div. 2)部分题解
- elasticsearch实战三部曲之三:搜索操作
- 视频需求超平常数 10 倍,却节省了 60% 的 IT 成本投入是一种什么样的体验?
- Cortex - M3 中断和异常的区别
- Windows 下安装mysql
- 自动设置图片的序号_巧用word软件中的题注功能,让你插入的对象自动排序与更新...
- 五、规则组织的衍生组织——纬山形组织数学模型的建立
- 论文浅尝 | 基于图匹配神经网络的跨语言知识图对齐 (ACL 2019)
- Memcached服务端自动启动(转载)
- 都啥年代了,求你别再说Redis是单线程了!
- scala学习-Description	Resource	Path	Location	Type value toDF is not a member of org.apache.spark.rdd.R
- linux下给源码安装好的php支持pdo_mysql
- Android: app不被系统kill掉
- vivo S10 PRO怎么解锁vivoS10解锁平台刷机教程屏幕锁激活手机锁不记得了可以用这个方法教程重装系统固件软件程序使用方法流程
- easyboot制作xp、win7、ubuntu、老毛桃winpe、深山红叶dos工具箱多合一光盘
- Qt QListView详解
- delphi中panel控件应用
- python中分号中没内容_被“嫌弃”的分号的一生:不要在Python中使用无用分号了...
- 为出海掘金创造更多可能 助力开发者触达全球用户
- Cpp Careless Tony
热门文章
- securecrt鼠标右键的配置
- 51nod 1102 【单调栈】
- Xcode8.0 删除插件路径
- UVA 10635 王子和公主
- F#学习之路(3) 如何组织程序(下)
- python asp.net requests_python3使用requests访问asp页面时出错,返回500
- python 动态_python实现动态创建类的方法分析
- Linux系统诊断实践-内存基础
- 程序员该知道的7个必经阶段
- 卸载vuecli3_卸载vue2.0并升级vue_cli3.0的实例讲解