【bzoj1924】[Sdoi2010]所驼门王的宝藏（tarjan+STL+dp）

题目：我是超链接

题意：

一个r*c的图中，有n个宫殿。 每个宫殿有一个类型。 

类型1:可以到达他所在的行的任意宫殿。 

类型2:可以到达他所在的列的任意宫殿。 

类型3:可以到达他四周八个格子的任意宫殿。 

现在你从任意一个宫殿开始，询问你最多访问多少个宫殿。 

题解：

不知道每一行每一列确实有多少个，那用vector真是再好不过了

map<int,int>mp[1e7];最多有1e7行，在本题是一个行的int对应一个列的int，赋值为一个编号，表示这个点。本题的目的是找8个方向的门是否存在（不存在就不用连了啊，反正只能访问这一个门）【map<int>mp[1e7] mp[x[i]]=y[i]可以去掉重复的门，这道题用不到

是每一行选定一个横天门A，向该行横天门B连双向边，其余门单向边，这样B与横天门A成强连通，就不用再给B连一边A连过的边了，反正都能到。纵列同理

空间太大的话会T哦，不记忆化搜索的话会WA哦（虽然不知道为什么

代码：

#include <cstdio>
#include <iostream>
#include <map>
#include <vector>
#include <cstring>
#include <algorithm>
#define R 1000005
#define N 100005
using namespace std;
map<int,int>mp[R];
vector<int>a[R],b[R];
int cc[8][2]={{1,0},{-1,0},{0,1},{0,-1},{1,1},{-1,-1},{-1,1},{1,-1}};
int r,c,n,ans,tot,nxt[R],point[N],v[R],id[N],x[N],y[N],Dfn[N],Low[N],NN,strack[N],tmp;
int num,belong[N],size[N],deep[N],in[N],tot1,nxt1[R],point1[N],v1[R];
bool vis[R];
void addline(int x,int y){if (x==y) return;++tot; nxt[tot]=point[x]; point[x]=tot; v[tot]=y;}
void addline1(int x,int y){++tot1; nxt1[tot1]=point1[x]; point1[x]=tot1; v1[tot1]=y;}
void build()
{
    int i,j;
    for (i=1;i<=r;i++)
    {
        int x=0,t=a[i].size();
        for (j=0;j<t;j++)
          if (id[a[i][j]]==1) {x=a[i][j]; break;}
        for (j=0;j<t;j++)//和这一行的所有点相连咯 
        {
            addline(x,a[i][j]);
            if (id[a[i][j]]==1) addline(a[i][j],x);
        }
    }
    for (i=1;i<=c;i++)
    {
        int x=0,t=b[i].size();
        for (j=0;j<t;j++)
          if (id[b[i][j]]==2) {x=b[i][j]; break;}
        for (j=0;j<t;j++)
        {
            addline(x,b[i][j]);
            if (id[b[i][j]]==2) addline(b[i][j],x);
        }
    }
    for (i=1;i<=n;i++)
      if (id[i]==3)
        for (j=0;j<8;j++)
        {
            int t=mp[x[i]+cc[j][0]][y[i]+cc[j][1]];
            if (t) addline(i,t);//如果此处有宫室 
        }
}
void tarjan(int now)
{
    Dfn[now]=Low[now]=++NN; vis[now]=1; strack[++tmp]=now;
    for (int i=point[now];i;i=nxt[i])
      if (!Dfn[v[i]])
      {
        tarjan(v[i]);
        Low[now]=min(Low[now],Low[v[i]]);
      }
      else if (vis[v[i]]) Low[now]=min(Low[now],Dfn[v[i]]);
    if (Low[now]==Dfn[now])
    {
        num++;
        while (strack[tmp]!=now)
        {
            belong[strack[tmp]]=num;
            vis[strack[tmp]]=0; tmp--;
            size[num]++;
        }
        belong[strack[tmp]]=num;
        vis[strack[tmp]]=0; tmp--;
        size[num]++;
    }
}
void rebuild()
{
    for (int i=1;i<=n;i++)
      for (int j=point[i];j;j=nxt[j])
        if (belong[v[j]]!=belong[i]) 
          addline1(belong[i],belong[v[j]]),in[belong[v[j]]]++;
}
void dp(int now)
{
	if (deep[now]) return;
    for (int i=point1[now];i;i=nxt1[i])
    {
        dp(v1[i]);
        deep[now]=max(deep[now],deep[v1[i]]);
    }
    deep[now]+=size[now];
    ans=max(ans,deep[now]);
}
int main()
{
    int i;
    scanf("%d%d%d",&n,&r,&c);
    for (i=1;i<=n;i++)
    {
        scanf("%d%d%d",&x[i],&y[i],&id[i]);
        mp[x[i]][y[i]]=i;
        a[x[i]].push_back(i);
        b[y[i]].push_back(i);
    }
    build();
    for (i=1;i<=n;i++) 
      if (!Dfn[i]) tarjan(i);
    rebuild();
    for (i=1;i<=num;i++)
      if (!in[i]) dp(i);
    printf("%d",ans);
}