离散数学:过河问题

过河问题(一)图论方法

农夫过河问题

问题描述

一个农夫带着一头狼、一头羊、一颗白菜过河。他面前只有一条船,只能容纳他和一件物品,只有农夫会划船。如果农夫不在场,狼会吃羊、羊会吃白菜,农夫在场则不会。求将所有物品运到对岸的方案。

解题思路

根据物品的位置定义状态,若在左岸记为1,右岸记为0,于是最终方案就是(1,1,1,1)-->(0,0,0,0)所经过的路径。

1、定义状态

2、列举所有状态(人、狼、羊、菜)

3、删除不合理的状态(狼和羊、羊和菜)

4、连边(模拟一次渡河)

5、寻找路径

寻找(1111)-->(0000)的边,可以用寻路算法如bfs、dfs,如果要求最短路可以用最短路算法如bfs、Dijsktra等,当然这里图很简单,可直接观察出来。

复制代码

(1111)-->(0101)-->(1101)-->(0001)-->(1011)-->(0010)-->(1010)-->(0000)(两条最短路之一)
       左岸                                 右岸
1、人  狼  羊 花                           空
2、狼  花                                 人 羊
3、人  狼 花                               羊
4、花                                     人 狼  羊
5、人  羊  花                              狼
6、羊                                     人 花  狼
7、人 羊                                  狼 花
8、空                                     狼 花 人 羊   

复制代码

传教士与吃人恶魔的问题

问题描述

有三个传教士和三个吃人恶魔要渡过一条河,河中有一条船,只能装下两个人。在任何地方(无论是岸边还是船上),如果吃人恶魔数量多于传教士数量,吃人恶魔就会吃掉传教士。问:怎么才能让这些都安全过河?

解题思路

1、定义状态

2、列举所有状态

3、删除不合理状态

4、连边(模拟依次渡河变化)

5、寻找路径

寻找(33 L 00)-->(00 R 33)的路径

复制代码

其中一条路径
(33 L 00)-->(31 R 01)-->(32 L 01)-->(30 R 03)-->(31 L 02)-->(11 R 22)-->(22 L 01)-->(02 R 31)-->(03 L 30)-->(01 R 32)-->(02 L 31)-->(00 R 33)
1、两个吃人恶魔过河
2、一个吃人恶魔回来
3、两个吃人恶魔过河
4、一个吃人恶魔回来
5、两个传教士过河
6、一个传教士和一个吃人恶魔回来
7、两个传教士回来
8、一个吃人恶魔回去
9、两个吃人恶魔过河
10、一个吃人恶魔回去
11、两个吃人恶魔过河

复制代码

四人过桥问题

问题描述

在一个漆黑的夜里,四位旅游者来到一座狭窄而没有护栏的桥边,如果不借助手电筒的话,大家是无论也不敢过去。不幸的是四个人中只有一只手电筒,而桥窄得只够两个人同时通过。如果各自单独过桥得话,四个人所需要的时间分别是1、2、5、10分钟,如果两个人同时过桥,所需要的时间是较慢的那个人单独行动时的时间。问:如何设计一个方案,让四个人尽快过桥。

解题思路

与前面两个相比,这次不仅要求方案,同时要求时间最短。 

同样需要定义状态,四个人+手电筒的位置

1、定义状态

2、建图

分为每次通过一个人和每次两个人,都是带权无向边。

(下面只连接了与(01111)的边)

3、寻找最短路

寻找(L 1111)-->(R 0000)的最短路,即最短路算法中(01111)-->(10000)的最短路,以下是利用Dijstra算法的解决方法。

最终答案为(2 + 1 + 10 + 2 + 2) = 17.

