jzoj3895. 数字对（B组——Day7）

jzoj3895. 数字对（B组——Day7）

题目

Description

小H是个善于思考的学生，现在她又在思考一个有关序列的问题。
她的面前浮现出一个长度为n的序列{ai}，她想找出一段区间[L, R](1 <= L <= R <= n)。
这个特殊区间满足，存在一个k(L <= k <= R)，并且对于任意的i(L <= i <= R)，ai都能被ak整除。这样的一个特殊区间 [L, R]价值为R - L。
小H想知道序列中所有特殊区间的最大价值是多少，而有多少个这样的区间呢？这些区间又分别是哪些呢？你能帮助她吧。

Input

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

Output

第一行两个整数，num和val,表示价值最大的特殊区间的个数以及最大价值。
第二行num个整数，按升序输出每个价值最大的特殊区间的L.

Sample Input

输入1：

5
4 6 9 3 6

输入2：

5
2 3 5 7 11

Sample Output

输出1：

1 3
2

输出2：

5 0
1 2 3 4 5

Data Constraint

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.

解析

jzoj的题解说明

30% ：

暴力枚举判断。O(n^4)。

60% ：

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

100%：

在 60%的基础上，最小值与 GCD 都使用 RMQ 算法来求解，对于这道题推荐使用ST 表。最大价值仍然使用二分。复杂度 O(nlogn)。

自己的思路

想不到正解，直接枚举1~n为a[k]，再以a[k]为中心，依次向左和右寻找，直到不能整除a[k]为止。每次寻找得出的结果R-L若大于之前的则将记录L的数组清零，重新开始记录；若等于则直接记录。最后输出前还要排序…
比赛的时候没排，结果暴蛋QWQ
注意：当L相同时，不用重复记录。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;

int n,maxn;
int f[500010];
long long a[500010];

int main(){
	scanf("%d",&n);
	for (int i=1;i<=n;i++)
		scanf("%lld",&a[i]);
	for (int i=1;i<=n;i++){
		int l=i,r=i;
		while (a[l-1]%a[i]==0 && l-1>0) l--;
		while (a[r+1]%a[i]==0 && r<n) r++;
		if (maxn<r-l){
			memset(f,0,sizeof(f));
			f[0]++;
			f[f[0]]=l;
			maxn=r-l;
		}
		else if (maxn==r-l && f[f[0]]!=l){
			f[0]++;
			f[f[0]]=l;
		}
	}
	sort(f+1,f+f[0]+1);
	printf("%d %d\n",f[0],maxn);
	for (int i=1;i<=f[0];i++)
		printf("%d ",f[i]);
	return 0;
}