BZOJ1924_所驼门王的宝藏_KEY

题目传送门

这道题苟了我好久，因为链表的内存问题，之后再细讲。

首先这是一道Tarjan+DAG上DP的题目。

有三种门，对于每种门可以和其他门相连。即连边。

使用链表快速查询连边。

建完图后可以进行Tarjan缩点。

然后做一遍DAG上DP就好了。（记搜）

然后因为建图时会有很多条边，而行列最多只有100000个，所以要分开定义。

不然会爆内存。

code：

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

int read()
{
    char c;while(c=getchar(),c<''||c>'');
    int x=c-'';while(c=getchar(),c>=''&&c<='')x=x*+c-'';
    return x;
}

int N,R,C,x,y,fx,fy,o,a[][];

struct list{
    int head[],nxt[],To[],cnt;
    list(){memset(head,-,sizeof head);memset(nxt,-,sizeof nxt);}

    void add(int x,int y){
        To[cnt]=y;
        nxt[cnt]=head[x];
        head[x]=cnt;
        cnt++;
    }
}M,NM;

struct list2{
    int head[],nxt[],To[],cnt;
    list2(){memset(head,-,sizeof head);memset(nxt,-,sizeof nxt);}

    void add(int x,int y){
        To[cnt]=y;
        nxt[cnt]=head[x];
        head[x]=cnt;
        cnt++;
    }
}A,B;

int DFN[],LOW[],stack[],top,cnt,vis[],fa[];
int otk[],f[],sk[],ans;

void tarjan(int now)
{
    stack[++top]=now;
    DFN[now]=LOW[now]=++cnt;
    vis[now]=;
        for(int i=M.head[now];i!=-;i=M.nxt[i]){
            if(!DFN[M.To[i]])tarjan(M.To[i]),LOW[now]=min(LOW[now],LOW[M.To[i]]);
            else if(vis[M.To[i]])LOW[now]=min(LOW[now],DFN[M.To[i]]);
        }
    if(DFN[now]==LOW[now]){
        while(stack[top]!=now)
            fa[stack[top]]=now,sk[now]++,vis[stack[top]]=,top--;
        fa[stack[top]]=now;
        sk[now]++;
        vis[stack[top--]]=;
    }
    return ;
}

int dfs(int now)
{
    if(vis[now])return f[now];
    vis[now]=;
        for(int i=NM.head[now];i!=-;i=NM.nxt[i]){
            f[now]=max(f[now],dfs(NM.To[i]));
        }
    f[now]+=sk[now];
    ans=max(ans,f[now]);
    return f[now];
}

int main()
{
    N=read(),R=read(),C=read();
    register int i,j,k;
        for(i=;i<=N;i++){
            a[i][]=read(),a[i][]=read(),a[i][]=read();
            A.add(a[i][],i);B.add(a[i][],i);
        }
        for(i=;i<=N;i++){
            x=a[i][],y=a[i][],o=a[i][];
            if(o==)
                for(j=A.head[x];j!=-;j=A.nxt[j])
                    if(A.To[j]!=i)M.add(i,A.To[j]);
            if(o==)
                for(j=B.head[y];j!=-;j=B.nxt[j])
                    if(B.To[j]!=i)M.add(i,B.To[j]);
            if(o==){
                for(k=-;k<;k++)
                for(j=A.head[x+k];j!=-;j=A.nxt[j]){
                    fx=a[A.To[j]][],fy=a[A.To[j]][];
                    if(fx<=x+&&fx>=x-&&fy<=y+&&fy>=y-&&A.To[j]!=i)
                        M.add(i,A.To[j]);
                }
            }
        }
        for(i=;i<=N;i++)if(!DFN[i])tarjan(i);
    memset(vis,,sizeof vis);
        for(i=;i<=N;i++){
            if(!vis[i]){vis[i]=;
                for(j=M.head[i];j!=-;j=M.nxt[j]){
                    if(fa[M.To[j]]!=fa[i])
                    otk[fa[M.To[j]]]++;
                }
            }
        }
        for(i=;i<=N;i++)
            for(j=M.head[i];j!=-;j=M.nxt[j])
                if(fa[i]!=fa[M.To[j]])NM.add(fa[i],fa[M.To[j]]);
    memset(vis,,sizeof vis);
        for(i=;i<=N;i++)
            if(!otk[fa[i]]&&!vis[fa[i]])
                dfs(i);
    printf("%d",ans);
    return ;
}

Vector版传送门