[TJOI2017]不勤劳的图书管理员

最新推荐文章于 2021-02-22 22:44:50 发布
最新推荐文章于 2021-02-22 22:44:50 发布
阅读量363 收藏
点赞数
分类专栏： 数据结构
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Richard_for_OI/article/details/81035771
版权

Solution1

分块+块内树状数组。可以踩掉标算，甚至碾压。
#include <cstdio>
#include <cstring>
#include <cmath>
#define rint register int
#define N 50010
#include <algorithm>
typedef long long ll;
using namespace std;
const ll mo = 1e9 + 7;
inline char gc() {
    static char now[1<<16], *S, *T;
    if(S == T) {T = (S = now) + fread(now, 1, 1<<16, stdin); if(S == T) return EOF;}
    return *S++;
}
inline int read() {
    int x = 0, f = 1; char c = gc();
    while(c < '0' || c > '9') {if(c == '-') f = -1; c = gc();}
    while(c >= '0' && c <= '9') {x = x * 10 + c - 48; c = gc();}
    return x * f;
}
inline void print(int x) {
    if(x == 0) {puts("0"); return ;}
    if(x < 0) putchar('-'), x = 0 - x;
    int d[30], now = -1;
    while(x) {d[++now] = x % 10; x/= 10;}
    for(int i = now; i >= 0; --i) putchar('0' + d[i]); puts("");
}
int n, m, block, num;
int belong[N], q[N], mp[N]; ll s[500][200];
struct node {int x, y;}a[N], b[500];
inline void build(int x) {
    for(rint i = b[x].x; i <= b[x].y; ++i) q[i] = a[i].x;
    sort(q+b[x].x, q+b[x].y+1); for(rint i = 1; i <= b[x].y - b[x].x + 1; ++i) s[x][i] = 0;
    for(rint i = b[x].x; i <= b[x].y; ++i) s[x][i - b[x].x + 1] = (s[x][i - b[x].x] + mp[q[i]]) % mo;
}
ll bit1[N], ans = 0, bit2[N];
inline void add1(int x, ll v) {for(rint i = x; i <= n; i+= i & -i) bit1[i]+= v, bit1[i]%= mo;}
inline ll query1(int x) {ll ret = 0; for(rint i = x; i; i-= i & -i) ret+= bit1[i], ret%= mo; return ret;}
inline void add2(int x, ll v) {for(rint i = x; i <= n; i+= i & -i) bit2[i]+= v, bit2[i]%= mo;}
inline ll query2(int x) {ll ret = 0; for(rint i = x; i; i-= i & -i) ret+= bit2[i], ret%= mo; return ret;}
inline void init() {
    for(rint i = 1; i <= n; ++i) bit1[i] = bit2[i] = 0;
    for(rint i = n; i >= 1; --i) {
        ans+= query1(a[i].x - 1) + query2(a[i].x) * a[i].y; ans%= mo; if(ans < 0) ans+= mo;
        add1(a[i].x, a[i].y); add2(a[i].x, 1);
    }
}
int main() {
    n = read(); m = read();
    for(rint i = 1; i <= n; ++i) {a[i].x = read(); a[i].y = read(); mp[a[i].x] = a[i].y;}
    block = (int)sqrt(n)/2; num = 0;
    for(rint i = 1; i <= n; i+= block) {
        b[++num].x = i; b[num].y = min(n, i + block - 1);
        for(rint j = i; j <= b[num].y; ++j) belong[j] = num;
    }
    for(rint i = 1; i <= num; ++i) build(i); init();
    for(rint i = 1; i <= m; ++i) {
        int x = read(), y = read(); if(x > y) swap(x, y); if(x == y) {print(ans); continue;}
        int b1 = belong[x], b2 = belong[y];
        for(rint j = b1 + 1; j <= b2 - 1; ++j) {
            int sz = b[j].y - b[j].x + 1;
            int p1 = lower_bound(q+b[j].x, q+b[j].y+1, a[x].x) - q;
            p1-= b[j].x; ll temp1 = s[j][p1]; if(temp1 < 0) temp1+= mo;
            ll ctb1 = (p1 * a[x].y + temp1) % mo;
            int p2 = lower_bound(q+b[j].x, q+b[j].y+1, a[y].x) - q;
            p2-= b[j].x; ll temp2 = s[j][p2]; if(temp2 < 0) temp2+= mo;
            ll ctb2 = (p2 * a[y].y + temp2) % mo;
            ans+= ctb2 - ctb1; if(ans < 0) ans+= mo;
            ll temp = s[j][sz]; if(temp < 0) temp+= mo;
            ctb1 = ((sz - p1) * a[x].y + temp - temp1) % mo;
            ctb2 = ((sz - p2) * a[y].y + temp - temp2) % mo;
            ans+= ctb1 - ctb2; if(ans < 0) ans+= mo;
        }
        ll ad = 0;
        for(rint j = x + 1; j <= min(y - 1, b[b1].y); ++j) {
            if(a[x].x > a[j].x) ad-= a[x].y + a[j].y, ad%= mo; else ad+= a[x].y + a[j].y, ad%= mo; if(ad < 0) ad+= mo;
            if(a[y].x > a[j].x) ad+= a[y].y + a[j].y, ad%= mo; else ad-= a[y].y + a[j].y, ad%= mo; if(ad < 0) ad+= mo;
        }
        ll ad1 = ad;
        for(rint j = max(x + 1, b[b2].x); j <= y - 1; ++j) {
            if(a[j].x > a[y].x) ad-= a[y].y + a[j].y, ad%= mo; else ad+= a[y].y + a[j].y, ad%= mo; if(ad < 0) ad+= mo;
            if(a[j].x > a[x].x) ad+= a[x].y + a[j].y, ad%= mo; else ad-= a[x].y + a[j].y, ad%= mo; if(ad < 0) ad+= mo;
        }
        if(b1 == b2) ans+= ad1; else ans+= ad; ans%= mo;
        swap(a[x], a[y]); build(b1); build(b2);
        if(a[x].x < a[y].x) ans-= a[y].y + a[x].y; else ans+= a[y].y + a[x].y; if(ans < 0) ans+= mo; ans%= mo;
        print(ans);
    }
    return 0;
}
洛谷(O2)评测结果：用时: 7676ms / 内存: 4085KB

