动态规划

动态规划算法适应于解的最优化问题，通常可按以下步骤求解：

     （1）找出最优解的性质，并刻画其结构特征。--即其最优解含有子问题的最优解。

     （2）递归的定义最优质。--即找出动态规划方程。

     （3）自底向上的方式求解最优值。

     （4）根据计算最优值时得到的信息，构造最优解。

一、LCS（输出LCS的长度和其中一个LCS）

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
using namespace std;
int main(){
    char a[100],b[100];
    int dp[100][100];
    while(cin>>(a+1)>>(b+1)){
        memset(dp,0,sizeof(dp));
        int len1=strlen(a+1),len2=strlen(b+1);
        for(int i=1;i<=len1;i++){
            for(int j=1;j<=len2;j++){
                if(a[i]==b[j]) dp[i][j]=dp[i-1][j-1]+1;
                else dp[i][j]=max(dp[i][j-1],dp[i-1][j]);
            }
        }
        cout<<"the LCS's length is:"<<dp[len1][len2]<<endl;
        string ans;
        int i=len1,j=len2;
        while(i>=1&&j>=1){
            if(dp[i][j]==dp[i-1][j-1]+1&&dp[i][j]>dp[i-1][j]&&dp[i][j]>dp[i][j-1]){
                ans+=a[i];
                i--;
                j--;
            }
            else if(dp[i][j]==dp[i][j-1]) j--;
            else if(dp[i][j]==dp[i-1][j]) i--;
        }
        cout<<"one of the LCSs is: "<<endl;
        for(int i=ans.size()-1;i>=0;i--)
            cout<<ans[i];
        cout<<endl;
    }
    return 0;
}

二、0-1背包问题（输出可以拿到的最大重量和所选物品）

#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
int Knapsack(int n,int m,int *w,int *v){
    int dp[100][100];
    memset(dp,0,sizeof(dp));
    for(int i=1;i<=n;i++){
        for(int j=w[i];j<=m;j++)
            dp[i][j]=max(dp[i-1][j],dp[i-1][j-w[i]]+v[i]);
    }
//    for(int i=1;i<=n;i++){
//        for(int j=0;j<=m;j++)
//            cout<<dp[i][j]<<' ';
//        cout<<endl;
//    }
    cout<<"被选中的物品有:"<<endl;
    int i=n,j=m;
    while(i>0&&j>0){
        if(dp[i][j]!=dp[i-1][j]){
            cout<<"第"<<i<<"件物品"<<endl;
            j-=w[i];
        }
        i--;
    }
    return dp[n][m];
}
int main(){
    int n,w[100],v[100],m;
    cout<<"请输入物品总数和背包容量:"<<endl;
    while(cin>>n>>m){
        cout<<"请输入物品的重量和价值:"<<endl;
        for(int i=1;i<=n;i++)
            cin>>w[i]>>v[i];
        cout<<"可以得到的最大价值为:"<<Knapsack(n,m,w,v)<<endl;
    }
    return 0;
}

三、矩阵连乘问题

#include <iostream>
using namespace std;
class P_MatrixChain{
private:
    int n;      //n个矩阵
    int *p;     //矩阵维数
    int **m;    //最少乘次数
    int **s;    //间断点
public:
    P_MatrixChain(int temp_n,int *temp_p){
        n=temp_n;
        p=temp_p;
        m=new int*[n+1];
        for(int i=0;i<=n;i++) m[i]=new int[n+1];
        s=new int*[n+1];
        for(int i=0;i<=n;i++) s[i]=new int[n+1];
    }
    void MatrixChain();
    void TraceBack(int,int);
};
void P_MatrixChain::MatrixChain(){
    for(int i=1;i<=n;i++) m[i][i]=0;
    for(int i=2;i<=n;i++){
        for(int l=1;l<=n-i+1;l++){
            int r=l+i-1;
            m[l][r]=p[l-1]*p[l]*p[r]+m[l+1][r];
            s[l][r]=l;
            for(int k=l+1;k<r;k++){
                int t=m[l][k]+m[k+1][r]+p[l-1]*p[k]*p[r];
                if(t<m[l][r]) {
                    m[l][r]=t;
                    s[l][r]=k;
                }
            }
        }
    }
}
void P_MatrixChain::TraceBack(int l,int r){
    if(l==r) return ;
    TraceBack(l,s[l][r]);
    TraceBack(s[l][r]+1,r);
    cout<<"Multiply A"<<l<<','<<s[l][r];
    cout<<" and A"<<s[l][r]+1<<','<<r<<endl;
}
int main(){
    int *p,n;
    cin>>n;
    p=new int[n+1];
    for(int i=0;i<=n;i++) cin>>p[i];
    P_MatrixChain mx(n,p);
    mx.MatrixChain();
    mx.TraceBack(1,n);
    return 0;
}