实验6 动态规划实验二

OJ练习
1*. Max Sum: http://acm.hdu.edu.cn/showproblem.php?pid=1003

2*. Super Jumping! Jumping! Jumping!: 【最长上升子序列求和】

http://acm.hdu.edu.cn/showproblem.php?pid=1087

3*. Constructing Roads In JGShining’s Kingdom:
http://acm.hdu.edu.cn/showproblem.php? pid=1025
4*. 最少拦截系统: http://acm.hdu.edu.cn/showproblem.php?pid=1257
5*.小明系列问题——小明序列:
http://acm.hdu.edu.cn/showproblem.php?pid=4521
6*.Removed Interval: http://acm.hdu.edu.cn/showproblem.php?pid=5489
7*.LIS again: http://acm.hdu.edu.cn/showproblem.php?pid=5141
8*. Bridging signals: http://acm.hdu.edu.cn/showproblem.php?pid=1950

实验内容
1. 最长递增子序列(Longest Increasing Subsequence, LIS)。给出一个序列a1,a2,a3,a4,a5,a6,a7…an，求它的一个子序列（设为s1,s2,…sn），使得这个子序列满足这样的性质：s1

#include<cstdio>
#include<cstring>
#include<queue>
#include<iostream>
#include<stack>
#include<algorithm>
#define maxn 100005
#define maxm 50005
#define INF 0x3f3f3f3f
#define ll long long
#define MOD 1000000007
using namespace std;

int dp[maxn];
int a[maxn];

int main()
{
    int n;
    while(cin>>n){
    memset(dp,0,sizeof(dp));
    for(int i=1;i<=n;i++){
        cin>>a[i];
    }
    int Max;
    dp[1]=1;
    for(int i=2;i<=n;i++){
        Max=0;
        for(int j=1;j<=i;j++){
            if(a[j]<a[i] && Max<dp[j]){
                Max=dp[j];
            }
        }
        dp[i]=Max+1;
    }
        int ans=0;
        for(int i=1;i<=n;i++)
            ans=max(ans,dp[i]);
        cout<<ans<<endl;
    }
}
/*
7
1 7 3 5 9 4 8
*/

2. 最长递增子序列(Longest Increasing Subsequence, LIS)。在练习1的基础上对代码进行改进，输出一条最长递增子序列，例如：输入{1, 7, 3, 5, 9, 4, 8}，输出{1, 3, 5, 9}或{1, 3, 4, 8}或{1, 3, 5, 8}。

源代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<stack>
#include<algorithm>
#define maxn 100005
#define maxm 50005
#define INF 0x3f3f3f3f
#define ll long long
#define MOD 1000000007
using namespace std;
int a[300],dp[500],pre[500];
using namespace std;
int id;
void printf()
{
    stack<int> s;
    int k=id;
    while(k>=0)
    {
        s.push(a[k]);
        k=pre[k];
    }
    while(!s.empty())
    {
        cout<<s.top()<<" ";
        s.pop();
    }
}
int main()
{
    int n;scanf("%d",&n);
    int ans=0;
    for(int i=0;i<n;i++){
            scanf("%d",&a[i]);
    }
    memset(dp,0,sizeof dp);
    id=1;
    pre[0]=-1;
    for(int i=1;i<=n;i++)
    {
        dp[i]=1;
        for(int j=i-1;j>=1;j--)
        if(a[i]>a[j] && dp[j]+1>dp[i]) {
                dp[i]=dp[j]+1;
                pre[i]=j;
         }
        if(dp[i]>ans)
        {
             ans=dp[i];
             id=i;
        }
    }
    printf();
}

1. 最大子段和。给定n个整数（可能是负数）组成的序列a[1], a[2], a[3], …, a[n]，求该序列的子段和如a[i]+a[i+1]+…+a[j]的最大值。
   源代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<stack>
#include<algorithm>
#define maxn 505
#define maxm 50005
#define INF 0x3f3f3f3f
#define ll long long
#define MOD 1000000007
using namespace std;
int a[maxn],dp[maxn],sum[maxn];
using namespace std;
int id;

int Max(int n)
{
    dp[0]=0;
    int ans=0;
    for(int i=1;i<=n;i++){
        if(dp[i-1]>0) dp[i]=dp[i-1]+a[i];
        else dp[i]=a[i];
        ans=max(ans,dp[i]);
    }
    return ans;
}

int main()
{
    int n;scanf("%d",&n);
    int ans=0;
    for(int i=0;i<n;i++){
        scanf("%d",&a[i]);
    }
    memset(dp,0,sizeof dp);
    printf("%d\n",Max(n));
}

1. 石子合并问题。在一条直线上有n堆石子，每堆有一定的数量，每次可以将两堆相邻的石子合并，合并后放在两堆的中间位置，合并的费用为两堆石子的总数。求把所有石子合并成一堆的最小花费。例如：输入{1，2，3，4，5}，输出33。【3+6+9+15=33】
   源代码：

#include<cstdio>
#include<cstring>
#include<queue>
#include<iostream>
#include<stack>
#define maxn 105
#define maxm 50005
#define INF 0x3f3f3f3f
#define ll long long
using namespace std;

int n;
int a[maxn];
int dp[maxn][maxn];
int sum[maxn];


int main()
{
    while(cin>>n)
    {
        memset(sum,0,sizeof(sum));
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=n;i++) {
            cin>>a[i];
            sum[i]=sum[i-1]+a[i];
        }
        for(int i=1;i<=n;i++)
            dp[i][i]=0;
        for(int r=2;r<=n;r++){
            for(int i=1;i<=n-r+1;i++){
                int j=i+r-1;
                dp[i][j]=INF;
                for(int k=i;k<j;k++){
                    dp[i][j]=min(dp[i][j],dp[i][k]+dp[k+1][j]+sum[j]-sum[i-1]);
                }
            }
        }
        printf("%d\n",dp[1][n]);
    }
    return 0;
}

5*. 设计一个算法，在O(nlogn)的时间内计算一个序列的最长递增子序列的长度。
源代码：

#include<cstdio>
#include<cstring>
#include<queue>
#include<iostream>
#include<stack>
#define maxn 50005
#define maxm 505
#define INF 0x3f3f3f3f
#define ll long long
using namespace std;


int a[maxn];
int dp[maxn];

int main()
{
    int n,tp;
    cin>>n;
    for(int i=0;i<n;i++){
        cin>>a[i];
    }
    int len=0;
    dp[0]=a[0];
    for(int i=1;i<n;i++)
    {
        if(dp[len]<a[i])
            dp[++len]=a[i];
        else{
            int temp = lower_bound(dp, dp+len, a[i]) - dp;
            dp[temp] = a[i];
        }
    }
    printf("%d\n",len+1);
}