【bzoj2081】[Poi2010]Beads

4 篇文章 0 订阅

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

题解
枚举长度Hash判重,为什么我找了好久的模数才能过,数据也太强了吧

#include<bits/stdc++.h>
#define N 500005
#define ll long long
#define inf 1000000009
#define mod 2010516623LL
using namespace std;
inline 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 n;
map<ll,bool>mp;
ll fac[200005],a[200005],sum1[200005],sum2[200005],ANS=0;
vector<int>q;
ll get1(int l,int r)
{
    ll ans=(sum1[r]-sum1[l-1]+mod)%mod;
    return ans*fac[n-r]%mod;
}
ll get2(int l,int r)
{
    ll ans=(sum2[r]-sum2[l-1]+mod)%mod;
    return ans*fac[n-r]%mod;
}
int main()
{
    n=read();for (int i=1;i<=n;i++) a[i]=read();
    fac[0]=1LL;
    for (int i=1;i<=n;i++)fac[i]=(ll)fac[i-1]*666623333LL%mod;
    for (int i=1;i<=n;i++)
    {
        sum1[i]=(sum1[i-1]+(ll)fac[i]*a[i]%mod)%mod;
        sum2[i]=(sum2[i-1]+(ll)fac[i]*a[n-i+1]%mod)%mod;
    }
    for (int i=1;i<=n;i++)
    {
        mp.clear();
        int ans=0;
        if((ll)ANS*i>n) break;
        for (int j=1;j+i-1<=n;j+=i)
        {
            int Hash1=get1(j,j+i-1);
            int Hash2=get2(n-(j+i-1)+1,n-j+1);
            int Hash=(ll)Hash1*Hash2%mod;
            if (!mp[Hash]){mp[Hash]=1;ans++;}
        }
        if (ans>ANS){ANS=ans;q.clear();q.push_back(i);}
        else if (ans==ANS) q.push_back(i);
    }
    printf("%lld %d\n",ANS,q.size());
    for (int i=0;i<q.size();i++)
    {
        printf("%d",q[i]);if (i<q.size()-1) printf(" ");
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值