记录编号 |
365910 |
评测结果 |
AAAAAAAAAA |
题目名称 |
[HZOI 2016] 偏序 |
最终得分 |
100 |
用户昵称 |
SiriusRen |
是否通过 |
通过 |
代码语言 |
C++ |
运行时间 |
9.378 s |
提交时间 |
2017-01-22 00:22:44 |
内存使用 |
150.22 MiB |
显示代码纯文本
//By SiriusRen
#include <bitset>
#include <cstdio>
#include <algorithm>
using namespace std;
int n;long long ans;
bitset<25000>b[50000],jy;
struct Node{int a,b,c,id;}node[50000];
bool cmpa(const Node &a,const Node &b){return a.a<b.a;}
bool cmpb(const Node &a,const Node &b){return a.b<b.b;}
bool cmpc(const Node &a,const Node &b){return a.c<b.c;}
inline int rd(){
int x=0;char p=getchar();
while(p<'0'||p>'9')p=getchar();
while(p>='0'&&p<='9')x=x*10+p-'0',p=getchar();
return x;
}
int main(){
freopen("partial_order.in","r",stdin);
freopen("partial_order.out","w",stdout);
scanf("%d",&n);
register int i;
for(i=0;i<n/2;i++)b[i]|=jy,jy.set(i);
for(i=n/2;i<n;i++)b[i]|=jy;
for(i=0;i<n;i++)node[i].a=rd();
for(i=0;i<n;i++)node[i].b=rd();
for(i=0;i<n;node[i].id=i,i++)node[i].c=rd();
jy.reset();sort(node,node+n,cmpa);
for(i=0;i<n;i++){b[node[i].id]&=jy;if(node[i].id<n/2)jy.set(node[i].id);}
jy.reset();sort(node,node+n,cmpb);
for(i=0;i<n;i++){b[node[i].id]&=jy;if(node[i].id<n/2)jy.set(node[i].id);}
jy.reset();sort(node,node+n,cmpc);
for(i=0;i<n;i++){b[node[i].id]&=jy;if(node[i].id<n/2)jy.set(node[i].id);}
for(i=0;i<n;i++)ans+=b[i].count();
jy.reset();
for(i=0;i<n;i++)b[i].reset();
for(i=n/2;i<n;i++)b[i]|=jy,jy.set(i-n/2);
jy.reset();sort(node,node+n,cmpa);
for(i=0;i<n;i++){b[node[i].id]&=jy;if(node[i].id>=n/2)jy.set(node[i].id-n/2);}
jy.reset();sort(node,node+n,cmpb);
for(i=0;i<n;i++){b[node[i].id]&=jy;if(node[i].id>=n/2)jy.set(node[i].id-n/2);}
jy.reset();sort(node,node+n,cmpc);
for(i=0;i<n;i++){b[node[i].id]&=jy;if(node[i].id>=n/2)jy.set(node[i].id-n/2);}
for(i=0;i<n;i++)ans+=b[i].count();
printf("%lld\n",ans);
}