若干道偏序问题(STL,分块)

找了4道题目

BZOJ陌上花开(权限题,提供洛谷链接)

Cogs2479偏序

Cogs2580偏序II

Cogs2639偏序++


作为一个正常人,肯定先看三维偏序

做法很多呀

首先,由于

智商不够数据结构来补 $ \ \ \ \ \ \ $--菊开

所以我们用最傻逼的数据结构来做这道题目

第一维:排序

第二维:树状数组

第三维:平衡树

于是乎,我们得到了一个复杂度为\(O(nlog^2n)\)的做法

并且常数巨大

这个做法到这里去看



第二种做法
CDQ分治
相信大家都会逆序对的求法
可以归并排序求逆序对
逆序对相当于是一个二维偏序
第一维就相当于它的位置
第二维就是值
每次归并排序的时候相当于分成了两部分
只有左侧的部分才能对右侧的部分产生贡献

这里同理,第一维排序

第二维类似于归并排序就行合并

同时一遍归并一边计算第三维

是不是想到了逆序对还可以用树状数组来求

所以这里一样的,

第三维用一个BIT来求就行了

这种做法看这里

恩,三维偏序没了


现在看Cogs的【偏序】

四维偏序模板题

可以顺着三维偏序来思考

既然多了一维

我就多套个CDQ分治呀

那不就是CDQ分治套树套树

或者CDQ分治套CDQ分治加树状数组

复杂度\(O(nlog^3n)\)

好像也可以诶。。。


好好好

我们继续

Cogs【偏序II】

不就是个五维偏序吗

CDQ套CDQ套CDQ

或者CDQ套CDQ套树套树

没毛病,时间复杂度\(O(nlog^4n)\)

可以做可以做,没问题的


那,我们继续

Cogs【偏序++】

七维偏序呀

CDQ套CDQ套CDQ套CDQ套CDQ

或者CDQ套CDQ套CDQ套CDQ套树套树

时间复杂度\(O(nlog^6n)\)

绝对跑不过的

而且给定的维数是K

你还得分类讨论,看着就是不可做的题目呀。。。

怎么办怎么办。。。


首先,请允许我orz大佬 ——中国翅王

这是他的PPT


FHR的做法的复杂度:\(O(n\sqrt n)\)

首先,我们知道

一个向量的偏序数量

等于它每一维的偏序数量的交集的大小

所以,按照每一维排序之后

我们就可以求出比每一维小的集合

然后求一个交集就行了

但是,暴力求,,,不现实吧。

而且交集也不好求呀

所以,来搞个\(bitset\)

STL大法好

但是,这样不就是\(O(n^2)\)的暴力吗???

没错,我们再来一发——分块大法好

分块之后就可以把一个\(n\)变成\(\sqrt n\)

于是复杂度就降到了\(O(n \sqrt n)\)

美滋滋

献上丑陋的代码(【偏序++】)

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<bitset>
#include<queue>
using namespace std;
#define MAX 45000
inline int read()
{
int x=0,t=1;char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=-1,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return x*t;
}
int n,blk;
bitset<45000> bs[15][250];
pair<int,int> val[15][MAX];
int bl[MAX],bll[MAX],blr[MAX];
int K,num;
long long Ans;
int f[15][MAX];
int find(int k,int x)
{
int l=1,r=n,ans=0;
while(l<=r)
{
int mid=(l+r)>>1;
if(val[k][mid].first<=x)ans=x,l=mid+1;
else r=mid-1;
}
return ans;
}
inline bitset<MAX> getbst(int p,int x)
{
bitset<MAX> ans;ans.reset();
int pp=find(p,x);
if(pp<=0)return ans;
int pre=(pp-1)/blk;
int st=pre*blk+1;
ans=bs[p][pre];
for(int i=st;i<=pp;++i)ans.set(val[p][i].second);
return ans;
}
void Solve()
{
bitset<MAX> ans;ans.reset();
for(int i=1;i<=n;++i)
{
ans.set();
for(int j=1;j<=K;++j)
ans&=getbst(j,f[j][i]);
Ans+=ans.count()-1;
}
}
int main()
{
freopen("partial_order_plus.in","r",stdin);
freopen("partial_order_plus.out","w",stdout);
n=read();blk=sqrt(n);K=read()+1;
for(int i=1;i<=n;++i)val[1][i]=make_pair(f[1][i]=i,i);
for(int j=2;j<=K;++j)
for(int i=1;i<=n;++i)
val[j][i]=make_pair(f[j][i]=read(),i);
for(int i=1;i<=K;++i)sort(&val[i][1],&val[i][n+1]);
for(int i=1;i<=n;++i)
bl[i]=(i-1)/blk+1;num=bl[n];
for(int i=1;i<=num;++i)
bll[i]=(i-1)*blk+1,blr[i]=i*blk;blr[num]=n;
for(int j=1;j<=K;++j)
for(int i=1;i<=num;++i)
{
bs[j][i]=bs[j][i-1];
for(int k=bll[i];k<=blr[i];++k)
bs[j][i][val[j][k].second]=1;
}
Solve();
printf("%lld\n",Ans);
return 0;
}

