贪心算法的设计技术
用贪心法 求问题的解
【例7-1】
问题分析
命题：
那么B是S’的一个最优解。若不然，假如S’有解B’，B’>B，那么用B’替换B以后得到解{i1 ,i2 , …,ik}∪B’，将比S的活动更多，这与S是最优解矛题得证。
计算模型
struct Active{
startTime s;//开始时间
endTime e;//结束时间
selectflag f;//选标识
}A[n];
（2）计算：
活动i 与 活动j 相容 -> A[ i ].s >= A[ j ].e 或 A[ j ].s >= A[ i ].e
用归并排序 或 其他任何高效的排序算法完成 以a[i].e的从小到大 的排序，形成排序A[1].e≤A[2].e ≤…≤A[n].e
【例7-2】数列极差。
问题分析
命题：
证明：
(1) 当k=3时，
数列a={a1 , a2 , a3}，不妨令a1< a2< a3 ，这样
——思考：
A={a1,a2},B={b1,b2,b3}, 设[A]=[a1 a2]=a1*a2+1 一种值在运算中{B}={b1,b2,b3}表示三个 顺序 ，即[b1 b2 b3], [b1 b3 b2], [b2 b3 b1]，若有A<a3<B<a4, 试比较 [[A],a3,{B},a4]与 [[A],a4,{B},a3]猜测：[[A],a3,{B},a4] > [[A],a4,{B},a3]
因为a3<a4 ，设a4=a3+k , a3=a4-k
[[A],a3,{B},a4]
=[ [A] ,a4-k ,{B} ,a4 ]
= [ [A]*(a4-k)+1 , {B},a4 ]
=[ [ A,a4]-[A]*k ,{B},a4 ] ——（1）B={b1,b2,b3} 取任意一个次序
[[A],a3,{B},a4] ——代入（1）式
=[ [A,a4]-[A]*k ,{ b1, b2,b3,b4},a4 ]
=[ [A,a4]-[A]*k ,b1, {b2,b3,b4},a4 ]
= [ ( [A,a4]-[A]*k )*b1+1 , {b2,b3,b4},a4 ]
=[ [A,a4,b1] -[A]*k*b1 , {b2,b3,b4},a4 ]
=.....
= [ [A,a4,b1,b2,b3]-[A]*k*b1*b2*b3 ,a4 ]
=[ [A,a4,b1,b2,b3] *a4 -[A]*k*b1*b2*b3*a4 +1]
=[ [A,a4,b1,b2,b3] *a4 +1 -[A]*k*b1*b2*b3*a4 ]从结果倒推，找[[A],a4,{B},a3] 形式
aj=ai+k ,a4=a3+k ,
[[A],a3,{B},a4]
=[ [A,a4,b1,b2,b3] *a4 +1 -[A]*k*b1*b2*b3*a4 ] (上文，代入a4=a3+k)
=[ [A,a4,b1,b2,b3] *(a3+k) +1 -[A]*k*b1*b2*b3*a4 ]
=[ [A,a4,b1,b2,b3] *a3+[A,a4,b1,b2,b3]*k +1 -[A]*k*b1*b2*b3*a4 ]
=[ [A,a4,b1,b2,b3,a3]+[A,a4,b1,b2,b3]*k -[A]*k*b1*b2*b3*a4 ]
=[ [A,a4,b1,b2,b3,a3]+k*( [A,a4,b1,b2,b3] -[A]*b1*b2*b3*a4 ) ]
= [ [A,a4,B,a3+k] +1 - [A]*k*b1*b2*b3*a4 ]
=[ [A,a4,B] *(a3+k)+1 - [A]*k*b1*b2*b3*a4 ]
=[ [A,a4,B] *a3+1+[A,a4,B] *k - [A]*k*b1*b2*b3*a4 ]
=[ [A,a4,B,a3] + [A,a4,B] *k - [A]*k*b1*b2*b3*a4 ]
= [ [A,a4,B,a3] + k* ( [A,a4,B] - [A]*b1*b2*b3*a4 ) ]
只需比较 [A,a4,B] - [A]*b1*b2*b3*a4 与0
( [A]*a4 +1 )*B+1 -[A]*b1*b2*b3*a4
=展开或者非常显然 >0因此[[A],a3,{B},a4] > [[A],a4,{B},a3]
设A={a1,a2...an}
[A,ai]
= [{a1,a2,...an} ,ai ]
=[a1,a2,a3...an,ai]
=[ a1*a2+1,a3...an,ai]
=[ ((a1*a2)+1 )*a3+1,....an,ai]
=[ [A],ai ]
(2)假设k=n时命题成立。
——引理
设有数列集合A和B，正整数a i , a j ，且有A<a i < a j <B，其中，B={b 1 , b 2 , …b m }，{B}表示B中元素的任意组合序列，则有[ [A] , a i , {B}, a j ]> [ [A], a j , {B}, a i ]。引理证明：第一个表达式 用第二个表达的出来，比较差别
∵ai< aj，∵不妨设aj = ai +k, ai = aj – k // （k>0）
ai = aj – k
[[A], ai, {B}, aj]
= [[A], aj – k, {B}, aj]
= [[A]*(aj – k)+1, {B}, aj]
= [[A]*aj –[A]*k+1, {B}, aj] = [ [A]*aj+1 – [A]*k, {B}, aj]
=[ [A, aj] –[A]*k, {B}, aj]将B={b1 , b2 , …bm }代入上式，并取B的[任意]一个次序
[[A], ai, B, aj]
= [[A, aj] –[A]*k, b1 ,{ b2 , …bm}, aj]
= [([A, aj] –[A]*k)*b1 +1,{ b2 , …bm}, aj]
=[ [A, aj] *b1 –[A]*k*b1 +1,{ b2 , …bm}, aj]
= [ [A, aj,b1 ] –[A]*k*b1 ,{ b2 , …bm}, aj]
= [ [A, aj,b1 ,b2 ,…bm] –[A]*k*b1 * b2 *…*bm, aj]
= [ ([A, aj,b1 ,b2 ,…bm]–[A]*k*b1 * b2 *…*bm) * aj +1]
= [ [A, aj,b1 ,b2 ,…bm] * aj–[A]*k*b1 * b2 *…*bm* aj +1]代入aj = ai +k，B= b1,b2,…bm
[[A], ai, B, aj]
= [ [A, aj, B] * (ai+k)+1 –[A]*k*b1* b2*…*bm* aj]
= [ [A, aj, B] * ai+1+[A, aj, B] *k –[A]*k*b1* b2*…*bm* aj]
= [ [A, aj, B, ai] +[A, aj, B] *k –[A]*k*b1* b2*…*bm* aj]
= [ [A, aj, B, ai] + k * ([A, aj, B] –[A]*b1* b2*…*bm* aj)]
//证明+的这个 k * ([A, aj, B] –[A]*b1* b2*…*bm* aj)]>0依据B次序的任意性，可得[[A], a i ,{B}, a j ] = [ [A, a j , {B}, a i ] + k * ([A, a j , B] –[A]*b 1 * b 2 *…*b m * a j )]因为[A, a j]=[A]* aj +1，所以[A, a j, B] >[A]*b1 * b2 *…*bm* a j ，//[ A aj B ]= ( [ A ]*aj +1 )*B+1 =[ A ] * aj*B +B +1则必定有 [[A], a i ,{B}, a j ] > [A, a j , {B}, a i ]∴引理成立。引理推广[{A}, a i ,{B}, a j ] > [{A}, a j , {B}, a i ] 当 [{A}]<a i < a j <[{B}]A的任意顺序引理：[ [A] , ai, {B}, aj ]> [ [A], aj, {B}, ai ]。
根据计算可以得出 [A , ai ] = [ [A], ai ]
但是[ai,A] ！= [ai,[A] ]
比如[ai ,a1,a2]= (ai*a1+1)*a2+1 =ai*a1*a2+ai*a1+a2+1
[ai,[A]]=ai*(a1*a2+1)=ai*a1*a2+ai
按照数列里最小值优先组合的规则，仍可以得到最大值。
任意性可知，加入一个新元素
按数列元素最小值优先组合的规则，仍可以得到最大值。则k=n+1时命题成立
计算模型
代码：
#include<bits/stdc++.h>
using namespace std;
//堆排序 -完全二叉树
/*
0
1 2
3 4 5 6
*/
void show(int a[],int n)
{
for(int i=0;i<n;i++)
{
cout<<a[i]<<"\t";
}
cout<<endl<<endl;
}
void min_h(int a[],int start,int end)
{
//建立父节点指针和子节点指针
int f=start;///父节点
int c=f*2+1;//左孩子
int t;//用于交换的临时变量
while(c<=end)//若子节点指针在范围内，进行比较
{
//找到最小值的下标
if(c+1<=end && a[c]>a[c+1])
//右孩子在范围中，左孩子>右孩子
{
c++;
}//选中左右孩纸中较小值
//如果父节点 < 子节点
//此时 父<右<左 调整完成
if(a[f]<a[c])return;
else//交换 父>右 左>右 ，右是最小的 换到父节点上去
{
t=a[f];a[f]=a[c];a[c]=t;
f=c;
c=f*2+1;
}
}
}
void min(int a[],int n)
{
//初始化 i从最后一个父节点开始调整
for(int i=n/2-1;i>=0;i--)
{
min_h(a,i,n-1);
}
for(int i=n-1;i>0;i--)
{
swap(a[0],a[i]);
min_h(a,0,i-1);
}
}
int cal(int a[], int n)
{
int v1=0,v2=1;
while (n > 2)
{
show(a,n);
min(a, n);//max
a[v1] = a[v1] * a[v2] + 1;
a[v2] = a[n-1];//下次还需排序，因此不用在意顺序，直接将最后一个向前移动
n = n - 1;//数目--
}
//此时n==2
return a[0] * a[1] + 1;
}
int main()
{
int n = 4;
int a[4] = { 8,2,1,3 };
cout<<"最大极差：" << cal(a, n);
return 0;
}
【例7-3】最优装载。
问题分析
贪心策略：轻者优先
| 算法设计与描述 | 算法分析 |
| 输入： | （1）问题输入规模为n |
| 输出： | |
| （2）选择集装箱是主要工作，时间渐进复杂度T(n)=n （3）排序最快的时间复杂度T(n)=O(nlogn) （4）上述（2）（3）是并列执行的，按照第2章的定理2.2可知时间渐进复杂度T(n)=O(nlogn) |
t[]数组存放集装箱的原始编号。通过t 找到x 下标变量
【例7-4】键盘输入
n1=“1 2 4 3 5 8 6 3”n2=“2 3 1 1 8 3 1”n3=“1 2 3 4 5 6 7”n4=“1 2 0 0 8 3”=》n1=“1 2 3 5 3”n2=“2 1 1 1” 应该是1131 需要回溯，只需向上回溯一位n3=“1 2 3 4 5 6 7” 应该直接划掉后三个 没删除， 有木有n4=“1 0 0 3” 没删够， 有木有
计算模型
代码：
#include<bits/stdc++.h>
using namespace std;
int n = 6;//长度
int s = 3;//要删除的数字个数
char a[] = { '1','2','0','0','8','3' };//数字
int data[3] = { 0 };//记录删除的元素在原数字中的位置
void show(char a[])
{
for(int i=0;i<strlen(a);i++)
cout<<a[i]<<'\t';
cout<<endl;
}
void del(char a[], int b, int k)//删除,物理覆盖
{
for (int i = b; i < strlen(a); i++)
{
a[i] = a[i + k];
}
}
void delD(char a[], int s)
{
int i = 0, j = 0, j1 = -1;//存在回溯
int len = strlen(a);//截止到\o
while (i < s && j < strlen(a))
{
show(a);
while (a[j] <= a[j + 1])//一直是前面小
{
j++;
}
//前面大于后面
if (j < strlen(a))
{
//删除a[j]
del(a, j, 1);
//删除次数 i
if (j >= j1)data[i] = j + i;//删除的第i个元素，在原数组中的位置是j+i
else data[i] = data[i - 1] - 1;//回溯，回溯到上一个删除位置 的前一个位置
//
j1 = j;
i++;
j--;
if (j < 0)j = 0;
}
}
//显示的时候高位的0不显示
while (a[0] == '0' && strlen(a))
del(a, 0, 1);
}
int main()
{
delD(a, s);
cout<<"位置"<<endl;
for(int i=0;i<s;i++)
{
cout<<data[i]<<"\t";
}
cout<<endl<<"新整数";
for (int i = 0; i < n-s; i++)
cout << a[i];
return 0;
}
1268
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