Solution2

树套树做法。标算。大概我天生常数大。
#include <cstdio>
#include <cstring>
#include <vector>
using namespace std;
inline char gc() {
    static char now[1<<16], *S, *T;
    if(S == T) {T = (S = now) + fread(now, 1, 1<<16, stdin); if(S == T) return EOF;}
    return *S++;
}
inline int read() {
    int x = 0; char c = gc();
    while(c < '0' || c > '9') c = gc();
    while(c >= '0' && c <= '9') {x = x * 10 + c - 48; c = gc();}
    return x;
}
void print(int x) {
    if(x/10) print(x/10);
    putchar('0' + x % 10);
}
const int N = 50010;
const int mod = 1e9 + 7;
typedef long long ll;
int n, m, a[N], val[N]; ll ans = 0;
ll bit[N];
inline void addbit(int x, int V) {for(; x <= n; x+= x & -x) bit[x]+= V, bit[x]%= mod;}
inline ll askbit(int x) {ll ret = 0; for(; x; x-= x & -x) ret+= bit[x], ret%= mod; return ret;}
inline void init() {
    for(int i = n; i; --i) {
        ans = (ans + askbit(a[i] - 1)) % mod;
        addbit(a[i], val[i]);
    }
    memset(bit, 0, sizeof(bit));
    for(int i = 1; i <= n; ++i) {
        ans = (ans + ((askbit(n) - askbit(a[i])) % mod + mod) % mod) % mod;
        addbit(a[i], val[i]);
    }
}
int rt[N], len = 0, L[N*256], R[N*256], cnt[N*256]; ll v[N*256]; vector<int> st;
inline int newnode() {
    int ret;
    if(st.begin() != st.end()) {ret = st.back(); st.pop_back();}
    else ret = ++len;
    return ret;
}
void add(int &p, int l, int r, int pos, int x) {
    if(!p) p = newnode();
    v[p]+= x; v[p]%= mod; if(v[p] < 0) v[p]+= mod;
    cnt[p]+= (x > 0) ? (1) : (-1);
    if(l < r) {
        int mid = (l + r)>>1;
        if(pos <= mid) add(L[p], l, mid, pos, x);
        else add(R[p], mid + 1, r, pos, x);
    }
    if(!cnt[p]) {st.push_back(p); cnt[p] = v[p] = L[p] = R[p] = 0; p = 0;}
}
vector<int> b1, b2;
inline void getBIT(int x, int y) {
    b1.clear(); b2.clear();
    for(; y; y-= y & -y) b1.push_back(rt[y]);
    --x; for(; x; x-= x & -x) b2.push_back(rt[x]);
}
ll getV(int p, int l, int r, int x, int y) {
    if(!p) return 0;
    if(x <= l && r <= y) return v[p];
    int mid = (l + r)>>1;
    ll sum = 0;
    if(x <= mid) sum+= getV(L[p], l, mid, x, y), sum%= mod;
    if(mid + 1 <= y) sum+= getV(R[p], mid + 1, r, x, y), sum%= mod;
    return sum;
}
ll getC(int p, int l, int r, int x, int y) {
    if(!p) return 0;
    if(x <= l && r <= y) return cnt[p];
    int mid = (l + r)>>1;
    int sum = 0;
    if(x <= mid) sum+= getC(L[p], l, mid, x, y);
    if(mid + 1 <= y) sum+= getC(R[p], mid + 1, r, x, y);
    return sum;
}
int main() {
    n = read(); m = read();
    for(int i = 1; i <= n; ++i) {
        a[i] = read(); val[i] = read();
        for(int j = i; j <= n; j+= j & -j) {
            add(rt[j], 1, n, a[i], val[i]);
        }
    }
    init();
    for(int i = 1; i <= m; ++i) {
        int x = read(), y = read(); if(x > y) swap(x, y); if(x == y) {print(ans); puts(""); continue;}
        if(a[x] < a[y]) ans = (ans + val[x] + val[y]) % mod; else ans = ((ans - val[x] - val[y]) % mod + mod) % mod;
        getBIT(x + 1, y - 1);
        for(vector<int>::iterator it = b1.begin(); it != b1.end(); ++it) {
            ans = ((ans - (getV(*it, 1, n, 1, a[x] - 1) + getC(*it, 1, n, 1, a[x] - 1) * val[x] % mod) % mod) % mod + mod) % mod;
            ans = (ans + (getV(*it, 1, n, a[x] + 1, n) + getC(*it, 1, n, a[x] + 1, n) * val[x] % mod) % mod) % mod;

            ans = ((ans - (getV(*it, 1, n, a[y] + 1, n) + getC(*it, 1, n, a[y] + 1, n) * val[y] % mod) % mod) % mod + mod) % mod;
            ans = (ans + (getV(*it, 1, n, 1, a[y] - 1) + getC(*it, 1, n, 1, a[y] - 1) * val[y] % mod) % mod) % mod;
        }
        for(vector<int>::iterator it = b2.begin(); it != b2.end(); ++it) {
            ans = ((ans - (getV(*it, 1, n, a[x] + 1, n) + getC(*it, 1, n, a[x] + 1, n) * val[x] % mod) % mod) % mod + mod) % mod;
            ans = (ans + (getV(*it, 1, n, 1, a[x] - 1) + getC(*it, 1, n, 1, a[x] - 1) * val[x] % mod) % mod) % mod;

            ans = ((ans - (getV(*it, 1, n, 1, a[y] - 1) + getC(*it, 1, n, 1, a[y] - 1) * val[y] % mod) % mod) % mod + mod) % mod;
            ans = (ans + (getV(*it, 1, n, a[y] + 1, n) + getC(*it, 1, n, a[y] + 1, n) * val[y] % mod) % mod) % mod;
        }
        print(ans); puts("");
        for(int j = x; j <= n; j+= j & -j) add(rt[j], 1, n, a[x], -val[x]);
        for(int j = y; j <= n; j+= j & -j) add(rt[j], 1, n, a[y], -val[y]);
        swap(a[x], a[y]); swap(val[x], val[y]);
        for(int j = x; j <= n; j+= j & -j) add(rt[j], 1, n, a[x], val[x]);
        for(int j = y; j <= n; j+= j & -j) add(rt[j], 1, n, a[y], val[y]);
    }
    return 0;
}
洛谷(O2)评测结果：用时: 14720ms / 内存: 75117KB

