一开始我用分块大法,分成$\sqrt{n}$块,每个块上维护一个Splay,然后balabala维护一下,时间复杂度是$O(n\sqrt{n}logn)$。后来对拍的时候发现比$O(n^2)$的暴力跑得还慢,xxy学长说是Splay常数太大2333333

考试的时候没想到可以在每个块上建一个$10^5$的数组来存储每个数字出现的次数,而是用了常数巨大且复杂度多了一个log的SplayQwQ,发现自己完全没有对空间复杂度的认识啊(┙>∧<)┙へ┻┻

标算是块状链表,什么balabala比较基础地维护,卡着空间开2333333

我把块的大小设为$[\frac{\sqrt{n}}{2},\sqrt{n}×2)$,在codevs上TLE,,,

后来把块的大小改成了$[\sqrt{n},\sqrt{n}×2)$,1s内能轻松跑过。

也许是因为某些玄学的原因吧,,,

#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N = 100000;
const int B = 316;
const int BB = 632; int in() {
int k = 0, fh = 1; char c = getchar();
for(; c < '0' || c > '9'; c = getchar())
if (c == '-') fh = -1;
for(; c >= '0' && c <= '9'; c = getchar())
k = (k << 3) + (k << 1) + c - '0';
return k * fh;
} struct BLOCK {
BLOCK() {
nxt = NULL;
len = 0;
memset(times, 0, sizeof(times));
}
BLOCK *nxt;
int a[BB + B + 3], times[N + 1], len;
}; int cnt = 0;
int n; namespace BlockList {
BLOCK *head;
void Build(BLOCK * t) {head = t;}
void split(BLOCK *r) {
int rlen = r->len / 2, tlen = r->len - rlen, to = rlen;
BLOCK *t = new BLOCK;
memcpy(t->a + 1, r->a + rlen + 1, sizeof(int) * tlen);
for(int i = 1; i <= tlen; ++i) {
++t->times[t->a[i]];
--r->times[t->a[i]];
}
t->len = tlen;
r->len = rlen;
t->nxt = r->nxt;
r->nxt = t;
}
void merge(BLOCK *r) {
BLOCK *t = r->nxt;
if (t == NULL) return;
int tlen = t->len, to = r->len;
memcpy(r->a + to + 1, t->a + 1, sizeof(int) * tlen);
for(int i = 1; i <= tlen; ++i)
++r->times[t->a[i]];
r->len += tlen;
t = t->nxt;
delete r->nxt;
r->nxt = t;
if (r->len >= BB) split(r);
}
BLOCK *find(int &k) {
BLOCK *t = head;
while (k - t->len > 0 && t != NULL) {
k -= t->len;
t = t->nxt;
}
return t;
}
void work(int l, int r) {
BLOCK *t = find(r);
int num = t->a[r];
--t->times[num];
int tlen = --t->len;
for(int i = r; i <= tlen; ++i)
t->a[i] = t->a[i + 1];
if (t->len < B) merge(t);
t = find(l);
for(int i = ++t->len; i > l; --i)
t->a[i] = t->a[i - 1];
t->a[l] = num;
++t->times[num];
if (t->len >= BB) split(t);
}
int query(int l, int r, int k) {
BLOCK *tl = find(l), *tr = find(r);
int ret = 0;
if (tl == tr) {
for(int i = l; i <= r; ++i)
if (tl->a[i] == k) ++ret;
return ret;
} else {
int lentl = tl->len;
for(int i = l; i <= lentl; ++i)
if (tl->a[i] == k) ++ret;
for(int i = 1; i <= r; ++i)
if (tr->a[i] == k) ++ret;
for(tl = tl->nxt; tl != tr && tl != NULL; ret += tl->times[k], tl = tl->nxt);
return ret;
}
}
} int main() {
n = in();
int c = 0, k;
BLOCK *tmp = new BLOCK;
BlockList::Build(tmp);
for(int i = 1; i <= n; ++i) {
++c;
if (c > B) {
c = 1;
tmp->len = B;
tmp->nxt = new BLOCK;
tmp = tmp->nxt;
}
k = in();
tmp->a[c] = k;
++tmp->times[k];
}
tmp->len = c;
int m = in(), op, l, r;
while (m--) {
op = in(); l = in(); r = in();
if (op == 1)
BlockList::work(l, r);
else {
k = in();
printf("%d\n", BlockList::query(l, r, k));
}
}
return 0;
}

继续颓文化课,期末考试Bless All!

