【NOIP 校内模拟】T1 优美的序列(二分+st表+卡常)

我是菜鸡 我是蒟蒻 我好菜

ak一定是区间最小的值,且是所有数(包括自己)的最大公约数

我没看出来 没救了 noip爆零了 回家养猪了 没学上了 怎么办

gcd有单调性 gcd有单调性 gcd有单调性 gcd有单调性 gcd有单调性 可以二分 可以二分 可以二分 可以二分 要说多少遍才记得到

我没想到 没救了 noip爆零了 回家养猪了 没学上了 怎么办

stl少用 stl少用 stl少用 stl少用 要说多少遍

我忘了 没救了 noip爆零了 回家养猪了 没学上了 怎么办

#include<bits/stdc++.h>
#define N 500005
#define re register 
using namespace std;
int n,a[N],gcd[N][21],maxd;
vector<int> ans;
inline int calc(int a,int b)
{
	re int t;
	while (b)
		t=a,a=b,b=t%b;	
	return a;
}
int LOG[N];
void init()
{
	LOG[0]=-1; 
	for(re int i=1;i<=n;++i)	gcd[i][0]=a[i],LOG[i]=LOG[i/2]+1;
	for(re int j=1;j<=21;++j)
	{
		for(re int i=1;i+(1<<j)-1<=n;++i)
			gcd[i][j]=calc(gcd[i][j-1],gcd[i+(1<<(j-1))][j-1]);
	}
}
inline int query(int l,int r)
{
	int k=LOG[r-l+1];
	return calc(gcd[l][k],gcd[r-(1<<k)+1][k]);
}
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(NULL);cout.tie(NULL);
	cin>>n;
	for(re int i=1;i<=n;++i)	cin>>a[i];
	init();
	for(re int i=1;i<=n;++i)	//枚举最小的 
	{
		re int l=1,r=i,mid;
		re int ans1=i,ans2=i;
		while(l<=r)
		{
			mid=(l+r)>>1;
			if(query(mid,i)==a[i])	ans1=mid,r=mid-1;
			else l=mid+1;
		}
		l=i,r=n;
		while(l<=r)
		{
			mid=(l+r)>>1;
			if(query(i,mid)==a[i])	ans2=mid,l=mid+1;
			else r=mid-1;
		}
		if(ans2-ans1==maxd)	ans.push_back(ans1);
		if(ans2-ans1>maxd)	ans.clear(),maxd=ans2-ans1,ans.push_back(ans1);
	}
	vector<int>::iterator it=unique(ans.begin(),ans.end());
	cout<<it-ans.begin()<<" "<<maxd<<'\n';
	for(re vector<int>::iterator i=ans.begin();i!=it;++i)	cout<<*i<<' ';
}
posted @ 2018-11-02 17:00  Patrickpwq  阅读(289)  评论(0编辑  收藏  举报