POJ 3171 区间最小花费覆盖 （DP+线段树

Cleaning Shifts

Time Limit: 1000MS		Memory Limit: 65536K
Total Submissions: 4245		Accepted: 1429

Description

Farmer John's cows, pampered since birth, have reached new heights of fastidiousness. They now require their barn to be immaculate. Farmer John, the most obliging of farmers, has no choice but hire some of the cows to clean the barn. 

Farmer John has N (1 <= N <= 10,000) cows who are willing to do some cleaning. Because dust falls continuously, the cows require that the farm be continuously cleaned during the workday, which runs from second number M to second number E during the day (0 <= M <= E <= 86,399). Note that the total number of seconds during which cleaning is to take place is E-M+1. During any given second M..E, at least one cow must be cleaning. 

Each cow has submitted a job application indicating her willingness to work during a certain interval T1..T2 (where M <= T1 <= T2 <= E) for a certain salary of S (where 0 <= S <= 500,000). Note that a cow who indicated the interval 10..20 would work for 11 seconds, not 10. Farmer John must either accept or reject each individual application; he may NOT ask a cow to work only a fraction of the time it indicated and receive a corresponding fraction of the salary. 

Find a schedule in which every second of the workday is covered by at least one cow and which minimizes the total salary that goes to the cows.

Input

Line 1: Three space-separated integers: N, M, and E. 

Lines 2..N+1: Line i+1 describes cow i's schedule with three space-separated integers: T1, T2, and S.

Output

Line 1: a single integer that is either the minimum total salary to get the barn cleaned or else -1 if it is impossible to clean the barn.

Sample Input

3 0 4
0 2 3
3 4 2
0 0 1

Sample Output

5

Hint

Explanation of the sample: 

FJ has three cows, and the barn needs to be cleaned from second 0 to second 4. The first cow is willing to work during seconds 0, 1, and 2 for a total salary of 3, etc. 

Farmer John can hire the first two cows.

Source

USACO 2005 December Silver

 

题意：给n个区间及其代价值，问要覆盖[M,E]区间至少要花费多少代价；

解法：这是一个dp问题，先列出方程。

F[i]表示取[0,i]这个区间的代价,初始化F[M-1]=0,答案就是F[E].

则方程为F[a[i].T2]=min(F[a[j].T2])+a[i].s (T1-1<=a[j].T2<T2),找min的过程用线段树实现。

将a[i]按T2从小到大排列,逐步更新最小值。

代码：

#include"bits/stdc++.h"

#define ll long long
#define vl vector<ll>
#define ci(x) scanf("%d",&x)
#define pi(x) printf("%d
",x)
#define pl(x) printf("%lld
",x)
#define rep(i, n) for(int i=0;i<n;i++)
using namespace std;
const int NN = 1e6 + 5;
int n,s,t;
struct P{int x,y,s;};
P a[NN];
bool cmp(P a,P b){
    return a.y<b.y;
}
const ll INF = 0x3fffffffffffffff;
struct SegMin {
    int N;
    vl is;vl mul;vl add;
    ll init;
    ll merge(ll a, ll b) {
        return min(a, b);
    }
    void push(int o, int L, int R, ll m, ll a) {
        is[o] = is[o] * m + a;
        mul[o] = mul[o] * m;
        add[o] = add[o] * m + a;
    }

    SegMin(int n, ll init=INF) {
        N = 1;
        while (N < n) N *= 2;
        this->init = init;
        is = vl(N * 4, init);
        mul = vl(N * 4, 1);
        add = vl(N * 4);
    }

    SegMin(vl a, ll init=INF) {
        int n = a.size();
        N = 1;
        while (N < n) N *= 2;
        this->init = init;
        is = vl(N * 2);
        mul = vl(N * 2, 1);
        add = vl(N * 2);
        copy(a.begin(), a.end(), is.begin() + N);
        for (int i = N - 1; i > 0; i--) {
            is[i] = merge(is[i << 1], is[i << 1 | 1]);
        }
    }

    void update(int l, int r, ll m, ll a) {
        if (l < r) update(1, 0, N, l, r, m, a);
    }

    void update(int o, int L, int R, int l, int r, ll m, ll a) {
        if (l <= L && R <= r) {
            push(o, L, R, m, a);
        } else {
            int M = (L + R) >> 1;
            push(o, L, M, R);
            if (l < M) update(o << 1, L, M, l, r, m, a);
            if (r > M) update(o << 1 | 1, M, R, l, r, m, a);
            is[o] = merge(is[o << 1], is[o << 1 | 1]);
        }
    }

    void push(int o, int L, int M, int R) {
        if (mul[o] != 1 || add[o] != 0) {
            push(o << 1, L, M, mul[o], add[o]);
            push(o << 1 | 1, M, R, mul[o], add[o]);
            mul[o] = 1;
            add[o] = 0;
        }
    }

    ll query(int l, int r) {
        if (l < r) return query(1, 0, N, l, r);
        return init;
    }

    ll query(int o, int L, int R, int l, int r) {
        if (l <= L && R <= r) {
            return is[o];
        } else {
            int M = (L + R) >> 1;
            push(o, L, M, R);
            ll res = init;
            if (l < M) res = merge(res, query(o << 1, L, M, l, r));
            if (r > M) res = merge(res, query(o << 1 | 1, M, R, l, r));
            is[o] = merge(is[o << 1], is[o << 1 | 1]);
            return res;
        }
    }
};

int main(){
    ci(n),ci(s),ci(t);//s从1开始
    s++,t++;
    int ma=0;
    for(int i=0;i<n;i++) ci(a[i].x),ci(a[i].y),ci(a[i].s);
    for(int i=0;i<n;i++) a[i].x++,a[i].y++,ma=max(ma,a[i].y);
    sort(a,a+n,cmp);
    SegMin seg(ma+1);
    seg.update(0,ma+1,0,INF);
    seg.update(0,s,0,0);

    for(int i=0;i<n;i++){
        if(a[i].y<s) continue;
        int L=a[i].x-1,R=a[i].y;
        ll res=seg.query(L,R)+a[i].s;
        res=min(seg.query(R,R+1),res);//与前面的最小值取min
        seg.update(R,R+1,0,res);
    }
    ll ans=seg.query(t,ma+1);
    if(ans>=INF) puts("-1");//未覆盖到
    else pl(ans);
    return 0;
}

转载于:https://www.cnblogs.com/mj-liylho/p/9502012.html