「BZOJ1924」「SDOI2010」 所驼门王的宝藏 tarjan + dp(DAG 最长路）

「BZOJ1924」[SDOI2010] 所驼门王的宝藏 tarjan + dp(DAG 最长路）

-----------------------------------------------------------------------------------------------------------------------------------------

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

-----------------------------------------------------------------------------------------------------------------------------------------

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

1、“横天门”:由该门可以传送到同行的任一宫室;

2、“纵寰门”:由该门可以传送到同列的任一宫室;

3、“自 由 门”:由该门可以传送到以该门所在宫室为中心周围8格中任一宫室(如果目标宫室存在的话)。

深谋远虑的 Henry 当然事先就搞到了所驼门王当年的宫殿招标册,书册上详细记录了每扇传送门所属宫室及类型。而且,虽然宫殿内外相隔,但他自行准备了一种便携式传送门,可将自己传送到殿内任意一间宫室开始寻宝,并在任意一间宫室结束后传送出宫。整座宫殿只许进出一次,且便携门无法进行宫室之间的传送。不过好在宫室内传送门的使用没有次数限制,每间宫室也可以多次出入。

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

-----------------------------------------------------------------------------------------------------------------------------------------
「输入格式」

第一行给出三个正整数 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,所有的传送门位置互不相同。

「输出格式」

只有一个正整数,表示你确定的路线所经过不同藏宝 宫室的最大数目。

-----------------------------------------------------------------------------------------------------------------------------------------
「样例输入输出」

输入

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

输出

9

-----------------------------------------------------------------------------------------------------------------------------------------
「数据规模和约定」

测试点编号 N R C
1 16 20 20
2 300 1,000 1,000
3 500 100,000 100,000
4 2,500 5,000 5,000
5 50,000 5,000 5,000
6 50,000 1,000,000 1,000,000
7 80,000 1,000,000 1,000,000
8 100,000 1,000,000 1,000,000
9 100,000 1,000,000 1,000,000
10 100,000 1,000,000 1,000,000

-----------------------------------------------------------------------------------------------------------------------------------------

题解

考虑到 N <= 100,000 矩阵上的点应该很稀疏
于是考虑直接建图做

但是连边的时候将一行中所有“横天门”连在一起 复杂度是 O(N^2）
将一列中所有“纵寰门”连在一起 复杂度也是 O(N^2)

于是考虑将一行中所有“横天门”缩成一个点（将每个“横天门”向该行第一个“横天门”连双向边）
将一列中所有“纵寰门”也缩成一个点（同理）

复杂度就可以降为线性的

再在建好的图上跑 tarjan 把图转为 DAG 再跑 有向无环图的最长路 （dp）

 

  1 #include<bits/stdc++.h>
  2 using namespace std;
  3 const int MAXN=1e5+5;
  4 const int MAXE=1e6+6;
  5 int dx[8]={0,0,1,1,1,-1,-1,-1},dy[8]={1,-1,0,1,-1,0,1,-1};
  6 int edge,nex[MAXE],head[MAXN],tail[MAXE];
  7 int redge,rnex[MAXE],rhead[MAXN],rtail[MAXE];
  8 int n,m,r,c;
  9 struct rec{ int x,y,t; } loc[MAXN];
 10 int f[MAXN];
 11 map<int,int> mp[MAXE];
 12 vector<int> a[MAXE],b[MAXE];
 13 int q[MAXN],top=0,scc=0;
 14 int use[MAXN],cmp[MAXN],low[MAXN],dfn[MAXN],num[MAXN],cnt,ans;
 15 void add(int u,int v){
 16     if (u==v) return;
 17     edge++,nex[edge]=head[u],head[u]=edge,tail[edge]=v;
 18 }
 19 void radd(int u,int v){
 20     redge++,rnex[redge]=rhead[u],rhead[u]=redge,rtail[redge]=v;
 21 }
 22 void clean(){ memset(use,0,sizeof(use)); }
 23 void build(){
 24     for (int i=1;i<=r;i++){
 25         int x=0;
 26         for (int j=0;j<a[i].size();j++) if (loc[a[i][j]].t==1){ x=a[i][j]; break; }
 27         for (int j=0;j<a[i].size();j++){
 28             add(x,a[i][j]);
 29             if (loc[a[i][j]].t==1) add(a[i][j],x);
 30         }
 31     }
 32     for (int i=1;i<=c;i++){
 33         int x=0;
 34         for (int j=0;j<b[i].size();j++) if (loc[b[i][j]].t==2){ x=b[i][j]; break; }
 35         for (int j=0;j<b[i].size();j++){
 36             add(x,b[i][j]);
 37             if (loc[b[i][j]].t==2) add(b[i][j],x);
 38         }
 39     }
 40     for (int i=1;i<=n;i++){
 41         if (loc[i].t!=3) continue;
 42         int x=loc[i].x,y=loc[i].y;
 43         for (int j=0;j<8;j++){
 44             if (mp[x+dx[j]][y+dy[j]]>0) add(i,mp[x+dx[j]][y+dy[j]]);
 45         }
 46     }
 47 }
 48 void tarjan(int u){
 49     low[u]=dfn[u]=++cnt;
 50     q[++top]=u,use[u]=1;
 51     for (int e=head[u];e;e=nex[e]){
 52         int v=tail[e];
 53         if (!dfn[v]){
 54             tarjan(v); low[u]=min(low[u],low[v]);
 55         }
 56         else if (use[v]) low[u]=min(low[u],dfn[v]);
 57     }
 58     if (low[u]==dfn[u]){
 59         int x=0; scc++;
 60         while (x!=u){
 61             x=q[top--];use[x]=0;
 62             cmp[x]=scc;num[scc]++;
 63         }
 64     }
 65 }
 66 void rebuild(){
 67     for (int u=1;u<=n;u++){
 68         for (int e=head[u];e;e=nex[e]){
 69             int v=tail[e];
 70             if (cmp[u]!=cmp[v]) radd(cmp[u],cmp[v]);
 71         }
 72     }
 73 }
 74 void dp(int u){
 75     use[u]=1;
 76     for (int e=rhead[u];e;e=rnex[e]){
 77         int v=rtail[e];
 78         if (!use[v]) dp(v);
 79         f[u]=max(f[u],f[v]);
 80     }
 81     f[u]+=num[u];
 82     ans=max(ans,f[u]);
 83 }
 84 int main(){
 85     scanf("%d%d%d",&n,&r,&c);
 86     for (int i=1;i<=n;i++){
 87         scanf("%d%d%d",&loc[i].x,&loc[i].y,&loc[i].t);
 88         a[loc[i].x].push_back(i);
 89         b[loc[i].y].push_back(i);
 90         mp[loc[i].x][loc[i].y]=i;
 91     } 
 92     build();
 93     clean();
 94     for (int i=1;i<=n;i++){
 95         if (!dfn[i]) tarjan(i);
 96     }
 97     rebuild();
 98     clean();
 99     for (int i=1;i<=scc;i++){
100         if (!use[i]) dp(i);
101     }
102     printf("%d\n",ans);
103     return 0;
104 } 

View Code

 

转载于:https://www.cnblogs.com/shjrd-dlb/p/9043847.html