CF932F Escape Through Leaf(DP,斜率优化)

SB 题。

写出 DP 方程:\(f_i\) 表示从 \(i\) 跳的最小值。

\(i\) 是叶子就是 \(0\),否则就是选个子树中的 \(v\),\(f_i=\min(f_v+a_ib_v)\)。

至于优化,求出每个子树中的凸包就行了。启发式合并保证复杂度。

复杂度 \(O(n\log^2 n)\)。

没错,我又用了回家路线那又臭又长的写法。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=100010;
#define FOR(i,a,b) for(int i=(a);i<=(b);i++)
#define ROF(i,a,b) for(int i=(a);i>=(b);i--)
#define MEM(x,v) memset(x,v,sizeof(x))
inline int read(){
    int x=0,f=0;char ch=getchar();
    while(ch<'0' || ch>'9') f|=ch=='-',ch=getchar();
    while(ch>='0' && ch<='9') x=x*10+ch-'0',ch=getchar();
    return f?-x:x;
}
struct line{
    int k;
    ll b;
    bool operator<(const line &l)const{
        if(k!=l.k) return k>l.k;
        return b<l.b;
    }
};
struct point{
    double x;
    int k;
    ll b;
    bool operator<(const point &p)const{return x<p.x;} 
};
int n,el,head[maxn],to[maxn*2],nxt[maxn*2],A[maxn],B[maxn];
ll f[maxn];
set<line> inter[maxn];
set<point> hull[maxn];
inline void add(int u,int v){to[++el]=v;nxt[el]=head[u];head[u]=el;}
double interx(line l1,line l2){
    return l1.k==l2.k?1e18:1.0*(l2.b-l1.b)/(l1.k-l2.k);
}
void remove(int id,set<line>::iterator it){
    set<line>::iterator it1=it,it2=it;it2++;
    if(it1!=inter[id].begin()){
        it1--;
        hull[id].erase((point){interx(*it,*it1),it1->k,it1->b});
        it1++;
    }
    if(it2!=inter[id].end()) hull[id].erase((point){interx(*it,*it2),it->k,it->b});
    if(it1!=inter[id].begin() && it2!=inter[id].end()){
        it1--;
        hull[id].insert((point){interx(*it1,*it2),it1->k,it1->b});
    }
    inter[id].erase(it);
}
void insert(int id,line l){
    set<line>::iterator it=inter[id].insert(l).first;
    set<line>::iterator it1=it,it2=it;it2++;
    if(it1!=inter[id].begin()){
        it1--;
        if(it1->k==it->k) return void(inter[id].erase(*it));
        it1++;
    }
    if(it1!=inter[id].begin() && it2!=inter[id].end()){
        it1--;
        if(interx(*it,*it1)>=interx(*it,*it2)) return void(inter[id].erase(*it));
        it1++;
    }
    if(it1!=inter[id].begin()){
        it1--;
        hull[id].insert((point){interx(*it,*it1),it1->k,it1->b});
        it1++;
    }
    if(it2!=inter[id].end()) hull[id].insert((point){interx(*it,*it2),it->k,it->b});
    if(it1!=inter[id].begin() && it2!=inter[id].end()){
        it1--;
        hull[id].erase((point){interx(*it1,*it2),it1->k,it1->b});
        it1++;
    }
    it=it1=inter[id].find(l);
    while(it1!=inter[id].begin()){
        it1--;
        if(it1==inter[id].begin()) break;
        it2=it1;it2--;
        if(interx(*it2,*it)<=interx(*it2,*it1)) remove(id,it1);
        else break;
        it=it1=inter[id].find(l);
    }
    it=it1=inter[id].find(l);it1++;
    while(it1!=inter[id].end()){
        it2=it1;it2++;
        if(it2!=inter[id].end() && interx(*it2,*it)>=interx(*it2,*it1) || it1->k==it->k) remove(id,it1);
        else break;
        it=it1=inter[id].find(l);it1++;
    }
}
void dfs(int u,int F){
    for(int i=head[u];i;i=nxt[i]){
        int v=to[i];
        if(v==F) continue;
        dfs(v,u);
        if(inter[u].size()<inter[v].size()) inter[u].swap(inter[v]),hull[u].swap(hull[v]);
        for(set<line>::iterator it=inter[v].begin();it!=inter[v].end();it++) insert(u,*it);
        inter[v].clear();hull[v].clear();
    }
    if(!inter[u].empty()){
        set<point>::iterator it=hull[u].lower_bound((point){A[u],0,0});
        int k;
        ll b;
        if(it==hull[u].end()){
            set<line>::iterator it=inter[u].end();it--;
            k=it->k;b=it->b;
        }
        else k=it->k,b=it->b;
        f[u]=1ll*k*A[u]+b;
    }
    insert(u,(line){B[u],f[u]});
}
int main(){
    n=read();
    FOR(i,1,n) A[i]=read();
    FOR(i,1,n) B[i]=read();
    FOR(i,1,n-1){
        int u=read(),v=read();
        add(u,v);add(v,u);
    }
    dfs(1,0);
    FOR(i,1,n) printf("%lld ",f[i]);
}
上一篇:Leaf for Mac (rss阅读器)


下一篇:两种开源分布式主键生成器简介