题目链接:https://nanti.jisuanke.com/t/26156
题目大意:给你一个n个点的图,图上只存在以下两种双向边:
1、E{1,i}(1<i<=n)
2、E{i,i+1}(1<i<n)
题目已经给出以上两种边的所有边长,接下来有q次询问,每次询问给出一个d,使得图中所有的边的权值都增加d,但增加的边长是一次性的(也就是每次询问的d对接下来的询问是没有影响的),每次询问要你输出点1到点n的最短路径。
题目思路:由题目所给出的条件我们可以对每次询问都推出如下的式子:
ans=min{sum[n]-sum[j]+e1[j]+(n-j+1)*d}(2<=j<=n)
其中sum[i]表示第二种边的边权的前缀和(sum[n]-sum[j]就表示点 j 到点 n的距离)
e1[i]表示点1到 i 这个的点直连边的边长;
接下来我们进行假设,对于k<j<=n,有点1选择从点 j 走到点n的的路径要小于点1选择从点 k 走到点n的的路径,
我们可以表达出如下的式子:
sum[n]-sum[j]+e1[j]+(n-j+1)*d < sum[n]-sum[k]+e1[k]+(n-k+1)*d
移项合并同类项之后得到
((sum[j]-e1[j])-(sum[k]-e1[k]))/(j-k) > -d;
令yj=sum[j]-e1[j] , yk=sum[k]-e1[k] , xj=j , xk = k;
那么上式就可以化为 (yj-yk)/(xj-xk) > -d,左边明显就是一个斜率公式,我们就可以用斜率优化的思想来解决这个问题了。
首先我们可以先将没有加d的原始边进行预处理,构造一个上凸包,
接着再将查询离线化,令每个d=-d,再按d的大小从大到小进行排序,最后再维护上凸包即可求得答案了。
具体实现看代码
#include <bits/stdc++.h>
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define fi first
#define se second
#define pb push_back
#define MP make_pair
#define FIN freopen("in.txt","r",stdin)
using namespace std;
typedef long long LL;
typedef pair<int,int>pii;
typedef pair<LL,LL>pll;
typedef double db;
const int MX=1e5+7;
struct point{
LL x,y;
point(){}
point(LL _x,LL _y){
x=_x;y=_y;
}
}q[MX],p[MX];
int n,m;
LL e1[MX],e2[MX],d[MX],ans[MX],sum[MX];
int id[MX];
bool cmp(int a,int b){
return d[a]>d[b];
}
double get_k(point a,point b){
return (double)(a.y-b.y)/(a.x-b.x);
}
int main(){
int T;cin>>T;
while(T--){
scanf("%d%d",&n,&m);
for(int i=2;i<=n;i++) scanf("%lld",&e1[i]);
sum[2]=0;
for(int i=2;i<n;i++){
scanf("%lld",&e2[i]);
sum[i+1]=sum[i]+e2[i];
}
for(int i=2;i<=n;i++) sum[i]-=e1[i];
int head=1,tail=0;
p[0]=point();
for(int i=2;i<=n;i++){
p[i]=point(i,sum[i]);
while(tail>=2 && get_k(p[i],q[tail])>get_k(q[tail],q[tail-1])) tail--;
q[++tail]=p[i];
}
for(int i=1;i<=m;i++){
scanf("%lld",&d[i]);
d[i]=-d[i];
id[i]=i;
}
sort(id+1,id+m+1,cmp);
for(int i=1;i<=m;i++){
LL D=d[id[i]];
while(head<tail && D<get_k(q[head+1],q[head])) head++;
ans[id[i]]=sum[n]-sum[q[head].x]+e1[n]+(n-q[head].x+1)*(-D);
}
for(int i=1;i<=m;i++)
printf("%lld%c",ans[i],i==m?'\n':' ');
}
return 0;
}