田忌赛马问题

郭佳奇：动态规划与数学软件的实现

田忌赛马问题

一 问题描述

田忌与齐王赛马，双方各有n 匹马参赛（n<=100），每场比赛赌注为1两黄金，现已知齐王与田忌的每匹马的速度，并且齐王肯定是按马的速度从快到慢出场，现要你写一个程序帮助田忌计算他最好的结果是赢多少两黄金（输用负数表示）。

算法思想

先排序，齐王的马的速度放在数组a 中，田忌的马的速度放在数组b 中。本问题应用的算法是动态规划和贪心算法相结合解决的。从两人的最弱的马入手：

若田忌的马快，就让这两匹马比赛；

若田忌的马慢，干脆就让他对付齐王最快的马；

若两匹马的速度相等，这时有两种选择方案，或者它俩比赛，或者对付齐王最快的马。 定义子问题：l(i ，j)为齐王的从第i 匹马开始的j 匹马与田忌的最快的j 匹马比赛，田忌所获得的最大收益。

则：⎪⎩

⎪⎨⎧->-+-+--+--+-<-+-=时当时＝当时当]1[]1[)1,1(]1[]1[)}1,(),1,1(max{

]1[]1[)1,(),(j b j i a j i l j b j i a j i l j i l j b j i a j i l j i l 程序具体实现时，为了适合c 数据从0开始，稍加变动，定义子问题：l(i ，j)为齐王的从第i 匹马开始到第i ＋j 匹马共j+1匹马与田忌的最快的j+1匹马比赛，田忌所获得的最大收益。初始化时：l[i][0]表示齐王的第i 匹马与田忌最快的马比赛的结果。

二 程序源代码

#include

void readdata();

void init();

int N,n,a[100],b[100],l[100][100];

void main()

{

int i,j,k;

scanf("%d",&N);//测试例子得个数

for(k=0;k

{

readdata();

init();

for(i=n-2;i>=0;i--)

for(j=1;j

if(a[i+j]

l[i][j]=l[i][j-1]+1;

else if(a[i+j]>b[j])

l[i][j]=l[i+1][j-1]-1;

else if(l[i+1][j-1]-1>l[i][j-1])

l[i][j]=l[i+1][j-1]-1;

运筹学上机报告

else

l[i][j]=l[i][j-1];

printf("%d\n",l[0][n-1]);

}

}

void readdata()

{

int i;

scanf("%d",&n);//马的个数：-

for(i=0;i

scanf("%d",&a[i]);//每只马的速度；

for(i=0;i

scanf("%d",&b[i]);//每只马的速度；

}

int* qsort(int a[100],int n)//对输入的马的速度的无序序列进行排序{

int i,j,t;

for(i=0;i

for(j=i+1;j

if(a[i]

{t=a[i];a[i]=a[j];a[j]=t;}

// for(i=0;i

// printf("%3d",a[i]);

// printf("\n");

return a;

}

void init()

{

int i;

qsort(a,n);

qsort(b,n);

for(i=0;i

{

if(a[i]

l[i][0]=1;

else if(a[i]==b[0])

l[i][0]=0;

else

l[i][0]=-1;

}

}

三动态规划的求解方法做出总结

用动态规划解决多阶段决策问题效率是很高的，而且思路清晰简便，同时易于实现，虽