BZOJ-2081-[Poi2010]Beads（hash+暴力）

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代表珠子的长度，（n<=200000），第二行是由空格分开的颜色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

 

题解

枚举长度，哈希暴力。

但是——————

被这道题坑了好久，模数怎么取交上去都WA，本来自己是正着和倒着都判一遍，后来参考了hzw学长的代码发现可以乘起来，这样正确率会高很多

这里有一个小优化：如果当前最大值乘上枚举的长度大于n的话，就可以直接break了

P.s.第一道权限题，用某人的号交的0.0

#include<bits/stdc++.h>
#define N 200005
#define mod 998532737 
#define ll long long
using namespace std;
const ll x=12138;
int n,Max,tot,num;
int a[N],t[N];
ll b[N],h1[N],h2[N];
map<ll,bool> mp;
int main(){
    scanf("%d",&n);
    for (int i=1;i<=n;i++) scanf("%d",&a[i]);
    b[0]=1;
    for (int i=1;i<=n;i++) b[i]=(b[i-1]*x)%mod;
    for (int i=1;i<=n;i++) h1[i]=(h1[i-1]*x+a[i])%mod;
    for (int i=n;i>=1;i--) h2[i]=(h2[i+1]*x+a[i])%mod;
    for (int k=1;k<=n;k++){
        mp.clear();
        if (n<Max*k) break;
        num=0;
        for (int i=1;i<=n/k;i++){
            ll s1=((h1[i*k]-h1[(i-1)*k]*b[k])%mod+mod)%mod,s2=((h2[(i-1)*k+1]-h2[i*k+1]*b[k])%mod+mod)%mod;
            if (!mp[s1*s2]){
                num++;
                mp[s1*s2]=true;
            }
        }
        if (num>Max) Max=num,tot=1,t[1]=k; else
        if (num==Max) t[++tot]=k;
    }
    printf("%d %d\n",Max,tot);
    for (int i=1;i<=tot-1;i++)
        printf("%d ",t[i]);
    printf("%d\n",t[tot]);
    return 0;
}