[动规] 旅行计划

最新推荐文章于 2024-08-18 08:31:43 发布

莓尊尊strawberry

最新推荐文章于 2024-08-18 08:31:43 发布

阅读量171

点赞数

分类专栏： -.- 散题 -.- 文章标签：算法图论动态规划

本文链接：https://blog.csdn.net/apigisyou/article/details/126515334

版权

-.- 散题 -.- 专栏收录该内容

22 篇文章 0 订阅

订阅专栏

旅行计划

题目描述

小明要去一个国家旅游。这个国家有# $N$ 个城市，编号为 $1$ 至 $N$ ，并且有 $M$ 条道路连接着，小明准备从其中一个城市出发，并只往东走到城市i停止。

所以他就需要选择最先到达的城市，并制定一条路线以城市i为终点，使得线路上除了第一个城市，每个城市都在路线前一个城市东面，并且满足这个前提下还希望游览的城市尽量多。

现在，你只知道每一条道路所连接的两个城市的相对位置关系，但并不知道所有城市具体的位置。现在对于所有的i，都需要你为小明制定一条路线，并求出以城市 $i$ 为终点最多能够游览多少个城市。

点击查看完整题目信息。

题解

状态转移方程： $d p (v) = ma x (d p (v), d p (u) + 1)$ 。

$d p (v)$ 表示：节点 v 能游览的最多城市，u 是一条边的起点，v 是一条边的终点。

题中的图一定是有向无环图，因为方向取决于小明行走的方向，小明总是向东，所以每两个点之间的边的方向都是向右（东）的。

代码（cpp）：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<cmath>
using namespace std;

int N, M;
int first[100001], curAdjEdge, in_deg[100001]; // 邻接表相关
int sortedPoint[100001], sortedCount; // 拓扑排序相关
int dp[100001];

struct EDGE // 相当于邻接表的 next 数组
{
  int end; // 线段终点
  int next; // 指向同一起点的下一条线段，nextEdge = edge[curEdge.next]
} edge[200001];

// 生成邻接表
void genAdjList(int start, int end)
{
  ++ curAdjEdge;
  edge[curAdjEdge].next = first[start];
  edge[curAdjEdge].end = end;
  first[start] = curAdjEdge;
}

// 拓扑排序
void topSort()
{
  int u, v; // 起点，终点
  int i;
  queue<int> q;
  for (i = 0; i < N; ++ i)
    if (in_deg[i] == 0)
      q.push(i);

  while(! q.empty())
  {
    u = q.front();
    q.pop();
    sortedPoint[++ sortedCount] = u;
    for (i = first[u]; i; i = edge[i].next) // 遍历邻接表，u 为起点的所有路线
    {
      v = edge[i].end;
      -- in_deg[v]; // 入度 -1
      if (in_deg[v] == 0)
        q.push(v);
    }
  }
}

int main()
{
  int u, v; // 起点，终点
  int i, j;
  scanf("%d %d", &N, &M);
  for (i = 0; i < M; ++ i)
  {
    scanf("%d %d", &u, &v);
    genAdjList(u, v);
    ++ in_deg[v]; // 入度 +1
  }
  topSort();
  for (i = 1; i < N + 1; ++ i) 
    dp[i] = 1;
  for (i = 1; i < N + 1; ++ i)
  {
    u = sortedPoint[i];
    for (j = first[u]; j; j = edge[j].next) // 遍历邻接表，u 为起点的所有路线
    {
      v = edge[j].end;
      dp[v] = max(dp[v], dp[u] + 1);
    }
  }
  for (i = 1; i < N + 1; ++ i)
    printf("%d\n", dp[i]);
  return 0;
}

拓扑排序

拓扑排序（Topological sorting）排的是有向无环图（DAG）的节点的顺序，例如有图1>2>3<4<5，排过序的节点就是1>5>2>4>3。同时也可以用拓扑排序判断图是有环还是无环的，如果有图1>2>3>1，就形成了环，这个图不能进行拓扑排序。

例如，“做完旅行计划”的前提是“掌握拓扑排序”，这是如果说“掌握拓扑排序”的前提是“做完旅行计划”，就形成了一个环，不能计算拓扑排序；但是如果只说“做完旅行计划”的前提是“掌握拓扑排序”，拓扑排序的结果就是：“掌握拓扑排序”>“做完旅行计划”。

伪代码（Kahn 算法）：

L ← 清空以存放排过序的节点
S ← 初始化，推入图中入度为 0 的节点
while S 不为空 do
    弹出 S 的一个节点 n
    向 L 推入节点 n
    for each 节点 m（m 是所有以节点 n 为起点的边 e 的终点） do // 遍历
        从图中移除边 e
        if 节点 m 入度为 0 then
            向 S 推入节点 m
if 图中有边 then
    return error (图是有环的)
else
    return 拓扑排序过的队列 L

时间复杂度：对于图 G = (V, E)，复杂度为 O(V + E)。V 表示点，E 表示边。

题解代码里的topSort函数是这段伪代码的一种实现。