HDU-5792 World is Exploding(树状数组)
题目大意:给一个整数序列,统计四元组(a,b,c,d)的个数,满足条件1:a<>b<>c<>d;条件2:<a,b>组成一个顺序对,<c,d>组成一个逆序对。(a、b、c、d均为下标)
代码如下:从所有四元组中减去不满足条件的四元组。用顺序对数乘以逆序对数得到只满足条件2的四元组数目sum,从sum减去不满足条件1的四元组数目便是答案。sum中只有四种不满足条件1的情况,即:a=c、a=d、b=c、b=d。四种情况的含义分别为顺序对的左端点与逆序对的左端点相同、顺序对的左端点与逆序对的右端点相同、顺序对的右端点与逆序对的左端点相同、顺序对的右端点与逆序对的右端点相同。维护四个数组A、B、C、D,数组中的位置 i 分别表示以 i 为右端点的顺序、逆序以及以 i 为左端点的顺序、逆序对数,便可计算出相应情况下多余的四元组数目。
这四个数组是通过维护树状数组获得的,但是通过维护线段树却超时啦。。。
代码如下:
# include<iostream>
# include<cstdio>
# include<map>
# include<vector>
# include<cstring>
# include<algorithm>
using namespace std;
# define LL long long const int N=50000; int a[N+5];
map<int,int>mp;
vector<int>v;
int vis[N+5];
int cnt[N+5];
int A[N+5];//以i为右端点的顺序对个数
int B[N+5];//以i为右端点的逆序对个数
int C[N+5];//以i为左端点的顺序对个数
int D[N+5];//以i为左端点的逆序对个数 int lowbit(int x)
{
return x&(-x);
} int query(int x)
{
int res=0;
while(x>=1){
res+=cnt[x];
x-=lowbit(x);
}
return res;
} void update(int x,int n)
{
while(x<=n){
++cnt[x];
x+=lowbit(x);
}
} int main()
{
int n;
while(~scanf("%d",&n))
{
v.clear();
mp.clear();
for(int i=1;i<=n;++i){
scanf("%d",a+i);
v.push_back(a[i]);
}
sort(v.begin(),v.end());
int m=unique(v.begin(),v.end())-v.begin();
for(int i=0;i<m;++i) mp[v[i]]=i+1; memset(cnt,0,sizeof(cnt));
memset(vis,0,sizeof(vis));
LL sum1=0,sum2=0;
for(int i=1;i<=n;++i){
A[i]=query(mp[a[i]]-1);
B[i]=i-1-A[i]-vis[mp[a[i]]];
update(mp[a[i]],m+1);
++vis[mp[a[i]]];
sum1+=A[i];
sum2+=B[i];
} memset(cnt,0,sizeof(cnt));
memset(vis,0,sizeof(vis));
for(int i=n;i>=1;--i){
D[i]=query(mp[a[i]]-1);
C[i]=n-i-D[i]-vis[mp[a[i]]];
update(mp[a[i]],m+1);
++vis[mp[a[i]]];
} LL sum3=0;
for(int i=1;i<=n;++i){
sum3+=(LL)C[i]*(LL)D[i];
sum3+=(LL)C[i]*(LL)B[i];
sum3+=(LL)A[i]*(LL)D[i];
sum3+=(LL)A[i]*(LL)B[i];
}
printf("%lld\n",sum1*sum2-sum3);
}
return 0;
}
线段树的超时代码:
# include<iostream>
# include<cstdio>
# include<map>
# include<vector>
# include<cstring>
# include<algorithm>
using namespace std;
# define LL long long
# define mid (l+(r-l)/2) const int N=50000; int a[N+5];
vector<int>v;
int A[N+5];//以i为右端点的顺序对个数
int B[N+5];//以i为右端点的逆序对个数
int C[N+5];//以i为左端点的顺序对个数
int D[N+5];//以i为左端点的逆序对个数
int tr[N*4+5];
map<int,int>mp; void pushUp(int rt)
{
tr[rt]=tr[rt<<1]+tr[rt<<1|1];
} void build(int rt,int l,int r)
{
tr[rt]=0;
if(l==r) return ;
build(rt<<1,l,mid);
build(rt<<1|1,mid+1,r);
} void update(int rt,int l,int r,int pos)
{
if(l==r)
++tr[rt];
else{
if(pos<=mid) update(rt<<1,l,mid,pos);
else update(rt<<1|1,mid+1,r,pos);
pushUp(rt);
}
} int query(int rt,int l,int r,int L,int R)
{
if(L<=l&&r<=R) return tr[rt];
int res=0;
if(L<=mid) res+=query(rt<<1,l,mid,L,R);
if(R>mid) res+=query(rt<<1|1,mid+1,r,L,R);
return res;
} int main()
{
int n;
while(~scanf("%d",&n))
{
mp.clear();
v.clear();
for(int i=1;i<=n;++i){
scanf("%d",a+i);
v.push_back(a[i]);
}
sort(v.begin(),v.end());
int len=unique(v.begin(),v.end())-v.begin();
for(int i=0;i<len;++i)
mp[v[i]]=i+1;
build(1,0,len+1);
LL sum1=0,sum2=0;
for(int i=1;i<=n;++i){
A[i]=query(1,0,len+1,0,mp[a[i]]-1);
B[i]=query(1,0,len+1,mp[a[i]]+1,len+1);
update(1,0,len+1,mp[a[i]]);
sum1+=A[i];
sum2+=B[i];
}
build(1,0,len+1);
for(int i=n;i>=1;--i){
C[i]=query(1,0,len+1,mp[a[i]]+1,len+1);
D[i]=query(1,0,len+1,0,mp[a[i]]-1);
update(1,0,len+1,mp[a[i]]);
}
LL sum3=0;
for(int i=1;i<=n;++i){
sum3+=(LL)C[i]*(LL)D[i];
sum3+=(LL)C[i]*(LL)B[i];
sum3+=(LL)A[i]*(LL)D[i];
sum3+=(LL)A[i]*(LL)B[i];
}
printf("%lld\n",sum1*sum2-sum3);
}
return 0;
}
HDU-5792 World is Exploding(树状数组)的更多相关文章
- HDU 5792 World is Exploding 树状数组+枚举
题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=5792 World is Exploding Time Limit: 2000/1000 MS (Ja ...
- hdu 5792 World is Exploding 树状数组
World is Exploding 题目连接: http://acm.hdu.edu.cn/showproblem.php?pid=5792 Description Given a sequence ...
- hdu 5792 World is Exploding 树状数组+离散化+容斥
World is Exploding Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Other ...
- HDU 5862 Counting Intersections(离散化+树状数组)
HDU 5862 Counting Intersections(离散化+树状数组) 题目链接http://acm.split.hdu.edu.cn/showproblem.php?pid=5862 D ...
- hdu 5517 Triple(二维树状数组)
Triple Time Limit: 12000/6000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)Total Sub ...
- 2016 Multi-University Training Contest 5 1012 World is Exploding 树状数组+离线化
http://acm.hdu.edu.cn/showproblem.php?pid=5792 1012 World is Exploding 题意:选四个数,满足a<b and A[a]< ...
- HDU 1394 Minimum Inversion Number ( 树状数组求逆序数 )
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1394 Minimum Inversion Number ...
- HDU 5862 Counting Intersections (树状数组)
Counting Intersections 题目链接: http://acm.split.hdu.edu.cn/showproblem.php?pid=5862 Description Given ...
- hdu 5592 ZYB's Game 树状数组
ZYB's Game Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://acm.hdu.edu.cn/showproblem.php?pid=55 ...
- HDU 1394 Minimum Inversion Number (树状数组求逆序对)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1394 题目让你求一个数组,这个数组可以不断把最前面的元素移到最后,让你求其中某个数组中的逆序对最小是多 ...
随机推荐
- 蓝桥杯 ALGO-4 结点选择 (树形动态规划)
问题描述 有一棵 n 个节点的树,树上每个节点都有一个正整数权值.如果一个点被选择了,那么在树上和它相邻的点都不能被选择.求选出的点的权值和最大是多少? 输入格式 第一行包含一个整数 n . 接下来的 ...
- Hibernate对象映射类型
Hibernate understands both the Java and JDBC representations of application data. The ability to rea ...
- Object Oriented Programming python
Object Oriented Programming python new concepts of the object oriented programming : class encapsula ...
- Asp.Net中Ajax实现登陆判断
Default.aspx: <head runat="server"> <title>无标题页</title> <script type= ...
- php大力力 [029节] 做PHP项目如何下载js文件:使用腾讯浏览器把网上案例页面存储到本地
php大力力 [029节] 做PHP项目如何下载js文件:使用腾讯浏览器把网上案例页面存储到本地 yeah,搞定啦 php大力力 [029节] 做PHP项目如何下载js文件:使用腾讯浏览器把网上案例页 ...
- php大力力 [005节] php大力力简单计算器001
2015-08-22 php大力力005. php大力力简单计算器001: 上网看视频,看了半天,敲击代码,如下: <html> <head> <title>简单计 ...
- 安装VMware Tools找不到内核头文件
http://blog.csdn.net/bobbat/article/details/38568885 安装VMware Tools,解决无法找到kernel header path的问题 安装 V ...
- CODEVS1380 没有上司的舞会 (树形DP)
f[i,0] 表示 第i个人不参加舞会 f[i,1] 表示 第i个人参加舞会 f[i,1]=sigma(f[j,0])+v[i] j 为 i 的孩子 f[i,1]=sigma(max(f[j,0] ...
- Unity3d各平台资源路径文件夹
之前一直是PC项目,公司终于考虑移动平台了,但是试验了几把,感觉移动平台资源管理路径还是有很多隐藏的注意事项. 比如在PC上可以做到随便读写,但是在移动平台就涉及到权限问题. 看到小伙伴的总结,还是要 ...
- mysql联合索引详解
联合索引又叫复合索引.对于复合索引:Mysql从左到右的使用索引中的字段,一个查询可以只使用索 引中的一部份,但只能是最左侧部分.例如索引是key index (a,b,c). 可以支持a | a,b ...