几道很Interesting的偏序问题的更多相关文章

  1. 两道很好的dp题目【4.29考试】

    A 问题描述: 对于一个排列,考虑相邻的两个元素,如果后面一个比前面一个大,表示这个位置是上升的,用I表示,反之这个位置是下降的,用D表示.如排列3,1,2,7,4,6,5可以表示为DIIDID. 现 ...

  2. 【学习笔记】使用 bitset 求解较高维偏序问题

    求解五维偏序 给定 \(n(\le 3\times 10^4)\) 个五元组,对于每个五元组 \((a_i, b_i, c_i, d_i, e_i)\),求存在多少个 \(1\le j\le n\) ...

  3. 高斯消元几道入门题总结POJ1222&&POJ1681&&POJ1830&&POJ2065&&POJ3185

    最近在搞高斯消元,反正这些题要么是我击败了它们,要么就是这些题把我给击败了.现在高斯消元专题部分还有很多题,先把几道很简单的入门题总结一下吧. 专题:http://acm.hust.edu.cn/vj ...

  4. 我也不知道什么是"莫比乌斯反演"和"杜教筛"

    我也不知道什么是"莫比乌斯反演"和"杜教筛" Part0 最近一直在搞这些东西 做了将近超过20道题目吧 也算是有感而发 写点东西记录一下自己的感受 如果您真的 ...

  5. Codeforces Round #415 (Div. 2) 翻车啦

    A. Straight «A» time limit per test 1 second memory limit per test 256 megabytes input standard inpu ...

  6. ZOJ Problem Set - 1006 Do the Untwist

    今天在ZOJ上做了道很简单的题目是关于加密解密问题的,此题的关键点就在于求余的逆运算: 比如假设都是正整数 A=(B-C)%D 则 B - C = D*n + A 其中 A < D 移项 B = ...

  7. wex5 实战 单页模式下的多页面数据同步

    在wex5官方教程中,关于多页模式与单页模式进行了对比.两者最大的区别在于: 1 web加载速度,单页模式快于多页模式 2  多页模式对加载机制进行了预加载,一次加载之后再次加载,就会加快. 但是,由 ...

  8. 【BZOJ-3306】树 线段树 + DFS序

    3306: 树 Time Limit: 10 Sec  Memory Limit: 256 MBSubmit: 792  Solved: 262[Submit][Status][Discuss] De ...

  9. hdu Remainder

    这道题是道很明显的bfs题.因为对数论没什么研究 ,所以这道题目里的两个关键点并不知道,看了别人的题解才知道 . 1.为避免取模后出现负数,采用:x%y=(x%y+y)%y 2.全部采用对m*k取模后 ...

随机推荐

  1. 2n皇后问题

    此题为蓝桥杯基础练习题. 问题描述 给定一个n*n的棋盘,棋盘中有一些位置不能放皇后.现在要向棋盘中放入n个黑皇后和n个白皇后,使任意的两个黑皇后都不在同一行.同一列或同一条对角线上,任意的两个白皇后 ...

  2. POJ 1023 The Fun Number System

    Description In a k bit 2's complement number, where the bits are indexed from 0 to k-1, the weight o ...

  3. Dynamics CRM 2015-Ribbon In Basic Home Tab

    前文中有说到关于Form上Ribbon的配置以及控制,而Ribbon Button还可以在其它地方的配置,今天就来说说在Basic Home Tab里面的配置,效果图如下: 具体配置Customiza ...

  4. Java经典编程题50道之四十三

    一个偶数总能表示为两个素数之和. public class Example43 {    public static void main(String[] args) {        f();   ...

  5. ch7复用类

    导出类的初始化是从基类开始向下扩展的,先初始化基类,再初始化由基类继承而来的类. 若类B需要类A中的一些甚至全部方法,但类B实际上不是并不是真正的类A,则可以通过代理的方式在B中实现所需要的A的方法, ...

  6. Yii2按需加载图片怎么做?

    按需加载图片应该用 jQuery LazyLoad 图片延迟加载按需加载文件夹应该用 Yii::import

  7. 关于UIButton嵌入到UIView点击无反应问题的解决方法和delegate的简单用法示例(转载)

    做项目封装UIView的时候碰到的问题,没想到有个哥们儿还写成博客,特此收藏! 问题是这样的,几个界面用到同一个自定义返回按钮,于是就想着把这个按钮单独封装起来,添加一个UIView类,在里面自定义U ...

  8. 硬件能力与智能AI-Zoomla!逐浪CMS2 x3.9.2正式发布

    北京时间2017年9月10日,领先的CMS网站内容管理系统与生产力软件研发厂商-Zoomla!逐浪CMS团队发布其年度重要产品:Zoomla!逐浪CMS2 x3.9.2,引领国内门户.移动.微商以及生 ...

  9. POJ - 2253 Frogger 单源最短路

    题意:给定n个点的坐标,问从第一个点到第二个点的最小跳跃范围.d(i)表示从第一个点到达第i个点的最小跳跃范围. AC代码 #include <cstdio> #include <c ...

  10. JavaScript的this和作用域

    本文主要讨论一下JS的作用域和this关键字.作用域,就是你的方法或者变量可访问的区域,是他们执行的上下文.如果你见过这样的代码: function someFunc() { var _this = ...