bzoj1924 [Sdoi2010]所驼门王的宝藏(tarjan缩点+拓扑排序+dp)

时间:2021-02-12 22:07:57

把每个点向它能到达的点连一条有向边,那么答案就是这个图的最长链。然后就是一个强连通分量缩点+拓扑排序dp求最长路的套路了。但是,边不能直接暴力连!比如 同一行可能会有好多横天门,然后 就爆掉了!!!其实如果一行有许多横天门,我们不需要对每个横天门都连许多边,因为同一行的横天门肯定是属于同一个强连通的,反正我们也要给他们缩点,所以我们只拿一个横天门出来做,让他去连边,对于同一行的横天门连双向边,其他门连单向边即可。
tips:洛谷数据极强,推荐大家一试。

#include <bits/stdc++.h>
using namespace std;
#define inf 0x3f3f3f3f
#define ll long long
#define N 100010
#define M 1000010
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,m,kk,num=0,h[N],h1[N],dfn[N],low[N],dfnum=0,bel[N],scc=0,sz[N],qq[N];
int dx[]={1,1,1,-1,-1,-1,0,0},dy[]={-1,0,1,-1,0,1,-1,1},in[N],dp[N],top=0,ans=0;
struct node{
    int x,y,op;
}p[N];
vector<int>a[M];
vector<int>b[M];
map<int,int>mp[N];
struct edge{
    int to,next; 
}data[M],data1[M];
inline void add(int x,int y){
    data[++num].to=y;data[num].next=h[x];h[x]=num;
}
inline void add1(int x,int y){
    data1[++num].to=y;data1[num].next=h1[x];h1[x]=num;in[y]++;
}
bool inq[N];
void tarjan(int x){
    dfn[x]=low[x]=++dfnum;qq[++top]=x;inq[x]=1;
    for(int i=h[x];i;i=data[i].next){
        int y=data[i].to;
        if(!dfn[y]) tarjan(y),low[x]=min(low[x],low[y]);
        else if(inq[y]) low[x]=min(low[x],dfn[y]);
    }
    if(low[x]==dfn[x]){
        ++scc;while(1){
            int y=qq[top--];inq[y]=0;
            bel[y]=scc;sz[scc]++;if(x==y) break;
        }
    }
}
inline void rebuild(){
    num=0;
    for(int x=1;x<=kk;++x){
        for(int i=h[x];i;i=data[i].next){
            int y=data[i].to;
            if(bel[x]!=bel[y]) add1(bel[x],bel[y]);
        }
    }
}
void Topodp(){
    queue<int>q;
    for(int i=1;i<=scc;++i)
        if(!in[i]) q.push(i),dp[i]=sz[i];
    while(!q.empty()){
        int x=q.front();q.pop();
        for(int i=h1[x];i;i=data1[i].next){
            int y=data1[i].to;
            dp[y]=max(dp[y],dp[x]+sz[y]);
            ans=max(ans,dp[y]);
            if(--in[y]==0) q.push(y);
        }
    }
}
int main(){
// freopen("a.in","r",stdin);
    kk=read();n=read();m=read();
    for(int i=1;i<=kk;++i){
        p[i].x=read(),p[i].y=read(),p[i].op=read();
        mp[p[i].x][p[i].y]=i;a[p[i].x].push_back(i);b[p[i].y].push_back(i);
    }
    for(int i=1;i<=n;++i){
        int x=0,t=a[i].size();
        while(x<t&&p[a[i][x]].op!=1) ++x;//找到第一个横天门
        if(x==t) continue;
        for(int j=0;j<t;++j){
            if(j==x) continue;
            add(a[i][x],a[i][j]);//连单向边
            if(p[a[i][j]].op==1) add(a[i][j],a[i][x]);//两个横天门之间连双向边
        }
    }
    for(int i=1;i<=m;++i){
        int x=0,t=b[i].size();
        while(x<t&&p[b[i][x]].op!=2) ++x;
        if(x==t) continue;
        for(int j=0;j<t;++j){
            if(j==x) continue;
            add(b[i][x],b[i][j]);
            if(p[b[i][j]].op==2) add(b[i][j],b[i][x]);
        }
    }
    for(int i=1;i<=kk;++i)
        if(p[i].op==3){
            for(int k=0;k<8;++k){
                int xx=p[i].x+dx[k],yy=p[i].y+dy[k];
                if(xx<1||xx>n||yy<1||yy>m||!mp[xx][yy]) continue;
                add(i,mp[xx][yy]);
            }
        }
    for(int i=1;i<=kk;++i) if(!dfn[i]) tarjan(i);
    rebuild();Topodp();
    printf("%d\n",ans);
    return 0;
}