Magic Bracelet POJ - 2888

题意：

不同项链个数，长度n，（n<=1e9）,颜色m种，m<=10 ,给定k个数对（a，b）表示a颜色和b颜色不能相邻

思路：

n过多的时候就自然的想到    有 euler（n/i）个   间隔为j的置换群，可以直接得到 循环节个数gcd（j，n）等于i的个数

统计n的约数就可以在O（sqrt（n））的复杂度完成，

我们知道  gcd 为循环节个数 ，也即是一个旋转的原串（个人的理解）

如图  gcd=5，只要这个串是合法的就可以是一种方案，

合法的方案是   第一个与第二个相邻合法，  第二个与第三个相邻合法，。。。。

直到最后一个和下一串的第一个合法，由于该方案是该置换群下的不动点

所以下一串的第一个也就是这一串的第一个

不经惭愧啊，我也是看了大佬们的博客才想起邻接矩阵的

a[i][j]=1表示从i到j是合法的相邻，a[[i][j]=0表示从i到j是不合法的相邻，

矩阵开n幂后，a[i][j]表示i经过n条路到j的路的条数

而答案需要的也就是 一种颜色出发经过gcd-1种颜色后回到自己的方案数

也就是邻接矩阵开n方后 的 sigma  a[ i ][ i ]

#include<iostream>
#include<string.h>
#include<stdio.h>
#define ll int
#define mod  9973
using namespace std;
struct node
{
	ll a[10][10];
};

node A;
node B;
int n,m;

node cheng(node a,node b)
{
	node c;
	memset(c.a,0,sizeof(c.a));
	for(int i=0;i<m;i++)
	for(int j=0;j<m;j++)
	for(int k=0;k<m;k++)
	if(a.a[i][k]&&b.a[k][j])
	c.a[i][j]=(c.a[i][j]+a.a[i][k]*b.a[k][j])%mod;
	return  c;
}

ll s(ll n)
{
	node x=A;
	node y=B;
	while(n)
	{
		if(n&1) y=cheng(y,x);
		x=cheng(x,x);
		n>>=1;
	}
   ll num=0;
   for(int i=0;i<m;i++)
   	 num=(num+y.a[i][i])%mod;
   
  return num;
	 
}

ll qpow(ll a,ll b)
{
	ll ans=1;
	a%=mod;
	while(b)
	{
		if(b&1) ans=(ans*a)%mod;;
		a=(a*a)%mod;
		b>>=1;
	}
	return ans;
}

ll euler(ll n)
{
	ll res=n;
	for(int i=2;i*i<=n;i++)
	{
		if(n%i==0)
		{
			res/=i;
			res*=(i-1);
			while(n%i==0) n/=i;
		}
	}
	if(n>1) 
	{
		res/=n;
			res*=(n-1);
	}
	return res%mod;
}


ll solve()
{

   ll ans=0;
   
   for(ll i=1;i*i<=n;i++)
   {
   		if(n%i==0)
   		{		
   		   ans=(ans+euler(n/i)*s(i))%mod;  //euler(n/i) gcd=i 的个数     gcd为循环节个数 
		   if(i*i!=n)
		   	ans=(ans+euler(i)*s(n/i))%mod;
		   	
   		}
   }

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

int main()
{

	int T;
	int k,a,b;
	cin>>T;
	memset(B.a,0,sizeof(B.a));
	for(int i=0;i<10;i++)
	B.a[i][i]=1;
	while(T--)
	{
   scanf("%d%d%d",&n,&m,&k);
	for(int i=0;i<m;i++)
	for(int j=0;j<m;j++)
	A.a[i][j]=1;
	
	for(int i=0;i<k;i++)
	{
	scanf("%d%d",&a,&b);	
	A.a[a-1][b-1]=A.a[b-1][a-1]=0;	
	}
	
	printf("%d\n",solve());	
		
	}


	return 0;
}