欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页

程序设计思维与实践 Week7 作业 (3/4/数据班)

程序员文章站 2022-03-18 13:41:19
...

A - TT 的魔法猫

题意:
众所周知,TT 有一只魔法猫。

这一天,TT 正在专心致志地玩《猫和老鼠》游戏,然而比赛还没开始,聪明的魔法猫便告诉了 TT 比赛的最终结果。TT 非常诧异,不仅诧异于他的小猫咪居然会说话,更诧异于这可爱的小不点为何有如此魔力?

魔法猫告诉 TT,它其实拥有一张游戏胜负表,上面有 N 个人以及 M 个胜负关系,每个胜负关系为 A B,表示 A 能胜过 B,且胜负关系具有传递性。即 A 胜过 B,B 胜过 C,则 A 也能胜过 C。

TT 不相信他的小猫咪什么比赛都能预测,因此他想知道有多少对选手的胜负无法预先得知,你能帮帮他吗?
思路:
分析题意,胜负关系具有传递性,所以使用floyed来做。

  • 其中用dis数组表示胜负关系,dis[a][b]=1时表示a强于b,dis[a][b]=0时表示胜负关系不知道。
  • 松弛条件:dis[i][j]=max(dis[i][j],dis[i][k]&dis[k][j]);

总结:
我第一次做的时候,没有剪枝,TE了,分析之后可以发现,当i与k的胜负关系未知时,
即dis[i][k]==0时,就不用考虑k和j的胜负关系,可以降低复杂度。
代码:

#include <iostream>
using namespace std;
const int N = 500 + 1;
int dis[N][N];
int main()
{
    int num;
    cin >> num;
    for (int c = 0; c < num; c++)
    {
        int n, m;
        cin >> n >> m;
        for (int j = 0; j < m; j++)
        {
            int x, y;
            cin >> x >> y;
            dis[x][y] = 1;
        }
        for (int k = 1; k <= n; k++)
            for (int i = 1; i <= n; i++)
                if (dis[i][k] != 0)
                    for (int j = 1; j <= n; j++)
                        dis[i][j] = max(dis[i][j], dis[i][k] & dis[k][j]);
        int ans = 0;
        for (int i = 1; i <= n; i++)
            for (int j = i + 1; j <= n; j++)
                if (dis[i][j] == 0 && dis[j][i] == 0)
                    ans++;
        cout << ans << endl;
        for (int i = 1; i <= n; i++)
            for (int j = 1; j <= n; j++)
                dis[i][j] = 0;
    }
    //system("pause");
}

B - TT 的旅行日记

题意:
众所周知,TT 有一只魔法猫。

今天他在 B 站上开启了一次旅行直播,记录他与魔法猫在喵星旅游时的奇遇。 TT 从家里出发,准备乘坐猫猫快线前往喵星机场。猫猫快线分为经济线和商业线两种,它们的速度与价钱都不同。当然啦,商业线要比经济线贵,TT 平常只能坐经济线,但是今天 TT 的魔法猫变出了一张商业线车票,可以坐一站商业线。假设 TT 换乘的时间忽略不计,请你帮 TT 找到一条去喵星机场最快的线路,不然就要误机了!
思路:
学长讲了两种思路,这里主要说一下第一种。
题意要求如果走商业线的话,只能走一条。如果商业线为(a,b),则我们可能的线路为
起点->a->b->终点 或 起点->终点 或 起点->b->a->终点。所以我们的方法为:

  1. 从起点求dijkstra最短路,记录dis1数组。
  2. 从终点求dijkstra最短路,记录dis2数组。
  3. 使用min求取最短路。

代码:

#include <iostream>
#include <cstring>
#include <climits>
#include <queue>
using namespace std;
const int N = 100000 + 1;
int vis[N], dis1[N], dis2[N], head[N], path1[N], path2[N], pathlong = INT_MAX;
int pos1, pos2;
bool flage = false;
vector<int> v;
priority_queue<pair<int, int>> q;
int n, s, e, tot = 0;
struct Edge
{
    int v, w, nxt;
} Edges[2 * N];
void add(int x, int y, int z)
{
    Edges[++tot].v = y, Edges[tot].w = z;
    Edges[tot].nxt = head[x];
    head[x] = tot;
}
void dijkstra(int s, int *dis, int *path)
{
    while (q.size())
        q.pop();
    dis[s] = 0;
    q.push(make_pair(0, s));
    while (q.size())
    {
        int x = q.top().second;
        q.pop();
        if (vis[x])
            continue;
        vis[x] = 1;
        for (int i = head[x]; i; i = Edges[i].nxt)
        {
            int y = Edges[i].v, w = Edges[i].w;
            if (dis[y] > dis[x] + w)
            {
                dis[y] = dis[x] + w;
                q.push(make_pair(-dis[y], y));
                path[y] = x;
            }
        }
    }
}
int main()
{
    while (scanf("%d%d%d", &n, &s, &e) == 3)
    {
        if (!flage)
            flage = true;
        else
            cout << endl;
        memset(head, 0, sizeof(head));
        memset(vis, 0, sizeof(vis));
        memset(path1, 0, sizeof(path1));
        memset(path2, 0, sizeof(path2));
        memset(dis1, 999999, sizeof(dis1));
        memset(dis2, 999999, sizeof(dis2));
        int m;
        cin >> m;
        for (int i = 0; i < m; i++)
        {
            int x, y, z;
            cin >> x >> y >> z;
            add(x, y, z);
            add(y, x, z);
        }
        dijkstra(s, dis1, path1);
        memset(vis, 0, sizeof(vis));
        dijkstra(e, dis2, path2);
        int k, ans = INT_MAX;
        cin >> k;
        for (int i = 0; i < k; i++)
        {
            int x, y, z;
            cin >> x >> y >> z;
            if (dis1[x] + dis2[y] + z < pathlong)
            {
                pos1 = x;
                pos2 = y;
                pathlong = dis1[x] + dis2[y] + z;
                ans = x;
            }
            if (dis2[x] + dis1[y] + z < pathlong)
            {
                pos1 = y;
                pos2 = x;
                pathlong = dis2[x] + dis1[y] + z;
                ans = y;
            }
        }
        if (dis1[e] < pathlong)
        {
            for (int i = e; i != s; i = path1[i])
                v.push_back(i);
            v.push_back(s);
            for (int i = v.size() - 1; i >= 1; i--)
                cout << v[i] << " ";
            cout << v[0] << endl
                 << "Ticket Not Used" << endl
                 << dis1[e] << endl;
        }
        else
        {
            for (int i = pos1; i != s; i = path1[i])
                v.push_back(i);
            v.push_back(s);
            for (int i = pos2; i != e; i = path2[i])
                v.insert(v.begin(), i);
            v.insert(v.begin(), e);
            for (int i = v.size() - 1; i >= 1; i--)
                cout << v[i] << " ";
            cout << v[0] << endl
                 << ans << endl
                 << pathlong << endl;
        }
        v.clear();
        ans = 0;
        tot = 0;
        pathlong = INT_MAX;
    }
    //system("pause");
}