然而bzoj就不那么友好了。

分块：3448 kb 41460 ms。
树套树：251920 kb 100000+ ms TLE
Richard_for_OI
关注
专栏目录
洛谷P3759 [TJOI2017]不勤劳的图书管理员(二维线段树)
p_b_p_b的博客
06-01 495
裸的动态逆序对,但维护的东西很多写的一不是一般的麻烦采用二维线段树,外层区间,内层权值query的返回值是pair（或用全局变量）代码里有注释// P3759 [TJOI2017]不勤劳的图书管理员 //二维线段树 #include&lt;bits/stdc++.h&gt; #define Pli pair&lt;ll,int&gt; #define mp make_pair #define f...
【洛谷3759】[TJOI2017] 不勤劳的图书管理员（树套树）
CXR的博客
10-27 344
点此看题面 大致题意： 给定一个序列，每个元素有两个属性aia_iai​和viv_ivi​，每次操作改变两个元素的位置，求每次操作后∑vi+vj[i&amp;amp;lt;j,ai&amp;amp;gt;aj]\sum{v_i+v_j}[i&amp;amp;lt;j,a_i&amp;amp;gt;a_j]∑vi​+vj​[i&amp;lt;j,ai​&amp;gt;aj​]。 关于题意的解读 其实，题目差不多就是让我们求逆序对（只不过每个逆序...
bzoj4889 [TJOI2017]不勤劳的图书管理员
Elijahqi
04-20 285
http://www.elijahqi.win/archives/3142 题目描述 加里敦大学有个帝国图书馆，小豆是图书馆阅览室的一个书籍管理员。他的任务是把书排成有序的，所以无序的书让他产生厌烦，两本乱序的书会让小豆产生这两本书页数的和的厌烦度。现在有n本被打乱顺序的书，在接下来m天中每天都会因为读者的阅览导致书籍顺序改变位置。因为小豆被要求在接下来的m天中至少要整理一次图书。小豆想知道，...
[BZOJ4889][Tjoi2017]不勤劳的图书管理员-分块-树状数组
「嫌われ者のフィロソフィ」
01-01 506
不勤劳的图书管理员Description加里敦大学有个帝国图书馆，小豆是图书馆阅览室的一个书籍管理员。他的任务是把书排成有序的，所以无序的书让他产生厌烦，两本乱序的书会让小豆产生这两本书页数的和的厌烦度。现在有n本被打乱顺序的书，在接下来m天中每天都会因为读者的阅览导致书籍顺序改变位置。因为小豆被要求在接下来的m天中至少要整理一次图书。小豆想知道，如果他前i天不去整理，第i天他的厌烦度是多少，这样他
bzoj4889 [Tjoi2017]不勤劳的图书管理员（线段树套权值线段树）
Icefox的博客
04-17 395
相当于动态维护带权逆序对。 交换x，y位置的数，只会对他们之间的数产生影响。分类讨论即可。我写的线段树套权值线段树，卡空间，需要写垃圾回收qaq 复杂度O(nlog2n)O(nlog^2n)
BZOJ4889: [TJOI2017]不勤劳的图书管理员
weixin_30917213的博客
09-13 82
BZOJ4889: [TJOI2017]不勤劳的图书管理员 额，这个题由于被卡评测了，于是变成了权限题。。。 本蒟蒻表示没钱氪金。。。 当然，可以去洛谷/$LOJ$搞搞事。。。 洛谷P3759 [TJOI2017]不勤劳的图书管理员 LOJ#2639. 「TJOI2017」不勤劳的图书管理员 这里附上洛谷的题面。 题目描述 加里敦大学有个帝国图书馆，小豆是图书馆阅览室的一个书籍管理...
NOI2017天津市选拔赛(TJOI2017)题目及测试数据
07-04
【标题】"NOI2017天津市选拔赛(TJOI2017)题目及测试数据"涉及的是中国信息学奥林匹克竞赛（National Olympiad in Informatics, NOI）的一个地方选拔赛，即2017年天津赛区的比赛。NOI是由中国计算机学会主办的一项...
洛谷P3758 [TJOI2017]可乐
qq_35827379的博客
02-22 251
原题链接 题目描述 加里敦星球的人们特别喜欢喝可乐。因而，他们的敌对星球研发出了一个可乐机器人，并且放在了加里敦星球的 11 号城市上。这个可乐机器人有三种行为： 停在原地，去下一个相邻的城市，自爆。它每一秒都会随机触发一种行为。现在给加里敦星球城市图，在第 00 秒时可乐机器人在 11 号城市，问经过了 tt 秒，可乐机器人的行为方案数是多少？ 输入格式 第一行输入两个正整数 NN，MM。NN 表示城市个数，MM 表示道路个数。 接下来 MM 行每行两个整数 uu，vv，表示 uu，vv 之间有一条道路。
TJOI 2017 (Pollard_rho,MR素数测试）
WedsonLin的博客
12-13 165
TJOI 2017 (Pollard_rho,MR素数测试） 题目链接 题意是求分数BA\frac{B}{A}AB​在模M下的逆元，没有则输出-1。 本题的两大核心是：Pollard_rho算法，和MR素数测试算法。 BA\frac{B}{A}AB​在模M下有逆元当且仅当gcd(A,M)=1gcd(A,M)=1gcd(A,M)=1，因为分数不是最简的，所以只要分数可以化简到B′A′\frac{B'}{A'}A′B′​，其中gcd(A′,M)=1gcd(A',M)=1gcd(A′,M)=1即有逆元。 我一开始
python画龙舟_BZOJ4891 TJOI2017龙舟（Polllard-Rho）
weixin_39688750的博客
12-08 759
对给定模数分解质因数后约分即可。依然常数巨大过不了。#include#include#include#include#include#includeusing namespace std;#define ll long long#define N 10010char getc(){char c=getchar();while ((c'Z')&&(c'z')&&(c'...
【loj2639】[Tjoi2017]不勤劳的图书管理员
weixin_34220963的博客
10-16 123
#2639. 「TJOI2017」不勤劳的图书管理员 题目描述 加里敦大学有个帝国图书馆，小豆是图书馆阅览室的一个书籍管理员。他的任务是把书排成有序的，所以无序的书让他产生厌烦，两本乱序的书会让小豆产生这两本书页数的和的厌烦度。现在有nnn本被打乱顺序的书，在接下来mmm天中每天都会因为读者的阅览导致书籍顺序改变位置。因为小豆被要求在接下来的mmm天中至少要整理一次...
bzoj4889 [Tjoi2017]不勤劳的图书管理员
a6t2007的博客
04-17 110
题目描述 题解： 很考验读题能力的一道题…… 首先我们知道要求的是一个类似逆序对的东西，统计时要加上两端权值。 考虑先把初始$ans$求出来，后续操作只会改变$[l,r]$内部关系，这里分块处理。 对于$ans$有影响的有：与端点有关的逆序对形成/破坏次数，已经中间节点与两端关系。 块内用树状数组维护小于等于$x$的有几个，还有小于等于$x$的权值和是多少。 再维护一下块内...
不勤劳的图书管理员
StaroForgin的博客
12-09 205
不勤劳的图书管理员 题解 很明显的一道动态逆序对的题，我们可以通过在每次交换时维护它的两点对交换的贡献得出答案。 不过如果我们只用的方法暴力去维护的话就只能得到20分，于是我们开始想优化。 有些大佬们用分块+CDQ分治秒了这道题，比笔者的树套树快了好多，不过笔者也只会打树状树组套线段树。 好了，讲一下树套树怎么打这道题。 我们内层用线段树维护它的权值，外层用树状数组维护区间。每次查询时...
[BZOJ4889][洛谷P3759][TJOI2017]不勤劳的图书管理员 分块+树状数组
TYB的博客
10-17 449
分块+树状数组
[TJOI2017]不勤劳的图书管理员 数列分块+二分
Exception
04-26 156
原题:https://loj.ac/problem/2639 题解:就是求动态逆序对的和。考虑暴力的方法，对于区间[L,R] 把[L+1.R-1]的元素和L,R比，其实就是要求[L+1,R-1]中所有大于L的和，所有小于L的和，很容易想到数列分块，把每一个块排序，二分查找求和。要注意不要用STL,注意常数优化。 #include<bits/stdc++.h> #define ll...
[TJOI2017] 不勤劳的图书管理员
weixin_30378623的博客
08-14 92
题目描述 加里敦大学有个帝国图书馆，小豆是图书馆阅览室的一个书籍管理员。他的任务是把书排成有序的，所以无序的书让他产生厌烦，两本乱序的书会让小豆产生这两本书页数的和的厌烦度。现在有n本被打乱顺序的书，在接下来m天中每天都会因为读者的阅览导致书籍顺序改变位置。因为小豆被要求在接下来的m天中至少要整理一次图书。小豆想知道，如果他前i天不去整理，第i天他的厌烦度是多少，这样他好选择厌烦度最小...
luogu3759 [TJOI2017]不勤劳的图书管理员
dianjiaxian1205的博客
03-16 107
分块+权值逆序对 #include <algorithm> #include <iostream> #include <cstdio> #include <cmath> using namespace std; typedef long long ll; int n, m, uu, vv, blc, bel[50005], ans, ccn...
P3853 [TJOI2007] 路标设置
最新发布
07-16
题目描述 在一条公路上，有 $n$ 个路标，第 $i$ 个路标的位置为 $x_i$，每个路标上都有一个数字 $a_i$。现在需要在这条公路上设置 $m$ 个路障，每个路障的位置为 $y_i$，每个路障会覆盖一段长度为 $l_i$ 的区间 $[y_i-l_i,y_i+l_i]$，其中 $l_i$ 是正整数。 对于每个路障，它会覆盖的区间内的路标的数字 $a_i$ 的和为 $s_i$，即 $s_i=\sum_{j=1}^n[a_j\in[y_i-l_i,y_i+l_i]]a_j$。 现在需要在不超过 $k$ 个路障的前提下，最大化所有路障覆盖的数字和的最小值。 输入格式 第一行包含三个整数 $n,m,k$。 接下来 $n$ 行，每行包含两个整数 $x_i,a_i$。 接下来 $m$ 行，每行包含两个整数 $y_i,l_i$。 输出格式 输出一个实数，表示最大化所有路障覆盖的数字和的最小值，保留 $2$ 位小数。 数据范围 $1\leq n,m\leq 200$, $1\leq k\leq 10$, $\leq x_i,y_i\leq 10^9$, $1\leq a_i,l_i\leq 10^6$ 输入样例1 5 2 1 1 2 2 3 3 4 4 5 5 6 3 1 7 1 输出样例1 3.50 输入样例2 5 2 2 1 2 2 3 3 4 4 5 5 6 3 1 7 1 输出样例2 2.50 算法1 (二分答案+贪心) $O(m^2\log w)$ 二分答案+贪心 时间复杂度 参考文献 python3 代码 C++ 代码 算法2 (暴力枚举) $O(n^2m^2)$ blablabla 时间复杂度 参考文献 C++ 代码
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值