商仆过河问题数学建模文档格式.docx
- 文档编号:19309159
- 上传时间:2023-01-05
- 格式:DOCX
- 页数:16
- 大小:123.19KB
商仆过河问题数学建模文档格式.docx
《商仆过河问题数学建模文档格式.docx》由会员分享,可在线阅读,更多相关《商仆过河问题数学建模文档格式.docx(16页珍藏版)》请在冰豆网上搜索。
可得如下的递归树:
(转下页)
(注解:
当K为奇数时,船在B岸;
当K为偶数时,船在A岸。
)
通过分析该递归树,知道求解关键在于正确地写出基本的状态转移关系式和恰当的边界条件。
因为k为奇数时,船是从A岸驶向B岸,k为偶数时。
船是由B岸驶回A岸。
所以状态SK随决策DK变化的规律是
SK+1=SK+(-1)KDK,k=l,2,?
,
称之为状态转移律,这样,制定过河方案就归结为如下的多步决策问题:
每一步,船由A岸驶向B岸或B岸驶回A岸,都要对船上的人员(商人UK,随从VK各几人)作出决策,在保证安全的前提下即两岸的商人数XK都不比随从数YK少,用有限步使人员全部过河.用状态(变量)SK表示某一岸的人员状况,决策(变量)DK表示船上的人员状况,可以找出状态SK随决策DK变化的规律.这样安全过河问题就转化为:
求决策DK∈D(k=1,2,……,n),使得状态SK∈S,按照状态转移律,由初始状态S1=(3,3),经有限步n到达状态SK+1=(O,O)。
模型建立:
SK+1=SK+(-1)KDK,k=l,2,3,其中DK∈D={(UK,VK)|UK+VK=l,2},{其中SK∈(XK,YK)|(XK=0,YK=1,2,3);
(XK=3,YK=0,1,2,3);
(XK=YK=1,2)},Sn+1=(0,0)
这就是三个商人的过河问题模型。
六、模型求解:
穷举法:
计算机编程(见附)
先建立编程的基本过程,然后考虑模型,再编写程序。
然后就可以得出结果了。
主程序流程图
图解法:
状态s=(x,y)16个格点
允许状态10个●点
允许决策移动1或2格;
k奇,左下移;
k偶,右上移.
总共需要11步
可以得出经过11步的渡河就能达到安全渡河的目标及满足渡河的次数尽量少的条件。
这11步的渡河方案就是上面程序运行结果中船上下面的一列。
八、模型的检验
用2名商人和2名随从的过河问题的解决思路,检验3名商人和3名随从的过河问题。
九、模型的拓展和延伸
通过三名商人和三名随从的过河问题的解决方案,可以进一步计算四名商人和四名随从的过河问题,通过计算机编程可以设计m名商人和n名随从的过河问题。
十、总结
这是通过数学分析的方法解决实用问题,经过问题提出、问题假设、问题分析、模型建立、模型求解、模型检验的过程,解决商人过河问题。
然后扩展延伸到n个商人的问题。
学习数学建模以来,重新认识了学习数学的乐趣,也重新认识了数学,本以为数学是单调的,枯燥的,学习了之后,发现数学是普遍存在我们生活之中的。
解决现实中的问题,很多都需要数学。
沉浸在数学的世界里,发现学习是有趣的;
相比于机械的认识各个组织器官,建立一个数学模型解决问题是十分有趣的。
参考文献:
(1)傅清祥.《数据结构与算法》.王晓东.北京:
电子工业出版社1998.
(2)姜启瑟.《数学建模》(第二版).北京:
高等教育出版社,2000.
(3)运筹学教材编写组.《运筹学》(修订版).北京:
清华大学出版社。
2001.
附:
商仆过河的C程序及运行截屏:
#include<
iostream>
usingnamespacestd;
structNode
{intnMer;
intnSer;
intlength;
};
classA
{
public:
A();
~A();
voidTspt();
//过河的动作
voiddoLeft(intnhead,intntail,intnlength);
private:
boolislegal(intnm,intns);
//判断是否满足约束条件,满足为true
Node*funTspt(intnm,intns,boolflag);
//添加STEP[head]可以向后延伸的节点
boolnoRepeat(intnm,intns);
//没有重复返回TRUE
voidfunshow(inta[][2],intntail);
boolfunLeft(Nodend,intb1,intb2,intn);
voidshow(ints[],intp[][2],int&
top,int&
count,inta[]);
inthead;
inttail;
intn;
//商仆的对数
intnB;
//船最多的载人数目
Node*STEP;
A:
:
~A()
free(STEP);
}
A()
cout<
<
"
请输入商仆的对数S="
;
F:
cin>
>
n;
if(n==1)
{
nB=2;
船最多载人的数目K="
nB;
}
elseif(n==2)
cout<
船最多载人的数目可以取:
for(intx=n;
x<
=2*n;
x++)
{
cout<
、"
}
endl;
请输入船最多载人的数目K="
elseif(n==3)
for(intx=n-1;
elseif(n==4)
elseif(n>
=5&
&
n<
=100)
for(intx=4;
elseif(n<
1||n>
100)
本程序仅在S=(0…100)以内保证其正确性"
请重新输入商仆的对数S="
gotoF;
STEP=(Node*)malloc(sizeof(Node)*10000);
memset(STEP,0,sizeof(Node)*10000);
head=tail=0;
STEP[0].nMer=STEP[0].nSer=n;
intmain()
问题描述:
S个商人各带一个随从乘船过河,一只小船只能容纳K人,由他们自己划船。
Aa;
a.Tspt();
return0;
voidA:
show(ints[],intp[][2],int&
count,inta[])
if(top==-1)
return;
//已找到目标状态需,输出数据
if(top==STEP[head].length)
***********"
++count<
***********"
funshow(p,top+1);
B:
top--;
if(top==-1)
return;
C:
s[top]--;
if(STEP[(s[top])].length!
=top)//退过了
s[top]=a[top];
gotoB;
if(funLeft(STEP[(s[top])],p[top-1][0],p[top-1][1],top-1)==false)
gotoC;
p[top][0]=STEP[(s[top])].nMer;
p[top][1]=STEP[(s[top])].nSer;
show(s,p,top,count,a);
//在中间加入节点STEP[(s[top+1])]
if(funLeft(STEP[(s[top+1])],p[top][0],p[top][1],top)==true)//符合条件
{
top++;
else//不符合条件
E:
s[top+1]--;
if(STEP[(s[top+1])].length==top)//退过了,到了下一层
s[top+1]=a[top+1];
D:
s[top]--;
if(STEP[(s[top])].length!
=top)//退过了,到了下一层
{
for(inti=top;
i<
=STEP[head].length;
i++)
s[i]=a[i];
top--;
if(top==-1)
return;
gotoD;
}
if(top==0)
return;
if(funLeft(STEP[(s[top])],p[top-1][0],p[top-1][1],top-1)==false)
p[top][0]=STEP[(s[top])].nMer;
p[top][1]=STEP[(s[top])].nSer;
show(s,p,top,count,a);
if(funLeft(STEP[(s[top+1])],p[top][0],p[top][1],top)==false)
gotoE;
}
doLeft(intnhead,intntail,intnlength)
inta[1000];
inta1[1000];
intsp[1000][2];
boolflag=false;
memset(a,0xff,4000);
memset(a1,0xff,4000);
memset(sp,0xff,8000);
if(STEP[head].length%2==0)
flag=true;
while(STEP[head].length==nlength-1)
{
funTspt(STEP[head].nMer,STEP[head].nSer,flag);
head++;
for(inti=0;
head+1;
a[(STEP[i].length)]=i;
a1[(STEP[i].length)]=i;
sp[0][0]=sp[0][1]=n;
STEP[head].nMer=STEP[head].nSer=0;
inttop=0;
intcount=0;
show(a1,sp,top,count,a);
boolA:
funLeft(Nodend,intb1,intb2,intn)
boolflag=abs(nd.nMer-b1)+abs(nd.nSer-b2)<
nB+1
&
abs(nd.nMer-b1)+abs(nd.nSer-b2)>
0;
if(flag==false)
returnfalse;
if(n%2==0&
b1>
=nd.nMer&
b2>
=nd.nSer)
returntrue;
if(n%2==1&
b1<
b2<
returnfalse;
Tspt()
{
Node*temp=newNode;
temp=NULL;
while(head<
=tail)
if(STEP[head].length%2==0)
flag=true;
else
flag=false;
temp=funTspt(STEP[head].nMer,STEP[head].nSer,flag);
if(NULL!
=temp)
break;
if(head>
tail)
此问题无解!
exit
(1);
doLeft(temp->
nMer,temp->
nSer,temp->
length);
//temp->
nMer表示head
deletetemp;
Node*A:
funTspt(intnm,intns,boolflag)
{//flag==true向对岸运输
Node*nd=NULL;
inttemp=1;
inttM=STEP[head].nMer;
//可供运输的商人数
inttS=STEP[head].nSer;
//可供运输的仆人数
if(flag==false)//向此岸运输
tM=n-STEP[head].nMer;
tS=n-STEP[head].nSer;
temp=-1;
tM+1&
nB+1;
i++)//i表示运输的商人数
for(intj=0;
j<
tS+1&
nB-i+1;
j++)//j表示运输的仆人数
if(i+j==0)
continue;
intp=STEP[head].nMer-temp*i;
intq=STEP[head].nSer-temp*j;
if(islegal(p,q)==true&
noRepeat(p,q)==true)
if(p==0&
q==0)
{
tail++;
STEP[tail].length=STEP[head].length+1;
STEP[tail].nMer=p;
STEP[tail].nSer=q;
nd=(Node*)malloc(sizeof(Node));
nd->
length=STEP[head].length+1;
nMer=head;
nSer=tail;
returnnd;
}
tail++;
STEP[tail].length=STEP[head].length+1;
STEP[tail].nMer=p;
STEP[tail].nSer=q;
}
}
returnnd;
noRepeat(intnm,intns)
intj1=0;
j1=1;
for(inti=j1;
tail+1;
if(STEP[i].length%2==j1&
nm==STEP[i].nMer&
ns==STEP[i].nSer)
returnfalse;
returntrue;
islegal(intnm,intns)
{//商人数少于仆人数或者商人数为0
if((nm==0)||(nm==n)||(nm==ns))
funshow(inta[][2],intntail)
商人数仆人数"
ntail;
第"
i<
次"
a[i][0]<
"
a[i][1]<
if(i!
=ntail-1&
i%2==0)
-->
("
abs(a[i+1][0]-a[i][0])<
"
<
abs(a[i+1][1]-a[i][1])<
)"
elseif(i!
i%2==1)
<
--("
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 过河 问题 数学 建模