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 题目让你求一个数组,这个数组可以不断把最前面的元素移到最后,让你求其中某个数组中的逆序对最小是多 ...
随机推荐
- 最大子段和-Program A
最大子段和 Time Limit:1000MS Memory Limit:32768KB 64bit IO Format:%I64d & %I64u Description G ...
- Program C--二分
My birthday is coming up and traditionally I’m serving pie. Not just one pie, no, I have a number N ...
- const的全面理解
const关键字用来作甚?const是一个类型修饰符.常见的类型修饰符有哪些? short long unsigned signed static autoextern register 定义一个变量 ...
- Thread IsBcakgroud
C#中,Thread类有一个IsBackground 的属性.MSDN上对它的解释是:获取或设置一个值,该值指示某个线程是否为后台线程.个人感觉这样的解释等于没有解释. .Net中的线程,可以分为后台 ...
- SharePoint 2013 开发——发布SharePoint应用程序
博客地址:http://blog.csdn.net/FoxDave 前几篇我们介绍了开发.部署和调试SharePoint应用程序的基础,本篇介绍更实用的操作,当我们开发一个SharePoint应用 ...
- PHP Array 函数
PHP Array 简介 array 函数允许您对数组进行操作. PHP 支持单维和多维的数组.同时提供了用数据库查询结果来构造数组的函数. 安装 array 函数是 PHP 核心的组成部分.无需安装 ...
- LINQ基础 之 LINQ TO SQL (二)
配置LINQ TO SQL 首先添加一个Linq to sql文件,以.dbml结尾的文件.无法把表拖拽到.dbml文件中,提示“所选对象使用不支持的数据提供程序” 解决方案 在服务器资源管理器中右键 ...
- poj3126 筛素数+bfs
//Accepted 212 KB 16 ms //筛素数+bfs #include <cstdio> #include <cstring> #include <iost ...
- (转)Tomcat的目录结构
原文:http://hi.baidu.com/qinyougen/item/beeb506abb3e1d08a1cf0ffb Tomcat的目录结构 一.TOMCAT的目录结构 /bin:存放wind ...
- 2016- 1- 16 NSThread 的学习
一:NSThread的概念: 二:NSThread的使用: 1.创建一个Thread 1.1第一种方法: - (void)test1{ NSString *str = @"zhengli&q ...