项湫

导航

装载问题

1. 问题

 

 

 

整数规划问题,0-1 背包问题

 

2. 解析

思路:轻者先装,直到再装任何集装箱将使轮船载重量超过 C 时停止。

定理:对于任何正整数 k,算法(轻者先装)对 k 个集装箱的实例得到

最优解。

证明:(数学归纳法)

1)k=1,只有 1 个集装箱,其重量小于 C。任何装法都只有一种方式,

因此都是最优解,因此轻者先装也是最优解。

2)归纳假设:假设算法对于规模为 k 的输入都能得到最优解。

考虑规模为 k+1 的输入, N={1,2,...,k+1},W={w1,w2,...,wk+1}是集装箱的重量,w1,w2,...,wk+1 。

N 中拿掉最轻的集装箱,得到 k 规模的输入:

N’=N-{1}={2,3,4,...,k+1}

W’=W-{w1}

C’=C-w1

根据归纳假设,对于 k 个输入,N’、W’、C’的最优解为 I’,即

I’为 N’,不含 1 的最优解(归纳假设得),

I=I’U{1}

那么 I 必然是 N 的最优解,这也是算法对于 N,W,C 的解。

证明:I 必然是 N 的最优解,采用反证法,即假设 I 不是 N 的最优解。

1)构建最优解 I*(N,含 1):假设 I 不是 N 的最优解。则必然存在

最优解 I*,如果 I*中没有 1,用 1 替代 I*中的第一个集装箱标号得到的解也是

最优解(个数不变,因此也是最优解),使得 I*为包含 1 的关于 N 的最优解

|I*|>|I|。

2) 构建最优解 I*’( N’,不含 1):因为 I*为包含 1 的关于 N

的最优解,构建的 I*’=I*-{1}是不包含 1 的最优解(待证明),即关于 N’、

W’、C’的最优解N’、W’、C’不包含 1)。

证明:I*’=I*-{1}是不包含 1 的最优解。

证(反证法):如果 I*’=I*-{1}不是不包含 1 的最优解,那么存在最优

I*’’ ( N’,不含 1,使得|I*’|<|I*’’|,

|I*’+{1}|<|I*’’ +{1}|,与 I*’+{1}= I*为最优矛盾,因此

I*’=I*-{1}是不包含 1 的最优解。

 

3)构建的最优解 I*’与归纳假设的最优解 I’比较:由(1)

|I*|>|I|得,| I*’|= |I*-{1}|>|I-{1}|=|I’|,与 I’的最优性矛盾(最优解 I*-

{1}大于最优解 I’)。

 

3. 设计

4. 分析

O(nlogn)

 

5. 源码

#include<stdio.h>

#define M 100

int n=3;

int x[]={0,0,0};

int c1=50;

int w[]={10,40,40};

int x1[]={0,0,0};

 

int Constrain(int t){

 int i;

 int sum=0;

 static int sum1=0;

 

 

 for(i=0;i<=t;i++){

 if(x[i]==1)

  sum+=w[i];

 }

 printf("--------------->%d\n",sum);

 for(i=0;i<n;i++){

   printf("%d",x1[i]);

   }

   printf("\n");

 

 if(sum>c1){

 x[t]=0;

 return 0;

 }

 

  else {

   if(sum1<sum){

    for(i=0;i<=t;i++){

    if(x[i]==1)  x1[i]=1;

    else x1[i]=0;

   }

   }

  return 1;

  }

}

 

int  Bound(int t){

 if(t<n)

  return 1;

 else return 0;

}

void Backtrack(int t){

 

 int i;

 

 if(t<n){

  for(i=0;i<=1;i++){

  x[t]=i;

  if(Constrain(t)&&Bound(t))

   Backtrack(t+1);

  }

 }

}

int main(){

 int i;

 int sum=0;

Backtrack(0);

 for(i=0;i<n;i++){

 if(x1[i]==0)

sum+=w[i];

 }  

if(sum>c1){

printf("无解");

return 0;

}

  for(i=0;i<n;i++){

 if(x1[i]==1)

  printf("%d  ",w[i]);  

 }  

 

}

posted on 2021-05-17 15:40  项湫  阅读(170)  评论(0编辑  收藏  举报