简介:分块算法主要是把区间划分成sqrt(n)块,从而降低暴力的复杂度,

其实这算是一种优化的暴力吧,复杂度O(n*sqrt(n))

题意:给定一个数列:a[i]    (1<= i <= n)    K[j]表示 在区间 [l,r]中j出现的次数。

有t个查询,每个查询l,r,对区间内所有a[i],求sigma(K[a[i]]^2*a[i])

思路:离线+分块处理

分块和离线处理:

将n个数分成sqrt(n)块,设每块有bsize个数, 并且我们计算出每个询问的左端点所在的块号(q[i].b = q[i].l/ bsize)。

对所有询问进行排序:

先按块号排序(块号小的在前),如果块号相等就要右端点排序(右端点小的在前)

解法:每次跟上次的询问区间比较,把多出来的减掉,把少的加上去。 当然第一个询问直接算。

如果一个数已经出现了x次,那么需要累加(2*x+1)*a[i],因为(x+1)^2*a[i] = (x^2 +2*x + 1)*a[i],

x^2*a[i]是出现x次的结果,(x+1)^2 * a[i]是出现x+1次的结果。就是暴力的处理。

复杂度分析:

处理左端点的复杂度:

对于相邻询问左端点的差值不会超过sqrt(n), 所以t个询问的总体复杂度为O(t*sqrt(n))。

处理右端点的复杂度:

对于每个块内的几个查询,因为right是单调递增的,所以极限复杂度为O(n),  而且一共有sqrt(n)个块

所以总体复杂度位O(n*sqrt(n));

因此总的时间复杂度为O(t*sqrt(n)  +  n*sqrt(n))。

为什么选择sqrt(n)分块:

我们从上面复杂度分析中可以得知, 左右断点的复杂度是独立的,

当块数少了,左边复杂度加大,右边复杂度减少,

反之 当块数多了,左边复杂度减少,右边复杂度加大,

块数选择sqrt(n)是为了总体复杂度的最小。

代码:

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
const int maxn = 200005;
typedef long long LL;
LL a[maxn], cnt[maxn * 5], ans[maxn], res;
int L, R; struct node {
int l, r, b, id;
bool operator <(const node &t) const {
if (b == t.b)
return r < t.r;
return b < t.b;
}
} q[maxn]; LL query(int x, int y, int flag) {
if (flag) {
for (int i = x; i < L; i++) {
res += ((cnt[a[i]] << 1) + 1) * a[i];
cnt[a[i]]++;
}
for (int i = L; i < x; i++) {
cnt[a[i]]--;
res -= ((cnt[a[i]] << 1) + 1) * a[i];
}
for (int i = y + 1; i <= R; i++) {
cnt[a[i]]--;
res -= ((cnt[a[i]] << 1) + 1) * a[i];
}
for (int i = R + 1; i <= y; i++) {
res += ((cnt[a[i]] << 1) + 1) * a[i];
cnt[a[i]]++;
} } else {
for (int i = x; i <= y; i++) {
res += ((cnt[a[i]] << 1) + 1) * a[i];
cnt[a[i]]++;
}
}
L = x, R = y;
return res;
}
int n, t;
int main() {
int i; scanf("%d%d", &n, &t);
for (i = 1; i <= n; i++)
scanf("%I64d", &a[i]);
int bsize = sqrt(n + 0.5); for (i = 0; i < t; i++) {
scanf("%d%d", &q[i].l, &q[i].r);
q[i].b = q[i].l / bsize;
q[i].id = i;
} sort(q, q + t);
memset(cnt, 0, sizeof(cnt));
res = 0;
for (i = 0; i < t; i++)
ans[q[i].id] = query(q[i].l, q[i].r, i); for (i = 0; i < t; i++)
printf("%I64d\n", ans[i]); return 0;
}

