博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
BZOJ 3744: Gty的妹子序列 [分块]
阅读量:4692 次
发布时间:2019-06-09

本文共 2561 字,大约阅读时间需要 8 分钟。

题意:询问区间内逆序对数


 

感觉这种题都成套路题了

两个预处理$f[i][j]$块i到j的逆序对数,$s[i][j]$前i块$\le j$的有多少个

f我直接处理成到元素j,方便一点

用个树状数组就行了

预处理和查询都带$log$所以还是开根号n比较科学吧

#include 
#include
#include
#include
#include
using namespace std;typedef long long ll;const int N=5e4+5, M=300;inline int read(){ char c=getchar();int x=0,f=1; while(c<'0'||c>'9'){
if(c=='-')f=-1; c=getchar();} while(c>='0'&&c<='9'){x=x*10+c-'0'; c=getchar();} return x*f;}int n, a[N], Q, mp[N], l, r;struct meow{
int l,r;} b[M];int block, m, pos[N];inline void ini(){ block=sqrt(n); m=(n-1)/block+1; for(int i=1;i<=n;i++) pos[i]=(i-1)/block+1; for(int i=1;i<=m;i++) b[i].l=(i-1)*block+1, b[i].r=i*block; b[m].r=n;}struct Block{ int c[N]; inline void add(int p,int v) {
for(;p<=*mp;p+=(p&-p)) c[p]+=v;} inline int sum(int p) {
int re=0; for(;p;p-=(p&-p)) re+=c[p]; return re;} int f[M][N], s[M][N]; void Set(int x){ for(int i=b[x].l; i<=n; i++) f[x][i]= f[x][i-1]+sum(*mp)-sum(a[i]), add(a[i], 1);// printf("f %d %d %d\n",x,i,f[x][i]); memset(c, 0, sizeof(c)); for(int i=b[x].l; i<=b[x].r; i++) s[x][a[i]]++; for(int i=1; i<=*mp; i++) s[x][i]+= s[x][i-1]; for(int i=1; i<=*mp; i++) s[x][i]+= s[x-1][i];// printf("s %d %d %d\n",x,i,s[x][i]);; } int Que(int l,int r){ int pl=pos[l], pr=pos[r]; int ans=0; if(pl==pr){ for(int i=l;i<=r;i++) ans+= sum(*mp)-sum(a[i]), add(a[i], 1); for(int i=l;i<=r;i++) add(a[i], -1); }else{ ans= f[pl+1][r]; //printf("ans %d \n",ans); for(int i=b[pr].l; i<=r; i++) add(a[i], 1); for(int i=b[pl].r; i>=l; i--) ans+= sum(a[i]-1) + s[pr-1][a[i]-1] - s[pl][a[i]-1], add(a[i], 1); for(int i=b[pr].l; i<=r; i++) add(a[i], -1); for(int i=l; i<=b[pl].r; i++) add(a[i], -1); } return ans; }}B;int main(){ freopen("in","r",stdin); n=read(); for(int i=1;i<=n;i++) a[i]=mp[i]=read(); sort(mp+1, mp+1+n); *mp=unique(mp+1, mp+1+n) - mp - 1; for(int i=1;i<=n;i++) a[i]=lower_bound(mp+1, mp+1+*mp, a[i]) - mp; //for(int i=1;i<=n;i++) printf("%d ",a[i]); puts(""); ini(); for(int i=1;i<=m;i++) B.Set(i); Q=read(); int lastans=0; while(Q--){ //l=read(), r=read(); l=read()^lastans, r=read()^lastans; if(l>r) swap(l, r); lastans=B.Que(l, r); printf("%d\n",lastans); }}

 

转载于:https://www.cnblogs.com/candy99/p/6579556.html

你可能感兴趣的文章
[整理]在命令行执行 UIAutomation
查看>>
Overlapping Rectangles
查看>>
新浪IP库地址
查看>>
centos7 eclispe 编译C++遇到的问题总结
查看>>
1-ML的基本概念
查看>>
自动附加Sqlserver数据库
查看>>
02.设计模式
查看>>
LeetCode 371. Sum of Two Integers
查看>>
ZBrush功能特性之法线贴图
查看>>
virtualbox虚拟机ubuntu操作系统,设置网络互通、访问,能访问虚拟机swoole的http服务...
查看>>
磁盘inode节点被占满的解决方法
查看>>
js 中对象属性的特性
查看>>
周四的工作 张孝祖
查看>>
怎么教软件工程课的看法
查看>>
mac下mysql安装及配置启动
查看>>
PAT 甲级 1020 Tree Traversals (二叉树遍历)
查看>>
安全增强 Linux (SELinux) 剖析
查看>>
防御性编程
查看>>
软件开发工程师必须知道的20个常识
查看>>
OC各种遍历方法的效率比较
查看>>