洛谷P4887 第十四分块(前体)(二次离线莫队)
题面
题解
lxl大毒瘤
我们考虑莫队,在移动端点的时候相当于我们需要快速计算一个区间内和当前数字异或和中\(1\)的个数为\(k\)的数有几个,而这个显然是可以差分的,也就是\([l,r]\)的询问可以拆成\([1,r]-[1,l-1]\)
我们考虑莫队移动指针的过程,以\([l,r]\)移动左指针到\(p\)为例,要减去的答案是\(l\)和\([1,r]-[1,l-1]\),\(l+1\)和\([1,r]-[1,l]\),...,总的来说,我们我们要对于\([1,r]\)这个前缀计算\([l,p-1]\)的答案,并对每一个\(i\)计算出它和\([1,i-1]\)的答案并做个前缀和
对于前面的,因为前缀是固定的,我们可以在\(r\)上面开一个\(vector\),然后对于每一个\(r\),枚举所有的\([l,p-1]\)暴力加上贡献。所以我们现在需要维护一个\(O(n)\)次插入和\(O(m\sqrt{m})\)次询问的数据结构,扫描线的话可以做到\(O({14\choose 7})\)插入和\(O(1)\)查询,复杂度足够了
然后是后面的,用上面那个数据结构就可以了
//minamoto
#include<bits/stdc++.h>
#define R register
#define pb push_back
#define inline __inline__ __attribute__((always_inline))
#define fp(i,a,b) for(R int i=(a),I=(b)+1;i<I;++i)
#define fd(i,a,b) for(R int i=(a),I=(b)-1;i>I;--i)
#define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)
template<class T>inline bool cmax(T&a,const T&b){return a<b?a=b,1:0;}
template<class T>inline bool cmin(T&a,const T&b){return a>b?a=b,1:0;}
using namespace std;
typedef long long ll;
char buf[1<<21],*p1=buf,*p2=buf;
inline char getc(){return p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++;}
int read(){
R int res,f=1;R char ch;
while((ch=getc())>'9'||ch<'0')(ch=='-')&&(f=-1);
for(res=ch-'0';(ch=getc())>='0'&&ch<='9';res=res*10+ch-'0');
return res*f;
}
char sr[1<<21],z[20];int K=-1,Z=0;
inline void Ot(){fwrite(sr,1,K+1,stdout),K=-1;}
void print(R ll x){
if(K>1<<20)Ot();if(x<0)sr[++K]='-',x=-x;
while(z[++Z]=x%10+48,x/=10);
while(sr[++K]=z[Z],--Z);sr[++K]='\n';
}
const int N=2e5+5;
int a[N],bl[N],st[N],blo,top,n,m,k;ll s1[N],s2[N],s[N],ret[N],ans[N],res;
struct node{
int l,r,id;
inline node(){}
inline node(R int li,R int ri,R int ii):l(li),r(ri),id(ii){}
inline bool operator <(const node &b)const{return bl[l]==bl[b.l]?r<b.r:l<b.l;}
}q[N];vector<node>Q[N];
int main(){
// freopen("testdata.in","r",stdin);
n=read(),m=read(),k=read(),blo=250;
for(R int i=0,c,x;i<16384;++i){
c=0;
for(x=i;x;x^=x&-x)++c;
if(c==k)st[++top]=i;
}
fp(i,1,n)a[i]=read(),bl[i]=(i-1)/blo+1;
fp(i,1,m)q[i].l=read(),q[i].r=read(),q[i].id=i;
sort(q+1,q+1+m);
for(R int i=1,l=q[1].r+1,r=q[1].r;i<=m;++i){
if(l<q[i].l)Q[r].pb(node(l,q[i].l-1,q[i].id<<1));
else if(l>q[i].l)Q[r].pb(node(q[i].l,l-1,q[i].id<<1));
l=q[i].l;
if(r<q[i].r)Q[l-1].pb(node(r+1,q[i].r,q[i].id<<1|1));
else if(r>q[i].r)Q[l-1].pb(node(q[i].r+1,r,q[i].id<<1|1));
r=q[i].r;
}
fp(i,1,n){
s1[i]=s1[i-1]+s[a[i]];
fp(k,1,top)++s[a[i]^st[k]];
s2[i]=s2[i-1]+s[a[i]];
for(vector<node>::iterator it=Q[i].begin();it!=Q[i].end();++it)
fp(k,it->l,it->r)ret[it->id]+=s[a[k]];
}
for(R int i=1,l=q[1].r+1,r=q[1].r;i<=m;++i){
if(l<q[i].l)res-=ret[q[i].id<<1]-s2[q[i].l-1]+s2[l-1];
else if(l>q[i].l)res+=ret[q[i].id<<1]-s2[l-1]+s2[q[i].l-1];
l=q[i].l;
if(r<q[i].r)res+=s1[q[i].r]-s1[r]-ret[q[i].id<<1|1];
else if(r>q[i].r)res-=s1[r]-s1[q[i].r]-ret[q[i].id<<1|1];
r=q[i].r,ans[q[i].id]=res;
}
fp(i,1,m)print(ans[i]);
return Ot(),0;
}
洛谷P4887 第十四分块(前体)(二次离线莫队)的更多相关文章
- [洛谷P4887]第十四分块(前体)
题目大意: 给定一个长度为\(n\)的序列\(a\),\(k\),和\(m\)次询问. 每次询问给定区间\([l,r]\),求满足\(l\leqslant i< j\leqslant r\)且\ ...
- P4887 第十四分块(前体) 莫队
题意: 给你一个序列,每次询问l,r问多少个a[i]^a[j]有k个1,k固定. 序列长度1e5,a[i]<=2^14 时限1s,空间40M 题解: 个人其实开始没什么思路,看了题解也好久,题解 ...
- 洛谷P4689 [Ynoi2016]这是我自己的发明(莫队,树的dfn序,map,容斥原理)
洛谷题目传送门 具体思路看别的题解吧.这里只提两个可能对常数和代码长度有优化的处理方法. I 把一个询问拆成\(9\)个甚至\(16\)个莫队询问实在是有点珂怕. 发现询问的一边要么是一个区间,要么是 ...
- 【洛谷5398】[Ynoi2018]GOSICK(二次离线莫队)
题目: 洛谷 5398 当我刚学莫队的时候,他们告诉我莫队能解决几乎所有区间问题: 现在,当我发现一个区间问题似乎难以用我所了解的莫队解决的时候,他们就把这题的正解叫做 XXX 莫队.--题记 (以上 ...
- 洛谷 P4887 -【模板】莫队二次离线(第十四分块(前体))(莫队二次离线)
题面传送门 莫队二次离线 mol ban tea,大概是这道题让我第一次听说有这东西? 首先看到这类数数对的问题可以考虑莫队,记 \(S\) 为二进制下有 \(k\) 个 \(1\) 的数集,我们实时 ...
- Luogu4887 第十四分块(前体)
sto \(lxl\) orz 考虑莫队,每次移动端点,我们都要询问区间内和当前数字异或有 \(k\) 个 \(1\) 的数字个数 询问 \([l,r]\) 可以再次离线,拆成询问 \([1,l-1] ...
- 【LuoguP4887】第十四分块(前体)
题目链接 题意 区间两数异或在二进制下有 \(k\) 个 \(1\) 的对数. Sol 普通莫队的话,如果要实时维护好区间内的答案需要支持区间对一个数求答案. 直接做不是很好做,容易发现其实这也就是一 ...
- 【洛谷1494】[国家集训队] 小Z的袜子(莫队)
点此看题面 大致题意: 有\(N\)只从\(1\sim N\)编号的袜子,告诉你每只袜子的颜色,\(M\)组询问,每组询问给你一个区间\([L\sim R]\),让你求出小Z随机抽出\(2\)只袜子时 ...
- 洛谷P5398 [Ynoi2018]GOSICK(二次离线莫队)
题面 传送门 题解 维包一生推 首先请确保您会二次离线莫队 那么我们现在的问题就是怎么转移了,对于\(i\)和前缀\([1,r]\)的贡献,我们拆成\(b_i\)和\(c_i\)两部分,其中\(b_i ...
随机推荐
- [luoguP1962] 斐波那契数列(矩阵快速幂)
传送门 解析详见julao博客连接 http://worldframe.top/2017/05/10/清单-数学方法-——-矩阵/ ——代码 #include <cstdio> #incl ...
- hdu - 2660 Accepted Necklace (二维费用的背包问题)
http://acm.hdu.edu.cn/showproblem.php?pid=2660 f[v][u]=max(f[v][u],f[v-1][u-w[i]]+v[i]; 注意中间一层必须逆序循环 ...
- [bzoj2561]最小生成树_网络流_最小割_最小生成树
最小生成树 bzoj-2561 题目大意:题目链接. 注释:略. 想法: 我们发现: 如果一条权值为$L$的边想加入到最小生成树上的话,需要满足一下条件. 就是求出原图的最小生成树之后,这个边当做非树 ...
- POJ3728 The merchant解题报告
Description There are N cities in a country, and there is one and only one simple path between each ...
- JSP中HTTP状态码
以下内容引用自http://wiki.jikexueyuan.com/project/jsp/http-status-codes.html: HTTP请求格式和HTTP响应消息的格式一样,都有以下结构 ...
- 七牛云上传图片到cdn,cdn返回的是一个只有图片的url
如:
- 移动硬盘/U盘上装Windows 7旗舰版(VHD版)
真正的移动版WIN7,在移动硬盘/U盘上运行的WIN7 工具准备 - 联想Y450本本,已安装Windows 7旗舰版(或者WINPE3.0版),用来给移动WIN7做引导 -Win7.vhd,15G, ...
- C#.NET 无法直接启动带有类库输出类型的项目怎么办
我把Driver.cs文件去掉了一行注释,发现报错 右击这个解决方案,选择属性,然后再启动项目中改成MySample
- vux 全局使用 loading / toast / alert
1.入口文件 main.js import { LoadingPlugin, ToastPlugin, AlertPlugin } from 'vux' Vue.use(LoadingPlugin); ...
- NYOJ 298-点的变换(经典矩阵解决点平移、缩放、翻转和旋转)
题目地址:NYOJ 298 思路:该题假设用对每一个点模拟的操作.时间复杂度为O(n+m),结果肯定超时.然而利用矩阵乘法能够在O(m)的时间内把全部的操作合并为一个矩阵,然后每一个点与该矩阵相乘能够 ...