标签：include int ybtoj ans 珠子 while 项链 例题
题目
给定 m m m 种颜色的珠子，每种颜色珠子的个数均不限，将这些珠子做成长度为 n n n 的项链。
问能做成多少种不重复的项链。
两条项链相同，当且仅当两条项链通过旋转或是沿中轴线翻转后能重合在一起，且对应珠子的颜色相同。
多组数据
n<32,m<32
思路
polya计数定理。
先考虑旋转。
旋转共n种置换，第i种的轮换个数为gcd(n,i)
接下来是对称。
对称分奇偶性讨论
奇数部分有
n
n
n种置换，轮换个数均为
(
n
+
1
)
/
2
(n+1)/2
(n+1)/2。
偶数部分又分2部分
第一部分以中线对称，共
n
/
2
n/2
n/2种，每种轮换均为
n
/
2
+
1
n/2+1
n/2+1。
第二部分以对点连线对称，共
n
/
2
n/2
n/2种，轮换数量均为
n
/
2
n/2
n/2。
所以旋转的公式为
∑
i
=
1
n
m
g
c
d
(
i
,
n
)
\sum_{i=1}^n m^{gcd(i,n)}
∑i=1nmgcd(i,n)
奇数公式为
n
∗
m
n
/
2
+
n
m
o
d
2
n*m^{n/2+n\bmod2}
n∗mn/2+nmod2
偶数公式为
n
/
2
∗
(
m
n
/
2
+
1
+
m
n
/
2
)
n/2*(m^{n/2+1}+m^{n/2})
n/2∗(mn/2+1+mn/2)
记得最后除
2
∗
n
2*n
2∗n（
∣
G
∣
|G|
∣G∣）
code：
#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
//When I wrote this code,God and I unterstood what was I doing
inline int read()
{
int ret,c,f=1;
while (((c=getchar())> '9'||c< '0')&&c!='-');
if (c=='-') f=-1,ret=0;
else ret=c-'0';
while ((c=getchar())>='0'&&c<='9') ret=(ret<<1)+(ret<<3)+(c^48);
return ret*f;
}
int n,m;
int ksm(int x,int y)
{
int ans=1;
while (y)
{
if (y&1) ans=ans*x;
x=x*x;
y>>=1;
}
return ans;
}
int gcd(int x,int y)
{
int r=x%y;
while (r)
{
x=y,y=r,r=x%y;
}
return y;
}
int main()
{
n=read(),m=read();
while (n!=0||m!=0)
{
int ans=0;
for (int i=1;i<=m;i++) ans+=ksm(n,gcd(i,m));
if (m&1) ans+=m*ksm(n,m/2+m%2);
else ans+=m/2*(ksm(n,m/2+1)+ksm(n,m/2));
cout<<ans/2/m<<endl;
n=read(),m=read();
}
return 0;
}
//Now,only God know
//但行好事，莫问前程
《论提前开金牌导航的感觉》
标签：include,int,ybtoj,ans,珠子,while,项链,例题 来源： https://blog.csdn.net/weixin_49843717/article/details/117933247
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。