https://www.luogu.org/problemnew/show/P4125

题解

首先一开始所有的线段互不相交。

那么对于第二问来说,一定存在一种方法使得所有线段都朝着一个方向动。

比如说我们要让所有线段从上往下走。

那么上面的线段得向下面的线段连边。

这是一个\(DAG\),考虑怎么建出来。

我们可以先用扫描线,还是因为线段互不相交,所以在扫描线移动的过程中,当前所有线段的相对位置是不变的,所以我们可以把每条线段用斜截式表示,然后用set维护他们的关系,每次插入一条线段就和前驱后继连边就可以了。

然后因为一些恶心的原因得把坐标表示为\((x+1,y)\)。

然后再去考虑第一问,如果直接做的话,感觉需要用数据结构去维护\(DAG\)。

然后我们考虑把操作序列反过来,这样的话删线段就变成了加线段。然后考虑什么情况下加入是不合法的。

比如我们从下往上加线段,设这个线段的横坐标区间为\((l,r)\),那么如果不合法,那么当且仅当这个区间内有在它下面的线段,用上我们刚刚求的拓扑序就是这个线段的拓扑序比较小。

判断的话就线段树维护区间最小值就好了。

因为有四个方向,所以我们还得横纵各做一次扫描线,线段树维护最大值和最小值。

代码

