线段树+树状数组

感觉这东西就是每棵线段树管的区间变了,以前学的时候线段树总是只管一个点或者管(1-i)这些点,但是这东西如果加上树状数组的思想,每棵线段树管( i-(i&-i)+1 ~ i )这些区间,那么动态修改单点就不用nlog修改只用改这个点在那log棵线段树改log*log次;查询的话变慢了点,之前的 log变成了 log * log;
精品博客 讲的真的很好 感谢大佬:
博客园 TaylorSwift13 的 树状数组套权值线段树
动态逆序对
代码

Dynamic Rankings
代码

#include <iostream>
#include <algorithm>
#include <cstring>
#include <sstream>
using namespace std;
typedef long long LL;
const int N = 5e5 + 4;
int n, m;
int a[N], f[N];
int rt[N], ls[N<<6], rs[N<<6], val[N<<6], idx = 0;namespace seg
{void modify(int &u, int l, int r, int x, int f){if(!u) u = ++idx;val[u] += f;if(l==r)return ;int mid = l+r>>1;x <= mid ? modify(ls[u], l, mid, x, f) : modify(rs[u], mid+1, r, x, f);}int quary(int u, int l, int r, int L, int R){if(!u||L>r||R<l)return 0;if(L<=l&&R>=r)return val[u];int mid = l+r>>1;return quary(ls[u], l, mid, L, R) + quary(rs[u], mid+1, r, L, R);}
};namespace tree
{int quary(int u, int l, int r){int sum = 0;for(int x = u;x >= 1;x -= (x&-x))sum += seg::quary(rt[x], 0, n+1, l, r);return sum;}void modify(int u,int g, int f){for(int x = u;x <= n;x += (x&-x))seg::modify(rt[x], 0, n+1, g, f);}
};namespace tr{int tr[N]={};void modify(int u,int f){while(u<=n+1)tr[u]+=f,u+=(u&-u);}int quary(int u){int sum=0;while(u)sum+=tr[u],u-=(u&-u);return sum;}
}int main() 
{scanf("%d%d",&n, &m);for(int i = 1;i <= n;i ++)scanf("%d",a + i);LL ans = 0;for(int i = 1;i <= n;i ++){ans += tree::quary(i-1, a[i]+1, n+1);tree::modify(i, a[i], 1);tr::modify(a[i], 1);f[a[i]] = i;}while(m --){int x;scanf("%d", &x);	cout<<ans<<endl;ans -= tree::quary(f[x]-1, x+1, n+1);ans -= tr::quary(x-1) - tree::quary(f[x]-1, 0, x-1);tree::modify(f[x], x, -1);tr::modify(x, -1);}return 0;
}
#include <iostream>
#include <algorithm>
#include <cstring>
#include <sstream>
using namespace std;typedef long long ll;
const int N = 1e5 + 4;
int n, m;
int rt[N], ls[N<<9], rs[N<<9], val[N<<9], idx;
int x[N], cntx, y[N], cnty, a[N];int quary(int l, int r, int k)
{	if(l == r)return l;int sum = 0;for(int i = 0;i < cntx ;i ++)sum -= val[ls[x[i]]];for(int i = 0;i < cnty ;i ++)sum += val[ls[y[i]]];if(sum >= k){for(int i = 0;i < cntx ;i ++)x[i] = ls[x[i]];for(int i = 0;i < cnty ;i ++)y[i] = ls[y[i]];return quary(l, l+r>>1, k);}else {for(int i = 0;i < cntx ;i ++)x[i] = rs[x[i]];for(int i = 0;i < cnty ;i ++)y[i] = rs[y[i]];return quary((l+r>>1) + 1, r, k - sum);}
}void  modify(int &u, int l, int r, int x, int f)
{if(!u)u = ++idx;val[u] += f;if(l == r)return ;int mid = l+r>>1;x <= mid ? modify(ls[u], l, mid, x, f) : modify(rs[u], mid+1, r, x, f);
}int main() 
{	scanf("%d%d", &n, &m);for(int i = 1;i <= n;i ++){scanf("%d", a+i);for(int u = i;u <= n;u +=(u&-u))modify(rt[u], 0, 1e9, a[i], 1);}while(m --){char c[2];int l, r, k;scanf("%s%d%d", c, &l, &r);if(c[0]=='Q'){scanf("%d", &k);cntx = cnty = 0;for(int j = l-1;j >= 1;j -= (j&-j)) x[cntx++] = rt[j];for(int u = r;u >= 1; u -= (u&-u)) y[cnty++] = rt[u];cout<<quary(0, 1e9, k)<<endl;}else{for(int u = l;u <= n;u +=(u&-u))modify(rt[u], 0, 1e9, a[l], -1);a[l] = r;for(int u = l;u <= n;u +=(u&-u))modify(rt[u], 0, 1e9, a[l], 1);}}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/316183.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

cf1009F. Dominant Indices

cf1009F. Dominant Indices 题意&#xff1a; 1号节点为根&#xff0c;问对于以每个点为根的子树种&#xff0c;求某个深度节点最多的层数。如果多个相等&#xff0c;输出深度小的。 题解&#xff1a; 直接dsu就完事了&#xff0c;相当于是求子树的众数&#xff0c;但是注意…

Asp.Net Core Web应用程序—探索

前言作为一个Windows系统下的开发者&#xff0c;我对于Core的使用机会几乎为0&#xff0c;但是考虑到微软的战略规划&#xff0c;我觉得&#xff0c;Core还是有先了解起来的必要。因为&#xff0c;目前微软已经搞出了两个框架了&#xff0c;一个是Net标准(.NetFramework)&#…

智能优化算法应用:基于跳蛛算法3D无线传感器网络(WSN)覆盖优化 - 附代码

智能优化算法应用&#xff1a;基于跳蛛算法3D无线传感器网络(WSN)覆盖优化 - 附代码 文章目录 智能优化算法应用&#xff1a;基于跳蛛算法3D无线传感器网络(WSN)覆盖优化 - 附代码1.无线传感网络节点模型2.覆盖数学模型及分析3.跳蛛算法4.实验参数设定5.算法结果6.参考文献7.MA…

Transformation HDU - 4578

思路&#xff1a;三种改变的顺序&#xff0c;先变 后乘 再 加。 #include <iostream> #include <algorithm> #include <cstring> #include <sstream> using namespace std; typedef long long LL; const int N 5e5 4, mod 1e4 7; int n, m;struct…

cf375D. Tree and Queries

cf375D. Tree and Queries 题意&#xff1a; 给你一颗有根树&#xff0c;每个点都有一个颜色&#xff0c;有m次询问&#xff0c;问以u为根的子树中&#xff0c;相同颜色数量超过k的有多少种颜色&#xff1f; 题解&#xff1a; 这个题做法很多&#xff0c;有莫队分块&#x…

让你的ASP.NET Core应用程序更安全

对于ASP.NET Core应用程序&#xff0c;除了提供认证和授权机制来保证服务的安全性&#xff0c;还需要考虑下面的一些安全因素&#xff1a;CSRF强制HTTPS安全的HTTP HeadersCSRFASP.NET Core通过AntiForgeryToken来阻止CSRF攻击&#xff0c;一般来说&#xff0c;当用户做表单提交…

201409-5 拼图

七十分超时代码&#xff0c;具体思路可以仿照acwing蒙德里安的梦想 就是枚举第i层放 a &#xff0c;第i-1层放b的合法方案其中b和a都是同一块&#xff1b; #include <iostream> #include <algorithm> #include <cstring> #include <sstream> using nam…

.net core 中间件管道底层剖析

.net core 管道&#xff08;Pipeline&#xff09;是什么&#xff1f;由上图可以看出&#xff0c;.net core 管道是请求抵达服务器到响应结果返回的中间的一系列的处理过程&#xff0c;如果我们简化一下成下图来看的话&#xff0c;.net core 的管道其实就是中间件的部分。微软中…

Codeforces Round #741 (Div. 2)

Codeforces Round #741 (Div. 2) 题号题目知识点AThe Miracle and the SleeperBScenes From a MemoryCRingsD1Two Hundred Twenty One (easy version)D2Two Hundred Twenty One (hard version)ERescue Niwen!FTubular Bells

201609-5 祭坛

超时暴力60分 #include <iostream> #include <algorithm> #include <cstring> #include <sstream> using namespace std; typedef long long LL; const int N 5e5 4, mod 1e9 7;vector<int> v[N]; int y[N], sizy, x[N], sizx, numx[N], a[…

cf1562E. Rescue Niwen!

cf1562E. Rescue Niwen! 题意&#xff1a; 我们定义一个字符串s1s2s3…sn的展开为:s1,s1s2,s1s2…sn,s2,s2s3,s2s3…sn,…,sn 找到字符串s“扩展”的最大递增子序列的大小(根据字典序大小比较) 题解&#xff1a; 第一感觉就是求最长上升子序列的变形 按照字典序大小比较规则…

架构为什么要以领域为核心

很久以前, 人们以为地球是太阳系的中心.然后一位聪明人, 哥白尼, 他改变了我们对太阳系的看法. 他认为太阳是太阳系的中心:这是对太阳更好的一个解释, 更简单也更具说服力. 事实上, 以太阳为中心的模型确实是更优雅的.上面这件事也发生在软件开发里. 下面这个就是很多开发者惯用…

201403-5 任务调度

哇&#xff0c;ccf csp认证考试 历年真题解&#xff08;一本书&#xff09;真厉害。 #include<iostream> #include<cstdio> #include<algorithm> #include<cstring>using namespace std;typedef long long LL; typedef pair<int,int> PII; con…

cf1562D Two Hundred Twenty One

cf1562D Two Hundred Twenty One 题意&#xff1a; 定义一个前缀和公式&#xff1a;a1−a2a3−a4..∑i1n(−1)i−1∗aia_{1}-a_{2}a_{3}-a_{4}..\sum_{i1}^{n}(-1)^{i-1}*a_{i}a1​−a2​a3​−a4​..∑i1n​(−1)i−1∗ai​ 然后给你一个长度为n的序列&#xff0c;只包含{-1…

C# - 为引用类型重定义相等性 - 继承相关

派生类这是上面Citizen类的一个子类&#xff1a;下面我重写object.Equals() 方法&#xff1a;大部分逻辑都在base.Equals()方法里了&#xff0c;首先如果父类的Equals()方法返回false&#xff0c;那么下面也就不用做啥了。但是如果父类Equals()认为这两个实例是相等的&#xff…

201612-5 卡牌游戏

**根据题目样例解释得到每种卡牌拥有状态之间的关系&#xff0c;然后转换成等式&#xff0c;高斯消元是2^(3n) ** 80分超时代码&#xff1a; #include <iostream> #include <algorithm> #include <cstring> #include <vector> #include <queue>…

cf1562 C. Rings

cf1562 C. Rings 题意&#xff1a; 给你一个长度为n的01串&#xff0c;在01串选两个连续区间S和T&#xff0c;要求区间长度>⌊n2⌋\lfloor \frac{n}{2} \rfloor⌊2n​⌋。 现在定义一个函数f(S):将S01串以2二进制转化成10进制&#xff0c;要求f(S)是f(T)的倍数 题解&…

微软百名员工签名力挺996.ICU

中国程序员上传到 GitHub 的 996.ICU repo 火速在互联网广泛传播时&#xff0c;996 工作制引起了全球的广泛关注&#xff0c;Python 之父直指这是不人道的行为&#xff0c;事情经过不断发酵&#xff0c;中国官方媒体也接连发声表态要警惕「996 工作制」。就在今日&#xff0c;微…

201512-5 矩阵

只会暴力&#xff0c;答案没看懂&#xff0c;呜呜呜。乘的时候还乱七八糟的。 #include <iostream> #include <algorithm> #include <cstring> #include <vector> #include <queue> #include <bitset> #define ls (u<<1) #define …

P1174 打砖块

P1174 打砖块 题意&#xff1a; 题解&#xff1a; 参考题解&#xff1a; I_AM_HelloWord danxmz2006 这两个博客结合看&#xff0c;大致就能理解 我们只在N处转移&#xff0c;面对Y类的块无需决策&#xff0c;因为Y类的块可以一直打 不同的打砖块的顺序&#xff0c;决定了我…