【bzoj1924】【SDOI2010】所驼门王的宝藏

最新推荐文章于 2020-08-10 20:03:21 发布
最新推荐文章于 2020-08-10 20:03:21 发布
阅读量306 收藏 1
点赞数
分类专栏： 题解 dp 拓扑排序 文章标签： c++ bzoj 省选
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_38064038/article/details/78158391
版权
题解 同时被 3 个专栏收录
56 篇文章 0 订阅
订阅专栏
dp
14 篇文章 0 订阅
订阅专栏
拓扑排序
2 篇文章 0 订阅
订阅专栏

题意

  一个 r∗c 的矩阵中，有 n 个特殊点，每个特殊点有一个传送门（单向），共三种：
  ①.到达同一行的特殊点（横天门）
  ②.到达同一列的特殊点（纵寰门）
  ③.到达同一九宫格的特殊点（自由门）
  求最多可以到达多少个特殊点

解法

Tarjan缩点+拓扑排序+图上 DP ：
  和之前的考试的 teach 一题非常相似，不过这一题不能够使用 spfa ，而必须使用 DP
  而为了消除后效性，所以可以按照拓扑序进行 DP
  设 fu 表示到达 u 点（Tarjan缩点之后的编号）能获得的最大点数，那么有：
   fson=max（fson，fu+sizson） ， son 为 u 的孩子，sizson为 son 的联通块的点数
  另外，本题的建图是一个重点：
  容易想到，如果在第 i 个藏宝宫室可以传送到第j个藏宝宫室，就连一条 i−>j 的有向边， Tarjan 强连通分量缩点之后，在新图上找出一条路径，使得这条路径上的所有点对应的强连通分量大小之和最大，这个可以按照拓扑序进行递推来实现。
  但是，在建图上存在一个问题，从横天门或纵寰门引出的边可能特别多。在极端情况下， 105 个横天门在同一行里出现，这样时空复杂度都是无法承受的。考虑这一点进行优化：由于在同一行的横天门一定属于同一个强连通分量，所以在建边时，只需要对在同一行的横天门构建出一个环即可，不需要两两进行连边。而此行内的其他宫室，只需要从这个环中的任意一点向这个宫室连边即可。对于纵寰门也是一样。
  此外，在建图的实现上有一些小技巧，具体见代码。

复杂度

O（ knlogn ）， k <script type="math/tex" id="MathJax-Element-23">k</script>为常数

