bzoj 1835: [ZJOI2010]base 基站选址

题意：

用最少代价，按题目要求覆盖所有点。

题解：

不怎么会做，orz状态表示，知道$f[i][j]$表示前i个，建j个塔i号点必建并且全部覆盖的最小代价。
因为只考虑i前的，所以统计答案就可以多加一个无限远的点，转移到那里就是答案了。
于是有$f[i][j]=min(f[k][j-1]+cost(k,i))+c[i]$
$cost(k,i)$是k到i中需要补偿的点。
先上暴力dp的代码：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#define LL long long
using namespace std;
const int inf=(1<<28);
int f[20010][110],n,k,sum[20010],s[20010],c[20010],w[20010];
void cmin(int &x,int y)
{
    if(x==-1) x=y;
    x=min(x,y);
}
int cost(int l,int r)
{
    int ans=0;
    for(int i=l+1;i<r;i++)
        if(sum[i]+s[i]<sum[r]&&(l==0||sum[l]<sum[i]-s[i])) ans+=w[i];
    return ans;
}
int main()
{
    scanf("%d %d",&n,&k);
    sum[1]=0;
    for(int i=2;i<=n;i++) scanf("%d",&sum[i]);//距离 
    for(int i=1;i<=n;i++) scanf("%d",&c[i]);//建基站费用 
    for(int i=1;i<=n;i++) scanf("%d",&s[i]);//接受范围 
    for(int i=1;i<=n;i++) scanf("%d",&w[i]);//补偿费用 
    n++;sum[n]=inf;c[n]=0;k++;
    memset(f,-1,sizeof(f));
    f[0][0]=0;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=min(i,k);j++)
        {
            for(int k=0;k<i;k++)
                if(f[k][j-1]!=-1) cmin(f[i][j],f[k][j-1]+cost(k,i));
            f[i][j]+=c[i];
        }
    int ans=-1;
    for(int i=0;i<=k;i++) cmin(ans,f[n][i]);
    printf("%d",ans);
}

