用m种颜色给n个点的一个环染色,有一些限制条件,表示颜色a和颜色b不能相邻,考虑旋转同构,求方案数。(m<=10,n<=1e9)
n很大,只能
O(n−√)
枚举循环节个数k,枚举素数计算欧拉函数。因为颜色的限制,我们现在的问题在于要计算:给定一些颜色不能相邻,把k个位置染色的方案数。(这里准确的说,不仅仅是给k个位置染色,还要保证头尾颜色也要合法),因此其实就是求从颜色i走K步回到颜色i的方案数。我们可以构建矩阵来计算。复杂度
O(m3logn)
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define inf 0x3f3f3f3f
#define ll long long
#define mod 9973
#define N 100010
inline char gc(){
static char buf[1<<16],*S,*T;
if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
return *S++;
}
inline int read(){
int x=0,f=1;char ch=gc();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
return x*f;
}
int n,m,prime[N],tot=0,phi[N];
bool notprime[N],mark[12];
struct Matrix{
int a[11][11];
Matrix(){}
Matrix(bool t){memset(a,0,sizeof(a));if(t) for(int i=1;i<=m;++i) a[i][i]=1;}
friend Matrix operator*(Matrix a,Matrix b){
Matrix res(0);
for(int i=1;i<=m;++i)
for(int j=1;j<=m;++j)
for(int k=1;k<=m;++k) res.a[i][j]=(res.a[i][j]+a.a[i][k]*b.a[k][j])%mod;
return res;
}friend Matrix operator^(Matrix a,int k){
Matrix res(1);for(;k;k>>=1,a=a*a) if(k&1) res=res*a;return res;
}
}a;
inline void init(){
notprime[1]=1;phi[1]=1;
for(int i=2;i<=100000;++i){
if(!notprime[i]) prime[++tot]=i,phi[i]=i-1;
for(int j=1;prime[j]*i<=100000;++j){
notprime[prime[j]*i]=1;
if(i%prime[j]==0){phi[i*prime[j]]=phi[i]*prime[j];break;}
phi[i*prime[j]]=phi[i]*phi[prime[j]];
}
}
}
inline int ph(int x){
if(x<=100000) return phi[x]%mod;
int res=x,xx=x;
for(int j=1;prime[j]*prime[j]<=x&&xx!=1;++j)
if(xx%prime[j]==0){
res-=res/prime[j];
while(xx%prime[j]==0) xx/=prime[j];
}if(xx!=1) res-=res/xx;return res%mod;
}
inline int calc(int x){
Matrix ans(0);
ans=a^x;int res=0;
for(int i=1;i<=m;++i) res+=ans.a[i][i];return res%mod;
}
inline int ksm(int x,int k){
x%=mod;int res=1;for(;k;k>>=1,x=x*x%mod) if(k&1) res=res*x%mod;return res;
}
int main(){
// freopen("a.in","r",stdin);
int tst=read();init();
while(tst--){
n=read();m=read();int owo=read(),ans=0;memset(mark,0,sizeof(mark));
for(int i=1;i<=m;++i)
for(int j=1;j<=m;++j) a.a[i][j]=1;
while(owo--){
int x=read(),y=read();a.a[x][y]=a.a[y][x]=0;
}int x=1;
for(;x*x<n;++x) if(n%x==0) ans+=calc(x)*ph(n/x)+calc(n/x)*ph(x),ans%=mod;
if(x*x==n) ans+=calc(x)*ph(x),ans%=mod;
printf("%d\n",ans*ksm(n,mod-2)%mod);
}return 0;
}