题目描述
有一条奶牛冲出了围栏,来到了一处圣地(对于奶牛来说),上面用牛语写着一段文字。
现用汉语翻译为:
有 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;
}