hdu 4187 Alphabet Soup

这题的主要就是找循环节数，这里用找字符串最小覆盖来实现，也就是n-next[n],证明在这http://blog.csdn.net/fjsd155/article/details/6866991

#include<iostream>

#include<stdio.h>

#include<algorithm>

#include<iomanip>

#include<cmath>

#include<string>

#include<vector>

using namespace std;

__int64 mod=100000007;

int dis[360005],data[360005],next[360005];

__int64 euler(__int64 n)

{

    int i;

    __int64 ans=1;

    for(i=2;i*i<=n;i++)

        if(n%i==0)

        {

            ans*=i-1;

            n/=i;

            while(n%i==0)

            {

                ans*=i;

                n/=i;

            }

        }

    if(n>1) ans*=n-1;

    return ans%mod;

}

int get_next(int n)

{

    int i=0,j=-1;

    next[0]=-1;

    while(i<n)

    {

        if(j==-1||dis[i]==dis[j])

        {

            i++;

            j++;

            next[i]=j;

        }

        else j=next[j];

    }

    i=n-j;

    if(n%i)

        return n;

    return i;

}

__int64 pows(__int64 a,__int64 b)

{

    __int64 ans=1;

    a%=mod;

    while(b)

    {

        if(b&1) ans=(ans*a)%mod;

        b>>=1;

        a=(a*a)%mod;

    }

    return ans%mod;

}

__int64 polya(__int64 m,__int64 n)

{

    int i,j;

    __int64 ans=0;

    for(i=1;i*i<=n;i++)

    {

        if(n%i==0)

        {

            ans=(ans+pows(m,n/i)*euler(i)%mod)%mod;

            if(i*i!=n)

                ans=(ans+pows(m,i)*euler(n/i)%mod)%mod;

        }

    }

    return (ans*pows(n,mod-2))%mod;

}

int main()

{

    int n,i,j,k,t,m,s,p;

    while(cin>>s>>p)

    {

        if(s==-1&&p==-1) break;

        for(i=0;i<p;i++) cin>>data[i];

        sort(data,data+p);

        for(i=1;i<p;i++)

            dis[i]=data[i]-data[i-1];

        dis[0]=360000-data[p-1]+data[0];

        int len=get_next(p);

        printf("%I64d\n",polya(pows(s,len),p/len));

    }

    return 0;

}

View Code