然后可以发现能用线段树优化。
还要开个堆，维护当前接受范围最短的点。
然后二分出不能接受的点的影响范围，线段树随便搞搞用好了。
口胡能力不足，具体代码。
各种SB错误荒废一晚上。
code:

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<queue>
#include<algorithm>
#define LL long long
using namespace std;
struct node{
    LL x,i;
    bool operator < (node a) const{return a.x<x;}
};
priority_queue<node> b[110];
const LL inf=(1<<29);
LL f[20010][110],n,m,sum[20010],s[20010],c[20010],w[20010];
struct trnode{
    LL lc,rc,c;
    LL u;
}tr[40010];LL tot=0;
void cmin(LL &x,LL y)
{
    if(x==-1) x=y;
    x=min(x,y);
}
void update(LL x)
{
    LL c=tr[x].u,lc=tr[x].lc,rc=tr[x].rc;
    //if(!lc) lc=tr[x].lc=++tot;
    //if(!rc) rc=tr[x].rc=++tot;
    tr[lc].c+=c;tr[rc].c+=c;
    tr[lc].u+=c;tr[rc].u+=c;
    tr[x].u=0;
}
void change(LL x,LL l,LL r,LL k,LL c)
{
    //if(!x) x=++tot;
    if(l==r){tr[x].c=c;return;}
    LL mid=(l+r)/2;
    if(tr[x].u!=0) update(x);
    if(k<=mid) change(tr[x].lc,l,mid,k,c);
    else change(tr[x].rc,mid+1,r,k,c);
    tr[x].c=min(tr[tr[x].lc].c,tr[tr[x].rc].c);
}
void add(LL x,LL l,LL r,LL fl,LL fr,LL c)
{
    //if(!x) x=++tot;
    if(l==fl&&r==fr)
    {
        tr[x].u+=c;
        tr[x].c+=c;
        return;
    }
    if(tr[x].u!=0) update(x);
    LL mid=(l+r)/2;
    if(fr<=mid) add(tr[x].lc,l,mid,fl,fr,c);
    else if(fl>mid) add(tr[x].rc,mid+1,r,fl,fr,c);
    else add(tr[x].lc,l,mid,fl,mid,c),add(tr[x].rc,mid+1,r,mid+1,fr,c);
    tr[x].c=min(tr[tr[x].lc].c,tr[tr[x].rc].c);
}
LL findans(LL x,LL l,LL r,LL fl,LL fr)
{
    //if(!x) return inf;
    if(l==fl&&r==fr) return tr[x].c;
    LL mid=(l+r)/2;
    if(tr[x].u!=0) update(x);
    if(fr<=mid) return findans(tr[x].lc,l,mid,fl,fr);
    if(fl>mid) return findans(tr[x].rc,mid+1,r,fl,fr);
    return min(findans(tr[x].lc,l,mid,fl,mid),findans(tr[x].rc,mid+1,r,mid+1,fr));
}
void solve(LL L,LL i)
{
    if(sum[L]>=sum[i]-s[i]) return;
    LL l=L,r=i,ans;
    while(l<=r)
    {
        LL mid=(l+r)/2;
        if(sum[mid]<sum[i]-s[i]) l=mid+1,ans=mid;
        else r=mid-1;
    }
    add(1,1,n,L,ans,w[i]);
}
LL bt(LL l,LL r)
{
    //printf("%lld %lld %lld\n",tot,l,r);
    LL x=++tot;
    tr[x].lc=tr[x].rc=0;tr[x].c=inf;
    tr[x].u=0;
    if(l!=r)
    {
        LL mid=(l+r)/2;
        tr[x].lc=bt(l,mid);
        tr[x].rc=bt(mid+1,r);
    }
    return x;
}
int main()
{
    //freopen("a.in","r",stdin);
    scanf("%lld %lld",&n,&m);
    sum[1]=0;
    for(LL i=2;i<=n;i++) scanf("%lld",&sum[i]);//距离 
    for(LL i=1;i<=n;i++) scanf("%lld",&c[i]);//建基站费用 
    for(LL i=1;i<=n;i++) scanf("%lld",&s[i]);//接受范围 
    for(LL i=1;i<=n;i++) scanf("%lld",&w[i]);//补偿费用 
    n++;sum[n]=inf;c[n]=0;m++;sum[0]=-inf;c[0]=0;
    memset(f,-1,sizeof(f));
    tot=0;
    for(LL i=0;i<=n;i++) tot+=w[i],f[i][0]=tot;
    while(!b[1].empty()) b[1].pop();
    tot=0;
    for(LL i=1;i<=n;i++)
    {
        while(!b[1].empty()&&b[1].top().x<sum[i]) tot+=w[b[1].top().i],b[1].pop();//printf("\n");
        f[i][1]=tot+c[i];
        node tmp;tmp.x=sum[i]+s[i];tmp.i=i;
        b[1].push(tmp);
    }
    for(LL j=2;j<=m;j++)
    {
        //memset(tr,0,sizeof(tr));
        tot=0;bt(1,n);
        change(1,1,n,j-1,f[j-1][j-1]);
        while(!b[j].empty()) b[j].pop();
        for(LL i=j;i<=n;i++)
        {
            while(!b[j].empty()&&b[j].top().x<sum[i]) solve(j-1,b[j].top().i),b[j].pop();//printf("\n");
            f[i][j]=findans(1,1,n,j-1,i-1)+c[i];
            change(1,1,n,i,f[i][j-1]);
            node tmp;tmp.x=sum[i]+s[i];tmp.i=i;
            b[j].push(tmp);
        }
    }
    LL ans=-1;
    for(LL i=0;i<=m;i++) cmin(ans,f[n][i]);//,printf("%lld ",f[n][i]);printf("\n");
    printf("%lld",ans);
}