省选专练SDOI2010所驼门王的宝藏

这是一个奇怪的题

首先，其思路和APIO2009-3掠夺计划相比难得多，数据变大了啊。

于是两两建边不可行，考虑用拉链判重。

从横天门或纵寰门引出的边可能特别多。在极端情况下，10^5个横天门在同一行里出现，这样时空复杂度都是无法承受的。考虑这一点进行优化：由于在同一行的横天门一定属于同一个强连通分量，所以在建边时，只需要对在同一行的横天门构建出一个环即可，不需要两两进行连边。而此行内的其他宫室，只需要从这个环中的任意一点向这个宫室连边即可。对于纵寰门也是一样。

然后缩了点想干嘛就干嘛top还是SPFA都可以

#include<iostream>#include<cstdio>#include<stack>#include<algorithm>#include<cmath>#include<cstring>#include<queue>#include<map>using namespace std;inline void read(int &x){    int f=1;    x=0;    char ch=getchar();    while(ch<'0'||ch>'9'){        if(ch=='-')            f=-1;        ch=getchar();    }    while(ch<='9'&&ch>='0'){        x=x*10+ch-'0';        ch=getchar();    }    x*=f;}int dx[]={0,1,0,-1,1,1,-1,-1};int dy[]={1,0,-1,0,-1,1,-1,1};int n,r,c;struct Front_star{    int u,v,w,nxt;}edge[5000000],e[5000000],ee[5000000];int cnt=0;int first[500000]={0};int head[500000]={0};struct Node{    int x,y,t;}pre[2000007]/*first array*/,nxt[2000007],a[2000007];void addedge(int u,int v,int w){    cnt++;    e[cnt].u=u;    e[cnt].v=v;    e[cnt].w=w;    e[cnt].nxt=first[u];    first[u]=cnt;}void adde(int u,int v,int w){    cnt++;    ee[cnt].u=u;    ee[cnt].v=v;    ee[cnt].w=w;    ee[cnt].nxt=head[u];    head[u]=cnt;}void add(){    for(int i=1;i<=n;i++){        if(a[i].t==1){            for(int j=pre[a[i].x].x;j;j=nxt[j].x){                if(i==j){                    continue;                }                addedge(i,j,1);            }        }        if(a[i].t==2){            for(int j=pre[a[i].y].y;j;j=nxt[j].y){                if(i==j){                    continue;                }                addedge(i,j,1);            }        }        if(a[i].t==3){            for(int k=0;k<8;k++){                int xx=a[i].x+dx[k];                for(int j=pre[xx].x;j;j=nxt[j].x){                    if(i==j)                        continue;                    if(abs(a[j].y-a[i].y)<=1)                    addedge(i,j,1);                }             }        }    }}int tot=0;//--------------int low[500000]={0};int dfn[500000]={0};int color[500000]={0};int siz[500000]={0};int inqueue[500000]={0};int scc=0;stack<int> S;void tarjan(int u){    tot++;    low[u]=dfn[u]=tot;    S.push(u);    inqueue[u]=1;    for(int i=first[u];i;i=e[i].nxt){        int v=e[i].v;        if(!dfn[v]){            tarjan(v);            low[u]=min(low[u],low[v]);        }        else{            if(inqueue[v]){                low[u]=min(low[u],dfn[v]);            }        }    }    if(low[u]==dfn[u]){        int x;        scc++;        do{            x=S.top();            S.pop();            color[x]=scc;            siz[scc]++;            inqueue[x]=0;        }while(x!=u);    }}int rd[500001]={0};int cd[500001]={0};int dis[500001]={0};void SPFA(){    queue<int>q;    for(int i=1;i<=scc;i++){        if(rd[i]==0){            q.push(i);            dis[i]=siz[i];        }    }    while(!q.empty()){        int x=q.front();        q.pop();        for(int i=head[x];i;i=ee[i].nxt){            int v=ee[i].v;            if(rd[v]==0)                continue;            rd[v]--;            if(dis[v]<dis[x]+siz[v]){                dis[v]=dis[x]+siz[v];            }            if(!rd[v]){                q.push(v);            }        }    }}int main(){    read(n);    read(r);    read(c);    for(int i=1;i<=n;i++){        int x,y,t;        read(x);        read(y);        read(t);        a[i].x=x;        a[i].y=y;        a[i].t=t;        nxt[i].x=pre[x].x;        pre[x].x=i;        nxt[i].y=pre[y].y;        pre[y].y=i;    }    add();    for(int i=1;i<=n;i++){        if(!dfn[i]){            tarjan(i);        }    }//cnt=0;//memset(first,0,sizeof(first));    for(int i=1;i<=n;i++){        for(int j=first[i];j;j=e[j].nxt){            int v=e[j].v;            int u=i;            if(color[u]!=color[v]){                adde(color[u],color[v],1);                rd[color[v]]++;                cd[color[u]]++;            }        }    }//for(int i=1;i<=scc;i++)//cout<<rd[i]<<" "<<siz[i]<<endl;    SPFA();    int ans=-1;//cout<<"4q382436939486724389279834"<<endl;    for(int i=1;i<=scc;i++){        ans=max(ans,dis[i]);//cout<<siz[i]<<" "<<cd[i]<<endl;    }    cout<<ans;}