bzoj 2081: [Poi2010]Beads

2081: [Poi2010]Beads

Time Limit: 10 Sec   Memory Limit: 259 MB
Submit: 346   Solved: 135
[ Submit][ Status][ Discuss]

Description

Zxl有一次决定制造一条项链，她以非常便宜的价格买了一长条鲜艳的珊瑚珠子，她现在也有一个机器，能把这条珠子切成很多块（子串），每块有k（k>0）个珠子，如果这条珠子的长度不是k的倍数，最后一块小于k的就不要拉（nc真浪费），保证珠子的长度为正整数。 Zxl喜欢多样的项链，为她应该怎样选择数字k来尽可能得到更多的不同的子串感到好奇，子串都是可以反转的，换句话说，子串（1，2，3）和（3,2,1）是一样的。写一个程序，为Zxl决定最适合的k从而获得最多不同的子串。 例如：这一串珠子是： (1,1,1,2,2,2,3,3,3,1,2,3,3,1,2,2,1,3,3,2,1)， k=1的时候，我们得到3个不同的子串： (1),(2),(3) k=2的时候，我们得到6个不同的子串： (1,1),(1,2),(2,2),(3,3),(3,1),(2,3) k=3的时候，我们得到5个不同的子串： (1,1,1),(2,2,2),(3,3,3),(1,2,3),(3,1,2) k=4的时候，我们得到5个不同的子串： (1,1,1,2),(2,2,3,3),(3,1,2,3),(3,1,2,2),(1,3,3,2)

Input

共有两行，第一行一个整数n代表珠子的长度，（），第二行是由空格分开的颜色ai（1<=ai<=n）。

Output

也有两行，第一行两个整数，第一个整数代表能获得的最大不同的子串个数，第二个整数代表能获得最大值的k的个数，第二行输出所有的k（中间有空格）。

Sample Input

21
1 1 1 2 2 2 3 3 3 1 2 3 3 1 2 2 1 3 3 2 1

Sample Output

6 1
2

HINT

Source

by poi

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<map>
#define N 1000003
#define p 200191
#define ll unsigned long long 
using namespace std;
int n,m,a[N],maxn;
ll s1[N],s2[N],mul[N];
vector<int> num;
map<ll,int> mp;
ll gethash(int l,int r)
{
  if (l<=r)  return s1[r]-s1[l-1]*mul[r-l+1];
  else return s2[r]-s2[l+1]*mul[l-r+1];
}
void solve (int x)
{
  if (maxn*x>n) return;
  int ans=0;
  mp.clear();
  for (int i=1;i<=n;i+=x)
  if (i+x-1<=n)
  {
  	ll t=gethash(i,i+x-1)*gethash(i+x-1,i);
  	if (!mp[t]) ans++,mp[t]=1; 
  }
  if (ans>maxn) maxn=ans,num.clear();
  if (ans==maxn) num.push_back(x);
}
int main()
{
 scanf("%d",&n);
 for (int i=1;i<=n;i++) scanf("%d",&a[i]);
 mul[0]=1;
 for (int i=1;i<=n;i++) mul[i]=mul[i-1]*p;
 for (int i=1;i<=n;i++)  s1[i]=s1[i-1]*p+a[i];
 for (int i=n;i>=1;i--)	 s2[i]=s2[i+1]*p+a[i];
 for (int i=1;i<=n;i++) solve(i);
 printf("%d %d\n",maxn,num.size());
 for(int i=0;i<num.size();i++)
 if (i!=num.size()-1) printf("%d ",num[i]);
 else printf("%d\n",num[i]);
}