【CodeVS 5032】【省队集训2016 Day5 T1】Play with array的更多相关文章

  1. [2018湖南省队集训] 6.24 T1 marshland

    题面在这里! 一开始感觉像一个类似二分图的最小割,于是成功跑偏2333333 很容易发现一个关键性质,'L'的两个角落在的偶数格 的行(或者列)的奇偶性一定不同.... 于是我们再把偶数格按照行(或者 ...

  2. JS省队集训记

    不知不觉省队集训已经结束,离noi也越来越近了呢 论考前实战训练的重要性,让我随便总结一下这几天的考试 Day 1 T1 唉,感觉跟xj测试很像啊?meet in middle,不过这种题不多测是什么 ...

  3. HN2018省队集训

    HN2018省队集训 Day1 今天的题目来自于雅礼的高二学长\(dy0607\). 压缩包下载 密码: 27n7 流水账 震惊!穿着该校校服竟然在四大名校畅通无阻?霸主地位已定? \(7:10\)从 ...

  4. 2018HN省队集训

    HNOI2018省队集训 Day 1 流水账 T1 tree 换根+求\(lca\)+求子树和,一脸bzoj3083遥远的国度的既视感.子树和讨论一下就好了,\(lca\)?也是大力讨论一波. 先写了 ...

  5. 2017FJ省队集训 游记

    2017FJ省队集训 游记 又是一篇流水账 Day 1 今天是省队集训的第一天.早上骑车去八中,到的时候汗流太多浑身湿透被杨哥哥和runzhe2000 d了,一个说我去游泳了一个说我打球了...流完汗 ...

  6. UOJ #274. 【清华集训2016】温暖会指引我们前行 [lct]

    #274. [清华集训2016]温暖会指引我们前行 题意比较巧妙 裸lct维护最大生成树 #include <iostream> #include <cstdio> #incl ...

  7. UOJ_274_[清华集训2016]温暖会指引我们前行_LCT

    UOJ_274_[清华集训2016]温暖会指引我们前行_LCT 任务描述:http://uoj.ac/problem/274 本题中的字典序不同在于空串的字典序最大. 并且题中要求排序后字典序最大. ...

  8. UOJ 275. 【清华集训2016】组合数问题

    UOJ 275. [清华集训2016]组合数问题 组合数 $C_n^m $表示的是从 \(n\) 个物品中选出 \(m\) 个物品的方案数.举个例子,从$ (1,2,3)(1,2,3)$ 三个物品中选 ...

  9. UOJ #269. 【清华集训2016】如何优雅地求和

    UOJ #269. [清华集训2016]如何优雅地求和 题目链接 给定一个\(m\)次多项式\(f(x)\)的\(m+1\)个点值:\(f(0)\)到\(f(m)\). 然后求: \[ Q(f,n,x ...

随机推荐

  1. POJ1062昂贵的聘礼[最短路建模]

    昂贵的聘礼 Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 45892   Accepted: 13614 Descripti ...

  2. USACO GCD Extreme(II)

    题目大意:求gcd(1,2)+gcd(1,3)+gcd(2,3)+...+gcd(n-1,n) ---------------------------------------------------- ...

  3. Git/Github + TortoiseGit 使用教程

    前言 Git是一个开源的分布式版本控制系统,用以有效.高速的处理从很小到非常大的项目版本管理. 在github上有很多优秀的项目,一个伟大的学习宝库.本文分享使用tortoisegit对github/ ...

  4. lock关键字只不过是C#提供的语法糖

    lock关键字只不过是C#提供的语法糖, 最终使用的还是Monitor类. Monitor类的Enter方法要求传入的参数不为null, 否则会有ArgumentNullException excep ...

  5. BeanShell Assertion in Jmeter

    以下为几个beanshell assertion的栗子: if (ResponseCode != null && ResponseCode.equals ("200" ...

  6. win10的安装与下载

    1.下载介质创建工具 https://www.microsoft.com/zh-cn/software-download/windows10 2. 下载iso https://www.microsof ...

  7. Docker简明教程(以安装wget程序为例)

    本文计划: 一.安装Docker(Centos) 二.注册Docker官网帐号 三.下载基础centos镜像,安装需要的软件和环境后,push到自己的repository 一.安装Docker(Cen ...

  8. iis 发布MVC HTTP错误 403.14

    1. 在web.config增加<modules runAllManagedModulesForAllRequests="true" /> 如:问题依然不能解决.报错貌 ...

  9. SQL 按特定字段值排序

    SQL 按特定字段值排序的代码,有需要的朋友可以参考下. id, name shandong01 name1 shandong02 name2 shandong03 name3 beijing01 n ...

  10. maven spring profile 协同

    请在网上查相关的使用情景,这里直接上要点.另外,可能不只一种方法,但这里只有一种. 1.POM.XML片段,使web.xml文件中有关活跃spring profile的内容可以被maven自动替换 & ...