70分spfa代码

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<vector>
#include<queue>
#include<set>
#include<map>
using namespace std;
const int MAXN=100010;
struct point
{
    int x,y;
    void gi() { scanf("%d%d",&x,&y); }
    bool operator < (const point a) const
    {
        return x==a.x ? y<a.y : x<a.x;
    }
}p[MAXN];
struct node
{
    int next,to;
}t[MAXN*2],e[MAXN*12];
int head[MAXN],Head[MAXN],num;
int type[MAXN],siz[MAXN];
int dfn[MAXN],low[MAXN];
int vis[MAXN],dis[MAXN];
int be[MAXN],s[MAXN];
int cnt,top,tot;
int n,r,c,ans;
int dx[9]={0,1,1,1,0,-1,-1,-1,0};
int dy[9]={0,1,0,-1,-1,-1,0,1,1};
vector<int> H[MAXN*10],L[MAXN*10];
vector<int> N[MAXN*10],M[MAXN*10];
map<int,int> id;
queue<int> q;
set<point> P;
void add(int u,int v)
{
    if( u==v )   return ;
    t[++num]=(node){ head[u],v };
    head[u]=num;
}
void Add(int u,int v)
{
    e[++num]=(node){ Head[u],v };
    Head[u]=num;
}
void Tarjan(int k)
{
    int x;
    vis[k]=1,s[++top]=k;
    dfn[k]=low[k]=++cnt;
    for(int i=head[k]; i ;i=t[i].next)
    {
        x=t[i].to;
        if( !dfn[x] )
            Tarjan( x ),low[k]=min( low[k],low[x] );
        else
            if( vis[x] )   low[k]=min( low[k],dfn[x] );
    }
    if( dfn[k]==low[k] )
    {
        tot++;
        do
        {
            x=s[top--],vis[x]=0;
            be[x]=tot,siz[tot]++;
            for(int i=head[x]; i ;i=t[i].next)   if( !be[t[i].to] || be[t[i].to]!=tot )   Add( tot,t[i].to );
        }while( x!=k );
    }
}
void spfa(int s)
{
    int tmp;
    dis[s]=siz[s];
    q.push( s );
    while( !q.empty() )
    {
        tmp=q.front(),q.pop();
        vis[tmp]=0;
        for(int i=Head[tmp],x; i ;i=e[i].next)
        {
            x=be[e[i].to];
            if( x==tmp )   continue ;
            if( dis[x]<dis[tmp]+siz[x] )
            {
                dis[x]=dis[tmp]+siz[x];
                if( !vis[x] )   vis[x]=1,q.push( x );
            }
        }
    }
}
int main()
{
    set<point>::iterator   l;
    int ex,ey;
    scanf("%d%d%d",&n,&r,&c);
    for(int i=1;i<=n;i++)
    {
        p[i].gi(),scanf("%d",&type[i]);
        if( type[i]==1 )   H[p[i].x].push_back( i ),M[p[i].y].push_back( i );
        if( type[i]==2 )   L[p[i].y].push_back( i ),N[p[i].x].push_back( i );
        if( type[i]==3 )   N[p[i].x].push_back( i ),M[p[i].y].push_back( i );
        id[c*(p[i].x-1)+p[i].y]=i;
        P.insert( p[i] );
    }
    for(int i=1,x;i<=r;i++)
    {
        x=H[i].size();
        if( x<2 )   continue ;
        for(int j=1;j<=x-1;j++)   add( H[i][j-1],H[i][j] );
        add( H[i][x-1],H[i][0] );
    }
    for(int i=1,x;i<=c;i++)
    {
        x=L[i].size();
        if( x<2 )   continue ;
        for(int j=1;j<=x-1;j++)   add( L[i][j-1],L[i][j] );
        add( L[i][x-1],L[i][0] );
    }
    for(int i=1,x;i<=r;i++)
    {
        x=N[i].size()-1;
        if( x<0 || !H[i].size() )   continue ;
        for(int j=0;j<=x;j++)   add( H[i][0],N[i][j] );
    }
    for(int i=1,x;i<=c;i++)
    {
        x=M[i].size()-1;
        if( x<0 || !L[i].size() )   continue ;
        for(int j=0;j<=x;j++)   add( L[i][0],M[i][j] );
    }
    for(int i=1,x;i<=n;i++)
    {
        if( type[i]!=3 )   continue ;
        for(int k=1;k<=8;k++)
        {
            ex=p[i].x+dx[k],ey=p[i].y+dy[k];
            if( ex<1 || ex>r || ey<1 || ey>c )   continue ;
            l=P.find( (point){ ex,ey } );
            if( l==P.end() )   continue ;
            add( i,id[c*(ex-1)+ey] );
        }
    }
    num=0;
    for(int i=1;i<=n;i++)   if( !dfn[i] )   Tarjan( i );
    for(int i=1;i<=tot;i++)   vis[i]=0;
    for(int i=1;i<=tot;i++)   if( !dis[i] )   spfa( i );
    for(int i=1;i<=tot;i++)   ans=max( ans,dis[i] );
    printf("%d\n",ans);
    return 0;
}

