题目描述

加里敦大学有个帝国图书馆,小豆是图书馆阅览室的一个书籍管理员。他的任务是把书排成有序的,所以无序的书让他产生厌烦,两本乱序的书会让小豆产生这两本书页数的和的厌烦度。现在有n本被打乱顺序的书,在接下来m天中每天都会因为读者的阅览导致书籍顺序改变位置。因为小豆被要求在接下来的m天中至少要整理一次图书。小豆想知道,如果他前i天不去整理,第i天他的厌烦度是多少,这样他好选择厌烦度最小的那天去整理。
## 输入输出格式

输入格式:

第一行会有两个数,n,m分别表示有n本书,m天

接下来n行,每行两个数,ai和vi,分别表示第i本书本来应该放在ai的位置,这本书有vi页,保证不会有放置同一个位置的书

接下来m行,每行两个数,xj和yj,表示在第j天的第xj本书会和第yj本书会因为读者阅读交换位置

输出格式:

一共m行,每行一个数,第i行表示前i天不去整理,第i天小豆的厌烦度,因为这个数可能很大,所以将结果模10^9 +7后输出.

输入输出样例

输入样例#1: 复制

5 5

1 1

2 2

3 3

4 4

5 5

1 5

1 5

2 4

5 3

1 3

输出样例#1: 复制

42

0

18

28

48

说明

对于20%的数据,1 ≤ ai; xj; yj ≤ n ≤ 5000, m ≤ 5000, vi ≤ 10^5

对于100%的数据,1 ≤ ai; xj; yj ≤ n ≤ 50000, m ≤ 50000, vi ≤ 10^5

Solution

看一波数据范围,50000,5s。猜一波时间复杂度,\(nlogn^2\)?\(nlogn^3\)?\(nlogn^2\sqrt n\)?

想了一波树套树,发现可以做。这里我用的是树状数组套主席树。

题目大意就是给你一些数对,每个数有一个权值,也有一个厌烦度。

当一对数为逆序对,然后算总厌烦度。可以发现这题和[CQOI2011]动态逆序对很像,只不过这里有一个厌烦度,动态逆序对那道题相当于权值为1,所以这道题是那道题的升级版。

可以理解,题目中所说的本来位置,就是指判断当前这对数是否是逆序对中的标准,也就是逆序对中数的权值。

我们每次只需统计前面比他小的数字厌烦度之和是多少,并且存下有多少个比他小,设当前数字的位置为\(x\),它的本来位置为\(a[x]\)(也就是判断逆序对的权值),每个数字的厌烦度为\(w[i]\),所以当前这个数与其他数产生的厌烦度为

\[\sum_{i=1,a[x]<a[i]}^{x-1}w[i]+\sum_{i=x+1,a[x]>a[i]}^{n}w[i]
\]

但想到若两个交换的数字为逆序对,那么需要特殊判断。

还要注意细节,本题要开\(long long\),而且要注意取模和空间问题。

主席树的节点都开\(long long\)就会爆空间。

代码如下:

