洛谷 P4125 [WC2012]记忆中的水杉树【扫描线+set+线段树】
我没有找到能在bzojAC的代码……当然我也WA了……但是我在洛谷过了，那就假装过了吧
minmax线段树一开始写的只能用min更新min，max更新max，实际上是可以互相更新的……
首先看第二问，注意到因为没有相交，所以可以全都按某种顺序像同一个方向移动来完成游戏，这个顺序通过扫描线找到，用set维护经过当前x坐标的的线段的y坐标，具体用y=kx+b的形式维护，因为不相交，所以随着x的变化，线段的y点的大小关系不会改变；每次插入一个点的时候，找一下它的前驱后继，分别连边（这里的前驱后继就相当于向上和向下第一个撞上的线段），然后得到一张有向图，拓扑序就是游戏的操作顺序
然后看第一问，对横向和纵向都做一遍拓扑，记录点在拓扑序中的出现的时间，然后用两个线段树分别维护对于横向/纵向点的出现时间的最大最小值，然后倒着做，每次查询线段树里就相当于查询游戏中剩下的线段的出现时间，判一下合法即可
#include<iostream>
#include<cstdio>
#include<set>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int N=300005;
int n,p[N],q[N],d[N],r[2][N],hs[2][N],ht[2][N],c[N],wi,h[N],cnt,tp[N],ans;
struct dian
{
int x,y;
}a[N],b[N];
struct bian
{
int ne,to;
}e[N<<1];
struct xds
{
int l,r,mn,mx,lmn,lmx;
void add(int v)
{
if(v<mn)
mn=v;
if(v>mx)
mx=v;
if(v<lmn)
lmn=v;
if(v>lmx)
lmx=v;
}
};
struct qwe
{
int x,id,f;
double k,b;
qwe(int X=0,int ID=0,int F=0,double K=0,double B=0)
{
x=X,id=ID,f=F,k=K,b=B;
}
bool operator < (const qwe &z) const
{
return k*wi+b<z.k*wi+z.b;
}
bool operator == (const qwe &z) const
{
return id==z.id;
}
}f[N<<1];
bool cmp(const qwe &a,const qwe &b)
{
return a.x<b.x||(a.x==b.x&&a.f>b.f);
}
int read()
{
int r=0,f=1;
char p=getchar();
while(p>'9'||p<'0')
{
if(p=='-')
f=-1;
p=getchar();
}
while(p>='0'&&p<='9')
{
r=r*10+p-48;
p=getchar();
}
return r*f;
}
void add(int u,int v)
{//cerr<<u<<" "<<v<<endl;
cnt++;
e[cnt].ne=h[u];
e[cnt].to=v;
h[u]=cnt;
d[v]++;
}
void wk(int hs[],int ht[],int r[])
{
int tot=0,has=0;
for(int i=1;i<=n;i++)
c[++tot]=a[i].x+1,c[++tot]=b[i].x;
sort(c+1,c+1+tot);
for(int i=1;i<=tot;i++)
if(i==1||c[i]!=c[i-1])
c[++has]=c[i];
tot=0;
for(int i=1;i<=n;i++)
{
hs[i]=lower_bound(c+1,c+1+has,a[i].x+1)-c,ht[i]=lower_bound(c+1,c+1+has,b[i].x)-c;
f[++tot]=qwe(hs[i],i,1,(double)(b[i].y-a[i].y)/(double)(b[i].x-a[i].x),a[i].y-a[i].x*(double)(b[i].y-a[i].y)/(double)(b[i].x-a[i].x));
f[++tot]=qwe(ht[i],i,-1,(double)(b[i].y-a[i].y)/(double)(b[i].x-a[i].x),a[i].y-a[i].x*(double)(b[i].y-a[i].y)/(double)(b[i].x-a[i].x));
}
sort(f+1,f+1+tot,cmp);
memset(d,0,sizeof(d));
memset(h,0,sizeof(h));
cnt=0;
set<qwe>s;
set<qwe>::iterator it;
for(int i=1,w=1;i<=has;i++)
{
wi=c[i];
while(w<=tot&&f[w].x==i)
{
if(f[w].f==1)
{
it=s.lower_bound(f[w]);
if(it!=s.end())
add(f[w].id,it->id);
if(it!=s.begin())
add((--it)->id,f[w].id);
s.insert(f[w]);
}
else
s.erase(f[w]);
w++;
}
}
tot=0;
queue<int>q;
for(int i=1;i<=n;i++)
if(!d[i])
q.push(i);
while(!q.empty())
{
int u=q.front();
r[u]=++tot;
q.pop();
for(int i=h[u];i;i=e[i].ne)
if(!(--d[e[i].to]))
q.push(e[i].to);
}
}
struct clc
{
xds t[N<<2];
void build(int ro,int l,int r)
{
t[ro].l=l,t[ro].r=r,t[ro].mx=t[ro].lmx=-1e9,t[ro].mn=t[ro].lmn=1e9;
if(l==r)
return;
int mid=(l+r)>>1;
build(ro<<1,l,mid);
build(ro<<1|1,mid+1,r);
}
void pd(int ro)
{
if(t[ro].lmx!=-1e9)
{
t[ro<<1].add(t[ro].lmx);
t[ro<<1|1].add(t[ro].lmx);
t[ro].lmx=-1e9;
}
if(t[ro].lmn!=1e9)
{
t[ro<<1].add(t[ro].lmn);
t[ro<<1|1].add(t[ro].lmn);
t[ro].lmn=1e9;
}
}
void update(int ro,int l,int r,int v)
{
if(t[ro].l==l&&t[ro].r==r)
{
t[ro].add(v);
return;
}
pd(ro);
int mid=(t[ro].l+t[ro].r)>>1;
if(r<=mid)
update(ro<<1,l,r,v);
else if(l>mid)
update(ro<<1|1,l,r,v);
else
update(ro<<1,l,mid,v),update(ro<<1|1,mid+1,r,v);
t[ro].mn=min(t[ro<<1].mn,t[ro<<1|1].mn);
t[ro].mx=max(t[ro<<1].mx,t[ro<<1|1].mx);
}
pair<int,int>ques(int ro,int l,int r)
{
if(t[ro].l==l&&t[ro].r==r)
return make_pair(t[ro].mn,t[ro].mx);
pd(ro);
int mid=(t[ro].l+t[ro].r)>>1;
if(r<=mid)
return ques(ro<<1,l,r);
else if(l>mid)
return ques(ro<<1|1,l,r);
else
{
pair<int,int>x=ques(ro<<1,l,mid),y=ques(ro<<1|1,mid+1,r);
return make_pair(min(x.first,y.first),max(x.second,y.second));
}
}
}t0,t1;
int main()
{
n=read();
for(int i=1;i<=n;i++)
{
a[i].x=read(),a[i].y=read(),b[i].x=read(),b[i].y=read();
if(a[i].x>b[i].x)
swap(a[i],b[i]);
}
wk(hs[0],ht[0],r[0]);//cerr<<"OK"<<endl;
for(int i=1;i<=n;i++)
{
swap(a[i].x,a[i].y),swap(b[i].x,b[i].y);
if(a[i].x>b[i].x)
swap(a[i],b[i]);
}
wk(hs[1],ht[1],r[1]);
// for(int i=1;i<=n;i++)
// cerr<<r[0][i]<<" ";cerr<<endl;
// for(int i=1;i<=n;i++)
// cerr<<r[1][i]<<" ";cerr<<endl;
t0.build(1,1,2*n);
t1.build(1,1,2*n);
for(int i=1;i<=n;i++)
p[i]=read(),q[i]=read();
for(int i=n;i>=1;i--)
{
if(q[i]==1||q[i]==3)
{
pair<int,int>nw=t0.ques(1,hs[0][p[i]],ht[0][p[i]]);
if((q[i]==1&&nw.second>r[0][p[i]])||(q[i]==3&&nw.first<r[0][p[i]]))
ans=i;
}
else
{
pair<int,int>nw=t1.ques(1,hs[1][p[i]],ht[1][p[i]]);//cerr<<nw.first<<" "<<nw.second<<endl;
if((q[i]==2&&nw.second>r[1][p[i]])||(q[i]==0&&nw.first<r[1][p[i]]))
ans=i;
}
t0.update(1,hs[0][p[i]],ht[0][p[i]],r[0][p[i]]);
t1.update(1,hs[1][p[i]],ht[1][p[i]],r[1][p[i]]);
}
printf("%d\n",ans);
for(int i=1;i<=n;i++)
tp[r[0][i]]=i;
for(int i=1;i<=n;i++)
printf("%d 3\n",tp[i]);
return 0;
}
