[bzoj2506]calc_分块处理
calc bzoj-2506
题目大意:给一个长度为n的非负整数序列A1,A2,…,An。现有m个询问,每次询问给出l,r,p,k,问满足l<=i<=r且Ai mod p = k的值i的个数。
注释:$0\le n,m\le 10^5$,$1\le p\le 10^4$。
想法:
直接把询问离线,就变成了求前$i$个数中$\% p=k$的数的个数。
显然直接做肯定做不了。
我们考虑按照mod数分块。
因为上限是$10^4$
如果$p\le 100$,我们用一个$f[i][j]$表示枚举到当前数除以$i$余$j$的数的个数。更新的复杂度为$O(\sqrt{p})$
反之$p<100$,$g[i]$表示枚举到当前数值为$i$的数的个数。更新的复杂度为$O(1)$
考虑如何查询答案:
如果$p\le 100$,直接输出,复杂度为$O(1)$。
反之,我们需要查询$g[k]+g[k+p]+g[k+2p]+...$。时间复杂度为$O(\sqrt{q})$。
故此总时间复杂度为$O(mlogm+n\sqrt{p})$。
Code:
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define N 100010
using namespace std;
int f1[110][110],f2[N];
int ans[2][N],a[N],cnt;
inline char nc() {static char *p1,*p2,buf[100000]; return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;}
int rd() {int x=0; char c=nc(); while(!isdigit(c)) c=nc(); while(isdigit(c)) x=(x<<3)+(x<<1)+(c^48),c=nc(); return x;}
struct Node {int dic,id,f,p,k;}q[N<<1];
inline bool cmp(const Node &a,const Node &b) {return a.dic<b.dic;}
inline void update(int x)
{
for(int i=1;i<=100;i++) f1[i][x%i]++;
f2[x]++;
}
int main()
{
// freopen("bzoj2506.in","r",stdin);
int mx=0;
int n=rd(),m=rd(); for(int i=1;i<=n;i++) a[i]=rd(),mx=max(mx,a[i]);
for(int i=1;i<=m;i++)
{
int l=rd(),r=rd(),p=rd(),k=rd();
q[++cnt].dic=l-1,q[cnt].p=p,q[cnt].k=k,q[cnt].f=0,q[cnt].id=i;
q[++cnt].dic= r; q[cnt].p=p,q[cnt].k=k,q[cnt].f=1,q[cnt].id=i;
}
sort(q+1,q+cnt+1,cmp);
int pnt=0;
for(int i=1;i<=cnt;i++)
{
while(pnt<=q[i].dic) update(a[pnt]),pnt++;
int p=q[i].p,k=q[i].k;
if(p<=100)
{
ans[q[i].f][q[i].id]=f1[p][k];
}
else for(int j=k;j<=mx;j+=p)
{
ans[q[i].f][q[i].id]+=f2[j];
}
}
for(int i=1;i<=m;i++) printf("%d\n",ans[1][i]-ans[0][i]);
return 0;
}
小结:这种根据当前范围思考不同的解决策略的分类讨论思想是极其重要的。
[bzoj2506]calc_分块处理的更多相关文章
- PHP搭建大文件切割分块上传功能
背景 在网站开发中,文件上传是很常见的一个功能.相信很多人都会遇到这种情况,想传一个文件上去,然后网页提示"该文件过大".因为一般情况下,我们都需要对上传的文件大小做限制,防止出现 ...
- POJ2104 K-th Number [分块做法]
传送:主席树做法http://www.cnblogs.com/candy99/p/6160704.html 做那倒带修改的主席树时就发现分块可以做,然后就试了试 思想和教主的魔法差不多,只不过那个是求 ...
- HDU 4467 分块
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4467 题意:给定n个点m条边的无向图,点被染色(黑0/白1),边带边权.然后q个询问.询问分为两种: ...
- 2016 ACM/ICPC Asia Regional Dalian Online 1010 Weak Pair dfs序+分块
Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 262144/262144 K (Java/Others)Total Submissio ...
- CC countari & 分块+FFT
题意: 求一个序列中顺序的长度为3的等差数列. SOL: 对于这种计数问题都是用个数的卷积来进行统计.然而对于这个题有顺序的限制,不好直接统计,于是竟然可以分块?惊为天人... 考虑分块以后的序列: ...
- bzoj2002弹(dan)飞绵羊 分块水过
据说是道lct求深度的题 但是在小猫大的指点下用分块就n^1.5水过了 = =数据忘记加强系列 代码极其不美观,原因是一开始是听小猫大讲的题意,还以为是弹到最前面... #include <cs ...
- BZOJ2506: calc
Description 给一个长度为n的非负整数序列A1,A2,…,An.现有m个询问,每次询问给出l,r,p,k,问满足l<=i<=r且Ai mod p = k的值 ...
- C语言两种查找方式(分块查找,二分法)
二分法(必须要保证数据是有序排列的): 分块查找(数据有如下特点:块间有序,块内无序):
- [New Portal]Windows Azure Storage (14) 使用Azure Blob的PutBlock方法,实现文件的分块、离线上传
<Windows Azure Platform 系列文章目录> 相关内容 Windows Azure Platform (二十二) Windows Azure Storage Servic ...
随机推荐
- webpack 报错:Module build failed: Unknown word (1:1)
解决方法:一是确保css配置里的"style-loader"必须在"css-loader"之前,二是将整个css配置注释掉,如下图:
- 4 Visual Effects 视觉效果 读书笔记 第四章
4 Visual Effects 视觉效果 读书笔记 第四章 Well, circles and ovals are good, but how about drawing r ...
- PPTP的搭建
一.准备 1.检查是否支持pptp modprobe ppp-compress-18 && echo yes yes支持 2.是否开启tun cat /dev/net/tun 返回ca ...
- 【ARM开发板】迅为IMX6开发板QT下LVDS和HDMI双屏异显
本文转自迅为论坛:http://www.topeetboard.com 平台:迅为-IMX6开发板 首先开发板分别连接9.7寸屏和HDMI显示器,然后使用MfgTool工具烧写QT系统,然后拨码开关设 ...
- Oracle中的执行计划
使用autotrace sqlplus系统参数:SQL> set autotrace trace onSQL> select * from dual;DUM---XExecution Pl ...
- 解决docker pull镜像速度慢的问题
直接下载Docker镜像时,由于种种原因,经常下载失败,即使连接成功也是很慢,怎么办呢 目前我知道可以提升速度的办法:DaoCloud 提供Docker Hub Mirror服务 用户可以进入Dock ...
- 使用Caliburn.Micro系列2:Convention
CM中实现一个比较有意思的特性,就是智能匹配. 通常使用MVVM的写法:在匹配 View和ViewModel时会使用DataContext,在匹配数据属性时使用Binding,在匹配事件命令时使用Co ...
- Java 类执行顺序
1.如果父类有静态成员赋值或者静态初始化块,执行静态成员赋值和静态初始化块2.如果类有静态成员赋值或者静态初始化块,执行静态成员赋值和静态初始化块3.将类的成员赋予初值(原始类型的成员的值为规定值,例 ...
- Vue之x-template(1)
今天,我们来讲一个比较有趣的一个功能吧 先来看一段代码示例: <html> <head> <meta charset="utf-8"> < ...
- Debug:This kind of launch is configured to openthe debug perspective when it解决办法
http://blog.sina.com.cn/s/blog_7ca3aa020100zlha.html 启动tomcat时,myeclipse报错: This kind of launch is c ...