算法分析与设计——LCS最长公共子串问题和完全背包问题

---

一、问题

LCS（最长公共子序列问题）：
给定序列X=<x1,x2,…，xm>, Y=<y1,y2,…，yj>,求X和Y的最长公共子序列

基础背包问题：
一个旅行者随身携带一个背包，可以放入n种物品，其中物品j的重量和价值分别是wj，vj，j=1，2，…，n，如果背包最大重量限制是b，怎么选择放入背包的物品使得背包的价值最大？

二、解析

LCS（最长公共子序列问题）：
综上，状态转移方程如下：
如果xi=yi，c[i][j] c[i-1][j-1]+1，因为第i位和第j位相同，所以最长序列长度要加上这个相等

如果xi!=yi，c[i][j]= c[i-1][j], c[i][j-1]，当第i位和第j位不相同时，则这个长度要从y的前一位j-1与x的i位或者x的前一位i-1与y的j位中能选出的最长子序列长度中选

背包问题：
设Fk(y)表示只允许装前k种物品时，背包总重不超过y时背包的最大价值
则在得到Fk(y)之前的状态为Fk-1(y1)，在此状态下想得到Fk(y)的策略是：

1. 继续装k-1号物品
2. 装k号物品
   
   Fk(y)表示只允许装前 k 种物品，背包总重不超过 y（重量）时背包的最大价值

不装第 k 种物品： Fk-1(y)

至少装 1 件 k 中物品：仅仅装一件物品 k，排除一件物品 k 后(减去重量,加上其价值)递归(可能再减一件物品 k，或者 k-1 物品),物品 k 的重量和价值分别为 wk,vk;Fk(y-wk)表示在(y-wk)重量中,(已排除一件 k 物品的前提下）只允许装前 k 种物品时的最大价值。

F0(y)表示背包不装物品的价值
Fk(0)背包重量限制为 0 时的最大价值
F1(y)只装第一种物品，[y/w1]个数
Fk(y)=-∞用于比较

转移方程：
F[i][j]=max(F[i-1][j],F[i][j-w[i]]+v[i])，此方程还可以压缩成一维数组

举例列表：

三、设计

LCS：

void lcs()
{
    // 先定义初始状态，x0与yi以及xi与y0因为有一个序列长度为0，所以最长子序列一定为0
    for (int i = 0; i <= m; ++i) c[i][0] = 0;
    for (int i = 0; i <= n; ++i) c[0][i] = 0;

    //开始用状态转移方程计算
    for (int i = 1; i <= m; ++i)
    {
        for (int j = 1; j <= n; ++j)
        {
            if (x[i] == y[j]){
                res[i][j] = 0;
                c[i][j] = c[i - 1][j - 1] + 1;
            }
            else
            {
                if (c[i - 1][j] > c[i][j - 1])
                {
                    res[i][j] = 1;
                    c[i][j] = c[i - 1][j];
                }
                else
                {
                    res[i][j] = 2;
                    c[i][j] = c[i][j - 1];
                }
            }
        }
    }
}

string getlcs()
{
    int i = m, j = n;
    string ans = "";
    while (i > 0 && j > 0)
    {
        if (res[i][j] == 0)
        {
            ans = x[i] + ans;
            i--, j--;
        }
        else if (res[i][j] == 1)
            i--;
        else if (res[i][j] == 2)
            j--;
    }
    return ans;
}

完全背包：

    for (int i = 1; i <= n; i++) 
    {
		for (int j = 1; j <= b; j++) 
        {
			if (w[i] > j) f[i][j] = f[i - 1][j];
			else f[i][j] = max(f[i - 1][j], f[i][j - w[i]] + v[i]);
		}
	}

四、分析

LCS需要两重循环DP，第一重为串A的长度M，第二重为串B的长度N，所以时间复杂度为O(MN)

完全背包需要两重循环DP，第一重为物品数量N，第二重为背包容量W，所以时间复杂度为O(N*W)

五、代码

LCS:

#include<iostream>
#include<vector>
#include<cmath>
#include<time.h>
#include<cstring>
#include<algorithm>
 
using namespace std;

int m, n;
int c[1000][1000];//记录当取到xi位和yj位时最长子序列的长度
string x, y;//x,y两个子序列
int res[1000][1000];//标记数组，记录判断xi？=yj的状态

void lcs()
{
    // 先定义初始状态，x0与yi以及xi与y0因为有一个序列长度为0，所以最长子序列一定为0
    for (int i = 0; i <= m; ++i) c[i][0] = 0;
    for (int i = 0; i <= n; ++i) c[0][i] = 0;

    //开始用状态转移方程计算
    for (int i = 1; i <= m; ++i)
    {
        for (int j = 1; j <= n; ++j)
        {
            if (x[i] == y[j]){
                res[i][j] = 0;
                c[i][j] = c[i - 1][j - 1] + 1;
            }
            else
            {
                if (c[i - 1][j] > c[i][j - 1])
                {
                    res[i][j] = 1;
                    c[i][j] = c[i - 1][j];
                }
                else
                {
                    res[i][j] = 2;
                    c[i][j] = c[i][j - 1];
                }
            }
        }
    }
}

string getlcs()
{
    int i = m, j = n;
    string ans = "";
    while (i > 0 && j > 0)
    {
        if (res[i][j] == 0)
        {
            ans = x[i] + ans;
            i--, j--;
        }
        else if (res[i][j] == 1)
            i--;
        else if (res[i][j] == 2)
            j--;
    }
    return ans;
}

int main()
{
    cin >> x >> y;
    m = x.size(), n = y.size();
    //cout << m << " " << n << endl;
    x = " " + x;
    y = " " + y;
    lcs();
    int ans = c[m][n];
    cout << "最长公共子串长度： " << ans << endl;
    cout << "最长公共子串： " << getlcs() << endl;
    system("pause");
    return 0;
}

/*
测试案例
abcdfge
abcdefg
*/

完全背包:

#include<iostream>
#include<vector>
#include<cmath>
#include<time.h>
#include<cstring>
#include<algorithm>
 
using namespace std;

#define N 20
int f[N+2][N+2];//全局变量，自动初始化为0  
int w[N];  
int v[N];  
int main(){  
    int n,b;
    cout << "请输入物品种类个数" << endl;
    cin>>n;//物品种类
    cout << "请输入背包最大容量" << endl;
    cin>>b;//背包容量
    cout << "请输入每个种类物品的重量和价值" << endl;
    for (int i=1;i<=n; i++)  
    {  
        cin >> w[i] >> v[i];  
    }  
    for (int i = 1; i <= n; i++) 
    {
		for (int j = 1; j <= b; j++) 
        {
			if (w[i] > j) f[i][j] = f[i - 1][j];
			else f[i][j] = max(f[i - 1][j], f[i][j - w[i]] + v[i]);
		}
	}
    cout << "当装"<<n<< "种物品时的最大价值是" <<f[n][b] <<endl;//输出最优解
    system("pause");
    return 0;
}

/*
测试案例
4
10
2 1
3 3
4 5
7 9
*/