1. 程式人生 > >bzoj2369/bzoj2687(貪心+單調決策性DP)

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]