JZOJ5705. 【gdsoi2018 day3】谁是冠军

题目描述

你猜
http滑稽://172.1滑稽6.0.132/sen滑稽ior/#contest滑稽/show/237滑稽2/8

30%

暴力？
反正我写的暴力1分都拿不了

60%

不知道

80%

不知道

100%

其实可以直接连边+Tarjan，但是慢在连边$O(n^2)$
于是考虑用线段树优化。

因为Tarjan找的是强连通分量，也就是说a–>b和a–>c–>b是等价的

而题目要求的是abc中的两项，所以可以分开来做三遍
按照一维排序，另一维丢到权值线段树里

把线段树的每个区间都当作一个点，插入一个点时把所有包含它的区间跟它连边，再把它跟所有它包含的区间连边
显然两者的个数都是$\log n$的
那么一共就连了$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;
}