组合数学中的项链计数

给c种不同颜色宝石能穿成多少种长度为s的宝石项链（本质不同）

Burnside定理的应用：

当n为奇数时，有n种翻转，每种翻转都是以一个顶点和该顶点对边的中点对称。有k^(n/2+1)*n种。

当n为偶数时，有n种翻转，其中一半是以两个对应顶点，另一半是以两条对边对称。有k^(n/2+1)*n/2+k^(n/2)*n/2种。

考虑旋转：枚举旋转角度360/n*i，(0<i<=n)，也就是一个置换。经过该置换，颜色仍保持不变的着色方案有k^GCD(n,i)种。

 

一个长度为n的环，每i个上同一种颜色，可以上多少种颜色。

假设起点在x，则x，x+i，x+2*i，……，x+k*i，……

假设在第t次，第一次回到起点，则x=(x+t*i)%n => t*i%n=0 => t=LCM(i,n)/i=n*i/GCD(n,i)/i=n/GCD(n,i)。

那么可以上n/t种颜色，即n/(n/GCD(n,i))种，所以旋转的着色方案有k^GCD(n,i)种。

 
 

   
   

    
    

   
   
   
   

    
    

     
     #include<cstdio>
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
     #include<iostream>
    
    
   
   

   
   

    
    

   
   
   
   

    
    
 
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
     using 
     
     namespace 
     
     std;
    
    
   
   

   
   

    
    

   
   
   
   

    
    
 
    
    
   
   

   
   

    
    

   
   
   
   

    
    
 
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
     typedef 
     
     long 
     
     long ll;
    
    
   
   

   
   

    
    

   
   
   
   

    
    
 
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
     ll power(ll a,ll b)
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
     {
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
         ll ans=
     
     1ll;
    
    
   
   

   
   

    
    

   
   
   
   

    
    
    
     
     while(b)
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
         {
    
    
   
   

   
   

    
    

   
   
   
   

    
    
        
     
     if(b&
     
     1)
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
             ans=ans*a;
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
             a=a*a;
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
             b=b>>
     
     1;
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
         }
    
    
   
   

   
   

    
    

   
   
   
   

    
    
    
     
     return ans;
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
     }
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
     ll gcd(ll a,ll b)
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
     {
    
    
   
   

   
   

    
    

   
   
   
   

    
    
    
     
     return b ? gcd(b,a%b)  :  a;
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
     }
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
     int main()
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
     {
    
    
   
   

   
   

    
    

   
   
   
   

    
    
    
     
     int c,s;
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
         ll ans;
    
    
   
   

   
   

    
    

   
   
   
   

    
    
    
     
     while(
     
     scanf(
     
     "%d%d",&c,&s)!=EOF)
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
         {
    
    
   
   

   
   

    
    

   
   
   
   

    
    
        
     
     if(s&
     
     1)
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
             ans=power(c,s/
     
     2+
     
     1)*s;
    
    
   
   

   
   

    
    

   
   
   
   

    
    
        
     
     else
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
             ans=power(c,s/
     
     2)*(s/
     
     2)+power(c,s/
     
     2+
     
     1)*(s/
     
     2);
    
    
   
   

   
   

    
    

   
   
   
   

    
    
        
     
     for(
     
     int i=
     
     1;i<=s;i++)
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
             ans+=power(c,gcd(s,i));
    
    
   
   

   
   

    
    

   
   
   
   

    
    
        
     
     printf(
     
     "%lld\n",(ans/
     
     2)/s);
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
         }
    
    
   
   

   
   

    
    

   
   
   
   

    
    
    
     
     return 
     
     0;
    
    
   
   

   
   

    
    

   
   
   
   

    
    

     
     }