[2018雅礼3-27]subset 分类讨论+三维偏序+二维偏序

首先有个显然的结论,对于一个集合 S S ,当ax,bx,cx均不为最大值时从 S S 中删去x,那么 |S|3 | S | ≤ 3 且一个 S S 唯一对应一个答案。
我们先用O(nlog2n)的三维偏序预处理出对于每个 x x ,有多少个y满足 ax>ay,bx>by,cx>cy a x > a y , b x > b y , c x > c y ,记作 dx d x
类似的,用 O(nlogn) O ( n log ⁡ n ) 的二维偏序预处理出有两维大于,另一维不限的个数,记作 dabx,dbcx,dacx d a b x , d b c x , d a c x
于是讨论 |S| | S |
1. |S|=1 | S | = 1 ,显然所有单独的列都为一种合法的 S S ,贡献n
2. |S|=2 | S | = 2 ,那么要从所有二元集合里去除某一列比另一列三维都大的情况,剩下的均为合法 S S ,那么就是(n2)i=1ndi
3. |S|=3 | S | = 3 ,先考虑某一列三维均为最大的情况,记作 A A ;再考虑某一列至少有两个最大值的情况,记作B,那么 |A|=ni=1(di2),|B|=ni=1((dabi2)+(dbci2)+(daci2)) | A | = ∑ i = 1 n ( d i 2 ) , | B | = ∑ i = 1 n ( ( d a b i 2 ) + ( d b c i 2 ) + ( d a c i 2 ) ) ,但是 A A 会在B中算重三次,所以最后贡献就是 (n3)|B|+2|A| ( n 3 ) − | B | + 2 | A |

代码:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 100010
#define ll long long
using namespace std;
int n,c[N];
ll ans;
int read()
{
    int x=0,f=1;char ch=getchar();
    for(;ch<'0'||ch>'9';ch=getchar()) if(ch=='-') f=-1;
    for(;ch>='0'&&ch<='9';ch=getchar()) x=x*10+ch-'0';
    return x*f;
}
struct node
{
    int x,y,z,d,dx,dy,dz;
}a[N],b[N];
void add(int x,int d){for(;x<=n;x+=(x&-x)) c[x]+=d;}
int qry(int x){int r=0;for(;x;x-=(x&-x)) r+=c[x];return r;}
bool cmpx(node p,node q)
{
    return p.x<q.x;
}
bool cmpy(node p,node q)
{
    return p.y<q.y;
}
ll G(int x)
{
    return (ll)x*(x-1)/2;
}
void solve(int L,int R,int l,int r)
{
    if(L>=R||l>=r) return ;
    int mid=(l+r>>1);
    for(int i=L;i<=R;i++)
        if(a[i].y<=mid) add(a[i].z,1);
        else a[i].d+=qry(a[i].z);   
    int top=L-1,tmp;
    for(int i=L;i<=R;i++)
        if(a[i].y<=mid) b[++top]=a[i],add(a[i].z,-1);
    tmp=top;
    for(int i=L;i<=R;i++)
        if(a[i].y>mid) b[++top]=a[i];
    for(int i=L;i<=R;i++)
        a[i]=b[i];
    solve(L,tmp,l,mid);
    solve(tmp+1,R,mid+1,r);             
}
int main()
{
    n=read();
    for(int i=1;i<=n;i++)
        a[i].x=read();
    for(int i=1;i<=n;i++)
        a[i].y=read();
    for(int i=1;i<=n;i++)
        a[i].z=read();
    sort(a+1,a+n+1,cmpx);   
    solve(1,n,1,n);
    sort(a+1,a+n+1,cmpx);
    memset(c,0,sizeof(c));
    for(int i=1;i<=n;i++)
        a[i].dz+=qry(a[i].y),add(a[i].y,1);
    memset(c,0,sizeof(c));
    for(int i=1;i<=n;i++)
        a[i].dy+=qry(a[i].z),add(a[i].z,1);
    memset(c,0,sizeof(c));  
    sort(a+1,a+n+1,cmpy);
    for(int i=1;i<=n;i++)
        a[i].dx+=qry(a[i].z),add(a[i].z,1); 
    ans=G(n)+n; 
    for(int i=1;i<=n;i++)
        ans-=a[i].d;
    ans+=G(n)*(n-2)/3;
    for(int i=1;i<=n;i++)
    {
        ll A=G(a[i].d),B=G(a[i].dx)+G(a[i].dy)+G(a[i].dz);
        ans-=(B-2*A);
    }
    printf("%lld",ans);
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
SQLAlchemy 是一个 SQL 工具包和对象关系映射(ORM)库,用于 Python 编程语言。它提供了一个高级的 SQL 工具和对象关系映射工具,允许开发者以 Python 类和对象的形式操作数据库,而无需编写大量的 SQL 语句。SQLAlchemy 建立在 DBAPI 之上,支持多种数据库后端,如 SQLite, MySQL, PostgreSQL 等。 SQLAlchemy 的核心功能: 对象关系映射(ORM): SQLAlchemy 允许开发者使用 Python 类来表示数据库表,使用类的实例表示表中的行。 开发者可以定义类之间的关系(如一对多、多对多),SQLAlchemy 会自动处理这些关系在数据库中的映射。 通过 ORM,开发者可以像操作 Python 对象一样操作数据库,这大大简化了数据库操作的复杂性。 表达式语言: SQLAlchemy 提供了一个丰富的 SQL 表达式语言,允许开发者以 Python 表达式的方式编写复杂的 SQL 查询。 表达式语言提供了对 SQL 语句的灵活控制,同时保持了代码的可读性和可维护性。 数据库引擎和连接池: SQLAlchemy 支持多种数据库后端,并且为每种后端提供了对应的数据库引擎。 它还提供了连接池管理功能,以优化数据库连接的创建、使用和释放。 会话管理: SQLAlchemy 使用会话(Session)来管理对象的持久化状态。 会话提供了一个工作单元(unit of work)和身份映射(identity map)的概念,使得对象的状态管理和查询更加高效。 事件系统: SQLAlchemy 提供了一个事件系统,允许开发者在 ORM 的各个生命周期阶段插入自定义的钩子函数。 这使得开发者可以在对象加载、修改、删除等操作时执行额外的逻辑。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值