【LOJ576】「LibreOJ NOI Round #2」签到游戏（最小生成树）（线段树）

传送门

---

题解：

首先这是一个非常常见以及SB的生成树模型。

我们把一个端点固定的边的代价函数拎出来看一看。发现是从中间向两边单调不增，显然是gcd的性质。

那么很显然地，只会有两种边$(1,i),(i,n)$。

再想一想发现这两种边的代价函数一定在中间会交起来，最小生成树一定存在一个$p$，使得所有边都是$[1,i](i>p),[i,n](i\le p),[1,n]$

线段树维护区间gcd，然后线段树上二分找到这个$p$，再线段树上找到前缀$gcd$和后缀gcd变化的$log$个位置算一下答案就行了。

---

官方题解的一点小锅：

其实官方题解里面这句话可能是写的时候没想清楚：

直接暴力gcd复杂度是对的。

我们把取模后值改变的取模称为有效取模，否则为无效取模。

显然无效取模只可能在取$gcd$的第一次执行，即无效取模总数量为$O(\log n\log B)$

而每一次有效取模会导致值变为原来的一半，而线段树一共会递归到最多$O(\log B)$个叶子。所以均摊分析一下，有效取模的数量最多为$O({\log }^{2}B)$。

然后剩下的部分求$gcd$本身都是$O(1)$，不存在什么复杂度的问题。

这个故事告诉我们均摊分析是多么的重要

---

代码：

#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const

namespace IO{
	inline char gc(){
		static cs int Rlen=1<<22|1;
		static char buf[Rlen],*p1,*p2;
		return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,Rlen,stdin),p1==p2)?EOF:*p1++;
	}
	
	template<typename T>
	inline T get(){
		char c;T num;
		while(!isdigit(c=gc()));num=c^48;
		while(isdigit(c=gc()))num=(num+(num<<2)<<1)+(c^48);
		return num;
	}
	inline int gi(){return get<int>();}
}
using namespace IO;

using std::cerr;
using std::cout;

cs int N=1e5+7;

int g[N<<2];

#define lc u<<1
#define rc u<<1|1
inline void build(int u,int l,int r){
	if(l==r){g[u]=gi();return ;}int mid=l+r>>1;
	build(lc,l,mid);build(rc,mid+1,r);g[u]=std::__gcd(g[lc],g[rc]);
}

inline void modify(int u,int l,int r,int p){
	if(l==r){g[u]=gi();return ;}int mid=l+r>>1;
	(p<=mid)?modify(lc,l,mid,p):modify(rc,mid+1,r,p);g[u]=std::__gcd(g[lc],g[rc]);
}

inline int find(int u,int l,int r,int a,int b){
	if(l==r)return l;int mid=l+r>>1;
	int pre=std::__gcd(a,g[lc]),suf=std::__gcd(b,g[rc]);
	return (pre<=suf)?find(lc,l,mid,a,suf):find(rc,mid+1,r,pre,b); 
}

inline ll get_pre(int u,int l,int r,int p,int val){
	if(l==r)return std::__gcd(val,g[u]);
	int mid=l+r>>1,pre=std::__gcd(val,g[rc]);
	if(p<=mid)return get_pre(lc,l,mid,p,pre);
	else if(g[lc]%pre==0)return (mid-l+1ll)*pre+get_pre(rc,mid+1,r,p,val);
	return get_pre(lc,l,mid,p,pre)+get_pre(rc,mid+1,r,p,val);
}

inline ll get_suf(int u,int l,int r,int p,int val){
	if(l==r)return std::__gcd(val,g[u]);
	int mid=l+r>>1,suf=std::__gcd(val,g[lc]);
	if(mid<p)return get_suf(rc,mid+1,r,p,suf);
	else if(g[rc]%suf==0)return (r-mid+0ll)*suf+get_suf(lc,l,mid,p,val);
	return get_suf(rc,mid+1,r,p,suf)+get_suf(lc,l,mid,p,val);
}

int n,q;
signed main(){
#ifdef zxyoi
	freopen("sign.in","r",stdin);
#endif
	n=gi(),q=gi();build(1,1,n);
	while(q--){
		modify(1,1,n,gi());int p=find(1,1,n,0,0);
		cout<<get_pre(1,1,n,p,0)+get_suf(1,1,n,p,0)-g[1]<<"\n";
	}
	return 0;
}