1.猴子吃桃问题(比较简单的一个例子)
(1)案例提出
猴子第一天摘下若干个桃子,当即吃了一半,还不过瘾,又多吃了一个。第二天又将剩下的桃子吃掉了一半,又多吃了一个。以后每天将前一天剩下的桃子吃掉一半,再多吃一个。直到第十天只剩下一个桃子了,求第一天共摘了多少个桃子?
(2)递推设计
设置递推变量:f(k),
递推关系式:f(k+1)=f(k)/2-m,f(k)=2(f(k+1)+m)
初始条件:f(10)=1
约束条件:9–1
#include<stdio.h>
int main()
{
int k,f(10)=1;
for(k=9;k>=1;k--)
{
f(k)=2*(f(k+1)+m);
}
printf("%d\n", f(1));
return 0;
}
(3)递归设计
#include<stdio.h>
int fun(int d)
{
if(d>=10)
{
return 1;
}
else
{
return (fun(d+1)+1)*2;
}
}
int main()
{
int m;
m=fun(1);
printf("%d\n", m);
return 0;
}
运行结果:
2.插入乘号问题(相对复杂)
(1)案例提出
在一个由n个数字组成的数字串中插入r个乘号(1<=r<n<16),将它分成r+1个整数,找出一种乘号的插入方法,使得r+1个整数的乘积最大。(这是书上的一道典型例子)
(2)穷举设计
#include<iostream>
#include<string>
using namespace std;
string str;
int s[10];
int t[10],f[10];
int main()
{
int i,j,count=1;
_int64 d,y,ans=0;
cin>>str;
int k=str.length();
for(i=0;i<k;i++)
s[i+1]=str[i]-'0';
t[0]=0;t[4]=7;
for(t[1]=1;t[1]<=4;t[1]++)
for(t[2]=t[1]+1;t[2]<=5;t[2]++)
for(t[3]=t[2]+1;t[3]<=6;t[3]++)
{
y=1;
for(i=0;i<=3;i++)
{
d=0;
for(j=t[i]+1;j<=t[i+1];j++)
d=d*10+s[j];
y*=d;
}
if(ans<y)
{
ans=y;
for(i=1;i<4;i++)
f[i]=t[i];
}
}
printf("分别插入到%d %d %d之后!!!最大值为%lld\n",f[1],f[2],f[3],ans);
return 0;
}
程序运行:
(3)动态规划求解
1.建立递推关系
设f(i,k)表示在前i位数中插入k个乘号所得乘积的最大值,a(i,j)表示从第i个数字到第j个数字所组成的j-i+1(i<=j)位整数值。
为了寻求递推关系,先看一个实例:对给定的847313926,如何插入r=5个乘号,使其乘积最大?我们的目标是为了求取最优值f(9,5)。
①设前8个数字中已插入4个乘号,则最大乘积为f(8,4)*6;
②设前7个数字中已插入4个乘号,则最大乘积为f(7,4)*26;
③设前6个数字中已插入4个乘号,则最大乘积为f(6,4)*926;
④设前5个数字中已插入4个乘号,则最大乘积为f(5,4)*3926;
比较以上四个数值的最大值即为f(9,5)。
依此类推,为了求f(8,4):
①设前7个数字中已插入3个乘号,则最大乘积为f(7,3)*2;
②设前6个数字中已插入3个乘号,则最大乘积为f(6,3)*92;
③设前5个数字中已插入3个乘号,则最大乘积为f(5,3)*392;
④设前4个数字中已插入3个乘号,则最大乘积为f(4,3)*1392;
比较以上4个数值的最大值即为f(8,4)。
一般地,为了求取f(i,k),考察数字串的前i个数字,设前j(k<=j<i)个数字中已插入k-1个乘号的基础上,在第j个数字后插入第t个乘号,显然此时的最大乘积为f(j,k-1)*a(j+1,i)。于是可以得递推关系式:
f(i,k)=max(f(j,k-1)*a(j+1, i)) (k<=j<i)
前j个数字没有插入乘号时的值显然为前j个数字组成的整数,因而得边界值为:
f(j,0)=a(1,j) (1<=j<=i)
2.递推计算最优值
为简单计,在程序设计中省略a数组,用变量d替代。
for (d = 0, j = 1; j <= n; j++) {
d = d * 10 + b[j - 1]; // 输入数字串每一位赋值给b数组
f[j][0] = d; // 计算初始值f[j][0]
}
for (k = 1; k <= r; k++)
{
for (i = k + 1; i <= n; i++)
for (j = k; j < i; j++) {
for (d = 0, u = j + 1; u <= i; u++)
d = d * 10 + b[u - 1]; // 计算d即为a(j+1,i)
if (f[i][k] < f[j][k - 1] * d) // 递推求取f[i][k]
f[i][k] = f[j][k - 1] * d;
}
printf("最优值为%.0f",f[n][r]);
}
3.构造最优解
为了能打印相应的插入乘号的乘积式,设置标注位置的数组t(k)与c(i,k),其中c(i,k)为相应的f(i,k)的第k个乘号的位置,而t(k)标明第k个乘号“”的位置,例如,t(2)=3,表示第2个“”号在第3个数字后面。
当给数组元素赋值f(i,k)=f(j,k-1)*d时,作相应赋值c(i,k)=j,表示f(i,k)的第k个乘号的位置是j。在求得f(n,r)的第r个乘号位置t(r)=c(n,r)=j的基础上,其他t(k) (1<=k<=r-1)可应用下式逆推产生:
t(k)=c(t(k+1),k)
根据t数组的值,可直接按字符形式打印表面积所求得的插入乘号的乘积式。
4.程序实现代码:
#include <stdio.h>
#include <string.h>
int main()
{
char numStr[16];
int i, j, k, len, u, r, b[16], t[16], c[16][16];
double f[17][17], d;
printf("请输入整数:");
scanf("%s", numStr);
printf("请输入插入的乘号个数r: ");
scanf("%d", &r);
len = strlen(numStr);
if (len <= r) {
printf("输入的整数位够数不够或r太大!/n");
return 0;
}
printf("在整数%s中插入%d个乘号,使乘积最大:\n", numStr, r);
for (d = 0, j = 0; j <= len - 1; j++)
b[j] = numStr[j] - '0';
for (d = 0, j = 1; j <= len; j++) {
d = d * 10 + b[j - 1]; // 把b数组的一个字符转化为数值
f[j][0] = d; // f[j][0]赋初始值
}
for (k = 1; k <= r; k++)
for (i = k + 1; i <= len; i++)
for (j = k; j < i; j++) {
for (d = 0, u = j + 1; u <= i; u++)
d = d * 10 + b[u - 1];
if (f[i][k] < f[j][k - 1] * d) {
f[i][k] = f[j][k - 1] * d;
c[i][k] = j;
}
}
t[r] = c[len][r];
for (k = r - 1; k >= 1; k--)
t[k] = c[t[k + 1]][k]; // 逆推出第k个*号的位置
t[0] = 0; t[r + 1] = len;
for (k = 1; k <= r + 1; k++) {
for (u = t[k - 1] + 1; u <= t[k]; u++)
putchar(numStr[u - 1]); // 输出最优解
if (k < r + 1)
putchar('*');
}
printf("=%.0f\n", f[len][r]); // 输出最优值
return 0;
}
程序运行示例:
这也是我写的算法报告作业,个人感觉算法还是挺有趣的,但是要深入研究算法可能就很难了,在学算法时基本都是用c语言编写程序,推荐一个练习编程题目的网站:南阳理工ACM