题目描述
有N个村庄坐落在一条直线上，第i(i>1)个村庄距离第1个村庄的距离为Di。需要在这些村庄中建立不超过K个通讯基站，在第i个村庄建立基站的费用为Ci。如果在距离第i个村庄不超过Si的范围内建立了一个通讯基站，那么就村庄被基站覆盖了。如果第i个村庄没有被覆盖，则需要向他们补偿，费用为Wi。现在的问题是，选择基站的位置，使得总费用最小。
输入输出格式
输入格式：
输入文件的第一行包含两个整数N,K，含义如上所述。
第二行包含N-1个整数，分别表示D2,D3,…,DN ，这N-1个数是递增的。
第三行包含N个整数，表示C1,C2,…CN。
第四行包含N个整数，表示S1,S2,…,SN。
第五行包含N个整数，表示W1,W2,…,WN。
输出格式：
输出文件中仅包含一个整数，表示最小的总费用。
输入输出样例
输入样例#1：
3 2
1 2
2 3 2
1 1 0
10 20 30
输出样例#1：
4
说明
40%的数据中，N<=500；
100%的数据中，K<=N，K<=100，N<=20,000，Di<=1000000000，Ci<=10000，Si<=1000000000，Wi<=10000。
不会主席树+决策单调，otz
蒟蒻只能
dp+线段树优化
f[i][j]表示在第i个村庄修建第j个基站
首先对数据进行预处理详见代码有注释
用线段树维护区间最小值与区间增大时花费的增加
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=,maxk=;
struct A {
int l,r,w;
bool operator<(const A &a)const {
return r<a.r;
}
} a[maxn];
int mn[maxn<<],lazy[maxn<<],d[maxn],c[maxn],s[maxn],w[maxn];
int n,m,k,f[maxn][maxk],ans=;
void build(int l,int r,int rt)
{
lazy[rt]=;
if(l==r)
{
mn[rt]=f[l][k-];
return;
}
int mid=l+r>>;
build(l,mid,rt<<);
build(mid+,r,rt<<|);
mn[rt]=min(mn[rt<<],mn[rt<<|]);
}
void modify(int L,int R,int d,int l,int r,int rt)
{
if(L<=l&&R>=r) {
mn[rt]+=d;
lazy[rt]+=d;
return;
}
int mid=(l+r)>>;
if(L<=mid)modify(L,R,d,l,mid,rt<<);
if(R>mid)modify(L,R,d,mid+,r,rt<<|);
mn[rt]=min(mn[rt<<],mn[rt<<|])+lazy[rt];
}
int query(int L,int R,int l,int r,int rt)
{
if(L<=l&&R>=r) return mn[rt];
int mid=(l+r)>>;
int ans=0x7fffffff;
if(L<=mid)ans=min(query(L,R,l,mid,rt<<)+lazy[rt],ans);
if(R>mid) ans=min(ans,query(L,R,mid+,r,rt<<|)+lazy[rt]);
return ans;
}
//]f[i][j]表示在第i个村庄修建第j个基站
inline void init()
{
for(int i=;i<=n;i++)scanf("%d",&d[i]);
for(int i=;i<=n;i++)scanf("%d",&c[i]);
for(int i=;i<=n;i++)scanf("%d",&s[i]);
for(int i=;i<=n;i++)scanf("%d",&a[i].w);
for(int i=;i<=n;i++)
{
a[i].l=lower_bound(d+,d+n+,d[i]-s[i])-d;//可被覆盖范围，右
a[i].r=upper_bound(d+,d+n+,d[i]+s[i])-d-;//左
f[a[i].r+][]+=a[i].w;
w[a[i].l-]+=a[i].w;
ans+=a[i].w;//printf("%d %d %d\n",a[i].l,a[i].r,a[i].w);
}
}
int main() {
scanf("%d%d",&n,&m);
init();
sort(a+,a+n+);//对右端点排序
for(int i=;i<=n;i++)f[i][]+=f[i-][];// 预处理花费
for(int i=n;i;i--)w[i]+=w[i+];//第二位枚举到前i个基站，后边基站不被覆盖（不考虑）自然要加上这些花费
for(int i=;i<=n;i++)
{
f[i][]+=c[i];
ans=min(ans,f[i][]+w[i]);
}
f[][]=c[];
for(k=;k<=m;k++) //枚举建第几个基站
{
int cur=;
build(,n,);
for(int i=k;i<=n;i++)
{
while(cur<=n&&a[cur].r<i) // 计算建基站花费
{
if(a[cur].l>=k)
modify(k-,a[cur].l-,a[cur].w,,n,);
cur++;
}
f[i][k]=query(k-,i-,,n,)+c[i];
ans=min(ans,f[i][k]+w[i]);
}
}
printf("%d\n",ans);
return ;
}