P1868 饥饿的奶牛

题目描述

有一条奶牛冲出了围栏，来到了一处圣地（对于奶牛来说），上面用牛语写着一段文字。

现用汉语翻译为：

有 NN 个区间，每个区间 x,yx,y 表示提供的 x\sim yx∼y 共 y-x+1y−x+1 堆优质牧草。你可以选择任意区间但不能有重复的部分。

对于奶牛来说，自然是吃的越多越好，然而奶牛智商有限，现在请你帮助他。

输入格式

第一行一个整数 NN。

接下来 NN 行，每行两个数 x,yx,y，描述一个区间。

输出格式

输出最多能吃到的牧草堆数。

输入输出样例

输入 #1复制

3
1 3
7 8
3 4

输出 #1复制

5

说明/提示

1 \leq n \leq 1.5 \times 10^51≤n≤1.5×105，0 \leq x \leq y \leq 3 \times 10^60≤x≤y≤3×106。

线性动态规划。

思路1：91分，超时1个点  时间复杂度O(N*N)

dp[i]:表示以第i个点为结尾获取的最长区间长度。

注意：既然第i个点为结尾，那么该点的y值应该是最大的，所以所有点应该按照y值从小到大排序,防止后效性。

dp[i]:表示以第i个点为结尾获取的最长区间长度。

转移方程：dp[i]=max(dp[j])+a[i].len  其中1<=j<i，a[j].y<a[i].x

#include<bits/stdc++.h>
using namespace std;
struct node{
    int x,y,len;
}a[150005];
int dp[150005],ans,k;

bool cmp(const node&a, const node&b){
    return a.y<b.y;
}
int main()
{
    int n;
    cin>>n;
    for(int i=1; i<=n; i++){
        cin>>a[i].x>>a[i].y;
        a[i].len=a[i].y-a[i].x+1;
    }
    sort(a+1,a+1+n,cmp);//必须先排序
    ans=dp[1]=a[1].len;
    for(int i=2; i<=n; i++){
       k=0;
       for(int j=1; j<i; j++){//每次在i前面找一个最大的
          if(a[i].x>a[j].y){
             k=max(k,dp[j]);
          }
       }
       dp[i]=a[i].len+k;//最大的加上自己的长度
       ans=max(ans,dp[i]);
    }
    cout<<ans<<endl;
    return 0;
}

 

思路二：DP+二分

要用二分优化，重点要理解二分的前提：数列存在单调性。

1.还是要根据y值从小到大排序。

2.状态: dp[i]:前i个区间能得到的最大长度。

  对于第i堆：

    (1)不取  dp[i]=dp[i-1];

    (2)取：MAX(dp[j]+a[i].len)   1<=j<i,且a[j].y<a[i].x

   综上可得 dp[i]=max(dp[j]+a[i].len,dp[i-1]);  1<=j<i,且a[j].y<a[i].x

通过上面的状态转移方程，可以得到dp[i]是单调不降的。

证明：

  (1)若不取  dp[i]=dp[i-1];  不降,成立

  (2)取：MAX(dp[j]+a[i].len)   1<=j<i,且a[j].y<a[i].x

      若存在j,则 dp[j]<dp[i]  j<i;  递增，成立

     若不存在j，则dp[i]=max(dp[i-1],a[i].len);即要么只取第i个（这种情况dp[i-1]<a[i].len=dp[i], dp[i]>dp[i-1].成立），要么不取第i个（dp[i]=dp[i-1]）

故dp[i]单调不降。

 

3.有了2 的结论，那么不需要通过循环来一个个找j了，二分找到离i最近的符合a[j].y<a[i].x的值即可，时间复杂度O(n*logn);

#include<bits/stdc++.h>
using namespace std;
struct node{
    int x,y,len;
}a[150005];
int dp[150005],ans,k;

bool cmp(const node&a, const node&b){
    return a.y<b.y;
}

int solve(int l, int r,int key)//在(1,i-1)找到最后一个a[j].y<a[i]
{//注意这里的l，r是寻找区间（1，i-1）左端点1的左边一个0，右端点的右边一个i
    int m,ans=0;
    while(l+1!=r){
        m=(l+r)>>1;
        if(a[m].y<key)
        {
            l=m;
            ans=l;
        }
        else r=m;
    }
    return ans;
}
int main()
{
    int n;
    cin>>n;
    for(int i=1; i<=n; i++){
        cin>>a[i].x>>a[i].y;
        a[i].len=a[i].y-a[i].x+1;
    }
    sort(a+1,a+1+n,cmp);
    dp[0]=0;
    for(int i=1; i<=n; i++){
       dp[i]=max(dp[i-1],max(a[i].len,a[i].len+dp[solve(0,i,a[i].x)]));//所以这里的搜索区间是（0，i）
    }
    cout<<dp[n];
    return 0;
}