题解 POJ3171 【Cleaning Shifts】

题目链接：Link

Problem

Solution

一个很直观的想法是设 $ f(x) $ 表示覆盖[L,x]需要花费的最小代价。
为避免后效性，可将奶牛右端点排序。设当前奶牛为 $ [L,R],c $ ，则状态转移方程为：
\[ f(R) = \min\limits_{L-1 \le x \le R} \{ f(x) \} + c \]
观察到这个dp过程事实上是在不断地在f数组里求区间最小值、单点修改，用线段树维护即可。

Code

#include<cstdio>
struct char_reader
{
    FILE* f; char *buf,*p1,*p2; int size;
    char_reader(FILE* fin,int bufsize=1000000) { f=fin; size=bufsize; p1=p2=0; buf=new char[bufsize]; }
    inline int operator()() { return p1==p2&&(p2=(p1=buf)+fread(buf,1,size,f),p1==p2)?EOF:*p1++; }
};
struct char_writer
{
    FILE* f; char *buf,*p,*end;
    char_writer(FILE* fout,int bufsize=1000000) { f=fout; buf=new char[bufsize]; p=buf; end=buf+bufsize; }
    ~char_writer() { fwrite(buf,p-buf,1,f); }
    inline char operator()(char ch) { return end==p&&(fwrite(buf,end-buf,1,f),p=buf),*p++=ch; }
};
char_reader gch(stdin);
char_writer wch(stdout);
template<typename T> inline void read(T& t)
{
    t=0; bool f=false; char ch;
    while(ch=gch(),!((ch>='0'&&ch<='9')||ch=='-'));
    if(ch=='-') f=true,ch=gch();
    t=ch-'0';
    while(ch=gch(),ch>='0'&&ch<='9') t=t*10+ch-'0';
    if(f) t=-t;
}
template<typename T> inline void write(T t)
{
    int stk[20],cnt=0;
    if(t==0) { wch('0'); return; }
    if(t<0) { wch('-'); t=-t; }
    while(t>0) { stk[cnt++]=t%10; t/=10; }
    while(cnt) wch(stk[--cnt]+'0');
}
inline void write(char t) { wch(t); }
inline void write(const char *&s) { while(*s) wch(*s++); }
#if __cplusplus >= 201103L
template<typename T,typename... Args> inline void read(T& t,Args&... args) { read(t); read(args...); }
template<typename T,typename... Args> inline void write(T t,Args... args) { write(t); write(args...); }
#else
template<typename A_t,typename B_t> inline void read(A_t &a,B_t &b) { read(a); read(b); }
template<typename A_t,typename B_t,typename C_t> inline void read(A_t &a,B_t &b,C_t &c) { read(a); read(b); read(c); }
template<typename A_t,typename B_t,typename C_t,typename D_t> inline void read(A_t &a,B_t &b,C_t &c,D_t &d) { read(a); read(b); read(c); read(d); }
template<typename A_t,typename B_t> inline void write(A_t a,B_t b) { write(a); write(b); }
template<typename A_t,typename B_t,typename C_t> inline void write(A_t a,B_t b,C_t c) { write(a); write(b); write(c); }
template<typename A_t,typename B_t,typename C_t,typename D_t> inline void write(A_t a,B_t b,C_t c,D_t d) { write(a); write(b); write(c); write(d); }
#endif
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=25005;
const int maxm=1000005;
const int oo=0x3f3f3f3f;
int n,L,R;
struct Node { int a,b,c; };
inline bool operator<(const Node &a,const Node &b) { return a.b<b.b; }
Node v[maxn];
int f[maxm<<2],pl,pr,val;
#define Lc (o<<1)
#define Rc (o<<1|1)
void build(int o,int L,int R)
{
    if(L==R) { f[o]=oo; return; }
    int M=(L+R)>>1;
    build(Lc,L,M); build(Rc,M+1,R);
    f[o]=min(f[Lc],f[Rc]);
}
int query(int o,int L,int R)
{
    if(pl<=L&&R<=pr) return f[o];
    int M=(L+R)>>1,res=oo;
    if(pl<=M) res=min(res,query(Lc,L,M));
    if(pr>M) res=min(res,query(Rc,M+1,R));
    return res;
}
void update(int o,int L,int R)
{
    if(L==R) { f[o]=val; return; }
    int M=(L+R)>>1;
    if(pr<=M) update(Lc,L,M);
    else update(Rc,M+1,R);
    f[o]=min(f[Lc],f[Rc]);
}
int main()
{
    #ifdef local
    freopen("pro.in","r",stdin);
    #endif
    read(n,L,R); L++; R++;
    for(int i=1;i<=n;i++) { read(v[i].a,v[i].b,v[i].c); v[i].a++; v[i].b++; }
    sort(v+1,v+1+n);
    build(1,L-1,R); val=0; pl=pr=L-1; update(1,L-1,R);
    for(int i=1;i<=n;i++)
    {
        pl=v[i].a-1; pr=v[i].b;
        val=query(1,L-1,R)+v[i].c;
        pl=pr=v[i].b;
        if(query(1,L-1,R)>val) update(1,L-1,R);
    }
    pl=pr=R;
    val=query(1,L-1,R);
    write(val!=oo?val:-1,'\n');
    return 0;
}

转载于:https://www.cnblogs.com/happyZYM/p/11392590.html