CSP-S模拟赛20190921 T3 御神渡【分治求异色点之间的凸包】

题目描述：

$N\le 10^5$

题目分析：

边$(i,j)$的权值是$C_i+C_j-A_i\ast A_j$，本质上是最小生成树问题。但是$n\le 10^5$。

不难用反证法证明一个集合向外的所有边中最小的边一定在最小生成树中，基于这个基本性质，我们考虑与一个点相连的所有边中的最小值。

令$z=C_i+C_j-A_i\ast A_j$，对于点$i$来说，$C_j=A_i\ast A_j+z-C_i$，可看做$y=A_{i}x+z-C_i$，要截距最小，容易用斜率优化解决最小值问题。

此时有一个大致的想法就是每次找每个点的最小边，加入最小生成树，连通块缩成一个点，直到只剩一个点

出题人题解：

这个分治非常地强 XD

Code：

#include<cstdio>
#include<cmath>
#include<algorithm>
#define LL long long
#define maxn 100005
using namespace std;
int n,m,scc[maxn],N,res[maxn],id[maxn],A[maxn],C[maxn],f[maxn];
LL ans;
int q[maxn],t,tmp[maxn];
bool cmp(int i,int j){return scc[i]<scc[j];}
inline int find(int x){return !f[x]?x:f[x]=find(f[x]);}
inline int LL calc(int i,int j){return C[i]+C[j]-1ll*A[i]*A[j];}
void solve(int lc,int rc,int l,int r){
	if(lc==rc) return;
	int mc=(lc+rc)>>1,mid;
	for(int i=r;;i--) if(scc[id[i]]==mc) {mid=i;break;}
	solve(lc,mc,l,mid),solve(mc+1,rc,mid+1,r);
	t=0;
	for(int i=l;i<=mid;q[++t]=id[i],i++)
		while(t>1&&1ll*(A[q[t]]-A[q[t-1]])*(C[id[i]]-C[q[t-1]])<=1ll*(A[id[i]]-A[q[t-1]])*(C[q[t]]-C[q[t-1]])) t--;
	for(int i=mid+1,h=1;i<=r;i++){
		while(h<t&&C[q[h+1]]-C[q[h]]<=1ll*(A[q[h+1]]-A[q[h]])*A[id[i]]) h++;
		if(!res[id[i]]||calc(id[i],res[id[i]])>calc(id[i],q[h])) res[id[i]]=q[h];
	}
	t=0;
	for(int i=mid+1;i<=r;q[++t]=id[i],i++)
		while(t>1&&1ll*(A[q[t]]-A[q[t-1]])*(C[id[i]]-C[q[t-1]])<=1ll*(A[id[i]]-A[q[t-1]])*(C[q[t]]-C[q[t-1]])) t--;
	for(int i=l,h=1;i<=mid;i++){
		while(h<t&&C[q[h+1]]-C[q[h]]<=1ll*(A[q[h+1]]-A[q[h]])*A[id[i]]) h++;
		if(!res[id[i]]||calc(id[i],res[id[i]])>calc(id[i],q[h])) res[id[i]]=q[h];
	}
	for(int i=l,j=mid+1,k=l;k<=r;k++)
		tmp[k]=i>mid?id[j++]:j>r?id[i++]:A[id[i]]<A[id[j]]?id[i++]:id[j++];
	for(int i=l;i<=r;i++) id[i]=tmp[i];
}
int main()
{
	//freopen("gokami.in","r",stdin);
	//freopen("gokami.out","w",stdout);
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%d",&A[i]);
	for(int i=1;i<=n;i++) scanf("%d",&C[i]);
	while(1){
		N=0;
		for(int i=1;i<=n;i++) if(!f[i]) scc[i]=++N;
		if(N==1) break;
		for(int i=1;i<=n;i++) if(f[i]) scc[i]=scc[find(f[i])];
		for(int i=1;i<=n;i++) id[i]=i,res[i]=0;
		sort(id+1,id+1+n,cmp);
		solve(1,N,1,n);
		for(int i=1;i<=n;i++) if(find(i)!=find(res[i])) ans+=calc(i,res[i]),f[find(i)]=find(res[i]);
	}
	printf("%lld\n",ans);
}