Jzoj3780 Magical GCD

给你一个序列s,求区间[l,r]使得gcd(s[l]~s[r])*(r-l+1)最大,输出这个最大值,n<=200000,cases<=20

显然,我们暴力是不可以过的,那么我们考虑一下的做法:

我们维护一个g数组,g[i]表示gcd(s[i]~s[r])

那么我们每次移动r指针,让后重新更新g数组,显然每个g[i]可以O(lg k)更新,但是这样的复杂度依然是n^2lgk

我们发现,对于g[i]=g[j]而且i<j,我们可以把j扔掉不管,这样的话可以将复杂度优化到nlgklgk

为什么?因为显然,g数组是单调的,而且g[i]|g[i+1],那么这就导致g[i+1]>=2*g[i],所以整个g数组长度不会超过lg k

#pragma GCC opitmize("O3")
#pragma G++ opitmize("O3")
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define L long long
using namespace std;
inline L gcd(L a,L b){
	for(L c;b;a=b,b=c) c=a%b;
	return a;
}
struct Num{ int id; L x; } g[100];
int n,t=0; L A=0,x;
int F(){
	scanf("%d",&n); t=A=0;
	*g=(Num){0,0};
	for(int i=1,k;i<=n;++i){
		scanf("%lld",&x);
		g[++t]=(Num){i,x};
		for(int j=k=1;j<=t;++j){
			g[j].x=gcd(g[j].x,x); 
			A=max(A,g[j].x*(i-g[j].id+1));
			if(g[j].x!=g[j-1].x) g[k++]=g[j];
		}
		t=--k;
	}
	printf("%lld\n",A);
}
int main(){
	int T; 
	for(scanf("%d",&T);T--;F());
}


posted @ 2017-09-28 16:26  扩展的灰(Extended_Ash)  阅读(97)  评论(0编辑  收藏  举报