这题说的给了100000个数有100000次操作 询问 L和R 区间内 在D位上为P的个数,用树状数组存 要开[10][10][100000]的int 开不了但是能开 这么大的unsign short 这样我们将这个树状数组一分为二 50000 个位前面 50000 为后面 我们知道unshort 范围在60000多显然这样是可以开的了

#include <iostream>
#include <cstdio>
#include <string.h>
using namespace std;
typedef __int64 ll;
const int maxn=;
unsigned short C[][][][];
int n;
int lowbit(int x){
return x&(-x);
}
void add(int op, int id, int floor, int x,int v){
int ma=op==?min(,n):n-;
while(x<=ma){
C[op][id][floor][x]+=v;
x+=lowbit(x);
}
}
int sum( int op, int id, int floor, int x){
int ans=;
while(x>){
ans+=C[op][id][floor][x];
x-=lowbit(x);
}
return ans;
}
int A[maxn];
int solve(int L, int R, int D, int P){
int s=;
if(R>){
s=sum(,P,D,R-)+sum(,P,D,);
}else s= sum(,P,D,R);
if(L>){
s=s-sum(,P,D,L--)-sum(,P,D,);
}else s-=sum(,P,D,L-);
return s;
}
int main()
{
int cas,m;
scanf("%d",&cas);
while(cas--){
scanf("%d%d",&n,&m);
memset(C,,sizeof(C));
for(int i=; i<=n; ++i)
{
scanf("%d",&A[i]);
int t =A[i];
int loc=;
int ko=i>?i-:i;
int op=i>?:;
while(loc<){
add(op,t%,loc,ko,);
loc++;
t/=;
}
}
char st[];
int L,R,D,P;
for(int i=; i<m; ++i){
scanf("%s",st);
if(st[]=='Q'){
scanf("%d%d%d%d",&L,&R,&D,&P);
int ans=solve(L,R,D-,P);
printf("%d\n",ans);
}else{
int x,y;
scanf("%d%d",&x,&y);
int t=A[x];
int loc=;
int lc=x;
int op;
if(x>) lc-=,op=;
else op=;
while(loc<){
add(op,t%,loc,lc,-);
loc++;
t/=;
}
A[x]=y;
t = y;
loc=;
while(loc<){
add(op,t%,loc,lc,);
loc++;
t/=;
}
}
}
}
return ;
}

还有一种操作 就是离线,将每个操作记下来, 10000次操作 做10次 先做在第一位 开一个【100000】【10】的数组,然后先存每个数的第一位,然后按照操作顺序回答所有有关第一位的询问,然后同样的操作第二位 第三位... 意思做了10次这样的整套查询操作

#include <iostream>
#include <cstdio>
#include <string.h>
#include <vector>
using namespace std;
typedef __int64 ll;
const int maxn=;
int C[][maxn];
int n,m;
int eps[];
int chose(int V, int t){
return V / eps[t] % ;
}
int lowbit(int x){
return x&(-x);
}
void add( int x, int id , int v){
while(x<=n){
C[id][x]+=v;
x+=lowbit(x);
}
}
int sum( int x, int id ){
int ans=;
while(x>){
ans+=C[id][x];
x-=lowbit(x);
}
return ans;
}
int A[maxn];
int B[maxn];
int ans[maxn];
int op[maxn],L[maxn],R[maxn],D[maxn],P[maxn];
bool use[];
void inti(int lo){
memset( C , , sizeof( C ) );
for(int i=; i<=n; ++i){
B[i]=A[ i ];
int id=chose( B[ i ] , lo );
add(i,id,);
}
}
void solve(int los){
for(int i =; i<m; ++i){
if(op[i]==){
int x=L[i];
add(x,chose(B[x],los),-);
B[x]=R[ i ];
add( x , chose(B[x],los) , );
}else if( op[i]== && D[i]==los ){
ans[i]=sum(R[i],P[i])-sum(L[i]-,P[i]);
}
}
}
int main()
{
eps[]=;
for(int i=; i<=; ++i)
eps[i]=eps[i-]*;
int cas;
scanf("%d",&cas);
while(cas--){
scanf("%d%d",&n,&m);
for(int i=; i<=n; ++ i){
scanf("%d",&A[i]);
A[i];
}
memset(use,false,sizeof(use));
char ch[];
for(int i=; i<m; ++i ){
scanf("%s",ch);
if(ch[]=='S'){
op[i]=;
scanf("%d%d",&L[i],&R[i]);
}else{
op[i]=;
scanf("%d%d%d%d",&L[i],&R[i],&D[i],&P[i]);
use[ D[i] ] = true ;
}
}
for(int i=; i<=; ++i){
if(use[i]==false)continue;
inti(i);
solve(i);
}
for(int i=; i<m; ++i){
if(op[i]==) printf("%d\n",ans[i]);
}
}
return ;
}

