HDU-5868-Different Circle Permutation-Burnside定理和polay计数

题目链接

题意：等价于在一个环上，有n个点，可以染成白色和黑色，问在黑色不可以相邻，在旋转同构意义下的方案数有多少？

环上的计数问题，基本都可以用burnside定理和polay计数来计算。

思路：首先先确定有n种置换，即旋转n个角度。对于每一种置换i,其循环节长度为$gcd(i,n)$；对于每个长度循环节要求$calc(i)$其实就是$f(gcd(i,n))$ 其中$f(x)$表示没有旋转同构下，环的个数。而$f(n)$的递推式子很好求。

$f(n)=f(n-1)+f(n-2)$用个矩阵快速幂就很快能求出来了。

接着公式就出来了 $$ans=\sum^{n}_{i=1}{f(gcd(i,n))}$$

把这个公式优化一下就好了 $$ans\sum^{}_{d|n}{f(d)*euler(n/d)}$$

ps注意特判n=1

#include<bits/stdc++.h>
using namespace std;
const long long mod=1000000007;
long long extend_gcd(long long a,long long b,long long &x,long long &y)
{
    if(a==0&&b==0) return -1;
    if(b==0){x=1;y=0;return a;}
    long long d=extend_gcd(b,a%b,y,x);
    y-=a/b*x;
    return d;
}
long long mod_rev(long long a,long long n)
{
    long long x,y;
    long long d=extend_gcd(a,n,x,y);
    if(d==1) return (x%n+n)%n;
    else return -1;
}
int euler(int n) {
    int ans=n;
    for(int i=2;i*i<=n;i++) {
        if(n%i==0) {
            ans=ans/i*(i-1);
            while(n%i==0) n/=i;
        }
    }
    if(n>1) ans=ans/n*(n-1);
    return ans;
}
struct MAT{
    long long a[2][2];
};
MAT mul(MAT x,MAT y) {
    MAT tt;
    for(int i=0;i<=1;i++)
        for(int j=0;j<=1;j++) {
            tt.a[i][j]=0;
            for(int k=0;k<=1;k++) {
                tt.a[i][j]+=x.a[i][k]*y.a[k][j]%mod;
                tt.a[i][j]%=mod;
            }
        }
    return tt;
}
long long fast_mat(int x) {
    MAT ans,tt;
    ans.a[0][0]=ans.a[1][1]=1;
    ans.a[0][1]=ans.a[1][0]=0;
    tt.a[0][0]=tt.a[0][1]=tt.a[1][0]=1;
    tt.a[1][1]=0;
    x-=3;
    while(x) {
        if(x&1) {
            ans=mul(ans,tt);
        }
        tt=mul(tt,tt);
        x/=2;
    }
    return (ans.a[0][0]*4+ans.a[0][1]*3)%mod;
}
long long calc(int x) {
    if(x==1) return 1;
    if(x==2) return 3;
    if(x==3) return 4;
    return fast_mat(x);
}
long long sov(int n) {
    if(n==1) return 2;
    long long ans=0;
    for(int i=1;i*i<=n;i++) {
        if(n%i==0) {
            if(i*i==n) {
                ans=(ans+calc(i)*euler(i)%mod)%mod;
            }
            else {
                ans=(ans+calc(i)*euler(n/i)%mod+calc(n/i)*euler(i)%mod)%mod;
            }
        }
    }
    ans=ans*mod_rev(n,mod)%mod;

    return ans;
}
int main()
{
    int n;
    while(scanf("%d",&n)!=EOF) {
        long long ans=sov(n);
        printf("%lld\n",ans);
    }
    return 0;
}