100分DP代码

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<vector>
#include<queue>
#include<set>
#include<map>
using namespace std;
const int MAXN=100010;
struct point
{
    int x,y;
    void gi() { scanf("%d%d",&x,&y); }
    bool operator < (const point a) const
    {
        return x==a.x ? y<a.y : x<a.x;
    }
}p[MAXN];
struct node
{
    int next,to;
}t[MAXN*2],e[MAXN*12];
int head[MAXN],Head[MAXN],num;
int dfnm[MAXN],deg[MAXN];
int type[MAXN],siz[MAXN];
int dfn[MAXN],low[MAXN];
int vis[MAXN],dis[MAXN];
int be[MAXN],s[MAXN];
int cnt,top,tot;
int n,r,c,ans;
int dx[9]={0,1,1,1,0,-1,-1,-1,0};
int dy[9]={0,1,0,-1,-1,-1,0,1,1};
vector<int> H[MAXN*10],L[MAXN*10];
vector<int> N[MAXN*10],M[MAXN*10];
map<int,int> id;
queue<int> q;
set<point> P;
void add(int u,int v)
{
    if( u==v )   return ;
    t[++num]=(node){ head[u],v };
    head[u]=num;
}
void Add(int u,int v)
{
    e[++num]=(node){ Head[u],v };
    Head[u]=num;
}
void Tarjan(int k)
{
    int x;
    vis[k]=1,s[++top]=k;
    dfn[k]=low[k]=++cnt;
    for(int i=head[k]; i ;i=t[i].next)
    {
        x=t[i].to;
        if( !dfn[x] )
            Tarjan( x ),low[k]=min( low[k],low[x] );
        else
            if( vis[x] )   low[k]=min( low[k],dfn[x] );
    }
    if( dfn[k]==low[k] )
    {
        tot++;
        do
        {
            x=s[top--],vis[x]=0;
            be[x]=tot,siz[tot]++;
            for(int i=head[x]; i ;i=t[i].next)   if( !be[t[i].to] || be[t[i].to]!=tot )   Add( tot,t[i].to );
        }while( x!=k );
    }
}
void topology()
{
    int tmp;
    for(int i=1;i<=tot;i++)   if( !deg[i] )   q.push( i );
    while( !q.empty() )
    {
        tmp=q.front(),q.pop();
        dfn[tmp]=++cnt,dfnm[cnt]=tmp;
        for(int i=Head[tmp],x; i ;i=e[i].next)
        {
            x=be[e[i].to],deg[x]--;
            if( !deg[x] )   q.push( x );
        }
    }
}
int main()
{
    set<point>::iterator   l;
    int ex,ey;
    scanf("%d%d%d",&n,&r,&c);
    for(int i=1;i<=n;i++)
    {
        p[i].gi(),scanf("%d",&type[i]);
        if( type[i]==1 )   H[p[i].x].push_back( i ),M[p[i].y].push_back( i );
        if( type[i]==2 )   L[p[i].y].push_back( i ),N[p[i].x].push_back( i );
        if( type[i]==3 )   N[p[i].x].push_back( i ),M[p[i].y].push_back( i );
        id[c*(p[i].x-1)+p[i].y]=i;
        P.insert( p[i] );
    }
    for(int i=1,x;i<=r;i++)
    {
        x=H[i].size();
        if( x<2 )   continue ;
        for(int j=1;j<=x-1;j++)   add( H[i][j-1],H[i][j] );
        add( H[i][x-1],H[i][0] );
    }
    for(int i=1,x;i<=c;i++)
    {
        x=L[i].size();
        if( x<2 )   continue ;
        for(int j=1;j<=x-1;j++)   add( L[i][j-1],L[i][j] );
        add( L[i][x-1],L[i][0] );
    }
    for(int i=1,x;i<=r;i++)
    {
        x=N[i].size()-1;
        if( x<0 || !H[i].size() )   continue ;
        for(int j=0;j<=x;j++)   add( H[i][0],N[i][j] );
    }
    for(int i=1,x;i<=c;i++)
    {
        x=M[i].size()-1;
        if( x<0 || !L[i].size() )   continue ;
        for(int j=0;j<=x;j++)   add( L[i][0],M[i][j] );
    }
    for(int i=1;i<=n;i++)
    {
        if( type[i]!=3 )   continue ;
        for(int k=1;k<=8;k++)
        {
            ex=p[i].x+dx[k],ey=p[i].y+dy[k];
            if( ex<1 || ex>r || ey<1 || ey>c )   continue ;
            l=P.find( (point){ ex,ey } );
            if( l==P.end() )   continue ;
            add( i,id[c*(ex-1)+ey] );
        }
    }
    num=0;
    for(int i=1;i<=n;i++)   if( !dfn[i] )   Tarjan( i );
    for(int i=1;i<=tot;i++)   dfn[i]=0;
    cnt=0;
    for(int i=1;i<=tot;i++)
        for(int j=Head[i]; j ;j=e[j].next)   if( be[e[j].to]!=i )   deg[be[e[j].to]]++;
    topology();
    for(int i=1;i<=tot;i++)   dis[i]=siz[i];
    for(int i=1,x;i<=tot;i++)
    {
        x=dfnm[i];
        for(int j=Head[x]; j ;j=e[j].next)
            if( be[e[j].to]!=x )
                dis[be[e[j].to]]=max( dis[be[e[j].to]],dis[x]+siz[be[e[j].to]] );
    }
    for(int i=1;i<=tot;i++)   ans=max( ans,dis[i] );
    printf("%d\n",ans);
    return 0;
}
清疚
关注
专栏目录
bzoj1924 所驼门王的宝藏 有向图最长链
lych的博客
02-09 1718
思路还是简单的，代码就。。。。反正bzoj上比我快的都比我长（3k+）     显然，如果把每个点向它能到达的点连一条有向边，那么答案就是这个图的最长链（废话）。然后就是一个景点的强连通分量缩点+dp最长路。但这道题目不能直接暴力连，因为对于同一行的横天门，如果直接暴力连的话是O(N^2)的。实际上同一行的横天门，只需要连向一个切连双向边即可；另外，同一行的横天门如果直接向该行的其他种类的门连边
bzoj2190 SDOI2008仪仗队
老臣
04-10 327
分析：我觉得这题肯定有公式，但是没推出来。。 正解是求解欧拉函数，因为所有能被看见的点横纵坐标必定互质。。#include<iostream> #include<cstring> #include<cstdio> using namespace std; int N,ans; int p[80000],phi[80000],prime[80000]; inline void getphi(){
[SDOI2010]所驼门王的宝藏 --tarjan缩点+最长路
weixin_30718391的博客
08-24 132
[SDOI2010]所驼门王的宝藏 题目描述 在宽广的非洲荒漠中，生活着一群勤劳勇敢的羊驼家族。被族人恭称为“先知”的Alpaca L. Sotomon是这个家族的领袖，外人也称其为“所驼门王”。所驼门王毕生致力于维护家族的安定与和谐，他曾亲自率军粉碎河蟹帝国主义的野蛮侵略，为族人立下赫赫战功。所驼门王一生财宝无数，但...
[SDOI2010] 所驼门王的宝藏
weixin_30546933的博客
06-01 101
Description 在宽广的非洲荒漠中，生活着一群勤劳勇敢的羊驼家族。被族人恭称为“先知”的Alpaca L. Sotomon是这个家族的领袖，外人也称其为“所驼门王”。所驼门王毕生致力于维护家族的安定与和谐，他曾亲自率军粉碎河蟹帝国主义的野蛮侵略，为族人立下赫赫战功。所驼门王一生财宝无数，但因其生性节俭低调，他将财宝埋藏在自己设计的地下宫殿里，这也是今天Henry Curtis故事的起点。H...
BZOJ1924: [Sdoi2010]所驼门王的宝藏
wxh010910
12-30 462
Description Input 第一行给出三个正整数 N, R, C。 以下 N 行，每行给出一扇传送门的信息，包含三个正整数xi, yi, Ti，表示该传送门设在位于第 xi行第yi列的藏宝宫室，类型为 Ti。Ti是一个1~3间的整数， 1表示可以传送到第 xi行任意一列的“横天门”，2表示可以传送到任意一行第 yi列的“纵寰门”，3表示可以传送到周围 8格宫室的“自由门
【洛谷2403】[SDOI2010] 所驼门王的宝藏（Tarjan+dfs遍历）
CXR的博客
07-28 568
点此看题面 大致题意：一个由R∗CR∗CR*C间矩形宫室组成的宫殿中的NNN间宫室里埋藏着宝藏。由一间宫室到达另一间宫室只能通过传送门，且只有埋有宝藏的宫室才有传送门。传送门分为3种，分别可以到达同行的任一宫室（横天门）、同列的任一宫室（纵寰门）和以该宫室为中心周围8个的任一宫室（自由门）。现在你可以从任一宫室开始寻宝，并可以在任一宫室结束寻宝，请求出最多可获得的宝藏数目（每个宝藏只能获得一次）...
[bzoj1924][强联通][最长路]所驼门王的宝藏
Rose_max的博客
02-09 411
Description 整座宫殿呈矩阵状,由 R×C 间矩形宫室组成,其中有 N 间宫室里埋藏着宝藏,称作藏宝宫室。宫殿里外、相邻宫室间都由坚硬的实体墙阻隔,由一间宫室到达另一间只能通过所驼门王独创的移动方式——传送门。所驼门王为这 N 间藏宝宫室每间都架设了一扇传送门,没有宝藏的宫室不设传送门,所有的宫室传送门分为三种: “横天门”:由该门可以传送到同行的任一宫室;...
BZOJ 1925 [Sdoi2010]地精部落
无尽
09-14 477
抖动序列DP思维难度好大- -，我就发一个比较详细的题解吧记f[i][j]表示长度为i（可以理解为序列是1~j的排列），首项取值为[1,j]的第一段为下降抖动序列方案数。（这个方程还真是很难想出来，要考虑到长度为i的序列可以由长度为i-1的序列通过首项的关系转移过来）方程：f[i][j]=f[i][j-1]+f[i-1][i-j]f[i][j-1]好理解，这可以直接求出首项为[1,j-1]的方案数，
BZOJ 3994 Sdoi2015 约数个数和 莫比乌斯反演
世界
04-16 4791
题目大意：求∑ni=1∑mj=1d(ij)\sum_{i=1}^n\sum_{j=1}^md(ij) 首先我们有一个很神的结论： ∑ni=1∑mj=1d(ij)=∑ni=1∑mj=1⌊ni⌋⌊mj⌋[gcd(i,j)==1]\sum_{i=1}^n\sum_{j=1}^md(ij)=\sum_{i=1}^n\sum_{j=1}^m\lfloor\frac n i\rfloor\lfloor\fr
BZOJ4518 Sdoi2016 征途 【斜率优化DP】 *
Dream_maker的博客
09-01 234
BZOJ4518 Sdoi2016 征途 Description Pine开始了从S地到T地的征途。 从S地到T地的路可以划分成n段，相邻两段路的分界点设有休息站。 Pine计划用m天到达T地。除第m天外，每一天晚上Pine都必须在休息站过夜。所以，一段路必须在同一天中走完。 Pine希望每一天走的路长度尽可能相近，所以他希望每一天走的路的长度的方差尽可能小。 帮助Pine求出...
【SDOI2010】【BZOJ1924】所驼门王的宝藏
CreationAugust is 14 years old forever
09-22 1933
Description Input第一行给出三个正整数 N, R, C。 以下 N 行，每行给出一扇传送门的信息，包含三个正整数xi, yi, Ti，表示该传送门设在位于第 xi行第yi列的藏宝宫室，类型为 Ti。Ti是一个1~3间的整数， 1表示可以传送到第 xi行任意一列的“横天门”，2表示可以传送到任意一行第 yi列的“纵寰门”，3表示可以传送到周围 8格宫室的“自由门”。 保证 1≤xi≤
bzoj1924: [Sdoi2010]所驼门王的宝藏
zyy大蒟蒻的博客
07-12 450
传送门 我们首先可以用O（NlogN）的时间建出来所有的边。 然后我们大力tarjan求联通块 之后大力拓扑排序一发就可以了。#include<cstdio> #include<iostream> #include<cstdlib> #include<cmath> #include<cstring> #include<algorithm> #include<vector> #include<
洛谷P2403 [SDOI2010]所驼门王的宝藏
tylon2006的博客
08-10 182
Description 在宽广的非洲荒漠中,生活着一群勤劳勇敢的羊驼家族。被族人恭称为“先知”的 Alpaca L. Sotomon 是这个家族的领袖,外人也称其为“所驼门王”。所驼门王毕生致力于维护家族的安定与和谐,他曾亲自率军粉碎河蟹帝国主义的野蛮侵略,为族人立下赫赫战功。所驼门王一生财宝无数,但因其生性节俭低调, 他将财宝埋藏在自己设计的地下宫殿里,这也是今天 Henry Curtis 故事的起点。Henry 是一个爱财如命的贪婪家伙,而又非常聪明,他费尽心机谋划了这次盗窃行动,破解重重机关后来到这座
【SDOI2010】所驼门王的宝藏
weixin_30362233的博客
07-31 104
题面 https://www.luogu.org/problem/P2403 题解 我写的是让同一行、同一列的点形成一个环，表示它们互相可达。 $aysn$写的是对于同一行列建一个新点表示互相可达。 看了窝还是太弱了，还须多练。 #include<cstdio> #include<iostream> #include<cstring&gt...
bzoj1924 [Sdoi2010]所驼门王的宝藏
wumingshi
05-28 491
目标在我们遥不可及的地方，我们难道就只能翘首仰望么？
BZOJ 1924 所驼门王的宝藏【思维建图+强连通缩点+DAG图上做最长路】这题好尼玛综合........
mengxiang000000
02-08 631
1924: [Sdoi2010]所驼门王的宝藏 Time Limit: 5 Sec  Memory Limit: 128 MB Submit: 913  Solved: 379 [Submit][Status][Discuss] Description Input 第一行给出三个正整数 N, R, C。 以下 N 行，每行给出一扇传送门的信息，包含三个正整数xi, yi
1924: [Sdoi2010]所驼门王的宝藏
OZY的博客
09-09 315
dfs+缩点
【题解】SDOI2010所驼门王的宝藏(强连通分量+优化建图)
weixin_30482383的博客
08-09 118
【题解】SDOI2010所驼门王的宝藏(强连通分量+优化建图) 最开始我想写线段树优化建图的说，数据结构学傻了233 虽然矩阵很大，但是没什么用，真正有用的是那些关键点 考虑关键点的类型： 横走型 竖走型 八连通型 本质上只有两种类型(走一大串/走八连通)，我们考虑这样一种建图方法： 对于每一行每一列建立一个点(点权为\(0\)) 对于关键点建立一个点(点权为\(1\)) 然后考...
[SDOI2010]BZOJ 1924所驼门王的宝藏-强连通分量-缩点-拓扑排序-dp
Mys_C_K的博客
04-25 511
题目链接：右转进入题目 题目大意：自行参考题目 题解：tarjan缩点后拓排一下dp乱搞即可。 代码：//BZOJ 1924 //SDOI 2010 #include #include #include #include #include #include #include #define MAXN 100010 #define MAXRC 1000010 using namespace s
洛谷 P1640 BZOJ 1854 [SCOI2010]连续攻击游戏
最新发布
05-14
这道题目需要进行博弈论的思考。我们可以使用动态规划来解决这道题目。 首先，我们可以定义 $f_{i,j}$ 表示当前剩下 $i$ 到 $j$ 这些数时，当前玩家最多能获得的分数。我们可以发现，在这个状态中，当前玩家只能从 ...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值