插花问题的“动态规划法”算法

// Layout.cpp : Defines the entry point for the console application.
/*
作者：成晓旭
时间：2001年10月11日(11:35:38-12:35:00)
内容：完成插花问题的“动态规划法”算法及注解
*/
#include "stdafx.h"
#include "string.h"
#define MAX(A,B) ((A) > (B) ? (A):(B))
//--------------------鲜花问题--------------------
#define F 100
#define V 100
/*
插花问题描述:
将f束鲜花插入v个花瓶中,使达到最徍的视觉效果,
问题相关约定及插花要求:
鲜花被编号为1--f,花瓶被编号为1--v,花瓶按从小到
大顺序排列,一只花瓶只能插一支花,鲜花i插入花瓶j中的
视觉效果效果值已知,编号小的鲜花所放入的花瓶编号也小
问题求解思路:
花瓶j(1<=j<=v)中插入鲜花的可能编号为[1..j](编号
小的鲜花所放入的花瓶编号也小);
设数组p[i][j]表示鲜花i插入花瓶j的好看程度,数组
q[i][j]表示[1..i]束鲜花插入[1..j]个花瓶所能得到的最大
好看程度,初始化q[0][0] = 0;q[0][j]=0(1<=j<=v),则q[f][v]
是问题的解.
特别地,j束鲜花插入到前面的j只花瓶中,所得到的好看
程度是q[j][j] = p[1][1]+p[2][2]+...+[j][j].现将插花过
程按花瓶排列顺序划分成不同阶段,则在第j阶段,第i束鲜花
若放入第j号花瓶,最大好看程度是q[i-1][j-1]+p[i][j];第i束鲜
花若放入前j-1个花瓶中的某一个,所得的好看程度是q[i][j-1],
那么在第j阶段,插入第i束鲜花所能得到的最大好看程度为:
q[i][j] = MAX(q[i-1][j-1]+p[i][j],q[i][j-1]),要使q[i][j]
最大,应使q[i-1][j-1]和q[i][j-1]也最大
*/
//初始化函数
void Initialize(int *f,int *v,int p[][V])
{
int i,j;
printf("输入鲜花数量及花瓶个数:");
scanf("%d%d",f,v);
printf("顺序输入各鲜花插入各花瓶的好看程度:/n");
for(i=1;i<=*f;i++)
for(j=1;j<=*v;j++)
p[i][j] = i+j;
//scanf("%d",&p[i][j]);
}
//鲜花问题处理函数
int Sove(int p[][V],int f,int v,int *way)
{
int i,j,newv,q[F][V];
q[0][0] = 0;
/*设置v个花瓶分别被插入v束鲜花时各号花瓶对应的(初始)最大好看程度*/
for(j=1;j<=v;j++)
{
q[0][j] = 0;
/*设置第j束鲜花放入第j号花瓶中的最大好看程度*/
q[j][j] = q[j-1][j-1]+p[j][j];
}
for(j=1;j<=v;j++)
for(i=1;i<j;i++)
q[i][j] = MAX(q[i-1][j-1]+p[i][j],q[i][j-1]);
newv = v;
for(i=f;i>0;i--)
{
while(q[i-1][newv-1]+p[i][newv] < q[i][newv])
newv--;
//确定鲜花i插在花瓶newv中,并准备考虑前一只花瓶　
way[i] = newv--;
}
return(q[f][v]);
}
//--------------------鲜花问题--------------------
//--------------------最长子串问题--------------------
#define N 100
char a[N],b[N],str[N];
//计算两个序列最长公共子序列的长度
int Get_LongSubStr_Len(char *a,char *b,int c[][N])
{
int m=strlen(a),n=strlen(b),//两个序列的长度
i,j;//循环变量
for(i=0;i<=m;i++)c[i][0] = 0;
for(i=1;i<=n;i++)c[0][i] = 0;
for(i=1;i<=m;i++)
for(j=1;j<=n;j++)
if(a[i-1]==b[j-1])
c[i][j] = c[i-1][j-1]+1;
else
c[i][j] = MAX(c[i-1][j],c[i][j-1]);
/*
if(c[i-1][j]>=c[i][j-1])
c[i][j] = c[i-1][j];
else
c[i][j] = c[i][j-1];
*/
return(c[m][n]);
}
//构造最长公共子序列
char *Build_LongSubStr(char s[],char *a,char *b)
{
int i=strlen(a),j=strlen(b),
k,c[N][N];
k = Get_LongSubStr_Len(a,b,c);
s[k] = '/0';
while(k>0)
{
if(c[i][j]==c[i-1][j])
i--;
else
{
if(c[i][j]==c[i][j-1])
j--;
else
{
s[--k]=a[i-1];
i--;
j--;
}
}
}
return(s);
}
//--------------------最长子串问题--------------------
int main(int argc, char* argv[])
{
int i,f,v,p[F][V],way[F];
//-----------------------------------
/*
Initialize(&f,&v,p);
printf("最大好看程度为%d/n",Sove(p,f,v,way));
printf("插有鲜花的花瓶是:/n");
for(i=1;i<=f;i++)
printf("%4d",way[i]);
*/
//-----------------------------------
printf("输入两个字符串(长度<%d):/n",N);
scanf("%s%s",a,b);
printf("两个串的最长公共子序列是:%s/n",Build_LongSubStr(str,a,b));
//-----------------------------------
printf("/n/n应用程序正在运行....../n");
return 0;
}