首页 > 代码库 > hdu 4630 树状数组+离线操作+GCD
hdu 4630 树状数组+离线操作+GCD
http://acm.hdu.edu.cn/showproblem.php?pid=4630
重新认识了树状数组。
首先要记住那个树形的图,然后+或-lowbit(i)是自己根据具体问题设定的,不要死于+或者-,
树状数组的特点:
1、+lowbit(i)可以到达包含结点i的上一层父节点 所以用于值的更改
2、-lowbit(i)可以到达不包含i所代表区间的上一层父节点 所以用于值的求和---每个不相交的段加起来
3、C[i]的含义也是根据具体问题去做设定的,但是c[i]覆盖了a[i-2^lowbit(i)+1]...a[i]这个长为lowbit(i)的区间
然后谈本题:
pre[i]见注释
#include <cstdio> #include <cstring> #include <algorithm> #include <string> #include <iostream> #include <iomanip> #include <cmath> #include <map> #include <set> #include <queue> using namespace std; #define ls(rt) rt*2 #define rs(rt) rt*2+1 #define ll long long #define ull unsigned long long #define rep(i,s,e) for(int i=s;i<e;i++) #define repe(i,s,e) for(int i=s;i<=e;i++) #define CL(a,b) memset(a,b,sizeof(a)) #define IN(s) freopen(s,"r",stdin) #define OUT(s) freopen(s,"w",stdout) const int MAXN = 50000+100; struct Query{ int l,r; int id; bool operator < (const Query &t)const{return r<t.r;} }q[MAXN]; int N,c[MAXN*2],pre[MAXN],num[MAXN],ans[MAXN];// inline int lowbit(int x){return x&(-x);} void update(int i, int v) { while(i>0) { c[i]=max(c[i],v); i-=lowbit(i); } } int query(int x) { int ret=0; while(x<=N) { ret=max(c[x],ret); x+=lowbit(x); } return ret; } vector<int>divs[MAXN]; void caldivs() { for(int i=1;i<MAXN;i++) for(int j=i;j<MAXN;j+=i) divs[j].push_back(i); } void init() { CL(c,0); CL(pre,0); } int main() { //IN("hdu4630.txt"); int ncase; caldivs(); scanf("%d",&ncase); while(ncase--) { init(); scanf("%d",&N); for(int i=1;i<=N;i++) scanf("%d",&num[i]); int Q; scanf("%d",&Q); for(int i=1;i<=Q;i++) { scanf("%d%d",&q[i].l,&q[i].r); q[i].id=i; } sort(q+1,q+1+Q); for(int i=1,k=1;i<=N;i++) { int x=num[i]; for(int j=0;j<divs[x].size();j++) { int y=divs[x][j]; if(pre[y])//说明现在y是第二次出现 前一次出现位置是pre[y],此次是i update(pre[y],y);//维护c[i]为pre[y]到i的gcd pre[y]=i; } while(q[k].r==i && k<=Q) { ans[q[k].id]=query(q[k].l); k++; } } for(int i=1;i<=Q;i++) { printf("%d\n",ans[i]); } } return 0; }
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。