hdu 5381 The sum of gcd(线段树+gcd)
将查询离线处理,依照r排序,然后从左向右处理每一个A[i],碰到查询时处理。用线段树维护。每一个节点表示从[l,i]中以l为起始的区间gcd总和。所以每次改动时须要处理[1,i-1]与i的gcd值。可是由于gcd值是递减的,成log级,对于每一个gcd值记录其区间就可以。然后用线段树段改动,可是是改动一个等差数列。
#include <cstdio>
#include <cstring>
#include <vector>
#include <algorithm> using namespace std;
typedef long long ll;
const int maxn = 1e4 + 5;
#define lson(x) (x<<1)
#define rson(x) ((x<<1)|1) namespace SegTree {
int lc[maxn << 2], rc[maxn << 2];
ll A[maxn << 2], D[maxn << 2], S[maxn << 2]; void maintain (int u, ll a, ll d) {
A[u] += a;
D[u] += d;
int n = rc[u] - lc[u] + 1;
S[u] += a * n + d * (n-1) * n / 2;
} void pushup(int u) {
S[u] = S[lson(u)] + S[rson(u)];
} void pushdown (int u) {
if (A[u] || D[u]) {
int mid = ((lc[u] + rc[u]) >> 1) + 1;
maintain(lson(u), A[u], D[u]);
maintain(rson(u), A[u] + (mid - lc[u]) * D[u], D[u]);
A[u] = D[u] = 0;
}
} void build (int u, int l, int r) {
lc[u] = l;
rc[u] = r;
A[u] = D[u] = S[u] = 0; if (l == r) {
return;
} int mid = (l + r) >> 1;
build (lson(u), l, mid);
build (rson(u), mid + 1, r);
pushup(u);
} void modify (int u, int l, int r, ll a, ll d) {
if (l <= lc[u] && rc[u] <= r) {
maintain(u, a + (lc[u] - l) * d, d);
return;
} pushdown(u);
int mid = (lc[u] + rc[u]) >> 1;
if (l <= mid)
modify(lson(u), l, r, a, d);
if (r > mid)
modify(rson(u), l, r, a, d);
pushup(u);
} ll query (int u, int l, int r) {
if (l <= lc[u] && rc[u] <= r)
return S[u]; pushdown(u);
ll ret = 0;
int mid = (lc[u] + rc[u]) >> 1;
if (l <= mid)
ret += query(lson(u), l, r);
if (r > mid)
ret += query(rson(u), l, r);
return ret;
}
}; int gcd (int a, int b) {
return b == 0 ? a : gcd(b, a%b);
} struct State {
int l, r, idx;
State (int l = 0, int r = 0, int idx = 0): l(l), r(r), idx(idx) {}
bool operator < (const State& u) const { return r < u.r; }
}S[maxn], G[maxn];
int N, Q, A[maxn];
ll R[maxn]; void init () {
scanf("%d", &N);
for (int i = 1; i <= N; i++)
scanf("%d", &A[i]);
scanf("%d", &Q);
for (int i = 1; i <= Q; i++) {
scanf("%d%d", &S[i].l, &S[i].r);
S[i].idx = i;
}
sort(S + 1, S + 1 + Q); SegTree::build(1, 1, N);
} void solve () {
int n = 0, p = 1;
for (int i = 1; i <= N; i++) {
for (int j = 0; j < n; j++)
G[j].idx = gcd(G[j].idx, A[i]);
G[n++] = State(i, i, A[i]);
int mv = 0;
for (int j = 1; j < n; j++) {
if (G[mv].idx == G[j].idx)
G[mv] = State(G[mv].l, G[j].r, G[j].idx);
else
G[++mv] = G[j];
}
n = mv + 1;
//printf("%d:%d\n", i, n);
for (int i = 0; i < n; i++) {
// printf("%d %d %d\n", G[i].l, G[i].r, G[i].idx);
int k = G[i].r - G[i].l + 1;
SegTree::modify(1, G[i].l, G[i].r, 1LL * G[i].idx * k, -G[i].idx);
if (G[i].l > 1)
SegTree::modify(1, 1, G[i].l - 1, 1LL * G[i].idx * k, 0);
} while (p <= N && S[p].r == i) {
R[S[p].idx] = SegTree::query(1, S[p].l, S[p].l);
p++;
}
}
} int main () {
int cas;
scanf("%d", &cas);
while (cas--) {
init ();
solve ();
for (int i = 1; i <= Q; i++)
printf("%lld\n", R[i]);
}
return 0;
}
hdu 5381 The sum of gcd(线段树+gcd)的更多相关文章
- HDU 1394 Minimum Inversion Number(线段树求最小逆序数对)
HDU 1394 Minimum Inversion Number(线段树求最小逆序数对) ACM 题目地址:HDU 1394 Minimum Inversion Number 题意: 给一个序列由 ...
- codeforces 1217E E. Sum Queries? (线段树
codeforces 1217E E. Sum Queries? (线段树 传送门:https://codeforces.com/contest/1217/problem/E 题意: n个数,m次询问 ...
- hdu 5274 Dylans loves tree(LCA + 线段树)
Dylans loves tree Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 131072/131072 K (Java/Othe ...
- HDU 3074.Multiply game-区间乘法-线段树(单点更新、区间查询),上推标记取模
Multiply game Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Tot ...
- Bash and a Tough Math Puzzle CodeForces 914D 线段树+gcd数论
Bash and a Tough Math Puzzle CodeForces 914D 线段树+gcd数论 题意 给你一段数,然后小明去猜某一区间内的gcd,这里不一定是准确值,如果在这个区间内改变 ...
- 2015 Multi-University Training Contest 8 hdu 5381 The sum of gcd
The sum of gcd Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)To ...
- hdu 5381 The sum of gcd
知道对于一个数列,如果以x为左(右)端点,往右走,则最多会有log(a[x])个不同的gcd,并且有递减性 所以会分成log段,每一段的gcd相同 那我们可以预处理出对于每一个位置,以这个位置为左端点 ...
- HDU 5381 The sum of gcd (技巧,莫队算法)
题意:有一个含n个元素的序列,接下来有q个询问区间,对每个询问区间输出其 f(L,R) 值. 思路: 天真单纯地以为是道超级水题,不管多少个询问,计算量顶多就是O(n2) ,就是暴力穷举每个区间,再直 ...
- HDU 5875 Function (2016年大连网络赛 H 线段树+gcd)
很简单的一个题的,结果后台数据有误,自己又太傻卡了3个小时... 题意:给你一串数a再给你一些区间(lef,rig),求出a[lef]%a[lef+1]...%a[rig] 题解:我们可以发现数字a对 ...
随机推荐
- mysql数据库增、删、改、查等基本命令
测试环境:windows7 64位 mysql.exe.Navicat Lite for MySQL.mysql 5.0.18 mysql数据库的基本结构: 数据库(database)包含多个表(ta ...
- 字符串匹配算法(在字符串T中查找是否有与字符串P相同的子串)
T称为目标串(Target)或主串 ,P称为模式串(Pattren) 或子串 1. 简单字符串模式匹配算法 原理:用字符串P的字符依次与字符串T中的字符进行比较,首先将字符串P从第0个位置起与主串T的 ...
- Struts2,get/set 自动获取/设置数据ActionSupport 类
主页:http://struts.apache.org/在用户请求和模块化处理方面以及页面的展现这块,Struts2 发挥了强大的作用:相对于传统的Jsp+Servlet 模式,Struts2 更适合 ...
- Python 字典值相加
#字典值相加 def union_dic(*objs): _keys = set(sum([obj.keys() for obj in objs],[])) _total = {} for _key ...
- 【bzoj1408】[Noi2002]Robot 数论+dp
题目描述 输入 输出 样例输入 3 2 1 3 2 5 1 样例输出 8 6 75 题解 语文题+数论+dp 花了大段讲述什么叫mu,什么叫phi,只是新定义的mu将2看作有平方因子,新定义的phi( ...
- CSA Round 84 Mahattan Center
题目 题目大意 给定平面上的 $n$ 个点和常数 $k$,求 $x$ 轴上的点 $p$ 到 $n$ 个点中距其最近的 $k$ 个点的距离之和的最小值.两点之间的距离定义为曼哈顿距离. 数据范围 $1\ ...
- javasript 按值传递
现象总结如下: 1.JS的基本类型,是按值传递的.2.对于对象而言:分两种情况(a).如果传递给函数的参数是对象,并且修改了这个对象的属性(某些字段的值),那么奇妙的问题就来了.原参数就被修改了.(b ...
- 急速安装Ubuntu/windows双操作系统
本文出自:http://www.cnblogs.com/svitter FAQ 因为很多人都不看FAQ,比如像我,所以直接把FAQ写在最前面,然后把正文卸载最后面逼你看- - 常用软件下载(官网) d ...
- pat 甲级 1078. Hashing (25)
1078. Hashing (25) 时间限制 100 ms 内存限制 65536 kB 代码长度限制 16000 B 判题程序 Standard 作者 CHEN, Yue The task of t ...
- js回到顶部
<!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...