题意:等价于在一个环上,有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=∑i=1nf(gcd(i,n))
把这个公式优化一下就好了
ans∑d|nf(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;
}