【noip模拟赛】数字对

http://hzwer.com/5015.html

好题一道，二分的单调性很隐秘，还介绍了ST算法的拓展

【题目描述】

小H是个善于思考的学生，现在她又在思考一个有关序列的问题。

她的面前浮现出一个长度为n的序列{ai}，她想找出一段区间[L, R](1 <= L <= R <= n)。

这个特殊区间满足，存在一个k(L <= k <= R)，并且对于任意的i(L <= i <= R)，ai都能被ak整除。这样的一个特殊区间 [L, R]价值为R - L。

小H想知道序列中所有特殊区间的最大价值是多少，而有多少个这样的区间呢？这些区间又分别是哪些呢？你能帮助她吧。

【输入格式】

第一行，一个整数n.

第二行，n个整数，代表ai.

【输出格式】

第一行两个整数，num和val,表示价值最大的特殊区间的个数以及最大价值。

第二行num个整数，按升序输出每个价值最大的特殊区间的L.

【样例输入1】

5

4 6 9 3 6

【样例输出1】

1 3

2

【样例输入2】

5

2 3 5 7 11

【样例输出2】

5 0

1 2 3 4 5

【数据范围】

30%: 1 <= n <= 30 , 1 <= ai <= 32.

60%: 1 <= n <= 3000 , 1 <= ai <= 1024.

80%: 1 <= n <= 300000 , 1 <= ai <= 1048576.

100%: 1 <= n <= 500000 , 1 <= ai < 2 ^ 31.

题解

30% ：暴力枚举判断。O(n^4) or O(n^3)。

60% ：特殊区间的特点实际上就是区间最小值等于这个区间的GCD，于是暴力或递推算出每个区间最小值与GCD。而对于最大价值，可以通过二分来进行求解。复杂度O(n ^ 2)。

100%：在60%的基础上，最小值与GCD都使用RMQ算法来求解，对于这道题推荐使用ST表，线段树会被艹飞。复杂度O(nlogn)。

实际上此题数据水机智的暴力踩标程

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#define ll long long
#define inf 2147483647
using namespace std;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int bin[20];
int n,top,res;
int a[500005],L[500005],q[500005],ans[500005];
int mn[500005][19],g[500005][19];
int gcd(int x,int y)
{
	return y==0?x:gcd(y,x%y);
}
void pre()
{
	for(int i=1;i<=n;i++)
		mn[i][0]=g[i][0]=a[i];
	for(int j=1;j<=18;j++)
		for(int i=1;i<=n;i++)
		{
			if(i+bin[j]-1<=n)
			{
				mn[i][j]=min(mn[i][j-1],mn[i+bin[j-1]][j-1]);
				g[i][j]=gcd(g[i][j-1],g[i+bin[j-1]][j-1]);
			}
			else break;
		}
}
bool query(int l,int r)
{
	int t=L[r-l+1],R=r-bin[t]+1;
	return min(mn[l][t],mn[R][t])==gcd(g[l][t],g[R][t]);
}
void solve(int x)
{
	for(int i=1;i<=n;i++)
		if(i+x-1<=n)
			if(query(i,i+x-1))
				ans[++top]=i;
}
int main()
{
	//freopen("pair.in","r",stdin);
	//freopen("pair.out","w",stdout);
	bin[0]=1;for(int i=1;i<20;i++)bin[i]=bin[i-1]<<1;
	L[0]=-1;for(int i=1;i<=500000;i++)L[i]=L[i>>1]+1;
	n=read();
	for(int i=1;i<=n;i++)a[i]=read();
	pre();
	int l=1,r=500000;
	while(l<=r)
	{
		int mid=(l+r)>>1;
		top=0;
		solve(mid);
		if(top)res=mid,l=mid+1;
		else r=mid-1;
	}
	top=0;
	solve(res);
	printf("%d %d\n",top,res-1);
	for(int i=1;i<=top;i++)
		printf("%d ",ans[i]);
	return 0;
}