Contest
时间限制:1秒 空间限制:131072K

题目描述

n支队伍一共参加了三场比赛。
一支队伍x认为自己比另一支队伍y强当且仅当x在至少一场比赛中比y的排名高。
求有多少组(x,y),使得x自己觉得比y强,y自己也觉得比x强。
(x, y), (y, x)算一组。
 

输入描述:

第一行一个整数n,表示队伍数; 接下来n行,每行三个整数a[i], b[i], c[i],分别表示i在第一场、第二场和第三场比赛中的名次;n 最大不超过200000

输出描述:

输出一个整数表示满足条件的(x,y)数;64bit请用lld
示例1

输入

4
1 3 1
2 2 4
4 1 2
3 4 3

输出

5

求三维偏序对,在二维的基础上可以进行,这个可以BIT+CDQ啊,但是超时妥妥的

超时代码

#include<stdio.h>
#include<bits/stdc++.h>
using namespace std;
#define maxn 200002
typedef long long ll;
int T,n,tot,id[maxn];
inline char _getchar()
{
static const int BUFSIZE = ;
static char buf[BUFSIZE] ;
static char *psta=buf, *pend=buf ;
if (psta >= pend){
psta = buf;
pend = buf + fread(buf,,BUFSIZE,stdin);
if (psta >= pend)
return -;
}
return *psta++;
}
inline void read(int &x)
{
char c;
int ans=;
for(c=getchar(); c<''||c>''; c=getchar());
while(c>=''&&c<='')
ans=(ans<<)+(ans<<)+(c-''),c=getchar();
x=ans;
}
struct node
{
int x,y,z,cnt,ans,id;
} p[maxn];
int cmpx(node a,node b)
{
if(a.x!=b.x)return a.x<b.x;
if(a.y!=b.y)return a.y<b.y;
return a.z<b.z;
}
int cmpy(node a,node b)
{
if(a.y!=b.y)return a.y<b.y;
return a.z<b.z;
}
int cmpid(node a,node b)
{
return a.id<b.id;
}
struct BIT
{
#define lowbit(x) (x&(-x))
int b[maxn];
void init()
{
memset(b,,sizeof(b));
}
void update(int x,int v)
{
while(x<=maxn)
{
b[x]+=v;
x+=lowbit(x);
}
}
int query(int x)
{
int ans=;
while(x)
{
ans+=b[x];
x-=lowbit(x);
}
return ans;
}
void clear(int x)
{
while(x<=maxn)
{
b[x]=;
x+=lowbit(x);
}
}
} bit;
void CDQ(int l,int r)
{
if(l==r)
{
p[l].ans+=p[l].cnt-;
return ;
}
int mid=(l+r)>>;
CDQ(l,mid);
CDQ(mid+,r);
sort(p+l,p+mid+,cmpy);
sort(p+mid+,p+r+,cmpy);
int j=l;
for(int i=mid+; i<=r; i++)
{
for(; j<=mid&&p[j].y<=p[i].y; j++)
bit.update(p[j].z,p[j].cnt);
p[i].ans+=bit.query(p[i].z);
}
for(int i=l; i<j; i++)bit.update(p[i].z,-p[i].cnt);
sort(p+l,p+r+,cmpy);
}
int main()
{
read(n);
bit.init();
tot=;
for(int i=; i<=n; i++)
{
read(p[i].x);
read(p[i].y);
read(p[i].z);
p[i].ans=,p[i].id=i;
}
sort(p+,p+n+,cmpx);
for(int i=; i<=n; i++)
{
if(i!=&&p[i-].x==p[i].x&&p[i-].y==p[i].y&&p[i-].z==p[i].z)
p[tot].cnt++;
else p[++tot]=p[i],p[tot].cnt=;
id[p[i].id]=tot;
p[tot].id=tot;
}
CDQ(,tot);
sort(p+,p+tot+,cmpid);
ll ans=1ll*n*(n-)/;
for(int i=; i<=n; i++)
ans-=p[id[i]].ans;
printf("%lld",ans);
return ;
}

为什么要CDQ啊,严重拖慢了时间啊,但是扔到CF上是可以AC的,谁让CF的机器跑的快呢

附上AC代码

#include<bits/stdc++.h>
using namespace std;
#define LL long long
const int maxn=;
const int K = ;
struct Point
{
int x,y,z;
bool operator < (const Point& rhs) const {return z<rhs.z;}
}p[maxn];
int n,c[maxn/K+][maxn/K+];
inline int lowbit(int x) {return x&(-x);}
int query2D(int qx,int qy)
{
int ret=;
while(qy)
{
ret+=c[qx][qy];
qy-=lowbit(qy);
}
return ret;
}
int query1D(int qx,int qy)
{
int ret=;
while(qx)
{
ret+=query2D(qx,qy);
qx-=lowbit(qx);
}
return ret;
}
void modify2D(int px,int py,int v)
{
while(py<=(n/K))
{
c[px][py]+=v;
py+=lowbit(py);
}
}
void modify1D(int px,int py,int v)
{
while(px<=(n/K))
{
modify2D(px,py,v);
px+=lowbit(px);
}
}
int X[maxn],Y[maxn];
int main()
{
scanf("%d",&n);
for(int i=;i<n;i++) scanf("%d%d%d",&p[i].x,&p[i].y,&p[i].z);
for(int i=;i<=n;i++) X[i]=n+,Y[i]=n+;
sort(p,p+n);
LL ans=(LL)n*(n-)/2LL;
for(int i=;i<n;i++)
{
ans-=query1D(p[i].x/K,p[i].y/K);
for(int j=p[i].x;j>p[i].x-(p[i].x%K);j--) if(Y[j]<p[i].y) ans--;
for(int j=p[i].y;j>p[i].y-(p[i].y%K);j--) if(X[j]<=p[i].x-(p[i].x%K)) ans--;
X[p[i].y]=p[i].x,Y[p[i].x]=p[i].y;
modify1D((p[i].x+K-)/K,(p[i].y+K-)/K,);
}
printf("%lld\n",ans);
}

