bzoj2369/bzoj2687(贪心+单调决策性DP)

显然取比较大的区间比较好,而且从均值不等式的角度能感觉到取的区间越多价值越小。。所以就猜想取2个是最优的。。

设已经取2个较大的区间,再加入一个较小的区间,那么模拟一下可以发现,最理想的情况是交增加的和并减少的是相同的。。

然后还是根据均值不等式可以知道,这个价值会减小,所以去2个区间必然是最优的。。

然后既然要取较大的区间,那么包含的区间也得去了。。而被包含区间所能产生的贡献,只能是被包含的区间和包含区间产生的价值。假设被包含区间和其他区间能产生更大的价值,那么显然包含区间和这个区间产生的价值会更大。。

然后就变成几个交叉区间之间的问题了。。

设d[i]为第i个区间和前i-1个区间产生的最大价值

d[i]=max{(r[i]-l[j])*(r[j]-l[i])}=max{r[i]r[j]+l[i]l[j]+r[j]l[j]}-l[i]r[i]

由于l[i]和r[i]都是递增的

那么对2个决策点j<k<i,如果r[i]r[j]+l[i]l[j]+r[j]l[j]-l[i]r[i]<r[i]r[k]+l[i]l[k]+r[k]l[k]-l[i]r[i]

那么由于l[j]<l[k],r[j]<r[k],此后k的值必定j优,所以这样这个dp具有决策单调性。。

然后用单调队列维护+二分就可以把复杂度降到O(nlogn)

 

 

 

/*
 *        ┏┓    ┏┓
 *        ┏┛┗━━━━━━━┛┗━━━┓
 *        ┃       ┃  
 *        ┃   ━    ┃
 *        ┃ >   < ┃
 *        ┃       ┃
 *        ┃... ⌒ ...  ┃
 *        ┃       ┃
 *        ┗━┓   ┏━┛
 *          ┃   ┃ Code is far away from bug with the animal protecting          
 *          ┃   ┃   神兽保佑,代码无bug
 *          ┃   ┃           
 *          ┃   ┃        
 *          ┃   ┃
 *          ┃   ┃           
 *          ┃   ┗━━━┓
 *          ┃       ┣┓
 *          ┃       ┏┛
 *          ┗┓┓┏━┳┓┏┛
 *           ┃┫┫ ┃┫┫
 *           ┗┻┛ ┗┻┛
 */
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#include<cmath>
#include<map>
#include<stack>
#include<set>
#include<bitset>
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,l,r) for(int i=l;i>=r;i--)
#define link(x) for(edge *j=h[x];j;j=j->next)
#define mem(a) memset(a,0,sizeof(a))
#define ll long long
#define eps 1e-8
#define succ(x) (1LL<<(x))
#define lowbit(x) (x&(-x))
#define sqr(x) ((x)*(x))
#define mid (x+y>>1)
#define NM 1000005
#define nm 200005
#define pi 3.1415926535897931
using namespace std;
const ll inf=1e9+7;
ll read(){
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return f*x;
}
 
 
 
 
struct tmp{
    int l,r;
    bool operator<(const tmp&o)const{return l<o.l||(l==o.l&&r>o.r);}
    ll operator*(const tmp&o)const{return (ll)(o.r-l)*(r-o.l);}
}a[NM];
int n,q[NM],qh,qt,f[NM],v[NM];
ll d[NM],ans;


int main(){
    n=read();
    inc(i,1,n)a[i].l=read(),a[i].r=read();
    sort(a+1,a+1+n);
    inc(i,1,n)if(a[i].l<inf){
	for(int j=i+1;j<=n&&a[i].r>=a[j].r;j++)ans=max(ans,1ll*(a[i].r-a[i].l)*(a[j].r-a[j].l)),a[j].l=inf,a[j].r=inf;
    }
    sort(a+1,a+1+n);
    while(a[n].l==inf)n--;
    q[qh=qt=1]=1;f[1]=1;
    inc(i,2,n){
	f[i]=max(f[i],f[i-1]);
	while(qh<=qt&&f[i]!=q[qh])qh++;
	if(qh<=qt)d[i]=a[q[qh]]*a[i];
	int s=i+1;
	while(qh<=qt){
	    s=n+1;
	    for(int x=i+1,y=n;x<=y;)
		if(a[q[qt]]*a[mid]<a[i]*a[mid])s=mid,y=mid-1;else x=mid+1;
	    if(s<=v[q[qt]])qt--;else break;
	}
	if(s<=n)v[i]=s,f[s]=i,q[++qt]=i;
    }
    //inc(i,1,n)printf("%d ",f[i]);putchar('\n');
    //inc(i,1,n)printf("%lld ",d[i]);putchar('\n');
    inc(i,1,n)ans=max(ans,d[i]);
    return 0*printf("%lld\n",ans);
}

 

 

 

 

2369: 区间

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 175  Solved: 82
[Submit][Status][Discuss]

Description

对于一个区间集合

{A1,A2……Ak}(K>1,Ai不等于Aj(i不等于J),定义其权值

 

S=|A1∪A2∪……AK|*|A1∩A2……∩Ak|

即它们的交区间的长度乘上它们并区间的长度。

显然,如果这些区间没有交集则权值为0。

Your Task

给定你若干互不相等的区间,选出若干区间使其权值最大。

Input

第一行n表示区间的个数

接下来n行每行两个整数l r描述一个区间[l,r]

Output

 

在一行中输出最大权值

Sample Input

4
1 6
4 8
2 7
3 5

Sample Output

24

HINT

 



样例解释

    选择[1,6]和[2,7]是最优的。

数据约定

100%:1<N<=10^6,1<=L<R<=10^6

 

Source

 

[Submit][Status][Discuss]



  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值