题目描述
你猜
http滑稽://172.1滑稽6.0.132/sen滑稽ior/#contest滑稽/show/237滑稽2/8
30%
暴力?
反正我写的暴力1分都拿不了
60%
不知道
80%
不知道
100%
其实可以直接连边+Tarjan,但是慢在连边
O(n2)
O
(
n
2
)
于是考虑用线段树优化。
因为Tarjan找的是强连通分量,也就是说a–>b和a–>c–>b是等价的
而题目要求的是abc中的两项,所以可以分开来做三遍
按照一维排序,另一维丢到权值线段树里
把线段树的每个区间都当作一个点,插入一个点时把所有包含它的区间跟它连边,再把它跟所有它包含的区间连边
显然两者的个数都是
logn
log
n
的
那么一共就连了
nlogn
n
log
n
条边
但是因为修改有时效性,前面的点不能包括后面的点
所以要用主席树
最后Tarjan就行了
(还有可以先离散化一下在搞,这样会方便&快一些)
我不会告诉你我没开人工栈就A了
尤其是本机跑不过但OJ过了
code
#include <algorithm>
#include <iostream>
#include <cstdlib>
#include <cstdio>
#define fo(a,b,c) for (a=b; a<=c; a++)
#define fd(a,b,c) for (a=b; a>=c; a--)
#define min(a,b) (a<b?a:b)
using namespace std;
int tr[10000001][2];
int a[20000001][2];
int ls[10000001];
struct AA{int a,b,c,id;} A[100001];
int ans[10000001];
int Ans[100001];
int dfn[10000001];
int low[10000001];
int d[10000001];
bool bz[10000001];
int Id[100001];
int n,m,i,j,k,l,len;
void swap(int &x,int &y) {int z=x;x=y;y=z;}
bool cmp(AA a,AA b) {return a.a<b.a;}
bool Cmp(int a,int b) {return a<b;}
void _New(int x,int y)
{
m++;
a[m][0]=y;
a[m][1]=ls[x];
ls[x]=m;
}
void New(int t,int x)
{
len++;
_New(len,tr[t][x]);
tr[len][0]=tr[tr[t][x]][0];
tr[len][1]=tr[tr[t][x]][1];
tr[t][x]=len;
}
void scan()
{
int i;
scanf("%d",&n);
fo(i,1,n)
scanf("%d%d%d",&A[i].a,&A[i].b,&A[i].c),A[i].id=i;
len=n;
}
void init()
{
int i,j;
fo(i,1,3)
{
sort(A+1,A+n+1,cmp);
fo(j,1,n) A[j].a=j,swap(A[j].a,A[j].b),swap(A[j].b,A[j].c);
}
}
void find(int I,int t,int l,int r,int y)
{
int mid=(l+r)/2;
if (r<=y)
{
if (l<r)
_New(I,t);
else
_New(I,Id[l]);
return;
}
if (tr[t][0])
find(I,tr[t][0],l,mid,y);
if (tr[t][1] && mid<y)
find(I,tr[t][1],mid+1,r,y);
}
void change(int I,int t,int l,int r,int x)
{
int mid=(l+r)/2;
if (l==r) return;
_New(t,I);
if (x<=mid)
{
if (l<mid)
New(t,0);
else
tr[t][0]=Id[l];
change(I,tr[t][0],l,mid,x);
}
else
{
if (mid+1<r)
New(t,1);
else
tr[t][1]=Id[r];
change(I,tr[t][1],mid+1,r,x);
}
}
void tj(int t)
{
int i;
bz[t]=1;
dfn[t]=++j;
low[t]=j;
d[++k]=t;
for (i=ls[t]; i; i=a[i][1])
{
if (!dfn[a[i][0]])
tj(a[i][0]);
if (bz[a[i][0]])
low[t]=min(low[t],low[a[i][0]]);
}
if (dfn[t]==low[t]) for (ans[0]=0; d[k+1]!=t; bz[d[k]]=0,ans[++ans[0]]=d[k--]);
}
void work()
{
int i,j,root;
sort(A+1,A+n+1,cmp);
fo(i,1,n)
Id[A[i].b]=A[i].id;
root=++len;
fo(i,1,n)
{
find(A[i].id,root,1,n,A[i].b-1);
len++;
tr[len][0]=tr[root][0];
tr[len][1]=tr[root][1];
root=len;
change(A[i].id,root,1,n,A[i].b);
}
fo(j,1,n) swap(A[j].a,A[j].b),swap(A[j].b,A[j].c);
}
int main()
{
freopen("champion.in","r",stdin);
freopen("champion.out","w",stdout);
scan();
init();
fo(i,1,3) work();
j=0;
k=0;
fo(i,1,n)
if (!dfn[i])
tj(i),k=0;
fo(i,1,ans[0])
if (ans[i]<=n)
Ans[++Ans[0]]=ans[i];
sort(Ans+1,Ans+Ans[0]+1,Cmp);
fo(i,1,Ans[0])
printf("%d\n",Ans[i]);
fclose(stdin);
fclose(stdout);
return 0;
}