复制代码

  1 #include<stdio.h>
  2 #include<iostream>
  3 #include<string>
  4 #include<queue>
  5 #include<algorithm>
  6 using namespace std;
  7 
  8 //定义图中结点
  9 struct Node
 10 {
 11     int u, d;            //该节点的编号与距离
 12     bool operator < (const Node x) const
 13     {
 14         return  d > x.d;
 15     }
 16 };
 17 
 18 //边结构体的定义
 19 struct Edge
 20 {
 21     int to;     
 22     int w;       
 23     int next;  
 24 };
 25 
 26 const int INF = 0x3f3f3f3f;
 27 const int V = 32 + 10;
 28 const int E = 32 * 32 + 10;
 29 int dis[V];            //源到各顶点的最短距离
 30 int vis[V];            //记录是否被收录,用来代替集合S
 31 int head[V];          //head[i]表示顶点i的第一条边的数组下标,"-1"表示顶点i没有边
 32 Edge edge[E];
 33 
 34 inline void AddEdge(int a, int b, int w, int id)
 35 {
 36     edge[id].to = b;
 37     edge[id].w = w;
 38     edge[id].next = head[a];  
 39     head[a] = id;             
 40     return;
 41 }
 42 
 43 //s为起点
 44 void Dijsktra(int s)
 45 {
 46     priority_queue<Node>q;            //取出集合T中的最小值
 47     memset(vis, 0, sizeof(vis));
 48     memset(dis, INF, sizeof(dis));    
 49 
 50     dis[s] = 0;
 51     q.push(Node{ s, dis[s] });
 52     while (!q.empty())
 53     {
 54         Node x = q.top(); q.pop();
 55         int u = x.u;
 56 
 57         if (vis[u])    continue;
 58 
 59         vis[u] = true;
 60         for (int i = head[u]; i != -1; i = edge[i].next)    //松弛与u直接相邻的顶点
 61         {
 62             int v = edge[i].to;
 63             int w = edge[i].w;
 64             if (!vis[v] && dis[u] + w < dis[v])
 65             {
 66                 dis[v] = dis[u] + w;
 67                 q.push(Node{ v,dis[v] });
 68             }
 69         }
 70     }
 71 }
 72 
 73 const int score[] = { 1,2,5,10 };   //每个人单独行动的时间
 74 
 75 int main()
 76 {
 77     //建图
 78     memset(head, -1, sizeof(head));
 79     int id = 0;
 80     for (int i = 0; i < (1 << 4); i++)
 81     {
 82         int bits[4];
 83         for (int j = 0; j < 4; j++)  bits[j] = (i >> j) & 1;
 84         //一次走一个人
 85         for (int j = 0; j < 4; j++)  if (bits[j])
 86         {
 87             int tmp = i - (1 << j) + 16;
 88             int w = score[j];
 89             AddEdge(i, tmp, w, id++);
 90             AddEdge(tmp, i, w, id++);
 91         }
 92         //一次走两个人
 93         for(int j = 0;j < 3;j++)
 94             for (int k = j + 1; k < 4; k++)   if (bits[j] && bits[k])
 95             {
 96                 int tmp = i - (1 << j) - (1 << k) + 16;
 97                 int w = max(score[j],score[k]);
 98                 AddEdge(i, tmp, w, id++);
 99                 AddEdge(tmp, i, w, id++);
100             }
101     }
102     Dijsktra(15);
103     printf("%d\n", dis[16]);
104 
105     return 0;
106 }

复制代码

此类问题很多,但大多可用图论的思想做(虽然不一定是速度最快的),后续在补充吧,有问题直接留言!

 

过河问题(二)

问题描述

POJ 1700

在一个漆黑的夜里,N为旅行者来到狭窄而没有护栏的桥边。如果不借助手电筒的话,大家是无论如何也不敢过桥去的。不幸的是,N个人一共只带了一只手电筒,而桥窄得只够让两个人同时过。如果各自单独过桥的话,N人所需要的时间已知;而如果两人同时过桥,所需要的时间就是走得比较慢的那个人单独行动时所需的时间。问题是,如何设计一个方案,让这N人尽快过桥。

解题思路

很容易想到使用贪心策略,但具体怎么做呢?

假设共5个人T1、T2、T3、T4、T5,时间一次增加,如何用最少的时间使得5人都过河?

考虑先将T4和T5过河(最慢的两个安排好了,无后效性),经过简单比较,有两种较优方法。

1、T1和T2过河,T1回来,T4和T5过河,T2回来(总时间:T2 + T1 + T5 + T2)

2、T1和T4过河,T1回来,T1和T5过河,T1回来(总时间:T4 + T1 + T5 + T1)