Wannafly模拟赛2的更多相关文章

  1. Wannafly模拟赛5 A Split 暴力 二分加速

    你有一个大小为

  2. Wannafly模拟赛 A.矩阵(二分答案+hash)

    矩阵 时间限制:1秒 空间限制:131072K 题目描述 给出一个n * m的矩阵.让你从中发现一个最大的正方形.使得这样子的正方形在矩阵中出现了至少两次.输出最大正方形的边长. 输入描述: 第一行两 ...

  3. Wannafly模拟赛

    题目描述 给出一个n * m的矩阵.让你从中发现一个最大的正方形.使得这样子的正方形在矩阵中出现了至少两次.输出最大正方形的边长. 输入描述: 第一行两个整数n, m代表矩阵的长和宽: 接下来n行,每 ...

  4. Wannafly模拟赛3

    喜欢想飞啊,毕竟这里的题目还都算挺好的中文题,不过做英文题久了还是感觉英文题更好理解一点 反蝴蝶效应 时间限制:1秒 空间限制:65536K 题目描述 一只南美洲亚马孙河流域热带雨林中的蝴蝶,偶尔扇动 ...

  5. 牛客网Wannafly模拟赛

    A矩阵 时间限制:1秒 空间限制:131072K 题目描述 给出一个n * m的矩阵.让你从中发现一个最大的正方形.使得这样子的正方形在矩阵中出现了至少两次.输出最大正方形的边长. 输入描述: 第一行 ...

  6. Wannafly模拟赛2 B river(拉格朗日乘数法)

    题目 https://www.nowcoder.com/acm/contest/4/B题意 有n条南北流向的河并列排着,水流速度是v,现在你需要从西岸游到东岸,总共T个时间,你的游泳速度是u,问东岸的 ...

  7. Wannafly模拟赛2 C alliances(dfs序+二分)

    题目 https://www.nowcoder.com/acm/contest/4/C 题意 由n个点组成一个树,有m个帮派,每个帮派由一些个点组成,这些点以及它们两两路径上的所有点都属于该帮派的管辖 ...

  8. Wannafly交流赛1_B_硬币【数学】

    Wannafly交流赛1_B_硬币[数学] 链接:https://www.nowcoder.com/acm/contest/69/B 来源:牛客网 题目描述 蜥蜴的生日快到了,就在这个月底! 今年,蜥 ...

  9. NOIP模拟赛20161022

    NOIP模拟赛2016-10-22 题目名 东风谷早苗 西行寺幽幽子 琪露诺 上白泽慧音 源文件 robot.cpp/c/pas spring.cpp/c/pas iceroad.cpp/c/pas ...

随机推荐

  1. 微软Coco Blockchain Framework:一键解决企业级区块链三大难题

    近年来,异军突起的“区块链”受到全行业的广泛关注,众多企业级用户在积极拥抱新技术的过程中却面临三大难题:性能.隐私和组织管理.如果不能很好地解决这些“顽固分子”,区块链技术就相对局限,很难发挥出应有的 ...

  2. 【Python图像特征的音乐序列生成】关于数据库到底在哪里下载

    毕竟原网站一个是14年前的一个是16年前的…… 1,http://ifdo.ca/~seymour/nottingham/nottingham.html 这个网站可以下载zip包. 2,https:/ ...

  3. 三维GIS-室内寻径功能实现

    期末,要交一个大作业,正巧之前跑国图借书的时候,晕头转向的,国图内居然没有导航!!!借这个机会做一个室内导航的demo,只是半成品,还需要加入室内定位,匹配一下坐标才能在实际中使用. demo:利用蜂 ...

  4. VC 对话框设置背景颜色和图片

    改变背景颜色,有两种方法: 1.在app的初始化函数中调用:void SetDialogBkColor( COLORREF clrCtlBk = RGB(192, 192, 192), COLORRE ...

  5. UVA439 knightMoves (A*启发搜索)

    第一个A*,纪念下. A*要保证最短路一定要估价函数小于等于实际值,越接近越好 估价函数取Manhattan距离除以二. //Rey #include<cstdio> #include&l ...

  6. [torch] torch.contiguous

    torch.contiguous 作用 连续存储,因为view的操作要求的是连续的内容. 详细 考虑下面的操作,transpose操作只是改变了stride,而实际数组存储的内容并没有得到任何改变,即 ...

  7. 使用js将div高度设置为100%

      在开发的工程中使用到了一些开源的bootstrap模板进行开发,在遇到一些需要替换的内容部分部分时,经常出现高度设置100%无法生效的问题,这里来用js强行设置一下.   思路:js监听窗口的缩放 ...

  8. 扫雷游戏 NOIP(入门)

    题目描述: 扫雷游戏是一款十分经典的单机小游戏.它的精髓在于,通过已翻开格子所提示的周围格地雷数,来判断未翻开格子里是否是地雷. 现在给出n行m列的雷区中的地雷分布,要求计算出每个非地雷格的周围格地雷 ...

  9. http 基础与通讯原理

    目录 http 基础与通讯原理 Internet 与中国 1990年10月 注册CN顶级域名 1993年3月2日 接入第一根专线 1994年4月20日 实现与互联网的全功能连接 1994年5月21日 ...

  10. Python数据分析实战视频教程【小蚊子数据分析实战课程】

    点击了解更多Python课程>>> Python数据分析实战视频教程[小蚊子数据分析实战课程] [课程概述] Python数据分析实战' 适用人群:适合需提升竞争力.提升工作效率.喜 ...