// luogu-judger-enable-o2
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const ll mod=1e9+7;
struct TREE
{
int ln,rn,shu;
ll zhi;
}t[15001000];
int root[210000];
ll tot,a[101000],hash[101001],n,m,id[101000],w[109101];
ll qll[20101],qrr[20100],q1,q2,ge,ff=0;
ll lowbit(ll x) {return ((x)&(-x));}
void gai(int &node,int l,int r,ll hs,ll v,ll kk)
{
if(!node) node=++tot;
t[node].zhi+=v;
t[node].shu+=kk;
if(l==r) return;
ll mid=(l+r)/2;
if(hs<=mid) gai(t[node].ln,l,mid,hs,v,kk);
else gai(t[node].rn,mid+1,r,hs,v,kk);
}
void add(ll p,ll v,ll kk){for(ll i=p;i<=n;i+=lowbit(i)) gai(root[i],1,n,hash[p],v,kk);}
ll SUM()
{
ll ans1=0,ans2=0;ff=0;
for(ll i=1;i<=q1;i++) ans1+=t[t[qrr[i]].ln].zhi,ff+=t[t[qrr[i]].ln].shu;
for(ll i=1;i<=q2;i++) ans2+=t[t[qll[i]].ln].zhi,ff-=t[t[qll[i]].ln].shu;
return ans1-ans2;
}
ll rk(ll qr,ll ql,ll l,ll r,ll k)
{
q1=0;q2=0;ll ans=0;ge=0;
for(ll i=qr;i>=1;i-=lowbit(i)) qrr[++q1]=root[i];
for(ll i=ql;i>=1;i-=lowbit(i)) qll[++q2]=root[i];
while(l<r)
{
ll mid=(l+r)/2;
if(k<=a[mid])
{
for(ll i=1;i<=q1;i++) qrr[i]=t[qrr[i]].ln;
for(ll i=1;i<=q2;i++) qll[i]=t[qll[i]].ln;
r=mid;
}
else
{
ll lsiz=SUM();
for(ll i=1;i<=q1;i++) qrr[i]=t[qrr[i]].rn;
for(ll i=1;i<=q2;i++) qll[i]=t[qll[i]].rn;
l=mid+1;ans+=lsiz;ge+=ff;
}
}
return ans;
}
ll SUM1()
{
ll ans1=0,ans2=0;ff=0;
for(ll i=1;i<=q1;i++) ans1+=t[t[qrr[i]].rn].zhi,ff+=t[t[qrr[i]].rn].shu;
for(ll i=1;i<=q2;i++) ans2+=t[t[qll[i]].rn].zhi,ff-=t[t[qll[i]].rn].shu;
return ans1-ans2;
}
ll rkk(ll qr,ll ql,ll l,ll r,ll k)
{
q1=0;q2=0;ll ans=0;ge=0;
for(ll i=qr;i>=1;i-=lowbit(i)) qrr[++q1]=root[i];
for(ll i=ql;i>=1;i-=lowbit(i)) qll[++q2]=root[i];
while(l<r)
{
ll mid=(l+r)/2;
if(k<=a[mid])
{
ll rsiz=SUM1();
for(ll i=1;i<=q1;i++) qrr[i]=t[qrr[i]].ln;
for(ll i=1;i<=q2;i++) qll[i]=t[qll[i]].ln;
r=mid;ans+=rsiz;ge+=ff;
}
else
{
for(ll i=1;i<=q1;i++) qrr[i]=t[qrr[i]].rn;
for(ll i=1;i<=q2;i++) qll[i]=t[qll[i]].rn;
l=mid+1;
}
}
return ans;
}
int main()
{
ll x,y;
long long ans=0;
cin>>n>>m;
for(ll i=1;i<=n;i++)
scanf("%lld%lld",&a[i],&w[i]),hash[i]=a[i],id[a[i]]=i;
sort(a+1,a+1+n);
for(ll i=1;i<=n;i++)
add(i,w[i],1);
for(ll i=1;i<=n;i++)
{
ans+=rkk(i-1,0,1,n,hash[i]);
ans+=rk(n,i,1,n,hash[i]);
}
for(ll i=1;i<=m;i++)
{
ll x1,y1;
scanf("%lld%lld",&x,&y);
x1=x;y1=y;
ans-=(rkk(x1-1,0,1,n,hash[x1])+ge*w[x1]);
ans-=(rk(n,x1,1,n,hash[x1])+ge*w[x1]);
ans-=(rkk(y1-1,0,1,n,hash[y1])+ge*w[y1]);
ans-=(rk(n,y1,1,n,hash[y1])+ge*w[y1]);
if((x1<y1&&hash[x1]>hash[y1])||(x1>y1&&hash[x1]<hash[y1]))
ans+=w[x1]+w[y1];
add(x1,-w[x1],-1);add(y1,-w[y1],-1);
swap(hash[x1],hash[y1]);swap(w[x1],w[y1]);
add(x1,w[x1],1);add(y1,w[y1],1);
x1=x;y1=y;
ans+=(rkk(x1-1,0,1,n,hash[x1])+ge*w[x1]);
ans+=(rk(n,x1,1,n,hash[x1])+ge*w[x1]);
ans+=(rkk(y1-1,0,1,n,hash[y1])+ge*w[y1]);
ans+=(rk(n,y1,1,n,hash[y1])+ge*w[y1]);
if((x1<y1&&hash[x1]>hash[y1])||(x1>y1&&hash[x1]<hash[y1]))
ans-=(w[x1]+w[y1]);
printf("%lld\n",ans%mod);
}
}

博主蒟蒻,可以随意转载,但必须附上原文链接k-z-j

