码迷,mamicode.com
首页 > 其他好文 > 详细

[ACM] POJ 2409 Let it Bead (Polya计数)

时间:2018-09-04 13:50:40      阅读:147      评论:0      收藏:0      [点我收藏+]

标签:div   参考   poj   pac   tar   name   while   amp   stream   

参考:https://blog.csdn.net/sr_19930829/article/details/38108871

 1 #include <iostream>
 2 #include <string.h>
 3 using namespace std;
 4 int c,n;//c种颜色,n个珠子
 5 int ans;
 6 
 7 int gcd(int a,int b)
 8 {
 9     return b==0?a:gcd(b,a%b);
10 }
11 
12 int power(int p,int n)//快速幂求公式里的k的nc(g)次方
13 {
14     int ans=1;
15     while(n)
16     {
17         if(n&1)
18             ans*=p;
19         p*=p;
20         n/=2;
21     }
22     return ans;
23 }
24 
25 int main()
26 {
27     while(cin>>c>>n&&(c||n))
28     {
29         ans=0;
30         for(int i=1;i<=n;i++)
31             ans+=power(c,gcd(n,i));//相当于求图3-1中左边的循环群染色数
32         if(n&1)//是奇数,有n个包含(n/2+1)个循环节的循环群
33             ans+=n*power(c,n/2+1);
34         else
35             ans+=(power(c,n/2+1)+power(c,n/2))*(n/2);
36         ans/=2*n;//别忘了除以置换群的总个数
37         cout<<ans<<endl;
38     }
39     return 0;
40 }

 

技术分享图片

技术分享图片

技术分享图片

技术分享图片

 

[ACM] POJ 2409 Let it Bead (Polya计数)

标签:div   参考   poj   pac   tar   name   while   amp   stream   

原文地址:https://www.cnblogs.com/hemeiwolong/p/9583595.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!