hdu5057 分块处理,当数值大于数据范围时树状数组 真是巧 将大数据分为小数据来处理的更多相关文章

  1. 【BZOJ2141】排队 树状数组+分块

    [BZOJ2141]排队 Description 排排坐,吃果果,生果甜嗦嗦,大家笑呵呵.你一个,我一个,大的分给你,小的留给我,吃完果果唱支歌,大家乐和和.红星幼儿园的小朋友们排起了长长地队伍,准备 ...

  2. 【BZOJ 3529】 [Sdoi2014]数表 (莫比乌斯+分块+离线+树状数组)

    3529: [Sdoi2014]数表 Description 有一张N×m的数表,其第i行第j列(1 < =i < =礼,1 < =j < =m)的数值为能同时整除i和j的所有 ...

  3. 【BZOJ 3295】动态逆序对 - 分块+树状数组

    题目描述 给定一个1~n的序列,然后m次删除元素,每次删除之前询问逆序对的个数. 分析:分块+树状数组 (PS:本题的CDQ分治解法见下一篇) 首先将序列分成T块,每一块开一个树状数组,并且先把最初的 ...

  4. 沈阳网络赛J-Ka Chang【分块】【树状数组】【dfs序】

    Given a rooted tree ( the root is node 11 ) of NN nodes. Initially, each node has zero point. Then, ...

  5. Bzoj 3295: [Cqoi2011]动态逆序对 分块,树状数组,逆序对

    3295: [Cqoi2011]动态逆序对 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 2886  Solved: 924[Submit][Stat ...

  6. Bzoj 2141: 排队 分块,逆序对,树状数组

    2141: 排队 Time Limit: 4 Sec  Memory Limit: 259 MBSubmit: 1310  Solved: 517[Submit][Status][Discuss] D ...

  7. HDU 4638 Group 【树状数组,分块乱搞(莫队算法?)】

    根据题目意思,很容易得出,一个区间里面连续的段数即为最少的group数. 题解上面给的是用树状数组维护的. 询问一个区间的时候,可以一个一个的向里面添加,只需要判断a[i]-1 和 a[i]+1是否已 ...

  8. bzoj 4765 普通计算姬(树状数组 + 分块)

    http://www.lydsy.com/JudgeOnline/problem.php?id=4765 很nice的一道题啊(可能是因为卡了n久终于做出来了 题意就是给你一棵带点权的有根树,sum( ...

  9. 牛客小白月赛12 F 华华开始学信息学 (分块+树状数组)

    链接:https://ac.nowcoder.com/acm/contest/392/F来源:牛客网 时间限制:C/C++ 2秒,其他语言4秒 空间限制:C/C++ 32768K,其他语言65536K ...

随机推荐

  1. GIS-009-Cesium 使用

    //加载ArcGIS 发布的地图服务MapServervar url='http://Jason:6080/arcgis/rest/services/SampleWorldCities/MapServ ...

  2. List自定义排序

    List自定义排序我习惯根据Collections.sort重载方法来实现,下面我只实现自己习惯方式.还有一种就是实现Comparable接口. 挺简单的,直接上代码吧. package com.so ...

  3. Python 处理输入输出

    sys.stdin.read() 用于接收标准输入,也就是让用户通过键盘进行输入sys.stdout.write() 用于打印标准输出,也就是把输入的数据输出到屏幕sys.stderr.write() ...

  4. securecrt 的安装

    http://bbs.feng.com/read-htm-tid-6939481.html ssh  -t  ip地址@用户名  -p 22

  5. PHP Web 木马扫描器代码

    scanner.php:<?php/**************PHP Web木马扫描器************************//* [+] 作者: alibaba *//* [+] ...

  6. Chisel常用命令总结

    Chisel简介 Chisel是Facebook开源的一款lldb调试工具,其实就是对系统lldb命令的封装,开发者可以通过简化的命令更方便的进行调试工作.开源地址:https://github.co ...

  7. php第二例

    参考: http://www.php.cn/code/3645.html 前言 由于navicat在linux平台不能很好的支持, PHP的学习转到windows平台. php IDE: PhpSto ...

  8. JS 构造图片Image对象

    var image=new Image(); image.src=""; console.log(image.width);

  9. 微信accesstoken回调

    errcode=-1的时候,开发文档中说明是系统异常,至于具体原因不明 不过有一种原因是AppID以及AppSecret错误 其它可能原因还待发现

  10. TCP协议的基本规则和在Java中的使用

    TCP协议是面向连接的,相对于UDP协议来说效率较低,但是比较安全,数据不容易丢失.TCP协议类似打电话的过程,在一端拨号时必须等待对方回应,确定两端建立了连接通道才能传送信息. 在Java中TCP被 ...