[TJOI2017] 不勤劳的图书管理员的更多相关文章

  1. 【bzoj4889】: [Tjoi2017]不勤劳的图书管理员 分块-BIT

    [bzoj4889]: [Tjoi2017]不勤劳的图书管理员 题目大意:给定一个序列(n<=50000),每个数有一个编码ai(ai<=50000)和权值vi(vi<=100000 ...

  2. 【BZOJ4889】[Tjoi2017]不勤劳的图书管理员 分块+树状数组

    [BZOJ4889][Tjoi2017]不勤劳的图书管理员 题目描述 加里敦大学有个帝国图书馆,小豆是图书馆阅览室的一个书籍管理员.他的任务是把书排成有序的,所以无序的书让他产生厌烦,两本乱序的书会让 ...

  3. 洛谷P3759 - [TJOI2017]不勤劳的图书管理员

    Portal Description 给出一个\(1..n(n\leq5\times10^4)\)的排列\(\{a_n\}\)和数列\(\{w_n\}(w_i\leq10^5)\),进行\(m(m\l ...

  4. 【loj2639】[Tjoi2017]不勤劳的图书管理员

    #2639. 「TJOI2017」不勤劳的图书管理员 题目描述 加里敦大学有个帝国图书馆,小豆是图书馆阅览室的一个书籍管理员.他的任务是把书排成有序的,所以无序的书让他产生厌烦,两本乱序的书会让小豆产 ...

  5. [P3759][TJOI2017]不勤劳的图书管理员(分块+树状数组)

    题目描述 加里敦大学有个帝国图书馆,小豆是图书馆阅览室的一个书籍管理员.他的任务是把书排成有序的,所以无序的书让他产生厌烦,两本乱序的书会让小豆产生 这两本书页数的和的厌烦度.现在有n本被打乱顺序的书 ...

  6. BZOJ4889 & 洛谷3759:[TJOI2017]不勤劳的图书管理员——题解

    https://www.lydsy.com/JudgeOnline/problem.php?id=4889 https://www.luogu.org/problemnew/show/P3759 加里 ...

  7. [bzoj4889] [Tjoi2017]不勤劳的图书管理员

    Description 加里敦大学有个帝国图书馆,小豆是图书馆阅览室的一个书籍管理员.他的任务是把书排成有序的,所以无序的书让他产生厌烦,两本乱序的书会让小豆产生这两本书页数的和的厌烦度.现在有n本被 ...

  8. 【bzoj4889】[Tjoi2017]不勤劳的图书管理员 树状数组+分块+二分

    题目描述(转自洛谷) 加里敦大学有个帝国图书馆,小豆是图书馆阅览室的一个书籍管理员.他的任务是把书排成有序的,所以无序的书让他产生厌烦,两本乱序的书会让小豆产生这两本书页数的和的厌烦度.现在有n本被打 ...

  9. 4889: [Tjoi2017]不勤劳的图书管理员 树套树

    国际惯例的题面(Bzoj没有,洛谷找的):动态加权逆序对,一眼树套树.256MB内存,5e4范围,不虚不虚.首先把交换改成两个插入和两个删除.考虑插入和删除的贡献,就是统计前面比这个值大的数的数值和, ...

随机推荐

  1. 性能学习之--loaderunner中run-time setting常用功能

    先打开run-time setting界面 一.Run Logic   设置迭代次数,只控制action的迭代次数,init和end只执行一次. 如果迭代次数设置10,10个并发用户,那么init和e ...

  2. BZOJ 1012 [JSOI2008]最大数maxnumber【线段树】

    水题,每次记录一下当前有多少个数,然后按照题目所指示的那样模拟就行,每次向线段树末尾插入(其实是修改)题目中指定的数,然后询问当前的个数到前面Q个数中最大值是多少结果就是,好久不碰线段树了,用数组模拟 ...

  3. poj3180 The Cow Prom

    The Cow Prom Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 2373   Accepted: 1402 Desc ...

  4. Chapter 4-5

    1.切片对象  sequence[起始索引:结束索引:步进值] 对象身份的比较 is /is not 2.eval()参数是一个字符串, 可以把这个字符串当成表达式来求值. >>>x ...

  5. calc BZOJ 2655

    calc [问题描述] 一个序列a1,...,an是合法的,当且仅当: 长度为给定的n. a1,...,an都是[1,A]中的整数. a1,...,an互不相等. 一个序列的值定义为它里面所有数的乘积 ...

  6. 说说icon图标

    咳咳,其实我是想copy过来的,然而,他竟然是用代码写的图标... (正经脸)话说icon图标是一种网页中常用图标的一种,网络上有各式各样的应用案例,在此就不多啰嗦了.其实我也不过是用着现成的而已,所 ...

  7. StoryBoard中使用xib

    转自:http://blog.csdn.net/li6185377/article/details/8131042 一般自定义View       代码方式 有       在初始化的时候添加 子Vi ...

  8. D3拖动效果

    <!DOCTYPE html> <html> <head> <meta charset="UTF-8"> <title> ...

  9. .net core webapi jwt 更为清爽的认证 ,续期很简单

    我的方式非主流,控制却可以更加灵活,喜欢的朋友,不妨花一点时间学习一下 jwt认证分为两部分,第一部分是加密解密,第二部分是灵活的应用于中间件,我的处理方式是将获取token放到api的一个具体的co ...

  10. .NET Core 3.0之创建基于Consul的Configuration扩展组件

    写在前面 经过前面三篇关于.NET Core Configuration的文章之后,本篇文章主要讨论如何扩展一个Configuration组件出来.如果前面三篇文章没有看到,可以点击如下地址访问 .N ...