美梦1（JSOI2014SC)

【问题描述】

　　这天晚上，约翰做了个奇怪的美梦。他拥有了分别分布在N座高高低低的山上的N个池塘，N座山连成一条直线，从左往右第i座山的高度是Hi。池塘中的鱼都是他请专家运用科学的方法专门养殖的，为了保护每个池塘的生态环境，他现在要在这N座山上建造若干个看护点。约翰是个很节约的人，在第i座山建造看护点的花费为Ci。假设在第i座山建造一个看护点，则往左或者往右第一座不比这座山低的山将挡住看护的视线。譬如说： {Hi} = {1 4 4 5 7 2}表示第一座山高度为1，第二座山高度为4。。。 如果在第1座山建造一个看护点，则可以看护第1,2两个池塘。如果在第5座山上建造一个看护点，则左右的池塘都能被看护到。如果在第3座山上建造一个看护点，则能够看护到第2，3，4个池塘。

Problem Task:

　　要求能够看护到所有的池塘，建造看护点的最小代价是多少。即建造看护点的山对应的花费Ci之和最小。

Problem Input:

　　第一行包含一个正整数N，N满足1<=N<=1000000。

　　第二行包含N个正整数，第i个正整数Hi满足1<=Hi<=10^9，表示第i座山的高度

　　第三行包含N个正整数，第i个正整数Ci满足1<=Ci<=10^9，表示在第i座山建造看护点的代价为Ci

Problem Output:

　　一行包含一个正整数C，表示最小的代价。

Problem Input Example:

3

1 1 1

2 2 2

Problem output Example:

2

这题是个很裸的动归  ， 不过数据范围有些怕人... 

Way1：设f[i] 表示 看护1~i座山所需最小花费，则显然有f[i] = min{f[k] + C_t} (Lt<=k<=Rt = i)

不难发现可以用线段树维护最优值，单点修改，区间查询。

Way2: 设f[i] 表示 看护1~i座山，且第i座山上修建看护点 ，则显然有 f[i] = min{f[j]} + C_i (R_j >= L_i )

不难发现也可以偶那个线段树维护，单点查询，区间修改

Codes:

#include<set>
#include<queue>
#include<cstdio>
#include<vector>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
using namespace std;
const int N = 1000010;
typedef long long lld;
const lld inf = lld(21474836470000000);
#define Ch1 (i<<1)
#define Ch2 (Ch1|1)
#define For(i,n) for(int i=1;i<=n;i++)
#define Rep(i,l,r) for(int i=l;i<=r;i++)
#define Down(i,r,l) for(int i=r;i>=l;i--)
struct tnode{
    int l,r,mid;
    lld min;
}T[N<<2];
int Left[N],Right[N],top,n,H[N],C[N],st[N];
lld f[N];

void read(int &v){
    int num = 0; char ch = getchar();
    while(ch>'9'||ch<'0') ch = getchar();
    while(ch>='0'&&ch<='9'){
        num = num * 10 + ch - '0';
        ch = getchar();
    }
    v = num;
}

void init(){
    read(n);
    For(i,n) read(H[i]);
    For(i,n) read(C[i]);
    st[top=1] = 0;H[0] = H[n+1] = 2147483647;
    For(i,n){
        while(H[st[top]]<H[i]) top--;
        Left[i] = max(st[top],1);
        st[++top] = i;
    }
    st[top=1] = n+1;
    Down(i,n,1){
        while(H[st[top]]<H[i]) top--;
        Right[i] = min(st[top],n);
        st[++top] = i;
    }
}

void Build(int l,int r,int i){
    T[i].l = l;T[i].r = r;T[i].mid = (l+r)>>1;
    T[i].min = inf;
    if(l==r) return;
    Build(l,T[i].mid,Ch1); Build(T[i].mid+1,r,Ch2);
}
vector< pair<int,int> > G[N];
//f[i] = min{f[t] + cost[k]}

void Modify(int i,int x,lld delta){
    if(T[i].l==T[i].r){
        T[i].min = delta;
        return;
    }
    if(x<=T[i].mid)  Modify(Ch1,x,delta);
    else             Modify(Ch2,x,delta);
    T[i].min = min(T[Ch1].min,T[Ch2].min);
}

lld query(int l,int r,int i){
    if(l<=T[i].l&&T[i].r<=r) return T[i].min;
    if(r<=T[i].mid)  return query(l,r,Ch1);
    if(l>T[i].mid)   return query(l,r,Ch2);
    return min(query(l,T[i].mid,Ch1) , query(T[i].mid+1,r,Ch2));
}

int main(){
    freopen("fish1.in","r",stdin);
    freopen("fish1.out","w",stdout);
    init();
    Build(0,n,1);
    Modify(1,0,0);
    For(i,n){
        G[Right[i]].push_back(make_pair(Left[i],C[i]));
        G[i].push_back(make_pair(Left[i],C[i]));
    }
    For(i,n){
        f[i] = inf;
        for(int j = 0;j<G[i].size();j++) 
           f[i] = min(f[i],query(G[i][j].first-1,i,1) + G[i][j].second);
        Modify(1,i,f[i]);
    }
    cout<<f[n]<<endl;
    return 0;
}