[置顶] CF 86D Powerful array 分块算法入门,n*sqrt(n)的更多相关文章

  1. CF 86D Powerful array 【分块算法,n*sqrt(n)】

    给定一个数列:A1, A2,……,An,定义Ks为区间(l,r)中s出现的次数. t个查询,每个查询l,r,对区间内所有a[i],求sigma(K^2*a[i]) 离线+分块 将n个数分成sqrt(n ...

  2. CF 86D Powerful array

    离线+分块 将n个数分成sqrt(n)块. 对所有询问进行排序,排序标准:       1. Q[i].left /block_size < Q[j].left / block_size (块号 ...

  3. Codeforces 86D Powerful array (莫队算法)

    题目链接 Powerful array 给你n个数,m次询问,Ks为区间内s的数目,求区间[L,R]之间所有Ks*Ks*s的和. $1<=n,m<=200000,   1<=s< ...

  4. CodeForces 86D Powerful array(莫队+优化)

    D. Powerful array time limit per test 5 seconds memory limit per test 256 megabytes input standard i ...

  5. Codeforces 86D Powerful array (莫队)

    D. Powerful array time limit per test 5 seconds memory limit per test 256 megabytes input standard i ...

  6. [置顶] hdu2815 扩展Baby step,Giant step入门

    题意:求满足a^x=b(mod n)的最小的整数x. 分析:很多地方写到n是素数的时候可以用Baby step,Giant step, 其实研究过Baby step,Giant step算法以后,你会 ...

  7. Codeforces#86D Powerful array(分块暴力)

    Description An array of positive integers a1, a2, ..., an is given. Let us consider its arbitrary su ...

  8. Codeforces 86D Powerful array(莫队算法)

    和BZOJ2038差不多..复习一下. #include<cstdio> #include<cmath> #include<algorithm> using nam ...

  9. Codeforces 86D - Powerful array(莫队算法)

    题目链接:http://codeforces.com/problemset/problem/86/D 题目大意:给定一个数组,每次询问一个区间[l,r],设cnt[i]为数字i在该区间内的出现次数,求 ...

随机推荐

  1. C#高级编程零散知识点

    1.206-实现单链表的添加和插入 207-实现单链表的其他功能和 3.209-Lambda表达式 4.301-栈的介绍和BCL中的栈 4.501-进程和线程的概念[00_12_06][2015122 ...

  2. Oracle Split 函数

    为了让 PL/SQL 函数返回数据的多个行,必须通过返回一个 REF CURSOR 或一个数据集合来完成.REF CURSOR 的这种情况局限于可以从查询中选择的数据,而整个集合在可以返回前,必须进行 ...

  3. Datamatrix码

    DataMatrix二维条码原名Datacode,由美国国际资料公司(International Data Matrix, 简称ID Matrix)於1989年发明.DataMatrix二维条码是一种 ...

  4. 开机时进入 grub rescue>的解决方法

    本机是centOS7和win8的双系统 之前在win上把一个空的磁盘空间释放了 可能造成了grub的一些问题 具体还没有研究过 开机后无法正常进入grub引导画面 而是跳出一串英文+ grub res ...

  5. Phalcon框架中的另类使用

    不像传统的PHP框架,假设框架想被还有一个框架使用仅仅能通过rpc或是引入文件等的方式.Phalcon能够在其他框架中直接使用.这是因为Phalcon是以扩展的形式存在的,在server载入时会直接载 ...

  6. 将其它图片格式转为.eps格式

    假设是用origin的话就不存在这个问题. 倘若你是用excel绘图的话. 1.先将excel导出为.pdf文件 2.用Adobe Acrobat Pro 打开pdf,用其它的pdf软件貌似不行 3. ...

  7. Ext JS学习第四天 我们所熟悉的javascript(三)

    此文用来记录学习笔记: •javascript之函数 •this关键字的使用 –this关键字总是指向调用者,谁调用函数,this就指向谁 •call.apply的使用 –call和apply主要应用 ...

  8. for循环语句之求和,阶乘,求偶,求n次篮球蹦起高度

    for循环语句格式: ;;/*循环条件*/i++/*状态改变*/) { //循环体,执行代码:(break;跳出循环体) } for 穷举法用循环把各种可能的情况都走一遍,然后用if条件把满足要求的结 ...

  9. python学习(一)

    1 python一切皆为对象,因为现实 包含了一系列的数据和操作这些数据的方法的一个整体,就叫作对象. 自行车 属性:手刹车,轮胎,脚踏板方法:如何前进的方法,控制停止的方法,控制方向 实际内容 男人 ...

  10. chapter3习题

    // 2013年11月4日21:47:21 # include <stdio.h> # include <math.h> int main() { int n; double ...