然后问题简化为只有3个人T1、T2、T3,n≤3的情况直接观察。

推广到n个人(n ≥ 4),只需取(T1 + 2*T2 + Tn)和(2*T1 + Tn-1 + Tn)中的较小者,问题化为前n-2个人的情况。

代码实现

复制代码

 1 #include<stdio.h>
 2 #include<iostream>
 3 #include<algorithm>
 4 using namespace std;
 5 
 6 const int maxn = 1000 + 10;
 7 int n,times[maxn];
 8 
 9 int main()
10 {
11     int T;
12     scanf("%d", &T);
13     while (T--)
14     {
15         scanf("%d", &n);
16         for (int i = 0; i < n; i++)  scanf("%d", &times[i]);
17         sort(times, times + n);
18         int sum = 0;
19         while (n > 3)
20         {
21             sum += min((times[0] + times[1] * 2 + times[n - 1]), (times[0] * 2 + times[n - 2] + times[n - 1]));
22             n -= 2;
23         }
24         if (n == 3)  sum += (times[0] + times[1] + times[2]);
25         if (n == 2)  sum += times[1];
26         if (n == 1)  sum += times[0];
27 
28         printf("%d\n", sum);
29     }
30     return 0;
31 }

复制代码

 参考链接:

https://www.cnblogs.com/ShiChaoPeng/p/5858690.html

https://blog.csdn.net/qq_40788630/article/details/79332774?utm_source=blogxgwz1

参考链接:中国大学mooc 刘铎老师  离散数学

  • 4
    点赞
  • 41
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
商人过河问题可以使用数学建模来解决,而MATLAB是一个非常适合进行数学建模和计算的工具。下面是一个使用MATLAB进行商人过河问题数学建模的示例: 首先,我们可以使用二进制向量来表示商人、野人和船的位置状态,其中0表示左岸,1表示右岸。假设商人和野人的数量分别为M和N。 下面是MATLAB代码示例: % 商人和野人的数量 M = 3; N = 3; % 初始状态(左岸) initial_state = [ones(1, M+N), 0]; % 目标状态(右岸) target_state = [zeros(1, M+N), 1]; % 状态转移规则函数 state_transition = @(state, action) state + action; % 判断状态是否合法 is_valid_state = @(state) ... all(state(1:M) >= state(M+1:end) | state(1:M) == 0) && ... all(state(M+1:end) >= state(1:M) | state(M+1:end) == 0); % 判断是否达到目标状态 is_goal_state = @(state) all(state == target_state); % 使用递归函数解决问题 solution = recursive_solve(initial_state, []); % 递归函数 function solution = recursive_solve(state, path) % 如果已经达到目标状态,则返回路径 if is_goal_state(state) solution = path; return; end % 遍历所有可能的动作 actions = generate_actions(state); for i = 1:size(actions, 1) action = actions(i, :); % 计算新状态 new_state = state_transition(state, action); % 如果新状态合法,则继续递归求解 if is_valid_state(new_state) solution = recursive_solve(new_state, [path; action]); % 如果找到解,则返回 if ~isempty(solution) return; end end end % 如果找不到解,则返回空 solution = []; end % 生成所有可能的动作 function actions = generate_actions(state) M = sum(state(1:end-1) == 1); N = sum(state(1:end-1) == 0); % 动作格式:[商人移动数 野人移动数 船移动方向] % 商人独自移动 actions = [-1 0 -1; -2 0 -1; 1 0 1; 2 0 1]; % 野人独自移动 actions = [actions; 0 -1 -1; 0 -2 -1; 0 1 1; 0 2 1]; % 商人和野人一起移动 for i = 1:M for j = 1:N if i + j <= 2 actions = [actions; -i -j -1; i j 1]; end end end end 这个示例代码使用了递归求解的方法来找到商人过河问题的解。代码中的state_transition函数定义了状态转移规则,is_valid_state函数判断状态是否合法,is_goal_state函数判断是否达到目标状态。generate_actions函数生成所有可能的动作。 注意:这只是商人过河问题的一个简单数学建模示例,实际问题可能需要更多的约束和复杂的规则。你可以根据具体需求进行修改和扩展。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值