bzoj1924&luogu2403 所驼门王的宝藏

http://www.elijahqi.win/archives/1010
题目描述

在宽广的非洲荒漠中，生活着一群勤劳勇敢的羊驼家族。被族人恭称为“先知”的Alpaca L. Sotomon是这个家族的领袖，外人也称其为“所驼门王”。所驼门王毕生致力于维护家族的安定与和谐，他曾亲自率军粉碎河蟹帝国主义的野蛮侵略，为族人立下赫赫战功。所驼门王一生财宝无数，但因其生性节俭低调，他将财宝埋藏在自己设计的地下宫殿里，这也是今天Henry Curtis故事的起点。Henry是一个爱财如命的贪婪家伙，而又非常聪明，他费尽心机谋划了这次盗窃行动，破解重重机关后来到这座地下宫殿前。

整座宫殿呈矩阵状，由R×C间矩形宫室组成，其中有N间宫室里埋藏着宝藏，称作藏宝宫室。宫殿里外、相邻宫室间都由坚硬的实体墙阻隔，由一间宫室到达另一间只能通过所驼门王独创的移动方式——传送门。所驼门王为这N间藏宝宫室每间都架设了一扇传送门，没有宝藏的宫室不设传送门，所有的宫室传送门分为三种：

“横天门”：由该门可以传送到同行的任一宫室；
“纵寰门”：由该门可以传送到同列的任一宫室；
“自由门”：由该门可以传送到以该门所在宫室为中心周围8格中任一宫室（如果目标宫室存在的话）。
深谋远虑的Henry当然事先就搞到了所驼门王当年的宫殿招标册，书册上详细记录了每扇传送门所属宫室及类型。而且，虽然宫殿内外相隔，但他自行准备了一种便携式传送门，可将自己传送到殿内任意一间宫室开始寻宝，并在任意一间宫室结束后传送出宫。整座宫殿只许进出一次，且便携门无法进行宫室之间的传送。不过好在宫室内传送门的使用没有次数限制，每间宫室也可以多次出入。

现在Henry已经打开了便携门，即将选择一间宫室进入。为得到尽多宝藏，他希望安排一条路线，使走过的不同藏宝宫室尽可能多。请你告诉Henry这条路线最多行经不同藏宝宫室的数目。

输入输出格式

输入格式：

输入文件sotomon.in第一行给出三个正整数N, R, C。

以下N行，每行给出一扇传送门的信息，包含三个正整数xi, yi, Ti，表示该传送门设在位于第xi行第yi列的藏宝宫室，类型为Ti。Ti是一个1~3间的整数，1表示可以传送到第xi行任意一列的“横天门”，2表示可以传送到任意一行第yi列的“纵寰门”，3表示可以传送到周围8格宫室的“自由门”。

保证1≤xi≤R，1≤yi≤C，所有的传送门位置互不相同。

输出格式：

输出文件sotomon.out只有一个正整数，表示你确定的路线所经过不同藏宝宫室的最大数目。

输入输出样例

输入样例#1：

10 7 7
2 2 1
2 4 2
1 7 2
2 7 3
4 2 2
4 4 1
6 7 3
7 7 1
7 5 2
5 2 1
输出样例#1：

9
说明

数据规模和约定：

在洛谷map会被卡内存 把map开到10^5就能水过 不知为何

这题和上一题一样，就是建图非常困难

建完图，缩点之后成为DAG 然后跑拓扑排序 使得其无后效性然后递推求出最大值

建图的时候 每一行每一列当op==1||op==2的时候 我们因为无论如何也是会缩成一个点所以我们只需要单向连接即可，这样时间满足条件，空间满足条件

然后如果再次出现op==1||op==2的边我们则往回连一条边，这样可以使得tarjan的时候缩点成为一个点 bzoj 正常过

#include<cstdio>
#include<map>
#include<vector>
#include<queue>
#include<stack>
#define N 110000
using namespace std;
int x1[10]={0,-1,-1,-1,0,1,1,1,0};
int y1[10]={0,-1,0,1,1,1,0,-1,-1};
inline int read(){
    int x=0;char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while (ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=getchar();}
    return x;
}
map<int,int> mm[10*N];
vector<int> q1[10*N],q2[10*N];
stack<int> st;queue<int>q;
struct node{
    int x,y,next;
}data[N*10],data1[N*10];
int dfn[N],low[N],h[N],s,num,n,r,c,x[N],y[N],op[N],b[N],size[N],h1[N],in[N],f[N];bool stackf[N];
inline void insert1(int x,int y){
    if (x==y) return;
    data[++num].y=y;data[num].next=h[x];h[x]=num;data[num].x=x;
}
void tarjan(int x){
    dfn[x]=low[x]=++num;stackf[x]=true;st.push(x);
    for (int i=h[x];i;i=data[i].next){
        int y=data[i].y;
        if (!dfn[y]) tarjan(y),low[x]=min(low[x],low[y]);else if (stackf[y]) low[x]=min(low[x],dfn[y]);
    }
    if (dfn[x]==low[x]){
        ++s;int y;
        do{
            y=st.top();st.pop();b[y]=s;size[s]++;stackf[y]=false;
        }while(y!=x);
    }
}
inline void insert2(int x,int y){
    data1[++num].y=y;data1[num].next=h1[x];h1[x]=num;
}
int main(){
    //freopen("bzoj1924.in","r",stdin);
    n=read();r=read();c=read();
    for (int i=1;i<=n;++i) {
        x[i]=read();y[i]=read();op[i]=read();
        mm[x[i]][y[i]]=i;q1[x[i]].push_back(i);q2[y[i]].push_back(i);
    }
    for (int i=1;i<=r;++i){
        int nm=q1[i].size(),now=0;
        for (int j=0;j<nm;++j) if (op[q1[i][j]]==1) {now=q1[i][j];break;}
        for (int j=0;j<nm;++j){insert1(now,q1[i][j]);if (op[q1[i][j]]==1&&now!=q1[i][j]) insert1(q1[i][j],now);}    
    }
    for (int i=1;i<=c;++i){
        int nm=q2[i].size(),now=0;
        for (int j=0;j<nm;++j) if (op[q2[i][j]]==2) {now=q2[i][j];break;}
        for (int j=0;j<nm;++j){insert1(now,q2[i][j]);if (op[q2[i][j]]==2&&now!=q2[i][j]) insert1(q2[i][j],now);}    
    }
    for (int i=1;i<=n;++i)
        if (op[i]==3)
            for (int j=1;j<=8;++j)
                if (mm[x[i]+x1[j]][y[i]+y1[j]]) insert1(i,mm[x[i]+x1[j]][y[i]+y1[j]]);
    num=0;        
    for (int i=1;i<=n;++i) if(!dfn[i]) tarjan(i);num=0;
    //for (int i=1;i<=n;++i) printf("%d ",b[i]);
    for (int x=1;x<=n;++x){
        for (int j=h[x];j;j=data[j].next){
            int y=data[j].y;
            if (b[x]!=b[y]) insert2(b[x],b[y]),in[b[y]]++;
        }
    }
    int ans=0;
    for (int i=1;i<=s;++i) if (!in[i]) q.push(i),f[i]=size[i];
    while (!q.empty()){
        int x=q.front();q.pop();
        for (int i=h1[x];i;i=data1[i].next){
            int y=data1[i].y;
            f[y]=max(f[y],f[x]+size[y]);
            ans=max(ans,f[y]);
            if (!--in[y]) q.push(y);
        }
    }
    printf("%d\n",ans);
    return 0;
}