C - TT 的美梦

题意:
这一晚,TT 做了个美梦!

在梦中,TT 的愿望成真了,他成为了喵星的统领!喵星上有 N 个商业城市,编号 1 ~ N,其中 1 号城市是 TT 所在的城市,即首都。

喵星上共有 M 条有向道路供商业城市相互往来。但是随着喵星商业的日渐繁荣,有些道路变得非常拥挤。正在 TT 为之苦恼之时,他的魔法小猫咪提出了一个解决方案!TT 欣然接受并针对该方案颁布了一项新的政策。

具体政策如下:对每一个商业城市标记一个正整数,表示其繁荣程度,当每一只喵沿道路从一个商业城市走到另一个商业城市时,TT 都会收取它们(目的地繁荣程度 - 出发地繁荣程度)^ 3 的税。

TT 打算测试一下这项政策是否合理,因此他想知道从首都出发,走到其他城市至少要交多少的税,如果总金额小于 3 或者无法到达请悄咪咪地打出 ‘?’。
思路:
分析题意后,可以知道这个题是含有负权边的有向图单源最短路问题。
当出现负环时,交税为负数,需要着重考虑。
使用spfa顺利解决,具体步骤为:

  1. 从一个城市开始搜索
  2. 到达某一个城市经过的边数超过n-1,则出现负环
  3. 从该城市搜索,到达的路径都标记上,不符合题意。

代码:

#include <iostream>
#include <cstring>
#include <climits>
#include <math.h>
#include <queue>
using namespace std;
const int N = 100000 + 5;
const int M = 200 + 5;
int vis[M], dis[M], cnt[M], head[N], a[N];
bool flag[N];
int n, tot = 0;
struct Edge
{
    int v, w, nxt;
} Edges[N];
void add(int x, int y, int z)
{
    Edges[++tot].v = y, Edges[tot].w = z;
    Edges[tot].nxt = head[x];
    head[x] = tot;
}
void bfs(int s)
{
    queue<int> q;
    while (q.size())
        q.pop();
    q.push(s);
    flag[s] = true;
    while (!q.empty())
    {
        int x = q.front();
        q.pop();
        for (int i = head[x]; i; i = Edges[i].nxt)
        {
            int y = Edges[i].v;
            if (!flag[y])
            {
                q.push(y);
                flag[y] = true;
            }
        }
    }
}
void spfa(int s)
{
    queue<int> q;
    while (q.size())
        q.pop();
    for (int i = 1; i <= n; i++)
    {
        vis[i] = cnt[i] = 0;
        dis[i] = 99999;
        flag[i] = false;
    }
    dis[s] = 0;
    vis[s] = 1;
    q.push(s);
    while (!q.empty())
    {
        int x = q.front();
        q.pop();
        vis[x] = 0;
        for (int i = head[x]; i; i = Edges[i].nxt)
        {
            int y = Edges[i].v;
            if (dis[y] > dis[x] + Edges[i].w)
            {
                cnt[y] = cnt[x] + 1;
                if (cnt[y] >= n)
                    bfs(y);
                dis[y] = dis[x] + Edges[i].w;
                if (vis[y] == 0 && !flag[y])
                {
                    q.push(y);
                    vis[y] = 1;
                }
            }
        }
    }
}
int main()
{
    
    int t;
    scanf("%d", &t);
    for (int k = 1; k <= t; k++)
    {
        tot = 0;
        memset(a, 0, sizeof(a));
        memset(head, 0, sizeof(head));
        scanf("%d", &n);
        for (int i = 1; i <= n; i++)
            scanf("%d", &a[i]);
        int m;
        scanf("%d", &m);
        for (int i = 1; i <= m; i++)
        {
            int a1, b1;
            scanf("%d%d", &a1, &b1);
            int c1 = pow(a[b1] - a[a1], 3);
            add(a1, b1, c1);
        }
        spfa(1);
        printf("Case %d:\n", k);
        int Q;
        scanf("%d", &Q);
        for (int i = 1; i <= Q; i++)
        {
            int p;
            scanf("%d", &p);
            if (flag[p] || dis[p] < 3 || dis[p] == 99999)
                printf("?\n");
            else
                printf("%d\n", dis[p]);
        }
    }
    //system("pause");
}

总结:
这一周主要学习了图的遍历方法,其中spfa判断负环的方法较难,需要重点掌握。

相关标签: 实验报告