#include<bits/stdc++.h>
#define N 200009
#define inf 2e9
using namespace std;
typedef long long ll;
queue<int>q;
int du[N],tot,head[N],n,xx[N],yy[N],ans[N],_xx[N],_yy[N],tag1[N],tag2[N],c[N],ans2,ans1[N];
ll nowx,b[N<<1];
inline ll rd(){
ll x=0;char c=getchar();bool f=0;
while(!isdigit(c)){if(c=='-')f=1;c=getchar();}
while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);c=getchar();}
return f?-x:x;
}
struct node{int a,b;}mu[N];
struct edge{int n,to;}e[N<<2];
inline void add(int u,int v){e[++tot].n=head[u];e[tot].to=v;head[u]=tot;du[v]++;}
struct point{double x,y;};
struct line{
point x,y;int id;
double k,b;
inline bool operator ==(const line &b)const{
return id==b.id;
}
inline bool operator <(const line &other)const{
return nowx*k+b<nowx*other.k+other.b;
}
}a[N];
struct segment_tree_t{
int mi[N<<2],ma[N<<2],trmi[N<<2],trma[N<<2];
segment_tree_t(){memset(mi,0x3f,sizeof(mi));memset(trmi,0x3f,sizeof(trmi));}
inline void pushdown(int cnt){
mi[cnt<<1]=min(mi[cnt<<1],mi[cnt]);
ma[cnt<<1]=max(ma[cnt<<1],ma[cnt]);
trmi[cnt<<1]=min(trmi[cnt<<1],mi[cnt]);
trma[cnt<<1]=max(trma[cnt<<1],ma[cnt]);
mi[cnt<<1|1]=min(mi[cnt<<1|1],mi[cnt]);
ma[cnt<<1|1]=max(ma[cnt<<1|1],ma[cnt]);
trmi[cnt<<1|1]=min(trmi[cnt<<1|1],mi[cnt]);
trma[cnt<<1|1]=max(trma[cnt<<1|1],ma[cnt]);
ma[cnt]=0;mi[cnt]=inf;
}
inline void ins(int cnt,int l,int r,int L,int R,int x){
if(l>=L&&r<=R){
ma[cnt]=max(ma[cnt],x);
mi[cnt]=min(mi[cnt],x);
trma[cnt]=max(trma[cnt],x);
trmi[cnt]=min(trmi[cnt],x);
return;
}
int mid=(l+r)>>1;
pushdown(cnt);
if(mid>=L)ins(cnt<<1,l,mid,L,R,x);
if(mid<R)ins(cnt<<1|1,mid+1,r,L,R,x);
trmi[cnt]=min(trmi[cnt<<1],trmi[cnt<<1|1]);
trma[cnt]=min(trma[cnt<<1],trma[cnt<<1|1]);
}
inline int querymin(int cnt,int l,int r,int L,int R){
if(l>=L&&r<=R)return trmi[cnt];
int mid=(l+r)>>1,ans=inf;
pushdown(cnt);
if(mid>=L)ans=min(ans,querymin(cnt<<1,l,mid,L,R));
if(mid<R)ans=min(ans,querymin(cnt<<1|1,mid+1,r,L,R));
return ans;
}
inline int querymax(int cnt,int l,int r,int L,int R){
if(l>=L&&r<=R)return trma[cnt];
int mid=(l+r)>>1,ans=0;
pushdown(cnt);
if(mid>=L)ans=max(ans,querymax(cnt<<1,l,mid,L,R));
if(mid<R)ans=max(ans,querymax(cnt<<1|1,mid+1,r,L,R));
return ans;
}
}X,Y;
set<line>s;
set<line>::iterator its;
vector<int>vec[N],anti_vec[N];
vector<int>::iterator it;
int main(){
n=rd();
for(int i=1;i<=n;++i){
a[i].x.x=rd();a[i].x.y=rd();a[i].y.x=rd();a[i].y.y=rd();
if(a[i].x.x>a[i].y.x)swap(a[i].x,a[i].y);
if(a[i].x.x!=a[i].y.x){
a[i].k=(a[i].y.y-a[i].x.y)/(a[i].y.x-a[i].x.x);
a[i].b=a[i].y.y-a[i].y.x*a[i].k;
}
else a[i].b=a[i].y.y;
b[++b[0]]=a[i].x.x+1;b[++b[0]]=a[i].y.x;
}
sort(b+1,b+b[0]+1);
b[0]=unique(b+1,b+b[0]+1)-b-1;
for(int i=1;i<=n;++i){
xx[i]=lower_bound(b+1,b+b[0]+1,a[i].x.x+1)-b;
yy[i]=lower_bound(b+1,b+b[0]+1,a[i].y.x)-b;
a[i].id=i;
vec[xx[i]].push_back(i);
anti_vec[yy[i]].push_back(i);
}
for(int i=1;i<=b[0];++i){
nowx=b[i];
for(it=vec[i].begin();it!=vec[i].end();++it){
line now=a[*it];
its=s.insert(now).first;
++its;
if(its!=s.end()){
add(now.id,its->id);
}
--its;
if(its!=s.begin()){
its--;
add(its->id,now.id);
}
}
for(it=anti_vec[i].begin();it!=anti_vec[i].end();++it){
line now=a[*it];
s.erase(now);
}
vec[i].clear();anti_vec[i].clear();
}
for(int i=1;i<=n;++i)if(!du[i])q.push(i);
while(!q.empty()){
int u=q.front();q.pop();ans[++ans[0]]=u;tag1[u]=ans[0];
for(int i=head[u];i;i=e[i].n){
int v=e[i].to;
if(!--du[v])q.push(v);
}
}
memset(head,0,sizeof(head));tot=0;
for(int i=1;i<=n;++i){
if(a[i].x.y>a[i].y.y)swap(a[i].x,a[i].y);
if(a[i].x.y!=a[i].y.y){
a[i].k=(a[i].y.x-a[i].x.x)/(a[i].y.y-a[i].x.y);
a[i].b=a[i].y.x-a[i].y.y*a[i].k;
}
else a[i].k=0,a[i].b=a[i].y.x;
c[++c[0]]=a[i].x.y+1;;c[++c[0]]=a[i].y.y;
}
sort(c+1,c+c[0]+1);
c[0]=unique(c+1,c+c[0]+1)-c-1;
for(int i=1;i<=n;++i){
_xx[i]=lower_bound(c+1,c+c[0]+1,a[i].x.y+1)-c;
_yy[i]=lower_bound(c+1,c+c[0]+1,a[i].y.y)-c;
a[i].id=i;
vec[_xx[i]].push_back(i);
anti_vec[_yy[i]].push_back(i);
}
for(int i=1;i<=c[0];++i){
nowx=c[i];
for(it=vec[i].begin();it!=vec[i].end();++it){
line now=a[*it];
its=s.insert(now).first;
++its;
if(its!=s.end()){
add(now.id,its->id);
}
--its;
if(its!=s.begin()){
its--;
add(its->id,now.id);
}
}
for(it=anti_vec[i].begin();it!=anti_vec[i].end();++it){
line now=a[*it];
s.erase(now);
}
}
for(int i=1;i<=n;++i)if(!du[i])q.push(i);
while(!q.empty()){
int u=q.front();q.pop();ans1[++ans1[0]]=u;tag2[u]=ans1[0];
for(int i=head[u];i;i=e[i].n){
int v=e[i].to;
if(!--du[v])q.push(v);
}
}
for(int i=1;i<=n;++i){
mu[i].a=rd();mu[i].b=rd();
}
for(int i=n;i>=1;--i){
int id=mu[i].a;
if(mu[i].b==0){
if(Y.querymin(1,1,c[0],_xx[id],_yy[id])<tag2[id]){
ans2=i;
}
}
if(mu[i].b==1){
if(X.querymax(1,1,b[0],xx[id],yy[id])>tag1[id]){
ans2=i;
}
}
if(mu[i].b==2){
if(Y.querymax(1,1,c[0],_xx[id],_yy[id])>tag2[id]){
ans2=i;
}
}
if(mu[i].b==3){
if(X.querymin(1,1,b[0],xx[id],yy[id])<tag1[id]){
ans2=i;
}
}
X.ins(1,1,b[0],xx[id],yy[id],tag1[id]);
Y.ins(1,1,c[0],_xx[id],_yy[id],tag2[id]);
}
printf("%d\n",ans2);
for(int i=1;i<=n;++i)printf("%d %d\n",ans[i],3);
return 0;
}

