HDU4630(树状数组,离线)
2015-10-28 18:14
204 查看
题意是给你一个1-n的排列,给q个询问,求每个询问区间里面的最大gcd。
先预处理每一个数的因数。然后按照询问的右区间排序,记录到当前询问为止的最左边界,还有某个因子的最右,次右位置,每次这个因子的位置能更新时都用树状数组维护下。
先预处理每一个数的因数。然后按照询问的右区间排序,记录到当前询问为止的最左边界,还有某个因子的最右,次右位置,每次这个因子的位置能更新时都用树状数组维护下。
#include <bits/stdc++.h> using namespace std; #define maxn 51111 vector <int> fac[maxn]; //每个数的因子 int num[maxn], c[maxn], a[maxn], n; int pre[maxn][2]; //有i因子的最大,次大位置 struct node { int l, r, ans, id; bool operator < (const node &a) const { return r < a.r || (r == a.r && l > a.l); } }qu[maxn]; bool cmp (const node &a, const node &b) { return a.id < b.id; } int lowbit (int x) { return x&(-x); } void add (int pos, int num) { for (int i = pos; i > 0; i -= lowbit (i)) { c[i] = max (c[i], num); } } int query (int pos) { int ans = 0; for (int i = pos; i <= n; i+= lowbit (i)) { ans = max (ans , c[i]); } return ans; } void solve (int pos) { int kk = a[pos]; for (int i = 0; i < fac[kk].size (); i++) { int v = fac[kk][i]; //++num[v]; if (pre[v][0] < pos) { pre[v][1] = pre[v][0]; pre[v][0] = pos; } else if (pre[v][1] < pos) { pre[v][1] = pos; } if (pre[v][1] != -1) { add (pre[v][1], v); } } } int main () { //freopen ("in", "r", stdin); ios::sync_with_stdio(false); int t; cin >> t; for (int i = 1; i <= 50000; i++) { fac[i].clear (); } for (int i = 1; i <= 50000; i++) { for (int j = i; j <= 50000; j+= i) fac[j].push_back (i); } while (t--) { cin >> n; for (int i = 1; i <= n; i++) { cin >> a[i]; } int q; cin >> q; for (int i = 1; i <= q; i++) { cin >> qu[i].l >> qu[i].r; qu[i].id = i; } sort (qu+1, qu+1+q); memset (c, 0, sizeof c); memset (num, 0, sizeof num); memset (pre, -1, sizeof pre); int pre_l, pre_r; for (int i = 1; i <= q; i++) { if (i == 1) { pre_r = qu[i].r, pre_l = qu[i].l; for (int pos = qu[i].r; pos >= qu[i].l; pos--) { solve (pos); } } else { for (int j = qu[i].r; j >= pre_r+1; j--) { solve (j); } for (int j = pre_l-1; j >= qu[i].l; j--) { solve (j); } } qu[i].ans = query (qu[i].l); pre_l = min (pre_l, qu[i].l); pre_r = qu[i].r; } sort (qu+1, qu+1+q, cmp); for (int i = 1; i <= q; i++) cout << qu[i].ans << endl; } return 0; }
相关文章推荐
- [转]在Eclipse中使用JUnit4进行单元测试(初级篇)
- 射影几何:基本形有几种?
- Curator Framework的基本使用方法
- git命令总结
- linux uevent分析
- 场景二:synchronized修饰方法
- Could not find file 'Microsoft.Windows.CommonLanguageRuntime, Version=2.0.50727.0' 问题
- andorid serviceconnection
- 封装Android文件基本操作方法
- Python *与** 参数问题
- STM32软件复位(基于库文件V3.5)
- 使用python进行数据转码
- CSS基础学习一:CSS概述
- 【安卓】ExpandableListView中的position
- C++编程 (三)--- 深入C++后台开发
- 【未解决】html中的a标签id与class冲突
- 河北民间管理系统利益相关者分析
- php xml转数组
- Linux CentOS 6.5nginx +tomcat +redis分布式web应用的session共享配置
- python函数参数带星号*