• BZOJ2795: [Poi2012]A Horrible Poem


    题解:

    首先判断长度len是否是l-r的循环节可以做到O(1)

    只需判断 l-l+len-1 和 r-len+1-r 是否相同即可,这一步可以用hash

    然后循环节一定是总长度的约数。

    这样可以做到n√n,但是会T。

    一个鬼畜的优化是 枚举在每块内每个字母的个数,这个个数必须是总个数的约数,然后我们取一下gcd,再枚举约数。

    实测可过。

    另一种做法可以戳这里:http://blog.csdn.net/zeyu_king/article/details/41989279

    代码:

     1 #include<cstdio>
     2 #include<cstdlib>
     3 #include<cmath>
     4 #include<cstring>
     5 #include<algorithm>
     6 #include<iostream>
     7 #include<vector>
     8 #include<map>
     9 #include<set>
    10 #include<queue>
    11 #include<string>
    12 #define inf 1000000000
    13 #define maxn 600000
    14 #define maxm 500+100
    15 #define eps 1e-10
    16 #define ull unsigned long long
    17 #define pa pair<int,int>
    18 #define for0(i,n) for(int i=0;i<=(n);i++)
    19 #define for1(i,n) for(int i=1;i<=(n);i++)
    20 #define for2(i,x,y) for(int i=(x);i<=(y);i++)
    21 #define for3(i,x,y) for(int i=(x);i>=(y);i--)
    22 #define mod 1000000007
    23 using namespace std;
    24 inline int read()
    25 {
    26     int x=0,f=1;char ch=getchar();
    27     while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    28     while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}
    29     return x*f;
    30 }
    31 int n,m,a[maxn],c[maxn],s[maxn][30];
    32 ull b[maxn],power[maxn];
    33 const ull base=13131;
    34 inline ull get(int x,int y){return b[x]-b[x+y]*power[y];}
    35 inline int gcd(int x,int y){return y?gcd(y,x%y):x;}
    36 inline int solve()
    37 {
    38     int x=read(),y=read(),z=y-x+1,len=y-x+1,t=0;
    39     for1(i,26)z=gcd(z,s[y][i]-s[x-1][i]);
    40     if(z==1)return len;
    41     for(int i=1;i*i<=z;i++)if(z%i==0)
    42     {
    43        int tmp=len/z*i;
    44        if(get(x,len-tmp)==get(y-(len-tmp)+1,len-tmp))return tmp;
    45        c[++t]=len/i;
    46     }
    47     for3(i,t,1)if(get(x,len-c[i])==get(y-(len-c[i])+1,len-c[i]))return c[i];
    48     return len;
    49 }
    50 int main()
    51 {
    52     freopen("input.txt","r",stdin);
    53     freopen("output.txt","w",stdout);
    54     n=read();n=0;
    55     char ch=getchar();
    56     while(ch<'a'||ch>'z')ch=getchar();
    57     for(;ch>='a'&&ch<='z';ch=getchar())a[++n]=ch-'a'+1;
    58     for3(i,n,1)b[i]=b[i+1]*base+a[i];
    59     for1(i,n)for1(j,26)s[i][j]=s[i-1][j]+(a[i]==j);
    60     power[0]=1;
    61     for1(i,n)power[i]=power[i-1]*base;
    62     m=read();
    63     while(m--)printf("%d
    ",solve());
    64     return 0;
    65 }
    View Code

    2795: [Poi2012]A Horrible Poem

    Time Limit: 50 Sec  Memory Limit: 128 MB
    Submit: 121  Solved: 72
    [Submit][Status]

    Description


    给出一个由小写英文字母组成的字符串S,再给出q个询问,要求回答S某个子串的最短循环节。
    如果字符串B是字符串A的循环节,那么A可以由B重复若干次得到。

    Input

    第一行一个正整数n (n<=500,000),表示S的长度。
    第二行n个小写英文字母,表示字符串S。
    第三行一个正整数q (q<=2,000,000),表示询问个数。
    下面q行每行两个正整数a,b (1<=a<=b<=n),表示询问字符串S[a..b]的最短循环节长度。

    Output

    依次输出q行正整数,第i行的正整数对应第i个询问的答案。

    Sample Input

    8
    aaabcabc
    3
    1 3
    3 8
    4 8

    Sample Output

    1
    3
    5

    HINT

     

    Source

  • 相关阅读:
    Codeforces Round #581 (Div. 2)
    Codeforces Round #605 (Div. 3)
    cin,cin.get(),cin.getline(),getline()
    容斥原理原理
    词法分析器
    Web学习开始。
    Convex hull凸包
    对max_flow做一个总结
    Vue中动画封装
    Vue中的动画特效
  • 原文地址:https://www.cnblogs.com/zyfzyf/p/4198545.html
Copyright © 2020-2023  润新知