[WC2012]记忆中的水杉树的更多相关文章

  1. 洛谷 P4125 [WC2012]记忆中的水杉树【扫描线+set+线段树】

    我没有找到能在bzojAC的代码--当然我也WA了--但是我在洛谷过了,那就假装过了吧 minmax线段树一开始写的只能用min更新min,max更新max,实际上是可以互相更新的-- 首先看第二问, ...

  2. 某寺庙,有小和尚、老和尚若干。有一水缸,由小和尚用水桶从井中提水入缸,老和尚用水桶从缸里取水饮用。水缸可容10桶水,水取自同一井中。水井径窄,每次只能容一个水桶取水。水桶总数为3个。每次入、取缸水仅为1桶,且不可以同时进行。试用P、V操作给出小和尚、老和尚动作的算法描述。

    寺庙和尚打水 设信号量mutex_gang, mutex_jing, gang_empty, gang_full, count分别表示使用缸互斥, 使用井互斥, 缸空, 缸满, 水桶总个数 semap ...

  3. CityEngine中动态水的实现

    地址:http://pan.baidu.com/share/link?shareid=3871210059&uk=3492170216 密码:am5b 在今年Esri全球用户大会和Esri中国 ...

  4. 聊聊Android5.0中的水波纹效果

    水波纹效果已经不是什么稀罕的东西了,用过5.0新控件的小伙伴都知道这个效果,可是如果使用一个TextView或者Button或者其它普通控件的话,你是否知道如何给它设置水波纹效果呢?OK,我们今天就来 ...

  5. 记忆中的像素块褪色了吗?用开源的体素编辑器重新做个 3D 的吧!

    本文适合对图形表现有兴趣的美术或者开发人员 本文作者:HelloGitHub-Joey 早期的的显示设备像素颗粒较大,使得显示内容的颗粒感严重,像是由一堆方块组成的.比较好的例子就是 GBA 上的游戏 ...

  6. 借助Algorithmia网站API:用AI给黑白照片上色,复现记忆中的旧时光

    先看DEMOhttps://demos.algorithmia.com/colorize-photos/ 了解ColorfulImageColorizationhttps://algorithmia. ...

  7. BZOJ 4971: [Lydsy1708月赛]记忆中的背包

    神仙构造 分成x个1和一堆>=w-x的大物品 (x<=20 w>=50) 则拼成w的方案中有且仅有一个大物品 若最终序列中有x个1,有一个大物品为w-k,可以提供C(x,k)种方案 ...

  8. poj1163 【记忆化搜索·水】

    题意: 一个这样的三角形,他可以往下的左或者往下的右走.求一个在最后一行的最大. 思路: 额...就是搜一下..记录一下...肯定有重合的情况. code- //#include <bits/s ...

  9. [bzoj4971]记忆中的背包

    为了使得方案的形式较为单一,不妨强制物品体积为1或$\ge \lceil\frac{w}{2}\rceil$,那么假设最终有$x$个1且$\ge \lceil\frac{w}{2}\rceil$的物品 ...

随机推荐

  1. 从DT时代云栖大会聊聊恒生电子

    从IT到DT,除了HOMS和配资,本文结合互联网金融的背景,帮助读者对恒生电子600570了解更多. 马云在2015杭州·云栖大会Computing Conference发表致辞时多次强调DT(Dat ...

  2. Map m=new HashMap()

    Map<String,String> m=new HashMap<String,String>() 等于 HashMap<String,String> hashMa ...

  3. vue封装element中table组件

    后台系统,table被用的次数比较多,所以决定提出来作为组件 1.新建一个Table.vue文件 <!--region 封装的分页 table--> <template>  & ...

  4. Java ——日期时间 日期时间相关类 随机数 定义类属性时建议使用引用数据类型

    本节重点思维导图 Date对象创建 两个构造函数: Date() ----使用当前日期和时间来初始化对象 Date(long millisec) -----接收一个参数,该参数是从1970年1月1日起 ...

  5. python接口自动化:调试接口的代码(无token情况下)

    实现代码如下: #接口调试 import requests,time class api_demo1: def __init__(self,RequestWay,url,data): self.s=r ...

  6. JAVA调用shell脚本利用ansible修改多节点上的redis参数

    创建hosts文件 创建ansible-playbook执行时所用到的hosts文件,例如 /etc/redis/hosts 利用shell命令根据传入的host名和地址写入hosts文件: #set ...

  7. python字典-基础

    一.解释 像列表一样,“字典”是许多值的集合.但不像列表的下标,字典的索引可以 使用许多不同数据类型,不只是整数.字典的索引被称为“键”,键及其关联的值 称为“键-值”对. 二.列表创建方式 1. I ...

  8. C++学习笔记(五)--指针、NULL、引用

    1. C++中已经定义了NULL为0:#define NULL 0 指针p可以指向空值NULL即 p = NULL;表示该指针变量不指向任何变量. 注意:指针未初始化与指针为NULL不同, p = N ...

  9. 从SVN下检出项目内容【步骤】

    1.新创建一个新的工作环境,然后new--->other--->SVN 2.点击Next,然后进行检出项目的操作,如下图所示: 3.再点击Next,进行输入指定的url地址,从指定的url ...

  10. [Codeforces 1205B]Shortest Cycle(最小环)

    [Codeforces 1205B]Shortest Cycle(最小环) 题面 给出n个正整数\(a_i\),若\(a_i \& a_j \neq 0\),则连边\((i,j)\)(注意i- ...